CN108457644B - 一种元素俘获能谱测井的伽马能谱解谱方法及装置 - Google Patents

Abstract

本发明提供了一种元素俘获能谱测井的伽马能谱解谱方法及装置。该解谱方法首先对元素俘获能谱测井仪获取的数据进行预处理；然后通过对不同元素谱形、特征峰道及本底之间的理论分析和数值计算，按照元素间相互影响程度构建主元素群和辅助元素群；进而在主、辅助元素群构建的基础上利用最小二乘法进行解谱；最后根据最小二乘法解谱获得的各元素产额通过理论进行能谱重构，将测量的伽马能谱与重构伽马能谱对比进行误差控制，提高解谱精度。该解谱方法具有高精度、可操作性强等优点，对不同类型储层均具有较好的普适性，同时为其他类型仪器的中子伽马能谱数据的处理提供了重要的参考和借鉴意义。

Description

一种元素俘获能谱测井的伽马能谱解谱方法及装置

技术领域

本发明涉及一种元素俘获能谱测井的伽马能谱解谱方法及装置，属于油气勘探领域中的测井技术领域。

背景技术

元素俘获能谱(Elemental Capture Spectroscopy，简称ECS)测井仪是由斯伦贝谢公司向中国市场推出的新型地层元素测井仪器，该仪器利用中子与地层原子核碰撞发生核反应的原理得到伽马能谱，该伽马能谱包含了地层中主要造岩元素Si、Ca、Fe、Al、S、Ti、H、Gd等元素的信息，为复杂储层岩性精细评价提供了一种新途径。元素俘获能谱测井解释的核心是利用各元素标准谱对地层原始测量谱进行刻度，通过解谱获得各种元素的产额。在元素俘获能谱仪器获取的测量谱中，既包含了俘获伽马信息又包含了非弹性散射信息，而对油田现场解释评价最重要的是俘获伽马信息，因此在对测量谱进行解谱处理时必须考虑非弹性散射信息对解谱精度的影响。与此同时，不同地层元素之间的相互影响也制约着中子俘获伽马能谱的解谱精度。

以往学者虽然对元素俘获能谱测井的解谱方法开展了研究，但主要侧重于数值模拟和理论阐述上，没有针对实际井下数据提供有效的技术方法，因此未能在实际生产中进行应用。

发明内容

为解决上述技术问题，本发明的目的在于一种元素俘获能谱测井的伽马能谱解谱方法及装置。该方法是针对元素俘获能谱测井仪器得到的原始测量谱进行解谱的方法，具有高精度、可操作性强等优点。

为了实现上述目的，本发明提供了一种元素俘获能谱测井的伽马能谱解谱方法，其包括以下步骤：

步骤1：获取、分析研究区资料，所述资料至少包括元素俘获伽马能谱测井数据(在本发明中也称为测量伽马能谱、测量谱、测量混合伽马能谱、总元素伽马能谱测量谱等)和全岩氧化物分析数据，确定该区块主要元素类型；

步骤2：对元素俘获伽马能谱测井数据进行预处理，所述预处理包括能窗的选择、能谱平滑滤波、归一化处理以及非弹性散射信息扣除，得到预处理后的元素俘获伽马能谱；

步骤3：构建由不同元素组成的主元素群和辅助元素群，以明确不同元素的解谱先后顺序；

步骤4：根据步骤3得到的主元素群和辅助元素群，在经步骤2预处理后的元素俘获伽马能谱中，利用最小二乘法先对主元素进行解谱，然后从预处理后的总元素伽马能谱(即经步骤2预处理后的元素俘获伽马能谱)中扣除全部主元素所占贡献，再利用最小二乘法对辅助元素进行解谱，从而得到各元素的相对产额；

步骤5：根据步骤4得到的各元素的相对产额，结合归一化处理后的单元素的元素俘获标准伽马能谱(在本发明中也称为单元素标准谱)，进行元素俘获伽马能谱的重构，并将重构的伽马能谱与测量伽马能谱进行对比，确定解谱结果是否可靠。

根据本发明的具体实施方式，优选地，在上述解谱方法的步骤1中，所获取、分析的研究区资料还包括常规测井数据、录井数据和地质资料等中的一种或几种的组合。

根据本发明的具体实施方式，优选地，在上述解谱方法的步骤1中，研究区主要元素类型是通过以下方式确定的：通过分析研究区的全岩氧化物分析数据、以及可选择地，常规测井数据、录井数据和地质资料中的一种或几种的组合，得出不同种类氧化物所占岩石的重量百分比，依据各氧化物重量百分比确定该区块主要元素类型。其中，所述地质资料主要包括：岩心薄片和/或地质报告等。全岩氧化物分析数据是确定研究区主要元素类型所必需的资料，同时可以参考常规测井数据、录井数据和地质资料中的一种或几种资料。

