共查询到20条相似文献，搜索用时 583 毫秒
1.
Regularization methods are used to recover a unique and stable solution in ill-posed geophysical inverse problems. Due to the connection of homogeneous operators that arise in many geophysical inverse problems to the Fourier basis, for these operators classical regularization methods possess some limitations that one may try to circumvent by wavelet techniques.
In this paper, we introduce a two-step wavelet-based regularization method that combines classical regularization methods with wavelet transform to solve ill-posed linear inverse problems in geophysics. The power of the two-step wavelet-based regularization for linear inversion is twofold. First, regularization parameter choice is straightforward; it is obtained from a priori estimate of data variance. Second, in two-step wavelet-based regularization the basis can simultaneously diagonalize both the operator and the prior information about the model to be recovered. The latter is performed by wavelet-vaguelette decomposition using orthogonal symmetric fractional B-spline wavelets.
In the two-step wavelet-based regularization method, at the first step where fully classical tools are used, data is inverted for the Moore-Penrose solution of the problem, which is subsequently used as a preliminary input model for the second step. Also in this step, a model-independent estimate of data variance is made using nonparametric estimation and L-curve analysis. At the second step, wavelet-based regularization is used to partially recover the smoothness properties of the exact model from the oscillatory preliminary model.
We illustrated the efficiency of the method by applying on a synthetic vertical seismic profiling data. The results indicate that a simple non-linear operation of weighting and thresholding of wavelet coefficients can consistently outperform classical linear inverse methods. 相似文献
In this paper, we introduce a two-step wavelet-based regularization method that combines classical regularization methods with wavelet transform to solve ill-posed linear inverse problems in geophysics. The power of the two-step wavelet-based regularization for linear inversion is twofold. First, regularization parameter choice is straightforward; it is obtained from a priori estimate of data variance. Second, in two-step wavelet-based regularization the basis can simultaneously diagonalize both the operator and the prior information about the model to be recovered. The latter is performed by wavelet-vaguelette decomposition using orthogonal symmetric fractional B-spline wavelets.
In the two-step wavelet-based regularization method, at the first step where fully classical tools are used, data is inverted for the Moore-Penrose solution of the problem, which is subsequently used as a preliminary input model for the second step. Also in this step, a model-independent estimate of data variance is made using nonparametric estimation and L-curve analysis. At the second step, wavelet-based regularization is used to partially recover the smoothness properties of the exact model from the oscillatory preliminary model.
We illustrated the efficiency of the method by applying on a synthetic vertical seismic profiling data. The results indicate that a simple non-linear operation of weighting and thresholding of wavelet coefficients can consistently outperform classical linear inverse methods. 相似文献
2.
The technique of seismic amplitude-versus-angle inversion has been widely used to estimate lithology and fluid properties in seismic exploration. The amplitude-versus-angle inversion problem is intrinsically ill-posed and generally stabilized by the use of L2-norm regularization methods but with drawback of smoothing important boundaries between adjacent layers. In this study, we propose a sparse Bayesian linearized solution for amplitude-versus-angle inversion problem to preserve the sharp geological interfaces. In this regard, a priori constraint term with two regularization functions is presented: the sparse constraint regularization and the low-frequency model information. In addition, to obtain high-resolution reflectivity estimation, the model parameters decorrelation technique combined with dipole decomposition method is employed. We validate the applicability of the presented method by both synthetic and real seismic data from the Gulf of Mexico. The accuracy improvement of the presented method is also confirmed by comparing the results with the commonly used Bayesian linearized amplitude-versus-angle inversion. 相似文献
3.
Iteratively re‐weighted and refined least squares algorithm for robust inversion of geophysical data 
下载免费PDF全文
下载免费PDF全文 A robust metric of data misfit such as the ?1‐norm is required for geophysical parameter estimation when the data are contaminated by erratic noise. Recently, the iteratively re‐weighted and refined least‐squares algorithm was introduced for efficient solution of geophysical inverse problems in the presence of additive Gaussian noise in the data. We extend the algorithm in two practically important directions to make it applicable to data with non‐Gaussian noise and to make its regularisation parameter tuning more efficient and automatic. The regularisation parameter in iteratively reweighted and refined least‐squares algorithm varies with iteration, allowing the efficient solution of constrained problems. A technique is proposed based on the secant method for root finding to concentrate on finding a solution that satisfies the constraint, either fitting to a target misfit (if a bound on the noise is available) or having a target size (if a bound on the solution is available). This technique leads to an automatic update of the regularisation parameter at each and every iteration. We further propose a simple and efficient scheme that tunes the regularisation parameter without requiring target bounds. This is of great importance for the field data inversion where there is no information about the size of the noise and the solution. Numerical examples from non‐stationary seismic deconvolution and velocity‐stack inversion show that the proposed algorithm is efficient, stable, and robust and outperforms the conventional and state‐of‐the‐art methods. 相似文献
4.
The anisotropy of the land surface can be best described by the bidirectional reflectance distribution function (BRDF). As
the field of multiangular remote sensing advances, it is increasingly probable that BRDF models can be inverted to estimate
the important biological or climatological parameters of the earth surface such as leaf area index and albedo. The state-of-the-art
of BRDF is the use of the linear kernel-driven models, mathematically described as the linear combination of the isotropic
kernel, volume scattering kernel and geometric optics kernel. The computational stability is characterized by the algebraic
operator spectrum of the kernel-matrix and the observation errors. Therefore, the retrieval of the model coefficients is of
great importance for computation of the land surface albedos. We first consider the smoothing solution method of the kernel-driven
BRDF models for retrieval of land surface albedos. This is known as an ill-posed inverse problem. The ill-posedness arises
from that the linear kernel driven BRDF model is usually underdetermined if there are too few looks or poor directional ranges,
or the observations are highly dependent. For example, a single angular observation may lead to an under-determined system
whose solution is infinite (the null space of the kernel operator contains nonzero vectors) or no solution (the rank of the
coefficient matrix is not equal to the augmented matrix). Therefore, some smoothing or regularization technique should be
applied to suppress the ill-posedness. So far, least squares error methods with a priori knowledge, QR decomposition method for inversion of the BRDF model and regularization theories for ill-posed inversion were
developed. In this paper, we emphasize on imposing a priori information in different spaces. We first propose a general a priori imposed regularization model problem, and then address two forms of regularization scheme. The first one is a regularized
singular value decomposition method, and then we propose a retrieval method in I
1 space. We show that the proposed method is suitable for solving land surface parameter retrieval problem if the sampling
data are poor. Numerical experiments are also given to show the efficiency of the proposed methods.
Supported by National Natural Science Foundation of China (Grant Nos. 10501051, 10871191), and Key Project of Chinese National
Programs for Fundamental Research and Development (Grant Nos. 2007CB714400, 2005CB422104) 相似文献
5.
地球物理反演是获取地球信息的重要手段,其求解具有严重的不适定性.为获得稳定的反问题结果,通常需要在目标泛函中加入正则化约束项.正确地估计正则化参数一直是地球物理反问题中的难点.目前存在的选取方法需要根据大量的试验来确定正则化参数,工作量十分巨大,并且存在很大的经验性,很难得到最优的正则化参数.针对这个问题,本文提出了一种基于广义Stein无偏风险估计的正则化参数求取方法.该方法的具体思路是通过求解模型参数均方误差的广义Stein无偏风险估计函数,在反问题求解过程中自动求取正则化参数.本文模型测试结果表明,相比于目前常用的方法,通过该方法得到的正则化参数是最优的. 相似文献
6.
A Trade-Off Solution between Model Resolution and Covariance in Surface-Wave Inversion 总被引：1，自引：0，他引：1
Jianghai Xia Yixian Xu Richard D. Miller Chong Zeng 《Pure and Applied Geophysics》2010,167(12):1537-1547
Regularization is necessary for inversion of ill-posed geophysical problems. Appraisal of inverse models is essential for meaningful interpretation of these models. Because uncertainties are associated with regularization parameters, extra conditions are usually required to determine proper parameters for assessing inverse models. Commonly used techniques for assessment of a geophysical inverse model derived (generally iteratively) from a linear system are based on calculating the model resolution and the model covariance matrices. Because the model resolution and the model covariance matrices of the regularized solutions are controlled by the regularization parameter, direct assessment of inverse models using only the covariance matrix may provide incorrect results. To assess an inverted model, we use the concept of a trade-off between model resolution and covariance to find a proper regularization parameter with singular values calculated in the last iteration. We plot the singular values from large to small to form a singular value plot. A proper regularization parameter is normally the first singular value that approaches zero in the plot. With this regularization parameter, we obtain a trade-off solution between model resolution and model covariance in the vicinity of a regularized solution. The unit covariance matrix can then be used to calculate error bars of the inverse model at a resolution level determined by the regularization parameter. We demonstrate this approach with both synthetic and real surface-wave data. 相似文献
7.
A. J. W. DUIJNDAM 《Geophysical Prospecting》1988,36(8):878-898
This paper gives a review of Bayesian parameter estimation. The Bayesian approach is fundamental and applicable to all kinds of inverse problems. Its basic formulation is probabilistic. Information from data is combined with a priori information on model parameters. The result is called the a posteriori probability density function and it is the solution to the inverse problem. In practice an estimate of the parameters is obtained by taking its maximum. Well-known estimation procedures like least-squares inversion or l1 norm inversion result, depending on the type of noise and a priori information given. Due to the a priori information the maximum will be unique and the estimation procedures will be stable except (in theory) for the most pathological problems which are very unlikely to occur in practice. The approach of Tarantola and Valette can be derived within classical probability theory. The Bayesian approach allows a full resolution and uncertainty analysis which is discussed in Part II of the paper. 相似文献
8.
通过引入带有二阶正则算子的正则化项,建立了一种双参数混合正则化方法.为确定最佳正则化参数,这里主要应用L-曲线法、偏差原理和广义交叉校验准则的优化组合来确定.首先对理论模型进行了数值模拟,通过与截断奇异值分解法、共轭梯度法及标准Tikhonov正则化法的结果比较,表明该方法不仅精度高,而且对于数据的随机扰动具有稳定性.... 相似文献
9.
Joint inversion of Rayleigh‐wave dispersion data and vertical electric sounding data: synthetic tests on characteristic sub‐surface models 下载免费PDF全文
下载免费PDF全文 In the traditional inversion of the Rayleigh dispersion curve, layer thickness, which is the second most sensitive parameter of modelling the Rayleigh dispersion curve, is usually assumed as correct and is used as fixed a priori information. Because the knowledge of the layer thickness is typically not precise, the use of such a priori information may result in the traditional Rayleigh dispersion curve inversions getting trapped in some local minima and may show results that are far from the real solution. In this study, we try to avoid this issue by using a joint inversion of the Rayleigh dispersion curve data with vertical electric sounding data, where we use the common‐layer thickness to couple the two methods. The key idea of the proposed joint inversion scheme is to combine methods in one joint Jacobian matrix and to invert for layer S‐wave velocity, resistivity, and layer thickness as an additional parameter, in contrast with a traditional Rayleigh dispersion curve inversion. The proposed joint inversion approach is tested with noise‐free and Gaussian noise data on six characteristic, synthetic sub‐surface models: a model with a typical dispersion; a low‐velocity, half‐space model; a model with particularly stiff and soft layers, respectively; and a model reproduced from the stiff and soft layers for different layer‐resistivity propagation. In the joint inversion process, the non‐linear damped least squares method is used together with the singular value decomposition approach to find a proper damping value for each iteration. The proposed joint inversion scheme tests many damping values, and it chooses the one that best approximates the observed data in the current iteration. The quality of the joint inversion is checked with the relative distance measure. In addition, a sensitivity analysis is performed for the typical dispersive sub‐surface model to illustrate the benefits of the proposed joint scheme. The results of synthetic models revealed that the combination of the Rayleigh dispersion curve and vertical electric sounding methods in a joint scheme allows to provide reliable sub‐surface models even in complex and challenging situations and without using any a priori information. 相似文献
10.
— We discuss and illustrate graphically with simple 2-D problems, four common pitfalls in geophysical nonlinear inversion. The first one establishes that the Lagrange multiplier, used to incorporate a priori information in the geophysical inverse problem, should be the largest value still compatible with a reasonable data fitting. This procedure should be used only when the interpreter is sure about the importance of the a priori information used to stabilize the inverse problem relative to the geophysical observations. Because this is rarely the case, the user should use the smallest Lagrange multiplier still producing stable solutions. The second pitfall is an attempt to automatically estimate the Lagrange multiplier by decreasing it along the iterative process used to solve the nonlinear optimization problem. Consequently, at the last iteration, the Lagrange multiplier may be so small that the problem may become ill-posed and any computed solution in this case is meaningless. The third pitfall is related to the incorporation of a priori information by a technique known as “Jumping.” This formulation, from the viewpoint of the class of Acceptable Gradient Methods, is incomplete and may lead to a premature halt in the iteration, and, consequently, to solutions far from the true one. Finally, the fourth pitfall is an inadequate convergence criterion which stops the iteration when the data misfit drops just below the noise level, irrespective of the fact that the functional to be minimized may not have attained its minimum. This means that the a priori information has not been completely incorporated, so that this stopping criterion partially neutralizes the effect of the stabilizing functional, and opens the possibility of obtaining unstable, meaningless estimates. 相似文献
11.
12.
Different from the stacked seismic data, pre-stack data includes abundant information about shear wave and density. Through inversing the shear wave and density information from the pre-stack data, we can determine oil-bearing properties from different incident angles. The state-of-the-art inversion methods obtain either low vertical resolution or lateral discontinuities. However, the practical reservoir generally has sharp discontinuities between different layers in vertically direction and is horizontally smooth. Towards obtaining the practical model, we present an inversion method based on the regularized amplitude-versus-incidence angle (AVA) data to estimate the piecewise-smooth model from pre-stack seismic data. This method considers subsurface stratum as a combination of two parts: a piecewise smooth part and a constant part. To fix the ill-posedness in the inversion, we adopt four terms to define the AVA inversion misfit function: the data misfit itself, a total variation regularization term acting as a sparsing operator for the piecewise constant part, a Tikhonov regularization term acting as a smoothing operator for the smooth part, and the last term to smoothly incorporate a priori information for constraining the magnitude of the estimated model. The proposed method not only can incorporate structure information and a priori model constraint, but also is able to derive into a convex objective function that can be easily minimized using iterative approach. Compared with inversion results of TV and Tikhonov regularization methods, the inverted P-wave velocity, S-wave velocity and density of the proposed method can better delineate the piecewise-smooth characteristic of strata. 相似文献
13.
2D data modelling by electrical resistivity tomography for complex subsurface geology 总被引：3，自引：1，他引：3
A new tool for two‐dimensional apparent‐resistivity data modelling and inversion is presented. The study is developed according to the idea that the best way to deal with ill‐posedness of geoelectrical inverse problems lies in constructing algorithms which allow a flexible control of the physical and mathematical elements involved in the resolution. The forward problem is solved through a finite‐difference algorithm, whose main features are a versatile user‐defined discretization of the domain and a new approach to the solution of the inverse Fourier transform. The inversion procedure is based on an iterative smoothness‐constrained least‐squares algorithm. As mentioned, the code is constructed to ensure flexibility in resolution. This is first achieved by starting the inversion from an arbitrarily defined model. In our approach, a Jacobian matrix is calculated at each iteration, using a generalization of Cohn's network sensitivity theorem. Another versatile feature is the issue of introducing a priori information about the solution. Regions of the domain can be constrained to vary between two limits (the lower and upper bounds) by using inequality constraints. A second possibility is to include the starting model in the objective function used to determine an improved estimate of the unknown parameters and to constrain the solution to the above model. Furthermore, the possibility either of defining a discretization of the domain that exactly fits the underground structures or of refining the mesh of the grid certainly leads to more accurate solutions. Control on the mathematical elements in the inversion algorithm is also allowed. The smoothness matrix can be modified in order to penalize roughness in any one direction. An empirical way of assigning the regularization parameter (damping) is defined, but the user can also decide to assign it manually at each iteration. An appropriate tool was constructed with the purpose of handling the inversion results, for example to correct reconstructed models and to check the effects of such changes on the calculated apparent resistivity. Tests on synthetic and real data, in particular in handling indeterminate cases, show that the flexible approach is a good way to build a detailed picture of the prospected area. 相似文献
14.
Man's engineering activities are concentrated on the uppermost part of the earth's crust which is called engineering-geologic zone. This zone is characterized by a significant spatialtemporal variation of the physical properties status of rocks, and saturating waters. This variation determines the specificity of geophysical and, particularly, geoelectrical investigations. Planning of geoelectric investigations in the engineering-geologic zone and their subsequent interpretation requires a priori) geologic-geophysical information on the main peculiarities of the engineering-geologic and hydrogeologic conditions in the region under investigation. This information serves as a basis for the creation of an initial geoelectric model of the section. Following field investigations the model is used in interpretation. Formalization of this a priori) model can be achieved by the solution of direct geoelectric problems. An additional geologic-geophysical information realized in the model of the medium allows to diminish the effect of the “principle of equivalence” by introducing flexible limitations in the section's parameters. Further geophysical observations as well as the correlations between geophysical and engineering-geologic parameters of the section permit the following step in the specification of the geolectric model and its approximation to the real medium. Next correction of this model is made upon accumulation of additional information. The solution of inverse problems with the utilization of computer programs permits specification of the model in the general iterational cycle of interpretation. 相似文献
15.
In this study, we formulate an improved finite element model‐updating method to address the numerical difficulties associated with ill conditioning and rank deficiency. These complications are frequently encountered model‐updating problems, and occur when the identification of a larger number of physical parameters is attempted than that warranted by the information content of the experimental data. Based on the standard bounded variables least‐squares (BVLS) method, which incorporates the usual upper/lower‐bound constraints, the proposed method (henceforth referred to as BVLSrc) is equipped with novel sensitivity‐based relative constraints. The relative constraints are automatically constructed using the correlation coefficients between the sensitivity vectors of updating parameters. The veracity and effectiveness of BVLSrc is investigated through the simulated, yet realistic, forced‐vibration testing of a simple framed structure using its frequency response function as input data. By comparing the results of BVLSrc with those obtained via (the competing) pure BVLS and regularization methods, we show that BVLSrc and regularization methods yield approximate solutions with similar and sufficiently high accuracy, while pure BVLS method yields physically inadmissible solutions. We further demonstrate that BVLSrc is computationally more efficient, because, unlike regularization methods, it does not require the laborious a priori calculations to determine an optimal penalty parameter, and its results are far less sensitive to the initial estimates of the updating parameters. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
16.
17.
地球物理反问题线性化处理之后, 各种反演算法归结为对病态线性方程组的求解. 为了快速准确地计算出地球物理参数, 本文提出了一种全新的基于LSQR算法的混合差分进化算法(Hybrid Differential Evolution Algorithm, HDE). 该算法利用LSQR算法给出DE算法的初始种群, 提高DE算法的计算速度和稳定性. 在不同噪声水平下, 对四种正则化方法Tikhonov、TSVD、LSQR和HDE的反演结果进行详细比较. 理论模型和实际数据反演的结果都表明: 改进的HDE算法应用于地球物理反问题的求解是成功的: 反演结果与原设定模型具有较高的相关性, 在稳定性和准确性上较常规的反演算法都具有一定的优势; 而且不需要给定正则化参数, 具有更强的实用性. 相似文献
18.
目前瑞雷波多阶模式频散曲线反演中仅考虑数据的拟合,缺乏对模型的约束,不能很好地刻画地层间断面的问题,针对此问题,研究了瑞雷波多阶模式频散曲线稀疏正则化反演方法.正演模拟基于广义反射-透射系数法,数值计算上采用一种快速求根方法,与二等分方法相比,能够在很短的时间内达到最优的收敛效果;反演建模时采用L1范数正则化方法对模型... 相似文献
19.
The methods of anomaly transformations considered are based on a system of combined analysis of the geophysical field and a priori) information on the structure of a geological object. The methods involve calculation of a transformative polynomial (describing geophysical noise) which makes it possible to separate the residual field component related to the geological characteristic under study in a correlatively optimal way. The structure of the transformative polynomial is determined by the nature of the geophysical noise that is eliminated by the field transformation. Various correlation methods of anomaly transformations arise, depending on the structure of the transformative polynomial chosen. By way of example, the correlation method employed for separating the geophysical anomalies is shown to be highly effective in investigating the local geological structure. 相似文献
20.
A. J. W. DUIJNDAM 《Geophysical Prospecting》1988,36(8):899-918
A parameter estimation or inversion procedure is incomplete without an analysis of uncertainties in the results. In the fundamental approach of Bayesian parameter estimation, discussed in Part I of this paper, the a posteriori probability density function (pdf) is the solution to the inverse problem. It is the product of the a priori pdf, containing a priori information on the parameters, and the likelihood function, which represents the information from the data. The maximum of the a posteriori pdf is usually taken as a point estimate of the parameters. The shape of this pdf, however, gives the full picture of uncertainty in the parameters. Uncertainty analysis is strictly a problem of information reduction. This can be achieved in several stages. Standard deviations can be computed as overall uncertainty measures of the parameters, when the shape of the a posteriori pdf is not too far from Gaussian. Covariance and related matrices give more detailed information. An eigenvalue or principle component analysis allows the inspection of essential linear combinations of the parameters. The relative contributions of a priori information and data to the solution can be elegantly studied. Results in this paper are especially worked out for the non-linear Gaussian case. Comparisons with other approaches are given. The procedures are illustrated with a simple two-parameter inverse problem. 相似文献