[发明专利]一种模拟岩心孔隙空间流体分布的方法及系统

权利要求书

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

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

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

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

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

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

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

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

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

fik(x+eiΔt,t+Δt)-fik(x,t)=-1τk[fik(x,t)-fik(eq)(x,t)],]]>其中，τ^k为第K相流体的松弛时间，为第K相流体在时刻t、位置x处的节点上第i方向的分布函数，为下一个时间步长Δt的分布函数，为基于DnQb模型的第K相流体的平衡态分布函数fik(eq)=ρkωi[1+ei·uk(eq)Cs2+(ei·uk(eq))22Cs4-uk(eq)·uk(eq)2Cs2],]]>K为正整数，C_s为玻尔兹曼系数，e_i为速度矢量，ω_i为权系数，u为格子声速。