手機直線：
咨詢熱線：0571-56211150
探地雷達(GPR)是采用高頻電磁波在地下電性界面(目標體)上的反射與透射信息來探測地質目標體的一種地球物理方法。具有高效、快速、無損、抗干擾能力強的優點，廣泛應用于工程與環境地球物理探測的各個領域，成為淺部勘探的重要方法技術。GPR正演模擬對數據解釋具有重要的指導作用。通過對GPR模型的正演模擬，可以加深對GPR探測剖面的認識，提高解釋精度。時域有限差分法(FDTD)以其算法簡單、靈活，程序易實現，而得到了廣泛的應用，但是FDTD法對于物性參數分布復雜或場域幾何特征不規則的模型，適應性差。為此，引入有限元分析法(FEM)，該算法能適用于求解各種復雜問題、求解過程規范，計算程序通用性強等優點具有廣泛的應用性。它在GPR正演模擬中應用起步較晚，主要是借鑒聲波、彈性波的有限元求解理論。如：將聲波方程代替電磁波的傳播方程進行了二維介質的正演模擬；底青云等進行了帶衰減項GPR剖面的FEM仿真模擬及偏移處理；謝輝等利用二十節點等參單元對Pulse探地雷達進行了有限單元正演模擬計算，建立了電磁波在Pulse探地雷達天線中的全反射模型；Arias等將GPR有限單元模擬技術應用到考古成像中。以往利用FEM進行GPR正演都是基于矩形網格剖分，對復雜GPR模型擬合效果較差，難以實現復雜地電模型的正演模擬。為此，本研究采用三角形剖分方式對復雜GPR地電模型進行FEM正演模擬，并把該結果與基于矩形剖分的FEM進行了對比。
設計了起伏界面、“V”字形界面2種GPR地電模型，應用基于三角形網格剖分的FEM算法結合透射邊界條件對設定的模型進行了正演模擬，并把該計算結果與基于矩形網格剖分的FEM算法進行了模擬效果對比。
起伏界面模型，模型為一個10 m×5 m的矩形區域，分為上下兩層介質，中間是起伏分界面，上層介質的相對介電常數為 5.0，電導率為0.001 s/m，下層介質的相對介電常數為50.0，電導率為0.01 s/m。空間網格步長為0.1 m，網格總數為100×50。GPR波脈沖激勵源的中心頻率為100 MHz，時間步長為0.25 ns，時窗長度為100 ns。分別應用基于三角形和矩形網格剖分的FEM對該模型進行正演模擬，其模擬所得的GPR正演模擬剖面如圖所示。由圖可見：在剖面中40 ns左右時刻，存在一能量較強反射界面，通過時深轉換計算可知，該反射界面與圖模型中的上、下兩層分界面位置相吻合，并且圖中變化劇烈的角點、界面起伏波動大、欠光滑的斷棱處，產生了角點繞射波形。但通過對比圖剖面可知，由于矩形剖分擬合起伏界面主要是以階梯狀模型擬合，三角形剖分能以三角單元的某條邊來更近似的擬合起伏的界面，因此三角形剖分擬合起伏分界面的能力強于矩形剖分，導致圖中產生的類似于角點繞射的波形能量明顯比圖的要弱。
專業從事機械產品設計│有限元分析│強度分析│結構優化│技術服務與解決方案
杭州納泰科技咨詢有限公司
本文出自杭州納泰科技咨詢有限公司www.hnyunhua.com，轉載請注明出處和相關鏈接！