CN103616494A - 一种模拟岩心孔隙空间流体分布的方法及系统 - Google Patents

Abstract

本发明实施例公开了一种模拟岩心孔隙空间流体分布的方法及系统，通过天然岩心样本，建立数字岩心模拟区域，该数字岩心模拟区域中包括第一区域和第二区域，第一区域为岩石骨架，第二区域为孔隙空间，根据天然岩心样本的含水饱和度及流体的类型，设定流体在第二区域中的初始分布和Shan-Chen模型的格子玻尔兹曼方法的演化方程的演化初始条件，基于Shan-Chen模型的格子玻尔兹曼方法的演化方程，并根据演化初始条件和流体在第二区域中的初始分布，对流体在第二区域中的分布情况进行演化计算，当演化计算达到稳定时，获取流体在第二区域中的分布情况。该方法提高了流体分布情况模拟的精确度，有利于岩石物理属性的定量研究。

Description

一种模拟岩心孔隙空间流体分布的方法及系统

技术领域

本发明涉及石油测井技术领域，特别涉及一种模拟岩心孔隙空间流体分布的方法及系统。

背景技术

由于对油气资源的需求不断增加，油气田的勘探和开发逐步由常规储层转向非常规储层，比如页岩、致密砂岩、碳酸盐岩、油砂等。岩石的物理研究在油气储层评价中所占据的地位也越来越重要，针对以上这些复杂储层，岩石物理研究遇到了诸多困难，例如：低孔、低渗岩石驱替困难，裂缝发育的碳酸盐岩难以取到代表性的岩心，页岩和油砂很难开展岩石物理实验等。而岩石物理实验无法定量研究储层微观参数对岩石宏观物理属性的影响，亟需一种储层微观参数的研究方法。

发明内容

本发明实施例中提供的一种模拟岩心孔隙空间流体分布的方法及系统，实现了可通过格子玻尔兹曼方法模拟岩石孔隙中的流体的分布情况，提供一条定量研究储层微观参数对岩石宏观物理属性影响的途径。

本发明实施例的目的是通过以下技术方案实现的：

一种模拟岩心孔隙空间流体分布的方法，包括：

根据天然岩心样本，建立数字岩心模拟区域，所述数字岩心模拟区域中包括第一区域和第二区域，所述第一区域为岩石骨架，所述第二区域为孔隙空间，其中，所述第一区域为所述数字岩心模拟区域的边界，所述第一区域和所述第二区域的体积大小为整数格子单位；

根据所述天然岩心样本的含水饱和度及流体的类型，设定流体在所述第二区域中的初始分布和Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程的演化初始条件；

基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程，并根据所述演化初始条件和流体在所述第二区域中的初始分布，对流体在所述第二区域中的分布情况进行演化计算；

当所述演化计算达到稳定时，获取流体在所述第二区域中的分布情况。

一种模拟岩心孔隙空间流体分布的系统，包括：CT扫描仪、计算机系统和图像服务器，其中，

所述CT扫描仪用于对所述天然岩心样本进行X射线扫描，以获得所述天然岩心样本的结构图像；

所述图像服务器用于根据所述CT扫描仪获得的所述天然岩心样本的结构图像，进行三维重构，生成数字岩心模拟区域的图像；

所述计算机系统用于根据所述数字岩心模拟区域的图像，建立数字岩心模拟区域，所述数字岩心模拟区域中包括第一区域和第二区域，所述第一区域为岩石骨架，所述第二区域为孔隙空间，其中，所述第一区域为所述数字岩心模拟区域的边界，所述第一区域和第二区域的体积大小为整数格子单位；根据所述天然岩心样本的含水饱和度及流体的类型，设定流体在所述第二区域中的初始分布和Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程的演化初始条件；基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程，并根据所述演化初始条件和流体在所述第二区域中的初始分布，对流体在所述第二区域中的分布情况进行演化计算；当所述演化计算达到稳定时，获取流体在所述第二区域中的分布情况；

