CN104329075B - 元素测井中元素俘获标准谱的获取方法 - Google Patents

元素测井中元素俘获标准谱的获取方法 Download PDF

Info

Publication number
CN104329075B
CN104329075B CN201410453115.6A CN201410453115A CN104329075B CN 104329075 B CN104329075 B CN 104329075B CN 201410453115 A CN201410453115 A CN 201410453115A CN 104329075 B CN104329075 B CN 104329075B
Authority
CN
China
Prior art keywords
model
card
capture
test pit
test
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
CN201410453115.6A
Other languages
English (en)
Other versions
CN104329075A (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.)
Xi'an Austrian Electronic Instrument Ltd By Share Ltd
CNPC Bohai Drilling Engineering Co Ltd
Original Assignee
Xi'an Austrian Electronic Instrument Ltd By Share Ltd
CNPC Bohai Drilling Engineering Co Ltd
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 Xi'an Austrian Electronic Instrument Ltd By Share Ltd, CNPC Bohai Drilling Engineering Co Ltd filed Critical Xi'an Austrian Electronic Instrument Ltd By Share Ltd
Priority to CN201410453115.6A priority Critical patent/CN104329075B/zh
Publication of CN104329075A publication Critical patent/CN104329075A/zh
Application granted granted Critical
Publication of CN104329075B publication Critical patent/CN104329075B/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
    • E21B47/00Survey of boreholes or wells
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mining & Mineral Resources (AREA)
  • Geology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Fluid Mechanics (AREA)
  • Geochemistry & Mineralogy (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Geophysics (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明提供一种元素测井中元素俘获标准谱的获取方法,建立模型、校验模型等步骤,其核心是按照碳氧比饱和度刻度井的条件建立仿真模型,然后使用建立好的仿真模型来计算得到元素俘获标准谱。本发明的初始仿真模型是根据刻度井以及测井仪器的设计尺寸建立的,同时对测井仪器的某些结构进行了合理优化,减小了建模的难。保温瓶的内壳、外壳中间是真空,我们将保温瓶的内壳、外壳合并到一起,将仪器的真空部分合在一起。建立模型需要的MCNP输入卡包括栅元卡、曲面卡、材料卡、计数卡四个部分。本发明无需建立大量的刻度井,只需在少量的刻度井内进行常规的刻度实验,成本很低。

Description

元素测井中元素俘获标准谱的获取方法
技术领域
本发明涉及元素测井中一种元素俘获标准谱的获取方法。
背景技术
目前,元素测井中使用的元素俘获标准谱大多是采用仪器在填充纯净物质或者是某种元素的氧化物的刻度井中测量得到。测量元素俘获标准谱需要使用中子源将刻度井中的待测量物质活化,使其放出伽马射线,使用仪器测量这些物质的俘获伽马能谱作为元素测井中使用的元素俘获标准谱。元素测井使用的元素俘获标准谱包括硅、钙、硫、铁、钛、钆、氢、氯等多种元素,用测量的方法获取这些元素俘获标准谱需要建立大量的刻度井,建设周期长,建设成本昂贵,具有一定的危险性,不利于元素测井的推广。
发明内容
本发明提出一种使用计算机仿真软件来获取元素测井中元素俘获标准谱的方法,其按照刻度井实际的条件(测井仪器、井眼条件、地层条件等)来建立仿真模型,使用蒙特卡罗数值模拟软件MCNP5来计算硅、钙、硫、铁、钛、钆、氢、氯等多种元素的标准谱。
本发明的技术方案:
该元素测井中元素俘获标准谱的获取方法,包括以下步骤:
1)建立模型:
根据刻度井以及测井仪器的设计尺寸建立初始仿真模型,同时将整个模型中的空气部分置成真空;
上述初始仿真模型的MCNP输入卡包括栅元卡、曲面卡、材料卡、计数卡四个部分;
上述栅元卡描述仪器的各个部件之间的空隙、井液、套管、水泥环、地层、各个部件的结构信息,各个部件的结构信息包括晶体类型、探头封装、保温瓶、仪器外壳;
上述曲面卡描述围成各个部件的面的位置以及尺寸信息;
上述材料卡记录仪器的各个部件的材料成分;
上述计数卡记录模拟计算中使用的放射源的类型、能量、强度以及最终的统计类型;
2)校验模型:
上述校验模型包括校正和验证的步骤;
上述校正是指在刻度井中进行实验来对初始模型进行校正,用于减小测井仪器的实际尺寸与设计尺寸的差异;上述刻度井为水+砂岩刻度井;
上述验证是指对于校正好的模型还需要在多口刻度井中进行实验,利用实验数据与仿真得到的数据进行对比,从而验证校正过的模型是否正确;上述多口刻度井分别为油+砂岩刻度井、水+灰岩刻度井、纯水刻度井;
上述校正或验证的具体步骤如下:
(1)测井仪器下放到刻度井里贴一侧井壁,使用定点测量的方式以保证中子源处于刻度井的中心位置;
(2)测井仪器通电检查,仪器正常后开始测量;
(3)将测量得到的伽马俘获能谱与建立好的初始仿真模型计算得到的伽马俘获能谱放在一起对比,对仿真模型进行校验;
若实测的伽马俘获能谱与仿真的俘获能谱的铁峰不一致,则模型中的铁的厚度需要调整;
若实测伽马俘获能谱的氢峰与仿真的俘获伽马能谱的氢峰不一致,则模型中的孔隙度需要调整;
验证的标准是:在感兴趣的能量段内,实测谱与计算谱的残差绝对值的和△小于设定值;
3)计算仿真的元素俘获标准谱:
使用已经校验过的仿真模型,计算元素俘获标准谱;
4)计算仿真其它元素的俘获标准谱对模型中的探测器、屏蔽体不做改动,地层材料使用元素俘获标准谱中的某种元素的单质来替换,其他部分的材料用单质氧来替换,模型的尺寸、密度不要改变,依此模型来计算其它某种元素俘获标准谱。
上述设定值为0.1。
本发明的优点:
1、成本低。本发明无需建立大量的刻度井,只需在少量的刻度井内进行常规的刻度实验,成本很低。
2、速度快。本发明可采用多台计算机并行计算,大大缩减了标准谱的获取周期。
3、安全性高。本发明只需要进行少量的中子源实验,对环境、实验人员的伤害小。
附图说明
图1是本发明的流程图;
图2是在砂岩刻度井中实测的俘获伽马能谱;
图3是校正之后仿真模型计算得到的俘获伽马能谱;
图4是两者的对比图,图中的横坐标是伽马射线的能量单位是MeV,纵坐标是伽马射线计数率经过归一化之后的结果;其中,带●的曲线为实测砂岩的能谱,带▲的曲线为模拟计算得到的砂岩能谱。
具体实施方式
本发明的核心是按照碳氧比饱和度刻度井的条件建立仿真模型,然后使用建立好的仿真模型来计算得到元素俘获标准谱。
1、建立模型。本发明的初始仿真模型是根据刻度井以及测井仪器的设计尺寸建立的,同时对测井仪器的某些结构进行了合理优化,减小了建模的难。例如高能中子、光子在空气中运动轨迹、能量几乎不发生改变,所以我们将整个模型中的空气部分置成真空,这样能节约建模的时间同时减少MCNP运行的机时;保温瓶的内壳、外壳中间是真空,我们将保温瓶的内壳、外壳合并到一起,将仪器的真空部分合在一起。建立模型需要的MCNP输入卡包括栅元卡、曲面卡、材料卡、计数卡四个部分。具体如下:
C栅元卡
C曲面卡
C计数卡+材料卡
MODE N P
IMP:N 1 12r 0
IMP:P 1 12r 0
SDEF erg=14par=1pos=0 0-27.5axs=0 0 1rad=d1ext=d2cel=8
si1 0 1.97
si2 5
$空气
2、校验模型。校验模型分为校正和验证两个部分,由于刻度井在建设的过程中,不能保证与图纸设计的尺寸完全一致,包括测井仪器的实际尺寸总会与设计尺寸略有差异,虽然差异不大,但是为了保证仿真得到的元素俘获标准谱的精度,我们通过在刻度井中进行实验来对初始模型进行校正来减小这些差异;对于校正好的模型我们还需要在多口刻度井中进行实验,利用实验数据与仿真得到的数据进行对比,从而验证校正过的模型是否正确。本发明共选用了四口刻度井来进行实验,分别为水+砂岩井、油+砂岩井、水+灰岩井、纯水井,水砂岩刻度井用来校正模型,其他三口刻度井用来验证仿真模型。
实验步骤如下:
(1)测井仪器下放到刻度井里贴一侧井壁,由于刻度井较浅,需要使用定点测量的方式,保证中子源处于刻度井的中心位置。
(2)测井仪器通电检查,仪器正常后开始测量,测量是注意观察俘获伽马能谱的谱型,保证测量结果有足够的统计精度。
将测量得到的伽马俘获能谱与建立好的初始仿真模型计算得到的伽马俘获能谱放在一起对比,对仿真模型进行校验。仪器的尺寸与设计尺寸的差异会导致实测的伽马俘获能谱与仿真的俘获能谱的铁峰不一致,刻度井的地层孔隙度的差异会导致实测伽马俘获能谱的氢峰与仿真的俘获伽马能谱的氢峰不一致,通过这两个比较明显的峰来校正仿真模型。铁峰不一致则模型中的铁的厚度需要调整,氢峰不一致则模型中的孔隙度需要调整。而验证的标准是:在感兴趣的能量段内,实测谱与计算谱的残差绝对值的和(△)小于0.1,用校正好的仿真模型分别计算俘获伽马能谱并与实测的能谱做比较,计算△,下表为各不同刻度井的结果
刻度井
砂岩刻度井 0.072094612
黄铁矿刻度井 0.051872332
灰岩刻度井 0.060769614
硬石膏刻度井 0.073072631
表1四个刻度井的△
图2是实测的砂岩刻度井内的测量俘获伽马能谱结果与仿真的俘获伽马能谱。
计算仿真的元素俘获标准谱。
使用已经校验过的仿真模型,计算元素俘获标准谱。对模型进行如下处理:模型中的探测器、屏蔽体不做改动,地层材料使用元素俘获标准谱中的某种元素的单质来替换,其他部分的材料用单质氧来替换,模型的尺寸、密度不要改变,依此模型来计算元素俘获标准谱即可。

