激发极化法参数的有限元数值模拟.docVIP

激发极化法参数的有限元数值模拟 　　[摘要]在三维地电模型条件下，本文用有限元方法模拟了不同模型下频率域激电参数的响应。用不同频率下Cole-Cole模型所得到的复电阻率值替换大地模型中围岩和体极化异常体的电阻率，用有限单元法对大地模型进行六面体和四面体剖分，得到各个节点上的复电位，根据激电参数的定义得到中梯装置和偶极装置下的Fs和Фs，通过对比正演模型可知该方法是可行和合理的。 　　[关键词]激发极化法 有限单元法 幅频率Fs 相频率Фs 　　[中图分类号] P631.3+24 [文献码] B[文章编号] 1000-405X（2014）-5-225-3 　　1引言 　　激发极化法引入我国已经有将近半个世纪了，从最初的以观测交变总场电位差幅值为基础的“变频激电法”，发展到后来的“奇次谐波法”，和由何继善院士发明的“双频激电法”、“伪随机信号方案”[1]。对于激发极化法的数值模拟国内外已经有很多人做了研究，Coggon[2] 是电法勘探领域里最早发表论文阐述有限单元法进行正演模拟的人，他从电磁场总能量最小原理出发，实现了二维地电断面有限元法正演计算，我国地球物理学家李大潜是最早将有限单元法引入到我国的人，后来徐世浙[3] 、罗延钟[4] 、阮百尧[5] [7] 等国内地球物理学家对有限元法在地球物理中的发展进行了开拓性的工作。 　　目前激电法的数值模拟主要以“等效电阻率法”的时间域模拟和以计算Cole-Cole模型的频率域激电法模拟为主[6] 。频率域激电法中的激电参数很多，本文通过用有限单元方法来模拟其中的主要参数，对分析和减少地球物理中的多解性做一些有益的探索和尝试。 　　2数值模拟的基本原理 　　由于频率域激电法中所选用的频率一般在10-2～102 Hz之间，所以可以忽略位移电流的影响。为了简化计算在模拟计算的时候不考虑电磁效应（矢量势A不随时间变化），在体极化条件下，将三维模型中的大地和极化异常体的电阻率换成用Cole-Cole模型表示的复电阻率，通过模拟计算就可以得到不同频率下地表的复电位响应[8] 。三维构造条件下双点电源电场模型的电位边值问题为[9] 。 　　与上述边值问题（1）等价的变分问题是[3] 　　令泛函数方程 的变分为零，就得到一组线性方程组。 　　对三维地电模型通常的剖分方式是六面体单元剖分，但它的精度比较低。在不增加节点数又能尽量多地剖分单元，将一个六面体单元首先剖分成两个三棱柱单元，再对三棱柱单元进行右手螺旋规则的四面体剖分[11] （如图1，2所示）。 　　3频率域激电法中的参数 　　频率域激电法中的参数很多，有视幅频率 、相位差 和电位差的虚实分量等，参照幅频率的定义又出现了视相频率 ，这么多参数中只有一个量是独立的。视幅频率的表达式可表示为： 　　式中ΔVL、ΔVH分别表示供低频和高频电流时测量电极MN间的电位差。 　　相频率是描述相位相对于频率的变化关系，或其频散特性，借鉴频散率与幅频率的定义来描述激电相位与频率的关系，相频率的计算公式为： 　　其中φL、φH分别为低频和高频时的激电绝对相位。 　　4正演模拟 　　4.1模型1 　　从图4看，当异常体缓慢向下移动时，在地表得到的幅频率极大值也会缓慢变小，这符合实际情况，当改变频率对时，幅频率的变化曲线没有大的改变，说明频率对异常情况的影响不大（图5、6），时间常数是区分矿物属性的重要参数，图7中可以看出不同时间常数下幅频率曲线的变化情况。 　　4.2模型2 　　根据频率域激电参数幅频率的定义形式，可以利用相位参数对正演模型进行研究，假设有一个直立硫化物矿体，它的时间常数 ，极化率 ，其它参数与模型1相似。可以得到偶极装置下频率对为2，2/13Hz时的相频率断面图8。 　　从图8中可以看出在断面图中心位置存在相互对称的两个相频率极大值区域，根据正演模型可知它们的对称中心位置是异常体的位置。在区分相同充电率和电阻率，不同时间常数的体极化异常体时，相频率参数可以作为区分它们的重要参数。 　　5结论 　　（1）根据模拟结果，可以看出模拟得到的视幅频率能够清楚地反映异常体的位置和极化强度情况，以及幅频率随充电率、时间常数的变化情况，与实际情况是符合的，也说明频率域激电和时间域激电在反映激电异常的能力上是等效的。在电性界面复杂的情况下，为用有限单元法进行正演模拟提供了很好的应用前景。 　　（2）用有限单元法来模拟频率域中的激电参数，方法简单、过程明晰，适合于电性复杂分布的地电模型，此外，相频率参数Фs作为区分不同时间常数异常体的做法也是值得关注和研究的，为减少地球物理中的多解性做一些有益的探索和尝试。 　　The finite element numerical simulation for indu