王國富 龍琴 王小紅 張法全 葉金才

摘 要： 針對傳統(tǒng)的礦井瞬變電磁探測大多以二維視電阻率斷面等值線圖顯示，難以對低阻異常的空間位置及分布范圍做出準(zhǔn)確的判斷，提出一種二維視電阻率斷面等值線圖的三維重建方法。首先利用閾值分割原理將電阻率斷面圖的低阻異常區(qū)與其他區(qū)域分開，然后通過體繪制方法實現(xiàn)低阻異常區(qū)域的三維圖像顯示。實驗結(jié)果表明，三維重建能夠有效顯示垂向視電阻率的變化特征，提高了對低阻異常體空間位置的判定。

關(guān)鍵詞： 視電阻率斷面圖； 礦井瞬變電磁探測； 閾值分割； 三維重建； 等值線； 體繪制

中圖分類號： TN911.7?34 文獻標(biāo)識碼： A 文章編號： 1004?373X（2018）07?0033?03

Three?dimensional reconstruction of apparent resistivity section diagram

WANG Guofu， LONG Qin， WANG Xiaohong， ZHANG Faquan， YE Jincai

（School of Information and Communication， Guilin University of Electronic Technology， Guilin 541004， China）

Abstract： The traditional mine transient electromagnetic detection results are displayed in the mode of isogram of two?dimensional apparent resistivity section diagram， which is difficult to accurately judge the spatial position and distribution range of low resistivity anomaly. For this problem， a three?dimensional reconstruction method of two?dimensional apparent resistivity section isogram is proposed. The threshold segmentation principle is used to separate the low resistivity anomaly area from other areas， and then the volume rendering method is used to realize the three?dimensional image display of low resistivity anomaly area. The experimental results show that the three?dimensional reconstruction method can display the variation characteristics of vertical apparent resistivity effectively， and improve the judgment accuracy for the spatial position with low resistivity abnormal body.

Keywords： apparent resistivity section diagram； mine transient electromagnetic detection； threshold segmentation； 3D reconstruction； isogram； volume rendering

0 引 言

近些年發(fā)生的煤礦安全事故大都造成眾多的工人傷亡及巨大的經(jīng)濟損失[1]，其中含水低阻結(jié)構(gòu)體的存在是造成礦井災(zāi)害的主要因素。礦井瞬變電磁法[2]因其對低阻體敏感，具有較好的探測方向，易攜帶等特性而被廣泛地用于探測礦井中的含水結(jié)構(gòu)體[3?5]。然而，傳統(tǒng)的礦井瞬變電磁探測結(jié)果大部分以二維視電阻率等值線斷面圖的形式解釋，二維等值線圖能夠較清楚地顯示橫向斷面的阻值特征，但無法清晰地顯示垂向延展面上的視電阻率特征分布情況[6]，難以準(zhǔn)確地判斷出異常體的空間位置及分布范圍，特別是垂直方向上的異常體變化特征，因此，三維立體的顯示出低阻體不僅有利于確定礦井中含水體結(jié)構(gòu)的垂直特征，而且有利于物探人員的觀察。

三維重建技術(shù)在醫(yī)學(xué)及工業(yè)領(lǐng)域已得到廣泛應(yīng)用，醫(yī)學(xué)中常用三維重建技術(shù)對人體掃描得到的CT圖像進行處理[7]，由此可以立體地觀察人體器官、骨骼、血管等，可以幫助醫(yī)務(wù)人員準(zhǔn)確判斷病情，快速找出病因，工業(yè)中常用三維重建獲取模型的立體圖形[8]，取得了較好的效果，在此提出將三維重建技術(shù)應(yīng)用到礦井安全防護領(lǐng)域，結(jié)合實際地質(zhì)資料及探測數(shù)據(jù)，對礦井瞬變電磁探測數(shù)據(jù)進行預(yù)處理，采用三維重建技術(shù)解釋礦井瞬變電磁探測成果，實現(xiàn)垂向視電阻率變化的可視化，使探測結(jié)果顯示更直觀，為探測含水構(gòu)造提供了依據(jù)。

1 礦井瞬變電磁原理

礦井瞬變電磁法在井下巷道空間鋪設(shè)發(fā)射線圈并給線圈通以一定的電流，由此發(fā)射回線的周圍空間會形成一個穩(wěn)定的磁場。在某一時刻斷開電流，空間磁場由于一次磁場在巖層中激發(fā)的感應(yīng)電流得以維持而不會立即消失，在一次電流斷開時，處于發(fā)射線圈周圍空間中的礦井巷道頂、底板和側(cè)幫最先感應(yīng)到電流斷開瞬間一次磁場所產(chǎn)生的劇烈變化，因此，由一次磁場激發(fā)的感應(yīng)電流會局限在巷道附近巖層中，而且?guī)r層中各處的感應(yīng)電流分布不均，隨著時間的推移，巖層中感應(yīng)電流強度變?nèi)醪⒙呌诰鶆颉?/p>

研究結(jié)果表明，任何時刻導(dǎo)電巖層中渦旋電流在巷道內(nèi)產(chǎn)生的磁場可以用一個簡單的水平電流環(huán)來等效，該電流環(huán)由發(fā)射線向外擴散的過程像從發(fā)射線圈中“吹”出來的一系列“煙圈”，將此擴散過程稱為“煙圈”效應(yīng)[9]。在礦井巷道有限的空間探測范圍內(nèi)，利用接收線圈接收的二次感應(yīng)場反映發(fā)射線圈周圍所有地質(zhì)體的導(dǎo)電性分布，因此視電阻率計算公式可表示為：

式中：[C]為全空間響應(yīng)系數(shù)；[t]為二次場衰減時間；[N]為線圈匝數(shù)；[S]為接收線圈的面積；[VI]為接收線圈歸一化二次場電位值。

2 閾值分割

閾值分割是三維重建中非常重要的環(huán)節(jié)，能夠有效分割存在灰度差異的圖像，簡化了后續(xù)的分析和處理[10]。為了得到圖像的指定區(qū)域，把圖像看作灰度值不同的目標(biāo)區(qū)與背景區(qū)，目標(biāo)不同，相應(yīng)的灰度值也會不一樣，為此可以選取一個合適的閾值，以閾值為分界線，劃分圖像中的各個像素點為背景區(qū)域和目標(biāo)區(qū)域，從而突出顯示需要重建的目標(biāo)，減少非目標(biāo)區(qū)域?qū)χ亟ㄐЧ挠绊?，提高重建效率?/p>

假定閾值[t，]Otsu最大類間方法利用閾值將圖像分為前景與背景[11]，計算前景與背景圖像的方差，當(dāng)方差達到最大值時，[t]為最佳閾值，圖像中的像素灰度級高于[t]時，保持原灰度級不變，其他像素灰度級都設(shè)置為0，其原理可表示為：

式中：[G（x，y）]為閾值分割后的灰度值；[F（x，y）]為原始灰度值。

圖1為原始二維視電阻率斷面等值線圖，黑色部分為低阻區(qū)域，將原始圖像經(jīng)過閾值分割后得到圖2所示的目標(biāo)區(qū)域圖，與圖1相比，圖2去除了非異常區(qū)域圖像，保留了需要三維重建的異常區(qū)目標(biāo)圖像，減少了背景區(qū)域?qū)δ繕?biāo)區(qū)域重建的影響。

3 三維重建

三維重建實現(xiàn)了二維視電阻率斷面圖的三維模擬顯示，為了充分利用反演得到的斷面圖，運用體繪制算法實現(xiàn)斷面圖像的三維顯示，體繪制首先按照位置和角度信息對斷面圖的數(shù)據(jù)進行規(guī)則化處理，然后在三維空間中形成一個規(guī)則的數(shù)據(jù)場，體素為規(guī)則數(shù)據(jù)場網(wǎng)格中的節(jié)點，描述對象的密度等屬性信息，數(shù)據(jù)場中的體素可當(dāng)做接收或發(fā)出光線的點，依據(jù)每個體素的屬性不同分配相應(yīng)的透明度與顏色，最后得到具有一定顏色與透明度的三維圖像，圖3為體繪制處理流程。

圖4為低阻異常區(qū)的三維重建結(jié)果圖，圖形呈現(xiàn)半透明的狀態(tài)，圖中[x]坐標(biāo)為水平方向延展距離，[y]坐標(biāo)為垂直軸線方向距離，[z]坐標(biāo)為垂直軸線豎直方向的延展距離，通過設(shè)置view（az，el）函數(shù)中的參數(shù)來旋轉(zhuǎn)三維立體圖，其中az表示方位角，el表示仰角，由此可從不同角度觀察視電阻率的空間分布及形狀。圖4a）～圖4c）分別從前向、側(cè)向及斜上方三個視角對圖像進行觀察，由圖中可看出重建圖像中有三處異常區(qū)域，主要集中在深度80～130 m處，其中一處較大的異常區(qū)域從測點號50一直延伸到測點號120，測點號100左右低阻體范圍突然變小，之后繼續(xù)擴大，另外兩處異常區(qū)相對較小，可以判斷出其中一處較小的含水體在測號10～40，深度為20～35 m，低阻體范圍隨著測點號的增大逐漸變寬，另一含水體分布在測號0～20，深度為100～135 m。

綜合前向、側(cè)向與斜上方三個方位的觀察，深度35～80 m為探測安全區(qū)，沒有低阻異常體，在測點50～140，采掘面到深度80 m沒有低阻體出現(xiàn)，結(jié)合探測結(jié)果可以得出，重建得到的三維圖像符合人們的視覺習(xí)慣，與傳統(tǒng)的二維等值線圖解釋相比，能夠清晰地顯示垂向電阻率變化特征。

4 結(jié) 論

