CN111220639B - 基于核磁共振的气驱水时岩心含气饱和度测定方法及装置 - Google Patents
基于核磁共振的气驱水时岩心含气饱和度测定方法及装置 Download PDFInfo
- Publication number
- CN111220639B CN111220639B CN201811416441.4A CN201811416441A CN111220639B CN 111220639 B CN111220639 B CN 111220639B CN 201811416441 A CN201811416441 A CN 201811416441A CN 111220639 B CN111220639 B CN 111220639B
- Authority
- CN
- China
- Prior art keywords
- core
- nuclear magnetic
- water
- gas
- 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
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 102
- 238000000034 method Methods 0.000 title claims abstract description 72
- 239000011435 rock Substances 0.000 title abstract description 21
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 257
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 84
- 238000001228 spectrum Methods 0.000 claims abstract description 58
- 238000001514 detection method Methods 0.000 claims abstract description 34
- 230000008569 process Effects 0.000 claims description 28
- 239000012530 fluid Substances 0.000 claims description 19
- 229920001971 elastomer Polymers 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 9
- -1 polytetrafluoroethylene Polymers 0.000 claims description 9
- 229920001343 polytetrafluoroethylene Polymers 0.000 claims description 9
- 239000004810 polytetrafluoroethylene Substances 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 abstract description 35
- 238000005303 weighing Methods 0.000 description 11
- 238000013461 design Methods 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 7
- 230000009286 beneficial effect Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 230000002708 enhancing effect Effects 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 4
- 238000005065 mining Methods 0.000 description 4
- 239000003208 petroleum Substances 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000001035 drying Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000006073 displacement reaction Methods 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
- 239000011148 porous material Substances 0.000 description 2
- 238000000685 Carr-Purcell-Meiboom-Gill pulse sequence Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
Abstract
本发明提供了一种基于核磁共振的气驱水时岩心含气饱和度测定方法及装置，其中，该方法包括：采集低场核磁共振方式下饱和水状态的待测岩心T2谱；根据饱和水状态的待测岩心T2谱，确定饱和水状态的待测岩心中水核磁信号峰面积；采集低场核磁共振方式下气驱水时的待测岩心T2谱；根据气驱水时的待测岩心T2谱，确定气驱水时的待测岩心中水核磁信号峰面积；根据饱和水状态的待测岩心中水核磁信号峰面积和气驱水时的待测岩心中水核磁信号峰面积，确定气驱水时岩心含气饱和度。上述技术方案实现了气驱水时岩心含气饱和度的在线测定，提高了含气饱和度的测定精度，同时还可以实现气驱水过程中岩心含气饱和度的连续检测。
Description
技术领域
本发明涉及油藏勘探技术领域，特别涉及一种基于核磁共振的气驱水时岩心含气饱和度测定方法及装置。
背景技术
含气饱和度是储层评价的一个重要参数，准确地测定气驱水过程中岩心含气饱和度是成藏过程物理模拟的重要任务之一。目前常用的气驱水过程中岩心含气饱和度测定是结合称重测量来完成的，有很大的局限性。主要体现在以下方面：岩心取出称重过程中会有部分水因暴露于外界环境中而损失，导致含水量测量值偏小，含气饱和度偏大，测试结果不精确；只能进行气驱水过程中设计点的含气饱和度测定，无法实现气驱水过程中含气饱和度的连续检测。
发明内容
本发明实施例提供了一种基于核磁共振的气驱水时岩心含气饱和度测定方法，用以实现气驱水时岩心含气饱和度的在线测定，提高含气饱和度的测定精度，同时实现气驱水过程中岩心含气饱和度的连续检测，该方法包括：
采集低场核磁共振方式下饱和水状态的待测岩心T2谱；
根据所述饱和水状态的待测岩心T2谱，确定饱和水状态的待测岩心中水核磁信号峰面积；
采集低场核磁共振方式下气驱水时的待测岩心T2谱；
根据所述气驱水时的待测岩心T2谱，确定气驱水时的待测岩心中水核磁信号峰面积；
根据饱和水状态的待测岩心中水核磁信号峰面积和气驱水时的待测岩心中水核磁信号峰面积，确定气驱水时岩心含气饱和度；
其中，根据饱和水状态的待测岩心中水核磁信号峰面积和气驱水时的待测岩心中水核磁信号峰面积，确定气驱水时岩心含气饱和度，包括：
根据饱和水状态的待测岩心中水核磁信号峰面积和气驱水时的待测岩心中水核磁信号峰面积，确定气驱水时待测岩心含水饱和度；
根据所述气驱水时待测岩心含水饱和度，确定气驱水时待测岩心含气饱和度。
本发明实施例还提供了一种基于核磁共振的气驱水时岩心含气饱和度测定装置，用以实现气驱水时岩心含气饱和度的在线测定，提高含气饱和度的测定精度，同时实现气驱水过程中岩心含气饱和度的连续检测，该装置包括：
第一采集单元，用于采集低场核磁共振方式下饱和水状态的待测岩心T2谱；
饱和水状态峰面积确定单元，用于根据所述饱和水状态的待测岩心T2谱，确定饱和水状态的待测岩心中水核磁信号峰面积；
第二采集单元，用于采集低场核磁共振方式下气驱水时的待测岩心T2谱；
气驱水时峰面积确定单元，用于根据所述气驱水时的待测岩心T2谱，确定气驱水时的待测岩心中水核磁信号峰面积；
含气饱和度确定单元，用于根据饱和水状态的待测岩心中水核磁信号峰面积和气驱水时的待测岩心中水核磁信号峰面积，确定气驱水时岩心含气饱和度；
其中，所述含气饱和度确定单元具体用于：
根据饱和水状态的待测岩心中水核磁信号峰面积和气驱水时的待测岩心中水核磁信号峰面积，确定气驱水时待测岩心含水饱和度；
根据所述气驱水时待测岩心含水饱和度，确定气驱水时待测岩心含气饱和度。
本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述基于核磁共振的气驱水时岩心含气饱和度测定方法。
本发明实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有执行基于核磁共振的气驱水时岩心含气饱和度测定方法的计算机程序。
本发明实施例提供的技术方案通过：采集低场核磁共振方式下饱和水状态的待测岩心T2谱；根据饱和水状态的待测岩心T2谱，确定饱和水状态的待测岩心中水核磁信号峰面积；采集低场核磁共振方式下气驱水时的待测岩心T2谱；根据气驱水时的待测岩心T2谱，确定气驱水时的待测岩心中水核磁信号峰面积；根据饱和水状态的待测岩心中水核磁信号峰面积和气驱水时的待测岩心中水核磁信号峰面积，确定气驱水时岩心含气饱和度，实现了采用低场核磁共振的方式检测气驱水时岩心中水的核磁信号，实现气驱水时岩心含气饱和度的在线测定，可避免称重法在测定含水量时因样品暴露在环境中的损失量，使含气饱和度的测定更加准确，同时还可以实现气驱水过程中岩心含气饱和度的连续检测。
本发明实施例提供的技术方案适用于石油、地质、矿业领域在实验室进行气驱水实验时岩心含气饱和度的测定，有助于增强成藏物理模拟实验的效果，可以服务于储层评价和油气成藏机理研究。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1是本发明实施例中基于核磁共振的气驱水时岩心含气饱和度测定方法的流程示意图；
图2是本发明实施例中低场核磁共振方式检测到的气驱水过程中岩心T2谱的示意图；
图3是利用本发明实施例中基于核磁共振的气驱水时岩心含气饱和度测定方法测试结果与现有称重法测试结果的对比示意图；
图4是本发明实施例中基于核磁共振的气驱水时岩心含气饱和度测定装置的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并不作为对本发明的限定。
由于考虑到现有称重法在测定含水量时因样品暴露在环境中的损失量，使含气饱和度的测定不精确，同时无法实现气驱水过程中岩心含气饱和度的连续检测的技术问题，发明人提出了一种基于核磁共振的气驱水时岩心含气饱和度在线测定方案，该方案采用低场核磁共振的方式检测气驱水时岩心中水的核磁信号强度，能够实现气驱水时岩心含气饱和度的在线测定，可避免称重法在测定含水量时因样品暴露在环境中的损失量，使含气饱和度的测定更加准确，同时还可以实现气驱水过程中岩心含气饱和度的连续检测。从而本发明的方案更适用于石油、地质、矿业领域在实验室进行气驱水实验时岩心含气饱和度的测定，有助于增强成藏物理模拟实验的效果，可以服务于储层评价和油气成藏机理研究。下面对该基于核磁共振的气驱水时岩心含气饱和度测定方案进行详细介绍如下。
图1是本发明实施例中基于核磁共振的气驱水时岩心含气饱和度测定方法的流程示意图，如图1所示，该方法包括如下步骤：
步骤101：采集低场核磁共振方式下饱和水状态的待测岩心T2谱；
步骤102：根据所述饱和水状态的待测岩心T2谱，确定饱和水状态的待测岩心中水核磁信号峰面积；
步骤103：采集低场核磁共振方式下气驱水时的待测岩心T2谱；
步骤104：根据所述气驱水时的待测岩心T2谱，确定气驱水时的待测岩心中水核磁信号峰面积；
步骤105：根据饱和水状态的待测岩心中水核磁信号峰面积和气驱水时的待测岩心中水核磁信号峰面积，确定气驱水时岩心含气饱和度；
其中，根据饱和水状态的待测岩心中水核磁信号峰面积和气驱水时的待测岩心中水核磁信号峰面积，确定气驱水时岩心含气饱和度，包括：
根据饱和水状态的待测岩心中水核磁信号峰面积和气驱水时的待测岩心中水核磁信号峰面积，确定气驱水时待测岩心含水饱和度；
根据所述气驱水时待测岩心含水饱和度，确定气驱水时待测岩心含气饱和度。
本发明实施例提供的技术方案通过：采集低场核磁共振方式下饱和水状态的待测岩心T2谱；根据饱和水状态的待测岩心T2谱，确定饱和水状态的待测岩心中水核磁信号峰面积；采集低场核磁共振方式下气驱水时的待测岩心T2谱；根据气驱水时的待测岩心T2谱，确定气驱水时的待测岩心中水核磁信号峰面积；根据饱和水状态的待测岩心中水核磁信号峰面积和气驱水时的待测岩心中水核磁信号峰面积，确定气驱水时岩心含气饱和度，实现了采用低场核磁的方式检测气驱水时岩心中水的核磁信号，实现气驱水时岩心含气饱和度的在线测定，可避免称重法在测定含水量时因样品暴露在环境中的损失量，使含气饱和度的测定更加准确，同时还可以实现气驱水过程中岩心含气饱和度的连续检测。
本发明实施例提供的技术方案适用于石油、地质、矿业领域在实验室进行气驱水实验时岩心含气饱和度的测定，有助于增强成藏物理模拟实验的效果，可以服务于储层评价和油气成藏机理研究。
下面对本发明实施例提供的基于核磁共振的气驱水时岩心含气饱和度测定方法的各个步骤进行详细介绍如下。
一、首先，介绍步骤101。
具体实施时，在进行采集T2谱之前，首先对待测岩心依次进行烘干、抽真空、加压饱和等操作，使岩心孔隙中完全饱和水。然后，将岩心装于无磁夹持器中，并置于低场核磁共振设备检测范围内。
具体实施时，上述烘干、抽真空、加压饱和等操作可以参考国家标准GB/T29172-2012《岩心分析方法》进行，提高了含气饱和度测定的准确度。
在一实施例中，所述饱和水状态的待测岩心可以被至于无磁夹持器内部；所述无磁夹持器内部的胶套可以为聚四氟乙烯材料制胶套；
具体实施时，无磁夹持器以及无磁夹持器内部胶套由聚四氟乙烯材料构成，也提高了含气饱和度测定的准确度。
在一实施例中，采集低场核磁共振方式下饱和水状态的待测岩心T2谱，可以包括：
配置低场核磁共振方式下的低场核磁共振检测参数与待测岩心物性特征及流体性质相匹配；
在低场核磁共振检测参数与待测岩心物性特征及流体性质相匹配后，采集低场核磁共振方式下饱和水状态的待测岩心T2谱。
具体实施时，低场核磁共振检测参数的设置与岩心物性特征及流体性质相匹配，提高了含气饱和度测定的准确度。
二、其次，介绍上述步骤102。
在一实施例中，根据所述饱和水状态的待测岩心T2谱(可以如图2所示)，确定饱和水状态的待测岩心中水核磁信号峰面积，包括按照如下公式确定饱和水状态的待测岩心中水核磁信号峰面积：
其中，a代表饱和水状态的待测岩心中水核磁信号峰面积，T2代表饱和水状态的横向弛豫时间，t为饱和水状态的横向弛豫时间最大值，A(T2)代表饱和水状态的横向弛豫时间为T2时刻的水信号强度。
具体实施时，采用低场核磁共振方式对岩心进行检测，采集饱和水状态岩心的T2谱，计算岩心中水核磁信号峰面积其中，a代表饱和水状态的待测岩心中水核磁信号峰面积，T2代表饱和水状态的横向弛豫时间，t为饱和水状态的横向弛豫时间最大值，A(T2)代表饱和水状态的横向弛豫时间为T2时刻的水信号强度。采用低场核磁共振的方式，检测饱和水状态岩心中水的核磁信号强度，能够实现气驱水时岩心含气饱和度的在线测定，可避免称重法在测定含水量时因样品暴露在环境中的损失量，使含气饱和度的测定更加准确，同时还可以实现气驱水过程中岩心含气饱和度的连续检测。
三、接着，介绍上述步骤103。
在一个实施例中，所述气驱水时的待测岩心可以被至于无磁夹持器内部；所述无磁夹持器内部的胶套可以为聚四氟乙烯材料制胶套。
具体实施时，无磁夹持器以及无磁夹持器内部胶套由聚四氟乙烯材料构成，也提高了含气饱和度测定的准确度。
在一个实施例中，采集低场核磁共振方式下气驱水时的待测岩心T2谱，可以包括：
配置低场核磁共振方式下的低场核磁共振检测参数与待测岩心物性特征及流体性质相匹配；
在低场核磁共振检测参数与待测岩心物性特征及流体性质相匹配后，采集低场核磁共振方式下气驱水时的待测岩心T2谱。
具体实施时，低场核磁共振检测参数的设置与岩心物性特征及流体性质相匹配，提高了含气饱和度测定的准确度。
四、接着，介绍上述步骤104。
在一实施例中，根据所述气驱水时的待测岩心T2谱(可以如图2所示)，确定气驱水时的待测岩心中水核磁信号峰面积，包括按照如下公式确定气驱水时的待测岩心中水核磁信号峰面积：
其中，b代表气驱水时的待测岩心中水核磁信号峰面积，T2'代表气驱水时的横向弛豫时间，t'为气驱水时的横向弛豫时间最大值，B(T2')代表气驱水时的横向弛豫时间为T2'时刻的水信号强度。
具体实施时，采用低场核磁共振方式对气驱水时的岩心进行检测，采集驱替过程中岩心的T2谱，计算气驱水时岩心中水核磁信号峰面积其中，b代表气驱水时的待测岩心中水核磁信号峰面积，T2'代表气驱水时的横向弛豫时间，t'为气驱水时的横向弛豫时间最大值，B(T2')代表气驱水时的横向弛豫时间为T2'时刻的水信号强度。采用低场核磁共振的方式，检测气驱水时岩心中水的核磁信号强度，能够实现气驱水时岩心含气饱和度的在线测定，可避免称重法在测定含水量时因样品暴露在环境中的损失量，使含气饱和度的测定更加准确，同时还可以实现气驱水过程中岩心含气饱和度的连续检测。
具体实施时，可以进行不同设计点的气驱水实验，即针对不同注入压力点的气驱水实验，可以得到不同压力下的气驱水时岩心含气饱和度。
五、接着，介绍上述步骤105。
具体实施时，计算气驱水时岩心含水饱和度那么，气驱水时岩心含气饱和度/>
下面举一实例以说明本发明如何实施。
按照以下操作步骤进行基于核磁共振的气驱水时岩心含气饱和度在线测定：
步骤(1)：选取砂岩样品1(样品直径2.5cm，长度5.7cm，孔隙度11.5％，渗透率0.437mD)，参考国家标准GB/T29172-2012《岩心分析方法》对岩心依次进行烘干、抽真空、加压饱和，使岩心孔隙中完全饱和水；
步骤(2)：将岩心装于内部胶套由聚四氟乙烯材料构成的无磁夹持器中置于低场核磁共振设备检测范围内；
步骤(3)：采用低场核磁共振方式对岩心进行检测，其中核磁设备主频12.5MHz，低场核磁共振设备检测参数的设置与岩心物性特征及流体性质相匹配，核磁检测采用CPMG序列，参数设置：SF＝12MHZ；O1＝0.578MHZ；TW＝40000ms；NS＝64；TE＝0.300ms；NECH＝15000，采集饱和水状态岩心的T2谱，通过公式(1)计算岩心中水核磁信号峰面积a(16158)；
步骤(4)：进行2.5MPa压力点的气驱水实验；
步骤(5)：采用低场核磁共振方式对气驱水时的岩心进行检测，其中低场核磁共振检测参数的设置与步骤(2)中所设置的参数相同。采集驱替压力2.5MPa时岩心的T2谱，通过公式(2)计算2.5MPa时岩心中水核磁信号峰面积b(8770)；
步骤(6)：通过公式(3)、(4)计算气驱水时岩心含水饱和度、含气饱和度(Φw：54.3％；Φg：45.7％)
公式说明：
符号说明：
a：饱和水状态的待测岩心中水核磁信号峰面积，无量纲；
b：气驱水时待测岩心中水核磁信号峰面积，无量纲；
t：饱和水状态的横向弛豫时间最大值，ms；
t'：气驱水时的横向弛豫时间最大值，ms；
T2：代表饱和水状态的横向弛豫时间，ms；
T2'：代表气驱水时的横向弛豫时间，ms；
A(T2)：横向弛豫时间为T2时刻的饱和水岩心中水信号强度，无量纲；
B(T2')：横向弛豫时间为T2'时刻的气驱水时岩心中水信号强度，无量纲；
Φw：岩心含水饱和度，％；
Φg：岩心含气饱和度，％。
按照上述实施例的方法，分别对三个设计点(2.5MPa、5MPa、12MPa)进行基于核磁共振的气驱水时岩心含气饱和度在线测定，并对每个设计点进行常规称重法含气饱和度测定，与采用本发明的方法进行含气饱和度测定的结果进行分析比较，测定数据请参见表1，并将结果绘制成图3，可见两种方法所测结果的差别，本发明提供的技术方案提高了含气饱和度测定的精度。
表1
| 含气饱和度值 | 设计点1(2.5MPa) | 设计点2(5MPa) | 设计点3(12MPa) |
| 称重法 | 40.1％ | 46.8％ | 54.3％ |
| 本发明的方法 | 45.7％ | 53.0％ | 61.6％ |
基于同一发明构思，本发明实施例中还提供了一种基于核磁共振的气驱水时岩心含气饱和度测定装置，如下面的实施例。由于基于核磁共振的气驱水时岩心含气饱和度测定装置解决问题的原理与上述基于核磁共振的气驱水时岩心含气饱和度测定方法相似，因此基于核磁共振的气驱水时岩心含气饱和度测定装置的实施可以参考上述基于核磁共振的气驱水时岩心含气饱和度测定方法的实施，重复之处不再赘述。以下所使用的，术语“模块”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
图4是本发明实施例中基于核磁共振的气驱水时岩心含气饱和度测定装置的结构示意图，如图4所示，该装置包括：
第一采集单元01，用于采集低场核磁共振方式下饱和水状态的待测岩心T2谱；
饱和水状态峰面积确定单元02，用于根据所述饱和水状态的待测岩心T2谱，确定饱和水状态的待测岩心中水核磁信号峰面积；
第二采集单元03，用于采集低场核磁共振方式下气驱水时的待测岩心T2谱；
气驱水时峰面积确定单元04，用于根据所述气驱水时的待测岩心T2谱，确定气驱水时的待测岩心中水核磁信号峰面积；
含气饱和度确定单元05，用于根据饱和水状态的待测岩心中水核磁信号峰面积和气驱水时的待测岩心中水核磁信号峰面积，确定气驱水时岩心含气饱和度；
其中，所述含气饱和度确定单元具体用于：
根据饱和水状态的待测岩心中水核磁信号峰面积和气驱水时的待测岩心中水核磁信号峰面积，确定气驱水时待测岩心含水饱和度；
根据所述气驱水时待测岩心含水饱和度，确定气驱水时待测岩心含气饱和度。
在一个实例中，所述饱和水状态的待测岩心被至于无磁夹持器内部；所述无磁夹持器内部的胶套为聚四氟乙烯材料制胶套；
所述气驱水时的待测岩心被至于无磁夹持器内部；所述无磁夹持器内部的胶套为聚四氟乙烯材料制胶套。
在一个实例中，所述第一采集单元具体用于：
配置低场核磁共振方式下的低场核磁共振检测参数与待测岩心物性特征及流体性质相匹配；
在低场核磁共振检测参数与待测岩心物性特征及流体性质相匹配后，采集低场核磁共振方式下饱和水状态的待测岩心T2谱；
所述第二采集单元具体用于：
配置低场核磁共振方式下的低场核磁共振检测参数与待测岩心物性特征及流体性质相匹配；
在低场核磁共振检测参数与待测岩心物性特征及流体性质相匹配后，采集低场核磁共振方式下气驱水时的待测岩心T2谱。
本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述基于核磁共振的气驱水时岩心含气饱和度测定方法。
本发明实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有执行基于核磁共振的气驱水时岩心含气饱和度测定方法的计算机程序。
本发明实施提供的技术方案的有益技术效果为：
采用低场核磁共振的方式检测气驱水时岩心中水的核磁信号强度，能够实现气驱水时岩心含气饱和度的在线测定，可避免称重法在测定含水量时因样品暴露在环境中的损失量，使含气饱和度的测定更加准确，同时还可以实现气驱水过程中岩心含气饱和度的连续检测。
本发明的方法适用于石油、地质、矿业领域在实验室进行气驱水实验时岩心含气饱和度的测定，有助于增强成藏物理模拟实验的效果，可以服务于储层评价和油气成藏机理研究。
显然，本领域的技术人员应该明白，上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明实施例不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于核磁共振的气驱水时岩心含气饱和度测定方法，其特征在于，包括：
采集低场核磁共振方式下饱和水状态的待测岩心T2谱；
根据所述饱和水状态的待测岩心T2谱，确定饱和水状态的待测岩心中水核磁信号峰面积；
采集低场核磁共振方式下气驱水时的待测岩心T2谱；
根据所述气驱水时的待测岩心T2谱，确定气驱水时的待测岩心中水核磁信号峰面积；
根据饱和水状态的待测岩心中水核磁信号峰面积和气驱水时的待测岩心中水核磁信号峰面积，确定气驱水时岩心含气饱和度；
其中，根据饱和水状态的待测岩心中水核磁信号峰面积和气驱水时的待测岩心中水核磁信号峰面积，确定气驱水时岩心含气饱和度，包括：
将气驱水时的待测岩心中水核磁信号峰面积占饱和水状态的待测岩心中水核磁信号峰面积的百分数作为气驱水时待测岩心含水饱和度；
将百分之一百与所述气驱水时待测岩心含水饱和度的差值作为气驱水时待测岩心含气饱和度。
2.如权利要求1所述的基于核磁共振的气驱水时岩心含气饱和度测定方法，其特征在于，根据所述饱和水状态的待测岩心T2谱，确定饱和水状态的待测岩心中水核磁信号峰面积，包括按照如下公式确定饱和水状态的待测岩心中水核磁信号峰面积：
其中，a代表饱和水状态的待测岩心中水核磁信号峰面积，T2代表饱和水状态的横向弛豫时间，t为饱和水状态的横向弛豫时间最大值，A(T2)代表饱和水状态的横向弛豫时间为T2时刻的水信号强度。
3.如权利要求1所述的基于核磁共振的气驱水时岩心含气饱和度测定方法，其特征在于，根据所述气驱水时的待测岩心T2谱，确定气驱水时的待测岩心中水核磁信号峰面积，包括按照如下公式确定气驱水时的待测岩心中水核磁信号峰面积：
其中，b代表气驱水时的待测岩心中水核磁信号峰面积，T2'代表气驱水时的横向弛豫时间，t'为气驱水时的横向弛豫时间最大值，B(T2')代表气驱水时的横向弛豫时间为T2'时刻的水信号强度。
4.如权利要求1所述的基于核磁共振的气驱水时岩心含气饱和度测定方法，所述饱和水状态的待测岩心被至于无磁夹持器内部；所述气驱水时的待测岩心被至于无磁夹持器内部；所述无磁夹持器内部的胶套为聚四氟乙烯材料制胶套。
5.如权利要求1所述的基于核磁共振的气驱水时岩心含气饱和度测定方法，采集低场核磁共振方式下饱和水状态的待测岩心T2谱，包括：
配置低场核磁共振方式下的低场核磁共振检测参数与待测岩心物性特征及流体性质相匹配；
在低场核磁共振检测参数与待测岩心物性特征及流体性质相匹配后，采集低场核磁共振方式下饱和水状态的待测岩心T2谱；
采集低场核磁共振方式下气驱水时的待测岩心T2谱，包括：
配置低场核磁共振方式下的低场核磁共振检测参数与待测岩心物性特征及流体性质相匹配；
在低场核磁共振检测参数与待测岩心物性特征及流体性质相匹配后，采集低场核磁共振方式下气驱水时的待测岩心T2谱。
6.一种基于核磁共振的气驱水时岩心含气饱和度测定装置，其特征在于，包括：
第一采集单元，用于采集低场核磁共振方式下饱和水状态的待测岩心T2谱；
饱和水状态峰面积确定单元，用于根据所述饱和水状态的待测岩心T2谱，确定饱和水状态的待测岩心中水核磁信号峰面积；
第二采集单元，用于采集低场核磁共振方式下气驱水时的待测岩心T2谱；
气驱水时峰面积确定单元，用于根据所述气驱水时的待测岩心T2谱，确定气驱水时的待测岩心中水核磁信号峰面积；
含气饱和度确定单元，用于根据饱和水状态的待测岩心中水核磁信号峰面积和气驱水时的待测岩心中水核磁信号峰面积，确定气驱水时岩心含气饱和度；
其中，所述含气饱和度确定单元具体用于：
将气驱水时的待测岩心中水核磁信号峰面积占饱和水状态的待测岩心中水核磁信号峰面积的百分数作为气驱水时待测岩心含水饱和度；
将百分之一百与所述气驱水时待测岩心含水饱和度的差值作为气驱水时待测岩心含气饱和度。
7.如权利要求6所述的基于核磁共振的气驱水时岩心含气饱和度测定装置，所述饱和水状态的待测岩心被至于无磁夹持器内部；所述气驱水时的待测岩心被至于无磁夹持器内部；所述无磁夹持器内部的胶套为聚四氟乙烯材料制胶套。
8.如权利要求6所述的基于核磁共振的气驱水时岩心含气饱和度测定装置，所述第一采集单元具体用于：
配置低场核磁共振方式下的低场核磁共振检测参数与待测岩心物性特征及流体性质相匹配；
在低场核磁共振检测参数与待测岩心物性特征及流体性质相匹配后，采集低场核磁共振方式下饱和水状态的待测岩心T2谱；
所述第二采集单元具体用于：
配置低场核磁共振方式下的低场核磁共振检测参数与待测岩心物性特征及流体性质相匹配；
在低场核磁共振检测参数与待测岩心物性特征及流体性质相匹配后，采集低场核磁共振方式下气驱水时的待测岩心T2谱。
9.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5任一所述方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至5任一所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811416441.4A CN111220639B (zh) | 2018-11-26 | 2018-11-26 | 基于核磁共振的气驱水时岩心含气饱和度测定方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811416441.4A CN111220639B (zh) | 2018-11-26 | 2018-11-26 | 基于核磁共振的气驱水时岩心含气饱和度测定方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111220639A CN111220639A (zh) | 2020-06-02 |
| CN111220639B true CN111220639B (zh) | 2024-03-26 |
Family
ID=70825460
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811416441.4A Active CN111220639B (zh) | 2018-11-26 | 2018-11-26 | 基于核磁共振的气驱水时岩心含气饱和度测定方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111220639B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112485281B (zh) * | 2020-10-19 | 2022-04-01 | 中国地质大学（北京） | 动态测定多孔介质中气体水合物饱和度和渗透率的方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1043564A (zh) * | 1988-12-19 | 1990-07-04 | 大庆石油管理局勘探开发研究院 | 用于测定油藏岩心油、气、水三相相对渗透率实验方法 |
| CN104316554A (zh) * | 2014-10-27 | 2015-01-28 | 中国石油天然气股份有限公司 | 储层孔隙水可动性测试方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| BR112012017778A2 (pt) * | 2010-01-22 | 2018-08-14 | Prad Res & Development Ltd | método para determinar capacidade de umectação de formações de rocha utilizando medições de ressonância magnética nuclear, método para determinar uma relaxatividade de superfície de uma formação de rocha de subsuperfície usando medições de ressonância magnética nuclear feitas de dentro de um furo de poço penetrando a formação de rocha, método para determinar uma relaxatividade de superfície de uma formação de rocha de subsuperfície, e método para determinar saturação de água e de hidrocarbonetos em uma formação de rocha de subsuperfície usando medições de tempo de relaxação de ressonância magnética nuclear (nmr) e medições de constante de difusão |
-
2018
- 2018-11-26 CN CN201811416441.4A patent/CN111220639B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1043564A (zh) * | 1988-12-19 | 1990-07-04 | 大庆石油管理局勘探开发研究院 | 用于测定油藏岩心油、气、水三相相对渗透率实验方法 |
| CN104316554A (zh) * | 2014-10-27 | 2015-01-28 | 中国石油天然气股份有限公司 | 储层孔隙水可动性测试方法 |
Non-Patent Citations (2)
| Title |
|---|
| 基于核磁共振技术的致密砂岩气运移特征研究;房涛等;《甘肃科学学报》;第29卷(第4期);第22-27页 * |
| 用核磁共振技术测量低渗含水气藏中的束缚水饱和度;杨正明等;《石油钻采工艺》;第30卷(第3期);第56-59页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111220639A (zh) | 2020-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10161891B1 (en) | Method for characterizing rock physical characteristics of deeply buried carbonate rocks | |
| WO2018164953A1 (en) | Absolute porosity and pore size determination of pore types in media with varying pore sizes using nmr | |
| CN106153662A (zh) | 岩心应力敏感性的测量方法 | |
| RU2432587C2 (ru) | Способ анализа данных с общими и различающимися свойствами | |
| CN104297280B (zh) | 利用核磁共振技术定量评价岩心洗油效果的方法 | |
| NO20110236A1 (no) | Monte Carlo-metoden for Laplace-inversjonen av NMR-data | |
| CN106442599B (zh) | 岩石含油量测定方法及装置 | |
| US10534055B2 (en) | NMR method for determining non-oil volume of a rock sample | |
| CN111337408A (zh) | 一种利用低场核磁共振设备测试岩石裂缝孔隙度的方法 | |
| CN106897531B (zh) | 一种低渗透石灰岩储层渗透率的定量评价方法 | |
| CN107782705A (zh) | 岩石含油量测定方法及装置 | |
| US20240027379A1 (en) | Method for quantitative evaluation on sensitivity of shale oil and gas reservoir to injected fluids | |
| CN113189129A (zh) | 一种岩石裂缝孔隙度检测工艺 | |
| CN111220639B (zh) | 基于核磁共振的气驱水时岩心含气饱和度测定方法及装置 | |
| CN106644879A (zh) | 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置 | |
| CN106290103A (zh) | 一种页岩气储层中粘土微孔孔隙度的测定方法 | |
| CN105004747A (zh) | 一种核磁共振测量煤芯平均孔隙压缩系数的方法 | |
| CN105241912A (zh) | 低场核磁共振测量页岩有机质含量的方法及装置 | |
| CN105758883B (zh) | 含油污泥水洗过程中含油率检测方法及装置 | |
| CN109901238A (zh) | 一种基于应力差电阻率实验的高应力地层电阻率校正方法 | |
| CN106066494A (zh) | 一种火成岩核磁共振孔隙度校正方法和t2分布校正方法 | |
| CN108918574A (zh) | 核磁共振测量原油含氢指数的方法 | |
| CN110410058B (zh) | 一种校正岩心实验结果刻度二维核磁测井的方法 | |
| EP2437078A1 (en) | NMR method for porous media | |
| CN110609332B (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 |