CN112360428B - 一种针对中子寿命测井的俘获截面计算方法及装置 - Google Patents

一种针对中子寿命测井的俘获截面计算方法及装置 Download PDF

Info

Publication number
CN112360428B
CN112360428B CN202011091508.9A CN202011091508A CN112360428B CN 112360428 B CN112360428 B CN 112360428B CN 202011091508 A CN202011091508 A CN 202011091508A CN 112360428 B CN112360428 B CN 112360428B
Authority
CN
China
Prior art keywords
capture section
value
preset value
neutron
preset
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011091508.9A
Other languages
English (en)
Other versions
CN112360428A (zh
Inventor
张泉滢
张宫
邓瑞
刘国斌
刘佳怡
罗晓芸
朱可馨
廖绘梅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yangtze University
Original Assignee
Yangtze University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Yangtze University filed Critical Yangtze University
Priority to CN202011091508.9A priority Critical patent/CN112360428B/zh
Publication of CN112360428A publication Critical patent/CN112360428A/zh
Application granted granted Critical
Publication of CN112360428B publication Critical patent/CN112360428B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00—Survey of boreholes or wells
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V5/00—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
    • G01V5/04—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging
    • G01V5/08—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays
    • G01V5/10—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays using neutron sources
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12—Simultaneous equations, e.g. systems of linear equations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/30—Nuclear fission reactors

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Geology (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Mining & Mineral Resources (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Algebra (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geophysics (AREA)
  • Fluid Mechanics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Computing Systems (AREA)
  • Operations Research (AREA)
  • Measurement Of Radiation (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

本发明涉及一种针对中子寿命测井的俘获截面计算方法，包括以下步骤：设置初始俘获截面预设值数列，根据俘获截面预设值数列生成对应的中子寿命预设值数列；根据中子寿命预设值数列构建多指数衰减方程；从时间谱衰减区域中选取多个数据点代入多指数衰减方程，构建多阶多元线性超定方程组；采用正则化方法对非负最小二乘法的约束条件进行修正，得到线性超定方程组的等价方程组；求解等价方程组得到衰减项幅度值，进而构建俘获截面预设值的贡献分布图以获取俘获截面近似值；根据俘获截面近似值对俘获截面预设值数列进行优化，基于优化后的俘获截面预设值数列构建优化贡献分布图以获取最终的俘获截面值。本发明具有计算精度高的技术效果。

Description

一种针对中子寿命测井的俘获截面计算方法及装置
技术领域
本发明涉及中子寿命测井技术领域，尤其涉及一种针对中子寿命测井的俘获截面计算方法、装置及计算机存储介质。
背景技术
中子寿命测井通过记录脉冲中子源放出的高能快中子在地层中产生的热中子或俘获伽马时间谱来进行地层宏观俘获截面Σ和中子寿命τ测量，在储层含水饱和度计算、油水和气水界面划分等方面有广泛的应用，对油田剩余油动态监测和提高油气采收率具有非常重要的意义。
目前，中子寿命测井计算地层宏观俘获截面的主要方法有时间门法和双指数拟合法。时间门法通过在时间谱中选取两个固定时间窗，利用窗计数来计算时间谱曲线的指数衰减速率，进而求得地层的宏观俘获截面；双指数拟合法通过在时间谱中选取特定时间窗的信息进行函数拟合，得到地层宏观俘获截面。
然而，常规时间门法易受井眼信息的干扰，双指数拟合法对数据统计性要求较高，因此，现有方法对时间窗的选择要求极高，靠前或靠后的时间窗都会引入较大井眼干扰或者统计误差。此外，在复杂的测井条件下，中子和伽马时间谱形态和统计性会产生剧烈变化，基于固定时间窗的时间门法和双指数拟合法会产生较大的误差，严重降低了中子寿命测井在复杂测井条件下的精度。
发明内容
有鉴于此，有必要提供一种针对中子寿命测井的俘获截面计算方法、装置及计算机存储介质，用以解决现有时间谱处理方法在计算俘获截面时，计算精度受限于时间窗选择和测井条件的问题。
本发明提供一种针对中子寿命测井的俘获截面计算方法，包括以下步骤：
设置初始的俘获截面预设值数列，根据所述俘获截面预设值数列生成对应的中子寿命预设值数列；
根据所述中子寿命预设值数列，基于热中子和伽马时间谱指数衰减规律，构建描述时间谱计数率随时间衰减规律的多指数衰减方程；
从时间谱衰减区域中选取多个数据点，代入所述多指数衰减方程，构建多阶多元线性超定方程组；
采用正则化方法对非负最小二乘法的约束条件进行修正，基于修正后的约束条件，得到所述线性超定方程组的等价方程组；
采用非负约束最小二乘法求解所述等价方程组，得到衰减项幅度值，根据所述衰减项幅度值构建俘获截面预设值的贡献分布图，根据所述贡献分布图获取俘获截面近似值；
根据所述俘获截面近似值对所述俘获截面预设值数列进行二次优化，基于优化后的俘获截面预设值数列重新构建优化贡献分布图，根据所述优化贡献分布图获取俘获截面值。
进一步的，设置初始的俘获截面预设值数列，具体为：
设定地层中子宏观俘获截面初始范围，根据所述俘获截面初始范围生成一组地层宏观俘获截面预设值数列。
进一步的，根据所述俘获截面预设值数列生成对应的中子寿命预设值数列具体为：
τi＝4545.5/Σi
其中，∑i表示俘获截面预设值数列的第i个值，τi为中子寿命预设值数列的第i个值，i＝1,2,…,m，俘获截面预设值数列和中子寿命预设值数列的项数均为m。
进一步的，根据所述中子寿命预设值数列，基于热中子和伽马时间谱指数衰减规律，构建描述时间谱计数率随时间衰减规律的多指数衰减方程，具体为：
其中，c(t)为t时刻的时间谱计数率，τi为中子寿命预设值数列的第i个值，i＝1,2,…,m，m为中子寿命预设值数列的项数，ai速率1/τi对应的指数衰减项幅度值。
进一步的，从时间谱衰减区域中选取多个数据点，代入所述多指数衰减方程，构建多阶多元线性超定方程组，具体为：
其中，ck表示tk时刻的数据点的时间谱计数率，k＝1,2,…,n，n为选取的数据点的数量，τi为中子寿命预设值数列的第i个值，ai为速率1/τi对应的指数衰减项幅度值，i＝1,2,…,m，m为中子寿命预设值数列的项数；
当n＞m时，所述线性超定方程组为n阶m元超定方程组，其矩阵形式为：
其中，A为衰减项幅度值矩阵，H为系数矩阵，C为时间谱计数率矩阵。
进一步的，采用正则化方法对非负最小二乘法的约束条件进行修正，具体为：
其中，λ为正则化因子，hik＝exp(-tk/τi)，即系数矩阵H第i行第k列的元素，ai为速率1/τi对应的指数衰减项幅度值，即衰减项幅度值矩阵A的第i个元素，ck表示tk时刻的数据点的时间谱计数率，即时间谱计数率矩阵C的第k个元素；
基于修正后的约束条件，得到所述线性超定方程组的等价方程组，具体为：
其中，I为m×m阶单位对角矩阵。
进一步的，根据所述贡献分布图获取俘获截面近似值，具体为：
获取所述贡献分布图的第一个贡献峰峰值对应的俘获截面预设值作为俘获截面近似值；
根据所述优化贡献分布图获取俘获截面值，具体为：
获取优化贡献分布图中第一个贡献峰峰值对应的俘获截面预设值作为俘获截面值。
进一步的，根据所述俘获截面近似值对所述俘获截面预设值数列进行二次优化，具体为：
以所述俘获截面近似值作为中心点设定加密区间，设定加密间隔，所述加密间隔小于所述俘获截面预设值数列的公差；
基于所述加密间隔所述俘获截面预设值数列中位于所述加密区间内的项进行加密处理，得到优化后的俘获截面预设值数列。
本发明还提供一种针对中子寿命测井的俘获截面计算装置，包括处理器以及存储器，所述存储器上存储有计算机程序，所述计算机程序被所述处理器执行时，实现所述针对中子寿命测井的俘获截面计算方法。
本发明还提供一种计算机存储介质，其上存储有计算机程序，所述计算机该程序被处理器执行时，实现所述针对中子寿命测井的俘获截面计算方法。
有益效果：本发明基于俘获截面预设值数列和多指数衰减方程构建了多阶多元线性超定方程组，通过求解方程组构建了俘获截面预设值贡献分布图，成功排除了井眼信息对计算结果的影响。同时，为了消除时间谱统计性对计算结果的影响，本发明没有直接采用非负约束最小二乘法对超定方程组进行求解，而是采用正则化方法对非负约束最小二乘法的约束条件进行了修正，基于修正后的约束条件得到线性超定方程组的等价方程组，然后再进行等价方程组的求解，最后得到俘获截面值。本发明在计算地层宏观俘获截面过程中成功排除了井眼信息和时间谱统计性对计算结果的干扰，计算结果不再受时间窗选择的限制，大大提高了中子寿命测井在复杂地层条件下的精度。
附图说明
图1为本发明提供的针对中子寿命测井的俘获截面计算方法第一实施例的方法流程图；
图2为图1步骤S3中用于选取数据点的热中子时间谱图；
图3为图1步骤S5中贡献分布图；
图4为图1步骤S6中优化贡献分布图。
具体实施方式
下面结合附图来具体描述本发明的优选实施例，其中，附图构成本申请一部分，并与本发明的实施例一起用于阐释本发明的原理，并非用于限定本发明的范围。
实施例1
如图1所示，本发明的实施例1提供了针对中子寿命测井的俘获截面计算方法，以下简称本方法，包括以下步骤：
S1、设置初始的俘获截面预设值数列，根据所述俘获截面预设值数列生成对应的中子寿命预设值数列；
S2、根据所述中子寿命预设值数列，基于热中子和伽马时间谱指数衰减规律，构建描述时间谱计数率随时间衰减规律的多指数衰减方程；
S3、从时间谱衰减区域中选取多个数据点，代入所述多指数衰减方程，构建多阶多元线性超定方程组；
S4、采用正则化方法对非负最小二乘法的约束条件进行修正，基于修正后的约束条件，得到所述线性超定方程组的等价方程组；
S5、采用非负约束最小二乘法求解所述等价方程组，得到衰减项幅度值，根据所述衰减项幅度值构建俘获截面预设值的贡献分布图，根据所述贡献分布图获取俘获截面近似值；
S6、根据所述俘获截面近似值对所述俘获截面预设值数列进行二次优化，基于优化后的俘获截面预设值数列重新构建优化贡献分布图，根据所述优化贡献分布图获取俘获截面值。
本发明公开了一种基于正则化最小二乘法的中子寿命测井俘获截面计算方法；首先，根据常见地层俘获截面范围预先设置一组宏观俘获截面和中子寿命数组；然后，利用给定的中子寿命数组构建多指数衰减方程描述中子(或伽马)时间谱衰减规律；其次，选用中子(或伽马)时间谱数据构建关于指数衰减项幅度值的超定方程组；为了消除时间谱统计性对计算结果的影响，利用正则化方法对超定方程组求解条件进行修正，形成基于正则化非负最小二乘法的等价方程组；最后，利用非负最小二乘法对等价方程组进行求解，得到关于地层俘获截面预设值的贡献分布，进而得到地层宏观俘获截面和中子寿命的计算结果。该方法与现有的时间谱处理方法相比较，在时间谱数据处理方法不再受时间窗选择的影响，并很好地解决了井眼和计数统计性问题。
现有方法计算结果受限于时间窗选择的原因，是因为方法本身处理过程不能排除井眼或者时间谱统计性的影响；时间谱的前半段数据统计性好但井眼信息比重大，后半段数据计数统计性差但地层信息比重大；因此，在选择数据的时候需要选择井眼信息比重小且计数统计性好的数据点来进行处理，这就是需要选择时间窗的原因。
而本方法通过求解多元超定方程组得到中子俘获截面预设值贡献分布，成功分离了时间谱中井眼和地层俘获截面的贡献，从而排除了井眼干扰；同时，采用正则化方法对非负最小二乘法修正，排除了时间谱数据统计性对计算结果的干扰。由于本方法具有自动排除井眼和时间谱统计性干扰的优势，在选择时间谱数据方面不用考虑井眼信息比重和计数统计性好坏，摆脱了时间窗选择的限制，可以直接对整个时间谱衰减区域进行处理，因此，本方法适用于任何复杂测井条件的时间谱。
与现有方法相比，本方法在计算热中子寿命和宏观俘获截面过程中不再受时间窗选择的影响，并很好地排除了井眼和计数统计性对计算结果的影响，大大提高了中子寿命测井在复杂测井条件下的应用效果。
优选的，设置初始的俘获截面预设值数列，具体为：
设定地层中子宏观俘获截面初始范围，根据所述俘获截面初始范围生成一组地层宏观俘获截面预设值数列。
首先根据常见地层中子宏观俘获截面范围，设定地层中子宏观俘获截面初始范围(∑min,∑max)，其中，∑min为俘获截面最小值，∑max为俘获截面最大值；具体的，本实施例中俘获截面初始范围选取1cu到200cu；
根据俘获截面初始范围生成一组地层宏观俘获截面预设值数列
其中，∑i表示俘获截面预设值数列的第i个值，d为俘获截面预设值数列的公差，m为俘获截面预设值数列的项数；
本实施例中m取101，∑i＝10+(i-1)d,d＝(200-1)/(101-1)；(i＝1,2,…,101)。
优选的，根据所述俘获截面预设值数列生成对应的中子寿命预设值数列具体为：
τi＝4545.5/Σi
其中，∑i表示俘获截面预设值数列的第i个值，τi为中子寿命预设值数列的第i个值，i＝1,2,…,m，俘获截面预设值数列和中子寿命预设值数列的项数均为m，m＝101。
根据地层宏观俘获截面和中子寿命的数学关系，生成对应的中子寿命预设值数列。
优选的，根据所述中子寿命预设值数列，基于热中子和伽马时间谱指数衰减规律，构建描述时间谱计数率随时间衰减规律的多指数衰减方程，具体为：
其中，c(t)为t时刻的时间谱计数率，ai为速率1/τi对应的指数衰减项幅度值，ai为非负值，其大小可以反映对应的衰减项在中子衰减过程中的贡献。
优选的，从时间谱衰减区域中选取多个数据点，代入所述多指数衰减方程，构建多阶多元线性超定方程组，具体为：
其中，ck表示tk时刻的数据点的时间谱计数率，k＝1,2,…,n，n为选取的数据点的数量，τi为中子寿命预设值数列的第i个值，ai为速率1/τi对应的指数衰减项幅度值，i＝1,2,…,m，m为中子寿命预设值数列的项数。
当n＞m时，所述线性超定方程组为n阶m元超定方程组，其矩阵形式为：
其中，A为衰减项幅度值矩阵，H为系数矩阵，C为时间谱计数率矩阵；
Am×1＝(a1,a2,…,am-1,am)T；
Cn×1＝(c1,c2,…,cn-1,cn)T；
具体的，本实施例以图2中40p.u.饱含水石灰岩中的热中子时间谱(中子俘获截面参考值为12.88cu)为例，图2为不同地层条件下的热中子时间谱，记录时间为0-2000μs，时间道数为200，时间道间隔为10μs；其中，在时间谱100-500μs时间段内的数据井眼信息比重较大；在矿化度为100kppm和300kppm的情况下，1200-1800μs时间段的数据存在计数统计性较差的问题；图2中中子计数率即时间谱计数率。从图2的热中子时间谱中选取数据点构建超定方程组，在时间谱衰减部分选取数据点，本实施例选取100μs到1800μs时间段内的数据，共计171个数据点，利用构造的多指数衰减方程，得到一个171阶101元超定方程组，
优选的，采用正则化方法对非负最小二乘法的约束条件进行修正，具体为：
修正前的约束条件为：
修正后的约束条件为：
其中，λ为正则化因子，hik＝exp(-tk/τi)，即系数矩阵H第i行第k列的元素，ai为速率1/τi对应的指数衰减项幅度值，即衰减项幅度值矩阵A的第i个元素，ck表示tk时刻的数据点的时间谱计数率，即时间谱计数率矩阵C的第k个元素；
基于修正后的约束条件，得到所述线性超定方程组的等价方程组，具体为：
其中，I为m×m阶单位对角矩阵。
具体的，本实施例中：
A101×1＝(a1,a2,...,a100,a101)T；
C171×1＝(c1,c2,...,c170,c171)T；
优选的，根据所述贡献分布图获取俘获截面近似值，具体为：
获取所述贡献分布图的第一个贡献峰峰值对应的俘获截面预设值作为俘获截面近似值；
根据所述优化贡献分布图获取俘获截面值，具体为：
获取优化贡献分布图中第一个贡献峰峰值对应的俘获截面预设值作为俘获截面值。
具体的，选取合适的正则化因子λ，本实施例中λ取0.001，采用非负约束最小二乘法获取等价方程组的最优解(即m个指数衰减项对应的幅度值)。根据衰减项幅度值，得到关于中子俘获截面预设值的贡献分布图，如图3所示。根据中子(伽马)时间谱衰减规律可知，在整个时间谱衰减区域中，地层部分对应的时间谱衰减速率是最低且最稳定的。因此，在所有存在贡献的衰减项中，地层衰减项对应的俘获截面预设值应该是最小的。鉴于上述分析，从中子俘获截面贡献分布图中选取第一个贡献峰峰值所对应的地层宏观俘获截面预设值作为地层宏观俘获截面的近似值∑a，图3中第一个贡献峰峰值所对应的地层宏观俘获截面预设值为13.06cu。
优选的，根据所述俘获截面近似值对所述俘获截面预设值数列进行二次优化，具体为：
以所述俘获截面近似值作为中心点设定加密区间，设定加密间隔，所述加密间隔小于所述俘获截面预设值数列的公差；
基于所述加密间隔所述俘获截面预设值数列中位于所述加密区间内的项进行加密处理，得到优化后的俘获截面预设值数列。
以宏观俘获截面近似值∑a为参照，对地层宏观俘获截面预设值数列重新优化，对近似值所在的加密区间范围进行预设值加密处理。本实施例选取13.06cu±2cu区间进行加密，加密间隔0.25cu，得到加密优化后的俘获截面预设值数列，然后基于优化后的俘获截面预设值数列重新获取贡献分布图，得到优化贡献分布图，如图4所示。优化贡献分布图的获取方法与贡献分布图的获取方法一致。选取优化贡献分布图中的地层宏观俘获截面的近似值作为地层宏观俘获截面值，进而可以得到对应的中子寿命值，图4中第一个贡献峰峰值所对应的地层宏观俘获截面预设值为12.81c.u.，即俘获截面值为12.81c.u.。
利用本方法，采用图2中不同时间段(例如200μs-1800μs和300μs-1400μs)的数据对对矿化度为100kppm以及300kppm的饱含水石灰岩地层时间谱进行处理，得到地层宏观俘获截面值进行对比，如表1所示。
表1、采用不同时间段的数据得到中子俘获截面结果
从表1可以看出，两个时间段的数据在井眼影响和数据统计性方面存在很大的差异，但是本方法利用上述两组数据得到的俘获截面结果与地层参考值基本一致。此外，在时间谱统计性较差的情况下(例如图2中40p.u.的100kppm和300kppm饱含水地层)，本方法计算依然能够得到准确的计算结果。这表明本发明在时间谱处理过程中不再受时间窗选择的限制，并很好地解决了井眼和计数统计性问题。
本方法在计算热中子寿命和宏观俘获截面过程中不再受时间窗选择的影响，降低了井眼和时间谱统计性对计算结果的影响，提高了中子寿命测井在复杂地层条件下的精度。
实施例2
本发明的实施例2提供了针对中子寿命测井的俘获截面计算装置，包括处理器以及存储器，所述存储器上存储有计算机程序，所述计算机程序被所述处理器执行时，实现实施例1提供的针对中子寿命测井的俘获截面计算方法。
本发明实施例提供的针对中子寿命测井的俘获截面计算装置，用于实现针对中子寿命测井的俘获截面计算方法，因此，针对中子寿命测井的俘获截面计算方法所具备的技术效果，针对中子寿命测井的俘获截面计算装置同样具备，在此不再赘述。
实施例3
本发明的实施例3提供了计算机存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时，实现实施例1提供的针对中子寿命测井的俘获截面计算方法。
本发明实施例提供的计算机存储介质，用于实现针对中子寿命测井的俘获截面计算方法，因此，针对中子寿命测井的俘获截面计算方法所具备的技术效果，计算机存储介质同样具备，在此不再赘述。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。

Claims (8)

1.一种针对中子寿命测井的俘获截面计算方法， 其特征在于，包括以下步骤：
设置初始的俘获截面预设值数列，根据所述俘获截面预设值数列生成对应的中子寿命预设值数列；
根据所述中子寿命预设值数列，基于热中子和伽马时间谱指数衰减规律，构建描述时间谱计数率随时间衰减规律的多指数衰减方程；
从时间谱衰减区域中选取多个数据点，代入所述多指数衰减方程，构建多阶多元线性超定方程组；
采用正则化方法对非负最小二乘法的约束条件进行修正，基于修正后的约束条件，得到所述线性超定方程组的等价方程组；
采用非负约束最小二乘法求解所述等价方程组，得到衰减项幅度值，根据所述衰减项幅度值构建俘获截面预设值的贡献分布图，根据所述贡献分布图获取俘获截面近似值；
根据所述俘获截面近似值对所述俘获截面预设值数列进行二次优化，基于优化后的俘获截面预设值数列重新构建优化贡献分布图，根据所述优化贡献分布图获取俘获截面值；
从时间谱衰减区域中选取多个数据点，代入所述多指数衰减方程，构建多阶多元线性超定方程组，具体为：
；
其中，表示/>时刻的数据点的时间谱计数率，/>，/>为选取的数据点的数量，/>为中子寿命预设值数列的第/>个值，/>为速率/>对应的指数衰减项幅度值，/>，/>为中子寿命预设值数列的项数；
当时，所述线性超定方程组为/>阶/>元超定方程组，其矩阵形式为：
；
其中，为衰减项幅度值矩阵，/>为系数矩阵，/>为时间谱计数率矩阵；
采用正则化方法对非负最小二乘法的约束条件进行修正，具体为：
其中，为正则化因子，/>，即系数矩阵/>第/>行第/>列的元素，/>为速率对应的指数衰减项幅度值，即衰减项幅度值矩阵/>的第/>个元素，/>表示/>时刻的数据点的时间谱计数率，即时间谱计数率矩阵/>的第/>个元素；
基于修正后的约束条件，得到所述线性超定方程组的等价方程组，具体为：
其中，为/>阶单位对角矩阵。
2.根据权利要求1所述的针对中子寿命测井的俘获截面计算方法，其特征在于，设置初始的俘获截面预设值数列，具体为：
设定地层中子宏观俘获截面初始范围，根据所述俘获截面初始范围生成一组地层宏观俘获截面预设值数列。
3.根据权利要求1所述的针对中子寿命测井的俘获截面计算方法，其特征在于，根据所述俘获截面预设值数列生成对应的中子寿命预设值数列，具体为：
其中，表示俘获截面预设值数列的第/>个值，/>为中子寿命预设值数列的第/>个值，，俘获截面预设值数列和中子寿命预设值数列的项数均为/>。
4.根据权利要求1所述的针对中子寿命测井的俘获截面计算方法，其特征在于，根据所述中子寿命预设值数列，基于热中子和伽马时间谱指数衰减规律，构建描述时间谱计数率随时间衰减规律的多指数衰减方程，具体为：
其中，为/>时刻的时间谱计数率，/>为中子寿命预设值数列的第/>个值，/>，为中子寿命预设值数列的项数，/>为速率/>对应的指数衰减项幅度值。
5.根据权利要求1所述的针对中子寿命测井的俘获截面计算方法，其特征在于，根据所述贡献分布图获取俘获截面近似值，具体为：
获取所述贡献分布图的第一个贡献峰峰值对应的俘获截面预设值作为俘获截面近似值；
根据所述优化贡献分布图获取俘获截面值，具体为：
获取优化贡献分布图中第一个贡献峰峰值对应的俘获截面预设值作为俘获截面值。
6.根据权利要求1所述的针对中子寿命测井的俘获截面计算方法，其特征在于，根据所述俘获截面近似值对所述俘获截面预设值数列进行二次优化，具体为：
以所述俘获截面近似值作为中心点设定加密区间，设定加密间隔，所述加密间隔小于所述俘获截面预设值数列的公差；
基于所述加密间隔所述俘获截面预设值数列中位于所述加密区间内的项进行加密处理，得到优化后的俘获截面预设值数列。
7.一种针对中子寿命测井的俘获截面计算装置，其特征在于，包括处理器以及存储器，所述存储器上存储有计算机程序，所述计算机程序被所述处理器执行时，实现如权利要求1-6任一项所述的针对中子寿命测井的俘获截面计算方法。
8.一种计算机存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时，实现如权利要求1-6任一项所述的针对中子寿命测井的俘获截面计算方法。
CN202011091508.9A 2020-10-13 2020-10-13 一种针对中子寿命测井的俘获截面计算方法及装置 Active CN112360428B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011091508.9A CN112360428B (zh) 2020-10-13 2020-10-13 一种针对中子寿命测井的俘获截面计算方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011091508.9A CN112360428B (zh) 2020-10-13 2020-10-13 一种针对中子寿命测井的俘获截面计算方法及装置

Publications (2)

Publication Number Publication Date
CN112360428A CN112360428A (zh) 2021-02-12
CN112360428B true CN112360428B (zh) 2024-02-20

Family

ID=74507913

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011091508.9A Active CN112360428B (zh) 2020-10-13 2020-10-13 一种针对中子寿命测井的俘获截面计算方法及装置

Country Status (1)

Country Link
CN (1) CN112360428B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113669048B (zh) * 2021-07-30 2024-05-24 长江大学 热中子时间衰减测井宏观俘获截面值提取方法
CN115098823A (zh) * 2022-06-29 2022-09-23 西南石油大学 一种中子寿命测井地层宏观俘获截面提取方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5374823A (en) * 1993-10-28 1994-12-20 Computalog U.S.A., Inc. Pulsed neutron decay tool for measuring gamma radiation energy spectra for fast neutron inelastic collisions and thermal neutron capture events
CN101078775A (zh) * 2007-06-29 2007-11-28 西安奥华电子仪器有限责任公司 脉冲中子双谱饱和度测井方法
US7408150B1 (en) * 2007-06-25 2008-08-05 Schlumberger Technology Corporation Well logging method for determining formation characteristics using pulsed neutron capture measurements
CN104114811A (zh) * 2011-09-08 2014-10-22 卡博陶粒有限公司 定位诱导的地下岩层裂缝中的示踪支撑剂的与岩性和井眼状况无关的方法
CN111032997A (zh) * 2017-12-26 2020-04-17 中国石油天然气股份有限公司 一种元素俘获能谱测井的伽马能谱解谱方法及装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5374823A (en) * 1993-10-28 1994-12-20 Computalog U.S.A., Inc. Pulsed neutron decay tool for measuring gamma radiation energy spectra for fast neutron inelastic collisions and thermal neutron capture events
US7408150B1 (en) * 2007-06-25 2008-08-05 Schlumberger Technology Corporation Well logging method for determining formation characteristics using pulsed neutron capture measurements
CN101078775A (zh) * 2007-06-29 2007-11-28 西安奥华电子仪器有限责任公司 脉冲中子双谱饱和度测井方法
CN104114811A (zh) * 2011-09-08 2014-10-22 卡博陶粒有限公司 定位诱导的地下岩层裂缝中的示踪支撑剂的与岩性和井眼状况无关的方法
CN111032997A (zh) * 2017-12-26 2020-04-17 中国石油天然气股份有限公司 一种元素俘获能谱测井的伽马能谱解谱方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
PNN测井热中子时间谱的多尺度滤波研究;房文静;范宜仁;张峰;;核电子学与探测技术(03);第385-389页 *

Also Published As

Publication number Publication date
CN112360428A (zh) 2021-02-12

Similar Documents

Publication Publication Date Title
CN112360428B (zh) 一种针对中子寿命测井的俘获截面计算方法及装置
US4390783A (en) Spectroscopic analysis with background compensation
US4721853A (en) Thermal decay time logging method and apparatus
US6737642B2 (en) High dynamic range analog-to-digital converter
CN109799541A (zh) 一种γ能谱测量谱漂及积分非线性校正方法
US20200292717A1 (en) Spectral gamma ray downhole logging tool
CN100498377C (zh) 模拟到数字转换的移位误差校正
Radzevičius et al. Upper Wenlock δ 13 C chemostratigraphy, conodont biostratigraphy and palaeoecological dynamics in the Ledai-179 drill core (Eastern Lithuania).
US4388529A (en) System for simultaneous measurement of borehole and formation neutron lifetimes
US6124590A (en) Method for determining thermal neutron capture cross-section of earth formations using measurements from multiple capture gamma ray detectors
US6327549B1 (en) Differential correction method and apparatus
KR920003040A (ko) 랜덤 및 펄스형 특성을 지니거나 펄스형으로 변환될 수 있는 물리량의 측정방법 및 그의 감마선 분광 측정법에의 적용
US5491342A (en) Apparatus and method for nuclear camera calibration
CN108646285B (zh) 伽马射线能谱测量方法、装置及系统
DE2941535A1 (de) Verfahren und vorrichtung fuer die optimierung von messungen der daempfungscharakteristiken thermischer neutronen
JP3268251B2 (ja) 対象核種の濃度の連続測定方法
Gin et al. Development of a technique for high-speed γ-ray spectrometry
Ameli et al. Optical background measurement in a potential site for the NEMO KM undersea neutrino telescope
Wouters et al. Beta-decay energies and masses of In 1 0 3− 1 0 5
EP0045686A1 (en) Spectroscopic analysis with background compensation
Dremin et al. Charm in cosmic rays (The long-flying component of EAS cores)
WO2024076861A2 (en) Spectral alignment method for induced gamma ray logging
CN106324692A (zh) 一种地面伽玛能谱测量铀含量数据调平校正方法
Stoyanova et al. A study on the transition between seniority-type and collective excitations in 204Po and 206Po
Gamba et al. Fast-Timing measurements in 100zr using labr3 (ce) detectors coupled with gammasphere

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant