Joint inversion of Rayleigh group and phase velocities for S-wave velocity structure of the 2021 M_(S)6.0 Luxian earthquake source area,China

On September 16,2021,a MS6.0 earthquake struck Luxian County,one of the shale gas blocks in the Southeastern Sichuan Basin,China.To understand the seismogenic environment and its mechanism,we inverted a fine three-dimensional S-wave velocity model from ambient noise tomography using data from a newly deployed dense seismic array around the epicenter,by extracting and jointly inverting the Rayleigh phase and group velocities in the period of 1.6–7.2 s.The results showed that the velocity model varied significantly beneath different geological units.The Yujiasi syncline is characterized by low velocity at depths of~3.0–4.0 km,corresponding to the stable sedimentary layer in the Sichuan Basin.The eastern and western branches of the Huayingshan fault belt generally exhibit high velocities in the NE-SW direction,with a few local low-velocity zones.The Luxian MS6.0 earthquake epicenter is located at the boundary between the high-and low-velocity zones,and the earthquake sequences expand eastward from the epicenter at depths of 3.0–5.0 km.Integrated with the velocity variations around the epicenter,distribution of aftershock sequences,and focal mechanism solution,it is speculated that the seismogenic mechanism of the main shock might be interpreted as the reactivation of pre-existing faults by hydraulic fracturing.

Three dimensional shear wave velocity structure of crust and upper mantle in South China Sea and its adjacent regions by surface waveform inversion

