[发明专利]一种层状油藏非均质表皮系数计算方法

说明书

本发明公开了一种层状油藏非均质表皮系数计算方法，根据渗流理论，推导出层状油藏非均质表皮系数S_vk的计算公式，即式中N是总层数，本发明的给出了具体的层状油藏非均质表皮系数计算方法以及推导过程，证明了层状油藏非均质产生的表皮系数与各层的物性有关，层间非均质性越强，产生的表皮系数越大，但影响是呈对数关系，非均质程度的变化对产生的表皮系数影响愈来愈小。

技术领域

本发明涉及一种计算方法，具体是一种层状油藏非均质表皮系数计算方法。

背景技术

储层损害评价是试井分析与储层评价的重要内容。除储层损害产生表皮系数外，其它偏离理想渗流的因素也会产生表皮系数。有关各种非理想滤流因素产生的表皮系数计算方法有许多报道。斜井、部分打开地层、油藏边界、非达西渗流、射孔、非牛顿幕律流体、径向非均质、不渗透直线边界、各向异性、渗透率间断、定压边界、变产量、地层压力敏感的表皮系数已有计算方法。但，现在还没有针对层状油藏非均质表皮系数计算方法。

发明内容

本发明的目的在于提供一种层状油藏非均质表皮系数计算方法，以解决上述背景技术中提出的问题。

为实现上述目的，本发明提供如下技术方案：

一种层状油藏非均质表皮系数计算方法，根据渗流理论，推导出层状油藏非均质表皮系数S_vk的计算公式，即

式中N是总层数，φ——孔隙度，

其具体推导过程如下：

S1、建立物理模型，物理模型为无窜流二层模型，根据物理模型得到数学模型：

P_j|_t＝0＝P_i j＝1,2

式中：B——地层体积系数，m³/m³；c_t——综合压缩系数，1/MPa；h——地层厚度，m；k——地层渗透率，μm²；P——地层压力，MPa；P_i——原始地层压力，MPa；q——井产量，m³/d；r——极坐标，m；t——生产时间，h；μ——流体粘度，mpa·s；

上述数学模型(1)的解为：

当生产时间足够长，流体作径向流时，式(2)可简化为：

式中r_w——井半径，m；由式(3)可得：

井产量为：

式中h＝h₁+h₂，S_vk为纵向非均质表皮系数；

q＝q₁+q₂ (6)

把式(4)、式(5)代入式(6)，并取得

故纵向非均质表皮系数为

上式(7)表明：纵向非均质表皮系数与各层物性有关；

相比渗透率，孔隙度与综合压缩系数变化较小，当忽略层间孔隙度与综合压缩系数变化的影响时，则式(7)可简化为

S2、多层层状油藏非均质表皮系数