[發明專利]地震數據時差、相位差自動識別及校正方法

有效

• 申請號: 201910572834.2
• 申請日: 2019-06-27
• 公開（公告）號: CN110261911B
• 公開（公告）日: 2021-07-20
• 發明（設計）人: 王東凱;李繼光;韓宏偉;趙愛國;韓智;金昌昆;張麗
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勝利油田分公司物探研究院
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 濟南日新專利代理事務所(普通合伙) 37224
• 代理人: 崔曉艷

1.地震數據時差、相位差自動識別及校正方法，其特征在于，該地震數據時差、相位差自動識別及校正方法包括：

步驟1，輸入標準地震道數據與待校正地震道數據；

步驟2，確定步長，在指定相位區間對待校正地震道數據進行相位校正掃描；

步驟3，對應每次掃描，對標準地震道數據和校正后的地震道數據做互相關處理；

步驟4，更新互相關值最大時對應的時差、相位差校正量；

步驟5，重復步驟2-4，完成在指定區間內地震道數據的相位掃描；

步驟6，完成時差、相位差校正量的自動識別；

步驟7，對地震道數據進行定量時差、相位差校正；

步驟8，輸出校正后地震道數據；

在步驟1中，輸入待連片的兩塊地震資料，確定拼接處的標準地震道數據A與待校正地震道數據B；

在步驟2中，以0.1度為相位掃描步長，[-π,π]為相位掃描區間，對地震道數據B信號S_B(t)進行相位校正掃描；對于指定的相位校正量θ，校正后的地震道數據B信號S_B(θ，t)表示為：

S_B(θ，t)＝S_B(t)·cosθ-H(S_B(t))·sinθ (式1)

其中H(S_B(t))為地震道數據B進行希爾伯特變換后的地震道數據；

在步驟3中，對標準地震道數據A信號S_A(t)和校正后的S_B(θ，t)做互相關處理，得到互相關函數R_AB(τ)的最大值R_max，具體過程表示為：

其中τ為相位校正后的地震道數據B在時間刻度上的平移量，對于長度相等且為N的兩道離散地震信號的互相關，τ的取值范圍為[-(N-1),N-1]；

在步驟4中，更新掃描過程中的互相關最大值R_max，并記錄對應的相位校正量θ和時間平移量τ；

在步驟5中，重復步驟2-4，計算下一個相位掃描間隔時地震道數據B的相位校正結果以及該結果與標準地震道數據A的互相關函數，記錄并更新互相關最大值R_max、相位校正量θ和時間平移量τ，直至全部計算完成指定區間內的相位掃描；

在步驟6中，完成時差、相位差校正量的自動識別，掃描過程結束后，更新停止時互相關最大值R_max對應的時間平移量τ和相位校正量θ就是地震道數據B到標準地震道數據A的時差和相位差；

在步驟7中，對地震道數據B所在的待拼接地震數據中的所有地震道進行定量的相位差和時差校正，校正量分別為θ和τ；

在步驟8中，將標準地震道數據A所在的待拼接地震數據與校正后的地震道數據B所在的待拼接地震數據進行連片處理，輸出一致性校正處理后的連片地震數據。

2.地震數據時差、相位差自動識別及校正方法，其特征在于，該地震數據時差、相位差自動識別及校正方法包括：

步驟1，輸入標準地震道數據與待校正地震道數據；

步驟2，確定步長，在指定相位區間對待校正地震道數據進行相位校正掃描；

步驟3，對應每次掃描，對標準地震道數據和校正后的地震道數據做互相關處理；

步驟4，更新互相關值最大時對應的時差、相位差校正量；

步驟5，重復步驟2-4，完成在指定區間內地震道數據的相位掃描；

步驟6，完成時差、相位差校正量的自動識別；

步驟7，對地震道數據進行定量時差、相位差校正；

步驟8，輸出校正后地震道數據；

在步驟1中，輸入多分量地震資料，包括加速度分量和速度分量，確定加速度分量中的標準地震道數據C與速度分量中的待匹配地震道數據D；

在步驟2中，以0.1度為相位掃描步長，為相位掃描區間，對待匹配地震道數據D信號S_D(t)進行相位校正掃描；對于指定的相位校正量θ，校正后的地震道數據D信號S_D(θ，t)表示為：

S_D(θ，t)＝S_D(t)·cosθ-H(S_D(t))·sinθ (式3)

其中H(S_D(t))為地震道數據D進行希爾伯特變換后的地震道數據；

在步驟3中，對標準地震道數據C信號S_C(t)和校正后的S_D(θ，t)做互相關處理，得到互相關函數R_CD(τ)的最大值R_max，具體過程表示為：

其中τ為相位校正后的地震道數據D在時間刻度上的平移量，對于長度相等且為N的兩道離散地震信號的互相關，τ的取值范圍為[-(N-1),N-1]；

在步驟4中，更新掃描過程中的互相關最大值，并記錄對應的相位校正量和時間平移量；

在步驟5中，重復步驟2-4，計算下一個相位掃描間隔時地震道數據D的相位校正結果以及該結果與標準地震道數據C的互相關函數，記錄并更新互相關最大值、相位校正量和時間平移量，直至全部計算完成指定區間內的相位掃描；

在步驟6中，完成時差、相位差校正量的自動識別，掃描過程結束后，更新停止時互相關最大值對應的時間平移量和相位校正量就是地震道數據D到標準地震道數據C匹配的時差和相位差；

在步驟7中，對地震道數據D所在的速度分量數據中的所有地震道進行步驟6中求取得到的定量相位差和時差校正；

在步驟8中，輸出校正處理后的速度分量數據，為后續與加速度分量數據的雙檢合成處理做準備。