雷东记 周猛 孟慧 赵晨光

摘要：实验室煤体复电阻率测试旨在提供可靠的复电阻率数据，为进一步现场应用基础研究。影响煤体复电阻率测量精度的因素众多，其中电极板为最重要的因素之一，极板的厚度和大小都会对煤体复电阻率测量产生影响，因此电极板的变化对测量参数的准确性有重要意义。通过选用导电性较好的紫铜材料作为极板，对3种不同厚度的极板材料进行了复电阻率测量，分析极板厚度对激发极化效应的影响，测试了直径为1，2，3，4，5 cm 5种大小极板对煤体复电阻测量的影响，并使用Cole-Cole模型进行数据拟合。得出结论如下：1）极板本身的极化效应和极板与煤样端面贴合程度是影响煤样复电阻测量误差的主要因素，3种不同厚度的极板材料中导电铜纸误差最小;2）对同一个煤样而言极板面积变化仅对复电阻率幅值有较大影响，对煤样激发极化的影响可以忽略。

关键词：复电阻率;极板材料;Cole-Cole模型

中图分类号：TD12

文献标志码：A文章编号：1000-582X（2022）02-103-11

Abstract：   The purpose of complex resistivity measurement of coal in laboratory is to provide reliable data of complex resistivity for further field application. The measurement accuracy of coal complex resistivity is affected by many factors， and electrode plate is one of the most important factors. The thickness and the size of electrode plate will affect the result of coal complex resistivity measurement， so the change of electrode plate is of great significance to the accuracy of measurement parameters. In this paper， the complex resistivity of the electrode material with three different thicknesses is measured by using the copper material with better conductivity as the electrode plate. The influence of the thickness of electrode plate on the induced polarization effect is analyzed. The influence of different sizes of electrode plates（diameter of 1 cm， 2 cm， 3 cm， 4 cm and 5 cm） on the measurement of the complex resistivity of coal is tested， and the Cole-Cole model is used for data fitting. The results show that 1） the polarization effect of the electrode plate itself and the bonding degree between the electrode plate and the end face of the coal sample are the main factors affecting the measurement error of the complex resistivity of the coal sample， and the error of the conductive copper paper is the smallest among the three electrode materials. 2） For the same coal sample， the change of the area of the electrode plate only has a great influence on the amplitude of the complex resistivity， and the influence of the area on the induced polarization of the coal sample can be ignored.

Keywords： complex resistivity; electrode material; Cole-Cole model

復电阻率法（CR， complex resistivity）作为一种新兴的电法勘探技术，其原理为利用岩石在外电场激励下产生电压响应滞后的激发极化效应来达到地球物理勘探的目的[1]。复电阻率法最早应用于金属矿产资源勘探[2]和地下水资源探测[3]，并取得了较好的效果。1978年，Pelton等[4]将Cole-Cole模型应用到激发极化理论中，极大推动了复电阻率法的发展。80年代以来，中国开展了复电阻率法在油气勘探中的理论研究与应用[5-7]，其测井技术在区分含油水层、评价岩石湿润性和渗透性等方面进行了大量研究[8-9]，何继善[10]开发了双频激电法，将装备笨重的时间域激电仪改为便携性高的频率域激电仪，并获得更好的测量效果。如今，复电阻率法应用广泛，其“非侵入性”的特点对于地下特征研究具有重要价值，例如土壤污染物监测[11]，地下水修复监测[12]，以及细菌活性监测[13]等，展现了复电阻率法广阔的应用前景。

煤體复电阻率实验是研究煤体激发极化现象与理论的基础，也是复电阻率法测井技术应用的前提，但目前仅有少量的相关研究见于报道[14-15]。贾将等[16]设计了一套复电阻率测量装置，提高了测量信号信噪比，可实现时域和频域测量，推动了煤岩复电阻率实验的发展。魏光华等[17]将电极测量结构改为四电极，相较于传统的二电极中极化电压的影响，大大提高了测量精度，但却忽略了电极板本身的极化影响;杨迪等[18]为避免电极板的极化效应，在样品和测量电极之间使用浸泡过饱和硫酸铜溶液的海绵隔开，Khajehnouri等[19]通过复电阻率法测量混凝土样件的内部损伤，利用导电凝胶将铜极板与样品结合，达到了较高的测量精度。关于极板本身的极化和极板大小对实验测量结果的影响不可忽略，但目前还鲜有报道。

笔者选用导电性较好的紫铜极板，分别测量3种厚度极板的复电性参数，研究极板本身极化的影响;通过测量不同接触面积的煤样复电阻率实验，研究接触面积对复电阻率测量的影响，相关研究对煤体复电阻率实验和测井技术发展具有重要意义。

1 实验系统与样品

1.1 实验测量系统

