CN112065339B - 多储层气藏合采产能预测方法及装置 - Google Patents
多储层气藏合采产能预测方法及装置 Download PDFInfo
- Publication number
- CN112065339B CN112065339B CN202010907844.XA CN202010907844A CN112065339B CN 112065339 B CN112065339 B CN 112065339B CN 202010907844 A CN202010907844 A CN 202010907844A CN 112065339 B CN112065339 B CN 112065339B
- Authority
- CN
- China
- Prior art keywords
- gas
- reservoir
- representing
- flow
- phase
- 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.)
- Active
Links
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 123
- 238000000034 method Methods 0.000 title claims abstract description 45
- 239000007788 liquid Substances 0.000 claims abstract description 56
- 239000011435 rock Substances 0.000 claims abstract description 54
- 230000005465 channeling Effects 0.000 claims abstract description 34
- 230000005514 two-phase flow Effects 0.000 claims abstract description 30
- 239000003245 coal Substances 0.000 claims description 67
- 239000010410 layer Substances 0.000 claims description 37
- 239000012071 phase Substances 0.000 claims description 35
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 31
- 239000007791 liquid phase Substances 0.000 claims description 24
- 230000035699 permeability Effects 0.000 claims description 18
- 239000011148 porous material Substances 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 12
- 239000011229 interlayer Substances 0.000 claims description 12
- 230000015572 biosynthetic process Effects 0.000 claims description 11
- 239000004576 sand Substances 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 6
- 239000008398 formation water Substances 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 239000011800 void material Substances 0.000 claims description 5
- 230000000630 rising effect Effects 0.000 claims description 4
- 230000004927 fusion Effects 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 2
- 230000007704 transition Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 12
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 12
- 230000008859 change Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000001179 sorption measurement Methods 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 230000002706 hydrostatic effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000003795 desorption Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/14—Obtaining from a multiple-zone well
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Primary Health Care (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Human Resources & Organizations (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Animal Husbandry (AREA)
- Agronomy & Crop Science (AREA)
- Filling Or Discharging Of Gas Storage Vessels (AREA)
Abstract
本发明公开了一种多储层气藏合采产能预测方法及装置,该方法包括:获取待开采气藏各储层的岩石属性信息,其中,待开采气藏包括:相邻的第一储层和第二储层,第一储层和第二储层的产气方式不同;根据第一储层的岩石属性信息,确定第一储层的第一初始含气量;根据第二储层的岩石属性信息,确定第二储层的第二初始含气量;根据第一储层和第二储层的岩石属性信息,确定第一储层与第二储层之间的气体窜流量;确定合采井筒内气液两相流压降;根据第一初始含气量、第二初始含气量、气体窜流量和合采井筒内气液两相流压降,确定待开采气藏的产气量。本发明实现了对多储层气藏合采产能的准确定量预测。
Description
技术领域
本发明涉及气藏开发领域,尤其涉及一种多储层气藏合采产能预测方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
油气藏数值模拟是油气田科学开发的基础,是油田生产指标计算、剩余油分布预测的必备技术。目前,市面上的一些油气藏数值模拟软件(例如,Eclips、CMG、Commet等)在模拟单一油气藏时比较成熟,但在模拟多储层气藏时有限制。现有技术大多基于定性分析以及现场经验,来对多储层气藏开发过程中的产气规律进行分析,唯一可借鉴的基础理论模型是常规多层油气藏合采的研究成果。
与常规油气藏相比,多储层气藏产气机理具有多样性(例如,含有砂岩气的多煤层气藏中,煤层有吸附气,也有游离气;而砂岩夹层只有游离气)。在多储层气藏开采过程中,产气方式不同的两个储层会相互影响(例如,煤层与砂岩层之间存在气体窜流),使得常规多层油气藏的研究成果不能直接用来解释多煤层气藏中气/水的运移规律。
因而,如何提供一种多储层气藏合采产能的定量预测方法,对于多储层气藏的开发具有十分重要的意义。
发明内容
本发明实施例提供一种多储层气藏合采产能预测方法,用以解决现有常规多储层油气藏合采产能预测方法无法对多储层气藏合采产能进行准确预测的技术问题,该方法包括:获取待开采气藏各储层的岩石属性信息,其中,待开采气藏包括:相邻的第一储层和第二储层,第一储层和第二储层的产气方式不同;根据第一储层的岩石属性信息,确定第一储层的第一初始含气量;根据第二储层的岩石属性信息,确定第二储层的第二初始含气量;根据第一储层和第二储层的岩石属性信息,确定第一储层与第二储层之间的气体窜流量;确定合采井筒内气液两相流压降;根据第一初始含气量、第二初始含气量、气体窜流量和合采井筒内气液两相流压降,确定待开采气藏的产气量。
本发明实施例还提供一种多储层气藏合采产能预测装置,用以解决现有常规多储层油气藏合采产能预测方法无法对多储层气藏合采产能进行准确预测的技术问题,该装置包括:岩石属性信息获取模块,用于获取待开采气藏各储层的岩石属性信息,其中,待开采气藏包括:相邻的第一储层和第二储层,第一储层和第二储层的产气方式不同;第一储层含气量确定模块,用于根据第一储层的岩石属性信息,确定第一储层的第一初始含气量;第二储层含气量确定模块,用于根据第二储层的岩石属性信息,确定第二储层的第二初始含气量;气体窜流量确定模块,用于根据第一储层和第二储层的岩石属性信息,确定第一储层与第二储层之间的气体窜流量;井筒气液两相流压降确定模块,用于确定合采井筒内气液两相流压降;合采产气量确定模块,用于根据第一初始含气量、第二初始含气量、气体窜流量和合采井筒内气液两相流压降,确定待开采气藏的产气量。
本发明实施例还提供一种计算机设备,用以解决现有常规多储层油气藏合采产能预测方法无法对多储层气藏合采产能进行准确预测的技术问题,该计算机设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时实现上述多储层气藏合采产能预测方法。
本发明实施例还提供一种计算机可读存储介质,用以解决现有常规多储层油气藏合采产能预测方法无法对多储层气藏合采产能进行准确预测的技术问题,计算机可读存储介质存储有执行上述多储层气藏合采产能预测方法的计算机程序。
本发明实施例中,在对多储层气藏合采产能进行预测的时候,首先根据各储层的岩石属性信息,确定各储层的初始含气量,然后根据相邻两个储层的岩石属性信息,确定相邻两个储层之间的气体窜流量,结合合采井筒在各储层位置处的井底流压,来确定多储层气藏合采产能。由于多储层气藏中不同储层的产气方式可能不同,且产气方式不同的相邻两个储层之间会有气体窜流影响,本发明实施例综合考虑各储层产气量、相邻两个储层之间的气体窜流量,以及合采井筒在各储层位置处的井底流压来确定多储层气藏合采产能,实现了对多储层气藏合采产能的准确定量预测。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例中提供的一种多储层气藏合采产能预测方法流程图;
图2为本发明实施例中提供的一种多煤层合采垂直井筒内气液两相流型示意图;
图3为本发明实施例中提供的一种确定各段截面含气率和密度的方法流程图;
图4为本发明实施例中提供的多储层气藏合采产能预测装置示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
在本说明书的描述中,所使用的“包含”、“包括”、“具有”、“含有”等,均为开放性的用语,即意指包含但不限于。参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”等的描述意指结合该实施例或示例描述的具体特征、结构或者特点包含于本申请的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。各实施例中涉及的步骤顺序用于示意性说明本申请的实施,其中的步骤顺序不作限定,可根据需要作适当调整。
本发明实施例中提供了一种多储层气藏合采产能预测方法,图1为本发明实施例中提供的一种多储层气藏合采产能预测方法流程图,如图1所示,该方法可以包括如下步骤:
S101,获取待开采气藏各储层的岩石属性信息,其中,待开采气藏包括:相邻的第一储层和第二储层,第一储层和第二储层的产气方式不同。
需要说明的是,上述S101中的待开采气藏可以是任意一种包含多个储层的多储层气藏;上述岩石属性信息包括但不限于第一储层的岩石属性信息和第二储层的岩石属性信息;上述第一储层和第二储层可以是多储层气藏中任意两个相邻的储层。
作为一种可选的实施例,本发明实施例中的待开采气藏可以是含有砂岩夹层的多煤层气藏。下面以含有砂岩夹层的多煤层气藏为例来对本发明实施例进行说明,假设第一储层煤层;第二储层为砂岩层。
S102,根据第一储层的岩石属性信息,确定第一储层的第一初始含气量。
由于煤层中的含气量由吸附气和游离气组成,作为一种可选的实施方式,当第一储层为煤层的情况下,上述S102具体可以通过如下公式,计算第一储层的第一初始含气量:
其中,V1表示第一储层的第一初始含气量,m3;Vc表示煤层体积,m3;φc表示煤层孔隙度;Scg表示煤层含气饱和度;Scw表示煤层含水饱和度;Bg表示煤层气体体积分数;Bw表示煤层水体积分数;Rcw表示煤层中液相溶解气液比;VL表示朗格缪尔体积;PL表示朗格缪尔压力;p表示煤层孔隙压力,MPa。
需要注意的是,Langmuir等温吸附公式被广泛用于描述甲烷在煤层中的吸附和解吸,因此使用了Langmuir等温吸附来描述甲烷在煤层中的吸附情况。
S103,根据第二储层的岩石属性信息,确定第二储层的第二初始含气量。
由于砂岩中的含气量只有游离气,因而,作为一种可选的实施方式,当第二储层为砂岩层的情况下,上述S103具体可以通过如下公式,计算第二储层的第二初始含气量:
其中,V2表示第二储层的第二初始含气量,m3;Vs表示砂岩层体积,m3;φs表示砂岩层孔隙度;Ssg表示砂岩层含气饱和度;Ssw表示砂岩层含水饱和度;Bg表示砂岩层气体体积分数;Bw表示砂岩层水体积分数;Rsw表示砂岩层中液相溶解气液比。
S104,根据第一储层和第二储层的岩石属性信息,确定第一储层与第二储层之间的气体窜流量。
当含有砂岩夹层的多煤层气藏开采过程中,煤层与砂岩层之间会发生气体窜流,因而,作为一种可选的实施方式,当第一储层为煤层、第二储层为砂岩层的情况下,上述S104具体可以通过如下公式,计算第一储层与第二储层之间的气体窜流量:
其中,qgcross表示煤层与砂岩层之间的气体窜流量,m3/s;pcg表示煤岩层孔隙气相压力;psg表示砂岩层孔隙气相压力;μcg表示煤层中气体的动力粘度MPa·s;μsg表示砂岩层中气体的动力粘度;ρcg表示煤层条件下气体的密度,kg/m3;ρsg表示砂岩层条件下气体的密度,kg/m3;kcrg表示煤岩层气体相对渗透率;ksrg表示砂岩层气体相对渗透率;hc表示煤层的厚度,m;hs表示砂岩层的厚度,m;kcz表示煤岩层垂向渗透率;ksz表示砂岩层垂向渗透率;β1和β2表示基于熔合界面的层间窜流阻力系数;ρsc表示地面条件下气体的密度,kg/m3;Rs表示单位体积水中溶解气量;pcw表示煤岩层孔隙水相压力;psw表示砂岩层孔隙水相压力;μcw表示煤层中水相动力粘度;kcrw表示煤岩层水相对渗透率;μsw表示砂岩层水相粘度;ksrw表示砂岩层水相对渗透率。
S105,确定合采井筒内气液两相流压降。
作为一种可选的实施方式,上述S105具体可以通过如下公式确定合采井筒内气液两相流压降:
其中,
ρwb=(1-α)ρw (5)
其中,pwf2表示合采井筒中第二储层位置处的第二井底流压;pwf1表示合采井筒中第一储层位置处的第一井底流压;ρwb表示合采井筒中不同位置处气液两相密度;g表示重力加速度,9.8m/s2;z1表示第一储层的垂直深度,m;z2表示第二储层的垂直深度,m;ρw表示水密度;α表示含气截面率。
需要说明的是,含有砂岩夹层的多煤层气藏井产气阶段,垂直井筒内为气液两相(例如,气水两相)流动,最下层煤层至动液面,由于压力变化和上层产气补给,井筒内气液两相流动相关参数甚至流型是发生变化的。因而,为了准确预测各层井底流压,对垂直井筒内的相关流动情况(气液流型和截面含气率)进行准确预测至关重要。根据气液流型图和煤层气产量数据分析可以确定,合采井筒内的气液流型主要是泡状流和段塞流,如图2所示,因此本发明实施例只对这两种流型的截面含气率以及两者的转变准则进行阐述。
首先,合采井筒内气液两相流的气液流型由泡状流向段塞流的转变满足如下条件:
气泡在静液中的上升速度表示如下:
泡状流的截面含气率可以通过如下公式计算:
其中,σ表示气液界面张力,10-3N·m-1;ρl表示液相密度,kg·m-3;ρg表示气相密度,kg·m-3;g表示重力加速度,m·s-2;αg表示泡状流的截面含气率,无量纲;vsg表示气体表观流速,m·s-1;vsl表示液相表观流速,m·s-1。
需要说明的是,充分发展的段塞流平均截面含气率,则段塞流液段截面含气率可以通过如下公式计算:
其中,
αlgvlg+(1-αlg)vll=vm (10)
其中,αlg表示液段截面含气率,无量纲;σ表示气液界面张力;ρl表示液相密度;ρg表示气相密度;g表示重力加速度;Cl表示与液相雷诺数有关的一个系数;D表示井筒内径;μl表示液相粘度,Pa·s;vm表示气液两相的混合流速,m·s-1;vlg表示液段气相流速,m·s-1;vll表示液段液相流速,m·s-1;n表示与雷诺数有关的一个常数,无量纲。
Taylor气泡的速度:
vTB=1.2vll+0.35ZTB (12)
其中,vTB为Taylor气泡的速度,m·s-1;zTB为Taylor气泡在静液中的上升速度,m·s-1。
Taylor气泡段含气率和液膜流速:
αTBgvTB-(1-αTBg)vf=vm (14)
其中,αTBg为Taylor气泡段含气率,无量纲;vf为液膜流速,m·s-1;X为洛-马参数,无量纲;Y为管路倾角对流型的影响参数。
Taylor气泡段长度与段塞单元段长度比值:
其中,LTB表示泰勒气泡段长度,L表示段塞单元长度;αTBg表示Taylor气泡段含气率;alg表示液段截面含气率;vsg表示气体表观流速;vlg表示液段气相流速。
段塞流平均截面含气率:
其中,αTBg为泰勒气泡段含气率,无量纲;vf为液膜流速,m·s-1;X为洛-马参数,无量纲;Y为管路倾角对流型的影响参数;vsl表示液相表观流速。
S106,根据第一初始含气量、第二初始含气量、气体窜流量和合采井筒内气液两相流压降,确定待开采气藏的产气量。
需要说明的是,上述S102、S103、S104、S105的执行顺序仅为示例,原则上可以任意交换,或并列分别执行。
由上可知,在本发明实施例提供的多储层气藏合采井产气量预测方法,在获取到多储层气藏各储层的岩石属性信息后,根据各储层的岩石属性信息,确定各储层的初始含气量,然后根据相邻两个储层的岩石属性信息,确定相邻两个储层之间的气体窜流量,结合合采井筒在各储层位置处的井底流压,来确定多储层气藏合采产能。由于多储层气藏中不同储层的产气方式可能不同,且产气方式不同的相邻两个储层之间会有气体窜流影响,本发明实施例综合考虑各储层产气量、相邻两个储层之间的气体窜流量,以及合采井筒在各储层位置处的井底流压来确定多储层气藏合采产能,实现了对多储层气藏合采产能的准确定量预测。
图3为本发明实施例中提供的一种确定各段截面含气率和密度的方法流程图,如图3所示,假设某一含有砂岩夹层的多煤层气藏包括煤层1、砂岩层和煤层2。通过图3所示的流程,计算各段截面含气率和密度,直到根据第n+1次的管流密度与第n次的管流密度之差的绝对值与第n次的管流密度的比值小于一个阈值(例如,0.01),结束流程。
本发明实施例中还提供了一种多储层气藏合采井产气量预测装置,如下面的实施例所述。由于该装置实施例解决问题的原理与多储层气藏合采井产气量预测方法相似,因此该装置实施例的实施可以参见方法的实施,重复之处不再赘述。
图4为本发明实施例中提供的多储层气藏合采产能预测装置示意图,如图4所示,该装置包括:岩石属性信息获取模块41、第一储层含气量确定模块42、第二储层含气量确定模块43、气体窜流量确定模块44、井筒气液两相流压降确定模块45、合采产气量确定模块46。
其中,岩石属性信息获取模块41,用于获取待开采气藏各储层的岩石属性信息,其中,待开采气藏包括:相邻的第一储层和第二储层,第一储层和第二储层的产气方式不同;第一储层含气量确定模块42,用于根据第一储层的岩石属性信息,确定第一储层的第一初始含气量;第二储层含气量确定模块43,用于根据第二储层的岩石属性信息,确定第二储层的第二初始含气量;气体窜流量确定模块44,用于根据第一储层和第二储层的岩石属性信息,确定第一储层与第二储层之间的气体窜流量;井筒气液两相流压降确定模块45,用于确定合采井筒内气液两相流压降;合采产气量确定模块46,用于根据第一初始含气量、第二初始含气量、气体窜流量和合采井筒内气液两相流压降,确定待开采气藏的产气量。
由上可知,在本发明实施例提供的多储层气藏合采井产气量预测方法,通过岩石属性信息获取模块41获取多储层气藏各储层的岩石属性信息,通过第一储层含气量确定模块42根据第一储层的岩石属性信息,确定第一储层的初始含气量;通过第二储层含气量确定模块43根据第二储层的岩石属性信息,确定第二储层的初始含气量;通过气体窜流量确定模块44根据相邻两个储层的岩石属性信息,确定相邻两个储层之间的气体窜流量,通过井筒气液两相流压降确定模块45确定合采井筒内气液两相流压降;通过合采产气量确定模块46根据第一储层和第二储层的初始含气量,以及第一储层和第二储层之间的气体窜流量,结合合采井筒内气液两相流压降,来确定多储层气藏合采产能。由于多储层气藏中不同储层的产气方式可能不同,且产气方式不同的相邻两个储层之间会有气体窜流影响,本发明实施例综合考虑各储层产气量、相邻两个储层之间的气体窜流量,以及合采井筒在各储层位置处的井底流压来确定多储层气藏合采产能,实现了对多储层气藏合采产能的准确定量预测。
本发明实施例还提供一种计算机设备,用以解决现有常规多储层油气藏合采产能预测方法无法对多储层气藏合采产能进行准确预测的技术问题,该计算机设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时实现上述多储层气藏合采产能预测方法。
本发明实施例还提供一种计算机可读存储介质,用以解决现有常规多储层油气藏合采产能预测方法无法对多储层气藏合采产能进行准确预测的技术问题,计算机可读存储介质存储有执行上述多储层气藏合采产能预测方法的计算机程序。
综上所述,本发明实施例基于煤岩双孔-单渗模型和砂岩单孔-单渗模型,以及煤岩与砂岩层间窜流模型和井筒气-水两相管流模型,构建了考虑层间窜流的多煤层气藏全过程气-水两相流动数值模型;采用数值的方法对模型进行求解,获得多煤层气合采井产气量和产水量以及层间窜流量随时间的变化,分析了多煤层气藏合采井产气动态变化规律,井筒气-水两相管流对合采的影响,以及窜流位置对合采井产气量的影响。通过本发明实施例提供的多储层气藏合采产能预测方法,能够实现对多储层气藏合采产能的准确定量预测。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种多储层气藏合采产能预测方法,其特征在于,包括:
获取待开采气藏各储层的岩石属性信息,其中,所述待开采气藏包括:相邻的第一储层和第二储层,第一储层和第二储层的产气方式不同;
根据第一储层的岩石属性信息,确定第一储层的第一初始含气量;
根据第二储层的岩石属性信息,确定第二储层的第二初始含气量;
根据第一储层和第二储层的岩石属性信息,确定第一储层与第二储层之间的气体窜流量;
确定合采井筒内气液两相流压降;
根据所述第一初始含气量、所述第二初始含气量、所述气体窜流量和所述合采井筒内气液两相流压降,确定待开采气藏的产气量;
其中,第一储层为煤层,第二储层为砂岩层,其中,根据第一储层和第二储层的岩石属性信息,确定第一储层与第二储层之间的气体窜流量,包括:
通过如下公式,计算第一储层与第二储层之间的气体窜流量:
其中,qgcross表示煤层与砂岩层之间的气体窜流量;pcg表示煤岩层孔隙气相压力;psg表示砂岩层孔隙气相压力;μcg表示煤层中气体的动力粘度;μsg表示砂岩层中气体的动力粘度;ρcg表示煤层条件下气体的密度;ρsg表示砂岩层条件下气体的密度;kcrg表示煤岩层气体相对渗透率;ksrg表示砂岩层气体相对渗透率;hc表示煤层的厚度;hs表示砂岩层的厚度;kcz表示煤岩层垂向渗透率;ksz表示砂岩层垂向渗透率;β1和β2表示基于熔合界面的层间窜流阻力系数;ρsc表示地面条件下气体的密度;Rs表示单位体积水中溶解气量;pcw表示煤岩层孔隙水相压力;psw表示砂岩层孔隙水相压力;μcw表示煤层中水相动力粘度;kcrw表示煤岩层水相对渗透率;μsw表示砂岩层水相粘度;ksrw表示砂岩层水相对渗透率。
5.如权利要求1所述的方法,其特征在于,所述合采井筒内气液两相流的气液流型包括:泡状流和段塞流,所述方法还包括:
通过如下公式计算泡状流的截面含气率:
其中,σ表示气液界面张力;ρl表示液相密度;ρg表示气相密度;g表示重力加速度;αg表示泡状流的截面含气率;vsg表示气体表观流速;vsl表示液相表观流速;
通过如下公式计算段塞流的截面含气率:
其中,αlgvlg+(1-αlg)vll=vm;
其中,αlg表示液段截面含气率;σ表示气液界面张力;ρl表示液相密度;ρg表示气相密度;g表示重力加速度;Cl表示与液相雷诺数有关的一个系数;D表示井筒内径;μl表示液相粘度;vm表示气液两相的混合流速;vlg表示液段气相流速;vll表示液段液相流速;n表示与雷诺数有关的一个常数。
7.一种多储层气藏合采产能预测装置,其特征在于,包括:
岩石属性信息获取模块,用于获取待开采气藏各储层的岩石属性信息,其中,所述待开采气藏包括:相邻的第一储层和第二储层,第一储层和第二储层的产气方式不同;
第一储层含气量确定模块,用于根据第一储层的岩石属性信息,确定第一储层的第一初始含气量;
第二储层含气量确定模块,用于根据第二储层的岩石属性信息,确定第二储层的第二初始含气量;
气体窜流量确定模块,用于根据第一储层和第二储层的岩石属性信息,确定第一储层与第二储层之间的气体窜流量;
井筒气液两相流压降确定模块,用于确定合采井筒内气液两相流压降;
合采产气量确定模块,用于根据所述第一初始含气量、所述第二初始含气量、所述气体窜流量和所述合采井筒内气液两相流压降,确定待开采气藏的产气量;
其中,第一储层为煤层,第二储层为砂岩层;
所述气体窜流量确定模块还用于通过如下公式,计算第一储层与第二储层之间的气体窜流量:
其中,qgcross表示煤层与砂岩层之间的气体窜流量;pcg表示煤岩层孔隙气相压力;psg表示砂岩层孔隙气相压力;μcg表示煤层中气体的动力粘度;μsg表示砂岩层中气体的动力粘度;ρcg表示煤层条件下气体的密度;ρsg表示砂岩层条件下气体的密度;kcrg表示煤岩层气体相对渗透率;ksrg表示砂岩层气体相对渗透率;hc表示煤层的厚度;hs表示砂岩层的厚度;kcz表示煤岩层垂向渗透率;ksz表示砂岩层垂向渗透率;β1和β2表示基于熔合界面的层间窜流阻力系数;ρsc表示地面条件下气体的密度;Rs表示单位体积水中溶解气量;pcw表示煤岩层孔隙水相压力;psw表示砂岩层孔隙水相压力;μcw表示煤层中水相动力粘度;kcrw表示煤岩层水相对渗透率;μsw表示砂岩层水相粘度;ksrw表示砂岩层水相对渗透率。
8.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至6任一项所述多储层气藏合采产能预测方法。
9.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至6任一项所述多储层气藏合采产能预测方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010907844.XA CN112065339B (zh) | 2020-09-02 | 2020-09-02 | 多储层气藏合采产能预测方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010907844.XA CN112065339B (zh) | 2020-09-02 | 2020-09-02 | 多储层气藏合采产能预测方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112065339A CN112065339A (zh) | 2020-12-11 |
CN112065339B true CN112065339B (zh) | 2021-10-26 |
Family
ID=73665188
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010907844.XA Active CN112065339B (zh) | 2020-09-02 | 2020-09-02 | 多储层气藏合采产能预测方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112065339B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113297803B (zh) * | 2021-06-17 | 2022-12-06 | 东北石油大学 | 一种油气产量智能模拟预估方法和系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6112817A (en) * | 1997-05-06 | 2000-09-05 | Baker Hughes Incorporated | Flow control apparatus and methods |
CN105649578A (zh) * | 2016-01-29 | 2016-06-08 | 西南石油大学 | 煤层气直井双层分压合采方法及其装置 |
CN106481332A (zh) * | 2015-08-31 | 2017-03-08 | 中国石油化工股份有限公司 | 用于确定页岩气多段压裂水平井内外区动态储量的方法 |
CN109826621A (zh) * | 2019-01-17 | 2019-05-31 | 西安科技大学 | 一种煤层气多层合采气水两相渗流实验装置及测试方法 |
CN110929971A (zh) * | 2018-09-20 | 2020-03-27 | 中国石油化工股份有限公司 | 一种多层叠置致密气藏开发方法及井型选择方法 |
CN111209646A (zh) * | 2018-11-20 | 2020-05-29 | 中国石油化工股份有限公司 | 一种致密低渗砂岩气藏多层合采井产气量劈分方法和装置 |
-
2020
- 2020-09-02 CN CN202010907844.XA patent/CN112065339B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6112817A (en) * | 1997-05-06 | 2000-09-05 | Baker Hughes Incorporated | Flow control apparatus and methods |
CN106481332A (zh) * | 2015-08-31 | 2017-03-08 | 中国石油化工股份有限公司 | 用于确定页岩气多段压裂水平井内外区动态储量的方法 |
CN105649578A (zh) * | 2016-01-29 | 2016-06-08 | 西南石油大学 | 煤层气直井双层分压合采方法及其装置 |
CN110929971A (zh) * | 2018-09-20 | 2020-03-27 | 中国石油化工股份有限公司 | 一种多层叠置致密气藏开发方法及井型选择方法 |
CN111209646A (zh) * | 2018-11-20 | 2020-05-29 | 中国石油化工股份有限公司 | 一种致密低渗砂岩气藏多层合采井产气量劈分方法和装置 |
CN109826621A (zh) * | 2019-01-17 | 2019-05-31 | 西安科技大学 | 一种煤层气多层合采气水两相渗流实验装置及测试方法 |
Non-Patent Citations (1)
Title |
---|
煤层气/砂岩气混合气藏合采产能预测及排采优化;郭肖等;《科学技术与工程》;20190630;第19卷(第17期);162-167 * |
Also Published As
Publication number | Publication date |
---|---|
CN112065339A (zh) | 2020-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Javadpour et al. | Nanoscale gas flow in shale gas sediments | |
US20100154514A1 (en) | Method of determining the evolution of petrophysical properties of a rock during diagenesis | |
Alfi et al. | How to improve our understanding of gas and oil production mechanisms in liquid-rich shale | |
Radwan et al. | Gas adsorption and reserve estimation for conventional and unconventional gas resources | |
WO2016028564A1 (en) | Methods for monitoring fluid flow and transport in shale gas reservoirs | |
Yanze et al. | The role of diffusion for nonequilibrium gas injection into a fractured reservoir | |
Esene et al. | A modeling strategy to investigate carbonated water injection for EOR and CO2 sequestration | |
Mjaavatten et al. | A model for gas coning and rate-dependent gas/oil ratio in an oil-rim reservoir | |
CN112065339B (zh) | 多储层气藏合采产能预测方法及装置 | |
Huang et al. | Effects of water invasion law on gas wells in high temperature and high pressure gas reservoir with a large accumulation of water-soluble gas | |
Wang et al. | Experimental investigation of reservoir fluid interlayer crossflow through fracture during the drainage stage of coal measure gas well | |
Christman et al. | Comparison of laboratory-and field-observed CO2 tertiary injectivity | |
Haghshenas et al. | Simulation of enhanced recovery using CO2 in a liquid-rich Western Canadian unconventional reservoir: accounting for reservoir fluid adsorption and compositional heterogeneity | |
Sang et al. | Effects of the laminated-structure and mixed wettability on the oil/water relative permeabilities and oil productions in shale oil formations | |
Bera et al. | Multiphysics gas transport in nanoporous unconventional reservoirs: Challenges of mathematical modelling | |
Bi et al. | Compositional Simulation of Cyclic Gas Injection in Liquid-Rich Shale Reservoirs Using Existing Simulators with a Framework for Incorporating Nanopores | |
CN108798658A (zh) | 一种石油运移参数的确定方法及装置 | |
Ming et al. | A new model for improving the prediction of liquid loading in horizontal gas wells | |
Lopez et al. | Sorption-dependent permeability of shales | |
Fazelipour | Predicting asphaltene precipitation in oilfields via the technology of compositional reservoir simulation | |
Salimi et al. | The influence of heterogeneity, wetting, and viscosity ratio on oil recovery from vertically fractured reservoirs | |
Izadi et al. | An investigation of mechanistic foam modeling for optimum field development of CO2 foam eor application | |
Li et al. | A generalized model for gas flow prediction in shale matrix with deduced coupling coefficients and its macroscopic form based on real shale pore size distribution experiments | |
Alagoz et al. | Calculation of Bottomhole Pressure in Two-Phase Wells Using Beggs and Brill Method: Sensitivity Analysis | |
Cho et al. | Effects of fluid flow in a porous cuttings-bed on cuttings transport efficiency and hydraulics |
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 |