地层元素测井仪器解谱分析确定产额的方法
【专利摘要】本发明提供地层元素测井仪器解谱分析确定产额的方法,包括确定段区数n及各区段左右边界kn和kn+1、计算标准俘获谱归一化响应矩阵A、计算待解谱在各段区总计数值ci、计算对角加权矩阵W、计算加权最小二乘法解谱得产额yj及利用残差平方和的平均数Δ2判定拟合效果。本发明采用卡段法解谱方法。卡段法:有选择性的选择元素的特征能谱所在的道区解谱。通过MC数值模拟发现该方法对解谱没有太大的影响。
【专利说明】地层元素测井仪器解谱分析确定产额的方法
【技术领域】
[0001] 本发明涉及一种对地层元素测井仪器进行解谱以分析确定产额的方法。
【背景技术】
[0002] 对地层元素测井而言,一般都是采用最小二乘法或者最小二乘法的改进方法进 行解谱。这种解谱方法计算出来的地层元素种类不会出错,但是对于某些含量较少或者是 灵敏度比较低的元素计算的过程中误差较大。
【发明内容】
[0003] 本发明目的是提供一种地层元素测井仪器解谱分析确定产额的方法,其基本上能 反映出地层元素的种类和相对含量。
[0004] 一种地层元素测井仪器解谱分析确定产额的方法,其特殊之处在于:包括以下步 骤:
[0005] 1)确定段区数n及各区段左右边界kn和kn+1 ;所述kn和kn+1分别表示第n个段区 的左右边界;n为不大于256的整数
[0006] 2)计算标准俘获谱归一化响应矩阵A
[0007] 响应矩阵A的计算方法如下;
[0008] 设au是响应矩阵A的(i,j)元,即A的第i行第j列的元素;
[0009]贝1J
[0010]
【权利要求】
1. 一种地层元素测井仪器解谱分析确定产额的方法,其特征在于:包括以下步骤: 1) 确定段区数η及各区段左右边界kn和kn+1 ;所述kn和kn+1分别表示第η个段区的左 右边界;η为不大于256的整数 2) 计算标准俘获谱归一化响应矩阵A 响应矩阵A的计算方法如下; 设%是响应矩阵A的(i,j)元,即A的第i行第j列的元素; 则 (1 ?=Σ? ; kn和kn+1分别表示第η个段区的左右边界,Nw是准谱归一化之后的第k道的计数; 3) 计算待解谱在各段区总计数值Ci 设在实际测量中得到的256道γ能谱的第i道计数为Ci ; 则 Cfh n = \, 2,... η' k-K 其中,匕和、+1分别表示第n个段区的左右边界,NkM是实测地层谱归一化后的第k道 计数; 4) 计算对角加权矩阵W 设对角加权阵为W,Wu是W的(i,i)元,令权系数r = Wii,则 1 wu = W) = - c; 则: wjj Ο Ο Ο Ο n Ο Ο W= Ο Ο ··· Ο Ο ο ο Wnn_ 5) 计算加权最小二乘法解谱得产额yj 地层测量谱的第i道的计数可以表示为,地层中所有元素在第i道内产生的Y射线即 标准谱的线性叠加; 1? c- = a^yj + εi /=1、2、3......m >ι 即叠加的系数为该地层中元素的产额; 其中是由元素的标准谱得到的mXm,的响应矩阵A的(i,j)元,y」表示第j种元 素的百分含量,即产额,表示误差,m为元素总数,m,为总道数,采用加权最小二乘法求 解,得: y = [ATWA]_1[ATWC] 其中 C = (c" c2, · · · · cn),y = (y" y2, · · · ym); 6)利用残差平方和的平均数Λ2判定拟合效果; 设残差平方和的平均数为 Δ =^£(c,-少 11 i=\ j=l 若确定的精度Λ2〈 ε,则满足要求,否则返回步骤1)重新进行解谱。
【文档编号】G01N33/24GK104239710SQ201410453182
【公开日】2014年12月24日 申请日期:2014年9月5日 优先权日:2014年9月5日
【发明者】石文夫, 张政彬, 谭红建 申请人:西安奥华电子仪器有限责任公司