实验采用四极法（图1所示）测量煤岩复电阻率，克服了二级法中电极本身对煤岩样品复电阻参数的误差影响，四极法是将测量电极和供电电极分开，其中AB两个电极为供电电极，MN两个电极为测量电极，两个电极互不影响。选用的测量仪为IM3533-01LCR四端子阻抗测试仪，如图2所示，其中HCUR端子为电流发生端子，HPOT端子为高侧电压检测端子，LCUR端子为电流检测端子，LPOT端子为低侧电压检测端子，四端子结构可提高测量精度。该阻抗测试仪具有速度快、范围广、精度高等优点。实验选取的测量参数为复电阻实部R和虚部X，测量频段为0.01～100 kHz。因在0.01～0.1 kHz频率，由于频率过低，容易受到外界干扰造成测量结果波动较大。当频率大于0.1 Hz，测量结果趋于稳定，在10～100 kHz有明显频散特征[20]。由于电感作用极板开始出现明显的频散现象，为了保证测量精度的要求，在每次测量前需对仪器进行30 min以上的预热，并且减少外界噪声的干扰。

1.2 测试流程

首先将仪器放置在干净隔音的房屋中，开机预热1 h以上，连接电脑，选择测量参数为Rs和X，测量频率范围为0.1～100 kHz，选择测量速度为SLOW，测量速度越慢，测量准确度越高，选择平均值为3，即每个点的测量次数为3次，取平均值。设置完参数之后要进行测量补偿，分别进行线路补偿、开路补偿和短路补偿，其中外界噪声对补偿效果影响很大，要在静音环境下进行。补偿成功后，将煤柱放置在极板中间，使极板与煤柱接触良好，开始测量，连续测量3次，测量结束后保存实验数据并关机。

1.3 样品制备

本实验选取了来自焦作矿区的无烟煤和平顶山矿区的肥煤进行分析，煤样的煤质特征如表1所示。根据标准试件尺寸要求，将煤加工成50 mm×100 mm煤柱，为减少端面效应，上下两端面打磨光滑，要求平整度不大于0.02%，加工好的样品如图3所示。

2 实验测量与结果

2.1 极板厚度复电阻测量实验

极板选取为导电较好的紫铜，其中极板厚度为2.00，6.00，0.06 mm，其中0.06 mm为铜箔纸。将极板材料分别加工成60 mm×60 mm正方形，以至于完全覆盖煤样截面。导电极板实物图如图4所示。

实验对3种厚度的铜极板进行复电阻率实部和虚部测量，测量结果如图5所示。在0.1～10 kHz以内，3种材料实部结果相似，皆在0 Ω·m附近，而在10 kHz以后，厚极板实部和虚部数值发生明显的频散现象，造成此现象的原因是厚极板随着频率的增高，即电流发生迅速变化时，就会产生电磁辐射现象，导致明显的电磁干扰现象，对实验造成很大的误差，而导电铜纸和薄铜极板的数值以及走势规律相似。为进一步选择，选取实验煤样分别用导电铜纸和薄极板进行复电阻测量。测量结果如图6所示，无烟煤测量结果中薄极板的实部和虚部较导电铜纸结果规律性较差，尤其是虚部在频率较低时，出现逐渐下降阶段，与激发极化理论不相符，这是薄极板自身的极化效应以及极板与煤样间无法紧密贴合，留有空隙造成的端面效应共同导致的，属于实验测量误差，而导电铜纸可以与煤样端面进行良好的贴合;肥煤测量结果中导电纸所测的实部和虚部具有明显的频散特征，虚部曲线极值点，而薄电极板测量的虚部没有拐点，不利于进一步的分析。所以综合考虑，选择导电铜纸作为极板材料。

2.2 极板面积复电阻测量实验

在实验室测量时，通常将极板完全覆盖煤样端面，而现场测量是不可能实现的，所以为了更贴近实际应用情况，研究不同面积的导电纸对复电阻率测量数据的影响，实验采用5种1，2，3，4，5 cm不同直径的导电铜纸作为电极材料，如图7所示，并对无烟煤进行复电阻率测量，进行对比分析。

测量结果如图8所示，随着测量面积的增加，煤样的复电阻率实部和虚部都呈现逐渐减少的趋势，且不同直径的实部拐点和虚部的极值点大致在1 000 Hz频率段中。

3 分析与讨论

3.1 复电阻率频散机理

由基本电磁理论可知，交变电磁场中存在2种电流，一是传导电流，二是位移电流。传导电流是由自由电子或者其他带电粒子的定向运动所产生;而位移电流为极化分子定向排列（介电极化）所致，两者间相位相差π/2。煤体复电阻率频散特性的微观机理是煤岩孔隙中大量带正电荷的离子游离于溶液中，煤体介质中存在大量带负电荷的成分，且因正负电荷间存在引力，故形成导体溶液双电层平衡结构，交变电流场引起煤体介质中电荷的分化和转移，在外电流的激励下，正负电荷发生电化学反应开始重新分布，致使孔隙中离子双电层形状发生改变重新达到新的动态平衡，形成激发极化效应。

