稠油油藏SAGD产油速度预测方法与流程

稠油油藏sagd产油速度预测方法
技术领域
1.本发明涉及稠油油藏开发技术领域，尤其涉及稠油油藏sagd产油速度预测方法。


背景技术：

2.本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
3.世界范围内稠油资源十分丰富，占油气资源总量的三分之一以上。在常规油气勘探开发难度越来越大，能源矛盾日益突出的当下，大力发展稠油开发技术符合需求。
4.蒸汽辅助重力泄油(sagd)是高效开发稠油资源的前沿技术之一，在稠油油藏开发中得到广泛应用。sagd技术的原理是两口水平井置于油藏下方，生产井置于油藏底部，注入井在生产井上方，生产井上方形成蒸汽腔。通过注入井持续向蒸汽腔内注入蒸汽，热量通过热传导进入油藏内，被加热的原油以沿冷凝界面几乎平行的方式流向生产井，油不断通过生产井产出。
5.sagd生产过程中蒸汽腔分为上升、横向发育和下降三个不同阶段，sagd蒸汽腔在各个阶段的发育状况对稠油油藏的产量有着至关重要的影响。如何计算每个阶段的产油速度，准确表征每个阶段产量的变化特征是预测稠油油藏sagd开发产能，制定合理开发对策和提高采收率的关键。
6.随着sagd技术的不断发展，油田对预测sagd产油速度提出了更高的要求。目前，主要通过bulter模型预测sagd产油速度，但该模型只能计算蒸汽腔在部分阶段的蒸汽腔发育状态和产量，无法描述三个阶段的蒸汽腔发育状态和产油速度，而且该模型无法计算sagd过程中含油饱和度变化。因此，急需能够预测sagd上升、横向发育和下降阶段产油速度的预测方法。


技术实现要素：

