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.

Accurate simulations of pure-viscoacoustic wave propagation in tilted transversely isotropic media

Forward modeling of seismic wave propagation is crucial for the realization of reverse time migration(RTM) and full waveform inversion(FWI) in attenuating transversely isotropic media. To describe the attenuation and anisotropy properties of subsurface media, the pure-viscoacoustic anisotropic wave equations are established for wavefield simulations, because they can provide clear and stable wavefields. However, due to the use of several approximations in deriving the wave equation and the introduction of a fractional Laplacian approximation in solving the derived equation, the wavefields simulated by the previous pure-viscoacoustic tilted transversely isotropic(TTI) wave equations has low accuracy. To accurately simulate wavefields in media with velocity anisotropy and attenuation anisotropy, we first derive a new pure-viscoacoustic TTI wave equation from the exact complex-valued dispersion formula in viscoelastic vertical transversely isotropic(VTI) media. Then, we present the hybrid finite-difference and low-rank decomposition(HFDLRD) method to accurately solve our proposed pure-viscoacoustic TTI wave equation. Theoretical analysis and numerical examples suggest that our pure-viscoacoustic TTI wave equation has higher accuracy than previous pure-viscoacoustic TTI wave equations in describing q P-wave kinematic and attenuation characteristics. Additionally, the numerical experiment in a simple two-layer model shows that the HFDLRD technique outperforms the hybrid finite-difference and pseudo-spectral(HFDPS) method in terms of accuracy of wavefield modeling.

Evaluation of 3D crustal seismic velocity models in southwest China:Model performance,limitation,and prospects

Southwest China is a tectonically and seismically active region,witnessing strong deformation due to the collision between the Indian and Eurasian plates.Constraining the subsurface velocity structure of this region is thus important in understanding the tectonics and geodynamic processes of continental collision and in mitigating seismic hazards.Numerous studies have provided various 3D seismic velocity models in southwest China.However,discrepancies exist among these models,and less effort has been made to quantify the reliability and accuracy of these existing velocity models.In this study,we use regional 3D waveform simulation to evaluate the performance of various regional crustal 3D velocity models in reproducing observed seismograms.We particularly focus on two recent earthquake sequence in the region,the 2021 Yunnan Yangbi MS_(6.4) earthquake sequence and the 2022 Sichuan Luding MS_(6.8) earthquake sequence.The tested 3D velocity models include the Southwest China Community Velocity model V1.0,the Unified Seismic Tomography Models for Continental China Lithosphere V2.0,the adjoint full waveform tomography model of the crustal and upper mantle beneath Eastern Tibetan Plateau,and the shallow seismic structure model beneath continental China.Our results show that the tested 3D velocity models generally capture well long-period(<0.2 Hz)waveforms,indicating that the 3D models adequately resolve overall large-scale subsurface structures.However,the 3D synthetics show discrepancies in higher frequencies(0.05–0.3 Hz)and the performance of the 3D velocity models varies from region to region,suggesting that smaller scale heterogeneities are not well constrained.Including shallow velocity structures(<10 km)can improve the waveform fitting,emphasizing the importance of incorporating shallow structures in waveform modeling.The full-waveform tomography model shows a slighter better performance than the other models,especially for the body-waves,highlighting the advantages of full-waveform method in achieving sub-wavelength resolution despite the usage of very long-period waveforms.In light of these comparison results of model performance,we identify the advantages and limitations of different seismic tomography models and methods,and we propose to incorporate different tomography methods and datasets to better constrain subsurface structures.While our target region in this study is southwest China,the analysis that we have conducted can be applied to other regions of various scales and tectonic settings for quantitative seismic model evaluation.

Seismic fragility assessment of RC frame structure designed according to modern Chinese code for seismic design of buildings

Following several damaging earthquakes in China, research has been devoted to find the causes of the collapse of reinforced concrete （RC） building sand studying the vulnerability of existing buildings. The Chinese Code for Seismic Design of Buildings （CCSDB） has evolved over time, however, there is still reported earthquake induced damage of newly designed RC buildings. Thus, to investigate modern Chinese seismic design code, three low-, mid- and high-rise RC frames were designed according to the 2010 CCSDB and the corresponding vulnerability curves were derived by computing a probabilistic seismic demand model （PSDM）.The PSDM was computed by carrying out nonlinear time history analysis using thirty ground motions obtained from the Pacific Earthquake Engineering Research Center. Finally, the PSDM was used to generate fragility curves for immediate occupancy, significant damage, and collapse prevention damage levels. Results of the vulnerability assessment indicate that the seismic demands on the three different frames designed according to the 2010 CCSDB meet the seismic requirements and are almost in the same safety level.

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.

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.

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.

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.

