礦井瞬變電磁法超前探測合成孔徑成像研究

程久龍， 陳丁， 薛國強， 邱浩， 周項通

1 中國礦業(yè)大學(xué)(北京)煤炭資源與安全開采國家重點實驗室， 北京　100083 2 中國礦業(yè)大學(xué)(北京) 地球科學(xué)與測繪工程學(xué)院， 北京　100083 3 中國科學(xué)院礦產(chǎn)資源研究重點實驗室，中國科學(xué)院地質(zhì)與地球物理研究所， 北京　100029 4 煤炭科學(xué)研究總院礦山安全技術(shù)研究分院， 北京　100013

?

礦井瞬變電磁法超前探測合成孔徑成像研究

程久龍1,2， 陳丁2， 薛國強3， 邱浩4， 周項通2

1 中國礦業(yè)大學(xué)(北京)煤炭資源與安全開采國家重點實驗室， 北京100083 2 中國礦業(yè)大學(xué)(北京) 地球科學(xué)與測繪工程學(xué)院， 北京100083 3 中國科學(xué)院礦產(chǎn)資源研究重點實驗室，中國科學(xué)院地質(zhì)與地球物理研究所， 北京100029 4 煤炭科學(xué)研究總院礦山安全技術(shù)研究分院， 北京100013

摘要在礦井瞬變電磁超前探測數(shù)據(jù)處理和解釋中，為了更好地識別巷道前方含水體電性分界面，利用擴散場與波動場間的函數(shù)關(guān)系，將礦井瞬變電磁法超前探測中瞬變電磁波場轉(zhuǎn)換為擬地震波場，實現(xiàn)波場轉(zhuǎn)換.同時，為了解決波場轉(zhuǎn)換后信號幅值偏弱，以增加信噪比利于解釋，利用雷達的合成孔徑方法對波場轉(zhuǎn)換后的擬地震信號進行相關(guān)疊加處理.將礦井瞬變電磁法巷道超前探測的同一測點不同探測方向的數(shù)據(jù)作為孔徑數(shù)據(jù)進行合成孔徑成像，以突出弱異常而提高信噪比，結(jié)合應(yīng)用實例驗證了方法的效果.研究表明，礦井瞬變電磁法合成孔徑成像方法能夠提高電性界面信息的信噪比和界面幾何分辨率，對礦井瞬變電磁法精細超前探測含水構(gòu)造具有重要的理論意義和應(yīng)用價值.關(guān)鍵詞礦井瞬變電磁法； 波場轉(zhuǎn)換； 擬地震信號； 合成孔徑成像； 電性界面

1引言

近幾年來，礦井瞬變電磁法在礦井工作面頂板、底板巖層富水性探測和掘進巷道超前探測中發(fā)揮了重要作用(程久龍等, 2014).由于礦井瞬變電磁法勘探是體積勘探，受體積效應(yīng)的制約，其分辨率難以滿足生產(chǎn)實際要求(Cheng et al.,2015)，并且受井下施工環(huán)境和金屬物體的影響，采集到的數(shù)據(jù)中往往存在一定的噪聲，致使信噪比降低，而且，礦井瞬變電磁法探測成果大多是以電阻率等值線斷面圖的形式表達，不能準確地分辨實際地質(zhì)情況中最為重要的地層電性界面.基于波場變換的瞬變電磁擬地震成像技術(shù)因能夠相對突出異常體的物性和幾何參數(shù)而越來越受到人們的重視(李貅和薛國強, 2013)，然而由于巷道超前探測方向不可能很密，加上采樣點相對稀疏等問題，所以波場變換成像的信噪比和分辨率有待進一步提高，合成孔徑雷達及航空瞬變電磁中合成孔徑成像的思想為礦井瞬變電磁法合成孔徑成像提供了新的思路.

