[發明專利]基于方位彈性阻抗傅里葉級數展開的裂縫預測方法及系統

有效

• 申請號: 201811024016.0
• 申請日: 2018-09-04
• 公開（公告）號: CN110873897B
• 公開（公告）日: 2021-07-27
• 發明（設計）人: 郭旭升;劉曉晶;陳超;王威;彭嫦姿;肖秋紅;石美璟;劉苗苗
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勘探分公司
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京思創大成知識產權代理有限公司 11614
• 代理人: 董曉盈

1.一種基于方位彈性阻抗傅里葉級數展開的裂縫預測方法，包括：

根據原始道集，獲得包含方位角信息的角道集；

根據方位角與入射角將所述角道集進行角度劃分，獲得不同方位角、不同入射角的部分疊加地震數據；

根據測井的縱波速度、橫波速度、密度、各向異性信息，計算與所述部分疊加地震數據對應的每一個方位角、每一個入射角的方位彈性阻抗曲線；

根據每一個方位角、每一個入射角的方位彈性阻抗曲線開展正演，與對應的部分疊加地震數據進行匹配標定，提取每一個方位角、每一個入射角對應的地震子波，進而建立每一個方位角、每一個入射角對應的彈性阻抗模型；

根據所述部分疊加地震資料、所述彈性阻抗模型、所述地震子波，計算每一個方位角、每一個入射角對應的方位彈性阻抗體，并對方位彈性阻抗體進行傅里葉級數展開，實現裂縫預測。

2.根據權利要求1所述的基于方位彈性阻抗傅里葉級數展開的裂縫預測方法，其中，通過公式(1)計算方位彈性阻抗曲線：

其中，EI(θ，φ)為方位彈性阻抗曲線，g為橫波速度與縱波速度的比值的平方，V_p為縱波速度，V_s為橫波速度，ρ為密度，ε^(v)、δ^(v)、γ分別表示為Thomsen弱各向異性參數，θ表示地震波入射角，φ表示為測線方向與裂縫傾向的夾角，V_p0、V_s0、ρ₀分別表示為縱波速度、橫波速度、密度的均值。

3.根據權利要求1所述的基于方位彈性阻抗傅里葉級數展開的裂縫預測方法，其中，通過公式(2)計算方位彈性阻抗體：

其中，F(m)為目標函數，G為正演算子，m為模型參數向量，d為觀測數據向量，i＝1,2，…，N，N為采樣點數量，X_n為地震數據的協方差矩陣，C為積分矩陣，ξ為相對彈性阻抗，表示第i組模型參數向量，表示第i組模型參數向量的協方差矩陣，X_Cm-ξ表示模型誤差的協方差矩陣。

4.根據權利要求1所述的基于方位彈性阻抗傅里葉級數展開的裂縫預測方法，其中，通過公式(3)進行傅里葉級數展開：

ln(EI_A(θ，φ))＝A₀+A₂cos2φ+A₄cos4φ (3)

其中，A₀、A₂與A₄分別為零階、二階、四階傅里葉系數，表達式分別為：

其中，EI_A(θ,φ)表示歸一化后的方位彈性阻抗，g為橫波速度與縱波速度的比值的平方，V_p為縱波速度，V_s為橫波速度，ρ為密度，ε^(v)、δ^(v)、分別表示為Thomsen弱各向異性參數，φ表示為測線方向與裂縫傾向的夾角，θ表示地震波入射角。

5.根據權利要求1所述的基于方位彈性阻抗傅里葉級數展開的裂縫預測方法，其中，所述部分疊加地震數據的方位角與入射角分別為劃分區域的方位角與入射角的平均值。