一种数字岩心建模方法技术

本发明专利技术公开了一种数字岩心建模方法。步骤包括：S1、取样；S2、确认X射线分辨率起始值；S3、对圆柱形岩心样品三轴加压；S4、开启CT设备，通过X射线源对载物台上经过三轴加压的圆柱形岩心样品进行旋转扫描最终形成初始二维扫描灰度图像；S5、进行二值化分割获取二值灰度图像，进而形成三维数字岩心图像。本发明专利技术的有益效果是：通过核磁共振检测，确认X射线分辨率起始值，提高圆柱形岩心样品CT扫描使的便捷性，并且扫描时X射线分辨率大于该起始值，保证CT扫描层析成像孔隙精度，进而保证后续数字岩心建模精度。模精度。模精度。

【技术实现步骤摘要】
一种数字岩心建模方法


[0001]本专利技术涉及数字岩心
，特别涉及一种数字岩心建模方法。

技术介绍

[0002]现代化生产中，通过对地下岩层结构、成份等方面检测来获取石油、天然气等能源信息，对岩层内部结构常用检测方法为数字岩心技术，数字岩心技术在油藏工程以及地学领域得到广泛应用和推广，其中包括数字岩心建模及流动模拟技术等。目前国内外对于数字岩心建模方法多采用X射线扫描方法，即通过X射线扫描岩心样品，利用图像处理算法将二维图像构建为三维数字岩心图像。
[0003]三维数字岩心建模方法由物理法、数学法和CT扫描法，物理法成像结果准确，但是费用较高。数学方法虽然建立的数字岩心孔隙连通性较好，但是只适用于成岩过程简单岩石的数字岩心的重建，不适用于成岩过程复杂的岩石数字岩心重建。因此目前，CT扫描法是常用数字岩心建模方法，但是怎样确定CT扫描时X射线分辨率，以及怎样提高数字岩心构建精度，是本领域人员急需解决的问题。

技术实现思路

