基于混合深度生成模型的岩心三维图像重建方法技术

本发明专利技术公开了一种基于混合深度生成模型的岩心三维图像重建方法，包括以下步骤：(1)对二值三维岩心CT图像进行预处理，构建训练样本集；(2)基于变分自编码器与生成对抗网络的融合网络，构建三维重建混合网络模型；(3)定义基于孔隙度的约束函数；(4)设计网络优化训练策略；(5)基于上述模型与训练策略，完成对岩心CT图像的三维重建。本发明专利技术可提高网络训练的稳定性，对均质、非均质岩心CT图像均有较好的重建效果，在石油地质领域有重要的应用价值。

【技术实现步骤摘要】
基于混合深度生成模型的岩心三维图像重建方法
本专利技术涉及一种微观结构三维重建方法，尤其涉及一种岩心微观结构三维图像重建方法，属于微观结构三维重建

技术介绍
在石油地质分析过程中，通常需要获得三维岩心结构来定量地研究渗流的微观机理，分析储集层渗流性质和运移规律。随着CT(computedtomography，电子计算机X射线断层扫描技术)分辨率的提高，CT已较多应用到岩石三维成像中，但CT的分辨率与样本尺寸是相矛盾的。为了得到准确的分析结果，需要获取高分辨率的孔隙结构图像，但扫描的样本尺寸受到很大限制，使岩样的代表性有所欠缺。因此很难得到可供分析的岩心高分辨率三维孔隙结构图像。基于数值模拟的三维多孔介质材料特征分析需要大量、多样化的三维孔隙样本以及重复的模拟实验。通过CT等物理方式直接获取岩心三维模型不仅耗时，及其经济性也不强。基于模型的随机三维重建可以有效地解决这一问题。与通过CT等硬件直接扫描重建不同，基于模型的随机重建利用有限的二维图像来重建三维结构。更具体地说，它通常使用包含在二维图像中的有限信息完成三维重建。目前研究比较深入的算法主要有：高斯随机场算法、模拟退火算法、多点统计算法等。传统的三维重建算法存在重建效率低、重建尺寸与时间复杂度呈几何倍增关系、生成样本多样性不足等问题，难以将其应用于实际工业场景。随时技术的不断发展，基于深度学习的三维重建算法越来越受到人们关注。但将深度学习技术应用于三维岩心CT图像三维重建的研究还较少，大部分还是针对二维图像进行二维图像修复。Lukas(Reconstructionofthree-dimensionalporousmediausinggenerativeadversarialneuralnetworks.PhysicalReviewE,96,043–309,2017)于2017年提出一种基于生成对抗网络(GANs)的多孔介质三维重建算法。生成对抗网络是一种深度学习模型，是近年来复杂分布上无监督学习最具前景的方法之一。生成对抗网络通过生成器(Generator)与判别器(Discriminator)的相互博弈，使两者的性能都获得提升，让生成器学习到训练样本的数据分布，从而生成与训练样本数据分布相似的样本。由于生成对抗网络的训练过程是在高维空间寻找生成器和鉴别器的纳什均衡点，所以训练过程存在不稳定，全局收敛困难等问题。变分自编码器是一种经典的生成模型，它由编码器和解码器组成。编码器通过学习训练样本数据集的内在特征，将其从二维图像空间映射到低维隐层空间。解码器的工作原理正好与编码器相反，它将隐层空间向量Z重新投射回高维图像空间，生成新的样本。但其主要缺陷是生成的样本较模糊，生成样本质量较低。同时，变分自编码器和生成对抗网络也存在各自的优点：变分自编码器性能较稳定，可以学习到训练样本的数据特征并映射到隐层空间，使隐层空间向量Z和训练样本数据分布相关联，不容易出现模式坍塌等问题。生成对抗网络可以生成清晰锐利的图片，且其生成器可以隐式无监督地学习到训练样本的数据分布，具有较强的拟合能力。另外，生成对抗网络有较好的网络兼容性，可以与经典的生成模型进行融合。将生成对抗网络单独应用于CT岩心图像的三维重建时，其稳定性，收敛效率，重建结果的多样性都存在极大的不确定性，严重影响了三维重建性能。同样地，变分自编码器虽然训练过程比较稳定，但受限于其简单的网络结构，其很难准确学习到CT岩心图像的三维特征。对于如何利用不种类型的生成网络，完成对CT岩心图像三维重建仍有许多问题需要研究解决。
技术实现思路
为解决上述问题，本专利技术将变分自编码器中的编码器与生成对抗网络进行融合，提出一种新的混合深度生成模型。该模型充分利用编码器和生成对抗网络各自的优点，提高CT岩心图像的三维重建的稳定性与生成样本的多样性，重建过程主要包括以下五个步骤：(1)对二值三维岩心CT图像进行预处理，构建训练样本集；(2)基于变分自编码器与生成对抗网络，构建三维重建混合网络模型；(3)定义基于孔隙度的约束函数；(4)设计网络优化训练策略；(5)基于上述模型与训练策略，完成对二值岩心CT图像的三维重建。附图说明图1是本专利技术中基于混合深度生成模型的岩心三维重建方法的流程框图；图2是本专利技术中用于岩心CT图像三维重建的融合网络结构框图；图3是本专利技术中编码器与生成器的融合网络结构框图；图4是本专利技术实施例中真实均质孔隙样本与生成样本的对比图，其中(a)为输入图像，(b)为真实样本的三维均质孔隙结构图及XYZ三方向截面图，(c)是为生成样本的三维均质孔隙结构图及XYZ三方向截面图；图5是本专利技术实施例中真实均质孔隙样本与生成均质孔隙样本的统计参数比较；图6是本专利技术实施例中真实非均质孔隙样本与生成非均质孔隙样本的对比图，其中(a)为输入图像，(b)为真实样本的非均质孔隙三维结构图，(c)是为生成样本的三维非均质孔隙结构图；图7是本专利技术实施例中真实非均质孔隙样本与生成非均质孔隙样本在X,Y,Z三个方向的统计参数比较。具体实施方式下面参照附图更详细地描述本专利技术的实施方式。虽然附图及下文给出了本专利技术的实施方式，但是本专利技术可以多种形式实现，而不被附图及下文描述的实施方式所限制。提供附图及下文描述的实施方式是为了能够让本专利技术更加完整及准确地被本领域技术人员所理解。图1中，基于混合深度生成模型的岩心三维重建方法，具体可以分为以下五个步骤：(1)对二值三维岩心CT图像进行预处理，构建训练样本集；(2)基于变分自编码器与生成对抗网络，构建三维重建混合网络模型；(3)定义基于孔隙度的约束函数；(4)设计网络优化训练策略；(5)基于上述模型与训练策略，完成对二值岩心CT图像的三维重建。具体地，所述步骤(1)中，由于卷积神经网络对二值数据不敏感，且训练图像尺寸决定其能否代表原始岩心CT图像的三维形态与统计特征，则所述二值三维岩心CT图像预处理是指首先对原始的二值岩心CT图像加入标准差为5的标准高斯噪声，再将其按本专利确定的尺寸约束进行缩放处理；最后，沿X,Y,Z三个方向将加噪后的原始岩心CT图像按本专利确定的递进步长切分为小尺寸的训练图像。记原始岩心CT图像的三维尺寸为N×N×N，训练图像的三维尺寸为n×n×n，为了使训练图像能够保持原始图像的局部空间结构信息，需要对原始CT图像及训练图像大小进行约束，设原始图像缩放后的大小为则约束记缩放因子为ρ，则同样地，记X,Y,Z三方向递进步长为σ个像素，则σ与n的约束关系为n＝4×σ。本专利要求训练图像大小n为2的幂次方，n≥64。在本专利技术中，作为实施实例，本专利技术的原始岩心CT图像大小为10003，规定训练图像大小为1283。首先对原始数据集加入标准差为5的标准高斯噪声；接着对原始岩心CT图像进行缩放操作，ρ＝0.5，在构建训练图像数据集时，由n＝4×σ，则σ＝32，即以32个像素为步长，沿X,Y,Z三个方向遍历，构建体素大小本文档来自技高网...

【技术保护点】
1.基于混合深度生成模型的岩心三维图像重建方法，其特征在于包括以下步骤：/n(1)对二值三维岩心CT图像进行预处理，构建训练样本集；/n(2)基于变分自编码器与生成对抗网络的融合网络，构建三维重建混合网络模型；/n(3)定义基于孔隙度的约束函数；/n(4)设计优化训练策略；/n(5)基于上述模型与训练策略，完成对岩心CT图像的三维重建。/n

【技术特征摘要】
1.基于混合深度生成模型的岩心三维图像重建方法，其特征在于包括以下步骤：
(1)对二值三维岩心CT图像进行预处理，构建训练样本集；
(2)基于变分自编码器与生成对抗网络的融合网络，构建三维重建混合网络模型；
(3)定义基于孔隙度的约束函数；
(4)设计优化训练策略；
(5)基于上述模型与训练策略，完成对岩心CT图像的三维重建。


2.根据权利要求1所述的基于混合深度生成模型的岩心三维图像重建方法，其特征在于步骤(1)所述的构建训练样本集的方式，由于卷积神经网络对二值数据不敏感，且训练图像尺寸决定其能否代表原始岩心CT图像的三维形态与统计特征，则步骤(1)构建训练样本集的特征在于首先对原始二值岩心CT图像加入标准差为5的标准高斯噪声；再将其按本专利确定的尺寸约束进行缩放处理；最后，沿X,Y,Z三个方向将加噪后的原始岩心CT图像按本专利确定的递进步长切分为小尺寸的训练图像；记原始岩心CT图像的三维尺寸为N×N×N，训练图像的三维尺寸为n×n×n，为了使训练图像能够保持原始图像的局部空间结构信息，需要对原始CT图像及训练图像大小进行约束，设原始图像缩放后的大小为则约束记缩放因子为ρ，则同样地，记X,Y,Z三方向递进步长为σ个像素，则σ与n的约束关系为n＝4×σ；本专利要求训练图像大小n为2的幂次方，n≥64。


3.根据权利要求1所述的基于混合深度生成模型的岩心三维图像重建方法，其特征在于步骤(2)所述的融合网络构建，与经典生成对抗网络不同，本发明所提出的融合模型不直接为生成器输入服从标准正态分布的隐向量，而是将训练样本中第一层图像作为输入至编码器，图像经编码器处理后，将输入图像映射至高维隐向量Zvae；另外，为保证生成样本的多样性，将编码器输出的高维隐向量Zvae与低维标准正态分布Znoise连接合并，再将整合后的隐向量整体输入至生成器,即：Ztotal＝Zvae+Znoise。

【专利技术属性】
技术研发人员：张帆，何小海，滕奇志，陈洪刚，吴小强，王正勇，
申请(专利权)人：四川大学，
类型：发明
国别省市：四川;51