Fractured reservoir modeling by discrete fracture network and seismic modeling in the Tarim Basin,China

Fractured reservoirs are an important target for oil and gas exploration in the Tarim Basin and the prediction of this type of reservoir is challenging.Due to the complicated fracture system in the Tarim Basin,the conventional AVO inversion method based on HTI theory to predict fracture development will result in some errors.Thus,an integrated research concept for fractured reservoir prediction is put forward in this paper.Seismic modeling plays a bridging role in this concept,and the establishment of an anisotropic fracture model by Discrete Fracture Network （DFN） is the key part.Because the fracture system in the Tarim Basin shows complex anisotropic characteristics,it is vital to build an effective anisotropic model.Based on geological,well logging and seismic data,an effective anisotropic model of complex fracture systems can be set up with the DFN method.The effective elastic coefficients,and the input data for seismic modeling can be calculated.Then seismic modeling based on this model is performed,and the seismic response characteristics are analyzed.The modeling results can be used in the following AVO inversion for fracture detection.

Optimization of a precise integration method for seismic modeling based on graphic processing unit

General purpose graphic processing unit （GPU） calculation technology is gradually widely used in various fields. Its mode of single instruction, multiple threads is capable of seismic numerical simulation which has a huge quantity of data and calculation steps. In this study, we introduce a GPU-based parallel calculation method of a precise integration method （PIM） for seismic forward modeling. Compared with CPU single-core calculation, GPU parallel calculating perfectly keeps the features of PIM, which has small bandwidth, high accuracy and capability of modeling complex substructures, and GPU calculation brings high computational efficiency, which means that high-performing GPU parallel calculation can make seismic forward modeling closer to real seismic records.

Migrated Exploding Reflectors in Evaluation of Finite Difference Solution for Inhomogeneous Seismic Models

Earth is inhomogeneous, which means its elastic characteristics change with depth. The seismic method employs the propagation of waves throughout the earth to locate different structures and stratigraphy. Understanding the wave propagation is an important matter in exploration seismology;therefore modeling of seismic wave is an important tool. To validate the interpreted earth model out of the seismic data, seismic synthetic seismograms should be generated in a process named “seismic forward modeling”. Finite difference method is used as one of the most common numerical modeling techniques. In this paper the accuracy of finite difference method in seismic section modeling is explored on different modeled data set of heterogeneous earth. It is shown that finite difference method completes with migration to reposition the events in their correct location. Two different migration methods are used and various velocities are also tested to determine an appropriate migration velocity. Finally the validly of finite difference modeling is examined using a 2D structural similarity index technique.

Review of Artificial Intelligence for Oil and Gas Exploration: Convolutional Neural Network Approaches and the U-Net 3D Model

Deep learning, especially through convolutional neural networks (CNN) such as the U-Net 3D model, has revolutionized fault identification from seismic data, representing a significant leap over traditional methods. Our review traces the evolution of CNN, emphasizing the adaptation and capabilities of the U-Net 3D model in automating seismic fault delineation with unprecedented accuracy. We find: 1) The transition from basic neural networks to sophisticated CNN has enabled remarkable advancements in image recognition, which are directly applicable to analyzing seismic data. The U-Net 3D model, with its innovative architecture, exemplifies this progress by providing a method for detailed and accurate fault detection with reduced manual interpretation bias. 2) The U-Net 3D model has demonstrated its superiority over traditional fault identification methods in several key areas: it has enhanced interpretation accuracy, increased operational efficiency, and reduced the subjectivity of manual methods. 3) Despite these achievements, challenges such as the need for effective data preprocessing, acquisition of high-quality annotated datasets, and achieving model generalization across different geological conditions remain. Future research should therefore focus on developing more complex network architectures and innovative training strategies to refine fault identification performance further. Our findings confirm the transformative potential of deep learning, particularly CNN like the U-Net 3D model, in geosciences, advocating for its broader integration to revolutionize geological exploration and seismic analysis.

Amplitude Variation with Offset Responses Modeling Study of Walkaway Vertical Seismic Profile Data at CO_2 Geological Storage Site,Ketzin,Germany

An important component of any CO_2 sequestration project is seismic monitoring for tracking changes in subsurface physical properties,such as velocity and density.Different reservoirs have different amplitude variation with offset(AVO) responses,which can define underground conditions. In the present paper we investigate walkaway vertical seismic profile(VSP) AVO response to CO_2 injection at the Ketzin site,the first European onshore CO_2 sequestration pilot study dealing with research on geological storage of CO_2.First,we performed rock physics analysis to evaluate the effect of injected CO_2 on seismic velocity using the Biot-Gassmann equation.On the basis of this model,the seismic response for different CO_2 injection saturation was studied using ray tracing modeling.We then created synthetic walkaway VSP data,which we then processed.In contrast,synthetic seismic traces were created from borehole data.Finally,we found that the amplitude of CO_2 injected sand layer with different gas saturations were increased with the offset when compared with the original brine target layer.This is the typical classⅢAVO anomaly for gas sand layer.The AVO responses matched the synthetic seismic traces very well.Therefore,walkaway VSP AVO response can monitor CO_2 distribution in the Ketzin area.

