謝建林 許家林 李曉林

(1.中國礦業(yè)大學(xué)礦業(yè)工程學(xué)院,江蘇省徐州市,221008;2.中國礦業(yè)大學(xué)煤炭資源與安全開采國家重點(diǎn)實(shí)驗(yàn)室,江蘇省徐州市,221008;3.中國礦業(yè)大學(xué)徐海學(xué)院,江蘇省徐州市,221008)

頂板離層檢測(cè)的地質(zhì)雷達(dá)物理模擬衰減分析＊

謝建林1,2,3許家林1,2李曉林1,2

(1.中國礦業(yè)大學(xué)礦業(yè)工程學(xué)院,江蘇省徐州市,221008;2.中國礦業(yè)大學(xué)煤炭資源與安全開采國家重點(diǎn)實(shí)驗(yàn)室,江蘇省徐州市,221008;3.中國礦業(yè)大學(xué)徐海學(xué)院,江蘇省徐州市,221008)

利用地質(zhì)雷達(dá)(GPR)對(duì)頂板離層物理模型進(jìn)行檢測(cè),根據(jù)菲涅耳公式在邊值關(guān)系條件下求出入射波、反射波和折射波的振幅關(guān)系,理論推導(dǎo)得出離層區(qū)的衰減大于非離層區(qū)的結(jié)論。通過物理模擬實(shí)驗(yàn)并對(duì)輸出的地質(zhì)雷達(dá)時(shí)間剖面圖像進(jìn)行特征分析,得出在離層界面的反射信號(hào)強(qiáng)烈造成了更大的衰減的結(jié)論。同時(shí)對(duì)接收的電磁波信號(hào)進(jìn)行低通濾波,通過對(duì)信號(hào)的振幅分析得出非離層區(qū)反射信號(hào)能量大于離層區(qū)反射信號(hào)能量的結(jié)論,驗(yàn)證了理論推導(dǎo)的正確性,同時(shí)也證明地質(zhì)雷達(dá)對(duì)頂板離層進(jìn)行檢測(cè)是可行的。

頂板離層 離層監(jiān)測(cè) 地質(zhì)雷達(dá) 物理模擬 衰減分析

目前人們認(rèn)為要做到準(zhǔn)確預(yù)防頂板冒頂事故,一是對(duì)冒頂機(jī)理進(jìn)行深入的研究;二是進(jìn)行合理的頂板控制設(shè)計(jì);三是進(jìn)行支護(hù)質(zhì)量與頂板動(dòng)態(tài)監(jiān)測(cè)。通過對(duì)支護(hù)質(zhì)量與頂板動(dòng)態(tài)監(jiān)測(cè),可以避免絕大部分冒頂事故的發(fā)生。但我國學(xué)者在研究我國近年發(fā)生較多的復(fù)合頂板推垮型事故時(shí),總結(jié)出這類冒頂事故發(fā)生前沒有明顯征兆,壓力不大,活柱下縮與頂?shù)装逑鲁烈膊幻黠@,可見這類事故通過支護(hù)質(zhì)量與頂板動(dòng)態(tài)監(jiān)測(cè)是很難預(yù)見的。造成這類頂板事故的主要因素是復(fù)合頂板形成的離層,因此了解頂板離層實(shí)際發(fā)展的狀況,從而判斷由離層導(dǎo)致的頂板事故發(fā)生的可能性,并為采取有效措施控制或避免事故發(fā)生提供有效地判斷依據(jù),確保頂板安全性很有必要。

本文對(duì)頂板離層進(jìn)行地質(zhì)雷達(dá)物理模擬,通過對(duì)輸出的地質(zhì)雷達(dá)時(shí)間剖面圖像以及信號(hào)振幅進(jìn)行分析,為地質(zhì)雷達(dá)對(duì)頂板離層進(jìn)行檢測(cè)提供理論依據(jù),并驗(yàn)證其方法的可行性。

1 地質(zhì)雷達(dá)模擬參數(shù)概述

地質(zhì)雷達(dá)方法是一種用于確定地下介質(zhì)分布的廣譜電磁技術(shù)。地質(zhì)雷達(dá)利用一個(gè)天線發(fā)射高頻寬頻帶電磁波,另一天線接收來自地下介質(zhì)界面的反射波。電磁波在介質(zhì)中傳播時(shí),其路徑、電磁場(chǎng)強(qiáng)度與波形將隨通過介質(zhì)的電性質(zhì)及幾何形態(tài)而變化。因此,根據(jù)接收到波的旅行時(shí)間、幅度與波形資料可推斷介質(zhì)的結(jié)構(gòu)。

地質(zhì)雷達(dá)利用寬帶的電磁波以脈沖形式來檢測(cè)地表之下或確定不可視的物體或結(jié)構(gòu)。經(jīng)過多年的發(fā)展,地質(zhì)雷達(dá)已成為工程和環(huán)境地球物理領(lǐng)域中的一種重要的方法,解決諸如公路路面層和基底結(jié)構(gòu)、高層建筑基底形態(tài)、近地表土壤層結(jié)構(gòu)、地下水位面和巖石分層、地下水污染評(píng)價(jià)等問題。

