Comparison of different inversion methods of D-T_(2)two-dimensional nuclear magnetic resonance logging and applicability analysis

D-T_(2)two-dimensional nuclear magnetic resonance(2D NMR)logging technology can distinguish pore fluid types intuitively,and it is widely used in oil and gas exploration.Many 2D NMR inversion methods(e.g.,truncated singular value decomposition(TSVD),Butler-Reds-Dawson(BRD),LM-norm smoothing,and TIST-L1 regularization methods)have been proposed successively,but most are limited to numerical simulations.This study focused on the applicability of different inversion methods for NMR logging data of various acquisition sequences,from which the optimal inversion method was selected based on the comparative analysis.First,the two-dimensional NMR logging principle was studied.Then,these inversion methods were studied in detail,and the precision and computational efficiency of CPMG and diffusion editing(DE)sequences obtained from oil-water and gas-water models were compared,respectively.The inversion results and calculation time of truncated singular value decomposition(TSVD),Butler-Reds-Dawson(BRD),LM-norm smoothing,and TIST-L1 regularization were compared and analyzed through numerical simulations.The inversion method was optimized to process SP mode logging data from the MR Scanner instrument.The results showed that the TIST-regularization and LM-norm smoothing methods were more accurate for the CPMG and DE sequence echo trains of the oil-water and gas-water models.However,the LM-norm smoothing method was less time-consuming,making it more suitable for logging data processing.A case study in well A25 showed that the processing results by the LM-norm smoothing method were consistent with GEOLOG software.This demonstrates that the LM-norm smoothing method is applicable in practical NMR logging processing.

Inversion method of deflection of the vertical based on SWOT wide-swath altimeter data

The deflection of the vertical(DOV)is the key information in the study of ocean gravity field.However,in most areas,the precision of the prime component of DoV is significantly lower than that of the meridian component.To obtain higher accuracy and resolution of ocean gravity information,researchers have proposed a novel altimeter called the wide-swath altimeter.This altimeter allows for the simultaneous acquisition of high-precision and high-resolution two-dimensional measurements of sea surface height(SSH).In this paper,the Surface Water and Ocean Topography(SWOT)mission with a wide-swath altimeter on board is selected for research.One cycle of swoT sea surface height data is simulated to inverse the DOV in the Arabian Sea(45°E—80°E,0°-30°N),and the inversion results are compared with those of conventional altimeter data.The results demonstrate that the difference between the meridian and prime components derived from the inversion of swoT wide-swath data is minimal,significantly outperforming the inversion results of conventional nadir altimeter data.The advantage of swoT wide-swath altimeter lies in its ability to use the multi-directional geoid slope at any sea surface measurement point to invert the components in the meridian and prime directions.To investigate the impact of this advantage on inversion precision,this paper employs a method to calculate the gradient of the geoid in multiple directions to invert DoV components.The improvement effect of calculating the gradient of the geoid in multiple directions on the precision of DoV component is analyzed.It is found that the accuracy of DoV inversion has significantly improved with the increase of geodetic gradient calculation direction.In addition,the effects of various errors and grid spacing in SwoT wide sea surface height data on the precision of Dov inversion are also analyzed.

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

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

A phase inversion based sponge-like polysulfonamide/SiO_2 composite separator for high performance lithium-ion batteries

In this work,a sponge-like polysulfonamide(PSA)/SiO_2 composite membrane is unprecedentedly prepared by the phase inversion method,and successfully demonstrated as a novel separator of lithium-ion batteries(LIBs).Compared to the commercial polypropylene(PP) separator,the sponge-like PSA/SiO_2 composite possesses better physical and electrochemical properties,such as higher porosity,ionic conductivity,thermal stability and flame retarding ability.The LiCoO_2/Li half-cells using the sponge-like composite separator demonstrate superior rate capability and cyclability over those using the commercial PP separator.Moreover,the sponge-like composite separator can ensure the normal operation of LiCoO_2/Li half-cell at an extremely high temperature of 90 °C,while the commercial PP separator cannot.All these encouraging results suggest that this phase inversion based sponge-like PSA/SiO_2 composite separator is really a promising separator for high performance LIBs.

Detection performance and inversion processing of logging-while-drilling extra-deep azimuthal resistivity measurements

