CN105134196A - 一种缝洞型碳酸盐岩凝析气井开发指标分析方法及装置 - Google Patents
一种缝洞型碳酸盐岩凝析气井开发指标分析方法及装置 Download PDFInfo
- Publication number
- CN105134196A CN105134196A CN201510556736.1A CN201510556736A CN105134196A CN 105134196 A CN105134196 A CN 105134196A CN 201510556736 A CN201510556736 A CN 201510556736A CN 105134196 A CN105134196 A CN 105134196A
- Authority
- CN
- China
- Prior art keywords
- gas
- condensate
- production
- pressure
- place
- 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
Abstract
本发明提供一种缝洞型碳酸盐岩凝析气井开发指标分析方法及装置,该方法包括:获取单井中的标准参数,根据标准参数计算单井动态储量;根据预设的天然气产量分段递减规律确定天然气产量以及各生产时间点的天然气累计产量;根据PVT实验数据确定两相偏差系数、气油比数据;建立物质平衡方程,将天然气累计产量及两相偏差系数代入方程,计算各生产时间点的地层压力理论值,并根据实测静压数据校正生成地层压力校正值;根据地层压力校正值及气油比数据预测对应的凝析油产量;以废弃地层压力为约束,根据物质平衡方程计算凝析油累计产量;再结合单井动态储量分别计算凝析油及天然气采出程度。
Description
技术领域
本发明属于凝析气藏开发技术领域,具体地,是一种关于缝洞型碳酸盐岩凝析气井开发指标分析的方法及装置。
背景技术
对于缝洞型碳酸盐岩凝析气藏开发来说,具有三方面的典型特征:1)储层性质复杂,发育不同程度的孔、洞、缝等储集介质,储层具有极强的非均质性;2)流体性质复杂,油、气状态随储层压力而发生相态变化;3)生产规律复杂,凝析油产量、天然气产量、地层压力递减快,不能按照定井底流压或井口定产的方式进行生产。总体上看,该类凝析气藏储层非均质性极强、井间连通性差、单井控制范围有限,普遍表现为“一井一藏”的模式,且储层内存在较为复杂的流体相态变化,每个储渗单元都具有不同的生产特征。
国内外针对缝洞型碳酸盐岩凝析气藏开发的研究相对较少,常规的分析方法对于该类气藏适用性较差,并具有很大的局限性,且预测结果准确性低。例如:若用数值模拟方法预测缝洞型凝析气井的开发指标,除地质建模困难、难以提供相对可靠的地质模型外,采用常规渗流数学模型描述大洞、大孔及大缝等大尺度介质内流体的流动规律,其适用性还存在较大分歧。综上所述,对于缝洞型凝析气藏关键开发指标(如凝析油产量、天然气产量、地层压力等)的准确预测目前仍是一个亟待解决的技术难题。
发明内容
本发明实施例的主要目的在于提供一种缝洞型碳酸盐岩凝析气井开发指标分析方法及装置,以克服现有技术中存在的上述问题,从而实现对缝洞型凝析气藏关键开发指标的准确分析。
为了实现上述目的,本发明实施例提供一种缝洞型碳酸盐岩凝析气井开发指标分析方法,所述的分析方法包括:获取单井中的标准参数,根据所述标准参数计算单井动态储量,所述单井动态储量包括:天然气地质储量和凝析油地质储量;根据预设的天然气产量分段递减规律预测天然气产量,并预测各生产时间点的天然气累计产量;根据PVT实验数据确定两相偏差系数、气油比数据以及两相偏差系数、气油比数据分别与地层压力的关系曲线;根据地层平均压力与露点压力的关系建立物质平衡方程,将所述的天然气累计产量及两相偏差系数代入所述物质平衡方程,计算各生产时间点的地层压力理论值,并根据实测静压数据校正地层压力理论值,生成地层压力校正值;根据所述的地层压力校正值及气油比数据生成对应的凝析油产量;以废弃地层压力为约束,根据所述物质平衡方程计算凝析油累计产量;根据所述的单井动态储量、凝析油累计产量及天然气累计产量分别计算凝析油采出程度及天然气采出程度。
在一实施例中,上述的获取单井中的标准参数,根据所述标准参数计算单井动态储量,包括:步骤a:获取凝析油产量及天然气产量,并根据所述凝析油产量计算气体当量;步骤b:根据所述天然气产量及气体当量计算凝析气产量;步骤c:根据区块平均单井天然气地质储量确定估算初值;步骤d:根据规整化拟时间公式及规整化拟压力公式分别计算规整化拟时间tca及规整化拟压力并生成关系曲线;步骤e:对所述关系曲线进行数据点线性回归,确定凝析气地质储量回归值;步骤f:判断所述估算初值与凝析气地质储量回归值之间的误差是否小于一预设值,如果是,则将所述凝析气地质储量回归值确定为所述的单井动态储量;否则返回执行所述步骤c~步骤f,直至所述估算初值与凝析气地质储量回归值之间的误差小于所述预设值。
在一实施例中,上述的的物质平衡方程为:
其中,G为气井凝析气地质储量,单位为108m3;Gp为已采出凝析气累计产量,单位为108m3;p为目前平均地层压力,单位为MPa;pi为原始地层压力,单位为MPa;Z为目前平均地层压力下的凝析气偏差系数,无量纲;Zi为原始地层压力下的凝析气偏差系数,无量纲;ZT为目前平均地层压力下的两相凝析气偏差系数,无量纲;ZTi为露点压力压力下的凝析气两相偏差系数,无量纲。
在一实施例中,上述的以废弃地层压力为约束,根据所述物质平衡方程计算凝析油累计产量,包括:根据井口外输压力下限确定井底流压,并通过以下公式确定所述废弃地层压力:其中,p为废弃地层压力,单位为MPa;pwf为井底流压,单位为MPa;q为气井废弃产量,单位为104m3;A、B为产能方程系数;以所述废弃地层压力为约束,根据所述物质平衡方程计算气井停止生产时间;根据所述气井停止时间及凝析油产量计算凝析油累计产量。
在一实施例中,计算所述凝析油采出程度的公式为:凝析油采出程度=凝析油累计产量/凝析油地质储量;计算所述天然气采出程度的公式为:天然气采出程度=天然气累计产量/天然气地质储量。
本发明实施例还提供一种缝洞型碳酸盐岩凝析气井开发指标分析装置,所述的分析装置包括:单井动态储量计算单元,用于获取单井中的标准参数,根据所述标准参数计算单井动态储量,所述单井动态储量包括:天然气地质储量和凝析油地质储量;天然气累计产量生成单元,用于根据预设的天然气产量分段递减规律预测天然气产量,并生成各生产时间点的天然气累计产量;PVT实验数据生成单元,用于根据PVT实验数据确定两相偏差系数、气油比数据以及两相偏差系数、气油比数据分别与地层压力的关系曲线;地层压力校正值生成单元,用于根据地层平均压力与露点压力的关系建立物质平衡方程,将所述的天然气累计产量及两相偏差系数代入所述物质平衡方程,计算各生产时间点的地层压力理论值,并根据实测静压数据生成地层压力校正值;凝析油产量生成单元,用于根据所述的地层压力校正值及气油比数据生成不同生产时间对应的凝析油产量;凝析油累计产量计算单元,用于以废弃地层压力为约束,根据所述物质平衡方程计算凝析油累计产量;采出程度计算单元,用于根据所述的单井动态储量、凝析油累计产量及天然气累计产量分别计算凝析油采出程度及天然气采出程度。
在一实施例中,上述的单井动态储量计算单元具体用于:步骤a:获取凝析油产量及天然气产量,并根据所述凝析油产量计算气体当量;步骤b:根据所述天然气产量及气体当量计算凝析气产量;步骤c:根据区块平均单井天然气地质储量确定估算初值;步骤d:根据规整化拟时间公式及规整化拟压力公式分别计算规整化拟时间tca及规整化拟压力并生成关系曲线;步骤e:对所述关系曲线进行数据点线性回归,确定凝析气地质储量回归值;步骤f:判断所述估算初值与凝析气地质储量回归值之间的误差是否小于一预设值,如果是,则将所述凝析气地质储量回归值确定为所述的单井动态储量;否则返回执行所述步骤c~步骤f,直至所述估算初值与凝析气地质储量回归值之间的误差小于所述预设值。
在一实施例中,上述的物质平衡方程为:
其中,G为气井凝析气地质储量,单位为108m3;Gp为已采出凝析气累计产量,单位为108m3;p为目前平均地层压力,单位为MPa;pi为原始地层压力,单位为MPa;Z为目前平均地层压力下的凝析气偏差系数,无量纲;Zi为原始地层压力下的凝析气偏差系数,无量纲;ZT为目前平均地层压力下的两相凝析气偏差系数,无量纲;ZTi为露点压力压力下的凝析气两相偏差系数,无量纲。
在一实施例中,上述的凝析油累计产量计算单元包括:废弃地层压力确定模块,用于根据井口外输压力下限确定井底流压,并通过以下公式确定所述废弃地层压力:其中,p为废弃地层压力,单位为MPa;pwf为井底流压,单位为MPa;q为气井废弃产量,单位为104m3;A、B为产能方程系数;气井停止生产时间计算模块,用于以所述废弃地层压力为约束,根据所述物质平衡方程计算气井停止生产时间;凝析油累计产量计算模块,用于根据所述气井停止时间及凝析油产量计算凝析油累计产量。
在一实施例中,计算所述凝析油采出程度的公式为:凝析油采出程度=凝析油累计产量/凝析油地质储量;计算所述天然气采出程度的公式为:天然气采出程度=天然气累计产量/天然气地质储量。
本发明实施例的有益效果在于,通过本发明,解决了常规分析方法难以处理的缝洞型碳酸盐岩凝析气井生产动态预测难题,可以准确预测凝析气井的凝析油、天然气产量,评价凝析油、天然气可采储量及采收率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例中的缝洞型碳酸盐岩凝析气井开发指标分析方法的流程图;
图2为本发明实施例中的单井动态储量计算曲线;
图3为本发明实施例中的单井天然气产量递减分析图;
图4为本发明实施例中的气油比、两相偏差系数与地层压力的关系曲线;
图5为本发明实施例中地层压力分析曲线;
图6为本发明实施例中凝析油分析曲线;
图7为本发明实施例中采出程度分析图;
图8为本发明实施例中的缝洞型碳酸盐岩凝析气井开发指标分析装置的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供一种缝洞型碳酸盐岩凝析气井开发指标分析方法及装置。以下结合附图对本发明进行详细说明。
本发明实施例提供一种缝洞型碳酸盐岩凝析气井开发指标分析方法,如图1所示,该缝洞型碳酸盐岩凝析气井开发指标分析方法主要包括以下各步骤:
步骤S101:获取单井中的标准参数,根据标准参数计算单井动态储量,单井动态储量包括:天然气地质储量和凝析油地质储量;
步骤S102:根据预设的天然气产量分段递减规律预测天然气产量,并预测未来的天然气累计产量;
步骤S103:根据PVT实验数据确定两相偏差系数、气油比数据以及两相偏差系数、气油比数据分别与地层压力的关系曲线;
步骤S104:根据地层平均压力与露点压力的关系建立物质平衡方程,将天然气累计产量及两相偏差系数代入物质平衡方程,计算各时间点的地层压力理论值,并根据实测静压数据校正地层压力理论值,生成地层压力校正值;
步骤S105:根据地层压力校正值及气油比数据生成对应的凝析油产量;
步骤S106:以废弃地层压力为约束,根据物质平衡方程计算凝析油累计产量;
步骤S107:根据单井动态储量、凝析油累计产量及天然气累计产量分别计算凝析油及天然气采出程度。
通过以上步骤S101至步骤S107,首先根据单井生产数据确定单井动态储量、确定天然气产量的递减规律,预测天然气的产量;然后根据凝析气藏物质平衡方程和PVT实验数据,预测凝析油的产量以及地层压力;进而可以得到可采储量和采出程度。本发明实施例的缝洞型碳酸盐岩凝析气井开发指标分析方法解决了常规分析方法难以处理的缝洞型碳酸盐岩凝析气井生产动态预测难题,可以准确预测凝析气井的凝析油、天然气产量,评价凝析油、天然气可采储量及采收率。
以下结合上述各步骤,对本发明实施例的缝洞型碳酸盐岩凝析气井开发指标分析方法做详细说明。
上述的步骤S101,获取单井中的标准参数,根据标准参数计算单井动态储量,单井动态储量包括:天然气地质储量和凝析油地质储量。
具体地,在该步骤S101中,首先获取凝析油产量和天然气产量等标准参数,然后,将凝析油产量代入公式(1),将凝析油产量折算为凝析油的气体当量:
其中,Ge为凝析油的气体当量,单位为m3/m3;γo为凝析油的相对密度,无量纲;TSC为标准状态下的温度,单位为K;pSC为标准状态下的压力,单位为MPa;Mo为凝析油的分子量,单位为kg/kmol。
进一步地,将天然气产量及公式(1)计算生成的气体当量Ge代入公式(2),计算凝析气产量:
其中,Rs为油气比,单位为m3/m3;qg为日产气量,单位为m3/d。
然后,按如下步骤计算单井动态储量:
1)根据区块平均单井天然气地质储量(OGIP)情况,确定估算初值;
2)计算规整化拟时间tca、规整化拟压力及绘图参数团即:
其中,Cti为原始条件下综合压缩系数,单位为MPa–1;q为日产量,单位为m3/d;Δpp为规整化拟压力差,单位为MPa;tca为规整化拟时间,单位为小时(hour);G为凝析气地质储量,单位为108m3。
3)根据上一步骤计算的规整化拟时间tca、规整化拟压力绘制曲线;
4)结合该曲线数据点线性回归,获得最佳拟合,直线段外推与x轴的交点即为凝析气地质储量G;
5)重复上述步骤1)至步骤4),直到估算初值与凝析气地质储量之间的误差小于一设定数值,则凝析气地质储量收敛,此时的凝析气地质储量即为所要求出的单井动态储量,计算结果如图2所示。此时可根据气油比将单井动态储量分解为天然气地质储量和凝析油地质储量。
上述的步骤S102:根据预设的天然气产量分段递减规律预测天然气产量,并生成各生产时间点的天然气累计产量。
在进行天然气产量统计时,可根据生产趋势分段确定递减规律,进而确定单井天然气产量的变化趋势,如图3所示。具体地,生产趋势可根据凝析气井的生产时间的长短来进行划分。如果凝析气井生产时间较长(生产时间大于330天),可用实际生产数据确定分段递减规律(如指数式或幂函数形式等);若凝析气井生产时间较短或未投产,可参考地质条件类似的邻井的递减规律来确定天然气产量,并结合生产时间来计算天然气累计产量。
上述步骤S103,根据PVT实验数据确定两相偏差系数、气油比数据以及两相偏差系数、气油比数据分别与地层压力的关系曲线。具体地,是根据定容衰竭实验,确定两相偏差系数;分析流体压力与体积关系及井流物组分,结合原始压力下偏差系数及测试压力,校正定容衰竭实验累计采出数据,从而得到气油比数据,并进一步确定两相偏差系数、气油比数据分别与地层压力的关系曲线,如图4所示。
上述步骤S104是对地层压力进行分析的步骤。当地层平均压力高于露点压力时,气藏中流体为单相气;当地层平均压力低于露点压力时,气藏中气液两相共存。利用累计产量等数据,进行物质平衡分析,确定地层压力。
其中,当地层平均压力高于露点压力时,可直接应用定容气藏物质平衡方程;当地层平均压力低于露点压力时,将定容气藏物质平衡方程中单相体积系数和偏差系数考虑为两相体积系数和两相偏差系数(ZT),该参数可以通过实验室确定。具体地,该物质平衡方程可通过公式(3)表示:
其中,G为气井凝析气地质储量,单位为108m3;Gp为已采出天然气累计产量,单位为108m3;p为目前平均地层压力,单位为MPa;pi为原始地层压力,单位为MPa;Z为目前平均地层压力下的天然气偏差系数,无量纲;Zi为原始地层压力下的天然气偏差系数,无量纲;ZT为目前平均地层压力下的两相天然气偏差系数,无量纲;ZTi为露点压力压力下的天然气两相偏差系数,无量纲。
确定了物质平衡方程后,将步骤S102中生成的天然气累计产量以及步骤S103中PVT实验提供的两相偏差系数代入该物质平衡方程,计算不同生产时间对应的地层压力(地层压力理论值);再根据实测静压数据对地层压力(理论值)进行校正,从而得到各生产时间点的地层压力(校正值),该值受实测地层压力约束,较为可靠,计算结果如图5所示。
在通过步骤S104计算出地层压力校正值后,执行上述步骤S105,根据该地层压力校正值以及步骤S103中生成的气油比数据预测未来的凝析油产量,如图6所示。
上述的步骤S106,以废弃地层压力为约束,根据物质平衡方程计算凝析油累计产量。具体地,根据井口外输压力下限,将井口压力通过井筒多相流理论折算为井底流压,进而根据公式(4)确定废弃地层压力:
其中,p为废弃地层压力,单位为MPa;pwf为井底流压,单位为MPa;q为气井废弃产量,单位为104m3;A、B为产能方程系数。
然后,以废弃地层压力为约束,根据上述的物质平衡方程计算气井停止生产的时间,进而根据该气井停止时间及凝析油产量计算凝析油累计产量。
上述的步骤S107,根据单井动态储量、凝析油累计产量及天然气累计产量分别计算凝析油采出程度及天然气采出程度。
在步骤S107中,以步骤S101计算得到的单井凝析油和天然气的动态储量为基础,根据步骤S102中预测生成的天然气累计产量,以及步骤S106中计算生成的凝析油累计产量,分别计算天然气采出程度和凝析油采出程度,如图7所示。具体地,计算凝析油采出程度的公式为:凝析油采出程度=凝析油累计产量/凝析油地质储量;计算天然气采出程度的公式为:天然气采出程度=天然气累计产量/天然气地质储量。
综上所述,根据本发明实施例的缝洞型碳酸盐岩凝析气井开发指标分析方法,能够分别分析得出缝洞型碳酸盐岩凝析气井凝析油产量、天然气产量、地层压力、可采储量及采收率等关键开发指标,从而为解决该类气藏开发指标分析难题提供有效途径。
本发明实施例还提供一种缝洞型碳酸盐岩凝析气井开发指标分析装置,如图8所示,该缝洞型碳酸盐岩凝析气井开发指标分析装置主要包括以下各部分:单井动态储量计算单元1、天然气累计产量生成单元2、PVT实验数据生成单元3、地层压力校正值生成单元4、凝析油产量生成单元5、凝析油累计产量计算单元6及采出程度计算单元7等。
其中,上述的单井动态储量计算单元1用于获取单井中的标准参数,根据标准参数计算单井动态储量,单井动态储量包括:天然气地质储量和凝析油地质储量;天然气累计产量生成单元2用于根据预设的天然气产量分段递减规律预测天然气产量,并生成各生产时间点的天然气累计产量;PVT实验数据生成单元3用于根据PVT实验数据确定两相偏差系数、气油比数据以及两相偏差系数、气油比数据分别与地层压力的关系曲线;地层压力校正值生成单元4用于根据地层平均压力与露点压力的关系建立物质平衡方程,将天然气累计产量及两相偏差系数代入物质平衡方程,计算各生产时间点的地层压力理论值,并根据实测静压数据生成地层压力校正值;凝析油产量生成单元5用于根据地层压力校正值及气油比数据生成不同生产时间对应的凝析油产量;凝析油累计产量计算单元6用于以废弃地层压力为约束,根据物质平衡方程计算凝析油累计产量;采出程度计算单元7用于根据单井动态储量、凝析油累计产量及天然气累计产量分别计算凝析油采出程度及天然气采出程度。
综上所述,首先根据单井生产数据确定单井动态储量、确定天然气产量的递减规律,预测未来天然气产量;然后根据凝析气藏物质平衡方程和PVT实验数据,预测凝析油的产量以及地层压力;进而可以得到可采储量和采收率。本发明实施例的缝洞型碳酸盐岩凝析气井开发指标分析装置解决了常规分析方法难以处理的缝洞型碳酸盐岩凝析气井生产动态预测难题,可以准确预测凝析气井的凝析油、天然气产量,评价凝析油、天然气可采储量及采收率。
以下结合上述各个部分,对本发明实施例的缝洞型碳酸盐岩凝析气井开发指标分析装置做详细说明。
上述的单井动态储量计算单元1,获取单井中的标准参数,根据标准参数计算单井动态储量,单井动态储量包括:天然气地质储量和凝析油地质储量。
具体地,该单井动态储量计算单元1首先获取凝析油产量和天然气产量等标准参数,然后,将凝析油产量代入公式(5),将凝析油产量折算为凝析油的气体当量:
其中,Ge为凝析油的气体当量,单位为m3/m3;γo为凝析油的相对密度,无量纲;TSC为标准状态下的温度,单位为K;pSC为标准状态下的压力,单位为MPa;Mo为凝析油的分子量,单位为kg/kmol。
进一步地,单井动态储量计算单元1将天然气产量及公式(5)计算生成的气体当量Ge代入公式(6),计算凝析气产量:
其中,Rs为油气比,单位为m3/m3;qg为日产气量,单位为m3/d。
然后,按如下步骤计算单井动态储量:
1)根据区块平均单井天然气地质储量(OGIP)情况,确定估算初值;
2)计算规整化拟时间tca、规整化拟压力及绘图参数团即:
其中,Cti为原始条件下综合压缩系数,单位为MPa–1;q为日产量,单位为m3/d;Δpp为规整化拟压力差,单位为MPa;tca为规整化拟时间,单位为小时(hour);G为凝析气地质储量,单位为108m3。
3)根据上一步骤计算的规整化拟时间tca、规整化拟压力绘制曲线;
4)结合该曲线数据点线性回归,获得最佳拟合,直线段外推与x轴的交点即为凝析气地质储量G;
5)重复上述步骤1)至步骤4),直到估算初值与凝析气地质储量之间的误差小于一设定数值,则凝析气地质储量收敛,此时的凝析气地质储量即为所要求出的单井动态储量,计算结果如图2所示。此时可根据气油比将单井动态储量分解为天然气地质储量和凝析油地质储量。
上述的天然气累计产量生成单元2用于根据预设的天然气产量分段递减规律预测天然气产量,并生成各生产时间点的天然气累计产量。
在进行天然气产量统计时,天然气累计产量生成单元2可根据生产趋势分段确定递减规律,进而确定单井天然气产量的变化趋势,如图3所示。具体地,生产趋势可根据凝析气井的生产时间的长短来进行划分。如果凝析气井生产时间较长(生产时间大于330天),可用实际生产数据确定分段递减规律(如指数式或幂函数形式等);若凝析气井生产较短或未投产,可参考地质条件类似的邻井的递减规律来确定天然气产量,并结合生产时间来计算天然气累计产量。
上述PVT实验数据生成单元3则根据PVT实验数据确定两相偏差系数、气油比数据以及两相偏差系数、气油比数据分别与地层压力的关系曲线。具体地,是根据定容衰竭实验,确定两相偏差系数;分析流体压力与体积关系及井流物组分,结合原始压力下偏差系数及测试压力,校正定容衰竭实验累计采出数据,从而得到气油比数据,并进一步确定两相偏差系数、气油比数据分别与地层压力的关系曲线,如图4所示。
上述地层压力校正值生成单元4是用以对地层压力进行分析。当地层平均压力高于露点压力时,气藏中流体为单相气;当地层平均压力低于露点压力时,气藏中气液两相共存。利用累计产量等数据,进行物质平衡分析确定地层压力。
其中,当地层平均压力高于露点压力时,可直接应用定容气藏物质平衡方程;当地层平均压力低于露点压力时,将定容气藏物质平衡方程中单相体积系数和偏差系数考虑为两相体积系数和两相偏差系数(ZT),该参数可以通过实验室确定。具体地,该物质平衡方程可通过公式(7)表示:
其中,G为气井凝析气地质储量,单位为108m3;Gp为已采出天然气累计产量,单位为108m3;p为目前平均地层压力,单位为MPa;pi为原始地层压力,单位为MPa;Z为目前平均地层压力下的天然气偏差系数,无量纲;Zi为原始地层压力下的天然气偏差系数,无量纲;ZT为目前平均地层压力下的两相天然气偏差系数,无量纲;ZTi为露点压力压力下的天然气两相偏差系数,无量纲。
确定了物质平衡方程后,将通过天然气累计产量生成单元2生成的天然气累计产量以及PVT实验数据生成单元3中PVT实验提供的两相偏差系数代入该物质平衡方程,计算不同生产时间对应的地层压力(地层压力理论值);再根据实测静压数据对地层压力(理论值)进行校正,从而得到各生产时间点的地层压力(校正值),该值受实测地层压力约束,较为可靠,计算结果如图5所示。
在通过地层压力校正值生成单元4计算出地层压力校正值后,执行上述凝析油产量生成单元5,根据该地层压力校正值以及步骤S103中生成的气油比数据预测未来的凝析油产量,如图6所示。
上述的累积产量计算单元6,以废弃地层压力为约束,根据物质平衡方程计算凝析油累计产量。具体地,根据井口外输压力下限,将井口压力通过井筒多相流理论折算为井底流压,进而根据公式(8)确定废弃地层压力:
其中,p为废弃地层压力,单位为MPa;pwf为井底流压,单位为MPa;q为气井废弃产量,单位为104m3;A、B为产能方程系数。
然后,以废弃地层压力为约束,根据上述的物质平衡方程计算气井停止生产的时间,进而根据该气井停止时间及凝析油产量计算凝析油累计产量。
上述的采出程度计算单元7,根据单井动态储量、凝析油累计产量及天然气累计产量分别计算凝析油采出程度及天然气采出程度。
采出程度计算单元7,以单井动态储量计算单元1计算得到的单井凝析油和天然气的动态储量为基础,根据天然气累计产量生成单元2中预测生成的天然气累计产量,以及累积产量计算单元6中计算生成的凝析油累计产量,分别计算天然气采出程度和凝析油采出程度,如图7所示。具体地,计算凝析油采出程度的公式为:凝析油采出程度=凝析油累计产量/凝析油地质储量;计算天然气采出程度的公式为:天然气采出程度=天然气累计产量/天然气地质储量。
综上所述,根据本发明实施例的缝洞型碳酸盐岩凝析气井开发指标分析装置,能够分析得出缝洞型碳酸盐岩凝析气井凝析油产量、天然气产量、地层压力、可采储量及采收率等关键开发指标,从而为解决该类气藏开发指标分析难题提供有效途径。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤,可以通过程序来指令相关的硬件来完成,该程序可以存储于一计算机可读取存储介质中,比如ROM/RAM、磁碟、光盘等。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种缝洞型碳酸盐岩凝析气井开发指标分析方法,其特征在于,所述的分析方法包括:
获取单井中的标准参数,根据所述标准参数计算单井动态储量,所述单井动态储量包括:天然气地质储量和凝析油地质储量;
根据预设的天然气产量分段递减规律确定天然气产量以及各生产时间点的天然气累计产量;
根据PVT实验数据确定两相偏差系数、气油比数据以及两相偏差系数、气油比数据分别与地层压力的关系曲线;
根据地层平均压力与露点压力的关系建立物质平衡方程,将所述的天然气累计产量及两相偏差系数代入所述物质平衡方程,计算各生产时间点的地层压力理论值,并根据实测静压数据校正计算结果及所述的地层压力理论值生成地层压力校正值;
根据所述的地层压力校正值及气油比数据生成对应的凝析油产量;
以废弃地层压力为约束,根据所述物质平衡方程计算凝析油累计产量;
根据所述的单井动态储量、凝析油累计产量及天然气累计产量分别计算凝析油采出程度及天然气采出程度。
2.根据权利要求1所述的缝洞型碳酸盐岩凝析气井开发指标分析方法,其特征在于,获取单井中的标准参数,根据所述标准参数计算单井动态储量,包括:
步骤a:获取凝析油产量及天然气产量,并根据所述凝析油产量计算气体当量;
步骤b:根据所述天然气产量及气体当量计算凝析气产量;
步骤c:根据区块平均单井天然气地质储量确定估算初值;
步骤d:根据规整化拟时间公式及规整化拟压力公式分别计算规整化拟时间tca及规整化拟压力并生成关系曲线;
步骤e:对所述关系曲线进行数据点线性回归,确定凝析气地质储量回归值;
步骤f:判断所述估算初值与凝析气地质储量回归值之间的误差是否小于一预设值,如果是,则将所述凝析气地质储量回归值确定为所述的单井动态储量;否则返回执行所述步骤c~步骤f,直至所述估算初值与凝析气地质储量回归值之间的误差小于所述预设值。
3.根据权利要求1所述的缝洞型碳酸盐岩凝析气井开发指标分析方法,其特征在于,所述的物质平衡方程为:
其中,G为气井凝析气地质储量,单位为108m3;Gp为已采出凝析气累计产量,单位为108m3;p为目前平均地层压力,单位为MPa;pi为原始地层压力,单位为MPa;Z为目前平均地层压力下的凝析气偏差系数,无量纲;Zi为原始地层压力下的凝析气偏差系数,无量纲;ZT为目前平均地层压力下的两相凝析气偏差系数,无量纲;ZTi为露点压力压力下的凝析气两相偏差系数,无量纲。
4.根据权利要求1所述的缝洞型碳酸盐岩凝析气井开发指标分析方法,其特征在于,以废弃地层压力为约束,根据所述物质平衡方程计算凝析油累计产量,包括:
根据井口外输压力下限确定井底流压,并通过以下公式确定所述废弃地层压力:
其中,p为废弃地层压力,单位为MPa;pwf为井底流压,单位为MPa;q为气井废弃产量,单位为104m3;A、B为产能方程系数;
以所述废弃地层压力为约束,根据所述物质平衡方程计算气井停止生产时间;
根据所述气井停止时间及凝析油产量计算凝析油累计产量。
5.根据权利要求1所述的缝洞型碳酸盐岩凝析气井开发指标分析方法,其特征在于,计算所述凝析油采出程度的公式为:凝析油采出程度=凝析油累计产量/凝析油地质储量;
计算所述天然气采出程度的公式为:天然气采出程度=天然气累计产量/天然气地质储量。
6.一种缝洞型碳酸盐岩凝析气井开发指标分析装置,其特征在于,所述的分析装置包括:
单井动态储量计算单元,用于获取单井中的标准参数,根据所述标准参数计算单井动态储量,所述单井动态储量包括:天然气地质储量和凝析油地质储量;
天然气累计产量生成单元,用于根据预设的天然气产量分段递减规律确定天然气产量以及各生产时间点的天然气累计产量;
PVT实验数据生成单元,用于根据PVT实验数据确定两相偏差系数、气油比数据以及两相偏差系数、气油比数据分别与地层压力的关系曲线;
地层压力校正值生成单元,用于根据地层平均压力与露点压力的关系建立物质平衡方程,将所述的天然气累计产量及两相偏差系数代入所述物质平衡方程,计算各生产时间点的地层压力理论值,并根据实测静压数据生成地层压力校正值;
凝析油产量生成单元,用于根据所述的地层压力校正值及气油比数据生成不同生产时间对应的凝析油产量;
凝析油累计产量计算单元,用于以废弃地层压力为约束,根据所述物质平衡方程计算凝析油累计产量;
采出程度计算单元,用于根据所述的单井动态储量、凝析油累计产量及天然气累计产量分别计算凝析油采出程度及天然气采出程度。
7.根据权利要求6所述的缝洞型碳酸盐岩凝析气井开发指标分析装置,其特征在于,所述的单井动态储量计算单元具体用于:
步骤a:获取凝析油产量及天然气产量,并根据所述凝析油产量计算气体当量;
步骤b:根据所述天然气产量及气体当量计算凝析气产量;
步骤c:根据区块平均单井天然气地质储量确定估算初值;
步骤d:根据规整化拟时间公式及规整化拟压力公式分别计算规整化拟时间tca及规整化拟压力并生成关系曲线;
步骤e:对所述关系曲线进行数据点线性回归,确定凝析气地质储量回归值;
步骤f:判断所述估算初值与凝析气地质储量回归值之间的误差是否小于一预设值,如果是,则将所述凝析气地质储量回归值确定为所述的单井动态储量;否则返回执行所述步骤c~步骤f,直至所述估算初值与凝析气地质储量回归值之间的误差小于所述预设值。
8.根据权利要求6所述的缝洞型碳酸盐岩凝析气井开发指标分析装置,其特征在于,所述的物质平衡方程为:
其中,G为气井凝析气地质储量,单位为108m3;Gp为已采出凝析气累计产量,单位为108m3;p为目前平均地层压力,单位为MPa;pi为原始地层压力,单位为MPa;Z为目前平均地层压力下的凝析气偏差系数,无量纲;Zi为原始地层压力下的凝析气偏差系数,无量纲;ZT为目前平均地层压力下的两相凝析气偏差系数,无量纲;ZTi为露点压力压力下的凝析气两相偏差系数,无量纲。
9.根据权利要求6所述的缝洞型碳酸盐岩凝析气井开发指标分析装置,其特征在于,所述的凝析油累计产量计算单元包括:
废弃地层压力确定模块,用于根据井口外输压力下限确定井底流压,并通过以下公式确定所述废弃地层压力:
其中,p为废弃地层压力,单位为MPa;pwf为井底流压,单位为MPa;q为气井废弃产量,单位为104m3;A、B为产能方程系数;
气井停止生产时间计算模块,用于以所述废弃地层压力为约束,根据所述物质平衡方程计算气井停止生产时间;
凝析油累计产量计算模块,用于根据所述气井停止时间及凝析油产量计算凝析油累计产量。
10.根据权利要求6所述的缝洞型碳酸盐岩凝析气井开发指标分析装置,其特征在于,计算所述凝析油采出程度的公式为:凝析油采出程度=凝析油累计产量/凝析油地质储量;
计算所述天然气采出程度的公式为:天然气采出程度=天然气累计产量/天然气地质储量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510556736.1A CN105134196B (zh) | 2015-09-02 | 2015-09-02 | 一种缝洞型碳酸盐岩凝析气井开发指标分析方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510556736.1A CN105134196B (zh) | 2015-09-02 | 2015-09-02 | 一种缝洞型碳酸盐岩凝析气井开发指标分析方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105134196A true CN105134196A (zh) | 2015-12-09 |
CN105134196B CN105134196B (zh) | 2018-05-04 |
Family
ID=54719733
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510556736.1A Active CN105134196B (zh) | 2015-09-02 | 2015-09-02 | 一种缝洞型碳酸盐岩凝析气井开发指标分析方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105134196B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108661631A (zh) * | 2017-04-01 | 2018-10-16 | 中国石油化工股份有限公司 | 一种页岩气井产量预测方法 |
CN108664679A (zh) * | 2017-04-01 | 2018-10-16 | 中国石油化工股份有限公司 | 一种油气井生产数据分析方法 |
CN110058000A (zh) * | 2018-01-15 | 2019-07-26 | 中国石油天然气股份有限公司 | 原油体积系数的获取方法及装置 |
CN110905484A (zh) * | 2018-09-14 | 2020-03-24 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏井间连通程度计算方法 |
CN111324938A (zh) * | 2018-12-13 | 2020-06-23 | 中国石油天然气股份有限公司 | 基质供气能力的计算方法及系统 |
CN111648768A (zh) * | 2020-06-16 | 2020-09-11 | 陕西延长石油(集团)有限责任公司 | 一种确定致密气藏产水气井合理产量制度的方法 |
CN111680814A (zh) * | 2019-03-11 | 2020-09-18 | 中国石油化工股份有限公司 | 一种碳酸盐岩缝洞型油藏弹性驱动合理配产优化方法 |
CN111852462A (zh) * | 2019-04-29 | 2020-10-30 | 中国石油天然气股份有限公司 | 油井动态储量获取方法和装置 |
CN112324420A (zh) * | 2019-07-30 | 2021-02-05 | 中国石油天然气股份有限公司 | 井控地质储量的预测方法及系统 |
CN112343573A (zh) * | 2019-08-07 | 2021-02-09 | 中国石油天然气股份有限公司 | 凝析气井产量的预测方法 |
CN112395724A (zh) * | 2019-07-30 | 2021-02-23 | 中国石油天然气股份有限公司 | 预测水力压裂地层效果的方法及系统 |
CN113323656A (zh) * | 2021-06-16 | 2021-08-31 | 中海石油(中国)有限公司 | 封闭凝析气藏的开发指标预测方法及计算机可读存储介质 |
CN113622902A (zh) * | 2020-04-21 | 2021-11-09 | 中国石油天然气股份有限公司 | 凝析气藏地层天然气密度确定方法及装置 |
CN116127874A (zh) * | 2023-02-15 | 2023-05-16 | 西南石油大学 | 一种凝析井地层压力及单井控制储量联合一体的评价方法 |
CN116383573A (zh) * | 2023-03-20 | 2023-07-04 | 中海石油(中国)有限公司海南分公司 | 一种基于多区相变传质渗流耦合的凝析气产能评价方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0482911A2 (en) * | 1990-10-26 | 1992-04-29 | Halliburton Company | Evaluation of fluid loss for subsurface fracturing |
CN101942984A (zh) * | 2010-01-27 | 2011-01-12 | 中国石油化工股份有限公司 | 碳酸盐岩缝洞型油藏注水替油开采方法 |
CN102508317A (zh) * | 2011-11-09 | 2012-06-20 | 中国石油天然气股份有限公司 | 一种识别缝洞型碳酸盐岩储层流体性质的方法 |
CN103413030A (zh) * | 2013-07-24 | 2013-11-27 | 中国石油天然气股份有限公司 | 一种缝洞型碳酸盐岩气藏动态分析方法及系统 |
CN104699954A (zh) * | 2015-01-30 | 2015-06-10 | 中国地质大学(武汉) | 一种岩溶缝洞型碳酸盐岩油藏采收率的标定方法 |
-
2015
- 2015-09-02 CN CN201510556736.1A patent/CN105134196B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0482911A2 (en) * | 1990-10-26 | 1992-04-29 | Halliburton Company | Evaluation of fluid loss for subsurface fracturing |
CN101942984A (zh) * | 2010-01-27 | 2011-01-12 | 中国石油化工股份有限公司 | 碳酸盐岩缝洞型油藏注水替油开采方法 |
CN102508317A (zh) * | 2011-11-09 | 2012-06-20 | 中国石油天然气股份有限公司 | 一种识别缝洞型碳酸盐岩储层流体性质的方法 |
CN103413030A (zh) * | 2013-07-24 | 2013-11-27 | 中国石油天然气股份有限公司 | 一种缝洞型碳酸盐岩气藏动态分析方法及系统 |
CN104699954A (zh) * | 2015-01-30 | 2015-06-10 | 中国地质大学(武汉) | 一种岩溶缝洞型碳酸盐岩油藏采收率的标定方法 |
Non-Patent Citations (1)
Title |
---|
郭平等: "缝洞型凝析气藏衰竭开采动态实验研究", 《石油钻采工艺》 * |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108664679B (zh) * | 2017-04-01 | 2021-07-27 | 中国石油化工股份有限公司 | 一种油气井生产数据分析方法 |
CN108664679A (zh) * | 2017-04-01 | 2018-10-16 | 中国石油化工股份有限公司 | 一种油气井生产数据分析方法 |
CN108661631A (zh) * | 2017-04-01 | 2018-10-16 | 中国石油化工股份有限公司 | 一种页岩气井产量预测方法 |
CN110058000A (zh) * | 2018-01-15 | 2019-07-26 | 中国石油天然气股份有限公司 | 原油体积系数的获取方法及装置 |
CN110058000B (zh) * | 2018-01-15 | 2022-03-01 | 中国石油天然气股份有限公司 | 原油体积系数的获取方法及装置 |
CN110905484A (zh) * | 2018-09-14 | 2020-03-24 | 中国石油化工股份有限公司 | 一种缝洞型碳酸盐岩油藏井间连通程度计算方法 |
CN111324938A (zh) * | 2018-12-13 | 2020-06-23 | 中国石油天然气股份有限公司 | 基质供气能力的计算方法及系统 |
CN111324938B (zh) * | 2018-12-13 | 2022-11-04 | 中国石油天然气股份有限公司 | 基质供气能力的计算方法及系统 |
CN111680814B (zh) * | 2019-03-11 | 2023-04-07 | 中国石油化工股份有限公司 | 一种碳酸盐岩缝洞型油藏弹性驱动合理配产优化方法 |
CN111680814A (zh) * | 2019-03-11 | 2020-09-18 | 中国石油化工股份有限公司 | 一种碳酸盐岩缝洞型油藏弹性驱动合理配产优化方法 |
CN111852462A (zh) * | 2019-04-29 | 2020-10-30 | 中国石油天然气股份有限公司 | 油井动态储量获取方法和装置 |
CN112395724A (zh) * | 2019-07-30 | 2021-02-23 | 中国石油天然气股份有限公司 | 预测水力压裂地层效果的方法及系统 |
CN112324420B (zh) * | 2019-07-30 | 2023-07-25 | 中国石油天然气股份有限公司 | 井控地质储量的预测方法及系统 |
CN112324420A (zh) * | 2019-07-30 | 2021-02-05 | 中国石油天然气股份有限公司 | 井控地质储量的预测方法及系统 |
CN112343573A (zh) * | 2019-08-07 | 2021-02-09 | 中国石油天然气股份有限公司 | 凝析气井产量的预测方法 |
CN112343573B (zh) * | 2019-08-07 | 2023-09-26 | 中国石油天然气股份有限公司 | 凝析气井产量的预测方法 |
CN113622902A (zh) * | 2020-04-21 | 2021-11-09 | 中国石油天然气股份有限公司 | 凝析气藏地层天然气密度确定方法及装置 |
CN113622902B (zh) * | 2020-04-21 | 2023-09-26 | 中国石油天然气股份有限公司 | 凝析气藏地层天然气密度确定方法及装置 |
CN111648768B (zh) * | 2020-06-16 | 2023-03-17 | 陕西延长石油(集团)有限责任公司 | 一种确定致密气藏产水气井合理产量制度的方法 |
CN111648768A (zh) * | 2020-06-16 | 2020-09-11 | 陕西延长石油(集团)有限责任公司 | 一种确定致密气藏产水气井合理产量制度的方法 |
CN113323656A (zh) * | 2021-06-16 | 2021-08-31 | 中海石油(中国)有限公司 | 封闭凝析气藏的开发指标预测方法及计算机可读存储介质 |
CN116127874A (zh) * | 2023-02-15 | 2023-05-16 | 西南石油大学 | 一种凝析井地层压力及单井控制储量联合一体的评价方法 |
CN116127874B (zh) * | 2023-02-15 | 2023-12-19 | 西南石油大学 | 一种凝析井地层压力及单井控制储量联合一体的评价方法 |
CN116383573A (zh) * | 2023-03-20 | 2023-07-04 | 中海石油(中国)有限公司海南分公司 | 一种基于多区相变传质渗流耦合的凝析气产能评价方法 |
CN116383573B (zh) * | 2023-03-20 | 2023-10-10 | 中海石油(中国)有限公司海南分公司 | 一种基于多区相变传质渗流耦合的凝析气产能评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105134196B (zh) | 2018-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105134196A (zh) | 一种缝洞型碳酸盐岩凝析气井开发指标分析方法及装置 | |
O'sullivan et al. | Fluid and heat flow in gas-rich geothermal reservoirs | |
Fevang | Gas condensate flow behavior and sampling | |
CN105464652A (zh) | 缝洞型碳酸盐岩油藏缝洞单元动态储量计算方法和系统 | |
CN105569646A (zh) | 一种油气井技术可采储量预测方法 | |
CN107563899A (zh) | 油气井产能预测方法及装置 | |
Andersen | Analytical modeling and correction of steady state relative permeability experiments with capillary end effects–An improved intercept method, scaling and general capillary numbers | |
Payne | Material-balance calculations in tight-gas reservoirs: the pitfalls of P/Z plots and a more accurate technique | |
CN107130955B (zh) | 井底流压的确定方法及储集体天然能量的确定方法 | |
Czarnota et al. | Semianalytical horizontal well length optimization under pseudosteady-state conditions | |
Douglas et al. | A data assimilation enabled model for coupling dual porosity flow with free flow | |
CN106481315A (zh) | 陆上砂岩油藏单井可采储量快速确定模型及建立方法 | |
Thompson et al. | Multiphase flowing material balance without relative permeability curves | |
Mott | Calculating well deliverability in gas condensate reservoirs | |
RU2569522C1 (ru) | Способ определения давления в скважине | |
Ojo et al. | Dynamic material balance equation and solution technique using production and PVT data | |
Gunn et al. | Principles for wellbore simulator validation and calibration using matching analysis—I. Analytical techniques | |
Ballin et al. | Cupiagua: modeling of a complex fractured reservoir using compositional upscaling | |
Khalid et al. | Inflow performance relationship for horizontal wells producing from multi-layered heterogeneous solution gas-drive reservoirs | |
Ballin et al. | Cupiagua: A complex full-field fractured reservoir study using compositional upscaling | |
Vardcharragosad et al. | Linear vs. Radial boundary-dominated flow: Implications for gas-well-decline analysis | |
Franz et al. | Calculation of geothermal stored heat from a numerical model for reserve estimation | |
Baig et al. | Productivity Increase Estimation for Multi Stage Fracturing in Horizontal Wells for Tight Oil Reservoirs | |
Zheng | Integrated production data analysis of horizontal fractured well in unconventional reservoir | |
CN113496301B (zh) | 油气田资产评价方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |