共查询到10条相似文献，搜索用时 31 毫秒
1.
海洋可控源电磁法(MCSEM)三维正反演理论现如今已经成为地球物理学研究的热点和难点之一,准确、高效、稳定的正演计算是实现快速反演计算的基础.三维正演数值模拟技术的发展已相对成熟,一些学者已将研究如何提高正演计算效率的目光转移到研究如何提高线性方场组的计算速度.为了提高MCSEM的三维正演问题的计算效率,本文首先从频域三维海洋电磁控制方程出发,然后利用Yee氏交错网格有限体积法在三维空间离散方程组,并施以第一类Dirichlet边界条件获得大型稀疏复系数线性方程组,最后引入3种不同几何多重网格迭代算法求解该线性方程组.为了检验GMG算法的正确性,通过建立一维层状油气模型,将3种GMG算法计算结果与Kerry Key等开发的二维开源程序MARE2DEM计算结果进行对比,两种程序求解电场分布的曲线能够很好的吻合,表明GMG算法能正确求解海洋电磁正演问题,且两种程序求解的相对误差数量级在1以下,表明GMG算法具有较高的求解精度.为分析GMG算法的计算效率,我们首先想要模拟出一个更加真实的海洋地下环境,将沉积(背景)层电阻率设计为三轴各向异性,然后在此环境中建立三维海洋油气油气模型,实现MCSEM三维正演计算.通过改变网格数,实施3种GMG迭代算法与GCROT迭代算法求解,结果表明:GMG算法求解三维海洋可控源电磁正演问题算法稳定,计算效率高.GMG算法作为Krylov子空间迭代算法的预条件器求解三维海洋可控源电磁正演问题,不仅能加快求解速度,而且能提高算法的稳定性. 相似文献
2.
《应用地球物理》2017,(1)
为加快直流电阻率法三维正演模拟的计算速度,本文引入一种新型的代数多重网格算法一聚集代数多重网格算法(AGMG)。首先从直流电阻率法满足的电位二次场微分方程出发,采用七点有限差分格式进行离散,结合混合边界条件形成大型稀疏求解线性方程组;然后详细给出AGMG法聚集粗化的成对聚集算法及技术流程,并采用V循环AGMG预处理共轭梯度(CG)算法(AGMG-CG)求解线性方程组,最终实现直流电阻率三维正演模拟。通过典型地电模型数值模拟研究,并与成熟的直流电阻率三维正演模拟程序(3DDCXH)结果及解析结果对比验证了本文给出算法可行性和准确性。另外通过对不同剖分网格和不同模型的数值模拟,并与传统迭代算法(ILU-BCGSTAB、ILU-GCR、SSOR-CG)对比表明,AGMG-CG算法不论从迭代次数还是迭代时间上都有显著优势,同时具有近乎线性快速下降、迭代次数随网格大小增加而缓慢增加等优点。因此,本文算法具有收敛精度高、收敛快、迭代稳定等优点,为提高直流电阻率法三维正演模拟的计算效率提供了可能。 相似文献
3.
4.
为了计算带任意地形的各向异性介质中二维大地电磁响应,本文在非结构化网格的基础上,采用有限体积法,开发了二维大地电磁各向异性正演模拟的新算法.首先,从Maxwell方程出发,推导二维各向异性介质中大地电磁场的边值问题;然后,采用三角网格自动生成技术对求解区域进行非结构化网格剖分,进而构建节点中心控制体积单元,利用有限体积方法,得到求解边值问题的大型稀疏线性方程组;最后,利用Pardiso精确地计算了大地电磁响应值.三个各向异性模型的计算结果表明,本文开发的有限体积算法,不仅能够高精度求解带任意地形的大地电磁电导率各向异性问题,而且对于同一模型,该方法的计算消耗和精度都与有限单元法相当.因此,有限体积法是处理电磁法各向异性问题的一种有效方法. 相似文献
5.
6.
高精度、快速有效的正演模拟算法是三维电磁正反演的前提.为了提高海洋电磁三维数值模拟的精度和效率,本文提出利用基于Gauss-Lobatto-Chebyshev（GLC）基函数的谱元法进行海洋可控源三维电磁正演模拟.谱元法结合有限元法和谱方法的优点.我们通过应用伽辽金加权残差法离散二次电场矢量亥姆赫兹方程,在单元内选择混合阶GLC多项式的张量积作为高阶矢量插值基函数,在求解大型稀疏线性方程组时利用直接求解器进行快速求解,从而实现了三维海洋可控源电磁快速高精度正演模拟.一维和三维模型正演结果验证了本文算法的有效性和准确性.典型模型的数值结果表明谱元法是一种有效的三维海洋可控源电磁正演数值方法,能在稀疏网格剖分情况下获得精确的海洋电磁正演模拟响应. 相似文献
7.
为了克服空气层和地表耦合以及避免一次场计算,开发适合不同类型场源、不同应用范围的频率域三维正演模拟统一平台,本文从麦克斯韦基本方程出发,推导基于Lorenz规范条件的磁矢势和标势耦合方程;通过将不同类型场源分解成一系列短导线(电性)源组合,采用交错网格采样和有限体积技术对方程进行离散得到对称大型稀疏线性方程组,并采用Jacobi迭代预处理QMR(Quasi-Minimum-Residual,拟最小残差)算法进行求解,我们成功实现不同类型场源、不同应用范围的频率域电磁法三维正演模拟.通过层状模型下大地电磁法以及有限长接地导线和大回线磁性源激发下的电磁场响应模拟,并与一维解析解对比验证算法的有效性.进而,我们利用该算法平台的模拟结果对典型地电模型在不同场源激发下频率域电磁法响应特征进行对比分析.本文算法研究及实现为建立频率域电磁法三维正反演统一框架打下基础. 相似文献
8.
正演是电磁法勘探野外工作参数选取、室内资料处理与解释的基础,精确、稳定、高效的三维正演算法尤为重要.本文采取先求解拉普拉斯域电场、再由Gaver-Stehfest算法获得时间域磁场的思路,基于电场异常场Helmholtz方程实现了交错网格有限差分法和有限体积法对回线源瞬变电磁法的三维正演.通过对比低阻块状体的积分方程法、时域有限差分法、矢量有限单元法和SLDM法的数值解,验证了交错网格有限差分法和有限体积法的正确性.由于交错网格有限差分法、有限体积法和基于矩形块单元的矢量有限单元法将待求电场均定义在矩形块单元棱边上,因此三种数值算法可采用相同方法进行电场待求量编码、计算背景场和后处理.然而,与矢量有限单元法相比,交错网格有限差分法和有限体积法的系数矩阵更加稀疏,求解效率更高.通过对水平低阻板状体三维模型的数值模拟,我们发现本研究中交错网格有限差分法比有限体积法精度更高;再利用一维解析法求解相应三层层状地电模型的感应电动势,我们还发现两种数值算法和一维解析法计算的感应电动势等值线形状吻合程度高,只是数值范围略有差异. 相似文献
9.
针对由多环状天线槽构成的方位随钻电磁波测井仪器真实结构,基于柱坐标系下耦合势Helmholtz方程和三维有限体积法研究建立一套相应的三维电磁响应数值模拟算法.首先,采用柱坐标系下三维不规则Yee氏交错网格对计算区域进行剖分,保证了钻铤表面附近的剖分网格与非光滑刻槽钻铤表面同形;并利用电场延拓边界条件逼近钻铤表面的理想电导体边界条件,保证非光滑刻槽钻铤表面耦合势方程离散精度.在此基础上,应用标准均质化技术确定非均质网格的等效电导率和磁导率,并借助三维有限体积法对柱坐标系中耦合势方程以及环状天线槽内的发射天线进行离散,得到一个大型非对称稀疏代数方程,通过PARDISO直接求解器求解离散方程.最后,通过与柱状均匀地层中模式匹配算法所得的数值结果的对比验证该三维算法的有效性,并进一步通过数值结果研究考察金属钻铤和天线槽对方位随钻电磁响应的影响,以及层状TI地层模型中该仪器的响应特征. 相似文献
10.
面波频散反演的数值模拟 总被引：1，自引：0，他引：1
大量的数值模拟表明,ＳＶＤ和ＬＳＱＲＤ在面波频散网格反演两步法中的应用效率都很好。但ＳＶＤ可以用分辨矩阵、信息矩阵和协方差矩阵对解估计进行数学上客观有效的评价,而对于大型稀疏方程的求解,ＬＳＱＲＤ确是一种内存注小,计算速度快以及分辨抗噪能力都较强的算法。在现有计算机运算速度较快、内存可以扩弃较大的条件下,实测数据量不很大时,应采用ＳＶＤ算法进行线性反演。 相似文献