所述图像服务器还用于根据所述计算机系统获取的流体在所述第二区域中的分布情况，生成流体在所述第二区域中的二维分布图像或生成流体在所述第二区域中的三维分布图像。

本发明实施例中，通过天然岩心样本，建立数字岩心模拟区域，所述数字岩心模拟区域中包括第一区域和第二区域，所述第一区域为岩石骨架，所述第二区域为孔隙空间，根据所述天然岩心样本的含水饱和度及流体的类型，设定流体在所述第二区域中的初始分布和Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程的演化初始条件，基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程，并根据所述演化初始条件和流体在所述第二区域中的初始分布，对流体在所述第二区域中的分布情况进行演化计算，当所述演化计算达到稳定时，获取流体在所述第二区域中的分布情况。实现了充分考虑流体自身的物理特性、流体与岩石骨架之间和流体与流体之间的相互作用的基础上，模拟岩石孔隙中的流体分布情况，能真实再现表面张力、润湿性等基本界面现象，提高了模拟精确度，有利于对岩石物理属性进行精确研究。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1是本发明实施例中数字岩心模拟区域示意图；

图2是本发明实施例岩心孔隙空间初始流体分布的等效示意图；

图3是本发明实施例中提供的一种模拟岩心孔隙空间流体分布的方法流程图。

图4是本发明实施例中提供的数字岩心内部流体分布的三维图像；

图5是本发明实施例中提供的不同时间演化步长下，数字岩心模拟区域中流体分布横截面图；

图6是本发明实施例中提供的一种模拟岩心孔隙空间流体分布的系统示意图。

具体实施方式

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

本发明实施例中，需要根据天然岩心样本，进行建模。建模需要利用计算机重构天然岩心样本的三维数字岩心，包括岩石所具有的骨架、孔隙和喉道的几何形状、大小、分布、相互连通情况等，如图1所示，模拟特定尺寸的数字岩心（例如：200*200*200），在模拟区域四周设置整数个格子单位的岩石骨架（如图1中的白色区域），岩石骨架内部为孔隙结构（如图1中的黑色区域），如图2所示为岩心孔隙空间初始流体分布的等效示意图，黑色边框为岩石骨架，黑色的球体表示孔隙，黑色的管道表示喉道。在此需要说明，图2只是示意性的说明，实际应用中的图示取决于天然岩心样本的实际孔隙结构。本发明实施例以这一模型为例，说明本发明实施例的一种模拟岩心孔隙空间流体分布的方法，流程图如图3所示，步骤如下：

步骤301、根据天然岩心样本，建立数字岩心模拟区域；

其中，为了建立实际的数字岩心模型，需要对天然岩心样本进行X射线扫描，并通过计算机进行三维重构，生成数字岩心模拟区域。所述数字岩心模拟区域中包括第一区域和第二区域，所述第一区域为岩石骨架，所述第二区域为孔隙结构，其中，所述第一区域为所述数字岩心模拟区域的边界，所述第一区域和第二区域的体积大小为整数格子单位。利用X射线计算机层析(CT)扫描仪,对天然岩心进行扫描，并进行三维重构，可获得实际的数字岩心模型，如图1所示为岩石的骨架结构，内部存在复杂的孔隙结构。为了示意性的说明岩心内部的孔隙结构，如图2展示了岩心孔隙结构的等效图，其中，空白部分的球体和管道分别表示孔隙和喉道，喉道将不同大小的孔隙进行连通。

步骤302、根据所述天然岩心样本的含水饱和度及流体的类型，设定流体在所述第二区域中的初始分布和Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程的演化初始条件；

