CN107229080B - 一种元素测井中子俘获伽马能谱的获取方法 - Google Patents

一种元素测井中子俘获伽马能谱的获取方法 Download PDF

Info

Publication number
CN107229080B
CN107229080B CN201710369388.6A CN201710369388A CN107229080B CN 107229080 B CN107229080 B CN 107229080B CN 201710369388 A CN201710369388 A CN 201710369388A CN 107229080 B CN107229080 B CN 107229080B
Authority
CN
China
Prior art keywords
gamma spectra
gamma
single element
data
neutron absorption
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.)
Expired - Fee Related
Application number
CN201710369388.6A
Other languages
English (en)
Other versions
CN107229080A (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.)
Lanzhou University
Original Assignee
Lanzhou University
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 Lanzhou University filed Critical Lanzhou University
Priority to CN201710369388.6A priority Critical patent/CN107229080B/zh
Publication of CN107229080A publication Critical patent/CN107229080A/zh
Application granted granted Critical
Publication of CN107229080B publication Critical patent/CN107229080B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V5/00Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
    • G01V5/04Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging
    • G01V5/08Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays
    • G01V5/10Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays using neutron sources
    • G01V5/101Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays using neutron sources and detecting the secondary Y-rays produced in the surrounding layers of the bore hole
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/36Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V13/00Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00

