本發(fā)明涉及一種適用于復(fù)雜地形條件下的高密度地電阻率測量系統(tǒng)及方法。

背景技術(shù)：

大地電阻率的測量被廣泛用于尋找地下水源，解決人畜用水及工農(nóng)業(yè)用水問題；用于水文工程，環(huán)境的地質(zhì)勘探及高分辨電阻率法工程地質(zhì)勘探；用于金屬與非金屬礦產(chǎn)資源勘探，能源勘探，城市物探，鐵道及橋梁工程勘探，并用于找地?zé)?，確定水庫壩基和防洪大堤隱患位置。

目前，現(xiàn)有的高密度地電阻率測量儀基本是多功能直流電法儀配合多路電極轉(zhuǎn)換器的方式。每根電極既做供電電極又做測量電極，供電后極化長時間不穩(wěn)定，很難達(dá)到快速高精度檢測；此外，測量場地周圍存在空間電磁干擾時會影響測量數(shù)據(jù)甚至是采集到錯誤數(shù)據(jù)，對后期數(shù)據(jù)反演時帶來反演誤差，甚至是得到錯誤結(jié)果；除此之外，現(xiàn)在的地電阻率測量儀器在考慮溫度、濕度測量影響方面，要么是非高密度儀器具有溫度、濕度測量功能，在高密度類型的地電阻率測量儀器中往往未考慮溫度、濕度對地電阻率的影響。對此，特別是針對復(fù)雜地形而言，這兩個環(huán)境參數(shù)是影響地電阻率的重要影響因素，因此急需一種針對復(fù)雜地形條件的地電阻率測量儀器。

技術(shù)實現(xiàn)要素：

本發(fā)明為了解決現(xiàn)有的高密度地電阻率測量儀不能適應(yīng)復(fù)雜地質(zhì)條件下的低電阻率測量的問題，提出了一種適用于復(fù)雜地形條件下的高密度地電阻率測量系統(tǒng)及方法，本發(fā)明抗電磁干擾能力強，測量快速、靈敏度高。

為了實現(xiàn)上述目的，本發(fā)明采用如下技術(shù)方案：

一種適用于復(fù)雜地形條件下的高密度地電阻率測量系統(tǒng)，包括主控單元、電極轉(zhuǎn)換單元和全橋電路，所述主控單元通過電極轉(zhuǎn)換單元控制多個電極模塊的高壓供電與測量狀態(tài)，全橋電路產(chǎn)生交流信號，并通過激發(fā)電極模塊注入電流到大地；

所述主控單元通過電極轉(zhuǎn)換單元控制電極模塊的排列方式、極距和測點掃描，對電極模塊采集的電流、電壓、溫度、濕度和經(jīng)緯度信號進行圓滑處理，構(gòu)建反演圖，得到復(fù)雜地形條件的地址情況。

所述主控單元，包括控制模塊、通信模塊和數(shù)據(jù)采集模塊，所述控制模塊通過通信模塊與電極轉(zhuǎn)換單元連接并通信，發(fā)送電極模塊的控制命令，所處數(shù)據(jù)采集模塊包括兩條模數(shù)轉(zhuǎn)換通道，分別采集電壓與電流數(shù)據(jù)。

所述電極轉(zhuǎn)換單元，包括解碼模塊、串轉(zhuǎn)并模塊、驅(qū)動模塊、繼電器陣列和電纜接口，繼電器陣列由多個單刀雙擲開關(guān)繼電器組成，每4個繼電器為一組，控制一個電極模塊，電纜接口為具有多個插孔的航空插座，每一個插孔通過導(dǎo)線連接到繼電器陣列的一組繼電器；解碼模塊接收控制模塊發(fā)送的信號，將其解碼后，利用SPI總線將解碼后數(shù)據(jù)發(fā)送給串轉(zhuǎn)并模塊，串轉(zhuǎn)并模塊對數(shù)據(jù)進行轉(zhuǎn)化輸出，輸出信號送入驅(qū)動模塊以驅(qū)動繼電器陣列中的繼電器打開與關(guān)閉，選擇各個電極模塊的接通與斷開。

