The AVO Effect of Formation Pressure on Time-Lapse Seismic Monitoring in Marine Carbon Dioxide Storage

The phase change of CO_(2) has a significant bearing on the siting, injection, and monitoring of storage. The phase state of CO_(2) is closely related to pressure. In the process of seismic exploration, the information of formation pressure can be response in the seismic data. Therefore, it is possible to monitor the formation pressure using time-lapse seismic method. Apart from formation pressure, the information of porosity and CO_(2) saturation can be reflected in the seismic data. Here, based on the actual situation of the work area, a rockphysical model is proposed to address the feasibility of time-lapse seismic monitoring during CO_(2) storage in the anisotropic formation. The model takes into account the formation pressure, variety minerals composition, fracture, fluid inhomogeneous distribution, and anisotropy caused by horizontal layering of rock layers(or oriented alignment of minerals). From the proposed rockphysical model and the well-logging, cores and geological data at the target layer, the variation of P-wave and S-wave velocity with formation pressure after CO_(2) injection is calculated. And so are the effects of porosity and CO_(2) saturation. Finally, from anisotropic exact reflection coefficient equation, the reflection coefficients under different formation pressures are calculated. It is proved that the reflection coefficient varies with pressure. Compared with CO_(2) saturation, the pressure has a greater effect on the reflection coefficient. Through the convolution model, the seismic record is calculated. The seismic record shows the difference with different formation pressure. At present, in the marine CO_(2) sequestration monitoring domain, there is no study involving the effect of formation pressure changes on seismic records in seafloor anisotropic formation. This study can provide a basis for the inversion of reservoir parameters in anisotropic seafloor CO_(2) reservoirs.

Hcable for Time-Lapse Seismic Monitoring of Marine Carbon Capture and Storage

To ensure project safety and secure public support, an integrated and comprehensive monitoring program is needed within a carbon capture and storage(CCS) project. Monitoring can be done using many well-established techniques from various fields, and the seismic method proves to be the crucial one. This method is widely used to determine the CO_(2) distribution, image the plume development, and quantitatively estimate the concentration. Because both the CO_(2) distribution and the potential migration pathway can be spatially small scale, high resolution for seismic imaging is demanded. However, obtaining a high-resolution image of a subsurface structure in marine settings is difficult. Herein, we introduce the novel Hcable(Harrow-like cable system) technique, which may be applied to offshore CCS monitoring. This technique uses a highfrequency source(the dominant frequency>100 Hz) to generate seismic waves and a combination of a long cable and several short streamers to receive seismic waves. Ultrahigh-frequency seismic images are achieved through the processing of Hcable seismic data. Hcable is then applied in a case study to demonstrate its detailed characterization for small-scale structures. This work reveals that Hcable is a promising tool for timelapse seismic monitoring of oceanic CCS.

Monitoring the change in horizontal stress with multi-wave time-lapse seismic response based on nonlinear elasticity theory

Monitoring the change in horizontal stress from the geophysical data is a tough challenge, and it has a crucial impact on broad practical scenarios which involve reservoir exploration and development, carbon dioxide (CO_(2)) injection and storage, shallow surface prospecting and deep-earth structure description. The change in in-situ stress induced by hydrocarbon production and localized tectonic movements causes the changes in rock mechanic properties (e.g. wave velocities, density and anisotropy) and further causes the changes in seismic amplitudes, phases and travel times. In this study, the nonlinear elasticity theory that regards the rock skeleton (solid phase) and pore fluid as an effective whole is used to characterize the effect of horizontal principal stress on rock overall elastic properties and the stress-dependent anisotropy parameters are therefore formulated. Then the approximate P-wave, SV-wave and SH-wave angle-dependent reflection coefficient equations for the horizontal-stress-induced anisotropic media are proposed. It is shown that, on the different reflectors, the stress-induced relative changes in reflectivities (i.e., relative difference) of elastic parameters (i.e., P- and S-wave velocities and density) are much less than the changes in contrasts of anisotropy parameters. Therefore, the effects of stress change on the reflectivities of three elastic parameters are reasonably neglected to further propose an AVO inversion approach incorporating P-, SH- and SV-wave information to estimate the change in horizontal principal stress from the corresponding time-lapse seismic data. Compared with the existing methods, our method eliminates the need for man-made rock-physical or fitting parameters, providing more stable predictive power. 1D test illustrates that the estimated result from time-lapse P-wave reflection data shows the most reasonable agreement with the real model, while the estimated result from SH-wave reflection data shows the largest bias. 2D test illustrates the feasibility of the proposed inversion method for estimating the change in horizontal stress from P-wave time-lapse seismic data.

