一种结合中轴点和最大内切球的孔隙网络模型构建方法技术

本发明专利技术公开一种结合中轴点和最大内切球的孔隙网络模型构建方法，包括如下步骤：步骤1.基于岩心CT扫描二值化图像采用MC算法构建岩心孔隙空间表面模型；步骤2.采用LKC算法提取孔隙空间表面模型的中轴点；步骤3.采用PowerCrust算法提取孔隙空间最大内切球；步骤4.结合中轴点和最大内切球提取中轴内切球；步骤5.识别中轴内切球中的主球、从球和公共从球以定义孔隙和喉道，最终完成孔隙网络模型的构建。建。建。

【技术实现步骤摘要】
一种结合中轴点和最大内切球的孔隙网络模型构建方法


[0001]本公开属于数字岩心建模
，具体地说是涉及一种结合中轴点和最大内切球的孔隙网络模型构建方法。

技术介绍

[0002]随着社会经济的迅速发展，人类对石油和天然气的需求量日渐增多，为了提高油气藏中石油和天然气的采收率，除了在宏观尺度进行理论研究和技术开发外，还需要进行微观尺度的研究，深入其主要存储空间—岩心孔隙空间内部，因此，研究一种可以体现数字岩心孔隙空间拓扑特征的，由球体和柱体简化孔隙和喉道的几何模型构建方法是非常必要的。
[0003]现有孔隙网络模型可分为规则拓扑孔隙网络模型和真实拓扑孔隙网络模型。
[0004]规则拓扑模型的特点是孔隙和喉道分布规则，具有结构简单、容易构建等优点，但同时该模型连通性较差、无法还原真实岩心结构。而真实拓扑孔隙网络模型是指基于真实岩心二维图像构建的孔隙网络模型。构建方法主要包括多向扫描法、Voronoi多面体法、居中轴线法及最大球算法。
[0005]其中，多向扫描法通过对岩心孔隙空间进行各个方向切片、扫描来寻找孔隙和喉道，具有操作简单、计算结果较为准确的优点，但是该方法对孔隙的定义并不明确。
[0006]Voronoi多面体法通过记录其顶点为孔隙，顶点间的连线为喉道，统计孔喉相关参数获得岩心孔隙网络模型。但是，此方法只能构建以成岩过程模拟方法建立的岩心模型，适用范围较小。
[0007]居中轴线法通过燃烧算法找到孔隙空间的居中轴线，并将其表示为喉道，居中轴线的交点表示为孔隙，从而近似得到等价拓扑结构的孔隙网络模型。此方法虽然可以较准确地保留整个岩心孔隙空间的基本拓扑结构特征和形态特征，但由于生成的居中轴线节点过多，使孔隙识别结果存在不确定性。
[0008]最大球算法通过寻找孔隙空间中的任意体素点的最大内切球，并用该局部最大内切球来表示孔隙，连接此局部最大内切球的所有最小球体表示喉道。此算法构建的孔隙网络模型可有效进行孔隙的准确识别，但其生成的喉道长度比实际喉道长度偏长，因此也具有一定的局限性。

技术实现思路

[0009]针对
技术介绍
中存在的问题,为克服现有技术的不足,本专利技术公开一种结合中轴点与最大内切球的孔隙网络模型构建方法，有效反映了孔隙空间几何拓扑特性,渗流模拟效率更高,准确性更好。
[0010]结合中轴点和最大内切球的孔隙网络模型构建方法包括如下步骤：
[0011]步骤1.基于岩心CT扫描序列图像，首先经过图像处理得到二值化图像后形成数字岩心，然后采用MC算法构建孔隙空间表面模型。
[0012]步骤2.基于步骤1中所构建的孔隙空间表面模型，采用LKC算法提取孔隙空间中轴点。
[0013]步骤3.基于步骤1中所构建的孔隙空间表面模型，采用PowerCrust算法提取孔隙空间最大内切球。
[0014]步骤4.结合中轴点和最大内切球完成孔隙网络模型的构建。
[0015]进一步地，所述步骤1中采用MC算法构建孔隙空间表面模型过程为：首先，读入数字岩心所有的体数据放入集合P；其次，取集合P中相邻上下两层相近的8个点作为立方体的顶点，并将其灰度值依次与等值面阈值进行比较，如果某条边两个顶点的灰度值分别小于和大于给定阈值，则可确定等值点；然后，确定等值面与立方体的相交方式；最后，采用线性插值法计算等值面的位置坐标和法向量，并得到孔隙空间表面模型。
[0016]进一步地，所述步骤2中孔隙空间表面模型中轴点提取过程为：首先，用离散的体素点来表示岩心空间，具体包括孔隙空间体素点和固体基质体素点；其次，确定边界点。若体素点的6连通邻域至少有一个是固体基质体素点，则确定该体素点为边界点；然后，划分边界点种类。按照上、下、东、南、西、北6个方向，将边界点划分为U、B、E、S、W、N类；最后，按照U、B、E、S、W、N六个方向的顺序搜索每个方向上满足特定条件的边界点并删除，直到搜索不到满足该条件的边界点，说明中轴点已生成，否则，继续标记新的边界点。
[0017]进一步地，所述步骤3中提取孔隙空间最大内切球的方法为：首先基于步骤1构建的孔隙空间表面模型的三角网格顶点提取采样点集并进行Delaunay三角剖分和构建Voronoi图；然后，计算采样点集的极点并标记内极点和外极点；最后，基于内极点和其到对应采样点距离提取岩心空间最大内切球。
[0018]进一步地，所述步骤4中孔隙网络模型的构建方法如下：
[0019](1)提取中轴内切球；
[0020](2)识别主球、从球和公共从球；
[0021](3)定义孔隙和喉道，并完成孔隙网络模型构建。
[0022]至此，即可得到基于结合中轴点与最大内切球构建的岩心孔隙网络模型。
[0023]本专利技术的有益效果是：
[0024]本专利技术结合了居中轴点和最大内切球法的优点，准确反映了孔隙空间的拓扑特征及几何特性，适用于复杂岩心的孔隙结构特征研究，并为下一步对其进行数值分析与渗流模拟奠定了基础。
附图说明
[0025]为了更清楚地说明本专利技术实例或现有技术中的技术方案，下面将对实例或现有技术描述中所需要使用的附图作简单地介绍，显而易见，下面描述中的附图是本专利技术的一些实例，对于本领域的普通技术人员来讲，在不付出创造性劳动的前提下，还可根据这些附图获得其他附图。
[0026]图1为孔隙网络模型构建流程图；
[0027]图2为构建的孔隙空间表面模型；
[0028]图3白色部分为提取的表面模型中轴点；
[0029]图4为提取的孔隙空间最大内切球；
[0030]图5为提取的中轴内切球；
[0031]图6为中轴内切球相交的7种处理办法；
[0032]图7为结合中轴点和内切球构建的岩心孔隙网络模型。
具体实施方式
[0033]应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属
的普通技术人员通常理解的相同含义。
[0034]需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和或“包括”时，其指明存在特征、步骤、操作、器件、组件和或它们的组合。
[0035]本专利技术提出一种结合中轴点与内切球的孔隙网络模型构建方法，图1为具体操作流程。
[0036]首先，基于数字岩心采用MC算法构建如图2所示的孔隙空间表面模型。
[0037]其次，基于孔隙空间表面模型采用LKC算法通过逐层删除岩心孔隙空间体素点中满足条件的边界点以提取岩心孔隙空间的中轴点，结果如图3所示；同时，以MC算法构建的岩心孔隙空间表面模型的三角网格顶点作为采样点集，采用PowerCrust算法提取岩心孔隙空间最大内切球，结果如图4所示。
[0038]然后，结合中本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.结合中轴点和最大内切球的孔隙网络模型构建方法，其特征在于，包括如下步骤：步骤1.基于岩心CT扫描二值化图像采用MC算法构建岩心孔隙空间表面模型；步骤2.采用LKC算法提取孔隙空间表面模型的中轴点；步骤3.采用PowerCrust算法提取孔隙空间最大内切球；步骤4.结合中轴点和最大内切球提取中轴内切球；步骤5.识别中轴内切球中的主球、从球和公共从球以定义孔隙和喉道，最终完成孔隙网络模型的构建。2.如权利要求1所述的孔隙网络模型构建方法，其特征在于，所述步骤中3中最大内切球的建立过程为：(1)以步骤1中构建的孔隙空间表面模型的三角格网顶点作为采样点集进行Delaunay三角剖分；(2)基于Delaunay三角网格构建Voronoi图；(3)计算采样点集的极点，并标记内极点和外极点；(4)以内极点为球心，内极点到对应采样点的距离为半径构建内极点球，内极点球集合即为建立的岩心最大内切球。3.如权利要求1所述的孔隙网络模型构建方法，其特征在于，所述步骤5中孔隙网络模型构建方法如下：(1)中轴内切球提取；(2)主球、从球和公共从球的识别；(3)定义孔隙和喉道；(4)按照(2)处理所有的中轴内切球，直至处理完最小半径的中轴内切球时算法结束。4.如权利要求3所述的孔隙网络模型构建方法，其特征在于，所述主球、从球和公共从球的识别采用如下方法：(1)对所有中轴内切球按照半径从大到小的顺序排序；(2)标记半径最大的中轴内切球为主球，所有与该球相交且小于该球半径的球为从球；...

【专利技术属性】
技术研发人员：牛文杰，余志燕，赵元坤，宋新猛，刘雨，巩愉，杨珊珊，
申请(专利权)人：中国石油大学华东，
类型：发明
国别省市：