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

    熱堆中釷鈾轉(zhuǎn)化規(guī)律

    2015-12-02 11:39:40張海青曹長(zhǎng)青朱天寶朱智勇
    核技術(shù) 2015年5期
    關(guān)鍵詞:熱中子半衰期核素

    張海青 林 俊 曹長(zhǎng)青 朱天寶 朱智勇

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

    熱堆中釷鈾轉(zhuǎn)化規(guī)律

    張海青 林 俊 曹長(zhǎng)青 朱天寶 朱智勇

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

    釷鈾燃料循環(huán)以其優(yōu)異的物理和化學(xué)特性,受到核能界的廣泛關(guān)注。本文利用單群的點(diǎn)燃耗計(jì)算程序ORIGEN,分別研究了釷燃料在沸水堆(Boiling Water Reactor, BWR)、壓水堆(Pressurized Water Reactor, PWR)和加拿大重水鈾反應(yīng)堆(Canada Deuterium Oxide Uranium, CANDU,又稱坎杜堆)能譜中輻照時(shí),232Th、233Th、233Pa、233U等核素生成量隨中子注量率和中子能譜的變化規(guī)律,并探索了多次“輻照-冷卻”循環(huán)對(duì)釷鈾轉(zhuǎn)化率的影響。計(jì)算結(jié)果表明,能譜相同時(shí),233Th和233Pa存量的最大值與注量率有關(guān);233U存量的最大值與注量率無(wú)關(guān),大概在注量(注量率×?xí)r間)為4×1016n·cm?2左右;注量率相同時(shí),能譜越硬,233U存量的最大值越大。采取循環(huán)“輻照-冷卻”可以提高233Th-233U的轉(zhuǎn)化率,對(duì)于相同的總輻照時(shí)間,每次循環(huán)周期內(nèi)的輻照時(shí)間越短,相對(duì)于總輻照時(shí)間相同的單次輻照,轉(zhuǎn)化率增量提高越明顯;當(dāng)總輻照時(shí)間超過(guò)兩個(gè)月時(shí),循環(huán)輻照對(duì)轉(zhuǎn)化率增量的作用較小,與單次輻照不冷卻相比,轉(zhuǎn)化率相對(duì)增量不超過(guò)1倍。

    中子注量率,中子能譜,釷鈾轉(zhuǎn)化率

    232Th的天然同位素豐度接近100%,它與238U一樣,都是可增殖核材料。近年來(lái),隨著鈾資源消耗量的急劇增加,釷鈾燃料循環(huán)研究越來(lái)越為世界核能界所關(guān)注。我國(guó)釷資源相對(duì)豐富,據(jù)統(tǒng)計(jì),國(guó)內(nèi)釷資源約為鈾資源的6倍[1]。所以,中國(guó)科學(xué)院上海應(yīng)用物理研究所啟動(dòng)了“未來(lái)先進(jìn)核裂變能-釷基熔鹽核能系統(tǒng)”(Thorium Molten Salt Reactor, TMSR)固態(tài)燃料釷基熔鹽堆的建設(shè)。該反應(yīng)堆將采用高溫氣冷堆的球形燃料元件,堆芯由釷球和鈾球混合組成。

    近年來(lái)世界各國(guó)很重視釷的利用[2]。美國(guó)、法國(guó)、德國(guó)、加拿大等都開(kāi)展了不少研究工作,取得一定進(jìn)展,其中包括釷在動(dòng)力堆、非傳統(tǒng)動(dòng)力堆[3?4]中的應(yīng)用,釷的應(yīng)用方法[5?6]和燃料循環(huán)技術(shù)等[7?9]。但都是在特定能譜下進(jìn)行的研究。由于中子能譜、中子注量率不同,釷鈾轉(zhuǎn)化情況就不同。本文對(duì)二氧化釷燃料球在燃燒過(guò)程中相關(guān)核素的釷鈾轉(zhuǎn)化規(guī)律進(jìn)行了研究,包括不同熱堆模型下,中子能譜、中子注量率、輻照方式等對(duì)釷鈾轉(zhuǎn)化過(guò)程中相關(guān)核素生成規(guī)律的影響,為將來(lái)釷球燃料元件在堆中的放置方法、增殖周期、添料換料和后處理周期的確定提供了參考數(shù)據(jù)。

    1 計(jì)算工具和模型

    計(jì)算工具為ORIGEN[10],它是由美國(guó)橡樹(shù)嶺國(guó)家實(shí)驗(yàn)室(Oak Ridge National Laboratory, ORNL)研發(fā)的一個(gè)多功能點(diǎn)燃耗及放射性衰變計(jì)算機(jī)程序。該程序能模擬核燃料循環(huán)過(guò)程中燃料與中子反應(yīng)后,各放射性物質(zhì)的積累、衰變,核素或元素的添加和元素的連續(xù)移除等各種處理過(guò)程,其中包含1119種獨(dú)立裂變產(chǎn)物,129種錒系核素,698種結(jié)構(gòu)材料或活化材料數(shù)據(jù),包含衰變庫(kù)和光子數(shù)據(jù)庫(kù)。本文所選的數(shù)據(jù)庫(kù)分別是ORIGEN中的能量為0.0253 eV的熱中子數(shù)據(jù)庫(kù)(THERMAL.LIB)、壓水堆數(shù)據(jù)庫(kù)(PWRU.LIB)、沸水堆(BWRU.LIB)和重水堆(CANDUSEU.LIB)核素中子反應(yīng)截面數(shù)據(jù)庫(kù)。其中PWRU和BWRU的燃料為低富集度鈾,卸載燃耗為33.0 MWd?kgHM?1和27.5 MWd?kgHM?1,CANDUSEU的燃料為1.2 wt%低富集度濃縮鈾(Slightly Enriched Uranium, SEU)。

    計(jì)算模型是直徑6 cm的釷燃料球[11?12](含釷量為6.5 g,232Th富集度為100%),該燃料球由燃料區(qū)和非燃料區(qū)組成,燃料區(qū)半徑2.5 cm,是由釷燃料的三結(jié)構(gòu)同向性型(Tri-structural Iso-tropic, TRISO)包覆顆粒均勻分布在石墨中組成的;非燃料區(qū)的內(nèi)、外半徑分別為2.5 cm和3 cm,材料為石墨。TRISO包覆顆粒由核心和4層包覆層組成,其幾何參數(shù)和成分如表1所示。

    表1 TRISO包覆燃料顆粒設(shè)計(jì)參數(shù)Table 1 Design parameters of TRISO coated fuel particles.

    2 結(jié)果與討論

    2.1 中子注量率對(duì)釷鈾轉(zhuǎn)化相關(guān)核素存量的影響

    為研究注量率對(duì)釷鈾轉(zhuǎn)化的影響,本文選定4種注量率:1×1012n·cm?2·s?1、1×1013n·cm?2·s?1、5×1013n·cm?2·s?1、1×1014n·cm?2·s?1;運(yùn)用同一截面庫(kù)THERMAL.LIB分別計(jì)算了在這4種不同注量率下釷球燃料元件中232Th、233Th、233Pa、233U存量的變化情況。

    圖1給出了232Th的核反應(yīng)鏈,其反應(yīng)主要為(n,γ)、(n,2n)、α衰變。其貝特曼方程可寫(xiě)為:

    在熱中子輻照下,(n,2n)反應(yīng)可忽略不計(jì),由于232Th的半衰期非常長(zhǎng),所以α衰變也可忽略不計(jì)。最后解方程(1),可得232Th隨時(shí)間的變化關(guān)系式(2)。即232Th隨時(shí)間和注量率按指數(shù)關(guān)系下降。圖2為ORIGEN計(jì)算的232Th質(zhì)量隨時(shí)間的變化情況。

    圖1 232Th的核反應(yīng)鏈Fig.1 Nuclear chain reaction of 232Th.

    圖2 不同注量率下232Th質(zhì)量隨時(shí)間的變化Fig.2 Inventory of 232Th vs. time at different neutron fluence rates.

    圖3為233Th的存量占最初232Th (232Th0)摩爾數(shù)的份額隨時(shí)間的變化。由圖3可知,233Th在120 min左右達(dá)到最大值,其最大值主要決定于中子注量率的大小。而達(dá)到最大值的時(shí)間與注量率大小無(wú)關(guān)。這主要是因?yàn)椋?33Th的生成主要由232Th的(n,γ)反應(yīng)而來(lái),所以注量率大時(shí)生成量就大。其消耗主要為自身的β衰變,與注量率大小無(wú)關(guān),所以注量率大時(shí),峰值就大。但半衰期非常短,僅為22.3 min,累積質(zhì)量受時(shí)間的影響很小。所以不同注量率下達(dá)到最大值的時(shí)間基本相同。達(dá)到最大值后,隨著232Th的消耗,峰值也緩慢地減小。

    圖4是233Pa的存量占最初232Th(232Th0)摩爾數(shù)的份額隨時(shí)間的變化。由圖4可知,233Pa與233Th的變化趨勢(shì)相同,這也是因?yàn)槠浒胨テ诒容^短的原因(T1/2=26.95 d)。當(dāng)中子注量率大時(shí),233Pa峰值也高。由圖1可見(jiàn),233Pa主要來(lái)自233Th的β衰變,233Th含量大時(shí),233Pa生成量也較多。對(duì)于消耗部分,由于233Pa的自身β衰變比其(n,γ)反應(yīng)截面要大的多,所以233Pa的消耗主要決定于自身衰變,其與注量率大小也基本無(wú)關(guān)。

    圖3 不同注量率下233ThInv/232Th0的變化Fig.3 233ThInv/232Th0 vs. time at different neutron fluence rates.

    圖4 不同注量率下233PaInv/232Th0的變化Fig.4 233PaInv/232Th0 vs. time at different neutron fluence rates.

    由圖5可見(jiàn),注量率越高時(shí),233U達(dá)到飽和的時(shí)間越短。但是同一能譜下,233U峰值大小受注量率大小影響很小。為證明這個(gè)結(jié)論,本文專(zhuān)門(mén)做了在與圖5相同的4種注量率下,中子注量相同時(shí)233U的生成情況,如圖6所示,可見(jiàn)不同中子注量率時(shí)233U生成量幾乎在同一注量下達(dá)到峰值。

    圖5 不同注量率下233UInv/232Th0的變化Fig.5 233UInv/232Th0 vs. time at different neutron fluence rates.

    圖6 相同輻照注量下233UInv/232Th0的變化Fig.6 233UInv/232Th0 vs. time at the same fluencies.

    綜上所述,在反應(yīng)堆中利用釷,233Th和233Pa生成量的最大值與注量率有關(guān),而233U的最大值與注量率無(wú)關(guān),大概在4×1016n·cm?2左右。這也是想要提取最大質(zhì)量233U時(shí)的熱中子注量。

    2.2 中子能譜對(duì)釷鈾轉(zhuǎn)化相關(guān)核素存量的影響

    在中子注量率4×1014n·cm?2·s?1下,采用不同中子截面數(shù)據(jù)庫(kù),計(jì)算233Pa和233U的存量占初始232Th(232Th0)摩爾數(shù)的份額隨時(shí)間的變化,結(jié)果如圖7所示。

    在不同能譜下,233Pa生成量與233U生成規(guī)律不同。233Pa在全熱中子譜下,所達(dá)到的飽和值最大;而能譜較硬的壓水堆和沸水堆,飽和值較小。不同數(shù)據(jù)庫(kù)中,232Th的中子俘獲反應(yīng)截面不同,而233Th半衰期只有22.3 min,很快能達(dá)到飽和值,所以232Th的中子俘獲截面(n,γ)大小是影響233Pa最大質(zhì)量的主要因素。由表2可知,能譜越軟,232Th的(n,γ)反應(yīng)截面越大。當(dāng)233Pa質(zhì)量達(dá)到最大后,使233Pa生成速率減小,233Pa飽和值也慢慢減小。

    在如圖7(b)所示的注量率下,熱中子數(shù)據(jù)庫(kù)和重水堆數(shù)據(jù)庫(kù)計(jì)算結(jié)果表明,233U生成量在300 d左右達(dá)到最大值,其中重水堆和熱中子庫(kù)的最大值分別約為1.5%和1.2%。而壓水堆和沸水堆的233U質(zhì)量在1 200 d左右達(dá)到最大值,取值在3.5%?4%。這里,采用熱中子截面庫(kù)計(jì)算的233U飽和值最小,因233U與熱中子發(fā)生(n,f)反應(yīng)截面(表2)最大,比CANDUSEU.LIB大20 b左右,而生成233U的232Th(n,γ)反應(yīng),THERMAL.LIB庫(kù)截面只比CANDUSEU大2 b左右,致使CANDUSEU.LIB比THERMAL.LIB數(shù)據(jù)庫(kù)下233U的飽和值稍大。

    圖7 不同能譜時(shí)233PaInv/232Th0 (a)和233UInv/232Th0 (b)的變化Fig.7 233PaInv/232Th0 (a) and 233UInv/232Th0 (b) vs. time at different neutron spectra.

    表2 不同截面數(shù)據(jù)庫(kù)相關(guān)核素與中子反應(yīng)的有效截面Table 2 Effective cross-section of neutron-induced reactions in different database.

    由此可見(jiàn),中子能譜越軟,233U/232Th0在越短時(shí)間內(nèi)達(dá)到最大,且最大值較小。要想提取最大質(zhì)量的233U,需要在能譜較硬的地方輻照到飽和,而要想使更多的233U為有效增值因子(keff)做貢獻(xiàn),即消耗較多的233U,則需要在較軟的地方輻照。

    2.3 循環(huán)輻照對(duì)釷鈾轉(zhuǎn)化率的影響

    釷鈾轉(zhuǎn)化率主要是針對(duì)232Th轉(zhuǎn)化生成233U的總量來(lái)說(shuō)的,包括其存量和消耗量,定義如式(3):

    運(yùn)用同一中子能譜(THERMAL.LIB),采用不同中子注量率,連續(xù)輻照下232Th-233U轉(zhuǎn)化率的計(jì)算結(jié)果見(jiàn)圖8。從圖8中可以看出,連續(xù)輻照相同時(shí)間,在相同能譜下,232Th-233U的轉(zhuǎn)化率與注量率幾乎成正比。注量在1.90×1022n·cm?2內(nèi),其轉(zhuǎn)化率不超過(guò)12.7%,這表明一次通過(guò)式燃料輻照方案釷的利用率很低,所以需要尋求一種增加釷-鈾轉(zhuǎn)化率的輻照方案。

    圖8 232Th-233U的轉(zhuǎn)化率Fig.8 Conversion rate of 232Th to 233U.

    由圖1可知,233U由233Pa衰變而來(lái),如果在輻照適當(dāng)時(shí)間后,將燃料從堆中拿出來(lái),就減少了233Pa與中子發(fā)生的其它反應(yīng)路徑的消耗,233Pa將可以完全衰變?yōu)?33U,從而增加釷-鈾轉(zhuǎn)化率。所以本節(jié)采用燃料輻照模式為“輻照-冷卻-輻照-冷卻……”的開(kāi)循環(huán)模式,對(duì)循環(huán)過(guò)程中轉(zhuǎn)化率的影響因素進(jìn)行分析,包括注量率、冷卻時(shí)間、輻照時(shí)間、循環(huán)次數(shù)等,以求得到最佳的輻照循環(huán)方案。

    2.3.1 注量率大小的影響

    以釷燃料球每次輻照兩個(gè)月,冷卻233Pa的5個(gè)半衰期的循環(huán)方式循環(huán)6次為例,將不同注量率下的轉(zhuǎn)化率與單次輻照相同時(shí)間時(shí)的轉(zhuǎn)化率相比,結(jié)果值見(jiàn)表3。在第6次循環(huán)結(jié)束時(shí),若注量率相差4倍(5×1013n·cm?2·s?1和2×1014n·cm?2·s?1),轉(zhuǎn)化率和轉(zhuǎn)化率的增量也相差4倍。所以,注量率對(duì)轉(zhuǎn)化率及轉(zhuǎn)化率增量(循環(huán)相對(duì)于不循環(huán)時(shí)的增量)的影響與注量率的變化成正比。

    表3 循環(huán)輻照轉(zhuǎn)化率及其相對(duì)于單次輻照的轉(zhuǎn)化率增量Table 3 Conversion rate of cyclic irradiation, in comparison with that of continuous irradiation (n·cm?2·s?1).

    2.3.2 冷卻時(shí)間的影響

    冷卻時(shí)間越長(zhǎng),由233Pa轉(zhuǎn)化到233U的量就越大。但從實(shí)用性考慮,不能無(wú)限長(zhǎng),要選擇一個(gè)合適的冷卻時(shí)間。冷卻4個(gè)半衰期的時(shí)候,核素的衰變率已達(dá)到96.9%,于是我們選擇4個(gè)、5個(gè)、7個(gè)半衰期做比較,以注量率為5×1013n·cm?2·s?1為例,輻照兩個(gè)月后,分別求出衰變4個(gè)、5個(gè)半衰期時(shí),232Th-233U的轉(zhuǎn)化率與衰變7個(gè)半衰期時(shí)的差別,將該數(shù)值列于表4。冷卻時(shí)間選4或5個(gè)半衰期即可達(dá)到與最大值非常相近的值。所以,可根據(jù)實(shí)際情況選擇冷卻4個(gè)或5個(gè)半衰期。

    表4 每次循環(huán)的衰變時(shí)間對(duì)轉(zhuǎn)化率的影響Table 4 Effect of decay time in each cycle on conversion rate.

    2.3.3 輻照時(shí)間、循環(huán)次數(shù)的影響

    圖4中,233Pa在連續(xù)輻照200 d時(shí)存量達(dá)到峰值,在前30 d的時(shí)候,濃度上升非???;30?60 d上升較快;60?100 d上升速度減慢;100?200 d曲線變得比較平緩,直到達(dá)到最大值。所以,本文分別選輻照一個(gè)月、兩個(gè)月、6個(gè)月和200 d來(lái)研究循環(huán)輻照中輻照時(shí)間對(duì)轉(zhuǎn)化率的影響。

    圖9給出注量率為5×1013n·cm?2·s?1、總輻照時(shí)間相同時(shí)(180 d),采用以下前三種不同循環(huán)方式下的轉(zhuǎn)化率:(1) 輻照一個(gè)月,冷卻5個(gè)半衰期,循環(huán)6次,共用時(shí)33個(gè)月;(2) 輻照兩個(gè)月,冷卻5個(gè)半衰期,循環(huán)三次,共用時(shí)26個(gè)月;(3) 輻照6個(gè)月,冷卻5個(gè)半衰期,循環(huán)1次,共用10.5個(gè)月。

    由圖9可知,即使循環(huán)方式不同,在圖9的增量很小,經(jīng)計(jì)算在0.01%?0.02%。這主要是因?yàn)樵趫D示注量范圍內(nèi),轉(zhuǎn)化率本來(lái)就很小,不超過(guò)0.6%。

    表5給出注量率為5×1013n·cm?2·s?1時(shí),5種循環(huán)方式下,輻照一段時(shí)間再冷卻后的值與只輻照不冷卻相比,各個(gè)輻照階段轉(zhuǎn)化率的相對(duì)增量。相對(duì)增量的計(jì)算式如下:

    圖9 總輻照時(shí)間相同時(shí)不同循環(huán)方式的轉(zhuǎn)化率—:輻照過(guò)程,---:冷卻過(guò)程Fig.9 Conversion rate of different cyclic manner at the same total irradiation time.—: Irradiation process, ---: Decay process

    表5 注量率為5×1013 n·cm?2·s?1時(shí),循環(huán)輻照相對(duì)于單次輻照232Th-233U轉(zhuǎn)化率的相對(duì)增量Table 5 Conversion rate increment of 232Th to 233U compare different cyclic manner to single irradiation with the same neutron fluence rate 5×1013 n·cm?2·s?1.

    由表5可看出,循環(huán)輻照中,每次循環(huán)的輻照時(shí)間越短,相對(duì)于相同總輻照時(shí)間的連續(xù)輻照,轉(zhuǎn)化率的增量越大,每次循環(huán)的輻照時(shí)間越長(zhǎng),轉(zhuǎn)化率相對(duì)增量越不明顯。這是因?yàn)椋?33Pa半衰期比233U短得多,輻照時(shí)間越長(zhǎng),前邊生成的233Pa已經(jīng)都衰變成233U,達(dá)到了與冷卻同樣的效果,所以轉(zhuǎn)化率增量就比短期輻照小。

    可見(jiàn),采用循環(huán)輻照冷卻的方法可以提高釷鈾轉(zhuǎn)化率,對(duì)于相同的總輻照時(shí)間,每次循環(huán)的輻照時(shí)間越短,相對(duì)于連續(xù)輻照,轉(zhuǎn)化率增量提高的倍數(shù)越多,作用越明顯。但由于圖9和表5條件下絕對(duì)轉(zhuǎn)化率本來(lái)就很小,所以采用循環(huán)輻照,與連續(xù)輻照相比絕對(duì)增量提高很小。因此,是否采取循環(huán)輻照,要根據(jù)實(shí)際情況來(lái)考慮。如果輻照時(shí)間比較長(zhǎng)(大于兩個(gè)月),就沒(méi)必要。如果輻照時(shí)間較短,循環(huán)輻照不僅可以提高轉(zhuǎn)化率,還可以節(jié)省中子注量、降低成本,所以是不錯(cuò)的選擇。

    3 結(jié)語(yǔ)

    通過(guò)用ORIGEN點(diǎn)燃耗計(jì)算程序,分析不同條件下釷鈾轉(zhuǎn)化規(guī)律,可得以下結(jié)論:

    (1) 同一能譜,注量率越大,233U生成速率越快;但233U的最大值與注量率無(wú)關(guān)。注量相同時(shí),最大值幾乎相同(圖6),大概在4×1016n·cm?2。

    (2) 與233U不同,能譜相同時(shí),233Th在120 min左右達(dá)到最大值,其最大值主要決定于中子注量率的大小。而達(dá)到最大值的時(shí)間與注量率大小無(wú)關(guān)。233Pa與233Th的變化趨勢(shì)相同,即最大值的大小決定于注量率,達(dá)到最大值的時(shí)間與注量率大小無(wú)關(guān)。

    (3) 在不同能譜下,能譜越硬,233U飽和值越大(圖7),如果要提取最大質(zhì)量的233U,應(yīng)該在能譜較硬的地方輻照;能譜越軟,消耗掉的233U越多,在反應(yīng)堆中利用233U做裂變?nèi)剂?,?yīng)該選用較軟的能譜。

    (4) 在不同能譜下,233Pa生成量與233U生成規(guī)律不同。233Pa在全熱中子譜下,所達(dá)到的飽和值最大;能譜越硬,飽和值越小。

    (5) 采取循環(huán)輻照冷卻可以提高233Th-233U的轉(zhuǎn)化率,每次循環(huán)周期內(nèi)輻照時(shí)間越短,相對(duì)增量提高的越明顯,當(dāng)輻照時(shí)間超過(guò)兩個(gè)月時(shí),循環(huán)輻照對(duì)轉(zhuǎn)化率增量的作用較小,與單次輻照不冷卻相比,轉(zhuǎn)化率相對(duì)增量不超過(guò)1倍。

    1 張家驊, 包伯榮. 我國(guó)釷資源調(diào)查及釷-鈾燃料循環(huán)研究[A]. 張家驊. 釷-鈾核燃料循環(huán)研究[C]. 上海: 中國(guó)科學(xué)院上海應(yīng)用物理研究所, 2009: 22?26

    ZHANG Jiahua, BAO Borong. Investigation of thorium resources in China and study of thorium-uranium fuel cycle[A]. ZHANG Jiahua. Study of thorium-uranium nuclear fuel cycle[C]. Shanghai: Shanghai Institute of Applied Physics, Chinese Academy of Sciences, 2009: 22?26

    2 Troyanov M F. Thorium fuel utilization: options and trends[R]. Vienna: IAEA-TECDOC-1319, 2002

    3 Furukawa K, Lecocq A, Kato Y, et al. Thorium molten-salt nuclear energy synergetics[J]. Journal of Nuclear Science and Technology, 1990, 27(12): 1157?1178

    4 Furukawa K, Mitachi K, Kato Y, et al. Small molten-salt reactors with a rational thorium fuel-cycle[J]. Nuclear Engineering and Design, 1992, 136(1): 157?165

    5 Murogov V M, Zinin A J, Ilyunin V G, et al. Fast reactors with different fuels in uranium-plutonium and combined fuel cycle[G]. IPPE-1920, 1998, 23(1): 203?206

    6 Rutten H J, Teuchert E. Advanced safety features of pebble bed HTR's with Th-utilization[R]. Julich: KFA, 1993

    7 Yu J Y, Wang K, Sollychin R, et al. Thorium fuel cycle of a thorium-based advanced nuclear energy system[J]. Progress in Nuclear Energy, 2004, 45(1): 71?84

    8 Chang J J, Chang J P, Won I K. Dynamic analysis of a thorium fuel cycle in CANDU reactors[J]. Annals of Nuclear Energy, 2008, 35(10): 1842–1848

    9 Yamamoto T, Suwarno H, Kayano H, et al. Development of new reactor fuel materials[J]. Journal of Nuclear Science and Technology, 1995, 32(3): 260?262

    10 Croff A G. A user's manual for the ORIGEN2 computer code[M]. USA: Oak Ridge National Laboratory, 1980

    11 熊文綱, 李文新, 王敏. 基于熱堆的釷鈾轉(zhuǎn)換過(guò)程中232U生成的模擬計(jì)算[J]. 核技術(shù), 2012, 35(5): 395?400

    XIONG Wengang, LI Wenxin, WANG Min. Simulative calculation of232U productions in thorium-uranium transform process based on thermal reactor[J]. Nuclear Techniques, 2012, 35(5): 395?400

    12 尹文靜, 張海青, 曹長(zhǎng)青, 等. TRISO釷鈾包覆燃料顆粒裂變氣體生成規(guī)律研究[J]. 核技術(shù), 2014, 37(1): 010601

    YIN Wenjing, ZHANG Haiqing, CAO Changqing, et al. Investigation of fission product gas formation in TRISO Th/U coated fuel particles[J]. Nuclear Techniques, 2014, 37(1): 010601

    CLC TL249

    Simulation study on232Th-233U conversion in thermal reactors

    ZHANG Haiqing LIN Jun CAO Changqing ZHU Tianbao ZHU Zhiyong

    (Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Jiading Campus, Shanghai 201800, China)

    Background: Thorium-uranium fuel cycle is attracting more and more attention because of its unique physical and chemical characteristics. Purpose: Realizing the utilization of thorium fuel in thermal reactors can save the valuable natural uranium resources and produce more fissile fuel, which is conducive to the nuclear energy sustainable development. Methods: ORIGEN code was used to simulate the build-up characteristics of232Th,233Th,233Pa and233U in neutron fluency rates and neutron spectra of typical Boiling Water Reactor (BWR), Pressurized Water Reactor (PWR) and Canada Deuterium Oxide Uranium (CANDU) reactor, and the effects of multiple“irradiation-cooling” cycles on the thorium uranium conversion rate under various irradiation time were analyzed by numerical comparison. Results: Simulation results showed that the maximum inventory of233Th and233Pa is irrelevant to the neutron fluency rates when neutron spectrum is fixed. When the neutron fluency rates are preset, the harder the neutron spectrum is, the larger the maximum inventory of233U will be obtained. Multiple“irradiation-cooling” cycles can improve the conversion rate of232Th-233U, but the relative increment is becoming smaller and smaller, compared to the continuous irradiation conversion rate. Conclusion: This work provides relevant theoretical basis for thermal reactors thorium-uranium fuel cycle research.

    Neutron fluence rate, Neutron spectrum, Conversion rate of232Th-233U

    TL249

    10.11889/j.0253-3219.2015.hjs.38.050601

    中國(guó)科學(xué)院戰(zhàn)略性先導(dǎo)科技專(zhuān)項(xiàng)(No.XDA02030200)資助

    張海青,女,1981年出生,2011年于中國(guó)原子能科學(xué)研究院獲博士學(xué)位,從事釷鈾燃料技術(shù)研究,E-mail: zhanghaiqing@sinap.ac.cn

    2014-12-09,

    2015-03-05

    猜你喜歡
    熱中子半衰期核素
    熱中子透射成像轉(zhuǎn)換屏物理設(shè)計(jì)研究
    核素分類(lèi)開(kāi)始部分的6種7核素小片分布
    核素分類(lèi)的4量子數(shù)
    基于語(yǔ)言學(xué)中文學(xué)術(shù)圖書(shū)的半衰期分析研究*
    單晶硅受照熱中子注量率的雙箔活化法測(cè)量研究
    基于引用半衰期的我國(guó)五官學(xué)期刊文獻(xiàn)老化研究
    基于JCR?的國(guó)外臨床醫(yī)學(xué)學(xué)科半衰期
    脈沖中子-裂變中子探測(cè)鈾黃餅的MCNP模擬
    植物對(duì)核素鍶的吸附與富集作用研究現(xiàn)狀
    核素分布的三循環(huán)
    91久久精品电影网| 男女视频在线观看网站免费| 非洲黑人性xxxx精品又粗又长| 国产精品永久免费网站| 日日摸夜夜添夜夜添av毛片| 国产精品综合久久久久久久免费| 亚州av有码| 国产成人精品一,二区| 蜜桃亚洲精品一区二区三区| 永久免费av网站大全| 国产熟女欧美一区二区| 你懂的网址亚洲精品在线观看 | 亚洲自拍偷在线| 国产精品福利在线免费观看| 国产一级毛片七仙女欲春2| 亚洲精品久久久久久婷婷小说 | 亚洲aⅴ乱码一区二区在线播放| 国产精华一区二区三区| 国产精品美女特级片免费视频播放器| 精品人妻熟女av久视频| 1024手机看黄色片| 国产在视频线精品| 亚洲av不卡在线观看| 国产午夜精品一二区理论片| 亚洲精品乱码久久久久久按摩| 亚洲欧美成人综合另类久久久 | .国产精品久久| 国产真实伦视频高清在线观看| a级毛片免费高清观看在线播放| 亚洲精品国产av成人精品| 极品教师在线视频| 久久久久国产网址| 2021天堂中文幕一二区在线观| 国产av不卡久久| 日韩欧美精品免费久久| 精品人妻一区二区三区麻豆| 美女脱内裤让男人舔精品视频| 免费播放大片免费观看视频在线观看 | 一边亲一边摸免费视频| 午夜精品在线福利| 国产免费又黄又爽又色| 白带黄色成豆腐渣| 在线观看一区二区三区| 99国产精品一区二区蜜桃av| 亚洲精品一区蜜桃| 在线天堂最新版资源| 草草在线视频免费看| 成年女人永久免费观看视频| 三级毛片av免费| 免费观看精品视频网站| 国内少妇人妻偷人精品xxx网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日本色播在线视频| 亚洲最大成人av| 人妻少妇偷人精品九色| 秋霞伦理黄片| 波多野结衣高清无吗| 高清av免费在线| 免费观看的影片在线观看| 老师上课跳d突然被开到最大视频| .国产精品久久| 精品一区二区免费观看| 亚洲综合精品二区| 久久精品熟女亚洲av麻豆精品 | 欧美日韩在线观看h| 欧美性感艳星| 中文字幕免费在线视频6| 日韩,欧美,国产一区二区三区 | 色尼玛亚洲综合影院| 国产精品人妻久久久久久| 一个人看视频在线观看www免费| 男女国产视频网站| 久久6这里有精品| 18禁裸乳无遮挡免费网站照片| 亚洲欧美清纯卡通| 中文字幕免费在线视频6| 国产免费男女视频| 国产精品一区二区在线观看99 | 久久婷婷人人爽人人干人人爱| 日韩欧美精品v在线| 九九久久精品国产亚洲av麻豆| 国内揄拍国产精品人妻在线| 有码 亚洲区| 久久午夜福利片| 秋霞在线观看毛片| 国产麻豆成人av免费视频| 久久久久久久国产电影| 青春草视频在线免费观看| 久久久久久久久中文| 岛国毛片在线播放| 国产午夜福利久久久久久| 身体一侧抽搐| 久久亚洲精品不卡| 国产伦精品一区二区三区视频9| 男女下面进入的视频免费午夜| 精品午夜福利在线看| 视频中文字幕在线观看| 国产亚洲精品久久久com| 免费看日本二区| 免费人成在线观看视频色| 一个人看的www免费观看视频| 亚洲激情五月婷婷啪啪| 国产精品久久久久久精品电影| 99九九线精品视频在线观看视频| 99久国产av精品国产电影| 中文字幕亚洲精品专区| 久久这里有精品视频免费| 波野结衣二区三区在线| 成人三级黄色视频| 久久久精品94久久精品| 国产亚洲精品久久久com| 亚洲国产精品久久男人天堂| 日本黄大片高清| 国产黄a三级三级三级人| 在线免费观看不下载黄p国产| 中文字幕av在线有码专区| 日本五十路高清| 免费看美女性在线毛片视频| 亚洲久久久久久中文字幕| 内射极品少妇av片p| 91aial.com中文字幕在线观看| 亚洲av免费在线观看| 国产精品久久久久久精品电影小说 | 国产91av在线免费观看| 国产av码专区亚洲av| 天天躁夜夜躁狠狠久久av| av在线天堂中文字幕| 少妇人妻精品综合一区二区| 天堂av国产一区二区熟女人妻| 欧美一区二区精品小视频在线| 中文字幕熟女人妻在线| 蜜桃亚洲精品一区二区三区| 大话2 男鬼变身卡| 97超视频在线观看视频| 成年免费大片在线观看| 黄色配什么色好看| 国产黄a三级三级三级人| 在现免费观看毛片| 人妻制服诱惑在线中文字幕| 国产国拍精品亚洲av在线观看| 性色avwww在线观看| videos熟女内射| 最近中文字幕2019免费版| 亚洲国产欧洲综合997久久,| 午夜激情福利司机影院| 欧美激情国产日韩精品一区| 长腿黑丝高跟| 三级国产精品片| 久久久久久久久久久免费av| 亚州av有码| 在线观看美女被高潮喷水网站| av线在线观看网站| 亚洲第一区二区三区不卡| 婷婷六月久久综合丁香| 51国产日韩欧美| 麻豆一二三区av精品| 精品一区二区三区人妻视频| 男女边吃奶边做爰视频| 国产视频内射| 国产熟女欧美一区二区| 丰满乱子伦码专区| 青春草国产在线视频| 欧美不卡视频在线免费观看| 色吧在线观看| 国产麻豆成人av免费视频| 日韩,欧美,国产一区二区三区 | 午夜福利成人在线免费观看| 人人妻人人澡人人爽人人夜夜 | 欧美日韩精品成人综合77777| 日本爱情动作片www.在线观看| 99久国产av精品| 亚洲最大成人中文| 国产私拍福利视频在线观看| 午夜老司机福利剧场| 青春草视频在线免费观看| 国产高潮美女av| 欧美性感艳星| 夜夜爽夜夜爽视频| 自拍偷自拍亚洲精品老妇| 国产精品久久视频播放| 亚洲在线观看片| 日日啪夜夜撸| 日日摸夜夜添夜夜爱| 啦啦啦观看免费观看视频高清| 日产精品乱码卡一卡2卡三| 我要搜黄色片| 午夜精品一区二区三区免费看| 永久免费av网站大全| 天美传媒精品一区二区| 国产亚洲av片在线观看秒播厂 | av在线观看视频网站免费| 美女内射精品一级片tv| 国产淫语在线视频| 亚洲欧美精品专区久久| 两个人视频免费观看高清| 精品久久久久久成人av| 国产精品一区二区在线观看99 | 26uuu在线亚洲综合色| 99久久无色码亚洲精品果冻| av专区在线播放| 久久精品久久久久久噜噜老黄 | 欧美不卡视频在线免费观看| 国产精品人妻久久久影院| 久久99热6这里只有精品| 一个人看视频在线观看www免费| www日本黄色视频网| 亚洲国产精品成人综合色| 免费看日本二区| 中文字幕av在线有码专区| 99热全是精品| 国产精品久久久久久av不卡| 亚洲va在线va天堂va国产| 免费看日本二区| 看黄色毛片网站| 国产成人精品久久久久久| 99热网站在线观看| 性色avwww在线观看| 直男gayav资源| 免费一级毛片在线播放高清视频| 老司机影院毛片| 男女啪啪激烈高潮av片| 中文字幕熟女人妻在线| 欧美xxxx性猛交bbbb| 国产成人午夜福利电影在线观看| ponron亚洲| 99在线人妻在线中文字幕| 免费黄色在线免费观看| 亚洲国产日韩欧美精品在线观看| 人妻夜夜爽99麻豆av| 国产国拍精品亚洲av在线观看| 一级毛片我不卡| 级片在线观看| 国语自产精品视频在线第100页| 黄片无遮挡物在线观看| 一边亲一边摸免费视频| 免费观看的影片在线观看| 又爽又黄a免费视频| 亚洲性久久影院| 日日啪夜夜撸| 97在线视频观看| 中文乱码字字幕精品一区二区三区 | 婷婷色麻豆天堂久久 | 如何舔出高潮| 日韩视频在线欧美| 99视频精品全部免费 在线| 男人舔女人下体高潮全视频| 免费大片18禁| 亚洲欧美中文字幕日韩二区| 国产黄色视频一区二区在线观看 | 日本免费在线观看一区| 精品人妻偷拍中文字幕| 婷婷六月久久综合丁香| 最近的中文字幕免费完整| 久久鲁丝午夜福利片| 欧美一区二区国产精品久久精品| 亚洲不卡免费看| 2022亚洲国产成人精品| 亚洲高清免费不卡视频| 色综合色国产| 视频中文字幕在线观看| 免费人成在线观看视频色| 99久久精品热视频| 成年免费大片在线观看| 久久久午夜欧美精品| 日本一二三区视频观看| 成人高潮视频无遮挡免费网站| 久久久精品94久久精品| 国产又黄又爽又无遮挡在线| 搡女人真爽免费视频火全软件| 亚洲综合精品二区| 欧美日韩综合久久久久久| 一级毛片久久久久久久久女| 亚洲在线自拍视频| 搡女人真爽免费视频火全软件| 国产精品久久久久久精品电影小说 | 精品熟女少妇av免费看| 乱人视频在线观看| 精品国内亚洲2022精品成人| 亚洲成人中文字幕在线播放| 日韩亚洲欧美综合| 看片在线看免费视频| 亚洲精品日韩av片在线观看| 欧美色视频一区免费| 色网站视频免费| 床上黄色一级片| 国产v大片淫在线免费观看| 成年女人看的毛片在线观看| 毛片一级片免费看久久久久| 青春草视频在线免费观看| 欧美日本亚洲视频在线播放| 亚洲av免费高清在线观看| 日本一本二区三区精品| 国产成人免费观看mmmm| 99热6这里只有精品| 黄色配什么色好看| 国产高清国产精品国产三级 | 日韩精品青青久久久久久| 成人午夜高清在线视频| 久久鲁丝午夜福利片| 久久人人爽人人爽人人片va| 欧美极品一区二区三区四区| 国产三级中文精品| 亚洲国产色片| 在线a可以看的网站| 久久99热这里只频精品6学生 | 国产伦精品一区二区三区四那| av免费在线看不卡| 一边摸一边抽搐一进一小说| 欧美不卡视频在线免费观看| 欧美性猛交╳xxx乱大交人| 啦啦啦韩国在线观看视频| 亚洲精品日韩av片在线观看| 久久草成人影院| 日产精品乱码卡一卡2卡三| 一区二区三区乱码不卡18| 看免费成人av毛片| 国产亚洲最大av| 大又大粗又爽又黄少妇毛片口| 日韩大片免费观看网站 | 午夜福利网站1000一区二区三区| 天天躁日日操中文字幕| 国产精品1区2区在线观看.| 嫩草影院精品99| 国产又黄又爽又无遮挡在线| 男女下面进入的视频免费午夜| 国产一级毛片在线| 亚洲精品乱码久久久久久按摩| 精品一区二区三区视频在线| 人人妻人人澡人人爽人人夜夜 | 九色成人免费人妻av| 久久久a久久爽久久v久久| 国产精品无大码| 少妇熟女欧美另类| 国产免费福利视频在线观看| 国产精品嫩草影院av在线观看| 精品酒店卫生间| 一卡2卡三卡四卡精品乱码亚洲| 亚洲电影在线观看av| 菩萨蛮人人尽说江南好唐韦庄 | 日本免费在线观看一区| 女人十人毛片免费观看3o分钟| 欧美+日韩+精品| 欧美精品国产亚洲| 亚洲av电影不卡..在线观看| 欧美又色又爽又黄视频| 成年免费大片在线观看| 日韩av在线免费看完整版不卡| 国产成人a∨麻豆精品| 在线播放无遮挡| 久久久久久大精品| 国产色婷婷99| 久久久久久久久大av| 伦精品一区二区三区| 免费一级毛片在线播放高清视频| 久久精品久久精品一区二区三区| 欧美极品一区二区三区四区| 日本黄色视频三级网站网址| 成人特级av手机在线观看| 在现免费观看毛片| 国产午夜精品论理片| 婷婷色综合大香蕉| 亚洲aⅴ乱码一区二区在线播放| 日韩欧美国产在线观看| 只有这里有精品99| 国产高清有码在线观看视频| 大又大粗又爽又黄少妇毛片口| 亚洲av中文av极速乱| 美女内射精品一级片tv| 久久久久久久久中文| 精品人妻视频免费看| 国产成人精品婷婷| 国产精品伦人一区二区| 亚洲欧美日韩东京热| 精品一区二区免费观看| 老司机影院成人| 久久久久国产网址| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品无人区乱码1区二区| 亚洲,欧美,日韩| 亚洲欧洲国产日韩| 日韩亚洲欧美综合| 欧美97在线视频| 3wmmmm亚洲av在线观看| 亚洲五月天丁香| 成人美女网站在线观看视频| 91久久精品电影网| 青青草视频在线视频观看| 99久国产av精品| 国产真实伦视频高清在线观看| 51国产日韩欧美| 中文亚洲av片在线观看爽| 国产黄片视频在线免费观看| kizo精华| 午夜亚洲福利在线播放| 国产精品人妻久久久久久| 成人漫画全彩无遮挡| a级一级毛片免费在线观看| 亚洲av福利一区| 十八禁国产超污无遮挡网站| 男女下面进入的视频免费午夜| 国产视频内射| 日本免费一区二区三区高清不卡| 在线观看美女被高潮喷水网站| 欧美+日韩+精品| 97在线视频观看| 亚洲av成人精品一区久久| 久久精品夜夜夜夜夜久久蜜豆| 99热全是精品| 成人三级黄色视频| 少妇熟女欧美另类| 欧美一区二区精品小视频在线| 亚洲国产最新在线播放| 亚洲成色77777| 亚洲婷婷狠狠爱综合网| 亚洲av二区三区四区| 国产国拍精品亚洲av在线观看| 天堂av国产一区二区熟女人妻| 亚洲av成人av| 小蜜桃在线观看免费完整版高清| 一级毛片电影观看 | 日韩精品有码人妻一区| 国产又黄又爽又无遮挡在线| 99久久中文字幕三级久久日本| 国产一区有黄有色的免费视频 | 精品国产三级普通话版| 乱码一卡2卡4卡精品| 国产成人精品婷婷| 中文字幕久久专区| 综合色av麻豆| 伊人久久精品亚洲午夜| 欧美成人精品欧美一级黄| 亚洲精品自拍成人| 熟女人妻精品中文字幕| 久久精品91蜜桃| 国产真实乱freesex| 日本黄色视频三级网站网址| videossex国产| 99久久九九国产精品国产免费| 乱人视频在线观看| 国内精品一区二区在线观看| 日本av手机在线免费观看| 国产伦在线观看视频一区| av在线播放精品| 亚洲性久久影院| 国产精品爽爽va在线观看网站| 一级毛片aaaaaa免费看小| 亚洲av电影在线观看一区二区三区 | 中国国产av一级| 国产激情偷乱视频一区二区| 久久精品影院6| 丝袜喷水一区| 亚洲最大成人av| 国产在线男女| 国国产精品蜜臀av免费| 久久这里有精品视频免费| 超碰av人人做人人爽久久| 欧美日本视频| 青春草亚洲视频在线观看| 国产毛片a区久久久久| 亚洲aⅴ乱码一区二区在线播放| 美女脱内裤让男人舔精品视频| 国产在视频线精品| 夫妻性生交免费视频一级片| 亚洲国产色片| 久久99蜜桃精品久久| 又爽又黄无遮挡网站| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲真实伦在线观看| 三级毛片av免费| 97人妻精品一区二区三区麻豆| 99久久成人亚洲精品观看| 亚洲av男天堂| 日韩欧美在线乱码| 国产精品三级大全| 久久精品国产亚洲av涩爱| 免费观看性生交大片5| 成人高潮视频无遮挡免费网站| 国语对白做爰xxxⅹ性视频网站| 夫妻性生交免费视频一级片| 亚洲av免费高清在线观看| 日产精品乱码卡一卡2卡三| 三级国产精品片| 草草在线视频免费看| 国产精品永久免费网站| 亚洲熟妇中文字幕五十中出| 波多野结衣巨乳人妻| av在线天堂中文字幕| 久久这里有精品视频免费| 国产黄片视频在线免费观看| 欧美不卡视频在线免费观看| 国产v大片淫在线免费观看| 蜜臀久久99精品久久宅男| 老司机福利观看| 精华霜和精华液先用哪个| 老司机福利观看| 免费观看a级毛片全部| 最近最新中文字幕大全电影3| av又黄又爽大尺度在线免费看 | 在线观看一区二区三区| 亚洲欧美精品综合久久99| 精品久久久久久久久亚洲| 国产91av在线免费观看| 亚洲综合色惰| 夜夜爽夜夜爽视频| 少妇裸体淫交视频免费看高清| 特级一级黄色大片| 一本久久精品| 久久久久久九九精品二区国产| 国产精品久久久久久av不卡| 国产成人精品久久久久久| 欧美成人一区二区免费高清观看| 亚洲中文字幕日韩| 久久久久性生活片| 2021天堂中文幕一二区在线观| 99热这里只有是精品在线观看| 久久久午夜欧美精品| 直男gayav资源| 国产精品一及| 日本一本二区三区精品| 国产高清视频在线观看网站| 欧美+日韩+精品| av在线播放精品| 99久久九九国产精品国产免费| 国产淫片久久久久久久久| 午夜福利高清视频| 国产亚洲一区二区精品| 午夜a级毛片| 久久久久久久国产电影| 精品久久久久久久久av| 精品酒店卫生间| 插阴视频在线观看视频| 亚洲精品,欧美精品| 嘟嘟电影网在线观看| 国产男人的电影天堂91| 国产精品人妻久久久久久| 少妇的逼好多水| 成人欧美大片| 久久这里有精品视频免费| 久久精品久久久久久噜噜老黄 | 亚洲人成网站在线观看播放| 久久综合国产亚洲精品| 成人美女网站在线观看视频| 只有这里有精品99| 国产精品蜜桃在线观看| 成人三级黄色视频| 成年av动漫网址| 国产成人免费观看mmmm| 亚洲人成网站在线播| 成人毛片60女人毛片免费| 99久久精品国产国产毛片| 国产精品国产高清国产av| av卡一久久| 大香蕉久久网| 亚洲婷婷狠狠爱综合网| 国产毛片a区久久久久| 我的老师免费观看完整版| 美女黄网站色视频| 中文字幕亚洲精品专区| 亚洲国产精品sss在线观看| 午夜福利成人在线免费观看| 日产精品乱码卡一卡2卡三| 水蜜桃什么品种好| 2022亚洲国产成人精品| 国产人妻一区二区三区在| 欧美最新免费一区二区三区| 2021少妇久久久久久久久久久| 日本黄大片高清| 中文字幕精品亚洲无线码一区| 日韩一本色道免费dvd| 六月丁香七月| 成人高潮视频无遮挡免费网站| 天天躁日日操中文字幕| 国产精品一区二区三区四区免费观看| 欧美精品国产亚洲| 国产色爽女视频免费观看| 亚洲真实伦在线观看| 欧美成人免费av一区二区三区| 波多野结衣高清无吗| 欧美xxxx黑人xx丫x性爽| 免费看光身美女| 久久久久久大精品| 国产色爽女视频免费观看| 最新中文字幕久久久久| 麻豆精品久久久久久蜜桃| av国产免费在线观看| 久久草成人影院| 亚洲内射少妇av| 麻豆成人av视频| 激情 狠狠 欧美| 亚洲欧美精品自产自拍| 青春草国产在线视频| 国产男人的电影天堂91| 亚洲,欧美,日韩| 成人性生交大片免费视频hd| 久久婷婷人人爽人人干人人爱| 国产精品电影一区二区三区| 精品久久久久久久末码| 国产免费一级a男人的天堂| 日本三级黄在线观看| 一级二级三级毛片免费看| 美女脱内裤让男人舔精品视频| 嘟嘟电影网在线观看| 美女黄网站色视频| 黄色一级大片看看| 最近的中文字幕免费完整| 免费播放大片免费观看视频在线观看 | 精品久久久久久成人av| 久久精品影院6| 一级黄色大片毛片| 免费看日本二区| 美女黄网站色视频| 国产精品一区二区三区四区免费观看| 麻豆乱淫一区二区| 久久久成人免费电影| 麻豆成人av视频|