Purposeless repeated acquisition time-lapse seismic data processing

In China, most oil fields are continental sedimentation with strong heterogeneity, which on one side makes the reservoir prospecting and development more difficult, but on the other side provides more space for searching residual oil in matured fields. Time-lapse seismic reservoir monitoring technique is one of most important techniques to define residual oil distribution. According to the demand for and development of time-lapse seismic reservoir monitoring in China, purposeless repeated acquisition time-lapse seismic data processing was studied. The four key steps in purposeless repeated acquisition time-lapse seismic data processing, including amplitude-preserved processing with relative consistency, rebinning, match filtering and difference calculation, were analyzed by combining theory and real seismic data processing. Meanwhile, quality control during real time-lapse seismic processing was emphasized.

Quantitative prediction of oil saturation of unconsolidated sandstone reservoir based on time-lapse seismic “relative difference method”: Taking Zeta oil field in West Africa as an example

In view of the disadvantage that the absolute difference of time-lapse seismic(the difference between monitoring data and base data) is not only related to the change of oil saturation, but also closely related to the thickness of reservoir, a time-lapse seismic "relative difference method"(the ratio of monitoring data to base data) not affected by the thickness of reservoir but only related to the change of fluid saturation, is proposed through seismic forward modeling after fluid displacement simulation. Given the same change of fluid saturation, the absolute difference of time-lapse seismic conforms to the law of "tuning effect" and seismic reflection of "thin bed", and the remaining oil prediction method based on absolute difference of time-lapse seismic is only applicable to the reservoirs with uniform thickness smaller than the tuning thickness or with thickness greater than the tuning thickness. The relative difference of time-lapse seismic is not affected by reservoir thickness, but only related to the change of fluid saturation. It is applicable to all the deep-sea unconsolidated sandstone reservoirs which can exclude the effect of pressure, temperature, pore type and porosity on seismic. Therefore, the relation between the relative difference of time-lapse seismic and the change of fluid saturation, which is obtained from seismic forward modeling after Gassmann fluid displacement simulation, can be used to quantitatively predict the change of reservoir water saturation and then the distribution of the remaining oil. The application of this method in deep sea Zeta oil field in west Africa shows that it is reasonable and effective.

A NEW INVERSION METHOD OF TIME-LAPSE SEISMIC

Time-Lapse Seismic improves oil recovery ratio by dynamic reservoir monitoring. Because of the large number of seismic explorations in the process of time-lapse seismic inversion, traditional methods need plenty of inversion calculations which cost high computational works. The method is therefore inefficient. In this paper, in order to reduce the repeating computations in traditional, a new time-lapse seismic inversion method is put forward. Firstly a homotopy-regularization method is proposed for the first time inversion. Secondly, with the first time inversion results as the initial value of following model, a model of the second time inversion is rebuilt by analyzing the characters of time-lapse seismic and localized inversion method is designed by using the model. Finally, through simulation, the comparison between traditional method and the new scheme is given. Our simulation results show that the new scheme could save the algorithm computations greatly.

Time-lapse培养与桌面台式培养箱胚胎培养对胚胎发育的影响

目的探讨Time-lapse培养与桌面台式培养箱胚胎培养对胚胎发育的影响。方法选取2021年3月—2023年3月在本院生殖中心行IVF-ET/ICSI-ET治疗的613个周期的胚胎发育数据,随机分配胚胎培养方式。分为对照组与研究组,其中对照组为桌面台式培养箱方式(360个周期),研究组为Time-lapse培养与桌面台式培养箱方式(253个周期),统计两组胚胎发育质量。结果研究组中正常受精数、可移植胚胎数、形成囊胚数、优质胚胎数、高评分囊胚数、受精率、冷冻胚胎率、胚囊形成率、优质胚胎率、高评分囊胚率高于对照组,组间差异有统计学意义(P<0.05)。结论Time-lapse培养与桌面台式培养箱胚胎培养方式优于桌面台式培养箱方式,可为胚胎培养提供适宜的环境、温度、湿度等,有助于胚胎发育,具有较强的实用性和推广性,值得推广应用。

