[發明專利]一種利用地震資料預測鉆井漏失的方法

在審

• 申請號: 201811117109.8
• 申請日: 2018-09-25
• 公開（公告）號: CN110941010A
• 公開（公告）日: 2020-03-31
• 發明（設計）人: 齊晴;孫振濤;田建華;陳勇;楊勤林;曹少蕾;董清源;朱博華;梁碩博
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京聿宏知識產權代理有限公司 11372
• 代理人: 吳大建;張杰

1.一種利用地震資料預測鉆井漏失的方法，其特征在于，包括以下步驟：

S1：對原始地震數據體進行降噪濾波處理，得到濾波數據體；

S2：對所述濾波數據體進行相干處理得到相干體；

S3：將所述相干體與所述原始地震數據體進行融合得到融合數據體；

S4：根據所述融合數據體的剖面確定相干體的分布區域，根據所述相干體的分布區域預測裂縫發育帶和/或溶洞發育帶的位置；

S5：根據所述裂縫發育帶和/或溶洞發育帶的位置預測鉆井漏失層段。

2.根據權利要求1所述的利用地震資料預測鉆井漏失的方法，其特征在于，所述降噪濾波處理包括以下步驟：

對原始地震數據體進行高斯平滑處理。

3.根據權利要求2所述的利用地震資料預測鉆井漏失的方法，其特征在于，對原始地震數據體進行高斯平滑處理具體包括：

對原始地震數據體中當前采樣點的數據u_σ進行x、y方向上的偏導計算，公式如下：

根據高斯平滑公式分別對原始地震數據體中當前采樣點的數據u_σ在x、y方向上的偏導進行尺度為ρ的高斯平滑，得到結構矩陣的結構張量J_ρ，所述高斯平滑公式為u_σ⁰＝u_σ*G_ρ；其中，G_ρ為高斯核參數，表示以G_ρ為參數的高斯核；

所述結構矩陣的結構張量J_ρ為：

其中，表示原始地震數據體中當前采樣點的數據的梯度，j₁₁表示u_σ在x方向上的偏導進行尺度為ρ的高斯平滑后的結果，j₁₂表示u_σ在x、y方向上的偏導進行尺度為ρ的高斯平滑后的結果，j₂₂表示u_σ在y方向上的偏導進行尺度為ρ的高斯平滑后的結果。

4.根據權利要求3所述的利用地震資料預測鉆井漏失的方法，其特征在于，所述降噪濾波處理還包括：根據矩陣特征分解定理對所述結構張量J_ρ進行分解得到

其中，λ₁、λ₂為J_ρ的特征值，λ₁、λ₂對應的正交歸一化特征向量為：

其中，

v₁表示平行于梯度的方向，變化率為v₂表示垂直于梯度的方向，變化率為θ為梯度方向的角度；

根據和判斷原始地震數據體當前采樣點根據結構張量提取的局部特征，其擴散系數μ₁和μ₂的表達式為：

其中，α為常數；

根據矩陣特征分解定理和擴散系數得到構造擴散張量，其表達式為：

其中，a、b、c分別為構造擴散張量D的分量，具體表示為：

5.根據權利要求4所述的利用地震資料預測鉆井漏失的方法，其特征在于，所述降噪濾波處理還包括：

根據所述構造擴散張量D得到代表擴散濾波結果的張量擴散方程：

t代表擴散時間，代表當前采樣點的數據的散度；

在正交坐標系(x,y)中，對張量擴散方程進行離散化，得到濾波迭代公式：

u^k和u^k+1分別表示原始地震數據體在kΔt和(k+1)Δt時刻的濾波結果，Δt為迭代一次的擴散時間，k代表采樣點的個數；

通過濾波迭代公式得到濾波數據體。