所述電極模塊，包含不銹鋼電極、溫濕度測量模塊和GPS模塊，不銹鋼電極直接與大地接觸，用于向大地注入電流和測量兩個電極之間的電壓，溫濕度模塊用于測量當(dāng)前狀態(tài)的環(huán)境溫度與濕度，GPS模塊用于記錄測量電極的經(jīng)緯度。

每一次測量過程有4個電極模塊處于工作狀態(tài)，其中兩個電極模塊向大地注入電流，兩個電極模塊測量兩點之間的電位差。

所述主控單元通過電纜與各個電極模塊連接，所述電纜為內(nèi)含多根導(dǎo)線的電纜和多個針頭的航空插頭；

所述電纜的導(dǎo)線數(shù)量與繼電器數(shù)量一致，針頭數(shù)量與插孔數(shù)量一致。

基于上述系統(tǒng)的工作方法，包括以下步驟：

(1)根據(jù)初始測量值計算接地電阻和偏置電壓的值，并根據(jù)信號大小選擇放大倍數(shù)與偏置電壓；

(2)采集數(shù)據(jù)，對信號過沖和隨機噪聲干擾進行抑制和剔除，提取有效信號，轉(zhuǎn)換成相應(yīng)的二維數(shù)據(jù)反演格式；

(3)加入地形及不均勻體數(shù)據(jù)，加入現(xiàn)場溫濕度數(shù)據(jù)修正，利用二維數(shù)據(jù)反演軟件進行二維數(shù)據(jù)反演，使用最小二乘法等迭代方法計算，直到滿足設(shè)定值誤差停止；

(4)對反演圖結(jié)合地質(zhì)情況進行分析，分析各種地質(zhì)情況。

所述步驟(1)中，電極模塊按照等間距間隔插入大地，通過電纜連接主控單元和電極模塊。

所述步驟(1)中，第一次測量分兩步進行，第一步測量接地電阻，第二步測量偏置電壓；在第一步測量完成后，根據(jù)電壓的值大小選擇放大倍數(shù)，接地電阻值的計算是通過計算測量兩點之間的電位差的電極模塊的供電電壓和供電電流的關(guān)系得到，第二步測量完后，根據(jù)測得電壓的最大值和最小值計算得出信號的偏置電壓。

所述步驟(2)中，進行數(shù)字平均累加和信號過沖剔除的數(shù)據(jù)處理。

所述步驟(3)中地形數(shù)據(jù)，通過室內(nèi)試驗，模擬計算近地表局部電性不均體對視電阻率影響規(guī)律，并進行導(dǎo)電紙試驗，在對數(shù)據(jù)進行仔細(xì)對比分析后，發(fā)現(xiàn)采用溫納裝置和偶極裝置在不同地形條件下可能出現(xiàn)相差較大的數(shù)值，而溫納裝置和斯倫貝爾裝置所測得的數(shù)據(jù)趨勢則相同。在精度方面，該發(fā)明高密度測量系統(tǒng)溫納裝置的精度最高，且非常穩(wěn)定，其次是斯倫貝爾裝置，偶極裝置的精度較差并且不穩(wěn)定。故在野外復(fù)雜地電條件的測量中，建議使用溫納裝置進行測量。

所述步驟(3)中溫濕度數(shù)據(jù)，通過儀器自帶的溫濕度探頭測量，針對溫濕度數(shù)據(jù)的反演也進行了室內(nèi)實驗以及野外試驗。

溫度變化對土壤電阻率的影響研究：由于電阻率的影響因素存在耦合效應(yīng)，整體研究相對困難，所以，儀器先選擇溫度這一單一影響因素來觀察其對電阻率的影響，通過理論分析與實驗室物理模擬，研究溫度變化對土壤電阻率的影響規(guī)律。使用控制變量法對溫度影響因素進行物理模擬試驗。針對溫度影響因素的特點，選擇適當(dāng)含水率的不同土壤進行試驗，改變試驗土壤的溫度，測量出對應(yīng)的電阻率隨溫度的變化曲線，并對曲線加以分析、擬合，得到電阻率受溫度的影響規(guī)律從而將其應(yīng)用到實際之中。

