[發明專利]一種基于常規測井的致密砂礫巖儲層孔隙結構評價方法

有效

• 申請號: 202110389007.7
• 申請日: 2021-04-12
• 公開（公告）號: CN113065286B
• 公開（公告）日: 2022-03-18
• 發明（設計）人: 李瑋;劉向君;施雷庭;張玉龍
• 申請（專利權）人: 西南石油大學
• 主分類號: G06F30/27
• 分類號: G06F30/27;G06K9/62;E21B49/00
• 代理公司: 北京中索知識產權代理有限公司 11640
• 代理人: 唐亭

1.一種基于常規測井的致密砂礫巖儲層孔隙結構評價方法，其特征在于，步驟如下：

S1、交會分析

將9條常規測井曲線兩兩交會，共計36種交會情況，得到36張交會圖；采用中心距離法計算每張交會圖上三類孔隙結構數據點群的總類間距；在36張交會圖中選取總類間距較大的多張交會圖，將選取的交會圖中涉及的測井參數作為敏感參數帶入后續步驟；

S2、高維輸入空間降維

基于所選取的測井參數作為輸入，采用主成分分析方法得到能夠表征高維輸入空間的第一、第二主成分，實現降維；

S3、明確聚類中心

首先對各類孔隙結構數據點群計算其初始聚類中心，并計算各點距離其聚類中心的歐式距離；將距離最遠的點剔除并重新計算聚類中心，重復此過程，直至聚類中心的移動量小于首次剔出點后移動量的十分之一時，停止此過程；此時對應的聚類中心為最終聚類中心；依次完成各類點集合的聚類中心計算；

S4、常規測井孔隙結構劃分

將實際地層常規測井數據中的敏感參數值通過主成分分析方法得到第一主成分值F1和第二主成分值F2，然后將得到的F1、F2投影到第一、第二主成分交會圖中，然后計算并比較F1、F2對應的點與各聚類中心之間的距離，該點距離哪類點的聚類中心最近就屬于哪一類。

2.如權利要求1所述的基于常規測井的致密砂礫巖儲層孔隙結構評價方法，其特征在于，所述9條常規測井曲線為井徑、自然伽馬、自然電位、補償密度、補償中子、補償聲波、深側向電阻率、淺側向電阻率、微球形聚焦電阻率。

3.如權利要求1所述的基于常規測井的致密砂礫巖儲層孔隙結構評價方法，其特征在于，所述步驟S1中，計算每張交會圖上三類孔隙結構數據點群的總類間距，具體方法是：首先對每類孔隙結構數據點群做平均，得到該類數據點群的平均點作為中心點；然后分別計算各類數據點群的中心點之間的歐式距離作為各類數據點群之間類間距，求和得到總類間距。

4.如權利要求3所述的基于常規測井的致密砂礫巖儲層孔隙結構評價方法，其特征在于，步驟S1中，總類間距的計算過程如下：

首先，計算各類數據點群的中心點之間的歐式距離，公式如下：

式中，i、j＝1，2，3為類標號，D_ij為i類孔隙結構點群的中心點與j類孔隙結構點群的中心點之間的距離；x_i和x_j為中心點的橫坐標，x_max和x_min分別為橫坐標的最大值和最小值；y_i和y_j為中心點的縱坐標，y_max和y_min分別為縱坐標的最大值和最小值；

總類間距的計算公式為：

D＝D₁₂+D₂₃+D₁₃

式中，D為總類間距；D₁₂是一類孔隙結構點群的中心點與二類孔隙結構點群的中心點之間的距離；D₁₃是一類孔隙結構點群的中心點與三類孔隙結構點群的中心點之間的距離；D₂₃是二類孔隙結構點群的中心點與三類孔隙結構點群的中心點之間的距離。

5.如權利要求3所述的基于常規測井的致密砂礫巖儲層孔隙結構評價方法，其特征在于，所述步驟S1中，計算出每張交會圖的總類間距后，按照總類間距由大到小的順序將36張交會圖排序，選取總類間距最大的前4-5張交會圖。

6.如權利要求1所述的基于常規測井的致密砂礫巖儲層孔隙結構評價方法，其特征在于，所述步驟S3中，首先對主成分交會圖上各類孔隙結構點群求平均值，得到初始聚類中心。