Landscapes

  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Geophysics (AREA)
  • Molecular Biology (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Measurement Of Radiation (AREA)

Abstract

本发明公开了一种元素测井中子俘获伽马能谱的获取方法,在实验室条件下测量单元素的伽马能谱并进行分析,确定探测器的能量分辨率;针对地层元素俘获伽马能谱测井仪器参数建立模型;采用蒙特卡洛方法模拟地层元素俘获伽马能谱测井仪器井下的热中子俘获伽马能谱;对模拟得到的单元素伽马能谱进行处理分析,从而确定某一种元素的单元素模拟谱。本发明的有益效果是针对地层元素测井技术中ECS测井仪建立模型,得到某一种元素的单元素模拟谱,为ECS测井仪器的数据处理提供重要技术支撑。

Description

一种元素测井中子俘获伽马能谱的获取方法
技术领域
本发明属于地球物理测井技术领域,涉及一种元素测井中子俘获伽马能谱的获取方法。
背景技术
元素俘获能谱测井仪器(简称ECS)采用的是241Am-9Be中子源,放射出中子的能谱范围为0~10MeV。在中子与地层物质的相互作用过程中,主要可分为中子慢化阶段和中子吸收阶段两个阶段。数据探测系统探测这两个阶段发出的γ射线。在中子慢化阶段,只有非弹性散射能发出γ射线。在中子吸收阶段,热中子与原子核发生俘获反应可以放射出γ射线。当发生俘获反应时,中子被靶核吸收形成复合核,而后复合核通过发射特征γ光子回到基态,该俘获反应阶段的中子能量较低。元素俘获能谱测井主要是中子与地层发生非弹性散射和俘获这两个核反应发出γ光子,探测器到探测γ射线。最后对探测到的γ射线进行分析、解谱和计算出地层中各种元素的含量。
在对元素俘获伽马能谱测井获得的原始测量谱进行解谱时需要用到单元素标准谱。所谓单元素标准谱就是一种元素与中子发生作用发出γ光子,并被探测器记录的谱线。国外斯伦贝谢公司确定单元素标准伽马能谱的方法是使用标准井实验确定,建立一个标准实验井,若需要测量某种元素的单元素谱,则将含有该元素的化合物填充在标准井内并且用探测仪测量该井得到伽马能谱,最后剥去其他干扰元素的伽马能谱就得到该种元素的单元素伽马能谱。国内目前还没有针对单元素建立的标准井,因而无法通过实际仪器测量获得标准谱,从而制约了国内相关仪器的研发和油田现场的实际应用。
发明内容
本发明的目的在于提供一种元素测井中子俘获伽马能谱的获取方法,解决了国内目前还没有针对单元素建立的标准井,因而无法通过实际仪器测量获得标准谱,从而制约了国内相关仪器的研发和油田现场的实际应用的问题。
本发明按照以下步骤进行:
步骤一:在实验室条件下测量单元素的伽马能谱并进行分析,确定探测器的能量分辨率;
步骤二:针对地层元素俘获伽马能谱测井仪器参数建立模型;
步骤三:采用蒙特卡洛方法模拟地层元素俘获伽马能谱测井仪器井下的热中子俘获伽马能谱;
步骤四:对模拟得到的单元素伽马能谱进行处理分析,从而确定某一种元素的单元素模拟谱。
进一步,步骤一中,通过搭建实验平台得到一些地层单元素样品伽马能谱,分析不同元素的伽马能谱并计算出BGO探测器的能量分辨率;
能量分辨率与能量以及半高宽的关系公式:
能谱半高宽和能谱关系式:
E是粒子能量,a的单位是MeV,b的单位是c无单位。
进一步,步骤二中模型选择241Am(ɑ,n)9Be同位素中子源;选择MCNP中的通用源模型且将其设置为体积源;测井仪模型部分是根据ECS测井仪所使用的测井仪器的几何尺寸,部件构成和其他一些相关参数建立;地层模型部分外径为160cm,内径为20cm,高为240cm的同心圆柱体。
进一步,步骤三中,根据步骤二获得的中子俘获伽马能谱数据,处理分析各元素的能量窗内的模拟伽马能谱计数,包括:1)对所述中子俘获伽马能谱数据进行预处理,获得预处理后的单元素中子俘获伽马能谱模拟数据;2)对所述预处理后的单元素中子俘获伽马能谱进行归一化处理,获得归一化处理后的单元素中子俘获伽马能谱数据;3)根据所述归一化处理后的单元素中子俘获伽马能谱数据,获取所述反映各元素的能量窗内的伽马能谱计数。
进一步,步骤四中,对模拟得到的单元素伽马能谱进行处理分析包括:对所述中子俘获伽马能谱数据进行预处理,为了得到比较平滑的数据,必须将这些毛刺去除,对所述预处理后的单元素中子俘获伽马能谱进行归一化处理,获得归一化处理后的单元素中子俘获伽马能谱数据;根据所述归一化处理后的单元素中子俘获伽马能谱数据,获取所述反映各元素的能量窗内的伽马能谱计数;将步骤二中的模拟模型地层部分加入某种元素单质或氧化物,通过步骤三、步骤四后,最终得到某一种元素的单元素模拟谱,利用不同的单元素样品重复步骤二、步骤三、步骤四,进而可以得到一套较为完整的元素俘获伽马谱。
进一步,毛刺去除采用五点滤波法,五点三次平滑公式如下:
为了消除随机误差的影响,提高实验数据的质量,取相邻的5个数据点,拟合出一条3次曲线来,然后用3次曲线上相应的位置的数据值作为滤波后结果,获得预处理后的单元素中子俘获伽马能谱模拟数据。
本发明的有益效果是针对地层元素测井技术中ECS测井仪建立模型,得到某一种元素的单元素模拟谱,为ECS测井仪器的数据处理提供重要技术支撑。
附图说明
图1是本申请实施例中地层元素测井中子俘获伽马能谱获取方法流程图;
图2是本申请一实施例中蒙特卡罗模拟模型的示意图;
图3是本申请一实施例中一套单元素伽马能谱图。
图中:301-镅铍同位素中子源,302-BGO晶体探测器,303-模拟井眼,304-地层元素测井中的地层部分。
具体实施方式
为了使本技术领域的技术人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述。下面结合附图,对本申请实施例的具体实施方式作进一步的详细说明。
参考图1,本申请实施例提供了一种地层元素俘获伽马能谱测井单元素中子俘获伽马能谱获取方法,所述方法包括:
S1、在实验室条件下测量单元素的伽马能谱并进行分析,确定探测器的能量分辨率;
本申请实施例中,地层元素测井中BGO探测器的分辨率可以通过搭建一个标准的实验平台获取,该实验平台由中子源、慢化剂、探测系统、数据采集系统和能谱分析系统组成。通过该实验平台得到一些样品伽马能谱,通过对单元素获混合样品的俘获伽马能谱分析可以计算出BGO探测器的能量分辨率,能量分辨率与能量以及半高宽的关系公式:
能谱半高宽和能谱关系式:
E是粒子能量。a的单位是MeV,b的单位是c无单位
S2、针对ECS测井仪器建立模型;
选择241Am(ɑ,n)9Be同位素中子源,该中子源具有体积小,放射性活度大,半衰期长等优点,其中还有一个最重要的优点是中子能量范围大。基于MCNP提供的几种源模型和实际测井中源的类型,本发明选择MCNP中的通用源模型且将其设置为体积源。模型主要包括两个部分,其一是测井仪部分,其二地层部分。测井仪部分是根据ECS测井仪所使用的测井仪器的几何尺寸,部件构成和其他一些相关参数建立的。地层部分根据地层的普遍参数,其外径为160cm,内径为20cm,高为240cm的同心圆柱体。这样建立的模型与实际情况相近,能够很好地模拟实际测井。
S3、采用蒙特卡洛方法模拟ECS仪器井下的热中子俘获伽马能谱;
中子与地层元素发生相互作用情况复杂(包括弹性散射﹑快中子非弹性散射以及热中子俘获反应等),一般数值方法难以解决,然而随着计算机的普及和迅速发展,蒙特卡罗能够比较逼真的描述粒子的输运过程。
根据步骤二获得的中子俘获伽马能谱数据,处理分析各元素的能量窗内的模拟伽马能谱计数,包括:1)对所述中子俘获伽马能谱数据进行预处理,获得预处理后的单元素中子俘获伽马能谱模拟数据;2)对所述预处理后的单元素中子俘获伽马能谱进行归一化处理,获得归一化处理后的单元素中子俘获伽马能谱数据;3)根据所述归一化处理后的单元素中子俘获伽马能谱数据,获取所述反映各元素的能量窗内的伽马能谱计数。
S4、对模拟得到的单元素伽马能谱进行处理分析,最终得到某一种元素的单元素模拟谱。
对模拟得到的单元素伽马能谱进行处理分析包括:
对所述中子俘获伽马能谱数据进行预处理,最主要平滑处理,为了得到比较平滑的数据,就必须将这些毛刺去除。可以用五点滤波法去除毛刺,五点三次平滑公式如下:
为了消除随机误差的影响,提高实验数据的质量。取相邻的5个数据点,可以拟合出一条3次曲线来,然后用3次曲线上相应的位置的数据值作为滤波后结果,获得预处理后的单元素中子俘获伽马能谱模拟数据;
对所述预处理后的单元素中子俘获伽马能谱进行归一化处理,归一化公式如下:
需要对单元素标准谱的积分面积为1,在离散情况下是道址所对应的纵坐标值的和为1,获得归一化处理后的单元素中子俘获伽马能谱数据;
根据所述归一化处理后的单元素中子俘获伽马能谱数据,获取所述反映各元素的能量窗内的伽马能谱计数。
将步骤S2中的模拟模型地层部分加入某种元素单质(单质为固体的元素)或氧化物(针对单质不是固体的元素),通过S3、S4步骤后,最终得到某一种元素的单元素模拟谱。利用不同的单元素样品重复S2、S3、S4步骤,进而可以得到一套较为完整的元素俘获伽马谱。
为了清楚的说明本申请实施例的有益效果,下面举例说明:
表2为本申请实施例中地层元素测井中在实验平台下计算出的BGO探测器的分辨率。
表2
图2为本申请实施例中地层元素测井蒙特卡罗模拟模型示意图,模拟计算条件为:镅铍同位素中子源301,该中子源具有体积小,放射性活度大,半衰期长等优点。BGO晶体探测器302,其密度ρ=7.31g/cm3对粒子有高的阻止能力,对γ射线的吸收系数比Na(Tl)大,是目前探测效率最高的一种闪烁体,荧光光谱范围在350~650nm。BGO比较突出的优点是对γ射线有很高的探测效率和阻止能力。模拟井眼303,其为直径为20cm,高为240cm圆柱体,根据斯伦贝谢公司ECS测井仪,井眼中填充淡水。地层元素测井中的地层部分304,其其外径为160cm,内径为20cm,高为240cm的同心圆柱体。为了得到某单元素中子俘获伽马能谱,地层部分填充某种元素单质(单质为固体的元素)或氧化物(针对单质不是固体的元素)。
图3为本申请实施例中地层元素测井蒙特卡罗模拟中子俘获伽马能谱图。将步骤S2中的模拟模型地层部分加入某种元素单质(单质为固体的元素)或氧化物(针对单质不是固体的元素),通过S3、S4步骤后,最终得到某一种元素的单元素模拟谱。重复S2、S3、S4步骤,得到一套较为完整的元素俘获伽马谱。能窗为1.2MeV~8.4MeV,对应道址选取的能窗范围是30~210道。
本发明在实验室条件下测量单元素的伽马能谱并进行分析,确定探测器的能量分辨率;针对地层元素测井技术中ECS测井仪建立模型;采用蒙特卡洛方法模拟ECS井下的热中子俘获伽马能谱;对模拟得到的单元素伽马能谱进行处理分析,最终得到某一种元素的单元素模拟谱,为后续伽马能谱的解谱提供重要技术支撑,为ECS测井仪器的数据处理提供重要技术支撑。
以上所述仅是对本发明的较佳实施方式而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施方式所做的任何简单修改,等同变化与修饰,均属于本发明技术方案的范围内。

