论文部分内容阅读
由于中深度工程地球物理勘探的迫切需要,以美国EH-4电导率成像系统为代表的高频率大地电磁法在我国地球物理勘探行业应用越来越广泛。高频大地电磁法属于采集天然场信号的被动源电磁方法,采集的信号频率范围为10Hz~100KHz,研究深度从地下的十几米至上千米。论文研究的基于矢量有限元的高频率的被动源电磁测深法,紧密结合当前我国中深度地球物理勘探中电磁法大量应用的实际情况,不仅具有学术意义,对指导工程实践也有一定的价值。论文研究的重点是以标量有限元和矢量有限元为核心的、以高频率电磁波为对象的正演数值模拟方法。主要包括四部分:建立与高频大地电磁法边值问题相等价的变分问题,为有限元数值计算奠定数理基础;标量/矢量有限元方法的研究;加速标量/矢量有限元矩阵系统方程迭代求解技术的研究;实现高频大地电磁参数精确快速正演模拟的研究。希望给后续的正反演研究工作提供科学的依据和参考,在减少地球物理多解性方面做一些有益的探索和实践。本文的主要内容和创新成果如下:1.在前人工作的基础上,从麦克斯韦电磁场双旋度方程出发,分别利用广义变分原理和加权余量法推导了高频大地电磁场在有耗介质中边值问题的稳定泛函,建立了变分方程,为有限元的计算的稳定性和准确性提供了数理保证。2.传统的标量有限元在解决矢量电磁场边值问题时,需要将未知量转化为标量场问题,然后进行求解。这种基于标量基函数处理矢量电磁场问题时,会造成非物理解或伪解问题的出现、异常体表面强加边界条件的不方便以及处理介质或导体边缘及角的困难性。本文利用了一种新型的矢量插值基函数来近似未知函数,将自由度赋给单元的棱边而不是节点,避免了传统标量有限元的困难,所做的工作证明,矢量有限元在高频大地电磁数值模拟中应用的非常成功,得到了较好的结果。3.论文研究了一种新型的大型病态线性方程组的求解方法——改进的威尔金森方法。应用这种方法进行求解不仅加速了有限元线性方程组的迭代收敛速度,而且大大降低了矩阵方程的求解时间,特别是在求解大地电磁场矩阵方程的情况下,此时矩阵方程的性态极差,一般的数值求解方法很难收敛到比较满意的截断门限,但该方法依然能够取得较好的效果。4.以典型的模型,计算并详细分析了三维情况下,标量有限元和矢量有限元高频大地电磁的响应特征、精度和速度,分别探讨了山谷地形、山脊地形和复杂地形条件下电磁异常的特点和变化规律。在本文的结论部分还指出了一些不足和今后工作的建议。