CN114622882B - 稠油油藏sagd产油速度预测方法 - Google Patents
稠油油藏sagd产油速度预测方法 Download PDFInfo
- Publication number
- CN114622882B CN114622882B CN202011433241.7A CN202011433241A CN114622882B CN 114622882 B CN114622882 B CN 114622882B CN 202011433241 A CN202011433241 A CN 202011433241A CN 114622882 B CN114622882 B CN 114622882B
- Authority
- CN
- China
- Prior art keywords
- oil
- steam cavity
- oil production
- stage
- sagd
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000003921 oil Substances 0.000 title claims abstract description 454
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 254
- 238000010796 Steam-assisted gravity drainage Methods 0.000 title claims abstract description 116
- 238000000034 method Methods 0.000 title claims abstract description 66
- 239000000295 fuel oil Substances 0.000 title claims abstract description 59
- 239000010410 layer Substances 0.000 claims abstract description 121
- 238000011161 development Methods 0.000 claims abstract description 67
- 239000010779 crude oil Substances 0.000 claims abstract description 41
- 230000008859 change Effects 0.000 claims abstract description 33
- 238000009792 diffusion process Methods 0.000 claims abstract description 23
- 230000035699 permeability Effects 0.000 claims abstract description 22
- 230000001174 ascending effect Effects 0.000 claims abstract description 17
- 239000011229 interlayer Substances 0.000 claims abstract description 17
- 230000000630 rising effect Effects 0.000 claims description 38
- 238000004364 calculation method Methods 0.000 claims description 27
- 238000004590 computer program Methods 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 238000003860 storage Methods 0.000 claims description 7
- 230000005484 gravity Effects 0.000 claims description 6
- 239000003027 oil sand Substances 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 11
- 230000000694 effects Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 239000010734 process oil Substances 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 241000237503 Pectinidae Species 0.000 description 1
- 238000010795 Steam Flooding Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 230000005494 condensation Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 235000020637 scallop Nutrition 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
- E21B43/24—Enhanced recovery methods for obtaining hydrocarbons using heat, e.g. steam injection
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/14—Obtaining from a multiple-zone well
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Environmental & Geological Engineering (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Computer Hardware Design (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Geometry (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种稠油油藏SAGD产油速度预测方法,该方法包括:根据油藏尺寸、蒸汽腔的扇形泄油区域最大蒸汽腔夹角、小层间蒸汽腔夹角变化、蒸汽腔的侧面泄油区域小层间夹角变化、油藏厚度和蒸汽腔前缘距离变化值将油藏分成多个小层;当蒸汽腔未发育到生产井处时,计算预热结束后SAGD初始阶段蒸汽腔最大扩散速度;确定蒸汽腔上升阶段、横向发育阶段、下降阶段的t时刻产油速度;分别改变油藏厚度、油相渗透率和原油粘度,重复上述,确定油藏厚度、油相渗透率和原油粘度对SAGD产油速度的影响规律。本发明提高了预测的准确性。
Description
技术领域
本发明涉及稠油油藏开发技术领域,尤其涉及稠油油藏SAGD产油速度预测方法。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
世界范围内稠油资源十分丰富,占油气资源总量的三分之一以上。在常规油气勘探开发难度越来越大,能源矛盾日益突出的当下,大力发展稠油开发技术符合需求。
蒸汽辅助重力泄油(SAGD)是高效开发稠油资源的前沿技术之一,在稠油油藏开发中得到广泛应用。SAGD技术的原理是两口水平井置于油藏下方,生产井置于油藏底部,注入井在生产井上方,生产井上方形成蒸汽腔。通过注入井持续向蒸汽腔内注入蒸汽,热量通过热传导进入油藏内,被加热的原油以沿冷凝界面几乎平行的方式流向生产井,油不断通过生产井产出。
SAGD生产过程中蒸汽腔分为上升、横向发育和下降三个不同阶段,SAGD蒸汽腔在各个阶段的发育状况对稠油油藏的产量有着至关重要的影响。如何计算每个阶段的产油速度,准确表征每个阶段产量的变化特征是预测稠油油藏SAGD开发产能,制定合理开发对策和提高采收率的关键。
随着SAGD技术的不断发展,油田对预测SAGD产油速度提出了更高的要求。目前,主要通过Bulter模型预测SAGD产油速度,但该模型只能计算蒸汽腔在部分阶段的蒸汽腔发育状态和产量,无法描述三个阶段的蒸汽腔发育状态和产油速度,而且该模型无法计算SAGD过程中含油饱和度变化。因此,急需能够预测SAGD上升、横向发育和下降阶段产油速度的预测方法。
发明内容
本发明实施例提供一种稠油油藏SAGD产油速度预测方法,用以提高预测的准确性,该方法包括:
根据油藏尺寸、蒸汽腔的扇形泄油区域最大蒸汽腔夹角、小层间蒸汽腔夹角变化、蒸汽腔的侧面泄油区域小层间夹角变化、油藏厚度和蒸汽腔前缘距离变化值将油藏分成多个小层;
当蒸汽腔未发育到生产井处时,计算预热结束后SAGD初始阶段蒸汽腔最大扩散速度;
当蒸汽腔处于上升阶段时,计算扇形泄油区域和侧面泄油区域t时刻产油速度,根据扇形泄油区域和侧面泄油区域t时刻产油速度确定蒸汽腔上升阶段t时刻产油速度;
当蒸汽腔处于横向发育阶段时,计算蒸汽腔横向发育阶段t时刻产油速度;
当蒸汽腔处于下降阶段时,计算蒸汽腔下降阶段t时刻产油速度;
分别改变油藏厚度、油相渗透率和原油粘度,重复上述步骤,根据预热结束后SAGD初始阶段蒸汽腔最大扩散速度、蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度和蒸汽腔下降阶段t时刻产油速度确定油藏厚度、油相渗透率和原油粘度对SAGD产油速度的影响规律。
本发明实施例还提供一种稠油油藏SAGD产油速度预测装置,用以提高预测的准确性,该装置包括:
油藏分层模块,用于根据油藏尺寸、蒸汽腔的扇形泄油区域最大蒸汽腔夹角、小层间蒸汽腔夹角变化、蒸汽腔的侧面泄油区域小层间夹角变化、油藏厚度和蒸汽腔前缘距离变化值将油藏分成多个小层;
扩散速度计算模块,用于当蒸汽腔未发育到生产井处时,计算预热结束后SAGD初始阶段蒸汽腔最大扩散速度;
上升阶段产油速度计算模块,用于当蒸汽腔处于上升阶段时,计算扇形泄油区域和侧面泄油区域t时刻产油速度,根据扇形泄油区域和侧面泄油区域t时刻产油速度确定蒸汽腔上升阶段t时刻产油速度;
横向发育阶段产油速度计算模块,用于当蒸汽腔处于横向发育阶段时,计算蒸汽腔横向发育阶段t时刻产油速度;
下降阶段产油速度计算模块,用于当蒸汽腔处于下降阶段时,计算蒸汽腔下降阶段t时刻产油速度;
产油速度的影响规律分析模块,用于分别改变油藏厚度、油相渗透率和原油粘度,重复上述步骤,根据预热结束后SAGD初始阶段蒸汽腔最大扩散速度、蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度和蒸汽腔下降阶段t时刻产油速度确定油藏厚度、油相渗透率和原油粘度对SAGD产油速度的影响规律。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述稠油油藏SAGD产油速度预测方法。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述稠油油藏SAGD产油速度预测方法的计算机程序。
本发明实施例中,根据油藏尺寸、蒸汽腔的扇形泄油区域最大蒸汽腔夹角、小层间蒸汽腔夹角变化、蒸汽腔的侧面泄油区域小层间夹角变化、油藏厚度和蒸汽腔前缘距离变化值将油藏分成多个小层,基于分层理论对稠油油藏SAGD产油速度进行预测,可以同时预测SAGD蒸汽腔上升、横向发育和下降三个阶段的产油速度,实现稠油油藏SAGD全过程产油速度预测。与现有方法相比,本发明极大地降低了预测SAGD产油速度所产生的人力和时间成本,提高了预测的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例中稠油油藏SAGD产油速度预测方法流程图一;
图2为本发明实施例中稠油油藏SAGD产油速度预测方法流程图二;
图3为本发明实施例中稠油油藏SAGD产油速度预测方法流程图三;
图4为本发明实施例中上升阶段小层划分示意图;
图5是本发明实施例中横向发育和下降阶段小层划分示意图;
图6是本发明实施例中例1中上升阶段计算结果和实验结果拟合图;
图7是本发明实施例中例1中横向发育阶段计算结果和实验结果拟合图;
图8是本发明实施例中例1中下降阶段计算结果和实验结果拟合图;
图9是本发明实施例中例2中上升阶段计算结果和实验结果拟合图;
图10是本发明实施例中例2中横向阶段计算结果和实验结果拟合图;
图11是本发明实施例中例2中下降阶段计算结果和实验结果拟合图;
图12是本发明实施例中例3中油藏厚度对产油速率影响规律图;
图13是本发明实施例中例3中原油粘度对产油速率影响规律图;
图14为本发明实施例中稠油油藏SAGD产油速度预测装置结构框图一;
图15为本发明实施例中稠油油藏SAGD产油速度预测装置结构框图二。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
图1为本发明实施例中稠油油藏SAGD产油速度预测方法流程图一,如图1所示,该方法包括:
步骤101:SAGD过程中,蒸汽腔存在扇形和侧面泄油区域。根据油藏尺寸(指的是油藏或者室内实验油藏模型)、蒸汽腔的扇形泄油区域最大蒸汽腔夹角、小层间蒸汽腔夹角变化、蒸汽腔的侧面泄油区域小层间夹角变化、油藏厚度和蒸汽腔前缘距离变化值将油藏分成多个小层;
步骤102:当蒸汽腔未发育到生产井处时,计算预热结束后SAGD初始阶段蒸汽腔最大扩散速度;
步骤103:当蒸汽腔处于上升阶段时,计算扇形泄油区域和侧面泄油区域t时刻产油速度,根据扇形泄油区域和侧面泄油区域t时刻产油速度确定蒸汽腔上升阶段t时刻产油速度;
步骤104:当蒸汽腔处于横向发育阶段时,计算蒸汽腔横向发育阶段t时刻产油速度;
步骤105:当蒸汽腔处于下降阶段时,计算蒸汽腔下降阶段t时刻产油速度;
步骤106:分别改变油藏厚度、油相渗透率和原油粘度,重复上述步骤,根据预热结束后SAGD初始阶段蒸汽腔最大扩散速度、蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度和蒸汽腔下降阶段t时刻产油速度确定油藏厚度、油相渗透率和原油粘度对SAGD产油速度的影响规律。
该预测方法可以同时预测SAGD蒸汽腔上升、横向发育和下降三个阶段的产油速度,实现稠油油藏SAGD全过程产油速度预测。与现有方法相比,本发明极大地降低了预测SAGD产油速度所产生的人力和时间成本,提高了预测的准确性。
在本发明实施例中,针对现有技术中存在的“该模型无法计算SAGD过程中含油饱和度变化”这一缺点,如图2所示,该方法还可以包括:
步骤107:根据蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度和蒸汽腔下降阶段t时刻产油速度,计算蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间Δt;
步骤108:当第一个小层的含油饱和度达到残余油饱和度时,根据蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度、蒸汽腔下降阶段t时刻产油速度和蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间Δt,计算每个小层t+Δt时刻的含油饱和度。
具体的,确定油藏厚度、油相渗透率和原油粘度对SAGD产油速度的影响规律还需要重复这两个步骤。
在本发明实施例中,通过以下公式计算预热结束后SAGD初始阶段蒸汽腔最大扩散速度:
其中,Umax,initial是预热结束后SAGD初始阶段蒸汽腔最大扩散速度,m/s;ko为油相渗透率,m2;ρo是原油密度,kg/m3;a为经验常数,值为0.4;m为常数,值在3-5之间;g是重力加速度,m/s2;φ是孔隙度,%;ΔSo是初始含油饱和度,%;Dstart是蒸汽腔的初始直径,m;α是油砂热扩散系数,m2/s;μos是原油粘度,m2/s。
在本发明实施例中,蒸汽腔上升阶段,根据以下公式计算扇形泄油区域t时刻产油速度:
其中,qs为扇形泄油区域t时刻产油速度,m3/s;dqs,i为t时刻第i小层的产油速度,m3/s;L是井的长度,m;h是油藏厚度,m;kro为油相相对渗透率,m2;是流动势函数;Ui是第i小层蒸汽腔移动速度,m/s;Xi是第i小层距蒸汽腔前缘距离,m;μ(Xi)为蒸汽腔前缘Xi处油的粘度,m2/s;θ为蒸汽腔夹角,°;dA是过流面积,m2;n为小层的总数。
在本发明实施例中,蒸汽腔上升阶段,通过以下公式计算侧面泄油区域t时刻产油速度:
其中,qc为侧面泄油区域t时刻产油速度,m3/s;dqc,i为t时刻第i小层的产油速度,m3/s;μ(Xi)为第i小层距离蒸汽腔前缘距离处的原油粘度,m2/s;γ为泄油区域任意一小层的夹角,°;Xi为第i小层距离蒸汽腔前缘的距离。
在本发明实施例中,根据公式(6)和(10)计算蒸汽腔上升阶段t时刻产油速度:
q1=qs+qc (11)
其中,q1为上升阶段t时刻产油速度,m3/s。
在本发明实施例中,横向发育阶段,通过以下公式计算t时刻产油速度:
其中,q2为横向发育阶段t时刻产油速度,m3/s;dq2,i为横向发育阶段t时刻第i小层的产油速度,m3/s。
在本发明实施例中,下降阶段,通过以下公式计算t时刻产油速度:
其中,q3为下降阶段t时刻产油速度,m3/s;dq3,i为下降阶段t时刻第i小层的产油速度,m3/s。
在本发明实施例中,根据质量守恒定律,蒸汽腔从第i个小层移动到第i+1个时,第i小层的可采油已被采完。因此,蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间:
其中,Δt是蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间;So,1是第一个小层的含油饱和度,%;dqj,1是蒸汽腔第一个小层的产油速度,m3/s,So,R是残余油饱和度,%;ΔV为第一小层的体积,m3;j=1,2,3;φ是孔隙度,%。
在本发明实施例中,当第一个小层的含油饱和度达到残余油饱和度时,其他小层的含油饱和度也会相应的降低,每个小层t+Δt时刻含油饱和度计算公式如下:
其中,为第i小层t+Δt时刻含油饱和度;/>为第i小层t时刻含油饱和度;ΔA是每个小层的横截面积,m2。
在本发明实施例中,所述步骤103中的侧面泄油区域t时刻产油速度中的Xi和步骤104中蒸汽腔横向发育阶段t时刻产油速度中的Xi,通过以下公式计算第i小层距离蒸汽腔前缘的距离:
Xi=h(tanγ-tanθ) (18)
在本发明实施例中,所述步骤105中蒸汽腔下降阶段t时刻产油速度中的Xi,通过以下公式计算第i小层距离蒸汽腔前缘的距离:
其中,W为油藏宽度的一半,m。
在本发明实施例中,所述步骤103、所述步骤104和所述步骤105中,蒸汽腔移动到第i个小层时的扩展速度:
其中,ΔXi是与第i小层相邻两个小层之间的距离,m。
在本发明实施例中,所述步骤103和所述步骤105中,第i个小层距离蒸汽腔前缘距离处的原油粘度:
下面举例说明本发明预测方法。
实例1
具体的流程图如图3所示,包括如下步骤:
步骤1.本实例1中,以室内实验油藏模型为基础(长度、高度和厚度尺寸为35×21×3cm),蒸汽腔上升阶段扇形泄油区域最大蒸汽腔夹角θmax为15°,小层间蒸汽腔夹角变化Δθ为1°,侧面泄油区域小层间夹角变化Δγ为1°,油藏厚度h为0.21m,蒸汽腔前缘距离变化值ΔX为0.01m,根据上述参数将整个油藏分为n个小层。蒸汽驱上升阶段见图4所示,横向发育和下降阶段如图5所示。
步骤2.当蒸汽腔未发育到生产井处时,通过公式(1)计算预热结束后SAGD初始阶段蒸汽腔最大扩散速度Umax,initial。
本实例1中,Ko为9.44×10-10m2,ρo为980kg/m3,μos为104×10-6m2/s,Dstart为0.08m,孔隙度φ为0.3,ΔSo为0.95,α为5.87×10-7m2/s,计算所得Umax,initial为4.2787×10-5m/s。
步骤3.蒸汽腔上升阶段,根据公式(2)-(6)计算扇形泄油区域t时刻产油速度qs。
步骤4.蒸汽腔上升阶段,通过公式(7)-(10)计算侧面泄油区域t时刻产油速度qc。
步骤5.根据公式(6)和(10),计算蒸汽腔上升阶段t时刻产油速度q1,即采用公式(11)。
本实例1中,根据公式(2)~(11)计算所得上升阶段产油速度如图6所示。由图6可知,随着产油速度随时间不断增加,这是由于随着蒸汽腔不断向上发育,泄油高度不断增加,重力泄油能力不断增强,蒸汽冷凝成水释放热量,原油不断被加热,所以产油速度逐渐增加,在蒸汽腔到达模型顶部后产油速度达到最大值。此外,由图6可知,计算所得产油速度结果和实验结果基本吻合,相对误差仅为4.1%,计算精度较高,表明该方法可以较为准确的预测SAGD蒸汽腔上升阶段产油速度。
步骤6.横向发育阶段,通过公式(12)-(13)计算t时刻产油速度q2。
本实例1中,根据公式(12)~(13)计算所得横向发育阶段产油速度如图7所示。由图7可知,该阶段产油速度基本保持稳定,这是由于蒸汽腔缓慢的向两侧扩展,该阶段泄油高度已经达到最大值,重力泄油能力保持稳定,所以产油速率基本保持稳定。该计算结果与实验现象相吻合,此外,由图7可知,计算所得产油速度结果更接近与实验结果,相对误差为4.5%,比Bulter模型计算结果更加准确。由此表明,本发明所提方法较Bulter模型预测结果更为准确。
步骤7.下降阶段,通过公式(14)-(15)计算t时刻产油速度q3。
本实例1中,根据公式(14)~(15)计算所得下降阶段产油速度如图8所示。由图8可知,产油速度随着时间逐渐降低,这是由于该阶段蒸汽腔不断下降,重力泄油能力不断降低,模型上部的原油已经大量采出,所以该阶段的产油速率不断下降。此外,由图8可知,计算所得产油速度结果和实验结果基本吻合,相对误差仅为0.99%,满足工程应用的预测精度。
步骤8.根据质量守恒定律,蒸汽腔从第i个小层移动到第i+1个时,第i小层的可采油已被采完。因此,蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间Δt可以由公式(16)计算。
步骤9.当第一个小层的含油饱和度达到残余油饱和度时,其他小层的含油饱和度也会相应的降低,每i小层t+Δt时刻含油饱和度公式(17)如下。
根据本发明优选的,分别改变油藏厚度、油相渗透率和原油粘度,重复步骤1~9,从而确定上述参数对SAGD产油速度的影响规律。
实施例2、
与实施例1区别在于,以室内实验所用油藏模型长度、高度和厚度尺寸为62.5×25×5cm,Ko为9.58×10-10m2,ρo为980kg/m3,μos为104×10-6m2/s,Dstart为0.10m,φ为0.35,ΔSo为0.9,α为5.66×10-7m2/s,计算步骤与实施例1相同。计算所得Umax,initial为2.6731×10- 5m/s。
根据公式(2)~(21)计算得到上升、横向发育和下降阶段产油速度如图9~图11所示。由图9~图11可知,计算结果和实验结果基本吻合,三个阶段计算相对误差分别为4.1%、4.5%和5.0%。由此可知本发明所述预测方法可以准确预测稠油油藏SAGD三个阶段产油速度。
实施例3、
与实施例1区别在于,分别改变油藏高度为15cm(或21cm),原油粘度为实施例1中原油粘度的1.2倍,研究油藏高度和原油粘度对SAGD产油速度的影响规律。计算步骤与实施例1相同。计算所得Umax,initial为4.2580×10-5m/s。
油藏高度对SAGD产油速度的影响如图12所示。由图12可知,随着油藏厚度的增加,蒸汽腔上升阶段向上发育的时间越长,到达盖层的时间也越靠后。此外,随着油藏厚度的增加,横向发育和下降阶段SAGD产油越大,最大产油速度也越大。但是随着油藏厚度的增加,蒸汽腔上升阶段SAGD产油速度变化不大,说明油藏厚度对蒸汽腔上升阶段的SAGD产油速度影响不大,对横向发育和下降阶段的SAGD产油速度影响较大。由此可知,油藏高度越大,SAGD开发效果越好。
原油粘度对SAGD产油速度的影响如图13所示。由图13可知,随着原油粘度的减小,上升、横向发育和下降阶段SAGD产油速度增大,最大产油速度也增大,横向阶段的时间越短,由此可知,对于原油粘度较小的油藏,SAGD开发效果越好。
本发明实施例中还提供了一种稠油油藏SAGD产油速度预测装置,如下面的实施例所述。由于该装置解决问题的原理与稠油油藏SAGD产油速度预测方法相似,因此该装置的实施可以参见稠油油藏SAGD产油速度预测方法的实施,重复之处不再赘述。
图14为本发明实施例中稠油油藏SAGD产油速度预测装置结构框图一,如图14所示,该装置包括:
油藏分层模块01,用于根据油藏尺寸、蒸汽腔的扇形泄油区域最大蒸汽腔夹角、小层间蒸汽腔夹角变化、蒸汽腔的侧面泄油区域小层间夹角变化、油藏厚度和蒸汽腔前缘距离变化值将油藏分成多个小层;
扩散速度计算模块02,用于当蒸汽腔未发育到生产井处时,计算预热结束后SAGD初始阶段蒸汽腔最大扩散速度;
上升阶段产油速度计算模块03,用于当蒸汽腔处于上升阶段时,计算扇形泄油区域和侧面泄油区域t时刻产油速度,根据扇形泄油区域和侧面泄油区域t时刻产油速度确定蒸汽腔上升阶段t时刻产油速度;
横向发育阶段产油速度计算模块04,用于当蒸汽腔处于横向发育阶段时,计算蒸汽腔横向发育阶段t时刻产油速度;
下降阶段产油速度计算模块05,用于当蒸汽腔处于下降阶段时,计算蒸汽腔下降阶段t时刻产油速度;
产油速度的影响规律分析模块06,用于分别改变油藏厚度、油相渗透率和原油粘度,重复上述步骤,根据预热结束后SAGD初始阶段蒸汽腔最大扩散速度、蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度和蒸汽腔下降阶段t时刻产油速度确定油藏厚度、油相渗透率和原油粘度对SAGD产油速度的影响规律。
在本发明实施例中,如图15所示,该装置还可以包括:
时间Δt计算模块07,用于根据蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度和蒸汽腔下降阶段t时刻产油速度,计算蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间Δt;
含油饱和度计算模块08,用于当第一个小层的含油饱和度达到残余油饱和度时,根据蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度、蒸汽腔下降阶段t时刻产油速度和蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间Δt,计算每个小层t+Δt时刻的含油饱和度。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述稠油油藏SAGD产油速度预测方法。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述稠油油藏SAGD产油速度预测方法的计算机程序。
本发明实施例中,根据油藏尺寸、蒸汽腔的扇形泄油区域最大蒸汽腔夹角、小层间蒸汽腔夹角变化、蒸汽腔的侧面泄油区域小层间夹角变化、油藏厚度和蒸汽腔前缘距离变化值将油藏分成多个小层,基于分层理论对稠油油藏SAGD产油速度进行预测,可以同时预测SAGD蒸汽腔上升、横向发育和下降三个阶段的产油速度,实现稠油油藏SAGD全过程产油速度预测。与现有方法相比,本发明极大地降低了预测SAGD产油速度所产生的人力和时间成本,提高了预测的准确性。
另外,本发明还具有如下有益效果:
1、利用分层方法对油藏进行分层能够模拟蒸汽腔的动态发育过程,能够计算蒸汽腔扩展速率和含油饱和度分布。
2、本发明所述预测方法所需基础参数较少,预测过程简单、易懂,可以满足现场技术人员在参数较少情况下快速预测的需求。
3、本发明所述预测方法可用于预测不同原油性质和油藏条件下的SAGD产油速度,满足了现场技术人员全面评估SAGD产油速度的需求。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (30)
1.一种稠油油藏SAGD产油速度预测方法,其特征在于,包括:
根据油藏尺寸、蒸汽腔的扇形泄油区域最大蒸汽腔夹角、小层间蒸汽腔夹角变化、蒸汽腔的侧面泄油区域小层间夹角变化、油藏厚度和蒸汽腔前缘距离变化值将油藏分成多个小层;
当蒸汽腔未发育到生产井处时,计算预热结束后SAGD初始阶段蒸汽腔最大扩散速度;
当蒸汽腔处于上升阶段时,计算扇形泄油区域和侧面泄油区域t时刻产油速度,根据扇形泄油区域和侧面泄油区域t时刻产油速度确定蒸汽腔上升阶段t时刻产油速度;
当蒸汽腔处于横向发育阶段时,计算蒸汽腔横向发育阶段t时刻产油速度;
当蒸汽腔处于下降阶段时,计算蒸汽腔下降阶段t时刻产油速度;
分别改变油藏厚度、油相渗透率和原油粘度,重复上述步骤,根据预热结束后SAGD初始阶段蒸汽腔最大扩散速度、蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度和蒸汽腔下降阶段t时刻产油速度确定油藏厚度、油相渗透率和原油粘度对SAGD产油速度的影响规律。
2.如权利要求1所述的稠油油藏SAGD产油速度预测方法,其特征在于,还包括:
根据蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度和蒸汽腔下降阶段t时刻产油速度,计算蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间Δt;
当第一个小层的含油饱和度达到残余油饱和度时,根据蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度、蒸汽腔下降阶段t时刻产油速度和蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间Δt,计算每个小层t+Δt时刻的含油饱和度。
3.如权利要求1所述的稠油油藏SAGD产油速度预测方法,其特征在于,通过以下公式计算预热结束后SAGD初始阶段蒸汽腔最大扩散速度:
其中,Umax,initial是预热结束后SAGD初始阶段蒸汽腔最大扩散速度,m/s;ko为油相渗透率,m2;ρo是原油密度,kg/m3;a为经验常数,值为0.4;m为常数,值在3-5之间;g是重力加速度,m/s2;φ是孔隙度,%;ΔSo是初始含油饱和度,%;Dstart是蒸汽腔的初始直径,m;α是油砂热扩散系数,m2/s;νos是原油粘度,m2/s。
4.如权利要求3所述的稠油油藏SAGD产油速度预测方法,其特征在于,蒸汽腔上升阶段,根据以下公式计算扇形泄油区域t时刻产油速度:
其中,qs为扇形泄油区域t时刻产油速度,m3/s;dqs,i为t时刻第i小层的产油速度,m3/s;L是井的长度,m;h是油藏厚度,m;kro为油相相对渗透率,m2;是流动势函数;Ui是第i小层蒸汽腔移动速度,m/s;Xi为第i小层距离蒸汽腔前缘的距离,m;ν(Xi)为蒸汽腔前缘Xi处油的粘度,m2/s;θ为蒸汽腔夹角,°;dA是过流面积,m2;n为小层的总数。
5.如权利要求4所述的稠油油藏SAGD产油速度预测方法,其特征在于,蒸汽腔上升阶段,通过以下公式计算侧面泄油区域t时刻产油速度:
其中,qc为侧面泄油区域t时刻产油速度,m3/s;dqc,i为t时刻第i小层的产油速度,m3/s;μ(Xi)为第i小层距离蒸汽腔前缘距离Xi处的原油粘度,m2/s;γ为泄油区域任意一小层的夹角,°;Xi为第i小层距离蒸汽腔前缘的距离。
6.如权利要求5所述的稠油油藏SAGD产油速度预测方法,其特征在于,按照如下公式计算蒸汽腔上升阶段t时刻产油速度:
q1=qs+qc;
其中,q1为上升阶段t时刻产油速度,m3/s。
7.如权利要求6所述的稠油油藏SAGD产油速度预测方法,其特征在于,横向发育阶段,通过以下公式计算t时刻产油速度:
其中,q2为横向发育阶段t时刻产油速度,m3/s;dq2,i为横向发育阶段t时刻第i小层的产油速度,m3/s。
8.如权利要求7所述的稠油油藏SAGD产油速度预测方法,其特征在于,下降阶段,通过以下公式计算t时刻产油速度:
其中,q3为下降阶段t时刻产油速度,m3/s;dq3,i为下降阶段t时刻第i小层的产油速度,m3/s。
9.如权利要求2所述的稠油油藏SAGD产油速度预测方法,其特征在于,蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间:
其中,Δt是蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间;So,1是第一个小层的含油饱和度,%;dqj,1是蒸汽腔第一个小层的产油速度,m3/s,So,R是残余油饱和度,%;ΔV为第一小层的体积,m3;j=1,2,3;φ是孔隙度,%。
10.如权利要求9所述的稠油油藏SAGD产油速度预测方法,其特征在于,每个小层t+Δt时刻含油饱和度计算公式如下:
其中,为第i小层t+Δt时刻含油饱和度;/>为第i小层t时刻含油饱和度;ΔA是每个小层的横截面积,m2。
11.如权利要求5或8所述的稠油油藏SAGD产油速度预测方法,其特征在于,第i小层距离蒸汽腔前缘的距离通过以下公式计算:
Xi=h(tanγ-tanθ);
其中,h是油藏厚度,m;γ为泄油区域任意一小层的夹角,°;θ为蒸汽腔夹角,°。
12.如权利要求7所述的稠油油藏SAGD产油速度预测方法,其特征在于,第i小层距离蒸汽腔前缘的距离通过以下公式计算:
其中,W为油藏宽度的一半,m。
13.如权利要求8所述的稠油油藏SAGD产油速度预测方法,其特征在于,蒸汽腔移动到第i个小层时的扩展速度:
其中,ΔXi是与第i小层相邻两个小层之间的距离,m。
14.如权利要求8所述的稠油油藏SAGD产油速度预测方法,其特征在于,第i个小层距离蒸汽腔前缘距离处的原油粘度:
其中,μos是原油粘度,m2/s;m为常数,值在3-5之间;Ui是第i小层蒸汽腔移动速度,m/s;Xi为第i小层距离蒸汽腔前缘的距离;α是油砂热扩散系数,m2/s。
15.一种稠油油藏SAGD产油速度预测装置,其特征在于,包括:
油藏分层模块,用于根据油藏尺寸、蒸汽腔的扇形泄油区域最大蒸汽腔夹角、小层间蒸汽腔夹角变化、蒸汽腔的侧面泄油区域小层间夹角变化、油藏厚度和蒸汽腔前缘距离变化值将油藏分成多个小层;
扩散速度计算模块,用于当蒸汽腔未发育到生产井处时,计算预热结束后SAGD初始阶段蒸汽腔最大扩散速度;
上升阶段产油速度计算模块,用于当蒸汽腔处于上升阶段时,计算扇形泄油区域和侧面泄油区域t时刻产油速度,根据扇形泄油区域和侧面泄油区域t时刻产油速度确定蒸汽腔上升阶段t时刻产油速度;
横向发育阶段产油速度计算模块,用于当蒸汽腔处于横向发育阶段时,计算蒸汽腔横向发育阶段t时刻产油速度;
下降阶段产油速度计算模块,用于当蒸汽腔处于下降阶段时,计算蒸汽腔下降阶段t时刻产油速度;
产油速度的影响规律分析模块,用于分别改变油藏厚度、油相渗透率和原油粘度,重复上述步骤,根据预热结束后SAGD初始阶段蒸汽腔最大扩散速度、蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度和蒸汽腔下降阶段t时刻产油速度确定油藏厚度、油相渗透率和原油粘度对SAGD产油速度的影响规律。
16.如权利要求15所述的稠油油藏SAGD产油速度预测装置,其特征在于,还包括:
时间Δt计算模块,用于根据蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度和蒸汽腔下降阶段t时刻产油速度,计算蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间Δt;
含油饱和度计算模块,用于当第一个小层的含油饱和度达到残余油饱和度时,根据蒸汽腔上升阶段t时刻产油速度、蒸汽腔横向发育阶段t时刻产油速度、蒸汽腔下降阶段t时刻产油速度和蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间Δt,计算每个小层t+Δt时刻的含油饱和度。
17.如权利要求15所述的稠油油藏SAGD产油速度预测装置,其特征在于,通过以下公式计算预热结束后SAGD初始阶段蒸汽腔最大扩散速度:
其中,Umax,initial是预热结束后SAGD初始阶段蒸汽腔最大扩散速度,m/s;ko为油相渗透率,m2;ρo是原油密度,kg/m3;a为经验常数,值为0.4;m为常数,值在3-5之间;g是重力加速度,m/s2;φ是孔隙度,%;ΔSo是初始含油饱和度,%;Dstart是蒸汽腔的初始直径,m;α是油砂热扩散系数,m2/s;μos是原油粘度,m2/s。
18.如权利要求17所述的稠油油藏SAGD产油速度预测装置,其特征在于,蒸汽腔上升阶段,根据以下公式计算扇形泄油区域t时刻产油速度:
其中,qs为扇形泄油区域t时刻产油速度,m3/s;dqs,i为t时刻第i小层的产油速度,m3/s;L是井的长度,m;h是油藏厚度,m;kro为油相相对渗透率,m2;是流动势函数;Ui是第i小层蒸汽腔移动速度,m/s;Xi是第i小层距蒸汽腔前缘距离,m;μ(Xi)为蒸汽腔前缘Xi处油的粘度,m2/s;θ为蒸汽腔夹角,°;dA是过流面积,m2;n为小层的总数。
19.如权利要求18所述的稠油油藏SAGD产油速度预测装置,其特征在于,蒸汽腔上升阶段,通过以下公式计算侧面泄油区域t时刻产油速度:
其中,qc为侧面泄油区域t时刻产油速度,m3/s;dqc,i为t时刻第i小层的产油速度,m3/s;μ(Xi)为第i小层距离蒸汽腔前缘距离处的原油粘度,m2/s;γ为泄油区域任意一小层的夹角,°;Xi为第i小层距离蒸汽腔前缘的距离。
20.如权利要求19所述的稠油油藏SAGD产油速度预测装置,其特征在于,按照如下公式计算蒸汽腔上升阶段t时刻产油速度:
q1=qs+qc;
其中,q1为上升阶段t时刻产油速度,m3/s。
21.如权利要求20所述的稠油油藏SAGD产油速度预测装置,其特征在于,横向发育阶段,通过以下公式计算t时刻产油速度:
其中,q2为横向发育阶段t时刻产油速度,m3/s;dq2,i为横向发育阶段t时刻第i小层的产油速度,m3/s。
22.如权利要求21所述的稠油油藏SAGD产油速度预测装置,其特征在于,下降阶段,通过以下公式计算t时刻产油速度:
其中,q3为下降阶段t时刻产油速度,m3/s;dq3,i为下降阶段t时刻第i小层的产油速度,m3/s。
23.如权利要求16所述的稠油油藏SAGD产油速度预测装置,其特征在于,蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间:
其中,Δt是蒸汽腔前缘第一个小层含油饱和度到达残余油饱和度的时间;So,1是第一个小层的含油饱和度,%;dqj,1是蒸汽腔第一个小层的产油速度,m3/s,So,R是残余油饱和度,%;ΔV为第一小层的体积,m3;j=1,2,3;φ是孔隙度,%。
24.如权利要求23所述的稠油油藏SAGD产油速度预测装置,其特征在于,每个小层t+Δt时刻含油饱和度计算公式如下:
其中,为第i小层t+Δt时刻含油饱和度;/>为第i小层t时刻含油饱和度;ΔA是每个小层的横截面积,m2。
25.如权利要求19或23所述的稠油油藏SAGD产油速度预测装置,其特征在于,第i小层距离蒸汽腔前缘的距离通过以下公式计算:
Xi=h(tanγ-tanθ);
其中,h是油藏厚度,m;γ为泄油区域任意一小层的夹角,°;θ为蒸汽腔夹角,°。
26.如权利要求22所述的稠油油藏SAGD产油速度预测装置,其特征在于,第i小层距离蒸汽腔前缘的距离通过以下公式计算:
其中,W为油藏宽度的一半,m。
27.如权利要求23所述的稠油油藏SAGD产油速度预测装置,其特征在于,蒸汽腔移动到第i个小层时的扩展速度:
其中,ΔXi是与第i小层相邻两个小层之间的距离,m。
28.如权利要求23所述的稠油油藏SAGD产油速度预测装置,其特征在于,第i个小层距离蒸汽腔前缘距离处的原油粘度:
其中,μos是原油粘度,m2/s;m为常数,值在3-5之间;Ui是第i小层蒸汽腔移动速度,m/s;Xi为第i小层距离蒸汽腔前缘的距离;α是油砂热扩散系数,m2/s。
29.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至14任一所述方法。
30.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至14任一所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011433241.7A CN114622882B (zh) | 2020-12-10 | 2020-12-10 | 稠油油藏sagd产油速度预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011433241.7A CN114622882B (zh) | 2020-12-10 | 2020-12-10 | 稠油油藏sagd产油速度预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114622882A CN114622882A (zh) | 2022-06-14 |
CN114622882B true CN114622882B (zh) | 2024-03-26 |
Family
ID=81896502
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011433241.7A Active CN114622882B (zh) | 2020-12-10 | 2020-12-10 | 稠油油藏sagd产油速度预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114622882B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102758603A (zh) * | 2012-07-10 | 2012-10-31 | 中国石油天然气股份有限公司 | 一种超稠油油藏sagd开采后期注空气开采方法 |
CA2898897A1 (en) * | 2015-07-29 | 2017-01-29 | Alberta Innovates - Technology Futures | Partial height steam chamber sagd |
CN106593368A (zh) * | 2016-12-07 | 2017-04-26 | 中国石油天然气股份有限公司 | 一种改善sagd开发效果的预处理方法 |
CN106951649A (zh) * | 2017-03-27 | 2017-07-14 | 中国石油大学(华东) | 一种测定水平井sagd蒸汽腔扩展速度的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017131850A1 (en) * | 2015-12-01 | 2017-08-03 | Conocophillips Company | Single well cross steam and gravity drainage (sw-xsagd) |
-
2020
- 2020-12-10 CN CN202011433241.7A patent/CN114622882B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102758603A (zh) * | 2012-07-10 | 2012-10-31 | 中国石油天然气股份有限公司 | 一种超稠油油藏sagd开采后期注空气开采方法 |
CA2898897A1 (en) * | 2015-07-29 | 2017-01-29 | Alberta Innovates - Technology Futures | Partial height steam chamber sagd |
CN106593368A (zh) * | 2016-12-07 | 2017-04-26 | 中国石油天然气股份有限公司 | 一种改善sagd开发效果的预处理方法 |
CN106951649A (zh) * | 2017-03-27 | 2017-07-14 | 中国石油大学(华东) | 一种测定水平井sagd蒸汽腔扩展速度的方法 |
Non-Patent Citations (1)
Title |
---|
提高直井与水平井组合SAGD泄油速率技术研究;李巍;刘永建;;石油钻探技术(第02期);第87-91页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114622882A (zh) | 2022-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110334431B (zh) | 一种低渗透致密气藏单井控制储量计算及剩余气分析方法 | |
CN106777651B (zh) | 注采平衡原理的油水井产量劈分方法 | |
CN107506566A (zh) | 一种新型泥石流动力学数值模拟分析方法及系统 | |
He et al. | Pore-scale experimental investigation of the effect of supercritical CO2 injection rate and surface wettability on salt precipitation | |
Rasmusson et al. | Residual trapping of carbon dioxide during geological storage—Insight gained through a pore-network modeling approach | |
Jackson et al. | Dynamic-wetting effects in finite-mobility-ratio Hele-Shaw flow | |
WO2012109191A1 (en) | A quantitative method of determining safe steam injection pressure for enhanced oil recovery operations | |
CN114622882B (zh) | 稠油油藏sagd产油速度预测方法 | |
Zhu et al. | Numerical simulation of displacement characteristics of CO2 injected in pore-scale porous media | |
Moghanloo et al. | Applying fractional flow theory to evaluate CO2 storage capacity of an aquifer | |
CN106056460B (zh) | 一种确定化学驱提高采收率贡献的计算方法 | |
CN104991992B (zh) | 孔板水流压力恢复长度计算方法 | |
CN111155974B (zh) | 缝洞型油藏注水替油的油水置换状态预测方法 | |
Bagudu et al. | Pore-to-core-scale network modelling of CO 2 migration in porous media | |
Guo et al. | Vertically integrated dual-porosity and dual-permeability models for CO2 sequestration in fractured geological formation | |
CN104392131A (zh) | 一种水驱沙过程中破碎岩石渗流场计算方法 | |
Sherafati et al. | Dynamic relative permeability and Simulation of WAG injection processes | |
CN112112639B (zh) | 一种凝析气藏循环注气条件下地层压力确定方法及系统 | |
Zhou et al. | Influence of wettability in immiscible displacements with lattice Boltzmann method | |
CN111502639A (zh) | 一种救援井最小压井排量的确定方法 | |
CN116205162B (zh) | 一种高温高压气藏结垢动态沉积堵塞预测模型构建方法 | |
CN114592838B (zh) | 一种层状超稠油油藏蒸汽驱潜力评价方法及系统 | |
RU2669647C1 (ru) | Способ разработки залежи высоковязкой и сверхвязкой нефти тепловыми методами на поздней стадии разработки | |
CN106529199A (zh) | 一种砾岩油藏化学驱井距的确定方法 | |
Ghazvini et al. | A new mathematical model for force gravity drainage in fractured porous media |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |