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 effic...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 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 con...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.展开更多
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 ...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.展开更多
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 migr...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.展开更多
Existing studies indicate that gas hydrate-bearing formations exhibit notable seismic velocity dispersion and attenuation. The Shenhu area of the South China Sea hold significant gas hydrate resource potential;however...Existing studies indicate that gas hydrate-bearing formations exhibit notable seismic velocity dispersion and attenuation. The Shenhu area of the South China Sea hold significant gas hydrate resource potential;however, the relationship between seismic velocity dispersion, attenuation properties, and gas-hydrate saturation remains insufficiently understood. Furthermore, a significant mismatch exists between the real seismic angle gather near a well and the synthetic angle gather generated using the convolution method, and this discrepancy may arise from the seismic velocity dispersion and attenuation characteristics of the gas hydrate-bearing formations. In this paper, we develop a rock physics model that integrates White's and Dvorkin's models, accounting for varied types of gas-hydrate occurrence states,specifically tailored to the gas hydrate-bearing formations in the Shenhu area. This model is calibrated with well log data and employed to investigate how gas-hydrate saturation influences seismic velocity dispersion and attenuation. Numerical analysis reveals the coexistence of two types of gas-hydrate occurrence states in the region: high gas-hydrate saturation formations are dominated by loadbearing-type gas hydrate, and formations containing both gas hydrate and free gas may exhibit either load-bearing or pore-filling types. The seismic velocity dispersion and attenuation properties vary significantly depending on the gas-hydrate occurrence state. We further apply the proposed model to generate seismic velocity and attenuation logs at various frequencies. These logs are used in seismic forward modeling employing both the convolution method and the propagator matrix method. Well tie analysis indicates that the synthetic angle gather incorporating attenuation via the propagator matrix method aligns more closely with the real seismic angle gather than the convolution method. This study provides valuable insights into frequency-dependent amplitude versus offset(AVO) analysis and the seismic interpretation of gas hydrate-bearing formations in the South China Sea.展开更多
The Longmenshan(LMS)fault zone is located at the junction of the eastern Tibetan Plateau and the Sichuan Basin and is of great significance for studying regional tectonics and earthquake hazards.Although regional velo...The Longmenshan(LMS)fault zone is located at the junction of the eastern Tibetan Plateau and the Sichuan Basin and is of great significance for studying regional tectonics and earthquake hazards.Although regional velocity models are available for the LMS fault zone,high-resolution velocity models are lacking.Therefore,a dense array of 240 short-period seismometers was deployed around the central segment of the LMS fault zone for approximately 30 days to monitor earthquakes and characterize fine structures of the fault zone.Considering the large quantity of observed seismic data,the data processing workflow consisted of deep learning-based automatic earthquake detection,phase arrival picking,and association.Compared with the earthquake catalog released by the China Earthquake Administration,many more earthquakes were detected by the dense array.Double-difference seismic tomography was adopted to determine V_(p),V_(s),and V_(p)/V_(s)models as well as earthquake locations.The checkerboard test showed that the velocity models have spatial resolutions of approximately 5 km in the horizontal directions and 2 km at depth.To the west of the Yingxiu–Beichuan Fault(YBF),the Precambrian Pengguan complex,where most of earthquakes occurred,is characterized by high velocity and low V_(p)/V_(s)values.In comparison,to the east of the YBF,the Upper Paleozoic to Jurassic sediments,where few earthquakes occurred,show low velocity and high V_(p)/V_(s)values.Our results suggest that the earthquake activity in the LMS fault zone is controlled by the strength of the rock compositions.When the high-resolution velocity models were combined with the relocated earthquakes,we were also able to delineate the fault geometry for different faults in the LMS fault zone.展开更多
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 st...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.展开更多
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 variatio...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.展开更多
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 sei...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.展开更多
： 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 modeli...： 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.展开更多
The perfectly matched layer（PML） was first introduced by Berenger as an absorbing boundary condition for electromagnetic wave propagation.In this article,a method is developed to ex-tend the PML to simulating seismi...The perfectly matched layer（PML） was first introduced by Berenger as an absorbing boundary condition for electromagnetic wave propagation.In this article,a method is developed to ex-tend the PML to simulating seismic wave propagation in fluid-saturated porous medium.This non-physical boundary is used at the computational edge of a Forsyte polynomial convolutional differenti-ator（FPCD） algorithm as an absorbing boundary condition to truncate unbounded media.The incor-poration of PML in Biot＇s equations is given.Numerical results show that the PML absorbing bound-ary condition attenuates the outgoing waves effectively and eliminates the reflections adequately.展开更多
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, transduce...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.展开更多
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 formulatio...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.展开更多
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 diff...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%.展开更多
Fractured reservoirs have always been a big favorable area for oil and gas reservoirs,so prediction of fractures is also a research hotspot in recent years.Due to the diversity of fracture development and the unclear ...Fractured reservoirs have always been a big favorable area for oil and gas reservoirs,so prediction of fractures is also a research hotspot in recent years.Due to the diversity of fracture development and the unclear development mechanism,fracture prediction has always been a major problem.Simple numerical simulation In this paper,seismic attribute is combined with numerical simulation,logging data and actual seismic profile are used as constraints,inversion impedance value and coherent attribute are combined,and finally a property model more in line with the actual geological conditions is established.The wave equation calculation and migration processing were used to obtain the numerical simulation profile,and the actual seismic profile,fracture detection profile and numerical simulation profile were combined for analysis:①The numerical simulation section under this modeling method can greatly correspond to the actual seismic section,and the reflected results can better reflect the changes of response characteristics.②The reliability and applicability of the fracture detection technology can be determined by comparing the forward simulation profile with the fracture detection profile.展开更多
To study the ground motion intensity measures(IMs)suitable for the design of seismic performance with a focus on longitudinal resistance in tunnel structures,21 different seismic intensity parameters are selected for ...To study the ground motion intensity measures(IMs)suitable for the design of seismic performance with a focus on longitudinal resistance in tunnel structures,21 different seismic intensity parameters are selected for nonlinear calculation and analysis of tunnel structures,in order to determine the optimal IM for the longitudinal seismic performance of tunnel structures under different site conditions.An improved nonlinear beam-spring model is developed to calculate the longitudinal seismic response of tunnels.The PQ-Fiber model is used to simulate the longitudinal nonlinear behavior of tunnel structures and the tangential interactions between the tunnel and the soil is realized by load in the form of moment.Five different site types are considered and 21 IMs is evaluated against four criteria:effectiveness,practicality,usefulness,and sufficiency.The results indicate that the optimal IMs are significantly influenced by the site conditions.Specifically,sustained maximum velocity(V_(SM))emerges as the optimal IM for circular tunnels in soft soil conditions(CaseⅠsites),peak ground velocity(V PG)is best suited for CaseⅡsites,sustained maximum acceleration(A_(SM))is ideal for both CaseⅢand CaseⅤsites,and peak ground acceleration(A PG)for CaseⅣsites.As site conditions transition from CaseⅠto CaseⅤ,from soft to hard,the applicability of acceleration-type intensity parameters gradually decreases,while the applicability of velocity-type intensity parameters gradually increases.展开更多
Large-scale 3D physical models of complex structures can be used to simulate hydrocarbon exploration areas. The high-fidelity simulation of actual structures poses challenges to model building and quality control. Suc...Large-scale 3D physical models of complex structures can be used to simulate hydrocarbon exploration areas. The high-fidelity simulation of actual structures poses challenges to model building and quality control. Such models can be used to collect wideazimuth, multi-azimuth, and full-azimuth seismic data that can be used to verify various 3D processing and interpretation methods. Faced with nonideal imaging problems owing to the extensive complex surface conditions and subsurface structures in the oil-rich foreland basins of western China, we designed and built the KS physical model based on the complex subsurface structure. This is the largest and most complex 3D physical model built to date. The physical modeling technology advancements mainly involve 1） the model design method, 2） the model casting flow, and 3） data acquisition. A 3D velocity model of the physical model was obtained for the first time, and the model building precision was quantitatively analyzed. The absolute error was less than 3 mm, which satisfies the experimental requirements. The 3D velocity model obtained from 3D measurements of the model layers is the basis for testing various imaging methods. Furthermore, the model is considered a standard in seismic physical modeling technology.展开更多
The pseudospectral method has been applied to the simulation of seismic wave propagation in 2-D global Earth model. When a whole Earth model is considered, the center of the Earth is included in the model and then sin...The pseudospectral method has been applied to the simulation of seismic wave propagation in 2-D global Earth model. When a whole Earth model is considered, the center of the Earth is included in the model and then singularity arises at the center of the Earth where r=0 since the 1/r term appears in the wave equations. In this paper, we extended the global seismic wavefield simulation algorithm for regular grid mesh to staggered grid configuration and developed a scheme to solve the numerical problems associated with the above singularity for a 2-D global Earth model defined on staggered grid using pseudospectral method. This scheme uses a coordinate transformation at the center of the model, in which the field variables at the center are calculated in Cartesian coordinates from the values on the grids around the center. It allows wave propagation through the center and hence the wavefield at the center can be stably calculated. Validity and accuracy of the scheme was tested by compared with the discrete wavenumber method. This scheme could also be suitable for other numerical methods or models parameterized in cylindrical or spherical coordinates when singularity arises at the center of the model.展开更多
Strong near-fault ground motion, usually caused by the fault-rupture and characterized by a pulse-like velocity- wave form, often causes dramatic instantaneous seismic energy （Jadhav and Jangid 2006）. Some reinforce...Strong near-fault ground motion, usually caused by the fault-rupture and characterized by a pulse-like velocity- wave form, often causes dramatic instantaneous seismic energy （Jadhav and Jangid 2006）. Some reinforced concrete （RC） bridge columns, even those built according to ductile design principles, were damaged in the 1999 Chi-Chi earthquake. Thus, it is very important to evaluate the seismic response of a RC bridge column to improve its seismic design and prevent future damage. Nonlinear time history analysis using step-by-step integration is capable of tracing the dynamic response of a structure during the entire vibration period and is able to accommodate the pulsing wave form. However, the accuracy of the numerical results is very sensitive to the modeling of the nonlinear load-deformation relationship of the structural member. FEMA 273 and ATC-40 provide the modeling parameters for structural nonlinear analyses of RC beams and RC columns. They use three parameters to define the plastic rotation angles and a residual strength ratio to describe the nonlinear load- deformation relationship of an RC member. Structural nonlinear analyses are performed based on these parameters. This method provides a convenient way to obtain the nonlinear seismic responses of RC structures. However, the accuracy of the numerical solutions might be further improved. For this purpose, results from a previous study on modeling of the static pushover analyses for RC bridge columns （Sung et al. 2005） is adopted for the nonlinear time history analysis presented herein to evaluate the structural responses excited by a near-fault ground motion. To ensure the reliability of this approach, the numerical results were compared to experimental results. The results confirm that the proposed approach is valid.展开更多
The Chilean Pampean flat slab subduction segment is characterized by the nearly horizontal subduction of the Nazca Plate within the depth range of 100-120 km.Numerous seismic tomography studies have been conducted to ...The Chilean Pampean flat slab subduction segment is characterized by the nearly horizontal subduction of the Nazca Plate within the depth range of 100-120 km.Numerous seismic tomography studies have been conducted to investigate its velocity structure;however,they have used only seismic body wave data or surface wave data.As a result,the existing velocity models in the region may have relatively large uncertainties.In this study,we use body wave arrival times from earthquakes occurring in central Chile between 2014 and 2019,as well as Rayleigh wave phase velocity maps at periods of 5-80 s from ambient noise empirical Green’s functions in Chile.By jointly using body wave arrival times and surface wave dispersion data,we refine the VS model and improve earthquake locations in the central Chile subduction zone.Compared with previous velocity models,our velocity model better reveals an eastward-dipping high-velocity plate representing the subducting Nazca Plate,which is 40-50 km thick and is more consistent with the slab thickness estimated by receiver function imaging and thermal modeling.Overall,the intraslab seismicity distribution spatially correlates well with the slab high-velocity anomalies except along the subduction paths of the CopiapóRidge and Juan Fernández Ridge.Additionally,parallel low-velocity stripes are imaged beneath the subducting plate,which are likely associated with the accumulated melts.The joint inversion velocity model also resolves widespread low-velocity anomalies in the crust beneath the Central Volcanic Zone of the central Andes,likely representing crustal magma chambers for various volcanoes.展开更多
基金supported by the National Natural Science Foundation of China (grant numbers: 41874160 and 92055213)。
文摘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.
基金co-supported by the National Basic Research Program of China(Grant No.2011CB201103)the National Science and Technology Major Project(GrantNo.2011ZX05004003)
文摘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.
基金supported by the National Natural Science Foundation of China (Nos 40974066 and 40821062)National Basic Research Program of China (No 2007CB209602)
文摘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.
基金supported by the National Natural Science Foundation of China(Grant No.42025403)the Youth Innovation Promotion Association,Chinese Academy of Sciences(Grant No.2023074).
文摘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.
基金supported by National Natural Science Foundation of China(W2431028,42122029)SINOPEC Fundamental Research Program(P24258)CNPC Investigations on fundamental experiments and advanced theoretical methods in geophysical prospecting applications(2022DQ0604-02).
文摘Existing studies indicate that gas hydrate-bearing formations exhibit notable seismic velocity dispersion and attenuation. The Shenhu area of the South China Sea hold significant gas hydrate resource potential;however, the relationship between seismic velocity dispersion, attenuation properties, and gas-hydrate saturation remains insufficiently understood. Furthermore, a significant mismatch exists between the real seismic angle gather near a well and the synthetic angle gather generated using the convolution method, and this discrepancy may arise from the seismic velocity dispersion and attenuation characteristics of the gas hydrate-bearing formations. In this paper, we develop a rock physics model that integrates White's and Dvorkin's models, accounting for varied types of gas-hydrate occurrence states,specifically tailored to the gas hydrate-bearing formations in the Shenhu area. This model is calibrated with well log data and employed to investigate how gas-hydrate saturation influences seismic velocity dispersion and attenuation. Numerical analysis reveals the coexistence of two types of gas-hydrate occurrence states in the region: high gas-hydrate saturation formations are dominated by loadbearing-type gas hydrate, and formations containing both gas hydrate and free gas may exhibit either load-bearing or pore-filling types. The seismic velocity dispersion and attenuation properties vary significantly depending on the gas-hydrate occurrence state. We further apply the proposed model to generate seismic velocity and attenuation logs at various frequencies. These logs are used in seismic forward modeling employing both the convolution method and the propagator matrix method. Well tie analysis indicates that the synthetic angle gather incorporating attenuation via the propagator matrix method aligns more closely with the real seismic angle gather than the convolution method. This study provides valuable insights into frequency-dependent amplitude versus offset(AVO) analysis and the seismic interpretation of gas hydrate-bearing formations in the South China Sea.
基金supported by the Scientific Research Foundation for High-level Talents of Anhui University of Science and Technology under Grant 2024yjrc64the National Key R&D Program of China under Grant 2018YFC1504102。
文摘The Longmenshan(LMS)fault zone is located at the junction of the eastern Tibetan Plateau and the Sichuan Basin and is of great significance for studying regional tectonics and earthquake hazards.Although regional velocity models are available for the LMS fault zone,high-resolution velocity models are lacking.Therefore,a dense array of 240 short-period seismometers was deployed around the central segment of the LMS fault zone for approximately 30 days to monitor earthquakes and characterize fine structures of the fault zone.Considering the large quantity of observed seismic data,the data processing workflow consisted of deep learning-based automatic earthquake detection,phase arrival picking,and association.Compared with the earthquake catalog released by the China Earthquake Administration,many more earthquakes were detected by the dense array.Double-difference seismic tomography was adopted to determine V_(p),V_(s),and V_(p)/V_(s)models as well as earthquake locations.The checkerboard test showed that the velocity models have spatial resolutions of approximately 5 km in the horizontal directions and 2 km at depth.To the west of the Yingxiu–Beichuan Fault(YBF),the Precambrian Pengguan complex,where most of earthquakes occurred,is characterized by high velocity and low V_(p)/V_(s)values.In comparison,to the east of the YBF,the Upper Paleozoic to Jurassic sediments,where few earthquakes occurred,show low velocity and high V_(p)/V_(s)values.Our results suggest that the earthquake activity in the LMS fault zone is controlled by the strength of the rock compositions.When the high-resolution velocity models were combined with the relocated earthquakes,we were also able to delineate the fault geometry for different faults in the LMS fault zone.
基金supported by National Natural Science Foundation of China(41504109,41404099)the Natural Science Foundation of Shandong Province(BS2015HZ008)the project of "Distinguished Professor of Jiangsu Province"
文摘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 European Commission,German Federal Ministry of Education and Research,German Federal Ministry of Economics and Technology as well as Research Institute and Industry are gratefully acknowledged for funding and supporting CO2 Storage by Injection into a Natural Storage Site CO2SINK(Project No.502599)
文摘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.
基金This work was funded by National Natural Science Foundation of China （No. 40474044）.
文摘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.
文摘： 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.
基金supported by the National Natural ScienceFoundation of China (No. 40804008)
文摘The perfectly matched layer（PML） was first introduced by Berenger as an absorbing boundary condition for electromagnetic wave propagation.In this article,a method is developed to ex-tend the PML to simulating seismic wave propagation in fluid-saturated porous medium.This non-physical boundary is used at the computational edge of a Forsyte polynomial convolutional differenti-ator（FPCD） algorithm as an absorbing boundary condition to truncate unbounded media.The incor-poration of PML in Biot＇s equations is given.Numerical results show that the PML absorbing bound-ary condition attenuates the outgoing waves effectively and eliminates the reflections adequately.
文摘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.
基金Supported by the National Natural Science Foundation of China(Nos. 41206043, 40930845)the Open Foundation of Key Laboratory of Marine Geology and Environment of Chinese Academy of Sciences(No. MGE2011KG07)+1 种基金the Knowledge Innovation Program of Chinese Academy of Sciences (No. KZCX2-YW-229)the National Basic Research Program of China (973 Program) (No. 2009CB219505)
文摘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.
基金supported by the National Nature Science Foundation of China(No.41474112)the National Science and Technology Major Project(No.2017ZX05005-004)
文摘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%.
文摘Fractured reservoirs have always been a big favorable area for oil and gas reservoirs,so prediction of fractures is also a research hotspot in recent years.Due to the diversity of fracture development and the unclear development mechanism,fracture prediction has always been a major problem.Simple numerical simulation In this paper,seismic attribute is combined with numerical simulation,logging data and actual seismic profile are used as constraints,inversion impedance value and coherent attribute are combined,and finally a property model more in line with the actual geological conditions is established.The wave equation calculation and migration processing were used to obtain the numerical simulation profile,and the actual seismic profile,fracture detection profile and numerical simulation profile were combined for analysis:①The numerical simulation section under this modeling method can greatly correspond to the actual seismic section,and the reflected results can better reflect the changes of response characteristics.②The reliability and applicability of the fracture detection technology can be determined by comparing the forward simulation profile with the fracture detection profile.
基金National Key Research and Development Program of China(No.2022YFC3004300)the National Natural Science Foundation of China(No.52378475).
文摘To study the ground motion intensity measures(IMs)suitable for the design of seismic performance with a focus on longitudinal resistance in tunnel structures,21 different seismic intensity parameters are selected for nonlinear calculation and analysis of tunnel structures,in order to determine the optimal IM for the longitudinal seismic performance of tunnel structures under different site conditions.An improved nonlinear beam-spring model is developed to calculate the longitudinal seismic response of tunnels.The PQ-Fiber model is used to simulate the longitudinal nonlinear behavior of tunnel structures and the tangential interactions between the tunnel and the soil is realized by load in the form of moment.Five different site types are considered and 21 IMs is evaluated against four criteria:effectiveness,practicality,usefulness,and sufficiency.The results indicate that the optimal IMs are significantly influenced by the site conditions.Specifically,sustained maximum velocity(V_(SM))emerges as the optimal IM for circular tunnels in soft soil conditions(CaseⅠsites),peak ground velocity(V PG)is best suited for CaseⅡsites,sustained maximum acceleration(A_(SM))is ideal for both CaseⅢand CaseⅤsites,and peak ground acceleration(A PG)for CaseⅣsites.As site conditions transition from CaseⅠto CaseⅤ,from soft to hard,the applicability of acceleration-type intensity parameters gradually decreases,while the applicability of velocity-type intensity parameters gradually increases.
基金sponsored by National Science and Technology Major Project(2011ZX05046-001)
文摘Large-scale 3D physical models of complex structures can be used to simulate hydrocarbon exploration areas. The high-fidelity simulation of actual structures poses challenges to model building and quality control. Such models can be used to collect wideazimuth, multi-azimuth, and full-azimuth seismic data that can be used to verify various 3D processing and interpretation methods. Faced with nonideal imaging problems owing to the extensive complex surface conditions and subsurface structures in the oil-rich foreland basins of western China, we designed and built the KS physical model based on the complex subsurface structure. This is the largest and most complex 3D physical model built to date. The physical modeling technology advancements mainly involve 1） the model design method, 2） the model casting flow, and 3） data acquisition. A 3D velocity model of the physical model was obtained for the first time, and the model building precision was quantitatively analyzed. The absolute error was less than 3 mm, which satisfies the experimental requirements. The 3D velocity model obtained from 3D measurements of the model layers is the basis for testing various imaging methods. Furthermore, the model is considered a standard in seismic physical modeling technology.
基金supported by the National Natural Science Foundation of China under grant Nos.40474012,40874020 and 40821062
文摘The pseudospectral method has been applied to the simulation of seismic wave propagation in 2-D global Earth model. When a whole Earth model is considered, the center of the Earth is included in the model and then singularity arises at the center of the Earth where r=0 since the 1/r term appears in the wave equations. In this paper, we extended the global seismic wavefield simulation algorithm for regular grid mesh to staggered grid configuration and developed a scheme to solve the numerical problems associated with the above singularity for a 2-D global Earth model defined on staggered grid using pseudospectral method. This scheme uses a coordinate transformation at the center of the model, in which the field variables at the center are calculated in Cartesian coordinates from the values on the grids around the center. It allows wave propagation through the center and hence the wavefield at the center can be stably calculated. Validity and accuracy of the scheme was tested by compared with the discrete wavenumber method. This scheme could also be suitable for other numerical methods or models parameterized in cylindrical or spherical coordinates when singularity arises at the center of the model.
基金Science Council, Chinese Taipei 106, Under Grant No. NSC-95-2625-Z-027-004
文摘Strong near-fault ground motion, usually caused by the fault-rupture and characterized by a pulse-like velocity- wave form, often causes dramatic instantaneous seismic energy （Jadhav and Jangid 2006）. Some reinforced concrete （RC） bridge columns, even those built according to ductile design principles, were damaged in the 1999 Chi-Chi earthquake. Thus, it is very important to evaluate the seismic response of a RC bridge column to improve its seismic design and prevent future damage. Nonlinear time history analysis using step-by-step integration is capable of tracing the dynamic response of a structure during the entire vibration period and is able to accommodate the pulsing wave form. However, the accuracy of the numerical results is very sensitive to the modeling of the nonlinear load-deformation relationship of the structural member. FEMA 273 and ATC-40 provide the modeling parameters for structural nonlinear analyses of RC beams and RC columns. They use three parameters to define the plastic rotation angles and a residual strength ratio to describe the nonlinear load- deformation relationship of an RC member. Structural nonlinear analyses are performed based on these parameters. This method provides a convenient way to obtain the nonlinear seismic responses of RC structures. However, the accuracy of the numerical solutions might be further improved. For this purpose, results from a previous study on modeling of the static pushover analyses for RC bridge columns （Sung et al. 2005） is adopted for the nonlinear time history analysis presented herein to evaluate the structural responses excited by a near-fault ground motion. To ensure the reliability of this approach, the numerical results were compared to experimental results. The results confirm that the proposed approach is valid.
基金funded by the National Natural Science Foundation of China under Grant 42230101the National Agency for Research and Development of Chile(ANID)by Projects AFB180004 and AFB220002the ANID Programa de Cooperación Internacional(PCI)Grant PII-180003.
文摘The Chilean Pampean flat slab subduction segment is characterized by the nearly horizontal subduction of the Nazca Plate within the depth range of 100-120 km.Numerous seismic tomography studies have been conducted to investigate its velocity structure;however,they have used only seismic body wave data or surface wave data.As a result,the existing velocity models in the region may have relatively large uncertainties.In this study,we use body wave arrival times from earthquakes occurring in central Chile between 2014 and 2019,as well as Rayleigh wave phase velocity maps at periods of 5-80 s from ambient noise empirical Green’s functions in Chile.By jointly using body wave arrival times and surface wave dispersion data,we refine the VS model and improve earthquake locations in the central Chile subduction zone.Compared with previous velocity models,our velocity model better reveals an eastward-dipping high-velocity plate representing the subducting Nazca Plate,which is 40-50 km thick and is more consistent with the slab thickness estimated by receiver function imaging and thermal modeling.Overall,the intraslab seismicity distribution spatially correlates well with the slab high-velocity anomalies except along the subduction paths of the CopiapóRidge and Juan Fernández Ridge.Additionally,parallel low-velocity stripes are imaged beneath the subducting plate,which are likely associated with the accumulated melts.The joint inversion velocity model also resolves widespread low-velocity anomalies in the crust beneath the Central Volcanic Zone of the central Andes,likely representing crustal magma chambers for various volcanoes.
