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

    基于氘氚中子源硼中子俘獲治療的中子慢化整形研究

    2022-01-20 00:56:52朱益楠林作康盧林遠(yuǎn)郁海燕陳金根謝雷東
    核技術(shù) 2022年1期
    關(guān)鍵詞:中子通量快中子反射層

    朱益楠林作康盧林遠(yuǎn)郁海燕陳金根謝雷東

    1(中國(guó)科學(xué)院上海應(yīng)用物理研究所上海201800)

    2(中國(guó)科學(xué)院大學(xué)北京100049)

    硼中子俘獲治療(Boron Neutron Capture Therapy,BNCT)是一種用于治療惡性腫瘤的二元靶向放療技術(shù),具有很大的發(fā)展前景。其主要原理為將含有10B的硼攜帶劑注射至患者體內(nèi),使之與腫瘤細(xì)胞特異性結(jié)合,隨后用熱中子束或超熱中子束在病患區(qū)進(jìn)行照射,在腫瘤細(xì)胞內(nèi)發(fā)生10B(n,α)7Li反應(yīng),產(chǎn)生的α粒子和7Li離子具有很高的線性能量傳輸(Linear Energy Transfer,LET)值,在體內(nèi)的射程分別為9 μm和5 μm,相當(dāng)于一個(gè)細(xì)胞的直徑尺度,對(duì)周?chē)<?xì)胞的損傷可忽略不計(jì)[1]。

    世界上首次BNCT臨床實(shí)驗(yàn)是1951年由美國(guó)布魯克海文國(guó)家實(shí)驗(yàn)室同馬薩諸塞州總醫(yī)院通過(guò)核反應(yīng)堆導(dǎo)出的熱中子束進(jìn)行的[2]。由于核反應(yīng)堆多建于偏遠(yuǎn)地區(qū),運(yùn)營(yíng)程序復(fù)雜,潛在危險(xiǎn)性大,且體積大、維護(hù)成本高,多數(shù)早期BNCT反應(yīng)堆因經(jīng)費(fèi)與政策等因素相繼關(guān)停,反應(yīng)堆BNCT中子源逐漸退出歷史舞臺(tái)。

    近年來(lái),基于加速器驅(qū)動(dòng)的BNCT中子源正在火熱地發(fā)展起來(lái),其原理是通過(guò)經(jīng)高能加速器加速的質(zhì)子轟擊靶材料,產(chǎn)生高通量的中子。常見(jiàn)的反應(yīng)如:7Li(p,n)7Be、9Be(p,n)9B。與反應(yīng)堆中子源相比,加速器中子源成本較低、設(shè)備運(yùn)行及操作較反應(yīng)堆簡(jiǎn)便、固有安全性高、占地面積小、公眾接受程度高、適合建在人群密集的醫(yī)院等優(yōu)勢(shì)[3]。日本在BNCT誕生初期就開(kāi)展了針對(duì)惡性腦膠質(zhì)瘤和惡性黑色素瘤的臨床實(shí)驗(yàn)及相應(yīng)的基礎(chǔ)性研究,在此方面具有較豐富的經(jīng)驗(yàn),與此同時(shí),基于加速器的BNCT治療系統(tǒng)也在不斷開(kāi)發(fā)中[4-7]。2020年3月,住友重工的BNCT治療系統(tǒng)獲得日本厚生勞動(dòng)省批準(zhǔn),正式成為癌癥治療設(shè)備,同年5月,首位患者在日本南東北綜合醫(yī)院BNCT研究中心接受治療,日本也因此成為了世界上首個(gè)推動(dòng)BNCT走向商業(yè)化的國(guó)家。

    對(duì)于加速器驅(qū)動(dòng)BNCT所用靶材,不同材料具有各自的優(yōu)缺點(diǎn)。如鋰靶的熔點(diǎn)較低,散熱性較差,在質(zhì)子轟擊靶的過(guò)程中難以保持其原有形狀[8]。鈹靶有劇毒,更換靶材難度大,以及發(fā)生(p,n)反應(yīng)的閾能高,需要使用高能加速器對(duì)質(zhì)子加速至2.5~4 MeV才能發(fā)生反應(yīng)。因此,針對(duì)高性能BNCT中子源研發(fā)所需的一些關(guān)鍵技術(shù)仍在探索中。

    自1998年以來(lái),意大利熱那亞大學(xué)開(kāi)始對(duì)氘氚中子源應(yīng)用于BNCT的項(xiàng)目進(jìn)行了研究,同一時(shí)期,美國(guó)勞倫斯伯克利國(guó)家實(shí)驗(yàn)室也開(kāi)始對(duì)緊湊型、低功率、高通量的氘氚中子源進(jìn)行研發(fā)[9]。這種基于氘氚反應(yīng)的緊湊型中子源,相較于傳統(tǒng)的高能加速器中子源,具有更高的安全性、更高的中子產(chǎn)額、更低的入射粒子能量、更低的成本以及更小的尺寸,有望成為未來(lái)可供選擇的BNCT中子源之一[10]。然而,由于氘氚中子源產(chǎn)生的中子能量為14.1 MeV,不能直接用于BNCT,因此,需要設(shè)計(jì)一種束流整形組件(Beam Shaping Assembly,BSA),將產(chǎn)生的快中子慢化至BNCT所需的超熱中子能區(qū)。此外,中子束的輸出端參數(shù)應(yīng)滿足國(guó)際原子能機(jī)構(gòu)(International Atomic Energy Agency,IAEA)所給出的建議值(表1)。

    表1 IAEA推薦的治療用中子束指標(biāo)[11]Table 1 Targeted values of therapeutic neutron beam recommended by IAEA[11]

    本文主要針對(duì)緊湊型氘氚中子源,設(shè)計(jì)一種基于14.1 MeV快中子下的BSA,通過(guò)蒙特卡羅程序(Monte Carlo N Particle Transport Code,version 5,MCNP5)模擬計(jì)算,評(píng)估中子束輸出端的參數(shù),對(duì)氘氚中子源用于BNCT的可行性進(jìn)行研究。

    1 材料與方法

    1.1 D-T中子源

    緊湊型氘氚中子源是一種利用射頻放電產(chǎn)生氘離子和氚離子,并經(jīng)過(guò)100~300 keV的電場(chǎng)加速后轟擊至含鈦、鋁氚靶,在靶上發(fā)生氘氚反應(yīng)產(chǎn)生中子的新型中子源,其核反應(yīng)式:

    式中:Q=17.586 MeV,產(chǎn)生的單能中子能量為14.1 MeV。由于其能量單一、準(zhǔn)直性好,中子產(chǎn)額可達(dá)到1012~1014n·s-1,是合適的BNCT中子源[12-13]。但由于氘氚反應(yīng)產(chǎn)生的中子相對(duì)于其他中子源來(lái)說(shuō),中子能量高,慢化難度大,在確保將快中子慢化至超熱中子能區(qū)的同時(shí)還要保證超熱中子的通量足夠大,這成為了目前為止氘氚中子源得到實(shí)際應(yīng)用于BNCT的重要挑戰(zhàn)。

    本文利用緊湊型氘氚中子源開(kāi)展對(duì)束流整形組件的物理設(shè)計(jì)及分析,探討將該中子源用作BNCT中子源的可行性,為未來(lái)BNCT中子源候選的緊湊型氘氚中子源提供理論基礎(chǔ)。

    1.2 蒙特卡羅模擬及計(jì)算模型

    本文的計(jì)算使用蒙特卡羅模擬程序MCNP5是由美國(guó)洛斯阿拉莫斯國(guó)家實(shí)驗(yàn)室開(kāi)發(fā)的一套用于模擬中子、電子及光子或其耦合態(tài)在物質(zhì)中輸運(yùn)過(guò)程的計(jì)算程序。BSA的作用是將高能中子慢化至可用于進(jìn)行治療的超熱中子,因此,必須降低快中子、熱中子的份額以及減少γ射線的污染。本文設(shè)計(jì)的BSA主要由快中子倍增層、慢化層、反射層、熱中子吸收層以及γ屏蔽層等五部分組成,整體形狀為圓柱體,BSA模型如圖1所示。

    圖1 束流整形體(BSA)示意圖Fig.1 Diagram of beam shaping assembly

    在模擬計(jì)算過(guò)程中,中子源定義為半徑為0.1 cm的單向圓形面源,方向?yàn)閤軸正向,中子能量定義為14.1 MeV,強(qiáng)度設(shè)為1.0×1014n·s-1,通過(guò)調(diào)用JEFF3.1核數(shù)據(jù)庫(kù),使用F1、F2、F4及F6卡對(duì)模擬過(guò)程中的中子通量和劑量率等有關(guān)參數(shù)進(jìn)行計(jì)算,為了使整體計(jì)算結(jié)果的相對(duì)誤差小于1%,運(yùn)行源粒子數(shù)大于108個(gè)。

    2 結(jié)果與討論

    2.1 中子倍增層

    中子在經(jīng)過(guò)BSA的不同種材料時(shí),由于散射和俘獲吸收,不可避免地會(huì)損失一部分中子。為了使經(jīng)BSA慢化后的超熱中子通量滿足IAEA的建議值,設(shè)計(jì)一層與中子源相鄰的中子倍增層,以增加入射中子數(shù),從而部分補(bǔ)償在慢化期間損失的中子。D-T中子源能量高達(dá)14.1 MeV,在此能量下可利用多種材料進(jìn)行中子倍增。

    基于評(píng)價(jià)核截面數(shù)據(jù)庫(kù)(Evaluated Nuclear Data File,ENDF),選擇以下7種在14.1 MeV能量下具有相對(duì)較高的(n,2n)或(n,3n)反應(yīng)截面的材料,即鈹、鉍、釷、鉛、238U、天然鈾(238U豐度為99.27%)以及20%富集度的濃縮鈾,幾何形狀為球體,作為中子倍增層候選材料。通過(guò)統(tǒng)計(jì)經(jīng)由倍增材料后的中子數(shù),結(jié)果如圖2所示。圖2中針對(duì)天然鈾、238U和20%富集度的濃縮鈾對(duì)比可知,隨著235U比例的增加,中子倍增效果也會(huì)明顯增加。雖然20%富集度的濃縮鈾的中子倍增性能遠(yuǎn)優(yōu)于天然鈾,但是相較于20%富集度的濃縮鈾而言,天然鈾的成本更低廉,加工更方便,可獲取性更高,更適用于未來(lái)醫(yī)院等場(chǎng)所的環(huán)境。因此,本文最終選擇半徑為14 cm的天然鈾球作為中子倍增層。

    圖2 不同倍增材料下的中子數(shù)與倍增球半徑的關(guān)系Fig.2 Relationship between neutron number and the radius of neutron multiplier consisted of different materials

    圖3為中子經(jīng)過(guò)半徑為14 cm天然鈾球倍增材料后相應(yīng)的能譜。由圖3可知,快中子在經(jīng)過(guò)倍增后發(fā)生核反應(yīng),能譜會(huì)往左移動(dòng),但仍有一部分14.1 MeV的入射中子沒(méi)有參與核反應(yīng),直接從倍增層透射出去,透射率為12.4%。結(jié)合圖2分析可知,增加了中子倍增層后,不僅使得中子數(shù)得到增加,并且能夠降低中子能量,起到一定的慢化作用,使得D-T中子的慢化難度在一定程度上得到減輕。

    圖3 經(jīng)14 cm半徑天然鈾球倍增后的中子能譜Fig.3 Neutron spectrum after multiplying with 14 cm radius of natural uranium

    2.2 慢化層

    在針對(duì)慢化層的設(shè)計(jì)時(shí),候選材料應(yīng)滿足在快中子能區(qū)具有較高的彈性散射截面和較高的非彈性散射截面,以及較低的吸收截面。同時(shí),所選材料需在超熱中子能區(qū)具有較低的散射截面和吸收截面,以減少超熱中子在到達(dá)患者體表之前的損失,同時(shí)也應(yīng)考慮慢化材料在吸收中子后導(dǎo)致的γ射線污染。

    本文設(shè)計(jì)的慢化層為多層慢化組合結(jié)構(gòu),對(duì)于上一小節(jié)中經(jīng)過(guò)倍增層的中子來(lái)說(shuō),能量大于10 MeV的快中子仍占有一定的比重,平均中子能量仍然較高,因此,需要尋找一種在快中子能區(qū)散射截面較大的材料作為第一慢化材料。在快中子能區(qū),由于中子能量已達(dá)到發(fā)生非彈性散射的閾能,在此能區(qū)重核的慢化效果優(yōu)于輕核[14],如:Pb、Bi等元素,這些元素不僅有較大的彈性散射截面和非彈性散射截面,而且吸收截面非常小,因此,選用含Pb、Bi等重元素的材料作為第一慢化材料,再選用利于調(diào)節(jié)中子能量至超熱能區(qū)的較輕元素,如:F、Mg、Al、Ti等元素作為第二慢化材料。而快中子在含氟材料中慢化至超熱中子的宏觀散射截面與宏觀吸收截面的比值最大,因而具有較好的慢化性能[15]。因此,本文最終選用PbF2、PbF4、BiF3作為第一慢化材料的候選材料,MgF2、AlF3、TiF3作為第二慢化材料的候選材料,根據(jù)現(xiàn)有的BSA設(shè)計(jì)經(jīng)驗(yàn),為控制中子經(jīng)過(guò)慢化體的損失及慢化體尺寸,將兩層慢化體的厚度總長(zhǎng)設(shè)置為60 cm,半徑設(shè)置為25 cm,調(diào)節(jié)慢化體長(zhǎng)度組合進(jìn)行模擬,篩選出慢化效果最佳的組合,圖4和圖5分別給出了在兩層慢化體總長(zhǎng)度為60 cm時(shí),不同組合下的中子慢化效果。

    在本文所設(shè)計(jì)的慢化體中,由于D-T中子的能量較高,因此在前兩層慢化體的設(shè)計(jì)中,首要任務(wù)是將較高能的快中子能量降低,設(shè)計(jì)過(guò)程中主要以快中子占總中子的份額以及超熱中子與快中子的比值為依據(jù),盡量使得較多的快中子能慢化下來(lái)。

    由圖4可知,當(dāng)選擇25 cm厚的BiF3與35 cm厚的MgF2以及50 cm厚的BiF3與10 cm厚的TiF3這兩種組合慢化構(gòu)型時(shí),快中子通量與總中子通量的比值均達(dá)到最小值,且兩者的數(shù)值較為接近,另結(jié)合圖5分析可知,在上述兩種組合慢化構(gòu)型下,超熱中子通量與快中子通量的比值也達(dá)到最大值,是較為理想的兩種慢化設(shè)計(jì)方案。為了進(jìn)一步比較上述兩種組合慢化構(gòu)型,分別模擬計(jì)算了中子經(jīng)過(guò)這兩種組合慢化層后的能譜,結(jié)果如圖6所示。

    圖4 不同組合情況下的快中子通量與總中子通量的比值Fig.4 The ratio of fast neutron flux to total neutron flux at different combinations

    圖5 不同組合情況下的超熱中子通量與快中子通量的比值Fig.5 The ratio of epithermal neutron flux to fast neutron flux at different combinations

    圖6 中子經(jīng)過(guò)兩種不同組合慢化層后的能譜Fig.6 The spectra of neutrons passing through two different combinations of moderating layers

    由 圖6可 知,采 用50 cm厚 的BiF3與10 cm厚TiF3這一組合慢化結(jié)構(gòu)能使中子能譜更好地集中在超熱能區(qū)的范圍之內(nèi),而選擇25 cm厚的BiF3與35 cm厚的MgF2會(huì)使得中子能譜過(guò)熱,導(dǎo)致能譜在熱中子能區(qū)有一個(gè)小峰出現(xiàn),而在BNCT的過(guò)程中并不希望有過(guò)多的熱中子。綜合以上考慮,本文選擇50 cm厚的BiF3與10 cm厚的TiF3的組合結(jié)構(gòu)作為慢化層。

    2.3 補(bǔ)充慢化層兼中子引出通道

    中子引出通道的主要作用是將高能中子經(jīng)慢化后輸出的超熱中子束精確聚焦至患者腫瘤所在位置,以便于使照射野達(dá)到最佳尺寸,從而達(dá)到更好的治療效果。另外,考慮到高能中子經(jīng)過(guò)此前所采用的組合慢化結(jié)構(gòu)仍未被完全慢化至超熱中子能段,因此,補(bǔ)充慢化層選用Al2O3、AlF3、Fluental[16]、MgO這4種具有進(jìn)一步慢化快中子能力的材料,并進(jìn)行模擬計(jì)算。其中,F(xiàn)luental是一種由芬蘭國(guó)家技術(shù)研究中心開(kāi)發(fā)的專利材料,主要成分為69%質(zhì)量比的AlF3、30%的Al和1%的LiF。在本文所設(shè)計(jì)的中子引出通道中,超熱中子束出口處直徑設(shè)置為15 cm,以滿足大部分腫瘤治療的情況。補(bǔ)充慢化層的初始厚度設(shè)置為10 cm,并按步長(zhǎng)2 cm遞增至24 cm,圖7和圖8分別為不同厚度下,中子經(jīng)補(bǔ)充慢化層后出口端的相關(guān)參數(shù)。

    由圖8可知,當(dāng)MgO在19 cm處時(shí),超熱中子與快中子的比值最高,但是考慮到MgO在起到慢化作用的同時(shí),吸收了大量超熱中子,導(dǎo)致超熱中子數(shù)降低,同時(shí)會(huì)產(chǎn)生大量的γ放射性污染。Al2O3有著與MgO相似的結(jié)果,而氟化物則表現(xiàn)出更優(yōu)異的慢化性能,這與前文提到過(guò)的含氟材料的快中子慢化至超熱中子的能力相對(duì)更好這一結(jié)論相符合。當(dāng)AlF3厚度處于17 cm與19 cm時(shí),計(jì)算得到的超熱中子通量與熱中子通量之比接近一致,隨著AlF3厚度的進(jìn)一步增加,超熱中子數(shù)會(huì)隨之減少。結(jié)合圖7的分析結(jié)果,選用17 cm厚的AlF3作為補(bǔ)充慢化材料是較合理的選擇。

    圖7 幾種不同厚度補(bǔ)充慢化層下的超熱中子通量Fig.7 Epithermal neutron flux under several supplementary moderators with different thicknesses

    圖8 幾種不同厚度補(bǔ)充慢化層下的超熱中子通量與快中子通量的比值Fig.8 Ratio of epithermal neutron flux to fast neutron flux under several supplementary moderators with different thicknesses

    2.4 反射層

    反射層應(yīng)選擇在超熱中子能區(qū)具有高彈性散射截面和低吸收截面的材料,其包覆在整個(gè)慢化體的最外層,使一部分散射出去的中子能夠反射回來(lái),重新進(jìn)入慢化層,從而可以起到提高中子經(jīng)濟(jì)性的作用。本文對(duì)Pb和Bi兩種效果相對(duì)較好的反射層材料進(jìn)行不同厚度下的模擬計(jì)算,計(jì)算結(jié)果如圖9和圖10所示。

    通過(guò)圖9和圖10可以明顯看出,在增加了反射層之后,超熱中子數(shù)明顯增加,超熱中子通量與快中子通量的比值也隨之增加,且Pb作為反射層的性能要優(yōu)于Bi。由圖10可得,在Pb反射層厚度大于10 cm后,超熱中子通量與快中子通量的比值增加趨于平緩,且超熱中子通量與快中子通量的比值已經(jīng)大于20,滿足了IAEA的官方建議值??紤]到BSA的整體質(zhì)量和尺寸大小,可選用厚度為10 cm的Pb作為反射層。

    圖9 不同厚度鉍、鉛反射層下的超熱中子通量Fig.9 Epithermal neutron flux under different thicknesses of Bi and Pb reflectors

    圖10 不同厚度鉍、鉛反射層下的超熱中子通量與快中子通量的比值Fig.10 The ratio of epithermal neutron flux to fast neutron flux under different thicknesses of Bi and Pb reflectors

    2.5 熱中子吸收層及γ屏蔽層

    中子在慢化過(guò)程中,不可避免會(huì)產(chǎn)生熱中子,而在BNCT治療過(guò)程中并不希望有過(guò)多的熱中子照射到患者體表。因此,中子束照射至病患處之前,應(yīng)盡可能把多余的熱中子吸收掉。鎘(Cd)的天然同位素中的113Cd自然豐度為12.22%,在熱中子能區(qū)具有極大的吸收截面,僅需0.2 mm的鎘(Cd)就可將超熱中子通量與熱中子通量的比值提高到100以上,滿足IAEA的建議值[17]。但鎘(Cd)在吸收了熱中子后會(huì)產(chǎn)生γ射線,為了屏蔽掉這部分γ射線,本文選用常用的γ過(guò)濾材料Pb放置于熱中子吸收層之后,經(jīng)模擬計(jì)算得到的快中子劑量率與γ劑量率隨Pb屏蔽層的厚度變化的曲線,結(jié)果如圖11和圖12所示。由圖11可以看出,當(dāng)Pb屏蔽層從0.5 cm遞增到5 cm時(shí),快中子劑量率與超熱中子通量的比值始終小于2×10-13Gy·cm2,可滿足IAEA的建議標(biāo)準(zhǔn)。結(jié)合圖12,當(dāng)Pb過(guò)濾層大于3.5 cm時(shí),γ劑量率與超熱中子通量的比值小于2×10-13Gy·cm2,亦符合IAEA的建議標(biāo)準(zhǔn)。因此,考慮整個(gè)BSA的大小以及質(zhì)量,本文最終選擇3.5 cm厚的Pb作為γ屏蔽層。

    圖11 不同厚度鉛屏蔽層下的快中子劑量率與超熱中子通量的比值Fig.11 Ratio of fast neutron dose rate to epithermal neutron flux under different thicknesses of Pb filter

    圖12 不同厚度鉛屏蔽層下的γ劑量率與超熱中子通量的比值Fig.12 Ratio of gamma dose rate to epithermal neutron flux under different thicknesses of Pb filter

    3 結(jié)語(yǔ)

    綜上所述,本文選用了半徑為14 cm的天然鈾球作為中子倍增層,50 cm厚的BiF3和10 cm厚的TiF3作為組合慢化層,17 cm厚的AlF3作為補(bǔ)充慢化層兼中子引出通道,0.2 mm厚的Cd作為熱中子吸收層,3.5 cm厚的Pb作為γ屏蔽層,以及10 cm厚的Pb作為反射層,整體的BSA模型如圖13所示。另外,計(jì)算了中子束出口處的中子單向流量與通量的比值,即J/Φ,該值可以很好反映出中子束的方向性。強(qiáng)度為1.0×1014n·s-1的D-T中子源經(jīng)過(guò)該BSA慢化,最終得到的中子能譜和空氣端出口處參數(shù)如圖14和表2所示。

    表2 D-T中子經(jīng)BSA慢化整形后得到的空氣端出口參數(shù)Table 2 The air exit parameters of D-T neutron after BSA moderation and shaping

    圖13 BSA模型圖Fig.13 Diagram of BSA model

    圖14 D-T中子經(jīng)BSA慢化整形后的中子能譜圖Fig.14 Neutron energy spectrum of D-T neutron moderated and shaped by BSA

    因此,D-T中子在經(jīng)過(guò)優(yōu)化的BSA之后所得到的參數(shù)均能很好地滿足IAEA的建議值,說(shuō)明D-T中子源在應(yīng)用于BNCT上具有可行性。

    作者貢獻(xiàn)聲明朱益楠負(fù)責(zé)文章的起草及最終版本修訂;盧林遠(yuǎn)負(fù)責(zé)資料的收集;陳金根負(fù)責(zé)論文的修改;林作康、郁海燕負(fù)責(zé)資料的搜集及整理;謝雷東負(fù)責(zé)研究的提出及設(shè)計(jì)。

    猜你喜歡
    中子通量快中子反射層
    反射層對(duì)NaI(Tl)閃爍體探測(cè)器探測(cè)效率的影響
    科技視界(2021年10期)2021-05-20 01:46:42
    基于有限元某壓水堆反射層組件模態(tài)分析的建模方法
    四川建筑(2020年3期)2020-07-18 01:56:24
    基于協(xié)同進(jìn)化的航空高度單粒子翻轉(zhuǎn)故障生成方法研究
    導(dǎo)向管快中子注量計(jì)算方法研究*
    49—2堆孔道快中子注量率的測(cè)定
    商情(2016年44期)2017-03-05 02:46:22
    隨鉆雙感應(yīng)測(cè)井儀背景影響研究
    燃料元件中235U豐度的快中子誘發(fā)裂變測(cè)量方法研究
    修正快中子通量以提高碳氧測(cè)量精度的研究
    某型“三代”核電機(jī)組與M310機(jī)組堆芯測(cè)量系統(tǒng)
    煤炭?jī)?nèi)中子通量與元素含量關(guān)系
    色尼玛亚洲综合影院| 人人妻,人人澡人人爽秒播| 亚洲第一青青草原| 国产亚洲av嫩草精品影院| 国产v大片淫在线免费观看| 中文字幕久久专区| 大型黄色视频在线免费观看| 国产精品二区激情视频| 亚洲成av人片免费观看| 高清在线国产一区| 午夜两性在线视频| 宅男免费午夜| 桃红色精品国产亚洲av| 我的亚洲天堂| 在线观看www视频免费| 草草在线视频免费看| 两个人视频免费观看高清| 99久久国产精品久久久| 91在线观看av| 国产亚洲精品第一综合不卡| 国产伦人伦偷精品视频| 老司机福利观看| 一本综合久久免费| 97人妻精品一区二区三区麻豆 | 久久人妻av系列| 亚洲欧美精品综合一区二区三区| 久久99热这里只有精品18| 亚洲七黄色美女视频| 国产区一区二久久| 不卡av一区二区三区| av中文乱码字幕在线| 久久香蕉国产精品| 国产日本99.免费观看| 久久久久久久久中文| 日韩免费av在线播放| 老汉色av国产亚洲站长工具| 亚洲成人久久爱视频| 欧美一级a爱片免费观看看 | 欧美日韩中文字幕国产精品一区二区三区| x7x7x7水蜜桃| 搞女人的毛片| 久久精品aⅴ一区二区三区四区| 国产成人精品无人区| 成年女人毛片免费观看观看9| 国产1区2区3区精品| 老司机靠b影院| 精品乱码久久久久久99久播| 欧美日韩一级在线毛片| 国产亚洲精品综合一区在线观看 | 黄色视频,在线免费观看| 三级毛片av免费| 国产激情欧美一区二区| 一进一出抽搐gif免费好疼| 亚洲熟妇中文字幕五十中出| 热99re8久久精品国产| 国产1区2区3区精品| 91字幕亚洲| 国产黄a三级三级三级人| 久久午夜亚洲精品久久| 一级毛片高清免费大全| 国产亚洲精品第一综合不卡| 亚洲一码二码三码区别大吗| 色av中文字幕| 精品少妇一区二区三区视频日本电影| 一二三四社区在线视频社区8| 一本一本综合久久| 人妻丰满熟妇av一区二区三区| 亚洲精品中文字幕在线视频| 久久久久久大精品| 两人在一起打扑克的视频| 久久精品国产综合久久久| 国产精品久久久久久亚洲av鲁大| 18禁美女被吸乳视频| 久久久久国产一级毛片高清牌| 99国产精品一区二区三区| 成人三级黄色视频| 精品久久久久久久末码| 久久青草综合色| 国产精品综合久久久久久久免费| 成人免费观看视频高清| 精品欧美一区二区三区在线| 国产亚洲精品久久久久久毛片| 国产一卡二卡三卡精品| 男人的好看免费观看在线视频 | 999精品在线视频| 欧美激情极品国产一区二区三区| www国产在线视频色| 欧美 亚洲 国产 日韩一| 丁香欧美五月| 18禁裸乳无遮挡免费网站照片 | 一区二区日韩欧美中文字幕| 精品欧美国产一区二区三| 一级毛片女人18水好多| 男人舔女人的私密视频| 欧美日韩乱码在线| 最新美女视频免费是黄的| 亚洲电影在线观看av| 天堂动漫精品| 欧美丝袜亚洲另类 | 色综合婷婷激情| 嫩草影视91久久| 亚洲人成77777在线视频| 最新在线观看一区二区三区| 一区二区三区精品91| 国产成人av激情在线播放| 99久久久亚洲精品蜜臀av| 国产精品九九99| 最近在线观看免费完整版| 看片在线看免费视频| 国产精品影院久久| 日本 欧美在线| 国产高清videossex| 熟女电影av网| 露出奶头的视频| 久久久国产成人精品二区| 久久久久国产精品人妻aⅴ院| 一区福利在线观看| 99精品久久久久人妻精品| 可以在线观看毛片的网站| 此物有八面人人有两片| 亚洲aⅴ乱码一区二区在线播放 | 男人舔女人的私密视频| 亚洲黑人精品在线| 成人欧美大片| 精品久久久久久,| 99国产精品一区二区蜜桃av| 久久久久精品国产欧美久久久| 日本黄色视频三级网站网址| 淫秽高清视频在线观看| 亚洲男人的天堂狠狠| 亚洲欧美精品综合一区二区三区| 欧美激情 高清一区二区三区| 欧美色视频一区免费| 免费一级毛片在线播放高清视频| 日本五十路高清| 久久久久久久精品吃奶| 可以在线观看毛片的网站| 男女午夜视频在线观看| 欧美在线一区亚洲| 国产激情偷乱视频一区二区| 国产精品久久久久久亚洲av鲁大| 国产精品 国内视频| 黄色a级毛片大全视频| 好男人电影高清在线观看| 日韩欧美免费精品| 高清在线国产一区| 男女午夜视频在线观看| 大香蕉久久成人网| 国产成人精品久久二区二区免费| 欧美av亚洲av综合av国产av| 一本一本综合久久| 婷婷亚洲欧美| 亚洲精品国产精品久久久不卡| 成年免费大片在线观看| 久久久久久免费高清国产稀缺| 成人免费观看视频高清| 亚洲欧洲精品一区二区精品久久久| 悠悠久久av| 免费观看人在逋| www.熟女人妻精品国产| 亚洲av成人不卡在线观看播放网| 国产伦人伦偷精品视频| 国内精品久久久久精免费| 亚洲av成人不卡在线观看播放网| 成人午夜高清在线视频 | av在线天堂中文字幕| 色综合欧美亚洲国产小说| 中文亚洲av片在线观看爽| 国产午夜精品久久久久久| 亚洲成国产人片在线观看| 一级黄色大片毛片| 天天添夜夜摸| 国产又黄又爽又无遮挡在线| 亚洲最大成人中文| 久久精品aⅴ一区二区三区四区| 亚洲精品久久国产高清桃花| 免费看日本二区| 又黄又爽又免费观看的视频| 亚洲久久久国产精品| 国产视频内射| 精品少妇一区二区三区视频日本电影| 欧美日韩亚洲综合一区二区三区_| 国产久久久一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 亚洲熟女毛片儿| 国产亚洲欧美在线一区二区| 亚洲av电影不卡..在线观看| 在线观看www视频免费| 欧美性猛交╳xxx乱大交人| 国产成人啪精品午夜网站| 成人av一区二区三区在线看| 性色av乱码一区二区三区2| 久久久久精品国产欧美久久久| 男人舔女人的私密视频| 嫁个100分男人电影在线观看| 精品高清国产在线一区| 国产精品免费一区二区三区在线| 中文字幕精品亚洲无线码一区 | 日韩欧美免费精品| 可以在线观看毛片的网站| 丁香欧美五月| 久久久精品国产亚洲av高清涩受| 亚洲av五月六月丁香网| 天天添夜夜摸| 久久人妻av系列| 久久人人精品亚洲av| 十八禁人妻一区二区| 国产成人精品久久二区二区免费| 亚洲精品久久成人aⅴ小说| 一级作爱视频免费观看| 亚洲专区国产一区二区| 国产亚洲精品久久久久5区| 欧美久久黑人一区二区| 男女下面进入的视频免费午夜 | 一边摸一边做爽爽视频免费| 亚洲狠狠婷婷综合久久图片| 高清在线国产一区| 2021天堂中文幕一二区在线观 | 久久欧美精品欧美久久欧美| 日本 欧美在线| 看黄色毛片网站| 婷婷六月久久综合丁香| 老司机午夜福利在线观看视频| 亚洲精品国产区一区二| 中亚洲国语对白在线视频| 看片在线看免费视频| 男女下面进入的视频免费午夜 | 亚洲中文字幕一区二区三区有码在线看 | 色在线成人网| tocl精华| 国产乱人伦免费视频| 色综合站精品国产| 精品一区二区三区av网在线观看| 97人妻精品一区二区三区麻豆 | 久久香蕉精品热| 国产精品av久久久久免费| 九色国产91popny在线| 国产成人系列免费观看| 男人操女人黄网站| 又紧又爽又黄一区二区| 国内毛片毛片毛片毛片毛片| 男女那种视频在线观看| 国产高清激情床上av| 久久久精品国产亚洲av高清涩受| 一边摸一边做爽爽视频免费| 日韩欧美免费精品| 人人妻人人看人人澡| 免费av毛片视频| 国产精品亚洲av一区麻豆| 身体一侧抽搐| 日本五十路高清| 香蕉国产在线看| 欧美国产精品va在线观看不卡| 日韩三级视频一区二区三区| 法律面前人人平等表现在哪些方面| 男人舔奶头视频| 嫁个100分男人电影在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 天天添夜夜摸| 真人做人爱边吃奶动态| 长腿黑丝高跟| 久久午夜综合久久蜜桃| 日韩欧美一区视频在线观看| 日本撒尿小便嘘嘘汇集6| 国产一区二区三区视频了| 啪啪无遮挡十八禁网站| 老司机深夜福利视频在线观看| 欧美激情极品国产一区二区三区| 国产高清videossex| 成人三级黄色视频| 国产区一区二久久| 人人妻,人人澡人人爽秒播| 香蕉国产在线看| 手机成人av网站| 欧美性猛交黑人性爽| xxx96com| 男人舔女人下体高潮全视频| 中文字幕精品亚洲无线码一区 | 18美女黄网站色大片免费观看| netflix在线观看网站| 99热只有精品国产| 国产激情久久老熟女| av天堂在线播放| 久久午夜亚洲精品久久| 欧美三级亚洲精品| 久久精品成人免费网站| 亚洲av五月六月丁香网| 99国产精品一区二区蜜桃av| 韩国av一区二区三区四区| 哪里可以看免费的av片| 天堂影院成人在线观看| 久久婷婷人人爽人人干人人爱| e午夜精品久久久久久久| 每晚都被弄得嗷嗷叫到高潮| 在线观看午夜福利视频| 激情在线观看视频在线高清| 久久人妻福利社区极品人妻图片| 国产精品电影一区二区三区| 丰满的人妻完整版| 国产av一区在线观看免费| 巨乳人妻的诱惑在线观看| 在线看三级毛片| 一区二区三区高清视频在线| 亚洲五月婷婷丁香| 熟女电影av网| 一本久久中文字幕| 亚洲国产高清在线一区二区三 | 亚洲自偷自拍图片 自拍| 亚洲av日韩精品久久久久久密| 色综合欧美亚洲国产小说| 亚洲真实伦在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲一区中文字幕在线| 久久精品国产亚洲av高清一级| 亚洲国产日韩欧美精品在线观看 | 国产av一区在线观看免费| 啦啦啦免费观看视频1| 桃色一区二区三区在线观看| 伊人久久大香线蕉亚洲五| 一本一本综合久久| 色哟哟哟哟哟哟| 丝袜人妻中文字幕| 黑人巨大精品欧美一区二区mp4| 欧美zozozo另类| 亚洲专区字幕在线| 国产亚洲精品久久久久5区| 精品卡一卡二卡四卡免费| 欧美激情极品国产一区二区三区| 我的亚洲天堂| 在线天堂中文资源库| 成人18禁高潮啪啪吃奶动态图| 国产高清激情床上av| 欧美成人免费av一区二区三区| 村上凉子中文字幕在线| 亚洲成国产人片在线观看| 99国产极品粉嫩在线观看| 国产国语露脸激情在线看| 成人av一区二区三区在线看| 精品欧美国产一区二区三| 亚洲电影在线观看av| 69av精品久久久久久| 超碰成人久久| 国产成人精品无人区| 欧美 亚洲 国产 日韩一| 香蕉丝袜av| 色综合婷婷激情| videosex国产| 国产亚洲精品av在线| 亚洲中文av在线| 久久青草综合色| 免费高清视频大片| 亚洲精品av麻豆狂野| 十八禁网站免费在线| 亚洲在线自拍视频| 变态另类丝袜制服| 丁香六月欧美| 91成年电影在线观看| 亚洲精品av麻豆狂野| 日韩三级视频一区二区三区| 99精品在免费线老司机午夜| 中文字幕人妻丝袜一区二区| 日本一本二区三区精品| 美女高潮喷水抽搐中文字幕| 中文在线观看免费www的网站 | 男人的好看免费观看在线视频 | 亚洲久久久国产精品| 91老司机精品| 亚洲精品粉嫩美女一区| 真人做人爱边吃奶动态| 黄色视频,在线免费观看| 91九色精品人成在线观看| 亚洲狠狠婷婷综合久久图片| 欧美日本视频| 精品免费久久久久久久清纯| 亚洲最大成人中文| 亚洲五月天丁香| 国语自产精品视频在线第100页| 免费搜索国产男女视频| 变态另类成人亚洲欧美熟女| 一边摸一边抽搐一进一小说| 日本成人三级电影网站| 在线观看免费日韩欧美大片| 妹子高潮喷水视频| 一进一出抽搐动态| 亚洲美女黄片视频| 免费看日本二区| 久久国产精品男人的天堂亚洲| 久久精品国产99精品国产亚洲性色| 成人三级黄色视频| 法律面前人人平等表现在哪些方面| 欧美黑人精品巨大| 国产又黄又爽又无遮挡在线| 男女床上黄色一级片免费看| 欧美最黄视频在线播放免费| 女性生殖器流出的白浆| 搡老熟女国产l中国老女人| 中文字幕另类日韩欧美亚洲嫩草| 不卡av一区二区三区| 国产成人精品久久二区二区91| 美女国产高潮福利片在线看| 亚洲黑人精品在线| 黄网站色视频无遮挡免费观看| 国产色视频综合| 亚洲国产欧美日韩在线播放| bbb黄色大片| av天堂在线播放| 亚洲成国产人片在线观看| www日本黄色视频网| 老司机福利观看| 一区二区日韩欧美中文字幕| 亚洲天堂国产精品一区在线| 欧美黄色片欧美黄色片| 亚洲精品av麻豆狂野| 国产一区二区三区在线臀色熟女| 好看av亚洲va欧美ⅴa在| 51午夜福利影视在线观看| 久久久久国产精品人妻aⅴ院| 一级片免费观看大全| 麻豆成人av在线观看| 给我免费播放毛片高清在线观看| 亚洲成人久久性| 久久久精品国产亚洲av高清涩受| 免费高清在线观看日韩| 中文资源天堂在线| 哪里可以看免费的av片| 中文资源天堂在线| 啪啪无遮挡十八禁网站| 精品一区二区三区视频在线观看免费| 操出白浆在线播放| 亚洲国产精品999在线| 成人亚洲精品av一区二区| 国产精品香港三级国产av潘金莲| 欧美国产日韩亚洲一区| 亚洲精品美女久久久久99蜜臀| 麻豆一二三区av精品| 国产精品爽爽va在线观看网站 | 成人亚洲精品一区在线观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品美女久久久久99蜜臀| 一级毛片高清免费大全| а√天堂www在线а√下载| 免费无遮挡裸体视频| 一本久久中文字幕| 97碰自拍视频| 色综合欧美亚洲国产小说| 在线永久观看黄色视频| 脱女人内裤的视频| 91av网站免费观看| 侵犯人妻中文字幕一二三四区| 成人手机av| 亚洲全国av大片| 亚洲成国产人片在线观看| 久久精品aⅴ一区二区三区四区| 国产在线精品亚洲第一网站| 免费高清在线观看日韩| 曰老女人黄片| 制服诱惑二区| 国产野战对白在线观看| 日本在线视频免费播放| 变态另类丝袜制服| 超碰成人久久| 亚洲五月婷婷丁香| 人人妻,人人澡人人爽秒播| 人妻久久中文字幕网| 宅男免费午夜| 色播在线永久视频| 亚洲一码二码三码区别大吗| 国产午夜福利久久久久久| 国产国语露脸激情在线看| 亚洲国产精品sss在线观看| 在线播放国产精品三级| 一二三四社区在线视频社区8| 亚洲午夜理论影院| 免费搜索国产男女视频| 亚洲国产欧美日韩在线播放| 欧美中文日本在线观看视频| 亚洲成人国产一区在线观看| 色播在线永久视频| 日本 欧美在线| 黄色 视频免费看| 国产视频内射| 国产精品久久久久久精品电影 | 成人午夜高清在线视频 | e午夜精品久久久久久久| 国产高清激情床上av| 怎么达到女性高潮| 一夜夜www| 女同久久另类99精品国产91| 日本一区二区免费在线视频| 18禁黄网站禁片免费观看直播| 成人18禁高潮啪啪吃奶动态图| 亚洲无线在线观看| 曰老女人黄片| 在线观看日韩欧美| 搡老岳熟女国产| 九色国产91popny在线| 国产精品爽爽va在线观看网站 | 国产亚洲精品久久久久5区| 日本成人三级电影网站| 中文字幕精品亚洲无线码一区 | 性欧美人与动物交配| 精品国产美女av久久久久小说| 国产激情偷乱视频一区二区| 亚洲国产精品成人综合色| 欧美亚洲日本最大视频资源| 欧美中文日本在线观看视频| 亚洲男人天堂网一区| 国产熟女xx| 18美女黄网站色大片免费观看| 亚洲专区中文字幕在线| 婷婷丁香在线五月| 国产欧美日韩精品亚洲av| 中文字幕精品免费在线观看视频| 大香蕉久久成人网| 精品久久久久久,| 久久人妻av系列| 欧美性猛交╳xxx乱大交人| 久久国产精品影院| 午夜激情福利司机影院| 午夜福利一区二区在线看| 免费一级毛片在线播放高清视频| 日韩av在线大香蕉| 亚洲全国av大片| www.精华液| 久久香蕉国产精品| 国产乱人伦免费视频| 操出白浆在线播放| 久久精品aⅴ一区二区三区四区| 成人欧美大片| 啪啪无遮挡十八禁网站| 91字幕亚洲| 午夜成年电影在线免费观看| 亚洲无线在线观看| 侵犯人妻中文字幕一二三四区| 非洲黑人性xxxx精品又粗又长| 久热这里只有精品99| 男男h啪啪无遮挡| 国产精品 国内视频| 伊人久久大香线蕉亚洲五| √禁漫天堂资源中文www| 亚洲精品粉嫩美女一区| 精品福利观看| 国产成人一区二区三区免费视频网站| 黄色女人牲交| 精品国产国语对白av| 这个男人来自地球电影免费观看| 免费看a级黄色片| 最新美女视频免费是黄的| 香蕉丝袜av| 国产一级毛片七仙女欲春2 | 男女床上黄色一级片免费看| 国产精品乱码一区二三区的特点| 制服诱惑二区| 日韩精品免费视频一区二区三区| 一级毛片精品| a级毛片在线看网站| 亚洲精品在线观看二区| 久久午夜综合久久蜜桃| 国产亚洲精品久久久久5区| 亚洲一区二区三区色噜噜| 国产亚洲欧美98| 亚洲av美国av| www日本在线高清视频| 亚洲人成电影免费在线| 欧美黑人精品巨大| 黑人操中国人逼视频| 日韩大码丰满熟妇| 亚洲成人精品中文字幕电影| 怎么达到女性高潮| 免费搜索国产男女视频| 欧美又色又爽又黄视频| 久久午夜综合久久蜜桃| 性色av乱码一区二区三区2| 少妇粗大呻吟视频| 国产精品 欧美亚洲| 听说在线观看完整版免费高清| 在线观看www视频免费| 婷婷丁香在线五月| 琪琪午夜伦伦电影理论片6080| 变态另类成人亚洲欧美熟女| 最新在线观看一区二区三区| 长腿黑丝高跟| 国产精品98久久久久久宅男小说| e午夜精品久久久久久久| 亚洲国产欧洲综合997久久, | 一级作爱视频免费观看| 两个人免费观看高清视频| 亚洲国产日韩欧美精品在线观看 | 亚洲自拍偷在线| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产精品成人综合色| 久久欧美精品欧美久久欧美| xxx96com| 欧美在线一区亚洲| 午夜免费成人在线视频| 国内精品久久久久精免费| 男女视频在线观看网站免费 | 欧美日韩乱码在线| 国产精品综合久久久久久久免费| 国产av一区在线观看免费| 日韩欧美免费精品| 久久久久久九九精品二区国产 | 一区二区三区激情视频| 老司机福利观看| 亚洲狠狠婷婷综合久久图片| 亚洲精品中文字幕在线视频| 日本三级黄在线观看| 国内毛片毛片毛片毛片毛片| 一进一出抽搐gif免费好疼| 精品久久久久久久久久久久久 | 99久久99久久久精品蜜桃| 亚洲成av人片免费观看| 久久精品aⅴ一区二区三区四区| 亚洲一区二区三区不卡视频|