共查询到18条相似文献，搜索用时 201 毫秒
1.
积分方程法三维模拟井地电法并行算法研究 总被引：1，自引：1，他引：1
在积分方程法模拟三维电磁响应原理的基础上,对井地电法三维正演模拟并行算法做了研究。同时,还利用Fortran 6.5和消息传递接口(MPI)开发了井地电法三维正演模拟并行计算程序。在求解大型线性方程组时,采用的是稳定的双共轭梯度(BICGSTAB)算法,单机正演程序与并行正演程序的计算结果有着很好的一致性,这说明了并行算法的正确性,为进一步开发井地电法三维反演并行程序奠定了基础,也为井地电法三维正、反演程序实用化提供了技术支持。 相似文献
2.
ZTEM(Z轴倾子电磁法)是一种天然场源的频率域航空电磁法,其特点是磁场垂直分量在空中机载平台测量,磁场水平分量在地面的固定基站测量,具有勘探深度大、速度快、成本低、覆盖面积大等技术优势。本文实现了ZTEM二维有限差分正演和二维非线性共轭梯度(NLCG)反演算法。研究对象是倾子资料,反演过程中通过解"拟正演"问题来避开雅克比矩阵的直接计算。通过理论模型合成数据反演试算,验证了ZTEM倾子资料二维NLCG反演算法的稳定性与可靠性。与大地电磁(MT) TE模式阻抗资料反演结果进行对比,发现在异常体横向边界的约束方面,ZTEM倾子反演比MT阻抗反演更具优越性。 相似文献
3.
随着地球物理设备和探测技术的不断发展,快速处理大规模地球物理数据的需求也随之增长。为了解决三维重力数据密度反演的耗时问题,提出一种并行的预处理共轭梯度算法来提高计算效率。本文分别采用两种不同的预处理算子通过组合模型数据反演进行测试比较,并利用迭代残差和计算用时共同评价其加速效果。结果表明：对称逐次超松弛预处理方法比对角预处理方法反演计算速度快,密度结果更贴近实际模型;与传统串行的共轭梯度算法相比,本文并行预处理快速算法可以获得近19倍的加速比。将该算法应用于美国Vinton盐丘的实测重力数据中,反演结果能够很好地圈定出岩体的位置,验证了本文并行预处理共轭梯度法在三维重力数据快速反演中的高效性和可行性。 相似文献
4.
电阻率三维复杂结构的快速反演 总被引：3，自引：1，他引：2
运用共轭梯度迭代算法解三维反问题中的线性方程组,并结合求Jacobi矩阵 G 的Rodi算法,则每次反演迭代仅需一次正演计算,大大加快了计算速度,实现了直流电阻率三维快速反演。另外,由于避免了存储 G 和 G T G 所需的庞大存储量,以及在三维反演中加入光滑约束,有利于精细网格下的复杂模型反演。 相似文献
5.
近年来,瑞利波波形反演技术因其避开了常规频散曲线计算,直接进行波场计算和反演不再受水平层状介质理论假设的限制,得到广大学者的高度重视。但瑞利波波形反演过程中需要不断进行波场正演和逆推计算。另外,由于浅地表速度较小,模拟计算时需要较小的网格间距才能避免数值频散,这无疑大大增加了正演模拟的计算量。对于这一问题,通常采用并行化设计来提高正演模拟的计算效率。本文基于消息传递接口（MPI）并行有限差分算法,以区域分解思路将模型区间分解成若干子区域,各区域互相通信,共同完成对模型的正演计算。并详细给出了区域分解、坐标转换、区域通信、波场合并等并行方案中的具体实现方法和实现步骤。通过对弹性模型、Kelvin黏弹性模型和标准线弹性固体（SLS）黏弹性模型不同并行方案的计算结果进行分析,验证了本文并行方案的可行性和有效性。并行计算结果表明,与单处理器计算时间相比,增加处理器数目可以明显减少计算时间,但随着处理器数目的增加,不同处理器之间的通信时间也增大;因此,并行时需要选择合适的处理器数目。对于黏弹性介质模型,SLS黏弹性模型的并行计算效率优于Kelvin黏弹性模型。 相似文献
6.
三维电阻率法对反演的精度和速度的要求越来越高,而正演是反演的基础,因此直流电阻率三维正演计算的速度和精度是三维电阻率反演实用化的关键。这里利用对称超松弛预条件共轭梯度法(SSOR-PCG),求解有限差分法离散生成的大型稀疏线性方程组,预条件矩阵的选择大大降低了系数矩阵的条件数,结合矩阵的一维非零元素压缩存储模式,使得正演计算速度得以提高,而内存占用量明显减小。在直流电阻率三维正演中采用异常场法,提高了电源点附近的解的精度。利用编制的有限差分正演程序,对两层模型、垂直接触带模型和低阻异常体模型进行了数值模拟,计算结果表明该算法是可行的,且可以明显提高正演计算的速度和精度。 相似文献
7.
8.
设计了一个简单的三维模型:三维低阻异常体和高阻异常体位于一维层状介质模型中,以"十"字型和"米"字型观测剖面方式作为模拟方式,开展模型一维、二维、三维反演技术有效性对比试验,其中一维反演计算采用自适应正则化(ARIA)反演,二维反演计算采用非线性共轭梯度(NLCG)反演,三维反演计算采用REBOCC三维反演,将剖面下方不同的反演结果与原始模型进行对比分析,发现:在该模型条件下,一维和二维反演都能得到反映模型真实信息的结果,REBOCC三维反演更偏向重建原始模型的宏观电性特征;在一维反演结果中,对于深部的电性结构特征,通常TE极化模式的反演结果好于TM极化模式,TE/TM几何平均值反演结果介于前两者之间;在二维反演结果中,通常TM极化模式的反演结果好于TE极化模式,TE/TM联合模式反演结果与TM极化模型相当,甚至更好;模型REBOCC三维反演相比一维、二维反演更易受反演结果多解性影响,REBOCC三维反演结果偏向重建原始模型的宏观电性特征。 相似文献
9.
《物探化探计算技术》2015,(4)
近年来2.5维复电阻率反演取得了一些理论和应用方面的成果,但计算过程复杂、计算量巨大,特别是正演计算过程中计算量大导致了反演速度较慢。这里在2.5维复电阻率反演计算基础上,将反演算法中的正演过程的排列和频率循环分割成多个规模较小的程序段,并发执行正演过程中的每个排列和频率循环的求解计算过程,并根据并行计算的节点数目或者计算核心数目,动态调整每个节点或者计算核心的计算量,从而提高了算法的计算效率和适应性。 相似文献
10.
为探究三维地质体对二维反演结果地影响以及二维反演中极化模式的选择问题,设计了一系列二维、三维地电模型进行正演计算,并对正演结果进行二维反演计算。研究结果表明:相位受到三维畸变的影响较小,视电阻率较大;当三维异常体模型走向延伸较小时,TE模式的反演结果误差较大,TM模式的反演结果相对较好,但是反演得到了虚假的下覆构造,TE+TM联合反演结合了TE和TM模式反演的优点,可以较好地反映异常体的赋存位置及电阻率值;当三维异常体模型走向延伸长度增大时,三维数据的二维反演结果是可靠的,用TM和TE+TM模式反演都比较合理。 相似文献
11.
傅氏变换有限差分(即FFD)混合波场延拓法波动方程正演模拟精度高,但三维FFD因其计算复杂、数据存储巨大而难以推广。这里提出三维FFD一步法正演模拟并行算法,把一个串行的大规模频率循环分割成并行的多个小规模频率循环,并发执行多个波场延拓来提高运算速度,并通过由给定的资源设定并行数,以及适时重排三维数据,大大增强了算法的适应性和时效性。通过对比并行和非并行正演实验,验证了算法的有效性。 相似文献
12.
大地电磁场三维地形影响的矢量有限元数值模拟 总被引：1，自引：0，他引：1
地形对大地电磁场的影响规律研究对于实际大地电磁测深反演解释具有重要的意义。基于矢量有限元法进行了大地电磁三维地形响应数值模拟。首先与二维地形模拟结果比较,验证了算法的正确性;接着,分别对1个山峰地形和1个山谷地形进行了模拟;最后,比较了二维模拟与三维模拟结果。结果表明：三维阻抗张量的Zxy模式和Zyx模式都会受到地形影响;用二维地形模拟结果去近似三维地形响应,二维TE模式与真实三维响应结果差别很大,TM模式接近真实三维响应结果;如采用二维反演处理带地形的三维MT数据,TM极化模式反演更为合理。 相似文献
13.
14.
利用一维正则化反演进行大地电磁测深数据拟二维反演解释 总被引：1，自引：0，他引：1
大地电磁测深的反演问题是不适定的,其反演结果不稳定,且具有非唯一性。通过在目标函数中采用正则化方法,可以使得不适定反演问题具有稳定的反演结果,并改善解的稳定性和非唯一性问题。为了提高野外大地电磁测深数据的处理效率和初步解释的精度,提出了大地电磁测深数据的一维正则化反演进行拟二维反演解释方法。这里所述的大地电磁测深一维反演解释,与以往的解释方法不同,其思路首先用Bostick反演的深度来控制层参数,使反演计算的模型参数仅存在电阻率;最后采用阻尼高斯-牛顿算法进行反演计算,并将Bostick反演结果作为反演计算的初始模型。通过模型试算,结果表明其处理速度快、解释直观,对野外大地电磁测深数据进行初步反演解释是可行的。 相似文献
15.
目前，大地电磁测深的资料处理与正、反演软件都是FORTRAN语言编写的，由于FORTRAN语言本身的局限性，使得这个软件存在很多的缺陷，但是因为用FORTRAN语言编写的算法已经比较成熟，将它全部用其他语言移植也比较费时费力，为了避免资料的浪费，提高大地电磁测深的资料处理与交互式系统的开发效率，并且更加完善其功能，在开发大地电磁测深二维交互式正反演系统时，对大地电磁测深中的计算部分仍用现有的FORTRAN程序，对有关的界面与图形诉开发部分用VC＋＋语言来编写，以增强和体现其交互式的功能。 相似文献
16.
地球化学场及其在隐伏矿体三维预测中的作用 总被引：1，自引：1，他引：0
本文从物理学"场论"的角度介绍了地球化学场的概念,并从场源、场作用及地球化学指标分布三个基本要素出发阐述了地球化学场的扩散、对流-扩散的动力学机制。首先,按照采样介质,将地球化学场分为原生的岩石地球化学场,以及次生的土壤、水系沉积物、水文和气体地球化学场,描述了组成各类地球化学场的要素。其次,将地球化学场的分析方法概括为静态的空间结构分析和动态的时空结构分析两种,提出要以场的动力学机制为基础,利用地球物理学中的正、反演理论来研究四维时空中地球化学场的发展和演化。最后,探讨了地球化学场时空结构分析与三维地学模拟两者之间彼此补充和相互验证的关系,三维地学模拟构建的场源及空间介质模型,为地球化学场的正、反演提供了初始条件,地球化学场反演的结果又可用来修正三维地质模型;探讨了地球化学场与大数据分析技术间的关系,即采用大数据的"数据驱动"的思路来挖掘其与多元地学数据之间的隐性联系,探索其与成矿过程的关联性。地球化学场与三维地学模拟、大数据分析技术的结合将为隐伏矿体三维预测中地下成矿物质的分布和演化提供依据。 相似文献
17.
在三维电阻率的正演计算中往往涉及到快速、准确求解大型线性方程纽Ax=b的问题。通过采用有限差分法来构造出求解点电源三维地电场的大型稀疏对称线性方程组。并引入Lanczos迭代技术，构造出三对角阵方程组，然后采用正交分解法进行求解，它是Krylov子空间方法中的一种。与传统迭代算法相比，它占用内存少，收敛速度快且稳定。针对大型稀疏矩阵及MATLAB语言的特点，采用简单记录矩阵的非零元素值及其所在行、列值的方法存储大型稀疏矩阵，可大大节省机器内存，提高运算速度。理论分析和计算实例显示，此算法是地电三维正演计算的有效方法，为下一步的反演计算打好基础。 相似文献
18.
对理论模型合成的稀疏测线数据和面积性数据以及实测单一测线数据进行了三维反演。通过对比分析合成稀疏测线数据和面积性数据的三维反演结果发现,使用4个张量阻抗元素作为反演数据的稀疏测线数据的三维反演结果可以较好地反映三维目标体,虽然效果可能不如面积性数据的三维反演结果。实测单一测线数据的三维反演结果也较好地揭示测线周围地下基岩界面的三维起伏情况。合成数据和实测资料的三维反演结果表明:可以采用三维反演解释的方法对大地电磁稀疏测线资料进行反演解释。同时,从合成数据的三维反演结果分析可以看出,对大地电磁稀疏测线资料进行三维反演,应尽可能反演所有张量数据;对大地电磁面积性资料进行三维反演,可以只使用两个非对角元素数据,但使用所有张量数据的效果更好。 相似文献