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.展开更多
The classical Gauss-Jordan method for matrix inversion involves augmenting the matrix with a unit matrix and requires a workspace twice as large as the original matrix as well as computational operations to be perform...The classical Gauss-Jordan method for matrix inversion involves augmenting the matrix with a unit matrix and requires a workspace twice as large as the original matrix as well as computational operations to be performed on both the original and the unit matrix. A modified version of the method for performing the inversion without explicitly generating the unit matrix by replicating its functionality within the original matrix space for more efficient utilization of computational resources is presented in this article. Although the algorithm described here picks the pivots solely from the diagonal which, therefore, may not contain a zero, it did not pose any problem for the author because he used it to invert structural stiffness matrices which met this requirement. Techniques such as row/column swapping to handle off-diagonal pivots are also applicable to this method but are beyond the scope of this article.展开更多
A modified version of the Gauss-Jordan algorithm for performing In-Place matrix inversion without using an augmenting unit matrix was described in a previous article by the author. He had also developed several Struct...A modified version of the Gauss-Jordan algorithm for performing In-Place matrix inversion without using an augmenting unit matrix was described in a previous article by the author. He had also developed several Structural Engineering softwares during his career using that method as their analysis engine. He chose matrix inversion because it was suitable for in-core solution of large numbers of vectors for the same set of equations as encountered in structural analysis of moving, dynamic and seismic loadings. The purpose of this article is to provide its readers with its theoretical background and detailed computations of an In-Place matrix inversion task as well as a Visual Basic routine of the algorithm for direct incorporation into Visual Basic 6TM softwares and Visual Basic for ApplicationsTM macros in MS-ExcelTM spreadsheets to save them time and effort of software development.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
In this study,an inverse design framework was established to find lightweight honeycomb structures(HCSs)with high impact resistance.The hybrid HCS,composed of re-entrant(RE)and elliptical annular re-entrant(EARE)honey...In this study,an inverse design framework was established to find lightweight honeycomb structures(HCSs)with high impact resistance.The hybrid HCS,composed of re-entrant(RE)and elliptical annular re-entrant(EARE)honeycomb cells,was created by constructing arrangement matrices to achieve structural lightweight.The machine learning(ML)framework consisted of a neural network(NN)forward regression model for predicting impact resistance and a multi-objective optimization algorithm for generating high-performance designs.The surrogate of the local design space was initially realized by establishing the NN in the small sample dataset,and the active learning strategy was used to continuously extended the local optimal design until the model converged in the global space.The results indicated that the active learning strategy significantly improved the inference capability of the NN model in unknown design domains.By guiding the iteration direction of the optimization algorithm,lightweight designs with high impact resistance were identified.The energy absorption capacity of the optimal design reached 94.98%of the EARE honeycomb,while the initial peak stress and mass decreased by 28.85%and 19.91%,respectively.Furthermore,Shapley Additive Explanations(SHAP)for global explanation of the NN indicated a strong correlation between the arrangement mode of HCS and its impact resistance.By reducing the stiffness of the cells at the top boundary of the structure,the initial impact damage sustained by the structure can be significantly improved.Overall,this study proposed a general lightweight design method for array structures under impact loads,which is beneficial for the widespread application of honeycomb-based protective structures.展开更多
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.展开更多
基金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.
文摘The classical Gauss-Jordan method for matrix inversion involves augmenting the matrix with a unit matrix and requires a workspace twice as large as the original matrix as well as computational operations to be performed on both the original and the unit matrix. A modified version of the method for performing the inversion without explicitly generating the unit matrix by replicating its functionality within the original matrix space for more efficient utilization of computational resources is presented in this article. Although the algorithm described here picks the pivots solely from the diagonal which, therefore, may not contain a zero, it did not pose any problem for the author because he used it to invert structural stiffness matrices which met this requirement. Techniques such as row/column swapping to handle off-diagonal pivots are also applicable to this method but are beyond the scope of this article.
文摘A modified version of the Gauss-Jordan algorithm for performing In-Place matrix inversion without using an augmenting unit matrix was described in a previous article by the author. He had also developed several Structural Engineering softwares during his career using that method as their analysis engine. He chose matrix inversion because it was suitable for in-core solution of large numbers of vectors for the same set of equations as encountered in structural analysis of moving, dynamic and seismic loadings. The purpose of this article is to provide its readers with its theoretical background and detailed computations of an In-Place matrix inversion task as well as a Visual Basic routine of the algorithm for direct incorporation into Visual Basic 6TM softwares and Visual Basic for ApplicationsTM macros in MS-ExcelTM spreadsheets to save them time and effort of software development.
基金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 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.
基金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.
基金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.
基金the financial supports from National Key R&D Program for Young Scientists of China(Grant No.2022YFC3080900)National Natural Science Foundation of China(Grant No.52374181)+1 种基金BIT Research and Innovation Promoting Project(Grant No.2024YCXZ017)supported by Science and Technology Innovation Program of Beijing institute of technology under Grant No.2022CX01025。
文摘In this study,an inverse design framework was established to find lightweight honeycomb structures(HCSs)with high impact resistance.The hybrid HCS,composed of re-entrant(RE)and elliptical annular re-entrant(EARE)honeycomb cells,was created by constructing arrangement matrices to achieve structural lightweight.The machine learning(ML)framework consisted of a neural network(NN)forward regression model for predicting impact resistance and a multi-objective optimization algorithm for generating high-performance designs.The surrogate of the local design space was initially realized by establishing the NN in the small sample dataset,and the active learning strategy was used to continuously extended the local optimal design until the model converged in the global space.The results indicated that the active learning strategy significantly improved the inference capability of the NN model in unknown design domains.By guiding the iteration direction of the optimization algorithm,lightweight designs with high impact resistance were identified.The energy absorption capacity of the optimal design reached 94.98%of the EARE honeycomb,while the initial peak stress and mass decreased by 28.85%and 19.91%,respectively.Furthermore,Shapley Additive Explanations(SHAP)for global explanation of the NN indicated a strong correlation between the arrangement mode of HCS and its impact resistance.By reducing the stiffness of the cells at the top boundary of the structure,the initial impact damage sustained by the structure can be significantly improved.Overall,this study proposed a general lightweight design method for array structures under impact loads,which is beneficial for the widespread application of honeycomb-based protective structures.
基金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.
