CN115860266A - 一种页岩气/煤层气井产能评价方法、***及电子设备 - Google Patents
一种页岩气/煤层气井产能评价方法、***及电子设备 Download PDFInfo
- Publication number
- CN115860266A CN115860266A CN202310119422.XA CN202310119422A CN115860266A CN 115860266 A CN115860266 A CN 115860266A CN 202310119422 A CN202310119422 A CN 202310119422A CN 115860266 A CN115860266 A CN 115860266A
- Authority
- CN
- China
- Prior art keywords
- gas
- zone
- fracture
- matrix
- well
- 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
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种页岩气/煤层气井产能评价方法、***及电子设备,涉及页岩气/煤层气产能评价技术领域,该方法包括根据目标井的储层参数和钻完井参数建立目标水平井仿真几何模型;根据上述模型构建控制方程集合及其对应的初始条件和边界条件,并在初始条件和边界条件的限制下,求解控制方程集合,计算目标井在生产过程中的产能评价参数;将实测产能评价参数与计算得到的产能评价参数进行拟合,并将最优拟合结果对应的关键参数和产能评价参数分别确定为目标井在生产过程中的最优关键参数和最终产能评价参数。本发明能够准确获取页岩气/煤层气井的最终可采储量、日产气中吸附气/游离气比例、裂缝导流能力、基质扩散系数和吸附能力参数。
Description
技术领域
本发明涉及页岩气/煤层气产能评价技术领域,特别是涉及一种基于同位素分馏的页岩气/煤层气井产能评价方法、***及电子设备。
背景技术
由于页岩气储层和深层煤层气储层低渗致密的特性,通常不存在天然产能,需要借助水平钻井和水力压裂等措施来实现页岩储层的经济高效开采。因此,页岩气/煤层气井产能评价和预测以及生产过程吸附气/游离气比例动用规律,由于事关页岩气/煤层气开发方案设计、压裂评价、生产动态分析、开发方案的优化和调整,一直以来受到了学术界和工业界的广泛关注和重视。
传统的产能分析方法包括经验产量递减分析方法、解析求解方法和数值模拟方法。经验产量递减分析方法包括Arps递减、改进双曲递减、扩展指数递减、Duong递减、改进Duong递减、幂指数递减和改进幂指数递减等产量递减分析方法,但不同经验模型均有各自的适用流态,只能对井底流动压力恒定或近似恒定条件下的生产数据进行分析,且要求生产数据连续稳定避免长时间关井。目前中国页岩气/煤层气勘探开发正处于起步阶段,受地表条件及地面管网限制,气井难以实现井底压力的恒定生产方式,多数采用前期变产量变压力的生产方式,因此经验产量递减分析方法在国内页岩气/煤层气井中的应用受到了很大程度的限制。
相较于经验产量递减分析方法,解析求解方法和数值模拟方法具备更加科学、严谨的数学物理基础。解析求解方法是基于基质孔隙和裂缝耦合的气体渗流机理,通过简化求解得到储层压力或产气量解析解的方法,而数值模拟方法则是考虑储层多种流动机理(渗流、扩散、吸附-解吸)建立控制方程,通过商业化或自研的数值模拟软件来进行数值模拟和求解的方法。解析求解方法和数值模拟方法具有更好的拟合效果,但缺点在于操作复杂,待标定的模型参数较多且多解性强,模型准确性大多依赖于历史拟合确定的模型参数是否准确。目前,用于标定的数据主要为产量或生产压力等较为单一的数据,缺乏更多元、对生产状态更敏感的指标来进行约束或验证。此外,现有产能评价方法的主要功能为产量和最终可采储量预测,无法对生产过程吸附气、游离气的动用比例进行定量解析。
近年来的研究发现,页岩气/煤层气井生产过程中存在显著的甲烷碳同位素分馏,并且同位素分馏特征与采收率、总含气量、吸附气/游离气比例、岩石渗透率等因素密切相关,这为利用同位素指标评价页岩气/煤层气井产能以及生产过程吸附气/游离气动用比例提供了全新的思路和途径。然而,目前有关页岩气/煤层气产出过程同位素分馏的研究多停留在定性认识和半定量分析阶段,定量化表征模型相对较少且主要针对均质岩心的气体解析过程,难以推广应用于多尺度多机理耦合的页岩气/煤层气生产过程,这很大程度上限制了对同位素分馏机理的深入认识,以及同位素分馏方法在页岩气/煤层气勘探开发定量评价技术领域的应用。
发明内容
鉴于此,本发明的目的是提供一种页岩气/煤层气井产能评价方法、***及电子设备。
为实现上述目的,本发明提供了如下方案:
第一方面,本发明提供了一种页岩气/煤层气井产能评价方法,包括:
获取目标井的储层参数和钻完井参数,并根据所述储层参数和所述钻完井参数,建立目标水平井仿真几何模型;所述目标为页岩气或煤层气,所述储层参数包括岩石密度、地层温度、储层初始气体压力、初始气体同位素值和含气孔隙度;所述钻完井参数包括水平井长度、压裂段数、总簇数、压裂缝间距、半缝长、压裂缝高度和压裂缝缝宽;所述目标水平井仿真几何模型包括裂缝区和基质区;所述基质区包括压裂改造区以及未压裂改造区;
根据所述目标水平井仿真几何模型,构建控制方程集合;所述控制方程集合包括多个控制方程组,分别为第一控制方程组、第二控制方程组和第三控制方程组;所述第一控制方程组包括目标井生产过程中裂缝区内12CH4气体流动控制方程和13CH4气体流动控制方程;所述第二控制方程组包括能够表征目标井生产过程中压裂改造区基质孔隙内游离态12CH4、游离态13CH4、吸附态12CH4和吸附态13CH4耦合流动的控制方程;所述第三控制方程组包括能够表征目标井生产过程中未压裂改造区基质孔隙内游离态12CH4、游离态13CH4、吸附态12CH4和吸附态13CH4耦合流动的控制方程;所述控制方程集合用于根据目标井的关键参数确定目标井的产能评价参数;所述关键参数包括裂缝导流能力、压裂改造区基质孔隙内12CH4的气体扩散系数、压裂改造区基质孔隙内13CH4的气体扩散系数、未压裂改造区基质孔隙内12CH4的气体扩散系数、未压裂改造区基质孔隙内13CH4的气体扩散系数、基质朗格缪尔体积和基质朗格缪尔压力;所述产能评价参数包括游离气产气量、吸附气产气量、日产气量、吸附气比例、游离气比例以及日产气表观同位素值;
确定所述控制方程集合的初始条件和边界条件;所述控制方程集合的初始条件是根据储层初始气体压力和初始气体同位素值确定的,所述控制方程集合的边界条件是根据实测井底流压和日产气同位素值确定的;所述控制方程集合的初始条件包括第一初始条件、第二初始条件和第三初始条件;所述第一初始条件为裂缝区内12CH4气体流动和13CH4气体流动的初始条件;所述第二初始条件为压裂改造区内12CH4气体流动和13CH4气体流动的初始条件;所述第三初始条件为未压裂改造区内12CH4气体流动和13CH4气体流动的初始条件;所述控制方程集合的边界条件包括控制边界条件、连续性边界条件和绝缘边界条件;
在所述控制方程集合的初始条件和边界条件的限制下,求解所述控制方程集合,计算得到目标井在生产过程中任意时刻的产能评价参数;
获取目标井在生产过程中任意时刻实测的日产气量和日产气表观同位素值,并将实测的日产气量和日产气表观同位素值,与计算得到的日产气量和日产气表观同位素值进行拟合,得到拟合结果;
将最优拟合结果对应的关键参数确定为目标井在生产过程中的最优关键参数,将所述最优拟合结果对应的产能评价参数确定为目标井在生产过程中最终的产能评价参数;所述最优拟合结果为拟合误差最小的拟合结果。
第二方面,一种页岩气/煤层气井产能评价***,包括:
目标水平井仿真几何模型建立模块,用于获取目标井的储层参数和钻完井参数,并根据所述储层参数和所述钻完井参数,建立目标水平井仿真几何模型;所述目标为页岩气或煤层气,所述储层参数包括岩石密度、地层温度、储层初始气体压力、初始气体同位素值和含气孔隙度;所述钻完井参数包括水平井长度、压裂段数、总簇数、压裂缝间距、半缝长、压裂缝高度和压裂缝缝宽;所述目标水平井仿真几何模型包括裂缝区和基质区;所述基质区包括压裂改造区以及未压裂改造区;
控制方程集合构建模块,用于根据所述目标水平井仿真几何模型,构建控制方程集合;所述控制方程集合包括多个控制方程组,分别为第一控制方程组、第二控制方程组和第三控制方程组;所述第一控制方程组包括目标井生产过程中裂缝区内12CH4气体流动控制方程和13CH4气体流动控制方程;所述第二控制方程组包括能够表征目标井生产过程中压裂改造区基质孔隙内游离态12CH4、游离态13CH4、吸附态12CH4和吸附态13CH4耦合流动的控制方程;所述第三控制方程组包括能够表征目标井生产过程中未压裂改造区基质孔隙内游离态12CH4、游离态13CH4、吸附态12CH4和吸附态13CH4耦合流动的控制方程;所述控制方程集合用于根据目标井的关键参数确定目标井的产能评价参数;所述关键参数包括裂缝导流能力、压裂改造区基质孔隙内12CH4的气体扩散系数、压裂改造区基质孔隙内13CH4的气体扩散系数、未压裂改造区基质孔隙内12CH4的气体扩散系数、未压裂改造区基质孔隙内13CH4的气体扩散系数、基质朗格缪尔体积和基质朗格缪尔压力;所述产能评价参数包括游离气产气量、吸附气产气量、日产气量、吸附气比例、游离气比例以及日产气表观同位素值;
初始条件和边界条件确定模块,用于确定所述控制方程集合的初始条件和边界条件;所述控制方程集合的初始条件是根据储层初始气体压力和初始气体同位素值确定的,所述控制方程集合的边界条件是根据实测井底流压和日产气同位素值确定的;所述控制方程集合的初始条件包括第一初始条件、第二初始条件和第三初始条件;所述第一初始条件为裂缝区内12CH4气体流动和13CH4气体流动的初始条件;所述第二初始条件为压裂改造区内12CH4气体流动和13CH4气体流动的初始条件;所述第三初始条件为未压裂改造区内12CH4气体流动和13CH4气体流动的初始条件;所述控制方程集合的边界条件包括控制边界条件、连续性边界条件和绝缘边界条件;
任意时刻产能评价参数计算模块,用于在所述控制方程集合的初始条件和边界条件的限制下,求解所述控制方程集合,计算得到目标井在生产过程中任意时刻的产能评价参数;
拟合结果确定模块,用于获取目标井在生产过程中任意时刻实测的日产气量和日产气表观同位素值,并将实测的日产气量和日产气表观同位素值,与计算得到的日产气量和日产气表观同位素值进行拟合,得到拟合结果;
最优关键参数和最终产能评价参数确定模块,用于将最优拟合结果对应的关键参数确定为目标井在生产过程中的最优关键参数,将所述最优拟合结果对应的产能评价参数确定为目标井在生产过程中最终的产能评价参数;所述最优拟合结果为拟合误差最小的拟合结果。
第三方面,本发明提供了一种电子设备,包括存储器及处理器,所述存储器用于存储计算机程序,所述处理器运行所述计算机程序以使所述电子设备执行根据第一方面所述的一种页岩气/煤层气井产能评价方法。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供了一种基于同位素分馏的页岩气/煤层气井产能评价方法、***及电子设备,通过建立页岩气/煤层气井生产过程多尺度(人工裂缝、天然裂缝、基质孔隙)多机理(渗流、扩散和吸附-解吸)耦合的同位素分馏定量表征模型,提出单一机理的解耦方法,为揭示页岩气/煤层气井产出过程同位素分馏特征、影响因素及机理提供理论支撑。同时,通过标定页岩气/煤层气井实测生产数据,优化同位素分馏模型参数,进而为评价和预测页岩气/煤层气井产能、确定页岩气/煤层气井生产过程吸附气/游离气动用规律及最终动用比例提供技术手段。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的页岩气/煤层气井产能评价方法的流程示意图;
图2为本发明实施例提供的页岩气/煤层气井产能评价***的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例一
如图1所示,本发明实施例提供了一种页岩气/煤层气井产能评价方法,包括:
步骤100:获取目标井的储层参数和钻完井参数,并根据所述储层参数和所述钻完井参数,建立目标水平井仿真几何模型;所述目标为页岩气或煤层气,所述储层参数包括岩石密度、地层温度、储层初始气体压力、初始气体同位素值和含气孔隙度;所述钻完井参数包括水平井长度、压裂段数、总簇数、压裂缝间距、半缝长、压裂缝高度和压裂缝缝宽;所述目标水平井仿真几何模型包括裂缝区和基质区;所述基质区包括压裂改造区以及未压裂改造区。
在本发明实施例中,选取特定页岩气/煤层气井,即目标井,开展生产过程中产气量和井口气甲烷碳同位素的连续监测,汇总页岩气/煤层气井的储层参数和钻完井参数;根据上述页岩气/煤层气井的储层参数和钻完井参数,建立同时包含(人工/天然)裂缝区、压裂改造区、未压裂改造区三个区域的页岩气/煤层气水平井仿真几何模型,即目标水平井仿真几何模型。
步骤200:根据所述目标水平井仿真几何模型,构建控制方程集合;所述控制方程集合包括第一控制方程组、第二控制方程组和第三控制方程组;所述第一控制方程组包括目标井生产过程中裂缝区内12CH4气体流动控制方程和13CH4气体流动控制方程;所述第二控制方程组包括能够表征目标井生产过程中压裂改造区基质孔隙内游离态12CH4、游离态13CH4、吸附态12CH4和吸附态13CH4耦合流动的控制方程;所述第三控制方程组包括能够表征目标井生产过程中未压裂改造区基质孔隙内游离态12CH4、游离态13CH4、吸附态12CH4和吸附态13CH4耦合流动的控制方程;所述控制方程集合用于根据目标井的关键参数确定目标井的产能评价参数;所述关键参数包括裂缝导流能力、压裂改造区基质孔隙内12CH4的气体扩散系数、压裂改造区基质孔隙内13CH4的气体扩散系数、未压裂改造区基质孔隙内12CH4的气体扩散系数、未压裂改造区基质孔隙内13CH4的气体扩散系数、基质朗格缪尔体积和基质朗格缪尔压力;所述产能评价参数包括游离气产气量、吸附气产气量、日产气量、吸附气比例、游离气比例以及日产气表观同位素值。
在本发明实施例中,步骤200具体包括:
根据所述目标水平井仿真几何模型,建立页岩气/煤层气井生产过程中裂缝区内12CH4和13CH4气体流动的控制方程组。
其中,所述第一控制方程组为:
根据所述目标水平井仿真几何模型,建立能够有效表征页岩气/煤层气井生产过程中压裂改造区基质孔隙内游离态与吸附态12CH4、13CH4耦合流动的控制方程组。
其中,基于菲克扩散第二定律、多组分竞争吸附理论和朗格缪尔等温吸附方程推导得到的第二控制方程组,具体公式为:
式中,B1,B2,C1和C2分别表示为:
P为12CH4的气体压力;P*为13CH4的气体压力;t为生产时间;为压裂改造区基质孔隙内12CH4的气体扩散系数;/>为压裂改造区基质孔隙内13CH4的气体扩散系数;/>为压裂改造区内的孔隙度;/>为压裂改造区质量平衡比参数,/>,/>为压裂改造区内的基质岩石视密度,VL为压裂改造区内基质的朗格缪尔(英文为Langmiur)体积,R为理想气体常数,T为储层温度,Vm为气体标况摩尔体积常数,Z为气体压缩因子;/>为压裂改造区基质内12CH4的朗格缪尔常数,/>,PL为压裂改造区内的基质朗格缪尔压力;/>,/>为压裂改造区基质内13CH4的朗格缪尔常数,α为吸附相和游离相之间的平衡分馏因子。
根据所述目标水平井仿真几何模型,建立能够有效表征页岩气/煤层气井生产过程未压裂改造区基质孔隙内游离态与吸附态12CH4、13CH4耦合流动的控制方程组。
其中,基于Fick扩散第二定律、多组分竞争吸附理论和朗格缪尔等温吸附方程推导得到的控制方程组为:
式中,B1,B2,C1和C2分别表示为:
P为12CH4的气体压力;P*为13CH4的气体压力;t为生产时间;为未压裂改造区基质孔隙内12CH4的气体扩散系数;/>为未压裂改造区基质孔隙内13CH4的气体扩散系数;/>为未压裂改造区的孔隙度;/>为未压裂改造区质量平衡比参数,/>,/>为未压裂改造区内的基质岩石视密度,VL为未压裂改造区内的基质朗格缪尔体积,R为理想气体常数,T为储层温度,Vm为气体标况摩尔体积常数,Z为气体压缩因子;/>为未压裂改造区基质内12CH4的朗格缪尔常数,/>,PL为未压裂改造区内的基质朗格缪尔压力;/>,/>为未压裂改造区基质内13CH4的朗格缪尔常数,α为吸附相和游离相之间的平衡分馏因子。
因为未压裂改造区和压裂改造区只会导致两个区域的孔隙结构参数(如孔隙度、岩石视密度)不一致,而不会改变岩石的吸附能力(如基质朗格缪尔体积、基质内12CH4的朗格缪尔常数、基质朗格缪尔压力和基质内13CH4的朗格缪尔常数等)。
步骤300:确定所述控制方程集合的初始条件和边界条件;所述控制方程集合的初始条件是根据储层初始气体压力和初始气体同位素值确定的,所述控制方程集合的边界条件是根据实测井底流压和日产气同位素值确定的;所述控制方程集合的初始条件包括第一初始条件、第二初始条件和第三初始条件;所述第一初始条件为裂缝区内12CH4气体流动和13CH4气体流动的初始条件;所述第二初始条件为压裂改造区内12CH4气体流动和13CH4气体流动的初始条件;所述第三初始条件为未压裂改造区内12CH4气体流动和13CH4气体流动的初始条件;所述控制方程集合的边界条件包括控制边界条件、连续性边界条件和绝缘边界条件。
其中,所述控制方程集合的初始条件的具体公式为:
式中,t为生产时间;P0为储层初始气体压力,x, y, z分别所述目标水平井仿真几何模型所在空间坐标体系中一位置点的x, y, z轴坐标。
所述控制边界条件的具体公式为:
所述连续性边界条件的具体公式为:
式中,为基质区域与裂缝区域交界面且位于基质区域一侧的12CH4的气体压力梯度,/>为基质区域与裂缝区域交界面且位于裂缝区域一侧的12CH4的气体压力梯度,/>为基质区域与裂缝区域交界面且位于基质区域一侧的13CH4的气体压力梯度,/>为基质区域与裂缝区域交界面且位于裂缝区域一侧的13CH4的气体压力梯度,/>为压裂改造区与未压裂改造区交界面且位于未压裂改造区一侧的12CH4的气体压力梯度,/>为压裂改造区与未压裂改造区交界面且位于位于压裂改造区一侧的12CH4的气体压力梯度,/>为压裂改造区与未压裂改造区交界面且位于未压裂改造区一侧的13CH4的气体压力梯度,/>为压裂改造区与未压裂改造区交界面且位于位于压裂改造区一侧的13CH4的气体压力梯度。两项相等表示通过从一侧流出交界面的气体通量与从交界面流出到另一侧的气体通量一致。/>
所述绝缘边界条件的具体公式为:
式中,为目标水平井仿真几何模型的外边界面处的12CH4的气体压力梯度,/>为目标水平井仿真几何模型的对称面处的12CH4的气体压力梯度,为目标水平井仿真几何模型的外边界面处的13CH4的气体压力梯度,为目标水平井仿真几何模型的对称面处的13CH4的气体压力梯度。
步骤400:在所述控制方程集合的初始条件和边界条件的限制下,求解所述控制方程集合,计算得到目标井在生产过程中任意时刻的产能评价参数;所述产能评价参数包括游离气产气量、吸附气产气量、日产气量、吸附气比例、游离气比例以及日产气表观同位素值。
在本发明实施例中,步骤400具体包括:
在所述控制方程集合的初始条件和边界条件的限制下,求解所述控制方程集合,得到第一分布、第二分布以及第三分布;所述第一分布为裂缝区内12CH4、13CH4压力场的时间-空间分布;所述第二分布为压裂改造区基质孔隙内12CH4、13CH4压力场的时间-空间分布;所述第三分布为未压裂改造区基质孔隙内12CH4、13CH4压力场的时间-空间分布。
根据所述第一分布、所述第二分布以及所述第三分布,对储层内的气体压力和气体覆盖率进行体积积分,计算得到目标井在生产过程中任意时刻的气体产气量;所述气体产气量包括游离态12CH4的产气量、游离态13CH4的产气量、吸附态12CH4的产气量和吸附态13CH4的产气量。
其中,体积积分方法的具体公式为:
式中,、/>、/>和/>分别为ti-1~ti时间段(即为在生产过程中)内游离态12CH4的产气量、游离态13CH4的产气量、吸附态12CH4的产气量、吸附态13CH4的产气量,θ和θ*分别为吸附态12CH4的气体覆盖率、吸附态13CH4的气体覆盖率,其计算公式为:/>
根据计算得到的目标井在生产过程中任意时刻的气体产气量,计算得到目标井在生产过程中任意时刻的产能评价参数。
所述游离气产气量、所述吸附气产气量、所述日产气量的计算公式为:
所述吸附气比例、所述游离气比例的计算公式为:
所述日产气表观同位素值的计算公式为:
步骤500:获取目标井在生产过程中任意时刻实测的日产气量和日产气表观同位素值,并将实测的日产气量和日产气表观同位素值,与计算得到的日产气量和日产气表观同位素值进行拟合,得到拟合结果。
步骤600:将最优拟合结果对应的关键参数确定为目标井在生产过程中的最优关键参数,将所述最优拟合结果对应的产能评价参数确定为目标井在生产过程中最终的产能评价参数;所述最优拟合结果为拟合误差最小的拟合结果。
进一步地,本发明实施例提供的方法还包括:
根据最优关键参数,预测未来时段目标井的产能评价参数。
一个示例为:制定页岩气/煤层气井未来的生产方案,设定相应的井底流压(即公式(7)中的控制边界条件),延长正演计算的时间(时间足够长以保证累积产气量不再明显增加),即可预测得到特定页岩气/煤层气井在指定生产方式下的最终可采收量(EUR),评价页岩气/煤层气井生产过程中吸附气/游离气的动用规律和最终动用比例。
实施例二
为了执行上述实施例一对应的方法,以实现相应的功能和技术效果,下面提供一种页岩气/煤层气井产能评价***。
如图2所示,本发明实施例提供的一种页岩气/煤层气井产能评价***,包括:
目标水平井仿真几何模型建立模块1,用于获取目标井的储层参数和钻完井参数,并根据所述储层参数和所述钻完井参数,建立目标水平井仿真几何模型;所述目标为页岩气或煤层气,所述储层参数包括岩石密度、地层温度、储层初始气体压力、初始气体同位素值和含气孔隙度;所述钻完井参数包括水平井长度、压裂段数、总簇数、压裂缝间距、半缝长、压裂缝高度和压裂缝缝宽;所述目标水平井仿真几何模型包括裂缝区和基质区;所述基质区包括压裂改造区以及未压裂改造区。
控制方程集合构建模块2,用于根据所述目标水平井仿真几何模型,构建控制方程集合;所述控制方程集合包括多个控制方程组,分别为第一控制方程组、第二控制方程组和第三控制方程组;所述第一控制方程组包括目标井生产过程中裂缝区内12CH4气体流动控制方程和13CH4气体流动控制方程;所述第二控制方程组包括能够表征目标井生产过程中压裂改造区基质孔隙内游离态12CH4、游离态13CH4、吸附态12CH4和吸附态13CH4耦合流动的控制方程;所述第三控制方程组包括能够表征目标井生产过程中未压裂改造区基质孔隙内游离态12CH4、游离态13CH4、吸附态12CH4和吸附态13CH4耦合流动的控制方程;所述控制方程集合用于根据目标井的关键参数确定目标井的产能评价参数;所述关键参数包括裂缝导流能力、压裂改造区基质孔隙内12CH4的气体扩散系数、压裂改造区基质孔隙内13CH4的气体扩散系数、未压裂改造区基质孔隙内12CH4的气体扩散系数、未压裂改造区基质孔隙内13CH4的气体扩散系数、基质朗格缪尔体积和基质朗格缪尔压力;所述产能评价参数包括游离气产气量、吸附气产气量、日产气量、吸附气比例、游离气比例以及日产气表观同位素值。
初始条件和边界条件确定模块3,用于确定所述控制方程集合的初始条件和边界条件;所述控制方程集合的初始条件是根据储层初始气体压力和初始气体同位素值确定的,所述控制方程集合的边界条件是根据实测井底流压和日产气同位素值确定的;所述控制方程集合的初始条件包括第一初始条件、第二初始条件和第三初始条件;所述第一初始条件为裂缝区内12CH4气体流动和13CH4气体流动的初始条件;所述第二初始条件为压裂改造区内12CH4气体流动和13CH4气体流动的初始条件;所述第三初始条件为未压裂改造区内12CH4气体流动和13CH4气体流动的初始条件;所述控制方程集合的边界条件包括控制边界条件、连续性边界条件和绝缘边界条件。
任意时刻产能评价参数计算模块4,用于在所述控制方程集合的初始条件和边界条件的限制下,求解所述控制方程集合,计算得到目标井在生产过程中任意时刻的产能评价参数。
拟合结果确定模块5,用于获取目标井在生产过程中任意时刻实测的日产气量和日产气表观同位素值,并将实测的日产气量和日产气表观同位素值,与计算得到的日产气量和日产气表观同位素值进行拟合,得到拟合结果。
最优关键参数和最终产能评价参数确定模块6,用于将最优拟合结果对应的关键参数确定为目标井在生产过程中的最优关键参数,将所述最优拟合结果对应的产能评价参数确定为目标井在生产过程中最终的产能评价参数;所述最优拟合结果为拟合误差最小的拟合结果。
实施例三
本发明实施例提供一种电子设备包括存储器及处理器,该存储器用于存储计算机程序,该处理器运行计算机程序以使电子设备执行实施例一的页岩气/煤层气井产能评价方法。
可选地,上述电子设备可以是服务器。
另外,本发明实施例还提供一种计算机可读存储介质,其存储有计算机程序,该计算机程序被处理器执行时实现实施例一的页岩气/煤层气井产能评价方法。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的***而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (7)
1.一种页岩气/煤层气井产能评价方法,其特征在于,包括:
获取目标井的储层参数和钻完井参数,并根据所述储层参数和所述钻完井参数,建立目标水平井仿真几何模型;所述目标为页岩气或煤层气,所述储层参数包括岩石密度、地层温度、储层初始气体压力、初始气体同位素值和含气孔隙度;所述钻完井参数包括水平井长度、压裂段数、总簇数、压裂缝间距、半缝长、压裂缝高度和压裂缝缝宽;所述目标水平井仿真几何模型包括裂缝区和基质区;所述基质区包括压裂改造区以及未压裂改造区;
根据所述目标水平井仿真几何模型,构建控制方程集合;所述控制方程集合包括多个控制方程组,分别为第一控制方程组、第二控制方程组和第三控制方程组;所述第一控制方程组包括目标井生产过程中裂缝区内12CH4气体流动控制方程和13CH4气体流动控制方程;所述第二控制方程组包括能够表征目标井生产过程中压裂改造区基质孔隙内游离态12CH4、游离态13CH4、吸附态12CH4和吸附态13CH4耦合流动的控制方程;所述第三控制方程组包括能够表征目标井生产过程中未压裂改造区基质孔隙内游离态12CH4、游离态13CH4、吸附态12CH4和吸附态13CH4耦合流动的控制方程;所述控制方程集合用于根据目标井的关键参数确定目标井的产能评价参数;所述关键参数包括裂缝导流能力、压裂改造区基质孔隙内12CH4的气体扩散系数、压裂改造区基质孔隙内13CH4的气体扩散系数、未压裂改造区基质孔隙内12CH4的气体扩散系数、未压裂改造区基质孔隙内13CH4的气体扩散系数、基质朗格缪尔体积和基质朗格缪尔压力;所述产能评价参数包括游离气产气量、吸附气产气量、日产气量、吸附气比例、游离气比例以及日产气表观同位素值;
确定所述控制方程集合的初始条件和边界条件;所述控制方程集合的初始条件是根据储层初始气体压力和初始气体同位素值确定的,所述控制方程集合的边界条件是根据实测井底流压和日产气同位素值确定的;所述控制方程集合的初始条件包括第一初始条件、第二初始条件和第三初始条件;所述第一初始条件为裂缝区内12CH4气体流动和13CH4气体流动的初始条件;所述第二初始条件为压裂改造区内12CH4气体流动和13CH4气体流动的初始条件;所述第三初始条件为未压裂改造区内12CH4气体流动和13CH4气体流动的初始条件;所述控制方程集合的边界条件包括控制边界条件、连续性边界条件和绝缘边界条件;
在所述控制方程集合的初始条件和边界条件的限制下,求解所述控制方程集合,计算得到目标井在生产过程中任意时刻的产能评价参数;
获取目标井在生产过程中任意时刻实测的日产气量和日产气表观同位素值,并将实测的日产气量和日产气表观同位素值,与计算得到的日产气量和日产气表观同位素值进行拟合,得到拟合结果;
将最优拟合结果对应的关键参数确定为目标井在生产过程中的最优关键参数,将所述最优拟合结果对应的产能评价参数确定为目标井在生产过程中最终的产能评价参数;所述最优拟合结果为拟合误差最小的拟合结果;
所述第一控制方程组为:
所述第二控制方程组为:
其中,B1,B2,C1和C2分别表示为:
P为12CH4的气体压力;P*为13CH4的气体压力;t为生产时间;为压裂改造区基质孔隙内12CH4的气体扩散系数;/>为压裂改造区基质孔隙内13CH4的气体扩散系数;/>为压裂改造区内的孔隙度;/>为压裂改造区质量平衡比参数,/>,/>为压裂改造区内的基质岩石视密度,VL为压裂改造区内的基质朗格缪尔体积,R为理想气体常数,T为储层温度,Vm为气体标况摩尔体积常数,Z为气体压缩因子;/>为压裂改造区基质内12CH4的朗格缪尔常数,/>,PL为压裂改造区内的基质朗格缪尔压力;/>,/>为压裂改造区基质内13CH4的朗格缪尔常数,α为吸附相和游离相之间的平衡分馏因子;
所述第三控制方程组为:
其中,B1,B2,C1和C2分别表示为:
3.根据权利要求1所述的一种页岩气/煤层气井产能评价方法,其特征在于,所述控制边界条件为:
所述连续性边界条件为:
式中,为基质区域与裂缝区域交界面且位于基质区域一侧的12CH4的气体压力梯度,/>为基质区域与裂缝区域交界面且位于裂缝区域一侧的12CH4的气体压力梯度,/>为基质区域与裂缝区域交界面且位于基质区域一侧的13CH4的气体压力梯度,/>为基质区域与裂缝区域交界面且位于裂缝区域一侧的13CH4的气体压力梯度,/>为压裂改造区与未压裂改造区交界面且位于未压裂改造区一侧的12CH4的气体压力梯度,/>为压裂改造区与未压裂改造区交界面且位于位于压裂改造区一侧的12CH4的气体压力梯度,/>为压裂改造区与未压裂改造区交界面且位于未压裂改造区一侧的13CH4的气体压力梯度,/>为压裂改造区与未压裂改造区交界面且位于位于压裂改造区一侧的13CH4的气体压力梯度;
所述绝缘边界条件为:
4.根据权利要求1所述的一种页岩气/煤层气井产能评价方法,其特征在于,所述在所述控制方程集合的初始条件和边界条件的限制下,求解所述控制方程集合,计算得到目标井在生产过程中任意时刻的产能评价参数,具体包括:
在所述控制方程集合的初始条件和边界条件的限制下,求解所述控制方程集合,得到第一分布、第二分布以及第三分布;所述第一分布为裂缝区内12CH4、13CH4压力场的时间-空间分布;所述第二分布为压裂改造区基质孔隙内12CH4、13CH4压力场的时间-空间分布;所述第三分布为未压裂改造区基质孔隙内12CH4、13CH4压力场的时间-空间分布;
根据所述第一分布、所述第二分布以及所述第三分布,对储层内的气体压力和气体覆盖率进行体积积分,计算得到目标井在生产过程中任意时刻的气体产气量;所述气体产气量包括游离态12CH4的产气量、游离态13CH4的产气量、吸附态12CH4的产气量和吸附态13CH4的产气量;
根据计算得到的目标井在生产过程中任意时刻的气体产气量,计算得到目标井在生产过程中任意时刻的产能评价参数。
5.根据权利要求1所述的一种页岩气/煤层气井产能评价方法,其特征在于,还包括:
根据最优关键参数,预测未来时段目标井的产能评价参数。
6.一种页岩气/煤层气井产能评价***,其特征在于,包括:
目标水平井仿真几何模型建立模块,用于获取目标井的储层参数和钻完井参数,并根据所述储层参数和所述钻完井参数,建立目标水平井仿真几何模型;所述目标为页岩气或煤层气,所述储层参数包括岩石密度、地层温度、储层初始气体压力、初始气体同位素值和含气孔隙度;所述钻完井参数包括水平井长度、压裂段数、总簇数、压裂缝间距、半缝长、压裂缝高度和压裂缝缝宽;所述目标水平井仿真几何模型包括裂缝区和基质区;所述基质区包括压裂改造区以及未压裂改造区;
控制方程集合构建模块,用于根据所述目标水平井仿真几何模型,构建控制方程集合;所述控制方程集合包括多个控制方程组,分别为第一控制方程组、第二控制方程组和第三控制方程组;所述第一控制方程组包括目标井生产过程中裂缝区内12CH4气体流动控制方程和13CH4气体流动控制方程;所述第二控制方程组包括能够表征目标井生产过程中压裂改造区基质孔隙内游离态12CH4、游离态13CH4、吸附态12CH4和吸附态13CH4耦合流动的控制方程;所述第三控制方程组包括能够表征目标井生产过程中未压裂改造区基质孔隙内游离态12CH4、游离态13CH4、吸附态12CH4和吸附态13CH4耦合流动的控制方程;所述控制方程集合用于根据目标井的关键参数确定目标井的产能评价参数;所述关键参数包括裂缝导流能力、压裂改造区基质孔隙内12CH4的气体扩散系数、压裂改造区基质孔隙内13CH4的气体扩散系数、未压裂改造区基质孔隙内12CH4的气体扩散系数、未压裂改造区基质孔隙内13CH4的气体扩散系数、基质朗格缪尔体积和基质朗格缪尔压力;所述产能评价参数包括游离气产气量、吸附气产气量、日产气量、吸附气比例、游离气比例以及日产气表观同位素值;
初始条件和边界条件确定模块,用于确定所述控制方程集合的初始条件和边界条件;所述控制方程集合的初始条件是根据储层初始气体压力和初始气体同位素值确定的,所述控制方程集合的边界条件是根据实测井底流压和日产气同位素值确定的;所述控制方程集合的初始条件包括第一初始条件、第二初始条件和第三初始条件;所述第一初始条件为裂缝区内12CH4气体流动和13CH4气体流动的初始条件;所述第二初始条件为压裂改造区内12CH4气体流动和13CH4气体流动的初始条件;所述第三初始条件为未压裂改造区内12CH4气体流动和13CH4气体流动的初始条件;所述控制方程集合的边界条件包括控制边界条件、连续性边界条件和绝缘边界条件;
任意时刻产能评价参数计算模块,用于在所述控制方程集合的初始条件和边界条件的限制下,求解所述控制方程集合,计算得到目标井在生产过程中任意时刻的产能评价参数;
拟合结果确定模块,用于获取目标井在生产过程中任意时刻实测的日产气量和日产气表观同位素值,并将实测的日产气量和日产气表观同位素值,与计算得到的日产气量和日产气表观同位素值进行拟合,得到拟合结果;
最优关键参数和最终产能评价参数确定模块,用于将最优拟合结果对应的关键参数确定为目标井在生产过程中的最优关键参数,将所述最优拟合结果对应的产能评价参数确定为目标井在生产过程中最终的产能评价参数;所述最优拟合结果为拟合误差最小的拟合结果;
所述第一控制方程组为:
所述第二控制方程组为:
其中,B1,B2,C1和C2分别表示为:
P为12CH4的气体压力;P*为13CH4的气体压力;t为生产时间;为压裂改造区基质孔隙内12CH4的气体扩散系数;/>为压裂改造区基质孔隙内13CH4的气体扩散系数;/>为压裂改造区内的孔隙度;/>为压裂改造区质量平衡比参数,/>,/>为压裂改造区内的基质岩石视密度,VL为压裂改造区内的基质朗格缪尔体积,R为理想气体常数,T为储层温度,Vm为气体标况摩尔体积常数,Z为气体压缩因子;/>为压裂改造区基质内12CH4的朗格缪尔常数,/>,PL为压裂改造区内的基质朗格缪尔压力;/>,/>为压裂改造区基质内13CH4的朗格缪尔常数,α为吸附相和游离相之间的平衡分馏因子;
所述第三控制方程组为:
其中,B1,B2,C1和C2分别表示为:
7.一种电子设备,其特征在于,包括存储器及处理器,所述存储器用于存储计算机程序,所述处理器运行所述计算机程序以使所述电子设备执行根据权利要求1至5中任一项所述的一种页岩气/煤层气井产能评价方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310119422.XA CN115860266B (zh) | 2023-02-16 | 2023-02-16 | 一种页岩气/煤层气井产能评价方法、***及电子设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310119422.XA CN115860266B (zh) | 2023-02-16 | 2023-02-16 | 一种页岩气/煤层气井产能评价方法、***及电子设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115860266A true CN115860266A (zh) | 2023-03-28 |
CN115860266B CN115860266B (zh) | 2023-05-09 |
Family
ID=85658144
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310119422.XA Active CN115860266B (zh) | 2023-02-16 | 2023-02-16 | 一种页岩气/煤层气井产能评价方法、***及电子设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115860266B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117147318A (zh) * | 2023-10-12 | 2023-12-01 | 东北石油大学三亚海洋油气研究院 | 一种页岩储层多尺度孔隙应力敏感性评价方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110054795A1 (en) * | 2009-08-27 | 2011-03-03 | Conocophillips Company | Petrophysical Evaluation of Subterranean Formations |
US20120312530A1 (en) * | 2001-04-11 | 2012-12-13 | Pope John M | In-Situ Detection and Analysis of Methane in Coal Bed Methane Formations with Spectrometers |
CN104389594A (zh) * | 2014-10-13 | 2015-03-04 | 成都创源油气技术开发有限公司 | 页岩气井产能评价预测方法 |
CN105046006A (zh) * | 2015-07-29 | 2015-11-11 | 中国石油天然气股份有限公司 | 一种页岩气藏水平井多段压裂产能预测方法及装置 |
WO2016041189A1 (zh) * | 2014-09-19 | 2016-03-24 | 杨顺伟 | 一种评价页岩气储层及寻找甜点区的方法 |
CN107133373A (zh) * | 2016-02-29 | 2017-09-05 | 中国石油化工股份有限公司 | 一种页岩气藏、井筒及地面管网的耦合模拟方法 |
US20180016896A1 (en) * | 2016-07-15 | 2018-01-18 | Geoisochem Corporation | Assessing Permeability |
CN108665086A (zh) * | 2017-04-01 | 2018-10-16 | 中国石油化工股份有限公司 | 一种页岩气井产量预测计算方法 |
CN109594968A (zh) * | 2017-09-28 | 2019-04-09 | 中国石油化工股份有限公司 | 一种页岩气多段压裂水平井压后裂缝参数评价方法及*** |
US20220049604A1 (en) * | 2020-08-12 | 2022-02-17 | Southwest Petroleum University | Shale gas well dynamic production allocating method |
-
2023
- 2023-02-16 CN CN202310119422.XA patent/CN115860266B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120312530A1 (en) * | 2001-04-11 | 2012-12-13 | Pope John M | In-Situ Detection and Analysis of Methane in Coal Bed Methane Formations with Spectrometers |
US20110054795A1 (en) * | 2009-08-27 | 2011-03-03 | Conocophillips Company | Petrophysical Evaluation of Subterranean Formations |
WO2016041189A1 (zh) * | 2014-09-19 | 2016-03-24 | 杨顺伟 | 一种评价页岩气储层及寻找甜点区的方法 |
CN104389594A (zh) * | 2014-10-13 | 2015-03-04 | 成都创源油气技术开发有限公司 | 页岩气井产能评价预测方法 |
CN105046006A (zh) * | 2015-07-29 | 2015-11-11 | 中国石油天然气股份有限公司 | 一种页岩气藏水平井多段压裂产能预测方法及装置 |
CN107133373A (zh) * | 2016-02-29 | 2017-09-05 | 中国石油化工股份有限公司 | 一种页岩气藏、井筒及地面管网的耦合模拟方法 |
US20180016896A1 (en) * | 2016-07-15 | 2018-01-18 | Geoisochem Corporation | Assessing Permeability |
CN108665086A (zh) * | 2017-04-01 | 2018-10-16 | 中国石油化工股份有限公司 | 一种页岩气井产量预测计算方法 |
CN109594968A (zh) * | 2017-09-28 | 2019-04-09 | 中国石油化工股份有限公司 | 一种页岩气多段压裂水平井压后裂缝参数评价方法及*** |
US20220049604A1 (en) * | 2020-08-12 | 2022-02-17 | Southwest Petroleum University | Shale gas well dynamic production allocating method |
Non-Patent Citations (1)
Title |
---|
李文镖 等: "页岩气/煤层气运移过程中的同位素分馏研究进展" * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117147318A (zh) * | 2023-10-12 | 2023-12-01 | 东北石油大学三亚海洋油气研究院 | 一种页岩储层多尺度孔隙应力敏感性评价方法 |
CN117147318B (zh) * | 2023-10-12 | 2024-05-07 | 东北石油大学三亚海洋油气研究院 | 一种页岩储层多尺度孔隙应力敏感性评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115860266B (zh) | 2023-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109441422B (zh) | 一种页岩气井间距优化开采方法 | |
CN110472276B (zh) | 一种裂缝型油气储层加密井改造优化方法 | |
CN106351651B (zh) | 页岩气井产能的预测方法及装置 | |
CN104750896B (zh) | 一种缝洞型碳酸盐岩油藏数值模拟方法 | |
CN110210157A (zh) | 一种页岩气藏压裂水平井非稳态扩散下产能计算方法 | |
CN104948163A (zh) | 一种页岩气井产能测定方法 | |
CN108661631B (zh) | 一种产量预测方法 | |
CN105089615B (zh) | 一种基于油藏模型的测井数据历史回归处理方法 | |
US20170212276A1 (en) | Method for characterizing and exploiting a subterranean formation comprising a network of fractures | |
CN112282744B (zh) | 一种非常规油气藏井网部署优化方法和装置 | |
CN110162808B (zh) | 一种确定页岩气井吸附气和游离气产出贡献的方法 | |
CN111927417A (zh) | 一种页岩气分段压裂水平井组储量动用状况评价方法 | |
CN115860266A (zh) | 一种页岩气/煤层气井产能评价方法、***及电子设备 | |
CN115330060A (zh) | 基于储层与工程因素分析的页岩油水平井产能计算方法 | |
Huang et al. | A review of flow mechanism and inversion methods of fracture network in shale gas reservoirs | |
CN111950112B (zh) | 一种适用于底部封闭的碳酸盐岩储层动态分析方法 | |
CN107451671B (zh) | 用于预测页岩地层压裂后初始产能的方法及*** | |
CN111950111B (zh) | 一种适用于底部开放的碳酸盐岩储层动态分析方法 | |
CN111997581B (zh) | 非均质油藏开发方法、装置及电子设备 | |
CN112381259A (zh) | 一种致密砾岩储层的产能动态预测计算方法 | |
CN112035993B (zh) | 一种底部定压的碳酸盐岩储层测试评价方法 | |
CN107704646B (zh) | 一种致密储层体积改造后的建模方法 | |
CN108664678B (zh) | 一种产量预测方法 | |
Ballin et al. | New reservoir dynamic connectivity measurement for efficient well placement strategy analysis under depletion | |
CN110486005B (zh) | 一种用于识别页岩气井甜点层的方法 |
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 |