• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于SVM的放射源快速定位技術(shù)研究

    2023-09-21 09:20:16李子錕樊海軍張濟(jì)鵬
    核技術(shù) 2023年9期
    關(guān)鍵詞:放射源能譜射線

    管 弦 魏 星 李子錕,2 樊海軍 張濟(jì)鵬 孫 濤

    1(國民核生化災(zāi)害防護(hù)國家重點實驗室 北京 102205)

    2(成都理工大學(xué) 地學(xué)核技術(shù)四川省重點實驗室 成都 610059)

    隨著放射源在能源、工業(yè)、農(nóng)業(yè)、醫(yī)療等領(lǐng)域的廣泛應(yīng)用,輻射安全要求也日益嚴(yán)格,截至2021年,國內(nèi)生產(chǎn)、銷售、使用放射性同位素和射線裝置的單位有92 453家,目前,在用放射源共156 539枚,當(dāng)年共發(fā)生5 起輻射事故,其中放射源丟失事故就占3起[1]。遺失放射源對社會與公眾的生命財產(chǎn)安全造成嚴(yán)重威脅,而常見放射源通常較小,發(fā)生丟失事故后只能使用輻射探測儀定位。

    傳統(tǒng)放射源定位方法是通過固定式監(jiān)測節(jié)點或者便攜式輻射測量儀遍歷待測區(qū)域,耗時長且效率低,需要專業(yè)人員進(jìn)入未知輻射水平區(qū)域內(nèi)操作,難以保障工作人員的生命健康安全。目前已有較多放射源快速定位相關(guān)研究,張振朝等[2]設(shè)計了一種由NaI、CsI、鍺酸鉍(Bi4Ge3O12,BGO)3種晶體與鉛耦合組成的γ 射線方向探測器,利用各晶體性能差異特點,快速確定γ射線入射方向。周楊等[3]設(shè)計了一套基于3 種NaI(Tl)晶體的γ 射線源定位陣列系統(tǒng),通過建立各探測器全能峰占比與入射光子角度間關(guān)系函數(shù)實現(xiàn)對放射源快速定位。楊靜遠(yuǎn)等[4]設(shè)計了多探測器放射源定位系統(tǒng),基于核電廠環(huán)境監(jiān)測數(shù)據(jù),建立指紋定位方法。Willis 等[5]設(shè)計一套四陣列式NaI(Tl)放射源定位裝置,估計入射射線角度。Durbin 等[6]使用機(jī)器學(xué)習(xí)算法GLAM,基于“十字形”陣列NaI(Tl)探測器對放射源進(jìn)行快速定位。Fragkos等[7]利用小型CZT探測網(wǎng)絡(luò),實現(xiàn)了對無屏蔽放射源的精準(zhǔn)定位。眾多研究成果表明,利用陣列探測器可實現(xiàn)放射源的快速、精準(zhǔn)定位。

    本文針對快速尋源定位問題,設(shè)計了一套NaI(Tl)探測陣列放射源定位裝置,利用蒙特卡羅方法,模擬4 個探測器計數(shù)率隨入射射線角度的變化曲線,基于支持向量機(jī)SVM(Support Vector Machine)回歸擬合獲得計數(shù)率-角度響應(yīng)模型,使用模擬數(shù)據(jù)初步驗證模型的有效性,并開展了放射源角度定位與位置定位物理實驗,進(jìn)一步驗證了探測裝置和定位模型可精準(zhǔn)、快速定位放射源位置。

    1 探測陣列定位原理

    探測陣列放射源定位裝置由4個10 cm×10 cm×10 cm 的NaI(Tl)閃爍探測器組成,裝置結(jié)構(gòu)如圖1所示。4個探測器呈“田”字型緊靠排列擺放在水平面上,探測器之間。在放射源定位過程中,源與探測裝置的距離通常較遠(yuǎn),遠(yuǎn)大于10 倍探測裝置尺寸(2 m),此時放射源對探測裝置所展開的空間角較小,因此,可近似認(rèn)為到達(dá)探測裝置表面的源γ射線為平行束[8-9]。方向為點源到探測裝置中心的連線向量,此時各探測器對源射線的響應(yīng)只與入射角度和能量相關(guān)。

    圖1 NaI(Tl)探測陣列放射源定位裝置結(jié)構(gòu)示意圖Fig.1 Structure diagram of array NaI (Tl) detector device for radiation source location

    為探究入射角度和探測器響應(yīng)之間的關(guān)系,利用蒙特卡羅模擬軟件SuperMC,開展平行束入射數(shù)值模擬,模擬模型如圖2 所示。以陣列探測器中心為原點建立直角坐標(biāo)系,以象限分別對探測器進(jìn)行編號。為提升采樣效率,降低模擬方差,源項設(shè)置為半徑0.5 m的圓形面源,面源中心與原點間的連線為源平面中心軸法線,連線長度為50 cm,取法線反方向與y軸正方向的夾角為入射射線參考角度θ,初級射線為發(fā)射角度(θ)的γ 平行束,模擬能量為0.661 MeV,發(fā)射粒子總數(shù)為109個,面源與探測器間為標(biāo)準(zhǔn)大氣壓空氣。圖3分別展示了θ分別為60°和210°時,γ射線在各探測器內(nèi)的碰撞分布。

    圖2 平行束入射蒙卡模擬模型示意圖Fig.2 Diagram of Monte Caro simulation model with parallel beam γ-ray

    圖3 γ射線在陣列探測器內(nèi)的碰撞分布 (a) θ=60°,(b) θ=210°Fig.3 Collision distribution of particles in array detector (a) θ = 60°,(b) θ = 210°

    從圖3 中可得出,平行束入射時,4 個探測器之間互相屏蔽,當(dāng)初級γ 射線入射角度發(fā)生變化后,γ射線在各探測器內(nèi)發(fā)生碰撞概率存在明顯變化??梢?,4 個探測器所測得能譜計數(shù)與入射射線角度存在一定關(guān)系,因此,模擬入射角度θ范圍[10°~360°],每隔10°取一發(fā)射角,4個探測器的全譜總計數(shù)隨角度分布如圖4所示。

    圖4 陣列探測器模擬總計數(shù)隨入射射線角度分布曲線Fig.4 Distribution curves of MC simulation gross counts with γ-ray incidence angles

    圖4縱坐標(biāo)為模擬各探測器發(fā)生能量沉積的光子數(shù),從中可得出,對所有探測器,其在固定角度范圍內(nèi)呈較為平滑的變化趨勢,4 個探測器的總計數(shù)與入射射線角度存在唯一對應(yīng)關(guān)系,且總計數(shù)與入射射線角度θ呈一定函數(shù)關(guān)系。但當(dāng)θ變化至90°、180°、270°、360°處時,探測器總計數(shù)變化趨勢發(fā)生明顯變化,這是由陰影效應(yīng)所導(dǎo)致,當(dāng)θ處于[90°,180°]范圍內(nèi),進(jìn)入1號探測器的部分初級γ射線被4號探測器所屏蔽,進(jìn)入4 號探測器的初級射線則沒有被屏蔽,而當(dāng)θ處于[0°,90°]范圍內(nèi)時則相反。致使在二維空間內(nèi)難以通過單個函數(shù)描述總計數(shù)和入射角度對應(yīng)關(guān)系,為解決這一問題,本文使用支持向量線性回歸方法對入射射線角度進(jìn)行預(yù)測。

    2 SVM線性回歸原理

    2.1 線性回歸理論基礎(chǔ)

    支持向量機(jī)(Support Vector Machine,SVM)是一種基于統(tǒng)計學(xué)習(xí)理論的前饋網(wǎng)絡(luò)學(xué)習(xí)方法,其通過核函數(shù),將數(shù)據(jù)從原始低維空間映射至高維歐式空間,建立超平面帶,在高維空間內(nèi)實現(xiàn)對數(shù)據(jù)線性回歸,從而解決非線性回歸問題[10]。

    假定現(xiàn)有一組測量數(shù)據(jù),包含各探測器總計數(shù)c與對應(yīng)入射射線角度θ,以改組數(shù)據(jù)為訓(xùn)練集,其如式(1)所示。

    式中:R為實數(shù)空間;R4為四維實數(shù)空間,對應(yīng)4個探測器計數(shù);l為測量點數(shù)量。

    從圖4 中可看出,此數(shù)據(jù)集無法使用線性函數(shù)描述輸入和輸出的對應(yīng)關(guān)系,因此,將原始數(shù)據(jù)映射到高維空間,具體如式(2)所示[11]。

    經(jīng)過映射Φ(xi)后,數(shù)據(jù)集TΦ在高維空間H內(nèi)為近似線性數(shù)據(jù)。

    給定高維空間的H一個超平面,其表達(dá)式為[11-12]:

    式中:ω為權(quán)重向量;b為偏置。將該超平面沿Y軸方向上下平移ε形成一個超平面帶,使得所有的數(shù)據(jù)點都被包含在內(nèi),將非線性回歸問題轉(zhuǎn)化為凸優(yōu)化問題,如式(4)所示[11-12]。

    式中:加和為ε-不敏感訓(xùn)練誤差的總和;C為懲罰因子,用于權(quán)衡訓(xùn)練誤差和懲罰項‖ω‖2;ξi和ξi' 為兩個非負(fù)松弛變量,其代表允許數(shù)據(jù)點越過超平面帶的程度[12-13]。

    為求解該最優(yōu)化問題,引入拉格朗日函數(shù),如式(5)所示[12-13]。

    式(5)中αi、a'i、ξi、ξ'i為拉格朗日乘子,將上式轉(zhuǎn)化為對偶問題,如式(6)所示[13]。

    在[0,C]范圍內(nèi)尋得最優(yōu)解,則凸二次規(guī)劃問題的最優(yōu)解如式(7)所示[14]。

    角度回歸預(yù)測函數(shù)如式(8)所示。

    為降低運算維度,使用高斯核函數(shù)代替映射內(nèi)積,則角度回歸預(yù)測函數(shù)如下所示[14-15]。

    2.2 入射角度線性回歸預(yù)測模型

    使用入射角度在[10°,180°]的模擬各探測器總計數(shù)作為訓(xùn)練集,采用高斯核函數(shù),使用留一驗證法獲得最佳核函數(shù)參數(shù),建立線性回歸模型。

    采用各探測器能譜計數(shù)占比作為預(yù)測模型輸入,真實角度作為輸出,如式(10)所示:

    式(10)中,cji為第i個測量點的j號探測器能譜計數(shù)。由于陣列探測器的尺寸相對于離源距離較小,可認(rèn)為在源與探測器的距離遠(yuǎn)大于探測器尺寸時,各探測器對放射源的張角以及距離一致,此時能譜計數(shù)占比只與入射射線角度相關(guān)。則計數(shù)占比誤差為:

    式中:ΔN為能譜計數(shù)標(biāo)準(zhǔn)誤差。

    將計數(shù)占比與真實角度作為訓(xùn)練集,其表達(dá)式如式(12):

    式中:n為訓(xùn)練集內(nèi)的測量點數(shù)量;Tx為訓(xùn)練集輸入,其是由4個探測器能譜計數(shù)之比構(gòu)建的矩陣;Ty為訓(xùn)練集輸出,其是由各測量點對應(yīng)真實角度組成的列向量。

    建立模型過程中,從訓(xùn)練集內(nèi)隨機(jī)抽取一個測量點的計數(shù)占比pi和真實角度θi為測試集,用于測試預(yù)測模型的輸出與真實輸出Ty間的誤差,多次迭代直至達(dá)到誤差最小,模型公式如式(13)所示:

    式中:Tx,i為第i個測量點4個探測器能譜計數(shù)之比組成的列向量;θ為預(yù)測模型輸出角度;μ為高斯核函數(shù)參考向量;σ為高斯核函數(shù)參數(shù);αˉ為最優(yōu)權(quán)重向量;bˉ為最優(yōu)偏置。μ、σ、αˉ、bˉ均通過求解凸優(yōu)化問題得到[9]。通過此回歸模型即可對射線入射角度進(jìn)行預(yù)測。

    使用蒙特卡羅模擬[10°~180°]范圍內(nèi)的模擬各探測器總計數(shù)能譜計數(shù)之比作為訓(xùn)練集獲得模型最優(yōu)參數(shù),在[190°~360°]范圍內(nèi)的模擬各探測器能譜計數(shù)之比作為測試集,測試入射角度預(yù)測模型,測試結(jié)果如表1所示。

    表1 MC模擬放射源角度預(yù)測結(jié)果(190°~360°)Table 1 Angle predictions of MC simulation radioactive source

    從表1 中可得出,基于SVM 角度回歸預(yù)測模型,最大角度絕對偏差為2.53°,平均絕對偏差為0.52°,對于模擬數(shù)據(jù)可精確估計入射射線角度。

    3 放射源定位物理實驗

    為驗證角度回歸預(yù)測模型定位效果,開展放射源定位物理實驗,探測裝置由4 個10 cm×10 cm×10 cm的NaI(Tl)探測器陣列組合而成。探測器安裝于一旋轉(zhuǎn)平臺之上,平臺可在水平方向360°旋轉(zhuǎn),放射源固定在一線性軌道之上,軌道長10 m,可遠(yuǎn)程操控放射源至軌道任意位置。實驗裝置實物如圖5所示。

    圖5 探測器裝置實物圖(a) NaI(Tl)陣列探測器,(b) 陣列旋轉(zhuǎn)臺Fig.5 Physical photos of array NaI(Tl) detector radiation source location device(a) Array NaI(Tl) detector, (b) Array rotation platform

    圖6 線性軌道實物圖(a) 軌道整體,(b) 放射源托架平臺Fig.6 Physical diagram of linear orbits(a) Orbits, (b) Radioactive source bracket platform

    3.1 放射源角度響應(yīng)實驗

    為獲取陣列探測器對入射射線角度的響應(yīng)關(guān)系,開展放射源角度響應(yīng)實驗,源位于軌道中心5 m處,探測器與放射源的連線與軌道垂直,為確保入射射線為平行束,兩者直線距離為4 m,調(diào)整探測器與源在同一水平面高度,通過旋轉(zhuǎn)平臺改變陣列探測器與放射源的相對角度。實驗裝置擺放以及測量過程如圖7所示。

    實驗過程中,記錄4個探測器的總計數(shù),單次測量時間為3 min,每次測量完畢后將平臺沿逆時針方向旋轉(zhuǎn)10°,圖7 展示了θ=0°時各探測器位置分布。實驗結(jié)果如圖8所示。

    圖8 放射源角度響應(yīng)實驗結(jié)果Fig.8 Experiment results of radioactive angle response

    從圖8 中可得出,實驗測得各探測器總計數(shù)隨入射射線角度的變化趨勢與模擬結(jié)果基本一致,然而,由于實驗測量總計數(shù)受放射源活度和放射源離探測器距離影響,難以直接通過探測器能譜總計數(shù)建立線性回歸預(yù)測模型,考慮到輻射事件的獨立性,對于固定能量的γ 平行束,各探測器總計數(shù)與所有探測器凈計數(shù)和之比與射入探測器的γ射線注量無關(guān),因此實驗采用的訓(xùn)練集n如式(14):

    式中:Ni為第i號探測器的能譜總計數(shù),則ni的誤差可表示為:

    訓(xùn)練集n隨入射射線角度分布如圖9所示。

    圖9 角度響應(yīng)實驗?zāi)茏V計數(shù)占比Fig.9 Proportions of gross count from energy spectra of radioactive angle response experiment

    同樣將[10°~180°]范圍內(nèi)的實驗測得各探測器能譜計數(shù)之比和真實角度分別作為訓(xùn)練集輸入和輸出,獲得模型最優(yōu)參數(shù),并建立角度預(yù)測回歸模型,在[190°~360°]范圍內(nèi)的實驗測得各探測器能譜計數(shù)之比作為測試集作為測試集輸入獲得預(yù)測角度,結(jié)果如表2所示。

    表2 角度響應(yīng)實驗角度預(yù)測結(jié)果Table 2 Prediction of angle response experiment results

    從表2 中可知,角度預(yù)測結(jié)果最大偏差出現(xiàn)在360°處,為6.5°,最小偏差出現(xiàn)在300°處,為0.18°,平均偏差為3.35°,表明在源與探測器距離固定情況下,預(yù)測模型能較為準(zhǔn)確預(yù)測入射射線角度。

    3.2 放射源位置響應(yīng)實驗

    為驗證角度預(yù)測回歸模型是否適用于不同距離下的場景,開展放射源位置響應(yīng)實驗,如圖10所示,實驗過程中,以軌道最左端為零點,遠(yuǎn)程操控137Cs點源從軌道左端起點(20 cm)處,移動至軌道右端終點(940 cm)處至軌道中心處,每隔40 cm 取一坐標(biāo)點進(jìn)行測量,單次測量時間為3 min,探測器與放射源平行,記錄4個探測器的全譜計數(shù),實驗過程中分別改變陣列探測器離軌道距離S,實驗結(jié)果如圖11所示。

    圖10 放射源位置響應(yīng)實驗示意圖Fig.10 Diagram of radioactive location response experiment

    將位置響應(yīng)實驗結(jié)果依照式(14)構(gòu)建測試集輸入,通過模型獲得預(yù)測角度,各距離實驗的角度預(yù)測結(jié)果如圖12。

    圖12 位置響應(yīng)實驗角度預(yù)測平均偏差Fig.12 Predicted mean angle deviation of location response experiment

    從圖12 中可得出,隨著探測器離軌道距離增大,角度預(yù)測平均偏差也隨之降低,這是由于當(dāng)源距離探測器較近時,點源對探測器的張角較大,此時入射射線不可視為平行束,因此誤差較大,當(dāng)源與探測器距離相對較遠(yuǎn)時(≥3 m),角度預(yù)測平均誤差低于3°,表明平行束入射時,模型能精確預(yù)測入射射線角度。

    4 定位影響因素測試

    4.1 相對高度影響測試

    在實際放射源定位場景中,放射源與定位探測器往往不處于同一水平面,為評估探測器與放射源相對高度差對定位效果的影響,開展探測陣列高度響應(yīng)實驗,實驗過程中,探測陣列與軌道的直線距離調(diào)整為4 m,探測陣列中軸線與軌道平行,放射源移動過程與位置響應(yīng)實驗一致,單次測量時間為3 min,通過操縱陣列旋轉(zhuǎn)臺,將探測陣列相對軌道的高度分別調(diào)至0.3 m、0.4 m、0.5 m、0.6 m,記錄各探測器的全譜計數(shù),實驗結(jié)果如圖13所示。

    圖13 高度響應(yīng)實驗結(jié)果 (a) 相對高度0.3 m,(b) 相對高度0.4 m,(c) 相對高度0.5 m,(d) 相對高度0.6 mFig.13 Results of relative height response experiment(a) Relative height 0.3 m, (b) Relative height 0.4 m, (c) Relative height 0.5 m, (d) Relative height 0.6 m

    從圖13 和表3 中可得出,不同高度下的探測器全譜計數(shù)響應(yīng)曲線未發(fā)生明顯變化,隨著探測器高度提升,雖然探測器與點源距離增大,但各探測器能譜計數(shù)反而提升,這可能由于γ 射線在垂直方向上不同角度入射方形探測器時,探測器對射線束的探測截面積也隨之增大,接收到初級射線的數(shù)量發(fā)生改變,致使對點源探測效率發(fā)生改變,導(dǎo)致隨探測器高度增高,能譜計數(shù)增大。使用角度預(yù)測模型對實驗數(shù)據(jù)進(jìn)行處理,定位平均偏差如圖14所示。

    表3 高度響應(yīng)實驗結(jié)果Table 3 Results of relative height response experiment

    圖14 高度響應(yīng)實驗角度預(yù)測平均偏差Fig.14 Predicted mean angle deviation of relative height response experiment

    從圖14中可得出,探測陣列與放射源相對高度為0 m 時定位偏差最小,隨著相對高度增加偏差逐漸增大,但變化幅度較小,偏差均小于5°,表明在一定相對高度范圍內(nèi),角度預(yù)測模型能夠較為精確地對放射源進(jìn)行定位。

    4.2 放射源能量與強(qiáng)度影響測試

    在實際定位場景中,考慮到探測器不同活度和不同能量的放射源的探測能力存在一定差異,定位能力必然存在一定不同,為探究對不同能量射線源的定位能力,開展蒙特卡羅模擬,模擬模型與圖1中的模型基本相同,分別設(shè)置射線能量為:50 keV、100 keV、200 keV、300 keV、500 keV、750 keV、1 000 eV、1 500 keV、2 000 keV、3 000 keV,模擬粒子數(shù)為109個,得到各探測器的歸一化全譜計數(shù)。使用角度預(yù)測模型定位放射源,結(jié)果如圖15所示。

    圖15 各能量射線入射角度預(yù)測平均偏差Fig.15 Predicted mean angle deviation with various energy γ-ray

    從圖15 中可得出,低能射線源(50~1 000 keV)的定位角度偏差相對較低,最大值為0.41°,高能射線源(1 500~3 000 keV)的定位角度偏差相對較高,最大值為0.77°,角度偏差整體隨能量而增大,這是由于兩個方面的原因所導(dǎo)致:一是當(dāng)射線能量越高,探測器的探測效率越低;二是當(dāng)射線能量高于1.022 MeV 后,發(fā)生電子對效應(yīng)產(chǎn)生的逃逸光子有可能進(jìn)入相鄰探測并發(fā)生能量沉積,致使響應(yīng)規(guī)律發(fā)生變化,但由于本文設(shè)計的定位裝置晶體體積較大,因此這兩個因素對定位結(jié)果影響較小,整體定位角度偏差較小。

    在以上模擬的基礎(chǔ)上,探究探測陣列與定位模型可定位不同能量射線源的最小活度。物理實驗過程中探測器的本底計數(shù)率約為1 600,重構(gòu)后不同活度下探測器全譜計數(shù)如式(16)所示:

    式中:N為探測器全譜計數(shù);I為放射源射線強(qiáng)度,為模擬歸一化全譜計數(shù);Poisson(1 600)為服從泊松分布的本底計數(shù)。設(shè)置射線強(qiáng)度范圍為[104,1010],使用不同強(qiáng)度下的探測器全譜計數(shù)作為輸入,獲得定位偏差變化規(guī)律。圖16 為能量為50 keV 時定位偏差隨射線強(qiáng)度的變化曲線。

    圖16 定位平均偏差隨γ射線源強(qiáng)度變化Fig.16 Varition of mean localization deviation with the primary γ-ray intensity

    從圖16 中可得出,隨著射線強(qiáng)度的增大,定位偏差也逐漸降低,定義定位偏差為5°時的射線強(qiáng)度為最小可定位強(qiáng)度,不同能量的最小可定位強(qiáng)度如圖17所示。

    圖17 放射源最小可定位強(qiáng)度Fig.17 Minimum localizable intensity of radiation source

    在探測陣列距離放射源4 m 時,源強(qiáng)度達(dá)到最小可定位強(qiáng)度即可確保射入探測陣列的光子所產(chǎn)生的全譜計數(shù)能夠?qū)崿F(xiàn)精準(zhǔn)定位,而探測陣列表面的射線強(qiáng)度與源離探測陣列的距離的平方成反比。因此,基于圖17中最小可定位強(qiáng)度可獲得不同距離下的最小可定位強(qiáng)度。

    5 結(jié)語

    針對丟失放射源快速定位問題,提出一種基于SVM的放射源角度回歸預(yù)測模型,該模型通過的各探測器能譜計數(shù)之比,快速預(yù)測放射源與探測器間的相對角度。放射源角度響應(yīng)實驗角度預(yù)測結(jié)果顯示,模型預(yù)測角度最大偏差為6.5°,最小偏差為0.18°,平均偏差為3.35°,表明在放射源與探測器距離固定的場景下,回歸模型能夠準(zhǔn)確預(yù)測相對角度;放射源位置響應(yīng)實驗角度預(yù)測結(jié)果顯示,隨著探測器與放射源的距離增大,平均偏差逐漸降低,最大偏差為9.21°,最小偏差為1.77°,說明當(dāng)入射射線越趨近于平行束時,回歸模型預(yù)測精度越高;通過相對高度測試實驗,證明定位模型在一定相對高度差內(nèi)均可實現(xiàn)對放射源的精準(zhǔn)定位;基于不同能量射線入射蒙卡模擬結(jié)果,獲得了不同能量下的最小可定位放射源強(qiáng)度,該方法可快速、精確預(yù)測放射源與陣列探測器的相對角度。

    作者貢獻(xiàn)聲明管弦負(fù)責(zé)算法設(shè)計、實驗數(shù)據(jù)收集與處理、文章的起草與修訂;孫濤負(fù)責(zé)研究的提出與設(shè)計和文章的最終修訂;魏星負(fù)責(zé)實驗方案設(shè)計;李子錕負(fù)責(zé)開展實驗與蒙特卡羅模擬;樊海軍負(fù)責(zé)資料調(diào)研;張濟(jì)鵬負(fù)責(zé)算法實現(xiàn)與測試。

    猜你喜歡
    放射源能譜射線
    寧夏銥-192放射源輻射事故調(diào)查及分析
    核安全(2022年3期)2022-06-29 09:17:50
    一起銫-137放射源失控事故應(yīng)急監(jiān)測探討
    核安全(2022年3期)2022-06-29 09:17:34
    能譜CT在術(shù)前預(yù)測胰腺癌淋巴結(jié)轉(zhuǎn)移的價值
    “直線、射線、線段”檢測題
    『直線、射線、線段』檢測題
    放射源在線監(jiān)控技術(shù)在醫(yī)院的應(yīng)用分析
    電子測試(2018年10期)2018-06-26 05:54:18
    赤石脂X-射線衍射指紋圖譜
    中成藥(2017年3期)2017-05-17 06:09:16
    M87的多波段輻射過程及其能譜擬合
    電子材料分析中的能譜干擾峰
    準(zhǔn)單色X射線機(jī)替代241Am放射源的測厚應(yīng)用研究
    同位素(2014年2期)2014-04-16 04:57:21
    国产精品人妻久久久久久| 亚洲欧美成人精品一区二区| 内射极品少妇av片p| videos熟女内射| 免费无遮挡裸体视频| 不卡视频在线观看欧美| 高清av免费在线| 中文精品一卡2卡3卡4更新| 久久精品影院6| 亚洲综合精品二区| 亚洲高清免费不卡视频| 免费在线观看成人毛片| 中文资源天堂在线| 亚洲欧美精品专区久久| 美女xxoo啪啪120秒动态图| 国产人妻一区二区三区在| 一级爰片在线观看| 日韩欧美 国产精品| 在线免费观看的www视频| 亚洲av电影在线观看一区二区三区 | 亚洲性久久影院| 久久久久免费精品人妻一区二区| 黄色一级大片看看| 在线天堂最新版资源| 永久免费av网站大全| 国内少妇人妻偷人精品xxx网站| 婷婷色麻豆天堂久久 | 桃色一区二区三区在线观看| 中文字幕亚洲精品专区| 啦啦啦韩国在线观看视频| 精品久久久久久久人妻蜜臀av| 国产精品乱码一区二三区的特点| 精品国产一区二区三区久久久樱花 | 秋霞伦理黄片| 18+在线观看网站| h日本视频在线播放| 久久精品人妻少妇| 最新中文字幕久久久久| 国产精品,欧美在线| 久久久精品欧美日韩精品| 久久久久精品久久久久真实原创| 国产精品国产高清国产av| 大香蕉久久网| 插逼视频在线观看| 国产一区二区三区av在线| 亚洲伊人久久精品综合 | av又黄又爽大尺度在线免费看 | 男女边吃奶边做爰视频| 免费看美女性在线毛片视频| 欧美成人一区二区免费高清观看| 老司机影院成人| 国产黄色视频一区二区在线观看 | 国产精品国产三级专区第一集| 身体一侧抽搐| 自拍偷自拍亚洲精品老妇| 久久99精品国语久久久| 色视频www国产| 欧美丝袜亚洲另类| 在线免费十八禁| 欧美3d第一页| 亚洲精品一区蜜桃| 赤兔流量卡办理| 国产毛片a区久久久久| av天堂中文字幕网| 成年女人看的毛片在线观看| 日韩人妻高清精品专区| 亚洲欧美成人综合另类久久久 | 免费看日本二区| 久久久a久久爽久久v久久| 日韩在线高清观看一区二区三区| 人体艺术视频欧美日本| 国产精品无大码| 天堂√8在线中文| 18禁在线无遮挡免费观看视频| 亚洲精品一区蜜桃| 日本黄大片高清| АⅤ资源中文在线天堂| a级毛色黄片| 伦精品一区二区三区| 国产精品无大码| 两个人视频免费观看高清| 99热这里只有精品一区| 免费观看在线日韩| 色综合站精品国产| 久久久久性生活片| 91狼人影院| 亚洲成人精品中文字幕电影| av国产久精品久网站免费入址| 天堂影院成人在线观看| 国产免费一级a男人的天堂| 亚洲欧洲日产国产| 亚洲av不卡在线观看| 精品酒店卫生间| 亚洲av男天堂| 成人午夜高清在线视频| 欧美zozozo另类| 久久国产乱子免费精品| 亚洲怡红院男人天堂| 国国产精品蜜臀av免费| 两个人视频免费观看高清| 麻豆成人午夜福利视频| 99久久人妻综合| 91狼人影院| 寂寞人妻少妇视频99o| www.色视频.com| 一级二级三级毛片免费看| 大香蕉97超碰在线| 晚上一个人看的免费电影| 可以在线观看毛片的网站| 国产高清视频在线观看网站| 免费看美女性在线毛片视频| 久久国内精品自在自线图片| 女的被弄到高潮叫床怎么办| 夫妻性生交免费视频一级片| 国产国拍精品亚洲av在线观看| 水蜜桃什么品种好| 欧美激情久久久久久爽电影| 国产精品,欧美在线| 久久久国产成人免费| 国产伦精品一区二区三区四那| 久久欧美精品欧美久久欧美| 丝袜美腿在线中文| 国产成人a区在线观看| 天天躁夜夜躁狠狠久久av| 欧美一区二区亚洲| 亚洲天堂国产精品一区在线| 久久久午夜欧美精品| 九九热线精品视视频播放| 男人的好看免费观看在线视频| 中文欧美无线码| 久久久久久伊人网av| 亚洲高清免费不卡视频| 18禁动态无遮挡网站| 国产伦在线观看视频一区| 国产乱人视频| 久久久久久国产a免费观看| 亚洲最大成人手机在线| 国产精品久久电影中文字幕| 看片在线看免费视频| 亚洲经典国产精华液单| 国产极品天堂在线| 久久精品影院6| 久久午夜福利片| 国产91av在线免费观看| 中国美白少妇内射xxxbb| 欧美精品一区二区大全| 国产欧美另类精品又又久久亚洲欧美| 一二三四中文在线观看免费高清| 久久精品夜色国产| 午夜精品一区二区三区免费看| www.av在线官网国产| 亚洲真实伦在线观看| 在线免费观看的www视频| 少妇猛男粗大的猛烈进出视频 | 最近视频中文字幕2019在线8| 免费看日本二区| 国产高清不卡午夜福利| 国国产精品蜜臀av免费| 亚洲国产精品国产精品| 日产精品乱码卡一卡2卡三| 两个人视频免费观看高清| 乱码一卡2卡4卡精品| 一区二区三区高清视频在线| 91在线精品国自产拍蜜月| 欧美高清性xxxxhd video| 伊人久久精品亚洲午夜| 人人妻人人澡人人爽人人夜夜 | 在线观看美女被高潮喷水网站| 男插女下体视频免费在线播放| 国产91av在线免费观看| 高清视频免费观看一区二区 | 日韩一本色道免费dvd| 久久久久久伊人网av| 99国产精品一区二区蜜桃av| 欧美成人免费av一区二区三区| 尾随美女入室| 久久精品人妻少妇| 欧美精品国产亚洲| 在线免费十八禁| 国产精品麻豆人妻色哟哟久久 | 国内少妇人妻偷人精品xxx网站| 高清av免费在线| 国产精品熟女久久久久浪| 中文在线观看免费www的网站| 人体艺术视频欧美日本| 国产高清视频在线观看网站| 一卡2卡三卡四卡精品乱码亚洲| 亚洲第一区二区三区不卡| 我要看日韩黄色一级片| 内地一区二区视频在线| 亚洲激情五月婷婷啪啪| 午夜老司机福利剧场| 国产午夜精品一二区理论片| 日韩三级伦理在线观看| 秋霞在线观看毛片| 永久网站在线| 国产免费视频播放在线视频 | 亚洲精品亚洲一区二区| 亚洲va在线va天堂va国产| 人人妻人人看人人澡| 免费播放大片免费观看视频在线观看 | 一卡2卡三卡四卡精品乱码亚洲| 亚洲中文字幕一区二区三区有码在线看| 久久久久精品久久久久真实原创| 秋霞在线观看毛片| 日本黄色视频三级网站网址| 免费大片18禁| 亚洲av中文字字幕乱码综合| 我的女老师完整版在线观看| 中文字幕免费在线视频6| 国产女主播在线喷水免费视频网站 | 亚洲精品影视一区二区三区av| 欧美潮喷喷水| 日韩欧美三级三区| 亚洲精品乱码久久久v下载方式| 国产乱人视频| 国产精品熟女久久久久浪| 国产探花极品一区二区| 久久精品综合一区二区三区| 国产精品蜜桃在线观看| 精品人妻偷拍中文字幕| 插阴视频在线观看视频| 黄色配什么色好看| 久久综合国产亚洲精品| 国产欧美日韩精品一区二区| 99热网站在线观看| 精品99又大又爽又粗少妇毛片| 日本免费一区二区三区高清不卡| 插阴视频在线观看视频| 久久久久久久久久久免费av| 一级二级三级毛片免费看| 真实男女啪啪啪动态图| 色噜噜av男人的天堂激情| 精品人妻熟女av久视频| 欧美日韩综合久久久久久| 偷拍熟女少妇极品色| 久久久精品大字幕| 亚洲怡红院男人天堂| 久久久久性生活片| 日日撸夜夜添| 国产亚洲精品久久久com| 欧美人与善性xxx| 久久人人爽人人爽人人片va| 日本三级黄在线观看| av在线观看视频网站免费| 亚洲欧美清纯卡通| 国产精品一区二区在线观看99 | 国产精品久久视频播放| 亚洲国产精品合色在线| 国产亚洲最大av| 天堂影院成人在线观看| 美女大奶头视频| 国产人妻一区二区三区在| 亚洲成色77777| 美女xxoo啪啪120秒动态图| 黄色欧美视频在线观看| 日日撸夜夜添| 99久国产av精品| 男女啪啪激烈高潮av片| 亚洲国产色片| 在线a可以看的网站| 一级二级三级毛片免费看| 99热这里只有是精品50| 亚洲伊人久久精品综合 | 日韩精品青青久久久久久| 国产在视频线精品| 白带黄色成豆腐渣| 91在线精品国自产拍蜜月| 国产午夜精品一二区理论片| 床上黄色一级片| 日韩欧美国产在线观看| 26uuu在线亚洲综合色| 中文字幕精品亚洲无线码一区| 插逼视频在线观看| av国产久精品久网站免费入址| 成年av动漫网址| 自拍偷自拍亚洲精品老妇| 在线免费十八禁| 熟女电影av网| 能在线免费看毛片的网站| 日本熟妇午夜| 国产视频首页在线观看| 小蜜桃在线观看免费完整版高清| 国产国拍精品亚洲av在线观看| 尤物成人国产欧美一区二区三区| 最新中文字幕久久久久| 一卡2卡三卡四卡精品乱码亚洲| 久久精品人妻少妇| 国内精品宾馆在线| 熟女电影av网| 精品久久久久久久末码| 久久久久久久久中文| 国产午夜精品久久久久久一区二区三区| 国产av不卡久久| 精品熟女少妇av免费看| 综合色丁香网| 国产高清三级在线| 日产精品乱码卡一卡2卡三| 中文字幕熟女人妻在线| 国产乱人偷精品视频| 22中文网久久字幕| 日日摸夜夜添夜夜爱| 国产免费又黄又爽又色| 男的添女的下面高潮视频| 久久久久久久午夜电影| 久久精品国产99精品国产亚洲性色| 日韩,欧美,国产一区二区三区 | 一个人观看的视频www高清免费观看| 亚洲成人久久爱视频| 99久久精品一区二区三区| 亚洲国产精品sss在线观看| 免费观看的影片在线观看| 中文字幕av在线有码专区| 91午夜精品亚洲一区二区三区| 午夜亚洲福利在线播放| 国产精品嫩草影院av在线观看| 亚洲在线观看片| 青春草视频在线免费观看| 欧美日韩在线观看h| av国产久精品久网站免费入址| 国产又黄又爽又无遮挡在线| 国产精品久久久久久久久免| 波多野结衣巨乳人妻| 亚洲成人精品中文字幕电影| 综合色丁香网| 又黄又爽又刺激的免费视频.| 国内精品一区二区在线观看| 国产精品国产三级专区第一集| 日韩欧美三级三区| 男人狂女人下面高潮的视频| 国产av码专区亚洲av| 精品国产三级普通话版| 男插女下体视频免费在线播放| 搞女人的毛片| 久久精品国产亚洲av涩爱| 亚洲av.av天堂| 国产高清视频在线观看网站| 一区二区三区乱码不卡18| 韩国高清视频一区二区三区| 欧美+日韩+精品| 黄片无遮挡物在线观看| 18禁裸乳无遮挡免费网站照片| a级毛色黄片| 又黄又爽又刺激的免费视频.| 成年版毛片免费区| 精品久久久久久久久av| 99热全是精品| 国产成人精品一,二区| 国产精品国产三级专区第一集| 一区二区三区乱码不卡18| 99久久成人亚洲精品观看| 日产精品乱码卡一卡2卡三| 午夜老司机福利剧场| 中国美白少妇内射xxxbb| АⅤ资源中文在线天堂| 男人狂女人下面高潮的视频| 国产女主播在线喷水免费视频网站 | 久久婷婷人人爽人人干人人爱| 三级毛片av免费| 内射极品少妇av片p| 变态另类丝袜制服| 91久久精品国产一区二区三区| 在现免费观看毛片| 久久久久网色| 国产免费男女视频| av又黄又爽大尺度在线免费看 | 日韩欧美国产在线观看| 亚洲精品日韩在线中文字幕| 最近手机中文字幕大全| 中文资源天堂在线| 亚洲精品色激情综合| 亚洲国产成人一精品久久久| 亚洲精品aⅴ在线观看| 69人妻影院| 99热精品在线国产| 国内揄拍国产精品人妻在线| 熟女人妻精品中文字幕| 国产毛片a区久久久久| 91久久精品国产一区二区三区| 精品人妻一区二区三区麻豆| 国产精品无大码| 日韩 亚洲 欧美在线| 亚洲电影在线观看av| 搡老妇女老女人老熟妇| 欧美丝袜亚洲另类| 简卡轻食公司| 日韩,欧美,国产一区二区三区 | 国产在视频线在精品| 久久精品国产自在天天线| 国产中年淑女户外野战色| 欧美精品一区二区大全| 一区二区三区高清视频在线| 欧美色视频一区免费| 少妇高潮的动态图| 搡老妇女老女人老熟妇| 国产精品99久久久久久久久| 国产精品不卡视频一区二区| 色综合站精品国产| 别揉我奶头 嗯啊视频| 国产精品无大码| 一级黄片播放器| 两个人的视频大全免费| 亚洲最大成人手机在线| 欧美高清成人免费视频www| 国产一区二区三区av在线| 青春草国产在线视频| 亚洲国产高清在线一区二区三| or卡值多少钱| 亚洲在线观看片| 免费不卡的大黄色大毛片视频在线观看 | 国产高清有码在线观看视频| 午夜视频国产福利| 黑人高潮一二区| 免费黄网站久久成人精品| 黄片无遮挡物在线观看| 人妻少妇偷人精品九色| 精品不卡国产一区二区三区| 免费观看性生交大片5| 国产亚洲精品av在线| 长腿黑丝高跟| 十八禁国产超污无遮挡网站| 蜜臀久久99精品久久宅男| 在线免费十八禁| 插逼视频在线观看| 免费搜索国产男女视频| 亚洲性久久影院| 久久久久久久午夜电影| 十八禁国产超污无遮挡网站| 精品国产三级普通话版| 一级av片app| 亚洲精品自拍成人| 欧美色视频一区免费| 夜夜看夜夜爽夜夜摸| 国产精品久久久久久av不卡| 高清视频免费观看一区二区 | 欧美xxxx性猛交bbbb| 国产在视频线在精品| 中文精品一卡2卡3卡4更新| 亚洲在线观看片| 久久久久久久午夜电影| 一级毛片aaaaaa免费看小| 嫩草影院精品99| 看免费成人av毛片| 男插女下体视频免费在线播放| 日韩欧美 国产精品| 一级黄片播放器| 欧美成人午夜免费资源| 久久国产乱子免费精品| 99久久成人亚洲精品观看| 男人舔女人下体高潮全视频| 国产亚洲av片在线观看秒播厂 | 狂野欧美激情性xxxx在线观看| 欧美+日韩+精品| av卡一久久| 亚洲欧美精品自产自拍| 国产人妻一区二区三区在| 亚洲欧美精品专区久久| 99久久成人亚洲精品观看| 国产av不卡久久| 日日摸夜夜添夜夜添av毛片| 嫩草影院入口| 欧美一区二区精品小视频在线| 人妻系列 视频| 最新中文字幕久久久久| 啦啦啦观看免费观看视频高清| 麻豆国产97在线/欧美| 天堂av国产一区二区熟女人妻| 有码 亚洲区| 亚洲av中文字字幕乱码综合| 天堂√8在线中文| 18禁裸乳无遮挡免费网站照片| 波多野结衣巨乳人妻| 亚洲国产欧美人成| 亚洲av免费高清在线观看| 国产av码专区亚洲av| 欧美成人一区二区免费高清观看| 国产精品一区二区在线观看99 | 老司机影院毛片| 精品久久久久久久久av| 日韩欧美精品免费久久| 精品不卡国产一区二区三区| 欧美xxxx黑人xx丫x性爽| 婷婷色av中文字幕| 精品国产三级普通话版| 国产大屁股一区二区在线视频| 亚洲美女视频黄频| 91精品国产九色| 永久免费av网站大全| 午夜激情欧美在线| 村上凉子中文字幕在线| 又爽又黄无遮挡网站| 国产精品国产三级国产专区5o | 看黄色毛片网站| 国产免费一级a男人的天堂| 麻豆国产97在线/欧美| 熟女人妻精品中文字幕| 亚洲无线观看免费| 热99在线观看视频| 亚洲国产最新在线播放| 最后的刺客免费高清国语| 中文字幕久久专区| 日韩精品有码人妻一区| 国产精品嫩草影院av在线观看| 老女人水多毛片| 在线播放无遮挡| 日本五十路高清| 成人漫画全彩无遮挡| 欧美成人免费av一区二区三区| 国产女主播在线喷水免费视频网站 | 久久久久久久久久久丰满| 国产午夜福利久久久久久| 91精品伊人久久大香线蕉| 中文字幕久久专区| 日日撸夜夜添| 免费看日本二区| 日日啪夜夜撸| 男人的好看免费观看在线视频| 一个人看视频在线观看www免费| 欧美精品国产亚洲| 欧美高清性xxxxhd video| 亚洲精品,欧美精品| 好男人在线观看高清免费视频| 夫妻性生交免费视频一级片| 男女国产视频网站| 99热这里只有是精品50| 亚洲美女视频黄频| 青春草国产在线视频| 亚洲中文字幕一区二区三区有码在线看| 女人十人毛片免费观看3o分钟| 日日摸夜夜添夜夜添av毛片| 亚洲熟妇中文字幕五十中出| 麻豆av噜噜一区二区三区| 一区二区三区高清视频在线| 亚洲国产色片| 免费av不卡在线播放| 国产黄色小视频在线观看| 可以在线观看毛片的网站| 久久久久久伊人网av| 久久草成人影院| 亚洲精品日韩在线中文字幕| 好男人在线观看高清免费视频| 精品国产露脸久久av麻豆 | 黄片无遮挡物在线观看| 久久久久网色| 免费av不卡在线播放| 亚洲精品影视一区二区三区av| 国产成人freesex在线| 亚洲最大成人av| 亚洲人成网站在线观看播放| 丰满人妻一区二区三区视频av| 精品久久久久久久久久久久久| 美女国产视频在线观看| 99久久人妻综合| 午夜a级毛片| 日韩 亚洲 欧美在线| 亚洲欧美成人精品一区二区| 久久久精品94久久精品| 中文资源天堂在线| av视频在线观看入口| 久久精品国产亚洲av涩爱| 乱人视频在线观看| 国产综合懂色| 亚洲av男天堂| 午夜精品国产一区二区电影 | 乱人视频在线观看| 亚洲性久久影院| 欧美一区二区国产精品久久精品| 波多野结衣巨乳人妻| 听说在线观看完整版免费高清| a级一级毛片免费在线观看| 国产亚洲91精品色在线| 国产亚洲av片在线观看秒播厂 | 亚洲色图av天堂| 亚洲av一区综合| 欧美人与善性xxx| 噜噜噜噜噜久久久久久91| 欧美变态另类bdsm刘玥| 色综合站精品国产| av专区在线播放| 久久久久久久午夜电影| 看片在线看免费视频| 少妇的逼好多水| av专区在线播放| 丰满乱子伦码专区| 天堂av国产一区二区熟女人妻| 激情 狠狠 欧美| 国产伦理片在线播放av一区| 少妇裸体淫交视频免费看高清| 久久亚洲国产成人精品v| 搡老妇女老女人老熟妇| 看免费成人av毛片| 两个人视频免费观看高清| 少妇熟女aⅴ在线视频| 日本av手机在线免费观看| 国产精品99久久久久久久久| 日本免费a在线| 搞女人的毛片| 91在线精品国自产拍蜜月| 国产毛片a区久久久久| 日日干狠狠操夜夜爽| 国产精品99久久久久久久久| 欧美激情在线99| 免费av不卡在线播放| 亚洲精品国产av成人精品| 三级国产精品欧美在线观看| 午夜爱爱视频在线播放| 日韩av在线免费看完整版不卡| 最近2019中文字幕mv第一页| 午夜福利在线在线| 99久国产av精品| 亚洲av成人精品一二三区| 高清在线视频一区二区三区 | 日本欧美国产在线视频| 成人毛片a级毛片在线播放| 国产精品久久久久久精品电影小说 |