在本发明的解谱方法中，步骤2对伽马能谱进行预处理，能够消除仪器在测量过程中受井眼环境、仪器固有特性等因素的影响，提高解谱的精度。

在能窗的选择方面，元素俘获能谱测井探测到的伽马射线能量一般集中在0-10MeV，对应道址为0-255，总共256道。低能区段容易受到康普顿平台效应以及仪器的噪声等因素影响而产生很高的计数，即大量干扰信号，应在解谱之前予以剔除。而在高能区段，能谱的计数往往非常低，在解谱的过程中对最终的求解结果影响不大，也应在解谱之前予以剔除。如果没有完全去掉带有不利影响的数据，会对解谱产生影响；但是，如果去掉过多的数据会去掉有用信息，也会对解谱产生影响。因此，选取适当的能量范围是解谱前必须解决的关键问题。

根据本发明的具体实施方式，优选地，在上述解谱方法的步骤2中，所选择的能窗范围是30-210道。本发明提出的选取解谱能窗的基本原则是不去掉元素全能峰所对应的能量。本发明的发明人通过对国际原子能机构公布的核数据进行查询，查出解谱涉及的各种元素的特征峰位，根据峰位所对应的能量信息，确定解谱的能窗范围。各元素特征峰位按能量从低到高排序为1.38，1.808，1.924，1.951，2.073，2.092，2.223，2.282，2.379，2.828，3.033，3.22，3.539，3.587，3.691，3.981，4.419，4.733，4.737，4.869，4.933，5.42，5.9，5.92，6.018，6.11，6.36，6.379，6.395，6.418，6.42，6.76，7.414，7.631，7.646，7.724，7.769，7.79，8.153(MeV)。根据这些数据将能窗的能量范围取为1.2MeV～8.4MeV，其对应的能窗范围是30-210道，如图1所示。

根据本发明的具体实施方式，优选地，在上述解谱方法的步骤2中，能谱平滑滤波的方式是采用Savitzky-Golay滤波器，并用Savitzky-Golay五点滤波法，对所选择的能窗范围内的元素俘获伽马能谱测井数据进行滤波，Savitzky-Golay五点滤波法的公式如下式1所示：

其中，

表示某一道址上滤波后的计数，y_i表示该道址上的计数，y_i-1表示该道址的前一道址上的计数，y_i-2表示该道址的前二道址上的计数，y_i+1表示该道址的后一道址上的计数，y_i+2表示该道址的后二道址上的计数。

元素俘获伽马能谱测井仪器的测量过程中，存在固有的统计涨落，能谱平滑滤波能够消除测量伽马能谱中存在的假峰、弱峰等问题。本发明采用的Savitzky-Golay滤波器，在滤波的窗长内，利用多项式拟合测量伽马能谱的数据，得到滤波后的结果，即拟合多项式在滤波点对应的值。通过对原始测量伽马能谱的分析，本发明提出的Savitzky-Golay五点滤波法可以消除毛刺。本发明的滤波方法的实质是对测量数据做加权平均，使得能谱数据中每个数据波动都做了平均，这种平均对数据趋势影响非常小，而波动的平均值为零，可有效消除统计涨落的影响。

根据本发明的具体实施方式，优选地，在上述解谱方法的步骤2中，归一化处理采用的方法是：对经过能窗范围选择、能谱平滑滤波之后的元素俘获伽马能谱测井数据，使181道(逐道)的能谱数据之和为10，得到归一化处理后的元素俘获伽马能谱，所采用的公式如下式2所示：

其中，N_Gkj为j深度点对应的归一化处理后的元素俘获伽马能谱在第k道的计数，N_kj为j深度点对应的经过能窗范围选择、能谱平滑滤波之后的元素俘获伽马能谱在第k道的计数，j为地层下的深度点。

在本发明的具体实施方式中，归一化处理是将整个能量段的能谱计数累加，再将每道计数除以能窗范围内所有道址对应的计数之和，从而使归一化后的数据在0～10之间，整个能量段的能谱数据之和为10。在本发明中，归一化处理是对滤波后的能窗范围内的伽马能谱数据进行归一化，获取归一化后的伽马能谱数据，使归一化处理后的俘获元素定性计算值可以准确反映整个井段的变化，因此归一化处理可以为元素定性计算和分离的方法提供保证，也可以准确地反应能窗范围内的地层变化。