The effect of pore fluid on seismicity: a computer model

The influence of fluid on seismicity of a computerized system is analyzed in this paper. The diffusion equation of fluid in a crustal fault area is developed and used in the calculation of a spring slide damper model. With mirror imagin boundary condition and three initial conditions, the equation is solved for a dynamic model that consists of six seismic belts and eight seismogenous sources in each belt with both explicit algorithm and implicit algorithm. The analysis of the model with water sources shows that the implicit algorithm is better to be used to calculate the model. Taking a constant proportion of the pore pressure of a broken element to that of its neighboring elements, the seismicity of the model is calculated with mirror boundary condition and no water source initial condition. The results shows that the frequency and magnitude of shocks are both higher than those in the model with no water pore pressure, which provides more complexity to earthquake prediction.

Ultrasonic sensitivity-improved fiber-optic Fabry-Perot interferometer using a beam collimator and its application for ultrasonic imaging of seismic physical models

An ultrasonic sensitivity-improved fiber-optic Fabry-Perot interferometer （FPI） is proposed and employed for ultra- sonic imaging of seismic physical models （SPMs）. The FPI comprises a flexible ultra-thin gold film and the end face of a graded-index multimode fiber （MMF）, both of which are enclosed in a ceramic tube. The MMF in a specified length can collimate the diverged light beam and compensate for the light loss inside the air cavity, leading to an increased spectral fringe visibility and thus a steeper spectral slope. By using the spectral sideband filtering technique, the collimated FP1 shows an improved ultrasonic response. Moreover, two-dimensional images of two SPMs are achieved in air by recon- structing the pulse-echo signals through using the time-of-flight approach. The proposed sensor with easy fabrication and compact size can be a good candidate for high-sensitivity and high-precision nondestructive testing of SPMs.

A GIS-based time-dependent seismic source modeling of Northern Iran

The first step in any seismic hazard study is the definition of seismogenic sources and the estimation of magnitude-frequency relationships for each source. There is as yet no standard methodology for source modeling and many researchers have worked on this topic. This study is an effort to define linear and area seismic sources for Northern Iran. The linear or fault sources are developed based on tectonic features and characteristic earthquakes while the area sources are developed based on spatial distribution of small to moderate earthquakes. Time-dependent recurrence relationships are developed for fault sources using renewal approach while time-independent frequency-magnitude relationships are proposed for area sources based on Poisson process. GIS functionalities are used in this study to introduce and incorporate spatial- temporal and geostatistical indices in delineating area seismic sources. The proposed methodology is used to model seismic sources for an area of about 500 by 400 square kilometers around Tehran. Previous researches and reports are studied to compile an earthquake/fault catalog that is as complete as possible. All events are transformed to uniform magnitude scale; duplicate events and dependent shocks are removed. Completeness and time distribution of the compiled catalog is taken into account. The proposed area and linear seismic sources in conjunction with defined recurrence relationships can be used to develop time-dependent probabilistic seismic hazard analysis of Northern Iran.

Seismicity acceleration model and its application to several earthquake regions in China

With the theory of subcritical crack growth, we can deduce the fundamental equation of regional seismicity acceleration model. Applying this model to intraplate earthquake regions, we select three earthquake subplates: North China Subplate, Chuan Dian Block and Xinjiang Subplate, and divide the three subplates into seven researched regions by the difference of seismicity and tectonic conditions. With the modified equation given by Sornette and Sammis (1995), we analysis the seismicity of each region. To those strong earthquakes already occurred in these region, the model can give close fitting of magnitude and occurrence time, and the result in this article indicates that the seismicity acceleration model can also be used for describing the seismicity of intraplate. In the article, we give the magnitude and occurrence time of possible strong earthquakes in Shanxi, Ordos, Bole Tuokexun, Ayinke Wuqia earthquake regions. In the same subplate or block, the earthquake periods for each earthquake region are similar in time interval. The constant α in model can be used to describe the intensity of regional seismicity, and for the Chinese Mainland, α is 0.4 generally. To the seismicity in Taiwan and other regions with complex tectonic conditions, the model does not fit well at present.

Seismic Physical Modeling Technology and Its Applications

