Acoustic Velocity-Based Inversion of the Physical Properties of Sea Ice in the Central Arctic Region

Studying the Arctic sea ice contributes to a comprehensive understanding of the climate system in polar regions and offers valuable insights into the interplay between polar climate change and the global climate and environment.One of the key research aspects is the investigation of the temperature,salinity,and density parameters of sea ice to obtain essential insights.During the 11th Chinese National Arctic Research Expedition,acoustic velocity was measured on an ice core at a short-term ice station,however,temperature,salinity,and density were not measured.In the present work,we utilized a genetic algorithm to invert these obtained acoustic velocity data to sea ice temperature,salinity,and density parameters on the basis of the relationship between acoustic velocity and the physical properties of Arctic summer sea ice.We validated the effectiveness of this inversion procedure by comparing its findings with those of other researchers.The results indicate that within the normalized depth range of 0.43-0.94,the ranges for temperature,salinity,and density are -0.48--0.29℃,1.63-3.35,and 793.1-904.1 kg m^(-3),respectively.

Shear wave velocity model using HVSR inversion beneath Bandar Lampung City

The horizontal-to-vertical spectral ratio(HVSR)method has been used to characterize site-effect parameters that are indispensable in seismic hazard and risk-reduction studies in urban areas and rapid land-use planning.This method is widely used because it is the cheapest and simplest geophysical method for the acquisition and processing stages.In subsequent developments,the HVSR method has been widely used to determine elastic rock parameters,particularly shear wave velocity(v_(S)),through the HVSR curve inversion process.Furthermore,the v_(S)structural model can be used to delineate the presence of complex geological structures,particularly faults and sedimentary basins.Bandar Lampung is a city in Lampung Province with many fault structures and groundwater basins to the south.There are 83 HVSR measurement points around Bandar Lampung for delineating the presence of fault structures and groundwater basins.We produced the HVSR curve from the measurement results and then performed an inversion process using the particle swarm optimization algorithm to obtain v_(S)for the depth profile.Subsequently,from this profile,we produced a two-dimensional(2D)lateral and vertical model.The mean v_(S)value was calculated from all the measurement points,and we found stiff soil layers reaching depths of approximately 5 m,with a value of v_(S)<330 m/s.A bedrock layer with a velocity exceeding 1250 m/s was visible at a depth of 100 m.Based on the 2D model,the v_(S)structure shows that the city of Bandar Lampung is divided into two zones,with a NW-SE boundary.The north-middle-eastern part of the city consists of harder rocks.This harder rock is characterized by extremely high v_(S)values,starting from a depth of 50 m.In contrast,the south-middle-west exhibits a low-moderate v_(S)anomaly associated with groundwater basins SW of the city.From the 2D v_(S)structural model,fault structures can be found along the city,characterized by a contrast of v_(S)values from low to medium and from medium to high.

Inversion of walkaway VSP data in the presence of lateral velocity heterogeneity

Multi-azimuth walkaway vertical seismic profiling is an established technique for the estimation of in situ slowness surfaces and inferring anisotropy parameters.Normally,this technique requires the assumption of lateral homogeneity,which makes the horizontal slowness components at depths of downhole receivers equal to those measured at the surface.Any violations of this assumption,such as lateral heterogeneity or nonzero dip of intermediate interfaces,lead to distortions in reconstructed slowness surfaces and,consequently,to errors in estimated anisotropic parameters.In this work,we relax the assumption of lateral homogeneity and discuss how to correct vertical seismic profile data for weak lateral heterogeneity.We describe a procedure of downward continuation of recorded traveltimes that accounts for the presence of both vertical inhomogeneity and weak lateral heterogeneity,which produces correct slowness surfaces at depths of downhole receivers,noticing that sufficiently dense receiver coverage along a borehole is required to separate influences of vertical and lateral heterogeneity on measured traveltimes and obtain accurate estimates of the slowness surfaces.Once the slowness surfaces are found and a desired type of anisotropic model to be inverted is selected,the corresponding anisotropic parameters,providing the best fit to the estimated slownesses,can be obtained.We invert the slowness surfaces of P-waves for parameters of the simplest anisotropic model describing dipping fractures(transversely isotropic medium with a tilted symmetry axis).Five parameters of this model,namely,the P-wave velocity V0 in the direction of the symmetry axis,Thomsen's anisotropic coefficients e and d,the tilt n,and the azimuth b of the symmetry axis,can be estimated in a stable manner when maximum source offset is greater than half of receiver depth.

Crustal and uppermost mantle structure of the northeastern Qinghai-Xizang Plateau from joint inversion of surface wave dispersions and receiver functions with P velocity constraints