7.本发明实施例提供一种稠油油藏sagd产油速度预测方法，用以提高预测的准确性，该方法包括：
8.根据油藏尺寸、蒸汽腔的扇形泄油区域最大蒸汽腔夹角、小层间蒸汽腔夹角变化、蒸汽腔的侧面泄油区域小层间夹角变化、油藏厚度和蒸汽腔前缘距离变化值将油藏分成多个小层；
9.当蒸汽腔未发育到生产井处时，计算预热结束后sagd初始阶段蒸汽腔最大扩散速度；
10.当蒸汽腔处于上升阶段时，计算扇形泄油区域和侧面泄油区域t时刻产油速度，根据扇形泄油区域和侧面泄油区域t时刻产油速度确定蒸汽腔上升阶段t时刻产油速度；
11.当蒸汽腔处于横向发育阶段时，计算蒸汽腔横向发育阶段t时刻产油速度；
12.当蒸汽腔处于下降阶段时，计算蒸汽腔下降阶段t时刻产油速度；
13.分别改变油藏厚度、油相渗透率和原油粘度，重复上述，根据预热结束后sagd初始
阶段蒸汽腔最大扩散速度、蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度和蒸汽腔下降阶段t时刻产油速度确定油藏厚度、油相渗透率和原油粘度对sagd产油速度的影响规律。
14.本发明实施例还提供一种稠油油藏sagd产油速度预测装置，用以提高预测的准确性，该装置包括：
15.油藏分层模块，用于根据油藏尺寸、蒸汽腔的扇形泄油区域最大蒸汽腔夹角、小层间蒸汽腔夹角变化、蒸汽腔的侧面泄油区域小层间夹角变化、油藏厚度和蒸汽腔前缘距离变化值将油藏分成多个小层；
16.扩散速度计算模块，用于当蒸汽腔未发育到生产井处时，计算预热结束后sagd初始阶段蒸汽腔最大扩散速度；
17.上升阶段产油速度计算模块，用于当蒸汽腔处于上升阶段时，计算扇形泄油区域和侧面泄油区域t时刻产油速度，根据扇形泄油区域和侧面泄油区域t时刻产油速度确定蒸汽腔上升阶段t时刻产油速度；
18.横向发育阶段产油速度计算模块，用于当蒸汽腔处于横向发育阶段时，计算蒸汽腔横向发育阶段t时刻产油速度；
19.下降阶段产油速度计算模块，用于当蒸汽腔处于下降阶段时，计算蒸汽腔下降阶段t时刻产油速度；
20.产油速度的影响规律分析模块，用于分别改变油藏厚度、油相渗透率和原油粘度，重复上述，根据预热结束后sagd初始阶段蒸汽腔最大扩散速度、蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度和蒸汽腔下降阶段t时刻产油速度确定油藏厚度、油相渗透率和原油粘度对sagd产油速度的影响规律。
21.本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述稠油油藏sagd产油速度预测方法。
22.本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述稠油油藏sagd产油速度预测方法的计算机程序。
23.本发明实施例中，根据油藏尺寸、蒸汽腔的扇形泄油区域最大蒸汽腔夹角、小层间蒸汽腔夹角变化、蒸汽腔的侧面泄油区域小层间夹角变化、油藏厚度和蒸汽腔前缘距离变化值将油藏分成多个小层，基于分层理论对稠油油藏sagd产油速度进行预测，可以同时预测sagd蒸汽腔上升、横向发育和下降三个阶段的产油速度，实现稠油油藏sagd全过程产油速度预测。与现有方法相比，本发明极大地降低了预测sagd产油速度所产生的人力和时间成本，提高了预测的准确性。
附图说明
24.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
25.图1为本发明实施例中稠油油藏sagd产油速度预测方法流程图一；
26.图2为本发明实施例中稠油油藏sagd产油速度预测方法流程图二；
27.图3为本发明实施例中稠油油藏sagd产油速度预测方法流程图三；
28.图4为本发明实施例中上升阶段小层划分示意图；
29.图5是本发明实施例中横向发育和下降阶段小层划分示意图；
30.图6是本发明实施例中例1中上升阶段计算结果和实验结果拟合图；
31.图7是本发明实施例中例1中横向发育阶段计算结果和实验结果拟合图；
32.图8是本发明实施例中例1中下降阶段计算结果和实验结果拟合图；
33.图9是本发明实施例中例2中上升阶段计算结果和实验结果拟合图；
34.图10是本发明实施例中例2中横向阶段计算结果和实验结果拟合图；
35.图11是本发明实施例中例2中下降阶段计算结果和实验结果拟合图；
36.图12是本发明实施例中例3中油藏厚度对产油速率影响规律图；
37.图13是本发明实施例中例3中原油粘度对产油速率影响规律图；
38.图14为本发明实施例中稠油油藏sagd产油速度预测装置结构框图一；
39.图15为本发明实施例中稠油油藏sagd产油速度预测装置结构框图二。
具体实施方式
40.为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
41.图1为本发明实施例中稠油油藏sagd产油速度预测方法流程图一，如图1所示，该方法包括：
42.步骤101：sagd过程中，蒸汽腔存在扇形和侧面泄油区域。根据油藏尺寸(指的是油藏或者室内实验油藏模型)、蒸汽腔的扇形泄油区域最大蒸汽腔夹角、小层间蒸汽腔夹角变化、蒸汽腔的侧面泄油区域小层间夹角变化、油藏厚度和蒸汽腔前缘距离变化值将油藏分成多个小层；
43.步骤102：当蒸汽腔未发育到生产井处时，计算预热结束后sagd初始阶段蒸汽腔最大扩散速度；
44.步骤103：当蒸汽腔处于上升阶段时，计算扇形泄油区域和侧面泄油区域t时刻产油速度，根据扇形泄油区域和侧面泄油区域t时刻产油速度确定蒸汽腔上升阶段t时刻产油速度；
45.步骤104：当蒸汽腔处于横向发育阶段时，计算蒸汽腔横向发育阶段t时刻产油速度；
46.步骤105：当蒸汽腔处于下降阶段时，计算蒸汽腔下降阶段t时刻产油速度；
47.步骤106：分别改变油藏厚度、油相渗透率和原油粘度，重复上述，根据预热结束后sagd初始阶段蒸汽腔最大扩散速度、蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度和蒸汽腔下降阶段t时刻产油速度确定油藏厚度、油相渗透率和原油粘度对sagd产油速度的影响规律。
48.该预测方法可以同时预测sagd蒸汽腔上升、横向发育和下降三个阶段的产油速度，实现稠油油藏sagd全过程产油速度预测。与现有方法相比，本发明极大地降低了预测
sagd产油速度所产生的人力和时间成本，提高了预测的准确性。
49.在本发明实施例中，针对现有技术中存在的“该模型无法计算sagd过程中含油饱和度变化”这一缺点，如图2所示，该方法还可以包括：
50.步骤107：根据蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度和蒸汽腔下降阶段t时刻产油速度，计算蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间δt；
51.步骤108：当第一个小层的含油饱和度达到残余油饱和度时，根据蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度、蒸汽腔下降阶段t时刻产油速度和蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间δt，计算每个小层t+δt时刻的含油饱和度。
52.具体的，确定油藏厚度、油相渗透率和原油粘度对sagd产油速度的影响规律还需要重复这两个步骤。
53.在本发明实施例中，通过以下公式计算预热结束后sagd初始阶段蒸汽腔最大扩散速度：
[0054][0055]
其中，u
max,initial
是预热结束后sagd初始阶段蒸汽腔最大扩散速度，m/s；ko为油相渗透率，m2；ρo是原油密度，kg/m3；a为经验常数，值为0.4；m为常数，值在3-5之间；g是重力加速度，m/s2；φ是孔隙度，％；δso是初始含油饱和度，％；d
start
是蒸汽腔的初始直径，m；α是油砂热扩散系数，m2/s；μ
os
是原油粘度，m2/s。
[0056]
在本发明实施例中，蒸汽腔上升阶段，根据以下公式计算扇形泄油区域t时刻产油速度：
[0057][0058][0059][0060][0061][0062]
其中，qs为扇形泄油区域t时刻产油速度，m3/s；dq
s,i
为t时刻第i小层的产油速度，m3/s；l是井的长度，m；h是油藏厚度，m；k
ro
为油相相对渗透率，m2；是流动势函数；ui是第i小层蒸汽腔移动速度，m/s；xi是第i小层距蒸汽腔前缘距离，m；μ(xi)为蒸汽腔前缘xi处油的粘度，m2/s；θ为蒸汽腔夹角，
°
；da是过流面积，m2；n为小层的总数。
[0063]
在本发明实施例中，蒸汽腔上升阶段，通过以下公式计算侧面泄油区域t时刻产油速度：
[0064][0065][0066][0067][0068]
其中，qc为侧面泄油区域t时刻产油速度，m3/s；dq
c,i
为t时刻第i小层的产油速度，m3/s；μ(xi)为第i小层距离蒸汽腔前缘距离处的原油粘度，m2/s；γ为泄油区域任意一小层的夹角，
°
；xi为第i小层距离蒸汽腔前缘的距离。
[0069]
在本发明实施例中，根据公式(6)和(10)计算蒸汽腔上升阶段t时刻产油速度：
[0070]
q1＝qs+qcꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(11)
[0071]
其中，q1为上升阶段t时刻产油速度，m3/s。
[0072]
在本发明实施例中，横向发育阶段，通过以下公式计算t时刻产油速度：
[0073][0074][0075]
其中，q2为横向发育阶段t时刻产油速度，m3/s；dq
2,i
为横向发育阶段t时刻第i小层的产油速度，m3/s。
[0076]
在本发明实施例中，下降阶段，通过以下公式计算t时刻产油速度：
[0077][0078][0079]
其中，q3为下降阶段t时刻产油速度，m3/s；dq
3,i
为下降阶段t时刻第i小层的产油速度，m3/s。
[0080]
在本发明实施例中，根据质量守恒定律，蒸汽腔从第i个小层移动到第i+1个时，第i小层的可采油已被采完。因此，蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间：
[0081][0082]
其中，δt是蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间；s
o,1
是
第一个小层的含油饱和度，％；dq
j,1
是蒸汽腔第一个小层的产油速度，m3/s，s
o,r
是残余油饱和度，％；δv为第一小层的体积，m3；j＝1,2,3；φ是孔隙度，％。
[0083]
在本发明实施例中，当第一个小层的含油饱和度达到残余油饱和度时，其他小层的含油饱和度也会相应的降低，每个小层t+δt时刻含油饱和度计算公式如下：
[0084][0085]
其中，为第i小层t+δt时刻含油饱和度；为第i小层t+δt时刻含油饱和度；δa是每个小层的横截面积，m2。
[0086]
在本发明实施例中，所述步骤103中的侧面泄油区域t时刻产油速度中的xi和步骤104中蒸汽腔横向发育阶段t时刻产油速度中的xi，通过以下公式计算第i小层距离蒸汽腔前缘的距离：
[0087]
xi＝h(tanγ-tanθ)
ꢀꢀꢀꢀꢀꢀ
(18)
[0088]
在本发明实施例中，所述步骤105中蒸汽腔下降阶段t时刻产油速度中的xi，通过以下公式计算第i小层距离蒸汽腔前缘的距离：
[0089][0090]
其中，w为油藏宽度的一半，m。
[0091]
在本发明实施例中，所述步骤103、所述步骤104和所述步骤105中，蒸汽腔移动到第i个小层时的扩展速度：
[0092][0093]
其中，δxi是与第i小层相邻两个小层之间的距离，m。
[0094]
在本发明实施例中，所述步骤103和所述步骤105中，第i个小层距离蒸汽腔前缘距离处的原油粘度：
[0095][0096]
下面举例说明本发明预测方法。
[0097]
实例1
[0098]
具体的流程图如图3所示，包括如下步骤：
[0099]
步骤1.本实例1中，以室内实验油藏模型为基础(长度、高度和厚度尺寸为35
×
21
×
3cm)，蒸汽腔上升阶段扇形泄油区域最大蒸汽腔夹角θ
max
为15
°
，小层间蒸汽腔夹角变化δθ为1
°
，侧面泄油区域小层间夹角变化δγ为1
°
，油藏厚度h为0.21m，蒸汽腔前缘距离变化值δx为0.01m，根据上述参数将整个油藏分为n个小层。蒸汽驱上升阶段见图4所示，横向发育和下降阶段如图5所示。
[0100]
步骤2.当蒸汽腔未发育到生产井处时，通过公式(1)计算预热结束后sagd初始阶段蒸汽腔最大扩散速度u
max,initial
。
[0101]
本实例1中，ko为9.44
×
10-10
m2，ρo为980kg/m3，μ
os
为104
×
10-6
m2/s，d
start
为0.08m，
孔隙度φ为0.3，δso为0.95，α为5.87
×
10-7
m2/s，计算所得u
max,initial
为4.2787
×
10-5
m/s。
[0102]
步骤3.蒸汽腔上升阶段，根据公式(2)-(6)计算扇形泄油区域t时刻产油速度qs。
[0103]
步骤4.蒸汽腔上升阶段，通过公式(7)-(10)计算侧面泄油区域t时刻产油速度qc。
[0104]
步骤5.根据公式(6)和(10)，计算蒸汽腔上升阶段t时刻产油速度q1，即采用公式(11)。
[0105]
本实例1中，根据公式(2)～(11)计算所得上升阶段产油速度如图6所示。由图6可知，随着产油速度随时间不断增加，这是由于随着蒸汽腔不断向上发育，泄油高度不断增加，重力泄油能力不断增强，蒸汽冷凝成水释放热量，原油不断被加热，所以产油速度逐渐增加，在蒸汽腔到达模型顶部后产油速度达到最大值。此外，由图6可知，计算所得产油速度结果和实验结果基本吻合，相对误差仅为4.1％，计算精度较高，表明该方法可以较为准确的预测sagd蒸汽腔上升阶段产油速度。
[0106]
步骤6.横向发育阶段，通过公式(12)-(13)计算t时刻产油速度q2。
[0107]
本实例1中，根据公式(12)～(13)计算所得横向发育阶段产油速度如图7所示。由图7可知，该阶段产油速度基本保持稳定，这是由于蒸汽腔缓慢的向两侧扩展，该阶段泄油高度已经达到最大值，重力泄油能力保持稳定，所以产油速率基本保持稳定。该计算结果与实验现象相吻合，此外，由图7可知，计算所得产油速度结果更接近与实验结果，相对误差为4.5％，比bulter模型计算结果更加准确。由此表明，本发明所提方法较bulter模型预测结果更为准确。
[0108]
步骤7.下降阶段，通过公式(14)-(15)计算t时刻产油速度q3。
[0109]
本实例1中，根据公式(14)～(15)计算所得下降阶段产油速度如图8所示。由图8可知，产油速度随着时间逐渐降低，这是由于该阶段蒸汽腔不断下降，重力泄油能力不断降低，模型上部的原油已经大量采出，所以该阶段的产油速率不断下降。此外，由图8可知，计算所得产油速度结果和实验结果基本吻合，相对误差仅为0.99％，满足工程应用的预测精度。
[0110]
步骤8.根据质量守恒定律，蒸汽腔从第i个小层移动到第i+1个时，第i小层的可采油已被采完。因此，蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间δt可以由公式(16)计算。
[0111]
步骤9.当第一个小层的含油饱和度达到残余油饱和度时，其他小层的含油饱和度也会相应的降低，每i小层t+δt时刻含油饱和度公式(17)如下。
[0112]
根据本发明优选的，分别改变油藏厚度、油相渗透率和原油粘度，重复步骤1～9，从而确定上述参数对sagd产油速度的影响规律。
[0113]
实施例2、
[0114]
与实施例1区别在于，以室内实验所用油藏模型长度、高度和厚度尺寸为62.5
×
25
×
5cm，ko为9.58
×
10-10
m2，ρo为980kg/m3，μ
os
为104
×
10-6m2/s，d
start
为0.10m，φ为0.35，δso为0.9，α为5.66
×
10-7
m2/s，计算步骤与实施例1相同。计算所得u
max,initial
为2.6731
×
10-5
m/s。
[0115]
根据公式(2)～(21)计算得到上升、横向发育和下降阶段产油速度如图9～图11所示。由图9～图11可知，计算结果和实验结果基本吻合，三个阶段计算相对误差分别为4.1％、4.5％和5.0％。由此可知本发明所述预测方法可以准确预测稠油油藏sagd三个阶段
产油速度。
[0116]
实施例3、
[0117]
与实施例1区别在于，分别改变油藏高度为15cm(或21cm)，原油粘度为实施例1中原油粘度的1.2倍，研究油藏高度和原油粘度对sagd产油速度的影响规律。计算步骤与实施例1相同。计算所得u
max,initial
为4.2580
×
10-5
m/s。
[0118]
油藏高度对sagd产油速度的影响如图12所示。由图12可知，随着油藏厚度的增加，蒸汽腔上升阶段向上发育的时间越长，到达盖层的时间也越靠后。此外，随着油藏厚度的增加，横向发育和下降阶段sagd产油越大，最大产油速度也越大。但是随着油藏厚度的增加，蒸汽腔上升阶段sagd产油速度变化不大，说明油藏厚度对蒸汽腔上升阶段的sagd产油速度影响不大，对横向发育和下降阶段的sagd产油速度影响较大。由此可知，油藏高度越大，sagd开发效果越好。
[0119]
原油粘度对sagd产油速度的影响如图13所示。由图13可知，随着原油粘度的减小，上升、横向发育和下降阶段sagd产油速度增大，最大产油速度也增大，横向阶段的时间越短，由此可知，对于原油粘度较小的油藏，sagd开发效果越好。
[0120]
本发明实施例中还提供了一种稠油油藏sagd产油速度预测装置，如下面的实施例所述。由于该装置解决问题的原理与稠油油藏sagd产油速度预测方法相似，因此该装置的实施可以参见稠油油藏sagd产油速度预测方法的实施，重复之处不再赘述。
[0121]
图14为本发明实施例中稠油油藏sagd产油速度预测装置结构框图一，如图14所示，该装置包括：
[0122]
油藏分层模块01，用于根据油藏尺寸、蒸汽腔的扇形泄油区域最大蒸汽腔夹角、小层间蒸汽腔夹角变化、蒸汽腔的侧面泄油区域小层间夹角变化、油藏厚度和蒸汽腔前缘距离变化值将油藏分成多个小层；
[0123]
扩散速度计算模块02，用于当蒸汽腔未发育到生产井处时，计算预热结束后sagd初始阶段蒸汽腔最大扩散速度；
[0124]
上升阶段产油速度计算模块03，用于当蒸汽腔处于上升阶段时，计算扇形泄油区域和侧面泄油区域t时刻产油速度，根据扇形泄油区域和侧面泄油区域t时刻产油速度确定蒸汽腔上升阶段t时刻产油速度；
[0125]
横向发育阶段产油速度计算模块04，用于当蒸汽腔处于横向发育阶段时，计算蒸汽腔横向发育阶段t时刻产油速度；
[0126]
下降阶段产油速度计算模块05，用于当蒸汽腔处于下降阶段时，计算蒸汽腔下降阶段t时刻产油速度；
[0127]
产油速度的影响规律分析模块06，用于分别改变油藏厚度、油相渗透率和原油粘度，重复上述，根据预热结束后sagd初始阶段蒸汽腔最大扩散速度、蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度和蒸汽腔下降阶段t时刻产油速度确定油藏厚度、油相渗透率和原油粘度对sagd产油速度的影响规律。
[0128]
在本发明实施例中，如图15所示，该装置还可以包括：
[0129]
时间δt计算模块07，用于根据蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度和蒸汽腔下降阶段t时刻产油速度，计算蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间δt；
[0130]
含油饱和度计算模块08，用于当第一个小层的含油饱和度达到残余油饱和度时，根据蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度、蒸汽腔下降阶段t时刻产油速度和蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间δt，计算每个小层t+δt时刻的含油饱和度。
[0131]
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述稠油油藏sagd产油速度预测方法。
[0132]
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述稠油油藏sagd产油速度预测方法的计算机程序。
[0133]
本发明实施例中，根据油藏尺寸、蒸汽腔的扇形泄油区域最大蒸汽腔夹角、小层间蒸汽腔夹角变化、蒸汽腔的侧面泄油区域小层间夹角变化、油藏厚度和蒸汽腔前缘距离变化值将油藏分成多个小层，基于分层理论对稠油油藏sagd产油速度进行预测，可以同时预测sagd蒸汽腔上升、横向发育和下降三个阶段的产油速度，实现稠油油藏sagd全过程产油速度预测。与现有方法相比，本发明极大地降低了预测sagd产油速度所产生的人力和时间成本，提高了预测的准确性。
[0134]
另外，本发明还具有如下有益效果：
[0135]
1、利用分层方法对油藏进行分层能够模拟蒸汽腔的动态发育过程，能够计算蒸汽腔扩展速率和含油饱和度分布。
[0136]
2、本发明所述预测方法所需基础参数较少，预测过程简单、易懂，可以满足现场技术人员在参数较少情况下快速预测的需求。
[0137]
3、本发明所述预测方法可用于预测不同原油性质和油藏条件下的sagd产油速度，满足了现场技术人员全面评估sagd产油速度的需求。
[0138]
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0139]
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0140]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0141]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或
其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0142]
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。