根据本发明的具体实施方式，优选地，在上述解谱方法的步骤2中，非弹性散射信息扣除的方法是：对归一化处理后的元素俘获伽马能谱选取30-54道、55-75道和76-210道三段道址的计数分别用0.9、0.7和0.8作为扣除系数进行非弹性信息扣除。由于元素俘获能谱测井获取的伽马能谱中既包含了俘获伽马信息又包含了非弹性散射信息，而对油田现场解释评价最重要的是俘获伽马信息，因此在对测量谱进行剥谱处理时必须扣除非弹性散射信息对总元素伽马能谱测量谱的贡献，可以有效提高后期解谱的精度。在本发明的具体实施方式中，对测量谱中非弹性散射信息扣除的方法采用的是分段扣除方法，根据不同元素对不同能量段总计数的贡献，在不同能量段分别采用不同的扣除系数消除测量谱中非弹性散射信息的影响。

根据本发明的具体实施方式，优选地，在上述解谱方法的步骤3中，所构建的主元素群包括：Si、Ca、S、H、Cl、Ti、Fe、Na、Ba和Gd等元素中的一种或几种元素，这些元素的贡献占总元素伽马能谱的80％以上；所构建的辅助元素群包括：Mg、K、Cr、Ni、I、Tb和Al等元素中的一种或几种元素。单元素标准伽马能谱在谱形、特征峰位的差异性是元素俘获能谱测井解谱的核心和关键，但实际上，不同元素之间在谱形及特征峰位分布上会存在一定的相似性和重叠。如图2所示，Na元素和Mg元素在谱形上存在一定的相似性，在特征能峰所在道址上存在重叠，因此，在实际解谱过程中如果不考虑元素之间的相互影响因素，则会降低解谱的精度，从而影响油气储层评价的准确性。本发明提出和构建的元素群是在考虑不同元素单元素标准伽马能谱相互影响的基础上建立的；主元素群主要包括构成岩石骨架部分的元素；辅助元素群的构建是为了消除或者辅助提高其它元素的解谱精度而建立的，辅助元素Mg和K主要用于分析判断S的产额，辅助元素Mg、Cr、Ni用于分析判断Si、Ca和Fe元素的产额，辅助元素Cr、Ni、I、Tb和Al用于分析判断H、Gd和Ti素的产额，进而提高解谱的精度。

根据本发明的具体实施方式，优选地，在上述解谱方法的步骤4中，利用最小二乘法进行解谱是按照如下式3所示的算法进行的，以得到各元素的相对产额：

其中，y_j为第j种元素的相对产额(即该元素对总计数的贡献)，a_ij为第j种元素的归一化处理后的单元素的元素俘获标准伽马能谱在第i道的计数，c_i为预处理后的元素俘获伽马能谱在第i道的计数，ε_i为校正系数，ε_i≤0.1。

本发明利用上述最小二乘法公式在对主元素群进行解谱后，从预处理后的总元素伽马能谱中对主元素所占贡献进行扣除，再利用上述最小二乘法公式对辅助元素群进行再次解谱，逐步得到主、辅元素群中的各元素产额信息，从而有效提高解谱精度。

根据本发明的具体实施方式，优选地，在上述解谱方法的步骤5中，所述归一化处理后的单元素的元素俘获标准伽马能谱是通过以下方法得到的：对于单元素的元素俘获标准伽马能谱，使256道(逐道)的能谱数据之和为10，所采用的公式如下式4所示：

其中，N_Gkj为j深度点对应的归一化处理后的单元素的元素俘获标准伽马能谱在第k道的计数，N_kj为j深度点对应的归一化处理前的单元素的元素俘获标准伽马能谱在第k道的计数，j为地层下的深度点。

根据本发明的具体实施方式，优选地，在上述解谱方法的步骤5中，元素俘获伽马能谱的重构是按照下式5计算重构伽马能谱在第i道的计数(即所有解谱元素在第i道的计数)，得到各道的计数后，绘制重构伽马能谱：

其中，X_i为重构伽马能谱在第i道的计数，y_j为第j种元素的相对产额(由式3计算得到)，a_ij为第j种元素的归一化处理后的单元素的元素俘获标准伽马能谱在第i道的计数(由式4计算得到)，ε_i为校正系数，ε_i≤0.1。

根据本发明的具体实施方式，优选地，在上述解谱方法的步骤5中，将重构的伽马能谱与测量伽马能谱进行对比是按照式6进行：

︱C_i－X_i︱《ε (式6)

c_i为测量伽马能谱在第i道的计数，X_i为重构的伽马能谱在第i道的计数，ε为相对误差；

当相对误差小于或等于5％时，认为解谱结果是可靠的；当相对误差大于5％，则认为解谱结果不可靠。

在本发明的具体实施方式中，伽马能谱的重构方式是参考剥谱分析方法，认为测量混合伽马能谱是单元素的线性组合，在确定地层元素含量时以获取的单元素标准谱数据库为基础，根据已知俘获和非弹性散射元素产额结果进行重构，然后将实际测量的伽马能谱与重构伽马能谱进行对比，确定解谱结果可靠与否，进而控制解谱精度。