合成孔徑成像技術(shù)最早是應(yīng)用于雷達方面，國內(nèi)外學(xué)者做了一系列的研究，研究內(nèi)容主要包括算法、應(yīng)用及成像技術(shù)的對比等(Cetin and Karl, 2001;王治樂等, 2004；Samadi et al.,2009)，應(yīng)用領(lǐng)域較為廣泛.在超聲波探測方面，Jensen等(2006)利用合成孔徑成像技術(shù)來提高醫(yī)學(xué)超聲波檢測圖的成像分辨率，Ganguli等(2012)研究了混凝土超聲波探傷的合成孔徑聚焦成像技術(shù).在電磁法探測方面，Mason等(2001)研究了利用寬頻井中雷達合成孔徑干涉法來繪制3D礦體圖，張春城和周正歐(2004)研究了基于stolt偏移的探地雷達合成孔徑成像.在全空間瞬變電磁數(shù)值模擬方面，李展輝和黃清華(2014)研究了將復(fù)頻率參數(shù)完全匹配層吸收邊界應(yīng)用到均勻全空間模型和包含巷道及異常體的全空間模型.上述研究沒有涉及異常體界面的定位.

在瞬變電磁勘探方面，國內(nèi)學(xué)者鑒于合成孔徑雷達和航空瞬變電磁之間的相似性首先把合成孔徑成像的思想應(yīng)用于航空瞬變電磁，研究了利用合成孔徑方法對航空瞬變電磁數(shù)據(jù)進行成像，以提高分辨率，增加勘探深度(李貅等, 2010).之所以可以把合成孔徑雷達的思想引入到地面瞬變電磁數(shù)據(jù)處理中，是因為在不同測點所接收到的由同一個發(fā)射源所激勵的二次場數(shù)據(jù)等效于合成孔徑雷達單點發(fā)射源時的陣列數(shù)據(jù)，而地面瞬變電磁所有的發(fā)射源則構(gòu)成了動源觀測.由于井下瞬變電磁法探測裝置及數(shù)據(jù)采集方式與地面基本上相似，所以，像地面瞬變電磁合成孔徑成像方法一樣，對井下瞬變電磁數(shù)據(jù)做類似的處理是可行的.本文利用合成孔徑成像技術(shù)對礦井瞬變電磁波場變換后的擬地震信號進行處理，結(jié)合超前探測的工程應(yīng)用實例對方法效果進行對比分析和驗證，實現(xiàn)突出弱異常且提高界面分辨率的目的.

2 波場變換原理

瞬變電磁場與其相對應(yīng)的虛擬波場之間的數(shù)學(xué)積分形式為(Lee et al., 1989)

(1)

將公式(1)改寫為數(shù)值積分形式為(李貅和薛國強, 2013)

(i=1,2,…,m,j=1,2,…,n)，

(2)

其中m為采樣個數(shù)，n為數(shù)值積分時積分區(qū)間剖分的個數(shù)，hj為積分步長.

將公式(2)寫成矩陣形式為

(3)

需要注意的是方程組(3)中的系數(shù)矩陣A常呈病態(tài).為了獲得穩(wěn)定的解，可采用阻尼法最小二乘法(楊文采, 1997)，建立目標函數(shù)為

(4)

其中ε2為阻尼因子，令公式(4)的梯度為0可得

(5)

進而有

(6)

其中(ATA+ε2I)+為(ATA+ε2I)的Moore-Penrose(M-P)廣義逆.對于廣義逆的求解，可采用IMqrginv算法(Ataei, 2014).阻尼因子需在解估計的分辨率和方差間折衷并通過實驗選取(楊文采, 1997; 陳本池等, 1999).礦井瞬變電磁超前探測數(shù)據(jù)是在同一掘進工作面處不同方向進行采集的，在數(shù)據(jù)采集段內(nèi)，誤差水平是一定的，所以對同一個掘進工作面處的超前探測數(shù)據(jù)，其波場變換可采用同一個阻尼因子.

波場變換時，瞬變電磁場中的誤差會給變換后的波場造成較大的影響，特別是當(dāng)解估計的單位協(xié)方差矩陣(楊文采, 1997)較大時，這種影響會更大.為此，在波場變換前需對所采集到的數(shù)據(jù)進行濾波處理，根據(jù)實驗選取小波變換法(Toetal., 2009)進行濾波效果較好.

