測線分類在坑透反演成圖中的應用

郭方，李培根，齊順

(1.福建工程學院信息科學與工程學院，福州 350108; 2.福州華虹智能科技開發有限公司，福州 350001)

測線分類在坑透反演成圖中的應用

郭方1，李培根2，齊順2

(1.福建工程學院信息科學與工程學院，福州 350108; 2.福州華虹智能科技開發有限公司，福州 350001)

針對目前井下坑透時干擾因素較多而導致接收場強噪聲大的問題，提出一種新的接收場強濾波方法。探索了測線長度分類在提高坑透反演成圖質量中的應用。首先將工作面內發射－接收距離相同射線的接收場強數據經過分類劃為一組，然后利用小波分析方法對各組測線的接收場強進行降噪處理，經過降噪處理的數據用來進行坑透反演成圖。理論仿真和實際數據分析表明:基于測線分類的濾波方法能在一定程度上改善坑透反演成圖的質量，較好地突出異常區的位置，對坑透反演成像有一定的指導意義。

坑透;測線分類;反演成圖

無線電磁波坑透技術(坑透法)是利用探測目標與周圍介質之間的電性差異來研究、確定目標體位置、形態、大小及物性參數的一種物探方法。當前，坑透法與瞬變電磁法、電法等一起構成煤炭綜采的高安全保障技術手段［1－4］。利用此技術，通過數據反演，可以方便地了解綜采工作面的復雜地質構造，判斷地質內的異常區域［5－10］。在綜采數據反演成像過程中，接收端場強作為工作面異常體的間接表征，其數據質量對成像效果有著重要的影響。但是在實際工作中，由于周邊環境和儀器本身的干擾，接收場強噪聲很大［11］。針對這種情況，一些學者提出相應的去噪方法，總結出一些求取初始場強的方法［12－13］。這些方法通過多次接收求平均來消除儀器等造成的隨機誤差，但井下風管、水管、電纜、金屬導體及動力設備所引起的輕微場強誤差一般很難去除，從而在較大程度上影響坑透反演成圖的質量。

基于此，本文在有關學者已有工作的基礎上，根據綜采面測線分布特點和規律，提出一種新的接收場強濾波方式。由于一個綜采面有很多條發射－接收點距離相同的測線，這些線并不全部經過異常點，一般情況下，只有很少的一部分經過異常點，而那些沒經過異常點的測線場強值就代表了綜采面的背景接收值，理論上這些值應該相同或相近。因此，可以把這些相同長度測線的場強值分為一組，將各組場強值進行濾波處理后，用濾波后的新值代替原接收場強值進行坑透反演。

1 坑透理論

坑透法以無線電磁波發射為基礎。根據電磁波傳播理論，電磁波在均勻介質中的輻射場隨距離呈一定的規律變化［12］。

假設原點A在輻射源(天線軸)中心，在無限均勻、各向同性的介質中，觀測點P到A點的距離為r(觀測點在輻射場內，見圖1)。在此條件下求解導電介質中的波動方程，即可求得P點的電磁場強度H。式(1)中:H0為初始場強，取決于發射機功率及周圍介質;β為煤層對電磁波能量的吸收系數;r為測點到輻射源原點的距離;sinθ為方向因子，θ是天線軸線與觀測方向的夾角，一般取90°。其中，β又可以用以下公式表示:

式(2)中:ω為角頻率;μ為導磁率;ε為介電常數; σ為導電率。

當頻率ω=2πf不變時，β=F(ε、μ、σ)，即β是ε、μ、σ的函數。當電磁波穿越不同電性的煤巖層時，μ差別不大。所以ε、σ的改變會引起β與H值的明顯變化。

圖1 偶極子天線輻射場

當電磁波在地下巖層中傳播時，由于各種巖、礦石電性(電阻率ρ和介電常數ε)不同，對電磁波能量的吸收不同，于是可根據這些不同結構的吸收參數來判斷工作面內的異常區。低阻巖層對電磁波具有較強的吸收作用，當波前進方向遇到斷裂構造所出現的界面時，電磁波將在界面上產生反射和折射作用，造成能量的損耗，使接收巷道中的電磁波信號十分微弱甚至接收不到透射信號，形成所謂的透射異常(又稱陰影異常)，見圖2。研究采區煤層、各種構造及地質體對電磁波的影響所造成的各種無線電波透視異常，從而進行地質推斷和解釋，這就是坑透法的基本原理。

圖2 無線坑透接收場強

為計算方便，實際應用中，式(1)又可以表示為