Lithospheric structure beneath the northeastern Qinghai-Xizang Plateau is of vital significance for studying the geodynamic processes of crustal thickening and expansion of the Qinghai-Xizang Plateau. We conducted a joint inversion of receiver functions and surface wave dispersions with P-wave velocity constraints using data from the Chin Array Ⅱ temporary stations deployed across the Qinghai-Xizang Plateau. Prior to joint inversion, we applied the H-κ-c method(Li JT et al., 2019) to the receiver function data in order to correct for the back-azimuthal variations in the arrival times of Ps phases and crustal multiples caused by crustal anisotropy and dipping interfaces. High-resolution images of vS, crustal thickness, and vP/vSstructures in the Qinghai-Xizang Plateau were simultaneously derived from the joint inversion. The seismic images reveal that crustal thickness decreases outward from the Qinghai-Xizang Plateau. The stable interiors of the Ordos and Alxa blocks exhibited higher velocities and lower crustal vP/vSratios. While, lower velocities and higher vP/vSratios were observed beneath the Qilian Orogen and Songpan-Ganzi terrane(SPGZ), which are geologically active and mechanically weak, especially in the mid-lower crust.Delamination or thermal erosion of the lithosphere triggered by hot asthenospheric flow contributes to the observed uppermost mantle low-velocity zones(LVZs) in the SPGZ. The crustal thickness, vS, and vP/vSratios suggest that whole lithospheric shortening is a plausible mechanism for crustal thickening in the Qinghai-Xizang Plateau, supporting the idea of coupled lithospheric-scale deformation in this region.

Converted wave AVO inversion for average velocity ratio and shear wave reflection coefficient

Based on the empirical Gardner equation describing the relationship between density and compressional wave velocity, the converted wave reflection coefficient extrema attributes for AVO analysis are proposed and the relations between the extrema position and amplitude, average velocity ratio across the interface, and shear wave reflection coefficient are derived. The extrema position is a monotonically decreasing function of average velocity ratio, and the extrema amplitude is a function of average velocity ratio and shear wave reflection coefficient. For theoretical models, the average velocity ratio and shear wave reflection coefficient are inverted from the extrema position and amplitude obtained from fitting a power function to converted wave AVO curves. Shear wave reflection coefficient sections have clearer physical meaning than conventional converted wave stacked sections and establish the theoretical foundation for geological structural interpretation and event correlation. ＂The method of inverting average velocity ratio and shear wave reflection coefficient from the extrema position and amplitude obtained from fitting a power function is applied to real CCP gathers. The inverted average velocity ratios are consistent with those computed from compressional and shear wave well logs.

Tomographic velocity inversion for ADCIGs in areas with a rugged surface

Pre-stack depth migration velocity analysis is one of the key techniques influencing image quality. As for areas with a rugged surface and complex subsurface, conventional prestack depth migration velocity analysis corrects the rugged surface to a known datum or designed surface velocity model on which to perform migration and update the velocity. We propose a rugged surface tomographic velocity inversion method based on angle-domain common image gathers by which the velocity field can be updated directly from the rugged surface without static correction for pre-stack data and improve inversion precision and efficiency. First, we introduce a method to acquire angle-domain common image gathers （ADCIGs） in rugged surface areas and then perform rugged surface tornographic velocity inversion. Tests with model and field data prove the method to be correct and effective.

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.

Wave equation tomographic velocity inversion method based on the Born/Rytov approximation

This paper discusses Born/Rytov approximation tomographic velocity inversion methods constrained by the Fresnel zone. Calculations of the sensitivity kernel function and traveltime residuals are critical in tomographic velocity inversion. Based on the Bom/Rytov approximation of the frequency-domain wave equation, we derive the traveltime sensitivity kemels of the wave equation on the band-limited wave field and simultaneously obtain the traveltime residuals based on the Rytov approximation. In contrast to single-ray tomography, the modified velocity inversion method improves the inversion stability. Tests of the near- surface velocity model and field data prove that the proposed method has higher accuracy and Computational efficiency than ray theory tomography and full waveform inversion methods.

Rock critical porosity inversion and S-wave velocity prediction

A critical porosity model is often used to calculate the dry frame elastic modulus by the rock critical porosity value which is affected by many factors. In practice it is hard for us to obtain an accurate critical porosity value and we can generally take only an empirical critical porosity value which often causes errors. In this paper, we propose a method to obtain the rock critical porosity value by inverting P-wave velocity and applying it to predict S-wave velocity. The applications of experiment and log data both show that the critical porosity inversion method can reduce the uncertainty resulting from using an empirical value in the past and provide the accurate critical porosity value for predicting S-wave velocity which significantly improves the prediction accuracy.

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.

Seismic velocity inversion based on CNN-LSTM fusion deep neural network

Based on the CNN-LSTM fusion deep neural network,this paper proposes a seismic velocity model building method that can simultaneously estimate the root mean square(RMS)velocity and interval velocity from the common-midpoint(CMP)gather.In the proposed method,a convolutional neural network(CNN)Encoder and two long short-term memory networks(LSTMs)are used to extract spatial and temporal features from seismic signals,respectively,and a CNN Decoder is used to recover RMS velocity and interval velocity of underground media from various feature vectors.To address the problems of unstable gradients and easily fall into a local minimum in the deep neural network training process,we propose to use Kaiming normal initialization with zero negative slopes of rectifi ed units and to adjust the network learning process by optimizing the mean square error(MSE)loss function with the introduction of a freezing factor.The experiments on testing dataset show that CNN-LSTM fusion deep neural network can predict RMS velocity as well as interval velocity more accurately,and its inversion accuracy is superior to that of single neural network models.The predictions on the complex structures and Marmousi model are consistent with the true velocity variation trends,and the predictions on fi eld data can eff ectively correct the phase axis,improve the lateral continuity of phase axis and quality of stack section,indicating the eff ectiveness and decent generalization capability of the proposed method.

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.

