有限元-无限元耦合法在三维直流电和电磁数值模拟中的应用.pdf

摘要 目前在三维直流电阻率法有限元正演数值模拟中存在以下两个 问题：第一个问题是，三维地电模型的快速建立及属性添加还存在困 难。目前国内多数学者仍然依靠在文本文件中手工输入网格信息的方 式进行模型建立及前处理，此方法耗时巨大、容易出错、且无法建立 复杂的三维模型。第二个问题是，传统的混合边界条件虽然可以得到 较高精度的解，但是由于系统矩阵与电源位置相关，对于电源位置多 次改变的装置，如偶极一偶极装置，每改变一次电源位置就需重新形 成系数矩阵及解方程组，以此为基础的反演计算将十分耗时。目前广 条件，即在半无限边界上强制电位为零或其法相导数为零，或假定混 合边界条件中的距离项为常数，以上方案要求有限元离散区域必须取 得非常大以减小边界条件造成的误差，这必将增加节点数和计算量。 在三维可控源电磁法有限元正演模拟中，数据量大和解方程耗时 是主要问题。由于电磁法研究的问题尺度较大，目标体覆盖区域一般 为几百米至几千米，造成有限元网格数量众多。Dirichlet外边界条件 是目前易于实现且效果较好的外边界条件，但其要求外边界远离场源 和目标体，一般范围都要在数万米，这无疑在我们不关心的区域增加 了许多有限元网格。考虑到电磁法数值模拟中每个有限元结点包含多 个自由度，且待求解的未知数均为复数，使得最终形成的系统矩阵规 模庞大，方程组求解困难。 在本文中，我们通过对通用三维前后处理软件GiD进行简单的 二次开发，实现了可视化的快速建立复杂三维地电模型，并且实现了 前处理、计算和后处理的整合；通过引入无限单元与有限元相结合， 形成了有限元一无限元偶和算法，替代了传统的人工边界条件，解决 了截断边界条件误差大、节点多的问题。对GiD的二次开发仅需要 用简单的脚本语言编写用户自定义“问题类型，即可在图形化界面 中建立并输出可用于已有有限元计算程序的初始模型。Astley波包映 射无限元被用来将电位分布或电磁场分布延伸到无限元处，并使其在 无限远处衰减为零。我们还提出了一种全新的无限元形函数，它保持 了系数矩阵的稀疏对称性，并且与多种常见无限元形函数的对比显示 在精度和时间消耗上均占优势。最后，在三维直流电法和三维电磁法 正演模拟中，通过若干不同的模型计算验证了本文提出的有限元一无 限元耦合算法的正确性和实用性。 本文提出的三维模型的建立方法快速、方便，且具备通用性；提 出的有限元一无限元耦合算法可以有效减少有限元网格剖分区域的范 围，节省结点数，加快计算速度。在三维直流电法数值模拟中，在仅 包含测区的计算范围内即可得到与混合边界条件相当的计算精度；在 三维可控源电磁法数值模拟中，本方法可以将各个方向的边界范围均 缩小到几千米，计算结果与数万米边界范围下施加Dirichlet边界条件 的传统有限元法相差无几。 关键词模型，前处理，可视化，有限元，无限元 lI ABSTRACT thatbame aretwo the3Ddirectcurrent Currently,thereproblems finiteelementforward firstishowto (DC)resistivity modeling．The build 3D model and efficiently complexgeoelectrical effectively．Most arestillbuild3Dmodelswithina．txtfileand domesticscholars by manual istime andnot pron