致密岩心的油水毛管压力曲线的测试方法和系统.pdf

《致密岩心的油水毛管压力曲线的测试方法和系统.pdf》由会员分享，可在线阅读，更多相关《致密岩心的油水毛管压力曲线的测试方法和系统.pdf（19页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 202010547211.2 (22)申请日 2020.06.16 (71)申请人 北京大学 地址 100871 北京市海淀区颐和园路5号 申请人 北京斯派克科技发展有限公司 (72)发明人 王为民马德铭朱涛涛 (74)专利代理机构 北京万象新悦知识产权代理 有限公司 11360 代理人 黄凤茹 (51)Int.Cl. G01N 24/08(2006.01) (54)发明名称 致密岩心的油水毛管压力曲线的测试方法 和系统 (57)摘要 本发明公布了一种致密岩心的油水毛管压 力曲线。

2、的测试方法和系统， 包括： 磁体、 梯度线 圈、 样品移动平台、 核磁共振射频线圈、 核磁共振 谱仪、 射频功率放大器、 梯度放大器、 电源、 水冷 机、 包括主机和显示器的计算机、 前置放大器； 梯 度线圈内配备有水冷系统； 通过设计具有较短回 波时间的纯相位编码脉冲序列， 实现采集致密岩 心内的快衰减信号， 结合离心设备获得致密岩心 中的油水两相共存时流体饱和度分布， 用于测试 致密岩心的油水毛管压力曲线， 更适合油藏开发 现场使用。 权利要求书3页 说明书9页 附图6页 CN 111650234 A 2020.09.11 CN 111650234 A 1.一种致密岩心的油水毛管压力曲线的。

3、测试方法， 通过设计具有较短回波时间的纯相 位编码脉冲序列采集致密岩心内的快衰减信号， 结合离心设备获得致密岩心中的油水两相 共存时流体饱和度分布， 用于测试致密岩心的油水毛管压力曲线； 包括如下步骤： 1)将致密岩心进行洗油、 烘干处理， 并测试氦气孔隙度； 2)将致密岩心抽真空后饱和常规模拟地层水， 并进行加压饱和处理， 确保流体饱和度 大于95， 得到饱和常规模拟地层水状态的致密岩心样品； 3)使用岩心标准样对核磁共振设备进行初始化调试： 核磁共振设备的主磁体具有防涡流层， 核磁共振谱仪配有多组数学波形预增强参数， 对短时间内快速变化的梯度电流进行涡流补偿， 使得在纯相位编码序列快速编码。

4、期间， 梯 度线圈无涡流产生； 在核磁共振设备的梯度线圈内布管进行水冷降温， 针对线圈不同部位的发热量不同进 行布管， 对线圈磁场影响小， 散热效率高， 避免用于纯相位序列编码的梯度电流施加在梯度 线圈上聚集热量； 将岩心标准样置入核磁共振设备的核磁共振射频线圈中； 射频线圈在显示与操作平台测试核磁共振的FID信号， 对FID信号进行傅里叶变换得到 FID信号的中心频率， 并将获得的中心频率置入核磁共振参数主频中， 同时通过改变射频脉 冲的能量获得核磁共振设备的90 翻转角参数， 置入核磁共振设备相应参数中； 即完成核磁 共振设备的初始化调试； 4)将步骤2)得到的饱和常规模拟地层水状态的致密。

5、岩心样品使用不产生磁共振信号 的四氟乙烯密封带包裹后， 放入试管； 5)将步骤4)得到的装有岩心样品的试管垂直置入核磁共振射频线圈内； 6)设计纯相位编码序列 -SPI， 将核磁共振设备加载纯相位编码序列； 包括： 激励脉冲采用小角度硬脉冲， 多次小角度硬脉冲激励和信号采集均在阶梯相位编码的 情况下进行； 全部相位编码过程由一维k空间中心变化至正最大编码梯度， 等待弛豫到稳态等待的 时间RD后， 编码梯度再次由k空间中心变化至负最大编码梯度， 完成由双半一维K空间组成1 维K空间填充的过程； 纯相位编码序列的最短回波时间TE表示为公式1： 其中， TRF为小角度硬脉冲激励的时间， tp为相位编。

6、码时间； 全部相位编码过程中的磁共振信号衰减表示为公式2： 其中， M0为平衡态时的磁化矢量， tp为相位编码时间， T2*为有效横向弛豫时间， 为翻转 角度， M0与流体量成正比； 将核磁共振设备加载纯相位编码序列， 在显示与操作平台输入测试样品轴向测试范 围， 参考步骤3)得到的90 翻转角参数， 输入小角度硬脉冲激励的时间参数TRF、 相位编码时 权利要求书 1/3 页 2 CN 111650234 A 2 间参数tp； 核磁共振设备运行序列 -SPRI， 对得到的信号进行傅里叶变换， 得到饱和状态致密岩 心的流体饱和度分布； 7)将步骤5)中的饱和常规模拟地层水状态的致密岩心样品去除四。

