一种基于铁非弹性散射伽马的随钻气层识别装置及方法

1.本发明涉及测井技术领域，尤其是涉及一种基于铁非弹性散射伽马的随钻气层识别装置及方法。


背景技术：

2.随着石油勘探开发的不断深入，脉冲中子测井技术在随钻气层识别和划分方面发挥着越来越重要作用；基于天然气和油水等孔隙流体在中子减速和俘获能力方面的巨大差异，可以利用中子寿命测井、中子孔隙度测井以及中子俘获测井等方法进行气层识别，但由于热中子和俘获伽马射线在地层输运过程较为复杂，气层识别结果容易受地层泥质含量和地层水矿化度等因素影响。
3.高能快中子在地层输运过程简单，受地层环境因素影响较小，对气层识别更有利，但由于衰减快和探测效率低，很难满足测井需要；因此，通常采用与高能快中子分布相关的总非弹性散射伽马信息进行气层识别，但是容易受到地层密度衰减的影响。


技术实现要素：

4.有鉴于此，有必要提供一种基于铁非弹性散射伽马的随钻气层识别装置及方法，用以解决采用总非弹性散射伽马信息进行气层识别容易受到地层密度衰减的影响的技术问题。
5.为了实现上述目的，第一方面，本发明提供了一种基于铁非弹性散射伽马的随钻气层识别装置，包括开槽钻铤、可控中子源、源强探测器及伽马探测器，所述可控中子源、所述源强探测器及所述伽马探测器均固定于所述开槽钻铤上，所述源强探测器位于所述可控中子源与所述伽马探测器之间且用于获取所述可控中子源释放的高能中子通量，所述源强探测器被由钨镍铁材料制成的壳体所包裹，所述伽马探测器用于获取非弹性散射伽马能谱中的铁非弹性散射伽马信息，以进行气层识别。
6.进一步地，所述伽马探测器与所述可控中子源之间的距离为70cm。
7.第二方面，本发明还提供了一种基于铁非弹性散射伽马的随钻气层识别方法，适用于所述基于铁非弹性散射伽马的随钻气层识别装置，且包括：获取地层孔隙度，通过所述源强探测器获取源强探测器计数，并通过所述伽马探测器获取非弹性散射伽马能谱；
8.在所述非弹性散射伽马能谱上确定铁非弹性散射伽马峰的能量范围，并获取纯净铁非弹性散射伽马峰计数；
9.通过所述地层孔隙度、所述源强探测器计数及所述纯净铁非弹性散射伽马峰计数确定地层的含气饱和度。
10.进一步地，获取地层孔隙度，通过所述源强探测器获取源强探测器计数，并通过所述伽马探测器获取非弹性散射伽马能谱，具体包括：
11.将所述基于铁非弹性散射伽马的随钻气层识别装置下入到实际地层的钻孔的一预设位置；
12.通过所述可控中子源持续释放中子，通过所述源强探测器获取源强探测器计数，通过所述伽马探测器获取非弹性散射伽马能谱；
13.获取所述钻孔在所述预设位置处的孔隙度。
14.进一步地，在所述非弹性散射伽马能谱上确定铁非弹性散射伽马峰的能量范围，并获取纯净铁非弹性散射伽马峰计数，具体包括：
15.在所述非弹性散射伽马能谱上确定铁非弹性散射伽马峰的能量范围，得到所述铁非弹性散射伽马峰的起点和终点；
16.将所述铁非弹性散射伽马峰的起点和终点进行连线，确定一条分界线方程；
17.获取所述铁非弹性散射伽马峰位于所述分界线上方的面积，作为纯净铁非弹性散射伽马峰计数。
18.进一步地，在所述非弹性散射伽马能谱上确定铁非弹性散射伽马峰的能量范围的具体方法是：在所述非弹性散射伽马能谱上选取能量峰值为0.84mev的非弹性散射伽马峰，作为铁非弹性散射伽马峰，并确定其能量范围。
19.进一步地，通过所述地层孔隙度、所述源强探测器计数及所述纯净铁非弹性散射伽马峰计数确定地层的含气饱和度，具体包括：获取r
fe
与孔隙度及含气饱和度的刻度关系，其中，r
fe
为源强探测器计数与纯净铁非弹性散射伽马峰计数的比值；
20.根据所述源强探测器计数及所述纯净铁非弹性散射伽马峰计数确定地层的源强探测器计数与纯净铁非弹性散射伽马峰计数的比值r
fe
；
21.根据获取到的r
fe
与孔隙度及含气饱和度的刻度关系、地层孔隙度及地层的r
fe
，确定地层的含气饱和度。
22.进一步地，获取r
fe
与孔隙度及含气饱和度的刻度关系，其中，r
fe
为源强探测器计数与纯净铁非弹性散射伽马峰计数的比值，具体包括：将所述基于铁非弹性散射伽马的随钻气层识别装置下入到刻度井内，所述刻度井的不同位置处的孔隙度及含气饱和度均为已知；获取刻度井内不同位置处的孔隙度、含气饱和度及r
fe
，其中，r
fe
为源强探测器计数与纯净铁非弹性散射伽马峰计数的比值；建立孔隙度、含气饱和度及r
fe
的刻度关系。
23.进一步地，获取刻度井内不同位置处的孔隙度、含气饱和度及r
fe
的步骤中，获取刻度井内不同位置处的r
fe
的具体方法包括：通过所述可控中子源持续释放中子，通过所述源强探测器获取刻度井内不同位置处的源强探测器计数，并通过所述伽马探测器获取刻度井内不同位置处的非弹性散射伽马能谱；在所述非弹性散射伽马能谱上确定铁非弹性散射伽马峰的能量范围，并获取纯净铁非弹性散射伽马峰计数；根据所述源强探测器计数及所述纯净铁非弹性散射伽马峰计数确定地层的源强探测器计数与纯净铁非弹性散射伽马峰计数的比值r
fe
。
24.进一步地，建立孔隙度、含气饱和度及r
fe
的刻度关系，具体为：获取不同孔隙度及不同含气饱和度条件下的源强探测器计数与纯净铁非弹性散射伽马峰计数的比值r
fe
，建立表达r
fe
与孔隙度及含气饱和度sg的关系的图版。
25.与现有技术相比，本发明提出的技术方案的有益效果是：通过采用基于铁非弹性散射伽马信息进行气层识别的技术方案，即利用非弹性散射伽马能谱中的纯净铁非弹性散射伽马峰计数代替高能快中子信息进行气层识别，不仅保留了高能快中子受地层因素影响较小的优势，同时也克服了总非弹性散射伽马信息受地层密度影响的缺点，具有更高的气
层灵敏度。
附图说明
26.图1是本发明提供的基于铁非弹性散射伽马的随钻气层识别装置的一实施例的结构示意图；
27.图2是图1中剖面a
‑
a的剖视图；
28.图3是本发明提供的基于铁非弹性散射伽马的随钻气层识别方法的一实施例的流程示意图；
29.图4是伽马探测器获取的非弹性散射伽马能谱；
30.图5是图4中的非弹性散射伽马能谱中提取纯净铁非弹性散射伽马峰计数的方法的示意图；
31.图6是距离可控中子源70cm处的纯净铁非弹性散射伽马峰计数与高能快中子通量的关系；
32.图7是图3中步骤s1的流程示意图；
33.图8是图3中步骤s2的流程示意图；
34.图9是图3中步骤s3的流程示意图；
35.图10是理想条件下r
fe
与孔隙度及含气饱和度sg的关系的图版；
36.图中：1
‑
开槽钻铤、2
‑
可控中子源、3
‑
源强探测器、4
‑
伽马探测器、5
‑
壳体、6
‑
仪器外壳、11
‑
泥浆倒流通道、sg
‑
含气饱和度、r
fe
‑
源强探测器计数与纯净铁非弹性散射伽马峰计数的比值。
具体实施方式
37.下面结合附图来具体描述本发明的优选实施例，其中，附图构成本申请一部分，并与本发明的实施例一起用于阐释本发明的原理，并非用于限定本发明的范围。
38.实施例1：
39.请参照图1和图2，本发明提供了一种基于铁非弹性散射伽马的随钻气层识别装置，包括开槽钻铤1、可控中子源2、源强探测器3及伽马探测器4，开槽钻铤1内具有泥浆倒流通道11，所述可控中子源2、所述源强探测器3及所述伽马探测器4均固定于所述开槽钻铤1上，所述源强探测器3位于所述可控中子源2与所述伽马探测器4之间且用于获取所述可控中子源2释放的高能中子通量(即源强探测器计数)，所述源强探测器3被由钨镍铁材料制成的壳体5所包裹，源强探测器3采用钨镍铁屏蔽的原因是为了阻止地层中的中子进入源强探测器3，源强探测器3仅记录来自可控中子源2直射的中子信息，从而方便用来检测中子源波动，所述伽马探测器4用于获取非弹性散射伽马能谱中的铁非弹性散射伽马信息，以进行气层识别，可控中子源2、源强探测器3及伽马探测器4均收容于仪器外壳6内，仪器外壳6由钢材制成且内嵌于开槽钻铤1内。
40.具体地，请参照图1，所述可控中子源2为d
‑
t中子源，采用重复脉冲工作方式，所述可控中子源2的脉冲宽度介于10μs与30μs之间，所述可控中子源2的工作周期不小于100μs。
41.优选地，请参照图1，所述可控中子源2的脉冲宽度为10μs，所述可控中子源2的工作周期为200μs。
42.优选地，请参照图1，所述源强探测器为氦
‑
4探测器。
43.优选地，请参照图1，所述伽马探测器4与所述可控中子源2之间的距离为70cm。在这个距离下，伽马探测器4距离可控中子源2较远，非弹性散射伽马能谱中的纯净铁非弹性散射伽马峰计数与探测器附近的高能快中子通量近似成正比，几乎不受地层密度衰减影响，因此可以代替高能快中子信息进行气层识别。
44.优选地，请参照图1，所述伽马探测器4采用高能量分辨率和探测效率的溴化镧闪烁晶体。所述溴化镧闪烁晶体的长度为20cm。
45.实施例2：
46.请参照图3，本发明还提供了一种基于铁非弹性散射伽马的随钻气层识别方法，适用于所述基于铁非弹性散射伽马的随钻气层识别装置，且包括：
47.s1、获取地层孔隙度，通过所述源强探测器获取源强探测器计数，并通过所述伽马探测器获取非弹性散射伽马能谱(见图4)。
48.s2、在所述非弹性散射伽马能谱上确定铁非弹性散射伽马峰的能量范围(见图5)，并获取纯净铁非弹性散射伽马峰计数。
49.s3、通过所述地层孔隙度、所述源强探测器计数及所述纯净铁非弹性散射伽马峰计数确定地层的含气饱和度。
50.铁非弹性散射伽马射线是指由高能快中子与铁元素非弹性碰撞产生的伽马射线，地层中几乎没有铁，而开槽钻铤及仪器外壳几乎是由铁构成的，所以铁非弹性散射伽马射线的主要来源是开槽钻铤及仪器外壳铁非弹性散射伽马射线对非弹性散射伽马能谱的贡献主要表现在两个方面：(1)非弹性散射伽马能谱计数整体增加；(2)非弹性散射伽马能谱中出现铁非弹性散射伽马峰。
51.在伽马探测器距离中子源较远的情况下(本实例中伽马探测器与可控中子源之间的距离为70cm)，伽马探测器附近产生的铁非弹性散射伽马射线几乎不受地层密度衰减作用的影响，直接被伽马探测器记录，是能谱中铁非弹性散射伽马峰计数的主要来源。假定不考虑非弹性散射伽马能谱本底对铁非弹性散射伽马峰的贡献，那么纯净铁非弹性散射伽马峰计数可以表示如下
52.n
fe
≈nσ
in
‑
fe
n
f
53.其中，n
f
为伽马探测器附近的高能中子通量，在距离较远的情况下，高能中子通量难以通过仪器测量得到，n为一个快中子与fe原子发生非弹性碰撞释放的平均伽马光子数，σ
in
‑
fe
为一个快中子与开槽钻铤中的fe原子发生非弹性碰撞的概率。在开槽钻铤和伽马探测器参数不变的情况，n和σ
in
‑
fe
可以看做常数，纯净铁非弹性散射伽马峰计数与伽马探测器附近的高能快中子通量近似存在正比例关系(如图6所示)，因此，可以提取纯净铁非弹性散射伽马峰计数代替高能快中子通量进行气层识别，这样既可以保留高能快中子受地层因素影响较小的优势，又克服了在距离较远的情况下，高能中子通量难以通过仪器测量得到的弊端。
54.从图6还可以看出，纯净铁非弹性散射伽马峰计数还不受地层密度衰减影响，而现有技术中，采用总非弹性散射伽马信息代替高能快中子通量进行气层识别会受到地层密度衰减的影响，因此，采用非弹性散射伽马能谱中的纯净铁非弹性散射伽马峰计数代替高能快中子信息进行气层识别，克服了总非弹性散射伽马信息受地层密度影响的缺点，具有更
高的气层灵敏度。
55.需要指出的是：图6中的伽马探测器附近的高能快中子通量是通过模拟和理论分析得出的。
56.具体地，请参照图7，所述步骤s1具体包括：
57.s11、将所述基于铁非弹性散射伽马的随钻气层识别装置下入到实际地层的钻孔的一预设位置；
58.s12、通过所述可控中子源持续释放中子，通过所述源强探测器获取源强探测器计数，通过所述伽马探测器获取非弹性散射伽马能谱；
59.s13、获取所述钻孔在所述预设位置处的孔隙度，可以通过其他测井方法获取地层孔隙度，此为现有技术，对此不再赘述。
60.具体地，请参照图8，所述步骤s2具体包括：
61.s21、在所述非弹性散射伽马能谱上确定铁非弹性散射伽马峰的能量范围。
62.具体方法是：在所述非弹性散射伽马能谱上选取能量峰值为0.84mev的非弹性散射伽马峰，作为铁非弹性散射伽马峰。由理论分析可以得出：理想情况下，铁非弹性散射伽马射线的能量为0.84mev，从而只需要在非弹性散射伽马能谱上选取能量峰值为0.84mev的非弹性散射伽马峰，即为铁非弹性散射伽马峰。
63.本实施例中，请参照图4和图5，铁非弹性散射伽马峰的能量范围为0.78mev～1.00mev，选取能量0.78mev和1.00mev对应的伽马能谱数据点作为铁非弹性散射伽马峰的起点和终点。
64.s22、将所述铁非弹性散射伽马峰的起点和终点进行连线，得到分界线方程。所述分界线将非弹性散射伽马能谱中铁非弹性散射伽马峰分为上下两部分，分界线以上的峰面计数为纯净铁非弹性散射伽马峰计数，分界线以下的峰面计数为非弹性散射伽马能谱本底计数。
65.s23、获取所述铁非弹性散射伽马峰位于所述分界线上方的面积，得到纯净铁非弹性散射伽马峰计数，本实施例中，利用0.78mev～1.00mev范围内非弹性散射伽马能谱计数之和减去分界线下方的伽马能谱本底计数之和，得到纯净铁非弹性散射伽马峰计数。
66.具体地，请参照图9，所述步骤s3具体包括：
67.s31、获取r
fe
与孔隙度及含气饱和度的刻度关系，其中，r
fe
为源强探测器计数与纯净铁非弹性散射伽马峰计数的比值。
68.r
fe
与孔隙度及含气饱和度的刻度关系的获取方法包括；
69.s311、将所述基于铁非弹性散射伽马的随钻气层识别装置下入到刻度井内，所述刻度井的不同位置处的孔隙度及含气饱和度均为已知；
70.s312、获取刻度井内不同位置处的孔隙度、含气饱和度及r
fe
，其中，r
fe
为源强探测器计数与纯净铁非弹性散射伽马峰计数的比值，获取刻度井内不同位置处的r
fe
的具体方法包括：(1)通过所述可控中子源持续释放中子，通过所述源强探测器获取刻度井内不同位置处的源强探测器计数，并通过所述伽马探测器获取刻度井内不同位置处的非弹性散射伽马能谱；(2)在所述非弹性散射伽马能谱上确定铁非弹性散射伽马峰的能量范围，并获取纯净铁非弹性散射伽马峰计数；(3)根据所述源强探测器计数及所述纯净铁非弹性散射伽马峰计数确定地层的源强探测器计数与纯净铁非弹性散射伽马峰计数的比值r
fe
。
71.s313、建立孔隙度、含气饱和度及r
fe
的刻度关系。
72.具体为：获取不同孔隙度及不同含气饱和度条件下的源强探测器计数与纯净铁非弹性散射伽马峰计数的比值r
fe
，建立表达r
fe
与孔隙度及含气饱和度sg的关系的图版(见图10)。
73.s32、根据所述源强探测器计数及所述纯净铁非弹性散射伽马峰计数确定地层的源强探测器计数与纯净铁非弹性散射伽马峰计数的比值r
fe
。
74.为了避免中子源强对非弹性散射伽马测量的影响，采用源强探测器计数与纯净铁非弹性散射伽马峰计数的比值r
fe
进行气层划分，
[0075][0076]
其中，n
s
为源强探测器计数，n
fe
为纯净铁非弹性散射伽马峰计数。
[0077]
s33、根据获取到的r
fe
与孔隙度及含气饱和度的关系、地层孔隙度及地层的r
fe
，确定地层的含气饱和度。具体为：
[0078]
将获取到的源强探测器计数与纯净铁非弹性散射伽马峰计数的比值r
fe
及孔隙度数据投点到表达r
fe
与孔隙度及含气饱和度的关系的图版中，得到所述钻孔在所述预设位置处的含气饱和度。
[0079]
通过采用上述技术方案，相对于现有的采用总非弹性散射伽马信息进行气层识别的技术方案来说，本技术方案能够避免总非弹性散射伽马信息受地层密度衰减的影响，提高了对气层识别的灵敏度。如表1所示，纯净铁非弹性散射伽马峰计数在气水层的相对变化量要比总非弹性散射伽马计数信息高5％
‑
18％，这表明基于铁非弹性散射伽马信息的气层识别方法对气层识别更加灵敏。
[0080]
表1不同孔隙度条件下的气水层相对变化量
[0081]
[0082]
综上所述，通过采用基于铁非弹性散射伽马信息进行气层识别的技术方案，即利用非弹性散射伽马能谱中的纯净铁非弹性散射伽马峰计数代替高能快中子信息进行气层识别，不仅保留了高能快中子受地层因素影响较小的优势，同时也克服了总非弹性散射伽马信息受地层密度影响的缺点，具有更高的气层灵敏度。
[0083]
以上所述仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。