一种地层元素测井仪器解谱分析确定产额的方法
【专利摘要】本发明提供一种地层元素测井仪器解谱分析确定产额的方法，包括计算标准俘获谱归一化响应矩阵A、计算实测地层俘获谱计数ci、计算对角加权矩阵W、计算加权最小二乘法解谱得产额yj、利用残差平方和的平均数Δ2判定拟合效果等步骤，本发明采用逐道法解谱方法。采用完整的256道能窗解谱。通过M C数值模拟发现该方法对解谱没有太大的影响。利用逐道法对未知地层进行解谱，发现解析的结果误差还较大，但是基本上能反映出地层元素的种类和相对含量。
【专利说明】一种地层元素测井仪器解谱分析确定产额的方法
【技术领域】
[0001] 本发明涉及一种对地层元素测井仪器进行解谱以分析确定产额的方法。
【背景技术】
[0002] 对地层元素测井而言，一般都是采用最小二乘法或者最小二乘法的改进方法进行 解谱。这种解谱方法计算出来的地层元素种类不会出错，但是对于某些含量较少或者是灵 敏度比较低的元素计算的过程中误差较大。
【发明内容】
[0003] 本发明目的是提供一种地层元素测井仪器解谱分析确定产额的方法，其基本上能 反映出地层元素的种类和相对含量。
[0004] 一种地层元素测井仪器解谱分析确定产额的方法，其特殊之处在于：包括以下步 骤：
[0005] 1)计算标准俘获谱归一化响应矩阵A
[0006] 响应矩阵A是用加权最小二乘法解谱时的设计矩阵；
[0007] 设au是响应矩阵A的（i，j)元，即A的第i行第j列的元素；
[0008] 贝 1J
[0009] aij=Njji = 1,2,... 256
[0010] 其中A」是标准谱归一化后第j种元素的的第i道计数；
[0011] 2)计算实测地层俘获谱计数Ci
[0012] 设在实际测量中得到的256道y能谱的第i道计数为Ci ;
[0013] 贝丨J
[0014] Cj = NiM i = 1,2,... 256
[0015] 其中N iM是实测地层谱归一化后的第i道计数；总道数i = 256 ;
[0016] 3)计算对角加权矩阵W
[0017] 设对角加权阵为W, w"是W的（i, i)元，令权系数r= wn，则
【权利要求】
1. 一种地层元素测井仪器解谱分析确定产额的方法，其特征在于：包括以下步骤： 1) 计算标准俘获谱归一化响应矩阵A 响应矩阵A是用加权最小二乘法解谱时的设计矩阵； 设％是响应矩阵A的（i，j)元，即A的第i行第j列的元素； 则 aij = N^· i = 1，2, · · · 256 其中是标准谱归一化后第j种元素的的第i道计数； 2) 计算实测地层俘获谱计数Ci 设在实际测量中得到的256道γ能谱的第i道计数为Ci ; 则 ci = Nm i = 1，2, ... 256 其中NiM是实测地层谱归一化后的第i道计数；总道数i = 256 ; 3) 计算对角加权矩阵W 设对角加权阵为W，Wu是W的（i，i)元，令权系数r = Wii，则 =叼=7 Vi 则 wf, 0 0 0 0 1? 0 0 W- 0 0 ··· 0 〇 0 〇 Wm_ 4) 计算加权最小二乘法解谱得产额yj 地层测量谱的第i道的计数可以表示为，地层中所有元素在第i道内产生的Y射线即 标准谱的线性叠加； m Cy - ^ /- \ Λ 2 Λ 3....//7 Μ 即叠加的系数为该地层中元素的产额； 其中ay是由元素的标准谱得到的mXm'的响应矩阵Α的（i，j)元，y」表示第j种元 素的百分含量，即产额，ε i表示误差，m为元素总数，m'为总道数； 采用加权最小二乘法求解，得： y = [ATWA]_1[ATWC] 其中 C = (Cp c2，· · · · cn)，y = (yp y2，· · · ym); 5) 利用残差平方和的平均数△2判定拟合效果； 设残差平方和的平均数为 1. m a2 = ?χ，/): ηΤ? μ 若确定的精度Λ2 < ε，则满足要求，否则返回步骤1)重新进行解谱。
【文档编号】G06F19/00GK104239709SQ201410453114
【公开日】2014年12月24日 申请日期:2014年9月5日 优先权日:2014年9月5日
【发明者】石文夫, 张政彬, 谭红建, 王志勇, 李国英, 张崇军, 于之深, 嵇成高 申请人:西安奥华电子仪器有限责任公司, 中国石油集团渤海钻探工程有限公司