7、氟乙烯密封带后烘 干至恒重， 饱和使用氧化氘配置的模拟地层水； 然后将岩心样品置入充满模拟油的离心杯 中， 采用离心机进行油驱水状态的离心操作； 8)将油驱水状态的岩心样品从离心杯中取出， 重复步骤5)6)， 获得油驱水状态致密 岩心沿轴向的油相饱和度分布； 9)根据油水含烃因子的不同， 将致密岩心油驱水状态的油相轴向饱和度分布转换为离 心后水相轴向饱和度分布结果； 10)分别计算得到岩心上不同位置毛细管压力值； 11)以不同位置含水饱和度为横轴， 对应毛细管压力为纵轴， 绘制致密岩心内油水两相 共存时毛管压力曲线； 通过上述步骤， 实现采用核磁共振纯相位编码对致密岩心的油水毛管压力进行测试。。

8、 2.如权利要求1所述致密岩心的油水毛管压力曲线的测试方法， 其特征是， 核磁共振谱 仪具体配有12组数学波形预增强参数对短时间内快速变化的梯度电流进行涡流补偿。 3.如权利要求1所述致密岩心的油水毛管压力曲线的测试方法， 其特征是， 在核磁共振 设备的梯度线圈内布局铜质水管进行水冷降温， 针对线圈不同部位的发热量不同进行布 管。 4.如权利要求1所述致密岩心的油水毛管压力曲线的测试方法， 其特征是， 岩心标准样 的长度大于等于待测致密岩心的长度； 岩心标准样的孔隙度大于等于20。 5.如权利要求1所述致密岩心的油水毛管压力曲线的测试方法， 其特征是， 步骤10)具 体通过公式3计算得到岩心上。

9、不同位置毛细管压力值： 其中， Pc(r)为致密岩心样品不同位置的油水毛管压力， r为致密岩心样品某一位置距离 离心机转轴轴心的径向距离， 为步骤7)中使用的模拟地层水与模拟油之间的密度差， 为离心机转速， r2为致密岩心样品距离离心机转轴轴心最远处的径向距离。 6.一种实现权利要求15所述致密岩心的油水毛管压力曲线的测试方法的核磁共振 纯相位编码测试系统， 包括： 磁体、 梯度线圈、 样品移动平台、 核磁共振射频线圈、 核磁共振 谱仪、 射频功率放大器、 梯度放大器、 电源、 水冷机、 包括主机和显示器的计算机、 前置放大 器； 梯度线圈内配备有水冷系统； 梯度线圈固定在磁体上； 样品移动平。

10、台固定在磁体上； 核磁共振射频线圈固定在样品 移动平台上； 被测试样品置于核磁共振射频线圈中； 水冷机连接到梯度线圈的水冷系统； 梯 度放大器与梯度线圈相连接； 计算机与核磁共振谱仪相连接； 射频功率放大器与核磁共振 射频线圈相连接； 核磁共振射频线圈与前置放大器相连接； 前置放大器与核磁共振谱仪相 连接； 磁体用于提供B0静磁场； 梯度线圈用于提供XYZ三个方向的梯度场； 核磁共振射频线 圈用于产生磁共振B1射频场， 激发并接收样品产生的磁共振信号； 光纤核磁共振谱仪用于 权利要求书 2/3 页 3 CN 111650234 A 3 控制脉冲序列所需的射频发射、 梯度发射、 数据采集滤波、 。

11、增益放大、 正交解调、 K空间排布； 射频功率放大器用于提供脉冲序列所需的射频能量； 梯度放大器用于提供脉冲序列所 需的梯度电流， 最大输出电流150A,最快爬升时间60us,具有过流过压自动保护功能； 电源 用于为梯度放大器提供直流电压； 水冷机连接到梯度线圈的水冷系统， 通过制冷机带走循 环水的温度； 计算机主机用于安装脉冲序列编辑程序、 脉冲序列测试程序、 数据采集分析显 示程序； 计算机显示器用于显示测试信息； 前置放大器用于对射频线圈中样品产生的核磁 共振信号进行放大。 7.如权利要求6所述的核磁共振纯相位编码测试系统， 其特征是， 所述梯度放大器为X、 Y、 Z方向梯度放大器； 和。

12、/或， 所述磁体为U型结构， 均匀区为长300mm*直径250mm的圆柱体， 均匀度30ppm。 8.如权利要求6所述的核磁共振纯相位编码测试系统， 其特征是， 还可包括电子机柜和 操作台， 电子机柜用于放置并固定电子部件， 提供各电子部件的良好接地； 操作台用于放置 计算机， 提供核磁共振操作人员的工作环境。 9.如权利要求6所述的核磁共振纯相位编码测试系统， 其特征是， 所述梯度线圈用于为 测试系统提供三个方向的梯度场； 通过150安培电流时， 产生梯度场130mT/m。 10.如权利要求6所述的核磁共振纯相位编码测试系统， 其特征是， 所述射频功率放大 器的工作频率范围为1MHz-350。

13、MHz， 最大功率350W； 和/或， 所述前置放大器的增益30dB， 工作频率范围为1MHz-60MHz， 噪声系数0.1dB。 权利要求书 3/3 页 4 CN 111650234 A 4 致密岩心的油水毛管压力曲线的测试方法和系统 技术领域 0001 本发明属于核磁共振技术领域， 涉及一种致密岩心的油水毛管压力曲线的测试方 法和系统。 背景技术 0002 致密油作为非常规油气资源， 储量丰富， 大大改变了我国能源结构， 具有非常重要 的战略意义及商业开发前景， 其储层岩心物性参数的测试分析显得尤其重要， 毛管压力曲 线作为岩心的重要物性参数， 在渗流机理研究、 储层开发评价、 储层成藏勘。

14、探等领域起到非 常重要的作用。 0003 目前毛管压力曲线核磁共振的测试方法中是通过常规医用核磁共振成像技术 (Magnetic Resonance Image)结合离心设备对样品进行流体饱和度分布的测试， 并根据离 心力计算公式， 获得样品毛管压力曲线的分布； 此方法仅适用于高孔隙度、 高渗透率的岩 心， 由于高孔高渗的岩心内部流体丰度高， 核磁共振信号较强且弛豫时间长， 核磁共振成像 信号容易采集， 采用常规医用核磁共振成像技术即可满足应用需求。 0004 致密岩心由于其孔隙度小(流体丰度低)、 渗透率低、 黏土含量高， 多以小孔径为主 等物性特征， 导致其被激发出的核磁共振信号强度低， 。

15、衰减快(即弛豫时间短)， 10ms以上的 弛豫组份少， 而常规医用核磁共振成像技术的重要参数回波时间TE(Time to Echo， 回波时 间)较长， 以0.3T核磁共振成像设备为例， 其一维成像脉冲序列见图1， 回波时间TE-3由图中 90 软脉冲-1、 180 软脉冲-2、 采样时间-4、 读出梯度补偿梯度-5共同决定， 最短回波时间 TE-3在根本上受频率编码方法的限制， 即使提高硬件水平， 由于TE与分辨率正相关， 也就决 定了该方法高分辨率与较短TE不能同时兼顾， 。 一维成像时最短回波时间约为15ms， 多用于 肌肉、 脂肪、 组织液等长弛豫特征的核磁共振成像分析， 而致密岩心为。

16、短弛豫， 横向弛豫时 间T2均小于1ms， 故目前已存在的毛管压力曲线核磁共振的测试技术不适用于致密岩心毛 管压力曲线的核磁共振分析， 且目前已存在的毛管压力曲线核磁共振的测试技术难以对石 油工业实际开发生产中岩心内油水两相共存时的毛管压力曲线进行测试分析。 发明内容 0005 为了克服上述现有毛管压力曲线核磁共振测试技术的不足， 本发明提供一种致密 岩心的油水毛管压力曲线的方法和系统， 通过设计具有较短回波时间的纯相位编码脉冲序 列， 实现采集致密岩心内的快衰减信号， 结合离心设备获得致密岩心中的油水两相共存时 流体饱和度分布， 再根据离心力公式， 计算得到油水毛管压力曲线。 0006 本发。

17、明的核心是： 0007 为了使得核磁共振设备能够采集到致密储层岩芯的信号， 本发明采用纯相位编码 技术， 设计具有短回波时间的纯相位编码序列 -SPI( -single-point imaging， 角度的单 点采集序列， 也称为纯相位编码序列)， 如图3。 采用硬脉冲设计以缩短激发时间； 激励脉冲 采用小角度硬脉冲， 多次小角度硬脉冲激励和信号采集均是在阶梯相位编码的情况下进 说明书 1/9 页 5 CN 111650234 A 5 行， 全部相位编码过程由一维k空间中心变化至正最大编码梯度， 等待RD(Relaxation Delay， 弛豫到稳态等待的时间)后编码梯度再次由k空间中心变化。

18、至负最大编码梯度， 完成 由双半一维K空间组成1维K空间填充的过程。 0008 纯相位编码序列的最短回波时间TE由以下公式1给出： 0009 0010 其中， TRF为小角度硬脉冲激励的时间， tp为相位编码时间(Time of Phase Coding)。 使得在较短TE时间内(约100微秒)可以采集到岩心内流体产生的核磁共振信号； 0011 全部相位编码过程中磁共振信号衰减见公式2： 0012 0013 其中， M0为平衡态时的磁化矢量， tp为相位编码时间， T2*为有效横向弛豫时间， 为 翻转角度， M0与流体量成正比； 0014 整个测试过程中， 测试样品呈现为质子密度加权， 测试样。