由于礦井瞬變電磁超前探測位置是在礦井巷道掘進工作面處，故采集的信號是前后雙向回波，可借助于半空間和全空間瞬變響應(yīng)關(guān)系將全空間瞬變響應(yīng)校正為半空間瞬變響應(yīng)，也即相當(dāng)于僅有超前探測時的前方回波信號.半空間和全空間瞬變響應(yīng)磁場間的關(guān)系(KaufmanandEaton, 2001)為

(7)

3礦井瞬變電磁合成孔徑成像方法

3.1雷達和地面瞬變電磁合成孔徑成像原理

合成孔徑雷達(SAR)就是利用雷達-目標的相對運動，以單元陣在不同相對空間位置上接收的時間采樣信號替代大尺寸源的空間采樣信號.由雷達-目標相對運動形成的軌跡構(gòu)成一個合成孔徑來代替龐大的真實孔徑，把尺寸較小的真實天線孔徑用數(shù)據(jù)處理的方法合成一較大的等效天線孔徑的雷達.合成孔徑雷達的特點是分辨率高，能有效地識別目標物(Chan and Koo, 2008)，合成孔徑成像原理如圖1a所示.

地面瞬變電磁法是在地表布設(shè)一個矩形回線發(fā)射框，在發(fā)射框的內(nèi)部和外部接收瞬變電磁信號.地面瞬變電磁法也是在移動方式中觀測瞬變電磁信號，并且不同發(fā)射源間的信號滿足類相干性.因此地面瞬變電磁法也可以使用合成孔徑成像技術(shù)(Guo et al., 2012；李貅等, 2012).地面瞬變電磁合成孔徑原理如圖1b所示.

地面瞬變電磁合成孔徑成像公式為(李貅等, 2012)

(j=1,2,…,m)，

(8)

圖1　雷達合成孔徑和地面電磁法合成孔徑(a) 雷達合成孔徑；(b)地面瞬變電磁合成孔徑.Fig.1　Schematic diagram of SAI of SAR (a) and ground TEM (b)

從以上雷達和地面瞬變電磁合成孔徑成像技術(shù)中可以看出，合成孔徑成像方法的本質(zhì)是不同測點信號間的相關(guān)疊加.

3.2礦井瞬變電磁超前探測合成孔徑成像原理和方法

礦井瞬變電磁法超前探測的地點是在巷道掘進工作面處，如圖2a所示，圖中矩形小回線源緊靠掘進工作面掌子面上，虛線圓圈代表掘進工作面前方瞬變電磁一次場不同時刻的“煙圈”.常采用扇形布置，即在掘進工作面位置放置小線框天線，以巷道正前方為中心軸，按一定的角度(10°或15°)往兩幫分別進行不同方向探測，最終形成扇形探測斷面，如圖2b所示，圖中箭頭表示同一斷面上不同的探測方向(發(fā)射線框的法線方向)，虛線框表示掘進工作面前方的異常體，設(shè)掘進工作面正前方為z軸，z軸所指方向為0°，用與z軸所夾的角度θ來定義探測方向(z軸左側(cè)的角度為負，右側(cè)的角度為正)，即：U(-θi)表示z軸左側(cè)與z軸成θi角度方向上的實測信號值，U(+θi)表示z軸右側(cè)與z軸成θi角度方向上的實測信號值，U(θ0)表示z軸方向上的實測信號值.

從圖2礦井瞬變電磁超前探測的觀測方式中可以看出：同一個掘進工作面處其周圍全空間范圍內(nèi)的地質(zhì)條件在觀測期間是一定的，不同探測方向(也即不同發(fā)射源，相當(dāng)于地面瞬變電磁法的不同測點)觀測到的瞬變電磁信號間具有類相干性.因此，礦井瞬變電磁數(shù)據(jù)具備相干疊加的條件，也即可以采用合成孔徑成像技術(shù).

在礦井瞬變電磁擬地震波場變換、反褶積和時深轉(zhuǎn)換等的基礎(chǔ)上(程久龍等, 2013)，根據(jù)合成孔徑成像的思想，本文提出合成孔徑成像公式為