根据所述天然岩心样本的含水饱和度在数字岩心模拟区域的孔隙结构中随机选取像素点，并把选取的像素点设定为流体，例如：如图2所示，根据所述天然岩心样本的含水饱和度在第二区域（孔隙结构）中随机选取像素点（图中孔隙中的黑点），并把选取的像素点设为液相（设定为液滴），孔隙中的其它像素点便是气相。

Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程的演化初始条件包括：数字岩心模拟区域边界的周期性边界条件，流体与固体之间的无滑移反弹边界条件，例如：假定岩石是亲水岩石，数值模拟中岩石的含水饱和度为78%、流体与固体间作用系数分别取

（水与岩石骨架），

（气与岩石骨架），流体之间的作用系数取g_f＝-0.26。

本步骤中采用LBGK模型中的DnQb模型模拟流体运动和分布情况，将流体看成多个粒子微团组成，这些粒子微团可以向若干方向任意流动。DnQb（n维空间,b个离散速度）模型是目前LB（Lattice Boltzmann，格子玻尔兹曼）方法中的常用的LBGK计算模型，DnQb模型以简单松弛过程，刻画粒子在短时间内的局部运动规律，经过反复演化逐步形成复杂多变的宏观现象，DnQb模型根据n和b的不同取值可以有多种模型例如：D1Q3、D1Q5、D2Q9、D3Q15和D3Q19等，本步骤中可以采用D3Q19模型对流体分布区域进行模拟演化。

为了引入微观分子间作用力，能够反应出流体力学的本质，本实施例中，结合Shan-Chen模型（简称SC模型）反映微观分子间作用力用格子Boltzmann方法模拟两相流，可以通过shan-chen模型引入微观分子间作用力，能够反应出流体力学的本质。根据模拟区域的流体情况，可以分别采用单组份Shan-Chen模型（例如：只包括气和水两相的情况），多组份Shan-Chen模型（例如：包括油、气和水多相的情况）。

步骤303、基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程，并根据所述演化初始条件和流体在所述第二区域中的初始分布，对流体在所述第二区域中的分布情况进行演化计算；

其中，基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程为 $f_i^k(x+e_i\mathrm{\Δt},t+\mathrm{\Δt})-f_i^k(x,t)=-\frac{1}{{\mathrm{\τ}}^k}[f_i^k(x,t)-f_i^{k\left(\mathrm{eq}\right)}(x,t)],$ τ^k为第K相流体的松弛时间，

为第K相流体在时刻t、位置x处的节点上第i方向的分布函数，

为下一个时间步长的分布函数，

为第K相流体的平衡态分布函数，根据DnQb模型，第K相流体的平衡态分布函数可以表示为 $f_i^{k\left(\mathrm{eq}\right)}={\mathrm{\ρ}}^k{\mathrm{\ω}}_i[1+\frac{e_i\·u^{k\left(\mathrm{eq}\right)}}{{C_s}^2}+\frac{{(e_i\·u^{k\left(\mathrm{eq}\right)})}^2}{{{2C}_s}^4}-\frac{u^{k\left(\mathrm{eq}\right)}\·u^{k\left(\mathrm{eq}\right)}}{{{2C}_s}^2}],$ C_s为玻尔兹曼系数，e_i为速度矢量，ω_i为权系数，u为格子声速；对于本实施例中用到的D3Q19模型，C_s为玻尔兹曼系数等于第K相流体的平衡态分布函数可以表示为 $f_i^{k\left(\mathrm{eq}\right)}={\mathrm{\ρ}}^k{\mathrm{\ω}}_i[1+{3e}_i\·u^{k\left(\mathrm{eq}\right)}+\frac{9}{2}{(e_i\·u^{k\left(\mathrm{eq}\right)})}^2-\frac{3}{2}{u}^{k\left(\mathrm{eq}\right)}\·u^{k\left(\mathrm{eq}\right)}].$

第K相流体的平衡态分布函数

中的u^k(eq)为u、

和

矢量之和，其中，u为格子声速，