濕度變化對土壤電阻率的影響研究：濕度對土壤電阻率的影響研究也是利用控制變量法，保持其他因素不變，僅僅改變土壤的濕度并對其電阻率進行測量。選擇適當(dāng)?shù)姆椒▽y量結(jié)果進行擬合，分析濕度變化對土壤電阻率的影響規(guī)律。

所述步驟(3)中地形及局部不均勻體對視電阻率的影響研究：地形及不均勻體對視電阻率的影響研究采用的是二維導(dǎo)電紙模擬試驗。將導(dǎo)電紙裁剪成不同的地形形狀或者在水平地下制作不均勻體，并使用成熟商業(yè)儀器進行數(shù)據(jù)的采集工作，得到不同測量裝置對應(yīng)的高密度剖面，對視電阻率剖面等值線圖進行分析研究，得到地形和不均勻體對視電阻率的影響規(guī)律。

所述步驟(3)中二維數(shù)據(jù)反演軟件進行二維數(shù)據(jù)反演，包括直流電法和低頻交流電法二維正反演軟件設(shè)計的原理推導(dǎo)、算法研究、程序設(shè)計和模型實驗，并對程序的運算能力和程序的計算精度作了檢測和評價。通過正演研究，求得地面或地下的電位分布，由于地電條件的復(fù)雜性，一般場論中的經(jīng)典方法難于求得這些問題的解，除了地電條件比較規(guī)則簡單的少數(shù)問題以外，大量的實際電法探測問題目前尚未得出準(zhǔn)確的解析解。因此，該發(fā)明從電場滿足的微分方程出發(fā)，用近似方法求其數(shù)值解。在電阻率正演計算中，實際的地質(zhì)構(gòu)造都是三維的，電流在地下流動的規(guī)律也是三維的，三維正演模擬更能反映復(fù)雜地質(zhì)構(gòu)造的視電阻率分布特征，但需要消耗大量的計算資源。當(dāng)假設(shè)地質(zhì)體有一個維度是無限延伸的，將具有一定走向的三維穩(wěn)定電流場轉(zhuǎn)化為2.5D問題，在走向上進行傅里葉變換，使得原來的空間域(x,y,z)問題轉(zhuǎn)化為波數(shù)域(kx,y,z)問題，三維問題就可以轉(zhuǎn)化為二維問題求解，計算速度將大幅度提高。本發(fā)明反演過程中，要求解的問題通常是病態(tài)的，正則化參數(shù)的應(yīng)用能有效改善方程組的病態(tài)問題。為提高反演質(zhì)量，確定比較合理的模型，如何選取合適的正則化因子往往是地球物理反演工作研究的重要內(nèi)容。利用正則化方法是減少反演的多解性和提高反演過程的穩(wěn)定性最為有效的手段，正則化因子值的選取直接影響反演結(jié)果的精度和反演過程的穩(wěn)定性。

本發(fā)明的有益效果為：

(1)本發(fā)明利用低頻交流供電技術(shù)，針對具體的復(fù)雜環(huán)境選擇不同的激勵信號頻點，對當(dāng)前地電阻率進行測量，為在強空間電磁干擾的環(huán)境下進行測量提供了一種新型的有效方法。

(2)本發(fā)明能根據(jù)需要選擇不同測量方法，確保測量數(shù)據(jù)的有效性和可靠性；電極轉(zhuǎn)換裝置可根據(jù)需要設(shè)置多個頻點，根據(jù)工作需要進行靈活選用，實現(xiàn)電阻率和交流激電功能；

(3)本發(fā)明適用于復(fù)雜地質(zhì)條件下地電阻率測量：利用GPS可測量出每一個電極的經(jīng)緯度和高度，并將各電極的相對高度信息融入的反演算法中，可進行復(fù)雜地形地電阻率測量反演；

