凝析气藏动态产能预测方法技术

公开了一种凝析气藏动态产能预测方法。该方法可以包括：步骤1：设定初始压力、起始时间tstart、结束时间与时间间隔△t，获取物性参数；步骤2：计算油相渗流方程与气相渗流方程，获得产油气速度表达式；步骤3：将地层划分为N个等间距的环形地层；步骤4：分别计算每一个环形地层在tstart的产油与产气速度；步骤5：分别计算每一个环形地层在时间t到时间t+△t内的油气含量变化量；步骤6：计算在时间t+△t每一个环形地层的地层压力，替换初始压力pstart，计算每一个环形地层在时间t+△t的产油速度与产气速度。本发明专利技术明确毛管力与相态变化对凝析气藏产能的影响，在此基础上准确预测凝析气藏在不同时刻的油气产量及对应井流物组成。

Dynamic Productivity Prediction Method for Condensate Gas Reservoir

【技术实现步骤摘要】
凝析气藏动态产能预测方法
本专利技术涉及油气田开发领域，更具体地，涉及一种凝析气藏动态产能预测方法。
技术介绍
凝析气在原始地层压力下一般呈单一气相状态。但在开发过程中，随着气体不断采出，地层压力下降，在上露点之下开始析出凝析油，形成油、气两相。伴随压力与相态的变化，地层油气流动过程十分复杂，常规气井产能预测与评价方法并不适用，如何建立一种新的适用于凝析气藏的动态产能预测方法，关键在于如何将油气相态随压力的变化融于油气渗流模型当中。1950年Fevang和Whitson在《ModelingGas-CondensateWellDeliverability》一文中利用数值模拟研究了凝析气藏的生产动态。通过分析生产气油比(GOR)随时间的变化，考虑流动过渡区的影响，提出三区流动模型:①内区，靠近井筒的区域，油气均可流动；②中间区，凝析油聚集的区域，仅有气相流动；③外区，地层压力高于露点，没有凝析油析出的单相气流区域。1985年Jones和Raghavan在《InterpretationofFlowingWellResponseinGasCondensateWells》中提出用油藏积分和井底积分两相拟压力函数方法研究凝析气藏不稳定渗流问题。2004年康晓东等在《凝析气藏考虑毛管数和非达西效应的渗流特征》中开展了考虑高速非达西与毛管数效应对油气两相流动的影响，在定解条件下，得到了拟稳态形式的流入动态方程。2006年石德佩等在《考虑相变的凝析气井产能方程》一文中运用稳态理论和拟压力方法,建立了考虑相态变化的凝析气井多相流产能方程，并比较了二项式和指数式产能方程的结果。2014年卢德堂等在《考虑闪蒸的凝析气井动态产能计算》中运用SRK状态方程进行闪蒸计算，明确压力与饱和度关系，并通过拟压力方法建立凝析气井产能方程。2014年张安刚等在《一种新的凝析气藏无因次产能评价方程》文中考虑气体非达西流动，将拟压力函数进行泰勒展开，通过稳态理论和多项式回归得到凝析气井无因次产能公式。目前国内外对凝析气藏产能预测方法的研究主要采用拟压力的方法，其结果较为准确可靠。但是还存在以下问题：①多相渗流微分方程的线性化过程中忽略了许多非线性项；②拟压力方法得到的产能需要转化后才能得到实际产能随时间的变化；③低渗透储层中毛管力对反凝析堵塞的影响较大，但是拟压力方法普遍忽略了毛管力对油气渗流以及反凝析伤害的影响。因此，有必要开发一种凝析气藏动态产能预测方法。公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
本专利技术提出了一种凝析气藏动态产能预测方法，其能够明确毛管力与相态变化对凝析气藏产能的影响，并且在此基础上准确预测凝析气藏在不同时刻的油气产量及对应井流物组成。所述方法可以包括：步骤1：设定初始压力pstart、起始时间tstart、结束时间tend与时间间隔△t，根据相平衡计算，获取物性参数；步骤2：根据所述物性参数，计算油相渗流方程与气相渗流方程，进而获得产油速度表达式与产气速度表达式；步骤3：以井筒中心为圆心，将地层划分为N个等间距的环形地层；步骤4：根据初始压力pstart、所述产油速度表达式与所述产气速度表达式，分别计算每一个环形地层在tstart的产油速度与产气速度；步骤5：针对每一个环形地层，根据所述环形地层在时间t的产油速度与产气速度，分别计算每一个环形地层在时间t到时间t+△t内的油含量变化量与气含量变化量；步骤6：针对每一个环形地层，根据所述环形地层在时间t到时间t+△t内的油气含量变化量，计算在时间t+△t每一个环形地层的地层压力，以所述地层压力替换初始压力pstart，计算每一个环形地层在时间t+△t的产油速度与产气速度；其中，时间t＝tstart+i△t，i＝0,1，…，优选地，所述油相渗流方程为：其中，k为储层渗透率，r为与井中心的径向距离，po为油相压力，μo为油相粘度，kro为油相的相对渗透率，vo为油相的渗流速度。优选地，所述产油速度表达式为：其中，qo表示产油速度，h表示油层厚度。优选地，所述环形地层On的产油速度为：其中，qon表示产油速度，kron表示环形地层On中的油相相对渗透率，pn表示环形地层On的地层压力，rn表示环形地层On的半径，μon为环形地层On中的油相粘度，n表示环形地层序号，n＝1,2,…,N。优选地，所述环形地层On在时间t到时间t+△t内的油含量变化量为：其中，表示环形地层On在时间t到时间t+△t内的油含量变化量，n＝1,2,…,N。优选地，所述气相渗流方程为：其中，pg为气相的压力，μg为气相粘度，krg为气相的相对渗透率，vg为气相的渗流速度，ρg为气相密度，βg为气相非达西系数，其中，根据公式(3)计算气相非达西系数：其中，φ表示孔隙度，Sg表示气相饱和度。优选地，所述产气速度表达式为：其中，qg表示产气速度。优选地，所述环形地层On的产气速度为：其中，qgn表示环形地层On的产气速度，krgn表示环形地层On中的气相相对渗透率，ρgn表示环形地层On中的气相密度，μgn为环形地层On中的气相粘度，βgn为环形地层On的气相非达西系数。优选地，所述环形地层On在时间t到时间t+△t内的气含量变化量为：其中，表示环形地层On在时间t到时间t+△t内的气含量变化量，Mgn为环形地层On的气相摩尔质量。优选地，t+Δt时间后环形地层On的地层压力为：其中，R为理想气体常数，T为温度，Zgn为环形地层On的气相偏差系数。本专利技术具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本专利技术的特定原理。附图说明通过结合附图对本专利技术示例性实施例进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施例中，相同的参考标号通常代表相同部件。图1示出了根据本专利技术的凝析气藏动态产能预测方法的步骤的流程图。图2示出了根据本专利技术的一个实施例的地层离散化的示意图。图3示出了根据本专利技术的一个实施例的油气相对渗透率曲线的示意图。图4示出了根据本专利技术的一个实施例的产油/产气速度对比的示意图。图5示出了根据本专利技术的一个实施例的累积产油/产气量对比的示意图。图6示出了根据本专利技术的一个实施例的生产气油比对比的示意图。图7示出了根据本专利技术的一个实施例的压力场径向分布的示意图。图8示出了根据本专利技术的一个实施例的井流物组成变化的示意图。具体实施方式下面将参照附图更详细地描述本专利技术。虽然附图中显示了本专利技术的优选实施例，然而应该理解，可以以各种形式实现本专利技术而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本专利技术更加透彻和完整，并且能够将本专利技术的范围完整地传达给本领域的技术人员。图1示出了根据本专利技术的凝析气藏动态产能预测方法的步骤的流程图。在该实施例中，根据本专利技术的凝析气藏动态产能预测方法可以包括：步骤1：设定初始压力pstart、起始时间tstart、结束时间tend与时间间隔△t，根据相平衡计算，获取物性参数；步骤2：根本文档来自技高网...

