CN106951649B - 一种测定水平井sagd蒸汽腔扩展速度的方法 - Google Patents

一种测定水平井sagd蒸汽腔扩展速度的方法 Download PDF

Info

Publication number
CN106951649B
CN106951649B CN201710188041.1A CN201710188041A CN106951649B CN 106951649 B CN106951649 B CN 106951649B CN 201710188041 A CN201710188041 A CN 201710188041A CN 106951649 B CN106951649 B CN 106951649B
Authority
CN
China
Prior art keywords
vapor chamber
formula
depth
expansion rate
well
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
Application number
CN201710188041.1A
Other languages
English (en)
Other versions
CN106951649A (zh
Inventor
李兆敏
武守亚
王宏远
鹿腾
杨建平
魏耀
王诗中
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201710188041.1A priority Critical patent/CN106951649B/zh
Publication of CN106951649A publication Critical patent/CN106951649A/zh
Application granted granted Critical
Publication of CN106951649B publication Critical patent/CN106951649B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • 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
    • E21B43/2406—Steam assisted gravity drainage [SAGD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mining & Mineral Resources (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Geology (AREA)
  • Geochemistry & Mineralogy (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Fluid Mechanics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本发明涉及一种测定水平井SAGD蒸汽腔扩展速度的方法，基于一定的假设条件、蒸汽腔扩展理论以及传热导热微分方程，推导出蒸汽腔外缘温度分布函数以及蒸汽腔前缘扩展速度预测模型。根据数学模型并结合现场检测资料，测定水平井组不同时刻蒸汽腔边缘倾斜角以及扩展速度。本发明可直接采用监测井监测温度数据测定蒸汽腔扩展速度，与现有预测蒸汽腔扩展速度方法比较，节省现有方法繁琐过程和冗长时间，节约人力物力财力，快捷准确的获得扩展速度值。

Description

一种测定水平井SAGD蒸汽腔扩展速度的方法
技术领域
本发明涉及一种测定水平井SAGD蒸汽腔扩展速度的方法，属于油气田开发工程技术领域。
背景技术
随着原油开采技术的日渐成熟，稠油和超稠油油藏的开发技术逐渐受到人们的重视。相比于常规油藏，稠油油藏中的原油粘度较大，流动性较差，常规采油技术难以达到较好的采油效果。我国拥有丰富的稠油和超稠油资源，主要的几大油区有：克拉玛依油区、辽河油区、胜利油区、塔里木油区等。目前蒸汽辅助重力泄油(SAGD)技术是比较前沿的热力采油技术，最早由Butler博士在1978年通过分析研究盐井注水采盐的原理，通过类比方法应用到稠油热采技术过程中，提出SAGD技术基本概念。
SAGD技术在加拿大经过20多年的推广已经得到广泛的应用，已是一套工业化成熟的技术。利用重力原理，SAGD技术的采收率最高可以达到70％以上，是常规蒸汽吞吐的2～3倍。SAGD是将流体热对流与热传导相结合，以蒸汽作为加热介质，依靠原油的重力作用进行开发的稠油热采技术，因此SAGD过程中蒸汽腔的扩展发育对SAGD开发效果具有重要影响。同时在SAGD开发方式中，蒸汽腔的形成特点及变化扩展趋势对生产效果具有极其重要的影响。
SAGD技术应用开发过程如下：超稠油经过油层预热形成热连通后，上部注汽井注入蒸汽向上超覆在地层中形成蒸汽腔，蒸汽腔向上及侧向扩展，注入蒸汽与原油发生热交换，加热后的原油和蒸汽冷凝水在重力作用下泄流至下部生产井产出，随着原油的不断采出，蒸汽腔将逐渐扩大。蒸汽腔的扩展可划分为泄油初期(蒸汽腔上升阶段)、泄油高峰期(蒸汽腔沿油层顶部向外扩展阶段)、泄油末期(蒸汽腔扩展到油藏侧边界或井组控制边界时，蒸汽腔沿边界下降阶段)3个阶段，蒸汽腔扩展过程如图1所示。
蒸汽腔波及规律研究是SAGD跟踪调查的一项必不可少的研究工作，而蒸汽腔扩展速度是研究蒸汽驱波及规律的重要参数，其表征SAGD技术在油藏条件下的实施效果。预测求解现场实践中SAGD蒸汽腔扩展速度大小对SAGD工艺具有重要参考价值，因此对现场指导有重要意义。目前对于蒸汽腔扩展速度的求解主要是采用油藏数值模拟方法：由现场工艺及生产数据，建立油藏数值模拟模型并进行历史拟合，获得实际油藏近似模型，从而模拟油藏SAGD开发，预测蒸汽腔扩展速度。此方法计算量较大，不适用于快速预测求解蒸汽腔扩展速度。
发明内容
针对现有计算蒸汽腔扩展速度方法的不足，本发明提供一种测定水平井SAGD蒸汽腔扩展速度的方法。本发明基于一定的假设条件、蒸汽腔扩展理论以及传热导热微分方程，推导出蒸汽腔外缘温度分布函数以及蒸汽腔前缘扩展速度预测模型。根据数学模型并结合现场监测资料，计算水平井组不同时刻蒸汽腔边缘倾斜角以及扩展速度。
本发明的技术方案如下：
一种测定水平井SAGD蒸汽腔扩展速度的方法，包括步骤如下：
S1、蒸汽腔简化模型
在蒸汽腔沿油层顶部向外扩展阶段，为便于研究蒸汽腔扩展速度，做以下假设：
1)蒸汽腔已达油藏顶层，向两侧扩展；
2)只考虑蒸汽腔前缘法线方向的一维驱替过程；
3)传热方向仅为垂直于蒸汽腔外缘的方向，即一维传热过程；
4)蒸汽与冷油之间以平板前缘和恒定速度；
5)只考虑垂直于蒸汽腔外缘方向的传导，不考虑对流；
6)在一定时刻系统处于准稳态过程，即蒸汽腔沿边缘法线方向以固定速度推进；
7)忽略稠油流动过程中的热损失；
8)油层热容不随温度变化；
9)导热系数随温度线性递减。
因此，蒸汽腔扩展过程可简化为如图2所示形式。
S2、蒸汽腔扩展简化模型
S2.1温度分布函数
基本传热微分方程为：
式(1)中，K表示油藏导热系数；T表示油藏温度；ρc表示冷凝液密度；cpc表示冷凝液比热容；ρr表示砂岩密度；cpr表示砂岩比热容；表示x、y、z方向上达西速率；表示净热流通量。
引入变量表观距离ξ：
式(2)中，Ux表示蒸汽腔边缘沿x方向上的运移速度；基于S1小节的准稳态等假设条件，并将式(2)带入式(1)，可得：
设定边界条件为：
式(4)中，Tst表示蒸汽注入温度；Tr表示油藏原始温度。
并设定导热系数是关于温度的线性函数：
K＝A+BT (5)
式(5)中A、B表示油藏导热系数函数系数。
结合式(1)～(4)可推导得蒸汽腔边缘温度分布函数：
式(6)中x表示油藏位置距蒸汽腔边缘距离；t表示蒸汽腔边缘推进时间。因此式(6)可有效描述蒸汽腔一维扩展过程中某时刻某位置处的温度。
S2.2扩展速度模型
由式(6)可知距蒸汽腔边缘距离x处，在已知蒸汽注入温度Tst、油藏原始温度Tr的条件下，可推导出蒸汽腔扩展速度。在同一监测井下，可分别测得三处不同深度下的温度，如图3所示。
式(6)可简化表示为：
式(7)中：
c1，c2表示扩展速度函数系数。
在同一时刻测得的两处深度、温度：
其中x可由深度表示：
x0、x1、x2分别表示同一时刻蒸汽腔边缘外同一监测位置下三个监测点x方向坐标，h0、h1、h2表示所对应深度。由式(9)(10)化简可得：
式(11)可写为：
因此，蒸汽腔水平扩展速度U为：
根据式(13)可在测得监测深度h及对应监测温度T的条件下求得此刻蒸汽腔扩展速度。
当油藏导热系数K视为常数时，蒸汽腔水平扩展速度U为：
S3、扩展速度求解
S3.1饱和岩石导热系数求解
由岩石、原油、水等物质的导热系数，经过理论公式求得：
油藏导热系数K可由下式求得：
式中：λd—干燥砂岩的导热系数，W/(m·℃)，取值0.877W/(m·℃)；
λw—水的导热系数，W/(m·℃)，取值0.675W/(m·℃)；
λa—空气的导热系数，W/(m·℃)，取值0.040W/(m·℃)；
ρl—饱和砂岩的密度，g/cm3；
ρd—干燥砂岩的密度，g/cm3，取值2.08g/cm3；
m—胶结系数，无因次，取值1.8。
热容量M是指单位体积的油藏岩石，温度升高1℃所需的热量，单位为kJ/(m3·℃)，与密度ρ、比热容C的关系为：
M＝ρgC (16)
油藏热容量可由下式求得：
式(17)中Bi表示体积系数，下标r,o,w,g分别表示砂岩、原油、水与气体；S表示饱和度。
S3.2倾斜角求解
由公式(7)和(13)可知，在注入井深度和蒸汽腔顶层深度范围内，深度差h-h1与温度函数f(T)-f(T1)成正比。
测定监测温度数据(T、T1)、可求得深度差与温度函数的回归参数c，从而可求得当温度T＝Tst时的假想深度h0。
根据几何关系求得蒸汽腔外缘倾斜角度θ：
式中H2表示水平生产井深度；S表示水平生产井与监测井水平距离。
S3.3蒸汽腔扩展速度求解
由蒸汽腔速度扩展理论公式(13)，结合所求物性参数及蒸汽腔边缘倾斜角，可获得蒸汽腔水平扩展速度。
S3.4蒸汽腔扩展位置求解
由生产井位置及蒸汽腔边缘倾斜角，可获得蒸汽腔在盖层扩展位置：
s＝(hcap-h0)/tanθ (20)
式(20)中s表示蒸汽腔扩展边缘在盖层处的位置与生产井水平距离；hcap表示盖层深度。由生产井深度即可求得蒸汽腔扩展位置。
本发明的有益效果：
本发明采用上述方法可直接采用监测井监测温度数据测定蒸汽腔扩展速度，与现有预测蒸汽腔扩展速度方法比较，节省现有方法繁琐过程和冗长时间，节约人力物力财力，快捷准确的获得扩展速度值。
附图说明
图1为本发明蒸汽腔扩展示意图。
图2为本发明蒸汽腔扩展原理图。
图3为本发明蒸汽腔温度监测图。
具体实施方式
下面通过具体实施例对本发明作进一步说明，但不限于此。
实施例1
以某井区某水平井组注蒸汽开采工程为例，通过对应观察井可检测到距水平注入井13.9米处的垂向温度分布。设定蒸汽腔温度为250℃，油藏原始温度为15℃，生产井深度350米。
根据观察井资料，选取下表数据(三个时间节点，每个时间节点六点监测点数据)：
表1观察井深度-温度资料
(1)物性参数
水平井组所对应油层孔隙度为29.7％，渗透率为1078.1mD，含油饱和度为70％。原油与岩石物性如表2所示。
表2观察井深度-温度资料
将数据带入式(15)可求得油藏导热系数为2.1808W/(m·℃)。
将数据带入式(17)可求得油藏热容量为1.8574×106J/m3·℃。
(2)倾斜角
导热系数的回归参数A比B的数量级高出很多，因此式(13)中简化系数c1近似等于0。因此由式(7)(13)可知，在注入井深度和蒸汽腔顶层深度范围内，深度差h-h1与温度函数f(T)-f(T1)成正比。
由表1所列时间1监测数据可求得回归参数c＝5.4903。从而可求得当温度T＝Tst时的深度h0＝320.4809m。
根据式(19)可求得蒸汽腔外缘倾斜角度，即θ＝64.7851°(tanθ＝2.1237)。同理可求得其他两个监测时间点处的倾斜角，如表3所示。
表3倾斜角
时间1 时间2 时间3
θ 64.7851 64.1849 64.4188
(3)蒸汽腔扩展速度
由蒸汽腔速度扩展理论公式，结合所求物性参数及蒸汽腔边缘倾斜角，可获得蒸汽腔水平扩展速度。蒸汽腔水平扩展速度理论公式为：
监测井三个监测时间点所对应的扩展速度如表4所示。
表4蒸汽腔水平扩展速度
时间1 时间2 时间3
U m/月 0.9548 1.2618 0.7911
结合工程实践以及国内外参考文献可知，蒸汽腔扩展速度范围为0.7～1.5m/月，对比计算结果表明该方法对预测蒸汽腔扩展速度有一定的准确性。
以上实施方式仅用于说明本发明，而并非对本发明的限制，有关技术领域的普通技术人员，在不脱离本发明的精神和范围的情况下，还可以做出各种变化和变型，因此所有等同的技术方案也属于本发明的范畴，本发明的专利保护范围应由权利要求限定。

Claims (1)

1.一种测定水平井SAGD蒸汽腔扩展速度的方法，包括步骤如下：
1)求解导热系数K
式中：λd—干燥砂岩的导热系数，W/(m·℃)，取值0.877W/(m·℃)；
λw—水的导热系数，W/(m·℃)，取值0.675W/(m·℃)；
λa—空气的导热系数，W/(m·℃)，取值0.040W/(m·℃)；
ρl—饱和砂岩的密度，g/cm3；
ρd—干燥砂岩的密度，g/cm3，取值2.08g/cm3；
m—胶结系数，无因次，取值1.8；
2)求解油藏热容量M
热容量M与密度ρ、比热容C的关系为：
M＝ρgC (16)
油藏热容量计算式为：
式(17)中B表示体积系数，下标r,o,w,g分别表示砂岩、原油、水与气体；S表示饱和度；
3)求解倾斜角θ
在注入井深度和蒸汽腔顶层深度范围内，深度差h-h1与温度函数f(T)-f(T1)成正比，如式(18)所示；
c表示深度差与温度函数的回归参数，θ蒸汽腔外缘倾斜角度；
测定监测井监测温度数据T、T1，求得深度差与温度函数的回归参数c，从而求得当温度T＝Tst时的深度h0；
根据几何关系求得蒸汽腔外缘倾斜角度θ：
式中H2表示水平生产井深度；X表示水平生产井与监测井水平距离；
4)求解蒸汽腔水平扩展速度U
测定监测井监测温度数据T1、T2和深度数据h1、h2，由蒸汽腔扩展速度公式(22)，结合公式(16)和物性参数M、K及蒸汽腔边缘倾斜角θ，求得蒸汽腔水平扩展速度U；
CN201710188041.1A 2017-03-27 2017-03-27 一种测定水平井sagd蒸汽腔扩展速度的方法 Active CN106951649B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710188041.1A CN106951649B (zh) 2017-03-27 2017-03-27 一种测定水平井sagd蒸汽腔扩展速度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710188041.1A CN106951649B (zh) 2017-03-27 2017-03-27 一种测定水平井sagd蒸汽腔扩展速度的方法