2 工作面測線分類

2.1 工作面測線分布

無線電波透視觀測一般在兩巷道間進行，其中一巷道布置發射點，向煤層中發射某一頻率的電磁波，接受點布置在另一巷道進行相應接收。為提高探測效率，目前普遍采用定點法方式［13］。發射機在一定的時間內相對固定位置，接收機在一定范圍內逐點觀測其場強值，即定點發射多點等間距(扇形)接收。

在實際情況下，綜采工作面寬度大都在150 m左右，長度為500 m以上。測線采用扇形結構，各接收點的間距在10 m左右，一般接收點為11 m (圖3)。因此各測線長度差別不很大，即理論接收場強值相差很小。

圖3 無線坑透測線分布

2.2 工作區內測線分類

在實際坑透過程中，由于一個工作面往往有多個相同的扇形測區，若每個測區都采用對稱分布，則一個測區就有兩條長度相同的測線。對于n個發射測點的綜采工作面，就有2n個相同的測線。因此可把這2n條相同測線的場強接收值化為一組進行信號處理。這樣，整個工作面測線可分為m組相同長度的測線。對于每組測線，根據式(1)，理論上它們的接收值應該相等，但實際上由于測試誤差、現場干擾和異常區的影響，這些測線的值并不完全相等。考慮到這些測線的長度相同，因此可以對每組測線進行濾波處理。

近年來，創造力研究領域更具包容性地關注創造性的情感．許多研究顯示，積極的情感影響個人的創造過程．什么是積極的創造性情感，以下研究能給出一定的啟示．

2.3 小波濾波去噪

在信號濾波的眾多方法中，由于小波變換具有良好的時頻局部化和多分辨率等特性，因此可很好地刻畫信號邊緣、尖峰、斷點等非平穩特性，在去除噪聲對信號影響的同時，還可以保留信號的突變位置和邊緣信息，具有更好的濾波效果。所以，以小波變換為基礎的濾波方法廣泛用于各種信號處理過程［14－15］，其基本定義如下:

對任一函數

則它的連續小波變換為

由于每組接收場強的坑透數據具有非平穩特性，往往每組數據既有噪聲數據，又有異常點的突變數據，因此較適合于小波處理。

小波降噪的一般處理過程為:①選擇小波包進行小波正變換;②閾值處理;③小波反變換。

3 實驗結果分析

為測試算法效果，本文根據式(3)建立均勻介質中的異常構造幾何模型(如圖4所示)，并取背景吸收系數為0.003 db/m，異常區吸收系數為0.005 db/m。3塊異常區大小都是30 m×20 m，模型區域網格劃分為100×20。

圖4 綜采面仿真模型

取初始場強為120 db進行系統仿真。為反映實際測試環境，仿真后工作區內各接收場強值加入幅值為±0.5 db的環境噪聲，其場強曲線如圖5所示。

從圖5可以看出:仿真工作面共有200 m、200.25 m、201.00 m、202.24 m、203.96 m、206.16 m 6種測線，因此整個工作面測線被分為6組。各組測線和小波濾波后的場強值曲線如圖6所示。

圖5 仿真場強分布

圖6 各組仿真場強值和濾波后的場強值

將原始場強值和濾波后的場強值分別代入式(3)進行工作面成圖反演，迭代500次得到圖7的反演成圖。對比圖4的仿真模型可以看出，濾波后的異常點2和4更加突出。

圖7 仿真數據反演成圖

3.2 綜采數據分析

為測試實際算法效果，本文對某礦采集的數據進行分析。礦體綜采面長度為1 200 m，寬度為185 m。按照圖3測線分布測得工作面的場強分布曲線如圖8所示。

圖8 實測場強分布

在本工作面測量中，共有185.00 m、185.27 m、186.08 m、187.42 m、189.28 m、191.64 m共6種測線。整個工作面測線被分為6組，各組測線和小波濾波后的場強值曲線如圖9所示。

圖9 各組測線場強值和濾波后的場強值

將濾波后的場強接收值代入式(3)，利用SIRT算法求解此方程組，迭代500次后得到工作面反演圖如圖10所示。從圖10可看出，工作面共有4處異常區。對比圖10(a)和(b)可知:使用小波濾波后的反演成圖中，異常區2、3、4的位置信息更為具體。現場實際回采發現，工作面內的4個異常區的位置與圖10(b)成圖結果較為接近。

圖10 不同接收場強反演成圖

4 結束語

