Bayesian seismic multi-scale inversion in complex Laplace mixed domains

Seismic inversion performed in the time or frequency domain cannot always recover the long-wavelength background of subsurface parameters due to the lack of low-frequency seismic records. Since the low-frequency response becomes much richer in the Laplace mixed domains, one novel Bayesian impedance inversion approach in the complex Laplace mixed domains is established in this study to solve the model dependency problem. The derivation of a Laplace mixed-domain formula of the Robinson convolution is the first step in our work. With this formula, the Laplace seismic spectrum, the wavelet spectrum and time-domain reflectivity are joined together. Next, to improve inversion stability, the object inversion function accompanied by the initial constraint of the linear increment model is launched under a Bayesian framework. The likelihood function and prior probability distribution can be combined together by Bayesian formula to calculate the posterior probability distribution of subsurface parameters. By achieving the optimal solution corresponding to maximum posterior probability distribution, the low-frequency background of subsurface parameters can be obtained successfully. Then, with the regularization constraint of estimated low frequency in the Laplace mixed domains, multi-scale Bayesian inversion inthe pure frequency domain is exploited to obtain the absolute model parameters. The effectiveness, anti-noise capability and lateral continuity of Laplace mixed-domain inversion are illustrated by synthetic tests. Furthermore,one field case in the east of China is discussed carefully with different input frequency components and different inversion algorithms. This provides adequate proof to illustrate the reliability improvement in low-frequency estimation and resolution enhancement of subsurface parameters, in comparison with conventional Bayesian inversion in the frequency domain.

Data inversion of multi-dimensional magnetic resonance in porous media

Since its inception in the 1970s,multi-dimensional magnetic resonance(MR)has emerged as a powerful tool for non-invasive investigations of structures and molecular interactions.MR spectroscopy beyond one dimension allows the study of the correlation,exchange processes,and separation of overlapping spectral information.The multi-dimensional concept has been re-implemented over the last two decades to explore molecular motion and spin dynamics in porous media.Apart from Fourier transform,methods have been developed for processing the multi-dimensional time-domain data,identifying the fluid components,and estimating pore surface permeability via joint relaxation and diffusion spectra.Through the resolution of spectroscopic signals with spatial encoding gradients,multi-dimensional MR imaging has been widely used to investigate the microscopic environment of living tissues and distinguish diseases.Signals in each voxel are usually expressed as multi-exponential decay,representing microstructures or environments along multiple pore scales.The separation of contributions from different environments is a common ill-posed problem,which can be resolved numerically.Moreover,the inversion methods and experimental parameters determine the resolution of multi-dimensional spectra.This paper reviews the algorithms that have been proposed to process multidimensional MR datasets in different scenarios.Detailed information at the microscopic level,such as tissue components,fluid types and food structures in multi-disciplinary sciences,could be revealed through multi-dimensional MR.

JACOBI POLYNOMIALS USED TO INVERT THE LAPLACE TRANSFORM

Given the Laplace transform F(s) of a function f(t), we develop a new algorithm to find on approximation to f(t) by the use of the dassical Jacobi polynomials. The main contribution of our work is the development of a new and very effective method to determine the coefficients in the finite series ex-pansion that approximation f(t) in terms of Jacobi polynomials. Some numerical examples are illustrated.

SOLVING COUPLED PSEUDO-PARABOLIC EQUATION USING A MODIFIED DOUBLE LAPLACE DECOMPOSITION METHOD

In this paper, the modification of double Laplace decomposition method is pro- posed for the analytical approximation solution of a coupled system of pseudo-parabolic equation with initial conditions. Some examples are given to support our presented method. In addition, we prove the convergence of double Laplace transform decomposition method applied to our problems.

A Regularization Method for Approximating the Inverse Laplace Transform

A new method for approximating the inerse Laplace transform is presented. We first change our Laplace transform equation into a convolution type integral equation, where Tikhonov regularization techniques and the Fourier transformation are easily applied. We finally obtain a regularized approximation to the inverse Laplace transform as finite sum

