[发明专利]岩心孔隙度参数场的获取方法有效
申请号: | 201610884265.1 | 申请日: | 2016-10-10 |
公开(公告)号: | CN106525680B | 公开(公告)日: | 2020-01-07 |
发明(设计)人: | 李玥洋;赵益;徐昌海;杨舒然 | 申请(专利权)人: | 中国石油天然气股份有限公司 |
主分类号: | G01N15/08 | 分类号: | G01N15/08 |
代理公司: | 11138 北京三高永信知识产权代理有限责任公司 | 代理人: | 董亚军 |
地址: | 100007 北京市*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明公开了一种岩心孔隙度参数场的获取方法,属于石油地球物理勘探领域。该方法包括:在k张岩心二值化图像上切割边长为L的方形截面,并将其对应的多个二维矩阵按三维顺序叠加,形成三维矩阵并封装至单元数组的各个元胞数组中,获取每个元胞数组对应的孔隙度,进而获得单元数组的孔隙度参数场。根据单元数组的孔隙度参数场计算三维切割样品的平均孔隙度模拟值,测定岩心样品的平均孔隙度实测值,当平均孔隙度模拟值和平均孔隙度实测值的误差小于5%时,确定L和k的取值为切割参数,并且所获取的单元数组的孔隙度参数场作为岩心样品的岩心孔隙度参数场。该方法能保证获得的岩心孔隙度参数场真实反映岩心样品孔隙度,且减少数据量。 | ||
搜索关键词: | 参数场 孔隙度 平均孔隙 岩心孔隙度 单元数组 岩心样品 数组 元胞 实测 切割 三维 石油地球物理勘探 二值化图像 二维矩阵 方形截面 获得单元 切割参数 三维矩阵 顺序叠加 岩心 数据量 边长 封装 保证 | ||
【主权项】:
1.一种岩心孔隙度参数场的获取方法,包括:利用X射线沿岩心样品的轴向间隔预定距离对所述岩心样品进行CT扫描,得到k张岩心CT图像;对所述k张岩心CT图像依次进行灰度化处理、对比度增强处理、二值化处理,得到k张岩心二值化图像,其特征在于,所述方法还包括:/n步骤a、在k张所述岩心二值化图像上分别切割边长为L的方形截面,并将k张所述方形截面所对应的多个二维矩阵按三维顺序叠加,形成三维矩阵;/n步骤b、将所述三维矩阵封装至单元数组的各个元胞数组中,通过对每个所述元胞数组中的孔隙度参数进行局部反演,获得每个所述元胞数组对应的孔隙度,进而获得所述单元数组的孔隙度参数场;/n步骤c、根据所述单元数组的孔隙度参数场计算三维切割样品的平均孔隙度模拟值,测定所述岩心样品的平均孔隙度实测值,当所述平均孔隙度模拟值和所述平均孔隙度实测值之间的误差小于5%时,确定所述L和所述k的取值为切割参数,并且将基于所述切割参数所获取的单元数组的孔隙度参数场作为所述岩心样品的岩心孔隙度参数场,/n其中,所述根据所述单元数组的孔隙度参数场计算三维切割样品的平均孔隙度模拟值包括:/n对每个所述元胞数组中的孔隙度参数进行局部反演,通过第一公式计算所述每个所述元胞数组对应的孔隙度;/n根据每个所述元胞数组对应的孔隙度,通过第二公式计算所述单元数组的平均孔隙度;/n根据所述单元数组的平均孔隙度,通过第三公式计算所述三维切割样品的平均孔隙度模拟值;/n所述第一公式如下所示:/n
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油天然气股份有限公司,未经中国石油天然气股份有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201610884265.1/,转载请声明来源钻瓜专利网。
- 上一篇:智能碳刷旋头打弧机
- 下一篇:致密砂岩稳态法气水相对渗透率曲线测定方法