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

    基于白光中子源的核石墨硼當(dāng)量測量

    2020-11-25 13:25:40王小鶴胡繼峰陳金根蔡翔舟王納秀王宏偉韓建龍
    原子能科學(xué)技術(shù) 2020年11期
    關(guān)鍵詞:中子源中子當(dāng)量

    王小鶴,胡繼峰,陳金根,蔡翔舟,王納秀,王宏偉,韓建龍,*

    (1.中國科學(xué)院 上海應(yīng)用物理研究所,上海 201800;2.中國科學(xué)院 先進(jìn)核能創(chuàng)新研究院,上海 201800;3.中國科學(xué)院大學(xué),北京 100049;4.中國科學(xué)院 上海高等研究院,上海 201204)

    核石墨具有耐高溫、熱中子散射截面大、吸收截面小等特性,被廣泛用作慢化劑或反射層材料,目前全球已有百余座核反應(yīng)堆使用核石墨作中子慢化劑和反射層[1]。中國科學(xué)院上海應(yīng)用物理研究所承擔(dān)的釷基熔鹽堆項目[2]中也采用石墨作為慢化劑、反射層材料。雖然石墨的熱中子吸收截面很小,但核石墨在反應(yīng)堆中用量巨大,故截面對反應(yīng)堆的設(shè)計和安全性具有重大影響。由于原料產(chǎn)地、生產(chǎn)工藝的不同,使得石墨中含有不同的雜質(zhì),尤其是Li、B等熱中子吸收截面大的雜質(zhì),它們吸收中子后會引起堆芯反應(yīng)性變化[3-4]。因此,反應(yīng)堆設(shè)計對核石墨硼當(dāng)量的數(shù)值提出限制性要求,材料中雜質(zhì)元素對熱中子的吸收以相當(dāng)于硼的量來表示,稱為熱中子吸收的硼當(dāng)量[5],如美國熔鹽堆MSRE要求石墨硼當(dāng)量<80 ppm[6],高溫氣冷堆HTR-10要求石墨硼當(dāng)量<1.3 ppm[7],在釷基熔鹽實驗堆中,也要求石墨的硼當(dāng)量<2 ppm。硼當(dāng)量是反映核材料純度的重要指標(biāo)之一。

    硼當(dāng)量的測量方法分為兩種,一種是元素分析法,如電感耦合等離子體光譜法(ICP-AES)[8]、電感耦合等離子體質(zhì)譜法(ICP-MS)[9]、激光剝蝕-電感耦合等離子體質(zhì)譜法(LA-ICP-MS)[10]、輝光放電質(zhì)譜法(GD-MS)[11]等,這些方法必須測定雜質(zhì)的含量,根據(jù)硼當(dāng)量因子即可計算出硼當(dāng)量,所需時間較長,存在取樣、溶樣、進(jìn)樣過程中引入不確定性等問題,且核反應(yīng)堆中石墨用量多、體積大,元素分析法難以體現(xiàn)大塊材料的平均硼當(dāng)量;另一種是熱中子宏觀吸收截面法,此方法進(jìn)行硼當(dāng)量測量需強(qiáng)中子源及復(fù)雜的探測器裝置,相關(guān)研究較少。1950年,Mireille初次使用Ra-Be中子源進(jìn)行了石墨宏觀吸收截面的測量實驗[12],當(dāng)測試樣品體積為0.13 m3、質(zhì)量為200 kg時,將宏觀吸收截面轉(zhuǎn)化為硼當(dāng)量后相對不確定度約為8%,但樣品體積太過巨大,無法用于日常測試。2005年,Bolewski等使用252Cf中子源進(jìn)行了核石墨硼當(dāng)量測量[13],當(dāng)樣品體積為25 cm3時,在僅考慮統(tǒng)計誤差的前提下,其測量不確定約為3.1 ppm,如果考慮其他誤差來源,實際不確定度遠(yuǎn)高于此值。這些實驗均采用同位素中子源,源強(qiáng)低且易受散射效應(yīng)影響,導(dǎo)致測量精度低、時間長。

    利用白光中子源(photoneutron source,PNS)進(jìn)行硼當(dāng)量測量具有較大優(yōu)勢:首先白光中子源的強(qiáng)度較高,有益于測量精度及測量效率的提高;其次,采用透射法進(jìn)行測量,可單獨(dú)處理樣品,降低實驗中不確定度來源,還可避免具有輻射、腐蝕等特性的材料對實驗裝置造成污染;白光中子源可獲得與熱堆較為接近的熱譜,中子能譜峰值能量約0.06 eV[14],可得到在反應(yīng)堆環(huán)境下的等效硼當(dāng)量。本文基于15 MeV電子加速器驅(qū)動的白光中子源[15],建立利用透射法進(jìn)行核石墨硼當(dāng)量測量的理論方法,并開展核石墨硼當(dāng)量的測量研究。

    1 理論介紹

    白光中子源產(chǎn)生平行中子束,基于中子透射率TR來推導(dǎo)硼當(dāng)量的計算公式[16]:

    (1)

    其中:I0為入射中子束強(qiáng)度;I為經(jīng)過厚度為x的靶后未碰撞的中子束強(qiáng)度;σ為靶核的微觀截面;N為單位體積內(nèi)的原子核數(shù);Σ為宏觀截面。當(dāng)樣品含有雜質(zhì)時,Σ=ΣC+ΣB,ΣC為基底石墨靶核本身的宏觀截面(散射和吸收截面),ΣB為雜質(zhì)元素的宏觀截面,包含散射和吸收截面。雜質(zhì)元素含量很低,其核子數(shù)密度NB遠(yuǎn)低于被測石墨的核子數(shù)密度NC,可認(rèn)為ΣB主要來自吸收截面,即ΣB=σBNB。因此,有:

    (2)

    假設(shè)第1個樣品為純樣品(不包含雜質(zhì)),即ΣB=0。當(dāng)測量第n個幾何形狀、材質(zhì)相同的樣品時,若NB?NC,可認(rèn)為各樣品的密度相同,基體樣品本身的總截面相同,即ΣC不變。ln TR的變化由式(2)中第2項引起。將兩個樣品的ln TR相減即可得到第n個樣品中的硼當(dāng)量BE:

    (3)

    其中:A為刻度參數(shù),A=σBxNAρ/MB;NA為阿伏伽德羅常數(shù);ρ為材料密度。根據(jù)上式,測量刻度樣品的中子透射率、厚度、密度,即可得到測量系統(tǒng)的刻度參數(shù),建立測量硼當(dāng)量理論方法。

    當(dāng)測量樣品的密度、厚度與刻度樣品一致時,只需獲得測量樣品的Δln TR*,即可得到測量樣品的硼當(dāng)量BE*:

    (4)

    中子源、樣品、探測器的物理參數(shù)均對核石墨硼當(dāng)量測量有重要影響,但白光中子源的布局及運(yùn)行參數(shù)相對固定,因此僅利用MCNP模擬分析樣品的厚度、直徑、密度、均勻度的影響。模擬計算采用5塊刻度樣品,硼當(dāng)量分別為0.11、42.68、85.26、170.29、340.13 ppm,厚度均為1.3 cm。其中,硼當(dāng)量為0.11 ppm的石墨樣品作為基體樣品進(jìn)行相對扣除,其他為刻度樣品。模擬計算采用CENDL-TMSR-V1數(shù)據(jù)庫[17],使用分段銜接與減方差結(jié)合的方法提高計算的可靠性[15],計算中考慮了石墨熱散射效應(yīng)。

    2 模擬分析

    圖1示出了樣品厚度為13、15 mm(a)以及樣品直徑為 4、5、6 cm(b)時的模擬結(jié)果,擬合度均大于0.999。結(jié)果顯示,隨著厚度增加,刻度參數(shù)也在逐步增加,樣品厚度每增加1 mm,刻度參數(shù)增大約8%,影響較大。與直徑5 cm的刻度參數(shù)相比,直徑為4 cm時刻度參數(shù)偏低約10.9%,導(dǎo)致硼當(dāng)量被高估,直徑6 cm時兩者刻度參數(shù)基本一致,相對偏差為0.37%。因此,當(dāng)樣品直徑大于或等于束流直徑(5 cm)時,對測量結(jié)果的影響較小。

    圖2示出了樣品密度分別為1.80、1.82、1.84 g/cm3時模擬計算的結(jié)果。結(jié)果顯示,刻度參數(shù)分別為7.49×10-5、7.64×10-5、7.76×10-5,擬合度均大于0.999。隨著密度的增加,刻度參數(shù)也會增加,密度每變化1%,可導(dǎo)致刻度參數(shù)增大約1%,可見密度的影響也較大。

    對于元素分析法,樣品的不均勻性會導(dǎo)致測量結(jié)果與實際結(jié)果差異較大,基于白光中子源的硼當(dāng)量測量方法也需分析樣品不均勻性的影響。采用上述刻度樣品,其中硼當(dāng)量為170.29 ppm的樣品采用4種分布模式。該樣品分為4個區(qū)域(圖3),每個區(qū)域?qū)?yīng)1個硼當(dāng)量,每個區(qū)域中硼分布均勻,樣品總硼當(dāng)量不變。第1種分布方式為均勻分布,即每個區(qū)域的硼當(dāng)量均為170.29 ppm,其他3種分布方式各區(qū)域的硼當(dāng)量與第1種分布方式的硼當(dāng)量的倍數(shù)關(guān)系列于表1,計算結(jié)果示于圖4。結(jié)果顯示,A分別為7.54×10-5、7.60×10-5、7.53×10-5、7.58×10-5,擬合度均大于0.998。與第1種分布方式相比,其他3種分布的相對偏差分別為0.8%、0.1%、0.5%,可見樣品的不均勻性對該方法的影響遠(yuǎn)小于對元素分析法的影響。

    圖1 樣品厚度(a)、直徑(b)對刻度參數(shù)的影響Fig.1 Influence of sample thickness (a) and diameter (b) on A

    圖2 樣品密度對刻度參數(shù)的影響Fig.2 Influence of sample density on A

    圖3 樣品分區(qū)示意圖Fig.3 Sketch map of sample subarea

    表1 樣品中各區(qū)域的硼當(dāng)量分布Table 1 Boron equivalent distribution in sample

    3 實驗測量

    3.1 白光中子源

    TMSR白光中子源為電子直線加速器驅(qū)動的中子源[14],電子能量最大達(dá)20 MeV,設(shè)計最大平均電子束功率為1.5 kW。慢化前初級中子能譜的峰值能量約1 MeV,慢化后可得到連續(xù)的熱能譜。

    圖4 樣品均勻度對刻度參數(shù)的影響Fig.4 Influence of sample homogeneity on A

    實驗采用6LiF/ZnS(Ag)探測器,型號為EJ426HD2,直徑50 mm,厚度0.48 mm。數(shù)據(jù)獲取系統(tǒng)采用CAEN公司生產(chǎn)的波形采樣器DT5720以及多功能數(shù)據(jù)獲取軟件[18],可實現(xiàn)硬件控制、數(shù)據(jù)獲取、波形在線顯示等功能,并與實驗靶控制軟件耦合,實現(xiàn)靶位自動定時更換、多參數(shù)獲取、自動報警等功能。

    3.2 樣品制備

    為獲得滿足要求的核石墨刻度樣品,利用核級石墨粉末與硼粉混合并進(jìn)行高溫?zé)Y(jié)的方法來制備刻度樣品。采用高精度的電子天平(精度為0.01 mg)進(jìn)行稱量,利用行星式球磨機(jī)按比例多次混合,每次混合2 h,以保證樣品中硼分布相對均勻。利用真空熱壓爐對其進(jìn)行熱壓成片,5 h內(nèi)升至1 850 ℃,保溫保壓5 h,然后自然冷卻。根據(jù)配比率,可得出樣品中硼當(dāng)量,所得刻度樣品參數(shù)列于表2。

    表2 石墨刻度樣品參數(shù)Table 2 Parameter of graphite calibration sample

    3.3 實驗測量

    實驗分兩輪進(jìn)行,共測試5塊樣品。每輪測量5個靶,包括3個石墨刻度樣品、本底靶以及空靶。其中,N靶為純石墨樣品,其余兩個樣品靶為硼當(dāng)量相對較大的刻度樣品,本底靶為12.8 cm的含硼5%的聚乙烯塊[19]。第1輪實驗的刻度樣品依次為I、N、G,第2輪實驗的刻度樣品依次為J、N、H。每輪實驗的電子束流打靶功率均為250 W,電子能量為16.5 MeV。每個靶測量時間為5 min,每個循環(huán)共25 min。

    4 結(jié)果討論

    4.1 束流穩(wěn)定性監(jiān)測

    在束流60°方向設(shè)置監(jiān)測探測器,用來監(jiān)測束流的穩(wěn)定性及中子注量率的變化。圖5示出了兩輪實驗中空靶時監(jiān)測探測器與TOF探測器每5 min內(nèi)的中子計數(shù)比。圖5顯示,大部分監(jiān)測探測器和TOF探測器的中子計數(shù)比是穩(wěn)定的,出現(xiàn)異常的位置是因為加速器出現(xiàn)掉束或真空異常的情況,這也表明監(jiān)測探測器的中子計數(shù)可體現(xiàn)束流的變化,可用于對TOF探測器中子計數(shù)的歸一。

    圖5 空靶時監(jiān)測探測器與TOF探測器中子計數(shù)比Fig.5 Ratio of neutron countings for monitoring detector and TOF detector in blank target

    4.2 數(shù)據(jù)檢驗

    根據(jù)戈羅貝斯舍棄標(biāo)準(zhǔn)對監(jiān)測探測器獲取數(shù)據(jù)進(jìn)行檢驗[20],將個別測量值與平均值相差較大的數(shù)據(jù)去除。兩輪實驗分別生成115、168組數(shù)據(jù)。經(jīng)過檢驗后分別選擇了91、144組數(shù)據(jù),所有樣品的TOF探測器的中子計數(shù)示于圖6。

    圖6 TOF探測器中子計數(shù)Fig.6 Neutron counting for each target obtained by TOF detector

    4.3 數(shù)據(jù)修正

    硼當(dāng)量測量針對的是雜質(zhì)對熱中子的吸收,可根據(jù)指定能區(qū)內(nèi)的總中子計數(shù)變化體現(xiàn)雜質(zhì)對熱中子的吸收。結(jié)合白光中子源的能譜及測量精度[15],對比了4×10-9~1×10-5MeV之間總中子計數(shù)變化。圖7示出了各刻度樣品5 min內(nèi)原始總中子計數(shù)。

    圖7 各測量樣品5 min內(nèi)中子計數(shù)Fig.7 Neutron counting for each sample in 5 min

    對各刻度樣品的數(shù)據(jù)進(jìn)行本底扣除,并利用監(jiān)測探測器的中子計數(shù)進(jìn)行歸一,以降低束流穩(wěn)定性帶來的影響,計算公式如下:

    (5)

    其中:In、I0、IB分別為刻度樣品、空靶、本底靶的中子計數(shù);Mn、M0分別為刻度樣品、空靶的監(jiān)測探測器中子計數(shù)。各樣品修正后歸一到5 min內(nèi)的中子計數(shù)列于表3。

    表3 修正后歸一到5 min內(nèi)各樣品的中子計數(shù)Table 3 Neutron counting for each sample in 5 min after modification

    模擬結(jié)果顯示樣品的厚度、密度對測量影響較大,而所用刻度樣品的厚度與密度也難以保持一致,需對其進(jìn)行修正。式(2)可變?yōu)榕c厚度、密度相關(guān)的公式:

    (6)

    將公式兩邊同時除x即可得到單位厚度下的ln TR。假設(shè)純樣品密度為ρ0,第n塊純樣品密度為在測量條件一致時,根據(jù)上式可獲得對應(yīng)的用于不同密度下的基體效應(yīng)扣除,式(3)可表示為:

    (7)

    根據(jù)式(7)修正后的計算結(jié)果示于圖8。結(jié)果顯示,修正前,Δln TR明顯偏高;修正后,Δln TR明顯降低,且基本呈線性。

    4.4 實驗測量結(jié)果

    圖9示出了石墨硼當(dāng)量刻度實驗的測量結(jié)果。結(jié)果顯示,實驗測量與模擬計算的模擬刻度參數(shù)分別為6.35×10-5、7.61×10-5,擬合度均大于0.998,兩者的相對偏差為19.8%。

    差異主要來自兩個方面:1) 模擬計算中,束斑的中子呈高斯分布,但實驗中,其分布并非是嚴(yán)格的高斯分布,且由于束流的不穩(wěn)定也會產(chǎn)生輕微的變化,由此會導(dǎo)致實際測量中子數(shù)低于模擬計算值,從而對穿透率產(chǎn)生影響;2) 在熱區(qū)峰值處的實際測量的熱中子數(shù)低于模擬計算值,且隨著能量的增加,探測器探測效率急劇下降,也會對結(jié)果產(chǎn)生影響。整體而言,實驗測量結(jié)果與理論計算結(jié)果趨勢一致,根據(jù)該刻度曲線,即可進(jìn)行測試樣品的硼當(dāng)量測量。

    圖8 修正前、后結(jié)果對比Fig.8 Comparison of results with and without modification

    圖9 實驗結(jié)果與MCNP模擬結(jié)果對比Fig.9 Comparison of measurement and simulation results

    4.5 不確定度分析

    基于白光中子源利用透射法進(jìn)行硼當(dāng)量測量實驗中所涉及的測量量為中子計數(shù)、樣品厚度、密度等,在進(jìn)行不確定度分析時,除考慮中子計數(shù)的統(tǒng)計誤差,還需考慮樣品厚度、密度、樣品均勻度、束流穩(wěn)定性帶來的影響。

    根據(jù)式(7),給出Δln TR的相對不確定度計算公式:

    (8)

    厚度的不確定度主要來自多次測量及測量儀器本身的誤差,密度的不確定度主要來自厚度、直徑及質(zhì)量的測量誤差。經(jīng)過計算,N、J、I、H、G靶的厚度相對不確定度分別為0.53%、1.47%、1.04%、1.17%、1.25%,密度的相對不確定度分別為0.57%、1.49%、1.15%、1.19%、1.4%。

    對TR而言,相對不確定度主要來自中子統(tǒng)計誤差、本底以及束流穩(wěn)定性,表達(dá)式如下:

    (9)

    經(jīng)計算,TR相對不確定度均約0.21%。

    根據(jù)式(8)~(9),可給出各刻度樣品Δln TR的相對不確定度為2.3%、1.76%、1.86%、2.05%。當(dāng)刻度樣品的硼當(dāng)量一定時,刻度參數(shù)的相對不確定度為2.1%。考慮雜質(zhì)分布不均勻帶來的影響,刻度參數(shù)的相對不確定度為2.9%。

    根據(jù)式(4),可給出測量樣品的相對不確定度計算公式為:

    (10)

    假設(shè)第1塊石墨刻度樣品為測量樣品,根據(jù)硼當(dāng)量刻度參數(shù),計算出硼當(dāng)量為41.78 ppm,相對不確定度為3.01%,高于同等條件下已有實驗的硼當(dāng)量測量精度[13]。

    5 總結(jié)

    本文根據(jù)透射法原理,建立了基于電子加速器驅(qū)動的白光中子源進(jìn)行核石墨硼當(dāng)量測量的實驗方法。利用5個刻度樣品完成刻度實驗,通過對實驗數(shù)據(jù)的一系列修正,建立硼當(dāng)量相關(guān)的刻度曲線,通過該刻度曲線及測量樣品的中子計數(shù)即可獲得測量樣品的硼當(dāng)量。計算結(jié)果顯示,基于該方法獲得測量樣品硼當(dāng)量的相對不確定度為3.01%,高于同條件下已有實驗的硼當(dāng)量測量精度,可用于工程中對核石墨的硼當(dāng)量檢測。

    猜你喜歡
    中子源中子當(dāng)量
    “華龍一號”取消二次中子源研究
    “超級顯微鏡”
    少兒科技(2020年3期)2020-05-13 13:15:25
    3D打印抗中子輻照鋼研究取得新進(jìn)展
    中國散裂中子源項目冷源系統(tǒng)設(shè)計
    中國散裂中子源首次打靶成功獲得中子束流
    基于PLC控制的中子束窗更換維護(hù)系統(tǒng)開發(fā)與研究
    黃河之聲(2016年24期)2016-02-03 09:01:52
    DORT 程序進(jìn)行RPV 中子注量率計算的可靠性驗證
    超壓測試方法對炸藥TNT當(dāng)量計算結(jié)果的影響
    中子深度定量分析的相對分析法
    計算物理(2014年2期)2014-03-11 17:01:27
    亚洲情色 制服丝袜| 黄色一级大片看看| 成人亚洲精品一区在线观看| 久久久久久久久久久免费av| 国产成人精品无人区| 22中文网久久字幕| 日本午夜av视频| 亚洲精品乱久久久久久| 国产男人的电影天堂91| 综合色丁香网| 99久国产av精品国产电影| 日韩欧美精品免费久久| 亚洲欧美中文字幕日韩二区| tube8黄色片| videosex国产| 日本午夜av视频| 欧美日韩视频高清一区二区三区二| 在线免费观看不下载黄p国产| 国产免费又黄又爽又色| 少妇高潮的动态图| 亚洲天堂av无毛| 黑人猛操日本美女一级片| 男人爽女人下面视频在线观看| 一个人免费看片子| 欧美日韩综合久久久久久| 久久精品久久久久久噜噜老黄| xxxhd国产人妻xxx| 激情视频va一区二区三区| 国产成人午夜福利电影在线观看| 一级片免费观看大全| 亚洲内射少妇av| 五月天丁香电影| 国产一区有黄有色的免费视频| 成人亚洲欧美一区二区av| 日韩精品有码人妻一区| 91在线精品国自产拍蜜月| 国产精品久久久久久精品古装| 欧美激情 高清一区二区三区| 亚洲精品视频女| 黑人猛操日本美女一级片| 97人妻天天添夜夜摸| 亚洲国产看品久久| 男女免费视频国产| 777米奇影视久久| videossex国产| 亚洲伊人色综图| 久久精品人人爽人人爽视色| 久久久久精品久久久久真实原创| 中文字幕av电影在线播放| 如何舔出高潮| 久久精品国产亚洲av涩爱| 久久99蜜桃精品久久| 亚洲国产精品一区三区| 天堂俺去俺来也www色官网| 熟女av电影| 欧美日韩成人在线一区二区| 久久午夜综合久久蜜桃| 一边摸一边做爽爽视频免费| 考比视频在线观看| 男女边摸边吃奶| 国产乱人偷精品视频| 伦理电影大哥的女人| 国产极品粉嫩免费观看在线| 一二三四中文在线观看免费高清| 国产欧美日韩一区二区三区在线| 国产一区二区激情短视频 | 国产精品偷伦视频观看了| 日本av免费视频播放| 99热这里只有是精品在线观看| 亚洲四区av| 午夜激情av网站| 韩国av在线不卡| 久久久精品区二区三区| 自线自在国产av| 国产又爽黄色视频| 久久久久久久久久成人| 亚洲中文av在线| 一级片免费观看大全| 一区二区av电影网| 九色成人免费人妻av| 免费观看无遮挡的男女| a级毛片黄视频| 欧美日韩亚洲高清精品| 亚洲欧美色中文字幕在线| 宅男免费午夜| 夫妻性生交免费视频一级片| 亚洲国产日韩一区二区| 狠狠婷婷综合久久久久久88av| 肉色欧美久久久久久久蜜桃| 久久久国产精品麻豆| 中文天堂在线官网| 国产国语露脸激情在线看| 精品一区二区免费观看| 激情五月婷婷亚洲| 久久综合国产亚洲精品| 性色avwww在线观看| 我要看黄色一级片免费的| 国产亚洲一区二区精品| 999精品在线视频| 高清毛片免费看| 午夜免费观看性视频| 赤兔流量卡办理| 亚洲av电影在线进入| 亚洲av国产av综合av卡| 久久热在线av| 国产精品蜜桃在线观看| 欧美激情极品国产一区二区三区 | 狂野欧美激情性bbbbbb| 日韩在线高清观看一区二区三区| 日本av免费视频播放| 久久国产亚洲av麻豆专区| 精品99又大又爽又粗少妇毛片| 免费高清在线观看日韩| 亚洲成国产人片在线观看| 亚洲国产色片| 久久久欧美国产精品| 国产一区二区三区av在线| 免费观看性生交大片5| 黑人巨大精品欧美一区二区蜜桃 | 97精品久久久久久久久久精品| 久久亚洲国产成人精品v| 日本午夜av视频| 亚洲婷婷狠狠爱综合网| 亚洲国产毛片av蜜桃av| 亚洲成人av在线免费| 国产一区二区三区av在线| 亚洲内射少妇av| 久久热在线av| 天美传媒精品一区二区| 午夜免费鲁丝| 久久精品国产综合久久久 | 精品人妻一区二区三区麻豆| 亚洲伊人色综图| 蜜桃在线观看..| 免费av不卡在线播放| 精品国产露脸久久av麻豆| 男女啪啪激烈高潮av片| 女人精品久久久久毛片| 午夜免费鲁丝| 人妻少妇偷人精品九色| 99久久人妻综合| 久久精品国产亚洲av涩爱| 天天躁夜夜躁狠狠久久av| www.色视频.com| 校园人妻丝袜中文字幕| 水蜜桃什么品种好| 日韩制服丝袜自拍偷拍| 熟女av电影| 久久久久久久久久久久大奶| 亚洲成国产人片在线观看| 最近2019中文字幕mv第一页| 99re6热这里在线精品视频| 亚洲国产精品专区欧美| 国产老妇伦熟女老妇高清| 国产精品不卡视频一区二区| 成人国产麻豆网| 在现免费观看毛片| 日韩欧美一区视频在线观看| 哪个播放器可以免费观看大片| 91在线精品国自产拍蜜月| 在线观看www视频免费| 亚洲国产毛片av蜜桃av| 国产免费又黄又爽又色| 中国美白少妇内射xxxbb| www.av在线官网国产| 91精品三级在线观看| 成人亚洲精品一区在线观看| 边亲边吃奶的免费视频| 日韩不卡一区二区三区视频在线| 久久热在线av| 91精品伊人久久大香线蕉| 久久人人爽人人片av| 成年人免费黄色播放视频| 七月丁香在线播放| 欧美精品亚洲一区二区| 色婷婷久久久亚洲欧美| 熟女av电影| 亚洲伊人久久精品综合| a 毛片基地| 十八禁网站网址无遮挡| 一级毛片 在线播放| 97精品久久久久久久久久精品| 大香蕉久久网| 国产成人精品婷婷| 久久99热6这里只有精品| 免费在线观看完整版高清| 国产免费福利视频在线观看| 午夜日本视频在线| 免费久久久久久久精品成人欧美视频 | 黄色视频在线播放观看不卡| 在线天堂中文资源库| 精品午夜福利在线看| 丰满饥渴人妻一区二区三| 亚洲成人av在线免费| 多毛熟女@视频| 韩国av在线不卡| 国产乱人偷精品视频| 国产在线视频一区二区| 国产不卡av网站在线观看| 亚洲经典国产精华液单| 18在线观看网站| 成人国产av品久久久| 免费女性裸体啪啪无遮挡网站| 国产福利在线免费观看视频| 99热国产这里只有精品6| 亚洲情色 制服丝袜| 少妇人妻精品综合一区二区| 国产av一区二区精品久久| 中文天堂在线官网| 国产亚洲最大av| 日韩制服丝袜自拍偷拍| 亚洲精品自拍成人| 丰满少妇做爰视频| 日韩在线高清观看一区二区三区| av片东京热男人的天堂| 一边摸一边做爽爽视频免费| 色94色欧美一区二区| 大码成人一级视频| 亚洲精品国产av成人精品| 欧美成人精品欧美一级黄| 一区二区日韩欧美中文字幕 | 午夜福利在线观看免费完整高清在| 22中文网久久字幕| 少妇人妻 视频| 欧美少妇被猛烈插入视频| 最近手机中文字幕大全| 涩涩av久久男人的天堂| 免费看光身美女| 中文字幕另类日韩欧美亚洲嫩草| 日韩不卡一区二区三区视频在线| 亚洲精品av麻豆狂野| 人妻少妇偷人精品九色| 国产免费又黄又爽又色| 天天操日日干夜夜撸| 国产综合精华液| 国产一区有黄有色的免费视频| 日本色播在线视频| 丝袜脚勾引网站| 色94色欧美一区二区| 亚洲欧洲精品一区二区精品久久久 | 亚洲av成人精品一二三区| 国产成人aa在线观看| 夫妻性生交免费视频一级片| 精品第一国产精品| 国产色婷婷99| 少妇人妻 视频| 国产乱来视频区| 久久女婷五月综合色啪小说| 亚洲精品乱久久久久久| 国产老妇伦熟女老妇高清| 涩涩av久久男人的天堂| 在线天堂中文资源库| 在线看a的网站| 国产精品人妻久久久影院| 一二三四中文在线观看免费高清| 少妇熟女欧美另类| 男女无遮挡免费网站观看| 免费观看在线日韩| 18禁裸乳无遮挡动漫免费视频| 午夜老司机福利剧场| 久久久久久久国产电影| 日韩成人伦理影院| 久久久久久久久久久久大奶| 一级毛片电影观看| 精品一区二区三区视频在线| 免费少妇av软件| 亚洲精品乱码久久久久久按摩| 人人妻人人添人人爽欧美一区卜| 久久精品熟女亚洲av麻豆精品| 精品一区二区免费观看| 婷婷色av中文字幕| 久久久久久久久久久久大奶| 最近手机中文字幕大全| 精品酒店卫生间| 在线观看国产h片| 欧美精品高潮呻吟av久久| 少妇被粗大猛烈的视频| 人妻系列 视频| 成年美女黄网站色视频大全免费| 一区二区三区四区激情视频| 国产精品国产av在线观看| 成人国产av品久久久| 久久精品久久精品一区二区三区| 久久久欧美国产精品| 五月玫瑰六月丁香| 老熟女久久久| 一区在线观看完整版| 亚洲五月色婷婷综合| 精品熟女少妇av免费看| 亚洲成人手机| 亚洲国产日韩一区二区| 国产日韩一区二区三区精品不卡| 久久精品熟女亚洲av麻豆精品| 一级毛片 在线播放| 中文字幕免费在线视频6| 18在线观看网站| 美女国产高潮福利片在线看| 欧美日韩一区二区视频在线观看视频在线| 久久韩国三级中文字幕| 在线精品无人区一区二区三| av天堂久久9| 久久狼人影院| 少妇的丰满在线观看| 国产深夜福利视频在线观看| 欧美日韩精品成人综合77777| 少妇熟女欧美另类| 亚洲av男天堂| 韩国高清视频一区二区三区| 一级,二级,三级黄色视频| 大香蕉97超碰在线| 国产精品一区二区在线观看99| 亚洲精品久久午夜乱码| 久久av网站| 亚洲精品国产色婷婷电影| 美女中出高潮动态图| 一本色道久久久久久精品综合| 蜜桃国产av成人99| 日产精品乱码卡一卡2卡三| a 毛片基地| 久久精品人人爽人人爽视色| 国产一区亚洲一区在线观看| 欧美精品亚洲一区二区| 日韩欧美一区视频在线观看| av播播在线观看一区| a 毛片基地| 插逼视频在线观看| 亚洲,一卡二卡三卡| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国国产精品蜜臀av免费| 国产免费又黄又爽又色| 黑人高潮一二区| 男女午夜视频在线观看 | 亚洲丝袜综合中文字幕| 亚洲综合色网址| 在线观看人妻少妇| 国产精品秋霞免费鲁丝片| 在线亚洲精品国产二区图片欧美| 如日韩欧美国产精品一区二区三区| 少妇 在线观看| 伊人亚洲综合成人网| 啦啦啦视频在线资源免费观看| 国产xxxxx性猛交| 精品亚洲成国产av| 人体艺术视频欧美日本| 久久人人爽av亚洲精品天堂| 午夜影院在线不卡| 最新中文字幕久久久久| 51国产日韩欧美| 一级毛片黄色毛片免费观看视频| 日本av手机在线免费观看| av视频免费观看在线观看| 99热网站在线观看| 亚洲欧美一区二区三区国产| 亚洲国产日韩一区二区| 国产精品三级大全| 熟妇人妻不卡中文字幕| 国产精品 国内视频| 99久国产av精品国产电影| 国产成人精品福利久久| 午夜免费鲁丝| 日韩在线高清观看一区二区三区| 一边亲一边摸免费视频| av在线app专区| 国产女主播在线喷水免费视频网站| av国产久精品久网站免费入址| 尾随美女入室| 又粗又硬又长又爽又黄的视频| 免费黄网站久久成人精品| 国产 一区精品| 国产成人aa在线观看| 欧美日韩综合久久久久久| 综合色丁香网| av视频免费观看在线观看| 日本与韩国留学比较| 亚洲四区av| 国产成人一区二区在线| xxx大片免费视频| 女人被躁到高潮嗷嗷叫费观| 成人毛片60女人毛片免费| 久久午夜福利片| 日韩成人伦理影院| 国国产精品蜜臀av免费| 高清在线视频一区二区三区| 亚洲av中文av极速乱| 亚洲高清免费不卡视频| 午夜视频国产福利| 99热网站在线观看| 国产精品久久久久久久久免| 观看av在线不卡| 18禁观看日本| 成人国产麻豆网| 亚洲国产精品国产精品| 亚洲成人av在线免费| 亚洲综合精品二区| 寂寞人妻少妇视频99o| 亚洲国产精品成人久久小说| 91午夜精品亚洲一区二区三区| 亚洲欧美日韩另类电影网站| 蜜臀久久99精品久久宅男| 久久精品久久久久久久性| 免费观看a级毛片全部| 女的被弄到高潮叫床怎么办| 亚洲av成人精品一二三区| 日日啪夜夜爽| 精品午夜福利在线看| 午夜福利乱码中文字幕| 久久ye,这里只有精品| 欧美激情 高清一区二区三区| 久久人人爽av亚洲精品天堂| 成年人免费黄色播放视频| av卡一久久| 在线天堂最新版资源| 99久国产av精品国产电影| 中文天堂在线官网| 麻豆精品久久久久久蜜桃| 黄片无遮挡物在线观看| 久久久久久久大尺度免费视频| 久久这里只有精品19| 一个人免费看片子| 男女边吃奶边做爰视频| 国产毛片在线视频| 一二三四中文在线观看免费高清| 精品国产一区二区三区久久久樱花| 国产亚洲精品久久久com| 中文天堂在线官网| 久久99热6这里只有精品| 日韩中字成人| 乱码一卡2卡4卡精品| 午夜激情久久久久久久| 午夜免费男女啪啪视频观看| 免费大片18禁| 亚洲经典国产精华液单| 2021少妇久久久久久久久久久| 亚洲av欧美aⅴ国产| 欧美成人午夜精品| 少妇被粗大猛烈的视频| 在线观看人妻少妇| 色哟哟·www| 国产毛片在线视频| 毛片一级片免费看久久久久| 欧美xxⅹ黑人| 韩国av在线不卡| 哪个播放器可以免费观看大片| 黑人巨大精品欧美一区二区蜜桃 | 久久久久久人妻| 桃花免费在线播放| 日韩一区二区视频免费看| 9色porny在线观看| 少妇高潮的动态图| 国产亚洲精品久久久com| 日韩一本色道免费dvd| 婷婷色av中文字幕| 日日撸夜夜添| 成人午夜精彩视频在线观看| 国产精品女同一区二区软件| av有码第一页| 免费日韩欧美在线观看| 视频在线观看一区二区三区| 久久久久久久大尺度免费视频| 七月丁香在线播放| 久久久精品免费免费高清| 日韩成人av中文字幕在线观看| 日韩人妻精品一区2区三区| 老司机影院成人| 这个男人来自地球电影免费观看 | 国产精品一二三区在线看| 国产欧美日韩综合在线一区二区| 亚洲第一av免费看| 黑人欧美特级aaaaaa片| 色网站视频免费| 肉色欧美久久久久久久蜜桃| 秋霞在线观看毛片| 少妇人妻 视频| 亚洲第一av免费看| 黄片播放在线免费| 99精国产麻豆久久婷婷| 久久精品夜色国产| 视频区图区小说| 亚洲国产毛片av蜜桃av| 精品第一国产精品| 99精国产麻豆久久婷婷| 久久久久国产精品人妻一区二区| 国产极品粉嫩免费观看在线| 欧美97在线视频| 免费大片18禁| 国产av码专区亚洲av| 亚洲人与动物交配视频| 黄片无遮挡物在线观看| 国产亚洲av片在线观看秒播厂| 一级毛片电影观看| 久久精品人人爽人人爽视色| 午夜激情av网站| 五月伊人婷婷丁香| 在线观看一区二区三区激情| 在线观看美女被高潮喷水网站| 大陆偷拍与自拍| 最近2019中文字幕mv第一页| 免费播放大片免费观看视频在线观看| 国产淫语在线视频| 亚洲,一卡二卡三卡| 国产免费福利视频在线观看| 国产无遮挡羞羞视频在线观看| 免费观看a级毛片全部| av在线老鸭窝| 免费大片18禁| 免费观看av网站的网址| 欧美丝袜亚洲另类| 爱豆传媒免费全集在线观看| 欧美人与善性xxx| 精品久久久精品久久久| 国产精品 国内视频| 狠狠精品人妻久久久久久综合| 国产成人精品婷婷| 91精品三级在线观看| 一本大道久久a久久精品| 久久99精品国语久久久| 精品人妻一区二区三区麻豆| 久久久久人妻精品一区果冻| 人体艺术视频欧美日本| 建设人人有责人人尽责人人享有的| 国产成人91sexporn| 热99久久久久精品小说推荐| 久久精品国产综合久久久 | 日韩制服丝袜自拍偷拍| 一级毛片电影观看| 国产男女超爽视频在线观看| 秋霞在线观看毛片| av免费在线看不卡| 免费在线观看完整版高清| 一本久久精品| av播播在线观看一区| 久久久a久久爽久久v久久| 多毛熟女@视频| 国产成人一区二区在线| 亚洲欧美日韩卡通动漫| 亚洲,欧美精品.| 久久这里有精品视频免费| 久久人人爽av亚洲精品天堂| 久久鲁丝午夜福利片| 免费观看无遮挡的男女| 久久女婷五月综合色啪小说| 99久久精品国产国产毛片| 1024视频免费在线观看| 欧美日韩亚洲高清精品| 亚洲伊人久久精品综合| 日韩精品免费视频一区二区三区 | 亚洲av免费高清在线观看| 免费观看av网站的网址| 热99国产精品久久久久久7| 人体艺术视频欧美日本| 我要看黄色一级片免费的| 欧美精品一区二区大全| 自线自在国产av| 国产片内射在线| 99国产综合亚洲精品| 在线观看www视频免费| 五月玫瑰六月丁香| 欧美激情极品国产一区二区三区 | 美女主播在线视频| 久久精品aⅴ一区二区三区四区 | 一二三四在线观看免费中文在 | 晚上一个人看的免费电影| 亚洲内射少妇av| 色94色欧美一区二区| 国产日韩欧美视频二区| 黑人高潮一二区| 国产欧美亚洲国产| 最近的中文字幕免费完整| 国产成人免费观看mmmm| 久久久久久久久久成人| 9热在线视频观看99| 亚洲国产av新网站| 欧美老熟妇乱子伦牲交| 久久国内精品自在自线图片| 欧美日韩视频精品一区| 日韩一区二区视频免费看| 久久午夜福利片| 成人手机av| 亚洲国产看品久久| 亚洲高清免费不卡视频| 久久99蜜桃精品久久| 日韩视频在线欧美| 热re99久久国产66热| 亚洲精品,欧美精品| 国产又爽黄色视频| av卡一久久| 91国产中文字幕| 美女中出高潮动态图| 王馨瑶露胸无遮挡在线观看| 亚洲美女黄色视频免费看| 精品国产一区二区三区四区第35| 在线观看免费高清a一片| 汤姆久久久久久久影院中文字幕| 久久久a久久爽久久v久久| 亚洲,一卡二卡三卡| 亚洲精品日韩在线中文字幕| 欧美精品高潮呻吟av久久| 国产成人免费无遮挡视频| 超色免费av| 亚洲三级黄色毛片| 精品国产国语对白av| 午夜av观看不卡| 考比视频在线观看| 2021少妇久久久久久久久久久| 免费播放大片免费观看视频在线观看| a级片在线免费高清观看视频| 欧美 亚洲 国产 日韩一| 少妇人妻 视频| 欧美+日韩+精品| 最近2019中文字幕mv第一页| 亚洲成国产人片在线观看| 性高湖久久久久久久久免费观看| 三级国产精品片| 2018国产大陆天天弄谢| 99国产综合亚洲精品|