CN114526056A - 一种井下节流气井井筒积液高度计算方法 - Google Patents
一种井下节流气井井筒积液高度计算方法 Download PDFInfo
- Publication number
- CN114526056A CN114526056A CN202011225231.4A CN202011225231A CN114526056A CN 114526056 A CN114526056 A CN 114526056A CN 202011225231 A CN202011225231 A CN 202011225231A CN 114526056 A CN114526056 A CN 114526056A
- Authority
- CN
- China
- Prior art keywords
- height
- liquid
- pressure
- choke
- gas
- 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
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
-
- 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Mining
Abstract
本发明涉及一种井下节流气井井筒积液高度计算方法,属于采气工程技术领域,根据油管压降公式及油套环空压降公式计算的井底流压的相等关系,以及油管内的液柱高度与油套环空内的液柱高度具有等比关系,建立了井筒积液高度的计算模型,然后通过采集井下节流器入口压力值和节流器位置环空压力值;将节流器入口压力值和节流器位置环空压力值代入井筒积液高度计算模型,计算油管内的液柱高度和油套环空内的液柱高度,解决现有低压集输气井的井筒内积液高度计算方式准确度低、不可靠的问题。
Description
技术领域
本发明涉及一种井下节流气井井筒积液高度计算方法,属于采气工程技术领域。
背景技术
含水气藏配套排水采气工艺的一个重要指标就是井筒积液程度,根据井筒积液程度,能够优选出经济有效的排水采气工艺。因此,准确的计算出井筒积液高度,对于含水气藏的高效稳定开发具有重要的意义。
含水天然气藏开发方式主要有高压集气与低压集气。对于高压集气气井,目前有多种方法可判断井筒内的积液高度,例如,利用井口高压气产生的次声波来探测环空积液面;或者是利用回声仪来监测环空液面高度;又或者是利用油套压差法来计算井筒内积液高度。例如,申请号为201810836603.3、名称为一种高压气井井筒积液面自动监测装置及方法的中国发明专利申请,公开了利用井口高压气产生的次声波来探测环空积液面的内容;申请号为201511016002.0、名称为一种实时测试气井井筒积液量的方法的中国发明专利申请中,公开了利用回声仪来监测环空液面高度的方法。
但对于低压集输气井,由于井下节流器的存在,采用上述的几种积液高度计算方法计算得到的井筒积液高度并不准确。
发明内容
本发明的目的是提供一种井下节流气井井筒积液高度计算方法,用以解决现有低压集输气井的井筒内积液高度计算方式准确度低、不可靠的问题。
本发明提供了一种井下节流气井井筒积液高度计算方法,包括如下步骤:
采集井下节流器入口压力值和节流器位置环空压力值;
将节流器入口压力值和节流器位置环空压力值代入井筒积液高度计算模型进行计算,所述井筒积液高度计算模型根据以下关系建立:
基于油管压降公式计算得到的井底流压与基于油套环空压降公式计算得到的井底流压是相等的;
油管内的液柱高度与油套环空内的液柱高度成设定比例。
进一步的,所述井筒积液高度计算模型为:
p1+pf1+ρgH1=p2+pf2+ρgH2
H2=kH1
ρg(1-k)H1=(p2-p1)+(pf2-pf1)
式中:H1为油管内的液柱高度,m;p1为节流器入口压力,MPa;pf1为油管内液面深度到节流器深度的流动气水混合物压损,MPa;H2为油套环空内的液柱高度,m;p2为节流器位置环空压力值,MPa;pf2为环空内液面深度到节流器深度的静止气水混合物压损,MPa;ρ为液体密度,kg/m3;g为重力加速度,m/s2,k为环空液柱高度与油管液柱高度的比值。
进一步的,所述计算方法为迭代计算:
1)设定(pf2-pf1)的差值为0;
2)根据相关计算方法,计算节流器入口压力值(p1)、节流器位置环空压力值(p2)、环空液柱高度与油管液柱高度的比值(k);
3)代入井筒积液高度计算模型计算油管积液高度,记为H1-1;
4)反算(pf2_pf1)的值,记为(pf2-pf1)-1;
5)将p1、p2、k、(pf2-pf1)-1再次代入计算模型计算油管积液高度,记为H1-2;
6)若H1-2与H1-1的差值大于设定精度阈值,则(pf2_pf1)初值设为(pf2_pf1)-1,再次重复1)-6)的步骤,直到H1-2与H1-1的差值小于设定精度阈值。
进一步的,节流器入口压力值的采用以下方法获取:
1)构建油管内液柱高度为0m的井筒压力剖面曲线;
2)选取井底流压值,根据管流模型、气液两相嘴流模型,计算得到节流器嘴径;
3)计算所述节流器嘴径与实际节流器嘴径的偏差,若偏差大于设定精度阈值,调整所述井底流压值,返回步骤2);若偏差小于设定精度阈值,根据此时对应的井底流压和管流模型计算得到对应的节流器入口压力。
进一步的,节流器位置环空压力值采用以下静气柱方程计算得到:
进一步的,环空液柱高度与油管液柱高度的比值k的计算公式为:
式中:pts为井口套压,MPa;WGR为液气比,m3/104m3;erf为高斯误差函数。
进一步的,当已知油管积液面深度、节流器位置深度,则油管内液面深度到节流器深度的流动气水混合物压损(pf1)由管流模型计算得到。
进一步的,当已知油套环空积液面深度、节流器位置深度,则环空内液面深度到节流器深度的静止气水混合物压损(pf2)由静气柱方程计算得到。
进一步的,还包括判断计算得到井筒积液高度是否有效的步骤:将计算得到的井筒积液高度与节流器位置的高度进行比较,如果井筒积液高度低于节流器位置的高度,则计算结果有效。
本发明充分结合井下节流气井的生产特点,计算方法简单、可靠地定量确定了井筒的积液程度,有利于优选排水采气工艺,细化工艺制度参数。
附图说明
图1是本发明井筒积液高度计算方法实施例1中的计算流程图;
图2是本发明应用示例计算结果柱状图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明,即所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
以下结合实施例对本发明的特征和性能作进一步的详细描述。
井筒积液高度计算方法实施例1:
本实施例中提供了一种井下节流气井井筒积液高度计算方法,根据井下节流气井的生产特点,从理论研究入手,建立了井筒积液高度的计算模型,用来指导井下节流气井的排水采气工作,保障气井的稳定生产。
为了满足计算需要,本实施例中结合井下节流气井的生产特点,建立了一种具体的井筒积液高度的计算模型。
以井下节流器为起点进行节点分析,根据油管特性和油套环空特性,可将井底流压的环空表达式与井底流压的油管表达式联立,得到节流器入口压力、压损与液柱高度等的关系表达式。
流压的油管表达式:pwf=p1+pf1+ρgH1
流压的环空表达式:pwf=p2+pf2+ρgH2
联立表达式:p1+pf1+ρgH1=p2+pf2+ρgH2
式中:pwf为井底流压,单位MPa;H1为油管内的液柱高度,单位m;p1为节流器入口压力,单位MPa;pf1为油管内液面深度到节流器深度的流动气水混合物压损,单位MPa;H2为油套环空内的液柱高度,单位m;p2为节流器位置环空压力值,单位MPa;pf2为环空内液面深度到节流器深度的静止气水混合物压损,单位MPa;ρ为液体密度,单位kg/m3;g为重力加速度,单位m/s2。
本实施例中,该环空与油管的液柱高度成设定比例,即
H2=kH1
通过将该等比关系与上述油管与环空中压力的相关关系式联立,组成了本实施例中井筒积液高度的计算模型:
ρg(1-k)H1=(p2-p1)+(pf2-pf1)
H2=kH1
其中,k为环空液柱高度与油管液柱高度的比值,无量纲。
基于上述井筒积液高度的计算公式,如图1所示,本实施例中的井下节流气井井筒积液高度计算方法主要包括以下步骤:
1)设定油管动气柱压损与环空静气柱的压损差值为0;
1)采集井下节流器入口压力值和节流器位置环空压力值;
2)将节流器入口压力值和节流器位置环空压力值代入上述井筒积液高度计算模型,计算油管内的液柱高度和/或油套环空内的液柱高度;
3)根据油管内的液柱高度、油套环空内的液柱高度,计算油管动气柱压损、环空静气柱的压损;
4)将节流器入口压力值、节流器位置环空压力值、油管动气柱压损、环空静气柱的压损值代入上述井筒积液高度计算模型,再次计算油管内的液柱高度和/或油套环空内的液柱高度;
5)对比前后两次油管高度值是否满足设定精度阈值,若不满足重复上述计算步骤。
井筒积液高度计算方法实施例2:
本实施例,与上述井筒积液高度计算方法实施例1的区别在于,本实施例中,为了提升计算精度,还根据历史测量数据,总结出了环空液柱高度与油管液柱高度比值的理论计算公式。
根据实测的井筒环空液柱高度、油管液柱高度、井口套压、液气比等数据,建立了计算环空液柱高度与油管液柱高度比值的计算公式:
式中:pts为井口套压,MPa;WGR为液气比,单位m3/104m3;erf为高斯误差函数。
井筒积液高度计算方法实施例3:
上述实施例中,通常直接采集井下节流器入口压力值,本实施例中,还提出了可以根据气液两相管流模型、气液两相嘴流模型,进行节流器入口压力值的拟合,主要包括以下步骤:
(1)构建一条油管内液柱高度为0m的井筒压力剖面,该剖面具有如下特征,即:①产气量、产液量、井筒结构、管流模型、嘴流模型等参数与原井筒完全一致;②因为条件一致,在油管积液面上部井筒,构建的井筒压力剖面曲线与原井筒压力剖面曲线呈平行管线,二者之间的距离与井底流压有关;③特别的,存在一个井底流压值,使得在油管积液面上部井筒,构建的井筒压力剖面与原井筒压力剖面曲线完全重合,二者具有相同的节流器入口压力;
(2)选取一个井底流压值,由管流模型计算出节流器入口压力,并根据井口实际油压,由管流模型计算出节流器出口压力;
(3)根据气液两相嘴流模型,带入构建的油管内液柱高度为0m的井筒压力剖面的节流器入口压力、节流器出口压力,计算出节流器嘴径,并与实际的节流器嘴径对比,计算二者的误差;
(4)以误差精度1%为控制标准,不断调整井底流压计算节流器嘴径,当满足误差控制标准时,判断构建的井筒压力剖面与原井筒压力剖面曲线完全重合,此时得到的节流器入口压力即是真实井筒压力剖面的节流器入口压力。
井筒积液高度计算方法实施例4:
上述实施例中,主要是通过直接采集节流器位置环空压力值的方式获取对应的输入参数,例如通过设置压力传感器/压力计直接采集。
本实施例中,还给出了通过井口套压等参数,直接通过理论公式计算得到节流器位置环空压力值,具体采用以下公式进行计算:
式中:p2为节流器位置环空压力值,单位MPa;pts为井口套压,单位MPa;γg为天然气相对密度,无量纲;H3为井口到节流器位置的垂直深度,单位m;为井筒静气柱平均压缩因子,无量纲;为井筒静气柱平均温度,单位K。由于压缩因子中隐含有所求井底流压,故无法用显示表示,需要采用迭代法求解。
井筒积液高度计算方法实施例5:
本实施例与上述实施例的区别在于,本实施例中,还包括判断计算得到井筒积液高度是否有效的步骤:将计算得到的井筒积液高度与节流器位置的高度进行比较,如果井筒积液高度低于节流器位置的高度,则计算结果有效,可用来指导井下节流气井的排水采气工作,保障气井的稳定生产。
基于上述井筒积液高度计算方法实施例中的计算方法,提供了一种具体应用示例,如图2:
具体的,以东胜气田A井为例,通过本实施例中的井筒积液高度计算方法,分析对应井的井筒积液程度:
步骤1:根据上述建立的计算环空液柱高度与油管液柱高度比值的经验公式,代入气井的套压数据、液气比数据,计算得到环空液柱高度与油管液柱高度比值。
具体的,选取该井2018年6月13日至2019年12月6日之间11个时间节点的气井套压数据、液气比数据,带入经验公式,得到11个时间节点的环空液柱高度与油管液柱高度比值(见表1),从计算结果来看,A井的环空液柱高度与油管液柱高度比值在0.58-0.65之间波动,平均0.62。
步骤2:根据上述节流器入口压力值的拟合过程,获取气井油管管径参数、实际节流器嘴径及深度参数、井口油压参数、日产气量参数、日产液量参数、天然气组分参数,并不断调整井底流压,通过拟合曲线的过程获取节流器入口压力。
具体的,根据拟合步骤,计算了2018年6月13日至2019年12月6日之间11个时间节点的节流器入口压力(见表1)。从计算结果看,该井的节流器入口压力由11.46MPa降至了2.6MPa。
步骤3:根据建立的计算节流器位置环空压力值的计算公式,输入气井节流器深度参数、井口套压参数、天然气组分参数,计算节流器位置环空压力值。
具体的,计算了2018年6月13日至2019年12月6日之间11个时间节点的节流器位置环空压力值(见表1)。从计算结果看,A井节流器位置环空压力由11.55MPa降至3.81MPa。
步骤4:根据建立的井下节流气井井筒积液高度计算公式,按照计算步骤,输入计算得到环空液柱高度与油管液柱高度比值、节流器入口压力值、节流器位置环空压力值,计算出油管积液高度、环空积液高度。
具体的,计算了2018年6月13日至2019年12月6日之间11个时间节点的井筒积液高度,见表1。
步骤5:反算油管动气柱压损、环空静气柱的压损,进而再次计算油管积液高度、环空积液高度,并根据计算精度值,确定计算轮次。
具体的,计算了2018年6月13日至2019年12月6日之间11个时间节点油管动气柱压损与环空静气柱的压损,并再次计算了油管积液高度、套管积液高度。从计算结果看,有7个时间节点第一轮与第二轮计算的井筒积液高度误差值小于5%,满足计算精度要求;有4个节点经过3-5轮的计算,井筒积液高度的误差值满足了计算精度。
步骤6:判断油管内积液面是否在节流器下方,以确定计算结果是否有效。
具体的,A井的地层中部垂深为3170m,节流器垂深为1849m,计算的油管液面深度均大于2320m,液面位于节流器下放,计算结果有效。
以上所述,仅为本发明的较佳实施例,并不用以限制本发明,本发明的专利保护范围以权利要求书为准,凡是运用本发明的说明书及附图内容所作的等同结构变化,同理均应包含在本发明的保护范围内。
表1:
Claims (9)
1.一种井下节流气井井筒积液高度计算方法,其特征在于,包括如下步骤:
采集井下节流器入口压力值和节流器位置环空压力值;
将节流器入口压力值和节流器位置环空压力值代入井筒积液高度计算模型进行计算,所述井筒积液高度计算模型根据以下关系建立:
基于油管压降公式计算得到的井底流压与基于油套环空压降公式计算得到的井底流压是相等的;
油管内的液柱高度与油套环空内的液柱高度成设定比例。
2.根据权利要求1所述的井下节流气井井筒积液高度计算方法,其特征在于,所述井筒积液高度计算模型为:
p1+pf1+ρgH1=p2+pf2+ρgH2
H2=kH1
ρg(1-k)H1=(p2-p1)+(pf2-pf1)
式中:H1为油管内的液柱高度,m;p1为节流器入口压力,MPa;pf1为油管内液面深度到节流器深度的流动气水混合物压损,MPa;H2为油套环空内的液柱高度,m;p2为节流器位置环空压力值,MPa;pf2为环空内液面深度到节流器深度的静止气水混合物压损,MPa;ρ为液体密度,kg/m3;g为重力加速度,m/s2,k为环空液柱高度与油管液柱高度的比值。
3.根据权利要求2所述的井下节流气井井筒积液高度计算方法,其特征在于,所述井筒积液高度计算模型的计算方法为迭代计算:
1)设定(pf2-pf1)的差值为0;
2)根据相关计算方法,计算节流器入口压力值(p1)、节流器位置环空压力值(p2)、环空液柱高度与油管液柱高度的比值(k);
3)代入井筒积液高度计算模型计算油管积液高度,记为H1-1;
4)反算(pf2-pf1)的值,记为(pf2-pf1)-1;
5)将p1、p2、k、(pf2-pf1)-1再次代入计算模型计算油管积液高度,记为H1-2;
6)若H1-2与H1-1的差值大于设定精度阈值,则(pf2-pf1)初值设为(pf2-pf1)-1,再次重复1)-6)的步骤,直到H1-2与H1-1的差值小于设定精度阈值。
4.根据权利要求1所述的井下节流气井井筒积液高度计算方法,其特征在于,节流器入口压力值的采用以下方法获取:
1)构建油管内液柱高度为0m的井筒压力剖面曲线;
2)选取井底流压值,根据管流模型、气液两相嘴流模型,计算得到节流器嘴径;
3)计算所述节流器嘴径与实际节流器嘴径的偏差,若偏差大于设定精度阈值,调整所述井底流压值,返回步骤2);若偏差小于设定精度阈值,根据此时对应的井底流压和管流模型计算得到对应的节流器入口压力。
7.根据权利要求2所述的井下节流气井井筒积液高度计算方法,其特征在于,当已知油管积液面深度、节流器位置深度,则油管内液面深度到节流器深度的流动气水混合物压损(pf1)由管流模型计算得到。
8.根据权利要求2所述的井下节流气井井筒积液高度计算方法,其特征在于,当已知油套环空积液面深度、节流器位置深度,则环空内液面深度到节流器深度的静止气水混合物压损(pf2)由静气柱方程计算得到。
9.根据权利要求2所述的井下节流气井井筒积液高度计算方法,其特征在于,还包括判断计算得到井筒积液高度是否有效的步骤:将计算得到的井筒积液高度与节流器位置的高度进行比较,如果井筒积液高度低于节流器位置的高度,则计算结果有效。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011225231.4A CN114526056A (zh) | 2020-11-05 | 2020-11-05 | 一种井下节流气井井筒积液高度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011225231.4A CN114526056A (zh) | 2020-11-05 | 2020-11-05 | 一种井下节流气井井筒积液高度计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114526056A true CN114526056A (zh) | 2022-05-24 |
Family
ID=81618491
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011225231.4A Pending CN114526056A (zh) | 2020-11-05 | 2020-11-05 | 一种井下节流气井井筒积液高度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114526056A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114996662A (zh) * | 2022-08-08 | 2022-09-02 | 成都英沃信科技有限公司 | 一种柱塞井井筒积液量的确定方法 |
-
2020
- 2020-11-05 CN CN202011225231.4A patent/CN114526056A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114996662A (zh) * | 2022-08-08 | 2022-09-02 | 成都英沃信科技有限公司 | 一种柱塞井井筒积液量的确定方法 |
CN114996662B (zh) * | 2022-08-08 | 2022-11-08 | 西南石油大学 | 一种柱塞井井筒积液量的确定方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101253402B (zh) | 评价欠饱和煤层气储层的方法 | |
CN105089639B (zh) | 一种煤层气井井筒流动动态预测方法 | |
CN104504604B (zh) | 一种定性气井井筒积液的方法 | |
CN110847894B (zh) | 一种井下节流气井流压的确定方法 | |
CN105178943A (zh) | 一种实时校正井筒压力的方法 | |
EP3707345B1 (en) | Determining wellbore leak crossflow rate between formations in an injection well | |
CN106401570A (zh) | 页岩气井产水的确定方法、积液的确定方法及排液方法 | |
CN110735633A (zh) | 一种低渗碳酸盐岩气藏气井井筒积液早期判识方法 | |
CN114526056A (zh) | 一种井下节流气井井筒积液高度计算方法 | |
CN113338915B (zh) | 一种判断气井是否积液及预测积液高度的方法 | |
CN111963151A (zh) | 一种通过抽油机悬点静载荷确定地层压力的方法 | |
CN112796725B (zh) | 一种分段压裂页岩气井压裂段产气贡献率确定方法及系统 | |
CN107725035B (zh) | 一种用于非均匀产液水平井判别来水方向的压力监测方法 | |
CN109522579B (zh) | 水平井压裂施工破裂压力预测方法 | |
CN107339096A (zh) | 一种水平井快速找水方法 | |
CN114060009B (zh) | 一种井下节流天然气井积液判断方法 | |
CN111206919B (zh) | 长井段高产气井储层段井筒压力计算方法 | |
Yao et al. | Case study on diagnosis and identify the degree of bottom hole liquid accumulation in double-branch horizontal wells in PCOC | |
CN107066768A (zh) | 油气上窜速度计算方法 | |
CN114427444A (zh) | 一种自喷期采油井的井底压力预测方法 | |
CN113803054B (zh) | 油水界面深度确定方法及防止生产井水窜的预警方法 | |
CN109869144A (zh) | 一种辨别油藏高压物性参数真伪的方法 | |
CN117436319B (zh) | 一种基于地面示功图的抽油机井生产气油比计算方法 | |
CN110188374B (zh) | 一种井内有气体条件下连续油管井下压力模拟方法 | |
CN117988816A (zh) | 一种井下节流气井积液诊断方法 |
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 |