• <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无| 人人妻人人澡人人爽人人夜夜| 久久久久国产精品人妻一区二区| 五月玫瑰六月丁香| 大片电影免费在线观看免费| 观看美女的网站| 久久99热这里只频精品6学生| 久久99一区二区三区| 亚洲三级黄色毛片| 色婷婷久久久亚洲欧美| 日韩中文字幕视频在线看片| 一级a做视频免费观看| 在线观看美女被高潮喷水网站| 99久国产av精品国产电影| 最近手机中文字幕大全| 国产老妇伦熟女老妇高清| 亚洲四区av| 国产老妇伦熟女老妇高清| 亚洲国产成人一精品久久久| av黄色大香蕉| 一区二区日韩欧美中文字幕 | 亚洲,欧美,日韩| 欧美日韩综合久久久久久| 夜夜爽夜夜爽视频| 又黄又爽又刺激的免费视频.| 在线观看免费视频网站a站| 国产一区二区三区综合在线观看 | 国产免费视频播放在线视频| 五月天丁香电影| 建设人人有责人人尽责人人享有的| 99热全是精品| 亚洲国产欧美日韩在线播放| 乱码一卡2卡4卡精品| 在线观看三级黄色| 国国产精品蜜臀av免费| 校园人妻丝袜中文字幕| 18禁在线无遮挡免费观看视频| 在现免费观看毛片| 亚洲精品成人av观看孕妇| 在线亚洲精品国产二区图片欧美| 午夜影院在线不卡| 99re6热这里在线精品视频| 免费观看av网站的网址| 两性夫妻黄色片 | 91久久精品国产一区二区三区| 精品少妇内射三级| 我要看黄色一级片免费的| 水蜜桃什么品种好| 如何舔出高潮| 国产亚洲欧美精品永久| 一本—道久久a久久精品蜜桃钙片| 18禁裸乳无遮挡动漫免费视频| 亚洲一级一片aⅴ在线观看| 桃花免费在线播放| 七月丁香在线播放| 男女边吃奶边做爰视频| 美女xxoo啪啪120秒动态图| 国产成人精品久久久久久| 久久99精品国语久久久| 99久久中文字幕三级久久日本| 久久热在线av| 在线 av 中文字幕| 精品99又大又爽又粗少妇毛片| 成年动漫av网址| 国产白丝娇喘喷水9色精品| 超色免费av| 大码成人一级视频| 亚洲欧美日韩卡通动漫| 国产精品一二三区在线看| 在线观看人妻少妇| 91精品三级在线观看| 久久国产精品大桥未久av| 国产精品.久久久| 午夜福利在线观看免费完整高清在| 国产成人精品在线电影| 精品人妻熟女毛片av久久网站| 欧美激情国产日韩精品一区| 少妇被粗大的猛进出69影院 | 乱人伦中国视频| 国产亚洲精品第一综合不卡 | 美女国产视频在线观看| 最近最新中文字幕免费大全7| 一本色道久久久久久精品综合| 97在线视频观看| 夜夜骑夜夜射夜夜干| 纯流量卡能插随身wifi吗| 男人添女人高潮全过程视频| 亚洲性久久影院| 国产成人精品无人区| 只有这里有精品99| 大码成人一级视频| 精品少妇内射三级| 在线免费观看不下载黄p国产| 国产精品蜜桃在线观看| 免费少妇av软件| 亚洲精品美女久久久久99蜜臀 | 亚洲av电影在线观看一区二区三区| 久久久精品区二区三区| 成人免费观看视频高清| av又黄又爽大尺度在线免费看| 中文乱码字字幕精品一区二区三区| 夜夜爽夜夜爽视频| 男人操女人黄网站| 精品卡一卡二卡四卡免费| 新久久久久国产一级毛片| 亚洲第一区二区三区不卡| 中文精品一卡2卡3卡4更新| 男人爽女人下面视频在线观看| h视频一区二区三区| 国产极品天堂在线| 美女内射精品一级片tv| 一个人免费看片子| 日韩av在线免费看完整版不卡| 国产亚洲最大av| 免费看av在线观看网站| 国产av精品麻豆| 亚洲av男天堂| 欧美 亚洲 国产 日韩一| 少妇精品久久久久久久| 18在线观看网站| 视频中文字幕在线观看| 亚洲欧美色中文字幕在线| 精品福利永久在线观看| 黄色怎么调成土黄色| 成年av动漫网址| 熟女电影av网| 久久午夜综合久久蜜桃| 十八禁高潮呻吟视频| 永久网站在线| 91久久精品国产一区二区三区| 亚洲精品中文字幕在线视频| 亚洲第一区二区三区不卡| 夫妻午夜视频| 中文字幕免费在线视频6| 波野结衣二区三区在线| 男女国产视频网站| av在线app专区| 精品国产露脸久久av麻豆| 亚洲av福利一区| a级毛片黄视频| 国产不卡av网站在线观看| av免费观看日本| 午夜老司机福利剧场| 国产免费又黄又爽又色| 2022亚洲国产成人精品| 高清不卡的av网站| 一个人免费看片子| 国产精品成人在线| 国产精品熟女久久久久浪| 赤兔流量卡办理| 91aial.com中文字幕在线观看| 午夜福利视频精品| 赤兔流量卡办理| 欧美日韩精品成人综合77777| 午夜91福利影院| 免费人成在线观看视频色| 亚洲精品,欧美精品| 久久 成人 亚洲| 国产精品麻豆人妻色哟哟久久| 母亲3免费完整高清在线观看 | 国产成人免费无遮挡视频| 国产爽快片一区二区三区| 国产av精品麻豆| 欧美另类一区| 欧美成人精品欧美一级黄| 色网站视频免费| 久久人人爽人人爽人人片va| 激情视频va一区二区三区| 看十八女毛片水多多多| 高清av免费在线| 久久精品夜色国产| 国产高清不卡午夜福利| 男女边吃奶边做爰视频| 三上悠亚av全集在线观看| 精品熟女少妇av免费看| 91在线精品国自产拍蜜月| 国产色爽女视频免费观看| 亚洲激情五月婷婷啪啪| 亚洲欧美清纯卡通| 欧美激情国产日韩精品一区| 日韩av免费高清视频| 超碰97精品在线观看| 成人亚洲欧美一区二区av| 九色成人免费人妻av| 国内精品宾馆在线| 香蕉精品网在线| 国产午夜精品一二区理论片| 色5月婷婷丁香| 永久网站在线| 七月丁香在线播放| av国产精品久久久久影院| 秋霞在线观看毛片| 青春草视频在线免费观看| 一级毛片黄色毛片免费观看视频| 免费播放大片免费观看视频在线观看| 欧美性感艳星| 欧美日韩av久久| 国产成人精品一,二区| 国产一区二区激情短视频 | 国产一区有黄有色的免费视频| 老司机影院毛片| 高清av免费在线| 飞空精品影院首页| 亚洲精品久久久久久婷婷小说| 热re99久久国产66热| 日韩精品免费视频一区二区三区 | 一级黄片播放器| 成人免费观看视频高清| 91久久精品国产一区二区三区| 久久精品国产a三级三级三级| 91午夜精品亚洲一区二区三区| 日日爽夜夜爽网站| 欧美激情 高清一区二区三区| 精品久久久久久电影网| 久热久热在线精品观看| 视频在线观看一区二区三区| 高清欧美精品videossex| a级片在线免费高清观看视频| 极品人妻少妇av视频| 秋霞在线观看毛片| 天天操日日干夜夜撸| 国精品久久久久久国模美| 男人添女人高潮全过程视频| 亚洲av.av天堂| 麻豆乱淫一区二区| 国产免费一级a男人的天堂| av国产精品久久久久影院| 成年动漫av网址| 视频区图区小说| 99热6这里只有精品| 九九在线视频观看精品| 久久97久久精品| 精品人妻在线不人妻| 日韩制服骚丝袜av| 女人被躁到高潮嗷嗷叫费观| 在线观看免费日韩欧美大片| 美女国产高潮福利片在线看| 色网站视频免费| 欧美日韩视频高清一区二区三区二| 在线观看人妻少妇| 美女主播在线视频| 青春草国产在线视频| 免费女性裸体啪啪无遮挡网站| xxxhd国产人妻xxx| 久久久精品免费免费高清| 日韩不卡一区二区三区视频在线| 九色亚洲精品在线播放| 国产av码专区亚洲av| 黑人欧美特级aaaaaa片| 精品人妻在线不人妻| 99热网站在线观看| 自线自在国产av| 视频在线观看一区二区三区| 国产精品久久久久久精品电影小说| 国产白丝娇喘喷水9色精品| 成人影院久久| h视频一区二区三区| 久久久久久伊人网av| 久久久久国产精品人妻一区二区| 99国产精品免费福利视频| 国产不卡av网站在线观看| 99国产精品免费福利视频| 一区二区三区精品91| 亚洲欧美日韩卡通动漫| 狠狠精品人妻久久久久久综合| 久久鲁丝午夜福利片| 99热网站在线观看| 大片电影免费在线观看免费| 亚洲美女视频黄频| 亚洲一级一片aⅴ在线观看| 亚洲欧美日韩卡通动漫| 久久久久久人妻| 久久女婷五月综合色啪小说| 国产一区二区三区av在线| 丰满迷人的少妇在线观看| 欧美丝袜亚洲另类| 99国产精品免费福利视频| 看非洲黑人一级黄片| 国产免费视频播放在线视频| 狠狠精品人妻久久久久久综合| 夫妻性生交免费视频一级片| 精品一区二区免费观看| 美国免费a级毛片| 免费黄频网站在线观看国产| 成年美女黄网站色视频大全免费| 啦啦啦在线观看免费高清www| 51国产日韩欧美| 国产av精品麻豆| 一区二区三区四区激情视频| a级片在线免费高清观看视频| 又粗又硬又长又爽又黄的视频| av线在线观看网站| 国产一区二区在线观看日韩| 欧美3d第一页| 91aial.com中文字幕在线观看| 欧美日韩av久久| 国产熟女欧美一区二区| 久久午夜福利片| 欧美激情极品国产一区二区三区 | 久久久a久久爽久久v久久| 国产国语露脸激情在线看| 国产无遮挡羞羞视频在线观看| 母亲3免费完整高清在线观看 | 久久99热6这里只有精品| 久久99一区二区三区| 人妻人人澡人人爽人人| 国产av码专区亚洲av| 国产成人a∨麻豆精品| 亚洲成av片中文字幕在线观看 | 少妇猛男粗大的猛烈进出视频| 人妻一区二区av| 欧美少妇被猛烈插入视频| 一级爰片在线观看| 国产精品人妻久久久影院| 亚洲少妇的诱惑av| 日本vs欧美在线观看视频| 亚洲国产欧美日韩在线播放| 五月伊人婷婷丁香| 18禁动态无遮挡网站| 欧美精品一区二区免费开放| 欧美3d第一页| 欧美日韩视频精品一区| 国产精品久久久久久久久免| 欧美日韩精品成人综合77777| 久久99热6这里只有精品| 亚洲精品乱码久久久久久按摩| 99久国产av精品国产电影| 人体艺术视频欧美日本| 精品久久久精品久久久| av.在线天堂| 丰满少妇做爰视频| 欧美日韩综合久久久久久| 国产精品久久久久久av不卡| 伦理电影大哥的女人| av国产精品久久久久影院| 在线精品无人区一区二区三| 有码 亚洲区| 久热久热在线精品观看| 热99久久久久精品小说推荐| 精品亚洲成a人片在线观看| 欧美日韩一区二区视频在线观看视频在线| 中文精品一卡2卡3卡4更新| 人妻系列 视频| 黄色配什么色好看| 精品午夜福利在线看| 在线观看www视频免费| 精品视频人人做人人爽| 欧美精品人与动牲交sv欧美| 欧美性感艳星| 乱码一卡2卡4卡精品| 大香蕉久久成人网| 亚洲av福利一区| 下体分泌物呈黄色| 欧美老熟妇乱子伦牲交| 亚洲激情五月婷婷啪啪| 最近中文字幕高清免费大全6| 国产xxxxx性猛交| 美女主播在线视频| 五月天丁香电影| 久久久久国产网址| 少妇的逼好多水| 飞空精品影院首页| 久久免费观看电影| 午夜91福利影院| 免费高清在线观看日韩| 日韩免费高清中文字幕av| 美女脱内裤让男人舔精品视频| 午夜91福利影院| 精品国产一区二区三区久久久樱花| 99九九在线精品视频| 国产精品一区二区在线观看99| 少妇猛男粗大的猛烈进出视频| 亚洲国产日韩一区二区| 国产成人免费无遮挡视频| 女人久久www免费人成看片| 91在线精品国自产拍蜜月| 日本黄色日本黄色录像| av卡一久久| 国产精品.久久久| 在线观看免费高清a一片| 五月开心婷婷网| 啦啦啦中文免费视频观看日本| 亚洲精品美女久久久久99蜜臀 | 伦精品一区二区三区| 黑人欧美特级aaaaaa片| 久久久久精品人妻al黑| 国产成人a∨麻豆精品| 国产精品一区二区在线观看99| 王馨瑶露胸无遮挡在线观看| 九九在线视频观看精品| 久久久久网色| 婷婷色麻豆天堂久久| 99热这里只有是精品在线观看| 国产精品偷伦视频观看了| 免费看av在线观看网站| 欧美日韩精品成人综合77777| 三上悠亚av全集在线观看| 搡女人真爽免费视频火全软件| 国产深夜福利视频在线观看| 久久久久视频综合| 一区二区日韩欧美中文字幕 | 成人综合一区亚洲| 国产精品三级大全| 免费在线观看完整版高清| 久久人人爽av亚洲精品天堂| 欧美激情极品国产一区二区三区 | 一级片免费观看大全| 亚洲国产精品成人久久小说| 人人妻人人澡人人爽人人夜夜| 七月丁香在线播放| 成人黄色视频免费在线看| 最近手机中文字幕大全| 国产伦理片在线播放av一区| 亚洲av成人精品一二三区| 国产精品国产三级国产专区5o| 亚洲欧美精品自产自拍| 国产成人a∨麻豆精品| 熟女av电影| 爱豆传媒免费全集在线观看| 午夜精品国产一区二区电影| 日韩人妻精品一区2区三区| 日韩欧美精品免费久久| 欧美激情国产日韩精品一区| 国产亚洲午夜精品一区二区久久| 日日啪夜夜爽| 五月伊人婷婷丁香| 三级国产精品片| 中文天堂在线官网| 又黄又粗又硬又大视频| 色94色欧美一区二区| xxxhd国产人妻xxx| 亚洲av在线观看美女高潮| 亚洲美女搞黄在线观看| 边亲边吃奶的免费视频| 黄片播放在线免费| tube8黄色片| 国产精品偷伦视频观看了| 美女国产高潮福利片在线看| 久久久久久久久久久免费av| 丰满迷人的少妇在线观看| 日韩精品免费视频一区二区三区 | 丰满乱子伦码专区| 欧美97在线视频| 日韩不卡一区二区三区视频在线| 99国产综合亚洲精品| 啦啦啦视频在线资源免费观看| 老司机影院毛片| 在现免费观看毛片| 婷婷色麻豆天堂久久| 亚洲国产欧美日韩在线播放| 精品久久久久久电影网| 亚洲丝袜综合中文字幕| 黑人欧美特级aaaaaa片| 免费av不卡在线播放| 五月天丁香电影| 18在线观看网站| 国产精品蜜桃在线观看| 黑人高潮一二区| 精品人妻在线不人妻| 国产综合精华液| 国产免费现黄频在线看| 美女国产高潮福利片在线看| 亚洲美女黄色视频免费看| av线在线观看网站| 有码 亚洲区| 深夜精品福利| 18禁在线无遮挡免费观看视频| 寂寞人妻少妇视频99o| 亚洲欧美日韩另类电影网站| 制服诱惑二区| 高清在线视频一区二区三区| 国产激情久久老熟女| 亚洲久久久国产精品| 99九九在线精品视频| 交换朋友夫妻互换小说| 热re99久久精品国产66热6| 久久这里只有精品19| 一级爰片在线观看| 久久人人爽av亚洲精品天堂| 99国产综合亚洲精品| 九色亚洲精品在线播放| 免费在线观看黄色视频的| 一级爰片在线观看| 久久综合国产亚洲精品| 九九在线视频观看精品| 在线免费观看不下载黄p国产| 中文字幕人妻熟女乱码| 夜夜爽夜夜爽视频| 黄色怎么调成土黄色| 熟女电影av网| 麻豆精品久久久久久蜜桃| 亚洲av成人精品一二三区| 男女无遮挡免费网站观看| 黑人欧美特级aaaaaa片| 国产熟女欧美一区二区| 青春草国产在线视频| 日韩制服骚丝袜av| 另类亚洲欧美激情| 国产1区2区3区精品| 高清黄色对白视频在线免费看| 性色av一级| 亚洲av男天堂| 18禁国产床啪视频网站| 亚洲国产最新在线播放| 精品久久久精品久久久| 色婷婷久久久亚洲欧美| 亚洲国产日韩一区二区| 国产成人av激情在线播放| 久久久久久伊人网av| 色94色欧美一区二区| 久久 成人 亚洲| 精品人妻熟女毛片av久久网站| 久久久国产一区二区| 夜夜爽夜夜爽视频| 少妇的逼好多水| 精品一品国产午夜福利视频| 午夜久久久在线观看| 午夜福利乱码中文字幕| 一级毛片电影观看| 18+在线观看网站| 制服人妻中文乱码| 亚洲精品一二三| 国产亚洲最大av| 欧美老熟妇乱子伦牲交| 两个人看的免费小视频| 亚洲久久久国产精品| 啦啦啦中文免费视频观看日本| 黄色一级大片看看| av国产精品久久久久影院| 国产精品.久久久| 亚洲国产最新在线播放| 香蕉丝袜av| 国产老妇伦熟女老妇高清| 母亲3免费完整高清在线观看 | 国产精品三级大全| 亚洲精华国产精华液的使用体验| 亚洲av电影在线观看一区二区三区| 亚洲国产精品999| 一区二区三区乱码不卡18| 超色免费av| 日本91视频免费播放| 搡老乐熟女国产| 美女脱内裤让男人舔精品视频| 青春草国产在线视频| 18禁裸乳无遮挡动漫免费视频| 国产毛片在线视频| 日本欧美国产在线视频| 日本午夜av视频| 成年女人在线观看亚洲视频| 黄色毛片三级朝国网站| 中文天堂在线官网| av线在线观看网站| 成人18禁高潮啪啪吃奶动态图| 插逼视频在线观看| 国产精品久久久久久精品古装| 精品久久蜜臀av无| 岛国毛片在线播放| 国产 精品1| 日韩制服骚丝袜av| 欧美激情 高清一区二区三区| 婷婷色综合www| 制服诱惑二区| 日韩精品免费视频一区二区三区 | 亚洲,欧美,日韩| 亚洲av免费高清在线观看| 婷婷色综合大香蕉| 两个人免费观看高清视频| 毛片一级片免费看久久久久| 2022亚洲国产成人精品| 丝袜喷水一区| 看免费av毛片| 日韩人妻精品一区2区三区| 伊人久久国产一区二区| 成人漫画全彩无遮挡| 永久网站在线| 男人舔女人的私密视频| 久久ye,这里只有精品| 欧美日韩一区二区视频在线观看视频在线| 国产精品久久久久久精品古装| 国产亚洲欧美精品永久| 国产精品一区二区在线观看99| 建设人人有责人人尽责人人享有的| 亚洲精品色激情综合| 精品酒店卫生间| 少妇人妻精品综合一区二区| 一本—道久久a久久精品蜜桃钙片| 国产探花极品一区二区| av女优亚洲男人天堂| 精品久久久久久电影网| 亚洲,欧美,日韩| 下体分泌物呈黄色| 久久久久久伊人网av| 精品人妻熟女毛片av久久网站| 欧美日韩av久久| 美女国产视频在线观看| 亚洲在久久综合| 国产黄色视频一区二区在线观看| 国产成人免费观看mmmm| 在线精品无人区一区二区三| 黄色视频在线播放观看不卡| 免费人妻精品一区二区三区视频| 女的被弄到高潮叫床怎么办| 男的添女的下面高潮视频| 精品久久国产蜜桃| 丝袜喷水一区| 国国产精品蜜臀av免费| 18+在线观看网站| 国精品久久久久久国模美| 99精国产麻豆久久婷婷| 色婷婷久久久亚洲欧美| 久久99精品国语久久久| 国产深夜福利视频在线观看| 一级片'在线观看视频|