The major storage space types in the carbonate reservoir in the Ordovician in the TZ45 area are secondary dissolution caves.For the prediction of caved carbonate reservoir,post-stack methods are commonly used in the o...The major storage space types in the carbonate reservoir in the Ordovician in the TZ45 area are secondary dissolution caves.For the prediction of caved carbonate reservoir,post-stack methods are commonly used in the oilfield at present since pre-stack inversion is always limited by poor seismic data quality and insufficient logging data.In this paper,based on amplitude preserved seismic data processing and rock-physics analysis,pre-stack inversion is employed to predict the caved carbonate reservoir in TZ45 area by seriously controlling the quality of inversion procedures.These procedures mainly include angle-gather conversion,partial stack,wavelet estimation,low-frequency model building and inversion residual analysis.The amplitude-preserved data processing method can achieve high quality data based on the principle that they are very consistent with the synthetics.Besides,the foundation of pre-stack inversion and reservoir prediction criterion can be established by the connection between reservoir property and seismic reflection through rock-physics analysis.Finally,the inversion result is consistent with drilling wells in most cases.It is concluded that integrated with amplitude-preserved processing and rock-physics,pre-stack inversion can be effectively applied in the caved carbonate reservoir prediction.展开更多
A three-dimensional(3D)electromagnetic(EM)inversion algorithm based on the nonlinear conjugate gradient(NLCG)method and a two-color plane Gauss-Seidel(GS)multigrid(MG)forward solver is developed to improve inversion e...A three-dimensional(3D)electromagnetic(EM)inversion algorithm based on the nonlinear conjugate gradient(NLCG)method and a two-color plane Gauss-Seidel(GS)multigrid(MG)forward solver is developed to improve inversion efficiency.The results indicate that the computational efficiency of each inversion can be improved by approximately a factor of three by using the proposed MG solver.First,the accuracy of the MG solver is validated through a test on a synthetic model.Next,the numerical performance of the inversion algorithm is evaluated using this model.Finally,the inversion algorithm is applied to a field EM data collected at the Beiya gold polymetallic ore district.A 3D resistivity model is obtained,and the formation process of the metal ore is analyzed.展开更多
Co-assembling chiral molecules with achiral compounds via non-covalent interactions like areneperfluoroarene(AP) interactions offers an effective approach for fabricating chiral functional materials.Herein,chiral mole...Co-assembling chiral molecules with achiral compounds via non-covalent interactions like areneperfluoroarene(AP) interactions offers an effective approach for fabricating chiral functional materials.Herein,chiral molecules L/D-PF1 and L/D-PF2 with pyrene groups were synthesized and its chiroptical properties upon co-assembly with achiral compound octafluoronaphthalene(OFN) through AP interaction were systemically studied.The co-assembly of L/D-PF1/OFN and L/D-PF2/OFN exhibited distinct chiroptical properties such as circular dichroism(CD) and circularly polarized luminescence(CPL) signals.Chirality transfer from the chirality center of L/D-PF1 and L/D-PF2 to the achiral OFN and chiral amplification were successfully achieved.Besides,no significant CPL signal was observed in the self-assembly of L/DPF1 or L/D-PF2 while co-assembly with OFN exhibited obvious CPL amplification induced by AP interaction.Notably,a reversal CD signal and CPL signal could be observed in L/D-PF2/OFN when the molar ratio changed from 1:1 to 1:2 while not found in L/D-PF1/OFN,indicating that that minor structural changes of molecules could cause large changes in assembly.In addition,a series of computational calculations were conducted to verify the AP interaction between L-PF1/L-PF2 and OFN.This work demonstrated that arene-perfluoroarene interaction could drive chiral transfer,chiral amplification and chiral inversion and provided a new method for the preparation of chiroptical materials.展开更多
To address the challenges of complex fluvial sandbody distribution and difficult remaining oil recovery in mature continental oilfields,this study focuses on key issues in reservoir identification such as ambiguous na...To address the challenges of complex fluvial sandbody distribution and difficult remaining oil recovery in mature continental oilfields,this study focuses on key issues in reservoir identification such as ambiguous narrow-channel boundaries and subdivision of multi-stage superimposed sandbodies.Taking the Upper Cretaceous continental sandstone in the Sazhong Oilfield of the Daqing Placanticline as an example,a technical system integrating OVT high-resolution processing,multi-attribute fusion,and varible-scale inversion was developed to establish a complete workflow from seismic processing to reservoir prediction and remaining oil recovery.The following results are obtained.First,the Offset Vector Tile(OVT)seismic processing technology is extended,for the first time,from fracture imaging to sandbody prediction,in order to address the weak seismic responses from boundaries of narrow and thin sandbodies.A geology-oriented OVT partitioning method is developed to significantly improve the imaging accuracy,enabling identification of channel sandbodies as narrow as 50 m.Second,an amplitude-coherence dual-attribute fusion method is proposed for predicting narrow channel boundaries between wells.Constrained by a sedimentary unit-level sequence chronostratigraphic framework,this method accurately delineates 800-2000 m long subaqueous distributary channels with bifurcation-convergence features.Third,considering the superimposition of multi-stage channels,a three-level variable-scale stratigraphic model(sandstone groups,sublayers,sedimentary units)is constructed to overcome single-scale modeling limitations,successfully characterizing key sedimentary features like meandering river“cut-offs”through 3D seismic inversion.Based on these advances,a direct link between seismic prediction and remaining oil recovery is established.The horizontal wells deployed using narrow-channel predictions encountered oil-bearing sandstones in the horizontal section by 97%,and achieved initial daily production of 12.5 t per well.Precise identification of individual channel boundaries within 17 composite sandbodies guided recovery processes in 135 wells,yielding an average daily increase of 2.8 t per well and a cumulative increase of 13.6×10^(4)t.展开更多
The forward model of optical fiber strain induced by fractures,together with the associated model resolution matrix,is used to demonstrate the interpretability of fracture parameters once the fracture intersects the f...The forward model of optical fiber strain induced by fractures,together with the associated model resolution matrix,is used to demonstrate the interpretability of fracture parameters once the fracture intersects the fiber.A regularized inversion framework for fracture parameters is established to evaluate the influence of measured data quality on the accuracy of iterative regularized inversion.An interpretation approach for both fracture width and height is proposed,and the synthetic forward data with measurement error and field examples are employed to validate the accuracy of the simultaneous inversion of fracture width and height.The results indicate that,after the fracture contacts the fiber,the strain response is strongly sensitive only to the fracture parameters at the intersection location,whereas the interpretability of parameters at other locations remains limited.The iterative regularized inversion method effectively suppresses the impact of measurement error and exhibits high computational efficiency,showing clear advantages for inversion applications.When incorporating the first-order regularization with a Neumann boundary constraint on the tip width,the inverted fracture-width distribution becomes highly sensitive to fracture height;thus,combined with a bisection strategy,simultaneous inversion of fracture width and height can be achieved.Examination using the model resolution matrix,noisy synthetic data,and field data confirms that the iterative regularized inversion model for fracture width and height provides high interpretive accuracy and can be applied to the calculation and analysis of fracture width,fracture height,net pressure and other parameters.展开更多
Young's modulus and Poisson's ratio are crucial parameters for reservoir characterization and rock brittleness evaluation.Conventional methods often rely on indirect computation or approximations of the Zoeppr...Young's modulus and Poisson's ratio are crucial parameters for reservoir characterization and rock brittleness evaluation.Conventional methods often rely on indirect computation or approximations of the Zoeppritz equations to estimate Young's modulus,which can introduce cumulative errors and reduce the accuracy of inversion results.To address these issues,this paper introduces the analytical solution of the Zoeppritz equation into the inversion process.The equation is re-derived and expressed in terms of Young's modulus,Poisson's ratio,and density.Within the Bayesian framework,we construct an objective function for the joint inversion of PP and PS waves.Traditional gradient-based algorithms often suffer from low precision and the computational complexity.In this study,we address limitations of conventional approaches related to low precision and complicated code by using Circle chaotic mapping,Levy flights,and Gaussian mutation to optimize the quantum particle swarm optimization(QPSO),named improved quantum particle swarm optimization(IQPSO).The IQPSO demonstrates superior global optimization capabilities.We test the proposed inversion method with both synthetic and field data.The test results demonstrate the proposed method's feasibility and effectiveness,indicating an improvement in inversion accuracy over traditional methods.展开更多
With the development of exploration of oil and gas resources,the requirements for seismic inversion results are getting more accurate.In particular,it is hoped that the distribution patterns of oil and gas reservoirs ...With the development of exploration of oil and gas resources,the requirements for seismic inversion results are getting more accurate.In particular,it is hoped that the distribution patterns of oil and gas reservoirs can be finely characterized,and the seismic inversion results can clearly characterize the location of stratigraphic boundaries and meet the needs of accurate geological description.Specifically,for pre-stack AVO inversion,it is required to be able to distinguish smaller geological targets in the depth or time domain,and clearly depict the vertical boundaries of the geological objects.In response to the above requirements,we introduce the preprocessing regularization of the adaptive edge-preserving smooth filter into the pre-stack AVO elastic parameter inversion to clearly invert the position of layer boundary and improve the accuracy of the inversion results.展开更多
The extensive application of pre-stack depth migration has produced huge volumes of seismic data,which allows for the possibility of developing seismic inversions of reservoir properties from seismic data in the depth...The extensive application of pre-stack depth migration has produced huge volumes of seismic data,which allows for the possibility of developing seismic inversions of reservoir properties from seismic data in the depth domain.It is difficult to estimate seismic wavelets directly from seismic data due to the nonstationarity of the data in the depth domain.We conduct a velocity transformation of seismic data to make the seismic data stationary and then apply the ridge regression method to estimate a constant seismic wavelet.The estimated constant seismic wavelet is constructed as a set of space-variant seismic wavelets dominated by velocities at different spatial locations.Incorporating the weighted superposition principle,a synthetic seismogram is generated by directly employing the space-variant seismic wavelets in the depth domain.An inversion workflow based on the model-driven method is developed in the depth domain by incorporating the nonlinear conjugate gradient algorithm,which avoids additional data conversions between the time and depth domains.The impedance inversions of the synthetic and field seismic data in the depth domain show good results,which demonstrates that seismic inversion in the depth domain is feasible.The approach provides an alternative for forward numerical analyses and elastic property inversions of depth-domain seismic data.It is advantageous for further studies concerning the stability,accuracy,and efficiency of seismic inversions in the depth domain.展开更多
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 cha...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.展开更多
For a typical marine shale reservoir in the Jiaoshiba area, Sichuan Basin of China, P-impedance is sensitive for identifying lithology but not suitable for indicating good shale reservoirs. In comparison, density is a...For a typical marine shale reservoir in the Jiaoshiba area, Sichuan Basin of China, P-impedance is sensitive for identifying lithology but not suitable for indicating good shale reservoirs. In comparison, density is an important quantity, which is sensitive for identifying the organic-rich mud shale from non-organic-rich mud shale. Due to the poor data quality and incidence angle range, density cannot be easily inverted by directly solving the ill-posed pre-stack seismic inversion in this area. Meanwhile, the traditional density regularizations implemented by directly using the more robust P-impedance inversion tend to be inaccurate for recovering density for this shale reservoir. In this paper, we combine the P-impedance and the minus uranium to construct the pseudo-P-impedance（PIp） at well locations. The PIp is observed to be sensitive for identifying organic-rich mud shale and has a good correlation with density in this area. We employ the PIp–density relation into the pre-stack inversion framework to estimate density. Three types of regularization are tested on both numerical and field data： These are no regularization, traditional regularization and the proposed approach. It is observed that the proposed method is better for recovering the density of organic-rich mud shale in the Jiaoshiba area.展开更多
Incremental Nonlinear Dynamic Inversion(INDI)is a control approach that has gained popularity in flight control over the past decade.Besides the INDI law,several common additional components complement an INDI-based c...Incremental Nonlinear Dynamic Inversion(INDI)is a control approach that has gained popularity in flight control over the past decade.Besides the INDI law,several common additional components complement an INDI-based controller.This paper,the second part of a two-part series of surveys on INDI,aims to summarize the modern trends in INDI and its related components.Besides a comprehensive components specification,it addresses their most common challenges,compares different variants,and discusses proposed advances.Further important aspects of INDI are gain design,stability,and robustness.This paper also provides an overview of research conducted concerning these aspects.This paper is written in a tutorial style to familiarize researchers with the essential specifics and pitfalls of INDI and its components.At the same time,it can also serve as a reference for readers already familiar with INDI.展开更多
An innovative gradient inversion approach employing the natural element method within the framework of least square regularization was proposed to enhance the quantitative interpretation of self-potential(SP)data orig...An innovative gradient inversion approach employing the natural element method within the framework of least square regularization was proposed to enhance the quantitative interpretation of self-potential(SP)data originating from mineral polarization.The results indicated that the natural element method effectively addressed the challenge of subdividing complex resistivity models and aided in the accurate forward calculation of SP.By applying this approach to synthetic SP data and lab-measured SP data associated with redox electrochemical half-cell reactions of iron−copper metal blocks within the geobattery model,the 3D fine structure of buried orebody models was successfully reconstructed and the spatial distribution of SP current sources was mapped.This study significantly contributes to understanding the quantitative relationship between the polarization process of metal deposits and their corresponding SP responses and provides a valuable reference for delineating metal deposits in both terrestrial and marine environments through SP surveys.展开更多
The brittleness index(BI)is crucial for predicting engineering sweet spots and designing fracturing operations in shale oil reservoir exploration and development.Seismic amplitude variation with offset(AVO)inversion i...The brittleness index(BI)is crucial for predicting engineering sweet spots and designing fracturing operations in shale oil reservoir exploration and development.Seismic amplitude variation with offset(AVO)inversion is commonly used to obtain the BI.Traditionally,velocity,density,and other parameters are firstly inverted,and the BI is then calculated,which often leads to accumulated errors.Moreover,due to the limited of well-log data in field work areas,AVO inversion typically faces the challenge of limited information,resulting in not high accuracy of BI derived by existing AVO inversion methods.To address these issues,we first derive an AVO forward approximation equation that directly characterizes the BI in P-wave reflection coefficients.Based on this,an intelligent AVO inversion method,which combines the advantages of traditional and intelligent approaches,for directly obtaining the BI is proposed.A TransUnet model is constructed to establish the strong nonlinear mapping relationship between seismic data and the BI.By incorporating a combined objective function that is constrained by both low-frequency parameters and training samples,the challenge of limited samples is effectively addressed,and the direct inversion of the BI is stably achieved.Tests on model data and applications on field data demonstrate the feasibility,advancement,and practicality of the proposed method.展开更多
As geological exploration conditions become increasingly complex, meeting the requirements of precise geological exploration necessitates the development of a controlled-source audio magnetotelluric (CSAMT) inversion ...As geological exploration conditions become increasingly complex, meeting the requirements of precise geological exploration necessitates the development of a controlled-source audio magnetotelluric (CSAMT) inversion method that considers anisotropy to improve the effectiveness of inversion accuracy and interpretation accuracy of data. This study is based on the 3D fi nite-diff erence forward modeling of axis anisotropy using the reciprocity theorem to calculate the Jacobian matrix by applying the search method to automatically search for the Lagrange operator. The aim is to establish inversion iteration equations to achieve the axis anisotropic Occam's 3D inversion of tensor CSAMT in data space. Further, we obtain an underground axis anisotropic 3D geoelectric model by inverting the impedance data of tensor CSAMT. Two synthetic data examples show that using the isotropic tensor CSAMT algorithm to directly invert data in anisotropic media can generate false anomalies, leading to incorrect geological interpretations. Meanwhile, the proposed anisotropic inversion algorithm can eff ectively improve the accuracy of data inversion in anisotropic media. Further, the inversion examples verify the eff ectiveness and stability of the algorithm.展开更多
The Pamir Plateau,at the northwestern margin of the Tibetan Plateau,is a key region for investigating continental collision and plateau uplifting.To probe its deep structure,we collected seismic data from 263 stations...The Pamir Plateau,at the northwestern margin of the Tibetan Plateau,is a key region for investigating continental collision and plateau uplifting.To probe its deep structure,we collected seismic data from 263 stations across 11 research projects.We applied cross-correlation to noise data and extracted surface wave dispersion data from cross-correlation functions.The extracted dispersion data were subsequently inverted using a 3-D transdimensional Bayesian inversion method(rj-3 DMcMC).The inversion result reveals several crustal low-velocity zones(LVZs).Their formation is likely related to crustal thickening,the exposure of gneiss domes,and thicker sedimentary sequences compared to surrounding areas.In the lower crust and upper mantle,the LVZs in southern Pamir and southeastern Karakoram evolve into high-velocity zones,which expand northeastward with increasing depth.This suggests northward underthrusting of the Indian Plate.We also analyzed the Moho using both the standard deviation of S-wave velocity and the S-wave velocity structure.Results show that significant variations in velocity standard deviation reliably delineate the Moho interface.展开更多
Ultrasound computed tomography(USCT)is a noninvasive biomedical imaging modality that offers insights into acoustic properties such as the sound speed(SS)and acoustic attenuation(AA)of the human body,enhancing diagnos...Ultrasound computed tomography(USCT)is a noninvasive biomedical imaging modality that offers insights into acoustic properties such as the sound speed(SS)and acoustic attenuation(AA)of the human body,enhancing diagnostic accuracy and therapy planning.Full waveform inversion(FWI)is a promising USCT image reconstruction method that optimizes the parameter fields of a wave propagation model via gradient-based optimization.However,twodimensional FWI methods are limited by their inability to account for three-dimensional wave propagation in the elevation direction,resulting in image artifacts.To address this problem,we propose a three-dimensional time-domain full waveform inversion algorithm to reconstruct the SS and AA distributions on the basis of a fractional Laplacian wave equation,adjoint field formulation,and gradient descent optimization.Validated by two sets of simulations,the proposed algorithm has potential for generating high-resolution and quantitative SS and AA distributions.This approach holds promise for clinical USCT applications,assisting early disease detection,precise abnormality localization,and optimized treatment planning,thus contributing to better healthcare outcomes.展开更多
Classical multi-channel technology can significantly reduce the pre-stack seismic inversion uncertainty, especially for complex geology such as high dipping structures. However, due to the consideration of complex str...Classical multi-channel technology can significantly reduce the pre-stack seismic inversion uncertainty, especially for complex geology such as high dipping structures. However, due to the consideration of complex structure or reflection features, the existing multi-channel inversion methods have to adopt the highly time-consuming strategy of arranging seismic data trace-by-trace, limiting its wide application in pre-stack inversion. A fast pre-stack multi-channel inversion constrained by seismic reflection features has been proposed to address this issue. The key to our method is to re-characterize the reflection features to directly constrain the pre-stack inversion through a Hadamard product operator without rearranging the seismic data. The seismic reflection features can reflect the distribution of the stratum reflection interface, and we obtained them from the post-stack profile by searching the shortest local Euclidean distance between adjacent seismic traces. Instead of directly constructing a large-size reflection features constraint operator advocated by the conventional methods, through decomposing the reflection features along the vertical and horizontal direction at a particular sampling point, we have constructed a computationally well-behaved constraint operator represented by the vertical and horizontal partial derivatives. Based on the Alternating Direction Method of Multipliers (ADMM) optimization, we have derived a fast algorithm for solving the objective function, including Hadamard product operators. Compared with the conventional reflection features constrained inversion, the proposed method is more efficient and accurate, proved on the Overthrust model and a field data set.展开更多
Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this cha...Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this challenge,nonlinear stress boundaries for a numerical model are determined through regression analysis of a series of nonlinear coefficient matrices,which are derived from the bubbling method.Considering the randomness and flexibility of the bubbling method,a parametric study is conducted to determine recommended ranges for these parameters,including the standard deviation(σb)of bubble radii,the non-uniform coefficient matrix number(λ)for nonlinear stress boundaries,and the number(m)and positions of in situ stress measurement points.A model case study provides a reference for the selection of these parameters.Additionally,when the nonlinear in situ stress inversion method is employed,stress distortion inevitably occurs near model boundaries,aligning with the Saint Venant's principle.Two strategies are proposed accordingly:employing a systematic reduction of nonlinear coefficients to achieve high inversion accuracy while minimizing significant stress distortion,and excluding regions with severe stress distortion near the model edges while utilizing the central part of the model for subsequent simulations.These two strategies have been successfully implemented in the nonlinear in situ stress inversion of the Xincheng Gold Mine and have achieved higher inversion accuracy than the linear method.Specifically,the linear and nonlinear inversion methods yield root mean square errors(RMSE)of 4.15 and 3.2,and inversion relative errors(δAve)of 22.08%and 17.55%,respectively.Therefore,the nonlinear inversion method outperforms the traditional multiple linear regression method,even in the presence of a systematic reduction in the nonlinear stress boundaries.展开更多
Based on waveform fitting,full waveform inversion(FWI)is an important inversion method with the ability to reconstruct multi-parameter models in high precision.However,the strong nonlinear equation used in FWI present...Based on waveform fitting,full waveform inversion(FWI)is an important inversion method with the ability to reconstruct multi-parameter models in high precision.However,the strong nonlinear equation used in FWI presents the following challenges,such as low convergence efficiency,high dependence on the initial model,and the energy imbalance in deep region of the inverted model.To solve these inherent problems,we develop a timedomain elastic FWI method based on gradient preconditioning with the following details:(1)the limited memory Broyden Fletcher Goldfarb Shanno method with faster convergence is adopted to im-prove the inversion stability;(2)a multi-scaled inversion strategy is used to alleviate the nonlinear inversion instead of falling into the local minimum;(3)in addition,the pseudo-Hessian preconditioned illumination operator is involved for preconditioning the parameter gradients to improve the illumination equilibrium degree of deep structures.Based on the programming implementation of the new method,a deep depression model with five diffractors is used for testing.Compared with the conventional elastic FWI method,the technique proposed by this study has better effectiveness and accuracy on the inversion effect and con-vergence,respectively.展开更多
Subsurface reservoirs commonly exhibit layered structures.Conventional amplitude variation with angle(AVA)inversion,which relies on the Zoeppritz equation and its approximations,often fails to accurately estimate elas...Subsurface reservoirs commonly exhibit layered structures.Conventional amplitude variation with angle(AVA)inversion,which relies on the Zoeppritz equation and its approximations,often fails to accurately estimate elastic parameters because it assumes single-interface models and ignores multiple reflections and transmission losses.To address these limitations,this study proposes a novel prestack time-frequency domain joint inversion method that utilizes the reflection matrix method(RMM)as the forward operator.The RMM accurately simulates wave propagation in layered media,while the joint inversion framework minimizes the misfit between observed and synthetic data in both the time and frequency domains.By incorporating Bayesian theory to optimize the inversion process,the method effectively balances contributions from both time-domain waveforms and frequency-domain spectral information through a weighting factor.Tests on both synthetic data and field data demonstrate that the proposed method outperforms conventional AVA inversion and time-domain waveform inversion in accuracy and robustness.Furthermore,the method demonstrates good robustness against variations in initial models,random noise,and coherent noise interference.This study provides a practical and effective approach for high-precision reservoir characterization,with potential applications in complex layered media.展开更多
基金supported by National Basic Research Program(2006CB202304)of Chinaco-supported by the National Basic Research Program of China(Grant No.2011CB201103)the National Science and Technology Major Project of China(Grant No.2011ZX05004003)
文摘The major storage space types in the carbonate reservoir in the Ordovician in the TZ45 area are secondary dissolution caves.For the prediction of caved carbonate reservoir,post-stack methods are commonly used in the oilfield at present since pre-stack inversion is always limited by poor seismic data quality and insufficient logging data.In this paper,based on amplitude preserved seismic data processing and rock-physics analysis,pre-stack inversion is employed to predict the caved carbonate reservoir in TZ45 area by seriously controlling the quality of inversion procedures.These procedures mainly include angle-gather conversion,partial stack,wavelet estimation,low-frequency model building and inversion residual analysis.The amplitude-preserved data processing method can achieve high quality data based on the principle that they are very consistent with the synthetics.Besides,the foundation of pre-stack inversion and reservoir prediction criterion can be established by the connection between reservoir property and seismic reflection through rock-physics analysis.Finally,the inversion result is consistent with drilling wells in most cases.It is concluded that integrated with amplitude-preserved processing and rock-physics,pre-stack inversion can be effectively applied in the caved carbonate reservoir prediction.
基金financially supported by the National Science and Technology Major Project,China(No.2024ZD1002100)the National Natural Science Foundation of China(Nos.42330801,42474112,42504062)+1 种基金the China Postdoctoral Science Foundation(No.2024M761704)Shuimu Tsinghua Scholar Program of Tsinghua University,China(No.2024SM114)。
文摘A three-dimensional(3D)electromagnetic(EM)inversion algorithm based on the nonlinear conjugate gradient(NLCG)method and a two-color plane Gauss-Seidel(GS)multigrid(MG)forward solver is developed to improve inversion efficiency.The results indicate that the computational efficiency of each inversion can be improved by approximately a factor of three by using the proposed MG solver.First,the accuracy of the MG solver is validated through a test on a synthetic model.Next,the numerical performance of the inversion algorithm is evaluated using this model.Finally,the inversion algorithm is applied to a field EM data collected at the Beiya gold polymetallic ore district.A 3D resistivity model is obtained,and the formation process of the metal ore is analyzed.
基金financially supported by the National Natural Science Foundation of China (Nos.22171165 and 22371170)Natural Science Foundation of Shandong Province (No.ZR2022MB080)Scientific and Technological Frontiers in Project of Henan Province(No.242102110192)。
文摘Co-assembling chiral molecules with achiral compounds via non-covalent interactions like areneperfluoroarene(AP) interactions offers an effective approach for fabricating chiral functional materials.Herein,chiral molecules L/D-PF1 and L/D-PF2 with pyrene groups were synthesized and its chiroptical properties upon co-assembly with achiral compound octafluoronaphthalene(OFN) through AP interaction were systemically studied.The co-assembly of L/D-PF1/OFN and L/D-PF2/OFN exhibited distinct chiroptical properties such as circular dichroism(CD) and circularly polarized luminescence(CPL) signals.Chirality transfer from the chirality center of L/D-PF1 and L/D-PF2 to the achiral OFN and chiral amplification were successfully achieved.Besides,no significant CPL signal was observed in the self-assembly of L/DPF1 or L/D-PF2 while co-assembly with OFN exhibited obvious CPL amplification induced by AP interaction.Notably,a reversal CD signal and CPL signal could be observed in L/D-PF2/OFN when the molar ratio changed from 1:1 to 1:2 while not found in L/D-PF1/OFN,indicating that that minor structural changes of molecules could cause large changes in assembly.In addition,a series of computational calculations were conducted to verify the AP interaction between L-PF1/L-PF2 and OFN.This work demonstrated that arene-perfluoroarene interaction could drive chiral transfer,chiral amplification and chiral inversion and provided a new method for the preparation of chiroptical materials.
基金Supported by the China National Science and Technology Major Project(2025ZD1407000)PetroChina Science and Technology Major Project(2023ZZ22)。
文摘To address the challenges of complex fluvial sandbody distribution and difficult remaining oil recovery in mature continental oilfields,this study focuses on key issues in reservoir identification such as ambiguous narrow-channel boundaries and subdivision of multi-stage superimposed sandbodies.Taking the Upper Cretaceous continental sandstone in the Sazhong Oilfield of the Daqing Placanticline as an example,a technical system integrating OVT high-resolution processing,multi-attribute fusion,and varible-scale inversion was developed to establish a complete workflow from seismic processing to reservoir prediction and remaining oil recovery.The following results are obtained.First,the Offset Vector Tile(OVT)seismic processing technology is extended,for the first time,from fracture imaging to sandbody prediction,in order to address the weak seismic responses from boundaries of narrow and thin sandbodies.A geology-oriented OVT partitioning method is developed to significantly improve the imaging accuracy,enabling identification of channel sandbodies as narrow as 50 m.Second,an amplitude-coherence dual-attribute fusion method is proposed for predicting narrow channel boundaries between wells.Constrained by a sedimentary unit-level sequence chronostratigraphic framework,this method accurately delineates 800-2000 m long subaqueous distributary channels with bifurcation-convergence features.Third,considering the superimposition of multi-stage channels,a three-level variable-scale stratigraphic model(sandstone groups,sublayers,sedimentary units)is constructed to overcome single-scale modeling limitations,successfully characterizing key sedimentary features like meandering river“cut-offs”through 3D seismic inversion.Based on these advances,a direct link between seismic prediction and remaining oil recovery is established.The horizontal wells deployed using narrow-channel predictions encountered oil-bearing sandstones in the horizontal section by 97%,and achieved initial daily production of 12.5 t per well.Precise identification of individual channel boundaries within 17 composite sandbodies guided recovery processes in 135 wells,yielding an average daily increase of 2.8 t per well and a cumulative increase of 13.6×10^(4)t.
基金Supported by the Ministry of Education U40 Program(ZYGXONJSKYCXNLZCXM-E19)National Natural Science Foundation of China(52574078)。
文摘The forward model of optical fiber strain induced by fractures,together with the associated model resolution matrix,is used to demonstrate the interpretability of fracture parameters once the fracture intersects the fiber.A regularized inversion framework for fracture parameters is established to evaluate the influence of measured data quality on the accuracy of iterative regularized inversion.An interpretation approach for both fracture width and height is proposed,and the synthetic forward data with measurement error and field examples are employed to validate the accuracy of the simultaneous inversion of fracture width and height.The results indicate that,after the fracture contacts the fiber,the strain response is strongly sensitive only to the fracture parameters at the intersection location,whereas the interpretability of parameters at other locations remains limited.The iterative regularized inversion method effectively suppresses the impact of measurement error and exhibits high computational efficiency,showing clear advantages for inversion applications.When incorporating the first-order regularization with a Neumann boundary constraint on the tip width,the inverted fracture-width distribution becomes highly sensitive to fracture height;thus,combined with a bisection strategy,simultaneous inversion of fracture width and height can be achieved.Examination using the model resolution matrix,noisy synthetic data,and field data confirms that the iterative regularized inversion model for fracture width and height provides high interpretive accuracy and can be applied to the calculation and analysis of fracture width,fracture height,net pressure and other parameters.
基金supported by Fundamental Research Funds for the Central Universities,CHD300102264715National Key Research and Development Program of China under Grant 2021YFA0716902Natural Science Basic Research Program of Shaanxi 2024JCYBMS-199。
文摘Young's modulus and Poisson's ratio are crucial parameters for reservoir characterization and rock brittleness evaluation.Conventional methods often rely on indirect computation or approximations of the Zoeppritz equations to estimate Young's modulus,which can introduce cumulative errors and reduce the accuracy of inversion results.To address these issues,this paper introduces the analytical solution of the Zoeppritz equation into the inversion process.The equation is re-derived and expressed in terms of Young's modulus,Poisson's ratio,and density.Within the Bayesian framework,we construct an objective function for the joint inversion of PP and PS waves.Traditional gradient-based algorithms often suffer from low precision and the computational complexity.In this study,we address limitations of conventional approaches related to low precision and complicated code by using Circle chaotic mapping,Levy flights,and Gaussian mutation to optimize the quantum particle swarm optimization(QPSO),named improved quantum particle swarm optimization(IQPSO).The IQPSO demonstrates superior global optimization capabilities.We test the proposed inversion method with both synthetic and field data.The test results demonstrate the proposed method's feasibility and effectiveness,indicating an improvement in inversion accuracy over traditional methods.
基金support of China national key project 41904130 and key research project 041020080060.
文摘With the development of exploration of oil and gas resources,the requirements for seismic inversion results are getting more accurate.In particular,it is hoped that the distribution patterns of oil and gas reservoirs can be finely characterized,and the seismic inversion results can clearly characterize the location of stratigraphic boundaries and meet the needs of accurate geological description.Specifically,for pre-stack AVO inversion,it is required to be able to distinguish smaller geological targets in the depth or time domain,and clearly depict the vertical boundaries of the geological objects.In response to the above requirements,we introduce the preprocessing regularization of the adaptive edge-preserving smooth filter into the pre-stack AVO elastic parameter inversion to clearly invert the position of layer boundary and improve the accuracy of the inversion results.
基金supported by the National Natural Science Foundation of China(No.41574130,41874143 and 41374134)the National Science and Technology Major Project of China(No.2016ZX05014-001-009)the Sichuan Provincial Youth Science&Technology Innovative Research Group Fund(No.2016TD0023)
文摘The extensive application of pre-stack depth migration has produced huge volumes of seismic data,which allows for the possibility of developing seismic inversions of reservoir properties from seismic data in the depth domain.It is difficult to estimate seismic wavelets directly from seismic data due to the nonstationarity of the data in the depth domain.We conduct a velocity transformation of seismic data to make the seismic data stationary and then apply the ridge regression method to estimate a constant seismic wavelet.The estimated constant seismic wavelet is constructed as a set of space-variant seismic wavelets dominated by velocities at different spatial locations.Incorporating the weighted superposition principle,a synthetic seismogram is generated by directly employing the space-variant seismic wavelets in the depth domain.An inversion workflow based on the model-driven method is developed in the depth domain by incorporating the nonlinear conjugate gradient algorithm,which avoids additional data conversions between the time and depth domains.The impedance inversions of the synthetic and field seismic data in the depth domain show good results,which demonstrates that seismic inversion in the depth domain is feasible.The approach provides an alternative for forward numerical analyses and elastic property inversions of depth-domain seismic data.It is advantageous for further studies concerning the stability,accuracy,and efficiency of seismic inversions in the depth domain.
基金supported by the National Science and Technology Major Project(No.2016ZX05047-002-001)
文摘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.
基金NSFC and Sinopec Joint Key Project (U1663207)the China Geology Survey Project (DD20160195)+2 种基金973 Program (2014CB239104)National Key S&T Projects (2017ZX05049002)China Postdoctoral Science Foundation for the financial support
文摘For a typical marine shale reservoir in the Jiaoshiba area, Sichuan Basin of China, P-impedance is sensitive for identifying lithology but not suitable for indicating good shale reservoirs. In comparison, density is an important quantity, which is sensitive for identifying the organic-rich mud shale from non-organic-rich mud shale. Due to the poor data quality and incidence angle range, density cannot be easily inverted by directly solving the ill-posed pre-stack seismic inversion in this area. Meanwhile, the traditional density regularizations implemented by directly using the more robust P-impedance inversion tend to be inaccurate for recovering density for this shale reservoir. In this paper, we combine the P-impedance and the minus uranium to construct the pseudo-P-impedance（PIp） at well locations. The PIp is observed to be sensitive for identifying organic-rich mud shale and has a good correlation with density in this area. We employ the PIp–density relation into the pre-stack inversion framework to estimate density. Three types of regularization are tested on both numerical and field data： These are no regularization, traditional regularization and the proposed approach. It is observed that the proposed method is better for recovering the density of organic-rich mud shale in the Jiaoshiba area.
文摘Incremental Nonlinear Dynamic Inversion(INDI)is a control approach that has gained popularity in flight control over the past decade.Besides the INDI law,several common additional components complement an INDI-based controller.This paper,the second part of a two-part series of surveys on INDI,aims to summarize the modern trends in INDI and its related components.Besides a comprehensive components specification,it addresses their most common challenges,compares different variants,and discusses proposed advances.Further important aspects of INDI are gain design,stability,and robustness.This paper also provides an overview of research conducted concerning these aspects.This paper is written in a tutorial style to familiarize researchers with the essential specifics and pitfalls of INDI and its components.At the same time,it can also serve as a reference for readers already familiar with INDI.
基金supported by the National Natural Science Foundation of China(No.42174170)。
文摘An innovative gradient inversion approach employing the natural element method within the framework of least square regularization was proposed to enhance the quantitative interpretation of self-potential(SP)data originating from mineral polarization.The results indicated that the natural element method effectively addressed the challenge of subdividing complex resistivity models and aided in the accurate forward calculation of SP.By applying this approach to synthetic SP data and lab-measured SP data associated with redox electrochemical half-cell reactions of iron−copper metal blocks within the geobattery model,the 3D fine structure of buried orebody models was successfully reconstructed and the spatial distribution of SP current sources was mapped.This study significantly contributes to understanding the quantitative relationship between the polarization process of metal deposits and their corresponding SP responses and provides a valuable reference for delineating metal deposits in both terrestrial and marine environments through SP surveys.
基金supposed by the National Nature Science Foundation of China(Grant No.42304131)the Natural Science Foundation of Heilongjiang Province(Grant No.LH2023D012)+1 种基金the Heilongjiang Postdoctoral Fund(Grant No.LBH-Z22092)the Basic Research Fund for Universities in Xinjiang Uygur Autonomous Region(Grant No.XJEDU2023P166)。
文摘The brittleness index(BI)is crucial for predicting engineering sweet spots and designing fracturing operations in shale oil reservoir exploration and development.Seismic amplitude variation with offset(AVO)inversion is commonly used to obtain the BI.Traditionally,velocity,density,and other parameters are firstly inverted,and the BI is then calculated,which often leads to accumulated errors.Moreover,due to the limited of well-log data in field work areas,AVO inversion typically faces the challenge of limited information,resulting in not high accuracy of BI derived by existing AVO inversion methods.To address these issues,we first derive an AVO forward approximation equation that directly characterizes the BI in P-wave reflection coefficients.Based on this,an intelligent AVO inversion method,which combines the advantages of traditional and intelligent approaches,for directly obtaining the BI is proposed.A TransUnet model is constructed to establish the strong nonlinear mapping relationship between seismic data and the BI.By incorporating a combined objective function that is constrained by both low-frequency parameters and training samples,the challenge of limited samples is effectively addressed,and the direct inversion of the BI is stably achieved.Tests on model data and applications on field data demonstrate the feasibility,advancement,and practicality of the proposed method.
基金supported by Heilongjiang Province Basic Research Business Expenses for Universities Heilongjiang University Special Fund Project (Grant No. 2023-KYYWF-1494)the Natural Science Foundation of Jiangxi Province (Grant No. 20212BAB213023)。
文摘As geological exploration conditions become increasingly complex, meeting the requirements of precise geological exploration necessitates the development of a controlled-source audio magnetotelluric (CSAMT) inversion method that considers anisotropy to improve the effectiveness of inversion accuracy and interpretation accuracy of data. This study is based on the 3D fi nite-diff erence forward modeling of axis anisotropy using the reciprocity theorem to calculate the Jacobian matrix by applying the search method to automatically search for the Lagrange operator. The aim is to establish inversion iteration equations to achieve the axis anisotropic Occam's 3D inversion of tensor CSAMT in data space. Further, we obtain an underground axis anisotropic 3D geoelectric model by inverting the impedance data of tensor CSAMT. Two synthetic data examples show that using the isotropic tensor CSAMT algorithm to directly invert data in anisotropic media can generate false anomalies, leading to incorrect geological interpretations. Meanwhile, the proposed anisotropic inversion algorithm can eff ectively improve the accuracy of data inversion in anisotropic media. Further, the inversion examples verify the eff ectiveness and stability of the algorithm.
基金supported by the National Natural Science Foundation of China(Grant No.42174126)the Alliance of International Science Organizations(ANSO)Project(Grant No.ANSO-CR-PP2022-04)+1 种基金the Deep Earth Probe and Mineral Resources Exploration National Science and Technology Major Project(Grant Nos.2024ZD1002206,2024ZD1002201)Key R&D Program of Xinjiang Uyghur Autonomous Region(Grant No.2024B03013-2)。
文摘The Pamir Plateau,at the northwestern margin of the Tibetan Plateau,is a key region for investigating continental collision and plateau uplifting.To probe its deep structure,we collected seismic data from 263 stations across 11 research projects.We applied cross-correlation to noise data and extracted surface wave dispersion data from cross-correlation functions.The extracted dispersion data were subsequently inverted using a 3-D transdimensional Bayesian inversion method(rj-3 DMcMC).The inversion result reveals several crustal low-velocity zones(LVZs).Their formation is likely related to crustal thickening,the exposure of gneiss domes,and thicker sedimentary sequences compared to surrounding areas.In the lower crust and upper mantle,the LVZs in southern Pamir and southeastern Karakoram evolve into high-velocity zones,which expand northeastward with increasing depth.This suggests northward underthrusting of the Indian Plate.We also analyzed the Moho using both the standard deviation of S-wave velocity and the S-wave velocity structure.Results show that significant variations in velocity standard deviation reliably delineate the Moho interface.
基金supported by the National Key Research and Development Program of China(2022YFA1404400)the National Natural Science Foundation of China(62122072,12174368,61705216,62405306)+4 种基金Anhui Provincial Department of Science and Technology(202203a07020020,18030801138)Anhui Provincial Natural Science Foundation(2308085QA21,2408085QF187)the USTC Research Funds of the Double First-Class Initiative(YD2090002015)the Institute of Artificial Intelligence at Hefei Comprehensive National Science Center(23YGXT005)the Fundamental Research Funds for the Central Universities(WK2090000083).
文摘Ultrasound computed tomography(USCT)is a noninvasive biomedical imaging modality that offers insights into acoustic properties such as the sound speed(SS)and acoustic attenuation(AA)of the human body,enhancing diagnostic accuracy and therapy planning.Full waveform inversion(FWI)is a promising USCT image reconstruction method that optimizes the parameter fields of a wave propagation model via gradient-based optimization.However,twodimensional FWI methods are limited by their inability to account for three-dimensional wave propagation in the elevation direction,resulting in image artifacts.To address this problem,we propose a three-dimensional time-domain full waveform inversion algorithm to reconstruct the SS and AA distributions on the basis of a fractional Laplacian wave equation,adjoint field formulation,and gradient descent optimization.Validated by two sets of simulations,the proposed algorithm has potential for generating high-resolution and quantitative SS and AA distributions.This approach holds promise for clinical USCT applications,assisting early disease detection,precise abnormality localization,and optimized treatment planning,thus contributing to better healthcare outcomes.
基金We would like to acknowledge the sponsorship of the National Natural Science Foundation of China(42004092,42030103,41974119)Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(Qingdao)(Grant No.2021QNLM020001-6)Young Elite Scientists Sponsorship Program by CAST(2021QNRC001).
文摘Classical multi-channel technology can significantly reduce the pre-stack seismic inversion uncertainty, especially for complex geology such as high dipping structures. However, due to the consideration of complex structure or reflection features, the existing multi-channel inversion methods have to adopt the highly time-consuming strategy of arranging seismic data trace-by-trace, limiting its wide application in pre-stack inversion. A fast pre-stack multi-channel inversion constrained by seismic reflection features has been proposed to address this issue. The key to our method is to re-characterize the reflection features to directly constrain the pre-stack inversion through a Hadamard product operator without rearranging the seismic data. The seismic reflection features can reflect the distribution of the stratum reflection interface, and we obtained them from the post-stack profile by searching the shortest local Euclidean distance between adjacent seismic traces. Instead of directly constructing a large-size reflection features constraint operator advocated by the conventional methods, through decomposing the reflection features along the vertical and horizontal direction at a particular sampling point, we have constructed a computationally well-behaved constraint operator represented by the vertical and horizontal partial derivatives. Based on the Alternating Direction Method of Multipliers (ADMM) optimization, we have derived a fast algorithm for solving the objective function, including Hadamard product operators. Compared with the conventional reflection features constrained inversion, the proposed method is more efficient and accurate, proved on the Overthrust model and a field data set.
基金funded by the National Key R&D Program of China(Grant No.2022YFC2903904)the National Natural Science Foundation of China(Grant Nos.51904057 and U1906208).
文摘Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this challenge,nonlinear stress boundaries for a numerical model are determined through regression analysis of a series of nonlinear coefficient matrices,which are derived from the bubbling method.Considering the randomness and flexibility of the bubbling method,a parametric study is conducted to determine recommended ranges for these parameters,including the standard deviation(σb)of bubble radii,the non-uniform coefficient matrix number(λ)for nonlinear stress boundaries,and the number(m)and positions of in situ stress measurement points.A model case study provides a reference for the selection of these parameters.Additionally,when the nonlinear in situ stress inversion method is employed,stress distortion inevitably occurs near model boundaries,aligning with the Saint Venant's principle.Two strategies are proposed accordingly:employing a systematic reduction of nonlinear coefficients to achieve high inversion accuracy while minimizing significant stress distortion,and excluding regions with severe stress distortion near the model edges while utilizing the central part of the model for subsequent simulations.These two strategies have been successfully implemented in the nonlinear in situ stress inversion of the Xincheng Gold Mine and have achieved higher inversion accuracy than the linear method.Specifically,the linear and nonlinear inversion methods yield root mean square errors(RMSE)of 4.15 and 3.2,and inversion relative errors(δAve)of 22.08%and 17.55%,respectively.Therefore,the nonlinear inversion method outperforms the traditional multiple linear regression method,even in the presence of a systematic reduction in the nonlinear stress boundaries.
基金supported by the Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(Qingdao)(Grant No.2021QNLM020001)the National Key R&D Program of China(Grant No.2019YFC0605503C)+2 种基金the Major Scientific and Technological Projects of China National Petroleum Corporation(CNPC)(Grant No.ZD2019-183-003)the National Outstanding Youth Science Foundation(Grant No.41922028)the National Innovation Group Project(Grant No.41821002).
文摘Based on waveform fitting,full waveform inversion(FWI)is an important inversion method with the ability to reconstruct multi-parameter models in high precision.However,the strong nonlinear equation used in FWI presents the following challenges,such as low convergence efficiency,high dependence on the initial model,and the energy imbalance in deep region of the inverted model.To solve these inherent problems,we develop a timedomain elastic FWI method based on gradient preconditioning with the following details:(1)the limited memory Broyden Fletcher Goldfarb Shanno method with faster convergence is adopted to im-prove the inversion stability;(2)a multi-scaled inversion strategy is used to alleviate the nonlinear inversion instead of falling into the local minimum;(3)in addition,the pseudo-Hessian preconditioned illumination operator is involved for preconditioning the parameter gradients to improve the illumination equilibrium degree of deep structures.Based on the programming implementation of the new method,a deep depression model with five diffractors is used for testing.Compared with the conventional elastic FWI method,the technique proposed by this study has better effectiveness and accuracy on the inversion effect and con-vergence,respectively.
基金the sponsorship of National Natural Science Foundation of China(42325403)Deep Earth Probe and Mineral Resources Exploration-National Science and Technology Major Project of China(2024ZD1004201)。
文摘Subsurface reservoirs commonly exhibit layered structures.Conventional amplitude variation with angle(AVA)inversion,which relies on the Zoeppritz equation and its approximations,often fails to accurately estimate elastic parameters because it assumes single-interface models and ignores multiple reflections and transmission losses.To address these limitations,this study proposes a novel prestack time-frequency domain joint inversion method that utilizes the reflection matrix method(RMM)as the forward operator.The RMM accurately simulates wave propagation in layered media,while the joint inversion framework minimizes the misfit between observed and synthetic data in both the time and frequency domains.By incorporating Bayesian theory to optimize the inversion process,the method effectively balances contributions from both time-domain waveforms and frequency-domain spectral information through a weighting factor.Tests on both synthetic data and field data demonstrate that the proposed method outperforms conventional AVA inversion and time-domain waveform inversion in accuracy and robustness.Furthermore,the method demonstrates good robustness against variations in initial models,random noise,and coherent noise interference.This study provides a practical and effective approach for high-precision reservoir characterization,with potential applications in complex layered media.
