基于有限元法的二维各向异性大地电磁场数值模拟研究.docxVIP

基于有限元法的二维各向异性大地电磁场数值模拟研究

一、引言

1.1研究背景与意义

在地球物理勘探领域，大地电磁法凭借其独特优势占据着至关重要的地位。该方法利用地球内部与大气层间的电磁场变化，依据电磁场在介质中传播时的反射、折射、散射等特性，结合地下储层矿物、岩石、水体等介质的电磁响应，能够获取地下介质的物理、化学和结构特征信息，广泛应用于寻找矿床、地下水、构造裂缝等地下信息的勘探工作，对推动资源勘探和地质研究有着关键作用。

实际的大地电磁场呈现出强烈的二维各向异性，即电阻率在水平和垂直方向上具有不同性质。这种特性在诸多地质构造和地下介质中普遍存在，例如片麻岩和其它类型的岩石具有明显的电导率各向异性，在断层破碎带和裂隙发育带上，呈脉状或条带状分布的岩体常常呈现出强烈的宏观电导率各向异性。在实际勘探中，若忽视这种二维各向异性，将导致对地下介质电磁响应的模拟出现偏差，进而影响对地下地质结构和资源分布的准确判断。因此，对二维各向异性进行精确建模和计算，是实现精准模拟地下介质电磁响应的核心需求。

当前，有限元法已成为大地电磁场数值模拟的重要方法之一。它能够将连续的物理问题区域划分为一系列较小的、有限的、并且相互连接的子域（即有限元），每个有限元内通过一组特定的近似函数来表示场变量的分布，进而建立关于未知量的代数方程组求解物理量分布。然而，针对二维各向异性大地电磁场的有限元数值模拟研究仍存在诸多局限。尤其是在面对复杂地质构造和近地表高阻抗等情况时，数值模拟的精度和计算效率问题亟待解决。高精度的数值模拟结果可以为地球物理勘探提供更为可靠的数据支持，助力勘探人员更准确地判断地下地质结构和资源分布情况，提升勘探效率和准确性，降低勘探成本和风险；高效的计算方法能够在更短时间内处理大量数据，满足实际勘探工作中的时效性需求。此外，对二维各向异性大地电磁场的深入研究，还能进一步推动大地电磁场数值模拟方法和技术的发展，为地球物理领域的理论研究和实际应用开辟新的道路。

1.2国内外研究现状

国外在二维各向异性大地电磁场有限元数值模拟领域起步较早，取得了一系列具有影响力的成果。Reddy等首次提出了二维情况下方位各向异性介质正演模拟的简单有限元方法，为后续研究奠定了重要基础；Pek等提出的任意双轴各向异性的有限体积算法，因对复杂各向异性情况的良好适应性，受到广泛认可和应用。在复杂地质构造模拟方面，一些研究通过构建精细的地质模型，利用有限元法模拟电磁场传播，分析各向异性对电磁场响应的影响，在油气勘探等领域取得了实际应用效果。

国内学者在该领域也积极开展研究并取得显著进展。黄一凡等利用基于非结构网格的自适应有限元法求解对角各向异性情况下大地电磁正演控制方程，既有利于起伏地形和复杂地质构造的建模，又便于模拟包含多种尺寸结构的地电模型，通过对经典模型的正演模拟试算和复杂海底起伏各向异性模型的反演试算，验证了方法的准确性和有效性；霍光谱等从麦克斯韦方程出发，引入张量电导率，使用有限差分法求解关于平行走向的电场分量和磁场分量的二阶偏微分方程，研究带地形的二维电性各向异性结构的正演模拟，分析地形因素对大地电磁响应的影响。

尽管国内外在该领域已取得一定成果，但仍存在一些不足之处。部分数值模拟方法在处理复杂地质构造和近地表高阻抗问题时，精度和计算效率难以兼顾，导致模拟结果与实际情况存在偏差；不同数值求解方法在精度和计算效率方面的对比研究还不够系统全面，缺乏针对特定复杂场景的优化策略；对模拟结果在地球物理勘探实际应用中的分析和验证还需进一步加强，以提高模拟结果的实用性和指导价值。

1.3研究内容与方法

本研究围绕二维各向异性大地电磁场的有限元数值模拟展开，具体内容涵盖以下几个方面：

构建包含二维各向异性的大地电磁场模型：基于国产有限元软件，深入研究二维各向异性的特性，建立精确的大地电磁场数学模型，并对模型进行合理的网格划分和边界条件设置，确保模型能够准确反映实际地质情况和电磁场传播条件。

探究二维各向异性电磁场传播特性：深入分析在地质构造变化和近地表电阻率变化等复杂情况下，二维各向异性电磁场的传播特性，研究其对电场和磁场响应的具体影响，揭示其中的内在规律。

比较不同数值求解方法：全面比较不同数值求解方法在二维各向异性大地电磁场数值模拟中的精度和计算效率差异，针对复杂场景下的数值模拟问题进行深入探讨，提出切实可行的改进措施和优化方案。

模拟结果应用分析：充分利用数值模拟结果，深入分析其在地球物理勘探应用中的问题，如对矿体分布、地下水和油气分布、地质构造等方面的影响，并积极探讨如何优化电磁场勘探方案，以提升勘探效率和准确性，为实际勘探工作提供有力的理论支持和技术指导。

本研究主要采用有限元方法对二维各向异性大地电磁场进行数值模拟。首先，建立包含二维各向异性的大地电磁