劉子強，齊興華，王有龍

（1.新疆煤田地質(zhì)局一五六勘探隊 新疆 烏魯木齊 830023；2.新疆工程學院 新疆 烏魯木齊 830022）

高密度電阻率法（HDR）是一種網(wǎng)絡勘探方法。根據(jù)巖石和土壤導電性的差異，研究了人工穩(wěn)定電流場作用下土壤中導電電流的分布。現(xiàn)場測量時，所有電極（幾十到幾百個）都可以放置在觀測段的測量點上。然后利用程控電極轉(zhuǎn)換裝置和微機工程機電測量儀實現(xiàn)數(shù)據(jù)的快速自動采集。將測量結(jié)果發(fā)送到微機后，還可以對數(shù)據(jù)進行處理，并給出地電剖面分布的圖形結(jié)果[1-2]。

電阻率采集過程中，電極系系數(shù)是將測得的電阻值轉(zhuǎn)換為電阻率值的關(guān)鍵系數(shù)。一般認為，該系數(shù)主要由電極間的相互位置決定，其中α排列的電極系系數(shù)理論計算公式如下：

式中：K為電極系系數(shù)，單位為m；AM、AN、BM、BN為電極間距，單位為m。

如公式（1）所示，電極系系數(shù)K是一個僅僅與各電極之間的距離相關(guān)的數(shù)值，但當電阻率測量裝置形式和極距確定，K值便可以通過公式進行計算[3-4]。但是在實際使用過程中，所使用的電極材料和大小不同的電阻率剖面儀制造商不同和地層的深度是不同的，由于人為因素和環(huán)境因素，有時可能使用不同的供電方法或電極并非完全插入到地層[5]。這些都可能導致實際電極系系數(shù)與理論電極系系數(shù)之間存在差異，導致所計算出的視電阻率與真電阻率之間的誤差較大。

為了使計算得出的視電阻率與真電阻率之間的誤差變小，有必要對各種因素對實際電極系系數(shù)的影響開展研究，為校正影響奠定基礎(chǔ)。因此，提出了一種基于COMSOL的模擬計算方法[6-8]，利用有限元法建立一種電阻率法的三維正演模型對實際電極系系數(shù)進行模擬計算，重點分析了不同電極半徑、電極插入地層深度和不同電極間距的條件下對實際電極系系數(shù)的影響，從而加深從業(yè)者對電極系系數(shù)的理解，對在生產(chǎn)用采用高密度電法及與之類似的其他電阻率測量方法時能夠提高探測精度，有較好的參考價值。

1 模型建立

文章從3個方面進行了論述。第一，在保持電極材料、電極間距、電極埋入地層深度等因素不變的情況下，通過改變電極半徑來探討電極半徑對實際電極系系數(shù)的影響。第二，在保持電極材料、電極間距、電極半徑等各因素不變的情況下，只改變電極插入地層的深度，以探討電極插入地層深度對實際電極系系數(shù)的影響。第三，在保持電極材料、電極插入地層深度、電極半徑等各因素不變的情況下，通過調(diào)整電極間距來探討電極間距對實際電極系系數(shù)的影響。

在解決簡單模型的模擬問題時，使用偏微分方程即可滿足需求，但是當幾何結(jié)構(gòu)或者物理場復雜時，偏微分方程便無法取得解析解。將復雜模型進行離散化從而構(gòu)建近似方程并獲取近似解以達到對復雜模型模擬的目的即是有限元法的本質(zhì)。COMSOL Multiphysics是一款基于有限元法的模擬軟件，近年來被各領(lǐng)域?qū)W者廣泛使用，取得了良好的效果，本文即采用該軟件開展上述研究。

1.1 理論基礎(chǔ)

在地球物理勘探中，經(jīng)常通過對目標地質(zhì)體電阻率的測量來完成對其大小、形狀、礦物種類、含量等多方面性質(zhì)的判斷。其中，經(jīng)常用到的一種方法就是電阻率剖面法。電阻率剖面法是指A、M、N、B電極距保持不變，同時沿—定剖面方向逐點觀測視電阻率（ρs），研究剖面方向地下一定深度的巖、礦石電阻率變化的一組方法。

在地下介質(zhì)均勻各向同性的條件下，電阻率計算式：

式中：R為電阻率，單位為Ω·m；

U為電壓值，單位為V；

I為電流值，單位為A。

由公式（2）可知在模型中只要能夠獲取到除K以外的物理量則可求出實際的K值。

對于時變電場和電流，在磁效應可以忽略不計的情況下，物理場應遵循電流守恒方程：

式中：J為電流密度，單位A/m2；

ρ為電荷密度，單位C/m3；

t為時間，單位為s。

其中

式中：σ為電導率，單位為S/m；

V為電壓，單位為V（伏特）；

E為電場強度，單位為V/m。

1.2 模型參數(shù)

模型由地層和電極棒組成，如下圖1所示，其中正方體為地層，其長寬高均為100 m，其電導率設(shè)定為1 S/m，即電阻率為1 Ω·m。其上4個小圓柱為銅電極，電極直徑初始0.004 m，電極初始間距2 m，自左起依次為A、M、N、B電極，其中A電極和B電極分別為正負供電電極，向地層中供直流電流，M電極和N電極為測量電極，測量所處位置的電壓。模型中設(shè)定A電極與B電極之間的電流恒定為1A，并通過探針的形式測量M點電壓及N點電壓。

圖1 模型示意圖