(9)

其中U′(θm+N，tk)為合成孔徑成像后θm+N探測方向上的信號值，θm是指合成孔徑成像的起始探測方向，而θm+2N為合成孔徑成像的結(jié)束探測方向，m和m+2N分別為起始探測方向編號和結(jié)束探測方向編號，2N+1為合成孔徑的長度，M是指同一探測方向上重復(fù)測量次數(shù)，Wij是權(quán)重函數(shù)，U(θi，tj)是礦井瞬變電磁波場變換后θi探測方向上第j次測量的擬地震波振幅值，τijk為權(quán)函數(shù)Wij取極值時的時移量.

權(quán)函數(shù)Wij為相關(guān)系數(shù)，表達式為

(10)

圖2　巷道瞬變電磁超前探測(a) 巷道瞬變電磁超前探測立體； (b) 巷道瞬變電磁超前探測平面.Fig.2　Schematic diagram of roadway advanced detection using MTEM(a) 3D view； (b) Planar view.

其中i依次取m、m+1、…、m+2N，j取1、2、…M，τijk為公式(10)取極大值時的相對時移量，K為波場變換后某一探測方向上的類時間門數(shù).礦井瞬變電磁超前探測合成孔徑成像原理如圖3所示.

相關(guān)系數(shù)表征了兩個探測方向上波場變換后信號值間的相關(guān)程度.相關(guān)系數(shù)越接近于1或-1，相關(guān)性越強；反之，相關(guān)系數(shù)越接近于0，相關(guān)性越弱.在通常情況下，通過以下取值范圍判斷信號間的相關(guān)性(表1).

表1　相關(guān)系數(shù)與相關(guān)性對照表

圖3　礦井瞬變電磁超前探測合成孔徑原理(a) 方向-θi上的單次探測； (b) 方向θ0上的單次探測； (c) 方向θi上的單次探測.Fig.3　Principle ofsynthetic aperture schematicof advanced detection in roadway(a) Detection in the direction of -θi; (b) Detection in the direction of θ0; (c) Detection in the direction of θi.

4工程應(yīng)用

某礦巷道掘進至G22+5 m位置，需要查清掘進迎頭前方是否存在含水地質(zhì)體，對迎頭前方100 m范圍內(nèi)砂巖富水情況、斷層帶富水情況進行總體控制，對含水異常區(qū)進行地質(zhì)解釋.選擇礦井瞬變電磁法進行超前探測，裝置方式為中心回線，回線邊長為2 m×2 m，電流1.8 A，關(guān)斷時間約為110 μs，實測數(shù)據(jù)的時間區(qū)間為[116.8 μs, 7088 μs]，均在關(guān)斷時間之外.巷道圍巖為砂巖，電阻率大于20 Ωm.

圖4是G22+5 m探測點水平斷面11個方向(以巷道掘進方向為0°，往左幫每隔10°一個探測方向，到左幫50°，往右?guī)兔扛?0°一個探測方向，到右?guī)?0°)的實測感應(yīng)電動勢曲線圖，圖中曲線表示不同探測方向上相同時間門間的連線.

圖4　G22+5 m位置處TEM超前探測水平斷面實測感應(yīng)電動勢Fig.4　Measured induction electromotive force on horizontal sectionby MTEM at G22+5 m

對圖4實測數(shù)據(jù)分別進行濾波和全空間校正后，進行波場轉(zhuǎn)換和時深轉(zhuǎn)換等數(shù)據(jù)處理(程久龍等, 2013)，可以得到該測點水平斷面上11個探測方向(從左幫-50°至右?guī)?50°，每隔10°依次探測)的擬地震信號，橫軸表示深度，縱軸代表探測方向角，如圖5所示.從圖中可以看出，擬地震信號有明顯的兩個同相軸，但總體偏弱.第一個相位在30～40 m范圍內(nèi)，能量相對較強，反映的是高阻到低阻的界面特征；第二個相位在65～80 m范圍內(nèi)，能量相對較弱，反映的是低阻到高阻的界面特征(程久龍等, 2013).

