文章快速检索     高级检索
  地质与资源 2020, Vol. 29 Issue (1): 76-84  
0

引用本文
商木元. 岩石电容——电介质电容之一[J]. 地质与资源, 2020, 29(1): 76-84.  
SHANG Mu-yuan. ROCK CAPACITANCE: ONE OF THE DIELECTRIC CAPACITANCE[J]. Geology and Resources, 2020, 29(1): 76-84.  

岩石电容——电介质电容之一
商木元     
天津地质研究院, 天津 300061
摘要:岩石是电介质之一.发现了岩石具有电容性,就应该探讨电介质是否具有电容性.本文根据岩石充放电一次电压和二次电压得到的岩石电容性给出了电介质电性模型.依据介质为电容电阻串联的新的模型,利用跨学科的前人对黄铜矿、黄铁矿、玻璃、云母、琥珀以及聚乙烯胶囊的时间域和频率域实验数据进行验证,并进一步与介电常数经典理论对比验证.验证结果表明,新模型的电性参数关系与黄铜矿、黄铁矿、玻璃、云母、琥珀及聚乙烯胶囊实测电性参数关系非常吻合;新模型适用于各种时间的弛豫以及直流和交流环境;新模型也完全能表达各种介电常数和电损耗的关系.主要结论有,介质是漏电的含有内阻的电容;它可以看作若干含内阻电容和漏电电阻的并联电路;介质的阻抗及等效电容取决于介质的电容、电容内电阻、介质电阻及频率;岩石电容是电介质电容的一种.
关键词岩石电容    内阻    电介质    弛豫时间    吸收电流    
中图分类号:P589.1            文献标志码:A             文章编号:1671-1947(2020)01-0076-09
ROCK CAPACITANCE: ONE OF THE DIELECTRIC CAPACITANCE
SHANG Mu-yuan     
Tianjin Institute of Geology, Tianjin 300061, China
Abstract: Rock is a kind of dielectric. Since the rock is found to be capacitive, it should be discussed about the capacitive characters of dielectric. The electrical model of dielectric is established based on the rock capacitance obtained by primary voltage and secondary voltage during charging and discharging. According to the new model with capacitor and resistor as the dielectric, verification is carried out using the time domain and frequency domain experimental data of chalcopyrite, pyrite, glass, mica, amber and polyethylene capsules by interdisciplinary predecessors. It is further verified by comparison with the classical theory of permittivity. The results show that the electrical parameters of the new model are perfectly consistent with those of chalcopyrite, pyrite, glass, mica, amber and polyethylene capsules. The new model is suitable for all kinds of relaxation time, direct current and alternating current environments, and can also fully express the relations between various permittivity and electrical loss. The main conclusions are as follows:Dielectric is a leaky internal resistance-containing capacitor; It can be regarded as a number of parallel circuits containing internal resistance capacitor and leakage resistor; The impedance and equivalent capacitance of dielectric depend on the capacitance, internal resistance, dielectric resistance and frequency; Rock capacitance is one of the dielectric capacitance.
Key words: rock capacitance    internal resistance    dielectric    relaxation time    absorption current    

0 引言

20世纪70年代,在地球物理勘探领域,Pelton给出了岩石阻抗的经验公式[1]. 21世纪初,有人给出了岩石具有电容的结论[2].在电介质领域,针对细胞膜相电导率研究,有人提出了高频和低频下的阻抗公式[3]. Debye给出了介电常数和弛豫时间及频率关系的经验公式[4],根据该公式可以得到介电常数和损耗的关系为一个半圆. Cole等对该公式做了一些发展,增加了一个圆弧[5-6]. Davidson、Havriliak等进一步发展,增加了一个斜歪的圆弧[7].现代对非均匀体的研究表明,介电常数和损耗的关系可以为多个圆弧.上述关系都有相应实验数据的支持[8-10].电介质研究表明其电性和岩石电性完全一致.

1 模型的建立