本发明还提供了一种元素俘获能谱测井的伽马能谱解谱装置，其包括：

数据采集及主要元素类型确定模块，用于获取、分析研究区资料，所述资料至少包括元素俘获伽马能谱测井数据和全岩氧化物分析数据，确定该区块主要元素类型；

元素俘获伽马能谱测井数据预处理模块，用于对元素俘获伽马能谱测井数据进行预处理，所述预处理包括能窗的选择、能谱平滑滤波、归一化处理以及非弹性散射信息扣除，得到预处理后的元素俘获伽马能谱；

主元素群和辅助元素群构建模块，用于构建由不同元素组成的主元素群和辅助元素群，以明确不同元素的解谱先后顺序；

基于元素群的最小二乘法解谱模块，用于根据主元素群和辅助元素群，在预处理后的元素俘获伽马能谱中，利用最小二乘法先对主元素进行解谱，然后从预处理后的总元素伽马能谱(即预处理后的元素俘获伽马能谱)中扣除全部主元素所占贡献，再利用最小二乘法对辅助元素进行解谱，从而得到各元素的相对产额；

伽马能谱重构及误差控制模块，用于根据各元素的相对产额，结合归一化处理后的单元素的元素俘获标准伽马能谱，进行元素俘获伽马能谱的重构，并将重构的伽马能谱与测量伽马能谱进行对比，确定解谱结果是否可靠。

在上述的装置中，优选地，在数据采集及主要元素类型确定模块中，所获取、分析的研究区资料还包括常规测井数据、录井数据和地质资料等中的一种或几种的组合。

在上述的装置中，优选地，在数据采集及主要元素类型确定模块中，研究区主要元素类型是通过以下方式确定的：通过分析研究区的全岩氧化物分析数据、以及可选择地，常规测井数据、录井数据和地质资料等中的一种或几种的组合，得出不同种类氧化物所占岩石的重量百分比，依据各氧化物重量百分比从而确定该区块主要元素类型。

在上述的装置中，优选地，在元素俘获伽马能谱测井数据预处理模块中，所选择的能窗范围是30-210道。

在上述的装置中，优选地，在元素俘获伽马能谱测井数据预处理模块中，能谱平滑滤波的方式是采用Savitzky-Golay滤波器，并用Savitzky-Golay五点滤波法，对所选择的能窗范围内的元素俘获伽马能谱测井数据进行滤波，Savitzky-Golay五点滤波法的公式如下式1所示：

其中，

表示某一道址上滤波后的计数，y_i表示该道址上的计数，y_i-1表示该道址的前一道址上的计数，y_i-2表示该道址的前二道址上的计数，y_i+1表示该道址的后一道址上的计数，y_i+2表示该道址的后二道址上的计数。

在上述的装置中，优选地，在元素俘获伽马能谱测井数据预处理模块中，归一化处理采用的方法是：对经过能窗范围选择、能谱平滑滤波之后的元素俘获伽马能谱测井数据，使181道(逐道)的能谱数据之和为10，得到归一化处理后的元素俘获伽马能谱，所采用的公式如下式2所示：

其中，N_Gkj为j深度点对应的归一化处理后的元素俘获伽马能谱在第k道的计数，N_kj为j深度点对应的经过能窗范围选择、能谱平滑滤波之后的元素俘获伽马能谱在第k道的计数，j为地层下的深度点。

在上述的装置中，优选地，在元素俘获伽马能谱测井数据预处理模块中，非弹性散射信息扣除的方法是：对归一化处理后的元素俘获伽马能谱选取30-54道、55-75道和76-210道三段道址的计数分别用0.9、0.7和0.8作为扣除系数进行非弹性信息扣除

在上述的装置中，优选地，在主元素群和辅助元素群构建模块中，所构建的主元素群包括：Si、Ca、S、H、Cl、Ti、Fe、Na、Ba和Gd等元素中的一种或几种元素；所构建的辅助元素群包括：Mg、K、Cr、Ni、I、Tb和Al等元素中的一种或几种元素。

在上述的装置中，优选地，在基于元素群的最小二乘法解谱模块中，利用最小二乘法进行解谱是按照如下式3所示的算法进行的，以得到各元素的相对产额：

其中，y_j为第j种元素的相对产额(即该元素对总计数的贡献)，a_ij为第j种元素的归一化处理后的单元素的元素俘获标准伽马能谱在第i道的计数，c_i为预处理后的元素俘获伽马能谱在第i道的计数，ε_i为校正系数，ε_i≤0.1。

在上述的装置中，优选地，在伽马能谱重构及误差控制模块中，所述归一化处理后的单元素的元素俘获标准伽马能谱是通过以下方法得到的：对单元素的元素俘获标准伽马能谱，使256道(逐道)的能谱数据之和为10，所采用的公式如下式4所示：

