CN103590818A - 一种井筒内温度分布半解析确定方法及装置 - Google Patents
一种井筒内温度分布半解析确定方法及装置 Download PDFInfo
- Publication number
- CN103590818A CN103590818A CN201310495190.4A CN201310495190A CN103590818A CN 103590818 A CN103590818 A CN 103590818A CN 201310495190 A CN201310495190 A CN 201310495190A CN 103590818 A CN103590818 A CN 103590818A
- Authority
- CN
- China
- Prior art keywords
- section
- pit shaft
- formula
- temperature
- temperature distribution
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000002834 transmittance Methods 0.000 claims description 66
- 239000012530 fluid Substances 0.000 claims description 32
- DWPVVZZGGGCRRM-UHFFFAOYSA-N (4-methoxyphenyl)-(4-methylpiperazin-1-yl)methanone Chemical compound C1=CC(OC)=CC=C1C(=O)N1CCN(C)CC1 DWPVVZZGGGCRRM-UHFFFAOYSA-N 0.000 claims description 14
- 238000004519 manufacturing process Methods 0.000 claims description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 13
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 abstract description 3
- 238000004134 energy conservation Methods 0.000 abstract 2
- 239000003921 oil Substances 0.000 description 34
- 238000010586 diagram Methods 0.000 description 8
- 238000011084 recovery Methods 0.000 description 8
- 239000010779 crude oil Substances 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 239000007788 liquid Substances 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 230000005514 two-phase flow Effects 0.000 description 2
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明提供了一种井筒内温度分布半解析确定方法及装置,所述方法包括:将井筒分成N段,每段井筒内外之间的传热系数为常数;同时,基于能量平衡原理得到每段井筒内外之间的能量守恒表达式,基于传热学得到每段油管内外之间的传热系数与井筒内温度的线性关系;根据每段井筒内外之间的能量守恒表达式、传热系数与温度的线性关系,利用井底已知条件由下往上依次逐段获取井筒内温度分布。
Description
技术领域
本发明涉及油气藏开采领域,特别涉及一种利用地层静态数据和油气井动态数据确定井筒内温度分布半解析方法及装置。
背景技术
原油中的溶解气受温度的影响,在采油过程中,如果井筒内温度过低,气体会从原油中分离出来,一方面井筒中形成两相流,严重影响井的采油效率,另一方面由于气体的分离,原油的粘度增加,增大了井筒中流体流动的阻力。因此井筒温度场的分布对研究井筒流体的流动状态和采油工程有着至关重要的作用。
1、现有技术一的技术方案
该方法综合考虑油层条件、地温梯度、油井产量、井筒管柱结构和掺入液量等因素的影响,基于能量平衡方程得到井筒内及环空内温度分布表达式:
式中:W为地面产出混合液的水当量,W/℃;W1为油层产出液水当量,W/℃;W2为井筒注入液体水当量,W/℃;t为沿井深任一点处注入液体的温度,℃;θ为沿井深任一点处混合液的温度,℃;D为由井口算起沿井筒的深度,m;kl1为油管内流体与环空中流体之间的传热系数,W/m·℃;kl2为环空流体与地层之间的传热系数,W/m·℃;t0为地表年平均温度,℃;m为地温梯度,℃/100m;C1,C2,r1,r2为与边界条件有关的常数。
2、现有技术二的技术方案
该方法基于传热学与能量平衡原理,建立如下方程组:
利用数值方法提出环空掺稀油井筒内温度计算模型。
式中:W为油管内混合液的水当量,W/℃;W2为套管内掺入稀油的水当量,W/℃;T为油管内混合液的温度,℃;t为套管内稀油的温度,℃;l为由井口算起沿井筒的深度,m;kt为油管内流体与环空中流体之间的传热系数,W/(m·℃);kc为环空流体与地层之间的传热系数,W/(m·℃);t0为地表年平均温度,℃;m为地温梯度,℃/100m。
上述现有的井筒内温度分布求解方法,主要存在以下两个问题,一是要求参数较多,求解过程较繁琐;二是求解过程中认为物质的传热系数是个常数,不随物体的温度影响。实际上对于求解井筒内温度的分布规律只需要将井筒以外的作为一个整体来看待即可,同时根据传热学知识,物体的传热系数是温度的函数关系式,而井底和井口内温度相差较大,如果整个井筒的传热系数视为常数来计算井筒内温度分布就存在很大误差。
发明内容
为解决上述问题,本发明提出一种井筒内温度分布半解析确定方法及装置,建立精确计算井筒温度场分布模型,准确预测稠油开采过程中井筒内流体的温度分布能优化采油工艺、节约资源、提高生产效率。
为实现上述目的,本发明提供了一种井筒内温度分布半解析确定方法,所述方法包括:
将井筒分成N段,每段井筒内外之间的传热系数为常数;同时,基于能量平衡原理得到每段井筒内外之间的能量守恒表达式,基于传热学得到每段油管内外之间的传热系数与井筒内温度的线性关系;
根据每段井筒内外之间的能量守恒表达式、传热系数与温度的线性关系,利用井底已知条件由下往上依次逐段获取井筒内温度分布。
可选的,在本发明一实施例中,所述线性关系的表达式为:
k=at+b
式中,k为油管内外之间的传热系数,t为井筒内温度,a、b为已知常数。
可选的,在本发明一实施例中,所述每段井筒内外之间的能量守恒表达式为:
Widti=ki(ti-tei)dhi
式中,i为井筒段的序号,i=1,2,…,N,第i段井筒的区间为(hi-1,hi);Wi为第i井筒段产出液的水当量,W/℃;ti为第i井筒段内流体温度,℃;ki为第i井筒段内外之间的传热系数;W/(m·℃),tei为地层温度,℃,根据地面温度tes和地温梯度m并利用tei=tes+mhi计算得到。
可选的,在本发明一实施例中,所述获取井筒内温度分布的步骤包括:
利用井筒内外之间能量守恒表达式获取关于hi的每段井筒内温度分布的通解表达式;其中,所述通解表达式中含有常数系数式Ci;
利用第N段井筒的已知边界条件为:hN=H,tN=tes+mH,根据所述传热系数与温度的线性关系获取第N段井筒的传热系数kN;其中,整个井筒长度为H;
利用所述第N段井筒的传热系数和所述每段井筒内温度分布的通解获取第N段井筒常数系数式CN;
根据每段井筒内温度分布的通解表达式,利用第N段井筒的已知边界条件为:hi=H,tN=tes+mH以及第N段井筒常数系数式CN获取第N段井筒内温度分布解析式tNh;其中,所述第N段井筒内温度分布式为:
式中,hN-1≤h≤hN;当h=hN-1时,由第N段井筒内温度分布式得到上一段井筒下端点的温度;
按照获取所述第N段井筒内温度分布式的步骤可得到对应井筒段内温度分布式,即根据第i段井筒内温度ti的通解表达式及上一段井筒下端点的温度可得到第i段内温度分布式,由下往上依次求解即得到井筒内的温度分布。
可选的,在本发明一实施例中,所述第i段井筒内温度ti的通解表达式:
可选的,在本发明一实施例中,所述第N段井筒常数系数式CN表达式为:
为实现上述目的,本发明还提供了一种井筒内温度分布半解析确定装置,所述装置包括:
井筒预处理单元,用于将井筒分成N段,每段井筒内外之间的传热系数为常数;同时,基于能量平衡原理得到每段井筒内外之间的能量守恒表达式,基于传热学得到每段油管内外之间的传热系数与井筒内温度的线性关系;
井筒内温度分布获取单元,用于根据每段井筒内外之间的能量守恒表达式、传热系数与温度的线性关系,利用井底已知条件由下往上依次逐段获取井筒内温度分布。
可选的,在本发明一实施例中,所述井筒预处理单元中采用的线性关系的表达式为:
k=at+b
式中,k为油管内外之间的传热系数,t为井筒内温度,a、b为已知常数。
可选的,在本发明一实施例中,所述井筒预处理单元中采用的每段井筒内外之间的能量守恒表达式为:
Widti=ki(ti-tei)dhi
式中,i为井筒段的序号,i=1,2,?,N,第i段井筒的区间为(hi-1,hi);Wi为第i井筒段产出液的水当量,W/℃;ti为第i井筒段内流体温度,℃;ki为第i井筒段内外之间的传热系数;W/(m·℃),tei为地层温度,℃,根据地面温度tes和地温梯度m并利用tei=tes+mhi计算得到。
可选的,在本发明一实施例中,所述井筒内温度分布获取单元包括:
每段井筒内温度通解表达式获取模块,用于利用井筒内外之间能量守恒表达式获取关于hi的每段井筒内温度分布的通解表达式;其中,所述通解表达式中含有常数系数式Ci;
第N段井筒的传热系数获取模块,用于利用第N段井筒的已知边界条件为:hi=H,tN=tes+mH,根据所述传热系数与温度的线性关系获取第N段井筒的传热系数kN;其中,整个井筒长度为H;
第N段井筒常数系数式获取模块,用于利用所述第N段井筒的传热系数和所述每段井筒内温度分布的通解获取第N段井筒常数系数式CN;
第N段井筒内温度分布式获取模块,用于根据每段井筒内温度分布的通解表达式,利用第N段井筒的已知边界条件为:hi=H,tN=tes+mH以及第N段井筒常数系数式CN获取第N段井筒内温度分布解析式tNh;其中,所述第N段井筒内温度分布式为:
式中,hN-1≤h≤hN;当h=hN-1时,由第N段井筒内温度分布式得到上一段井筒下端点的温度;
循环模块,用于按照获取所述第N段井筒内温度分布式的步骤可得到对应井筒段内温度分布式,即根据第i段井筒内温度ti的通解表达式及上一段井筒下端点的温度可得到第i段内温度分布式,由下往上依次求解即得到井筒内的温度分布。
可选的,在本发明一实施例中,所述每段井筒内温度通解表达式获取模块得到的每段井筒内温度分布的通解表达式为:
可选的,在本发明一实施例中,所述第N段井筒常数系数式获取模块获取的第N段井筒常数系数式CN表达式为:
上述技术方案具有如下有益效果:将井筒内外视为两个独立的热力系统,尤其是将井筒外视为一个热力系统,减少了目前求取井筒内温度分布所需的参数数量和简化了求解过程;两个热力系统之间的传热系数随温度的变化而变化,同时利用系统节点思想将非线性方程转化为线性方程求解,通过从井底往上逐段求取各区段的温度解析表达式,从而得到整个井筒内温度分布值。在理论上能更加准确的计算预测出井筒内温度分布规律。进一步根据金属材料的传热性,考虑油管的传热系数随井筒内温度的不同而变化,使得求解结果能更接近和反映真实的传热物理现象,有利于优化采油工艺和提高生产效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提出的一种井筒内温度分布半解析确定方法流程图;
图2为本发明提出的一种井筒内温度分布半解析确定装置框图;
图3为本发明提出的一种井筒内温度分布半解析确定装置中井筒内温度分布获取单元框图;
图4为实施例中开采井截面示意图;
图5为实施例中应用本技术方案得到的井筒温度分布与实际采油井温度测量值对比图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的技术方案的工作原理为:针对目前求解井筒内温度分布存在的问题,基于传热学和能量守恒原理,将井筒内和井筒外看成两个热力系统,井筒内流体通过油管向外系统传递热量,并考虑传热系数为温度的应变函数。同时利用节点思想沿井筒方向将研究的井筒系统分成足够小的若干段,此时每段的两端点的温度相差不大,因此可以将每段的传热系数视为该段温度下的常数。然后由井底已知条件由下往上逐段利用解析方法求解温度的分布值。
如图1所示,为本发明提出的一种井筒内温度分布半解析确定方法流程图。所述方法包括:
步骤101):将井筒分成N段,每段井筒内外之间的传热系数为常数;同时,基于能量平衡原理得到每段井筒内外之间的能量守恒表达式,基于传热学得到每段油管内外之间的传热系数与井筒内温度的线性关系;
在本技术方案中,将井筒分成N段可以有多种方式,比如:平均分段。还有,可以随机分段。但是,在工程实践中,随机分段时,每一段的端点的位置是可以已知的。
步骤102):根据每段井筒内外之间的能量守恒表达式、传热系数与温度的线性关系,利用井底已知条件由下向上逐段获取井筒内温度分布。
可选的,在本发明一实施例中,所述线性关系的表达式为:
k=at+b
式中,k为油管内外之间的传热系数,t为井筒内温度,a、b为已知常数。
可选的,在本发明一实施例中,所述每段井筒内外之间的能量守恒表达式为:
Widti=ki(ti-tei)dhi
式中,i为井筒段的序号,i=1,2,…,N,第i段井筒的区间为(hi-1,hi);Wi为第i井筒段产出液的水当量,W/℃;ti为第i井筒段内流体温度,℃;ki为第i井筒段内外之间的传热系数;W/(m·℃),tei为地层温度,℃,根据地面温度tes和地温梯度m并利用tei=tes+mhi计算得到。
可选的,在本发明一实施例中,所述获取井筒内温度分布的步骤包括:
利用井筒内外之间能量守恒表达式获取关于hi的每段井筒内温度分布的通解表达式;其中,所述通解表达式中含有常数系数式Ci;
利用第N段井筒的已知边界条件为:hi=H,tN=tes+mH,根据所述传热系数与温度的线性关系获取第N段井筒的传热系数kN;其中,整个井筒长度为H;
利用所述第N段井筒的传热系数和所述每段井筒内温度分布的通解获取第N段井筒常数系数式CN;
根据每段井筒内温度分布的通解表达式,利用第N段井筒的已知边界条件为:hi=H,tN=tes+mH以及第N段井筒常数系数式CN获取第N段井筒内温度分布解析式tNh;其中,所述第N段井筒内温度分布式为:
式中,hN-1≤h≤hN;当h=hN-1时,由第N段井筒内温度分布式得到上一段井筒下端点的温度;
按照获取所述第N段井筒内温度分布式的步骤可得到对应井筒段内温度分布式,即根据第i段井筒内温度ti的通解表达式及上一段井筒下端点的温度可得到第i段内温度分布式,由下往上依次求解即得到井筒内的温度分布。
可选的,在本发明一实施例中,所述第i段井筒内温度ti的通解表达式:
可选的,在本发明一实施例中,所述第N段井筒常数系数式CN表达式为:
如图2所示,为本发明提出的一种井筒内温度分布半解析确定装置框图。所述装置包括:
井筒预处理单元201,用于将井筒分成N段,每段井筒内外之间的传热系数为常数;同时,基于能量平衡原理得到每段井筒内外之间的能量守恒表达式,基于传热学得到每段油管内外之间的传热系数与井筒内温度的线性关系;
井筒内温度分布获取单元202,用于根据每段井筒内外之间的能量守恒表达式、传热系数与温度的线性关系,利用井底已知条件由下向上逐段获取井筒内温度分布。
可选的,在本发明一实施例中,所述井筒预处理单元201中采用的线性关系的表达式为:
k=at+b
式中,k为油管内外之间的传热系数,t为井筒内温度,a、b为已知常数。
可选的,在本发明一实施例中,所述井筒预处理单元201中采用的每段井筒内外之间的能量守恒表达式为:
Widti=ki(ti-tei)dhi
式中,i为井筒段的序号,i=1,2,…,N,第i段井筒的区间为(hi-1,hi);Wi为第i井筒段产出液的水当量,W/℃;ti为第i井筒段内流体温度,℃;ki为第i井筒段内外之间的传热系数;W/(m·℃),tei为地层温度,℃,根据地面温度tes和地温梯度m并利用tei=tes+mhi计算得到。
如图3所示,为本发明提出的一种井筒内温度分布半解析确定装置中井筒内温度分布获取单元框图。所述井筒内温度分布获取单元202包括:
每段井筒内温度通解表达式获取模块2021,用于利用井筒内外之间能量守恒表达式获取关于hi的每段井筒内温度分布的通解表达式;其中,所述通解表达式中含有常数系数式Ci;
第N段井筒的传热系数获取模块2022,用于利用第N段井筒的已知边界条件为:hi=H,tN=tes+mH,根据所述传热系数与温度的线性关系获取第N段井筒的传热系数kN;其中,整个井筒长度为H,
第N段井筒常数系数式获取模块2023,用于利用所述第N段井筒的传热系数和所述每段井筒内温度分布的通解获取第N段井筒常数系数式CN;
第N段井筒内温度分布式获取模块2024,用于根据每段井筒内温度分布的通解表达式,利用第N段井筒的已知边界条件为:hi=H,tN=tes+mH以及第N段井筒常数系数式CN获取第N段井筒内温度分布解析式tNh;其中,所述第N段井筒内温度分布式为:
式中,hN-1≤h≤hN;当h=hN-1时,由第N段井筒内温度分布式得到上一段井筒下端点的温度;
循环模块2025,用于按照获取所述第N段井筒内温度分布式的步骤可得到对应井筒段内温度分布式,即根据第i段井筒内温度ti的通解表达式及上一段井筒下端点的温度可得到第i段内温度分布式,由下往上依次求解即得到井筒内的温度分布。
可选的,在本发明一实施例中,所述每段井筒内温度通解表达式获取模块2021得到的每段井筒内温度分布的通解表达式为:
可选的,在本发明一实施例中,所述第N段井筒常数系数式获取模块2023获取的第N段井筒常数系数式CN表达式为:
实施例:
本实施例技术方案实现的具体步骤包括:
1)基础参数:根据钻完井资料得到井的深度H,m;通过井的测试资料获得地温梯度m和地面温度tes,℃;通过生产数据得到产出液的水当量W,W/℃;基于传热学知识得到油管内外之间的传热系数k与温度t的线性关系,k=at+b。
2)如图4所示,为实施例中开采井截面示意图。利用节点思想,沿井筒将研究系统分成N段,共有N+1个节点。在第i段,即在区间(hi-1,hi)内,根据能量平衡有:
Widti=ki(ti-tei)dhi (1)
式中,i为井筒段的序号,i=1,2,…,N,第i段井筒的区间为(hi-1,hi);Wi为第i井筒段产出液的水当量,W/℃;ti为第i井筒段内流体温度,℃;ki为第i井筒段内外之间的传热系数;W/(m·℃),tei为地层温度,℃,根据地面温度tes和地温梯度m并利用式(2)计算得到。
tei=tes+mhi (2)
3)根据微分方程求解方法对公式(1)进行解方程,得到式(1)的通解为:
式(3)是一个非线性方程,当系统划分的节点数足够多时,每个区间的传热系数可以视为常数,此时即将非线性方程(3)化为线性方程。
当i=N时,即区间(hN-1,hN)内有边界条件:hN=H,tN=tes+mH,此时第N区间的传热系数kN由步骤1)中的线性关系式k=at+b可求得,带入式(3)可求得常数表达式为:
将已知的边界条件及常数表达式(4)带入式(3)得第N区间温度分布解析表达式:
式中:tNh为第N区间井筒内温度,℃;h为深度,m,hN-1≤h≤hN。当h=hN-1时,从式(5)得到温度值即为第N-1区间的下端点值。
4)从井底逐段往上按照步骤3)依次可以求解得到研究系统每段的温度分布值,即井筒内各点的温度分布。
如图4所示,沿井筒方向流体从井底向井口流动。垂直于井筒方向,由于井筒流体的温度要比同一深度处地层温度高,井筒内流体的部分热量通过油管向油管外系统传播。如图5所示,为实施例中应用本技术方案得到的井筒温度分布与实际采油井温度测量值对比图。曲线C为随深度变化地层温度的变化值。倒三角曲线A为传热系数不随温度变化而计算得到的温度值;正方形散点为实际测试值,曲线B为本发明计算值。从图中可以看出,曲线B比曲线A更与实际测试值一致,能够更好地满足生产和研究精度要求。
原油中的溶解气受温度的影响,在采油过程中,如果井筒内温度过低,气体会从原油中分离出来,一方面井筒中形成两相流,严重影响井的采油效率,另一方面由于气体的分离,原油的粘度增加,增大了井筒中流体流动的阻力。因此井筒温度场的分布对研究井筒流体的流动状态和采油工程有着至关重要的作用。特别是稠油粘度受温度的影响很大,准确预测稠油开采过程中井筒内流体的温度分布能优化采油工艺、节约资源、提高生产效率,因此,本申请的技术方案提出了一种建立精确计算井筒温度场分布模型。
最后应说明的是:上述仅用以说明本发明而并非限制本发明所描述的技术方案;尽管本说明书对本发明已进行了详细的说明,但是,本领域的技术人员仍然可以对本发明进行修改或等同替换,一切不脱离本发明的精神和范围的技术方案及其改进,其均应涵盖在本发明的权利要求范围中。
Claims (12)
1.一种井筒内温度分布半解析确定方法,其特征在于,所述方法包括:
将井筒分成N段,每段井筒内外之间的传热系数为常数;同时,基于能量平衡原理得到每段井筒内外之间的能量守恒表达式,基于传热学得到每段油管内外之间的传热系数与井筒内温度的线性关系;
根据每段井筒内外之间的能量守恒表达式、传热系数与温度的线性关系,利用井底已知条件由下往上依次逐段获取井筒内温度分布。
2.如权利要求1所述的方法,其特征在于,所述线性关系的表达式为:
k=at+b
式中,k为油管内外之间的传热系数,t为井筒内温度,a、b为已知常数。
3.如权利要求2所述的方法,其特征在于,所述每段井筒内外之间的能量守恒表达式为:
Widti=ki(ti-tei)dhi
式中,i为井筒段的序号,i=1,2,…,N,第i段井筒的区间为(hi-1,hi);Wi为第i井筒段产出液的水当量,W/℃;ti为第i井筒段内流体温度,℃;ki为第i井筒段内外之间的传热系数;W/(m·℃),tei为地层温度,℃,根据地面温度tes和地温梯度m并利用tei=tes+mhi计算得到。
4.如权利要求3所述的方法,其特征在于,所述获取井筒内温度分布的步骤包括:
利用井筒内外之间能量守恒表达式获取关于hi的每段井筒内温度分布的通解表达式;其中,所述通解表达式中含有常数系数式Ci;
利用第N段井筒的已知边界条件为:hN=H,tN=tes+mH,根据所述传热系数与温度的线性关系获取第N段井筒的传热系数kN;其中,整个井筒长度为H;
利用所述第N段井筒的传热系数和所述每段井筒内温度分布的通解获取第N段井筒常数系数式CN;
根据每段井筒内温度分布的通解表达式,利用第N段井筒的已知边界条件为:hi=H,tN=tes+mH以及第N段井筒常数系数式CN获取第N段井筒内温度分布解析式tNh;其中,所述第N段井筒内温度分布式为:
式中,hN-1≤h≤hN;当h=hN-1时,由第N段井筒内温度分布式得到上一段井筒下端点的温度;
按照获取所述第N段井筒内温度分布式的过程可得到对应井筒段内温度分布式,即根据第i段井筒内温度ti的通解表达式及上一段井筒下端点的温度可得到第i段内温度分布式,由下往上依次求解即得到井筒内的温度分布。
5.如权利要求4所述的方法,其特征在于,所述第i段井筒内温度ti的通解表达式:
6.如权利要求4所述的方法,其特征在于,所述第N段井筒常数系数式CN表达式为:
7.一种井筒内温度分布半解析确定装置,其特征在于,所述装置包括:
井筒预处理单元,用于将井筒分成N段,每段井筒内外之间的传热系数为常数;同时,基于能量平衡原理得到每段井筒内外之间的能量守恒表达式,基于传热学得到每段油管内外之间的传热系数与井筒内温度的线性关系;
井筒内温度分布获取单元,用于根据每段井筒内外之间的能量守恒表达式、传热系数与温度的线性关系,利用井底已知条件由下往上依次逐段获取井筒内温度分布。
8.如权利要求7所述的装置,其特征在于,所述井筒预处理单元中采用的线性关系的表达式为:
k=at+b
式中,k为油管内外之间的传热系数,t为井筒内温度,a、b为已知常数。
9.如权利要求8所述的装置,其特征在于,所述井筒预处理单元中采用的每段井筒内外之间的能量守恒表达式为:
Widti=ki(ti-tei)dhi
式中,i为井筒段的序号,i=1,2,…,N,第i段井筒的区间为(hi-1,hi);Wi为第i井筒段产出液的水当量,W/℃;ti为第i井筒段内流体温度,℃;ki为第i井筒段内外之间的传热系数;W/(m·℃),tei为地层温度,℃,根据地面温度tes和地温梯度m并利用tei=tes+mhi计算得到。
10.如权利要求9所述的装置,其特征在于,所述井筒内温度分布获取单元包括:
每段井筒内温度通解表达式获取模块,用于利用井筒内外之间能量守恒表达式获取关于hi的每段井筒内温度分布的通解表达式;其中,所述通解表达式中含有常数系数式Ci;
第N段井筒的传热系数获取模块,用于利用第N段井筒的已知边界条件为:hi=H,tN=tes+mH,根据所述传热系数与温度的线性关系获取第N段井筒的传热系数kN;其中,整个井筒长度为H;
第N段井筒常数系数式获取模块,用于利用所述第N段井筒的传热系数和所述每段井筒内温度分布的通解获取第N段井筒常数系数式CN;
第N段井筒内温度分布式获取模块,用于根据每段井筒内温度分布的通解表达式,利用第N段井筒的已知边界条件为:hi=H,tN=tes+mH以及第N段井筒常数系数式CN获取第N段井筒内温度分布解析式tNh;其中,所述第N段井筒内温度分布式为:
式中,hN-1≤h≤hN;当h=hN-1时,由第N段井筒内温度分布式得到上一段井筒下端点的温度;
循环模块,用于按照获取所述第N段井筒内温度分布式的步骤可得到对应井筒段内温度分布式,即根据第i段井筒内温度ti的通解表达式及上一段井筒下端点的温度可得到第i段内温度分布式,由下往上依次求解即得到井筒内的温度分布。
11.如权利要求10所述的装置,其特征在于,所述每段井筒内温度通解表达式获取模块得到的每段井筒内温度分布的通解表达式为:
12.如权利要求10所述的装置,其特征在于,所述第N段井筒常数系数式获取模块获取的第N段井筒常数系数式CN表达式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310495190.4A CN103590818B (zh) | 2013-10-21 | 2013-10-21 | 一种井筒内温度分布半解析确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310495190.4A CN103590818B (zh) | 2013-10-21 | 2013-10-21 | 一种井筒内温度分布半解析确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103590818A true CN103590818A (zh) | 2014-02-19 |
CN103590818B CN103590818B (zh) | 2016-01-06 |
Family
ID=50081114
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310495190.4A Active CN103590818B (zh) | 2013-10-21 | 2013-10-21 | 一种井筒内温度分布半解析确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103590818B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103867174A (zh) * | 2014-02-26 | 2014-06-18 | 中国石油天然气股份有限公司 | 一种稠油水平井吸汽状况分析方法及系统 |
CN104376188A (zh) * | 2014-12-08 | 2015-02-25 | 电子科技大学 | 一种地热井热损失的计算方法 |
CN104389589A (zh) * | 2014-10-11 | 2015-03-04 | 中国石油天然气股份有限公司 | 基于空心杆确定井筒温度场分布的方法以及系统 |
CN104453861A (zh) * | 2014-11-11 | 2015-03-25 | 中国石油天然气股份有限公司 | 一种高压气井井筒温度分布的确定方法以及系统 |
CN104806230A (zh) * | 2015-02-16 | 2015-07-29 | 中国石油天然气股份有限公司 | 超临界注汽井的井筒温度场及压力场分布的计算方法 |
CN105003238A (zh) * | 2015-07-24 | 2015-10-28 | 中国石油化工股份有限公司 | 利用井筒压力温度剖面分析井下蒸汽干度方法 |
CN105114060A (zh) * | 2015-08-26 | 2015-12-02 | 中国石油天然气股份有限公司 | 一种计算注气井的井筒温度分布的方法及装置 |
CN105370266A (zh) * | 2015-12-01 | 2016-03-02 | 中国石油天然气股份有限公司 | 火烧油层分层电点火注气井井筒温度分布确定方法及装置 |
CN105422084A (zh) * | 2015-11-26 | 2016-03-23 | 中国石油天然气股份有限公司 | 一种热水循环加热降粘工艺的井筒温度场获取方法及装置 |
CN104373117B (zh) * | 2014-10-11 | 2017-01-04 | 中国石油天然气股份有限公司 | 基于空心杆电加热确定井筒温度场分布的方法以及系统 |
CN106503299A (zh) * | 2016-09-27 | 2017-03-15 | 中国石油天然气股份有限公司 | 注蒸汽井井筒总传热系数计算方法 |
CN106934106A (zh) * | 2017-02-17 | 2017-07-07 | 中国海洋石油总公司 | 一种获取固井循环温度的方法及装置 |
CN112211618A (zh) * | 2019-06-25 | 2021-01-12 | 中国石油天然气股份有限公司 | 一种油柱高度的获取方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090032259A1 (en) * | 2007-05-08 | 2009-02-05 | Baker Hughes Incorporated | Downhole applications of composites having aligned nanotubes for heat transport |
CN102587887A (zh) * | 2012-04-05 | 2012-07-18 | 中国石油化工股份有限公司 | 一种气井井况预测方法 |
CN103226641A (zh) * | 2013-05-10 | 2013-07-31 | 中国石油大学(华东) | 深水气液两相流循环温度压力耦合计算方法 |
-
2013
- 2013-10-21 CN CN201310495190.4A patent/CN103590818B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090032259A1 (en) * | 2007-05-08 | 2009-02-05 | Baker Hughes Incorporated | Downhole applications of composites having aligned nanotubes for heat transport |
CN102587887A (zh) * | 2012-04-05 | 2012-07-18 | 中国石油化工股份有限公司 | 一种气井井况预测方法 |
CN103226641A (zh) * | 2013-05-10 | 2013-07-31 | 中国石油大学(华东) | 深水气液两相流循环温度压力耦合计算方法 |
Non-Patent Citations (3)
Title |
---|
曾祥林等: "预测井筒压力及温度分布的机理模型", 《西南石油学院学报(自然科学版)》 * |
杨雄文等: "预测井筒流动温度分布的新方法", 《大庆石油地质与开发》 * |
毛伟等: "计算气井井筒温度分布的新方法", 《西南石油学院学报》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103867174B (zh) * | 2014-02-26 | 2016-08-17 | 中国石油天然气股份有限公司 | 一种稠油水平井吸汽状况分析方法及系统 |
CN103867174A (zh) * | 2014-02-26 | 2014-06-18 | 中国石油天然气股份有限公司 | 一种稠油水平井吸汽状况分析方法及系统 |
CN104389589A (zh) * | 2014-10-11 | 2015-03-04 | 中国石油天然气股份有限公司 | 基于空心杆确定井筒温度场分布的方法以及系统 |
CN104389589B (zh) * | 2014-10-11 | 2017-02-15 | 中国石油天然气股份有限公司 | 基于空心杆确定井筒温度场分布的方法以及系统 |
CN104373117B (zh) * | 2014-10-11 | 2017-01-04 | 中国石油天然气股份有限公司 | 基于空心杆电加热确定井筒温度场分布的方法以及系统 |
CN104453861A (zh) * | 2014-11-11 | 2015-03-25 | 中国石油天然气股份有限公司 | 一种高压气井井筒温度分布的确定方法以及系统 |
CN104453861B (zh) * | 2014-11-11 | 2017-03-08 | 中国石油天然气股份有限公司 | 一种高压气井井筒温度分布的确定方法以及系统 |
CN104376188A (zh) * | 2014-12-08 | 2015-02-25 | 电子科技大学 | 一种地热井热损失的计算方法 |
CN104806230B (zh) * | 2015-02-16 | 2018-01-05 | 中国石油天然气股份有限公司 | 超临界注汽井的井筒温度场及压力场分布的计算方法 |
CN104806230A (zh) * | 2015-02-16 | 2015-07-29 | 中国石油天然气股份有限公司 | 超临界注汽井的井筒温度场及压力场分布的计算方法 |
CN105003238A (zh) * | 2015-07-24 | 2015-10-28 | 中国石油化工股份有限公司 | 利用井筒压力温度剖面分析井下蒸汽干度方法 |
CN105114060A (zh) * | 2015-08-26 | 2015-12-02 | 中国石油天然气股份有限公司 | 一种计算注气井的井筒温度分布的方法及装置 |
CN105114060B (zh) * | 2015-08-26 | 2018-02-02 | 中国石油天然气股份有限公司 | 一种计算注气井的井筒温度分布的方法及装置 |
CN105422084A (zh) * | 2015-11-26 | 2016-03-23 | 中国石油天然气股份有限公司 | 一种热水循环加热降粘工艺的井筒温度场获取方法及装置 |
CN105422084B (zh) * | 2015-11-26 | 2018-07-13 | 中国石油天然气股份有限公司 | 一种热水循环加热降粘工艺的井筒温度场获取方法及装置 |
CN105370266A (zh) * | 2015-12-01 | 2016-03-02 | 中国石油天然气股份有限公司 | 火烧油层分层电点火注气井井筒温度分布确定方法及装置 |
CN106503299A (zh) * | 2016-09-27 | 2017-03-15 | 中国石油天然气股份有限公司 | 注蒸汽井井筒总传热系数计算方法 |
CN106934106A (zh) * | 2017-02-17 | 2017-07-07 | 中国海洋石油总公司 | 一种获取固井循环温度的方法及装置 |
CN106934106B (zh) * | 2017-02-17 | 2021-08-10 | 中国海洋石油总公司 | 一种获取固井循环温度的方法及装置 |
CN112211618A (zh) * | 2019-06-25 | 2021-01-12 | 中国石油天然气股份有限公司 | 一种油柱高度的获取方法及装置 |
CN112211618B (zh) * | 2019-06-25 | 2023-11-28 | 中国石油天然气股份有限公司 | 一种油柱高度的获取方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103590818B (zh) | 2016-01-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103590818B (zh) | 一种井筒内温度分布半解析确定方法及装置 | |
CN104594872B (zh) | 一种优化致密气藏压裂水平井裂缝导流能力的方法 | |
CN110210157A (zh) | 一种页岩气藏压裂水平井非稳态扩散下产能计算方法 | |
Davarpanah et al. | A simulation study to control the oil production rate of oil-rim reservoir under different injectivity scenarios | |
Zhao et al. | A simulation method for modified isochronal well testing to determine shale gas well productivity | |
CN104314547A (zh) | Hafd水平井压裂酸化分段优化方法及系统 | |
CN105201484A (zh) | 一种直井分层压裂层段优选及施工参数优化设计方法 | |
CN106202746A (zh) | 模拟多孔介质中水流达西速度的Yeh‑多尺度有限元方法 | |
CN107191169A (zh) | 用于低产低效老油井的重复压裂方法 | |
CN106021659A (zh) | 一种冲蚀-二氧化碳腐蚀耦合作用下的天然气注采井管柱腐蚀速率的确定方法 | |
CN103939066A (zh) | 一种一注多采井组定注水量确定油井产液量的方法 | |
CN110397430A (zh) | 一种预测碳酸盐岩前置液酸压裂缝导流能力分布的方法 | |
CN107313759A (zh) | 低渗稠油油藏直井热采压力分布预测方法及系统 | |
CN106547973A (zh) | 一种稠油油藏注蒸汽过程汽窜通道的定量描述方法 | |
CN106469333B (zh) | 一种低渗稠油油藏水平井热采压力分布预测方法 | |
CN104060985A (zh) | 一种层状油藏调剖堵水堵剂进入深度测试方法及系统 | |
CN112541287A (zh) | 疏松砂岩压裂充填防砂增产调剖一体化设计方法 | |
CN105527198A (zh) | 一种渗流运动规律的测量方法 | |
CN106285600B (zh) | 获取油气压裂增产中井底温度场的方法及装置 | |
CN104615806A (zh) | 一种凝胶与化学剂交替注入驱油数值模拟研究方法 | |
Shen et al. | Spontaneous imbibition in asymmetric branch-like throat structures in unconventional reservoirs | |
CN106097118A (zh) | 一种相对渗透率曲线确定稠油油藏聚合物驱时机范围的方法 | |
CN104343429A (zh) | 一种确定聚合物驱粘度比的方法 | |
CN106050204B (zh) | 一种考虑岩石孔隙压缩系数温敏效应的稠油热采模拟方法 | |
CN106777663A (zh) | 一种考虑天然裂缝的压裂液滤失速度计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |