二維探地雷達正演

1、引言地質(zhì)雷達正演模擬是分析探測問題、研究電磁波在介質(zhì)中傳播規(guī)律的有效手段，在地質(zhì)雷達的實際應(yīng)用及技術(shù)研究中具有十分重要的意義。通過正演模擬各類地質(zhì)模型的數(shù)據(jù)，在很大程度上方便了對地質(zhì)雷達技術(shù)的研究，而不必花費大量時間去實地采集，大大節(jié)約了時間和成本。使用合理的正演方法能夠得到正確的、實際的雷達反射信號，能讓我們更加客觀的了解地下介質(zhì)的分布、物理特性等重要信息?？傊刭|(zhì)雷達正演模擬技術(shù)在地質(zhì)雷達領(lǐng)域中發(fā)揮著重要的作用。關(guān)于GPGMAX的介紹GPRMAX是一種探地雷達正演模擬軟件，其基礎(chǔ)是時域有限差分（FDTD）,由AntonisGiannopoulos教授研發(fā)。利用GPRMAX做探地雷達正演數(shù)

2、值模擬，前提是需要編寫關(guān)于要模擬的模型數(shù)據(jù)的輸入文件，輸入文件中的命令可分為以下四個步驟：1）描述性命令語句。描述GPR所要掃描的模型，即模型的大小和掃描的步長。2）邊界吸收相關(guān)命令語句。指明邊界吸收條件及公式。3）所要模擬的介質(zhì)及幾何體的相關(guān)命令語句。說明模型中所含的各種介質(zhì)的一些特性及埋于介質(zhì)中的一些簡單的不同尺寸的幾何體。4）GPR掃描及輸出的相關(guān)命令。在模型中放置觸發(fā)點，輸出點及掃描的道數(shù)。上述四個步驟完成后，在安裝目錄下單擊GprMax2D.exe文件，啟動后在光標(biāo)處輸入文件的路徑名（如.in或者.txt）,如圖1。若輸入文件的格式正確，軟件才會開始數(shù)值計算，反之則會自行終止。最后輸

3、出文件會自動生成，分別是.geo及.out類型。圖1圖1GprMax2D軟件啟動界面圖形顯示計算完.geo和.out類型文件之后，下一步就是利用matlab軟件顯示數(shù)據(jù)圖像。.geo輸出的是模型數(shù)據(jù)，.out輸出的是模型電磁波正演數(shù)據(jù)。利用matlab程序，可以數(shù)據(jù)圖形化。具體程序如下圖。其中每個程序包含一個調(diào)用函數(shù)，也就是子程序，該程序在下載軟件時別人就編輯好了，我們可以直接利用。fil&g4o=bre4.geo：meshdata5headerjmedia=gpr)nax2g(filegeo);figured);M,NN=M,NN=size(meshdata):imagesc(1:MN)*h

4、eader.dx2axis(equal);(1header,dy,meshdata)xlabelCx(m);ylabelCy(jr);圖2,geo文件matlab圖形顯示程序fileout=bre4.out;|HeaderFields=gprm.ax(fileout):N=l:Header.NSteps;%移動次數(shù)Position=Header.dx*Header.(M-L)*(Header.dx*Header.TxStepX);%天線每次斫在位矍Data(:,:)=Fields.ez(:,1,:%轉(zhuǎn)換額組格式figure(2):iiuagesc(Position,Fields.t*le9,D

5、ata):%畫圖colorbarxlabelCAntenaPosition(m)；ylabelCt(ns):圖3.out文件matlab圖形顯示程序GPRMAX的理論基礎(chǔ)GprMax2D基于FDTD算法,FDTD法是K.S.Yee在1966年提出來的。該方法直接求解依賴于時間變量的Maxwell旋度方程組：AXH卷+JXE=-掌-J*(2)dt，XB=p*(3)XD=p其中E,H,D和B分別是電場強度，磁場強度,電位移矢量和磁感應(yīng)強度。J和J*分別是電流密度和磁流密度。p和p*分別是電荷密度和磁荷密度。簡而言之即是將旋度方程化為一組電場和磁場各分量的偏微分方程，然后將電場和磁場各分量智交叉取樣

6、，利用二階精度的中心差近似將這一組偏微分算符轉(zhuǎn)換為差分形式,這樣達到在一定空間和一段時間上對邊界電磁場的數(shù)據(jù)抽樣，在時域?qū)﹄姶抛饔眠^程進行直接模擬。旋度方程的差分化運用于每一個FDTD網(wǎng)格，而數(shù)值模擬的結(jié)果將直接通過這些網(wǎng)格內(nèi)方程的差分化得到，并且是多次重復(fù)計算的結(jié)果。在每一次的重復(fù)過程中，電磁波都傳遞到FDTD網(wǎng)格，每一次所耗費的時間即是to因此如果給定重發(fā)的次數(shù)(我們把它成為掃描的道數(shù))，時間窗(掃描一道所需的時間),我們就能知道FDTD解答器模擬指定的一塊區(qū)域的耗時。FDTD法對整個計算空間劃分網(wǎng)格。為保證計算精確度，通常每波長至少用10個以上網(wǎng)格。時間步長的確定則利用數(shù)值穩(wěn)定性條件確定

