CN108121894A - 一种垂直气井井筒多液滴夹带临界气流量计算方法 - Google Patents
一种垂直气井井筒多液滴夹带临界气流量计算方法 Download PDFInfo
- Publication number
- CN108121894A CN108121894A CN201711387764.0A CN201711387764A CN108121894A CN 108121894 A CN108121894 A CN 108121894A CN 201711387764 A CN201711387764 A CN 201711387764A CN 108121894 A CN108121894 A CN 108121894A
- Authority
- CN
- China
- Prior art keywords
- liquid
- critical
- gas
- drop
- droplet
- 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
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Abstract
本发明涉及一种考虑气井井筒压力、温度和液流量条件影响的垂直气井井筒多液滴夹带临界气流量的计算方法，属于气井采气工艺技术领域，气井雾状流场中的液滴不计其数，而现有连续携液液滴夹带临界气流量计算方法未考虑井筒压力、温度及液流速等因素对雾状多液滴流场中液滴的动力学特征影响，也未考虑雾状流场中大量液滴相互作用对液滴受力情况的影响，导致适用条件受限。本发明基于垂直气井井筒雾状流场中多液滴动力学特征，提出了一种垂直气井井筒多液滴夹带临界气流量计算方法，所述方法的中的液滴的平均直径、曳力系数等关键参数与气井井筒流动条件联系紧密，具有较强的理论基础。该方法适用范围更宽泛，准确性更好。
Description
技术领域
本发明属于气井采气工艺技术领域，具体涉及一种垂直气井井筒雾状流场中多液滴夹带的临界携液气流量的计算方法。
背景技术
气井开始积液的产气量为气井连续携液临界气流量。当气井产量小于井筒压力及温度条件下的连续携液临界气流量时，气井开始积液。井筒积液将增加井底回压，降低气井产能，井筒积液量太大可使低压气井完全停喷。连续携液气井井筒流型是环雾流，气芯中的大量液滴被气流连续稳定夹带。气井连续携液生产的基础是气流量大于多液滴夹带的临界气流量。因此，准确预测气井井筒多液滴夹带的临界气流量对于优化气井配产气量具有十分重要的意义。
现有连续携液液滴夹带临界气流量计算方法，如Turner圆球模型、Coleman圆球模型、李闽椭球模型、王毅忠球帽模型，是通过对气流中的最大液滴进行受力分析得到，是建立在单液滴受力特征基础之上的，没有考虑井筒压力、温度及液流速等因素对雾状多液滴流场中液滴的动力学特征影响，没有考虑雾状流场中大量液滴相互作用对液滴受力情况的影响，导致现有单液滴模型适用条件受限。为此，本发明基于气井井筒雾状流场中多液滴动力学特征，提出了一种垂直气井井筒多液滴夹带临界气流量的计算方法。
发明内容
气井雾状流场中的液滴不计其数。现有连续携液液滴模型是通过对气流中的最大液滴进行受力分析得到，是建立在单液滴受力特征基础之上的，没有考虑井筒流动条件，如液流速、界面张力、气体密度、粘度等因素对雾状多液滴流场中液滴的动力学特征影响，没有考虑雾状流场中大量液滴相互作用对液滴受力情况的影响，导致适用条件受限。为此，本发明基于气井井筒雾状流场中多液滴动力学特征，提出了一种垂直气井井筒多液滴夹带临界携液气流量的计算方法。
本申请提供的一种垂直气井井筒多液滴夹带临界携液气流量的计算方法是这样实现的：
步骤一：假设多液滴夹带临界气流量qSC,Crit,new；
步骤二：根据压力p、温度T、管径D、液流量qL计算流体的物性参数和流动参数；
步骤三：计算雾状多液滴流场中的持液率HLC；
步骤四：计算雾状多液滴流场中液滴平均索特直径d32；
步骤五：计算雾状多液滴流场中液滴平均索特直径对应的韦伯数Wed32；
步骤六：计算雾状多液滴流场中液滴间横向及纵向距离do；
步骤七：计算雾状多液滴流场中液滴的曳力系数CD；
步骤八：计算雾状多液滴流场中液滴夹带的特征参数Cm,d32；
步骤九：计算雾状多液滴流场中液滴夹带的临界气流速uCrit；
步骤十：计算雾状多液滴流场中液滴夹带的临界携液气流量qSC,Crit；
步骤十一：若步骤一假设的气流量qSC,Crit,new与步骤十一计算的临界携液气流量qSC,Crit误差满足一定的精度，所假设的临界携液气流量qSC,Crit,new为该流动条件下的临界携液气流量；否则重复骤一至步骤十。
参照后文的说明和附图，详细公开了本申请的特定实施方式，指明了本申请相关参数的计算方式。应该理解本申请的实施方式在范围上并不因而受到限制。在所附权利要求的精神和条款的范围内，本申请的实施方式包括许多改变、修改和等同。
本发明的有益效果在于：
本发明考虑了井筒流动条件，如液流速、界面张力、气体密度、粘度等因素对雾状多液滴流场中液滴的动力学特征影响，同时考虑了雾状流场中大量液滴相互作用对液滴受力情况的影响，是基于气井井筒雾状流场中多液滴动力学特征提出的一种气井临界携液气流量计算方法。所述方法中的关键参数具有较强的理论基础，与现有的单液滴模型相比，是适用范围更宽泛，准确性更好。
附图说明
所包括的附图用来提供对本申请实施例的进一步的理解，其构成了说明书的一部分，用于例示本申请的实施方式，并与文字描述一起来阐释本申请的原理。显而易见地，下面描述中的附图仅仅是本申请的一些实施例。
图1为假象的多液滴单元，液滴排列成正六面体，分别在六面体的六个端部，液滴为圆球状，液滴尺寸均等。图中do为液滴间横向及纵向距离，d32为液滴平均索特直径。图1中每个液滴被8个单元共享，每个单元实际有效液滴个数为1。
图2是多液滴夹带临界气流量计算流程，图中ε为设定的偏差。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本申请一部分实施例，而不是全部的实施例。
本实施例提供一种垂直气井井筒多液滴夹带临界携液气流量的计算方法。包括如下步骤：
步骤一：假设多液滴夹带临界携液气流量qSC,Crit,new；
所假设的多液滴夹带临界携液气流量qSC,Crit,new是根据前一步计算的多液滴夹带临界携液气流量qSC,Crit进行调整。本专利采用牛顿迭代法求取该循环步多液滴夹带临界携液气流量的假设值。
（1）
式中 qSC,Crit,new为此循环新假设的多液滴夹带临界携液气流量；qSC,Crit,old为上一次循环假设的多液滴夹带临界携液气流量；f(qSC,Crit,old)为上一次循环后计算的多液滴夹带临界携液气流量，即步骤十一的qSC,crit；f'(qSC,Crit,old)为上一次循环后计算的多液滴夹带临界携液气流量随所假设值的梯度，即△qSC,Crit/△qSC,Crit,old。
步骤二：根据压力p、温度T、管径D、液流量qL计算流体的物性参数(气体偏差系数z、气体的密度ρG、液体密度ρL、气体的粘度μG、气液界面张力σ)和流动参数(表观气流速uSG、表观液流速uSL)；
步骤三：由下式计算雾状多液滴流场中的持液率HLC，即液滴所占体积份额。
（2）
式中HLC为持液率；uSG为表观气流速，m/s；uSL为表观液流速，m/s。
步骤四：由Azzopardi提出的半经验半理论关系式计算雾状多液滴流场中的液滴索特直径d32：
（3）
式中d32为液滴平均索特直径，m；λ为中间参数，其定义式为λ=(σ/ρLg)1/2；WeL为液相的韦伯数，其定义式为WeL=ρLuSG 2λ/σ；GLE为液相质量流量通量，计算式为GLE=uSLρL=qLρL/21600πD2；ρL为液体的密度，kg/m3；uSG为表观气流速，m/s；uSL为液流速，m/s；σ为气液界面张力，N/m；D为管径，m；g为重力加速度，9.8。
步骤五：由下式计算雾状多液滴流场中液滴平均索特直径对应的临界韦伯数Wed32。
（4）
式中Wed32为液滴平均索特直径对应的临界韦伯数；uSG为表观气流速，m/s；λ为中间参数，其定义式为(λ=σ/ρLg)；WeL为液相的韦伯数，其定义式为WeL=ρLuSG 2λ/σ；GLE为液相质量流量通量，计算式为GLE= uSLρL=qLρL/21600πD2；ρL为液体的密度，kg/m3；ρG为气体的密度，kg/m3；uSL为表观液流速，m/s；σ为气液界面张力，N/m；D为管径，m；g为重力加速度，9.8。
步骤六：由下式计算所述的雾状多液滴流场中液滴间横向及纵向距离do：
（5）
式中 do为液滴间距离，m；HLC为持液率；d32为液滴平均索特直径，m。
步骤七：由下式计算雾状多液滴流场中液滴的曳力系数CD：
（6）
式中CD为雾状多液滴相互作用下的曳力系数；0.42为单液滴的曳力系数；1-0.8(d0/d32)-0.45为考虑液滴相互影响对曳力系数提出的修正系数；do为液滴间距离，m；d32为液滴平均索特直径，m。
步骤八：由下式计算所述雾状多液滴流场中多液滴夹带的特征参数Cm,d32。
（7）
式中Cm,d32为雾状多液滴流场中多液滴夹带的特征参数；CD为雾状多液滴流场中液滴的曳力系数；Wed32为液滴平均索特直径对应的韦伯数；g为重力加速度，9.8。
步骤九：由下式计算雾状多液滴流场中多液滴夹带的临界气流速uCrit。
（8）
式中uCrit为雾状多液滴流场中多液滴夹带的临界气流速；Cm,d32为多液滴夹带的特征参数；ρL为液体的密度，kg/m3；ρG为气体的密度，kg/m3；σ为气液界面张力，N/m。
步骤十：由下式计算雾状多液滴流场中多液滴夹带的临界携液气流量qSC,Crit：
（9）
式中qSC,Crit为多液滴夹带的临界携液气流量，m3/d；uCrit为多液滴夹带的临界气流速，m/s；p为压力，MPa；A为油管过流截面积，m2；z为气体的偏差系数；T为温度，K。
实施例：某井深3000m，油管尺寸62mm，产液量为1m3/d，产液中不含油，产气量为3×104m3/d，测试井底流压为11.58MPa，测试井底温度为375.15K，天然气相对密度为0.59。以上述实例井说明所述多液滴夹带临界气流速计算方法的过程。
第一步：流体物性参数及流动参数计算：液体的密度ρL为960kg/m3，气体的密度ρG为68.28kg/m3，气体的偏差系数z取0.95，气体的粘度μG为0.015mPa.s，液体的粘度μG为0.28mPa.s，气液界面张力σ为0.047N/m，表观液流速uSL为0.0038m/s。
步骤二：假设多液滴夹带的临界气流量qSC,Crit,new为2×104m3/d，计算表观气流速uSG为0.852m/s；
步骤四：由式（3）计算雾状多液滴流场中的液滴索特直径d32为4.53mm；
步骤五：由式（4）计算雾状多液滴流场中液滴平均索特直径对应的临界韦伯数Wed32为4.82；
步骤六：由式（5）计算雾状多液滴流场中液滴间横向及纵向距离do为21.8mm；
步骤七：由式（6）计算雾状多液滴流场中液滴的曳力系数CD为0.25；
步骤八：由式（7）计算雾状多液滴流场中多液滴夹带的特征参数Cm,d32为3.96；
步骤九：由式（8）计算雾状多液滴流场中多液滴夹带的临界气流速uCrit为1.22 m/s；
步骤十：由式（9）计算雾状多液滴流场中多液滴夹带的临界携液气流量qSC,Crit为2.8×104m3/d；
步骤十一：比较判断，|qSC,Crit,new - qSC,Crit|/ qSC,Crit=0.4>0.005（计算精度）；说明所假设的雾状多液滴流场中多液滴夹带的临界气流量qSC,Crit,new不满足计算要求。
根据前一步假设的多液滴夹带的临界气流量qSC,Crit,new和计算的多液滴夹带的临界气流量qSC,Crit按照式（2）假设新的多液滴夹带的临界气流量qSC,Crit,new，重复步骤四到步骤十。
当假设的多液滴夹带的临界气流量qSC,Crit,new为3.02×104m3/d时，新计算的多液滴夹带的临界气流量qSC,Crit为3.05×104m3/d。所假设的多液滴夹带的临界气流量满足计算精度要求，计算结束。
同时以该井为例计算了不同压力及产液量条件下的多液滴夹带的临界气流量，如表1所示。从表1可知，随压力增加，多液滴夹带临界的气流量增加；随产液量增加，多液滴夹带的临界气流量增加。
表1气井不同产液量及井底压力下的多液滴夹带临界气流量。
Claims (11)
1.一种垂直气井井筒多液滴夹带临界气流量计算方法是建立在以下四点假设基础之上：①气井产出液体中全部以液滴形式夹带；②液滴在气芯中分布成正六面体排列；③井筒流场中液滴为圆球状；④液滴尺寸相等且等于平均液滴索特直径d32。
2.一种垂直气井井筒多液滴夹带临界气流量计算方法，计算步骤如下
步骤一：假设多液滴夹带临界气流量qSC,Crit,new；
步骤二：根据压力p、温度T、管径D、液流量qL计算流体的物性参数和流动参数；
步骤三：计算雾状多液滴流场中的持液率HLC；
步骤四：计算雾状多液滴流场中液滴平均索特直径d32；
步骤五：计算雾状多液滴流场中液滴平均索特直径对应的韦伯数Wed32；
步骤六：计算雾状多液滴流场中液滴间横向及纵向距离do；
步骤七：计算雾状多液滴流场中液滴的曳力系数CD；
步骤八：计算雾状多液滴流场中液滴夹带的特征参数Cm,d32；
步骤九：计算雾状多液滴流场中多液滴夹带的临界气流速uCrit；
步骤十：计算雾状多液滴流场中多液滴夹带的临界携液气流量qSC,Crit；
步骤十一：若步骤一假设的多液滴夹带临界气流量qSC,Crit,new与步骤十一计算的多液滴夹带临界气流量qSC,Crit误差满足一定的精度，所假设的qSC,Crit,new为该流动条件下的多液滴夹带临界气流量；否则重复骤一至步骤十。
3.如权利要求2所述的一种垂直气井井筒多液滴夹带临界气流量计算方法，其特征在于，所述的临界携液气流量qSC,Crit,new是根据前一步计算的临界携液气流量qSC,Crit进行调整，本专利采用牛顿迭代法求取该循环步的临界携液气流量的假设值：
<math display = 'block'>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>SC</mi>
<mo>,</mo>
<mi>Crit</mi>
<mo>,</mo>
<mi>new</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>SC</mi>
<mo>,</mo>
<mi>Crit</mi>
<mo>,</mo>
<mi>old</mi>
</mrow>
</msub>
<mo>&minus;</mo>
<mfrac>
<mrow>
<mi>f</mi>
<mo stretchy='false'>(</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>SC</mi>
<mo>,</mo>
<mi>Crit</mi>
<mo>,</mo>
<mi>old</mi>
</mrow>
</msub>
<mo stretchy='false'>)</mo>
</mrow>
<mrow>
<msup>
<mi>f</mi>
<mo>&#8242;</mo>
</msup>
<mo stretchy='false'>(</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>SC</mi>
<mo>,</mo>
<mi>Crit</mi>
<mo>,</mo>
<mi>old</mi>
</mrow>
</msub>
<mo stretchy='false'>)</mo>
</mrow>
</mfrac>
</mrow>
</math> （1）
式中 qSC,Crit,new为此循环新假设的携液气流量；qSC,Crit,old为上一次循环假设的携液气流量；f(qSC,Crit,old)为上一次循环后计算的携液临界气流量，即步骤十一的qSC,crit；f'(qSC,Crit,old)为上一次循环后计算的临界携液气流量随所假设临界携液气流量的梯度，即△qSC,Crit/△qSC,Crit,old。
4.如权利要求1所述的一种垂直气井井筒多液滴夹带临界气流量计算方法，其特征在于，所述的雾状多液滴流场中的持液率HLC由下式计算：
<math display = 'block'>
<mrow>
<msub>
<mi>H</mi>
<mi>LC</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>u</mi>
<mi>SL</mi>
</msub>
<mrow>
<msub>
<mi>u</mi>
<mi>SG</mi>
</msub>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>E</mi>
</msub>
<msub>
<mi>u</mi>
<mi>SL</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</math> （2）
式中 HLC为持液率；uSG为表观气流速，m/s；uSL为表观液流速，m/s。
5.如权利要求1所述的一种垂直气井井筒多液滴夹带临界气流量计算方法，其特征在于，所述的雾状多液滴流场中液滴平均索特直径d32由下式计算：
<math display = 'block'>
<mrow>
<msub>
<mi>d</mi>
<mn>32</mn>
</msub>
<mo>=</mo>
<mi>&lambda;</mi>
<mo stretchy='false'>(</mo>
<mfrac>
<mn>15.4</mn>
<msubsup>
<mi>We</mi>
<mi>L</mi>
<mn>0.58</mn>
</msubsup>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mn>3.5</mn>
<msub>
<mi>G</mi>
<mi>LE</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>u</mi>
<mi>SG</mi>
</msub>
</mrow>
</mfrac>
<mo stretchy='false'>)</mo>
</mrow>
</math> （3）
式中d32为液滴平均索特直径，m；λ为中间参数，其定义式为(λ=σ/ρLg)；WeL为液相的韦伯数，其定义式为WeL=ρLuSG 2λ/σ；GLE为液相质量流量通量，计算式为GLE=uSLρL=qLρL/21600πD2；ρL为液体的密度，kg/m3；uSG为表观气流速，m/s；uSL为液流速，m/s；σ为气液界面张力，N/m；D为管径，m；g为重力加速度，9.8。
6.如权利要求2所述的一种垂直气井井筒多液滴夹带临界气流量计算方法，其特征在于，雾状多液滴流场中液滴平均索特直径对应的临界韦伯数Wed32由下式计算：
<math display = 'block'>
<mrow>
<msub>
<mi>We</mi>
<mrow>
<mi>d</mi>
<mn>32</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>G</mi>
</msub>
<msubsup>
<mi>u</mi>
<mi>SG</mi>
<mn>2</mn>
</msubsup>
<mi>&lambda;</mi>
</mrow>
<mi>&sigma;</mi>
</mfrac>
<mo stretchy='false'>(</mo>
<mfrac>
<mn>15.4</mn>
<msubsup>
<mi>We</mi>
<mi>L</mi>
<mn>0.58</mn>
</msubsup>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mn>3.5</mn>
<msub>
<mi>G</mi>
<mi>LE</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<msub>
<mi>u</mi>
<mi>SG</mi>
</msub>
</mrow>
</mfrac>
<mo stretchy='false'>)</mo>
</mrow>
</math> （4）
式中Wed32为液滴平均索特直径对应的临界韦伯数；uSG为表观气流速，m/s；λ为中间参数，其定义式为(λ=σ/ρLg)；WeL为液相的韦伯数，其定义式为WeL=ρLuSG 2λ/σ；GLE为液相质量流量通量，计算式为GLE= uSLρL=qLρL/21600πD2；ρL为液体的密度，kg/m3；ρG为气体的密度，kg/m3；uSL为表观液流速，m/s；σ为气液界面张力，N/m；D为管径，m；g为重力加速度，9.8。
7.如权利要求1所述的一种垂直气井井筒多液滴夹带临界气流量计算方法，其特征在于，所述的雾状多液滴流场中液滴间横向及纵向距离do由下式计算：
<math display = 'block'>
<mrow>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msub>
<mi>d</mi>
<mn>32</mn>
</msub>
<msup>
<mrow>
<mo stretchy='false'>(</mo>
<mfrac>
<mi>&pi;</mi>
<mrow>
<mn>6</mn>
<msub>
<mi>H</mi>
<mi>LC</mi>
</msub>
</mrow>
</mfrac>
<mo stretchy='false'>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>3</mn>
</mrow>
</msup>
</mrow>
</math> （5）
式中 do为液滴间距离，m；HLC为持液率；d32为液滴平均索特直径，m。
8.如权利要求1所述的一种垂直气井井筒多液滴夹带临界气流量计算方法，其特征在于，所述的雾状多液滴流场中液滴的曳力系数CD由下式计算：
<math display = 'block'>
<mrow>
<msub>
<mi>C</mi>
<mi>D</mi>
</msub>
<mo>=</mo>
<mn>0.42</mn>
<mo stretchy='false'>(</mo>
<mn>1</mn>
<mo>&minus;</mo>
<mn>0.8</mn>
<mo stretchy='false'>(</mo>
<msup>
<mrow>
<mrow>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
<mo>/</mo>
<msub>
<mi>d</mi>
<mn>32</mn>
</msub>
</mrow>
<mo stretchy='false'>)</mo>
</mrow>
<mrow>
<mo>&minus;</mo>
<mn>0.45</mn>
</mrow>
</msup>
<mo stretchy='false'>)</mo>
</mrow>
</math> （6）
式中CD为雾状多液滴相互作用下的曳力系数；0.42为单液滴的曳力系数；1-0.8(d0/d32)-0.45为考虑液滴相互影响对曳力系数提出的修正系数；do为液滴间距离，m；d32为液滴平均索特直径，m。
9.如权利要求1所述的一种垂直气井井筒多液滴夹带临界气流量计算方法，其特征在于，所述雾状多液滴流场中多液滴夹带的特征参数Cm,d32由下式计算：
<math display = 'block'>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mi>d32</mi>
</mrow>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo stretchy='false'>(</mo>
<mfrac>
<mrow>
<mn>4</mn>
<mi>g</mi>
<msub>
<mi>We</mi>
<mi>d32</mi>
</msub>
</mrow>
<mrow>
<mn>3</mn>
<msub>
<mi>C</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<mo stretchy='false'>)</mo>
</mrow>
<mn>0.25</mn>
</msup>
</mrow>
</math> （7）
式中Cm,d32为雾状多液滴流场中多液滴夹带的特征参数；CD为雾状多液滴流场中液滴的曳力系数；Wed32为液滴平均索特直径对应的韦伯数；g为重力加速度，9.8。
10.如权利要求1所述的一种垂直气井井筒多液滴夹带临界气流量计算方法，其特征在于，所述雾状多液滴流场中多液滴夹带的临界气流速uCrit由下式计算：
<math display = 'block'>
<mrow>
<msub>
<mi>u</mi>
<mi>Crit</mi>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mi>d32</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo stretchy='false'>(</mo>
<mfrac>
<mrow>
<mi>&sigma;</mi>
<mo stretchy='false'>(</mo>
<msub>
<mi>&rho;</mi>
<mi>L</mi>
</msub>
<mo>&minus;</mo>
<msub>
<mi>&rho;</mi>
<mi>G</mi>
</msub>
<mo stretchy='false'>)</mo>
</mrow>
<msubsup>
<mi>&rho;</mi>
<mi>G</mi>
<mn>2</mn>
</msubsup>
</mfrac>
<mo stretchy='false'>)</mo>
</mrow>
<mn>0.25</mn>
</msup>
</mrow>
</math> （8）
式中uCrit为雾状多液滴流场中多液滴夹带的临界气流速；Cm,d32为流液滴夹带的特征参数；ρL为液体的密度，kg/m3；ρG为气体的密度，kg/m3；σ为气液界面张力，N/m。
11.如权利要求1所述的一种垂直气井井筒多液滴夹带临界气流量计算方法，其特征在于，所述雾状多液滴流场中多液滴夹带的临界携液气流量qSC,Crit由下式计算：
<math display = 'block'>
<mrow>
<msub>
<mi>q</mi>
<mrow>
<mi>SC</mi>
<mo>,</mo>
<mi>Crit</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>2.5</mn>
<mo>&times;</mo>
<msup>
<mn>10</mn>
<mn>8</mn>
</msup>
<mfrac>
<mrow>
<mi>pA</mi>
<msub>
<mi>u</mi>
<mi>Crit</mi>
</msub>
</mrow>
<mi>zT</mi>
</mfrac>
</mrow>
</math> （9）
式中qSC,Crit为雾状多液滴流场中多液滴夹带临界携液气流量，m3/d；uCrit为雾状多液滴流场中多液滴夹带的临界气流速，m/s；p为压力，MPa；A为油管流动截面积，m2；z为气体的偏差系数；T为温度，K。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711387764.0A CN108121894B (zh) | 2017-12-20 | 2017-12-20 | 一种垂直气井井筒多液滴夹带临界气流量计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711387764.0A CN108121894B (zh) | 2017-12-20 | 2017-12-20 | 一种垂直气井井筒多液滴夹带临界气流量计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108121894A true CN108121894A (zh) | 2018-06-05 |
| CN108121894B CN108121894B (zh) | 2019-09-10 |
Family
ID=62230763
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711387764.0A Active CN108121894B (zh) | 2017-12-20 | 2017-12-20 | 一种垂直气井井筒多液滴夹带临界气流量计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108121894B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109765265A (zh) * | 2019-01-24 | 2019-05-17 | 中国石油大学(华东) | 测量深水气井环空测试液保温性能的装置及方法 |
| CN111400978A (zh) * | 2020-06-08 | 2020-07-10 | 西南石油大学 | 一种考虑液滴形变和多参数影响的临界携液流量计算方法 |
| CN112345420A (zh) * | 2020-11-09 | 2021-02-09 | 中北大学 | 一种旋转填料床中液滴粒径及粒径分布的测试装置和方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104533388A (zh) * | 2014-12-23 | 2015-04-22 | 中国石油天然气股份有限公司 | 一种测量气井临界携液流量的方法 |
-
2017
- 2017-12-20 CN CN201711387764.0A patent/CN108121894B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104533388A (zh) * | 2014-12-23 | 2015-04-22 | 中国石油天然气股份有限公司 | 一种测量气井临界携液流量的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 周德胜: "气井携液多液滴模型研究", 《水动力学研究与进展》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109765265A (zh) * | 2019-01-24 | 2019-05-17 | 中国石油大学(华东) | 测量深水气井环空测试液保温性能的装置及方法 |
| CN111400978A (zh) * | 2020-06-08 | 2020-07-10 | 西南石油大学 | 一种考虑液滴形变和多参数影响的临界携液流量计算方法 |
| CN111400978B (zh) * | 2020-06-08 | 2020-09-29 | 西南石油大学 | 一种考虑液滴形变和多参数影响的临界携液流量计算方法 |
| CN112345420A (zh) * | 2020-11-09 | 2021-02-09 | 中北大学 | 一种旋转填料床中液滴粒径及粒径分布的测试装置和方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108121894B (zh) | 2019-09-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108121894A (zh) | 一种垂直气井井筒多液滴夹带临界气流量计算方法 | |
| CN110175340B (zh) | 一种适用于不确定工况的离心压气机叶轮气动稳健优化设计方法 | |
| CN113392618B (zh) | 一种印刷电路板掩膜线宽补偿量预测方法 | |
| CN108829960B (zh) | 激光增材再制造同轴送粉喷嘴粉流流场的数值模拟方法 | |
| CN108256262B (zh) | 一种轴对称射流稳压腔参数设计的数值模拟方法 | |
| CN110991760A (zh) | 一种高气液比产水气井临界携液流速预测方法 | |
| CN116401897B (zh) | 一种基于近似解和能量方程的半解析裂缝扩展模拟方法 | |
| Zhang et al. | Computational fluid dynamics (CFD) modeling of proppant transport in a plug and perf completion with different perforation phasing | |
| Jiang et al. | Discrete element method–computational fluid dynamics analyses of flexible fibre fluidization | |
| CN106595787B (zh) | 一种低产气井利用持气率计算气相流量的方法 | |
| Yuan et al. | Influence of droplet coalescence and breakup on the separation process in wave‐plate separators | |
| CN113987714A (zh) | 一种柱塞泵球面配流副油膜动态多场求解方法及系统 | |
| CN111025898B (zh) | 一种用于流程工业大规模过程控制的降维辨识方法 | |
| CN108090321A (zh) | 一种水平气井倾斜管液膜携带临界气流量计算方法 | |
| Demir et al. | Determination of the hydraulic properties of a flat type drip emitter using computational fluid dynamics | |
| CN112632767B (zh) | 一种均衡堵塞与冲蚀的适度挡砂精度设计方法 | |
| CN111967143B (zh) | 一种适于深水气井的临界携液流量预测方法 | |
| CN108629105A (zh) | 一种计算等离子喷涂液滴铺展过程流场的方法 | |
| CN105404749B (zh) | 一种空排止回阀的建模方法 | |
| Chen et al. | Interactions of spout jets in a multiple‐spouted bed | |
| CN110032818B (zh) | 一种无边界游泳池水跌流量计算方法 | |
| CN106640004A (zh) | 注汽锅炉出口的蒸汽热力参数的计算方法及其装置 | |
| WB et al. | Theoretical analysis of bubble formation in a co-flowing liquid | |
| Matsuzaki et al. | Numerical study on particle motions in swirling flows in a cyclone separator | |
| Hui et al. | Investigation on numerical computational methodology for the flow characteristics of aviation piston pump |
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 |