CN108956417B - 分析岩石孔隙无效吸水量的同位素核磁方法 - Google Patents
分析岩石孔隙无效吸水量的同位素核磁方法 Download PDFInfo
- Publication number
- CN108956417B CN108956417B CN201810590383.0A CN201810590383A CN108956417B CN 108956417 B CN108956417 B CN 108956417B CN 201810590383 A CN201810590383 A CN 201810590383A CN 108956417 B CN108956417 B CN 108956417B
- Authority
- CN
- China
- Prior art keywords
- water
- rock
- spectrum
- core
- sample
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 137
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000010521 absorption reaction Methods 0.000 title claims abstract description 18
- 230000005389 magnetism Effects 0.000 title claims abstract description 8
- 239000011435 rock Substances 0.000 claims abstract description 73
- XLYOFNOQVPJJNP-ZSJDYOACSA-N Heavy water Chemical compound [2H]O[2H] XLYOFNOQVPJJNP-ZSJDYOACSA-N 0.000 claims abstract description 64
- 238000001228 spectrum Methods 0.000 claims abstract description 49
- 239000012530 fluid Substances 0.000 claims abstract description 29
- 230000005311 nuclear magnetism Effects 0.000 claims abstract description 23
- 238000006073 displacement reaction Methods 0.000 claims abstract description 21
- 238000002347 injection Methods 0.000 claims abstract description 20
- 239000007924 injection Substances 0.000 claims abstract description 20
- YZCKVEUIGOORGS-IGMARMGPSA-N Protium Chemical compound [1H] YZCKVEUIGOORGS-IGMARMGPSA-N 0.000 claims abstract description 19
- 241000720974 Protium Species 0.000 claims abstract description 19
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 18
- JEGUKCSWCFPDGT-UHFFFAOYSA-N h2o hydrate Chemical compound O.O JEGUKCSWCFPDGT-UHFFFAOYSA-N 0.000 claims abstract description 16
- 230000008569 process Effects 0.000 claims abstract description 16
- 238000009736 wetting Methods 0.000 claims abstract description 11
- 230000002745 absorbent Effects 0.000 claims abstract description 6
- 239000002250 absorbent Substances 0.000 claims abstract description 6
- 238000009826 distribution Methods 0.000 claims abstract description 6
- 239000011148 porous material Substances 0.000 claims description 17
- 239000007788 liquid Substances 0.000 claims description 16
- 238000004140 cleaning Methods 0.000 claims description 9
- 238000005481 NMR spectroscopy Methods 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000002360 preparation method Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- YZCKVEUIGOORGS-OUBTZVSYSA-N Deuterium Chemical compound [2H] YZCKVEUIGOORGS-OUBTZVSYSA-N 0.000 claims description 4
- 230000004323 axial length Effects 0.000 claims description 4
- 229910052805 deuterium Inorganic materials 0.000 claims description 4
- 238000002474 experimental method Methods 0.000 claims description 4
- 239000001307 helium Substances 0.000 claims description 3
- 229910052734 helium Inorganic materials 0.000 claims description 3
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 claims description 3
- 125000001967 indiganyl group Chemical group [H][In]([H])[*] 0.000 claims description 3
- 239000003960 organic solvent Substances 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 238000005406 washing Methods 0.000 claims description 3
- 238000001035 drying Methods 0.000 claims 1
- 238000002372 labelling Methods 0.000 abstract description 2
- 230000004044 response Effects 0.000 abstract description 2
- 238000010183 spectrum analysis Methods 0.000 abstract description 2
- 239000000523 sample Substances 0.000 description 46
- 238000004458 analytical method Methods 0.000 description 6
- 239000000047 product Substances 0.000 description 5
- 238000005259 measurement Methods 0.000 description 4
- 238000004445 quantitative analysis Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000009792 diffusion process Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 239000007789 gas Substances 0.000 description 3
- 238000004043 dyeing Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 230000005415 magnetization Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 241000256602 Isoptera Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000011549 displacement method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 238000011049 filling Methods 0.000 description 1
- 238000003682 fluorination reaction Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000000608 laser ablation Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011946 reduction process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007921 spray Substances 0.000 description 1
- 230000004936 stimulating effect Effects 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 239000011800 void material Substances 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/0806—Details, e.g. sample holders, mounting samples for testing
-
- 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
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
-
- 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
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
- G01N33/246—Earth materials for water content
-
- 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
- G01N2015/0813—Measuring intrusion, e.g. of mercury
-
- 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
- G01N2015/0866—Sorption
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Immunology (AREA)
- Engineering & Computer Science (AREA)
- Dispersion Chemistry (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明属于油田注水开发技术领域，具体涉及分析岩石孔隙无效吸水量的同位素核磁方法，岩心样品的岩心孔隙中饱和润湿相流体普通水H2O，得到饱和普通水后的岩心T2谱，重新让岩心孔隙中饱和润湿相流体重水D2O，得到岩石T2谱作为基线；向饱和重水的岩心样品中注入氟化油，注入普通水H2O，模拟注水过程，再注入氟化油，分析不动水的含量，并得到残余T2谱；取残余T2谱与岩石的基线T2谱之差，所围范围即为无效孔隙吸水部分，得到无效吸水的孔隙的孔径分布情况。本发明利用普通水和重水的核磁响应信号不同，将驱替、同位素标记和核磁T2谱分析联合在一起，确定岩石不动水的来源，以及无效注水量，有效减少了样本中氕元素的影响。
Description
技术领域
本发明属于油田注水开发技术领域，具体涉及分析岩石孔隙无效吸水量的同位素核磁方法。
背景技术
在在油田注水开发过程中，注入的水进入储层后，一部分水进入微观孔隙内，变成不动水，这部分水为孔隙无效吸水。
油气开发过程中，随着油气的不断采出，地层压力逐渐下降，油气自喷流量逐渐下降。为了补充地层能量，常规需要采用注水方法。在注水开发实际应用中，需要对储层孔隙无效吸水进行定量分析，即注入水多少变成了束缚的不动水。无效吸水量是储层注水开发配注量设计的重要参数。
现在尚无对孔隙无效吸水定量分析的成熟方案。
目前分析储层含水饱和度分布的方案有电阻法，因为含水饱和度不同，岩层的电阻不同，根据电阻和含水饱和度之间的关系，可以根据电阻反算储层含水饱和度。电阻法的缺点：可以确定含水饱和度，但无法确定束缚水，即不动水的饱和度，不能用于分析无效吸水量。
还可以用激光刻蚀或者化学腐蚀方法，在玻璃上腐蚀，形成二维多孔介质。然后使用染色后的水进行驱替分析，同时进行视频拍摄。之后，利用图像分析，确定染色水的分布。这种人工岩心法的缺点：仅能实现二维多孔结构，与实际三维岩心有区别；同时视频图像无法观察到微米级孔隙中的流体。
发明内容
为了解决上述技术问题，针对岩石孔隙微观无效吸水定量分析问题，本发明提供一种分析岩石孔隙无效吸水量的同位素核磁方法，可以定量给出驱替水注入岩石后，转变为不动的无效水的量，解决当前常规方法，无法分析无效吸水具体含量的问题。
具体的技术方案为：
分析岩石孔隙无效吸水量的同位素核磁方法，岩心样品的岩心孔隙中饱和润湿相流体普通水H2O，得到饱和普通水后的岩心T2谱，得到普通水体积与T2谱面积之间的转换关系；重新让岩心孔隙中饱和润湿相流体重水D2O，得到岩石T2谱作为基线；向饱和重水的岩心样品中注入氟化油，恢复了地层原始含水饱和度和含油饱和度状态；注入普通水H2O，模拟注水过程，再注入氟化油，分析不动水的含量，并得到残余T2谱；取残余T2谱与岩石的基线T2谱之差，所围范围即为无效孔隙吸水部分，根据普通水体积与T2谱面积转换系数，得到无效吸水量。
具体的，包括以下步骤：
(1)岩石样品清洗制备
钻取一个圆柱形岩心样品，岩心样品直径D厘米，岩心样品轴长L厘米，体积为V，并对岩心样品两侧端面进行机械抛光；用有机溶剂对岩心样品进行洗油处理，然后对岩心样品进行烘干，完成岩心样品制备；利用氦气方法，测定岩心样品孔隙度为孔隙体积为
(2)饱和普通水H2O
将清洗完毕的岩心样品放在普通水H2O液体中，静置不少于24小时，让岩心孔隙中饱和润湿相流体普通水H2O；
(3)流体体积和T2谱面积比例关系标定
将步骤(2)中饱和普通水H2O后的岩心样品放在核磁驱替仪器磁体仓岩心夹持器中，核磁共振仪器外磁场B0调为0.6T，交变磁场B1的频率调整为氕的共振频率，此时得到饱和普通水后的岩心T2谱，其面积为TS；根据TS和岩心孔隙体积VP，确定核磁信号面积和普通水体积之间的转换系数T2W＝VP/TS；
(4)饱和重水D2O
将岩心样品按照步骤(1)再次进行清洗处理，清洗完毕的岩心样品放在重水D2O液体中，静置不少于24小时，让岩心孔隙中饱和润湿相流体重水D2O；
(5)T2谱基线标定
将步骤(4)中饱和重水后的样品放在核磁驱替仪器磁体仓岩心夹持器中，核磁共振仪器外磁场B0调为0.6T，交变磁场B1的频率调整为氕的共振频率，此时得到的岩石T2谱作为基线；由于氘在氕的共振频率下，不会出现T2信号，因此此时接受到的T2谱主要来自于岩石固体部分中氕发出的信号；
(6)注入氟化油
以油藏充注时的压差Pcd向饱和重水的岩心样品中注入氟化油，同时在出口安置一个计量器，计量出口被氟化油驱出的重水体积VD，当VD不再改变时，停止注入氟化油；此时即认为恢复了地层原始含水饱和度和含油饱和度状态；
(7)注入普通水H2O，模拟注水驱替过程
步骤(6)之后，紧接以恒速V1向岩心样品中注入液体，所述的注入液体为普通水H2O，模拟注水过程，同时在出口安置一个计量器，计量出口被普通水H2O驱除的氟化油的体积VO，当注入普通水体积为NPV倍于岩石孔隙总体积时，停止注入普通水H2O；NPV＝注入液体体积/岩石孔隙总体积；
(8)注入氟化油，分析不动水的含量
步骤(7)之后，紧接以恒速V1向岩心样品中注入氟化油，同时在出口安置一个计量器，计量出口被氟化油驱出的普通水H2O体积VH，当VH不再改变时，停止注入氟化油；此时即认为岩心中剩余的水为不动的水；
(9)岩心孔隙无效吸水量确定
步骤(8)实验后岩心样品中残留的普通水即是无效注水；步骤(8)之后，将核磁共振仪器外磁场B0调为0.6T，交变磁场B1的频率调整为氕的共振频率；由于重水在氕的共振频率条件下，无核磁回波信号，因此得到岩心注油步骤(8)后的T2谱信号为残余普通水的信号，称之为残余T2谱；将残余T2谱与步骤(5)中岩石的基线T2谱之差，得到岩心注水驱替无效吸水T2谱图，所围范围即为无效孔隙吸水部分，该部分谱面积为WTS；根据步骤(3)流体体积和T2谱面积比例关系标定得到的比例系数T2W，得到注入的普通水转变为不动无效水的具体体积Vfa＝WTS*T2W；同时根据T2谱弛豫时间和孔隙半径之间的转换关系，从岩心注水驱替无效吸水T2谱图中，得到无效吸水的孔隙的孔径分布情况。
本发明提供的分析岩石孔隙无效吸水量的同位素核磁方法，利用普通水和重水的核磁响应信号不同，将驱替、同位素标记和核磁T2谱分析联合在一起，确定岩石不动水的来源，以及无效注水量；先测定样品基线T2谱，将残余T2谱与基线T2谱的差值作为最终无效吸水的T2谱部分，有效减少了样本中氕元素的影响。
附图说明
图1为岩石核磁T2谱示意图；
图2为核磁驱替仪器***；
图3为实施例的岩心孔隙饱和普通水后测得的T2谱
图4为实施例的岩石T2谱线基线示意图；
图5为实施例的岩石注水驱替无效吸水T2谱示意图。
具体实施方式
结合实施例说明本发明的具体技术方案。
分析大庆扶余油层砂岩的无效吸水量。岩心样品加工为2.5cm直径，6cm轴长的样品。
本方案中所使用的核磁驱替仪器装置如图2，包括：空气压缩机、恒速恒压双缸泵、驱替液仓、氟化油仓、磁体箱、射频装置等。磁体箱提供稳恒外场B0，射频装置包括交变信号B1，以及核磁共振回波信号检测探头。
该仪器使用步骤如下：
(1)将岩心样品放在磁体箱的样品仓岩心夹持器中。
(2)调节磁体箱高温高压控制***，岩心夹持器温度和压力达到地层温压力条件。
(3)按照恒温或者恒压方式，向岩心中压入液体。
(4)打开射频装置，收集岩心中样品的T2谱信号。
本发明的分析岩石孔隙无效吸水量的同位素核磁方法，利用氢元素有两个同位素氕H和氘D，并且在相同的外磁场强度和激发频率下，两种同位素水分子H2O和D2O的核磁信号响应不同的特点，联合驱替方法，定量分析岩石孔隙无效吸水的含量。
其中，核磁共振原理为：原子核是带正电荷的粒子，不能自旋的核没有磁矩，能自旋的核有循环的电流，会产生磁场，形成磁矩(μ)。当自旋核处于磁场强度为B0的外磁场中时，除自旋外，还会绕B0运动，这种运动情况与陀螺的运动情况十分相象，称为拉莫尔进动。自旋核进动的角速度ω0与外磁场强度B0成正比，比例常数即为磁旋比γ：ω0＝2πv0＝γB0，式中v0是进动频率。
微观磁矩在外磁场中的取向是量子化的，自旋量子数为I的原子核在外磁场作用下只可能有2I+1个取向，每一个取向都可以用一个自旋磁量子数m来表示，m与I之间的关系是：
m＝I，I-1，I-2…-I
原子核的每一种取向都代表了核在该磁场中的一种能量状态。
让处于外磁场B0中的自旋核接受垂直于B0的交变磁场B1的激发，当B1辐射的能量恰好等于自旋核两种不同取向的能量差时，处于低能态E1的自旋核吸收电磁辐射能跃迁到高能态E2。这种现象称为核磁共振，简称NMR。
T1弛豫时间和T2弛豫时间：当激励交变磁场B1撤掉后，自旋核在B0的作用下，会从高能态E2的核会自动跃迁到低能态E1，这个过程称之为弛豫。在这个过程中，低能态E1上的核数目逐渐增多，这个增多的过程，称之为纵向弛豫，其增加过程所花费的弛豫时间为T1。同时高能态E2上的核自旋方向也会发生横向变化，但能量不变的过程称之为横向弛豫，横向弛豫导致激励磁场散相，其减少过程所花费的弛豫时间为T2。其中横向是指垂直于磁场B0的所有水平方向。纵上弛豫是一种能量弛豫现象，是能量熵增加过程；横向弛豫是一种几何相位弛豫现象，是一种几何熵增加的过程。
多孔介质：内部包含复杂孔喉结构的介质，如砂岩、页岩等。
多孔介质中流体弛豫时间T2与孔径r的关系：
多孔介质中流体原子核的横向相位弛豫时间T2的决定机制如下：
其中为横向体积弛豫项，T2B为流体自由弛豫时间；为横向表面弛豫项，S为孔隙表面积，V为孔隙体积，ρ横向弛豫强度因子；为扩散弛豫项，G为磁场梯度，TE为回波间隔，γ为磁旋比，D为扩散系数。
当孔隙中仅有一种流体时，T2B要远远大于T2，因此横向体积弛豫项可以忽略；当磁场足够均匀，磁场梯度非常小，回波间隔TE很小的时候，扩散弛豫项也可以忽略。上述条件下，横向弛豫时间只与横向表面弛豫相关：其中F为形状因子，r为孔隙半径。上述表明，小孔径中的流体，相位发生横向弥散的速度越快，横向弛豫时间越短，这部分孔隙中的流体的磁化强度首先减少；大孔径中的流体，相位发生横向弥散的速度越慢，横向弛豫时间越长，这部分孔隙中的流体的磁化强度最后变弱。
岩石核磁T2谱：由于多孔介质的孔隙半径存在一个分部F(r),而不同的半径r的弛豫时间T2不同，因此多孔介质的横向弛豫时间不是一个单值，而是一个分部—T2分部。核磁脉冲序列测量记录得到的自旋回波串按指数衰减，即各单指数衰减叠加：
其中Aj为测量得到的第j个回波信号幅度，ΔTCE是测量回波间隔，T2i为第i种孔径流体的横向弛豫时间常数，Ei为第i种孔隙的横向弛豫强度。回波信号的多指数拟合，就是预设多个T2(T2i)数值，然后找多个Ei，使其可以拟合实测回波信号，T2i和Ei组成核磁T2谱，如图1，其中T2i对应孔隙半径，Ei对应该半径孔隙的体积大小。T2谱图下面的面积对应于岩石孔隙中流体总体积。
具体步骤如下：
(1)岩石样品清洗制备
钻取一个圆柱形岩心样品，岩心样品直径D厘米，岩心样品轴长L厘米，体积为V，并对岩心样品两侧端面进行机械抛光；用有机溶剂对岩心样品进行洗油处理，然后对岩心样品进行烘干，完成岩心样品制备；利用氦气方法，测定岩心样品孔隙度为孔隙体积为
(2)饱和普通水H2O
将清洗完毕的岩心样品放在普通水H2O液体中，静置不少于24小时，让岩心孔隙中饱和润湿相流体普通水H2O；
(3)流体体积和T2谱面积比例关系标定
将步骤(2)中饱和普通水H2O后的岩心样品放在核磁驱替仪器，如图2所示，磁体仓岩心夹持器中，核磁共振仪器外磁场B0调为0.6T，交变磁场B1的频率调整为氕的共振频率，此时得到饱和普通水后的岩心T2谱，如图3所示，其面积为TS；根据TS和岩心孔隙体积VP，确定核磁信号面积和普通水体积之间的转换系数T2W＝VP/TS；
(4)饱和重水D2O
将岩心样品按照步骤(1)再次进行清洗处理，清洗完毕的岩心样品放在重水D2O液体中，静置不少于24小时，让岩心孔隙中饱和润湿相流体重水D2O；
(5)T2谱基线标定
将步骤(4)中饱和重水后的样品放在核磁驱替仪器，磁体仓岩心夹持器中，核磁共振仪器外磁场B0调为0.6T，交变磁场B1的频率调整为氕的共振频率，此时得到的岩石T2谱作为基线，如图4所示；由于氘在氕的共振频率下，不会出现T2信号，因此此时接受到的T2谱主要来自于岩石固体部分中氕发出的信号；
(6)注入氟化油
以油藏充注时的压差Pcd向饱和重水的岩心样品中注入氟化油，同时在出口安置一个计量器，计量出口被氟化油驱出的重水体积VD，当VD不再改变时，停止注入氟化油；此时即认为恢复了地层原始含水饱和度和含油饱和度状态；
(7)注入普通水H2O，模拟注水驱替过程
步骤(6)之后，紧接以恒速V1向岩心样品中注入液体，所述的注入液体为普通水H2O，模拟注水过程，同时在出口安置一个计量器，计量出口被普通水H2O驱除的氟化油的体积VO，当注入普通水体积为NPV倍于岩石孔隙总体积时，停止注入普通水H2O；NPV＝注入液体体积/岩石孔隙总体积；
(8)注入氟化油，分析不动水的含量
步骤(7)之后，紧接以恒速V1向岩心样品中注入氟化油，同时在出口安置一个计量器，计量出口被氟化油驱出的普通水H2O体积VH，当VH不再改变时，停止注入氟化油；此时即认为岩心中剩余的水为不动的水；
(9)岩心孔隙无效吸水量确定
步骤(8)实验后岩心样品中残留的普通水即是无效注水；步骤(8)之后，将核磁共振仪器外磁场B0调为0.6T，交变磁场B1的频率调整为氕的共振频率；由于重水在氕的共振频率条件下，无核磁回波信号，因此得到岩心注油步骤(8)后的T2谱信号为残余普通水的信号，称之为残余T2谱；将残余T2谱与步骤(5)中岩石的基线T2谱之差，得到岩心注水驱替无效吸水T2谱图，如图5所示，所围范围即为无效孔隙吸水部分，该部分谱面积为WTS；根据步骤(3)流体体积和T2谱面积比例关系标定得到的比例系数T2W，得到注入的普通水转变为不动无效水的具体体积Vfa＝WTS*T2W；同时根据T2谱弛豫时间和孔隙半径之间的转换关系，从图5岩心注水驱替无效吸水T2谱图中，得到无效吸水的孔隙的孔径分布情况。
大庆扶余油层砂岩的注入水转换为无效吸水的量为0.4ml。而常规电阻法和人造岩心实验方法，无法得到无效吸水量。
Claims (1)
1.分析岩石孔隙无效吸水量的同位素核磁方法，岩心样品的岩心孔隙中饱和润湿相流体普通水H2O，得到饱和普通水后的岩心T2谱，得到普通水体积与T2谱面积之间的转换关系；重新让岩心孔隙中饱和润湿相流体重水D2O，得到岩石T2谱作为基线；向饱和重水的岩心样品中注入氟化油，恢复了地层原始含水饱和度和含油饱和度状态；注入普通水H2O，模拟注水过程，再注入氟化油，分析不动水的含量，并得到残余T2谱；取残余T2谱与岩石的基线T2谱之差，所围范围即为无效孔隙吸水部分，根据普通水体积与T2谱面积转换系数，得到无效吸水量；其特征在于，包括以下步骤：
(1)岩石样品清洗制备
钻取一个圆柱形岩心样品，岩心样品直径D厘米，岩心样品轴长L厘米，体积为V，并对岩心样品两侧端面进行机械抛光；用有机溶剂对岩心样品进行洗油处理，然后对岩心样品进行烘干，完成岩心样品制备；利用氦气方法，测定岩心样品孔隙度为孔隙体积为
(2)饱和普通水H2O
将清洗完毕的岩心样品放在普通水H2O液体中，静置不少于24小时，让岩心孔隙中饱和润湿相流体普通水H2O；
(3)流体体积和T2谱面积比例关系标定
将步骤(2)中饱和普通水H2O后的岩心样品放在核磁驱替仪器磁体仓岩心夹持器中，核磁共振仪器外磁场B0调为0.6T，交变磁场B1的频率调整为氕的共振频率，此时得到饱和普通水后的岩心T2谱，其面积为TS；根据TS和岩心孔隙体积VP，确定核磁信号面积和普通水体积之间的转换系数T2W＝VP/TS；
(4)饱和重水D2O
将岩心样品按照步骤(1)再次进行清洗处理，清洗完毕的岩心样品放在重水D2O液体中，静置不少于24小时，让岩心孔隙中饱和润湿相流体重水D2O；
(5)T2谱基线标定
将步骤(4)中饱和重水后的样品放在核磁驱替仪器磁体仓岩心夹持器中，核磁共振仪器外磁场B0调为0.6T，交变磁场B1的频率调整为氕的共振频率，此时得到的岩石T2谱作为基线；由于氘在氕的共振频率下，不会出现T2信号，因此此时接受到的T2谱主要来自于岩石固体部分中氕发出的信号；
(6)注入氟化油
以油藏充注时的压差Pcd向饱和重水的岩心样品中注入氟化油，同时在出口安置一个计量器，计量出口被氟化油驱出的重水体积VD，当VD不再改变时，停止注入氟化油；此时即认为恢复了地层原始含水饱和度和含油饱和度状态；
(7)注入普通水H2O，模拟注水驱替过程
步骤(6)之后，紧接以恒速V1向岩心样品中注入液体，所述的注入液体为普通水H2O，模拟注水过程，同时在出口安置一个计量器，计量出口被普通水H2O驱除的氟化油的体积VO，当注入普通水体积为NPV倍于岩石孔隙总体积时，停止注入普通水H2O；NPV＝注入液体体积/岩石孔隙总体积；
(8)注入氟化油，分析不动水的含量
步骤(7)之后，紧接以恒速V1向岩心样品中注入氟化油，同时在出口安置一个计量器，计量出口被氟化油驱出的普通水H2O体积VH，当VH不再改变时，停止注入氟化油；此时即认为岩心中剩余的水为不动的水；
(9)岩心孔隙无效吸水量确定
步骤(8)实验后岩心样品中残留的普通水即是无效注水；步骤(8)之后，将核磁共振仪器外磁场B0调为0.6T，交变磁场B1的频率调整为氕的共振频率；由于重水在氕的共振频率条件下，无核磁回波信号，因此得到岩心注油步骤(8)后的T2谱信号为残余普通水的信号，称之为残余T2谱；将残余T2谱与步骤(5)中岩石的基线T2谱之差，得到岩心注水驱替无效吸水T2谱图，所围范围即为无效孔隙吸水部分，该部分谱面积为WTS；根据步骤(3)流体体积和T2谱面积比例关系标定得到的比例系数T2W，得到注入的普通水转变为不动无效水的具体体积Vfa＝WTS*T2W；同时根据T2谱弛豫时间和孔隙半径之间的转换关系，从岩心注水驱替无效吸水T2谱图中，得到无效吸水的孔隙的孔径分布情况。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810590383.0A CN108956417B (zh) | 2018-06-09 | 2018-06-09 | 分析岩石孔隙无效吸水量的同位素核磁方法 |
| PCT/CN2018/109269 WO2019233000A1 (zh) | 2018-06-09 | 2018-10-08 | 分析岩石孔隙无效吸水量的同位素核磁方法 |
| US16/398,033 US10801941B2 (en) | 2018-06-09 | 2019-04-29 | Isotope nuclear magnetic method for analyzing ineffective water absorption of rock pores |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810590383.0A CN108956417B (zh) | 2018-06-09 | 2018-06-09 | 分析岩石孔隙无效吸水量的同位素核磁方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108956417A CN108956417A (zh) | 2018-12-07 |
| CN108956417B true CN108956417B (zh) | 2019-09-24 |
Family
ID=64488122
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810590383.0A Active CN108956417B (zh) | 2018-06-09 | 2018-06-09 | 分析岩石孔隙无效吸水量的同位素核磁方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US10801941B2 (zh) |
| CN (1) | CN108956417B (zh) |
| WO (1) | WO2019233000A1 (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2758777C1 (ru) * | 2021-03-22 | 2021-11-01 | Федеральное государственное бюджетное учреждение науки Институт оптики атмосферы им. В.Е. Зуева Сибирского отделения Российской академии наук | Способ измерения размера пор гидрофильных материалов |
Families Citing this family (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10495589B2 (en) * | 2018-04-17 | 2019-12-03 | Saudi Arabian Oil Company | Determining permeability of porous media based on nuclear magnetic resonance measurement |
| CN110410058B (zh) * | 2019-06-20 | 2023-01-31 | 中国石油化工集团有限公司 | 一种校正岩心实验结果刻度二维核磁测井的方法 |
| CN110487835B (zh) * | 2019-09-18 | 2022-04-05 | 西南石油大学 | 一种计算致密油气藏储层饱和度指数的新方法 |
| CN112782477B (zh) * | 2019-11-11 | 2024-05-14 | 中国石油化工股份有限公司 | 一种测定岩心不同润湿状态下电性响应特征的方法及*** |
| CN111537541B (zh) * | 2020-05-26 | 2023-04-07 | 西安石油大学 | 致密油藏co2驱储层动用特征评价方法 |
| CN114062217A (zh) * | 2020-08-07 | 2022-02-18 | 中国石油化工股份有限公司 | 一种饱和岩心带压渗吸实验装置及实验方法 |
| CN112326710B (zh) * | 2020-11-10 | 2022-06-03 | 中国石油大学(华东) | 一种泥页岩饱和水后各离心状态核磁共振t2谱校正方法 |
| CN112858133B (zh) * | 2021-01-12 | 2023-09-12 | 西安石油大学 | 致密油藏裂缝动态渗吸排驱规律评价方法 |
| CN112505085B (zh) * | 2021-02-05 | 2021-04-09 | 西南石油大学 | 基于核磁共振的孔隙度有效应力系数测定方法 |
| CN113533156B (zh) * | 2021-06-30 | 2024-02-02 | 西安石油大学 | 页岩油储集层微观孔隙结构特征及多类孔隙流体识别方法 |
| US11614417B2 (en) | 2021-07-06 | 2023-03-28 | Saudi Arabian Oil Company | Determining saturation in low resistivity pay zones |
| CN113916745B (zh) * | 2021-10-12 | 2024-05-10 | 延安大学 | 一种无损测量水驱气藏微观孔隙结构变化规律的实验方法 |
| CN114088585B (zh) * | 2021-11-15 | 2022-08-12 | 中国石油大学(北京) | 一种页岩润湿角预测模型建立的方法、装置、设备及介质 |
| CN114235641B (zh) * | 2021-12-22 | 2023-10-10 | 常州工学院 | 一种核磁共振测量致密岩石润湿性的方法 |
| CN114417593B (zh) * | 2022-01-13 | 2023-06-06 | 东北石油大学 | 一种评价孔隙连通性对页岩采收率影响的实验方法 |
| CN114486978A (zh) * | 2022-01-26 | 2022-05-13 | 东北石油大学 | 一种定量分析岩心剩余油的方法及*** |
| CN114778588B (zh) * | 2022-04-27 | 2024-05-14 | 长江大学 | 基于nmr表征致密砂岩不同类型流体分布的方法 |
| CN115326645B (zh) * | 2022-08-19 | 2024-04-23 | 中国石油大学(北京) | 基于岩石的基质扩散过程实时动态监测方法 |
Family Cites Families (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2530026A1 (fr) * | 1982-07-09 | 1984-01-13 | Bengsch Eberhard | Procede pour identifier et pour determiner les origines bio- et/ou techno-synthetiques de substances organiques |
| MXPA03010645A (es) * | 2002-12-03 | 2004-06-14 | Schlumberger Technology Bv | Metodo y aparato que utiliza mediciones de resonancia magnetica nuclear (rmn) para reunir informacion sobre una propiedad de la formacion terrestres que rodea un pozo de sondeo. |
| WO2005091018A1 (en) * | 2004-03-18 | 2005-09-29 | Baker Hughes Incorporated | Rock properties prediction, categorization, and recognition from nmr echo-trains using linear and nonlinear regression |
| CN101004134A (zh) * | 2007-01-10 | 2007-07-25 | 辽河石油勘探局 | 利用t2弛豫谱解释稠油储层的方法 |
| EP2019311A1 (en) * | 2007-07-23 | 2009-01-28 | F.Hoffmann-La Roche Ag | A method for detecting a target substance by nuclear magnetic resonance |
| CN101915716B (zh) * | 2010-07-06 | 2012-01-04 | 中国石油天然气股份有限公司 | 一种判断储层岩石的润湿性的方法 |
| CN102944571B (zh) * | 2012-10-17 | 2016-02-03 | 中国地质大学（北京） | 一种测量煤中不同状态水分含量的方法 |
| EP3069127B1 (en) * | 2013-11-13 | 2023-07-26 | NanoNord A/S | A method for quantitative determination of nitrogen in an aqueous fluid |
| CN104316554B (zh) * | 2014-10-27 | 2016-07-06 | 中国石油天然气股份有限公司 | 储层孔隙水可动性测试方法 |
| CN104697915B (zh) * | 2015-03-20 | 2017-12-29 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 一种页岩微观孔隙大小及流体分布的分析方法 |
| CN105223116B (zh) * | 2015-08-28 | 2017-09-01 | 中国石油天然气集团公司 | 一种基于核磁共振谱系数法计算束缚水饱和度的方法 |
| CN106442600B (zh) * | 2016-11-23 | 2017-10-13 | 中国石油大学(华东) | 页岩束缚水含量确定方法 |
| CN107894386B (zh) * | 2017-11-14 | 2019-09-13 | 西安石油大学 | 超临界二氧化碳注入对低渗透砂岩油藏孔喉结构影响的定量评价方法 |
-
2018
- 2018-06-09 CN CN201810590383.0A patent/CN108956417B/zh active Active
- 2018-10-08 WO PCT/CN2018/109269 patent/WO2019233000A1/zh active Application Filing
-
2019
- 2019-04-29 US US16/398,033 patent/US10801941B2/en active Active
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2758777C1 (ru) * | 2021-03-22 | 2021-11-01 | Федеральное государственное бюджетное учреждение науки Институт оптики атмосферы им. В.Е. Зуева Сибирского отделения Российской академии наук | Способ измерения размера пор гидрофильных материалов |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108956417A (zh) | 2018-12-07 |
| US20190257733A1 (en) | 2019-08-22 |
| US10801941B2 (en) | 2020-10-13 |
| WO2019233000A1 (zh) | 2019-12-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108956417B (zh) | 分析岩石孔隙无效吸水量的同位素核磁方法 | |
| Dai et al. | Oil migration in nanometer to micrometer sized pores of tight oil sandstone during dynamic surfactant imbibition with online NMR | |
| Xue et al. | Study of drainage and percolation of nitrogen–water flooding in tight coal by NMR imaging | |
| Guo et al. | Advances in low-field nuclear magnetic resonance (NMR) technologies applied for characterization of pore space inside rocks: a critical review | |
| Elsayed et al. | A review on the applications of nuclear magnetic resonance (NMR) in the oil and gas industry: laboratory and field-scale measurements | |
| US4424487A (en) | Dispersion coefficient determination | |
| Cai et al. | Monitoring oil displacement and CO2 trapping in low-permeability media using NMR: A comparison of miscible and immiscible flooding | |
| CN107561112B (zh) | 一种获取岩石渗透率剖面的核磁共振方法 | |
| US20060116828A1 (en) | Methods and apparatus for measuring capillary pressure in a sample | |
| Lebon et al. | Experimental measurement of dispersion processes at short times using a pulsed field gradient NMR technique | |
| Elsayed et al. | The effect of clay content on the spin–spin NMR relaxation time measured in porous media | |
| Zheng et al. | Microscopic mechanism of fracturing fluid imbibition in stimulated tight oil reservoir | |
| Maojin et al. | Nuclear magnetic resonance simulations of nano-scale cores and microscopic mechanisms of oil shale | |
| CN105940185A (zh) | 使用核磁共振进行的井下流体监测 | |
| Elsayed et al. | New technique for evaluating fracture geometry and preferential orientation using pulsed field gradient nuclear magnetic resonance | |
| Lev et al. | Opportunities for characterizing geological flows using magnetic resonance imaging | |
| Amin et al. | Study of infiltration into a heterogeneous soil using magnetic resonance imaging | |
| Wang et al. | Ionic strength-dependent pre-asymptoic diffusion coefficient distribution in porous media-determination through the pulsed field gradient technique | |
| Geng et al. | A novel method for quantitative analysis of brine and oil saturation in a sealed low permeable core based on one-dimensional (1D) NMR technique | |
| Fleury | Measurement of interfacial area from NMR time dependent diffusion and relaxation measurements | |
| Montague et al. | Visualization and simulation of density driven convection in porous media using magnetic resonance imaging | |
| Zhao et al. | Experimental investigation on spontaneous counter‐current imbibition in water‐wet natural reservoir sandstone core using MRI | |
| Chen et al. | Quantitative magnetic resonance imaging methods for core analysis | |
| Gidcumb | Characterization of methane storage in gas shales using NMR | |
| CN107807078B (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 |