CN103235350B - 放射性测井仪器稳定性检测及刻度方法与装置 - Google Patents
放射性测井仪器稳定性检测及刻度方法与装置 Download PDFInfo
- Publication number
- CN103235350B CN103235350B CN201310126466.1A CN201310126466A CN103235350B CN 103235350 B CN103235350 B CN 103235350B CN 201310126466 A CN201310126466 A CN 201310126466A CN 103235350 B CN103235350 B CN 103235350B
- Authority
- CN
- China
- Prior art keywords
- well logging
- logging instrument
- collection result
- normal distribution
- radioactivity well
- 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
Links
Landscapes
- Measurement Of Radiation (AREA)
Abstract
本发明公开了一种放射性测井仪器稳定性检测及刻度方法与装置,克服目前放射性测井仪器稳定性检测效率较低且可靠性较差的不足,该放射性测井仪器稳定性检测方法包括:放射性测井仪器位于套管中的第一深度时,采集符合正态分布的p次计数获得第一采集结果;计算第一采集结果的均值、方差及变异系数;放射性测井仪器位于套管中的第二深度时,采集符合正态分布的p次计数,获得第二采集结果;计算第二采集结果的均值、方差及变异系数;根据第一采集结果的均值、方差及变异系数,以及第二采集结果的均值、方差及变异系数,确定放射性测井仪器的稳定性。与现有技术相比,本申请的实施例可以定量评价放射性测井仪器的质量和性能,效率高、可靠性高。
Description
技术领域
本发明涉及放射性测井仪器的检测技术,尤其涉及一种放射性测井仪器稳定性检测及刻度方法与装置。
背景技术
放射性测井指在钻孔中利用岩石天然放射性、γ射线与物质的相互作用,以及中子与物质的相互作用等一系列效应,来研究地层情况的一组测井方法。
放射性测井主要包括自然伽马测井、伽马测井以及中子测井等。这些测井得到的数据都是等时间间隔内的计数,均符合放射性衰变的随机性、统计性。
在放射性测量中,如果测量仪器不稳定,存在系统误差,则多次测量的计数率误差将不符合泊松分布。如果仪器稳定,不存在仪器系统误差,则多次测量的计数率误差来自放射性涨落,在平均值较小时,服从泊松分布;平均值较大时,服从高斯分布。一般情况下,放射性测井中的计数率较大,因此,服从高斯分布。
稳定性是仪表在规定工作条件内,某些性能随时间保持不变的能力。在放射性测井中,仪器的稳定性对测量结果的准确性至关重要。
在放射性测井仪器进行刻度或测井时,仪器加电后需要一段时间才能稳定测量。而这个时间主要靠操作人员的经验来判断,往往造成仪器尚未达到稳定状态就开始进行测量的情况,导致仪器刻度不能通过或刻度不够准确。
另外,在仪器质量检验中,生产人员往往通过是否刻度通过来判断仪器是否合格,无法定量判断仪器稳定时的状态以及仪器能稳定工作的最长时间。
发明内容
本发明所要解决的技术问题是克服目前放射性测井仪器稳定性检测效率较低且可靠性较差的不足。
为了解决上述技术问题,本发明提供了一种放射性测井仪器稳定性检测方法,所述放射性测井仪器中设置有放射源;该方法包括:
所述放射性测井仪器位于套管中的第一深度时,采集所述放射源符合正态分布的p次计数,获得第一采集结果;其中p大于等于8;
计算所述第一采集结果的均值、方差及变异系数;
所述放射性测井仪器位于所述套管中的第二深度时,采集所述放射源符合正态分布的p次计数,获得第二采集结果;
计算所述第二采集结果的均值、方差及变异系数;
根据所述第一采集结果的均值、方差及变异系数,以及所述第二采集结果的均值、方差及变异系数,确定所述放射性测井仪器的稳定性。
优选地,采集所述放射源符合正态分布的p次计数,获得第一采集结果,包括:
在所述第一深度时采集所述放射源的至少一轮p次计数,利用X2检验方法对所述第一深度时采集到的至少一轮p次计数是否符合正态分布进行判断,将其中符合正态分布的一轮p次计数作为所述第一采集结果;
采集所述放射源符合正态分布的p次计数,获得第二采集结果,包括:
在所述第二深度时采集所述放射源的至少一轮p次计数,利用X2检验方法对所述第二深度时采集到的至少一轮p次计数是否符合正态分布进行判断,将其中符合正态分布的一轮p次计数作为所述第二采集结果。
优选地,根据所述第一采集结果的均值、方差及变异系数,以及所述第二采集结果的均值、方差及变异系数,确定所述放射性测井仪器的稳定性,包括:
利用T检验和F检验判断出所述第一采集结果的均值与所述第二采集结果的均值、所述第一采集结果的方差与所述第二采集结果的方差以及所述第一采集结果的变异系数与所述第二采集结果的变异系数是否有统计上的显著差异时,确认所述放射性测井仪器具备稳定性。
本申请还提供了一种放射性测井仪器刻度方法,该方法包括:
获得所述放射性测井仪器在多块刻度器中各自符合正态分布的p次测量数据;其中,p大于等于8;
所述放射性测井仪器在第i+1块刻度器中获得的符合正态分布的p次测量数据N[i+1,p]与所述放射性测井仪器在第一块刻度器中获得的符合正态分布的p次测量数据N[1,p]没有统计上的显著差异时,根据所述测量数据N[i+1,p]进行刻度参数的计算,完成所述放射性测井仪器在所述第i+1块刻度器中的刻度;其中,i大于等于1。
优选地,利用T检验和F检验判断所述放射性测井仪器在第i块刻度器中获得的符合正态分布的p次测量数据N[i+1,p]与所述放射性测井仪器在第一块刻度器中获得的符合正态分布的p次测量数据N[1,p]是否具有统计上的显著差异。
本申请还提供了一种放射性测井仪器稳定性检测方法,该方法包括:
采集所述放射性测井仪器位于井中或刻度器中时符合正态分布的p次测量数据N[1,p]并保存,记录保存时刻为第一时刻;其中p大于等于8;
采集第i+1次所述放射性测井仪器位于井中或刻度器中时符合正态分布的p次测量数据N[i+1,p]并保存,记录保存时刻为第i+1时刻;其中,i大于等于1;
在从N[1,p]到N[i+1,p]这(i+1)×p个数据符合正态分布,且所述p次测量数据N[i+1,p]与所述p次测量数据N[1,p]没有统计上的显著差异时,将所述第i+1时刻与所述第一时刻的差值作为所述放射性测井仪器的最大稳定时间。
优选地,利用T检验和F检验判断N[1,p]与N[i+1,p]是否具有统计上的显著差异。
本申请还提供了一种放射性测井仪器稳定性检测装置,所述放射性测井仪器中设置有放射源;该装置包括:
采集模块,设置为所述放射性测井仪器位于套管中的第一深度时,采集所述放射源符合正态分布的p次计数,获得第一采集结果;所述放射性测井仪器位于所述套管中的第二深度时,采集所述放射源符合正态分布的p次计数,获得第二采集结果;其中p大于等于8;
计算模块,设置为计算所述第一采集结果的均值、方差及变异系数;计算所述第二采集结果的均值、方差及变异系数;
判断模块,设置为根据所述第一采集结果的均值、方差及变异系数,以及所述第二采集结果的均值、方差及变异系数,确定所述放射性测井仪器的稳定性。
本申请还提供了一种放射性测井仪器刻度装置,该装置包括:
获取模块,设置为获得所述放射性测井仪器在多块刻度器中各自符合正态分布的p次测量数据;其中,p大于等于8;
刻度模块,设置为所述放射性测井仪器在第i+1块刻度器中获得的符合正态分布的p次测量数据N[i+1,p]与所述放射性测井仪器在第一块刻度器中获得的符合正态分布的p次测量数据N[1,p]没有统计上的显著差异时,根据所述测量数据N[i+1,p]进行刻度参数的计算,完成所述放射性测井仪器在所述第i+1块刻度器中的刻度;其中,i大于等于1。
本申请还提供了一种放射性测井仪器稳定性检测装置,该装置包括:
采集模块,设置为采集所述放射性测井仪器位于井中或刻度器中时符合正态分布的p次测量数据N[1,p]并保存,记录保存时刻为第一时刻;采集第i+1次所述放射性测井仪器位于井中或刻度器中时符合正态分布的p次测量数据N[i+1,p]并保存,记录保存时刻为第i+1时刻;其中,i大于等于1,p大于等于8;
确定模块,设置为在从N[1,p]到N[i+1,p]这(i+1)×p个数据符合正态分布,且所述p次测量数据N[i+1,p]与所述p次测量数据N[1,p]没有统计上的显著差异时,将所述第i+1时刻与所述第一时刻的差值作为所述放射性测井仪器的最大稳定时间。
与现有技术相比,本申请的实施例可以定量评价放射性测井仪器的质量和性能,效率高、可靠性高。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
附图用来提供对本发明技术方案的进一步理解,并且构成说明书的一部分,与本申请的实施例一起用于解释本发明的技术方案,并不构成对本发明技术方案的限制。
图1为本申请实施例放射性测井仪器稳定性检测方法的流程示意图。
图2为本申请实施例放射性测井仪器刻度方法的流程示意图。
图3为本申请实施例放射性测井仪器稳定性检测方法的流程示意图。
图4为本申请实施例放射性测井仪器稳定性检测装置的构造示意图。
图5为本申请实施例放射性测井仪器刻度装置的构造示意图。
图6为本申请实施例放射性测井仪器稳定性检测装置的构造示意图。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式,借此对本发明如何应用技术手段来解决技术问题,并达成技术效果的实现过程能充分理解并据以实施。本申请实施例以及实施例中的各个特征在不相冲突前提下的相互结合,均在本发明的保护范围之内。
高斯分布又称为正态分布,对于服从高斯分布的总体分布来说,其概率密度定义如下:
其中σ为标准差,m为期望。
有限次测量的标准差用它们的标准偏差s来表示:
k为测量次数,Δi为第i次测量值与平均值之间的偏差。
变异系数等于标准差与期望的绝对值之比σ/m。
χ2分布检验方法中,设有n个独立变量e,每个随机变量都服从标准化的高斯分布(期望值为0,方差为1),则统计量χ2=也是一个随机变量,并服从自由度为v(独立变量数减去约束条件)的χ2分布:
其中:为Γ函数。
χ2值大于的概率为:
对于给定的自由度v,概率a与值一一对应。
对于如下的正态分布:
令则原正态分布就转变为标准化的正态分布进而可以转化为χ2分布。
假定一组放射性测量的变量ni(i=1,2,3,...,k)服从同一个高斯分布,期望值为m,标准误差为则服从自由度为k-1的χ2分布。
由于k个变量属于同一个分布,故只有k-1个变量是独立的,所以χ2分布的自由度为k-1。
χ2取值应与自由度为k-1的χ2分布相对应。
如果χ2取值与χ2分布不一致,说明ni服从同一个高斯分布的假设有问题。
χ2检验的具体步骤如下:
取一个小概率值α1=0.05,计算相应自由度下的χ2分布概率积分值,当测量值时,ni服从同一个高斯分布的假设正确;当时,ni服从同一个高斯分布的假设不成立。
因为χ2值出现在大于值的概率只有0.05,属于小概率事件。
同样取α2=0.95,计算对应的当测量值时,ni服从同一个高斯分布的假设不成立;当测量值ni服从同一个高斯分布的假设正确。
故只有当时,ni服从同一个高斯分布的假设才是正确的。
对放射性测井仪器测量的两组数据(X1,X2,…,Xn1)和(Y1,Y2,...,Yn2),分别求得对应的平均值和及标准偏差Sx和Sy,和中包含了系统偏离的信息,而Sx和Sy中包含了波动的信息。
可以通过T检验来判断两次测量的均值是否存在显著差异,通过F检验判断两次测量的标准差是否存在显著差异。
计算两组数据的Sx 2和Sy 2,用Sx 2和Sy 2的大的作分子,小的为分母,求出F值:
然后计算置信概率95%的F0.05值进行比较。
若F>F0.05,说明两组数据差异显著,说明两次测量时仪器波动很大,不必继续t检验;若F<F0.05,说明两组数据无明显差异,还须进一步用t检验法检验和之间是否存在显著性差异。
t检验法是在Sx和Sy无明显差异时,用于继续检验和之间有无显著差异的方法。具体作法是依下式计算t值:
再根据总自由度(f=n1+n2-2)计算置信概率95%的t0.05值。
比较t与t0.05,若t<t0.05,则认为和之间无明显差异。
若t>t0.05,则和有显著差异,两次测量存在系统偏差。
本申请实施例的检测放射性测井仪器稳定性的方法,可以校验放射性测井仪器在测井前后的状态是否稳定。如图1所示,本申请实施例的放射性测井仪器稳定性检测方法,主要包括如下内容。
步骤S111,在放射性测井仪器测井中,如果需要装源则装源,然后将放射性测井仪器放入井中,加电。如果放射源是中子,则采集中子计数。如果放射源是伽马源,则采集伽马计数。将仪器下放至套管中的第一深度,保持仪器贴紧刻度器的状态不动。其中,仪器可以水平置于刻度器内或者垂直贴紧刻度器。
步骤S112,采集p次中子计数或者伽马计数,获得第一采集结果N[1,p],其中p大于等于8。
步骤S113,利用X2检验方法判断第一采集结果N[1,p]是否符合正态分布,如果符合正态分布则转步骤S114;如果不符合正态分布,则放弃N[1,p]并转步骤S112继续采集下一组p个计数。如果多次检验结果均得出第一采集结果N[1,p]不符合正态分布,则可以认为仪器有问题,此时停止测量。
步骤S114,计算第一采集结果N[1,p]的均值、方差及变异系数。
步骤S115,开始测井,并在测井结束后,将仪器移动至套管中的第二深度,保持仪器不动。
步骤S116,采集p次中子计数或者伽马计数,获得第二采集结果N[2,p]。
步骤S117,利用X2检验方法判断第二采集结果N[2,p]是否符合正态分布,如果符合正态分布则转步骤S118,否则转步骤S116继续采集下一组p个计数N[2,p]。如果多次检验结果均得出第二采集结果N[2,p]不符合正态分布,则可以认为仪器有问题,此时停止测量。
步骤S118,计算第二采集结果N[2,p]的均值、方差及变异系数。
步骤S119,利用T检验和F检验判断第一采集结果N[1,p]的均值、方差及变异系数和第二采集结果N[2,p]的均值、方差及变异系数有无统计上的显著差异。结束测井。
本申请的实施例中,放射性测井仪器位于套管的第一深度时,采集放射源的p次计数至少会进行一轮。同样地,放射性测井仪器位于套管的第二深度时,采集放射源的p次计数也至少会进行一轮。
本申请的实施例中,如果第一采集结果N[1,p]的均值和第二采集结果N[2,p]的均值有统计上的显著差异,或者第一采集结果N[1,p]的方差和第二采集结果N[2,p]的方差有统计上的显著差异,或者第一采集结果N[1,p]的变异系数和第二采集结果N[2,p]的变异系数有统计上的显著差异,则说明放射性测井仪器在测井前后的状态不够稳定。如果第一采集结果N[1,p]及第二采集结果N[2,p]在均值、方差及变异系数这三方面均没有统计上的显著差异,则说明放射性测井仪器在测井前后的状态是稳定的。
本申请的实施例中,本领域的技术人员对于第一采集结果N[1,p]的均值、方差及变异系数和第二采集结果N[2,p]的均值、方差及变异系数有无统计上的显著差异,是可以理解和分辨的。
如图2所示,本申请实施例的放射性测井仪器刻度方法主要包括如下内容。
步骤S211,如果需要装源则装源,将放射性测井仪器放入刻度器,加电。如果放射源是中子,则采集中子计数。如果放射源是伽马源,则采集伽马计数。
步骤S212,保持仪器贴紧刻度器的状态不动,对中子计数或者伽马计数进行实时处理,利用X2检验方法判断p次测量的计数是否符合正态分布,如果是,则转至步骤S213,否则重复执行步骤S211。其中,仪器可以水平置于刻度器内或者垂直贴紧刻度器。
步骤S213,仪器放入第1块刻度器中开始刻度,采集p次计数,记为N[1,p],p≥8。
步骤S214,利用X2检验方法判断p次测量数据N[1,p]是否符合正态分布,如果是转至步骤S215,否则重复执行步骤S213。
步骤S215,保存N[1,p]数据,并根据N[1,p]数据进行刻度参数的计算,完成放射性测井仪器在第1块刻度器中的刻度。
上述刻度参数,是根据测量的计数率与刻度器的物理值对应关系,计算出来的仪器参数,或者在计算仪器参数的过程中所产生的中间值。
步骤S216,将仪器放入第i+1块刻度器中。其中,i大于等于1。采集获得p次计数,记为N[i+1,p]。
步骤S217,利用X2检验方法判断p次测量数据N[i+1,p]是否符合正态分布,如果是转至步骤S218,否则重复执行步骤S217。
步骤S218,保存N[i+1,p],并利用T检验和F检验判断N[1,p]与N[i+1,p]有无统计上的显著差异,如果没有统计上的显著差异则转至步骤S219,否则转至步骤S213,重新刻度。如果经过多次的循环,N[1,p]与N[i+1,p]总是有统计上的显著差异,则可以认为仪器不稳定,可以得出仪器不稳定的结论。
步骤S219,放射性测井仪器在第i+1块刻度器中的刻度结束。
其中,对应刻度是否结束,本领域的普通技术人员根据具体的应用,是完全知晓刻度在何时结束。比如对于密度的刻度,一共有四个步骤;这四个步骤完成,密度的刻度即完成。又比如对于孔隙度的刻度,是一个步骤;这一个步骤完成,孔隙度的刻度即完成。
本申请实施例的检测放射性测井仪器稳定性的方法,可以检验放射性测井仪器的质量。如图3所示,本申请实施例的放射性测井仪器稳定性检测方法,主要包括如下内容。
步骤S311,如果需要装源则装源,将放射性测井仪器放入井中或刻度器中,加电。如果放射源是中子,则采集中子计数;采集如果放射源是伽马源,则采集伽马计数;记录仪器开始预热时间。
步骤S312,保持仪器贴紧井壁或者刻度器的状态不动,对中子计数或者伽马计数进行实时处理,采集p次计数,记为测量数据N[1,p]。
步骤S313,利用X2检验方法判断p次测量数据N[1,p]是否符合正态分布。如果是,则转至步骤S314,否则返回执行步骤S312。
步骤S314,计算获得N[1,p]的均值、方差以及变异系数,保存N[1,p],并将其均值、方差以及变异系数作为仪器测试开始时的稳定性指标,记录此刻时间time1作为第一时刻。
步骤S315,采集第i+1次的p次计数,记为N[i+1,p];其中,i大于等于1。
步骤S316,利用X2检验方法确定N[1,p],N[2,p],...N[i+1,p]的(i+1)×p个数据符合正态分布。
步骤S317,在N[1,p],N[2,p],...N[i+1,p]这(i+1)×p个数据符合正态分布时,计算获得N[i+1,p]的均值、方差以及变异系数,保存N[i+1,p],并将N[i+1,p]的均值、方差以及变异系数作为仪器测试结束时的稳定性指标,记录此时时间time(i+1)作为第i+1时刻。
步骤S318,利用T检验和F检验判断N[1,p]与N[i+1,p]有无统计上的显著差异。具体地,判断N[1,p]与N[i+1,p]的均值是否具有统计上的显著差异,判断N[1,p]与N[i+1,p]的方差是否具有统计上的显著差异,以及判断N[1,p]与N[i+1,p]的变异系数是否具有统计上的显著差异。如果判断处N[1,p]与N[i+1,p]具有统计上的显著差异,则i=i+1,转步骤S315,否则转至步骤S319。
步骤S319,根据第一时刻以及第i+1时刻得到仪器最大稳定时间Tmax=time(i+1)-time1,结束稳定性测试。
本申请的实施例判断放射性测井仪器在测井前后以及刻度时是否状态稳定。本申请的实施例在仪器生产检验时,能够自动判断仪器稳定所需的时间,以及仪器能够稳定工作的最长时间。本申请的实施例能够给出仪器稳定性指标。
如图4所示,本申请实施例的放射性测井仪器稳定性检测装置,主要包括采集模块410、计算模块420以及判断模块430。所述放射性测井仪器中设置有放射源。
采集模块410,设置为所述放射性测井仪器位于套管中的第一深度时,采集所述放射源符合正态分布的p次计数,获得第一采集结果;所述放射性测井仪器位于所述套管中的第二深度时,采集所述放射源符合正态分布的p次计数,获得第二采集结果;其中p大于等于8;
计算模块420,与采集模块410相连,设置为计算所述第一采集结果的均值、方差及变异系数;计算所述第二采集结果的均值、方差及变异系数;
判断模块430,与计算模块420相连,设置为根据所述第一采集结果的均值、方差及变异系数,以及所述第二采集结果的均值、方差及变异系数,确定所述放射性测井仪器的稳定性。
所述采集模块410设置为在所述第一深度时采集所述放射源的至少一轮p次计数,利用X2检验方法对所述第一深度时采集到的至少一轮p次计数是否符合正态分布进行判断,将其中符合正态分布的一轮p次计数作为所述第一采集结果;在所述第二深度时采集所述放射源的至少一轮p次计数,利用X2检验方法对所述第二深度时采集到的至少一轮p次计数是否符合正态分布进行判断,将其中符合正态分布的一轮p次计数作为所述第二采集结果。
如图5所示,本申请实施例的放射性测井仪器刻度装置,主要包括获取模块510以及刻度模块520。
获取模块510,设置为获得所述放射性测井仪器在多块刻度器中各自符合正态分布的p次测量数据;其中,p大于等于8;
刻度模块520,与获取模块510相连,设置为所述放射性测井仪器在第i+1块刻度器中获得的符合正态分布的p次测量数据N[i+1,p]与所述放射性测井仪器在第一块刻度器中获得的符合正态分布的p次测量数据N[1,p]没有统计上的显著差异时,根据所述测量数据N[i+1,p]进行刻度参数的计算,完成所述放射性测井仪器在所述第i+1块刻度器中的刻度;其中,i大于等于1。
如图6所示,本申请实施例的放射性测井仪器稳定性检测装置,主要包括采集模块610以及确定模块620。
采集模块610,设置为采集所述放射性测井仪器位于井中或刻度器中时符合正态分布的p次测量数据N[1,p]并保存,记录保存时刻为第一时刻;采集第i+1次所述放射性测井仪器位于井中或刻度器中时符合正态分布的p次测量数据N[i+1,p]并保存,记录保存时刻为第i+1时刻;其中,i大于等于1,p大于等于8。
确定模块620,与采集模块610相连,设置为在从N[1,p]到N[i+1,p]这(i+1)×p个数据符合正态分布,且所述p次测量数据N[i+1,p]与所述p次测量数据N[1,p]没有统计上的显著差异时,将所述第i+1时刻与所述第一时刻的差值作为所述放射性测井仪器的最大稳定时间。
采用本申请的实施例,无需测前测后校验,即可定量检测放射性测井仪器在测井前后的稳定性,简化了测井程序。将本申请的实施例应用于放射性测井仪器刻度时,可以实时检测仪器的稳定性,减少人为判断产生的刻度误差,提高放射性测井仪的刻度、测量准确度。本申请的实施例,在仪器质量检验中,可用来判断仪器的起始稳定时间和结束稳定时间,以及对应的稳定性指标,为定量判断仪器性能提高依据。
虽然本发明所揭露的实施方式如上,但所述的内容仅为便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属领域内的技术人员,在不脱离本发明所揭露的精神和范围的前提下,可以在实施的形式及细节上进行任何的修改与变化,但本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (9)
1.一种放射性测井仪器稳定性检测方法,该方法包括:
所述放射性测井仪器位于套管中的第一深度时,采集符合正态分布的p次计数,获得第一采集结果;其中p大于等于8;
计算所述第一采集结果的均值、方差及变异系数;
所述放射性测井仪器位于所述套管中的第二深度时,采集符合正态分布的p次计数,获得第二采集结果;
计算所述第二采集结果的均值、方差及变异系数;
根据所述第一采集结果的均值、方差及变异系数,以及所述第二采集结果的均值、方差及变异系数,确定所述放射性测井仪器的稳定性,
具体包括:利用T检验和F检验判断出所述第一采集结果的均值与所述第二采集结果的均值、所述第一采集结果的方差与所述第二采集结果的方差以及所述第一采集结果的变异系数与所述第二采集结果的变异系数没有统计上的显著差异时,确认所述放射性测井仪器具备稳定性。
2.根据权利要求1所述的方法,其中,
采集符合正态分布的p次计数,获得第一采集结果,包括:
在所述第一深度时采集至少一轮p次计数,利用χ2分布检验方法对所述第一深度时采集到的至少一轮p次计数是否符合正态分布进行判断,将其中符合正态分布的一轮p次计数作为所述第一采集结果;
采集符合正态分布的p次计数,获得第二采集结果,包括:
在所述第二深度时采集至少一轮p次计数,利用χ2分布检验方法对所述第二深度时采集到的至少一轮p次计数是否符合正态分布进行判断,将其中符合正态分布的一轮p次计数作为所述第二采集结果。
3.一种放射性测井仪器刻度方法,该方法包括:
获得所述放射性测井仪器在多块刻度器中各自符合正态分布的p次测量数据;其中,p大于等于8;
所述放射性测井仪器在第i+1块刻度器中获得的符合正态分布的p次测量数据N[i+1,p]与所述放射性测井仪器在第一块刻度器中获得的符合正态分布的p次测量数据N[1,p]没有统计上的显著差异时,根据所述测量数据N[i+1,p]进行刻度参数的计算,完成所述放射性测井仪器在所述第i+1块刻度器中的刻度;其中,i大于等于1。
4.根据权利要求3所述的方法,其中:
利用T检验和F检验判断所述放射性测井仪器在第i块刻度器中获得的符合正态分布的p次测量数据N[i+1,p]与所述放射性测井仪器在第一块刻度器中获得的符合正态分布的p次测量数据N[1,p]是否具有统计上的显著差异。
5.一种放射性测井仪器稳定性检测方法,该方法包括:
采集所述放射性测井仪器位于井中或刻度器中时符合正态分布的p次测量数据N[1,p]并保存,记录保存时刻为第一时刻;其中p大于等于8;
采集第i+1次所述放射性测井仪器位于井中或刻度器中时符合正态分布的p次测量数据N[i+1,p]并保存,记录保存时刻为第i+1时刻;其中,i大于等于1;
在从N[1,p]到N[i+1,p]这(i+1)×p个数据符合正态分布,且所述p次测量数据N[i+1,p]与所述p次测量数据N[1,p]没有统计上的显著差异时,将所述第i+1时刻与所述第一时刻的差值作为所述放射性测井仪器的最大稳定时间。
6.根据权利要求5所述的方法,其中:
利用T检验和F检验判断N[1,p]与N[i+1,p]是否具有统计上的显著差异。
7.一种放射性测井仪器稳定性检测装置,该装置包括:
采集模块,设置为所述放射性测井仪器位于套管中的第一深度时,采集符合正态分布的p次计数,获得第一采集结果;所述放射性测井仪器位于所述套管中的第二深度时,采集符合正态分布的p次计数,获得第二采集结果;其中p大于等于8;
计算模块,设置为计算所述第一采集结果的均值、方差及变异系数;计算所述第二采集结果的均值、方差及变异系数;
判断模块,设置为根据所述第一采集结果的均值、方差及变异系数,以及所述第二采集结果的均值、方差及变异系数,确定所述放射性测井仪器的稳定性,具体包括:利用T检验和F检验判断出所述第一采集结果的均值与所述第二采集结果的均值、所述第一采集结果的方差与所述第二采集结果的方差以及所述第一采集结果的变异系数与所述第二采集结果的变异系数没有统计上的显著差异时,确认所述放射性测井仪器具备稳定性。
8.一种放射性测井仪器刻度装置,该装置包括:
获取模块,设置为获得所述放射性测井仪器在多块刻度器中各自符合正态分布的p次测量数据;其中,p大于等于8;
刻度模块,设置为所述放射性测井仪器在第i+1块刻度器中获得的符合正态分布的p次测量数据N[i+1,p]与所述放射性测井仪器在第一块刻度器中获得的符合正态分布的p次测量数据N[1,p]没有统计上的显著差异时,根据所述测量数据N[i+1,p]进行刻度参数的计算,完成所述放射性测井仪器在所述第i+1块刻度器中的刻度;其中,i大于等于1。
9.一种放射性测井仪器稳定性检测装置,该装置包括:
采集模块,设置为采集所述放射性测井仪器位于井中或刻度器中时符合正态分布的p次测量数据N[1,p]并保存,记录保存时刻为第一时刻;采集第i+1次所述放射性测井仪器位于井中或刻度器中时符合正态分布的p次测量数据N[i+1,p]并保存,记录保存时刻为第i+1时刻;其中,i大于等于1,p大于等于8;
确定模块,设置为在从N[1,p]到N[i+1,p]这(i+1)×p个数据符合正态分布,且所述p次测量数据N[i+1,p]与所述p次测量数据N[1,p]没有统计上的显著差异时,将所述第i+1时刻与所述第一时刻的差值作为所述放射性测井仪器的最大稳定时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310126466.1A CN103235350B (zh) | 2013-04-12 | 2013-04-12 | 放射性测井仪器稳定性检测及刻度方法与装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310126466.1A CN103235350B (zh) | 2013-04-12 | 2013-04-12 | 放射性测井仪器稳定性检测及刻度方法与装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103235350A CN103235350A (zh) | 2013-08-07 |
CN103235350B true CN103235350B (zh) | 2016-01-27 |
Family
ID=48883401
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310126466.1A Active CN103235350B (zh) | 2013-04-12 | 2013-04-12 | 放射性测井仪器稳定性检测及刻度方法与装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103235350B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105158610B (zh) * | 2015-09-14 | 2018-12-11 | 广西电网有限责任公司电力科学研究院 | 一种变压器状态预警数据可疑值的筛查处理方法 |
CN106706017B (zh) * | 2016-12-06 | 2019-04-26 | 深圳怡化电脑股份有限公司 | 一种器件稳定性测试方法及装置 |
CN106813703B (zh) * | 2016-12-20 | 2019-06-28 | 歌尔科技有限公司 | 一种测试产品功能的方法和装置 |
CN108049865B (zh) * | 2017-09-12 | 2021-08-03 | 中国石油天然气集团公司 | 随钻感应测井电场刻度方法 |
CN110487315B (zh) * | 2019-08-28 | 2021-09-28 | 广东核电合营有限公司 | 一种仪表漂移的分析系统及方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5196698A (en) * | 1991-11-01 | 1993-03-23 | Baker Hughes Corporation | Method and apparatus for nuclear logging using lithium detector assemblies |
CA2820922A1 (en) * | 2010-12-10 | 2012-06-14 | Schlumberger Canada Limited | Method for estimation of borehole and formation properties from nuclear logging measurements |
CN102562047A (zh) * | 2010-12-31 | 2012-07-11 | 中国石油天然气集团公司 | 阵列感应仪器居中情况下的井眼环境校正系统及校正方法 |
-
2013
- 2013-04-12 CN CN201310126466.1A patent/CN103235350B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5196698A (en) * | 1991-11-01 | 1993-03-23 | Baker Hughes Corporation | Method and apparatus for nuclear logging using lithium detector assemblies |
CA2820922A1 (en) * | 2010-12-10 | 2012-06-14 | Schlumberger Canada Limited | Method for estimation of borehole and formation properties from nuclear logging measurements |
CN102562047A (zh) * | 2010-12-31 | 2012-07-11 | 中国石油天然气集团公司 | 阵列感应仪器居中情况下的井眼环境校正系统及校正方法 |
Non-Patent Citations (2)
Title |
---|
提高521补偿密度测井仪器刻度的准确性;王相彬;《石油仪器》;20121031;第26卷(第5期);第38-40页 * |
测井分析中的精度评价问题;侯哲国;《石油勘探与开发》;19811231(第6期);第34-43页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103235350A (zh) | 2013-08-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103235350B (zh) | 放射性测井仪器稳定性检测及刻度方法与装置 | |
CN102608649B (zh) | 统计分布式γ或X射线能谱解谱方法 | |
CN104570047B (zh) | 伽马能谱测井仪自稳谱装置和方法 | |
CN106782710B (zh) | 一种反应堆测瞬发中子衰减参数的数据处理方法 | |
CN103513287B (zh) | 一种利用直流可控中子源计算地层密度的测井方法 | |
CN102084271B (zh) | 由核谱学确定绝对元素浓度 | |
CN107229080B (zh) | 一种元素测井中子俘获伽马能谱的获取方法 | |
CN103678778A (zh) | 一种放射性物化探信息综合方法 | |
CN106597560B (zh) | 一种利用快中子场分布表征的中子伽马密度测井方法 | |
CN1151385C (zh) | 确定岩层中油饱和度的方法 | |
CN106855990A (zh) | 核电机组仪表通道测量误差论证方法 | |
CN109581468B (zh) | 一种环境本底下弱γ放射源识别方法 | |
CN107288629A (zh) | 一种基于新型n‑γ双粒子探测器的中子伽马密度测井方法 | |
CN115331852B (zh) | 一种次临界反应堆控制棒反应性价值测量方法 | |
CN115826070A (zh) | 利用岩心编录仪确定待探测的白岗岩型铀矿中具有开采价值位置的方法 | |
CN115012920A (zh) | 基于双clyc双粒子探测器的可控中子源多谱测井仪及测井方法 | |
CN108222927A (zh) | 一种基于x射线源的密度测井方法 | |
CN113466919A (zh) | 一种大气中氡气浓度检测系统及方法 | |
US5272629A (en) | Method for determining the slowing down length and the porosity of a formation surrounding a borehole | |
CN115614022B (zh) | 铀裂变瞬发中子测井的井眼水层原位修正方法 | |
CN105673006B (zh) | 一种用于确定剩余油饱和度的碳氧比值计算方法 | |
CN115248463B (zh) | 一种d-t源中子孔隙度测井非弹性散射影响的校正方法 | |
CN112016215B (zh) | 一种基于牛顿迭代法的密度测井六参数反演计算方法 | |
CN113902348B (zh) | 一种基于孔内伽玛测试快速评价隧道工程放射性风险等级的方法 | |
RU2754993C1 (ru) | Способ калибровки счетного канала реактиметра |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CP01 | Change in the name or title of a patent holder |
Address after: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing Co-patentee after: China Oilfield Services Limited Patentee after: China Offshore Oil Group Co., Ltd. Address before: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing Co-patentee before: China Oilfield Services Limited Patentee before: China National Offshore Oil Corporation |