圖5　G22+5 m位置處TEM超前探測水平斷面波場變換后擬地震信號Fig.5　Pseudo-seismic wave-field signalafter transformon horizontal section by MTEM at G22+5 m

根據(jù)公式(10)，對水平斷面11個探測方向的擬地震數(shù)據(jù)分別計算相鄰兩探測方向之間的相關(guān)系數(shù)，判斷相關(guān)性，見表2.從表中可以看出，水平斷面上相鄰兩測點數(shù)據(jù)之間最小的相關(guān)系數(shù)為0.12，最大的相關(guān)系數(shù)為0.9726.

孔徑長度的選擇，即中心探測方向確定以后，具

體應(yīng)該選取幾個探測方向上的數(shù)據(jù)進行疊加而且疊加后會有好的效果，一般情況下應(yīng)通過試驗來選取.根據(jù)礦井瞬變電磁超前探測的特點，水平斷面上只有11個方向的探測數(shù)據(jù)，探測方向數(shù)較少，選擇了3個方向和4個方向上的數(shù)據(jù)分別進行合成孔徑試驗，最終確定選擇3個方向進行合成孔徑成像.

圖6為巷道迎頭G22+5m位置處礦井瞬變電磁超前探測合成孔徑成像后擬地震信號，可以看出，波場變換曲線的極值位置沒有改變，深度方向未發(fā)生偏移，即合成孔徑處理不影響波場曲線極值點與深度的對應(yīng)關(guān)系.

對比該斷面合成孔徑前(圖5)與合成孔徑后(圖6)擬地震信號波形圖，合成孔徑之后的波形圖在相關(guān)性較好的探測方向較合成孔徑之前的信號得到了明顯加強，如在相關(guān)性較好的-20°探測方向到+20°探測方向間的探測數(shù)據(jù)，合成孔徑前的兩處曲線幅值較為平緩，但合成孔徑后，波形幅值增大，信號得到了顯著增強，并且波形寬度未發(fā)生變化；而在相關(guān)性相對較差的-50°探測方向以及50°探測方向，合成孔徑前曲線幅值較為平緩，在合成孔徑后，波形幅值被削弱；其他探測方向合成孔徑后波場曲線幅值得到了不同程度的加強.整體而言，各個探測方向上有意義的弱異常得到了加強，信噪比得到了提高，合成孔徑效果比較顯著.

表2　G22+5 m位置處水平面內(nèi)相鄰探測方向擬地震數(shù)據(jù)間相關(guān)系數(shù)與相關(guān)性列表

由以上分析可以看出：合成孔徑處理不影響擬地震波場曲線極值點與深度的對應(yīng)關(guān)系.相關(guān)性較強的兩個探測方向上，其擬地震信號極值點在合成孔徑后得到了明顯增強；而相關(guān)性較弱的兩個探測方向上，其擬地震信號極值點在合成孔徑后被削弱.這說明合成孔徑成像在礦井瞬變電磁超前探測中突出弱異常、壓制假異常進而提高分辨率的有效性.

將圖6的擬地震波場曲線進行歸一化處理，并輔以實際斷面扇形坐標，標定電性界面位置，以利于解釋.圖7是G22+5 m位置瞬變電磁超前探測合成孔徑成像斷面圖與常規(guī)視電阻率斷面圖對比，可以看出：常規(guī)視電阻率斷面圖上有一明顯的低阻異常，位于巷道前方約20～60 m范圍，中心在34 m附近；合成孔徑成像斷面圖有兩個電性界面，第一電性界面位于前方35 m處，信號能量較強，應(yīng)該是高阻到低阻的突變界面；第二個電性界面位于前方73 m處，信號能量相對較弱，應(yīng)該是低阻到高阻的突變界面.兩處界面大體上為直線，并與探測方向近垂直.巷道掘進實際揭露了探測位置為前方33～70 m范圍砂巖含水，為弱富水，與合成孔徑成像解釋成果吻合較好，而常規(guī)視電阻率斷面圖所反映的異常與實際揭露存在一定的偏差.這進一步說明，僅依據(jù)斷面圖視電阻率的變化來確定電性界面位置存在較大誤差，而合成孔徑成像確定電性界面具有明顯優(yōu)勢.

