[發(fā)明專(zhuān)利]一種用于核電設(shè)備可靠性Gamma-Poisson模型分布參數(shù)的估算方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201210127023.X
• 申請(qǐng)日: 2012-04-27
• 公開(kāi)（公告）號(hào): CN103106326A
• 公開(kāi)（公告）日: 2013-05-15
• 發(fā)明（設(shè)計(jì)）人: 何劼;張彬彬;何建東;劉海濱;胡軍濤;張琴芳
• 申請(qǐng)（專(zhuān)利權(quán)）人: 上海核工程研究設(shè)計(jì)院
• 主分類(lèi)號(hào): G06F19/00
• 分類(lèi)號(hào): G06F19/00
• 代理公司: 核工業(yè)專(zhuān)利中心 11007
• 代理人: 高尚梅

技術(shù)領(lǐng)域

本發(fā)明屬于核電廠概率安全評(píng)價(jià)技術(shù)領(lǐng)域，具體涉及一種用于核電設(shè)備可靠性Gamma-Poisson模型分布參數(shù)的估算方法。

背景技術(shù)

核電設(shè)備可靠性Gamma-Poisson模型通常用以描述核電設(shè)備的運(yùn)行失效，例如核電廠的泵在持續(xù)運(yùn)行一段時(shí)間內(nèi)發(fā)生故障的次數(shù)。在Gamma-Poisson模型中，Gamma(α，β)分布的分布參數(shù)α和β值往往很難估算。如果α和β值不夠準(zhǔn)確或者缺失，則核電廠概率安全評(píng)價(jià)(PSA)的準(zhǔn)確性將受到影響，而且日后無(wú)法對(duì)特定電廠的設(shè)備可靠性參數(shù)進(jìn)行貝葉斯更新。

發(fā)明內(nèi)容

本發(fā)明的目的在于提供一種用于核電設(shè)備可靠性Gamma-Poisson模型分布參數(shù)的估算方法，可以利用核電設(shè)備運(yùn)行的歷史數(shù)據(jù)，解決計(jì)算分布參數(shù)困難的問(wèn)題。

本發(fā)明的技術(shù)方案如下：一種用于核電設(shè)備可靠性Gamma-Poisson模型分布參數(shù)的估算方法，該方法具體包括如下步驟：

步驟1、采集核電設(shè)備的運(yùn)行數(shù)據(jù)，利用矩法估計(jì)獲得分布參數(shù)的初步估計(jì)值；

采集核電設(shè)備的運(yùn)行數(shù)據(jù)，獲取核電設(shè)備失效率λ的均值和方差利用矩法估計(jì)獲得分布參數(shù)α、β的初步估計(jì)值α₀和β₀；

步驟2、根據(jù)步驟1獲取的分布參數(shù)α的估計(jì)值α₀，以k為步長(zhǎng)，獲取分布參數(shù)α值的一組備選值α_j(j＝1，2，…)，利用似然方程獲取一組與備選值α_j(j＝1，2，…)相對(duì)應(yīng)的分布參數(shù)β值的一組備選值β_j(j＝1，2，…)；

步驟3、根據(jù)步驟2所獲得的分布參數(shù)備選值(α_j，β_j)，獲取最終的分布參數(shù)估計(jì)值；

將步驟2中所獲得的分布參數(shù)備選值(α_j，β_j)代入檢驗(yàn)函數(shù)f(α，β)中，計(jì)算獲得使檢驗(yàn)函數(shù)f(α，β)達(dá)到最小值f_m(α，β)時(shí)所對(duì)應(yīng)的(α_m，β_m)值，其中，

檢驗(yàn)函數(shù)f(α，β)為：

f(α,β)=|Σi=1n1β[(α+xi)tiβ+ti-xi]|]]>

式中，x_i為第i個(gè)樣本的失效次數(shù)；t_i，為第i個(gè)樣本的運(yùn)行時(shí)間；n為采集樣本容量；

檢驗(yàn)函數(shù)f(α，β)達(dá)到最小值f_m(α，β)時(shí)所對(duì)應(yīng)的(α_m，β_m)值即為所要獲得的分布參數(shù)α、β的最終估計(jì)值

所述的步驟3中檢驗(yàn)函數(shù)f(α，β)達(dá)到最小值f_m(α，β)時(shí)所對(duì)應(yīng)的(α_m，β_m)值中，若α_m等于步驟2中獲取一組備選值α_j(j＝1，2，…)范圍的邊界時(shí)，則以步驟1所獲得的分布參數(shù)的初步估計(jì)值α₀和β₀作為分布參數(shù)α、β的最終估計(jì)值。