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

    含在線處理的熔鹽堆源項(xiàng)計(jì)算

    2014-02-17 01:37:50張志宏夏曉彬朱興望
    核技術(shù) 2014年2期
    關(guān)鍵詞:熔鹽核素堆芯

    張志宏 夏曉彬 朱興望 蔡 軍

    1(中國科學(xué)院上海應(yīng)用物理研究所 嘉定園區(qū) 上海 201800)

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

    含在線處理的熔鹽堆源項(xiàng)計(jì)算

    張志宏1,2夏曉彬1朱興望1蔡 軍1

    1(中國科學(xué)院上海應(yīng)用物理研究所 嘉定園區(qū) 上海 201800)

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

    熔鹽堆是第四代國際核能論壇推薦的6種先進(jìn)四代堆之一,具有良好的中子經(jīng)濟(jì)性、固有安全性、可在線后處理等特點(diǎn)。基于熔鹽燃料和冷卻劑的流動(dòng)性,熔鹽堆可對(duì)堆內(nèi)產(chǎn)生的放射性氣體進(jìn)行在線移除,周期性地從熔鹽泵中除去135Xe和85Kr等中子毒物,以提高燃料循環(huán)過程中的中子經(jīng)濟(jì)性。目前主要的源項(xiàng)計(jì)算軟件并不完全適用于熔鹽堆在線除氣情況下的源項(xiàng)計(jì)算,為了解熔鹽堆在線處理情況下的源項(xiàng)分布,為熔鹽堆輻射防護(hù)設(shè)計(jì)提供參考依據(jù),并驗(yàn)證熔鹽堆在線處理情況下源項(xiàng)計(jì)算方法的可行性和計(jì)算結(jié)果的可靠性,采用實(shí)際運(yùn)行過的美國橡樹嶺國家實(shí)驗(yàn)室熔鹽實(shí)驗(yàn)堆(Molten Salt Reactor Experiment, MSRE)參數(shù),應(yīng)用自主開發(fā)的PostTRITON程序計(jì)算了含在線處理情況下的源項(xiàng),分析了熔鹽堆裂變產(chǎn)物產(chǎn)額及氚的產(chǎn)生情況,并應(yīng)用SCALE5.1軟件包中TRITON模塊計(jì)算了無在線處理情況下的源項(xiàng),分析了數(shù)據(jù)偏差的原因,并與橡樹嶺國家實(shí)驗(yàn)室的計(jì)算數(shù)據(jù)進(jìn)行了比較,計(jì)算結(jié)果在±8%的范圍內(nèi)符合一致。通過本研究,驗(yàn)證了在線處理下熔鹽堆源項(xiàng)計(jì)算方法的可行性和結(jié)果的可靠性,可為熔鹽堆的輻射防護(hù)設(shè)計(jì)提供重要的參考。

    SCALE5.1,熔鹽實(shí)驗(yàn)堆(Molten Salt Reactor Experiment, MSRE),源項(xiàng)計(jì)算,在線處理

    熔鹽堆(Molten Salt Reactor, MSR)是第四代國際核能論壇(GIF)推薦的6種先進(jìn)四代堆候選堆型之一,采用流動(dòng)的熔鹽作為燃料和冷卻劑,具有資源的可持續(xù)性、高度的安全性、良好的經(jīng)濟(jì)性和可靠的防擴(kuò)散性等特點(diǎn),滿足核能可持續(xù)發(fā)展的需要[1]。熔鹽實(shí)驗(yàn)堆 (Molten Salt Reactor Experiment, MSRE)是美國橡樹嶺國家實(shí)驗(yàn)室(Oak Ridge National Laboratory, ORNL)為了驗(yàn)證熔鹽堆的可行性而建立的實(shí)驗(yàn)堆。根據(jù)文獻(xiàn)[2],MSRE在運(yùn)行過程中分別采用了235U和233U作為燃料并進(jìn)行在線除氣,其運(yùn)行歷史為熔鹽堆后續(xù)研究及相關(guān)數(shù)據(jù)庫的建立提供了寶貴的經(jīng)驗(yàn)。為保證乏燃料儲(chǔ)存及處理的安全,Bell M J 曾對(duì)ORIGEN軟件進(jìn)行修改并計(jì)算了MSRE的源項(xiàng),發(fā)布了計(jì)算結(jié)果[2]。

    反應(yīng)堆源項(xiàng)計(jì)算給出的核素濃度、裂變產(chǎn)物活度及衰變熱等結(jié)果,為反應(yīng)堆屏蔽設(shè)計(jì)、輻射監(jiān)測(cè)儀表量程設(shè)計(jì)以及事故后果分析等提供重要的數(shù)據(jù)依據(jù)。熔鹽流動(dòng)性和在線處理影響了熔鹽堆的源項(xiàng)分布,目前主流的源項(xiàng)計(jì)算軟件并不具備熔鹽堆液態(tài)燃料流動(dòng)和在線處理情況下源項(xiàng)計(jì)算的功能,為了解熔鹽堆在線處理情況下的源項(xiàng)分布,為熔鹽堆輻射防護(hù)設(shè)計(jì)提供參考依據(jù),本文應(yīng)用自主開發(fā)編寫的具有在線移除計(jì)算功能的擴(kuò)展程序PostTRITON以實(shí)現(xiàn)含在線處理的熔鹽堆源項(xiàng)計(jì)算,并采用MSRE熔鹽堆參數(shù),計(jì)算了MSRE熔鹽堆在線處理情況下的裂變產(chǎn)物產(chǎn)額及氚的產(chǎn)生情況。同時(shí)應(yīng)用SCALE5.1軟件包中的TRITON (Transport Rigor Implemented with Time-dependent Operation for Neutronic depletion)模塊進(jìn)行了無在線處理情況下的源項(xiàng)計(jì)算,分析了數(shù)據(jù)偏差的原因,并與橡樹嶺國家實(shí)驗(yàn)室的計(jì)算數(shù)據(jù)進(jìn)行了比較。通過本研究,可驗(yàn)證熔鹽堆在線處理情況下源項(xiàng)計(jì)算方法的可行性和可靠性,同時(shí)為熔鹽堆的輻射防護(hù)設(shè)計(jì)提供參考依據(jù)。

    1 計(jì)算方法

    1.1 MSRE堆芯結(jié)構(gòu)

    MSRE堆芯由熔鹽、石墨、哈氏合金組成,直徑為140 cm,高度為163 cm[3],共有513根完整的石墨棒和104塊的石墨塊。石墨棒為5.08 cm× 5.08cm×160 cm的結(jié)構(gòu),并四側(cè)鏤空,兩兩拼接形成3.048 cm×1.016 cm的熔鹽通道[4],總共拼出1140個(gè)熔鹽通道、3個(gè)控制棒通道和1個(gè)樣品通道。石墨棒結(jié)構(gòu)及布局如圖1所示。

    本文選取MSRE運(yùn)行初期的235U燃料展開源項(xiàng)計(jì)算,熔鹽燃料LiF-BeF2-TrF4-UF4(65-29.2-5-0.8 mol%)堆芯入口溫度約635 °C,出口溫度約663 °C,其中鈾的原子濃度為235U 35%、238U 64.4%、234U 0.3%、236U 0.3%[5]。

    圖1 MSRE石墨棒結(jié)構(gòu)Fig.1 Arrangement of MSRE graphite stringer.

    1.2 計(jì)算方法

    根據(jù)MSRE堆芯結(jié)構(gòu)建立圖2所示的計(jì)算模型,模型的建立需要考慮以不同的幾何形狀在石墨棒上拼接出4個(gè)燃料半通道,并以陣列的方法排列513根石墨棒構(gòu)建出1140個(gè)燃料通道,堆芯圍桶(Core can)和反應(yīng)堆容器(Vessel)則采用哈氏合金材料[3]。在計(jì)算模型的邊界條件設(shè)置方面,堆芯圍桶設(shè)置為White邊界,該邊界條件適用于圓柱等非直角幾何[6];反應(yīng)堆容器則應(yīng)設(shè)置為鏡面反射邊界,以保證堆內(nèi)中子不損失。整個(gè)計(jì)算過程中假定熔鹽燃料在堆芯內(nèi)均勻分布,采用238群中子截面庫。計(jì)算過程中所需要的輻照時(shí)間、輻照功率參考MSRE的實(shí)際運(yùn)行時(shí)間功率[2],見表1。

    圖2 TRITON建立的MSRE堆芯模型Fig.2 Model of MSRE core in TRITON.

    表1 MSRE運(yùn)行時(shí)間及對(duì)應(yīng)功率[2]Table1 Approximate power history of MSRE fuel salt[2].

    選用SCALE5.1中的TRITON控制模塊[7]進(jìn)行無在線處理情況下的熔鹽堆源項(xiàng)分析。SCALE5.1程序的全稱是Standardized Computer Analyses for Licensing Evaluation,由美國橡樹嶺國家實(shí)驗(yàn)室開發(fā)并維護(hù),主要用于核反應(yīng)堆物理分析、臨界安全分析、輻射屏蔽計(jì)算等方面,并可以對(duì)核設(shè)施的設(shè)計(jì)進(jìn)行安全評(píng)估。TRITON控制模塊主要提供一種集成自動(dòng)化的多群共振截面處理、中子輸運(yùn)計(jì)算和源項(xiàng)計(jì)算功能。TRITON主要集成BONAMI、NITAWL和CENTRM模塊進(jìn)行核素的共振處理,其中BONAMI利用Bondarenko方法計(jì)算不可分辨共振區(qū)的共振自屏蔽截面,NITAWL通過利用Nordheim積分方法求解燃料-慢化劑兩區(qū)中子慢化方程對(duì)分辨共振能區(qū)截面進(jìn)行處理,CENTRM通過利用離散縱標(biāo)(SN)方法求解一維或均勻問題的連續(xù)能量方程,對(duì)可分辨共振能區(qū)截面進(jìn)行處理。中子輸運(yùn)方面,可采用一維離散縱標(biāo)法XSDRNPM模塊、二維離散縱標(biāo)法NEWT模塊,以及三維蒙特卡洛方法KENO V.a和KENO-VI模塊[8-9];源項(xiàng)計(jì)算方面,主要采用ORIGEN-S程序[10]進(jìn)行源項(xiàng)計(jì)算。

    圖3是典型的TRITON計(jì)算流程圖[11]。由圖3,TRITON程序的計(jì)算流程為:(1) 運(yùn)行TRITON/KENO輸入文件,程序BONAMI、NITAWL和CENTRM模塊對(duì)核素進(jìn)行共振截面處理;(2) 共振處理后的數(shù)據(jù)傳遞給三維輸運(yùn)計(jì)算KENO進(jìn)行輸運(yùn)計(jì)算;(3) KMART程序?qū)ψ⒘柯蔬M(jìn)行再處理;(4) COUPLE程序?qū)ENO輸運(yùn)計(jì)算的多群截面數(shù)據(jù)以及注量率數(shù)據(jù)傳遞給ORIGEN-S進(jìn)行計(jì)算;(5) ORIGEN-S計(jì)算規(guī)定時(shí)間步長(zhǎng)內(nèi)生成或消失的核素量;(6) 如未完成所有的時(shí)間步長(zhǎng),計(jì)算的材料組份再傳遞給共振處理模塊,依次循環(huán)計(jì)算;如果完成所有時(shí)間步長(zhǎng),OPUS功能模塊將源項(xiàng)計(jì)算的結(jié)果輸出[12]。

    圖3 TRITON源項(xiàng)計(jì)算流程圖[11]Fig.3 Flowchart of TRITON[11].

    1.3 在線除氣下的源項(xiàng)計(jì)算

    SCALE5.1中TRITON模塊在耦合ORIGEN-S的時(shí)候只簡(jiǎn)單地調(diào)用了其源項(xiàng)計(jì)算功能,并未調(diào)用相應(yīng)的在線移除計(jì)算功能[7],因此不具備在線除氣下的源項(xiàng)計(jì)算能力。在線除氣對(duì)熔鹽堆源項(xiàng)的影響分析主要應(yīng)用自主開發(fā)的具有在線處理情況下源項(xiàng)計(jì)算能力的PostTRITON程序開展。

    在TRITON模塊中輸運(yùn)計(jì)算和截面處理過程通常占了較大的時(shí)間比例,一次輸運(yùn)計(jì)算平均耗時(shí)7min多。根據(jù)ORNL的技術(shù)報(bào)告[2],MSRE運(yùn)行期間每487 s移除38%的惰性氣體Kr和Xe。按照TRITON的計(jì)算流程,若單純只進(jìn)行在線除氣計(jì)算功能的添加,則每進(jìn)行487 s的源項(xiàng)計(jì)算就需要進(jìn)行一次輸運(yùn)計(jì)算和截面處理??紤]到移除周期與反應(yīng)堆運(yùn)行時(shí)間的差別,為了節(jié)省計(jì)算成本,并不需要在每個(gè)移除周期都進(jìn)行輸運(yùn)計(jì)算,可以選擇性地完成這一過程。

    根據(jù)這一思路開發(fā)建立PostTRITON程序以實(shí)現(xiàn)含在線處理情況下的熔鹽堆源項(xiàng)計(jì)算。PostTRITON程序是以TRITON為基礎(chǔ)進(jìn)行擴(kuò)展的計(jì)算程序,擴(kuò)展的內(nèi)容主要是重新組織計(jì)算序列,在TRITON進(jìn)行截面處理后調(diào)用ORIGRN-S源項(xiàng)計(jì)算和在線移除功能[10]。計(jì)算流程與TRITON大致相同,在圖3所示的步驟(5)中添加了ORIGEN-S的在線除氣功能進(jìn)行在線除氣下的源項(xiàng)計(jì)算,并在步驟(5)完成的時(shí)候?qū)Ψ磻?yīng)堆功率進(jìn)行判斷:若反應(yīng)堆運(yùn)行功率發(fā)生改變,則先進(jìn)行步驟(1)-(4);若反應(yīng)堆穩(wěn)功率運(yùn)行,則程序按照用戶設(shè)定的頻率在一段特定時(shí)間后進(jìn)行步驟(1)-(4),其余時(shí)間段只進(jìn)行步驟(5)。例如反應(yīng)堆功率為1.76 MW運(yùn)行185 d,若設(shè)定該時(shí)間段內(nèi)更新截面庫的次數(shù)為5次,則在程序運(yùn)行中的第0 d、37 d、74 d、111 d、148 d進(jìn)行截面更新,其余時(shí)間則只進(jìn)行源項(xiàng)計(jì)算。

    2 計(jì)算結(jié)果與分析

    2.1 核素濃度

    ORNL依據(jù)表1中的MSRE運(yùn)行歷史數(shù)據(jù)計(jì)算了MSRE運(yùn)行714 d后停堆時(shí)的核素產(chǎn)額,并發(fā)布了相關(guān)的計(jì)算結(jié)果[2]。為分析在線處理對(duì)于源項(xiàng)的影響并驗(yàn)證計(jì)算結(jié)果的可靠性,本工作開展了兩類計(jì)算:(I) 應(yīng)用SCALE5.1軟件包中TRITON模塊開展的無在線除氣情況下的源項(xiàng);(II) 應(yīng)用PostTRITON進(jìn)行的在線移除Kr、Xe情況下的源項(xiàng)。計(jì)算結(jié)果見圖4。

    圖4 MSRE停堆后燃料中的源項(xiàng)情況(a) 錒系核素,(b) 裂變產(chǎn)物Fig.4 Inventory of MSRE fuel salt after shutdown. (a) Actinides, (b) Fission products

    從圖4(a)中可以看出,無論是計(jì)算I的結(jié)果還是計(jì)算II的結(jié)果,主要的錒系核素產(chǎn)額均與ORNL的計(jì)算數(shù)據(jù)[2]符合得較好。圖4(b)顯示,計(jì)算II中主要的裂變產(chǎn)物產(chǎn)額與ORNL的計(jì)算數(shù)據(jù)有較好的符合度;而計(jì)算I中部分核素如87Rb、88Sr、89Y、131Xe、135Cs、137Cs和137Ba等,則與ORNL計(jì)算數(shù)據(jù)存在一定的差異,最大差異在兩個(gè)數(shù)量級(jí)左右。導(dǎo)致這一差異的原因,主要是這些核素都是Kr和Xe的衰變相關(guān)的核素,這些核素與Kr和Xe存在如下關(guān)系:

    在線移除惰性氣體改變了這些衰變鏈的平衡關(guān)系,這部分核素在兩種計(jì)算中產(chǎn)額的變化即是在線移除Kr和Xe對(duì)于反應(yīng)堆內(nèi)源項(xiàng)影響的體現(xiàn)。

    2.2 裂變產(chǎn)物比活度及衰變熱

    確定反應(yīng)堆運(yùn)行期間的裂變產(chǎn)物化學(xué)組分及其放射性可為事故情況下放射性物質(zhì)釋放和分布的分析提供參考依據(jù)。反應(yīng)堆中的裂變產(chǎn)物包括近40種不同元素的200余種不同核素,其中多數(shù)帶有放射性,并通過發(fā)射α粒子或β、γ射線而衰變,衰變子體往往也是放射性的。裂變產(chǎn)物所具有的放射性以及衰變過程中產(chǎn)生的熱量是反應(yīng)堆輻射防護(hù)設(shè)計(jì)的重要參考,也是燃料離堆后處理過程中的重要參數(shù)。核能行業(yè)標(biāo)準(zhǔn)指出,乏燃料離堆后需經(jīng)五年以上的冷卻方能進(jìn)行貯存[13]。應(yīng)用TRITON和PostTRITON分別計(jì)算得到了無在線處理和有在線處理兩種情況下停堆冷卻295 d后(總第1009 d)裂變產(chǎn)物的比活度和衰變熱情況,如表2所示。

    可見在線除氣對(duì)于乏燃料中裂變產(chǎn)物的比活度基本不會(huì)產(chǎn)生影響,但是在個(gè)別核素的衰變熱方面則會(huì)存在量級(jí)上的差別,計(jì)算結(jié)果可為后續(xù)乏燃料儲(chǔ)存和處理過程的輻射安全設(shè)計(jì)提供參考。

    2.3 氚的產(chǎn)量

    MSRE使用FLiBe作為燃料載體,其中Li和F等主要成分通過中子反應(yīng)產(chǎn)生相當(dāng)數(shù)量的氚,并以具有強(qiáng)腐蝕性、強(qiáng)擴(kuò)散能力的氣態(tài)氚HT和氟化氚TF等形式存在。熔鹽堆運(yùn)行過程中氚的產(chǎn)生途徑主要為[14]:(1) 快中子與7Li 反應(yīng)產(chǎn)生,閾能為2.82MeV;(2) 熱中子與6Li反應(yīng)產(chǎn)生,反應(yīng)截面比快中子與7Li的反應(yīng)截面大四個(gè)數(shù)量級(jí)[15];(3) 通過三分核裂變反應(yīng)產(chǎn)生。

    氚具有較長(zhǎng)的半衰期(12.43 a),隨著高溫熔鹽的流動(dòng)HT和TF會(huì)進(jìn)入熔鹽回路,導(dǎo)致管道的腐蝕和損傷,影響材料的力學(xué)性能和使用壽命,并經(jīng)管道設(shè)備滲透而向環(huán)境釋放。氚釋放到自然環(huán)境后,會(huì)通過吸入和攝入等途徑進(jìn)入人體。因此,計(jì)算堆內(nèi)氚的產(chǎn)量可為氚的防護(hù)設(shè)計(jì)提供必要的依據(jù)。

    應(yīng)用TRITON建立的圖2模型,分別計(jì)算了經(jīng)同位素分離后7Li富集度為99.90%依次遞增到99.99%、99.993%依次遞增到99.997%情況下氚的產(chǎn)量情況,計(jì)算結(jié)果如圖5所示。由圖5可知,經(jīng)同位素分離到7Li富集度99.90%之后氚產(chǎn)量依然高于1 TBq·(day·MW)-1;分離到7Li富集度99.96%以上時(shí)氚的產(chǎn)量可控制在1 TBq·(day·MW)-1以下;7Li富集度為99.99%時(shí)氚的產(chǎn)量已可控制在0.1TBq·(day·MW)-1左右。ORNL在源項(xiàng)計(jì)算時(shí)采用的LiF中Li的同位素比為99.997%的7Li、0.003%的6Li[3],由此計(jì)算而得氚的產(chǎn)量為0.0616TBq·(day·MW)-1[16]。本計(jì)算中,當(dāng)7Li富集度為99.997%時(shí),計(jì)算所得的氚產(chǎn)量為0.0619TBq·(day·MW)-1,與ORNL計(jì)算結(jié)果在誤差±0.4%內(nèi)符合。

    圖5 氚產(chǎn)率隨鋰富集度的變化Fig.5 Dependence of tritium production on the enrichment of lithium.

    3 結(jié)語

    本文應(yīng)用自主開發(fā)的在線處理計(jì)算擴(kuò)展程序PostTRITON綜合控制TRITON的截面處理功能和ORIGRN-S在線處理下的源項(xiàng)計(jì)算功能,實(shí)現(xiàn)了含在線除氣過程的熔鹽堆源項(xiàng)計(jì)算;同時(shí)應(yīng)用Scale5.1軟件包中的TRITON進(jìn)行了熔鹽堆無在線處理情況下的源項(xiàng)計(jì)算,并與橡樹嶺國家實(shí)驗(yàn)室的計(jì)算數(shù)據(jù)進(jìn)行了比較。計(jì)算結(jié)果表明:

    (1) 無在線處理的計(jì)算中Kr和Xe及其衰變子體的產(chǎn)額數(shù)據(jù)與ORNL計(jì)算數(shù)據(jù)[2]相比差異較大,最大差異在兩個(gè)數(shù)量級(jí)左右;而含在線處理的計(jì)算中這些核素的數(shù)據(jù)與ORNL的計(jì)算數(shù)據(jù)[2]相比在±8%誤差范圍內(nèi)相符,這些結(jié)果表明本研究提供的在線除氣下熔鹽堆源項(xiàng)計(jì)算方法是可行的。

    (2) 明確了熔鹽堆中氚的產(chǎn)量與7Li富集度的定量關(guān)系,可為燃料制備過程中鋰的純化和熔鹽堆輻射安全設(shè)計(jì)提供重要的數(shù)據(jù)參考。

    1 GIF-002-00. A technology roadmap for Generation IV nuclear energy systems[R]. USA: DOE, 2002

    2 Bell M J. Calculated radioactivity of MSRE fuel salt[R]. ORNL-TM-2970, USA: ORNL, 1970

    3 Robertson R C. MSRE design and operations report Part I: Description of reactor design[R]. ORNL-TM-0728, USA: ORNL, 1965

    4 McCoy H E, McNabb M E. Post-irradiation examination of materials from the MSRE[R]. ORNL-TM-4174, USA: ORNL, 1977

    5 Haubenreich P N, Engle J R, Prince B E. MSRE design and operations report Part III: Nuclear Analysis[R]. ORNL-TM-0730, USA: ORNL, 1964

    6 O’Dell R D, Alcouffe R E. Transport calculations for nuclear analyses: theory and guidelines for effective use of transport codes[M]. LA-10983, USA: Los Alamos National Laboratory, 1987

    7 DeHart M D. TRITON: A two-dimensional transport and depletion module for characterization of spent nuclear fuel[M]. ORNL/TM-2005/39, Revision 5.1, Vol. I, Book 3, Sect. T1, USA: ORNL, 2006

    8 Petrie L M, Landers N F, Hollenbach D F, et al. KENO V. a: AN improved Monte Carlo criticality program[M]. ORNL/TM-2005/39, Version 5.1, Vol. II, Book 2, Sect. F11, USA: ORNL, 2006

    9 Hollenbach D F, Petrie L M, Landers N F. KENO-VI: A general quadratic version of the KENO program[M]. ORNL/TM-2005/39, Version 5.1, Vol. II, Book 3, Sect. F17, USA: ORNL, 2006

    10 Gauld I C, Hermann O W, Westfall R M. ORIGEN-S: SCALE system module to calculate fuel depletion, actinide transmutation, fission product buildup and decay, and associated radiation[M]. ORNL/TM-2005/39, Version 5.1, Vol. II, Book 1, Sect. F7, USA: ORNL, 2006

    11 DeHart M D, Bowman S M. Reactor physics methods and analysis capabilities in SCALE[J]. Nuclear Technology, 2011, 174(2): 196-213

    12 熊文綱. 熱堆中釷鈾燃料循環(huán)耗分析[D]. 上海: 中國科學(xué)院上海應(yīng)用物理研究所, 2012

    XIONG Wengang. Thorium-uranium fuel cycle burnup analysis for thermal reactor[D]. Shanghai: Shanghai Institute of Applied Physics, Chinese Academy of Sciences, 2012

    13 EJ/T 878-2011. 乏燃料離堆貯存水池安全設(shè)計(jì)準(zhǔn)則[S].中國核工業(yè)總公司, 2011

    EJ/T 878-2011. Guidelines for the safety design of the storage of spent fuel pit[S]. China National Nuclear Corporation, 2011

    14 黃豫, 劉衛(wèi), 肖德濤, 等. 熔鹽堆中氚的控制和監(jiān)測(cè)[J].核技術(shù), 2011, 34(8): 632- 636

    HUANG Yu, LIU Wei, XIAO Detao, et al. Control and monitoring of tritium in molten salt reactor[J]. Nuclear Techniques, 2011, 34(8): 632-636

    15 蔣國強(qiáng), 羅德禮, 陸光達(dá), 等. 氚與氚的工程技術(shù)[M].北京: 國防工業(yè)出版社, 2007: 130-133

    JIANG Guoqiang, LUO Deli, LU Guangda, et al. Tritium and industry techniques of tritium[M]. Beijing: National Defense Industry Press, 2007: 130-133

    16 Mays G T, Smith A N, Engle J R. Distribution and behavior of tritium in the coolant-salt technology facility[R]. ORNL/TM-5759, USA: ORNL, 1977

    CLCTL364

    Source terms calculation for the MSRE with on-line removing radioactive gases

    ZHANG Zhihong1,2XIA Xiaobin1ZHU Xingwang1CAI Jun1
    1(Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Jiading Campus, Shanghai 201800, China)
    2(University of Chinese Academy of Sciences, Beijing 100049, China)

    Background: Molten salt reactor (MSR), a candidate of the Generation IV reactors, has the advantages of inherent safety, on-line refueling and good neutron economy. Neutrons absorption of the poisons such as135Xe and85Kr would affect the reactivity of the reactor. As a result, on-line processing is required to remove gases when the reactor is run. Most of the computation codes that are available for calculation of the source terms have not taken liquid fuel flow and online reprocessing into account. Purpose: This paper attempts to analyze the impacts of on-line removing radioactive gases on the source terms and provide reliable source terms for the radiation protection of MSR. Methods: Source terms of the Molten Salt Reactor Experiment (MSRE), - originally developed at the Oak Ridge National Laboratory (ORNL) of USA, with/without on-line processing are calculated by using self-developed code PostTRITON and the TRITON module of the SCALE5.1 package. Results: The productions of the radioactivity of fission products and tritium were calculated for MSRE, and compared with the data given by ORNL. The results were consistent within the error of ±8%. Conclusion: The calculation method of the source terms for the MSRE with on-line removing radioactive gases is feasible and reliable to provide important data for the radiation protection of MSRs.

    SCALE5.1, Molten Salt Reactor Experiment (MSRE), Source terms, Online processing

    TL364

    10.11889/j.0253-3219.2014.hjs.37.020603

    中國科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)(XDA02005002、XDA02005004)資助

    張志宏,男,1986年出生,2009年畢業(yè)于山東大學(xué),現(xiàn)為中國科學(xué)院上海應(yīng)用物理研究所博士研究生,反應(yīng)堆安全

    夏曉彬,E-mail: xiaxiaobin@sinap.ac.cn

    2013-10-18,

    2013-12-19

    猜你喜歡
    熔鹽核素堆芯
    核素分類開始部分的6種7核素小片分布
    熔鹽在片堿生產(chǎn)中的應(yīng)用
    核素分類的4量子數(shù)
    NaF-KF熔鹽體系制備Ti2CTx材料的研究
    應(yīng)用CDAG方法進(jìn)行EPR機(jī)組的嚴(yán)重事故堆芯損傷研究
    純鈦的熔鹽滲硼
    基于Hoogenboom基準(zhǔn)模型的SuperMC全堆芯計(jì)算能力校驗(yàn)
    核技術(shù)(2016年4期)2016-08-22 09:05:32
    大型燃?xì)馊埯}爐的研發(fā)和工藝控制
    工業(yè)爐(2016年1期)2016-02-27 12:34:11
    壓水堆堆芯中應(yīng)用可燃毒物的兩個(gè)重要實(shí)驗(yàn)
    基于SOP規(guī)程的大亞灣堆芯冷卻監(jiān)測(cè)系統(tǒng)改造
    久久精品亚洲av国产电影网| 国产熟女午夜一区二区三区| 国产成人啪精品午夜网站| 亚洲精品一二三| 老司机深夜福利视频在线观看| 亚洲人成电影免费在线| 男女午夜视频在线观看| 成人18禁在线播放| 交换朋友夫妻互换小说| 精品国产一区二区三区四区第35| 99久久99久久久精品蜜桃| 99热网站在线观看| 国产精品av久久久久免费| 成人av一区二区三区在线看| 搡老乐熟女国产| 交换朋友夫妻互换小说| 久久天堂一区二区三区四区| 黄网站色视频无遮挡免费观看| 搡老乐熟女国产| 欧美日韩福利视频一区二区| 老司机影院毛片| 国产日韩欧美在线精品| 国产精品熟女久久久久浪| 亚洲国产精品一区二区三区在线| 伊人久久大香线蕉亚洲五| 99国产精品一区二区三区| 熟女少妇亚洲综合色aaa.| 欧美日韩福利视频一区二区| 久久人妻福利社区极品人妻图片| 国产又色又爽无遮挡免费看| 韩国精品一区二区三区| 亚洲欧美日韩另类电影网站| 亚洲国产欧美在线一区| 大香蕉久久成人网| 女人久久www免费人成看片| 99国产极品粉嫩在线观看| 一本综合久久免费| 搡老乐熟女国产| 国产97色在线日韩免费| 啪啪无遮挡十八禁网站| 国产一区二区激情短视频| 三上悠亚av全集在线观看| 在线十欧美十亚洲十日本专区| 欧美日韩亚洲国产一区二区在线观看 | 久久国产精品影院| 国产人伦9x9x在线观看| 亚洲人成电影免费在线| 人成视频在线观看免费观看| 宅男免费午夜| 亚洲精品中文字幕在线视频| 美女福利国产在线| 国产单亲对白刺激| 久久久久国内视频| 在线天堂中文资源库| 日韩免费av在线播放| 我要看黄色一级片免费的| 成人国语在线视频| 99riav亚洲国产免费| 久久中文字幕人妻熟女| 亚洲国产毛片av蜜桃av| 在线观看66精品国产| 咕卡用的链子| 脱女人内裤的视频| √禁漫天堂资源中文www| 91麻豆av在线| 欧美精品一区二区免费开放| 国产亚洲精品久久久久5区| 久久精品国产亚洲av香蕉五月 | 亚洲va日本ⅴa欧美va伊人久久| 精品国产国语对白av| 在线观看免费高清a一片| 午夜91福利影院| 亚洲欧美精品综合一区二区三区| 丰满饥渴人妻一区二区三| 激情在线观看视频在线高清 | 人人妻人人澡人人爽人人夜夜| 国产在线免费精品| 亚洲一码二码三码区别大吗| 久久精品国产亚洲av高清一级| 国产免费福利视频在线观看| 啦啦啦 在线观看视频| 在线观看一区二区三区激情| 欧美亚洲日本最大视频资源| a级毛片黄视频| 亚洲午夜理论影院| 免费av中文字幕在线| 亚洲伊人久久精品综合| 国产极品粉嫩免费观看在线| 亚洲av欧美aⅴ国产| 亚洲精品中文字幕一二三四区 | 正在播放国产对白刺激| 可以免费在线观看a视频的电影网站| 久9热在线精品视频| 国产成人精品久久二区二区91| 国产一区二区三区在线臀色熟女 | 免费不卡黄色视频| 免费观看a级毛片全部| 汤姆久久久久久久影院中文字幕| 久久中文字幕一级| 少妇裸体淫交视频免费看高清 | 久久久久久人人人人人| 欧美精品啪啪一区二区三区| 麻豆乱淫一区二区| 搡老岳熟女国产| h视频一区二区三区| 深夜精品福利| 十八禁网站免费在线| 国产一区二区三区综合在线观看| 国产精品免费一区二区三区在线 | 50天的宝宝边吃奶边哭怎么回事| 老司机福利观看| 777米奇影视久久| 亚洲人成伊人成综合网2020| 精品一区二区三区四区五区乱码| 99热国产这里只有精品6| 一夜夜www| 少妇裸体淫交视频免费看高清 | 色播在线永久视频| 国产亚洲欧美精品永久| 99热国产这里只有精品6| 嫁个100分男人电影在线观看| 日韩欧美免费精品| 18禁国产床啪视频网站| 久久精品国产亚洲av香蕉五月 | 亚洲国产欧美网| a在线观看视频网站| 日本精品一区二区三区蜜桃| 精品亚洲乱码少妇综合久久| 亚洲成人免费电影在线观看| 国产成+人综合+亚洲专区| 中文字幕最新亚洲高清| av又黄又爽大尺度在线免费看| 日本av免费视频播放| 久久精品国产亚洲av高清一级| 1024视频免费在线观看| 亚洲精品美女久久久久99蜜臀| 亚洲第一av免费看| 国产精品免费视频内射| 老司机深夜福利视频在线观看| 精品卡一卡二卡四卡免费| 久久ye,这里只有精品| 夜夜夜夜夜久久久久| 男女边摸边吃奶| 叶爱在线成人免费视频播放| 日韩三级视频一区二区三区| 亚洲美女黄片视频| 三上悠亚av全集在线观看| 曰老女人黄片| 水蜜桃什么品种好| 色播在线永久视频| 精品卡一卡二卡四卡免费| 人人妻人人添人人爽欧美一区卜| 免费久久久久久久精品成人欧美视频| 99久久精品国产亚洲精品| 黄色片一级片一级黄色片| 1024视频免费在线观看| 性色av乱码一区二区三区2| 天天躁夜夜躁狠狠躁躁| 欧美亚洲日本最大视频资源| 一边摸一边做爽爽视频免费| 亚洲三区欧美一区| 国产精品亚洲一级av第二区| 老司机在亚洲福利影院| 侵犯人妻中文字幕一二三四区| 亚洲va日本ⅴa欧美va伊人久久| 99精品在免费线老司机午夜| 多毛熟女@视频| 熟女少妇亚洲综合色aaa.| 91麻豆av在线| 欧美日韩亚洲国产一区二区在线观看 | 啦啦啦在线免费观看视频4| 国产成+人综合+亚洲专区| 精品亚洲乱码少妇综合久久| 日本精品一区二区三区蜜桃| 久久天躁狠狠躁夜夜2o2o| 欧美黄色淫秽网站| 啦啦啦中文免费视频观看日本| 99国产精品一区二区三区| 婷婷成人精品国产| 国产日韩欧美在线精品| 亚洲熟妇熟女久久| 大陆偷拍与自拍| 无人区码免费观看不卡 | 亚洲欧美色中文字幕在线| 国产精品一区二区在线观看99| 亚洲熟女精品中文字幕| 50天的宝宝边吃奶边哭怎么回事| 水蜜桃什么品种好| 精品一区二区三卡| 黄色视频,在线免费观看| 国产熟女午夜一区二区三区| 久久99热这里只频精品6学生| 日本黄色视频三级网站网址 | 宅男免费午夜| 亚洲av第一区精品v没综合| 十分钟在线观看高清视频www| 中文字幕高清在线视频| 国产高清视频在线播放一区| 午夜福利乱码中文字幕| 一夜夜www| 成人国产av品久久久| 亚洲第一av免费看| 日韩欧美免费精品| 极品人妻少妇av视频| 日日摸夜夜添夜夜添小说| 中文亚洲av片在线观看爽 | 国产91精品成人一区二区三区 | 国产1区2区3区精品| 精品国产亚洲在线| 亚洲,欧美精品.| 亚洲精品自拍成人| 欧美黄色淫秽网站| 精品久久久精品久久久| 亚洲国产欧美日韩在线播放| 国产一卡二卡三卡精品| 麻豆av在线久日| 国产高清videossex| 下体分泌物呈黄色| 亚洲第一av免费看| 亚洲七黄色美女视频| 久久久久国产一级毛片高清牌| 两性午夜刺激爽爽歪歪视频在线观看 | 人人妻人人爽人人添夜夜欢视频| 亚洲国产中文字幕在线视频| 少妇精品久久久久久久| 日本一区二区免费在线视频| 另类精品久久| 一夜夜www| 国产精品电影一区二区三区 | 男女边摸边吃奶| 精品国产国语对白av| 宅男免费午夜| 欧美一级毛片孕妇| 三上悠亚av全集在线观看| 久久久久久久久久久久大奶| 欧美 亚洲 国产 日韩一| 啦啦啦 在线观看视频| 午夜成年电影在线免费观看| 亚洲美女黄片视频| 美女主播在线视频| 侵犯人妻中文字幕一二三四区| 99国产综合亚洲精品| 久久国产精品人妻蜜桃| 日韩制服丝袜自拍偷拍| 欧美日韩亚洲综合一区二区三区_| 国产精品美女特级片免费视频播放器 | 亚洲五月婷婷丁香| 搡老乐熟女国产| 欧美亚洲日本最大视频资源| 99re6热这里在线精品视频| 国产av又大| 国产成人精品久久二区二区免费| 欧美 日韩 精品 国产| 国产精品国产av在线观看| 亚洲人成电影观看| 91老司机精品| 成人国产一区最新在线观看| 日韩视频一区二区在线观看| 国产精品亚洲av一区麻豆| 色综合婷婷激情| 三上悠亚av全集在线观看| kizo精华| 在线观看免费日韩欧美大片| 久久午夜综合久久蜜桃| 妹子高潮喷水视频| 成年人免费黄色播放视频| 女警被强在线播放| 国产亚洲精品一区二区www | www.熟女人妻精品国产| 成人18禁高潮啪啪吃奶动态图| 黄网站色视频无遮挡免费观看| 亚洲欧美色中文字幕在线| 国产精品av久久久久免费| 三级毛片av免费| 欧美国产精品va在线观看不卡| 巨乳人妻的诱惑在线观看| 狠狠婷婷综合久久久久久88av| 亚洲伊人久久精品综合| 亚洲色图综合在线观看| 国产男靠女视频免费网站| 久久精品91无色码中文字幕| 窝窝影院91人妻| www.自偷自拍.com| 99精国产麻豆久久婷婷| 熟女少妇亚洲综合色aaa.| 一级毛片精品| 国产在线观看jvid| 1024视频免费在线观看| 啦啦啦中文免费视频观看日本| 亚洲精品在线观看二区| av超薄肉色丝袜交足视频| 久久ye,这里只有精品| 不卡一级毛片| 亚洲人成77777在线视频| 国产精品亚洲av一区麻豆| 精品少妇内射三级| 激情视频va一区二区三区| 老汉色av国产亚洲站长工具| 亚洲熟女精品中文字幕| 久久久久久久久免费视频了| 水蜜桃什么品种好| 国产成人系列免费观看| 国产精品.久久久| 久久久精品免费免费高清| 欧美黄色片欧美黄色片| 电影成人av| 国产激情久久老熟女| 又黄又粗又硬又大视频| av线在线观看网站| 国产午夜精品久久久久久| 大香蕉久久成人网| 9色porny在线观看| 亚洲精品中文字幕在线视频| 精品国产一区二区三区四区第35| 免费观看人在逋| av片东京热男人的天堂| 黄色毛片三级朝国网站| 国产一区二区三区视频了| 一区二区日韩欧美中文字幕| 高清黄色对白视频在线免费看| 国产高清国产精品国产三级| 一边摸一边做爽爽视频免费| 日韩欧美一区二区三区在线观看 | 午夜视频精品福利| 亚洲av美国av| 色精品久久人妻99蜜桃| av网站在线播放免费| 高清欧美精品videossex| 丁香欧美五月| 久久久国产精品麻豆| 99香蕉大伊视频| av天堂久久9| 变态另类成人亚洲欧美熟女 | 99re在线观看精品视频| 国产男女超爽视频在线观看| av福利片在线| 窝窝影院91人妻| 999久久久精品免费观看国产| 久久久精品免费免费高清| 最近最新中文字幕大全免费视频| 18禁裸乳无遮挡动漫免费视频| 王馨瑶露胸无遮挡在线观看| 日韩免费av在线播放| 日本黄色日本黄色录像| 在线观看免费日韩欧美大片| 欧美精品av麻豆av| 亚洲五月色婷婷综合| av欧美777| 色精品久久人妻99蜜桃| 丰满人妻熟妇乱又伦精品不卡| 中国美女看黄片| 日韩一卡2卡3卡4卡2021年| 少妇的丰满在线观看| 亚洲av成人一区二区三| 久久精品亚洲熟妇少妇任你| 在线观看一区二区三区激情| 18禁美女被吸乳视频| 一本色道久久久久久精品综合| 国产在视频线精品| 丁香六月天网| 精品国产乱码久久久久久男人| 视频区图区小说| 国产精品二区激情视频| 女人被躁到高潮嗷嗷叫费观| 亚洲av国产av综合av卡| 国产精品国产高清国产av | a级毛片在线看网站| 亚洲一区中文字幕在线| 电影成人av| 一边摸一边抽搐一进一出视频| 亚洲精品国产精品久久久不卡| 免费在线观看日本一区| 国产精品九九99| 看免费av毛片| 成年人免费黄色播放视频| 日本av手机在线免费观看| 天堂动漫精品| 色视频在线一区二区三区| 国产精品二区激情视频| 国产精品电影一区二区三区 | 成年人黄色毛片网站| 国产日韩欧美视频二区| 中文欧美无线码| 热re99久久精品国产66热6| 亚洲五月婷婷丁香| 老司机在亚洲福利影院| 操出白浆在线播放| av线在线观看网站| 国产精品亚洲av一区麻豆| 精品亚洲成a人片在线观看| 国产精品美女特级片免费视频播放器 | 国产三级黄色录像| 亚洲av成人一区二区三| 国产亚洲精品久久久久5区| 国产精品久久久久成人av| 欧美激情高清一区二区三区| 在线 av 中文字幕| av网站免费在线观看视频| 国产精品99久久99久久久不卡| 久久亚洲真实| 97在线人人人人妻| 日韩有码中文字幕| 亚洲少妇的诱惑av| 可以免费在线观看a视频的电影网站| 亚洲国产毛片av蜜桃av| 美女扒开内裤让男人捅视频| 欧美午夜高清在线| 午夜免费鲁丝| 成人免费观看视频高清| 在线观看免费高清a一片| 国产精品一区二区在线观看99| 欧美日韩国产mv在线观看视频| 亚洲人成77777在线视频| 欧美日韩国产mv在线观看视频| 国产精品av久久久久免费| 欧美人与性动交α欧美软件| 亚洲黑人精品在线| 成人特级黄色片久久久久久久 | 久久精品亚洲熟妇少妇任你| 一级片免费观看大全| 女同久久另类99精品国产91| 制服诱惑二区| 韩国精品一区二区三区| 91国产中文字幕| 50天的宝宝边吃奶边哭怎么回事| 国产日韩欧美视频二区| 精品一区二区三区视频在线观看免费 | 动漫黄色视频在线观看| 欧美日本中文国产一区发布| 亚洲精品乱久久久久久| 欧美日韩黄片免| 免费女性裸体啪啪无遮挡网站| 99国产精品一区二区三区| 亚洲专区字幕在线| av天堂在线播放| 人人妻人人澡人人爽人人夜夜| 男女免费视频国产| 久久这里只有精品19| 色综合婷婷激情| 91字幕亚洲| 性高湖久久久久久久久免费观看| a在线观看视频网站| 国产精品久久电影中文字幕 | 午夜福利一区二区在线看| 啦啦啦 在线观看视频| 99精品久久久久人妻精品| 久久国产精品人妻蜜桃| 制服人妻中文乱码| 午夜免费鲁丝| 啦啦啦中文免费视频观看日本| 精品久久久精品久久久| 成年版毛片免费区| 另类精品久久| 亚洲精品中文字幕在线视频| 变态另类成人亚洲欧美熟女 | 在线观看免费日韩欧美大片| 黑人猛操日本美女一级片| 色播在线永久视频| 成人永久免费在线观看视频 | 在线永久观看黄色视频| 国产伦人伦偷精品视频| 一区二区三区精品91| 亚洲精品乱久久久久久| 日本av手机在线免费观看| 99九九在线精品视频| 在线 av 中文字幕| www日本在线高清视频| 久久精品国产亚洲av香蕉五月 | 丝袜在线中文字幕| 自线自在国产av| 久久av网站| 老鸭窝网址在线观看| 欧美精品一区二区免费开放| 久久久久久亚洲精品国产蜜桃av| 色播在线永久视频| 新久久久久国产一级毛片| 黄色视频不卡| 午夜福利在线观看吧| 欧美国产精品一级二级三级| 亚洲va日本ⅴa欧美va伊人久久| 精品亚洲成国产av| 精品乱码久久久久久99久播| 久久免费观看电影| 亚洲avbb在线观看| 又紧又爽又黄一区二区| 日韩精品免费视频一区二区三区| 99精品欧美一区二区三区四区| 成人亚洲精品一区在线观看| 国产免费福利视频在线观看| 国产精品欧美亚洲77777| 久久久久网色| 欧美老熟妇乱子伦牲交| 最近最新中文字幕大全免费视频| 久久久精品94久久精品| 啦啦啦在线免费观看视频4| 成人影院久久| 在线观看66精品国产| 日韩欧美一区视频在线观看| 精品国产乱码久久久久久男人| 又大又爽又粗| 国产一区二区在线观看av| 大码成人一级视频| 中文字幕av电影在线播放| 中文字幕精品免费在线观看视频| 自线自在国产av| 大香蕉久久网| 狂野欧美激情性xxxx| 丰满迷人的少妇在线观看| 国产精品久久久久久精品电影小说| 天天躁夜夜躁狠狠躁躁| 久久久久精品人妻al黑| 中文字幕色久视频| 不卡一级毛片| 香蕉丝袜av| 亚洲国产中文字幕在线视频| 久久久久久久大尺度免费视频| 这个男人来自地球电影免费观看| netflix在线观看网站| 美女高潮喷水抽搐中文字幕| 精品午夜福利视频在线观看一区 | 人妻 亚洲 视频| 国产精品九九99| 国产精品亚洲av一区麻豆| 亚洲av成人一区二区三| 日韩欧美一区视频在线观看| 纵有疾风起免费观看全集完整版| 后天国语完整版免费观看| 99九九在线精品视频| 国产亚洲精品第一综合不卡| 美女高潮喷水抽搐中文字幕| 免费人妻精品一区二区三区视频| 一级毛片女人18水好多| 又紧又爽又黄一区二区| 午夜免费成人在线视频| 十八禁网站网址无遮挡| 久9热在线精品视频| 丝袜美足系列| 他把我摸到了高潮在线观看 | netflix在线观看网站| 国产欧美日韩一区二区三| 三级毛片av免费| 日韩一区二区三区影片| 视频在线观看一区二区三区| 亚洲色图综合在线观看| 亚洲视频免费观看视频| 波多野结衣一区麻豆| 我要看黄色一级片免费的| 十八禁网站免费在线| 欧美一级毛片孕妇| 丝袜在线中文字幕| 天堂中文最新版在线下载| 国产一区二区 视频在线| 新久久久久国产一级毛片| 18在线观看网站| 国产在视频线精品| 欧美老熟妇乱子伦牲交| 热99re8久久精品国产| 欧美国产精品va在线观看不卡| 18禁国产床啪视频网站| aaaaa片日本免费| 国产精品99久久99久久久不卡| 国产在线视频一区二区| 精品久久蜜臀av无| 飞空精品影院首页| 亚洲欧美一区二区三区黑人| 黑人欧美特级aaaaaa片| 精品亚洲成国产av| 国产黄频视频在线观看| 一本综合久久免费| 亚洲精品粉嫩美女一区| 老司机午夜福利在线观看视频 | 亚洲精品成人av观看孕妇| 久久久国产精品麻豆| 精品少妇内射三级| 国产伦理片在线播放av一区| 久久久精品国产亚洲av高清涩受| 欧美日韩av久久| 老司机影院毛片| 国产日韩一区二区三区精品不卡| 在线观看人妻少妇| 如日韩欧美国产精品一区二区三区| 在线观看舔阴道视频| 一区二区日韩欧美中文字幕| 免费不卡黄色视频| 麻豆乱淫一区二区| 天堂中文最新版在线下载| 精品国产乱子伦一区二区三区| 欧美精品高潮呻吟av久久| 欧美精品一区二区大全| 亚洲精品av麻豆狂野| 男女午夜视频在线观看| 又大又爽又粗| 美女福利国产在线| 久久国产精品男人的天堂亚洲| 麻豆国产av国片精品| 我要看黄色一级片免费的| kizo精华| 成人国产av品久久久| 欧美日韩亚洲高清精品| 国产色视频综合| 中文亚洲av片在线观看爽 | kizo精华| 亚洲成国产人片在线观看| 97人妻天天添夜夜摸| 视频区欧美日本亚洲| 亚洲中文日韩欧美视频| 老司机影院毛片| 91九色精品人成在线观看| 日韩视频一区二区在线观看| 搡老熟女国产l中国老女人| 肉色欧美久久久久久久蜜桃| 欧美性长视频在线观看| 人成视频在线观看免费观看| 久久国产精品影院| 欧美成狂野欧美在线观看|