圖6　G22+5 m位置處TEM超前探測合成孔徑成像后擬地震信號Fig.6　Pseudo-seismic wave-field signalafter SRIon horizontal section by MTEM at G22+5 m

圖7　合成孔徑成像斷面與常規(guī)視電阻率扇形斷面對比(a) 常規(guī)視電阻率斷面； (b) 合成孔徑成像斷面.Fig.7　Comparison ofgeneral resistivity section (a) and SAI section (b)

5結(jié)論

論文在將礦井瞬變電磁波場轉(zhuǎn)換為擬地震波場的基礎(chǔ)之上，實現(xiàn)了礦井瞬變電磁合成孔徑成像，并結(jié)合工程應(yīng)用實例對方法的有效性和實用性進行了探討和驗證，主要結(jié)論如下：

(1) 通過對比分析合成孔徑前后波形圖可知，合成孔徑之后的波形圖在相關(guān)性較好的探測方向上較合成之前的信號得到了加強，且波形寬度未發(fā)生變化；而在相關(guān)性相對較差的方向上合成孔徑后較合成孔徑前由噪聲引起的小幅值異常進一步減小，降低了隨機噪聲，提高了信噪比.

(2) 合成孔徑成像的效果與各探測方向之間的相關(guān)系數(shù)或相關(guān)性密切相關(guān)，較大的相關(guān)系數(shù)會獲得較好的合成孔徑成像效果.

(3) 對實際資料的處理和分析可知，低阻異常體在常規(guī)的視電阻率等值線圖中表現(xiàn)為一個放大的低阻區(qū)域，而在合成孔徑成像后的成果圖上則表現(xiàn)為更為精確的兩個電性界面.合成孔徑成像處理后電性界面的分辨率顯著提高.一般來說，在高阻到低阻的突變界面處，波場變換后信號較強，對應(yīng)界面合成孔徑成像后的能量較強；而在低阻到高阻的突變界面處，波場變換后信號較弱，對應(yīng)界面合成孔徑成像后的能量也相對較弱.對于后者，增強波場變換后信號有待進一步研究.

(4) 合成孔徑成像方法能夠突出瞬變電磁數(shù)據(jù)中所包含的電性界面信息，彌補了礦井瞬變電磁法在異常體深度確定方面的不足，提高信噪比，突出弱異常進而提高分辨率及勘探精度，發(fā)展了礦井瞬變電磁法對低阻異常電性界面的精細探測.

References

Ataei A. 2014. Improved Qrginv algorithm for computing moore-penrose inverse matrices.ISRNAppliedMathematics, 2014: Article ID 641706.

Cetin M, Karl W C. 2001. Feature-enhanced synthetic aperture radar image formation based on nonquadratic regularization.IEEETransactionsonImageProcessing, 10(4): 623-631.

Chan Y K, Koo V C. 2008. An introduction to synthetic aperture Radar(SAR).ProgressinElectromagneticsResearchB, 2: 27-60. Chen B C, Li J M, Zhou F T. 1999. Wave-field conversion method for transient electromagnetic field.OGP(in Chinese), 34(5): 539-545.

Cheng J L, Li F, Peng S P, et al. 2014. Research progress and development direction on advanced detection in mine roadway working face using geophysical methods.JournalofChinaCoalSociety(in Chinese), 39(8): 1742-1750.

Cheng J L, Li F, Peng S P, et al. 2015. Joint inversion of TEM and DC in roadway advanced detection based on particle swarm optimization.JournalofAppliedgeophysics,123: 30-35.

Cheng J L, Qiu H, Ye Y T, et al, 2013. Research on wave-field transformation and data processing of the mine transient electromagnetic method.JournalofChinaCoalSociety(in Chinese), 38(9): 1646-1650. Ganguli A, Rappaport C M, Abramo M, et al. 2012. Synthetic aperture imaging for flaw detection in a concrete medium.NDT&EInternational, 45(1): 79-90.