19、品的磁共振信号不再受 到纵向时间T1和TR的影响， 同时纯相位编码技术在很大程度上不受磁化率引起的磁场变 化、 化学位移和顺磁杂质引起的磁共振信号的影响； 在本发明采用的纯相位编码序列中， 用 于编码的梯度电流在短时间内快速变化， 梯度线圈会产生涡流， 进而干扰样品产生的核磁 共振信号。 为此本发明使用的核磁共振设备的主磁体具有防涡流层， 核磁共振谱仪配有12 组数学波形预增强参数对短时间内快速变化的梯度电流进行涡流补偿， 使得在纯相位编码 序列快速编码期间， 梯度线圈无涡流产生， 不会影响样品产生的核磁共振信号； 0015 在本发明采用的纯相位编码序列中， 用于编码的梯度电流长时间施加在梯度。

20、线圈 上， 会导致梯度线圈积聚热量。 为了避免热量聚集， 温度升高影响信号采集， 本发明在核磁 共振设备的梯度线圈内布局铜质水管进行水冷降温， 针对线圈不同部位的发热量不同进行 布管， 对线圈磁场影响小， 散热效率高； 0016 本发明的测试样品为致密岩心， 核磁共振信号信噪比低， 为了避免信号传输过程 中的信号进一步损耗， 核磁共振谱仪采用光纤传输模式， 致密岩心样品产生的核磁共振信 号经过前置放大器放大后由电信号转换为光信号， 经光纤传输至核磁共振谱仪中； 0017 为了区分测试样品致密岩心中的油与水， 本发明使用氧化氘(D2O)配置的模拟地 层水作为被测试致密岩心的润湿相， 同时使用模拟。

21、油作为被测试致密岩心的非润湿相， 致 密岩心中油水两相共存时， 仅有油相产生核磁共信号。 通过油水不同的含烃因子将离心后 岩心中油相流体饱和度分布换算成水相流体饱和度分布， 进而得到致密岩心内油水毛管压 力曲线， 更适合油藏开发现场使用。 0018 本发明提供的技术方案是： 0019 一种致密岩心的油水毛管压力曲线的测试方法致密岩心， 通过设计具有较短回波 时间的纯相位编码脉冲序列采集致密岩心内的快衰减信号， 结合离心设备获得致密岩心中 的油水两相共存时流体饱和度分布， 用于测试致密岩心的油水毛管压力曲线， 包括如下步 骤： 0020 1)将致密岩心进行洗油、 烘干处理， 并测试氦气孔隙度。 。

22、0021 2)将致密岩心抽真空后饱和常规模拟地层水， 并进行加压饱和处理， 确保流体饱 说明书 2/9 页 6 CN 111650234 A 6 和度95， 得到饱和常规模拟地层水状态的致密岩心样品备用。 0022 3)使用岩心标准样对核磁共振设备进行初始化调试： 0023 首先将标准样置入核磁共振射频线圈中， 要求标准样的长度大于等于待测致密岩 心的长度， 标准样的孔隙度大于等于20。 射频线圈在显示与操作平台测试核磁共振的FID (free induction decays)信号， 对FID信号进行傅里叶变换得到FID信号的中心频率， 并将 获得的中心频率置入核磁共振参数主频中， 同时通过。

23、改变射频脉冲的能量获得核磁共振设 备的90 翻转角参数， 置入核磁共振设备相应参数中。 核磁共振设备初始化调试结束。 0024 本发明使用的核磁共振设备的主磁体具有防涡流层， 核磁共振谱仪配有12组数学 波形预增强参数对短时间内快速变化的梯度电流进行涡流补偿核， 使得在纯相位编码序列 快速编码期间， 梯度线圈无涡流产生。 在本发明采用的纯相位编码序列中， 用于编码的梯度 电流长时间施加在梯度线圈上， 会导致梯度线圈积聚热量。 为了避免热量聚集， 温度升高影 响信号采集， 本发明在核磁共振设备的梯度线圈内布局铜质水管进行水冷降温， 针对线圈 不同部位的发热量不同进行布管， 对线圈磁场影响小， 散。

24、热效率高。 0025 4)将步骤2)得到的饱和常规模拟地层水状态的致密岩心样品使用不产生磁共振 信号的四氟乙烯密封带包裹后， 放入试管。 0026 5)将步骤4)得到的装有样品的试管垂直置入核磁共振射频线圈内。 0027 6)设计纯相位编码序列； 0028 核磁共振设备加载纯相位编码序列， 在显示与操作平台输入测试样品轴向测试范 围(大于标准样的长度)， 参考步骤3)得到的90 翻转角参数， 输入TRF、 tp参数。 核磁共振设备 运行 -SPI序列， 对得到的信号进行傅里叶变换， 得到饱和状态致密岩心的流体饱和度分 布。 0029 其中， TRF为小角度硬脉冲激励的时间， tp为相位编码时间。

25、。 如公式1中， 核磁共振设 备在较短TE时间内(约100微秒)可以采集到测试样品的核磁共振信号； 这样足以采集致密 岩心内的快衰减信号。 0030 7)将步骤5)中的饱和常规模拟地层水状态的致密岩心样品去除四氟乙烯密封带 后烘干至恒重， 饱和使用氧化氘(D2O)配置的模拟地层水， 然后将岩心样品置入充满模拟油 的离心杯中， 离心机进行油驱水状态的离心操作。 0031 8)将油驱水状态的岩心样品从离心杯中取出， 重复步骤5)、 步骤6)， 获得油驱水状 态致密岩心沿轴向的油相饱和度分布。 0032 9)根据油水含烃因子的不同， 将致密岩心油驱水状态的油相轴向饱和度分布转换 为离心后水相轴向饱和。

