Multi-scale data joint inversion of minerals and porosity in altered igneous reservoirs—A case study in the South China Sea

There are abundant igneous gas reservoirs in the South China Sea with significant value of research,and lithology classification,mineral analysis and porosity inversion are important links in reservoir evaluation.However,affected by the diverse lithology,complicated mineral and widespread alteration,conventional logging lithology classification and mineral inversion become considerably difficult.At the same time,owing to the limitation of the wireline log response equation,the quantity and accuracy of minerals can hardly meet the exploration requirements of igneous formations.To overcome those issues,this study takes the South China Sea as an example,and combines multi-scale data such as micro rock slices,petrophysical experiments,wireline log and element cutting log to establish a set of joint inversion methods for minerals and porosity of altered igneous rocks.Specifically,we define the lithology and mineral characteristics through core slices and mineral data,and establish an igneous multi-mineral volumetric model.Then we determine element cutting log correction method based on core element data,and combine wireline log and corrected element cutting log to perform the lithology classification and joint inversion of minerals and porosity.However,it is always difficult to determine the elemental eigenvalues of different minerals in inversion.This paper uses multiple linear regression methods to solve this problem.Finally,an integrated inversion technique for altered igneous formations was developed.The results show that the corrected element cutting log are in good agreement with the core element data,and the mineral and porosity results obtained from the joint inversion based on the wireline log and corrected element cutting log are also in good agreement with the core data from X-ray diffraction.The results demonstrate that the inversion technique is applicable and this study provides a new direction for the mineral inversion research of altered igneous formations.

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.

An empirical method for joint inversion of wave and wind parameters based on SAR and wave spectrometer data

Synthetic aperture radar(SAR)and wave spectrometers,crucial in microwave remote sensing,play an essential role in monitoring sea surface wind and wave conditions.However,they face inherent limitations in observing sea surface phenomena.SAR systems,for instance,are hindered by an azimuth cut-off phenomenon in sea surface wind field observation.Wave spectrometers,while unaffected by the azimuth cutoff phenomenon,struggle with low azimuth resolution,impacting the capture of detailed wave and wind field data.This study utilizes SAR and surface wave investigation and monitoring(SWIM)data to initially extract key feature parameters,which are then prioritized using the extreme gradient boosting(XGBoost)algorithm.The research further addresses feature collinearity through a combined analysis of feature importance and correlation,leading to the development of an inversion model for wave and wind parameters based on XGBoost.A comparative analysis of this model with ERA5 reanalysis and buoy data for of significant wave height,mean wave period,wind direction,and wind speed reveals root mean square errors of 0.212 m,0.525 s,27.446°,and 1.092 m/s,compared to 0.314 m,0.888 s,27.698°,and 1.315 m/s from buoy data,respectively.These results demonstrate the model’s effective retrieval of wave and wind parameters.Finally,the model,incorporating altimeter and scatterometer data,is evaluated against SAR/SWIM single and dual payload inversion methods across different wind speeds.This comparison highlights the model’s superior inversion accuracy over other methods.

Fast 3D joint inversion of gravity and magnetic data based on cross gradient constraint

The gravity and magnetic data can be adopted to interpret the internal structure of the Earth.To improve the calculation efficiency during the inversion process and the accuracy and reliability of the reconstructed physical property models,the triple strategy is adopted in this paper to develop a fast cross-gradient joint inversion for gravity and magnetic data.The cross-gradient constraint contains solving the gradients of the physical property models and performing the cross-product calculation of their gradients.The sparse matrices are first obtained by calculating the gradients of the physical property models derived from the first-order finite difference.Then,the triple method is applied to optimize the storages and the calculations related to the gradients of the physical property models.Therefore,the storage compression amount of the calculations related to the gradients of the physical property models and the cross-gradient constraint are reduced to one-fold of the number of grid cells at least,and the compression ratio increases with the increase of the number of grid cells.The test results from the synthetic data and field data prove that the structural coupling is achieved by using the fast cross-gradient joint inversion method to effectively reduce the multiplicity of solutions and improve the computing efficiency.

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.

Joint inversion of gravity and vertical gradient data based on modified structural similarity index for the structural and petrophysical consistency constraint

