[發(fā)明專利]地表一致性和非一致性聯(lián)合靜校正方法

在審

• 申請?zhí)枺?/td>: 202011102822.2
• 申請日: 2020-10-15
• 公開（公告）號: CN112269210A
• 公開（公告）日: 2021-01-26
• 發(fā)明（設計）人: 楊華臣;張建中;任凱
• 申請（專利權(quán)）人: 中國海洋大學
• 主分類號: G01V1/36
• 分類號: G01V1/36;G01V1/32;G01V1/28
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.地表一致性和非一致性聯(lián)合靜校正方法，其特征在于，包括如下步驟：

步驟一：對陸上地震資料進行預處理，預處理包括振幅均衡、去除面波、壓制隨機噪聲；

步驟二：利用初至走時數(shù)據(jù)建立近地表速度模型，近地表速度模型采用相同大小的正方形網(wǎng)格離散；

步驟三：計算地表平均高程；

步驟四：對近地表速度模型中地表以上的空氣層以地表速度進行填充；

步驟五：計算炮點到地表平均高程的地表一致性靜校正量；炮點到地表平均高程的地表一致性靜校正量計算公式為：

，（1）

式中，t_s1是炮點到地表平均高程的地表一致性靜校正量，i是炮點與地表平均高程間離散網(wǎng)格的序號，n₁是炮點與地表平均高程間離散網(wǎng)格的最大序號，dz是離散網(wǎng)格的厚度，v_i是第i個離散網(wǎng)格的速度；

步驟六：計算檢波點到地表平均高程的地表一致性靜校正量；檢波點到地表平均高程的地表一致性靜校正量計算公式為：

，（2）

式中，t_r1是檢波點到地表平均高程的地表一致性靜校正量，j是檢波點與地表平均高程間離散網(wǎng)格的序號，m₁是檢波點與地表平均高程間離散網(wǎng)格的最大序號，dz是離散網(wǎng)格的厚度，v_j是第j個離散網(wǎng)格的速度；

步驟七：對預處理后的地震記錄進行地表一致性靜校正；地表一致性靜校正的計算公式為：

，（3）

式中，d₁是地表一致性靜校正后的地震記錄，x_s是炮點x軸坐標，x_r是檢波點x軸坐標，t_sr¹是地表一致性靜校正后的反射波走時，d是預處理后的地震記錄，t_sr是地表一致性靜校正前的反射波走時，t_s1是炮點到地表平均高程的地表一致性靜校正量，t_r1是檢波點到地表平均高程的地表一致性靜校正量；

步驟八：將地表一致性靜校正后的地震資料從時間-空間域變換到τ-p域；采用下式進行變換：

，（4）

式中，m是變換后的地震資料，τ是截距時間，p是射線參數(shù)，d₁是地表一致性靜校正后的地震記錄，x_s是炮點x軸坐標，x_r是檢波點x軸坐標，t_sr¹是地表一致性靜校正后的反射波走時；

步驟九：計算射線參數(shù)對應的截距時間校正量；采用下式計算截距時間校正量：

，（5）

式中，Δτ是截距時間校正量，k是射線經(jīng)過的近地表速度模型離散網(wǎng)格序號，n_z是射線經(jīng)過的近地表速度模型離散網(wǎng)格序號的最大值，dz是離散網(wǎng)格的厚度，q_k是射線經(jīng)過第k個網(wǎng)格時的垂直慢度；

步驟十：對變換后的地震資料進行截距時間校正；采用下式進行截距時間校正：

，（6）

式中，m’是截距時間校正后的地震資料，τ’是校正后的截距時間，p是射線參數(shù)，m是變換后的地震資料，τ是截距時間，Δτ是截距時間校正量；

步驟十一：將截距時間校正后的地震資料從τ-p域變換到時間-空間域；采用下式進行變換：

，（7）

式中，d’是變換后時間-空間域的地震資料，x_s是炮點x軸坐標，x_r是檢波點x軸坐標，t_sr²是地表非一致性靜校正后的反射波走時，m’是截距時間校正后的地震資料，τ’是校正后的截距時間，p是射線參數(shù)。