This paper introduces the seismic physical modeling technology in the CNPC Key Lab of Geophysical Exploration. It includes the seismic physical model positioning system, the data acquisition system, sources, transducers, model materials, model building techniques, precision measurements of model geometry, the basic principles of the seismic physical modeling and experimental methods, and two physical model examples.

Accuracy of the staggered-grid finite-difference method of the acoustic wave equation for marine seismic reflection modeling

Seismic wave modeling is a cornerstone of geophysical data acquisition, processing, and interpretation, for which finite-difference methods are often applied. In this paper, we extend the velocity- pressure formulation of the acoustic wave equation to marine seismic modeling using the staggered-grid finite-difference method. The scheme is developed using a fourth-order spatial and a second-order temporal operator. Then, we define a stability coefficient (SC) and calculate its maximum value under the stability condition. Based on the dispersion relationship, we conduct a detailed dispersion analysis for submarine sediments in terms of the phase and group velocity over a range of angles, stability coefficients, and orders. We also compare the numerical solution with the exact solution for a P-wave line source in a homogeneous submarine model. Additionally, the numerical results determined by a Marmousi2 model with a rugged seafloor indicate that this method is sufficient for modeling complex submarine structures.

One-way wave equation seismic prestack forward modeling with irregular surfaces

数学地震检波器(MG )(et ) 并且相等时间的叠原则被用来实现向前用单程的声学的波浪方程与不规则的表面当模特儿的地震 prestack。这个方法从表面下的使用受到地震主要思考一套虚拟 MG。接收装置能位于在任何地方不规则的观察表面。而且， ETS 方法利用单程的声学的波浪方程到容易并且快速想象并且外推地震思考数据。方法用普通射击集合由二个简单模型，的数字前面的建模计算了的高单个噪音的比率被说明有扁平的表面的和有不规则的表面的，和复杂正常指责模型。为不规则的表面地形学的一个 prestack 深度迁居方法被用来与高精确性复制正常差错模型。

Research on percolation model and criticality of seismicity

Making use of modern nonlinear physics theory and earthquake focus theory, combined with seismicity characteristics, the percolation model of earthquake activity is given in this paper. We take the seismogenic process of alarge earthquake as a phase transition process of percolation and apply the renormalization method to phase transition of percolation. The critical property of the system, which is like percolation probability exponential andcorrelative length exponential, etc, can be calculated under the fixed point as which in the renormalization transformation infinite correlative length in percolation phase transition is taken. The percolation phase transition process of two large earthquakes, which are Haicheng and Tangshan event occurred in 1975 and 1976 respectively, hasbeen discussed by means of seismicity data before and after two shocks.

Fault activation and induced seismicity in geological carbon storage--Lessons learned from recent modeling studies

In the light of current concerns related to induced seismicity associated with geological carbon sequestration(GCS),this paper summarizes lessons learned from recent modeling studies on fault activation,induced seismicity,and potential for leakage associated with deep underground carbon dioxide(CO2) injection.Model simulations demonstrate that seismic events large enough to be felt by humans require brittle fault properties and continuous fault permeability allowing pressure to be distributed over a large fault patch to be ruptured at once.Heterogeneous fault properties,which are commonly encountered in faults intersecting multilayered shale/sandstone sequences,effectively reduce the likelihood of inducing felt seismicity and also effectively impede upward CO2leakage.A number of simulations show that even a sizable seismic event that could be felt may not be capable of opening a new flow path across the entire thickness of an overlying caprock and it is very unlikely to cross a system of multiple overlying caprock units.Site-specific model simulations of the In Salah CO2storage demonstration site showed that deep fractured zone responses and associated microseismicity occurred in the brittle fractured sandstone reservoir,but at a very substantial reservoir overpressure close to the magnitude of the least principal stress.We conclude by emphasizing the importance of site investigation to characterize rock properties and if at all possible to avoid brittle rock such as proximity of crystalline basement or sites in hard and brittle sedimentary sequences that are more prone to injection-induced seismicity and permanent damage.

Improvements in seismic event locations in a deep western U.S. coal mine using tomographic velocity models and an evolutionary search algorithm

Methods of improving seismic event locations were investigated as part of a research study aimed at reducing ground control safety hazards. Seismic event waveforms collected with a 23-station three-dimensional sensor array during longwall coal mining provide the data set used in the analyses. A spatially variable seismic velocity model is constructed using seismic event sources in a passive tomographic method. The resulting three-dimensional velocity model is used to relocate seismic event positions. An evolutionary optimization algorithm is implemented and used in both the velocity model development and in seeking improved event location solutions. Results obtained using the different velocity models are compared. The combination of the tomographic velocity model development and evolutionary search algorithm provides improvement to the event locations.

