共查询到20条相似文献,搜索用时 46 毫秒
1.
We develop an inversion procedure using the total variation (TV) regularization method as a stabilizing function to invert surface gravity data to retrieve 3-D density models of geologic structures with sharp boundaries. The developed inversion procedure combines several effective algorithms to solve the TV regularized problem. First, a matrix form of the gradient vector is designed using the Kronecker product to numerically approximate the 3-D TV function. The piecewise polynomial truncated singular value decomposition (PP-TSVD) algorithm is then used to solve the TV regularized inverse problem. To obtain a density model with depth resolution, we use a sensitivity-based depth weighting function. Finally, we apply the Genetic Algorithm (GA) to select the best combination of the PP-TSVD algorithm and the depth weighting function parameters. 3-D simulations conducted with synthetic data show that this approach produces sub-surface images in which the structures are well separated in terms of sharp boundaries, without the need of a priori detailed density model. The method applied to a real dataset from a micro-gravimetry survey of Gotvand Dam, southwestern Iran, clearly delineates subsurface cavities starting from a depth of 40 m within the area of the dam reservoir. 相似文献
2.
Seismic anisotropy provides important constraints on deformation patterns of Earth's material. Rayleigh wave dispersion data with azimuthal anisotropy can be used to invert for depth-dependent shear wavespeed azimuthal anisotropy, therefore reflecting depth-varying deformation patterns in the crust and upper mantle. In this study, we propose a two-step method that uses the Neighborhood Algorithm(NA) for the point-wise inversion of depth-dependent shear wavespeeds and azimuthal anisotropy from Rayleigh wave azimuthally anisotropic dispersion data. The first step employs the NA to estimate depthdependent VSV(or the elastic parameter L) as well as their uncertainties from the isotropic part Rayleigh wave dispersion data. In the second step, we first adopt a difference scheme to compute approximate Rayleigh-wave phase velocity sensitivity kernels to azimuthally anisotropic parameters with respect to the velocity model obtained in the first step. Then we perform the NA to estimate the azimuthally anisotropic parameters Gc/L and Gs/L at depths separately from the corresponding cosine and sine terms of the azimuthally anisotropic dispersion data. Finally, we compute the depth-dependent magnitude and fast polarization azimuth of shear wavespeed azimuthal anisotropy. The use of the global search NA and Bayesian analysis allows for more reliable estimates of depth-dependent shear wavespeeds and azimuthal anisotropy as well as their uncertainties.We illustrate the inversion method using the azimuthally anisotropic dispersion data in SE Tibet, where we find apparent changes of fast axes of shear wavespeed azimuthal anisotropy between the crust and uppermost mantle. 相似文献
3.
Borehole-wall imaging is currently the most reliable means of mapping discontinuities within boreholes. As these imaging techniques are expensive and thus not always included in a logging run, a method of predicting fracture frequency directly from traditional logging tool responses would be very useful and cost effective. Artificial neural networks (ANNs) show great potential in this area. ANNs are computational systems that attempt to mimic natural biological neural networks. They have the ability to recognize patterns and develop their own generalizations about a given data set. Neural networks are trained on data sets for which the solution is known and tested on data not previously seen in order to validate the network result. We show that artificial neural networks, due to their pattern recognition capabilities, are able to assess the signal strength of fracture-related heterogeneity in a borehole log and thus fracture frequency within a borehole. A combination of wireline logs (neutron porosity, bulk density, P-sonic, S-sonic, deep resistivity and shallow resistivity) were used as input parameters to the ANN. Fracture frequency calculated from borehole televiewer data was used as the single output parameter. The ANN was trained using a back-propagation algorithm with a momentum learning function. In addition to fracture frequency within a single borehole, an ANN trained on a subset of boreholes in an area could be used for prediction over the entire set of boreholes, thus allowing the lateral correlation of fracture zones. 相似文献
4.
Based on the long-wavelength approximation, a set of parallel vertical fractures embedded in periodic thin interbeds can be regarded as an equivalent orthorhombic medium. Rock physics is the basis for constructing the relationship between fracture parameters and seismic response. Seismic scattering is an effective way to inverse anisotropic parameters. In this study, we propose a reliable method for predicting the Thomsen’s weak anisotropic parameters and fracture weaknesses in an orthorhombic fractured reservoir using azimuthal pre-stack seismic data. First, considering the influence of fluid substitution in mineral matrix, porosity, fractures and anisotropic rocks, we estimate the orthorhombic anisotropic stiffness coefficients by constructing an equivalent rock physics model for fractured rocks. Further, we predict the logging elastic parameters, Thomsen’s weak parameters, and fracture weaknesses to provide the initial model constraints for the seismic inversion. Then, we derive the P-wave reflection coefficient equation for the inversion of Thomsen’s weak anisotropic parameters and fracture weaknesses. Cauchy-sparse and smoothing-model constraint regularization taken into account in a Bayesian framework, we finally develop a method of amplitude variation with angles of incidence and azimuth (AVAZ) inversion for Thomsen’s weak anisotropic parameters and fracture weaknesses, and the model parameters are estimated by using the nonlinear iteratively reweighted least squares (IRLS) strategy. Both synthetic and real examples show that the method can directly estimate the orthorhombic characteristic parameters from the azimuthally pre-stack seismic data, which provides a reliable seismic inversion method for predicting Thomsen’s weak anisotropic parameters and fracture weaknesses. 相似文献
5.
在长波长假设条件下,水平层状地层中发育一组垂直排列的裂缝构成了等效正交各向异性介质.各向异性参数与裂缝弱度参数的估算有助于非均匀各向异性介质的各向异性特征描述,而弹性逆散射理论是非均匀介质参数反演的有效途径.基于地震散射理论,我们首先推导了非均匀正交介质中纵波散射系数方程,并通过引入正交各向异性特征参数,提出了一种新颖的正交各向异性方位弹性阻抗参数化方法.为了提高反演的稳定性与横向连续性,我们发展了贝叶斯框架下的正交各向异性方位弹性阻抗反演方法,同时考虑了柯西稀疏约束正则化和平滑模型约束正则化,最终使用非线性的迭代重加权最小二乘策略实现了各向异性特征参数的稳定估算.模型和实际资料处理表明,反演结果与测井解释数据相吻合,证明了该方法能够稳定可靠地从方位叠前地震资料中获取各向异性特征参数,减小参数估算的不确定性,为非均匀正交介质的各向异性预测提供了一种高可靠性的地震反演方法. 相似文献
7.
The analytical solution and algorithm for simulating the electric potential in an arbitrarily anisotropic multilayered medium produced by a point DC source is here proposed. The solution is presented as a combination of Hankel transforms of integer order and Fourier transforms based on the analytical recurrent equations obtained for the potential spectrum. For the conversion of the potential spectrum into the space domain, we have applied the algorithm of the Fast Fourier Transform for logarithmically spaced points. A comparison of the modelling results with the power‐series solution for two‐layered anisotropic structures demonstrated the high accuracy and computing‐time efficiency of the method proposed. The results of the apparent‐resistivity calculation for both traditional pole‐pole and tensor arrays above three‐layered sequence with an azimuthally anisotropic second layer are presented. The numerical simulations show that both arrays have the same sensitivity to the anisotropy parameters. This sensitivity depends significantly on the resistivity ratio between anisotropic and adjacent layers and increases for the models with a conductive second layer. 相似文献
8.
重力数据的密度界面反演是位场数据解释中的一项主要工作, 在区域构造演化、深部莫霍面确定等领域的研究中发挥重要作用.近年来, 数据驱动的深度学习方法广泛地应用在地球物理数据处理与反演中, 本文提出一种基于深度学习U-net网络的重力数据密度界面反演方法.首先, 对半椭球体界面模型进行随机抽取和组合进而形成地下起伏界面数据集, 并基于Parker正演理论对界面数据集进行重力异常正演计算, 为深度学习网络模型的训练提供特征完备的数据源; 其次, 设计了基于U-net网络模型的深度学习界面反演算法, 在传统的损失函数基础上增加光滑损失项和过拟合抑制项, 提高重力界面反演结果的光滑性和收敛效率; 最后通过测试样本集进行反演预测, 验证建立深度学习网络模型的泛化性.本文通过理论模型和实际数据试验分析了本文方法在密度界面反演中的有效性和实用性, 基于改进损失函数约束的深度学习界面反演方法有效地提高了密度界面反演的收敛效率和计算稳定性. 相似文献
9.
Traditional traveltime inversion for anisotropic medium is, in general, based on a “weak” assumption in the anisotropic property, which simplifies both the forward part (ray tracing is performed once only) and the inversion part (a linear inversion solver is possible). But for some real applications, a general (both “weak” and “strong”) anisotropic medium should be considered. In such cases, one has to develop a ray tracing algorithm to handle with the general (including “strong”) anisotropic medium and also to design a non-linear inversion solver for later tomography. Meanwhile, it is constructive to investigate how much the tomographic resolution can be improved by introducing the later arrivals. For this motivation, we incorporated our newly developed ray tracing algorithm (multistage irregular shortest-path method) for general anisotropic media with a non-linear inversion solver (a damped minimum norm, constrained least squares problem with a conjugate gradient approach) to formulate a non-linear inversion solver for anisotropic medium. This anisotropic traveltime inversion procedure is able to combine the later (reflected) arrival times. Both 2-D/3-D synthetic inversion experiments and comparison tests show that (1) the proposed anisotropic traveltime inversion scheme is able to recover the high contrast anomalies and (2) it is possible to improve the tomographic resolution by introducing the later (reflected) arrivals, but not as expected in the isotropic medium, because the different velocity (qP, qSV and qSH) sensitivities (or derivatives) respective to the different elastic parameters are not the same but are also dependent on the inclination angle. 相似文献
10.
推导了各向异性介质中由弹性系数表示的方位动校NMO速度的具体表达式,表明各向异性介质中方位NMO速度程椭圆形状,并分别对具水平对称轴的横向各向同性介质(HTI)、正交介质和单斜各向异性介质及在不同的裂隙填充物的性质下方位NMO速度进行了计算,结果表明裂隙的存在对NMO速度的影响不仅与裂隙密度有关,还取决于裂隙填充物的性质.同时,研究表明对于裂隙型单斜各向异性介质,其方位NMO速度椭圆轴向并不象HTI介质和正交介质中的那样与自然坐标系的坐标轴一致,而是发生了一定角度的偏离,其大小与裂隙填充物的性质、两组裂隙密度的比值及裂隙间的夹角等因素有关,研究结果为进一步区分裂隙介质的类型及裂隙填充物的性质提供依据. 相似文献
11.
位场向下延拓是位场数据处理和反演中的重要运算,但是它的不稳定性影响了它在许多处理和反演方法技术中的应用.本文通过把位场向下延拓视为向上延拓的反问题,得到向下延拓的褶积型线性积分方程,再利用Fourier变换矩阵的正交对称特性,并结合矩阵的奇异值分解和广义逆原理,提出了一种稳定的不需要进行求逆运算的位场向下延拓广义逆方法——波数域广义逆算法,解决了位场大深度向下延拓的不稳定性问题.把这种方法用于三维理论模型数据和实际磁场数据的向下延拓获得了理想的结果. 相似文献
12.
对于包含有垂向裂缝的横向各向同性地层或含有多组正交裂缝的各向同性地层,正交各向异性介质模型是最简单的与实际地层相符的方位各向异性模型.本文对单层水平反射界面正交各向异性模型采用射线追踪法计算了全方位角变化的P波走时,时距曲线表现出强方位各向异性.采用小生境遗传算法,对三条成一定角度的测线的走时信息进行速度和各向异性参数反演.模型算例表明,此方法可以得到高精度的裂缝方位角、P波垂直速度和较高精度的Thomsen各向异性参数. 相似文献
13.
The work develops the approximation approach to solving the inverse MTS problem with the use of neural networks. The inverse
problem is considered in model classes of parametrized geoelectric structures, whose electric conductivity is controlled by
a few hundreds of macroparameters ( N ∼ 300). An approximate inverse operator of the problem is constructed for each model class as a neural network, whose coefficients
are determined in the process of training on a representative sample of standard examples of forward problem solutions. The
problem of determination of the model class of geolectric structures corresponding to the presented input MT data is solved
with the use of the neural network classifier constructed for the available set of model classes of structures. Regularizing
factors and errors of the neural network method are analyzed. The operation of the algorithm is illustrated by examples of
the 2-D inversion of synthetic MT data. 相似文献
14.
Anisotropy is often observed due to the thin layering or aligned micro‐structures, like small fractures. At the scale of cross‐well tomography, the anisotropic effects cannot be neglected. In this paper, we propose a method of full‐wave inversion for transversely isotropic media and we test its robustness against structured noisy data. Optimization inversion techniques based on a least‐square formalism are used. In this framework, analytical expressions of the misfit function gradient, based on the adjoint technique in the time domain, allow one to solve the inverse problem with a high number of parameters and for a completely heterogeneous medium. The wave propagation equation for transversely isotropic media with vertical symmetry axis is solved using the finite difference method on the cylindrical system of coordinates. This system allows one to model the 3D propagation in a 2D medium with a revolution symmetry. In case of approximately horizontal layering, this approximation is sufficient. The full‐wave inversion method is applied to a crosswell synthetic 2‐component (radial and vertical) dataset generated using a 2D model with three different anisotropic regions. Complex noise has been added to these synthetic observed data. This noise is Gaussian and has the same amplitude f? k spectrum as the data. Part of the noise is localized as a coda of arrivals, the other part is not localized. Five parameter fields are estimated, (vertical) P‐wave velocity, (vertical) S‐wave velocity, volumetric mass and the Thomsen anisotropic parameters epsilon and delta. Horizontal exponential correlations have been used. The results show that the full‐wave inversion of cross‐well data is relatively robust for high‐level noise even for second‐order parameters such as Thomsen epsilon and delta anisotropic parameters. 相似文献
15.
为提高水平层状介质中三维电磁波散射和逆散射数值模拟的效率,在对角张量近似(DTA)的基础上根据不同回代方式得到了求解积分方程的DTA1和DTA2两种近似. 这两种近似可以作为计算积分方程稳定型双共轭梯度快速Fourier变换(BCGS-FFT)算法的初始猜测值和预条件因子,从而形成效率更高的混合DTA-BCGS算法. 散射实例说明了DTA2的高精度和混合DTA-BCGS算法尤其是混合DTA2-BCGS算法的高效率. 由于DTA2近似程度更高,将DTA2与变型Born迭代反演方法(DBIM)相结合形成了一种对三维异常体进行重构的快速电磁波逆散射技术. 文中的逆散射实例说明所开发的逆散射技术对重构水平层状介质中的任意三维异常体是非常有效的. 相似文献
16.
地震预警是地震减灾工作的重要途径,而震级预估是整个地震紧急预警系统中重要且较为困难的一个环节.目前,世界上多个国家和地区都已建立了各自的地震预警系统,并且形成了特征频率(τp和τc等)相关和特征振幅(Pd等)相关的两类震级紧急预警的方法,但各有局限性.本文在已有的方法和理论基础上,运用机器学习算法,将日本KIK和KNET台网从2015年至2017年所记录到的843条地震目录,55426条记录作为全数据集,设计、训练出一套用于常见震级范围的机器学习震级预估模型.与已有方法的预估结果相比,机器学习方法不仅使预估的整体误差和方差下降,同时多台联合评估单一地震事件的截面方差也更低.本研究的结果显示了机器学习算法在震级紧急预估问题上具有较广阔的应用前景,同时也为较为复杂的深度学习类算法框架下端到端模型提供了实践基础和研究可能. 相似文献
17.
由所建立的三维qP波相速度表示式出发,导出并解析求解各向异性介质中的频散方程,得到三维各向异性介质中的相移算子,进而将以相移算子为基础的对称非平稳相移方法推广到各向异性介质,发展了一个三维各向异性介质的深度偏移方法. 文中使用的各向异性介质的速度模型与现行的各向异性构造的速度估计方法一致,将各向同性、弱各向异性及强各向异性统一在一个模型中. 所建立的各向异性介质对称非平稳相移波场延拓算子可以同时适应速度及各向异性参数横向变化;文中给出的算例虽然是针对二维VTI介质的,但所提出的算法同样适用于三维TI介质. 相似文献
18.
利用差异演化(Differential Evolution)非线性全局优化算法,设计了一种反演地壳速度模型和进行地震定位的方法,并给出了反演结果的具体分析.利用有限差分算法计算速度模型的走时场,可以节省大量的计算量,加快计算过程.反演得到的地壳速度模型和地震的震源参数可以直接用于地震层析成像研究,还可以利用地震层析成像得到的三维速度结构对地震重新定位,从而得到较为精确的震源参数.地壳速度模型的反演方法也可以用于三维速度结构的反演. 相似文献
19.
Despite impressive progress in the development and application of electromagnetic (EM) deterministic inverse schemes to map the 3-D distribution of electrical conductivity within the Earth, there is one question which remains poorly addressed—uncertainty quantification of the recovered conductivity models. Apparently, only an inversion based on a statistical approach provides a systematic framework to quantify such uncertainties. The Metropolis–Hastings (M–H) algorithm is the most popular technique for sampling the posterior probability distribution that describes the solution of the statistical inverse problem. However, all statistical inverse schemes require an enormous amount of forward simulations and thus appear to be extremely demanding computationally, if not prohibitive, if a 3-D set up is invoked. This urges development of fast and scalable 3-D modelling codes which can run large-scale 3-D models of practical interest for fractions of a second on high-performance multi-core platforms. But, even with these codes, the challenge for M–H methods is to construct proposal functions that simultaneously provide a good approximation of the target density function while being inexpensive to be sampled. In this paper we address both of these issues. First we introduce a variant of the M–H method which uses information about the local gradient and Hessian of the penalty function. This, in particular, allows us to exploit adjoint-based machinery that has been instrumental for the fast solution of deterministic inverse problems. We explain why this modification of M–H significantly accelerates sampling of the posterior probability distribution. In addition we show how Hessian handling (inverse, square root) can be made practicable by a low-rank approximation using the Lanczos algorithm. Ultimately we discuss uncertainty analysis based on stochastic inversion results. In addition, we demonstrate how this analysis can be performed within a deterministic approach. In the second part, we summarize modern trends in the development of efficient 3-D EM forward modelling schemes with special emphasis on recent advances in the integral equation approach. 相似文献
20.
当地震信号通过复杂地球介质时,地层除了表现为各向异性,还表现为内在的粘弹性特征.因此,为准确描述地震波在地球介质中的传播特征,理想的地球介质模型应该能够模拟岩石的各向异性特征和衰减特征.本文给出了各向异性粘弹性介质模型的波动方程及其差分格式,并利用有限差分法实现了地震波波场数值模拟.结果表明了该介质模型中地震波场特征与各向异性主轴方位和介质的粘滞性参数之间的关系. 相似文献
|