Joint inversion is one of the most effective methods for reducing non-uniqueness for geophysical inversion.The current joint inversion methods can be divided into the structural consistency constraint and petrophysical consistency constraint methods,which are mutually independent.Currently,there is a need for joint inversion methods that can comprehensively consider the structural consistency constraints and petrophysical consistency constraints.This paper develops the structural similarity index(SSIM)as a new structural and petrophysical consistency constraint for the joint inversion of gravity and vertical gradient data.The SSIM constraint is in the form of a fraction,which may have analytical singularities.Therefore,converting the fractional form to the subtractive form can solve the problem of analytic singularity and finally form a modified structural consistency index of the joint inversion,which enhances the stability of the SSIM constraint applied to the joint inversion.Compared to the reconstructed results from the cross-gradient inversion,the proposed method presents good performance and stability.The SSIM algorithm is a new joint inversion method for petrophysical and structural constraints.It can promote the consistency of the recovered models from the distribution and the structure of the physical property values.Then,applications to synthetic data illustrate that the algorithm proposed in this paper can well process the synthetic data and acquire good reconstructed results.

Curupira V1.0: Joint Inversion of VES and TEM for Environmental and Mass Movements Studies

An innovative inversion code, named “Curupira v1.0”, has been developed using Matlab to determine the vertical distribution of resistivity beneath the subsoil. The program integrates Vertical Electrical Sounding (VES), successful in shallow subsurface exploration and Time Domain Electromagnetic (TEM) techniques, better suited for deeper exploration, both of which are widely employed in geophysical exploration. These methodologies involve calculating subsurface resistivity through appropriate inversion processes. To address the ill-posed nature of inverse problems in geophysics, a joint inversion scheme combining VES and TEM data has been incorporated into Curupira v1.0. The software has been tested on both synthetic and real-world data, the latter of which was acquired from the Parana sedimentary basin which we summarise here. The results indicate that the joint inversion of VES and TEM techniques offers improved recovery of simulated models and demonstrates significant potential for hydrogeological studies.

Joint PP and PS AVO inversion based on Bayes theorem

Based on the Aki-Richards approximate equations for reflection coefficients and Bayes theorem, we developed an inversion method to estimate P- and S-wave velocity contrasts and density contrast from combined PP and PS data. This method assumes that the parameters satisfy a normal distribution and introduces the covariance matrix to describe the degree of correlation between the parameters and thus to improve the inversion stability. Then, we suppose that the parameter sequence is subject to the Cauchy distribution and employs another matrix Q to describe the parameter sequence sparseness to improve the inversion result resolution. Tests on both synthetic and real multi-component data prove that this method is valid, efficient, more stable, and more accurate compared to methods using PP data only.

Generalized linear joint PP-PS inversion based on two constraints

Conventional joint PP-PS inversion is based on approximations of the Zoeppritz equations and assumes constant VP/VS;therefore,the inversion precision and stability cannot satisfy current exploration requirements.We propose a joint PP-PS inversion method based on the exact Zoeppritz equations that combines Bayesian statistics and generalized linear inversion.A forward model based on the exact Zoeppritz equations is built to minimize the error of the approximations in the large-angle data,the prior distribution of the model parameters is added as a regularization item to decrease the ill-posed nature of the inversion,low-frequency constraints are introduced to stabilize the low-frequency data and improve robustness,and a fast algorithm is used to solve the objective function while minimizing the computational load.The proposed method has superior antinoising properties and well reproduces real data.

Joint inversion of gravity and magnetic data for a two-layer model

Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose a joint gravity and magnetic inversion methodfor two-layer models by concentrating on the relationship between the change of thicknessI and position of the middle layer and anomaly and discuss the effects of the key parameters. Model tests and application to field data show the validity of this method.

AVO approximation for PS-wave and its application in PP/PS joint inversion

Multi-component exploration has many advantages over ordinary P-wave exploration. PP/PS joint AVO analysis and inversion are useful and powerful methods to discriminate between reservoir and non-productive lithology. In this paper, we derive a new PS-wave reflection coefficient approximation equation which is more accurate at larger incidence angles. The equation is simplified for small incidence angles, which makes AVO analysis clearer and easier for angles less than 30 degrees. Based on this approximation, a PP/PS joint inversion is introduced. A real data example shows that oil sands, brine sands and shales can be differentiated based on the P- to S-wave velocity ratio from the PP/PS joint inversion. Fluid factors and Poisson＇s ratio also indicate an anomaly in the target zone at the oil well location.