(4)本發(fā)明適用于復(fù)雜地質(zhì)條件下地電阻率測量：利用溫濕度數(shù)據(jù)得到經(jīng)驗公式加入正反演算法中，修正原始數(shù)據(jù)，并可以根據(jù)不同條件反演推斷相應(yīng)的電阻率；

(5)本發(fā)明抗干擾能力強：利用H橋產(chǎn)生交流電流注入大地，通過兩點大地之間的交流電壓信號，利用交流電阻率測量方法計算視電阻率，比直流測量方法具有更強的干擾能力；

(6)本發(fā)明精度高：地電阻率的大小與溫度、濕度均有較大關(guān)系，可記錄每一個電極所處位置的環(huán)境溫度、濕度情況，將這兩個參數(shù)融入到反演方法中，可得到更為準(zhǔn)確的反演結(jié)果；除此之外，儀器中利用24位模數(shù)轉(zhuǎn)換器進行數(shù)據(jù)采集，在無干擾條件下測量精度可達(dá)到10uV；

(7)本發(fā)明具有結(jié)構(gòu)簡單、成本低、功耗低等特點，可長時間進行野外工作，同時可靠性好，不受降雨、大霧等天氣因素的影響，易于維護。

附圖說明

圖1為本發(fā)明的整體結(jié)構(gòu)示意圖；

圖2為本發(fā)明的主控采集模塊示意圖；

圖3為本發(fā)明的電極控制模塊示意圖；

圖4為有限單元法正演程序流程圖；

圖5阻尼型高斯-牛頓法反演流程。

其中，1、主控模塊，2、電極控制模塊，3、全橋模塊，4、電極模塊，5、電纜，6、控制器模塊，7、RS422通信模塊，8、數(shù)據(jù)采集模塊,9、解碼模塊，10、串轉(zhuǎn)并模塊,11、驅(qū)動模塊,12、繼電器陣列,13、電纜接口。

具體實施方式：

下面結(jié)合附圖與實施例對本發(fā)明作進一步說明。

如圖1所示，一種用于復(fù)雜地質(zhì)條件下低電阻率測量的高密度電阻率測量儀，包括主控模塊1、電極控制模塊2、全橋模塊3、電極模塊4、電纜5五部分組成，主控模塊為核心，控制其它模塊協(xié)調(diào)工作。

所述主控模塊1，由控制器模塊6、RS422通信模塊7、數(shù)據(jù)采集模塊8組成，控制器模塊6采用STM32單片機為主控單片機，控制器模塊6通過RS422通信模塊7與電極控制模塊2連接并通信，用于將電極控制命令發(fā)送到電極控制模塊2，數(shù)據(jù)采集模塊8是由24位模數(shù)轉(zhuǎn)換芯片AD7710設(shè)計的兩路模數(shù)轉(zhuǎn)換器通道，分別用于電壓、電流采集，控制器模塊6通過SPI接口控制數(shù)據(jù)采集模塊8進行數(shù)據(jù)采集。

所述電極控制模塊2，由解碼模塊9、串轉(zhuǎn)并模塊10、驅(qū)動模塊11、繼電器陣列12和電纜接口13構(gòu)成，串轉(zhuǎn)并模塊10由16片74HC595芯片串聯(lián)而成，驅(qū)動模塊11由16片ULN2803構(gòu)成，繼電器陣列12由128個單刀雙擲開關(guān)繼電器組成，每4個繼電器控制電極模塊4中的1個電極，電纜接口13由1個32插孔的航空插座組成，每一個插孔用1根導(dǎo)線連接到繼電器陣列12中的4個繼電器；解碼模塊9接收主控模塊1發(fā)送的信號，將其解碼后，利用SPI總線將解碼后數(shù)據(jù)發(fā)送給串轉(zhuǎn)并模塊10；串轉(zhuǎn)并模塊10可將16字節(jié)數(shù)據(jù)轉(zhuǎn)換成128位二進制數(shù)據(jù)輸出，該輸出信號送入驅(qū)動模塊11后，用于驅(qū)動繼電器陣列12中的繼電器打開與關(guān)閉，達(dá)到選擇電極模塊4中32個電極的接通與斷開的目的。

