一种储层孔隙孔径分布确定方法和装置制造方法及图纸

本文提供一种储层孔隙孔径分布确定方法和装置，方法包括获取储层样品；对储层样品进行核磁共振得到T2分布谱；对储层样品进行孔隙提取得到多个孔隙，并根据各个孔隙的形状参数计算得到储层样品的综合形状因子，形状参数至少包括孔隙的最大直径和面积；对储层样品进行X射线衍射分析确定储层样品的表面弛豫率；根据综合形状因子和表面弛豫率转换T2分布谱，得到储层样品的孔径分布。本文提供的方法能够根据孔隙的最大直径和面积计算得到储层样品的综合形状因子以及根据X射线衍射结果得到储层样品的表面弛豫率，以对T2分布谱进行转换得到孔隙孔径分布，提高了所得到的孔径分布的准确性，且极大地降低了成本，为储层油气资源的勘探与开发提供支撑。探与开发提供支撑。探与开发提供支撑。

【技术实现步骤摘要】
一种储层孔隙孔径分布确定方法和装置


[0001]本文涉及石油天然气勘探
，尤其是一种储层孔隙孔径分布确定方法和装置。

技术介绍

[0002]孔隙结构是制约油气勘探的重要因素，其不仅影响着储层储集油气的能力，而且影响着油气等流体在储层中流动的难易程度。近年来，对致密油气和页岩油气的勘探越来越广泛，孔隙结构对致密油气和页岩油气这些储层的勘探开发的影响也更加明显，因此储层孔隙结构的定量表征对油气勘探和开发至关重要。
[0003]在目前的对储层孔隙结构的分析方法中，核磁共振技术是其中一种非常重要的分析方法。核磁共振技术可以在不损坏样品的基础上，通过获取样品中充满的流体的核磁信号进而获取样品的孔隙结构，因此，相对于其他的储层孔隙结构的分析方法，核磁共振技术具有无损的优势。由于核磁共振技术分析储层孔隙结构主要是通过将饱含水样品的核磁T2谱分布转化为孔隙的孔径分布曲线进而定量评价孔隙结构，其转化过程如下式所示：
[0004]r＝C
×
T2；
[0005]其中，r为孔隙直径，T2为表面弛豫时间，其通过核磁共振实验获取；C为转换系数。可以看出，若要将核磁共振实验的T2谱分布准确转化为储层孔隙的孔径分布，就需要对C的取值进行准确的判定。但现有的测定转换系数的方法，包括根据储层表面弛豫率和孔隙空间形状来计算的方法、以及直接根据储层岩性的不同对其进行赋值的方法，均存在着准确性较低的问题，因此，导致现有的核磁共振技术仍无法实现对储层孔隙结构的准确分析，制约了页岩油气和致密砂岩等储层的油气的深入勘探开发。
[0006]有鉴于此，本文旨在提供一种储层孔隙孔径分布确定方法和装置。

技术实现思路

