探地雷達因其無損探測的特性被廣泛應用。然而,無遮蔽雷達天線激發(fā)的電磁波可以在空氣中傳播,當遇到地面異常體(如樹木和石塊)時將產生繞射回波(圖1),從而對地下反射波和繞射波造成干擾,必須設法壓制。現(xiàn)有的地表繞射波分離方法效果不佳,對有效信號損傷較為明顯。
圖1 探地雷達觀測系統(tǒng)與地表繞射波示意圖
中科院地質與地球物理所地球與行星物理院重點實驗室的博士生李超與導師張金海研究員等合作,開發(fā)了一種基于迭代Stolt偏移/反偏移的地表繞射波分離方法。他們首先利用真空中的光速開展Stolt聚焦偏移以識別繞射波,并對識別出的偏移域繞射波頂點進行反偏移,從而生成未偏移域的窗函數(shù),并對輸入的探地雷達數(shù)據(jù)進行分割,形成以該繞射波為主體的數(shù)據(jù)體(圖2c),剩余的數(shù)據(jù)體(圖2b)則以地下反射波和繞射波為主。利用真空中的光速對窗內的數(shù)據(jù)進行Stolt偏移,從而使地表繞射波聚焦在偏移域中的局部區(qū)域而免受窗外信號對其形成干擾(圖2d)。相比之下,原本被地表繞射波所掩蓋的有效信號(即窗口內的有效信號)經過偏移后會被拉伸,從而與聚焦的地表繞射波進一步分離,只有極小部分仍被覆蓋(圖2d)。此時,在偏移域利用新的窗函數(shù)再一次分離地表繞射波和有效信號,并且分別對二者進行反偏移,就能得到第一次迭代恢復的有效信號(圖2h)和地表繞射波(圖2g)。最后,分離出的地表繞射波中仍殘存一些微弱的地下有效信號(圖2g),將其作為下一輪迭代的輸入數(shù)據(jù),但采用亞光速和超光速進行完全一致的二次乃至三次分離,就可以進一步分選出更多殘存的地下有效信號(圖3)。
圖2 基于光速的聚焦偏移及地表繞射波分離過程及結果
圖3 基于多速度迭代偏移的地表繞射波分離過程及結果
探地雷達野外數(shù)據(jù)測試結果表明(圖4):新方法能夠有效去除地表繞射波,且?guī)缀醪粨p傷地下的反射波和繞射波。該方法對于裸露天線探地雷達的數(shù)據(jù)處理和地下結構成像具有重要意義,對于地外天體雷達探測也具有實際意義。
圖4 實際數(shù)據(jù)地表繞射波分離結果(左側為孤立地表繞射波,右側為重疊地表繞射波)
研究成果發(fā)表于勘探地震學領域權威期刊Geophysics (李超, 林羿, 呂文敏, 張金海*. Eliminating above-surface diffractions from ground-penetrating radar data using iterative Stolt migration[J]. Geophysics, 2021, 86(1): H1-H11. DOI: 10.1190/GEO2019-0796.1)



