CN113484222B - 致密气储层渗透率瓶颈区确定方法及相关设备 - Google Patents
致密气储层渗透率瓶颈区确定方法及相关设备 Download PDFInfo
- Publication number
- CN113484222B CN113484222B CN202110899572.8A CN202110899572A CN113484222B CN 113484222 B CN113484222 B CN 113484222B CN 202110899572 A CN202110899572 A CN 202110899572A CN 113484222 B CN113484222 B CN 113484222B
- Authority
- CN
- China
- Prior art keywords
- permeability
- gas reservoir
- sandstone
- pore
- tight 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.)
- Active
Links
- 230000035699 permeability Effects 0.000 title claims abstract description 133
- 238000000034 method Methods 0.000 title claims abstract description 46
- 239000011148 porous material Substances 0.000 claims abstract description 132
- 238000006073 displacement reaction Methods 0.000 claims abstract description 27
- 239000007789 gas Substances 0.000 claims description 122
- 238000002474 experimental method Methods 0.000 claims description 51
- 239000012071 phase Substances 0.000 claims description 47
- 238000004088 simulation Methods 0.000 claims description 31
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 26
- 239000011435 rock Substances 0.000 claims description 25
- 238000004590 computer program Methods 0.000 claims description 22
- 238000012360 testing method Methods 0.000 claims description 14
- 239000007788 liquid Substances 0.000 claims description 13
- 229910052757 nitrogen Inorganic materials 0.000 claims description 13
- 238000001179 sorption measurement Methods 0.000 claims description 13
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 12
- 229910052753 mercury Inorganic materials 0.000 claims description 12
- 238000003860 storage Methods 0.000 claims description 12
- 238000005481 NMR spectroscopy Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 10
- 239000011800 void material Substances 0.000 claims description 10
- 238000005119 centrifugation Methods 0.000 claims description 7
- 239000008346 aqueous phase Substances 0.000 claims description 5
- 238000009736 wetting Methods 0.000 claims description 5
- 238000002347 injection Methods 0.000 claims description 2
- 239000007924 injection Substances 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 abstract description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 23
- 238000010586 diagram Methods 0.000 description 13
- 238000005325 percolation Methods 0.000 description 8
- 239000000243 solution Substances 0.000 description 8
- 238000010219 correlation analysis Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 241001415846 Procellariidae Species 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000013211 curve analysis Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 239000012153 distilled water Substances 0.000 description 1
- 238000011234 economic evaluation Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 230000005514 two-phase flow Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
-
- 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
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种致密气储层渗透率瓶颈区确定方法及相关设备。该方法包括:获取所述致密气储层的砂岩微观孔隙结构特征,获取所述致密气储层的渗透率瓶颈区边界条件,基于所述砂岩微观孔隙结构特征和瓶颈区边界条件确定所述致密气储层的渗透率瓶颈区。通过上述方法能够用确定利用砂岩的比表面积、平均孔喉半径、特征结构参数、排驱压力确定致密气储层渗透率瓶颈区,对指导油气生产具有重要意义。
Description
技术领域
本说明书涉及油气领域,更具体地说,本发明涉及一种致密气储层渗透率瓶颈区确定方法及相关设备。
背景技术
在致密砂岩气藏的相对渗透率曲线中,存在气、水两相均无法有效流动的渗透率瓶颈区,在该区域中对气体和水的相对渗透率都非常低,渗透率瓶颈区的确认对指导油气生产具有重要意义。目前国内外对渗透率瓶颈区的研究较少,缺少根据生产实际的渗透率瓶颈区表征。当前确定致密气储层渗透率瓶颈区时,需要工程人员到油田现场进行气水驱替实验确定渗透率瓶颈区,实验较为复杂,影响确定效率。
发明内容
在发明内容部分中引入了一系列简化形式的概念,这将在具体实施方式部分中进一步详细说明。本发明的发明内容部分并不意味着要试图限定出所要求保护的技术方案的关键特征和必要技术特征,更不意味着试图确定所要求保护的技术方案的保护范围。
为至少部分地解决上述问题,第一方面,本发明提出一种致密气储层渗透率瓶颈区确定方法,其特征在于,包括:
获取上述致密气储层的砂岩微观孔隙结构特征;
获取上述致密气储层的渗透率瓶颈区边界条件,其中,上述渗透率瓶颈区边界条件由数值模拟仿真得到的;
基于上述砂岩微观孔隙结构特征和上述渗透率瓶颈区边界条件确定上述致密气储层的渗透率瓶颈区。
可选的,上述获取上述致密气储层的砂岩微观孔隙结构特征,包括:
获取所述致密气储层的砂岩样本;
对上述砂岩样本进行目标实验,得到上述砂岩微观空隙结构特征,其中上述目标实验包括低温液氮吸附实验、核磁共振实验、离心实验和高压压汞实验,上述砂岩微观孔隙结构特征包括孔隙类型、比表面积、总孔体积、特征结构参数、平均孔喉半径、孔喉大小和孔隙连通性。
可选的,上述获取上述致密气储层的砂岩微观孔隙结构特征,包括:
基于上述低温液氮吸附实验获取砂岩的上述孔隙类型、上述比表面积、上述总孔体积和上述特征结构参数;
基于高压压汞试验获取上述砂岩的上述平均孔喉半径和上述孔喉大小;
基于上述离心实验和上述核磁共振实验获取上述孔隙连通性。
可选的,上述方法包括:
用排驱压力表示上述孔喉大小,其中,上述排驱压力是指非润湿相刚进入上述砂岩样本最大孔喉时给予的压力。
可选的,上述获取上述致密气储层的瓶颈区边界条件,包括:
获取多条致密气储层相渗透率曲线;
基于上述致密气储层相渗透率曲线进行数值模拟操作;
基于上述数值模拟操作的结果,确定上述瓶颈区边界条件。
可选的,上述基于上述致密气储层相渗透率曲线进行数值模拟操作,包括:
基于水相在岩样中的相对渗透率、气相在岩样中的相对渗透率及水相和气相的干扰程度将上述致密气储层相渗透率曲线分组,其中,每个分组至少包含一个曲线;
将分组的曲线分别进行均值处理,得到特征曲线;
基于上述特征曲线进行数值模拟。
可选的,上述基于上述砂岩微观孔隙结构特征和瓶颈区边界条件确定上述致密气储层的渗透率瓶颈区,包括:
基于砂岩的上述比表面积、上述平均孔喉半径、上述特征结构参数、上述排驱压力确定上述致密气储层的渗透率瓶颈区。
第二方面,本发明还提出一种致密气储层渗透率瓶颈区确定装置,其特征在于,包括:
第一获取单元,用于获取上述致密气储层的砂岩微观孔隙结构特征;
第二获取单元,用于获取上述致密气储层的渗透率瓶颈区边界条件;
瓶颈区确定单元,用于基于上述砂岩微观孔隙结构特征和上述渗透率瓶颈区边界条件确定上述致密气储层的渗透率瓶颈区。
第三方面,本发明还提出一种电子设备,包括:储存器、处理器以及存储在上述存储器中并可在上述处理器上运行的计算机程序,其特征在于,上述处理器用于执行存储器中存储的计算机程序时实现上述第一方面任一种致密气储层渗透率瓶颈区确定方法的步骤。
第四方面,本发明还提出一种计算机可读存储介质,其上存储有计算机程序,可选的,上述计算机程序被处理器执行时实现上述第一方面任一种致密气储层渗透率瓶颈区确定方法。
综上,通过实验获取致密气储层的砂岩微观孔隙结构特征,对致密气储层相渗透曲线进行分类,并将分类后数据用于仿真确定致密气储层的边界条件,基于微观孔隙结构与致密气储层渗透率瓶颈区相关性分析,用砂岩微观孔隙结构特征确定所述致密气储层的渗透率瓶颈区,通过上述方法确定能够确定致密气储层渗透率瓶颈区,上述方法均可在实验室完成,无需工程人员到油田现场作业,这种方法简单有效,经济可靠,对指导油气生产具有重要意义。
本发明的一种致密气储层渗透率瓶颈区确定方法,本发明的其它优点、目标和特征将部分通过下面的说明体现,部分还将通过对本发明的研究和实践而为本领域的技术人员所理解。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本说明书的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1为本申请实施例提供的一种可能的致密气储层渗透率瓶颈区确定方法流程示意图;
图2为本申请实施例提供的一种可能的直井渗透率分布图;
图3为本申请实施例提供的一种可能的水平井日产量统计图;
图4为本申请实施例提供的一种可能的直井日产量统计图;
图5为本申请实施例提供的一种可能的第Ⅰ类相渗曲线示意图;
图6为本申请实施例提供的一种可能的第II类相渗曲线示意图;
图7为本申请实施例提供的一种可能的第III类相渗曲线示意图;
图8为本申请实施例提供的一种可能的第IV类相渗曲线示意图;
图9为本申请实施例提供的一种可能的第V类相渗曲线示意图;
图10为本申请实施例提供的一种可能的第VI类相渗曲线示意图;
图11为本申请实施例提供的一种可能的致密气储层渗透率瓶颈区确定装置结构示意图;
图12为本申请实施例提供的一种可能的致密气储层渗透率瓶颈区确定电子设备结构示意图。
具体实施方式
本申请实施例提供了一种致密气储层渗透率瓶颈区确定方法及相关设备,可以用砂岩微观空隙结构特征确定致密气储层的瓶颈区。
本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”、“第三”、“第四”等(如果存在)是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的实施例能够以除了在这里图示或描述的内容以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。
请参阅图1,为本申请实施例提供的一种可能的及一种致密气储层渗透率瓶颈区确定方法流程示意图,具体可以包括:
S110,获取上述致密气储层的砂岩微观孔隙结构特征;
具体的,获取上述砂岩样本,对砂岩样本进行低温液氮吸附实验、核磁共振实验、离心实验和高压压汞实验等试验,来获取上述砂岩微观孔隙结构特征包括孔隙类型、比表面积、总孔体积、特征结构参数、平均孔喉半径、孔喉大小和孔隙连通性。
S120,获取上述致密气储层的渗透率瓶颈区边界条件,其中,上述渗透率瓶颈区边界条件由数值模拟结果确定;
具体的,渗透率瓶颈区最基本的概念是在致密岩石中存在一个饱和度区域,在这个区域中水相和气相的相对渗透率都很低,以至于气相和水相都没有流动能力。利用仿真软件进行仿真,获取不同曲线数据对应的油气产量,根据上油气产量得出上述渗透率瓶颈区边界条件。
S130,基于上述砂岩微观孔隙结构特征和瓶颈区边界条件确定上述渗透率致密气储层的渗透率瓶颈区。
具体的,用砂岩的比表面积、平均孔喉半径、特征结构参数、排驱压力确定所述致密气储层的渗透率瓶颈区。
综上,上述实施例提供的致密气储层渗透率瓶颈区确定方法,通过上述目标试验获取砂岩微观孔隙结构特征,通过上述数值模拟确定渗透率瓶颈区的边界条件,经过相关性分析,最终可用砂岩的比表面积、平均孔喉半径、特征结构参数、排驱压力确定所述致密气储层的渗透率瓶颈区,这种瓶颈区确定方法简单有效,对指导油气生产具有重大意义。
在一些示例中,上述获取上述致密气储层的砂岩微观孔隙结构特征,包括:
获取上述致密气储层的砂岩样本;
对上述砂岩样本进行目标实验,得到上述砂岩微观空隙结构特征,其中上述目标实验包括低温液氮吸附实验、核磁共振实验、离心实验和高压压汞实验,上述砂岩微观孔隙结构特征包括孔隙类型、比表面积、总孔体积、特征结构参数、平均孔喉半径、孔喉大小和孔隙连通性。
具体的,通过低温液氮吸附实验、核磁共振实验、离心实验和高压压汞实验这些常规性的试验得到试验数据并将这些数据进行整理分析,得出砂岩样本的孔隙类型、比表面积、总孔体积、特征结构参数、平均孔喉半径、孔喉大小和孔隙连通性。
综上,通过上述的温液氮吸附等实验,可以测得砂岩样本的孔隙类型等结构特征,上述试验均可在实验室内完成相比于现有确定致密气储层渗透率瓶颈区在油井现场进行气水驱替实验,上述实验更为简便,实验环境更好,实验成本更低。
在一些示例中,上述获取上述致密气储层的砂岩微观孔隙结构特征,包括:
基于上述低温液氮吸附实验获取砂岩的上述孔隙类型、上述比表面积、上述总孔体积和上述特征结构参数;
基于高压压汞试验获取上述砂岩的上述平均孔喉半径和上述孔喉大小;
基于上述离心实验和上述核磁共振实验获取上述孔隙连通性。
具体地,低温液氮吸附实验是利用液氮的等温吸附特性,运用低温液氮吸附方法测定实验岩石样本孔隙的孔隙类型、比表面积、总孔体积和特征结构参数。将上述砂岩样本泡入蒸馏水24小时,并将上述砂岩样本放入离心机,使其成为束缚水状态。利用高压压汞实验获得进退汞曲线,通过曲线分析计算可以得到上述砂岩样本的孔喉特征,并对各结构特征进行定量分析,可得出岩样的平均孔喉半径和孔喉大小。核磁共振实验实验根据氢原子核本身存在的磁性同外加磁场之间的作用来进行观测,通过上述砂岩样本中流体存在的氢原子核生成的核磁共振弛豫信号变化和弛豫速率测出岩心样本孔隙连通性。
综上,岩石样本孔隙的孔隙类型、比表面积、总孔体积和特征结构参数都可以通过上述的液氮吸附实验等单一实验或两个实验的组合获取,上述实验均为常规实验,与现有确定渗透率瓶颈区的气水驱替实验比实验费用低,实验环境好。
在一些示例中,用排驱压力表示上述孔喉大小,其中,上述排驱压力是指非润湿相刚进入上述砂岩样本最大孔喉时给予的压力。
具体地,排驱压力是高压压汞实验中的一项数据,是指的非润湿相刚进入上述砂岩样本最大孔喉时给予的压力,排驱压力越大表示最大孔喉越小,孔喉的质量越差,排驱压力越小最大孔喉越大,孔喉的质量越好。
综上,可以用高压压汞试验中的排驱压力来表示孔喉的大小,孔喉大小影响孔喉的质量,孔喉的质量影响水、气在砂岩中的流动性,通过排驱压力来表述孔喉的大小,排驱压力获得比较简单,并可用其来确定渗透率瓶颈区。
请参阅图2-图4,在一些示例中,上述获取上述致密气储层的瓶颈区边界条件,包括:
获取多条致密气储层相渗透率曲线;
基于上述致密气储层相渗透率曲线进行数值模拟操作;
基于上述数值模拟的结果,确定上述瓶颈区边界条件。
具体地,根据获取的渗透率曲线数据,利用Petrel构建的骨架和属性场分布,气水相渗,流体高压物性数据、生产动态数据建立完善的数值模型。建立模型时,在气藏构造图上进行模拟网格的划分,选取矩形网格系统。纵向分为14个模拟层,网格方向、尺寸与井位相适应。共两种不同的井位图,其中一个共22口直井,井距设置为800×800;另一个模型为13口水平井,水平井水平段长度为1000m,井距设置为1000×600。将Petrel的孔隙度场、上述渗透率曲线数据导入数值模型,利用气水相渗、流体高压物性参数对模型饱和度场、压力场进行初始化,得到初始化后的数值模型,直井渗透率分布图如图2所示。六条归一化的相渗曲线对应共六个方案,开采时间为十年,单井日产量设置为6000m3。仿真后结果如图3和图4所示,从直井、水平井的产量和经济评价来判断,第III类曲线和第IV类相渗曲线之间差距较大。且IV、V、VI类曲线的产量较差。所以,IV、V、VI类曲线存在广义的渗透率瓶颈区,第IV类曲线对应的渗透率曲线为上述瓶颈区边界条件。
综上,根据获取的渗透率曲线数据,通过设置合理的参数,进行仿真模拟便可确定致密气储层的边界条件,上述方法相比现有技术简单有效,所需费用低。
在一些示例中,上述基于上述致密气储层相渗透率曲线进行数值模拟操作,包括:
基于水相在岩样中的相对渗透率、气相在岩样中的相对渗透率及水相和气相的干扰程度将上述致密气储层相渗透率曲线分组,其中,每个分组至少包含一个曲线;
将分组的曲线分别进行均值处理,得到特征曲线;
基于上述特征曲线进行数值模拟。
具体地,将气水相对渗透率曲线进行整理和分类。图5-图10共包括共108条致密砂岩储层岩样的气水相对渗透率曲线,将这些曲线进行总结和分类,并利用平均法对相渗曲线进行归一化处理。
表1相渗曲线分类
K<sub>rg</sub>(S<sub>wi</sub>) | K<sub>rw</sub>(S<sub>gr</sub>) | K<sub>x</sub>(S<sub>x</sub>) | |
I | K<sub>rg</sub>(S<sub>wi</sub>)>0.6 | K<sub>rw</sub>(S<sub>gr</sub>)>0.6 | >0.16 |
II | K<sub>rg</sub>(S<sub>wi</sub>)>0.4 | K<sub>rw</sub>(S<sub>gr</sub>)>0.6 | 0.1-0.16 |
III | K<sub>rg</sub>(S<sub>wi</sub>)>0.4 | K<sub>rw</sub>(S<sub>gr</sub>)<0.4 | <0.1 |
IV | K<sub>rg</sub>(S<sub>wi</sub>)>0.4 | K<sub>rw</sub>(S<sub>gr</sub>)>0.4 | <0.1 |
V | K<sub>rg</sub>(S<sub>wi</sub>)<0.4 | K<sub>rw</sub>(S<sub>gr</sub>)>0.6 | <0.1 |
VI | K<sub>rg</sub>(S<sub>wi</sub>)<0.3 | K<sub>rw</sub>(S<sub>gr</sub>)<0.4 | <0.05 |
表中Krg(Swi)表示水相在岩样中的相对渗透率、Krw(Sgr)表示气相在岩样中的相对渗透率、Kx(Sx)表示水相和气相的干扰程度。
为了将相渗曲线的进行量化分类,根据气水相渗曲线的Krg(Swi)、Krw(Sgr)、Kx(Sx)分类,将气-水相渗曲线分为六类,并利用平均法对收集的六类相渗曲线进行归一化处理。利用相对渗透率曲线分类归一化处理后的数据进行仿真模拟。
综上,采用大量的砂岩储层岩样的气、水相对渗透率曲线,使得测试的方法更接近于真实值,同时根据水相在岩样中的相对渗透率、气相在岩样中的相对渗透率和水相和气相的干扰程度进行分组,并将分组后的数据归一化处理后用于仿真,在不影响准确性的前提下节省了仿真资源。
在一些示例中,基于上述砂岩微观孔隙结构特征和瓶颈区边界条件确定上述致密气储层的渗透率瓶颈区,包括:
基于砂岩的上述比表面积、上述平均孔喉半径、上述特征结构参数、上述排驱压力确定上述致密气储层的渗透率瓶颈区。
具体地,根据数值模拟得到的致密气储层渗透率瓶颈区边界条件,并根据上述实验得到多个岩样的微观孔隙结构特征,为了确定哪些结构特征能够确定致密气储层,本实施例选取的是开放型楔V形孔孔隙类型的岩样,并针对岩样的比表面积、总孔体积、特征结构参数、平均孔喉半径和孔喉大小这些结构特征和由数值模拟确定的渗透率瓶颈区边界条件进行相关性分析,具体结果详见表2。
表2相关性分析
表中,Krg(Swi)表示水相在岩样中的相对渗透率,Krw(Sgr)表示气相在岩样中的相对渗透率,Kx(Sx)表示水相和气相的干扰程度,R2表示相关性,其中R2最大为1,数值越大说明相关性越高。因此根据表2的数据分析可知,比表面积同Krg(Swi)、Krw(Sgr)和Kx(Sx)关联性较好。总孔体积同所有特征点数据的相关性很低,R2都在0.2以下,说明总孔体积对气水两相流动几乎没有影响。平均孔喉半径与Krg(Swi)、Kx(Sx)的相关性好,Krg(Swi)、Kx(Sx)随着平均孔直径的减小而减小。特征结构参数则同Krg(Swi)、Krw(Sgr)和Kx(Sx)的相对性好,反映孔喉分选程度和孔喉连通程度,特征结构参数越大,气水渗流能力越好,进而Krg(Swi)、Krw(Sgr)越大。排驱压力同Krg(Swi)、Krw(Sgr)和Kx(Sx)相关性较好。排驱压力越高,说明储层孔隙结构差和连通性差,特征点参数就越小,岩样的渗流能力越差。根据相渗曲线的数值模拟结果,同对应岩样的微观孔隙结构参数相对应,对相渗曲线渗透率瓶颈区的范围进行表征。综上可以用比表面、平均孔喉半径、孔隙类型、排驱压力和特征结构参数对渗透率瓶颈区进行确定。经上述方法,确定了当孔隙类型为开放型的似片状颗粒组成的裂缝型孔,且比表面大于1.96m2/g,平均孔喉半径小于0.2um,特征结构参数小于3.815、排驱压力大于0.875MPa的岩样存在渗透率瓶颈区现象。
综上,通过上述相关性分析确定了相关性大的比表面、平均孔喉半径、特征结构参数和排驱压力用来确定渗透率瓶颈区,获取上述结构特征的试验方法简单经济,因此基于砂岩的比表面积、平均孔喉半径、特征结构参数、排驱压力确定上述致密气储层的渗透率瓶颈区方法与现有技术相比具有显著进步。
请参阅图11,本申请实施例中致密气储层渗透率瓶颈区确定装置的一个实施例,可以包括:
第一获取单元21,用于获取上述致密气储层的砂岩微观孔隙结构特征;
第二获取单元22,获取上述致密气储层的渗透率瓶颈区边界条件;
瓶颈区确定单元23,用于基于上述砂岩微观孔隙结构特征和上述渗透率瓶颈区边界条件确定上述致密气储层的渗透率瓶颈区。
综上,上述实施例提供的致密气储层渗透率瓶颈区确定装置,可对上述实验数据进行统计,对致密气储层相渗透曲线进行分类,并将分类后数据用于仿真确定致密气储层的边界条件,基于微观孔隙结构与致密气储层渗透率瓶颈区相关性分析,用所述比表面积、所述平均孔喉半径、所述特征结构参数、所述排驱压力确定所述致密气储层的渗透率瓶颈区,通过上述装置能够确定致密气储层渗透率瓶颈区,对指导油气生产具有重要意义。
请参阅图12,图12为本申请实施例提供的电子设备的实施例示意图。
如图12所示,本申请实施例还提供一种电子设备,包括存储器310、处理器320及存储在存储器320上并可在处理器上运行的计算机程序311,处理器320执行计算机程序311时实现上述致密气储层渗透率瓶颈区确定方法的任一步骤。
综上,上述实施例提供的致密气储层渗透率瓶颈区电子设备,其内部的计算机程序可对上述实验数据进行统计,对致密气储层相渗透曲线进行分类,并将分类后数据用于仿真确定致密气储层的边界条件,基于微观孔隙结构与致密气储层渗透率瓶颈区相关性分析,用所述比表面积、所述平均孔喉半径、所述特征结构参数、所述排驱压力确定所述致密气储层的渗透率瓶颈区,通过上述电子设备能够确定致密气储层渗透率瓶颈区,对指导油气生产具有重要意义。
由于本实施例所介绍的电子设备为实施本申请实施例中一种致密气储层渗透率瓶颈区确定装置所采用的设备,故而基于本申请实施例中所介绍的方法,本领域所属技术人员能够了解本实施例的电子设备的具体实施方式以及其各种变化形式,所以在此对于该电子设备如何实现本申请实施例中的方法不再详细介绍,只要本领域所属技术人员实施本申请实施例中的方法所采用的设备,都属于本申请所欲保护的范围。
在具体实施过程中,该计算机程序311被处理器执行时可以实现图1对应的实施例中任一实施方式。
需要说明的是,在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详细描述的部分,可以参见其它实施例的相关描述。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式计算机或者其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本申请实施例还提供了一种计算机程序产品,该计算机程序产品包括计算机软件指令,当计算机软件指令在处理设备上运行时,使得处理设备执行如图1对应实施例中的致密气储层渗透率瓶颈区确定中的流程。
计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行计算机程序指令时,全部或部分地产生按照本申请实施例的流程或功能。计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。计算机指令可以存储在计算机可读存储介质中,或者从一个计算机可读存储介质向另一计算机可读存储介质传输,例如,计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线(digital subscriber line,DSL))或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输。计算机可读存储介质可以是计算机能够存储的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。可用介质可以是磁性介质,(例如,软盘、硬盘、磁带)、光介质(例如,DVD)、或者半导体介质(例如固态硬盘(solid state disk,SSD))等。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的系统,装置和单元的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在本申请所提供的几个实施例中,应该理解到,所揭露的系统,装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本申请各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(Random Access Memory,RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
以上,以上实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修该,或者对其中部分技术特征进行等同替换;而这些修该或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。
Claims (8)
1.一种致密气储层渗透率瓶颈区确定方法,其特征在于,包括:
获取所述致密气储层的砂岩微观孔隙结构特征,其中,所述砂岩微观孔隙结构特征包括孔隙类型、比表面积、总孔体积、特征结构参数、平均孔喉半径、孔喉大小和孔隙连通性;
获取所述致密气储层的渗透率瓶颈区边界条件,其中,所述渗透率瓶颈区边界条件由数值模拟仿真得到的;
基于所述砂岩微观孔隙结构特征和所述渗透率瓶颈区边界条件确定所述致密气储层的渗透率瓶颈区;
其中,所述获取所述致密气储层的瓶颈区边界条件,包括:
获取多条致密气储层相渗透率曲线;
基于所述致密气储层相渗透率曲线进行数值模拟操作;
基于所述数值模拟操作的结果,确定所述瓶颈区边界条件;
其中,所述基于所述砂岩微观孔隙结构特征和瓶颈区边界条件确定所述致密气储层的渗透率瓶颈区,包括:
基于砂岩的比表面积、平均孔喉半径、特征结构参数、排驱压力确定所述致密气储层的渗透率瓶颈区。
2.如权利要求1所述的方法,其特征在于,所述获取所述致密气储层的砂岩微观孔隙结构特征,包括:
获取所述致密气储层的砂岩样本;
对所述砂岩样本进行目标实验,得到所述砂岩微观空隙结构特征,其中所述目标实验包括低温液氮吸附实验、核磁共振实验、离心实验和高压压汞实验,所述砂岩微观孔隙结构特征包括孔隙类型、比表面积、总孔体积、特征结构参数、平均孔喉半径、孔喉大小和孔隙连通性。
3.如权利要求2所述的方法,其特征在于,所述对所述砂岩样本进行目标实验,得到所述砂岩微观空隙结构特征,包括:
基于所述低温液氮吸附实验获取砂岩的所述孔隙类型、所述比表面积、所述总孔体积和所述特征结构参数;
基于高压压汞试验获取所述砂岩的所述平均孔喉半径和所述孔喉大小;
基于所述离心实验和所述核磁共振实验获取所述孔隙连通性。
4.如权利要求3所述的方法,其特征在于,包括:
用排驱压力表示所述孔喉大小,其中,所述排驱压力是指非润湿相刚进入所述砂岩样本最大孔喉时给予的压力。
5.如权利要求1所述的方法,其特征在于,所述基于所述致密气储层相渗透率曲线进行数值模拟操作,包括:
基于水相在岩样中的相对渗透率、气相在岩样中的相对渗透率及水相和气相的干扰程度将所述致密气储层相渗透率曲线分组,其中,每个分组至少包含一个曲线;
将分组的曲线分别进行均值处理,得到特征曲线;
基于所述特征曲线进行数值模拟。
6.一种致密气储层渗透率瓶颈区确定装置,其特征在于,包括:
第一获取单元,用于获取所述致密气储层的砂岩微观孔隙结构特征,其中,所述砂岩微观孔隙结构特征包括孔隙类型、比表面积、总孔体积、特征结构参数、平均孔喉半径、孔喉大小和孔隙连通性;
第二获取单元,用于获取所述致密气储层的渗透率瓶颈区边界条件;
瓶颈区确定单元,用于基于所述砂岩微观孔隙结构特征和所述渗透率瓶颈区边界条件确定所述致密气储层的渗透率瓶颈区;
其中,所述获取所述致密气储层的瓶颈区边界条件,包括:
获取多条致密气储层相渗透率曲线;
基于所述致密气储层相渗透率曲线进行数值模拟操作;
基于所述数值模拟操作的结果,确定所述瓶颈区边界条件;
其中,所述基于所述砂岩微观孔隙结构特征和瓶颈区边界条件确定所述致密气储层的渗透率瓶颈区,包括:
基于砂岩的比表面积、平均孔喉半径、特征结构参数、排驱压力确定所述致密气储层的渗透率瓶颈区。
7.一种电子设备,包括:储存器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器用于执行存储器中存储的计算机程序时实现如权利要求1-5中任一项所述的致密气储层渗透率瓶颈区确定方法的步骤。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于:所述计算机程序被处理器执行时实现如权利要求1-5中任一项所述的致密气储层渗透率瓶颈区确定方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110899572.8A CN113484222B (zh) | 2021-08-06 | 2021-08-06 | 致密气储层渗透率瓶颈区确定方法及相关设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110899572.8A CN113484222B (zh) | 2021-08-06 | 2021-08-06 | 致密气储层渗透率瓶颈区确定方法及相关设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113484222A CN113484222A (zh) | 2021-10-08 |
CN113484222B true CN113484222B (zh) | 2022-09-13 |
Family
ID=77945636
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110899572.8A Active CN113484222B (zh) | 2021-08-06 | 2021-08-06 | 致密气储层渗透率瓶颈区确定方法及相关设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113484222B (zh) |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103279647B (zh) * | 2013-05-03 | 2016-08-03 | 中国石油天然气集团公司 | 一种基于孔隙特征参数计算致密砂岩储层渗透率的方法 |
CN103335928B (zh) * | 2013-05-30 | 2015-03-04 | 中国石油天然气集团公司 | 一种测量孔隙岩石渗透率的方法和装置 |
CN104833789B (zh) * | 2015-04-09 | 2016-08-24 | 中国石油大学(北京) | 利用致密砂岩微观孔隙结构确定气水关系的方法 |
CN107133452B (zh) * | 2017-04-18 | 2019-12-03 | 中国石油大学(北京) | 油藏渗流数值模拟方法及装置 |
CN107390268B (zh) * | 2017-07-05 | 2019-06-11 | 中国石油天然气股份有限公司 | 低渗透砂砾岩成岩圈闭的识别方法及装置 |
US11774622B2 (en) * | 2017-11-06 | 2023-10-03 | Schlumberger Technology Corporation | System and method for reducing a feasible region of solutions in a relative permeability and capillary pressure curve |
CN108918396A (zh) * | 2018-09-13 | 2018-11-30 | 中国地质大学(北京) | 一种用三轴应力测试致密储层岩石渗透率的方法 |
CN109736794A (zh) * | 2019-01-24 | 2019-05-10 | 中国石油天然气股份有限公司 | 裂缝性孔隙型致密砂岩储层特征确定方法、装置及系统 |
CN111521532B (zh) * | 2019-02-01 | 2023-02-10 | 中国石油天然气股份有限公司 | 一种致密储层含气饱和度的测定方法 |
CN110206522B (zh) * | 2019-06-10 | 2019-12-20 | 西南石油大学 | 一种页岩气藏压裂水平井压裂液返排模拟方法 |
-
2021
- 2021-08-06 CN CN202110899572.8A patent/CN113484222B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113484222A (zh) | 2021-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106525688B (zh) | 一种饱和页岩孔隙流体分离及饱和度计算的实验方法 | |
CN105447762B (zh) | 一种流体替换的低渗透油藏水淹信息的计算方法 | |
CN110927035A (zh) | 一种低渗致密砂岩束缚水饱和度计算方法 | |
CN109030311A (zh) | 基于核磁共振t2谱敏感参数的孔隙结构分类与识别方法 | |
CN106442600A (zh) | 页岩束缚水含量确定方法 | |
CN106840993B (zh) | 一种确定岩心可动流体空间的方法及装置 | |
CN111521539B (zh) | 一种致密砂砾岩储层中不同类型孔隙结构定量表征方法 | |
CN104374685B (zh) | 一种数字岩心三孔隙组分求取储层渗透率的方法 | |
CN105240006B (zh) | 一种适用于火山岩储层的油水层识别方法 | |
CN107132171A (zh) | 一种基于压汞‑氮吸附联测数据确定致密储层孔径分布的方法 | |
CN103778328A (zh) | 基于标准偏差分析的储层敏感孔喉提取方法 | |
CN114370269B (zh) | 深层碳酸盐岩气藏有效储层物性下限综合确定方法 | |
CN113484222B (zh) | 致密气储层渗透率瓶颈区确定方法及相关设备 | |
CN111507986A (zh) | 一种孔隙内波及情况的确定方法和装置 | |
CN110222981B (zh) | 一种基于参数二次选取的储层分类评价方法 | |
CN116821786A (zh) | 一种低渗储层分类方法、装置及介质 | |
Li et al. | Pore size distribution of a tight sandstone reservoir and its effect on micro pore‐throat structure: a case study of the Chang 7 member of the Xin'anbian Block, Ordos Basin, China | |
CN112049623A (zh) | 一种中基性火山岩束缚水出水的判定方法和装置 | |
CN110593864B (zh) | 一种利用核磁-声波流体识别的方法和装置 | |
CN110656924A (zh) | 一种超低渗透油藏分类方法 | |
CN113376069A (zh) | 致密储层物性下限的确定方法及装置 | |
CN114577693B (zh) | 泥页岩储层确定方法、装置、计算机设备及存储介质 | |
CN114062190B (zh) | 表征页岩油微观赋存特征的方法、装置、终端及存储介质 | |
CN113642656B (zh) | 一种低渗砂岩储层开采方式确定方法及相关装置 | |
CN113642157B (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 |