Claims (1)

1.一种元素测井中子俘获伽马能谱的获取方法,其特征在于按照以下步骤进行:
步骤一:在实验室条件下测量单元素的伽马能谱并进行分析,确定探测器的能量分辨率;通过搭建实验平台得到一些地层单元素样品伽马能谱,分析不同元素的伽马能谱并计算出BGO探测器的能量分辨率;
能量分辨率与能量以及半高宽的关系公式:
能谱半高宽和能谱关系式:
E是粒子能量,a的单位是MeV,b的单位是c无单位;
步骤二:针对地层元素俘获伽马能谱测井仪器参数建立模型;
模型选择241Am(ɑ,n)9Be同位素中子源;选择MCNP中的通用源模型且将其设置为体积源;测井仪模型部分是根据ECS测井仪所使用的测井仪器的几何尺寸,部件构成和其他一些相关参数建立;地层模型部分外径为160cm,内径为20cm,高为240cm的同心圆柱体;
步骤三:采用蒙特卡洛方法模拟地层元素俘获伽马能谱测井仪器井下的热中子俘获伽马能谱;根据步骤二获得的中子俘获伽马能谱数据,处理分析各元素的能量窗内的模拟伽马能谱计数,包括:1)对所述中子俘获伽马能谱数据进行预处理,获得预处理后的单元素中子俘获伽马能谱模拟数据;2)对所述预处理后的单元素中子俘获伽马能谱进行归一化处理,获得归一化处理后的单元素中子俘获伽马能谱数据;3)根据所述归一化处理后的单元素中子俘获伽马能谱数据,获取所述反映各元素的能量窗内的伽马能谱计数;
步骤四:对模拟得到的单元素伽马能谱进行处理分析,从而确定某一种元素的单元素模拟谱,对模拟得到的单元素伽马能谱进行处理分析包括:对所述中子俘获伽马能谱数据进行预处理,为了得到比较平滑的数据,必须将这些毛刺去除,对所述预处理后的单元素中子俘获伽马能谱进行归一化处理,获得归一化处理后的单元素中子俘获伽马能谱数据;根据所述归一化处理后的单元素中子俘获伽马能谱数据,获取所述反映各元素的能量窗内的伽马能谱计数;将步骤二中的模拟模型地层部分加入某种元素单质或氧化物,通过步骤三、步骤四后,最终得到某一种元素的单元素模拟谱,利用不同的单元素样品重复步骤二、步骤三、步骤四,进而可以得到一套较为完整的元素俘获伽马谱;
所述毛刺去除采用五点滤波法去除毛刺,五点三次平滑公式如下:
为了消除随机误差的影响,提高实验数据的质量,取相邻的5个数据点,拟合出一条3次曲线来,然后用3次曲线上相应的位置的数据值作为滤波后结果,获得预处理后的单元素中子俘获伽马能谱模拟数据。
CN201710369388.6A 2017-05-23 2017-05-23 一种元素测井中子俘获伽马能谱的获取方法 Expired - Fee Related CN107229080B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710369388.6A CN107229080B (zh) 2017-05-23 2017-05-23 一种元素测井中子俘获伽马能谱的获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710369388.6A CN107229080B (zh) 2017-05-23 2017-05-23 一种元素测井中子俘获伽马能谱的获取方法