表示流体与流体之间的相互作用力，表示流体与固体之间的相互作用力，τ^k为第K相流体的松弛时间，内聚力 $F_{f-f}^k(x,t)=-{\mathrm{G\ψ}}^k(x,t)\underset{i}{\mathrm{\Σ}}{\mathrm{\ω}}_i{\mathrm{\ψ}}^{k^{\′}}(x+e_i,t)e_i$ 表示第K相流体与第K’相流体之间的相互作用力，黏附力 $F_{f-s}^k(x,t)=-{G_{\mathrm{fs}}^k\mathrm{\ψ}}^k(x,t)\underset{i}{\mathrm{\Σ}}{\mathrm{\ω}}_{i}s(x+e_i)e_i$ 表示流体间的相互作用力大小，G为流体间的相互作用力系数，

为第K相流体与固体间作用力系数，ψ(x,t)为与密度有关的势函数。

也可以表示流体的润湿性，如果

表示第K相流体为非润湿相，如果

表示第K相流体为润湿相。

另外，当模拟的流体为油水系统时，由于油水密度差异不大，因此在模拟中流体重力产生的影响可忽略不计，所述u^k(eq)可以表示为u、

和

的矢量之和，当模拟的流体为气液系统时，由于气体压缩比较大、气液密度差异大，所以要考虑流体重力影响。所述第K相流体的平衡态分布函数中的所述u^k(eq)可以表示为u、

和τ^kg的矢量之和。

本步骤中，基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程进行演化计算之前，可以对演化方程中的变量增加微小变换量Δ，以触发演化方程进行演化计算。

步骤304、当所述演化计算达到稳定时，获取流体在所述第二区域中的分布情况。

其中，当所述基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程的迭代预设时间步长和/或所述数字岩心模拟区域中流体分布不发生变化和/或所述数字岩心模拟区域的含水饱和度不变时，根据迭代之后的所述基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程，生成所述第二区域中流体的分布情况。为了直观表达流体的分布情况可以通过数字图像处理技术，生成流体在所述第二区域中的二维分布图像或生成流体在所述第二区域中的三维分布图像。如图4所示，为通过计算机构建的三维流体分布图像清楚的显示孔隙空间的流体分布，图中岩石骨架设为透明，黑色区域为水和气，在实际显示过程中可以采用颜色来区分水和气，例如：蓝色代表水，红色代表气。

本发明实施例中，根据天然岩心样本，建立数字岩心模拟区域，所述数字岩心模拟区域中包括第一区域和第二区域，所述第一区域为岩石骨架，所述第二区域为孔隙空间，根据所述天然岩心样本的含水饱和度及流体的类型，设定流体在所述第二区域中的初始分布和Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程的演化初始条件，基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程，并根据所述演化初始条件和流体在所述第二区域中的初始分布，对流体在所述第二区域中的分布情况进行演化计算，当所述演化计算达到稳定时，获取流体在所述第二区域中的分布情况。实现了充分考虑流体自身的物理特性、流体与岩石骨架之间和流体与流体之间的相互作用的基础上，模拟岩石孔隙中的流体分布情况，能真实再现表面张力、润湿性等基本界面现象，提高了模拟精确度，有利于对岩石物理属性进行精确研究。

本发明实施例中，基于三维数字岩心进行岩石物理属性研究，引入格子玻尔兹曼方法以流体流动的微观机理为基础，采用粒子分布函数的演化模拟流体的特性，可以准确合理地描述物理系统的微观动力学机制，从而确定岩心孔隙空间的流体分布。将实验室利用X射线CT扫描观察得到的气和水在岩石孔隙空间的实际分布情况与图4展示的分布情况比较，发现数值模拟得到的气和水分布规律和实验得到的气和水分布规律是吻合的，说明利用本实施例中提供的技术方案确定岩石孔隙空间的气和水分布是可行的。

