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

    上海光源光束線站BL18U軔致輻射劑量的計算與測量

    2014-03-06 09:07:36徐加強夏曉彬盛尹祥子秦銳芳費平安劉嗣鑫
    核技術 2014年7期
    關鍵詞:棚屋劑量率模擬計算

    徐加強 夏曉彬 盛尹祥子 秦銳芳 張 冬 費平安 姚 明 劉嗣鑫

    上海光源光束線站BL18U軔致輻射劑量的計算與測量

    徐加強1,2夏曉彬1盛尹祥子1秦銳芳1張 冬1費平安1姚 明1劉嗣鑫1

    1(中國科學院上海應用物理研究所 嘉定園區(qū) 上海 201800)
    2(中國科學院大學 北京 100049)

    軔致輻射是第三代先進同步輻射光源光束線站輻射屏蔽的重點。上海光源(Shanghai Synchrotron Radiation Facility, SSRF)是國際上正在運行的第三代先進同步輻射光源之一。上海光源蛋白質微晶結構光束線(BL18U)線站是正在試運行中的光束線站。本文采用FLUKA軟件模擬計算了蛋白質微晶結構光束線(BL18U)線站因軔致輻射的散射和光中子引起的輻射劑量率分布。分析了因狹縫開孔尺寸、儲存環(huán)內(nèi)流強的變化所引起的輻射劑量變化。用高靈敏度的光子和中子探測器測量了BL18U光學棚屋外的光子和光中子劑量率。測量結果表明模擬計算的可靠性。本研究采用的模擬和測量方法可用于其它線站的輻射劑量水平評估,并為上海光源后續(xù)線站的屏蔽設計提供參考。

    同步輻射,光束線站,蒙特卡羅方法,F(xiàn)LUKA,軔致輻射,輻射測量

    高能電子在同步光源儲存環(huán)內(nèi)旋轉時,會與真空管內(nèi)殘余氣體分子的庫侖場發(fā)生非彈性碰撞,產(chǎn)生氣體軔致輻射。在儲存環(huán)偏轉磁鐵之間的直線節(jié)段,產(chǎn)生的氣體軔致輻射將沿電子的軌道不斷增強,形成張角為m0c2/E0大小尖銳前沖,光子能譜范圍從0-E0連續(xù)呈1/E0分布的軔致輻射源項,這里m0c2為靜止電子能量,E0為初始電子能量[1-2]。氣體軔致輻射伴隨插入件產(chǎn)生的同步輻射光一起傳輸?shù)较掠喂馐€站,并與光束線站上的各類光學器件相互作用后形成電磁級聯(lián)簇射。光子能量超過器件光核反應閾值(大于5-15 MeV)的軔致輻射,通過巨核偶極共振衰變、準氘核光致蛻變和核內(nèi)崩潰三種機制產(chǎn)生光核反應中子(以下簡稱“光中子”)[3]。軔致輻射的散射光子和光中子是光束線站屏蔽的主要輻射源項[4-5]。

    光束線站的輻射源項和屏蔽設計已有許多研究成果,如軔致輻射、光中子能譜、輻射劑量的測量與計算[3,6]、次級軔致輻射的屏蔽計算等[4-5,7-9]。這些研究都局限于測量和計算某個散射靶,而缺少對整條光束線站的輻射劑量計算和測量。本研究采用FLUKA[10-11]軟件開展蒙特卡羅模擬計算氣體軔致輻射在光束線站上產(chǎn)生的輻射劑量。FLUKA軟件能夠模擬計算包括軔致輻射與各類光學器件相互作用產(chǎn)生的電磁級聯(lián)簇射,以及通過光核反應產(chǎn)生的光中子,并可以將注量轉換為周圍劑量當量[12]得到各類粒子所致的輻射劑量。研究對象選擇上海光源(Shanghai Synchrotron Radiation Facility, SSRF)新建的典型插入件光束線站——蛋白質微晶結構光束線(BL18U)。為了模擬計算結果與測量結果比較,對BL18U棚屋外的輻射劑量分布進行了實驗測量。

    1 BL18U線站布局簡介

    上海光源的光束線站一般由一個光學棚屋和一個實驗棚屋組成,BL18U的布局特點是具有比較長的光學棚屋,以保證入射光高通量的前提下,照射至樣品處的光斑尺寸盡量小。為減少建設代價,在其布局上分為兩個光學棚屋。圖1給出了沿光束線BL18U光學棚屋布局示意圖,兩個光學棚屋分別為圖1中所示的OE1和OE2。這兩個光學棚屋之間通過一段7 m長的光束線輸運管道連接。OE1左側為儲存環(huán)屏蔽墻,右側為實驗大廳。光束線輸運管管道和OE2兩側部分空間布置了機柜,其余空間也是實驗大廳。

    圖1給出了BL18U光束線上對軔致輻射有較大影響的關鍵光學設備。這些設備包括:光學前端區(qū)的安全光閘(Safety Shutter)、在OE1內(nèi)的水冷狹縫(Slits)和準直吸收器1(Collimator 1)、在OE2內(nèi)的單色儀(Monochromator)和吸收器2以及安全光閘2 (Safety Shutter 2)。

    圖1 BL18U線站平面布局及監(jiān)測點布置示意圖Fig.1 Layout of BL18U plane view and arrangement of measurement points.

    2 劑量率的模擬計算

    2.1 計算模型的建立

    采用FLUKA軟件模擬計算了氣體軔致輻射在直線節(jié)內(nèi)產(chǎn)生后,在光束線站BL18U光學棚屋內(nèi)的傳輸過程。計算模型根據(jù)BL18U布局圖(見圖1)和BL18U的射線追跡圖(見圖2)建立。如圖2所示,3.5 GeV能量的電子在長12.5 m的標準直線節(jié)真空管內(nèi)傳輸,在直線節(jié)末端電子通過1.27 T的偏轉磁場轉移到儲存環(huán)的下一個單元;產(chǎn)生的氣體軔致輻射通過打開的前端安全光閘和儲存環(huán)的鋸齒屏蔽墻開孔,進入到BL18U的光學棚屋;然后與OE1內(nèi)的水冷狹縫和準直吸收器相互作用,透過狹縫和準直吸收器開孔的軔致輻射又通過光束輸運管道進入到OE2內(nèi),與單色儀作用后最終損失在最后一個準直吸收器上。氣體軔致輻射被器件散射后入射到光學棚屋和管道上,產(chǎn)生的光子和光中子經(jīng)棚屋屏蔽層的衰減,最后到達棚屋外。BL18U插入件所在直線節(jié)處的真空度平均值約99.32 nPa,其真空管內(nèi)殘余氣體有效Z值為3.42,其主要成分依據(jù)真空測試數(shù)據(jù)如表1所示。

    表1 SSRF儲存環(huán)直線節(jié)真空管內(nèi)的殘余氣體成分Table 1 Compositions of the residual gas in the electron storage ring of SSRF.

    圖2 BL18U軔致輻射射線追跡示意圖Fig.2 Bremsstrahlung ray trace of BL18U.

    2.2 劑量分布

    用FLUKA軟件模擬了BL18U線站在儲存環(huán)工作在恒流注入模式,束流強度為220mA和300 mA時,氣體軔致輻射與各類光學器件相互作用后產(chǎn)生的所有粒子在光學棚屋區(qū)域的劑量率分布。束流強度為220 mA時的結果如圖3所示。從圖3可看出輻射劑量最大點集中在水冷狹縫和單色儀周圍。

    圖3 BL18U線站光學棚屋的劑量率分布Fig.3 Dose rate distribution at BL18U optical enclosure.

    光束線站的軔致輻射劑量會因線站狹縫開口大小、儲存環(huán)束流和直線節(jié)真空狀態(tài)變化。沿圖1所示的1-20個測量參考點連線(與束流中心同高度1.3 m),計算得到不同狹縫開孔下的光子和光中子的劑量率分布如圖4所示。BL18U在調(diào)試時狹縫開口為0.4mm×0.4mm,在正常工作狀態(tài)時則為3.5mm×8.0mm。從圖4可以看到,狹縫開孔較小時,軔致輻射主要損失在水冷狹縫所在的OE1內(nèi),狹縫開孔較大時,OE2周圍的劑量率開始增大。

    圖4 不同狹縫開口大小,棚屋外的光子(a)和光中子(b)劑量率分布Fig.4 Dose rate distribution of photons (a) and photoneutrons (b) out of optical enclosure for different size of slits.

    儲存環(huán)流強為300 mA,直線節(jié)真空度為133.32nPa時,氣體軔致輻射引起的光學棚屋外輻射劑量率分布如圖5所示。從計算結果可知,經(jīng)過棚屋屏蔽層的衰減吸收,光子最大劑量率為0.3μSv·h-1,光中子的最大劑量率為0.55 μSv·h-1。水冷狹縫中心距離鋸齒墻約2 m,而棚屋右側墻距離束流中心1.2 m。棚屋外的中子劑量率的最大點對應在圖1中的○2測量點,是水冷狹縫側向約90°的位置處;而光子劑量率最大點在在圖1中的○3號測量點,約在水冷狹縫靶點前方側向的45°位置處。

    圖5 儲存環(huán)300 mA流強,直線節(jié)133.32 nPa真空條件下棚屋外的劑量率分布Fig.5 Dose rate distribution at the outer OE1 (In case of beam current of 300 mA in storage ring and vacuum of 133.32 nPa in straight section).

    2.3 光子和光中子能譜

    通過模擬計算統(tǒng)計棚屋內(nèi)距離水冷狹縫側向0.5 m處和棚屋外距離水冷狹縫側向1.2 m處的光子和光中子能譜如圖6和圖7所示。由于軔致輻射能譜近似呈1/E分布,經(jīng)狹縫、單色器和準直吸收器散射后的光子在經(jīng)過16 mm以上的鉛屏蔽墻后,低能部分被衰減吸收,能量最可幾分布在1 MeV左右,見圖6。光核反應中,巨共振光核反應的截面較大,因此散射到側向的中子主要來自于各項同性的巨共振中子,平均能量約2 MeV,如圖7所示。

    圖6 OE1側向屏蔽墻內(nèi)外的光子能譜Fig.6 Photon spectrum at the inner and outer OE1.

    圖7 OE1內(nèi)外光中子能譜 (a) 所有中子,(b) 巨共振中子Fig.7 Photoneutron spectra at the inner and outer OE1. (a) All neutrons, (b) Giant resonance neutrons

    3 劑量率的測量

    棚屋外輻射劑量的測量選用高靈敏度環(huán)境輻射水平γ監(jiān)測儀(Area Gamma Monitor, AGM5-12-001)和高靈敏度環(huán)境輻射水平中子監(jiān)測儀(Environmental Neutron Monitor, ENM5-12-001)。γ監(jiān)測儀AGM采用球形高氣壓電離室作為γ探測器,其直徑為250 mm,不銹鋼壁厚1.8 mm,容積8.5 L,內(nèi)充3.333 MPa純氬氣。其靈敏度好于0.15nSv·Plus-1,測量范圍0.01-100 μSv·h-1,能量響應 50keV-3MeV。對照圖6所示棚屋外的光子能譜,AGM的能量響應范圍涵蓋了光子注量的主要區(qū)域,未涵蓋住的部分占光子總注量的20%左右。

    高靈敏度中子監(jiān)測儀ENM的BF3正比計數(shù)管為?50mm×350mm的柱體,充氣壓力79.992 kPa。能量響應范圍在0.25 eV-16 MeV,能量響應范圍能夠涵蓋如圖7(b)所示的巨共振中子能譜。其靈敏度較高,可以監(jiān)測低水平的中子劑量,探測下限在1nSv·h-1。該探測器對650 μSv·h-1以下的γ光子劑量不靈敏,總的不確定度小于15%[13]。

    測量時,儲存環(huán)運行在恒流注入模式,電子量3.5 GeV,儲存環(huán)流強220 mA,束流平均壽命16 h。BL18U光束線站工作在正常工作狀態(tài),此時水冷狹縫開孔為3.5 mm×8.0 mm。在光學棚屋外側,選擇了如圖1所示的20個測量點,用AGM和ENM進行輻射劑量的測量。測量點儀器中心高度1.3 m,與束流中心同高,測量點之間間距一般在1 m。每個地點測量約5 min,測量結果取平均值并扣除本底,其中光子本底劑量率平均為0.069 μSv·h-1,中子本底平均劑量率為0.010 μSv·h-1。

    4 測量結果比較分析

    用ENM測量OE1、輸運管道和OE2外的光核反應產(chǎn)生的中子劑量率。測量結果如圖8(a)所示,橫坐標為每個測量點距離鋸齒墻的距離。這些測量結果與相同位置處模擬計算得到的周圍當量劑量一起表示在圖8(a)中。如圖8(a)所示,ENM的測量結果與只統(tǒng)計巨共振中子的計算結果符合較好。其中測量得到的光中子的最大劑量率,在扣除本底后為0.24μSv·h-1。模擬計算中只統(tǒng)計巨共振中子時,棚屋外最大劑量率為0.23 μSv·h-1,兩者符合得很好。用AGM測量光學棚屋外20個測量點處光子的劑量率結果如圖8(b)所示。由圖8(b)看出,棚屋外光子的劑量水平較低,計算得到棚屋外最大劑量率為0.14 μSv·h-1,接近實驗大廳內(nèi)的輻射本底水平。受輻射本底漲落的影響,測量結果誤差較大。在63%的誤差范圍內(nèi),光子劑量率的模擬計算結果和實驗測量結果符合。

    圖8 光學棚屋外輻射劑量率計算與測量結果的比較 (a) 中子,(b) 光中子Fig.8 Dose rate comparisons of calculated and measured results at outside optical enclosure neutrons (a) and photons (b).

    5 結語

    本文用蒙特卡羅模擬軟件FLUKA,通過對上海光源光束線站BL18U的布局和射線追跡分析建立計算模型,模擬計算了上海光源儲存環(huán)運行在不同狹縫開口、束流強度和真空條件下該線站的氣體軔致輻射劑量率分布、光子和光中子能譜、以及線站光學棚屋外的輻射劑量。比較分析了儲存環(huán)工作在220 mA恒流注入模式下BL18U光學棚屋和輸運管道外的輻射劑量的計算與實驗測量結果。光中子劑量率的計算結果與測量結果符合得很好,而光子劑量率在63%的誤差范圍內(nèi)一致。測量結果表明模擬計算的可靠性。BL18U是上海光源16個標準直線節(jié)插入件光束線站之一,本研究采用的模擬和測量方法可用于其它線站的輻射劑量水平評估,并為上海光源后續(xù)線站的屏蔽設計提供參考。

    1 Rindi A. Gas bremsstrahlung from electron storage rings[J]. Health Physics, 1982, 42(2): 187-193

    2 Pisharody M, Job P K, Magill S, et al. Measurement of gas bremsstrahlung from electron storage rings[J]. Nuclear Instruments and Methods in Physics Research A, 1997, 401(2-3): 442-462

    3 Pisharody M, Semones E, Job P K. Dose measurements of bremsstrahlung-produced neutrons from thick targets[J]. Nuclear Instruments and Methods in Physics Research A, 1999, 430: 542-558

    4 Asai J, Hirayama H. Radiation shielding from secondary gas bremsstrahlung for the first optics enclosure at the Canadian Light Source[J]. Nuclear Instruments and Methods in Physics Research A, 2004, 527(3): 264-272

    5 Asai J, Hirayama H. Radiological design study of first optics enclosure against secondary gas bremsstrahlung for the Canadian Light Source[J]. Nuclear Instruments and Methods in Physics Research A, 2004, 517(1-3): 9-18

    6 Sheu R J, Wang J P, Sheu R D, et al. Gas bremsstrahlung and induced photoneutrons in the NSRRC's electron storage ring: a comparison of measurements and simulations[J]. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, 2004, 217(4): 555-563

    7 Asai J. Radiation shielding considerations against gas bremsstrahlung for the transverse wall at the Canadian Light Source[J]. Nuclear Instruments and Methods in Physics Research Section B, 2006, 251(2): 461-469

    8 Asai J, Hirayama H. Shielding design study against gas bremsstrahlung for the CLS first optics enclosures[J]. Radiation Measurements, 2006, 41: 186-193

    9 Job P K, Casey W R. Shielding calculations for the National Synchrotron Light Source-II experimental beamlines[J]. Nuclear Instruments and Methods in Physics Research A, 2013, 698: 196-201

    10 Fasso A, Ferrari A, Ranft J, et al. FLUKA: status and prospective for hadronic applications[M]. Kling A, Barao F, Nakagawa M, et al. ed. Proceedings of the Monte Carlo 2000 Conference, Lisbon, Oct 23-26, 2000, Berlin: Springer-Verlag, 2001: 955-960

    11 Fasso A, Ferrari A, Sala P R. Electron-photon transport in FLUKA: status[M]. Kling A, Barao F, Nakagawa M, et al. ed. Proceedings of the Monte Carlo 2000 Conference, Lisbon, Oct 23-26, 2000, Berlin: Springer-Verlag, 2001: 159-164

    12 Pelliccioni M. Overview of fluence-to-effective dose and fluence-to-ambient dose equivalent conversion coefficients for high energy radiation calculated using the FLUKA code[J]. Radiation Protection Dosimetry, 2000, 88(4): 279-297

    13 Li J P, Tang Y L, Li S D, et al. Neutron energy response of a modified Andersson-Braun Rem Counter[J]. Radiation Protection Dosimetry, 1996, 67(3): 179-185

    CLC TL72, R142

    Dose calculation and measurement for bremsstrahlung at BL18U beamline of SSRF

    XU Jiaqiang1,2XIA Xiaobin1SHENG Yinxiangzi1QIN Ruifang1ZHANG Dong1FEI Ping’an1YAO Ming1LIU Sixin1
    1(Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Jiading Campus, Shanghai 201800, China)
    2(University of Chinese Academy of Sciences, Beijing 100049, China)

    Background: Gas bremsstrahlung is one of the most important radiation sources that needs to be taken into consideration for shielding design of beamlines at the third generation synchrotron radiation light source. Shanghai Synchrotron Radiation Facility (SSRF) is one of the third generation synchrotron radiation light source in the world. The Protein Micro-crystallography Beamline (BL18U) is one of the commissioning beamlines and is a representative insertion device beamline at SSRF. Purpose: Estimation of radiation dose induced by scattering bremsstrahlung and photoneutrons at BL18U. Methods: Dose rate distribution induced by scattering bremsstrahlung and photoneutrons at BL18U are performed by Monte Carlo simulation code FLUKA. The radiation dose was analyzed with the variation of slits size, beam current at storage ring and the vacuum. Dose rate of photons and photoneutrons at the outside of the optical enclosure of BL18U were measured by using high sensitivity photon and neutron monitors. Results: The measurement results show that the reliability of the simulation. Conclusion: The simulation and measurement methods presented in this study can be applied to evaluate the dose rate level of other beamline stations at SSRF, and provide references to the shielding design for the following beamlines at SSRF in the near future.

    Synchrotron radiation, Beamline station, Monte Carlo method, FLUKA, Bremsstrahlung, Radiation measurement

    TL72,R142

    10.11889/j.0253-3219.2014.hjs.37.070101

    徐加強,男,1974年出生,2001年于上海大學獲碩士學位,現(xiàn)為博士研究生,研究領域為核技術及應用

    夏曉彬,E-mail: xiaxiaobin@sinap.ac.cn

    2014-04-16,

    2014-05-12

    猜你喜歡
    棚屋劑量率模擬計算
    R1234ze PVTx熱物性模擬計算
    能源工程(2022年1期)2022-03-29 01:06:26
    尋找荒野棚屋
    看世界(2022年2期)2022-02-14 22:08:03
    澳一家九口與世隔絕住棚屋過田園生活
    奇聞怪事(2020年4期)2020-12-21 03:57:55
    甲狀腺乳頭狀癌患者術后首次131Ⅰ治療后輻射劑量率的影響因素及出院時間的探討
    擠出發(fā)泡片材褶皺分析及模擬計算
    中國塑料(2016年9期)2016-06-13 03:18:54
    X線照射劑量率對A549肺癌細胞周期的影響
    癌癥進展(2016年10期)2016-03-20 13:15:43
    ITER極向場線圈PF4維修區(qū)域停堆劑量率分析
    如何有效調(diào)整ELEKTA PRECISE加速器的劑量率
    實際發(fā)射工況下底排藥柱結構完整性的模擬計算
    火炸藥學報(2014年3期)2014-03-20 13:17:43
    馬克斯·穆恩的瘋狂割草機的秘密
    美女被艹到高潮喷水动态| 五月伊人婷婷丁香| 人人妻人人看人人澡| 好男人在线观看高清免费视频| 亚洲美女视频黄频| 精品午夜福利视频在线观看一区| 国产不卡一卡二| 国产熟女xx| 麻豆国产av国片精品| 日日夜夜操网爽| 在线观看美女被高潮喷水网站 | 精品久久久久久,| 成熟少妇高潮喷水视频| 午夜福利免费观看在线| 国产一级毛片七仙女欲春2| av黄色大香蕉| 国产一区二区激情短视频| 在线十欧美十亚洲十日本专区| 夜夜夜夜夜久久久久| 九色成人免费人妻av| 夜夜看夜夜爽夜夜摸| 午夜久久久久精精品| 日本a在线网址| 欧美区成人在线视频| 国产成人欧美在线观看| 亚洲精品国产精品久久久不卡| 亚洲专区国产一区二区| 日韩欧美在线二视频| 亚洲精品久久国产高清桃花| 国产主播在线观看一区二区| 国产精品99久久99久久久不卡| xxx96com| 久9热在线精品视频| 国产一区二区激情短视频| 日本免费一区二区三区高清不卡| 午夜免费男女啪啪视频观看 | 在线观看av片永久免费下载| 精品久久久久久久久久免费视频| 18禁裸乳无遮挡免费网站照片| 久久精品国产清高在天天线| 在线天堂最新版资源| 国产成人a区在线观看| 宅男免费午夜| 亚洲va日本ⅴa欧美va伊人久久| 精品99又大又爽又粗少妇毛片 | 最近最新中文字幕大全免费视频| 美女免费视频网站| 亚洲精品一卡2卡三卡4卡5卡| 国产精品电影一区二区三区| 一级作爱视频免费观看| 色综合站精品国产| 18禁黄网站禁片免费观看直播| 一区二区三区激情视频| 国产精品野战在线观看| 人妻久久中文字幕网| 日韩欧美精品免费久久 | 97碰自拍视频| 韩国av一区二区三区四区| 色综合欧美亚洲国产小说| 丝袜美腿在线中文| 久久久久久国产a免费观看| 亚洲美女黄片视频| 91九色精品人成在线观看| 黑人欧美特级aaaaaa片| 男人和女人高潮做爰伦理| 在线看三级毛片| 欧美日韩中文字幕国产精品一区二区三区| 18美女黄网站色大片免费观看| 色尼玛亚洲综合影院| 亚洲精品影视一区二区三区av| 欧美性猛交黑人性爽| 97超视频在线观看视频| 国产精品久久久久久亚洲av鲁大| av天堂中文字幕网| 人妻丰满熟妇av一区二区三区| 人妻久久中文字幕网| 国产探花在线观看一区二区| 国产一区二区三区在线臀色熟女| 黄色女人牲交| 少妇高潮的动态图| 日日摸夜夜添夜夜添小说| 亚洲av不卡在线观看| 亚洲五月天丁香| 国产精品,欧美在线| 久久久久亚洲av毛片大全| 成年版毛片免费区| 日韩欧美在线乱码| 欧美大码av| 亚洲avbb在线观看| 亚洲人与动物交配视频| 久久精品91无色码中文字幕| 成人特级av手机在线观看| 黄色女人牲交| 欧美av亚洲av综合av国产av| av福利片在线观看| 免费在线观看日本一区| 欧美日韩亚洲国产一区二区在线观看| 最近在线观看免费完整版| 国模一区二区三区四区视频| 舔av片在线| 久久精品影院6| 他把我摸到了高潮在线观看| 欧美bdsm另类| 2021天堂中文幕一二区在线观| 不卡一级毛片| 少妇熟女aⅴ在线视频| 脱女人内裤的视频| 一卡2卡三卡四卡精品乱码亚洲| 麻豆国产av国片精品| 免费看十八禁软件| 长腿黑丝高跟| 黄色女人牲交| 美女免费视频网站| 免费看a级黄色片| 精品不卡国产一区二区三区| 欧美日本视频| 成年女人毛片免费观看观看9| 最近最新免费中文字幕在线| 黄色视频,在线免费观看| 国产一区二区亚洲精品在线观看| 日本与韩国留学比较| 日韩欧美三级三区| 国产爱豆传媒在线观看| 精品一区二区三区av网在线观看| 少妇裸体淫交视频免费看高清| 国产伦一二天堂av在线观看| 中文字幕精品亚洲无线码一区| 亚洲av免费在线观看| 91在线观看av| 天堂√8在线中文| av视频在线观看入口| 精品久久久久久久久久免费视频| 最后的刺客免费高清国语| 琪琪午夜伦伦电影理论片6080| 欧美激情久久久久久爽电影| 又紧又爽又黄一区二区| 女警被强在线播放| 日韩欧美精品免费久久 | 国产精品国产高清国产av| 一卡2卡三卡四卡精品乱码亚洲| 1024手机看黄色片| 长腿黑丝高跟| 精品乱码久久久久久99久播| 久久精品国产自在天天线| 国产男靠女视频免费网站| 成年女人看的毛片在线观看| 成年版毛片免费区| 伊人久久精品亚洲午夜| 亚洲一区二区三区不卡视频| 少妇的逼水好多| 婷婷丁香在线五月| 精品无人区乱码1区二区| 午夜精品在线福利| 国产高清有码在线观看视频| 一个人免费在线观看的高清视频| 在线免费观看的www视频| 搡老岳熟女国产| 欧美日本亚洲视频在线播放| ponron亚洲| 亚洲av二区三区四区| h日本视频在线播放| 91九色精品人成在线观看| 看片在线看免费视频| 国产成人啪精品午夜网站| 97超级碰碰碰精品色视频在线观看| 国产精华一区二区三区| 国产三级黄色录像| 97人妻精品一区二区三区麻豆| 熟妇人妻久久中文字幕3abv| 嫁个100分男人电影在线观看| 热99在线观看视频| 99国产精品一区二区三区| 99热6这里只有精品| 99在线人妻在线中文字幕| 国产精品亚洲一级av第二区| 美女 人体艺术 gogo| 一进一出好大好爽视频| 精品福利观看| 成熟少妇高潮喷水视频| 欧美一区二区国产精品久久精品| 日本成人三级电影网站| 特大巨黑吊av在线直播| 我要搜黄色片| 亚洲精品在线观看二区| 国产亚洲精品综合一区在线观看| 99精品久久久久人妻精品| 九九久久精品国产亚洲av麻豆| 国内精品一区二区在线观看| 久久久久久久亚洲中文字幕 | 桃红色精品国产亚洲av| 亚洲五月婷婷丁香| 国产欧美日韩一区二区精品| 麻豆成人午夜福利视频| 午夜福利在线观看吧| 国模一区二区三区四区视频| 久久九九热精品免费| 亚洲最大成人手机在线| 别揉我奶头~嗯~啊~动态视频| 麻豆国产av国片精品| 国产美女午夜福利| 三级毛片av免费| 最新在线观看一区二区三区| 色综合欧美亚洲国产小说| 最近视频中文字幕2019在线8| 床上黄色一级片| 亚洲精品国产精品久久久不卡| 一个人免费在线观看的高清视频| 嫩草影院精品99| 久久久久久久午夜电影| 偷拍熟女少妇极品色| 91在线精品国自产拍蜜月| 91精品伊人久久大香线蕉| 能在线免费看毛片的网站| 午夜老司机福利剧场| 春色校园在线视频观看| 国产美女午夜福利| 日韩一区二区三区影片| 男人狂女人下面高潮的视频| 精品99又大又爽又粗少妇毛片| av在线老鸭窝| 国产免费又黄又爽又色| 国产色婷婷99| 亚洲无线观看免费| 岛国毛片在线播放| 大香蕉久久网| 岛国毛片在线播放| 免费观看av网站的网址| 国产一区有黄有色的免费视频 | 一级毛片aaaaaa免费看小| 日韩在线高清观看一区二区三区| 日韩成人伦理影院| 网址你懂的国产日韩在线| av.在线天堂| 蜜臀久久99精品久久宅男| 黄片无遮挡物在线观看| 亚洲伊人久久精品综合| 久久精品综合一区二区三区| 日日摸夜夜添夜夜爱| 最近2019中文字幕mv第一页| 别揉我奶头 嗯啊视频| 国产激情偷乱视频一区二区| 欧美激情久久久久久爽电影| 午夜福利网站1000一区二区三区| 深夜a级毛片| 最近最新中文字幕免费大全7| 色哟哟·www| 两个人视频免费观看高清| 精品欧美国产一区二区三| 免费看日本二区| 观看免费一级毛片| 伊人久久精品亚洲午夜| 91久久精品电影网| 九九久久精品国产亚洲av麻豆| 国产黄a三级三级三级人| 五月天丁香电影| 精品一区二区三区视频在线| 非洲黑人性xxxx精品又粗又长| 亚洲图色成人| av在线亚洲专区| 亚洲欧美一区二区三区黑人 | 亚洲精品成人av观看孕妇| 中文字幕制服av| 午夜精品一区二区三区免费看| 亚洲最大成人中文| 你懂的网址亚洲精品在线观看| 午夜福利在线观看免费完整高清在| 国产精品久久久久久久电影| 免费大片黄手机在线观看| 亚洲欧美日韩卡通动漫| 五月伊人婷婷丁香| 啦啦啦韩国在线观看视频| 91精品伊人久久大香线蕉| 亚洲精品影视一区二区三区av| 婷婷色综合www| 美女内射精品一级片tv| 能在线免费看毛片的网站| 亚洲欧美清纯卡通| 麻豆成人午夜福利视频| 久久精品夜色国产| 别揉我奶头 嗯啊视频| 久久久久免费精品人妻一区二区| 中文字幕av成人在线电影| 中文精品一卡2卡3卡4更新| 久久久久久国产a免费观看| 国产精品一区二区三区四区久久| 亚洲欧洲国产日韩| 2021天堂中文幕一二区在线观| 国产伦精品一区二区三区四那| 午夜精品一区二区三区免费看| 久久精品国产亚洲av天美| 久久草成人影院| 麻豆av噜噜一区二区三区| 亚洲四区av| 乱系列少妇在线播放| 熟女人妻精品中文字幕| 国产人妻一区二区三区在| 午夜精品在线福利| 免费看日本二区| 亚洲一级一片aⅴ在线观看| 五月天丁香电影| 99久久人妻综合| 亚洲av在线观看美女高潮| 高清在线视频一区二区三区| 丝袜喷水一区| 国产 一区 欧美 日韩| 日韩欧美一区视频在线观看 | 国产有黄有色有爽视频| 国产毛片a区久久久久| 全区人妻精品视频| 91午夜精品亚洲一区二区三区| 亚洲真实伦在线观看| 精华霜和精华液先用哪个| 国产精品av视频在线免费观看| 亚洲自偷自拍三级| 久久久久久久午夜电影| 激情 狠狠 欧美| 深爱激情五月婷婷| 69av精品久久久久久| 久久国内精品自在自线图片| 一本久久精品| 少妇人妻精品综合一区二区| 久热久热在线精品观看| 国产在线一区二区三区精| 日韩视频在线欧美| 欧美人与善性xxx| 精品国产露脸久久av麻豆 | 中文字幕久久专区| 亚洲欧美成人精品一区二区| 成年版毛片免费区| 国产熟女欧美一区二区| 人妻制服诱惑在线中文字幕| 亚洲性久久影院| 18+在线观看网站| 黄色一级大片看看| 国产熟女欧美一区二区| 免费观看av网站的网址| 全区人妻精品视频| 亚洲av日韩在线播放| 成年女人看的毛片在线观看| 神马国产精品三级电影在线观看| 97人妻精品一区二区三区麻豆| 日本黄大片高清| 国产永久视频网站| 最新中文字幕久久久久| 亚洲综合色惰| 中国国产av一级| 性色avwww在线观看| 亚洲天堂国产精品一区在线| 国产成人freesex在线| 精品国产露脸久久av麻豆 | 亚洲欧美成人精品一区二区| 在线a可以看的网站| freevideosex欧美| 国产成人91sexporn| 久久久久久久久久黄片| 亚洲欧洲日产国产| 午夜福利视频精品| 99热6这里只有精品| 免费少妇av软件| 不卡视频在线观看欧美| 纵有疾风起免费观看全集完整版 | 国产 一区 欧美 日韩| 亚洲性久久影院| 黄色日韩在线| 男女啪啪激烈高潮av片| 黄片wwwwww| 亚洲国产欧美在线一区| av一本久久久久| 一级毛片久久久久久久久女| 亚洲欧美清纯卡通| 亚洲精品第二区| 极品少妇高潮喷水抽搐| 国产精品一区二区在线观看99 | 欧美最新免费一区二区三区| 日韩av不卡免费在线播放| 亚洲国产精品sss在线观看| 综合色丁香网| 在现免费观看毛片| 成年av动漫网址| 高清视频免费观看一区二区 | 少妇裸体淫交视频免费看高清| 亚洲国产精品成人综合色| 成人性生交大片免费视频hd| 精品亚洲乱码少妇综合久久| 亚洲精品456在线播放app| 国产成人精品一,二区| 99视频精品全部免费 在线| 日韩av不卡免费在线播放| 高清日韩中文字幕在线| 国产高清有码在线观看视频| 麻豆成人午夜福利视频| 久久久久久久大尺度免费视频| 免费看日本二区| 少妇的逼水好多| 老司机影院毛片| 国产高潮美女av| 26uuu在线亚洲综合色| 中国国产av一级| 性色avwww在线观看| 免费黄色在线免费观看| 男人爽女人下面视频在线观看| 亚洲四区av| 国产乱来视频区| 蜜桃久久精品国产亚洲av| 亚洲av中文字字幕乱码综合| 免费av不卡在线播放| 国产爱豆传媒在线观看| 日本色播在线视频| 国产精品日韩av在线免费观看| 丰满人妻一区二区三区视频av| 最近最新中文字幕免费大全7| 精品一区在线观看国产| 男女国产视频网站| 亚洲va在线va天堂va国产| 久久99热6这里只有精品| 欧美97在线视频| 亚洲色图av天堂| 久久精品人妻少妇| 久久精品久久久久久噜噜老黄| or卡值多少钱| 好男人在线观看高清免费视频| 国产91av在线免费观看| 麻豆成人av视频| 日韩av在线免费看完整版不卡| 国产精品一及| 天美传媒精品一区二区| 久99久视频精品免费| 国产综合懂色| 色综合色国产| 国产精品福利在线免费观看| 麻豆国产97在线/欧美| 一边亲一边摸免费视频| 欧美zozozo另类| 性色avwww在线观看| 国产91av在线免费观看| 熟妇人妻久久中文字幕3abv| 天堂中文最新版在线下载 | 久久综合国产亚洲精品| 久久99热这里只频精品6学生| 欧美97在线视频| 麻豆成人午夜福利视频| 国产在视频线在精品| 国产午夜精品论理片| 日韩不卡一区二区三区视频在线| 亚洲精品国产av成人精品| 久久这里有精品视频免费| 国产乱人视频| 综合色av麻豆| 18禁动态无遮挡网站| 亚洲色图av天堂| www.色视频.com| 成人性生交大片免费视频hd| 亚洲精品乱码久久久v下载方式| 夜夜爽夜夜爽视频| 特级一级黄色大片| 亚洲成色77777| 国产黄频视频在线观看| 婷婷六月久久综合丁香| 精品久久国产蜜桃| 丰满乱子伦码专区| av在线天堂中文字幕| 日日摸夜夜添夜夜爱| 美女cb高潮喷水在线观看| 国产精品久久久久久久久免| 熟妇人妻不卡中文字幕| videos熟女内射| 小蜜桃在线观看免费完整版高清| 亚洲av中文字字幕乱码综合| 三级男女做爰猛烈吃奶摸视频| 天堂中文最新版在线下载 | 精品酒店卫生间| 99久久九九国产精品国产免费| 人妻夜夜爽99麻豆av| 亚洲av中文av极速乱| 亚洲欧美成人综合另类久久久| 有码 亚洲区| 超碰av人人做人人爽久久| 亚洲综合精品二区| 日韩欧美国产在线观看| 国产av码专区亚洲av| 一级毛片 在线播放| 国内精品宾馆在线| 国产探花在线观看一区二区| 成人毛片60女人毛片免费| 亚洲熟妇中文字幕五十中出| 亚洲四区av| 亚洲av成人精品一区久久| 国产精品一区二区三区四区免费观看| av免费观看日本| 美女xxoo啪啪120秒动态图| 男插女下体视频免费在线播放| 亚洲人成网站高清观看| 成人毛片60女人毛片免费| 日韩av免费高清视频| 最近2019中文字幕mv第一页| 91精品一卡2卡3卡4卡| 久久精品人妻少妇| 亚洲av电影在线观看一区二区三区 | 国产亚洲5aaaaa淫片| 淫秽高清视频在线观看| a级毛色黄片| 伊人久久精品亚洲午夜| 搡女人真爽免费视频火全软件| 久久97久久精品| 亚洲国产精品国产精品| 熟女电影av网| 亚洲欧洲国产日韩| 午夜福利视频精品| 真实男女啪啪啪动态图| 男的添女的下面高潮视频| 国产一区二区亚洲精品在线观看| 国产精品国产三级国产专区5o| 色5月婷婷丁香| 国产精品福利在线免费观看| 99九九线精品视频在线观看视频| 欧美成人a在线观看| 国产人妻一区二区三区在| 久久久久网色| 大香蕉97超碰在线| 国产精品美女特级片免费视频播放器| 欧美zozozo另类| 日本欧美国产在线视频| 成人综合一区亚洲| 亚洲电影在线观看av| 国产亚洲午夜精品一区二区久久 | 美女国产视频在线观看| 午夜福利视频1000在线观看| 成人午夜精彩视频在线观看| 亚洲成人精品中文字幕电影| 听说在线观看完整版免费高清| 欧美3d第一页| 99久久九九国产精品国产免费| 亚洲18禁久久av| 日本黄色片子视频| 大陆偷拍与自拍| 日韩制服骚丝袜av| 国产精品一区www在线观看| 十八禁网站网址无遮挡 | 亚洲av成人精品一区久久| 久久精品久久久久久噜噜老黄| 一本久久精品| 亚洲经典国产精华液单| 一级片'在线观看视频| 亚洲精品中文字幕在线视频 | 国产亚洲最大av| 亚洲天堂国产精品一区在线| 国产女主播在线喷水免费视频网站 | 免费大片18禁| 亚洲在线观看片| 男人舔女人下体高潮全视频| 午夜激情欧美在线| 成年女人看的毛片在线观看| 观看美女的网站| 国产综合懂色| 我的老师免费观看完整版| 99九九线精品视频在线观看视频| 亚洲av中文字字幕乱码综合| 久久这里只有精品中国| 狂野欧美白嫩少妇大欣赏| 久久久久久久久久久免费av| 国产三级在线视频| 国语对白做爰xxxⅹ性视频网站| 免费av观看视频| 久久人人爽人人爽人人片va| 干丝袜人妻中文字幕| 欧美xxⅹ黑人| 精品一区二区三区视频在线| 一个人看的www免费观看视频| 床上黄色一级片| 乱码一卡2卡4卡精品| 免费高清在线观看视频在线观看| 毛片一级片免费看久久久久| 国产伦理片在线播放av一区| 国产精品国产三级国产av玫瑰| a级毛色黄片| 国产一区二区三区综合在线观看 | 嫩草影院精品99| 日韩欧美国产在线观看| av在线观看视频网站免费| 精品久久久精品久久久| 99re6热这里在线精品视频| 欧美性感艳星| 日日摸夜夜添夜夜添av毛片| 乱人视频在线观看| 可以在线观看毛片的网站| 国产精品一区二区三区四区久久| 亚洲在线观看片| 亚洲av福利一区| 一级毛片aaaaaa免费看小| 一级片'在线观看视频| 亚洲三级黄色毛片| 国产成人精品婷婷| 成人综合一区亚洲| av免费在线看不卡| 久久久a久久爽久久v久久| av免费观看日本| 特大巨黑吊av在线直播| 深夜a级毛片| 成人亚洲精品av一区二区| 精品国产露脸久久av麻豆 | 亚洲18禁久久av| 久久国内精品自在自线图片| 国产精品一区二区三区四区免费观看| 嫩草影院新地址| 小蜜桃在线观看免费完整版高清| av播播在线观看一区| 人妻少妇偷人精品九色| 免费无遮挡裸体视频| 男女国产视频网站| 免费黄频网站在线观看国产| 日日撸夜夜添| 亚洲一区高清亚洲精品| 日本免费在线观看一区| 成人无遮挡网站| 欧美不卡视频在线免费观看|