Publications (2)

Publication Number Publication Date
CN107229080A CN107229080A (zh) 2017-10-03
CN107229080B true CN107229080B (zh) 2018-07-20

Family

ID=59934652

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710369388.6A Expired - Fee Related CN107229080B (zh) 2017-05-23 2017-05-23 一种元素测井中子俘获伽马能谱的获取方法

Country Status (1)

Country Link
CN (1) CN107229080B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107907901B (zh) * 2017-11-06 2020-03-17 中广核工程有限公司 核电站中子、伽马能谱与剂量的测量方法和系统
CN108333637B (zh) * 2017-12-14 2020-02-04 兰州大学 一种提高元素测井技术确定元素含量准确度的方法
CN108457644B (zh) 2017-12-26 2020-04-10 中国石油天然气股份有限公司 一种元素俘获能谱测井的伽马能谱解谱方法及装置
CN111553110B (zh) * 2020-04-22 2023-03-24 西北核技术研究院 一种基于蒙特卡罗模拟的核电磁脉冲电流源项计算方法
CN112114382A (zh) * 2020-08-19 2020-12-22 中国石油天然气集团有限公司 一种地层元素测井中子俘获伽马标准谱的校正方法
CN112258432B (zh) * 2020-10-23 2021-09-21 成都大学 基于复合型结构元素数学形态学的中子-伽马甄别方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010002796A2 (en) * 2008-06-25 2010-01-07 Schlumberger Canada Limited Absolute elemental concentrations from nuclear spectroscopy
CN105093343B (zh) * 2015-08-06 2018-04-06 中国石油天然气股份有限公司 一种地层元素俘获伽马能谱的测量装置

