共查询到20条相似文献，搜索用时 31 毫秒
1.
Improved Monte Carlo inversion of surface wave data 总被引：2，自引：0，他引：2
Inversion of surface wave data suffers from solution non‐uniqueness and is hence strongly biased by the initial model. The Monte Carlo approach can handle this non‐uniqueness by evidencing the local minima but it is inefficient for high dimensionality problems and makes use of subjective criteria, such as misfit thresholds, to interpret the results. If a smart sampling of the model parameter space, which exploits scale properties of the modal curves, is introduced the method becomes more efficient and with respect to traditional global search methods it avoids the subjective use of control parameters that are barely related to the physical problem. The results are interpreted drawing inference by means of a statistical test that selects an ensemble of feasible shear wave velocity models according to data quality and model parameterization. Tests on synthetic data demonstrate that the application of scale properties concentrates the sampling of model parameter space in high probability density zones and makes it poorly sensitive to the initial boundary of the model parameters. Tests on synthetic and field data, where boreholes are available, prove that the statistical test selects final results that are consistent with the true model and which are sensitive to data quality. The implemented strategies make the Monte Carlo inversion efficient for practical applications and able to effectively retrieve subsoil models even in complex and challenging situations such as velocity inversions. 相似文献
2.
A. Ghorbani C. Camerlynck N. Florsch P. Cosenza A. Revil 《Geophysical Prospecting》2007,55(4):589-605
The inversion of induced‐polarization parameters is important in the characterization of the frequency electrical response of porous rocks. A Bayesian approach is developed to invert these parameters assuming the electrical response is described by a Cole–Cole model in the time or frequency domain. We show that the Bayesian approach provides a better analysis of the uncertainty associated with the parameters of the Cole–Cole model compared with more conventional methods based on the minimization of a cost function using the least‐squares criterion. This is due to the strong non‐linearity of the inverse problem and non‐uniqueness of the solution in the time domain. The Bayesian approach consists of propagating the information provided by the measurements through the model and combining this information with a priori knowledge of the data. Our analysis demonstrates that the uncertainty in estimating the Cole–Cole model parameters from induced‐polarization data is much higher for measurements performed in the time domain than in the frequency domain. Our conclusion is that it is very difficult, if not impossible, to retrieve the correct value of the Cole–Cole parameters from time‐domain induced‐polarization data using standard least‐squares methods. In contrast, the Cole–Cole parameters can be more correctly inverted in the frequency domain. These results are also valid for other models describing the induced‐polarization spectral response, such as the Cole–Davidson or power law models. 相似文献
3.
Prestack seismic inversion plays an important role in estimating elastic parameters that are sensitive to reservoirs and fluid underground. In this paper, a simultaneous inversion method named FMR‐AVA (Fluid Factor, Mu (Shear modulus), Rho (Density)‐Amplitude Variation with Angle) is proposed based on partial angle stack seismic gathers. This method can be used for direct inversion for the fluid factor, shear modulus and density of heterogeneous reservoirs. Firstly, an FMR approximation equation of a reflection coefficient is derived based on poroelasticity with P‐ and S‐wave moduli. Secondly, a stable simultaneous AVA inversion approach is presented in a Bayesian scheme. This approach has little dependence on initial models. Furthermore, it can be applied in heterogeneous reservoirs whose initial models for inversion are not easy to establish. Finally, a model test shows the superiority of this FMR‐AVA inversion method in stability and independence of initial models. We obtain a reasonable fluid factor, shear modulus and density even with smooth initial models and moderate Gaussian noise. A real data case example shows that the inverted fluid factor, shear modulus and density fit nicely with well log interpretation results, which verifies the effectiveness of the proposed method. 相似文献
4.
We develop a two‐dimensional full waveform inversion approach for the simultaneous determination of S‐wave velocity and density models from SH ‐ and Love‐wave data. We illustrate the advantages of the SH/Love full waveform inversion with a simple synthetic example and demonstrate the method's applicability to a near‐surface dataset, recorded in the village ?achtice in Northwestern Slovakia. Goal of the survey was to map remains of historical building foundations in a highly heterogeneous subsurface. The seismic survey comprises two parallel SH‐profiles with maximum offsets of 24 m and covers a frequency range from 5 Hz to 80 Hz with high signal‐to‐noise ratio well suited for full waveform inversion. Using the Wiechert–Herglotz method, we determined a one‐dimensional gradient velocity model as a starting model for full waveform inversion. The two‐dimensional waveform inversion approach uses the global correlation norm as objective function in combination with a sequential inversion of low‐pass filtered field data. This mitigates the non‐linearity of the multi‐parameter inverse problem. Test computations show that the influence of visco‐elastic effects on the waveform inversion result is rather small. Further tests using a mono‐parameter shear modulus inversion reveal that the inversion of the density model has no significant impact on the final data fit. The final full waveform inversion S‐wave velocity and density models show a prominent low‐velocity weathering layer. Below this layer, the subsurface is highly heterogeneous. Minimum anomaly sizes correspond to approximately half of the dominant Love‐wavelength. The results demonstrate the ability of two‐dimensional SH waveform inversion to image shallow small‐scale soil structure. However, they do not show any evidence of foundation walls. 相似文献
5.
Assessing uncertainty in refraction seismic traveltime inversion using a global inversion strategy 
下载免费PDF全文
下载免费PDF全文 J. Tronicke 《Geophysical Prospecting》2015,63(5):1188-1197
To analyse and invert refraction seismic travel time data, different approaches and techniques have been proposed. One common approach is to invert first‐break travel times employing local optimization approaches. However, these approaches result in a single velocity model, and it is difficult to assess the quality and to quantify uncertainties and non‐uniqueness of the found solution. To address these problems, we propose an inversion strategy relying on a global optimization approach known as particle swarm optimization. With this approach we generate an ensemble of acceptable velocity models, i.e., models explaining our data equally well. We test and evaluate our approach using synthetic seismic travel times and field data collected across a creeping hillslope in the Austrian Alps. Our synthetic study mimics a layered near‐surface environment, including a sharp velocity increase with depth and complex refractor topography. Analysing the generated ensemble of acceptable solutions using different statistical measures demonstrates that our inversion strategy is able to reconstruct the input velocity model, including reasonable, quantitative estimates of uncertainty. Our field data set is inverted, employing the same strategy, and we further compare our results with the velocity model obtained by a standard local optimization approach and the information from a nearby borehole. This comparison shows that both inversion strategies result in geologically reasonable models (in agreement with the borehole information). However, analysing the model variability of the ensemble generated using our global approach indicates that the result of the local optimization approach is part of this model ensemble. Our results show the benefit of employing a global inversion strategy to generate near‐surface velocity models from refraction seismic data sets, especially in cases where no detailed a priori information regarding subsurface structures and velocity variations is available. 相似文献
6.
The problem of equivalence in direct current (DC) resistivity and electromagnetic methods for a thin resistive and conducting layer is well‐known. Attempts have been made in the past to resolve this problem through joint inversion. However, equivalence still remains an unresolved problem. In the present study, an effort is made to reduce non‐uniqueness due to equivalence using global optimization and joint inversion by successive refinement of the model space. A number of solutions derived for DC resistivity data using very fast simulated annealing global inversion that fits the observations equally well, follow the equivalence principle and show a definite trend. For a thin conductive layer, the quotient between resistivity and thickness is constant, while for a resistive one, the product between these magnitudes is constant. Three approaches to obtain very fast simulated annealing solutions are tested. In the first one, layer resistivities and thicknesses are optimized in a linear domain. In the second, layer resistivities are optimized in the logarithmic domain and thicknesses in the linear domain. Lastly, both layer resistivities and thicknesses are optimized in the logarithmic domain. Only model data from the mean models, corresponding to very fast simulated annealing solutions obtained for approach three, always fit the observations. The mean model defined by multiple very fast simulated annealing solutions shows extremely large uncertainty (almost 100%) in the final solution after inversion of individual DC resistivity or electromagnetic (EM) data sets. Uncertainty associated with the intermediate resistive and conducting layers after global optimization and joint inversion is still large. In order to reduce the large uncertainty associated with the intermediate layer, global optimization is performed over several iterations by reducing and redefining the search limits of model parameters according to the uncertainty in the solution. The new minimum and maximum limits are obtained from the uncertainty in the previous iteration. Though the misfit error reduces in the solution after successive refinement of the model space in individual inversion, it is observed that the mean model drifts away from the actual model. However, successive refinement of the model space using global optimization and joint inversion reduces uncertainty to a very low level in 4–5 iterations. This approach works very well in resolving the problem of equivalence for resistive as well as for conducting layers. The efficacy of the approach has been demonstrated using DC resistivity and EM data, however, it can be applied to any geophysical data to solve the inherent ambiguities in the interpretations. 相似文献
7.
本文针对地层横波慢度径向分层模型,分析了地层横波慢度的径向非均匀性对弯曲波频散的影响.基于径向非均匀与均匀模型之间弯曲波频散的差异,结合微扰法和Backus-Gilbert(BG)理论建立了反演横波慢度径向分布的方程,求取了地层横波慢度的径向分布.在无噪声和参数误差时,反演结果较好地反映了实际地层横波慢度的径向分布,当井孔流体或井外地层纵波慢度的选取误差在10%内变化时,反演结果基本保持不变;存在信噪比(SNR)为20 dB(信号的功率为噪声的100倍)或10 dB(信号的功率为噪声的10倍)噪声时,反演结果没有发生明显的改变,其相对误差基本控制在10%以内,可见噪声对反演结果的影响不大.以上反演结果说明,本文采用的结合微扰法和BG理论的反演方法来估测地层横波慢度的径向分布时,具有很好的鲁棒性,可以被用于现场了解井壁周围的地层性质. 相似文献
8.
The discussion deals with the effect of shear wave velocity uncertainties on 1D seismic ground response analysis. In particular, the paper refers to uncertainties deriving from the solution of the inverse problem in surface wave methods. We address some issues related to the evaluation of “equivalent” profiles from surface wave data, the inversion strategy and the numerical simulation of seismic site response. The pitfalls in the analyses point out the need for more refined studies to draw general conclusions on the subject. 相似文献
9.
A nonlinear method for multiparameter inversion of pre‐stack seismic data based on anisotropic Markov random field 下载免费PDF全文
下载免费PDF全文 Hongbing Zhang Qiang Guo Lifeng Liang Chenghao Cao Zuoping Shang 《Geophysical Prospecting》2018,66(3):461-477
Multiparameter inversion for pre‐stack seismic data plays a significant role in quantitative estimation of subsurface petrophysical properties. However, it remains a complicated problem due to the non‐unique results and unstable nature of the processing; the pre‐stack seismic inversion problem is ill‐posed and band‐limited. Combining the full Zoeppritz equation and additional assumptions with edge‐preserving regularisation can help to alleviate these problems. To achieve this, we developed an inversion method by constructing a new objective function that includes edge‐preserving regularisation and soft constraints based on anisotropic Markov random fields and is intended especially for layered formations. We applied a fast simulated annealing algorithm to solve the nonlinear optimisation problem. The method directly obtains reflectivity RPP values using the full Zoeppritz equation instead of its approximations and effectively controls the stability of the multiparameter inversion by assuming a sectionally constant S‐ and P‐wave velocity ratio and using the generalised Gardner equation. We substituted the inverted parameters, i.e., the P‐wave velocity, the fitting deviation of S‐wave velocity, and the density were inverted instead of the P‐wave velocity, the S‐wave velocity, and the density, and the generalised Gardner equation was applied as a constraint. Test results on two‐dimensional synthetic data indicated that our substitution obtained improved results for multiparameter inversion. The inverted results could be improved by utilising high‐order anisotropic Markov random field neighbourhoods at early stages and low‐order anisotropic Markov random field neighbourhoods in the later stages. Moreover, for layered formations, using a large horizontal weighting coefficient can preserve the lateral continuity of layers, and using a small vertical weighting coefficient allows for large longitudinal gradients of the interlayers. The inverted results of the field data revealed more detailed information about the layers and matched the logging curves at the wells acceptably over most parts of the curves. 相似文献
10.
Joint PP‐ and PSV‐wave amplitudes versus offset and azimuth inversion for fracture compliances in horizontal transversely isotropic media 下载免费PDF全文
下载免费PDF全文 We propose a robust approach for the joint inversion of PP‐ and PSV‐wave angle gathers along different azimuths for the elastic properties of the homogeneous isotropic host rock and excess compliances due to the presence of fractures. Motivated by the expression of fluid content indicator in fractured reservoirs and the sensitivity of Lamé impedances to fluid type, we derive PP‐ and PSV‐wave reflection coefficients in terms of Lamé impedances, density, and fracture compliances for an interface separating two horizontal transversely isotropic media. Following a Bayesian framework, we construct an objective function that includes initial models. We employ the iteratively reweighted least‐squares algorithm to solve the inversion problem to estimate unknown parameters (i.e., Lamé impedances, density, and fracture compliances) from PP‐ and PSV‐wave angle gathers along different azimuths. Synthetic tests reveal that the unknown parameters estimated using the joint inversion approach match true values better than those estimated using a PP‐wave amplitude inversion only. A real data test indicates that reasonable results for subsurface fracture detection are obtained from the joint inversion approach. 相似文献
11.
Surface‐wave inversion for a P‐velocity profile with a constant depth gradient of the squared slowness 下载免费PDF全文
下载免费PDF全文 Surface waves are often used to estimate a near‐surface shear‐velocity profile. The inverse problem is solved for the locally one‐dimensional problem of a set of homogeneous horizontal elastic layers. The result is a set of shear velocities, one for each layer. To obtain a P‐wave velocity profile, the P‐guided waves should be included in the inversion scheme. As an alternative to a multi‐layered model, we consider a simple smooth acoustic constant‐density velocity model, which has a negative constant vertical depth gradient of the squared P‐wave slowness and is bounded by a free surface at the top and a homogeneous half‐space at the bottom. The exact solution involves Airy functions and provides an analytical expression for the dispersion equation. If the ratio is sufficiently small, the dispersion curves can be picked from the seismic data and inverted for the continuous P‐wave velocity profile. The potential advantages of our model are its low computational cost and the fact that the result can serve as a smooth starting model for full‐waveform inversion. For the latter, a smooth initial model is often preferred over a rough one. We test the inversion approach on synthetic elastic data computed for a single‐layer P‐wave model and on field data, both with a small ratio. We find that a single‐layer model can recover either the shallow or deeper part of the profile but not both, when compared with the result of a multi‐layer inversion that we use as a reference. An extension of our analytic model to two layers above a homogeneous half‐space, each with a constant vertical gradient of the squared P‐wave slowness and connected in a continuous manner, improves the fit of the picked dispersion curves. The resulting profile resembles a smooth approximation of the multi‐layered one but contains, of course, less detail. As it turns out, our method does not degrade as gracefully as, for instance, diving‐wave tomography, and we can only hope to fit a subset of the dispersion curves. Therefore, the applicability of the method is limited to cases where the ratio is small and the profile is sufficiently simple. A further extension of the two‐layer model to more layers, each with a constant depth gradient of the squared slowness, might improve the fit of the modal structure but at an increased cost. 相似文献
12.
Tijmen Jan Moser 《Geophysical Prospecting》2014,62(3):475-490
We study the stability of source mechanisms inverted from data acquired at surface and near‐surface monitoring arrays. The study is focused on P‐wave data acquired on vertical components, as this is the most common type of acquisition. We apply ray modelling on three models: a fully homogeneous isotropic model, a laterally homogeneous isotropic model and a laterally homogeneous anisotropic model to simulate three commonly used models in inversion. We use geometries of real arrays, one consisting in surface receivers and one consisting in ‘buried’ geophones at the near‐surface. Stability was tested for two of the frequently observed source mechanisms: strike‐slip and dip‐slip and was evaluated by comparing the parameters of correct and inverted mechanisms. We assume these double‐couple source mechanisms and use quantitatively the inversion allowing non‐double‐couple components to measure stability of the inversion. To test the robustness we inverted synthetic amplitudes computed for a laterally homogeneous isotropic model and contaminated with noise using a fully homogeneous model in the inversion. Analogously amplitudes computed in a laterally homogeneous anisotropic model were inverted in all three models. We show that a star‐like surface acquisition array provides very stable inversion up to a very high level of noise in data. Furthermore, we reveal that strike‐slip inversion is more stable than dip‐slip inversion for the receiver geometries considered here. We show that noise and an incorrect velocity model may result in narrow bands of source mechanisms in Hudson's plots. 相似文献
13.
Multistep inversion workflow for 3D long‐offset damped elastic waves in the Fourier domain 下载免费PDF全文
下载免费PDF全文 We present a new workflow for imaging damped three‐dimensional elastic wavefields in the Fourier domain. The workflow employs a multiscale imaging approach, in which offset lengths are laddered, where frequency content and damping of the data are changed cyclically. Thus, the inversion process is launched using short‐offset and low‐frequency data to recover the long spatial wavelength of the image at a shallow depth. Increasing frequency and offset length leads to the recovery of the fine‐scale features of the model at greater depths. For the fixed offset, we employ (in the imaging process) a few discrete frequencies with a set of Laplace damping parameters. The forward problem is solved with a finite‐difference frequency‐domain method based on a massively parallel iterative solver. The inversion code is based upon the solution of a least squares optimisation problem and is solved using a nonlinear gradient method. It is fully parallelised for distributed memory computational platforms. Our full‐waveform inversion workflow is applied to the 3D Marmousi‐2 and SEG/EAGE Salt models with long‐offset data. The maximum inverted frequencies are 6 Hz for the Marmousi model and 2 Hz for the SEG/EAGE Salt model. The detailed structures are imaged successfully up to the depth approximately equal to one‐third of the maximum offset length at a resolution consistent with the inverted frequencies. 相似文献
14.
At present, near-surface shear wave velocities are mainly calculated through Rayleigh wave dispersion-curve inversions in engineering surface investigations, but the required calculations pose a highly nonlinear global optimization problem. In order to alleviate the risk of falling into a local optimal solution, this paper introduces a new global optimization method, the shuffle frog-leaping algorithm (SFLA), into the Rayleigh wave dispersion-curve inversion process. SFLA is a swarm-intelligence-based algorithm that simulates a group of frogs searching for food. It uses a few parameters, achieves rapid convergence, and is capability of effective global searching. In order to test the reliability and calculation performance of SFLA, noise-free and noisy synthetic datasets were inverted. We conducted a comparative analysis with other established algorithms using the noise-free dataset, and then tested the ability of SFLA to cope with data noise. Finally, we inverted a real-world example to examine the applicability of SFLA. Results from both synthetic and field data demonstrated the effectiveness of SFLA in the interpretation of Rayleigh wave dispersion curves. We found that SFLA is superior to the established methods in terms of both reliability and computational efficiency, so it offers great potential to improve our ability to solve geophysical inversion problems. 相似文献
15.
2.5D modelling approximates 3D wave propagation in the dip‐direction of a 2D geological model. Attention is restricted to raypaths for waves propagating in a plane. In this way, fast inversion or migration can be performed. For velocity analysis, this reduction of the problem is particularly useful. We review 2.5D modelling for Born volume scattering and Born–Helmholtz surface scattering. The amplitudes are corrected for 3D wave propagation, taking into account both in‐plane and out‐of‐plane geometrical spreading. We also derive some new inversion/migration results. An AVA‐compensated migration routine is presented that is simplified compared with earlier results. This formula can be used to create common‐image gathers for use in velocity analysis by studying the residual moveout. We also give a migration formula for the energy‐flux‐normalized plane‐wave reflection coefficient that models large contrast in the medium parameters not treated by the Born and the Born–Helmholtz equation results. All results are derived using the generalized Radon transform (GRT) directly in the natural coordinate system characterized by scattering angle and migration dip. Consequently, no Jacobians are needed in their calculation. Inversion and migration in an orthorhombic medium or a transversely isotropic (TI) medium with tilted symmetry axis are the lowest symmetries for practical purposes (symmetry axis is in the plane). We give an analysis, using derived methods, of the parameters for these two types of media used in velocity analysis, inversion and migration. The kinematics of the two media involve the same parameters, hence there is no distinction when carrying out velocity analysis. The in‐plane scattering coefficient, used in the inversion and migration, also depends on the same parameters for both media. The out‐of‐plane geometrical spreading, necessary for amplitude‐preserving computations, for the TI medium is dependent on the same parameters that govern in‐plane kinematics. For orthorhombic media, information on additional parameters is required that is not needed for in‐plane kinematics and the scattering coefficients. Resolution analysis of the scattering coefficient suggests that direct inversion by GRT yields unreliable parameter estimates. A more practical approach to inversion is amplitude‐preserving migration followed by AVA analysis. SYMBOLS AND NOTATION A list of symbols and notation is given in Appendix D . 相似文献
16.
反演瑞雷波频散曲线能有效获取地层横波速度和厚度.但由于其高度的非线性、多参数、多极值等特点,传统的全局搜索方法易出现收敛速度慢、早熟收敛及搜索精度低的问题.鉴于此,本文提出并测试了基于萤火虫优化算法（FA）和带惯性权重的蝙蝠优化算法（WBA）的新的瑞雷波频散曲线反演策略.在瑞雷波频散曲线反演中,FA全局搜索能力强,但后期搜索精度低,而WBA局部搜索能力强,搜索精度高,但易出现早熟收敛.故本文将二者结合,提出了一种新的优化策略,称其为WFBA,即在反演前期使用FA,后期使用WBA,很好地解决了FA后期搜索精度低及WBA早熟收敛的问题.本文首先反演了三个典型理论模型的无噪声、含噪声的数据,验证了WFBA对瑞雷波数据反演的有效性与稳定性.然后将WFBA与WBA、FA单独反演以及不含惯性权重的FBA和粒子群优化算法（PSO）反演的结果进行了对比,说明了WFBA相对于WBA、FA、FBA和PSO具有更稳定、收敛速度更快、求解精度更高等优点.最后,反演了来自美国怀俄明地区的实测资料,检验了WFBA对瑞雷波数据反演的实用性.理论模型试算和实测资料分析表明,WFBA很适用于瑞雷波频散曲线的定量解释,具有很高的实用性价值. 相似文献
17.
Analysis of different parameterisations of waveform inversion of compressional body waves in an elastic transverse isotropic Earth with a vertical axis of symmetry 下载免费PDF全文
下载免费PDF全文 In a multi‐parameter waveform inversion, the choice of the parameterisation influences the results and their interpretations because leakages and the tradeoff between parameters can cause artefacts. We review the parameterisation selection when the inversion focuses on the recovery of the intermediate‐to‐long wavenumbers of the compressional velocities from the compressional body (P) waves. Assuming a transverse isotropic medium with a vertical axis of symmetry and weak anisotropy, analytical formulas for the radiation patterns are developed to quantify the tradeoff between the shear velocity and the anisotropic parameters and the effects of setting to zero the shear velocity in the acoustic approach. Because, in an anisotropic medium, the radiation patterns depend on the angle of the incident wave with respect to the vertical axis, two particular patterns are discussed: a transmission pattern when the ingoing and outgoing slowness vectors are parallel and a reflection pattern when the ingoing and outgoing slowness vectors satisfy Snell's law. When the inversion aims at recovering the long‐to‐intermediate wavenumbers of the compressional velocities from the P‐waves, we propose to base the parameterisation choice on the transmission patterns. Since the P‐wave events in surface seismic data do not constrain the background (smooth) vertical velocity due to the velocity/depth ambiguity, the preferred parameterisation contains a parameter that has a transmission pattern concentrated along the vertical axis. This parameter can be fixed during the inversion which reduces the size of the model space. The review of several parameterisations shows that the vertical velocity, the Thomsen parameter δ, or the Thomsen parameter ε have a transmission pattern along the vertical axis depending on the parameterisation choice. The review of the reflection patterns of those selected parameterisations should be done in the elastic context. Indeed, when reflection data are also inverted, there are potential leakages of the shear parameter at intermediate angles when we carry out acoustic inversion. 相似文献
18.
Derecke Palmer 《Geophysical Prospecting》2010,58(4):561-575
The tau‐p inversion algorithm is widely employed to generate starting models with many computer programs that implement refraction tomography. However, this algorithm can frequently fail to detect even major lateral variations in seismic velocities, such as a 50 m wide shear zone, which is the subject of this study. By contrast, the shear zone is successfully defined with the inversion algorithms of the generalized reciprocal method. The shear zone is confirmed with a 2D analysis of the head wave amplitudes, a spectral analysis of the refraction convolution section and with numerous closely spaced orthogonal seismic profiles recorded for a later 3D refraction investigation. Further improvements in resolution, which facilitate the recognition of additional zones with moderate reductions in seismic velocity, are achieved with a novel application of the Hilbert transform to the refractor velocity analysis algorithm. However, the improved resolution also requires the use of a lower average vertical seismic velocity, which accommodates a velocity reversal in the weathering. The lower seismic velocity is derived with the generalized reciprocal method, whereas most refraction tomography programs assume vertical velocity gradients as the default. Although all of the tomograms are consistent with the traveltime data, the resolution of each tomogram is comparable only with that of the starting model. Therefore, it is essential to employ inversion algorithms that can generate detailed starting models, where detailed lateral resolution is the objective. Non‐uniqueness can often be readily resolved with head wave amplitudes, attribute processing of the refraction convolution section and additional seismic traverses, prior to the acquisition of any borehole data. It is concluded that, unless specific measures are taken to address non‐uniqueness, the production of a single refraction tomogram that fits the traveltime data to sufficient accuracy does not necessarily demonstrate that the result is either correct, or even the most probable. 相似文献
19.
Uncertainty analysis and probabilistic segmentation of electrical resistivity images: the 2D inverse problem 下载免费PDF全文
下载免费PDF全文 Juan Luis Fernández‐Martínez Shan Xu Colette Sirieix Zulima Fernández‐Muniz Joëlle Riss 《Geophysical Prospecting》2017,65(Z1):112-130
In this paper, we present the uncertainty analysis of the 2D electrical tomography inverse problem using model reduction and performing the sampling via an explorative member of the Particle Swarm Optimization family, called the Regressive‐Regressive Particle Swarm Optimization. The procedure begins with a local inversion to find a good resistivity model located in the nonlinear equivalence region of the set of plausible solutions. The dimension of this geophysical model is then reduced using spectral decomposition, and the uncertainty space is explored via Particle Swarm Optimization. Using this approach, we show that it is possible to sample the uncertainty space of the electrical tomography inverse problem. We illustrate this methodology with the application to a synthetic and a real dataset coming from a karstic geological set‐up. By computing the uncertainty of the inverse solution, it is possible to perform the segmentation of the resistivity images issued from inversion. This segmentation is based on the set of equivalent models that have been sampled, and makes it possible to answer geophysical questions in a probabilistic way, performing risk analysis. 相似文献
20.
目前瑞雷波多阶模式频散曲线反演中仅考虑数据的拟合,缺乏对模型的约束,不能很好地刻画地层间断面的问题,针对此问题,研究了瑞雷波多阶模式频散曲线稀疏正则化反演方法.正演模拟基于广义反射-透射系数法,数值计算上采用一种快速求根方法,与二等分方法相比,能够在很短的时间内达到最优的收敛效果;反演建模时采用L1范数正则化方法对模型... 相似文献