水驱采收率低的原因.doc

线性系统内的驱替 3.1 引言 本章主要论述了前缘推进理论驱替机理1-3。主要目的是通过简单数学模型介绍EOR基本概念，这些简单模型与之后提到的复杂模型有基本特征相同。同时也研究了前缘推进理论预测线性渗流系统水驱动态应用4。前缘驱替理论还被用于粘性水驱和化学驱，比如聚合物驱和表面活性剂驱。本章介绍了注化学段塞。最后描述了一种流体驱替另一种流体时产生的扩散现象，混相流，粘性指进及其对驱替的影响。 3.2 水驱动态——前缘推进方程 3.2.1 前缘推进及相关方程。 这部分论述了线性系统水驱油过程中前缘推进理论的应用4。其中，假定岩石均质，孔隙度φ，渗透率k，长L，横截面积A，束缚水饱合度为Siw，束缚水饱合度是指不可动水饱合度，即，水相相对渗透率，krw为零。含气饱合度也为零。当注入水以前缘推进所需流速注入线性系统时，系统中含水饱合度Sw随时间变化率如方程3.1前缘推进方程或称Buckley-Leverett方程，是由Willhite4得到的。 (3.1) 其中，xSw为含水饱合度Sw所对应坐标，A为横截面积，φ为孔隙度，qt为注入量，fw为含水率，t为注入时间。 当系统水平并且忽略重力和毛管力时，含水率可通过方程3.2计算： (3.2) 其中kw为水相渗透率，ko为油相渗透率，μw为水的粘度，μo为油的粘度，fw是无因次含水率。 图3.1是含水率与含水饱合度特征曲线。含水率对含水饱合度的导数，可通过作相应饱合度下fw-Sw曲线的切线作图得到，若已知kro(Sw)和 krw(Sw)可通过计算得到。对于不同的流体岩石物性，前缘推进方程驱替前缘的饱合度不连续，含水饱合度由Siw变为Swf(驱替前缘含水饱合度)。不连续性的出现是由于含水饱合度较低(Swf)处的流速小于驱替前缘含水饱合度对应的流速，并被后者取代。 通过作Siw处含水率曲线的切线得到驱替前缘含水饱合度，如图3.1所示。切线斜率由方程3.3得到： (3.3) 其中，fwf 为驱替前缘含水饱合度对应含水率；fwf 为束缚水饱合度下的含水率；Swf为驱替前缘饱合度；Siw为束缚水饱合度。 驱替前缘或者饱合度间断点与速度的关系为： (3.4) 其中，束缚水不流动，如图3.1所示。fiw=0且vwf由下式得到： (3.5) 当含水饱合度小于Swf时，均以驱替前缘速度推进。 方程3.3也说明当驱替前缘通过多孔介质时在含水饱合度断点质量守恒，通常称为“激波理论”。详见3.3.1。 确定某特定含水饱合度对应位置是方程3.1对时间积分得到的： (3.6) 由于每个含水饱合度面对应推进速度是一定的，含水饱合度面位置与时间关系曲线是从原点出发的一系列直线。通过将下列无因次量代入方程3.6可得到含水饱合度面位置与时间的无因次关系曲线。 (3.7) 其中，xD为与注入端的无因次距离。 （3.8） 其中，tD为无因次时间（与Qi、注入流体的PV数相对应），因此公式3.6变为： （3.9） 图3.2是通过前缘推进方程预测的不同含水饱和度下无因次距离与时间的关系曲线。饱和度以相同的速度前进且位于驱替前缘（flood-front path）。驱替前缘前面的区域拥有相同的饱和度。如图3.2中xD与tD曲线中递减的斜率所示，大于Swf的饱和度前进速度逐渐减缓。有时候称这个区域为扩散波。 我们可以通过时间与距离的图像制作截面来构建饱和度剖面或者饱和度历史。一个饱和度剖面是某一固定时间下沿某一个截面所有饱和度位置的图像，如图3.2虚线所示，其中时间tD=0.15。图3.3是根据图3.2得到的tD=0.15时刻饱和度剖面。饱和度历史是指在某一xD值下，饱和度与时间的关系图。图3.4显示了在xD=1时，含水饱和度与无因次时间tD的关系；该图表明了线性系统末端含水饱和度的变化（illustrates the arrival of water saturations at the end of the linear system）。 我们通过确定水驱替区域的平均含水饱和度来确定驱替动态。在水突破（驱替前缘含水饱和度到达xD=1处）之前，只有油产出且产油量等于注水量。当水突破时时间记为tDbt，此时产油体积用PV数表示。 （3.10） 在水