根據坑透工作面數據分布特點，提出了一種新的數據處理方法。該方法將工作面內長度相同的測線接收場強值歸為一組，在此基礎上對每組數據進行小波降噪處理。濾波后的數據被用來進行坑透反演成圖。理論和實踐證明:將測線進行歸類濾波可在一定程度上壓制環境背景噪聲，提高反演成圖的清晰度和分辨率，從而與實際更為接近，對現場應用和進一步的算法研究具有一定的指導意義。

［1］寧書年，張紹紅.無線電波層析成像技術及在礦井坑透中的應用［J］.煤炭學報，2001，26(5):468－472.

［2］郭文波，李貅，薛國強，等.瞬變電磁快速成像解釋系統研究［J］.地球物理學報，2005 48(6):1400－1405.

［3］劉樹才，劉鑫明，姜志海，等.煤層底板導水裂隙演化規律的電法探測研究［J］.巖石力學與工程學報，2009，28(2):348－356.

［4］韓連生，范廣海，賈建軍，等.地震映像法在鐵礦隱伏空區探測中的應用研究［J］.工程地球物理學報，2012，9(4):469－472.

［5］湯友誼，陳江峰，彭立世.無線電波坑道透視構造煤的研究［J］.煤炭學報，2002，27(3):254－258.

［6］劉煥新，王利宏，劉樹才.坑透高精度CT層析成像測量方法試驗研究［J］.山東煤炭科技，2010，(1):103－105.

［7］劉加利.坑透技術在復雜地質條件下的應用［J］.河北煤炭，2005，(4):39－40.

［8］孫四清，湯友誼，陳江峰.工作面坑透資料在瓦斯地質預測中應用的研究［J］.煤礦安全，2005，36(5):1－3.

［9］王均雙，薄夫利，馬沖.坑透CT成像技術在工作面地質構造探測中的應用［J］.煤炭科學技術，2008，36 (10):93－96.

［10］侯恩科，郝珠成.坑透層析成像技術在金礦定位預測中的應用［J］.地質與勘探，2001，37(4):43－45.

［11］冀學宏.論影響無線電波坑透的干擾因素及排除方法［J］.西山科技，2001，6:18－22.

［12］吳燕清.地下電磁波探測及應用研究［D］.長沙:中南大學，2002.

［13］肖玉林.煤礦綜采工作面無線電波透視技術研究［D］.淮南:安徽理工大學，2010.

［14］蔡強富，陳鵬，李固，等.電磁超聲回波信號的改進閾值函數降噪法研究［J］.國外電子測量技術，2011，30 (12):48－51.

［15］李越囡，孫迎.聽覺誘發電位信號的小波消噪方法研究［J］.儀器儀表學報，2010，31(3):541－545.

(責任編輯 楊黎麗)

Application of Classified Probing Lines in Inverting Image with Tunnel Perspective

GUO Fang1，LI Pei－gen2，QI Shun2
(1.Department of Information Science and Engineering，Fujian University of Technology，Fuzhou 350108，China;
2.Fuzhou Huahong Intelligent Technology Development Co.，Ltd，Fuzhou 350001，China)

Aiming for the issue of receiving field intensity with much noise，a novel method to filter the noise of receiving field intensity is proposed，and the technology of enhancing the imaging quality of radio electromagnetic wave tunnel perspective with classified probing lines is explored initially.Firstly，the

field intensity values of the same probing lines are classified as a group.Then the each group data is filtered with wavelet.Finally，the optimized data is used to invert the image of working plane.The attenuation formula of electromagnetic are established，and the background absorption coefficient and initial field intensity is solved by LM optimization algorithm.Theoretical analysis and experimental results indicate that the method could prompt the imaging quality to some extent，and the location of abnormal area could be enhanced.The proposed method is significant to guide the inverting and imaging of the working plane in the process of radio electromagnetic wave tunnel perspective.

tunnel perspective;probing lines classification;inverting image

TD15

A

1674－8425(2014)03－0122－05

10.3969/j.issn.1674－8425(z).2014.03.022

2013－11－28

福建省教育廳資助項目(JB13142);福建工程學院科研啟動基金資助項目(GY－Z12078)

郭方(1977—)，男，安徽人，講師，主要從事光電信息處理方面研究。

郭方，李培根，齊順.測線分類在坑透反演成圖中的應用［J］.重慶理工大學學報:自然科學版，2014(3): 122－126.

format:GUO Fang，LI Pei－gen，QI Shun.Application of Classified Probing Lines in Inverting Image with Tunnel Perspective［J］.Journal of Chongqing University of Technology:Natural Science，2014(3):122－126.