为更好的说明本发明实施例，基于Shan-Chen模型的格子玻尔兹曼LBGK方法进行模拟演化的过程，图5展示了不同时间演化步长下岩石孔隙空间中气和水两相密度分布图。图5展示三维数字岩心垂直方向的截面图，图中给出了不同迭代时间步长（T=100、T=6000、T=20000），岩心孔隙空间的气和水密度分布。从演化过程可以看到，气和水在表面张力作用下逐渐分离，气体逐渐融合在一块，在润湿性的作用下，气体与岩石骨架分离，聚集到孔隙中间，而水吸附在岩石孔隙壁面，形成水膜或占据孔隙空间的一些小角隅。为了把最后模拟得到的岩心用于电性研究，需通过数字图像处理技术从密度分布图像中分割出表征骨架、水以及气体的像素点，

如图6所示，本发明实施例中还提供了一种模拟岩心孔隙空间流体分布的系统，包括：CT扫描仪610、计算机系统620和图像服务器630，其中，

所述CT扫描仪610，用于对所述天然岩心样本进行X射线扫描，以获得所述天然岩心样本的结构图像；

所述图像服务器630，用于根据所述CT扫描仪获得的所述天然岩心样本的结构图像，进行三维重构，生成数字岩心模拟区域的图像；

所述计算机系统620用于根据所述数字岩心模拟区域的图像，建立数字岩心模拟区域，所述数字岩心模拟区域中包括第一区域和第二区域，所述第一区域为岩石骨架，所述第二区域为孔隙空间，其中，所述第一区域为所述数字岩心模拟区域的边界，所述第一区域和第二区域的体积大小为整数格子单位；根据所述天然岩心样本的含水饱和度及流体的类型，设定流体在所述第二区域中的初始分布和Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程的演化初始条件；基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程，并根据所述演化初始条件和流体在所述第二区域中的初始分布，对流体在所述第二区域中的分布情况进行演化计算；当所述演化计算达到稳定时，获取流体在所述第二区域中的分布情况；

所述图像服务器630还用于根据所述计算机系统620获取的流体在所述第二区域中的分布情况，生成流体在所述第二区域中的二维分布图像或生成流体在所述第二区域中的三维分布图像。

其中，所述计算机系统620可以为高性能的计算机工作站、计算机集群系统或分布式计算系统等。

所述图像服务器630为安装有三维可视化软件的服务器。

通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到本发明可借助软件加必需的硬件平台的方式来实现，当然也可以全部通过硬件来实施，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案对背景技术做出贡献的全部或者部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行本发明各个实施例或者实施例的某些部分所述的方法。

以上对本发明进行了详细介绍，本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (10)

1.一种模拟岩心孔隙空间流体分布的方法，其特征在于，包括：

根据天然岩心样本，建立数字岩心模拟区域，所述数字岩心模拟区域中包括第一区域和第二区域，所述第一区域为岩石骨架，所述第二区域为孔隙空间，其中，所述第一区域为所述数字岩心模拟区域的边界，所述第一区域和所述第二区域的体积大小为整数格子单位；

根据所述天然岩心样本的含水饱和度及流体的类型，设定流体在所述第二区域中的初始分布和Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程的演化初始条件；

基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程，并根据所述演化初始条件和流体在所述第二区域中的初始分布，对流体在所述第二区域中的分布情况进行演化计算；

当所述演化计算达到稳定时，获取流体在所述第二区域中的分布情况。

2.根据权利要求1所述的方法，其特征在于，所述演化初始条件包括流体与固体间作用系数、流体间的作用系数、模拟区域边界的周期性边界条件和流体与固体之间的无滑移反弹边界条件。

3.根据权利要求1所述的方法，其特征在于，所述根据天然岩心样本，建立数字岩心模拟区域，包括：

对所述天然岩心样本进行X射线CT扫描，并通过计算机进行三维重构，生成所述数字岩心模拟区域。

4.根据权利要求1所述的方法，其特征在于，所述基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程为

$f_i^k(x+e_i\mathrm{\Δt},t+\mathrm{\Δt})-f_i^k(x,t)=-\frac{1}{{\mathrm{\τ}}^k}[f_i^k(x,t)-f_i^{k\left(\mathrm{eq}\right)}(x,t)],$ 其中，τ^k为第K相流体的松弛时间，