學者们为描述岩石复电阻率频散特征，依据等效电路法建立了各种模型，如CPA模型[21]等。目前描述激发极化的模型中应用最广的是Cole-Cole模型，该图9中（a）是岩（矿）石激发极化单元的Co1e-Cole模型，图9（b）是它的等效电路。

3.2 Cole-Cole模型拟合和数据分析

无烟煤使用导电铜纸测量时虚部为单极值曲线，使用薄极板测量时虚部为双极值曲线，在之前的研究中[20]笔者发现Cole-Cole模型对双极值曲线拟合效果较差，这里不再研究无烟煤。运用上述Cole-Cole模型推导出各分量公式对肥煤数据（实部电阻率和虚部电阻率）进行拟合。如图10所示，拟合得到模型参数如表2所示。

从图9和表2可以看出Cole-Cole模型对导电铜纸和薄极板测量的肥煤虚部数值拟合效果都很好，其R2达到0.99。但从拟合参数的数值来看，除c值变化不大，其余参数都有较大的变化，其中薄极板所测的m值为10.28，与m值的取值范围0～1有较大偏差，而导电纸所测的m值为1，在m值取值范围内。所以，通过模型拟合参数也可证明导电纸测量的数据具有较好的效果和可信度，与前文实验结果相印证。

为研究不同接触面积对煤样复电阻率结果的影响，采用Cole-Cole模型对实验测量的复电阻率虚部进行拟合，拟合曲线如图11所示。

4 结 论

通过测量3种极板的复电性参数，研究极板本身极化的影响，并开展了不同接触面积的煤样复电阻率实验研究，得出结论如下：

1）极板本身的极化效应和极板与煤样端面贴合程度是影响煤样复电阻测量误差的主要因素，3种不同厚度的极板中测量结果：厚极板的频散曲线干扰严重，而薄极板谱和导电纸谱彼此接近，更为稳定。对于同一样品，在不同厚度电极板的条件下采集的复电阻率频谱完全不同，厚极板为电极板进行测量会产生误差，不能用于数据解释;薄极板无法与煤样端面紧密贴合，且自身存在极化效应;以导电纸为电极板的测量结果与激发极化现象更为吻合，且经模型参数验证其误差也是最小，即以导电纸为电极板进行煤体复电阻率测量可提供更为可靠的复电阻率频谱。

2）同一个煤样，极板面积变化仅对复电阻率幅值有较大影响：即随着测量面积的增加，煤样的复电阻率实部和虚部绝对值都呈现逐渐减少的趋势，且不同直径的实部拐点和虚部的极值点大致在1 000 Hz频率段中。依据复电阻率频散机理和模型参数对比可知，参数ρ0随着煤样的有效导电面积增加而减小，极板面积变化对参数m、c、τ的影响较小，即对煤样激发极化的影响较小，在实际测量中可忽略。

相关研究为提高实验室复电阻率测量精度提供依据，为进一步复电阻率法的现场应用打下了基础。

参考文献：

[1] Vinegar H J， Waxman M H， Best M H， et al. Induced polarization logging： borehole modeling， tool design and field tests[C]∥SPWLA 26th Annual Logging Symposium， 17-20 June， Dallas， Texas. Society of Petrophysicists & Well Log Analysts， 1985： 62.

[2] 何继善， 鲍光淑. 频率域激发极化法的观测系统[J]. 中南矿冶学院学报， 1986， 17（4）： 1-9.

He J S， Bao G S. On induced polarization surveying system[J]. Journal of Central-South Institute of Mining and Metallurgy， 1986， 17（4）： 1-9. （in Chinese）

[3] Chongo M， Christiansen A V， Fiandaca G， et al. Mapping localised freshwater anomalies in the brackish paleo-lake sediments of the Machile-Zambezi Basin with transient electromagnetic sounding， geoelectrical imaging and induced polarisation[J]. Journal of Applied Geophysics， 2015， 123： 81-92.

[4] Pelton W H， Ward S H， Hallof P G， et al. Mineral discrimination and removal of inductive coupling with multifrequency IP[J]. Geophysics， 1978， 43（3）： 588-609.

[5] 罗延钟， 吴之训. 谱激电法中频率相关系数的应用[J]. 地球物理学报， 1992， 35（4）： 490-500.

Luo Y Z， Wu Z X. The application of frequency dependent factor in spectral induced polarization method[J]. Chinese Journal of Geophysics， 1992， 35（4）： 490-500. （in Chinese）