Nonlinear joint PP-PS AVO inversion based on improved Bayesian inference and LSSVM

Multiwave seismic technology promotes the application of joint PP–PS amplitude versus offset (AVO) inversion;however conventional joint PP–PS AVO inversioan is linear based on approximations of the Zoeppritz equations for multiple iterations. Therefore the inversion results of P-wave, S-wave velocity and density exhibit low precision in the faroffset;thus, the joint PP–PS AVO inversion is nonlinear. Herein, we propose a nonlinear joint inversion method based on exact Zoeppritz equations that combines improved Bayesian inference and a least squares support vector machine (LSSVM) to solve the nonlinear inversion problem. The initial parameters of Bayesian inference are optimized via particle swarm optimization (PSO). In improved Bayesian inference, the optimal parameter of the LSSVM is obtained by maximizing the posterior probability of the hyperparameters, thus improving the learning and generalization abilities of LSSVM. Then, an optimal nonlinear LSSVM model that defi nes the relationship between seismic refl ection amplitude and elastic parameters is established to improve the precision of the joint PP–PS AVO inversion. Further, the nonlinear problem of joint inversion can be solved through a single training of the nonlinear inversion model. The results of the synthetic data suggest that the precision of the estimated parameters is higher than that obtained via Bayesian linear inversion with PP-wave data and via approximations of the Zoeppritz equations. In addition, results using synthetic data with added noise show that the proposed method has superior anti-noising properties. Real-world application shows the feasibility and superiority of the proposed method, as compared with Bayesian linear inversion.

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.

Rupture process of the 2011 Tohoku earthquake from the joint inversion of teleseismic and GPS data

Teleseismic and GPS data were jointly inverted for the rupture process of the 2011 Tohoku earthquake. The inversion results show that it is a bilateral rupture event with an average rupture velocity less than 2.0 km/s along the fault strike direction. The source rupture process consists of three sub-events, the first oc- curred near the hypocenter and the rest two ruptured along the up-dip direction and broke the sea bed, causing a maximum slip of about 30 m. The large-scale sea bed breakage may account for the tremendous tsunami disaster which resulted in most of the death and missing in this mega earthquake.

2D joint inversion of CSAMT and magnetic data based on cross-gradient theory

A two-dimensional forward and backward algorithm for the controlled-source audio-frequency magnetotelluric （CSAMT） method is developed to invert data in the entire region （near, transition, and far） and deal with the effects of artificial sources. First, a regularization factor is introduced in the 2D magnetic inversion, and the magnetic susceptibility is updated in logarithmic form so that the inversion magnetic susceptibility is always positive. Second, the joint inversion of the CSAMT and magnetic methods is completed with the introduction of the cross gradient. By searching for the weight of the cross-gradient term in the objective function, the mutual influence between two different physical properties at different locations are avoided. Model tests show that the joint inversion based on cross-gradient theory offers better results than the single-method inversion. The 2D forward and inverse algorithm for CSAMT with source can effectively deal with artificial sources and ensures the reliability of the final joint inversion algorithm.

Source model of the 11th July 2004 Zhongba earthquake revealed from the joint inversion of InSAR and seismological data

We use interferometric synthetic aperture radar （InSAR） and broadband seismic waveform data to estimate a source model of the 11th July, 2004 M W 6.2 Zhongba earthquake, Tibet of China. This event occurred within the seismically active zone of southwestern Tibetan Plateau where the east-west extension of the upper crust is observed. Because of limitations in one pair of InSAR data available, there are trade-offs among centroid depth, rupture area and amount of slip. Available seismic data tightly constrain the focal mechanism and centroid depth of the earthquake but not the horizontal location. Together, two complementary data sets can be used to identify the actual fault plane, better constrain the slip model and event location. We first use regional seismic waveform to estimate point source mechanism, then InSAR data is used to obtain better location. Finally, a joint inversion of teleseismic P-waves and InSAR data is performed to obtain a distributed model. Our preferred point source mechanism indicates a seismic moment of ～2.2×10 18 N·m （～M W 6.2）, a fault plane solution of 171° （342 ° ）/42 ° （48 ° ）/-83 ° （-97 ° ）, corresponding to strike/dip/rake, and a depth of 11 km. The fault plane with strike of 171 ？ and dip of 42° is identified as the ruptured fault with the aid of InSAR data. The preferred source model features compact area of slips between depth of 5–11 km and 10 km along strike with maximum slip amplitude of about 1.5 m.