We present systematic investigations on the physics,detection performance and inversion of logging-while-drilling extradeep azimuthal resistivity measurements(EDARM).First,the definitions of EDRAM measurements are discussed,followed by the derivation of the attenuation and phase-shift geometrical factors to illustrate the relative contributions of formation units to the observed signals.Then,a new definition of detection depth,which considers the uncertainty of inversion results caused by the data noise,is proposed to quantify the detection capability of ED ARM.Finally,the B ayesian theory associated with Markov chain Monte Carlo sampling is introduced for fast processing of EDARM data.Numerical results show that ED ARM is capable of detecting the azimuth and distance of remote bed boundaries,and the detection capability increases with increasing spacing and resistivity contrast.The EDARM tool can accommodate a large range of formation resistivity and is able to provide the resistivity anisotropy at arbitrary relative dipping angles.In addition,multiple bed boundaries and reservoir images near the borehole are readily obtained by using the Bayesian inversion.

A time domain induced polarization relaxation time spectrum inversion method based on a damping factor and residual correction

Relaxation time spectra （RTS） derived from time domain induced polarization data （TDIP） are helpful to assess oil reservoir pore structures. However, due to the sensitivity to the signal-to-noise ratio （SNR）, the inversion accuracy of the traditional singular value decomposition （SVD） inversion method reduces with a decrease of SNR. In order to enhance the inversion accuracy and improve robustness of the inversion method to the SNR, an improved inversion method, based on damping factor and spectrum component residual correction, is proposed in this study. The numerical inversion results show that the oscillation of the RTS derived from the SVD method increased with a decrease of SNR, which makes it impossible to get accurate inversion components. However, the SNR has little influence on inversion components of the improved method, and the RTS has high inversion accuracy and robustness. Moreover, RTS derived from core sample data is basically in accord with the pore-size distribution curve, and the RTS derived from the actual induced polarization logging data is smooth and continuous, which indicates that the improved method is practicable.

Low-field NMR inversion based on low-rank and sparsity restraint of relaxation spectra

In this paper,we proposed a novel method for low-field nuclear magnetic resonance(NMR)inversion based on low-rank and sparsity restraint(LRSR)of relaxation spectra,with which high quality construction is made possible for one-and two-dimensional low-field and low signal to noise ratio NMR data.In this method,the low-rank and sparsity restraints are introduced into the objective function instead of the smoothing term.The low-rank features in relaxation spectra are extracted to ensure the local characteristics and morphology of spectra.The sparsity and residual term are contributed to the resolution and precision of spectra,with the elimination of the redundant relaxation components.Optimization process of the objective function is designed with alternating direction method of multiples,in which the objective function is decomposed into three subproblems to be independently solved.The optimum solution can be obtained by alternating iteration and updating process.At first,numerical simulations are conducted on synthetic echo data with different signal-to-noise ratios,to optimize the desirable regularization parameters and verify the feasibility and effectiveness of proposed method.Then,NMR experiments on solutions and artificial sandstone samples are conducted and analyzed,which validates the robustness and reliability of the proposed method.The results from simulations and experiments have demonstrated that the suggested method has unique advantages for improving the resolution of relaxation spectra and enhancing the ability of fluid quantitative identification.

An improved data space inversion method to predict reservoir state fields via observed production data

A data-space inversion(DSI)method has been recently proposed and successfully applied to the history matching and production prediction of reservoirs.Based on Bayesian theory,DSI can directly and effectively obtain good posterior flow predictions without inversion of geological parameters of reservoir model.This paper presents an improved DSI method to fast predict reservoir state fields(e.g.saturation and pressure profiles)via observed production data.Firstly,a large number of production curves and state data are generated by reservoir model simulation to expand the data space of original DSI.Then,efficient history matching only on the observed production data is carried out via the original DSI to obtain related parameters which reflects the weight of the real reservoir model relative to prior reservoir models.Finally,those parameters are used to predict the oil saturation and pressure profiles of the real reservoir model by combining large amounts of state data of prior reservoir models.Two examples including conventional heterogeneous and unconventional fractured reservoir are implemented to test the performances of predicting saturation and pressure profiles of this improved DSI method.Besides,this method is also tested in a real field and the obtained results show the high computational efficiency and high accuracy of the practical application of this method.

Transport properties of a binary mixture of CO_2-N_2 from the pair potential energy functions based on a semi-empirical inversion method

