CN113533156B - 页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法 - Google Patents
页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法 Download PDFInfo
- Publication number
- CN113533156B CN113533156B CN202110740441.5A CN202110740441A CN113533156B CN 113533156 B CN113533156 B CN 113533156B CN 202110740441 A CN202110740441 A CN 202110740441A CN 113533156 B CN113533156 B CN 113533156B
- Authority
- CN
- China
- Prior art keywords
- core
- shale
- water
- pore
- magnetic resonance
- 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
- 239000011148 porous material Substances 0.000 title claims abstract description 144
- 239000012530 fluid Substances 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 25
- 239000003079 shale oil Substances 0.000 title claims abstract description 24
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 112
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 46
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 41
- 238000001228 spectrum Methods 0.000 claims abstract description 36
- 238000009826 distribution Methods 0.000 claims abstract description 31
- 238000010438 heat treatment Methods 0.000 claims abstract description 24
- 230000003595 spectral effect Effects 0.000 claims abstract description 16
- 230000001186 cumulative effect Effects 0.000 claims description 37
- 238000005119 centrifugation Methods 0.000 claims description 20
- LYCAIKOWRPUZTN-UHFFFAOYSA-N Ethylene glycol Chemical compound OCCO LYCAIKOWRPUZTN-UHFFFAOYSA-N 0.000 claims description 12
- 238000001035 drying Methods 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 8
- 230000000704 physical effect Effects 0.000 claims description 8
- 239000000203 mixture Substances 0.000 claims description 7
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 6
- 239000011707 mineral Substances 0.000 claims description 6
- 238000000926 separation method Methods 0.000 claims description 5
- CAHQGWAXKLQREW-UHFFFAOYSA-N Benzal chloride Chemical compound ClC(Cl)C1=CC=CC=C1 CAHQGWAXKLQREW-UHFFFAOYSA-N 0.000 claims description 4
- 239000002904 solvent Substances 0.000 claims description 4
- 238000004140 cleaning Methods 0.000 claims description 3
- 238000003825 pressing Methods 0.000 claims description 2
- 238000012512 characterization method Methods 0.000 claims 5
- 239000011435 rock Substances 0.000 abstract description 7
- 230000008569 process Effects 0.000 abstract description 5
- 238000013508 migration Methods 0.000 abstract description 4
- 230000005012 migration Effects 0.000 abstract description 4
- 239000008398 formation water Substances 0.000 abstract 1
- 238000007781 pre-processing Methods 0.000 abstract 1
- 238000011161 development Methods 0.000 description 8
- 238000002474 experimental method Methods 0.000 description 8
- 239000002734 clay mineral Substances 0.000 description 7
- 230000009471 action Effects 0.000 description 6
- 239000004927 clay Substances 0.000 description 6
- 230000035699 permeability Effects 0.000 description 6
- 230000007423 decrease Effects 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 3
- 239000013078 crystal Substances 0.000 description 3
- 239000007789 gas Substances 0.000 description 3
- 230000036571 hydration Effects 0.000 description 3
- 238000006703 hydration reaction Methods 0.000 description 3
- 229910052739 hydrogen Inorganic materials 0.000 description 3
- 239000001257 hydrogen Substances 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000003628 erosive effect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 150000001768 cations Chemical class 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000013467 fragmentation Methods 0.000 description 1
- 238000006062 fragmentation reaction Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000004094 surface-active agent Substances 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
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
Landscapes
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Dispersion Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
Abstract
本发明公开了一种页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法，包括：1)对岩心预处理；2)高压注入KCl溶液，再常压注入地层水，得到完全饱和水状态下的核磁共振T2谱；3)根据完全饱和水状态下的核磁共振T2谱分布的形态，将孔隙结构划分为三种类型，4)依次在不同转速条件下离心，并测量不同离心状态下岩心的核磁共振T2谱，确定最佳离心力；5)重复步骤2)，将最佳离心力离心后的岩心进行热处理，分别测量不同温度后岩心的核磁共振T2谱；6)根据上述的T2谱分布，计算得到可动水与束缚水的T2截止值T2C1，毛管束缚水与不可采出水的T2截止值T2C2。采用本发明方法可以明确页岩油储集层的孔隙结构特征及其孔隙流体赋存和运移规律。
Description
技术领域
本发明属于油气开发领域，涉及一种页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法。
背景技术
我国页岩油气资源丰富，具有巨大的开发潜力，但页岩油储集层呈现低孔、特低渗、孔隙结构差、非均质性严重等特点，给大规模商业开采带来了极大困难，如何优选出高品质储集层是决定页岩油能否增产稳产的关键。
页岩复杂的矿物组成和孔隙结构特征导致孔隙流体类型也千差万别，准确划分不同类型孔隙流体是优选高品质页岩油储集层的关键。通常，页岩孔隙流体类型可划分为可动流体和束缚流体，其中束缚流体主要包括毛管束缚流体和黏土束缚流体。毛管束缚流体主要指细孔微喉中赋存流体在巨大毛管压力作用下而无法参与流动的流体，但当加入表面活性剂或升温后界面张力降低引起毛管压力降低，这部分束缚流体又能重新参与流动。黏土束缚流体是指黏土晶层在面与面堆积形成晶体过程中，带有负电荷的黏土晶层吸附水合阳离子并形成水化膜，这部分流体又称为不可采出流体。
为识别核磁共振T2谱分布中不同类型孔隙流体，定义T2C1和T2C2分别为可动流体与束缚流体的弛豫时间截止值及毛管束缚流体与黏土束缚流体的弛豫时间截止值。目前本领域技术人员在研究页岩核磁共振特性的时候都是采用常规砂岩的评价方法，然而由于相比普通砂岩，页岩孔喉致密、有机碳及黏土矿物含量高导致孔隙结构特征必然存在明显差异，由于客观认识及实验方法上的限制，使得对T2C2的研究进展相对缓慢，而T2C2值的确定对评价页岩油后期生产动态具有重要意义。
发明内容
为明确页岩油储集层的孔隙结构特征及其孔隙流体赋存和运移规律，本发明提供一种页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法，采用核磁共振实验方法，在对目标页岩孔隙结构分类的基础上，开展离心试验与热处理实验，识别划分了不同孔隙结构页岩中3类孔隙流体的T2弛豫时间界限，定量给出了目标页岩的全孔径分布特征及多类流体识别方法，取得的成果有助于进一步认识页岩储层及识别划分页岩孔隙流体类型。
本发明所采用的技术方案是：
一种页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法，该方法具体按照如下步骤实施：
1)选取目标页岩储层中若干岩心，先后采用二氯甲苯和乙二醇溶剂反复清洗岩心，烘干后测量每块岩心的基础物性及矿物组成；
2)将岩心放入烘箱中进行高温烘干，以去除岩心中的水蒸气，然后对烘干后的岩心进行核磁共振扫描，获取岩心基础信号；
3)将岩心放入高压容器中抽真空，再以高压向容器中注入的KCl溶液避免样品发生水敏反应，高压饱和岩心后在常压下对饱和水的岩心进行核磁共振扫描，得到完全饱和水状态下的核磁共振T2谱；
4)根据完全饱和水状态下的核磁共振T2谱分布的形态，将目标储层页岩的孔隙结构划分为Ⅰ、Ⅱ和Ⅲ三种类型，三种类型对应的T2谱分布范围和波峰峰值依次降低，平均孔喉半径和孔隙中流体赋存量依次减小，孔喉连通性及渗流能力不断变差；
5)依次在不同转速条件下对岩心进行离心，并测量不同离心状态下岩心的核磁共振T2谱，确定最佳离心速度；
6)重复步骤3)，将最佳离心速度离心后的岩心放入不同温度烘箱中进行热处理，分别测量不同温度后岩心的核磁共振T2谱，并确定阈值温度；
7)根据步骤3)的页岩完全饱和水、步骤5)的最佳离心速度离心后、步骤6)的阈值温度热处理后三种状态下的T2谱分布，绘制对应的3条累积振幅随弛豫时间的变化曲线，然后反向延长最佳离心力离心后的累积振幅不变时的直线段，与完全饱和水状态下的累积振幅曲线的交点对应的弛豫时间即为可动水与束缚水的T2截止值(T2C1)，同理可以确定毛管束缚水与不可采出水的T2截止值(T2C2)；
8)将完全饱和水下和最佳离心力离心后的累积振幅之差再除以完全饱和水下的累积振幅即为可动水饱和度；同理，最佳离心力离心后和阈值温度热处理后的累积振幅之差再除以完全饱和水状态下的累积振幅即为毛管束缚水饱和度。
优选地，上述步骤2)中，放入烘箱中高温120℃烘干24h。
优选地，上述步骤3)中，将岩心放入高压容器中抽真空，再以30MPa高压向容器中注入1mol/L浓度的KCl溶液(此浓度能避免页岩发生水敏反应)，高压饱和岩心后在常压下对饱和水的岩心进行核磁共振扫描。
优选地，步骤5)中确定最佳离心力具体为：随着离心转速逐渐变大，岩心累积振幅曲线由下降速度较快到下降速度明显变缓，曲线上累积振幅明显变化的拐点对应的离心速度作为目标储层页岩的最佳离心转速。
优选地，上述步骤5)中，采用Optima L-100XP型超速离心机依次在转速6000、8000、10000、12000和14000r/min(对应离心力分别为1.18、2.76、4.11、5.36和71.8MPa)条件下对岩心进行离心，并测量不同离心状态下岩心的核磁共振T2谱。
优选地，步骤6)中，确定阈值温度具体为：随着温度逐渐升高，岩心累积振幅曲线由下降速度较快到下降速度明显变缓，曲线上累积振幅明显变化的拐点对应的离心速度作为目标储层页岩的最佳离心转速。
优选地，上述步骤6)中，将最佳离心力离心后的岩心放入不同温度烘箱中进行热处理，逐渐升高烘箱温度至40、60、80、100、120和140℃(当温度超过140℃后会对页岩孔隙中黏土矿物胶结方式产生破坏，固最高温度只到140℃)，分别测量不同温度后岩心的核磁共振T2谱。
本发明的有益效果：
本发明采用核磁共振实验方法，在对目标页岩孔隙结构分类的基础上，开展离心试验与热处理实验，识别划分了不同孔隙结构页岩中3类孔隙流体的T2弛豫时间界限，定量给出了目标页岩的全孔径分布特征及多类流体识别方法；采用发明方法可以明确页岩油储集层的孔隙结构特征及其孔隙流体赋存和运移规律，有助于进一步认识页岩储层及识别划分页岩孔隙流体类型，为页岩储层后续的开发提供了基础。
附图说明
图1是3种不同类型孔隙结构页岩的T2谱分布；
图2是3类孔隙结构代表性页岩在不同离心作用下的T2谱分布；
图3是3类代表性页岩累积振幅(归一化)随离心转速的变化规律；
图4是3类孔隙结构代表性页岩在不同温度热处理后的T2谱变化；
图5是3类孔隙结构代表性页岩累积振幅(归一化)随离心转速的变化规律；
图6是3类孔隙结构代表性页岩的双T2截止值。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明提供一种页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法，包括如下步骤：
1)选取目标页岩储层中若干岩心，先后采用二氯甲苯和乙二醇溶剂反复清洗岩心，烘干后测量每块岩心的基础物性及矿物组成；
2)将岩心放入烘箱中进行高温烘干，以去除岩心中的水蒸气，然后对烘干后的岩心进行核磁共振扫描，获取岩心基础信号；
3)将岩心放入高压容器中抽真空，再以高压向容器中注入的KCl溶液，高压饱和岩心后在常压下对饱和水的岩心进行核磁共振扫描，得到完全饱和水状态下的核磁共振T2谱；
4)根据完全饱和水状态下的核磁共振T2谱分布的形态，将目标储层页岩的孔隙结构划分为Ⅰ、Ⅱ和Ⅲ三种类型，三种类型对应的T2谱分布范围和波峰峰值依次降低，平均孔喉半径和孔隙中流体赋存量依次减小，孔喉连通性及渗流能力不断变差；
5)依次在不同转速条件下对岩心进行离心，并测量不同离心状态下岩心的核磁共振T2谱，确定最佳离心速度；
6)重复步骤3)，将最佳离心速度离心后的岩心放入不同温度烘箱中进行热处理，分别测量不同温度后岩心的核磁共振T2谱，并确定阈值温度；
7)根据步骤3)的页岩完全饱和水、步骤5)的最佳离心力离心后、步骤6)的阈值温度热处理后三种状态下的T2谱分布，绘制对应的3条累积振幅随弛豫时间的变化曲线，然后反向延长最佳离心力离心后的累积振幅不变时的直线段，与完全饱和水状态下的累积振幅曲线的交点对应的弛豫时间即为可动水与束缚水的T2截止值(T2C1)，同理可以确定毛管束缚水与不可采出水的T2截止值(T2C2)；
8)将完全饱和水下和最佳离心力离心后的累积振幅之差再除以完全饱和水下的累积振幅即为可动水饱和度；同理，最佳离心力离心后和阈值温度热处理后的累积振幅之差再除以完全饱和水状态下的累积振幅即为毛管束缚水饱和度。
进一步地，上述步骤5)中确定最佳离心力具体为：随着离心转速逐渐变大，岩心累积振幅曲线由下降速度较快到下降速度明显变缓，曲线上累积振幅明显变化的拐点对应的离心速度作为目标储层页岩的最佳离心转速。
进一步地，上述步骤6)中，确定阈值温度具体为：随着温度逐渐升高，岩心累积振幅曲线由下降速度较快到下降速度明显变缓，曲线上累积振幅明显变化的拐点对应的离心速度作为目标储层页岩的最佳离心转速。
进一步地，上述步骤2)中，需要把岩心放入烘箱中高温120℃烘干24h。
进一步地，上述步骤3)中，将岩心放入高压容器中抽真空，再以30MPa高压向容器中注入1mol/L浓度的KCl溶液(此浓度能避免页岩发生水敏反应)，高压饱和岩心后在常压下对饱和水的岩心进行核磁共振扫描。
进一步地，上述步骤5)中，采用Optima L-100XP型超速离心机依次在转速6000、8000、10000、12000和14000r/min(对应离心力分别为1.18、2.76、4.11、5.36和71.8MPa)条件下对岩心进行离心，并测量不同离心状态下岩心的核磁共振T2谱。
进一步地，上述步骤6)中，将最佳离心力离心后的岩心放入不同温度烘箱中进行热处理，逐渐升高烘箱温度至40、60、80、100、120和140℃(当温度超过140℃后会对页岩孔隙中黏土矿物胶结方式产生破坏，固最高温度只到140℃)，分别测量不同温度后岩心的核磁共振T2谱。
实施例：
对吴起油田某区位的页岩储层的微观孔隙结构特征及多类孔隙流体进行识别，以如下步骤进行实施：
步骤1：取吴起油田某区位井中的页岩，取样深度位于2684～2711m的长7油层组。从岩心柱上钻取若干块直径为2.5cm、长度为8cm左右的小岩心柱，然后分别测定岩心的基础物性和矿物组成及含量。目标页岩孔隙度介于4.4～14.6％，渗透率介于0.003～0.182×10-3μm2，有机碳含量(TOC)介于1.22～4.16％，石英含量介于9.5～41.9％，黏土矿物介于2.8～13.7％。实验页岩样品先后采用二氯甲苯和乙二醇溶剂反复清洗小岩心柱，烘干后测量每块岩心的基础物性及矿物组成。为避免切割过程中岩心发生碎裂，从所有岩心中选取15块结构完整且不含裂缝的岩心，再将每块岩心横向切割成3段(长度分别为1、2.5和4.5cm)，再次烘干后分别进行核磁共振实验。
步骤2：低场核磁共振实验
基于低场岩心核磁共振理论可知，岩石孔隙流体中的氢核会在外加磁场的作用下发生共振，通过测量氢核共振信号的幅度和弛豫速率可以间接表征岩石孔隙结构特征。通常，氢核的弛豫时间T2与孔隙半径成正比，共振信号的幅度与孔隙中流体赋存量成正比。因而采用核磁共振T2谱不但能够描述岩心的孔隙结构，还能定量评价孔隙中流体的分布特征。
本次核磁共振实验考虑到页岩纳米孔隙发育，T2谱中短弛豫孔隙占比高，因此采用0.069ms的回波间隔进行测量，其他参数分别为等待时间1.5s，扫描次数64，回波个数8192。
步骤2具体实验步骤如下：
步骤2.1：选取长度为4.5cm的岩心柱，放入烘箱中高温120℃烘干24h，以去除岩心中的水蒸气，然后对烘干后的岩心进行核磁共振扫描，获取岩心基础信号；
步骤2.2：对装有岩心的高压容器抽真空后，再以30MPa高压向容器中注入1mol/L浓度的KCl溶液(此浓度能避免页岩发生水敏反应)，高压饱和岩心后在常压下对饱和水的岩心进行核磁共振扫描；
步骤2.3：采用Optima L-100XP型超速离心机依次在转速6000、8000、10000、12000和14000r/min(对应离心力分别为1.18、2.76、4.11、5.36和71.8MPa)条件下对岩心进行离心，并测量不同离心状态下岩心的核磁共振T2谱；
步骤2.4：重复步骤2.2，将最佳离心力离心后的页岩放入不同温度烘箱中进行热处理，逐渐升高烘箱温度至40、60、80、100、120和140℃(当温度超过140℃后会对页岩孔隙中黏土矿物胶结方式产生破坏，固最高温度只到140℃)，分别测量不同温度后岩心的核磁共振T2谱。
步骤3：根据步骤2.2中页岩完全饱和水T2谱分布，确定页岩微观孔隙结构分类由15块页岩在完全饱和水状态下的T2谱分布(图1)可知，目标页岩具有较强的非均质性，根据T2谱分布的形态可以进一步将目标储层页岩的孔隙结构划分为3类。
Ⅰ类孔隙结构页岩的T2谱分布为左峰高于右峰的连续型双峰形态，弛豫时间分布范围广(0.01～2000ms)，平均孔隙度为12.8％，平均渗透率为0.141×10-3μm2，说明此类页岩储层物性相对较好，孔喉组合多样化，左峰代表的微小孔隙发育程度高于右峰代表的大孔隙和微裂缝发育程度，且微小孔隙与大孔隙(或微裂缝)之间连通性较好。
Ⅱ类孔隙结构页岩的T2谱分布属于左峰远大于右峰的不连续型双峰形态，弛豫时间范围变窄(0.01～200ms)，两波峰对应振幅明显降低，且两峰之间连通性变差，平均渗透率降至0.038×10-3μm2，说明此类页岩储层物性变差，平均孔隙半径减小，孔隙中流体赋存量也在不断降低。
Ⅲ类孔隙结构页岩的T2谱分布虽然以多峰分布为主，但主峰的分布范围大幅变窄，对应弛豫时间在0.01～1ms，平均孔隙度和渗透率均大幅下降，分别仅为5.3％和0.008×10-3μm2，说明此类页岩物性很差，孔隙类型单一，纳米孔隙极其发育，而大孔隙或微裂缝发育程度很低，孔隙中流体赋存量大幅降低。
步骤3：不同类型孔隙流体识别方法
步骤3.1：饱和-离心T2谱对比，得出各页岩岩心的最佳离心速度
图2为3类孔隙结构页岩中代表性页岩在不同离心力作用下的T2谱分布。从图2中可以看出，3块页岩T2谱中左右两峰对应振幅均随离心转速(离心力)的增加而不断降低，当离心转速由12000r/min(5.36MPa)增大至14000r/min(7.18MPa)时，Ⅰ类1#页岩和Ⅱ类9#页岩的T2振幅的下降幅度明显减小(Ⅲ类15#页岩由于不同离心力作用下振幅变化很小，未看出明显变化，需要进一步分析)。说明左峰代表的微小孔隙和右峰代表的大孔隙和微裂缝中均含有可动水，且随着离心力的增加，可动水在克服毛管压力后会优先被排出。当离心力增大至5.36MPa时，页岩孔隙中能够被离心出的可动水已基本达到最大，继续增大离心力非但不能大幅增加可动水的产出，还有可能导致页岩在水化作用下发生碎裂，而毛管束缚水与不可采出水仍滞留在孔隙中无法排出。
基于低场核磁共振测试中“核磁共振信号振幅的大小与页岩孔隙中流体的赋存量成正比”这一原理，可以采用T2谱中所有弛豫时间对应的振幅之和来表征页岩中流体的总赋存量。图3为不同离心转速下归一化后的累积振幅(归一化是指不同离心力下的累积振幅除以完全饱和水下的累积振幅)变化规律。从图3中可以看出，当离心转速由6000r/min增至12000r/min的过程中，3块页岩累积振幅的下降速度较快，而当离心转速由12000r/min增至14000r/min时，累积振幅的下降速度明显变缓。严格意义上而言，12000r/min并不是累积振幅明显变化的拐点，但却可以近似认为12000r/min是目标储层页岩的最佳离心转速。此外，Ⅰ类1#页岩累积振幅的降低幅度明显大于另外2块页岩，而Ⅲ类15#页岩累积振幅的降低幅度最小，说明1#页岩孔隙中可动水饱和度最大，而15#页岩中可动水饱和度最小，造成可动水饱和度差异的主要原因一方面是由于孔隙半径变小，孔喉组合方式变得单一，孔喉中的毛管压力增大；另一方面由于黏土矿物含量增大，更多的水被黏土矿物颗粒吸附，导致黏土束缚水饱和度增加，可动水饱和度大幅降低。
步骤3.2：饱和-热处理T2谱对比，确定阈值温度
离心试验仅能区分页岩孔隙中的可动水与束缚水，而无法有效识别束缚水中的毛管束缚水与不可采出水。毛管束缚水主要是指孔隙中水相在巨大毛管压力作用下被锁定，在最佳离心力作用下无法参与流动的水，但当页岩受热升温时，气水界面张力大幅减小，毛管压力的作用减弱，原来被毛细管力束缚的水相能够蒸发逃逸出来。因而可以采用热处理升温的方式来识别页岩中的毛管束缚水和不可采出水。
图4为3类孔隙结构页岩中代表性页岩在不同温度热处理后的T2谱分布。从图4中可以看出，在最佳离心力离心后，随着热处理温度的升高，3块页岩的T2振幅将会进一步大幅降低，主峰波峰对应的弛豫时间向左移动。结合图5可以进一步看出，1#页岩累积振幅的降低幅度明显大于另外2块页岩，而15#页岩的累积振幅降低幅度最小，这一现象与离心过程中累积振幅的变化规律相似，这主要因为Ⅰ类页岩孔喉半径大，孔喉连通性好，毛管束缚水受热蒸发排出过程中需要克服的毛管压力比另外2类页岩小。此外，还可以看出，当温度由40℃升高至100℃过程中，累积振幅的下降幅度较大，且下降速度较快，而当温度由100℃继续升高至140℃过程中，累积振幅的下降幅度明显变缓，并逐渐出现了拐点。与最佳离心力的选取方式一样，虽然100℃并不是严格意义上的阈值温度，但却可以近似选取100℃作为阈值温度。
步骤3.3：确定不同孔隙流体类型T2截止值
为准确识别页岩孔隙中可动水与束缚水、毛管束缚水与不可采出水之间的弛豫时间界限(即T2截止值)，首先根据页岩完全饱和水、最佳离心力离心后及阈值温度热处理后三种状态下的T2谱分布，绘制对应的3条累积振幅随弛豫时间的变化曲线(图6虚线)，然后反向延长最佳离心力离心后的累积振幅不变时的直线段(蓝色虚线)，与完全饱和水状态下的累积振幅曲线(黑色虚线)的交点对应的弛豫时间即为可动水与束缚水的T2截止值(T2C1)，同理可以确定毛管束缚水与不可采出水的T2截止值(T2C2)。此外，完全饱和水下和最佳离心力离心后的累积振幅之差再除以完全饱和水下的累积振幅即为可动水饱和度。同理，最佳离心力离心后和阈值温度热处理后的累积振幅之差再除以完全饱和水状态下的累积振幅即为毛管束缚水饱和度，而剩余饱和度即为不可采出水饱和度。
由图6可知，Ⅰ类孔隙结构页岩的T2C1在0.87～1.22ms，T2C2在0.08～0.26ms，对应的平均可动水饱和度、毛管束缚水饱和度和不可采出水饱和度分别为30.5％、34.2％和35.3％。Ⅱ类孔隙结构页岩的T2C1在0.98～1.52ms，T2C2在0.24～0.42ms，对应的平均可动水饱和度、毛管束缚水饱和度和不可采出水饱和度分别为26.6％、32.7％和40.7％。Ⅲ类孔隙结构页岩的T2C1在1.64～2.27ms，T2C2在0.41～0.72ms，对应的平均可动水饱和度、毛管束缚水饱和度和不可采出水饱和度分别为17.2％、32.7％和50.1％。综上可以看出，随着页岩孔隙结构不断变差，可动水和不可采出水的T2截止值不断增大，可动水饱和度不断降低、不可采出水饱和度不断增大，而毛管束缚水饱和度变化较小。
综上所述，根据本发明页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法，对吴起油田某区位的页岩储层进行分析得到如下结论：
(1)基于核磁共振T2谱分布形态，目标页岩孔隙结构可划分为Ⅰ、Ⅱ和Ⅲ三种类型，其对应的T2谱分布范围和波峰峰值依次降低，平均孔喉半径和孔隙中流体赋存量依次减小，孔喉连通性及渗流能力不断变差。
(2)Ⅰ类页岩孔隙以残余粒间孔、粒间溶蚀孔和微裂缝为主，喉道以孔隙缩小型为主；Ⅱ类页岩孔隙以粒间和粒内溶蚀孔为主，喉道以短导管状和缩颈型为主；Ⅲ类页岩孔隙以孔径小于100nm的晶间孔为主，孔喉间界限模糊，以树形结构连通方式为主。
(3)采用离心试验与热处理相结合的方法能够对页岩孔隙中的可动流体、毛管束缚流体与不可采出流体进行识别，Ⅰ、Ⅱ和Ⅲ类页岩的可动流体截止值T2C1分别为1.1、1.24和1.92ms，不可采出流体截止值T2C2分别为0.2、0.3和0.54ms。
(4)目标页岩全孔径分布特征为Ⅰ类页岩中的可动流体赋存于孔径大于24nm的中、大孔隙，平均饱和度在30.5％；不可采出流体则赋存于孔径小于4.5nm的微孔，平均饱和度在35.3％。Ⅱ类页岩中的可动流体赋存于孔径大于41.9nm的大孔隙中，平均饱和度在26.6％；不可采出流体赋存于孔径小于10.3nm的微、小孔隙，平均饱和度在40.7％。Ⅲ类页岩中的可动流体赋存于孔径大于93.6nm的大孔隙，平均饱和度在17.2％；不可采出流体则赋存于孔径小于26.3nm的中、小孔隙中，平均饱和度达到50.1％。
复杂的页岩孔隙结构特征导致多样的孔隙流体类型，上述实施例的研究结果明确了不同孔隙结构页岩中各类孔隙流体的赋存和运移特征，其为页岩全孔径孔隙流体划分及页岩油高效开发提供了依据。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。
Claims (6)
1.一种页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法，其特征在于，包括以下步骤：
1)选取目标页岩储层中若干岩心，先后采用二氯甲苯和乙二醇溶剂反复清洗岩心，烘干后测量每块岩心的基础物性及矿物组成；
2)将岩心放入烘箱中进行高温烘干，以去除岩心中的水蒸气；
3)将岩心放入高压容器中抽真空，再以高压向容器中注入的KCl溶液，高压饱和岩心后在常压下对饱和水的岩心进行核磁共振扫描，得到完全饱和水状态下的核磁共振T2谱；
4)根据完全饱和水状态下的核磁共振T2谱分布的形态，将目标储层页岩的孔隙结构划分为Ⅰ、Ⅱ和Ⅲ三种类型，三种类型对应的T2谱分布范围和波峰峰值依次降低，平均孔喉半径和孔隙中流体赋存量依次减小，孔喉连通性及渗流能力不断变差；
5)依次在不同转速条件下对岩心进行离心，并测量不同离心状态下岩心的核磁共振T2谱，确定最佳离心力；确定最佳离心力具体为：随着离心转速逐渐变大，岩心累积振幅曲线由下降速度较快到下降速度明显变缓，曲线上累积振幅明显变化的拐点对应的离心速度作为目标储层页岩的最佳离心转速；
6)重复步骤3)，将最佳离心速度离心后的岩心放入不同温度烘箱中进行热处理，分别测量不同温度后岩心的核磁共振T2谱，并确定阈值温度；
7)根据步骤3)的页岩完全饱和水、步骤5)的最佳离心速度离心后、步骤6)的阈值温度热处理后三种状态下的T2谱分布，绘制对应的3条累积振幅随弛豫时间的变化曲线，然后反向延长最佳离心速度离心后的累积振幅不变时的直线段，与完全饱和水状态下的累积振幅曲线的交点对应的弛豫时间即为可动水与束缚水的T2截止值T2C1，同理可以确定毛管束缚水与不可采出水的T2截止值T2C2。
2.根据权利要求1所述的一种页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法，其特征在于：
步骤2)中，将岩心放入烘箱中高温120℃烘干24h。
3.根据权利要求2所述的一种页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法，其特征在于：
步骤3)中，将岩心放入高压容器中抽真空，再以30MPa高压向容器中注入1mol/L浓度的KCl溶液，高压饱和岩心后在常压下对饱和水的岩心进行核磁共振扫描。
4.根据权利要求1所述的一种页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法，其特征在于：
步骤6)中，确定阈值温度具体为：随着温度逐渐升高，岩心累积振幅曲线由下降速度较快到下降速度明显变缓，曲线上累积振幅明显变化的拐点对应的离心速度作为目标储层页岩的最佳离心转速。
5.根据权利要求1所述的一种页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法，其特征在于：
步骤5)中，采用Optima L-100XP型超速离心机依次在转速6000、8000、10000、12000和14000r/min，对应离心力分别为1.18、2.76、4.11、5.36和71.8MPa条件下对岩心进行离心，并测量不同离心状态下岩心的核磁共振T2谱。
6.根据权利要求1所述的一种页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法，其特征在于：
步骤6)中，将最佳离心速度离心后的岩心放入不同温度烘箱中进行热处理，逐渐升高烘箱温度至40、60、80、100、120和140℃，分别测量不同温度后岩心的核磁共振T2谱。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110740441.5A CN113533156B (zh) | 2021-06-30 | 2021-06-30 | 页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110740441.5A CN113533156B (zh) | 2021-06-30 | 2021-06-30 | 页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113533156A CN113533156A (zh) | 2021-10-22 |
| CN113533156B true CN113533156B (zh) | 2024-02-02 |
Family
ID=78126419
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110740441.5A Active CN113533156B (zh) | 2021-06-30 | 2021-06-30 | 页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113533156B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115788373B (zh) * | 2021-11-26 | 2024-04-30 | 中国石油天然气集团有限公司 | 一种co2与页岩油藏油气置换规律评价方法 |
| CN114778588B (zh) * | 2022-04-27 | 2024-05-14 | 长江大学 | 基于nmr表征致密砂岩不同类型流体分布的方法 |
| CN116519731B (zh) * | 2023-07-03 | 2023-08-25 | 中国石油大学(华东) | 一种基于分子动力学模拟的页岩油可动用界限确定方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103760082A (zh) * | 2014-01-08 | 2014-04-30 | 中国石油天然气股份有限公司 | 页岩层系致密储层原油有效可动空间确定方法和装置 |
| CN104697915A (zh) * | 2015-03-20 | 2015-06-10 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 一种页岩微观孔隙大小及流体分布的分析方法 |
| CN105987925A (zh) * | 2015-02-28 | 2016-10-05 | 中国石油化工股份有限公司 | 油基钻屑的含油率和/或含油量的测量方法及其应用 |
| CN106525688A (zh) * | 2016-11-21 | 2017-03-22 | 中国石油大学(华东) | 一种饱和页岩孔隙流体分离及饱和度计算的实验方法 |
| CN108169099A (zh) * | 2018-01-17 | 2018-06-15 | 西南石油大学 | 一种基于核磁共振的页岩气储层孔隙结构定量计算方法 |
| CN109030311A (zh) * | 2018-07-16 | 2018-12-18 | 西南石油大学 | 基于核磁共振t2谱敏感参数的孔隙结构分类与识别方法 |
| CN110160934A (zh) * | 2019-06-06 | 2019-08-23 | 西安石油大学 | 基于核磁共振技术评价致密储层孔隙结构分类标准的方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10401313B2 (en) * | 2016-05-24 | 2019-09-03 | Schlumberger Technology Corporation | Estimating downhole fluid volumes using multi-dimensional nuclear magnetic resonance measurements |
| EP3593159B1 (en) * | 2017-03-07 | 2022-06-01 | Saudi Arabian Oil Company | Absolute porosity and pore size determination of pore types in media with varying pore sizes using nmr |
| CN108458960B (zh) * | 2018-03-27 | 2019-10-29 | 中国石油大学(华东) | 富有机质泥页岩含氢组分、孔隙度及孔径的评价方法 |
| CN108956417B (zh) * | 2018-06-09 | 2019-09-24 | 中国石油大学(华东) | 分析岩石孔隙无效吸水量的同位素核磁方法 |
-
2021
- 2021-06-30 CN CN202110740441.5A patent/CN113533156B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103760082A (zh) * | 2014-01-08 | 2014-04-30 | 中国石油天然气股份有限公司 | 页岩层系致密储层原油有效可动空间确定方法和装置 |
| CN105987925A (zh) * | 2015-02-28 | 2016-10-05 | 中国石油化工股份有限公司 | 油基钻屑的含油率和/或含油量的测量方法及其应用 |
| CN104697915A (zh) * | 2015-03-20 | 2015-06-10 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 一种页岩微观孔隙大小及流体分布的分析方法 |
| CN106525688A (zh) * | 2016-11-21 | 2017-03-22 | 中国石油大学(华东) | 一种饱和页岩孔隙流体分离及饱和度计算的实验方法 |
| CN108169099A (zh) * | 2018-01-17 | 2018-06-15 | 西南石油大学 | 一种基于核磁共振的页岩气储层孔隙结构定量计算方法 |
| CN109030311A (zh) * | 2018-07-16 | 2018-12-18 | 西南石油大学 | 基于核磁共振t2谱敏感参数的孔隙结构分类与识别方法 |
| CN110160934A (zh) * | 2019-06-06 | 2019-08-23 | 西安石油大学 | 基于核磁共振技术评价致密储层孔隙结构分类标准的方法 |
Non-Patent Citations (6)
| Title |
|---|
| Considerations for the acquisition and inversion of NMR T2 data in shales;M. Nadia Testamanti;Journal of Petroleum Science and Engineering;第174卷;全文 * |
| Experimental study of the NMR characteristics of shale reservoir rock;Sun Junchang Etal.;Science & Technology Review;第30卷(第14期);全文 * |
| 中国南方海相页岩储层可动流体及T_2截止值核磁共振研究;周尚文;刘洪林;闫刚;薛华庆;郭伟;;石油与天然气地质(第04期);全文 * |
| 基于核磁共振弛豫谱技术的页岩储层物性与流体特征研究;姚艳斌等;煤炭学报;第43卷(第1期);全文 * |
| 基于核磁共振驱替技术的页岩含油性价评价;何梦卿等;石油化工应用;第36卷(第11期);全文 * |
| 核磁共振录井可动流体t_2截止值确定方法;宋明会;;录井工程(03);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113533156A (zh) | 2021-10-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113533156B (zh) | 页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法 | |
| CN112858133A (zh) | 致密油藏裂缝动态渗吸排驱规律评价方法 | |
| CN111537541B (zh) | 致密油藏co2驱储层动用特征评价方法 | |
| CN107941838B (zh) | 一种二氧化碳驱油过程中沥青质沉淀对孔喉分布影响的定量评价方法 | |
| Humphry et al. | Impact of wettability on residual oil saturation and capillary desaturation curves | |
| Bai et al. | Effect of fracture on production characteristics and oil distribution during CO2 huff-n-puff under tight and low-permeability conditions | |
| CN109030311A (zh) | 基于核磁共振t2谱敏感参数的孔隙结构分类与识别方法 | |
| CN104075959A (zh) | 泥页岩油水饱和度精确测试方法 | |
| CN110541692A (zh) | 一种适用于超强水敏稠油油藏的化学驱油方法 | |
| CN112051190A (zh) | 一种复合渗透酸渗透能力评价方法 | |
| Huang et al. | The influence of CO2 huff and puff in tight oil reservoirs on pore structure characteristics and oil production from the microscopic scale | |
| CN113075102A (zh) | 一种建立多孔介质自发渗吸量与时间关系数学模型的方法 | |
| Wang et al. | Evaporation process and pore size distribution in tight sandstones: a study using NMR and MICP | |
| CN112414918B (zh) | 一种低渗储层注采井间砂体连通性的判别方法及装置 | |
| CN113933333A (zh) | 致密储层co2驱沥青质沉积对岩石润湿性变化特征评价方法 | |
| Fei et al. | Electrical responses and classification of complex waterflooded layers in carbonate reservoirs: A case study of Zananor Oilfield, Kazakhstan | |
| CN113176186A (zh) | 碳酸盐岩原油注入压力和孔喉下限的判断方法 | |
| Luo et al. | Simultaneous capillary pressure and wettability determination for tight bakken cores using an ultra-high-speed centrifuge | |
| Lu et al. | Trapping damage analysis of water locking/aqueous phase of tight sandstone gas reservoir in Ke Keya structure | |
| CN109899043B (zh) | 一种周期注聚提高原油采收率幅度的定量预测方法 | |
| CN114687730A (zh) | 一种水锁损害条件下的生产压差确定方法 | |
| Wang et al. | Study on movable fluid of low permeability reservoir with NMR technology | |
| CN114778588B (zh) | 基于nmr表征致密砂岩不同类型流体分布的方法 | |
| CN112782477A (zh) | 一种测定岩心不同润湿状态下电性响应特征的方法及系统 | |
| CN113916926A (zh) | 致密储层co2驱沥青质沉积作用下孔喉堵塞特征评价方法 |
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 |