共查询到18条相似文献，搜索用时 171 毫秒
1.
可控源音频大地电磁法（CSAMT）正演具有计算耗时、精度低的特点;常规的反演方法边界处理简单,占用内存大,搜索方向也难以精准,严重影响反演效果。为了实现CSAMT快速高精度正演,改善边界处理方法以减少内存占用,优化搜索方向以减少迭代反演次数,提出基于边界约束有限内存的拟牛顿CSAMT一维反演方法。首先,在CSAMT正演中基于高斯积分与并行算法实现了一维介质模型条件下有限长导线源激发的水平电场快速高精度计算。其次,在反演方法中,分别构建了基于相对误差和绝对误差的两种目标函数,并引入光滑模型约束,采取自适应正则化策略更新正则化因子,并采用基于边界约束的有限内存拟牛顿算法（LBFGS_B算法）,实现了CSAMT一维准确快速反演。以含低、高阻薄层的多层地电模型的水平电场Ex振幅为反演数据进行直接反演和模型检验。结果表明,该方法能够反映出地层电阻率随深度变化的趋势,并对高、低阻薄层有较好的分辨能力,低阻薄层分辨效果更佳。实际CSAMT资料的反演剖面的分辨率优于常规的连续介质反演方法,地质效果明显,表明该方法具有良好的应用前景。 相似文献
2.
ZTEM(Z轴倾子电磁法)是一种天然场源的频率域航空电磁法,其特点是磁场垂直分量在空中机载平台测量,磁场水平分量在地面的固定基站测量,具有勘探深度大、速度快、成本低、覆盖面积大等技术优势。本文实现了ZTEM二维有限差分正演和二维非线性共轭梯度(NLCG)反演算法。研究对象是倾子资料,反演过程中通过解"拟正演"问题来避开雅克比矩阵的直接计算。通过理论模型合成数据反演试算,验证了ZTEM倾子资料二维NLCG反演算法的稳定性与可靠性。与大地电磁(MT) TE模式阻抗资料反演结果进行对比,发现在异常体横向边界的约束方面,ZTEM倾子反演比MT阻抗反演更具优越性。 相似文献
3.
将无网格局部Petrov-Galerkin算法用于大地电磁二维正演。介绍了该方法的基本原理;从大地电磁二维边值问题出发,利用子域法详细推导了与之对应的局部Petrov-Galerkin弱式方程,并用高斯积分法将其离散化。论述了无网格局部Petrov-Galerkin法较无单元Galerkin法及有限元法的优缺点,最后通过二维模型的计算验证了算法的有效性。 相似文献
4.
基于有限单元法的二维／三维大地电磁正演模拟策略 总被引：1，自引：0，他引：1
对于二维和三维大地电磁正演问题,有限单元法最后形成了一个线性方程组KX=p。方程组中的K是大型稀疏的带状对称复系数矩阵,其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难。不完全LU分解处理的BICGSTAB算法,可用于该线性方程组的求解,并且具有速度快,精度高,稳定性好等优点。为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下,设计了非均匀网格剖分。在程序编制中,因只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间。通过对二维模型和三维模型电磁响应的计算,验证了该算法的正确性。 相似文献
5.
交错网格有限差分算法以其高效、精确、实用等优点在地震波数值模拟中得到广泛应用。目前交错网格有限差分的精度已达到时间4阶、空间2N阶;然而在求空间三次导数时,差分格式实际上并未达到所谓的2N阶精度,而是采用了低阶的差分格式及差分系数,这样有利于提高大尺度空间正演时的计算效率;但从计算精度的角度考虑,有必要推导出准确的满足2N阶精度的交错网格有限差分格式及差分系数,以得到更高精度的正演结果。笔者利用Taylor公式展开首次推导出了可导函数任意次导数的任意偶数阶精度的差分近似式及相应的差分系数,从而完善了常规高精度交错网格有限差分算法。采用新推导的交错网格有限差分格式得到的正演波形与解析解进行了对比,证明了新推导的差分格式的正确性,并与常规差分格式的正演波形进行了比较,结果显示,新推导出的交错网格有限差分格式模拟结果稳定性好,精度更高。 相似文献
6.
本文研究了利用阻抗信息进行可控源电磁勘探有限内存拟牛顿法三维反演的技术。首先用理论模型来验证有限内存拟牛顿法反演的准确性和可行性。观测参数为复阻抗Zxy分量,采用交错网格有限差分方法计算模型响应,反演采用有限内存拟牛顿法。数值模拟结果表明：有限内存拟牛顿法反演迭代速度较快,每4 min迭代一次,拟合差由146.00下降到1.78,收敛稳定;异常体的位置与理论模型吻合较好,有效地验证了有限内存拟牛顿法可控源三维反演的正确性。为了进一步验证该方法的实用性,将其应用到隐伏钼矿可控源电磁勘探工作中。工区的反演结果显示：在工区北西段深部存在高阻异常,其上为低阻异常。截取过钻孔的3号测线发现,-500~-100 m的位置表现为低阻,东侧存在向上涌起的高阻。推测此低阻为矿化蚀变带,延伸较深。该异常与钻探资料揭示的钼矿脉一致,证明了反演结果的准确性。因此,利用可控源观测得到的阻抗信息进行有限内存拟牛顿法三维反演,可以获得可靠的三维电阻率分布。 相似文献
7.
8.
为了加快大地电磁三维正演的求解速度,本文将一种新型的代数多重网格算法——聚集多重网格（aggregation-based algebraic multigrid, AGMG）算法引入大地电磁三维正演模拟中。首先从准静态条件下的麦克斯韦方程出发,利用交错网格有限体积法进行离散,并采用第一类Dirichlet边界条件形成大型稀疏复线性方程组;然后阐述AGMG算法的粗化策略和套迭代技术,并实施3种不同的AGMG求解算法：1）传统的V循环AGMG算法;2）AGMG预处理共轭梯度（AGMG-CG）;3）AGMG预处理广义共轭残差法（AGMG-GCR）。最终实现大地电磁法三维正演模拟。对典型地电模型进行正演模拟,并与已有的大地电磁三维正反演程序（ModEM）进行结果对比,以验证本文算法的准确性。另外,不同剖分网格和极化方式正演模拟结果与准残量最小化（QMR）迭代算法的对比表明,AGMG预处理求解算法（AGMG-CG、AGMG-GCR）不仅能够改善算法的稳定性,而且能够快速有效地求解正演问题;其中AGMG-GCR迭代次数更少,求解速度更快,误差衰减曲线更光滑,在144×152×104网格剖分情况下,相对于现有ModEM程序能够提高十几倍的计算速度,尤其适合大规模大地电磁三维正演问题。 相似文献
9.
《物探与化探》2015,(3)
介绍了有限差分算法在大地电磁测深法中的应用,推导了一维层状介质和连续介质模型下的有限差分算法,提出了一种新的网格剖分方法,并通过两层和三层介质模型的正演计算,验证了算法的正确性以及网格剖分方法的合理性。设计了一种连续介质模型,与两种不同程度近似的多层层状介质模型进行了正演结果的对比,指出了研究连续介质模型正演的必要性。对于连续介质模型,计算中剖分得越详细,其结果越接近真实值,而随着误差的逐渐减小,精确度的提高不再明显,但资源与时间的消耗将大幅增加。因此,无限制地减小网格步长是不必要的。同时,为了兼顾误差大小与网格步长的关系,采用了一种将一次网格剖分进行二次差分计算的方法得到新的结果,使计算结果明显得到了改善。 相似文献
10.
利用有限元方法进行大地电磁正演数值模拟时,由于是在有限网格区域上的数值计算,模拟计算时的网格边界为截断边界,而有限元数值模拟时的大地电磁场边界条件需要在足够远处才能够满足,所以截断边界的存在可能会使大地电磁正演模拟的边界条件无法满足,致使对计算结果和计算精度产生影响。利用有限元二维正演程序,在网格边界处加载一维情况下的大地电磁场,然后固定研究区域的网格剖分,并对一维地电模型和二维地电模型在改变有限元网格边界大小的情况下进行计算。在对一维模型进行模拟计算时,截断边界对边界条件没有影响,边界条件自然满足。而对二维模型进行模拟计算时,截断边界的存在对计算结果有较大影响。利用趋肤深度作为有限元网格边界变化的量度,通过改变网格边界大小,对不同的二维地电模型进行计算比较,总结出适合大地电磁有限元正演模拟的参考网格边界。 相似文献
11.
大地电磁正演的精度对后续反演意义重大。决定大地电磁有限元正演精度的因素,主要有插值函数、方程组的求解误差,以及网格的剖分程度。而网格剖分的合适与否,为有限元的求解提供了先决条件。因此,从边界条件和横向网格以及纵向网格三方面,讨论了网格剖分对大地电磁正演精度的影响。经研究表明,对于一维和左右对称的二维介质,左右和下边界都无需放置在无穷远处,并且横向网格对精度影响不大。而由于纵向网格与近地表的最小趋肤深度密切相关,所以地表网格不得大于1/4个趋肤深度。 相似文献
12.
13.
采用非结构三角形网格结合快速Occam算法对海洋可控源电磁(CSEM)与大地电磁(MT)数据进行二维联合反演试验研究。非结构三角形网格能够准确地模拟起伏地形和复杂地质构造,对反演目标区域采用精细网格剖分,其他区域采用粗网格剖分,在满足精度的前提下减少了不必要的计算量。将CSEM与MT数据加入到同一反演数据集中,通过对联合反演数据权重公式进行推导,构建CSEM及MT数据相关权重因子,控制不同数据的拟合权重来实现联合反演。最后对不同模型进行了反演计算,结果表明联合反演比单一反演对海底构造与异常体的还原度更高,验证了联合反演算法的可靠性。 相似文献
14.
基于一阶速度-应力波动方程,采用高阶交错网格有限差分数值模拟方法,对弹性及黏弹性TTI介质进行正演数值模拟。模拟时采用完全匹配层吸收边界条件(PML)消除边界反射。同时设计了层状介质模型、断层模型,通过模型的正演计算,得到了不同时刻的地震波波场快照及合成地震记录,分析其波场运动学及动力学特征。模拟结果表明,交错网格有限差分法可以很好地完成对复杂介质的波场模拟,具有较高的精度和可靠性。 相似文献
15.
16.
为提高电磁波场数值计算精度,对时间域有限差分法与时间域伪谱法进行了对比研究。时间域有限差分法是一种目前流行的电磁场时域数值计算方法,已被广泛应用于求解与时间有关的偏导数方程。对于大规模数值计算,时间域有限差分法需要较多的内存空间。时间域伪谱法基于时间域有限差分法,该方法使用快速傅里叶度换来计算麦克斯韦方程中的空间导数。由于傅里叶变换的准确性,时间域伪谱法使计算精度提高,数值计算时所需的格子数大大减少,这极大地节省了计算机内存空间,适合于大规模正反演问题的数值计算。 相似文献
17.
频率-空间域正演模拟是频率域及Laplace-Fourier域全波形反演的基础,起伏地表条件下波形反演算法的关键是正演算法中考虑起伏地表的影响。基于带PML吸收边界的声波波动方程,在已有最优9点有限差分正演算法的基础上构建了起伏地表条件下频率-空间域正演算法。通过应用变网格技术,进一步提高算法的计算效率、降低内存开销,使得大规模起伏地表模型的频率域正反演问题成为可能。理论分析及数值测试表明：通过对近地表区域进行局部网格加密,可有效地压制由于矩形网格离散引起的角点散射;结合变网格技术可较易获得5倍以上计算效率的提高及内存占用的降低,且随着模型尺度的增加及地表起伏高程差的减小,倍数将显著增加;在细网格与粗网格交界处产生的虚假反射振幅幅值控制在原始波场的2%以内,满足地震波场正反演的需求。 相似文献
18.
三维电阻率法对反演的精度和速度的要求越来越高,而正演是反演的基础,因此直流电阻率三维正演计算的速度和精度是三维电阻率反演实用化的关键。这里利用对称超松弛预条件共轭梯度法(SSOR-PCG),求解有限差分法离散生成的大型稀疏线性方程组,预条件矩阵的选择大大降低了系数矩阵的条件数,结合矩阵的一维非零元素压缩存储模式,使得正演计算速度得以提高,而内存占用量明显减小。在直流电阻率三维正演中采用异常场法,提高了电源点附近的解的精度。利用编制的有限差分正演程序,对两层模型、垂直接触带模型和低阻异常体模型进行了数值模拟,计算结果表明该算法是可行的,且可以明显提高正演计算的速度和精度。 相似文献