其中，N_Gkj为j深度点对应的归一化处理后的单元素的元素俘获标准伽马能谱在第k道的计数，N_kj为j深度点对应的归一化处理前的单元素的元素俘获标准伽马能谱在第k道的计数，j为地层下的深度点。

在上述的装置中，优选地，在伽马能谱重构及误差控制模块中，元素俘获伽马能谱的重构是按照下式5计算重构伽马能谱在第i道的计数(即所有解谱元素在第i道的计数)，得到各道的计数后，绘制重构伽马能谱：

其中，X_i为重构伽马能谱在第i道的计数，y_j为第j种元素的相对产额(由式3计算得到)，a_ij为第j种元素的归一化处理后的单元素的元素俘获标准伽马能谱在第i道的计数(由式4计算得到)，ε_i为校正系数，ε_i≤0.1。

在上述的装置中，优选地，在伽马能谱重构及误差控制模块中，将重构的伽马能谱与测量伽马能谱进行对比是按照式6进行：

︱C_i－X_i︱《ε (式6)

c_i为测量伽马能谱在第i道的计数，X_i为重构的伽马能谱在第i道的计数，ε为相对误差；

当相对误差小于或等于5％时，认为解谱结果是可靠的；当相对误差大于5％，则认为解谱结果不可靠。

本发明对元素俘获伽马能谱测井资料的处理方法进行了深入研究之后，提供了一种高精度的、具有可操作性的针对元素俘获能谱测井仪器得到的伽马能谱的解谱方法及装置，为准确开展元素俘获能谱的解谱提供了优秀的技术手段。本发明首先对元素俘获能谱测井仪器获取的伽马能谱数据进行预处理，消除了仪器在测量过程中受井眼环境、仪器固有特性等因素影响。然后首次在考虑不同元素间相互影响的基础上提出了基于元素群的中子俘获伽马能谱解谱思想和具体实现方法；通过对不同元素谱形、特征峰道及本底之间的理论分析和数值计算，按照元素间相互影响程度构建了主元素群和辅助元素群，并首次明确了不同元素解谱的先后顺序。之后首次基于最小二乘法实现了基于元素群的中子俘获伽马能谱解谱的方法，并根据最小二乘法解谱获得的各元素产额通过理论进行能谱重构，将实际测量的伽马能谱与重构伽马能谱对比进行误差控制，有效提高了解谱精度。油田现场应用表明，本发明的解谱方法及装置对不同类型储层均具有较好的普适性，同时为其他类型仪器的中子伽马能谱数据处理提供了重要的参考和借鉴意义。

附图说明

图1为解谱能窗的选择示意图。

图2为Na元素和Mg元素单元素的元素俘获标准伽马能谱对比图。

图3为DAS3井2403.1956m的各元素解谱效果。

图4为DAS3井2403.1956m的重构伽马能谱和测量伽马能谱对比图。

图5为DAS3井2403.1956m的预处理后的元素俘获伽马能谱和测量伽马能谱对比图。

图6为DAS3井2403.1956m的主元素群重构伽马能谱和测量伽马能谱对比图。

具体实施方式

为了对本发明的技术特征、目的和有益效果有更加清楚的理解，结合油田现场实际应用实例，现对本发明的技术方案进行以下详细说明，但不能理解为对本发明的可实施范围的限定。

实施例1

本实施例提供了一种元素俘获能谱测井的伽马能谱解谱方法，其可以包括以下步骤：

1.总元素伽马能谱测量谱的获取以及研究区主要元素分布种类的确定

利用元素俘获能谱测井仪器对对大庆油田DAS3井储层段3230-3270m进行测井，获得各深度点的地层元素总的伽马能谱，即元素俘获伽马能谱测井数据，每个深度点对应一张谱图。获取、分析大庆油田DAS3井的全岩氧化物分析实验结果、录井数据，发现在储层段3230m-3270m主要包括Si、Ca、Al、Fe、S、K、Na、Mg、Ti、Gd、H和Ba元素，这些元素的氧化物占岩石重量百分含量的98.9％。

2.对元素俘获伽马能谱数据进行预处理

对DAS3井储层段3230-3270m各深度点的测量谱分别进行预处理，分别得到各深度点对应的预处理后的元素俘获伽马能谱，所述预处理包括能窗的选择、能谱平滑滤波、归一化处理以及非弹性散射信息扣除。

其中，所选择的能窗范围是30-210道，如图1所示。

能谱平滑滤波的方式是采用Savitzky-Golay滤波器，并用Savitzky-Golay五点滤波法，对所选择的能窗范围内的测量谱进行滤波，Savitzky-Golay五点滤波法的公式如下式1所示：

其中，