为第K相流体在时刻t、位置x处的节点上第i方向的分布函数，

为下一个时间步长Δt的分布函数，

为基于DnQb模型的第K相流体的平衡态分布函数 $f_i^{k\left(\mathrm{eq}\right)}={\mathrm{\ρ}}^k{\mathrm{\ω}}_i[1+\frac{e_i\·u^{k\left(\mathrm{eq}\right)}}{{C_s}^2}+\frac{{(e_i\·u^{k\left(\mathrm{eq}\right)})}^2}{{{2C}_s}^4}-\frac{u^{k\left(\mathrm{eq}\right)}\·u^{k\left(\mathrm{eq}\right)}}{{{2C}_s}^2}],$ K为正整数，C_s为玻尔兹曼系数，e_i为速度矢量，ω_i为权系数，u为格子声速。

5.根据权利要求4所述的方法，其特征在于，所述第K相流体的平衡态分布函数中的所述u^k(eq)为u、

和矢量之和，其中，u为格子声速，内聚力 $F_{f-f}^k(x,t)=-{\mathrm{G\ψ}}^k(x,t)\underset{i}{\mathrm{\Σ}}{\mathrm{\ω}}_i{\mathrm{\ψ}}^{k^{\′}}(x+e_i,t)e_i$ 表示第K相流体与第K’相流体之间的相互作用力，黏附力表示流体间的相互作用力大小，G为流体间的相互作用力系数，

为第K相流体与固体间作用力系数，ψ(x,t)为与密度有关的势函数。

6.根据权利要求4所述的方法，其特征在于，当模拟的流体为气体和液体的系统时，所述第K相流体的平衡态分布函数中的所述u^k(eq)可以表示为u、

和τ^kg的矢量之和。

7.根据权利要求4所述的方法，其特征在于，所述DnQb模型具体为D3Q19模型。

8.根据权利要求1所述的方法，其特征在于，所述方法，还包括：

通过对所述基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程中的变量增加微小变换量Δ，触发所述基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程进行演化计算。

9.根据权利要求1所述的方法，其特征在于，所述当所述演化计算达到稳定时，获取流体在所述第二区域中的分布情况，包括：

当所述基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程的迭代预设时间步长和/或所述数字岩心模拟区域中流体分布不发生变化和/或所述数字岩心模拟区域的含水饱和度不变时，根据迭代之后的所述基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程，生成所述第二区域中流体的分布情况。

10.一种模拟岩心孔隙空间流体分布的系统，其特征在于，包括：CT扫描仪、计算机系统和图像服务器，其中，

所述CT扫描仪用于对所述天然岩心样本进行X射线扫描，以获得所述天然岩心样本的结构图像；

所述图像服务器用于根据所述CT扫描仪获得的所述天然岩心样本的结构图像，进行三维重构，生成数字岩心模拟区域的图像；

所述计算机系统用于根据所述数字岩心模拟区域的图像，建立数字岩心模拟区域，所述数字岩心模拟区域中包括第一区域和第二区域，所述第一区域为岩石骨架，所述第二区域为孔隙空间，其中，所述第一区域为所述数字岩心模拟区域的边界，所述第一区域和第二区域的体积大小为整数格子单位；根据所述天然岩心样本的含水饱和度及流体的类型，设定流体在所述第二区域中的初始分布和Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程的演化初始条件；基于Shan-Chen模型的格子玻尔兹曼LBGK方法的演化方程，并根据所述演化初始条件和流体在所述第二区域中的初始分布，对流体在所述第二区域中的分布情况进行演化计算；当所述演化计算达到稳定时，获取流体在所述第二区域中的分布情况；

所述图像服务器还用于根据所述计算机系统获取的流体在所述第二区域中的分布情况，生成流体在所述第二区域中的二维分布图像或生成流体在所述第二区域中的三维分布图像。

Images