1.一种基于CT扫描图像与电成像图像融合构建多尺度数字岩心方法,其特征在于,包括以下步骤:S1、由电成像图像中提取大尺度的孔隙与岩石骨架,计算其局部孔隙度分布和孔隙尺寸分布,并以此构建基于电成像图像的径向方向大尺度三维数字岩心;
S2、以CT扫描图像提取小尺度的孔隙与岩石骨架,计算其局部孔隙度分布和孔隙尺寸分布,并以此构建基于CT图像的深度纵向方向小尺度三维数字岩心;
S3、分别选取大尺度三维数字岩心和小尺度三维数字岩心的深度相同的水平切面,融合生成一个多尺度二维水平切面,并以所述多尺度二维水平切面构建多尺度三维数字岩心。
2.根据权利要求1所述的基于CT扫描图像与电成像图像融合构建多尺度数字岩心方法,其特征在于,所述步骤S1包括以下子步骤:S11、加载电成像测井资料,进行电成像测井资料预处理;
S12、采用图像修复方法生成全井眼电成像图像;
S13、在电成像测井图像上,采用图像分割技术自动提取大尺度的孔隙与岩石骨架;
S14、所述二维电成像图像通过多点地质统计方法逐层构建基于电成像图像径向方向大尺度三维数字岩心。
3.根据权利要求2所述的基于CT扫描图像与电成像图像融合构建多尺度数字岩心方法,其特征在于,所述步骤S12包括以下子步骤:S121、从待修复区域的边界点中,以优先度量准则,找到一个优先的修复点;
S122、以修复点为中心形成一个待修复块,在已知信息区域中按照待修复块与匹配块之间相似性的匹配规则,寻找最佳匹配块;
S123、将待修复块中的空白信息用最佳匹配块中对应的已知信息进行填充,更新边界,直至全部图像修复完毕,得到全井眼电成像图像。
4.根据权利要求1所述的基于CT扫描图像与电成像图像融合构建多尺度数字岩心方法,其特征在于,所述步骤S3包括以下子步骤:S31、基于CT图像,选择小尺度三维数字岩心的一个水平切面,计算其局部孔隙度分布和孔隙尺寸分布;
S32、基于电成像图像,选择大尺度三维数字岩心的一个水平切面,其深度与小尺度三维数字岩心的水平切面的深度相对,计算其局部孔隙度分布和孔隙尺寸分布;
S33、融合小尺度三维数字岩心与大尺度三维数字岩心的水平切面,生成的多尺度二维水平切面;
S34、采用多点地质统计方法,对不同深度的融合生产的多尺度二维水平切面逐层生成多尺度三维数字岩心。
5.根据权利要求4所述的基于CT扫描图像与电成像图像融合构建多尺度数字岩心方法,其特征在于,所述步骤S34包括以下子步骤:S341、选择多尺度二维水平切面作为训练图像,并对训练图像进行水平与垂直两种方向多种滤波,得到滤波值;
S342、将训练图像依据滤波值进行模式分类;
S343、如果待模拟区域中待生成点的周边像素与训练图像中对应位置的像素相似,根据训练图像的对应点像素计算待生成点处的特征值;
S344、从所有模式分类中找到一个具有相近或相似的特征值的最佳分类,从最佳分类中随机选择一个数据事件,对待生成图像进行粘贴,依据下层二维训练图像生成上层仿真图像,保证上下两层图像具有相同的统计特征;
S345、生成的仿真图像作为新的训练图像,重复递增的生成仿真图像,直至产生三维数字岩心所有的图像,最终得到多尺度三维数字岩心。
6.根据权利要求5所述的基于CT扫描图像与电成像图像融合构建多尺度数字岩心方法,其特征在于,所述多种滤波包括水平与垂直方向均值、一阶导数、二阶导数滤波。
7.根据权利要求1所述的基于CT扫描图像与电成像图像融合构建多尺度数字岩心方法,其特征在于,所述局部孔隙度分布的计算流程如下:定义局部孔隙度集合,用于保存各个像素所对应的局部孔隙度;
选择图像中的一个像素点,搜索该点邻域中的每一点,如果搜索到的点为孔隙类型,则孔隙计数加一,最终得到邻域内的孔隙类型像素点的个数;
直至图像中所有像素点全部选择完毕,则统计图像中孔隙度大小与对应的个数,得到局部孔隙度分布显示。
8.根据权利要求7所述的基于CT扫描图像与电成像图像融合构建多尺度数字岩心方法,其特征在于,局部孔隙度分布的函数定义如下:定义K(r,L)为多孔介质内部以向量r末端为中心、边长为L的立方体区域;
测量单元K(r,L)的孔隙度定义为:
φ(r,L)=V[C∩K(r,L)]/V[K(r,L)]其中,V(G)为某集合G∈Rd的体积,C表示孔隙空间;
局部孔隙度分布函数定义为:
μ(φ,L)=∑δ[φ–φ(r,L)]/m
其中,m为系统中测量单元K(r,L)的个数,δ(x)为狄拉克函数,μ(φ,L)表示边长为L。
9.根据权利要求1所述的基于CT扫描图像与电成像图像融合构建多尺度数字岩心方法,其特征在于,所述孔隙尺寸分布的计算流程如下:定义孔隙尺寸集合,用于保存各个像素的孔隙尺寸;
选择图像中的一个像素点,搜索该点邻域中的每一点,如果搜索到的点为非孔隙类型,则计算两点间的距离,并保存到孔隙尺寸集合中;
直至图像中所有像素点全部选择完毕,统计孔隙尺寸集合中每个孔隙尺寸对应的个数,得到孔隙尺寸分布显示。