CSAMT二维正演研究

发布时间：2017-05-02 01:00

  本文关键词：CSAMT二维正演研究，由笔耕文化传播整理发布。

【摘要】： 目前国内大多数CSAMT资料处理都用一维反演方法，即假定地球介质是水平层状地电结构，这对于横向变化较大的介质结构的反演则是不真实的。为了合理的模拟可控信号源的三维特征应该使用三维有限元方法。受计算机内存和计算速度的限制，三维方法的实用化受到限制。在许多情况下，地电结构沿走向变化很小，只沿倾向发生变化。这种地电结构是二维的，而人工源是三维的，因此CSAMT资料的观测可用2.5维有限元方法进行数值模拟。本文从麦克斯韦方程组出发，建立了2.5维有限元CSAMT数值模拟方法，其核心是把地电参数变化小的走向方向通过傅立叶变换转化成波数域，用一系列波数模拟二维波数域电磁场的特征。用横向均匀的三层地电结构模型展示了2.5维数值模拟的特征，并和一维模拟结果进行了比较，证实了有限元方法2.5维CSAMT数值模拟的可靠性。在此基础上对地电结构已知的实际模型进行了2.5维数值模拟，计算其电磁场值，视电阻率和视阻抗相位，从而判断此2.5维数值模拟的正确性。为2.5维反演打下基础。 本文首先调研了CSAMT国内外的研究现状，确定研究的突破点和创新点在于通过反傅立叶变换，求出侧面剖面十公里范围内电磁场的变化特征，并与一维模型的解析解进行对比，了解二维模拟的精度。在第二章中，研究一维地电模型条件水平电偶源的理论公式和解析解，利用快速汉克变换滤波系数法，做了一维正演计算，并与前人相同地电模型下的解析解进行对比，确定了一维正演程序的可靠性和正确性，进而为验证二维有限单元算法做基础。第三章，具体讨论了沿构造方向地电结构不变或变化较小的条件下，在该方向做傅立叶变换，变换到波数域，通过麦克斯韦方程组及电磁场边界条件，建立电偶源情况下电磁场满足的偏微分方程组，，然后通过泛函建立与偏微分方程组等价的变分问题。第四章，详细介绍变分问题下有限单元数值计算方法，包括三角单元剖分，线性方程组的建立，下三角的存储和线性方程组的简化，电磁场其它分量的求取。由于在有限单元算法中采用了二次场算法，线性方程组的求解需要波数域一次场的电磁场值，所以第五章详细推导了波数域电磁场一次场的求解公式。求解到波数域电磁场场值后，需要做反傅立叶变换才能得到常用空间域的电磁场值，第六章则详细介绍了反傅立叶变换的理论和方法。为了检验有限单元算法的正确性，需要做一维地电条件下有限单元数值解和解析解的对比，在第七章，本文详细列出了H型、K型、KH型地电结构下有限单元数值解和解析解的对比图及误差图，对比的场值包括Ex、Hy和Hz，对比剖面包括1公里处、2公里处、4公里处、6公里处、8公里处、10公里处、12公里处；频率范围为10~(-2)～10~4Hz。为进一步验证有限单元二维正演的有效性，在第八章做了中间低阻体和中间高阻体两种二维模型的二维正演计算，绘出了这两种模型条件下的视电阻率图和视相位图。 由于可控音频大地电磁采用的是人工源，所以比起天然源有很高的信躁比，即使在有工业干扰的地区，也可以获得较好的信躁比。近年来，CSAMT方法在国内外打开了应用局面，并取得了显著的效果。该方法广泛应用在地下水资源勘察、地热勘察、工程物探、金属矿勘查、煤田地质灾害预测等领域。 本文利用有限单元算法，实现了电偶源二维正演程序。其创新点在于，在二维地电结构下，x方向布源，通过逆傅立叶变换算出侧面剖面12公里以内的电磁场特征值，并做了一维地电结构条件下有限单元算法结果和解析解的比较，从而验证该有限单元算法的准确性。另外，本文还做了两种异常体条件下的电磁场正演结果。 通过有限单元数值模拟的结果和解析解的对比，认为该正演结果在10公里以内可以达到较满意的精度。通过对异常体的正演模拟，得到了反映异常的效果，因此认为该正演程序有效，在低频情况下做适当改进方可用来做静校正，近场校正和二维反演。 不足之处在于，低频(频率小于10Hz)时，有限单元算法和解析解结果相差较大，需要做理论上及算法的探讨和研究。随着侧面剖面y值的增大，有限单元算法和解析解结果的误差随之增大，这是由于逆傅立叶变换波数选取的范围不够，以及在波数域中电磁场的变化曲线拟合程度不够。在靠近有限单元网格的边缘，数值模拟结果和解析解结果误差也较大，这就要求网格的剖分要合适，边界网格要够大。靠近源的测点电磁场值也算不准确，因为在有源的地方，电磁场趋近无穷大，电磁场边界条件不满足。
【关键词】：2.5D CSAMT 有限元法 数值模拟
【学位授予单位】：中国地质大学
【学位级别】：硕士
【学位授予年份】：2007
【分类号】：P631.3
【目录】：

• 摘要6-8
• ABSTRACT8-11
• 第一章 前言11-14
• 第二章 CSAMT一维正演14-18
• 2.1 水平层状大地上电偶极子的电磁场及其快速汉克变换算法14-15
• 2.2 一维模型对比15-18
• 第三章 水平电偶源有限元数值模拟基本原理18-27
• 3.1 变分问题的建立18-21
• 3.2 偏微分方程与变分问题的等价行证明21-27
• 第四章 有限单元数值计算方法27-54
• 4.1 单元剖分和线性方程组建立27-43
• 4.1.1 单元剖分27-35
• 4.1.2 线性方程组的建立35-42
• 4.1.3 直接求异常的线性方程组42-43
• 4.2 下三角刚度阵的求取及线性方程组的简化43-51
• 4.3 其它分量的求取51-54
• 第五章 水平电偶源波数域一次场的计算54-59
• 第六章 傅立叶逆变换求空间域电磁场值59-61
• 第七章 有限单元算法结果与解析解的对比61-71
• 7.1 H型断面对比图61-64
• 7.2 K型断面对比图64-67
• 7.3 KH型断而对比图67-71
• 第八章 二维地电模型数值模拟71-75
• 8.1 柱状低阻体模拟71-73
• 8.2 柱状高阻体模拟73-75
• 第九章 结论与建议75-76
• 致谢76-77
• 参考文献77-79

【引证文献】

中国期刊全文数据库 前3条

1 吴桂桔;胡祥云;刘慧;;CSAMT三维正演数值模拟研究进展[J];地球物理学进展;2010年05期

2 吴艳辉;王璞s

本文编号：339945