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.

Four-channel catalytic micro-reactor based on alumina hollow fiber membrane for efficient catalytic oxidation of CO

The traditional automotive catalytic converter using commercial ceramic honeycomb carriers has many problems such as high back pressure,low engine efficiency,and high usage of precious metals.This study proposes a four-channel catalytic micro-reactor based on alumina hollow fiber membrane,which uses phase inversion method for structural molding and regulation.Due to the advantages of its carrier,it can achieve lower ignition temperature under low noble metal loading.With Pd/CeO_(2) at a loading rate of 2.3%(mass),the result showed that the reaction ignition temperature is even less than 160℃,which is more than 90℃ lower than the data of commercial ceramic substrates under similar catalyst loading and airspeed conditions.The technology in turn significantly reduces the energy consumption of the reaction.And stability tests were conducted under constant conditions for 1000 h,which proved that this catalytic converter has high catalytic efficiency and stability,providing prospects for the design of innovative catalytic converters in the future.

Inverse reliability analysis and design for tunnel face stability considering soil spatial variability

The traditional deterministic analysis for tunnel face stability neglects the uncertainties of geotechnical parameters,while the simplified reliability analysis which models the potential uncertainties by means of random variables usually fails to account for soil spatial variability.To overcome these limitations,this study proposes an efficient framework for conducting reliability analysis and reliability-based design(RBD)of tunnel face stability in spatially variable soil strata.The three-dimensional(3D)rotational failure mechanism of the tunnel face is extended to account for the soil spatial variability,and a probabilistic framework is established by coupling the extended mechanism with the improved Hasofer-Lind-Rackwits-Fiessler recursive algorithm(iHLRF)as well as its inverse analysis formulation.The proposed framework allows for rapid and precise reliability analysis and RBD of tunnel face stability.To demonstrate the feasibility and efficacy of the proposed framework,an illustrative case of tunnelling in frictional soils is presented,where the soil's cohesion and friction angle are modelled as two anisotropic cross-correlated lognormal random fields.The results show that the proposed method can accurately estimate the failure probability(or reliability index)regarding the tunnel face stability and can efficiently determine the required supporting pressure for a target reliability index with soil spatial variability being taken into account.Furthermore,this study reveals the impact of various factors on the support pressure,including coefficient of variation,cross-correlation between cohesion and friction angle,as well as autocorrelation distance of spatially variable soil strata.The results also demonstrate the feasibility of using the forward and/or inverse first-order reliability method(FORM)in high-dimensional stochastic problems.It is hoped that this study may provide a practical and reliable framework for determining the stability of tunnels in complex soil strata.

Changes of coastline and tidal flat and its implication for ecological protection under human activities: Take China’s Bohai Bay as an example

The change processes and trends of shoreline and tidal flat forced by human activities are essential issues for the sustainability of coastal area,which is also of great significance for understanding coastal ecological environment changes and even global changes.Based on field measurements,combined with Linear Regression(LR)model and Inverse Distance Weighing(IDW)method,this paper presents detailed analysis on the change history and trend of the shoreline and tidal flat in Bohai Bay.The shoreline faces a high erosion chance under the action of natural factors,while the tidal flat faces a different erosion and deposition patterns in Bohai Bay due to the impact of human activities.The implication of change rule for ecological protection and recovery is also discussed.Measures should be taken to protect the coastal ecological environment.The models used in this paper show a high correlation coefficient between observed and modeling data,which means that this method can be used to predict the changing trend of shoreline and tidal flat.The research results of present study can provide scientific supports for future coastal protection and management.

A Two-step Inverse Method for Quantitative Damage Evaluation of Plate-like Structures Using Vibration Approach

This study presents a novel two-step approach to assess plate-like structural laminar damages,particularly for delamination damage detection of composite structures.Firstly,a 2-D continuous wavelet transform is employed to identify the damage location and sizes from vibration curvature data.An inverse method is subsequently then used to determine the bending stiffness reduction ratio along a specified direction,enabling the quantification of the delamination severity.The method employed in this study is an extension of the one-dimensional inverse method developed in a previous work of the authors.The applicability of the two-step inverse approach is demonstrated in a simulation analysis and by an experimental study on a cantilever composite plate containing a single delamination.The inverse method is shown to have the capacity to reveal the detailed damage information of delamination within a constrained searching space and can be used to determine the effective flexural stiffness of composite plate structures,even in cases of complex delamination damage.

Experimental Procedure and Hardening Model to Consider Forming and Baking Effects in Simulation

This paper investigated an experimental method for bake hardening properties, a technique for deriving the true stress-strain curves after reaching the maximum load, and a constitutive equation considering both work hardening and bake hardening in order to apply the work hardening occurring in the forming process of parts and the bake hardening induced in the baking process to an automotive crash simulation. A general bake hardening test is that a pre-tensioned specimen is baked and then the same specimen is tensioned again without any further treatment. For a bake hardening test of automotive steel with a tensile strength of 1.2 GPa or more, fractures often occur in curvature section outside, an extensometer due to the difference in the material strength caused by non-uniform bake hardening. This causes a problem in that the bake hardening properties cannot be obtained. In this paper, to prevent curvature fracture, tensile specimens were re-machined in the uniformly deformed region of large specimens subjected to pre-strain, and the re-machined specimens with uniform strength in all regions were re-tensioned. In the bake hardening test of ultra-high strength steels with a tensile strength of 1 GPa or more, shear band fractures occur when the pre-strain is large. This makes it impossible to obtain a true stress-strain curve because there is no uniformly deformed region under a tensile test. To overcome this problem, a new method to calculate the true stress-strain curve by comparing experimental results and the load calculated by the local strain obtained from digital images was developed. This method can be applied not only where shear band deformation occurs, but also in necking deformation, and true stress-strain curves for strains up to 2 - 3 times the uniform elongation can be obtained. A new constitutive equation was developed since an appropriate hardening model is required to simultaneously apply the work hardening and the bake hardening to the simulation. For the newly developed model, the user material subroutine of LS-Dyna was configured, and the simulation was performed on the single hat specimens with pre-strain. When both work hardening and bake hardening were considered, there was a significant increase in absorbed energy compared to when only work hardening was considered. This means that both work hardening and bake hardening should be considered in the car crash simulations to enhance the accuracy of the simulation.

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.