二维油藏数值模拟.pptVIP

第四章 二维油藏的数值模拟方法 二维单相流的数值模拟方法 二维油水两相流的数值模拟方法 第一节 二维单相流的数值模拟方法 一 、数学模型 1. 假设条件 1) 符合达西渗流定律 2) 等温渗流 3) 单相流体流动 4) 二维(x,y)平面流动 5) 岩石不可压缩，流体微可压缩 6) 油藏非均质和各向异性 7) 不考虑重力 2. 质量守恒方程 由三维单相质量守恒方程逐步简化到上述假设条件。 1) 三维单相非均质油藏可压缩流体和岩石 2) 二维(x,y)单相非均质油藏各向异性，可压缩流体和岩石 3) 二维(x,y)单相非均质油藏各向异性，岩石不可压缩，流体微可压缩 假设矩形或正方形边界的油藏中有一口或若干口井生产或注入，求油藏中的压力分布。 3. 初始条件和边界条件 B.C 1) 外边界 2) 内边界 式中?—点源函数 ?＝0 网格块中无井 ?＝1 网格块中有井 Pwf或Piwf 可求以下问题： 1. 定压外边界条件下 1) 内边界定产，求不同时间的油藏内压力分布和井底流压。 2)内边界定流压，求不同时间的油藏内压力分布和产量。 2. 封闭外边界条件下 1) 内边界定产，求不同时间的油藏内压力分布和井底流压。 2) 内边界定流压，求不同时间的油藏内压力分布和产量。 二、差分方程组的建立 1. 对（1）式左端第一项进行二阶差商，取不等距网格 2. 对（1）式左端第二项进行二阶差商 同（2）式可得 3. 对（1）式右端进行一阶差商 则差分方程组为 （2）＋（3）＋qv＝（4） 两端乘以 令 得隐式差分方程组为 若采用行标准排列 则（5）式可写为 （6）式可写为 式中 A— 五对角系数矩阵， 带宽 三、不同内外边界条件下的压力线性代数方程组 1. 外边界定压、内边界定产 假设取点中心网格。x方向为Nx个网格,y方向为Ny个网格. 1) 由于外边界定压，则边界压力可以不进行计算，只要计算内部各个节点（Nx-2）（Ny-2）的压力 2) 计算 3) 计算 k从内部节点按行标准排列开始算。 4) 列出系数矩阵方程，解五对角系数矩阵方程可得P。 5) 已知Qv，可利用产量方程求出Pwf。若各向同性，正方形网格，其产量公式为： 2. 外边界定压，内边界定流动压力。 网格取法同1。 1) 由于外边界定压，则边界压力可以不进行计算，只要计算内部各个节点（Nx-2）（Ny-2）的压力 2) 计算 Pi,j有两种处理方法 用n时刻的数值，即 为已知，称为显式处理，Qv 为 已知。 用n+1时刻的值，即 为 未知值，称为隐式处理。需将产量部分系数的未知数放到方程左边，已知数放在方程右边。 令 3) 计算 4) 列出系数矩阵方程，解五对角系数矩阵方程可得P。 5) 已知Pi,j和Pwf，可利用产量公式计算Qv。 3. 外边界封闭，内边界定产 假设取块中心网格。x方向为Nx个网格, y方向为Ny个网格。 1) 外边界封闭处理方法 根据达西定律，网格块之间的流动速度 因此有两种方法处理封闭边界。 为此在边界外虚拟一圈网格，并令其与边界网格压力 值相等。 加上虚拟网格后，总节点数为(Nx+2)(Ny+2) ，需求 的节点数为Nx?Ny 。 在虚拟网格中假设 当网格间用调和平均取 时，则为零，因此 2) 计算 3) 计算 4) 列出系数矩阵方程，解五对角系数矩阵方程可得P。 5) 已知Qv，可利用产量方程求出Pwf。 4. 外边界封闭，内边界定流压 网格取法和边界条件表达式同3。 1) 外边界封闭的处理方法同3。 2) 计算 3) 计算 4) 列出系数矩阵方程，解五对角系数矩阵方程可得P。 5) 已知Pi,j和Pwf，可利用产量公式计算Qv。 四、计算框图 第二节 二维油水两相流的数值模拟方法 一 、数学模型 1. 假设条件 1) 符合达西渗流定律 2) 等温渗流 3) 油水两相流动，油水两组分，油组分仅在油相中，水组分仅在水相中 4)