Joint Inversion of the 3D P Wave Velocity Structure of the Crust and Upper Mantle under the Southeastern Margin of the Tibetan Plateau Using Regional Earthquake and Teleseismic Data

The special seismic tectonic environment and frequent seismicity in the southeastern margin of the Qinghai-Tibet Plateau show that this area is an ideal location to study the present tectonic movement and background of strong earthquakes in China's Mainland and to predict future strong earthquake risk zones. Studies of the structural environment and physical characteristics of the deep structure in this area are helpful to explore deep dynamic effects and deformation field characteristics, to strengthen our understanding of the roles of anisotropy and tectonic deformation and to study the deep tectonic background of the seismic origin of the block＇s interior. In this paper, the three-dimensional （3D） P-wave velocity structure of the crust and upper mantle under the southeastern margin of the Qinghai-Tibet Plateau is obtained via observational data from 224 permanent seismic stations in the regional digital seismic network of Yunnan and Sichuan Provinces and from 356 mobile China seismic arrays in the southern section of the north-south seismic belt using a joint inversion method of the regional earthquake and teleseismic data. The results indicate that the spatial distribution of the P-wave velocity anomalies in the shallow upper crust is closely related to the surface geological structure, terrain and lithology. Baoxing and Kangding, with their basic volcanic rocks and volcanic clastic rocks, present obvious high-velocity anomalies. The Chengdu Basin shows low-velocity anomalies associated with the Quaternary sediments. The Xichang Mesozoic Basin and the Butuo Basin are characterised by low- velocity anomalies related to very thick sedimentary layers. The upper and middle crust beneath the Chuan-Dian and Songpan-Ganzi Blocks has apparent lateral heterogeneities, including low-velocity zones of different sizes. There is a large range of low-velocity layers in the Songpan-Ganzi Block and the sub-block northwest of Sichuan Province, showing that the middle and lower crust is relatively weak. The Sichuan Basin, which is located in the western margin of the Yangtze platform, shows high-velocity characteristics. The results also reveal that there are continuous low-velocity layer distributions in the middle and lower crust of the Daliangshan Block and that the distribution direction of the low-velocity anomaly is nearly SN, which is consistent with the trend of the Daliangshan fault. The existence of the low-velocity layer in the crust also provides a deep source for the deep dynamic deformation and seismic activity of the Daliangshan Block and its boundary faults. The results of the 3D P-wave velocity structure show that an anomalous distribution of high-density, strong-magnetic and high-wave velocity exists inside the crust in the Panxi region. This is likely related to late Paleozoic mantle plume activity that led to a large number of mafic and ultra-mafic intrusions into the crust. In the crustal doming process, the massive intrusion of mantle-derived material enhanced the mechanical strength of the crustal medium. The P-wave velocity structure also revealed that the upper mantle contains a low-velocity layer at a depth of 80-120 km in the Panxi region. The existence of deep faults in the Panxi region, which provide conditions for transporting mantle thermal material into the crust, is the deep tectonic background for the area＇s strong earthquake activity.

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.

Seismic traveltime inversion of 3D velocity model with triangulated interfaces

Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks,which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block,and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization,which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models.The damped least squares method is employed in seismic traveltime inversion,which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model.

S-wave velocity structure inferred from receiver function inversion in Tengchong volcanic area

Tengchong volcanic area is located near the impinging and underthrust margin of India and Eurasia plates. The volcanic activity is closely related to the tectonic environment. The deep structure characteristics are inferred from the receiver function inversion with the teleseismic records in the paper. The results show that the low velocity zone is influenced by the NE-trending Dayingjiang fault. The S-wave low velocity structure occurs obviously in the southern part of the fault, but unobviously in its northern part. There are low velocity zones in the shallow po-sition, which coincides with the seismicity. It also demonstrates that the low velocity zone is directly related to the thermal activity in the volcanic area. Therefore, we consider that the volcano may be alive again.

Full-waveform Velocity Inversion Based on the Acoustic Wave Equation

Full-waveform velocity inversion based on the acoustic wave equation in the time domain is investigated in this paper. The inversion is the iterative minimization of the misfit between observed data and synthetic data obtained by a numerical solution of the wave equation. Two inversion algorithms in combination with the CG method and the BFGS method are described respectively. Numerical computations for two models including the benchmark Marmousi model with complex structure are implemented. The inversion results show that the BFGS-based algorithm behaves better in inversion than the CG-based algorithm does. Moreover, the good inversion result for Marmousi model with the BFGS-based algorithm suggests the quasi-Newton methods can provide an important tool for large-scale velocity inversion. More computations demonstrate the correctness and effectives of our inversion algorithms and code.

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

Using the data of P-wave network and Zhejiang and travel time recorded at the Shanxi-reservoir seismological 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 lowvelocity anomaly zone. This may be related to water penetration.

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.