[發(fā)明專利]核電廠點(diǎn)源面源組合的復(fù)合輻射源強(qiáng)逆推方法及系統(tǒng)

有效

• 申請(qǐng)?zhí)枺?/td>: 201610038853.3
• 申請(qǐng)日: 2016-01-20
• 公開（公告）號(hào): CN106991621B
• 公開（公告）日: 2019-10-29
• 發(fā)明（設(shè)計(jì)）人: 陳義學(xué);賀淑相;臧啟勇;張涵;王夢(mèng)琪
• 申請(qǐng)（專利權(quán)）人: 華北電力大學(xué)
• 主分類號(hào): G06Q50/06
• 分類號(hào): G06Q50/06
• 代理公司: 北京康思博達(dá)知識(shí)產(chǎn)權(quán)代理事務(wù)所(普通合伙) 11426
• 代理人: 劉冬梅;路永斌

1.一種核電廠點(diǎn)源面源組合的復(fù)合輻射源強(qiáng)逆推方法，其特征在于，該方法包括如下步驟：

步驟一，用探測(cè)器探測(cè)核電廠內(nèi)的劑量率D₁,D₂,D₃…D_i；

步驟二，根據(jù)探測(cè)到的劑量率信息，建立如下式(一)所示的含有輻射源強(qiáng)度的超定方程組，

其中，所述超定方程組的系數(shù)矩陣a_i,j通過下式(二)和(三)得到，

在式(三)中，當(dāng)輻射源為點(diǎn)源時(shí)，p＝0，M＝1，N＝1；當(dāng)輻射源為面源時(shí)，p＝1；

步驟三，通過最小二乘法處理步驟二中的超定方程組，得到如下式(四)所示輻射源強(qiáng)度信息，

S_j,0＝(a_j,i·a_i,j)^-1·a_j,i·D_i (四)；

其中，D_i表示第i個(gè)探測(cè)器探測(cè)得到的劑量率；j表示輻射源的個(gè)數(shù)；m表示輻射源個(gè)數(shù)能達(dá)到的最大值；S_j表示第j個(gè)輻射源的強(qiáng)度；S_j,0表示初始計(jì)算未進(jìn)行迭代的第j個(gè)輻射源的強(qiáng)度；a_i,j表示系數(shù)矩陣，是第j個(gè)輻射源對(duì)第i個(gè)探測(cè)器的劑量響應(yīng)系數(shù)；BD(E_l,L(E_l,r₀→r_p)表示積累因子，是E_l和L(E_l,r₀→r_p)的函數(shù)；L(E_l,r₀→r_p)表示光學(xué)距離，是E_l和r₀→r_p的函數(shù)；μ表示截面/線性衰減系數(shù)；r₀→r_p表示輻射源到探測(cè)點(diǎn)的距離；C(E_l)表示通量-劑量轉(zhuǎn)換因子，是E_l的函數(shù)；E_l表示能量，是核電廠中輻射源發(fā)出的伽瑪射線的平均能量；表示離散源強(qiáng)；M和N分別表示面源在二維坐標(biāo)上離散后兩個(gè)坐標(biāo)軸上的離散標(biāo)號(hào)；

在步驟三之后，所述方法還包括如下步驟，

步驟四，根據(jù)步驟三中得到的輻射源強(qiáng)度信息計(jì)算探測(cè)器位置處的劑量率，D′₁,D′₂,D′₃…D′_i；

步驟五，對(duì)探測(cè)器探測(cè)到的劑量率信息和計(jì)算得到的探測(cè)器位置處的劑量率信息進(jìn)行線性擬合，得到擬合后的兩者關(guān)系的線性方程，進(jìn)而得到擬合參數(shù)，所述擬合參數(shù)包括：平均不確定度、擬合優(yōu)度和對(duì)應(yīng)的權(quán)重矩陣；

步驟六，將步驟五中得到的權(quán)重矩陣迭代至步驟二中的超定方程組，得到加權(quán)的超定方程，進(jìn)而重復(fù)步驟二、步驟三和步驟四，直至獲得期望的輻射源強(qiáng)度信息；

其中，D′_i表示計(jì)算出的第i個(gè)探測(cè)器位置處的劑量率。