摘要：在工程实际应用中发现，桥塞坐封后，由于地层因素和人为操作，上下端流体的温度和压力会产生变化，使桥塞的坐封失效。其中，二级桥塞受影响失效比例最大。因此，二级桥塞坐封后的段塞区内流体状态参数的获取及相关影响因素的分析对工程应用十分重要，然而国内外学者对注采井筒温度压力场研究颇多，还未涉及针对桥塞坐封方面的温度压力场变化的研究。因此，基于瞬态温度场模型，结合流体状态方程，运用MATLAB软件编程，计算得段塞区温度和压力随时间及井深变化的关系。针对井下5 km处二级桥塞进行计算，设初始温度60℃，结果显示：(1)随着时间推移，井筒回温使段塞区温度升高，压力随之升高，20 h后变化趋势停止，即当井筒温度与地层相等时，温度和压力上升趋势停止；(2)进入稳态后，段塞区深度对其压力及温度有影响，且温度压力随深度的增大而呈正比例增加；(3)进入稳态后，段塞区间隔距离对段塞区温度及压力无影响。计算结果可为工程实际应用提供一定参考价值。

　　关键词：瞬态；温度压力场；桥塞段塞区；二级桥塞

　　0引言

　　在中浅层试油施工中出现的干层、水层、气层及异常高压等特殊层位[1]，为方便后续试油，封堵废弃层位，通常采用桥塞进行封层[2-3]。同时，对于部分短期无开发计划的试油结束井也采用桥塞封井[4-6]。桥塞坐封时通常采用电缆或管柱将其输送到井筒预定位置，通过火药爆破、液压坐封或者机械坐封工具产生的压力作用于上卡瓦，拉力作用于张力棒，通过上下锥体对密封胶筒施以上压下拉两个力，当拉力达到一定值时，张力棒断裂，坐封工具与桥塞脱离。此时桥塞中心管上的锁紧装置发挥效能，上下卡瓦破碎并镶嵌在套管内壁上，胶筒膨胀并密封，完成坐封。此外，桥塞也用于深层气井的已经试井层封堵，为上返测试、压裂改造等工艺技术的成功实施提供保障[7]。但在高温高压深井中，由于地层压力较大，如图1所示的单级桥塞额定压差往往无法满足正常施工需求。因此，在工程实际中，通常采用双级桥塞利用逐级降压的方式使桥塞的工作压差满足工作条件。国内外学者对于非稳态注采问题[8-10]和桥塞结构[11-12]的研究关注度非常高，但对于温度场压力场理论计算应用于桥塞坐封效果及影响机理方面几乎没有涉及。但是，双级桥塞井下试油作业工况复杂，桥塞工作状态难以明确，亟需通过温度场、压力场相关理论探明坐封位置[13-14]、坐封间距等因素对桥塞坐封后上下端面压力的影响规律，为后续作业提供理论指导。

　　本文针对二级桥塞坐封期间，建立瞬态温度压力场模型。进行以下计算模拟研究：(1)对井下5 000 m深度处，二级桥塞的段塞区温度压力变化进行研究；(2)利用有限元分析软件对井筒流体进行模拟仿真分析，验证理论的正确性；(3)探讨进入稳态后，段塞区深度对段塞区温度及压力的影响规律；(4)探讨进入稳态后，段塞区间隔距离对段塞区温度及压力的影响规律。

　　1瞬态温度压力场模型创建

　　二级桥塞坐封后，在一级桥塞与二级桥塞之间会形成一段封闭的空间——段塞区。段塞区内流体在地层温度的作用下温度升高，流体体积受热膨胀。同时，段塞区部位套管受注入流体的影响，温度下降，发生径向收缩。在以上影响的作用下，必然导致段塞区压力有所变化，应当采用温度场理论对其进行分析。

　　已知，桥塞坐封后，由于地层与井筒内部存在较大温差。由热力学第二定律，井筒内温度会按照一定的规律上升至与地层温度相同。这一过程为非稳态导热过程，因此，根据非稳态热传导规律[15-16]，导热微分方程及其定解条件为：

　　式中：t为温度，℃；τ为时间常数，s；a为热扩散率，a=λ/ρcp，m2/s；x为径向位置，m；h为表面传热系数，W/(m2·K)；λ为导热系数，W/(m·K)；R为井筒半径，m。

　　2二级桥塞段塞区温度和压力变化规律

　　针对瞬态计算模型，对井筒状态参数做出合理假设，并进行MATLAB编程，获得二级桥塞段塞区温度和压力变化规律。

　　2.1时间对二级桥塞段塞区温度和压力影响规律分析

　　由式(16)可知，段塞区压力与温度按照一定规律增加。

　　另，针对井筒工况，建立了瞬态计算模型，对井筒状态参数做出合理假设，并进行MATLAB编程如下：

　　同时，使用有限元分析软件对井筒流体进行模拟仿真分析，以验证计算结果。首先，建立井筒流体模型，定义其内部均为液态水，初始温度为60℃。其次，对模型划分网格，设置网格大小为5 mm，如图3所示。然后，设置上下端面条件为绝热，设置侧面对流换热系数为10-3 W/(mm2⋅℃)，设置环境温度为170℃。最后，设置10个子步，每个子步10 000 s，总计时间为100 000 s。图4所示为初始时刻与达到稳态时热通量示意图，其结果与编程计算结果一致，故本文所得结论可靠。

　　由上述结果可知：5 000 m深处，在桥塞坐封后，由于地温的影响，段塞区的温度回升，压力也随之升高。
　　2.2进入稳态后段塞区深度对段塞区温度及压力的影响规律分析

