发明名称 |
岩心孔隙度参数场的获取方法 |
摘要 |
本发明公开了一种岩心孔隙度参数场的获取方法,属于石油地球物理勘探领域。该方法包括:在k张岩心二值化图像上切割边长为L的方形截面,并将其对应的多个二维矩阵按三维顺序叠加,形成三维矩阵并封装至单元数组的各个元胞数组中,获取每个元胞数组对应的孔隙度,进而获得单元数组的孔隙度参数场。根据单元数组的孔隙度参数场计算三维切割样品的平均孔隙度模拟值,测定岩心样品的平均孔隙度实测值,当平均孔隙度模拟值和平均孔隙度实测值的误差小于5%时,确定L和k的取值为切割参数,并且所获取的单元数组的孔隙度参数场作为岩心样品的岩心孔隙度参数场。该方法能保证获得的岩心孔隙度参数场真实反映岩心样品孔隙度,且减少数据量。 |
申请公布号 |
CN106525680A |
申请公布日期 |
2017.03.22 |
申请号 |
CN201610884265.1 |
申请日期 |
2016.10.10 |
申请人 |
中国石油天然气股份有限公司 |
发明人 |
李玥洋;赵益;徐昌海;杨舒然 |
分类号 |
G01N15/08(2006.01)I |
主分类号 |
G01N15/08(2006.01)I |
代理机构 |
北京三高永信知识产权代理有限责任公司 11138 |
代理人 |
董亚军 |
主权项 |
一种岩心孔隙度参数场的获取方法,包括:利用X射线沿岩心样品的轴向间隔预定距离对所述岩心样品进行CT扫描,得到k张岩心CT图像;对所述k张岩心CT图像依次进行灰度化处理、对比度增强处理、二值化处理,得到k张岩心二值化图像,其特征在于,所述方法还包括:步骤a、在k张所述岩心二值化图像上分别切割边长为L的方形截面,并将k张所述方形截面所对应的多个二维矩阵按三维顺序叠加,形成三维矩阵;步骤b、将所述三维矩阵封装至单元数组的各个元胞数组中,通过对每个所述元胞数组中的孔隙度参数进行局部反演,获得每个所述元胞数组对应的孔隙度,进而获得所述单元数组的孔隙度参数场;步骤c、根据所述单元数组的孔隙度参数场计算三维切割样品的平均孔隙度模拟值,测定所述岩心样品的平均孔隙度实测值,当所述平均孔隙度模拟值和所述平均孔隙度实测值之间的误差小于5%时,确定所述L和所述k的取值为切割参数,并且将基于所述切割参数所获取的单元数组的孔隙度参数场作为所述岩心样品的岩心孔隙度参数场。 |
地址 |
100007 北京市东城区东直门北大街9号中国石油大厦 |