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

    基于加速器7Li(p,n)反應(yīng)的BNCT中子源及慢化體的優(yōu)化設(shè)計

    2022-06-09 01:51:40郭志琦劉昌奇張煒忠戶志鳴徐大鵬姚澤恩
    核技術(shù) 2022年5期
    關(guān)鍵詞:中子通量熱中子產(chǎn)額

    郭志琦 劉昌奇 張煒忠 戶志鳴 徐大鵬 姚澤恩,2

    1(蘭州大學(xué)核科學(xué)與技術(shù)學(xué)院 蘭州 730000)

    2(蘭州大學(xué)教育部中子應(yīng)用技術(shù)工程研究中心 蘭州 730000)

    3(蘭州大學(xué)生命科學(xué)學(xué)院 蘭州 730000)

    中子作為一種高傳能線密度(Linear Energy Transfer,LET)射線,具有較高的輻射生物效應(yīng),所以廣泛應(yīng)用于生命及醫(yī)學(xué)的各個領(lǐng)域[1-2]。其中硼中子俘獲療法(Boron Neutorn Capture Therapy,BNCT)是一種利用中子束和放射療法治療惡性腫瘤的方法,它在常規(guī)治療方法難以治療的惡性腫瘤(如腦膠質(zhì)腫瘤)上[3-5]有著非常廣泛的應(yīng)用前景。含10B藥物進入人體后,會在腫瘤細胞內(nèi)富集,當(dāng)中子束照射腫瘤時,10B(n,α)7Li反應(yīng)生成高LET的α粒子和7Li核,它們能殺死≤10 μm的腫瘤細胞[6],實現(xiàn)在細胞水平上靶向選擇腫瘤細胞,這是其他常規(guī)治療方法無法比擬的,在原理上優(yōu)于目前外科手術(shù)、放射療法、化學(xué)療法、免疫療法和基因療法。BNCT發(fā)展初期,反應(yīng)堆源是主要且唯一的中子源,這限制了BNCT的推廣和發(fā)展[7-9];20世紀(jì)90年代,基于加速器的硼中子俘獲治療(AB-BNCT)提出以后,許多國家都瞄準(zhǔn)了AB-BNCT的研究。2009年,日本成功研制出世界上第一臺基于加速器的硼中子俘獲治療的裝置。近年來美國、芬蘭、英國等十幾個國家也在開展基于加速器的硼中子俘獲治療的研究[10]。中國的散裂中子源于2017年8月開始研發(fā)適用于醫(yī)院的BNCT裝置,在原有的強流質(zhì)子加速器上研制,有望制成國內(nèi)首臺用于BNCT的加速器裝置[11]。相較于核反應(yīng)堆中子源,加速器中子源在keV能區(qū)擁有更高的中子產(chǎn)額,且伴隨放出的γ射線較少,另外還有設(shè)備簡單、成本低及安全性高等優(yōu)勢,故其更適用于醫(yī)用BNCT。加速器中子源一般通過A(a,n)B反應(yīng)產(chǎn)生中子,入射粒子a一般為α等帶電粒子,而靶核一般選用2D、7Li、9Be等中子結(jié)合能低、能級間距大的輕核[12-14]。當(dāng)入射質(zhì)子能量在2.0~3.5 MeV時,基于7Li(p,n)7Be反應(yīng)的中子源的產(chǎn)額更高,更適用于BNCT治療。

    本文主要采用蒙特卡羅程序Geant4模擬計算不同能量的質(zhì)子轟擊不同厚度的鋰靶產(chǎn)生的中子束的能量、角分布及產(chǎn)額等數(shù)據(jù),從而得出合理的基于加速器的硼中子俘獲治療中子靶;在此基礎(chǔ)上設(shè)計中子束流整形體(Beam Shaping Assembly,BSA)的模型,對BSA進行優(yōu)化模擬,計算出最優(yōu)的材料組合及厚度,為AB-BNCT的設(shè)計提供數(shù)據(jù)參考。

    1 7Li(p,n)7Be反應(yīng)中子分布特性研究

    1.1 質(zhì)子打靶模型

    基于低能質(zhì)子加速器的中子源一般采用Li靶或Be靶,利用核反應(yīng)7Li(p,n)和9Be(p,n)反應(yīng)產(chǎn)生中子。當(dāng)質(zhì)子能量低于3.5 MeV,鋰靶的中子產(chǎn)額要 高 于 鈹 靶[1],考 慮 到7Li(p,n)反 應(yīng) 閾 能 為1.88 MeV,因此本文的設(shè)計基于7Li(p,n)反應(yīng),主要研究2.0~3.5 MeV的入射質(zhì)子轟擊鋰靶所產(chǎn)生的中子能譜。質(zhì)子打靶的模型如圖1所示,鋰靶的半徑為1 cm。在Geant4(版本10.06)模擬計算中,物理過程使用的是QGSP_BIC_AIIHP包,QGSP_BIC_AIIHP是一個物理過程集成包,它描述了質(zhì)子和中子與物質(zhì)相互作用的過程。其中,包括電磁物理過程、衰變過程、離子輸運以及強子物理過程等。另外,數(shù)據(jù)庫則是采用Talys程序[15]得到的雙微分截面數(shù)據(jù)庫,即TENDLE1.3.2庫。

    圖1 質(zhì)子打靶模型示意圖Fig.1 Schematic diagram of proton target model

    1.2 7Li(p,n)7Be中子源的設(shè)計結(jié)果與分析

    經(jīng)模擬計算,2.5 MeV的質(zhì)子束轟擊0.2 cm的鋰靶,中子產(chǎn)額為8.446×1011n·mC-1,中子最高能量為0.785 MeV,平均能量為0.308 MeV,最大中子出射角為180°,平均中子出射角為63.23°。以上參數(shù)已有的研究數(shù)據(jù)[16]為:中子產(chǎn)額為8.83×1011n·mC-1,中子最高能量為0.787 MeV,平均能量為0.326 MeV,最大中子出射角為180°,平均中子出射角為62.9°。經(jīng)比對,本文使用Geant4計算得到的數(shù)據(jù)與已有的研究數(shù)據(jù)吻合度較高,故而證明使用Geant4模擬計算是可行的。

    為了得出7Li(p,n)反應(yīng)產(chǎn)生中子的能譜特性,本文采用蒙特卡羅程序Geant4對不同能量的質(zhì)子轟擊不同厚度的鋰靶的中子產(chǎn)額、出射中子的平均能量、角分布等進行了模擬計算。由于質(zhì)子轟擊鋰的反應(yīng)閾能為1.88 MeV,所以選用能量為2 MeV、2.5 MeV、3 MeV和3.5 MeV的質(zhì)子分別轟擊不同厚度的鋰靶,設(shè)計鋰靶的厚度分別為0.04 mm、0.06 mm、0.08 mm、0.1 mm、0.15 mm、0.2 mm、0.3 mm和0.5 mm。本文將從平均能量及中子產(chǎn)額進行分析,得出合適的基于加速器7Li(p,n)反應(yīng)的BNCT中子源,結(jié)果如圖2所示。入射質(zhì)子的能量為2 MeV,鋰靶厚度為0.04~0.5 mm時,質(zhì)子轟擊鋰靶產(chǎn)生的中子產(chǎn)額以及平均能量基本不變。質(zhì)子能量大于等于2.5 MeV,鋰靶厚度大于0.1 mm時,中子產(chǎn)額相較于2 MeV時提升了一個數(shù)量級。如圖2所示,質(zhì)子能量為2.5 MeV、3 MeV和3.5 MeV,對應(yīng)鋰靶厚度分別為0.1 mm、0.2 mm和0.3 mm時,它們的中子產(chǎn)額和平均中子能量都達到了穩(wěn)定值,之后再增加鋰靶厚度,中子產(chǎn)額和中子平均能量不再變化。入射質(zhì)子能量為2 MeV時,平均能量為0.074 MeV,接近適合BNCT的超熱中子能量,但其中子產(chǎn)額偏低。當(dāng)鋰靶厚度是0.2 mm時,入射質(zhì)子能量3 MeV的中子產(chǎn)額是入射質(zhì)子能量2.5 MeV的兩倍,而其中子平均能量為0.47 MeV,只比2.5 MeV的入射質(zhì)子高0.16 MeV。入射質(zhì)子能量為3.5 MeV的中子產(chǎn)額雖然最高,但是中子平均能量達0.706 MeV,這將給后面的慢化工作造成一定的困難。因此,7Li(p,n)7Be反應(yīng)中的鋰靶厚度建議設(shè)計為0.2 mm,入射質(zhì)子能量選用3 MeV,對應(yīng)此時的中子產(chǎn)額為2.859×10-4n·s-1,中子平均能量為0.47 MeV。

    圖2 中子平均能量隨鋰靶厚度的變化趨勢(a)和中子產(chǎn)額隨鋰靶厚度的變化趨勢(b)Fig.2 Variation trends of average neutron energy(a)and neutron yield(b)with lithium target thickness

    1.3 7Li(p,n)7Be能譜分布及角分布

    上述§1.2中通過分析得出使用3 MeV的質(zhì)子轟擊0.2 mm的鋰靶產(chǎn)生的中子束適合作為基于加速器7Li(p,n)反應(yīng)的BNCT中子源。圖3(a)為不同能量的質(zhì)子轟擊0.2 mm鋰靶的能譜圖。當(dāng)入射質(zhì)子能量為2.5 MeV和3 MeV時,中子能譜出現(xiàn)兩個峰,其峰位分別是0.09 MeV和0.46 MeV,第二個峰出現(xiàn)是因為7Li(p,n)反應(yīng)在2.25 MeV處有個共振峰,圖3(b)[14]為7Li(p,n)反應(yīng)截面圖;當(dāng)入射中子能量為3.5 MeV時,由于入射質(zhì)子能量較高,鋰靶較薄,質(zhì)子在鋰靶中的能量沉積較低,從而導(dǎo)致出射中子能譜呈現(xiàn)單峰。

    圖3 不同能量質(zhì)子入射0.2 mm的鋰靶的出射中子能譜圖(a)和7Li(p,n)反應(yīng)截面(b)Fig.3 Outgoing neutron spectra of different energies of protons incident into a 0.2 mm lithium target(a)and7Li(p,n)reaction cross section(b)

    質(zhì)子轟擊鋰靶產(chǎn)生的中子能量較高,不符合BNCT設(shè)計要求,因此有必要對其開展慢化、準(zhǔn)直工作以達到實際指標(biāo)要求。其中,中子慢化層設(shè)計需要對中子角分布進行分析,以實現(xiàn)合適能區(qū)中子產(chǎn)額最大化利用。圖4給出了出射中子立體角(dY/dΩ)分布以及空間弧度角(dY/dθ)分布計算結(jié)果。由圖4(a)可以看出,中子角分布在0°~90°范圍內(nèi)具有前傾性,當(dāng)單位角度為50°時中子產(chǎn)額達到最大,另外,出射中子在0°~90°范圍內(nèi)中子積分產(chǎn)額達到了總產(chǎn)額77%以上。因此,基于中子前傾性質(zhì),采用質(zhì)子入射方向與BSA軸線方向夾角呈0°的設(shè)計對靶出射中子進行慢化、準(zhǔn)直。圖5為3 MeV入射質(zhì)子轟擊0.2 mm鋰靶產(chǎn)生的出射中子能譜-角度二維分布計算結(jié)果??梢钥闯?,中子能量集中在0.09~0.46 MeV,與圖3(a)能譜所給出結(jié)論一致,另外,出射中子角分布具有明顯前傾性。該7Li(p,n)反應(yīng)出射中子能譜-角度分布將作為后續(xù)BSA慢化程序輸入中子源,通過外耦合方式開展對基于BNCT中子源BSA設(shè)計的研究。

    圖4 不同能量質(zhì)子轟擊0.2 mm鋰靶出射中子的立體角分布(a)和空間弧度角分布(b)Fig.4 Solid angle distribution(a)and radian angle distribution(b)of the outgoing neutron of 0.2 mm lithium target bombarded by protons with different energies

    圖5 3 MeV質(zhì)子轟擊0.2 mm鋰靶能譜-角度二維圖Fig.5 Energy spectrum-angular two-dimensional diagram of a 0.2-mm lithium target bombarded by 3-MeV protons

    2 BSA的設(shè)計與評價

    2.1 BSA的模型設(shè)計

    根據(jù)IAEA-TECDOC-1223報告推薦的指數(shù)[15],用于BNCT的中子束必須滿足以下條件:1)超熱中子通量Φepi>1.0×109n·cm-2·s-1;2)Φepi與熱中子通量Φth的比值Φepi/Φth>100;3)Φepi與快中子通量Φfast的比值Φepi/Φfast>20;4)γ射線劑量Dγ與超熱中子通量Φepi的 比 值Dγ/Φepi<2×10-13Gy·cm2·n-1(epi);詳見表1。

    表1 BSA的設(shè)計參數(shù)要求Table 1 BSA design parameter requirements

    為了滿足國際原子能機構(gòu)IAEA-TECDOC-1223報告推薦的指數(shù),本文采用Geant4對BSA進行了設(shè)計。BSA主要由中子通道、反射體、慢化層、熱中子吸收層、γ吸收層以及準(zhǔn)直器組成,通過模擬優(yōu)化,給出合適的慢化層、熱中子吸收層以及γ吸收層的材料和厚度。

    2.2 單種材料作為慢化層的設(shè)計

    BSA主要由中子通道、慢化層、熱中子吸收層、γ吸收層以及準(zhǔn)直器組成。常見的中子慢化材料有氟化鎂(MgF2)、氟化鋰(LiF)、氟化鈦(TiF3)、氟化鋁(AlF3)、Fluental(復(fù)合材料,69% AlF3+30% Al+1%LiF)、聚乙烯(PE)、重水(D2O)、石墨(C)等。為了評價上述材料的慢化性能,利用Geant4對單獨的每種材料的慢化性能做了模擬計算。慢化效果主要與慢化層的厚度相關(guān),選擇材料厚度區(qū)間為10~35 cm,以5 cm為步長,模擬上述8種材料的Φepi、Φepi/Φth、Φepi/Φfast、Dγ/Φepi數(shù)值,如圖6(a)所示。由圖6(a)可知,隨著慢化層厚度的增加,超熱中子通量呈先增大后減小的趨勢(PE和D2O除外),在15 cm左右達到最大值。3 MeV的質(zhì)子束轟擊鋰靶產(chǎn)生的中子的平均能量只有0.47 MeV,使用PE和D2O做慢化材料會出現(xiàn)過慢化的情況,導(dǎo)致有一部分能量較低的中子會直接被吸收,超熱中子通量會降得很低,所以PE和D2O不適合做3 MeV的質(zhì)子轟擊鋰靶產(chǎn)生的中子束的慢化材料。

    圖6 超熱中子通量與材料厚度(a)及超熱中子占比(b)的關(guān)系Fig.6 Relationship between epithermal neutron flux and material thickness(a)and proportion of epithermal neutrons(b)

    理想的用于BNCT治療的中子束,除了Φepi達到所要求的參數(shù)值外,最重要的是提高超熱中子所占比例,盡可能減少熱中子和快中子,提高中子束的“質(zhì)量”。所以將超熱中子通量在總中子通量(Φtotal)所占比例Φepi/Φtotal作為橫坐標(biāo),以超熱中子通量Φepi為縱坐標(biāo)作圖,選出當(dāng)Φepi/Φtotal最高時,Φepi也最大時的材料,該材料就是單種材料做慢化體時的最佳材料。如圖6(b)所示,單獨使用一種材料做慢化層時,F(xiàn)luental和AlF3兩種材料在Φepi/Φtotal相等時,超熱中子通量Φepi遠低于其他幾種材料,所以這兩種材料可以排除,不予考慮。MgF2、LiF、TiF3、C 4種材料中,當(dāng)Φepi/Φtotal相等時,MgF2的超熱中子通量最高。所以,單獨使用一種材料作為慢化層時,MgF2是最好的選擇。

    2.3 不同種材料混合作為慢化層的設(shè)計

    為了找出最合適的慢化層組合方式,提高超熱中子的通量,本文通過模擬嘗試了上述幾種材料的組合。如圖6(a)所示,當(dāng)慢化體厚度大于15 cm,在相同的厚度下,TiF3的超熱中子通量最高,從上面的分析中得出單獨使用一種材料時,MgF2是最合適的材料。MgF2在10 keV很多個共振峰[17-18],可以使快中子快速失去能量,所以將MgF2放在臨近中子源的第一層,可以快速降低快中子能量并減少慢化層的總體厚度;TiF3中的19F的第一、第二激發(fā)態(tài)分別為109.9 keV和197.1 keV,在100 keV以上有較高的非彈性散射截面,可以用來降低經(jīng)過第一層MgF2慢化后能量還未減小的快中子,TiF3中的48Ti是一個原子序數(shù)較大的元素,對慢化后的超熱中子吸收較小,所以將TiF3設(shè)置為第二層慢化材料。通過前面的模擬發(fā)現(xiàn),將3 MeV質(zhì)子束轟擊0.2 mm的鋰靶產(chǎn)生的中子慢化至適合用于BNCT治療的中子束,所需要的慢化層厚度大約為30~36 cm。通過大量的模擬優(yōu)化后發(fā)現(xiàn),MgF2的厚度在24~29 cm、TiF3的厚度在6~12 cm時,經(jīng)過慢化后的中子束Φepi/Φth、Φepi/Φfast、Dγ/Φepi等參數(shù)在滿足IAEA的標(biāo)準(zhǔn)的條件下,超熱中子通量可以達到一個較高的水平。如圖7所示,隨著MgF2和TiF3厚度的增加,Φepi呈一個下降的趨勢。從圖7可以看出,并不是所有的組合方式超熱中子通量Φepi都是大于1×109n·cm-2·s-1(入射質(zhì)子流強為15 mA),為了找出最合適的組合方式,將圖7中超熱中子通量Φepi大于1×109n·cm-2·s-1且其他參數(shù)符合表1要求的組合方式列于表2。

    圖7 MgF2+TiF3組合作慢化層時的Φepi結(jié)果Fig.7 Φepi results of MgF2+TiF3 cooperative moderating layer

    如表2所示,滿足BSA出口參數(shù)要求的組合方式一共有4個,雖然“25 cm MgF2+11 cm TiF3”組合擁有更高的Φepi/Φth值,但是它的Dγ/Φepi值相對偏高;“26 cm MgF2+10 cm TiF3”的組合Φepi/Φth值略低,但是其他參數(shù)都更符合BSA出口參數(shù)的要求,所以“26 cm MgF2+10 cm TiF3”的組合是更好的選擇。

    表2 4種“MgF2+TiF3”組合的能譜參數(shù)Table 2 Energy spectrum parameters of four MgF2+TiF3 combinations

    2.4 熱中子吸收層的設(shè)計

    經(jīng)過慢化后的中子束不可避免地會產(chǎn)生一些熱中子,而熱中子主要會沉積在人的皮膚層,對人皮膚造成損傷,所以需要在慢化體之后增加一層熱中子吸收材料來降低熱中子的通量。6Li是一種良好的熱中子吸收材料,它的熱中子吸收截面很高(940 b),且6Li(n,α)核反應(yīng)產(chǎn)生的γ光子更少,所以本文選擇的熱中子吸收材料是6Li,質(zhì)量密度為0.11 g·cm-3。我們將慢化體的厚度定為“26 cm MgF2+10 cm TiF3”的組合,對熱中子吸收層的厚度進行優(yōu)化模擬。如圖8所示,隨著6Li的厚度的增加,Φepi/Φth呈上升趨勢,當(dāng)6Li的厚度為0.08 cm時,Φepi/Φth=125.23,符合BSA參數(shù)設(shè)計要求,因此,熱中子吸收層厚度為0.08 cm。

    圖8 經(jīng)過慢化層后的中子束在不同6Li厚度下的Φepi/Φth值Fig.8 Φepi/Φth values of the neutron beam after moderating layer at different6Li thicknesses

    2.5 γ吸收層的設(shè)計

    在BSA設(shè)計過程中,中子穿透慢化材料時,必然會產(chǎn)生γ射線,而多余的γ劑量會對人體產(chǎn)生危害,因此需要增加一層吸收材料吸收γ射線。最常見的的γ屏蔽材料是鉛,但是鉛對生物體有一定的毒性,所以我們選擇了具有相同γ射線屏蔽能力,但更綠色無污染的材料鉍(Bi)。通過改變鉍厚度,對整個BNCT系統(tǒng)進行了模擬計算,以選取最優(yōu)材料厚度參數(shù)。結(jié)果如圖9所示,隨著鉍厚度的增加,Dγ/Φepi呈下降趨勢,當(dāng)鉍的厚度為1.5 cm時,Dγ/Φepi=1.82×10-13Gy·cm2·n-1(epi),此時處于較低劑量水平,滿足BSA出口束流參數(shù)的要求。因此,γ吸收層鉍的厚度設(shè)計為1.5 cm。

    圖9 Dγ/Φepi與Bi厚度的關(guān)系Fig.9 Relationship between Dγ/Φepi and Bi thickness

    2.6 BSA的設(shè)計結(jié)果

    經(jīng)過前面的分析討論,本文最終給出了一種“26 cm MgF2+10 cm TiF3”的慢化層組合,慢化層半徑20 cm;反射體材料選用Teflon,包裹在鋰靶和慢化層的外圍,半徑為40 cm,長度為62 cm;熱中子吸收層位于慢化層下游一側(cè),材料選用6Li,厚度0.08 cm;γ屏蔽層的材料為鉍,厚度1.5 cm;準(zhǔn)直器的材料選用含硼聚乙烯,厚度8 cm,在底面設(shè)置一中子出口,半徑為7.5 cm。慢化體的整體示意圖如圖10所示。慢化后的能譜如圖11所示,由圖11可以看出,經(jīng)過慢化后的中子在超熱能區(qū)(0.5 eV~10 keV)具有很高的通量,這進一步驗證了我們設(shè)計的合理性。

    圖10 BSA整體設(shè)計示意圖1:半徑40 cm的反射體;2:質(zhì)子通道;3:0.02 cm7Li;4:26 cm MgF2;5:10 cm TiF3;6:0.08 cm6Li;7:1.5 cm Bi;8:8 cm含硼聚乙烯;9:中子出口Fig.10 Schematic diagram of BSA overall design1:Reflector with radius of 40 cm;2:Proton channel;3:0.02 cm7Li;4:26 cm MgF2;5:10 cm TiF3;6:0.08 cm6Li;7:1.5 cm Bi;8:8 cm boron-polyethylene;9:Neutron outlet

    圖11 慢化后的中子能譜圖Fig.11 Neutron energy spectrum after smoderation

    3 結(jié)語

    本工作利用Geant4研究了基于加速器7Li(p,n)反應(yīng)的BNCT中子源的能譜角分布,設(shè)計了3 MeV的質(zhì)子轟擊0.2 mm鋰靶的方案,其產(chǎn)生的中子束的平均能量為0.47 MeV,中子產(chǎn)額為2.859×10-4n·s-1。另外,通過外耦合方式將BNCT中子源能譜角分布與BSA設(shè)計工作結(jié)合,開展了基于加速器7Li(p,n)反應(yīng)的BNCT中子源BSA設(shè)計研究,提出了“26 cm MgF2+10 cm TiF3”的組合作為慢化層,0.08 cm的6Li作為熱中子吸收層,1.5 cm的Bi作為γ射線的吸收層。模擬計算結(jié)果表明:入射質(zhì)子流為15 mA時,Φepi=1.16×109n·cm-2·s-1,快 中子 與 超 熱 中子 比 率Φepi/Φfast=21.64,熱中子與超熱中子比率Φepi/Φth=125.23,γ成 分Dγ/Φepi=1.82×10-13Gy·cm2·n-1(epi)。本文的研究結(jié)果可以為AB-BNCT的設(shè)計提供重要的數(shù)據(jù)參考,為推動我國BNCT的發(fā)展提供重要的理論依據(jù)。

    作者貢獻聲明郭志琦:負責(zé)文章的數(shù)據(jù)模擬,起草撰寫以及最終版本的修訂;劉昌奇:負責(zé)文章模擬程序的編寫及修改;張煒忠:負責(zé)文章資料的查閱及整理;戶志鳴:負責(zé)文章資料的查閱和整理;徐大鵬:負責(zé)文章的修改;姚澤恩:負責(zé)文章的修改和整體把握。

    猜你喜歡
    中子通量熱中子產(chǎn)額
    熱中子透射成像轉(zhuǎn)換屏物理設(shè)計研究
    一個可靠和準(zhǔn)確的光電產(chǎn)額譜模型及應(yīng)用*
    基于貝葉斯更新方法的235U熱中子獨立裂變產(chǎn)額協(xié)方差估計
    針對裂變產(chǎn)額和半衰期的燃耗計算靈敏度和不確定度分析方法
    基于協(xié)同進化的航空高度單粒子翻轉(zhuǎn)故障生成方法研究
    單晶硅受照熱中子注量率的雙箔活化法測量研究
    裂變產(chǎn)物活度計算通用程序開發(fā)
    修正快中子通量以提高碳氧測量精度的研究
    某型“三代”核電機組與M310機組堆芯測量系統(tǒng)
    脈沖中子-裂變中子探測鈾黃餅的MCNP模擬
    一级a爱视频在线免费观看| 交换朋友夫妻互换小说| 亚洲欧美成人综合另类久久久| 50天的宝宝边吃奶边哭怎么回事| 国产亚洲av高清不卡| 国产在线观看jvid| 久久精品国产亚洲av高清一级| 国产麻豆69| 一区二区三区精品91| e午夜精品久久久久久久| 啦啦啦中文免费视频观看日本| 亚洲精品第二区| h视频一区二区三区| 黄频高清免费视频| 爱豆传媒免费全集在线观看| 男人爽女人下面视频在线观看| 免费黄频网站在线观看国产| 亚洲av日韩在线播放| 一本—道久久a久久精品蜜桃钙片| 亚洲中文日韩欧美视频| 亚洲国产成人一精品久久久| 免费高清在线观看日韩| 少妇人妻久久综合中文| 国产成+人综合+亚洲专区| 国产精品影院久久| 黄频高清免费视频| 多毛熟女@视频| 亚洲一区中文字幕在线| av视频免费观看在线观看| av网站免费在线观看视频| 中文精品一卡2卡3卡4更新| 美国免费a级毛片| 丝瓜视频免费看黄片| 岛国毛片在线播放| 99久久精品国产亚洲精品| 精品少妇内射三级| 极品少妇高潮喷水抽搐| 男女之事视频高清在线观看| 精品一区二区三区四区五区乱码| 日韩免费高清中文字幕av| 97在线人人人人妻| 男女免费视频国产| 丰满人妻熟妇乱又伦精品不卡| 国产精品久久久久久精品电影小说| 波多野结衣av一区二区av| 俄罗斯特黄特色一大片| 精品欧美一区二区三区在线| 日韩 亚洲 欧美在线| 午夜福利乱码中文字幕| 亚洲精品自拍成人| 久久久久精品国产欧美久久久 | 波多野结衣av一区二区av| 美女脱内裤让男人舔精品视频| 极品人妻少妇av视频| 丁香六月欧美| 日韩欧美国产一区二区入口| 大陆偷拍与自拍| 精品久久久久久电影网| 亚洲欧美清纯卡通| 国产精品香港三级国产av潘金莲| av网站免费在线观看视频| 亚洲国产中文字幕在线视频| 亚洲中文字幕日韩| 国产精品九九99| 99热国产这里只有精品6| 日日摸夜夜添夜夜添小说| 老熟妇仑乱视频hdxx| 国产一区二区三区综合在线观看| 99久久精品国产亚洲精品| 国产成人一区二区三区免费视频网站| 黄色毛片三级朝国网站| 国产野战对白在线观看| 精品福利永久在线观看| 高清在线国产一区| 超碰97精品在线观看| 99re6热这里在线精品视频| 亚洲色图 男人天堂 中文字幕| 黄片播放在线免费| 国产深夜福利视频在线观看| 亚洲一码二码三码区别大吗| 成人三级做爰电影| www.精华液| av不卡在线播放| 巨乳人妻的诱惑在线观看| kizo精华| 两性夫妻黄色片| 日日夜夜操网爽| 伊人久久大香线蕉亚洲五| 久久女婷五月综合色啪小说| 十分钟在线观看高清视频www| 日韩免费高清中文字幕av| 日本欧美视频一区| 人妻人人澡人人爽人人| 日本欧美视频一区| 亚洲精品久久久久久婷婷小说| svipshipincom国产片| 国产视频一区二区在线看| 亚洲成人免费av在线播放| 成年动漫av网址| 真人做人爱边吃奶动态| 午夜福利视频精品| 91九色精品人成在线观看| 99国产精品一区二区三区| 欧美 亚洲 国产 日韩一| 国产一级毛片在线| 在线精品无人区一区二区三| 成年人黄色毛片网站| 亚洲精品国产色婷婷电影| 欧美人与性动交α欧美软件| 亚洲国产精品成人久久小说| 极品人妻少妇av视频| 亚洲免费av在线视频| 亚洲精品国产色婷婷电影| 国产1区2区3区精品| 999精品在线视频| 久热这里只有精品99| 国产精品一区二区在线观看99| 国产精品成人在线| 2018国产大陆天天弄谢| 国产高清videossex| 午夜免费鲁丝| 一本综合久久免费| 久久久久久免费高清国产稀缺| 国产淫语在线视频| 大陆偷拍与自拍| 免费在线观看黄色视频的| 久久ye,这里只有精品| 亚洲自偷自拍图片 自拍| 纯流量卡能插随身wifi吗| 国产又色又爽无遮挡免| 成年美女黄网站色视频大全免费| 亚洲人成电影免费在线| 欧美少妇被猛烈插入视频| 亚洲欧美清纯卡通| 999久久久国产精品视频| 国产精品九九99| 欧美黄色淫秽网站| 国产亚洲一区二区精品| 国产精品一区二区在线观看99| 大码成人一级视频| 欧美亚洲日本最大视频资源| 欧美xxⅹ黑人| 高潮久久久久久久久久久不卡| 伊人久久大香线蕉亚洲五| 亚洲欧美精品自产自拍| 久久99热这里只频精品6学生| 极品少妇高潮喷水抽搐| 中亚洲国语对白在线视频| 91老司机精品| 免费少妇av软件| 久久久久国内视频| 欧美乱码精品一区二区三区| a级片在线免费高清观看视频| 午夜精品久久久久久毛片777| av电影中文网址| 精品人妻在线不人妻| 亚洲美女黄色视频免费看| 国产成人精品久久二区二区免费| 久热这里只有精品99| 亚洲色图 男人天堂 中文字幕| 王馨瑶露胸无遮挡在线观看| 久久毛片免费看一区二区三区| 99精品久久久久人妻精品| 美国免费a级毛片| 久久毛片免费看一区二区三区| 中文字幕人妻熟女乱码| 国产精品1区2区在线观看. | 又黄又粗又硬又大视频| 国产不卡av网站在线观看| 99热国产这里只有精品6| 热99国产精品久久久久久7| 国产91精品成人一区二区三区 | 黄色毛片三级朝国网站| 另类亚洲欧美激情| 亚洲av男天堂| 99热全是精品| 欧美xxⅹ黑人| 成人18禁高潮啪啪吃奶动态图| 飞空精品影院首页| a级毛片在线看网站| av在线播放精品| 国产又爽黄色视频| 99热国产这里只有精品6| 女性被躁到高潮视频| 99精品欧美一区二区三区四区| 日本vs欧美在线观看视频| 老司机在亚洲福利影院| 亚洲精品国产精品久久久不卡| 亚洲第一欧美日韩一区二区三区 | 亚洲精品一区蜜桃| 中文欧美无线码| 国产成+人综合+亚洲专区| 宅男免费午夜| 久久久久久亚洲精品国产蜜桃av| 婷婷色av中文字幕| 老鸭窝网址在线观看| 一级毛片女人18水好多| 两个人免费观看高清视频| bbb黄色大片| 欧美日韩av久久| 黄色视频不卡| 91大片在线观看| 国产精品影院久久| 欧美人与性动交α欧美精品济南到| 夜夜骑夜夜射夜夜干| 黄片小视频在线播放| 夜夜夜夜夜久久久久| 秋霞在线观看毛片| 满18在线观看网站| 9191精品国产免费久久| 亚洲成人国产一区在线观看| 亚洲精品国产一区二区精华液| 欧美 日韩 精品 国产| 男女免费视频国产| 咕卡用的链子| 黄色视频不卡| 91大片在线观看| 一个人免费看片子| 99热国产这里只有精品6| 97人妻天天添夜夜摸| 美女高潮到喷水免费观看| 精品久久蜜臀av无| 一区在线观看完整版| 午夜91福利影院| 精品乱码久久久久久99久播| 国产成人啪精品午夜网站| 999久久久国产精品视频| 国产精品 国内视频| 少妇猛男粗大的猛烈进出视频| 午夜免费观看性视频| 国产成人av教育| 热99国产精品久久久久久7| 老汉色av国产亚洲站长工具| 欧美精品av麻豆av| 日韩电影二区| 少妇的丰满在线观看| 亚洲欧美激情在线| 久久人妻福利社区极品人妻图片| 国产免费视频播放在线视频| 日韩欧美一区视频在线观看| 国产成人欧美在线观看 | 亚洲精品第二区| 久久中文字幕一级| av在线播放精品| 国产av一区二区精品久久| 菩萨蛮人人尽说江南好唐韦庄| 国产亚洲av高清不卡| 免费女性裸体啪啪无遮挡网站| 看免费av毛片| 久久亚洲国产成人精品v| 久久久久视频综合| 97在线人人人人妻| 三级毛片av免费| 麻豆乱淫一区二区| 亚洲av欧美aⅴ国产| 午夜91福利影院| 中文字幕人妻丝袜一区二区| 91大片在线观看| 国产真人三级小视频在线观看| 欧美人与性动交α欧美精品济南到| 欧美日韩精品网址| 男女之事视频高清在线观看| 搡老熟女国产l中国老女人| 中文字幕人妻丝袜一区二区| 国产成人啪精品午夜网站| 女人高潮潮喷娇喘18禁视频| 日本一区二区免费在线视频| 国产一级毛片在线| 中文欧美无线码| 9热在线视频观看99| 精品久久久久久久毛片微露脸 | 两个人免费观看高清视频| 桃花免费在线播放| 97在线人人人人妻| 国产一区二区在线观看av| 欧美日韩视频精品一区| 91九色精品人成在线观看| www.999成人在线观看| 日本猛色少妇xxxxx猛交久久| 精品人妻一区二区三区麻豆| 69精品国产乱码久久久| 动漫黄色视频在线观看| 黄片大片在线免费观看| 欧美成人午夜精品| a级毛片在线看网站| 亚洲自偷自拍图片 自拍| 久久精品人人爽人人爽视色| 免费一级毛片在线播放高清视频 | 一区二区三区精品91| 亚洲午夜精品一区,二区,三区| 俄罗斯特黄特色一大片| 丰满少妇做爰视频| 国产男女超爽视频在线观看| av国产精品久久久久影院| kizo精华| 久久午夜综合久久蜜桃| 国产不卡av网站在线观看| 国产亚洲欧美精品永久| 国产在线观看jvid| av又黄又爽大尺度在线免费看| 欧美+亚洲+日韩+国产| 国产av又大| 国产精品免费视频内射| 最近最新免费中文字幕在线| 精品福利观看| 国产精品久久久久久人妻精品电影 | 亚洲精品国产色婷婷电影| 美女午夜性视频免费| 国产国语露脸激情在线看| 啦啦啦啦在线视频资源| 国产精品偷伦视频观看了| 丝袜脚勾引网站| 中文字幕人妻丝袜制服| 日韩免费高清中文字幕av| 久久女婷五月综合色啪小说| 亚洲国产欧美网| 国产一区二区三区综合在线观看| 脱女人内裤的视频| 亚洲国产日韩一区二区| 色视频在线一区二区三区| 亚洲成人国产一区在线观看| 国产日韩一区二区三区精品不卡| 国产91精品成人一区二区三区 | 亚洲精品成人av观看孕妇| 精品福利永久在线观看| 777久久人妻少妇嫩草av网站| 免费观看av网站的网址| 又黄又粗又硬又大视频| 狠狠婷婷综合久久久久久88av| 大片免费播放器 马上看| 国产欧美亚洲国产| 在线观看免费高清a一片| 捣出白浆h1v1| 窝窝影院91人妻| 十八禁网站免费在线| 久久久久久久久免费视频了| 亚洲精品美女久久久久99蜜臀| 欧美精品av麻豆av| 久久中文看片网| 久久ye,这里只有精品| av视频免费观看在线观看| 免费久久久久久久精品成人欧美视频| 久久亚洲精品不卡| 丰满人妻熟妇乱又伦精品不卡| 精品一品国产午夜福利视频| 国产精品 欧美亚洲| 精品一区二区三卡| 国产亚洲午夜精品一区二区久久| 久热这里只有精品99| 熟女少妇亚洲综合色aaa.| 国产真人三级小视频在线观看| 狂野欧美激情性xxxx| 丝袜喷水一区| 91麻豆av在线| 国产主播在线观看一区二区| 正在播放国产对白刺激| 美女高潮喷水抽搐中文字幕| 脱女人内裤的视频| 成人亚洲精品一区在线观看| 99九九在线精品视频| 日韩视频一区二区在线观看| 一区二区三区乱码不卡18| 精品国产一区二区三区久久久樱花| 亚洲国产欧美网| 国产男女超爽视频在线观看| 亚洲专区字幕在线| 久久99一区二区三区| 热99国产精品久久久久久7| 性色av乱码一区二区三区2| 美国免费a级毛片| 久久人妻福利社区极品人妻图片| www日本在线高清视频| 性色av乱码一区二区三区2| 我要看黄色一级片免费的| 91av网站免费观看| 老司机影院成人| 久久99一区二区三区| 成年人午夜在线观看视频| 永久免费av网站大全| 9热在线视频观看99| 母亲3免费完整高清在线观看| 下体分泌物呈黄色| a在线观看视频网站| 大片电影免费在线观看免费| 男女无遮挡免费网站观看| av网站在线播放免费| 一区福利在线观看| 国产日韩一区二区三区精品不卡| 亚洲精品美女久久av网站| xxxhd国产人妻xxx| 一个人免费在线观看的高清视频 | 国产一区二区三区综合在线观看| 黄色a级毛片大全视频| 久久久久精品人妻al黑| 国产成人a∨麻豆精品| 国产男女超爽视频在线观看| 欧美性长视频在线观看| 一级片免费观看大全| av免费在线观看网站| 汤姆久久久久久久影院中文字幕| 悠悠久久av| 久久九九热精品免费| 亚洲精品在线美女| 国产成人av教育| 97在线人人人人妻| 一级片'在线观看视频| 久久精品aⅴ一区二区三区四区| 国产片内射在线| 黄片小视频在线播放| 黄色怎么调成土黄色| 久久久久精品人妻al黑| 精品国产国语对白av| 成在线人永久免费视频| 色老头精品视频在线观看| 久久天躁狠狠躁夜夜2o2o| 日韩一区二区三区影片| av福利片在线| 精品亚洲成国产av| 高清av免费在线| 精品视频人人做人人爽| 精品第一国产精品| 国产一区二区激情短视频 | 色精品久久人妻99蜜桃| 999精品在线视频| 一进一出抽搐动态| 国产日韩欧美在线精品| 亚洲欧洲精品一区二区精品久久久| 午夜福利,免费看| 精品少妇久久久久久888优播| 国产亚洲欧美在线一区二区| 老司机深夜福利视频在线观看 | 丝袜在线中文字幕| 两人在一起打扑克的视频| 国产精品久久久av美女十八| 男女高潮啪啪啪动态图| 日韩欧美一区二区三区在线观看 | 黄色 视频免费看| 日本av免费视频播放| 男人舔女人的私密视频| 日韩中文字幕视频在线看片| 国产精品久久久人人做人人爽| 欧美大码av| 欧美日本中文国产一区发布| 99九九在线精品视频| 国产福利在线免费观看视频| 建设人人有责人人尽责人人享有的| 香蕉国产在线看| 好男人电影高清在线观看| 国产精品亚洲av一区麻豆| 亚洲精品久久成人aⅴ小说| 一级黄色大片毛片| 超色免费av| 12—13女人毛片做爰片一| 极品人妻少妇av视频| 国产真人三级小视频在线观看| 欧美精品av麻豆av| 女人久久www免费人成看片| 亚洲专区字幕在线| 最新的欧美精品一区二区| 12—13女人毛片做爰片一| 国产成人免费观看mmmm| 欧美变态另类bdsm刘玥| 啦啦啦视频在线资源免费观看| 黄色视频,在线免费观看| 美女国产高潮福利片在线看| 香蕉国产在线看| 99久久人妻综合| 午夜激情久久久久久久| 在线精品无人区一区二区三| 美国免费a级毛片| a级毛片在线看网站| 国产成人欧美| 可以免费在线观看a视频的电影网站| av视频免费观看在线观看| 啦啦啦中文免费视频观看日本| 叶爱在线成人免费视频播放| 97精品久久久久久久久久精品| 美女主播在线视频| 免费黄频网站在线观看国产| av天堂久久9| 久久久国产精品麻豆| 亚洲va日本ⅴa欧美va伊人久久 | www日本在线高清视频| 免费高清在线观看视频在线观看| 99久久人妻综合| 日本猛色少妇xxxxx猛交久久| 十八禁网站免费在线| 国产亚洲精品一区二区www | 亚洲av欧美aⅴ国产| 蜜桃国产av成人99| 日韩大片免费观看网站| 久久亚洲国产成人精品v| 大片免费播放器 马上看| 交换朋友夫妻互换小说| 少妇粗大呻吟视频| 可以免费在线观看a视频的电影网站| kizo精华| 咕卡用的链子| 亚洲精品久久久久久婷婷小说| 好男人电影高清在线观看| 高清av免费在线| 国产在线视频一区二区| 色婷婷久久久亚洲欧美| 大香蕉久久成人网| 国产日韩欧美亚洲二区| 国产片内射在线| 国产免费现黄频在线看| 精品福利观看| 国产国语露脸激情在线看| 久久狼人影院| 啦啦啦中文免费视频观看日本| 成年美女黄网站色视频大全免费| 不卡一级毛片| 国产成人欧美在线观看 | 久久久精品免费免费高清| 丁香六月欧美| 精品一区二区三卡| 50天的宝宝边吃奶边哭怎么回事| 免费观看av网站的网址| 国产无遮挡羞羞视频在线观看| 女人精品久久久久毛片| 国产伦理片在线播放av一区| 午夜福利影视在线免费观看| 一区二区三区精品91| 新久久久久国产一级毛片| 搡老乐熟女国产| 精品久久久久久久毛片微露脸 | 男女之事视频高清在线观看| 欧美日韩黄片免| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品999| 人妻人人澡人人爽人人| 国产黄频视频在线观看| 欧美av亚洲av综合av国产av| 精品一区二区三区av网在线观看 | 日韩大码丰满熟妇| 男女午夜视频在线观看| 久久久久久久久免费视频了| 亚洲欧美精品自产自拍| 飞空精品影院首页| 首页视频小说图片口味搜索| 日韩视频在线欧美| 高清av免费在线| av视频免费观看在线观看| 欧美精品一区二区免费开放| 一级片'在线观看视频| 午夜福利乱码中文字幕| 久久综合国产亚洲精品| 女人久久www免费人成看片| 午夜福利影视在线免费观看| 欧美精品人与动牲交sv欧美| 久久人人爽人人片av| 欧美日韩一级在线毛片| 老熟女久久久| 自拍欧美九色日韩亚洲蝌蚪91| 久久毛片免费看一区二区三区| 我的亚洲天堂| 国产男女超爽视频在线观看| av不卡在线播放| 又大又爽又粗| 黄色 视频免费看| 99热网站在线观看| 天堂8中文在线网| 欧美黄色片欧美黄色片| 欧美性长视频在线观看| 久久中文看片网| 国产亚洲午夜精品一区二区久久| 国产老妇伦熟女老妇高清| 日韩欧美一区视频在线观看| 国产人伦9x9x在线观看| 亚洲第一青青草原| 50天的宝宝边吃奶边哭怎么回事| 啦啦啦 在线观看视频| 免费黄频网站在线观看国产| 国产区一区二久久| 免费少妇av软件| 国产亚洲午夜精品一区二区久久| 亚洲成人国产一区在线观看| 中文字幕av电影在线播放| 国产精品一区二区在线观看99| 精品久久久精品久久久| 丝袜美腿诱惑在线| 日韩免费高清中文字幕av| 久久99一区二区三区| 黑人巨大精品欧美一区二区mp4| 欧美成狂野欧美在线观看| www.自偷自拍.com| 大码成人一级视频| 国产又色又爽无遮挡免| 一级毛片精品| 久久人妻福利社区极品人妻图片| 手机成人av网站| 欧美成狂野欧美在线观看| 黑人巨大精品欧美一区二区蜜桃| 首页视频小说图片口味搜索| 日韩欧美一区视频在线观看| 国产成人av激情在线播放| 久久精品人人爽人人爽视色| 色婷婷久久久亚洲欧美| 欧美日韩成人在线一区二区| 狠狠精品人妻久久久久久综合| 日韩大片免费观看网站| 午夜福利一区二区在线看| 色视频在线一区二区三区| 99精国产麻豆久久婷婷| 精品福利观看| 欧美人与性动交α欧美精品济南到| 新久久久久国产一级毛片| √禁漫天堂资源中文www| 久久av网站| 黑人操中国人逼视频| 人成视频在线观看免费观看| 精品人妻一区二区三区麻豆| 国产99久久九九免费精品| 18在线观看网站|