[0007]针对现有技术的上述问题，本文的目的在于，提供一种储层孔隙孔径分布确定方法和装置，以解决现有技术中无法准确获得转化系数导致对储层孔隙孔径分布分析判定不准确的问题。
[0008]为了解决上述技术问题，本文的具体技术方案如下：
[0009]第一方面，本文提供一种储层孔隙孔径分布确定方法，所述方法包括：
[0010]获取储层样品；
[0011]对所述储层样品进行核磁共振，得到所述储层样品的T2分布谱；
[0012]对所述储层样品进行孔隙提取得到多个孔隙，并根据各个孔隙的形状参数计算得到所述储层样品的综合形状因子，所述形状参数至少包括孔隙的最大直径和面积；
[0013]对所述储层样品进行X射线衍射分析，并根据X射线衍射分析结果确定所述储层样品的表面弛豫率；
[0014]根据所述综合形状因子和所述表面弛豫率转换所述T2分布谱，得到所述储层样品
的孔径分布。
[0015]进一步地，对所述储层样品进行孔隙提取得到多个孔隙，并根据各个孔隙的形状参数计算得到所述储层样品的综合形状因子，包括：
[0016]通过扫描电镜获取所述储层样品的图像；
[0017]根据预设的灰度阈值识别出所述图像中的多个孔隙及其形状参数；
[0018]根据各个孔隙的形状参数，计算得到各个孔隙的形状因子；
[0019]根据各个孔隙的形状因子和各个孔隙的面积，计算得到所述储层样品的综合形状因子。
[0020]具体地，所述各个孔隙的形状因子通过如下公式计算得到：
[0021][0022]其中，F
i
为第i个孔隙的形状因子；D
imax
为第i个孔隙的最大直径；S
i
为第i个孔隙的面积。
[0023]优选地，所述形状参数还包括孔隙的最小直径和周长，所述各个孔隙的形状因子通过如下公式计算得到：
[0024][0025]其中，F
i
为第i个孔隙的形状因子；D
imax
为第i个孔隙的最大直径；D
imin
为第i个孔隙的最小直径；L
i
为第i个孔隙的周长；S
i
为第i个孔隙的面积。
[0026]具体地，所述储层样品的综合形状因子通过如下公式计算得到：
[0027][0028]其中，F为综合形状因子；F
i
为第i个孔隙的形状因子；S
i
为第i个孔隙的面积；i的取值范围为1至n，n为识别出的孔隙的数量。
[0029]优选地，在根据各个孔隙的形状因子和各个孔隙的面积，计算得到所述储层样品的综合形状因子之前，所述方法还包括：
[0030]根据预设的面积阈值对识别出的孔隙进行筛选，去除孔隙面积小于所述面积阈值的孔隙。
[0031]具体地，根据X射线衍射分析结果确定所述储层样品的表面弛豫率，包括：
[0032]根据X射线衍射结果，确定所述储层样品中各矿物的含量，所述矿物至少包括黄铁矿、石英、钾长石、斜长石、尖晶石、碳酸盐岩矿物、黏土含量和菱铁矿中的一种或几种的组合；
[0033]根据各矿物的含量、各矿物的权重系数以及常系数，计算得到所述储层样品的表面弛豫率。
[0034]进一步地，所述储层样品的表面弛豫率通过如下公式计算得到：
[0035][0036]其中，ρ为储层样品的表面弛豫率，单位为μm/s；a0为常系数；x
j
为第j种矿物；a
j
为第j种矿物的权重系数；j的取值为1至N，N为确定出的矿物的数量；
[0037]其中，a0为15.3；x1为黄铁矿，a1为
‑
1.72；x2为石英，a2为
‑
2.16；x3为钾长石，a3为0.18；x4为斜长石，a4为1.36；x5为尖晶石，a5为1.93；x6为碳酸盐岩矿物，a6为1.90；x7为黏土，a7为
‑
0.37；以及x8为菱铁矿，a8为1.35。
[0038]优选地，在对所述储层样品进行孔隙提取之前，所述方法还包括：
[0039]对所述储层样品进行抽提以去除其中的可溶有机质；以及
[0040]对抽提处理后的所述储层样品进行氩离子抛光处理。
[0041]第二方面，本文还提供一种储层孔隙孔径分布确定装置，所述装置包括：
[0042]获取模块，用于获取储层样品；
[0043]T2分布谱获取模块，用于对所述储层样品进行核磁共振，得到所述储层样品的T2分布谱；
[0044]综合形状因子计算模块，用于对所述储层样品进行孔隙提取得到多个孔隙，并根据各个孔隙的形状参数计算得到所述储层样品的综合形状因子，所述形状参数至少包括孔隙的最大直径和面积；
[0045]表面弛豫率确定模块，用于对所述储层样品进行X射线衍射分析，并根据X射线衍射分析结果确定所述储层样品的表面弛豫率；
[0046]孔径分布获取模块，用于根据所述综合形状因子和所述表面弛豫率转换所述T2分布谱，得到所述储层样品的孔径分布。
[0047]采用上述技术方案，本文提供的一种储层孔隙孔径分布确定方法和装置，根据识别出的孔隙的最大直径和面积计算得本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种储层孔隙孔径分布确定方法，其特征在于，包括：获取储层样品；对所述储层样品进行核磁共振，得到所述储层样品的T2分布谱；对所述储层样品进行孔隙提取得到多个孔隙，并根据各个孔隙的形状参数计算得到所述储层样品的综合形状因子，所述形状参数至少包括孔隙的最大直径和面积；对所述储层样品进行X射线衍射分析，并根据X射线衍射分析结果确定所述储层样品的表面弛豫率；根据所述综合形状因子和所述表面弛豫率转换所述T2分布谱，得到所述储层样品的孔径分布。2.根据权利要求1所述的方法，其特征在于，对所述储层样品进行孔隙提取得到多个孔隙，并根据各个孔隙的形状参数计算得到所述储层样品的综合形状因子，进一步为：通过扫描电镜获取所述储层样品的图像；根据预设的灰度阈值识别出所述图像中的多个孔隙及其形状参数；根据各个孔隙的形状参数，计算得到各个孔隙的形状因子；根据各个孔隙的形状因子和各个孔隙的面积，计算得到所述储层样品的综合形状因子。3.根据权利要求2所述的方法，其特征在于，所述各个孔隙的形状因子通过如下公式计算得到：其中，F
i
为第i个孔隙的形状因子；D
imax
为第i个孔隙的最大直径；S
i
为第i个孔隙的面积。4.根据权利要求2所述的方法，其特征在于，所述形状参数还包括孔隙的最小直径和周长，所述各个孔隙的形状因子通过如下公式计算得到：其中，F
i
为第i个孔隙的形状因子；D
imax
为第i个孔隙的最大直径；D
imin
为第i个孔隙的最小直径；L
i
为第i个孔隙的周长；S
i
为第i个孔隙的面积。5.根据权利要求2所述的方法，其特征在于，所述储层样品的综合形状因子通过如下公式计算得到：其中，F为综合形状因子；F
i
为第i个孔隙的形状因子；S
i
为第i个孔隙的面积；i的取值范围为1至n，n为识别出的孔隙的数量。6.根据权利要求2所述的方法，其特征在于，在根据各个孔隙的形状因子和各个孔隙的面积，计算得到所述储层样品的综合形状...

【专利技术属性】
技术研发人员：姜福杰，陈迪，张晨曦，霍丽娜，庞雄奇，陈冬霞，张雪，胡涛，宋泽章，朱晨曦，庞宏，陈君青，黄任达，
申请(专利权)人：中国石油大学北京，
类型：发明
国别省市：