Target-Oriented Inversion of Time-Lapse Seismic Waveform Data

Full waveform inversion of time-lapse seismic data can be used as a means of estimating the reservoir changes due to the production.Since the repeated computa-tions for the monitor surveys lead to a large computational cost,time-lapse full wave-form inversion is still considered to be a challenging task.To address this problem,we present an efficient target-oriented inversion scheme for time-lapse seismic data using an integral equation formulation with Gaussian beam based Green’s function approach.The proposed time-lapse approach allows one to perform a local inversion within a small region of interest(e.g.a reservoir under production)for the monitor survey.We have verified that the T-matrix approach is indeed naturally target-oriented,which was mentioned by Jakobsen and Ursin[24]and allows one to reduce the compu-tational cost of time-lapse inversion by focusing the inversion on the target-area only.This method is based on a new version of the distorted Born iterative T-matrix inverse scattering method.The Gaussian beam and T-matrix are used in this approach to perform the wavefield computation for the time-lapse inversion in the baseline model from the survey surface to the target region.We have provided target-oriented inversion results of the synthetic time-lapse waveform data,which shows that the proposed scheme reduces the computational cost significantly.

Numerical Modeling of Anisotropic Elastic-Wave Sensitivity Propagation for Optimal Design of Time-Lapse Seismic Surveys

Reliable subsurface time-lapse seismic monitoring is crucial for many geo-physical applications,such as enhanced geothermal system characterization,geologic carbon utilization and storage,and conventional and unconventional oil/gas reservoir characterization,etc.We develop an elastic-wave sensitivity propagation method for optimal design of cost-effective time-lapse seismic surveys considering the fact that most of subsurface geologic layers and fractured reservoirs are anisotropic instead of isotropic.For anisotropic media,we define monitoring criteria using qP-and qS-wave sensitivity energies after decomposing qP-and qS-wave components from the total elastic-wave sensitivity wavefield using a hybrid time-and frequency-domain approach.Geophones should therefore be placed at locations with significant qP-and qS-wave sensitivity energies for cost-effective time-lapse seismic monitoring in an anisotropic geology setting.Our numerical modeling results for a modified anisotropic Hess model demonstrate that,compared with the isotropic case,subsurface anisotropy changes the spatial distributions of elastic-wave sensitivity energies.Consequently,it is necessary to consider subsurface anisotropies when designing the spatial distri-bution of geophones for cost-effective time-lapse seismic monitoring.This finding suggests that it is essential to use our new anisotropic elastic-wave sensitivity modeling method for optimal design of time-lapse seismic surveys to reliably monitor the changes in subsurface reservoirs,fracture zones or target monitoring regions.

Introduction to a recently released dataset entitled CSNCD:A Comprehensive Dataset of Chinese Seismic Network

China and its adjacent regions are heavily affected by earthquake disasters,which has led to strong demand and vigorous development of the Chinese Seismic Network(CSN).Since 2004,China Earthquake Networks Center(CENC)has been authorized by China Earthquake Administration(CEA)to take charge of CSN management,mainly observing product output and quality control.

Probabilistic analysis of tunnel face seismic stability in layered rock masses using Polynomial Chaos Kriging metamodel