表示某一道址上滤波后的计数，y_i表示该道址上的计数，y_i-1表示该道址的前一道址上的计数，y_i-2表示该道址的前二道址上的计数，y_i+1表示该道址的后一道址上的计数，y_i+2表示该道址的后二道址上的计数。

归一化处理采用的方法是：对经过能窗范围选择、能谱平滑滤波之后的测量谱，使181道(逐道)的能谱数据之和为10，得到归一化处理后的测量谱，所采用的公式如下式2所示：

其中，N_Gkj为j深度点对应的归一化处理后的测量谱在第k道的计数，N_kj为j深度点对应的经过能窗范围选择、能谱平滑滤波之后的测量谱在第k道的计数，j为地层下的深度点。

非弹性散射信息扣除的方法是：对归一化处理后的测量谱选取30-54道、55-75道和76-210道三段道址的计数分别用0.9、0.7和0.8作为扣除系数进行非弹性信息扣除。

3.构建主元素群和辅助元素群

根据步骤1知道，DAS3井储层段3230m-3270m的主要元素种类为Si、Ca、Al、Fe、S、K、Na、Mg、Ti、Gd、H和Ba。据此构建的主元素群为：Si、Ca、S、Ti、Fe、Na、Gd；考虑到仪器本底的影响，在构建辅助元素群时增加了Tb元素，因此构建的辅助元素群为：Mg、K、Al和Tb。

4.基于元素群的最小二乘法解谱

根据步骤3得到的主元素群和辅助元素群，在经步骤2预处理后的元素俘获伽马能谱中，利用最小二乘法先对主元素(即：Si、Ca、S、Ti、Fe、Na、H、Ba、和Gd)进行解谱，然后从预处理后的总元素伽马能谱(即经步骤2预处理后的元素俘获伽马能谱)中扣除全部主元素所占贡献，再利用最小二乘法对辅助元素(即：Mg、K、Al和Tb)进行解谱，从而得到各元素的相对产额；

利用最小二乘法进行解谱是按照如下式3所示的算法进行的，以得到各元素的相对产额：

其中，y_j为第j种元素的相对产额(即该元素对总计数的贡献)，a_ij为第j种元素的归一化处理后的单元素的元素俘获标准伽马能谱在第i道的计数，c_i为预处理后的元素俘获伽马能谱在第i道的计数，ε_i为校正系数，ε_i≤0.1。

以DAS3井2403.1956m深度点获取的元素俘获伽马能谱为例，解谱得到的各元素的相对产额如图3所示。由于Td为本底产额，没有参考价值，故在图3中没有展示。

5.重构伽马能谱并控制误差

据步骤4得到的各元素的相对产额，结合归一化处理后的单元素的元素俘获标准伽马能谱，进行元素俘获伽马能谱的重构，并将重构的伽马能谱与测量伽马能谱进行对比，确定解谱结果是否可靠。

其中，所述归一化处理后的单元素的元素俘获标准伽马能谱是通过以下方法得到的：对于单元素的元素俘获标准伽马能谱，使256道(逐道)的能谱数据之和为10，所采用的公式如下式4所示：

其中，N_Gkj为j深度点对应的归一化处理后的单元素的元素俘获标准伽马能谱在第k道的计数，N_kj为j深度点对应的归一化处理前的单元素的元素俘获标准伽马能谱在第k道的计数，j为地层下的深度点。

考虑到前二道址的计数没有实际意义，因此也可以采用下式4-1对单元素的元素俘获标准伽马能谱进行归一化处理：

元素俘获伽马能谱的重构是按照下式5计算重构伽马能谱在第i道的计数(即所有解谱元素在第i道的计数)，得到各道的计数后，绘制重构伽马能谱：

其中，X_i为重构伽马能谱在第i道的计数，y_j为第j种元素的相对产额(由式3计算得到)，a_ij为第j种元素的归一化处理后的单元素的元素俘获标准伽马能谱在第i道的计数(由式4或式4-1计算得到)，ε_i为校正系数，ε_i≤0.1。

将重构的伽马能谱与测量伽马能谱进行对比是按照式6进行：

︱C_i－X_i︱《ε (式6)

c_i为测量伽马能谱在第i道的计数，X_i为重构的伽马能谱在第i道的计数，ε为相对误差；

当相对误差小于或等于5％时，认为解谱结果是可靠的；当相对误差大于5％，则认为解谱结果不可靠。

以DAS3井2403.1956m深度点各元素解谱结果为例，分别对重构伽马能谱和测量伽马能谱(图4)、预处理后的元素俘获伽马能谱和测量伽马能谱(图5)、主元素群重构伽马能谱(即按照式5计算全部主元素在第i道的计数，得到各道的计数后，重构的伽马能谱)和测量伽马能谱(图6)进行了对比，由图4-6可以看出，本实施例的解谱结果可靠，并且解谱方法精确、有效。