地球物理勘探发现,在向大地和岩矿标本直流充电和放电时有一次和二次电压现象,有人据此得到了岩石具有电容的结论[2];交流充电有相位移和阻抗随频率变化现象[11-13].前者是电容电阻串联电路的表现,后者也是电容的表现形式.而岩矿能够自然放电表明岩矿是导电的.岩矿也是一种介质.因此,可以用上述岩矿性质表达介质性质.考虑到介质组成的不均一性,介质电容电路等效于若干串联电容电阻的并联电路.介质的导纳Y和角频率ωω=2πf)的关系可表示为

$ Y=j \omega \frac{C_{1}}{1+j R_{1} C_{1} \omega}+j \omega \frac{C_{2}}{1+j R_{2} C_{2} \omega}+\cdots+j \omega \frac{C_{n}}{1+j R_{n} C_{n} \omega}+\frac{1}{R_{0}} $ (1)

式中$\dot j $为单位向量,ω为频率,C1C2Cn为介质内的并联电容,R1R2Rn为介质内相应电容的内阻(电容的串联电阻,电容内电导的倒数),R0为介质的漏电电阻(介质电导的倒数).

式(1)还可以写为

$ \begin{array}{l} Y=j \omega\left(\frac{C_{1}}{1+\left(R_{1} C_{1} \omega\right)^{2}}+\frac{C_{2}}{1+\left(R_{2} C_{2} \omega\right)^{2}}+\cdots+\frac{C_{n}}{1+\left(R_{n} C_{n} \omega\right)^{2}}\right) \\ +\frac{R_{1} C_{1}^{2} \omega^{2}}{1+\left(R_{1} C_{1} \omega\right)^{2}}+\frac{R_{2} C_{2}^{2} \omega^{2}}{1+\left(R_{2} C_{2} \omega\right)^{2}}+\cdots+\frac{R_{n} C_{n}^{2} \omega^{2}}{1+\left(R_{n} C_{n} \omega\right)^{2}}+\frac{1}{R_{0}} \end{array} $ (2)

式(2)中虚部和实部分别为容纳和电导,容纳是等效电容Ce和频率ω的乘积.容纳和电导比值为电流电压相位差的正切[14-15],也是介质电损耗角的余切.电导的倒数为介质的阻抗.

一定电压U下,并联的串联电容电阻的充放电I和时间t的关系可写为(正负符号分别充电和放电)

$ I=\pm U\left(\frac{1}{R_{1}}\left(\mathrm{e}^{\frac{-t}{R_{1} C_{1}}}\right)+\frac{1}{R_{2}}\left(\mathrm{e}^{\frac{-t}{R_{2}-C_{2}}}\right)+\cdots+\frac{1}{R_{n}}\left(\mathrm{e}^{\frac{-t}{R_{n} C_{n}}}\right)\right) $ (3)

含介质平行板电容短接充放电相当于平行板为导线接点的介质充放电.真空电容短接充放电是瞬时完成的,没有弛豫时间,其电容不随频率变化.当在平行板电容内加入介质后,电容增加,产生弛豫时间,电容随频率变化.充电过程中有吸收电流,短接弛豫放电电流和吸收电流镜像对应.平行板电容短接充放电瞬间完成,因此,可以认为产生弛豫的充放电是介质的缘故.含介质的平行板电容短接能够储存和释放电量并具有弛豫时间,表明介质具有电容.电容电阻串联充放电有一时间常数,这和弛豫时间相对应.这表明介质电容存在内阻.平行板电容加入介质后等效于平板电容并联上若干串联电容电阻.平行板电容不随频率变化.加入介质后电容的导纳Y可写为平行板电容容纳与介质电容导纳的和,有

$ \begin{aligned} Y &=j \omega C_{\mathrm{P}}+j \omega \frac{C_{1}}{1+j R_{1} C_{1} \omega}+j \omega \frac{C_{2}}{1+j R_{2} C_{2} \omega}+\cdots+\\ & j \omega \frac{C_{n}}{1+j R_{n} C_{n} \omega}+\frac{1}{R_{0}} \end{aligned} $ (4)

