共查询到18条相似文献，搜索用时 851 毫秒
1.
多道瞬变电磁法(MTEM)是近年来发展起来的一种新的地球物理勘查技术,其数据采集方式与地震法类似,因此,采用瞬变电磁拟地震解释方法对MTEM数据进行处理解释具有一定的优越性.研究MTEM虚拟波场有效提取方法具有重要意义.在以往奇异值分解法、预条件正则化共轭梯度法两种不同的常规提取方法研究基础上,本文提出采用相关叠加法提取MTEM虚拟波场信息.首先采用预条件正则化共轭梯度法对全时段MTEM数据进行虚拟波场提取,然后采用同种方法对划分的各时间段数据进行虚拟波场提取,最后对全时段提取结果与各个时间段提取结果进行相关性叠加,叠加结果作为最终的提取结果.实测MTEM数据虚拟波场提取结果表明,采用相关叠加法可以得到稳定、光滑的虚拟波场波形曲线,且抗干扰能力强. 相似文献
2.
瞬变电磁扩散场到虚拟波场的精细积分变换算法 总被引：1，自引：0，他引：1
现有的瞬变电磁解释与反演方法主要对电阻率参数进行解释,较难给出准确的地质构造信息.为了能够较多获得地质构造信息,可将瞬变电磁扩散场信号转换为虚拟波场信号,提高瞬变电磁法的分辨率.本文采用精细积分法,实现了瞬变电磁扩散场到虚拟波场的转换.将一个高度病态的线性方程组转换成求积分的过程,大大降低了解决病态问题的难度.积分步长以2的指数增加,在保证精度的前提下,具有极高的计算效率;设置合理的终止迭代条件,进一步提高本文方法的适用性.对典型地电模型的计算,表明本文方法具有较高的精度;减小两个波峰之间的距离,证明了本文方法具有较好的分辨能力;对含有噪声信号的扩散场进行波场反变换,证明了本文方法具有较好的抗噪性.最后计算三维模型和实测数据的虚拟波场,可以很好地反映出三维地质体的界面信息,证明本文方法稳定、可靠,可以对地下地质目标体的界面进行有效识别.本文方法相比预条件正则化共轭梯度法(PRCG),在计算效率方面提升了 4倍左右;本文方法的相对误差在5％以内,而PRCG的最大相对误差可达40％～50％,精度最高可提升10倍;在分辨率与稳定性方面,也表现出较好的效果. 相似文献
3.
通过波场反变换从瞬变电磁信号中提取虚拟波场是定位电性界面和刻画地下导电目标形态的有效手段.波场反变换过程中,降低虚拟波场速度有利于改善虚拟波场对电性界面的分辨能力.然而,波场变换核函数随时间、虚拟时间以及虚拟波场速度变化的动态范围极大,且随虚拟波场降速而进一步扩大,造成波场反变换问题高度不适定,制约了虚拟波场提取的精度和效果.对此,在虚拟波场降速的条件下,本文以核函数在无穷区间上积分的解析解为准绳,通过最优化与数据拟合得到了经验公式,确定了波场反变换所需的最优积分区间;在适应采集时间与虚拟波场速度的同时,最大限度缩短了积分区间长度,降低了核函数变化的动态范围,压制了波场反变换方程的病态程度.此外,通过引入精细积分方法求解波场反变换方程,提高了波场反变换方程的求解精度.为验证方法的有效性与可靠性,对正演模拟的电性源瞬变电磁垂直磁场分量及其时间导数分量进行了波场反变换;结果回代后的拟合数据与正演数据的相对误差小于5%.在所提取的虚拟波场记录中观察到了合理的运动学特征,如实反映了低阻目标层深度、厚度和二次场观测时间范围的变化.在计算条件相同的前提下,通过降低虚拟波场速度,改善了波场反变换的精... 相似文献
4.
根据瞬变电磁场所满足的扩散方程与波动场所满足的波动方程之间的数学关系,运用扫时波场变换技术,得到稳定的虚拟波场.在虚拟波场的基础上,借用地震勘探中的Born近似逆散射技术实现对电性界面进行成像的目的.在理论推导中,直接从时间域波动方程出发,运用格林定理、引入均匀半空间的Green函数建立散射场关于速度扰动量的表达式;在数值求解中,运用奇异值分解的方法实现对第一类Fredholm积分方程的求解.通过对理论模型的计算以及实例应用,均表明瞬变电磁虚拟波场Born近似算法对地下电性界面具有良好的分辨能力. 相似文献
5.
在波场的正变换中,保证瞬变电磁场的计算精度的同时,应用两步最优化算法,成功地控制和减少了积分系数个数和离散数字积分的采样点个数,解决了在波场反变换式中,由于积分系数过多而产生的欠定方程组的问题,同时改善了第一类算子方程的不适定性.在波场反变换中,将正则化算法用于逆变换过程的计算中,通过采用偏差原理和Newton迭代格式选出最优的正则化参数,使得反变换所得到的波场稳定、可靠.通过对数值计算结果与已知波场函数对比,证明了该方法的有效性和实用性. 相似文献
6.
7.
8.
目前,瞬变电磁数据的主要处理结果是给出地下介质的(视)电阻率信息,较难给出清晰的地质界面信息.本文基于有限差分偏移成像原理,实现了电性源瞬变电磁虚拟波场的偏移成像技术.根据瞬变电磁扩散场与虚拟波场之间的数学关系,本文使用精细积分法对第一类Fredholm积分方程求解,然后消除直达波,得到只含有地层反射波的虚拟波场.在此基础上,本文从时间域波动方程出发,对上行波的一级近似方程反向外推,得到时间-空间域二维与三维的波场递推公式.通过求解地下空间的波场,并使用反褶积方法削弱虚拟波场的波形展宽效应,从而提高对地质界面的纵向分辨率,获得较为准确的地质界面信息.本文分别采用均匀半空间模型、D型模型、HK型模型以及两层含水采空区模型对本文方法进行验证,并对实测数据进行处理,偏移成像的结果与地质资料相吻合,说明瞬变电磁有限差分偏移成像方法在识别电性界面以及地层特征方面具有可行性. 相似文献
9.
10.
11.
A pseudo-seismic interpretation method is an alternative way to process and explain transient electromagnetic (TEM) data, and has become a popular research field in recent years. TEM signals which satisfy the diffusion equation can be converted by means of a mathematical transformation into ones which obey the wave equation. For an ill-posed problem of this kind of transformation, a sub-regularization algorithm is developed in this paper to extract a virtual wavelet of the TEM field. According to the conventional designation of TEM recordings, the entire integration period is divided into seven time intervals. In order to avoid low accuracy in the calculations, high-density wavefield data has been calculated based on the former sub-division. Therefore, the virtual wavelet can be extracted successfully by using an optimized algorithm to obtain high-density integral coefficients for all time windows, and a satisfactory condition number of the coefficient matrix while taking a different channel number in each time period. The Tikhonov regularization inversion scheme is used to determine the optimal parameters based on minimizing a least squares misfit, and the Newton iterative formula is used to obtain optimal regularization parameters. Both synthetic model simulations and a real data interpretation example indicate that the proposed pseudo-seismic wavefield method is a suitable alternative way to interpret TEM data. 相似文献
12.
13.
In this study, we apply Tikhonov’s regularization algorithm for a 3-D density inversion from the gravity-gradiometry data. To reduce the non-uniqueness of the inverse solution (carried out without additional information from geological evidence), we implement the depth-weighting empirical function. However, the application of an empirical function in the inversion equation brings the bias problem of the regularization factor when a traditional Tikhonov’s algorithm is applied. To solve the bias problem of regularization factor selection, we present a standardized solution that comprises two parts for solving a 3-D constrained inversion equation, specifically the linear matrix transformation and Tikhonov’s regularization algorithm. Since traditional regularization techniques become numerically inefficient when dealing with large number of data, we further apply methods which include the Simultaneous Iterative Reconstruction Technique (SIRT) and the wavelet compression combined with Least Squares QR-decomposition (LSQR). In our simulation study, we demonstrate that SIRT as well as the wavelet compression plus LSQR algorithm improve the computation efficiency, while provide results which closely agree with that obtained from applying Tikhonov’s regularization. In particular, the algorithm of wavelet compression plus LSQR shows the best computing efficiency, because it combines the advantages of coefficients compression of big matrix and fast solution of sparse matrix. Similar findings are confirmed from the vertical gravity gradient data inversion for detecting potential deposits at the Kauring (near Perth, Western Australia) testing site. 相似文献
14.
波形反演是一种利用全波场信息,通过最小化预测波场和实际波场的残差来揭示地下岩性和构造信息的方法.本文首先简述了常规拟牛顿算法的原理,之后利用一种新的拟牛顿公式对Davidon-Fletcher-Powell(DFP)和Broyden-Fletcher-Goldfarb-Shanno(BFGS)算法进行了修正,改进后的BFGS算法在近似Hessian矩阵逆矩阵时,不仅考虑了梯度和模型信息,还加入了目标函数本身的信息,而且对于每次迭代,基本没有增加计算量.数值试验表明,相对常规拟牛顿方法,修正BFGS算法在保证反演精度的同时,明显提高了反演效率. 相似文献
15.
通过对波场的时间二阶积分运算以增强地震数据中的低频成分,提出了一种可有效减小对初始速度模型依赖性的地震数据全波形反演方法—时间二阶积分波场的全波形反演方法.根据散射理论中的散射波场传播方程,推导出时间二阶积分散射波场的传播方程,再利用一阶Born近似对时间二阶积分散射波场传播方程进行线性化.在时间二阶积分散射波场传播方程的基础上,利用散射波场反演地下散射源分布,再利用波场模拟的方法构建地下入射波场,然后根据时间二阶积分散射波场线性传播方程中散射波场与入射波场、速度扰动间的线性关系,应用类似偏移成像的公式得到速度扰动的估计,以此建立时间二阶积分波场的全波形迭代反演方法.最后把时间二阶积分波场的全波形反演结果作为常规全波形反演的初始模型可有效地减小地震波场全波形反演对初始模型的依赖性.应用于Marmousi模型的全频带合成数据和缺失4Hz以下频谱成分的缺低频合成数据验证所提出的全波形反演方法的正确性和有效性,数值试验显示缺失4Hz以下频谱成分数据的反演结果与全频带数据的反演结果没有明显差异. 相似文献
16.
全波形反演是一种高精度的反演方法,其目标函数是一个强非线性函数,易受局部极值影响,而且反演过程计算量较大.波场重构反演是近几年提出的一种改进的全波形反演理论.该反演方法通过将波动方程作为惩罚项引入到目标函数中,通过拓宽解的寻找空间减弱了局部极小值的影响,而且反演过程不需要计算伴随波场,提高了计算效率.但该反演方法一直缺少准确的惩罚因子算法,直接影响到该方法的准确度.本文将波场重构反演拓展到时间域并利用梯度法进行波场重构.频率域的惩罚因子用来加强波动方程的约束,而时间域惩罚因子表现为调节模拟波场和实际波场的权重因子.为此,我们根据约束优化理论,在波动方程准确以及重构波场与反演参数解耦的假设下,提出以波动方程为目标函数的新的惩罚因子算法.根据波形反演在应用时普遍存在的噪音干扰、子波错误和低频信息缺失的情况下,应用部分Sigsbee2A模型合成数据对本文提出的算法进行实验.数值实验结果表明：基于新的惩罚因子算法,在其他信息不准确的情况下,波场重构反演可以给出高精度的反演结果. 相似文献
17.
准确模拟TTI介质中弹性波的传播是研究地震各向异性、AVO反演的基础. 在二维加权近似解析离散化(WNAD)算法的基础上, 本文发展的并行WNAD算法是一种研究三维横向各向同性(TI)介质中弹性波传播的、快速高效的数值模拟方法. 我们首先介绍三维WNAD方法的构造过程, 然后与经典的差分格式--交错网格(SG)算法进行了比较. 理论分析和数值算例表明, WNAD算法比交错网格算法更适合在高性能计算机上进行大规模弹性波场模拟. 同时, 本文利用并行的WNAD方法研究了弹性波在TTI介质中的传播规律, 观测了TI介质中弹性波传播的重要特征:横波分离、体波耦合和速度各向异性等. 在TTI介质分界面处, 弹性波产生更加复杂的折射、反射和波型转化, 使得波场非常复杂, 研究和辨别不同类型的波能够加深我们对由裂隙诱导的各向异性介质的认识. 相似文献
18.
叠前深度偏移理论及方法一直是地震数据成像中研究的热点问题.业界对单程波叠前深度偏移方法和逆时深度偏移开展了深入的研究,但对双程波方程波场深度延拓理论及成像方法的研究还鲜有报道.本文以地表记录的波场值为基础,利用单程波传播算子估计波场对深度的偏导数,为在深度域求解双程波方程提供充分的边界条件,并提出利用矩阵分解理论实现双程波方程的波场深度外推.通过对强速度变化介质中传播波场的计算,与传统的单程波偏移方法相比,本文提出的偏移方法计算的波场与常规有限差分技术计算的波场相一致,证明了本方法计算的准确性.通过对SEAM模型的成像,在相同的成像参数下,与传统的单程波偏移算法和逆时深度偏移算法方法相比,本文提出的偏移方法能够提供更少的虚假成像和更清晰的成像结果.本文所提偏移算法具有深度偏移和双程波偏移的双重特色,推动和发展了双程波叠前深度偏移的理论和实践. 相似文献