Claims (2)

1.一种元素测井中元素俘获标准谱的获取方法,其特征在于:包括以下步骤:
1)建立模型:
根据刻度井以及测井仪器的设计尺寸建立初始仿真模型,同时将整个模型中的空气部分置成真空;
所述初始仿真模型的MCNP输入卡包括栅元卡、曲面卡、材料卡、计数卡四个部分;
所述栅元卡描述仪器的各个部件之间的空隙、井液、套管、水泥环、地层、各个部件的结构信息,各个部件的结构信息包括晶体类型、探头封装、保温瓶、仪器外壳;
所述曲面卡描述围成各个部件的面的位置以及尺寸信息;
所述材料卡记录仪器的各个部件的材料成分;
所述计数卡记录模拟计算中使用的放射源的类型、能量、强度以及最终的统计类型;
2)校验模型:
所述校验模型包括校正和验证的步骤;
所述校正是指在刻度井中进行实验来对初始模型进行校正,用于减小测井仪器的实际尺寸与设计尺寸的差异;所述刻度井为水+砂岩刻度井;
所述验证是指对于校正好的模型还需要在多口刻度井中进行实验,利用实验数据与仿真得到的数据进行对比,从而验证校正过的模型是否正确;所述多口刻度井分别为油+砂岩刻度井、水+灰岩刻度井、纯水刻度井;
所述校正或验证的具体步骤如下:
(1)测井仪器下放到刻度井里贴一侧井壁,使用定点测量的方式以保证中子源处于刻度井的中心位置;
(2)测井仪器通电检查,仪器正常后开始测量;
(3)将测量得到的伽马俘获能谱与建立好的初始仿真模型计算得到的伽马俘获能谱放在一起对比,对仿真模型进行校验;
若实测的伽马俘获能谱与仿真的俘获能谱的铁峰不一致,则模型中的铁的厚度需要调整;
若实测伽马俘获能谱的氢峰与仿真的俘获伽马能谱的氢峰不一致,则模型中的孔隙度需要调整;
验证的标准是:在感兴趣的能量段内,实测谱与计算谱的残差绝对值的和△小于设定值;
3)计算仿真的元素俘获标准谱:
使用已经校验过的仿真模型,计算元素俘获标准谱;
4)计算仿真其它元素的俘获标准谱
对模型中的探测器、屏蔽体不做改动,地层材料使用元素俘获标准谱中的某种元素的单质来替换,其他部分的材料用单质氧来替换,模型的尺寸、密度不要改变,依此模型来计算其它某种元素俘获标准谱。
2.根据权利要求1所述元素测井中元素俘获标准谱的获取方法,其特征在于:所述设定值为0.1。
CN201410453115.6A 2014-09-05 2014-09-05 元素测井中元素俘获标准谱的获取方法 Active CN104329075B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410453115.6A CN104329075B (zh) 2014-09-05 2014-09-05 元素测井中元素俘获标准谱的获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410453115.6A CN104329075B (zh) 2014-09-05 2014-09-05 元素测井中元素俘获标准谱的获取方法