Face stability is an essential issue in tunnel design and construction.Layered rock masses are typical and ubiquitous;uncertainties in rock properties always exist.In view of this,a comprehensive method,which combines the Upper bound Limit analysis of Tunnel face stability,the Polynomial Chaos Kriging,the Monte-Carlo Simulation and Analysis of Covariance method(ULT-PCK-MA),is proposed to investigate the seismic stability of tunnel faces.A two-dimensional analytical model of ULT is developed to evaluate the virtual support force based on the upper bound limit analysis.An efficient probabilistic analysis method PCK-MA based on the adaptive Polynomial Chaos Kriging metamodel is then implemented to investigate the parameter uncertainty effects.Ten input parameters,including geological strength indices,uniaxial compressive strengths and constants for three rock formations,and the horizontal seismic coefficients,are treated as random variables.The effects of these parameter uncertainties on the failure probability and sensitivity indices are discussed.In addition,the effects of weak layer position,the middle layer thickness and quality,the tunnel diameter,the parameters correlation,and the seismic loadings are investigated,respectively.The results show that the layer distributions significantly influence the tunnel face probabilistic stability,particularly when the weak rock is present in the bottom layer.The efficiency of the proposed ULT-PCK-MA is validated,which is expected to facilitate the engineering design and construction.

An adaptive physics-informed deep learning method for pore pressure prediction using seismic data

Accurate prediction of formation pore pressure is essential to predict fluid flow and manage hydrocarbon production in petroleum engineering.Recent deep learning technique has been receiving more interest due to the great potential to deal with pore pressure prediction.However,most of the traditional deep learning models are less efficient to address generalization problems.To fill this technical gap,in this work,we developed a new adaptive physics-informed deep learning model with high generalization capability to predict pore pressure values directly from seismic data.Specifically,the new model,named CGP-NN,consists of a novel parametric features extraction approach(1DCPP),a stacked multilayer gated recurrent model(multilayer GRU),and an adaptive physics-informed loss function.Through machine training,the developed model can automatically select the optimal physical model to constrain the results for each pore pressure prediction.The CGP-NN model has the best generalization when the physicsrelated metricλ=0.5.A hybrid approach combining Eaton and Bowers methods is also proposed to build machine-learnable labels for solving the problem of few labels.To validate the developed model and methodology,a case study on a complex reservoir in Tarim Basin was further performed to demonstrate the high accuracy on the pore pressure prediction of new wells along with the strong generalization ability.The adaptive physics-informed deep learning approach presented here has potential application in the prediction of pore pressures coupled with multiple genesis mechanisms using seismic data.

Seismic performance evaluation of water supply pipes installed in a full-scale RC frame structure based on a shaking table test

As an important part of nonstructural components,the seismic response of indoor water supply pipes deserves much attention.This paper presents shaking table test research on water supply pipes installed in a full-scale reinforced concrete(RC)frame structure.Different material pipes and different methods for penetrating the reinforced concrete floors are combined to evaluate the difference in seismic performance.Floor response spectra and pipe acceleration amplification factors based on test data are discussed and compared with code provisions.A seismic fragility study of displacement demand is conducted based on numerical simulation.The acceleration response and displacement response of different combinations are compared.The results show that the combination of different pipe materials and different passing-through methods can cause obvious differences in the seismic response of indoor riser pipes.

Time-domain dynamic constitutive model suitable for mucky soil site seismic response

Soil nonlinear behavior displays noticeable effects on the site seismic response.This study proposes a new functional expression of the skeleton curve to replace the hyperbolic skeleton curve.By integrating shear modulus and combining the dynamic skeleton curve and the damping degradation coefficient,the constitutive equation of the logarithmic dynamic skeleton can be obtained,which considers the damping effect in a soil dynamics problem.Based on the finite difference method and the multi-transmitting boundary condition,a 1D site seismic response analysis program called Soilresp1D has been developed herein and used to analyze the time-domain seismic response in three types of sites.At the same time,this study also provides numerical simulation results based on the hyperbolic constitutive model and the equivalent linear method.The results verify the rationality of the new soil dynamic constitutive model.It can analyze the mucky soil site nonlinear seismic response,reflecting the deformation characteristics and damping effect of the silty soil.The hysteresis loop area is more extensive,and the residual strain is evident.

Monitoring seismicity in the southern Sichuan Basin using a machine learning workflow

