1、第30卷第1期合川气田井下节流气井动态分析新方法徐建军1，李云杰1，杨香艳2，杜怡鹤3（1.中国石油大庆油田有限责任公司重庆分公司，重庆402660；2.重庆工商大学外国语学院，重庆400067；3.西南石油大学石油与天然气工程学院，四川 成都610500）摘要合川气田属于我国典型的低渗透、低产气田，为提高开发效益，合川气田气井绝大多数已应用了井下节流工艺。但是，井下节流阻断了井口油压与井底流压的自然联系，给合川气田气井生产动态分析带来了极大挑战。为实现对井下节流气井井口动态监测资料的有效分析，文中以单井开采系统为研究对象，根据气藏、气水两相井筒管流和油嘴节流理论模型及动态分析特征化方法，建立
2、了新型井下节流气井系统综合生产动态分析方法。该方法将传统动态特征化分析法扩展到井下节流气井，并通过对合川气田实际井下节流气井生产数据的拟合分析，验证了新方法的可靠性和实用性。关键词低渗透；气水两相；井下节流；井筒管流；动态分析中图分类号：TE33+2文献标志码：A收稿日期：20220506；改回日期：20221111。第一作者：徐建军，男，1975年生，工程硕士，从事天然气开发管理工作。E-mail：。引用格式：徐建军，李云杰，杨香艳，等.合川气田井下节流气井动态分析新方法J.断块油气田，2023，30（1）：143148.XU Jianjun，LI Yunjie，YANG Xiangyan，
3、et al.A new method for performance analysis of downhole throttling gas wells in Hechuan gas fieldJ.Fault-Block Oil&Gas Field，2023，30（1）：143148.A new method for performance analysis of downhole throttling gas wells in Hechuan gas fieldXU Jianjun1,LI Yunjie1,YANG Xiangyan2,DU Yihe3(1.Chongqing Branch
4、of Daqing Oilfield Co.,Ltd.,PetroChina,Chongqing 402660,China;2.School of Foreign Languages,Chongqing Industrial and Commercial University,Chongqing 400067,China;3.Petroleum Engineering School,Southwest PetroleumUniversity,Chengdu 610500,China)Abstract:Hechuanisatypicallowpermeabilityandlowproductio
5、ngasfieldinChina.Inordertoimprovethedevelopmentefficiencyofthe gas field,most gas wells in Hechuan block have applied downhole throttling technology.However,the natural connection betweenwellhead oil pressure and bottom hole flow pressure is blocked by downhole throttling,which brings great challeng
6、es to the productionperformance analysis of gas wells in Hechuan block.In order to realize the analysis of wellhead dynamic monitoring data of downholethrottling gas wells,taking the single well production system as the research object,and combining the theoretical model and dynamicanalysis characte
7、rization method of gas reservoir,gas-water two-phase wellbore pipe flow and choke throttling,a new comprehensiveproduction dynamic analysis method of downhole throttling gas well system is established.This method extends the traditional dynamiccharacteristic analysis method into downhole throttling
8、wells,and verifies the reliability and practicability of new method by fitting andanalyzingtheproductiondataofactualdownholethrottlinggaswellsinHechuanblock.Key words:low permeability;gas-water two-phase;downhole throttling;wellbore pipe flow;performance analysis0引言合川气田属于我国典型的低渗透、低产气田，为降低开发成本，提高开发效益
9、，合川气田气井绝大多数应用了井下节流工艺15。通过井下节流，可以减少节流器上方系统的压力，并降低水合物形成温度，避免井口、井筒和地面管道中发生水合物堵塞，从而改善节流后段流速和携液量。但是，井下节流阻断了井口油压与井底流压的自然联系，这给合川气田气井生产动态分析带来了极大挑战67。井下节流工艺在低渗透气田中已应用近20 a，并取得了良好的应用效果。针对气井井下节流动态分析，目前已经开展了大量的研究。王宇等8应用节点系统分析方法，以井下节流器为节点，建立了气井井下节流过程井筒压力、温度动态预测模型，该模型充分考虑了流入动态和井筒的径向换热机制，可对井下节流压力、温度等重要参数进行预测。孙蓉等9针
10、对须家河气藏井下节流生产井的特点，开发出适合非酸性气井使用断块油气田FAULT-BLOCK OIL GAS FIELDdoi:10.6056/dkyqt2023010202023年1月断块油气田2023年1月的直井型电子脱挂装置，有效地解决了井下节流气井动态监控的技术难题。刘顺等10推导了考虑节流效应影响的压力公式，根据节流器上流油压和套管压力的变化规律，先计算出积液深度，然后从套管环空出发，采用井深迭代法，得到了计算井底压力的新方法。齐彩婷11通过节点系统分析方法，利用PIPESIM软件分析节流前后的井筒压力、温度分布，同时对节流下入深度、节流嘴尺寸、气液比相关节流参数进行敏感性分析，进而分
11、析不同情况排液的效果。虽然近年来国内外许多学者在油嘴节流模型、井下节流温度、压力场预测模型等方面开展了大量理论研究，但鲜有考虑气藏井筒油嘴节流耦合的产量动态分析方法，目前采用传统动态分析方法对井下节流气井进行分析存在很大缺陷。为此，本文以单井开采系统为研究对象，结合气藏、井筒和油嘴节流理论模型，建立了井下节流气井系统综合生产动态分析新方法，并通过对实际生产数据拟合，分析验证了该新方法的可靠性和实用性。1气藏井筒油嘴耦合动态分析模型1.1井下节流工艺原理气井的井下节流技术，是将井下节流装置坐放在井筒内的适当位置来达到冷却、降低压力的目的，从而使得节流部位的地层温度高于水合物生成的临界温度，防止井
12、筒中水合物的生成，进而减少地面节流及保温措施的一种采气技术。井下节流嘴相当于截面突缩装置，天然气流经节流嘴时流速会迅速增加，压力明显下降（见图1。图中：d为油嘴直径，T1，T2（T1T2）分别为流体通过油嘴前、后的温度，p1，p2（p1p2）分别为流体通过油嘴前、后的压力，红色箭头代表气体流动方向）。图1气井井下节流示意Fig.1Downhole throttling diagram of gas well1.2气藏模型为提高气井产能，合川气田大部分井采用压裂改造工艺。为此，以Ozkan&Raghavan点源函数法理论为基础，应用瞬时源函数理论和叠加原理，得到圆形封闭无限导流压裂井Laplac
13、e空间无因次压力响应解12：軍D軍 軍p=12s1-1乙K0rDs姨D軍+I0rDs姨D軍K1reDs姨D軍I1reDs姨D軍D軍d（1）式中：軍D为Laplace空间无因次压力；p为压力，Pa；s为Laplace空间无因次时间；K0，K1分别为修正0阶、1阶二类贝塞尔函数；rD为无因次半径；I0，I1分别为修正0阶、1阶一类贝塞尔函数；reD为无因次井控半径；为无因次距离。在无限导流压裂井压力响应解获得的基础上，Pratikno等13提出了定产圆形封闭条件有限导流压裂井无因次压力响应解：軍wD-12kj=1q軈fDD軍sxDj+1xDj乙K0 xD-xD軍+K0 xD+xD軍D軍dx+12s
14、1fD軍s s姨K1reDs姨D軍I1reDs姨D軍s姨1-xDD 軍0乙I0D 軍 d+s姨1+xDD 軍0乙I0D 軍 dD軍-sfq軈fDxD，D軍s=sFcDxD-sxD0乙s0乙q軈fDu，D軍s dudD軍s（2）式中：軍wD为Laplace空间无因次井底压力；q軈fD为Laplace无因次裂缝流量；FcD为无因次裂缝导流能力；sf为压裂井表皮系数；xD为无因次裂缝位置坐标；x为无因次离散段位置；xDj，xDj+1分别为裂缝离散段上、下位置无因次坐标；u为积分函数中的无因次位置；j为离散段个数；k为离散段时间步数。根据杜哈美叠加原理，变流动压力条件下的流量响应等于定流压的流量响应与
15、变流动压力的褶积：qDtDD 軍=tD0乙qDtD-D軍鄣D鄣D軍d（3）式中：qD为无因次流量；tD为无因次时间；为积分函数中的无因次时间；D为无因次压力。通过式（3）可以实现井底流压与产量之间的转化，结合现代产量递减特征化理论1415，可获得产量动态分析特征化图版（见图2。图中：qDd表示无因次规整化产量，qDdi表示无因次规整化产量积分，qDdid表示无因次规整化产量积分导数，tDd表示无因次物质平衡拟时间）。1.3井筒管流模型HagedornBrown两相管流模型在计算直井井筒压144第30卷第1期其中：=fg+fgCvg+fwCvwNNMZR（8）1=gfww（9）式中：为相对气体质
16、量分数；pr为临界压力比；N为气体多变膨胀指数；1为气体相对密度比；A1，A2分别为油嘴上游油管横截面积、油嘴喉部横截面积，m2；fg，fw分别为气、水的质量分数；Cvg，Cvw分别为气、水的定容比热容，J/(kgK)；M为气体分子摩尔质量，g/mol；Z为气体偏差因子；R为理想气体常数；g，w分别为混合物气、水相密度，kg/m3。根据上述气水相节流理论模型，建立了气水两相节流后临界压力/节流前压力（pc/p1）与相对气体质量分数的关系（见图3）。从图3可看出，当气井生产过程中出现气水两相流动后，临界压力/节流前压力明显下降。这也说明了井下节流油、套压差值大主要是气水两相节流导致的。图3pc/p1与的关系Fig.3The relationship between pc/p1and 1.5气藏井筒油嘴耦合分析模型基于上述气藏、井筒、油嘴压力计算模型，采用迭代法求解气井流动耦合系统。计算步骤（步骤中t为当前时间，n为时间步，i为迭代次数，为收敛精度）为：1）由给定的油嘴直径d及油嘴的下游压力pchd，求出最小的上游压力pwh，min，取piwh，n=pwh，min，i=0。力场方面具有准确
