一種地層元素測井儀器解譜分析確定產(chǎn)額的方法
【專利摘要】本發(fā)明提供一種地層元素測井儀器解譜分析確定產(chǎn)額的方法，包括計算標(biāo)準(zhǔn)俘獲譜歸一化響應(yīng)矩陣A、計算實測地層俘獲譜計數(shù)ci、計算對角加權(quán)矩陣W、計算加權(quán)最小二乘法解譜得產(chǎn)額yj、利用殘差平方和的平均數(shù)Δ2判定擬合效果等步驟，本發(fā)明采用逐道法解譜方法。采用完整的256道能窗解譜。通過M C數(shù)值模擬發(fā)現(xiàn)該方法對解譜沒有太大的影響。利用逐道法對未知地層進(jìn)行解譜，發(fā)現(xiàn)解析的結(jié)果誤差還較大，但是基本上能反映出地層元素的種類和相對含量。
【專利說明】一種地層元素測井儀器解譜分析確定產(chǎn)額的方法

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

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


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

[0003] 本發(fā)明目的是提供一種地層元素測井儀器解譜分析確定產(chǎn)額的方法，其基本上能 反映出地層元素的種類和相對含量。
[0004] 一種地層元素測井儀器解譜分析確定產(chǎn)額的方法，其特殊之處在于：包括以下步 驟：
[0005] 1)計算標(biāo)準(zhǔn)俘獲譜歸一化響應(yīng)矩陣A
[0006] 響應(yīng)矩陣A是用加權(quán)最小二乘法解譜時的設(shè)計矩陣；
[0007] 設(shè)au是響應(yīng)矩陣A的（i，j)元，即A的第i行第j列的元素；
[0008] 貝 1J
[0009] aij=Njji = 1,2,... 256
[0010] 其中A」是標(biāo)準(zhǔn)譜歸一化后第j種元素的的第i道計數(shù)；
[0011] 2)計算實測地層俘獲譜計數(shù)Ci
[0012] 設(shè)在實際測量中得到的256道y能譜的第i道計數(shù)為Ci ;
[0013] 貝丨J
[0014] Cj = NiM i = 1,2,... 256
[0015] 其中N iM是實測地層譜歸一化后的第i道計數(shù)；總道數(shù)i = 256 ;
[0016] 3)計算對角加權(quán)矩陣W
[0017] 設(shè)對角加權(quán)陣為W, w"是W的（i, i)元，令權(quán)系數(shù)r= wn，則

【權(quán)利要求】
1. 一種地層元素測井儀器解譜分析確定產(chǎn)額的方法，其特征在于：包括以下步驟： 1) 計算標(biāo)準(zhǔn)俘獲譜歸一化響應(yīng)矩陣A 響應(yīng)矩陣A是用加權(quán)最小二乘法解譜時的設(shè)計矩陣； 設(shè)％是響應(yīng)矩陣A的（i，j)元，即A的第i行第j列的元素； 則 aij = N^· i = 1，2, · · · 256 其中是標(biāo)準(zhǔn)譜歸一化后第j種元素的的第i道計數(shù)； 2) 計算實測地層俘獲譜計數(shù)Ci 設(shè)在實際測量中得到的256道γ能譜的第i道計數(shù)為Ci ; 則 ci = Nm i = 1，2, ... 256 其中NiM是實測地層譜歸一化后的第i道計數(shù)；總道數(shù)i = 256 ; 3) 計算對角加權(quán)矩陣W 設(shè)對角加權(quán)陣為W，Wu是W的（i，i)元，令權(quán)系數(shù)r = Wii，則 =叼=7 Vi 則 wf, 0 0 0 0 1? 0 0 W- 0 0 ··· 0 〇 0 〇 Wm_ 4) 計算加權(quán)最小二乘法解譜得產(chǎn)額yj 地層測量譜的第i道的計數(shù)可以表示為，地層中所有元素在第i道內(nèi)產(chǎn)生的Y射線即 標(biāo)準(zhǔn)譜的線性疊加； m Cy - ^ /- \ Λ 2 Λ 3....//7 Μ 即疊加的系數(shù)為該地層中元素的產(chǎn)額； 其中ay是由元素的標(biāo)準(zhǔn)譜得到的mXm'的響應(yīng)矩陣Α的（i，j)元，y」表示第j種元 素的百分含量，即產(chǎn)額，ε i表示誤差，m為元素總數(shù)，m'為總道數(shù)； 采用加權(quán)最小二乘法求解，得： y = [ATWA]_1[ATWC] 其中 C = (Cp c2，· · · · cn)，y = (yp y2，· · · ym); 5) 利用殘差平方和的平均數(shù)△2判定擬合效果； 設(shè)殘差平方和的平均數(shù)為 1. m a2 = ?χ，/): ηΤ? μ 若確定的精度Λ2 < ε，則滿足要求，否則返回步驟1)重新進(jìn)行解譜。
【文檔編號】G06F19/00GK104239709SQ201410453114
【公開日】2014年12月24日 申請日期:2014年9月5日 優(yōu)先權(quán)日:2014年9月5日
【發(fā)明者】石文夫, 張政彬, 譚紅建, 王志勇, 李國英, 張崇軍, 于之深, 嵇成高 申請人:西安奧華電子儀器有限責(zé)任公司, 中國石油集團(tuán)渤海鉆探工程有限公司