倾斜有水气藏产水气井产能的计算方法技术

本发明专利技术公开了一种倾斜有水气藏产水气井产能的计算方法，本发明专利技术从渗流理论入手，充分考虑有水气藏中储层的倾角对气井产能的影响，进而建立产水气井的产能计算方程，为倾斜有水气藏气井的产能预测和合理工作制度的制定提供合理的理论依据。本发明专利技术的优越性表现在：充分考虑了气藏的储层倾角；考虑了气藏具有边底水的情况，建立了气井产水对气井产能的影响关系；考虑了气水两相高速非达西渗流；建立了应用相对渗透率曲线，并考虑去除凝析水的影响下，绘制出产水率与含水饱和度的关系曲线。

【技术实现步骤摘要】

本专利技术涉及一种产水气井产能的计算方法，具体涉及一种倾斜有水气藏产水气井产能的计算方法。
技术介绍
随着气藏采出程度的增加，地层中存在的可动水会不断在井底聚集，造成气井生产水气比上升，严重影响气井产量和产能，并影响整个气田的储量评估及开发方案的实施。因此产水气井产能的确定和液相伤害程度对制定合理气井开采措施而言显得尤为重要。在气井产水对产能的影响方面前人做了大量研究，主要表现在：(1)考虑真实气体PVT参数随压力变化的高速非达西渗流单井模型的气井产能[1-3]；(2)采用井周渗透性变化对气井产能的影响推导气井产能方程[4]；(3)基于气井稳定和拟稳定流动状态产能方程推导产水气井的产能方程[5-13]；(4)考虑应力敏感的产水气井的产能方程，并通过实验和理论综合分析对气井产能的影响[14-15]。然而绝大部分气藏均为背斜气藏，且带有一定的倾角，前人研究建立的相关产水气井的产能方程都没有考虑倾斜有水气藏储层角度对产能的影响。
技术实现思路
针对现有技术中的缺陷，本专利技术提供了一种倾斜有水气藏产水气井产能的计算方法，该计算方法从渗流理论入手充分考虑到有水气藏中储层的倾角对气井产能的影响。本专利技术采取的技术方案如下：一种倾斜有水气藏产水气井产能的计算方法，所述倾斜有水气藏具备以下特征：具有一定储层倾角θ；气水彼此不互溶；储层全部参与生产，地层中的气和水径向流入井内；地层中的气和水微可压缩，且压缩系数为常数；流体黏度为常数，考虑气水两相高速非达西渗流而不考虑启动压力梯度；忽略毛管力的影响；流体为等温流动；产能计算方程为：式中，pe：地层压力，单位MPa；pwf：井底流动压力，单位MPa；ψ(pe)：压力为pe时的气水两相拟压力，单位MPa；ψ(pwf)：压力为pwf时的气水两相拟压力，单位MPa；A：产能方程达西系数；qsc：标准状况(温度为0℃、压力为1标准大气压)下气井的气体体积流量，单位m3/s；B：产能方程非达西系数；re：气藏控制半径，单位m；rw：井筒半径，单位m；Krw、Krg分别为水相和气相的相对渗透率，无因次；ρw、ρg分别为水和气体的密度，单位kg/m3；g为重力加速度，m/s2；θ为储层倾角，0°≤θ≤90°；μw、μg分别为水相和气相的黏度，单位mPa·s。优选的，式中，μw、μg分别为水相和气相的黏度，单位mPa·s；Krw、Krg分别为水相和气相的相对渗透率，无因次；a为水气质量比，单位kg/kg；ρsc为标准状况(温度为0℃、压力为1标准大气压)下气体的密度，单位kg/m3；h为油层厚度，单位m；re为气藏控制半径，单位m；rw为井筒半径，单位m；δ为常数7.644×1010；K为气藏渗透率，单位10-3μm2；r为气体渗流半径，单位m，rw≤r≤re；表皮系数为S，无因次。优选的，a为mw/mg，mw和mg分别为气、水的质量流量，单位kg/s；气体质量流量mg＝qscρsc，ρsc为标准状况(温度为0℃、压力为1标准大气压)下气体的密度，单位kg/m3；qsc为标准状况(温度为0℃、压力为1标准大气压)下气体的体积流量，单位m3/s。优选的，通过气水两相拟压力函数的定义：得到p＝pwf，p＝pe，pe为地层压力，单位MPa；pwf为井底流动压力，单位MPa。优选的，考虑气水两相高速非达西渗流和含有地层倾角的运动方程为：式中，θ为储层倾角，单位°；K为气藏渗透率，单位10-3μm2；Krw、Krg分别为水相和气相的相对渗透率，无因次；pw、pg分别为水相和气相的压力，单位MPa；Vw、Vg分别为水相和气相的速度，单位m/s；μw、μg分别为水相和气相的黏度，单位mPa·s；βw、βg分别为水相和气相的速度系数，单位m-1；ρw、ρg分别为水和气体的密度，单位kg/m3；水相和气相的速度系数为βw＝δ/Kw1.5，βg＝δ/Kg1.5，δ为常数7.644×1010，Kg、Kw分别为水相和气相的渗透率，单位10-3μm2；g为重力加速度，单位m/s2；r为气体渗流半径，单位m。优选的，由于忽略毛管力的影响，则pw＝pg＝p；式中，pw、pg分别为水相和气相压力，p为气藏压力，单位MPa。优选的，Vw和Vg由以下计算方法得到：式中，mg、mw分别为气、水的质量流量，单位kg/s；h为油层厚度，单位m；ρw、ρg分别为水和气体的密度，单位kg/m3；r为气体渗流半径，单位m，rw≤r≤re。优选的，气体密度ρg根据ρg＝m/v＝PMg/RT计算获得；其中，P表示绝对压力，单位MPa；R为摩尔气体常数0.008471；T表示绝对温度，单位K；m表示气体质量，单位kg；Mg表示气体相对平均分子量，计算公式为yi表示气体组分i的摩尔分数；Mi表示气体组分i的相对分子量；n表示气体的组分数量。优选的，根据天然气的组分数据，计算μg与p的关系曲线，并获得pe、pwf值下的μg，式中，μg为气相的黏度，单位mPa·s；Mg表示气体相对平均分子量；T表示绝对温度，单位K；ρg为气体的密度，单位kg/m3。优选的，Krw、Krg由以下方法获得：(1)根据含水率公式利用相对渗透率曲线绘制含水率与含水饱和度的关系曲线；式中，WGR表示生产水气比，单位m3/104m3；Rwgr表示凝析水水气比，单位m3/104m3；fw为含水率，％；(2)根据含水率的定义采用实际气藏中气井的水气比a，计算对应的含水率fw，然后根据(1)中的曲线，查找出对应含水率下的含水饱和度Sw，进而在相对渗透率曲线上查得Sw对应的Krw、Krg。本专利技术的有益效果在于：本专利技术从渗流理论入手，充分考虑有水气藏中储层的倾角对气井产能的影响，进而建立产水气井的产能模型，为倾斜有水气藏气井的产能预测和合理工作制度的制定提供合理的理论依据。本专利技术的优越性表现在：(1)本专利技术所述计算方法充分考虑了气藏的储层倾角，角度可以从0-90°，更加符合气藏的真实特性，当倾角为0°时，与常规气藏产能一致；(2)本专利技术所述计算方法考虑了气藏具有边底水的情况，建立了气井产水对气井产能的影响关系；(3)本专利技术所述计算方法考虑了气水两相高速非达西渗流；(4)本专利技术所述计算方法建立了应用相对渗透率曲线，并考虑去除凝析水的影响下，绘制出产水率与含水饱和度的关系曲线。附图说明图1倾斜气藏运动方程示意图，θ为储层倾角，g为重力加速度。图2含水率fw与含水饱和度Sw关系示意图。图3水气两相相对渗透率Krw、Krg曲线示意图。具体实施方式现结合附图对本专利技术的优选实施例进行详细描述。如图1所示，所述倾斜有水气藏具备以下特征：具有一定储层倾角θ；气水彼此不互溶；储层全部参与生产，地层中的气和水径向流入井内；地层中的气和水微可压缩，且压缩系数为常数；流体黏度为常数，考虑气水两相高速非达西渗流而不考虑启动压力梯度；忽略毛管力的影响；流体为等温流动。以下式中，θ为储层倾角，0°≤θ≤90°；K为气藏渗透率，单位10-3μm2；Krw、Krg分别为水相和气相的相对渗透率，无因次；pw、pg分别为水相和气相的压力，单位MPa；p为气藏压力，单位MPa；Vw、Vg分别为水相和气相的速度，单位m/s；μw、μg分别为水相和气相的黏度，单位mPa·s；βw本文档来自技高网...

