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

    炸藥化學(xué)爆炸環(huán)境中模擬钚氣溶膠源項實驗研究

    2016-10-14 02:15:04劉文杰馬慶鵬田立智王鵬來
    兵工學(xué)報 2016年3期
    關(guān)鍵詞:氣溶膠炸藥容器

    劉文杰,馬慶鵬,田立智,王鵬來

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

    炸藥化學(xué)爆炸環(huán)境中模擬钚氣溶膠源項實驗研究

    劉文杰,馬慶鵬,田立智,王鵬來

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

    采用高密封爆炸容器為實驗平臺,以銀、鎢、錫、鈰、釩5種金屬為钚的替代材料,開展爆轟加載條件下金屬材料氣溶膠源項實驗,對炸藥化學(xué)爆炸事故條件下钚氣溶膠源項開展實驗?zāi)M研究。利用8級沖擊式氣溶膠采樣器對金屬氣溶膠進(jìn)行間隔采樣,采用電感耦合等離子體質(zhì)譜分析方法對空氣動力學(xué)直徑10 μm以下的5種金屬氣溶膠源項分布進(jìn)行分析,并與美國钚氣溶膠外場擴(kuò)散實驗結(jié)果進(jìn)行對比。研究結(jié)果表明:金屬銀在實驗室條件下可實現(xiàn)對钚氣溶膠源項的模擬,其歸一化積累質(zhì)量分布與钚材料具有較高的一致性,其他4種材料均與钚差距較大。進(jìn)一步對銀氣溶膠進(jìn)行時間間隔采樣,并通過能譜分析及掃描電鏡等表征手段研究了銀氣溶膠源項隨時間的演化特性。結(jié)果表明:爆轟后初次采樣時,1.1~2.1 μm粒徑范圍的銀氣溶膠含量最高;二次采樣時,銀氣溶膠發(fā)生了顯著的凝并和沉降。

    兵器科學(xué)與技術(shù);钚氣溶膠;銀;源項;歸一化積累質(zhì)量分布;粒徑分布

    0 引言

    钚是一種化學(xué)性質(zhì)極為活潑的高毒性放射性核素,在其眾多的同位素中,239Pu是各類核裝置最常用的裂變材料。當(dāng)核裝置遭遇意外化學(xué)爆炸或燃燒等事故時,裝置中未達(dá)到臨界條件的钚材料(特指239Pu,下同)將在高溫高壓環(huán)境中液化并發(fā)生劇烈氧化反應(yīng),反應(yīng)產(chǎn)生的钚氧化物將以氣溶膠形式懸浮在空氣中并向外界擴(kuò)散[1-2]。由于钚氣溶膠粒徑極?。▉單⒚?、微米量級)極易被人體吸入,因此與常態(tài)钚相比可對人體產(chǎn)生更大的危害[3-4]。美國能源部及各武器實驗室的評估認(rèn)定空氣動力學(xué)直徑(AD)在2 μm以下的钚氣溶膠均可被吸入到肺深部,亞微米級微??赏高^毛細(xì)血管進(jìn)入血液循環(huán)系統(tǒng)并在肝臟和骨骼中沉積,并將钚氣溶膠顆粒造成生物體可吸入危害的AD上限定為10 μm,人體允許吸入钚的最高劑量為0.3 μg,急性致死量為4.3 mg[5].

    20世紀(jì)60年代前后,美國發(fā)生多起武器級核裝置因高能炸藥化爆引起的钚氣溶膠污染事故,促使其開展了一系列炸藥化學(xué)爆炸環(huán)境中钚氣溶膠源項及擴(kuò)散的研究工作。1963年5月~6月期間,美國和英國聯(lián)合開展了 Operation Roller Coaster (ORC)外場試驗[6],獲得了化爆條件下钚氣溶膠源項及其擴(kuò)散過程系統(tǒng)、詳細(xì)的原始數(shù)據(jù)。在這次實驗結(jié)束后,美國禁止了一切可能造成钚材料擴(kuò)散至大氣環(huán)境的外場實驗,ORC實驗為日后的钚氣溶膠擴(kuò)散數(shù)值模擬提供了最重要的源項數(shù)據(jù)。而后,相關(guān)的研究工作也轉(zhuǎn)入了實驗室實驗和數(shù)值模擬。在實驗室研究方面,美國桑迪亞國家實驗室(SNL)的Sagartz[7]通過PBX-9404炸藥加載銀球殼實驗?zāi)M核裝置在事故條件下钚氣溶膠源項的生成規(guī)律,研究了不同炸藥構(gòu)型條件下銀樣品氣溶膠化的轉(zhuǎn)化率,分析了炸藥爆能對銀氣溶膠轉(zhuǎn)化率的影響,但未說明使用金屬銀作為钚替代材料的依據(jù),也并未將銀氣溶膠源項數(shù)據(jù)與钚進(jìn)行對比。郝樊華等[8]開展了爆轟條件下銀氣溶膠的源項實驗,但其研究主要側(cè)重于銀氣溶膠的粒徑質(zhì)量分布(SMD)及分凝效應(yīng),并未與钚氣溶膠相聯(lián)系。目前,氣溶膠擴(kuò)散的數(shù)值計算,主要以美國洛斯阿拉莫斯國家實驗室(LANL)、勞倫斯利弗摩爾國家實驗室 (LLNL)、SNL三大國家實驗室開發(fā)的數(shù)值計算模型為研究手段,包括HOTSPOT[9]模型、DIFOUT[10]、MACCS2[11]和ERAD[12]模型等,前3種是高斯煙云模型,后一種是采用蒙特卡洛方法的三維數(shù)值模型。

    在模擬核事故條件下钚氣溶膠的擴(kuò)散機(jī)制時,所有模型的輸入?yún)?shù)均依賴于事故產(chǎn)生的源項數(shù)據(jù)。然而出于安全及政治方面的考慮,直接開展钚材料化學(xué)爆炸實驗研究缺乏現(xiàn)實可行性,通過實驗室實驗,尋找一種能夠有效模擬钚氣溶膠源項的常規(guī)金屬替代材料是較為可行的一種研究方法。

    郝樊華等[8]曾開展過氣溶膠源項的探索研究,本文基于其研究思路,建立了一種實驗室密封環(huán)境條件下,通過爆轟加載手段研究金屬氣溶膠源項的實驗方法,開展多種金屬材料在化學(xué)爆炸條件下產(chǎn)生的氣溶膠源項實驗,對AD小于10 μm(PM10)的金屬氣溶膠源項分布特性進(jìn)行分析,首次采用常規(guī)金屬替代材料對钚材料氣溶膠歸一化積累質(zhì)量分布(CMD)的模擬。該研究結(jié)果對意外核事故的預(yù)防及應(yīng)急響應(yīng)措施的制定具有重大的意義,同時為钚氣溶膠擴(kuò)散數(shù)值模擬提供了輸入?yún)?shù)。

    1 實驗方法

    1.1模擬材料選取

    本文選用銀、鎢、錫、鈰和釩5種金屬作為氣溶膠實驗的樣品材料。其中:錫和銀均為低熔點高沸點金屬,除密度相對較低外二者基本的表觀物理性質(zhì)與钚有相似性;鎢與钚的密度較為接近;金屬鈰與金屬钚都具有低熔點高沸點、非對稱性晶體結(jié)構(gòu)、多種同素異形體共存、發(fā)生相變時都發(fā)生體積變化等特性,諸多物相性質(zhì)上的相似性很大程度上與兩種金屬相近的亞殼層電子結(jié)構(gòu)有關(guān),钚與鈰均處于元素周期表中f電子向定域態(tài)轉(zhuǎn)變的位置;金屬釩與钚的沸點亦較為接近。上述5種金屬及钚的基本物理性質(zhì)如表1[13]所示。

    表1 氣溶膠實驗相關(guān)金屬物理性質(zhì)[13]Tab.1 Physical properties of metals used in experiment[13]

    1.2實驗方法及工況

    實驗系統(tǒng)由2.5 kg梯恩梯當(dāng)量高密封爆炸容器、氣溶膠采集系統(tǒng)(多級沖擊采樣器及質(zhì)量流量計)和抽氣系統(tǒng)(真空泵、角閥及波紋管道)組成,如圖1所示。

    圖1 爆轟加載金屬銀氣溶膠實驗裝置示意圖Fig.1 Simplified scheme of experimental set-up

    氣溶膠采樣器采用安德森TE-20-800型8級沖擊式采樣器,對AD在0~10 μm范圍內(nèi)的氣溶膠顆粒進(jìn)行粒徑分段采集,圖2(a)[14]為該采樣器的實物圖。該采樣器具有一個10.0 μm粒徑的前置切割級,用來切割、分離所有粒徑大于等于10.0 μm的顆粒,如圖2(b)[14]所示。中間8級切割粒徑分別為9.0 μm、5.8 μm、4.7 μm、3.3 μm、2.1 μm、1.1 μm、0.7 μm、0.4 μm,對應(yīng)的沖擊板編號為0~7,分別采集9.0~10.0 μm至0.4~0.7 μm粒徑的氣溶膠顆粒;后置級“Filter”用來收集所有粒徑小于0.4 μm的氣溶膠顆粒。實驗中選用對金屬顆粒物吸附性能最優(yōu)的高密度石英纖維濾膜(美國Tisch公司產(chǎn))作為氣溶膠樣品采集基底。實驗過程中,通過真空抽氣系統(tǒng)的可調(diào)角閥控制抽速,并通過質(zhì)量流量計監(jiān)測采樣流量,將其限制為28.3±1.5 L/min,以保證采樣器的粒徑切割精度。

    圖2 安德森TE-20-800型8級沖擊式采樣器[14]Fig.2 Anderson 8 stage non-viable cascade impactor[14]

    實驗主要采用圖3所示的平面對碰加載模型對銀片進(jìn)行爆轟加載,該構(gòu)型與美國ORC-DT钚氣溶膠外場擴(kuò)散實驗加載構(gòu)型類似。

    在美國ORC-DT實驗中,采用平面加載方式,裝置由1.3 kg钚材料及47.6 kg高能炸藥組成,放置在離地高度8 feet的鋼板上方1 foot高度[13]。本文實驗采用的炸藥單片尺寸為φ100 mm×20 mm,在炸藥中間φ70 mm范圍內(nèi)放置要研究的金屬材料,通過厚度控制夾在炸藥之間的金屬質(zhì)量最大不超過50 g(錫密度最小,厚度接近5 mm;鎢密度最大,厚約0.7 mm);起爆裝置使用BL-21雷管,傳爆藥為φ32 mm×11 mm規(guī)格的泰安;裝置外用厚5 mm的有機(jī)玻璃外殼進(jìn)行約束。在炸藥使用方面,主要采用RHT-901高能炸藥分別對5種材料進(jìn)行爆轟加載,對于銀材料采用RHT-901和奧克托今(HMX)兩種炸藥對其進(jìn)行爆轟加載,具體實驗參數(shù)如表2所示。

    圖3 平面對碰爆轟加載裝置示意圖Fig.3 Schematic diagram of detonative device

    表2 實驗參數(shù)Tab.2 Experimental conditions

    每發(fā)次實驗前,對容器內(nèi)壁進(jìn)行清洗,實驗狀態(tài)為常壓,不進(jìn)行抽真空處理。實驗中,盡量減少輔助裝置,并避免使用可能對實驗結(jié)果產(chǎn)生較大影響的輔助材料(加載裝置采用棉繩懸掛,雷管套及樣品約束件均采用有機(jī)玻璃)。

    爆轟后,目標(biāo)收集物為銀的氧化物。對于濾膜樣品中銀含量的分析,首先,對濾膜及其附著物進(jìn)行電化學(xué)消解,然后,采用電感耦合等離子體質(zhì)譜分析(ICP-MS)方法對銀元素含量進(jìn)行測定,測量精度可達(dá)0.01 μg.

    2 實驗結(jié)果與討論

    2.1實驗重復(fù)性驗證

    圖4 銀氣溶膠CMDFig.4 CMD of silver aerosol

    對比3發(fā)實驗結(jié)果可知,第1發(fā)、第2發(fā)實驗在3.3 μm以上粒徑段的 CMD基本吻合,在小于2.1 μm粒徑段有一定的差異,由于該段氣溶膠顆粒粒徑較小,易受實驗環(huán)境的微小差異及氣流擾動的影響,故存在一定的差異。第3發(fā)實驗結(jié)果與其他兩發(fā)稍有差異,表現(xiàn)在小粒徑氣溶膠含量更高,但整體分布趨勢相同。綜上可知,本文采用的實驗方法具有較好的重復(fù)性,采樣結(jié)果具有一定的隨機(jī)性,但在可接受范圍內(nèi)。

    2.2氣溶膠源項分析

    在裝置起爆2 min之內(nèi),開啟氣溶膠采集系統(tǒng)對容器內(nèi)氣體進(jìn)行取樣,每次取樣10 min,單次取樣抽氣量為 283 L,約占容器總體積(7 300 L)的3.877%.電化學(xué)消解后能譜分析給出容器內(nèi)氣相金屬氣溶膠SMD范圍,利用原子吸收光譜方法給出氣相金屬氣溶膠采集總量,進(jìn)而推算出爆轟后容器內(nèi)被測金屬氣溶膠的PM10總量(假定爆轟后容器內(nèi)氣相爆轟產(chǎn)物均勻分布)。每發(fā)次實驗均進(jìn)行了多次采樣,各采樣起始時間如表3所示。

    表3 各發(fā)次實驗采樣時間Tab.3 Sampling time

    圖5 金屬氣溶膠CMDFig.5 CMD of metal aerosol

    對實驗得到的各金屬氣溶膠CMD進(jìn)行分析,并與美國ORC-DT钚氣溶膠外場擴(kuò)散實驗結(jié)果[15]進(jìn)行對比。

    圖5為各發(fā)次實驗采樣得到的金屬氣溶膠CMD結(jié)果。由圖5(a)可知:5種金屬材料爆轟后第1次采樣得到的CMD結(jié)果與美國ORC-DT钚氣溶膠實驗結(jié)果均呈現(xiàn)較大的偏差,金屬鈰、錫及鎢氣溶膠CMD明顯高于钚氣溶膠;在RHT-901炸藥加載條件下,金屬銀及釩氣溶膠CMD結(jié)果與钚氣溶膠較為接近,但整體水平同樣偏高;HMX炸藥加載下,銀氣溶膠CMD結(jié)果與RHT-901加載時的結(jié)果存在較大差異,在各粒徑范圍的數(shù)值前者均高于后者,這表明炸藥性能對金屬氣溶膠源項分布具有較大的影響。進(jìn)一步對各材料實驗的多次采樣結(jié)果進(jìn)行分析,圖5(b)為銀氣溶膠4次采樣得到的CMD結(jié)果與ORC-DT钚氣溶膠實驗結(jié)果對比,可看出,HMX炸藥加載條件下兩次采樣結(jié)果均與钚氣溶膠差距較大。2.1~10.0 μm粒徑范圍內(nèi),爆轟后1.5 min的結(jié)果偏低,而爆轟后106.5 min的結(jié)果偏高,由于CMD為積分質(zhì)量比率,因此可以斷定,在這段時間內(nèi)大粒徑銀氣溶膠發(fā)生了明顯的沉降,導(dǎo)致較小粒徑氣溶膠相對含量上升。同時,由圖5(b)可知:RHT-901炸藥加載條件下,2.1~10.0 μm粒徑范圍內(nèi),第2次采樣結(jié)果(爆轟后106.5 min)較第1次采樣(爆轟后1.5 min)更為接近钚氣溶膠結(jié)果;而在0~2.1 μm粒徑范圍內(nèi),第1次采樣結(jié)果更為接近钚氣溶膠結(jié)果。對此,可得到兩個結(jié)論:1)在兩次采樣時間范圍內(nèi),銀氣溶膠發(fā)生了一定程度的凝并,即較小粒徑氣溶膠顆粒聚合形成較大顆粒,導(dǎo)致第2次采樣結(jié)果小粒徑氣溶膠含量降低。2)通過對比分析可知,除銀以外,其他幾種金屬在爆轟加載條件下的氣溶膠CMD與ORC-DT實驗獲得的钚氣溶膠CMD差異較大;而在RHT-901炸藥對碰加載條件下,對銀氣溶膠進(jìn)行不同時段的采樣,可能實現(xiàn)對ORC-DT實驗钚氣溶膠CMD的模擬,進(jìn)而可為外場擴(kuò)散數(shù)值模擬提供有效的源項參數(shù)。對此,將在下一階段的工作中繼續(xù)開展相關(guān)工作予以確認(rèn)。對于金屬釩,圖5(c)給出了兩次采樣得到的CMD結(jié)果,可看出,釩氣溶膠兩次采樣結(jié)果均不能對钚氣溶膠源項進(jìn)行有效模擬,同樣,由圖5(d)可知,其他3種金屬—鎢、錫、鈰共計9次采樣得到的結(jié)果均與钚氣溶膠源項差異較大。

    在氣溶膠可吸入釋放率(RRF,PM10氣溶膠總量占金屬樣品總質(zhì)量比率)方面,考慮到裝置在容器中心爆轟,可認(rèn)為氣溶膠顆粒在容器內(nèi)呈均勻分布,通過抽樣體積與容器容積之間的比例關(guān)系(3.877%)可推算出爆轟后容器內(nèi)氣溶膠的總量,進(jìn)而與樣品質(zhì)量相比得到各金屬材料的RRF,如表4所示(均為RHT-901炸藥加載)。對比5種材料RRF可知,在相同爆轟加載條件下,銀、鎢、錫、鈰4種金屬RRF均在10-4數(shù)量級左右。金屬釩明顯高于其他4種金屬,這主要因為實驗中金屬釩初始狀態(tài)為粉末狀,較易氣溶膠化。美國ORC-DT實驗得到炸藥化爆條件下钚材料的RRF約為20%,因此可知,本研究采用的5種替代金屬材料在RRF方面均無法對钚材料進(jìn)行模擬,同時,對5種金屬材料熔點及RRF可知,二者并未呈現(xiàn)出明顯的關(guān)聯(lián)。

    綜上所述,本文所采用的5種金屬材料,唯有銀氣溶膠可以對钚氣溶膠CMD實現(xiàn)有效模擬;除釩粉外,其他4種金屬的氣溶膠RRF均遠(yuǎn)遠(yuǎn)低于钚,且與被加載金屬熔點無明顯關(guān)聯(lián)。此外,對于銀氣溶膠在多次采樣期間發(fā)生的凝并及沉降現(xiàn)象,還需要開展進(jìn)一步研究。

    表4 RHT-901加載下金屬氣溶膠RRFTab.4 RRF of metal aerosol under detonation of RHT-901

    2.2氣溶膠凝并與沉降分析

    前文分析得知,銀氣溶膠在0~2.1 μm粒徑范圍內(nèi)出現(xiàn)較為明顯的凝并現(xiàn)象,而在2.1~10.0 μm粒徑范圍內(nèi)有明顯的沉降效應(yīng)。對此,可從銀氣溶膠的SMD及沉降物的掃描電鏡結(jié)果予以分析。

    圖6為兩種炸藥加載條件下銀氣溶膠4次采樣得到的SMD.

    圖6 銀氣溶膠SMDFig.6 SMD of Ag aerosol

    由圖6可看出:兩種加載條件下,爆轟后第1次采樣中1.1~2.1 μm粒徑范圍內(nèi)的氣溶膠含量均為最高,該粒徑范圍內(nèi)的氣溶膠也是最容易被人體吸入的,對人體危害最大;第2次采樣時,0~2.1 μm粒徑范圍內(nèi)氣溶膠含量較爆轟后初次采樣均有明顯下降,考慮到此粒徑范圍內(nèi)氣溶膠質(zhì)量較輕不易沉降,因此可以判定該粒徑范圍內(nèi)的氣溶膠發(fā)生了較為明顯的凝并效應(yīng);在RHT-901炸藥加載狀態(tài)下,初次采樣中,0~2.1 μm粒徑范圍內(nèi)氣溶膠含量較高,后期凝并效應(yīng)也更強(qiáng)。此外,由圖6可知,在小粒徑氣溶膠凝并的同時,較大粒徑(2.1~10.0 μm)氣溶膠含量并未上升,反而同樣呈現(xiàn)出下降趨勢,這表明較大粒徑氣溶膠沉降效應(yīng)強(qiáng)于凝并效應(yīng),粒徑越大,沉降效應(yīng)越強(qiáng),9.0~10.0 μm粒徑范圍內(nèi)氣溶膠沉降趨勢最為明顯。

    進(jìn)一步對RHT-901炸藥加載條件下容器底部固體沉降物進(jìn)行分析,實驗后24 h,打開爆炸容器入孔門,收集容器內(nèi)的爆轟殘留顆粒,通過機(jī)械共振篩樣機(jī)(美國西蒙公司Mol 20型,38~125 μm共分8級)完成顆粒物分級篩樣,圖7給出了顆粒尺寸在53~63 μm及0~38 μm范圍內(nèi)的殘留物掃描電鏡圖像。

    圖7 爆轟殘留顆粒掃描電鏡圖像Fig.7 Scanning electron microscope results of residual particles at the bottom of the confinement vessel

    由圖7(a)~圖7(c)可知,顆粒尺寸在53~63 μm范圍內(nèi),銀氣溶膠多數(shù)以單一的大顆粒狀態(tài)存在,顆粒表面僅有少數(shù)小粒徑氣溶膠顆粒附著;由圖7(d)~圖7(f)可知,顆粒尺寸在0~38 μm范圍內(nèi),銀氣溶膠顆粒表面狀況較為復(fù)雜,由5 000倍放大圖(見圖7(f))可看出,在較大粒徑的銀氣溶膠顆粒上附著有大量的小尺寸顆粒。這表明,爆轟后有大量的小尺寸銀氣溶膠與較大尺寸的氣溶膠顆粒發(fā)生凝并,繼而隨之沉降到容器底部,而更大粒徑的銀氣溶膠由于重力較大快速沉降,其表面并無較多的小尺寸氣溶膠附著。

    3 結(jié)論

    本文建立了一種密封環(huán)境條件下對金屬材料進(jìn)行爆轟加載的氣溶膠源項實驗方法,開展了對開放環(huán)境下钚氣溶膠源項分布進(jìn)行模擬的氣溶膠源項實驗。通過對銀、鎢、錫、鈰、釩等5種金屬進(jìn)行爆轟加載,獲得了其PM10氣溶膠源項分布,并與文獻(xiàn)中钚氣溶膠源項進(jìn)行了對比分析,具體得到以下結(jié)論:

    1)采用高密封環(huán)境下對常規(guī)金屬進(jìn)行爆轟加載的氣溶膠源項實驗方法,在開放環(huán)境下钚氣溶膠源項模擬方面具有一定的可行性。

    2)采用兩種炸藥開展了5種金屬材料、共計6發(fā)次爆轟加載實驗,結(jié)果表明,金屬銀在RHT-901加載條件下產(chǎn)生的氣溶膠源項,其CMD與美國ORC-DT實驗中钚氣溶膠CMD具有較高的一致性。金屬氣溶膠的RRF與其熔點并無關(guān)聯(lián),5種金屬材料均不能對钚氣溶膠RRF進(jìn)行有效模擬。

    3)爆轟加載條件下銀氣溶膠源項實驗結(jié)果表明,爆轟后初次采樣得到1.1~2.1 μm粒徑范圍內(nèi)的氣溶膠含量最高,隨著時間推移,該段粒徑范圍內(nèi)氣溶膠含量迅速降低,PM10銀氣溶膠存在明顯的凝并及沉降效應(yīng)。小粒徑氣溶膠含量越高,凝并越顯著;粒徑越大,沉降效應(yīng)越強(qiáng)。

    (References)

    [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/2):90-103.

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

    [3] 劉文杰,胡八一,李慶忠.核事故條件下钚氣溶膠源項研究綜述[J].安全與環(huán)境學(xué)報,2011,11(5):259-263. LIU Wen-jie,HU Ba-yi,LI Qing-zhong.Review of plutonium aerosol source-term in nuclear accident[J].Journal of Safety and Environment,2011,11(5):259-263.(in Chinese)

    [4] Stewart K,Thomas D C,Terry J L,et al.Atomic weapons research establishment,AWRE 0-29/65[R].Abingdon,Oxfordshire:United Kingdom Atomic Energy Authority,1965.

    [5] Condit R H.Plutonium dispersal in fires:summary of what is known,UCRL-ID-114164[R].CA,US:Lawrence Livermore National Laboratory,1993.

    [6] Titus R W.Operation roller coaster:project officers report-project 2.4,POR-2504[R].US:United States Weather Bureau,1965.

    [7] Sagartz M J.Violent reaction source term study,SAND-943252 [R].Albuquerque,NM,US:Sandia National Laboratory,1995.

    [8] 郝樊華,劉曉亞,肖成建,等.爆轟條件下的氣溶膠粒徑分布及Ag質(zhì)量分?jǐn)?shù)與粒徑的關(guān)系[J].核技術(shù),2009,32(5):343-346. HAO Fan-hua,LIU Xiao-ya,XIAO Cheng-jian,et al.Size distributing of explosion-generated particles in a sealed container[J]. Nuclear Techniques,2009,32(5):343-346.(in Chinese)

    [9] Homann S G,Wilson D V.HOTSPOT health physics code:Version 8.0,UCRL-MA-118617[R].CA,US:Lawrence Livermore National Laboratory,1995.

    [10] Shreve J D.Operation roller coaster:scientific directors report,DASA-1644[R].Washington DC:Department of Defense,1965.

    [11] Steele C M,Wald T L,Chanin D I.Plutonium explosive dispersal modeling using the MACCS2 computer code,LA-UR-98-1901 [R].NM,US:Los Alamos National Laboratory,1998.

    [12] Boughton B A,Delaurentis J M.Description and validation of ERAD:an atmospheric dispersion model for high explosive detonations,SAND 92-2069-UC-702[R].Albuquerque,NM:Sandia National Laboratory,1992.

    [13] 劉光啟,馬連湘,項曙光.化學(xué)化工物性數(shù)據(jù)手冊:無機(jī)卷[M].北京:化學(xué)工業(yè)出版社,2013. LIU Guang-qi,MA Lian-xiang,XIANG Shu-guang.Handbook of chemistry physical properties data[M].Beijing:Chemical Industry Press,2013.(in Chinese)

    [14] Tisch Environmental Inc.Model 20-800 ambient cascade impactor(non-viable)operations manual[M].Miami:Tisch Environmental Inc,1999.

    [15] Stephens D R.Source terms for plutonium aerosolization from nuclear weapon accidents,UCRL-ID-119303[R].CA,US:Lawrence Livermore National Laboratory,1995.

    Experimental Simulation of Aerosol Source Term of Surrogate Plutonium Produced by High-explosive Detonation

    LIU Wen-jie,MA Qing-peng,TIAN Li-zhi,WANG Peng-lai
    (Institute of Fluid Physics,China Academy of Engineering Physics,Mianyang 621900,Sichuan,China)

    An experimental method is developed to simulate the plutonium aerosol source term in the circumstance of high-explosive detonation with ordinary surrogate metals on an aerosol test facility consisting of a confinement vessel and a non-viable aerosol sizing sampler.Five kinds of metals,including silver (Ag),tungsten(W),tin(Sn),cerium(Ce)and vanadium(V),are detonated in confinement vessel.Aerosol samples with aerodynamic diameter of smaller than 10 μm are collected.The cumulative mass and size distribution of each metal aerosolized in experiment are investigated and compared to the plutonium aerosol source term obtained from operation roller coaster-double tracks experiment of USA. The results show that the silver aerosol can simulate the plutonium aerosol source term in the laboratory environment,of which the cumulative mass distribution is in good agreement with that of plutonium,while the aerosol source terms of other four metals are quite different from that of plutonium.Two samples of the silver aerosol in confinement vessel are collected at different time after detonation,which are compared based on energy spectrum analysis and scanning electron microscope characterization.A solid sam-ple on the bottom of confinement vessel is also collected,and the scanning electron microscope method is used to analyze the size distribution of the sediments.The results show that the AD of silver aerosol obtained just after detonation is in the range from 1.1 μm to 2.1 μm,while the aerosol size is larger for the late sample,suggesting that settlement and combination of the silver aerosol take place with time in confinement vessel.

    ordnance science and technology;plutonium aerosol;silver;source term;cumulative mass distribution;size distribution

    O648.18

    A

    1000-1093(2016)03-0447-08

    10.3969/j.issn.1000-1093.2016.03.009

    2015-04-22

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

    劉文杰(1982—),男,副研究員。E-mail:wjliu@mail.ustc.edu.cn;馬慶鵬(1988—),男,助理研究員。E-mail:maqingpeng@caep.cn

    猜你喜歡
    氣溶膠炸藥容器
    “炸藥”驚魂
    Different Containers不同的容器
    議論火炸藥數(shù)字化制造
    氣溶膠傳播之謎
    難以置信的事情
    氣溶膠中210Po測定的不確定度評定
    四川盆地秋季氣溶膠與云的相關(guān)分析
    取米
    大氣氣溶膠成核監(jiān)測
    Al粉對炸藥爆炸加速能力的影響
    特大巨黑吊av在线直播| 久久久久久伊人网av| 国产成人一区二区在线| 在线免费观看不下载黄p国产| 男女国产视频网站| 自拍欧美九色日韩亚洲蝌蚪91 | 免费播放大片免费观看视频在线观看| 国产一区二区三区综合在线观看 | 精品亚洲乱码少妇综合久久| 99九九线精品视频在线观看视频| 最新中文字幕久久久久| 天天操日日干夜夜撸| 亚洲欧美成人精品一区二区| 亚洲不卡免费看| 亚洲欧美日韩另类电影网站| 性色avwww在线观看| 最近最新中文字幕免费大全7| 中文欧美无线码| 亚洲av欧美aⅴ国产| 99精国产麻豆久久婷婷| 欧美日本中文国产一区发布| 午夜福利网站1000一区二区三区| 久久国产乱子免费精品| 少妇的逼好多水| 在线观看美女被高潮喷水网站| 午夜福利视频精品| 日韩一区二区三区影片| 日韩视频在线欧美| .国产精品久久| 成人毛片60女人毛片免费| 国产精品成人在线| 亚洲精品乱久久久久久| 一级毛片电影观看| 国产亚洲最大av| 精品国产乱码久久久久久小说| 蜜桃久久精品国产亚洲av| 色婷婷久久久亚洲欧美| 国产精品久久久久久精品电影小说| 久久久久网色| 一区二区三区乱码不卡18| 亚洲四区av| 麻豆成人午夜福利视频| 乱系列少妇在线播放| 亚洲人与动物交配视频| 欧美另类一区| 国产男人的电影天堂91| 国产精品三级大全| 黑人猛操日本美女一级片| 18禁裸乳无遮挡动漫免费视频| 日韩中字成人| 免费高清在线观看视频在线观看| av.在线天堂| 成人免费观看视频高清| 亚洲国产日韩一区二区| 久久久久国产网址| 亚洲精华国产精华液的使用体验| 少妇被粗大的猛进出69影院 | 少妇被粗大猛烈的视频| av.在线天堂| 国产成人精品一,二区| 欧美日韩av久久| 久久青草综合色| 久久午夜综合久久蜜桃| 亚洲国产欧美在线一区| 寂寞人妻少妇视频99o| 在线播放无遮挡| 国产日韩欧美在线精品| 日韩一区二区三区影片| 免费黄色在线免费观看| 精品少妇久久久久久888优播| 一边亲一边摸免费视频| av在线观看视频网站免费| 不卡视频在线观看欧美| 成人美女网站在线观看视频| 久久久久久久大尺度免费视频| 五月天丁香电影| 最后的刺客免费高清国语| 在线观看国产h片| 3wmmmm亚洲av在线观看| 精品一区在线观看国产| 久久青草综合色| 日本欧美国产在线视频| 久久这里有精品视频免费| 亚洲av国产av综合av卡| 国产精品三级大全| 国产精品久久久久久久久免| 亚洲欧美成人精品一区二区| a级毛片在线看网站| 美女脱内裤让男人舔精品视频| 国产成人精品福利久久| 亚洲精品自拍成人| 色哟哟·www| 久久午夜福利片| 午夜激情福利司机影院| 男人添女人高潮全过程视频| 成年人免费黄色播放视频 | 亚洲国产精品一区二区三区在线| 两个人的视频大全免费| 成人国产av品久久久| 26uuu在线亚洲综合色| 欧美最新免费一区二区三区| 人人妻人人爽人人添夜夜欢视频 | 18禁动态无遮挡网站| 丰满饥渴人妻一区二区三| 亚洲欧美日韩另类电影网站| 国产精品99久久99久久久不卡 | 国产黄色视频一区二区在线观看| 国产精品熟女久久久久浪| 一区二区三区四区激情视频| 搡老乐熟女国产| 亚洲av电影在线观看一区二区三区| 91精品国产国语对白视频| 久久国内精品自在自线图片| 欧美 亚洲 国产 日韩一| 亚洲人与动物交配视频| 国产在线男女| 欧美xxⅹ黑人| 精品久久久久久电影网| 嫩草影院新地址| 晚上一个人看的免费电影| 久热这里只有精品99| 高清午夜精品一区二区三区| 日韩中字成人| 久久久久久久精品精品| 麻豆乱淫一区二区| 在线亚洲精品国产二区图片欧美 | av在线app专区| 久久精品熟女亚洲av麻豆精品| 成人免费观看视频高清| 国产黄色免费在线视频| 久久久久视频综合| 婷婷色综合www| 狠狠精品人妻久久久久久综合| 狂野欧美激情性xxxx在线观看| 一本久久精品| 国产一区有黄有色的免费视频| 国产精品人妻久久久久久| av免费观看日本| 亚洲国产精品一区三区| 黑人巨大精品欧美一区二区蜜桃 | xxx大片免费视频| 久久精品久久精品一区二区三区| 18禁动态无遮挡网站| 内射极品少妇av片p| 一级,二级,三级黄色视频| 一级片'在线观看视频| 夜夜骑夜夜射夜夜干| 亚洲人成网站在线播| 久久av网站| 婷婷色麻豆天堂久久| 亚洲av.av天堂| 啦啦啦视频在线资源免费观看| a级毛色黄片| 精品一品国产午夜福利视频| 日本黄大片高清| 亚洲怡红院男人天堂| 久久6这里有精品| 亚洲欧美日韩另类电影网站| av线在线观看网站| av福利片在线观看| 精品一区二区免费观看| av国产久精品久网站免费入址| 成年人午夜在线观看视频| 在线观看人妻少妇| 制服丝袜香蕉在线| 久久久久视频综合| 国产成人aa在线观看| 黄色欧美视频在线观看| 永久免费av网站大全| 欧美日韩视频高清一区二区三区二| 国产精品嫩草影院av在线观看| 91精品国产九色| 亚洲国产精品999| 男女边摸边吃奶| 国产亚洲91精品色在线| 王馨瑶露胸无遮挡在线观看| 国产精品一区二区在线不卡| 日本vs欧美在线观看视频 | 狠狠精品人妻久久久久久综合| 亚洲无线观看免费| 欧美bdsm另类| 中文字幕av电影在线播放| 免费黄网站久久成人精品| 在线精品无人区一区二区三| 五月伊人婷婷丁香| 卡戴珊不雅视频在线播放| 午夜精品国产一区二区电影| 免费大片黄手机在线观看| 欧美日韩在线观看h| 成人漫画全彩无遮挡| 日韩 亚洲 欧美在线| 美女中出高潮动态图| 成人毛片60女人毛片免费| 亚洲一级一片aⅴ在线观看| 亚洲国产最新在线播放| av在线app专区| 九九在线视频观看精品| 亚洲欧美精品自产自拍| 国产精品国产三级国产av玫瑰| 丰满少妇做爰视频| 中文在线观看免费www的网站| 少妇高潮的动态图| 亚洲中文av在线| 水蜜桃什么品种好| 一级av片app| 中文在线观看免费www的网站| 伦理电影大哥的女人| 午夜福利网站1000一区二区三区| 亚洲av国产av综合av卡| 九九久久精品国产亚洲av麻豆| 国产欧美另类精品又又久久亚洲欧美| 久久青草综合色| 久久精品夜色国产| 狠狠精品人妻久久久久久综合| 少妇人妻精品综合一区二区| 成人毛片60女人毛片免费| 女人久久www免费人成看片| 天天躁夜夜躁狠狠久久av| 国产欧美日韩综合在线一区二区 | 色94色欧美一区二区| 欧美精品一区二区大全| 99久久中文字幕三级久久日本| 搡老乐熟女国产| 最新中文字幕久久久久| 亚洲,欧美,日韩| 一级黄片播放器| 国产黄色视频一区二区在线观看| 欧美成人精品欧美一级黄| 亚洲国产精品国产精品| 男人添女人高潮全过程视频| 少妇精品久久久久久久| 久久国产精品男人的天堂亚洲 | 国产成人a∨麻豆精品| 在线 av 中文字幕| 性色av一级| 精品亚洲乱码少妇综合久久| 国产精品一区二区在线观看99| 少妇人妻 视频| 99热全是精品| 伊人亚洲综合成人网| 国产免费视频播放在线视频| 久久ye,这里只有精品| 成年av动漫网址| 男女啪啪激烈高潮av片| 一本一本综合久久| av天堂中文字幕网| 啦啦啦视频在线资源免费观看| 两个人免费观看高清视频 | tube8黄色片| 中文欧美无线码| 国产片特级美女逼逼视频| 成人免费观看视频高清| 日本爱情动作片www.在线观看| 欧美97在线视频| 久久精品久久久久久噜噜老黄| 少妇 在线观看| 大陆偷拍与自拍| 校园人妻丝袜中文字幕| 少妇被粗大的猛进出69影院 | 精品久久久久久久久亚洲| 美女xxoo啪啪120秒动态图| 免费在线观看成人毛片| 国产一区二区三区综合在线观看 | 亚洲av福利一区| 国产乱来视频区| 免费看不卡的av| 国产精品偷伦视频观看了| 成人影院久久| 狠狠精品人妻久久久久久综合| 亚洲av二区三区四区| 国产熟女午夜一区二区三区 | 女性生殖器流出的白浆| 国内揄拍国产精品人妻在线| 久久国产亚洲av麻豆专区| 欧美成人午夜免费资源| 在线天堂最新版资源| av又黄又爽大尺度在线免费看| 国产亚洲精品久久久com| 成人综合一区亚洲| 免费人妻精品一区二区三区视频| 亚洲精品久久午夜乱码| 精品99又大又爽又粗少妇毛片| 国产在线视频一区二区| 精品久久久久久久久亚洲| 亚洲欧洲日产国产| 免费黄色在线免费观看| 国产精品人妻久久久久久| 你懂的网址亚洲精品在线观看| 各种免费的搞黄视频| 老司机亚洲免费影院| 国产午夜精品一二区理论片| 69精品国产乱码久久久| 97在线视频观看| 建设人人有责人人尽责人人享有的| 在线精品无人区一区二区三| 一本—道久久a久久精品蜜桃钙片| 黄色毛片三级朝国网站 | 99热这里只有精品一区| 亚洲电影在线观看av| 亚洲精品色激情综合| 99久国产av精品国产电影| 国产视频内射| 色94色欧美一区二区| 亚洲国产av新网站| 热99国产精品久久久久久7| 久热这里只有精品99| 欧美精品一区二区免费开放| 99九九在线精品视频 | 免费看光身美女| 3wmmmm亚洲av在线观看| 99久国产av精品国产电影| 中文字幕久久专区| 亚洲国产精品一区二区三区在线| 成人无遮挡网站| 嘟嘟电影网在线观看| 伊人久久精品亚洲午夜| 97超视频在线观看视频| 中文乱码字字幕精品一区二区三区| 在线观看一区二区三区激情| 男女啪啪激烈高潮av片| 亚洲精品久久午夜乱码| 国产精品熟女久久久久浪| 18禁动态无遮挡网站| 国产欧美日韩精品一区二区| 国产一区二区三区综合在线观看 | 80岁老熟妇乱子伦牲交| 成人18禁高潮啪啪吃奶动态图 | 成人特级av手机在线观看| 午夜免费鲁丝| 建设人人有责人人尽责人人享有的| 欧美最新免费一区二区三区| 观看美女的网站| 黑人巨大精品欧美一区二区蜜桃 | 69精品国产乱码久久久| 大香蕉97超碰在线| 欧美+日韩+精品| 精品国产一区二区久久| 日本黄色片子视频| 三级经典国产精品| 男女边吃奶边做爰视频| 亚洲不卡免费看| 日韩欧美 国产精品| 美女xxoo啪啪120秒动态图| 99久久人妻综合| 99国产精品免费福利视频| 国产成人91sexporn| av女优亚洲男人天堂| 九草在线视频观看| xxx大片免费视频| 亚洲欧美日韩另类电影网站| 国产黄片美女视频| 久久久国产一区二区| 午夜激情福利司机影院| 性高湖久久久久久久久免费观看| xxx大片免费视频| 久久精品久久久久久久性| 国产精品一区二区在线观看99| 欧美少妇被猛烈插入视频| 国产乱来视频区| 在线天堂最新版资源| 精品一区二区免费观看| 国产无遮挡羞羞视频在线观看| 全区人妻精品视频| 亚洲精品国产av成人精品| 国产精品久久久久久久电影| 插阴视频在线观看视频| 在线播放无遮挡| 视频区图区小说| 老司机影院成人| 日韩在线高清观看一区二区三区| 视频区图区小说| 欧美精品亚洲一区二区| 人体艺术视频欧美日本| 国产日韩一区二区三区精品不卡 | 欧美精品一区二区大全| 国产免费视频播放在线视频| 久久婷婷青草| 免费观看a级毛片全部| 精品卡一卡二卡四卡免费| 国产成人精品一,二区| 高清在线视频一区二区三区| 欧美老熟妇乱子伦牲交| 日韩精品有码人妻一区| 欧美性感艳星| 亚洲图色成人| 在线天堂最新版资源| 亚洲av男天堂| 日韩av在线免费看完整版不卡| 欧美xxⅹ黑人| 一级二级三级毛片免费看| 日本与韩国留学比较| 国产一区亚洲一区在线观看| 久久精品久久久久久久性| 免费观看在线日韩| 免费观看性生交大片5| 美女中出高潮动态图| 最新中文字幕久久久久| 韩国av在线不卡| 国产免费福利视频在线观看| 六月丁香七月| 在线精品无人区一区二区三| 大片电影免费在线观看免费| 久热久热在线精品观看| 色94色欧美一区二区| 亚州av有码| 18禁动态无遮挡网站| 你懂的网址亚洲精品在线观看| 日本av手机在线免费观看| 看免费成人av毛片| 99九九在线精品视频 | 亚洲欧美日韩卡通动漫| 免费av不卡在线播放| 亚洲欧美清纯卡通| 日韩人妻高清精品专区| 免费看av在线观看网站| 国产欧美日韩综合在线一区二区 | 日本色播在线视频| 内地一区二区视频在线| 十分钟在线观看高清视频www | 精品一区二区免费观看| 成人综合一区亚洲| 中国国产av一级| 国产成人午夜福利电影在线观看| 美女xxoo啪啪120秒动态图| 国产成人freesex在线| 精品少妇内射三级| 免费观看av网站的网址| 老司机影院成人| 日韩av不卡免费在线播放| 国产亚洲一区二区精品| 看非洲黑人一级黄片| 伊人久久精品亚洲午夜| av卡一久久| 肉色欧美久久久久久久蜜桃| 欧美日韩av久久| 欧美亚洲 丝袜 人妻 在线| 嘟嘟电影网在线观看| 一区二区av电影网| 成人综合一区亚洲| 美女xxoo啪啪120秒动态图| 久久久久国产网址| 久久久久国产精品人妻一区二区| 美女大奶头黄色视频| 天堂俺去俺来也www色官网| 久久人人爽人人爽人人片va| 亚洲精品乱久久久久久| 亚洲国产色片| 老司机影院毛片| 国产成人精品婷婷| 国产成人精品婷婷| 免费黄网站久久成人精品| 欧美精品一区二区免费开放| 国产淫片久久久久久久久| 夫妻性生交免费视频一级片| 亚洲欧美日韩卡通动漫| 久久久久久久久久成人| 亚洲精品国产色婷婷电影| 一个人免费看片子| 亚洲伊人久久精品综合| 国产精品蜜桃在线观看| 午夜av观看不卡| 麻豆精品久久久久久蜜桃| 两个人的视频大全免费| 一本一本综合久久| 免费观看av网站的网址| 99久久综合免费| 久久鲁丝午夜福利片| 人人妻人人添人人爽欧美一区卜| 国产精品不卡视频一区二区| 国产男女内射视频| 国产日韩欧美在线精品| 少妇人妻久久综合中文| 亚洲不卡免费看| 99热国产这里只有精品6| 国产精品久久久久久精品古装| 久久久午夜欧美精品| 国产老妇伦熟女老妇高清| 午夜福利在线观看免费完整高清在| 欧美日韩国产mv在线观看视频| 久久鲁丝午夜福利片| 97超视频在线观看视频| 免费久久久久久久精品成人欧美视频 | 狂野欧美激情性xxxx在线观看| 大又大粗又爽又黄少妇毛片口| 女人精品久久久久毛片| 十八禁高潮呻吟视频 | 精品亚洲乱码少妇综合久久| 亚洲av成人精品一二三区| 成人二区视频| 日产精品乱码卡一卡2卡三| 日韩一区二区三区影片| 日韩av在线免费看完整版不卡| 水蜜桃什么品种好| 精品人妻偷拍中文字幕| 中文乱码字字幕精品一区二区三区| 寂寞人妻少妇视频99o| 久久人人爽av亚洲精品天堂| 午夜免费男女啪啪视频观看| 国产亚洲最大av| 人人妻人人看人人澡| 亚洲欧洲精品一区二区精品久久久 | 国产 一区精品| videos熟女内射| 久久99热这里只频精品6学生| 免费看不卡的av| 美女福利国产在线| 高清黄色对白视频在线免费看 | 亚洲怡红院男人天堂| 精品一品国产午夜福利视频| 丰满迷人的少妇在线观看| 国产伦在线观看视频一区| 色婷婷av一区二区三区视频| 免费观看a级毛片全部| 一本大道久久a久久精品| 欧美老熟妇乱子伦牲交| 大片电影免费在线观看免费| 免费人成在线观看视频色| 精品一区二区免费观看| 亚洲国产精品国产精品| 亚洲综合色惰| 丝袜喷水一区| 亚洲自偷自拍三级| 成人综合一区亚洲| 午夜激情久久久久久久| 亚洲第一av免费看| 亚洲精品色激情综合| 中文字幕久久专区| 多毛熟女@视频| 蜜桃在线观看..| 欧美性感艳星| 99re6热这里在线精品视频| 中国国产av一级| 日本黄色日本黄色录像| 在线观看免费高清a一片| 老司机影院毛片| av不卡在线播放| 日韩欧美 国产精品| 午夜激情福利司机影院| 日本与韩国留学比较| 亚洲精品456在线播放app| av播播在线观看一区| 欧美 日韩 精品 国产| 国产伦精品一区二区三区视频9| 91成人精品电影| 一边亲一边摸免费视频| 亚洲精华国产精华液的使用体验| 青青草视频在线视频观看| 国产成人freesex在线| 王馨瑶露胸无遮挡在线观看| 在线看a的网站| 久久99热6这里只有精品| 日韩精品有码人妻一区| 老女人水多毛片| 亚洲自偷自拍三级| 亚洲欧美日韩东京热| 丝袜脚勾引网站| 国产欧美日韩综合在线一区二区 | 国产黄频视频在线观看| 午夜影院在线不卡| 国精品久久久久久国模美| 美女大奶头黄色视频| av线在线观看网站| 18禁裸乳无遮挡动漫免费视频| 精品一品国产午夜福利视频| 久久久久人妻精品一区果冻| 日韩,欧美,国产一区二区三区| 六月丁香七月| 一个人看视频在线观看www免费| 在线观看免费高清a一片| 中文字幕人妻丝袜制服| 国产精品99久久99久久久不卡 | 亚洲av二区三区四区| 中文字幕精品免费在线观看视频 | av黄色大香蕉| 国产亚洲91精品色在线| 久久毛片免费看一区二区三区| 久久综合国产亚洲精品| av有码第一页| 春色校园在线视频观看| 22中文网久久字幕| 国产精品国产三级国产专区5o| 久久精品国产自在天天线| 伊人久久国产一区二区| 国产国拍精品亚洲av在线观看| 欧美老熟妇乱子伦牲交| 中文字幕av电影在线播放| 最近中文字幕2019免费版| 男人爽女人下面视频在线观看| 久久99精品国语久久久| 晚上一个人看的免费电影| 亚洲伊人久久精品综合| 久久婷婷青草| 黄色怎么调成土黄色| 国产精品久久久久久久电影| 三级国产精品片| 久久久久久久久大av| 国产精品国产av在线观看| 久久97久久精品| 国产永久视频网站| 国产在线一区二区三区精| 亚洲av电影在线观看一区二区三区| 亚洲av日韩在线播放| 国产深夜福利视频在线观看| 偷拍熟女少妇极品色| 交换朋友夫妻互换小说| 涩涩av久久男人的天堂| 一区二区三区四区激情视频| 男人舔奶头视频| av在线老鸭窝| 亚洲欧美日韩卡通动漫| 国产精品国产三级专区第一集| 国产色婷婷99| av国产久精品久网站免费入址| 亚洲精品日韩av片在线观看| 亚洲人成网站在线观看播放|