地層元素測(cè)井儀器解譜分析確定產(chǎn)額的方法
【專(zhuān)利摘要】本發(fā)明提供地層元素測(cè)井儀器解譜分析確定產(chǎn)額的方法，包括確定段區(qū)數(shù)n及各區(qū)段左右邊界kn和kn+1、計(jì)算標(biāo)準(zhǔn)俘獲譜歸一化響應(yīng)矩陣A、計(jì)算待解譜在各段區(qū)總計(jì)數(shù)值ci、計(jì)算對(duì)角加權(quán)矩陣W、計(jì)算加權(quán)最小二乘法解譜得產(chǎn)額yj及利用殘差平方和的平均數(shù)Δ2判定擬合效果。本發(fā)明采用卡段法解譜方法。卡段法：有選擇性的選擇元素的特征能譜所在的道區(qū)解譜。通過(guò)MC數(shù)值模擬發(fā)現(xiàn)該方法對(duì)解譜沒(méi)有太大的影響。
【專(zhuān)利說(shuō)明】地層元素測(cè)井儀器解譜分析確定產(chǎn)額的方法

【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及一種對(duì)地層元素測(cè)井儀器進(jìn)行解譜以分析確定產(chǎn)額的方法。

【背景技術(shù)】
[0002] 對(duì)地層元素測(cè)井而言，一般都是采用最小二乘法或者最小二乘法的改進(jìn)方法進(jìn) 行解譜。這種解譜方法計(jì)算出來(lái)的地層元素種類(lèi)不會(huì)出錯(cuò)，但是對(duì)于某些含量較少或者是 靈敏度比較低的元素計(jì)算的過(guò)程中誤差較大。


【發(fā)明內(nèi)容】

[0003] 本發(fā)明目的是提供一種地層元素測(cè)井儀器解譜分析確定產(chǎn)額的方法，其基本上能 反映出地層元素的種類(lèi)和相對(duì)含量。
[0004] 一種地層元素測(cè)井儀器解譜分析確定產(chǎn)額的方法，其特殊之處在于：包括以下步 驟：
[0005] 1)確定段區(qū)數(shù)n及各區(qū)段左右邊界kn和kn+1 ;所述kn和kn+1分別表示第n個(gè)段區(qū) 的左右邊界；n為不大于256的整數(shù)
[0006] 2)計(jì)算標(biāo)準(zhǔn)俘獲譜歸一化響應(yīng)矩陣A
[0007] 響應(yīng)矩陣A的計(jì)算方法如下；
[0008] 設(shè)au是響應(yīng)矩陣A的（i，j)元，即A的第i行第j列的元素；
[0009]貝1J
[0010]

【權(quán)利要求】
1. 一種地層元素測(cè)井儀器解譜分析確定產(chǎn)額的方法，其特征在于：包括以下步驟： 1) 確定段區(qū)數(shù)η及各區(qū)段左右邊界kn和kn+1 ;所述kn和kn+1分別表示第η個(gè)段區(qū)的左 右邊界；η為不大于256的整數(shù) 2) 計(jì)算標(biāo)準(zhǔn)俘獲譜歸一化響應(yīng)矩陣A 響應(yīng)矩陣A的計(jì)算方法如下； 設(shè)％是響應(yīng)矩陣A的（i，j)元，即A的第i行第j列的元素； 則 (1 ?=Σ? ； kn和kn+1分別表示第η個(gè)段區(qū)的左右邊界，Nw是準(zhǔn)譜歸一化之后的第k道的計(jì)數(shù)； 3) 計(jì)算待解譜在各段區(qū)總計(jì)數(shù)值Ci 設(shè)在實(shí)際測(cè)量中得到的256道γ能譜的第i道計(jì)數(shù)為Ci ; 則 Cfh n = \, 2,... η' k-K 其中，匕和、+1分別表示第n個(gè)段區(qū)的左右邊界，NkM是實(shí)測(cè)地層譜歸一化后的第k道 計(jì)數(shù)； 4) 計(jì)算對(duì)角加權(quán)矩陣W 設(shè)對(duì)角加權(quán)陣為W，Wu是W的（i，i)元，令權(quán)系數(shù)r = Wii，則 1 wu = W) = - c； 則： wjj Ο Ο Ο Ο n Ο Ο W= Ο Ο ··· Ο Ο ο ο Wnn_ 5) 計(jì)算加權(quán)最小二乘法解譜得產(chǎn)額yj 地層測(cè)量譜的第i道的計(jì)數(shù)可以表示為，地層中所有元素在第i道內(nèi)產(chǎn)生的Y射線(xiàn)即 標(biāo)準(zhǔn)譜的線(xiàn)性疊加； 1? c- = a^yj + εi /=1、2、3......m >ι 即疊加的系數(shù)為該地層中元素的產(chǎn)額； 其中是由元素的標(biāo)準(zhǔn)譜得到的mXm，的響應(yīng)矩陣A的（i，j)元，y」表示第j種元 素的百分含量，即產(chǎn)額，表示誤差，m為元素總數(shù)，m，為總道數(shù)，采用加權(quán)最小二乘法求 解，得： y = [ATWA]_1[ATWC] 其中 C = (c" c2, · · · · cn)，y = (y" y2, · · · ym); 6)利用殘差平方和的平均數(shù)Λ2判定擬合效果； 設(shè)殘差平方和的平均數(shù)為 Δ =^￡(c,-少 11 i=\ j=l 若確定的精度Λ2〈 ε，則滿(mǎn)足要求，否則返回步驟1)重新進(jìn)行解譜。
【文檔編號(hào)】G01N33/24GK104239710SQ201410453182
【公開(kāi)日】2014年12月24日 申請(qǐng)日期:2014年9月5日 優(yōu)先權(quán)日:2014年9月5日
【發(fā)明者】石文夫, 張政彬, 譚紅建 申請(qǐng)人:西安奧華電子儀器有限責(zé)任公司