【技术保护点】
一种倾斜有水气藏产水气井产能的计算方法，其特征在于，所述倾斜有水气藏具备以下特征：具有一定储层倾角θ；气水彼此不互溶；储层全部参与生产，地层中的气和水径向流入井内；地层中的气和水微可压缩，且压缩系数为常数；流体黏度为常数，考虑气水两相高速非达西渗流而不考虑启动压力梯度；忽略毛管力的影响；流体为等温流动；产能计算方程为：ψ(pe)-ψ(pwf)=Aqsc2+Bqsc+(re-rw)·(Krwρw2gsinθμw+Krgρg2gsinθμg),]]>式中，pe：地层压力，单位MPa；pwf：井底流动压力，单位MPa；ψ(pe)：压力为pe时的气水两相拟压力，单位MPa；ψ(pwf)：压力为pwf时的气水两相拟压力，单位MPa；A：产能方程达西系数；qsc：标准状况(温度为0℃、压力为1标准大气压)下气井的气体体积流量，单位m3/s；B：产能方程非达西系数；re：气藏控制半径，单位m；rw：井筒半径，单位m；Krw、Krg分别为水相和气相的相对渗透率，无因次；ρw、ρg分别为水和气体的密度，单位kg/m3；g为重力加速度，m/s2；θ为储层倾角，0°≤θ≤90°；μw、μg分别为水相和气相的黏度，单位mPa·s。...

