基于非结构化网格的大地电磁2D自适应有限元正演模拟

有限元求解大地电磁正演问题时,对研究区域的剖分常规做法是采用规则化的网格。但规则化网格在剖分地形、断层、褶皱等复杂模型会产生较大的几何离散误差。针对上述情况,这里采用非结构化的四边形网格对二维地电模型进行网格剖分,并与自适应有限元相结合,由剖分的粗网格出发,利用每个频点下网格单元的后验误差估计值指导网格的局部加密,优化网格质量和数量,从而提高正演模拟的精度。通过一维K型地电模型利用本文方法算出数值解与解析解进行对比分析,验证了自适应有限元法在求解大地电磁二维正演模拟中的有效性。并通过对断层、褶皱模型的正演模拟,分析了其大地电磁正演响应特征。     

一种电磁法三维数值模拟的建模与网格剖分方法

为满足电磁法三维数值模拟解释对交互建模与可视化网格剖分的需求,利用计算机图形学、人机交互、拓扑关系学等技术,设计了起伏地形下三维矢量交互建模、四面体和六面体网格剖分流程,并开发了软件模块。该软件可提供起伏地形下电磁法三维数值模拟的四面体、六面体网格剖分的实用化工具,提高了电磁法三维在数据模拟与解释方面的可视化程度和计算效率。     

基于非结构三角网格的海洋CSEM和MT二维联合反演研究

采用非结构三角形网格结合快速Occam算法对海洋可控源电磁(CSEM)与大地电磁(MT)数据进行二维联合反演试验研究.非结构三角形网格能够准确地模拟起伏地形和复杂地质构造,对反演目标区域采用精细网格剖分,其他区域采用粗网格剖分,在满足精度的前提下减少了不必要的计算量.将CSEM与MT数据加入到同一反演数据集中,通过对联...     

频率域航空电磁法人机交互式二维正演研究

研究了频率域航空电磁法的电磁场理论、边界条件、有限元单元法求解等内容,重点介绍了网格剖分原理,实现了人机交互网格剖分软件以及二维有限单元数值模拟算法。通过二维数值模拟计算,了解了圆柱体模型的响应曲线特征,提高了频率域航空电磁数据的解释水平。     

适合地形起伏的二维有限单元数值模拟技术

在有限单元法矩形～三角形单元剖分的基础上,提出了适合复杂地形条件下的网格剖分方法:平行四边形～三角形剖分方法。计算结果表明,该方法具有快速、高效的特点。只要对原有程序稍加修改,即可实现对复杂地形条件下的复杂地电断面进行自动剖分和网格节点自动编号。因此,对于复杂地形条件下的地电模型,该方法能快速地实现较为精确的数值模拟。     

基于非结构化网格的瞬变电磁2．5维有限元正演模拟

利用Delaunay三角化这种网格非结构化方法。通过编程实现了二维模型的非结构化三角形网格剖分,并编写了中心回线法瞬变电磁2．5维有限元正演程序。与前人计算结果对比,在取得相同计算精度的情况下,与结构化网格相比,非结构化网格所需网格和节点数量大大减少,计算效率更高。通过将非结构化网格法引入到瞬变电磁2．5维正演模拟中,实现了对复杂二维地电模型的有限元数值模拟,提高了现有有限元算法的应用范围。     

任意四边形剖分二维MT有限元正演模拟

为了模拟实际的起伏地形,在前人研究的基础上,单元网格设计为任意四边形网格,单元内的场值双线性变化,采用了高斯数值积分法计算单元系数矩阵,并给出二维大地电场的辅助场表达式。通过模型算例表明,模拟结果与解析解的均方误差在1%以内,地形模型与前人的模拟结果相吻合。分析对比了两个不同坡度模型对视电阻率和相位的影响。采用任意四边形网格剖分法,降低了编程难度,可以方便地适应野外地形的起伏变化。     

MT有限元模拟中截断边界的影响