對(duì)于給定的介質(zhì),在一定的物理?xiàng)l件下(溫度、密度等),介質(zhì)的介電常數(shù)是定值,介電常數(shù)反映了處于電場(chǎng)中的介質(zhì)存儲(chǔ)電荷的能力。高頻電磁波在介質(zhì)中的傳播速度主要取決于介質(zhì)的介電常數(shù),其速度v(單位為m/ns)為:

式中:c——真空中的電磁波傳播速度,取0.3 m/ns;

εr——介質(zhì)的相對(duì)介電常數(shù),F/m。

如果要求的空間分辨率為x(單位為m),檢測(cè)介質(zhì)的相對(duì)介電常數(shù)為εr,則天線中心頻率f(單位為MHz)可由下式初步選定:

時(shí)窗大小W(單位為ns)主要取決于最大檢測(cè)深度H(單位為m)與地層電磁波的傳播速度v,應(yīng)滿足:

2 地質(zhì)雷達(dá)物理模擬理論分析

現(xiàn)假設(shè)模型一見圖1,介質(zhì)從上到下依次分別為空氣、潮濕砂巖與干沙質(zhì)土壤,假設(shè)其折射率依次分別為n1、n2與n3,電磁波的入射方向?yàn)閺目諝膺M(jìn)入砂巖。

圖1 模型一信號(hào)反射和折射示意圖

介質(zhì)的折射率和介質(zhì)的相對(duì)介電常數(shù)有關(guān),具體關(guān)系為:

空氣、潮濕砂巖與干沙質(zhì)土壤三種介質(zhì)的相對(duì)介電常數(shù)大小依次為1、8與5,由此可以計(jì)算得出n2最大,n3次之,n1最小。在圖1中,E1為發(fā)射源發(fā)出的入射電磁波的電場(chǎng)強(qiáng)度,E2為在空氣和砂巖界面反射波的電場(chǎng)強(qiáng)度,E3為在砂巖中折射波的電場(chǎng)強(qiáng)度,E4為在砂巖和土壤界面反射波的電場(chǎng)強(qiáng)度,E5為經(jīng)砂巖折射入空氣折射波的電場(chǎng)強(qiáng)度。我們接收的反射波主要是E2和E5,在這里為了計(jì)算方便將其他信號(hào)進(jìn)行忽略。

那么根據(jù)菲聶耳垂直入射公式,就有:

R1為圖1中接收到的電場(chǎng)強(qiáng)度的值。

現(xiàn)假設(shè)模型二見圖2,介質(zhì)從上到下依次分別為空氣、潮濕砂巖、空氣與干沙質(zhì)土壤,假設(shè)其折射率依次分別為n1、n2與n3,電磁波的入射方向?yàn)閺目諝膺M(jìn)入砂巖。與圖3相比,就是砂巖介質(zhì)與土壤介質(zhì)中存在空氣層,把這一層看成是離層區(qū)。

圖2 模型二信號(hào)反射和折射示意圖

從圖2可以看出,E10為發(fā)射源發(fā)出的入射電磁波的電場(chǎng)強(qiáng)度,E20為在空氣和砂巖界面反射波的電場(chǎng)強(qiáng)度,E30為在砂巖中折射波的電場(chǎng)強(qiáng)度,E40為在砂巖和空氣界面反射波的電場(chǎng)強(qiáng)度,E50為經(jīng)砂巖折射入空氣折射波的電場(chǎng)強(qiáng)度。我們接收的反射波主要是E20和E50,在這里為了計(jì)算方便將其他信號(hào)進(jìn)行忽略。

R2為圖2中接收到的電場(chǎng)強(qiáng)度的值。R1與R2相比,E2和E20是相同的,而后半部分E5和E50有差異,其中:

所以得出E5

3 地質(zhì)雷達(dá)物理模擬實(shí)驗(yàn)

3.1 地質(zhì)雷達(dá)物理模型

本次物理模擬實(shí)驗(yàn)的模型剖面圖見圖3,模型整體高度為45 cm,離層區(qū)厚度為15 cm。圖3中空氣介質(zhì)區(qū)域假設(shè)為離層區(qū)。

實(shí)驗(yàn)采用500 MHz地質(zhì)雷達(dá)屏蔽天線,檢測(cè)方向?yàn)閺哪Ｐ妥蠖碎_始向離層區(qū)方向移動(dòng),檢測(cè)步長(zhǎng)為2 cm,時(shí)窗大小為95 ns,檢測(cè)總道數(shù)為165,從第44道開始進(jìn)入離層區(qū)。

圖3 模型剖面圖

3.2 地質(zhì)雷達(dá)物理模擬實(shí)驗(yàn)結(jié)果

