本发明涉及地质学和油气田勘探，尤其涉及一种基于岩石铸体薄片照片计算岩石孔隙度的方法。
背景技术：
1、岩石铸体薄片资料在地质研究及油气勘探中充当重要作用，岩石铸体薄片孔隙度的求取是确定岩石物性参数及储集性能评价的重要基础。
2、传统对岩石铸体薄片的孔隙度求取主要通过人为经验的目估法，岩石通过铸体磨制薄片后，岩石中孔隙呈现蓝色，通过人为的铸体薄片鉴定，根据个人经验目估在显微镜下的视域中蓝色区域占比多少，大致的确定岩石铸体薄片的孔隙度大小。该类方法存在较大的人为误差，与铸体薄片鉴定人员的经验密切相关，不同人员针对铸体薄片鉴定的孔隙度存在差异。尽管，岩石铸体薄片中孔隙呈现蓝色，在显微镜的单偏光下与岩石矿物颜色差别较大，极易区分矿物组分与孔隙，但孔隙度(孔隙面积占比)的估算仍然存在人为误差和经验误差。另外，在储集空间及储集物性研究中，岩石铸体薄片孔隙度的确定十分关键，迫切需要寻求新的方法准确定量计算孔隙度。
3、图像识别技术能将图片网格化以及确定图片相对坐标、也能确定单位网格的颜色rgb，并根据颜色变化赋予数值等，所以基于图像识别技术在岩石铸体薄片孔隙的识别与统计面积大小具有明显优势。综上所述，为获得精度高、准确的岩石铸体薄片孔隙度大小，开展基于岩石铸体薄片照片计算岩石孔隙度的研究是十分有必要的。
技术实现思路
1、本发明的目的在于提供一种基于岩石铸体薄片照片计算岩石孔隙度的方法，旨在解决现有技术中的人为定性确定孔隙度已满足不了现在研究的需求的精度及准确性的技术问题。
2、为实现上述目的，本发明采用的一种基于岩石铸体薄片照片计算岩石孔隙度的方法，包括如下步骤：
3、获取岩石铸体薄片照片，并对图像进行初步处理；
4、确定岩石铸体薄片照片的相对坐标后，根据岩石铸体薄片照片像素大小及数量，对岩石铸体薄片照片进行网格化处理；
5、确定岩石铸体薄片照片rgb分布数据，并将rgb数据转换成xyz数据；
6、基于直方图阈值分割原理，在xyz数据中，根据孔隙空间范围内z值的阈值，开展孔隙分布区域的拾取，计算孔隙分布区域的面积；
7、根据岩石铸体薄片照片相对面积计算岩石铸体薄片孔隙度；
8、获取岩石铸体薄片的孔隙度。
9、其中，在获取岩石铸体薄片照片，并对图像进行初步处理的步骤中：
10、在单偏光下对岩石铸体薄片进行拍摄，获取岩石铸体薄片照片；
11、将岩石铸体薄片照片输入计算机，进行图像初步处理。
12、其中，在将岩石铸体薄片照片输入计算机，进行图像初步处理的步骤中：
13、初步处理包括边缘裁剪、污点处理、全变分图像去噪。
14、其中，在确定岩石铸体薄片照片的相对坐标后，根据岩石铸体薄片照片像素大小及数量，对岩石铸体薄片照片进行网格化处理的步骤中：
15、确定岩石铸体薄片照片的相对坐标；
16、根据岩石铸体薄片照片像素大小及数量，确定岩石铸体薄片照片纵横向的网格数量以及网格步长。
17、其中，在确定岩石铸体薄片照片的相对坐标的步骤中：
18、坐标为根据岩石铸体薄片照片的像素大小和数量确定的相对坐标。
19、其中，在确定岩石铸体薄片照片rgb分布数据，并将rgb数据转换成xyz数据的步骤中：
20、进行岩石铸体薄片照片的图像识别与处理，确定岩石铸体薄片照片的rgb平面数据分布；
21、将rgb数据转换成xyz数据。
22、其中，在基于直方图阈值分割原理，在xyz数据中，根据孔隙空间范围内z值的阈值，开展孔隙分布区域的拾取，计算孔隙分布区域的面积的步骤中：
23、确定岩石铸体薄片照片xyz数据中孔隙空间分布区域的z值的阈值区间；
24、根据z值的阈值区间，截取孔隙分布区域，计算孔隙空间分布的面积。
25、其中，在根据岩石铸体薄片照片相对面积计算岩石铸体薄片孔隙度的步骤中：
26、根据岩石铸体薄片照片确定的相对坐标计算图像的总面积；
27、根据孔隙分布面积大小，通过孔隙分布面积与岩石铸体薄片照片总面积比，计算出岩石铸体薄片的孔隙度。
28、本发明的一种基于岩石铸体薄片照片计算岩石孔隙度的方法，通过获取岩石铸体薄片照片，并对图像进行初步处理；确定岩石铸体薄片照片的相对坐标后，根据岩石铸体薄片照片像素大小及数量，对岩石铸体薄片照片进行网格化处理；将rgb数据转换成xyz数据，并确定岩石铸体薄片照片的xyz数据中，孔隙空间分布的z值的阈值区间；基于直方图阈值分割，计算孔隙分布区域的面积；根据岩石铸体薄片照片相对面积计算岩石铸体薄片孔隙度；获取岩石铸体薄片的孔隙度；通过上述方式获得精度高、准确的、定量的确定岩石铸体薄片照片中孔隙度的大小，解决了传统通过人为薄片鉴定粗略估算岩石铸体薄片照片中孔隙度大小，避免因人为经验造成误差的问题。
1.一种基于岩石铸体薄片照片计算岩石孔隙度的方法，其特征在于，包括如下步骤：
2.如权利要求1所述的基于岩石铸体薄片照片计算岩石孔隙度的方法，其特征在于，在获取岩石铸体薄片照片，并对图像进行初步处理的步骤中：
3.如权利要求2所述的基于岩石铸体薄片照片计算岩石孔隙度的方法，其特征在于，在将岩石铸体薄片照片输入计算机，进行图像初步处理的步骤中：
4.如权利要求1所述的基于岩石铸体薄片照片计算岩石孔隙度的方法，其特征在于，在确定岩石铸体薄片照片的相对坐标后，根据岩石铸体薄片照片像素大小及数量，对岩石铸体薄片照片进行网格化处理的步骤中：
5.如权利要求4所述的基于岩石铸体薄片照片计算岩石孔隙度的方法，其特征在于，在确定岩石铸体薄片照片的相对坐标的步骤中：
6.如权利要求1所述的基于岩石铸体薄片照片计算岩石孔隙度的方法，其特征在于，在确定岩石铸体薄片照片rgb分布数据，并将rgb数据转换成xyz数据的步骤中：
7.如权利要求6所述的基于岩石铸体薄片照片计算岩石孔隙度的方法，其特征在于，在基于直方图阈值分割原理，在xyz数据中，根据孔隙空间范围内z值的阈值，开展孔隙分布区域的拾取，计算孔隙分布区域的面积的步骤中：
8.如权利要求1所述的基于岩石铸体薄片照片计算岩石孔隙度的方法，其特征在于，在根据岩石铸体薄片照片相对面积计算岩石铸体薄片孔隙度的步骤中：
