共查询到20条相似文献，搜索用时 287 毫秒
1.
《地球物理学进展》2020,(3)
瞬变电磁场数字滤波算法计算时间短,精度高,但难以实现三维正演计算;有限差分算法可以实现三维正演计算,但计算时间长,当模型网格剖分数量不足时计算精度较低.结合三维有限差分正演算法和一维数字滤波正演算法各自优势,提出双模型三维正演方法,可以减少计算时间,提高计算精度.首先,给出了瞬变电磁场三维有限差分正演算法和一维数字滤波正演算法.然后,通过理论分析和公式推导,提出了通过三维有限差分正演算法计算异常场,通过一维数字滤波算法计算背景场,然后叠加得到总场的双模型方法,并给出了具体计算公式和流程.最后,建立了经典的Newman均匀介质含低阻长方体模型,以及Commer层状介质含垂直接触带模型,分别采用双模型方法和常规三维有限差分算法进行了正演计算,对计算结果进行了对比分析.结果表明:在相同模型网格数量条件下双模型方法计算精度显著高于常规三维有限差分算法;双模型方法可以通过较少的模型网格数量取得高精度的计算结果,从而可以在保证计算精度的前提下显著提高计算效率. 相似文献
2.
采用非线性最小二乘法结合蒙特卡罗法,实现宽带高频电磁场椭圆极化率数据的精确反演,确定地下层状介质的真实电阻率和介电常数.反演结果表明,对于均匀半空间和二层介质模型,最小二乘法能够很好地实现反演,而对于三层或更多层的介质,首先利用蒙特卡罗法确定拟合初始模型,再进行最小二乘反演,能够避免收敛到局部极小值,提高了反演的稳定性.为了加速正演响应函数的计算和迭代的速度,采用高密度采样的线性滤波算法,大大加快了该精确反演方法的速度.针对如覆盖区地质填图和土壤调查等大面积确定地质体性质的应用,本文还给出了一种近似反演方法（相位矢量图法）,能够快速获取视电阻率和视介电常数,不仅可以为应用提供有用的基础信息,而且可作为精确反演方法的初始模型. 相似文献
3.
标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论:采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。 相似文献
4.
考虑到可控源音频大地电磁法(CSAMT)电偶极发射源与地下介质的三维结构特点,本文采用非结构化网格剖分技术,开展了三维CSAMT方法有限元数值模拟研究,将三维电磁场的背景场和异常场分别求解,避免了电偶极发射源的奇异性问题,并减小了计算区域.推导了三维异常电场遵循的有限元方程,加入散度条件进行约束以消除电场伪解;对非结构化网格单元采用高斯加权平均算法,得到了精度较高的异常磁场.针对层状介质模型,与积分方程法对比,验证了有限元算法的正确性;计算分析了典型三维地质模型的电磁响应,异常体反映明显.结果表明本文算法正确、可靠,适用于三维地质模型的CSAMT方法正反演研究. 相似文献
5.
G-S变换的快速算法 总被引：18，自引：8，他引：10
在电磁场瞬变响应的数值计算中 ,常采用G S变换法作逆拉氏变换 .它是纯实数运算 ,而且只需对较少的拉氏变换变量s值作计算 （通常对每一采样时间选用 1 2个s值 ） ,因而是一种计算速度较快的算法 .但是 ,要对大量采样时间作计算 ,其计算量仍太大 .本文基于拉氏变换的延迟定理 ,建立了一种新的G S变换算法 .数值检验结果表明 ,新算法可成级次地减少对大量采样时间作G S变换的计算量 ,显著提高电磁场瞬变响应的计算速度 . 相似文献
6.
复电阻率法在矿产、油气勘探调查中发挥着重要作用.为了认识复杂构造的复电阻率法电磁场的变化规律,本文基于自适应有限元方法,采用非结构化网格,引入Cole-Cole模型,实现了电偶源2.5D复电阻率法电磁场正演,可以模拟复杂地形和地电结构,正演结果更符合野外实际地质情况.通过将本文的计算结果与半空间模型解析解、层状介质和起伏模型结果进行对比,验证了本文算法的正确性.最后,基于复杂地电模型,通过正演模拟,系统分析了地形、激电参数、复杂构造对复电阻率法电磁场的影响特征. 相似文献
7.
井间电磁测量的2．5维层析成像方法 总被引：1，自引：0，他引：1
利用正则化最小二乘反演方法实现了井间电磁测量数据的层析成像,对井间地层电阻率进行了重建。在成像算法中,我们假设了井间电磁的激发与接收采用电磁偶极子源,井间介质仅在二维（xoz）平面内变化。在数值模拟中,通过对构造走向（y方向）的Fourier变换,将三维电磁场问题转化为一系列二维问题,用等参有限元方法在波数域求解,使实际地层模型的处理得以实现。对于波数域中每个波数对应的电磁场方程采用等参有限元求解,并用高斯积分将波数域解变换为空间域电磁场。利用源与接收器电磁场的互易原理,实现了电磁场响应对电导率分布灵敏度的快速计算。针对正演模拟中源点的奇异性,我们采用具有一定面积的伪艿函数表达源电流分布,使数值解精度得到提高。用层状介质的解析解与数值计算结果的对比,验证了模拟算法的精度。用介质扰动产生的电磁场变化检验互易性定理计算灵敏度的有效性。对简单块状模型、斜向裂缝带模型及“大”字模型的模拟数据成像结果表明,本文介绍的层析成像方法是正确有效的。 相似文献
8.
本文发展了积分方程法用于层状介质中三维不均匀体的大地电磁模拟算法（简称MT）.分为二个步骤:第一步,异常体用等效的散射电流代替,通过层状介质中的格林函数,建立以散射电流为未知参数的积分方程;第二步,把求得的散射电流乘上相应的格林函数,即得地面上的二次电磁场,由此而计算出各种MT响应. 文中采用了数值滤波与插值、群变换以及格林矩阵带状化三个方面的数值处理方法,提高了计算效率.通过与已发表的三维MT计算结果对比及格林函数互易性检验,表明了该算法的正确性.在此基础上,进行了数值模拟,初步讨论了三维MT曲线的畸变特点. 相似文献
9.
《地球物理学进展》2016,(1)
本文对瞬变电磁法2.5维有限元正演进行了研究.从频率域麦克斯韦方程组出发,经过傅里叶变换推导出了走向y方向频率域电磁场响应的变分问题,然后运用频率域和时间域的转换公式求解出时间域瞬变电磁场的解.在求解频率域电磁场响应时为提高精度采用了基于二次插值的高阶有限元的算法,即单元网格插值为二次函数,同时推导出了经有限元离散后的泛函问题;在求解时间域电磁响应采用了正余弦变换的数字滤波算法.通过基本模型的正演,验证了算法的可行性.同时,也对比了基于G-S变换的线性有限元算法的数值结果,结果表明,本文采用的算法精度更高,层状模型最大延迟采样时间提高到了100 ms以上. 相似文献
10.
非规则采样快速傅里叶变换（NFFT）主要用于快速计算非规则采样数据的频谱及重建.该方法为非规则采样数据频谱重建技术的核心算法.在实现NFFT算法时,高速度和高精度计算是其应用的前提和关键.本文针对二维NFFT计算效率,应用表驱动思路进行改进,将Gauss褶积算子由矩形改进为椭圆以减少计算量,将e指数计算改进为乘法以加快计算速度,并建表解决NFFT算法在地震资料处理中的应用问题.本文同时给出了非规则采样地震数据NFFT谱重建方法.最后本文给出算例验证提出方法的计算速度和精度,和非规则采样地震资料重建结果. 相似文献
11.
针对基于互相关成像条件的探地雷达（GPR）逆时偏移计算效率低、存储量大及易产生低频假象的不足,本文将激发振幅成像条件应用于GPR逆时偏移成像中.通过在源点电磁波场正向传播过程计算每个网格点的能量密度,并保存最大能量密度的时刻和相应的电磁波场值;在接收点电磁波场逆向传播过程提取每个网格点最大能量密度时刻及对应的电磁波场值,并利用保存的最大能量源点电磁波场及走时做归一化,从而获得了依赖反射系数成像剖面,避免了源点正向传播电磁波场的存储和重建.此外,为了提高电磁波场的模拟精度,采用了基于三角形剖分的时间域有限单元法（FETD）计算电磁波正向和逆向传播过程.最后通过模型试算表明：激发振幅成像条件相比于归一化互相关成像条件,成像结果低频噪声更弱,空间分辨率更高,计算效率提高了近2倍. 相似文献
12.
本文针对电导率横向同性地层中三轴发射-接收线圈系的井间电磁响应,开展三维正演算法及模拟研究.首先将发射线圈简化为三轴正交磁偶极子源.为克服低频电磁产生的数值迭代收敛缓慢问题,利用低感应数预处理方法将频域Maxwell方程组转化为基于矢势与标势的Helmholtz方程.采用非等间距的Yee氏交错网格、电导率的体积加权平均技术以及积分形式的有限差分格式实现对控制方程的离散化.借助于不完全LU分解预处理的稳定双共轭梯度法计算井间电磁响应的数值解,并通过与有限元软件COMSOL的数值结果对比验证本文模型与算法的有效性.数值计算结果表明：在横向同性地层的垂直井模型中,测井响应横向分量xx和yy能反映出地层电导率各向异性信息,但横向探测能力较差.轴向分量zz虽只反映地层水平电阻率信息,但具有较好的横向探测能力,且对地层中异常体非常敏感,而交叉分量xz、zx则具备较强的层边界识别能力. 相似文献
13.
作为一种基于节点的计算方法,无网格法具有构造高阶导数方便,自适应分析便利的优点,特别适合复杂地质构造的数值模拟.本文针对实际地球物理勘探中存在的起伏地形和各向异性的地电结构,提出用无网格法来模拟大地电磁响应,采用复合二次径向基函数构造了形函数,推导了大地电磁无网格法等价线性方程组,研究了系数矩阵的压缩存储方法以及大型稀疏复线性方程组快速求解算法,实现了起伏地形下各向异性的2D大地电磁高精度数值模拟.基于层状模型验证了算法的正确性,计算结果表明:无论是TM模式还是TE模式,计算相对误差均小于1%;通过对地垒和地堑模型的模拟,得出了起伏地形对视电阻率和相位的影响规律;对起伏地形下含有不同各向异性系数异常体的模型进行了数值模拟,为开展复杂地质结构的电磁场特征研究以及地形校正奠定了理论基础.随着计算科学的快速发展,无网格法必将成为新的高精度电磁场数值模拟方法. 相似文献
14.
Hui Sun Jian-Guo Sun Zhang-Qing Sun Fu-Xing Han Zhi-Qiang Liu Ming-Chen Liu Zheng-Hui Gao Xiu-Lin Shi 《应用地球物理》2017,14(1):56-63
3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method (FMM) is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method (WFC), which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models. 相似文献
15.
为了克服空气层和地表耦合以及避免一次场计算,开发适合不同类型场源、不同应用范围的频率域三维正演模拟统一平台,本文从麦克斯韦基本方程出发,推导基于Lorenz规范条件的磁矢势和标势耦合方程;通过将不同类型场源分解成一系列短导线(电性)源组合,采用交错网格采样和有限体积技术对方程进行离散得到对称大型稀疏线性方程组,并采用Jacobi迭代预处理QMR(Quasi-Minimum-Residual,拟最小残差)算法进行求解,我们成功实现不同类型场源、不同应用范围的频率域电磁法三维正演模拟.通过层状模型下大地电磁法以及有限长接地导线和大回线磁性源激发下的电磁场响应模拟,并与一维解析解对比验证算法的有效性.进而,我们利用该算法平台的模拟结果对典型地电模型在不同场源激发下频率域电磁法响应特征进行对比分析.本文算法研究及实现为建立频率域电磁法三维正反演统一框架打下基础. 相似文献
16.
利用积分方程法实现了均匀导电半空间三维大地电磁响应的数值模拟。求取张量格林函数积分时,采用二次剖分算法解决计算中奇异值问题,对于含有贝塞尔函数的积分项,利用结合连分式展开的高斯求积代替常规的快速汉克尔变换方法,确保了张量格林函数的正确计算并提高了计算精度。最后通过数值模拟结果的对比及模型试算验证了算法的正确性,所实现的三维大地电磁数值模拟算法为理论研究三维地电构造的大地电磁响应的分布规律提供了有效的工具,也为研究三维反演算法奠定了基础。 相似文献
17.
本文从瞬变电磁均匀半空间二次磁场响应公式出发,提出了一种基于神经网络的视电阻率快速计算方法.以中心回线为例,根据瞬变响应公式的特点,简化网络结构,选用三层BP神经网络和误差训练算法,用均匀半空间样本数据进行训练,确定了收敛快、误差小的一步正割法和隐含单元数,得到基于不同采样时窗的一组网络参数.用本文方法与二分法、牛顿迭代法做模型计算比较,及最后的实验计算,说明算法的快速,准确.本文方法不依赖初始模型,避开了复杂的电磁场数值计算,实现了视电阻率的快速计算,对瞬变电磁法资料的快速解释有一定的参考价值. 相似文献
18.
针对线性走时插值算法(LTI)不能正确追踪逆向传播射线的问题, 目前已提出多种改进算法, 如扩张收缩LTI算法、 循环计算LTI算法、 动态网络最短路径射线追踪算法等, 但这些算法的计算效率普遍偏低. 在分析各种改进LTI算法的优劣后, 本文提出了改进动态网络最短路径射线追踪算法. 该改进算法依据波的传播规律以及LTI算法的基本方程, 排除动态网络最短路径射线追踪算法中大量冗余节点计算, 并采用传统的二叉树堆排序算法对波前阵列节点进行管理. 数值算例表明, 本文提出的改进算法具有较高的计算效率, 其计算效率是动态网络最短路径射线追踪算法的4.5—30倍, 是原始LTI算法的2—6.5倍; 当动态网络最短路径射线追踪算法采用堆排序算法时, 改进算法的计算效率是其3.5—15倍. 相似文献
19.
高频电磁波测井同时求解视电导率和视介电常数 总被引：6，自引：0，他引：6
利用 Ｌａｍｂｅｒｔ Ｗ函数的性质,针对单发单收和单发双收两种电磁波测井仪器,给出 了一种求解视介电常数和视电导率的方法．理论分析和数值实验表明,该方法在计算速度和 精度方面,明显优于传统的求解复杂非线性方程组的方法．与查表方法相比,该方法不需要 预先计算大量图版,且精度较高． 相似文献