The potential energy snrface of a CO2-N2 mixture is determined by using an inversion method, together with a new collision integral correlation [J. Phys. Chem. R@ Data 19 1179 （1990）]. With the new invert potential, the transport properties of CO2-N2 mixture are presented in a temperature range front 273.15 K to 3273.15 K at low density by employing the Chapman-Enskog scheme and the Wang Chang-Uhlenbeck de Boer theory, consisting of a viscosity coefficient, a thermal conductivity coefficient, a binary diffusion coefficient, and a thermal diffusion factor. The accuracy of the predicted results is estimated to be 2% for viscosity, 5% for thermal conductivity, and 10% for binary diffusion coefficient.

Inverse Construction Methods of Heterogeneous NURBS Object Based on Additive Manufacturing

According to the requirement of heterogeneous object modeling in additive manufacturing(AM),the Non-Uniform Rational B-Spline(NURBS)method has been applied to the digital representation of heterogeneous object in this paper.By putting forward the NURBS material data structure and establishing heterogeneous NURBS object model,the accurate mathematical unified representation of analytical and free heterogeneous objects have been realized.With the inverse modeling of heterogeneous NURBS objects,the geometry and material distribution can be better designed to meet the actual needs.Radical Basis Function(RBF)method based on global surface reconstruction and the tensor product surface interpolation method are combined to RBF-NURBS inverse construction method.The geometric and/or material information of regular mesh points is obtained by RBF interpolation of scattered data,and the heterogeneous NURBS surface or object model is obtained by tensor product interpolation.The examples have shown that the heterogeneous objects fitting to scattered data points can be generated effectively by the inverse construction methods in this paper and 3D CAD models for additive manufacturing can be provided.

Seismic fluid identification using a nonlinear elastic impedance inversion method based on a fast Markov chain Monte Carlo method

Elastic impedance inversion with high efficiency and high stability has become one of the main directions of seismic pre-stack inversion. The nonlinear elastic impedance inversion method based on a fast Markov chain Monte Carlo （MCMC） method is proposed in this paper, combining conventional MCMC method based on global optimization with a preconditioned conjugate gradient （PCG） algorithm based on local optimization, so this method does not depend strongly on the initial model. It converges to the global optimum quickly and efficiently on the condition that effi- ciency and stability of inversion are both taken into consid- eration at the same time. The test data verify the feasibility and robustness of the method, and based on this method, we extract the effective pore-fluid bulk modulus, which is applied to reservoir fluid identification and detection, and consequently, a better result has been achieved.

Magnetic interface forward and inversion method based on Padé approximation

The magnetic interface forward and inversion method is realized using the Taylor series expansion to linearize the Fourier transform of the exponential function. With a large expansion step and unbounded neighborhood, the Taylor series is not convergent, and therefore, this paper presents the magnetic interface forward and inversion method based on Pade approximation instead of the Taylor series expansion. Compared with the Taylor series, Pade＇s expansion＇s convergence is more stable and its approximation more accurate. Model tests show the validity of the magnetic forward modeling and inversion of Pade approximation proposed in the paper, and when this inversion method is applied to the measured data of the Matagami area in Canada, a stable and reasonable distribution of underground interface is obtained.

Calculation of Viscosity and Diffusion Coefficients for Two Binary Gaseous Mixtures Using the Semi-empirical Inversion Method

The accurate reduced potential energies for two binary gas mixtures including benzene- methanol and methane-tetrafluoromethane at low density have been obtained by direct inversion of the viscosity collision integral equations. The kinetic theory along with the extended principle of corresponding-states has been used to calculate the viscosity and dif- fusion coefficients over a wide range of temperature and composition. Good agreements between calculated and experimental data are obtained.

Sensitivity and inverse analysis methods for parameter intervals

This paper proposes a sensitivity analysis method for engineering parameters using interval analyses.This method substantially extends the application of interval analysis method.In this scheme,parameter intervals and decision-making target intervals are determined using the interval analysis method.As an example,an inverse analysis method for uncertainty is presented.The intervals of unknown parameters can be obtained by sampling measured data.Even for limited measured data,robust results can also be obtained with the inverse analysis method,which can be intuitively evaluated by the uncertainty expressed in terms of an interval.For complex nonlinear problems,an iteratively optimized inverse analysis model is proposed.In a given set of loose parameter intervals,all the unknown parameter intervals that satisfy the measured information can be obtained by an iteratively optimized inverse analysis model.The influences of measured precisions and the number of parameters on the results of the inverse analysis are evaluated.Finally,the uniqueness of the interval inverse analysis method is discussed.