A limit equilibrium fracture zone model to investigate seismicity in coal mines

This paper explores possible synergies between techniques used to minimise seismicity in deep South African gold mines and their applicability to control coal bumps. The paper gives a summary of the techniques used in the deep gold mines and a critical appraisal if these are useful in coal mines. The techniques typically include control of mining rate, preconditioning, optimisation of extraction sequences and centralised blasting. Of particular interest to the coal bump problem is an experimental limit equilibrium fracture zone model implemented in a displacement discontinuity code. This was recently developed for the gold mines to enable the interactive analysis of complex tabular mine layout extraction sequences. The model specifically accommodates energy dissipation computations in the developing fracture zone near the edges of these excavations. This allows the released energy to be used as a surrogate measure of ongoing seismic activity and addresses a number of the weaknesses in the traditional usage of this quantity as a criterion for the design of seismically active layouts. This paper investigates the application of the model to a hypothetical coal longwall layout and the specific problem of coal bumps.

Seismic fragility analysis of bridges by relevance vector machine based demand prediction model

A relevance vector machine(RVM)based demand prediction model is explored for efficient seismic fragility analysis(SFA)of a bridge structure.The proposed RVM model integrates both record-to-record variations of ground motions and uncertainties of parameters characterizing the bridge model.For efficient fragility computation,ground motion intensity is included as an added dimension to the demand prediction model.To incorporate different sources of uncertainty,random realizations of different structural parameters are generated using Latin hypercube sampling technique.Mean fragility,along with its dispersions,is estimated based on the log-normal fragility model for different critical components of a bridge.The effectiveness of the proposed RVM model-based SFA of a bridge structure is elucidated numerically by comparing it with fragility results obtained by the commonly used SFA approaches,while considering the most accurate direct Monte Carlo simulation-based fragility estimates as the benchmark.The proposed RVM model provides a more accurate estimate of fragility than conventional approaches,with significantly less computational effort.In addition,the proposed model provides a measure of uncertainty in fragility estimates by constructing confidence intervals for the fragility curves.

Numerical Modeling of Neotectonic Movements and State of Stresses in the Central Seismic Gap Region,Garhwal Himalaya

This paper presents finite element modeling(FEM) to simulate the present-day stress field and crustal deformation using NE-SW structural section in the central Seismic Gap region of the Garhwal Himalaya.Our study deals with the effect of geometrical characteristics and rock layer parameters on the upper crust.Modeling results show that two types of tectonic regimes developed in the central Seismic Gap region:the geotectonics of the northern part has been controlled by regional compression,whereas southern part is characterized by regional extension.Correspondingly,thrust faults are induced in the northern part and normal faults are extensively developed in the southern front.Those evidences noticeably indicate that the compressive tectonic environment of the Himalaya becomes change into the extensional tectonic regime in its front.The computed shear stress accumulation along the northern flat of Main Himalayan Thrust(MHT) implies that considerable amount of interseismic stress is building up along the MHT system in the Himalaya,which ultimately release through the possible future great Himalayan earthquake(M>8).The comparison between our modeled stress field,faulting pattern and horizontal shortening rate with the distribution of the microseismic events,focal mechanism solutions,active faulting and GPS data in the central Seismic Gap region shows good agreement.

An improved multidirectional velocity model for micro-seismic monitoring in rock engineering

An improved multidirectional velocity model was proposed for more accurately locating micro-seismic events in rock engineering. It was assumed that the stress wave propagation velocities from a micro-seismic source to three nearest monitoring sensors in a sensor's array arrangement were the same. Since the defined objective function does not require pre-measurement of the stress wave propagation velocity in the field, errors from the velocity measurement can be avoided in comparison to three traditional velocity models. By analyzing 24 different cases, the proposed multidirectional velocity model iterated by the Simplex method is found to be the best option no matter the source is within the region of the sensor's array or not. The proposed model and the adopted iterative algorithm are verified by field data and it is concluded that it can significantly reduce the error of the estimated source location.

Parameter Fitting of Seismic Intensity Attenuation Model in Xinjiang

By using the existing historical earthquake investigation data in Xinjiang,this paper obtained the envelope curves of isoseismal maps of 103 destructive earthquakes occurring from 1716 to 2010 after digitization of the data. The author summarized the seismic intensity attenuation laws in the Xinjiang region with the multiple regression fitting method. The intensity attenuation function of the elliptical model was provided and the fitting results in different periods and areas were compared. Finally, the intensity attenuation relationship in the Xinjiang region was obtained by the method of constraining the start and end of the attenuation curves.