[发明专利]一种用于地层元素测井的元素产额计算方法

说明书

本发明公开了一种应用于地层元素测井的元素产额计算方法，发明的元素产额计算方法主要包括以下几个步骤：(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最小值的元素产额为最终元素产额。该发明方法与常用的加权最小二乘元素产额计算方法相比较，降低了噪音干扰及和测量自身统计波动对元素产额计算准确度的影响，提高了元素产额计算精度；并且能够实现元素种类的自动确定。

技术领域

本发明涉及石油天然气开发领域，具体涉及一种针对地层元素测井的元素产额计算的方法。

背景技术

地层元素测井是指通过利用同位素中子源或是可控脉冲中子源，利用中子源发射的中子与测井仪器周围物质作用产生伽马射线；处理伽马探测器测量得到的伽马能谱数据，评价地层骨架元素组成的测井方法。元素测井仪器通常采用一个伽马晶体探测器，在井眼内进行测量，基于建立的获取的元素标准能谱数据，通过伽马能谱解谱，获取元素产额，进而利用氧化物闭合模型将元素产额转换为元素的干重含量。目前针对地层元素测井的元素产额计算方法主要是加权最小二乘方法。

放射性衰变是一个随机过程，因此针对放射性衰变的伽马射线的测量会在一定程度上受到统计波动性的影响。这些内在的波动不可避免地对测量伽马能谱产生影响，造成相同测量条件及测量时间的情况下，测量得到的伽马能谱的能道计数不完全相同。地层元素测井测量伽马能谱的这种自身统计特性，会对元素产额的计算产生影响。另外由于电路信号、自然伽马本底等因素产生的能谱噪音也会影响元素产额的计算准确度。除此之外，目前常用的加权最小二乘元素产额计算方法，需要人为确定元素的种类组合，且受测量能谱统计波动及噪音干扰影响较大。

基于上述地层元素测井元素产额计算中存在的问题，本发明提出一种应用于地层元素测井的元素产额计算，降低能谱噪音及自身统计波动性对元素产额计算的影响，提高元素产额计算的准确度。

发明内容

针对上述元素测井元素产额计算存在的问题，本发明提出了一种元素产额计算方法，降低噪音干扰与自身统计性对元素产额计算的影响，实现元素种类与产额的自动识别与准确计算。

其技术解决方案包括：

一种地层元素测井元素产额的计算方法，依次包括以下步骤：

(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中的矩阵元素，λ为约束项系数。

(d)如果|y^(k+1)-y^(k)|＜T或k＞k_T，终止迭代；否则，重复步骤(c)。