1.一种基于岩芯的石油储层孔隙度测算方法,其特征在于,所述方法包括以下步骤:步骤1,钻取岩芯,制作铸体薄片;
步骤2,扫描铸体薄片获得扫描照片,计算孔隙度;
步骤3,根据孔缝关系对孔隙度进行校正;
步骤4,获得校正后的孔隙度;
其中,步骤2中,扫描铸体薄片获得扫描照片,计算孔隙度的子步骤为:使用激光共聚焦显微镜扫描岩芯的铸体薄片,获得扫描片;对扫描片进行灰度化处理;
通过边缘检测算子获得边缘曲线,以边缘曲线将扫描片分割成多个封闭区域,设置封闭区域面积的阈值作为缝隙大小阈值THRS,将各个封闭区域中面积大于缝隙大小阈值的区域构成第一区域集合;
记第一区域集合中所有封闭区域的总面积为第一面积S1,使岩芯面积S0为当前扫描片的视域面积,计得孔隙度Φ=S1/S0;如果S1大于0则使第一区域集合中面积最大的封闭区域为MAXR,MAXR的几何重心点为MAXR0;
步骤3中,根据孔缝关系对孔隙度进行校正的子步骤为:以扫描片的各个封闭区域中不属于第一区域集合的封闭区域构成孔缝集合MSet,对MSet中的封闭区域根据大小进行降序排序,以MSeti为集合中第i个封闭区域,i为正整数;
初始化变量i为2,孔缝集合MSet的大小为N,i∈[1,N];设置一个变量作为校正比CD;
步骤3.1,以扫描片中面积最大的封闭区域的几何重心点为基准点PA,使i的值增加1;
步骤3.2,封闭区域MSeti进行角点检测,如果角点数量少于3则跳转步骤3.2.1,否则跳转步骤3.2.2;
步骤3.2.1,对封闭区域MSeti进行圆识别检测到圆,检测到的圆被封闭区域所包围且半径最大化;所述圆的圆心为C和半径为R,记当前封闭区域的几何重心点为O,计算第一缝隙偏移距离CJ;
;
如果R小于CJ跳转步骤3.3,否则跳转步骤3.4且使CD的值为R÷CJ;
式中,exp()为以自然对数为底的指数函数,L()为2点的距离,A()为封闭区域的面积,点P1,P2分别是封闭区域MSeti的边缘上距离O最小和最大的点,MAXR0为MAXR的几何重心点;
步骤3.2.2,取与封闭区域MSeti的几何重心点O距离最小的角点为P1,与封闭区域MSeti的几何重心点O距离最大的角点为P2,P1和P2构造线段PL,计算第二缝隙偏移距离CH:,
式中,exp()为以自然对数为底的指数函数,L()为2点的距离,A(MSeti)为封闭区域MSeti的面积,R(MSeti)为封闭区域MSeti的外轮廓的长度;
如果第二缝隙偏移距离CH的值大于或等于线段PL的长度则跳转步骤3.3,否则跳转步骤3.4且使CD的值为 L(P1,P2)÷CH;
步骤3.3,i如果小于N则使i的值增加1并跳转步骤3.2,否则跳转步骤3.5;
步骤3.4,获得校正缝隙面积:
,
步骤3.4.1,把Area的值加入第一面积S1,i如果小于N则使i的值增加1并跳转步骤3.2,否则跳转步骤3.5;
步骤3.5,输出校正后的第一面积S1。
2.根据权利要求1所述的一种基于岩芯的石油储层孔隙度测算方法,其特征在于,步骤
1中,钻取岩芯,制作铸体薄片的子步骤为:
钻取直径为25mm的标准岩心柱样,将环氧树脂浸染剂注入标准岩心柱样,待环氧树脂浸染剂固化后磨制成岩石薄片。
3.根据权利要求1所述的一种基于岩芯的石油储层孔隙度测算方法,其特征在于,步骤
4中,获得校正后的孔隙度的子步骤为:根据孔隙度Φ=S1/S0重新计算孔隙度获得校正孔隙度,输出校正孔隙度。
4.一种基于岩芯的石油储层孔隙度测算系统,其特征在于,所述系统包括:图像获取模块:用于获取铸体薄片的激光共聚焦扫描片;
图像处理模块:处理激光共聚焦扫描获得的扫描片,输出多个封闭区域;
数据处理模块:执行如权利要求 1 3任一项权利要求所述的一种基于岩芯的石油储层~孔隙度测算方法的步骤获得孔隙度。
5.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现权利要求1 3中任一项所述一种基于岩芯的石油储层孔隙度测算方法的步骤。
~
6.一种电子设备,其特征在于,包括:存储器,其上存储有计算机程序;处理器,用于执行所述存储器中的所述计算机程序,以实现权利要求1 3中任一项所述一种基于岩芯的石~油储层孔隙度测算方法的步骤。