Publications (2)

Publication Number Publication Date
CN104329075A CN104329075A (zh) 2015-02-04
CN104329075B true CN104329075B (zh) 2017-01-18

Family

ID=52403869

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410453115.6A Active CN104329075B (zh) 2014-09-05 2014-09-05 元素测井中元素俘获标准谱的获取方法

Country Status (1)

Country Link
CN (1) CN104329075B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105093343B (zh) 2015-08-06 2018-04-06 中国石油天然气股份有限公司 一种地层元素俘获伽马能谱的测量装置
CN106958443B (zh) * 2017-05-02 2020-10-09 长江大学 一种碳氧比测井钻井液侵入影响因素的校正方法及***
CN107238856B (zh) * 2017-05-26 2020-09-25 兰州大学 一种高通量氘-氚中子发生器中子平均能量的确定方法
CN108333637B (zh) * 2017-12-14 2020-02-04 兰州大学 一种提高元素测井技术确定元素含量准确度的方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5406078A (en) * 1992-05-28 1995-04-11 Halliburton Logging Services, Inc. Induced gamma ray spectorscopy well logging system
US5608215A (en) * 1994-09-16 1997-03-04 Schlumberger Technology Corporation Method and apparatus for determining density of earth formations
CN1756967A (zh) * 2002-12-31 2006-04-05 施蓝姆伯格海外股份有限公司 井下频谱处理方法和设备
CN101680965A (zh) * 2007-12-31 2010-03-24 普拉德研究及开发股份有限公司 用于使用锶测量勘探地质地层的方法和设备
CN201835828U (zh) * 2010-07-23 2011-05-18 中国石油化工集团公司 双脉冲中子伽马谱测井仪
CN102084271A (zh) * 2008-06-25 2011-06-01 普拉德研究及开发股份有限公司 由核谱学确定绝对元素浓度
CN103696765A (zh) * 2013-11-06 2014-04-02 中国石油大学(华东) 基于可控中子源的双LaBr3探测器元素能谱测井仪及测井方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8476584B2 (en) * 2011-08-31 2013-07-02 Baker Hughes Incorporated Methods for sourceless density downhole measurement using pulsed neutron generator

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5406078A (en) * 1992-05-28 1995-04-11 Halliburton Logging Services, Inc. Induced gamma ray spectorscopy well logging system
US5608215A (en) * 1994-09-16 1997-03-04 Schlumberger Technology Corporation Method and apparatus for determining density of earth formations
CN1756967A (zh) * 2002-12-31 2006-04-05 施蓝姆伯格海外股份有限公司 井下频谱处理方法和设备
CN101680965A (zh) * 2007-12-31 2010-03-24 普拉德研究及开发股份有限公司 用于使用锶测量勘探地质地层的方法和设备
CN102084271A (zh) * 2008-06-25 2011-06-01 普拉德研究及开发股份有限公司 由核谱学确定绝对元素浓度
CN201835828U (zh) * 2010-07-23 2011-05-18 中国石油化工集团公司 双脉冲中子伽马谱测井仪
CN103696765A (zh) * 2013-11-06 2014-04-02 中国石油大学(华东) 基于可控中子源的双LaBr3探测器元素能谱测井仪及测井方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
元素俘获能谱测井标准谱测量与模拟;武宏亮 等;《中国石油学会第十七届测井年会论文集 》;20110901;第266-275页 *
地层元素能谱测井标准伽马谱响应特性研究;张锋 等;《第十八届测井年会论文集》;20130901;第229-236页 *

