一种考虑凹凸性的岩石孔隙网络模型的孔喉截面构造方法
【技术领域】
[0001] 本发明涉及图像处理技术领域，尤其涉及一种考虑凹凸性的岩石孔隙网络模型的 孔喉截面构造方法，适用于基于孔隙网络模型进行的多孔介质重建或微观渗流模拟。
【背景技术】
[0002] 多孔介质孔隙空间的形状和连通性极不规则，拓扑结构复杂，其微观结构和物理 特性决定了许多宏观渗流性质。为了深入了解渗流机制和规律，人们通常从孔隙水平甚至 更微观的水平入手研究孔隙间的复杂渗流现象。孔隙网络模型作为现在常用的孔隙级建模 方法，具有可重复、可控制、计算速度快的优点，定量研究渗流规律简便易行。
[0003] 利用基于岩心CT切片可以对多孔介质进行三维重建，较完善的保留真实孔喉的 截面形状、拓扑结构等信息。进而可以对三维岩心做适当简化，提取其拓扑信息，并将多孔 介质抽象为具有理想几何形状的孔隙空间，建立孔隙网络模型。孔隙网络模型由喉道及其 相连的孔隙体构成，每个孔喉都有固定的截面形状，最常用的形状截面是C-T-S设置，即圆 形、三角形、正方形（PatzekTff,SilinDB.Shapefactorandhydraulicconductance innoncircularcapillaries:I.One-phasecreepingflow.Journalofcolloidand interfacescience，2001，236(2) :295-304.)。这种经典的截面形状设置保证了截面形状 因子（周长的平方与面积的比）相等，但是只能利用凸多边形进行形状等价，不能精确表征 多孔介质中大量存在的凹形截面。同时该经典的截面设置方案不能完全保证截面的周长和 面积与真实孔喉截面相等，造成水力半径（面积与周长的比）不准确。C-T-S截面设置作为 一种简单近似，可能丢失大量角隅信息，影响渗流机制的进行和渗流参数的计算。
【发明内容】
[0004] 本发明的目的是提高利用孔隙网络模型表征真实多孔介质时的精确性，从而快 速、准确地预测多孔介质渗流参数。为了达到以上目的本发明提供了一种考虑凹凸性的岩 石孔隙网络模型的孔喉截面构造方法，该方法充分考虑了多孔介质截面的凹凸性、水力半 径的等价性。
[0005] 本发明技术方案具体步骤如下：
[0006] (1)利用CT成像技术对岩心进行扫描，然后基于CT切片图像利用移动立方体法对 岩心进行三维重建。
[0007] (2)利用细化算法对三维数字岩心进行孔喉分割，将所有孔喉编号并排序，依次为 1，2, 3 · · ·N_，同时初始化N使得N= 1 ;统计每个孔喉截面周长P、面积A、形状因子G及 孔喉半径rin。
[0008] (3)利用数值实验随机构建包括凹四边形及凸四边形在内的Μ个四边形，Μ取值一 般为5000;画出Μ个四边形最大内角和形状因子的散点关系图，其轮廓曲线称之为截面最 大内角极值曲线，拟合曲线得到函数表达式。
[0009] (4)采用无因次水力半径Η判断第Ν个截面的凹凸性，根据最大内角极值曲线表达 式确定一个四边形最大内角β4。
[0010] (5)利用一条从角β4顶点出发的对角线1。将β4分为角ai和角α2两部分，同 时根据公式确定该对角线长度。
[0011] (6)建立以四条边长为未知数的截面非线性特征方程。利用牛顿迭代法求解四边 形参数，其中每条边长的迭代初值选为截面周长的四分之一，迭代结束条件为迭代前后绝 对值之差小于10 5或者迭代次数大于1000。
[0012] (7)检验截面非线性特征方程的解，若方程无正解或内角和不等于360°，返回步 骤（5)。若方程有符合物理意义的解且Ν〈Ν_，执行Ν=Ν+1，并转至步骤（4)，否则构造结 束。
[0013] 其中所述步骤（3)中的数值实验方法为：
[0014] 在直角坐标系中以原点为圆心构造一个单位圆，在单位圆范围内的四个象限分别 随机选取一点，依次连接各个点，组成一个四边形。
[0015] 所述步骤⑷中判断截面凹凸性的方法为：
[0016] H = rh/rin,
[0017] 其中rh为水力半径，rιη是孔隙半径。如果Η多0.5,则该截面形状是凸边形，否则 为凹边形。
[0018] 所述步骤（4)中的四边形最大内角公式为：
[0019]
[0020] 其中β4_表示四边形形状因子为G时最大内角的最小值，β4_表示四边形形状 因子为G时最大内角的最大值。
[0021] 所述步骤（5)中的对角线长度确定公式为：
[0022] 1。= 5A/P+k(0. 5Ρ-5Α/Ρ)，k e (〇, 1)
[0023] 其中k为符合均匀分布的随机数。
[0024] 所述步骤(6)中的截面非线性特征方程为：
[0025]
[0026] 共τ丄！、丄2、丄3、丄4ατ別衣不四边形四求边的边农。
[0027] 本发明具有以下有益效果及优点：
[0028](1)构造过程考虑了真实岩心中的凹形截面，使得多孔介质表征更准确。
[0029] (2)所构造截面不仅保证了与真实岩心截面形状因子等价，而且保证了水力半径 等价。
[0030] (3)精细的孔喉截面构造方法使得孔隙网络模型模拟得到的渗流参数更加精确。
【附图说明】
[0031] 图1为本发明的步骤流程图。
[0032] 图2为CT切片、孔隙岩石分割结果及三维重建结果。
[0033]图3为四边形数值实验构造方法及形状因子与最大内角的散点关系图。
[0034] 图4为凸形及凹形孔喉截面示意图。
[0035] 图5为截面构造示意图。
[0036] 图6为Η值不同时的截面构造结果。
【具体实施方式】
[0037] 结合附图及实施例对本发明作进一步说明：
[0038] 如图1所示，一种考虑凹凸性的岩石孔隙网络模型的孔喉截面构造方法，其步骤 如下：
[0039] (1)利用CT成像技术对岩心进行扫描，基于图像分割技术和移动立方体法对真实 岩心进行三维重建。图2依次为CT切片、孔隙岩石分割结果及三维重建结果。
[0040] ⑵利用细化算法对三维数字岩心进行孔喉分割，将所有孔喉编号并排序，依次为 1，2, 3 · · ·Ν_，同时初始化Ν使得Ν= 1 ;统计每个孔喉截面周长Ρ、面积Α、形状因子G及 孔喉半径rin。
[0041] (3)如图3(a)所示，在直角坐标系中以原点为圆心构造一个单位圆，在单位圆范 围内的四个象限分别随机选取一点A、B、C、D，依次连接各个点可组成一个四边形。如此往 复，构建5000个四边形，画出最大内角和形状因子的散点关系图如图3 (b)所示。在形状因 子G-定的