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

    Geant4模擬1.6 GeV質(zhì)子的輸運(yùn)過程

    2014-03-22 02:42:46姚志明宋顧周黑東煒馬繼明段寶軍韓長材
    核技術(shù) 2014年3期
    關(guān)鍵詞:面密度庫侖原子核

    姚志明 宋顧周 黑東煒 馬繼明 周 鳴 段寶軍 宋 巖 韓長材

    (西北核技術(shù)研究所 強(qiáng)脈沖輻射環(huán)境模擬與效應(yīng)國家重點(diǎn)實驗室 西安 710024)

    輻射成像技術(shù)是通過測量入射到被測物體上的射線束的參數(shù)變化來確定物體內(nèi)部幾何結(jié)構(gòu)和材料組分等信息的[1]。用于射線檢測的射線源包括X射線、γ射線、中子、電子、α粒子、質(zhì)子等。其中質(zhì)子是帶電強(qiáng)子,高能量的質(zhì)子具有穿透能力強(qiáng)、散射角分布與材料組分相關(guān)、源的單色性好和探測效率高等特點(diǎn)。美國科學(xué)家利用質(zhì)子作為探針做了大量的理論和實驗工作[1?4]。

    質(zhì)子的電離能損率與電子密度近似成正比,可以用于 CT技術(shù)[2];能量損失在射程末端達(dá)到極大值,可以用于治療癌癥[3];在不同材料中散射角不同,可以獲得高對比度的邊界圖像[1];與原子核碰撞后發(fā)生衰減,多次庫侖散射角與材料有關(guān),可以用于高能質(zhì)子照相技術(shù)[4]??傊|(zhì)子在醫(yī)學(xué)、材料科學(xué)、輻射成像等領(lǐng)域有著廣闊的應(yīng)用前景。

    然而,受制于質(zhì)子加速器的昂貴造價,除醫(yī)學(xué)上的應(yīng)用外,國內(nèi)的質(zhì)子照相技術(shù)尚處于理論分析和蒙卡模擬的起步階段。將于2018年建成的中國散裂中子源中的質(zhì)子加速器可以將質(zhì)子加速到1.6GeV[5]。本文將通過Geant4蒙卡模擬軟件對該能量下質(zhì)子在材料中的輸運(yùn)過程進(jìn)行模擬計算,分析個數(shù)衰減、能量損失、散射角分布等參數(shù)隨材料的面密度和組分的變化規(guī)律。

    1 質(zhì)子與物質(zhì)相互作用的理論公式

    質(zhì)子與物質(zhì)相互作用的物理過程包括電離能損、多次庫侖散射、與原子核的非彈性碰撞以及與原子核的彈性碰撞[4]。其中與原子核的彈性碰撞又可按碰撞后原子核處于基態(tài)或激發(fā)態(tài)劃分為彈性散射和準(zhǔn)彈性散射[6]。對于同一種物理過程不同的文獻(xiàn)中給出的經(jīng)驗公式有所差別,以下各選取了一種進(jìn)行論述。

    1.1 電離能損

    質(zhì)子與電子發(fā)生庫侖相互作用,能量發(fā)生衰減,能損率由Bethe-Bloch公式[6]計算:

    式 中, 能損率 的單位 是 MeV·g?1·cm2,K=0.307075 MeV·g?1·cm2;z為事件粒子電荷量(質(zhì)子z=1);Z為靶原子電荷量;A為靶原子原子質(zhì)量;b為質(zhì)子速度除以光速,即b=v/c;me為電子靜止質(zhì)量;γ2=1/(1?b2);Tmax是單次碰撞中質(zhì)子能夠傳遞給電子的最大動能;I是靶原子的平均激發(fā)能;C/Z是殼修正項;d/2是密度修正量。

    1.2 多次庫侖散射

    質(zhì)子受到原子核庫侖力的作用,發(fā)生多次庫侖散射(Multiple Coulomb Scattering, MCS)。每次散射后,質(zhì)子能量不變,運(yùn)動方向發(fā)生小的改變??臻g角分布可以用高斯分布來描述[4]:

    式中,θ是質(zhì)子偏離初始方向的角度;θ0是多次庫侖散射角的均方根值,可用下式估計:

    式中,p是質(zhì)子動量;β為質(zhì)子速度與光速的比值;li為面密度;Ri是材料的輻射長度,經(jīng)驗公式為:

    式中,Ri單位是g·cm?2;Z為靶原子電荷量;A為靶原子原子質(zhì)量。

    1.3 核反應(yīng)

    高能質(zhì)子與原子核中的質(zhì)子和中子發(fā)生非彈性碰撞。質(zhì)子個數(shù)以指數(shù)形式衰減[4]:

    式中,N為透射質(zhì)子個數(shù);N0為入射質(zhì)子個數(shù);li為面密度;λi是平均自由程,單位g·cm?2。λi可用下式估計:

    式中,A為相對原子質(zhì)量;NA為阿伏伽德羅常數(shù);A/NA為每個原子的質(zhì)量,g;σi是核反應(yīng)截面,可用下式估計:

    式中,ri=r0×A1/3,r0≈1.2 fm。

    1.4 彈性散射

    質(zhì)子與原子核發(fā)生碰撞,能量不發(fā)生損失或者損失非常小的部分,運(yùn)動方向發(fā)生改變。與電離和MCS不同,核反應(yīng)和彈性散射并不是每個質(zhì)子都會發(fā)生,而是有一定的概率會發(fā)生。

    如果碰撞時質(zhì)子能量不發(fā)生損失,原子核處于基態(tài),不產(chǎn)生粒子。微分截面可以由下式估計[6]:

    式中,p是質(zhì)子動量;h是普朗克常量;J1為一階貝塞爾函數(shù);R為原子核黑體有效半徑k=p/h。

    如果質(zhì)子能量損失非常小的部分,原子核將處于激發(fā)態(tài),或者發(fā)射出粒子。微分截面可以由下式估計[6]:

    2 Geant4的參數(shù)設(shè)置

    Geant4是CERN開發(fā)的一款用于模擬粒子輸運(yùn)過程的軟件工具包。與MCNP、EGS等蒙卡模擬軟件相比,Geant4具有源代碼完全開放的優(yōu)勢,用戶可以根據(jù)實際需要改進(jìn)和擴(kuò)展程序。Geant4用戶可以設(shè)置的參數(shù)主要包括物體的幾何結(jié)構(gòu)、粒子與物質(zhì)反應(yīng)的物理過程模型、粒子源的信息、粒子記錄的信息等[7]。

    Geant4的幾何結(jié)構(gòu)設(shè)置如圖1。整個“world”幾何體由真空填充,在“world”中心放置一定厚度的足夠大平板材料,距離平板5 mm位置放置足夠大的平板探測平面,用于記錄穿過物體后的粒子信息。粒子與物質(zhì)反應(yīng)的物理過程模型選取為G4hMultipleScattering、G4hIonisation、G4LElastic和G4ProtonInelasticProcess[8],分別對應(yīng)質(zhì)子與物質(zhì)的四種相互作用。能量為1.6 GeV的質(zhì)子在距離物體一定距離的位置沿z軸方向發(fā)射。探測平面記錄下透射粒子的能量、動量方向和位置信息。

    圖1 幾何結(jié)構(gòu)模型Fig.1 Geometry model.

    3 蒙卡模擬結(jié)果與討論

    3.1 電離能損

    3.1.1 單一材料中的能量損失

    Geant4中物理模型設(shè)置為G4hIonisation,計算單個質(zhì)子在足夠厚的鎢材料中的輸運(yùn)過程。圖2統(tǒng)計出每經(jīng)過1 cm厚鎢材料的能量沉積大小和剩余能量。

    可以看出,質(zhì)子在單一材料中的能量衰減規(guī)律是:在射程的前86%,單位長度的能量損失變化不大,而在射程的末端,能量損失率突然增大,稱為Bragg峰。利用這一特性,質(zhì)子被用于腫瘤治療[3]。

    由式(1),穿過單位面密度的能量損失率與材料的Z/A成正比。穿過一定面密度材料后的能量損失大小可以由下式給出:

    式中,ρ是體密度;S(E)是能損率中除去Z/A的部分;ηe是電子密度。

    進(jìn)而可以得到以下表達(dá)式:

    綜上所述,生長抑素輔助治療急性胰腺炎患者的療效顯著,可有效改善Try-2、CER及AMY水平,同時降低TNF-α、IL-6水平。具有較高的臨床推廣應(yīng)用價值。

    即測得入射和出射質(zhì)子能量,就可以重建電子密度,該方法被應(yīng)用于質(zhì)子CT技術(shù)[2]。

    圖2 質(zhì)子在鎢中的能量損失Fig.2 Proton energy loss in tungsten.

    3.1.2 不同材料中的能量損失

    對于不同面密度的鈹、銅和鎢材料,模擬計算了1.6 GeV質(zhì)子透射物體后的剩余能量,得到質(zhì)子的能量損失。三種材料的能量損失隨面密度的變化曲線如圖3所示。

    圖3 能量損失隨面密度變化曲線Fig.3 Energy loss vs. thickness.

    由圖 3,能量損失與面密度近似成正比關(guān)系,能量損失的大小可以反映已知材料的面密度信息;相同面密度下,不同材料的能量損失大小也不同,因而測得能量損失也可以用于區(qū)別不同的材料。

    3.2 多次庫侖散射

    3.2.1 單一材料中的散射分布

    Geant4物理模型設(shè)置為G4hMultipleScattering,選取面密度為5 g·cm?2、10 g·cm?2和15 g·cm?2的鎢材料,分別入射105個質(zhì)子,記錄下每個質(zhì)子偏離入射方向的角度,統(tǒng)計出角度分布情況,結(jié)果如圖4、5所示。

    圖4 多次庫侖散射角分布圖(dN/dW)Fig.4 MCS angular distribution (dN/dW).

    圖5 多次庫侖散射角分布圖(dN/dθ)Fig.5 MCS angular distribution (dN/dθ).

    這是高斯分布在空間立體角內(nèi)的積分,是高斯分布和正弦函數(shù)的乘積,縱坐標(biāo)中的個數(shù)除以角度的正弦值,就得到了高斯分布。

    式(3)表明散射角的均方根值與材料的輻射長度有關(guān)。對于不同面密度的物體,統(tǒng)計出了鈹、銅和鎢的多次庫侖散射角的均方根值,如圖6所示。

    圖6 θ0隨面密度變化曲線Fig.6 θ0 vs. thickness.

    可以看出,對于單一材料,θ0隨面密度而增大,測得θ0的大小就反映了面密度的大小。當(dāng)已知材料面密度時,不同材料的θ0的值不同,測得θ0的大小就反映了材料的組分信息。然而引起散射的因素還包括與原子核的彈性碰撞以及與原子核非彈性碰撞產(chǎn)生的次級質(zhì)子,θ0能否精確測量需要更多的探索。

    3.3 非彈性碰撞

    Geant4物理模型設(shè)置為G4ProtonInelasticProc,對于不同面密度的鈹、銅和鎢,統(tǒng)計出了能量未發(fā)生改變,即未發(fā)生非彈性碰撞的透射質(zhì)子個數(shù),圖7畫出了透射率隨面密度的變化曲線。

    圖7 透射率隨面密度變化曲線Fig.7 Transmission vs. thickness.

    由圖 7,對于單一材料,透射率隨面密度的增加呈指數(shù)衰減,與式(5)符合得較好。測得透射率大小就反映了面密度的大小。當(dāng)已知材料面密度時,不同材料的透射率的值不同,測得透射率的大小就反映了材料的組分信息。

    非彈性碰撞還會產(chǎn)生次級質(zhì)子,為了分析次級質(zhì)子的特性,選取面密度為5 g·cm?2、10 g·cm?2和15 g·cm?2的鎢材料進(jìn)行計算,物理模型設(shè)置為G4ProtonInelasticProcess,每次入射質(zhì)子個數(shù)105。模擬計算中發(fā)現(xiàn)有能量很低的質(zhì)子產(chǎn)生,而這部分質(zhì)子由于能量損失將不能透射物體。因而程序中同時加入G4ProtonInelasticProcess和G4hIonisation物理過程,統(tǒng)計出了透射質(zhì)子總個數(shù),減去單獨(dú)用G4ProtonInelasticProcess過程計算時能量未發(fā)生改變的直穿質(zhì)子的個數(shù),就得到了透射的次級質(zhì)子個數(shù)。對于次級質(zhì)子,圖8畫出了散射角分布圖,次級質(zhì)子個數(shù)較少,bin的劃分較大,為50 mrad。表1給出了產(chǎn)生的次級質(zhì)子個數(shù)和多次庫倫散射角均方根值大小。

    圖8 次級質(zhì)子的散射角分布Fig.8 Scattering angular distribution of secondary protons.

    三種面密度下散射角的大小要比MCS的散射角大很多,因而調(diào)整了橫坐標(biāo)的取值區(qū)間,分布函數(shù)都與高斯函數(shù)類似,具有中間多,周圍少的特點(diǎn)。此外,三種面密度下的分布沒有明顯的差別。

    3.4 彈性散射和準(zhǔn)彈性散射

    Geant4中物理模型設(shè)置為 G4LElastic,G4LElastic模型包括了彈性散射和非彈性散射過程。選取面密度為5 g·cm?2、10 g·cm?2和15 g·cm?2的鎢材料,分別入射 105個質(zhì)子,統(tǒng)計出發(fā)生彈性散射的比例(表1)。散射質(zhì)子的角度分布見圖9。

    表1 5 g·cm?2、10 g·cm?2和15 g·cm?2鎢中的散射角均方根值Table 1 RMS of scattering angle in 5 g·cm?2, 10 g·cm?2, 15 g·cm?2 tungsten.

    圖9 彈性散射的質(zhì)子角分布Fig.9 Elastic scattering angular distribution of protons.

    彈性碰撞的散射角也與高斯函數(shù)類似,具有中間多,周圍少的特點(diǎn),且三種面密度下的散射角分布沒有明顯的差別。

    3.5 討論

    質(zhì)子穿過物體過程中,能量損失、個數(shù)衰減、運(yùn)動方向發(fā)生偏轉(zhuǎn)。式(11)表明,通過測量能量損失可以給出電子的密度;在射程的前86%,能量損失與面密度近似成正比,且相同面密度的不同材料的能量損失大小不同。圖7表明,不發(fā)生非彈性碰撞的質(zhì)子個數(shù)隨面密度的增加呈指數(shù)衰減,且相同面密度的不同材料的個數(shù)損失不同。由表 1,引起質(zhì)子運(yùn)動方向發(fā)生偏轉(zhuǎn)的因素包括多次庫侖散射、非彈性碰撞產(chǎn)生的次級質(zhì)子以及與原子核的彈性散射。隨著面密度的增加,產(chǎn)生的次級質(zhì)子和發(fā)生彈性散射的質(zhì)子個數(shù)在逐漸增多;多次庫侖散射角逐漸增大,次級質(zhì)子的散射角逐漸減小,彈性散射角沒有明顯的變化規(guī)律。圖6表明,對于不同的材料,散射角的均方根值也不同。

    4 結(jié)語

    對于某種已知材料,能量損失的大小、個數(shù)損失的多少以及散射角的大小都可以反映材料的面密度。對于相同面密度的物體,不同材料的能量損失、個數(shù)損失和散射角的大小也不同,可以通過測量這些物理量來反映材料的組分信息。

    1 李家偉. 無損檢測手冊[M]. 北京: 機(jī)械工業(yè)出版社, 2002

    LI Jiawei. Non-destructive inspection handbook[M]. Beijing: China Machine Press, 2002

    2 Schulte R W, Bashkirov V, Loss Klock M C, et al. Density resolution of proton computed tomography[J]. Medical Physics, 2005, 24(4): 1035?1046

    3 樊明武. 用于醫(yī)學(xué)診斷和治療的質(zhì)子回旋加速器[J].中國工程科學(xué), 2000, 12(2): 9?15

    FAN Mingwu. Medical cyclotron used for diagnostic or therapy[J]. Engineering Science, 2000, 12(2): 9?15

    4 Morris C L, Ables E, Alrick K R, et al. Flash radiography with 24 GeV/c protons[J]. Journal of Applied Physics, 2011, 109(10): 104905_1?104905_10

    5 陳延偉. 中國散裂中子源(CSNS)[J]. 中國科學(xué)院院刊, 2011, 26(6): 726?729

    CHEN Yanwei. China spallation neutron source[J]. Bulletin of Chinese Academy of Sciences, 2011, 26(6): 726?729

    6 劉進(jìn), 章林文, 劉軍, 等. 快速高能質(zhì)子照相程序QMCPrad的研制[J]. 強(qiáng)激光與粒子束, 2012, 24(12): 2959?2964

    LIU Jin, Zhang Linwen, LIU Jun, et al. Development of code QMCPrad for fast high-energy proton radiography[J]. High Power Laser and Particle Beams, 2012, 24(12): 2959?2964

    7 Geant4 User’s Guide For Application Developers [EB/OL]. http://geant4.web.cern.ch/geant4/

    8 Physics Reference Manual[EB/OL]. http://geant4.web. cern.ch/geant4/

    猜你喜歡
    面密度庫侖原子核
    低面密度PE/PP雙組分紡粘非織造布的制備及其性能
    毛紡科技(2022年10期)2022-10-31 03:56:48
    1976年唐山強(qiáng)震群震后庫侖應(yīng)力演化及其與2020年古冶5.1級地震的關(guān)系
    地震研究(2021年1期)2021-04-13 01:04:46
    關(guān)于原子核結(jié)構(gòu)的討論
    纖維層面密度的光學(xué)新算法
    物質(zhì)構(gòu)成中的“一定”與“不一定”
    水刺產(chǎn)品面密度均勻度波動因素分析與自動控制
    基于粘彈庫侖應(yīng)力變化的后續(xù)最大地震震級估計及2008、2014年于田2次7.3級地震之間關(guān)系的討論
    中國地震(2015年1期)2015-11-08 11:11:18
    車身空腔用膨脹膠塊聲學(xué)性能影響因素分析
    粘接(2014年9期)2014-04-29 20:57:28
    一種周期庫侖作用勢優(yōu)化法的改進(jìn)
    計算物理(2014年1期)2014-03-11 17:01:03
    長程庫侖勢對高溫超導(dǎo)渦旋電荷的影響
    免费av观看视频| 最新中文字幕久久久久| 欧美中文综合在线视频| 精品日产1卡2卡| 黄色片一级片一级黄色片| 麻豆成人午夜福利视频| 色综合站精品国产| 性欧美人与动物交配| 超碰av人人做人人爽久久 | 校园春色视频在线观看| 亚洲欧美激情综合另类| 一级黄片播放器| 亚洲一区高清亚洲精品| 国产成人av教育| 国产真实乱freesex| 三级毛片av免费| 日本黄色片子视频| 国产精品亚洲av一区麻豆| 校园春色视频在线观看| xxx96com| 啦啦啦韩国在线观看视频| 午夜视频国产福利| 好看av亚洲va欧美ⅴa在| 老熟妇乱子伦视频在线观看| 欧美最新免费一区二区三区 | 久久精品国产99精品国产亚洲性色| 特大巨黑吊av在线直播| 99久国产av精品| 亚洲国产中文字幕在线视频| 蜜桃久久精品国产亚洲av| 国产aⅴ精品一区二区三区波| 国产色爽女视频免费观看| 熟女少妇亚洲综合色aaa.| 日韩欧美在线二视频| 老司机午夜福利在线观看视频| 欧美日韩一级在线毛片| 在线免费观看不下载黄p国产 | 久久国产精品影院| 老熟妇乱子伦视频在线观看| 亚洲国产高清在线一区二区三| 日本黄大片高清| 色在线成人网| 他把我摸到了高潮在线观看| 亚洲精品美女久久久久99蜜臀| 国产精品,欧美在线| 国产一区二区激情短视频| 国产亚洲欧美在线一区二区| 久久精品综合一区二区三区| 男女之事视频高清在线观看| 狠狠狠狠99中文字幕| 亚洲欧美一区二区三区黑人| 亚洲aⅴ乱码一区二区在线播放| 两个人视频免费观看高清| 欧美成人a在线观看| 久久精品91无色码中文字幕| 国产亚洲av嫩草精品影院| 亚洲国产色片| 午夜激情福利司机影院| 桃色一区二区三区在线观看| 色综合欧美亚洲国产小说| svipshipincom国产片| 国产在视频线在精品| 国产精品永久免费网站| 亚洲avbb在线观看| 亚洲第一电影网av| 免费在线观看日本一区| 亚洲国产日韩欧美精品在线观看 | 两个人视频免费观看高清| 夜夜爽天天搞| 国产亚洲精品久久久com| 熟女人妻精品中文字幕| 欧美黄色淫秽网站| 天堂影院成人在线观看| 中文字幕高清在线视频| 婷婷亚洲欧美| 老汉色∧v一级毛片| 2021天堂中文幕一二区在线观| 国产伦精品一区二区三区视频9 | 最近视频中文字幕2019在线8| 婷婷六月久久综合丁香| 国产精品1区2区在线观看.| 欧美高清成人免费视频www| 国产97色在线日韩免费| 给我免费播放毛片高清在线观看| 天美传媒精品一区二区| 久久精品国产清高在天天线| 国产爱豆传媒在线观看| 国产一区二区三区视频了| 亚洲av二区三区四区| 一个人看视频在线观看www免费 | 听说在线观看完整版免费高清| 一边摸一边抽搐一进一小说| 亚洲最大成人手机在线| 俺也久久电影网| 国产亚洲精品久久久com| 一级黄片播放器| 亚洲精品乱码久久久v下载方式 | 制服人妻中文乱码| svipshipincom国产片| 免费在线观看日本一区| 亚洲一区二区三区不卡视频| 国产毛片a区久久久久| 黄色成人免费大全| 91久久精品电影网| 亚洲欧美日韩无卡精品| 男女做爰动态图高潮gif福利片| 88av欧美| 亚洲性夜色夜夜综合| 中国美女看黄片| 欧美日韩一级在线毛片| 日韩欧美在线二视频| 窝窝影院91人妻| 久久精品亚洲精品国产色婷小说| 女同久久另类99精品国产91| 精品电影一区二区在线| 五月伊人婷婷丁香| 久久精品国产清高在天天线| 久久久国产精品麻豆| 亚洲av第一区精品v没综合| 男女视频在线观看网站免费| 一本久久中文字幕| 在线免费观看不下载黄p国产 | 午夜a级毛片| 美女黄网站色视频| a级毛片a级免费在线| 最近最新中文字幕大全免费视频| 18美女黄网站色大片免费观看| 男插女下体视频免费在线播放| 黄片大片在线免费观看| 国内精品久久久久久久电影| 精品久久久久久成人av| 亚洲中文日韩欧美视频| 欧美激情久久久久久爽电影| 国产欧美日韩精品一区二区| 啦啦啦韩国在线观看视频| 国产精品乱码一区二三区的特点| 我的老师免费观看完整版| 在线视频色国产色| 97超视频在线观看视频| 琪琪午夜伦伦电影理论片6080| 母亲3免费完整高清在线观看| 亚洲av免费高清在线观看| 99热6这里只有精品| 国产成人啪精品午夜网站| 精品一区二区三区av网在线观看| 99久久精品一区二区三区| 在线十欧美十亚洲十日本专区| 亚洲成人免费电影在线观看| 国产三级中文精品| 亚洲激情在线av| 国产一区二区激情短视频| 亚洲av一区综合| 成人性生交大片免费视频hd| 亚洲一区高清亚洲精品| 国产精品一区二区免费欧美| 熟妇人妻久久中文字幕3abv| 亚洲中文字幕日韩| 亚洲成人久久性| 搡老岳熟女国产| 88av欧美| 免费观看人在逋| 成人无遮挡网站| 国产精品亚洲av一区麻豆| 日韩大尺度精品在线看网址| 老熟妇乱子伦视频在线观看| 日韩人妻高清精品专区| 色噜噜av男人的天堂激情| 亚洲熟妇中文字幕五十中出| 男女视频在线观看网站免费| xxxwww97欧美| 99久久综合精品五月天人人| 亚洲人成网站高清观看| 特级一级黄色大片| 老熟妇仑乱视频hdxx| 久久久久久久久大av| 亚洲欧美一区二区三区黑人| 别揉我奶头~嗯~啊~动态视频| 亚洲精品456在线播放app | or卡值多少钱| 午夜福利在线观看吧| 久久草成人影院| 波多野结衣巨乳人妻| 熟女电影av网| eeuss影院久久| 超碰av人人做人人爽久久 | 麻豆一二三区av精品| 亚洲人成网站在线播放欧美日韩| 国产男靠女视频免费网站| 久久99热这里只有精品18| 亚洲av一区综合| 在线观看免费午夜福利视频| 人妻久久中文字幕网| 日韩欧美在线乱码| 国产99白浆流出| 丰满的人妻完整版| 日本免费a在线| 无人区码免费观看不卡| 国产aⅴ精品一区二区三区波| 一个人看视频在线观看www免费 | 美女 人体艺术 gogo| 日本一本二区三区精品| 在线观看66精品国产| 国产一区在线观看成人免费| 嫩草影院精品99| 亚洲国产色片| 国产高清三级在线| 国产精品98久久久久久宅男小说| 亚洲成人中文字幕在线播放| 国产精品爽爽va在线观看网站| 99精品久久久久人妻精品| 免费在线观看日本一区| 国产老妇女一区| 亚洲片人在线观看| 欧美日韩中文字幕国产精品一区二区三区| 国产亚洲av嫩草精品影院| 制服人妻中文乱码| 婷婷精品国产亚洲av| 精华霜和精华液先用哪个| 国产精品女同一区二区软件 | 国产不卡一卡二| а√天堂www在线а√下载| 午夜激情福利司机影院| 日日夜夜操网爽| 国产99白浆流出| 深夜精品福利| 亚洲五月婷婷丁香| 制服人妻中文乱码| 99视频精品全部免费 在线| 中文字幕人成人乱码亚洲影| 国内毛片毛片毛片毛片毛片| 成人性生交大片免费视频hd| 中文字幕人妻熟人妻熟丝袜美 | 亚洲av第一区精品v没综合| 美女高潮的动态| 中文资源天堂在线| 大型黄色视频在线免费观看| 国产精品三级大全| 亚洲av免费高清在线观看| 国产精品一区二区免费欧美| 亚洲中文字幕日韩| 丰满人妻一区二区三区视频av | 人妻丰满熟妇av一区二区三区| 免费av不卡在线播放| 五月伊人婷婷丁香| 两个人视频免费观看高清| 婷婷丁香在线五月| 国产亚洲精品一区二区www| 久久久久亚洲av毛片大全| 婷婷六月久久综合丁香| 91字幕亚洲| 国产欧美日韩一区二区三| 91久久精品国产一区二区成人 | 免费观看人在逋| 男女做爰动态图高潮gif福利片| 在线观看美女被高潮喷水网站 | 亚洲精品色激情综合| 在线观看午夜福利视频| 欧美极品一区二区三区四区| 欧美日韩瑟瑟在线播放| 在线视频色国产色| 国产伦人伦偷精品视频| 免费在线观看亚洲国产| 国产亚洲精品综合一区在线观看| 免费看a级黄色片| 国产精品精品国产色婷婷| 一级作爱视频免费观看| 精品久久久久久久久久免费视频| 熟女电影av网| 亚洲av美国av| 黄片大片在线免费观看| 免费在线观看日本一区| 日韩人妻高清精品专区| 国产av不卡久久| 精品熟女少妇八av免费久了| 一级毛片女人18水好多| 国内精品美女久久久久久| 午夜免费观看网址| 日韩高清综合在线| 久久久精品大字幕| 日本免费一区二区三区高清不卡| 成人特级av手机在线观看| aaaaa片日本免费| 欧美乱码精品一区二区三区| 亚洲成人精品中文字幕电影| 嫩草影院精品99| 精品无人区乱码1区二区| 欧美黄色淫秽网站| 高清在线国产一区| 午夜两性在线视频| 国产精品,欧美在线| 性色avwww在线观看| 欧美大码av| 日韩欧美 国产精品| 听说在线观看完整版免费高清| 成人鲁丝片一二三区免费| 两性午夜刺激爽爽歪歪视频在线观看| 成人三级黄色视频| 精品免费久久久久久久清纯| 蜜桃久久精品国产亚洲av| 女同久久另类99精品国产91| 国产精品一区二区三区四区久久| 真人一进一出gif抽搐免费| 99久久精品一区二区三区| 好男人电影高清在线观看| 亚洲狠狠婷婷综合久久图片| 亚洲人成伊人成综合网2020| 91九色精品人成在线观看| 搡老岳熟女国产| 少妇裸体淫交视频免费看高清| 91麻豆av在线| 精品福利观看| 宅男免费午夜| 日韩亚洲欧美综合| 麻豆一二三区av精品| 亚洲乱码一区二区免费版| 搡老妇女老女人老熟妇| 亚洲人成网站在线播| 久久香蕉精品热| 国产精品一区二区三区四区久久| 老司机午夜十八禁免费视频| 熟妇人妻久久中文字幕3abv| 国产主播在线观看一区二区| 色精品久久人妻99蜜桃| 日韩大尺度精品在线看网址| 岛国在线观看网站| 精品久久久久久,| 我要搜黄色片| 久久午夜亚洲精品久久| 免费人成在线观看视频色| 国产精品影院久久| 69人妻影院| www.熟女人妻精品国产| 亚洲成人免费电影在线观看| 欧美一区二区国产精品久久精品| 精品无人区乱码1区二区| 一本综合久久免费| 精华霜和精华液先用哪个| 欧美在线一区亚洲| av黄色大香蕉| 黄色成人免费大全| 99精品欧美一区二区三区四区| 亚洲欧美日韩高清专用| 国产高清三级在线| 国产高潮美女av| 日日干狠狠操夜夜爽| av黄色大香蕉| 一区二区三区高清视频在线| 在线观看美女被高潮喷水网站 | 国产一区二区三区在线臀色熟女| 成人鲁丝片一二三区免费| 欧美最新免费一区二区三区 | 国内精品美女久久久久久| 90打野战视频偷拍视频| 99久久综合精品五月天人人| 成人特级黄色片久久久久久久| 国产一区二区激情短视频| 女人被狂操c到高潮| 婷婷亚洲欧美| 99在线人妻在线中文字幕| 国产视频内射| 欧美日韩综合久久久久久 | 中出人妻视频一区二区| 又紧又爽又黄一区二区| 亚洲无线在线观看| 精品99又大又爽又粗少妇毛片 | 一级黄片播放器| 少妇高潮的动态图| 久久香蕉国产精品| 熟女少妇亚洲综合色aaa.| 国产精品98久久久久久宅男小说| 国产一区二区在线av高清观看| 午夜激情福利司机影院| 免费看光身美女| 在线免费观看的www视频| 黄色片一级片一级黄色片| 国产探花极品一区二区| 99久久精品热视频| 美女 人体艺术 gogo| 日韩国内少妇激情av| www.999成人在线观看| 国产成人欧美在线观看| 精品一区二区三区视频在线观看免费| 深爱激情五月婷婷| 国产高清三级在线| 免费看a级黄色片| 中国美女看黄片| 免费观看精品视频网站| 午夜福利18| 久久久久国内视频| 国产熟女xx| 岛国视频午夜一区免费看| 成年免费大片在线观看| 日本黄色视频三级网站网址| 男女之事视频高清在线观看| 国产探花在线观看一区二区| 欧美bdsm另类| 一进一出抽搐动态| 精品不卡国产一区二区三区| 99国产综合亚洲精品| 亚洲在线自拍视频| 成人三级黄色视频| 男人和女人高潮做爰伦理| 香蕉久久夜色| 88av欧美| 法律面前人人平等表现在哪些方面| 亚洲精品亚洲一区二区| 亚洲五月天丁香| 午夜免费男女啪啪视频观看 | 极品教师在线免费播放| 中文在线观看免费www的网站| 欧美日本亚洲视频在线播放| 中亚洲国语对白在线视频| 9191精品国产免费久久| 黄片小视频在线播放| 国产v大片淫在线免费观看| 日本免费一区二区三区高清不卡| 69人妻影院| 国产激情偷乱视频一区二区| 午夜福利欧美成人| 美女 人体艺术 gogo| 午夜两性在线视频| 亚洲av第一区精品v没综合| 啦啦啦韩国在线观看视频| 国产伦精品一区二区三区视频9 | 成年女人毛片免费观看观看9| 人妻夜夜爽99麻豆av| 最近最新免费中文字幕在线| 午夜日韩欧美国产| 亚洲精品色激情综合| 人妻久久中文字幕网| 此物有八面人人有两片| 日韩欧美精品免费久久 | 婷婷亚洲欧美| 精品久久久久久久毛片微露脸| 国产高清三级在线| 国产男靠女视频免费网站| 日本一二三区视频观看| 搞女人的毛片| 少妇的逼好多水| 婷婷丁香在线五月| 神马国产精品三级电影在线观看| 中国美女看黄片| 午夜视频国产福利| 高清日韩中文字幕在线| 99久国产av精品| 大型黄色视频在线免费观看| 亚洲五月天丁香| 少妇熟女aⅴ在线视频| 久久精品综合一区二区三区| 日本黄色视频三级网站网址| 日韩欧美一区二区三区在线观看| 欧美色欧美亚洲另类二区| 久久香蕉精品热| 一个人免费在线观看电影| 国内精品久久久久精免费| 国内揄拍国产精品人妻在线| av片东京热男人的天堂| 欧美不卡视频在线免费观看| 亚洲狠狠婷婷综合久久图片| 欧美一级a爱片免费观看看| 色噜噜av男人的天堂激情| 狂野欧美白嫩少妇大欣赏| 国产99白浆流出| 欧美成人a在线观看| 国产成+人综合+亚洲专区| 麻豆国产97在线/欧美| 尤物成人国产欧美一区二区三区| av国产免费在线观看| 欧美成狂野欧美在线观看| 黄色日韩在线| 性色av乱码一区二区三区2| 中文字幕高清在线视频| 免费高清视频大片| 制服人妻中文乱码| 久久精品国产亚洲av香蕉五月| 亚洲av不卡在线观看| 最好的美女福利视频网| 亚洲人成伊人成综合网2020| 亚洲中文字幕一区二区三区有码在线看| 91久久精品国产一区二区成人 | 国产精品国产高清国产av| 99久久精品热视频| 一卡2卡三卡四卡精品乱码亚洲| 村上凉子中文字幕在线| 亚洲人成伊人成综合网2020| 欧美日韩黄片免| 亚洲最大成人手机在线| 网址你懂的国产日韩在线| 国产黄a三级三级三级人| 男女之事视频高清在线观看| 一级毛片高清免费大全| 国产高清视频在线观看网站| 国产精品一区二区免费欧美| 欧美午夜高清在线| 亚洲人成网站在线播| www日本黄色视频网| 99热6这里只有精品| 国产精品综合久久久久久久免费| 国产亚洲精品久久久com| 成人高潮视频无遮挡免费网站| 深爱激情五月婷婷| 一级毛片高清免费大全| 国产欧美日韩精品亚洲av| 天堂网av新在线| 亚洲精品色激情综合| 中文字幕人妻熟人妻熟丝袜美 | 老司机午夜福利在线观看视频| 日韩人妻高清精品专区| 午夜精品在线福利| 亚洲 欧美 日韩 在线 免费| 日本一二三区视频观看| 少妇高潮的动态图| 超碰av人人做人人爽久久 | 色噜噜av男人的天堂激情| 99久国产av精品| 成人三级黄色视频| 两人在一起打扑克的视频| 国产精品久久久久久亚洲av鲁大| 成人性生交大片免费视频hd| 午夜福利在线观看免费完整高清在 | 国产成人av教育| 大型黄色视频在线免费观看| 久久久国产成人免费| 国产高清视频在线播放一区| 麻豆成人午夜福利视频| 三级国产精品欧美在线观看| 男人和女人高潮做爰伦理| 三级国产精品欧美在线观看| 亚洲精品久久国产高清桃花| 久久久久国产精品人妻aⅴ院| 两个人视频免费观看高清| 亚洲最大成人中文| 91在线精品国自产拍蜜月 | 欧美性猛交╳xxx乱大交人| 亚洲午夜理论影院| 久久久精品欧美日韩精品| 欧美乱码精品一区二区三区| 叶爱在线成人免费视频播放| 国产精品99久久久久久久久| 午夜亚洲福利在线播放| 亚洲一区二区三区不卡视频| 亚洲性夜色夜夜综合| 熟妇人妻久久中文字幕3abv| 观看美女的网站| 亚洲国产精品合色在线| 69av精品久久久久久| 日韩精品中文字幕看吧| 国产高清有码在线观看视频| 一个人看视频在线观看www免费 | 韩国av一区二区三区四区| 欧美一区二区亚洲| 老司机深夜福利视频在线观看| 狠狠狠狠99中文字幕| 国语自产精品视频在线第100页| 中亚洲国语对白在线视频| 久久国产精品影院| 毛片女人毛片| 两个人的视频大全免费| 99精品久久久久人妻精品| 91字幕亚洲| 国产亚洲欧美98| 精品久久久久久,| 99久国产av精品| 少妇丰满av| 黄色丝袜av网址大全| 国产精品久久久人人做人人爽| 国产av一区在线观看免费| 99久久精品热视频| 亚洲人成网站高清观看| 亚洲精品国产精品久久久不卡| 少妇的逼水好多| 99热精品在线国产| 亚洲中文日韩欧美视频| 岛国视频午夜一区免费看| 一级毛片女人18水好多| 精品久久久久久久毛片微露脸| 久久精品国产亚洲av涩爱 | 一级a爱片免费观看的视频| 亚洲av不卡在线观看| 免费电影在线观看免费观看| 国产aⅴ精品一区二区三区波| 岛国视频午夜一区免费看| 精品一区二区三区av网在线观看| 亚洲精品美女久久久久99蜜臀| 99热精品在线国产| 99国产综合亚洲精品| 欧美绝顶高潮抽搐喷水| 嫩草影院入口| 国内精品美女久久久久久| 无人区码免费观看不卡| 欧美又色又爽又黄视频| 欧美极品一区二区三区四区| 99久久精品热视频| av黄色大香蕉| 亚洲国产欧美人成| 国产成+人综合+亚洲专区| 日韩av在线大香蕉| 亚洲精品在线美女| 热99re8久久精品国产| 久久国产乱子伦精品免费另类| 亚洲人成网站在线播| 91麻豆精品激情在线观看国产| 夜夜躁狠狠躁天天躁| 天堂√8在线中文| 韩国av一区二区三区四区| 18禁美女被吸乳视频| 三级男女做爰猛烈吃奶摸视频| 禁无遮挡网站| 精品人妻1区二区| 亚洲男人的天堂狠狠| 波多野结衣高清作品| 亚洲国产日韩欧美精品在线观看 | 欧美+日韩+精品| 日韩av在线大香蕉| 精品国产亚洲在线|