一种井壁渗透率计算方法及装置制造方法及图纸

本申请实施例提供了一种井壁渗透率计算方法及装置。该方法包括：获取目标地层的电成像测井图像及测井数据；根据所述测井数据，计算所述电成像测井图像中每个像素点的孔隙度；根据所述电成像测井图像中每个像素点的孔隙度，计算所述电成像测井图像中每个像素点的渗透率；根据所述电成像测井图像中每个像素点的渗透率，获得所述目标地层全井段的地层渗透率值。本申请实施例的方法基于电成像测井图像，同时考虑了井壁中不同位置孔隙类型对渗透率的影响，方法简单实用，可以准确的获得碳酸盐地层井壁渗透率，为有效储层评价提供了基础。

【技术实现步骤摘要】

本专利技术涉及油田勘探测井
，尤其涉及一种井壁渗透率计算方法及装置。
技术介绍
渗透率参数是油田勘探开发过程中极其重要的参数之一，其计算精度的高低对油田开发效益具有重要影响。由于碳酸盐岩储层极强非均质性和各向异性的存在，碳酸盐岩渗透率的确定一直是测井评价的难题之一。目前测井计算渗透率主要有以下两种方法：一是通过岩心实验资料建立渗透率与孔隙度、粒度中值的拟合经验公式或多元回归等数学手段确定。该方法对于相对均质的碎屑岩储层来说，岩心分析的孔隙度和渗透率通常具有较好的相关性，统计得到的孔渗关系可以求取储层渗透率，但对于极强非均质性的碳酸盐岩储层，通过岩心实验建立的渗透率响应方程并不能代表整个储层的渗流能力，其计算精度和适用性差；二是利用双孔隙介质模型，将碳酸盐岩储层渗透率视为基块孔隙渗透率和裂缝渗透率共同组成，但由于裂缝渗透率一般比基块孔隙渗透率要大得多，故储层中有无裂缝，对其平均渗透率影响极大，因此该方法在油田现场实际应用中难以实现，实用性差。
技术实现思路
本申请实施例提供一种井壁渗透率计算方法及装置，以准确的获得碳酸盐岩地层的渗透率，为有效储层评价提供基础。为达到上述目的，一方面，本申请实施例提供了一种井壁渗透率计算方法，所述方法包括：获取目标地层的电成像测井图像及测井数据；根据所述测井数据，计算所述电成像测井图像中每个像素点的孔隙度；根据所述电成像测井图像中每个像素点的孔隙度，计算所述电成像测井图像中每个像素点的渗透率；根据所述电成像测井图像中每个像素点的渗透率，获得所述目标地层全井段的地层渗透率值。进一步地，在根据所述电成像测井图像中每个像素点的孔隙度，计算所述电成像测井图像中每个像素点的渗透率之前，还包括：利用图像分析技术，识别所述电成像测井图像中各图像特征单元对应的地质特征类型，并据此获得所述电成像测井图像中每个像素点的地质特征类型；根据所述电成像测井图像中每个像素点的地质特征类型，获取所述电成像测井图像中每个像素点的渗透性有效系数；对应的，根据所述电成像测井图像中像素点的孔隙度以及渗透性有效系数，计算所述电成像测井图像中每个像素点的渗透率。进一步地，采用以下公式计算所述电成像测井图像中每个像素点的渗透率：PERMi=kei*10ka*φi+kb]]>式中，PERMi表示电成像测井图像中像素点的渗透率；φi表示像素点的孔隙度；ka是孔隙度指数；kb是渗透率常数；kei是像素点所属的地质特征的渗透性有效系数。进一步地，所述利用图像分析技术，识别所述电成像测井图像中各图像特征单元对应的地质特征类型，包括：利用图像分割方法将所述电成像测井图像分割为三种类型的图像特征单元，所述三种类型的图像特征单元为高阻图像特征单元、过渡图像特征单元以及高导图像特征单元；定量计算所述图像特征单元的特征参数；将所述图像特征单元的特征参数与预设的地质特征类型的特征参数进行对比，获得所述图像特征单元对应的地质特征类型。进一步地，所述特征参数包括：特征灰度平均值、面积、周长、圆度、宽长比、粒径、主方向、对比度以及均匀性。进一步地，所述地质特征类型包括：低阻裂缝、溶孔、洞穴、泥质条带、泥质纹层、泥质厚层、诱导缝、缝合线、角砾以及高阻缝。进一步地，根据所述电成像测井图像中每个像素点的渗透率，获得所述目标地层全井段的地层渗透率值，包括：根据所述电成像测井图像中每个像素点的渗透率，获得井壁渗透率图像；对所述井壁渗透率图像进行逐深度窗统计，计算每个深度窗的渗透率，获得所述目标地层全井段的地层渗透率值。进一步地，采用以下公式计算每个深度窗的渗透率：PERM=Σi=0npixelPERMi]]>式中，PERM是对应深度窗的渗透率；npixel是深度窗内像素点数；PERMi像素点的渗透率。进一步地，采用阿尔奇公式计算所述电成像测井图像中每个像素点的孔隙度：φi=a·b·Rmf·Cim]]>其中：φi是计算的电成像测井图像中像素点的孔隙度；a、b是阿尔奇公式中岩性系数；m是地层胶结指数；Rmf是泥浆滤液电阻率；Ci是电成像测井图像中像素点的电导率值。另一方面，本申请实施例还提供了一种井壁渗透率计算装置，所述装置包括：获取单元，用于获取目标地层的电成像测井图像及测井数据；像素点孔隙度计算单元，用于根据所述测井数据，计算所述电成像测井图像中每个像素点的孔隙度；像素点渗透率计算单元，用于根据所述电成像测井图像中每个像素点的孔隙度，计算所述电成像测井图像中每个像素点的渗透率；地层渗透率获取单元，用于根据所述电成像测井图像中每个像素点的渗透率，获得所述目标地层全井段的地层渗透率值。进一步地，所述装置还包括：像素点地质特征类型获取单元，用于利用图像分析技术，识别所述电成像测井图像中各图像特征单元对应的地质特征类型，并据此获得所述电成像测井图像中每个像素点的地质特征类型；像素点渗透性有效系数获取单元，用于根据所述电成像测井图像中每个像素点的地质特征类型，获取所述电成像测井图像中每个像素点的渗透性有效系数；对应的，像素点渗透率计算单元用于根据所述电成像测井图像中每个像素点的孔隙度以及渗透性有效系数，计算所述电成像测井图像中每个像素点的渗透率。进一步地，所述像素点渗透率计算单元采用以下公式计算所述电成像测井图像中每个像素点的渗透率：PERMi=kei*10ka*φi+kb]]>式中，PERMi表示电成像测井图像中像素点的渗透率；φi表示像素点的孔隙度；ka是孔隙度指数；kb是渗透率常数；kei是像素点所属的地质特征的渗透性有效系数。进一步地，所述像素点地质特征类型获取单元包括：图像分割子单元，用于利用图像分割方法将所述电成像测井图像分割为三种类型的图像特征单元，所述三种类型的图像特征单元为高阻图像特征单元、过渡图像特征单元以及高导图像特征单元；特征参数计算子单元，用于定量计算所述图像特征单元的特征参数；地质特征类型获取子单元，用于将所述图像特征单元的特征参数与预设的地质特征类型的特征参数进行对比，获得所述图像特征单元对应的地质特征类型。进一步地，所述特征参数包括：特征灰度平均值、面积、周长、圆度、宽长比、粒径、主方向、对比度以及均匀性。进一步地，所述地质特征类型包括：低阻裂缝、溶孔、洞穴、泥质条带、泥质纹层、泥质厚层、诱导缝、缝合线、角砾以及高阻缝。进一步地，所述地层渗透率获取单元包括：井壁渗透率图像获取子单元，用于根据所述电成像测井图像中每个像素点的渗透率，获得井壁渗透率图像；深度窗渗透率计算子单元，用于对所述井壁渗透率图像进行逐深度窗统计，计算每个深度窗的渗透率，获得所述目标地层全井段的地层渗透率值。进一步地，所述深度窗渗透率计算子单元采用以下公式计算每个深度窗的渗透率：PERM=Σi=0npixelPERMi]]>式中，PERM是对应深度窗的渗透率；npixel是深度窗内像素点数；PERMi像素点的渗透率。进一步地，所述像素点孔隙度计算单元采用阿尔奇公式计算所述电成像测井图像中每个像素点的孔隙度：φi=a·b·Rmf·Cim]]>其中：φi是计算的电成像测井图像中像素本文档来自技高网...

