复电阻率法是以岩、矿石电阻率的频谱或时谱特性为基础,通过测量大地的视复电阻率频谱或视时变电阻率谱,寻找电性异常体,并根据复电阻率参数来评价电性异常体.
Pelton等 (1978)发表的文章“利用多频电测量作矿物区分和去除电磁耦合”被认为是复电阻率法的基础文献.19世纪40年代初期,Cole和Cole (1941)提出利用公式来描述激发极化效应,公式为
(1) |
该数学模型被称为Cole-Cole模型.其中,ρ0为零频电阻率, m为极化率,τ为时间常数,c为频率相关系数,这些参数统称为Cole-Cole模型或复电阻率 (频) 谱参数.Pelton等 (1978)通过对大量岩、矿石露头和标本的广频谱复电阻率测量,证实了Cole-Cole模型的正确性,并指出可以根据岩、矿石的Cole-Cole模型参数,按结构评价极化异常体.
Soininen (1984, 1985) 利用积分方程法模拟了均匀半空间中极化直立长方体的视复电阻率谱.Xiong等 (1986),罗延钟等 (1987)在同时考虑电磁效应和激电效应的情况下,实现了复杂围岩中三维极化体的数值模拟,并就其正演相应特征进行了深入的讨论.孟永良等 (2000)利用有限单元法实现了2.5维时间域频谱激电正演,其采用的时间谱电阻率法能对可极化和导电大地 (即同时包括IP和EM效应) 的电磁场瞬变响应作计算.蔡军涛等 (2007)以复电位偏微分方程为基础,利用有限元法实现了不考虑电磁效应的二维正演模拟.岳安平等 (2009)研究了包含激电效应的一维CSAMT电磁场响应.王大勇等 (2010)进行了三维复电阻率模型电磁场正演模拟研究,得出三维复电阻率体的电磁响应是由累积电荷、电磁感应和激电效应共同产生的.范翠松 (2013)进行了基于有限元法的复电阻率正、反演研究及应用,提出了电场振幅和电场相位联合反演的方案,实现了2.5维复电阻率法全区反演.
但是,目前仍然缺少对电磁场正演结果的深入分析,对电阻率和相位的变化范围及其影响因素也不甚了解.本文利用3D积分方程法和3D有限差分法对三维正演视电阻率和视相位结果进行了深入分析,得出激电效应和激发极化效应对正演结果的影响规律.
1 均匀半空间正演模型分析 1.1 模型描述与装置参数图 1为均匀半空间模型.大地背景电阻率为100 Ω·m,时间常数为0,极化率为0,频率系数为0.采用的观测方式为偶极-偶极,令AB为发射偶极,MN为接收偶极,O、O′分别为AB、MN中心点,发射偶极AB和接收偶极MN长度选用100 m,共观测了41个点,观测点坐标从 (-400, 0, 0) 至 (400, 0, 0).
图 2中 (a)(c) 和 (e) 是视电阻率曲线,(b)、(d) 和 (f) 是视相位曲线.对比 (a)、(c) 和 (e) 中三条视电阻率曲线,在Cole-Cole模型中,当频率为0.01 Hz时,视复电阻率值为大地背景电阻率100 Ω·m,当频率增加到128 Hz时,最大视复电阻率值达到90 Ω·m,可见复电阻率值随频率的增加而单调减小.激电效应可以看作是电容效应,所以,随着频率不断增高,激电效应逐渐减弱,复电阻率值逐渐减小.从空间上来看 (如图 2e),随着收发距增加,视电阻率值从90 Ω·m减小到45 Ω·m左右,可见视电阻率值随着收发距增加而不断减弱.
对比 (b)、(d) 和 (f) 中三条视相位曲线,视相位值恒为负.在0.01 Hz时 (如图b),视相位值为负的十分之几个毫弧度;在1 Hz时 (如图d),视相位值为负的几个到几十个毫弧度;在128 Hz时 (如图e),视相位值为负的几百个毫弧度.在高频时 (如图f),随着收发距变大视相位绝对值先变大后变小.
2 异常体空间正演模型分析 2.1 模型描述与装置参数图 3为低阻异常体空间模型.大地背景电阻率为100 Ω·m,异常体电阻率为10 Ω·m,时间常数为100,极化率为0.3,频率系数0.3,异常体中心点坐标为 (0, 0, 80),异常体的长、宽、高分别为160 m、80 m和80 m.采用的观测方式依然为偶极-偶极,发射偶极AB和接收偶极MN长度同样选用100 m,测线布置、采集频率与图 1中的均匀半空间模型相同.
图 4中 (a)、(c) 和 (e) 是视电阻率曲线,(b)、(d) 和 (f) 是视相位曲线.由 (a)、(c) 和 (e) 三条视电阻率曲线可以看出,视电阻率曲线在低阻异常体顶部出现小于正常背景电阻率的极小值,这是由于低阻异常体吸引电流的作用,使其顶部电流密度小于均匀半空间时的电流密度,从而出现视电阻率值减小的情况.而在其两侧,因为电流密度大于均匀半空间时的电流密度,从而出现视电阻率值大于背景电阻率值的情况,所以在低阻异常体两侧出现大于正常背景电阻率的极大值.
由 (b)、(d) 和 (f) 三条视相位曲线可以看出,在0.01 Hz时 (如图 4b),视相位值接近于零,在异常体中心点上方,视相位值为负的十几个毫弧度;在1 Hz时 (如图 2d),视相位值在负的几个到几十个毫弧度范围,在异常体上方,视相位值由负的十几个毫弧度迅速变化到负的二十几个毫弧度;在128 Hz时 (如图 4f),视相位值在负的二百到五百个毫弧度范围之间,视相位值由负的四百个毫弧度迅速变化到负的二百五十个毫弧度.
3 结论 3.1由均匀半空间模型的视相位曲线可以看出,视相位值恒为负的.且频率由低到高变化时,视相位的绝对值随之变大.而且在f=128 Hz时,视相位绝对值先变大后变小,曲线呈现出非单调曲线形态,这说明视相位与收发距的关系是非线性的.
3.2由低阻异常体模型的视相位曲线可以看出,频率由低到高,电磁效应增强,激发极化效应减弱.低频时,激发极化效应远大于电磁效应,激发极化效应使正演视相位变大;高频时,电磁效应远大于激发极化效应,电磁效应使正演视相位变小.
3.3电磁效应与激发极化效应都会影响实际数据的处理结果,在传统的激发极化法和电磁感应法中,电磁效应与激发极化效应互相干扰,因此在实际工作中必须予以考虑,否则会影响最终的结果.
致谢 感谢审稿专家提出的修改意见和编辑部的大力支持![] | Cai J T, Ruan B Y, Zhao G Z, et al. 2007. Two-dimensional modeling of complex resistivity using finite element method[J]. Chinese J. Geophys., 50(6): 1860–1876. DOI:10.3321/j.issn:0001-5733.2007.06.030 |
[] | Cole K S, Cole R H. 1941. Dispersion and absorption in dielectrics I. Alternating current characteristics[J]. The Journal of Chemical Physics, 9(4): 341–351. DOI:10.1063/1.1750906 |
[] | Fan C S. 2013. Research on complex resistivity forward and inversion with finite element method and its application (in Chinese)[Ph. D. thesis]. Changchun:Jilin University, 69-99. |
[] | Luo Y Z, Xiong Z H, Cui X W. 1987. Properties of the spectral IP anomaly in coexistence of IP and electeomagnetic effects[J]. Geophysical and Geochemical Exploration, 11(2): 106–114. |
[] | Meng Y L, Luo Y Z, Chang Y J. 2000. 2-Dimensional forward algorithm for time spectral resistivity[J]. Earth Science-Journal of China University of Geosciences, 25(6): 656–662. |
[] | Pelton W H, Ward S H, Hallof P G, et al. 1978. Mineral discrimination and removal of inductive coupling with multifrequency IP[J]. Geophysics, 43(3): 588–609. DOI:10.1190/1.1440839 |
[] | Soininen H. 1984. The behavior of the apparent resistivity phase spectrum in the case of a polarizable prism in an unpolarizable half-space[J]. Geophysics, 49(9): 1534–1540. DOI:10.1190/1.1441778 |
[] | Soininen H. 1985. The behavior of the apparent resistivity phase spectrum in the case of two polarizable media[J]. Geophysics, 50(5): 810–819. DOI:10.1190/1.1441955 |
[] | Wang D Y, Li T L, Li J P, et al. 2010. Forward simulation of a 3D complex resistivity model[J]. Progress in Geophysics, 25(1): 266–271. DOI:10.3969/j.issn.1004-2903.2010.01.035 |
[] | Xiong Z H, Luo Y Z, Wang S T, et al. 1986. Induced-polarization and electromagnetic modeling of a three-dimensional body buried in a two-layer anisotropic earth[J]. Geophysics, 51(12): 2235–2246. DOI:10.1190/1.1442077 |
[] | Yue A P, Di Q Y, Wang M Y, et al. 2009. 1-D forward modeling of the CSAMT signal incorporating IP effect[J]. Chinese J. Geophys., 52(7): 1937–1946. DOI:10.3969/j.issn.0001-5733.2009.07.030 |
[] | 蔡军涛, 阮百尧, 赵国泽, 等. 2007. 复电阻率法二维有限元数值模拟[J]. 地球物理学报, 50(6): 1860–1876. DOI:10.3321/j.issn:0001-5733.2007.06.030 |
[] | 范翠松. 2013. 基于有限元法的复电阻率正反演研究及应用[博士论文]. 长春: 吉林大学, 69-99. |
[] | 罗延钟, 熊宗厚, 崔先文. 1987. 同时存在激电和电磁效应的频谱激电异常性质[J]. 物探与化探, 11(2): 106–114. |
[] | 孟永良, 罗延钟, 昌彦君. 2000. 时间谱电阻率法的二维正演算法[J]. 地球科学-中国地质大学学报, 25(6): 656–662. |
[] | 王大勇, 李桐林, 李建平, 等. 2010. 三维复电阻率模型电磁场正演模拟研究[J]. 地球物理学进展, 25(1): 266–271. DOI:10.3969/j.issn.1004-2903.2010.01.035 |
[] | 岳安平, 底青云, 王妙月, 等. 2009. 含激电效应的CSAMT一维正演研究[J]. 地球物理学报, 52(7): 1937–1946. DOI:10.3969/j.issn.0001-5733.2009.07.030 |