特征线法激波理论油藏驱替.docVIP

####第一篇参考文献还是不错的特征线法[3]线性偏微分方程：……………….（1） 解为，对应常微分方程组;……………….(2)解为： 即为方程(1)的特征线，其中s是任意的曲线参数，例如可以令s=t。求的全微分得到： 把方程（2）代入上式得到：因为u是方程（1）的解，所以 ，即u沿特征线方向是常数。方程（1）和方程组（2）是等价的。有了上述理论，我们可以解决更一般的拟线性、非齐次的偏微分方程：…………(3)令u=z，引入新的变量z，z是方程（3）的解，所以 。假设把方程（3）写成下面的隐式形式：微分运算得到： 解得 ，即得到，代入方程（3）得到：此方程与方程（1）有着相同的形式，可以用特征线法求解，但是特征线上不变量是，不是u。当然方程还会有其他形式，在此不再讨论。激波理论在两相驱替中的应用这是一个拟线性方程，其中h是关于x，t的函数。对h求关于t的全微分得到： 如果令，则得到一开始的方程。对比得到由激波理论得到非连续情况下的间断连接条件或称为Rankkine-Hugoniot条件：其中q为流量项，u为密度项，加减号代表间断前后。在一维油水两相驱替过程中，根据连续性方程，可以得到下面的方程式：可以看出方程为一维拟线性方程，可以用特征线法求解。当 ，即含水饱和度不随时间变化时，此时满足为常数的 曲线即为特征线。特征速度为 ，即为B-L油水驱替饱和度面运动方程。含水饱和度波既具有发散波的性质，又具有收敛波性质，是一个混合波。在激波前缘处，由间断连接条件可以得到：即：激波理论在二氧化碳驱油中的应用碳化水驱油问题中：水先用CO2饱和，然后注入地层。进入原油的CO2将减少原油粘度，使其膨胀，从而提高采收率。数学模型；设C为水中CO2浓度， 为油中CO2浓度，则有：式中：Sw为含水饱和度，u为渗流速度， 为孔隙度。此式为CO2总组分浓度的对流方程。也可以写成：方程为一阶拟线性方程，可以用特征线法求解。当 时，即CO2浓度不随时间变化。 特征速度为 ，即为CO2等浓度面的推进速度方程。对上式积分就可以得到水中CO2浓度解C(x,t)。由于CO2的吸附作用，从而造成了CO2浓度波的前缘要滞后于水驱前缘。这样就存在两个激波，第一次发生在Sw2点，称为饱和度波（S-波）；第二次发生在Sw3点，称为CO2浓度波（C-波）。图1 碳化水驱油含水饱和度分布图与常规水驱油不同，碳化水驱油具有两条分流曲线，即常规水驱油分流曲线和CO2水驱油分流曲线。由于CO2提高了水的粘度，增加了原油的相对渗透率以及降低了原油粘度，即总体上减少了水的流度，而增加了油的流度，导致含水率下降，从而使分流曲线由原始位置向右移动。图2 碳化水驱油分流函数曲线图考虑指进影响的分数流量理论[2]考虑三组分驱替模型，并假设它们分布在油相和水相中，一维状况下的质量平衡方程： 其中 方程可以写为：定义浓度速度，即固定浓度值的移动速度。 方程的全微分满足一致性条件，可以看做有两个特征速速的特征根求解问题，解为： 其中，以此类推。代表了两个组分方向：fast和slow。一致性条件[1]：为了使波保持完整性，影响波的所有变量都要有相同的波速。如图4所示。波可分为压缩波、疏散波和平稳波，压缩波满足间断性条件即为激波。这里的浓度波为疏散波。图3[1]帮助理解两个速度波。图3 油藏驱替过程中含油饱和度变化、时间和距离原理图图4 组分路径原理图图5 流量随浓度变化曲线由于粘度带来的指进影响，对分数流量理论进行了如下的改进：改变了Fi的定义，其中fi2是i组分在油相中的分数流量。 类似于Kovel理论，我们修正上式： 其中为Kovel[4]因子，H为非均质系数。混相驱两维对流扩散方程求解[3]混相驱两维对流扩散方程： 如果我们把方程中的二阶项看成是关于x，y，t的函数，那么方程就变成了一阶方程。在第一小节中对于这种方程的解法已经做了介绍。可以推导出下面的常微分方程组：通过方程的转化使方程更容易求解。参考文献[1]Helfferich, F. G. (1981, February 1). Theory of Multicomponent, Multiphase Displacement in Porous Media. Society of Petroleum Engineers. doi:10.2118/8372-PA[2] Paul, G. W., Lake, L. W., Gould, T. L. (1984, January 1). A Simplified Predictive Model for CO2 Miscible Flooding. Society of Petroleum Engineers. doi:10.2118/13238-MS[3] Garder, A. O., Peac