We assembled approximately 328 seismic records. The data set was from 4 digitally recording long-period and broadband stations of CDSN. We carried out the inversion based on the partitioned waveform inversion (PWI). It partitions the large-scale optimization problem into a number of independent small-scale problems. We adopted surface waveform inversion with an equal block (2((2() discretization in order to acquire the images of shear velocity structure at different depths (from surface to 430 km) in the crust and upper-mantle. The resolution of all these anomalies has been established with (check-board( resolution tests. These results show significant difference in velocity, lithosphere and asthenosphere structure between South China Sea and its adjacent regions.

Surface wave inversion with unknown number of soil layers based on a hybrid learning procedure of deep learning and genetic algorithm

Surface wave inversion is a key step in the application of surface waves to soil velocity profiling.Currently,a common practice for the process of inversion is that the number of soil layers is assumed to be known before using heuristic search algorithms to compute the shear wave velocity profile or the number of soil layers is considered as an optimization variable.However,an improper selection of the number of layers may lead to an incorrect shear wave velocity profile.In this study,a deep learning and genetic algorithm hybrid learning procedure is proposed to perform the surface wave inversion without the need to assume the number of soil layers.First,a deep neural network is adapted to learn from a large number of synthetic dispersion curves for inferring the layer number.Then,the shear-wave velocity profile is determined by a genetic algorithm with the known layer number.By applying this procedure to both simulated and real-world cases,the results indicate that the proposed method is reliable and efficient for surface wave inversion.

Research on Multi-Wave Pore Pressure Prediction Method Based on Three Field Velocity Fusion

The optimization of velocity field is the core issue in reservoir seismic pressure prediction. For a long time, the seismic processing velocity analysis method has been used in the establishment of pressure prediction velocity field, which has a long research period and low resolution and restricts the accuracy of seismic pressure prediction;This paper proposed for the first time the use of machine learning algorithms, based on the feasibility analysis of wellbore logging pressure prediction, to integrate the CVI velocity inversion field, velocity sensitive post stack attribute field, and AVO P-wave and S-wave velocity reflectivity to obtain high-precision seismic P and S wave velocities. On this basis, high-resolution formation pore pressure and other parameters prediction based on multi waves is carried out. The pressure prediction accuracy is improved by more than 50% compared to the P-wave resolution of pore pressure prediction using only root mean square velocity. Practice has proven that the research method has certain reference significance for reservoir pore pressure prediction.

Stepwise joint inversion of surface wave dispersion,Rayleigh wave ZH ratio,and receiver function data for 1D crustal shear wave velocity structure

Accurate determination of seismic velocity of the crust is important for understanding regional tectonics and crustal evolution of the Earth. We propose a stepwise joint linearized inversion method using surface wave dispersion, Rayleigh wave ZH ratio （i.e., ellipticity）, and receiver function data to better resolve 1D crustal shear wave velocity （Vs） structure. Surface wave dispersion and Rayleigh wave ZH ratio data are more sensitive to absolute variations of shear wave speed at depths, but their sensi- tivity kernels to shear wave speeds are different and complimentary. However, receiver function data are more sensitive to sharp velocity contrast （e.g., due to the existence of crustal interfaces） and Vp/Vs ratios. The stepwise inversion method takes advantages of the complementary sensitivities of each dataset to better constrain the Vs model in the crust. We firstly invert surface wave dispersion and ZH ratio data to obtain a 1D smooth absolute vs model and then incorporate receiver function data in the joint inver- sion to obtain a finer Vs model with better constraints on interface structures. Through synthetic tests, Monte Carlo error analyses, and application to real data, we demonstrate that the proposed joint inversion method can resolve robust crustal Vs structures and with little initial model dependency.

Inversion of shear wave interval velocity on prestack converted wave data

Seismic velocity is important to migration of seismic data, interpretation of lithology and lithofacies as well as prediction of reservoir. The information of shear wave velocity is required to reduce the uncertainty for discriminating lithology, identifying fluid type in porous material and calculating gas saturation in reservoir prediction. Based on Zoeppritz equations, a numeral and scanning method was proposed in this paper. Shear wave velocity can be calculated with prestack converted wave data. The effects were demonstrated by inversion of theoretical and real seismic data.

Crustal velocity structure of central Gansu Province from regional seismic waveform inversion using firework algorithm

The firework algorithm（FWA） is a novel swarm intelligence-based method recently proposed for the optimization of multi-parameter, nonlinear functions. Numerical waveform inversion experiments using a synthetic model show that the FWA performs well in both solution quality and efficiency. We apply the FWA in this study to crustal velocity structure inversion using regional seismic waveform data of central Gansu on the northeastern margin of the Qinghai-Tibet plateau. Seismograms recorded from the moment magnitude（MW） 5.4 Minxian earthquake enable obtaining an average crustal velocity model for this region. We initially carried out a series of FWA robustness tests in regional waveform inversion at the same earthquake and station positions across the study region,inverting two velocity structure models, with and without a low-velocity crustal layer; the accuracy of our average inversion results and their standard deviations reveal the advantages of the FWA for the inversion of regional seismic waveforms. We applied the FWA across our study area using three component waveform data recorded by nine broadband permanent seismic stations with epicentral distances ranging between 146 and 437 km. These inversion results show that the average thickness of the crust in this region is 46.75 km, while thicknesses of the sedimentary layer, and the upper, middle, and lower crust are 3.15,15.69, 13.08, and 14.83 km, respectively. Results also show that the P-wave velocities of these layers and the upper mantle are 4.47, 6.07, 6.12, 6.87, and 8.18 km/s,respectively.

Inversion Method for Sound Velocity Profile of Eddy in Deep Ocean

The modal wave number tomography approach is used to obtain sound speed profile of water column in deep ocean. The approach consists of estimation of the local modal eigenvalues from complex pressure field and use of these data as input to modal perturbative inversion method for obtaining the local sound speed profile. The empirical orthonormal function (EOF) is applied to reduce the parameter search space. The ocean environment used for numerical simulations includes the Munk profile as the unperturbed background speed profile and a weak Gaussian eddy as the sound speed profile perturbation. The results of numerical simulations show the method is capable of monitoring the oceanic interior structure.

Two-step interface and velocity inversion—— Study o e of the Tangshan seismic region

This paper studies the computation method of two step inversion of interface and velocity in a region. The 3 D interface is described by a segmented incomplete polynomial; while the reconstruction of 3 D velocity is accomplished by the principle of least squares in functional space. The computation is carried out in two steps. The first step is to inverse the shape of 3 D interface; while the second step is to do 3 D velocity inversion by distributing the remaining residual errors of travel time in accordance with their weights. The data of seismic sounding in the Tangshan Luanxian seismic region are processed, from which the 3 D structural form in depth of the Tangshan seismic region and the 3 D velocity distribution in the crust below the Tangshan Luanxian seismic region are obtained. The result shows that the deep 3 D structure in the Tangshan seismic region trends NE on the whole and the structure sandwiched between the NE trending Fengtai Yejituo fault and the NE trending Tangshan fault is an uplifted zone of the Moho. In the 3 D velocity structure of middle lower crust below that region, there is an obvious belt of low velocity anomaly to exist along the NE trending Tangshan fault, the position of which tallies with that of the Tangshan seismicity belt. The larger block of low velocity anomaly near Shaheyi corresponds to a denser earthquake distribution. In that region, there is an NW trending belt of high velocity anomaly, probably a buried fault zone. The lower crust below the epicentral region of the Tangshan M S=7.8 earthquake is a place where the NE trending belt of low velocity anomaly meets the NW trending belt of high velocity anomaly. The two sets of structures had played an important role in controlling the preparation and occurrence of the M S=7.8 Tangshan earthquake.

Simultaneous Inversion of Earthquake Relocation and Velocity Structure in the Shanxi-reservoir,Wenzhou

Using the data of P-wave travel time recorded at the Shanxi-reservoir seismological network and Zhejiang and Fujian local networks, we implemented a simultaneous inversion of earthquake relocation and velocity structure and determined the new locations of earthquakes in the Shanxi-reservoir. The results show that: (1) the overall epicenter distribution is NW directed, and the Shanxi reservoir induced seismicity has a close relationship to the Shuangxi-Jiaoxiyang fault; (2) the focal depth of the Shanxi reservoir induced seismicity is 5.4km in average, less than the average focal depth in the South China earthquake zone; (3) the focal depth is shallower on the reservoir shore and deeper in the reservoir inundation area. At the beginning of the reservoir induced seismicity, the focal depth increased gradually. This may be due to the gradual penetration of water into a larger depth that induced deeper earthquakes; and (4) there is a low P-wave velocity anomaly in the study area, located at the intersection of multiple faults in the reservoir inundation area. The Shanxi reservoir induced seismicity mostly occurred in this low-velocity anomaly zone. This may be related to water penetration.

S-wave velocity structure in Tangshan earthquake region and its adjacent areas from joint inversion of receiver functions and surface wave dispersion

Using the seismic records of 83 temporary and 17 permanent broadband seismic stations deployed in Tangshan earthquake region and its adjacent areas(39°N–41.5°N,115.5°E–119.5°E),we conducted a nonlinear joint inversion of receiver functions and surface wave dispersion.We obtained some detailed information about the Tangshan earthquake region and its adjacent areas,including sedimentary thickness,Moho depth,and crustal and upper mantle S-wave velocity.Meanwhile,we also obtained the vP/vS structure along two sections across the Tangshan region.The results show that:(1)the Moho depth ranges from 30 km to 38 km,and it becomes shallower from Yanshan uplift area to North China basin;(2)the thickness of sedimentary layer ranges from 0 km to 3 km,and it thickens from Yanshan uplift region to North China basin;(3)the S-wave velocity structure shows that the velocity distribution of the upper crust has obvious correlation with the surface geological structure,while the velocity characteristics of the middle and lower crust are opposite to that of the upper crust.Compared with the upper crust,the heterogeneity of the middle and lower crust is more obvious;(4)the discontinuity of Moho on the two sides of Tangshan fault suggests that Tangshan fault cut the whole crust,and the low vS and high vP/vS beneath the Tangshan earthquake region may reflect the invasion of mantle thermal material through Tangshan fault.

A Quadratic precision generalized nonlinear global optimization migration velocity inversion method

为勘探地震学的一个重要研究话题是提供一个快精确速度模型从预先叠深度移植。瞄准如此的一个问题，我们建议二次的精确概括了非线性的全球优化迁居速度倒置。首先，我们丢弃在剩余深度和剩余速度之间有一种线性关系并且求婚的假设有使从每次重复的速度模型能尽可能快速走近真正的模特儿的二次的精确的一个速度模型修正方程。第二，我们使用概括非线性的倒置得到全球最佳的速度不安模型到所有踪迹。这个方法能帮助集中速度并且能也减少在倒置期间掉进本地最小的概率。合成数据和 Marmousi 数据例子证明我们的方法有更高的精确并且需要仅仅一些重复并且因而在复杂区域提高迁居速度分析(MVA ) 的有实行可能和精确性。

Forward prediction for tunnel geology and classification of surrounding rock based on seismic wave velocity layered tomography

Excavation under complex geological conditions requires effective and accurate geological forward-prospecting to detect the unfavorable geological structure and estimate the classification of surround-ing rock in front of the tunnel face.In this work,a forward-prediction method for tunnel geology and classification of surrounding rock is developed based on seismic wave velocity layered tomography.In particular,for the problem of strong multi-solution of wave velocity inversion caused by few ray paths in the narrow space of the tunnel,a layered inversion based on regularization is proposed.By reducing the inversion area of each iteration step and applying straight-line interface assumption,the convergence and accuracy of wave velocity inversion are effectively improved.Furthermore,a surrounding rock classification network based on autoencoder is constructed.The mapping relationship between wave velocity and classification of surrounding rock is established with density,Poisson’s ratio and elastic modulus as links.Two numerical examples with geological conditions similar to that in the field tunnel and a field case study in an urban subway tunnel verify the potential of the proposed method for practical application.

Petrophysical parameters inversion for heavy oil reservoir based on a laboratory-calibrated frequency-variant rock-physics model

Heavy oil has high density and viscosity, and exhibits viscoelasticity. Gassmann's theory is not suitable for materials saturated with viscoelastic fluids. Directly applying such model leads to unreliable results for seismic inversion of heavy oil reservoir. To describe the viscoelastic behavior of heavy oil, we modeled the elastic properties of heavy oil with varying viscosity and frequency using the Cole-Cole-Maxwell (CCM) model. Then, we used a CCoherent Potential Approximation (CPA) instead of the Gassmann equations to account for the fluid effect, by extending the single-phase fluid condition to two-phase fluid (heavy oil and water) condition, so that partial saturation of heavy oil can be considered. This rock physics model establishes the relationship between the elastic modulus of reservoir rock and viscosity, frequency and saturation. The viscosity of the heavy oil and the elastic moduli and porosity of typical reservoir rock samples were measured in laboratory, which were used for calibration of the rock physics model. The well-calibrated frequency-variant CPA model was applied to the prediction of the P- and S-wave velocities in the seismic frequency range (1–100 Hz) and the inversion of petrophysical parameters for a heavy oil reservoir. The pre-stack inversion results of elastic parameters are improved compared with those results using the CPA model in the sonic logging frequency (∼10 kHz), or conventional rock physics model such as the Xu-Payne model. In addition, the inversion of the porosity of the reservoir was conducted with the simulated annealing method, and the result fits reasonably well with the logging curve and depicts the location of the heavy oil reservoir on the time slice. The application of the laboratory-calibrated CPA model provides better results with the velocity dispersion correction, suggesting the important role of accurate frequency dependent rock physics models in the seismic prediction of heavy oil reservoirs.

A deep-learning-based approach for seismic surface-wave dispersion inversion(SfNet)with application to the Chinese mainland

Surface-wave tomography is an important and widely used method for imaging the crust and upper mantle velocity structure of the Earth.In this study,we proposed a deep learning(DL)method based on convolutional neural network(CNN),named SfNet,to derive the vS model from the Rayleigh wave phase and group velocity dispersion curves.Training a network model usually requires large amount of training datasets,which is labor-intensive and expensive to acquire.Here we relied on synthetics generated automatically from various spline-based vS models instead of directly using the existing vS models of an area to build the training dataset,which enhances the generalization of the DL method.In addition,we used a random sampling strategy of the dispersion periods in the training dataset,which alleviates the problem that the real data used must be sampled strictly according to the periods of training dataset.Tests using synthetic data demonstrate that the proposed method is much faster,and the results for the vS model are more accurate and robust than those of conventional methods.We applied our method to a dataset for the Chinese mainland and obtained a new reference velocity model of the Chinese continent(ChinaVs-DL1.0),which has smaller dispersion misfits than those from the traditional method.The high accuracy and efficiency of our DL approach makes it an important method for vS model inversions from large amounts of surface-wave dispersion data.

A two-step multi-frequency receiver function inversion method for shallow crustal S-wave velocity structure and its application across the basin-mountain range belts in Northeast China

A shallow crustal velocity structure(above 10 km depth) is essential for understanding the crustal structures and deformation and assessing the exploration prospect of natural resources, and also provides priori information for imaging deeper crustal and mantle structure. Passive-source seismic methods are cost-effective and advantageous for regional-scale imaging of shallow crustal structures compared to active-source methods. Among these passive methods, techniques utilizing receiver function waveforms and/or body-wave amplitude ratios have recently gained prominence due to their relatively high spatial resolution. However, in basin regions, reverberations caused by near-surface unconsolidated sedimentary layers often introduce strong non-uniqueness and uncertainty, limiting the applicability of such methods. To address these challenges, we propose a two-step inversion method that uses multi-frequency P-RF waveforms and P-RF horizontal-to-vertical amplitude ratios. Synthetic tests indicate that our two-step inversion method can mitigate the non-uniqueness of the inversion and enhance the stability of the results. Applying this method to teleseismic data from a linear seismic array across the sedimentary basins in Northeast China, we obtain a high-resolution image of the shallow crustal S-wave velocity structure along the array. Our results reveal significant differences between the basins and mountains. The identification of low-velocity anomalies(<2.8 km s^(-1)) at depths less than 1.0 km beneath the Erlian Basin and less than 2.5 km beneath the Songliao Basin suggests the existence of sedimentary layers. Moreover, the high-velocity anomalies(~3.4–3.8 km s^(-1)) occurring at depths greater than 7 km in the Songliao Basin may reflect mafic intrusions emplaced during the Early Cretaceous. Velocity anomaly distribution in our imaging result is consistent with the location of the major faults, uplifts, and sedimentary depressions, as well as active-source seismic results. This application further validates the effectiveness of our method in constraining the depth-dependent characteristics of the S-wave velocity in basins with unconsolidated sedimentary cover.

Study on Approach for Computer-Aided Design and Machining of General Cylindrical Cam Using Relative Velocity and Inverse Kinematics

Cylindrical Cam Mechanism which is one of the best eq uipments to accomplish an accurate motion transmission is widely used in the fie lds of industries, such as machine tool exchangers, textile machinery and automa tic transfer equipments. This paper proposes a new approach for the shape design and manufacturing of the cylindrical cam. The design approach uses the relative velocity concept and the manufacturing approach uses the inverse kinematics concept. For the shape desig n, the contact points between the cam and the follower roller are calculated bas ed on relative velocity of which the direction is on the common tangential line, and then the whole shape of cam is determined from transformation of the coordi nate system. For the manufacturing procedures, the location and the orientation of cutter path can be allocated corresponding to the designed shape data. The in tegral NC code for multi-axis CNC machining center is created using the inverse kinematics concept from the data of the location and the orientation of cutter path. As the advantages of the proposed approach, the machine tool is designed t o having an alternative size in fabricating the general cam, while the tool must be fitted to diameter size of the follower in the conventional approach. Finally, CAD/CAM program, "Cylindrical DAM", is developed on C++ lan guage. This program can perform shape design, manufacturing and kinematics simul ation, which can make integral NC code for multi-axis CNC machining center. The proposed method can be applied easily on fields of industries.

SOLVING THE VELOCITY PARAMETER OF THE 2-D WAVE INVERSE PROBLEMS WITH THE INTEGRATION-CHARACTERISTIC METHOD

For the 2-D wave inverse problems introduced from geophysical exploration, in this paper, the author presents integration-characteristic method to solve the velocity parameter, and then applies it to common shotpoint model data, in noise-free case. The accuracy is quite good.

Active source seismic imaging on near-surface granite body:case study of siting a geological disposal repository for high-level radioactive nuclear waste

In order to research whether it is suitable to set a geological disposal repository for high-level radioactive nuclear waste into one target granite body,two active source seismic profles were arranged near a small town named Tamusu,Western China.The study area is with complex surface conditions,thus the seismic exploration encountered a variettraveltimey of technical difculties such as crossing obstacles,de-noising harmful scattered waves,and building complex near-surface velocity models.In order to address those problems,techniques including cross-obstacle seismic geometry design,angle-domain harmful scattered noise removal,and an acoustic wave equation-based inversion method jointly utilizing both the and waveform of frst arrival waves were adopted.The fnal seismic images clearly exhibit the target rock’s unconformable contact boundary and its top interface beneath the sedimentary and weathered layers.On this basis,it could be confrmed that the target rock is not thin or has been transported by geological process from somewhere else,but a native and massive rock.There are a few small size fractures whose space distribution could be revealed by seismic images within the rock.The fractures should be kept away.Based on current research,it could be considered that active source seismic exploration is demanded during the sitting process of the geological disposal repository for nuclear waste.The seismic acquisition and processing techniques proposed in the present paper would ofer a good reference value for similar researches in the future.

Earthquake relocation and 3-dimensional crustal structure of P-wave velocity in cen-tral-western China

采用中国中西部地区(2l°～36°N，98°～112°E)193个地震台在1992～1999年间记录到的9988次地震的Pg和Sg震相走时的读数资料，用Roecker的SPHYPIT90程序，反演了该地区三维地壳P波速度结构，并用SPHREL3D90程序进行了地震的重新定位．反演结果揭示了中国中西部地区地震P波速度结构明显的横向不均匀性，这些不同深度上波速的横向变化多以该地区的活动断裂为分界线．可以看出活动断裂两侧存在明显的速度反差．通过重新定位，得到了6459次地震的震源参数，这些精确定位的地震震中明显沿该区活动断裂呈现条带状分布，其范围和尺度清晰地表示了这一地区地震活动与活动断裂的紧密关系．其中，82％重新精确定位的事件的震源深度在20km以内．这一结果与笔者用双差地震定位法得到的重新定位的震源深度分布相一致．