A scheme to treat the singularity in global seismic wavefield simulation using pseudospectral method with staggered grids

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.

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.

Analytical investigations of seismic responses for reinforced concrete bridge columns subjected to strong near-fault ground motion

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 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.

3-D seismic interpretation of stratigraphic and structural features in the Upper Jurassic to Lower Cretaceous sequence of the Gullfaks Field,Norwegian North Sea:A case study of reservoir development

The 3-D seismic dataset is a key tool to analyze and understand the mechanism of structural and stratigraphic hydrocarbon(HC)trapping in the subsurface.Conventionally used subsurface seismic characterization methods for fractures are based on the theory of effective anisotropy medium.The aim of this work is to improve the structural images with dense sampling of 3-D survey to evaluate structural and stratigraphic models for reservoir development to predict reservoir quality.The present study of the Gullfaks Field,located in the Norwegian North Sea Gullfaks sector,identifies the shallowest structural elements.The steepness of westward structural dip decreases eastward during the Upper Jurassic to Lower Cretaceous deposition.Reservoir sands consist of the Middle Jurassic deltaic deposits and Lower Jurassic fluvial channel and delta plain deposits.Sediment supply steadily prevails on sea-level rise and the succession displays a regressive trend indicated by a good continuous stacking pattern.The key factor for the development of reservoirs in the Gullfaks Field is fault transmissibility with spatially distributed pressure.The majority of mapped faults with sand-to-sand contacts are non-sealing,which provide restriction for the HC flow between the fault blocks.The traps for HC accumulation occur between the post-rift and syn-rift strata,i.e.antiform set by extensional system,unconformity trap at the top of syndeposition,and structural trap due to normal faults.Overall reservoir quality in the studied area is generally excellent with average 35%porosity and permeability in the Darcy range.Our findings are useful to better understand the development of siliciclastic reservoirs in similar geological settings worldwide.

ComprehensiveanalysesofseismicsourcelayerinXingtaiandTangshanseismicregionsandtheconditionsofmediaab

Statistic analyses of sectional focal depth frequency relation for 1739 M S≥4.0 events in North China, 541 M S≥4.0 events in Xingtai seismic region and 999 M S≥4.0 events in Tangshan seismic region have been carried out for the first step of the study. Then the sectional focal depth frequency diagrams (H N diagram) for each region were obtained. The study shows that: (1) The intraplate earthquakes in North China occur within the continental crust, and are mainly concentrated within a certain depth range of the middle and upper crust, making up a seismic source layer; (2) The seismic source layer consists of a main layer and an auxiliary layer adjacent to it. The tops and bottoms of the main and auxiliary layers correspond to the G and C, as well as to B and G geologic boundaries respectively; (3) The H N diagrams for Xingtai and Tangshan seismic regions are significantly different from each other, indicating the regionality and localization of the seismic source layers; (4) DSS and MTS observations show that the regionality and localization of the seismic source layer depend mainly on the environment of the deep structures of the region; (5) The seismic source layer is a substance of medium layer, which can be related to the brittle ductile transition zone within the crust. The seismic source layer is also a stress conducting layer. As the media above and below this layer are subjected to failure or flowage, then the stress will concentrate in the stress conducting layer, in which the earthquake occurs.

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

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

Forward modeling of fracture prediction based on seismic attribute modeling

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.

A drilling technology guided by well-seismic information integration

The predictions by drilling-related mechanical and geological models are in some degree inaccurate due to non-unique solution of seismic velocity model.To address this problem,a new drilling technology guided by well-seismic information integration is proposed which consists of seismic velocity update of drilled formations,seismic velocity prediction of the formation ahead of drilling bit,and the prediction of geological feature and drilling geological environmental factors ahead of bit.In this technology,real information(velocity,formation and depth)behind the drilling bit and local pre-stack seismic data around the wellbore being drilled are used to correct the primitive seismic velocity field for a re-migration of seismic data and to update geological features and drilling geological environmental factors ahead of the drilling bit.Field application shows that this technology can describe and predict the geological features,drilling geological environmental factors and complex drilling problems ahead of the bit timely and improve the prediction efficiency and accuracy greatly.These new updated results are able to provide scientific basis for optimizing drilling decisions.

Study on seismic numerical model on the Chinese mainland and its vicinity

This paper presents the stress incremental rate of each element by using finite element numerical simulation way.Combined with the practical continent boundary condition and medium parameters, this paper applies the stress increment rate to the more practical modified cell automaton model according to the cell automaton principle. The paper preliminarily gives the more practical seismic space-time repeated evolutive maps, and compares the imitative seismic energy release with the real seismic energy release. It is thought that this model has its advantageover the old one whether in the seismic imitative characteristics or in the comparison with the realsituation and is more practical model.

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.