CN114282387A - 一种基于dts的稠油油藏注蒸汽水平井综合评价方法 - Google Patents
一种基于dts的稠油油藏注蒸汽水平井综合评价方法 Download PDFInfo
- Publication number
- CN114282387A CN114282387A CN202111642620.1A CN202111642620A CN114282387A CN 114282387 A CN114282387 A CN 114282387A CN 202111642620 A CN202111642620 A CN 202111642620A CN 114282387 A CN114282387 A CN 114282387A
- Authority
- CN
- China
- Prior art keywords
- iteration
- inversion
- temperature
- matrix
- heavy oil
- 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.)
- Pending
Links
Images
Abstract
本发明公开了一种基于DTS的稠油油藏注蒸汽水平井综合评价方法,包括以下步骤:根据储层及井筒参数,设置初始赋值反演目标参数向量及阻尼因子,设置初始迭代步数;计算第k步反演迭代的温度剖面反演误差;计算第k步迭代的温度雅克比矩阵及对应的温度对角矩阵;根据温度对角矩阵计算反演目标参数增量;计算第k+1步的温度剖面和对应的k+1步的反演误差;将k步与k+1步中的反演误差进行比较;若满足迭代终止条件则完成迭代并输出第k+1步的反演参数向量作为稠油油藏注蒸汽水平井反演解释结果,若不满足则从第k步开始重复迭代,直至满足终止条件,本发明提供了一种对稠油油藏生产和量化评价的新思路和新手段,同时也能显著的提升稠油热采的经济效益。
Description
技术领域
本发明涉及一种基于DTS的稠油油藏注蒸汽水平井综合评价方法,属于油气藏开发技术领域。
背景技术
目前稠油油藏的开采有热开采和冷开采两种手段,其中稠油热采为主要的开采技术手段包括:蒸汽吞吐、蒸汽驱、蒸汽辅助重力泄油、电加热采油和超声波采油等。稠油热采技术种类繁多,在油田工程应用中蒸汽吞吐工艺是应用最成熟和使用最广泛的热采工艺,该工艺支撑了一半以上的稠油热采产量。
蒸汽吞吐过程包含了渗流力学、化学和热动力学等多种科学原理,是一个综合性很强的复杂过程,所以稠油热采工艺投资较高且开发风险也比普通油田大。而如何去认识稠油水平井的注汽剖面情况并解决水平井注汽剖面不均匀的问题,设定该油井最适宜的注汽量从而提升综合经济效益,成为了稠油开采的核心问题。
随着分布式光纤测温(Distributed Temperature Sensing)技术在石油领域的应用逐渐加深,采用DTS技术可以实现全井段实时监测,能提供准确的连续温度剖面数据。通过大量文献调研表明流体在井筒中引起的温度变化能够反映出流体的流动状态,完全可以将DTS技术运用于稠油油藏注蒸汽水平井综合评价的定量解释。
因此,通过建立稠油油藏注蒸汽水平井吸汽剖面的反演模型,来对流体流入量与温度变化对应的关系进行量化评价,完成通过温度信息来定量解释吸汽剖面,为水平井注汽效果评价及生产优化提供实际依据。该方法为稠油油藏注蒸汽水平井吸汽剖面的确定提供了一种全新的技术手段,同时也能显著地提升稠油热采的经济效益。
发明内容
本发明的目的在于克服现有技术的缺点,提供一种基于DTS的稠油油藏注蒸汽水平井综合评价方法,包括以下步骤:根据储层及井筒参数,设置初始赋值反演目标参数向量及阻尼因子,设置初始迭代步数;计算第k步反演迭代的温度剖面反演误差;计算第k步迭代的温度雅克比矩阵及对应的温度对角矩阵;根据温度对角矩阵计算反演目标参数增量;计算第k+1 步的温度剖面和对应的k+1步的反演误差;将k步与k+1步中的反演误差进行比较;若满足迭代终止条件则完成迭代并输出第k+1步的反演参数向量作为稠油油藏注蒸汽水平井反演解释结果,若不满足则从第k步开始重复迭代,直至满足终止条件,本发明提供了一种对稠油油藏生产和量化评价的新思路和新手段,同时也能显著的提升稠油热采的经济效益。
为实现以上技术效果,采用如下技术方案:
一种基于DTS的稠油油藏注蒸汽水平井综合评价方法,包括以下步骤:
进一步的,所述步骤S3和步骤S7中计算反演误差方法为:
式中,为反演目标参数向量,Ti为第i点处的井筒温度剖面计算值,N为反演目标参数的数量,δ表示微小波动变量,为第k步中反演目标参数在有N个目标参数时的值,为偏导数,反映其函数沿坐标轴的变化率,ej为单位向量中的第j分量,下标i=1,2…N,下标j=1,2…N。
进一步的,所述步骤S4中计算对应的温度对角矩阵Ωk的方法为:
式中,
进一步的,所述步骤S9中迭代终止条件为:
进一步的,所述反演目标参数为吸汽剖面、渗透率分布、蒸汽干度分布、加热半径和蒸汽波及范围。
本发明的有益效果为:
本发明提出一种基于DTS的稠油油藏注蒸汽水平井综合评价方法,具有如下有益效果:
1)本发明通过对稠油油藏注蒸汽水平井井筒温度数据进行反演,能够完成吸汽剖面、渗透率分布、蒸汽干度分布、加热半径、蒸汽波及范围和储层动用范围的定量解释;
2)克服了利用常规手段难以找到稠油油藏注汽剖面的问题,为稠油油藏注蒸汽吸汽剖面的确定提供了一种全新的技术手段,同时也能显著地提升稠油热采的经济效益;
3)本发明还能在渗透率分布、蒸汽干度分布、加热半径和蒸汽波及范围等进行反演解释,为本领域专业人员提供了对稠油油藏生产和量化评价的新思路。
本发明提供了一种对稠油油藏生产和量化评价的新思路和新手段,同时也能显著的提升稠油热采的经济效益。
附图说明
图1为本发明实施例中稠油油藏注蒸汽水平井吸汽剖面反演解释流程示意图;
图2为本发明实施例中稠油油藏水平井实测井筒温度剖面示意图;
图3为本发明实施例中反演模拟的井筒温度剖面与实测井筒温度剖面拟合示意图;
图4为本发明实施例中反演解释出的吸汽剖面分布示意图;
图5为本发明实施例中反演解释出的渗透率分布;
图6为本发明实施例中反演解释出的蒸汽干度分布;
图7为本发明实施例中反演解释出的加热半径;
图8为本发明实施例中反演解释出的储层动用范围。
具体实施方式
下面结合附图对本发明做进一步的描述,本发明的保护范围不局限于以下所述:
显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域技术人员在没有付出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1:
实例井X为位于渤海南部海域的水平井,该井采用套管固井完井,完钻井深2160m,表层套管下入深度384.35m,技术套管下入深度1850.21m,采用旋转导向造斜,造斜深度120m,最大井斜90.17°,最大位移1241.73m,最大全角变化率4.74°/30m,X井深2160m,裸眼水平段长311m,裸眼段采用筛管防砂完井,光纤长度2220.7m,井口预留80m,光纤入井最大深度2131.7m。
如图1所示,本发明的一种基于DTS的稠油油藏注蒸汽水平井综合评价方法,以稠油油藏作为目标水平井、以吸汽剖面作为反演目标参数为例阐述采用所述的方法进行稠油油藏注蒸汽水平井吸汽剖面反演解释的具体步骤,其他渗透率分布、蒸汽干度分布、加热半径和蒸汽波及范围等反演解释与实施例1中类似,不再赘述。
式中,为反演目标参数向量,Ti为第i点处的井筒温度剖面计算值,N为反演目标参数的数量,δ表示微小波动变量,为第k步中反演目标参数在有N个目标参数时的值,为偏导数,反映其函数沿坐标轴的变化率,ej为单位向量中的第j分量,下标i=1,2…N,下标j=1,2…N。
(9)判断是否满足如下迭代终止条件,若满足则完成迭代并输出结果作为稠油油藏注蒸汽水平井吸汽剖面反演解释结果,如图3,若不满足则从步骤(2)开始进行重复迭代,直至满足终止条件,输出产出剖面反演解释结果如图4所示。
其中反演迭代终止条件满足以下其一即可:
本发明所述的温度剖面模型如下:
注蒸汽垂直井段非等温流动模型:
注蒸汽水平井井段非等温流动模型:
注蒸汽垂直井段热损失模型:
注蒸汽垂水平井井段热损失模型:
垂直井段温度剖面预测模型:
水平井井段温度剖面预测模型:
式中:P为微元段压力,Pa;λ为流体在垂直井筒中流动时的阻力系数,无因次;D为垂直井筒的直径,m;g为重力加速度,9.8m/s2;ρm为注入流体的密度,kg/m3;vm为注入流体在井筒中的流动速度,m/s;vsg注入气体在井筒中的流动速度,m/s;is为蒸汽在井筒中的质量流量,kg/s;τc为套管内壁与蒸汽之间的摩擦力,N;Ah为微元段的横截面积,m2;T为微元段的温度,℃;R为传热过程中的总热阻,(m·K)/W;Ts为井筒中的蒸汽温度,℃; Th为水泥环外缘处的温度,℃;Q为微元段内的热损失,W;Te为原始地层温度,℃;M为井筒流体质量,kg;h1为饱和湿蒸汽中液体的焓,kJ/kg;Cp为饱和水蒸汽的热容,kJ/(kg·℃); dl为微元段的长度,m;dW为单位时间内摩擦力所做的功,W;vl为蒸汽在微元段内沿水平方向的流动速度,m/s;vr为微元段内的蒸汽沿径向流入油藏的速度,m/s;下标x表示在x 方向,下标y表示在y方向,下标z表示在z方向。
以上所述为本发明对稠油油藏注蒸汽中的吸汽剖面进行反演的步骤及方法。
但本发明可以但不限于对稠油油藏注蒸汽中的吸汽剖面,如图4,渗透率分布,如图5,蒸汽干度分布,如图6,加热半径,如图7,和储层动用范围,如图8,进行反演。其原理与步骤方法与实施例1中的吸汽剖面反演步骤和原理一致,在此不再赘述。
本发明公开了一种基于DTS的稠油油藏注蒸汽水平井综合评价方法,包括以下步骤:根据储层及井筒参数,设置初始赋值反演目标参数向量及阻尼因子,设置初始迭代步数;计算第k步反演迭代的温度剖面反演误差;计算第k步迭代的温度雅克比矩阵及对应的温度对角矩阵;根据温度对角矩阵计算反演目标参数增量;计算第k+1步的温度剖面和对应的k+1步的反演误差;将k步与k+1步中的反演误差进行比较;若满足迭代终止条件则完成迭代并输出第k+1步的反演参数向量作为稠油油藏注蒸汽水平井反演解释结果,若不满足则从第k步开始重复迭代,直至满足终止条件,本发明提供了一种对稠油油藏生产和量化评价的新思路和新手段,同时也能显著的提升稠油热采的经济效益。
至此,本领域技术人员认识到,虽然本文已详尽展示和描述了本发明的实施例,但是,在不脱离本发明精神和范围的情况下,仍可根据本发明公开的内容直接确定或推导符合本发明原理的许多其他变形或修改。因此,本发明的范围应被理解和认定为覆盖了所有这些其他变形或修改。
Claims (8)
1.一种基于DTS的稠油油藏注蒸汽水平井综合评价方法,其特征在于,所述评价方法包括以下步骤:
8.如权利要求1-7中任意一项所述的一种基于DTS的稠油油藏注蒸汽水平井综合评价方法,其特征在于,所述反演目标参数为吸汽剖面、渗透率分布、蒸汽干度分布、加热半径和蒸汽波及范围。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111642620.1A CN114282387A (zh) | 2021-12-29 | 2021-12-29 | 一种基于dts的稠油油藏注蒸汽水平井综合评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111642620.1A CN114282387A (zh) | 2021-12-29 | 2021-12-29 | 一种基于dts的稠油油藏注蒸汽水平井综合评价方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114282387A true CN114282387A (zh) | 2022-04-05 |
Family
ID=80878102
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111642620.1A Pending CN114282387A (zh) | 2021-12-29 | 2021-12-29 | 一种基于dts的稠油油藏注蒸汽水平井综合评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114282387A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115324546A (zh) * | 2022-08-28 | 2022-11-11 | 重庆科技学院 | Sagd蒸汽预热成效预测方法、模拟实验装置及其模拟方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102278103A (zh) * | 2011-08-25 | 2011-12-14 | 孙洪军 | 一种重力泄水辅助蒸汽驱提高深层超稠油油藏采收率方法 |
CN103527161A (zh) * | 2012-07-03 | 2014-01-22 | 中国石油化工股份有限公司 | 热采水平井堵水调剖的方法 |
CN112302607A (zh) * | 2020-07-07 | 2021-02-02 | 西南石油大学 | 一种致密气藏压裂水平井人工裂缝参数解释方法 |
CN112302606A (zh) * | 2020-07-07 | 2021-02-02 | 西南石油大学 | 一种低渗气藏压裂水平井产出剖面反演解释方法 |
-
2021
- 2021-12-29 CN CN202111642620.1A patent/CN114282387A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102278103A (zh) * | 2011-08-25 | 2011-12-14 | 孙洪军 | 一种重力泄水辅助蒸汽驱提高深层超稠油油藏采收率方法 |
CN103527161A (zh) * | 2012-07-03 | 2014-01-22 | 中国石油化工股份有限公司 | 热采水平井堵水调剖的方法 |
CN112302607A (zh) * | 2020-07-07 | 2021-02-02 | 西南石油大学 | 一种致密气藏压裂水平井人工裂缝参数解释方法 |
CN112302606A (zh) * | 2020-07-07 | 2021-02-02 | 西南石油大学 | 一种低渗气藏压裂水平井产出剖面反演解释方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115324546A (zh) * | 2022-08-28 | 2022-11-11 | 重庆科技学院 | Sagd蒸汽预热成效预测方法、模拟实验装置及其模拟方法 |
CN115324546B (zh) * | 2022-08-28 | 2023-09-08 | 重庆科技学院 | Sagd蒸汽预热成效预测方法、模拟实验装置及其模拟方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2002300917B2 (en) | Method of predicting formation temperature | |
CN109598099B (zh) | 一种考虑油藏与井筒耦合的双管sagd长水平井均匀注汽数值模拟方法 | |
CN112302606B (zh) | 一种低渗气藏压裂水平井产出剖面反演解释方法 | |
CN108266176B (zh) | 一种基于井筒模型的天然气井口流量计算方法 | |
CN104895560A (zh) | 一种深水测试井筒压力、温度场模拟及水合物预测方法 | |
CN111400893B (zh) | 一种套管式地埋管换热器流体温度场分析方法 | |
CN109948182A (zh) | 一种针对中深层地热井井间距的计算方法 | |
CN106770439A (zh) | 岩土层分层导热系数测试方法 | |
CN114282387A (zh) | 一种基于dts的稠油油藏注蒸汽水平井综合评价方法 | |
CN112780257A (zh) | 基于分布式光纤传感钻井液漏失监测系统及监测方法 | |
CN105134179A (zh) | 一种天然气井井筒压力和温度的分布数据的计算方法 | |
CN112302607B (zh) | 一种致密气藏压裂水平井人工裂缝参数解释方法 | |
CN114595504B (zh) | 一种页岩气藏压裂水平井产出剖面反演解释方法 | |
CN115130274A (zh) | 一种稠油油藏相似准则研究方法 | |
CN111894561B (zh) | 一种适用于欠平衡钻井的地层特性随钻解释方法 | |
CN115293066B (zh) | 考虑地层渗流传热效应的气井温度场计算方法 | |
CN104866681A (zh) | 高温高压油气斜井关井过程中温度压力数值模拟方法 | |
CN103615233A (zh) | 一种测定三维井眼中套管轴向力和侧向力的方法及设备 | |
RU2751528C1 (ru) | Способ определения скорости потока жидкости в скважине | |
CN115655769A (zh) | 一种中深层套管换热器换热实验系统及实验方法 | |
US6807854B2 (en) | Method of determining the thermal profile of a drilling fluid in a well | |
Abdelhafiz | Modelling and simulating the thermal interaction between the wellbore fluids and rock formation during drilling operations | |
CN109209357B (zh) | 一种生产测井解释方法 | |
CN113958299B (zh) | 一种基于温度测量的水平井吸水剖面反演方法 | |
RU2753129C1 (ru) | Способ определения поинтервальной скорости и расхода жидкости в скважине |
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 |