CN110702579B - 基于渗流实验的深层非均质滩相气藏流动单元划分方法 - Google Patents
基于渗流实验的深层非均质滩相气藏流动单元划分方法 Download PDFInfo
- Publication number
- CN110702579B CN110702579B CN201911000486.8A CN201911000486A CN110702579B CN 110702579 B CN110702579 B CN 110702579B CN 201911000486 A CN201911000486 A CN 201911000486A CN 110702579 B CN110702579 B CN 110702579B
- Authority
- CN
- China
- Prior art keywords
- flow
- reservoir
- seepage
- sample
- parameters
- 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
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000002474 experimental method Methods 0.000 title claims abstract description 44
- 230000035699 permeability Effects 0.000 claims abstract description 45
- 239000011148 porous material Substances 0.000 claims abstract description 39
- 238000000638 solvent extraction Methods 0.000 claims abstract description 31
- 238000005325 percolation Methods 0.000 claims abstract description 20
- 238000004458 analytical method Methods 0.000 claims abstract description 14
- 238000012360 testing method Methods 0.000 claims description 17
- 238000005192 partition Methods 0.000 claims description 10
- 239000012530 fluid Substances 0.000 claims description 7
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 6
- 229910052753 mercury Inorganic materials 0.000 claims description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 6
- 238000005266 casting Methods 0.000 claims description 4
- 238000007621 cluster analysis Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000001035 drying Methods 0.000 claims description 3
- 238000002360 preparation method Methods 0.000 claims description 3
- 238000010219 correlation analysis Methods 0.000 claims description 2
- 238000011068 loading method Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 238000003860 storage Methods 0.000 claims description 2
- 230000010287 polarization Effects 0.000 claims 1
- 238000004062 sedimentation Methods 0.000 claims 1
- 238000011161 development Methods 0.000 abstract description 5
- 239000013078 crystal Substances 0.000 description 16
- 238000004090 dissolution Methods 0.000 description 10
- 229910000514 dolomite Inorganic materials 0.000 description 8
- 239000010459 dolomite Substances 0.000 description 8
- 230000015572 biosynthetic process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 239000011435 rock Substances 0.000 description 6
- 208000035126 Facies Diseases 0.000 description 4
- 230000008901 benefit Effects 0.000 description 4
- 230000032823 cell division Effects 0.000 description 4
- 230000018109 developmental process Effects 0.000 description 4
- 239000002245 particle Substances 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000003628 erosive effect Effects 0.000 description 2
- 238000003703 image analysis method Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000000227 grinding Methods 0.000 description 1
- 238000004573 interface analysis Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000001308 synthesis method Methods 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
- G01N15/082—Investigating permeability by forcing a fluid through a sample
- G01N15/0826—Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
-
- 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
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Dispersion Chemistry (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Solid-Sorbent Or Filter-Aiding Compositions (AREA)
Abstract
本发明属于油气开发技术领域，公开了一种基于渗流实验的深层非均质滩相气藏流动单元划分方法。进行样品采集与分析：选取不同储层类型的样品；开展渗流特征实验：明确不同样品的渗流特征；分析样品孔隙度、渗透率、R35与渗流指数的相关性，优选流动单元划分关键参数；在渗流指数与渗透率、R35的交汇图中，将具有相近FI值的样品作为同一流动单元，并读取对应的渗透率及R35值，完成流动单元分类与阀门值的确定。本发明揭示了流动单元渗流特征差异的内涵，更加准确的区分流动单元；提供了一种划分流动单元的新方法，能够对强非均质孔隙型滩相气藏流动单元划分进行有效的支撑。
Description
技术领域
本发明属于油气开发领域中地质学研究技术领域，尤其涉及一种基于渗流实验的深层非均质滩相气藏流动单元划分方法。
背景技术
目前，最接近的现有技术：
我国海相碳酸盐岩滩相气藏分布广泛，多为孔隙型储层，具有孔隙类型多样、孔渗关系及孔隙结构复杂、高孔渗与低孔渗并存及非均质性强的特点，如普光气田飞仙关组、高磨气田龙王庙组。在开发过程中多采用长井段射孔的方式采气，进入中后期，气藏低孔低渗储层动用程度差的矛盾日益突出，影响气藏整体稳产。针对这一现状，开展流动单元划分，可以为细分气藏开发层系及合理动用剩余储量提供有价值的支撑。
流动单元的概念综合起来可以表述为:具有特定(相似)渗流特征，一定空间分布的、可以预测的储集体。流动单元划分经历多年发展，归纳起来主要包括四类方法：1)沉积学研究方法，如沉积界面分析法及沉积岩相约束法；2)基于多参数的储层物性研究方法，如旋回约束下多参数物性综合法、储层自然分段多参数聚类分析法；3)微观孔隙结构-岩石物理相研究方法，如孔喉半径法(R35)、岩石物理相法；4)基于储层层带指标(FZI)的研究方法。上述研究方法能静态反映储层的储集空间及孔喉结构的差异，能基本支撑油藏流动单元划分。但这些方法应用于气藏流动单元划分，尤其是孔渗关系复杂的非均质气藏，还存在不足：因为气体在储层中的渗流特征受温度压力影响较大，基于单纯的静态储层参数划分流动单元，不能真实反映气体在气藏中的渗流特征。
因此，形成一种渗流特征实验约束下的非均质滩相储层流动单元划分方法，它能真实反映气体在不同储层中渗流特征的差异，能揭示气藏流动单元划分的本质，也是对气藏流动单元划分的一种新的尝试。
综上所述，现有技术存在的问题是：
现有方法应用于气藏流动单元划分，尤其是孔渗关系复杂的非均质气藏，还存在不足：因为气体在储层中的渗流特征受温度压力影响较大，基于单纯的静态储层参数划分流动单元，不能真实反映气体在气藏中的渗流特征。
而且现有技术深层非均质滩相气藏流动单元划分中，区分流动单元的数据准确性差；不能够全面反映一般规律信息。
发明内容
针对现有技术存在的问题，本发明提供了一种基于渗流实验的深层非均质滩相气藏流动单元划分方法，能够对强非均质孔隙型气藏流动单元划分进行有效的支撑。
本发明是这样实现的，一种基于渗流实验的深层非均质滩相气藏流动单元划分方法，所述方法包括以下步骤：
步骤一，样品采集与分析：基于岩心观察及储层样品挑选，测试样品孔隙类型、孔隙度及渗透率，明确储层岩性及储集空间类型。
步骤二，开展高温高压渗流特征实验：选取不同储层类型的样品，开展渗流特征实验，明确不同样品的渗流特征。
步骤三，流动单元划分参数优选：分析样品孔隙度、渗透率、R35与渗流指数(FI)的相关性，优选流动单元划分关键参数。
步骤四，流动单元划分关及其关键参数阈值区间确定：在渗流指数(FI)与渗透率、R35的交汇图中，将具有相近FI值的样品作为同一流动单元，并读取对应的渗透率及R35值，完成流动单元分类与阀门值的确定。
进一步，步骤一中，所述样品采集与分析具体包括：
1)岩心观察与分析，采集具有不同岩性(相)、储集性的储层样品。
2)将样品备制成直径为2.5cm的柱塞样，完成孔隙度及渗透率测试。残样备制成铸体薄片，利用普通透射光偏关显微镜分析储层孔隙类型。
进一步，步骤二中，所述开展高温高压渗流特征实验具体包括：
采用高温高压恒温箱；设置实验温度：80℃，实验流体：99.99％甲烷，覆压(地层压力)：30Mpa。实验过程如下：
1)在步骤一基础上，优选能不同岩性/储层类型的柱塞样品，保证每类样品至少3个。
2)样品烘干，测量样品长度、直径及质量。
3)将样品装入夹持器，并提供一个稳定的覆(围)压P4；同时增加进气端压力P2和出气段回压P3至3Mpa，确保进气端压力与回压平衡。
4)4ml/min恒速增加进气段压力P2，观察末端流量计的变化，并记录时间(T,min)、P2、P3及流量(Q，ml/min)，直至流量稳定。
5)改变样品，重复上述步骤。
进一步，所述流动单元划分参数一般包括沉积特征参数(微相、岩性、粒度中值、储层厚度、泥质含量及净毛比等)、储层宏观参数(孔隙度、渗透率)、储层微观孔隙结构参数(孔喉半径)、流体物理参数(饱和度、粘度、密度)及综合参数(储层系数、流动层带指数)。
进一步，步骤三中，所述流动单元划分参数优选具体包括：
1)同一岩性有不同的孔隙类型、孔隙度及渗透率，因此沉积特征参数不能用于流动单元划分；选择能代表储层宏观及微观特征的孔隙度、渗透率及R35参数，结合FI参数，作为流动单元划分的基本参数。
基于渗流实验，获取稳定渗流状态下的流量(Q，单位ml/min)与压差(ΔP，单位Mpa)，提出渗流指数FI＝Q/ΔP，单位ml/(min.Pa)。
R35为进汞饱和度为35％时对应的喉道半径，单位为μm。该参数能直观反映孔喉半径及孔喉结构质量，它是影响渗透率的主要参数。R35的值可以通过压汞法或图像分析法直接通过测试数据获取，也可利用Winland(1)公式求取：
lgR35＝0.372+0.588lgK-0.684lgΦ (1)
2)整理样品的FI、渗透率、孔隙度及R35，开展交汇分析。提出与FI相关性强的参数以及FI参数作为流动单元划分的关键参数。
进一步，步骤四中，所述流动单元划分关及其关键参数阈值区间确定具体包括：
1)流动单元划分。在交汇图相关性分析基础上，分析交汇点的分异度。
2)流动单元关键参数阀值区间确定。
分为两个层次，第一个层次为基于样品测试结果的阀值区间确定，受样品数量限制，样品数越多，阀值区间越接近事实。
第二个层次，基于测井解释结果的阀值区间确定。在样品测试数据的基础上，利用测井解释结果，完成孔隙度、渗透率、R35及FI的数据准备(FI值，可以通过实验结果中FI与R35的关系式求取)。然后，在Sufer软件或SPSS软件中，用4类流动单元的认识约束，开展聚类分析，实现4类流动单元关键参数阀值区间的确定。
本发明的另一目的在于提供一种实施所述基于渗流实验的深层非均质滩相气藏流动单元划分方法的基于渗流实验的深层非均质滩相气藏流动单元划分系统。
本发明的另一目的在于提供一种实施所述基于渗流实验的深层非均质滩相气藏流动单元划分方法的信息数据处理终端。
本发明的另一目的在于提供一种计算机可读存储介质，包括指令，当其在计算机上运行时，使得计算机执行所述的基于渗流实验的深层非均质滩相气藏流动单元划分方法。
综上所述，本发明的优点及积极效果为：
本发明提供的基于渗流实验的深层非均质滩相孔隙型气藏流动单元划分方法，揭示流动单元渗流特征差异的内涵，从而更加准确的区分流动单元；提供了一种划分流动单元的新方法，能够对强非均质孔隙型储层流动单元划分进行有效的支撑。本发明在Sufer软件或SPSS软件中，将少样品点的分类拓展到多样品点的聚类分析，优点在于能够较为全面反映一般规律。
附图说明
图1是本发明实施例提供的基于渗流实验的深层非均质滩相气藏流动单元划分方法流程图。
图2是本发明实施例提供的基于渗流实验的深层非均质滩相气藏流动单元划分方法原理图。
图3是本发明实施例提供的实验装置示意图。
图4是本发明实施例提供的飞仙关组储层岩性及孔隙类型示意图。
图中：P1是残余鲕粒云岩，颗粒间溶蚀孔(G-Inter-DP)，蓝色为铸体；P2是残余鲕粒云岩，颗粒间溶蚀孔(G-Inter-DP)及颗粒内溶蚀孔(G-Intra-DP；P3是残余鲕粒云岩，鲕粒内溶蚀孔(G-Intra-DP)；P4是中晶白云岩，晶间溶孔(C-Inter-DP)；P5是中-细晶白云岩，晶间孔(C-Inter-P)及晶间溶孔(C-Inter-DP)；P6是细晶白云岩，晶间孔(C-Inter-P)。
图5是本发明实施例提供的渗透率与FI的关系示意图。
图6是本发明实施例提供的孔隙度与FI的关系示意图。
图7是本发明实施例提供的R35与FI的关系示意图。
图8是本发明实施例提供的基于FI与渗透率交汇图的流动单元划分图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
针对现有技术存在的问题，本发明提供了一种基于渗流实验的深层非均质滩相气藏流动单元划分方法，下面结合附图对本发明作详细的描述。
如图1所示，本发明实施例提供的一种基于渗流实验的深层非均质滩相气藏流动单元划分方法包括以下步骤：
S101，样品采集与分析：基于岩心观察及储层样品挑选，测试样品孔隙类型、孔隙度及渗透率，明确储层岩性及储集空间类型。
S102，开展高温高压渗流特征实验：选取不同储层类型的样品，开展渗流特征实验，明确不同样品的渗流特征。
S103，流动单元划分参数优选：分析样品孔隙度、渗透率、R35与渗流指数(FI)的相关性，优选流动单元划分关键参数。
S104，流动单元划分关及其关键参数阈值区间确定：在渗流指数(FI)与渗透率、R35的交汇图中，将具有相近FI值的样品作为同一流动单元，并读取对应的渗透率及R35值，完成流动单元分类与阀门值的确定。
在本发明实施例中，如图2所示，为本发明实施例提供的基于渗流实验的深层非均质滩相气藏流动单元划分方法原理。
下面结合实施例对本发明作进一步描述。
实施例
普光气田飞仙关组主要发育鲕粒滩相孔隙型储层，具有孔渗关系及孔隙结构复杂、高孔渗与低孔渗并存及非均质性强的特点。对此，本次实施提供了针对复杂非均质滩相气层流动单元划分方法，也是上述方法的进一步阐述。包括以下步骤。
步骤一：样品储层特征分析。
1).岩心观察与分析，采集具有不同岩性(相)、储集性的储层样品。其中，岩心样品的观察方法参考标准《SY/T 5336-1996岩心常规分析方法》。
飞仙关组储层主要为颗粒滩相，岩性为白云岩、残余鲕粒白云岩(如图4所示)。
2).将样品备制成直径为2.5cm的柱塞样，完成孔隙度及渗透率测试。残样备制成铸体薄片，利用普通透射光偏关显微镜分析储层孔隙类型。其中，孔隙度及渗透率测量，参考标准《SY/T 6385-2016覆压下岩石孔隙度及渗透率测定方法》；铸体薄片制作参考标准《SY/T 5913-2004岩石制片方法》。
飞仙关组孔隙类型主要为：粒间溶孔、粒内溶孔及粒间-粒内复合孔(如图4所示)。
步骤二：高温高压渗流特征实验结果。
我国滩相气藏一般具有埋藏深、地层温度高、地层压力大的特点。以普光气田为例，目前气藏的压力为30-50Mpa左右，地层温度1300c左右。为了尽可能模地层条件，客观真实反映渗流特征。本次设计采用高温高压恒温箱，实验设备如图3所示。设置实验温度：80℃；实验流体：99.99％甲烷；覆压(地层压力)：30Mpa。实验过程如下：
1).在步骤一基础上，优选能不同岩性/储层类型的柱塞样品，保证每类样品至少3个；
2).样品烘干，测量样品长度、直径及质量；
3).将样品装入夹持器，并提供一个稳定的覆(围)压P4(压力表4测定)；同时增加进气端压力P2(压力表P2测定)和出气段回压P3(压力表P3测定)至3Mpa，确保进气端压力与回压平衡。
4).恒速(4ml/min)增加进气段压力P2，观察末端流量计的变化，并记录时间(T,min)、P2、P3及流量(Q，ml/min)，直至流量稳定。
5).改变样品，重复上述步骤。
优选14个柱塞样品，开展高位高压渗流特征实验，记录实验结果(如表1所示)。
表1渗流特征实验结果及样品储层参数表
在本发明实施例中，图4是本发明实施例提供的飞仙关组储层岩性及孔隙类型；图中：P1是残余鲕粒云岩，颗粒间溶蚀孔(G-Inter-DP)，蓝色为铸体；P2是残余鲕粒云岩，颗粒间溶蚀孔(G-Inter-DP)及颗粒内溶蚀孔(G-Intra-DP；P3是残余鲕粒云岩，鲕粒内溶蚀孔(G-Intra-DP)；P4是中晶白云岩，晶间溶孔(C-Inter-DP)；P5是中-细晶白云岩，晶间孔(C-Inter-P)及晶间溶孔(C-Inter-DP)；P6是细晶白云岩，晶间孔(C-Inter-P)。
步骤三：流动单元划分参数优选。
流动单元划分参数一般包括沉积特征参数(如，微相、岩性、粒度中值、储层厚度、泥质含量及净毛比等)、储层宏观参数(如，孔隙度、渗透率)、储层微观孔隙结构参数(如孔喉半径)、流体物理参数(如，饱和度、粘度、密度)及综合参数(如，储层系数、流动层带指数)。
1).通过储层基本特征分析发现，同一岩性有不同的孔隙类型、孔隙度及渗透率，因此沉积特征参数不能用于流动单元划分。因为气藏的特殊性，流体物理参数及流动系数也不能支撑流动单元划分。基于此，选择能代表储层宏观及微观特征的孔隙度、渗透率及R35参数，结合FI参数，作为流动单元划分的基本参数。
基于渗流实验，获取稳定出气状态下的流量(Q，单位ml/min)与压差(ΔP，单位Mpa)，提出渗流指数FI＝Q/ΔP，单位ml/(min.Pa)。该指数能准确表达地层条件下气体通过储层的能力。
R35，为进汞饱和度为35％时对应的喉道半径，单位为μm。该参数能直观反映孔喉半径及孔喉结构质量，它是影响渗透率的主要参数。R35的值可以通过压汞法或图像分析法直接通过测试数据获取，也可利用Winland(1)公式求取：
lgR35＝0.372+0.588lgK-0.684lgΦ (1)
2).整理样品的FI、渗透率、孔隙度及R35，开展交汇分析。提出与FI相关性强的参数以及FI参数作为流动单元划分的关键参数。从图5、6、7可以看出，渗透率、R35两个参数与FI的拟合关系良好，说明渗透率及R35能支撑流动单元划分。孔隙度与FI拟合关系差，说明孔隙度不能支撑流动单元划分。因此，渗透率、R35及FI三个参数可以作为流动单元划分的关键参数。
步骤四：流动单元划分关及其关键参数阈值区间确定。
1).流动单元划分。从FI与渗透率交汇图(如图8所示)看出，明显存在5个分布区间，说明存在5种不同的渗流特征及对应的流动单元。
2).流动单元关键参数阀值区间确定。
基于样品测试结果的阀值区间确定。在FI与渗透率交汇图中(如图8所示)，通每类流动单元的中值线客运代表对应流动单元的阀值，进而确定5类流动单元FI、渗透率及R35的阀值区间(如表2所示)。该结果受样品数量影响，样品数越多，阀值区间越接近事实。
表2基于实验样品的流动单元分类及其阀值区间
基于测井解释结果的阀值区间确定。在样品测试数据的基础上，利用测井解释结果，完成批量孔隙度、渗透率、R35及FI的数据准备。FI可以通过图8获取的关系式求取。在SPSS软件中，用5类流动单元的认识约束，开展聚类分析，获取每类流动单元的阀值区间(如表3所示)。不同的流动单元的阀值区间值存在一定的重叠部分，而两种方法获取的5类流动单元的阀值区间基本上一致。该方法的优点在于样本点多，能够较为全面反映一般规律，缺点在于测井物性解释结果与实际结果存在误差。
表3基于测井解释样品的的流动单元分类及其阀值区间
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (9)
1.一种基于渗流实验的深层非均质滩相气藏流动单元划分方法，其特征在于，所述深层非均质滩相气藏流动单元划分方法包括以下步骤：
步骤一，样品采集与分析：测试挑选的储层样品孔隙类型、孔隙度及渗透率，分析储层岩性及储集空间类型；
步骤二，高温高压渗流特征实验：选取不同储层类型的样品，进行渗流特征实验，分析不同样品的渗流特征；
步骤三，流动单元划分参数优选：分析样品孔隙度、渗透率、R35与渗流指数FI的相关性，优选流动单元划分关键参数；
步骤三中，所述流动单元划分参数优选具体包括：
1）选择能代表储层宏观及微观特征的孔隙度、渗透率及R35参数，结合FI参数，作为流动单元划分的基本参数；
基于渗流实验，获取稳定渗流状态下的流量与压差，渗流指数FI= Q/∆P，单位ml/(min.Pa)；
R35为进汞饱和度为35%时对应的喉道半径，喉道半径是影响储层渗流能力的关键参数，单位为μm；R35的值通过压汞实验求取，或通过经验公式求取：
lgR35=0.732+0.588lgK-0.864lgΦ；
2）整理样品的FI、渗透率、孔隙度及R35，交汇分析；提出与FI相关性强的参数以及FI参数作为流动单元划分的关键参数；
步骤四，流动单元划分及关键参数阈值区间确定：在渗流指数FI与渗透率、R35的交汇图中，将具有相近FI值的样品作为同一流动单元，并读取对应的渗透率及R35值，完成流动单元分类与阀门值的确定。
2.如权利要求1所述的基于渗流实验的深层非均质滩相气藏流动单元划分方法，其特征在于，步骤一中，所述样品采集与分析具体包括：
1）岩心分析，采集具有不同岩性、储集性的储层样品；
2）将样品备制成直径为2.5cm的柱塞样，完成孔隙度及渗透率测试；残样备制成铸体薄片，利用普通透射光偏光显微镜分析储层孔隙类型。
3.如权利要求1所述的基于渗流实验的深层非均质滩相气藏流动单元划分方法，其特征在于，步骤二中，所述开展高温高压渗流特征实验中采用高温高压恒温箱；设置实验温度：80℃，实验流体：99.99%甲烷，覆压：30Mpa；
具体包括：
1）在步骤一基础上，优选不同岩性/储层类型的柱塞样品，每类样品至少3个；
2）样品烘干，测量样品长度、直径及质量；
3）将样品装入夹持器，并提供一个稳定的覆压P4；同时增加进气端压力P2和出气段回压P3至3Mpa，使进气端压力与回压平衡；
4）4ml/min恒速增加进气段压力P2，分析末端流量计的变化，并记录时间、P2、P3及流量，直至流量稳定；
5）改变样品，重复上述步骤1）-步骤4）。
4.如权利要求1所述的基于渗流实验的深层非均质滩相气藏流动单元划分方法，其特征在于，所述流动单元划分参数包括：
沉积特征参数、储层宏观参数、储层微观孔隙结构参数、流体物理参数及综合参数。
5.如权利要求1所述的基于渗流实验的深层非均质滩相气藏流动单元划分方法，其特征在于，步骤四中，所述流动单元划分及其关键参数阈值区间确定具体包括：
1）流动单元划分；在交汇图相关性分析基础上，分析交汇点的分异度；
2）分为两个层次确定流动单元关键参数阀值区间。
6.如权利要求5所述的基于渗流实验的深层非均质滩相气藏流动单元划分方法，其特征在于，分为两个层次确定流动单元关键参数阀值区间具体包括：
第一个层次为基于样品测试结果的阀值区间确定，受样品数量限制，样品数越多，阀值区间越接近事实；
第二个层次，基于测井解释结果的阀值区间确定；在样品测试数据的基础上，利用测井解释结果，完成孔隙度、渗透率、R35及FI的数据准备；在Sufer软件或SPSS软件中，用4类流动单元的认识约束，聚类分析，进行 4类流动单元关键参数阀值区间的确定。
7.一种实施权利要求1所述基于渗流实验的深层非均质滩相气藏流动单元划分方法的基于渗流实验的深层非均质滩相气藏流动单元划分系统。
8.一种实施权利要求1~6任意一项所述基于渗流实验的深层非均质滩相气藏流动单元划分方法的信息数据处理终端。
9.一种计算机可读存储介质，包括指令，当其在计算机上运行时，使得计算机执行如权利要求1-6任意一项所述的基于渗流实验的深层非均质滩相气藏流动单元划分方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911000486.8A CN110702579B (zh) | 2019-10-21 | 2019-10-21 | 基于渗流实验的深层非均质滩相气藏流动单元划分方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911000486.8A CN110702579B (zh) | 2019-10-21 | 2019-10-21 | 基于渗流实验的深层非均质滩相气藏流动单元划分方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110702579A CN110702579A (zh) | 2020-01-17 |
| CN110702579B true CN110702579B (zh) | 2022-05-20 |
Family
ID=69201963
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911000486.8A Active CN110702579B (zh) | 2019-10-21 | 2019-10-21 | 基于渗流实验的深层非均质滩相气藏流动单元划分方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110702579B (zh) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103902827A (zh) * | 2014-04-03 | 2014-07-02 | 刘鑫 | 一种碳酸盐岩水平井流动单元划分方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104153769B (zh) * | 2014-07-04 | 2017-01-04 | 中国石油大学(北京) | 一种缝洞型油藏流动单元的划分及评价方法 |
| CN106372446B (zh) * | 2016-10-12 | 2018-12-18 | 中国地质大学（北京） | 流动单元的划分方法和装置及井间对比方法和装置 |
| CN109958413B (zh) * | 2019-03-01 | 2021-03-30 | 中国石油大学(华东) | 一种特高含水期油藏动态流动单元划分方法 |
-
2019
- 2019-10-21 CN CN201911000486.8A patent/CN110702579B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103902827A (zh) * | 2014-04-03 | 2014-07-02 | 刘鑫 | 一种碳酸盐岩水平井流动单元划分方法 |
Non-Patent Citations (2)
| Title |
|---|
| 基于流动单元的致密砂岩气储层渗透率测井评价――以川中广安地区须家河组为例;陈志强等;《岩性油气藏》;20171201(第06期);79-86页 * |
| 榆林气田山_2~3亚段低渗砂岩气藏流动单元精细划分;杨国平等;《天然气勘探与开发》;20170325(第01期);48-55页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110702579A (zh) | 2020-01-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wenlian et al. | Evaluation of the stress sensitivity in tight reservoirs | |
| CN109030311B (zh) | 基于核磁共振t2谱敏感参数的孔隙结构分类与识别方法 | |
| Nabawy et al. | Pore-throat characterization in highly porous and permeable sandstones | |
| Cao et al. | Factors influencing microstructure and porosity in shales of the Wufeng-Longmaxi formations in northwestern Guizhou, China | |
| Aliakbardoust et al. | Effects of pore geometry and rock properties on water saturation of a carbonate reservoir | |
| Cheng et al. | Nuclear magnetic resonance features of low-permeability reservoirs with complex wettability | |
| Hosseini et al. | The effect of heterogeneity on NMR derived capillary pressure curves, case study of Dariyan tight carbonate reservoir in the central Persian Gulf | |
| CN110470584B (zh) | 一种评价渗吸和水锁综合效应的方法 | |
| Li et al. | A novel model for assessing the pore structure of tight sands and its application | |
| Pan et al. | Study on the pore structure, fluid mobility, and oiliness of the lacustrine organic-rich shale affected by volcanic ash from the Permian Lucaogou Formation in the Santanghu Basin, Northwest China | |
| Sneider et al. | Rock types, depositional history, and diagenetic effects, Ivishak reservoir, Prudhoe Bay field | |
| Landry et al. | Estimating mudrock oil-water relative permeability curves using digital rock physics | |
| Li et al. | Pore‐Throat Combination Types and Gas‐Water Relative Permeability Responses of Tight Gas Sandstone Reservoirs in the Zizhou Area of East Ordos Basin, China | |
| CN105464655B (zh) | 一种流体测井识别方法 | |
| Duan et al. | Characterization and classification of pore structures in deeply buried carbonate rocks based on mono-and multifractal methods | |
| CN110702579B (zh) | 基于渗流实验的深层非均质滩相气藏流动单元划分方法 | |
| CN113137225A (zh) | 一种基于砾石组分分析的砾岩储层评价方法 | |
| CN108802192B (zh) | 一种钙屑砂岩储层孔隙类型识别方法 | |
| Sneider et al. | Integrated rock-log calibration in the Elmworth field-Alberta, Canada: Reservoir rock detection and characterization: Part i | |
| CN112557277B (zh) | 一种致密岩石微孔隙连通性识别分类评价方法 | |
| CN115905917A (zh) | 综合静动态数据构建海域低渗气层分类曲线的方法 | |
| Li et al. | Application of nuclear magnetic resonance logging in the low-resistivity reservoir—Taking the XP area as an example | |
| Gao et al. | Influence of different pore structure types on the occurrence features of movable fluid in Chang-10 reservoir in Wuqi–Ansai Oilfield | |
| Wang et al. | Characterization of Movable Fluid Distribution in Tight Oil Reservoirs by NMR and Centrifugation | |
| Rahman et al. | Effects of microporosity on permeability and sonic velocity of Miocene carbonates and an approach to relate micrite microtextures with microporosity occurrences in Miocene carbonate reservoirs of offshore Sarawak, Malaysia |
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 |