26、度分布结果。 0033 10)分别计算岩心上不同位置毛细管压力值。 0034 具体计算公式为： 0035 0036 其中， Pc(r)为致密岩心样品不同位置的油水毛管压力， r为致密岩心样品某一位 置距离离心机转轴轴心的径向距离， 为步骤7)中使用的模拟地层水与模拟油之间的密 度差， 为离心机转速， r2为致密岩心样品距离离心机转轴轴心最远处的径向距离。 0037 11)以不同位置含水饱和度为横轴， 对应毛细管压力为纵轴， 绘制致密岩心内油水 说明书 3/9 页 7 CN 111650234 A 7 毛管压力曲线。 0038 通过上述步骤， 实现采用核磁共振纯相位编码对致密岩心的油水毛管压力进。

27、行测 试。 0039 本发明还提供一种核磁共振纯相位编码测试系统(见图二)， 包括： 磁体、 梯度线 圈、 样品移动平台、 核磁共振射频线圈、 核磁共振谱仪、 射频功率放大器、 梯度放大器(X、 Y、 Z 方向)、 电源、 水冷机、 计算机(包括主机和显示器)、 前置放大器； 整个梯度线圈被固定在磁 体上； 样品移动平台被固定在磁体上， 核磁共振射频线圈被固定样品移动平台上； 被测试样 品置于核磁共振射频线圈中； 水冷机连接到梯度线圈的水冷系统； 梯度放大器(X、 Y、 Z方向) 与梯度线圈相连接； 计算机与核磁共振谱仪相连接； 射频功率放大器与核磁共振射频线圈 相连接； 核磁共振射频线圈与前。

28、置放大器相连接； 前置放大器与核磁共振谱仪相连接； 0040 核磁共振纯相位编码测试系统还可包括电子机柜和操作台， 电子机柜用于放置并 固定电子部件， 提供各电子部件的良好接地； 操作台用于放置计算机(包括主机和显示器)， 提供核磁共振操作人员的工作环境； 0041 其中： 0042 磁体用于为测试系统提供B0静磁场， 磁体为U型结构， 均匀区为长300mm*直径 250mm的圆柱体， 均匀度30ppm。 0043 梯度线圈用于为测试系统提供XYZ三个方向的梯度场， 通过150安培电流时， 产 生梯度场130mT/m， 梯度线圈内配备有水冷系统； 0044 核磁共振射频线圈用于产生磁共振B1射。

29、频场， 激发并接收样品产生的磁共振信 号。 0045 光纤核磁共振谱仪用于控制整个系统中脉冲序列的所需射频发射、 梯度发射、 数 据采集滤波、 增益放大、 正交解调、 K空间排布。 0046 射频功率放大器用于提供脉冲序列所需的射频能量， 工作频率范围1MHz-350MHz， 最大功率350W。 0047 X、 Y、 Z方向的梯度放大器用于提供脉冲序列所需的梯度电流， 最大输出电流150A, 最快爬升时间60us,具有过流过压自动保护功能。 0048 电源为X、 Y、 Z梯度放大器提供150V直流电压。 0049 水冷机连接到梯度线圈的水冷系统， 通过制冷机带走循环水的温度。 0050 计算机。

30、主机安装有脉冲序列编辑程序、 脉冲序列测试程序、 数据采集分析显示程 序。 0051 计算机显示器用于显示测试信息。 0052 前置放大器用于对射频线圈中样品产生的核磁共振信号进行放大， 前置放大器增 益30dB,工作频率范围1MHz-60MHz， 噪声系数0.1dB。 0053 与现有技术相比， 本发明的有益效果是： 0054 针对现有技术无法满足致密岩心磁共振信号采集的问题， 本发明方法通过实现具 有较短回波时间的纯相位编码脉冲序列以及相关硬件配置， 解决了致密岩心毛管压力曲线 核磁共振测试的问题。 0055 对比现有技术， 本发明方法被测样品中流体的核磁共振信号呈现为质子密度加 权， 磁。

31、共振信号不被纵向时间T1和TR影响， 并且不受磁化率引起的磁场变化、 化学位移和顺 说明书 4/9 页 8 CN 111650234 A 8 磁杂质引起的磁共振信号的影响， 为致密岩心核磁共振流体饱和度分布测试较为理想的方 法。 0056 本发明方法依据氧化氘(D2O)不包含1H同位素， 无法产生核磁共振信号， 且物性与 普通水相同等特点， 使用氧化氘(D2O)配置模拟地层水溶液， 致密岩心内仅有油相产生核磁 共信号， 从而实现被测样品中油、 水核磁共振共信号的区分， 计算获得油、 水的流体饱和度 分布， 进而得到致密岩心内油水毛管压力曲线， 更适合油藏开发现场使用。 致密岩心。 附图说明 0。

