CN102052074A - 一种岩心数字化确定岩电参数的方法 - Google Patents

Abstract

本发明涉一种岩心数字化确定岩电参数m、n值的方法；对取岩心采用核磁快速测量及计算得到核磁T₂谱和孔隙度

岩心数字化后的数据体；分别通过公式

、

计算胶结指数m、饱和度指数n；利用建立的胶结指数

、饱和度指数

的函数关系式，分别代入未取心井段的测井计算地层孔隙度，得到未取心井段的胶结指数m值和未取心井段的饱和度指数n值；本发明能够较为迅速地求取m、n值，能够达到与实验室常规标准流程实验及计算获取的实验测试效果，克服了实验室在较短时间内无法获得岩电参数m、n值的难题，求取方便快捷，能够及时服务或满足实际生产的需要，具有较好的应用效果。

Description

一种岩心数字化确定岩电参数的方法

技术领域

本发明涉及石油测井技术领域，具体是一种岩心数字化后确定岩电参数m、n值的方法。

背景技术

在石油勘探和测井技术领域，通常储层岩电参数胶结指数m、饱和度指数n值的求取是通过选取能够代表本地区地质特征的系列岩心，按照《岩心分析方法(SY/T 5336-2006)》标准流程进行实验及计算，得到岩心的孔隙度φ，其中，岩心孔隙度φ为小数；按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算，得到岩心的胶结指数m值、饱和度指数n值，岩性系数a、b值。

岩电实验首先要求用地层水完全饱和岩心，然而要使岩心完全达到饱和或达到不同程度的含水饱和度，往往需要较长的时间周期，对于低孔隙度低渗透率岩心，岩电实验需要的实验周期会更长，由于受生产时间和条件的限制，岩电参数m、n值在较短时间内无法按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算得到，导致岩电参数m、n值无法及时服务和满足生产实际的需要。

发明内容

本发明能够提供一种岩心数字化后确定岩电参数m、n值的方法。

本发明采用下述步骤实现：

1)旋转式井壁取心或钻井取心后钻取岩心；

2)对步骤1)所述的系列岩心，采用核磁进行现场快速实验测量及计算得到核磁T₂谱及其孔隙度

等岩心数字化后的数据体；

3)对步骤1)所述的系列岩心，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算，分别得到胶结指数m、饱和度指数n值；

4)利用步骤2)建立的岩心数字化数据体以及步骤3)得到系列岩心的m、n值，分别通过公式(1)、公式(2)计算胶结指数m、饱和度指数n；

$m={\mathrm{\α}}_3\×{{\mathrm{\φ}}_{T_2}}^3+{\mathrm{\α}}_2\×{{\mathrm{\φ}}_{T_2}}^2+{\mathrm{\α}}_1\×{\mathrm{\φ}}_{T_2}+{\mathrm{\α}}_0---\left(1\right)$

$n={\mathrm{\β}}_1\×{\mathrm{Log}}_e\left({\mathrm{\φ}}_{T_2}\right)+{\mathrm{\β}}_0---\left(2\right)$

其中：

α₀、α₁、α₂、α₃为待定系数，利用最小二乘法拟合得到；β₀、β₁为待定系数，利用最小二乘法拟合得到；

5)利用步骤3)建立的胶结指数

饱和度指数

的函数关系式，分别代入未取心井段的测井计算地层孔隙度，即可得到未取心井段的胶结指数m值和未取心井段的饱和度指数n值。

本发明能够较为迅速地求取m、n值，能够达到与实验室常规标准流程实验及计算获取的实验测试效果，克服了实验室在较短时间内无法获得岩电参数m、n值的难题，求取方便快捷，能够及时服务或满足实际生产的需要，具有较好的应用效果。

附图说明

图1本发明函数公式计算m值与实验测定m值对比图；

图2本发明函数公式计算n值与实验测定n值对比图。

图3本发明函数公式计算求取的未取心井段的岩电参数m、n值曲线图。

具体实施方式

实施例：

1)利用旋转式井壁取心得到146块岩心样品；

2)对步骤1)所述的系列岩心，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算，分别得到胶结指数m、饱和度指数n值；

3)对步骤1)所述的系列岩心，采用核磁进行现场快速实验测量及计算得到岩心数字化后的核磁T₂谱及其孔隙度

4)对步骤2)得到的胶结指数m、对步骤3)得到的岩心核磁孔隙度

本发明胶结指数m计算公式