Evaluating gravity gradient components based on a reweighted inversion method

In gravity gradient inversion,to choose an appropriate component combination is very important,that needs to understand the function of each component of gravity gradient in the inversion.In this paper,based on the previous research on the characteristics of gravity gradient components,we propose a reweighted inversion method to evaluate the influence of single gravity gradient component on the inversion resolution The proposed method only adopts the misfit function of the regularized equation and introduce a depth weighting function to overcome skin effect produced in gravity gradient inversion.A comparison between different inversion results was undertaken to verify the influence of the depth weighting function on the inversion result resolution.To avoid the premise of introducing prior information,we select the depth weighting function based on the sensitivity matrix.The inversion results using the single-prism model and the complex model show that the influence of different components on the resolution of inversion results is different in different directions,however,the inversion results based on two kind of models with adding different levels of random noise are basically consistent with the results of inversion without noises.Finally,the method was applied to real data from the Vinton salt dome,Louisiana,USA.

Assessment of Three Methods for Estimating Abundance of Ungulate Populations

The accuracy of three methods (pellet group count, inverse line transect method and track count) for estimat-ing the density of ungulate population was assessed, based on a wapiti population with the known size in the enclosure of500 hm2 on Tonghe Hunting Farm, Heilongjiang province, China, during 1992-1993. The estimation of population densityusing the Pelld group count method was the closeSt to the actual value, whereas the density of the wapiti population wasunderestimated by the inverse line transect method and overestimated by the track count method. No signilicant differenceexisted the estimated and actual values for the pellet group method and the inverse line transed method, but thesignificant differen was found for the track count method. Therefore, it is suggested that the pellet group method should begiven the priority in the survey of ungulate population abundance in practice.

GENERALIZED INVERSE RATIONAL EXTRAPOLATION METHODS FOR MATRIX SEQUENCES

Assume that a convergent matrix sequence{An}:An→A（n→∞）, An,A∈C3×3.We want to form a new matrix sequence {Hn}, derived from {An}, which has also A aslimit and whose convergence is faster than the of {An}. Three rational extrapolation meth-ods for accelerating the convergence of matrix sequences {An} are presented in this paper.The underlying methods are based on the generalized inverse for matrices which is

Genetic algorithm-finite element method inversion of the factors determining the recent tectonic stress field of part of East Asia area

Genetic algorithm finite element method (GA FEM) is applied to the study of tectonic stress field of part of East Asia area. From the observed stress distribution, 2 D elastic plane stress inversion is made to deduce the boundary forces and investigate controlling factors. It is suggested that the continent continent collision is the dominant factor controlling the Chinese tectonic stress field. The ocean continent convergence along the subduction zone is an important factor. There exists tensile boundary force along the marginal sea.

Application of homotopy parameter inversion method in Miyun Reservoir

The large-scale convergence of homotopy parametric inversion method on the water quality model parameters calculated was used,with application in parametric inversion calculation of total phosphorus of Beijing Miyun Reservoir.Through calculated and compared the error of sedimentation rate by homotopy parametric inversion method and genetic inversion calculation method,the results indicate that homotopy parametric inversion method has good stability,calculating speed,and even if the initial selection away from the objective function,the solution still has a good convergence.

Estimation of Site Effects on Stations in the Capital Circle Region Using the Generalized Inversion Method

Generalized Inversion Method has been used to estimate the spatial variation of site effects,using the digital data of SH-waves recorded by 63 stations in the Capital Circle Region of China from 2001 to 2006.We gained the site effects of all stations participating in the calculation.We found that the site effect of rock was stabile and about 1.0 from 1.0Hz to 10.0Hz,while the site effect of deposit was high in low frequencies,about 3 ～ 7 from 1.0Hz to 8.0Hz,and the site effect was protuberant at about 5.0Hz,then fell as the frequency increased.The result shows the shape and intensity of station site effects are mainly influenced by the lithology below the station,and possibly also by the local geological structure.