DYNAMIC STRESS INTENSITY FACTORS AROUND TWO CRACKS NEAR AN INTERFACE OF TWO DISSIMILAR ELASTIC HALF－PLANES UNDER IN－PLANE SHEAR IMPACT LOAD

Transient stresses around two collinear cracks which lie in parallel with theinterface of the two dissimilar half-planes are studied in this article.The surfaces ofthe cracks are sheared suddenly. Application of the Fourier and Laplace transforms technique reduces the problem to that of solving dual integrai equations.To solvethese,the differences of.the crack surface displacements are expanded in a series offunctions which are automatically zero outside of the cracks. The unknown coefficients accompanied in the series are determined by the Schmidt method. The stress intensity .factors are defined in the Laplace transform domain and these are inverted numerically in the physical space .As an example ,the dynamic stress intensity factors around two cracks in a ceramic and steel bonded composite are numerically calculated.

Fractional-order generalized thermoelastic diffusion theory

The present work aims to establish a fractional-order generalized themoelastic diffusion theory for anisotropic and linearly thermoelastic diffusive media. To numerically handle the multi-physics problems expressed by a sequence of incomplete differential equations, particularly by a fractional equation, a generalized variational principle is obtained for the unified theory using a semi-inverse method. In numerical implementation, the dynamic response of a semi-infinite medium with one end subjected to a thermal shock and a chemical potential shock is investigated using the Laplace transform. Numerical results, i.e., non-dimensional temperature, chemical potential, and displacement, are presented graphically. The influence of the fractional order parameter on them is evaluated and discussed.

Calculating the transient behavior of grounding systems using inverse Laplace transform

This paper deals with a unified and novel approach for analyzing the frequency and time domain performance of grounding systems.The proposed procedure is based on solving the full set of Maxwell's equations in the frequency domain,and enables the exact computation of very near fields at the surface of the grounding grid,as well as far fields,by simple and accurate closed-form expressions for solving Sommerfeld integrals.In addition,the soil ionization is easily considered in the proposed method.The frequency domain responses are converted to the time domain by fast inverse Laplace transform.The results are validated and have shown acceptable accuracy.

Analytical solutions of citrate-phosphate coupled model of rice(Oryza sativa L.)roots

The citrate secreted by the rice(Oryza sativa L.)roots will promote the absorption of phosphate,and this process is described by the Kirk model.In our work,the Kirk model is divided into citrate sub-model and phosphate sub-model.In the citrate sub-model,we obtain the analytical solution of citrate with the Laplace transform,inverse Laplace transform and convolution theorem.The citrate solution is substituted into the phosphate sub-model,and the analytical solution of phosphate is obtained by the separation variable method.The existence of the solutions can be proved by the comparison test,the Weierstrass M-test and the Abel discriminating method.

Efficient Numerical Valuation of Continuous Installment Options

In this work we investigate the novel Kryzhnyi method for the numerical inverse Laplace transformation and apply it to the pricing problem of continuous installment options.We compare the results with the one obtained using other classical methods for the inverse Laplace transformation,like the Euler summation method or the Gaver-Stehfest method.

A Numerical Algorithm for Arbitrary Real-Order Hankel Transform

The Hankel transform is widely used to solve various engineering and physics problems,such as the representation of electromagnetic field components in the medium,the representation of dynamic stress intensity factors,vibration of axisymmetric infinite membrane and displacement intensity factors which all involve this type of integration.However,traditional numerical integration algorithms cannot be used due to the high oscillation characteristics of the Bessel function,so it is particularly important to propose a high precision and efficient numerical algorithm for calculating the integral of high oscillation.In this paper,the improved Gaver-Stehfest(G-S)inverse Laplace transform method for arbitrary real-order Bessel function integration is presented by using the asymptotic characteristics of the Bessel function and the accumulation of integration,and the optimized G-S coefficients are given.The effectiveness of the algorithm is verified by numerical examples.Compared with the linear transformation accelerated convergence algorithm,it shows that the G-S inverse Laplace transform method is suitable for arbitrary real order Hankel transform,and the time consumption is relatively stable and short,which provides a reliable calculation method for the study of electromagnetic mechanics,wave propagation,and fracture dynamics.