[發明專利]適合于致密砂巖地層的井震頻散校正方法

在審

• 申請號: 202110509081.8
• 申請日: 2021-05-11
• 公開（公告）號: CN115327640A
• 公開（公告）日: 2022-11-11
• 發明（設計）人: 張云濤;羅霞;李民龍;張云銀;劉秉建;張學芳;樊明;王雨潔;伊慧
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勝利油田分公司物探研究院
• 主分類號: G01V1/50
• 分類號: G01V1/50
• 代理公司: 濟南日新專利代理事務所(普通合伙) 37224
• 代理人: 崔曉艷

1.適合于致密砂巖地層的井震頻散校正方法，其特征在于，該適合于致密砂巖地層的井震頻散校正方法包括：

步驟1，選取目標儲層巖心，開展巖性、物性以及低頻巖石物理參數測量；

步驟2，構建全頻段致密砂巖巖石物理模型；

步驟3，根據實驗數據標定動態等效微觀尺度模型；

步驟4，進行井震速度頻散校正。

2.根據權利要求1所述的適合于致密砂巖地層的井震頻散校正方法，其特征在于，在步驟1中，選取目標地層典型塊巖心開展孔隙度、滲透率、密度、礦物組分這些參數測量，而后開展地層條件下低頻巖石物理參數測量，飽和狀態為干燥、飽水、飽油，得到巖心的孔隙度、滲透率、密度、礦物組分數據、以及地層條件下巖心的低頻縱橫波速度數據。

3.根據權利要求2所述的適合于致密砂巖地層的井震頻散校正方法，其特征在于，在步驟1中，開展地層條件下低頻巖石物理參數測量時，測量頻率為12-200Hz。

4.根據權利要求1所述的適合于致密砂巖地層的井震頻散校正方法，其特征在于，在步驟2中，動態等效微觀尺度模型的等效介質的體積模量K_eff和剪切模量u_eff可表示為：

其中，ω為角頻率，K為骨架體積模量、μ為剪切模量，ε為裂隙密度，為孔隙度，λ為拉梅系數，K_f為流體體模量，Kc、Kp、A、B均為中間參數沒有具體物理意義，A和B分別為

其中，

τ為時間尺度，控制著頻散特征出現的頻帶范圍，它與流體的粘度η成正比，與滲透率k成反比，并與微裂縫半徑a有關，其倒數稱為特征頻率；對于較小的縱橫比：

為顆粒尺度，v為固體顆粒的泊松比。

5.根據權利要求4所述的適合于致密砂巖地層的井震頻散校正方法，其特征在于，在步驟2中，根據Connell和Budiansky提出的公式，得到P波速度V_P和S波的速度V_S分別為：

式中，Re為求實部，ρ為密度。

6.根據權利要求1所述的適合于致密砂巖地層的井震頻散校正方法，其特征在于，在步驟3中，首先通過步驟1中得到的孔隙度、滲透率、密度、礦物組分、含水飽和度這些數據，基于動態等效微觀尺度模型計算全頻段縱橫波速度數據，而后調整裂隙密度與實驗室測量數據對比，驗證模型計算的精度。

7.根據權利要求1所述的適合于致密砂巖地層的井震頻散校正方法，其特征在于，在步驟4中，首先通過測井解釋成果數據獲得巖性組分數據，而后基于Voigt-Reuss-Hill模型計算骨架模量K和μ，繼而通過測井解釋成果數據獲得孔隙度數據，假設孔隙度全由硬孔隙組成，基于Berrymann模型計算干巖石彈性模量，然后通過測井解釋成果數據獲得流體組分數據，滲透率數據，基于動態等效微觀尺度模型計算測井頻段下縱波速度。

8.根據權利要求7所述的適合于致密砂巖地層的井震頻散校正方法，其特征在于，在步驟4中，如果計算縱波速度和實測縱波速度的差值小于誤差，那么這時可以直接計算得到地震頻段縱橫波速度，如果計算縱波速度和實測縱波速度的差值大于誤差，那么加入裂隙密度后再次按照以上流程進行計算，反復循環，直到計算縱波速度和實測縱波速度差值小于誤差。

9.根據權利要求8所述的適合于致密砂巖地層的井震頻散校正方法，其特征在于，步驟4包括：

4a，進行骨架模量計算；

4b，進行干巖石的彈性模量計算；

4c，進行含流體巖石全頻段縱橫波速度計算。