跨孔透地雷达断层扫描影像之解析度与精确度之研究.pdf

跨孔透地雷達斷層掃描影像之解析度與精確度之研究 張竝瑜 嘉南藥理科技大學環境資源管理系 摘要 本研究透過電磁波傳輸模擬跨孔透地雷達訊號，並運用阻尼最小平方法與商 業反演算程式分析斷層掃描資料，藉以比較不同方法之反演算影像解析度與準確 程度。驗證模型則採用一簡單砂泥互層模型與美國 STVZ 之穩態(Steady-state)與暫 態(Transient-state)之未飽和層地下水流動模型。結果顯示(1)在簡單砂泥互層模型 中，以阻尼最小平方法反演算之影像與運用程式 GeotomCG 準確度相當，均能顯 示砂質地層中所夾之泥層存在。(2)在穩態與暫態模型模擬試驗中，比較兩種不同 演算法，訊號衰減率反演算影像解析度均較波速反演算影像要差，顯示實際情況 下，雷達波傳遞可能有額外之能量損失，因此須於反演算的過程中加以考量。(3) 在暫態的試驗中，入滲試驗開始後第七日、以及第二十一日之反演算影像差值顯 示，由於阻尼最小平方法採用了阻尼以保證反演算之順利進行，ㄧ定程度減少了 反演算影像最高值與最低值間之誤差程度，解析度上要較 GeotomCG 來的差。但 比較所得之衰減率變化，阻尼最小平方法所得之值卻較 GeotomCG 接近原模型之 設定值，顯示以阻尼最小平方法所得之影像，要較 GeotomCG 能更接近模型原始 設定值，但卻犧牲其解析度，為在解析度與準確度間得到平衡，必須依在測試阻 尼值，以得到最佳之解析度與準確度組合。 前言 本研究透過電磁波傳輸模擬與掃描資料反演算分析，釐清造成透地雷達斷層 掃描之解析度與精確性之主要因子，以協助未來研究者利用最有效之透地雷達斷 層掃描，達成精確監測地層目標物(如地下污染物分布，未飽和層含水量分布， 或裂隙排列走向等)之目的。由於根據先前研究之發現，一般商用之為透地雷達 波速反演算開發之軟體，並不一定適用於訊號衰減率反演算之應用。因此本研究 採用自行開發撰寫之反演算程式，並採用簡單砂泥互層、穩態(Steady-state)與暫 態(Transient-state)之未飽和層地下水流動模型進行反演算，與目前之透地雷達掃 描反演算程式所反演之影像進行分析比較,以了解如何改進反演算之解析度與精 確度。本研究內容包含： (1) 透地雷達影像反演算程式之撰寫。 (2) 透地雷達斷層掃描反演算程式之檢討比較及反演算方法對掃描影像之解析度 探討。 (3)地層與地下水層物理狀況對透地雷達斷層掃描之誤差評估。 透地雷達斷層掃描(GPR Tomography)反演算(Inversion)原理 透地雷達斷層掃描方法(GPR Tomography)利用井間(Borehole-borehole)或地面- 井間(Surface-borehole)等不同之雷達天線排列，收集不同排列方式與不同天線排列 位置所得到之高密度雷達波訊號，經由反演算後推測介質中之物理性質分布狀 況。 對透地雷達斷層掃描資料反演算方法，一般均將發射與接收天線井間之區域 劃分為若干面積一致之小方格，假設雷達波由發射天線至接收天線行進的途徑可 視為垂直波前之直線進行，每一組天線排列位置所收集到之雷達波訊號資料如傳 遞時間或振幅衰減率可以以下式表示：(Sharma, 1997; Jackson Tweeton, 1996): M t (or α ) = ∑p d (2.1) i i j ij j =1 式 2.1 中，ti (αi ) 為某一組天線組合之雷達波行進路線 ith 所得到之雷達波傳 遞總時間與總訊號衰減率， p j 代表在某ㄧ方格 j th 內之平均波速或訊號衰減率， dij 為此雷達波行進路線 ith 在方格 j th 內之行進距離，而 M 則代表