利用有限元方法进行大地电磁正演数值模拟时,由于是在有限网格区域上的数值计算,模拟计算时的网格边界为截断边界,而有限元数值模拟时的大地电磁场边界条件需要在足够远处才能够满足,所以截断边界的存在可能会使大地电磁正演模拟的边界条件无法满足,致使对计算结果和计算精度产生影响。利用有限元二维正演程序,在网格边界处加载一维情况下的大地电磁场,然后固定研究区域的网格剖分,并对一维地电模型和二维地电模型在改变有限元网格边界大小的情况下进行计算。在对一维模型进行模拟计算时,截断边界对边界条件没有影响,边界条件自然满足。而对二维模型进行模拟计算时,截断边界的存在对计算结果有较大影响。利用趋肤深度作为有限元网格边界变化的量度,通过改变网格边界大小,对不同的二维地电模型进行计算比较,总结出适合大地电磁有限元正演模拟的参考网格边界。     

复杂电性结构大地电磁二维响应特征

为了改进计算区域离散化问题,本文利用自适应非结构化网格有限单元法求解二维地电结构下大地电磁场满足的加权余量表达式。在有限元求解电磁场的过程中,网格剖分越精细、计算精度越高,计算量也会越大。此外,结构化网格难以适应任意地形以及复杂地质构造。而自适应非结构化网格在电性变化剧烈的区域会自动加密,在电性缓变的区域则生成粗疏的网格,从而优化网格质量与数量。因此,文中引入COMSOL Multiphysics软件,以实现若干地电模型的构建及非结构化自由四边形单元网格化。将网格数据信息导入本文算法,计算大地电磁场响应,并与解析解及数值解对比。结果表明,基于非结构化网格的正演模拟精度高、适应性强,为计算区域网格化提供了新的方法。     

基于二次场二维起伏地形MT有限元数值模拟

通过计算二次场来进行二维大地电磁数值模拟;导出了二维大地电磁二次场的微分方程,利用有限单元法来解微分方程;对矩形网格进行对角线的二次剖分,更容易且真实地模拟起伏地形。对几个典型模型进行了试算,与前人总场法的计算结果做了比较,两者视电阻率曲线一致,证明本文算法是正确的;通过2个简单的算例说明复杂地表下2种极化模式的MT观测资料都有明显的异常,视电阻率在TM模式下比TE模式更易受地形影响,TE模式下视电阻率曲线形态与地形呈"正相关",TM模式下反之。     

MT二维正演的混合剖分算法及其在线源CSAMT正演计算中的应用

在矩形有限元和三角形有限元方法大地电磁测深（ＭＴ）二维正演计算的基础上，设计了矩形单元粗剖分、程序自动进行三角形单元细剖分的混合剖分算法程序，使得ＭＴ二维模型正演的计算速度大为提高。该方法同时被应用于线电流源变频电磁测深（ＣＳＡＭＴ）ＴＥ极化方式的二维正演计算。文章给出了在不同模型上计算的ＭＴ和线源ＣＳＡＭＴ的响应，分析了计算ＣＳＡＭＴ响应时，计算点离源的距离对视电阻率和相位值计算的影响，比较了ＭＴ和ＣＳＡＭＴ网格剖分的特点，提出了计算ＣＳＡＭＴ响应的改进的边界条件设计     

基于ILU-BICGSTAB算法的起伏地形MT二维正演模拟

针对起伏地形对异常体的探测效果的影响,利用有限单元法,基于MTLAB平台编写了大地电磁二维正演程序,为更好地模拟实际地形,采用矩形网格中进一步三角细化的剖分方式,并利用不完全LU分解的预条件BICGSTAB算法提高计算速度与精度,最后通过建立多种地电模型分析讨论了起伏地形对大地电磁场的影响,验证了该方法的实用性与准确性。     

基于非结构网格的带地形MT二维Occam反演及应用

本文对基于自适应非结构三角形网格的带地形MT数据二维Occam反演进行应用研究.自适应非结构三角形网格能够准确地模拟起伏地形和复杂地质构造,正演网格由有限元解的后验误差估计指导自动细化,保证了模型响应的准确性,反演网格在反演目标区域采用精细网格剖分,在模型边界区域采用粗网格剖分,在满足反演精度的前提下减少了不必要的反演...     

各向异性介质中大地电磁场的异常特征

