This technical note deals with the cauchy and related empirical transparent dispersion formulae to calculate the real n and imaginary k parts of the complex refractive index for a material. Highresolution threeterm avo inversion via a trivariate. This promoted the univariate cauchy probability distribution for high resolution radon transform and avo inversion downton and lines, 2004. Stat 451 lecture notes 0512 simulating random variables.
Therefore, using the cauchy distribution as a prior results in a regularization term that produces sparse solutions. Simulation lecture 2 inversion and transformation methods. The characteristic function of a random variable with the distribution n. Ir, denote any cumulative distribution function cdf continuous or. In this chapter, we present basic methods of generating random variables and.
Here we address some of the limitations of mvfsa by developing a new stochastic method, named greedy annealed importance sampling gais. We can generate from the cauchy distribution using inverse transform sampling, which says that if you plug random uniforms into the inverse cdf of a distribution, then what you get out has that distribution. A study on the method of joint inversion of multiscale. For a discrete distribution, the pdf is the probability that the variate takes the value x.
Invert cdf, return x i f 1 x u i i transformation method. Nonlinear multiparameter inversion using a hybrid global. Generating random variables consortium for mathematics and its. A bayesian impedance inversion method is proposed, under the assumption that the prior distribution is a gamma distribution. On the generalization of seismic tomography algorithms. Application of new monte carlo method for inversion of prestack seismic data yang xue1, mrinal k. The standard cauchy distribution has probability density function fx 1. In practice, computing the complex integral can be done by using the cauchy residue theorem. Conclusionsnonlinear inversion is the direction and trend of mt inversion. U has a cauchy distribution, eu, and all higher order. In this paper, we address these problems by using exact zoeppritz equations. Therefore a 2d cauchy probability density function pdf in the killing angles domain is defined 3,4.
The inversion test based on the model and the actual data prove that the inversion method described in this paper can improve the multisolution problem of inversion by using the information of. For the inversion, the computation of the covariance matrix and regularized terms requires prior estimation of model parameters, which makes the iterative inversion weakly nonlinear. Alemie and sacchi 2011, to enhance the inversion resolution. The scale parameter in the cauchy distribution is not very easy to estimate, and different methods are tested. Due to the positivity constraints 3 the killing angles lie in the domain. Whe cauchy distribution has heavier tails than the gaussian distribution.
Comparing the cauchy and gaussian normal density functions. We then combined the idea of generalized linear inversion with iterative reweighed leastsquares irls algorithm to solve the problem. With the sparsespike assumption of the reflectivity, taking logarithm transform to the poster pdf, the object function of the joint inversion of surface seismic data and borehole seismic data. The following is the plot of the cauchy inverse survival function. A blockyness constraint for seismic ava inversion core. Mar 07, 2016 this method is sometimes called the inversion method. The probability density function pdf, and the cumulative distribution function cdf, f, for the uniforma,b random. Theoretical model equation of cauchy transparent the earliest dispersion formula was established by cauchy 1836 who set up simple empirical dispersion.
The boxmuller transform starts with 2 random uniform numbers \u\ and \v\ generate an exponentially distributed variable \r2\ from \u\ using the inverse transform method this means that \r\ is an exponentially distributed variable on \0, \infty\ generate a variable \\theta\ uniformly. Prestack ava inversion of exact zoeppritz equations based. Characteristic functions and the central limit theorem. The inverse cdf method for simulating from a distribution. Normal random numbers can also be generated using the general inverse transform method e. First, the inverse problem was constructed and the modified cauchy distribution was introduced as the prior information by utilizing bayes theorem. Compared with commonly used gauss distribution and cauchy distribution, modified cauchy prior constraint fits the real data well, and can effectively protect weak, a study of multiscale seismic data joint inversion method. The cauchy distribution is used to model random variables with very heavy tails. In this method, cauchy distribution is utilized to the ratios of fluid factors, shear moduli and densities of perturbation media and background homogeneous media, respectively. The inversion method can also be applied to discrete.
The l 1norm constraint and the cauchy constraint can be applied to the leastsquares inversion, in order to generate a seismic reflectivity series with a longertailed distribution than the gaussian distribution, and thus to increase the effective bandwidth of the retrieved reflectivity series levy and fulagar 1981, sacchi et al 1998. Mar 01, 2017 alemie and sacchi 2011 proposed a highresolution threeterm avo inversion method based on the trivariate cauchy probability distribution by comparing it with the multivariate gaussian distribution. There is also a variant of boxmuller that does not require the use of expensive trigonometric. Examples of initialization of one or a batch of distributions. The standard cauchy distribution has pdf and cdf fx. Global optimization methods such as very fast simulated annealing vfsa and a multiple vfsa mvfsa inversion have been applied to seismic waveform inversion and uncertainty characterization respectively. The case where t 0 and s 1 is called the standard cauchy distribution. It is also known, especially among physicists, as the lorentz distribution after hendrik lorentz, cauchy lorentz distribution, lorentzian function, or breitwigner distribution. On the program implementation of a simple markov homogeneous. The standard form of 1 can be obtained by replacing a with 0 and b with 1 and is given by 2 1 0,1.
The equation for the standard cauchy distribution reduces to. Conventionally, the prior distribution in the inversion focuses on the gaussian distribution, huber distribution, cauchy or the improved cauchy distribution. For a general discussion of the cauchy distribution in statistics, the reader is referred to johnson and kotz 1970. The most notable difference between the gaussian and cauchy distributions is the longer and flatter tail of the latter. Stat 3701 notes on basic monte carlo fall 2020 september 11, 2020 contents 1 introduction 1 2 inversion method 2. The first general method that we present is called the inverse transform method. Suppose g has a geometric distribution, so the mass function is p. Seismic wave scattering inversion for fluid factor of. The results of the simulation analysis of the estimator properties have allowed conclusions to be drawn as concerning the application of the considered methods.
Posts inversion formula edit posts inversion formula for laplace transforms, named after emil post, 3 is a simplelooking but usually impractical formula for evaluating an inverse laplace transform. Sep 22, 20 prior information plays an important role in bayesian avo inversion. Cauchy prior distributionbased avo elastic parameter. When x has a continuous distribution then it has a probability density function fx. We could thus use the discrete inversetransform method, but of course it involves computing in advance pieces like k. At the same time, the relations among the model parameters are assumed. Blind inversion in nonlinear spacevariant imaging by using. Cauchy example we consider estimation of the location parameter in the cauchy location family f xx.
Theorem the inverse of a standard cauchy random variable x is. Pdf seismic parameter estimation from avo inversion. Unlike local optimization methods, sa method attempts to find the global minimum of the objective function. P, x2 p, which we wish to simulate i can nd a function. The existence of the moments of the cauchy distribution under. An alternative to the cauchy distribution sciencedirect. Index terms cauchy machine, helmholtz free energy, spacevariant imaging, blind inversion, sensor nonlinearities. Im trying to use the inverse cumulative distribution method to plot a histogram from the standard cauchy distribution and im getting a strange plot that doesnt look like the textbook standard cauchy. The cauchy distribution, named after augustin cauchy, is a continuous probability distribution. Pdf blind inversion in nonlinear spacevariant imaging.
Proof let the random variable x have the standard cauchy distribution. Cauchy priori distribution based bayesian avo reflectivity inversion may lead to sparse estimates that are sensitive to large reflectivities. However, not all the parameters in field data belong to one of the above distributions. Multichannel algorithms for seismic reflectivity inversion. The cauchy distribution is a member of the locationscale family, i. Transformation methods i suppose we i have a random variable y. With the sparsespike assumption of the reflectivity, taking logarithm transform to the.
Cauchy distribution is able to obtain strong reflections with large amplitudes and has superior noise handling performance. Cauchy distribution the cauchy distribution has pdf given by. In this paper, cauchy distribution is utilized for a priori probability distribution of the model parameters sacchi and ulrych 1995. Non linear multiparameter inversion 993 among local sampling of the model space. I know their corresponding probability density functions and distribution functions, and i think i have to calculate with hand on paper the inverse of fx for each one, right. Instead of building nmdimensional cauchy distribution, nmproduct of onedimensional cauchy distribution is applied sen and stoffa 1995. Chapter 2 inversion method, transformation methods and rejection sampling. The fact that the cauchy distribution has no moments, and therefore the law of large numbers does not apply, motivates researchers to generalize the cauchy. The assumptions above can be weakened to some extent. Mathematical and statistical functions for the cauchy distribution, which is commonly used in physics and finance.
We consider two methods to generate cauchy variate samples here. This triggered a series of articles in which cauchy developed a distribution, now called the cauchy distribution, which would force the method of ordinary least squares to fail. When its parameters correspond to a symmetric shape, the sortofmean is found by symmetry, and since the cauchy has no finite variance, that cant be used to match to a gaussian either. The inverse transform method of generating random variables. Stigler, 1974 he had discovered that the method of ordinary least squares gave the best linear unbiased estimator. Comparison of the leastsquarescriterion and the cauchy. I think i have my inverse function correct x tanpix 12 so i. Cauchy inver sion is also frequently used in the geophysical inversion as a robust optimization method 3. Prestack ava inversion of exact zoeppritz equations based on. When using statistical inversion based on bayes rule, the prior distribution is used to shape the solution, and the modified cauchy norm can help provide a solution with better focused layer boundaries. Boxmuller for generating normally distributed random numbers. Adaptiveregularization method with effective model constraint is superior to a simple data fit inversion, which can get a more stable and accurate solution.
Probabilities of intervals can be recovered from the characteristic function using the following inversion theorem. To sample the bernoulli distribution by inversion, take x 1 1 u6p. Cumulative distribution function the cumulative distribution function cdf is the probability that the variable takes a value less than or equal to x. The following algorithm is a direct consequence of the inversetransform method and the fact that cot. In this case, the cauchy distribution of data noise is assumed. Chapter 4 generating nonuniform random variables arizona math. Using inversion, we can sample the cauchy dis tribution by. Blind inversion in nonlinear spacevariant imaging using. Dec 27, 20 with this equation as forward solver, a prestack seismic bayesian inversion method is proposed to estimate the fluid factor of heterogeneous media. Through the test from the field and model data, this method is feasible and this model based on gamma distribution has a result with tall resolution. We can see that the cauchy function has narrower distributions at their peaks and approach zero more slowly away from the mean.
The inverse function method for modeling the cauchy distribution. The most basic distribution is the uniform distribution on 0. Avo inversion with tdistribution as priori constraint seg. To this end, we model the prior distribution of the avo parameters via a trivariate cauchy distribution. Other works relevant to robust inversion methods are beaton and tukey 1974, holland and welsch 1977. I think i have my inverse function correct x tanpix 12 so i would appreciate some help. X is a continuous random variable with distributed function f and density. Invited paper blind inversion in nonlinear spacevariant imaging by using cauchy machine ivica kopriva and harold szu digital media rf laboratory department of electrical and computer engineeringroom 308 george washington university 725 23rd street nw, washington dc 20052, usa email. The cauchy distribution has in nite mean and variance. The cauchy distribution is consistent with the distribution characteristic of the reflectivity series that is calculated from the field well logging curves shown in figure 1.
Use the inverse transform method to generate random variables with the cauchy. Paper open access bayesian acoustic impedance inversion with. Chapter pdf file introduction to probability, statistics, and. The standard cauchy distribution has pdf and cdf f x 1. The probability density function of a cauchy distribution cd with the location parameter, a and scale parameter b is given by,, 0 22 b f xab x b bxa.
This method is sometimes called the inversion method. The rst general method that we present is called the inverse transform method. In this case the usual inverse function need not be defined. Using the transformation technique, y gx 1x is a 11 transformation from.
Highresolution threeterm avo inversion by means of a. Avo inversion with tdistribution as priori constraint. Cauchy and related empirical dispersion formulae for. However, conventional prestack amplitude variation with incidence angle inversion method is strongly influenced by the accuracy of the approximate zoeppritz equations, which suppresses weak reflections coming from the commonly used prior distribution. Oct 17, 2010 this implies that it gives relatively higher probability for large parameters. A study of multiscale seismic data joint inversion method. The inversion approach uses two uniforms to generate one higherprecision uniform via the. Suppose g has a geometric distribution, so the mass function is pg g 1. Impedance joint inversion of borehole and surface seismic. Cauchy distribution the cauchy distribution, or the lorentzian distribution, is a continuous probability distribution that is the ratio of two independent normally distributed random variables if the denominator distribution has mean zero. We found an iterative algorithm to solve the bayesian inversion and, in addition, comparisons are provided with the classical inversion approach that uses a multivariate gaussian prior. The inversion method can also be applied to discrete random variables. Pdf appling simulated annealing algorithm in regularized. The cauchy distribution has the probability density function pdf.
169 988 887 448 129 1180 1445 196 727 1612 1318 1693 1619 736 1529 547 1574 1703 586 1382 657 811 1613 188 738 505 889 1153 800 907 1136 1198 983 1597 1723