Elastic modulus extraction based on generalized pre-stack PP–PS wave joint linear inversion

Joint PP–PS inversion offers better accuracy and resolution than conventional P-wave inversion. P-and S-wave elastic moduli determined through data inversions are key parameters for reservoir evaluation and fluid characterization. In this paper, starting with the exact Zoeppritz equation that relates P-and S-wave moduli, a coefficient that describes the reflections of P-and converted waves is established. This method effectively avoids error introduced by approximations or indirect calculations, thus improving the accuracy of the inversion results. Considering that the inversion problem is ill-posed and that the forward operator is nonlinear, prior constraints on the model parameters and modified low-frequency constraints are also introduced to the objective function to make the problem more tractable. This modified objective function is solved over many iterations to continuously optimize the background values of the velocity ratio, which increases the stability of the inversion process. Tests of various models show that the method effectively improves the accuracy and stability of extracting P and S-wave moduli from underdetermined data. This method can be applied to provide inferences for reservoir exploration and fluid extraction.

The distribution of deep source rocks in the GS Sag:joint MT-gravity modeling and constrained inversion

The coal-bearing strata of the deep Upper Paleozoic in the GS Sag have high hydrocarbon potential. Because of the absence of seismic data, we use electromagnetic （MT） and gravity data jointly to delineate the distribution of deep targets based on well logging and geological data. First, a preliminary geological model is established by using three-dimensional （3D） MT inversion results. Second, using the formation density and gravity anomalies, the preliminary geological model is modified by interactive inversion of the gravity data. Then, we conduct MT-constrained inversion based on the modified model to obtain an optimal geological model until the deviations at all stations are minimized. Finally, the geological model and a seismic profile in the middle of the sag is analysed. We determine that the deep reflections of the seismic profile correspond to the Upper Paleozoic that reaches thickness up to 800 m. The processing of field data suggests that the joint MT-gravity modeling and constrained inversion can reduce the multiple solutions for single geophysical data and thus improve the recognition of deep formations. The MT-constrained inversion is consistent with the geological features in the seismic section. This suggests that the joint MT and gravity modeling and constrained inversion can be used to delineate deep targets in similar basins.

Joint inversion of gravity and seismic data along a profile across the seismogenic fault of 2010 Yushu Ms7. 1 earthquake

Yushu Ms7.1 earthquake occurred on the Ganzi-Yushu fault zone, across which we carried out a joint relative-gravity and seismic-reflection survey, and then performed a gravity inversion constrained by the seismic-reflection result. Based on the data of complete Bouguer gravity anomaly and seismic reflection, we obtained a layered interface structure in deep crust down to Moho. Our study showed that the inversion could reveal the interfaces of strata along the survey profile and the directions of regional faults in two-dimension. From the characteristics of the observed topography of the Moho basement, we tentatively confirmed that the uplift of eastern edge of Qinghai-Tibet plateau was caused by the subduetion of the Indian plate.

3D inversion modeling of joint gravity and magnetic data based on a sinusoidal correlation constraint

Joint inversion based on a correlation constraint utilizes a linear correlation function as a structural constraint.The linear correlation function contains a denominator,which may result in a singularity as the objective function is optimized,leading to an unstable inversion calculation.To improve the robustness of this calculation,this paper proposes a new method in which a sinusoidal correlation function is employed as the structural constraint for joint inversion instead of the conventional linear correlation function.This structural constraint does not contain a denominator,thereby preventing a singularity.Compared with the joint inversion method based on a cross-gradient constraint,the joint inversion method based on a sinusoidal correlation constraint exhibits good performance.An application to actual data demonstrates that this method can process real data.