[發明專利]一種地層巖石骨架的測井識別方法

有效

• 申請號: 201310218843.4
• 申請日: 2013-06-04
• 公開（公告）號: CN104213899B
• 公開（公告）日: 2017-08-04
• 發明（設計）人: 李浩;魏修平;王丹丹;馮瓊;陳萍
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油勘探開發研究院
• 主分類號: E21B47/00
• 分類號: E21B47/00
• 代理公司: 北京思創畢升專利事務所11218
• 代理人: 郭韞

1.一種地層巖石骨架的測井識別方法，其特征在于：所述方法包括：

尋找電阻率趨于無窮大時的深度點，將該深度點的測井值作為巖石骨架值；所述巖石骨架值包括聲波時差、補償密度及補償中子測值；

對于有巖心分析實驗的探區，進行巖石骨架值的檢驗分析，得到通過巖心實驗驗證了的巖石骨架值；

將所述通過巖心實驗驗證了的巖石骨架值應用到同巖性未知探區，

其中，所述巖石骨架值的檢驗分析包括：

(1)、將所述巖石骨架值與利用巖心分析孔隙度與聲波或密度交會得到的巖石骨架值做對比分析，得到兩者的一致性；

(2)、將所述巖石骨架值輸入孔隙度計算公式，求取地層孔隙度；在已知取心井中，根據巖心分析的孔隙度對所求得的地層孔隙度的計算精度進行再檢驗，

其中，所述根據巖心分析的孔隙度對所求得的地層孔隙度的計算精度進行再檢驗是這樣實現的：將計算得到的地層孔隙度與巖心分析孔隙度做相對誤差分析，如果相對誤差值小于最低參考值，則所述巖石骨架值即為通過巖心實驗驗證了的巖石骨架值；

其中，所述最低參考值取10％；

其中，將所述通過巖心實驗驗證了的巖石骨架值應用到同巖性未知探區是這樣實現的：將所述通過巖心實驗驗證了的巖石骨架值代入所述孔隙度計算公式，得到同巖性未知探區的孔隙度；所述同巖性未知探區是指與所述有巖心分析實驗的探區的巖性相同的井，該探區的聲波時差骨架值或者密度骨架值或者中子骨架值是未知的；

其中，所述電阻率趨于無窮大是指地層電阻率大于2000歐姆米；

其中，選取質量合格的測井曲線，讀取所述深度點的測井值作為所述巖石骨架值。

2.根據權利要求1所述的地層巖石骨架的測井識別方法，其特征在于：所述利用巖心分析孔隙度與聲波或密度交會得到的巖石骨架值是這樣實現的：

把巖心分析孔隙度按巖性進行分類：在同一種巖性下，根據巖心分析孔隙度的深度值，在測井數據中找到該深度下相對應的聲波測井值或密度測井值；

把巖心分析孔隙度值和所述聲波測井值做二維交會圖，然后用直線擬合后，得到y＝kx+b的關系，其中y指代聲波時差，x指代孔隙度，當孔隙度x為零時，b就是聲波時差骨架值；

同樣，把巖心分析孔隙度值和密度測井值做二維交會圖，然后用直線擬合后，得到y＝kx+b的關系，其中y指代密度，x指代孔隙度，當孔隙度x為零時，b就是密度骨架值；

如果有巖心分析的視密度，則直接用視密度替代測井曲線數據中得到的密度值。

3.根據權利要求2所述的地層巖石骨架的測井識別方法，其特征在于：所述孔隙度計算公式如下：

其中，

式中：為巖石孔隙度，小數；Δt為測量的巖石聲波時差，單位為μs/m；Δt_ma為巖石骨架聲波時差，單位為μs/m；Δt_f為巖石孔隙流體的聲波時差，單位為μs/m；為密度孔隙度，小數；ρ_b為測量巖石的體積密度，單位為g/cm³；ρ_ma為巖石骨架密度，單位為g/cm³；ρ_f為巖石孔隙流體的體積密度，單位為g/cm³。

4.根據權利要求1所述的地層巖石骨架的測井識別方法，其特征在于：所述計算精度是指利用巖石骨架值計算得到的地層孔隙度與巖心分析孔隙度一致性的高低。