[6] 柯式镇， 冯启宁， 何亿成， 等. 电极法复电阻率测井研究[J]. 石油学报， 2006， 27（2）： 89-92.

Ke S Z， Feng Q N， He Y C， et al. Study on complex resistivity well logging with electrode antenna[J]. Acta Petrolei Sinica， 2006， 27（2）： 89-92. （in Chinese）

[7] 柯式镇. 井壁电成像测井全三维数值模拟与裂缝评价模型[J]. 中国科学（D辑：地球科学）， 2008， 38（S1）： 150-153.

Ke S Z. 3D numerical simulation and fracture evaluation model of borehole electrical Imaging logging[J]. Science in China（Series D： Earth Sciences）， 2008， 38（S1）： 150-153. （in Chinese）

[8] 肖占山， 曾志国， 朱世和， 等. 基于岩石电性参数频散特性评价润湿性的实验方法研究[J]. 地球物理学报， 2009， 52（5）： 1326-1332.

Xiao Z S， Zeng Z G， Zhu S H， et al. An experimental study of wettability evaluation based on frequency dispersion property of rock electric parameters[J]. Chinese Journal of Geophysics， 2009， 52（5）： 1326-1332.（in Chinese）

[9] Tong M S， Tao H G. Permeability estimating from complex resistivity measurement of shaly sand reservoir[J]. Geophysical Journal International， 2008， 173（2）： 733-739.

[10] 何继善. 双频激电法[M]. 北京： 高等教育出版社， 2006.

He J S. Dual frequency IP method[M]. Beijing： Higher Education Press， 2006. （in Chinese）

[11] Deng Y P， Shi X Q， Revil A， et al. Complex conductivity of oil-contaminated clayey soils[J]. Journal of Hydrology， 2018， 561： 930-942.

[12] Orozco A F， Velimirovic M， Tosco T， et al. Monitoring the injection of microscale zerovalent iron particles for groundwater remediation by means of complex electrical conductivity imaging[J]. Environmental Science & Technology， 2015， 49（9）： 5593-5600.

[13] Mellage A， Smeaton C M， Furman A， et al. Linking spectral induced polarization （SIP） and subsurface microbial processes： results from sand column incubation experiments[J]. Environmental Science & Technology， 2018， 52（4）： 2081-2090.

[14] Volkmann J， Klitzsch N. Wideband impedance spectroscopy from 1 mHz to 10 MHz by combination of four- and two-electrode methods[J]. Journal of Applied Geophysics， 2015， 114：191-201.

[15] Abdulsamad F， Florsch N， Schmutz M， et al. Assessing the high frequency behavior of non-polarizable electrodes for spectral induced polarization measurements[J]. Journal of Applied Geophysics， 2016， 135： 449-455.

[16] 賈将， 柯式镇， 张冰， 等. 岩心激发极化实验系统研究[J]. 中南大学学报（自然科学版）， 2019， 50（4）： 908-914.

Jia J， Ke S Z， Zhang B， et al. Research on core induced polarization experimental system[J]. Journal of Central South University（Science and Technology）， 2019，50（4）： 908-914. （in Chinese）

[17] 魏光华， 赵学亮， 李康， 等. 基于双极性电压激励信号的四电极电导率测量仪的设计[J]. 仪表技术与传感器， 2019（5）： 34-37， 49.

Wei G H， Zhao X L， Li K， et al. Design of four-electrode conductivity instrument based on bipolar voltage excitation[J]. Instrument Technique and Sensor， 2019（5）： 34-37， 49. （in Chinese）

[18] 杨迪. 天然岩矿石复电阻率测量及频谱曲线特征研究[D]. 北京： 中国地质大学（北京）， 2019.

Yang D. Research on complex resistivity measurement and frequency spectrum curve of nature rock[D]. Beijing： China University of Geosciences （Beijing）， 2019. （in Chinese）

[19] Khajehnouri Y， Chouteau M， Rivard P， et al. Measuring electrical properties of mortar and concrete samples using the spectral induced polarization method： laboratory set-up[J]. Construction and Building Materials， 2019， 210： 1-12.

[20] Lei D J， Zhao C G， Zhang Y G， et al. Experimental study on characteristics of water containing coal complex resistivity dispersion[J]. Journal of Loss Prevention in the Process Industries， 2018， 56： 328-333.

[21] 杨振威， 许江涛， 赵秋芳， 等. 复电阻率法（CR）发展现状与评述[J]. 地球物理学进展， 2015， 30（2）： 899-904.

Yang Z W， Xu J T， Zhao Q F， et al. Current situation and review of Complex Resistivity[J]. Progress in Geophysics， 2015， 30（2）： 899-904. （in Chinese）

（编 辑 郑 洁）