复杂条件下直流电阻率异常三维数值模拟研究.pdfVIP

21阢丌尧，村1．讲，徐世浙．I 999，}乜5I【半／激发极化法数据的一维反演程序：物探化抓汁算技术，2 l 6～I25 复杂条件下直流电阻率异常三维数值模拟研究 熊彬1，2阮百尧2罗延钟。 (1中国地质人学地球物理系武汉430074 2桂林I：学院资源与环境I．程系棒林541004) 1．引言 在山区进行电阻率法勘探时，复杂的地形条件是个重要的干扰因素。它往往会导致视 电阻率曲线发生严重的畸变，产生形似高、低阻体的假异常；有时甚至会掩没地质体的有 用异常；有时又歪曲有用异常，造成地质体异常形态的变更和位置的移动。而且，我国的 金属矿又多位于山区，因此对复杂条件下的电阻率异常规律进行研究具有重要的理论及实 际意义，是提高电阻率法地质效果的重要问题之一。 以往的文献中有不少用数值模拟的办法来研究复杂条件下的电场变化规律，但这其中 人都是对二维问题进行计算，或者是仅仅就纯地形的情况进行了讨论、分析，比如 研究是个难度很大的问题，上世纪八十年代以前，没有人涉及到这个问题。1984年，美国 地球物理学家Holcombe和Jirack试图用有限单元法来解决这一问题，然而需要求解含有 儿万个未知数的线性方程组，这只有在大型计算机上才能实现。同年，Oppliger用边界积 分方程研究这个问题，需解含1681个未知数的方程组。后来，徐世浙等(1985，1991) 开始用边界单元法研究这个问题，范国华等(1998)又采用有限差分的方法计算了三维地 形影响的电磁感应。 众所周知，有限单元法最适合于物性复杂分布的地球物理问题，而有限差分和边界元 法在这一领域又各有其局限性。本文将在前期层状介质数值模拟工作的基础上，给出对具 有复杂地形条件的地下异常体进行有限元数值模拟的具体办法。文中先通过对儿例典型的 角域纯地形问题进行研究，得到了与解析解拟合很好的视电阻率剖面及测深曲线。之后， 对复杂地形下赋存良导电矿体的情况进行分析，有限元数值计算的视电阻率曲线与土槽实 验结论一致。可见，随着计算机的迅速发展，用有限单元法求解三维电场问题是可行的、 高效的。 2．边值问题 三维构造中双点电源电场中电位的边值问题归纳如下： IV(oVv)=一218(A)+2：6(8)∈Q； {av／an=0∈t； (1) 【加／锄+l／(％一o)‘[‰／o ∈L cos(rA，)一rA／r．cos(rB，”)]v=0 式中·r为区域Q的地面边界，L为区域Q的地12／界，为边界的外法向方向 仃为介质的电导率，v为电位。 3．变分问题 与上述边值问题(1)等价的变分问题为： F(V)=f[I／20(VV)2—218(A)V+218(B)V】艘 n +l／2 (2) j叫(％一rA[rUr』eos(_，”)一_／％cos(‰，)】V2dF； r叩 8F(v)=0． 4．有限单元法 在前期工作的基础上，本文用四面体单元取代六面体单元，对研究区域进行剖分，在 每一网格单元内对电性参数(电导率、电位)进行三线性插值，较之以前更具针对性。在 野外实际情况下，地形往往是起伏不平的，而地形的影响又是不可忽视的，它直接影响到 异常的解释推断工作，而采用四面体单元对区域进行网格化正好解决了这一矛盾，原因是 它可以很好地拟合野外实际地形。 本研究由国家自然科学基金项目、广西自然科学基金匹配项目资助。 参考文献 P in modeled Wannamaker finite