【技术保护点】
1.一种凝析气藏动态产能预测方法，包括：步骤1：设定初始压力pstart、起始时间tstart、结束时间tend与时间间隔△t，根据相平衡计算，获取物性参数；步骤2：根据所述物性参数，计算油相渗流方程与气相渗流方程，进而获得产油速度表达式与产气速度表达式；步骤3：以井筒中心为圆心，将地层划分为N个等间距的环形地层；步骤4：根据初始压力pstart、所述产油速度表达式与所述产气速度表达式，分别计算每一个环形地层在tstart的产油速度与产气速度；步骤5：针对每一个环形地层，根据所述环形地层在时间t的产油速度与产气速度，分别计算每一个环形地层在时间t到时间t+△t内的油含量变化量与气含量变化量；步骤6：针对每一个环形地层，根据所述环形地层在时间t到时间t+△t内的油气含量变化量，计算在时间t+△t每一个环形地层的地层压力，以所述地层压力替换初始压力pstart，计算每一个环形地层在时间t+△t的产油速度与产气速度；其中，时间t＝tstart+i△t，i＝0,1，…，

【技术特征摘要】
1.一种凝析气藏动态产能预测方法，包括：步骤1：设定初始压力pstart、起始时间tstart、结束时间tend与时间间隔△t，根据相平衡计算，获取物性参数；步骤2：根据所述物性参数，计算油相渗流方程与气相渗流方程，进而获得产油速度表达式与产气速度表达式；步骤3：以井筒中心为圆心，将地层划分为N个等间距的环形地层；步骤4：根据初始压力pstart、所述产油速度表达式与所述产气速度表达式，分别计算每一个环形地层在tstart的产油速度与产气速度；步骤5：针对每一个环形地层，根据所述环形地层在时间t的产油速度与产气速度，分别计算每一个环形地层在时间t到时间t+△t内的油含量变化量与气含量变化量；步骤6：针对每一个环形地层，根据所述环形地层在时间t到时间t+△t内的油气含量变化量，计算在时间t+△t每一个环形地层的地层压力，以所述地层压力替换初始压力pstart，计算每一个环形地层在时间t+△t的产油速度与产气速度；其中，时间t＝tstart+i△t，i＝0,1，…，2.根据权利要求1所述的凝析气藏动态产能预测方法，其中，所述油相渗流方程为：其中，k为储层渗透率，r为与井中心的径向距离，po为油相压力，μo为油相粘度，kro为油相的相对渗透率，vo为油相的渗流速度。3.根据权利要求2所述的凝析气藏动态产能预测方法，其中，所述产油速度表达式为：其中，qo表示产油速度，h表示油层厚度。4.根据权利要求3所述的凝析气藏动态产能预测方法，其中，环形地层On的产油速度为：其中，qon表示产油速度，kron表示环形地层On中的油相相对渗透率，pn表...

【专利技术属性】
技术研发人员：杨阳，吕成远，伦增珉，王锐，胡伟，王欣，周霞，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：北京,11