Publications (2)

Publication Number Publication Date
CN106951649A CN106951649A (zh) 2017-07-14
CN106951649B true CN106951649B (zh) 2018-06-29

Family

ID=59473171

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710188041.1A Active CN106951649B (zh) 2017-03-27 2017-03-27 一种测定水平井sagd蒸汽腔扩展速度的方法

Country Status (1)

Country Link
CN (1) CN106951649B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109594969A (zh) * 2017-09-28 2019-04-09 中国石油天然气股份有限公司 蒸汽腔的解析方法
CN108920795B (zh) * 2018-06-22 2021-08-27 中国石油天然气股份有限公司 热采生产井井底蒸汽腔预测方法及装置
CN111119820B (zh) * 2018-10-30 2022-08-05 中国石油天然气股份有限公司 Sagd采油方法
CN114622882B (zh) * 2020-12-10 2024-03-26 中国石油天然气股份有限公司 稠油油藏sagd产油速度预测方法
CN112943194B (zh) * 2021-03-03 2023-01-06 中国石油天然气股份有限公司 一种预防sagd开发过程中边水下内侵的方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102272418A (zh) * 2008-11-28 2011-12-07 普拉德研究及开发股份有限公司 用于估计sagd过程特性的方法
CN104389568A (zh) * 2014-09-29 2015-03-04 中国石油大学(北京) 蒸汽辅助重力泄油过程中气体辅助用量的获取方法及装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160061014A1 (en) * 2014-08-28 2016-03-03 Cenovus Energy Inc. Hydraulically unitary well system and recovery process (huwsrp)

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102272418A (zh) * 2008-11-28 2011-12-07 普拉德研究及开发股份有限公司 用于估计sagd过程特性的方法
CN104389568A (zh) * 2014-09-29 2015-03-04 中国石油大学(北京) 蒸汽辅助重力泄油过程中气体辅助用量的获取方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
SAGD Process Monitoring in Heavy Oil Reservoir Using UWB Radar Techniques;Daniel Oloumi等;《IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES》;20160630;第64卷(第6期);第1884-1895页 *

