[發明專利]寬頻帶頻散曲線提取方法、裝置

有效

• 申請號: 202010978546.X
• 申請日: 2020-09-16
• 公開（公告）號: CN112083487B
• 公開（公告）日: 2021-12-14
• 發明（設計）人: 李俊倫;鄧寶;沈長明;劉江山;高級
• 申請（專利權）人: 中國科學技術大學
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 中科專利商標代理有限責任公司 11021
• 代理人: 周天宇

1.一種寬頻帶頻散曲線提取方法，應用于密集線性臺陣，其特征在于，包括：

利用所述密集線性臺陣內的各臺站采集地震面波數據；

分別以所述密集臺陣內的各個臺站為中心臺站，將所述密集臺陣按一定半徑劃分為多個陣內臺站和多個陣外臺站；

基于各所述中心臺站對應的多個陣內臺站采集的地震面波數據，計算各所述中心臺站對應的陣內相移頻散能量，以及，基于全部臺站采集的地震面波數據，計算各所述中心臺站對應的陣外相移頻散能量，包括：

基于各所述中心臺站對應的所述多個陣內臺站采集的所述地震面波數據，利用互相關函數，計算與各所述中心臺站對應的第一正支頻散能量、第一負支頻散能量和第一對稱疊加頻散能量；

從所述第一正支頻散能量、第一負支頻散能量和第一對稱疊加頻散能量中選出一種頻散能量，作為所述中心臺站的陣內相移頻散能量；

基于全部臺站采集的所述地震面波數據，利用互相關函數，計算與各所述中心臺站對應的第二正支頻散能量、第二負支頻散能量和第二對稱疊加頻散能量；

從所述第二正支頻散能量、第二負支頻散能量和第二對稱疊加頻散能量中選出一種頻散能量，作為所述中心臺站的陣外相移頻散能量；

根據各所述中心臺站對應的所述陣內相移頻散能量和所述陣外相移頻散能量，提取出對應的陣內頻散曲線和陣外頻散曲線；

給各所述中心臺站對應的所述陣內頻散曲線和所述陣外頻散曲線在各頻率點分別分配權重，并進行加權疊加計算，使所述陣內頻散曲線和所述陣外頻散曲線融合為一條寬頻帶頻散曲線。

2.根據權利要求1所述的方法，其特征在于，所述利用所述密集線性臺陣內的各臺站采集地震面波數據之后，方法還包括：

對所述地震面波數據分別進行單道預處理。

3.根據權利要求1所述的方法，其特征在于，所述分別以所述密集臺陣內的各個臺站為中心臺站，將所述密集臺陣按一定半徑劃分為多個陣內臺站和多個陣外臺站包括：

若所述密集線性臺陣為臺間距為d的N個臺站組成，當以第n個臺站為中心臺站時，以r＝md為半徑作圓，所述圓內部的多個臺站為所述多個陣內臺站，所述圓外部的多個臺站成為所述多個陣外臺站，n＝1，…，N。

4.根據權利要求1所述的方法，其特征在于，所述基于各所述中心臺站對應的所述多個陣內臺站采集的所述地震面波數據，利用互相關函數，計算與各所述中心臺站對應的第一正支頻散能量、第一負支頻散能量和第一對稱疊加頻散能量包括：

令表示所述第一正支頻散能量，表示所述第一負支頻散能量，表示所述第一對稱疊加頻散能量，f表示地震面波的頻率，C_T表示所述地震面波的估計相速度，x_kl表示第k個陣內臺站和第l個陣內臺站之間的距離，分別為所述k，l兩臺站記錄的地震面波數據的互相關函數正支和負支的傅里葉變換，所述多個陣內臺站為所述密集線性臺陣內的第n-m至n+m個臺站，則：

5.根據權利要求1所述的方法，其特征在于，所述基于全部臺站采集的所述地震面波數據，利用互相關函數，計算與各所述中心臺站對應的第二正支頻散能量、第二負支頻散能量和第二對稱疊加頻散能量包括：

令表示所述第二正支頻散能量，表示所述第二負支頻散能量，表示所述第二對稱疊加頻散能量，f表示所述地震面波的頻率，C_T表示所述地震面波的估計相速度，x_kl表示第k個所述陣外臺站和第l個陣內臺站之間的距離，分別為k，l兩臺站記錄的地震面波數據的互相關函數正支和負支的傅里葉變換，所述密集線性臺陣包括N個臺站，所述多個陣內臺站為所述密集線性臺陣內的第n-m至n+m個臺站，所述多個陣外臺站為所述密集線性臺陣內的除所述多個陣內臺站以外的臺站，則：