32、057 图1是现有常规0.3T核磁共振成像设备的一维成像脉冲序列； 0058 其中， 190度激发软脉冲； 2180度激发软脉冲； 3回波时间； 4回波数据采 集； 5读出梯度补偿梯度； 6读出梯度。 0059 图2是本发明实施例所提供的核磁共振岩心分析测试系统的结构框图； 0060 其中， 7永磁磁体； 8梯度线圈； 9样品移动平台； 10待测试岩心； 11核磁 共振射频线圈； 12核磁共振电子机柜； 13核磁共振谱仪； 14核磁共振射频功率放大 器； 15核磁共振X方向梯度功率放大器； 16核磁共振Y方向梯度功率放大器； 17核磁共 振Z方向梯度功率放大器； 18梯度放大器供电电源； 19。

33、水冷机； 20计算机； 21操作 台； 22前置放大器。 0061 图3是本发明实施例中使用的纯相位编码脉冲序列 -SPI时序图； 0062 其中， 23小角度激发硬脉冲； 24数据采集； 25正向梯度阶梯变化； 26重复 时间(Time of Repeatation)； 27相位编码时间(Time of Phase Coding)； 28正向编码 梯度变大值； 29弛豫时间(Relaxation Delay)。 0063 图4是本发明实施例中使用的离心杯结构图； 0064 其中， 30密封盖； 31岩心样品； 32穿孔聚四氟乙烯衬套； 33多孔挡板。 0065 图5是本发明实施例中使用的离心。

34、机的工作转盘示意图； 0066 其中， 34离心杯； 35转轴。 0067 图6是本发明实施例中致密岩心饱和水状态的纯相位编码投影测试结果。 0068 图7是本发明实施例中致密岩心离心后的水相换算投影结果。 0069 图8是本发明实施例中致密岩心离心后的氧化氘换算投影结果。 0070 图9是本发明实施例中致密岩心油水毛管压力曲线计算结果。 具体实施方式 0071 下面结合附图， 通过实施例进一步描述本发明， 但不以任何方式限制本发明的范 围。 0072 本发明提供一种致密岩心核磁共振纯相位编码测试油水毛管压力曲线的方法和 系统， 通过设计较短回波时间的纯相位编码脉冲序列采集致密岩心内的快衰减信。

35、号， 结合 离心设备获得致密岩心中的油水两相共存时流体饱和度分布， 再根据离心力公式， 计算得 到油水毛管压力曲线。 0073 为了在梯度编码作用下能够采集到致密岩心产生的快弛豫信号， 本发明硬件上包 括一台具备一维以上强梯度磁场和脉冲序列开发平台功能的核磁共振设备及一台离心机， 说明书 5/9 页 9 CN 111650234 A 9 该核磁共振设备的梯度磁场最高梯度强度不低于300mT/m， 约为常规医用磁共振成像系统 的10倍， 梯度最短爬升时间为10us。 核磁共振设备包括极环、 极板、 磁钢、 磁体、 磁共振谱仪。 0074 具体实施时， 致密岩心核磁共振纯相位编码测试油水毛管压力曲。

36、线的系统结构如 图2所示， 包括的部件具体为： 0075 A.7是核磁共振纯相位编码测试系统的磁体部分， 为测试系统提供B0静磁场， 磁体 为U型结构， 均匀区为长300mm*直径250mm的圆柱体， 均匀度30ppm。 0076 B.8是核磁共振纯相位编码测试系统的梯度线圈部分， 为测试系统提供XYZ三个 方向的梯度场， 通过150安培电流时， 产生梯度场130mT/m， 梯度线圈内配备有水冷系统， 整个梯度线圈被固定在磁体上。 0077 C.9是核磁共振纯相位编码测试系统的样品移动平台部分， 样品移动平台被固定 在磁体上， 核磁共振射频线圈被固定样品移动平台上。 0078 D.10是核磁共。

37、振纯相位编码测试系统的被测试样品， 被放置于核磁共振射频线圈 中。 0079 E.11是核磁共振纯相位编码测试系统的核磁共振射频线圈， 用于产生磁共振B1射 频场， 激发并接收样品产生的磁共振信号。 0080 F.12是核磁共振纯相位编码测试系统的电子机柜， 用于放置并固定电子部件， 提 供各电子部件的良好接地。 0081 G.13是核磁共振纯相位编码测试系统的光纤核磁共振谱仪， 用于控制整个系统中 脉冲序列的所需射频发射、 梯度发射、 数据采集滤波、 增益放大、 正交解调、 K空间排布。 0082 H.14是核磁共振纯相位编码测试系统的射频功率放大器， 用于提供脉冲序列所需 的射频能量， 工。

