CN102094639A - 重晶石泥浆井自然伽马能谱测井校正的方法 - Google Patents
重晶石泥浆井自然伽马能谱测井校正的方法 Download PDFInfo
- Publication number
- CN102094639A CN102094639A CN2009102470138A CN200910247013A CN102094639A CN 102094639 A CN102094639 A CN 102094639A CN 2009102470138 A CN2009102470138 A CN 2009102470138A CN 200910247013 A CN200910247013 A CN 200910247013A CN 102094639 A CN102094639 A CN 102094639A
- Authority
- CN
- China
- Prior art keywords
- energy window
- window
- energy
- 85api
- spectrum
- 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.)
- Granted
Links
Images
Landscapes
- Measurement Of Radiation (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
重晶石泥浆井自然伽马能谱测井校正的方法。目前国内外普遍采用二种方法解谱求地层钾、铀、钍含量,即,五能窗最小二乘法;三能窗剥谱法或三能窗矩阵法,其测量钾、铀、钍精度低。1)下放测井,当GR>85API,层厚>5m时,进行五能窗测量;2)进行三能窗解谱,利用三能窗解谱求获取二低能窗校正系数;3)上提测井,求取五能窗计数率,修正二能校窗校正系数;4)判定总伽马值GR是否大于85API?当总伽马值GR大于85API时,返回用第2)步;当总伽马值GR小于85API时,五能窗解谱用原系数;5)利用五能窗解谱,得出地层K、U、Th。本发明能正确判断地层岩性,测量钾、铀、钍精度高。
Description
技术领域
本发明涉及重晶石泥浆井自然伽马能谱测井校正的方法。
背景技术
目前国内外普遍采用二种方法解谱求地层钾、铀、钍含量,即,五能窗最小二乘法;三能窗剥谱法(或三能窗矩阵法)
1)五能窗最小二乘法解谱,
采用最小二乘法拟合,解出地层K、U、Th的含值,得到:
Y=(ATWA)-1ATWN…………………………(1)
式中:
地层K、U、Th的量值矩阵为:
各能窗标准计数率矩阵为:
各能窗标准计数率权重矩阵W为:
五能窗解谱的响应关系距阵
则,计算钾、铀、钍量值的误差为:
其中,V=(ATWA)-1
m------能窗数,对于五能窗,m=5
五能窗解谱比三能窗解谱多二个自由度,其解析误差降1.5~2.0倍,即五能谱解比三能窗解谱拥有更高精确度
2)能窗的能量区间的影响
(1)能窗的能量区间
(2)重晶石的影响
重晶石对五能窗中低能窗计数率降低;
重晶石对三能窗计数率影响小。
发明内容
本发明的目的是避免上述缺点,提供重晶石泥浆井自然伽马能谱测井校正的新方法,其测量精度高。
本发明的技术方案:重晶石泥浆井自然伽马能谱测井校正的方法,其特征在于:所述方法的步骤为:
1>、下放测井,当GR>85API,层厚>5m时,进行五能窗测量;
2>、进行三能窗解谱,利用三能窗解谱求获取二能窗校正系数;
3>、上提测井,求取五能窗计数率,修正二能校窗校正系数;
4>、判定总伽马值GR是否大于85API?
当总伽马值GR大于85API时,返回用第2>步;
当总伽马值GR小于85API时,五能窗解谱用原系数;
5>、利用五能窗解谱,得出地层K、U、Th含量。
所述的测量出五能窗的平均计数率N1、N2、N3、N4、N5,利用N3、N4、N5求出K、U、Th,然后求N1、N2能窗校正系数b1、b2。
所述的测量五能窗计数率N1、N2、N3、N4、N5,对N1、N2进行修正,N1=b1*N1;N2=b2*N2。
本发明的优点是测量精度高。
附图说明
以下通过优选实施例的附图和描述对本发明作进一步的说明。
图1为本发明***程序流程图;
具体实施方式
1)发明原理:采用三/五能窗综合解谱
假设,在一段井段井中重晶石含量一样.
(1)在泥岩地层段,具有高的钾、铀、钍含量,各能窗具有较高计数率,测量五能窗计数率N1、N2、N3、N4、N5。
(2)采用三能窗(N3、N4、N5)求出地层钾、铀、钍含量,
Y=(BTB)-1BTNt
三能窗解谱的响应关系距阵
三能窗标准计数率矩阵为:
利用三能窗解谱求出的K`、U`、Th`,反求五能窗二低能窗计数率n1,n2:
n1=a11*K`+a21*U`+a31*Th`
n2=a12*K`+a22*U`+a31*Th`
利用测量的五能窗中低能窗计数率N1,N2,求出校正系数b1,b2
b1=n1/N1;b2=n2/N2
(3)在全井段,采用五能窗解谱求出地层的K、U、Th,
Y=(ATW`A)-1ATW`NL,
其中,
各能窗计数率权重矩阵W`为:
2),本专利的五能窗计数率实时校正解谱方法,既消除了三能窗解谱精度低,又大大减少了重晶石泥浆井中重晶石对解谱的影响,提高了地层的K、U、Th测量精度。
参见图1,工作过程
(1)下放测井时,选择总伽马值GR大于85API的大段泥岩层,测量五能窗的平均计数率N1、N2、N3、N4、N5;
三能窗解谱求出K、U、Th,然后求N1、N2能窗校正系数b1,b2;
满足上述条件,不断修正b1,b2。
(2)上提测井时,测量五能窗计数率N1、N2、N3、N4、N5;对N1、N2进行修正,N1=b1*N1;N2=b2*N2;
对五能窗解谱求出K、U、Th;
(3)同时,判定总伽马值GR是否大于85API?
当,总伽马值GR大于85API时,返回用三能窗解谱求取新的校正系数b1,b2;
当,总伽马值GR小于85API时,五能窗解谱用原系数。
Claims (3)
1.重晶石泥浆井自然伽马能谱测井校正的方法,其特征在于:所述方法的步骤为:
1>、下放测井,当GR>85API,层厚>5m时,进行五能窗测量;
2>、进行三能窗解谱,利用三能窗解谱求获取二能窗校正系数;
3>、上提测井,求取五能窗计数率,修正二能校窗校正系数;
4>、判定总伽马值GR是否大于85API?
当总伽马值GR大于85API时,返回用第2>步;
当总伽马值GR小于85API时,五能窗解谱用原系数;
5>、利用五能窗解谱,得出地层K、U、Th含量。
2.根据权利要求1所述的重晶石泥浆井自然伽马能谱测井校正的方法,其特征在于:所述的测量出五能窗的平均计数率N1、N2、N3、N4、N5,利用N3、N4、N5求出K、U、Th,然后求N1、N2能窗校正系数b1、b2。
3.根据权利要求1所述的重晶石泥浆井自然伽马能谱测井校正的方法,其特征在于:所述的测量五能窗计数率N1、N2、N3、N4、N5,对N1、N2进行修正,N1=b1*N1;N2=b2*N2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102470138A CN102094639B (zh) | 2009-12-14 | 2009-12-14 | 重晶石泥浆井自然伽马能谱测井校正的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102470138A CN102094639B (zh) | 2009-12-14 | 2009-12-14 | 重晶石泥浆井自然伽马能谱测井校正的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102094639A true CN102094639A (zh) | 2011-06-15 |
CN102094639B CN102094639B (zh) | 2013-12-04 |
Family
ID=44127887
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009102470138A Expired - Fee Related CN102094639B (zh) | 2009-12-14 | 2009-12-14 | 重晶石泥浆井自然伽马能谱测井校正的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102094639B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103730173A (zh) * | 2013-11-26 | 2014-04-16 | 北京智崴科技发展有限公司 | 一种判断核电站反应堆下泄管道中裂变产物来源的方法 |
CN104121016A (zh) * | 2013-04-28 | 2014-10-29 | 核工业北京地质研究院 | 一种用于铀矿勘探的伽玛能谱测井仪校正方法 |
CN105181719A (zh) * | 2015-10-19 | 2015-12-23 | 核工业理化工程研究院 | 一种溶液中铀浓度在线测量方法 |
CN105425309A (zh) * | 2015-11-17 | 2016-03-23 | 山东科技大学 | 确定高伽马砂岩凝灰质含量的方法 |
CN110344817A (zh) * | 2019-02-21 | 2019-10-18 | 中陕核工业集团地质调查院有限公司 | 石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用 |
CN112649888A (zh) * | 2020-11-24 | 2021-04-13 | 东华理工大学 | 一种基于能谱测井特征谱峰的铀矿定量刻度系数求法 |
CN113250686A (zh) * | 2021-05-26 | 2021-08-13 | 核工业北京地质研究院 | 一种井下伽马能谱直接测铀方法及*** |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3858037A (en) * | 1973-01-18 | 1974-12-31 | Schlumberger Technology Corp | Well-logging measuring apparatus and method |
EP0070943A1 (en) * | 1981-07-31 | 1983-02-09 | Societe De Prospection Electrique Schlumberger | Natural gamma ray spectrum analysis and filtering technique |
GB2199139A (en) * | 1986-12-22 | 1988-06-29 | Exxon Production Research Co | Method for detecting drilling fluid in the annulus of a cased wellbore |
CN1136636A (zh) * | 1995-04-05 | 1996-11-27 | 江汉石油管理局测井研究所 | 激发极化和自然电位组合测井仪及解释方法 |
CN1243959A (zh) * | 1999-06-24 | 2000-02-09 | 佟广新 | 同位素吸水剖面测井方法 |
CN101037941A (zh) * | 2006-03-17 | 2007-09-19 | 陈为民 | 钻进式井壁取芯器 |
CN101126815A (zh) * | 2006-08-17 | 2008-02-20 | 中国石油天然气股份有限公司 | 一种利用地震岩性因子和岩性阻抗进行油气检测的方法 |
CN101414013A (zh) * | 2007-10-17 | 2009-04-22 | 中国石油天然气股份有限公司 | 一种利用地震资料确定地下流体的方法 |
-
2009
- 2009-12-14 CN CN2009102470138A patent/CN102094639B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3858037A (en) * | 1973-01-18 | 1974-12-31 | Schlumberger Technology Corp | Well-logging measuring apparatus and method |
EP0070943A1 (en) * | 1981-07-31 | 1983-02-09 | Societe De Prospection Electrique Schlumberger | Natural gamma ray spectrum analysis and filtering technique |
GB2199139A (en) * | 1986-12-22 | 1988-06-29 | Exxon Production Research Co | Method for detecting drilling fluid in the annulus of a cased wellbore |
CN1136636A (zh) * | 1995-04-05 | 1996-11-27 | 江汉石油管理局测井研究所 | 激发极化和自然电位组合测井仪及解释方法 |
CN1243959A (zh) * | 1999-06-24 | 2000-02-09 | 佟广新 | 同位素吸水剖面测井方法 |
CN101037941A (zh) * | 2006-03-17 | 2007-09-19 | 陈为民 | 钻进式井壁取芯器 |
CN101126815A (zh) * | 2006-08-17 | 2008-02-20 | 中国石油天然气股份有限公司 | 一种利用地震岩性因子和岩性阻抗进行油气检测的方法 |
CN101414013A (zh) * | 2007-10-17 | 2009-04-22 | 中国石油天然气股份有限公司 | 一种利用地震资料确定地下流体的方法 |
Non-Patent Citations (2)
Title |
---|
李传伟; 王江波; 慕德俊等: "基于最小二乘法的伽马能谱解谱方法及改进", 《测井技术》 * |
郭颖文: "自然伽马能谱数据解谱方法的研究", 《国外测井技术》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104121016A (zh) * | 2013-04-28 | 2014-10-29 | 核工业北京地质研究院 | 一种用于铀矿勘探的伽玛能谱测井仪校正方法 |
CN103730173A (zh) * | 2013-11-26 | 2014-04-16 | 北京智崴科技发展有限公司 | 一种判断核电站反应堆下泄管道中裂变产物来源的方法 |
CN103730173B (zh) * | 2013-11-26 | 2016-08-17 | 北京智崴科技发展有限公司 | 一种判断核电站反应堆下泄管道中裂变产物来源的方法 |
CN105181719A (zh) * | 2015-10-19 | 2015-12-23 | 核工业理化工程研究院 | 一种溶液中铀浓度在线测量方法 |
CN105425309A (zh) * | 2015-11-17 | 2016-03-23 | 山东科技大学 | 确定高伽马砂岩凝灰质含量的方法 |
CN110344817A (zh) * | 2019-02-21 | 2019-10-18 | 中陕核工业集团地质调查院有限公司 | 石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用 |
CN110344817B (zh) * | 2019-02-21 | 2020-10-16 | 中陕核工业集团地质调查院有限公司 | 石油孔固井水泥环属性参数在老井伽玛复测系数修正中的应用 |
CN112649888A (zh) * | 2020-11-24 | 2021-04-13 | 东华理工大学 | 一种基于能谱测井特征谱峰的铀矿定量刻度系数求法 |
CN113250686A (zh) * | 2021-05-26 | 2021-08-13 | 核工业北京地质研究院 | 一种井下伽马能谱直接测铀方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN102094639B (zh) | 2013-12-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102094639B (zh) | 重晶石泥浆井自然伽马能谱测井校正的方法 | |
US11313223B2 (en) | Gamma ray spectrum unfolding method for elemental capture spectroscopy logging and device therefor | |
CN104536056B (zh) | 小口径伽马能谱测井装置及数据采集传输和自稳方法 | |
CN104062307B (zh) | 一种车载伽玛能谱大气氡影响修正方法 | |
CN104121016A (zh) | 一种用于铀矿勘探的伽玛能谱测井仪校正方法 | |
CN104298883A (zh) | 一种油气资源评价中的烃源岩产烃率图版的建立方法 | |
CN107780923B (zh) | 一种基于泥质校正的含水饱和度模型的建立、仿真方法 | |
CN105545284B (zh) | 一种随钻伽马成像数据处理方法 | |
CN104318077A (zh) | 气候变化和人类活动对河川径流变化定量分析方法 | |
CN104615843B (zh) | 一种识别隐伏火山岩型铀矿化信息的物化探综合方法 | |
CN111045072A (zh) | 一种适用于CeBr3探测器的γ能谱迭代解谱算法 | |
CN104865614A (zh) | 基于可变骨架参数的复杂储层流体识别方法 | |
CN107657119B (zh) | 提高伽马能谱数据品质的测量方法 | |
TWI479176B (zh) | 可應用於碘化鈉偵檢器量測所得能譜且具高核種鑑別力之活度獲得方法 | |
CN106569252A (zh) | 一种校正伽玛总量测量型活性炭测氡仪一致性的方法 | |
CN104111482A (zh) | 一种双探测器x射线荧光测井探管及方法 | |
TWI495896B (zh) | 碘化鈉偵檢器量測所得能譜之活度計算方法 | |
CN100492055C (zh) | 一种氯能谱测井方法 | |
CN204116253U (zh) | 一种双探测器x射线荧光测井探管 | |
CN112987117A (zh) | 基于自然γ能谱测井多特征峰组合的铀矿定量换算系数求法 | |
Chruścińska et al. | Luminescence dating of bricks from the gothic Saint James Church in Toruń | |
CN110515115A (zh) | 一种用于伽马谱峰面积求解的评估方法 | |
CN115614022B (zh) | 铀裂变瞬发中子测井的井眼水层原位修正方法 | |
CN106324692A (zh) | 一种地面伽玛能谱测量铀含量数据调平校正方法 | |
CN105806852A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20131204 Termination date: 20181214 |
|
CF01 | Termination of patent right due to non-payment of annual fee |