• <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
    馬克斯·穆恩的瘋狂割草機的秘密
    亚洲av欧美aⅴ国产| 国产精品女同一区二区软件| a级毛片在线看网站| 成年人免费黄色播放视频| 精品国产一区二区三区四区第35| 美女视频免费永久观看网站| 麻豆精品久久久久久蜜桃| 黄网站色视频无遮挡免费观看| 侵犯人妻中文字幕一二三四区| 亚洲精品日本国产第一区| a级毛片黄视频| 成年人午夜在线观看视频| 男人爽女人下面视频在线观看| 夜夜骑夜夜射夜夜干| 久久久久久久久久久免费av| 色网站视频免费| 99香蕉大伊视频| 国产激情久久老熟女| 精品久久久精品久久久| 亚洲国产毛片av蜜桃av| 久久精品国产亚洲av高清一级| 精品一区二区免费观看| 亚洲国产最新在线播放| 国产野战对白在线观看| 99精国产麻豆久久婷婷| 中国三级夫妇交换| videossex国产| 成人二区视频| 18禁动态无遮挡网站| 国产深夜福利视频在线观看| 久久久久国产网址| 这个男人来自地球电影免费观看 | 亚洲人成77777在线视频| 一区在线观看完整版| 亚洲av电影在线观看一区二区三区| 午夜激情av网站| 在线观看一区二区三区激情| 制服丝袜香蕉在线| 一边亲一边摸免费视频| 在线免费观看不下载黄p国产| 男女无遮挡免费网站观看| 国产成人精品久久久久久| 免费不卡的大黄色大毛片视频在线观看| 国产免费一区二区三区四区乱码| 免费黄频网站在线观看国产| 叶爱在线成人免费视频播放| 一本—道久久a久久精品蜜桃钙片| 王馨瑶露胸无遮挡在线观看| 亚洲欧美色中文字幕在线| 国产精品av久久久久免费| 亚洲国产av影院在线观看| 在线观看免费高清a一片| 国产一级毛片在线| 亚洲经典国产精华液单| 日韩中字成人| av国产久精品久网站免费入址| 久久精品国产亚洲av天美| 国产白丝娇喘喷水9色精品| 99热国产这里只有精品6| 久久99热这里只频精品6学生| 婷婷色av中文字幕| 亚洲在久久综合| 亚洲美女黄色视频免费看| 国产精品女同一区二区软件| 欧美成人午夜免费资源| 国产乱人偷精品视频| 日日爽夜夜爽网站| 精品亚洲成国产av| 欧美亚洲日本最大视频资源| 欧美成人精品欧美一级黄| 麻豆乱淫一区二区| a级毛片黄视频| 在线免费观看不下载黄p国产| 国产精品99久久99久久久不卡 | 国产精品欧美亚洲77777| 可以免费在线观看a视频的电影网站 | 韩国精品一区二区三区| 午夜免费观看性视频| 色94色欧美一区二区| 丝袜人妻中文字幕| 18禁裸乳无遮挡动漫免费视频| 最近中文字幕高清免费大全6| 9热在线视频观看99| 免费看av在线观看网站| 亚洲精品国产av蜜桃| 成年女人在线观看亚洲视频| 深夜精品福利| 不卡视频在线观看欧美| 亚洲av日韩在线播放| 69精品国产乱码久久久| av网站在线播放免费| 啦啦啦视频在线资源免费观看| 狂野欧美激情性bbbbbb| 久久精品国产a三级三级三级| a级片在线免费高清观看视频| 久久国产精品大桥未久av| 麻豆av在线久日| av片东京热男人的天堂| 久久精品久久精品一区二区三区| 亚洲欧美日韩另类电影网站| 久久精品aⅴ一区二区三区四区 | 婷婷色综合大香蕉| 成人国产av品久久久| 老汉色∧v一级毛片| 日日啪夜夜爽| 免费在线观看完整版高清| 黄频高清免费视频| 国产精品一国产av| 各种免费的搞黄视频| 久久97久久精品| 国产精品免费视频内射| 久久精品国产鲁丝片午夜精品| 亚洲欧美成人精品一区二区| 男的添女的下面高潮视频| 一级毛片 在线播放| 久久 成人 亚洲| 美国免费a级毛片| 国产精品麻豆人妻色哟哟久久| 人成视频在线观看免费观看| 黄色怎么调成土黄色| 男人舔女人的私密视频| 国产av精品麻豆| 国产男人的电影天堂91| 久热久热在线精品观看| 久久 成人 亚洲| 亚洲精品一区蜜桃| 日韩av在线免费看完整版不卡| 国产综合精华液| 男人舔女人的私密视频| 电影成人av| av天堂久久9| 国产一级毛片在线| 中文字幕另类日韩欧美亚洲嫩草| 国产精品一区二区在线不卡| 涩涩av久久男人的天堂| 天天操日日干夜夜撸| 国产亚洲最大av| 蜜桃国产av成人99| 久久精品国产亚洲av高清一级| 国产片内射在线| 飞空精品影院首页| 寂寞人妻少妇视频99o| 国产人伦9x9x在线观看 | 成年女人毛片免费观看观看9 | 综合色丁香网| 亚洲第一av免费看| 午夜福利网站1000一区二区三区| 日韩成人av中文字幕在线观看| www.自偷自拍.com| 婷婷色麻豆天堂久久| 啦啦啦在线免费观看视频4| 女人高潮潮喷娇喘18禁视频| 2022亚洲国产成人精品| 嫩草影院入口| 欧美97在线视频| 亚洲精品日韩在线中文字幕| 好男人视频免费观看在线| 国产色婷婷99| 26uuu在线亚洲综合色| 一本大道久久a久久精品| 久久这里有精品视频免费| 大香蕉久久成人网| 高清欧美精品videossex| 久热这里只有精品99| 91aial.com中文字幕在线观看| 亚洲成人手机| 国产又色又爽无遮挡免| 亚洲国产av影院在线观看| 国产精品久久久久久精品电影小说| 热99久久久久精品小说推荐| 桃花免费在线播放| 侵犯人妻中文字幕一二三四区| 国产成人精品福利久久| 啦啦啦在线观看免费高清www| 97人妻天天添夜夜摸| 亚洲久久久国产精品| 菩萨蛮人人尽说江南好唐韦庄| 9热在线视频观看99| 久久久久精品人妻al黑| www.精华液| 香蕉精品网在线| 国产av码专区亚洲av| 欧美国产精品一级二级三级| 99久久人妻综合| 日韩av不卡免费在线播放| 亚洲精品av麻豆狂野| 欧美在线黄色| 搡女人真爽免费视频火全软件| 在现免费观看毛片| 国产精品一区二区在线观看99| 大片电影免费在线观看免费| 国产深夜福利视频在线观看| 中文欧美无线码| 精品国产国语对白av| 精品人妻在线不人妻| 最新中文字幕久久久久| 韩国av在线不卡| 亚洲av电影在线进入| 一级毛片电影观看| 国产精品麻豆人妻色哟哟久久| 久久精品亚洲av国产电影网| a级毛片黄视频| 日本黄色日本黄色录像| 一级毛片电影观看| 国产精品麻豆人妻色哟哟久久| 久久99蜜桃精品久久| 亚洲欧美一区二区三区久久| 午夜福利,免费看| 成人国语在线视频| 18禁观看日本| 亚洲人成电影观看| 老司机影院成人| 捣出白浆h1v1| 91久久精品国产一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 久久99蜜桃精品久久| 一区二区三区乱码不卡18| 国产精品不卡视频一区二区| 免费久久久久久久精品成人欧美视频| 乱人伦中国视频| 日韩av在线免费看完整版不卡| 嫩草影院入口| 欧美+日韩+精品| 亚洲av福利一区| 黄色 视频免费看| 纵有疾风起免费观看全集完整版| 精品久久久精品久久久| 欧美激情极品国产一区二区三区| av一本久久久久| 亚洲精品一区蜜桃| 看非洲黑人一级黄片| 国产精品久久久久久精品电影小说| 大香蕉久久网| 观看美女的网站| 亚洲国产看品久久| 七月丁香在线播放| 九九爱精品视频在线观看| 丝袜人妻中文字幕| 女性被躁到高潮视频| 久久av网站| av免费在线看不卡| 美女视频免费永久观看网站| 日韩一本色道免费dvd| 久久国产亚洲av麻豆专区| 免费黄频网站在线观看国产| 欧美bdsm另类| 亚洲五月色婷婷综合| 国产欧美亚洲国产| 国产人伦9x9x在线观看 | xxx大片免费视频| 99热全是精品| 欧美xxⅹ黑人| 一二三四中文在线观看免费高清| 国产精品久久久av美女十八| 少妇被粗大猛烈的视频| 免费黄色在线免费观看| 丝袜美腿诱惑在线| 午夜免费鲁丝| 国产午夜精品一二区理论片| 中国三级夫妇交换| 欧美成人午夜免费资源| 久久精品亚洲av国产电影网| 黄色视频在线播放观看不卡| 国产有黄有色有爽视频| 国产成人免费无遮挡视频| 秋霞伦理黄片| 日韩大片免费观看网站| 国产欧美日韩一区二区三区在线| 亚洲精品日韩在线中文字幕| 午夜老司机福利剧场| 赤兔流量卡办理| 日本wwww免费看| 一边摸一边做爽爽视频免费| 成年美女黄网站色视频大全免费| 9191精品国产免费久久| 日本黄色日本黄色录像| 中文字幕最新亚洲高清| 亚洲国产毛片av蜜桃av| 精品一区二区免费观看| 国产欧美日韩一区二区三区在线| 久久久久久久久久久免费av| 一本色道久久久久久精品综合| 免费看av在线观看网站| a级毛片在线看网站| 久久精品国产综合久久久| 伦精品一区二区三区| 久久久久国产一级毛片高清牌| 欧美亚洲日本最大视频资源| 亚洲一区二区三区欧美精品| 亚洲内射少妇av| 国产精品99久久99久久久不卡 | 免费在线观看完整版高清| 亚洲伊人久久精品综合| 欧美日韩亚洲高清精品| 18在线观看网站| 我要看黄色一级片免费的| 看免费成人av毛片| 午夜免费鲁丝| 国产毛片在线视频| 国产精品 欧美亚洲| 亚洲国产精品国产精品| 69精品国产乱码久久久| 激情视频va一区二区三区| 亚洲av欧美aⅴ国产| 日韩欧美精品免费久久| 伦理电影大哥的女人| 久久ye,这里只有精品| 国产成人精品一,二区| 精品一区二区三区四区五区乱码 | av国产精品久久久久影院| 午夜免费男女啪啪视频观看| 国产亚洲午夜精品一区二区久久| 国产免费又黄又爽又色| 中文精品一卡2卡3卡4更新| 极品少妇高潮喷水抽搐| 国产一区二区在线观看av| 晚上一个人看的免费电影| 91国产中文字幕| 在线免费观看不下载黄p国产| 午夜福利网站1000一区二区三区| 亚洲一区中文字幕在线| 精品福利永久在线观看| 欧美 亚洲 国产 日韩一| 大陆偷拍与自拍| 亚洲,一卡二卡三卡| 熟女av电影| 男女边吃奶边做爰视频| 国产免费又黄又爽又色| 一边亲一边摸免费视频| 日韩欧美一区视频在线观看| 国产精品免费大片| 午夜免费观看性视频| 国产成人精品婷婷| 国产又爽黄色视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产精品国产精品| 女人精品久久久久毛片| 国产深夜福利视频在线观看| 性少妇av在线| 日本91视频免费播放| 熟女少妇亚洲综合色aaa.| 国产在线一区二区三区精| 又大又黄又爽视频免费| 99久久人妻综合| 日韩不卡一区二区三区视频在线| 国产爽快片一区二区三区| 99香蕉大伊视频| 美女视频免费永久观看网站| 欧美bdsm另类| 老熟女久久久| 国产97色在线日韩免费| 国产日韩欧美在线精品| 大片电影免费在线观看免费| 人妻 亚洲 视频| 一级片'在线观看视频| 国产精品久久久av美女十八| 欧美精品国产亚洲| 一区二区日韩欧美中文字幕| 久久ye,这里只有精品| 久久久欧美国产精品| 最近最新中文字幕免费大全7| 在线观看www视频免费| 免费大片黄手机在线观看| 午夜福利,免费看| 日本wwww免费看| 国产成人aa在线观看| 国产在线免费精品| 久久久精品94久久精品| 晚上一个人看的免费电影| 边亲边吃奶的免费视频| 亚洲国产精品一区三区| www.自偷自拍.com| 男女下面插进去视频免费观看| 国产视频首页在线观看| 亚洲精品国产av蜜桃| 久久99一区二区三区| 亚洲成人手机| 在线天堂最新版资源| 午夜福利影视在线免费观看| 男女免费视频国产| 国产成人aa在线观看| 亚洲激情五月婷婷啪啪| 妹子高潮喷水视频| 免费播放大片免费观看视频在线观看| 99久久人妻综合| 日日撸夜夜添| 日本欧美视频一区| 激情五月婷婷亚洲| 99re6热这里在线精品视频| 日本爱情动作片www.在线观看| 亚洲成av片中文字幕在线观看 | 亚洲国产欧美日韩在线播放| 精品国产一区二区久久| 国产精品麻豆人妻色哟哟久久| 亚洲国产欧美在线一区| 亚洲精品久久成人aⅴ小说| www.精华液| 另类亚洲欧美激情| 9色porny在线观看| 日本黄色日本黄色录像| 亚洲婷婷狠狠爱综合网| 亚洲第一区二区三区不卡| 欧美国产精品va在线观看不卡| 欧美日韩一区二区视频在线观看视频在线| 在线 av 中文字幕| 汤姆久久久久久久影院中文字幕| 天堂俺去俺来也www色官网| 亚洲色图 男人天堂 中文字幕| 欧美成人精品欧美一级黄| 免费播放大片免费观看视频在线观看| 国产亚洲av片在线观看秒播厂| 一本—道久久a久久精品蜜桃钙片| 国产又爽黄色视频| 欧美精品av麻豆av| 日本黄色日本黄色录像| 日韩欧美精品免费久久| 欧美精品一区二区大全| 久久精品久久久久久久性| 观看av在线不卡| www.自偷自拍.com| 精品午夜福利在线看| 午夜福利乱码中文字幕| 亚洲av电影在线观看一区二区三区| 久久精品国产亚洲av涩爱| 卡戴珊不雅视频在线播放| 亚洲成色77777| av电影中文网址| 国产成人精品久久久久久| 777久久人妻少妇嫩草av网站| 午夜日本视频在线| 如何舔出高潮| 日韩大片免费观看网站| 成人午夜精彩视频在线观看| 国产精品女同一区二区软件| 成人亚洲精品一区在线观看| 老司机亚洲免费影院| 中文天堂在线官网| 婷婷色av中文字幕| 色婷婷av一区二区三区视频| 日本欧美视频一区| 久久影院123| 精品人妻熟女毛片av久久网站| 日日摸夜夜添夜夜爱| 国产男女超爽视频在线观看| 国产精品久久久久久精品电影小说| 日韩人妻精品一区2区三区| 亚洲三区欧美一区| 精品国产一区二区久久| 美女中出高潮动态图| www.熟女人妻精品国产| 黄色配什么色好看| 欧美 亚洲 国产 日韩一| 日韩av不卡免费在线播放| 超碰97精品在线观看| 亚洲欧美成人综合另类久久久| 新久久久久国产一级毛片| 亚洲内射少妇av| 午夜福利在线免费观看网站| 精品少妇内射三级| 边亲边吃奶的免费视频| 亚洲精品中文字幕在线视频| 亚洲av男天堂| 人妻 亚洲 视频| 日本av手机在线免费观看| 色94色欧美一区二区| 在现免费观看毛片| 91国产中文字幕| 日韩制服丝袜自拍偷拍| 大片电影免费在线观看免费| 伦精品一区二区三区| 国产免费一区二区三区四区乱码| 国产有黄有色有爽视频| 久久久久国产精品人妻一区二区| 青草久久国产| 国产日韩欧美亚洲二区| 男女边吃奶边做爰视频| 黄色配什么色好看| 国产精品嫩草影院av在线观看| 国产极品天堂在线| 九九爱精品视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 宅男免费午夜| 国产免费一区二区三区四区乱码| 电影成人av| 如何舔出高潮| 精品午夜福利在线看| 国产伦理片在线播放av一区| videossex国产| 国产女主播在线喷水免费视频网站| 少妇熟女欧美另类| 人人妻人人澡人人看| 一级爰片在线观看| 9色porny在线观看| 亚洲成人av在线免费| 国产精品久久久久久精品电影小说| 亚洲美女黄色视频免费看| 男女午夜视频在线观看| 国产人伦9x9x在线观看 | 曰老女人黄片| 欧美日韩精品成人综合77777| 少妇精品久久久久久久| 久久热在线av| 国产av码专区亚洲av| 亚洲五月色婷婷综合| 亚洲欧美成人综合另类久久久| 青春草视频在线免费观看| 综合色丁香网| 日本-黄色视频高清免费观看| 亚洲精品成人av观看孕妇| 亚洲欧美成人综合另类久久久| 男女无遮挡免费网站观看| 日本av手机在线免费观看| 肉色欧美久久久久久久蜜桃| 久久久精品区二区三区| 亚洲一级一片aⅴ在线观看| 日韩精品免费视频一区二区三区| 亚洲色图综合在线观看| 男女下面插进去视频免费观看| 999精品在线视频| 成人国产av品久久久| 免费观看在线日韩| 热re99久久精品国产66热6| 十分钟在线观看高清视频www| 国产成人午夜福利电影在线观看| 午夜日韩欧美国产| 这个男人来自地球电影免费观看 | 精品国产一区二区三区久久久樱花| 777久久人妻少妇嫩草av网站| 日韩伦理黄色片| 成年动漫av网址| 日日摸夜夜添夜夜爱| 交换朋友夫妻互换小说| 天天躁夜夜躁狠狠躁躁| 国产精品欧美亚洲77777| 国产白丝娇喘喷水9色精品| 欧美精品亚洲一区二区| 精品一区二区三区四区五区乱码 | 国产乱来视频区| 精品久久久久久电影网| 亚洲视频免费观看视频| 丝袜美足系列| 国产黄色免费在线视频| 涩涩av久久男人的天堂| 久久综合国产亚洲精品| 欧美日韩一区二区视频在线观看视频在线| 啦啦啦啦在线视频资源| 亚洲精品乱久久久久久| 久久久久精品人妻al黑| 欧美日韩av久久| 看非洲黑人一级黄片| 日韩一区二区视频免费看| 热re99久久国产66热| 亚洲精品日韩在线中文字幕| 久久精品国产a三级三级三级| 成人影院久久| 免费观看av网站的网址| 国产精品蜜桃在线观看| 多毛熟女@视频| 国产精品一区二区在线观看99| 在线免费观看不下载黄p国产| 久久久久久伊人网av| 人人妻人人爽人人添夜夜欢视频| 日本欧美国产在线视频| 日韩伦理黄色片| 赤兔流量卡办理| 性色avwww在线观看| 国产一区二区三区av在线| 国产免费又黄又爽又色| 卡戴珊不雅视频在线播放| 亚洲一区二区三区欧美精品| 国产在线一区二区三区精| 精品第一国产精品| 欧美xxⅹ黑人| 精品99又大又爽又粗少妇毛片| 丰满乱子伦码专区| 欧美日韩综合久久久久久| 麻豆乱淫一区二区| 日韩中文字幕视频在线看片| 国产一区二区三区综合在线观看| 亚洲欧美中文字幕日韩二区| 免费播放大片免费观看视频在线观看| 蜜桃在线观看..| 夜夜骑夜夜射夜夜干| 欧美日韩精品成人综合77777| 80岁老熟妇乱子伦牲交| 午夜福利影视在线免费观看| 观看av在线不卡| 赤兔流量卡办理| 夜夜骑夜夜射夜夜干| 中文字幕人妻熟女乱码| 午夜福利乱码中文字幕| 日韩欧美精品免费久久| 国产熟女午夜一区二区三区| 赤兔流量卡办理| 春色校园在线视频观看| 亚洲色图 男人天堂 中文字幕| h视频一区二区三区| 精品久久久精品久久久| 69精品国产乱码久久久| 国产精品国产av在线观看| 国产精品久久久久成人av| 一区在线观看完整版| 少妇被粗大的猛进出69影院| 亚洲欧美一区二区三区久久| 成人黄色视频免费在线看| 人妻 亚洲 视频| 又黄又粗又硬又大视频| 久久精品久久久久久噜噜老黄| 日韩在线高清观看一区二区三区| 国产精品熟女久久久久浪| 一区二区av电影网| 亚洲久久久国产精品| 一级a爱视频在线免费观看| 七月丁香在线播放| 黄色一级大片看看|