共查询到20条相似文献，搜索用时 93 毫秒
1.
Olu Agunloye 《Pure and Applied Geophysics》1981,119(5):1003-1023
A simple numerical inversion scheme for estimatingn-layer model parameters from observed geoelectrical resistivity data can be used in either the space or wavenumber domain. The technique utilizes Madden's Transmission Line Analogy to compute the resistivity transforms and linear filter theory to accomplish the excursions between the space and wavenumber domains. The inversion is effected by an iterative refinement scheme employing the stochastic inverse which is approximate to the generalized inverse. No singular decomposition analysis is required and the scheme is stable under ill conditions. The inversion scheme not only gives the desired estimates; it exposes redundant parameters and irrelevant data and is easily programmed on a desk-top mini computer. Examples of inverse modeling with hypothetical and field data are discussed. 相似文献
2.
F. AMINZADEH 《Geophysical Prospecting》1984,32(5):819-827
In this paper we develop a recursive algorithm to obtain the layer parameters of an elastic medium (density, P-wave velocity, S-wave velocity) from reflection coefficient matrices in terms of energy flux ratios for a non-normal incidence case. We define a layer impedance matrix, analogous to the impedance of an acoustic medium. Next we derive a matrix relationship between the layer impedance matrix of the n+ 1st layer and the reflection coefficient and parameter matrices of the nth layer. This relationship leads to recursively computing the parameters of the subsurface. We show that the elastic case—unlike the acoustic case—allows one to recover the layer parameters from the impedance matrix for non-normal incidence. The results of this work play a key role in the solution of the inverse problem with non-normal-incidence plane-wave seismic data when using a downward continuation technique. 相似文献
3.
Zhaohai Meng 《Geophysical Prospecting》2018,66(3):626-646
The quantitative explanation of the potential field data of three‐dimensional geological structures remains one of the most challenging issues in modern geophysical inversion. Obtaining a stable solution that can simultaneously resolve complicated geological structures is a critical inverse problem in the geophysics field. I have developed a new method for determining a three‐dimensional petrophysical property distribution, which produces a corresponding potential field anomaly. In contrast with the tradition inverse algorithm, my inversion method proposes a new model norm, which incorporates two important weighting functions. One is the L0 quasi norm (enforcing sparse constraints), and the other is depth‐weighting that counteracts the influence of source depth on the resulting potential field data of the solution. Sparseness constraints are imposed by using the L0 quasinorm on model parameters. To solve the representation problem, an L0 quasinorm minimisation model with different smooth approximations is proposed. Hence, the data space (N) method, which is much smaller than model space (M), combined with the gradient‐projected method, and the model space, combined with the modified Newton method for L0 quasinorm sparse constraints, leads to a computationally efficient method by using an N × N system versus an M × M one because N ? M. Tests on synthetic data and real datasets demonstrate the stability and validity of the L0 quasinorm spare norms inversion method. With the aim of obtaining the blocky results, the inversion method with the L0 quasinorm sparse constraints method performs better than the traditional L2 norm (standard Tikhonov regularisation). It can obtain the focus and sparse results easily. Then, the Bouguer anomaly survey data of the salt dome, offshore Louisiana, is considered as a real case study. The real inversion result shows that the inclusion the L0 quasinorm sparse constraints leads to a simpler and better resolved solution, and the density distribution is obtained in this area to reveal its geological structure. These results confirm the validity of the L0 quasinorm sparse constraints method and indicate its application for other potential field data inversions and the exploration of geological structures. 相似文献
4.
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. 相似文献
5.
6.
Some theoretical and practical limits to linear moment tensor inversion of surface waves are analyzed in detail, in particular when one or few stations are used for rapid determination of source parameters. The theory is briefly outlined and steps of preprocessing, especially corrections for source history and propagation, are discussed in order to guarantee the validity of the moment tensor inversion. The inverse problem is first studied from a theoretical point of view. Then the feasibility of the moment tensor inversion is tested from synthetic computations and the formalism is improved in order to obtain a better system conditioning. Finally, some cases of practical nonuniqueness of the solution are shown and possible restrictions to the applicability of the method are discussed. All steps are illustrated with the example of theM
s
=6.9 Erzincan (Turkey) earthquake of March, 13, 1992. 相似文献
7.
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. 相似文献
8.
9.
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. 相似文献
10.
A MAN/COMPUTER INTERPRETATION SYSTEM FOR RESISTIVITY SOUNDINGS OVER A HORIZONTALLY STRAFIFIED EARTH*
H. K. JOHANSEN 《Geophysical Prospecting》1977,25(4):667-691
The proposed system works as follows:
- 1 By a trial-and-error procedure using a graphic display terminal a geologically relevant layer sequence with parameters (ρj, dj) is adjusted to yield roughly the measured curve.
- 2 The resulting layer sequence is used as starting model for an iterative least squares procedure with singular value decomposition. Minimization of the sum of the squares of the logarithmic differences between measured and calculated values with respect to the logarithms of the resistivities and thicknesses as parameters linearizes the problem to a great extent, with two important implications:
- a) a considerable increase in speed (the number of iterations goes down), thus making it cheap to achieve the optimum solution;
- b) the confidence surfaces in parameter space are well approximated by the hyper-ellipsoids defined by the eigenvalues and eigenvectors of the normal equations.
- 3 It is possible to “freeze” any combination of parameters at predetermined values. Thus extra knowledge and/or hypotheses are easily incorporated and can be tested by rerunning step (2). The overall computing time for a practical case is of the order of 10 sec on a CDC 6400.
11.
A method is presented to estimate the elastic parameters and thickness of media that are locally laterally homogeneous using P‐wave and vertically polarized shear‐wave (SV‐wave) data. This method is a ‘layer‐stripping’ technique, and it uses many aspects of common focal point (CFP) technology. For each layer, a focusing operator is computed using a model of the elastic parameters with which a CFP gather can be constructed using the seismic data. Assuming local homogeneity, the resulting differential time shifts (DTSs) represent error in the model due to anisotropy and error in thickness. In the (τ?p) domain, DTSs are traveltimes Δτ that connect error in layer thickness z, vertical slowness q, and ray parameter p. Series expansion is used to linearize Δτ with respect to error in the elastic parameters and thickness, and least‐squares inversion is used to update the model. For stability, joint inversion of P and SV data is employed and, as pure SV data are relatively rare, the use of mode‐converted (PSV) data to represent SV in the joint inversion is proposed. Analytic and synthetic examples are used to demonstrate the utility and practicality of this inversion. 相似文献
12.
Transverse isotropy with a vertical axis of symmetry is a common form of anisotropy in sedimentary basins, and it has a significant influence on the seismic amplitude variation with offset. Although exact solutions and approximations of the PP-wave reflection coefficient for the transversely isotropic media with vertical axis of symmetry have been explicitly studied, it is difficult to apply these equations to amplitude inversion, because more than three parameters need to be estimated, and such an inverse problem is highly ill-posed. In this paper, we propose a seismic amplitude inversion method for the transversely isotropic media with a vertical axis of symmetry based on a modified approximation of the reflection coefficient. This new approximation consists of only three model parameters: attribute A, the impedance (vertical phase velocity multiplied by bulk density); attribute B, shear modulus proportional to an anellipticity parameter (Thomsen's parameter ε−δ); and attribute C, the approximate horizontal P-wave phase velocity, which can be well estimated by using a Bayesian-framework-based inversion method. Using numerical tests we show that the derived approximation has similar accuracy to the existing linear approximation and much higher accuracy than isotropic approximations, especially at large angles of incidence and for strong anisotropy. The new inversion method is validated by using both synthetic data and field seismic data. We show that the inverted attributes are robust for shale-gas reservoir characterization: the shale formation can be discriminated from surrounding formations by using the crossplot of the attributes A and C, and then the gas-bearing shale can be identified through the combination of the attributes A and B. We then propose a rock-physics-based method and a stepwise-inversion-based method to estimate the P-wave anisotropy parameter (Thomsen's parameter ε). The latter is more suitable when subsurface media are strongly heterogeneous. The stepwise inversion produces a stable and accurate Thomsen's parameter ε, which is proved by using both synthetic and field data. 相似文献
13.
Interval inversion of well-logging data for objective determination of textural parameters 总被引：1，自引：0，他引：1
A joint inversion method for the evaluation of well-logging data is presented, which is applicable to determine textural parameters,
i.e., cementation exponent, saturation exponent and tortuosity factor, simultaneously with conventional petrophysical properties.
The inversion techniques used today perform local interpretation. Since the number of unknowns is slightly lower than that
of the data estimated locally to one depth-point, a set of marginally overdetermined inverse problems has to be solved. For
preserving the overdetermination, textural parameters must be kept constant for longer depth intervals (i.e., 200–300 m), despite the fact that they seem to be varying faster with depth according to field experiences. An inversion
method was developed, which inverts data of a greater depth interval jointly in a highly overdetermined inversion procedure
and gives a better resolution (10 m or less) estimate for the textural parameters. In the paper, a set of inversion tests
on synthetic data as well as a field example prove the feasibility of the method. 相似文献
14.
In this paper, we describe a non‐linear constrained inversion technique for 2D interpretation of high resolution magnetic field data along flight lines using a simple dike model. We first estimate the strike direction of a quasi 2D structure based on the eigenvector corresponding to the minimum eigenvalue of the pseudogravity gradient tensor derived from gridded, low‐pass filtered magnetic field anomalies, assuming that the magnetization direction is known. Then the measured magnetic field can be transformed into the strike coordinate system and all magnetic dike parameters – horizontal position, depth to the top, dip angle, width and susceptibility contrast – can be estimated by non‐linear least squares inversion of the high resolution magnetic field data along the flight lines. We use the Levenberg‐Marquardt algorithm together with the trust‐region‐reflective method enabling users to define inequality constraints on model parameters such that the estimated parameters are always in a trust region. Assuming that the maximum of the calculated gzz (vertical gradient of the pseudogravity field) is approximately located above the causative body, data points enclosed by a window, along the profile, centred at the maximum of gzz are used in the inversion scheme for estimating the dike parameters. The size of the window is increased until it exceeds a predefined limit. Then the solution corresponding to the minimum data fit error is chosen as the most reliable one. Using synthetic data we study the effect of random noise and interfering sources on the estimated models and we apply our method to a new aeromagnetic data set from the Särna area, west central Sweden including constraints from laboratory measurements on rock samples from the area. 相似文献
15.
Common‐midpoint moveout of converted waves is generally asymmetric with respect to zero offset and cannot be described by the traveltime series t2(x2) conventionally used for pure modes. Here, we present concise parametric expressions for both common‐midpoint (CMP) and common‐conversion‐point (CCP) gathers of PS‐waves for arbitrary anisotropic, horizontally layered media above a plane dipping reflector. This analytic representation can be used to model 3D (multi‐azimuth) CMP gathers without time‐consuming two‐point ray tracing and to compute attributes of PS moveout such as the slope of the traveltime surface at zero offset and the coordinates of the moveout minimum. In addition to providing an efficient tool for forward modelling, our formalism helps to carry out joint inversion of P and PS data for transverse isotropy with a vertical symmetry axis (VTI media). If the medium above the reflector is laterally homogeneous, P‐wave reflection moveout cannot constrain the depth scale of the model needed for depth migration. Extending our previous results for a single VTI layer, we show that the interval vertical velocities of the P‐ and S‐waves (VP0 and VS0) and the Thomsen parameters ε and δ can be found from surface data alone by combining P‐wave moveout with the traveltimes of the converted PS(PSV)‐wave. If the data are acquired only on the dip line (i.e. in 2D), stable parameter estimation requires including the moveout of P‐ and PS‐waves from both a horizontal and a dipping interface. At the first stage of the velocity‐analysis procedure, we build an initial anisotropic model by applying a layer‐stripping algorithm to CMP moveout of P‐ and PS‐waves. To overcome the distorting influence of conversion‐point dispersal on CMP gathers, the interval VTI parameters are refined by collecting the PS data into CCP gathers and repeating the inversion. For 3D surveys with a sufficiently wide range of source–receiver azimuths, it is possible to estimate all four relevant parameters (VP0, VS0, ε and δ) using reflections from a single mildly dipping interface. In this case, the P‐wave NMO ellipse determined by 3D (azimuthal) velocity analysis is combined with azimuthally dependent traveltimes of the PS‐wave. On the whole, the joint inversion of P and PS data yields a VTI model suitable for depth migration of P‐waves, as well as processing (e.g. transformation to zero offset) of converted waves. 相似文献
16.
Abstract The structure and collapse of linear three-dimensional magnetic neutral points is studied by varying the four parameters (p, q,j|,j⊥ ) that define, in general, the linear field of a neutral point. The effect of these parameters on both the skeleton structure (i.e. the fan and spine) and the actual field line structure of the null is considered. It is found that one current component (j⊥ ) causes the skeleton structure of the null to fold up from its potential state, whereas the other current component (j |;) causes the field lines to bend. The two other parameters (p,q) determine the potential structure of the null and cause the null to transform from a three-dimensional null to a two-dimensional null and from a positive (type B) null to a negative (type A) null. To investigate the collapse of three-dimensional nulls, solutions to the linear, low-β ideal magnetohydrodynamic equations are found. It is found that three-dimensional null points can collapse if the field line foot-points are free and energy can propagate into the system. 相似文献
17.
G. C. Yetirmishli S. E. Kazimova I. E. Kazimov 《Izvestiya Physics of the Solid Earth》2011,47(9):847-856
We present the method for determining the velocity model of the Earth’s crust and the parameters of earthquakes in the Middle
Kura Depression from the data of network telemetry in Azerbaijan. Application of this method allowed us to recalculate the
main parameters of the hypocenters of the earthquake, to compute the corrections to the arrival times of P and S waves at the observation station, and to significantly improve the accuracy in determining the coordinates of the earthquakes.
The model was constructed using the VELEST program, which calculates one-dimensional minimal velocity models from the travel
times of seismic waves. 相似文献
18.
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. 相似文献
19.
A new method for the 2D inversion of induced polarization (IP) data in the time domain has been developed. The entire IP transients were observed and inverted into 2D Cole-Cole earth models, including resistivity, chargeability, relaxation time and the frequency constant. Firstly, a modified 1D time-domain electromagnetic algorithm was used to calculate the response of a layered polarizable ground. The transient signals were then inverted using the Marquardt method to derive the Cole-Cole parameters of each layer. However, model calculations showed that the EM effects could be neglected for the time range (>1 ms) and for the transmitter–receiver distances (<50 m) used in this study. Therefore, the induction effects were not considered for the solution of the 2D inverse problem and a DC solution was applied. An approximative forward algorithm was introduced in order to calculate the IP transients directly in the time domain and in order to speed up the inverse procedure. The approximation is highly accurate, and this is demonstrated by comparing the approximations with their exact solutions up to 3D. The inverse algorithm presented consists of two steps. The transient voltages of an array data set were inverted separately into a two-dimensional resistivity model for each time channel. The time-dependent resistivity of each cell was then interpreted as the response of a homogeneous half-space. In the 2D inversion algorithm, a 3D DC algorithm was used as a forward operator. The method only requires a standard 2D DC inversion and a homogenous half-space Cole-Cole inversion. The developed algorithm has been successfully applied to synthetic data sets and to a field data set obtained from a waste site situated close to Düren in Germany. 相似文献
20.
M. I. Shimelevich E. A. Obornev I. E. Obornev E. A. Rodionov 《Izvestiya Physics of the Solid Earth》2017,53(4):588-597
The iterative approximation neural network method for solving conditionally well-posed nonlinear inverse problems of geophysics is presented. The method is based on the neural network approximation of the inverse operator. The inverse problem is solved in the class of grid (block) models of the medium on a regularized parameterization grid. The construction principle of this grid relies on using the calculated values of the continuity modulus of the inverse operator and its modifications determining the degree of ambiguity of the solutions. The method provides approximate solutions of inverse problems with the maximal degree of detail given the specified degree of ambiguity with the total number of the sought parameters ~n × 103 of the medium. The a priori and a posteriori estimates of the degree of ambiguity of the approximated solutions are calculated. The work of the method is illustrated by the example of the three-dimensional (3D) inversion of the synthesized 2D areal geoelectrical (audio magnetotelluric sounding, AMTS) data corresponding to the schematic model of a kimberlite pipe. 相似文献