CP为平行板电容为介质影响电场后的电容,其值为含介质电容频率趋于无穷的电容值.这可以理解为频率趋于无穷时,介质含内电阻电容的等效电容趋于0,电路的电容就是平行板电容CP.式(4)可以还写为

$ \begin{aligned} Y=& j \omega\left(C_{\mathrm{P}}+\frac{C_{1}}{1+\left(R_{1} C_{1} \omega\right)^{2}}+\frac{C_{2}}{1+\left(R_{2} C_{2} \omega\right)^{2}}+\cdots+\right.\\ &\left.\frac{C_{n}}{1+\left(R_{n} C_{n} \omega\right)^{2}}\right)+\frac{R_{1} C_{1}^{2} \omega^{2}}{1+\left(R_{1} C_{1} \omega\right)^{2}}+\frac{R_{2} C_{2}^{2} \omega^{2}}{1+\left(R_{2} C_{2} \omega\right)^{2}} \\ &+\cdots+\frac{R_{n} C_{n}^{2} \omega^{2}}{1+\left(R_{n} C_{n} \omega\right)^{2}}+\frac{1}{R_{0}} \end{aligned} $ (5)

式(5)的虚部为容纳,实部为电导.在单位电压下,无功功率(Q)数值上等于容纳,有功功率(P)数值上等于电导.电导和频率ω的比值为单位周期电损耗/2π,取名为周期电耗Cε.容纳和电导的比值为电流电压位移角的正切,也是损耗角的余切.容纳是等效电容Ce和频率ω的乘积.式5还可以写成

$ \begin{aligned} Y=& j \omega\left[\left(C_{\mathrm{P}}+\frac{C_{1}}{1+\left(R_{1} C_{1} \omega\right)^{2}}+\frac{C_{2}}{1+\left(R_{2} C_{2} \omega\right)^{2}}+\cdots+\right.\right.\\ &\left.\frac{C_{n}}{1+\left(R_{n} C_{n} \omega\right)^{2}}\right)-j\left(\frac{R_{1} C_{1}^{2} \omega}{1+\left(R_{1} C_{1} \omega\right)^{2}}+\frac{R_{2} C_{2}^{2} \omega}{1+\left(R_{2} C_{2} \omega\right)^{2}}\right.\\ &\left.\left.+\cdots+\frac{R_{n} C_{n}^{2} \omega}{1+\left(R_{n} C_{n} \omega\right)^{2}}+\frac{1}{R_{0} \omega}\right)\right] \end{aligned} $ (6)

式(6)后面的表达式可以看作为复电容Ci表达式,那么有

$ \begin{aligned} C_{\rm i}=&\left(C_{\mathrm{P}}+\frac{C_{1}}{1+\left(R_{1} C_{1} \omega\right)^{2}}+\frac{C_{2}}{1+\left(R_{2} C_{2} \omega\right)^{2}}+\cdots+\right.\\ &\left.\frac{C_{n}}{1+\left(R_{n} C_{n} \omega\right)^{2}}\right)-j\left(\frac{R_{1} C_{1}^{2} \omega}{1+\left(R_{1} C_{1} \omega\right)^{2}}+\frac{R_{2} C_{2}^{2} \omega}{1+\left(R_{2} C_{2} \omega\right)^{2}}\right.\\ &\left.+\cdots+\frac{R_{n} C_{n}^{2} \omega}{1+\left(R_{n} C_{n} \omega\right)^{2}}+\frac{1}{\omega R_{0}}\right) \end{aligned} $ (7)
2 模型和实验数据对比 2.1 模型和岩矿数据对比

Pelton分别测量了黄铜矿、黄铁矿(含氯化钠溶液)、电阻电路电流电压相位角及阻抗对频率的关系[1]表 1).取一定电容、电阻值(表 2),利用式(1)得到的相位和阻抗能够与实测数据得到很好的吻合(图 12).

