本申请涉及页岩气领域，尤其涉及一种页岩气水平井积液位置判断方法。
背景技术：
页岩气是一种储量大、开采难的非常规油气资源，近年来，国内长水平井分段压裂理论、技术和装备的快速发展、成熟和完善，页岩气实现了工业化、商业化的开发。而页岩气一般具有初期产量高、中后期产量递减快、返排率低的特点。压裂液长期同气体从地层到井底再到井口，需要消耗更多的地层能量，后期地层压力不足时，将会有液体停滞在井底或水平段，对页岩气的排采稳产工作造成技术难题，尤其是积液位置的不确定性，对排采工艺的选择无指导性。
国内外页岩气开发时间较短，后期排采理论与技术一般均是借鉴常规天然气手段，无法同页岩气长水平段的特殊性建立联系，目前尚未见到页岩气水平井积液位置判断方法的公开报道。
技术实现要素：
本申请提供一种页岩气水平井积液位置判断方法，解决现有技术中页岩气水平井中无法准确且系统的判断出积液的位置问题。
本申请的技术方案是：
一种页岩气水平井积液位置判断方法，包括以下步骤：
s1.计算页岩气水平井沿井筒分布的压力温度剖面，判断沿所述井筒方向的流型分布，计算所述井筒不同位置处的气体流速；
s2.计算天然气偏差系数随压力温度的变化，已知天然气组分组成，计算出天然气拟临界压力、拟临界温度；
s3.计算页岩气水平井的水平段、倾斜段以及垂直段三处的临界携液流量；
s4.计算临界携液流量并对比不同位置处的气体真实流量，根据井眼轨迹分布绘制不同参数下产量与临界携液流量沿井眼轨迹对比变化图，根据所述对比变化图判断出积液的位置；若是对比变化图中的临界携液气量大于实际产气量，则所述页岩气水平井段中有积液；若是对比变化图中的临界携液气量小于实际产气量，则所述页岩气水平井段中没有积液。
作为本申请的一种技术方案，在所述步骤s1中，采用beggs-brill方法计算页岩气水平井沿井筒分布的压力温度剖面，所述beggs-brill方法包括以下公式：
式中：p-压力；z-沿井筒方向的长度；
3.根据权利要求2所述的页岩气水平井积液位置判断方法，其特征在于，在所述步骤s1中，包括以下步骤：
s11.采用beggs-brill方法的公式，根据基础数据，确定计算的起始点压力p、计算段数n、段深度△h；
s12.假设计算段内的压力差为△p，则计算段内的末端压力ph1，计算所述计算段内的平均压力
s13.计算所述计算段内的压力差δp′与末端压力p′h1，若满足|p′h1-ph1|＜0.0001，则将计算的末端压力作为下一段的起点压力；若不满足|p′h1-ph1|＜0.0001，则重新将p′h1替代ph1，继续计算，直到满足误差要求为止；
s14.重复步骤s11至s13来继续计算下一段压力分布，直到所有计算分割段计算完成为止。
作为本申请的一种技术方案，步骤s11中的所述基础数据包括液体密度、气体密度、液相体积流量、气相体积流量、液体表面张力、重力加速度、液体黏度、气体黏度、井筒倾斜角、油管内径、井口温度以及温度梯度。
作为本申请的一种技术方案，在步骤s2中，采用dranchuk-abu-kassem方法计算天然气偏差系数随压力温度的变化，所述dranchuk-abu-kassem方法包括以下公式：
采用牛顿迭代法对z进行迭代计算，式中：z-，a1＝0.3265，a2＝-1.0700，tpr-拟临界温度，a3＝-0.5339，a4＝0.01569，a5＝-0.05165，ρpr-无因次对比密度，a6＝0.5457，a7＝-0.7361，a8＝0.1844，a9＝0.1056，a10＝0.6134，a11＝0.7210，ppr-拟临界压力。
6.根据权利要求1所述的页岩气水平井积液位置判断方法，其特征在于，在步骤s3中，所述页岩气水平井的临界产气量为：
式中：qcr-临界产气量，m3/d；a-油管横截面积，m2，
作为本申请的一种技术方案，在步骤s3中，所述页岩气水平井的垂直段临界速度为：
式中：g-重力加速度；σ-气水界面张力，n/cm；ρl-液体密度，kg/m3；ρg-气体密度，kg/m3。
作为本申请的一种技术方案，在步骤s3中，所述页岩气水平井的倾斜段临界速度为：
式中：g-重力加速度；σ-气水界面张力，n/cm；ρl-液体密度，kg/m3；ρg-气体密度，kg/m3；α-井斜角。
作为本申请的一种技术方案，在步骤s3中，所述页岩气水平井的水平井段临界速度为：
式中：g-重力加速度；σ-气水界面张力，n/cm；ρl-液体密度，kg/m3；ρg-气体密度，kg/m3。
本申请的有益效果：
本申请通过利用dranchuk-abu-kassem方法计算天然气偏差系数随压力温度的变化，采用beggs-brill方法计算沿井筒分布的压力温度剖面，判断沿井筒方向的流型分布，计算不同位置处的气体流速，通过计算临界携液流量并对比不同位置出的气体真实流量，根据井眼轨迹分布绘制不同参数下产量与临界携液流量沿井眼轨迹对比变化图，根据该对比变化图来系统有效的判断出页岩气水平井段中的积液的位置；该方法科学系统，能够有效地判断出页岩气水平井段的积液的位置，方法科学，且使用效率高。
具体实施方式
实施例：
本申请实施例提供了一种页岩气水平井积液位置判断方法，其主要包括以下步骤：
s1.采用beggs-brill方法计算页岩气水平井沿井筒分布的压力温度剖面，判断沿井筒方向的流型分布，计算井筒不同位置处的气体流速；
s2.采用dranchuk-abu-kassem方法计算天然气偏差系数随压力温度的变化，已知天然气组分组成，计算出天然气拟临界压力、拟临界温度；
s3.计算页岩气水平井的水平段、倾斜段以及垂直段三处的临界携液流量；
s4.计算临界携液流量并对比不同位置处的气体真实流量，根据井眼轨迹分布绘制不同参数下产量与临界携液流量沿井眼轨迹对比变化图，根据对比变化图判断出积液的位置；若是对比变化图中的临界携液气量大于实际产气量，则页岩气水平井段中有积液；若是对比变化图中的临界携液气量小于实际产气量，则页岩气水平井段中没有积液。
需要说明的是，在本实施例中，beggs-brill方法包括以下公式：
式中：p-压力；z-沿井筒方向的长度；
需要说明的是，在本实施例中，dranchuk-abu-kassem方法包括以下公式：
采用牛顿迭代法对z进行迭代计算，式中：z-，a1＝0.3265，a2＝-1.0700，tpr-拟临界温度，a3＝-0.5339，a4＝0.01569，a5＝-0.05165，ρpr-无因次对比密度，a6＝0.5457，a7＝-0.7361，a8＝0.1844，a9＝0.1056，a10＝0.6134，a11＝0.7210，ppr-拟临界压力。需要说明的是，该方法适用于1.0≤tpr≤3.0，或0.7≤tpr≤1.0，ppr＜1.0的情况。
进一步地，在步骤s1中，其主要包括以下步骤：
s11.采用beggs-brill方法的公式，根据基础数据，确定计算的起始点压力p、计算段数n、段深度△h；
s12.假设计算段内的压力差为△p，则计算段内的末端压力ph1，计算计算段内的平均压力
s13.计算计算段内的压力差δp′与末端压力p′h1，若满足|p′h1-ph1|＜0.0001，则将计算的末端压力作为下一段的起点压力，若不满足|p′h1-ph1|＜0.0001，则重新将p′h1替代ph1，继续计算直到满足误差要求为止；
s14.重复步骤s11至s13来继续计算下一段压力分布，直到所有计算分割段计算完成。
需要说明的是，在本实施例中，步骤s11中的基础数据包括液体密度、气体密度、液相体积流量、气相体积流量、液体表面张力、重力加速度、液体黏度、气体黏度、井筒倾斜角、油管内径、井口温度以及温度梯度。
需要说明的是，在本实施例中，在步骤s3中，页岩气水平井的临界产气量为：
式中：qcr-临界产气量，m3/d；a-油管横截面积，m2，
需要说明的是，在本实施例中，在步骤s3中，页岩气水平井的垂直段临界速度为：
式中：g-重力加速度；σ-气水界面张力，n/cm；ρl-液体密度，kg/m3；ρg-气体密度，kg/m3。
需要说明的是，在本实施例中，在步骤s3中，页岩气水平井的倾斜段临界速度为：
式中：g-重力加速度；σ-气水界面张力，n/cm；ρl-液体密度，kg/m3；ρg-气体密度，kg/m3；α-井斜角。
需要说明的是，在本实施例中，在步骤s3中，页岩气水平井的水平井段临界速度为：
式中：g-重力加速度；σ-气水界面张力，n/cm；ρl-液体密度，kg/m3；ρg-气体密度，kg/m3。
需要说明的是，在本实施例中，本申请通过利用dranchuk-abu-kassem方法计算天然气偏差系数随压力温度的变化，采用beggs-brill方法计算沿井筒分布的压力温度剖面，判断沿井筒方向的流型分布，计算不同位置处的气体流速，通过计算临界携液流量并对比不同位置出的气体真实流量，根据井眼轨迹分布绘制不同参数下产量与临界携液流量沿井眼轨迹对比变化图，根据该对比变化图来系统有效的判断出页岩气水平井段中的积液的位置；该方法科学系统，能够有效地判断出页岩气水平井段的积液的位置，方法科学，且使用效率高。
以上所述仅为本申请的优选实施例而已，并不用于限制本申请，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，应包含在本申请的保护范围之内。