[0004]本专利技术的目的在于提供一种数字岩心建模方法，其通过核磁共振检测，确认X射线分辨率起始值，提高圆柱形岩心样品CT扫描使的便捷性，并且扫描时X射线分辨率大于该起始值，保证CT扫描层析成像孔隙精度，进而提高了数字岩心构建精度。
[0005]为实现上述目的，本专利技术采用如下技术方案，包括以下步骤：
[0006]S1、从岩心上截取圆柱形岩心样品；
[0007]S2、将所述圆柱形岩心样品放置在核磁共振检测设备上，通过核磁共振获取所述圆柱形岩心样品的T2图谱，计算T2图谱中T2曲线积分面积数值，通过高压压汞获取压汞孔喉半径分布与对应的T2图谱对比，根据岩心样品预设的定量关系将该计算结果转换得到孔喉半径分布频率；
[0008][0009]其中，T2为核磁横向弛豫时间，ms；ρ
t
为弛豫强度，Fs为孔隙形状因子，r
c
为孔隙半径；n为幂指数；
[0010]再根据该孔喉半径分布频率获取孔喉分布函数，对比饱和状态下T2谱换算的孔喉分布得到岩心样品的孔隙半径；
[0011][0012]其中，K为与孔喉比、形状因子相关的系数，r为喉道半径；
[0013]根据T2截止值范围确认X射线分辨率起始值为，当0<T2<10毫秒时，X射线分辨率为
r
±
0.01μm；当10<T2<100毫秒时，X射线分辨率为r
±
0.1μm；当T2＞100毫秒时，X射线分辨率为r
±
1μm；
[0014]S3、将所述圆柱形岩心样品由核磁共振检测设备取下，然后放在CT设备的载物台上，对圆柱形岩心样品三轴加压；
[0015]S4、开启CT设备，通过X射线源对载物台上经过三轴加压的圆柱形岩心样品进行旋转扫描，施加的X射线分辨率大于步骤S2确定的X射线分辨率起始值，通过该CT设备对圆柱形岩心样品进行层析成像，计算机记录灰度信息，最终形成初始二维扫描灰度图像；
[0016]S5、使用分水岭法对所述初始二维扫描灰度图像进行二值化分割，以分离图像中的孔隙空间和岩石骨架，获取二值灰度图像，确定孔隙和骨架边界，确定孔隙及网络骨架三维结构，形成三维数字岩心图像。
[0017]优选的是，步骤S1中所述圆柱形岩心样品的直径为4～26mm，高度为10～30mm。
[0018]优选的是，步骤S3中，所述三轴加压中，围压≥30Mpa，轴压≥30MPa。
[0019]优选的是，在所述CT设备的载物台上设有夹持器，在所述夹持器内设有加温系统，在所述圆柱形岩心样品进行层析成像时，开启所述加温系统对所述圆柱形岩心样品进行加温，使所述圆柱形岩心样品的温度与其取样处的地层温度一致。
[0020]优选的是，还包括：对步骤S4中获得的初始二维扫描灰度图像通过中值滤波法进行去躁处理。
[0021]优选的是，步骤S4中，施加的X射线分辨率＝分辨率起始值+1μm。
[0022]优选的是，步骤S5中，三维数字岩心图像采用“最大球”算法获取规则化的孔隙和孔喉模型来表征岩心样品的孔隙结构。
[0023]本专利技术的有益效果是：通过核磁共振检测，确认X射线分辨率起始值，提高圆柱形岩心样品CT扫描使的便捷性，并且扫描时X射线分辨率大于该起始值，保证CT扫描层析成像孔隙精度，进而提高了数字岩心构建精度。
附图说明
[0024]图1是专利技术一种数字岩心建模方法的流程图。
[0025]图2为本实施例所述岩心样品二维灰度图像。
[0026]图3为本实施例所述岩心样品分割孔隙后的二值化图像。
[0027]图4为本实施例所述岩心样品构建的三维模型。
具体实施方式
[0028]下面结合附图对专利技术做进一步的详细说明，以令本领域技术人员参照说明书文字能够据以实施。
[0029]应当理解，本文所使用的诸如“具有”、“包含”以及“包括”术语并不排除一个或其它元件或其组合的存在或添加。
[0030]如图1
‑
4所示，本专利技术的一种实现形式，为实现上述目的，本专利技术采用如下技术方案。
[0031]实施例1
[0032]S1、从岩心上截取圆柱形岩心样品；所述圆柱形岩心样品的直径为25mm，高度为
20mm。
[0033]S2、将所述圆柱形岩心样品放置在核磁共振检测设备上，通过核磁共振获取所述圆柱形岩心样品的T2图谱，计算T2图谱中T2曲线积分面积数值，通过高压压汞获取压汞孔喉半径分布与对应的T2图谱对比，根据岩心样品预设的定量关系将该计算结果转换得到孔喉半径分布频率；
[0034][0035]其中，T2为核磁横向弛豫时间，ms；ρ
t
为弛豫强度，Fs为孔隙形状因子，r
c
为孔隙半径；n为幂指数；
[0036]再根据该孔喉半径分布频率获取孔喉分布函数，对比饱和状态下T2谱换算的孔喉分布得到岩心样品的孔隙半径；
[0037][0038]其中，K为与孔喉比、形状因子相关的系数，r为喉道半径；
[0039]根据T2截止值范围确认X射线分辨率起始值为，当0<T2<10毫秒时，X射线分辨率为r
±
0.01μm；当10<T2<100毫秒时，X射线分辨率为r
±
0.1μm；当T2＞100毫秒时，X射线分辨率为r
±
1μm。
[0040]S3、将所述圆柱形岩心样品由核磁共振检测设备取下，然后放在CT设备的载物台上，对圆柱形岩心样品三轴加压；所述三轴加压中，围压＝50Mpa，轴压＝50MPa。
[0041]S4、开启CT设备，通过X射线源对载物台上经过三轴加压的圆柱形岩心样品进行旋转扫描，施加的X射线分辨率＝分辨率起始值
±
1μm。通过该CT设备对圆柱形岩本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种数字岩心建模方法，其特征在于，包括以下步骤：S1、从岩心上截取圆柱形岩心样品；S2、将所述圆柱形岩心样品放置在核磁共振检测设备上，通过核磁共振获取所述圆柱形岩心样品的T2图谱，计算T2图谱中T曲线积分面积数值，通过高压压汞获取压汞孔喉半径分布与对应的T2图谱对比，根据岩心样品预设的定量关系将该计算结果转换得到孔喉半径分布频率；其中，T2为核磁横向弛豫时间，ms；ρ
t
为弛豫强度，为孔隙形状因子，r
c
为孔隙半径；n为幂指数；再根据该孔喉半径分布频率获取孔喉分布函数，对比饱和状态下T2谱换算的孔喉分布得到岩心样品的孔隙半径；其中，K为与孔喉比、形状因子相关的系数，r为喉道半径；根据T2截止值范围确认X射线分辨率起始值，当0<T2<10毫秒时，X射线分辨率为r
±
0.01μm；当10<T2<100毫秒时，X射线分辨率为r
±
0.1μm；当T2＞100毫秒时，X射线分辨率为r
±
1μm；S3、将所述圆柱形岩心样品由核磁共振检测设备取下，然后放在CT设备的载物台上，对圆柱形岩心样品三轴加压；S4、开启CT设备，通过X射线源对载物台上经过三轴加压的圆柱形岩心样品进行旋转扫描，施加的X射线分辨率大于步骤S2确定的...

【专利技术属性】
技术研发人员：张铜耀，龚胜利，齐玉民，陈科，史长林，万欢，马良帅，何伟，张旭东，税蕾蕾，唐磊，杨发荣，何玉梅，黄向胜，
申请(专利权)人：中海油能源发展股份有限公司，
类型：发明
国别省市：