Also Published As

Publication number Publication date
CN106951649A (zh) 2017-07-14

Similar Documents

Publication Publication Date Title
CN106951649B (zh) 一种测定水平井sagd蒸汽腔扩展速度的方法
CN110206522B (zh) 一种页岩气藏压裂水平井压裂液返排模拟方法
CN107044277B (zh) 低渗透非均质油藏水平井重复压裂增产潜力评价方法
CN103161435B (zh) 一种稠油热采直井试井解释方法
CN105089595B (zh) 水平压裂裂缝导流作用下的油藏数值模拟方法及装置
CN104389568B (zh) 蒸汽辅助重力泄油过程中气体辅助用量的获取方法及装置
CN104594872A (zh) 一种优化致密气藏压裂水平井裂缝导流能力的方法
CN102272418A (zh) 用于估计sagd过程特性的方法
CN109522634A (zh) 一种致密气多段体积压裂水平井数值分析方法
CN103161436B (zh) 一种稠油热采水平井试井解释方法
CN104453834A (zh) 一种井组注采关系优化调整方法
CN106469333B (zh) 一种低渗稠油油藏水平井热采压力分布预测方法
CN114692472B (zh) 储层改造辅助降压开发天然气水合物藏的数值模拟方法
CN108868731A (zh) 一种裂缝性储层酸压动态综合滤失系数的计算方法
Zhang et al. Efficient flow rate profiling for multiphase flow in horizontal wells using downhole temperature measurement
CN111734394B (zh) 一种确定致密油藏压裂井不定常流井底压力的方法
CN108959767A (zh) 一种窄河道型气藏不同井型凝析油伤害数值模拟方法
CN103455667A (zh) 充气法治理承压含水层海水入侵的数值模拟方法
Chuanzhi et al. Identification and quantitative description of large pore path in unconsolidated sandstone reservoir during the ultra-high water-cut stage
Gao et al. Numerical simulation study of a novel horizontally layered enhanced geothermal system: A case study of the Qiabuqia geothermal area, Qinghai Province, China
Guo et al. Oil production rate predictions for steam assisted gravity drainage based on high-pressure experiments
Zhao et al. A field test data research based on a new hydraulic parameters quick test technology
CN104809283A (zh) 一种稠油油藏内的夹层识别方法
CN112989721B (zh) 一种致密油藏体积压裂水平井改造体积的快速计算方法
Hunt et al. Reinjection and gravity changes at Rotokawa geothermal field, New Zealand

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