三維重建技術(shù)中的閾值分割有效地將等值線斷面圖中的異常區(qū)與其他阻值區(qū)域區(qū)分，減少了數(shù)據(jù)處理量，與以往的二維視電阻率斷面圖等值線顯示相比，對異常體進行三維成像的結(jié)果顯示更加直觀，信息量也更豐富，與真實地質(zhì)體電性分布情況吻合較好，對三維重建結(jié)果圖進行旋轉(zhuǎn)，可實現(xiàn)從不同方向?qū)Ξ惓ｓw空間位置及分布范圍的觀察，對防治礦井水災(zāi)害具有重要的意義，方便非物探人員直觀了解瞬變電磁法探測結(jié)果。

參考文獻

[1] 程久龍，李飛，彭蘇萍，等.礦井巷道地球物理方法超前探測研究進展與展望[J].煤炭學(xué)報，2014，39（8）：1742?1750.

CHENG Jiulong， LI Fei， PENG Suping， et al. Research progress and development direction on advanced detection in mine roadway working face using geophysical methods [J]. Journal of China coal society， 2014， 39（8）： 1742?1750.

[2] 于景邨，胡兵，劉振慶，等.礦井瞬變電磁探測技術(shù)的應(yīng)用[J].物探與化探，2011，35（4）：532?535.

YU Jingcun， HU Bing， LIU Zhenqing， et al. Application of mine transient electromagnetic detection technology [J]. Geophysical and geochemical exploration， 2011， 35（4）： 532?535.

[3] 李富，廖國忠，劉鑫明.礦井瞬變電磁法在煤礦陷落區(qū)探測中的應(yīng)用[J].煤礦安全，2013，44（4）：155?158.

LI Fu， LIAO Guozhong， LIU Xinming. Application of mine transient electromagnetic method in coal mine collage area detection [J]. Safety in coal mines， 2013， 44（4）： 155?158.

[4] 王東偉，劉志新，武俊文，等.礦井瞬變電磁法在巷道迎頭超前探測中的應(yīng)用[J].工程地球物理學(xué)報，2011，8（4）：403?407.

WANG Dongwei， LIU Zhixin， WU Junwen， et al. Application of mine transient electromagnetic method to advanced detecting of roadway head [J]. Chinese journal of engineering geophysics， 2011， 8（4）： 403?407.

[5] 張衛(wèi)，林承灝，于仲，等.工作面頂板富水區(qū)的礦井瞬變電磁法探測[J].煤炭技術(shù)，2012，31（1）：164?166.

ZHANG Wei， LIN Chenghao， YU Zhong， et al. Detection of rich water areas about coal seam roof by mine transient electromagnetic method [J]. Coal technology， 2012， 31（1）： 164?166.

[6] 蔣宗霖，于景邨，孫偉濤.礦井瞬變電磁法低阻體的全空間響應(yīng)影響研究[J].煤炭科學(xué)技術(shù)，2012，40（8）：107?110.

JIANG Zonglin， YU Jingcun， SUN Weitao. Study on mine transient electromagnetic method affected to full space response of low resistance body [J]. Coal science and technology， 2012， 40（8）： 107?110.

[7] 蔡海洋，付俊輝，呂精巧.醫(yī)學(xué)圖像三維拼接重建優(yōu)化技術(shù)研究與仿真[J].計算機仿真，2012，29（7）：308?310.

CAI Haiyang， FU Junhui， L? Jingqiao. Research and simulation of optimization technology of 3D medical image mosaics and reconstruction [J]. Computer simulation， 2012， 29（7）： 308?310.

[8] 張愛東，李炬，陳發(fā)，等.工業(yè)CT圖像的三維重建[J].核電子學(xué)與探測技術(shù)，2005，25（4）：420?422.

ZHANG Aidong， LI Ju， CHEN Fa， et al. Three?dimensional ICT reconstruction [J]. Nuclear electronics and detection technology， 2005， 25（4）： 420?422.

[9] 張軍，趙瑩，李萍.礦井瞬變電磁法在超前探測中的應(yīng)用研究[J].工程地球物理學(xué)報，2012，9（1）：49?53.

ZHANG Jun， ZHAO Ying， LI Ping. Application research of mine transient electromagnetic method in advanced detection [J]. Chinese journal of engineering geophysics， 2012， 9（1）： 49?53.

[10] 范澤華，白鐵成.基于改進差分演化的高效圖像分割算法[J].現(xiàn)代電子技術(shù)，2016，39（14）：87?91.

FAN Zehua， BAI Teicheng. High efficient image segmentation algorithm based on improved differential evolution [J]. Modern electronics technique， 2016， 39（14）： 87?91.

[11] 劉紫燕，吳俊熊，毛攀，等.基于遺傳模擬退火算法的Otsu圖像分割研究[J].電視技術(shù)，2016，40（8）：15?18

LIU Ziyan， WU Junxiong， MAO Pan， et al. Image segmentation on genetic simulated annealing algorithm [J]. Video engineering， 2016， 40（8）： 15?18.