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

    LUNA全身伽馬刀輸出因子的蒙特卡羅模擬研究*

    2016-10-19 11:48:27袁繼龍程金生
    中國醫(yī)學(xué)裝備 2016年9期
    關(guān)鍵詞:伽馬刀蒙特卡羅模擬計(jì)算

    袁繼龍 程金生*

    LUNA全身伽馬刀輸出因子的蒙特卡羅模擬研究*

    袁繼龍①程金生①*

    目的:利用蒙特卡羅模擬方法,對LUNA全身伽馬刀重要?jiǎng)┝繉W(xué)參數(shù)輸出因子進(jìn)行研究。方法:采集廠商提供的技術(shù)數(shù)據(jù)文檔,建立LUNA全身伽馬刀蒙特卡羅模型,模型驗(yàn)證后結(jié)合廠商測量輸出因子所用探測器,設(shè)置蒙特卡羅模型相關(guān)參數(shù)進(jìn)行模擬計(jì)算,比較分析廠商提供的輸出因子與蒙特卡羅模擬計(jì)算出輸出因子的差異。結(jié)果:在3%偏差范圍內(nèi),不同準(zhǔn)直器條件下,廠商提供的輸出因子和蒙特卡羅模擬計(jì)算得到的輸出因子符合性好,并進(jìn)一步分析了相關(guān)差異的原因。結(jié)論:LUNA全身伽馬刀輸出因子的蒙特卡羅模擬研究為該設(shè)備的臨床使用和質(zhì)量控制提供了參考和依據(jù)。

    LUNA全身伽馬刀;蒙特卡羅模擬;輸出因子

    [First-author’s address] Key Laboratory, Institute for Radiological Protection and Nuclear Safety, Chinese Center for Disease Control and Prevention, Beijing 100088, China.

    LUNA全身伽馬刀是全身立體定向放射治療專用設(shè)備,目前全國已經(jīng)有40余臺(tái)LUNA全身伽馬刀應(yīng)用于臨床治療。LUNA全身伽馬刀有42顆60Co放射源在源體中沿著臨床X軸和Y軸平面呈扇形對稱分布,通過準(zhǔn)直器聚焦在等中心位置,LUNA全身伽瑪?shù)杜鋫?組不同規(guī)格,即6 mm×6 mm、8 mm×8 mm、14 mm×14 mm、14 mm×20 mm、14 mm×40 mm和14 mm×60 mm,分別對應(yīng)1號(hào)、2號(hào)、3號(hào)、4號(hào)、5號(hào)和6號(hào)準(zhǔn)直器的矩形孔準(zhǔn)直器。

    在臨床治療時(shí)放射源體連同準(zhǔn)直器可以繞著臨床Z軸進(jìn)行旋轉(zhuǎn)進(jìn)行非共面的聚焦照射,最大旋轉(zhuǎn)角度是以臨床Y軸為基準(zhǔn)的180o。6個(gè)型號(hào)的準(zhǔn)直器以滿足頭部、體部腫瘤的體積大小的立體定向放射治療照射需要進(jìn)行安裝使用,所有準(zhǔn)直器條件下在等中心焦點(diǎn)處形成的聚焦野,皆非傳統(tǒng)放射治療10 cm×10 cm的標(biāo)準(zhǔn)照射野,結(jié)合LUNA全身伽馬刀的臨床動(dòng)態(tài)旋轉(zhuǎn)照射治療特點(diǎn),其劑量學(xué)研究有著極其重要的臨床治療和質(zhì)量控制意義。由于LUNA全身伽馬刀主要在國內(nèi)使用,其相關(guān)劑量學(xué)研究報(bào)道很少[1-6]。為此,本研究結(jié)合蒙特卡羅模擬的優(yōu)勢進(jìn)行建模,比對廠商測試結(jié)果,對LUNA全身伽馬刀的輸出因子這一重要?jiǎng)┝繉W(xué)參數(shù)進(jìn)行研究,為LUNA全身伽馬刀放射治療的質(zhì)量控制提供參考和依據(jù)[7-10]。

    1 資料與方法

    1.1設(shè)備與材料

    Luna-260型LUNA全身伽馬刀(深圳一體公司生產(chǎn));EGSnrc蒙特卡羅模擬程序工具包(加拿大國家研究中心開發(fā))。

    1.2研究方法

    以LUNA全身伽馬刀為模型,根據(jù)廠商提供的技術(shù)文檔,將放射源、放置放射源的源體、預(yù)準(zhǔn)直器以及終準(zhǔn)直器等劑量學(xué)相關(guān)部件的幾何尺寸、相對空間位置分布、對應(yīng)材料及動(dòng)態(tài)旋轉(zhuǎn)模式作為模型的輸入?yún)?shù)進(jìn)行蒙特卡羅模擬計(jì)算。

    1.3蒙特卡羅模擬計(jì)算

    蒙特卡羅模擬計(jì)算使用通用EGSnrc蒙特卡羅模擬程序工具包,在BEAMnrc模塊3號(hào)源項(xiàng)中設(shè)置建模部分的相關(guān)部件幾何尺寸、相對空間位置分布和對應(yīng)材料,抽樣γ粒子數(shù)為2×109,獲得終準(zhǔn)直器出口位置的源相空間文件,其主要包括射線到達(dá)設(shè)定位置區(qū)域的位置、角動(dòng)量及對應(yīng)的能量等分布信息。

    (1)以此源相空間文件為輸入文件在DOSXYZnrc模塊中進(jìn)行設(shè)置,該部分設(shè)置分兩步進(jìn)行:①驗(yàn)證蒙特卡羅模型,對應(yīng)實(shí)際測量情況,在8號(hào)源項(xiàng)中設(shè)置相關(guān)源分布進(jìn)行靜態(tài)聚焦照射,模體體素為1 mm×1 mm×1 mm,抽樣γ粒子數(shù)為2×109,分析模擬結(jié)果3ddose文件,分別提取臨床X軸、Y軸和Z軸上的聚焦野劑量分布信息,計(jì)算劑量分布半高寬,與廠商技術(shù)文件提供相關(guān)信息進(jìn)行比對分析,驗(yàn)證模型[11];②在驗(yàn)證蒙特卡羅模型的基礎(chǔ)上,在20號(hào)源項(xiàng)中根據(jù)源體旋轉(zhuǎn)拉弧的角度和源的活度實(shí)際分布進(jìn)行設(shè)置。

    (2)在等中心處設(shè)置8 cm×8 cm×8 cm的有機(jī)玻璃模體,等中心位置體素根據(jù)不同準(zhǔn)直器條件下廠商測量輸出因子所用探測器靈敏體積進(jìn)行設(shè)置;①在3號(hào)、4號(hào)、5號(hào)及6號(hào)準(zhǔn)直器條件下,使用0.015 cc電離室進(jìn)行測量,其靈敏體積近似為2 mm×2 mm×5 mm,1號(hào)和2號(hào)準(zhǔn)直器條件下,使用Diode半導(dǎo)體探測器進(jìn)行測量;②其靈敏體積近似為1 mm×1 mm×1 mm,對這些感興趣區(qū)域體素和0o~180o的旋轉(zhuǎn)聚焦照射進(jìn)行編程設(shè)置;③在每個(gè)準(zhǔn)直器條件下,抽樣γ粒子數(shù)均為2×109,分析模擬結(jié)果3ddose文件,獲取等中心處劑量值,相應(yīng)劑量值(歸一到每一個(gè)粒子)偏差均控制在≤0.5%范圍內(nèi),偏差計(jì)算為公式1:

    [(MC模擬-廠商數(shù)據(jù))÷MC模擬]×100%(1)

    (3)以最大尺寸的6號(hào)準(zhǔn)直器為基準(zhǔn),分別計(jì)算其余準(zhǔn)直器條件下的輸出因子,并與廠商所提供技術(shù)文件中的相應(yīng)輸出因子進(jìn)行比較分析。

    2 結(jié)果

    (1)基于EGSnrc的蒙特卡羅模擬環(huán)境在Ubuntu12.04 32位系統(tǒng)下搭建,運(yùn)行CPU為4核InterCoreTM760 2.8 GHz,有效利用內(nèi)存大約為3.0 GB,每個(gè)準(zhǔn)直器條件下BEAMnrc模塊運(yùn)行為4.5 h左右,每個(gè)準(zhǔn)直器條件下DOSXYZnrc模塊運(yùn)行為8 h左右。LUNA全身伽馬刀靜態(tài)聚焦照射時(shí),等中心處聚焦野在X軸上的劑量分布半高寬見表1。

    表1 X軸上劑量分布半高寬(mm)

    LUNA全身伽馬刀靜態(tài)聚焦照射時(shí),等中心處聚焦野在Y軸上的劑量分布半高寬見表2。

    表2 Y軸上劑量分布半高寬(mm)

    LUNA全身伽馬刀靜態(tài)聚焦照射時(shí),等中心處聚焦野在Z軸上的劑量分布半高寬見表3。

    表3 Z軸上劑量分布半高寬(mm)

    (2)根據(jù)不同準(zhǔn)直器條件,分別進(jìn)行對應(yīng)的焦點(diǎn)處感興趣區(qū)域的體素設(shè)置,采用蒙特卡羅模擬計(jì)算,其等中心處感興趣區(qū)域劑量值結(jié)果歸一到每個(gè)抽樣模擬計(jì)算的粒子上,見表4。

    表4 蒙特卡羅模擬等中心處感興趣區(qū)域劑量值(×10-17cGy)

    (3)由蒙特卡羅模擬等中心處感興趣區(qū)域的劑量值計(jì)算LUNA全身伽馬刀不同準(zhǔn)直器條件下的輸出因子,并與廠商提供基于測量計(jì)算出的輸出因子進(jìn)行比較,其相應(yīng)偏差見表5。

    3 討論

    本研究目的是在蒙特卡羅模擬計(jì)算輔助下分析LUNA全身伽馬刀輸出因子,涉及的聚焦野焦點(diǎn)吸收劑量全部為小野模式,相較于傳統(tǒng)的放射治療射野劑量分布的測量,小野模式需要考慮靈敏體積較大的探測器,如0.6 cc及0.125 cc電離室,在聚焦野內(nèi)測量時(shí)存在側(cè)向電子欠平衡效應(yīng)及體積效應(yīng)問題,這兩個(gè)效應(yīng)會(huì)影響到測量結(jié)果的準(zhǔn)確性,而且LUNA全身伽馬刀所有準(zhǔn)直器條件下的聚焦野都是非標(biāo)準(zhǔn)野,直接用靈敏體積小的電離室如0.015 cc電離室和其他探測器如Diode半導(dǎo)體探測器和MOSFET場效應(yīng)管探測器測量前不但需要在標(biāo)準(zhǔn)輻射場下準(zhǔn)確刻度修正,而且需要預(yù)先了解所測聚焦野的劑量分布情況,以便在不同的聚焦野劑量分布條件下準(zhǔn)確使用合適的探測器。

    表5 廠商提供輸出因子與蒙特卡羅計(jì)算輸出因子比較

    本研究通過廠商技術(shù)文檔提供的相關(guān)參數(shù),建立LUNA全身伽馬刀蒙特卡羅模型,模擬計(jì)算得到的等中心焦點(diǎn)處三維劑量分布,進(jìn)一步提取出沿著臨床X軸、Y軸及Z軸上劑量分布半高寬和廠商提供的設(shè)計(jì)尺寸進(jìn)行對比,表1,表2和表3的結(jié)果對比表明,除最小準(zhǔn)直器條件下模擬結(jié)果大于設(shè)計(jì)允許偏差外,原因應(yīng)該是模體體素設(shè)置為1 mm×1 mm×1 mm,相較于5 mm的劑量分布半高寬相對較大,而且分析模擬計(jì)算結(jié)果時(shí),采用線性插值方法獲取模擬劑量分布半高寬的方法所造成,此差異在隨后研究中將采用更小的體素來分析該準(zhǔn)直器條件下聚焦野的劑量分布[12]。其余結(jié)果均在聚焦野設(shè)計(jì)允許偏差之內(nèi),因此本研究建立的LUNA全身伽馬刀蒙特卡羅模型得到驗(yàn)證。

    利用驗(yàn)證過的蒙特卡羅模型,參考蒙特卡羅模擬出劑量分布半高寬,可以選擇靈敏體積合適的探測器進(jìn)行輸出因子的測量。本研究首先為了驗(yàn)證廠商給出的輸出因子,按照廠商測試所采用探測器的靈敏體積在蒙特卡羅模型的模體中進(jìn)行焦點(diǎn)處感興趣區(qū)域設(shè)置進(jìn)行模擬計(jì)算[13-15]。表5顯示,廠商測試得到的輸出因子和本研究蒙特卡羅模型模擬計(jì)算得到的輸出因子的差異,其中4號(hào)、5號(hào)準(zhǔn)直器的輸出因子偏差很小,3號(hào)準(zhǔn)直器偏差較大,應(yīng)該是廠商測試所用0.015 cc電離室在此準(zhǔn)直器條件下的聚焦野內(nèi)出現(xiàn)體積平均效應(yīng),造成測量值偏低,測量結(jié)果和模擬結(jié)果出現(xiàn)如此偏差[16]。1號(hào)、2號(hào)準(zhǔn)直器條件下,輸出因子偏差相較于3號(hào)準(zhǔn)直器條件偏差小很多,但與5號(hào)和6號(hào)準(zhǔn)直器條件相比較偏差明顯大些,而且此偏差均表現(xiàn)為測試結(jié)果大于蒙特卡羅模擬結(jié)果,廠商測試1號(hào)和2號(hào)準(zhǔn)直器用的探測器是Diode半導(dǎo)體探測器,其靈敏體積的半導(dǎo)體有效原子序數(shù)大于測試及模擬所用有機(jī)玻璃模體,在測量過程中會(huì)有一定程度的過響應(yīng)效應(yīng)造成測量數(shù)據(jù)偏大,從而導(dǎo)致測試和模擬輸出因子的偏差。

    通過建立LUNA全身伽馬刀的蒙特卡羅模擬模型,可以更容易的分析LUNA全身伽馬刀的焦點(diǎn)處三維劑量分布情況,輔助解決現(xiàn)實(shí)中不容易實(shí)驗(yàn)測量的劑量學(xué)參數(shù)的問題,輔助實(shí)驗(yàn)測量中合適探測器的選擇,相應(yīng)復(fù)雜射野情況下修正因子的獲取和測試結(jié)果的科學(xué)合理性分析[17-18]。

    4 結(jié)語

    本研究介紹了建立LUNA全身伽馬刀蒙特卡羅模型的方法,在3%的偏差范圍內(nèi),廠商提供的輸出因子和蒙特卡羅模擬結(jié)果符合性好,并進(jìn)一步分析了輸出因子的測量和模擬結(jié)果在不同準(zhǔn)直器條件下偏差的原因。本研究為LUNA全身伽馬刀設(shè)備的臨床使用及質(zhì)量控制提供了參考和依據(jù)。同時(shí),下一步將進(jìn)行LUNA全身伽馬刀的劑量學(xué)研究。

    [1]張曉峰,喬寶龍.有關(guān)伽瑪?shù)兜募夹g(shù)問題探討[J].中國腫瘤,2007,16(7):582-584.

    [2]黃禹,端正澄,朱國力,等.全身伽瑪?shù)稒C(jī)械結(jié)構(gòu)設(shè)計(jì)與制造[J].機(jī)械制造,2006,44(7):28-30.

    [3]張重魁,莫軍揚(yáng).體部伽瑪?shù)杜R床應(yīng)用及進(jìn)展[J].實(shí)用癌癥進(jìn)展,2012,27(1):107-110.

    [4]康靜波.體部γ刀治療技術(shù)的進(jìn)展[J].武警醫(yī)學(xué),2009,20(6):485-488.

    [5]劉涵,郭吉豐,羅浩,等.伽瑪?shù)秳┝繄鰷y量分析及驗(yàn)證系統(tǒng)的設(shè)計(jì)和研究[J].中國醫(yī)療器械雜志,2009,29(6):402-405.

    [6]廖福錫,夏邦傳,莫莉,等.SGS-I立體定向伽瑪線全身放療系統(tǒng)源體定位系統(tǒng)的改進(jìn)[J].2010,17(5):12-13.

    [7]彭柳芬,孔令人,范善翔.Monte Carlo模擬在伽瑪?shù)吨委熤械膽?yīng)用[J].醫(yī)療設(shè)備信息,2006,26(1):9-11.

    [8]Giuseppe Battistoni,F(xiàn)abrizio Cappucci,Nicola Bertolino,et al.FLUKA Monte Carlo simulation for the Leksell Gamma Knife Perfexion radiosurgery system:Homogeneous media[J]. Phys Med,2013,29(6):656-661.

    [9]Battistoni G,Cappucci F,Bertolino N,et al.Monte Carlo calculated and experimentally determined output correction factors for small field detectors in Leksell Gamma Knife Perfexion beams[J]. Phys Med Biol,2015,60(10):3959-3973.

    [10]E Pappas,A Moutsatsos,E Pantelis,et al.An EGSnrc-based model of the gamma knife perfexion unit for small field dosimetric calculation[J].Physica Medica,2014,30:e108-e109.

    [11]Kim S,Yoshizumi TT,Yin FF,et al.Spiral computed tomography phase-space sourcemodel in the BEAMnrc/EGSnrc Monte Carlo system:implementation and validation[J].Phys Med Biol,2013,58(8):2609-2624.

    [12]吳愛東,吳宜燦.CT圖像的體元大小對EGSnrc蒙特卡羅劑量計(jì)算的影響[J].核技術(shù),2007,30(2):143-146.

    [13]Novotny J Jr,Bhatnagar JP,Quader MA,et al. Measurement of relative output factors for the 8 and 4 mm collimators of Leksell Gamma Knife Perfexion by film dosimetry[J].Med Phys,2009,36(5):1768-1774.

    [14]Ralston A,Liu P,Warrener K,et al.Small field diode correction factors derived using an air core fibre optic scintillation dosimeter and EBT2 film[J].Phys Med Biol,2012,57(9):2587-2602.

    [15]Mack A,Scheib SG,Major J,et al.Precision dosimetry for narrow photon beams used in radiosurgery-Determination of Gamma Knife?output factors[J].Med Phys,2002,29(9):2080-2089.

    [16]Tsai JS,Rivard MJ,Engler MJ,et al.Determination of the 4 mm Gamma Knife helmet relative output factor using a variety of detectors[J].Med Phys,2003,30(5):986-992.

    [17]Sauer OA,Wilbert J.Measurement of output factors for small photon beams[J].Med Phys,2007,34(6):1983-1988.

    [18]Novotny J,Bhatnagar JP,Niranjan A,et al. Dosimetric comparison of the Leksell Gamma Knife Perfexion and 4C[J].J Neurosurg,2008,109:8-14.

    Monte Carlo Simulation study of LUNA stereotactic body radiotherapy device’s output factors

    YUAN Ji-long, CHENG Jin-sheng


    China Medical Equipment,2016,13(9):10-13.

    Objective: Using Monte Carlo simulation method to study the important dosimetric parameter output factor of the LUNA stereotactic body radiotherapy device. Methods: Following the manufacturer’s technical documentation to establish the LUNA stereotactic body radiotherapy device’s Monte Carlo model and after model validation. After referring to the output factor measurement’s detector provided by manufacturer, set the relevant parameters in Monte Carlo model at last, compare and analyze the difference of output factor which provided between by manufacturer and calculated by Monte Carlo simulation. Results: In the range of 3% deviation,under the different collimator of condition, output factors which provided by manufacturers and calculated by Monte Carlo are in good agreement, and the reasons for the difference are also analyzed. Conclusion: LUNA stereotactic body radiotherapy device’s output factors of Monte Carlo simulated research work provides reference and basis for the clinical use and quality control.

    LUNA stereotactic body radiotherapy device; Monte carlo simulation; Output factor

    1672-8270(2016)09-0010-04

    R197.39

    A

    10.3969/J.ISSN.1672-8270.2016.09.003

    2016-06-02

    國際原子能機(jī)構(gòu)2014-2015國際合作項(xiàng)目(CPR6006)“加強(qiáng)現(xiàn)代核醫(yī)學(xué)、放射診斷和放射治療實(shí)踐的質(zhì)量保證”

    ①中國疾病預(yù)防控制中心輻射防護(hù)與核安全醫(yī)學(xué)所 輻射防護(hù)與核應(yīng)急中國疾病預(yù)防控制中心重點(diǎn)實(shí)驗(yàn)室 北京 100088

    chengjs3393@163.com

    袁繼龍,男,(1982- ),博士,助理研究員。中國疾病預(yù)防控制中心輻射防護(hù)與核安全醫(yī)學(xué)所輻射防護(hù)與核應(yīng)急中國疾病預(yù)防控制中心重點(diǎn)實(shí)驗(yàn)室,研究方向:輻射劑量與輻射防護(hù)。

    猜你喜歡
    伽馬刀蒙特卡羅模擬計(jì)算
    R1234ze PVTx熱物性模擬計(jì)算
    能源工程(2022年1期)2022-03-29 01:06:26
    復(fù)方丹參注射液聯(lián)合伽馬刀治療對原發(fā)性肝癌患者免疫功能的影響
    利用蒙特卡羅方法求解二重積分
    TACE與伽瑪?shù)吨委煾伟┎涣挤磻?yīng)發(fā)生情況、臨床治療效果分析
    伽馬刀治療多發(fā)腦轉(zhuǎn)移瘤的療效分析
    擠出發(fā)泡片材褶皺分析及模擬計(jì)算
    中國塑料(2016年9期)2016-06-13 03:18:54
    伽馬刀治療腦動(dòng)靜脈畸形的遠(yuǎn)期療效分析
    探討蒙特卡羅方法在解微分方程邊值問題中的應(yīng)用
    復(fù)合型種子源125I-103Pd劑量場分布的蒙特卡羅模擬與實(shí)驗(yàn)測定
    同位素(2014年2期)2014-04-16 04:57:20
    實(shí)際發(fā)射工況下底排藥柱結(jié)構(gòu)完整性的模擬計(jì)算
    亚洲精华国产精华液的使用体验| 午夜福利在线在线| 日韩制服骚丝袜av| 亚洲色图综合在线观看| 蜜桃亚洲精品一区二区三区| kizo精华| 亚洲欧美精品专区久久| 搡老乐熟女国产| 五月天丁香电影| 高清在线视频一区二区三区| 国产熟女欧美一区二区| 最新中文字幕久久久久| 国产黄片视频在线免费观看| 久久99热6这里只有精品| 成人午夜精彩视频在线观看| 91午夜精品亚洲一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 人妻少妇偷人精品九色| 中文字幕人妻熟人妻熟丝袜美| 精品少妇黑人巨大在线播放| 久久热精品热| 九色成人免费人妻av| videossex国产| 三级国产精品片| 国产成人精品一,二区| 免费看光身美女| 亚洲精品日韩在线中文字幕| 在线观看三级黄色| 男人和女人高潮做爰伦理| 熟妇人妻不卡中文字幕| 看十八女毛片水多多多| 日本av免费视频播放| 在线观看三级黄色| av女优亚洲男人天堂| 日日啪夜夜爽| 美女脱内裤让男人舔精品视频| 日韩中字成人| 特大巨黑吊av在线直播| 国产精品一区二区在线观看99| 深夜a级毛片| 欧美另类一区| 欧美精品亚洲一区二区| 国产淫语在线视频| 国产乱人视频| 丰满乱子伦码专区| a 毛片基地| 日本欧美视频一区| 欧美日韩一区二区视频在线观看视频在线| 亚洲伊人久久精品综合| 在线观看美女被高潮喷水网站| 亚洲四区av| 亚洲av在线观看美女高潮| 成人高潮视频无遮挡免费网站| 欧美xxxx性猛交bbbb| 久久精品国产a三级三级三级| 插阴视频在线观看视频| 国产一区二区三区av在线| 26uuu在线亚洲综合色| 我的女老师完整版在线观看| 国产成人精品婷婷| 成人一区二区视频在线观看| 欧美变态另类bdsm刘玥| 欧美激情极品国产一区二区三区 | 国产极品天堂在线| 色婷婷av一区二区三区视频| 91久久精品国产一区二区三区| 春色校园在线视频观看| av视频免费观看在线观看| 男人舔奶头视频| 黑人猛操日本美女一级片| 国产 精品1| 80岁老熟妇乱子伦牲交| 日本色播在线视频| av卡一久久| av线在线观看网站| 国产女主播在线喷水免费视频网站| 国产成人精品婷婷| 国产无遮挡羞羞视频在线观看| 人妻一区二区av| 国产一区二区三区av在线| av在线播放精品| 极品教师在线视频| 午夜视频国产福利| 91精品一卡2卡3卡4卡| 日韩强制内射视频| 亚洲中文av在线| 国产在线视频一区二区| 交换朋友夫妻互换小说| 亚洲欧美日韩无卡精品| 国产成人精品婷婷| 毛片一级片免费看久久久久| 日韩 亚洲 欧美在线| 久久精品国产亚洲av天美| 亚洲精华国产精华液的使用体验| 日本免费在线观看一区| 少妇人妻一区二区三区视频| 一本色道久久久久久精品综合| 色视频在线一区二区三区| 久久午夜福利片| 99热全是精品| 免费人成在线观看视频色| 亚洲色图av天堂| 三级国产精品欧美在线观看| 免费人妻精品一区二区三区视频| 亚洲国产av新网站| 久久ye,这里只有精品| 午夜福利影视在线免费观看| 一级av片app| 99久久精品一区二区三区| 波野结衣二区三区在线| 中文欧美无线码| 国产成人精品婷婷| 亚洲高清免费不卡视频| 麻豆精品久久久久久蜜桃| 国产高清三级在线| 看免费成人av毛片| 免费高清在线观看视频在线观看| 久久久久久久精品精品| tube8黄色片| 成年美女黄网站色视频大全免费 | 亚洲av国产av综合av卡| 精品久久久久久久末码| 国产精品麻豆人妻色哟哟久久| 18+在线观看网站| 亚洲av欧美aⅴ国产| 亚洲欧美成人综合另类久久久| 午夜免费鲁丝| 免费黄频网站在线观看国产| 街头女战士在线观看网站| 久久青草综合色| 99久久人妻综合| 永久网站在线| 日韩强制内射视频| 精品亚洲成国产av| 亚洲精品自拍成人| 天美传媒精品一区二区| 免费少妇av软件| 亚洲精品视频女| 国产综合精华液| 熟妇人妻不卡中文字幕| 两个人的视频大全免费| 久久6这里有精品| 老师上课跳d突然被开到最大视频| 在线看a的网站| 国产高清国产精品国产三级 | 亚洲精品乱久久久久久| 韩国av在线不卡| 人妻制服诱惑在线中文字幕| 又粗又硬又长又爽又黄的视频| 亚洲精品成人av观看孕妇| 国产在线视频一区二区| 亚洲内射少妇av| 超碰97精品在线观看| 波野结衣二区三区在线| 青春草亚洲视频在线观看| 中文字幕亚洲精品专区| 日韩av不卡免费在线播放| 免费观看性生交大片5| av在线app专区| 国产乱人视频| 亚洲色图综合在线观看| 国产欧美日韩精品一区二区| 国产精品av视频在线免费观看| 我要看黄色一级片免费的| 18+在线观看网站| 欧美bdsm另类| 国产亚洲精品久久久com| 日韩在线高清观看一区二区三区| 国产成人免费观看mmmm| 久久99热6这里只有精品| .国产精品久久| 亚洲自偷自拍三级| 九草在线视频观看| 男的添女的下面高潮视频| 亚洲自偷自拍三级| 丰满少妇做爰视频| 午夜精品国产一区二区电影| 91精品一卡2卡3卡4卡| 国产永久视频网站| 五月玫瑰六月丁香| 国产免费一区二区三区四区乱码| 亚洲成人一二三区av| 97超碰精品成人国产| 日本黄大片高清| 七月丁香在线播放| 国模一区二区三区四区视频| 纯流量卡能插随身wifi吗| 在线精品无人区一区二区三 | 久久韩国三级中文字幕| 国产v大片淫在线免费观看| 亚洲图色成人| 色综合色国产| 精品久久久久久久久亚洲| 97在线视频观看| 国产精品久久久久久久久免| 一个人看的www免费观看视频| 欧美bdsm另类| 国产高清有码在线观看视频| 国产精品久久久久久久久免| 国产精品一区二区三区四区免费观看| 久久精品久久久久久噜噜老黄| 免费观看的影片在线观看| 精品久久久噜噜| 这个男人来自地球电影免费观看 | 成人一区二区视频在线观看| 在线看a的网站| a级毛片免费高清观看在线播放| 一级av片app| 插逼视频在线观看| 色哟哟·www| 一边亲一边摸免费视频| 青春草国产在线视频| 精品酒店卫生间| 爱豆传媒免费全集在线观看| 在线免费观看不下载黄p国产| 十八禁网站网址无遮挡 | 久久久精品免费免费高清| 亚洲一区二区三区欧美精品| 激情 狠狠 欧美| 免费在线观看成人毛片| 尾随美女入室| 日韩中字成人| 我要看黄色一级片免费的| 欧美日本视频| 免费观看a级毛片全部| 女人十人毛片免费观看3o分钟| 丝瓜视频免费看黄片| 美女脱内裤让男人舔精品视频| 2018国产大陆天天弄谢| 欧美3d第一页| 黄色一级大片看看| 精品国产三级普通话版| 男人和女人高潮做爰伦理| 赤兔流量卡办理| 日本vs欧美在线观看视频 | 伊人久久精品亚洲午夜| 国产精品免费大片| 直男gayav资源| 这个男人来自地球电影免费观看 | 18禁在线无遮挡免费观看视频| 国产成人午夜福利电影在线观看| 日韩av免费高清视频| 女的被弄到高潮叫床怎么办| 日韩av在线免费看完整版不卡| 亚洲婷婷狠狠爱综合网| 欧美日本视频| 欧美少妇被猛烈插入视频| 中文字幕人妻熟人妻熟丝袜美| 中文欧美无线码| 国产片特级美女逼逼视频| 国产精品一区二区三区四区免费观看| 欧美国产精品一级二级三级 | 肉色欧美久久久久久久蜜桃| 精品久久久久久久末码| 人人妻人人爽人人添夜夜欢视频 | 久久久久久久精品精品| 中文资源天堂在线| 一级毛片电影观看| 免费人成在线观看视频色| 亚洲伊人久久精品综合| 99热这里只有精品一区| 2022亚洲国产成人精品| 黄色视频在线播放观看不卡| 制服丝袜香蕉在线| 中文字幕精品免费在线观看视频 | 久久久久网色| 一级毛片黄色毛片免费观看视频| 欧美xxxx黑人xx丫x性爽| 欧美老熟妇乱子伦牲交| 夫妻午夜视频| av专区在线播放| 日日撸夜夜添| 老熟女久久久| 搡女人真爽免费视频火全软件| 久久精品国产亚洲网站| 22中文网久久字幕| 日韩亚洲欧美综合| 午夜福利在线观看免费完整高清在| 国产av国产精品国产| 亚洲人成网站在线观看播放| 一级a做视频免费观看| 又爽又黄a免费视频| 亚洲熟女精品中文字幕| 美女高潮的动态| 午夜福利网站1000一区二区三区| 成人无遮挡网站| 九草在线视频观看| av国产精品久久久久影院| 国产成人午夜福利电影在线观看| 免费av中文字幕在线| av在线老鸭窝| 成人国产av品久久久| 我要看黄色一级片免费的| 亚洲色图av天堂| 色视频在线一区二区三区| 成人午夜精彩视频在线观看| 亚洲欧洲日产国产| 久久久久久久精品精品| 亚洲欧美日韩另类电影网站 | 久久热精品热| 高清av免费在线| 午夜激情福利司机影院| 国产精品人妻久久久久久| 精华霜和精华液先用哪个| 国产免费一区二区三区四区乱码| 免费少妇av软件| 国产黄色免费在线视频| 久久97久久精品| 中文在线观看免费www的网站| 亚洲av欧美aⅴ国产| 色视频在线一区二区三区| 婷婷色av中文字幕| av天堂中文字幕网| 久久久国产一区二区| 日韩一区二区三区影片| 亚洲国产日韩一区二区| 欧美性感艳星| 人妻一区二区av| 精品一区在线观看国产| 国产成人午夜福利电影在线观看| 久久ye,这里只有精品| 国产欧美日韩精品一区二区| 男女边吃奶边做爰视频| 久热这里只有精品99| 亚洲精品乱久久久久久| 丝袜喷水一区| 各种免费的搞黄视频| 直男gayav资源| 最新中文字幕久久久久| 丝瓜视频免费看黄片| 黄片wwwwww| 妹子高潮喷水视频| 街头女战士在线观看网站| 99热全是精品| 亚洲欧美精品专区久久| 亚洲av综合色区一区| 日日啪夜夜撸| 国产视频内射| 国产精品福利在线免费观看| 免费看日本二区| 欧美精品一区二区免费开放| 香蕉精品网在线| 亚洲美女黄色视频免费看| a级毛色黄片| 边亲边吃奶的免费视频| 人妻少妇偷人精品九色| 亚洲精品视频女| 51国产日韩欧美| 91久久精品电影网| 国产精品99久久99久久久不卡 | 精品国产一区二区三区久久久樱花 | 韩国av在线不卡| 亚洲精品亚洲一区二区| 青春草国产在线视频| 日产精品乱码卡一卡2卡三| 一级黄片播放器| 亚洲精品亚洲一区二区| 性色av一级| 国产黄频视频在线观看| 亚洲精品456在线播放app| 亚洲熟女精品中文字幕| 狠狠精品人妻久久久久久综合| 最黄视频免费看| 极品教师在线视频| 亚洲欧美成人综合另类久久久| 99热这里只有精品一区| 99热网站在线观看| 精品久久久噜噜| 在线观看三级黄色| 人妻 亚洲 视频| 国产男女内射视频| 蜜桃在线观看..| 亚洲av电影在线观看一区二区三区| 赤兔流量卡办理| 日韩强制内射视频| 99re6热这里在线精品视频| 综合色丁香网| 国产精品一区二区三区四区免费观看| 亚洲高清免费不卡视频| 亚洲电影在线观看av| 高清欧美精品videossex| 天堂中文最新版在线下载| 久久久色成人| 亚洲精品乱码久久久久久按摩| 亚洲精品第二区| 亚洲一区二区三区欧美精品| 婷婷色av中文字幕| 亚洲美女视频黄频| 亚洲精品一区蜜桃| 国产欧美另类精品又又久久亚洲欧美| 男人添女人高潮全过程视频| 成人特级av手机在线观看| 国产又色又爽无遮挡免| 视频中文字幕在线观看| 国产白丝娇喘喷水9色精品| av卡一久久| 美女脱内裤让男人舔精品视频| 春色校园在线视频观看| 欧美精品人与动牲交sv欧美| 麻豆乱淫一区二区| 3wmmmm亚洲av在线观看| 精品亚洲成a人片在线观看 | 不卡视频在线观看欧美| 亚洲婷婷狠狠爱综合网| 大话2 男鬼变身卡| 狂野欧美白嫩少妇大欣赏| 国产亚洲5aaaaa淫片| 国产69精品久久久久777片| 欧美精品一区二区免费开放| 高清毛片免费看| 国产精品一及| 午夜免费观看性视频| 国产亚洲最大av| 国产成人一区二区在线| 男男h啪啪无遮挡| av天堂中文字幕网| 在线 av 中文字幕| 日日啪夜夜撸| 欧美极品一区二区三区四区| 美女国产视频在线观看| 成年女人在线观看亚洲视频| 18禁在线播放成人免费| 成人亚洲精品一区在线观看 | 岛国毛片在线播放| 亚洲精品国产av成人精品| 亚洲欧美一区二区三区黑人 | 日韩成人伦理影院| 久久 成人 亚洲| 亚洲成人av在线免费| 亚洲欧美清纯卡通| 亚洲激情五月婷婷啪啪| 欧美日韩亚洲高清精品| 国产精品国产三级国产av玫瑰| 黄色欧美视频在线观看| 国产精品精品国产色婷婷| 国产欧美亚洲国产| 亚洲综合精品二区| 亚洲精品乱码久久久久久按摩| 啦啦啦啦在线视频资源| 亚洲伊人久久精品综合| 一级毛片aaaaaa免费看小| 青春草亚洲视频在线观看| 看免费成人av毛片| 免费观看的影片在线观看| 日韩人妻高清精品专区| 一级av片app| 99久久精品热视频| 美女中出高潮动态图| 我要看日韩黄色一级片| 黄色配什么色好看| 春色校园在线视频观看| freevideosex欧美| 熟妇人妻不卡中文字幕| 一级片'在线观看视频| 国产精品人妻久久久久久| 亚洲精品久久午夜乱码| 日本免费在线观看一区| 秋霞在线观看毛片| 亚洲av成人精品一二三区| 特大巨黑吊av在线直播| 伦精品一区二区三区| 五月天丁香电影| 黄色怎么调成土黄色| 国产成人精品一,二区| 又爽又黄a免费视频| 不卡视频在线观看欧美| 久久久午夜欧美精品| av免费观看日本| 少妇裸体淫交视频免费看高清| 国产视频首页在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲成人一二三区av| 中文字幕久久专区| 国产高清不卡午夜福利| 日日啪夜夜撸| 亚洲内射少妇av| 男人和女人高潮做爰伦理| 久久99热这里只频精品6学生| 91在线精品国自产拍蜜月| 欧美xxxx性猛交bbbb| 国产欧美日韩精品一区二区| 国产男女超爽视频在线观看| 欧美日韩在线观看h| 欧美一区二区亚洲| 妹子高潮喷水视频| 久久这里有精品视频免费| 一级毛片黄色毛片免费观看视频| 91aial.com中文字幕在线观看| 日本猛色少妇xxxxx猛交久久| 久久97久久精品| a 毛片基地| 内射极品少妇av片p| 免费观看a级毛片全部| 校园人妻丝袜中文字幕| 青春草国产在线视频| 赤兔流量卡办理| 纵有疾风起免费观看全集完整版| 亚洲av.av天堂| 亚洲aⅴ乱码一区二区在线播放| 午夜精品国产一区二区电影| 国产在线男女| 久久99蜜桃精品久久| 国语对白做爰xxxⅹ性视频网站| 国产精品久久久久久精品古装| 夜夜爽夜夜爽视频| 身体一侧抽搐| 亚洲国产欧美在线一区| 99久久综合免费| 免费黄网站久久成人精品| 国产精品一区二区三区四区免费观看| 自拍偷自拍亚洲精品老妇| 高清av免费在线| 国产精品嫩草影院av在线观看| 丝袜喷水一区| 久久精品久久久久久噜噜老黄| 一级黄片播放器| 亚洲一区二区三区欧美精品| 熟女人妻精品中文字幕| 人体艺术视频欧美日本| 国产高清有码在线观看视频| 欧美高清成人免费视频www| 久久人人爽人人爽人人片va| 99久久人妻综合| 免费观看a级毛片全部| 人妻系列 视频| 我的老师免费观看完整版| 毛片一级片免费看久久久久| 91狼人影院| 成人特级av手机在线观看| 免费黄网站久久成人精品| 国产精品偷伦视频观看了| 亚洲欧美中文字幕日韩二区| h日本视频在线播放| 综合色丁香网| 久久久欧美国产精品| xxx大片免费视频| 涩涩av久久男人的天堂| 又大又黄又爽视频免费| 91精品国产九色| 在线观看免费日韩欧美大片 | 全区人妻精品视频| 国产av国产精品国产| 人人妻人人爽人人添夜夜欢视频 | 日产精品乱码卡一卡2卡三| 在现免费观看毛片| 99精国产麻豆久久婷婷| 久久国产亚洲av麻豆专区| www.av在线官网国产| 99久久综合免费| 在线天堂最新版资源| 国产精品久久久久成人av| 亚洲av日韩在线播放| 国产美女午夜福利| 亚洲av国产av综合av卡| 日韩 亚洲 欧美在线| 男的添女的下面高潮视频| 色婷婷久久久亚洲欧美| 熟女电影av网| 性色avwww在线观看| 波野结衣二区三区在线| 日日摸夜夜添夜夜爱| 伊人久久精品亚洲午夜| 久久综合国产亚洲精品| 国产精品国产av在线观看| 大片电影免费在线观看免费| 丝袜喷水一区| 色哟哟·www| 国产大屁股一区二区在线视频| 精品国产三级普通话版| 在线亚洲精品国产二区图片欧美 | 欧美日韩亚洲高清精品| 精品人妻视频免费看| 热re99久久精品国产66热6| 国产 一区精品| 建设人人有责人人尽责人人享有的 | 亚洲人成网站高清观看| 精品午夜福利在线看| 久久久a久久爽久久v久久| 各种免费的搞黄视频| tube8黄色片| 日韩三级伦理在线观看| 成人免费观看视频高清| 免费观看无遮挡的男女| 伦理电影大哥的女人| 性色avwww在线观看| 身体一侧抽搐| 青春草国产在线视频| 最新中文字幕久久久久| 乱码一卡2卡4卡精品| 免费久久久久久久精品成人欧美视频 | 国产综合精华液| 精品少妇黑人巨大在线播放| 人人妻人人添人人爽欧美一区卜 | 免费高清在线观看视频在线观看| 精品人妻熟女av久视频| 日韩一本色道免费dvd| 久久久久性生活片| 蜜臀久久99精品久久宅男| a级毛色黄片| 国产精品成人在线| 欧美日韩亚洲高清精品| 亚洲综合色惰| 大又大粗又爽又黄少妇毛片口| 高清av免费在线| 小蜜桃在线观看免费完整版高清| 美女xxoo啪啪120秒动态图| 亚洲性久久影院| 精品久久久久久久末码| 午夜福利网站1000一区二区三区| 狂野欧美白嫩少妇大欣赏| 一区在线观看完整版| 亚洲人成网站在线播| 亚洲国产av新网站| 欧美性感艳星| 不卡视频在线观看欧美|