• <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
    馬克斯·穆恩的瘋狂割草機的秘密
    伦精品一区二区三区| 国内少妇人妻偷人精品xxx网站| 99精国产麻豆久久婷婷| 日韩一本色道免费dvd| 国产精品福利在线免费观看| 色吧在线观看| 少妇猛男粗大的猛烈进出视频| 久久精品久久久久久久性| 成人午夜精彩视频在线观看| 久久女婷五月综合色啪小说| 最新中文字幕久久久久| 一本色道久久久久久精品综合| 色哟哟·www| 校园人妻丝袜中文字幕| 18禁动态无遮挡网站| 人人妻人人爽人人添夜夜欢视频 | 嫩草影院新地址| 亚洲va在线va天堂va国产| 中文字幕免费在线视频6| 不卡视频在线观看欧美| 最新中文字幕久久久久| 91精品国产国语对白视频| 男人狂女人下面高潮的视频| 99国产精品免费福利视频| 亚洲精品久久久久久婷婷小说| 一级a做视频免费观看| 免费观看的影片在线观看| 国产日韩欧美亚洲二区| 亚洲av二区三区四区| 观看av在线不卡| 综合色丁香网| 精品99又大又爽又粗少妇毛片| 国产成人freesex在线| 亚洲欧美精品专区久久| 亚洲,一卡二卡三卡| 成人二区视频| 国产白丝娇喘喷水9色精品| 最黄视频免费看| 久久青草综合色| 国产精品久久久久久精品古装| kizo精华| 日本与韩国留学比较| 成人影院久久| 亚洲精品456在线播放app| 久久av网站| 免费人妻精品一区二区三区视频| 久久久久久久精品精品| 老熟女久久久| 国产成人精品久久久久久| 亚洲精品色激情综合| 欧美 亚洲 国产 日韩一| 日本91视频免费播放| 99视频精品全部免费 在线| 亚洲国产精品专区欧美| 日韩免费高清中文字幕av| 97在线人人人人妻| 天堂中文最新版在线下载| 亚州av有码| 亚洲av免费高清在线观看| 久久久久久久亚洲中文字幕| 亚洲激情五月婷婷啪啪| 黄色怎么调成土黄色| 国产精品伦人一区二区| 国产白丝娇喘喷水9色精品| 国产免费又黄又爽又色| 老司机影院毛片| 菩萨蛮人人尽说江南好唐韦庄| 美女内射精品一级片tv| 精品亚洲成国产av| 两个人的视频大全免费| 在线观看国产h片| 久久精品久久久久久久性| 在线精品无人区一区二区三| 国产精品福利在线免费观看| 99久久综合免费| 日日啪夜夜爽| 欧美精品一区二区大全| 国产精品一区二区在线不卡| 如何舔出高潮| 国产熟女欧美一区二区| 亚洲一级一片aⅴ在线观看| 婷婷色综合www| 黄色欧美视频在线观看| 亚洲精品日韩在线中文字幕| 久久ye,这里只有精品| 日韩一区二区三区影片| 成人亚洲精品一区在线观看| 国产一区二区在线观看av| 日韩欧美精品免费久久| 99久久精品一区二区三区| a级毛片在线看网站| 丰满乱子伦码专区| 中文字幕制服av| 亚洲av二区三区四区| 十分钟在线观看高清视频www | 国产在线免费精品| 国产欧美另类精品又又久久亚洲欧美| 国产美女午夜福利| 日本黄色日本黄色录像| 97超视频在线观看视频| 99热全是精品| 亚洲成色77777| 我的女老师完整版在线观看| 少妇人妻精品综合一区二区| 国产成人freesex在线| 亚洲欧洲国产日韩| av女优亚洲男人天堂| 日韩欧美一区视频在线观看 | 国产91av在线免费观看| 天堂中文最新版在线下载| 国产精品三级大全| 特大巨黑吊av在线直播| 欧美日韩综合久久久久久| 一本色道久久久久久精品综合| 精品一区二区三区视频在线| 日本黄色日本黄色录像| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美日韩精品成人综合77777| 美女福利国产在线| 亚洲丝袜综合中文字幕| av福利片在线| av女优亚洲男人天堂| 亚洲欧美成人精品一区二区| 精品一区二区免费观看| 久久人人爽人人片av| 亚洲精华国产精华液的使用体验| 性高湖久久久久久久久免费观看| 在线观看www视频免费| 国产无遮挡羞羞视频在线观看| 亚洲精品国产av蜜桃| 国产亚洲av片在线观看秒播厂| 欧美高清成人免费视频www| 精品酒店卫生间| 亚洲精品日韩在线中文字幕| av又黄又爽大尺度在线免费看| 观看美女的网站| 人妻少妇偷人精品九色| 中文在线观看免费www的网站| 黄色一级大片看看| 国产极品天堂在线| 日韩中文字幕视频在线看片| 少妇猛男粗大的猛烈进出视频| 建设人人有责人人尽责人人享有的| a级一级毛片免费在线观看| 夜夜看夜夜爽夜夜摸| 少妇人妻 视频| 国产综合精华液| 国产91av在线免费观看| 日韩av不卡免费在线播放| 人人妻人人添人人爽欧美一区卜| 一级爰片在线观看| 3wmmmm亚洲av在线观看| 中文乱码字字幕精品一区二区三区| 日韩在线高清观看一区二区三区| 一区二区av电影网| 麻豆乱淫一区二区| 久久久久久久久久久免费av| 国产精品一区二区在线不卡| 自线自在国产av| 制服丝袜香蕉在线| av在线播放精品| 极品少妇高潮喷水抽搐| 久久热精品热| videossex国产| 日韩一区二区视频免费看| 街头女战士在线观看网站| 男男h啪啪无遮挡| 日韩视频在线欧美| 26uuu在线亚洲综合色| 美女cb高潮喷水在线观看| 国产伦精品一区二区三区四那| 国产精品无大码| 精品国产国语对白av| 日韩精品免费视频一区二区三区 | 99热这里只有是精品50| 欧美日韩视频精品一区| 久久久久网色| 国产精品.久久久| 大码成人一级视频| 建设人人有责人人尽责人人享有的| 人人妻人人澡人人爽人人夜夜| 欧美精品国产亚洲| 丝瓜视频免费看黄片| 国产免费一级a男人的天堂| 欧美性感艳星| 精品久久久噜噜| 久久99蜜桃精品久久| 欧美日韩国产mv在线观看视频| av视频免费观看在线观看| 少妇丰满av| 国产精品久久久久久精品电影小说| 99re6热这里在线精品视频| 日日撸夜夜添| 特大巨黑吊av在线直播| 婷婷色麻豆天堂久久| 热re99久久国产66热| h视频一区二区三区| 哪个播放器可以免费观看大片| 大香蕉97超碰在线| 国产精品.久久久| 亚洲欧洲日产国产| 涩涩av久久男人的天堂| 国产 一区精品| 高清欧美精品videossex| 日韩成人伦理影院| 国产亚洲欧美精品永久| 纯流量卡能插随身wifi吗| 亚洲av二区三区四区| 黑人猛操日本美女一级片| 内地一区二区视频在线| 国产精品偷伦视频观看了| 丰满少妇做爰视频| 欧美激情极品国产一区二区三区 | 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美一区二区三区国产| 熟女电影av网| 国产成人a∨麻豆精品| 80岁老熟妇乱子伦牲交| 爱豆传媒免费全集在线观看| 18禁动态无遮挡网站| 亚洲av欧美aⅴ国产| 黄色毛片三级朝国网站 | 久久97久久精品| 成人美女网站在线观看视频| 精品一区二区免费观看| 久久久久久久大尺度免费视频| 久久99一区二区三区| 91在线精品国自产拍蜜月| 大又大粗又爽又黄少妇毛片口| 国产 一区精品| 中文乱码字字幕精品一区二区三区| 欧美激情极品国产一区二区三区 | 各种免费的搞黄视频| 妹子高潮喷水视频| 1024视频免费在线观看| 色老头精品视频在线观看| 成人免费观看视频高清| 淫妇啪啪啪对白视频 | 99热网站在线观看| 三上悠亚av全集在线观看| 国产精品 国内视频| 色播在线永久视频| www.av在线官网国产| 免费av中文字幕在线| 午夜免费成人在线视频| 别揉我奶头~嗯~啊~动态视频 | 精品国产乱码久久久久久小说| 水蜜桃什么品种好| 天堂中文最新版在线下载| 日本av手机在线免费观看| 国产淫语在线视频| 最新在线观看一区二区三区| 欧美+亚洲+日韩+国产| 亚洲精品久久久久久婷婷小说| 我要看黄色一级片免费的| 日本黄色日本黄色录像| 亚洲情色 制服丝袜| 国产真人三级小视频在线观看| 日韩有码中文字幕| 亚洲精品中文字幕在线视频| 亚洲精品自拍成人| 水蜜桃什么品种好| 69av精品久久久久久 | 91成人精品电影| 国产av又大| 午夜激情av网站| 久久中文字幕一级| 亚洲欧美一区二区三区久久| 国产精品.久久久| 中国国产av一级| 男人添女人高潮全过程视频| 老司机午夜福利在线观看视频 | 日韩 欧美 亚洲 中文字幕| 日韩有码中文字幕| 少妇粗大呻吟视频| 制服诱惑二区| 国产亚洲一区二区精品| 又大又爽又粗| 国产精品免费视频内射| cao死你这个sao货| 不卡av一区二区三区| 五月开心婷婷网| 91大片在线观看| 一区福利在线观看| 高清av免费在线| 91成人精品电影| bbb黄色大片| 日韩中文字幕欧美一区二区| 久久精品成人免费网站| 水蜜桃什么品种好| 日本a在线网址| 亚洲欧美一区二区三区久久| 老司机影院毛片| 国产av国产精品国产| 欧美午夜高清在线| 正在播放国产对白刺激| 两个人免费观看高清视频| 亚洲成国产人片在线观看| 久久人人97超碰香蕉20202| 青草久久国产| 韩国高清视频一区二区三区| 桃红色精品国产亚洲av| 亚洲精品第二区| 大型av网站在线播放| 91九色精品人成在线观看| 美女高潮到喷水免费观看| 两个人看的免费小视频| 黄片小视频在线播放| 黑人猛操日本美女一级片| 久久国产精品男人的天堂亚洲| 99久久精品国产亚洲精品| 国产高清视频在线播放一区 | 精品亚洲成a人片在线观看| 国产欧美日韩综合在线一区二区| 免费在线观看影片大全网站| 69精品国产乱码久久久| 97精品久久久久久久久久精品| 亚洲一码二码三码区别大吗| 欧美性长视频在线观看| 精品免费久久久久久久清纯 | 窝窝影院91人妻| 欧美黑人欧美精品刺激| 黄色怎么调成土黄色| 欧美激情 高清一区二区三区| 国产亚洲午夜精品一区二区久久| 午夜福利乱码中文字幕| svipshipincom国产片| 美女中出高潮动态图| 亚洲精品乱久久久久久| 国产成人系列免费观看| 亚洲国产欧美在线一区| 亚洲国产看品久久| www.精华液| 青春草视频在线免费观看| 各种免费的搞黄视频| 女人高潮潮喷娇喘18禁视频| 免费不卡黄色视频| 亚洲成人免费av在线播放| 菩萨蛮人人尽说江南好唐韦庄| 91字幕亚洲| 日韩有码中文字幕| 亚洲av成人不卡在线观看播放网 | 搡老岳熟女国产| 18禁裸乳无遮挡动漫免费视频| 亚洲精品国产av蜜桃| 久9热在线精品视频| 超碰97精品在线观看| 这个男人来自地球电影免费观看| 建设人人有责人人尽责人人享有的| 这个男人来自地球电影免费观看| 国产精品一二三区在线看| 一级,二级,三级黄色视频| 日韩 亚洲 欧美在线| 高清黄色对白视频在线免费看| 免费日韩欧美在线观看| 69av精品久久久久久 | av视频免费观看在线观看| 午夜成年电影在线免费观看| kizo精华| 黄色怎么调成土黄色| 亚洲成av片中文字幕在线观看| 黑人欧美特级aaaaaa片| 欧美黄色片欧美黄色片| 精品国产乱子伦一区二区三区 | 无限看片的www在线观看| 久久国产精品大桥未久av| 亚洲中文日韩欧美视频| 国产av精品麻豆| 午夜91福利影院| 国产精品自产拍在线观看55亚洲 | 日日夜夜操网爽| 黄色视频不卡| 精品卡一卡二卡四卡免费| 久久女婷五月综合色啪小说| 制服诱惑二区| 欧美亚洲日本最大视频资源| av福利片在线| 精品人妻1区二区| 桃花免费在线播放| 国产又爽黄色视频| 欧美另类一区| 亚洲男人天堂网一区| 欧美久久黑人一区二区| 黄色视频不卡| 搡老熟女国产l中国老女人| 国产精品av久久久久免费| 国产在视频线精品| 欧美精品一区二区大全| 国产精品一区二区精品视频观看| 在线观看免费日韩欧美大片| 久久久久久人人人人人| 亚洲伊人久久精品综合| 在线观看一区二区三区激情| 国产成人啪精品午夜网站| 欧美变态另类bdsm刘玥| 欧美少妇被猛烈插入视频| 国产又色又爽无遮挡免| 久久国产精品人妻蜜桃| 国产在线视频一区二区| 国产免费一区二区三区四区乱码| 男女之事视频高清在线观看| 国产精品99久久99久久久不卡| av超薄肉色丝袜交足视频| 侵犯人妻中文字幕一二三四区| 黄片大片在线免费观看| 狠狠精品人妻久久久久久综合| 精品国产一区二区久久| 天天影视国产精品| 国产高清videossex| 亚洲精品一二三| 黄片播放在线免费| 成人免费观看视频高清| 老司机深夜福利视频在线观看 | 亚洲国产欧美日韩在线播放| 亚洲国产av新网站| 久久国产精品男人的天堂亚洲| 欧美中文综合在线视频| 男女床上黄色一级片免费看| 色播在线永久视频| 高清黄色对白视频在线免费看| 亚洲精品一卡2卡三卡4卡5卡 | 97人妻天天添夜夜摸| 中亚洲国语对白在线视频| 日韩中文字幕欧美一区二区| bbb黄色大片| videos熟女内射| 老司机午夜福利在线观看视频 | 女人久久www免费人成看片| 99久久综合免费| 少妇裸体淫交视频免费看高清 | 久久99一区二区三区| 一二三四在线观看免费中文在| 久久精品亚洲av国产电影网| 自线自在国产av| 一本—道久久a久久精品蜜桃钙片| 久久狼人影院| 久久国产精品人妻蜜桃| 久久久国产精品麻豆| 美国免费a级毛片| 亚洲国产精品一区二区三区在线| 亚洲精品第二区| 久9热在线精品视频| 午夜免费成人在线视频| 美女扒开内裤让男人捅视频| 亚洲av美国av| netflix在线观看网站| 最近最新中文字幕大全免费视频| 久久狼人影院| 欧美大码av| 亚洲精品国产精品久久久不卡| 欧美日韩国产mv在线观看视频| 男女免费视频国产| 国产亚洲av高清不卡| 一区在线观看完整版| 女人高潮潮喷娇喘18禁视频| 夫妻午夜视频| 亚洲精品美女久久av网站| netflix在线观看网站| 国产成人精品久久二区二区91| 欧美成人午夜精品| 亚洲欧美成人综合另类久久久| av在线app专区| 欧美中文综合在线视频| 日本黄色日本黄色录像| 中文字幕精品免费在线观看视频| 最近最新免费中文字幕在线| 老熟妇仑乱视频hdxx| 久久这里只有精品19| 精品一区二区三卡| 国产成人av激情在线播放| 国产男女超爽视频在线观看| 欧美精品亚洲一区二区| 一区在线观看完整版| 丁香六月天网| 国产精品一区二区在线不卡| 色94色欧美一区二区| 女性被躁到高潮视频| 免费高清在线观看视频在线观看| 亚洲av成人不卡在线观看播放网 | 国产欧美日韩精品亚洲av| 在线观看舔阴道视频| 亚洲黑人精品在线| 欧美 日韩 精品 国产| 首页视频小说图片口味搜索| 亚洲成人国产一区在线观看| netflix在线观看网站| 美女大奶头黄色视频| 999精品在线视频| 成人国产av品久久久| 欧美日韩中文字幕国产精品一区二区三区 | 老司机影院毛片| 亚洲午夜精品一区,二区,三区| 欧美成狂野欧美在线观看| 精品人妻一区二区三区麻豆| 精品一区在线观看国产| 色94色欧美一区二区| 国产精品久久久久久精品电影小说| 日本av免费视频播放| 国产日韩欧美亚洲二区| 搡老乐熟女国产| av有码第一页| 黄色视频在线播放观看不卡| av福利片在线| avwww免费| 中文字幕最新亚洲高清| 色94色欧美一区二区| 丝袜脚勾引网站| 别揉我奶头~嗯~啊~动态视频 | 青青草视频在线视频观看| 99精国产麻豆久久婷婷| 亚洲精品成人av观看孕妇| 国产精品秋霞免费鲁丝片| 日韩精品免费视频一区二区三区| 精品国内亚洲2022精品成人 | 午夜免费观看性视频| 日日夜夜操网爽| 国产福利在线免费观看视频| 国产精品久久久久久精品古装| 国产精品亚洲av一区麻豆| 80岁老熟妇乱子伦牲交| 大片免费播放器 马上看| 最近最新免费中文字幕在线| 国产深夜福利视频在线观看| 国产片内射在线| 一个人免费在线观看的高清视频 | 久久人妻福利社区极品人妻图片| 黄网站色视频无遮挡免费观看| 免费不卡黄色视频| 黄片小视频在线播放| 一本久久精品| 超色免费av| 国产野战对白在线观看| 人人澡人人妻人| 中文字幕人妻丝袜制服| 亚洲精品国产区一区二| av网站在线播放免费| 19禁男女啪啪无遮挡网站| 亚洲国产精品成人久久小说| 女警被强在线播放| 少妇被粗大的猛进出69影院| 黑人操中国人逼视频| 亚洲精品日韩在线中文字幕| 97在线人人人人妻| 成人18禁高潮啪啪吃奶动态图| 国内毛片毛片毛片毛片毛片| 亚洲精品在线美女| 大码成人一级视频| 黄色视频不卡| 一二三四在线观看免费中文在| 国产麻豆69| 亚洲成人国产一区在线观看| 高清欧美精品videossex| 后天国语完整版免费观看| 美女午夜性视频免费| 窝窝影院91人妻| 日韩人妻精品一区2区三区| 嫩草影视91久久| 成人影院久久| 免费女性裸体啪啪无遮挡网站| 婷婷色av中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 中文字幕精品免费在线观看视频| 男人操女人黄网站| av欧美777| 午夜日韩欧美国产| 国产黄频视频在线观看| 久久国产精品男人的天堂亚洲| 亚洲精品av麻豆狂野| 久久免费观看电影| 国产亚洲精品第一综合不卡| 50天的宝宝边吃奶边哭怎么回事| 纵有疾风起免费观看全集完整版| 天天躁夜夜躁狠狠躁躁| 伦理电影免费视频| 国产精品av久久久久免费| 婷婷色av中文字幕| 国产激情久久老熟女| 午夜福利在线免费观看网站| 久久久精品区二区三区| 亚洲国产欧美在线一区| 视频区欧美日本亚洲| 国产精品国产三级国产专区5o| 丰满迷人的少妇在线观看| 欧美xxⅹ黑人| 亚洲 国产 在线| 久久影院123| 成年美女黄网站色视频大全免费| 久久久久久人人人人人| 正在播放国产对白刺激| 性色av乱码一区二区三区2| 在线亚洲精品国产二区图片欧美| 国产精品欧美亚洲77777| 亚洲九九香蕉| 悠悠久久av| 午夜成年电影在线免费观看| 美女脱内裤让男人舔精品视频| 丝袜在线中文字幕| 欧美少妇被猛烈插入视频| 免费观看人在逋| 天天躁夜夜躁狠狠躁躁| 久久国产精品人妻蜜桃| 永久免费av网站大全| 日韩欧美国产一区二区入口| 青青草视频在线视频观看| 啪啪无遮挡十八禁网站| av在线app专区| 一区二区三区四区激情视频| 这个男人来自地球电影免费观看| 在线观看一区二区三区激情| 男男h啪啪无遮挡| 亚洲美女黄色视频免费看| 午夜福利乱码中文字幕| 夜夜骑夜夜射夜夜干| 最黄视频免费看| 久久中文看片网|