為滿足使用有限元法進行模擬的精確度，在地層中所使用的網(wǎng)格最大為1 m，在電極上采用的網(wǎng)格最大為0.001 m，網(wǎng)格剖分后的結(jié)果如圖2所示。為保證模型計算精度的可靠，特采用了100 m尺寸的地層模型，保證地層的尺寸在研究范圍內(nèi)足夠大以減小電場邊界效應對模擬結(jié)果的影響，因此，對模型的邊界未進行虛擬域或其他邊界設(shè)置。

圖2 網(wǎng)格剖分結(jié)果

2 模擬結(jié)果及分析

初始狀態(tài)下的電場分布情況如圖3所示。由圖3可以看出，A電極電勢最大，B電極電勢最小，M電極和N電極電勢差別難以用肉眼觀察但軟件可通過探針進行準確探測。

圖3 電場分布圖

在此情況下理論電極系系數(shù)計算過程如下：

AM=2 m，AN=2 m，BM=2 m，BN=2 m代入公式（1）中，得出理論電極系系數(shù)K=12.56，單位為m。

實際電極系系數(shù)計算過程如下：

地層電阻率R為1 Ω·m，A電極與B電極之間的電流I為1 A，模擬結(jié)果顯示M電的電壓為1.8313 V，N點的電壓為1.754 7 V，因此M、N電極的電壓差U為0.076 6 V。將上述物理量帶入公式（2），可知此時的電極系系數(shù)為：

由上述計算過程可以看出，在初始情況下，實際電極系系數(shù)與理論電極系系數(shù)并不完全一致，存在由各種因素導致電極系系數(shù)有較大變化從而影響電阻率計算精度降低的可能性。

2.1 當電極半徑改變時的模擬結(jié)果

在其他條件均不改變時，電極半徑從0.002 m逐步增加到0.007 m時，通過仿真計算得到的實際電極系系數(shù)的變化，并與由（1）式得到的理論電極系系數(shù)進行比較，其結(jié)果見圖4。

圖4 電極半徑對電極系系數(shù)的影響變化模擬結(jié)果

從圖4中可以看出，實際電極系系數(shù)與理論電極系系數(shù)的差值隨著電極半徑的增大而增大。但是從變化的程度可以看出，電極的半徑變化對電極系系數(shù)的影響非常有限，因此在實際使用過程中，采用常用的電極基本不需要對電極系系數(shù)進行半徑方面的校正。

2.2 當電極插入地層深度改變時的模擬結(jié)果

當電極插入地層深度逐漸由0.1 m逐漸變?yōu)?.5 m且其他條件均不改變時，得出的實際電極系系數(shù)的變化并與由公式（1）得出的理論電極系系數(shù)作對比如圖5所示。

圖5 電極插入地層深度對電極系系數(shù)的影響變化模擬結(jié)果

從圖5中可以看出，實際電極系系數(shù)隨著電極插入地層深度的增大而增大，且與理論電極系系數(shù)的差值越來越大，因此，對電極插入地層深度為0.45 m時的電場分布情況（如圖6所示）開展了相關(guān)研究。

圖6 插入深度為0.45 m時的電場分布圖

通過圖6可看出，A電極供電是電勢最大，B電極電勢最小，M電極和N電極電勢近乎相等，但與圖3未改變地層深度時的電場分布情況相比，電極附近的電場分布情況發(fā)生了明顯的改變，電場強度最高的地方與電極的深入程度直接相關(guān)，因此，隨著電極插入程度的加深，電場的整體分布發(fā)生了與理想電極系情況下的點電極電場分布有了很明顯的區(qū)別，這種電場分布的變化最終導致了實際電極系系數(shù)與理論電極系系數(shù)之間的差別隨著插入深度的變化而發(fā)生明顯變化。

2.3 當電極間距改變時的模擬結(jié)果

在其他條件不變的情況下，在電極間距從2 m逐步增加到6 m時，得到了實際的電極系系數(shù)的變化，并將其與公式（1）中的理論電極系系數(shù)進行比較如圖7所示。

圖7 電極間距對電極系系數(shù)的影響變化模擬結(jié)果

從圖7中可以看出，當電極距增大時，實際電極系系數(shù)與理論電極系系數(shù)之間的差值逐漸減小，在達到一定的電極距之后二者基本趨于相同。這是由于在電極尺寸和插入地層深度不變的情況下，隨著電極距的增大，整個電極系的尺寸在不斷的增大，因此，電極的相對尺寸越來越小，即電極系的工作狀態(tài)越來越趨近于點電極時的工作狀態(tài)，從而導致實際電極系系數(shù)與理論電極系系數(shù)越來越接近。這也說明，在使用較大電極距時，電極系系數(shù)可以近似采用理論電極系系數(shù)，而在電極距較小時，則應對電極系系數(shù)進行相應的校正以降低視電阻率的計算誤差。

3 結(jié)論

本文通過使用COMSOL Multiphysics有限元模擬軟件分析在不同電極半徑、電極插入地層深度和不同電極間距的條件下對實際電極系系數(shù)的影響及模擬計算，得出的結(jié)論如下：

（1）電極半徑對實際電極系系數(shù)的影響非常小，使用常用電極時無須對電極半徑進行電極系系數(shù)的校正。

（2）隨著電極插入地層深度的增大，實際電極系系數(shù)與理論電極系系數(shù)之間的差值將明顯變大，體現(xiàn)出插入深度對實際電極系系數(shù)的影響較為明顯。因此，有必要在計算視電阻率的時候?qū)﹄姌O插入地層的深度進行實際電極系系數(shù)的校正。

（3）隨著電極間距的增加實際電極系系數(shù)與理論電極系系數(shù)之間的差值將越來越小并逐漸趨于一致。當電極插入深度不深且電極距較大時，可以直接采用理論電極系數(shù)而無須進行校正。