一种用于地层元素测井的元素产额计算方法技术

本发明专利技术公开了一种应用于地层元素测井的元素产额计算方法，发明专利技术的元素产额计算方法主要包括以下几个步骤：(1)首先自动选定一种元素种类组合。(2)然后利用奇异值分解方法计算对应元素产额初值y(0)，当计算得到的某种元素产额为负值时，赋值为0.001。(3)令k＝0，利用极大似然的约束迭代方法计算得到元素产额y(k+1)，令k＝k+1。(4)如果满足迭代截止条件，终止循环，然后计算AIC与BIC，及两者比值AIC/BIC。(5)重复步骤1‑3，直到所有元素种类组合计算完毕。(6)选择对应AIC/BIC最小值的元素产额为最终元素产额。该发明专利技术方法与常用的加权最小二乘元素产额计算方法相比较，降低了噪音干扰及和测量自身统计波动对元素产额计算准确度的影响，提高了元素产额计算精度；并且能够实现元素种类的自动确定。

【技术实现步骤摘要】
一种用于地层元素测井的元素产额计算方法
本专利技术涉及石油天然气开发领域，具体涉及一种针对地层元素测井的元素产额计算的方法。
技术介绍
地层元素测井是指通过利用同位素中子源或是可控脉冲中子源，利用中子源发射的中子与测井仪器周围物质作用产生伽马射线；处理伽马探测器测量得到的伽马能谱数据，评价地层骨架元素组成的测井方法。元素测井仪器通常采用一个伽马晶体探测器，在井眼内进行测量，基于建立的获取的元素标准能谱数据，通过伽马能谱解谱，获取元素产额，进而利用氧化物闭合模型将元素产额转换为元素的干重含量。目前针对地层元素测井的元素产额计算方法主要是加权最小二乘方法。放射性衰变是一个随机过程，因此针对放射性衰变的伽马射线的测量会在一定程度上受到统计波动性的影响。这些内在的波动不可避免地对测量伽马能谱产生影响，造成相同测量条件及测量时间的情况下，测量得到的伽马能谱的能道计数不完全相同。地层元素测井测量伽马能谱的这种自身统计特性，会对元素产额的计算产生影响。另外由于电路信号、自然伽马本底等因素产生的能谱噪音也会影响元素产额的计算准确度。除此之外，目前常用的加权最小二乘元素产额计算方法，需要人为确定元素的本文档来自技高网...

【技术保护点】
1.一种应用于地层元素测井的元素产额计算方法，其特征在于，依次包括以下步骤：(a)自动选取第一种元素种类组合。例如若地层最多含有9种元素，且最少含有7种元素，则元素种类的可能组合数一共有C(9,9)+C(9,8)+C(9,7)＝46种；若已知地层必定含有某一种或几种元素，则相应的元素种类组合数会相对减少。(b)选取一种元素种类组合之后，利用奇异值分解方法计算元素产额初值y(0)＝VS‑1UC，并令k＝0，若某种元素的产额为负值，将其赋值为0.001。式中：A＝U·S·VT，A为所选元素种类标准谱组成的m×n系数矩阵，其中m为选取的组合中的元素种类数，n为测量伽马能谱的能道数，U，S，V为矩阵A...

【技术特征摘要】
1.一种应用于地层元素测井的元素产额计算方法，其特征在于，依次包括以下步骤：(a)自动选取第一种元素种类组合。例如若地层最多含有9种元素，且最少含有7种元素，则元素种类的可能组合数一共有C(9,9)+C(9,8)+C(9,7)＝46种；若已知地层必定含有某一种或几种元素，则相应的元素种类组合数会相对减少。(b)选取一种元素种类组合之后，利用奇异值分解方法计算元素产额初值y(0)＝VS-1UC，并令k＝0，若某种元素的产额为负值，将其赋值为0.001。式中：A＝U·S·VT，A为所选元素种类标准谱组成的m×n系数矩阵，其中m为选取的组合中的元素种类数，n为测量伽马能谱的能道数，U，S，V为矩阵A的奇异值分解矩阵，C为测量得到的n×1伽马能谱矩...

【专利技术属性】
技术研发人员：刘军涛，张锋，刘树才，苗彬，
申请(专利权)人：中国矿业大学，
类型：发明
国别省市：江苏,32