王國富 龍琴 王小紅 張法全 葉金才
摘 要: 針對(duì)傳統(tǒng)的礦井瞬變電磁探測(cè)大多以二維視電阻率斷面等值線圖顯示,難以對(duì)低阻異常的空間位置及分布范圍做出準(zhǔn)確的判斷,提出一種二維視電阻率斷面等值線圖的三維重建方法。首先利用閾值分割原理將電阻率斷面圖的低阻異常區(qū)與其他區(qū)域分開,然后通過體繪制方法實(shí)現(xiàn)低阻異常區(qū)域的三維圖像顯示。實(shí)驗(yàn)結(jié)果表明,三維重建能夠有效顯示垂向視電阻率的變化特征,提高了對(duì)低阻異常體空間位置的判定。
關(guān)鍵詞: 視電阻率斷面圖; 礦井瞬變電磁探測(cè); 閾值分割; 三維重建; 等值線; 體繪制
中圖分類號(hào): TN911.7?34 文獻(xiàn)標(biāo)識(shí)碼: A 文章編號(hào): 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)濟(jì)損失[1],其中含水低阻結(jié)構(gòu)體的存在是造成礦井災(zāi)害的主要因素。礦井瞬變電磁法[2]因其對(duì)低阻體敏感,具有較好的探測(cè)方向,易攜帶等特性而被廣泛地用于探測(cè)礦井中的含水結(jié)構(gòu)體[3?5]。然而,傳統(tǒng)的礦井瞬變電磁探測(cè)結(jié)果大部分以二維視電阻率等值線斷面圖的形式解釋,二維等值線圖能夠較清楚地顯示橫向斷面的阻值特征,但無法清晰地顯示垂向延展面上的視電阻率特征分布情況[6],難以準(zhǔn)確地判斷出異常體的空間位置及分布范圍,特別是垂直方向上的異常體變化特征,因此,三維立體的顯示出低阻體不僅有利于確定礦井中含水體結(jié)構(gòu)的垂直特征,而且有利于物探人員的觀察。
三維重建技術(shù)在醫(yī)學(xué)及工業(yè)領(lǐng)域已得到廣泛應(yīng)用,醫(yī)學(xué)中常用三維重建技術(shù)對(duì)人體掃描得到的CT圖像進(jìn)行處理[7],由此可以立體地觀察人體器官、骨骼、血管等,可以幫助醫(yī)務(wù)人員準(zhǔn)確判斷病情,快速找出病因,工業(yè)中常用三維重建獲取模型的立體圖形[8],取得了較好的效果,在此提出將三維重建技術(shù)應(yīng)用到礦井安全防護(hù)領(lǐng)域,結(jié)合實(shí)際地質(zhì)資料及探測(cè)數(shù)據(jù),對(duì)礦井瞬變電磁探測(cè)數(shù)據(jù)進(jìn)行預(yù)處理,采用三維重建技術(shù)解釋礦井瞬變電磁探測(cè)成果,實(shí)現(xiàn)垂向視電阻率變化的可視化,使探測(cè)結(jié)果顯示更直觀,為探測(cè)含水構(gòu)造提供了依據(jù)。
1 礦井瞬變電磁原理
礦井瞬變電磁法在井下巷道空間鋪設(shè)發(fā)射線圈并給線圈通以一定的電流,由此發(fā)射回線的周圍空間會(huì)形成一個(gè)穩(wěn)定的磁場(chǎng)。在某一時(shí)刻斷開電流,空間磁場(chǎng)由于一次磁場(chǎng)在巖層中激發(fā)的感應(yīng)電流得以維持而不會(huì)立即消失,在一次電流斷開時(shí),處于發(fā)射線圈周圍空間中的礦井巷道頂、底板和側(cè)幫最先感應(yīng)到電流斷開瞬間一次磁場(chǎng)所產(chǎn)生的劇烈變化,因此,由一次磁場(chǎng)激發(fā)的感應(yīng)電流會(huì)局限在巷道附近巖層中,而且?guī)r層中各處的感應(yīng)電流分布不均,隨著時(shí)間的推移,巖層中感應(yīng)電流強(qiáng)度變?nèi)醪⒙呌诰鶆颉?/p>
研究結(jié)果表明,任何時(shí)刻導(dǎo)電巖層中渦旋電流在巷道內(nèi)產(chǎn)生的磁場(chǎng)可以用一個(gè)簡(jiǎn)單的水平電流環(huán)來等效,該電流環(huán)由發(fā)射線向外擴(kuò)散的過程像從發(fā)射線圈中“吹”出來的一系列“煙圈”,將此擴(kuò)散過程稱為“煙圈”效應(yīng)[9]。在礦井巷道有限的空間探測(cè)范圍內(nèi),利用接收線圈接收的二次感應(yīng)場(chǎng)反映發(fā)射線圈周圍所有地質(zhì)體的導(dǎo)電性分布,因此視電阻率計(jì)算公式可表示為:
式中:[C]為全空間響應(yīng)系數(shù);[t]為二次場(chǎng)衰減時(shí)間;[N]為線圈匝數(shù);[S]為接收線圈的面積;[VI]為接收線圈歸一化二次場(chǎng)電位值。
2 閾值分割
閾值分割是三維重建中非常重要的環(huán)節(jié),能夠有效分割存在灰度差異的圖像,簡(jiǎn)化了后續(xù)的分析和處理[10]。為了得到圖像的指定區(qū)域,把圖像看作灰度值不同的目標(biāo)區(qū)與背景區(qū),目標(biāo)不同,相應(yīng)的灰度值也會(huì)不一樣,為此可以選取一個(gè)合適的閾值,以閾值為分界線,劃分圖像中的各個(gè)像素點(diǎn)為背景區(qū)域和目標(biāo)區(qū)域,從而突出顯示需要重建的目標(biāo),減少非目標(biāo)區(qū)域?qū)χ亟ㄐЧ挠绊?,提高重建效率?/p>
假定閾值[t,]Otsu最大類間方法利用閾值將圖像分為前景與背景[11],計(jì)算前景與背景圖像的方差,當(dāng)方差達(dá)到最大值時(shí),[t]為最佳閾值,圖像中的像素灰度級(jí)高于[t]時(shí),保持原灰度級(jí)不變,其他像素灰度級(jí)都設(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 三維重建
三維重建實(shí)現(xiàn)了二維視電阻率斷面圖的三維模擬顯示,為了充分利用反演得到的斷面圖,運(yùn)用體繪制算法實(shí)現(xiàn)斷面圖像的三維顯示,體繪制首先按照位置和角度信息對(duì)斷面圖的數(shù)據(jù)進(jìn)行規(guī)則化處理,然后在三維空間中形成一個(gè)規(guī)則的數(shù)據(jù)場(chǎng),體素為規(guī)則數(shù)據(jù)場(chǎng)網(wǎng)格中的節(jié)點(diǎn),描述對(duì)象的密度等屬性信息,數(shù)據(jù)場(chǎng)中的體素可當(dāng)做接收或發(fā)出光線的點(diǎn),依據(jù)每個(gè)體素的屬性不同分配相應(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è)向及斜上方三個(gè)視角對(duì)圖像進(jìn)行觀察,由圖中可看出重建圖像中有三處異常區(qū)域,主要集中在深度80~130 m處,其中一處較大的異常區(qū)域從測(cè)點(diǎn)號(hào)50一直延伸到測(cè)點(diǎn)號(hào)120,測(cè)點(diǎn)號(hào)100左右低阻體范圍突然變小,之后繼續(xù)擴(kuò)大,另外兩處異常區(qū)相對(duì)較小,可以判斷出其中一處較小的含水體在測(cè)號(hào)10~40,深度為20~35 m,低阻體范圍隨著測(cè)點(diǎn)號(hào)的增大逐漸變寬,另一含水體分布在測(cè)號(hào)0~20,深度為100~135 m。
綜合前向、側(cè)向與斜上方三個(gè)方位的觀察,深度35~80 m為探測(cè)安全區(qū),沒有低阻異常體,在測(cè)點(diǎn)50~140,采掘面到深度80 m沒有低阻體出現(xiàn),結(jié)合探測(cè)結(jié)果可以得出,重建得到的三維圖像符合人們的視覺習(xí)慣,與傳統(tǒng)的二維等值線圖解釋相比,能夠清晰地顯示垂向電阻率變化特征。
4 結(jié) 論
三維重建技術(shù)中的閾值分割有效地將等值線斷面圖中的異常區(qū)與其他阻值區(qū)域區(qū)分,減少了數(shù)據(jù)處理量,與以往的二維視電阻率斷面圖等值線顯示相比,對(duì)異常體進(jìn)行三維成像的結(jié)果顯示更加直觀,信息量也更豐富,與真實(shí)地質(zhì)體電性分布情況吻合較好,對(duì)三維重建結(jié)果圖進(jìn)行旋轉(zhuǎn),可實(shí)現(xiàn)從不同方向?qū)Ξ惓sw空間位置及分布范圍的觀察,對(duì)防治礦井水災(zāi)害具有重要的意義,方便非物探人員直觀了解瞬變電磁法探測(cè)結(jié)果。
參考文獻(xiàn)
[1] 程久龍,李飛,彭蘇萍,等.礦井巷道地球物理方法超前探測(cè)研究進(jìn)展與展望[J].煤炭學(xué)報(bào),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] 于景邨,胡兵,劉振慶,等.礦井瞬變電磁探測(cè)技術(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ū)探測(cè)中的應(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] 王東偉,劉志新,武俊文,等.礦井瞬變電磁法在巷道迎頭超前探測(cè)中的應(yīng)用[J].工程地球物理學(xué)報(bào),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ū)的礦井瞬變電磁法探測(cè)[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].計(jì)算機(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é)與探測(cè)技術(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] 張軍,趙瑩,李萍.礦井瞬變電磁法在超前探測(cè)中的應(yīng)用研究[J].工程地球物理學(xué)報(bào),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ìn)差分演化的高效圖像分割算法[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.