[发明专利]一种多相流油藏数值模拟定压边界条件的处理方法

权利要求书

1.一种多相流油藏数值模拟定压边界条件的处理方法，其特征在于，包括以下步骤：

步骤S1、根据油藏中的流体流动状态，建立多相渗流压力控制方程和饱和度控制方程，并给出定压边界条件；

步骤S2、对建立的压力控制方程和饱和度控制方程采用对应的数值离散格式，得到压力控制方程的离散格式和饱和度控制方程的离散格式；

步骤S3、针对定压力边界条件，采用对角线扩大法进行处理，求解得到压力场的分布；

步骤S4、将求得的压力场代入到边界点上的压力离散方程中，利用方程的守恒性反求边界点上的流量；

所述边界点上的压力离散方程通过在压力控制方程的离散格式的线性方程组中，将定压力边界点上的控制方程单独列出得到；

S5.将求得的边界流量代入到饱和度控制方程的离散格式中求解饱和度场。

2.根据权利要求1所述的多相流油藏数值模拟定压边界条件的处理方法，其特征在于：所述步骤S1中，所述多相包括水相、气相和油相中的至少一种，具体通过以下方式实现：

S11、建立各相的质量守恒方程；

S12、将各相的质量守恒方程相加得到最终的压力控制方程；

S13、列出定压边界条件。

3.根据权利要求1所述的多相流油藏数值模拟定压边界条件的处理方法，其特征在于：所述步骤S2具体通过以下方式实现：

S21、将计算区域进行网格离散，得到网格系统，所述网格系统包括节点和单元；

S22、在网格系统上对压力控制方程进行离散，得到压力控制方程的离散格式，所述压力控制方程的离散格式由一组以压力为未知量的线性方程组构成；

S23、在网格系统上对各相饱和度控制方程进行离散，得到各相饱和度控制方程的离散格式，所述的饱和度的离散格式由一组以饱和度为未知量的线性方程组构成。

4.根据权利要求1所述的多相流油藏数值模拟定压边界条件的处理方法，其特征在于：所述步骤S3具体包括：

S31、在压力控制方程的离散格式中，利用对角线扩大法对压力控制方程的离散格式的线性方程组进行处理，处理的线性方程组的行即为边界节点所在的行，将节点对应的质量守恒方程消除；

S32、求解经过处理后的线性方程组，得到所有网格节点上的压力值，即得到压力场。

5.根据权利要求4所述的多相流油藏数值模拟定压边界条件的处理方法，其特征在于：所述步骤S4具体通过以下方式实现：

S41、在压力的离散格式的线性方程组中，将定压力边界点上的控制方程单独列出，所述定压力边界点上的控制方程为边界点的质量守恒方程，其中包含边界点及周围节点的压力；

S42、将S3计算得到的节点上的压力值代入到S41得到的边界点的质量守恒方程中，计算该方程的残差，即为边界点上的流量值。

6.根据权利要求1所述的多相流油藏数值模拟定压边界条件的处理方法，其特征在于：所述步骤S5包括：

S51、将S4计算得到的边界点上的流量值根据节点处的流度比进行流量剖分，得到各相的流量值；

S52、将各相的流量值代入到各相饱和度控制方程的离散格式中，求解得到各相的饱和度值。