【技术特征摘要】
1.一种倾斜有水气藏产水气井产能的计算方法，其特征在于，所述倾斜有水气藏具备以下特征：具有一定储层倾角θ；气水彼此不互溶；储层全部参与生产，地层中的气和水径向流入井内；地层中的气和水微可压缩，且压缩系数为常数；流体黏度为常数，考虑气水两相高速非达西渗流而不考虑启动压力梯度；忽略毛管力的影响；流体为等温流动；产能计算方程为：ψ(pe)-ψ(pwf)=Aqsc2+Bqsc+(re-rw)·(Krwρw2gsinθμw+Krgρg2gsinθμg),]]>式中，pe：地层压力，单位MPa；pwf：井底流动压力，单位MPa；ψ(pe)：压力为pe时的气水两相拟压力，单位MPa；ψ(pwf)：压力为pwf时的气水两相拟压力，单位MPa；A：产能方程达西系数；qsc：标准状况(温度为0℃、压力为1标准大气压)下气井的气体体积流量，单位m3/s；B：产能方程非达西系数；re：气藏控制半径，单位m；rw：井筒半径，单位m；Krw、Krg分别为水相和气相的相对渗透率，无因次；ρw、ρg分别为水和气体的密度，单位kg/m3；g为重力加速度，m/s2；θ为储层倾角，0°≤θ≤90°；μw、μg分别为水相和气相的黏度，单位mPa·s。2.根据权利要求1所述的倾斜有水气藏产水气井产能的计算方法，其特征在于，A=(μwKrw0.5+a2μgKrg0.5)ρsc2μgKrg0.5μwKrw0.54π2h2∫rwreδK1.5·1r2dr;]]>B=1K·2πhln(rerw+S)×(1+a)ρsc;]]>式中，μw、μg分别为水相和气相的黏度，单位mPa·s；Krw、Krg分别为水相和气相的相对渗透率，无因次；a为水气质量比，单位kg/kg；ρsc为标准状况(温度为0℃、压力为1标准大气压)下气体的密度，单位kg/m3；h为油层厚度，单位m；re为气藏控制半径，单位m；rw为井筒半径，单位m；δ为常数7.644×1010；K为气藏渗透率，单位10-3μm2；r为气体渗流半径，单位m，rw≤r≤re；表皮系数为S，无因次。3.根据权利要求2所述的倾斜有水气藏产水气井产能的计算方法，其特征在于，a为mw/mg，mw和mg分别为气、水的质量流量，单位kg/s；气体质量流量mg＝qscρsc，ρsc为标准状况(温度为0℃、压力为1标准大气压)下气体的密度，单位kg/m3；qsc为标准状况(温度为0℃、压力为1标准大气压)下气体的体积流量，单位m3/s。4.根据权利要求1所述的倾斜有水气藏产水气井产能的计算方法，其特征在于，通过气水两相拟压力函数的定义：得到p＝pwf，p＝pe，pe为地层压力，单位MPa；pwf为井底流动压力，单位MPa。5.根据权利要求1所述的倾斜有水气藏产水气井产能的计算方法，其特征在于，考虑气水两相高速非达西渗流和含有地层倾角的运动方程为：∂pw∂r=μwKKrwVw+βwρwVw2+ρwgsinθ---...

【专利技术属性】
技术研发人员：黄小亮，戚志林，严文德，雷登生，袁迎中，李继强，向祖平，肖前华，
申请(专利权)人：重庆科技学院，
类型：发明
国别省市：重庆;50