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

    爆轟加載下銀氣溶膠源項(xiàng)實(shí)驗(yàn)研究*

    2016-04-18 03:03:30劉文杰馬慶鵬王鵬來
    爆炸與沖擊 2016年6期
    關(guān)鍵詞:采樣器氣溶膠炸藥

    劉文杰,馬慶鵬,王鵬來

    (中國工程物理研究院流體物理研究所,四川 綿陽 621999)

    爆轟加載下銀氣溶膠源項(xiàng)實(shí)驗(yàn)研究*

    劉文杰,馬慶鵬,王鵬來

    (中國工程物理研究院流體物理研究所,四川 綿陽 621999)

    為研究钚氣溶膠生成規(guī)律及源項(xiàng)分布,提出了一種密封環(huán)境下的金屬氣溶膠源項(xiàng)實(shí)驗(yàn)方法,利用爆炸容器及氣溶膠采樣器等實(shí)驗(yàn)設(shè)備,開展了钚替代材料銀的氣溶膠源項(xiàng)生成實(shí)驗(yàn),分析了空氣動(dòng)力學(xué)直徑小于10 μm的銀氣溶膠源項(xiàng)分布特性,并與外場(chǎng)擴(kuò)散試驗(yàn)中钚氣溶膠數(shù)據(jù)進(jìn)行了對(duì)比。進(jìn)一步采用不同炸藥及加載裝置,對(duì)相同銀片樣品進(jìn)行加載實(shí)驗(yàn),研究了加載峰值壓力對(duì)銀氣溶膠源項(xiàng)分布的影響。研究結(jié)果表明,金屬銀可作為替代材料研究爆轟加載條件下钚氣溶膠的源項(xiàng)參數(shù),在一定實(shí)驗(yàn)條件下,二者歸一化積累質(zhì)量分布具有較高的一致性;爆轟加載產(chǎn)生的銀氣溶膠中,較小粒徑氣溶膠質(zhì)量含量較高。氣溶膠總量與加載峰值壓力呈二次函數(shù)分布,當(dāng)峰值壓力超過某個(gè)臨界值時(shí),氣溶膠總量將達(dá)到極值;此外,通過對(duì)同發(fā)次實(shí)驗(yàn)的3次不同采樣數(shù)據(jù)進(jìn)行分析,發(fā)現(xiàn)了較明顯的氣溶膠凝并和沉降現(xiàn)象。

    爆炸力學(xué);粒徑分布;歸一化積累質(zhì)量;銀氣溶膠;源項(xiàng);钚氣溶膠

    钚-239(下稱钚)材料一般通過反應(yīng)堆中產(chǎn)生的慢中子轟擊鈾-238人工產(chǎn)生,是武器級(jí)核裝置的重要裂變材料和核電站的初始加濃原料。作為高活性和高毒性的放射性核素,钚材料具有強(qiáng)烈地誘發(fā)骨腫瘤的生物毒性,且半衰期極長(zhǎng),約為24 400年。钚材料的化學(xué)性質(zhì)極為活潑,常見的非合金態(tài)α位相是一種室溫下堅(jiān)硬的金屬,但僅在小于122 ℃時(shí)穩(wěn)定。當(dāng)含钚材料的核裝置處于高能炸藥化學(xué)爆炸和燃料燃燒等事故環(huán)境中時(shí),钚材料將發(fā)生氧化反應(yīng)并氣溶膠化[1-2]。而钚氣溶膠擴(kuò)散至環(huán)境后極易因呼吸進(jìn)入生物肺部和血液循環(huán)系統(tǒng),使生物機(jī)體在存活期間遭受持續(xù)地內(nèi)照射[3-4]。R.H.Condit[5]將可能對(duì)人體產(chǎn)生吸入危害的钚氣溶膠顆粒的空氣動(dòng)力學(xué)直徑上限定為10 μm,急性致死劑量為4.3 mg。因此,掌握爆轟加載下钚氣溶膠的源項(xiàng)參數(shù)(產(chǎn)生氣溶膠的總質(zhì)量及粒徑分布),對(duì)意外核事故的預(yù)防及應(yīng)急響應(yīng)措施的制定具有重大的意義。

    對(duì)于核裝置化爆條件下钚氣溶膠的研究,僅有在內(nèi)華達(dá)核試驗(yàn)場(chǎng)進(jìn)行的大型外場(chǎng)安全試驗(yàn)ORC[4,6]。迄今為止,僅有一次系統(tǒng)完備地測(cè)量炸藥化爆條件下钚氣溶膠的擴(kuò)散試驗(yàn),ORC為幾十年間針對(duì)化爆條件下钚氣溶膠擴(kuò)散的數(shù)值模擬提供了重要且唯一的源項(xiàng)及擴(kuò)散分布數(shù)據(jù)[7-8]。但是,考慮對(duì)生物和環(huán)境的危害,之后對(duì)于钚氣溶膠的源項(xiàng)及擴(kuò)散研究轉(zhuǎn)入核事故、核試驗(yàn)場(chǎng)地的钚氣溶膠再懸浮[9-10]、沉積評(píng)估[11-12]以及實(shí)驗(yàn)室實(shí)驗(yàn)、數(shù)值模擬。M.J.Sagartz[13]通過PBX-9404炸藥加載銀球殼實(shí)驗(yàn),模擬核裝置在事故條件下钚氣溶膠源項(xiàng)生成規(guī)律,研究了不同炸藥構(gòu)型條件下銀樣品氣溶膠化的轉(zhuǎn)化率,分析了炸藥爆能對(duì)銀氣溶膠轉(zhuǎn)化率的影響規(guī)律,但未說明使用金屬銀作為钚替代材料的依據(jù)。伍懷龍等[14]對(duì)钚氣溶膠環(huán)境中放射性氣體的遷移及擴(kuò)散過程開展了實(shí)驗(yàn)研究,建立了描述放射性氣體遷移過程的物理模型,但對(duì)钚氣溶膠的源項(xiàng)特征未做明確描述。在數(shù)值模擬方面,相關(guān)研究機(jī)構(gòu)相繼開發(fā)了多種針對(duì)钚氣溶膠擴(kuò)散研究的數(shù)值計(jì)算模型,并對(duì)核事故擴(kuò)散場(chǎng)景進(jìn)行了數(shù)值計(jì)算。主要有Hotspot模型[15]、DIFOUT模型[16]、MACCS2模型[17]和ERAD模型[18],前3種是高斯煙云模型,后1種是采用蒙特卡洛方法的三維數(shù)值模型。

    在模擬核事故條件下钚氣溶膠的擴(kuò)散機(jī)制時(shí),所有模型的輸入?yún)?shù)均依賴于事故發(fā)生的源項(xiàng)數(shù)據(jù)。然而出于安全性等方面的考慮,直接開展钚材料化爆實(shí)驗(yàn)研究缺乏現(xiàn)實(shí)可行性。鑒于銀是一種低熔點(diǎn) (961.78 ℃)、高沸點(diǎn) (2 163 ℃)的金屬,除密度相對(duì)較低外,其基本物理化學(xué)性質(zhì)與钚具有一定的相似性,且不存在放射毒性等問題。因此,本文中,擬采用銀作為钚的替代材料,開展炸藥化爆條件下銀氣溶膠的源項(xiàng)實(shí)驗(yàn),研究粒徑10 μm以下 (PM10) 的銀氣溶膠源項(xiàng)分布特性,探索使用銀模擬钚材料氣溶膠源項(xiàng)分布的可行性,并進(jìn)一步利用多種爆轟模型及加載方式,研究加載峰值壓力對(duì)銀氣溶膠源項(xiàng)參數(shù)的影響特性。

    1 實(shí)驗(yàn)方法

    圖1 爆轟加載金屬銀氣溶膠實(shí)驗(yàn)裝置示意圖Fig.1 Illustration of experimental setup

    圖1為爆轟加載金屬銀的實(shí)驗(yàn)裝置示意圖,裝置由高密封爆炸容器、氣溶膠采集系統(tǒng)(多級(jí)沖擊采樣器、壓力傳感器、熱電偶及質(zhì)量流量計(jì))和抽氣系統(tǒng)(真空泵、角閥及波紋管道)組成。其中,爆炸容器為臥式柱形,兩端為球頭,設(shè)計(jì)TNT當(dāng)量為2.5 kg,容器采用高密封設(shè)計(jì),600 h總泄漏份額不大于10-8。氣溶膠采集系統(tǒng)的核心裝置是安德森八級(jí)撞擊式氣溶膠采樣器TE-20-800,可對(duì)0 ~ 10 μm粒徑范圍內(nèi)的氣溶膠進(jìn)行分段采集,圖2(a)為該采樣器的實(shí)物圖。該采樣器具有一個(gè)10 μm粒徑的前置切割級(jí),用來切割、分離所有粒徑大于等于10 μm的顆粒,如圖2(b)所示;中間八級(jí)切割粒徑分別為9.0、5.8、4.7、3.3、2.1、1.1、0.7和0.4 μm,對(duì)應(yīng)沖擊板為0~7,分別采集9~10 μm至0.4~0.7 μm粒徑的氣溶膠顆粒;后置級(jí)用來收集所有粒徑小于0.4 μm的氣溶膠顆粒。采樣器利用慣性原理對(duì)氣溶膠進(jìn)行粒徑切割,其切割原理如圖2(c)所示。當(dāng)含氣溶膠氣流由各級(jí)加速噴嘴噴出時(shí),受到?jīng)_擊板的阻擋,粒徑小、質(zhì)量輕的氣溶膠將隨氣流改變流動(dòng)方向流入下一級(jí),粒徑大、質(zhì)量重的氣溶膠將打在沖擊板上進(jìn)而被濾膜收集。利用該原理,每級(jí)濾膜收集到的氣溶膠粒徑分布在該級(jí)和上一級(jí)截?cái)嘀睆椒秶鷥?nèi),從而實(shí)現(xiàn)了氣溶膠多級(jí)采樣。選用對(duì)金屬顆粒物吸附性能最優(yōu)的高密度石英纖維濾膜作為氣溶膠樣品采集基底。通過壓力傳感器和熱電偶監(jiān)測(cè)容器環(huán)境壓力及溫度;通過真空抽氣系統(tǒng)的可調(diào)角閥控制抽速,并通過質(zhì)量流量計(jì)監(jiān)測(cè)采樣流量,將它限制為(28.3±1.5) L/min,以保證采樣器的粒徑切割精度。

    圖2 安德森八級(jí)撞擊式氣溶膠采樣器Fig.2 Anderson 8 stage non-viable cascade impactor

    圖3 平面對(duì)碰爆轟加載裝置示意圖Fig.3 Sketch of detonative device

    為了獲得較高的加載壓力,主要采用平面對(duì)碰加載模型(見圖3)對(duì)銀片進(jìn)行爆轟加載,該構(gòu)型與ORC-DT钚氣溶膠外場(chǎng)擴(kuò)散試驗(yàn)加載構(gòu)型[18]類似。單片炸藥尺寸為?100 mm×20 mm,銀片放置于兩片炸藥之間;起爆裝置使用BL-21雷管,傳爆藥為泰安?32 mm×11 mm;裝置外用厚度5 mm的有機(jī)玻璃外殼進(jìn)行約束;實(shí)驗(yàn)銀樣品均為高純銀片(含銀99.99%),質(zhì)量為45 g。

    實(shí)驗(yàn)狀態(tài)參數(shù)見表1,共開展7發(fā)次爆轟加載條件下的銀氣溶膠實(shí)驗(yàn)。其中前3發(fā)次為標(biāo)定實(shí)驗(yàn),用于評(píng)估采樣器的氣溶膠采樣誤差;后4發(fā)次為變加載壓力的銀氣溶膠源相實(shí)驗(yàn),是為了獲得多種加載壓力條件下銀氣溶膠生成總量及粒徑分布。加載分為單面和對(duì)碰兩種方式,單面加載實(shí)驗(yàn)裝置仍如圖3所示,但拆除了底部起爆雷管2和傳爆藥2。p為各加載條件下作用在銀樣品上的峰值壓力的理論計(jì)算結(jié)果。其中,平面對(duì)碰及單面加載兩種加載方式下作用在銀片上的峰值壓力以炸藥的CJ爆轟參數(shù)為輸入?yún)?shù),根據(jù)Rankine-Hugoniot曲線得到。

    表1 實(shí)驗(yàn)狀態(tài)參數(shù)Table 1 Experimental details

    實(shí)驗(yàn)前,先清洗容器內(nèi)壁,并盡可能減少容器內(nèi)可能影響氣溶膠樣品分析的實(shí)驗(yàn)輔助材料,如加載裝置采用棉繩懸掛而取消木支架,雷管套采用有機(jī)玻璃而不用鋁材料。

    2 結(jié)果及分析

    2.1 采樣數(shù)據(jù)匯總及誤差

    炸藥爆轟零時(shí)后,立即打開容器密封閥和真空抽氣系統(tǒng),利用采樣器開始對(duì)氣相爆轟產(chǎn)物進(jìn)行分級(jí)采集。每次采樣時(shí)間均為10 min,單次累積抽氣量為283 L,約為待測(cè)氣溶膠體系總量(容器總體積)的3.877 %。采樣結(jié)束后,回收采樣濾膜對(duì)沉積在各級(jí)濾膜上的氣溶膠顆粒物進(jìn)行物化分析:利用粒度分析儀或電化學(xué)消解后能譜分析,給出容器內(nèi)氣相銀氣溶膠粒徑分布范圍;利用原子吸收光譜方法,給出氣相銀氣溶膠采集總量。7發(fā)次實(shí)驗(yàn)的銀氣溶膠各粒徑的采樣分析結(jié)果如表2所示,其中實(shí)驗(yàn)3進(jìn)行了3次間隔采樣。

    對(duì)氣溶膠源項(xiàng)生成的分析主要從兩個(gè)方面開展:一是各粒徑氣溶膠質(zhì)量分布,二是氣溶膠的歸一化積累質(zhì)量分布(小于某一粒徑的氣溶膠質(zhì)量占采樣總質(zhì)量的百分比)。

    圖4為實(shí)驗(yàn)1~3標(biāo)定得到的銀氣溶膠歸一化積累質(zhì)量分布。對(duì)比采樣1、2可知,在3.3~10 μm粒徑范圍內(nèi),銀氣溶膠歸一化積累質(zhì)量分布基本吻合;在小于3.3 μm粒徑范圍內(nèi),因小粒徑顆粒懸浮運(yùn)動(dòng)具有一定隨機(jī)性且易受氣流擾動(dòng)的影響,3次采樣結(jié)果稍有差異,但變化趨勢(shì)一致。綜上,3次采樣銀氣溶膠歸一化積累質(zhì)量分布趨勢(shì)一致,僅在數(shù)值上稍有差異。因此,可以認(rèn)為,采用密封容器開展爆轟加載條件下銀氣溶膠源項(xiàng)的研究方法是可行的,具有較高的可重復(fù)性,其數(shù)據(jù)是合理、有效的。

    表2 銀氣溶膠采樣數(shù)據(jù)Table 2 Data of collected samples of silver aerosol

    圖4 爆轟加載下銀氣溶膠源項(xiàng)生成數(shù)據(jù)標(biāo)定Fig.4 Calibration of silver aerosol source-term after explosion

    表3 電化學(xué)消解銀標(biāo)樣試液濃度Table 3 Calibration error of samples

    2.2 銀氣溶膠源項(xiàng)分布及沉降、凝并

    圖5 銀氣溶膠源項(xiàng)分布Fig.5 Distribution of silver aerosol source-term

    圖5為7發(fā)次實(shí)驗(yàn)爆轟零時(shí)后采樣銀氣溶膠粒徑質(zhì)量分布,這7發(fā)次實(shí)驗(yàn)銀氣溶膠的粒徑質(zhì)量分布具有一定的相似性。各次采樣結(jié)果表明,爆轟加載產(chǎn)生的銀氣溶膠含量峰值均出現(xiàn)在1.1~3.3 μm粒徑范圍內(nèi),在不同的加載方式下,該峰值在1.1~2.1 μm和2.1~3.3 μm兩個(gè)粒徑范圍內(nèi)變化。爆轟零時(shí)后,采樣1、2、3a、4、5、6、7中,0.7~3.3 μm銀氣溶膠百分含量分別為66.8%、76.9%、91.0%、93.4%、79.8%、82.2%、60.1%。在整個(gè)采樣區(qū)間內(nèi),粒徑在0.7~3.3 μm范圍內(nèi)的氣溶膠含量占總量比例均高于60%。這表明,在對(duì)金屬銀進(jìn)行爆轟加載時(shí),產(chǎn)生較小粒徑 (0.7~3.3 μm) 氣溶膠含量較高,在RHT-901高猛炸藥加載下,這個(gè)粒徑范圍內(nèi)的氣溶膠含量甚至可達(dá)90%以上。

    圖6為實(shí)驗(yàn)3(RHT-901平面對(duì)碰加載)中3次間隔采樣銀氣溶膠歸一化積累質(zhì)量分布及粒徑質(zhì)量分布。由圖6(a)可知,爆轟加載零后2 h采樣3c的銀氣溶膠歸一化積累質(zhì)量分布,與ORC-DT外場(chǎng)擴(kuò)散試驗(yàn)钚氣溶膠數(shù)據(jù)[19-20]具有較高的一致性。這表明,采用銀作為替代材料對(duì)爆轟加載條件下钚氣溶膠源項(xiàng)的模擬是完全可行的。此外,對(duì)比3次采樣數(shù)據(jù)(見圖6(a))可知,三者在歸一化積累質(zhì)量分布方面有較大的差別,采樣時(shí)間越早,小粒徑氣溶膠含量越高。為此,對(duì)3次采樣得到的粒徑質(zhì)量分布(見圖6(b))做進(jìn)一步地分析。在0~2.1 μm粒徑范圍內(nèi),隨著采樣時(shí)間的推移,該范圍內(nèi)氣溶膠質(zhì)量呈現(xiàn)明顯的下降趨勢(shì);而在2.1~10 μm粒徑范圍內(nèi),第2次采樣得到的氣溶膠質(zhì)量比第1次有明顯下降,第3次采樣得到的氣溶膠質(zhì)量比第2次有明顯上升;在5.8~10 μm粒徑范圍內(nèi),第3次的結(jié)果甚至高于第1次采樣結(jié)果。這表明,隨著時(shí)間的推移,容器內(nèi)的銀氣溶膠發(fā)生了較顯著的凝并及沉降現(xiàn)象。

    圖6 實(shí)驗(yàn)3氣溶膠源項(xiàng)分布Fig.6 Source-term of silver aerosol of experiment 3

    綜上所述,在采用的實(shí)驗(yàn)密封環(huán)境實(shí)驗(yàn)條件下,采用RHT-901炸藥對(duì)銀金屬進(jìn)行平面對(duì)碰爆轟加載,產(chǎn)生的銀氣溶膠經(jīng)過2 h穩(wěn)定后,較接近開放環(huán)境化爆條件下钚氣溶膠源項(xiàng)分布。同時(shí),爆轟零時(shí)后產(chǎn)生的0.7~3.3 μm粒徑范圍內(nèi)銀氣溶膠占PM10銀氣溶膠比例較高,且隨著時(shí)間推移,容器內(nèi)小粒徑氣溶膠發(fā)生較明顯的凝并和沉降現(xiàn)象。

    2.3 加載峰值壓力的影響

    圖7 不同峰值壓力下采樣中銀氣溶膠總量Fig.7 Amount of silver aerosol under different pressures

    參照表1~2,得到加載峰值壓力與采樣中銀氣溶膠總量的對(duì)應(yīng)關(guān)系,并對(duì)數(shù)據(jù)點(diǎn)進(jìn)行二階多項(xiàng)式擬合,如圖7所示,其中RHT-901平面對(duì)碰加載實(shí)驗(yàn)1~3的氣溶膠總量取平均值。

    由圖7可知,隨著峰值壓力增大,銀氣溶膠總量呈先增大后減小的趨勢(shì),表明氣溶膠的濃度不是隨峰值壓力而漸增的。通過對(duì)比數(shù)值擬合曲線與實(shí)驗(yàn)采樣結(jié)果可知,二者一致性很高,這表明銀氣溶膠總量與加載峰值壓力存在典型的二次函數(shù)關(guān)系。依據(jù)擬合曲線可知,在45.85~83.11 GPa的峰值壓力區(qū)間內(nèi),存在一個(gè)銀氣溶膠產(chǎn)物總量極值,通過計(jì)算,得出該極值對(duì)應(yīng)的壓力約為68.69 GPa。對(duì)于這個(gè)結(jié)果,文獻(xiàn)[21]中提到,在對(duì)多個(gè)核爆試驗(yàn)場(chǎng)和次臨界試驗(yàn)場(chǎng)的擾動(dòng)再懸浮钚氣溶膠濃度測(cè)量時(shí)得出,核爆試驗(yàn)場(chǎng)的氣溶膠濃度明顯低于次臨界試驗(yàn)場(chǎng)。經(jīng)對(duì)核爆試驗(yàn)場(chǎng)爆心土壤進(jìn)行篩選、電鏡分析后得到大量的微小的琉璃化钚化合物。這表明,當(dāng)加載峰值壓力超過一定值時(shí),被加載金屬材料化合物可能會(huì)發(fā)生琉璃化現(xiàn)象,抑制小粒徑氣溶膠顆粒的形成,進(jìn)而導(dǎo)致氣溶膠濃度降低。

    3 結(jié) 論

    設(shè)計(jì)了爆炸容器內(nèi)金屬材料在爆轟加載條件下的氣溶膠源項(xiàng)實(shí)驗(yàn)方案,通過對(duì)銀材料爆轟加載獲得了PM10銀氣溶膠的源項(xiàng)參數(shù),實(shí)現(xiàn)了對(duì)ORC-DT試驗(yàn)中钚氣溶膠源項(xiàng)分布的模擬。具體獲得了以下結(jié)論。

    (1) 實(shí)驗(yàn)結(jié)果表明,在爆炸容器約束空間內(nèi),爆轟加載下產(chǎn)生的銀氣溶膠與钚氣溶膠的歸一化積累質(zhì)量分布具有較高的一致性。在一定條件下,可以使用銀作為替代材料研究钚氣溶膠的源項(xiàng)分布,為后續(xù)氣溶膠擴(kuò)散的數(shù)值計(jì)算提供輸入?yún)?shù)。

    (2) 在對(duì)金屬銀進(jìn)行爆轟加載時(shí),較小粒徑 (0.7~3.3 μm) 氣溶膠含量較高,在RHT-901高猛炸藥加載下,這個(gè)粒徑范圍內(nèi)的氣溶膠含量甚至可達(dá)總量的90%以上。

    (3) 結(jié)果表明,爆轟零時(shí)后,銀氣溶膠發(fā)生顯著的凝并和沉降現(xiàn)象:爆轟零時(shí)后1 h內(nèi),各粒徑范圍氣溶膠含量均有下降,較小粒徑氣溶膠含量下降更明顯;零后2 h內(nèi),較小粒徑氣溶膠含量進(jìn)一步減少,但較大粒徑氣溶膠含量不降反升。

    (4) 得到加載峰值壓力和生成PM10銀氣溶膠總量的變化規(guī)律,銀氣溶膠總量與加載峰值壓力呈二次函數(shù)分布,該現(xiàn)象發(fā)生的原因可能在于高加載壓力時(shí)銀材料發(fā)生了琉璃化。進(jìn)一步研究,在該實(shí)驗(yàn)裝置上,加載峰值壓力約為68.69 GPa時(shí),可能獲得的PM10銀氣溶膠質(zhì)量最大。對(duì)該極值,擬在后續(xù)的實(shí)驗(yàn)工作中予以確認(rèn)。

    [1] Martz J C, Haschke J M. A mechanism for combustive heating and explosive dispersal of plutonium[J]. Journal of Alloys and Compounds, 1998,266(1):90-103.

    [2] Haschke J M, Martz J C. Oxidation kinetics of plutonium in air from 500 to 3 500 ℃: Application to source terms for dispersal[J]. Journal of Alloys and Compounds, 1998,266(1):81-89.

    [3] 劉文杰,胡八一,李慶忠.核事故條件下钚氣溶膠源項(xiàng)研究綜述[J].安全與環(huán)境學(xué)報(bào),2011,11(5):259-263. Liu Wenjie, Hu Bayi, Li Qingzhong. Review of plutonium aerosol source-term in nuclear accident[J]. Journal of Safety and Environment, 2011,11(5):259-263.

    [4] Stewart K, Thomas D C, Terry J L, et al. Atomic weapons research establishment[R]. AWRE REPORT NO.0-29/65, 1965.

    [5] Condit R H. Plutonium dispersal in fires: Summary of what is known[R]. UCRL-ID-114164, 1993.

    [6] Titus R W. Operation roller coaster: Project officers report-project 2.4[R]. POR-2504, 1965.

    [7] Boughton B A, DeLaurentis J M. An integral model of plume rise from high explosive detonations[C]∥Davis L R, Catton I. Buoyant plumes: Proceedings of the 24th National Heat Transfer Conference. Pittsburgh, PA, USA: ASME, 1987:27-32.

    [8] Boughton B A, DeLaurentis J M, Dunn W E. A note on a stochastic model of turbulent dispersion in the atmosphere[J]. Boundary-Layer Meteorology, 1989,48(4):443-444.

    [9] Shinn J H, Homan D N, Hofmann C B. A summary of plutonium aerosol studies: Resuspension at the Nevada Test Site[R]. UCRL-90746, 1986.

    [10] Shinn J H. Complementary Pu resuspension study at Palomares, Spain[R]. UCRL-ID-150980, 2002.

    [11] Johansen M P, Kamboj S, Kuhne W W. Whole-organism concentration ratios for plutonium in wildlife from past US nuclear research data[J]. Journal of Environmental Radioactivity, 2013,126(6):412-419.

    [12] Johansen M P, Child D P, Davis E, et al. Plutonium in wildlife and soils at the Maralinga legacy site: Persistence over decadal time scales[J]. Journal of Environmental Radioactivity, 2014,131(2):72-80.

    [13] Sagartz M J. Violent reaction source term study[R]. SAND-943252, 1995.

    [14] 伍懷龍,田東風(fēng),郝樊華,等.钚氣溶膠環(huán)境中惰性氣體氪遷移過程研究[J].物理學(xué)報(bào),2011,60(3):162-165. Wu Huailong, Tian Dongfeng, Hao Fanhua, et al. Investigation on inert gases transport process in plutonium aerosol[J]. Acta Physica Sinica, 2011,60(3):162-165.

    [15] Homann S G, Wilson D V. Hotspot health physics code: Version8.0[R]. UCRL-MA-118617, 1995.

    [16] Shreve J D. Operation roller coaster: Scientific directors report[R]. DASA-1644, 1965.

    [17] Steele C M, Wald T L, Chanin D I. Plutonium explosive dispersal modeling using the MACCS2 computer vode[R]. LA-UR-98-1901, 1998.

    [18] Boughton B A, Delaurentis J M. Description and validation of ERAD: An atmospheric dispersion model for high explosive detonations[R]. SAND92-2069-UC-702, 1992.

    [19] Beasely R R. Operation roller coaster: Technical photography[R]. POR-2519, 1965.

    [20] Church H W. Cloud rise from high-explosives detonation[R]. SC-RR-68-903, 1969.

    [21] Shinn J H. Technical basis for air pathway assessment of resuspended radioactive aerosols: LLNL experiences at seven sites around the world[R]. UCRL-JC-115045, 1993.

    (責(zé)任編輯 丁 峰)

    Experimental research of silver aerosol source-term under explosive detonation

    Liu Wenjie, Ma Qingpeng, Wang Penglai

    (InstituteofFluidPhysics,ChinaAcademyofEngineeringPhysics,Mianyang621999,Sichuan,China)

    An experimental method in the confinement vessel for silver aerosol source-term investigation of metal under high explosive detonation was designed to simulate the plutonium aerosol source-term of nuclear devices under the circumstance of high explosive detonation. Seven silver source-term experiments under variable explosive pressure were conducted in the aerosol facility and the size distribution characteristics of the silver aerosol with an aerodynamic diameter less than 10 μm were examined and compared with the data from the field test. The results show that the silver aerosol source-term under specific conditions has good agreement with the plutonium aerosol source-term obtained by the U. S. government under test conditions simulating an open-air detonation accident of nuclear device. The amount of the silver aerosolized in the experiment is determined by the peak value of high explosive pressure that acted on the silver plate and exhibited a quadratic polynomial. The silver aerosol with an aerodynamic diameter ranging from 0.7 to 3.3 μm account for the major part of the silver aerosol induced by high explosive detonation. When the peak value exceeds a certain critical value, the total amount of aerosol reaches the maximum. Settlement and combination of aerosol were observed in the comparative analysis of the size-mass distribution of the samples collected at different times after the explosion.

    mechanics of explosion; size distribution; cumulative mass; silver aerosol; source-term; plutonium aerosol

    10.11883/1001-1455(2016)06-0774-07

    2015-04-09; < class="emphasis_bold">修回日期:2015-05-25

    2015-05-25

    國家自然科學(xué)基金項(xiàng)目(21307119)

    劉文杰(1982— ),男,博士,副研究員;

    馬慶鵬,maqingpeng@caep.cn。

    O381 <國標(biāo)學(xué)科代碼:1303510 class="emphasis_bold"> 國標(biāo)學(xué)科代碼:1303510 文獻(xiàn)標(biāo)志碼:A國標(biāo)學(xué)科代碼:1303510

    A

    猜你喜歡
    采樣器氣溶膠炸藥
    “炸藥”驚魂
    議論火炸藥數(shù)字化制造
    氣溶膠傳播之謎
    粉塵采樣器檢定和校準(zhǔn)證書中不確定度區(qū)別
    氣溶膠中210Po測(cè)定的不確定度評(píng)定
    四川盆地秋季氣溶膠與云的相關(guān)分析
    PM2.5小流量采樣器數(shù)據(jù)采集與導(dǎo)出系統(tǒng)優(yōu)化研究
    基于VB的PM2.5采樣器數(shù)據(jù)處理軟件的研究
    大氣氣溶膠成核監(jiān)測(cè)
    Al粉對(duì)炸藥爆炸加速能力的影響
    一本久久精品| 美女xxoo啪啪120秒动态图| 国产精品99久久久久久久久| 97超碰精品成人国产| 黄色日韩在线| 中文欧美无线码| 精品国产乱码久久久久久小说| 久久久久精品久久久久真实原创| 亚洲内射少妇av| 中文字幕av电影在线播放| 中文乱码字字幕精品一区二区三区| 国产精品久久久久久久久免| 亚洲av不卡在线观看| tube8黄色片| 伊人久久精品亚洲午夜| 欧美区成人在线视频| 一区二区三区精品91| 婷婷色麻豆天堂久久| 久久久久久久久久久久大奶| av免费观看日本| 日韩三级伦理在线观看| 亚洲精品国产成人久久av| 精品少妇内射三级| av卡一久久| 国语对白做爰xxxⅹ性视频网站| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美一区二区三区国产| 18禁动态无遮挡网站| 精品亚洲成国产av| 久久精品熟女亚洲av麻豆精品| kizo精华| 黄色配什么色好看| 免费观看在线日韩| 国产女主播在线喷水免费视频网站| 亚洲一区二区三区欧美精品| 3wmmmm亚洲av在线观看| 欧美精品国产亚洲| a 毛片基地| 欧美精品高潮呻吟av久久| 成人亚洲欧美一区二区av| 日韩精品有码人妻一区| 啦啦啦视频在线资源免费观看| 嘟嘟电影网在线观看| 麻豆成人av视频| 亚洲真实伦在线观看| 成人特级av手机在线观看| 久久久久精品性色| 亚洲欧美成人精品一区二区| 久久国产乱子免费精品| 黄色毛片三级朝国网站 | 色网站视频免费| 亚洲欧美成人综合另类久久久| 日本爱情动作片www.在线观看| 欧美丝袜亚洲另类| 亚洲无线观看免费| 日韩中字成人| 成人无遮挡网站| 国产乱人偷精品视频| 夫妻性生交免费视频一级片| 大码成人一级视频| 午夜激情福利司机影院| 久久久久久久久久成人| 日韩不卡一区二区三区视频在线| 一本久久精品| 国产av一区二区精品久久| 91久久精品电影网| 久久午夜综合久久蜜桃| 在线亚洲精品国产二区图片欧美 | 欧美日韩视频精品一区| 国产av国产精品国产| 纯流量卡能插随身wifi吗| 精品人妻熟女毛片av久久网站| 国产午夜精品久久久久久一区二区三区| 一级二级三级毛片免费看| 国产日韩欧美视频二区| 国产成人freesex在线| 日本午夜av视频| 最近中文字幕高清免费大全6| 国产精品久久久久久精品古装| 亚洲久久久国产精品| 深夜a级毛片| 天美传媒精品一区二区| 精品人妻一区二区三区麻豆| 永久免费av网站大全| 高清av免费在线| 亚洲欧美中文字幕日韩二区| 91精品一卡2卡3卡4卡| 国产黄频视频在线观看| 一区二区三区四区激情视频| 丝袜喷水一区| 我的老师免费观看完整版| 又大又黄又爽视频免费| 亚洲精华国产精华液的使用体验| 日本欧美视频一区| 久久婷婷青草| 性高湖久久久久久久久免费观看| 三级国产精品欧美在线观看| 美女主播在线视频| 欧美另类一区| 美女cb高潮喷水在线观看| 人妻少妇偷人精品九色| av线在线观看网站| 十八禁网站网址无遮挡 | 校园人妻丝袜中文字幕| 熟女人妻精品中文字幕| 国产一级毛片在线| 又爽又黄a免费视频| 亚洲久久久国产精品| 精品国产露脸久久av麻豆| 亚洲精品,欧美精品| 午夜激情久久久久久久| 久久午夜福利片| 老司机亚洲免费影院| 中文乱码字字幕精品一区二区三区| 黑人巨大精品欧美一区二区蜜桃 | 一级毛片电影观看| 亚州av有码| 观看美女的网站| 国产精品久久久久久精品电影小说| 麻豆乱淫一区二区| 欧美成人精品欧美一级黄| 国产精品国产av在线观看| 老司机影院毛片| 国产精品成人在线| 一本—道久久a久久精品蜜桃钙片| 国产精品一区二区在线不卡| 大香蕉97超碰在线| 女的被弄到高潮叫床怎么办| h视频一区二区三区| 全区人妻精品视频| 美女主播在线视频| av天堂久久9| 国产乱来视频区| 精品久久久噜噜| 亚洲激情五月婷婷啪啪| 午夜日韩欧美国产| 久久毛片免费看一区二区三区| tocl精华| 国产免费一区二区三区四区乱码| 在线观看一区二区三区激情| 69av精品久久久久久 | 9热在线视频观看99| 国产亚洲av片在线观看秒播厂| 久久久久国产精品人妻一区二区| 午夜福利视频精品| 狂野欧美激情性bbbbbb| 免费少妇av软件| 一个人免费在线观看的高清视频 | 久久精品成人免费网站| 久久久精品国产亚洲av高清涩受| 久久ye,这里只有精品| 无限看片的www在线观看| 中文字幕制服av| 国产真人三级小视频在线观看| 国产在线免费精品| 亚洲五月婷婷丁香| 欧美另类一区| 两性午夜刺激爽爽歪歪视频在线观看 | 男女下面插进去视频免费观看| 欧美黄色淫秽网站| 两人在一起打扑克的视频| 欧美激情高清一区二区三区| 久久精品久久久久久噜噜老黄| 男人操女人黄网站| 2018国产大陆天天弄谢| 欧美人与性动交α欧美软件| av线在线观看网站| 成人免费观看视频高清| 这个男人来自地球电影免费观看| 爱豆传媒免费全集在线观看| 亚洲 欧美一区二区三区| 九色亚洲精品在线播放| 波多野结衣av一区二区av| 国产精品av久久久久免费| 中文字幕人妻熟女乱码| 丰满迷人的少妇在线观看| 97精品久久久久久久久久精品| 精品一区二区三卡| 欧美xxⅹ黑人| 各种免费的搞黄视频| 在线观看舔阴道视频| 亚洲精品国产av蜜桃| 日本精品一区二区三区蜜桃| 久久免费观看电影| 精品欧美一区二区三区在线| 91九色精品人成在线观看| 狠狠精品人妻久久久久久综合| 久久亚洲国产成人精品v| 久久免费观看电影| √禁漫天堂资源中文www| 国产老妇伦熟女老妇高清| 久久国产精品影院| 欧美黑人精品巨大| 亚洲精品自拍成人| √禁漫天堂资源中文www| 老熟妇乱子伦视频在线观看 | 涩涩av久久男人的天堂| 国产精品.久久久| 99香蕉大伊视频| xxxhd国产人妻xxx| 亚洲欧美激情在线| 国产又色又爽无遮挡免| 伊人亚洲综合成人网| 爱豆传媒免费全集在线观看| 一区二区三区乱码不卡18| 成人国产一区最新在线观看| 五月开心婷婷网| 在线观看免费午夜福利视频| 天堂8中文在线网| 超碰97精品在线观看| 亚洲av美国av| 久9热在线精品视频| 日日爽夜夜爽网站| 18禁观看日本| 国产伦人伦偷精品视频| 少妇猛男粗大的猛烈进出视频| 久久精品aⅴ一区二区三区四区| av超薄肉色丝袜交足视频| 中文字幕av电影在线播放| 一区福利在线观看| 亚洲免费av在线视频| 亚洲专区国产一区二区| 99国产综合亚洲精品| 亚洲av电影在线观看一区二区三区| 涩涩av久久男人的天堂| 日本av免费视频播放| 老司机亚洲免费影院| 日本黄色日本黄色录像| 亚洲精品美女久久av网站| 日本vs欧美在线观看视频| 亚洲三区欧美一区| 国产精品 国内视频| 王馨瑶露胸无遮挡在线观看| 国产免费av片在线观看野外av| 视频区图区小说| 亚洲欧美日韩另类电影网站| 大型av网站在线播放| 国精品久久久久久国模美| avwww免费| 国产成人一区二区三区免费视频网站| 亚洲精品中文字幕在线视频| 汤姆久久久久久久影院中文字幕| 欧美精品啪啪一区二区三区 | 免费av中文字幕在线| 午夜免费成人在线视频| 免费在线观看黄色视频的| 一级a爱视频在线免费观看| 亚洲人成77777在线视频| 91大片在线观看| 天天影视国产精品| 视频区图区小说| 日韩一区二区三区影片| 纵有疾风起免费观看全集完整版| 精品高清国产在线一区| 亚洲精品美女久久av网站| 亚洲国产中文字幕在线视频| 一个人免费看片子| av福利片在线| 婷婷成人精品国产| 欧美亚洲 丝袜 人妻 在线| 中文字幕精品免费在线观看视频| 啦啦啦 在线观看视频| 国产亚洲午夜精品一区二区久久| 亚洲精品中文字幕在线视频| 午夜福利免费观看在线| 视频区图区小说| 午夜91福利影院| 男人添女人高潮全过程视频| 中国美女看黄片| 黑人欧美特级aaaaaa片| 亚洲精品国产区一区二| 亚洲久久久国产精品| 欧美成狂野欧美在线观看| 国产成人精品久久二区二区免费| 97在线人人人人妻| 91av网站免费观看| 一级片'在线观看视频| 亚洲精品第二区| 亚洲男人天堂网一区| bbb黄色大片| 搡老熟女国产l中国老女人| 欧美精品高潮呻吟av久久| 亚洲一区二区三区欧美精品| 精品少妇久久久久久888优播| 丝瓜视频免费看黄片| 狠狠婷婷综合久久久久久88av| 亚洲第一欧美日韩一区二区三区 | 老司机深夜福利视频在线观看 | 免费久久久久久久精品成人欧美视频| 97精品久久久久久久久久精品| 免费在线观看完整版高清| 欧美日韩福利视频一区二区| 午夜日韩欧美国产| 久久国产精品男人的天堂亚洲| av超薄肉色丝袜交足视频| 天堂8中文在线网| 99久久人妻综合| 大型av网站在线播放| 亚洲国产精品一区三区| 12—13女人毛片做爰片一| 黄色视频,在线免费观看| 欧美日韩亚洲国产一区二区在线观看 | 午夜免费成人在线视频| 日本a在线网址| 久久久精品94久久精品| 超碰97精品在线观看| 亚洲人成77777在线视频| 成人三级做爰电影| 建设人人有责人人尽责人人享有的| 十八禁网站网址无遮挡| 日韩一卡2卡3卡4卡2021年| a级毛片黄视频| 精品少妇黑人巨大在线播放| 性少妇av在线| 亚洲国产精品成人久久小说| 蜜桃在线观看..| 久久香蕉激情| 人妻人人澡人人爽人人| 日韩欧美一区视频在线观看| 久久性视频一级片| 99久久99久久久精品蜜桃| av欧美777| 美女主播在线视频| 岛国在线观看网站| 高清视频免费观看一区二区| 亚洲国产av影院在线观看| 国产精品99久久99久久久不卡| 人人妻人人爽人人添夜夜欢视频| 精品国产一区二区久久| 亚洲专区字幕在线| 考比视频在线观看| 久久久精品国产亚洲av高清涩受| 亚洲欧美清纯卡通| 午夜福利,免费看| 亚洲欧美成人综合另类久久久| 亚洲精品国产av成人精品| 亚洲成人国产一区在线观看| 国产成人免费观看mmmm| 久久国产精品影院| 人人妻人人澡人人爽人人夜夜| 精品久久蜜臀av无| 激情视频va一区二区三区| 一个人免费在线观看的高清视频 | 欧美国产精品va在线观看不卡| 久久天堂一区二区三区四区| 18在线观看网站| 日日夜夜操网爽| 国产主播在线观看一区二区| 国产熟女午夜一区二区三区| 不卡av一区二区三区| 国产一卡二卡三卡精品| 久久精品久久久久久噜噜老黄| 下体分泌物呈黄色| 岛国毛片在线播放| 色综合欧美亚洲国产小说| 热re99久久精品国产66热6| 久久国产精品大桥未久av| 精品欧美一区二区三区在线| 国产亚洲av高清不卡| 国产在线免费精品| 中文精品一卡2卡3卡4更新| 天天躁日日躁夜夜躁夜夜| 大片电影免费在线观看免费| 无限看片的www在线观看| 久久亚洲国产成人精品v| 国产一区二区在线观看av| 亚洲欧美日韩另类电影网站| 亚洲久久久国产精品| 操出白浆在线播放| 老司机深夜福利视频在线观看 | 91精品伊人久久大香线蕉| 男女床上黄色一级片免费看| 伦理电影免费视频| 国产视频一区二区在线看| 久久久水蜜桃国产精品网| 国产亚洲欧美精品永久| 亚洲中文字幕日韩| 肉色欧美久久久久久久蜜桃| 国产精品 欧美亚洲| 午夜视频精品福利| 日韩熟女老妇一区二区性免费视频| 高清视频免费观看一区二区| 操美女的视频在线观看| 精品国产超薄肉色丝袜足j| av有码第一页| 亚洲精品国产色婷婷电影| 国产亚洲精品第一综合不卡| 中国美女看黄片| 王馨瑶露胸无遮挡在线观看| 久热这里只有精品99| 精品久久久久久电影网| 美国免费a级毛片| kizo精华| 777米奇影视久久| 动漫黄色视频在线观看| 狠狠婷婷综合久久久久久88av| 成人18禁高潮啪啪吃奶动态图| videos熟女内射| 国精品久久久久久国模美| 国产精品熟女久久久久浪| 国产在线免费精品| 亚洲欧洲精品一区二区精品久久久| 午夜免费成人在线视频| 成人黄色视频免费在线看| 黄色视频在线播放观看不卡| 久久中文字幕一级| 日日摸夜夜添夜夜添小说| 国产免费av片在线观看野外av| 日本撒尿小便嘘嘘汇集6| 超碰97精品在线观看| av网站免费在线观看视频| 99九九在线精品视频| 人人妻人人爽人人添夜夜欢视频| 欧美一级毛片孕妇| 国产深夜福利视频在线观看| 国产精品成人在线| 97人妻天天添夜夜摸| 99国产综合亚洲精品| 亚洲色图综合在线观看| 热99久久久久精品小说推荐| 亚洲成国产人片在线观看| 每晚都被弄得嗷嗷叫到高潮| 久久久久国内视频| 亚洲一区中文字幕在线| 国产精品av久久久久免费| 一区二区三区四区激情视频| 国产一卡二卡三卡精品| 免费久久久久久久精品成人欧美视频| 精品人妻1区二区| 波多野结衣一区麻豆| 97精品久久久久久久久久精品| 天堂中文最新版在线下载| 国产免费视频播放在线视频| 欧美国产精品va在线观看不卡| a级片在线免费高清观看视频| av片东京热男人的天堂| 亚洲国产欧美日韩在线播放| 国产黄色免费在线视频| 久久精品亚洲熟妇少妇任你| 90打野战视频偷拍视频| xxxhd国产人妻xxx| 另类亚洲欧美激情| av又黄又爽大尺度在线免费看| 91字幕亚洲| 欧美成狂野欧美在线观看| 日韩 亚洲 欧美在线| 欧美激情高清一区二区三区| 少妇猛男粗大的猛烈进出视频| 99精品欧美一区二区三区四区| 正在播放国产对白刺激| 精品视频人人做人人爽| 久久久久视频综合| 考比视频在线观看| 啦啦啦在线免费观看视频4| 美女脱内裤让男人舔精品视频| 成年人午夜在线观看视频| 午夜福利在线观看吧| 美女视频免费永久观看网站| 99精品久久久久人妻精品| 男女午夜视频在线观看| 一级毛片电影观看| 亚洲久久久国产精品| 99久久99久久久精品蜜桃| 亚洲一码二码三码区别大吗| 国产无遮挡羞羞视频在线观看| 视频区图区小说| 午夜影院在线不卡| 两个人看的免费小视频| 黑丝袜美女国产一区| 欧美日韩av久久| 制服诱惑二区| 日韩制服丝袜自拍偷拍| 精品欧美一区二区三区在线| 亚洲国产精品999| 啦啦啦在线免费观看视频4| 中文字幕另类日韩欧美亚洲嫩草| 欧美精品亚洲一区二区| 免费观看a级毛片全部| 亚洲伊人色综图| 自拍欧美九色日韩亚洲蝌蚪91| 国产免费av片在线观看野外av| 国产99久久九九免费精品| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲av高清不卡| 国产国语露脸激情在线看| 啦啦啦免费观看视频1| 91大片在线观看| videosex国产| 777久久人妻少妇嫩草av网站| 国产欧美日韩一区二区精品| 国产一区有黄有色的免费视频| 欧美少妇被猛烈插入视频| 十八禁网站免费在线| 国产成+人综合+亚洲专区| 色精品久久人妻99蜜桃| 久久久久久久大尺度免费视频| 国产区一区二久久| 国产欧美日韩一区二区三 | 国产精品影院久久| 久久久精品94久久精品| videosex国产| 777久久人妻少妇嫩草av网站| 一区在线观看完整版| 在线观看人妻少妇| 男女国产视频网站| 五月开心婷婷网| 成年人午夜在线观看视频| tocl精华| 巨乳人妻的诱惑在线观看| 亚洲国产精品一区三区| 国产三级黄色录像| 丁香六月欧美| 国产成人精品在线电影| 丁香六月欧美| 亚洲精品美女久久av网站| 亚洲 欧美一区二区三区| 久久久久精品人妻al黑| 久久青草综合色| 捣出白浆h1v1| 一区二区三区四区激情视频| 在线天堂中文资源库| 精品国产乱子伦一区二区三区 | 高潮久久久久久久久久久不卡| 高清视频免费观看一区二区| 黑人欧美特级aaaaaa片| 高清视频免费观看一区二区| xxxhd国产人妻xxx| 女人久久www免费人成看片| 大片免费播放器 马上看| 在线观看免费视频网站a站| 欧美中文综合在线视频| 少妇被粗大的猛进出69影院| 久久久精品94久久精品| 亚洲国产精品999| 我的亚洲天堂| 女人久久www免费人成看片| 国产精品久久久久久精品电影小说| 别揉我奶头~嗯~啊~动态视频 | 久久久久久久精品精品| 亚洲av日韩在线播放| 12—13女人毛片做爰片一| 一级片'在线观看视频| 啦啦啦 在线观看视频| 岛国在线观看网站| 久久精品久久久久久噜噜老黄| 国产主播在线观看一区二区| 国产成人精品久久二区二区免费| 国产黄色免费在线视频| 精品亚洲成国产av| 色视频在线一区二区三区| 男女免费视频国产| 纯流量卡能插随身wifi吗| 多毛熟女@视频| 国产一区二区三区在线臀色熟女 | 老司机午夜十八禁免费视频| 蜜桃在线观看..| 人妻人人澡人人爽人人| 一区二区三区精品91| 国产精品麻豆人妻色哟哟久久| 黑人欧美特级aaaaaa片| 亚洲精品中文字幕一二三四区 | 欧美日韩中文字幕国产精品一区二区三区 | 国产免费一区二区三区四区乱码| 亚洲精品久久久久久婷婷小说| 国产野战对白在线观看| 高潮久久久久久久久久久不卡| 男女下面插进去视频免费观看| 亚洲精品国产av蜜桃| www.熟女人妻精品国产| 午夜福利在线免费观看网站| 丰满饥渴人妻一区二区三| 大码成人一级视频| 亚洲中文av在线| 一二三四在线观看免费中文在| 午夜老司机福利片| 国产伦理片在线播放av一区| 热99国产精品久久久久久7| 欧美日韩福利视频一区二区| 在线永久观看黄色视频| 一本大道久久a久久精品| 午夜精品久久久久久毛片777| 久久精品aⅴ一区二区三区四区| netflix在线观看网站| 亚洲美女黄色视频免费看| 国产淫语在线视频| 国产1区2区3区精品| 最近最新免费中文字幕在线| 丁香六月天网| av在线播放精品| 免费看十八禁软件| 国产精品久久久av美女十八| 免费在线观看日本一区| 色精品久久人妻99蜜桃| 一个人免费看片子| 丁香六月欧美| 亚洲人成电影观看| 男人添女人高潮全过程视频| 亚洲国产欧美网| 免费人妻精品一区二区三区视频| 日韩视频在线欧美| 99久久99久久久精品蜜桃| 欧美成狂野欧美在线观看| 久久亚洲精品不卡| 午夜福利影视在线免费观看| 最黄视频免费看| 蜜桃在线观看..| 在线十欧美十亚洲十日本专区| 国产成人一区二区三区免费视频网站| 亚洲av美国av| 黄色视频在线播放观看不卡| 免费在线观看黄色视频的| 亚洲欧美色中文字幕在线| av天堂久久9| 色94色欧美一区二区| 五月天丁香电影|