CN115492573B - 一种柱塞气举井地层流入动态确定方法 - Google Patents
一种柱塞气举井地层流入动态确定方法 Download PDFInfo
- Publication number
- CN115492573B CN115492573B CN202211455230.8A CN202211455230A CN115492573B CN 115492573 B CN115492573 B CN 115492573B CN 202211455230 A CN202211455230 A CN 202211455230A CN 115492573 B CN115492573 B CN 115492573B
- Authority
- CN
- China
- Prior art keywords
- pressure
- gas
- time
- oil pipe
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 42
- 238000004519 manufacturing process Methods 0.000 claims abstract description 36
- 238000009826 distribution Methods 0.000 claims abstract description 28
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 238000009499 grossing Methods 0.000 claims abstract description 12
- 230000008859 change Effects 0.000 claims abstract description 11
- 239000007788 liquid Substances 0.000 claims description 72
- 239000011541 reaction mixture Substances 0.000 claims description 17
- 230000015572 biosynthetic process Effects 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000006243 chemical reaction Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 8
- 238000006467 substitution reaction Methods 0.000 claims description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 4
- 238000004891 communication Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 239000008398 formation water Substances 0.000 claims description 3
- 238000009530 blood pressure measurement Methods 0.000 claims description 2
- 239000007791 liquid phase Substances 0.000 claims description 2
- 230000003068 static effect Effects 0.000 claims description 2
- 239000012530 fluid Substances 0.000 description 7
- 238000012545 processing Methods 0.000 description 3
- 230000001174 ascending effect Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 239000012071 phase Substances 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000003490 calendering Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
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
- 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
- E21B49/08—Obtaining fluid samples or testing fluids, in boreholes or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/203—Drawing of straight lines or curves
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- Databases & Information Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种柱塞气举井地层流入动态确定方法。它包括下述步骤:S1、测试关井阶段的井口油压和套压数据,并进行光滑处理;S2、根据光滑处理后的井口油压和套压数据得到关井阶段的t~t+1时刻的地层产气量qsc和井筒气水界面分布情况;S3:根据光滑处理后的井口油压和套压数据和井筒气水界面分布计算得到井底流压pwf;S4:根据地层产气量qsc和井底流压pwf做出qsc~pwf曲线,同时考虑近井带含水饱和度变化的产能方程形式拟合得到该井的产能方程,从而得到地层流入动态曲线。本发明可以准确认识柱塞气举井地层流入特征,能较为准确的确定柱塞气举井的地层流入动态曲线。
Description
技术领域
本发明涉及一种页岩气水平井套管,具体涉及一种柱塞气举井地层流入动态确定方法。
背景技术
非常规低渗、致密以及页岩气的产水气井进入到生产中后期,由于地层能量不足,难以带出地层中的产出水,导致井底积液,影响气井产能,严重的还会导致气井“压死”。柱塞气举由于其工艺简单,成本较低是目前较为广泛运用的一种排水采气方法。柱塞气举利用在油管中下入柱塞,周期性开关井带出井筒中的液体,增大生产压差,从而提升气井产气能力。
地层流入动态曲线是气井所有工作的基础,由于柱塞气举工艺和成本限制,难以通过下入井底压力计准确测试井底流压,因此,目前现场生产中,大多数柱塞气举井都未下井底压力计测试井底流压;同时,由于柱塞气举开井生产阶段井筒流动属于不稳定流动,测试的产量并不是地层流入井筒的产量,即井口产量并能不完全反映地层流入动态,只能通过井口压力变化反推地层流入动态。但是,在柱塞上行阶段中柱塞将井筒分隔成两个压力系统,计算井筒气量和气液分布较难,难以准确获取该阶段的地层流入动态,导致柱塞气举井的地层流入动态曲线(IPR曲线)难以确定,使得柱塞气举井缺少确定合理工作方式的主要依据,以及缺少分析油气井动态的基础。
本申请人发现现有技术至少存在以下技术问题:
现有技术中由于柱塞气举开井生产阶段井筒流动属于不稳定流动,测试的产量并不是地层流入井筒的产量,柱塞气举井的地层流入动态曲线(IPR曲线)难以确定。
发明内容
本发明的目的在于提供一种柱塞气举井地层流入动态确定方法,以解决现有技术由于柱塞气举开井生产阶段井筒流动属于不稳定流动,测试的产量并不是地层流入井筒的产量,柱塞气举井的地层流入动态曲线(IPR曲线)难以确定的技术问题。
为实现上述目的,本发明提供了以下技术方案:
本发明提供的一种柱塞气举井地层流入动态确定方法,利用关井阶段测试的井口油套压计算井底流压,包括下述步骤:
S1、输入测试关井阶段的井口油套压数据;
S2、根据井口油套压数据得到关井阶段的t~t+1时刻的地层产气量qsc;
S3:根据井口油套压数据计算得到t~t+1时刻井底流压pwf;
S4:根据地层产气量qsc和井底流压pwf绘制qsc~pwf曲线。
进一步的,还包括对所述油套压数据进行光滑处理的步骤。
进一步的,还包括井筒气水界面分布计算步骤。
进一步的,还包括含水饱和度计算步骤。
进一步的,还包括井筒气体参数计算步骤。
进一步的,所述油套压数据包括油管流压数据和套管流压数据。
基于上述技术方案,本发明实施例至少可以产生如下技术效果:
本发明提供的柱塞气举井地层流入动态确定方法,利用关井阶段测试的井口油套压计算井底流压,同时分别利用油管、环空体积和油、套压压力变化计算地层流入油套管中的气量,从而获得地层流入动态曲线;该方法能较为准确的确定柱塞气举井的地层流入动态曲线(IPR曲线),为准确认识柱塞气举井地层流入特征奠定基础,对柱塞在井筒中的运动模拟以及柱塞气举制度的制定和优化具有重要意义。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见的,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图示出的结构获得其他的附图。
图1是本发明各个阶段压力和产量变化曲线;
图2是本发明实施例中Z井测试的油套压和产气量曲线;
图3是本发明实施例中Z井油压的压力导数曲线;
图4是本发明实施例中Z井套压的压力导数曲线;
图5是本发明实施例中Z井油套压光滑前后的压力曲线对比;
图6是本发明实施例中Z井油压的压力导数曲线;
图7是本发明实施例中Z井套压的压力导数曲线;
图8是本发明实施例中在某一时刻下的井筒中的气液界面分布。
具体实施方式
应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。另外,各个实施例之间的技术方案可以相互结合,但是必须是以本领域普通技术人员能够实现为基础,当技术方案的结合出现相互矛盾或无法实现时应当人认为这种技术方案的结合不存在,也不在本发明要求的保护范围之内。
常规柱塞气举井可以分为三个周期:关井压力恢复、开井柱塞上行和续流阶段,各个阶段压力和产量变化曲线如图1所示;
本发明以关井时刻作为柱塞气举井一个周期的开始,关井压力恢复阶段可认为井筒中的柱塞不影响井筒压力变化,地层产出的气和水分别进入到油管和环空中,使得油套压逐渐上升。当油套压上升到一定值后,开井生产,柱塞在压差作用下从井底向井口举升,此时柱塞将自身上下分隔为两个独立的压力系统,井口产出气量和水量并不与此时的地层供给量相等。柱塞上部液段完全产出,柱塞到达井口捕捞装置后,继续开井生产的阶段为续流阶段。本发明以关井时刻作为柱塞气举井一个周期的开始,关井阶段井筒中气液运动较缓、压力稳定变化,可通过该阶段的井筒压力变化反推地层流入动态,以此作为该井的地层流入动态曲线。
本发明中所有的公式,t时刻和t+1时刻通用的公式未使用上标,上标为t的代表t时刻对应的数据,上标为t+1的代表t+1时刻对应的数据。
下述以Z井为例具体说明柱塞气举井地层流入动态的确定方法,Z井为四川长宁区块页岩气水平井,属于下倾型井,井深4350.12m,使用鱼骨柱塞进行排水采气,柱塞外径48.2mm,油管内径50.6mm。具体包括下述步骤:
S1、测试关井阶段的井口油压和套压数据,并进行光滑处理;
目前现场跟踪测试的井口油套压数据通常为30s测试一次,但是由于测试仪表精度和稳定性问题,测试的油套压数据存在一定的波动,虽然从油套压数值上反映不明显,但是转换成压力导数就能明显观察到这一波动现象,例如Z井测试的油套压曲线如图2所示,油套压的压力导数曲线如图3、图4所示。
从测试的油套压图2可以看出,关井阶段油套压曲线是连续光滑上升的,但是从压力导数曲线图3、图4可以看出,压力导数为点阵而不是光滑的曲线,说明测试的压力点并不连续光滑。为此需要对测试的油压和套压数据进行拟合光滑处理;对井口油压和套压数据进行光滑处理采用下式(1)进行:
式(1)中:y为光滑处理后的井口油压或套压数据,Pa;
x为步骤S1中测试的关井阶段的井口油压或套压数据,Pa;
n 1、n 2、n 3、n 4为拟合参数,使用测试的关井阶段的井口油压和套压数据,在数据处理软件Origin中,拟合得到n 1、n 2、n 3、n 4计算出光滑处理后的数据,光滑后的油套压曲线如图5所示。
Z井光滑使用的参数和光滑前后的相关系数R2表如表1,相关系数接近于1,拟合度很高。
表1 数据处理使用的参数及决定系数表
Z井光滑前后油套压的压力导数曲线对比如图6、图7所示,压力导数为光滑曲线。
S2、根据步骤S1中光滑处理后的井口油压和套压数据得到关井阶段的t~t+1时刻的地层产气量qsc和井筒气水界面分布情况;
S21、t~t+1时刻的地层产气量qsc的获得过程为:
①基于油套连通位置的压力相等原则,由于水是不可压缩流体,则可以认为当环空中的压力过高时,环空液体向油管转移,反之,当油管中压力过高油管中液体向环空中转移,参与转移的流体只有液体,气体则按分配直接进入环空和油管,在环空液面下降到油套连接处前均不会参与转移。在环空液面下降到油套连接处时,即环空无液面时,油、套管的压力变化直接反映了地层产出气在油管和环空中的分配。
任意时刻油管和环空中的总气量折算到地面条件下,由下式(2)得到:
式(2)中,ρ g为当前平均压力下的气体密度,kg/m3;ρ g-air为地面条件下的气体密度,kg/m3;A c为套管内径,m;A to为油管外径,m;A t为油管内径,m;V cg为换算到地面条件下的环空总气量,m3;V tg为换算到地面条件下的油管总气量,m3;H tw为油管内液柱高度,m;H ctw为环空液柱高度,m;Hs为油管鞋深度,m;
②那么从t时刻到t+1时刻的地层产气量,由下式(3)得到:
S22、t~t+1时刻井筒气水界面分布情况的获得过程为:
①在某一时刻下的井筒中的气液界面分布如下图8所示:
不论井筒中的气水分布如何,任意时刻都必须满足油套连通位置的压力相等,即满足下式(4):
②关井时期,井筒气体流动非常缓慢,油管和环空液面上的压力可通过气柱静压计算,如下式(5)表述:
式(5)中,为任意时刻环空液面上压力,Pa;为任意时刻油管液面上部压
力,Pa;为套压,Pa;环空气柱高度,m;为环空平均温度,K;为环空平均偏差
因子;为气体相对密度;为油压,Pa;油管气柱高度,m;为油管平均温度,K;
为油管平均偏差因子;由于确定平均压力下的偏差因子需要用到井筒压力分布,而使用式
(5)计算井筒压力分布,所以必须使用试算的方法计算井筒气体参数。
③阶段地层产水量由下式(6)表述:
④考虑地层的液相流入,从t时刻到t+1时刻对环空和油管中的液量具有下式(8)的统一物质平衡方程:
⑤气水界面分布与每个时间段下的阶段产水量有关,通过迭代的方法逐时刻计算气水界面分布情况,计算流程图如下所示:
对应的具体计算步骤为:
Step2:根据式(6)计算出t时刻的地层产水量V Lt;
S3:根据步骤S1中光滑处理后的井口油压和套压数据和步骤S2中的井筒气水界面分布计算得到井底流压pwf;井底流压pwf的获得过程为:
现场测压数据为30s一个测试点,t~t+1时刻的井底流压变化较小,使用t时刻的井底流压代表t~t+1时刻的井底流压,利用步骤S2得到的每一时刻的地层产气量,对应的井底流压由下式(9)表示:
式(9)中:为t时刻的井底流压,Pa;为t时刻环空液面上压力,Pa;
为t时刻的环空液面高度,m;为液体密度,kg/m3;g为重力系数,N/kg;为油管鞋到井
底的高度,m;其中, 可通过下式(5)t时刻套压计算:
式(5)中,为任意时刻环空液面上压力,Pa;为任意时刻油管液面上部压
力,Pa;为套压,MPa;H ctg环空气柱高度,m;为环空平均温度,K;为环空平均偏差
因子;为气体相对密度;为油压,Pa;油管气柱高度,m;为油管平均温度,K;
为油管平均偏差因子。
井筒中的气体密度、偏差因子分布受到井筒压力分布的影响,同时计算井筒压力分布时首先也需要明确井筒中的气体密度、偏差因子的分布,针对这一问题,采用试算法计算井筒中的压力、气体密度和偏差因子的分布。
式(9)中的气体密度和式(5)中的偏差因子通过1978年PR-Peneloux修正方程计算,计算偏差因子使用的平均压力为下式(12):
基于平均压力的井筒气体参数计算,计算环空压力则p代入p c,计算油管压力则p代入p t,
具体的计算步骤为:
步骤1:假设需要计算的液柱上部压力等于此时的井口压力,初始环空液柱上部压力p ctwu_0=p ct ,初始油管液柱上部压力p twu_0=p t ,
式(12)中,p cavg为环空平均压力,Pa;p tavg为油管平均压力,Pa;p c为套压,Pa;p t为油压,Pa;p ctwu为环空液面上压力,Pa;p twu为油管液面上部压力,Pa;
步骤6:气体密度和偏差因子通过1978年PR-Peneloux修正方程计算,将p ctwu和p c代入式(12)计算出平均压力;从而可以相应的计算出气体密度和偏差因子;
其中,下标_i指的是第i次迭代的参数,下标_o指的是假设的迭代初值;
计算油管气体的密度和偏差因子则代入油管相关参数,计算环空气体的密度和偏差因子则代入环空相关参数。
S4:根据步骤S2中得到的t~t+1时刻的地层产气量qsc和步骤S3中得到的井底流压pwf做出qsc~pwf曲线,同时产水气井引入含水饱和度的影响,以考虑近井带含水饱和度变化的产能方程形式拟合得到该井的产能方程,从而得到地层流入动态曲线;地层流入动态曲线的获得过程为:
①根据步骤S2中得到的t~t+1时刻的地层产气量qsc和步骤S3中得到的井底流压pwf做出qsc~pwf曲线,得到指数式产能方程,如下式(10):
式(10)中,C为系数;n为指数;q sc为t~t+1时刻的地层产气量;p wf为井底流压;
②引入含水饱和度的影响,得到产水气井的地层流入曲线,如下式(11):
式(11)中,c1、c2、c3、c4均为系数。
上式中的C,n,γ以及c1、c2、c3、c4均可以通过数据处理软件Origin拟合得到。
Z井地层流入动态方程拟合结果如表2所示;
表2 Z井产能曲线拟合结果
即该柱塞气举井的地层流入动态方程为:
上面对本发明的较佳实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域的普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护的范围由所附的权利要求书及其等效物界定。
Claims (4)
1.一种柱塞气举井地层流入动态确定方法,其特征在于,利用关井阶段测试的井口油套压计算井底流压,包括下述步骤:
S1、输入测试关井阶段的井口油套压数据;
S2、根据井口油套压数据得到关井阶段的t~t+1时刻的地层产气量qsc和井筒气水界面分布情况;
t~t+1时刻的地层产气量qsc的获得过程为:
任意时刻油管和环空中的总气量折算到地面条件下,由下式(2)得到:
式(2)中,ρ g为当前平均压力下的气体密度,kg/m3;ρ g-air为地面条件下的气体密度,kg/m3;A c为套管内径,m;A to为油管外径,m;A t为油管内径,m;V cg为换算到地面条件下的环空总气量,m3;V tg为换算到地面条件下的油管总气量,m3;H tw为油管内液柱高度,m;H ctw为环空液柱高度,m;Hs为油管鞋深度,m;
那么从t时刻到t+1时刻的地层产气量,由下式(3)得到:
t~t+1时刻井筒气水界面分布情况的获得过程为:
①不论井筒中的气水分布如何,任意时刻都必须满足油套连通位置的压力相等,即满足下式(4):
式(4)中,p ctwu为任意时刻环空液面上压力,Pa;p twu为任意时刻油管液面上部压力,Pa;Htw为任意时刻油管内液面高度,m;H ctw为任意时刻环空液面高度,m;ρ w为液体密度,kg/m3;g为重力系数,N/kg;
②关井时期,井筒气体流动非常缓慢,油管和环空液面上的压力通过气柱静压计算,如下式(5)表述:
式(5)中,p ctwu为任意时刻环空液面上压力,Pa;p twu为任意时刻油管液面上部压力,Pa;p c为套压,Pa;H ctg环空气柱高度,m;为环空平均温度,K;为环空平均偏差因子;γ g为气体相对密度;p t为油压,Pa;H tg油管气柱高度,m;为油管平均温度,K;为油管平均偏差因子;由于确定平均压力下的偏差因子需要用到井筒压力分布,而使用式(5)计算井筒压力分布,所以必须使用试算的方法计算井筒气体参数;
③阶段地层产水量由下式(6)表述:
式(6)中,VLt 为 t~t+1 时刻的阶段产水量,m3 ;qL 为地层日产水量,m3/d;Δt为 t~t+1 的时间,s;
④考虑地层的液相流入,从t时刻到t+1时刻对环空和油管中的液量具有下式(8)的统一物质平衡方程:
⑤气水界面分布与每个时间段下的阶段产水量有关,通过迭代的方法逐时刻计算气水界面分布情况;
S3:根据井口油套压数据和井筒气水界面分布计算得到t~t+1时刻井底流压pwf,井底流压数据的获得过程为:
现场测压数据为30s一个测试点,t~t+1时刻的井底流压变化较小,使用t时刻的井底流压代表t~t+1时刻的井底流压,利用步骤S2得到的每一时刻的地层产气量,对应的井底流压由下式(9)表示:
式(9)中:p t wf为t时刻的井底流压,Pa;p t ctwu为t时刻环空液面上压力,Pa;H t ctw为t时刻的环空液面高度,m;ρ w为液体密度,kg/m3;g为重力系数,N/kg;H x为油管鞋到井底的高度,m;其中,p t ctwu 通过下式(5)t时刻套压计算:
式(5)中,p ctwu为任意时刻环空液面上压力,Pa;p twu为任意时刻油管液面上部压力,Pa;p c为套压,MPa;H ctg环空气柱高度,m;为环空平均温度,K;为环空平均偏差因子;γ g为气体相对密度;p t为油压,Pa;H tg油管气柱高度,m;为油管平均温度,K;为油管平均偏差因子;
S4:根据地层产气量qsc和井底流压pwf绘制qsc~pwf曲线,同时产水气井引入含水饱和度的影响,以考虑近井带含水饱和度变化的产能方程形式拟合得到该井的产能方程,从而得到地层流入动态曲线;地层流入动态曲线的获得过程为:
①根据步骤S2中得到的t~t+1时刻的地层产气量qsc和步骤S3中得到的井底流压pwf做出qsc~pwf曲线,得到指数式产能方程,如下式(10):
式(10)中,C为系数;n为指数;q sc为t~t+1时刻的地层产气量;p wf为井底流压;
②引入含水饱和度的影响,得到产水气井的地层流入曲线,如下式(11):
式(11)中,c1、c2、c3、c4均为系数, C,n,γ以及c1、c2、c3、c4均通过数据处理软件Origin拟合得到。
3.根据权利要求1所述的柱塞气举井地层流入动态确定方法,其特征在于,还包括井筒气体参数计算步骤,具体的计算步骤为:
步骤1:假设需要计算的液柱上部压力等于此时的井口压力,初始环空液柱上部压力p ctwu_0=p ct ,初始油管液柱上部压力p twu_0=p t ,
步骤2:根据下式(12)计算出初始平均压力p cavg_0及p tavg_0;
式(12)中,p cavg为环空平均压力,Pa;p tavg为油管平均压力,Pa;p c为套压,Pa;p t为油压,Pa;p ctwu为环空液面上压力,Pa;p twu为油管液面上部压力,Pa;
步骤5:设置迭代精度a ct =|p ctwu_0-p ctwu_i|,如果a ct >1e-6那么令p ctwu_0=p ctwu_i重复步骤2-步骤4,如果a ct ≤1e-6那么令p ctwu=p ctwu_i退出循环;设置迭代精度a t =|p twu_0-p twu_i|,如果a t >1e-6那么令p twu_0=p twu_i重复步骤1- 步骤3,如果a t ≤1e-6那么令p twu=p twu_i退出循环,计算得到油管液柱上部压力;
步骤6:气体密度和偏差因子通过1978年PR-Peneloux修正方程计算,将p ctwu和p c代入式(12)计算出平均压力;从而相应的计算出气体密度和偏差因子;
其中,下标_i指的是第i次迭代的参数,下标_o指的是假设的迭代初值;
计算油管气体的密度和偏差因子则代入油管相关参数,计算环空气体的密度和偏差因子则代入环空相关参数。
4.根据权利要求1所述的柱塞气举井地层流入动态确定方法,其特征在于,所述油套压数据包括油管流压数据和套管流压数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211455230.8A CN115492573B (zh) | 2022-11-21 | 2022-11-21 | 一种柱塞气举井地层流入动态确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211455230.8A CN115492573B (zh) | 2022-11-21 | 2022-11-21 | 一种柱塞气举井地层流入动态确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115492573A CN115492573A (zh) | 2022-12-20 |
CN115492573B true CN115492573B (zh) | 2023-03-17 |
Family
ID=85114811
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211455230.8A Active CN115492573B (zh) | 2022-11-21 | 2022-11-21 | 一种柱塞气举井地层流入动态确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115492573B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104806229A (zh) * | 2015-04-01 | 2015-07-29 | 西南石油大学 | 柱塞气举排水采气生产测井系统及其控制方法 |
CN112199820A (zh) * | 2020-09-05 | 2021-01-08 | 黑龙江省荣泽石油设备有限公司 | 一种数字化条件下油井产能曲线测试方法 |
CN114329840A (zh) * | 2021-12-31 | 2022-04-12 | 四川轻化工大学 | 一种柱塞举升排水采气优化设计计算方法 |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7243730B2 (en) * | 2004-12-31 | 2007-07-17 | Casey Danny M | Well production optimizing system |
US20080240930A1 (en) * | 2005-10-13 | 2008-10-02 | Pumpwell Solution Ltd | Method and System for Optimizing Downhole Fluid Production |
US9500067B2 (en) * | 2011-10-27 | 2016-11-22 | Ambyint Inc. | System and method of improved fluid production from gaseous wells |
CN103590812B (zh) * | 2013-10-21 | 2016-06-08 | 中国石油天然气股份有限公司 | 一种气井积液量的计算方法、计算装置及确定方法 |
CN104790916B (zh) * | 2015-04-24 | 2017-04-19 | 中国石油集团川庆钻探工程有限公司长庆井下技术作业公司 | 一种利用油套压力平衡法排除气井积液方法 |
US9984180B2 (en) * | 2015-05-05 | 2018-05-29 | King Fahd University Of Petroleum And Minerals | Inflow performance relationship for multilateral wells |
CN106097122B (zh) * | 2016-06-16 | 2020-01-07 | 中国石油天然气股份有限公司 | 一种确定井下节流器合理打捞时机的方法 |
CN106869917A (zh) * | 2017-03-24 | 2017-06-20 | 中国石油集团川庆钻探工程有限公司 | 基于井底和井口流入动态曲线共同约束的产能评价方法 |
CN107387030B (zh) * | 2017-07-26 | 2020-04-10 | 中国海洋石油总公司 | 一种利用本井气的接力举升采油设计方法 |
CN108518204B (zh) * | 2018-03-20 | 2020-01-21 | 四川理工学院 | 一种间歇气举排液采气控制方法及系统 |
CN109162697A (zh) * | 2018-10-13 | 2019-01-08 | 中石化石油工程技术服务有限公司 | 一种天然气井油套环空流动动态描述方法 |
CN110439513B (zh) * | 2019-07-30 | 2021-08-31 | 北京雅丹石油技术开发有限公司 | 一种柱塞气举排液采气生产制度的优化方法 |
CN111489090B (zh) * | 2020-04-14 | 2023-06-06 | 中国石油大学(北京) | 井口回压对油井产量影响的分析方法、装置及存储介质 |
CN112630464B (zh) * | 2020-12-01 | 2022-02-11 | 西南石油大学 | 一种柱塞在不同斜度井中下落瞬时速度的测试装置及方法 |
CN115221666A (zh) * | 2021-04-21 | 2022-10-21 | 中国石油化工股份有限公司 | 一种排水采气工艺优选方法及系统 |
-
2022
- 2022-11-21 CN CN202211455230.8A patent/CN115492573B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104806229A (zh) * | 2015-04-01 | 2015-07-29 | 西南石油大学 | 柱塞气举排水采气生产测井系统及其控制方法 |
CN112199820A (zh) * | 2020-09-05 | 2021-01-08 | 黑龙江省荣泽石油设备有限公司 | 一种数字化条件下油井产能曲线测试方法 |
CN114329840A (zh) * | 2021-12-31 | 2022-04-12 | 四川轻化工大学 | 一种柱塞举升排水采气优化设计计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115492573A (zh) | 2022-12-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107506948B (zh) | 一种基于动态泄流体积的页岩油气综合产量分析方法 | |
CN110334431A (zh) | 一种低渗透致密气藏单井控制储量计算及剩余气分析方法 | |
CN108694254B (zh) | 一种变产变压生产气井产量经验递减曲线分析方法 | |
CN109594968A (zh) | 一种页岩气多段压裂水平井压后裂缝参数评价方法及系统 | |
CN111353205B (zh) | 用于致密气藏产水气井地层压力和动态产能的计算方法 | |
CN110344818B (zh) | 柱塞气举井柱塞上行阶段液体段塞界面跟踪方法 | |
CN111287740B (zh) | 基于真实应变下的异常高压气藏动态储量计算方法 | |
CN110847894B (zh) | 一种井下节流气井流压的确定方法 | |
CN105893679A (zh) | 低产水平井续流修正试井解释方法 | |
CN111894561B (zh) | 一种适用于欠平衡钻井的地层特性随钻解释方法 | |
CN115587674B (zh) | 油藏改建储气库扩容达产过程气井动态产能预测方法 | |
CN113236207A (zh) | 一种强非均质性储层中产水气井的定产量递减预测方法 | |
CN115492573B (zh) | 一种柱塞气举井地层流入动态确定方法 | |
CN108627436B (zh) | 一种基于施工动态监测数据的地下水封油库涌水量的预测方法 | |
CN110766562B (zh) | 一种油井间抽时间确定方法 | |
CN112836173B (zh) | 一种页岩油水平井产量计算方法 | |
CN111287741B (zh) | 一种致密油藏体积压裂改造区渗透率的快速计算方法 | |
CN113294147B (zh) | 一种考虑重力因素影响的单洞型断溶体储层试井解释方法 | |
CN114510847A (zh) | 低渗透油藏污染井产能计算方法、电子设备及存储介质 | |
CN112709547A (zh) | 一种产水气井堵水时机的判别方法 | |
CN115828787B (zh) | 一种改进的mdh试井分析方法 | |
CN109784003A (zh) | 一种基于水平井修井管柱的优化组合方法 | |
CN115163015B (zh) | 超稠油蒸汽驱后期产量的调控方法、装置及电子设备 | |
CN115126473B (zh) | 一种页岩气井标准化测试产量的计算方法 | |
CN115906428B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |