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

    瓦里安擴(kuò)充型動(dòng)態(tài)楔形板的蒙特卡羅模擬與驗(yàn)證

    2017-12-19 01:27:13戴振暉招什武楊耕朱琳張白霖朱遠(yuǎn)湖靳懷志王學(xué)濤張煜
    核技術(shù) 2017年12期
    關(guān)鍵詞:離軸射野電離室

    戴振暉 招什武 楊耕 朱琳 張白霖 朱遠(yuǎn)湖 靳懷志 王學(xué)濤 張煜

    ?

    瓦里安擴(kuò)充型動(dòng)態(tài)楔形板的蒙特卡羅模擬與驗(yàn)證

    戴振暉1招什武1楊耕1朱琳1張白霖1朱遠(yuǎn)湖1靳懷志1王學(xué)濤1張煜2

    1(廣東省中醫(yī)院 廣州 510006)2(南方醫(yī)科大學(xué) 廣州510515)

    基于蒙特卡羅方法模擬Varian直線加速器擴(kuò)充型動(dòng)態(tài)楔形板,建立動(dòng)態(tài)楔形板的質(zhì)量保證新工具。使用BEAMnrc的DYNJAWS組件模擬動(dòng)態(tài)楔形板,建立加速器源模型,用DOSXYZnrc程序計(jì)算體模內(nèi)的深度劑量和離軸劑量,獲得不同角度楔形板的楔形因子和離軸劑量分布曲線,與相同條件的測(cè)量結(jié)果進(jìn)行對(duì)比,驗(yàn)證動(dòng)態(tài)楔形板模型的精確性。蒙特卡羅模擬結(jié)果與測(cè)量結(jié)果取得了很好的一致,偏差小于2%。表明使用蒙特卡羅程序建立擴(kuò)充型動(dòng)態(tài)楔形板模型是可行的,能夠滿足臨床應(yīng)用的質(zhì)量保證需要。

    擴(kuò)充型動(dòng)態(tài)楔形板,蒙特卡羅模擬

    放療計(jì)劃設(shè)計(jì)過程中由于患者體表彎曲導(dǎo)致劑量分布不均勻,這種情況通常需要利用楔形板對(duì)射線束進(jìn)行修整,以獲得劑量均勻分布。擴(kuò)充型動(dòng)態(tài)楔形板(Enhanced Dynamic Wedge, EDW)不但具有傳統(tǒng)固定楔形板的功能,更具有其獨(dú)特的優(yōu)點(diǎn),EDW能夠提高靶區(qū)劑量的均勻性,同時(shí)降低正常組織的受照劑量,大大減少照射時(shí)間和跳數(shù),提高工作效率[1]。EDW射野產(chǎn)生的劑量分布精度關(guān)系到病人的療效,需要進(jìn)行專門的質(zhì)量保證[2],目前大部分放療單位使用Matrixx二維電離室矩陣進(jìn)行動(dòng)態(tài)楔形板質(zhì)量保證[3],本研究使用蒙特卡羅方法模擬瓦里安23EX直線加速器的EDW,首先采用BEAMnrc程序中的組件模塊DYNJAWS來模擬EDW,生成EDW源模型,然后用DOSXYZnrc程序計(jì)算體模內(nèi)劑量,并結(jié)合相同條件下實(shí)驗(yàn)測(cè)量結(jié)果分析楔形因子和離軸劑量分布曲線,對(duì)EDW模型進(jìn)行驗(yàn)證,開發(fā)一個(gè)獨(dú)立的EDW質(zhì)量保證工具。

    1 方法

    1.1 EDW工作原理

    EDW是通過加速器獨(dú)立準(zhǔn)直器運(yùn)動(dòng)和劑量率變化來實(shí)現(xiàn)特定劑量分布。瓦里安系統(tǒng)為每一檔能量的光子提供一套 60°動(dòng)態(tài)楔形板的Golden STT (Segmented Treatment Table)表,即最大野的加速器輸出跳數(shù)(Monitor Unit, MU)對(duì)應(yīng)準(zhǔn)直器運(yùn)動(dòng)擋塊位置的分配表。治療過程中STT表控制準(zhǔn)直器移動(dòng)速度和劑量率,生成相應(yīng)楔形角的EDW,用戶設(shè)定照射野大小和楔形角度后,系統(tǒng)從已存儲(chǔ)的Golden STT表中截取相應(yīng)的一段數(shù)據(jù),用它和相應(yīng)開野按一定比例合成所需角度的楔形野,最后結(jié)合所給加速器MU數(shù)生成治療用的STT表[4]。準(zhǔn)直器運(yùn)動(dòng)坐標(biāo)系統(tǒng)中靜止準(zhǔn)直器一邊為正,準(zhǔn)直器運(yùn)動(dòng)范圍為?20?9.5 cm,運(yùn)動(dòng)準(zhǔn)直器最后位置距離靜止準(zhǔn)直器0.5 cm。

    1.2 蒙特卡羅建模

    使用BEAMnrc模擬Varian23EX加速器6 MV光子束在二級(jí)準(zhǔn)直器上方生成相空間平面[5],該相空間平面作為EDW模擬的輸入源[6],使用BEAMnrc程序中的組件模塊DYNJAWS模擬EDW,在之前的研究中Kakakhel等[7]基于BEAMnrc程序模擬EDW,實(shí)現(xiàn)了EDW輸入文件的自動(dòng)化,探討了動(dòng)態(tài)和靜態(tài)兩種模式的優(yōu)劣。本研究直接使用更加精確的動(dòng)態(tài)模式模擬EDW,在BEAMnrc模擬過程中需要讀取包含EDW信息的輸入文件。輸入文件包括子野數(shù)量、子野概率、準(zhǔn)直器坐標(biāo)、-distance。準(zhǔn)直器的坐標(biāo):前表面正坐標(biāo)值(Front positive, FP)、后表面正坐標(biāo)值(Back positive, BP)、前表面負(fù)坐標(biāo)值(Front negative, FN)、后表面負(fù)坐標(biāo)值(Back negative, BN),準(zhǔn)直器固定。-distance包括-min和-max,-min表示源到準(zhǔn)直器前表面的距離;-max表示源到準(zhǔn)直器后表面的距離,如圖1所示。

    圖1 DYNJAWS坐標(biāo)

    子野概率用變量index描述,1是0?1之間生成的隨機(jī)數(shù),如果index(?1)<1≤index(),則使用第個(gè)子野。動(dòng)態(tài)模式中在第個(gè)和?1個(gè)子野之間進(jìn)行插值,表達(dá)連續(xù)的準(zhǔn)直器運(yùn)動(dòng)。選擇一個(gè)子野的概率是基于STT表計(jì)算得到,STT表定義了準(zhǔn)直器的位置和對(duì)應(yīng)的劑量,根據(jù)能量、楔形角和射野大小截取相關(guān)的STT表,從機(jī)器dynalog file中讀取準(zhǔn)直器位置,與STT表中準(zhǔn)直器的位置進(jìn)行比較,如果匹配則使用STT表中的劑量值作為選擇子野的概率;如果不匹配,則執(zhí)行插值,確定當(dāng)前子野的概率。在DYNJAWS模塊中讀取包含準(zhǔn)直器信息的SEQUENCE文件,完成不同角度的EDW的蒙卡建模[8]。

    BEAMnrc模擬含EDW的加速器機(jī)頭生成的相空間平面距離源44.5 cm,該相空間文件作為DosXYZnrc的輸入源文件,在DosXYZnrc程序的源選項(xiàng)中設(shè)置源到等中心距離(Distance from source to isocenter):55.5cm,最終源到體模表面的距離(Source Skin Distance, SSD)為100 cm。為提高計(jì)算效率,在BEAMnrc中設(shè)置電子截止能量CUT= 0.70MeV,光子截止能量CUT=0.01 MeV,電子射程截?cái)嗄芰縎AVE=2.0 MeV,采用方向韌致輻射分裂(Directional Bremsstrahlung Splitting, DBS)方差減少技術(shù),軔致輻射粒子分裂數(shù)BRSPL=1000,分裂野半徑30 cm,在DOSXYZnrc中設(shè)置s=40。BEAMnrc模擬粒子數(shù):5×108,DosXYZnrc模擬粒子數(shù):5×108;劑量計(jì)算體模規(guī)格:30 cm×30 cm× 20cm;單一體素大?。?.2 cm×0.2 cm×0.2 cm。通過DosXYZnrc計(jì)算體模內(nèi)的劑量,進(jìn)而計(jì)算不同大小對(duì)稱射野和非對(duì)稱射野的楔形因子,獲得10cm×10cm、20cm×20 cm射野,SSD 100 cm,水下10 cm處4個(gè)不同楔形角度的離軸劑量分布曲線和30°楔形角不同深度的離軸劑量分布曲線。

    1.3 實(shí)驗(yàn)測(cè)量

    楔形因子是在相同位置,照射相同跳數(shù),加和不加楔形板時(shí)射野中心軸上某一點(diǎn)劑量之比。楔形因子的測(cè)量條件:能量6 MV,SSD 100 cm。測(cè)量位置:射野中心軸水下10 cm處,測(cè)量對(duì)稱射野5cm×5 cm到20 cm×20 cm加EDW和不加EDW的劑量。測(cè)量過程使用FC65指型電離室,用Dose1劑量?jī)x接收測(cè)量數(shù)據(jù),動(dòng)態(tài)楔形板選擇1-IN方 向[9],電離室長(zhǎng)軸垂直于楔形方向和射野中心軸[10]。非對(duì)稱射野設(shè)計(jì):在楔形方向2固定于10 cm位置,1為0?20 cm,非楔形方向1:10 cm、2:10 cm,測(cè)量點(diǎn)位于射野的幾何中心0,0與1的位置關(guān)系如表1所示。

    表1 準(zhǔn)直器位置與射野中心的關(guān)系

    離軸Profiles測(cè)量條件:能量6 MV、SSD 100cm。射野大小:10 cm×10 cm、20 cm×20 cm,固體水下10 cm,動(dòng)態(tài)楔形板選擇2-OUT方向,用德國(guó)IBA (Ion Beam Applications)公司的Matrixx二維電離室矩陣(1020個(gè)電離室,電離室高:0.5 cm,直徑:0.45 cm,0.08 cm3的靈敏體積)測(cè)量不同楔形角的離軸劑量分布[11?12]。

    2 結(jié)果

    2.1 對(duì)稱射野的楔形因子

    對(duì)稱射野的楔形因子如圖2所示,楔形因子隨著射野大小增大而減小,10°?30°楔形板的楔形因子隨著射野增大楔形因子減小的幅度較大,45°和60°楔形板的楔形因子隨著射野增大楔形因子減小的幅度較小,5 cm×5 cm射野的楔形因子蒙特卡羅模擬值普遍大于測(cè)量值,模擬值與測(cè)量值平均偏差小于1%。

    圖2 對(duì)稱射野的楔形因子

    2.2 非對(duì)稱射野的楔形因子

    非對(duì)稱射野的楔形因子如圖3所示,非對(duì)稱射野的測(cè)量點(diǎn)位于射野幾何中心0,根據(jù)表1準(zhǔn)直器1在坐標(biāo)原點(diǎn)處,對(duì)應(yīng)的測(cè)量點(diǎn)0位于軸的5 cm處,非對(duì)稱射野的楔形因子隨著1準(zhǔn)直器的增大而減小,并且減小的幅度變小,模擬值與測(cè)量值平均偏差小于1%。

    圖3 非對(duì)稱射野的楔形因子

    2.3 離軸劑量分布曲線

    對(duì)比10 cm×10 cm、20 cm×20 cm的射野在水下10 cm處10°、25°、45°、60°楔形角的離軸劑量分布曲線,結(jié)果如圖4所示,模擬值和測(cè)量值取得了很好的一致,偏差小于2%,在高劑量區(qū)模擬值微低于測(cè)量值。相對(duì)劑量曲線以射野中心點(diǎn)的劑量歸一,其中蒙特卡羅劑量計(jì)算的誤差小于1%。

    圖4 不同角度EDW離軸劑量分布曲線 (a) 10 cm×10 cm,(b) 20 cm×20 cm

    對(duì)比10 cm×10 cm、20 cm×20 cm的射野在最大劑量深度(max)和10 cm深度處30°EDW的離軸劑量分布曲線結(jié)果如圖5所示,模擬值和測(cè)量值取得了很好的一致,偏差小于2%。相對(duì)劑量曲線以最大劑量深度射野中心點(diǎn)的劑量歸一,蒙特卡羅劑量計(jì)算的誤差小于1%。其中20 cm×20 cm射野的離軸劑量曲線在?10 ? 0區(qū)間內(nèi)劑量下降梯度較大,在0?10區(qū)間內(nèi)劑量下降梯度較小。

    圖5 30°EDW在不同深度的離軸劑量分布曲線 (a) 10 cm×10 cm,(b) 20 cm×20 cm

    3 討論

    本研究在BEAMnrc程序中使用DYNJAWS組件模塊模擬EDW,生成EDW源模型,使用DOSXYZnrc程序計(jì)算體模內(nèi)劑量,進(jìn)而計(jì)算楔形因子,獲得離軸劑量分布。設(shè)計(jì)幾種驗(yàn)證方法,用相同條件的測(cè)量結(jié)果驗(yàn)證EDW模型,模擬值與測(cè)量值取得了較好的一致。由于動(dòng)態(tài)楔形野各點(diǎn)的劑量由STT表決定,而且計(jì)算方程也是直接利用擋塊位置參數(shù),所以測(cè)量均在等中心平面進(jìn)行,便于保證測(cè)量點(diǎn)與計(jì)算點(diǎn)一致。對(duì)比蒙特卡羅模擬和測(cè)量的楔形因子,相同角度的動(dòng)態(tài)楔形板,楔形因子隨射野邊長(zhǎng)的增加平滑減小,且隨楔形角的增大這種減小更明顯,這主要?dú)w因于動(dòng)態(tài)楔形野的形成機(jī)制,動(dòng)態(tài)楔形板是利用加速器獨(dú)立準(zhǔn)直器運(yùn)動(dòng)和運(yùn)動(dòng)過程劑量率的變化來實(shí)現(xiàn)楔形劑量分布,對(duì)于某一角度楔形野先開野照射一定跳數(shù),然后準(zhǔn)直器逐漸關(guān)閉形成楔形效果,當(dāng)射野增大時(shí)必然會(huì)在楔形野關(guān)閉方向形成較高劑量,相應(yīng)射野中心軸上的累積劑量會(huì)降低,從而導(dǎo)致射野中心軸上楔形因子減小。對(duì)于小的照射野,蒙特卡羅模擬的楔形因子大于測(cè)量的楔形因子,這主要由于電離室測(cè)量存在誤差,在小野測(cè)量時(shí)電離室本身的影響就更加明顯[13]。相同射野時(shí)大楔形角EDW的離軸劑量分布曲線顯示蒙特卡羅模擬結(jié)果與測(cè)量結(jié)果偏差較大,尤其在高劑量區(qū)蒙特卡羅模擬結(jié)果普遍小于測(cè)量結(jié)果,由于楔形劑量分布是通過準(zhǔn)直器不斷關(guān)閉生成的,在高劑量區(qū)照射野變小,二維電離室矩陣對(duì)于小照射野的測(cè)量同樣存在誤差。相同楔形角大照射野的EDW離軸劑量分布曲線顯示蒙特卡羅模擬結(jié)果與測(cè)量結(jié)果偏差較大,這很可能因?yàn)樵诖笳丈湟吧尚ㄐ谓菚r(shí)準(zhǔn)直器步進(jìn)幅度變大,Matrixx內(nèi)部單個(gè)電離室的直徑是0.45 cm,電離室測(cè)量時(shí)取體積平均值,不同電離室之間采用了插值,導(dǎo)致測(cè)量結(jié)果誤差變大,而蒙特卡羅模擬過程取單一體素的劑量值,模擬過程體素大小與測(cè)量電離室大小不同,分辨率的差異會(huì)導(dǎo)致模擬值與測(cè)量值產(chǎn)生偏差,并且相應(yīng)楔形角度的STT表生成精確度對(duì)EDW的模擬精度也會(huì)產(chǎn)生影響。

    EDW在日常計(jì)劃設(shè)計(jì)過程中可以用來修正劑量分布,提高放療的效率,但需要定期執(zhí)行(Quality Assurance, QA)[14],對(duì)比蒙特卡羅模擬和實(shí)驗(yàn)測(cè)量的楔形因子和離軸劑量分布曲線,二者偏差小于2%,說明本文構(gòu)建的EDW蒙特卡羅模型是可靠的,為后續(xù)EDW放療計(jì)劃的蒙特卡羅驗(yàn)證提供了基礎(chǔ),可以作為獨(dú)立的EDW的QA工具。

    1 Koken P W, Heukelom S, Cuijpers J P,. On the practice of the clinical implementation of enhanced dynamic wedges[J]. Medical Dosimetry, 2003, 28(1): 13?19. DOI: 10.1016/S0958-3947(02)00140-1.

    2 Alaei P, Higgins P D. Performance evaluation and quality assurance of Varian enhanced dynamic wedges[J]. Journal of Applied Clinical Medical Physics, 2006, 7(1): 14?20. DOI: 10.1120/jacmp.2027.25367.

    3 于浪, 楊波, 李楠, 等. 應(yīng)用二維電離室矩陣進(jìn)行虛擬楔形板物理參數(shù)檢測(cè)[J]. 中華放射腫瘤學(xué)雜志, 2015, 24(3): 340?343. DOI: 10.3760/cma.j.issn.1004-4221.2015. 03.027. YU Lang, YANG Bo, LI Nan,. Dosimetry verification of virtual wedge using a two-dimensional ionization chamber array[J]. Chinese Journal of Radiation Oncology, 2015, 24(3): 340?343. DOI: 10.3760/cma.j. issn.1004-4221.2015.03.027.

    4 Varian Oncology Systems. C-series clinac enhanced dynamic wedge implementation guide[R]. Varian, 2002: 1?166.

    5 Rogers D W O, Faddegon B A, Ding G X,. BEAM: a Monte Carlo code to simulate radiotherapy treatment units[J]. Medical Physics, 1995, 22(5): 503?524. DOI: 10.1118/1.597552.

    6 熊綺麗, 石勇, 徐剛, 等. TrueBeam加速器多葉準(zhǔn)直器射野劑量學(xué)特性[J]. 核技術(shù), 2015, 38(10): 100203. DOI: 10.11889/j.0253-3219.2015.hjs.38.100203. XIONG Qili, SHI Yong, XU Gang,. Dosimetry characteristics of multi-leaf collimator field for TrueBeam[J]. Nuclear Techniques, 2015, 38(10): 100203. DOI: 10.11889/j.0253-3219.2015.hjs.38.100203.

    7 Kakakhel M B, Baveas E S, Fielding A L,. Validation and automation of the DYNJAWS component module of the BEAMnrc Monte Carlo code[J]. Australasian Physical and Engineering Sciences in Medicine, 2011, 34(1): 83?90. DOI: 10.1007/s13246-011-0060-x.

    8 Verhaegen F, Liu H H. Incorporating dynamic collimator motion in Monte Carlo simulations: an application in modelling a dynamic wedge[J]. Physics in Medicine and Biology, 2001, 46(2): 287?296. DOI: 10.1088/0031-9155/ 46/2/302.

    9 陳亮, 李玉民, 徐紅. 瓦里安加速器擴(kuò)充型動(dòng)態(tài)楔形板楔形因子特性研究[J]. 中華放射腫瘤學(xué)雜志, 2004, 13(1): 53?57. DOI: 10.3760/j.issn:1004-4221.2004.01. 013. CHEN Liang, LI Yumin, XU Hong. Characteristics of wedge factor of an enhanced dynamic wedge on Varian 600C accelerator[J]. Chinese Journal of Radiation Oncology, 2004, 13(1): 53?57. DOI: 10.3760/j.issn: 1004-4221.2004.01.013.

    10 Gibbons J P. Calculation of enhanced dynamic wedge factors for symmetric and asymmetric photon fields[J]. Medical Physics, 1998, 25(8): 1411?1418. DOI: 10.1118/ 1.598313.

    11 Oh S A, Kim S K, Kang M K,. Dosimetric verification of enhanced dynamic wedges by a 2D ion chamber array[J]. Journal of the Korean Physical Society, 2013, 63(11): 2215?2219. DOI: 10.1118/1. 4814791.

    12 Beavis A W, Weston S J, Whitton V J. Implementation of the Varian EDW into a commercial RTP system[J]. Physics in Medicine and Biology, 1996, 41(9): 1691?1704. DOI: 10.1088/0031-9155/41/9/009.

    13 Ahmad M, Deng J, Lund M W,. Clinical implementation of enhanced dynamic wedges into the pinnacle treatment planning system: Monte Carlo validation and patient-specific QA[J]. Physics in Medicine and Biology, 2009, 54(2): 447?465. DOI: 10.1088/0031-9155/54/2/018.

    14 Moeller J H, Leavitt D D, Klein E. The quality assurance of enhanced dynamic wedges[J]. Medical Dosimetry, 1997, 22(3): 241?246. DOI: 10.1016/S0958-3947(97) 00026-5.

    Monte Carlo simulation and validation of the enhanced dynamic wedge of Varian accelerator

    DAI Zhenhui1ZHAO Shiwu1YANG Geng1ZHU Lin1ZHANG Bailin1ZHU Yuanhu1JIN Huaizhi1WANG Xuetao1ZHANG Yu2

    1(Guangdong Provincial Hospital of Chinese Medicine, Guangzhou 510006, China) 2(Southern Medical University, Guangzhou 510515, China)

    Enhanced dynamic wedge (EDW) is commonly used to achieve dose uniformity in the target volume in radiotherapy. However, EDW must be verified before implementation in the treatment planning system.This paper is to simulate the EDW of Varian linear accelerator by Monte Carlo method, and to construct a new tool for quality assurance of EDW.EDW is simulated with the DYNJAWS module of BEAMnrc to build the source model. Depth dose and off-axis dose are calculated by use of DOSXYZnrc in a phantom to generate the EDW factor and profiles of EDW at different angles. The model is validated based on the comparison of Monte Carlo simulation and real measurement results.The simulation results are in good agreement with measurement results, and the deviation is less than 2%.It demonstrates that the EDW model based on Monte Carlo simulation is effective, which would achieve the requirement of quality assurance of EDW.

    EDW, Monte Carlo simulation

    DAI Zhenhui, male, born in 1984, graduated from Southern Medical University with a master’s degree in 2010, physicist, engineer, major in biomedical engineering

    WANG Xuetao, E-mail: wangxuetao0625@126.com

    2017-06-29,

    2017-08-23

    TL72

    10.11889/j.0253-3219.2017.hjs.40.120302

    戴振暉,男,1984年出生,2010年于南方醫(yī)科大學(xué)獲碩士學(xué)位,物理師,工程師,生物醫(yī)學(xué)工程專業(yè)

    王學(xué)濤,E-mail: wangxuetao0625@126.com

    2017-06-29,

    2017-08-23

    Supported by National Natural Science Foundation of China (No.61671230)

    國(guó)家自然科學(xué)基金(No.61671230)資助

    猜你喜歡
    離軸射野電離室
    離軸超構(gòu)透鏡設(shè)計(jì)與特性分析
    基于自由曲面雙波段離軸三反光學(xué)系統(tǒng)的優(yōu)化設(shè)計(jì)
    利用三維水箱測(cè)量的“環(huán)形機(jī)架”加速器“典型射線數(shù)據(jù)”驗(yàn)證研究
    PTW 729電離室矩陣不同驗(yàn)證方法用于宮頸癌術(shù)后調(diào)強(qiáng)放療計(jì)劃驗(yàn)證結(jié)果分析
    反射式紅外多波段準(zhǔn)直投影光學(xué)系統(tǒng)設(shè)計(jì)
    航空兵器(2019年2期)2019-05-30 00:00:00
    三維藍(lán)水箱(BPH)掃描測(cè)量系統(tǒng)在螺旋斷層加速器質(zhì)量控制檢測(cè)中的應(yīng)用
    寬譜段、動(dòng)態(tài)局部高分辨離軸主動(dòng)反射變焦系統(tǒng)
    DAVID系統(tǒng)探測(cè)MLC葉片位置誤差的能力測(cè)試與評(píng)估
    自動(dòng)多野排序執(zhí)行技術(shù)在適形調(diào)強(qiáng)放療中的應(yīng)用
    Precise直線加速器電離室的維護(hù)和故障排查
    精品久久久久久电影网| 天天躁狠狠躁夜夜躁狠狠躁| av一本久久久久| 少妇人妻 视频| 男女边摸边吃奶| 90打野战视频偷拍视频| 午夜免费男女啪啪视频观看| 国产av一区二区精品久久| 亚洲第一av免费看| 新久久久久国产一级毛片| 999精品在线视频| 一边摸一边做爽爽视频免费| 亚洲欧美精品综合一区二区三区| 国产在视频线精品| 久久久久人妻精品一区果冻| 女人精品久久久久毛片| av福利片在线| 这个男人来自地球电影免费观看 | 男人舔女人的私密视频| 亚洲av综合色区一区| 侵犯人妻中文字幕一二三四区| 大片免费播放器 马上看| 激情五月婷婷亚洲| 各种免费的搞黄视频| bbb黄色大片| 亚洲美女黄色视频免费看| 交换朋友夫妻互换小说| 亚洲第一青青草原| 精品一品国产午夜福利视频| 高清在线视频一区二区三区| 久久久久久免费高清国产稀缺| 国产精品嫩草影院av在线观看| 毛片一级片免费看久久久久| 午夜福利视频精品| 综合色丁香网| 丰满迷人的少妇在线观看| 大码成人一级视频| 亚洲国产av新网站| 一区二区三区精品91| 男女之事视频高清在线观看 | 麻豆乱淫一区二区| 欧美 亚洲 国产 日韩一| 色视频在线一区二区三区| 亚洲天堂av无毛| 国产精品无大码| 国产黄频视频在线观看| 超色免费av| 成人毛片60女人毛片免费| 日韩,欧美,国产一区二区三区| 日韩一区二区视频免费看| 一级爰片在线观看| 男女午夜视频在线观看| 国产亚洲精品第一综合不卡| 国产一级毛片在线| 国产一级毛片在线| 免费黄网站久久成人精品| 男人舔女人的私密视频| 国产精品蜜桃在线观看| 国产精品久久久久久精品古装| av女优亚洲男人天堂| 国产精品久久久久成人av| 午夜福利视频精品| 久久热在线av| 人人妻人人添人人爽欧美一区卜| av天堂久久9| 久久ye,这里只有精品| 成人漫画全彩无遮挡| 9色porny在线观看| 欧美精品av麻豆av| 精品一区在线观看国产| 中文字幕色久视频| 七月丁香在线播放| 国产麻豆69| 精品一区二区免费观看| 日韩不卡一区二区三区视频在线| 99热网站在线观看| 精品久久久久久电影网| 国产精品秋霞免费鲁丝片| 欧美日韩一级在线毛片| 日韩伦理黄色片| 99香蕉大伊视频| 欧美亚洲 丝袜 人妻 在线| 街头女战士在线观看网站| 欧美日韩视频精品一区| 免费人妻精品一区二区三区视频| 大香蕉久久成人网| 91国产中文字幕| 如日韩欧美国产精品一区二区三区| 99香蕉大伊视频| 七月丁香在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 大片电影免费在线观看免费| 青春草亚洲视频在线观看| av一本久久久久| 在线观看国产h片| 人人妻人人添人人爽欧美一区卜| 久久青草综合色| 国产成人av激情在线播放| 久久精品久久精品一区二区三区| 毛片一级片免费看久久久久| 久久精品国产综合久久久| 热99国产精品久久久久久7| 中文乱码字字幕精品一区二区三区| 亚洲国产最新在线播放| 王馨瑶露胸无遮挡在线观看| 午夜激情av网站| 亚洲一区中文字幕在线| 在线免费观看不下载黄p国产| 欧美精品亚洲一区二区| 啦啦啦啦在线视频资源| 在线免费观看不下载黄p国产| 久久久久国产精品人妻一区二区| 性色av一级| 丝袜在线中文字幕| 精品少妇黑人巨大在线播放| 午夜福利乱码中文字幕| 久热这里只有精品99| 老汉色av国产亚洲站长工具| 91精品三级在线观看| 国产一区二区激情短视频 | 国产精品亚洲av一区麻豆 | 丁香六月天网| 国产野战对白在线观看| 一区二区日韩欧美中文字幕| 中文字幕人妻丝袜一区二区 | 中国国产av一级| 久久影院123| 人人澡人人妻人| 国精品久久久久久国模美| 夜夜骑夜夜射夜夜干| 热99国产精品久久久久久7| 日本欧美国产在线视频| 久久人人爽人人片av| 妹子高潮喷水视频| 亚洲精品国产区一区二| 亚洲精品国产av蜜桃| 99久久人妻综合| 精品视频人人做人人爽| 免费看不卡的av| 五月开心婷婷网| 你懂的网址亚洲精品在线观看| 精品国产露脸久久av麻豆| 亚洲一区二区三区欧美精品| 午夜日本视频在线| 一区二区三区精品91| 美女主播在线视频| 日日撸夜夜添| 爱豆传媒免费全集在线观看| 久久ye,这里只有精品| 亚洲成人手机| 国产精品嫩草影院av在线观看| 中文乱码字字幕精品一区二区三区| 久久人人爽人人片av| 女的被弄到高潮叫床怎么办| 国产野战对白在线观看| 国产亚洲精品第一综合不卡| 久久久欧美国产精品| 精品少妇黑人巨大在线播放| 黑丝袜美女国产一区| 人成视频在线观看免费观看| 18在线观看网站| 亚洲国产av影院在线观看| 美女扒开内裤让男人捅视频| 亚洲,一卡二卡三卡| 亚洲综合精品二区| 亚洲五月色婷婷综合| 丰满饥渴人妻一区二区三| 久久久久久人妻| 丁香六月天网| 国产精品麻豆人妻色哟哟久久| 精品国产一区二区三区四区第35| 国产成人午夜福利电影在线观看| 久久精品国产亚洲av涩爱| 国产成人精品在线电影| 欧美日韩综合久久久久久| 国产在线免费精品| 免费观看av网站的网址| 18在线观看网站| 99精品久久久久人妻精品| 一本色道久久久久久精品综合| 自拍欧美九色日韩亚洲蝌蚪91| 91老司机精品| 在线观看免费高清a一片| 9热在线视频观看99| 午夜精品国产一区二区电影| 国产淫语在线视频| 宅男免费午夜| 伦理电影大哥的女人| 母亲3免费完整高清在线观看| 中文字幕人妻熟女乱码| 中国三级夫妇交换| 一级毛片电影观看| 久久青草综合色| 日韩 亚洲 欧美在线| 飞空精品影院首页| 久久精品熟女亚洲av麻豆精品| 久久av网站| 99国产综合亚洲精品| 久久性视频一级片| av福利片在线| 黑人猛操日本美女一级片| 黄色毛片三级朝国网站| 一级黄片播放器| 国产亚洲av片在线观看秒播厂| 赤兔流量卡办理| 国产成人精品在线电影| 一区二区三区精品91| 国产精品人妻久久久影院| 午夜免费鲁丝| 亚洲精品国产av成人精品| 婷婷色av中文字幕| 免费黄频网站在线观看国产| 1024香蕉在线观看| 黄色视频不卡| 国产女主播在线喷水免费视频网站| 成人国产麻豆网| 人人澡人人妻人| 婷婷色综合www| 亚洲精品国产av蜜桃| 日韩一区二区视频免费看| 精品酒店卫生间| 午夜福利一区二区在线看| 亚洲av日韩精品久久久久久密 | 成年av动漫网址| 考比视频在线观看| 99国产综合亚洲精品| 中文字幕人妻熟女乱码| 免费久久久久久久精品成人欧美视频| 观看美女的网站| 成年人午夜在线观看视频| 18禁动态无遮挡网站| 亚洲成人国产一区在线观看 | 国产免费现黄频在线看| 男女边摸边吃奶| 777米奇影视久久| 午夜福利视频精品| 欧美成人午夜精品| av在线老鸭窝| 亚洲久久久国产精品| 中文天堂在线官网| 久久久精品国产亚洲av高清涩受| 丰满迷人的少妇在线观看| 欧美日韩av久久| 一级爰片在线观看| 久久女婷五月综合色啪小说| 亚洲人成电影观看| 男女边吃奶边做爰视频| 亚洲自偷自拍图片 自拍| 久久久久久久国产电影| 欧美黑人精品巨大| 在线天堂最新版资源| 日韩一本色道免费dvd| 成人国语在线视频| 中文字幕av电影在线播放| 国产乱来视频区| 少妇精品久久久久久久| 看十八女毛片水多多多| 中文欧美无线码| 国产成人啪精品午夜网站| 嫩草影院入口| 狠狠精品人妻久久久久久综合| 久久久精品免费免费高清| 伦理电影大哥的女人| 久久99精品国语久久久| 欧美日韩亚洲国产一区二区在线观看 | 一级片免费观看大全| 亚洲国产av新网站| 久久热在线av| 老司机在亚洲福利影院| 新久久久久国产一级毛片| 久久鲁丝午夜福利片| 中文欧美无线码| 亚洲av电影在线观看一区二区三区| 超碰成人久久| 熟妇人妻不卡中文字幕| 国产精品香港三级国产av潘金莲 | 少妇精品久久久久久久| 黄色视频不卡| 汤姆久久久久久久影院中文字幕| 2018国产大陆天天弄谢| 国产精品熟女久久久久浪| 亚洲国产欧美日韩在线播放| 男女床上黄色一级片免费看| 大香蕉久久成人网| 国产亚洲午夜精品一区二区久久| 一级毛片我不卡| 中文字幕亚洲精品专区| 人人妻人人爽人人添夜夜欢视频| 天天躁夜夜躁狠狠久久av| 一边亲一边摸免费视频| 亚洲自偷自拍图片 自拍| 波多野结衣av一区二区av| 亚洲av日韩精品久久久久久密 | 岛国视频午夜一区免费看| 国产精品自产拍在线观看55亚洲| 18禁裸乳无遮挡免费网站照片 | 精品久久久久久久毛片微露脸| 丝袜人妻中文字幕| 性欧美人与动物交配| 啦啦啦韩国在线观看视频| 亚洲五月色婷婷综合| 国产xxxxx性猛交| 国内毛片毛片毛片毛片毛片| 无人区码免费观看不卡| 欧美一级a爱片免费观看看 | 这个男人来自地球电影免费观看| 老司机福利观看| 国产精品乱码一区二三区的特点 | 精品熟女少妇八av免费久了| 亚洲avbb在线观看| 最近最新中文字幕大全电影3 | 夜夜爽天天搞| 法律面前人人平等表现在哪些方面| 美女国产高潮福利片在线看| 久久国产精品男人的天堂亚洲| 国产伦人伦偷精品视频| 国产单亲对白刺激| 免费在线观看视频国产中文字幕亚洲| 久久久国产欧美日韩av| 国产精品精品国产色婷婷| 欧美黄色淫秽网站| 咕卡用的链子| 久久国产精品男人的天堂亚洲| cao死你这个sao货| 成人国语在线视频| 欧美 亚洲 国产 日韩一| 亚洲国产精品sss在线观看| 大型av网站在线播放| 美国免费a级毛片| 大型黄色视频在线免费观看| 国产av精品麻豆| 国产精品秋霞免费鲁丝片| 国产成+人综合+亚洲专区| 欧美在线黄色| 啪啪无遮挡十八禁网站| √禁漫天堂资源中文www| 手机成人av网站| 国产精品香港三级国产av潘金莲| 欧美成人一区二区免费高清观看 | 最近最新中文字幕大全电影3 | 九色国产91popny在线| 少妇 在线观看| 亚洲av成人av| 99久久国产精品久久久| 久久精品亚洲精品国产色婷小说| 亚洲美女黄片视频| 69av精品久久久久久| 久久久久精品国产欧美久久久| 少妇粗大呻吟视频| 亚洲人成77777在线视频| 熟女少妇亚洲综合色aaa.| 侵犯人妻中文字幕一二三四区| 天堂动漫精品| 免费在线观看亚洲国产| 一级毛片高清免费大全| 免费少妇av软件| 久久天躁狠狠躁夜夜2o2o| 欧美日韩亚洲国产一区二区在线观看| 我的亚洲天堂| 国产精品久久电影中文字幕| 女警被强在线播放| 好男人在线观看高清免费视频 | 香蕉国产在线看| 黄色毛片三级朝国网站| 午夜福利18| 大型av网站在线播放| av天堂久久9| 日韩精品青青久久久久久| 亚洲国产精品成人综合色| 亚洲情色 制服丝袜| 午夜福利高清视频| 免费人成视频x8x8入口观看| 国产人伦9x9x在线观看| 一区二区三区高清视频在线| 国产精品av久久久久免费| 国内久久婷婷六月综合欲色啪| 欧美日韩黄片免| 叶爱在线成人免费视频播放| 亚洲欧美精品综合久久99| 激情视频va一区二区三区| 精品久久久久久久毛片微露脸| 自拍欧美九色日韩亚洲蝌蚪91| 搞女人的毛片| 国产91精品成人一区二区三区| 三级毛片av免费| 亚洲欧美精品综合一区二区三区| 国产精品日韩av在线免费观看 | 欧美精品亚洲一区二区| 人人妻人人爽人人添夜夜欢视频| 女警被强在线播放| 熟妇人妻久久中文字幕3abv| 婷婷六月久久综合丁香| 真人一进一出gif抽搐免费| 一本大道久久a久久精品| 国产私拍福利视频在线观看| 久99久视频精品免费| 两性午夜刺激爽爽歪歪视频在线观看 | 天天添夜夜摸| 欧美另类亚洲清纯唯美| 99久久久亚洲精品蜜臀av| av视频在线观看入口| 午夜免费鲁丝| 国产三级在线视频| 18美女黄网站色大片免费观看| 99精品在免费线老司机午夜| 巨乳人妻的诱惑在线观看| 久久性视频一级片| 韩国av一区二区三区四区| 亚洲视频免费观看视频| 国内精品久久久久久久电影| 国产精品99久久99久久久不卡| 天堂动漫精品| 欧美一级毛片孕妇| 麻豆成人av在线观看| 女性生殖器流出的白浆| 欧美国产日韩亚洲一区| 精品人妻1区二区| 久久精品成人免费网站| 看片在线看免费视频| 两个人免费观看高清视频| aaaaa片日本免费| 亚洲熟女毛片儿| 丝袜美腿诱惑在线| 亚洲天堂国产精品一区在线| 亚洲成人精品中文字幕电影| 禁无遮挡网站| 熟妇人妻久久中文字幕3abv| 美女高潮喷水抽搐中文字幕| 国产欧美日韩综合在线一区二区| 搞女人的毛片| 国产成年人精品一区二区| 狂野欧美激情性xxxx| 欧美亚洲日本最大视频资源| 99国产综合亚洲精品| 国产高清videossex| 韩国精品一区二区三区| 99国产精品一区二区蜜桃av| 日日夜夜操网爽| 亚洲欧美日韩无卡精品| 黑丝袜美女国产一区| 亚洲avbb在线观看| 精品国产超薄肉色丝袜足j| 久久中文字幕一级| 多毛熟女@视频| 亚洲最大成人中文| 91麻豆精品激情在线观看国产| 此物有八面人人有两片| 久久性视频一级片| 精品国产国语对白av| 精品乱码久久久久久99久播| 成人国语在线视频| 国产成人欧美| 亚洲国产精品成人综合色| 丝袜美足系列| 国产又色又爽无遮挡免费看| 久久久久九九精品影院| av电影中文网址| 亚洲精品美女久久av网站| 91麻豆av在线| 一进一出抽搐gif免费好疼| 一级片免费观看大全| 91九色精品人成在线观看| 91精品三级在线观看| 欧美另类亚洲清纯唯美| 91字幕亚洲| 久久久久九九精品影院| 中文字幕精品免费在线观看视频| 老汉色av国产亚洲站长工具| 午夜久久久在线观看| 18禁观看日本| 18禁裸乳无遮挡免费网站照片 | 高清在线国产一区| 波多野结衣高清无吗| 欧美成人一区二区免费高清观看 | 亚洲av成人一区二区三| 深夜精品福利| 熟女少妇亚洲综合色aaa.| 最好的美女福利视频网| 精品欧美国产一区二区三| 国语自产精品视频在线第100页| 99国产综合亚洲精品| 看免费av毛片| 国产私拍福利视频在线观看| 性少妇av在线| 国产一区二区三区视频了| 久久精品91无色码中文字幕| netflix在线观看网站| 中文字幕人妻丝袜一区二区| 亚洲avbb在线观看| 国产av在哪里看| 女警被强在线播放| 亚洲中文字幕日韩| 18禁美女被吸乳视频| 久久精品国产综合久久久| 国产午夜精品久久久久久| 亚洲第一电影网av| 国产亚洲精品第一综合不卡| 在线观看日韩欧美| 日韩中文字幕欧美一区二区| 久久久久久久久久久久大奶| 精品久久久精品久久久| 精品乱码久久久久久99久播| 午夜成年电影在线免费观看| 久久人人精品亚洲av| 国产精品秋霞免费鲁丝片| 黄色 视频免费看| 嫁个100分男人电影在线观看| 777久久人妻少妇嫩草av网站| 久久精品国产亚洲av高清一级| 在线国产一区二区在线| 亚洲熟妇中文字幕五十中出| 天堂影院成人在线观看| 男人舔女人下体高潮全视频| 女人被狂操c到高潮| 神马国产精品三级电影在线观看 | 午夜视频精品福利| av网站免费在线观看视频| 国产免费男女视频| 19禁男女啪啪无遮挡网站| 久久久久久亚洲精品国产蜜桃av| 午夜激情av网站| 国产精品 国内视频| 一二三四在线观看免费中文在| 91av网站免费观看| 成人精品一区二区免费| 亚洲精品美女久久av网站| 亚洲在线自拍视频| 18禁国产床啪视频网站| 高清在线国产一区| 日韩欧美一区视频在线观看| 色综合欧美亚洲国产小说| 精品人妻1区二区| 热99re8久久精品国产| 嫩草影视91久久| 黄色片一级片一级黄色片| 韩国精品一区二区三区| 99香蕉大伊视频| 丝袜在线中文字幕| 好男人在线观看高清免费视频 | 九色亚洲精品在线播放| 久久久久久久精品吃奶| 无限看片的www在线观看| 欧美激情高清一区二区三区| 国产熟女xx| 大陆偷拍与自拍| 国产熟女xx| 久久中文字幕一级| 久久精品亚洲精品国产色婷小说| 丝袜美腿诱惑在线| 国产免费男女视频| 亚洲欧美一区二区三区黑人| 日韩国内少妇激情av| 亚洲av电影在线进入| 9191精品国产免费久久| 亚洲性夜色夜夜综合| 国产免费男女视频| 国内毛片毛片毛片毛片毛片| 91九色精品人成在线观看| 少妇裸体淫交视频免费看高清 | 日本欧美视频一区| 久久草成人影院| 精品久久久久久,| 91麻豆av在线| 精品国产一区二区三区四区第35| 性欧美人与动物交配| 亚洲国产毛片av蜜桃av| 亚洲第一欧美日韩一区二区三区| 在线观看一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 美女午夜性视频免费| 午夜成年电影在线免费观看| 日本欧美视频一区| 日本三级黄在线观看| 一区福利在线观看| 国产高清videossex| 免费在线观看完整版高清| 伦理电影免费视频| 国产私拍福利视频在线观看| 50天的宝宝边吃奶边哭怎么回事| 88av欧美| 看免费av毛片| 国产一区二区三区综合在线观看| 国产1区2区3区精品| 亚洲九九香蕉| av电影中文网址| 欧美日韩乱码在线| 国产精品乱码一区二三区的特点 | 亚洲午夜理论影院| 嫩草影院精品99| 757午夜福利合集在线观看| 欧美日韩乱码在线| 欧美绝顶高潮抽搐喷水| 99国产精品一区二区蜜桃av| 老熟妇乱子伦视频在线观看| 这个男人来自地球电影免费观看| 老司机在亚洲福利影院| 国产成人影院久久av| 黄网站色视频无遮挡免费观看| 青草久久国产| 久久人人97超碰香蕉20202| 两人在一起打扑克的视频| 久久国产精品男人的天堂亚洲| 少妇的丰满在线观看| 一级a爱片免费观看的视频| 亚洲精品国产一区二区精华液| 欧美午夜高清在线| 人妻丰满熟妇av一区二区三区| 成人免费观看视频高清| 国产熟女xx| cao死你这个sao货| 一边摸一边做爽爽视频免费| 给我免费播放毛片高清在线观看| 999久久久国产精品视频| 亚洲一区高清亚洲精品| 婷婷六月久久综合丁香|