$m={\mathrm{\α}}_3\×{{\mathrm{\φ}}_{T_2}}^3+{\mathrm{\α}}_2\×{{\mathrm{\φ}}_{T_2}}^2+{\mathrm{\α}}_1\×{\mathrm{\φ}}_{T_2}+{\mathrm{\α}}_0---\left(1\right)$

利用最小二乘法拟合得到待定系数α₃＝-0.0002、α₂＝0.0064、α₁＝-0.0439、α₀＝1.5193。

5)对步骤2)得到的饱和度指数n、对步骤3)得到的岩心孔隙度

本发明饱和度指数n计算公式

$n={\mathrm{\β}}_1\×{\mathrm{Log}}_e\left({\mathrm{\φ}}_{T_2}\right)+{\mathrm{\β}}_0---\left(2\right)$

利用最小二乘法拟合得到待定系数β₁＝-1.1378、β₀＝5.0032。

6)利用步骤4)得到的胶结指数m的函数关系式，见公式(3)，代入未取心井段的测井计算地层孔隙度，即可计算得到未取心井段的胶结指数m值。

$m=-0.0002\×{{\mathrm{\φ}}_{T_2}}^3+0.0064\×{{\mathrm{\φ}}_{T_2}}^2-0.0439\×{\mathrm{\φ}}_{T_2}+1.5193---\left(3\right)$

7)利用步骤5)得到的饱和度指数n的函数关系式，见公式(4)，代入未取心井段的测井计算地层孔隙度，即可计算得到未取心井段的饱和度指数n值。

$n=-1.1378\×{\mathrm{Log}}_e\left({\mathrm{\φ}}_{T_2}\right)+5.0032---\left(4\right)$

利用本发明能够达到与实验室岩心按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算获取的测试效果，而且求取方便快捷，能够及时服务或满足实际生产的需要，具有较好的应用效果。

利用本发明计算m值与岩心测定m值对比图见图1。

利用本发明计算n值与岩心测定n值对比图见图2。

利用本发明计算求取的未取心井段的岩电参数m、n值曲线图见图3。

图1是利用本发明计算m值与实验室按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算得到m值对比图，从图1可以看出，采用本发明计算m值与岩心实际测定m值均分布于45度线附近，应用效果较好。

图2是利用本发明计算n值与实验室按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算得到n值对比图，从图2可以看出，采用本发明计算n值与岩心实际测定n值均分布于45度线附近，应用效果明显。

图3是利用本发明计算求取的未取心井段的岩电参数m、n值曲线图，利用本发明计算求取的m、n值可用于地层含水饱和度的计算，能够达到与实验室按照《岩石电阻率参数实验室测量及计算方法(SY/T5385-1991)》标准流程进行实验及计算得到m、n值的效果，应用效果良好。

Claims (1)

1.一种岩心数字化确定岩电参数的方法，其特征在于：该方法采用下述步骤：

1)旋转式井壁取心或钻井取心后钻取岩心；

2)对步骤1)所述的系列岩心，采用核磁进行现场快速实验测量及计算得到核磁T₂谱及其孔隙度岩心数字化后的数据体；

3)对步骤1)所述的系列岩心，按照《岩石电阻率参数实验室测量及计算方法(SY/T 5385-1991)》标准流程进行实验及计算，分别得到胶结指数m、饱和度指数n值；

4)利用步骤2)建立的岩心数字化数据体以及步骤3)得到系列岩心的m、n值，分别通过公式(1)、公式(2)计算胶结指数m、饱和度指数n；

$m={\mathrm{\α}}_3\×{{\mathrm{\φ}}_{T_2}}^3+{\mathrm{\α}}_2\×{{\mathrm{\φ}}_{T_2}}^2+{\mathrm{\α}}_1\×{\mathrm{\φ}}_{T_2}+{\mathrm{\α}}_0---\left(1\right)$

$n={\mathrm{\β}}_1\×{\mathrm{Log}}_e\left({\mathrm{\φ}}_{T_2}\right)+{\mathrm{\β}}_0---\left(2\right)$

其中：

α₀、α₁、α₂、α₃为待定系数，利用最小二乘法拟合得到；

β₀、β₁为待定系数，利用最小二乘法拟合得到；

5)利用步骤4)建立的胶结指数

饱和度指数

的函数关系式，分别代入未取心井段的测井计算地层孔隙度，即可得到未取心井段的胶结指数m值和未取心井段的饱和度指数n值。

Images