Seismic modeling by combining the finite-difference scheme with the numerical dispersion suppression neural network

Seismic finite-difference(FD) modeling suffers from numerical dispersion including both the temporal and spatial dispersion, which can decrease the accuracy of the numerical modeling. To improve the accuracy and efficiency of the conventional numerical modeling, I develop a new seismic modeling method by combining the FD scheme with the numerical dispersion suppression neural network(NDSNN). This method involves the following steps. First, a training data set composed of a small number of wavefield snapshots is generated. The wavefield snapshots with the low-accuracy wavefield data and the high-accuracy wavefield data are paired, and the low-accuracy wavefield snapshots involve the obvious numerical dispersion including both the temporal and spatial dispersion. Second, the NDSNN is trained until the network converges to simultaneously suppress the temporal and spatial dispersion.Third, the entire set of low-accuracy wavefield data is computed quickly using FD modeling with the large time step and the coarse grid. Fourth, the NDSNN is applied to the entire set of low-accuracy wavefield data to suppress the numerical dispersion including the temporal and spatial dispersion.Numerical modeling examples verify the effectiveness of my proposed method in improving the computational accuracy and efficiency.

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.

Seismic modelling of tracked-vehicle signals for monitoring and verifi cation

To better understand characteristics of seismic signals of tracked vehicles measured when passing a sensor line,we numerically modelled force-pulse responses of a layered soil that is similar in its seismic properties to that found at the original measurement site.The vertical-force pulses from the road wheels rolling over the track elements are fitted to the measured ones.Single-pulse seismic waves vary with distance due to diff erent wave types,refl ections at layer boundaries,vehicle velocity and relative position of the left and right track elements.They are computed by a modelling program and superposed at sensor positions with the appropriate slant distance and time shift for each track element.These sum signals are in qualitative agreement with those from the original measurements.However,they are several magnitudes weaker and much smoother.Furthermore,higher frequencies are damped much less at larger distances.Due to the large variability of the sum signals,recognition of tracked-vehicle types exclusively through their seismic signals seems diffi cult.

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.

3D seismic forward modeling from the multiphysical inversion at the Ketzin CO_(2) storage site

From June 2008 to August 2013,approximately 67 kt of CO_(2) was injected into a deep saline formation at the Ketzin pilot CO_(2) storage site.During injection,3D seismic surveys have been performed to monitor the migration of sequestered CO_(2).Seismic monitoring results are limited by the acquisition and signal-to-noise ratio of the acquired data.The multiphysical reservoir simulation provides information regarding the CO_(2) fluid behavior,and the approximated model should be calibrated with the monitoring results.In this work,property models are delivered from the multiphysical model during 3D repeated seismic surveys.The simulated seismic data based on the models are compared with the real data,and the results validate the effectiveness of the multiphysical inversion method.Time-lapse analysis shows the trend of CO_(2) migration during and after injection.

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.

Seismic wavefield modeling based on time-domain symplectic and Fourier finite-difference method

Seismic wavefield modeling is important for improving seismic data processing and interpretation. Calculations of wavefield propagation are sometimes not stable when forward modeling of seismic wave uses large time steps for long times. Based on the Hamiltonian expression of the acoustic wave equation, we propose a structure-preserving method for seismic wavefield modeling by applying the symplectic finite-difference method on time grids and the Fourier finite-difference method on space grids to solve the acoustic wave equation. The proposed method is called the symplectic Fourier finite-difference （symplectic FFD） method, and offers high computational accuracy and improves the computational stability. Using acoustic approximation, we extend the method to anisotropic media. We discuss the calculations in the symplectic FFD method for seismic wavefield modeling of isotropic and anisotropic media, and use the BP salt model and BP TTI model to test the proposed method. The numerical examples suggest that the proposed method can be used in seismic modeling of strongly variable velocities, offering high computational accuracy and low numerical dispersion. The symplectic FFD method overcomes the residual qSV wave of seismic modeling in anisotropic media and maintains the stability of the wavefield propagation for large time steps.

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.

Brittleness index and seismic rock physics model for anisotropic tight-oil sandstone reservoirs

Brittleness analysis becomes important when looking for sweet spots in tightoil sandstone reservoirs. Hence, appropriate indices are required as accurate brittleness evaluation criteria. We construct a seismic rock physics model for tight-oil sandstone reservoirs with vertical fractures. Because of the complexities in lithology and pore structure and the anisotropic characteristics of tight-oil sandstone reservoirs, the proposed model is based on the solid components, pore connectivity, pore type, and fractures to better describe the sandstone reservoir microstructure. Using the model, we analyze the brittleness sensitivity of the elastic parameters in an anisotropic medium and establish a new brittleness index. We show the applicability of the proposed brittleness index for tight-oil sandstone reservoirs by considering the brittleness sensitivity, the rock physics response characteristics, and cross-plots. Compared with conventional brittleness indexes, the new brittleness index has high brittleness sensitivity and it is the highest in oil-bearing brittle zones with relatively high porosity. The results also suggest that the new brittleness index is much more sensitive to elastic properties variations, and thus can presumably better predict the brittleness characteristics of sweet spots in tight-oil sandstone reservoirs.

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

Mathematical geophone （MG） and equal-time stacking （ETS） principles are used to implement seismic prestack forward modeling with irregular surfaces using the oneway acoustic wave-equation. This method receives seismic primary reflections from the subsurface using a set of virtual MGs. The receivers can be located anywhere on an irregular observing surface. Moreover, the ETS method utilizes the one-way acoustic wave equation to easily and quickly image and extrapolate seismic reflection data. The method is illustrated using high single-noise ratio common shot gathers computed by numerical forward modeling of two simple models, one with a flat surface and one with an irregular surface, and a complex normal fault model. A prestack depth migration method for irregular surface topography was used to reoroduce the normal fault model with high accuracy.

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.

Seismic physical modeling and quality factor

Accurate Q parameter is hard to be obtained, but there is great difference between Q measurements from different measurement methods in seismic physical modelling. The influence factors, stability and accuracy of different methods are analyzed through standard sample experiment and the seismic physical modelling. Based on this, we proposed an improved method for improving accuracy of pulse transmission method, in which the samples with similar acoustic properties to the test sample are selected as the reference samples. We assess the stability and accuracy of the pulse transmission, pulse transmission insertion, and reflection wave methods for obtaining the quality factor Q using standard and reference samples and seismic physical modeling. The results suggest that the Q-values obtained by the pulse transmission method are strongly affected by diffraction and the error is 50% or greater, whereas the relative error of the improved pulse transmission method is about 10%. By using a theoretical diffraction correction method and the improved measurement method, the differences among the Q-measuring methods can be limited to within 10%.

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.

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.

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.