• <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)
    毛片女人毛片| 亚洲经典国产精华液单| 日韩欧美精品v在线| 亚洲av成人av| 简卡轻食公司| 啦啦啦韩国在线观看视频| 欧美日韩综合久久久久久| 人人妻人人澡欧美一区二区| 中文字幕熟女人妻在线| 最近的中文字幕免费完整| 日韩大尺度精品在线看网址| 亚洲欧美精品自产自拍| 久久久久久久久久黄片| 最近手机中文字幕大全| 日本在线视频免费播放| 国产精品国产高清国产av| 欧美日本亚洲视频在线播放| 成人漫画全彩无遮挡| 久久精品夜色国产| 欧美日韩国产亚洲二区| 亚洲欧美日韩无卡精品| 国国产精品蜜臀av免费| 麻豆成人午夜福利视频| 人人妻人人看人人澡| 寂寞人妻少妇视频99o| 日日撸夜夜添| 在现免费观看毛片| 日日摸夜夜添夜夜添小说| 免费在线观看影片大全网站| 成人特级av手机在线观看| 国产片特级美女逼逼视频| 国产又黄又爽又无遮挡在线| 身体一侧抽搐| 亚洲精品粉嫩美女一区| av女优亚洲男人天堂| 午夜精品一区二区三区免费看| 亚洲中文字幕日韩| 国产av不卡久久| 日韩亚洲欧美综合| 赤兔流量卡办理| ponron亚洲| 亚洲内射少妇av| 亚洲欧美日韩高清专用| 天堂av国产一区二区熟女人妻| 日本黄大片高清| 国产成年人精品一区二区| 午夜福利视频1000在线观看| 美女被艹到高潮喷水动态| 两个人的视频大全免费| 最近视频中文字幕2019在线8| 热99re8久久精品国产| 精品久久久久久久人妻蜜臀av| 91久久精品电影网| 激情 狠狠 欧美| 韩国av在线不卡| 日日摸夜夜添夜夜添av毛片| 村上凉子中文字幕在线| 色播亚洲综合网| 欧美另类亚洲清纯唯美| 国产老妇女一区| 亚洲综合色惰| 麻豆国产av国片精品| 国产色婷婷99| 精品人妻偷拍中文字幕| 偷拍熟女少妇极品色| 国产亚洲欧美98| 搡老熟女国产l中国老女人| 免费看a级黄色片| 18禁在线播放成人免费| 一卡2卡三卡四卡精品乱码亚洲| 亚洲欧美日韩东京热| 美女大奶头视频| 国产黄色小视频在线观看| 亚洲一区高清亚洲精品| 白带黄色成豆腐渣| 精品乱码久久久久久99久播| 国产美女午夜福利| 欧美+日韩+精品| 熟女人妻精品中文字幕| 日本一本二区三区精品| ponron亚洲| 亚洲精品国产成人久久av| 国产视频内射| 99久久九九国产精品国产免费| av国产免费在线观看| 少妇猛男粗大的猛烈进出视频 | 日本在线视频免费播放| 三级国产精品欧美在线观看| 18禁在线播放成人免费| 久久精品国产鲁丝片午夜精品| 国产精品电影一区二区三区| 亚洲国产高清在线一区二区三| 日韩欧美三级三区| 久久久国产成人免费| 国产精品野战在线观看| 91av网一区二区| 一区二区三区四区激情视频 | 三级男女做爰猛烈吃奶摸视频| 久久久久久伊人网av| 亚洲精品国产av成人精品 | 99九九线精品视频在线观看视频| 男人狂女人下面高潮的视频| 黄色日韩在线| 亚洲国产欧美人成| 欧美成人免费av一区二区三区| 网址你懂的国产日韩在线| 久久久久久伊人网av| 久久午夜亚洲精品久久| 欧美+日韩+精品| 色噜噜av男人的天堂激情| 一区二区三区免费毛片| 日韩一区二区视频免费看| 午夜福利在线在线| 亚洲高清免费不卡视频| 成人一区二区视频在线观看| av国产免费在线观看| 精品久久久久久久末码| 日韩一区二区视频免费看| 日本五十路高清| а√天堂www在线а√下载| 内地一区二区视频在线| 悠悠久久av| 最近视频中文字幕2019在线8| 99久国产av精品| 成人亚洲欧美一区二区av| 97超视频在线观看视频| 国产69精品久久久久777片| 国产精品乱码一区二三区的特点| 久久久久久久久大av| 国产精品嫩草影院av在线观看| 免费一级毛片在线播放高清视频| 精品久久久久久久久亚洲| 国产精品久久视频播放| 看黄色毛片网站| 国产精品无大码| 男女视频在线观看网站免费| av专区在线播放| 2021天堂中文幕一二区在线观| 欧美激情在线99| 十八禁国产超污无遮挡网站| 亚洲乱码一区二区免费版| 亚洲aⅴ乱码一区二区在线播放| 亚洲天堂国产精品一区在线| 国产爱豆传媒在线观看| 久久韩国三级中文字幕| 看片在线看免费视频| 神马国产精品三级电影在线观看| 一个人免费在线观看电影| 亚洲第一电影网av| 最好的美女福利视频网| 久久久久九九精品影院| 天堂av国产一区二区熟女人妻| 亚洲国产精品sss在线观看| 国产精品精品国产色婷婷| 亚洲高清免费不卡视频| 一a级毛片在线观看| 国国产精品蜜臀av免费| 在线免费十八禁| 国产又黄又爽又无遮挡在线| 国产在线男女| 九九在线视频观看精品| 大型黄色视频在线免费观看| avwww免费| 久久精品综合一区二区三区| 国产男人的电影天堂91| 俄罗斯特黄特色一大片| 日韩中字成人| 久久久午夜欧美精品| 日日啪夜夜撸| 免费电影在线观看免费观看| a级毛片a级免费在线| 免费搜索国产男女视频| 在线免费观看不下载黄p国产| 天堂影院成人在线观看| 国产精品亚洲美女久久久| 永久网站在线| 国产精品无大码| 亚洲综合色惰| 我的女老师完整版在线观看| 国产精品,欧美在线| 亚洲熟妇中文字幕五十中出| 观看美女的网站| 国产精品美女特级片免费视频播放器| 亚洲av电影不卡..在线观看| 99精品在免费线老司机午夜| 日韩,欧美,国产一区二区三区 | 男女之事视频高清在线观看| 男女视频在线观看网站免费| 久久午夜福利片| 一级毛片我不卡| 国产片特级美女逼逼视频| 亚洲av不卡在线观看| 久久久久久久久中文| 黄色一级大片看看| 亚洲成人久久性| 香蕉av资源在线| 日韩成人av中文字幕在线观看 | 不卡一级毛片| 欧美又色又爽又黄视频| 午夜福利18| 看非洲黑人一级黄片| 国产精品,欧美在线| 校园人妻丝袜中文字幕| 1024手机看黄色片| 精品日产1卡2卡| 真人做人爱边吃奶动态| 亚洲经典国产精华液单| 亚洲乱码一区二区免费版| 国产欧美日韩精品亚洲av| 午夜福利在线观看吧| 1024手机看黄色片| 久久精品国产亚洲网站| 如何舔出高潮| 国产aⅴ精品一区二区三区波| av在线亚洲专区| 99久久久亚洲精品蜜臀av| 亚洲av免费在线观看| 欧美xxxx黑人xx丫x性爽| 晚上一个人看的免费电影| 国产伦精品一区二区三区四那| 国语自产精品视频在线第100页| 亚洲自拍偷在线| 国产精品国产三级国产av玫瑰| 变态另类丝袜制服| 淫秽高清视频在线观看| 欧美日韩乱码在线| 亚洲人与动物交配视频| 99热这里只有是精品50| 国产精品99久久久久久久久| 久久久午夜欧美精品| 欧美性猛交黑人性爽| av女优亚洲男人天堂| 99热这里只有是精品50| 一个人看的www免费观看视频| 亚洲成人精品中文字幕电影| 一区福利在线观看| 国产精品人妻久久久影院| 丰满乱子伦码专区| 性色avwww在线观看| 校园春色视频在线观看| 久久久欧美国产精品| 寂寞人妻少妇视频99o| 一进一出好大好爽视频| 亚洲国产精品成人综合色| 国产精品一二三区在线看| 成人美女网站在线观看视频| 人人妻人人澡人人爽人人夜夜 | 色尼玛亚洲综合影院| 天天一区二区日本电影三级| 亚洲精品亚洲一区二区| 长腿黑丝高跟| 欧美xxxx黑人xx丫x性爽| 国产v大片淫在线免费观看| 18+在线观看网站| 精品一区二区三区人妻视频| 国产蜜桃级精品一区二区三区| 九色成人免费人妻av| 日本免费a在线| 一本久久中文字幕| 久久久久久久久久黄片| 舔av片在线| 美女 人体艺术 gogo| 成人欧美大片| 日韩欧美 国产精品| 久久中文看片网| 在线观看午夜福利视频| 亚洲av成人精品一区久久| 精品一区二区三区人妻视频| 国产国拍精品亚洲av在线观看| 简卡轻食公司| 亚洲无线观看免费| 3wmmmm亚洲av在线观看| 国产成人a∨麻豆精品| 久久久a久久爽久久v久久| 亚洲一区二区三区色噜噜| 变态另类丝袜制服| 中文字幕精品亚洲无线码一区| 午夜免费男女啪啪视频观看 | 成人性生交大片免费视频hd| 午夜日韩欧美国产| 最好的美女福利视频网| 欧美zozozo另类| 日韩欧美国产在线观看| 三级国产精品欧美在线观看| 综合色av麻豆| 欧美中文日本在线观看视频| 欧美不卡视频在线免费观看| 我的女老师完整版在线观看| 国产黄色视频一区二区在线观看 | 美女内射精品一级片tv| 最近最新中文字幕大全电影3| 寂寞人妻少妇视频99o| 男女边吃奶边做爰视频| 久久久久性生活片| 欧美+日韩+精品| 热99在线观看视频| 性色avwww在线观看| 人妻丰满熟妇av一区二区三区| 麻豆国产97在线/欧美| 性色avwww在线观看| 精品久久久久久久人妻蜜臀av| 日本五十路高清| av卡一久久| 中文字幕久久专区| 2021天堂中文幕一二区在线观| 69av精品久久久久久| 国产在视频线在精品| 91久久精品电影网| 亚洲性夜色夜夜综合| 如何舔出高潮| 国内精品美女久久久久久| 最近视频中文字幕2019在线8| 久久草成人影院| 三级经典国产精品| 亚洲国产欧洲综合997久久,| 少妇人妻精品综合一区二区 | 国产av一区在线观看免费| 成人av一区二区三区在线看| or卡值多少钱| 校园春色视频在线观看| 国产免费一级a男人的天堂| a级毛片a级免费在线| 欧洲精品卡2卡3卡4卡5卡区| 1000部很黄的大片| 舔av片在线| 午夜久久久久精精品| 嫩草影视91久久| 欧美zozozo另类| av专区在线播放| 美女免费视频网站| 一边摸一边抽搐一进一小说| 青春草视频在线免费观看| 国产真实伦视频高清在线观看| 欧美激情在线99| 联通29元200g的流量卡| 国产三级在线视频| 男人的好看免费观看在线视频| 看黄色毛片网站| 最新中文字幕久久久久| 亚洲四区av| 人人妻人人看人人澡| 欧美日韩精品成人综合77777| 成人性生交大片免费视频hd| 国产毛片a区久久久久| 亚洲av成人av| 成年女人毛片免费观看观看9| 全区人妻精品视频| 精品久久国产蜜桃| 中出人妻视频一区二区| 欧美人与善性xxx| 国产精品人妻久久久影院| 男人狂女人下面高潮的视频| 身体一侧抽搐| 久久6这里有精品| 国产精品野战在线观看| 男女做爰动态图高潮gif福利片| 能在线免费观看的黄片| 亚洲av五月六月丁香网| 亚洲av成人av| 免费黄网站久久成人精品| 国产高清视频在线播放一区| 亚洲欧美日韩卡通动漫| 真人做人爱边吃奶动态| 亚洲综合色惰| 免费看av在线观看网站| 国产精品一区二区性色av| 欧美在线一区亚洲| 久久久久久久久中文| 我的女老师完整版在线观看| 干丝袜人妻中文字幕| 永久网站在线| 欧美在线一区亚洲| 日韩欧美在线乱码| 免费不卡的大黄色大毛片视频在线观看 | 国产亚洲91精品色在线| 亚洲人成网站在线播| 久久久午夜欧美精品| 男人的好看免费观看在线视频| 亚洲,欧美,日韩| 国产成人福利小说| 亚洲久久久久久中文字幕| 1000部很黄的大片| 一级毛片久久久久久久久女| 欧美另类亚洲清纯唯美| 国产极品精品免费视频能看的| 嫩草影院精品99| 国产综合懂色| 五月伊人婷婷丁香| 亚洲人与动物交配视频| 你懂的网址亚洲精品在线观看 | 在线天堂最新版资源| 国产精品三级大全| 级片在线观看| 亚洲乱码一区二区免费版| 最近在线观看免费完整版| 日本-黄色视频高清免费观看| 少妇的逼好多水| 国内久久婷婷六月综合欲色啪| 久久精品国产亚洲av香蕉五月| 久久精品国产鲁丝片午夜精品| 久久精品国产清高在天天线| 国产成人a区在线观看| 国产精品人妻久久久影院| 99热精品在线国产| or卡值多少钱| 国内精品一区二区在线观看| 成人无遮挡网站| 亚洲电影在线观看av| 欧美色视频一区免费| 国产高清视频在线播放一区| 国产精品一二三区在线看| 亚洲精品日韩av片在线观看| 色尼玛亚洲综合影院| 美女黄网站色视频| 亚洲av中文字字幕乱码综合| 少妇裸体淫交视频免费看高清| 成人高潮视频无遮挡免费网站| 男女视频在线观看网站免费| 国产极品精品免费视频能看的| 国产真实乱freesex| 男人舔奶头视频| 免费在线观看成人毛片| 看非洲黑人一级黄片| 亚洲第一电影网av| 国产av不卡久久| 午夜日韩欧美国产| 成人午夜高清在线视频| 国产精品亚洲一级av第二区| 亚洲国产色片| 色在线成人网| 国内久久婷婷六月综合欲色啪| 久久久久久久久久黄片| 久久久久久久久中文| 少妇裸体淫交视频免费看高清| 黄色一级大片看看| 国产精华一区二区三区| 亚洲精品日韩在线中文字幕 | 看免费成人av毛片| 久久久久九九精品影院| 亚洲精品亚洲一区二区| 特大巨黑吊av在线直播| 丰满人妻一区二区三区视频av| 啦啦啦韩国在线观看视频| 嫩草影院精品99| 亚洲国产精品合色在线| 国产精华一区二区三区| 欧美又色又爽又黄视频| 如何舔出高潮| 久久久久久久久久黄片| 午夜视频国产福利| а√天堂www在线а√下载| 网址你懂的国产日韩在线| 日韩欧美精品免费久久| 12—13女人毛片做爰片一| 热99re8久久精品国产| 免费看a级黄色片| 亚洲久久久久久中文字幕| 日韩三级伦理在线观看| 免费av观看视频| 嫩草影院新地址| 欧美日韩综合久久久久久| 最近在线观看免费完整版| 日产精品乱码卡一卡2卡三| 中文字幕免费在线视频6| 97超视频在线观看视频| 亚洲欧美成人精品一区二区| 狠狠狠狠99中文字幕| 欧美日韩综合久久久久久| 在线免费十八禁| 日韩欧美精品v在线| 欧美高清性xxxxhd video| 丰满人妻一区二区三区视频av| 天堂av国产一区二区熟女人妻| 精品一区二区免费观看| 午夜久久久久精精品| 国产一区二区激情短视频| 内射极品少妇av片p| videossex国产| 老熟妇仑乱视频hdxx| 成人二区视频| 亚洲一级一片aⅴ在线观看| 波多野结衣巨乳人妻| 搡老妇女老女人老熟妇| 免费观看在线日韩| 色视频www国产| 久久婷婷人人爽人人干人人爱| 精品国内亚洲2022精品成人| 亚洲av第一区精品v没综合| 你懂的网址亚洲精品在线观看 | 精品一区二区三区视频在线| 深夜a级毛片| 国产aⅴ精品一区二区三区波| 欧美一区二区国产精品久久精品| 少妇高潮的动态图| 别揉我奶头~嗯~啊~动态视频| 深夜a级毛片| 国产综合懂色| 成人漫画全彩无遮挡| 亚洲专区国产一区二区| 在线免费观看不下载黄p国产| 最新中文字幕久久久久| 美女被艹到高潮喷水动态| 日日干狠狠操夜夜爽| 国产极品精品免费视频能看的| 桃色一区二区三区在线观看| 高清日韩中文字幕在线| 国产精品不卡视频一区二区| 国产精品久久视频播放| 97碰自拍视频| 亚洲第一电影网av| or卡值多少钱| h日本视频在线播放| 亚洲国产色片| 男女下面进入的视频免费午夜| 亚洲三级黄色毛片| 国产黄色小视频在线观看| 久久精品91蜜桃| 国产精品一二三区在线看| 91麻豆精品激情在线观看国产| 国产激情偷乱视频一区二区| 毛片一级片免费看久久久久| 日韩av不卡免费在线播放| 蜜桃久久精品国产亚洲av| 久久人人爽人人爽人人片va| 不卡视频在线观看欧美| 欧美色视频一区免费| 深夜a级毛片| 精品一区二区三区av网在线观看| 免费看美女性在线毛片视频| 国产亚洲精品久久久久久毛片| 在线免费观看不下载黄p国产| 男女下面进入的视频免费午夜| 亚洲丝袜综合中文字幕| 岛国在线免费视频观看| 成人毛片a级毛片在线播放| 午夜精品国产一区二区电影 | 亚洲激情五月婷婷啪啪| 桃色一区二区三区在线观看| 国产精品无大码| 老师上课跳d突然被开到最大视频| 国产精品三级大全| 少妇高潮的动态图| 最近的中文字幕免费完整| 久久久成人免费电影| 亚洲欧美成人精品一区二区| 一本久久中文字幕| 成人三级黄色视频| 欧美最黄视频在线播放免费| 欧美激情在线99| 久久热精品热| 国产乱人偷精品视频| 日韩强制内射视频| 全区人妻精品视频| 久久精品国产99精品国产亚洲性色| 亚洲av美国av| 天天一区二区日本电影三级| 国产真实乱freesex| 日韩中字成人| 波多野结衣巨乳人妻| 97热精品久久久久久| 99久久久亚洲精品蜜臀av| 欧美高清性xxxxhd video| 久久久久久大精品| 高清午夜精品一区二区三区 | 欧美激情久久久久久爽电影| 国产美女午夜福利| 99热网站在线观看| av在线亚洲专区| 一本一本综合久久| 18禁黄网站禁片免费观看直播| 免费人成在线观看视频色| 男女视频在线观看网站免费| 日产精品乱码卡一卡2卡三| 中文字幕人妻熟人妻熟丝袜美| 免费大片18禁| 99久久久亚洲精品蜜臀av| 欧洲精品卡2卡3卡4卡5卡区| 我要看日韩黄色一级片| 亚洲中文字幕一区二区三区有码在线看| 中出人妻视频一区二区| 国产一区二区亚洲精品在线观看| 亚洲专区国产一区二区| 亚洲一区高清亚洲精品| 九色成人免费人妻av| av福利片在线观看| 成人鲁丝片一二三区免费| 少妇的逼好多水| 波野结衣二区三区在线| 国产日本99.免费观看| 亚洲18禁久久av| 伊人久久精品亚洲午夜| 国产麻豆成人av免费视频| 精品一区二区免费观看| 熟妇人妻久久中文字幕3abv| 麻豆国产av国片精品| 久久99热这里只有精品18| 国产精品人妻久久久久久| 亚洲最大成人手机在线| avwww免费| 亚洲无线观看免费| 亚洲国产色片| 亚洲av.av天堂| 日韩,欧美,国产一区二区三区 | 国内精品宾馆在线| av在线观看视频网站免费| 中文字幕久久专区| 亚洲欧美日韩卡通动漫| 天天躁夜夜躁狠狠久久av| 如何舔出高潮| or卡值多少钱| 三级男女做爰猛烈吃奶摸视频| 国产单亲对白刺激| 成人特级av手机在线观看| 精品一区二区三区视频在线观看免费| 国产精品久久电影中文字幕|