[發明專利]核電廠點線面體組合的復合輻射源強逆推方法及系統

有效

• 申請號: 201610224072.3
• 申請日: 2016-04-12
• 公開（公告）號: CN107290770B
• 公開（公告）日: 2019-12-24
• 發明（設計）人: 陳義學;賀淑相;臧啟勇;王夢琪;張涵
• 申請（專利權）人: 華北電力大學
• 主分類號: G01T1/16
• 分類號: G01T1/16;G01T1/02
• 代理公司: 11426 北京康思博達知識產權代理事務所(普通合伙)
• 代理人: 劉冬梅;路永斌

1.一種核電廠點源線源面源體源組合的復合輻射源強逆推方法，其特征在于，該方法包括如下步驟：

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

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

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

在式(三)中，當輻射源為點源時，p＝0，L＝1，M＝1，N＝1；當輻射源為線源時，p＝1，M＝1，N＝1；當輻射源為面源時，p＝1，L＝1；當輻射源為體源時，p＝1；

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

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

其中，D_i表示第i個探測器探測得到的劑量率；j表示輻射源的個數；m表示輻射源個數能達到的最大值；S_j表示第j個輻射源的強度；S_j,0表示初始計算未進行迭代的第j個輻射源的強度；a_i,j表示系數矩陣，是第j個輻射源對第i個探測器的劑量響應系數；BD(E,L(μ(E),r₀→r_p)表示積累因子，是E和L(μ(E),r₀→r_p)的函數；L(μ(E),r₀→r_p)表示光學距離，是μ(E)和r₀→r_p的函數；μ(E)表示截面/線性衰減系數；r₀→r_p表示輻射源到探測點的距離；C(E)表示通量-劑量轉換因子，是E的函數；E表示能量，是核電廠中輻射源發出的伽瑪射線的平均能量；表示離散源強，其中當輻射源為體源時L、M和N分別表示體源在三維坐標上離散后三個坐標軸上的離散標號，當輻射源為面源時M和N分別表示面源在二維坐標上離散后二個坐標軸上的離散，當輻射源為線源時L表示線源在一維坐標上離散后坐標軸上的離散標號；

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

步驟四，根據步驟三中得到的輻射源強度信息計算探測器位置處的劑量率，D′₁,D′₂,D′₃…D′_i；

步驟五，對探測器探測到的劑量率信息和計算得到的探測器位置處的劑量率信息進行線性擬合，得到擬合后的兩者關系的線性方程，進而得到擬合參數，所述擬合參數包括：平均不確定度、擬合優度和對應的權重矩陣；

步驟六，將步驟五中得到的權重矩陣迭代至步驟二中的超定方程組，得到加權的超定方程，進而重復步驟二、步驟三和步驟四，直至獲得期望的輻射源強度信息；

其中，D′_i表示計算出的第i個探測器位置處的劑量率。