CN105545285A - 基于隔水管气液两相流识别的深水钻井气侵监测方法 - Google Patents
基于隔水管气液两相流识别的深水钻井气侵监测方法 Download PDFInfo
- Publication number
- CN105545285A CN105545285A CN201510689849.9A CN201510689849A CN105545285A CN 105545285 A CN105545285 A CN 105545285A CN 201510689849 A CN201510689849 A CN 201510689849A CN 105545285 A CN105545285 A CN 105545285A
- Authority
- CN
- China
- Prior art keywords
- gas
- marine riser
- gas cut
- cut
- void fraction
- 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
- 238000005553 drilling Methods 0.000 title claims abstract description 69
- 238000012544 monitoring process Methods 0.000 title claims abstract description 53
- 239000007788 liquid Substances 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000005514 two-phase flow Effects 0.000 title claims abstract description 16
- 239000011800 void material Substances 0.000 claims description 56
- 230000015572 biosynthetic process Effects 0.000 claims description 37
- 239000012530 fluid Substances 0.000 claims description 30
- 239000000523 sample Substances 0.000 claims description 30
- 239000012071 phase Substances 0.000 claims description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 11
- 230000000704 physical effect Effects 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 6
- 239000007791 liquid phase Substances 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 238000009533 lab test Methods 0.000 claims description 3
- 230000000149 penetrating effect Effects 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 239000013535 sea water Substances 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 2
- 230000006835 compression Effects 0.000 claims description 2
- 238000007906 compression Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 239000000203 mixture Substances 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 3
- 238000011002 quantification Methods 0.000 abstract 2
- 238000001514 detection method Methods 0.000 description 8
- 239000002245 particle Substances 0.000 description 3
- 239000007790 solid phase Substances 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000012806 monitoring device Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000011155 quantitative monitoring Methods 0.000 description 1
- 230000035484 reaction time Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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
- E21B47/00—Survey of boreholes or wells
- E21B47/04—Measuring depth or liquid level
-
- 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
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明属于深水油气钻井领域,具体地,涉及一种基于隔水管气液两相流识别的深水钻井气侵监测方法;在隔水管底部对位安装超声波探头监测气侵信号,根据建立的声波信号衰减与截面含气率的定量对应关系,结合深水钻井井筒气液两相流模型和含可信度地层压力预测方法,建立基于隔水管处气侵监测数据的井筒气侵程度反算方法,提前准确发现气侵的同时,实现对井筒气侵情况的实时定量描述,从而形成基于隔水管气液两相流识别的深水钻井气侵监测方法。本发明的所测数据更能真实、全面反映该截面的含气率情况,可以基于隔水管处的声波监测数据对井底预测原始压力进行修正,使得井筒气侵程度反算结果更接近真实情况,利于采取合理措施控制气侵。
Description
技术领域
本发明属于深水油气钻井领域,具体地,涉及一种基于隔水管气液两相流识别的深水钻井气侵监测方法。
背景技术
钻井过程中若钻遇高压气层会发生气侵,在气侵的初始阶段很难通过泥浆池液面检测、流量差溢流检测等方法或设备监测到,且随着井深的增加,提早发现气侵的难度增大,留给井控的反应时间明显减少,使得气侵的早期检测成为深井井控的重点和难点。而对于深水钻井,首先,随着水深的增加井涌余量减小,更容易发生气侵;其次,在深水高静压环境下,气泡的总体尺度较小,导致气体进入隔水管时泥浆池增量依然不会达到预警值,难以通过常规方法及时检测到气侵。因此,深水钻井中如何及时准确发现气侵是井控工作的重中之重。
近年来,有不少学者针对深水钻井气侵的早期检测进行了研究,包括深水钻井井筒气液两相溢流特征、深水钻井气液两相流温压场计算和深水钻井井筒气液两相流运移特点等基础研究以及泥浆池液面检测法、流量差溢流检测法、随钻声波气侵监测法等深水气侵监测方法的研究。其中声波气侵检测法是利用声波在气液两相中波速的明显差异对气侵进行早期监测的,该方法虽然可以实现气侵的早期监测,但其声波传播历经钻杆和环空内流体,沿程过长导致外界噪声等因素对监测结果产生较大影响,存在较大的误差和不稳定性,且很难分析井底的气侵程度。为了实现更准确的气侵早期监测,有学者提出采用超声波多普勒流量计在隔水管处对气侵进行监测,然而超声波多普勒流量计是基于超声波被流体中的固相颗粒或气泡反射产生频差来近似测量流体速度(认为固相颗粒或气泡与流体速度一致)的,测量结果反映的是固相颗粒和气泡的速度,且只能测量信息窗口区域内的流体速度,测量结果虽可定性检测气侵,但其精度不高,不能全面反映环空内的含气率情况,因此该方法有待进一步改进。同时,深水或超深水钻井对井控有了更高的要求,有必要在及时准确发现气侵的同时尽可能精确的确定井底气侵程度,以便采取合理的措施进行压井,达到安全高效井控的目的。
发明内容
为克服现有技术所存在的缺陷,本发明提供一种基于隔水管气液两相流识别的深水钻井气侵监测方法,可在气泡运移至隔水管底部时及时准确地检测到气侵,同时结合深水钻井井筒多相流模型和含可信度地层压力预测方法可以精确计算气侵时间并实时模拟井底气侵程度和井筒环空内的气体分布,为合理精确的控制气侵提供技术支持。
为实现上述目的,本发明所采用的技术方案如下:
基于隔水管气液两相流识别的深水钻井气侵监测方法,在隔水管底部对位安装超声波探头监测气侵信号,根据建立的声波信号衰减与截面含气率的定量对应关系,结合深水钻井井筒气液两相流模型和含可信度地层压力预测方法,建立基于隔水管处气侵监测数据的井筒气侵程度反算方法,提前准确发现气侵的同时,实现对井筒气侵情况的实时定量描述,从而形成基于隔水管气液两相流识别的深水钻井气侵监测方法。
相对于现有技术,本发明具有如下有益效果:
1、利用超声波在气液两相流中的衰减特性来反映截面含气率,较超声波多普勒测量法有如下优势:超声波穿透截面范围较大,所测数据更能真实、全面反映该截面的含气率情况;可以不考虑气泡对超声波的反射效应,只需通过分析探头接收到的信号即可建立声波信号与截面含气率的定量关系;可以通过多组实验,定量确定一定条件下(流体的类型、物性和流变性,流场和气泡大小等)的探头接收信号与截面含气率的对应关系,从而建立经验公式。
2、通过引入含可信度地层压力预测方法,可以基于隔水管处的声波监测数据对井底预测原始压力进行修正,使得井筒气侵程度反算结果更接近真实情况,利于采取合理措施控制气侵。
3、建立了基于隔水管处气侵监测数据的井筒气侵程度反算方法,在准确检测气侵的同时,还可计算气侵时间、气体到达井口所需的时间、总溢流量以及井筒环空中的气体分布等,实现了井筒气侵情况的实时定量描述,为合理高效井控提供技术支持。
附图说明
图1为隔水管超声波气侵监测示意图;
图2为超声波探头安装位置示意图;
图3为隔水管处截面含气率随时间的变化关系;
图4为气侵后总溢流量随时间的变化关系。
图中:1、隔水管,2、海底井口,3、海底电缆,4、平台监视器,5、超声波发射探头,6、超声波接收探头7、第一延迟块,8、第二延迟块;9、钻杆。
具体实施方式
如图1所示,基于隔水管气液两相流识别的深水钻井气侵监测方法,在隔水管底部对位安装超声波探头监测气侵信号,根据建立的声波信号衰减与截面含气率的定量对应关系,结合深水钻井井筒气液两相流模型和含可信度地层压力预测方法,建立基于隔水管处气侵监测数据的井筒气侵程度反算方法,提前准确发现气侵的同时,实现对井筒气侵情况的实时定量描述,从而形成基于隔水管气液两相流识别的深水钻井气侵监测方法;具体方法如下:
步骤一、获取深水井的参数
深水井的参数,包括:水深h、井深H、已固井段Hg、海水表面温度Tc、隔水管外径尺寸D、循环排量L、钻井液物性密度ρ、钻井液导热系数λl、钻井液稠度系数K、钻井液流性指数n、地层孔隙度Φ、渗透率k、地层导热系数λf、原始地层压力P0。
步骤二、确定实际工况下声波信号与截面含气率的定量表征关系
根据步骤一中深水井的实际参数,采用室内实验模拟该工况下超声波对气侵的响应情况,探头频率从0.2MHz开始以0.2MHz的步长逐步增加到2.0MHz分别进行实验,选取接收声信号幅值最大的探头频率作为最佳探头频率,并建立该频率探头超声波接收信号与截面含气率的定量关系;改变钻井液物性和排量进行实验,建立不同钻井液物性和排量条件的超声波接收信号与截面含气率的关系式,便于实际工况下通过声波监测数据实时计算截面含气率随时间的变化情况。
步骤三、在隔水管底部对位安装超声波探头并实施超声波监测
气侵气体到达隔水管下部时,深水的高静压使气体总体尺度较小,多数情况下,此时通过泥浆池液面很难及时发现并准确判断出气侵的发生,因此在隔水管底部安装气体探测装置可提前发现气侵。
如图1、图2所示,在隔水管1底部分别通过第一延迟块7和第二延迟块8对位安装超声波发射探头5、超声波接收探头6,第一延迟块7、第二延迟块8的外径与钻杆9的外径相切,保证声波传输路径更多的穿过隔水管环空,避免钻杆9对声波的反射,从而最大程度的全面反映隔水管环空的截面含气率情况。
超声波发生器采用900V高电压发生装置,保证不同频率的超声波穿透隔水管壁和气液两相流体,保证声波监测信号的质量。
如图1所示,超声波由超声波发射探头5发射,穿透隔水管1和钻杆9形成的环形空间到达超声波接收探头6,超声波声波监测信号经由海底电缆3实时传输至平台监视器4,平台人员可通过平台监视器4实时了解声波监测数据。
监测装置一旦发现气侵立即报警,同时计算截面含气率随时间的变化情况,便于分析井筒整体的气侵程度。
步骤四、根据声波监测数据实时描述井筒气侵程度
隔水管处一旦检测到气侵,便需要根据气侵情况判断是否需要采取措施进行井控,而判断的依据便是井底气侵时间和井筒气侵程度。
利用步骤二中建立的实际工况下的声波接收信号-截面含气率经验公式,基于步骤三中的隔水管处声波监测数据来计算截面含气率随时间的变化情况;然后利用深水钻井井筒环空气液两相流模型对井筒的气侵程度进行反算,进而确定井筒的气侵程度;具体步骤如下:
(1)、根据目标井的井身结构、钻井液性能、施工参数和地层参数等,建立井筒气液两相流计算模型,控制方程如下:
气相连续性方程:
液相连续性方程:
气、液两相混合运动方程:
气体上升速度经验方程:
vg=C0[vgEg+vlEl]+vrg(4)
气相状态方程:
另外,溢流量计算公式为:QY=A·Vlo(1-Eg)·Δt-Qlo·Δt(6)
总溢流量:
QZY(t)=QZY(t-1)+QZY(t)(7)
其中,A为环空横截面积,m2;Eg,El分别为气相,液相的体积分数,无量纲;ρg,ρgs,ρl分别为气相密度、气相相对密度、钻井液的密度,kg/m3;vg,vl分别为气相、液相的速度,m/s;Qg为单位时间单位厚度气层产出气体的质量,kg/(s·m);Vlo为井口流体流速,m/s;Qlo为泵流量,m3/s;QZY(t)为t时刻的总溢流量,m3;为沿程摩阻压降,Pa/m;为液柱压降,Pa/m;g为重力加速度,m2/s;C0为气相分布系数,无量纲;vrg为气体滑脱速度,m/s;Z为气体压缩因子,无量纲;T为温度,℃;P为压力,Pa。
(2)、引入含可信度地层压力预测方法(柯珂,管志川,周行.深水探井钻前含可信度的地层孔隙压力确立方法[J].中国石油大学学报(自然科学版),2009,33(5):61-67.)对目标井的地层压力进行预测。假设目标井井底原始地层压力为P0,通过计算确定目标井井底含可信度的地层压力取值区间为[P1,PN],对应的可信度为J;
(3)、采用步骤1建立的井筒气液两相流计算模型模拟不同井底原始地层压力条件下的气侵过程,得到一系列不同原始地层压力[P1,PN]条件下的隔水管底部截面含气率随时间的变化情况(变化趋势如图3所示),选取含气率开始非零的时间点[T01,T0N]和截面含气率达到EgM(根据实际需要在0-20%选取)的时间点[Tt1,TtN],计算隔水管处截面含气率(0-EgM)随时间的变化率:
则T0i和Ki即为隔水管处截面含气率反映井底气侵程度的特征参数;
(4)、假设实际工况下声波监测到Δt时间内隔水管底部截面含气率的变化率为Kc,检测到气泡的时刻为Tc,若满足:
|Kc-Kj|<ε(9)
则认为Kj对应的井底原始地层压力Pj为真实的地层压力。其中ε根据实际情况确定取值。
(5)、在步骤2中选取Pj对应的井筒气侵模拟结果,确定气泡到达隔水管底部所用时间,即监测到信号时气侵已经发生了T0j,则气侵发生的时刻为Tc-T0j,同时计算气体到达井口所需时间Tlim;总溢流量QZY根据公式(7)计算确定,从而实现根据隔水管处截面含气率实时确定气侵发生后总溢流量的大小。此外,还可计算气侵发生后每一时刻井筒环空中的气体分布,即不同井深处的截面含气率Eg(H,t)。
由此,根据气侵发生的时刻Tc-T0j、气体到达井口的剩余时间Tlim、任一时刻的总溢流量QZY以及任一时刻不同井深处的截面含气率Eg(H,t),不但可以对气侵进行精确的监测,还可对井筒气侵情况进行实时的定量描述,从而实现了基于隔水管气液两相流识别的深水钻井气侵监测。
为了更加详尽的阐述本发明的实施方式,下面结合具体实施例对本发明进一步进行解释和叙述:
一、获取深水井的参数
深水井相关参数:水深h=1524m,井深H=5030m,已固井段Hg=3153m,海水表面温度Tc=15℃,隔水管外径D=508mm,循环排量L=28L/S,钻井液密度ρ=1.29g/cm3,钻井液导热系数λl=1.73W/(m·℃),稠度系数K=0.27,流性指数n=0.4,地层孔隙度Φ=0.3%,渗透率K=50md,地层导热系数λf=2.25W/(m·℃),预测原始地层压力P0=70MPa。
二、确定实际工况下声波信号与截面含气率的定量表征关系
根据步骤一中深水井的实际参数,采用室内实验模拟该工况下超声波对气侵的响应情况,探头频率从0.2MHz开始以0.2MHz的步长逐步增加到2.0MHz分别进行实验,分析探头接收信号的质量,优选出1MHz为最佳探头频率,并建立该频率探头超声波接收信号与截面含气率的定量关系。改变钻井液物性和排量进行实验,建立不同钻井液物性和排量条件的超声波接收信号与截面含气率的关系式,便于实际工况下通过声波监测数据实时计算截面含气率随时间的变化情况。
三、在隔水管底部对位安装超声波探头并实施超声波监测
采用超声波透射检测方法,利用超声波在气液两相流体中的衰减特性进行气侵的定量监测。如图1、图2所示,在隔水管1底部通过第一延迟块7和第二延迟块8对位安装超声波发射探头5、超声波接收探头6,第一延迟块7、第二延迟块8的外径与钻杆9的外径相切,保证声波传输路径更多的穿过隔水管环空,避免钻杆9对声波的反射,从而最大程度的全面反映隔水管环空的截面含气率情况。
超声波发生器采用900V高电压发生装置,如图1所示,超声波由超声波发射探头5发射,穿透隔水管1和钻杆9形成的环形空间到达超声波接收探头6,超声波声波监测信号经由海底电缆3实时传输至平台监视器4,平台人员可通过平台监视器4实时了解声波监测数据。
监测装置一旦发现气侵立即报警,同时计算截面含气率随时间的变化情况,便于分析井筒整体的气侵程度。
步骤四、根据声波监测数据实时描述井筒气侵程度
具体步骤如下:
1)根据目标井的井身结构、钻井液性能、施工参数和地层参数等,建立井筒气液两相流计算模型;
2)根据地震资料,选取可信度90%,采用含可信度地层压力预测方法对目标井的地层压力进行预测,得到井底压力取值区间为[67.5,72.9];
3)采用步骤1建立的井筒气液两相流计算模型模拟不同井底原始地层压力条件下的气侵过程,得到一系列不同原始地层压力[67.5,72.9]条件下的隔水管底部截面含气率随时间的变化情况(变化趋势如图3所示),选取含气率开始非零的时间点[T01,T0N]和截面含气率达到EgM=10%的时间点[Tt1,TtN],利用公式(8)计算隔水管处截面含气率(0-10%)随时间的变化率:
Ki∈[0.082,0.142]
则T0i和Ki即为隔水管处截面含气率反映井底气侵程度的特征参数。
4)假设气侵发生时的时刻为0点,隔水管底部超声波探头监测到气侵的时刻为Tc,超声波探头监测到隔水管底部截面含气率(0-10%)的变化率Kc=0.09,取ε=0.001,则根据公式(9)可得:Kj=Kc=0.09对应的井底原始地层压力Pj=68.2MPa为真实的地层压力。
5)在步骤2中选取Pj=68.2MPa对应的井筒气侵模拟结果,确定气泡到达隔水管底部所用时间,即监测到信号时气侵已经发生了T0j=Tc=22.3min,同时计算气体到达井口所需时间Tlim=7.6min;总溢流量QZY根据公式(7)计算确定,其随时间的变化情况如图4所示,从而实现了根据隔水管处截面含气率实时确定气侵发生后总溢流量的大小。此外,还可利用步骤1建立的井筒气液两相流计算模型计算气侵发生后每一时刻井筒环空中的气体分布,即不同井深处的截面含气率Eg(H,t)。
由此,根据气侵已发生时间T0j=Tc=22.3min、气体到达井口的剩余时间Tlim=7.6min、任一时刻的总溢流量QZY(见图4)以及任一时刻不同井深处的截面含气率Eg(H,t),不但可以对气侵进行精确的监测,还可对井筒气侵情况进行实时的定量描述,从而实现了基于隔水管气液两相流识别的深水钻井气侵监测。
Claims (7)
1.一种基于隔水管气液两相流识别的深水钻井气侵监测方法,其特征在于:在隔水管底部对位安装超声波探头监测气侵信号,根据建立的声波信号衰减与截面含气率的定量对应关系,结合深水钻井井筒气液两相流模型和含可信度地层压力预测方法,建立基于隔水管处气侵监测数据的井筒气侵程度反算方法,提前准确发现气侵的同时,实现对井筒气侵情况的实时定量描述,从而形成基于隔水管气液两相流识别的深水钻井气侵监测方法。
2.根据权利要求1所述的基于隔水管气液两相流识别的深水钻井气侵监测方法,其特征在于,具体方法如下:
步骤一、获取深水井的参数
步骤二、确定实际工况下声波信号与截面含气率的定量表征关系
步骤三、在隔水管底部对位安装超声波探头并实施超声波监测
步骤四、根据声波监测数据实时描述井筒气侵程度。
3.根据权利要求1-2所述的基于隔水管气液两相流识别的深水钻井气侵监测方法,其特征在于,深水井的参数,包括:水深h、井深H、已固井段Hg、海水表面温度Tc、隔水管外径尺寸D、循环排量L、钻井液物性密度ρ、钻井液导热系数λl、钻井液稠度系数K、钻井液流性指数n、地层孔隙度Φ、渗透率k、地层导热系数λf、原始地层压力P0。
4.根据权利要求1-2所述的基于隔水管气液两相流识别的深水钻井气侵监测方法,其特征在于,步骤二的具体方法如下:根据步骤一中深水井的实际参数,采用室内实验模拟该工况下超声波对气侵的响应情况,探头频率从0.2MHz开始以0.2MHz的步长逐步增加到2.0MHz分别进行实验,选取接收声信号幅值最大的探头频率作为最佳探头频率,并建立该频率探头超声波接收信号与截面含气率的定量关系;改变钻井液物性和排量进行实验,建立不同钻井液物性和排量条件的超声波接收信号与截面含气率的关系式,便于实际工况下通过声波监测数据实时计算截面含气率随时间的变化情况。
5.根据权利要求1-4所述的基于隔水管气液两相流识别的深水钻井气侵监测方法,其特征在于,步骤三的具体方法如下:在隔水管底部分别通过第一延迟块和第二延迟块对位安装超声波发射探头、超声波接收探头,第一延迟块、第二延迟块的外径与钻杆的外径相切,超声波发生器采用900V高电压发生装置;超声波由超声波发射探头发射,穿透隔水管和钻杆形成的环形空间到达超声波接收探头,超声波声波监测信号经由海底电缆实时传输至平台监视器,平台人员可通过平台监视器实时了解声波监测数据。
6.根据权利要求1-5所述的基于隔水管气液两相流识别的深水钻井气侵监测方法,其特征在于,步骤四的具体方法如下:隔水管处一旦检测到气侵,便需要根据气侵情况判断是否需要采取措施进行井控,而判断的依据便是井底气侵时间和井筒气侵程度;利用步骤二中建立的实际工况下的声波接收信号-截面含气率经验公式,基于步骤三中的隔水管处声波监测数据来计算截面含气率随时间的变化情况;然后利用深水钻井井筒环空气液两相流模型对井筒的气侵程度进行反算,进而确定井筒的气侵程度。
7.根据权利要求6所述的基于隔水管气液两相流识别的深水钻井气侵监测方法,其特征在于,具体步骤如下:
(1)、根据目标井的井身结构、钻井液性能、施工参数和地层参数等,建立井筒气液两相流计算模型,控制方程如下:
气相连续性方程:
液相连续性方程:
气、液两相混合运动方程:
气体上升速度经验方程:
vg=C0[vgEg+vlEl]+vrg(4)
气相状态方程:
另外,溢流量计算公式为:QY=A·Vlo(1-Eg)·Δt-Qlo·Δt(6)
总溢流量:
QZY(t)=QZY(t-1)+QZY(t)(7)
其中,A为环空横截面积,m2;Eg,El分别为气相,液相的体积分数,无量纲;ρg,ρgs,ρl分别为气相密度、气相相对密度、钻井液的密度,kg/m3;vg,vl分别为气相、液相的速度,m/s;Qg为单位时间单位厚度气层产出气体的质量,kg/(s·m);Vlo为井口流体流速,m/s;Qlo为泵流量,m3/s;QZY(t)为t时刻的总溢流量,m3;为沿程摩阻压降,Pa/m;为液柱压降,Pa/m;g为重力加速度,m2/s;C0为气相分布系数,无量纲;vrg为气体滑脱速度,m/s;Z为气体压缩因子,无量纲;T为温度,℃;P为压力,Pa;
(2)、引入含可信度地层压力预测方法对目标井的地层压力进行预测,假设目标井井底原始地层压力为P0,通过计算确定目标井井底含可信度的地层压力取值区间为[P1,PN],对应的可信度为J;
(3)、采用步骤一建立的井筒气液两相流计算模型模拟不同井底原始地层压力条件下的气侵过程,得到一系列不同原始地层压力[P1,PN]条件下的隔水管底部截面含气率随时间的变化情况,选取含气率开始非零的时间点[T01,T0N]和截面含气率达到EgM的时间点[Tt1,TtN],,根据实际需要在0-20%选取,计算隔水管处截面含气率0-EgM随时间的变化率:
则T0i和Ki即为隔水管处截面含气率反映井底气侵程度的特征参数;
(4)、假设实际工况下声波监测到Δt时间内隔水管底部截面含气率的变化率为Kc,检测到气泡的时刻为Tc,若满足:
|Kc-Kj|<ε(9)
则认为Kj对应的井底原始地层压力Pj为真实的地层压力;其中ε根据实际情况确定取值;
(5)、在步骤2中选取Pj对应的井筒气侵模拟结果,确定气泡到达隔水管底部所用时间,即监测到信号时气侵已经发生了T0j,则气侵发生的时刻为Tc-T0j,同时计算气体到达井口所需时间Tlim;总溢流量QZY根据公式(7)计算确定,从而实现根据隔水管处截面含气率实时确定气侵发生后总溢流量的大小;此外,还可计算气侵发生后每一时刻井筒环空中的气体分布,即不同井深处的截面含气率Eg(H,t);根据气侵发生的时刻Tc-T0j、气体到达井口的剩余时间Tlim、任一时刻的总溢流量QZY以及任一时刻不同井深处的截面含气率Eg(H,t),不但可以对气侵进行精确的监测,还可对井筒气侵情况进行实时的定量描述,从而实现基于隔水管气液两相流识别的深水钻井气侵监测。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510689849.9A CN105545285B (zh) | 2015-10-22 | 2015-10-22 | 基于隔水管气液两相流识别的深水钻井气侵监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510689849.9A CN105545285B (zh) | 2015-10-22 | 2015-10-22 | 基于隔水管气液两相流识别的深水钻井气侵监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105545285A true CN105545285A (zh) | 2016-05-04 |
CN105545285B CN105545285B (zh) | 2018-08-10 |
Family
ID=55824790
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510689849.9A Active CN105545285B (zh) | 2015-10-22 | 2015-10-22 | 基于隔水管气液两相流识别的深水钻井气侵监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105545285B (zh) |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105956566A (zh) * | 2016-05-09 | 2016-09-21 | 中国海洋石油总公司 | 一种深海钻井气侵检测信号处理方法及处理电路 |
CN107060737A (zh) * | 2017-05-26 | 2017-08-18 | 中国石油天然气集团公司 | 一种随钻气侵模拟实验装置及实验方法 |
CN107060731A (zh) * | 2017-07-04 | 2017-08-18 | 中国石油集团钻井工程技术研究院 | 一种基于井涌余量的深水钻井套管下深修正方法 |
CN107389787A (zh) * | 2017-07-28 | 2017-11-24 | 电子科技大学 | 一种基于频域的气侵检测信号处理方法 |
CN107575212A (zh) * | 2017-10-18 | 2018-01-12 | 中国石油大学(北京) | 超声波随钻气侵监测装置及方法 |
CN108071384A (zh) * | 2017-12-07 | 2018-05-25 | 中国石油大学(北京) | 隔水管气侵监测系统及方法 |
CN109138988A (zh) * | 2018-11-01 | 2019-01-04 | 中国石油大学(华东) | 隔水管超声波气侵监测短节及监测方法 |
CN109164205A (zh) * | 2018-07-06 | 2019-01-08 | 覃楚倩 | 一种基于海底基盘的钻探钻井气体监测系统及其监测方法 |
CN109458172A (zh) * | 2018-11-01 | 2019-03-12 | 中国石油大学(华东) | 隔水管气侵监测工具及监测方法 |
CN110185433A (zh) * | 2019-05-16 | 2019-08-30 | 中国海洋石油集团有限公司 | 一种基于谱特征分析方法的隔水管气侵监测装置及方法 |
CN110312460A (zh) * | 2017-02-22 | 2019-10-08 | 奥林匹斯冬季和Ibe有限公司 | 准备设备和用于清洁和/或消毒医疗器械的准备设备的操作方法 |
CN111119864A (zh) * | 2018-11-01 | 2020-05-08 | 中国石油化工股份有限公司 | 基于气侵压力特征的溢流监测方法及系统 |
CN111364979A (zh) * | 2020-03-23 | 2020-07-03 | 中国石油大学(华东) | 一种基于超声波的井下气侵监测系统 |
CN111997586A (zh) * | 2020-08-07 | 2020-11-27 | 长江大学 | 一种利用阵列流量测井获取水平井持水率的方法 |
CN112483074A (zh) * | 2019-09-11 | 2021-03-12 | 中国石油化工股份有限公司 | 一种用于预测钻井过程气侵现象的方法及系统 |
CN112855122A (zh) * | 2020-12-31 | 2021-05-28 | 中国石油大学(华东) | 一种井下气液固三相流超声波气侵监测系统及实施方法 |
CN113006769A (zh) * | 2021-03-17 | 2021-06-22 | 中国石油大学(华东) | 一种复杂压力体系地层智能压井方法及装置 |
CN113323653A (zh) * | 2021-06-15 | 2021-08-31 | 中海油研究总院有限责任公司 | 一种深水钻井溢流早期预警方法和装置 |
CN113530524A (zh) * | 2021-07-12 | 2021-10-22 | 中国石油大学(华东) | 一种井筒流量监测系统及流量、含水率解释方法 |
CN113806919A (zh) * | 2021-08-24 | 2021-12-17 | 中国石油大学(华东) | 一种基于隔水管外含气率监测的深水地层参数预测方法 |
CN114109357A (zh) * | 2021-12-17 | 2022-03-01 | 中国石油大学(北京) | 深水气侵模拟实验装置及气侵判断方法 |
CN115060796A (zh) * | 2022-06-24 | 2022-09-16 | 西南石油大学 | 满管多相钻井液流速与含气量超声联合检测方法 |
CN114109357B (zh) * | 2021-12-17 | 2024-05-24 | 中国石油大学(北京) | 深水气侵模拟实验装置及气侵判断方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6993963B1 (en) * | 2000-09-22 | 2006-02-07 | Jon Steinar Gudmundsson | Method for determining pressure profiles in wellbores, flowlines and pipelines, and use of such method |
CN202788815U (zh) * | 2012-07-06 | 2013-03-13 | 中国石油天然气股份有限公司 | 油气侵监测报警系统 |
CN103032064A (zh) * | 2013-01-11 | 2013-04-10 | 西南石油大学 | 一种钻井过程中气侵位置检测的方法及装置 |
CN104153759A (zh) * | 2014-07-30 | 2014-11-19 | 中国石油集团钻井工程技术研究院 | 控压钻井气液两相流动模拟计算方法 |
CN104712320A (zh) * | 2015-01-29 | 2015-06-17 | 中国石油大学(华东) | 钻井过程中气侵早期监测装置及方法 |
-
2015
- 2015-10-22 CN CN201510689849.9A patent/CN105545285B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6993963B1 (en) * | 2000-09-22 | 2006-02-07 | Jon Steinar Gudmundsson | Method for determining pressure profiles in wellbores, flowlines and pipelines, and use of such method |
CN202788815U (zh) * | 2012-07-06 | 2013-03-13 | 中国石油天然气股份有限公司 | 油气侵监测报警系统 |
CN103032064A (zh) * | 2013-01-11 | 2013-04-10 | 西南石油大学 | 一种钻井过程中气侵位置检测的方法及装置 |
CN104153759A (zh) * | 2014-07-30 | 2014-11-19 | 中国石油集团钻井工程技术研究院 | 控压钻井气液两相流动模拟计算方法 |
CN104712320A (zh) * | 2015-01-29 | 2015-06-17 | 中国石油大学(华东) | 钻井过程中气侵早期监测装置及方法 |
Non-Patent Citations (7)
Title |
---|
任美鹏等: "深水钻井井筒气液两相溢流特征及其识别方法", 《工程热物理学报》 * |
任美鹏等: "钻井气液两相流体溢流与分布特征研究", 《工程热物理学报》 * |
朱红钧等: "气侵期间垂直环空气液两相流动模拟分析", 《油气田地面工程》 * |
林英松等: "垂直环空气液两相流流型的实验研究", 《石油大学学报(自然科学版)》 * |
柯珂等: "深水探井钻前含可信度的地层孔隙压力确立方法", 《中国石油大学学报(自然科学版)》 * |
陈平等: "深水钻井溢流早期监测技术研究现状", 《石油学报》 * |
隋秀香等: "基于多普勒测量技术的深水隔水管气侵早期监测研究", 《石油钻探技术》 * |
Cited By (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105956566A (zh) * | 2016-05-09 | 2016-09-21 | 中国海洋石油总公司 | 一种深海钻井气侵检测信号处理方法及处理电路 |
CN105956566B (zh) * | 2016-05-09 | 2021-10-15 | 中国海洋石油集团有限公司 | 一种深海钻井气侵检测信号处理方法及处理电路 |
CN110312460B (zh) * | 2017-02-22 | 2022-01-28 | 奥林匹斯冬季和Ibe有限公司 | 准备设备和用于清洁和/或消毒医疗器械的准备设备的操作方法 |
CN110312460A (zh) * | 2017-02-22 | 2019-10-08 | 奥林匹斯冬季和Ibe有限公司 | 准备设备和用于清洁和/或消毒医疗器械的准备设备的操作方法 |
CN107060737B (zh) * | 2017-05-26 | 2020-09-08 | 中国石油天然气集团公司 | 一种随钻气侵模拟实验装置及实验方法 |
CN107060737A (zh) * | 2017-05-26 | 2017-08-18 | 中国石油天然气集团公司 | 一种随钻气侵模拟实验装置及实验方法 |
CN107060731A (zh) * | 2017-07-04 | 2017-08-18 | 中国石油集团钻井工程技术研究院 | 一种基于井涌余量的深水钻井套管下深修正方法 |
CN107389787A (zh) * | 2017-07-28 | 2017-11-24 | 电子科技大学 | 一种基于频域的气侵检测信号处理方法 |
CN107389787B (zh) * | 2017-07-28 | 2019-09-24 | 电子科技大学 | 一种基于频域的气侵检测信号处理方法 |
CN107575212A (zh) * | 2017-10-18 | 2018-01-12 | 中国石油大学(北京) | 超声波随钻气侵监测装置及方法 |
CN107575212B (zh) * | 2017-10-18 | 2020-10-13 | 中国石油大学(北京) | 超声波随钻气侵监测装置及方法 |
CN108071384A (zh) * | 2017-12-07 | 2018-05-25 | 中国石油大学(北京) | 隔水管气侵监测系统及方法 |
CN109164205A (zh) * | 2018-07-06 | 2019-01-08 | 覃楚倩 | 一种基于海底基盘的钻探钻井气体监测系统及其监测方法 |
CN109458172A (zh) * | 2018-11-01 | 2019-03-12 | 中国石油大学(华东) | 隔水管气侵监测工具及监测方法 |
CN111119864A (zh) * | 2018-11-01 | 2020-05-08 | 中国石油化工股份有限公司 | 基于气侵压力特征的溢流监测方法及系统 |
CN109138988A (zh) * | 2018-11-01 | 2019-01-04 | 中国石油大学(华东) | 隔水管超声波气侵监测短节及监测方法 |
CN110185433A (zh) * | 2019-05-16 | 2019-08-30 | 中国海洋石油集团有限公司 | 一种基于谱特征分析方法的隔水管气侵监测装置及方法 |
CN112483074A (zh) * | 2019-09-11 | 2021-03-12 | 中国石油化工股份有限公司 | 一种用于预测钻井过程气侵现象的方法及系统 |
CN111364979B (zh) * | 2020-03-23 | 2023-05-23 | 中国石油大学(华东) | 一种基于超声波的井下气侵监测系统 |
CN111364979A (zh) * | 2020-03-23 | 2020-07-03 | 中国石油大学(华东) | 一种基于超声波的井下气侵监测系统 |
CN111997586A (zh) * | 2020-08-07 | 2020-11-27 | 长江大学 | 一种利用阵列流量测井获取水平井持水率的方法 |
CN111997586B (zh) * | 2020-08-07 | 2024-03-26 | 长江大学 | 一种利用阵列流量测井获取水平井持水率的方法 |
CN112855122A (zh) * | 2020-12-31 | 2021-05-28 | 中国石油大学(华东) | 一种井下气液固三相流超声波气侵监测系统及实施方法 |
CN113006769A (zh) * | 2021-03-17 | 2021-06-22 | 中国石油大学(华东) | 一种复杂压力体系地层智能压井方法及装置 |
CN113323653A (zh) * | 2021-06-15 | 2021-08-31 | 中海油研究总院有限责任公司 | 一种深水钻井溢流早期预警方法和装置 |
CN113323653B (zh) * | 2021-06-15 | 2024-05-03 | 中海油研究总院有限责任公司 | 一种深水钻井溢流早期预警方法和装置 |
CN113530524A (zh) * | 2021-07-12 | 2021-10-22 | 中国石油大学(华东) | 一种井筒流量监测系统及流量、含水率解释方法 |
CN113530524B (zh) * | 2021-07-12 | 2023-05-26 | 中国石油大学(华东) | 一种井筒流量监测系统及流量、含水率解释方法 |
CN113806919B (zh) * | 2021-08-24 | 2023-09-22 | 中国石油大学(华东) | 一种基于隔水管外含气率监测的深水地层参数预测方法 |
CN113806919A (zh) * | 2021-08-24 | 2021-12-17 | 中国石油大学(华东) | 一种基于隔水管外含气率监测的深水地层参数预测方法 |
CN114109357A (zh) * | 2021-12-17 | 2022-03-01 | 中国石油大学(北京) | 深水气侵模拟实验装置及气侵判断方法 |
CN114109357B (zh) * | 2021-12-17 | 2024-05-24 | 中国石油大学(北京) | 深水气侵模拟实验装置及气侵判断方法 |
CN115060796A (zh) * | 2022-06-24 | 2022-09-16 | 西南石油大学 | 满管多相钻井液流速与含气量超声联合检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105545285B (zh) | 2018-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105545285A (zh) | 基于隔水管气液两相流识别的深水钻井气侵监测方法 | |
RU2607826C2 (ru) | Скважинный инструмент для определения скорости потока | |
RU2374443C2 (ru) | Система оповещения о выбросе, использующая высокочастотный режим флюида в стволе скважины | |
US5561245A (en) | Method for determining flow regime in multiphase fluid flow in a wellbore | |
CN109386279B (zh) | 一种井筒气侵检测方法及系统 | |
CN103032064A (zh) | 一种钻井过程中气侵位置检测的方法及装置 | |
CN110344815B (zh) | 一种基于分布式光纤声音监测和分布式光纤温度监测的生产剖面监测方法 | |
US20190085686A1 (en) | Integrated logging tool method for identifying well damage | |
US20150338378A1 (en) | Ultrasonic signal time-frequency decomposition for borehole evaluation or pipeline inspection | |
CN101268251A (zh) | 声学流体分析仪 | |
Fu et al. | Development and testing of kick detection system at mud line in deepwater drilling | |
US20160291186A1 (en) | Seismic cable, system and method for acquiring information about seismic, microseismic and mechanical vibration incidents in a well | |
AU2016254907A1 (en) | A method of identifying a material and/or condition of a material in a borehole | |
Yuqiang et al. | The quantitative description of gas-cut degree in deepwater drilling | |
CN103926422A (zh) | 流体测量系统和方法 | |
CN105822264B (zh) | 天然气水合物层钻井井筒水合物动态分解位置检测方法 | |
Ünalmis | Sound speed in downhole flow measurement | |
CN103884625A (zh) | 一种基于数据融合的水库底泥密度探测方法和系统 | |
NL2022689B1 (en) | Method and system for non-intrusively determining cross-sectional variation for a fluidic channel | |
CN107191179A (zh) | 一种油气井动液面测试方法 | |
Tabjula et al. | Empirical correlations for predicting flow rates using distributed acoustic sensor measurements, validated with wellbore and flow loop data sets | |
Jia et al. | A review of dynamic fluid level detection for oil well | |
Shinyakov et al. | Analysis of methods for measuring the liquid level in the annular space of an oil well | |
WO2016086794A1 (zh) | 用于确定钻杆位置的系统和方法 | |
Wang et al. | Location estimation of subsurface fluid-filled fractures: Cepstral predominant peak analysis and numerical study |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |