一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法

1.本发明涉及油井勘探开采技术领域，特别是一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法。


背景技术：

2.室内研究和矿场实践表明，热力采油技术（如蒸汽吞吐、蒸汽驱、蒸汽辅助重力泄油等）是一项可以大幅度提高稠油油田采收率的有效技术，是目前陆上稠油开发的有效手段。渤海目前稠油热采处于规模化推广阶段，稠油热采油藏的粘度不断增加，油藏厚度不断扩大（地层条件原油粘度从50mpa
•
s到50000mpa
•
s，油藏单砂体厚度从6m扩展到40m），定量获取厚层稠油油藏注蒸汽吞吐或蒸汽驱纵向波及系数，有助于评价注蒸汽热采（包括蒸汽吞吐、蒸汽驱等）的有效加热范围、热采渗流模式，从而指导热采方案设计和注采参数优化。
3.目前关于注蒸汽的纵向波及系数，陆地油田通常在井间设置观察井，通过观察井适时监控地层中的温度，从而通过地层中的温度判断纵向上的动用情况。而海上稠油热采，由于开发成本高，通常采用大井距热采，且井间没有观察井，因此定量获取厚层稠油油藏注入蒸汽纵向波及系数成为难题。


技术实现要素：

4.本发明的目的在于克服现有技术的缺点，提供一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，解决了定量获取厚层稠油油藏注入蒸汽纵向波及系数的难题。
5.本发明的目的通过以下技术方案来实现：一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，包括以下步骤：s1、收集目标区块地质油藏基本参数；s2、计算目标区块蒸汽微元在油层垂直方向上的受力值；s3、计算目标区块蒸汽微元在油层水平方向上的受力值；s4、计算目标区块蒸汽微元在油藏垂直方向的合力和瞬时渗流速度；s5、计算目标区块蒸汽微元在油藏水平方向的合力和瞬时渗流速度；s6、计算垂向分流系数，并绘制垂向分流系数与井间距离的图版；s7、计算注入蒸汽的超覆程度，并绘制超覆程度与井间距离的图版；蒸汽超覆程度定义为累计运移到垂向上的蒸汽量与累计注入蒸汽量的比值，用d表示；s8、定量计算倾斜厚层稠油注蒸汽纵向波及系数。
6.进一步地，所述s1中的基本参数：是指收集目标区块的基础数据，包括油层厚度、地层倾角、注采井间距离、油的密度、水的密度、气体粘度、界面张力、接触角、孔隙半径、垂直方向渗透率、水平方向渗透率。为通过受力分析和渗流规律计算注入地层中蒸汽的分布做好数据准备。
7.进一步地，所述s2中计算目标区块蒸汽微元在油层垂直方向上的受力值，包括重力、浮力；
s21、重力的计算根据蒸汽微元的密度为、蒸汽微元的体积，再根据重力加速度g，计算得蒸汽微元的重力f重力；由于地层有倾角，以油层水平方向为x轴、油层垂直方向为y轴，建立直角坐标系，计算蒸汽微元的重力在x方向的分量f
x重力
、y方向的分量f
y重力
；s22、浮力的计算；在注入蒸汽突破之前，蒸汽微元受到的浮力f浮力看作是恒定的，方向为竖直向上；浮力f浮力，根据油相的密度、水相的密度、蒸汽微元的体积，重力加速度g，计算得到；由于地层有倾角，以油层水平方向为x轴、油层垂直方向为y轴，建立直角坐标系，计算浮力f
浮力
在x方向的分量f
x浮力
、y方向的分量f
y浮力
。
8.在步骤s2中，通过计算的重力和浮力值是决定注入蒸汽微元在垂向上运动的动力，由于考虑了浮力和地层倾角的影响，可使描述倾斜油藏垂直油层方向的蒸汽运动更加精确。
9.进一步地，所述s3中计算目标区块蒸汽微元在油层水平方向上的受力值，包括毛管力和驱替压力：s31、计算驱替压力，即计算注采井间主流线上的注采压力梯度，公式为其中，为生产井产量；为蒸汽注入量；为油层厚度；r为任意点到生产井的距离；p(r)为r处的压力；为r处的压力梯度；l为注采井间距离；为渗透率；为原油粘度；s32、计算毛管力对于蒸汽，受到毛管力为式中，为界面张力；为接触角；为孔隙半径；在水湿情况下，毛细管力是注入水运移的动力，取正值；在油湿情况下，毛管力是注入水运移的阻力，取负值；由于地层有倾角，根据该倾角，计算平行油藏方向的水平毛管力和垂直油藏方向的垂直毛管力。
10.在步骤s3中，通过注入蒸汽在水平方向上受力主要包括驱替力和毛管力，是决定蒸汽在水平方向上运动的动力，由于考虑了毛管力和地层倾角的影响，可使描述倾斜油藏油层水平方向的蒸汽运动更加精确。
11.进一步地，所述步骤s4计算目标区块蒸汽微元在油藏垂直方向的合力和瞬时渗流速度；s41、计算目标区块蒸汽微元在油藏垂直方向的合力
根据蒸汽微元在y方向的重力分量fy重力，蒸汽微元在y方向的重力分量fy浮力，垂直油藏方向的垂直毛细管力fy毛管力，计算得到蒸汽微元在油藏垂直方向的合力fy合力；s42、计算目标区块蒸汽微元在油藏垂直方向的瞬时渗流速度注入蒸汽沿油层斜上方渗流，可将速度分解为沿油层x方向的速度和y方向的速度；蒸汽在重力、浮力及毛管力的共同作用下，根据达西定律，计算得到y方向瞬时渗流速度，即在油藏垂直方向的瞬时渗流速度vx。
12.在步骤s4中，由于是倾斜油藏，并且s2、s3计算的各种力在油藏垂直方向均有分力，且作用方向不一致，因此需要求取s2和s3在垂直方向的合力，并根据合力和渗流理论，计算得到油藏垂直方向的瞬时渗流速度。由于考虑了多种作用力的合力和地层倾角的影响，让垂直反向的瞬时渗流速度计算更加精确。
13.进一步地，所述s5计算目标区块蒸汽微元在油藏水平方向的合力和瞬时渗流速度；s51、计算目标区块蒸汽微元在油藏水平方向的合力根据蒸汽微元在x方向的重力分量fx重力，蒸汽微元在x方向的重力分量fx浮力，驱替压力，垂直油藏方向的垂直毛细管力fx毛管力，计算得到蒸汽微元在油藏垂直方向的合力fx合力；s52、计算目标区块蒸汽微元在油藏水平方向的瞬时渗流速度根据达西定律，油滴在浮力及毛管力的共同作用下，计算得到x方向瞬时渗流速度为vy。
14.在步骤s5中，由于是倾斜油藏，第s2和s3计算的各种力在油藏水平方向均有分力，且作用方向不一致，因此需要求取s2和s3在水平方向的合力，并根据合力和渗流理论，计算得到油藏水平方向的瞬时渗流速度。由于考虑了多种作用力的合力和地层倾角的影响，让水平方向瞬时渗流速度计算更加精确。
15.进一步地，所述s6计算垂向分流系数，并绘制垂向分流系数与井间距离的图版；s61、计算垂向分流系数垂向分流系数为垂向瞬时渗流速度与总瞬时渗流速度的比值，计算方程为；s62、为纵坐标，井间距离l为横坐标，绘制垂向分流系数与井间距离l的图版。
16.在步骤s6中，根据s4、s5计算得到的垂直方向和水平方向瞬时渗流速度的计算结果，可以通过计算垂向分流系数，定量计算出注入蒸汽在垂向上的突进情况，解决了海上油田由于没有观察井，难以掌握注入蒸汽油层内分布情况的难题。
17.进一步地，所述s7计算注入蒸汽的超覆程度，并绘制超覆程度与井间距离的图版；
s71、计算注入蒸汽的超覆程度蒸汽超覆程度定义为累计运移到垂向上的蒸汽量与累计注入蒸汽量的比值，用d表示；将注采井间距离划分为n块，当蒸汽前缘运移到注采井间 i 块( 0 ＜ i ＜ n) 时，利用此处的径向蒸汽运移量与垂向流量比的乘积和第i-1块的超覆程度迭代，计算出此时的蒸汽超覆程度，计算公式为；s72、绘制超覆程度与井间距离的图版以为纵坐标，井间距离为横坐标，绘制超覆程度与注采井间距离l的图版。
18.由于s6计算的垂向分流量值是一个变化的量，在注采井间的不同位置具有不同的值，通过s7，可以描绘出注采井间不同位置处的蒸汽分布情况进一步地，所述s8定量计算倾斜厚层稠油注蒸汽纵向波及系数；s81、根据s7的图版，计算得到厚层稠油蒸汽吞吐初期不同分块的蒸汽超覆程度，计算公式为，其中，表示第块的垂向分流量；s82、通过，计算每个区块纵向波及无因次厚度；s83、用每个区块的长度乘以无因次厚度得到每个区块的面积；s84、用计算得到每个区块的每个区块面积厚度乘积；s85、计算得到面积加权后的纵向波及厚度，计算公式为；计算得到面积加权后的纵向波及系数vk，计算公式为。
19.s7计算的蒸汽超覆程度di，在注采井间的不同位置具有不同的值，通过s8，可以描绘出注采井间任何剖面的蒸汽分布面积，根据面积加权，可以定量计算出注采井间的蒸汽波及厚度。这样得到的波及厚度和波及系数值，考虑了注采井间不同位置的超覆程度差异，使得计算结果更加贴近矿场实际情况。
20.本发明具有以下优点：（1）基于注入蒸汽在地层中的受力分析，提出了通过计算注采井间距离处的纵向分流系数和不同位置的蒸汽超覆程度来预测倾斜厚层稠油油藏的纵向波及系数，能有效解决海上大井距热采且无观察井的条件下，定量预测倾斜厚层稠油注蒸汽吞吐或者蒸汽驱纵向波及系数的难题，为海上稠油油田的热采方案设计提供了支撑；（2）传统的注蒸汽的纵向波及系数，陆地油田通常在井间设置观察井，通过观察井适时监控地层中的温度，从而通过地层中的温度判断纵向上的动用情况；通俗地说，传统的注蒸汽的纵向波及系数，是通过经验情况进行实验的，最终获得相应系数；根据经验的方式，使得系数的结果非常不准确；若地质情况复杂，以经验的方式测得的数据不可靠；本发明提供的测量方式，是在对注入蒸汽在地层中受力分析的基础上，定量表征
出倾斜厚层稠油油藏注蒸汽开发不同方向上的瞬时渗流速度，根据距离注汽井点的不同距离计算出垂向分流系数，并计算得到注采井间不同位置的蒸汽超覆程度，从而定量计算出倾斜厚层稠油注入蒸汽的纵向波及系数；本发明提供的测量方式，通过计算的方式获得系数，不仅中间数据计算非常精确，而且最后经计算获得的结果也非常精确；为确定海上厚层稠油油藏注蒸汽吞吐和蒸汽驱方案设计、注采参数优化、水平井纵向位置优化等提供了非常可靠的依据；（3）本方案的测量方式，不凭借经验，而是实际计算出来的，操作性强，应用广泛。
附图说明
21.图1为定量预测倾斜厚层稠油注蒸汽纵向波及系数的主要步骤；图2为实施例1倾斜厚层稠油油藏（润湿性：油湿）地层中蒸汽微元受力示意图；图3为实施例1注采井间不同距离处的垂向分流量图版；图4为实施例1注采井间不同距离处的注入蒸汽超覆程度图版；图5为实施例2中典型油田注采井间不同距离处的垂向分流量图版；图6为实施例2中典型油田注采井间不同距离处的注入蒸汽超覆程度图版；表1为实施例1中a目标区块的基础数据；表2为实施例1中a目标区注采井间不同距离处的垂向分流量和蒸汽超覆程度计算表；表3为实施例2中b目标区块的基础数据；表4为实施2中b目标区块注采井间不同距离处的垂向分流量和蒸汽超覆程度计算表。
具体实施方式
22.下面结合附图对本发明做进一步的描述，但本发明的保护范围不局限于以下所述。
23.(实施例1)一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，主要包括以下步骤（图1）：s1，收集目标区块地质油藏基本参数等静态资料。收集目标区块的基础数据，包括油层厚度、地层原油粘度、渗透率等参数（表1），为通过受力分析和渗流规律计算注入地层中蒸汽的分布做好数据准备。表1中收集的为a目标区的基础数据。
24.表1a目标区块的基础数据s2，计算目标区块蒸汽微元在油层垂直方向上的重力和浮力值。计算的重力和浮力值是决定注入蒸汽微元在垂向上运动的动力，由于考虑了浮力和地层倾角的影响，可使
描述倾斜油藏垂直油层方向的蒸汽运动更加精确。如图2所示，厚层稠油油藏地层中某蒸汽微元主要受力有浮力、毛管力、重力和驱替压力。
25.对于蒸汽微元，受到的重力为：（a-1）式中：为蒸汽的密度，g /cm3；为蒸汽微元的体积，cm3；g 为重力加速度，m/s2。
26.由于地层具有一定倾角，以油层水平方向为x轴，油层垂直方向为y轴，建立直角坐标系统，计算重力在x方向和y方向的分量：（a-2）（a-3）在注入蒸汽突破以前，蒸汽微元受到的浮力可以近似看作是恒定的，方向为垂直向上，表达式为：（a-4）式中：为油相的密度，g /cm3；为水相的密度，g /cm3；为蒸汽微元的体积，cm3；g 为重力加速度，m/s2。
27.由于地层具有一定倾角，以油层水平方向为x轴，油层垂直方向为y轴，建立直角坐标系统，计算浮力在x方向和y方向的分量：（a-5）（a-6）根据表1的数据，油的密度取值0.95g /cm3，水的密度取值1.0g /cm3，重力加速度取值9.8 m/s2，地层倾角取值30
°
，蒸汽微元的体积取值1.0cm3。根据式子（a-2）、（a-3）、（a-4）、（a-5）计算得到单位体积（1.0cm3，）蒸汽所受到的x方向重力和浮力分别为：-4.655mn和-0.245mn；y方向的重力和浮力分别为8.063mn和-0.424mn。
28.s3，计算目标区块蒸汽微元在油层水平方向上的毛管力和驱替压力梯度值。注入蒸汽在水平方向上受力主要包括驱替力和毛管力，是决定蒸汽在水平方向上运动的动力，由于考虑了毛管力和地层倾角的影响，可使描述倾斜油藏油层水平方向的蒸汽运动更加精确。
29.用（a-7）计算注采井间主流线上的注采压力梯度。
30.（a-7）式中：为生产井产量，cm
3 /s；为蒸汽注入量，cm
3 /s；为油层厚度，m；r 为任意点到生产井的距离，cm；为 r 处的压力梯度，10-1
mpa/cm； l 为注采井间距离，cm；k 为渗透率，；为原油粘度，mpa
·
s。
31.对于蒸汽，受到的毛管力为：
（a-8）式中：为界面张力，mn/m；为接触角，
°
；为孔隙半径，mm。
32.在水湿情况下，毛细管力是注入水运移的动力，取正值；在油湿情况下，毛管力是注入水运移的阻力，取负值。毛管力的方向与注入水运移方向同向或反向，为便于计算与研究，分解为平行油藏方向的水平毛管力和垂直油藏方向的垂直毛管力。
33.（a-9）（a-10）根据表1的数据，界面张力取值0.015 n/m，接触角取值150
°
，孔隙半径取值0.00005m。由于接触角较大，可以判定为油湿，毛管力为注入水运移的阻力，取负值。根据式子（a-9）、（a-10）计算得到蒸汽所受到的x和y方向的毛管力分别为0.0015pa和-0.003pa。
34.s4，计算目标区块蒸汽微元在油藏垂直方向的合力和瞬时渗流速度。由于是倾斜油藏，s2和s3计算的各种力在油藏垂直方向均有分力，且作用方向不一致，因此需要求取s2和s3在垂直方向的合力，并根据合力和渗流理论，计算得到油藏垂直方向的瞬时渗流速度。由于考虑了多种作用力的合力和地层倾角的影响，让垂直反向的瞬时渗流速度计算更加精确。
35.联立式子（a-3）、（a-6）、（a-10），得到y方向合力计算方程：（a-11）根据表1的数据，联立（a-3）、（a-6）、（a-10），计算得到为7.635pa。
36.注入蒸汽沿油层斜上方渗流，可将速度分解为沿油层x方向的速度和y方向的速度。根据达西定律，蒸汽在重力、浮力及毛管力的共同作用下，y方向瞬时渗流速度为：（a-12）根据表1的数据，垂直方向渗透率取值500
ⅹ
10-3
µ
m2，气体粘度取值5 mpa
·
s，根据式子（a-11）、（a-12）计算得到蒸汽y方向的瞬时渗流速度为763.534
ꢀµ
m2/s。
37.s5，计算目标区块蒸汽微元在油藏水平方向的合力和瞬时渗流速度。由于是倾斜油藏，s2和s3计算的各种力在油藏水平方向均有分力，且作用方向不一致，因此需要求取s2和s3在水平方向的合力，并根据合力和渗流理论，计算得到油藏水平方向的瞬时渗流速度。由于考虑了多种作用力的合力和地层倾角的影响，让水平方向瞬时渗流速度计算更加精确。
38.联立式子（a-2）、（a-5）、（a-7）、（a-9），得到x方向合力计算方程：（a-13）将（a-2）、（a-5）、（a-5）计算结果代入（a-13），得到：（a-14）根据达西定律，油滴在浮力及毛管力的共同作用下，x方向瞬时渗流速度为：
（a-15）根据表1的数据，水平方向渗透率取值2000
ⅹ
10-3
µ
m2，气体粘度取值5mpa
·
s，根据式子（a-13）、（a-14）计算得到蒸汽x方向的瞬时渗流速度计算方程为：（a-16）s6，计算垂向分流系数，并绘制垂向分流系数与井间距离的图版。其特点在于，根据s4和s5计算得到的垂直方向和水平方向瞬时渗流速度的计算结果，可以通过计算垂向分流系数，定量计算出注入蒸汽在垂向上的突进情况，解决了海上油田由于没有观察井，难以掌握注入蒸汽油层内分布情况的难题。
39.垂向分流系数为垂向瞬时渗流速度与总瞬时渗流速度的比值，计算方程为（a-16）；以为纵坐标，井间距离l为横坐标，绘制垂向分流系数与井间距离l的图版（图3）。
40.（a-17）根据表1的数据和a区块的实际情况，为生产井产量，取值34722.22cm3/s（300m3/d）；为蒸汽注入量，取值34722.22cm3/s（300m3/d）；为油层厚度，取值40m；l为注采井间距离，取值15000cm；k为渗透率，取值；为蒸汽粘度，取值5mpa
·
s。将上述取值代入（a-17），计算方程变为：（a-18）s7，计算注入蒸汽的超覆程度，并绘制超覆程度与井间距离的图版。s6计算的垂向分流量值是一个变化的量，在注采井间的不同位置具有不同的值，通过s7，可以描绘出注采井间不同位置处的蒸汽分布情况。其特点在于，蒸汽超覆程度定义为累计运移到垂向上的蒸汽量与累计注入蒸汽量的比值，用d表示（表2）。具体计算方法为：将注采井间的距离划分为n块，当蒸汽前缘运移到注采井间第i块(0＜i＜n)时，利用此处的径向蒸汽运移量与垂向流量比的乘积和第i-1块的超覆程度di-1迭代，并用（a-19）计算出此时的蒸汽超覆程度di。并以为纵坐标，井间距离为横坐标，绘制超覆程度与注采井间距离l的图版（图4）。其中，表示第块的垂向分流量。
41.（a-19）表2a目标区注采井间不同距离处的垂向分流量和蒸汽超覆程度计算表
s8，定量计算倾斜厚层稠油注蒸汽纵向波及系数。s7计算的蒸汽超覆程度di，在注采井间的不同位置具有不同的值，通过s8，可以描绘出注采井间任何剖面的蒸汽分布面积，根据面积加权，可以定量计算出注采井间的蒸汽波及厚度。这样得到的波及厚度和波及系数值，考虑了注采井间不同位置的超覆程度差异，使得计算结果更加贴近矿场实际情况。
42.根据s7的图版，用（a-19）计算厚层稠油蒸汽吞吐初期不同分块的蒸汽超覆程度（表2第5列）；通过（1-di）计算每个区块纵向波及无因次厚度hi（表2第6列）；用每个区块的长度（10m）乘以无因次厚度1得到每个区块的面积ai（表2第7列）；用a
ihi
计算得到每个区块的每个区块面积厚度乘积（表2第8列）。
43.用（a-20）、（a-21），分别计算得到面积加权后的纵向波及厚度和纵向波及系数vk。
44.（a-20）（a-21）根据表2的数据，根据（a-20）计算得到面积加权后的无因次厚度为65.3（表2第8列求和）/150（表2第7列求和）=0.435；无因次厚度h取值为1，根据（a-20）计算得到纵向厚度波及系数为43.5%。
45.（实施例2）一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法，对某厚层倾斜稠油油藏b进行纵向波及系数预测，过程如下：s1，收集目标区块地质油藏基本参数等静态资料。收集目标区块的基础数据，包括油层厚度、地层原油粘度、渗透率等参数（表3）。表3中收集的为b目标区的基础数据。
46.表3b目标区的基础数据
s2，计算目标区块蒸汽微元在油层垂直方向上的受力值。对于蒸汽微元，受到的重力为：（a-1）式中：为蒸汽的密度，g /cm3；为蒸汽微元的体积，cm3；g 为重力加速度，m/s2。
47.由于地层具有一定倾角，以油层水平方向为x轴，油层垂直方向为y轴，建立直角坐标系统，计算重力在x方向和y方向的分量：（a-2）（a-3）在注入蒸汽突破以前，蒸汽微元受到的浮力可以近似看作是恒定的，方向为垂直向上，表达式为：（a-4）式中：为油相的密度，g /cm3；为水相的密度，g /cm3；为蒸汽微元的体积，cm3；g 为重力加速度，m/s2。
48.由于地层具有一定倾角，以油层水平方向为x轴，油层垂直方向为y轴，建立直角坐标系统，计算浮力在x方向和y方向的分量：（a-5）（a-6）根据表3的数据，油的密度取值0.95g /cm3，水的密度取值1.0g /cm3，重力加速度取值9.8 m/s2，地层倾角取值20
°
，蒸汽微元的体积取值1.0cm3。根据式子（a-2）、（a-3）、（a-4）、（a-5）计算得到单位体积（1.0cm3，）蒸汽所受到的x方向重力和浮力分别为：-3.184mn和-0.168mn；y方向的重力和浮力分别为8.749mn和-0.460mn。
49.s3，计算目标区块蒸汽微元在油层水平方向上的受力值。主要包括驱替力和毛管力，用（a-7）计算注采井间主流线上的注采压力梯度。
50.（a-7）式中：为生产井产量，cm3 /s；为蒸汽注入量，cm3 /s；为油层厚度，m；r 为任意点到生产井的距离，cm；为 r 处的压力梯度，10-1
mpa/cm； l 为注采井间距离，cm；k 为渗透率，；为原油粘度，mpa
·
s。
51.对于蒸汽，受到的毛管力为：
（a-8）式中：为界面张力，mn/m；为接触角，
°
；为孔隙半径，mm。
52.在水湿情况下，毛细管力是注入水运移的动力，取正值；在油湿情况下，毛管力是注入水运移的阻力，取负值。毛管力的方向与注入水运移方向同向或反向，为便于计算与研究，分解为平行油藏方向的水平毛管力和垂直油藏方向的垂直毛管力。
53.（a-9）（a-10）根据表3的数据，界面张力取值0.02 n/m，接触角取值100
°
，孔隙半径取值0.00005m。由于接触角较大，可以判定为油湿，毛管力为注入水运移的阻力，取负值。根据式子（a-9）、（a-10）计算得到蒸汽所受到的x和y方向的毛管力分别为0.001pa和-0.001pa。
54.s4，计算目标区块蒸汽微元在油藏垂直方向的合力和瞬时渗流速度。联立式子（a-3）、（a-6）、（a-10），得到y方向合力计算方程：（a-11）根据表3的数据，联立（a-3）、（a-6）、（a-10），计算得到为8.287pa。
55.注入蒸汽沿油层斜上方渗流，可将速度分解为沿油层x方向的速度和y方向的速度。根据达西定律，蒸汽在重力、浮力及毛管力的共同作用下，y方向瞬时渗流速度为：（a-12）根据表3的数据，垂直方向渗透率取值600
ⅹ
10-3
µ
m2，气体粘度取值5 mpa
·
s，根据式子（a-11）、（a-12）计算得到蒸汽y方向的瞬时渗流速度为994.40
ꢀµ
m2/s。
56.s5，计算目标区块蒸汽微元在油藏水平方向的合力和瞬时渗流速度。联立式子（a-2）、（a-5）、（a-7）、（a-9），得到x方向合力计算方程：（a-13）将（a-2）、（a-5）、（a-5）计算结果代入（a-13），得到：（a-14）根据达西定律，油滴在浮力及毛管力的共同作用下，x方向瞬时渗流速度为：（a-15）根据表3的数据，水平方向渗透率取值3000
ⅹ
10-3
µ
m2，气体粘度取值5 mpa
·
s，根据式子（a-13）、（a-14）计算得到蒸汽x方向的瞬时渗流速度计算方程为：（a-16）s6，计算垂向分流系数，并绘制垂向分流系数与井间距离的图版。垂向分流系数为垂向瞬时渗流速度与总瞬时渗流速度的比值，计算方程为（a-16）；以为纵坐标，井
间距离l为横坐标，绘制垂向分流系数与井间距离l的图版（图5）。
57.（a-17）根据表3的数据和a区块的实际情况，为生产井产量，取值34722.22cm3/s（300m3/d）；为蒸汽注入量，取值34722.22cm3/s（300m3/d）；为油层厚度，取值50m；l为注采井间距离，取值20000cm；k为渗透率，取值；为蒸汽粘度，取值5mpa
·
s。将上述取值代入（a-17），计算方程变为：（a-18）s7，计算注入蒸汽的超覆程度，并绘制超覆程度与井间距离的图版。蒸汽超覆程度定义为累计运移到垂向上的蒸汽量与累计注入蒸汽量的比值，用d表示（表4）。将注采井间距离划分为n块，当蒸汽前缘运移到注采井间i块(0＜i＜n)时，利用此处的径向蒸汽运移量与垂向流量比的乘积和第i-1块的超覆程度迭代，用（a-19）计算出此时的蒸汽超覆程度。并以为纵坐标，井间距离为横坐标，绘制超覆程度与注采井间距离l的图版（图6）。其中，表示第块的垂向分流量；。
58.（a-19）表4b目标区块注采井间不同距离处的垂向分流量和蒸汽超覆程度计算表s8，定量计算倾斜厚层稠油注蒸汽纵向波及系数。根据s7的图版，用（a-19）计算厚层稠油蒸汽吞吐初期不同分块的蒸汽超覆程度（表4第5列）；通过计算每个区块纵向波及无因次厚度（表4第6列）；用每个区块的长度（10m）乘以无因次厚度1得到每个区块的
面积（表4第7列）；用计算得到每个区块的每个区块面积厚度乘积（表4第8列）。
59.用（a-20）、（a-21），分别计算得到面积加权后的纵向波及厚度和纵向波及系数。
60.（a-20）（a-21）根据表4的数据，根据（a-20）计算得到面积加权后的无因次厚度为36.5（表4第8列求和）/200（表4第7列求和）=0.1825；无因次厚度h取值为1，根据（a-20）计算得到纵向厚度波及系数为18.25%。
61.上述实施例仅表达了较为优选的实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。