CN115012920B - 基于双clyc双粒子探测器的可控中子源多谱测井仪及测井方法 - Google Patents

基于双clyc双粒子探测器的可控中子源多谱测井仪及测井方法 Download PDF

Info

Publication number
CN115012920B
CN115012920B CN202210643863.5A CN202210643863A CN115012920B CN 115012920 B CN115012920 B CN 115012920B CN 202210643863 A CN202210643863 A CN 202210643863A CN 115012920 B CN115012920 B CN 115012920B
Authority
CN
China
Prior art keywords
double
clyc
particle detector
stratum
neutron
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
Application number
CN202210643863.5A
Other languages
English (en)
Other versions
CN115012920A (zh
Inventor
张锋
梁启轩
张慧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN202210643863.5A priority Critical patent/CN115012920B/zh
Publication of CN115012920A publication Critical patent/CN115012920A/zh
Application granted granted Critical
Publication of CN115012920B publication Critical patent/CN115012920B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Measurement Of Radiation (AREA)

Abstract

本发明公开了一种基于双CLYC双粒子探测器的可控中子源多谱测井仪,涉及油气勘探设备技术领域,基于双CLYC双粒子探测器的可控中子源多谱测井仪其仪器壳体内部由下而上依次设置有D‑T可控中子源、第一钨镍铁屏蔽体、近CLYC双粒子探测器、第二钨镍铁屏蔽体和远CLYC双粒子探测器,能够直接测量相同源距位置处地层热中子时间谱、伽马能谱和伽马时间谱,实现了对地层中子信息和伽马信息的精确探测。本发明还公开了一种基于双CLYC双粒子探测器可控中子源多谱测井仪的测井方法,实现了对地层中子孔隙度、中子伽马密度、含油饱和度和含气饱和度的同时测量以及对地层流体类型的准确划分,有利于地层的流体类型识别及含油气性评价。

Description

基于双CLYC双粒子探测器的可控中子源多谱测井仪及测井 方法
技术领域
本发明涉及油气勘探设备技术领域,尤其涉及一种基于双CLYC双粒子探测器的可控中子源多谱测井仪及测井方法。
背景技术
地层含油饱和度作为计算油气剩余储量、制定油气井开发计划的关键参数之一,对于油气井的勘探开发具有重要意义。目前,脉冲中子含油饱和度测井仪器设计中,常采用一个或多个伽马探测器测量地层伽马射线信息,利用伽马能谱信息或伽马时间谱信息评价地层的含油饱和度;或者采用多个中子探测器测量地层热中子信息,利用中子寿命差异评价地层的含油饱和度。在脉冲中子含油饱和度评价过程中,主要利用仪器测量的碳氧比或中子寿命信息,结合其他测井仪器测量的地层孔隙度信息对地层含油饱和度进行评价。由此可见,现阶段脉冲中子含油饱和度测井仪器的探测器仅能获取单一的中子信息或伽马信息,无法同时测量中子信息和伽马信息。
CLYC(Cs2LiYCl6:Ce3+)双粒子探测器具有同时测量中子信息和伽马信息的能力,其晶体密度为3.31g/cm3,伽马射线能量分辨率为4%,热中子能量分辨率为5.6%,伽马光产额为20000光子/MeV,中子光产额为70500光子/MeV,CLYC双粒子探测器能够在-30℃~180℃的温度区间内正常工作,具有优异的物理性能,能够适应井下严苛的工作环境。
因此,亟需将CLYC双粒子探测器应用于核地球物理勘探领域,为脉冲中子饱和度仪器的设计提供一种新思路。
发明内容
为解决上述技术问题,本发明公开了一种基于双CLYC(Cs2LiYCl6:Ce3+)双粒子探测器的可控中子源多谱测井仪及测井方法,本发明通过对双CLYC双粒子探测器可控中子源多谱测井仪的仪器结构、脉冲测量时序和仪器测量信息进行设计,利用双CLYC双粒子探测器同时获取地层的中子和伽马信息,仪器一次下井即可同时测量地层的中子孔隙度、中子伽马密度、含油饱和度和含气饱和度并准确确定地层的流体类型,有利于地层的流体类型识别及含油气性评价。
为实现上述目的,本发明采用下述技术方案:
一种基于双CLYC双粒子探测器的可控中子源多谱测井仪,所述可控中子源多谱测井仪的仪器壳体内部由下而上依次设置有D-T可控中子源、第一钨镍铁屏蔽体、近CLYC双粒子探测器、第二钨镍铁屏蔽体和远CLYC双粒子探测器,近CLYC双粒子探测器和远CLYC双粒子探测器,能够同时测量地层的热中子时间谱、伽马能谱和伽马时间谱,其中,伽马能谱包括碳氧比模式总伽马能谱和碳氧比模式俘获伽马能谱,伽马时间谱包括中子寿命模式伽马时间谱和本底伽马时间谱。
优选地,所述D-T可控中子源的产额不低于1×108/s;所述仪器壳体的直径设置为43~90mm;所述近CLYC双粒子探测器的直径设置为20~60mm、长度设置为40~60mm,与D-T可控中子源之间的距离设置为30~35cm;所述远CLYC双粒子探测器的直径设置为20~60mm、长度设置为100~140mm,与D-T可控中子源之间的距离设置为50~55cm;所述近CLYC双粒子探测器和远CLYC双粒子探测器采集的伽马能谱能量范围为0.1~10.3MeV,分为256道;
所述近CLYC双粒子探测器和远CLYC双粒子探测器均为Cs2LiYCl6:Ce3+闪烁晶体探测器,Cs2LiYCl6:Ce3+晶体的密度为3.31g/cm3,伽马能量分辨率为4~5%,热中子能量分辨率为4~6%。
优选地,所述第一钨镍铁屏蔽体的长度设置为7~12cm,第二钨镍铁屏蔽体的长度设置为1~2cm,第一钨镍铁屏蔽体和第二钨镍铁屏蔽体的直径均设置为40~90mm,第一钨镍铁屏蔽体和第二钨镍铁屏蔽体中均含有1%的硼元素。
优选地,所述可控中子源多谱测井仪脉冲测量时序的周期时长为100ms,包括时长为99ms的时间段T1和时长为1ms的时间段T2;时间段T1内包括4个记录周期,各记录周期内包括300个碳氧比模式记录周期和1个中子寿命模式记录周期,碳氧比模式记录周期的总时长为75μs,包括时间门G1和时间门G2,时间门G1的记录范围为第10~30μs,用于控制D-T可控中子源以脉冲形式发射快中子以及记录碳氧比模式总伽马能谱,时间门G2的记录范围为第35~75μs,用于记录碳氧比模式俘获伽马能谱和碳氧比模式伽马时间谱,中子寿命记录周期的总时长为2250μs,其中设置有时间门G3,时间门G3的记录范围为第75~2325μs,用于记录中子寿命模式伽马时间谱和热中子时间谱,时间段T2内设置有时间门G4,时间门G4的记录范围为2325~3325μs,用于记录本底伽马时间谱。
一种基于双CLYC双粒子探测器可控中子源多谱测井仪的测井方法,采用如上所述的基于双CLYC双粒子探测器的可控中子源多谱测井仪,利用基于双CLYC双粒子探测器的可控中子源多谱测井仪在井中贴井壁连续测量,其特征在于,同时测量地层的中子孔隙度、中子伽马密度、含油饱和度和含气饱和度,并对地层内的流体类型进行划分,综合评价地层的含油气性,具体包括以下步骤:
步骤1,将双CLYC双粒子探测器的可控中子源多谱测井仪在井中贴井壁连续测量,利用D-T可控中子源以脉冲形式发射快中子,近CLYC双粒子探测器与远CLYC双粒子探测器分别记录不同位置处的热中子时间谱、伽马能谱和伽马时间谱;
步骤2,根据近CLYC双粒子探测器与远CLYC双粒子探测器测量的热中子时间谱,获取近CLYC双粒子探测器与远CLYC双粒子探测器的热中子计数比,计算得到地层的中子孔隙度,如式(1)所示:
式中,φ为地层的中子孔隙度,Rth为近CLYC双粒子探测器与远CLYC双粒子探测器的热中子计数比,b0、b1、b2和b3均为中子孔隙度刻度系数;
步骤3,根据近CLYC双粒子探测器和远CLYC双粒子探测器测量的碳氧比模式总伽马能谱、碳氧比模式俘获伽马能谱,通过在碳氧比模式总伽马能谱中扣除碳氧比模式俘获伽马能谱得到非弹伽马能谱,分别得到近CLYC双粒子探测器与远CLYC双粒子探测器的非弹伽马能谱,确定近CLYC双粒子探测器和远CLYC双粒子探测器的非弹伽马计数比,计算得到地层的中子伽马密度,如式(2)所示:
ρ=a1ln(Rin)+a2ln(Rth)+a3 (2)
式中,ρ为地层的中子伽马密度,Rin为近CLYC双粒子探测器与远CLYC双粒子探测器的非弹伽马计数比,a1、a2、a3均为中子伽马密度刻度系数;
步骤4,根据近CLYC双粒子探测器的非弹伽马能谱获取近CLYC双粒子探测器的碳氧比值,再根据近CLYC双粒子探测器的碳氧比模式俘获伽马能谱获取近CLYC双粒子探测器的硅钙比值,计算得到地层含油饱和度,如式(3)所示:
式中,So为地层含油饱和度,COR为近CLYC双粒子探测器的碳氧比值,SCR为近CLYC双粒子探测器的硅钙比值,kw、bw、c1、c2和c3均为含油饱和度刻度系数;
步骤5,根据近CLYC双粒子探测器的非弹伽马能谱和热中子时间谱,确定近CLYC双粒子探测器的非弹热中子计数比值,再根据远CLYC双粒子探测器的非弹伽马能谱和热中子时间谱,确定远CLYC双粒子探测器的非弹热中子计数比值,计算得到地层含气饱和度,如式(4)所示:
式中,Sg为地层含气饱和度,RIT为近CLYC双粒子探测器非弹热中子计数比值与远CLYC双粒子探测器非弹热中子计数比值之间的差值,RITw为地层饱含水时近CLYC双粒子探测器非弹热中子计数比值与远CLYC双粒子探测器非弹热中子计数比值之间的差值,RITg为地层饱含气时近CLYC双粒子探测器非弹热中子计数比值与远CLYC双粒子探测器非弹热中子计数比值之间的差值,d1、d2、d3均为地层饱含水时双CLYC双粒子探测器的含气饱和度刻度系数,f1、f2、f3均为地层饱含气时双CLYC双粒子探测器的含气饱和度刻度系数;
步骤6,根据近CLYC双粒子探测器的伽马时间谱获取地层的宏观俘获截面值,结合近CLYC双粒子探测器的碳氧比值以及近、远CLYC双粒子探测器的热中子计数比,计算地层流体类型划分参数,如式(5)所示:
式中,D为地层流体类型划分参数,COR为近CLYC双粒子探测器的碳氧比值,Rth为近CLYC双粒子探测器与远CLYC双粒子探测器的热中子计数比,Σ为地层的宏观俘获截面值,Σma为地层骨架的宏观俘获截面值,A、B、C均为仪器刻度系数;
步骤7,根据地层的地层流体类型划分参数和中子孔隙度,基于MCNP数值模拟方法建立地层流体类型划分图版,确定地层内流体的类型,综合地层的含油饱和度、含气饱和度和流体类型,评价地层的含油气性。
优选地,所述步骤7中,具体包括以下步骤:
步骤7.1,根据双CLYC双粒子探测器可控中子源多谱测井仪的结构参数,结合双CLYC双粒子探测器可控中子源多谱测井仪在井中的实际探测情况,构建双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型;
步骤7.2,依次将双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型中地层设置为饱含水、饱含油和饱含气的砂岩地层,分别改变双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型中地层的孔隙度,模拟得到不同地层条件下的双CLYC双粒子探测器的热中子时间谱、伽马能谱和伽马时间谱,分别计算各地层条件所对应的地层流体类型划分参数;
步骤7.3,依次根据饱含砂岩地层、饱含油砂岩地层和饱含气砂岩地层中各地层孔隙度对应的地层流体类型划分参数,分别确定砂岩地层饱含水、饱含油和饱含气时双CLYC双粒子探测器可控中子源多谱测井仪地层流体类型划分参数的测量响应,绘制地层流体类型划分图版;
步骤7.4,根据地层流体类型划分参数和地层中子孔隙度,利用地层流体类型划分图版确定地层内流体的类型,综合地层的含油饱和度、含气饱和度和流体类型,评价地层的含油气性。
优选地,所述步骤7.2中,依次将双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型中的地层孔隙度设置为0、5%、10%、15%、20%、25%、30%。
本发明所带来的有益技术效果:
本发明基于CLYC(Cs2LiYCl6:Ce3+)双粒子探测器的探测特性,将CLYC双粒子探测器应用于脉冲中子饱和度测井仪器设计中,通过对仪器结构、脉冲测量时序和仪器测量信息进行研究,提出了一种基于双CLYC双粒子探测器的可控中子源多谱测井仪,该可控中子源多谱测井仪实现了对地层中子和伽马信息的同步测量,仪器一次下井即可获取地层的热中子时间谱、碳氧比模式总伽马能谱、碳氧比模式俘获伽马能谱、中子寿命模式伽马时间谱和地层本底伽马时间谱等多谱信息,即实现了对地层中子孔隙度、中子伽马密度、含油饱和度、含气饱和度的同时测量以及地层流体类型的准确获取,减少了仪器下井次数的同时有效提升了地层含油气性评价的准确性,为脉冲中子饱和度测井仪器的设计提供了一种新思路。
本发明还提出了一种基于双CLYC双粒子探测器可控中子源多谱测井仪的测井方法,利用可控中子源多谱测井仪的CLYC双粒子探测器记录相同源距位置处地层的中子场和伽马场,实现了对地层中子信息和伽马信息的同时测量,相比于常规利用不同源距中子探测器或伽马探测器组合探测的测井方法以及通过在探测器外包裹硼套的次生伽马中子记录方法,本发明方法中CLYC双粒子探测器记录的中子场和伽马场均具有明确的对应性和良好的统计性,实现了对相同源距位置处地层热中子时间谱、伽马能谱和伽马时间谱的直接测量,保证了地层流体类型划分的准确性,为地层的含油气性评价提供了新方法。
附图说明
图1为本发明基于双CLYC双粒子探测器的可控中子源多谱测井仪结构示意图。
图2为基于双CLYC双粒子探测器的可控中子源多谱测井仪的脉冲测量时序示意图。
图3为双CLYC双粒子探测器热中子计数比与地层中子孔隙度的刻度关系图。
图4为不同地层条件下中子伽马密度计算值与真实值的对比图。
图5为含油地层中近CLYC双粒子探测器碳氧比值的测量响应。
图6为含气地层中双CLYC双粒子探测器非弹热中子计数比差值的测量响应。
图7为地层流体类型划分图版。
图中,1、D-T可控中子源,2、第一钨镍铁屏蔽体,3、近CLYC双粒子探测器,4、第二钨镍铁屏蔽体,5、远CLYC双粒子探测器,6、地层,7、井壁,8、仪器壳体。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
一种基于双CLYC双粒子探测器的可控中子源多谱测井仪,如图1所示,可控中子源多谱测井仪的仪器壳体内部由下而上依次设置有D-T可控中子源1、第一钨镍铁屏蔽体2、近CLYC双粒子探测器3、第二钨镍铁屏蔽体4和远CLYC双粒子探测器5,近CLYC双粒子探测器3和远CLYC双粒子探测器5,能够同时测量地层的热中子时间谱、伽马能谱(包括碳氧比模式总伽马能谱和碳氧比模式俘获伽马能谱)和伽马时间谱(包括中子寿命模式伽马时间谱和本底伽马时间谱)。
本实施例中基于双CLYC双粒子探测器的可控中子源多谱测井仪,其仪器参数如下所示:
基于双CLYC双粒子探测器的可控中子源多谱测井仪,仪器壳体8的直径设置为43mm,D-T可控中子源1的产额不低于1×108/s。近CLYC双粒子探测器3和远CLYC双粒子探测器5共同构成双CLYC双粒子探测器,其中,近CLYC双粒子探测器3的直径设置为25.4mm、长度设置为50.8mm,与D-T可控中子源1之间的距离为30cm,远CLYC双粒子探测器5的直径设置为25.4mm、长度设置为101.6mm,与D-T可控中子源1之间的距离为55cm;近CLYC双粒子探测器3和远CLYC双粒子探测器5均为Cs2LiYCl6:Ce3+闪烁晶体探测器,近CLYC双粒子探测器3和远CLYC双粒子探测器5采集的伽马能谱能量范围为0.1~10.3MeV且分为256道。Cs2LiYCl6:Ce3 +晶体的密度为3.31g/cm3,伽马能量分辨率为4~5%,热中子能量分辨率为4~6%,伽马光产额20000光子/MeV,中子光产额70500光子/MeV,在-30℃~180℃的温度区间内正常工作,能够满足井下严苛探测环境的需求。
第一钨镍铁屏蔽体2设置于D-T可控中子源1与近CLYC双粒子探测器3之间,第二钨镍铁屏蔽体4设置于近CLYC双粒子探测器3与远CLYC双粒子探测器5之间,第一钨镍铁屏蔽体2和第二钨镍铁屏蔽体4均采用掺杂有1%硼元素的钨镍铁材料制成,第一钨镍铁屏蔽体2的直径设置为42mm、长度设置为800mm,第二钨镍铁屏蔽体4的直径设置为42mm、长度设置为20mm。
本实施例中基于双CLYC双粒子探测器的可控中子源多谱测井仪,其脉冲测量时序设计如图2所示,脉冲测量时序的周期时长为100ms,包括时长为99ms的时间段T1和时长为1ms的时间段T2;时间段T1内包括4个记录周期,各记录周期内包括300个碳氧比模式记录周期和1个中子寿命模式记录周期,碳氧比模式记录周期的总时长为75μs,包括时间门G1和时间门G2,时间门G1的记录范围为第10~30μs,用于控制D-T可控中子源以脉冲形式发射快中子以及记录碳氧比模式总伽马能谱,时间门G1记录的第10~15μs,D-T可控中子源处于预热状态,时间门G1记录的第15~30μs,D-T可控中子源以脉冲形式发生快中子;时间门G2的记录范围为第35~75μs,用于记录碳氧比模式俘获伽马能谱和碳氧比模式伽马时间谱;中子寿命记录周期的总时长为2250μs,其中设置有时间门G3,时间门G3的记录范围为第75~2325μs,用于记录中子寿命模式伽马时间谱和热中子时间谱,时间段T2内设置有时间门G4,时间门G4的记录范围为2325~3325μs,用于记录本底伽马时间谱。
本发明还提出了一种基于双CLYC双粒子探测器可控中子源多谱测井仪的测井方法,采用如上所述的基于双CLYC双粒子探测器的可控中子源多谱测井仪,能够同时测量地层的中子孔隙度、中子伽马密度、含油饱和度、含气饱和度以及确定地层的流体类型,综合评价地层的含油气性,具体包括以下步骤:
步骤1,将双CLYC双粒子探测器的可控中子源多谱测井仪在井中贴井壁7连续测量,利用D-T可控中子源以脉冲形式发射快中子,近CLYC双粒子探测器与远CLYC双粒子探测器分别记录不同位置处的热中子时间谱、伽马能谱和伽马时间谱,具体包括热中子时间谱、碳氧比模式总伽马能谱、碳氧比模式俘获伽马能谱、中子寿命模式伽马时间谱和本底伽马时间谱。
步骤2,根据近CLYC双粒子探测器与远CLYC双粒子探测器测量的热中子时间谱,获取近CLYC双粒子探测器与远CLYC双粒子探测器的热中子计数比,计算得到地层的中子孔隙度,如式(1)所示:
式中,φ为地层的中子孔隙度,Rth为近CLYC双粒子探测器与远CLYC双粒子探测器的热中子计数比,b0、b1、b2和b3均为中子孔隙度刻度系数,可通过标准刻度井刻度或MCNP数值模拟获得。
本实施例中中子孔隙度刻度系数通过MCNP数值模拟获得,根据双CLYC双粒子探测器可控中子源多谱测井仪的仪器参数和实际测量情况,利用MCNP数值模拟软件建立双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型,将建立双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型的地层6设置为饱含水的石灰岩地层,依次将建立双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型的地层孔隙度设置为0、5%、10%、15%、20%、25%和30%,分别模拟得到不同地层孔隙度条件下近CLYC双粒子探测器与远CLYC双粒子探测器测量的热中子时间谱,计算各地层条件下近CLYC双粒子探测器与远CLYC双粒子探测器的热中子计数比,得到双CLYC双粒子探测器热中子计数比与地层中子孔隙度的刻度关系,如图3所示,通过对双CLYC双粒子探测器热中子计数比与地层中子孔隙度的刻度关系进行拟合,确定中子孔隙度刻度系数b0、b1、b2和b3
步骤3,根据近CLYC双粒子探测器与远CLYC双粒子探测器测量的碳氧比模式总伽马能谱和碳氧比模式俘获伽马能谱,通过在碳氧比模式总伽马能谱中扣除碳氧比模式俘获伽马能谱得到非弹伽马能谱,分别得到近CLYC双粒子探测器与远CLYC双粒子探测器的非弹伽马能谱,确定近CLYC双粒子探测器和远CLYC双粒子探测器的非弹伽马计数比,计算得到地层的中子伽马密度,如式(2)所示:
ρ=a1ln(Rin)+a2ln(Rth)+a3 (2)
式中,ρ为地层的中子伽马密度,Rin为近CLYC双粒子探测器与远CLYC双粒子探测器的非弹伽马计数比,a1、a2、a3均为中子伽马密度刻度系数,通过将仪器置于标准刻度井中刻度后确定。
采用刻度后的双CLYC双粒子探测器可控中子源多谱测井仪对实测井进行测量,根据近CLYC双粒子探测器和远CLYC双粒子探测器的非弹伽马计数比计算得到地层中子伽马密度的计算值,并与地层的真实密度进行对比,如图4所示,发现采用本发明方法计算的地层中子伽马密度与地层真实密度一致,验证了本发明双CLYC双粒子探测器可控中子源多谱测井仪对地层中子伽马密度测量的准确性。
步骤4,根据近CLYC双粒子探测器的非弹伽马能谱获取近CLYC双粒子探测器的碳氧比值,再根据近CLYC双粒子探测器的碳氧比模式俘获伽马能谱获取近CLYC双粒子探测器的硅钙比值,计算得到地层含油饱和度,如式(3)所示:
式中,So为地层含油饱和度,COR为近CLYC双粒子探测器的碳氧比值,SCR为近CLYC双粒子探测器的硅钙比值,kw、bw、c1、c2和c3均为含油饱和度刻度系数。
本实施例中将双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型中的地层依次设置为饱含水砂岩地层和饱含水石灰岩地层,分别将地层孔隙度设置为0、5%、10%、15%、20%、25%和30%,模拟得到不同地层条件下近CLYC双粒子探测器测量的碳氧比模式总伽马能谱和碳氧比模式俘获伽马能谱,通过在碳氧比模式总伽马能谱中扣除碳氧比模式俘获伽马能谱得到非弹伽马能谱,得到不同地层条件下近CLYC双粒子探测器的非弹伽马能谱。
根据不同地层条件下近CLYC双粒子探测器的非弹伽马能谱获取近CLYC双粒子探测器的碳氧比值,再根据不同地层条件下近CLYC双粒子探测器的碳氧比模式俘获伽马能谱获取近CLYC双粒子探测器的硅钙比值,以近CLYC双粒子探测器的硅钙比值作为横坐标、碳氧比值作为纵坐标绘制碳氧硅钙交会图,通过分析碳氧硅钙交会图确定含油饱和度刻度系数kw和bw。
再利用双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型进行模拟,通过将双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型中的地层依次设置为砂岩地层,在地层饱含水和饱含油的条件下分别将地层孔隙度设置为0、5%、10%、15%、20%、25%和30%,模拟得到不同地层条件下近CLYC双粒子探测器的碳氧比值,如图5所示,计算相同孔隙度条件下饱含油砂岩地层碳氧比值与饱含水砂岩地层碳氧比值的差值,通过对不同孔隙度条件下饱含油砂岩地层碳氧比值与饱含水砂岩地层碳氧比值的差值曲线进行拟合,确定含油饱和度刻度系数c1、c2和c3
步骤5,根据近CLYC双粒子探测器的非弹伽马能谱和热中子时间谱,确定近CLYC双粒子探测器的非弹热中子计数比值,再根据远CLYC双粒子探测器的非弹伽马能谱和热中子时间谱,确定远CLYC双粒子探测器的非弹热中子计数比值,计算得到地层含气饱和度,如式(4)所示:
式中,Sg为地层含气饱和度,RIT为近CLYC双粒子探测器非弹热中子计数比值与远CLYC双粒子探测器非弹热中子计数比值之间的差值,RITw为地层饱含水时近CLYC双粒子探测器非弹热中子计数比值与远CLYC双粒子探测器非弹热中子计数比值之间的差值,RITg为地层饱含气时近CLYC双粒子探测器非弹热中子计数比值与远CLYC双粒子探测器非弹热中子计数比值之间的差值,d1、d2、d3均为地层饱含水时双CLYC双粒子探测器的含气饱和度刻度系数,f1、f2、f3均为地层饱含气时双CLYC双粒子探测器的含气饱和度刻度系数。
本实施例中利用MCNP数值模拟确定地层饱含水时双CLYC双粒子探测器的含气饱和度刻度系数d1、d2、d3以及地层饱含气时双CLYC双粒子探测器的含气饱和度刻度系数f1、f2、f3。将双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型中的地层设置为饱含水砂岩地层和饱含气砂岩地层,将地层孔隙度设置为0、5%、10%、15%、20%、25%和30%,模拟得到不同地层条件下双CLYC双粒子探测器的非弹伽马能谱和热中子时间谱,分别计算得到不同孔隙度条件下饱含水砂岩地层中近CLYC双粒子探测器非弹热中子计数比值与远CLYC双粒子探测器非弹热中子计数比值之间的差值以及不同孔隙度条件下饱含气砂岩地层中近CLYC双粒子探测器非弹热中子计数比值与远CLYC双粒子探测器非弹热中子计数比值之间的差值,如图6所示。根据饱含水砂岩地层中双CLYC双粒子探测器非弹热中子计数比差值随孔隙度的变化规律,拟合分析后得到地层饱含水时双CLYC双粒子探测器的含气饱和度刻度系数d1、d2和d3。再根据饱含气砂岩地层中双CLYC双粒子探测器非弹热中子计数比差值随孔隙度的变化规律,拟合分析后得到地层饱含气时双CLYC双粒子探测器的含气饱和度刻度系数f1、f2和f3
步骤6,根据近CLYC双粒子探测器的伽马时间谱获取地层的宏观俘获截面值,结合近CLYC双粒子探测器的碳氧比值以及近、远CLYC双粒子探测器的热中子计数比,计算地层流体类型划分参数,如式(5)所示:
式中,D为地层流体类型划分参数,COR为近CLYC双粒子探测器的碳氧比值,Rth为近CLYC双粒子探测器与远CLYC双粒子探测器的热中子计数比,Σ为地层的宏观俘获截面值,Σma为地层骨架的宏观俘获截面值,A、B、C均为仪器刻度系数。
步骤7,根据计算得到的地层流体类型划分参数,基于MCNP数值模拟方法建立地层流体类型划分图版,确定地层的流体类型,综合地层的含油饱和度、含气饱和度和流体类型,评价地层的含油气性,具体包括以下步骤:
步骤7.1,根据双CLYC双粒子探测器可控中子源多谱测井仪的结构参数,构建双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型,双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型中双CLYC双粒子探测器可控中子源多谱测井仪在井中贴井壁测量。
步骤7.2,依次将双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型中地层设置为饱含水、饱含油和饱含气的砂岩地层,分别将双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型中砂岩地层的孔隙度设置为0、5%、10%、15%、20%、25%、30%,模拟得到不同地层条件下的双CLYC双粒子探测器的热中子时间谱、伽马能谱和伽马时间谱,确定不同地层条件下近CLYC双粒子探测器的碳氧比值、近CLYC双粒子探测器与远CLYC双粒子探测器的热中子计数比、地层的宏观俘获截面值,分别计算各地层条件所对应的地层流体类型划分参数。
步骤7.3,依次根据不同孔隙度条件下饱含水砂岩地层、饱含油砂岩地层和饱含气砂岩地层所对应的地层流体类型划分参数,分别确定砂岩地层饱含水、饱含油和饱含气时双CLYC双粒子探测器可控中子源多谱测井仪地层流体类型划分参数的测量响应,绘制地层流体类型划分图版,如图7所示。
步骤7.4,根据地层流体类型划分参数和地层中子孔隙度,利用地层流体类型划分图版确定地层内流体的类型,综合地层的含油饱和度、含气饱和度和流体类型,评价地层的含油气性。
本发明方法能够准确获取地层的中子孔隙度、中子伽马密度、含油饱和度和含气饱和度,实现了对地层流体类型的准确划分,有利于地层的含油气性评价。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。

Claims (6)

1.一种基于双CLYC双粒子探测器可控中子源多谱测井仪的测井方法,采用基于双CLYC双粒子探测器的可控中子源多谱测井仪,利用基于双CLYC双粒子探测器的可控中子源多谱测井仪在井中贴井壁连续测量,所述可控中子源多谱测井仪的仪器壳体内部由下而上依次设置有D-T可控中子源、第一钨镍铁屏蔽体、近CLYC双粒子探测器、第二钨镍铁屏蔽体和远CLYC双粒子探测器,近CLYC双粒子探测器和远CLYC双粒子探测器,能够同时测量地层的热中子时间谱、伽马能谱和伽马时间谱,其中,伽马能谱包括碳氧比模式总伽马能谱和碳氧比模式俘获伽马能谱,伽马时间谱包括中子寿命模式伽马时间谱和本底伽马时间谱,其特征在于,同时测量地层的中子孔隙度、中子伽马密度、含油饱和度和含气饱和度,并对地层内的流体类型进行划分,综合评价地层的含油气性,具体包括以下步骤:
步骤1,将双CLYC双粒子探测器的可控中子源多谱测井仪在井中贴井壁连续测量,利用D-T可控中子源以脉冲形式发射快中子,近CLYC双粒子探测器与远CLYC双粒子探测器分别记录不同位置处的热中子时间谱、伽马能谱和伽马时间谱;
步骤2,根据近CLYC双粒子探测器与远CLYC双粒子探测器测量的热中子时间谱,获取近CLYC双粒子探测器与远CLYC双粒子探测器的热中子计数比,计算得到地层的中子孔隙度,如式(1)所示:
式中,φ为地层的中子孔隙度,Rth为近CLYC双粒子探测器与远CLYC双粒子探测器的热中子计数比,b0、b1、b2和b3均为中子孔隙度刻度系数;
步骤3,根据近CLYC双粒子探测器和远CLYC双粒子探测器测量的碳氧比模式总伽马能谱、碳氧比模式俘获伽马能谱,通过在碳氧比模式总伽马能谱中扣除碳氧比模式俘获伽马能谱得到非弹伽马能谱,分别得到近CLYC双粒子探测器与远CLYC双粒子探测器的非弹伽马能谱,确定近CLYC双粒子探测器和远CLYC双粒子探测器的非弹伽马计数比,计算得到地层的中子伽马密度,如式(2)所示:
ρ=a1ln(Rin)+a2ln(Rth)+a3 (2)
式中,ρ为地层的中子伽马密度,Rin为近CLYC双粒子探测器与远CLYC双粒子探测器的非弹伽马计数比,a1、a2、a3均为中子伽马密度刻度系数;
步骤4,根据近CLYC双粒子探测器的非弹伽马能谱获取近CLYC双粒子探测器的碳氧比值,再根据近CLYC双粒子探测器的碳氧比模式俘获伽马能谱获取近CLYC双粒子探测器的硅钙比值,计算得到地层含油饱和度,如式(3)所示:
式中,So为地层含油饱和度,COR为近CLYC双粒子探测器的碳氧比值,SCR为近CLYC双粒子探测器的硅钙比值,kw、bw、c1、c2和c3均为含油饱和度刻度系数;
步骤5,根据近CLYC双粒子探测器的非弹伽马能谱和热中子时间谱,确定近CLYC双粒子探测器的非弹热中子计数比值,再根据远CLYC双粒子探测器的非弹伽马能谱和热中子时间谱,确定远CLYC双粒子探测器的非弹热中子计数比值,计算得到地层含气饱和度,如式(4)所示:
式中,Sg为地层含气饱和度,RIT为近CLYC双粒子探测器非弹热中子计数比值与远CLYC双粒子探测器非弹热中子计数比值之间的差值,RITw为地层饱含水时近CLYC双粒子探测器非弹热中子计数比值与远CLYC双粒子探测器非弹热中子计数比值之间的差值,RITg为地层饱含气时近CLYC双粒子探测器非弹热中子计数比值与远CLYC双粒子探测器非弹热中子计数比值之间的差值,d1、d2、d3均为地层饱含水时双CLYC双粒子探测器的含气饱和度刻度系数,f1、f2、f3均为地层饱含气时双CLYC双粒子探测器的含气饱和度刻度系数;
步骤6,根据近CLYC双粒子探测器的伽马时间谱获取地层的宏观俘获截面值,结合近CLYC双粒子探测器的碳氧比值以及近、远CLYC双粒子探测器的热中子计数比,计算地层流体类型划分参数,如式(5)所示:
式中,D为地层流体类型划分参数,COR为近CLYC双粒子探测器的碳氧比值,Rth为近CLYC双粒子探测器与远CLYC双粒子探测器的热中子计数比,Σ为地层的宏观俘获截面值,Σma为地层骨架的宏观俘获截面值,A、B、C均为仪器刻度系数;
步骤7,根据地层的地层流体类型划分参数和中子孔隙度,基于MCNP数值模拟方法建立地层流体类型划分图版,确定地层内流体的类型,综合地层的含油饱和度、含气饱和度和流体类型,评价地层的含油气性。
2.根据权利要求1所述的一种基于双CLYC双粒子探测器可控中子源多谱测井仪的测井方法,其特征在于,所述步骤7中,具体包括以下步骤:
步骤7.1,根据双CLYC双粒子探测器可控中子源多谱测井仪的结构参数,结合双CLYC双粒子探测器可控中子源多谱测井仪在井中的实际探测情况,构建双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型;
步骤7.2,依次将双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型中地层设置为饱含水、饱含油和饱含气的砂岩地层,分别改变双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型中地层的孔隙度,模拟得到不同地层条件下的双CLYC双粒子探测器的热中子时间谱、伽马能谱和伽马时间谱,分别计算各地层条件所对应的地层流体类型划分参数;
步骤7.3,依次根据饱含砂岩地层、饱含油砂岩地层和饱含气砂岩地层中各地层孔隙度对应的地层流体类型划分参数,分别确定砂岩地层饱含水、饱含油和饱含气时双CLYC双粒子探测器可控中子源多谱测井仪地层流体类型划分参数的测量响应,绘制地层流体类型划分图版;
步骤7.4,根据地层流体类型划分参数和地层中子孔隙度,利用地层流体类型划分图版确定地层内流体的类型,综合地层的含油饱和度、含气饱和度和流体类型,评价地层的含油气性。
3.根据权利要求2所述的一种基于双CLYC双粒子探测器可控中子源多谱测井仪的测井方法,其特征在于,所述步骤7.2中,依次将双CLYC双粒子探测器可控中子源多谱测井仪数值模拟模型中的地层孔隙度设置为0、5%、10%、15%、20%、25%、30%。
4.根据权利要求1所述的一种基于双CLYC双粒子探测器可控中子源多谱测井仪的测井方法,其特征在于,所述D-T可控中子源的产额不低于1×108/s;所述仪器壳体的直径设置为43~90mm;所述近CLYC双粒子探测器的直径设置为20~60mm、长度设置为40~60mm,与D-T可控中子源之间的距离设置为30~35cm;所述远CLYC双粒子探测器的直径设置为20~60mm、长度设置为100~140mm,与D-T可控中子源之间的距离设置为50~55cm;所述近CLYC双粒子探测器和远CLYC双粒子探测器采集的伽马能谱能量范围为0.1~10.3MeV,分为256道;
所述近CLYC双粒子探测器和远CLYC双粒子探测器均为Cs2LiYCl6:Ce3+闪烁晶体探测器,Cs2LiYCl6:Ce3+晶体的密度为3.31g/cm3,伽马能量分辨率为4~5%,热中子能量分辨率为4~6%。
5.根据权利要求1所述的一种基于双CLYC双粒子探测器可控中子源多谱测井仪的测井方法,其特征在于,所述第一钨镍铁屏蔽体的长度设置为7~12cm,第二钨镍铁屏蔽体的长度设置为1~2cm,第一钨镍铁屏蔽体和第二钨镍铁屏蔽体的直径均设置为40~90mm,第一钨镍铁屏蔽体和第二钨镍铁屏蔽体中均含有1%的硼元素。
6.根据权利要求1所述的一种基于双CLYC双粒子探测器可控中子源多谱测井仪的测井方法,其特征在于,所述可控中子源多谱测井仪脉冲测量时序的周期时长为100ms,包括时长为99ms的时间段T1和时长为1ms的时间段T2;时间段T1内包括4个记录周期,各记录周期内包括300个碳氧比模式记录周期和1个中子寿命模式记录周期,碳氧比模式记录周期的总时长为75μs,包括时间门G1和时间门G2,时间门G1的记录范围为第10~30μs,用于控制D-T可控中子源以脉冲形式发射快中子以及记录碳氧比模式总伽马能谱,时间门G2的记录范围为第35~75μs,用于记录碳氧比模式俘获伽马能谱和碳氧比模式伽马时间谱,中子寿命记录周期的总时长为2250μs,其中设置有时间门G3,时间门G3的记录范围为第75~2325μs,用于记录中子寿命模式伽马时间谱和热中子时间谱,时间段T2内设置有时间门G4,时间门G4的记录范围为2325~3325μs,用于记录本底伽马时间谱。
CN202210643863.5A 2022-06-08 2022-06-08 基于双clyc双粒子探测器的可控中子源多谱测井仪及测井方法 Active CN115012920B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210643863.5A CN115012920B (zh) 2022-06-08 2022-06-08 基于双clyc双粒子探测器的可控中子源多谱测井仪及测井方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210643863.5A CN115012920B (zh) 2022-06-08 2022-06-08 基于双clyc双粒子探测器的可控中子源多谱测井仪及测井方法

Publications (2)

Publication Number Publication Date
CN115012920A CN115012920A (zh) 2022-09-06
CN115012920B true CN115012920B (zh) 2024-07-09

Family

ID=83072988

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210643863.5A Active CN115012920B (zh) 2022-06-08 2022-06-08 基于双clyc双粒子探测器的可控中子源多谱测井仪及测井方法

Country Status (1)

Country Link
CN (1) CN115012920B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115291288B (zh) * 2022-09-29 2022-12-30 中石化经纬有限公司 基于铁中子标记的随钻脉冲中子孔隙度智能处理方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201835828U (zh) * 2010-07-23 2011-05-18 中国石油化工集团公司 双脉冲中子伽马谱测井仪

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5900627A (en) * 1997-06-19 1999-05-04 Computalog Research, Inc. Formation density measurement utilizing pulse neutrons
CN101915093A (zh) * 2010-07-23 2010-12-15 中国石油化工集团公司 一种采集双脉冲中子伽马谱的时序设计方法
CN102518431B (zh) * 2011-12-26 2015-04-22 中国石油大学(华东) 基于可控中子源的随钻多参数测井方法
US10261214B2 (en) * 2016-05-02 2019-04-16 Schlumberger Technology Corporation Method and apparatus for separating gamma and neutron signals from a radiation detector and for gain-stabilizing the detector
CN107288629B (zh) * 2017-07-25 2020-07-14 中国石油大学(华东) 一种基于新型n-γ双粒子探测器的中子伽马密度测井方法
US10852283B2 (en) * 2017-12-14 2020-12-01 Schlumberger Technology Corporation Systems and methods for detecting a neutron generator operating in air
WO2019165452A1 (en) * 2018-02-26 2019-08-29 Starfire Industries Llc Neutron-gamma imaging system to render azimuthal images in oil and gas well structures
CN111335886B (zh) * 2020-02-06 2023-04-28 长江大学 一种中子伽马密度测井测量装置及方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201835828U (zh) * 2010-07-23 2011-05-18 中国石油化工集团公司 双脉冲中子伽马谱测井仪

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于CLYC探测器的可控源中子孔隙度测井数值模拟研究;周悦;《核技术》;20210415;第44卷(第4期);74-79 *

Also Published As

Publication number Publication date
CN115012920A (zh) 2022-09-06

Similar Documents

Publication Publication Date Title
US7365307B2 (en) Sigma/porosity tools with neutron monitors
CN108643890B (zh) 一种确定地层孔隙度的方法
CN103696765B (zh) 基于可控中子源的双LaBr3探测器元素能谱测井仪及测井方法
CN115012920B (zh) 基于双clyc双粒子探测器的可控中子源多谱测井仪及测井方法
CN103470252B (zh) 基于超热中子时间谱的瞬发中子测井及铀矿定量方法
CN112904436B (zh) 一种元素产额与热中子计数比组合的孔隙度测量方法
CN110454147B (zh) 一种可控源一体化核测井仪及测井方法
CN115267928B (zh) 一种用于随钻元素测井的能谱智能处理方法
CN112523741B (zh) 一种基于能谱测井交叉谱段的铀矿定量刻度系数求法
CN201137491Y (zh) 脉冲中子双谱流体饱和度测井仪
CN112392454B (zh) 一种基于中子活化分析自屏蔽校正的裂缝体密度定量计算方法
CN110469324B (zh) 一种基于脉冲中子测井的计算地层密度方法
CN115291288B (zh) 基于铁中子标记的随钻脉冲中子孔隙度智能处理方法
CN115267930A (zh) 一种基于d-t脉冲中子源的高灵敏度中子孔隙度测量方法
CN111123379B (zh) 一种基于双谱联合的纯净非弹伽马能谱获取方法
CN112523742B (zh) 一种基于能谱测井特征谱段求取自然伽马放射性元素含量的方法
CN112034524B (zh) 一种双探测器井中地层俘获截面计算方法
CN111638559B (zh) 一种基于脉冲中子测井的快中子散射截面表征方法
CN112630848A (zh) 一种基于能谱测井特征谱段的铀矿定量剥离系数求法
CN115248463B (zh) 一种d-t源中子孔隙度测井非弹性散射影响的校正方法
CN116025335A (zh) 基于钻铤次生非弹和俘获伽马信息交会的气层识别方法
CN113902348B (zh) 一种基于孔内伽玛测试快速评价隧道工程放射性风险等级的方法
CN115614022B (zh) 铀裂变瞬发中子测井的井眼水层原位修正方法
CN116357301A (zh) 一种基于有效反应截面的地层孔隙度计算方法和系统
CN215949458U (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