The multi-parameter inverse scattering problem of elastic waveequation with single fre- quency is investigated within Bornapproximation. By use of a wideband measuring scheme in which bothtransmitters and receivers sc...The multi-parameter inverse scattering problem of elastic waveequation with single fre- quency is investigated within Bornapproximation. By use of a wideband measuring scheme in which bothtransmitters and receivers scan over the half-space surface, theformula of the scattering field of elastic wave is derived. Fourtypes of mode conversion of elastic wave(P→P, P→S, S→P, S→S)areseparated from the scattering field. These components containsufficient information for usto recon- struct the configuration ofthe density and Lame parameters of the medium.展开更多
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.展开更多
Accurate assessment of snowpack volumetric liquid water content and bulk density is essential for understanding snow hydrology,avalanche risk management,and monitoring cryosphere changes.This study presents a novel du...Accurate assessment of snowpack volumetric liquid water content and bulk density is essential for understanding snow hydrology,avalanche risk management,and monitoring cryosphere changes.This study presents a novel dual-parameter inversion framework that integrates synthetic electromagnetic modelling,dimensionality reduction,and machine learning algorithms to extract relative permittivity and log-resistivity from ground-penetrating radar(GPR)data.Traditional snowpack measurements are invasive,labor-intensive,and limited to point observations.To overcome these limitations,we developed a non-invasive,scalable,and data-driven framework that uses synthetic GPR datasets representing diverse snowpack conditions with variable moisture and density profiles.Synthetic 1D time series reflections(A-scans)are generated using finite-difference time-domain simulations in the state-of-the-art electromagnetic simulator gprMax.Principal component analysis(PCA)is applied to compress each A-scan while preserving key features,which significantly improved and enhanced the model training efficiency.Four machine learning models,including random forest,neural network,support vector machine,and eXtreme gradient boosting,are trained on PCA-reduced features.Among these,the neural network model achieved the best performance,with R^(2)>0.97 for permittivity and R 2>0.92 for resistivity.Gaussian noise(signal-to-noise ratio of 6 dB)is introduced to the synthetic data,and then targeted domain adaptation is employed to enhance generalization to field data.The framework is validated on two contrasting GPR transects in the Altay Mountains of the Chinese mainland,representing moist(T750)and wet(G125)snowpack conditions.The neural network model predictions are most consistent with the GPR derived estimates,Snowfork measurements,and snow pit data,achieving volumetric liquid water content deviation of≤1.5% and bulk density error within the range of 30-84 kg m^(-3).The results demonstrate that machine learning-based inversion,supported by realistic simulations and data augmentation enables scalable,non-invasive snowpack characterization with significant applications in hydrological forecasting,snow monitoring,and water resource management.展开更多
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.展开更多
The parameter reconstruction of strong-scattering media is a challenge for conventional full waveform inversion(FWI).Direct envelope inversion(DEI)is an effective method for large-scale and strongscattering structures...The parameter reconstruction of strong-scattering media is a challenge for conventional full waveform inversion(FWI).Direct envelope inversion(DEI)is an effective method for large-scale and strongscattering structures imaging without the need of low-frequency seismic data.However,the current DEI methods are all based on the acoustic approximation.Whereas,in real cases,seismic records are the combined effects of the subsurface multi-parameters.Therefore,the study of DEI in elastic media is necessary for the accurate inversion of strong-scattering structures,such as salt domes.In this paper,we propose an elastic direct envelope inversion(EDEI)method based on wave mode decomposition.We define the objective function of EDEI using multi-component seismic data and derive its gradient formulation.To reduce the coupling effects of multi-parameters,we introduce the wave mode decomposition method into the gradient calculation of EDEI.The update of Vp is primarily the contributions of decomposed P-waves.Two approaches on Vs gradient calculation are proposed,i.e.using the petrophysical relation and wave mode decomposition method.Finally,we test the proposed method on a layered salt model and the SEG/EAGE salt model.The results show that the proposed EDEI method can reconstruct reliable large-scale Vp and Vs models of strong-scattering salt structures.The successive elastic FWI can obtain high-precision inversion results of the strong-scattering salt model.The proposed method also has a good anti-noise performance in the moderate noise level.展开更多
It is now common practice to perform simultaneous traveltime inversion for the velocity field and the reflector geometry in reflection/refraction tomography, or the velocity field and the hypocenter locations in regio...It is now common practice to perform simultaneous traveltime inversion for the velocity field and the reflector geometry in reflection/refraction tomography, or the velocity field and the hypocenter locations in regional earthquake tomography, but seldom are all three classes of model parameters updated simultaneously. This is mainly due to the trade-off between the different types of model parameters and the lack of different seismic phases to constrain the model parameters. Using a spherical-coordinate ray tracing algorithm for first and later（primary reflected） arrival tracing algorithm in combination with a popular linearized inversion solver, it is possible to simultaneously recover the three classes of model parameters in regional or global tomographic studies. In this paper we incorporate the multistage irregular shortest-path ray tracing algorithm（in a spherical coordinate system） with a subspace inversion solver to formulate a simultaneous inversion algorithm for triple model parameters updating using direct and later arrival time information.Comparison tests for two sets of data（noise free and added noise） indicate that the new triple-class parameter inversion algorithm is capable of obtaining nearly the same results as the double-class parameter inversion scheme. Furthermore,the proposed multi-parameter type inversion method is not sensitive to a modest level of picking error in the traveltime data, and also performs well with a relatively large uncertainty in earthquake hypocentral locations. This shows it to be a feasible and promising approach in regional or global tomographic applications.展开更多
Seismic migration and inversion are closely related techniques to portray subsurface images and identify hydrocarbon reservoirs.Seismic migration aims at obtaining structural images of subsurface geologic discontinuit...Seismic migration and inversion are closely related techniques to portray subsurface images and identify hydrocarbon reservoirs.Seismic migration aims at obtaining structural images of subsurface geologic discontinuities.More specifically,seismic migration estimates the reflectivity function(stacked average reflectivity or pre-stack angle-dependent reflectivity)from seismic reflection data.On the other hand,seismic inversion quantitatively estimates the intrinsic rock properties of subsurface formulations.Such seismic inversion methods are applicable to detect hydrocarbon reservoirs that may exhibit lateral variations in the inverted parameters.Although there exist many differences,pre-stack seismic migration is similar with the first iteration of the general linearized seismic inversion.Usually,seismic migration and inversion techniques assume an acoustic or isotropic elastic medium.Unconventional reservoirs such as shale and tight sand formation have notable anisotropic property.We present a linearized waveform inversion(LWI)scheme for weakly anisotropic elastic media with vertical transversely isotropic(VTI)symmetry.It is based on two-way anisotropic elastic wave equation and simultaneously inverts for the localized perturbations(ΔVp_(0)/Vp_(0)/Vs_(0)/Vs_(0)/,Δ∈,Δδ)from the long-wavelength reference model.Our proposed VTI-elastic LWI is an iterative method that requires a forward and an adjoint operator acting on vectors in each iteration.We derive the forward Born approximation operator by perturbation theory and adjoint operator via adjoint-state method.The inversion has improved the quality of the images and reduces the multi-parameter crosstalk comparing with the adjoint-based images.We have observed that the multi-parameter crosstalk problem is more prominent in the inversion images for Thomsen anisotropy parameters.Especially,the Thomsen parameter is the most difficult to resolve.We also analyze the multi-parameter crosstalk using scattering radiation patterns.The linearized waveform inversion for VTI-elastic media presented in this article provides quantitative information of the rock properties that has the potential to help identify hydrocarbon reservoirs.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
In this article,the multi-parameters Mittag-Leffler function is studied in detail.As a consequence,a series of novel results such as the integral representation,series representation and Mellin transform to the above ...In this article,the multi-parameters Mittag-Leffler function is studied in detail.As a consequence,a series of novel results such as the integral representation,series representation and Mellin transform to the above function,are obtained.Especially,we associate the multi-parameters Mittag-Leffler function with two special functions which are the generalized Wright hypergeometric and the Fox’s-H functions.Meanwhile,some interesting integral operators and derivative operators of this function,are also discussed.展开更多
Multi-parameter quantum estimation has attracted considerable attention due to its broad applications.Due to the complexity of quantum dynamics,existing research places significant emphasis on estimating parameters in...Multi-parameter quantum estimation has attracted considerable attention due to its broad applications.Due to the complexity of quantum dynamics,existing research places significant emphasis on estimating parameters in time-independent Hamiltonians.Here,our work makes an effort to explore multi-parameter estimation with time-dependent Hamiltonians.In particular,we focus on the discrimination of two close frequencies of a magnetic field by using a single qubit.We optimize the quantum controls by employing both traditional optimization methods and reinforcement learning to improve the precision for estimating the frequencies of the two magnetic fields.In addition to the estimation precision,we also evaluate the robustness of the optimization schemes against the shift of the control parameters.The results demonstrate that the hybrid reinforcement learning approach achieves the highest estimation precision,and exhibits superior robustness.Moreover,a fundamental challenge in multi-parameter quantum estimation stems from the incompatibility of the optimal control strategies for different parameters.We demonstrate that the hybrid control strategies derived through numerical optimization remain effective in enhancing the precision of multi-parameter estimation in spite of the incompatibilities,thereby mitigating incompatibilities between control strategies on the estimation precision.Finally,we investigate the trade-offs in estimation precision among different parameters for different scenarios,revealing the inherent challenges in balancing the optimization of multiple parameters simultaneously and providing insights into the fundamental distinction between quantum single-parameter estimation and multi-parameter estimation.展开更多
Battery safety has emerged as a critical challenge for achieving carbon neutrality,driven by the increasing frequency of thermal runaway incidents in electric vehicles(EVs)and stationary energy storage systems(ESSs).C...Battery safety has emerged as a critical challenge for achieving carbon neutrality,driven by the increasing frequency of thermal runaway incidents in electric vehicles(EVs)and stationary energy storage systems(ESSs).Conventional battery monitoring technologies struggle to track multiple physicochemical parameters in real time,hindering early hazard detection.Embedded optical fiber sensors have gained prominence as a transformative solution for next-generation smart battery sensing,owing to their micrometer size,multiplexing capability,and electromagnetic immunity.However,comprehensive reviews focusing on their advancements in operando multi-parameter monitoring remain scarce,despite their critical importance for ensuring battery safety.To address this gap,this review first introduces a classification and the fundamental principles of advanced battery-oriented optical fiber sensors.Subsequently,it summarizes recent developments in single-parameter battery monitoring using optical fiber sensors.Building on this foundation,this review presents the first comprehensive analysis of multifunctional optical fiber sensing platforms capable of simultaneously tracking temperature,strain,pressure,refractive index,and monitoring battery aging.Targeted strategies are proposed to facilitate the practical development of this technology,including optimization of sensor integration techniques,minimizing sensor invasiveness,resolving the cross-sensitivity of fiber Bragg grating(FBG)through structural innovation,enhancing techno-economics,and combining with artificial intelligence(AI).By aligning academic research with industry requirements,this review provides a methodological roadmap for developing robust optical sensing systems to ensure battery safety in decarbonization-driven applications.展开更多
Fluid identification and anisotropic parameters characterization are crucial for shale reservoir exploration and development.However,the anisotropic reflection coefficient equation,based on the transverse isotropy wit...Fluid identification and anisotropic parameters characterization are crucial for shale reservoir exploration and development.However,the anisotropic reflection coefficient equation,based on the transverse isotropy with a vertical axis of symmetry(VTI)medium assumption,involves numerous parameters to be inverted.This complexity reduces its stability and impacts the accuracy of seismic amplitude variation with offset(AVO)inversion results.In this study,a novel anisotropic equation that includes the fluid term and Thomsen anisotropic parameters is rewritten,which reduces the equation's dimensionality and increases its stability.Additionally,the traditional Markov Chain Monte Carlo(MCMC)inversion algorithm exhibits a high rejection rate for random samples and relies on known parameter distributions such as the Gaussian distribution,limiting the algorithm's convergence and sample randomness.To address these limitations and evaluate the uncertainty of AVO inversion,the IADR-Gibbs algorithm is proposed,which incorporates the Independent Adaptive Delayed Rejection(IADR)algorithm with the Gibbs sampling algorithm.Grounded in Bayesian theory,the new algorithm introduces support points to construct a proposal distribution of non-parametric distribution and reselects the rejected samples according to the Delayed Rejection(DR)strategy.Rejected samples are then added to the support points to update the proposal distribution function adaptively.The equation rewriting method and the IADR-Gibbs algorithm improve the accuracy and robustness of AVO inversion.The effectiveness and applicability of the proposed method are validated through synthetic gather tests and practical data applications.展开更多
In Rayleigh wave exploration,the inversion of dispersion curves is a crucial step for obtaining subsurface stratigraphic information,characterized by its multi-parameter and multi-extremum nature.Local optimization al...In Rayleigh wave exploration,the inversion of dispersion curves is a crucial step for obtaining subsurface stratigraphic information,characterized by its multi-parameter and multi-extremum nature.Local optimization algorithms used in dispersion curve inversion are highly dependent on the initial model and are prone to being trapped in local optima,while classical global optimization algorithms often suffer from slow convergence and low solution accuracy.To address these issues,this study introduces the Osprey Optimization Algorithm(OOA),known for its strong global search and local exploitation capabilities,into the inversion of dispersion curves to enhance inversion performance.In noiseless theoretical models,the OOA demonstrates excellent inversion accuracy and stability,accurately recovering model parameters.Even in noisy models,OOA maintains robust performance,achieving high inversion precision under high-noise conditions.In multimode dispersion curve tests,OOA effectively handles higher modes due to its efficient global and local search capabilities,and the inversion results show high consistency with theoretical values.Field data from the Wyoming region in the United States and a landfill site in Italy further verify the practical applicability of the OOA.Comprehensive test results indicate that the OOA outperforms the Particle Swarm Optimization(PSO)algorithm,providing a highly accurate and reliable inversion strategy for dispersion curve inversion.展开更多
Machine learning(ML)techniques have emerged as powerful tools for improving the predictive capabilities of Reynolds-averaged Navier-Stokes(RANS)turbulence models in separated flows.This improvement is achieved by leve...Machine learning(ML)techniques have emerged as powerful tools for improving the predictive capabilities of Reynolds-averaged Navier-Stokes(RANS)turbulence models in separated flows.This improvement is achieved by leveraging complex ML models,such as those developed using field inversion and machine learning(FIML),to dynamically adjust the constants within the baseline RANS model.However,the ML models often overlook the fundamental calibrations of the RANS turbulence model.Consequently,the basic calibration of the baseline RANS model is disrupted,leading to a degradation in the accuracy,particularly in basic wall-attached flows outside of the training set.To address this issue,a modified version of the Spalart-Allmaras(SA)turbulence model,known as Rubber-band SA(RBSA),has been proposed recently.This modification involves identifying and embedding constraints related to basic wall-attached flows directly into the model.It is shown that no matter how the parameters of the RBSA model are adjusted as constants throughout the flow field,its accuracy in wall-attached flows remains unaffected.In this paper,we propose a new constraint for the RBSA model,which better safeguards the law of wall in extreme conditions where the model parameter is adjusted dramatically.The resultant model is called the RBSA-poly model.We then show that when combined with FIML augmentation,the RBSA-poly model effectively preserves the accuracy of simple wall-attached flows,even when the adjusted parameters become functions of local flow variables rather than constants.A comparative analysis with the FIML-augmented original SA model reveals that the augmented RBSA-poly model reduces error in basic wall-attached flows by 50%while maintaining comparable accuracy in trained separated flows.These findings confirm the effectiveness of utilizing FIML in conjunction with the RBSA model,offering superior accuracy retention in cardinal flows.展开更多
基金Foundation of Ph.D Program of the State Education Commission of China
文摘The multi-parameter inverse scattering problem of elastic waveequation with single fre- quency is investigated within Bornapproximation. By use of a wideband measuring scheme in which bothtransmitters and receivers scan over the half-space surface, theformula of the scattering field of elastic wave is derived. Fourtypes of mode conversion of elastic wave(P→P, P→S, S→P, S→S)areseparated from the scattering field. These components containsufficient information for usto recon- struct the configuration ofthe density and Lame parameters of the medium.
基金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.
基金supported by the National Key R&D Program of China(Grant Nos.2023YFC3008300&2023YFC3008305)the National Natural Science Foundation of China(Grant No.42172320)+1 种基金the Key Laboratory of Mountain Hazards and Engineering Resilience,Institute of Mountain Hazards and Environment,Chinese Academy of Sciences(Grant Nos.KLMHER-Z06&KLMHER-T07)the Science and Technology Research Program of Institute of Mountain Hazards and Environment,Chinese Academy of Sciences(Grant No.IMHE-CXTD.04).
文摘Accurate assessment of snowpack volumetric liquid water content and bulk density is essential for understanding snow hydrology,avalanche risk management,and monitoring cryosphere changes.This study presents a novel dual-parameter inversion framework that integrates synthetic electromagnetic modelling,dimensionality reduction,and machine learning algorithms to extract relative permittivity and log-resistivity from ground-penetrating radar(GPR)data.Traditional snowpack measurements are invasive,labor-intensive,and limited to point observations.To overcome these limitations,we developed a non-invasive,scalable,and data-driven framework that uses synthetic GPR datasets representing diverse snowpack conditions with variable moisture and density profiles.Synthetic 1D time series reflections(A-scans)are generated using finite-difference time-domain simulations in the state-of-the-art electromagnetic simulator gprMax.Principal component analysis(PCA)is applied to compress each A-scan while preserving key features,which significantly improved and enhanced the model training efficiency.Four machine learning models,including random forest,neural network,support vector machine,and eXtreme gradient boosting,are trained on PCA-reduced features.Among these,the neural network model achieved the best performance,with R^(2)>0.97 for permittivity and R 2>0.92 for resistivity.Gaussian noise(signal-to-noise ratio of 6 dB)is introduced to the synthetic data,and then targeted domain adaptation is employed to enhance generalization to field data.The framework is validated on two contrasting GPR transects in the Altay Mountains of the Chinese mainland,representing moist(T750)and wet(G125)snowpack conditions.The neural network model predictions are most consistent with the GPR derived estimates,Snowfork measurements,and snow pit data,achieving volumetric liquid water content deviation of≤1.5% and bulk density error within the range of 30-84 kg m^(-3).The results demonstrate that machine learning-based inversion,supported by realistic simulations and data augmentation enables scalable,non-invasive snowpack characterization with significant applications in hydrological forecasting,snow monitoring,and water resource management.
基金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.
基金financial support jointly provided by the National Key R&D Program of China under contract number 2019YFC0605503Cthe Major Projects during the 14th Five-year Plan period under contract number 2021QNLM020001+2 种基金the National Outstanding Youth Science Foundation under contract number 41922028the Funds for Creative Research Groups of China under contract number 41821002the Major Projects of CNPC under contract number ZD2019-183-003。
文摘The parameter reconstruction of strong-scattering media is a challenge for conventional full waveform inversion(FWI).Direct envelope inversion(DEI)is an effective method for large-scale and strongscattering structures imaging without the need of low-frequency seismic data.However,the current DEI methods are all based on the acoustic approximation.Whereas,in real cases,seismic records are the combined effects of the subsurface multi-parameters.Therefore,the study of DEI in elastic media is necessary for the accurate inversion of strong-scattering structures,such as salt domes.In this paper,we propose an elastic direct envelope inversion(EDEI)method based on wave mode decomposition.We define the objective function of EDEI using multi-component seismic data and derive its gradient formulation.To reduce the coupling effects of multi-parameters,we introduce the wave mode decomposition method into the gradient calculation of EDEI.The update of Vp is primarily the contributions of decomposed P-waves.Two approaches on Vs gradient calculation are proposed,i.e.using the petrophysical relation and wave mode decomposition method.Finally,we test the proposed method on a layered salt model and the SEG/EAGE salt model.The results show that the proposed EDEI method can reconstruct reliable large-scale Vp and Vs models of strong-scattering salt structures.The successive elastic FWI can obtain high-precision inversion results of the strong-scattering salt model.The proposed method also has a good anti-noise performance in the moderate noise level.
基金partially supported by the Doctoral Programming Research Fund of Higher Education, Chinese Ministry of Education (No. 20110205110010)
文摘It is now common practice to perform simultaneous traveltime inversion for the velocity field and the reflector geometry in reflection/refraction tomography, or the velocity field and the hypocenter locations in regional earthquake tomography, but seldom are all three classes of model parameters updated simultaneously. This is mainly due to the trade-off between the different types of model parameters and the lack of different seismic phases to constrain the model parameters. Using a spherical-coordinate ray tracing algorithm for first and later（primary reflected） arrival tracing algorithm in combination with a popular linearized inversion solver, it is possible to simultaneously recover the three classes of model parameters in regional or global tomographic studies. In this paper we incorporate the multistage irregular shortest-path ray tracing algorithm（in a spherical coordinate system） with a subspace inversion solver to formulate a simultaneous inversion algorithm for triple model parameters updating using direct and later arrival time information.Comparison tests for two sets of data（noise free and added noise） indicate that the new triple-class parameter inversion algorithm is capable of obtaining nearly the same results as the double-class parameter inversion scheme. Furthermore,the proposed multi-parameter type inversion method is not sensitive to a modest level of picking error in the traveltime data, and also performs well with a relatively large uncertainty in earthquake hypocentral locations. This shows it to be a feasible and promising approach in regional or global tomographic applications.
文摘Seismic migration and inversion are closely related techniques to portray subsurface images and identify hydrocarbon reservoirs.Seismic migration aims at obtaining structural images of subsurface geologic discontinuities.More specifically,seismic migration estimates the reflectivity function(stacked average reflectivity or pre-stack angle-dependent reflectivity)from seismic reflection data.On the other hand,seismic inversion quantitatively estimates the intrinsic rock properties of subsurface formulations.Such seismic inversion methods are applicable to detect hydrocarbon reservoirs that may exhibit lateral variations in the inverted parameters.Although there exist many differences,pre-stack seismic migration is similar with the first iteration of the general linearized seismic inversion.Usually,seismic migration and inversion techniques assume an acoustic or isotropic elastic medium.Unconventional reservoirs such as shale and tight sand formation have notable anisotropic property.We present a linearized waveform inversion(LWI)scheme for weakly anisotropic elastic media with vertical transversely isotropic(VTI)symmetry.It is based on two-way anisotropic elastic wave equation and simultaneously inverts for the localized perturbations(ΔVp_(0)/Vp_(0)/Vs_(0)/Vs_(0)/,Δ∈,Δδ)from the long-wavelength reference model.Our proposed VTI-elastic LWI is an iterative method that requires a forward and an adjoint operator acting on vectors in each iteration.We derive the forward Born approximation operator by perturbation theory and adjoint operator via adjoint-state method.The inversion has improved the quality of the images and reduces the multi-parameter crosstalk comparing with the adjoint-based images.We have observed that the multi-parameter crosstalk problem is more prominent in the inversion images for Thomsen anisotropy parameters.Especially,the Thomsen parameter is the most difficult to resolve.We also analyze the multi-parameter crosstalk using scattering radiation patterns.The linearized waveform inversion for VTI-elastic media presented in this article provides quantitative information of the rock properties that has the potential to help identify hydrocarbon reservoirs.
基金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.
文摘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 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.
基金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.
基金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.
基金Supported by The National Undergraduate Innovation Training Program(Grant No.202310290069Z).
文摘In this article,the multi-parameters Mittag-Leffler function is studied in detail.As a consequence,a series of novel results such as the integral representation,series representation and Mellin transform to the above function,are obtained.Especially,we associate the multi-parameters Mittag-Leffler function with two special functions which are the generalized Wright hypergeometric and the Fox’s-H functions.Meanwhile,some interesting integral operators and derivative operators of this function,are also discussed.
基金supported by the National Natural Science Foundation of China(Grant No.12075323)the Innovation Program for Quantum Science and Technology(Grant No.2021ZD0300702).
文摘Multi-parameter quantum estimation has attracted considerable attention due to its broad applications.Due to the complexity of quantum dynamics,existing research places significant emphasis on estimating parameters in time-independent Hamiltonians.Here,our work makes an effort to explore multi-parameter estimation with time-dependent Hamiltonians.In particular,we focus on the discrimination of two close frequencies of a magnetic field by using a single qubit.We optimize the quantum controls by employing both traditional optimization methods and reinforcement learning to improve the precision for estimating the frequencies of the two magnetic fields.In addition to the estimation precision,we also evaluate the robustness of the optimization schemes against the shift of the control parameters.The results demonstrate that the hybrid reinforcement learning approach achieves the highest estimation precision,and exhibits superior robustness.Moreover,a fundamental challenge in multi-parameter quantum estimation stems from the incompatibility of the optimal control strategies for different parameters.We demonstrate that the hybrid control strategies derived through numerical optimization remain effective in enhancing the precision of multi-parameter estimation in spite of the incompatibilities,thereby mitigating incompatibilities between control strategies on the estimation precision.Finally,we investigate the trade-offs in estimation precision among different parameters for different scenarios,revealing the inherent challenges in balancing the optimization of multiple parameters simultaneously and providing insights into the fundamental distinction between quantum single-parameter estimation and multi-parameter estimation.
基金the financial supports of the National Natural Science Foundation of China(No.52372200)a project supported by the State Key Laboratory of Mechanics and Control for Aerospace Structures(No.MCAS-S-0324G01)。
文摘Battery safety has emerged as a critical challenge for achieving carbon neutrality,driven by the increasing frequency of thermal runaway incidents in electric vehicles(EVs)and stationary energy storage systems(ESSs).Conventional battery monitoring technologies struggle to track multiple physicochemical parameters in real time,hindering early hazard detection.Embedded optical fiber sensors have gained prominence as a transformative solution for next-generation smart battery sensing,owing to their micrometer size,multiplexing capability,and electromagnetic immunity.However,comprehensive reviews focusing on their advancements in operando multi-parameter monitoring remain scarce,despite their critical importance for ensuring battery safety.To address this gap,this review first introduces a classification and the fundamental principles of advanced battery-oriented optical fiber sensors.Subsequently,it summarizes recent developments in single-parameter battery monitoring using optical fiber sensors.Building on this foundation,this review presents the first comprehensive analysis of multifunctional optical fiber sensing platforms capable of simultaneously tracking temperature,strain,pressure,refractive index,and monitoring battery aging.Targeted strategies are proposed to facilitate the practical development of this technology,including optimization of sensor integration techniques,minimizing sensor invasiveness,resolving the cross-sensitivity of fiber Bragg grating(FBG)through structural innovation,enhancing techno-economics,and combining with artificial intelligence(AI).By aligning academic research with industry requirements,this review provides a methodological roadmap for developing robust optical sensing systems to ensure battery safety in decarbonization-driven applications.
基金the sponsorship of the Key Technology for Geophysical Prediction of Ultra-Deep Carbonate Reservoirs(P24240)the National Natural Science Foundation of China(U24B2020)the National Science and Technology Major Project of China for New Oil and Gas Exploration and Development(Grant No.2024ZD1400102)。
文摘Fluid identification and anisotropic parameters characterization are crucial for shale reservoir exploration and development.However,the anisotropic reflection coefficient equation,based on the transverse isotropy with a vertical axis of symmetry(VTI)medium assumption,involves numerous parameters to be inverted.This complexity reduces its stability and impacts the accuracy of seismic amplitude variation with offset(AVO)inversion results.In this study,a novel anisotropic equation that includes the fluid term and Thomsen anisotropic parameters is rewritten,which reduces the equation's dimensionality and increases its stability.Additionally,the traditional Markov Chain Monte Carlo(MCMC)inversion algorithm exhibits a high rejection rate for random samples and relies on known parameter distributions such as the Gaussian distribution,limiting the algorithm's convergence and sample randomness.To address these limitations and evaluate the uncertainty of AVO inversion,the IADR-Gibbs algorithm is proposed,which incorporates the Independent Adaptive Delayed Rejection(IADR)algorithm with the Gibbs sampling algorithm.Grounded in Bayesian theory,the new algorithm introduces support points to construct a proposal distribution of non-parametric distribution and reselects the rejected samples according to the Delayed Rejection(DR)strategy.Rejected samples are then added to the support points to update the proposal distribution function adaptively.The equation rewriting method and the IADR-Gibbs algorithm improve the accuracy and robustness of AVO inversion.The effectiveness and applicability of the proposed method are validated through synthetic gather tests and practical data applications.
基金sponsored by China Geological Survey Project(DD20243193 and DD20230206508).
文摘In Rayleigh wave exploration,the inversion of dispersion curves is a crucial step for obtaining subsurface stratigraphic information,characterized by its multi-parameter and multi-extremum nature.Local optimization algorithms used in dispersion curve inversion are highly dependent on the initial model and are prone to being trapped in local optima,while classical global optimization algorithms often suffer from slow convergence and low solution accuracy.To address these issues,this study introduces the Osprey Optimization Algorithm(OOA),known for its strong global search and local exploitation capabilities,into the inversion of dispersion curves to enhance inversion performance.In noiseless theoretical models,the OOA demonstrates excellent inversion accuracy and stability,accurately recovering model parameters.Even in noisy models,OOA maintains robust performance,achieving high inversion precision under high-noise conditions.In multimode dispersion curve tests,OOA effectively handles higher modes due to its efficient global and local search capabilities,and the inversion results show high consistency with theoretical values.Field data from the Wyoming region in the United States and a landfill site in Italy further verify the practical applicability of the OOA.Comprehensive test results indicate that the OOA outperforms the Particle Swarm Optimization(PSO)algorithm,providing a highly accurate and reliable inversion strategy for dispersion curve inversion.
基金supported by the National Natural Science Foundation of China(Grant Nos.12388101,12372288,U23A2069,and 92152301).
文摘Machine learning(ML)techniques have emerged as powerful tools for improving the predictive capabilities of Reynolds-averaged Navier-Stokes(RANS)turbulence models in separated flows.This improvement is achieved by leveraging complex ML models,such as those developed using field inversion and machine learning(FIML),to dynamically adjust the constants within the baseline RANS model.However,the ML models often overlook the fundamental calibrations of the RANS turbulence model.Consequently,the basic calibration of the baseline RANS model is disrupted,leading to a degradation in the accuracy,particularly in basic wall-attached flows outside of the training set.To address this issue,a modified version of the Spalart-Allmaras(SA)turbulence model,known as Rubber-band SA(RBSA),has been proposed recently.This modification involves identifying and embedding constraints related to basic wall-attached flows directly into the model.It is shown that no matter how the parameters of the RBSA model are adjusted as constants throughout the flow field,its accuracy in wall-attached flows remains unaffected.In this paper,we propose a new constraint for the RBSA model,which better safeguards the law of wall in extreme conditions where the model parameter is adjusted dramatically.The resultant model is called the RBSA-poly model.We then show that when combined with FIML augmentation,the RBSA-poly model effectively preserves the accuracy of simple wall-attached flows,even when the adjusted parameters become functions of local flow variables rather than constants.A comparative analysis with the FIML-augmented original SA model reveals that the augmented RBSA-poly model reduces error in basic wall-attached flows by 50%while maintaining comparable accuracy in trained separated flows.These findings confirm the effectiveness of utilizing FIML in conjunction with the RBSA model,offering superior accuracy retention in cardinal flows.