所述全橋模塊3，由2片IR4184芯片和4片NMOS管構(gòu)成H橋，用于產(chǎn)生交流信號，并通過電極注入電流到大地。

所述電極模塊4，一共有32個電極，每個電極包含不銹鋼電極、溫濕度測量模塊、GPS模塊，不銹鋼電極直接與大地接觸，用于向大地注入電流和測量兩個電極之間的電壓，溫濕度模塊用于測量當(dāng)前狀態(tài)的環(huán)境溫度與濕度，GPS模塊用于記錄測量電極的經(jīng)緯度；每一次測量過程有4個電極處于工作狀態(tài)，其中兩個用于向大地注入電流，兩個用于測量兩點之間的電位差。

所述電纜5，由一根內(nèi)含32根導(dǎo)線的電纜和32針的航空插頭組成，電纜每隔10米引出1根導(dǎo)線，用于連接電極模塊4中的電極。

本發(fā)明所提供的適用于復(fù)雜地質(zhì)條件下低電阻率測量的高密度電阻率測量儀測量方法是：

1、將電極模塊4中的32個電極按照等間距間隔插入大地，

2、將每個電極于電纜5中的32個抽頭連接，

3、將電纜5中的航空插頭接入電極控制模塊2中的電纜接口13，

4、利用8芯電纜連接將主控模塊1中的RS422通信模塊7和電極控制模塊2中的解碼模塊9連接；

5、在主控模塊中設(shè)置好儀器工作方式，則儀器可按照預(yù)先設(shè)定的方式進行工作。

首先根據(jù)實際情況要求在上層軟件對系統(tǒng)進行初始化設(shè)置，如：采集參數(shù)選擇、串口傳輸數(shù)據(jù)參數(shù)設(shè)置、算法設(shè)置等。

其次，進行第一次測量。根據(jù)此次測量值計算接地電阻和偏置電壓的值，并根據(jù)信號大小選擇放大倍數(shù)。接地電阻的值的計算主要是通過計算電極M、N的供電電壓和供給電極M、N的供電電流的關(guān)系得到。然后可以根據(jù)測得信號的最大值和最小值計算得出信號的偏置電壓。由上層軟件控制單片機將計算得到的偏置電壓加到最前端的儀用放大器上以完成校準(zhǔn)過程。

在做完以上工作以后，可以開始正常測量。將測得的數(shù)據(jù)通過串口傳到上層軟件上進行數(shù)據(jù)處理和數(shù)據(jù)存儲。在實際測量過程中，待測信號常常會受到工頻干擾、信號過沖和隨機噪聲干擾的影響，導(dǎo)致測量信號無法正常提取，因此對以上干擾的處理就顯得十分必要。工頻干擾在前端調(diào)理電路中已經(jīng)得到了抑制，所以就需要在上層軟件處理中對信號過沖和隨機噪聲干擾進行抑制和剔除。根據(jù)實際情況要求，數(shù)據(jù)處理模塊分為數(shù)字平均累加和信號過沖剔除兩部分。

將儀器采集的數(shù)據(jù)通過電法處理軟件對數(shù)據(jù)進行圓滑處理，再轉(zhuǎn)換成相應(yīng)的二維數(shù)據(jù)反演格式，并加入地形數(shù)據(jù)，利用二維數(shù)據(jù)反演軟件進行二維數(shù)據(jù)反演，使用最小二乘法的迭代方法，記錄多次迭代的誤差，誤差在小于10％的時候為滿足要求。最后對反演圖結(jié)合地質(zhì)情況進行分析，分析出各種地質(zhì)情況。

上述雖然結(jié)合附圖對本發(fā)明的具體實施方式進行了描述，但并非對本發(fā)明保護范圍的限制，所屬領(lǐng)域技術(shù)人員應(yīng)該明白，在本發(fā)明的技術(shù)方案的基礎(chǔ)上，本領(lǐng)域技術(shù)人員不需要付出創(chuàng)造性勞動即可做出的各種修改或變形仍在本發(fā)明的保護范圍以內(nèi)。