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.展开更多
Seismic inversion is one of the most widely used technologies for reservoir prediction. Many good results have been obtained but sometimes it fails to differentiate the lithologies and identify the fluids. However, se...Seismic inversion is one of the most widely used technologies for reservoir prediction. Many good results have been obtained but sometimes it fails to differentiate the lithologies and identify the fluids. However, seismic prestack elastic inversion based on rock physics modeling and analysis introduced in this paper is a significant method that can help seismic inversion and interpretation reach a new quantitative （or semi-quantitative） level from traditional qualitative interpretation. By doing rock physics modeling and forward perturbation analysis, we can quantitatively analyze the essential relationships between rock properties and seismic responses and try to find the sensitive elastic properties to the lithology, porosity, fluid type, and reservoir saturation. Finally, standard rock physics templates （RPT） can be built for specific reservoirs to guide seismic inversion interpretation results for reservoir characterization and fluids identification purpose. The gas sand distribution results of the case study in this paper proves that this method has unparalleled advantages over traditional post-stack methods, by which we can perform reservoir characterization and seismic data interpretation more quantitatively and efficiently.展开更多
In multi-component oil and gas exploration using ocean bottom nodes,converted wave data is rich in lithological and fracture information.One of the urgent problems to be solved is how to construct an accurate shear wa...In multi-component oil and gas exploration using ocean bottom nodes,converted wave data is rich in lithological and fracture information.One of the urgent problems to be solved is how to construct an accurate shear wave velocity model of the shallow sea bottom by leveraging the seismic wave information at the fluid-solid interface in the ocean,and improve the lateral resolution of marine converted wave data.Given that the dispersion characteristics of surface waves are sensitive to the S-wave velocity of subsurface media,and that Scholte surface waves,which propagate at the interface between liquid and solid media,exist in the data of marine oil and gas exploration,this paper proposes a Scholte wave inversion and modeling method based on oil and gas exploration using ocean bottom nodes.By using the method for calculating the Scholte wave dispersion spectrum based on the Bessel kernel function,the accuracy of dispersion spectrum analysis is improved,and more accurate dispersion curves are picked up.Through the adaptive weighted least squares Scholte wave dispersion inversion algorithm,the Scholte wave dispersion equation for liquid-solid media is solved,and the shear wave velocity model of the shallow sea bottom is calculated.Theoretical tests and applications of realdata have proven that this method can significantly improve the lateral resolution of converted wave data,provide high-quality data for subsequent inversion of marine multi-component oil and gas exploration data and reservoir reflection information,and contribute to the development of marine oil and gas exploration technology.展开更多
Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose...Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose a joint gravity and magnetic inversion methodfor two-layer models by concentrating on the relationship between the change of thicknessI and position of the middle layer and anomaly and discuss the effects of the key parameters. Model tests and application to field data show the validity of this method.展开更多
Surface-wave inversion is a powerful tool for revealing the Earth's internal structure.However,aside from shear-wave velocity(v_(S)),other parameters can influence the inversion outcomes,yet these have not been sy...Surface-wave inversion is a powerful tool for revealing the Earth's internal structure.However,aside from shear-wave velocity(v_(S)),other parameters can influence the inversion outcomes,yet these have not been systematically discussed.This study investigates the influence of various parameter assumptions on the results of surface-wave inversion,including the compressional and shear velocity ratio(v_(P)/v_(S)),shear-wave attenuation(Q_(S)),density(ρ),Moho interface,and sedimentary layer.We constructed synthetic models to generate dispersion data and compared the obtained results with different parameter assumptions with those of the true model.The results indicate that the v_(P)/v_(S) ratio,Q_(S),and density(ρ) have minimal effects on absolute velocity values and perturbation patterns in the inversion.Conversely,assumptions about the Moho interface and sedimentary layer significantly influenced absolute velocity values and perturbation patterns.Introducing an erroneous Mohointerface depth in the initial model of the inversion significantly affected the v_(S) model near that depth,while using a smooth initial model results in relatively minor deviations.The assumption on the sedimentary layer not only affects shallow structure results but also impacts the result at greater depths.Non-linear inversion methods outperform linear inversion methods,particularly for the assumptions of the Moho interface and sedimentary layer.Joint inversion with other data types,such as receiver functions or Rayleigh wave ellipticity,and using data from a broader period range or higher-mode surface waves,can mitigate these deviations.Furthermore,incorporating more accurate prior information can improve inversion results.展开更多
How to get the rapid and stable inversion results and reconstruct the clear subsurface resistivity structures is a focus problem in current magnetotelluric inversion. A stable solution of an ill-posed inverse problem ...How to get the rapid and stable inversion results and reconstruct the clear subsurface resistivity structures is a focus problem in current magnetotelluric inversion. A stable solution of an ill-posed inverse problem was obtained by the regularization methods in which some desired structures were imposed to stabilize the inverse problem. By the smoothness-constrained model and approximate sensitivity method, the stable subsurface resistivity structures were reconstructed. The synthetic examples show that the smoothness-constrained regularized inversion method is effective and can be reasonable to reconstruct three-dimensional subsurface resistivity structures.展开更多
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.展开更多
The interaction between a two-level atom and a single-mode field in the k-photon Jaynes-Cummings model （JCM） in the presence of the Stark shift and a Kerr medium is studied. All terms in the Hamiltonian, such as the...The interaction between a two-level atom and a single-mode field in the k-photon Jaynes-Cummings model （JCM） in the presence of the Stark shift and a Kerr medium is studied. All terms in the Hamiltonian, such as the single-mode field, its interaction with the atom, the contribution of the Stark shift and the Kerr medium effects are considered to be f-deformed. In particular, the effect of the initial state of the radiation field on the dynamical evolution of some physical properties such as atomic inversion and entropy squeezing are investigated by considering different initial field states （coherent, squeezed and thermal states）.展开更多
Theoretical analysis and practical observations show that fault dislocations can change the gravity field around the fault. Gravity changes which were caused by the repeated dislocations over a long period of time wer...Theoretical analysis and practical observations show that fault dislocations can change the gravity field around the fault. Gravity changes which were caused by the repeated dislocations over a long period of time were superimposed on the Bougeur gravity anomalies. These anomalies became the evidence of historical movement of fault as well as provide a way for the study of paleo earthquakes. This paper investigates inversion methods for the geological dislocation modeling of faults using the local Bouguer's gravity anomalies. To remove the effects of the irrelevant part of gravity anomalies to fault movements, we propose the robust nonlinear inversion method and set up the corresponding algorithm. Modeling examples indicate that the Marquardt's and Baye's least squares solutions depart from the true solution due to the attraction of gross errors in the data. The more seriously the data is contaminated, the more seriously the solutions are biased. In contrast, the proposed robust Marquardt's and Baye's inversion solutions can still maintain consistency with the solution without gross errors, even though 50 percent of the data is contaminated. This indicates that the proposed robust methods are effective. Using the proposed methods, we invert the geological dislocation models of the faults around the Erhai Lake in West Yunnan. The results show that the Northern Cangdong fault and the Erhai fault are normal dip slip faults with about 4 to 5 km dislocations; and that the Southern Cangdong fault has a less dip slip compared with the former two. A satisfactory fitting between the theoretical values of the inversion solution and the actual local gravity field is achievable.展开更多
The coal-bearing strata of the deep Upper Paleozoic in the GS Sag have high hydrocarbon potential. Because of the absence of seismic data, we use electromagnetic （MT） and gravity data jointly to delineate the distri...The coal-bearing strata of the deep Upper Paleozoic in the GS Sag have high hydrocarbon potential. Because of the absence of seismic data, we use electromagnetic （MT） and gravity data jointly to delineate the distribution of deep targets based on well logging and geological data. First, a preliminary geological model is established by using three-dimensional （3D） MT inversion results. Second, using the formation density and gravity anomalies, the preliminary geological model is modified by interactive inversion of the gravity data. Then, we conduct MT-constrained inversion based on the modified model to obtain an optimal geological model until the deviations at all stations are minimized. Finally, the geological model and a seismic profile in the middle of the sag is analysed. We determine that the deep reflections of the seismic profile correspond to the Upper Paleozoic that reaches thickness up to 800 m. The processing of field data suggests that the joint MT-gravity modeling and constrained inversion can reduce the multiple solutions for single geophysical data and thus improve the recognition of deep formations. The MT-constrained inversion is consistent with the geological features in the seismic section. This suggests that the joint MT and gravity modeling and constrained inversion can be used to delineate deep targets in similar basins.展开更多
Borehole-to-surface electrical imaging （BSEI） uses a line source and a point source to generate a stable electric field in the ground. In order to study the surface potential of anomalies, three-dimensional forward ...Borehole-to-surface electrical imaging （BSEI） uses a line source and a point source to generate a stable electric field in the ground. In order to study the surface potential of anomalies, three-dimensional forward modeling of point and line sources was conducted by using the finite-difference method and the incomplete Cholesky conjugate gradient （ICCG） method. Then, the damping least square method was used in the 3D inversion of the formation resistivity data. Several geological models were considered in the forward modeling and inversion. The forward modeling results suggest that the potentials generated by the two sources have different surface signatures. The inversion data suggest that the low- resistivity anomaly is outlined better than the high-resistivity anomaly. Moreover, when the point source is under the anomaly, the resistivity anomaly boundaries are better outlined than when using a line source.展开更多
As a primary parameter in the water quality model for shallow bays, the dispersion coefficient is traditionally determined with a trial-and-error method, which is time-consuming and requires much experience. In this p...As a primary parameter in the water quality model for shallow bays, the dispersion coefficient is traditionally determined with a trial-and-error method, which is time-consuming and requires much experience. In this paper, based on the measured data of chemical oxygen demand （COD）, the dispersion coefficient is calculated using an inversion method. In the process, the regularization method is applied to treat the ill-posedness, and an operator identity perturbation method is used to obtain the solu- tion. Using the model with an inverted dispersion coefficient, the distributions of COD, inorganic nitrogen （IN）, and inorganic phosphorus （IP） in Bohai Bay are predicted and compared with the measured data. The results indicate that the method is feasible and the inverted dispersion coefficient can be used to predict other pollutant distribution. This method may also be further extended to the inversion of other parameters in the water quality model.展开更多
Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well descri...Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks,which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block,and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization,which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models.The damped least squares method is employed in seismic traveltime inversion,which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model.展开更多
The Zhuxi tungsten deposit in Jiangxi Province,South China,contains a total W reserve of about 2.86 Mt at an average grade of 0.54 wt%WO3,representing the largest W deposit in the world.Numerous studies on the metallo...The Zhuxi tungsten deposit in Jiangxi Province,South China,contains a total W reserve of about 2.86 Mt at an average grade of 0.54 wt%WO3,representing the largest W deposit in the world.Numerous studies on the metallogeny of the deposit have included its timing,the ore-controlling structures and sedimentary host rocks and their implications for mineral exploration.However,the deep nappe structural style of Taqian-Fuchun metallogenic belt that hosts the W deposit,and the spatial shape and scale of deeply concealed intrusions and their sedimentary host rocks are still poorly defined,which seriously restricts the discovery of new deposits at depth and in surrounding areas of the W deposit.Modern 3 D geological modeling is an important tool for the exploration of concealed orebodies,especially in brownfield environments.There are obvious density contrast and weak magnetic contrast in the ore-controlling strata and granite at the periphery of the deposit,which lays a physical foundation for solving the 3 D spatial problems of the ore-controlling geological body in the deep part of the study area through gravity and magnetic modeling.Gravity data(1:50000)and aeromagnetic data(1:50000)from the latest geophysical surveys of 2016-2018 have been used,firstly,to carry out a potential field separation to obtain residual anomalies for gravity and magnetic interactive inversion.Then,on the basis of the analysis of the relationship between physical properties and lithology,under the constraints of surface geology and borehole data,human-computer interactive gravity and magnetic inversion for 18 cross-sections were completed.Finally,the 3 D geological model of the Zhuxi tungsten deposit and its periphery have been established through these 18 sections,and the spatial shape of the intrusions and strata with a depth of 5 km underground were obtained,initially realizing―transparency‖for ore-controlling bodies.According the analysis of the geophysical,geochemical,and geological characteristics of the Zhuxi tungsten deposit,we discern three principles for prospecting and prediction in the research area,and propose five new exploration targets in its periphery.展开更多
The maritime tropospheric duct is a low-altitude anomalous refractivity structure over the ocean surface,and it can significantly affect the performance of many shore-based/shipboard radar and communication systems.We...The maritime tropospheric duct is a low-altitude anomalous refractivity structure over the ocean surface,and it can significantly affect the performance of many shore-based/shipboard radar and communication systems.We propose the idea that maritime tropospheric ducts can be retrieved from ocean forward-scattered low-elevation global positioning system（GPS） signals.Retrieval is accomplished by matching the measured power patterns of the signals to those predicted by the forward propagation model as a function of the modified refractivity profile.On the basis of a parabolic equation method and bistatic radar equation,we develop such a forward model for computing the trapped propagation characteristics of an ocean forward-scattered GPS signal within a tropospheric duct.A new GPS scattering initial field is defined for this model to start the propagation modeling.A preliminary test on the performance of this model is conducted using measured data obtained from a 2009-experiment in the South China Sea.Results demonstrate that this model can predict GPS propagation characteristics within maritime tropospheric ducts and serve as a forward model for duct inversion.展开更多
Seismic impedance inversion is an important technique for structure identification and reservoir prediction.Model-based and data-driven impedance inversion are the commonly used inversion methods.In practice,the geoph...Seismic impedance inversion is an important technique for structure identification and reservoir prediction.Model-based and data-driven impedance inversion are the commonly used inversion methods.In practice,the geophysical inversion problem is essentially an ill-posedness problem,which means that there are many solutions corresponding to the same seismic data.Therefore,regularization schemes,which can provide stable and unique inversion results to some extent,have been introduced into the objective function as constrain terms.Among them,given a low-frequency initial impedance model is the most commonly used regularization method,which can provide a smooth and stable solution.However,this model-based inversion method relies heavily on the initial model and the inversion result is band limited to the effective frequency bandwidth of seismic data,which cannot effectively improve the seismic vertical resolution and is difficult to be applied to complex structural regions.Therefore,we propose a data-driven approach for high-resolution impedance inversion based on the bidirectional long short-term memory recurrent neural network,which regards seismic data as time-series rather than image-like patches.Compared with the model-based inversion method,the data-driven approach provides higher resolution inversion results,which demonstrates the effectiveness of the data-driven method for recovering the high-frequency components.However,judging from the inversion results for characterization the spatial distribution of thin-layer sands,the accuracy of high-frequency components is difficult to guarantee.Therefore,we add the model constraint to the objective function to overcome the shortages of relying only on the data-driven schemes.First,constructing the supervisor1 based on the bidirectional long short-term memory recurrent neural network,which provides the predicted impedance with higher resolution.Then,convolution constraint as supervisor2 is introduced into the objective function to guarantee the reliability and accuracy of the inversion results,which makes the synthetic seismic data obtained from the inversion result consistent with the input data.Finally,we test the proposed scheme based on the synthetic and field seismic data.Compared to model-based and purely data-driven impedance inversion methods,the proposed approach provides more accurate and reliable inversion results while with higher vertical resolution and better spatial continuity.The inversion results accurately characterize the spatial distribution relationship of thin sands.The model tests demonstrate that the model-constrained and data-driven impedance inversion scheme can effectively improve the thin-layer structure characterization based on the seismic data.Moreover,tests on the oil field data indicate the practicality and adaptability of the proposed method.展开更多
Flow units(FU)rock typing is a common technique for characterizing reservoir flow behavior,producing reliable porosity and permeability estimation even in complex geological settings.However,the lateral extrapolation ...Flow units(FU)rock typing is a common technique for characterizing reservoir flow behavior,producing reliable porosity and permeability estimation even in complex geological settings.However,the lateral extrapolation of FU away from the well into the whole reservoir grid is commonly a difficult task and using the seismic data as constraints is rarely a subject of study.This paper proposes a workflow to generate numerous possible 3D volumes of flow units,porosity and permeability below the seismic resolution limit,respecting the available seismic data at larger scales.The methodology is used in the Mero Field,a Brazilian presalt carbonate reservoir located in the Santos Basin,who presents a complex and heterogenic geological setting with different sedimentological processes and diagenetic history.We generated metric flow units using the conventional core analysis and transposed to the well log data.Then,given a Markov chain Monte Carlo algorithm,the seismic data and the well log statistics,we simulated acoustic impedance,decametric flow units(DFU),metric flow units(MFU),porosity and permeability volumes in the metric scale.The aim is to estimate a minimum amount of MFU able to calculate realistic scenarios porosity and permeability scenarios,without losing the seismic lateral control.In other words,every porosity and permeability volume simulated produces a synthetic seismic that match the real seismic of the area,even in the metric scale.The achieved 3D results represent a high-resolution fluid flow reservoir modelling considering the lateral control of the seismic during the process and can be directly incorporated in the dynamic characterization workflow.展开更多
Based on the pseudo-analytical equation of electromagnetic log for layered formation,an optimal boundary match method is proposed to adaptively truncate the encountered formation structures.An efficient integral metho...Based on the pseudo-analytical equation of electromagnetic log for layered formation,an optimal boundary match method is proposed to adaptively truncate the encountered formation structures.An efficient integral method is put forward to significantly accelerate the convergence of Sommerfeld integral.By asymptotically approximating and subtracting the first reflection/transmission waves from the scattered field,the new Sommerfeld integral method has addressed difficulties encountered by the traditional digital filtering method,such as low computational precision and limited operating range,and realized the acceleration of the computation speed of logging-while-drilling electromagnetic measurements(LWD EM).By making use of the priori information from the offset/pilot wells and interactively adjusting the formation model,the optimum initial guesses of the inversion model is determined in order to predict the nearby formation boundaries.The gradient optimization algorithm is developed and an interactive inversion system for the LWD EM data from the horizontal wells is established.The inverted results of field data demonstrated that the real-time interactive inversion method is capable of providing the accurate boundaries of layers around the wellbore from the LWD EM,and it will benefit the wellbore trajectory optimization and reservoir interpretation.展开更多
Research on quantitative models of suspended sediment concentration （SSC） using remote sensing technology is very important to understand the scouting and siltation variation in harbors and water channels. Based onl...Research on quantitative models of suspended sediment concentration （SSC） using remote sensing technology is very important to understand the scouting and siltation variation in harbors and water channels. Based onlaboratory study of the relationship between different suspended sediment concentrations and reflectance spectra measured synchronously, quantitative inversion models of SSC based on single factor, band ratio and sediment parameter were developed, which provides an effective method to retrieve the SSC from satellite images. Results show that the bl （430-500nm） and b3 （670-735nm） are the optimal wavelengths for the estimation of lower SSC and the b4 （780-835nm） is the optimal wavelength to estimate the higher SSC. Furthermore the band ratio B2/B3 can be used to simulate the variation of lower SSC better and the B4/B1 to estimate the higher SSC accurately. Also the inversion models developed by sediment parameters of higher and lower SSCs can get a relatively higher accuracy than the single factor and band ratio models.展开更多
There is not much correlation between sand content and P-wave impedance in faulted basins,and the prediction results cannot be guaranteed. Due to sedimentary facies control,there is inconsistency between sand content ...There is not much correlation between sand content and P-wave impedance in faulted basins,and the prediction results cannot be guaranteed. Due to sedimentary facies control,there is inconsistency between sand content and low frequency trend of P-wave impedance,causing problems for seismic inversion modeling,which directly affects final seismic inversion results. The acoustic impedance increases with burial depth. When the same layers of sand and shale formation endure different compaction,the acoustic impedance values will be different. Therefore,The seismic inversion modeling of faulted basin is different from that of conventional basins.The authors built an inversion model using uncompact sandstone percentage,directly compensating for the low frequency trend of the inversion model. In addition,the model's intermediate frequency is similar to P-wave impedance,ensuring that the inversion results are converged correctly. The final results of the inversion can be used directly as sandstone percentage. The aforementioned method was applied to Liaohe Beach in the Bijialing region and obtained optimistic lithology inversion effects: the results were correlated well with the depositional facies map and with lithology in the well borehole. The inversion results can be used to define the sand body horizontally and can separate sand bodies vertically,which is very difficult on conventional seismic section.Therefore the inversion results played an important role in reservoir prediction.展开更多
基金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.
文摘Seismic inversion is one of the most widely used technologies for reservoir prediction. Many good results have been obtained but sometimes it fails to differentiate the lithologies and identify the fluids. However, seismic prestack elastic inversion based on rock physics modeling and analysis introduced in this paper is a significant method that can help seismic inversion and interpretation reach a new quantitative （or semi-quantitative） level from traditional qualitative interpretation. By doing rock physics modeling and forward perturbation analysis, we can quantitatively analyze the essential relationships between rock properties and seismic responses and try to find the sensitive elastic properties to the lithology, porosity, fluid type, and reservoir saturation. Finally, standard rock physics templates （RPT） can be built for specific reservoirs to guide seismic inversion interpretation results for reservoir characterization and fluids identification purpose. The gas sand distribution results of the case study in this paper proves that this method has unparalleled advantages over traditional post-stack methods, by which we can perform reservoir characterization and seismic data interpretation more quantitatively and efficiently.
基金financially supported by the Scientific Research and Technology Development Project of China National Petroleum Corporation(No.2021ZG02)titled"Development of Seismic Data Processing Software for Ocean Nodes(OBN)"。
文摘In multi-component oil and gas exploration using ocean bottom nodes,converted wave data is rich in lithological and fracture information.One of the urgent problems to be solved is how to construct an accurate shear wave velocity model of the shallow sea bottom by leveraging the seismic wave information at the fluid-solid interface in the ocean,and improve the lateral resolution of marine converted wave data.Given that the dispersion characteristics of surface waves are sensitive to the S-wave velocity of subsurface media,and that Scholte surface waves,which propagate at the interface between liquid and solid media,exist in the data of marine oil and gas exploration,this paper proposes a Scholte wave inversion and modeling method based on oil and gas exploration using ocean bottom nodes.By using the method for calculating the Scholte wave dispersion spectrum based on the Bessel kernel function,the accuracy of dispersion spectrum analysis is improved,and more accurate dispersion curves are picked up.Through the adaptive weighted least squares Scholte wave dispersion inversion algorithm,the Scholte wave dispersion equation for liquid-solid media is solved,and the shear wave velocity model of the shallow sea bottom is calculated.Theoretical tests and applications of realdata have proven that this method can significantly improve the lateral resolution of converted wave data,provide high-quality data for subsequent inversion of marine multi-component oil and gas exploration data and reservoir reflection information,and contribute to the development of marine oil and gas exploration technology.
基金Supported by the National Natural Science Foundation of China(Grant No.40674063)National Hi-tech Research and Development Program of China(863Program)(Grant No.2006AA09Z311)
文摘Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose a joint gravity and magnetic inversion methodfor two-layer models by concentrating on the relationship between the change of thicknessI and position of the middle layer and anomaly and discuss the effects of the key parameters. Model tests and application to field data show the validity of this method.
基金supported by the Special Fund of the Institute of Geophysics, China Earthquake Administration (No. DQJB21B32)the National Key R&D Program of China (No. 2022YFF0800601)。
文摘Surface-wave inversion is a powerful tool for revealing the Earth's internal structure.However,aside from shear-wave velocity(v_(S)),other parameters can influence the inversion outcomes,yet these have not been systematically discussed.This study investigates the influence of various parameter assumptions on the results of surface-wave inversion,including the compressional and shear velocity ratio(v_(P)/v_(S)),shear-wave attenuation(Q_(S)),density(ρ),Moho interface,and sedimentary layer.We constructed synthetic models to generate dispersion data and compared the obtained results with different parameter assumptions with those of the true model.The results indicate that the v_(P)/v_(S) ratio,Q_(S),and density(ρ) have minimal effects on absolute velocity values and perturbation patterns in the inversion.Conversely,assumptions about the Moho interface and sedimentary layer significantly influenced absolute velocity values and perturbation patterns.Introducing an erroneous Mohointerface depth in the initial model of the inversion significantly affected the v_(S) model near that depth,while using a smooth initial model results in relatively minor deviations.The assumption on the sedimentary layer not only affects shallow structure results but also impacts the result at greater depths.Non-linear inversion methods outperform linear inversion methods,particularly for the assumptions of the Moho interface and sedimentary layer.Joint inversion with other data types,such as receiver functions or Rayleigh wave ellipticity,and using data from a broader period range or higher-mode surface waves,can mitigate these deviations.Furthermore,incorporating more accurate prior information can improve inversion results.
基金Project(20110162120064)supported by Higher School Doctor Subject Special Scientific Research Foundation of ChinaProject(10JJ6059)supported by the Natural Science Foundation of Hunan Province,China
文摘How to get the rapid and stable inversion results and reconstruct the clear subsurface resistivity structures is a focus problem in current magnetotelluric inversion. A stable solution of an ill-posed inverse problem was obtained by the regularization methods in which some desired structures were imposed to stabilize the inverse problem. By the smoothness-constrained model and approximate sensitivity method, the stable subsurface resistivity structures were reconstructed. The synthetic examples show that the smoothness-constrained regularized inversion method is effective and can be reasonable to reconstruct three-dimensional subsurface resistivity structures.
文摘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.
文摘The interaction between a two-level atom and a single-mode field in the k-photon Jaynes-Cummings model （JCM） in the presence of the Stark shift and a Kerr medium is studied. All terms in the Hamiltonian, such as the single-mode field, its interaction with the atom, the contribution of the Stark shift and the Kerr medium effects are considered to be f-deformed. In particular, the effect of the initial state of the radiation field on the dynamical evolution of some physical properties such as atomic inversion and entropy squeezing are investigated by considering different initial field states （coherent, squeezed and thermal states）.
文摘Theoretical analysis and practical observations show that fault dislocations can change the gravity field around the fault. Gravity changes which were caused by the repeated dislocations over a long period of time were superimposed on the Bougeur gravity anomalies. These anomalies became the evidence of historical movement of fault as well as provide a way for the study of paleo earthquakes. This paper investigates inversion methods for the geological dislocation modeling of faults using the local Bouguer's gravity anomalies. To remove the effects of the irrelevant part of gravity anomalies to fault movements, we propose the robust nonlinear inversion method and set up the corresponding algorithm. Modeling examples indicate that the Marquardt's and Baye's least squares solutions depart from the true solution due to the attraction of gross errors in the data. The more seriously the data is contaminated, the more seriously the solutions are biased. In contrast, the proposed robust Marquardt's and Baye's inversion solutions can still maintain consistency with the solution without gross errors, even though 50 percent of the data is contaminated. This indicates that the proposed robust methods are effective. Using the proposed methods, we invert the geological dislocation models of the faults around the Erhai Lake in West Yunnan. The results show that the Northern Cangdong fault and the Erhai fault are normal dip slip faults with about 4 to 5 km dislocations; and that the Southern Cangdong fault has a less dip slip compared with the former two. A satisfactory fitting between the theoretical values of the inversion solution and the actual local gravity field is achievable.
基金supported by the National Science and Technology Major Project(No.2016ZX05018006)the National Key Research Development Program(No.2016YFC0601104)the National Natural Science Foundation of China(No.41472136)
文摘The coal-bearing strata of the deep Upper Paleozoic in the GS Sag have high hydrocarbon potential. Because of the absence of seismic data, we use electromagnetic （MT） and gravity data jointly to delineate the distribution of deep targets based on well logging and geological data. First, a preliminary geological model is established by using three-dimensional （3D） MT inversion results. Second, using the formation density and gravity anomalies, the preliminary geological model is modified by interactive inversion of the gravity data. Then, we conduct MT-constrained inversion based on the modified model to obtain an optimal geological model until the deviations at all stations are minimized. Finally, the geological model and a seismic profile in the middle of the sag is analysed. We determine that the deep reflections of the seismic profile correspond to the Upper Paleozoic that reaches thickness up to 800 m. The processing of field data suggests that the joint MT-gravity modeling and constrained inversion can reduce the multiple solutions for single geophysical data and thus improve the recognition of deep formations. The MT-constrained inversion is consistent with the geological features in the seismic section. This suggests that the joint MT and gravity modeling and constrained inversion can be used to delineate deep targets in similar basins.
基金sponsored by the National Major Project(No.2016ZX05014-001)the National Natural Science Foundation of China(No.41172130 and U1403191)the Fundamental Research Funds for the Central Universities(No.2-9-2015-209)
文摘Borehole-to-surface electrical imaging （BSEI） uses a line source and a point source to generate a stable electric field in the ground. In order to study the surface potential of anomalies, three-dimensional forward modeling of point and line sources was conducted by using the finite-difference method and the incomplete Cholesky conjugate gradient （ICCG） method. Then, the damping least square method was used in the 3D inversion of the formation resistivity data. Several geological models were considered in the forward modeling and inversion. The forward modeling results suggest that the potentials generated by the two sources have different surface signatures. The inversion data suggest that the low- resistivity anomaly is outlined better than the high-resistivity anomaly. Moreover, when the point source is under the anomaly, the resistivity anomaly boundaries are better outlined than when using a line source.
基金supported by the National Natural Science Foundation of China (No. 10872144)the Global Environmental Foundation (No. TF053183)
文摘As a primary parameter in the water quality model for shallow bays, the dispersion coefficient is traditionally determined with a trial-and-error method, which is time-consuming and requires much experience. In this paper, based on the measured data of chemical oxygen demand （COD）, the dispersion coefficient is calculated using an inversion method. In the process, the regularization method is applied to treat the ill-posedness, and an operator identity perturbation method is used to obtain the solu- tion. Using the model with an inverted dispersion coefficient, the distributions of COD, inorganic nitrogen （IN）, and inorganic phosphorus （IP） in Bohai Bay are predicted and compared with the measured data. The results indicate that the method is feasible and the inverted dispersion coefficient can be used to predict other pollutant distribution. This method may also be further extended to the inversion of other parameters in the water quality model.
基金supported financially by the Ministry of Science and Technology of China(2011CB808904)the National Natural Science Foundation of China(Nos.41021063,41174075,41004034,41174043,and 41274090)
文摘Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks,which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block,and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization,which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models.The damped least squares method is employed in seismic traveltime inversion,which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model.
基金jointly supported by the National Key R&D Program of China(Grant No.2016YFC0600201)China Geological Survey project(Grant Nos.DD20190012,DD20160082)the National Natural Science Foundation of China(Grant Nos.92062108,41630320,41574133)。
文摘The Zhuxi tungsten deposit in Jiangxi Province,South China,contains a total W reserve of about 2.86 Mt at an average grade of 0.54 wt%WO3,representing the largest W deposit in the world.Numerous studies on the metallogeny of the deposit have included its timing,the ore-controlling structures and sedimentary host rocks and their implications for mineral exploration.However,the deep nappe structural style of Taqian-Fuchun metallogenic belt that hosts the W deposit,and the spatial shape and scale of deeply concealed intrusions and their sedimentary host rocks are still poorly defined,which seriously restricts the discovery of new deposits at depth and in surrounding areas of the W deposit.Modern 3 D geological modeling is an important tool for the exploration of concealed orebodies,especially in brownfield environments.There are obvious density contrast and weak magnetic contrast in the ore-controlling strata and granite at the periphery of the deposit,which lays a physical foundation for solving the 3 D spatial problems of the ore-controlling geological body in the deep part of the study area through gravity and magnetic modeling.Gravity data(1:50000)and aeromagnetic data(1:50000)from the latest geophysical surveys of 2016-2018 have been used,firstly,to carry out a potential field separation to obtain residual anomalies for gravity and magnetic interactive inversion.Then,on the basis of the analysis of the relationship between physical properties and lithology,under the constraints of surface geology and borehole data,human-computer interactive gravity and magnetic inversion for 18 cross-sections were completed.Finally,the 3 D geological model of the Zhuxi tungsten deposit and its periphery have been established through these 18 sections,and the spatial shape of the intrusions and strata with a depth of 5 km underground were obtained,initially realizing―transparency‖for ore-controlling bodies.According the analysis of the geophysical,geochemical,and geological characteristics of the Zhuxi tungsten deposit,we discern three principles for prospecting and prediction in the research area,and propose five new exploration targets in its periphery.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 61172031 and 41175012)
文摘The maritime tropospheric duct is a low-altitude anomalous refractivity structure over the ocean surface,and it can significantly affect the performance of many shore-based/shipboard radar and communication systems.We propose the idea that maritime tropospheric ducts can be retrieved from ocean forward-scattered low-elevation global positioning system（GPS） signals.Retrieval is accomplished by matching the measured power patterns of the signals to those predicted by the forward propagation model as a function of the modified refractivity profile.On the basis of a parabolic equation method and bistatic radar equation,we develop such a forward model for computing the trapped propagation characteristics of an ocean forward-scattered GPS signal within a tropospheric duct.A new GPS scattering initial field is defined for this model to start the propagation modeling.A preliminary test on the performance of this model is conducted using measured data obtained from a 2009-experiment in the South China Sea.Results demonstrate that this model can predict GPS propagation characteristics within maritime tropospheric ducts and serve as a forward model for duct inversion.
基金funded by R&D Department of China National Petroleum Corporation(2022DQ0604-04)the Strategic Cooperation Technology Projects of CNPC and CUPB(ZLZX2020-03)the Science Research and Technology Development of PetroChina(2021DJ1206).
文摘Seismic impedance inversion is an important technique for structure identification and reservoir prediction.Model-based and data-driven impedance inversion are the commonly used inversion methods.In practice,the geophysical inversion problem is essentially an ill-posedness problem,which means that there are many solutions corresponding to the same seismic data.Therefore,regularization schemes,which can provide stable and unique inversion results to some extent,have been introduced into the objective function as constrain terms.Among them,given a low-frequency initial impedance model is the most commonly used regularization method,which can provide a smooth and stable solution.However,this model-based inversion method relies heavily on the initial model and the inversion result is band limited to the effective frequency bandwidth of seismic data,which cannot effectively improve the seismic vertical resolution and is difficult to be applied to complex structural regions.Therefore,we propose a data-driven approach for high-resolution impedance inversion based on the bidirectional long short-term memory recurrent neural network,which regards seismic data as time-series rather than image-like patches.Compared with the model-based inversion method,the data-driven approach provides higher resolution inversion results,which demonstrates the effectiveness of the data-driven method for recovering the high-frequency components.However,judging from the inversion results for characterization the spatial distribution of thin-layer sands,the accuracy of high-frequency components is difficult to guarantee.Therefore,we add the model constraint to the objective function to overcome the shortages of relying only on the data-driven schemes.First,constructing the supervisor1 based on the bidirectional long short-term memory recurrent neural network,which provides the predicted impedance with higher resolution.Then,convolution constraint as supervisor2 is introduced into the objective function to guarantee the reliability and accuracy of the inversion results,which makes the synthetic seismic data obtained from the inversion result consistent with the input data.Finally,we test the proposed scheme based on the synthetic and field seismic data.Compared to model-based and purely data-driven impedance inversion methods,the proposed approach provides more accurate and reliable inversion results while with higher vertical resolution and better spatial continuity.The inversion results accurately characterize the spatial distribution relationship of thin sands.The model tests demonstrate that the model-constrained and data-driven impedance inversion scheme can effectively improve the thin-layer structure characterization based on the seismic data.Moreover,tests on the oil field data indicate the practicality and adaptability of the proposed method.
文摘Flow units(FU)rock typing is a common technique for characterizing reservoir flow behavior,producing reliable porosity and permeability estimation even in complex geological settings.However,the lateral extrapolation of FU away from the well into the whole reservoir grid is commonly a difficult task and using the seismic data as constraints is rarely a subject of study.This paper proposes a workflow to generate numerous possible 3D volumes of flow units,porosity and permeability below the seismic resolution limit,respecting the available seismic data at larger scales.The methodology is used in the Mero Field,a Brazilian presalt carbonate reservoir located in the Santos Basin,who presents a complex and heterogenic geological setting with different sedimentological processes and diagenetic history.We generated metric flow units using the conventional core analysis and transposed to the well log data.Then,given a Markov chain Monte Carlo algorithm,the seismic data and the well log statistics,we simulated acoustic impedance,decametric flow units(DFU),metric flow units(MFU),porosity and permeability volumes in the metric scale.The aim is to estimate a minimum amount of MFU able to calculate realistic scenarios porosity and permeability scenarios,without losing the seismic lateral control.In other words,every porosity and permeability volume simulated produces a synthetic seismic that match the real seismic of the area,even in the metric scale.The achieved 3D results represent a high-resolution fluid flow reservoir modelling considering the lateral control of the seismic during the process and can be directly incorporated in the dynamic characterization workflow.
基金Supported by the National Natural Science Foundation of China(41904109,41974146)National Science and Technology Major Project(2017ZX05019-005)+2 种基金China Postdoctoral Science Foundation(2018M640663)the Shandong Province Postdoctoral Innovation Projects(sdbh20180025)National Key Laboratory of Electromagnetic Environment Projects(6142403200307)。
文摘Based on the pseudo-analytical equation of electromagnetic log for layered formation,an optimal boundary match method is proposed to adaptively truncate the encountered formation structures.An efficient integral method is put forward to significantly accelerate the convergence of Sommerfeld integral.By asymptotically approximating and subtracting the first reflection/transmission waves from the scattered field,the new Sommerfeld integral method has addressed difficulties encountered by the traditional digital filtering method,such as low computational precision and limited operating range,and realized the acceleration of the computation speed of logging-while-drilling electromagnetic measurements(LWD EM).By making use of the priori information from the offset/pilot wells and interactively adjusting the formation model,the optimum initial guesses of the inversion model is determined in order to predict the nearby formation boundaries.The gradient optimization algorithm is developed and an interactive inversion system for the LWD EM data from the horizontal wells is established.The inverted results of field data demonstrated that the real-time interactive inversion method is capable of providing the accurate boundaries of layers around the wellbore from the LWD EM,and it will benefit the wellbore trajectory optimization and reservoir interpretation.
基金Under the auspices of the Key Program of National Natural Science Foundation of China(No.50339010)Huaihe Valley Open Fund Projects(No.Hx2007)
文摘Research on quantitative models of suspended sediment concentration （SSC） using remote sensing technology is very important to understand the scouting and siltation variation in harbors and water channels. Based onlaboratory study of the relationship between different suspended sediment concentrations and reflectance spectra measured synchronously, quantitative inversion models of SSC based on single factor, band ratio and sediment parameter were developed, which provides an effective method to retrieve the SSC from satellite images. Results show that the bl （430-500nm） and b3 （670-735nm） are the optimal wavelengths for the estimation of lower SSC and the b4 （780-835nm） is the optimal wavelength to estimate the higher SSC. Furthermore the band ratio B2/B3 can be used to simulate the variation of lower SSC better and the B4/B1 to estimate the higher SSC accurately. Also the inversion models developed by sediment parameters of higher and lower SSCs can get a relatively higher accuracy than the single factor and band ratio models.
基金Supported by Project of National Natural Science Foundation(No.41430322)
文摘There is not much correlation between sand content and P-wave impedance in faulted basins,and the prediction results cannot be guaranteed. Due to sedimentary facies control,there is inconsistency between sand content and low frequency trend of P-wave impedance,causing problems for seismic inversion modeling,which directly affects final seismic inversion results. The acoustic impedance increases with burial depth. When the same layers of sand and shale formation endure different compaction,the acoustic impedance values will be different. Therefore,The seismic inversion modeling of faulted basin is different from that of conventional basins.The authors built an inversion model using uncompact sandstone percentage,directly compensating for the low frequency trend of the inversion model. In addition,the model's intermediate frequency is similar to P-wave impedance,ensuring that the inversion results are converged correctly. The final results of the inversion can be used directly as sandstone percentage. The aforementioned method was applied to Liaohe Beach in the Bijialing region and obtained optimistic lithology inversion effects: the results were correlated well with the depositional facies map and with lithology in the well borehole. The inversion results can be used to define the sand body horizontally and can separate sand bodies vertically,which is very difficult on conventional seismic section.Therefore the inversion results played an important role in reservoir prediction.