表 1 黄铜矿、黄铁矿相位角阻抗测试数据 Table 1 Test data of phase angle and impedance for chalcopyrite and pyrite
表 2 黄铜矿、黄铁矿电容电阻取值表 Table 2 Capacitance and resistance values of chalcopyrite and pyrite
图 1 黄铜矿相位和阻抗曲线(双对数) Fig.1 Phase and impedance curves of chalcopyrite (double logarithm) 1—计算相位(calculated phase);2—计算阻抗(calculated impedance);3—实测相位(measured phase);4—实测阻抗(measured impedance)
图 2 黄铁矿相位和阻抗曲线(双对数) Fig.2 Phase and impedance curves of pyrite (double logarithm) 1—计算相位(calculating phase);2—计算阻抗(calculating impedance);3—实测相位(measured phase);4—实测阻抗(measured impedance)

由式(2)可以看到,当频率趋于0时,等效电容是并联电容的和,电阻为漏电电阻.当频率趋于无穷时,等效电容为0,阻抗倒数为内阻及漏电电阻倒数的和.这是电容频率趋于0时隔电,趋于无穷时相当于短路的特征. Pelton实测结果证实了这一点.

2.2 模型和玻璃、云母、琥珀数据对比

R. Kohlrauch测得了玻璃充电吸收电流数据,H. Neumann测得了琥珀的充放电数据(表 3[5],Voglis测得了云母的充放电数据(表 4[5].利用式(3),取一定数值(表 5~7),可以与实测数据得到很好的弥合(图 3~5).

表 3 玻璃、琥珀充放电测试数据 Table 3 Charge-discharge test data of glass and amber
表 4 云母放电测试数据 Table 4 Discharge test data of mica
表 5 玻璃充电时间常数τ、电流Ι取值表 Table 5 Charging time constant and current values of glass
表 6 云母放电时间常数τ、电流Ι取值表 Table 6 Discharge time constant and current values of mica
表 7 琥珀充放电时间常数τ、电流Ι取值表 Table 7 Charge-discharge time constant and current values of amber
图 3 玻璃充电曲线(双对数) Fig.3 Charging curve of glass (double logarithm) 1—计算充电损失(calculated charging loss);2—实测充电损失(measured charging loss)
图 4 琥珀充电和放电曲线(双对数) Fig.4 Charging and discharging curve of amber (double logarithm) 1—计算充电损失(calculated charging loss);2—实测放电损失(measured disharging loss);3—实测充电损失(measured charging loss)
图 5 云母放电曲线(双对数) Fig.5 Discharging curves of mica (double logarithm) 1—506 V计算曲线(calculated curve of 506 V);2—104 V计算曲线(calculated curve of 104 V);3—506 V放电实测(measured discharge of 506 V);4—104 V放电实测(measured discharge of 104 V)

Vollis用两个电压测量了云母的放电电流. 图 5中可以看到两条曲线形态完全一致.由于云母电阻电容不变,变化的只有电压,因此曲线完全相同.式(3)清楚地表明了这一点.注意到琥珀弛豫时间在30 min左右,这表明同一模型同样适用于短时间和长时间的弛豫现象.

2.3 模型和聚乙烯胶囊数据对比

Asami利用介电扫描显微镜测量了聚乙烯胶囊的介电常数、电损ε"、电导等数据(表 8[16].利用式(7)(5)取一定数值(表 9),能够很好地吻合这些数据(图 6~9).

表 8 聚乙烯胶囊测试数据 Table 8 Test data of polyethylene capsule
表 9 聚乙烯胶囊电容电阻取值表 Table 9 Capacitance and resistance values of polyethylene capsule
图 6 聚乙烯胶囊介电常数电损曲线 Fig.6 Permittivity-electrical loss curve of polyethylene capsule 1—电损耗计算(calculated electrical loss);2—电损耗实测(measured electrical loss)
图 7 聚乙烯胶囊电导容纳曲线 Fig.7 Conductivity and capacitance curve of polyethylene capsule 1—容纳率计算曲线(calculated capacitance curve);2—容纳率实测(measured capacitance)
图 8 聚乙烯胶囊介电常数频率曲线 Fig.8 Permittivity-frequency curve of polyethylene capsule 1—等效电容计算曲线(calculated equivalent capacitance);2—介电常数实测(measured permittivity)
图 9 聚乙烯胶囊电导频率曲线 Fig.9 Conductivity-frequency curve of polyethylene capsule 1—电导率计算曲线(calculated conductivity curve);2—电导率实测(measured conductivity)

图 79中电导最小值为介质电导(ω趋于0),电导最大值为介质电导与介质电容内导的和(ω趋于无穷).最小介电常数和最大介电常数分别是平行板电容CP、平行板电容与介质电容的和(CP+C)与真空介电常数的比值.这与式(5)(7)一致.

2.4 经典介电常数曲线和模型对比

利用式(7)(略去$\frac{{1}}{{ωR_0}}$),取一定数值(表 10),可以得到不同等效电容Ce和周期电耗Cε关系(图 10~13).

表 10 介质电容电阻取值表 Table 10 Capacitance and resistance values of dielectric
图 10 多峰等效电容Ce和周期电耗Cε Fig.10 Multimodal curve of equivalent capacitance vs. periodic electrical loss
图 11 斜歪等效电容Ce和周期电耗Cε曲线 Fig.11 Slope curve of equivalent capacitance vs. periodic electrical loss
图 12 弧状等效电容Ce和周期电耗Cε Fig.12 Arc curve of equivalent capacitance vs. periodic electrical loss
图 13 半圆等效电容Ce和周期电耗Cε Fig.13 Semicircle curve of equivalent capacitance vs. periodic electrical loss

结合表 10,可以看到,当电容、电阻值差距较大时,为多峰的曲线,随差距减小.曲线由斜歪到压扁的半圆,直至差距为0,关系曲线变为一个半圆.差距为0时,式(5)可以写成

$ Y=j \omega C_{\mathrm{P}}+j \omega \frac{C}{1+j R C \omega}+\frac{1}{R_{0}} $ (8)

$ C_{\rm i}=C_{\mathrm{P}}+\left(\frac{C}{1+j R C \omega}+\frac{1}{j \omega R_{0}}\right) $ (9)

C为介质的并联电容的和,R为电容内阻并联电路的电阻,R0很大时,可略去$\frac{{1}}{{jωR_0}}$.用时间常数τ替代RC,等式两边比上真空介电常数,式(9)实际上就变为Debye的介电常数的经验公式.参照式(5),单位电压下电损耗功率数值上等于导纳的实部.这说明Debye经验公式的电损ε"并不是实际电损耗,实际是单位周期电耗和2π的比值.对比表明,介质电容模型完全能解释经典介质阻抗、电容性实验数据.

3 温度、频率对等效电容影响的推论

介电常数会随温度变化[17].式(7)的实部为介质等效电容.可以看到电阻变大,等效电容减小.温度变化会引起电阻变化.对于温度系数为正数的介质,温度增加,等效电容减小;对于温度系数为负数的介质,温度增加,等效电容增加.另外,在特定频率下,可能由于共振等方面的影响,也可能会引起电阻的突变,从而导致等效电容的变化.这可能是介电谱在高频下产生跳跃的原因.

4 讨论

平行板电容CP在有介质的情况下大于真空下电容的讨论.实践知道,平板电容在加入导体后,电容会增加,一般认为导体相当于缩小了平行板间的距离,电场增加,电容变大.介质实际上也是导体,同样也会使电场变大,而使电容增加.另外,介质有一定折射率,一般会使电磁波的速度变慢.如果认为由于介质通过影响电场通导(场导)而使电磁波速度变慢的话,那么“场导”同样会影响平行板间的电场,它应该会使电容减小.

由式(6)知道,电路表现的电容为介质等效电容,它是平行板电容与介质等效电容的和.介电常数是这个等效电容和与真空介电常数的比值.由此可见,如果存在介电常数的话,它也是等效介电常数,并不是介电常数随频率变化.

目前认为,加入介质后电容增加是由于极化的作用.介质极化理论认为介质在外电场的作用产生极化电场,其方向与外电场方向相反,使总电场强度减小.但电容是和电场强度成正比的,电场强度减小不应引起电容增加.另外,当外电场消失后,与电场方向相反的介质电场应该产生由电容负极流向正极的短接放电电流,这与实际相反.再有,含介质电容在充电过程中有一吸收电流,并被认为是被介质吸收.极化理论并没有要求介质和平行板间有导体连接,绝缘的情况下,外电源电流不会流到介质里.介质也就不可能吸收电流.

5 结论

介质是存在电容的.没有平行板电容情况下,介质具有电容特性;在有平行板电容情况下,介质同样具有电容特性.

弛豫时间是电容时间常数决定的.介质成分均一,并联电容的电容、内阻都相等的情况下,时间常数就是弛豫时间.

介质电容内阻是可测量的.充电瞬间可以看作电容是短路的,其电导等于内电导和介质电导的和.测得漏电电导和充放电瞬时电导即可得到电容内电导.交流电情况下,当频率趋于无穷时,相当于电容短路.同样可以得到内电导.实际上,本文已经得到了黄铜矿、黄铁矿(含其他介质)及聚乙烯(含其他介质)的电容内阻.

吸收电流是介质电容充电电流.吸收电流和短接放电电流是镜像对应的,它反应了介质电容的电量储存能力.这个吸收电流产生无功功率,同时由于通过电容内阻而产生有功功率.

平行板会使介质电容的电容增加.平行板电极会使介质的外电场增加,因此介质电容会增加.

介质电损耗在单位电压下,其数值上等于介质电导与介质电容内导的和.

岩石电容是电介质电容的一种.

参考文献
[1]
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.
[2]
商木元. 岩石电容:物探"激发极化"的实质——对于物探"激发极化"一次和二次电压现象的讨论[J]. 地质与资源, 2019, 28(1): 98-108.
[3]
赵孔双. 稀薄细胞悬浊系模型的介电解析法[J]. 高等学校化学学报, 1997, 18(1): 107-111.
[4]
赵孔双. 介电谱方法及应用[M]. 北京: 化学工业出版社, 2008: 91.
[5]
Cole K S, Cole R H. Dispersion and absorption in dielectrics Ⅰ. Alternating current characteristics[J]. The Journal of Chemical Physics, 1941, 9(4): 341.
[6]
Cole K S, Cole R H. Dispersion and absorption in dielectrics Ⅱ. Direct current characteristics[J]. The Journal of Chemical Physics, 1942, 10(2): 98.
[7]
Havriliak S, Negami S. A complex plane representation of dielectric and mechanical relaxation processes in some polymers[J]. Polymer, 1967, 8: 161-210.
[8]
Cole R H. Evaluation of dielectric behavior by time domain spectroscopy. Ⅰ. Dielectric response by real time analysis[J]. The Journal of Physical Chemistry, 1975, 79(14): 1459-1469.
[9]
Cole R H. Evaluation of dielectric behavior by time domain spectroscopy. Ⅱ. Complex permittivity[J]. The Journal of Physical Chemistry, 1975, 79(14): 1469-1474.
[10]
David D W, Cole R H. Dielectric relaxation in glycerol, propylene glycol, and n-propanol[J]. The Journal of Chemical Physics, 1951, 19(12): 1484.
[11]
刘国兴. 电法勘探原理与方法[M]. 北京: 地质出版社, 2005: 108.
[12]
程志平. 电法勘探教程[M]. 北京: 冶金工业出版社, 2007: 89.
[13]
傅良魁. 电法勘探教程[M]. 北京: 地质出版社, 1983: 187.
[14]
邱关源. 电路[M]. 5版. 北京: 高等教育出版社, 2006: 202.
[15]
赵凯华, 陈熙谋. 电磁学[M]. 北京: 高等教育出版社, 1985: 155.
[16]
Asami K, Zhao K S. Dielectric measurement of a single sub-millimeter size microcapsule[J]. Colloid and Polymer Science, 1994, 272(1): 64-71.
[17]
陈季丹, 刘子玉. 电介质物理学[M]. 北京: 机械工业出版社, 1982: 132.