论文部分内容阅读
在磁法勘探的正演问题中，大多假设磁性体介质均匀，在均匀磁场中均匀磁化。这种简化的假设对了解各种磁性形体的磁场规律特征及反演解释有很火的作用。对大多数退磁作用可以忽略的非强磁性物体，只要其成分较均匀，这种假设基本符合实际情况。对强磁性物体，则需考虑退磁作用，如形体表面为二次曲面(三轴椭球体及其特例)，由于退磁场是均匀的，体内仍为均匀磁化，所以解析计算仍照样进行。若磁体表面为非二次曲面，则由于退磁作用的非均匀性，即使介质成分均匀，但退磁系数成为坐标的函数，磁性体内仍将形成较强的非均匀磁化，对磁异常的特征规律，必有影响；对于这种高磁化率矿体的磁异常正演却并不符合任何所有己知的近似，其中一个显著的困难就是需要解决磁场的边值问题。相比数值解来说，对这个问题的解析近似解只是按规则形体近似，从而进行退磁校正，对于磁异常的定量计算会产生非常大的误差，不利于解释工作。
为了解决这一问题，一般有三种可以使用的数值方法：积分方程法、有限单元法和边界单元法。积分方程法将边值问题转换成边界积分——第二类弗莱德霍姆(Fredholm)方程，但是在解方程的过程中会遇到“奇点”的难题；有限单元法则是将待求区域剖分，然后使用变分进行泛函求极值，利用理论物理学中的能量最小原理来进行求解；边界单元法的基本思想是将复杂形体异常的体积分通过高斯(Gauss)公式转化为面积分，再由格林(Green)公式转为线积分，而后累加求和得到整个形体的场值。
本文采用了边界单元法来计算这种高磁化率、非均匀磁化条件下的磁异常。相比于积分方程法，边界单元法不会产生奇点；相比于有限单元法在全区域的剖分，边界单元法只在研究区域的边界上进行单元的剖分。很明显，对比前两种方法，边界单元法要简单快捷的多。
通过边界单元法的计算，分别描述了均匀和非均匀磁化条件下的磁场特征，并和解析方法得到的结果进行了对比；同时，还研究了不同剩余磁化强度对磁异常的影响。另外，对于总磁场异常模量的特征也做了相应的说明。