【技术保护点】
一种井壁渗透率计算方法，其特征在于，所述方法包括：获取目标地层的电成像测井图像及测井数据；根据所述测井数据，计算所述电成像测井图像中每个像素点的孔隙度；根据所述电成像测井图像中每个像素点的孔隙度，计算所述电成像测井图像中每个像素点的渗透率；根据所述电成像测井图像中每个像素点的渗透率，获得所述目标地层全井段的地层渗透率值。

【技术特征摘要】
1.一种井壁渗透率计算方法，其特征在于，所述方法包括：获取目标地层的电成像测井图像及测井数据；根据所述测井数据，计算所述电成像测井图像中每个像素点的孔隙度；根据所述电成像测井图像中每个像素点的孔隙度，计算所述电成像测井图像中每个像素点的渗透率；根据所述电成像测井图像中每个像素点的渗透率，获得所述目标地层全井段的地层渗透率值。2.如权利要求1所述的井壁渗透率计算方法，其特征在于，在根据所述电成像测井图像中每个像素点的孔隙度，计算所述电成像测井图像中每个像素点的渗透率之前，还包括：利用图像分析技术，识别所述电成像测井图像中各图像特征单元对应的地质特征类型，并据此获得所述电成像测井图像中每个像素点的地质特征类型；根据所述电成像测井图像中每个像素点的地质特征类型，获取所述电成像测井图像中每个像素点的渗透性有效系数；对应的，根据所述电成像测井图像中每个像素点的孔隙度，计算所述电成像测井图像中每个像素点的渗透率，具体为：根据所述电成像测井图像中每个像素点的孔隙度以及渗透性有效系数，计算所述电成像测井图像中每个像素点的渗透率。3.如权利要求2所述的井壁渗透率计算方法，其特征在于，采用以下公式计算所述电成像测井图像中每个像素点的渗透率：PERMi=kei*10ka*φi+kb]]>式中，PERMi表示电成像测井图像中像素点的渗透率；φi表示像素点的孔隙度；ka是孔隙度指数；kb是渗透率常数；kei是像素点所属的地质特征的渗透性有效系数。4.如权利要求2所述的井壁渗透率计算方法，其特征在于，所述利用图像分析技术，识别所述电成像测井图像中各图像特征单元对应的地质特征类型，包括：利用图像分割方法将所述电成像测井图像分割为三种类型的图像特征单元，所述三种类型的图像特征单元为高阻图像特征单元、过渡图像特征单元以及高导图像特征单元；定量计算所述图像特征单元的特征参数；将所述图像特征单元的特征参数与预设的地质特征类型的特征参数进行对比，获得所述图像特征单元对应的地质特征类型。5.如权利要求4所述的井壁渗透率计算方法，其特征在于，所述特征参数包括：特征灰度平均值、面积、周长、圆度、宽长比、粒径、主方向、对比度以及均匀性。6.如权利要求4所述的井壁渗透率计算方法，其特征在于，所述地质特征类型包括：低阻裂缝、溶孔、洞穴、泥质条带、泥质纹层、泥质厚层、诱导缝、缝合线、角砾以及高阻缝。7.如权利要求1所述的井壁渗透率计算方法，其特征在于，根据所述电成像测井图像中每个像素点的渗透率，获得所述目标地层全井段的地层渗透率值，包括：根据所述电成像测井图像中每个像素点的渗透率，获得井壁渗透率图像；对所述井壁渗透率图像进行逐深度窗统计，计算每个深度窗的渗透率，获得所述目标地层全井段的地层渗透率值。8.如权利要求7所述的井壁渗透率计算方法，其特征在于，采用以下公式计算每个深度窗的渗透率：PERM=Σi=0npixelPERMi]]>式中，PERM是对应深度窗的渗透率；npixel是深度窗内像素点数；PERMi像素点的渗透率。9.如权利要求1所述的井壁渗透率计算方法，其特征在于，采用阿尔奇公式计算所述电成像测井图像中每个像素点的孔隙度：φi=a·b·Rmf·Cim]]>其中：φi是计算的电成像测井图像中像素点的孔隙度；a、b是阿尔奇公式中岩性系数；m是地层胶结指数；Rmf是泥浆滤液电阻率；Ci是电成像测井图像中像素点的电导率值。10.一种井壁渗透率计算装置，其特征在于，所述装置包括：获取...

【专利技术属性】
技术研发人员：李宁，冯周，武宏亮，王克文，冯庆付，张宫，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京;11