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

    虛擬單點(diǎn)源構(gòu)建醫(yī)用直線加速器束流模型的探討

    2021-04-20 07:57:50李惠娟焦勝修陳仲本劉小偉
    核技術(shù) 2021年4期
    關(guān)鍵詞:射野點(diǎn)源束流

    李惠娟 焦勝修 陳仲本 劉小偉

    1(中山大學(xué)新華學(xué)院 廣州 510520)

    2(中山大學(xué)物理學(xué)院 廣州 510275)

    放射治療是癌癥治療的三大主要手段之一。醫(yī)用直線加速器是目前放射治療的主流設(shè)備,其基本原理是通過(guò)高能電子打靶產(chǎn)生X射線對(duì)病灶區(qū)進(jìn)行照射,從而達(dá)到抑制和控制腫瘤的目的。為了達(dá)到該目的,吸收劑量計(jì)算是放射治療計(jì)劃制定的出發(fā)點(diǎn),而醫(yī)用直線加速器的束流模型是完成劑量計(jì)算所必須的基礎(chǔ)之一。因此對(duì)醫(yī)用直線加速器束流模型的研究在放射劑量計(jì)算的臨床應(yīng)用具有重要的意義。

    醫(yī)用直線加速器X射線的束流描述主要是蒙特卡羅模擬得到的相空間集合或由此得到的特性[1?2]和基于測(cè)量的劑量分布進(jìn)行參數(shù)擬合的虛源模型[3?8]。采用蒙特卡羅方法,根據(jù)實(shí)際加速器的結(jié)構(gòu),可模擬得到機(jī)頭出射粒子的相空間集合。雖然利用出射粒子的相空間集合更準(zhǔn)確地描述了加速器的束流,但由于各加速器之間存在差異,需針對(duì)具體的加速器進(jìn)行模擬,這需要廠家提供具體的加速器機(jī)頭的幾何及材料參數(shù)進(jìn)行模擬。虛源模型是采用虛擬的多個(gè)源對(duì)機(jī)頭出射的粒子進(jìn)行描述,各源的權(quán)重通過(guò)對(duì)實(shí)測(cè)的劑量分布(如百分深度劑量(Percentage Depth Dose,PDD)、離軸劑量比(Off Axis dose Ratio,OAR)、總散射因子(SCP)等),進(jìn)行擬合得到。醫(yī)用加速器的X射線可分為原射線和散射線,散射線主要來(lái)源于均整器,在虛源模型中多用面源進(jìn)行描述。但面源在束流方向及強(qiáng)度處理上較點(diǎn)源復(fù)雜,而采用點(diǎn)源將大大簡(jiǎn)化吸收劑量的計(jì)算。由于點(diǎn)源在射野中心點(diǎn)產(chǎn)生的注量不隨射野尺寸而變化,因而其不能正確描述散射源。采用高斯函數(shù)進(jìn)行卷積可使點(diǎn)源在射野中心點(diǎn)產(chǎn)生的注量隨射野尺寸而變化,從而得到面源的特征[9]。而本文嘗試采用另一種更簡(jiǎn)單的方法,通過(guò)設(shè)置點(diǎn)源的強(qiáng)度與射野的面積相關(guān)來(lái)描述散射源。

    本文初步探討了一種源強(qiáng)與射野的面積相關(guān)的單點(diǎn)源來(lái)描述散射線的醫(yī)用直線加速器X射線束流的虛源模型,并通對(duì)過(guò)實(shí)際測(cè)量的劑量分布進(jìn)行參數(shù)擬合來(lái)檢驗(yàn)?zāi)P偷倪m用性。

    1 材料與方法

    1.1 加速器X射線的單點(diǎn)源束流模型

    加速器的結(jié)構(gòu)示意圖如圖1 所示,其主要部件有 X 射 線 靶(Target)、初 級(jí) 準(zhǔn) 直 器(Primary Collimator)、均整器(Flatting Filter)、監(jiān)測(cè)電離室(Monitor Chamber)、反射鏡(Mirror)、次級(jí)準(zhǔn)直器(Secondary Collimator)及射野擋塊(Jaws)、多葉準(zhǔn)直器(Multi-leaf Collimator,MLC)等。均整器的作用是對(duì)X 射線束進(jìn)行調(diào)制。經(jīng)過(guò)均整器調(diào)制后,靠近中心軸的光子數(shù)減少,這樣就使得整體的光子注量分布沿半徑方向上較為平坦,從而使得它在模體中產(chǎn)生的劑量分布變得平坦。

    圖1 醫(yī)用直線加速器模型示意圖Fig.1 Sketch map of linear accelerator model

    單點(diǎn)源束流建模中,點(diǎn)源坐標(biāo)原點(diǎn)設(shè)在靶的中心點(diǎn),Z軸的方向?yàn)閳D1 中箭頭的指向。注量由式(1)描述:

    式中:Φ(r,z)為位置(x,y,z)處的注量,r2=x2+y2;ΦI(r,z)為打靶處產(chǎn)生射線貢獻(xiàn)份額,本文中稱為原射線;Φs(r,z)為均整器以及準(zhǔn)直器等部件產(chǎn)生散射線,該部分與射野有關(guān)的注量分布,本文稱為散射線。

    1.1.1 原射線建模

    由于電子打靶產(chǎn)生的韌致輻射光子發(fā)射不是各向同性的,同時(shí)考慮到均整器的調(diào)制作用,將原射線建模為各向異性的點(diǎn)源,射野內(nèi)射線的注量分布由式(2)描述:

    式中:Φ0可看成一種歸一化常數(shù);ɑ1、ɑ2、ɑ3為待擬合的參數(shù)。

    1.1.2 散射線建模

    散射線建模為強(qiáng)度與射野有關(guān)的點(diǎn)源。射野內(nèi)的散射線注量分布由式(3)描述:

    式中:F代表射野的面積;F0為參考射野的面積,實(shí)際處理時(shí)采用邊長(zhǎng)10 cm的方野的面積。

    1.1.3 射野外的注量

    由于本文主要關(guān)注射野內(nèi)的劑量,對(duì)于射野外的注量,僅采用簡(jiǎn)單的指數(shù)下降函數(shù)來(lái)描述:

    式中:Φout(d)代表射野外距射野邊緣的距離為d處的注量;ΦB代表射野邊緣處的注量;α為衰減系數(shù),通過(guò)調(diào)整其大小來(lái)符合野外的劑量分布。

    1.1.4 射線的能譜

    對(duì)醫(yī)用直線加速器X 射線能譜的研究,主要采用蒙特卡羅方法[2,10?11],如 Mohan 等[10]利用 EGS3 模擬計(jì)算了 Varian 不同機(jī)型的 X 射線能譜;Ding[2]利用EGS4 模擬計(jì)算了Varian2100EX 的X 射線能譜。在本研究中,劑量計(jì)算采用的X 射線能譜由BEAMnrc 對(duì) Varian TRILOGY 模擬得到,同時(shí)在計(jì)算中不考慮原射線和散射線能譜的差異。

    1.2 加速器劑量分布的測(cè)量

    測(cè)量采用0.13 cm3指形電離室(IBA,德國(guó))在三維水箱(Blue Phantom,IBA,德國(guó))測(cè)量 Varian TRILOGY 加速器6 MV 光子線在水中的劑量分布。該系統(tǒng)包括水箱、控制單元、探頭、儲(chǔ)水水箱、升降臺(tái)、計(jì)算機(jī)及控制軟件。水箱的體積為48 cm ×48 cm × 41cm,分 辨 率 為 0.1 mm,定 位 精 度為 ±0.1 mm,定 位 重 復(fù) 性 為 ±0.1 mm。 對(duì)3cm × 3cm、10 cm × 10 cm、40 cm × 40 cm 射野,分別測(cè)量了相應(yīng)的SCP、PDD 和5 cm、10 cm、20 cm 的OAR 用于模型擬合,并測(cè)量了6 cm × 6 cm、15cm ×15cm、20 cm × 20 cm 和 30 cm × 30 cm 射野的SCP對(duì)模型進(jìn)行檢驗(yàn)。測(cè)量中源皮距SSD= 100 cm,SCP測(cè)量的深度為5 cm。

    1.3 模型參數(shù)的確定

    對(duì)于某一個(gè)射野F,射野內(nèi)的各分項(xiàng)注量f0=對(duì)應(yīng)的 分 項(xiàng) 劑 量 分 布 分 別 為DF,0(r,z)、DF,1(r,z)、DF,2(r,z)和DF,3(r,z),則在此射野下的劑量為:

    式中:第一項(xiàng)包括了原射線的第一項(xiàng)和散射線項(xiàng)對(duì)劑量的貢獻(xiàn),其中模型有意義的是ɑ1、ɑ2、ɑ3、ɑ、b和= 1+c這6個(gè)參數(shù)。

    DF,0(r,z)、DF,1(r,z)、DF,2(r,z)和DF,3(r,z)可通過(guò)筆形束[12?13]、筒串卷積[14]或蒙特卡羅等方法計(jì)算得到,這里我們采用筆形束算法。

    筆形束算法的積分核由式(6)描述:

    其中:I(d)為積分核的深度分量,其徑向部分,除了權(quán)重因子wi(d)外,與深度d無(wú)關(guān)。分項(xiàng)劑量DF,k(x,y,d)可通過(guò)積分核與分項(xiàng)注量fk(x,y,d)的卷積得到:

    式中:fk(x',y',d)為深度d處對(duì)應(yīng)的第k項(xiàng)注量,k=0,1,2,3。在實(shí)際計(jì)算時(shí),我們采用一種發(fā)散型坐標(biāo)系(Fan-line coordinate)[13],在這種坐標(biāo)系下,積分核的深度分量I(d)和不同深度的徑向權(quán)重因子wi(d)是通過(guò)對(duì)EGSnrc 蒙特卡羅程序模擬得到的一系列不同半徑的錐束X 射線的劑量分布對(duì)式(6)進(jìn)行擬合而得到。采用發(fā)散型坐標(biāo)系,注量形式上與深度d無(wú)關(guān),因此劑量計(jì)算時(shí)不用對(duì)各個(gè)深度進(jìn)行卷積運(yùn)算,提高了計(jì)算效率。

    束流模型的參數(shù)通過(guò)對(duì)測(cè)量的劑量分布擬合得到,即通過(guò)調(diào)整參數(shù)使得目標(biāo)函數(shù)最小。目標(biāo)函數(shù)T為:

    式中:F代表射野為射野F下(ri,zj)處的測(cè)量劑量。當(dāng)設(shè)定某一個(gè)參考劑量時(shí),如10 cm ×10 cm射野最大劑量點(diǎn)處的劑量,它們可由測(cè)量得到的SCP、PDD 和OAR 確定??紤]到建成區(qū)的電子污染的影響,模型只對(duì)建成區(qū)之后的劑量進(jìn)行擬合。

    2 結(jié)果

    表1給出了模型參數(shù)擬合的結(jié)果。從擬合結(jié)果來(lái)看,當(dāng)射野從 3cm × 3cm 增加到 10 cm × 10 cm時(shí),式(5)中A的值從 0.896 增加到 0.998,增加了12%;而從10 cm × 10 cm增加到40 cm × 40 cm時(shí),A的值增加了8%,這表明小射野時(shí),散射的貢獻(xiàn)隨散射面積而較快增加。

    表1 模型參數(shù)的擬合結(jié)果Table 1 Fitting results of model parameters

    圖2給出了在不同射野下深度劑量分布的擬合值與測(cè)量值的比較。在建成區(qū)以后,擬合值與測(cè)量值在各個(gè)射野下都很好地符合。對(duì)于3cm × 3cm、10 cm × 10 cm、40 cm × 40 cm 射 野 ,從 水 下 深 度1.4 cm(最大劑量深度)到水下深度20 cm,劑量的擬合值與測(cè)量值的平均差異分別為1.3%、1.3%和1.1%;而從水下深度20 cm 后的平均差異分別為3.4%、2.0%和3.7%。水下深度20 cm后的差異較大的原因可能是擬合所用到的離軸劑量比對(duì)應(yīng)的深度局限于20 cm前,沒(méi)有對(duì)水下深度20 cm后的離軸劑量比進(jìn)行擬合。在建成區(qū)內(nèi),擬合值與測(cè)量值差異很大,這是因?yàn)榻ǔ蓞^(qū)內(nèi)存在電子污染,而模型不考慮污染電子。

    圖2 不同射野下深度劑量擬合值與測(cè)量值Fig.2 Comparison between fitting values and measured values of depth dose in the different fields

    圖3分別給出了水下深度5 cm和10 cm處的不同射野下離軸劑量的擬合值與測(cè)量值的比較。由圖3 可以看出,總體上擬合值與測(cè)量值差異不大。在5 cm 處 ,對(duì)于 3cm × 3cm、10 cm × 10 cm、40 cm ×40 cm射野,在OAR大于90%的區(qū)域,劑量的擬合值與測(cè)量值的平均差異分別為1.76%、1.30% 和0.65%;在10 cm處,分別為2.35%、0.75%和0.53%。

    圖3 不同射野水下5 cm(a)、10 cm(b)測(cè)量值與擬合值的離軸劑量比Fig.3 Comparison between fitting values and measured values of off-axis dose in the different fields at the depth of 5 cm(a)and 10 cm(b)underwater

    圖4 給出了不同射野的SCP的測(cè)量值與擬合值比較。由圖4 可以看出,總體上測(cè)量值與擬合值符合很好。3cm × 3cm射野的擬合值與測(cè)量值相差了1.42%,10 cm × 10 cm射野的擬合值與測(cè)量值相等,而40 cm × 40 cm 射野的擬合值與測(cè)量值相差0.55%。 對(duì) 于 作 為 檢 驗(yàn) 的 6 cm × 6 cm、15cm ×15cm、20 cm × 20 cm和30 cm × 30 cm射野,SCP的擬合值與測(cè)量值相差小于0.5%。

    3 討論

    X射線束流模型的建立應(yīng)該包含原射線和散射線。原射線通常設(shè)為點(diǎn)源。蒙特卡羅模擬表明散射光子最大的來(lái)源是均整器,其次是初級(jí)準(zhǔn)直器,而來(lái)自于次級(jí)準(zhǔn)直器(JAW、MLC)的散射光子數(shù)只占總散射光子數(shù)的0.2%[15]。因而對(duì)散射線建模時(shí),常將其設(shè)置為面源(均整器)和環(huán)源(初級(jí)準(zhǔn)直器),而將可運(yùn)動(dòng)并構(gòu)成不同射野的次級(jí)準(zhǔn)直部件(如JAWs和MLC)看成吸收體。但在計(jì)算束流注量時(shí)面源以及環(huán)源需將源劃分為點(diǎn)源計(jì)算后累加,這相對(duì)于點(diǎn)源大大增加了計(jì)算的時(shí)間??紤]到次級(jí)準(zhǔn)直部件對(duì)面源一些位置發(fā)射出光子的阻擋和吸收作用,將散射線建模為強(qiáng)度隨射野變化的點(diǎn)源可能是有益思路,它不僅具有面源對(duì)射野內(nèi)注量隨射野面積變化的特點(diǎn),同時(shí)減少劑量計(jì)算的時(shí)間。

    圖4 不同射野散射的總散射因子SCP的測(cè)量值與擬合值比較Fig.4 Comparison between fitted values and measured values of SCP in the different fields

    本文探討了一種簡(jiǎn)化的強(qiáng)度與射野的面積相關(guān)的單點(diǎn)源束流模型,此簡(jiǎn)化模型沒(méi)有考慮原射線和散射線能譜的差異、沒(méi)考慮原射線和散射線的點(diǎn)源位置的差異,對(duì)于射野外的注量,僅采用簡(jiǎn)單的指數(shù)下降函數(shù)來(lái)描述。盡管模型簡(jiǎn)單,但它卻很好地描述了射野內(nèi)的劑量分布,準(zhǔn)確地給出了SCP隨射野變化的關(guān)系,這說(shuō)明采用與射野面積相關(guān)的點(diǎn)源來(lái)描述散射線是可行的。

    4 結(jié)語(yǔ)

    本文初步探討了一個(gè)非常簡(jiǎn)單的單點(diǎn)源醫(yī)用加速器的X射線束流模型。該方法只是對(duì)采用單點(diǎn)源的模型初步探討,進(jìn)一步可以考慮原射線以及散射線的能譜差異,從而可進(jìn)一步提高模型的準(zhǔn)確度。而散射線的能譜,可通過(guò)蒙特卡羅模擬等方法來(lái)得到。在此模型中,單點(diǎn)源注量隨射野變化。結(jié)果表明:此射線束流模型可很好地重建體膜內(nèi)的三維劑量分布,給出正確的總散射因子。

    猜你喜歡
    射野點(diǎn)源束流
    利用三維水箱測(cè)量的“環(huán)形機(jī)架”加速器“典型射線數(shù)據(jù)”驗(yàn)證研究
    高能同步輻射光源低能束流輸運(yùn)線設(shè)計(jì)研究
    中國(guó)散裂中子源加速器注入束流損失調(diào)節(jié)研究
    關(guān)于脈沖積累對(duì)雙點(diǎn)源干擾影響研究
    靜止軌道閃電探測(cè)性能實(shí)驗(yàn)室驗(yàn)證技術(shù)研究
    三維藍(lán)水箱(BPH)掃描測(cè)量系統(tǒng)在螺旋斷層加速器質(zhì)量控制檢測(cè)中的應(yīng)用
    基于標(biāo)準(zhǔn)化點(diǎn)源敏感性的鏡面視寧度評(píng)價(jià)
    ECR中和器束流引出實(shí)驗(yàn)研究
    DAVID系統(tǒng)探測(cè)MLC葉片位置誤差的能力測(cè)試與評(píng)估
    自動(dòng)多野排序執(zhí)行技術(shù)在適形調(diào)強(qiáng)放療中的應(yīng)用
    男女无遮挡免费网站观看| 插逼视频在线观看| 免费观看在线日韩| 国产熟女午夜一区二区三区| 啦啦啦在线观看免费高清www| 亚洲 欧美一区二区三区| 在现免费观看毛片| 新久久久久国产一级毛片| 欧美激情 高清一区二区三区| 大话2 男鬼变身卡| 国产成人免费观看mmmm| 热re99久久国产66热| 美女视频免费永久观看网站| 欧美老熟妇乱子伦牲交| 1024视频免费在线观看| 少妇的逼好多水| 王馨瑶露胸无遮挡在线观看| 中文字幕免费在线视频6| 欧美老熟妇乱子伦牲交| 视频在线观看一区二区三区| 国产成人精品婷婷| 又黄又粗又硬又大视频| 国产无遮挡羞羞视频在线观看| 一区二区日韩欧美中文字幕 | av又黄又爽大尺度在线免费看| av卡一久久| 妹子高潮喷水视频| 中国美白少妇内射xxxbb| 日韩人妻精品一区2区三区| 七月丁香在线播放| 久久午夜福利片| 九九爱精品视频在线观看| 国产成人精品福利久久| 日本午夜av视频| 天堂8中文在线网| 日韩一区二区三区影片| 18+在线观看网站| 国产永久视频网站| 九色成人免费人妻av| 国产精品熟女久久久久浪| 日韩制服丝袜自拍偷拍| 国产免费现黄频在线看| 日韩精品有码人妻一区| 黄色毛片三级朝国网站| 久久久久久人人人人人| 亚洲国产日韩一区二区| 欧美日韩视频精品一区| 天堂8中文在线网| 在线观看www视频免费| a级片在线免费高清观看视频| 亚洲av中文av极速乱| 高清av免费在线| www日本在线高清视频| tube8黄色片| 国产av精品麻豆| 精品一区二区三区四区五区乱码 | 亚洲性久久影院| 免费人妻精品一区二区三区视频| 制服人妻中文乱码| a级毛色黄片| 免费看不卡的av| 1024视频免费在线观看| 制服人妻中文乱码| 韩国高清视频一区二区三区| 精品少妇黑人巨大在线播放| 亚洲熟女精品中文字幕| 综合色丁香网| 亚洲国产成人一精品久久久| 免费大片黄手机在线观看| 91国产中文字幕| 亚洲av电影在线观看一区二区三区| 99热6这里只有精品| 国产精品不卡视频一区二区| av播播在线观看一区| 两个人看的免费小视频| 丝袜在线中文字幕| 午夜福利网站1000一区二区三区| 在线观看三级黄色| 中文乱码字字幕精品一区二区三区| 日韩制服丝袜自拍偷拍| 在线观看一区二区三区激情| 欧美97在线视频| 丝袜喷水一区| kizo精华| 女人久久www免费人成看片| 最近最新中文字幕大全免费视频 | 精品国产一区二区久久| 日产精品乱码卡一卡2卡三| 午夜日本视频在线| 日韩三级伦理在线观看| 久久精品久久精品一区二区三区| 午夜91福利影院| 97人妻天天添夜夜摸| a级毛色黄片| 亚洲五月色婷婷综合| 国产深夜福利视频在线观看| 国产成人午夜福利电影在线观看| 精品熟女少妇av免费看| 国产精品成人在线| 少妇人妻 视频| 少妇人妻久久综合中文| 国产亚洲午夜精品一区二区久久| 99久久综合免费| 下体分泌物呈黄色| 激情视频va一区二区三区| 国产日韩欧美视频二区| 国产日韩欧美在线精品| 日韩熟女老妇一区二区性免费视频| 精品国产露脸久久av麻豆| 人人妻人人澡人人爽人人夜夜| 日产精品乱码卡一卡2卡三| 色吧在线观看| 免费观看性生交大片5| 日本av手机在线免费观看| av国产久精品久网站免费入址| 一区二区三区乱码不卡18| 午夜老司机福利剧场| 国产精品久久久久久精品古装| 黄色 视频免费看| 亚洲美女搞黄在线观看| 精品视频人人做人人爽| 黄色配什么色好看| 久久久久久久久久人人人人人人| 亚洲成人av在线免费| 男女下面插进去视频免费观看 | 国产精品一区二区在线不卡| 国产综合精华液| 国产 精品1| 国产在线一区二区三区精| 美女中出高潮动态图| 久热这里只有精品99| 看免费成人av毛片| 夫妻性生交免费视频一级片| 肉色欧美久久久久久久蜜桃| 汤姆久久久久久久影院中文字幕| 22中文网久久字幕| 国产一区二区在线观看av| 高清av免费在线| av国产久精品久网站免费入址| 欧美 亚洲 国产 日韩一| 热99国产精品久久久久久7| 亚洲av综合色区一区| videos熟女内射| 自线自在国产av| 老女人水多毛片| 午夜视频国产福利| videosex国产| 男女边摸边吃奶| 色婷婷av一区二区三区视频| 激情五月婷婷亚洲| 精品久久久精品久久久| 亚洲成人一二三区av| 男女无遮挡免费网站观看| 国产成人精品在线电影| 久久久亚洲精品成人影院| 亚洲av福利一区| 日韩欧美精品免费久久| 成年女人在线观看亚洲视频| 亚洲精品,欧美精品| 亚洲欧美日韩卡通动漫| 乱码一卡2卡4卡精品| 久久久亚洲精品成人影院| 91在线精品国自产拍蜜月| 十八禁高潮呻吟视频| 男的添女的下面高潮视频| 午夜免费鲁丝| 丁香六月天网| 日韩av在线免费看完整版不卡| xxxhd国产人妻xxx| 精品亚洲成国产av| 国产精品国产三级专区第一集| 亚洲色图综合在线观看| av在线老鸭窝| 超色免费av| 国产精品久久久久久av不卡| 欧美97在线视频| 日本av免费视频播放| 欧美日韩国产mv在线观看视频| 男女无遮挡免费网站观看| 欧美日韩视频精品一区| 国产亚洲精品久久久com| 国产淫语在线视频| 国产精品偷伦视频观看了| 国产极品粉嫩免费观看在线| 十八禁高潮呻吟视频| 成人午夜精彩视频在线观看| 成年美女黄网站色视频大全免费| 最黄视频免费看| 人人妻人人澡人人看| 亚洲综合色网址| 亚洲情色 制服丝袜| 母亲3免费完整高清在线观看 | 高清毛片免费看| 中文字幕精品免费在线观看视频 | 最新的欧美精品一区二区| 亚洲av欧美aⅴ国产| 大话2 男鬼变身卡| 免费高清在线观看视频在线观看| 日本wwww免费看| 黑人欧美特级aaaaaa片| 丰满少妇做爰视频| 侵犯人妻中文字幕一二三四区| 亚洲丝袜综合中文字幕| 最近中文字幕2019免费版| 婷婷色综合大香蕉| 久久ye,这里只有精品| 丰满少妇做爰视频| 高清欧美精品videossex| 老熟女久久久| 伦精品一区二区三区| 纵有疾风起免费观看全集完整版| xxx大片免费视频| 欧美日韩av久久| 精品亚洲乱码少妇综合久久| 久久久国产欧美日韩av| 亚洲丝袜综合中文字幕| 国产av精品麻豆| 亚洲天堂av无毛| 国产高清国产精品国产三级| 中文乱码字字幕精品一区二区三区| 国产乱来视频区| 日韩精品免费视频一区二区三区 | 最近中文字幕2019免费版| 亚洲一码二码三码区别大吗| 亚洲美女视频黄频| 国产精品三级大全| 我要看黄色一级片免费的| 亚洲少妇的诱惑av| 看十八女毛片水多多多| 草草在线视频免费看| 国产精品国产三级国产av玫瑰| 亚洲欧美日韩卡通动漫| videos熟女内射| 国产av精品麻豆| 99热网站在线观看| 国产 精品1| 国产成人午夜福利电影在线观看| 久久久亚洲精品成人影院| 成人18禁高潮啪啪吃奶动态图| 99热6这里只有精品| kizo精华| 最新中文字幕久久久久| 一级毛片我不卡| 午夜精品国产一区二区电影| 久久精品aⅴ一区二区三区四区 | 成人二区视频| 亚洲av在线观看美女高潮| 男人添女人高潮全过程视频| 亚洲经典国产精华液单| 日本与韩国留学比较| 丰满迷人的少妇在线观看| 人妻人人澡人人爽人人| 亚洲美女搞黄在线观看| 人体艺术视频欧美日本| 成人国产麻豆网| 成人无遮挡网站| 久久这里有精品视频免费| 一区二区三区精品91| 日日摸夜夜添夜夜爱| av黄色大香蕉| 边亲边吃奶的免费视频| 多毛熟女@视频| 久久久久网色| 晚上一个人看的免费电影| 熟女av电影| 成人毛片a级毛片在线播放| 少妇人妻 视频| 青春草亚洲视频在线观看| 麻豆精品久久久久久蜜桃| 99国产精品免费福利视频| 亚洲精品乱久久久久久| 午夜精品国产一区二区电影| 一区二区三区精品91| 成年人午夜在线观看视频| 久久久国产精品麻豆| 久久久久人妻精品一区果冻| 成人综合一区亚洲| 精品一区二区三区四区五区乱码 | 夫妻性生交免费视频一级片| 一区二区三区四区激情视频| 亚洲欧美成人精品一区二区| 大片电影免费在线观看免费| 精品亚洲乱码少妇综合久久| 婷婷色综合大香蕉| 韩国高清视频一区二区三区| 男女啪啪激烈高潮av片| 夜夜骑夜夜射夜夜干| 国产成人精品福利久久| 人人澡人人妻人| 男女无遮挡免费网站观看| av片东京热男人的天堂| 中文字幕亚洲精品专区| 91aial.com中文字幕在线观看| a 毛片基地| av.在线天堂| 美女福利国产在线| 亚洲,欧美,日韩| 久久久精品94久久精品| 国产一区有黄有色的免费视频| 欧美bdsm另类| 美女视频免费永久观看网站| 欧美日韩成人在线一区二区| 最近的中文字幕免费完整| 亚洲欧美一区二区三区黑人 | 亚洲av在线观看美女高潮| av免费观看日本| 黄网站色视频无遮挡免费观看| 交换朋友夫妻互换小说| 中文字幕精品免费在线观看视频 | 欧美日韩av久久| 精品熟女少妇av免费看| 一区二区三区乱码不卡18| 在线天堂中文资源库| 欧美性感艳星| 国产午夜精品一二区理论片| 少妇的逼好多水| 久久久久久久亚洲中文字幕| 日本vs欧美在线观看视频| 久久久久久久国产电影| 美女主播在线视频| 免费看av在线观看网站| 91国产中文字幕| 最近的中文字幕免费完整| 午夜福利网站1000一区二区三区| 国产在视频线精品| 亚洲av电影在线进入| 看免费av毛片| 高清视频免费观看一区二区| 国产精品一区二区在线不卡| 热99久久久久精品小说推荐| 欧美少妇被猛烈插入视频| 免费少妇av软件| 国产精品国产三级国产av玫瑰| 亚洲伊人色综图| 成人国产麻豆网| 在线亚洲精品国产二区图片欧美| 日本av免费视频播放| 黄网站色视频无遮挡免费观看| 日韩av在线免费看完整版不卡| 精品一区二区三区视频在线| 成人毛片60女人毛片免费| 午夜福利影视在线免费观看| 老司机影院毛片| 丝袜美足系列| 免费av不卡在线播放| 91精品伊人久久大香线蕉| 国产在线视频一区二区| 久久精品久久精品一区二区三区| 久久av网站| 多毛熟女@视频| 国产精品久久久av美女十八| 黑人高潮一二区| 免费黄频网站在线观看国产| 一本色道久久久久久精品综合| 精品人妻一区二区三区麻豆| 亚洲伊人久久精品综合| 涩涩av久久男人的天堂| 亚洲一级一片aⅴ在线观看| 精品亚洲成国产av| 中国三级夫妇交换| 精品亚洲乱码少妇综合久久| 七月丁香在线播放| 亚洲国产毛片av蜜桃av| 亚洲精品乱久久久久久| 国产xxxxx性猛交| 777米奇影视久久| 丝袜喷水一区| 国产男人的电影天堂91| 亚洲国产毛片av蜜桃av| 久久人人爽人人爽人人片va| 久久这里有精品视频免费| 国产精品欧美亚洲77777| 久久久精品免费免费高清| 国国产精品蜜臀av免费| 久久久欧美国产精品| 国产视频首页在线观看| 国产成人一区二区在线| 免费观看性生交大片5| 亚洲久久久国产精品| 午夜免费鲁丝| 少妇高潮的动态图| 丝袜美足系列| 欧美日韩一区二区视频在线观看视频在线| 99久久精品国产国产毛片| 你懂的网址亚洲精品在线观看| 一本大道久久a久久精品| 在线天堂最新版资源| 大香蕉97超碰在线| 精品福利永久在线观看| 亚洲欧美日韩另类电影网站| 国产精品成人在线| 18禁动态无遮挡网站| 又大又黄又爽视频免费| 99久久人妻综合| 亚洲国产av影院在线观看| 新久久久久国产一级毛片| 成人国语在线视频| 国产又爽黄色视频| 亚洲欧洲日产国产| 日韩在线高清观看一区二区三区| 欧美日韩视频高清一区二区三区二| 免费看光身美女| 国产1区2区3区精品| 国产一区二区激情短视频 | 亚洲伊人色综图| 国产国语露脸激情在线看| 自线自在国产av| 国产成人91sexporn| 久久久a久久爽久久v久久| 亚洲av福利一区| 一区二区日韩欧美中文字幕 | 亚洲国产看品久久| 男女边吃奶边做爰视频| 午夜福利,免费看| 日韩成人av中文字幕在线观看| 日韩,欧美,国产一区二区三区| 亚洲,一卡二卡三卡| 波野结衣二区三区在线| 夜夜骑夜夜射夜夜干| 国产免费一级a男人的天堂| 日韩 亚洲 欧美在线| 只有这里有精品99| 午夜精品国产一区二区电影| av免费在线看不卡| 亚洲精品久久成人aⅴ小说| 国产视频首页在线观看| 日韩,欧美,国产一区二区三区| 精品少妇内射三级| 肉色欧美久久久久久久蜜桃| 国产黄色视频一区二区在线观看| 制服诱惑二区| 高清毛片免费看| 爱豆传媒免费全集在线观看| 秋霞伦理黄片| 午夜免费鲁丝| 大陆偷拍与自拍| 色婷婷久久久亚洲欧美| 蜜臀久久99精品久久宅男| 热re99久久精品国产66热6| 久久精品久久久久久噜噜老黄| 亚洲少妇的诱惑av| 日本猛色少妇xxxxx猛交久久| 97超碰精品成人国产| 久久久久久久精品精品| 亚洲经典国产精华液单| 亚洲精品国产av成人精品| 国产女主播在线喷水免费视频网站| 少妇熟女欧美另类| 美女主播在线视频| av在线老鸭窝| 国产精品一区二区在线不卡| videosex国产| 欧美 亚洲 国产 日韩一| 日韩一本色道免费dvd| 亚洲av综合色区一区| 亚洲国产av新网站| 精品视频人人做人人爽| 午夜免费鲁丝| 香蕉丝袜av| 精品一区二区三区视频在线| 亚洲美女黄色视频免费看| 丝瓜视频免费看黄片| 国产一区亚洲一区在线观看| 国产精品国产av在线观看| 色婷婷久久久亚洲欧美| 欧美精品人与动牲交sv欧美| 大香蕉97超碰在线| 午夜视频国产福利| 一级a做视频免费观看| 少妇人妻精品综合一区二区| 夫妻性生交免费视频一级片| 丝袜美足系列| 性高湖久久久久久久久免费观看| 观看美女的网站| 国产伦理片在线播放av一区| 黑人欧美特级aaaaaa片| 国产伦理片在线播放av一区| 青春草亚洲视频在线观看| 国产不卡av网站在线观看| 成人影院久久| 80岁老熟妇乱子伦牲交| 国产精品免费大片| 夜夜爽夜夜爽视频| 免费久久久久久久精品成人欧美视频 | 亚洲av福利一区| 精品亚洲乱码少妇综合久久| 色吧在线观看| 午夜免费男女啪啪视频观看| 91精品伊人久久大香线蕉| 男女免费视频国产| 建设人人有责人人尽责人人享有的| 99热国产这里只有精品6| 久久久a久久爽久久v久久| 国产视频首页在线观看| 90打野战视频偷拍视频| xxxhd国产人妻xxx| 香蕉丝袜av| 秋霞在线观看毛片| 菩萨蛮人人尽说江南好唐韦庄| 青春草视频在线免费观看| 国产精品女同一区二区软件| 亚洲成人av在线免费| av电影中文网址| 这个男人来自地球电影免费观看 | 日韩一区二区三区影片| 精品一区二区免费观看| 日韩欧美一区视频在线观看| av国产精品久久久久影院| 天堂中文最新版在线下载| 天堂俺去俺来也www色官网| 日韩制服骚丝袜av| 国产又色又爽无遮挡免| 美女福利国产在线| 日韩av不卡免费在线播放| 成年av动漫网址| av又黄又爽大尺度在线免费看| 18禁国产床啪视频网站| 午夜老司机福利剧场| 最新的欧美精品一区二区| 亚洲av在线观看美女高潮| 丰满乱子伦码专区| 亚洲av电影在线进入| 狂野欧美激情性bbbbbb| 欧美日韩精品成人综合77777| 国产xxxxx性猛交| 亚洲精品久久午夜乱码| 永久网站在线| 免费看av在线观看网站| xxx大片免费视频| 久久久久精品人妻al黑| 亚洲精品日本国产第一区| 国产熟女午夜一区二区三区| 亚洲三级黄色毛片| 免费观看性生交大片5| 国产福利在线免费观看视频| 丝袜在线中文字幕| 国产亚洲欧美精品永久| 看免费成人av毛片| 黄色视频在线播放观看不卡| 人人妻人人澡人人看| 80岁老熟妇乱子伦牲交| 精品卡一卡二卡四卡免费| 亚洲激情五月婷婷啪啪| 亚洲少妇的诱惑av| 9色porny在线观看| 欧美日韩一区二区视频在线观看视频在线| 啦啦啦啦在线视频资源| 国产男人的电影天堂91| 大香蕉久久网| 亚洲欧美一区二区三区国产| 日本vs欧美在线观看视频| a级片在线免费高清观看视频| 一本久久精品| 波多野结衣一区麻豆| 男女无遮挡免费网站观看| 一二三四中文在线观看免费高清| 插逼视频在线观看| 97人妻天天添夜夜摸| 国产精品偷伦视频观看了| 999精品在线视频| 国产极品粉嫩免费观看在线| 国产亚洲av片在线观看秒播厂| 欧美精品国产亚洲| 高清视频免费观看一区二区| 欧美国产精品va在线观看不卡| 欧美日韩视频精品一区| 久久精品熟女亚洲av麻豆精品| 91精品伊人久久大香线蕉| 2022亚洲国产成人精品| 在线观看人妻少妇| 免费观看性生交大片5| 母亲3免费完整高清在线观看 | h视频一区二区三区| 一二三四中文在线观看免费高清| 1024视频免费在线观看| 午夜福利视频在线观看免费| 久久精品久久精品一区二区三区| 母亲3免费完整高清在线观看 | 久久久欧美国产精品| 最近手机中文字幕大全| 男人操女人黄网站| 久久精品aⅴ一区二区三区四区 | 秋霞在线观看毛片| 亚洲精品乱久久久久久| 九草在线视频观看| 亚洲五月色婷婷综合| 爱豆传媒免费全集在线观看| 岛国毛片在线播放| 一边摸一边做爽爽视频免费| 最近最新中文字幕免费大全7| 捣出白浆h1v1| 国产 精品1| 桃花免费在线播放| 不卡视频在线观看欧美| 亚洲人与动物交配视频| 十分钟在线观看高清视频www| 亚洲精品美女久久av网站| 一二三四在线观看免费中文在 | 免费黄网站久久成人精品| 久久久久人妻精品一区果冻| 国产免费福利视频在线观看| 九九爱精品视频在线观看| 国产日韩一区二区三区精品不卡| 在线看a的网站| 99久久人妻综合| 视频中文字幕在线观看| 国产精品 国内视频| tube8黄色片| 久久国产亚洲av麻豆专区| 只有这里有精品99| 日韩精品免费视频一区二区三区 | 大香蕉久久网| 青春草国产在线视频| 久久久国产一区二区| 波野结衣二区三区在线|