Also Published As

Publication number Publication date
CN107229080A (zh) 2017-10-03

Similar Documents

Publication Publication Date Title
CN107229080B (zh) 一种元素测井中子俘获伽马能谱的获取方法
CN105093343B (zh) 一种地层元素俘获伽马能谱的测量装置
Dong et al. Research on recognition of gas saturation in sandstone reservoir based on capture mode
CN105022084B (zh) 一种数字化中子谱仪
CN106250619B (zh) 一种确定地层矿物含量的方法和装置
CN108825219A (zh) 融合自然γ能谱与中子时间谱的铀矿测井刻度参数求法
Tang et al. A method for neutron-induced gamma spectra decomposition analysis based on Geant4 simulation
US6577697B2 (en) Field analysis of geological samples using delayed neutron activation analysis
Tian et al. Monte Carlo simulation of Cu, Ni and Fe grade determination in borehole by PGNAA technique
Wu et al. Standard spectrum measurement and simulation of elemental capture spectroscopy log
Abd El Gawad et al. Improving the analysis performance of gamma spectrometer using the Monte Carlo code for accurate measurements of uranium samples
Feng et al. Monte Carlo simulation of PGNAA system for determining element content in the rock sample
CN112523741B (zh) 一种基于能谱测井交叉谱段的铀矿定量刻度系数求法
Yu et al. Corrections of fast neutron inelastic scattering effects on DT neutron porosity logging
Shue et al. Thermal-neutron intensities in soils irradiated by fast neutrons from point sources
Xing et al. A LWD gas-reservoir recognition method based on the inelastic gamma information of the drill collar
Esmaeili-Sani et al. Gamma–gamma density and lithology tools simulation based on GEANT4 advanced low energy Compton scattering (GALECS) package
CN112523742A (zh) 一种基于能谱测井特征谱段的铀矿定量刻度系数求法
Cai et al. Design of a DT neutron source based PGNAA facility for element determination in aqueous solution
Liu et al. A method for improving the evaluation of elemental concentrations measured by geochemical well logging
CN106018231B (zh) 一种原位探测水蚀土壤中137Cs渗透深度的方法和系统
Wu et al. Calculating standard captured γ spectra of formation elements
Islam et al. Study of a laboratory-based gamma spectrometry for food and environmental samples
Wensheng et al. Quantitative analysis of neutron-capture gamma-ray energy spectra using direct demodulation
CN115614022B (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180720

Termination date: 20200523

CF01 Termination of patent right due to non-payment of annual fee