CN111855523A - 渗流实验的分析方法及其应用 - Google Patents
渗流实验的分析方法及其应用 Download PDFInfo
- Publication number
- CN111855523A CN111855523A CN201910347235.0A CN201910347235A CN111855523A CN 111855523 A CN111855523 A CN 111855523A CN 201910347235 A CN201910347235 A CN 201910347235A CN 111855523 A CN111855523 A CN 111855523A
- Authority
- CN
- China
- Prior art keywords
- fluid
- magnetic resonance
- nuclear magnetic
- sample
- saturation
- 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.)
- Pending
Links
- 238000002474 experimental method Methods 0.000 title claims abstract description 28
- 238000004458 analytical method Methods 0.000 title claims abstract description 19
- 239000012530 fluid Substances 0.000 claims abstract description 156
- 238000005481 NMR spectroscopy Methods 0.000 claims abstract description 89
- 238000001514 detection method Methods 0.000 claims abstract description 52
- 238000000034 method Methods 0.000 claims abstract description 48
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 44
- 230000008569 process Effects 0.000 claims abstract description 24
- 238000002347 injection Methods 0.000 claims abstract description 6
- 239000007924 injection Substances 0.000 claims abstract description 6
- 238000012545 processing Methods 0.000 claims abstract description 6
- 239000000523 sample Substances 0.000 claims description 100
- 238000009826 distribution Methods 0.000 claims description 44
- 238000001228 spectrum Methods 0.000 claims description 43
- 239000011435 rock Substances 0.000 claims description 15
- 238000000685 Carr-Purcell-Meiboom-Gill pulse sequence Methods 0.000 claims description 14
- 238000005070 sampling Methods 0.000 claims description 13
- 239000002689 soil Substances 0.000 claims description 11
- 238000005325 percolation Methods 0.000 claims description 9
- 238000009825 accumulation Methods 0.000 claims description 6
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 125000004435 hydrogen atom Chemical group [H]* 0.000 claims description 4
- 238000005520 cutting process Methods 0.000 claims description 3
- 239000003673 groundwater Substances 0.000 claims description 3
- 238000003895 groundwater pollution Methods 0.000 claims description 3
- 238000003672 processing method Methods 0.000 claims description 3
- 238000001764 infiltration Methods 0.000 claims 3
- 230000008595 infiltration Effects 0.000 claims 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 28
- 239000003921 oil Substances 0.000 description 16
- 239000000243 solution Substances 0.000 description 9
- 239000003350 kerosene Substances 0.000 description 8
- 239000007789 gas Substances 0.000 description 7
- 238000003911 water pollution Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 4
- XLYOFNOQVPJJNP-ZSJDYOACSA-N Heavy water Chemical compound [2H]O[2H] XLYOFNOQVPJJNP-ZSJDYOACSA-N 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000002592 echocardiography Methods 0.000 description 3
- 239000011148 porous material Substances 0.000 description 3
- 230000002902 bimodal effect Effects 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000003889 chemical engineering Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 239000013256 coordination polymer Substances 0.000 description 1
- 239000010779 crude oil Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000008398 formation water Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000003014 reinforcing effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000004753 textile Substances 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
-
- 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
- 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
Landscapes
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- High Energy & Nuclear Physics (AREA)
- Dispersion Chemistry (AREA)
- Environmental & Geological Engineering (AREA)
- Engineering & Computer Science (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种渗流实验的分析方法及其应用。该渗流实验的分析方法包括以下步骤：(a)使用孔隙度定标样品定标，确定核磁共振信号量与孔隙度的函数关系式；(b)向饱和第一流体的待测样品中注入第二流体，注入过程中，在所需时间节点上采集核磁共振信号量；其中，第一流体或第二流体中含有能够被核磁共振仪探测到的探测对象；(c)根据步骤(a)中的函数关系式以及步骤(b)中采集到的核磁共振信号量，计算出待测样品的孔隙度；(d)对步骤(b)中采集到的核磁共振信号量进行数据处理，计算得到所需时间节点的含探测对象流体的流体饱和度。该方法所得待测样品的孔隙度和流体饱和度精确度高，且能对渗流实验过程中流体的动态过程进行监测。
Description
技术领域
本发明涉及渗流特性研究领域，具体而言，涉及一种渗流实验的分析方法及其应用。
背景技术
渗流是指流体在孔隙介质中的流动，渗流在水利、地质、采矿、石油、环境保护、化工、生物、医疗等领域都有广泛的应用，如开发利用地下水资源、防止建筑物地基发生渗透变形、基坑排水等均需应用到渗流理论。流体饱和度是用来描述储层岩石孔隙中流体充满的程度，该参数影响油气藏储量的大小。当储层岩石孔隙中同时存在多种流体(原油、地层水或天然气)时，某种流体所占的体积百分数称为该种流体的饱和度。
传统渗流实验过程中，对样品内部流体饱和度的监测手段为，通过计量流体泵的注入量或者渗流尾端流出的流体量，间接计算样品内部的流体饱和度，但是该方法精确度较低。
核磁共振技术作为一种先进的无损检测手段，已经应用在医药、化工、生物、食品、纺织、能源地矿等多领域。然而现有的核磁共振技术仅能够检测到岩心等材料静态时的孔隙度和流体饱和度，很难对渗流实验过程中流体的动态过程进行监测。
有鉴于此，特提出本发明。
发明内容
本发明的第一目的在于提供一种渗流实验的分析方法，该方法采用核磁共振技术来测定待测样品的孔隙度和流体的饱和度，所得待测样品的孔隙度和流体饱和度精确度高，且能对渗流实验过程中流体的动态过程进行监测，得到渗流过程中流体饱和度的变化规律。
本发明的第二目的在于提供一种上述渗流实验的分析方法的应用。
为了实现本发明的上述目的，特采用以下技术方案：
第一方面，本发明提供了一种渗流实验的分析方法，包括以下步骤：
(a)使用孔隙度定标样品定标，确定核磁共振信号量与孔隙度的函数关系式；
(b)向饱和第一流体的待测样品中注入第二流体，注入过程中，在所需时间节点上采集核磁共振信号量；
其中，第一流体或第二流体中含有能够被核磁共振仪探测到的探测对象；
(c)根据步骤(a)中的函数关系式以及步骤(b)中采集到的核磁共振信号量，计算出待测样品的孔隙度；
(d)对步骤(b)中采集到的核磁共振信号量进行数据处理，计算得到所需时间节点的含探测对象流体的流体饱和度。
作为进一步优选的技术方案，在步骤(a)之前还依次包括校准核磁共振仪的基本参数、选择脉冲序列以及设置采样参数的步骤。
作为进一步优选的技术方案，所述基本参数包括中心频率和/或脉宽；
优选地，所述脉冲序列包括CPMG序列；
优选地，所述采样参数包括序列的回波间隔、回波个数、等待时间、增益或累加次数中的至少一种。
作为进一步优选的技术方案，步骤(b)和/或步骤(d)中，所述探测对象包括1H原子核。
作为进一步优选的技术方案，步骤(a)中，定标方法包括：采用核磁共振仪测量定标样品的核磁共振信号量，其中定标样品的孔隙度和体积已知，得到单位体积核磁共振信号量与孔隙度之间的关系式y＝ax+b，式中y为单位体积核磁共振信号量，x为定标样品孔隙度，a为斜率，b为截距。
作为进一步优选的技术方案，步骤(c)中的核磁共振信号量为待测样品饱和含探测对象的流体时所采集到的核磁共振信号量。
作为进一步优选的技术方案，步骤(c)中的计算方法为：x样＝(y饱和-b)/a，式中x样为待测样品孔隙度，y饱和为待测样品饱和含探测对象的流体时采集到的单位体积核磁共振信号量，a为步骤(a)函数关系式中的斜率，b为步骤(a)函数关系式中的截距。
作为进一步优选的技术方案，步骤(d)中，所述含探测对象流体的流体饱和度包括总含探测对象流体的流体饱和度、束缚含探测对象流体的流体饱和度和可动含探测对象流体的流体饱和度。
作为进一步优选的技术方案，步骤(d)中，所述数据处理的方法包括：对采集到的核磁共振信号量进行拉普拉斯反演，得到待测样品的T2分布谱，根据谱图形态确定T2截止值，得到含探测对象流体的流体饱和度，其中：
总含探测对象流体的流体饱和度＝渗流某一时间节点的T2分布谱总面积/饱和含探测对象流体的T2分布谱总面积，单位为百分比；
束缚含探测对象流体的流体饱和度＝渗流某一时间节点的T2分布谱T2小于等于T2截止值的面积/饱和含探测对象流体的T2分布谱总面积，单位为百分比；
可动含探测对象流体的流体饱和度＝渗流某一时间节点的T2分布谱T2大于T2截止值的面积/饱和含探测对象流体的T2分布谱总面积，单位为百分比。
第二方面，本发明提供了一种上述渗流实验的分析方法在油气开采、地下水开采、岩土体稳定性控制或地下水污染控制中的应用。
与现有技术相比，本发明的有益效果为：
本发明提供的渗流实验的分析方法利用定标样品确定出核磁共振信号量与孔隙度的函数关系式，然后测得待测样品饱和含探测对象的流体的核磁共振信号量，将该核磁共振信号量代入到上述函数关系式中得到待测样品的孔隙度，该方法仅需采集不同待测样品饱和含探测对象的流体的核磁共振信号量即可以方便地得到不同待测样品的孔隙度。同时，根据不同时间节点上的核磁共振信号量，经过数据处理后，得到所需时间节点的含探测对象流体的流体饱和度。
该方法采用核磁共振技术来测定待测样品的孔隙度和流体的饱和度，无损、方便、快捷，所得待测样品的孔隙度和流体饱和度精确度高，且能对渗流实验过程中流体的动态过程进行监测，得到渗流过程中流体饱和度的变化规律。
将上述渗流实验的分析方法应用于油气开采、地下水开采、岩土体稳定性控制或地下水污染控制中，能够提高油气开采和地下水开采的效率和精度，能够对岩土体稳定性控制和地下水污染控制提供有效参考，并提高岩土体稳定性控制和地下水污染控制的有效性和精确度。
附图说明
图1是本发明实施例1渗流过程中的T2分布谱；
图2是本发明实施例2渗流过程中的T2分布谱。
具体实施方式
下面将结合实施例对本发明的实施方案进行详细描述，但是本领域技术人员将会理解，下列实施例仅用于说明本发明，而不应视为限制本发明的范围。实施例中未注明具体条件者，按照常规条件或制造商建议的条件进行。
需要说明的是：
本发明中，如果没有特别的说明，本文所提到的所有实施方式以及优选实施方法可以相互组合形成新的技术方案。
本发明中，如果没有特别的说明，本文所提到的所有技术特征以及优选特征可以相互组合形成新的技术方案。
本发明中，除非另有说明，各个操作步骤可以顺序进行，也可以不按照顺序进行。优选地，本文中的操作步骤是顺序进行的。
除非另有说明，本文中所用的专业与科学术语与本领域熟练人员所熟悉的意义相同。此外，任何与所记载内容相似或均等的方法或材料也可应用于本发明中。
第一方面，在至少一个实施例中提供了一种渗流实验的分析方法，包括以下步骤：
(a)使用孔隙度定标样品定标，确定核磁共振信号量与孔隙度的函数关系式；
(b)向饱和第一流体的待测样品中注入第二流体，注入过程中，在所需时间节点上采集核磁共振信号量；
其中，第一流体或第二流体中含有能够被核磁共振仪探测到的探测对象；
(c)根据步骤(a)中的函数关系式以及步骤(b)中采集到的核磁共振信号量，计算出待测样品的孔隙度；
(d)对步骤(b)中采集到的核磁共振信号量进行数据处理，计算得到所需时间节点的含探测对象流体的流体饱和度。
需要说明的是：
上述“孔隙度定标样品”是指孔隙度数值及体积是已知并且准确的样品。
上述“定标”是指确定孔隙度定标样品的孔隙度所对应的核磁共振信号量参考值。
上述“饱和第一流体的待测样品”是指待测样品中的第一流体达到饱和状态，该饱和状态是指待测样品中能够容纳的第一流体的最大体积量。
上述“第一流体”和“第二流体”仅为表述方便的目的而设，并不意味着其具有重要性的差别，即不意味着第一流体更重要，也不意味着第二流体更重要。
本发明中所涉及的“核磁共振信号量”也可称为回波串数据、核磁共振数据或核磁共振信号。
“第一流体或第二流体中含有能够被核磁共振仪探测到的探测对象”是指第一流体或第二流体能够被核磁共振仪采集到核磁共振信号量，上述探测对象包括但不限于1H原子核、19F、23Na或31P等质子数与中子数不同时为偶数的原子核。
“含探测对象流体的流体饱和度”是指第一流体饱和度或第二流体饱和度。
待测样品的孔隙度是在待测样品饱和含探测对象的流体后得到的。
当第一流体中含有能够被核磁共振仪探测到的探测对象时，则饱和第一流体的待测样品在注入第二流体的时间节点为0时，即当还未注入第二流体时或渗流过程还未开始时，将此时采集到的核磁共振信号量代入步骤(a)中的函数关系式中，得到待测样品的孔隙度。相应的，得到的所需时间节点的饱和度为第一流体饱和度。
当第二流体中含有能够被核磁共振仪探测到的探测对象时，则饱和第一流体的待测样品在注入第二流体的时间节点为终止时间时，即当渗流过程结束时，将此时采集到的核磁共振信号量代入步骤(a)中的函数关系式中，得到待测样品的孔隙度。相应的，得到的所需时间节点的饱和度为第二流体饱和度。
上述渗流实验的分析方法利用定标样品确定出核磁共振信号量与孔隙度的函数关系式，然后测得待测样品饱和含探测对象的流体的核磁共振信号量，将该核磁共振信号量代入到上述函数关系式中得到待测样品的孔隙度，该方法仅需采集不同待测样品饱和含探测对象的流体的核磁共振信号量即可以方便地得到不同待测样品的孔隙度。同时，根据不同时间节点上的核磁共振信号量，经过数据处理后，得到所需时间节点的含探测对象流体的流体饱和度。
该方法采用核磁共振技术来测定待测样品的孔隙度和流体的饱和度，无损、方便、快捷，所得待测样品的孔隙度和流体饱和度精确度高，且由于能采集所需时间点的核磁共振信号量，因而能对渗流实验过程中流体的动态过程进行监测，得到渗流过程中流体饱和度的变化规律。
在一种优选的实施方式中，在步骤(a)之前还依次包括校准核磁共振仪的基本参数、选择脉冲序列以及设置采样参数的步骤。在依次经过上述步骤之后，即可使用孔隙度定标样品定标。
优选地，所述基本参数包括中心频率和/或脉宽；
优选地，所述脉冲序列包括CPMG序列。CPMG序列能够用于丈量T2值，是CP序列的改良，用Y轴位置的180°脉冲代替相位交替的180°±X脉冲，发生一系列的自旋回波，用CPMG序列发生的回波误差小且无积累。
优选地，所述采样参数包括序列的回波间隔、回波个数、等待时间、增益或累加次数中的至少一种。
在一种优选的实施方式中，步骤(b)和/或步骤(d)中，所述探测对象包括1H原子核。核磁共振信号量与含1H原子核的流体量成正比，监测和换算更加简单。
应当理解的是，定标样品定标时所采用的探测对象与步骤(b)和/或步骤(d)中的探测对象一致。因此，相应的，定标样品定标时，所采用的探测对象为1H原子核。
在一种优选的实施方式中，步骤(a)中，定标方法包括：采用核磁共振仪测量定标样品的核磁共振信号量，其中定标样品的孔隙度和体积已知，得到单位体积核磁共振信号量与孔隙度之间的关系式y＝ax+b，式中y为单位体积核磁共振信号量，x为定标样品孔隙度，a为斜率，b为截距。采用定标样品定标后，即可得到单位体积核磁共振信号量与孔隙度之间的函数关系式，当对不同待测样品进行渗流实验时，均采用上述关系式即可得到不同待测样品的孔隙度，因此，测定不同待测样品孔隙度的方法简便、效率高、误差小。
优选地，步骤(c)中的核磁共振信号量为待测样品饱和含探测对象的流体时所采集到的核磁共振信号量。上述“待测样品饱和含探测对象的流体”是指待测样品中的含探测对象的流体达到饱和状态，该饱和状态是指待测样品中能够容纳的含探测对象的流体的最大体积量。
优选地，步骤(c)中的计算方法为：x样＝(y饱和-b)/a，式中x样为待测样品孔隙度，y饱和为待测样品饱和含探测对象的流体时采集到的单位体积核磁共振信号量，a为步骤(a)函数关系式中的斜率，b为步骤(a)函数关系式中的截距。采用本优选实施方式中的方法即可反推出待测样品的孔隙度。
在一种优选的实施方式中，步骤(d)中，所述含探测对象流体的流体饱和度包括总含探测对象流体的流体饱和度、束缚含探测对象流体的流体饱和度和可动含探测对象流体的流体饱和度。
优选地，步骤(d)中，所述数据处理的方法包括：对采集到的核磁共振信号量进行拉普拉斯反演，得到待测样品的T2分布谱，根据谱图形态确定T2截止值，得到含探测对象流体的流体饱和度，其中：
总含探测对象流体的流体饱和度＝渗流某一时间节点的T2分布谱总面积/饱和含探测对象流体的T2分布谱总面积，单位为百分比；
束缚含探测对象流体的流体饱和度＝渗流某一时间节点的T2分布谱T2小于等于T2截止值的面积/饱和含探测对象流体的T2分布谱总面积，单位为百分比；
可动含探测对象流体的流体饱和度＝渗流某一时间节点的T2分布谱T2大于T2截止值的面积/饱和含探测对象流体的T2分布谱总面积，单位为百分比。
需要说明的是：
“拉普拉斯反演”是指拉普拉斯变换的反演。本发明不对上述“拉普拉斯反演”做特别限制，采用现有的任意一种可实现的反演方法即可。
根据谱图形态确定T2截止值的优点在于：能够快速低成本的得到T2截止值，并且相对于经验值，其精确度更高。
现有的通过计量流体泵的注入量或者渗流尾端流出流体量的方式，仅能间接计算待测样品内部的总体流体饱和度，无法确定待测样品内部束缚流体和可动流体饱和度。通过本发明上述优选方法能够准确得到待测样品的总含流体饱和度、束缚流体饱和度和可动流体饱和度。
优选地，当T2分布谱为单峰时，当分布谱形态为对称分布时，T2截止值与峰值位置一致；
当分布谱形态为非对称分布，且峰值形态左偏时，T2截止值位于峰值右侧斜坡段的半幅点处；
当分布谱形态为非对称分布，且峰值形态右偏时，T2截止值位于峰值左侧斜坡段的半幅点处；
优选地，当T2分布谱为双峰时，当双峰值明显，且存在明显波谷时，T2截止值位于波谷处；
当一个波峰明显，一个波峰不明显，有较为明显的平坦段，低峰/高峰的幅度比值大于0.5时，T2截止值位于斜坡段与平坦段的拐点处；
当可动峰明显，束缚峰不明显，有较为明显的斜坡段，低峰/高峰的幅度比值小于0.5时，T2截止值位于可动峰左侧斜坡段的半幅点处；
当束缚峰明显，可动峰不明显，有较为明显的斜坡段，低峰/高峰的幅度比值小于0.5时，T2截止值位于束缚峰右侧斜坡段顶部；
优选地，当T2分布谱为三峰时，当三峰明显，第1峰和第2峰在20ms左侧，第3峰在20ms右侧时，T2截止值位于第2峰和第3峰的波谷位置；
当三峰明显，第1峰在10ms左侧，第2峰在20-40ms之间，第3峰大于40ms时，T2截止值与第2峰峰值一致；
当第1峰在10ms左侧，第3峰在50ms右侧，第2峰峰值不明显，但与第3峰之间存在明显的拐点时，T2截止值位于第2峰和第3峰的拐点处；
当第1峰在20ms左侧，与第2峰存在明显的波谷，第2峰和第3峰在40ms右侧时，T2截止值位于第1峰和第3峰的波谷处。
第二方面，本发明提供了一种上述渗流实验的分析方法在油气开采、地下水开采、岩土体稳定性控制或地下水污染控制中的应用。
将上述渗流实验的分析方法应用于油气开采、地下水开采、岩土体稳定性控制或地下水污染控制中，能够提高油气开采和地下水开采的效率和精度，能够对岩土体稳定性控制和地下水污染控制提供有效参考，并提高岩土体稳定性控制和地下水污染控制的有效性和精确度。
需要说明的是：上述“油气开采”是指石油和/或天然气的开采；上述“岩土体稳定性控制”是指岩石和/或土体的稳定性控制，例如，加固建筑物地基或基坑排水等。
下面结合实施例对本发明做进一步的说明。
实施例1
一种氮气渗流饱水岩心实验的分析方法：
(1)启动核磁共振仪，放入校准油样，校准核磁共振仪的中心频率、脉宽等基本参数；
(2)选择CPMG脉冲序列，设置序列的回波间隔为0.2ms、回波个数为10000、等待时间为6000ms、增益15dB、累加次数为16；
(3)依次将孔隙度定标样品放入核磁共振仪使用CPMG序列进行采样，得到单位体积核磁共振信号量与孔隙度之间的函数关系式：y＝2.356x+0.618，式中x为孔隙度，y为单位体积核磁共振信号量；
(4)将饱水岩心放入核磁共振仪，使用CPMG序列进行采样，得到饱水岩心单位体积核磁共振信号量y饱和为59.6682，代入步骤(3)中的函数关系式计算得到待测样品孔隙度x样为25.06％；
(5)以0.5MPa的压力向饱水岩心中注入氮气，分别在1min、5min、10min、30min使用CPMG序列进行采样；
(6)对步骤(4)和(5)中得到的核磁共振数据进行一维拉普拉斯反演，得到不同时间节点的T2分布谱(如图1所示)，根据谱图形态确定T2截止值为47ms，根据以下公式计算出不同时间节点的各饱和度参数，结果列于表1：
含水饱和度＝渗流某一时间节点的T2分布谱总面积/饱水样的T2分布谱总面积；
束缚水饱和度＝渗流某一时间节点的T2分布谱T2小于等于T2截止值的面积/饱水样T2分布谱总面积；
可动水饱和度＝渗流某一时间节点的T2分布谱T2大于T2截止值的面积/饱水样T2分布谱总面积。
表1渗流过程中的饱和度计算结果
实施例2
一种重水渗流饱煤油岩心实验的分析方法：
(1)启动核磁共振仪，放入校准油样，校准核磁共振仪的中心频率、脉宽等基本参数；
(2)选择CPMG脉冲序列，设置序列的回波间隔为0.2ms、回波个数为10000、等待时间为6000ms、增益20dB、累加次数为32；
(3)依次将孔隙度定标样品放入核磁共振仪使用CPMG序列进行采样，得到单位体积核磁共振信号量与孔隙度之间的函数关系式：y＝431.93x+52.031，式中x为孔隙度，y为单位体积核磁共振信号量；
(4)将饱煤油岩心放入核磁共振仪，使用CPMG序列进行采样，得到饱煤油岩心单位体积核磁共振信号量y饱和为2450.129，代入步骤(3)中的函数关系式计算得到待测样品孔隙度x样为5.62％；
(5)以2mL/min的流速向饱煤油岩心中注入重水，分别在1min、5min、10min、30min使用CPMG序列进行采样；
(6)对步骤(4)和(5)中得到的核磁共振数据进行一维拉普拉斯反演，得到不同时间节点的T2分布谱(如图2所示)，根据谱图形态确定T2截止值为4.9ms，根据以下公式计算出不同时间节点的各饱和度参数，结果列于表2：
含油饱和度＝渗流某一时间节点的T2分布谱总面积/饱煤油样的T2分布谱总面积；
束缚油饱和度＝渗流某一时间节点的T2分布谱T2小于等于T2截止值的面积/饱煤油样T2分布谱总面积；
可动油饱和度＝渗流某一时间节点的T2分布谱T2大于T2截止值的面积/饱煤油样T2分布谱总面积。
表2渗流过程中的饱和度计算结果
| 编号 | 时间/min | 总含油饱和度 | 束缚油饱和度 | 可动油饱和度 |
| 1 | 1 | 93.96％ | 43.01％ | 50.95％ |
| 2 | 5 | 80.31％ | 42.10％ | 38.21％ |
| 3 | 10 | 68.92％ | 41.33％ | 27.59％ |
| 4 | 30 | 44.21％ | 35.70％ | 8.51％ |
对比例1
该对比例须与实施例1同时进行操作，增加步骤为：
1、分别测试岩心饱和水前后的质量m前、m饱和，饱和进入到岩心中的总水质量则为m水＝m饱和-m前；
2、在渗流尾端放置10ml量程量筒，用于接盛渗流流出的水；
3、在1min、5min、10min、30min使用CPMG序列采样的同时，读取量筒中水的量，记录于下表，并根据公式计算各时间节点的含水饱和度，注意该方法仅能得到总含水饱和度：
含水饱和度＝(饱和进入到岩心中的总水质量-渗流某一时间节点量筒中水质量)/饱和进入到岩心中的总水质量。
| 编号 | 时间/min | 水质量/g | 总含水饱和度 |
| 1 | 1 | 0 | 100％ |
| 2 | 5 | 1.1 | 82.03％ |
| 3 | 10 | 2.2 | 64.05％ |
| 4 | 30 | 4.1 | 33.01％ |
由于渗流出的水部分存在于管路中，且无法准确计量，所以导致该方法得到的含水饱和度存在偏差较大。
尽管已用具体实施例来说明和描述了本发明，然而应意识到，在不背离本发明的精神和范围的情况下可以作出许多其它的更改和修改。因此，这意味着在所附权利要求中包括属于本发明范围内的所有这些变化和修改。
Claims (10)
1.一种渗流实验的分析方法，其特征在于，包括以下步骤：
(a)使用孔隙度定标样品定标，确定核磁共振信号量与孔隙度的函数关系式；
(b)向饱和第一流体的待测样品中注入第二流体，注入过程中，在所需时间节点上采集核磁共振信号量；
其中，第一流体或第二流体中含有能够被核磁共振仪探测到的探测对象；
(c)根据步骤(a)中的函数关系式以及步骤(b)中采集到的核磁共振信号量，计算出待测样品的孔隙度；
(d)对步骤(b)中采集到的核磁共振信号量进行数据处理，计算得到所需时间节点的含探测对象流体的流体饱和度。
2.根据权利要求1所述的渗流实验的分析方法，其特征在于，在步骤(a)之前还依次包括校准核磁共振仪的基本参数、选择脉冲序列以及设置采样参数的步骤。
3.根据权利要求2所述的渗流实验的分析方法，其特征在于，所述基本参数包括中心频率和/或脉宽；
优选地，所述脉冲序列包括CPMG序列；
优选地，所述采样参数包括序列的回波间隔、回波个数、等待时间、增益或累加次数中的至少一种。
4.根据权利要求1所述的渗流实验的分析方法，其特征在于，步骤(b)和/或步骤(d)中，所述探测对象包括1H原子核。
5.根据权利要求1-4任一项所述的渗流实验的分析方法，其特征在于，步骤(a)中，定标方法包括：采用核磁共振仪测量定标样品的核磁共振信号量，其中定标样品的孔隙度和体积已知，得到单位体积核磁共振信号量与孔隙度之间的关系式y＝ax+b，式中y为单位体积核磁共振信号量，x为定标样品孔隙度，a为斜率，b为截距。
6.根据权利要求1-4任一项所述的渗流实验的分析方法，其特征在于，步骤(c)中的核磁共振信号量为待测样品饱和含探测对象的流体时所采集到的核磁共振信号量。
7.根据权利要求5所述的渗流实验的分析方法，其特征在于，步骤(c)中的计算方法为：x样＝(y饱和-b)/a，式中x样为待测样品孔隙度，y饱和为待测样品饱和含探测对象的流体时采集到的单位体积核磁共振信号量，a为步骤(a)函数关系式中的斜率，b为步骤(a)函数关系式中的截距。
8.根据权利要求1-4任一项所述的渗流实验的分析方法，其特征在于，步骤(d)中，所述含探测对象流体的流体饱和度包括总含探测对象流体的流体饱和度、束缚含探测对象流体的流体饱和度和可动含探测对象流体的流体饱和度。
9.根据权利要求8所述的渗流实验的分析方法，其特征在于，步骤(d)中，所述数据处理的方法包括：对采集到的核磁共振信号量进行拉普拉斯反演，得到待测样品的T2分布谱，根据谱图形态确定T2截止值，得到含探测对象流体的流体饱和度，其中：
总含探测对象流体的流体饱和度＝渗流某一时间节点的T2分布谱总面积/饱和含探测对象流体的T2分布谱总面积，单位为百分比；
束缚含探测对象流体的流体饱和度＝渗流某一时间节点的T2分布谱T2小于等于T2截止值的面积/饱和含探测对象流体的T2分布谱总面积，单位为百分比；
可动含探测对象流体的流体饱和度＝渗流某一时间节点的T2分布谱T2大于T2截止值的面积/饱和含探测对象流体的T2分布谱总面积，单位为百分比。
10.权利要求1-9任一项所述的渗流实验的分析方法在油气开采、地下水开采、岩土体稳定性控制或地下水污染控制中的应用。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910347235.0A CN111855523A (zh) | 2019-04-26 | 2019-04-26 | 渗流实验的分析方法及其应用 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910347235.0A CN111855523A (zh) | 2019-04-26 | 2019-04-26 | 渗流实验的分析方法及其应用 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111855523A true CN111855523A (zh) | 2020-10-30 |
Family
ID=72964879
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910347235.0A Pending CN111855523A (zh) | 2019-04-26 | 2019-04-26 | 渗流实验的分析方法及其应用 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111855523A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20200363352A1 (en) * | 2019-05-13 | 2020-11-19 | Exxonmobil Upstream Research Company | Advanced NMR Analysis of Porosity and Other Properties in Core Samples using Hydraulic Fluid Exchange |
| CN118243591A (zh) * | 2024-05-29 | 2024-06-25 | 苏州纽迈分析仪器股份有限公司 | 一种非规则形状岩心的孔隙度测试方法及装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6140817A (en) * | 1998-05-26 | 2000-10-31 | Schlumberger Technology Corporation | Magnetic resonance well logging method and apparatus |
| CN105181728A (zh) * | 2015-07-24 | 2015-12-23 | 中国石油天然气股份有限公司 | 核磁共振在线检测页岩气的方法 |
| CN105954170A (zh) * | 2016-05-04 | 2016-09-21 | 中国石油大学(华东) | 一种考虑背景噪声信号量的核磁孔隙度计算方法 |
| CN106290440A (zh) * | 2015-05-26 | 2017-01-04 | 中国石油天然气股份有限公司 | 驱油效果检测方法及装置 |
| CN108414560A (zh) * | 2018-03-06 | 2018-08-17 | 中国石油大学(华东) | 一种核磁-驱替联用装置评价致密油充注过程的方法 |
| CN109211748A (zh) * | 2017-07-07 | 2019-01-15 | 中国石油化工股份有限公司 | 致密砂岩储层孔隙油相动用率分析方法 |
-
2019
- 2019-04-26 CN CN201910347235.0A patent/CN111855523A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6140817A (en) * | 1998-05-26 | 2000-10-31 | Schlumberger Technology Corporation | Magnetic resonance well logging method and apparatus |
| CN106290440A (zh) * | 2015-05-26 | 2017-01-04 | 中国石油天然气股份有限公司 | 驱油效果检测方法及装置 |
| CN105181728A (zh) * | 2015-07-24 | 2015-12-23 | 中国石油天然气股份有限公司 | 核磁共振在线检测页岩气的方法 |
| CN105954170A (zh) * | 2016-05-04 | 2016-09-21 | 中国石油大学(华东) | 一种考虑背景噪声信号量的核磁孔隙度计算方法 |
| CN109211748A (zh) * | 2017-07-07 | 2019-01-15 | 中国石油化工股份有限公司 | 致密砂岩储层孔隙油相动用率分析方法 |
| CN108414560A (zh) * | 2018-03-06 | 2018-08-17 | 中国石油大学(华东) | 一种核磁-驱替联用装置评价致密油充注过程的方法 |
Non-Patent Citations (10)
| Title |
|---|
| 房涛等: "基于核磁共振技术的致密砂岩气运移特征研究", 《甘肃科学学报》 * |
| 房涛等: "基于核磁共振技术的致密砂岩气运移特征研究", 《甘肃科学学报》, no. 04, 25 August 2017 (2017-08-25), pages 24 - 29 * |
| 杨正明等: "用核磁共振技术测量低渗含水气藏中的束缚水饱和度", 《石油钻采工艺》 * |
| 杨正明等: "用核磁共振技术测量低渗含水气藏中的束缚水饱和度", 《石油钻采工艺》, no. 03, 15 June 2008 (2008-06-15), pages 62 - 65 * |
| 臧士宾等: "核磁共振技术在储层评价中的应用研究――以青海柴达木盆地为例", 《青海石油》 * |
| 臧士宾等: "核磁共振技术在储层评价中的应用研究――以青海柴达木盆地为例", 《青海石油》, no. 04, 15 December 2007 (2007-12-15), pages 10 - 14 * |
| 邵维志等: "利用T2谱形态确定T2截止值的方法探索", 《测井技术》, vol. 33, no. 5, pages 430 - 435 * |
| 霍迎冬: "核磁共振技术在致密油储层流体饱和度分析中的应用研究", 《西部探矿工程》 * |
| 霍迎冬: "核磁共振技术在致密油储层流体饱和度分析中的应用研究", 《西部探矿工程》, no. 06, 15 June 2018 (2018-06-15), pages 107 - 109 * |
| 霍迎冬;: "核磁共振技术在致密油储层流体饱和度分析中的应用研究", 西部探矿工程, no. 06, pages 107 - 109 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20200363352A1 (en) * | 2019-05-13 | 2020-11-19 | Exxonmobil Upstream Research Company | Advanced NMR Analysis of Porosity and Other Properties in Core Samples using Hydraulic Fluid Exchange |
| US11815481B2 (en) * | 2019-05-13 | 2023-11-14 | ExxonMobil Technology and Engineering Company | Advanced NMR analysis of porosity and other properties in core samples using hydraulic fluid exchange |
| CN118243591A (zh) * | 2024-05-29 | 2024-06-25 | 苏州纽迈分析仪器股份有限公司 | 一种非规则形状岩心的孔隙度测试方法及装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yan et al. | A robust NMR method to measure porosity of low porosity rocks | |
| CN109443867B (zh) | 一种对致密岩石的物性参数进行连续检测的方法 | |
| CN102998322B (zh) | 恒定梯度场核磁共振岩样分析方法及仪器 | |
| CN106050225B (zh) | 一种核磁共振测井t2谱100％纯水谱的确定方法 | |
| CN100349013C (zh) | 核磁共振测井t2谱t2截止值的确定方法 | |
| CN110296931B (zh) | 一种致密砂岩油水相对渗透率信息的表征方法及系统 | |
| US10247684B2 (en) | Nuclear magnetic resonance (NMR) distributions and pore information | |
| CN104777181B (zh) | 致密油核磁共振t2截止值及流体饱和度确定方法、装置 | |
| US20240027379A1 (en) | Method for quantitative evaluation on sensitivity of shale oil and gas reservoir to injected fluids | |
| SA519402506B1 (ar) | تحديد مسامية مطلقة وحجم مسامي لأنواع مسام ذات أحجام مسامية مختلفة | |
| CN106501144A (zh) | 一种基于核磁共振双截止值的致密砂岩渗透率计算方法 | |
| CN105804734A (zh) | 一种利用核磁共振测井识别稠油储层的方法 | |
| CN109682850B (zh) | 一种在线渗吸实验用核磁共振测试装置及实验方法 | |
| CN107655922B (zh) | 一种基于核磁共振的烃源岩含水量无损检测的方法 | |
| CN111855523A (zh) | 渗流实验的分析方法及其应用 | |
| CA3107041A1 (en) | Assessment of inaccessible pore volume for polymer flooding | |
| CN110687155A (zh) | 一种土壤薄膜水厚度的测试方法 | |
| CN104297278A (zh) | 改良式CPMG序列T2-G采集参数t0自适应方法 | |
| CN110161071A (zh) | 一种评价致密储层可动流体最小喉道半径的方法 | |
| CN103675722A (zh) | 岩石t2-g实验采集参数自动匹配方法 | |
| CN102590476A (zh) | 一种变性土膨胀系数测定装置及其应用 | |
| McPhee et al. | Nuclear magnetic resonance (NMR) | |
| CN104374679A (zh) | 一种确定储层渗透率贡献值的方法 | |
| EP2437078B1 (en) | NMR method for porous media | |
| CN115389387A (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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20201030 |