[发明专利]一种低渗砂砾岩储层数字岩心结构提取方法

权利要求书

1.一种低渗砂砾岩储层数字岩心结构提取方法，其特征在于，包括以下步骤：

步骤1：获取目标储层的岩心铸体薄片图像，结合比例尺统计矿物组分的粒径，基于矿物组分粒径大小及颗粒间的孔隙形态，确定储层岩心结构类型的个数；

步骤2：获取目标储层的数字岩心，对所述的数字岩心进行切片，绘制对应切片的数字岩心体素的灰度直方图，根据灰度直方图统计不同岩心结构类型的数字岩心体素灰度值范围；

步骤3：基于储层岩心结构类型的个数，确定数字岩心结构分析所需阈值的个数；基于不同岩心结构类型体素灰度值范围，确定不同岩心结构类型体素灰度值大小关系；

步骤4：利用最大类间距方法，对数字岩心切片进行多阈值优选计算，得到砂砾岩储层数字岩心最优阈值；

步骤5：基于最优阈值结果，进行数字岩心阈值分割，获取不同岩心结构类型的数字岩心体素点阵，完成数字岩心结构提取。

2.根据权利要求1所述的一种低渗砂砾岩储层数字岩心结构提取方法，其特征在于，步骤1所述的获取目标储层的岩心铸体薄片图像，具体是：从低渗砂砾岩储层的岩心获取岩心样品并制备铸体薄片，再通过图像收集系统获取铸体薄片不同位置的标准图像，所述的标准图像，是能够通过颜色差异区分矿物颗粒、孔隙裂缝信息。

3.根据权利要求1所述的一种低渗砂砾岩储层数字岩心结构提取方法，其特征在于，步骤1所述的结合比例尺统计矿物组分粒径，对所述铸体薄片图像中矿物组分颗粒绘制最大外接圆，测量全部外接圆直径数值，得到矿物组分粒径。

4.根据权利要求1所述的一种低渗砂砾岩储层数字岩心结构提取方法，其特征在于，步骤1所述的储层岩心结构类型，包括矿物颗粒、岩屑、填隙物、孔隙、裂缝。

5.根据权利要求1所述的一种低渗砂砾岩储层数字岩心结构提取方法，其特征在于，步骤2所述的获取目标储层的数字岩心，是对获取的岩心样品进行微米CT扫描，获得三维灰度值图像，构成数字岩心，所述的数字岩心能够通过灰度值差异表征三维岩心骨架分布情况。

6.根据权利要求1所述的一种低渗砂砾岩储层数字岩心结构提取方法，其特征在于，步骤2所述的对数字岩心进行切片，是对数字岩心任一方向进行切割后，取能够表征岩心结构特征截面图像。

7.根据权利要求1所述的一种低渗砂砾岩储层数字岩心结构提取方法，其特征在于，步骤2所述的根据灰度直方图统计不同岩心结构类型的数字岩心体素灰度值范围，是对所述数字岩心切片进行灰度值统计及分类，获得不同岩心结构类型的数字岩心体素灰度值范围，确定灰度直方图中不同岩心结构类型的数字岩心体素分布。

8.根据权利要求1所述的一种低渗砂砾岩储层数字岩心结构提取方法，其特征在于，步骤3所述的确定数字岩心结构分析所需阈值个数，是储层岩心结构类型的个数为n时对应阈值的个数为n-1。

9.根据权利要求1所述的一种低渗砂砾岩储层数字岩心结构提取方法，其特征在于，步骤4包括：根据数字岩心结构分析所需阈值个数，将全部数字岩心体素划分为n组集合，针对各集合数字岩心体素灰度值，计算组间方差和组内方差，当组间方差和组内方差两者之比达到最大值时，对应的阈值即为最优阈值；由于数字岩心结构分析所需阈值个数为n-1，因此需要结合不同岩心结构类型的数字岩心体素灰度值大小进行嵌套循环、迭代计算，获取最优阈值。

10.根据权利要求1所述的一种低渗砂砾岩储层数字岩心结构提取方法，其特征在于，步骤5包括：根据步骤4所获个数为n-1的最优阈值结果，将数字岩心体素按照灰度值大小分为n组体素点阵，根据步骤3所获的不同岩心结构类型体素灰度值大小关系确定各组体素点阵的岩心结构类型，实现数字岩心结构提取。