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

    新型涂硼正比計數器研制及其實驗測試和模擬計算

    2015-05-16 02:17:26張紫霞魏志勇方美華陳國云
    原子能科學技術 2015年3期
    關鍵詞:圓柱型中子源靈敏

    張紫霞,魏志勇,朱 立,方美華,陳國云,強 鵬

    (1.南京航空航天大學材料科學與技術學院,江蘇南京 210016;2.南京航空航天大學航天學院,江蘇南京 210016;3.中國船舶重工集團公司第七零五研究所,云南昆明 650118;4.南昌大學理學院,江西南昌 330031)

    新型涂硼正比計數器研制及其實驗測試和模擬計算

    張紫霞1,魏志勇2,*,朱 立3,方美華2,陳國云4,強 鵬2

    (1.南京航空航天大學材料科學與技術學院,江蘇南京 210016;2.南京航空航天大學航天學院,江蘇南京 210016;3.中國船舶重工集團公司第七零五研究所,云南昆明 650118;4.南昌大學理學院,江西南昌 330031)

    本文將10B粉與1,2-二氯乙烷溶劑、芳華樹脂粘結劑混合研制了中子靈敏層10B膜,以此制作一新型涂硼正比計數器。為增加探測效率,除管內壁涂硼外,在管內增置了14片雙面涂硼環(huán)氧薄片,并用241Am-Be源測試了其性能。活度為3.7×109Bq的241Am-Be源實驗表明:此新型涂硼正比計數器坪長約150V,坪斜為8.2%/100V(750~900V);工作電壓為800V時測得新型涂硼正比計數器的計數率為50s-1,靈敏度為0.71cm2。與管內壁僅涂硼的正比計數器相比,新型涂硼正比計數器中子靈敏面積增加到3.15倍后,探測器坪長從80V增至150V,坪斜從12.4%/100V改進到7.58%/100V,靈敏度提高到2.63倍。同時,以基于蒙特卡羅方法的Geant4平臺,模擬計算單能中子及241Am-Be源中子照射包裹高密度聚乙烯慢化材料的新型涂硼正比計數器的響應和探測效率。241Am-Be源模擬結果和實驗結果吻合較好。實驗測試及模擬結果表明,涂硼技術與增加靈敏面積相結合的正比計數器的設計較成功。

    中子;涂硼正比計數器;坪曲線;探測效率

    常用的氣體中子探測器有三氟化硼正比計數器、3He正比計數器、硼電離室等[1-2]。3He正比計數器被廣泛用于檢測中子。3He正比管制備工藝復雜,同時因為核反恐活動的活躍致使3He氣體供需緊張,價格飛漲了幾十倍[3]。而BF3由于其對環(huán)境存在潛在的危害,氣體分解會產生含氟自由基,導致其在對環(huán)境要求嚴格的航空航天領域的使用受到了限制[4]?;谏鲜銮闆r,研制一種不使用3He和BF3材料的中子探測器十分必要。

    國內外已有對涂硼正比計數器的研究,如北京核儀器廠ZJ1521型涂硼正比計數器和GE Reuter Stokes的涂硼正比計數器[5]等。但10B電極制作技術已成為中子探測器研制的瓶頸,傳統(tǒng)的電泳技術不僅成本很大,且無法推廣。本文將富集10B粉、樹脂和1,2-二氯乙烷的混合液刷涂在探測器材料上,形成一層對熱中子敏感的硼膜。在常規(guī)的圓柱型正比計數器的基礎上,在計數器管內部放置14片雙面涂硼環(huán)氧薄片,以達到提高正比計數器靈敏度的目的。并在241Am-Be中子源輻射場中對設計研制的新型涂硼正比計數器進行性能測試。

    1 新型涂硼正比計數器結構設計

    正比計數器的應用非常廣泛,與電離室相同,它是利用非自持放電工作的氣體電離探測器。氣體電離產生的離子對在電場的作用下向探測器的兩極運動,輸出電離產生的電信號。正比計數器工作于坪曲線的正比區(qū),在該區(qū)離子收集對電離產生的離子對具有正比放大作用。

    1.1 陽極絲結構

    正比計數器大多由同軸圓柱型電極構成,圓柱型金屬管作為陰極,陽極通常采用一直徑很細的鎢絲。圓柱型正比計數器的電場分布非均勻。正比管計數器對氣體放大須具備足夠高的電場。一般,正比計數器的工作電壓約1kV,要實現倍增的最小工作電壓Vmin[6-7]為:

    其中:a為陽極絲半徑;b為正比管半徑;K為對應氣體倍增所需的最小約化場強;p為管內填充氣體壓強。從式(1)可見,陽極絲直徑是影響正比計數器工作電壓的主要因素之一。正比計數器直徑為40mm時,陽極絲直徑與最小工作電壓的關系如圖1所示。陽極絲直徑越小,所需工作電壓越小。鍍金鎢絲是公認的理想陽極材料。本文選用英國Goodfellow公司的直徑20μm的鍍金鎢絲。

    1.2 整體結構

    本文設計的新型涂硼正比計數器剖面結構如圖2所示。新型涂硼正比計數器由厚度為2mm的圓柱型不銹鋼(φ40mm×270mm)組成,內充4×104Pa的P10(90%Ar+10%CH4)氣體,管壁為陰極,中心放置一直徑20μm的鍍金鎢絲(陽極絲),并加高壓。裝好的陽極絲需滿足:1)始終繃緊;2)所受張力不使其斷裂。本工作涂硼正比管的陽極絲安裝選用錫焊的方法固定[8-9],選用掛重物的方法使陽極絲繃緊。同時注意20μm鍍金鎢絲的懸掛重物不超過100g。

    圖1 最小工作電壓與陽極絲直徑的關系Fig.1 Relationship between minimum working voltage and diameter of anode wire

    圖2 新型涂硼正比計數器剖面結構Fig.2 Section structure of new type boron-coating proportional counter

    圖3 新型涂硼正比計數器實驗裝置Fig.3 Experiment set-up diagram for new type boron-coating proportional counter

    研制中的另一關鍵部分是熱中子靈敏10B膜,此硼膜厚度[10]為0.8mg/cm2。硼膜的制作采用浸酯涂硼的方法。將同位素豐度為94.5%的10B粉與1,2-二氯乙烷溶劑、芳華樹脂黏結劑混合溶液涂抹到涂面上,放入90℃烤箱中保溫2h[11]。在不改變涂硼正比計數器結構及大小的情況下,將14片環(huán)形環(huán)氧薄片(厚度為0.2mm)按正比計數器的軸向排列,置于計數器管內,通過增加中子靈敏面積,達到提高中子探測效率的目的。環(huán)氧薄片兩面均涂硼膜。

    2 新型涂硼正比計數器的性能測試

    采用241Am-Be中子源對新型涂硼中子正比計數器進行輻照,實驗裝置如圖3所示。新型涂硼正比計數器是熱中子靈敏的探測器。241Am-Be源是快中子源,實現對快中子的探測,需在熱中子敏感探測器外包裹一層快中子慢化材料。本工作中選擇高密度聚乙烯作為慢化MeV中子的材料。在新型涂硼正比計數器外包裹厚為7.5mm的高密度聚乙烯(密度為0.95g/cm3),整個探測器及其包裹體的外觀是尺寸為φ55mm×250mm的圓柱。

    2.1 脈沖信號堆積分析

    新型涂硼正比計數器分別用南京大學輻射研究室提供的活度為3.7×109Bq的241Am-Be中子源和中國原子能科學研究院提供的活度為2.8×107Bq的241Am-Be中子源進行測量。實驗測量系統(tǒng)示于圖3。制備新型涂硼正比計數器時,采取以下措施提高探測器的信噪比:選擇屏蔽好的電纜線作為探測器內部輸出信號線;探測器結構設計中加入保護電極;使用鋁箔包裹信號引出口與BNC接口處(因接口處暴露在空氣中),避免周圍環(huán)境噪聲對探測器核脈沖信號的干擾。

    新型涂硼正比計數器工作電壓為800V時,3.7×109Bq的241Am-Be中子源測得的中子信號電壓脈沖如圖4a所示。探測器輸出脈沖幅度為3.68V,上升時間tR=980ns,脈沖寬度tW=2.5μs,脈沖的下降時間tF=1.52μs,新型涂硼正比計數器計數率的限值為4× 105s-1。由此實驗測試結果可看出,新型涂硼正比計數器抗噪聲性能較好,輸出信號的噪聲較小。

    若核輻射探測器記錄核輻射事件時的計數率為n,Δt時間內N個脈沖出現的概率p[4,12]可寫為:

    根據測量脈沖的結果,得到新型涂硼正比計數器計數率為105s-1時,對應的Δt為980ns時間內新型涂硼正比計數器的脈沖堆積概率為9.34%。大部分線性放大器的電路設計時均內置了脈沖堆積抑制器,以使探測器的脈沖堆積概率更小。這表明本工作研制的新型涂硼正比計數器可用于高注量率中子輻射場的檢測。

    同時將新型涂硼正比計數器放置到活度為2.8×107Bq的241Am-Be中子源輻射場中進行實驗測量,脈沖輸出波形如圖4b所示。從圖4b可直觀看出脈沖波形主要有兩個,分別是能量為1.47MeVα粒子和能量為0.84MeV Li離子的貢獻。α粒子能量較高,其相應的脈沖波形也較高。

    2.2 高壓坪曲線

    涂硼正比計數器的高壓坪曲線的測量使用南京大學輻射研究室的241Am-Be中子源。中子源置于圓柱形石蠟慢化體內。將涂硼正比計數器置于石蠟體上方進行輻照測量(圖3),源與涂硼正比計數器的距離為25cm。新型涂硼正比計數器的坪曲線測量結果如圖5a所示。圖5b示出了僅在管內壁涂硼的正比計數器的高壓坪曲線。表1列出了新型涂硼正比計數器和僅在管內壁涂硼的正比計數器的坪曲線參數對比。

    圖4 新型涂硼正比計數器輸出的中子脈沖Fig.4 Output neutron pulse signal of new type boron-coating proportional counter

    圖5 涂硼正比計數器高壓坪曲線Fig.5 High voltage plateau curve of boron-coating proportional counter

    從圖5和表1可看出,管內增加14片雙面涂硼環(huán)氧片使探測器的靈敏面積增加了311.3cm2。與僅在管內壁涂硼的探測器相比,中子靈敏面積增加到3.15倍。從實驗測得的坪曲線、靈敏度參數可看出增加中子靈敏面積對涂硼正比計數器的性能有顯著提高。

    坪曲線改善的原因是在管壁內放置了雙面涂硼膜的PCB板。相對于管壁內敷的硼膜,PCB板敷的雙面10B膜更靠近陽極絲,打到PCB板上硼膜的中子通過10B(n,α)7Li核反應產生的粒子更靠近陽極絲。核相互作用產生的帶電粒子越靠近陽極絲,電場越強,復合損失越小,電子在強場區(qū)的收集越好,亂真放電越少。同時,導致探測器靈敏度提高的原因之一是:14片環(huán)形雙面涂硼環(huán)氧片是按正比計數器的軸向排列的,即環(huán)氧片上的硼膜垂直于正比計數器管內壁的硼膜。中子經慢化體充分慢化后,中子的分布才各向同性。中子經0.75cm厚的高密度聚乙烯有限慢化后,中子角分布有一定傾向性。進入新型涂硼正比計數器的徑向中子數大于軸向中子數。相比軸向中子,徑向中子與PCB版上的硼膜發(fā)生10B(n,α)7Li核反應的概率較小,對應的中子靈敏面積無顯著增加。因此導致探測器靈敏度未隨靈敏面積成倍增加,而略有降低。

    表1 正比計數器增加靈敏面積前、后的性能參數對比Table 1 Comparison of proportional counter performance parameters before and after increased neutron sensitive area

    3 Geant4模擬研究

    實際中,通過實驗的方法獲得新型涂硼正比計數器從熱中子到幾十MeV區(qū)間的能量響應具有一定的局限性。通過實驗與計算相結合的方法,用實驗驗證模擬計算結果,模擬計算的結果為實驗提供參考。采用Geant4模擬中子在新型涂硼正比計數器中的輸運過程。中子數據庫選擇Geant4軟件包中基于美國ENDF/B-Ⅵ中子評價數據庫[13]的G4NeutronHP數據庫。

    包裹7.5mm厚高密度聚乙烯的新型涂硼正比計數器是圓柱型快中子組合探測器。定義探測器的探測效率ε為與10B發(fā)生反應的進入P10氣體內的帶電粒子數n和入射到圓柱型快中子組合探測器的中子數nin之比,即ε=n/nin。同時,若能量為E的中子入射到圓柱型快中子組合探測器的中子注量為Φ(E),對應的探測器計數為n(E),其能量響應R(E)=n(E)/Φ(E)。

    3.1 Am-Be源模擬結果

    基于Geant4平臺仿真模擬圖3所示的實驗測量條件下,中子在圓柱型快中子組合探測器中的中子輸運過程,得到探測器的探測效率和靈敏度。模擬中幾何結構參數為:涂硼正比計數器的直徑為40mm,有效長度為200mm,管內填充0.4×105Pa的P10氣體;硼膜的厚度為0.8g/cm2;0.2mm厚的環(huán)氧片的內徑為13mm,外徑為34mm;高密度聚乙烯密度為0.95g/cm3,外徑為55mm,厚7.5mm;Am-Be中子源到探測器中心的距離為25cm。

    表2給出了Am-Be中子源照射包裹7.5mm厚的同軸圓柱高密度聚乙烯慢化體的新型涂硼正比計數器的Geant4模擬結果和表1中的實驗結果的對比。從表2可看出,模擬計算得到的探測器的探測效率和靈敏度較實驗值小。由于實驗環(huán)境復雜,且未考慮實驗測量中后續(xù)電子學設備的影響,導致存在這些誤差。同時也證明了基于Geant4平臺編寫的模擬程序的可靠性較高。

    表2 模擬結果和實驗結果對比Table 2 Comparison between results of simulation and experiment

    3.2 不同能量中子垂直入射新型涂硼正比計數器的響應

    圓柱型快中子組合探測器對不同能量中子的響應和探測效率可通過實驗和模擬獲得。但實驗方法上,對不同能量中子源的獲取較繁瑣,且受實驗條件限制。用Geant4蒙特卡羅方法模擬計算得到新型涂硼正比計數器對不同能量中子的響應和探測效率。新型涂硼正比計數器外包裹7.5mm厚的同軸圓柱高密度聚乙烯慢化材料,平行中子束垂直入射,如圖6所示。在10-9~20MeV能量區(qū)間劃分40個單能中子能點。利用Geant4蒙特卡羅方法模擬計算得到新型涂硼正比計數器對不同入射能量中子的響應和對應的探測效率,結果示于圖7。

    由動量守恒和能量守恒可知,中子和10B反應產生的α粒子和7Li離子的能量分別為Eα=1.47MeV、ELi=0.84MeV。根據SRIM2008軟件,計算得到α粒子和7Li離子在P10氣體中的射程分別為Rα=3.55μm、RLi=1.84μm。α粒子的能量和射程均較7Li粒子的高。因此,從圖7可看出,α粒子在P10氣體中產生的脈沖數較7Li粒子的多;涂硼電離室加一層高密度聚乙烯慢化材料后,在低能中子段,由于高密度聚乙烯對中子的慢化,進入中子靈敏層的慢中子增多,新型涂硼正比計數器響應和探測效率增加。隨中子能量的增加,聚乙烯對高能中子的慢化作用減弱,高能中子與10B的反應截面減小,涂硼正比計數器的能量響應和探測效率減小。

    圖6 新型涂硼正比計數器對不同能量中子的響應模擬Fig.6 Simulation of different energy neutron response for new type boron-coating proportional counter

    圖7 Geant4模擬計算的新型涂硼正比計數器的能量響應和探測效率Fig.7 Energy response and detection efficiency of new type boron-coating proportional counter by Geant4

    4 結論

    本文選用浸酯涂硼方法自制硼膜,并采用管內放置14片雙面涂硼環(huán)氧片增加靈敏面積的方法對正比計數器進行結構改進設計。該新型涂硼正比計數器能在較高溫度下工作,具有環(huán)保性和經濟性,正比特性優(yōu)于BF3氣體[14],輸出信號有很好的抗γ性能。在相同測試環(huán)境、中子源、工作電壓(800V)的條件下,增加了靈敏面積的新型涂硼正比計數器與僅管壁內涂硼的正比計數器相比,坪區(qū)從80V(760~840V)擴展到150V(750~900V),坪斜從12.4%/100V改進到7.58%/100V,探測器的計數率從19s-1增加到50s-1,靈敏度提高到了2.63倍。在800V工作電壓下,測得探測器計數率為50s-1,靈敏度為0.71cm2。目前國內的涂硼正比計數器的水平是坪長為100V,坪斜為40%/100V;國外的水平達到坪長180V,坪斜18%/100V,靈敏度達1.6cm2[15]。故本工作制備的涂硼正比計數器已接近國際水平。

    同時,通過模擬計算得到了本文研制的新型涂硼正比計數器的響應和探測效率。模擬計算Am-Be源照射新型涂硼正比計數器的結果與實驗結果吻合較好。響應曲線的峰值點對應的中子入射能量為6×10-8MeV。對應的響應曲線給出了探測器對不同能量中子的響應。測量更高能量的中子,可通過在新型涂硼正比計數器外圍包裹更厚的慢化體實現。

    感謝南京大學和中國原子能科學研究院為實驗測量提供Am-Be源。

    [1] 方美華,魏志勇,張紫霞,等.涂硼正比計數器中熱中子響應函數理論計算及模擬驗證[J].南京航空航天大學學報,2012,44(1):81-86.

    FANG Meihua,WEI Zhiyong,ZHANG Zixia,et al.Theoretical calculation and simulation of thermal neutron response function of boron-lined proportional counter[J].Journal of Nanjing University of Aeronautics &Astronautics,2012,44(1):81-86(in Chinese).

    [2] 張紫霞,魏志勇,方美華,等.空間環(huán)境中子效應及測量技術[J].裝備環(huán)境工程,2009,6(4):5-11.

    ZHANG Zixia,WEI Zhiyong,FANG Meihua,et al.Neutron radiation effects of space environment and measuring technique[J].Equipment Environmental Engineering,2009,6(4):5-11(in Chinese).

    [3] 陳國云,魏志勇,方美華,等.正比計數器漏電流研究[J].南京航空航天大學學報,2010,42(5):670-674.

    CHEN Guoyun,WEI Zhiyong,FANG Meihua,et al.Leakage current in proportional counters[J].Journal of Nanjing University of Aeronautics &Astronautics,2010,42(5):670-674(in Chinese).

    [4] MAZED D.Experimental gas amplification study in boron-lined proportional counter for neutron detection[J].Radiation Measurement,2007,42(2):245-250.

    [5] KOUZES R T,ELY J H,ERIKSON L E,et al.Neutron detection alternatives to3He for national security applications[J].Nuclear Instruments and Methods in Physics Research A,2010,623:1 035-1 045.

    [6] MITEV K,SEGUR P,ALKAA A,et al.Study of non-equilibrium electron avalanches,application to proportioanl counters[J].Nucl Insutrm Methods A,2005,538(1-3):672-685.

    [7] KISER R W.Characteristic parameters of gastube proportinal counters[J].Appl Sci Res B,1960,8(1):183-200.

    [8] HIRANO H.A high-resolution cylindrical drift chamber for measuring CP violation at the KEKB asymmetric e+e-collider[R].Tokoy:KEK,2007.

    [9] 王學仁,郝秀蘭.利用改進的掛重法測定絲的張力[J].計量學報,1993,14(3):209-213.

    WANG Xueren,HAO Xiulan.Measuring wire tension with an improved method of hanging weight[J].Acta Metrologica Scinica,1993,14(3):209-213(in Chinese).

    [10]陳國云,魏志勇,辛勇,等.涂硼電離室中子探測效率和靈敏度[J].原子能科學技術,2011,45(2):244-249.

    CHEN Guoyun,WEI Zhiyong,XIN Yong,et al.Detection efficiency and sensitivity of boronlined ionization chamber of neutron[J].Atomic Energy Science and Technology,2011,45(2):244-249(in Chinese).

    [11]陳國云,魏志勇,黃福成,等.浸酯涂硼:中子探測器靈敏層制作方法[J].原子能科學技術,2013,47(3):485-489.

    CHEN Guoyun,WEI Zhiyong,HUANG Fucheng,et al.Resin-dipped boron coating:Produciton method for sensitive layer of neutron detectors[J].Atomic Energy Science and Technology,2013,47(3):485-489(in Chinese).

    [12]王經謹,范天明,錢永康.核電子學[M].北京:原子能出版社,1983:50-52.

    [13]Geant4user’s guide for tookit developers[R/OL].[S.l.]:[s.n.],2007.http:∥geant4.web.cern.ch/geant4/UserDocumentation/Users-Guides/ForToolkitDeveloper/html/index.html.

    [14]汲長松.中子探測實驗方法[M].北京:原子能出版社,1998:82-86.

    [15]DIGHE P M.New cathode design boron lined proportional counters for neutron area monitoring application[J].Nucl Instrum Methods A,2007,575(3):461-465.

    Development of New Type Boron-coating Proportional Counter and Its Experimental Investigation and Simulation Calculation

    ZHANG Zi-xia1,WEI Zhi-yong2,*,ZHU Li3,FANG Mei-h(huán)ua2,CHEN Guo-yun4,QIANG Peng2
    (1.College of Material Science and Technology,Nanjing University of Aeronautics &Astronautics,Nanjing210016,China;
    2.College of Astronautics,Nanjing University of Aeronautics &Astronautics,Nanjing210016,China;
    3.No.705 Research and Development Institute,China State Shipbuilding Corporation,Kunming650118,China;4.School of Science,Nanchang University,Nanchang330031,China)

    Three materials wherein suitable proportion of isotope enriched10B powder,1,2-ethylene dichloride and formvar resin were blended to make10B neutron sensitive coating by oneself.A new type proportional counter that coated with10B neutron sensi-tive coating was made.Furthermore,in order to increase the character and quality of the neutron detector,a set of 14annulus epoxy sheets which were sided with10B film were placed in the tube.A series performance tests were done by241Am-Be neutron source.The tests of 3.7×109Bq241Am-Be neutron source show that the plateau length of detector is 150Vfrom 750Vto 900V,while the plateau slope is 8.2%/100V.When the working voltage is 800V,the count rate of new type boron-coating proportional counter is 50s-1.The level of sensitivity is 0.71cm2.Compared with the detector only coated with10B film in the inner walls of detector,neutron sensitity area of the new detector increases to 3.15times.The results show that the plateau length increases from 80Vto 150V,and the plateau slope is improved from 12.4%/100Vto 7.58%/100V,while the neutron sensitivity increases to 2.63times.Using Geant4 software based on Monte Carlo method,this paper presented the response and detection efficiency of new type boron-coating proportional counter,which was covered with φ55mm×250mm cylinder high density polyethylene moderator material.The simulation results of Geant4are in agreement with the results of241Am-Be neutron source experiment.It shows the reliability of simulation application.

    neutron;boron-coating proportional counter;plateau curve;detector efficiency

    TL816.3

    A

    1000-6931(2015)03-0545-07

    10.7538/yzk.2015.49.03.0545

    2013-12-02;

    2014-04-23

    江蘇省普通高校研究生科研創(chuàng)新計劃資助項目(CXLX11_0192);中央高校基本科研業(yè)務費專項資金資助項目(NS2012070)

    張紫霞(1986—),女,江蘇武進人,博士研究生,從事空間輻射、粒子探測、輻射劑量與測量研究

    *通信作者:魏志勇,E-mail:wzy_msc@nuaa.edu.cn

    猜你喜歡
    圓柱型中子源靈敏
    “華龍一號”取消二次中子源研究
    板式換熱器進口處液體分布器的數值模擬*
    “超級顯微鏡”
    少兒科技(2020年3期)2020-05-13 13:15:25
    低壓靈敏放大器研究
    電子測試(2018年18期)2018-11-14 02:30:26
    我有一個靈敏的鼻子
    發(fā)明解讀刀具夾持裝置以及包括這種刀具夾持裝置的手持式電動工具
    電動工具(2018年4期)2018-08-22 02:02:28
    中國散裂中子源項目冷源系統(tǒng)設計
    圓柱型功能梯度雙材料界面裂紋問題研究
    吳靈敏:愛折騰的農場主
    海峽姐妹(2018年4期)2018-05-19 02:12:43
    中國散裂中子源首次打靶成功獲得中子束流
    春色校园在线视频观看| 国产人妻一区二区三区在| 全区人妻精品视频| 少妇猛男粗大的猛烈进出视频 | 少妇高潮的动态图| 午夜免费鲁丝| 亚洲av.av天堂| 黑人高潮一二区| 丰满乱子伦码专区| 亚洲精品久久午夜乱码| 欧美日韩国产mv在线观看视频 | 国内少妇人妻偷人精品xxx网站| 一本一本综合久久| 两个人的视频大全免费| 日韩亚洲欧美综合| 国产爱豆传媒在线观看| 亚洲人与动物交配视频| 亚洲精品乱码久久久久久按摩| 欧美xxxx性猛交bbbb| 午夜日本视频在线| 韩国av在线不卡| 亚洲成人精品中文字幕电影| 久久99热6这里只有精品| 久久国产乱子免费精品| 干丝袜人妻中文字幕| av国产免费在线观看| 午夜精品国产一区二区电影 | 大香蕉97超碰在线| 三级国产精品欧美在线观看| 国产淫片久久久久久久久| 男人舔奶头视频| 国产v大片淫在线免费观看| 日韩不卡一区二区三区视频在线| 久久精品国产亚洲av天美| 精品一区二区免费观看| 观看免费一级毛片| 777米奇影视久久| a级一级毛片免费在线观看| 99九九线精品视频在线观看视频| 在线精品无人区一区二区三 | 国产精品国产三级国产av玫瑰| 免费少妇av软件| a级毛片免费高清观看在线播放| 人妻一区二区av| 色播亚洲综合网| 成年女人看的毛片在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲精品日韩av片在线观看| 别揉我奶头 嗯啊视频| 18禁裸乳无遮挡动漫免费视频 | 欧美日韩精品成人综合77777| 中文在线观看免费www的网站| 一边亲一边摸免费视频| 国产伦理片在线播放av一区| 女人十人毛片免费观看3o分钟| 国产熟女欧美一区二区| 亚洲欧美一区二区三区国产| 国产精品无大码| 国产在线一区二区三区精| 日韩在线高清观看一区二区三区| 看十八女毛片水多多多| av在线亚洲专区| 国产色爽女视频免费观看| 日本黄色片子视频| av国产久精品久网站免费入址| 久久精品久久精品一区二区三区| av在线蜜桃| 国产欧美日韩一区二区三区在线 | 人妻夜夜爽99麻豆av| 欧美3d第一页| 亚洲av福利一区| 久久鲁丝午夜福利片| 中文乱码字字幕精品一区二区三区| av播播在线观看一区| 男女下面进入的视频免费午夜| 五月伊人婷婷丁香| 成人毛片a级毛片在线播放| 网址你懂的国产日韩在线| 国产成人午夜福利电影在线观看| 亚洲不卡免费看| 国产色婷婷99| videossex国产| 国产真实伦视频高清在线观看| 国产成人精品一,二区| 嫩草影院新地址| 国内少妇人妻偷人精品xxx网站| 丝瓜视频免费看黄片| 国产男女内射视频| 网址你懂的国产日韩在线| 国产中年淑女户外野战色| 美女国产视频在线观看| 欧美老熟妇乱子伦牲交| 亚洲av免费高清在线观看| 97人妻精品一区二区三区麻豆| av在线app专区| 视频中文字幕在线观看| 纵有疾风起免费观看全集完整版| 亚洲人成网站在线播| 国产视频内射| 七月丁香在线播放| 国产成人一区二区在线| 亚洲四区av| 免费播放大片免费观看视频在线观看| 黄色欧美视频在线观看| 欧美激情在线99| 午夜日本视频在线| 国产一区有黄有色的免费视频| 免费观看无遮挡的男女| 亚洲精品日韩在线中文字幕| 日日啪夜夜爽| 黄色怎么调成土黄色| 久久久亚洲精品成人影院| 精品人妻一区二区三区麻豆| 在线观看美女被高潮喷水网站| av在线蜜桃| videos熟女内射| 欧美激情在线99| 亚洲真实伦在线观看| 欧美97在线视频| 久久午夜福利片| 欧美日韩视频精品一区| 欧美变态另类bdsm刘玥| 亚洲最大成人手机在线| 久久久成人免费电影| 天天一区二区日本电影三级| 免费在线观看成人毛片| 国产成人精品婷婷| 美女脱内裤让男人舔精品视频| 成人亚洲精品av一区二区| 嫩草影院精品99| 亚洲精品aⅴ在线观看| 乱码一卡2卡4卡精品| 嫩草影院入口| 看非洲黑人一级黄片| av在线亚洲专区| 韩国av在线不卡| 久久精品综合一区二区三区| 三级国产精品欧美在线观看| 乱码一卡2卡4卡精品| 亚洲高清免费不卡视频| 99热国产这里只有精品6| 日韩成人伦理影院| 麻豆乱淫一区二区| 欧美另类一区| 国产乱人视频| 青春草视频在线免费观看| 免费观看的影片在线观看| 七月丁香在线播放| 精品国产一区二区三区久久久樱花 | 日本欧美国产在线视频| 国产成人91sexporn| 夫妻性生交免费视频一级片| 中文字幕免费在线视频6| 嫩草影院新地址| av国产免费在线观看| 国产精品99久久久久久久久| 建设人人有责人人尽责人人享有的 | 99久国产av精品国产电影| 一级毛片 在线播放| 波野结衣二区三区在线| 国产在线一区二区三区精| 欧美老熟妇乱子伦牲交| 内射极品少妇av片p| 天天一区二区日本电影三级| 日本熟妇午夜| 伦理电影大哥的女人| 网址你懂的国产日韩在线| 久久久a久久爽久久v久久| 男人和女人高潮做爰伦理| 日本猛色少妇xxxxx猛交久久| 在线观看av片永久免费下载| 精品久久久久久久人妻蜜臀av| 免费黄色在线免费观看| 99九九线精品视频在线观看视频| 国产精品无大码| 国产精品99久久99久久久不卡 | 日韩中字成人| 91久久精品国产一区二区成人| 亚洲精品国产av蜜桃| 老司机影院毛片| 免费电影在线观看免费观看| 欧美成人a在线观看| 久久久久久久国产电影| 国产毛片a区久久久久| 精品人妻视频免费看| 黄色怎么调成土黄色| 直男gayav资源| 伊人久久精品亚洲午夜| 一级毛片黄色毛片免费观看视频| 精品久久久久久久久亚洲| 搡老乐熟女国产| 亚洲精品日本国产第一区| 乱码一卡2卡4卡精品| 成人鲁丝片一二三区免费| 欧美人与善性xxx| 中文欧美无线码| 日本爱情动作片www.在线观看| 国产黄色视频一区二区在线观看| 精品视频人人做人人爽| 制服丝袜香蕉在线| 亚洲国产欧美在线一区| 91精品伊人久久大香线蕉| 别揉我奶头 嗯啊视频| 亚洲成人av在线免费| 亚洲怡红院男人天堂| 热99国产精品久久久久久7| 波多野结衣巨乳人妻| 国产午夜精品久久久久久一区二区三区| 午夜免费鲁丝| 青春草国产在线视频| 大香蕉久久网| 亚州av有码| 黄色视频在线播放观看不卡| 久久久久久久久久成人| 视频区图区小说| 一级毛片我不卡| 国产欧美日韩精品一区二区| 精品熟女少妇av免费看| 久久人人爽人人爽人人片va| 五月天丁香电影| 91精品一卡2卡3卡4卡| 国产亚洲5aaaaa淫片| 一级毛片久久久久久久久女| 亚洲人与动物交配视频| 小蜜桃在线观看免费完整版高清| 边亲边吃奶的免费视频| 日本av手机在线免费观看| 狠狠精品人妻久久久久久综合| 国产熟女欧美一区二区| 亚洲国产欧美人成| 免费在线观看成人毛片| 国产精品一区二区在线观看99| 在线观看三级黄色| 日韩 亚洲 欧美在线| 国产 一区 欧美 日韩| 麻豆久久精品国产亚洲av| 在线观看国产h片| 免费黄网站久久成人精品| 身体一侧抽搐| 久久影院123| 亚洲国产精品国产精品| 国产免费福利视频在线观看| 99九九线精品视频在线观看视频| 久久精品久久久久久久性| 色综合色国产| 亚洲国产精品专区欧美| 国产美女午夜福利| 国产毛片a区久久久久| 精品午夜福利在线看| 免费在线观看成人毛片| 少妇高潮的动态图| 精品午夜福利在线看| 久久精品夜色国产| 久久久久国产网址| 成人国产av品久久久| 国产精品一区二区性色av| 亚洲国产欧美在线一区| 国产免费一区二区三区四区乱码| av线在线观看网站| 亚洲电影在线观看av| 国产成人免费观看mmmm| 插阴视频在线观看视频| 亚洲国产成人一精品久久久| 色综合色国产| 亚洲国产av新网站| 亚洲国产最新在线播放| 亚洲精品乱码久久久v下载方式| 国产老妇女一区| 亚洲美女视频黄频| 国产亚洲精品久久久com| 在线亚洲精品国产二区图片欧美 | 免费播放大片免费观看视频在线观看| 亚洲婷婷狠狠爱综合网| 99久久精品一区二区三区| 亚洲av男天堂| 亚洲精品一区蜜桃| 又大又黄又爽视频免费| 一级av片app| 午夜福利高清视频| 亚洲va在线va天堂va国产| 久久久久国产网址| 69人妻影院| 亚洲精华国产精华液的使用体验| 欧美少妇被猛烈插入视频| 真实男女啪啪啪动态图| 高清日韩中文字幕在线| 亚洲精品乱码久久久v下载方式| 国产精品福利在线免费观看| 三级国产精品片| av在线蜜桃| 一级毛片我不卡| 国产成人精品福利久久| 亚洲精品久久午夜乱码| 99re6热这里在线精品视频| 国国产精品蜜臀av免费| 乱码一卡2卡4卡精品| 精品国产乱码久久久久久小说| a级毛色黄片| 中文天堂在线官网| 亚洲精品日本国产第一区| 国产老妇女一区| 国产免费视频播放在线视频| 特级一级黄色大片| 天美传媒精品一区二区| 内射极品少妇av片p| 激情五月婷婷亚洲| 国产在线男女| 国产 一区 欧美 日韩| 久久精品久久久久久噜噜老黄| 免费在线观看成人毛片| 午夜福利在线在线| 久久精品国产亚洲网站| 汤姆久久久久久久影院中文字幕| 久久精品国产亚洲网站| 精品熟女少妇av免费看| 大香蕉97超碰在线| 日本-黄色视频高清免费观看| 美女内射精品一级片tv| 波野结衣二区三区在线| 一区二区三区免费毛片| 99精国产麻豆久久婷婷| 嫩草影院入口| 三级男女做爰猛烈吃奶摸视频| 日日摸夜夜添夜夜爱| 国产探花极品一区二区| 91在线精品国自产拍蜜月| 亚洲人成网站在线观看播放| 色哟哟·www| 久久久久久久亚洲中文字幕| 久久久久久伊人网av| 欧美成人午夜免费资源| 亚洲国产av新网站| 精品人妻一区二区三区麻豆| 我的女老师完整版在线观看| 男女那种视频在线观看| 爱豆传媒免费全集在线观看| tube8黄色片| 美女国产视频在线观看| 乱码一卡2卡4卡精品| 日韩制服骚丝袜av| 日日啪夜夜爽| 日韩大片免费观看网站| 在线观看人妻少妇| 国产一区亚洲一区在线观看| 免费少妇av软件| 爱豆传媒免费全集在线观看| 国产免费又黄又爽又色| 国产男女内射视频| 少妇的逼水好多| 男女边吃奶边做爰视频| 午夜免费男女啪啪视频观看| 成人漫画全彩无遮挡| 亚洲人成网站在线观看播放| 我要看日韩黄色一级片| 亚洲无线观看免费| 午夜免费鲁丝| 国产成人精品福利久久| 美女内射精品一级片tv| 在线天堂最新版资源| 久久久久九九精品影院| 亚洲电影在线观看av| 亚洲欧美清纯卡通| 国产午夜福利久久久久久| 男插女下体视频免费在线播放| 欧美区成人在线视频| 性色avwww在线观看| 美女脱内裤让男人舔精品视频| 成人无遮挡网站| 一个人看的www免费观看视频| 毛片一级片免费看久久久久| 最近中文字幕2019免费版| 成人无遮挡网站| 欧美区成人在线视频| 日韩欧美精品免费久久| 天堂中文最新版在线下载 | 欧美激情久久久久久爽电影| 丝袜美腿在线中文| 欧美日韩综合久久久久久| 能在线免费看毛片的网站| 免费av不卡在线播放| 一区二区三区免费毛片| 日本-黄色视频高清免费观看| 大香蕉97超碰在线| 欧美精品人与动牲交sv欧美| 51国产日韩欧美| 久久精品人妻少妇| 三级国产精品片| 日本与韩国留学比较| 性色av一级| 日韩欧美精品v在线| 欧美丝袜亚洲另类| av国产免费在线观看| 中文字幕制服av| 精品国产乱码久久久久久小说| 欧美高清性xxxxhd video| 国产精品秋霞免费鲁丝片| 别揉我奶头 嗯啊视频| 美女国产视频在线观看| 国产一区有黄有色的免费视频| 精品午夜福利在线看| 日韩电影二区| 久久影院123| 久久久久久久久久成人| 91精品伊人久久大香线蕉| 老司机影院毛片| 久久久久九九精品影院| 91久久精品电影网| 国产乱来视频区| 中文欧美无线码| 欧美日韩亚洲高清精品| 久久久久精品久久久久真实原创| videos熟女内射| 精品午夜福利在线看| 国产欧美日韩精品一区二区| 国产精品精品国产色婷婷| 欧美激情国产日韩精品一区| 女人十人毛片免费观看3o分钟| 亚洲经典国产精华液单| 一本久久精品| 亚洲精品aⅴ在线观看| 最近的中文字幕免费完整| 久久99精品国语久久久| 国产91av在线免费观看| 亚洲精品乱久久久久久| 午夜视频国产福利| 嫩草影院新地址| 亚洲av二区三区四区| 日韩欧美精品免费久久| 中文天堂在线官网| 啦啦啦啦在线视频资源| 久久久久精品久久久久真实原创| 国产高清三级在线| a级毛片免费高清观看在线播放| 尾随美女入室| 久久韩国三级中文字幕| 制服丝袜香蕉在线| 久久久久性生活片| 在线观看一区二区三区| 久久久久久久亚洲中文字幕| 国产高清国产精品国产三级 | 日韩伦理黄色片| av在线app专区| 国产亚洲91精品色在线| 中文资源天堂在线| 国产成人免费无遮挡视频| 99精国产麻豆久久婷婷| 精品人妻熟女av久视频| 亚洲成人精品中文字幕电影| 亚洲精品一区蜜桃| 一级a做视频免费观看| 亚洲无线观看免费| 好男人视频免费观看在线| 国产免费福利视频在线观看| 一区二区三区乱码不卡18| 久久久久久久久久久免费av| 午夜老司机福利剧场| 国产乱人偷精品视频| 黑人高潮一二区| 少妇的逼水好多| 久久久久久久精品精品| 老司机影院毛片| 中文乱码字字幕精品一区二区三区| freevideosex欧美| 内地一区二区视频在线| 一级片'在线观看视频| 午夜免费男女啪啪视频观看| 国产av国产精品国产| 欧美区成人在线视频| 国产精品一区www在线观看| 中文欧美无线码| 99久久人妻综合| 在线观看三级黄色| 人人妻人人爽人人添夜夜欢视频 | 精品国产三级普通话版| 午夜精品一区二区三区免费看| 激情 狠狠 欧美| 成人免费观看视频高清| 麻豆成人午夜福利视频| 亚洲熟女精品中文字幕| 免费黄频网站在线观看国产| 国产高清国产精品国产三级 | 国产伦精品一区二区三区视频9| 黄色欧美视频在线观看| 视频中文字幕在线观看| 国产精品久久久久久久电影| 亚洲欧美日韩另类电影网站 | 国产国拍精品亚洲av在线观看| 亚洲国产精品999| 亚洲精品影视一区二区三区av| 真实男女啪啪啪动态图| 激情 狠狠 欧美| 网址你懂的国产日韩在线| 亚洲精品日本国产第一区| 国产精品成人在线| 极品少妇高潮喷水抽搐| 男女国产视频网站| 三级经典国产精品| 熟女人妻精品中文字幕| av播播在线观看一区| 又爽又黄a免费视频| 久久久精品94久久精品| 白带黄色成豆腐渣| 精品人妻熟女av久视频| 免费播放大片免费观看视频在线观看| 男人添女人高潮全过程视频| 嫩草影院新地址| 成年女人看的毛片在线观看| 欧美丝袜亚洲另类| 亚洲av.av天堂| 午夜免费男女啪啪视频观看| 大片免费播放器 马上看| 欧美精品人与动牲交sv欧美| 亚洲精品国产av蜜桃| 啦啦啦在线观看免费高清www| 亚洲性久久影院| 一本久久精品| 国产精品成人在线| tube8黄色片| 国产精品女同一区二区软件| 精品国产露脸久久av麻豆| 免费黄色在线免费观看| 少妇人妻 视频| 秋霞伦理黄片| 国产成人a∨麻豆精品| 两个人的视频大全免费| 中文字幕亚洲精品专区| 午夜激情久久久久久久| 黄色配什么色好看| 亚洲国产精品成人久久小说| 在线观看一区二区三区| a级一级毛片免费在线观看| 国产免费视频播放在线视频| 日韩人妻高清精品专区| 亚洲图色成人| a级毛片免费高清观看在线播放| 成年版毛片免费区| 简卡轻食公司| 国产精品av视频在线免费观看| 国产欧美日韩精品一区二区| av卡一久久| 久久久亚洲精品成人影院| 亚洲美女搞黄在线观看| 国产一区二区亚洲精品在线观看| 黄色一级大片看看| 国产成人aa在线观看| 草草在线视频免费看| 国产精品国产av在线观看| 国模一区二区三区四区视频| 久久99精品国语久久久| 99热这里只有是精品50| 色视频在线一区二区三区| 日本黄大片高清| 香蕉精品网在线| 久久人人爽人人爽人人片va| 亚洲,欧美,日韩| 亚洲精品国产av成人精品| 91aial.com中文字幕在线观看| 免费黄频网站在线观看国产| 午夜免费鲁丝| 纵有疾风起免费观看全集完整版| 熟妇人妻不卡中文字幕| 精品一区在线观看国产| 国产av码专区亚洲av| 啦啦啦啦在线视频资源| 久久人人爽人人爽人人片va| 成年版毛片免费区| 久久精品国产亚洲网站| 爱豆传媒免费全集在线观看| 国产一区二区三区综合在线观看 | 美女被艹到高潮喷水动态| 丝袜美腿在线中文| 在线看a的网站| 在线免费观看不下载黄p国产| 国产一区二区三区av在线| 九九爱精品视频在线观看| 国产一区亚洲一区在线观看| 午夜免费男女啪啪视频观看| 精品少妇久久久久久888优播| 黄色一级大片看看| 最近手机中文字幕大全| 男女啪啪激烈高潮av片| 国产视频内射| 精品国产三级普通话版| 99热这里只有精品一区| 国产精品一及| 久热这里只有精品99| 国产综合精华液| 伦理电影大哥的女人| 最近手机中文字幕大全| 亚洲精品aⅴ在线观看| 蜜桃亚洲精品一区二区三区| 五月玫瑰六月丁香| 亚洲色图av天堂| 国产中年淑女户外野战色| 免费不卡的大黄色大毛片视频在线观看| 十八禁网站网址无遮挡 | 欧美成人一区二区免费高清观看| 黄片无遮挡物在线观看| 成人欧美大片| 一区二区三区乱码不卡18| 日本午夜av视频| 赤兔流量卡办理| 亚洲综合精品二区| 日韩中字成人| 在线播放无遮挡| 久久久久网色| 亚洲国产av新网站| 内射极品少妇av片p| 国产美女午夜福利| 97超视频在线观看视频| 黄色一级大片看看| 三级国产精品片| 亚洲内射少妇av| 国产熟女欧美一区二区| 国产欧美日韩一区二区三区在线 |