[發明專利]一種用于地層元素測井的元素產額計算方法

有效

• 申請號: 201811093260.2
• 申請日: 2018-09-19
• 公開（公告）號: CN109241644B
• 公開（公告）日: 2023-05-09
• 發明（設計）人: 劉軍濤;張鋒;劉樹才;苗彬
• 申請（專利權）人: 蘭州大學
• 主分類號: G06F17/18
• 分類號: G06F17/18;G06F17/16
• 代理公司: 青島潤集專利代理事務所(普通合伙) 37327
• 代理人: 趙以芳

1.一種應用于地層元素測井的元素產額計算方法，其特征在于，依次包括以下步驟：

(a)自動選取第一種元素種類組合，若地層最多含有9種元素，且最少含有7種元素，則元素種類的可能組合數一共有C(9,9)+C(9,8)+C(9,7)＝46種；若已知地層必定含有某一種或幾種元素，?則相應的元素種類組合數會相對減少；

(b)選取一種元素種類組合之后，利用奇異值分解方法計算元素產額初值y⁽⁰⁾＝VS^-1UC，并令k＝0，若某種元素的產額為負值，將其賦值為0.001，式中：A＝U·S·V^T，A為所選元素種類標準譜組成的m×n系數矩陣，其中m為選取的組合中的元素種類數，n為測量伽馬能譜的能道數，U，S，V為矩陣A的奇異值分解矩陣，C為測量得到的n×1伽馬能譜矩陣；

(c)令k＝k+1，然后計算：

；

式中：a(i,j)為矩陣A中的矩陣元素，λ為約束項系數；y^(k+1)(i)代表第k+1次迭代得到的第i種元素產額；y^(k)(i)代表第k次迭代得到的第i種元素產額；

(d)如果|y^(k+1)-y^(k)|＜T或k＞k^T，終止迭代；否則，重復步驟(c)；

式中：T與k^T分別為殘差及迭代次數截止值，然后計算AIC＝-20(y)+2·n及BIC＝-20(y)+ln(m)n；并計算兩者比值AIC/BIC，

式中：0(y)為擬合能譜與測量能譜的殘差；

(e)重復步驟(a)-(d)，直到所有元素種類組合的情況計算完畢；

(f)對比不同元素種類組合條件下，AIC與BIC的比值；當AIC/BIC具有最小值時，計算結果即為地層元素種類與元素產額。