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

    90o伴隨粒子法D-D中子產(chǎn)額測量修正因子計(jì)算

    2017-02-09 06:14:27李建一王俊潤黃智武盧小龍徐大鵬馬占文姚澤恩
    核技術(shù) 2017年1期
    關(guān)鍵詞:測量

    李建一 王俊潤,2 張 宇,2 黃智武 盧小龍,2 徐大鵬,2 張 杰 韋 崢,2 馬占文 姚澤恩,2

    1(蘭州大學(xué) 核科學(xué)與技術(shù)學(xué)院 蘭州 730000)

    2(蘭州大學(xué) 教育部中子應(yīng)用技術(shù)工程研究中心 蘭州 730000)

    90o伴隨粒子法D-D中子產(chǎn)額測量修正因子計(jì)算

    李建一1王俊潤1,2張 宇1,2黃智武1盧小龍1,2徐大鵬1,2張 杰1韋 崢1,2馬占文1姚澤恩1,2

    1(蘭州大學(xué) 核科學(xué)與技術(shù)學(xué)院 蘭州 730000)

    2(蘭州大學(xué) 教育部中子應(yīng)用技術(shù)工程研究中心 蘭州 730000)

    開展了90o伴隨粒子法D-D中子產(chǎn)額測量中的修正因子理論計(jì)算方法研究,基于MATLAB軟件平臺,開發(fā)了用于修正因子計(jì)算的計(jì)算機(jī)程序。計(jì)算給出了厚靶條件下,入射氘能量在20-700 keV范圍,90o伴隨粒子法D-D中子產(chǎn)額測量各向異性修正因子Rthick、中子和質(zhì)子產(chǎn)額比(Yd,n/Yd,p)thick及總修正因子RY,并與早先的研究結(jié)果進(jìn)行了對比,分析了計(jì)算結(jié)果的不確定度,總修正因子計(jì)算數(shù)據(jù)的不確定度約為2%。

    伴隨粒子法,D-D中子產(chǎn)額,中子發(fā)生器,修正因子

    緊湊型D-D中子發(fā)生器是一種可產(chǎn)生約2.5 MeV快中子的小型化加速器單能快中子源,與D-T中子發(fā)生器相比,不使用放射性氚靶,更具安全性,且存在氘離子自注入效應(yīng)、靶壽命長、運(yùn)行成本低、可代替同位素中子源等優(yōu)點(diǎn),在工業(yè)在線中子活化分析、中子照相、爆炸物及毒品檢測等方面有廣泛的應(yīng)用價(jià)值[1-2]。早先的緊湊型D-D中子發(fā)生器主要以密封中子管為主,密封中子管的優(yōu)點(diǎn)是尺寸很小,其缺點(diǎn)是使用微型潘寧離子源,D束流只有幾十微安,D-D中子產(chǎn)額低,且是一次性真空密封設(shè)備,壽命偏短,使用成本高。開發(fā)D束流大于1 mA,D-D中子產(chǎn)額大于108n·s-1,元件可更換的緊湊型長壽命D-D中子發(fā)生器已成為中子應(yīng)用技術(shù)開發(fā)的迫切需要。國際上一些重要實(shí)驗(yàn)室已發(fā)展了各種形式的緊湊型中子發(fā)生器[3-5],國內(nèi)此方面的研究相對滯后[6-7]。

    在蘭州大學(xué),一臺緊湊型D-D中子發(fā)生器正在研制中。其中,緊湊型D-D中子發(fā)生器的快中子產(chǎn)額的準(zhǔn)確測量是需解決的關(guān)鍵問題。伴隨粒子法是實(shí)現(xiàn)D-D中子發(fā)生器快中子產(chǎn)額絕對測量的手段之一,在低能D束下,D-D聚變反應(yīng)存在以下兩個反應(yīng)道:D(d,n)3He反應(yīng)和D(d,p)T反應(yīng)。其中,第一個反應(yīng)道產(chǎn)生快中子n和伴隨粒子3He,第二個反應(yīng)道產(chǎn)生質(zhì)子(p)和氚(T)粒子,因3He能量偏低(約0.78 MeV),難于準(zhǔn)確測量,一般通過測量第二個反應(yīng)道產(chǎn)生的能量較高的伴隨質(zhì)子(能量約3.05 MeV)來實(shí)現(xiàn)D-D中子產(chǎn)額測量[8]。

    在采用伴隨質(zhì)子法實(shí)現(xiàn)D-D中子產(chǎn)額測量中,有兩個參數(shù)需要通過理論計(jì)算來確定,即中子和質(zhì)子的產(chǎn)額比和出射中子的各向異性修正因子。早先的工作大多針對135o伴隨粒子法開展了修正因子計(jì)算[9-11],90o伴隨粒子法的相關(guān)研究較少,而且早先的工作采用了較早的D-D反應(yīng)截面數(shù)據(jù),計(jì)算結(jié)果存在較大的不確定度。本文針對緊湊型D-D中子發(fā)生器90o伴隨粒子法中子產(chǎn)額測量的需求,開展相關(guān)理論研究和修正因子計(jì)算方法研究,計(jì)算給出入射氘離子能量在20-700 keV范圍,90o伴隨粒子法D-D中子產(chǎn)額測量所需修正因子數(shù)據(jù),為緊湊型D-D中子發(fā)生器中子產(chǎn)額的測量提供基礎(chǔ)數(shù)據(jù)。數(shù)據(jù)也可用于其他類型D-D中子發(fā)生器伴隨粒子法D-D中子產(chǎn)額測量。

    1 緊湊型D-D中子發(fā)生器結(jié)構(gòu)及90o伴隨粒子法中子產(chǎn)額測量原理

    如圖1所示,緊湊型D-D中子發(fā)生器由離子源、離子束引出加速系統(tǒng)、靶系統(tǒng)、高壓接頭、真空腔和真空泵等部分裝配而成。離子束引出加速系統(tǒng)可將離子源引出D離子加速到100-120 keV,轟擊在D靶上,發(fā)生D-D聚變反應(yīng),產(chǎn)生能量約2.5 MeV的快中子。中子發(fā)生器結(jié)構(gòu)為圓筒形,長約98 cm,最大法蘭外徑為?28 cm。在90o方向有一管道用于安裝金硅面壘探測器,通過測量D-D反應(yīng)產(chǎn)生的伴隨質(zhì)子來實(shí)現(xiàn)D-D中子產(chǎn)額測量。

    圖1 緊湊型高產(chǎn)額D-D中子發(fā)生器結(jié)構(gòu)示意圖1. 離子源,2. 引出加速電極,3. 靶,4. 金硅面壘探測器,5. 靶冷卻管道接頭,6. 不銹鋼外殼,7. 絕緣環(huán),8. 冷卻管道,9. 高壓電纜,10. 真空抽氣口Fig.1 Schematic of compact high-yield D-D neutron generator.1. Ion source, 2. Accelerating electrode, 3. Target, 4. Au-Si surface barrier detector, 5. Cooling joints of target, 6. Stainless steel shell, 7. Insulating ring, 8. Cooling duct, 9. High voltage cables, 10. Vacuum suction

    在厚靶條件下,D-D中子產(chǎn)額可用式(1)描述[8]:

    式中:np為探測器測到的伴隨質(zhì)子計(jì)數(shù)率;Ωp為探測器對靶點(diǎn)所張立體角;Rthick為各向異性修正因子;(Yd,n/Yd,p)為D-D反應(yīng)中子與質(zhì)子的積分產(chǎn)額比。各向異性修正因子和中子與質(zhì)子的積分產(chǎn)額比可分別由式(2)、(3)描述[9-10],即:

    式中:Ed為入射氘束能量;σd,n(E)和σd,p(E)分別為實(shí)驗(yàn)室系中D(d,n)3He和D(d,p)T反應(yīng)的積分截面;(dσd,p/dω')(E,θ)為質(zhì)心系中D(d,p)T反應(yīng)對應(yīng)實(shí)驗(yàn)室系θ方向出射質(zhì)子的微分截面;(dE/dx)(E)為氘離子在氘鈦(DTi)靶中的阻止本領(lǐng);(dω'/dω)(E,θ)為質(zhì)心系和實(shí)驗(yàn)室系之間的立體角轉(zhuǎn)換因子。

    綜上所述,將Rthick和(Yd,n/Yd,p)thick的乘積記為總的修正因子RY,即:

    則D-D反應(yīng)的中子產(chǎn)額可表達(dá)為:

    由式(5)可見,要實(shí)現(xiàn)伴隨粒子法D-D反應(yīng)中子產(chǎn)額的絕對測量,必須根據(jù)D-D反應(yīng)截面數(shù)據(jù)計(jì)算出各向異性修正因子及中子與質(zhì)子的積分產(chǎn)額比。

    2 計(jì)算方法

    針對緊湊型D-D中子發(fā)生器90o伴隨粒子法中子產(chǎn)額測量,本研究將采用薄靶近似方法計(jì)算各向異性修正因子和中子與質(zhì)子的積分產(chǎn)額比,即將厚靶分割成很多足夠薄的薄層,每一薄層可被近似看作一個薄靶,則由式(2)、(3),各向異性修正因子和中子與質(zhì)子的積分產(chǎn)額比可近似轉(zhuǎn)化為下列求和形式,即:

    式中:Ei為氘束入射到第i層時的能量;ΔEi為氘離子在第i層中損失的能量。設(shè)E0=Ed,ΔEi=0,Ei和ΔEi可由式(8)、(9)計(jì)算,即:

    式中:Δxi-1為被分割的靶中第i-1薄層厚度。

    式中:m1、m2、m3、m4分別為射氘離子、靶核氘、出射質(zhì)子和氚的質(zhì)量;Q為D(d,p)T反應(yīng)的反應(yīng)能(Q=4.033 MeV);E'ip為質(zhì)心系中的質(zhì)子能量,可由式(13)-(15)計(jì)算[12]:

    式中:θ=90o為實(shí)驗(yàn)室系中質(zhì)子出射角;θc為質(zhì)心系中質(zhì)子出射角;Ei為入射氘離子能量;Eip為實(shí)驗(yàn)室系中θ方向的質(zhì)子能量,可由式(16)給出[13]:

    2.1 截面與阻止本領(lǐng)數(shù)據(jù)

    質(zhì)心系中,(dσd,p/dω')(Ei,θ)為D(d,p)T反應(yīng)對應(yīng)實(shí)驗(yàn)室系θ方向出射質(zhì)子的微分截面,可根據(jù)式(17)、(18)進(jìn)行計(jì)算[14-15],即:

    式中:a、b、c和A、B為常數(shù);dσd,p(90o)/dω'為質(zhì)心系中θc=90o的D(d,p)T反應(yīng)微分截面,常數(shù)和質(zhì)心系微分截面數(shù)據(jù)見文獻(xiàn)[14-15]。質(zhì)心系下(dσd,p/dω')(Ei,θ)的計(jì)算結(jié)果如圖2所示。

    圖2 質(zhì)心系下D(d,p)T反應(yīng)對應(yīng)實(shí)驗(yàn)室系θ=90o方向出射質(zhì)子的微分截面Fig.2 Differential cross sections for D(d,p)T in the centroid system corresponding to the laboratory system for the proton at θ=90o.

    實(shí)驗(yàn)室系D(d,n)3He和D(d,p)T反應(yīng)的積分截面σd,p和σd,n采用了ENDF-B-VI庫的推薦數(shù)據(jù),如圖3所示;式(6)、(7)所需的阻止本領(lǐng)數(shù)據(jù)由SRIM-2008程序計(jì)算給出[16]。

    圖3 實(shí)驗(yàn)室系D(d,n)3He和D(d,p)T反應(yīng)積分截面Fig.3 Integral reaction cross section of D(d,n)3He and D(d,p)T in the laboratory system.

    3 結(jié)果及討論

    根據(jù)上述計(jì)算方法及相關(guān)數(shù)據(jù),利用MATLAB軟件平臺,開發(fā)了一個計(jì)算程序。應(yīng)用所開發(fā)的程序,完成了入射氘束能量20-700 keV范圍,90o伴隨粒子法D-D中子產(chǎn)額測量所需的各向異性修正因子、中子和質(zhì)子產(chǎn)額比、總修正因子的理論計(jì)算,并將計(jì)算結(jié)果同早先Ruby的研究數(shù)據(jù)[9]進(jìn)行了比較,結(jié)果分別如圖4-6所示。

    圖4 各向異性因子Rthick計(jì)算結(jié)果與文獻(xiàn)[9]數(shù)據(jù)對比Fig.4 Result of anisotropy factor Rthickand the comparison with the Ref.[9] data.

    圖5 D(d,n)3He反應(yīng)中子和D(d,p)T反應(yīng)質(zhì)子積分產(chǎn)額比Fig.5 Yield ratio of the neutron for D(d,n)3He and the proton for D(d,p)T.

    圖6 總修正異性因子RY計(jì)算結(jié)果與文獻(xiàn)[9]數(shù)據(jù)對比Fig.6 Result of total correction factor RY and the comparison with the Ref.[9] data.

    由圖4可以看出,D束能量在100-500 keV能區(qū),本研究計(jì)算給出的各向異性修正因子Rthick結(jié)果與Ruby的數(shù)據(jù)[9]一致性較好,相對偏差小于0.8%;在D束流能量50 keV附近,兩者之間差異最大,最大偏差為13.6%(對應(yīng)D能量50 keV),這可能與Ruby使用了較早的截面數(shù)據(jù)有關(guān)。

    由圖5可以看出,在D束能量小于50 keV能區(qū),中子和質(zhì)子產(chǎn)額比小于1.0;在D束能量大于50 keV能區(qū),中子和質(zhì)子產(chǎn)額比大于1.0,并隨D束流能量而增大,這一特征是由圖3所示ENDF-B-VI庫給出的D(d,n)3He和D(d,p)T反應(yīng)積分截面新數(shù)據(jù)所導(dǎo)致,也是造成本研究計(jì)算給出的Rthick與Ruby給出的結(jié)果[9](Ruby使用了較早的反應(yīng)截面數(shù)據(jù))在D束能量50 keV附近存在較大差異的原因。

    由圖6可以看出,在D束能量50 keV能區(qū)附近,總修正因子RY的計(jì)算結(jié)果和Ruby給出的數(shù)據(jù)[9]之間也存在較大差異,最大偏差為16.2%;在D束能量200 keV 附近,差異較小,兩者之間的相對偏差約為4%,隨后隨著D束能量的增大,相對偏差逐漸增大,到500 keV附近,相對偏差約6%。上述偏差也與Ruby使用了較早的截面數(shù)據(jù)有關(guān)。

    4 結(jié)語

    對90o伴隨粒子法D-D中子產(chǎn)額測量中的修正因子計(jì)算方法進(jìn)行了研究,基于MATLAB軟件平臺,開發(fā)了用于修正因子計(jì)算的計(jì)算機(jī)程序,計(jì)算給出了厚靶條件下,入射氘能量在20-700 keV范圍,90o伴隨粒子法D-D中子產(chǎn)額測量各向異性修正因子Rthick、中子和質(zhì)子產(chǎn)額比(Yd,n/Yd,p)thick及總修正因子RY,并與早先的研究結(jié)果進(jìn)行了對比分析。就本工作的計(jì)算方法,只要厚靶被分割得足夠薄,薄靶近似所帶來的不確定度很小,可忽略。修正因子計(jì)算數(shù)據(jù)的不確定度主要來自于D離子在靶中的阻止本領(lǐng)數(shù)據(jù)和D-D反應(yīng)截面數(shù)據(jù)。其中,阻止本領(lǐng)計(jì)算的不確定度約為5%,由式(6)、(7)可以看出,方程中分子分母均包含阻止本領(lǐng)項(xiàng),阻止本領(lǐng)數(shù)據(jù)的不確定度相互抵消,故此項(xiàng)所帶來的不確定度也可以忽略不計(jì);ENDF-B-VI庫給出的D(d,n)3He和D(d,p)T反應(yīng)積分截面數(shù)據(jù)和文獻(xiàn)[14-15]給出的微分截面數(shù)據(jù)的不確定度約為2%,故修正因子計(jì)算數(shù)據(jù)的不確定度約為2%。

    1 魏國海, 韓松柏, 陳東風(fēng), 等. 中子照相技術(shù)在核燃料元件無損檢測中的應(yīng)用[J]. 核技術(shù), 2012, 35(11): 821-826. WEI Guohai, HAN Songbai, CHEN Dongfeng, et al. Application of neutron radiography for non-destructive testing nuclear fuel elements[J]. Nuclear Techniques, 2012, 35(11): 821-826.

    2 袁國軍, 肖才錦, 金象春, 等. 北京市場10種膳食魚微量元素的中子活化分析[J]. 核技術(shù), 2013, 36(9): 090207. DOI: 10.11889/j.0253-3219.2013.hjs.36.090207.

    YUAN Guojun, XIAO Caijin, JIN Xiangchun, et al. Study on trace elements of ten kinds of fish species at Beijing market by INAA[J]. Nuclear Techniques, 2013, 36(9): 090207. DOI: 10.11889/j.0253-3219.2013.hjs.36. 090207.

    3 Reijonen J, Gicquel F, Leung K N, et al. D-D neutron generator development at LBNL[J]. Applied Radiation and Isotopes, 2005, 63: 757-763. DOI: 10.1016/ j.apradiso.2005.05.024.

    4 Ludewigt B A, Wells R P, Reijonen J. High-yield D-T neutron generator[J]. Nuclear Instruments and Methods in Physics Research B, 2007, 261: 830-834. DOI: 10.1016/j.nimb.2007.04.246.

    5 Reijonen J. Neutron generators developed at LBNL for homeland security and imaging applications[J]. Nuclear Instruments and Methods in Physics Research B, 2007, 261: 272-276. DOI: 10.1016/j.nimb.2007.04.192.

    6 張國光, 沈冠仁, 豐樹強(qiáng), 等. 小型化中子發(fā)生器研制進(jìn)展[J]. 中國原子能科學(xué)研究院年報(bào), 2012, (1): 344-346.

    ZHANG Guoguang, SHEN Guanren, FENG Shuqiang, et al. The development of miniaturized neutron generator[J]. China Institute of Atomic Energy Annual Report, 2012, (1): 344-346.

    7 肖坤祥, 冉漢正, 曾清, 等. 高產(chǎn)額中子發(fā)生器研制[J].原子能科學(xué)技術(shù), 2012, 46(B09): 713-717. DOI: 1000-6931(2012)S0-0713-05.

    XIAO Kunxiang, RAN Hanzheng, ZENG Qing, et al. Development of high-yielding neutron generator[J]. Atomic Energy Science and Technology, 2012, 46(B09): 713-717. DOI: 1000-6931(2012)S0-0713-05.

    8 Yan Y, Lu X L, Yao Z E, et al. Study on the method with associated particle for measuring the neutron yield of D-D neutron generator[J]. Nuclear Instruments and Methods in Physics Research A, 2015, 773: 52-55. DOI: 10.1016/ j.nima.2014.11.023.

    9 Ruby L, Crawford R B. Anisotropy factors for determination of total neutron yield from the D(d,n)3He and T(d,n)4He reactions[J]. Nuclear Instruments and Methods, 1963, 24: 413-417.

    10 Jordanova J. Anisotropy factors for determination of neutron yield from the2H(d,n)3He and T(d,n)4He reaction[J]. Journal of Radio Analytical and Nuclear Chemistry, 1995, 200: 127-135.

    11 中國原子能科學(xué)研究院. 加速器單能中子源常用數(shù)據(jù)手冊[M]. 北京: 科學(xué)出版社, 1976: 215-217.

    China Institute of Atomic Energy. Accelerator monoenergetic neutron source common data manual[M]. Beijing: Science Press, 1976: 215-217.

    12 盧希庭. 原子核物理[M]. 北京: 原子能出版社, 2000: 235-236.

    LU Xiting. Nuclear physics[M]. Beijing: Atomic Energy Press, 2000: 235-236.

    13 Csikai G J. CRC handbook of fast neutron generators[M]. Vol.1. CRC Press Inc Boca Raton Fl, 1987.

    14 Brown R E, Jarmie N. Differential cross sections at low energies for2H(d,p)3H and2H(d,n)3He[J]. Physical Review C, 1990, 41(4): 1391-1400.

    15 Theus R B, McCarry W I, Beach L A. Angular distributions and cross-section ratios for the reactions2H(d,n)3He and2H(d,p)3H below 500 keV[J]. Nuclear Physics, 1966, 80: 273-288.

    16 James F Z. Beam interactions with materials and atoms[J]. Nuclear Instruments and Methods in Physics Research B, 2004, 219: 1027-1036.

    Calculation on the associated particle method for 90o to measure the neutron yield correction factors of D-D neutron generator

    LI Jianyi1WANG Junrun1,2ZHANG Yu1,2HUANG Zhiwu1LU Xiaolong1,2XU Dapeng1,2ZHANG Jie1WEI Zheng1,2MA Zhanwen1YAO Zeen1,2
    1(School of Nuclear Science and Technology, Lanzhou University, Lanzhou 730000, China)
    2(Engineering Research Center for Neutron Application, Ministry of Education, Lanzhou University, Lanzhou 730000, China)

    Background:In recent years, the neutron generator based on D-D reaction is widely investigated due to some potential applications and the advantage of the more security without the use of radioactive tritium target and long target lifetime with D ion self-injection effect in target. The accurate measurement of fast neutron yield of the D-D neutron generator is a key problem to be solved. For D-D neutron generator, the associated particle method is possibly one of the most accurate ways for the neutron source strength calibration. Purpose: In order to measure the neutron yield of D-D neutron generator using the associated particle method, the anisotropy correction factor and the yield ratio of neutron and proton need to be calculated. Methods: A method for calculating the anisotropy correction factor and the yield ratio and a computer program based on MATLAB software platform were developed. Results: In this work, the anisotropy correction factors, the yield ratio of neutron and proton and the total correction factors were presented for the thick-target and for 90° associated particle measurement in the deuterons energy region of20-700 keV. Conclusion: The uncertainty of the total correction factors (RY) comes mainly from the cross section data. The uncertainty of RYis estimated to be about 2%.

    LI Jianyi, male, born in 1990, graduated from Lanzhou University in 2014, master student, focusing on neutron physics and application technology

    YAO Zeen, E-mail: zeyao@lzu.edu.cn

    Associated particle method, D-D neutron yield, Neutron generator, Correction factors

    O57

    10.11889/j.0253-3219.2017.hjs.40.010201

    No.11375077)、國家重大科學(xué)儀器設(shè)備開發(fā)專項(xiàng)(No.2013YQ40861)、蘭州大學(xué)中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金(No.lzujbky-2016-bt08)資助

    李建一,男,1990年出生,2014年畢業(yè)于蘭州大學(xué),現(xiàn)為碩士研究生,主要研究領(lǐng)域?yàn)橹凶游锢砑皯?yīng)用技術(shù)

    姚澤恩,E-mail: zeyao@lzu.edu.cn

    2016-09-29,

    2016-11-16

    Supported by National Natural Science Foundation of China (No.11375077), National Key Scientific Instrument and Equipment Development Project (No.2013YQ40861), Fundamental Research Funds for the Central Universities (No.lzujbky-2016-bt08)

    Received date: 2016-09-29, accepted date: 2016-11-16

    猜你喜歡
    測量
    測量重量,測量長度……
    把握四個“三” 測量變簡單
    滑動摩擦力的測量和計(jì)算
    滑動摩擦力的測量與計(jì)算
    測量的樂趣
    二十四節(jié)氣簡易測量
    日出日落的觀察與測量
    滑動摩擦力的測量與計(jì)算
    測量
    測量水的多少……
    欧美97在线视频| 亚洲av福利一区| 在线观看免费午夜福利视频| 国产精品一区二区在线不卡| 侵犯人妻中文字幕一二三四区| 久久久久久人人人人人| 亚洲国产精品999| 日韩一本色道免费dvd| 国产精品一国产av| 看十八女毛片水多多多| 97在线人人人人妻| 无遮挡黄片免费观看| 日韩电影二区| kizo精华| 80岁老熟妇乱子伦牲交| 国产在线免费精品| 亚洲一码二码三码区别大吗| 色精品久久人妻99蜜桃| 国产一区二区三区av在线| 中文字幕人妻熟女乱码| 午夜激情久久久久久久| 成年女人毛片免费观看观看9 | 熟妇人妻不卡中文字幕| 国产一区有黄有色的免费视频| 免费观看性生交大片5| 国产精品一区二区在线不卡| 亚洲国产精品一区三区| 国产日韩一区二区三区精品不卡| 91精品伊人久久大香线蕉| 国产精品一区二区精品视频观看| 老司机靠b影院| 悠悠久久av| 纯流量卡能插随身wifi吗| 欧美精品高潮呻吟av久久| 国产日韩欧美在线精品| 人人妻人人爽人人添夜夜欢视频| 国产免费现黄频在线看| av网站免费在线观看视频| 哪个播放器可以免费观看大片| 成年女人毛片免费观看观看9 | 日韩一区二区视频免费看| 日韩欧美一区视频在线观看| 亚洲成人国产一区在线观看 | 久久精品国产亚洲av涩爱| 免费观看人在逋| 国产精品一国产av| 爱豆传媒免费全集在线观看| 亚洲国产精品一区三区| av卡一久久| 亚洲欧美中文字幕日韩二区| 久久久久久人人人人人| 成年动漫av网址| 国产毛片在线视频| 亚洲中文av在线| 欧美激情极品国产一区二区三区| 欧美黄色片欧美黄色片| 亚洲在久久综合| 美女高潮到喷水免费观看| 欧美最新免费一区二区三区| 天堂中文最新版在线下载| 国产一区二区三区av在线| 午夜老司机福利片| 大码成人一级视频| 成人18禁高潮啪啪吃奶动态图| 一区二区三区四区激情视频| 国产在视频线精品| 久久久精品区二区三区| 新久久久久国产一级毛片| 成人午夜精彩视频在线观看| 成人午夜精彩视频在线观看| 美女高潮到喷水免费观看| 成人毛片60女人毛片免费| 成人手机av| 国产成人免费观看mmmm| 亚洲国产精品一区三区| 又黄又粗又硬又大视频| 亚洲精品成人av观看孕妇| 黑人巨大精品欧美一区二区蜜桃| 巨乳人妻的诱惑在线观看| 国产一卡二卡三卡精品 | 久久久精品国产亚洲av高清涩受| 国产精品亚洲av一区麻豆 | 国产精品欧美亚洲77777| av在线观看视频网站免费| 99精品久久久久人妻精品| 免费在线观看黄色视频的| 日韩大片免费观看网站| 欧美久久黑人一区二区| 男男h啪啪无遮挡| 9热在线视频观看99| 中文字幕亚洲精品专区| 日韩 欧美 亚洲 中文字幕| 老汉色av国产亚洲站长工具| 美女脱内裤让男人舔精品视频| 亚洲av国产av综合av卡| 国产国语露脸激情在线看| 777米奇影视久久| 91老司机精品| 久久精品久久精品一区二区三区| 亚洲自偷自拍图片 自拍| 黄网站色视频无遮挡免费观看| av线在线观看网站| 国产又色又爽无遮挡免| 在线天堂最新版资源| 精品福利永久在线观看| 国产欧美日韩综合在线一区二区| 熟女少妇亚洲综合色aaa.| 校园人妻丝袜中文字幕| 女性被躁到高潮视频| 久久久久久人妻| 秋霞伦理黄片| 午夜福利在线免费观看网站| 婷婷色麻豆天堂久久| 亚洲欧洲国产日韩| 丰满乱子伦码专区| 亚洲综合精品二区| 丝袜在线中文字幕| 成人午夜精彩视频在线观看| 最近中文字幕2019免费版| 91国产中文字幕| 国产男人的电影天堂91| 啦啦啦在线免费观看视频4| 下体分泌物呈黄色| 亚洲国产成人一精品久久久| 亚洲精品国产av蜜桃| 欧美变态另类bdsm刘玥| 国产av码专区亚洲av| 超色免费av| 亚洲国产欧美在线一区| 国产视频首页在线观看| 99久国产av精品国产电影| 久久久久久久国产电影| 欧美老熟妇乱子伦牲交| 亚洲第一区二区三区不卡| 哪个播放器可以免费观看大片| 国产精品欧美亚洲77777| 又黄又粗又硬又大视频| 男人舔女人的私密视频| 人妻一区二区av| 亚洲国产欧美一区二区综合| 国产97色在线日韩免费| 日本爱情动作片www.在线观看| 水蜜桃什么品种好| 如何舔出高潮| 一本久久精品| 丝瓜视频免费看黄片| 日本黄色日本黄色录像| 男女边摸边吃奶| 少妇的丰满在线观看| 又粗又硬又长又爽又黄的视频| 久久性视频一级片| www.精华液| 久久精品亚洲熟妇少妇任你| 精品亚洲乱码少妇综合久久| 日韩,欧美,国产一区二区三区| 免费人妻精品一区二区三区视频| 啦啦啦在线观看免费高清www| 大片免费播放器 马上看| 麻豆精品久久久久久蜜桃| 午夜福利视频在线观看免费| 一本一本久久a久久精品综合妖精| 免费观看av网站的网址| 日韩制服骚丝袜av| 免费观看性生交大片5| 在线观看www视频免费| 考比视频在线观看| 精品福利永久在线观看| 国产男人的电影天堂91| 国产精品欧美亚洲77777| 精品一区二区三区av网在线观看 | 日本欧美国产在线视频| 国产日韩欧美亚洲二区| 性高湖久久久久久久久免费观看| 亚洲熟女精品中文字幕| 青青草视频在线视频观看| 国产精品国产三级国产专区5o| 久久精品aⅴ一区二区三区四区| 国产伦人伦偷精品视频| 日本vs欧美在线观看视频| 国产成人免费无遮挡视频| 999精品在线视频| 尾随美女入室| 亚洲人成电影观看| 久久午夜综合久久蜜桃| 丰满饥渴人妻一区二区三| 国产午夜精品一二区理论片| 人成视频在线观看免费观看| av又黄又爽大尺度在线免费看| 最近手机中文字幕大全| 18禁观看日本| 国产欧美亚洲国产| 看非洲黑人一级黄片| 最近中文字幕高清免费大全6| 国产无遮挡羞羞视频在线观看| 久久亚洲国产成人精品v| 免费高清在线观看视频在线观看| 两个人看的免费小视频| 青春草国产在线视频| 久久亚洲国产成人精品v| 免费人妻精品一区二区三区视频| av天堂久久9| 少妇 在线观看| h视频一区二区三区| 最近最新中文字幕大全免费视频 | 国产精品欧美亚洲77777| 高清欧美精品videossex| 新久久久久国产一级毛片| 天天添夜夜摸| 极品人妻少妇av视频| 99久久99久久久精品蜜桃| 成年av动漫网址| 伦理电影大哥的女人| 日本午夜av视频| 免费观看性生交大片5| 亚洲国产欧美日韩在线播放| 老汉色av国产亚洲站长工具| 一区在线观看完整版| 亚洲精品久久午夜乱码| 热re99久久精品国产66热6| 熟女av电影| 两性夫妻黄色片| 综合色丁香网| 中文字幕制服av| videosex国产| 日韩不卡一区二区三区视频在线| 99久久综合免费| 伊人久久国产一区二区| 精品国产露脸久久av麻豆| 国产一区亚洲一区在线观看| 热99国产精品久久久久久7| 秋霞伦理黄片| 日韩大码丰满熟妇| 国产一区二区 视频在线| 免费看av在线观看网站| 国产一区二区三区综合在线观看| 国产 一区精品| 毛片一级片免费看久久久久| 日韩制服丝袜自拍偷拍| 午夜福利在线免费观看网站| 99热网站在线观看| 天天躁日日躁夜夜躁夜夜| 免费黄频网站在线观看国产| 永久免费av网站大全| 搡老乐熟女国产| 悠悠久久av| 精品酒店卫生间| 成人国产麻豆网| 一级毛片我不卡| 日韩 欧美 亚洲 中文字幕| 麻豆av在线久日| 午夜免费鲁丝| 老司机靠b影院| 亚洲av日韩精品久久久久久密 | 日韩一本色道免费dvd| 搡老乐熟女国产| 女人被躁到高潮嗷嗷叫费观| videosex国产| 亚洲精品乱久久久久久| 在线 av 中文字幕| 欧美变态另类bdsm刘玥| 在线观看免费日韩欧美大片| 看免费av毛片| 韩国精品一区二区三区| 成人亚洲欧美一区二区av| 国产成人免费观看mmmm| 国产老妇伦熟女老妇高清| 两个人免费观看高清视频| 日韩大码丰满熟妇| 一区二区三区四区激情视频| 亚洲国产精品国产精品| 肉色欧美久久久久久久蜜桃| 熟女少妇亚洲综合色aaa.| 精品少妇一区二区三区视频日本电影 | 亚洲欧美日韩另类电影网站| 久久久欧美国产精品| 两性夫妻黄色片| 精品少妇内射三级| 久久婷婷青草| 如日韩欧美国产精品一区二区三区| 天美传媒精品一区二区| 韩国精品一区二区三区| 男女下面插进去视频免费观看| 久久久久国产一级毛片高清牌| 久久精品国产综合久久久| 日本91视频免费播放| 高清黄色对白视频在线免费看| 亚洲综合色网址| 一二三四在线观看免费中文在| 18禁国产床啪视频网站| 亚洲天堂av无毛| 精品卡一卡二卡四卡免费| 国产一区有黄有色的免费视频| 电影成人av| 久久久久久久精品精品| 91成人精品电影| 一级,二级,三级黄色视频| 久久青草综合色| 国产精品麻豆人妻色哟哟久久| 母亲3免费完整高清在线观看| 搡老乐熟女国产| 999久久久国产精品视频| 久久人妻熟女aⅴ| 国产亚洲最大av| 亚洲成人一二三区av| 90打野战视频偷拍视频| 成人国产av品久久久| 国产欧美日韩一区二区三区在线| 欧美日韩综合久久久久久| 国产精品三级大全| 精品国产超薄肉色丝袜足j| 午夜福利在线免费观看网站| 大片电影免费在线观看免费| 亚洲精品在线美女| 国产熟女欧美一区二区| 亚洲国产中文字幕在线视频| 老司机靠b影院| 国产熟女欧美一区二区| 丰满乱子伦码专区| 在线亚洲精品国产二区图片欧美| 亚洲第一区二区三区不卡| 成人影院久久| 自线自在国产av| 男女之事视频高清在线观看 | 婷婷色综合大香蕉| 国产免费一区二区三区四区乱码| 久久久久国产一级毛片高清牌| 日本wwww免费看| 最近中文字幕高清免费大全6| 高清在线视频一区二区三区| 国产探花极品一区二区| 成人国产av品久久久| 亚洲自偷自拍图片 自拍| 欧美在线黄色| 69精品国产乱码久久久| 色播在线永久视频| 午夜福利乱码中文字幕| 制服丝袜香蕉在线| 日韩人妻精品一区2区三区| 黑人欧美特级aaaaaa片| 亚洲av成人不卡在线观看播放网 | 热99国产精品久久久久久7| 亚洲男人天堂网一区| e午夜精品久久久久久久| 精品少妇一区二区三区视频日本电影 | 精品少妇黑人巨大在线播放| 久久久久久久久免费视频了| 尾随美女入室| 日韩成人av中文字幕在线观看| 国产成人一区二区在线| xxxhd国产人妻xxx| 欧美激情极品国产一区二区三区| 午夜福利视频精品| 免费高清在线观看日韩| 一本大道久久a久久精品| 欧美日韩亚洲综合一区二区三区_| 亚洲成人手机| 香蕉丝袜av| 亚洲国产最新在线播放| 黄频高清免费视频| 亚洲一区中文字幕在线| 电影成人av| 国产黄色视频一区二区在线观看| 这个男人来自地球电影免费观看 | 波多野结衣一区麻豆| 国产一级毛片在线| 成人国产麻豆网| www.自偷自拍.com| 亚洲美女视频黄频| 久久久国产欧美日韩av| 一级片'在线观看视频| 巨乳人妻的诱惑在线观看| 成人国语在线视频| 国产又爽黄色视频| 免费久久久久久久精品成人欧美视频| 一区二区日韩欧美中文字幕| 青春草亚洲视频在线观看| 天美传媒精品一区二区| 成人免费观看视频高清| 欧美日韩国产mv在线观看视频| 色网站视频免费| 一本久久精品| 天天操日日干夜夜撸| 女性被躁到高潮视频| 天天操日日干夜夜撸| 水蜜桃什么品种好| 欧美日韩成人在线一区二区| 如何舔出高潮| av福利片在线| 国产日韩欧美亚洲二区| 又粗又硬又长又爽又黄的视频| 免费少妇av软件| 亚洲国产欧美在线一区| 亚洲人成网站在线观看播放| 精品第一国产精品| 啦啦啦 在线观看视频| 免费观看性生交大片5| 亚洲中文av在线| 99久久综合免费| 亚洲成色77777| 久久97久久精品| 国产精品.久久久| av国产久精品久网站免费入址| 国产男女超爽视频在线观看| 一本大道久久a久久精品| 久久精品国产亚洲av高清一级| 欧美 亚洲 国产 日韩一| 你懂的网址亚洲精品在线观看| 午夜91福利影院| 黄色视频不卡| 亚洲,一卡二卡三卡| 国产在线免费精品| 最近手机中文字幕大全| 天天躁夜夜躁狠狠久久av| 亚洲精品第二区| 亚洲精品日韩在线中文字幕| 国产精品久久久久久精品古装| 婷婷色综合大香蕉| 久久精品国产a三级三级三级| 成人国产av品久久久| 成人黄色视频免费在线看| 国产欧美日韩一区二区三区在线| 一级毛片 在线播放| 日韩 欧美 亚洲 中文字幕| 黑人猛操日本美女一级片| 亚洲国产精品999| 色视频在线一区二区三区| 国产精品.久久久| av电影中文网址| 午夜福利网站1000一区二区三区| 建设人人有责人人尽责人人享有的| 如日韩欧美国产精品一区二区三区| 亚洲国产成人一精品久久久| 成人国产av品久久久| 搡老乐熟女国产| 国产探花极品一区二区| 亚洲欧美一区二区三区久久| 满18在线观看网站| 最近的中文字幕免费完整| 欧美亚洲 丝袜 人妻 在线| 婷婷成人精品国产| 99香蕉大伊视频| 乱人伦中国视频| 各种免费的搞黄视频| 欧美黑人精品巨大| 欧美日韩视频高清一区二区三区二| 伦理电影大哥的女人| 婷婷色综合www| 纵有疾风起免费观看全集完整版| 欧美成人午夜精品| 少妇被粗大猛烈的视频| 一个人免费看片子| 18禁动态无遮挡网站| 日韩伦理黄色片| 狠狠精品人妻久久久久久综合| 精品一区二区三区四区五区乱码 | 国产精品二区激情视频| 国产成人免费无遮挡视频| 在线天堂最新版资源| 在线 av 中文字幕| 两性夫妻黄色片| 亚洲av在线观看美女高潮| 精品午夜福利在线看| 久久天堂一区二区三区四区| 国产av国产精品国产| 人体艺术视频欧美日本| 丝袜脚勾引网站| 如日韩欧美国产精品一区二区三区| 99久久人妻综合| 大片电影免费在线观看免费| 99re6热这里在线精品视频| 男男h啪啪无遮挡| 满18在线观看网站| 欧美成人午夜精品| 亚洲国产av影院在线观看| 巨乳人妻的诱惑在线观看| 久久影院123| 大香蕉久久网| 亚洲一区中文字幕在线| 国产 一区精品| 精品亚洲成国产av| 久久午夜综合久久蜜桃| 看十八女毛片水多多多| 9热在线视频观看99| 亚洲美女黄色视频免费看| 久久精品国产a三级三级三级| 久久久久视频综合| 国产又爽黄色视频| 亚洲综合色网址| 欧美97在线视频| 成人国语在线视频| 精品人妻一区二区三区麻豆| 国产国语露脸激情在线看| 飞空精品影院首页| 国产极品粉嫩免费观看在线| 热re99久久精品国产66热6| 亚洲精品久久久久久婷婷小说| 国产精品免费视频内射| 亚洲精品自拍成人| 国产精品久久久久久精品电影小说| 国产av码专区亚洲av| 欧美人与性动交α欧美软件| 国产av精品麻豆| av卡一久久| 成人亚洲精品一区在线观看| 日日摸夜夜添夜夜爱| 日韩,欧美,国产一区二区三区| 黄片播放在线免费| 国产一卡二卡三卡精品 | 老汉色∧v一级毛片| 久久性视频一级片| 热re99久久精品国产66热6| 国产精品成人在线| 亚洲三区欧美一区| 免费高清在线观看视频在线观看| 青春草国产在线视频| 日本av免费视频播放| 国产一级毛片在线| 欧美人与性动交α欧美精品济南到| 老汉色∧v一级毛片| 亚洲一卡2卡3卡4卡5卡精品中文| 国产无遮挡羞羞视频在线观看| 欧美精品一区二区免费开放| 久久久精品区二区三区| 欧美精品一区二区免费开放| 蜜桃国产av成人99| 欧美人与性动交α欧美精品济南到| 悠悠久久av| 亚洲欧洲国产日韩| 亚洲美女搞黄在线观看| 日韩制服骚丝袜av| 国产亚洲av片在线观看秒播厂| 色婷婷av一区二区三区视频| 伊人久久国产一区二区| 18禁裸乳无遮挡动漫免费视频| 男女午夜视频在线观看| 亚洲av成人精品一二三区| 免费黄网站久久成人精品| 亚洲国产av新网站| 精品人妻一区二区三区麻豆| 18禁裸乳无遮挡动漫免费视频| 精品一区二区三卡| 激情五月婷婷亚洲| 免费在线观看完整版高清| 搡老乐熟女国产| 深夜精品福利| 校园人妻丝袜中文字幕| 啦啦啦 在线观看视频| 在线观看免费高清a一片| 午夜福利免费观看在线| 无限看片的www在线观看| 亚洲av日韩在线播放| 曰老女人黄片| 亚洲av国产av综合av卡| 妹子高潮喷水视频| 亚洲欧洲国产日韩| 丁香六月天网| 观看美女的网站| 亚洲七黄色美女视频| 女的被弄到高潮叫床怎么办| 国产亚洲午夜精品一区二区久久| 五月开心婷婷网| 在线免费观看不下载黄p国产| 日本wwww免费看| 欧美老熟妇乱子伦牲交| 免费在线观看视频国产中文字幕亚洲 | 久久久欧美国产精品| 丰满少妇做爰视频| 久久性视频一级片| 十八禁高潮呻吟视频| 久久午夜综合久久蜜桃| 免费女性裸体啪啪无遮挡网站| 国产成人欧美| 高清av免费在线| 亚洲成色77777| 国产精品一二三区在线看| 中文字幕色久视频| 久久鲁丝午夜福利片| 国产精品久久久久成人av| 男的添女的下面高潮视频| 青春草国产在线视频| 亚洲欧美成人综合另类久久久| 成年人午夜在线观看视频| 人妻一区二区av| 久久97久久精品| 大话2 男鬼变身卡| 久久久久国产精品人妻一区二区| 国产极品粉嫩免费观看在线| 国产野战对白在线观看| 国产成人啪精品午夜网站| 欧美黄色片欧美黄色片| 日日爽夜夜爽网站| 亚洲欧美一区二区三区黑人| 国产精品一区二区精品视频观看| 亚洲色图 男人天堂 中文字幕| 免费黄网站久久成人精品| 777米奇影视久久| 1024香蕉在线观看| 精品卡一卡二卡四卡免费| 美女大奶头黄色视频| 另类精品久久| 日日撸夜夜添| 亚洲欧美成人精品一区二区| 亚洲,一卡二卡三卡| 欧美激情高清一区二区三区 | 在线天堂中文资源库| 亚洲,一卡二卡三卡| 亚洲精品自拍成人| 国产精品久久久av美女十八| 久久97久久精品| 日韩大码丰满熟妇| 精品国产一区二区三区四区第35| 色播在线永久视频| 波多野结衣一区麻豆| 午夜免费男女啪啪视频观看| 又黄又粗又硬又大视频|