[發明專利]一種探地雷達數據三維時頻分析方法

在審

• 申請號: 201910312340.0
• 申請日: 2019-04-18
• 公開（公告）號: CN110109186A
• 公開（公告）日: 2019-08-09
• 發明（設計）人: 許軍才;朱先祥;任青文;沈振中
• 申請（專利權）人: 河海大學
• 主分類號: G01V3/12
• 分類號: G01V3/12;G01V3/38;G01S7/36
• 代理公司: 南京經緯專利商標代理有限公司 32200
• 代理人: 許方

1.一種實現探地雷達數據三維時頻譜分析方法，其特征在于，該方法的步驟如下：

(1)選定需要處理測線上的探地雷達回波信號，根據回波信號強度調整增益；

(2)根據干擾來源判斷干擾波形和其頻率范圍，采用濾波方法去除背景干擾；

(3)對勘測剖面上每道回波信號分別進行計算得到二維瞬時頻率譜；

(4)將剖面每道回波信號的二維瞬時頻率譜進行合成，構建出該測線的三維時頻譜數據體H(x,t,ω)；

(5)根據三維探地雷達時頻譜數據體，在位置、時間和頻率三個不同方向進行切片，得到分頻剖面、時頻剖面及層位剖面以此分析判斷變化特征異常位置和大小。

2.根據權利要求1所述的一種實現探地雷達數據三維時頻譜分析方法，其特征在于，步驟(3)中，對勘測剖面上每道回波信號分別進行計算得到二維瞬時頻率譜，方法如下：首先通過變分模態分解將待處理的探地雷達信號分解為若干特征模態函數，再將分離的特征模態函數經過希爾伯特變換得出各個特征模態函數的瞬時頻率，最后根據特征模態函數的瞬時頻率譜疊加得到測點處回波信號的二維瞬時頻率譜。

3.根據權利要求1或2所述的一種實現探地雷達數據三維時頻譜分析方法，其特征在于，步驟(5)中，分析判斷變化特征異常位置和大小方法如下：從測線各頻帶能量變化特征、固定測點處頻率隨時間的變化特征及沿反射界面頻率變化特征，根據各剖面中頻率能量譜極值區域范圍及分布判斷異常位置和大小。