Guo W B, Xue G Q, Li X, et al. 2012. Correlation analysis and imaging technique of TEM data.ExplorationGeophysics, 43(3): 137-148.

Jensen J A, Nikolov S I, Gammelmark K L, et al. 2006. Synthetic aperture ultrasound imaging.Ultrasonics, 44: e5-e16. Kaufman A A, Eaton P A. 2001. The Theory of Inductive Prospecting. New York: Elsevier,365-373. Lee K H, Liu G, Morrison H F. 1989. A new approach to modeling the electromagnetic response of conductive media.Geophysics, 54(9): 1180-1192.

Li X, Qi Z P, Liu Y A, et al. 2010. A research of synthetic aperture rapid imaging of ATEM. ∥26rdCGS, Chinese Geophysics 2010 (in Chinese). 639. Li X, Xue G Q. 2013. Study on Pseudo-Seismic Migration Imaging of Transient Electromagnetic Method. Beijing: Science Press.

Li X, Xue G Q, Liu Y A, et al. 2012. A research on TEM imaging method based on synthetic-aperture technology.ChineseJ.Geophys. (in Chinese), 55(1): 333-340, doi: 10.6038/j.issn.0001-5733.2012.01.034.

Li Z H, Huang Q H. 2014. Application of the complex frequency shifted perfectly matched layer absorbing boundary conditions in transient electromagnetic method modeling.ChineseJ.Geophys. (in Chinese), 57(4): 1292-1299, doi: 10.6038/cjg20140426. Mason I, Osman N, Liu Q, et al. 2001. Broadband synthetic aperture borehole radar interferometry.JournalofAppliedGeophysics,47(3-4): 299-308.

Qi Z P. 2011. 3D continuation imaging technology with a synthetic aperature based on the transient electromagnetic method (in Chinese). Xi′an: Chang′an University.

Samadi S, ?etin M, Masnadi-Shirazi M A. 2009. Multiple feature-enhanced synthetic aperture radar imaging.∥Proc. SPIE Algorithms Synthetic Aperture Radar Imagery, XVI.Orlando, Florida, USA. To A C, Moore J R, Glaser S D. 2009. Wavelet denoising techniques with applications to experimental geophysical data.SignalProcessing, 89(1): 144-160. Wang Z L, Zhou Y P, Zhang W, et al. 2004. Synthetic aperture imaging techniques.JournalofHarbinInstituteofTechnology(in Chinese), 36(3): 384-387.

Yang H Y, Deng J Z, Zhang H, et al. 2010. Research on full-space apparent resistivity interpretation technique in mine transient electromagnetic method.ChineseJ.Geophys.(in Chinese), 53(3): 651-656, doi: 10.3969/j.issn.0001-5733.2010.03.020.

Yang W C. 1997. Theory and Methods of Geophysical Inversion (in Chinese). Beijing: Geological Publishing House. 45-67.

Zhang C C, Zhou Z O. 2004. Ground penetrating radar synthetic aperture imaging based on Stolt migration.ChineseJournalofRadioScience(in Chinese), 19(3): 316-320.

附中文參考文獻

陳本池, 李金銘, 周鳳桐. 1999. 瞬變電磁場的波場變換算法. 石油地球物理勘探, 34(5): 539-545.

程久龍, 李飛, 彭蘇萍等. 2014. 礦井巷道地球物理方法超前探測研究進展與展望. 煤炭學(xué)報, 39(8): 1742-1750.

程久龍, 邱浩, 葉云濤等. 2013. 礦井瞬變電磁法波場變換與數(shù)據(jù)處理方法研究. 煤炭學(xué)報, 38(9): 1646-1650.

李貅, 戚志鵬, 劉銀愛等. 2010. 航空瞬變電磁合成孔徑快速成像方法研究.∥中國地球物理2010-中國地球物理學(xué)會第二十六屆年會論文集, 639.

李貅, 薛國強. 2013. 瞬變電磁法擬地震偏移成像研究. 北京: 科學(xué)出版社.

李貅, 薛國強, 劉銀愛等. 2012. 瞬變電磁合成孔徑成像方法研究. 地球物理學(xué)報, 55(1): 333-340, doi: 10.6038/j.issn.0001-5733.2012.01.034.