　　计算结果如图5所示，稳态后温度与地层温度一致，即温度随地层深度线性增加。由式(16)可知，在其他条件不变的情况下，压力也随地层深度呈线性增加趋势。

　　上述结果表明：在保持初始条件不变的情况下，随着桥塞下井深度增加，桥塞段的稳态时刻温度呈上升趋势。由于桥塞的段塞区坐封时体积大小被固定，温度的大幅提高势必会引起段塞区压力的明显升高。且初始时刻与稳态时的温差越大，压力升高就越明显。所以，可以根据这一结果给出合理桥塞应用的工程应用建议：建议在深层使用桥塞分层工艺时，需严格考虑其受温度场压力场变化的影响程度，针对不同情况选取合适的桥塞。
　　2.3进入稳态后段塞区间隔距离对段塞区温度及压力的影响规律分析
　　段塞区间隔距离变化，即体积变化。针对段塞区，已知状态方程(16)，由：

　　式中：P为压力，Pa；ρ为密度，kg/m3；M为摩尔质量，kg/mol；R为常数，8.31 J/(mol·K)；T为温度，K。

　　由式(18)可知：当达到稳态时，二级封隔器下端温度及压力变化与段塞区体积无关。所以，段塞区间隔对二级封隔器下端温度及压力无影响。

　　3结束语

　　针对工程实际应用中存在井筒内桥塞受温度压力影响失效的问题，本文考虑二级桥塞段塞区工作期间的实际情况，基于瞬态温度场压力场理论计算，建立分析模型，应用有限元软件模拟验证，对温度压力场变化对坐封效果的影响进行了探讨。经过对计算结果对比分析，得出以下结论。

　　(1)二级桥塞段塞区的温度对其压力变化有显著影响。随着时间推移，井筒回温使段塞区温度升高，压力随之升高。当井筒温度与地层相等时，温度和压力上升趋势停止。对于井深5 000 m处桥塞，20 h后达到稳态，温度升高与地层温度和初始温度的温差有关。

　　(2)进入稳态后，段塞区深度对其压力及温度有影响，且温度压力随深度的增大而呈正比例增加。因此，建议在深层使用桥塞分层工艺时，需严格考虑其受温度场压力场变化的影响程度，针对不同情况选取合适的桥塞。

　　(3)进入稳态后，考虑气体状态方程，得出结论为段塞区间隔距离对段塞区温度及压力无影响，故在实际工程应用中，可不考虑桥塞间隔对桥塞坐封效果的影响。

　　参考文献：

　　[1]马先林,周德胜,蔡文斌,等.基于可解释机器学习的水平井产能预测方法[J].西南石油大学学报(自然科学版),2022,44(4):81-90.

　　[2]费世祥,余浩杰,陈存良,等.致密砂岩气藏水平井开发关键技术——以苏里格气田上古生界为例[J].西安石油大学学报(自然科学版),2022,37(4):26-35.

　　[3]吴泽兵,席凯凯,王杰,等.水平井仿生磨鞋钻磨复合桥塞的设计与仿真[J].石油机械,2022,50(2):8-14.

　　[4]陈宁,邓凯,陈小荣,等.可固化堵漏技术在长庆油田的研究与应用[J].中国石油和化工标准与质量,2021,41(3):157-159.

　　[5]张凯权,刘飞,朱永国,等.非爆炸膨胀式过油管桥塞封堵技术[J].油气井测试,2020,29(3):32-37.

　　[6]宋海.临时封堵工艺在高压油气井中的应用[J].化工管理,2020(15):183-184.

　　[7]郑鑫,陈春明,贾德友.水平井桥塞分段压裂技术在Toson-Uul油田的应用评价[J].化学工程与装备,2022(4):90-92.

　　[8]WANG J,LI J,LIU G,et al.Development and application of well‐bore heat transfer model considering variable mass flow[J].Un‐derground Space,2021,6(3):316-328.

　　[9]LI X,ZHANG J,TANG X,et al.Study on wellbore temperature of riserless mud recovery system by CFD approach and numeri‐cal calculation[J].Petroleum,2020,6(2):163-169.

　　[10]WANG J,SUN B,CHEN W,et al.Calculation model of un‐steady temperature-pressure fields in wellbores and fractures of supercritical CO2 fracturing[J].Fuel,2019,253(1):1168-1183.

　　[11]平恩顺,王林,解同川,等.套管压裂井不压井完井工艺技术研究及应用[J].石油机械,2022,50(5):113-118.

　　[12]张帅,苗典远,江民盛,等.液氮冷冻暂堵技术在海上井控处置的可行性分析[J].钻采工艺,2022,45(2):28-32.

　　[13]ZHENG J,LI Z,DOU Y,et al.Study on Downhole Throttling Characteristics of High Water Content Gas[J].Geofluids,2022(1):1-14.

　　[14]窦益华,缑雅洁,郑杰,等.高温高压气井温度压力耦合分布研究[J].机械设计与制造工程,2022,51(3):103-107.

　　[15]冯剑,许博越,付建红,等.深水压井井筒瞬态传热及复杂流动行为研究[J].石油机械,2021,49(2):1-6.

　　[16]杨世铭,陶文铨.传热学(第五版)[M].北京:高等教育出版社,2019.

　　[17]Campo A.Rapid determination of spatio-temporal temperature an heat transfer in simple bodies cooled by convection:usage of caculators in lieu of Heisler-Grobercharts[J].IntComm Heat Mass Transfer,1997,24(4):553-564.