低渗透油藏非线性渗流机理及数值模拟方法研究.ppt.ppt

中国石油大学（华东） 3.4.2 井-网格流动方程 为了保证所提数学模型的适定性，井-网格压力方程中也必须考虑非线性渗流的影响。井所在的网格在数值模拟中处理为等效的平面径向流，考虑一个时间步内井与其所在网格间的稳定渗流，设等效半径为： 圆形油藏内为均匀介质，设该时间步内网格-井间p相的稳定流率为Qp，则任意半径处的流速为： 式中 当 时，对应于汇项。 时，对应于源项。 三、非线性渗流数模方法研究 * 中国石油大学（华东） 对上式进行积分并求取等效供液半径内的平均势，考虑到近井地带压力梯度大于临界启动压力梯度，可取a=0，并且，化简可得： 式中 表示井底流动势。 应该与网格p相流体的流动势。 由上式可得： 三、非线性渗流数模方法研究 * 中国石油大学（华东） 相等，可得： 上式即为模拟非线性渗流的井-网格方程，其中正负号对生产井取负值，注入井取正值。当用本模型模拟达西渗流时，井-网格方程与达西渗流一致。可以看出，井与网格间存在启动压力差，只有压差达到该值后才会发生流动，因此，低渗层的实际产出和注入量远低于按照达西流规律所得到的数值。 应该与网格p相流体的流动势 三、非线性渗流数模方法研究 * 中国石油大学（华东） 3.5 非线性数值模型解法研究 黑油模型中，在x方向上，网格块(i,j,k)与相邻网格块间的相传导率可定义如下： 三、非线性渗流数模方法研究 在黑油模型中， 为定值， 为弱非线性参数， 为强非线性参数。本文建立的渗流模型与黑油模型存在很大的不同，本文引入了非线性渗流修正系数修正绝对渗透率，该修正因子是与压力梯度有关的；本文还考虑了压敏效应，同样的是修正绝对渗透率，因此传导率的计算公式中的 项由黑油模型计算时的定值变成了强非线性的，普通的解法不能保证计算的稳定性和正确性，因此需采用全隐式解法。 * 中国石油大学（华东） 3.5 非线性数值模型解法研究 下面就上述提到的对绝对渗透率进行非线性渗流系数修正以及压敏效应修正的处理方式进行介绍如下： 利用有限差分法将非线性渗流微分方程离散化，其中绝对渗透率为压力梯度和压差的函数。在迭代计算中，先对方程中的压力和饱和度赋初值，然后根据赋的初值使系数线性化，迭代求解方程可得到压力和饱和度值，然后得到各个单元的压力梯度值，再由渗透率曲线得到各个单元的渗透率值。将此渗透率和刚计算出的饱和度再代入非线性渗流方程中又使其线性化。循环交替迭代，直至压力和饱和度值趋于稳定并满足计算精度要求，此值即为对应时段的压力和饱和度值，然后进入下一个时段的计算。 三、非线性渗流数模方法研究 * 中国石油大学（华东） 3.6 算例计算及对比分析 为了检验该数值模拟软件的正确性，采用该软件对不同井距的反九点面积井网概念模型进行了数值模拟计算。结合某特低渗透油藏实际储层物性和流体性质，基本油藏参数选取如下：孔隙度0.15；渗透率5×10-3μm2；a=0.8，b=14.0(1/Mpa)；地层原油相对密度0.8；粘度3；原始地层压力17.5；模型的顶面深度2000m；有效厚度15m；方向网格数为：Nx=41，Ny=41，Nz=1。工作制度以井底流压为约束条件，注水井取25MPa；采油井取10MPa；模拟时间为20年，最大时间步长为10天。应用该概念模型模拟预测了达西渗流模型、非线性渗流模型和拟启动压力梯度模型在不同井距下的油藏动态情况。 三、非线性渗流数模方法研究 * 中国石油大学（华东） 三、非线性渗流数模方法研究 a.井距150m b.井距200m c.井距250m 不同井距下压力梯度等值线图 渗流曲线上压力梯度弯曲段为 0.015-0.05MPa/m * 中国石油大学（华东） 渗流曲线上弯曲段为0.015-0.05MPa/m，可以看出，除井点处压力梯度较大外，油藏中大部分区域压力梯度比较低，处于非线性渗流曲线的弯曲段。非线性渗流模型更准确。 同时可以看出，随着井距的缩小，油藏动用范围增加，井距为150m时，压力梯度较大，边角井间都可以建立起连通关系；井距为250m时，注水井和角井间存在不流动区，只有边井见到注水效果。 三、非线性渗流数模方法研究 * 中国石油大学（华东） 3.6 算例计算及对比分析 左图可以看出，达西渗流模型的开发效果最好，在相同情况下，非线性渗流模型和拟启动压力梯度模型的开发效果比达西渗流模型差的多。200m达西渗流模型的采出程度与150m非线性渗流模型的相近。相同井距下，非线性渗流模型的采出程度比拟启动压力梯度模型的大，并