CN115372398A - 页岩油气储层对外来流体敏感性定量评价新方法 - Google Patents
页岩油气储层对外来流体敏感性定量评价新方法 Download PDFInfo
- Publication number
- CN115372398A CN115372398A CN202210845363.XA CN202210845363A CN115372398A CN 115372398 A CN115372398 A CN 115372398A CN 202210845363 A CN202210845363 A CN 202210845363A CN 115372398 A CN115372398 A CN 115372398A
- Authority
- CN
- China
- Prior art keywords
- shale
- rock sample
- porosity
- nuclear magnetic
- sensitivity
- 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
- 239000012530 fluid Substances 0.000 title claims abstract description 206
- 230000035945 sensitivity Effects 0.000 title claims abstract description 133
- 238000000034 method Methods 0.000 title claims abstract description 52
- 239000003079 shale oil Substances 0.000 title claims abstract description 28
- 239000011435 rock Substances 0.000 claims abstract description 222
- 230000035699 permeability Effects 0.000 claims abstract description 107
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 59
- 238000011156 evaluation Methods 0.000 claims abstract description 55
- 238000011158 quantitative evaluation Methods 0.000 claims abstract description 55
- 239000011159 matrix material Substances 0.000 claims abstract description 50
- 238000012360 testing method Methods 0.000 claims abstract description 48
- 238000005259 measurement Methods 0.000 claims abstract description 16
- 239000011148 porous material Substances 0.000 claims description 101
- 238000001035 drying Methods 0.000 claims description 99
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 94
- 230000008859 change Effects 0.000 claims description 61
- 238000002791 soaking Methods 0.000 claims description 46
- 238000002474 experimental method Methods 0.000 claims description 41
- 238000007654 immersion Methods 0.000 claims description 40
- 230000003595 spectral effect Effects 0.000 claims description 39
- 239000002245 particle Substances 0.000 claims description 38
- 229920006395 saturated elastomer Polymers 0.000 claims description 38
- 230000001186 cumulative effect Effects 0.000 claims description 37
- 238000001228 spectrum Methods 0.000 claims description 29
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 23
- 229910052753 mercury Inorganic materials 0.000 claims description 23
- 238000009826 distribution Methods 0.000 claims description 20
- 238000002347 injection Methods 0.000 claims description 20
- 239000007924 injection Substances 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 20
- 238000006243 chemical reaction Methods 0.000 claims description 16
- 230000015572 biosynthetic process Effects 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 11
- 239000000126 substance Substances 0.000 claims description 11
- 239000012153 distilled water Substances 0.000 claims description 9
- 238000009825 accumulation Methods 0.000 claims description 7
- 238000011068 loading method Methods 0.000 claims description 6
- 238000003825 pressing Methods 0.000 claims description 6
- 238000002360 preparation method Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 3
- 238000000227 grinding Methods 0.000 claims description 3
- 229930195733 hydrocarbon Natural products 0.000 claims description 3
- 150000002430 hydrocarbons Chemical class 0.000 claims description 3
- 238000001225 nuclear magnetic resonance method Methods 0.000 claims description 3
- 239000005416 organic matter Substances 0.000 claims description 3
- 238000005192 partition Methods 0.000 claims description 3
- 238000002459 porosimetry Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000002156 mixing Methods 0.000 claims description 2
- 238000003828 vacuum filtration Methods 0.000 claims description 2
- 239000004215 Carbon black (E152) Substances 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 6
- 238000011161 development Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000008595 infiltration Effects 0.000 description 2
- 238000001764 infiltration Methods 0.000 description 2
- 241001391944 Commicarpus scandens Species 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 238000000643 oven drying Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000010998 test 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
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
- G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/34—Purifying; Cleaning
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
- G01N1/44—Sample treatment involving radiation, e.g. heat
-
- 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
-
- 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
- G01N15/0886—Mercury 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N5/00—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid
- G01N5/04—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid by removing a component, e.g. by evaporation, and weighing the remainder
-
- 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
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- High Energy & Nuclear Physics (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Geochemistry & Mineralogy (AREA)
- Dispersion Chemistry (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Remote Sensing (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种页岩油气储层对外来流体敏感性定量评价新方法，该方法包括以下步骤：(一)、岩样制备；(二)、页岩孔隙度和层理缝渗透率对外来流体敏感性定量评价；(三)、页岩基质渗透率对外来流体敏感性定量评价。基于核磁共振测试技术，从孔隙度和渗透率这两个基础物性参数出发，创新提出了页岩油气储层对外来流体敏感性综合评价新方法，实现了页岩不同尺寸孔隙、基质和层理缝渗透率对外来流体敏感性的定量评价，解决了现有标准SY/T5358‑2010《储层敏感性流动实验评价方法》针对页岩岩样无法建立有效驱动体系、实验测量误差大和评价指标单一的问题，实现了页岩油气储层对外来流体敏感性定量精确评价。
Description
技术领域：
本发明涉及非常规页岩油气储层勘探开发领域，具体涉及一种页岩油气储层对外来流体敏感性评价新方法。
背景技术：
随着我国油气资源勘探开发技术日趋成熟，油气藏开发物性下限逐渐向非常规油气资源靠拢，页岩油气等非常规油气资源表现出巨大开发潜力，正逐渐成为我国重要接替战略资源。国内外目前普遍采用水平井体积压裂模式，同时结合衰竭式开采方式，来实现非常规页岩油气资源的有效开发。页岩油气储层压裂过程中注入大规模外来流体，不同类型外来流体对页岩储层孔隙结构和渗透率造成的伤害或改善程度需定量评价，因此有必要系统开展页岩油气储层对外来流体敏感性定量评价研究。
目前现有储层对外来流体敏感性定量评价方法主要是执行标准SY/T 5358-2010《储层敏感性流动实验评价方法》，该标准测量原理是基于传统达西定律，通过驱替方式测量岩样与外来流体接触作用前后渗透率的相对变化率来评价储层敏感性。
页岩油气储层属于超致密油气储层，微—纳米孔隙大量发育，基质渗透率远远低于常规砂岩储层。采用标准SY/T 5358-2010《储层敏感性流动实验评价方法》无法实现页岩油气储层对外来流体敏感性定量评价，主要原因是：
(1)针对页岩岩样无法建立有效驱动体系。实验过程中若采用大压差建立驱动体系，需要对岩样施加大数值围压，由于页岩岩样脆性较强，在大数值围压作用下，岩样容易破碎产生裂隙，进而影响实验测量精度。(2)页岩基质和层理缝渗透率非常低，需要较长时间达到稳定流动且测量时需要人为记录压差和流速等实验数据，误差相对较大。(3)现有标准SY/T 5358-2010《储层敏感性流动实验评价方法》仅通过测量外来流体作用前后岩样渗透率的相对变化评价储层对外来流体敏感性，评价指标单一。
针对现有标准的三个技术弊端，有必要建立一套适用于页岩油气储层敏感性定量评价新方法。本发明从孔隙度和渗透率这两个基础物性参数出发，基于核磁共振技术提出了一种页岩油气储层对外来流体敏感性评价新方法，实现了页岩这类超致密油气储层对外来流体敏感性的全面准确评价，为页岩油气储层对外来流体敏感性定量评价提供了一种全新的科学手段。
发明内容：
为了解决背景技术中存在的问题，本发明提出了一种页岩油气储层对外来流体敏感性评价新方法。本申请基于核磁共振测试技术，从孔隙度和渗透率这两个基础物性参数出发，创新提出了页岩油气储层对外来流体敏感性综合评价新方法，实现了页岩不同尺寸孔隙、基质和层理缝渗透率对外来流体敏感性的定量评价，解决了现有标准SY/T 5358-2010《储层敏感性流动实验评价方法》针对页岩岩样无法建立有效驱动体系、实验测量误差大和评价指标单一的问题，实现了页岩油气储层对外来流体敏感性定量精确评价。
本发明的技术方案是：一种页岩油气储层对外来流体敏感性定量评价新方法，该方法包括以下步骤：
(一)、岩样制备
选择来自同一取芯层段的页岩岩样3块，并将其标号为1#、2#和3#，其中1#和2#岩样用于确定目标烘干温度和时间，3#岩样用于完成孔隙度和渗透率对外来流体敏感性评价实验，3#页岩岩样要求：岩样为标准柱塞样，直径为2.5cm且长度应大于或等于5cm，实验前需对3#岩样进行预处理，将岩样切割成3块小岩样，并将其编号为4#、5#和6#，4#岩样用于高压压汞测试、5#岩样用于开展页岩孔隙度和层理缝渗透率敏感性评价实验、6#岩样研磨成10/20目颗粒用于页岩基质渗透率敏感性评价实验；
(二)、页岩孔隙度和层理缝渗透率对外来流体敏感性定量评价
页岩孔隙度和层理缝渗透率对外来流体敏感性定量评价过程包含以下3个步骤：
A、确定页岩岩样目标烘干温度和时间
在开展页岩孔隙度和层理缝渗透率对外来流体敏感性评价实验之前，需清除岩样中存有的烃类、水等可动流体，确保岩样清洁；
(1)选择来自同一取芯层段的页岩岩样2块，并将其标号为1#和2#；
(2)设定10个不同烘干温度Tt1、Tt2、Tt3、Tt4、Tt5、Tt6、Tt7、Tt8、Tt9、Tt10，温度按照20℃依次递增；
(3)在设定温度条件下烘干1#页岩岩样t0小时，测量1#页岩岩样在不同烘干温度条件下对应的岩样质量m1、m2、m3、m4、m5、m6、m7、m8、m9、m10并开展烘干后页岩二维核磁共振测试；
(4)绘制1#页岩岩样质量随烘干温度变化曲线；当页岩岩样内可动流体完全清除时，岩样质量随烘干温度变化曲线出现拐点，拐点对应的温度为目标烘干温度T0；
(5)设定10个不同烘干时间t1、t2、t3、t4、t5、t6、t7、t8、t9、t10，时间按照2小时依次递增；在确保实验准确度条件下，为减少测试工作量，可将t5或者t6时间设定与t0相同；
(6)在目标烘干温度T0条件下按照设定时间烘干2#页岩岩样，测量页岩岩样在不同烘干时间条件下对应的岩样质量m'1、m'2、m'3、m'4、m'5、m'6、m'7、m'8、m'9、m'10并开展烘干后页岩二维核磁共振测试；
(7)绘制2#页岩岩样质量随烘干时间变化曲线；当页岩岩样内可动流体完全清除时，岩样质量不随烘干时间变化，此时对应的时间为目标烘干时间t0；
(8)分析大于目标烘干时间t0之后的二维核磁共振T1和T2图谱，若小、中、大孔内的核磁信号量不发生改变，则可确定页岩岩样最终的目标烘干温度和烘干时间；若二维核磁共振T1和T2图谱中小、中、大孔内的核磁信号量发生改变，则需重新筛选岩样按照步骤(1)-(7)开展实验，再次确定页岩岩样目标烘干温度T0和时间t0；
B、开展页岩岩样高压压汞+外来流体浸泡+核磁共振联测实验
(1)页岩岩样高压压汞测试
将制备好的4#页岩岩样按照确定的目标烘干温度T0和时间t0进行烘干处理，将烘干后的岩样开展高压压汞测试，在半对数坐标条件下绘制孔隙分布频率随孔隙半径变化曲线，确定页岩岩样孔隙分布特征；
(2)页岩岩样与外来流体浸泡+核磁共振联测实验；
C、实验数据处理与分析
1)页岩层理缝渗透率对外来流体敏感性定量评价
(1)引入层理缝渗透率对外来流体敏感性定量评价指标，如公式(1)所示：
其中，RK—层理缝渗透率变化率，％；Kb—页岩干样层理缝渗透率，mD；Ka—页岩浸泡外来流体后层理缝渗透率，mD；
(2)按照公式(1)计算页岩浸泡外来流体前后层理缝渗透率变化率；
(3)参考SY/T 5358-2010《储层敏感性流动实验评价方法》敏感性影响程度评价指标，定量评价页岩层理缝渗透率对外来流体敏感性；
2)页岩总孔隙度对外来流体敏感性定量评价
(1)引入页岩总孔隙度对外来流体敏感性定量评价指标，如公式(2)所示：
其中，Rt—覆压孔渗仪法测量的总孔隙度变化率，％；φb—页岩干样孔隙度，％；φa—页岩浸泡外来流体后孔隙度，％；
(2)按照公式(2)计算页岩浸泡外来流体前后孔隙度变化率；
(3)参考SY/T 5358-2010《储层敏感性流动实验评价方法》敏感性影响程度评价指标，定量评价页岩总孔隙度对外来流体敏感性；
3)页岩不同尺寸孔隙对外来流体敏感性定量评价
(1)核磁共振弛豫时间T2值和孔隙尺寸时空转换系数计算；
(2)核磁共振T2谱信号量与页岩饱和水孔隙度标定关系建立
①基于重量法计算5#页岩岩样饱和水孔隙度，计算公式如下：
其中，φw—饱和水孔隙度，％；mw—饱和水后岩样重量，g；m0—烘干后岩样重量，g；ρw—水密度，g/cm3；d—岩样直径，cm；L—岩样长度，cm；
②基于核磁共振测试T2谱信号分量，按照公式(4)将页岩饱和水后T2谱核磁信号分量转化为孔隙度分量：
其中，
—饱和水后孔隙度分量，％；Sw—饱和水后T2谱核磁信号分量，PU；Sac,w—饱和水后T2谱核磁信号累计量，PU；Sac,d—饱和水T2谱第一个峰对应的核磁信号累计量，PU；
③基于步骤②饱和水后T2谱核磁信号分量与孔隙度分量转换原理，按照公式(5)将页岩浸泡外来流体后T2谱核磁信号分量转化为孔隙度分量：
其中，
—外来流体浸泡后孔隙度分量，％；Ss—外来流体浸泡后T2谱核磁信号分量，PU；
④基于步骤②和③得到的孔隙度分量，绘制页岩饱和水和浸泡外来流体后孔隙度分量和孔隙度分量累计值随T2弛豫时间变化关系曲线；需要指出的是，上述曲线需去除T2谱内第一个峰对应的核磁信号累计量；
(3)不同尺寸孔隙对外来流体敏感性定量评价
①基于孔隙度分量累计值随T2弛豫时间变化关系曲线，引入页岩总孔隙度和不同尺寸孔隙度对外来流体敏感性定量评价指标，如公式(6)、(7)、(8)、(9)和(10)所示：
其中，
—核磁共振法计算的总孔隙度变化率，％；
—外来流体浸泡后总孔隙度分量累计值，％；
—饱和水后总孔隙度分量累计值，％；
—微孔孔隙度变化率，％；
—外来流体浸泡后微孔孔隙度分量累计值，％；
—饱和水后微孔孔隙度分量累计值，％；
—小孔孔隙度变化率，％；
—外来流体浸泡后小孔孔隙度分量累计值，％；
—饱和水后小孔孔隙度分量累计值，％；
—中孔孔隙度变化率，％；
—外来流体浸泡后中孔孔隙度分量累计值，％；
—饱和水后中孔孔隙度分量累计值，％；
—大孔孔隙度变化率，％；
—外来流体浸泡后大孔孔隙度分量累计值，％；
—饱和水后大孔孔隙度分量累计值，％；
需要指出的是，针对页岩储层推荐不同尺寸孔隙划分标准为：孔隙半径＜0.01μm为微孔、0.01-0.1μm为小孔、0.1-1.0μm为中孔、孔隙半径＞1.0μm为大孔；
②基于时空转换系数k，计算微孔、小孔、中孔和大孔孔隙半径划分对应的T2弛豫时间值，即10nm、100nm和1000nm对应的T2弛豫时间值T2,10nm、T2,100nm和T2,1000nm；
③基于T2,10nm、T2,100nm和T2,1000nm值，利用公式(11)-(18)计算饱和水后微孔、小孔、中孔和大孔孔隙度
和外来流体浸泡后微孔、小孔、中孔和大孔孔隙度
其中，
—饱和水后微孔T2谱核磁信号累计值，PU；
—饱和水后小孔T2谱核磁信号累计值，PU；
—饱和水后中孔T2谱核磁信号累计值，PU；
—外来流体浸泡后微孔T2谱核磁信号累计值，PU；
—外来流体浸泡后小孔T2谱核磁信号累计值，PU；
—外来流体浸泡后中孔T2谱核磁信号累计值，PU；Sac,sd—外来流体浸泡后T2谱第一个峰对应的核磁信号累计量，PU；Sac,s—外来流体浸泡后T2谱核磁信号累计量，PU；
④基于公式(6)、(7)、(8)、(9)和(10)计算页岩岩样浸泡外来流体前后总孔隙度和不同尺寸孔隙度变化率，参考SY/T 5358-2010《储层敏感性流动实验评价方法》敏感性影响程度评价指标，给出页岩储层总孔隙度和不同尺寸孔隙对外来流体敏感性定量评价结果；
(三)、页岩基质渗透率对外来流体敏感性定量评价
页岩基质渗透率对外来流体敏感性定量评价包含以下步骤：
①样品制备：选取与孔隙度和层理缝渗透率对外来流体敏感性评价同一取芯层段的页岩岩样，将样品粉碎研磨至10/20目颗粒，质量大于30g；
②在目标烘干温度T0和时间t0条件下对颗粒样品进行烘干处理；
③基于压力脉冲衰减原理，利用SMP-200页岩基质渗透率仪测定6#样品研磨后颗粒的基质渗透率Kmb；
④将颗粒样品放置于活塞容器内，利用抽真空预处理系统对其抽真空，时间为48小时，为保护活塞容器，在其上部放置筛网；
⑤将预先配制好的外来流体放置于另一个活塞容器内，利用ISCO泵将外来流体注入到盛放颗粒样品的活塞容器内浸泡48小时，浸泡压力与地层压力相同；
⑥取出颗粒样品，将其在目标烘干温度T0和时间t0条件下烘干；
⑦利用SMP-200页岩基质渗透率仪测定外来流体浸泡后颗粒样品的基质渗透率Kma；
⑧基于公式(19)计算页岩颗粒样品浸泡外来流体前后基质渗透率变化率：
其中，RmK—页岩颗粒样品基质渗透率变化率，％；Kmb—页岩颗粒样品浸泡外来流体前(烘干后)基质渗透率，mD；Kma—页岩颗粒样品浸泡外来流体后基质渗透率，mD；
⑨参考SY/T 5358-2010《储层敏感性流动实验评价方法》敏感性影响程度评价指标，给出页岩基质渗透率对外来流体敏感性定量评价结果；
基于页岩浸泡外来流体前后总孔隙度、不同尺寸孔隙度、层理缝渗透率和基质渗透率变化规律，综合分析给出页岩油气储层对外来流体敏感性定量评价结果。
上述方案中的步骤B中第(2)步页岩岩样与外来流体浸泡+核磁共振联测实验包括以下步骤：
①将5#岩样按照确定的目标烘干温度T0和时间t0进行烘干处理并记录其长度L、直径d和质量m0；
②利用VINCI覆压孔渗仪测量烘干后岩样的孔隙度φb和层理缝渗透率Kb；
③将测量孔渗后的标准岩样放置于岩心夹持器内，利用手摇泵加载2MPa围压，之后利用抽真空预处理系统对岩样进行抽真空，时间为48小时；
④制备实验用蒸馏水并将其放置于活塞容器内，利用ISCO泵将活塞容器内蒸馏水以恒压方式注入岩样内，饱和过程中注入压力与围压呈阶梯方式递增且保持围压与注压差值为2MPa恒定，当注压达到地层压力时，停止饱和过程；岩样饱和过程在恒温箱内完成，温度与地层温度保持一致，总饱和时间不低于48小时；
⑤将饱和水后岩样从岩心夹持器内取出并记录其质量mw，测试岩样饱和水后核磁共振T2谱，为减小实验误差，需连续测量3次；
⑥将步骤⑤饱和水后岩样在目标烘干温度T0和时间t0进行烘干处理；
⑦将烘干后岩样放置于岩心夹持器内，利用手摇泵加载2MPa围压，之后利用抽真空预处理系统对岩样进行抽真空，时间为48小时；
⑧利用ISCO泵将活塞容器内外来流体以恒压方式注入岩样内，注入压力与围压呈阶梯方式递增且保持围压与注压差值为2MPa恒定，当注压达到地层压力时，停止饱和过程，总饱和时间不低于48小时；岩样饱和与浸泡外来流体过程在恒温箱内完成，温度与地层温度保持一致，浸泡时间为48小时；
⑨开展岩样浸泡外来流体后核磁共振T2谱测试，为减小实验误差，需连续测量3次；
⑩将步骤⑨浸泡外来流体后的页岩岩样在目标烘干温度T0和时间t0进行烘干处理；
上述方案中的步骤C中第3)步中的第(1)步核磁共振弛豫时间T2值和孔隙尺寸时空转换系数计算中包括以下步骤；
①提取5#岩样饱和水后核磁共振T2弛豫时间与核磁信号量数据，在半对数坐标条件下绘制核磁信号量随T2弛豫时间变化曲线；
②提取高压压汞测试结果中的孔隙半径和孔隙分布频率数据，在半对数坐标条件下绘制孔隙分布频率随孔隙半径变化曲线；
③将步骤①和②曲线数据整合在同一坐标系内，建立核磁共振与高压压汞时空转换曲线，其中X轴底部坐标为孔隙半径，Y轴主坐标为孔隙分布频率；X轴顶部坐标为T2弛豫时间，Y轴副坐标为核磁信号量；
④页岩岩样饱和水后核磁共振T2谱第一个峰代表有机质信号显示，因此以第二个峰对应的弛豫时间开始进行时空转换系数计算，当核磁信号峰值与高压压汞孔隙分布频率峰值逐一对应时，记录此时的核磁共振弛豫时间T2i和孔隙半径ri值；基于以上数据，即可计算页岩岩样核磁共振弛豫时间T2值和孔隙尺寸时空转换系数k＝(k1+k2+…+kn)/n,其中n代表核磁信号与高压压汞孔隙分布频率对应峰的数量。
本发明具有如下有益效果：本申请基于核磁共振测试技术，从孔隙度和渗透率这两个基础物性参数出发，创新提出了页岩油气储层对外来流体敏感性定量评价新方法，解决了现有标准SY/T 5358-2010《储层敏感性流动实验评价方法》针对页岩岩样无法建立有效驱动体系、实验测量时间长和误差大、敏感性评价指标单一的问题。本发明提出的高压压汞+外来流体浸泡+核磁共振联测实验方法可便捷、有效建立储层压力和温度模拟条件，具有实验测量周期短、精度高的优点。与此同时，本发明可分别给出页岩总孔隙、不同尺寸孔隙、基质和层理缝渗透率对外来流体敏感性的定量评价结果，实现了页岩油气储层对外来流体敏感性定量综合精确评价。
本发明从孔隙度和渗透率这两个基础物性参数出发，基于核磁共振技术提出了一种页岩油气储层对外来流体敏感性评价新方法，实现了页岩这类超致密油气储层对外来流体敏感性的全面准确评价，为页岩油气储层对外来流体敏感性定量评价提供了一种全新的科学手段。
附图说明：
附图1：实验用页岩岩样；
附图2：1#页岩烘干后二维核磁共振测试图谱(烘干温度110℃、烘干时间8小时)；
附图3：1#页岩岩样质量随烘干温度变化曲线图；
附图4：2#页岩烘干后二维核磁共振测试图谱(烘干温度110℃、烘干时间10小时)；
附图5：2#页岩岩样质量随烘干时间变化曲线图；
附图6：核磁共振弛豫时间T2值和孔隙尺寸时空转换曲线图；
附图7：孔隙度分量和孔隙度分量累计值随T2弛豫时间变化关系曲线图。
具体实施方式：
下面结合实例对本发明作进一步说明：
实施例：
一、岩样制备
选择来自松辽盆地青山口组青一段页岩岩样3块，并将其标号为1#、2#和3#。其中1#和2#岩样用于确定目标烘干温度和时间，3#岩样用于完成孔隙度和渗透率对外来流体敏感性评价实验。实验前对3#岩样进行预处理，将岩样切割成3块小岩样，并将其编号为4#、5#和6#，4#岩样用于高压压汞测试、5#岩样用于开展页岩孔隙度和层理缝渗透率敏感性评价实验、6#岩样研磨成10/20目颗粒用于页岩基质渗透率敏感性评价实验。实验用页岩岩样参见附图1。
二、页岩孔隙度和层理缝渗透率对外来流体敏感性定量评价
1、确定页岩岩样目标烘干温度和时间
(1)选择来自青山口组青一段页岩岩样2块，标号为1#和2#(见附图1)。
(2)设定烘干温度50、70、90、110、130、150、170、190、210、230℃。
(3)在设定温度条件下烘干1#页岩岩样8小时，测量其在不同烘干温度条件下对应的岩样质量并开展烘干后1#页岩二维核磁共振测试(见附图2)。
(4)绘制1#页岩岩样质量随烘干温度变化曲线(见附图3)。当页岩岩样内可动流体完全清除时，岩样质量随烘干温度变化曲线出现拐点，拐点对应的温度为目标烘干温度T0＝110℃。
(5)设定烘干时间2、4、6、8、10、12、14、16、18、20小时。
(6)在目标烘干温度110℃条件下按照设定时间烘干2#页岩岩样，测量页岩岩样在不同烘干时间条件下对应的岩样质量并开展烘干后页岩二维核磁共振测试(见附图4)。
(7)绘制2#页岩岩样质量随烘干时间变化曲线(见附图5)。当页岩岩样内可动流体完全清除时，岩样质量不随烘干时间变化，此时对应的时间为目标烘干时间t0＝8小时。
(8)由附图2和附图4对比分析可知：当烘干时间为10小时(大于目标烘干时间8小时)，二维核磁共振T1和T2图谱中小、中、大孔内的核磁信号量不发生改变，则可最终确定页岩岩样的目标烘干温度为110℃和烘干时间为8小时。
2、开展页岩岩样高压压汞+外来流体浸泡+核磁共振联测实验
(1)页岩岩样高压压汞测试
将4#页岩岩样在温度110℃条件下烘干8小时后开展高压压汞测试，在半对数坐标条件下绘制孔隙分布频率随孔隙半径变化曲线(见附图6)，确定页岩岩样孔隙分布特征。
(2)页岩岩样与外来流体浸泡+核磁共振联测实验
①将5#岩样在温度110℃条件下烘干8小时并测量其长度L＝3.02cm、直径d＝2.50cm和质量m0＝35.0320g；
②利用VINCI覆压孔渗仪测量烘干后页岩岩样的孔隙度φb＝6.854％和层理缝渗透率Kb＝0.184mD；
③将测量孔渗后的标准页岩岩样放置于岩心夹持器内，利用手摇泵加载2MPa围压，之后利用抽真空预处理系统对岩样进行抽真空，时间为48小时；
④制备实验用蒸馏水并将其放置于活塞容器内，利用ISCO泵将活塞容器内蒸馏水以恒压方式注入岩样内，饱和过程中注入压力与围压呈阶梯方式递增且保持围压与注压差值为2MPa恒定，当注压达到青一段地层压力时，停止饱和过程。岩样饱和过程在恒温箱内完成，温度与青一段地层温度保持一致，总饱和时间为48小时；
⑤将饱和水后岩样从岩心夹持器内取出并记录其质量mw＝35.8662g，测试岩样饱和水后核磁共振T2谱(见附图6)，该T2谱为三次测试结果平均值；
⑥将步骤⑤饱和水后页岩岩样在温度110℃条件下烘干8小时；
⑦将烘干后岩样放置于岩心夹持器内，利用手摇泵加载2MPa围压，之后利用抽真空预处理系统对岩样进行抽真空，时间为48小时；
⑧利用ISCO泵将活塞容器内外来流体(5％HCl)以恒压方式注入岩样内，注入压力与围压呈阶梯方式递增且保持围压与注压差值为2MPa恒定，当注压达到青一段地层压力时，停止饱和过程，总饱和时间为48小时。岩样饱和与浸泡外来流体(5％HCl)过程在恒温箱内完成，温度与青一段地层温度保持一致，浸泡时间为48小时；
⑨测试岩样浸泡外来流体(5％HCl)后核磁共振T2谱(见附图7)。需要指出的是，附图7内曲线(X轴为弛豫时间T2，Y轴为孔隙度分量)代表岩样浸泡5％HCl后孔隙度分量，T2谱核磁信号分量变化规律与孔隙度分量变化规律完全相同，二者数值上相差一个系数，因此浸泡外来流体(5％HCl)后的T2谱没有重复给出；
⑩将步骤⑨浸泡外来流体(5％HCl)后的页岩岩样在温度110℃条件下烘干8小时；
3、实验数据处理与分析
1)页岩层理缝渗透率对外来流体敏感性定量评价
(1)引入层理缝渗透率对外来流体敏感性定量评价指标，如公式(1)所示：
其中，RK—层理缝渗透率变化率，％；Kb—页岩干样层理缝渗透率，mD；Ka—页岩浸泡外来流体后层理缝渗透率，mD；
(2)按照公式(1)计算页岩浸泡外来流体前后层理缝渗透率变化率：Kb＝0.184mD、Ka＝0.196mD、RK＝6.52％。
(3)参考SY/T 5358-2010《储层敏感性流动实验评价方法》敏感性影响程度评价指标(见表1)，页岩层理缝渗透率对外来流体(5％HCl)敏感性为改善(弱)，页岩浸泡5％HCl后层理缝渗透率增加0.012mD。
表1敏感性影响程度评价指标
| 变化率R<sub>k</sub>(％) | 敏感性影响程度 |
| R<sub>k</sub>≤5 | 无 |
| 5＜R<sub>k</sub>≤30 | 弱 |
| 30＜R<sub>k</sub>≤50 | 中等偏弱 |
| 50＜R<sub>k</sub>≤70 | 中等偏强 |
| R<sub>k</sub>≥70 | 强 |
附注：若Rk为“+”值，则代表改善；若Rk为“-”值，则代表伤害。
2)页岩总孔隙度对外来流体敏感性定量评价
(1)引入页岩总孔隙度对外来流体敏感性定量评价指标，如公式(2)所示：
其中，Rt—覆压孔渗仪法测量的总孔隙度变化率，％；φb—页岩干样孔隙度，％；φa—页岩浸泡外来流体后孔隙度，％。
(2)按照公式(2)计算页岩浸泡外来流体前后孔隙度变化率：φb＝6.854％、φa＝7.571％、Rt＝10.46％。
(3)参考SY/T 5358-2010《储层敏感性流动实验评价方法》敏感性影响程度评价指标(见表1)，页岩总孔隙度对外来流体(5％HCl)敏感性为改善(弱)，页岩浸泡5％HCl后总孔隙度增加0.717％。
3)页岩不同尺寸孔隙对外来流体敏感性定量评价
(1)核磁共振弛豫时间T2值和孔隙尺寸时空转换系数计算
①提取5#岩样饱和水后核磁共振T2弛豫时间与核磁信号量数据，在半对数坐标条件下绘制核磁信号量随T2弛豫时间变化曲线(见附图6)。
②提取高压压汞测试结果中的孔隙半径和孔隙分布频率数据(见附图6)，在半对数坐标条件下绘制孔隙分布频率随孔隙半径变化曲线(见附图6)。
③将步骤①和②曲线数据整合在同一坐标系内，建立核磁共振与高压压汞时空转换曲线(见附图6)，其中X轴底部坐标为孔隙半径，Y轴主坐标为孔隙分布频率；X轴顶部坐标为T2弛豫时间，Y轴副坐标为核磁信号量。
④页岩岩样饱和水后核磁共振T2谱第一个峰代表有机质信号显示，因此以第二个峰对应的弛豫时间开始进行时空转换系数计算，当核磁信号峰值与高压压汞孔隙分布频率峰值逐一对应时，记录此时的核磁共振弛豫时间T2i和孔隙半径ri值(见表2)。基于表2数据，即可计算页岩岩样核磁共振弛豫时间T2值和孔隙尺寸时空转换系数k＝(9.49+9.00+7.07)/3＝8.52nm/ms。
表2孔隙半径与弛豫时间对照数据表
| r(nm) | T<sub>2</sub>(ms) | k(nm/ms) |
| r<sub>1</sub>＝4 | T<sub>21</sub>＝0.4217 | k<sub>1</sub>＝9.49 |
| r<sub>2</sub>＝16 | T<sub>22</sub>＝1.7783 | k<sub>2</sub>＝9.00 |
| r<sub>3</sub>＝63 | T<sub>23</sub>＝8.9125 | k<sub>3</sub>＝7.07 |
(2)核磁共振T2谱信号量与页岩饱和水孔隙度标定关系建立
①基于重量法计算5#页岩岩样饱和水孔隙度，计算公式如下：
其中，φw—饱和水孔隙度，％；mw—饱和水后岩样重量，g；m0—烘干后岩样重量，g；ρw—水密度，g/cm3；d—岩样直径，cm；L—岩样长度，cm。5#页岩岩样饱和水孔隙度
②基于核磁共振测试T2谱信号分量，按照公式(4)将页岩饱和水后T2谱核磁信号分量转化为孔隙度分量：
其中，
—饱和水后孔隙度分量(见附图7)，％；Sw—饱和水后T2谱核磁信号分量，PU；Sac,w—饱和水后T2谱核磁信号累计量，PU；Sac,d—饱和水T2谱第一个峰对应的核磁信号累计量，PU。
根据附图6内T2谱第一个峰对应的核磁信号累计量计算可得：Sac,d＝1670.7PU、Sac,w＝3349.1PU。
③基于步骤②饱和水后T2谱核磁信号分量与孔隙度分量转换原理，按照公式(5)将页岩浸泡外来流体(5％HCl)后T2谱核磁信号分量转化为孔隙度分量：
其中，
—外来流体浸泡后孔隙度分量(见附图7)，％；Ss—外来流体浸泡后T2谱核磁信号分量(见附图7)，PU。
④基于步骤②和③得到的孔隙度分量，绘制页岩饱和水和浸泡外来流体(5％HCl)后孔隙度分量和孔隙度分量累计值随T2弛豫时间变化关系曲线(见附图7)，需要指出的是，附图7内孔隙度分量曲线去除了核磁共振第一个峰包含的孔隙度分量值。
(3)不同尺寸孔隙对外来流体敏感性定量评价
①基于孔隙度分量累计值随T2弛豫时间变化关系曲线，引入页岩总孔隙度和不同尺寸孔隙度对外来流体敏感性定量评价指标，如公式(6)、(7)、(8)、(9)和(10)所示：
其中，
—核磁共振法计算的总孔隙度变化率，％；
—外来流体浸泡后总孔隙度分量累计值，％；
—饱和水后总孔隙度分量累计值，％；
—微孔孔隙度变化率，％；
—外来流体浸泡后微孔孔隙度分量累计值，％；
—饱和水后微孔孔隙度分量累计值，％；
—小孔孔隙度变化率，％；
—外来流体浸泡后小孔孔隙度分量累计值，％；
—饱和水后小孔孔隙度分量累计值，％；
—中孔孔隙度变化率，％；
—外来流体浸泡后中孔孔隙度分量累计值，％；
—饱和水后中孔孔隙度分量累计值，％；
—大孔孔隙度变化率，％；
—外来流体浸泡后大孔孔隙度分量累计值，％；
—饱和水后大孔孔隙度分量累计值，％。
需要指出的是，针对页岩储层推荐不同尺寸孔隙划分标准为：孔隙半径＜0.01μm为微孔、0.01-0.1μm为小孔、0.1-1.0μm为中孔、孔隙半径＞1.0μm为大孔。
②基于时空转换系数k，计算微孔、小孔、中孔和大孔孔隙半径划分对应的T2弛豫时间值，即10nm、100nm和1000nm对应的T2弛豫时间值T2,10nm、T2,100nm和T2,1000nm，
③基于T2,10nm、T2,100nm和T2,1000nm值，利用公式(11)-(18)计算饱和水后微孔、小孔、中孔和大孔孔隙度
和外来流体浸泡后微孔、小孔、中孔和大孔孔隙度
其中，
—饱和水后微孔T2谱核磁信号累计值，PU；
—饱和水后小孔T2谱核磁信号累计值，PU；
—饱和水后中孔T2谱核磁信号累计值，PU；
—外来流体浸泡后微孔T2谱核磁信号累计值，PU；
—外来流体浸泡后小孔T2谱核磁信号累计值，PU；
—外来流体浸泡后中孔T2谱核磁信号累计值，PU；Sac,sd—外来流体浸泡后T2谱第一个峰对应的核磁信号累计量，PU；；Sac,s—外来流体浸泡后T2谱核磁信号累计量，PU。
饱和水后微孔孔隙度
小孔孔隙度
中孔孔隙度
大孔孔隙度
外来流体(5％HCl)浸泡后，微孔孔隙度
小孔孔隙度
中孔孔隙度
大孔孔隙度
以
为例详细论述外来流体(5％HCl)浸泡后微孔孔隙度计算过程，其它不同尺寸孔隙度计算过程与之相类似，具体如下：
④5#页岩岩样饱和水后总孔隙度分量累计值与饱和水孔隙度相同，因此
由附图7可知：
基于公式(6)计算可知
参考SY/T 5358-2010《储层敏感性流动实验评价方法》敏感性影响程度评价指标(见表1)，页岩总孔隙度对外来流体(5％HCl)敏感性为改善(弱)，页岩浸泡5％HCl后核磁共振法计算的总孔隙度增加0.69％。
基于公式(7)、(8)、(9)和(10)计算页岩岩样浸泡外来流体(5％HCl)前后不同尺寸孔隙度变化率，参考SY/T 5358-2010《储层敏感性流动实验评价方法》敏感性影响程度评价指标(见表1)，给出页岩储层总孔隙和不同尺寸孔隙对外来流体(5％HCl)敏感性定量评价结果(见表3)。
表3总孔隙和不同尺寸孔隙对外来流体(5％HCl)敏感性定量评价结果
三、页岩基质渗透率对外来流体敏感性定量评价
松辽盆地青山口组青一段页岩大量发育微纳米级孔隙，层理缝渗透率一般在0.001-0.1mD，基质渗透率一般在10-9-10-4mD，采用SY/T5358-2010《储层敏感性流动实验评价方法》标准中传统的稳态法评估页岩基质渗透率对外来流体敏感性不具有现实意义。具体体现在以下两个方面：(1)页岩基质渗透率极低，需要较高的驱替压差和较长的测试时间才能达到稳定流动状态；(2)测量时需要人为记录实验数据，误差较大。基于上述两个原因，本发明基于非稳态法(压力脉冲衰减)建立了页岩基质渗透率对外来流体敏感性定量评价方法。
页岩基质渗透率对外来流体敏感性定量评价包含以下步骤：
①样品制备。选取与孔隙度和层理缝渗透率对外来流体(5％HCl)敏感性评价同一取芯层段的页岩岩样(6#岩样)，将样品粉碎研磨至10/20目颗粒，取颗粒质量为35g；
②将6#颗粒岩样在温度110℃条件下烘干8小时；
③基于压力脉冲衰减原理，利用SMP-200页岩基质渗透率仪测定6#样品研磨后颗粒的基质渗透率Kmb；
④将颗粒样品放置于活塞容器内，利用抽真空预处理系统对其抽真空，时间为48小时，为保护活塞容器，在其上部放置筛网；
⑤将预先配制好的外来流体(5％HCl)放置于另一个活塞容器内，利用ISCO泵将外来流体(5％HCl)注入到盛放颗粒样品的活塞容器内浸泡48小时，浸泡压力与青一段地层压力相同；
⑥取出颗粒样品，将其在温度110℃条件下烘干8小时；
⑦利用SMP-200页岩基质渗透率仪测定外来流体浸泡后颗粒样品的基质渗透率Kma；
⑧基于公式(19)计算页岩颗粒样品浸泡外来流体前后基质渗透率变化率：
其中，RmK—页岩颗粒样品基质渗透率变化率，％；Kmb—页岩颗粒样品浸泡外来流体前基质渗透率，mD；Kma—页岩颗粒样品浸泡外来流体后基质渗透率，mD。
⑨参考SY/T 5358-2010《储层敏感性流动实验评价方法》敏感性影响程度评价指标(见表1)，给出页岩基质渗透率对外来流体(5％HCl)敏感性定量评价结果(见表4)。
表4页岩基质渗透率对外来流体(5％HCl)敏感性定量评价结果
| 样品类型 | K<sub>mb</sub>(mD) | K<sub>ma</sub>(mD) | R<sub>mK</sub>(％) | 敏感性评价结果 |
| 页岩基质 | 1.49×10<sup>-4</sup> | 3.68×10<sup>-4</sup> | 146.98 | 改善(强) |
综上可知：页岩总孔隙度对外来流体(5％HCl)敏感性评价结果为改善(弱)、层理缝渗透率评价结果为改善(弱)、基质渗透率评价结果为改善(强)，综合分析给出实验用页岩样品对外来流体(5％HCl)敏感性评价结果为改善(弱)。
Claims (3)
1.一种页岩油气储层对外来流体敏感性定量评价新方法，该方法包括以下步骤：
(一)、岩样制备
选择来自同一取芯层段的页岩岩样3块，并将其标号为1#、2#和3#，其中1#和2#岩样用于确定目标烘干温度和时间，3#岩样用于完成孔隙度和渗透率对外来流体敏感性评价实验，3#页岩岩样要求：岩样为标准柱塞样，直径为2.5cm且长度应大于或等于5cm，实验前需对3#岩样进行预处理，将岩样切割成3块小岩样，并将其编号为4#、5#和6#，4#岩样用于高压压汞测试、5#岩样用于开展页岩孔隙度和层理缝渗透率敏感性评价实验、6#岩样研磨成10/20目颗粒用于页岩基质渗透率敏感性评价实验；
(二)、页岩孔隙度和层理缝渗透率对外来流体敏感性定量评价
页岩孔隙度和层理缝渗透率对外来流体敏感性定量评价过程包含以下3个步骤：
A、确定页岩岩样目标烘干温度和时间
在开展页岩孔隙度和层理缝渗透率对外来流体敏感性评价实验之前，需清除岩样中存有的烃类、水等可动流体，确保岩样清洁；
(1)选择来自同一取芯层段的页岩岩样2块，并将其标号为1#和2#；
(2)设定10个不同烘干温度Tt1、Tt2、Tt3、Tt4、Tt5、Tt6、Tt7、Tt8、Tt9、Tt10，温度按照20℃依次递增；
(3)在设定温度条件下烘干1#页岩岩样t0小时，测量1#页岩岩样在不同烘干温度条件下对应的岩样质量m1、m2、m3、m4、m5、m6、m7、m8、m9、m10并开展烘干后页岩二维核磁共振测试；
(4)绘制1#页岩岩样质量随烘干温度变化曲线；当页岩岩样内可动流体完全清除时，岩样质量随烘干温度变化曲线出现拐点，拐点对应的温度为目标烘干温度T0；
(5)设定10个不同烘干时间t1、t2、t3、t4、t5、t6、t7、t8、t9、t10，时间按照2小时依次递增；在确保实验准确度条件下，为减少测试工作量，可将t5或者t6时间设定与t0相同；
(6)在目标烘干温度T0条件下按照设定时间烘干2#页岩岩样，测量页岩岩样在不同烘干时间条件下对应的岩样质量m′1、m′2、m′3、m′4、m′5、m′6、m′7、m′8、m′9、m′10并开展烘干后页岩二维核磁共振测试；
(7)绘制2#页岩岩样质量随烘干时间变化曲线；当页岩岩样内可动流体完全清除时，岩样质量不随烘干时间变化，此时对应的时间为目标烘干时间t0；
(8)分析大于目标烘干时间t0之后的二维核磁共振T1和T2图谱，若小、中、大孔内的核磁信号量不发生改变，则可确定页岩岩样最终的目标烘干温度和烘干时间；若二维核磁共振T1和T2图谱中小、中、大孔内的核磁信号量发生改变，则需重新筛选岩样按照步骤(1)-(7)开展实验，再次确定页岩岩样目标烘干温度T0和时间t0；
B、开展页岩岩样高压压汞+外来流体浸泡+核磁共振联测实验
(1)页岩岩样高压压汞测试
将制备好的4#页岩岩样按照确定的目标烘干温度T0和时间t0进行烘干处理，将烘干后的岩样开展高压压汞测试，在半对数坐标条件下绘制孔隙分布频率随孔隙半径变化曲线，确定页岩岩样孔隙分布特征；
(2)页岩岩样与外来流体浸泡+核磁共振联测实验；
C、实验数据处理与分析
1)页岩层理缝渗透率对外来流体敏感性定量评价
(1)引入层理缝渗透率对外来流体敏感性定量评价指标，如公式(1)所示：
其中，RK—层理缝渗透率变化率，％；Kb—页岩干样层理缝渗透率，mD；Ka—页岩浸泡外来流体后层理缝渗透率，mD；
(2)按照公式(1)计算页岩浸泡外来流体前后层理缝渗透率变化率；
(3)参考SY/T 5358-2010《储层敏感性流动实验评价方法》敏感性影响程度评价指标，定量评价页岩层理缝渗透率对外来流体敏感性；
2)页岩总孔隙度对外来流体敏感性定量评价
(1)引入页岩总孔隙度对外来流体敏感性定量评价指标，如公式(2)所示：
其中，Rt—覆压孔渗仪法测量的总孔隙度变化率，％；φb—页岩干样孔隙度，％；φa—页岩浸泡外来流体后孔隙度，％；
(2)按照公式(2)计算页岩浸泡外来流体前后孔隙度变化率；
(3)参考SY/T 5358-2010《储层敏感性流动实验评价方法》敏感性影响程度评价指标，定量评价页岩总孔隙度对外来流体敏感性；
3)页岩不同尺寸孔隙对外来流体敏感性定量评价
(1)核磁共振弛豫时间T2值和孔隙尺寸时空转换系数计算；
(2)核磁共振T2谱信号量与页岩饱和水孔隙度标定关系建立
①基于重量法计算5#页岩岩样饱和水孔隙度，计算公式如下：
其中，φw—饱和水孔隙度，％；mw—饱和水后岩样重量，g；m0—烘干后岩样重量，g；ρw—水密度，g/cm3；d—岩样直径，cm；L—岩样长度，cm；
②基于核磁共振测试T2谱信号分量，按照公式(4)将页岩饱和水后T2谱核磁信号分量转化为孔隙度分量：
其中，
—饱和水后孔隙度分量，％；Sw—饱和水后T2谱核磁信号分量，PU；Sac,w—饱和水后T2谱核磁信号累计量，PU；Sac,d—饱和水T2谱第一个峰对应的核磁信号累计量，PU；
③基于步骤②饱和水后T2谱核磁信号分量与孔隙度分量转换原理，按照公式(5)将页岩浸泡外来流体后T2谱核磁信号分量转化为孔隙度分量：
其中，
—外来流体浸泡后孔隙度分量，％；Ss—外来流体浸泡后T2谱核磁信号分量，PU；
④基于步骤②和③得到的孔隙度分量，绘制页岩饱和水和浸泡外来流体后孔隙度分量和孔隙度分量累计值随T2弛豫时间变化关系曲线；需要指出的是，上述曲线需去除T2谱内第一个峰对应的核磁信号累计量；
(3)不同尺寸孔隙对外来流体敏感性定量评价
①基于孔隙度分量累计值随T2弛豫时间变化关系曲线，引入页岩总孔隙度和不同尺寸孔隙度对外来流体敏感性定量评价指标，如公式(6)、(7)、(8)、(9)和(10)所示：
其中，
—核磁共振法计算的总孔隙度变化率，％；
—外来流体浸泡后总孔隙度分量累计值，％；
—饱和水后总孔隙度分量累计值，％；
—微孔孔隙度变化率，％；
—外来流体浸泡后微孔孔隙度分量累计值，％；
—饱和水后微孔孔隙度分量累计值，％；
—小孔孔隙度变化率，％；
—外来流体浸泡后小孔孔隙度分量累计值，％；
—饱和水后小孔孔隙度分量累计值，％；
—中孔孔隙度变化率，％；
—外来流体浸泡后中孔孔隙度分量累计值，％；
—饱和水后中孔孔隙度分量累计值，％；
—大孔孔隙度变化率，％；
—外来流体浸泡后大孔孔隙度分量累计值，％；
—饱和水后大孔孔隙度分量累计值，％；
需要指出的是，针对页岩储层推荐不同尺寸孔隙划分标准为：孔隙半径＜0.01μm为微孔、0.01-0.1μm为小孔、0.1-1.0μm为中孔、孔隙半径＞1.0μm为大孔；
②基于时空转换系数k，计算微孔、小孔、中孔和大孔孔隙半径划分对应的T2弛豫时间值，即10nm、100nm和1000nm对应的T2弛豫时间值T2,10nm、T2,100nm和T2,1000nm；
③基于T2,10nm、T2,100nm和T2,1000nm值，利用公式(11)-(18)计算饱和水后微孔、小孔、中孔和大孔孔隙度
和外来流体浸泡后微孔、小孔、中孔和大孔孔隙度
其中，
—饱和水后微孔T2谱核磁信号累计值，PU；
—饱和水后小孔T2谱核磁信号累计值，PU；
—饱和水后中孔T2谱核磁信号累计值，PU；
—外来流体浸泡后微孔T2谱核磁信号累计值，PU；
—外来流体浸泡后小孔T2谱核磁信号累计值，PU；
—外来流体浸泡后中孔T2谱核磁信号累计值，PU；Sac,sd—外来流体浸泡后T2谱第一个峰对应的核磁信号累计量，PU；Sac,s—外来流体浸泡后T2谱核磁信号累计量，PU；
④基于公式(6)、(7)、(8)、(9)和(10)计算页岩岩样浸泡外来流体前后总孔隙度和不同尺寸孔隙度变化率，参考SY/T5358-2010《储层敏感性流动实验评价方法》敏感性影响程度评价指标，给出页岩储层总孔隙度和不同尺寸孔隙对外来流体敏感性定量评价结果；
(三)、页岩基质渗透率对外来流体敏感性定量评价
页岩基质渗透率对外来流体敏感性定量评价包含以下步骤：
①样品制备：选取与孔隙度和层理缝渗透率对外来流体敏感性评价同一取芯层段的页岩岩样，将样品粉碎研磨至10/20目颗粒，质量大于30g；
②在目标烘干温度T0和时间t0条件下对颗粒样品进行烘干处理；
③基于压力脉冲衰减原理，利用SMP-200页岩基质渗透率仪测定6#样品研磨后颗粒的基质渗透率Kmb；
④将颗粒样品放置于活塞容器内，利用抽真空预处理系统对其抽真空，时间为48小时，为保护活塞容器，在其上部放置筛网；
⑤将预先配制好的外来流体放置于另一个活塞容器内，利用ISCO泵将外来流体注入到盛放颗粒样品的活塞容器内浸泡48小时，浸泡压力与地层压力相同；
⑥取出颗粒样品，将其在目标烘干温度T0和时间t0条件下烘干；
⑦利用SMP-200页岩基质渗透率仪测定外来流体浸泡后颗粒样品的基质渗透率Kma；
⑧基于公式(19)计算页岩颗粒样品浸泡外来流体前后基质渗透率变化率：
其中，RmK—页岩颗粒样品基质渗透率变化率，％；Kmb—页岩颗粒样品浸泡外来流体前(烘干后)基质渗透率，mD；Kma—页岩颗粒样品浸泡外来流体后基质渗透率，mD；
⑨参考SY/T 5358-2010《储层敏感性流动实验评价方法》敏感性影响程度评价指标，给出页岩基质渗透率对外来流体敏感性定量评价结果；
基于页岩浸泡外来流体前后总孔隙度、不同尺寸孔隙度、层理缝渗透率和基质渗透率变化规律，综合分析给出页岩油气储层对外来流体敏感性定量评价结果。
2.根据权利要求1所述的页岩油气储层对外来流体敏感性定量评价新方法，其特征在于：步骤B中第(2)步页岩岩样与外来流体浸泡+核磁共振联测实验包括以下步骤：
①将5#岩样按照确定的目标烘干温度T0和时间t0进行烘干处理并记录其长度L、直径d和质量m0；
②利用VINCI覆压孔渗仪测量烘干后岩样的孔隙度φb和层理缝渗透率Kb；
③将测量孔渗后的标准岩样放置于岩心夹持器内，利用手摇泵加载2MPa围压，之后利用抽真空预处理系统对岩样进行抽真空，时间为48小时；
④制备实验用蒸馏水并将其放置于活塞容器内，利用ISCO泵将活塞容器内蒸馏水以恒压方式注入岩样内，饱和过程中注入压力与围压呈阶梯方式递增且保持围压与注压差值为2MPa恒定，当注压达到地层压力时，停止饱和过程；岩样饱和过程在恒温箱内完成，温度与地层温度保持一致，总饱和时间不低于48小时；
⑤将饱和水后岩样从岩心夹持器内取出并记录其质量mw，测试岩样饱和水后核磁共振T2谱，为减小实验误差，需连续测量3次；
⑥将步骤⑤饱和水后岩样在目标烘干温度T0和时间t0进行烘干处理；
⑦将烘干后岩样放置于岩心夹持器内，利用手摇泵加载2MPa围压，之后利用抽真空预处理系统对岩样进行抽真空，时间为48小时；
⑧利用ISCO泵将活塞容器内外来流体以恒压方式注入岩样内，注入压力与围压呈阶梯方式递增且保持围压与注压差值为2MPa恒定，当注压达到地层压力时，停止饱和过程，总饱和时间不低于48小时；岩样饱和与浸泡外来流体过程在恒温箱内完成，温度与地层温度保持一致，浸泡时间为48小时；
⑨开展岩样浸泡外来流体后核磁共振T2谱测试，为减小实验误差，需连续测量3次；
⑩将步骤⑨浸泡外来流体后的页岩岩样在目标烘干温度T0和时间t0进行烘干处理；
3.根据权利要求1所述的页岩油气储层对外来流体敏感性定量评价新方法，其特征在于：步骤C中第3)步中的第(1)步核磁共振弛豫时间T2值和孔隙尺寸时空转换系数计算中包括以下步骤：
①提取5#岩样饱和水后核磁共振T2弛豫时间与核磁信号量数据，在半对数坐标条件下绘制核磁信号量随T2弛豫时间变化曲线；
②提取高压压汞测试结果中的孔隙半径和孔隙分布频率数据，在半对数坐标条件下绘制孔隙分布频率随孔隙半径变化曲线；
③将步骤①和②曲线数据整合在同一坐标系内，建立核磁共振与高压压汞时空转换曲线，其中X轴底部坐标为孔隙半径，Y轴主坐标为孔隙分布频率；X轴顶部坐标为T2弛豫时间，Y轴副坐标为核磁信号量；
④页岩岩样饱和水后核磁共振T2谱第一个峰代表有机质信号显示，因此以第二个峰对应的弛豫时间开始进行时空转换系数计算，当核磁信号峰值与高压压汞孔隙分布频率峰值逐一对应时，记录此时的核磁共振弛豫时间T2i和孔隙半径ri值；基于以上数据，即可计算页岩岩样核磁共振弛豫时间T2值和孔隙尺寸时空转换系数k＝(k1+k2+…+kn)/n,其中n代表核磁信号与高压压汞孔隙分布频率对应峰的数量。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210845363.XA CN115372398B (zh) | 2022-07-19 | 2022-07-19 | 页岩油气储层对外来流体敏感性定量评价方法 |
| US18/353,553 US20240027379A1 (en) | 2022-07-19 | 2023-07-17 | Method for quantitative evaluation on sensitivity of shale oil and gas reservoir to injected fluids |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210845363.XA CN115372398B (zh) | 2022-07-19 | 2022-07-19 | 页岩油气储层对外来流体敏感性定量评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115372398A true CN115372398A (zh) | 2022-11-22 |
| CN115372398B CN115372398B (zh) | 2023-03-07 |
Family
ID=84061331
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210845363.XA Active CN115372398B (zh) | 2022-07-19 | 2022-07-19 | 页岩油气储层对外来流体敏感性定量评价方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US20240027379A1 (zh) |
| CN (1) | CN115372398B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116539815A (zh) * | 2023-06-07 | 2023-08-04 | 四川省科源工程技术测试中心有限责任公司 | 适用于油气藏工作液评价与优选的装置与方法 |
| CN116539655A (zh) * | 2023-07-06 | 2023-08-04 | 山东石油化工学院 | 一种基于核磁共振技术评价致密砂岩储层水敏性的方法 |
| CN117110349A (zh) * | 2023-06-28 | 2023-11-24 | 东北石油大学 | 一种页岩油储层岩石多尺度储集空间润湿性定量评价方法 |
| CN117147318A (zh) * | 2023-10-12 | 2023-12-01 | 东北石油大学三亚海洋油气研究院 | 一种页岩储层多尺度孔隙应力敏感性评价方法 |
| CN117760935A (zh) * | 2024-02-19 | 2024-03-26 | 东北石油大学三亚海洋油气研究院 | 页岩油储层不同尺度储集空间吞吐液量确定方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104990851A (zh) * | 2015-06-23 | 2015-10-21 | 西南石油大学 | 一种新的页岩敏感性实验研究方法 |
| US20180348396A1 (en) * | 2016-01-25 | 2018-12-06 | Halliburton Energy Services, Inc. | Determining Permeability In Subsurface Anisotropic Formations |
-
2022
- 2022-07-19 CN CN202210845363.XA patent/CN115372398B/zh active Active
-
2023
- 2023-07-17 US US18/353,553 patent/US20240027379A1/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104990851A (zh) * | 2015-06-23 | 2015-10-21 | 西南石油大学 | 一种新的页岩敏感性实验研究方法 |
| US20180348396A1 (en) * | 2016-01-25 | 2018-12-06 | Halliburton Energy Services, Inc. | Determining Permeability In Subsurface Anisotropic Formations |
Non-Patent Citations (4)
| Title |
|---|
| JIANGUANG WEI ET AL.: "Experimental study on oil recovery mechanism of CO2 associated enhancing oil recovery methods in low permeability reservoirs", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》 * |
| 张顺等: "东营凹陷页岩油储层孔隙演化", 《石油学报》 * |
| 朱维耀等: "层理缝对页岩渗透率的影响及表征", 《特种油气藏》 * |
| 罗超等: "超低渗岩心渗透率应力敏感性评价指标对比", 《石油地质与工程》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116539815A (zh) * | 2023-06-07 | 2023-08-04 | 四川省科源工程技术测试中心有限责任公司 | 适用于油气藏工作液评价与优选的装置与方法 |
| CN116539815B (zh) * | 2023-06-07 | 2024-03-19 | 四川省科源工程技术测试中心有限责任公司 | 适用于油气藏工作液评价与优选的装置与方法 |
| CN117110349A (zh) * | 2023-06-28 | 2023-11-24 | 东北石油大学 | 一种页岩油储层岩石多尺度储集空间润湿性定量评价方法 |
| CN117110349B (zh) * | 2023-06-28 | 2024-05-28 | 东北石油大学 | 一种页岩油储层岩石多尺度储集空间润湿性定量评价方法 |
| CN116539655A (zh) * | 2023-07-06 | 2023-08-04 | 山东石油化工学院 | 一种基于核磁共振技术评价致密砂岩储层水敏性的方法 |
| CN116539655B (zh) * | 2023-07-06 | 2023-10-10 | 山东石油化工学院 | 一种基于核磁共振技术评价致密砂岩储层水敏性的方法 |
| CN117147318A (zh) * | 2023-10-12 | 2023-12-01 | 东北石油大学三亚海洋油气研究院 | 一种页岩储层多尺度孔隙应力敏感性评价方法 |
| CN117147318B (zh) * | 2023-10-12 | 2024-05-07 | 东北石油大学三亚海洋油气研究院 | 一种页岩储层多尺度孔隙应力敏感性评价方法 |
| CN117760935A (zh) * | 2024-02-19 | 2024-03-26 | 东北石油大学三亚海洋油气研究院 | 页岩油储层不同尺度储集空间吞吐液量确定方法 |
| CN117760935B (zh) * | 2024-02-19 | 2024-05-03 | 东北石油大学三亚海洋油气研究院 | 页岩油储层不同尺度储集空间吞吐液量确定方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| US20240027379A1 (en) | 2024-01-25 |
| CN115372398B (zh) | 2023-03-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115372398B (zh) | 页岩油气储层对外来流体敏感性定量评价方法 | |
| CN108414560B (zh) | 一种核磁-驱替联用装置评价致密油充注过程的方法 | |
| CN106501144B (zh) | 一种基于核磁共振双截止值的致密砂岩渗透率计算方法 | |
| CN106525688B (zh) | 一种饱和页岩孔隙流体分离及饱和度计算的实验方法 | |
| CN103926186B (zh) | 一种定量评价注水对孔喉分布影响的方法 | |
| CN110595953B (zh) | 一种页岩混合润湿性的实验测试装置及方法 | |
| CN109884109B (zh) | 一种页岩孔隙水吸附量、游离量及其分布的评价方法 | |
| CN104075959B (zh) | 泥页岩油水饱和度精确测试方法 | |
| CN106153662A (zh) | 岩心应力敏感性的测量方法 | |
| CN104697915A (zh) | 一种页岩微观孔隙大小及流体分布的分析方法 | |
| CN110296931B (zh) | 一种致密砂岩油水相对渗透率信息的表征方法及系统 | |
| CN113075102B (zh) | 一种建立多孔介质自发渗吸量与时间关系数学模型的方法 | |
| CN106897531A (zh) | 一种低渗透石灰岩储层渗透率的定量评价方法 | |
| CN112946005A (zh) | 一种页岩微裂缝评价方法及其应用 | |
| CN112129802A (zh) | 一种水化页岩不同尺度孔隙体积增量的定量分析方法 | |
| RU2542998C1 (ru) | Способ лабораторного определения анизотропии абсолютной газопроницаемости на полноразмерном керне | |
| CN112487620B (zh) | 一种页岩油可动资源量的评价方法 | |
| Yin et al. | Development of an NMR workflow for determining nano-petrophysical properties of marine and lacustrine mudrocks | |
| CN116539655B (zh) | 一种基于核磁共振技术评价致密砂岩储层水敏性的方法 | |
| CN115479964A (zh) | 一种致密油储层压后返排的物理模拟方法 | |
| CN114594033B (zh) | 一种评价油气储层水敏、水锁和固相伤害的方法 | |
| CN115876659A (zh) | 一种获取页岩孔隙和裂缝体积的方法及系统 | |
| CN114076777A (zh) | 一种页岩含气饱和度的获取方法 | |
| CN117760935B (zh) | 页岩油储层不同尺度储集空间吞吐液量确定方法 | |
| CN116660306B (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 |