CN108827853B - 基于核磁共振的致密储层岩电测量装置及测量方法 - Google Patents
基于核磁共振的致密储层岩电测量装置及测量方法 Download PDFInfo
- Publication number
- CN108827853B CN108827853B CN201810492350.2A CN201810492350A CN108827853B CN 108827853 B CN108827853 B CN 108827853B CN 201810492350 A CN201810492350 A CN 201810492350A CN 108827853 B CN108827853 B CN 108827853B
- Authority
- CN
- China
- Prior art keywords
- rock sample
- pressure
- valve
- rock
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 198
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 55
- 238000005259 measurement Methods 0.000 title claims abstract description 17
- 238000000691 measurement method Methods 0.000 title abstract description 4
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims abstract description 34
- 229910052757 nitrogen Inorganic materials 0.000 claims abstract description 17
- 238000003860 storage Methods 0.000 claims abstract description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 94
- 229920006395 saturated elastomer Polymers 0.000 claims description 28
- 239000011148 porous material Substances 0.000 claims description 26
- 238000000034 method Methods 0.000 claims description 24
- 239000007789 gas Substances 0.000 claims description 17
- 238000001228 spectrum Methods 0.000 claims description 16
- 238000005303 weighing Methods 0.000 claims description 8
- 239000008398 formation water Substances 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 238000009738 saturating Methods 0.000 claims description 3
- 238000007789 sealing Methods 0.000 claims description 3
- 239000011800 void material Substances 0.000 claims description 3
- 238000009736 wetting Methods 0.000 claims description 3
- 230000005311 nuclear magnetism Effects 0.000 claims description 2
- 238000006073 displacement reaction Methods 0.000 description 9
- 238000009826 distribution Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000005192 partition Methods 0.000 description 5
- 238000003825 pressing Methods 0.000 description 4
- 238000013480 data collection Methods 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 231100000614 poison Toxicity 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000003440 toxic substance 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/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- 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
Abstract
本发明公开了一种基于核磁共振的致密储层岩电测量装置及测量方法，其中测量装置其包括的毛管压力电性联测仪包括高压氮气储罐、围压泵和用于装夹岩样的岩心夹持器，高压氮气储罐和围压泵均通过管道与岩心夹持器一端连接，岩心夹持器的另一端上连接的管道延伸至测量瓶内；岩心夹持器放置于核磁共振仪的测量腔内，高压氮气储罐与岩心夹持器之间的管道上设置有第一阀门和第一压力控制器，围压泵与岩心夹持器之间的管道上设置有第二阀门和第二压力控制器；岩心夹持器与测量瓶之间的管道上设置有第三阀门；岩心夹持器的两端分别通过一个电极与LCR数字电桥连接，第一压力控制器、第二压力控制器、LCR数字电桥和核磁共振仪均与数据采集控制台连接。
Description
技术领域
本发明涉及致密储层性能研究领域，具体涉及一种基于核磁共振的致密储层岩电测量装置及测量方法。
背景技术
目前，致密储层孔隙结构和岩石成分复杂、储集空间多样、非均质性强，故对其复杂孔喉系统的测井响应一直缺乏系统的研究，使得致密储层评价效果不好，复杂储层测井解释符合率低，难以准确获取含油气饱和度，使用传统的阿尔奇公式也显示出现了非阿尔奇特性。
常用的储层孔喉结构表征技术手段主要有铸体薄片、扫描电镜、毛管压力曲线法(压汞技术)、核磁共振及微纳米-CT扫描技术等。其中铸体薄片与扫描电镜都只能针对某个二维断面进行观察，经过后续图像处理提取有限的二维孔喉结构信息。
毛管压力曲线法最常用的是压汞技术，常规压汞技术不能直接测量吼道数量，只能给出不同吼道半径及对应吼道所控制的体积分布。恒速压汞技术受进汞压力限制，无法识别半径小于0.119μm的孔隙和喉道，且还涉及有毒物质的使用。微纳米-CT扫描法，扫描速度快，扫描覆盖范围大，提供孔喉结构定量参数，但测量方法复杂，且费用较高。
发明内容
针对现有技术中的上述不足，本发明提供的基于核磁共振的致密储层岩电测量装置及测量方法通过测量的参数能够计算致密储层的多个岩石电性参数。
为了达到上述发明目的，本发明采用的技术方案为：
第一方面，提供一种基于核磁共振的致密储层岩电测量装置，其包括毛管压力电性联测仪、核磁共振仪和数据采集控制台；毛管压力电性联测仪包括高压氮气储罐、围压泵和用于装夹岩样的岩心夹持器，高压氮气储罐和围压泵均通过管道与岩心夹持器一端连接，岩心夹持器的另一端上连接的管道延伸至放置在称重装置上的测量瓶内；
岩心夹持器放置于核磁共振仪的测量腔内，高压氮气储罐与岩心夹持器之间的管道上设置有第一阀门和第一压力控制器，围压泵与岩心夹持器之间的管道上设置有第二阀门和第二压力控制器；岩心夹持器与测量瓶之间的管道上设置有第三阀门；
岩心夹持器的两端分别通过一个电极与用于测量岩样电阻的LCR数字电桥连接，第一压力控制器、第二压力控制器、LCR数字电桥、核磁共振仪和称重装置均与数据采集控制台连接。
第二方面，提供一种采用基于核磁共振的致密储层岩电测量装置测量岩石电性参数的方法，其包括以下步骤：
S1，获取致密储层处的烘干岩样，并记录岩样的孔隙度、渗透率、长度、干重和直径，之后配置地层水，并在设定压力下对岩样进行饱和，并测量岩样的饱和重；
S2，当排空测量装置内的气体和水分后，关闭第一阀门、第二阀门和第三阀门，将饱和的岩样放置于岩心夹持器的密封腔内，开启核磁共振仪和LCR数字电桥，并记录岩样完全饱和水时的岩样温度、毛管压力P0、电阻和T2谱；
S3，根据记录的岩样温度、电阻和核磁共振弛豫时间，计算岩样完全饱和水时的电阻率、地层因素和孔隙半径：
rc＝ρFST2 (3)
其中，T为岩样完全饱和水时的岩样温度；Ro为岩样完全饱和水时的电阻率；r0为岩样完全饱和水时的电阻；c为电极系数；F为地层因素；Rw为地层水电阻率；φ为孔隙度；m为胶结指数；a为岩性系数；a1、b1为常数；T2为核磁共振横向弛豫时间；ρ为岩石横向表面弛豫率；FS为孔隙形状因子；rc为孔隙半径；
S4，开启第二阀门，采用围压泵给岩样施加设定围压后，开启第一阀门和第三阀门，当测量装置内的压力达到设置压力时，每隔设定时间记录一次岩样的温度、电阻和T2谱；
S5，当与测量瓶连接的管道无水流出时，关闭第二阀门和第三阀门，并记录岩样处于第一个点含水饱和度下测量瓶的重量及毛管压力、岩样温度、电阻和T2谱；
S6，打开第二阀门和第三阀门，继续采用围压泵给岩样施加设定围压，并每隔设定时间记录一次岩样的温度、电阻和核磁共振弛豫时间；
当与测量瓶连接的管道无水流出时，关闭第二阀门和第三阀门，记录岩样处于第二个点含水饱和度下测量瓶的重量及毛管压力、岩样温度、电阻和核磁共振弛豫时间；
S7，重复步骤S6，得到第三个点含水饱和度下测量瓶的重量及毛管压力、岩样温度、电阻和核磁共振弛豫时间，之后关闭第一阀门、第二阀门和第三阀门，取出岩样；
S8，分别计算岩样在第一个点含水饱和度、第二个点含水饱和度和第三个点含水饱和度时的电阻率和孔隙半径：
rcx＝ρFST2x (6)
其中，Rtx为岩样在第x个点含水饱和度下的电阻率；rx为岩样在第x个点含水饱和度下的电阻；Tx为岩样在第x个点含水饱和度下的温度；rcx为岩样在第x个点含水饱和度下的孔隙半径；T2x为岩样在第x个点含水饱和度下的核磁共振横向弛豫时间；Swpx为岩样的第x个点含水饱和度；Sw0为岩样的初始含水饱和度；mx为岩样在第x个点含水饱和度下的测量瓶的重量；ρw为水的密度；VP为岩样的空隙体积；m'1为岩样的饱和重；m'0为岩样的干重；
S9，选取多根不同的岩样，重复步骤S1至S3，采用不同岩样得到的地层因素计算公式和胶结指数的计算公式计算岩性系数a和胶结指数m；
S10，根据同一个岩样的多个点含水饱和度及该含水饱和度下的电阻率，计算岩样的饱和度指数和岩性系数：
其中，n为岩样的饱和度指数；b为岩样岩性系数；RI为电阻率增大系数。
本发明的有益效果为：本方案将饱和的岩样放置在岩心夹持器中，采用半渗隔板法和核磁共振仪相结合，通过高压氮气瓶和围压泵对饱和岩样施加一定围压和驱替压力，在无水流出时，通过LCR数字电桥和核磁共振仪可以测得驱替过程中不同含水饱和度下的岩心电阻和T2谱及毛管压力；
数据采集控制台通过岩心电阻和T2谱及毛管压力可以对整个测量装置进行监测，实时测量不同含水饱和度下岩样的电阻率、毛管压力和孔隙半径，从而对岩样的孔喉分布及含油气饱和度进行有效的评价及分析。
附图说明
图1为基于核磁共振的致密储层岩电测量装置的示意图。
其中，1、第一阀门；2、第二阀门；3、第三阀门；4、高压氮气储罐；5、第一压力控制器；6、岩心夹持器；7、核磁共振仪；8、LCR数字电桥；9、电极；11、温度采集模块；12、时间控制器；13、围压泵；14、第二压力控制器；15、测量瓶；16、称重装置；17、压力采集器；18、数据采集控制台；19、亲水隔板。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
如图1所示，该基于核磁共振的致密储层岩电测量装置包括毛管压力电性联测仪、核磁共振仪7和数据采集控制台18；毛管压力电性联测仪包括高压氮气储罐4、围压泵13和用于装夹岩样的岩心夹持器6，高压氮气储罐4和围压泵13均通过管道与岩心夹持器6一端连接，岩心夹持器6的另一端上连接的管道延伸至放置在称重装置16上的测量瓶15内。
岩心夹持器6放置于核磁共振仪7的测量腔内，高压氮气储罐4与岩心夹持器6之间的管道上设置有第一阀门1和第一压力控制器5，围压泵13与岩心夹持器6之间的管道上设置有第二阀门2和第二压力控制器14；岩心夹持器6与测量瓶15之间的管道上设置有第三阀门3。
其中的称重装置16为电子天平，通过岩样相邻饱和水状态时测得的重量差值可以得到岩样排出的水，通过排出水的质量可以快速计算出岩样此时的含水饱和度。
其中，第一压力控制器5和第二压力控制器14实则为一压力传感器，其可以选用型号为PT124G-128的压力传感器，数据采集控制台18可以为一台电脑，也可以选用控制芯片，其型号为TMS320DSC2X。
如图1所示，岩心夹持器6上设置有与数据采集控制台18连接的温度采集模块11，温度采集模块11可以选用温度计或温度传感器，当为温度传感器时，其可以选用型号为DS18B20数字温度传感器。
岩心夹持器6的两端分别通过一个电极9与用于测量岩样电阻的LCR数字电桥8连接，第一压力控制器5、第二压力控制器14、LCR数字电桥8、核磁共振仪7和称重装置16均与数据采集控制台18连接。
本方案提到的第一阀门1、第二阀门2和第三阀门3可以为普通的手动开启关闭的阀门，不过为了方便实现自动化控制，也可以选用可以自动调节的电磁阀，不过此时第一阀门1、第二阀门2和第三阀门3均需要与数据采集控制台18连接。
本测量装置在使用时，高压氮气储罐4和第一阀门1相配合用于给岩心夹持器6提供驱替压力，第一压力控制器5用于采集相应管道上的压力，数据采集控制台18通过其采集的压力数据判断给岩心夹持器6提供的驱替压力是否达到设定值。
围压泵13与第二阀门2相配合用户给岩心夹持器6提供围压，第二压力控制用于采集相应管道上的压力，数据采集控制台18通过其采集的压力数据判断给岩心夹持器6提供的围压是否达到设定值。
当与测量瓶15连接的管道无水流出时，表明岩样处于饱和水状态，此时第一压力控制器5采集的压力为毛管压力，LCR数字电桥8和核磁共振仪7可以随时将其采集的数据上传至数据采集控制台18，在岩样未处于饱和水状态时，数据采集控制台18可以每隔设定时间记录一次，当岩样处于饱和水状态时，数据采集控制台18需要记录该时刻的毛管压力、LCR数字电桥8采集的电阻、核磁共振仪7的T2谱及岩样的温度。
实施时，本方案优选岩心夹持器6内部放置岩心的密封腔的两端放置有与岩样两端接触的亲水隔板19，设置亲水隔板19后可以在突破一定喉道的毛管压力前，避免气体进入岩样中，以保证试验过程中测量的各个数据的准确性。
其中，隔板法的原理是在不超过某一压力的情况下亲水隔板只允许水通过，而不允许气通过。在驱替过程中，只有当外加的驱替压力等于或大于一定喉道的毛管压力时，非润湿相(气)才能通过喉道，进入孔隙中把湿相流体(水)排出。通过测量的出水量可以计算出岩芯含水饱和度，同时LCR数字电桥8可以测量出该饱和度下的岩样电阻，外加压力就相当于一定喉道的毛管压力)。
实施时，在测量瓶15与第三阀门3之间设置有与数据采集控制台18连接的压力采集器17，其可以选用型号为PT124G-128的压力传感器，该压力采集器17的设置，可以通过其采集的压力信号判断岩样中是否有水流出，避免人工观察出现误差引起后续岩石电性参数的计算不准确性。
再次参考图1，实施时，本方案优选基于核磁共振的致密储层岩电测量装置还包括分别与核磁共振仪7、LCR数字电桥8和数据采集控制台18连接的时间控制器12；设置时间控制器12后，可以通过时间控制器12控制核磁共振仪7和LCR数字电桥8每隔设定时间上传一次采集的数据。
本方案提供的测量装置采用半渗隔板法和核磁共振技术设备实时测量不同含水饱和度下岩样的电阻率、毛管压力和孔喉分布，其测量原理如下：
用半渗隔板法测定岩样的电阻率，首先配置地层水对岩样进行饱和，再通过高压氮气瓶和围压泵13对岩样施加一定围压和驱替压力，通过LCR数字电桥8和核磁共振仪7，记录驱替前100％饱和条件和驱替过程中不同含水饱和度下的岩心电阻和T2谱及毛管压力，通过换算公式(1)、(2)(5)和(6)得到不同含水饱和度下的孔隙半径和电阻率，过程中通过时间控制器12设定一定时间(半小时)实行实时监测，并通过数据采集装置进行整个实验装置进行监测，从而对岩样的孔喉分布及含油气饱和度进行有效的评价及分析。
至此已完成基于核磁共振的致密储层岩电测量装置的详细描述，下面接着采用测量装置测量岩石电性参数的方法进行详细说明。
该基于核磁共振的致密储层岩电测量装置测量岩石电性参数的方法包括步骤S1至步骤S10。
在步骤S1中，获取致密储层处的烘干岩样，并记录岩样的孔隙度、渗透率、长度、干重和直径，之后配置地层水，并在设定压力下对岩样进行饱和，并测量岩样的饱和重；
在步骤S2中，当排空测量装置内的气体和水分后，关闭第一阀门1、第二阀门2和第三阀门3，将饱和的岩样放置于岩心夹持器6的密封腔内，开启核磁共振仪7和LCR数字电桥8，并记录岩样完全饱和水时的岩样温度、毛管压力P0、电阻和T2谱。
在对测量装置内的气体和水分是否排空进行判断时，主要是通过第一压力控制器5、第二压力控制器14及压力采集器17相互配合实现，若第一压力控制器5、第二压力控制器14和压力采集器17均无信号输出时，表明测量装置内的气体和水分已排空。
在步骤S3中，根据记录的岩样温度、电阻和核磁共振弛豫时间，计算岩样完全饱和水时的电阻率、地层因素和孔隙半径：
rc＝ρFST2 (3)
其中，T为岩样完全饱和水时的岩样温度，℃；Ro为岩样完全饱和水时的电阻率，Ω·m；r0为岩样完全饱和水时的电阻，Ω；c为电极9系数(c＝1.072)；F为地层因素；Rw为地层水电阻率，Ω·m；φ为孔隙度；m为胶结指数；a为岩性系数；a1、b1为常数；T2为核磁共振横向弛豫时间；ρ为岩石横向表面弛豫率；S/V表示孔隙比表面；FS为孔隙形状因子(对球形孔隙，FS＝3；对柱状喉道，FS＝2)；rc为孔隙半径，μm。
在步骤S4中，开启第二阀门2，采用围压泵13给岩样施加设定围压后，开启第一阀门1和第三阀门3，当测量装置内的压力达到设置压力时，每隔设定时间记录一次岩样的温度、电阻和T2谱。
在步骤S5中，当与测量瓶15连接的管道无水流出时，关闭第二阀门2和第三阀门3，并记录岩样处于第一个点含水饱和度下测量瓶15的重量及毛管压力、岩样温度、电阻和T2谱。
其中，步骤S4和步骤S5中设定围压的判断主要是通过第二压力控制器14采集的气压进行判断，若是第二压力控制器14采集的气压达到设定围压，则关闭围压泵13；设置压力的判断主要通过第一压力控制器5采集的气压进行判断，若第一压力控制器5采集的气压达到设置压力，则关闭高压氮气储罐4。
在步骤S6中，打开第二阀门2和第三阀门3，继续采用围压泵13给岩样施加设定围压，并每隔设定时间记录一次岩样的温度、电阻和核磁共振弛豫时间；
当与测量瓶15连接的管道无水流出时，关闭第二阀门2和第三阀门3，记录岩样处于第二个点含水饱和度下测量瓶15的重量及毛管压力、岩样温度、电阻和核磁共振弛豫时间。
管道无水流出的判断主要是通过压力采集器17采集的信号进行确定，若是压力采集器17不输出信号，则无液体流过，表明管道无水流出。
在步骤S7中，重复步骤S6，得到第三个点含水饱和度下测量瓶15的重量及毛管压力、岩样温度、电阻和核磁共振弛豫时间，之后关闭第一阀门1、第二阀门2和第三阀门3，取出岩样；
在步骤S8中，分别计算岩样在第一个点含水饱和度、第二个点含水饱和度和第三个点含水饱和度时的电阻率和孔隙半径：
rcx＝ρFST2x (6)
其中，Rtx为岩样在第x个点含水饱和度下的电阻率，Ω·m；rx为岩样在第x个点含水饱和度下的电阻，Ω；Tx为岩样在第x个点含水饱和度下的温度；rcx为岩样在第x个点含水饱和度下的孔隙半径；T2x为岩样在第x个点含水饱和度下的核磁共振横向弛豫时间；Swpx为岩样的第x个点含水饱和度；Sw0为岩样的初始含水饱和度；mx为岩样在第x个点含水饱和度下的测量瓶15的重量；ρw为水的密度；VP为岩样的空隙体积；m'1为岩样的饱和重；m'0为岩样的干重；
在步骤S9中，选取多根不同的岩样，重复步骤S1至S3，采用不同岩样得到的地层因素计算公式和胶结指数的计算公式计算岩性系数a和胶结指数m；
在步骤S10中，根据同一个岩样的多个点含水饱和度及该含水饱和度下的电阻率，计算岩样的饱和度指数和岩性系数：
其中，n为岩样的饱和度指数；b为岩样岩性系数；RI为电阻率增大系数。
实施时，本方案优选当致密储层为含气储层或含油储层时，其含气饱和度或含油饱和度的计算公式为：
其中，Sqx为第x个点含气饱和度，Syx为第x个点含油饱和度。
实施时，本方案的测量岩石电性参数的方法还包括根据绘制的毛管压力与电阻率的曲线，计算电阻率幂指数β：
其中，Pcx为第x个点的含水饱和度下的毛管压力，MPa；
测量岩石电性参数的方法还包括构建毛管压力与孔隙半径之间的函数关系：
其中，σ为流体界面张力；θ为润湿接触角；
根据记录的T2谱和与其对应的电阻率，计算核磁拟合指数nt2：
其中，e为自然对数。
综上所述，本方案通过毛管压力电性联测仪及核磁共振仪7对岩样T2谱分布、毛管压力及岩样电阻进行高精度、高效率、易操作地实时检测，并通过测量岩石电性参数的方法快速地得到评价岩石性能的多个岩石电性参数，从而实现对致密储层含油气饱和度进行有效评价。
Claims (4)
1.一种基于核磁共振的致密储层岩电测量装置测量岩石电性参数的方法，其特征在于，基于核磁共振的致密储层岩电测量装置包括毛管压力电性联测仪、核磁共振仪和数据采集控制台；所述毛管压力电性联测仪包括高压氮气储罐、围压泵和用于装夹岩样的岩心夹持器，所述高压氮气储罐和围压泵均通过管道与岩心夹持器一端连接，所述岩心夹持器的另一端连接的管道延伸至放置在称重装置上的测量瓶内；
所述岩心夹持器放置于所述核磁共振仪的测量腔内，所述高压氮气储罐与岩心夹持器之间的管道上设置有第一阀门和第一压力控制器，所述围压泵与岩心夹持器之间的管道上设置有第二阀门和第二压力控制器；所述岩心夹持器与测量瓶之间的管道上设置有第三阀门；
所述岩心夹持器的两端分别通过一个电极与用于测量岩样电阻的LCR数字电桥连接，所述第一压力控制器、第二压力控制器、LCR数字电桥、核磁共振仪和称重装置均与所述数据采集控制台连接；
基于核磁共振的致密储层岩电测量装置测量岩石电性参数的方法包括以下步骤：
S1，获取致密储层处的烘干岩样，并记录岩样的孔隙度、渗透率、长度、干重和直径，之后配置地层水，并在设定压力下对岩样进行饱和，并测量岩样的饱和重；
S2，当排空测量装置内的气体和水分后，关闭第一阀门、第二阀门和第三阀门，将饱和的岩样放置于岩心夹持器的密封腔内，开启核磁共振仪和LCR数字电桥，并记录岩样完全饱和水时的岩样温度、毛管压力P0、电阻和T2谱；
S3，根据记录的岩样温度、电阻和核磁共振弛豫时间，计算岩样完全饱和水时的电阻率、地层因素和孔隙半径：
rc＝ρFST2 (3)
其中，T为岩样完全饱和水时的岩样温度；Ro为岩样完全饱和水时的电阻率；r0为岩样完全饱和水时的电阻；c为电极系数；F为地层因素；Rw为地层水电阻率；φ为孔隙度；m为胶结指数；a为岩性系数；a1、b1为常数；T2为核磁共振横向弛豫时间；ρ为岩石横向表面弛豫率；FS为孔隙形状因子；rc为孔隙半径；
S4，开启第二阀门，采用围压泵给岩样施加设定围压后，开启第一阀门和第三阀门，当测量装置内的压力达到设置压力时，每隔设定时间记录一次岩样的温度、电阻和T2谱；
S5，当与测量瓶连接的管道无水流出时，关闭第二阀门和第三阀门，并记录岩样处于第一个点含水饱和度下测量瓶的重量及毛管压力、岩样温度、电阻和T2谱；
S6，打开第二阀门和第三阀门，继续采用围压泵给岩样施加设定围压，并每隔设定时间记录一次岩样的温度、电阻和核磁共振弛豫时间；
当与测量瓶连接的管道无水流出时，关闭第二阀门和第三阀门，记录岩样处于第二个点含水饱和度下测量瓶的重量及毛管压力、岩样温度、电阻和核磁共振弛豫时间；
S7，重复步骤S6，得到第三个点含水饱和度下测量瓶的重量及毛管压力、岩样温度、电阻和核磁共振弛豫时间，之后关闭第一阀门、第二阀门和第三阀门，取出岩样；
S8，分别计算岩样在第一个点含水饱和度、第二个点含水饱和度、第三个点含水饱和度时的电阻率和孔隙半径及每个点的含水饱和度：
rcx＝ρFST2x (6)
其中，Rtx为岩样在第x个点含水饱和度下的电阻率；rx为岩样在第x个点含水饱和度下的电阻；Tx为岩样在第x个点含水饱和度下的温度；rcx为岩样在第x个点含水饱和度下的孔隙半径；T2x为岩样在第x个点含水饱和度下的核磁共振横向弛豫时间；Swpx为岩样的第x个点含水饱和度；Sw0为岩样的初始含水饱和度；mx为岩样在第x个点含水饱和度下的测量瓶的重量；ρw为水的密度；VP为岩样的空隙体积；m'1为岩样的饱和重；m'0为岩样的干重；
S9，选取多根岩样，重复步骤S1至S3，采用多根岩样得到的地层因素的计算公式和胶结指数的计算公式计算岩性系数a和胶结指数m；
S10，根据同一个岩样的多个点含水饱和度及该含水饱和度下的电阻率，计算岩样的饱和度指数和岩性系数：
其中，n为岩样的饱和度指数；b为岩样岩性系数；RI为电阻率增大系数。
2.根据权利要求1所述的方法，其特征在于，当致密储层为含气储层或含油储层时，其含气饱和度或含油饱和度的计算公式为：
其中，Sqx为第x个点含气饱和度，Syx为第x个点含油饱和度。
3.根据权利要求1所述的方法，其特征在于，还包括根据绘制的毛管压力与电阻率的曲线，计算电阻率幂指数β：
其中，Pcx为第x个点的含水饱和度下的毛管压力。
4.根据权利要求1所述的方法，其特征在于，还包括构建毛管压力与孔隙半径之间的函数关系：
其中，σ为流体界面张力；θ为润湿接触角；
根据记录的T2谱和与其对应的电阻率，计算核磁拟合指数nt2：
其中，e为自然对数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810492350.2A CN108827853B (zh) | 2018-05-22 | 2018-05-22 | 基于核磁共振的致密储层岩电测量装置及测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810492350.2A CN108827853B (zh) | 2018-05-22 | 2018-05-22 | 基于核磁共振的致密储层岩电测量装置及测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108827853A CN108827853A (zh) | 2018-11-16 |
| CN108827853B true CN108827853B (zh) | 2024-02-06 |
Family
ID=64148941
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810492350.2A Active CN108827853B (zh) | 2018-05-22 | 2018-05-22 | 基于核磁共振的致密储层岩电测量装置及测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108827853B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109932301A (zh) * | 2019-04-10 | 2019-06-25 | 西南石油大学 | 一种计算致密储层自发渗吸两相流体相对渗透率的方法 |
| CN110487835B (zh) * | 2019-09-18 | 2022-04-05 | 西南石油大学 | 一种计算致密油气藏储层饱和度指数的新方法 |
| CN110618170B (zh) * | 2019-09-26 | 2021-04-30 | 中国石油天然气股份有限公司 | 岩土的饱和度指数的确定方法及装置 |
| CN112782477A (zh) * | 2019-11-11 | 2021-05-11 | 中国石油化工股份有限公司 | 一种测定岩心不同润湿状态下电性响应特征的方法及系统 |
| CN112881472B (zh) * | 2021-01-18 | 2021-12-14 | 克拉玛依市昂科能源科技有限公司 | 流体驱替过程中岩样电阻率和核磁共振联测的方法与装置 |
| CN113358529B (zh) * | 2021-06-11 | 2022-04-05 | 西南石油大学 | 一种用于计算砂岩储层弥散系数的装置 |
| CN114441407A (zh) * | 2022-01-14 | 2022-05-06 | 合肥综合性国家科学中心能源研究院(安徽省能源实验室) | 低渗煤岩co2驱替过程动态可视化模拟试验系统与方法 |
Citations (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6008645A (en) * | 1997-03-11 | 1999-12-28 | Conoco Inc. | Prediction of permeability from capillary pressure curves derived from nuclear magnetic resonance pore size distributions |
| CN1914500A (zh) * | 2004-01-26 | 2007-02-14 | 皇家飞利浦电子股份有限公司 | 用于芯片上磁共振波谱分析的方法和器件 |
| CN101907586A (zh) * | 2010-06-11 | 2010-12-08 | 中国石油天然气股份有限公司 | 用于核磁共振测试岩心的高温高压夹持器 |
| CN202256095U (zh) * | 2011-10-26 | 2012-05-30 | 中国石油集团西部钻探工程有限公司 | 半透隔板岩心毛管压力电阻率测量装置 |
| CN102519831A (zh) * | 2011-10-26 | 2012-06-27 | 中国石油集团西部钻探工程有限公司 | 半透隔板岩心毛管压力电阻率测量装置 |
| CN102590284A (zh) * | 2012-02-02 | 2012-07-18 | 西南石油大学 | 一种测量高温高压岩心含水饱和度分布的装置 |
| CN102998322A (zh) * | 2011-09-14 | 2013-03-27 | 中国石油天然气股份有限公司 | 恒定梯度场核磁共振岩样分析方法及仪器 |
| CN103353462A (zh) * | 2013-06-17 | 2013-10-16 | 中国石油大学(华东) | 一种基于核磁共振成像的岩石非均质性定量评价方法 |
| CN103884633A (zh) * | 2014-03-05 | 2014-06-25 | 中国石油天然气股份有限公司 | 一种确定岩石渗透率的方法及装置 |
| CN105353419A (zh) * | 2015-11-25 | 2016-02-24 | 中国地质大学（北京） | 一种基于储集层分类的核磁毛管压力曲线构造方法 |
| CN105606517A (zh) * | 2016-03-25 | 2016-05-25 | 中国地质大学（北京） | 一种结合核磁共振测量低渗透储层相对渗透率的仪器 |
| CN105651963A (zh) * | 2015-12-30 | 2016-06-08 | 中国石油天然气股份有限公司 | 岩心自吸毛管压力曲线测试装置、工作方法及控制器 |
| CN106153662A (zh) * | 2016-06-17 | 2016-11-23 | 北京大学 | 岩心应力敏感性的测量方法 |
| CN106501144A (zh) * | 2016-09-13 | 2017-03-15 | 中国石油大学(华东) | 一种基于核磁共振双截止值的致密砂岩渗透率计算方法 |
| CN106814018A (zh) * | 2017-01-18 | 2017-06-09 | 西南石油大学 | 一种致密岩石气相相对渗透率测量装置及方法 |
| CN206410978U (zh) * | 2017-01-18 | 2017-08-15 | 西南石油大学 | 一种致密岩石气相相对渗透率测量装置 |
| CN107525720A (zh) * | 2017-08-22 | 2017-12-29 | 成都理工大学 | 一种测试致密储层敏感性的装置及方法 |
| CN108007954A (zh) * | 2017-11-29 | 2018-05-08 | 西南石油大学 | 一种致密岩心束缚水环境建立的实验装置及测试方法 |
| CN208155804U (zh) * | 2018-05-22 | 2018-11-27 | 西南石油大学 | 基于核磁共振的致密储层岩电测量装置 |
| CN209821099U (zh) * | 2019-03-14 | 2019-12-20 | 西南石油大学 | 基于核磁共振的多功能致密气储层动态参数联测装置 |
| CN210264648U (zh) * | 2019-06-10 | 2020-04-07 | 西南石油大学 | 一种多功能的岩心驱替装置 |
| CN112881472A (zh) * | 2021-01-18 | 2021-06-01 | 克拉玛依市昂科能源科技有限公司 | 流体驱替过程中岩样电阻率和核磁共振联测的方法与装置 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2844355B1 (fr) * | 2002-09-11 | 2005-03-25 | Inst Francais Du Petrole | Methode de mesure de la mouillabilite de roches par resonnance magnetique nucleaire |
| FR2920876B1 (fr) * | 2007-09-07 | 2009-12-04 | Inst Francais Du Petrole | Methode de mesure rapide de la saturation et de la resistivite d'un milieu poreux. |
| US9671525B2 (en) * | 2013-03-08 | 2017-06-06 | Saudi Arabian Oil Company | Determining continuous capillary pressure curves for subsurface earth formations using saturation and NMR log data |
| US10156531B2 (en) * | 2015-10-27 | 2018-12-18 | Saudi Arabian Oil Company | Measuring connectivity between different pore types in porous media |
| US10466381B2 (en) * | 2015-12-28 | 2019-11-05 | Baker Hughes, A Ge Company, Llc | NMR logging in formation with micro-porosity by using first echoes from multiple measurements |
| US10466153B2 (en) * | 2016-02-25 | 2019-11-05 | Exxonmobil Upstream Research Company | Coreflood testing system and methods for simultaneous measurement of key core properties |
-
2018
- 2018-05-22 CN CN201810492350.2A patent/CN108827853B/zh active Active
Patent Citations (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6008645A (en) * | 1997-03-11 | 1999-12-28 | Conoco Inc. | Prediction of permeability from capillary pressure curves derived from nuclear magnetic resonance pore size distributions |
| CN1914500A (zh) * | 2004-01-26 | 2007-02-14 | 皇家飞利浦电子股份有限公司 | 用于芯片上磁共振波谱分析的方法和器件 |
| CN101907586A (zh) * | 2010-06-11 | 2010-12-08 | 中国石油天然气股份有限公司 | 用于核磁共振测试岩心的高温高压夹持器 |
| CN102998322A (zh) * | 2011-09-14 | 2013-03-27 | 中国石油天然气股份有限公司 | 恒定梯度场核磁共振岩样分析方法及仪器 |
| CN202256095U (zh) * | 2011-10-26 | 2012-05-30 | 中国石油集团西部钻探工程有限公司 | 半透隔板岩心毛管压力电阻率测量装置 |
| CN102519831A (zh) * | 2011-10-26 | 2012-06-27 | 中国石油集团西部钻探工程有限公司 | 半透隔板岩心毛管压力电阻率测量装置 |
| CN102590284A (zh) * | 2012-02-02 | 2012-07-18 | 西南石油大学 | 一种测量高温高压岩心含水饱和度分布的装置 |
| CN103353462A (zh) * | 2013-06-17 | 2013-10-16 | 中国石油大学(华东) | 一种基于核磁共振成像的岩石非均质性定量评价方法 |
| CN103884633A (zh) * | 2014-03-05 | 2014-06-25 | 中国石油天然气股份有限公司 | 一种确定岩石渗透率的方法及装置 |
| CN105353419A (zh) * | 2015-11-25 | 2016-02-24 | 中国地质大学（北京） | 一种基于储集层分类的核磁毛管压力曲线构造方法 |
| CN105651963A (zh) * | 2015-12-30 | 2016-06-08 | 中国石油天然气股份有限公司 | 岩心自吸毛管压力曲线测试装置、工作方法及控制器 |
| CN105606517A (zh) * | 2016-03-25 | 2016-05-25 | 中国地质大学（北京） | 一种结合核磁共振测量低渗透储层相对渗透率的仪器 |
| CN106153662A (zh) * | 2016-06-17 | 2016-11-23 | 北京大学 | 岩心应力敏感性的测量方法 |
| CN106501144A (zh) * | 2016-09-13 | 2017-03-15 | 中国石油大学(华东) | 一种基于核磁共振双截止值的致密砂岩渗透率计算方法 |
| CN106814018A (zh) * | 2017-01-18 | 2017-06-09 | 西南石油大学 | 一种致密岩石气相相对渗透率测量装置及方法 |
| CN206410978U (zh) * | 2017-01-18 | 2017-08-15 | 西南石油大学 | 一种致密岩石气相相对渗透率测量装置 |
| CN107525720A (zh) * | 2017-08-22 | 2017-12-29 | 成都理工大学 | 一种测试致密储层敏感性的装置及方法 |
| CN108007954A (zh) * | 2017-11-29 | 2018-05-08 | 西南石油大学 | 一种致密岩心束缚水环境建立的实验装置及测试方法 |
| CN208155804U (zh) * | 2018-05-22 | 2018-11-27 | 西南石油大学 | 基于核磁共振的致密储层岩电测量装置 |
| CN209821099U (zh) * | 2019-03-14 | 2019-12-20 | 西南石油大学 | 基于核磁共振的多功能致密气储层动态参数联测装置 |
| CN210264648U (zh) * | 2019-06-10 | 2020-04-07 | 西南石油大学 | 一种多功能的岩心驱替装置 |
| CN112881472A (zh) * | 2021-01-18 | 2021-06-01 | 克拉玛依市昂科能源科技有限公司 | 流体驱替过程中岩样电阻率和核磁共振联测的方法与装置 |
Non-Patent Citations (5)
| Title |
|---|
| Study on the impacts of capillary number and initial water saturation on the residual gas distribution by NMR;Li T 等;《ENERGIES》;第12卷(第4期);2714 * |
| 基于核磁共振测井和毛管压力的储层渗透率计算方法;肖忠祥;肖亮;;原子能科学技术(10);6-9 * |
| 多孔介质中CO2与油（水）两相渗流的MRI研究;赵越超;《中国博士学位论文全文数据库基础科学辑》(第6期);A011-7 * |
| 致密砂岩储层流体参数核磁共振实验研究;郎东江;伦增珉;吕成远;潘伟义;孙爱军;;CT理论与应用研究(02);101-109 * |
| 致密砂岩毛细压力和电性参数及可动流体饱和度研究;雷浩;《，中国优秀硕士学位论文全文数据库工程科技I辑》(第11期);B019-280 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108827853A (zh) | 2018-11-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108827853B (zh) | 基于核磁共振的致密储层岩电测量装置及测量方法 | |
| US10161891B1 (en) | Method for characterizing rock physical characteristics of deeply buried carbonate rocks | |
| CN110346258B (zh) | 一种致密岩石油相相对渗透率测定的方法 | |
| CN208155804U (zh) | 基于核磁共振的致密储层岩电测量装置 | |
| CN208171813U (zh) | 一种多功能的渗透率测试装置 | |
| CN108195732B (zh) | 一种致密岩心渗吸实验装置及渗吸量测试方法 | |
| CN112881472B (zh) | 流体驱替过程中岩样电阻率和核磁共振联测的方法与装置 | |
| EP3593159A1 (en) | Absolute porosity and pore size determination of pore types in media with varying pore sizes using nmr | |
| CN110426336A (zh) | 一种路基土非饱和渗透系数测量系统及其相关测量方法 | |
| CN111650082B (zh) | 一种非饱和土土水特征曲线测量装置 | |
| CN110261280A (zh) | 一种高温高压岩心逆向渗吸在线监测实验装置及实验方法 | |
| CN105910941B (zh) | 基于压力板仪的冻土未冻水含量的测试方法 | |
| CN209821099U (zh) | 基于核磁共振的多功能致密气储层动态参数联测装置 | |
| Lu et al. | Constant flow method for concurrently measuring soil-water characteristic curve and hydraulic conductivity function | |
| Van Der Ploeg et al. | Polymer tensiometers with ceramic cones: Direct observations of matric pressures in drying soils | |
| CN106501151A (zh) | 一种基于渗吸和离子扩散特性的页岩孔径测量装置及方法 | |
| CN109781765A (zh) | 一种计算致密储层束缚水液膜厚度的新方法 | |
| CN108444890B (zh) | 测试中、高渗岩心液测渗透率的非稳态滴定装置及方法 | |
| CN109799177A (zh) | 一种多组岩样低速非达西渗流试验同时测量的装置和方法 | |
| Van Grinsven et al. | Automated in situ measurement of unsaturated soil water flux | |
| CN113218843A (zh) | 一种声电渗等多功能三轴实验系统及方法 | |
| CN208076347U (zh) | 岩石渗透率和压缩系数联测装置 | |
| CN110208497A (zh) | 一种便携式土壤给水度测试仪及测试方法 | |
| EP3612864B1 (en) | Method for determining a representative parameter of a porous sample and related assembly | |
| CN110082388A (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 |