[发明专利]按比例增大非结构化栅格(GRID)的渗透率的方法

权利要求书

1.一种把与代表多孔介质的细分度栅格的单元相关的渗透率按比例增大为与代表所述多孔介质的非结构化粗分度栅格的单元相关的渗透率的方法，包括下述步骤：

(a)通过把粗分度栅格用作计算栅格的起源，产生面内非结构化的Voronoi计算栅格，计算栅格的单元小于粗分度栅格的单元，并且所述计算栅格及所述粗分度栅格的每个单元具有一个节点；

(b)利用与细分度栅格相关的渗透率填充计算栅格；

(c)导出计算栅格的流动方程，求解所述流动方程，并计算所述计算栅格的节点间通量和压强梯度；

(d)利用在步骤(c)中计算的通量和压强梯度，计算与粗分度栅格相关的节点间平均通量和平均压强梯度；

(e)利用在步骤(d)中计算的平均通量和平均压强梯度，计算与粗分度栅格相关的按比例增大的渗透率。

2.按照权利要求1所述的方法，其中计算栅格单元的节点间连线平行于粗分度栅格的节点间连线。

3.按照权利要求1所述的方法，其中计算栅格的单元与细分度栅格的单元的大小近似相同。

4.按照权利要求1所述的方法，其中计算栅格的单元小于细分度栅格的单元。

5.按照权利要求1所述的方法，其中在步骤(b)中填充的渗透率被分配给计算栅格的节点，分配给计算栅格的指定节点的渗透率对应于细分度栅格的某一单元的预定渗透率，如果把计算栅格叠加在细分度栅格上，则所述某一单元会包含指定节点的位置。

6.按照权利要求5所述的方法，还包括通过求取形成指定节点间连线的两个节点处的渗透率的调和平均值，计算计算栅格的指定节点间连线的渗透率。

7.按照权利要求1所述的方法，其中在步骤(b)中填充的渗透率被分配给计算栅格的指定节点间连线的中点，分配给指定节点间连线的渗透率对应于细分度栅格的某一单元的预定渗透率，如果把计算栅格叠加在细分度栅格上，则所述某一单元会包含所述节点间连线的中点。

8.按照权利要求1所述的方法，其中通过只利用关于落入计算栅格预定子域内的计算栅格节点间连线计算的通量和压强梯度，确定步骤(d)中节点间平均通量和平均压强梯度的计算。

9.按照权利要求1所述的方法，其中通过只利用在步骤(c)中，关于计算栅格的平行于指定节点间连线的节点间连线计算的通量和压强梯度，确定步骤(d)中关于粗分度栅格的指定节点间连线的平均通量和平均压强梯度的计算。

10.按照权利要求1所述的方法，其中关于粗分度栅格的指定节点，确定在步骤(e)中计算的渗透率。

11.按照权利要求1所述的方法，其中通过计算在步骤(d)中，关于指定的节点间连线计算的平均通量与平均压强梯度的比值，关于粗分度栅格的指定节点间连线，确定在步骤(e)中计算的渗透率。

12.按照权利要求1所述的方法，其中细分度栅格的单元是结构化的。

13.按照权利要求1所述的方法，其中粗分度栅格是PEBI栅格。

14.按照权利要求1所述的方法，其中粗分度栅格和计算栅格都是PEBI栅格。

15.按照权利要求1所述的方法，其中粗分度栅格的节点间连线形成Delaunay三角形，并且在步骤(a)中产生的计算栅格含有相似的较小的Delaunay三角形，其数目等于n²与大分度栅格的Delaunay三角形数目的乘积，这里n是用于产生计算栅格的预定整数细化比。

16.按照权利要求1所述的方法，其中所有单元都是三维的。

17.按照权利要求16所述的方法，其中粗分度栅格和计算栅格都是面内非结构化，垂直方向结构化的栅格。

18.按照权利要求1所述的方法，还包括利用在步骤(e)中计算的渗透率，确定粗分度栅格的节点间连线传输率。

19.按照权利要求1所述的方法，其中步骤(c)的流动方程是单相、定态方程。