[发明专利]一种测井渗透率解释方法

说明书

技术领域

本发明涉及一种测井渗透率解释方法，属于采矿技术领域。

背景技术

渗透率是储层评价中非常重要的参数。目前，获得储层渗透率的途径主要有两种，一种是通过实验或工程手段直接测得，另一种是基于地球物理方法解释求取。前一种包括岩心实验测量法、电缆地层测试法、钻杆地层测试法等，后一种包括测井解释渗透率和地震解释渗透率等，由于获取手段不同，加之渗透率本身敏感性，其数值之间相差很大。

而校验渗透率是否精准主要有两个方面，一是与岩心实测的对比，但是由于样品质量、渗透率测量的精度影响，实测数据与地层真实渗透率有一定偏差；二是与生产动态相结合，除了人为工作制度的改变，储层质量是影响生产的根本原因。

渗透率作为储层物性研究的一个关键参数，其解释模型一直较难建立。利用测井资料进行渗透率解释是目前最常用做法，常规测井软件解释采用经验公式和理论模型来计算渗透率，由于公式或模型的经验性、人为因素影响、地区差异及缺乏生产动态验证，常导致最终的解释模型不能广泛应用或与实际结果不符。

发明内容

为了解决现有技术的不足，本发明提供了一种测井渗透率解释方法，该方法基于岩心、实验、测井、测试等资料，采用数学统计法及动静态相互校验法，建立渗透率多元线性回归测井解释模型，可以正确有效地计算渗透率。

本发明为解决其技术问题所采用的技术方案是：提供了一种测井渗透率解释方法，具体包括以下步骤：

(1)粒度解释：采用以下公式进行粒度解释：

M_d＝-0.346GR+18.775ρ+2.497γ+0.002GR²

-17.723ρ²+4.089ρ³-0.095γ²+0.001γ³

其中，M_d表示粒度，GR表示伽马测井，单位为API；ρ表示测井密度，单位为g/cm³；γ表示中子测井，单位为v/v；所述的伽马测井、测井密度和中子测井均为已知的电测井数值常数；

(2)孔隙度解释：对中子密度曲线进行质量控制和附加校正后，采用中子-密度交会方法得到地层孔隙度；

(3)渗透率解释：采用以下公式渗透率解释模型进行渗透率解释：

其中，K表示渗透率，单位为mD；Md表示步骤(1)计算得到的粒度，单位为Φ；表示步骤(2)计算得到的孔隙度，单位为f；P表示测井视电阻率，单位为Ω·m，所述的测井视电阻率为测井数据中的已知量。

步骤(3)进行渗透率解释后，通过以下步骤对渗透率解释模型进行校验：

(4-1)根据步骤(1)得到的粒度确定其中值，再根据实际采集到的岩样的粒度确定岩样粒度中值，对比粒度中值和岩样粒度中值，做其相关性散点图，若散点分布于45°线的两侧越集中，则粒度中值和岩样粒度中值的相关性越好；

(4-2)渗透率模型实测检验：根据步骤(3)的解释的渗透率作出测井解释渗透率散点图，将其与岩心样品实测渗透率散点图进行对比，若散点分布于45°线两侧越集中，则相关性越好。

(4-3)渗透率模型实测生产特征检验：对于单井，根据地层系数与生产测井的产量剖面对比检验渗透率解释模型，地层系数即地层有效厚度与有效渗透率的乘积，用于判定注水井各层位注入量，其中，有效厚度表示在现有经济技术条件下可以开采的油层厚度，油层厚度不包含差油层、油水同层、含油水层以及含水油层，物性越好，渗透率值越高，其地层系数越大，生产井产能好，注水井注入量高。

本发明基于其技术方案所具有的有益效果在于：该方法结合静态实验、测井及动态测试资料，利用数学统计的科学方法，经过基于地质条件的客观分析评价后对解释结果加以校验，模型结果客观真实且便于使用。

附图说明

图1是喉道分维数与孔隙度、渗透率、门槛压力、最大孔喉半径相关图，其中图1(1)是喉道分维数与孔隙度相关图，图1(2)是喉道分维数与渗透率相关图，图1(3)是喉道分维数与门槛压力相关图，图1(4)是喉道分维数与最大孔喉半径相关图。

图2是孔隙分维数与孔隙度、渗透率、门槛压力、最大半径相关图。其中图2(1)是孔隙分维数与孔隙度相关图，图2(2)是孔隙分维数与渗透率相关图，图2(3)是孔隙分维数与门槛压力相关图，图2(4)是孔隙分维数与最大孔喉半径相关图。

图3是岩样孔隙度分布直方图。