38、作频率范围1MHz-350MHz， 最大功率350W。 0083 I.15、 16、 17对应为核磁共振纯相位编码测试系统X、 Y、 Z方向的梯度放大器， 用于 提供脉冲序列所需的梯度电流， 最大输出电流150A,最快爬升时间60us,具有过流过压自动 保护功能。 0084 J.18是核磁共振纯相位编码测试系统的电源部分， 为X、 Y、 Z梯度放大器提供150V 直流电压。 0085 K.19是核磁共振纯相位编码测试系统的水冷机， 连接到梯度线圈的水冷系统， 通 过制冷机带走循环水的温度。 0086 L.20是核磁共振纯相位编码测试系统的计算机主机部分， 主机安装有脉冲序列编 辑程序、 脉冲序。

39、列测试程序、 数据采集分析显示程序。 0087 M.21是核磁共振纯相位编码测试系统的计算机显示器， 用于显示计算机程序。 0088 N.22是核磁共振纯相位编码测试系统的前置放大器部分， 用于射频线圈中样品产 生的核磁共振信号， 前置放大器增益30dB,工作频率范围1MHz-60MHz， 噪声系数0.1dB。 0089 致密岩心核磁共振纯相位编码测试中， 快速突变的相位编码梯度场会在核磁共振 设备的极环、 极板、 磁钢等金属部件中产生涡流， 涡流会使实际的梯度磁场的前后沿变缓和 幅度减小， 同时它又迟后于梯度场的消失， 会造成了磁共振信号的相位编码线性度混淆、 变 形,出现错误的测试结果， 。

40、本发明使用的核磁共振设备的磁体具备防涡流层， 且核磁共振谱 仪共有12组数学波形预增强参数对梯度快速切换时进行涡流补偿， 确保在梯度快速编码期 间无涡流影响样品产生的核磁共振信。 说明书 6/9 页 10 CN 111650234 A 10 0090 高分辨率短TE的脉冲序列运行期间会产生大量的热量， 对梯度线圈的稳定性要求 极高， 为保证样品测试过程中温度不变， 且磁体不会因为梯度线圈发热产生频率漂移带来 的测量误差， 本发明摒弃常规梯度线圈， 进行重新规划设计。 在梯度线圈内布局铜质水管进 行水冷降温， 针对线圈不同部位的发热量不同进行布管， 对线圈磁场影响小， 散热效率高。 图2为核磁共。

41、振设备示意图， 梯度线圈-8为水冷梯度线圈， 被固定在磁体极环上， 并通过滤 波器与梯度功率放大器输出端相连， 同时水冷梯度线圈内铜质水管与外部水冷循环系统相 连接， 带走样品测量过程中的梯度电流产生的热量。 0091 现有常规核磁共振设备多使用射频线缆进行信号的传输， 致密岩心核磁共振信号 弱， 射频线缆及射频转接头会引起一定的核磁共振信号的损耗， 且容易引入外界干扰， 本发 明使用的核磁共振谱仪13为光纤传输模式， 致密岩心样品产生的核磁共振信号经过前置放 大器22放大后由电信号转换为光信号， 经光纤传输至核磁共振谱仪中。 0092 本方法通过填充一维K空间， 并对一维K空间进行傅里叶变换。

42、获得岩心不同位置流 体饱和度的分布， 本方法核磁共振脉冲序列 -SPI时序见图3， RF-23为射频硬脉冲， TR(Time of Repeatation)-26为重复时间， tp-27为梯度编码时间， 0为k空间中心编码梯度值， Gmax- 28为k空间边缘编码梯度值， RD(Relaxation Delay)-29为弛豫到稳态等待的时间。 每个TR 内完成一次被测样品的磁共振信号激发、 一步相位编码和一次磁共振信号采集。 0093 全部相位编码过程由k空间中心变化至正最大编码梯度， 等待RD后编码梯度再次 由k空间中心变化至负最大编码梯度， 完成由双半K空间组成组成1维K空间填充的过程； 。

43、编 码过程中为缩短激发时间我们采用硬脉冲设计， 同时为了充分提高横向磁化矢量的使用效 率以及消除纵向弛豫T1的加权影响， 本方法激励脉冲采用小角度硬脉冲， 多次小角度硬脉 冲激励和信号采集均是在阶梯相位编码的情况下进行， 最短回波时间TE由以下公式1给出： 0094 0095 其中， TRF为小角度硬脉冲激励的时间， tp为相位编码时间(Time of Phase Coding)。 使得在较短TE时间内(约100微秒)可以采集到岩心内流体产生的核磁共振信号。 0096 全部相位编码过程中磁共振信号衰减见公式2： 0097 0098 其中， M0为平衡态时的磁化矢量， tp为相位编码时间， T2。

44、*为有效横向弛豫时间， 为 翻转角度， M0与流体量成正比， 整个测试过程中呈现为质子密度加权， 磁共振信号不再受到 纵向时间T1和TR的影响， 同时纯相位编码技术在很大程度上不受磁化率引起的磁场变化、 化学位移和顺磁杂质引起的磁共振信号的影响， 为致密岩心核磁共振流体饱和度分布测试 较为理想的方法。 0099 在水、 油毛管压力曲线的测量中， 需要对油水核磁共振信号进行识别。 由于氧化氘 (D2O)具有不包含1H同位素， 无法产生核磁共振信号， 且物性与普通水相同等特点。 本方法使 用氧化氘(D2O)配置的模拟地层水作为被测试致密岩心的润湿相， 同时使用模拟油作为被 测试致密岩心的非润湿相，。

