Electrical impedance tomography (EIT) aims to reconstruct the conductivity distribution using the boundary measured voltage potential. Traditional regularization based method would suffer from error propagation due to...Electrical impedance tomography (EIT) aims to reconstruct the conductivity distribution using the boundary measured voltage potential. Traditional regularization based method would suffer from error propagation due to the iteration process. The statistical inverse problem method uses statistical inference to estimate unknown parameters. In this article, we develop a nonlinear weighted anisotropic total variation (NWATV) prior density function based on the recently proposed NWATV regularization method. We calculate the corresponding posterior density function, i.e., the solution of the EIT inverse problem in the statistical sense, via a modified Markov chain Monte Carlo (MCMC) sampling. We do numerical experiment to validate the proposed approach.展开更多
Seismic amplitude variation with offset(AVO) inversion is an important approach for quantitative prediction of rock elasticity,lithology and fluid properties.With Biot-Gassmann's poroelasticity,an improved statist...Seismic amplitude variation with offset(AVO) inversion is an important approach for quantitative prediction of rock elasticity,lithology and fluid properties.With Biot-Gassmann's poroelasticity,an improved statistical AVO inversion approach is proposed.To distinguish the influence of rock porosity and pore fluid modulus on AVO reflection coefficients,the AVO equation of reflection coefficients parameterized by porosity,rock-matrix moduli,density and fluid modulus is initially derived from Gassmann equation and critical porosity model.From the analysis of the influences of model parameters on the proposed AVO equation,rock porosity has the greatest influences,followed by rock-matrix moduli and density,and fluid modulus has the least influences among these model parameters.Furthermore,a statistical AVO stepwise inversion method is implemented to the simultaneous estimation of rock porosity,rock-matrix modulus,density and fluid modulus.Besides,the Laplace probability model and differential evolution,Markov chain Monte Carlo algorithm is utilized for the stochastic simulation within Bayesian framework.Models and field data examples demonstrate that the simultaneous optimizations of multiple Markov chains can achieve the efficient simulation of the posterior probability density distribution of model parameters,which is helpful for the uncertainty analysis of the inversion and sets a theoretical fundament for reservoir characterization and fluid discrimination.展开更多
A review of ten-year's practice in developing the improved simultaneous physical retrieval method(ISPRM)is given in the hope that some creative ideas can be drawn from it.The improvement upon the SPRM is associate...A review of ten-year's practice in developing the improved simultaneous physical retrieval method(ISPRM)is given in the hope that some creative ideas can be drawn from it.The improvement upon the SPRM is associated with the under-determinedness of this ill-posed inverse problem.In our experiment,the precondition is observed that prior information must be independent of the satellite measurements.The well-posed retrieval theory has told us that the forward process is fundamental for the retrieval,and it is the bridge between the input of satellite radiance and the output of retrievals.In order to obtain a better result from the forward process. the full advantage of every prior information available must be taken.It is necessary to turn the ill- posed inverse problem into the well-posed one.Then by using the Ridge regression or Bayes algorithm to find the optimal combination among the first guess,the theoretical analogue information and the satellite observations,the impact of the under-determinedness of this inverse problem on the numerical solution is minimized.展开更多
文摘Electrical impedance tomography (EIT) aims to reconstruct the conductivity distribution using the boundary measured voltage potential. Traditional regularization based method would suffer from error propagation due to the iteration process. The statistical inverse problem method uses statistical inference to estimate unknown parameters. In this article, we develop a nonlinear weighted anisotropic total variation (NWATV) prior density function based on the recently proposed NWATV regularization method. We calculate the corresponding posterior density function, i.e., the solution of the EIT inverse problem in the statistical sense, via a modified Markov chain Monte Carlo (MCMC) sampling. We do numerical experiment to validate the proposed approach.
基金the sponsorship of National Grand Project for Science and Technology(2016ZX05024004,2017ZX05009001,2017ZX05032003)the Fundamental Research Funds for the Central Universities(20CX06036A)+1 种基金the Postdoctoral Applied Research Project of Qingdao(QDYY20190040)the Science Foundation from SINOPEC Key Laboratory of Geophysics(wtyjy-wx2019-01-04)。
文摘Seismic amplitude variation with offset(AVO) inversion is an important approach for quantitative prediction of rock elasticity,lithology and fluid properties.With Biot-Gassmann's poroelasticity,an improved statistical AVO inversion approach is proposed.To distinguish the influence of rock porosity and pore fluid modulus on AVO reflection coefficients,the AVO equation of reflection coefficients parameterized by porosity,rock-matrix moduli,density and fluid modulus is initially derived from Gassmann equation and critical porosity model.From the analysis of the influences of model parameters on the proposed AVO equation,rock porosity has the greatest influences,followed by rock-matrix moduli and density,and fluid modulus has the least influences among these model parameters.Furthermore,a statistical AVO stepwise inversion method is implemented to the simultaneous estimation of rock porosity,rock-matrix modulus,density and fluid modulus.Besides,the Laplace probability model and differential evolution,Markov chain Monte Carlo algorithm is utilized for the stochastic simulation within Bayesian framework.Models and field data examples demonstrate that the simultaneous optimizations of multiple Markov chains can achieve the efficient simulation of the posterior probability density distribution of model parameters,which is helpful for the uncertainty analysis of the inversion and sets a theoretical fundament for reservoir characterization and fluid discrimination.
基金Supported by NNSF of China under Grant(49794030#)National"973"Program No.4 (G1998040909#).
文摘A review of ten-year's practice in developing the improved simultaneous physical retrieval method(ISPRM)is given in the hope that some creative ideas can be drawn from it.The improvement upon the SPRM is associated with the under-determinedness of this ill-posed inverse problem.In our experiment,the precondition is observed that prior information must be independent of the satellite measurements.The well-posed retrieval theory has told us that the forward process is fundamental for the retrieval,and it is the bridge between the input of satellite radiance and the output of retrievals.In order to obtain a better result from the forward process. the full advantage of every prior information available must be taken.It is necessary to turn the ill- posed inverse problem into the well-posed one.Then by using the Ridge regression or Bayes algorithm to find the optimal combination among the first guess,the theoretical analogue information and the satellite observations,the impact of the under-determinedness of this inverse problem on the numerical solution is minimized.