[发明专利]一种岩层孔隙率分布计算方法

权利要求书

1.一种岩层孔隙率分布计算方法，其特征在于包括如下步骤：

S1.对岩层裂隙富集区域进行模拟，输出对应的地层运动原始三彩图像；

S2.对地层运动原始三彩图像进行预处理，生成对应的三维图像数字矩阵；

S3.通过分析三维图像数字矩阵，去除多余岩层间像素信息，得到岩层裂隙层间分布图；具体包括如下步骤：

A1.去除步骤S1输出的地层运动原始三彩图像image1中的多余上网格线；具体包括：以下部网格线沿上像素行，所有列为基准，设：

w₁＝I(i,j,1)/I(i,j,2)

w₂＝I(i,j,2)/I(i,j,3)

其中，w₁表示像素点红色通道对应的像素值与绿色通道像素值的比值；w₂表示像素点绿色通道对应的像素值与蓝色通道像素值的比值；I(i,j,1)表示像素点在红通道的像素值；I(i,j,2)表示像素点在绿通道的像素值；I(i,j,3)表示像素点在蓝通道的像素值；

若满足网格线沿上设定范围，则令I(i,j,1)＝I(i,j,2)＝I(i,j,3)＝0；网格线沿上设定范围为0.9＜w₁,w₂＜1.1；

A2.去除步骤S1输出的地层运动原始三彩图像image1中的多余下网格线；具体包括，以下部网络线沿下像素行，所有列为基准，设：

sum＝I(i,j,1)+I(i,j,2)+I(i,j,3)

满足网格线向下设定范围，则令I(i,j,1)＝I(i,j,2)＝I(i,j,3)＝0；I(i,j,1)表示像素点在红通道的像素值；I(i,j,2)表示像素点在绿通道的像素值；I(i,j,3)表示像素点在蓝通道的像素值；sum表示像素点在红、绿、蓝三通道的像素值之和；网格线向下设定范围为130＜sum＜170；

S4.沿岩层方向遍历岩层裂隙层间分布图，得到岩层裂隙代表像素点和岩层裂隙信息提取区域图；并绘制岩层裂隙层孔隙率对应的分布曲线图；

S5.设置全图像范围遍历窗口，对岩层裂隙层间分布图进行全图像范围遍历，全图像范围遍历窗口沿竖直和水平走向，生成岩层孔隙率分布云图。

2.根据权利要求1所述的岩层孔隙率分布计算方法，其特征在于所述的步骤S1具体包括：利用基于离散单元法，对岩层裂隙富集区域的地层进行模拟，对岩石基质和岩石接触进行参数化处理，生成岩层条件参数，根据岩层条件参数，输出对应的地层运动原始三彩图像image1；原彩图存储了以三种颜色[r,g,b]通道为颜色模式的像素值。

3.根据权利要求2所述的岩层孔隙率分布计算方法，其特征在于所述的岩层条件参数包括：块体的质量密度、弹性性能、塑性性能用于确定材料的力学特性；刚度、剪切刚度、抗拉强度、摩擦角、内聚力用于确定块体间接触的力学特性；地层块体尺寸。

4.根据权利要求3所述的岩层孔隙率分布计算方法，其特征在于所述的步骤S2中，预处理包括：将地层运动原始三彩图像image1生成对应的三维图像数字矩阵Array1，数字矩阵Array1采用三维形式反映三种颜色[r,g,b]通道，并分别在其对应通道以二维数组形式[[·]i]j显示像素位置，并通过像素值反映三种原彩色素值；其中，每一个像素位置包含三个颜色通道，红绿蓝三色像素值以[r,g,b]的形式，在像素位置中分布，同时r,g,b∈[0,255]。