Also Published As

Publication number Publication date
CN104329075A (zh) 2015-02-04

Similar Documents

Publication Publication Date Title
CN104329075B (zh) 元素测井中元素俘获标准谱的获取方法
Kang et al. In-situ stress measurements and stress distribution characteristics in underground coal mines in China
CN105093343B (zh) 一种地层元素俘获伽马能谱的测量装置
CN105426659B (zh) 获取不同燃耗下有效增殖因子对截面的灵敏度系数的方法
CN105628549B (zh) 一种煤层可解吸瓦斯含量的直接快速全自动测定方法
CN108457640A (zh) 融合瞬发中子时间谱修正自然γ总量的铀矿测井定量方法
CN107288629B (zh) 一种基于新型n-γ双粒子探测器的中子伽马密度测井方法
CN106503357A (zh) 基于穿层钻孔瓦斯涌出特征测定有效抽采半径的方法
Tian et al. A study of the principles and methods of quick validation of the outburst-prevention effect in the process of coal uncovering
Grieger et al. Neutron flux and spectrum in the Dresden Felsenkeller underground facility studied by moderated He 3 counters
CN105182422A (zh) 一种地层元素测井元素标准谱制作与应用方法
CN110145301A (zh) 一种随钻密度成像测井仪环境校正方法
CN108825219A (zh) 融合自然γ能谱与中子时间谱的铀矿测井刻度参数求法
Li et al. A study on the differences in radon exhalation of different lithologies at various depths and the factors influencing its distribution in northern Shaanxi, China
CN104074501A (zh) 一种裂变中子测井修正方法
CN112392454B (zh) 一种基于中子活化分析自屏蔽校正的裂缝体密度定量计算方法
CN109738937A (zh) 一种铀矿放射性测量仪器核查体源、其制备方法及应用
CN103711479B (zh) 基于超热中子与热中子比值的铀裂变瞬发中子测井技术
Zhou et al. Study on the influence of temperature and humidity on radon exhalation from a radon-containing solution
Doron et al. MCNP benchmarking of an inelastic neutron scattering system for soil carbon analysis
CN102536199B (zh) 一种可控源随钻核测井仪器刻度装置及刻度方法
Xie et al. Calibration and unfolding of the pulse height spectra of liquid scintillator-based neutron detectors using photon sources
CN209674006U (zh) 一种铀矿放射性测量仪器的核查体源装置
CN108222927A (zh) 一种基于x射线源的密度测井方法
CN107152277A (zh) 一种碳氧比测井计算剩余油饱和度的方法及***

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: 710061 Shaanxi province Xi'an Cuihua Road No. 254 Jialin building room 207

Applicant after: Xi'an Austrian electronic instrument Limited by Share Ltd

Applicant after: CNPC Bohai Drilling Engineering Co., Ltd.

Address before: 710061 Shaanxi province Xi'an Cuihua Road No. 254 Jialin building room 207

Applicant before: Xi'an Aohua Electronic Instrument Co., Ltd.

Applicant before: CNPC Bohai Drilling Engineering Co., Ltd.

COR Change of bibliographic data
C14 Grant of patent or utility model
GR01 Patent grant