给出各向异性介质二维地电断面大地电磁场的边值问题以及等价的变分问题。对计算区域采用矩形网格中进一步三角细化的剖分方式并在三角单元内进行线性插值,解出有限单元法数值解。通过典型模型的正演模拟,得到大地电磁测深曲线,并与前人的研究工作对比,验证了该方法的有效性。     

起伏地形条件下二维大地电磁各向异性正演

在大地电磁资料处理和解释中,大地电磁的各向异性正演一直是国内、外研究的前沿课题。这里首先从麦克斯韦方程组出发,得出各向异性介质二维大地电磁场的边值问题以及等价的变分问题,进而对计算区域进行完全非结构化三角形剖分,在单元内采用线性插值,将变分方程转化成线性方程组,求解出有限单元法数值解,获得各向异性介质下的电磁场值。讨论了电性主轴与笛卡尔坐标轴之间的夹角以及地形的起伏变化对视电阻率和相位值的影响,结果表明:起伏地形条件下,TE和TM两种极化模式的视电阻率和相位值都会受到影响,且TM模式受地形影响较大;由于介质的各向异性,视电阻率和相位断面图出现明显不同于各向同性的形态,导致横向分界面模糊。     

岩体结构三维无缝建模与四面体优化剖分

复杂地质特征给数值模拟的建模工作带来很大困难,针对该问题,提出并实现了岩体结构三维无缝建模方法。将直接法拓展到四面体网格剖分中,实现了基于岩体结构无缝模型的四面体优化剖分方法,从而为数值模拟提供既符合地质实际,又有较高网格质量的计算模型。利用该方法建立了某拱坝坝肩岩体三维地质模型,并进行四面体优化剖分。     

求解潜水含水层剖面二维稳定流问题探讨

本文在对现有求解潜水含水层剖面二维流的多种方法及其特点进行分析的基础上,探讨新的求解方法.提出的方法只需进行一次网格剖分,在网格节点编号保持不变的情况下,把潜水面与单元相交的形状概化为五边形、任意四边形和三角形,按五边形、任意四边形和三角形子程序准确地求出潜水面穿越的单元的渗透矩阵.同时不用计算完全在潜水面之上的单元的渗透矩阵.本文给出具体的思路、算法步骤,计算潜水面的位置,通过比较数值解与解析解,两者虽然存在一定的误差,但是该方法总体结果是比较令人满意的.     

电阻率层析成像及演中若干问题的讨论

以双异性点电流源电流场二维地电断面为例，着重讨论了电阻率展析成像反演中有限单元法数值模拟边界条件的选择、网格剖分以及模拟精度等问题，同时也讨论了成像反演中最优化及演方法及拟合参数的评价等。文中提出的新的网格剖分方法和选择的边界条件可有效地提高计算速度和数值模拟的精度。     

复杂区域四边形网格生成的一种改进方法

提出了一种改进的四边形网格生成算法。该法在行波法生成平面三角形网格的基础上，将一个三角形单元分解为三个四边形单元。经过优化处理，获得优良的计算域四边形网格。针对岩土工程结构问题，还提出了网格剖分中尖灭不连续面和锚杆及锚索的处理方案。地下洞室和边坡算例表明了这种方法的可行性和可靠性。     

航空大地电磁响应特征分析及反演研究

航空大地电磁法以其机动灵活、效率高等优点,适用于在地势崎岖的偏远环境中开展快速普查作业。在航空大地电磁探测中,起伏地形会对倾子响应产生影响,另外由于倾子是无量纲传输函数,只包含电性结构的空间变化信息,而与电阻率的数值水平无关,单独利用倾子数据进行反演时,反演多解性更大。为解决上述问题,首先采用四边形网格模拟地形起伏;其次在倾子数据的反演目标函数中加入少量地面大地电磁数据拟合项,克服只拟合倾子的多解性缺陷;最后利用二维地电模型,验证本文有限元算法的正演精度,并分析起伏地形对倾子响应数据的影响。针对倾子数据缺乏背景电阻率信息这一问题,开展了带地形的倾子反演和倾子与地面大地电磁数据反演研究。研究结果表明:向倾子数据中加入少量地面大地电磁数据进行反演,可有效降低单独倾子反演对初始模型的依赖。