[發明專利]含氣飽和度確定方法及裝置

有效

• 申請號: 201611263196.9
• 申請日: 2016-12-30
• 公開（公告）號: CN106842303B
• 公開（公告）日: 2019-07-05
• 發明（設計）人: 謝偉彪;司兆偉;莊東志;徐風;張建林;陳晶瑩;馬明明;殷秋麗
• 申請（專利權）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30
• 代理公司: 北京三高永信知識產權代理有限責任公司 11138
• 代理人: 董亞軍

1.一種含氣飽和度確定方法，其特征在于，所述方法包括：

分別獲取縱波和橫波在至少兩個巖樣中傳播的至少八個樣本波速參數，所述至少兩個巖樣中的每個巖樣對應至少兩個含氣飽和度，所述至少八個樣本波速參數包括至少四個橫波樣本參數和至少四個縱波樣本參數，每個巖樣對應至少兩個橫波樣本參數和至少兩個縱波樣本參數；

采用波速-時差轉換公式對所述至少八個樣本波速參數中的每個樣本波速參數進行轉換，得到至少八個樣本時差，所述波速-時差轉換公式為：

其中，DT為樣本時差，v為樣本波速參數；

根據所述至少四個縱波樣本時差，采用模型參數計算公式得到所述模型參數，所述模型參數計算公式為：

其中，ΔDT為所述模型參數，DTC為預設巖樣處于不同含水飽和度時對應的縱波時差，DTC_sg＝0為所述預設巖樣為第一巖樣時對應的第一縱波樣本時差，DTC_sg＝1為所述預設巖樣為第二巖樣時對應的第二縱波樣本時差，所述預設巖樣為所述至少一個巖樣中的任意一個巖樣；

根據每個巖樣對應的模型參數和每個巖樣對應的含氣飽和度，采用最小二乘法確定所述含氣飽和度計算公式，所述含氣飽和度計算公式為：

S_g∝ΔDT

其中，S_g為含氣飽和度，ΔDT為模型參數；

分別獲取橫波和縱波在待確定氣藏儲層中的傳播速度的參數，所述傳播速度的參數包括橫波速度參數和縱波速度參數；

根據所述橫波速度參數以及預先確定的轉換參數，分別確定第一縱波時差和第二縱波時差，所述第一縱波時差是縱波在含氣飽和度為0的待確定氣藏儲層中傳播時的縱波時差，所述第二縱波時差是縱波在含氣飽和度為1的待確定氣藏儲層中傳播時的縱波時差；

根據所述縱波速度參數、所述第一縱波時差、所述第二縱波時差和預先確定的模型參數，確定所述待確定氣藏儲層的含氣飽和度；

所述根據所述縱波速度參數、所述第一縱波時差、所述第二縱波時差和預先確定的模型參數，確定所述待確定氣藏儲層的含氣飽和度，包括：

將所述縱波速度參數、所述第一縱波時差、所述第二縱波時差和所述模型參數代入預先確定的含氣飽和度計算公式，得到所述待確定氣藏儲層的含氣飽和度。

2.根據權利要求1所述的方法，其特征在于，所述至少八個樣本時差包括與所述至少四個橫波樣本參數對應的至少四個橫波樣本時差，以及與所述至少四個縱波樣本參數對應的至少四個縱波樣本時差；

在所述采用波速-時差轉換公式對所述至少八個樣本波速參數中的每個樣本波速參數進行轉換之后，所述方法還包括：

根據所述至少四個橫波樣本時差和所述至少四個縱波樣本時差，采用最小二乘法確定所述轉換參數。

3.根據權利要求1所述的方法，其特征在于，每個巖樣包括含氣飽和度為0時的第一巖樣以及含氣飽和度為1時的第二巖樣，所述至少四個縱波樣本參數包括縱波在所述第一巖樣中傳播的至少兩個第一縱波樣本參數，以及在所述第二巖樣中傳播的至少兩個第二縱波樣本參數，所述至少八個樣本時差包括與所述至少四個縱波樣本參數對應的至少四個縱波樣本時差，所述至少四個縱波樣本時差包括分別與所述至少兩個第一縱波樣本參數對應的至少兩個第一縱波樣本時差，以及與所述至少兩個第二縱波樣本參數對應的至少兩個第二縱波樣本時差。

4.根據權利要求3所述的方法，其特征在于，在所述分別獲取縱波和橫波在至少一個巖樣中傳播的至少八個樣本波速參數之前，所述方法還包括：

獲取所述至少一個巖樣中每個巖樣對應的含水飽和度；

根據每個巖樣對應的含水飽和度，采用飽和度轉換公式計算得到每個巖樣的含氣飽和度，所述飽和度轉換公式為：

S_g＝1-S_w

其中，S_g為含氣飽和度，S_w為含水飽和度。