实施例2

本实施例提供了一种元素俘获能谱测井的伽马能谱解谱装置，其包括：

数据采集及主要元素类型确定模块，用于获取、分析研究区资料，所述资料至少包括元素俘获伽马能谱测井数据和全岩氧化物分析数据，确定该区块主要元素类型；

元素俘获伽马能谱测井数据预处理模块，用于对元素俘获伽马能谱测井数据进行预处理，所述预处理包括能窗的选择、能谱平滑滤波、归一化处理以及非弹性散射信息扣除，得到预处理后的元素俘获伽马能谱；

主元素群和辅助元素群构建模块，用于构建由不同元素组成的主元素群和辅助元素群，以明确不同元素的解谱先后顺序；

基于元素群的最小二乘法解谱模块，用于根据主元素群和辅助元素群，在预处理后的元素俘获伽马能谱中，利用最小二乘法先对主元素进行解谱，然后从预处理后的总元素伽马能谱(即预处理后的元素俘获伽马能谱)中扣除全部主元素所占贡献，再利用最小二乘法对辅助元素进行解谱，从而得到各元素的相对产额；

伽马能谱重构及误差控制模块，用于根据各元素的相对产额，结合归一化处理后的单元素的元素俘获标准伽马能谱，进行元素俘获伽马能谱的重构，并将重构的伽马能谱与测量伽马能谱进行对比，确定解谱结果是否可靠。

Claims (7)

1.一种元素俘获能谱测井的伽马能谱解谱方法，其包括以下步骤：

步骤1：获取、分析研究区资料，所述资料至少包括元素俘获伽马能谱测井数据和全岩氧化物分析数据，确定该研究区主要元素类型；

步骤2：对元素俘获伽马能谱测井数据进行预处理，所述预处理包括能窗的选择、能谱平滑滤波、归一化处理以及非弹性散射信息扣除，得到预处理后的元素俘获伽马能谱；

步骤3：构建由不同元素组成的主元素群和辅助元素群，以明确不同元素的解谱先后顺序；

步骤4：根据步骤3得到的主元素群和辅助元素群，在经步骤2预处理后的元素俘获伽马能谱中，利用最小二乘法先对主元素进行解谱，然后从经步骤2预处理后的元素俘获伽马能谱中扣除全部主元素所占贡献，再利用最小二乘法对辅助元素进行解谱，从而得到各元素的相对产额；

步骤5：根据步骤4得到的各元素的相对产额，结合归一化处理后的单元素的元素俘获标准伽马能谱，进行元素俘获伽马能谱的重构，并将重构的伽马能谱与测量伽马能谱进行对比，确定解谱结果是否可靠；

其中，在步骤2中，所选择的能窗范围是30-210道，归一化处理的方法是：对经过能窗范围选择、能谱平滑滤波之后的测量谱，使181道的能谱数据之和为10，得到归一化处理后的测量谱，所采用的公式如下式2所示：

其中，N_Gkj为j深度点对应的归一化处理后的测量谱在第k道的计数，N_kj为j深度点对应的经过能窗范围选择、能谱平滑滤波之后的测量谱在第k道的计数，k为道址，j为地层下的深度点；

非弹性散射信息扣除的方法是：对归一化处理后的测量谱选取30-54道、55-75道和76-210道三段道址的计数分别用0.9、0.7和0.8作为扣除系数进行非弹性信息扣除；

在步骤3中，所构建的主元素群包括：Si、Ca、S、H、Cl、Ti、Fe、Na、Ba和Gd中的一种或几种元素；所构建的辅助元素群包括：Mg、K、Cr、Ni、I、Tb和Al中的一种或几种元素，

在步骤5中，所述归一化处理后的单元素的元素俘获标准伽马能谱是通过以下方法得到的：对单元素的元素俘获标准伽马能谱，使256道的能谱数据之和为10，所采用的公式如下式4所示：

其中，N_Gkj为j深度点对应的归一化处理后的单元素的元素俘获标准伽马能谱在第k道的计数，N_kj为j深度点对应的归一化处理前的单元素的元素俘获标准伽马能谱在第k道的计数，k为道址，j为地层下的深度点。

2.根据权利要求1所述的解谱方法，其中，在步骤1中，所获取、分析的研究区资料还包括常规测井数据、录井数据和地质资料中的一种或几种的组合；研究区主要元素类型是通过以下方式确定的：通过分析研究区的全岩氧化物分析数据、以及可选择地，常规测井数据、录井数据和地质资料中的一种或几种的组合，得出不同种类氧化物所占岩石的重量百分比，依据各氧化物重量百分比确定该研究区主要元素类型。