45、 致密岩心中油水两相共存时， 仅有油相产生核磁共信号， 通过油 水不同的含烃因子将离心后岩心中油相流体饱和度分布换算成水相流体饱和度分布， 进而 得到致密岩心内油水毛管压力曲线， 更适合油藏开发现场使用。 说明书 7/9 页 11 CN 111650234 A 11 0100 实施例选用了一块致密砂岩， 岩心孔隙度为8.3， 渗透率为0.91毫达西， 岩心为 直径26mm*45mm圆柱体形状。 0101 1.使用有机溶剂对岩心进行洗油处理， 使岩心荧光级别小于3级。 0102 2.将岩心在116条件下烘干至恒重， 并放入干燥器中冷却至室温， 将岩心样品置 于真空条件下保持。 0103 3.根据。

46、SY/T5336-2006规定测试岩心氦气孔隙度为8.3， 渗透率为0.91毫达西， 岩心直径为26mm， 长度为45mm。 0104 4.使用蒸馏水(H2O)配置模拟地层水1， 使用蒸馏水(D2O)配置模拟地层水2， 然后将 模拟地层水1和2按照SY/T5336规定的方法抽真空保持。 0105 5.配置模拟油， 然后按照SY/T5336规定的方法抽真空保持。 0106 6.开启核磁共振岩心分析设备， 将孔隙度为25、 长度为50mm的标准样放入核磁 共振射频线圈中， 加载FID脉冲序列， 使用默认参数搜索当前磁体氢核对应的共振频率为 15.00788665MHz， 改变射频脉冲的能量获得被测。

47、样品的90度翻转角以及180度翻转角， P90 为6us， P180为12us。 加载 -SPI序列， 设置Z方向梯度纯相位编码的测试范围为70mm， 平均次 数NS4， 接收增益控制RG30， 获得高SNR的一维双半K空间结果， 对一维双半K空间结果进 行傅里叶变换， 得到标准样的核磁共振信号分布。 根据结果不断调整岩心样品轴向位置， 直 至岩心大致处于设备视野中心， 标记岩芯标准样中心位置。 0107 7.按照SY/T5336规定的方法， 对步骤3的岩芯样品进行加压饱和模拟地层水1处 理， 最终含水饱和度为98.5， 称重。 0108 8.岩芯饱和完成后， 使用不产生磁共振信号的四氟乙烯密。

48、封带包裹岩芯样品， 以 减少样品内部液体的蒸发， 将包裹好的放入试管并垂直置入磁共振射频线圈内， 参考步骤6 标记的位置。 0109 9.对样品进行质子(1H)一维中心扫描的 -SPI核磁共振测量， 加载 -SPI序列， 设 置Z方向梯度纯相位编码的测试范围为60mm， 设置平均次数NS64， 接收增益控制RG100， 获得高SNR的一维双半K空间结果， 对一维双半K空间结果进行傅里叶变换， 并根据岩心样品 的总水量进行了归一化处理。 得到了H2O在样品长度上的分布， 见图6。 0110 10.重复步骤2、 步骤7， 用模拟地层水2代替模拟地层水1。 0111 11.步骤10饱和完成后取出岩芯。

49、， 去除四氟乙烯密封带， 在岩心样品的圆柱表面覆 盖一热缩聚四氟乙烯管， 以保持岩心内的纵向流体流动。 采用穿孔聚四氟乙烯衬套支撑岩 心， 置入充满模拟油的离心杯中， 如图4(1是密封盖， 2是岩心样品， 3是穿孔聚四氟乙烯衬 套， 4是多孔挡板)。 对4个相同的离心杯进行配重处理后置入离心机的工作转盘中， 如图5(1 为离心杯， 2为转轴)。 0112 12.设置离心机温度为5， 设置离心机转速为20000RPM， 离心时间2小时， 进行油 驱水状态的离心操作。 0113 13.将油驱水状态的岩心样品从离心杯中取出， 称重， 重复步骤8、 9， 并进行了归一 化处理， 得到了油相在样品长度上。

50、的分布， 见图7。 0114 14.根据油水含烃因子将致密岩心油驱水状态油相轴向核磁共振信号分布换算为 对应水核磁共振信号沿岩心长度的分布， 由初始饱和模拟地层水1磁共振信号与换算的水 信号作差， 得到油驱水后D2O的核磁共振信号分布。 说明书 8/9 页 12 CN 111650234 A 12 0115 15.根据离心后D2O核磁共振信号分布除以初始饱和状态的H2O核磁共振信号分布， 得到离心后岩心样品中水相饱和度(Sw)沿岩心长度的分布， 见图8。 011616.根据公式分别计算岩心上不同位置的油水毛细管 压力值。 0117 17.以不同位置含水饱和度为横轴， 对应水的毛管压力为纵轴， 。