Monitoring seismicity in real time provides significant benefits for timely earthquake warning and analyses.In this study,we propose an automatic workflow based on machine learning(ML)to monitor seismicity in the southern Sichuan Basin of China.This workflow includes coherent event detection,phase picking,and earthquake location using three-component data from a seismic network.By combining Phase Net,we develop an ML-based earthquake location model called Phase Loc,to conduct real-time monitoring of the local seismicity.The approach allows us to use synthetic samples covering the entire study area to train Phase Loc,addressing the problems of insufficient data samples,imbalanced data distribution,and unreliable labels when training with observed data.We apply the trained model to observed data recorded in the southern Sichuan Basin,China,between September 2018 and March 2019.The results show that the average differences in latitude,longitude,and depth are 5.7 km,6.1 km,and 2 km,respectively,compared to the reference catalog.Phase Loc combines all available phase information to make fast and reliable predictions,even if only a few phases are detected and picked.The proposed workflow may help real-time seismic monitoring in other regions as well.

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.

Dynamically triggered seismicity on a tectonic scale:A review

Earthquakes triggered by dynamic disturbances have been confirmed by numerous observations and experiments.In the past several decades,earthquake triggering has attracted increasing attention of scholars in relation to exploring the mechanism of earthquake triggering,earthquake prediction,and the desire to use the mechanism of earthquake triggering to reduce,prevent,or trigger earthquakes.Natural earthquakes and large‐scale explosions are the most common sources of dynamic disturbances that trigger earthquakes.In the past several decades,some models have been developed,including static,dynamic,quasi‐static,and other models.Some reviews have been published,but explosiontriggered seismicity was not included.In recent years,some new results on earthquake triggering have emerged.Therefore,this paper presents a new review to reflect the new results and include the content of explosion‐triggered earthquakes for the reference of scholars in this area.Instead of a complete review of the relevant literature,this paper primarily focuses on the main aspects of dynamic earthquake triggering on a tectonic scale and makes some suggestions on issues that need to be resolved in this area in the future.

Potential of the Arkhangelsk seismic network for European Arctic monitoring

The Arkhangelsk Seismic Network(ASN)of the N.Laverov Federal Center for Integrated Arctic Research of the Ural Branch of the Russian Academy of Sciences,founded in 2003,includes 10 permanent seismic stations located on the coasts of the White,Barents,and Kara Seas and on the Arctic archipelagos of Novaya Zemlya,Franz Josef Land,and Severnaya Zemlya.The network is registered with the International Federation of Digital Seismograph Networks and the International Seismological Center.We used not only ASN data to process earthquakes but also the waveforms of various international seismic stations.The 13,000 seismic events were registered using ASN data for 2012-2022,and for 5,500 of them,we determined the parameters of the earthquake epicenters from the European Arctic.The spatial distribution of epicenters shows that the ASN monitors not only the main seismically active zones but also weak seismicity on the shelf of the Barents and Kara Seas.The representative magnitude of ASN was ML,rep=3.5.The level of microseismic noise has seasonal variations that affect the registration capabilities of each station included in the ASN and the overall sensitivity of the network as a whole.In summer,the sensitivity of the ASN decreased owing to the increasing microseismic and ambient noises,whereas in winter,the sensitivity of the ASN increased significantly because of the decrease.

The 2023 Turkey earthquake doublet: Earthquake relocation, seismic tomography, and stress field inversion

On February 6,2023,two earthquakes with magnitudes of M_(W) 7.8 and M_(W) 7.5 struck southeastern Turkey,causing significant casualties and economic losses.These seismic events occurred along the East Anatolian Fault Zone,a convergent boundary between the Arabian Plate and the Anatolian Subplate.In this study,we analyze the M_(W) 7.8 and M_(W) 7.5 earthquakes by comparing their aftershock relocations,tomographic images,and stress field inversions.The earthquakes were localized in the upper crust and exhibited steep dip angles.Furthermore,the aftershocks occurred either close to the boundaries of low and high P-wave velocity anomaly zones or within the low P-wave velocity anomaly zones.The East Anatolia Fault,associated with the M_(W) 7.8 earthquake,and the SürgüFault,related to the M_(W) 7.5 earthquake,predominantly experienced shear stress.However,their western sections experienced a combination of strike-slip and tensile stresses in addition to shear stress.The ruptures of the M_(W) 7.8 and M_(W) 7.5 earthquakes appear to have bridged a seismic gap that had seen sparse seismicity over the past 200 years prior to the 2023 Turkey earthquake sequence.