灰度反转图像。同时将Well矩阵中的数值也赋值为GRAY，Well矩阵将用于存储类型识别结果。
[0039]步骤六:以岩心照片每个象素的灰度数值高低作为砂岩、泥披和裂缝的判读依据，然后绘制三维灰度图。泥披的灰度范围一般为180-256 ;砂岩的灰度范围一般为50-180 ;裂缝的灰度范围一般为1-50。然后根据灰度数据，确定泥披、裂缝和砂岩。
[0040]步骤七:根据步骤六中裂缝的灰度范围(0-50)，确定裂缝区域的图像像素坐标。然后将Well矩阵中对应的数值设定为-9999，表示无效类型。
[0041]步骤八:根据步骤六中选取的砂岩和泥披的灰度数值范围，将Well矩阵中泥披区域的类型数值设定为1，砂岩区域和裂缝区域的类型数值设定为O。
[0042]步骤九:根据泥披识别二色图，计算泥质含量(ISH)曲线。
[0043]作为优选，所述步骤九中计算泥质含量曲线具体步骤为:
[0044]SI统计泥披二色图的像素行数目Row和列数目Col。
[0045]S2统计Well类型矩阵中每一列中等于I的泥披像素的数目，计为SumJ。
[0046]S3根据S2得的泥披像素点和，按照ISH = 100*Sum_j/ (Row)，得到每一列的泥质含量ISH。
[0047]S4依次计算每一列位置点的泥质含量曲线，输出总体泥质含量曲线。
[0048]实施例一:
[0049]根据本发明提出的基于岩心扫描的富泥披储层泥质含量计算方法，对阿尔伯达盆地麦凯河区块海相储层中的泥质含量进行了计算。具体实施步骤如下:
[0050]如图1所示，⑴取心井为Song-1井，取心段长度为1.5m。本次实施首先将取心段以500mm为间隔切割取心段，将取心段分割为3段。然后将3段岩心分别放在独立的510mm长度的3个岩心盒中。对3盒岩心进行切割，然后露出新鲜剖面，进行扫描拍照。每一盒岩心扫描照片的扫描精度为1_。扫描图片的尺寸为:宽度100个像素点，即100_ ;长度510个像素点，即510mm。
[0051]如图2所示，⑵将每一口井所有N段岩心照片全部拼接在一起。
[0052]如图3所示，(3)将岩心扫描图像Well转换为灰度图，灰度存储矩阵为GRAY。使用GRAY* (-1)+255计算公式对灰度矩阵GRAY进行处理，得到灰度反转图像。同时将Well矩阵中的数值也赋值为GRAY，Well矩阵将用于存储类型识别结果。
[0053]如图4所示，(4)以岩心照片每个象素的灰度数值高低作为砂岩、泥披和裂缝的判读依据，然后绘制三维灰度图。泥披的灰度范围为180-256，砂岩的灰度范围为60-180，裂缝的灰度范围为1-60。
[0054]如图5所示，(5)将Well矩阵中泥披区域的类型数值设定为1，砂岩区域和裂缝区域的类型数值设定为0，绘制泥披识别二色图。
[0055]如图6所示，(6)根据泥披识别二色图，计算泥质含量曲线ISH。
[0056]本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的实施方法，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。
【主权项】
1.一种基于岩心扫描的富泥披储层泥质含量计算方法，包括以下步骤: 步骤一:取心岩石样品，采用低温二氧化碳冷冻保存，转运到岩心切割实验室，每口井的长度为lm-20m的岩心样品，按深度顺序从上到下顺序，以500mm为间隔切割为N段，分别放置在510mm长的岩心盒中；每一个岩心盒标的侧面标上序号；然后在岩心盒上标明顶底深度，顶部深度为TOP，底部深度为B0TT0W，将每一个盒岩心放置到岩心切割仪上； 步骤二:将切割仪刀具移动到取心岩石样品中间，下刀切割；然后将切割开的岩石样品的新鲜面露出，采用扫描仪从顶到底扫描，扫描图像存储到二维矩阵ImageUndex}中； 步骤三:按照步骤一和步骤二，对每一盒岩心进行切割扫描，按照深度顺序，从顶到底，将顶底吻合的两张照片拼接在一起，使相邻的一张照片的底部与下一张照片的顶部吻合；依次遍历所有岩心照片，将每一口井所有N段岩心照片全部拼接在一起，然后将拼接后的图像存储到二维矩阵Well中； 步骤四:将岩心扫描图像Well转换为灰度图，灰度存储矩阵为GRAY ; 步骤五:使用GRAY*(-l)+255计算公式对步骤四得到的灰度矩阵GRAY进行处理，得到灰度反转图像；同时将WelI矩阵中的数值也赋值为GRAY，Well矩阵将用于存储类型识别结果; 步骤六:以岩心照片每个象素的灰度数值高低作为砂岩、泥披和裂缝的判读依据，绘制三维灰度图；泥披的灰度范围为180-256 ;砂岩的灰度范围一般为50-180 ;裂缝的灰度范围为1-50 ;根据灰度数据，确定泥披、裂缝和砂岩； 步骤七:根据步骤六中裂缝的灰度范围，确定裂缝区域的图像像素坐标，将Well矩阵中对应的数值设定为-9999，表示无效类型； 步骤八:根据步骤六中选取的砂岩和泥披的灰度数值范围，将Well矩阵中泥披区域的类型数值设定为1，砂岩区域和裂缝区域的类型数值设定为O ; 步骤九:根据泥披识别二色图，计算泥质含量曲线。2.根据权利要求1所述的一种基于岩心扫描的富泥披储层泥质含量计算方法，其特征在于:所述步骤二中每一盒岩心扫描照片的扫描精度为1mm，扫描图片的尺寸为:宽度100个像素点，即10mm ;长度510个像素点，即510mm。3.根据权利要求2所述的一种基于岩心扫描的富泥披储层泥质含量计算方法，其特征在于:所述步骤三中Well中存储的图片的尺寸为:宽度为100个像素点，长度为510*N个像素点。4.根据权利要求3所述的一种基于岩心扫描的富泥披储层泥质含量计算方法，其特征在于:所述步骤四中RAY矩阵中黑颜色的灰度数值为255，白颜色的灰度数值为O。5.根据权利要求4所述的一种基于岩心扫描的富泥披储层泥质含量计算方法，其特征在于:所述步骤九中计算泥质含量曲线具体步骤为: SI统计泥披二色图的像素行数目Row和列数目Col ; S2统计Well类型矩阵中每一列中等于I的泥披像素的数目，计为Sum_j ; S3根据S2得的泥披像素点和，按照ISH = 100*Sum_j/(ROW)，得到每一列的泥质含量ISH0
【专利摘要】本发明公开了一种基于岩心扫描的富泥披储层泥质含量计算方法，包括以下步骤：第一、取心岩石样品；第二、下刀切割；第三、切割扫描；第四、岩心扫描图像转换为灰度图；第五、得到灰度反转图像；第六、绘制三维灰度图；第七、确定裂缝区域的图像像素坐标；第八、将Well矩阵中泥披区域的类型数值设定为1，砂岩区域和裂缝区域的类型数值设定为0；第九、计算泥质含量曲线。本发明的有益效果如下：针对海相储层的岩心扫描结果，基于富泥披储层泥质含量计算方法，首次计算出了分辨率达到mm级的泥质含量曲线。解决常规自然伽马测井曲线对mm-cm基泥披层无响应，无法精确计算泥质含量的问题。
【IPC分类】G01N23/04
【公开号】CN105181721
【申请号】CN201510646484
【发明人】唐明明, 黄继新, 卢双舫, 陈和平, 梁宏儒, 燕贲惠, 马慧芳, 薛松, 刘爽
【申请人】中国石油大学(华东), 中国石油天然气股份有限公司勘探开发研究院
【公开日】2015年12月23日
【申请日】2015年10月8日