3.根据权利要求1所述的解谱方法，其中，在步骤2中，能谱平滑滤波的方式是采用Savitzky-Golay滤波器，并用Savitzky-Golay五点滤波法，对所选择的能窗范围内的元素俘获伽马能谱测井数据进行滤波，Savitzky-Golay五点滤波法的公式如下式1所示：

其中，

表示某一道址上滤波后的计数，i为道址编号，y_i表示该道址上的计数，y_i-1表示该道址的前一道址上的计数，y_i-2表示该道址的前二道址上的计数，y_i+1表示该道址的后一道址上的计数，y_i+2表示该道址的后二道址上的计数。

4.根据权利要求1所述的解谱方法，其中，在步骤4中，利用最小二乘法进行解谱是按照如下式3所示的算法进行的，以得到各元素的相对产额：

其中，y_j为第j种元素的相对产额，i为道址编号，a_ij为第j种元素的归一化处理后的单元素的元素俘获标准伽马能谱在第i道的计数，c_i为预处理后的元素俘获伽马能谱在第i道的计数，ε_i为校正系数，ε_i≤0.1。

5.根据权利要求1所述的解谱方法，其中，在步骤5中，元素俘获伽马能谱的重构是按照下式5计算重构伽马能谱在第i道的计数，得到各道的计数后，绘制重构伽马能谱：

其中，i为道址编号，X_i为重构伽马能谱在第i道的计数，y_j为第j种元素的相对产额，a_ij为第j种元素的归一化处理后的单元素的元素俘获标准伽马能谱在第i道的计数，ε_i为校正系数，ε_i≤0.1。

6.根据权利要求1所述的解谱方法，其中，在步骤5中，将重构的伽马能谱与测量伽马能谱进行对比是按照式6进行：

︱C_i－X_i︱＜＜ε (式6)

i为道址编号，c_i为测量伽马能谱在第i道的计数，X_i为重构的伽马能谱在第i道的计数，ε为相对误差；

当相对误差小于或等于5％时，认为解谱结果是可靠的；当相对误差大于5％，则认为解谱结果不可靠。

7.一种元素俘获能谱测井的伽马能谱解谱装置，其包括：

数据采集及主要元素类型确定模块，用于获取、分析研究区资料，所述资料至少包括元素俘获伽马能谱测井数据和全岩氧化物分析数据，确定该研究区主要元素类型；

元素俘获伽马能谱测井数据预处理模块，用于对元素俘获伽马能谱测井数据进行预处理，所述预处理包括能窗的选择、能谱平滑滤波、归一化处理以及非弹性散射信息扣除，得到预处理后的元素俘获伽马能谱，

其中，所选择的能窗范围是30-210道，归一化处理的方法是：对经过能窗范围选择、能谱平滑滤波之后的测量谱，使181道的能谱数据之和为10，得到归一化处理后的测量谱，所采用的公式如下式2所示：

其中，N_Gkj为j深度点对应的归一化处理后的测量谱在第k道的计数，N_kj为j深度点对应的经过能窗范围选择、能谱平滑滤波之后的测量谱在第k道的计数，k为道址，j为地层下的深度点，

非弹性散射信息扣除的方法是：对归一化处理后的测量谱选取30-54道、55-75道和76-210道三段道址的计数分别用0.9、0.7和0.8作为扣除系数进行非弹性信息扣除；

主元素群和辅助元素群构建模块，用于构建由不同元素组成的主元素群和辅助元素群，以明确不同元素的解谱先后顺序，所构建的主元素群包括：Si、Ca、S、H、Cl、Ti、Fe、Na、Ba和Gd中的一种或几种元素；所构建的辅助元素群包括：Mg、K、Cr、Ni、I、Tb和Al中的一种或几种元素；

基于元素群的最小二乘法解谱模块，用于根据主元素群和辅助元素群，在预处理后的元素俘获伽马能谱中，利用最小二乘法先对主元素进行解谱，然后从预处理后的元素俘获伽马能谱中扣除全部主元素所占贡献，再利用最小二乘法对辅助元素进行解谱，从而得到各元素的相对产额；

伽马能谱重构及误差控制模块，用于根据各元素的相对产额，结合归一化处理后的单元素的元素俘获标准伽马能谱，进行元素俘获伽马能谱的重构，并将重构的伽马能谱与测量伽马能谱进行对比，确定解谱结果是否可靠，

所述归一化处理后的单元素的元素俘获标准伽马能谱是通过以下方法得到的：对单元素的元素俘获标准伽马能谱，使256道的能谱数据之和为10，所采用的公式如下式4所示：

其中，N_Gkj为j深度点对应的归一化处理后的单元素的元素俘获标准伽马能谱在第k道的计数，N_kj为j深度点对应的归一化处理前的单元素的元素俘获标准伽马能谱在第k道的计数，k为道址，j为地层下的深度点。

Images