Matlab稳健回归法在FD-3019γ测井仪校准中的应用.pdf

《Matlab稳健回归法在FD-3019γ测井仪校准中的应用.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第32卷第4期铀矿地质Vo1．32NO．42016定7月UraniumGeologyJuly2016DOI：10．3969／j．issn．1000—0658．2016．04．008Matlab稳健回归法在FD一3019测井仪校准中的应用何涛(核工业280研究所，四川广汉618300)[摘要]简要介绍了FD一3019．y测井仪的校准方法，即Matlab稳健同归法及其抗差性。通过一个实例对比了Matlab稳健回归法和最小二乘法建立的一元线性回归方程，结果表明Matlab稳健回归法较之最小二乘法更为准确。[关键词]校准；FD

2、一3019)'测井仪；稳健回归；Matlab[文章编号]1000—0658(2016)04—0242一O4[中图分类号]P631．6[文献标志码]A在7测井中，每个月均要用固体标准镭方程。(1)求7照射量率换算系数：源对FD一3019y测井仪进行校准检查，以确I=：=(N～N)×k(1)定其7照射量率换算系数是否符合要求，通式中：1一距固体镭源为R处的7照射量率常用一元线性回归(最小二乘法)计算。用(nC·kg·h)；N一距固体镭源为R时最小二乘法进行回归分析时，数据偏离越大FD一3019仪的计数率(cps)；N。一F

3、D一3Ol9对回归方程的影响就越大。偏离较大的异常仪的底数(cps)；k一了照射量率换算系数数据比较容易剔除，而偏离达不到异常的数(nC·kg～·h一·cps)。据就不容易剔除，致使最小二乘法建立的回实际应用中，把式(1)简化为如下归方程的准确性降低，表明最小二乘法的抗方程：差性较差。稳健回归法采用M一估计作为目标I—A+k×N(2)函数，无需判别异常值，全部数据参与计算，A一～N×k(3)是目前应用最广泛的稳健估计方法，具有较式中：A一底数(nC·kg·h)。强的抗差性n]。用稳健回归法建立FD一3019依据不同的R

4、值，可以计算出相应的I仪校准回归方程能有效降低数据观测误差的值。以FD一3019仪器测量1O次的平均值作影响。为距固体镭源R处的N值。移除镭源，测量3O次的平均值作为测量底数N。值。在FD一1FD一3019"／测井仪校准方法3019仪器测程范围内均匀地给出不少于10个用固体标准镭源校准时，可以建立直线测量值的点，每个测量值的点上测量次数应[收稿日期]2015～O3一O9[改回日期]2016—03—24[作者简介]何涛(1965一)，男，高级工程师，1987年毕业于成都地质学院三系，主要从事地球物理勘探和数据处理工作。E

5、—mail：164O214830@qq．com第4期何涛：Matlab稳健回归法在FD一30197测井仪校准中的应用。243·不少于10个。2．2抗差性用测量值N作为式(2)中的N值进行回以y===1O一2x为例，X取1到1O的整数，归分析计算出A和k，利用式(3)计算出仪给Y增加随机误差，同时令Y中某一数据产生器的底数，记为N。为了便于区分和说明，较大偏离。执行如下代码：把用水面法或其他方法实测的仪器底数记为x一(1：10)’；形成数据XNo。那么，N与No的误差越小越好。y一1O一2*x+randn(10，1)；形

6、成数用(N—N)代替式(2)中的N进行回据Y归分析时，由于减去了仪器的底数，相当于y(10)一0；形成偏离数据式(3)中的N。等于0，此时计算出的A值应bls=regress(y，Eones(10，1)x])；最小为0。但实际情况并非如此，由于观测误差的二乘法回归影响计算出的A值不等于0，根据经验，此时brob—robustfit(x，Y，’logistic’，0．3)；％稳的A应当小于0．3，才能使计算出的铀含量误健回归差小于10×10～。scatter(x，Y，’filled’，’k’)；gridon；holdon

7、；绘制散点图2Matlab稳健回归y一10—2*x；形成直线Y一10—2x的在Matlab的统计工具箱中稳健回归函数Y数据为robustfit，通过查阅文件，其函数使用方plot(x，y，’k一，’LineWidth，2)；绘制y一法及作用如下。1O一2*X直线2．1函数说明plot(x，bls(1)+bls(2)*x，’k一’，常用的调用方式为：’Linewidth’，2)；％绘制最小二乘法回归直线d—robustfit(X，y，wfun，tune)plot(X，brob(1)+brob(2)X，’k一．’，其中：x

8、一列向量，直线方程Y—a+bx中的’LineWidth’，2)；绘制稳健回归直线自变量x；y一列向量，直线方程y—a+bx中legend(。数据点’，’最小二乘法’，’稳健回的因变量Y；d一返回值，向量，d(0)、归’)；％绘制图例d(1)分别表示直线方程y=a+bx的系数a和xlabel(’X’，’fontsize’，12，’f