7、。時域有限差分法的邊界吸收條件吸收邊界處理的好壞將直接影響包括計算精度和計算開銷在內(nèi)的時域有限差分法性能,因此吸收邊界條件始終是一個重點研究內(nèi)容。理想匹配層(PerfectlyMatchLayer:PML)是J.Be-renger于1994年首先提出,其基本思想是在計算區(qū)域邊界面附近引入虛擬各向異性有耗媒質(zhì),在一定的條件下,模擬空間與理想匹配層間、理想匹配層內(nèi)部層間完全匹配，模擬區(qū)域內(nèi)的外行電磁波可以無反射地進入有耗媒質(zhì)，并在有耗媒質(zhì)內(nèi)進行衰減，從而有效吸收模擬區(qū)域內(nèi)出射的外行波。直角坐標(biāo)中三維E場量在x方向和z方向的FDTD時域推進計算式，完整公式推導(dǎo)及其變量含義可參考文獻。模型舉例本文設(shè)計

8、了幾個模型參數(shù)，經(jīng)過GPRMAX正演計算以及matlab成圖，驗證了程序的正確性，通過整個過程的實踐，對探地雷達正演有了全新的認(rèn)識。另外通過改變模型參數(shù)，深刻理解了模型參數(shù)變化和正演曲線的相關(guān)關(guān)系，對探地雷達的學(xué)習(xí)大有裨益。這里設(shè)置了兩組模型參數(shù)，分別如下：Aiediuri：6.00.00.00.OOS1.00.0cccicmieilTediur:20.00.00.00.11.00.0wet.setdAdcmoln:2.60.2Odx.dy：0.00250.00251titre_wjrdow：ftbox：0.00.02.50.15wct_3cndAcylirder:0.150.10,0125p

9、ecflcylinder：0.3S0.0750.0125pecrtcylirvdcr：0.SS0.070.012Sgcdcylirder:0.7S0.050.0125pecftcylinder：0.260.10.C6T&free_SFacrrtcylir/Jer:1.1$6.0750.0$?54cylirrder:I.350.070.0375free_cpaceAcylirvdcr：1.550.050.0375frcjzpacertcylinder:1.750.10.CO5pec0cylinder：I.60.0750.0C6pecrtcylirvdcr：2.150.070.005pecUcyl

10、irxder:2.SS0.050.006pecrtline-source:1.090Oe6rickerlyLlneSourceAanalysis：115bre4,outbaix:0.0850.1525VyLineSourcc0.012e-flrx:0.11250.15250.02(X0iJrx.fitcpfi：0.020.0fteji_nalysis：bre4.geoAtltle：mNoisl4口rrcdiuB：6.00.00.00.0051.00.0cctietc時cdiuB：3.00.00.00.OG011.00.0dry_sandtfdnnain：2.50.650.00250.0025#

11、tim_昉n&M：12l9#box:0.00.02.50.6dry_sand#box:0.00.32.50.6ccncrete#box:1.750.21.95U.Sfree.spaceftcylindn*:0.650.10.0375free_space#linc_sourK：LO900e6rickcrKyLineSourccffanalyris:IISbrel.outbfiftx:0.08750.4525VlineSource0.012e-90.112S0.4S2S,#tsf_stps:0.020.0ffrx_stps:0.020.0#end_analysis:#georretry_file：

12、brel.geoflftitle：Model1Messages:y圖4模型1相關(guān)參數(shù)圖5模型2相關(guān)參數(shù)圖6模型1正演模擬模型12248102E60.511.5AntenaPosition(m)1000圖6模型1正演模擬模型12248102E60.511.5AntenaPosition(m)1000800600400200-28-409圖7二維模擬數(shù)據(jù)剖面圖通過對模擬結(jié)果進行分析，能得到這些結(jié)果，圓形不良地質(zhì)體在雷達圖像上呈雙曲線形態(tài)，雙曲線特征為圓形不良地質(zhì)體雷達波發(fā)射所致。在不同物質(zhì)界面處，由于界面兩側(cè)物質(zhì)彈性波性質(zhì)不同，可以看到強烈的雷達波反射，如圖中黃色直線條帶所示。x(m)圖8模型2正

13、演模擬模型x(m)0.511.52AntenaPositon(m)圖9二維模擬數(shù)據(jù)剖面圖00.511.52AntenaPositon(m)圖9二維模擬數(shù)據(jù)剖面圖02810S10005000-5001000通過對模擬結(jié)果進行分析，能得到這些結(jié)果，在不同物質(zhì)界面處，由于界面兩側(cè)物質(zhì)彈性波性質(zhì)不同，可以看到強烈的雷達波反射;如圖中黃色直線條帶所示。但是由于兩處異常體埋藏較深，所以電磁波的反射能量較小，在正演模擬剖面圖上幾乎看不到。結(jié)論通過本文上述模型的正演結(jié)果可知，結(jié)合GprMax2D與MATLAB進行GPR正演數(shù)值模擬是可行的，并且效果顯著，正演模擬圖像中能清晰的反映出典型不良地質(zhì)體的形態(tài)特征,具有良好的應(yīng)用前景。參考文獻1李亞飛.地質(zhì)雷達超前地質(zhì)預(yù)報正演模擬D.北京交通大學(xué),2011.2周奇才，周杰，范思遐，張根雷.高鐵隧道GPR技術(shù)研究與應(yīng)用幾物探與化探,2014,01:185-188.周奇才，李炳杰,鄭宇軒，何自強.基于GPRMax2D的探地雷達圖像正演模擬人工程地球物理學(xué)報,2008,04:396-399.尹光輝.基于GprMax的道路空洞探地雷達圖像