• <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ì)算
    丁香六月天网| 大片免费播放器 马上看| 免费在线观看完整版高清| 亚洲精品中文字幕在线视频| 50天的宝宝边吃奶边哭怎么回事| a在线观看视频网站| 91成年电影在线观看| 亚洲国产看品久久| 国产精品久久久久久精品古装| 女人爽到高潮嗷嗷叫在线视频| 最新美女视频免费是黄的| 国产麻豆69| 99国产精品一区二区蜜桃av | 精品少妇久久久久久888优播| 老司机午夜十八禁免费视频| cao死你这个sao货| 另类亚洲欧美激情| 日本黄色视频三级网站网址 | 午夜福利在线免费观看网站| 国产精品一区二区精品视频观看| 精品少妇久久久久久888优播| 欧美中文综合在线视频| 亚洲天堂av无毛| 欧美激情极品国产一区二区三区| 国产xxxxx性猛交| 久久久久久人人人人人| 岛国在线观看网站| 91精品国产国语对白视频| 亚洲国产欧美在线一区| 最黄视频免费看| 久久精品91无色码中文字幕| 国产精品av久久久久免费| 在线永久观看黄色视频| 每晚都被弄得嗷嗷叫到高潮| 日本黄色视频三级网站网址 | 午夜91福利影院| 亚洲第一欧美日韩一区二区三区 | 日韩欧美国产一区二区入口| 国产成+人综合+亚洲专区| 精品欧美一区二区三区在线| 亚洲九九香蕉| 蜜桃国产av成人99| 免费一级毛片在线播放高清视频 | 在线观看免费午夜福利视频| 日韩视频一区二区在线观看| av有码第一页| 久久久久精品人妻al黑| 国产精品久久久av美女十八| 久久久久久久大尺度免费视频| 欧美日韩黄片免| 欧美激情久久久久久爽电影 | 精品午夜福利视频在线观看一区 | 免费久久久久久久精品成人欧美视频| 如日韩欧美国产精品一区二区三区| 亚洲黑人精品在线| 一夜夜www| 在线天堂中文资源库| 啦啦啦免费观看视频1| 麻豆乱淫一区二区| 王馨瑶露胸无遮挡在线观看| 啦啦啦中文免费视频观看日本| 国产成人免费观看mmmm| 日本一区二区免费在线视频| tocl精华| 又大又爽又粗| 一边摸一边抽搐一进一小说 | 99国产精品免费福利视频| 成年动漫av网址| 一个人免费看片子| 97在线人人人人妻| 无遮挡黄片免费观看| 老熟女久久久| 午夜精品国产一区二区电影| 久久精品国产a三级三级三级| 国产欧美日韩一区二区三区在线| 建设人人有责人人尽责人人享有的| 国产精品熟女久久久久浪| 精品一区二区三区av网在线观看 | 久久精品亚洲熟妇少妇任你| 超色免费av| 免费看a级黄色片| 精品福利永久在线观看| 纯流量卡能插随身wifi吗| 欧美精品人与动牲交sv欧美| 美女午夜性视频免费| 亚洲国产欧美网| 最新在线观看一区二区三区| 一个人免费在线观看的高清视频| 日韩欧美一区二区三区在线观看 | 色精品久久人妻99蜜桃| 美女高潮到喷水免费观看| 亚洲精品久久成人aⅴ小说| 丝袜在线中文字幕| 亚洲专区字幕在线| 久久青草综合色| 18禁裸乳无遮挡动漫免费视频| 色精品久久人妻99蜜桃| 欧美精品人与动牲交sv欧美| www.999成人在线观看| 亚洲国产毛片av蜜桃av| www.自偷自拍.com| 亚洲精品国产区一区二| 变态另类成人亚洲欧美熟女 | www.精华液| 亚洲国产av影院在线观看| 又大又爽又粗| 一夜夜www| 久久国产精品大桥未久av| 欧美黑人欧美精品刺激| 亚洲欧美色中文字幕在线| 婷婷成人精品国产| 久热爱精品视频在线9| 又黄又粗又硬又大视频| 自拍欧美九色日韩亚洲蝌蚪91| 香蕉国产在线看| 成人黄色视频免费在线看| 国产欧美日韩一区二区三| 日韩欧美一区二区三区在线观看 | 欧美日韩成人在线一区二区| 男女下面插进去视频免费观看| 宅男免费午夜| 午夜福利一区二区在线看| 日韩欧美国产一区二区入口| 另类精品久久| 嫩草影视91久久| 12—13女人毛片做爰片一| 美女福利国产在线| 老司机福利观看| 在线观看免费午夜福利视频| 啪啪无遮挡十八禁网站| 国产亚洲欧美在线一区二区| 午夜福利乱码中文字幕| 男人舔女人的私密视频| 9热在线视频观看99| 青青草视频在线视频观看| 99久久精品国产亚洲精品| 麻豆av在线久日| 99精国产麻豆久久婷婷| 国产主播在线观看一区二区| 嫩草影视91久久| 黄片播放在线免费| 久久久国产成人免费| 国产成人精品久久二区二区91| 啦啦啦 在线观看视频| 日韩视频在线欧美| 一本大道久久a久久精品| 国产男女内射视频| 亚洲国产欧美在线一区| 午夜福利在线观看吧| 久久久精品免费免费高清| 色视频在线一区二区三区| 高潮久久久久久久久久久不卡| 久久ye,这里只有精品| 国产片内射在线| av视频免费观看在线观看| 一级毛片精品| 十分钟在线观看高清视频www| 欧美黄色淫秽网站| 国产成人欧美在线观看 | 男女床上黄色一级片免费看| 精品少妇黑人巨大在线播放| a级片在线免费高清观看视频| 国产不卡一卡二| 欧美av亚洲av综合av国产av| 国产深夜福利视频在线观看| 久久久久久久国产电影| 天天躁夜夜躁狠狠躁躁| 好男人电影高清在线观看| 久久婷婷成人综合色麻豆| 一边摸一边做爽爽视频免费| 国产精品免费视频内射| 一区二区三区激情视频| 一级毛片电影观看| 欧美日韩一级在线毛片| 极品人妻少妇av视频| 日本wwww免费看| 十八禁人妻一区二区| 亚洲av国产av综合av卡| 久久久久久久精品吃奶| a在线观看视频网站| 久久精品国产亚洲av高清一级| 狠狠狠狠99中文字幕| 国产精品偷伦视频观看了| 老汉色av国产亚洲站长工具| 欧美激情 高清一区二区三区| 人妻 亚洲 视频| 国产淫语在线视频| 国产高清激情床上av| 亚洲,欧美精品.| www.999成人在线观看| 十八禁网站免费在线| 国产在线视频一区二区| 久久精品aⅴ一区二区三区四区| 国产深夜福利视频在线观看| 中文字幕制服av| 精品亚洲乱码少妇综合久久| 伊人久久大香线蕉亚洲五| av天堂在线播放| 精品亚洲乱码少妇综合久久| 色在线成人网| 考比视频在线观看| 美女视频免费永久观看网站| www.自偷自拍.com| 日本欧美视频一区| 啦啦啦免费观看视频1| 国产主播在线观看一区二区| 午夜福利,免费看| 岛国毛片在线播放| 久久国产精品大桥未久av| 黄频高清免费视频| 纵有疾风起免费观看全集完整版| 久久免费观看电影| 99久久精品国产亚洲精品| 自拍欧美九色日韩亚洲蝌蚪91| 正在播放国产对白刺激| 一级毛片女人18水好多| 午夜福利视频精品| 亚洲色图综合在线观看| 男男h啪啪无遮挡| 一区二区三区乱码不卡18| 亚洲 欧美一区二区三区| 国精品久久久久久国模美| 日日爽夜夜爽网站| 欧美人与性动交α欧美软件| 久久青草综合色| xxxhd国产人妻xxx| h视频一区二区三区| 中文字幕制服av| 自线自在国产av| 国产亚洲精品第一综合不卡| 久久久精品免费免费高清| 亚洲av美国av| 高清视频免费观看一区二区| 99re6热这里在线精品视频| 国产亚洲欧美在线一区二区| 99国产综合亚洲精品| 国产欧美日韩一区二区三| 久久亚洲精品不卡| 99re在线观看精品视频| av欧美777| 美女高潮到喷水免费观看| 少妇裸体淫交视频免费看高清 | 脱女人内裤的视频| 午夜91福利影院| 9热在线视频观看99| 女人被躁到高潮嗷嗷叫费观| 美女视频免费永久观看网站| 亚洲国产成人一精品久久久| 日本欧美视频一区| 国产免费av片在线观看野外av| 欧美老熟妇乱子伦牲交| 亚洲av电影在线进入| 亚洲天堂av无毛| 亚洲人成伊人成综合网2020| 亚洲avbb在线观看| 久久99热这里只频精品6学生| 一本综合久久免费| 天堂俺去俺来也www色官网| 韩国精品一区二区三区| 一个人免费看片子| 日韩欧美一区二区三区在线观看 | 欧美激情久久久久久爽电影 | 国产精品久久久久久精品电影小说| 12—13女人毛片做爰片一| 亚洲国产毛片av蜜桃av| 日韩欧美三级三区| 18禁美女被吸乳视频| 电影成人av| 男女边摸边吃奶| 首页视频小说图片口味搜索| 亚洲综合色网址| 50天的宝宝边吃奶边哭怎么回事| 欧美日韩国产mv在线观看视频| 女性被躁到高潮视频| 黄色 视频免费看| 亚洲精品久久成人aⅴ小说| 国产在线免费精品| 最新的欧美精品一区二区| 日本wwww免费看| 777米奇影视久久| www.熟女人妻精品国产| 日日爽夜夜爽网站| 亚洲午夜精品一区,二区,三区| 999精品在线视频| 欧美在线黄色| 日韩熟女老妇一区二区性免费视频| 国产精品 国内视频| 免费不卡黄色视频| 黄色丝袜av网址大全| 99精品在免费线老司机午夜| www日本在线高清视频| 侵犯人妻中文字幕一二三四区| 涩涩av久久男人的天堂| 岛国毛片在线播放| 欧美日韩成人在线一区二区| 高清av免费在线| 欧美另类亚洲清纯唯美| 人人妻人人添人人爽欧美一区卜| 久久人妻av系列| 在线十欧美十亚洲十日本专区| 亚洲专区中文字幕在线| 蜜桃在线观看..| 91精品三级在线观看| 亚洲一码二码三码区别大吗| 视频区图区小说| 日日夜夜操网爽| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美日韩另类电影网站| 18禁观看日本| 久热这里只有精品99| 国产男女超爽视频在线观看| 桃红色精品国产亚洲av| 两人在一起打扑克的视频| 人人澡人人妻人| 超碰成人久久| 亚洲 国产 在线| 宅男免费午夜| 香蕉国产在线看| 一级毛片电影观看| 精品卡一卡二卡四卡免费| 飞空精品影院首页| 一区福利在线观看| 国产精品免费一区二区三区在线 | 国精品久久久久久国模美| 一级黄色大片毛片| 久久婷婷成人综合色麻豆| 精品亚洲乱码少妇综合久久| 我要看黄色一级片免费的| 精品少妇一区二区三区视频日本电影| 嫩草影视91久久| 91精品三级在线观看| a在线观看视频网站| 老司机福利观看| 久久 成人 亚洲| 丝袜在线中文字幕| 性高湖久久久久久久久免费观看| av福利片在线| 精品久久久久久电影网| 日本av手机在线免费观看| 一区二区三区激情视频| 国产不卡一卡二| 免费女性裸体啪啪无遮挡网站| 亚洲色图 男人天堂 中文字幕| 丰满人妻熟妇乱又伦精品不卡| 国产视频一区二区在线看| 久久中文看片网| 69精品国产乱码久久久| 午夜91福利影院| 精品乱码久久久久久99久播| 三级毛片av免费| 在线 av 中文字幕| 日本精品一区二区三区蜜桃| 成年人黄色毛片网站| 日日摸夜夜添夜夜添小说| 国产成人系列免费观看| cao死你这个sao货| 黄色丝袜av网址大全| 欧美乱妇无乱码| 丰满迷人的少妇在线观看| 国产精品免费视频内射| 久热这里只有精品99| 精品一区二区三区四区五区乱码| 97人妻天天添夜夜摸| 天堂动漫精品| 中文字幕人妻熟女乱码| 国产在线精品亚洲第一网站| 不卡一级毛片| 黄色视频不卡| 国产精品秋霞免费鲁丝片| 真人做人爱边吃奶动态| 中文字幕制服av| 一个人免费在线观看的高清视频| 热99re8久久精品国产| 久久热在线av| 一边摸一边抽搐一进一出视频| 日韩制服丝袜自拍偷拍| 国产精品国产av在线观看| 久久亚洲真实| 国产99久久九九免费精品| 国产亚洲精品第一综合不卡| 一夜夜www| 黄色丝袜av网址大全| 在线观看www视频免费| 另类亚洲欧美激情| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美一区二区三区久久| 免费看a级黄色片| 国产成人精品无人区| 丰满人妻熟妇乱又伦精品不卡| 建设人人有责人人尽责人人享有的| 国产精品二区激情视频| 制服诱惑二区| 精品高清国产在线一区| 好男人电影高清在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 精品福利永久在线观看| 午夜久久久在线观看| 一区二区三区激情视频| 啦啦啦 在线观看视频| 亚洲成a人片在线一区二区| 亚洲午夜精品一区,二区,三区| 精品视频人人做人人爽| 久久久久久久精品吃奶| 搡老乐熟女国产| 亚洲精品中文字幕在线视频| 久久av网站| 成人黄色视频免费在线看| 成人精品一区二区免费| 色在线成人网| 久久人妻熟女aⅴ| 国产精品一区二区在线不卡| 成人av一区二区三区在线看| 亚洲精品国产精品久久久不卡| 日本wwww免费看| 精品福利观看| 91精品三级在线观看| 亚洲人成伊人成综合网2020| 99国产精品一区二区三区| avwww免费| 久久久久国内视频| 日本av手机在线免费观看| 亚洲美女黄片视频| 最近最新免费中文字幕在线| 制服人妻中文乱码| 免费在线观看影片大全网站| 色94色欧美一区二区| 精品欧美一区二区三区在线| 亚洲精品一卡2卡三卡4卡5卡| 亚洲成av片中文字幕在线观看| 三级毛片av免费| 美女国产高潮福利片在线看| 欧美黑人欧美精品刺激| 日本撒尿小便嘘嘘汇集6| 变态另类成人亚洲欧美熟女 | 日韩欧美一区二区三区在线观看 | 国产真人三级小视频在线观看| 久久天躁狠狠躁夜夜2o2o| 国产日韩一区二区三区精品不卡| 好男人电影高清在线观看| 在线看a的网站| 怎么达到女性高潮| videos熟女内射| 免费看十八禁软件| 色综合婷婷激情| 两个人免费观看高清视频| 波多野结衣一区麻豆| 欧美乱码精品一区二区三区| 久久精品国产亚洲av香蕉五月 | 亚洲精品中文字幕一二三四区 | 久久国产精品男人的天堂亚洲| 十八禁人妻一区二区| 日日夜夜操网爽| 日韩视频在线欧美| 黄色片一级片一级黄色片| 亚洲av美国av| 国产av又大| 咕卡用的链子| 两个人免费观看高清视频| 久久午夜亚洲精品久久| 国产免费视频播放在线视频| 国产91精品成人一区二区三区 | 精品国产一区二区久久| 91精品国产国语对白视频| 国产不卡一卡二| 国产欧美日韩综合在线一区二区| 欧美日韩av久久| 欧美激情极品国产一区二区三区| www.精华液| av天堂在线播放| 久久精品亚洲av国产电影网| 不卡av一区二区三区| 精品少妇黑人巨大在线播放| 久久久久国产一级毛片高清牌| 欧美 日韩 精品 国产| 亚洲精华国产精华精| 免费一级毛片在线播放高清视频 | 亚洲国产毛片av蜜桃av| 久久久国产成人免费| 老司机在亚洲福利影院| 国产伦理片在线播放av一区| 久久久久久久精品吃奶| 国产成人免费无遮挡视频| 热99久久久久精品小说推荐| 成人国语在线视频| 悠悠久久av| 午夜福利在线观看吧| 91九色精品人成在线观看| 午夜福利乱码中文字幕| 欧美日韩一级在线毛片| a级毛片黄视频| 国产精品九九99| 成人影院久久| 国产在线观看jvid| 男人操女人黄网站| 国产精品免费视频内射| 久久中文字幕人妻熟女| 久久99热这里只频精品6学生| 国产真人三级小视频在线观看| 免费在线观看黄色视频的| 老熟妇乱子伦视频在线观看| 十八禁网站免费在线| 欧美国产精品一级二级三级| 18禁观看日本| 中国美女看黄片| 日韩免费高清中文字幕av| 1024香蕉在线观看| 国产精品亚洲av一区麻豆| 国产无遮挡羞羞视频在线观看| 黄色视频,在线免费观看| 色视频在线一区二区三区| 在线永久观看黄色视频| 女人精品久久久久毛片| 丝袜在线中文字幕| 欧美精品一区二区大全| 午夜日韩欧美国产| 黄色视频不卡| 黄色a级毛片大全视频| 国产精品香港三级国产av潘金莲| 操出白浆在线播放| 国产精品一区二区精品视频观看| 在线天堂中文资源库| 国产高清videossex| 免费黄频网站在线观看国产| 99热国产这里只有精品6| 久久久久久免费高清国产稀缺| 亚洲精品自拍成人| 免费日韩欧美在线观看| 国产午夜精品久久久久久| 欧美大码av| 国产成人系列免费观看| 欧美日韩一级在线毛片| 国产亚洲一区二区精品| 久久亚洲真实| 99国产精品一区二区三区| 国产一区二区激情短视频| 亚洲欧洲精品一区二区精品久久久| 乱人伦中国视频| 亚洲五月婷婷丁香| 一级a爱视频在线免费观看| 99精国产麻豆久久婷婷| 欧美大码av| 波多野结衣av一区二区av| 黄色片一级片一级黄色片| 午夜两性在线视频| 国产精品九九99| 欧美中文综合在线视频| 亚洲国产毛片av蜜桃av| 国产黄色免费在线视频| 在线av久久热| 757午夜福利合集在线观看| 精品人妻熟女毛片av久久网站| 777米奇影视久久| 操美女的视频在线观看| 国产精品熟女久久久久浪| videosex国产| √禁漫天堂资源中文www| 少妇精品久久久久久久| av超薄肉色丝袜交足视频| 99国产精品一区二区三区| 午夜福利一区二区在线看| 99热国产这里只有精品6| 国产亚洲精品一区二区www | 国产欧美日韩精品亚洲av| 法律面前人人平等表现在哪些方面| 久久婷婷成人综合色麻豆| 免费观看人在逋| www.精华液| 巨乳人妻的诱惑在线观看| 国产亚洲欧美在线一区二区| av片东京热男人的天堂| 在线十欧美十亚洲十日本专区| 天堂动漫精品| 中亚洲国语对白在线视频| √禁漫天堂资源中文www| 老司机靠b影院| 亚洲性夜色夜夜综合| 在线观看免费高清a一片| 国产主播在线观看一区二区| 熟女少妇亚洲综合色aaa.| 激情视频va一区二区三区| 亚洲色图 男人天堂 中文字幕| 一边摸一边抽搐一进一出视频| 狠狠狠狠99中文字幕| 丰满迷人的少妇在线观看| 国产免费av片在线观看野外av| 久久天躁狠狠躁夜夜2o2o| 国产精品一区二区精品视频观看| 人妻一区二区av| 久久ye,这里只有精品| 日韩制服丝袜自拍偷拍| 日韩有码中文字幕| 久热爱精品视频在线9| 国产欧美亚洲国产| 亚洲熟女精品中文字幕| 黄色毛片三级朝国网站| 丝袜喷水一区| 久久亚洲真实| 国产精品98久久久久久宅男小说| 免费女性裸体啪啪无遮挡网站| 我要看黄色一级片免费的| 久久人妻熟女aⅴ| 色婷婷久久久亚洲欧美| 国产亚洲精品久久久久5区| 欧美性长视频在线观看| 男女之事视频高清在线观看| 老司机亚洲免费影院| 狠狠狠狠99中文字幕| 精品福利观看| 丝袜人妻中文字幕| 日本黄色视频三级网站网址 | 国产午夜精品久久久久久| 男男h啪啪无遮挡| 国产成人精品在线电影| 欧美日韩国产mv在线观看视频| 视频区欧美日本亚洲|