對(duì)模型進(jìn)行檢測(cè)得到的同相軸圖進(jìn)行處理后見圖4。

圖4 同相軸圖

從圖4可以觀測(cè)得到在橫向44道后、縱向8 ns附近區(qū)域同相軸出現(xiàn)強(qiáng)烈的反射,根據(jù)地質(zhì)雷達(dá)檢測(cè)目的層深度的計(jì)算式可計(jì)算得出此區(qū)域?qū)?yīng)的正是模型中的離層區(qū)。

從全部信道中選擇第33道(見圖5)以及第100道(見圖6)離層區(qū)信號(hào)頻譜圖進(jìn)行對(duì)比分析。

對(duì)圖5和圖6兩道信號(hào)進(jìn)行低通濾波后,得到處理后的信號(hào)見圖7和圖8。

從圖7和圖8可以觀測(cè)得到離層區(qū)信號(hào)幅值略小于平均道信號(hào)振幅,而非離層區(qū)信號(hào)幅值大于平均道信號(hào)振幅,這驗(yàn)證了之前離層區(qū)的衰減大于非離層區(qū)的結(jié)論。

4 結(jié)論

本文通過物理模擬實(shí)驗(yàn)并對(duì)輸出的地質(zhì)雷達(dá)時(shí)間剖面圖像進(jìn)行特征分析,得出在離層界面的反射信號(hào)強(qiáng)烈造成了更大的衰減,從圖像上能夠直觀區(qū)分離層區(qū)與非離層區(qū)。同時(shí)對(duì)接收的電磁波信號(hào)進(jìn)行低通濾波,通過對(duì)接收信號(hào)的振幅分析得出非離層區(qū)反射信號(hào)能量大于離層區(qū)反射信號(hào)能量。綜上所述,驗(yàn)證了離層區(qū)的存在增加了信號(hào)衰減的理論推導(dǎo)正確性,同時(shí)也證明采用地質(zhì)雷達(dá)對(duì)煤礦頂板離層進(jìn)行檢測(cè)是可行的。

[1] 岑傳鴻.頂板災(zāi)害防治(第二版)[M].徐州:中國礦業(yè)大學(xué)出版社,1994

[2] 岑傳鴻,高存寶,杜計(jì)平.論預(yù)防冒頂事故的采場(chǎng)頂板控制[J].中國礦業(yè)大學(xué)學(xué)報(bào),1996(6)

[3] Qian Minggao.A study of the behaviour of overlying strata in longwall mining and its application to strata control[M].StrataMechanics,ElsevierScientific Publishing Company,1982

[4] 李大心編著.探地雷達(dá)方法與應(yīng)用[M].北京:地質(zhì)出版社,1994

[5] 劉四新,曾昭發(fā).頻散介質(zhì)中地質(zhì)雷達(dá)波傳播的數(shù)值模擬[J].地球物理學(xué)報(bào),2007(1)

[6] 鄧世坤.探地雷達(dá)技術(shù)野外工作參數(shù)選擇的基本原理[J].工程地球物理學(xué)報(bào),2005(5)

[7] 郭碩鴻編著.電動(dòng)力學(xué)(第二版)[M].北京:高等教育出版社,1997

Analysis of physically simulated attenuation of GPR in detecting mine roof separation

Xie Jianlin1,2,3,Xu Jialin1,2,Li Xiaolin1,2
(1.School of Mines,China University of Mining and Technology,Xuzhou,Jiangshu 221008,China;2.State Key Laboratory of Coal Resources and Safe Mining,China University of Mining
and Technology,Xuzhou,Jiangshu 221008,China;3.Xuhai College,China University of Mining and Technology,Xuzhou,Jiangshu 221008,China)

Theoretical derivation concludes that attenuation in the separation zone is higher than in areas without roof separation.The attenuation is detected by testing the physical model of roof separation using the Ground Penetrating Radar(GPR).The derivation is based on the relationship between the amplitudes of the incident wave,reflected wave and refractive wave calculated in the Fresnel formula under the boundary condition.Physical simulation and analysis of the characteristics of the time section as output from the GPR indicate that strong reflected signals in the interface of separation result in higher attenuation.Lowpass filtering of

electromagnetic signals is conducted,and analysis of the amplitudes of signals indicates that the energy of the reflected signals from separation-free areas is greater than those from areas with roof separation,which verifies the correctness of the theoretical derivation and the feasibility of using the GPR to detect roof separation.

roof separation,separation monitoring,GPR,physical simulation,attenuation analysis

TD325

A

國家自然科學(xué)基金項(xiàng)目(50974116);國家重點(diǎn)基礎(chǔ)研究計(jì)劃(973項(xiàng)目)(2007CB209403)

謝建林(1982-),男,江西泰和人,博士研究生,講師,從事巖層移動(dòng)方面的研究。

(責(zé)任編輯 張艷華)