[發明專利]一種分析儲層孔隙度的方法

在審

• 申請號: 201410058847.5
• 申請日: 2014-02-20
• 公開（公告）號: CN104863572A
• 公開（公告）日: 2015-08-26
• 發明（設計）人: 張軍;李軍;武清釗;蘇俊磊;南澤宇;王曉暢
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油勘探開發研究院
• 主分類號: E21B47/08
• 分類號: E21B47/08;E21B49/00
• 代理公司: 北京聿宏知識產權代理有限公司 11372
• 代理人: 吳大建;劉華聯

1.一種分析儲層孔隙度的方法，其包括：

步驟A，通過測井資料獲取儲層的測井參數；

步驟B，測量儲層中n個深度的巖心分析孔隙度；

步驟C，建立儲層的巖心分析孔隙度與測井參數的對應關系；

步驟D，利用步驟C中所建立的儲層的巖心分析孔隙度與測井參數的對應關系通過測井參數來確定待解釋儲層的孔隙度。

2.根據權利要求1所述的方法，其特征在于，所述測井參數包括測井實測值和/或測井計算值。

3.根據權利要求2所述的方法，其特征在于，所述測井實測值包括聲波時差測井值、密度測井值、中子測井值中的至少一個。

4.根據權利要求2所述的方法，其特征在于，所述測井計算值由聲波時差測井值、密度測井值、中子測井值中的至少一個計算得到。

5.根據權利要求1-4中任意一項所述的方法，其特征在于，步驟C包括：

步驟K，以測井參數為橫坐標，以巖心分析孔隙度為縱坐標，建立平面直角坐標系，并將n對巖心分析孔隙度和相同深度的測井參數所對應的n個數據點繪制于坐標系內，并且所述數據點的坐標為（x_j，y_j），

其中，x_j為測井參數值；

y_j為巖心分析孔隙度；

j=1、2、…、n-i；

i=0、1、2、…n-1；

步驟L，對n-i對巖心分析孔隙度和相同深度的測井參數所對應的n-i個數據點進行線性擬合獲得如式（Ⅰ）所示的擬合直線；

y=a_i·x+b_i??????????????式（Ⅰ）

其中，ai為所述擬合直線的斜率，bi為所述擬合直線的截距，ai和bi均為步驟L直線擬合過程中直接獲得；

步驟M，通過式（Ⅱ）計算各數據點與步驟L中的擬合直線的距離；

dj=|ai·xj-yj+bi|(ai2+1)]]>???????????式(II)

步驟N，依據各數據點與步驟L中的擬合直線的距離判斷是否存在異常點，并確定儲層的巖心分析孔隙度與測井參數的對應關系。

6.根據權利要求5所述的方法，其特征在于，在步驟N中，依據各數據點與步驟L中的擬合直線的距離對n-i個數據點進行排序，選出距離排序最大的數據點和距離排序中部的數據點，通過式（Ⅲ）計算距離排序最大的數據點和距離排序中部的數據點的距離倍數，然后將其與距離倍數界限值進行比較，并且，當所計算的距離倍數大于距離倍數界限值時，則判斷該距離排序最大的數據點為異常點，剔除該數據點，i值加1，重復步驟L到N，直到所計算的距離倍數大于距離倍數界限值，則判斷為無異常數據點，步驟L中所擬合的直線即為儲層的巖心分析孔隙度與測井參數的對應關系；

l=lenmax/len50??????????????式（Ⅲ）

其中，l為距離最大的數據點的距離倍數；

lenmax為距離最大的數據點的距離；

len50為排序在中部的數據點的距離，該數據點的排序為（n-i）/2。

7.根據權利要求6所述的方法，其特征在于，所述距離倍數界限值為2-5。