李展輝, 黃清華. 2014. 復(fù)頻率參數(shù)完全匹配層吸收邊界在瞬變電磁法正演中的應(yīng)用. 地球物理學(xué)報, 57(4): 1292-1299, doi: 10.6038/cjg20140426.

戚志鵬. 2011. 瞬變電磁三維合成孔徑延拓成像方法研究. 西安: 長安大學(xué).

王治樂, 周彥平, 張偉等. 2004. 合成孔徑成像技術(shù)對比研究. 哈爾濱工業(yè)大學(xué)學(xué)報, 36(3): 384-387.

楊海燕, 鄧居智, 張華等. 2010.礦井瞬變電磁法全空間視電阻率解釋方法研究. 地球物理學(xué)報, 53(3): 651-656, doi: 10.3969/j.issn.0001-5733.2010.03.020. 楊文采. 1997. 地球物理反演的理論與方法. 北京: 地質(zhì)出版社,45-67. 張春城, 周正歐. 2004. 基于Stolt偏移的探地雷達合成孔徑成像研究. 電波科學(xué)學(xué)報, 19(3): 316-320.

(本文編輯張正峰)

基金項目國家自然科學(xué)基金項目(51574250, 51174210)和高等學(xué)校博士學(xué)科點專項科研基金項目(20120023110014)聯(lián)合資助.

作者簡介程久龍，男，1965年生，教授，博士生導(dǎo)師，主要從事應(yīng)用地球物理科研與教學(xué)工作.E-mail：jlcheng@126.com

doi:10.6038/cjg20160230 中圖分類號P631

收稿日期2015-02-09，2015-09-28收修定稿

Synthetic aperture imaging in advanced detection of roadway using the mine transient electromagnetic method

CHENG Jiu-Long1,2, CHEN Ding2, XUE Guo-Qiang3, QIU Hao4, ZHOU Xiang-Tong2

1StateKeyLaboratoryofCoalResourcesandSafeMining,ChinaUniversityofMiningandTechnology,Beijing100083,China2CollegeofGeoscienceandSurveyingEngineering,ChinaUniversityofMiningandTechnology,Beijing100083,China3KeyLaboratoryofMineralResources,InstituteofGeologyandGeophysics,ChineseAcademyofSciences,Beijing100029,China4MineSafetyTechnologyBranchofChinaCoalResearchInstitute,Beijing100013,China

AbstractDuring the data processing and interpretation of the mine transient electromagnetic method (MTEM), it is an important goal to improve the identification of the geo-electrical interface between low resistivity abnormal body and surrounding rock in advanced detection. To achieve this end, the transient electromagnetic field is transformed into the pseudo-seismic wave-field based on the function relationship between spreading electromagnetic field and seismic wave-field in whole-space. Meanwhile, the signal of pseudo-seismic after transform is processed with correlative stack to strengthen the amplitude and improve signal-noise-ratio (SNR) using the synthetic aperture imaging(SAI). The data of different detection directions on one surveying point in the MTEM are imaged as aperture data, which would highlight weak anomalies and improve SNR. In addition, the application effect of SAI was testified with the in-site advanced detection in roadway using the MTEM.The results show that the SAI of the MTEM can improve SNR of the electrical interface information and highlight the geometrical resolution to increase exploration accuracy. The technique is significant theoretically and practically in accurate advanced detection of aquifer structure by the MTEM.

KeywordsMine transient electromagnetic method； Wave-field transformation； Pseudo-seismicsignal； Synthetic aperture imaging； Geo-electrical interface

程久龍， 陳丁， 薛國強等. 2016. 礦井瞬變電磁法超前探測合成孔徑成像研究.地球物理學(xué)報，59(2):731-738,doi:10.6038/cjg20160230.

Cheng J L, Chen D, Xue G Q, et al. 2016. Synthetic aperture imaging in advanced detection of roadway using the mine transient electromagnetic method.ChineseJ.Geophys. (in Chinese),59(2):731-738,doi:10.6038/cjg20160230.