二维极化介质的复电阻率正演的研究.pdfVIP

二维极化介质的复电阻率正演研究 蔡军涛阮百尧朱光萍 (桂林工学院资源与环境工程系桂林541004) 0～～ 激发极化的形成和衰减是～个比较缓慢的过程，在频率域中激电效应基本上只发生在超低频段上(，?·1 珂·10‘Hz)。对于这样缓慢变化的电场，通常可以忽略电磁感应和电磁辐射效应，因此复电阻率法的正演计算可 采用“等效电阻率法”，只要将地下各种地质体的真电阻率P，I／=1，2,3，…)换成相应的对给定频率厂’按柯尔一 柯尔模型计算的复电阻率 f 『 ， ]] ，一 ’¨’ pj(i∞)巩11--mj卜而寿”‘产l’2，…‘) 式中Po，、c，、f，和m，分别为第．，个地质体零频电阻率、频率相关系数、时间常数和充电率。则对无激电效 应的一次场电位作数值模拟的各种方法，皆可直接用来模拟计算体极化时频率域的总场电位值，进而计算给定电 极装置在该频率上的视复电阻率p。tia,J。 有一定走向的构造称为二维构造，但在复电阻率法中常用交流点电源供电，所以其产生的极化总场实质上是 三维的，可用傅立叶变换的方法，将这种极化场变换至波数域中进行求解。在波数域中，通过对无穷远边界的近 似处理，体极化二维构造中的交流点电源的极化总场电位的变分问题简化为： 即(训=』[·o-(2iao)f、V凇))2≯12巾)凇)2叫舻(伊(肋)卜 (2) 6F(U(ico))=O (2)式中U(fCOl为体极化总场的电位：II f∞)为交流点电源的电流强度；仃《fCO)为极化介质的复电导率，其 倒数为极化介质的复电阻率P(泐)，P(f∞}可由(1)式计算得到。 用有限单元法计算体极化总场电位的变分闯题。将整个二维区域划分为两个区域：目标区和边界区。目标区 为目标地质体赋存区域，按照需要以均匀网格剖分；边界区的网格单元步长呈指数递增，以模拟无穷远边界，边 界距离为整个区域的边界与目标区边界在X、Y方向上的距离。采用三角单元对积分区域(2进行剖分，在每一 个三角单元内对电性参数(电位【，If国)和复电导率盯If缈；)进行线性插值。计算出各个节点的波数域电位 U《ico}，采用傅立叶余弦逆变换求得三维空间的极化总电位“lf国l，然后对于给定的供电电极装置计算出视复 电阻率p。(j∞)。对不同测点的4个不同频率(4个频率分别为f=2。，厶=2¨‘，疋=2“2，无=2¨3， k=一3)的二维极化介质的视复电阻率p。lf国J，用阮百尧提出的递推反演方法计算4个柯尔一柯尔视谱参数 风一m。、r。和c一如此反复直至计算出所有测点的柯尔一柯尔视谱参数。 为了研究二维极化介质的复电阻率实、虚部和柯尔一柯尔视谱参数异常的形态特征，选取了几个二维极化模 型进行了数值模拟计算。采用偶极一偶极装置进行测量，偶极长度AB=I斛=4m，间隔系数n=l～8，供电电极移 164m。 通过对二维极化模型的数值模拟计算，发现二维极化介质相近频率的视复电阻率实、虚部断面图的异常形态 是基本类似的(只不过二维极化介质的视复电阻率虚部都为负值)。二维极化介质的4个柯尔一柯尔视谱参数与视 复电阻率的异常形态也比较相似，根据极化介质的柯尔一柯尔视谱参数断面图，可以大致判断极化介质的几何形态 和柯尔一柯尔真谱参数分布规律．但在利用柯尔一柯尔视谱参数断面图推断解释组合极化体时，必须注意其存在的 假异常的干扰。二维极化模型的视复电阻率和柯尔一柯尔视谱参数等值线断面图都有比较明显的异常反映，综合 两者可以对极化异常体的正确解释有很大的帮助。 本研究由广西留学回国人员基金项目(桂科回0448021)资助。 参考文献 …徐世浙．地球物理中的有限单元法．北京：科学出版社，1994． 12l阮百尧，罗润林．一种新的复电阻率频谱参数的递推反演方法．物探与化探计算技术，vol。25，NO．4，2003． · ·