CN117993228A - 缝洞型储气库注采相渗滞后分析方法、设备及存储介质 - Google Patents
缝洞型储气库注采相渗滞后分析方法、设备及存储介质 Download PDFInfo
- Publication number
- CN117993228A CN117993228A CN202410397681.3A CN202410397681A CN117993228A CN 117993228 A CN117993228 A CN 117993228A CN 202410397681 A CN202410397681 A CN 202410397681A CN 117993228 A CN117993228 A CN 117993228A
- Authority
- CN
- China
- Prior art keywords
- gas storage
- fracture
- reservoir
- injection
- type gas
- 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
- 238000003860 storage Methods 0.000 title claims abstract description 151
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 77
- 238000002347 injection Methods 0.000 title claims abstract description 52
- 239000007924 injection Substances 0.000 title claims abstract description 52
- 238000004458 analytical method Methods 0.000 title claims description 14
- 230000035699 permeability Effects 0.000 claims abstract description 57
- 238000000034 method Methods 0.000 claims abstract description 38
- 238000004088 simulation Methods 0.000 claims abstract description 31
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 20
- 238000005213 imbibition Methods 0.000 claims abstract description 15
- 238000006073 displacement reaction Methods 0.000 claims abstract description 14
- 238000001764 infiltration Methods 0.000 claims abstract description 14
- 230000008595 infiltration Effects 0.000 claims abstract description 14
- 238000002474 experimental method Methods 0.000 claims abstract description 13
- 238000009736 wetting Methods 0.000 claims description 23
- 239000004047 hole gas Substances 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 5
- 238000009826 distribution Methods 0.000 claims description 4
- 238000012986 modification Methods 0.000 claims description 4
- 230000004048 modification Effects 0.000 claims description 4
- 239000011435 rock Substances 0.000 claims description 4
- 239000008398 formation water Substances 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 239000007789 gas Substances 0.000 abstract description 213
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 abstract description 12
- 239000003345 natural gas Substances 0.000 abstract description 6
- 230000008569 process Effects 0.000 abstract description 5
- 239000003209 petroleum derivative Substances 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 7
- 230000006872 improvement Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 230000007423 decrease Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000000605 extraction Methods 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 241001415846 Procellariidae Species 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000003204 osmotic effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/21—Design, administration or maintenance of databases
- G06F16/211—Schema design and management
- G06F16/212—Schema design and management with details for data modelling support
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/203—Drawing of straight lines or curves
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Underground Structures, Protecting, Testing And Restoring Foundations (AREA)
Abstract
本发明涉及石油与天然气工程地下气体储存技术领域,公开了缝洞型储气库注采相渗滞后分析方法、设备及介质,缝洞型储气库注采相渗滞后分析方法包括:通过缝洞型储气库岩心多轮次气水相渗实验得到多周期渗吸驱替后的相对渗透率数据,用于生成缝洞型储气库的相渗滞后曲线;基于动态参数构建缝洞型储气库地质模型和数值模拟模型;计算相渗滞后曲线影响下多周期地下气体体积;评价相渗滞后对缝洞型储气库注采影响程度,调整缝洞型储气库运行方案中的单周期采气总量和各生产井配产。采用本方案能够解决目前缝洞型储气库乃至常规储气库运行模拟过程中忽略相渗滞后导致的误差问题,从而为储气库运行方案的制订提供指导。
Description
技术领域
本发明涉及石油与天然气工程地下气体储存技术领域,具体涉及一种缝洞型储气库注采相渗滞后分析方法、设备及存储介质。
背景技术
在一年的不同时期,天然气的供需关系也各有不同。建立天然气地下储气库能够起到调峰作用,是解决不同时期天然气供需不平衡的有效方法。缝洞型储气库具有储集能力强、采气速度快的特点;缝洞型储气库由于其多周期循环注采的特点,常出现相渗滞后等渗流机理,相渗滞后即交替进行渗吸驱替时由于毛管力等原因,更多的气体被困在孔隙中无法被采出,降低了缝洞型储气库的工作气量,进而影响到运行方案中制定的注采速度,因此在制定储气库运行方案时需要考虑相渗滞后对于注采的影响。
储气库数值模拟方法主要用于指导储气库运行方案的制订或计算某些关键参数。现有储气库数值模拟方法主要针对储气库储气空间计算或气井产气量和注气量计算,用以指导运行方案制订,但大多基于常规类型储气库中的基质型或基质-裂缝型储气库,未考虑裂缝和溶洞对储气库注采的影响,且在进行数值模拟时未考虑相渗滞后这一渗流机理,给储气库的运行模拟带来了很大的误差。缝洞型储气库相较于常规储气库,由于其裂缝和溶洞的存在,储层的平均孔隙度和平均渗透率都更高,储气空间和工作气量更大。因此在进行多周期注采时会有更多的气体滞留在储层中难以采出,导致忽略相渗滞后所产生的误差相较于常规储气库更为显著。因此,为了准确模拟缝洞型储气库的运行过程,需要厘清相渗滞后这一渗流机理对储气库注采的影响,从而为指导运行方案的设计提供技术支持。
发明内容
本发明意在提供缝洞型储气库注采相渗滞后分析方法、设备及存储介质,以解决目前缝洞型储气库乃至常规储气库运行模拟过程中忽略相渗滞后导致的误差问题。
为达到上述目的,本发明采用如下技术方案:
缝洞型储气库注采相渗滞后分析方法,包括:
S1,通过缝洞型储气库岩心多轮次气水相渗实验得到多周期渗吸驱替后的相对渗透率数据;
S2,获取建库目标气藏当前状态下的动态参数;
S3,基于所述的动态参数构建缝洞型储气库地质模型和数值模拟模型,并完成模型初始化,用于后续运行预测;
S4,根据实验测得的相对渗透率数据生成缝洞型储气库的相渗滞后曲线;
S5,模拟考虑相渗滞后和不考虑相渗滞后的缝洞型储气库运行,将相渗滞后曲线导入数值模拟模型中,一种方案使用导入的相渗滞后曲线,另一种方案使用原有相对渗透率曲线,分别进行缝洞型储气库多周期注采运行预测,计算出两种情况下缝洞型储气库多周期的地下气体体积;
S6,对所述地下气体体积进行分析处理,评价相渗滞后对缝洞型储气库注采影响程度。
本方案的原理及优点是:
1. 本发明通过在数值模拟过程中加入相渗滞后这一重要机理,降低了储气库运行模型过程中的误差,模拟结果更加准确。
2. 本发明将相对渗透率实验和数值模拟进行了结合,根据实验数据指导数值模拟参数的设置,提高了数值模拟相较于矿场实际情况的准确性,为储气库运行方案的制订提供了有利的支撑。
优选的,作为一种改进,所述动态参数包括:气藏建库区域范围、气藏圈闭线和断层线、气藏深度、气藏裂缝分布、气藏溶洞参数、气藏孔隙度场、气藏渗透率场、气藏饱和度场、气藏净毛比场、气藏温度场、气藏压力场、气藏气体组成、气藏地层水特征、气藏岩石特征以及气藏相对渗透率曲线。
技术效果:通过上述动态参数,能够使得地质模型和数值模拟模型更加真实的模拟现场实际情况,从而使得模拟的生产动态及提供的运行方案更加真实可靠。
优选的,作为一种改进,所述气水相渗实验使用高温高压油气水相渗测试仪完成,通过气水相渗实验仪器模拟储气库运行时多周期循环注采过程进行多轮次的两相流体相对渗透率测定。
技术效果:通过高温高压油气水相渗测试仪能够对真实岩心中的气水相对渗透率进行准确测试,测试精度高。
优选的,作为一种改进,所述S3中,采用软件Petrel构建缝洞型储气库的地质模型,采用数值模拟器INTERSECT构建缝洞型储气库的数值模拟模型。
技术效果:通过软件Petrel和数值模拟器INTERSECT,能够高精度的进行建模。
优选的,作为一种改进,所述S4中在处理相对渗透率滞后时,首先利用历史最大残余非润湿相饱和度关系得到渗吸曲线的终点;然后将该终点值带入预设方法中,计算得到渗吸曲线;
历史最大残余非润湿相饱和度关系及预设方法具体为:
其中,Sncrt为开始渗吸时达到的非润湿相残余饱和度,Sncrd为驱替过程中非润湿相临界流动饱和度,A为修改参数,Ct为Land常数,Snhy为驱替过程中非润湿相临界流动饱和度,b为Jerauld参数,Snmx为最大非润湿相饱和度,Sncri为临界非润湿相饱和度;
相对渗透率变量关系式为:
其中,krn为饱和度Sn对应的相对渗透率,krni为与渗吸毛细管压力曲线相对应的相对渗透率,krnd为与驱替毛细管压力曲线相对应的相对渗透率,Snorm为归一化气体饱和度。
技术效果:通过所述预设方法,所需数据少的特点,计算精度高,便于提高相渗滞后曲线生成精度。
优选的,作为一种改进,所述S6中对计算数据进行分析处理的目标包括:两种情况下对应的各周期储气库库容量、各周期储气库工作气量和各周期工作气比例。
技术效果:通过上述目标,便于对储气库的重要参数进行分析。
优选的,作为一种改进,所述S6中相渗滞后对缝洞型储气库注采影响程度包括相对渗透率滞后导致的储气库单周期注气能力和采气能力的下降程度和库容量下降程度。
技术效果:厘清单周期注气能力和采气能力有助于制订准确的运行方案保证储气库按计划运行,满足用气高峰期居民的用气需求。
一种电子设备,包括存储器和处理器,所述存储器用于存储计算机程序,该计算机程序被所述处理器执行时,执行所述的缝洞型储气库注采相渗滞后分析方法。
一种存储介质,该存储介质存储的计算机程序,能够被一个或多个处理器执行,能够用来实现所述的缝洞型储气库注采相渗滞后分析方法。
附图说明
图1为缝洞型储气库注采相渗滞后分析方法实施例的流程图;
图2为多轮次气水相渗实验流程示意图;
图3为缝洞型储气库2D地质模型示意图;
图4为缝洞型储气库3D地质模型示意图;
图5为缝洞型储气库数值模拟模型及圈闭范围示意;
图6为相渗滞后曲线图;
图7为两种算例地下气体总量图。
具体实施方式
下面通过具体实施方式进一步详细说明:
实施例基本如附图1所示:
缝洞型储气库注采相渗滞后分析方法,包括:
S1,通过缝洞型储气库岩心多轮次气水相渗实验得到多周期渗吸驱替后的相对渗透率数据。
所述相对渗透率数据包括岩心孔隙度、气相渗透率及对应含气饱和度、水相渗透率及对应的含水饱和度、等渗点饱和度和气水共渗区范围。
本实施例中,缝洞型储气库岩心多轮次气水相渗实验使用高温高压油气水相渗测试仪完成,具体型号为ISCO、HAG-250;该仪器能够模拟地层高温高压条件,对真实岩心中的气水相对渗透率进行准确测试,测试精度高。依据《GB/T 28912-2012岩石中两相流体相对渗透率测定方法》进行多轮次的两相流体相对渗透率测定,模拟储气库运行时多周期循环注采过程。多轮次气水相渗实验流程如图2所示。
S2,获取建库目标气藏当前状态下的动态参数。
本实施例中,建库目标气藏当前状态下的动态参数包括气藏建库区域范围、气藏圈闭线和断层线、气藏深度、气藏裂缝分布、气藏溶洞参数、气藏孔隙度场、气藏渗透率场、气藏饱和度场、气藏净毛比场、气藏温度场、气藏压力场、气藏气体组成、气藏地层水特征、气藏岩石特征以及气藏相对渗透率曲线。
S3,基于所述的动态参数构建缝洞型储气库地质模型和数值模拟模型,并完成模型初始化,用于后续运行预测。
将动态参数以文件输入或手动输入到数值模拟软件中,即能够生成对应的地质模型和数值模拟模型。所述地质模型和数值模拟模型通过历史拟合以接近现场实际情况,所模拟的生产动态及提供的运行方案更加真实可靠。相较于常规的基质型储气库或基质-裂缝型储气库,本实施例中额外考虑了溶洞对地质模型属性的影响,反映了缝洞型储气库溶洞发育的特点,溶洞是缝洞型储气库的重要储气空间,溶洞的存在使得缝洞型储气库能够储存更多的天然气。除此之外,溶洞还提高了储层的平均孔隙度和平均渗透率,增大了气体的流动速度,忽略溶洞会导致生成的模型存在巨大误差,进而影响储气库的运行方案制订。在生成用于计算的孔隙度场和渗透率场时,还考虑了溶洞参数的地质模型拥有更高的孔隙度和渗透率,在后续计算中,溶洞的存在还会进一步影响模型的饱和度分布。
本实施例中,构建缝洞型储气库地质模型和数值模拟模型时使用的软件分别为地质工程一体化软件Petrel、数值模拟器INTERSECT。Petrel用于国内外石油地质领域的地质建模及后续的数值模拟过程,作为成熟的商业软件,建模精度高,功能齐全,便于提高本方案中模型的精度。数值模拟器INTERSECT为Petrel配套数值模拟器,运算速度快,与Petrel建立的地质模型更加匹配。
获取实际开发区块的动态参数建立考虑裂缝、溶洞的缝洞型储气库地质模型,如图3、图4所示,数值模拟模型如图5所示。
S4,根据实验测得的相对渗透率数据生成缝洞型储气库的相渗滞后曲线。相渗滞后曲线特征为相对渗透率曲线中出现额外一根或多根气相曲线,且相较于原曲线有更高的残余气饱和度。
以某缝洞型气田改建储气库数值模拟为例,测定得到模拟多周期注采时每个周期的相对渗透率曲线,如图6所示,读取最后一个周期相渗曲线中气相曲线和x轴的交点取值,即考虑多周期注采相渗滞后时的最大残余气饱和度约为0.3,最大残余气饱和度是相对渗透率曲线和x坐标轴的交点取值。
在处理相对渗透率滞后时,首先利用历史最大残余非润湿相饱和度关系得到渗吸曲线的终点;然后将该终点值,即最大非润湿相饱和度带入预设方法对应公式中即可计算出渗吸曲线,即考虑了相渗滞后的相对渗透率曲线数据。
历史最大残余非润湿相饱和度关系及预设方法具体为:
其中,Sncrt为开始渗吸时达到的非润湿相残余饱和度,Sncrd为驱替过程中非润湿相临界流动饱和度,A为修改参数,Ct为Land常数,Snhy为驱替过程中非润湿相临界流动饱和度,b为Jerauld参数,Snmx为最大非润湿相饱和度,Sncri为临界非润湿相饱和度;
相对渗透率变量关系式为:
其中,krn为饱和度Sn对应的相对渗透率,krni为与渗吸毛细管压力曲线相对应的相对渗透率,krnd为与驱替毛细管压力曲线相对应的相对渗透率,Snorm为归一化气体饱和度。
多轮次的两相流体相对渗透率测得的最大残余气饱和度根据预设方法生成多周期相渗滞后曲线,如图6所示。
本实施例中,生成的缝洞型储气库的相渗滞后曲线使用的预设方法已被编写在数值模拟器INTERSECT中,在使用时调用即可。该预设方法具有所需数据少的特点,计算精度高,只需输入最大非润湿相饱和度即可进行下一步运算。
S5,模拟考虑相渗滞后和不考虑相渗滞后的缝洞型储气库运行,具体的,将相渗滞后曲线导入数值模拟模型中,一种方案使用导入的相渗滞后曲线,另一种方案使用原有相对渗透率曲线,分别进行缝洞型储气库多周期注采运行预测,计算出两种情况下缝洞型储气库多周期地下气体体积。其中使用了相渗滞后曲线的方案由于残余气饱和度更高,会导致更多的气无法被采出,在地下气体体积图中,曲线的峰值更低。如图7所示,考虑相渗滞后的曲线峰值更低。
设置工作制度为日产气量150万方/天,最低井底压力6.5MPa,分别进行考虑相渗滞后和不考虑相渗滞后的储气库多周期注采数值模拟,计算出的地下气体总量对比如图7所示。
S6,对所述地下气体体积进行分析处理,评价相渗滞后对缝洞型储气库注采影响程度,并根据评价结果调整缝洞型储气库运行方案中的单周期采气总量和各生产井配产。
本实施例中,对计算数据进行分析处理的目标包括:考虑相渗滞后和不考虑相渗滞后两种情况下对应的各周期储气库库容量、各周期储气库工作气量(单周期最大注气能力和单周期最大采气能力)和各周期工作气比例。其中,各周期储气库库容量为:地下气体体积曲线中对应的地下最大气体体积;各周期储气库工作气量为:地下气体体积曲线中每个周期的最大值减去最小值,其中单周期最大注气能力为地下气体体积最大值减去上一周期最小值,单周期最大采气能力为地下气体体积最大值减去下一周期最小值;各周期工作气比例为:地下气体体积曲线中每个周期的最大值减去最小值后除以地下气体体积曲线中对应的地下最大气体体积。
库容量代表了地下储存气体的体积,库容量越大,储气库一个运行周期所能采出和注入的气体就越多,是储气库运行的一个重要参数。工作气量代表了储气库一个运行周期所能采出和注入的气体,储气库运行方案中的注气速度及总量、采气速度及总量和单井配产都受到工作气量大小的影响。相渗滞后导致部分气体滞留在孔隙中难以被采出,在储气库下一个注气阶段,这部分气体被驱向地层远处,难以在下个采气阶段被采出从而变成了储气库垫底气,影响了储气库的经济效益和运行成本。工作气比例为工作气量和库容量的比值,能够侧面表征这部分新增垫底气的体积,从而有助于制订不同储气库运行周期的运行方案和评估经济效益和运行成本。
相渗滞后对缝洞型储气库注采影响程度包括相对渗透率滞后导致的储气库单周期注气能力和采气能力的下降程度和库容量下降程度。根据计算出的注气能力和采气能力的下降程度调整储气库运行方案中的注采速度,根据库容量的降低程度降低储气库运行方案中一个采气周期的计划采气总量。考虑到缝洞型储气库采气前期主要通过裂缝和溶洞供气,气体输运能力较强,而后期主要开采基质中的气体,采气速度和采气能力将出现大幅下降。除此之外,相较于常规储气库,缝洞型储气库由于其裂缝较多,在采气过程中存在更强的渗透率敏感,会使得渗透率随采气过程不断降低,进一步降低气井的采气能力。因此在调整储气库运行方案时,各生产井配产的降低程度应大于计算出的储气库单周期采气能力的下降程度,防止过早出现产量快速衰减,保证储气库的平稳运行。
库容量是储气库运行的一个重要参数,代表了地下储存气体的体积,库容量越大,储气库一个运行周期所能采出和注入的气体就越多。而储气库单周期注气能力和采气能力决定了储气库运行方案制订时的单井配产,厘清单周期注气能力和采气能力有助于制订准确的运行方案保证储气库按计划运行,满足用气高峰期居民的用气需求。现有的储气库运行数值模拟方法忽略了相渗滞后,计算出的采气能力偏大,在实际生产时会导致气井产量提前衰竭,不利于储气库的平稳运行。
一种电子设备,包括存储器和处理器,所述存储器用于存储计算机程序,该计算机程序被所述处理器执行时,执行所述的缝洞型储气库注采相渗滞后分析方法。
一种存储介质,该存储介质存储的计算机程序,能够被一个或多个处理器执行,能够用来实现所述的缝洞型储气库注采相渗滞后分析方法。
以上所述的仅是本发明的实施例,方案中公知的具体技术方案和/或特性等常识在此未作过多描述。应当指出,对于本领域的技术人员来说,在不脱离本发明技术方案的前提下,还可以作出若干变形和改进,这些也应该视为本发明的保护范围,这些都不会影响本发明实施的效果和专利的实用性。本申请要求的保护范围应当以其权利要求的内容为准,说明书中的具体实施方式等记载可以用于解释权利要求的内容。
Claims (9)
1.缝洞型储气库注采相渗滞后分析方法,其特征在于,包括:
S1,通过缝洞型储气库岩心多轮次气水相渗实验得到多周期渗吸驱替后的相对渗透率数据;
S2,获取建库目标气藏当前状态下的动态参数;
S3,基于所述的动态参数构建缝洞型储气库地质模型和数值模拟模型,并完成模型初始化,用于后续运行预测;
S4,根据实验测得的相对渗透率数据生成缝洞型储气库的相渗滞后曲线;
S5,模拟考虑相渗滞后和不考虑相渗滞后的缝洞型储气库运行,将相渗滞后曲线导入数值模拟模型中,一种方案使用导入的相渗滞后曲线,另一种方案使用原有相对渗透率曲线,分别进行缝洞型储气库多周期注采运行预测,计算出两种情况下缝洞型储气库多周期的地下气体体积;
S6,对所述地下气体体积进行分析处理,评价相渗滞后对缝洞型储气库注采影响程度。
2.根据权利要求1所述的缝洞型储气库注采相渗滞后分析方法,其特征在于,所述动态参数包括:气藏建库区域范围、气藏圈闭线和断层线、气藏深度、气藏裂缝分布、气藏溶洞参数、气藏孔隙度场、气藏渗透率场、气藏饱和度场、气藏净毛比场、气藏温度场、气藏压力场、气藏气体组成、气藏地层水特征、气藏岩石特征以及气藏相对渗透率曲线。
3.根据权利要求1所述的缝洞型储气库注采相渗滞后分析方法,其特征在于:所述气水相渗实验使用高温高压油气水相渗测试仪完成,通过气水相渗实验仪器模拟储气库运行时多周期循环注采过程进行多轮次的两相流体相对渗透率测定。
4.根据权利要求1所述的缝洞型储气库注采相渗滞后分析方法,其特征在于:所述S3中,采用软件Petrel构建缝洞型储气库的地质模型,采用数值模拟器INTERSECT构建缝洞型储气库的数值模拟模型。
5.根据权利要求1所述的缝洞型储气库注采相渗滞后分析方法,其特征在于:所述S4中在处理相对渗透率滞后时,首先利用历史最大残余非润湿相饱和度关系得到渗吸曲线的终点;然后将该终点值带入预设方法中,计算得到渗吸曲线;
历史最大残余非润湿相饱和度关系及预设方法具体为:
其中,Sncrt为开始渗吸时达到的非润湿相残余饱和度,Sncrd为驱替过程中非润湿相临界流动饱和度,A为修改参数,Ct为Land常数,Snhy为驱替过程中非润湿相临界流动饱和度,b为Jerauld参数,Snmx为最大非润湿相饱和度,Sncri为临界非润湿相饱和度;
相对渗透率变量关系式为:
其中,krn为饱和度Sn对应的相对渗透率,krni为与渗吸毛细管压力曲线相对应的相对渗透率,krnd为与驱替毛细管压力曲线相对应的相对渗透率,Snorm为归一化气体饱和度。
6.根据权利要求1所述的缝洞型储气库注采相渗滞后分析方法,其特征在于,所述S6中对计算数据进行分析处理的目标包括:两种情况下对应的各周期储气库库容量、各周期储气库工作气量和各周期工作气比例。
7.根据权利要求1所述的缝洞型储气库注采相渗滞后分析方法,其特征在于,所述S6中相渗滞后对缝洞型储气库注采影响程度包括相对渗透率滞后导致的储气库单周期注气能力和采气能力的下降程度和库容量下降程度。
8.一种电子设备,其特征在于:包括存储器和处理器,所述存储器用于存储计算机程序,该计算机程序被所述处理器执行时,执行如权利要求1-7任意一项所述的缝洞型储气库注采相渗滞后分析方法。
9.一种存储介质,其特征在于:该存储介质存储的计算机程序,能够被一个或多个处理器执行,能够用来实现如权利要求1至7任意一项所述的缝洞型储气库注采相渗滞后分析方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410397681.3A CN117993228B (zh) | 2024-04-03 | 2024-04-03 | 缝洞型储气库注采相渗滞后分析方法、设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410397681.3A CN117993228B (zh) | 2024-04-03 | 2024-04-03 | 缝洞型储气库注采相渗滞后分析方法、设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117993228A true CN117993228A (zh) | 2024-05-07 |
CN117993228B CN117993228B (zh) | 2024-05-31 |
Family
ID=90895445
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410397681.3A Active CN117993228B (zh) | 2024-04-03 | 2024-04-03 | 缝洞型储气库注采相渗滞后分析方法、设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117993228B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109356564A (zh) * | 2018-11-08 | 2019-02-19 | 中国海洋石油集团有限公司 | 一种基于三元动态参数的水驱油相渗时变规律表征方法 |
US20200202060A1 (en) * | 2017-05-15 | 2020-06-25 | University Of Houston System | Advanced technique for screening enhanced oil recovery and improved oil recovery methodologies for a petroleum reservoir |
CN112943176A (zh) * | 2021-02-01 | 2021-06-11 | 重庆科技学院 | 模拟油藏型储气库注采的气油相对渗透率曲线测定方法 |
CN113008753A (zh) * | 2021-02-01 | 2021-06-22 | 重庆科技学院 | 一种应用二维实验优化储气库气水相渗曲线的方法 |
CN114372352A (zh) * | 2021-12-23 | 2022-04-19 | 成都理工大学 | 渗流-温度双场耦合数值模拟预测复杂断块油藏储气库调峰能力方法 |
CN114880962A (zh) * | 2021-12-15 | 2022-08-09 | 成都理工大学 | 基于地层渗流-井筒管流耦合的致密气藏储气库单井注采能力一体化分析方法 |
CN114897436A (zh) * | 2022-06-07 | 2022-08-12 | 国家石油天然气管网集团有限公司 | 用于储气库的地层压力平衡方法和控制器 |
CN116306385A (zh) * | 2023-05-24 | 2023-06-23 | 中国石油大学(华东) | 一种油藏压裂渗吸增能数值模拟方法、系统、设备及介质 |
CN117129383A (zh) * | 2023-08-23 | 2023-11-28 | 成都理工大学 | 一种模拟储氢库储层中氢气水渗流滞后规律的装置及方法 |
-
2024
- 2024-04-03 CN CN202410397681.3A patent/CN117993228B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20200202060A1 (en) * | 2017-05-15 | 2020-06-25 | University Of Houston System | Advanced technique for screening enhanced oil recovery and improved oil recovery methodologies for a petroleum reservoir |
CN109356564A (zh) * | 2018-11-08 | 2019-02-19 | 中国海洋石油集团有限公司 | 一种基于三元动态参数的水驱油相渗时变规律表征方法 |
CN112943176A (zh) * | 2021-02-01 | 2021-06-11 | 重庆科技学院 | 模拟油藏型储气库注采的气油相对渗透率曲线测定方法 |
CN113008753A (zh) * | 2021-02-01 | 2021-06-22 | 重庆科技学院 | 一种应用二维实验优化储气库气水相渗曲线的方法 |
CN114880962A (zh) * | 2021-12-15 | 2022-08-09 | 成都理工大学 | 基于地层渗流-井筒管流耦合的致密气藏储气库单井注采能力一体化分析方法 |
CN114372352A (zh) * | 2021-12-23 | 2022-04-19 | 成都理工大学 | 渗流-温度双场耦合数值模拟预测复杂断块油藏储气库调峰能力方法 |
CN114897436A (zh) * | 2022-06-07 | 2022-08-12 | 国家石油天然气管网集团有限公司 | 用于储气库的地层压力平衡方法和控制器 |
CN116306385A (zh) * | 2023-05-24 | 2023-06-23 | 中国石油大学(华东) | 一种油藏压裂渗吸增能数值模拟方法、系统、设备及介质 |
CN117129383A (zh) * | 2023-08-23 | 2023-11-28 | 成都理工大学 | 一种模拟储氢库储层中氢气水渗流滞后规律的装置及方法 |
Non-Patent Citations (6)
Title |
---|
JINZI LIU等: "Coalbed Gas Reservoir Using Artificial Intelligent Algorithm of Calculating Relative Permeability Curve Method", 《2012 INTERNATIONAL CONFERENCE ON COMPUTER DISTRIBUTED CONTROL AND INTELLIGENT ENVIRONMENTAL MONITORING》, 5 April 2012 (2012-04-05), pages 191 - 194 * |
孙军昌等: "气藏型地下储气库建库注采机理与评价关键技术", 《天然气工业》, vol. 38, no. 04, 20 April 2018 (2018-04-20), pages 138 - 144 * |
杨羿: "多层系含硫储气库硫化氢淘洗规律数值模拟研究", 《第33届全国天然气学术年会(2023)论文集》, 1 June 2023 (2023-06-01) * |
汪周华等: "缝洞型碳酸盐岩气藏高温高压气水两相渗流特征", 《油气藏评价与开发》, vol. 7, no. 02, 26 April 2017 (2017-04-26), pages 47 - 52 * |
石磊等: "水侵气藏型储气库气水微观渗流规律", 《天然气勘探与开发》, vol. 43, no. 01, 25 March 2020 (2020-03-25), pages 58 - 63 * |
肖尊荣等: "基于双重介质嵌入式离散裂缝模型的致密油藏产能影响因素", 《科学技术与工程》, vol. 23, no. 25, 30 September 2023 (2023-09-30), pages 10780 - 10790 * |
Also Published As
Publication number | Publication date |
---|---|
CN117993228B (zh) | 2024-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11782741B2 (en) | Modeling of fluid introduction and/or fluid extraction elements in simulation of coreflood experiment | |
CN104750896B (zh) | 一种缝洞型碳酸盐岩油藏数值模拟方法 | |
CN110824556B (zh) | 一种非常规致密砂岩储层的岩石物理模型建立方法及应用 | |
CN101892837B (zh) | 地层因数确定方法及含油饱和度确定方法 | |
US20060069511A1 (en) | Gas reservoir evaluation and assessment tool method and apparatus and program storage device | |
CN104879103A (zh) | 一种分层注水效果分析方法 | |
RU2002122397A (ru) | Всеобъемлющая оптимизация пластового резервуара | |
CN105804737A (zh) | 一种基于迭代算法求解地层孔隙度的方法 | |
CN107808032A (zh) | 页岩气储层含水饱和度计算方法及系统 | |
CN111706318B (zh) | 一种确定低渗储层剩余油分布状况的方法 | |
CN114066084B (zh) | 一种基于机器学习的相渗曲线预测方法及系统 | |
CN107130952A (zh) | 一种三级复合石油增产方法 | |
CN114880962A (zh) | 基于地层渗流-井筒管流耦合的致密气藏储气库单井注采能力一体化分析方法 | |
Kolin et al. | Pressure build-up test analysis of the reservoir system with the multiphase flow | |
CN113933899A (zh) | 基于地震属性约束指导下的砂砾岩储层建模方法和系统 | |
CN117993228B (zh) | 缝洞型储气库注采相渗滞后分析方法、设备及存储介质 | |
CN111155980B (zh) | 一种水流优势通道识别方法及装置 | |
Lian et al. | Evaluation and applicability study on prediction methods of water inflow in mines | |
Zhang et al. | Storage ratio of CO2 miscible flooding in Chang 8 reservoir of H block in ordos basin under different injection methods | |
Lamas et al. | Application of a global interpretation model for assessment of the stress field for engineering purposes | |
Xiaoqing et al. | Prediction method of produced polymer concentration based on interwell connectivity | |
CN111119841A (zh) | 一种盐岩组分精细处理与岩性自动识别的评价方法 | |
Li et al. | Application of Petrel software in reserve estimation in Xinli-Xinbei oilfield, Jilin Province, China | |
CN115062552B (zh) | 一种缝洞油藏氮气吞吐效果预测方法及系统 | |
CN114810012B (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 |