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

    多源γ 射束立體定向放射治療計(jì)劃系統(tǒng)的劑量計(jì)算加速算法研究

    2012-12-05 09:37:48曹國(guó)剛
    中國(guó)醫(yī)療器械雜志 2012年1期
    關(guān)鍵詞:計(jì)算速度錐形線程

    【作 者】曹國(guó)剛

    上海世鵬實(shí)驗(yàn)室科技發(fā)展有限公司,上海,200072

    1951年,瑞典著名神經(jīng)外科專家Leksell教授提出立體定向放射外科(Stereotactic Radiosurgery)的概念。1968年,Leksell教授與Larsson、Rexed等核物理專家以及放射治療學(xué)專家、解剖學(xué)家通力合作,在瑞典的Uppsala大學(xué)首先研制成功世界上第一臺(tái)伽瑪?shù)叮_了立體定向放射神經(jīng)外科的新紀(jì)元[1-2]。立體定向放射治療是指采用立體定向手段進(jìn)行的高精度放射治療的技術(shù)。與傳統(tǒng)的放射治療相比,它在增加腫瘤照射劑量的同時(shí),分散了對(duì)健康組織的照射劑量,大大提高了腫瘤放射治療的臨床療效。經(jīng)過(guò)各國(guó)專家?guī)资甑呐?,立體定向放射治療現(xiàn)已形成了靜態(tài)式、旋轉(zhuǎn)式和弧形等中心三大類方式,出現(xiàn)了頭部伽瑪?shù)?、旋轉(zhuǎn)式頭部和體部伽瑪?shù)逗皖^部和體部X刀等多種裝置。

    放射治療的質(zhì)量取決于治療時(shí)的吸收劑量分布,而劑量計(jì)算是目前能全面給出吸收劑量分布的唯一方法,因此高精度地計(jì)算人體的三維吸收劑量分布就成為評(píng)價(jià)治療計(jì)劃優(yōu)劣的先決條件。放射治療計(jì)劃系統(tǒng)中,劑量計(jì)算總是在計(jì)算速度和計(jì)算精度間尋找平衡點(diǎn)?;贛onte Carlo粒子輸運(yùn)理論的劑量計(jì)算是目前最精確的算法,但計(jì)算所需的時(shí)間是驚人的,極少能用于臨床[3-4];基于Collapsed Cone卷積疊加模型的劑量計(jì)算方法,常用于X射線計(jì)劃系統(tǒng)[5-6];而查表插值法等近似模型則適合 γ - 射線的性質(zhì)及照射方法[7]。

    本文以陀螺刀(一種陀螺旋轉(zhuǎn)式多源 γ 射束立體定向放射治療系統(tǒng))治療計(jì)劃系統(tǒng)為例,以查表插值法近似模型為基礎(chǔ),介紹一種基于錐形坐標(biāo)系和多核并行化的劑量計(jì)算加速算法。該治療計(jì)劃系統(tǒng)由上海世鵬實(shí)驗(yàn)室科技發(fā)展有限公司自主研發(fā),已通過(guò)國(guó)家相關(guān)機(jī)構(gòu)檢驗(yàn)。

    1 多源 γ 射束立體定向放射治療計(jì)劃系統(tǒng)

    多源 γ 射束立體定向放射治療系統(tǒng)(Multi-Source Stereotactic Radiotherapy System with Gamma Beam),俗稱伽瑪?shù)?Gamma Knife),采用幾何聚焦的原理,利用立體定向裝置、CT模擬定位、計(jì)算機(jī)立體成像與三維治療計(jì)劃系統(tǒng)等技術(shù)手段,將由多個(gè)鈷-60放射源經(jīng)準(zhǔn)直器形成的很細(xì)的γ-射線在三維空間聚焦于一點(diǎn),要求該點(diǎn)與病灶完全重合,形成一個(gè)高劑量分布區(qū),稱之為靶點(diǎn)(shot)。由于靶點(diǎn)上的能量極高,可以在短時(shí)間內(nèi)將靶點(diǎn)擊毀;而靶點(diǎn)之外的健康組織所受的照射劑量卻很小,從而達(dá)到了比手術(shù)切除更好的治療效果[2,8]。

    治療計(jì)劃系統(tǒng)(treatment planning system,TPS)是伽瑪?shù)兜闹匾M成部分,是一種融合了諸如圖像處理和分析、數(shù)值計(jì)算和可視化、計(jì)算機(jī)圖形學(xué)、軟件工程、放射物理學(xué)以及腫瘤學(xué)等多種學(xué)科的復(fù)雜軟件系統(tǒng),醫(yī)生藉此可以方便迅速地制定出合理的治療計(jì)劃方案。治療計(jì)劃系統(tǒng)一般由以下幾個(gè)功能模塊組成:劑量學(xué)數(shù)據(jù)管理、病例數(shù)據(jù)管理、病例數(shù)據(jù)導(dǎo)入和導(dǎo)出、框架定位、圖像配準(zhǔn)與融合、輪廓勾畫、虛擬人體建模、計(jì)劃規(guī)劃與優(yōu)化、劑量計(jì)算與計(jì)劃評(píng)估、治療計(jì)劃輸出與存檔等[9]。

    目前的伽瑪?shù)吨委熡?jì)劃系統(tǒng)均采用正向設(shè)計(jì)、手工規(guī)劃的方法,計(jì)劃者負(fù)責(zé)選擇靶點(diǎn)相關(guān)的參數(shù),由計(jì)算機(jī)進(jìn)行劑量計(jì)算。這是一個(gè)反復(fù)試誤(try and error)的過(guò)程,在找到合適的方案之前要一直進(jìn)行下去。這種方法中劑量計(jì)算速度決定計(jì)劃設(shè)計(jì)時(shí)間[10],因此提高劑量計(jì)算速度成為計(jì)劃系統(tǒng)的關(guān)鍵。

    2 治療計(jì)劃系統(tǒng)中的劑量計(jì)算

    所謂劑量計(jì)算,就是計(jì)算粒子或光子射線照射到人體(或體模)中能量吸收的空間分布。由于粒子或光子線與人體組織的相互作用比較復(fù)雜,射線的種類、能量和照射方式不同,其形成的劑量分布會(huì)有很大的差別。在伽瑪?shù)吨委熡?jì)劃系統(tǒng)中,根據(jù) γ - 射線的性質(zhì)及照射方法,采用劑量計(jì)算近似模型,用實(shí)際測(cè)量的劑量學(xué)數(shù)據(jù)進(jìn)行查表插值的方法計(jì)算劑量分布的近似估計(jì)。該方法主要包括單源劑量計(jì)算、多源劑量分布疊加和多靶點(diǎn)劑量分布疊加。

    2.1 單源劑量計(jì)算方法

    單個(gè)靜態(tài) γ - 射線源對(duì)病人體內(nèi)某一點(diǎn)貢獻(xiàn)的劑量可以用下式表示:

    式中:D(x,y,z,p,n)為當(dāng)源在位置p,使用第n個(gè)準(zhǔn)直器照射時(shí),體內(nèi)某點(diǎn)(x,y,z)的劑量值;x,y,z為病人體內(nèi)某點(diǎn)的坐標(biāo);p為源位置;n為準(zhǔn)直器編號(hào);DF(n)為第n個(gè)準(zhǔn)直器的標(biāo)定劑量率;TMR(deff,n)為第n個(gè)準(zhǔn)直器,深度為deff時(shí)的組織最大劑量比(tissue maximum ratio,TMR);deff為計(jì)算點(diǎn)的有效深度(包含非均勻組織的電子密度校正);OAR(d,r,n)為第n個(gè)準(zhǔn)直器,離軸距離為r,深度為d的離軸比(off axis ratio,為深度校正幾何因子;SAD為源軸距,即源至等中心的距離;△d為計(jì)算點(diǎn)在中心軸上投影距離等中心的距離,遠(yuǎn)離源方向時(shí)取正,反之取負(fù)。

    2.2 多源劑量計(jì)算方法

    對(duì)于多個(gè)放射源或放射源為運(yùn)動(dòng)的情況,采用離散累加的方式計(jì)算劑量,所以計(jì)算點(diǎn)的實(shí)際劑量是多個(gè)放射源位置的劑量貢獻(xiàn)之和:

    其中,D(x,y,z,n)為使用第n個(gè)準(zhǔn)直器照射時(shí),體內(nèi)某點(diǎn)(x,y,z)的劑量值;Wp是源在位置p的權(quán)重。

    2.3 多靶點(diǎn)劑量計(jì)算方法

    如果有m個(gè)靶點(diǎn),病人體內(nèi)的劑量是多個(gè)靶點(diǎn)的加權(quán)和,即:

    其中,D(x,y,z)為體內(nèi)某點(diǎn)(x,y,z)的劑量值,wi是第i個(gè)靶點(diǎn)的權(quán)重,ni為第i個(gè)靶點(diǎn)的準(zhǔn)直器編號(hào)。

    3 劑量計(jì)算加速算法

    傳統(tǒng)劑量計(jì)算方法通常在直角坐標(biāo)系中進(jìn)行,不同計(jì)算點(diǎn)需重復(fù)計(jì)算深度、有效深度、離軸距離等中間結(jié)果多次,且這些中間結(jié)果無(wú)法重復(fù)使用,導(dǎo)致計(jì)算速度慢。因此,我們考慮建立一個(gè)“錐形”坐標(biāo)系,重復(fù)利用中間結(jié)果,提高計(jì)算速度。此外,劑量計(jì)算過(guò)程中,需多次計(jì)算源在不同位置的劑量分布,是一個(gè)典型的計(jì)算密集型過(guò)程,十分耗時(shí)。隨著多核處理器的可用性和性價(jià)比的不斷提高,越來(lái)越多的應(yīng)用將采用多核架構(gòu),多核并行化的劑量計(jì)算的研究更具有實(shí)用價(jià)值。

    3.1 基于錐形坐標(biāo)系的劑量計(jì)算

    首先建立錐形坐標(biāo)系。通過(guò)空間幾何變換將機(jī)器坐標(biāo)系變換到以源為原點(diǎn)的直角坐標(biāo)系(x,y,z),錐形坐標(biāo)系(d,φ,θ)與上述直角坐標(biāo)系關(guān)系如圖1。兩個(gè)坐標(biāo)系的轉(zhuǎn)換公式如下:

    其中,當(dāng)z=0時(shí),θ=0,φ=0;當(dāng)x=0,y>0時(shí),θ=90;當(dāng)x=0,y<0時(shí),θ=270。

    在錐形坐標(biāo)系中進(jìn)行劑量計(jì)算時(shí),可以減少對(duì)deff、d、r和△d等中間結(jié)果的大量重復(fù)計(jì)算,加快計(jì)算速度。在該坐標(biāo)系下進(jìn)行劑量計(jì)算后,只需將其結(jié)果轉(zhuǎn)換回直角坐標(biāo)系下即可。

    圖1 錐形坐標(biāo)系與直角坐標(biāo)系Fig.1 Cone coordingate system and cartesian coordinate system

    3.2 劑量計(jì)算的多核并行化

    多核技術(shù)是近年來(lái)全球計(jì)算機(jī)技術(shù)發(fā)展的重要內(nèi)容,自從英特爾在2006年底發(fā)布了全球第一款主流四核處理器至強(qiáng)5300后,計(jì)算機(jī)行業(yè)宣告正式進(jìn)入多核時(shí)代。多核計(jì)算將成為一種普及的計(jì)算模式,并將深刻影響企業(yè)和消費(fèi)者的使用模式[11-12]。為了提高劑量計(jì)算速度,本文提出一種多核并行化的劑量計(jì)算算法,以加快劑量計(jì)算速度。

    多核并行化將每個(gè)單源劑量計(jì)算分配給一個(gè)虛擬的處理器(核),所有處理器共享輸入數(shù)據(jù),形成數(shù)據(jù)并行的算法。每個(gè)虛擬的處理器處理對(duì)應(yīng)的計(jì)算過(guò)程,無(wú)需關(guān)心其他的計(jì)算過(guò)程,但需在對(duì)共享數(shù)據(jù)進(jìn)行讀寫時(shí),采用同步處理保證對(duì)共享數(shù)據(jù)的互斥訪問(wèn)。

    執(zhí)行共享數(shù)據(jù)的多核并行計(jì)算時(shí),線程不是孤立的,各線程間存在依賴關(guān)系,線程間對(duì)數(shù)據(jù)的共享可能會(huì)導(dǎo)致錯(cuò)誤的行為和結(jié)果。因此,應(yīng)用多線程實(shí)現(xiàn)并行計(jì)算必需引入同步機(jī)制協(xié)調(diào)各線程,以便安全共享數(shù)據(jù)。線程的同步機(jī)制包括臨界區(qū)(Critical sections)、互斥量(Mutex)、信號(hào)量(Semaphores)等。其中,臨界區(qū)和互斥量用來(lái)串行化多個(gè)線程共享數(shù)據(jù)的訪問(wèn),其思想是每個(gè)獨(dú)占性地訪問(wèn)一個(gè)資源的線程,可在訪問(wèn)那個(gè)資源之前鎖定臨界區(qū),訪問(wèn)之后解除鎖定。如果另一個(gè)線程也試圖鎖定該臨界區(qū),則該線程將被阻塞直到該臨界區(qū)空閑。信號(hào)量保存有代表可用資源數(shù)量的資源數(shù),鎖定信號(hào)量會(huì)減少資源數(shù),釋放信號(hào)量則增加資源數(shù)。只有在線程試圖鎖定資源數(shù)為0的信號(hào)量時(shí),線程才會(huì)被阻塞。信號(hào)量可以用來(lái)同步化同一進(jìn)程中的線程,也可同步化不同進(jìn)程中的線程。當(dāng)系統(tǒng)中有多個(gè)同種類型的資源共享時(shí),就可利用信號(hào)量在多個(gè)線程間分配這些資源。

    基于多核處理器的并行化劑量計(jì)算是共享輸入數(shù)據(jù)的并行化算法,其并行程序與串行程序十分相似,僅在同時(shí)計(jì)算累加劑量分布數(shù)據(jù)時(shí)用臨界區(qū)保證互斥訪問(wèn)。

    4 實(shí)驗(yàn)結(jié)果與分析

    4.1 實(shí)驗(yàn)條件

    實(shí)驗(yàn)采用直徑160 mm的有機(jī)玻璃球模,并將根據(jù)該球模產(chǎn)生的仿真CT數(shù)據(jù)輸入TPS;采用標(biāo)稱尺寸為12 mm,標(biāo)稱吸收劑量率為200cGy/min的聚焦野劑量學(xué)數(shù)據(jù);將靶點(diǎn)置于球模中心,滾筒旋轉(zhuǎn)120o,持續(xù)照射5分鐘。采用不同方法進(jìn)行劑量計(jì)算,比較不同方法的點(diǎn)劑量誤差、面積重合率和計(jì)算速度等參數(shù)。

    實(shí)驗(yàn)環(huán)境:HP xw6600圖形工作站;2顆英特爾至強(qiáng)E5405四核處理器,4G內(nèi)存;Windows XP Professional SP3。

    4.2 評(píng)價(jià)方式

    分別在直角坐標(biāo)系和錐形坐標(biāo)系下計(jì)算劑量分布,對(duì)比劑量分布中的點(diǎn)劑量誤差和面積重合率。用點(diǎn)劑量誤差評(píng)價(jià)某個(gè)固定點(diǎn)在兩種計(jì)算方法中絕對(duì)劑量的差異,而面積重合率評(píng)價(jià)兩種計(jì)算方法中相對(duì)劑量分布的差異。

    點(diǎn)劑量絕對(duì)誤差和相對(duì)誤差定義分別如下:

    其中,△D為某點(diǎn)的劑量絕對(duì)誤差,△DR為該點(diǎn)的劑量相對(duì)誤差,D1和D2分別為在直角坐標(biāo)系和錐形坐標(biāo)系下計(jì)算所得的該點(diǎn)劑量值。

    面積重合率定義如下:

    其中,△S為某指定等劑量線的面積重合率,S1和S2分別為在直角坐標(biāo)系和錐形坐標(biāo)系下計(jì)算所得的該等劑量線所包圍區(qū)域的面積,SC為S1和S2兩部分面積的重合部分。

    4.3 不同坐標(biāo)系下計(jì)算結(jié)果與分析

    根據(jù)過(guò)球模中心CT層面的劑量分布結(jié)果,統(tǒng)計(jì)點(diǎn)劑量誤差和面積重合率。兩種坐標(biāo)系下劑量計(jì)算結(jié)果的點(diǎn)劑量誤差和面積重合率如表1和表2,在兩種坐標(biāo)系下靶點(diǎn)計(jì)算量與計(jì)算時(shí)間結(jié)果如表3,表中各符號(hào)含義同4.2節(jié)。

    從上述數(shù)據(jù)可以看出,兩種計(jì)算方法在點(diǎn)劑量誤差和面積重合率兩個(gè)方面的差異均非常小,可以忽略不計(jì);且新方法可以極大提高計(jì)算速度,與原有方法相比,該方法可以節(jié)約65%的時(shí)間。

    表1 點(diǎn)劑量誤差Tab.1 Error of point dose

    表2 面積重合率Tab.2 Overlap rate of area

    表3 靶點(diǎn)劑量和計(jì)算時(shí)間Tab.3 Point dose of shot center and calculation time

    表4 多核并行計(jì)算時(shí)間與速度Tab.4 Time and speed in multi-core based parallelized dose calculation

    4.4 多核并行化劑量計(jì)算結(jié)果與分析

    表4列出了兩種計(jì)算方法使用不同數(shù)量核計(jì)算所需的時(shí)間和速度。其中“傳統(tǒng)”表示不采用多核并行的傳統(tǒng)方法,表中分別將兩種方法的傳統(tǒng)計(jì)算方法速度定義為1.0。

    圖2為計(jì)算速度與核數(shù)量之間的關(guān)系,圖中橫坐標(biāo)為劑量計(jì)算使用核的數(shù)量,縱坐標(biāo)為計(jì)算速度,虛線、點(diǎn)虛線、實(shí)線分別代表直角坐標(biāo)系、錐形坐標(biāo)系和理想情況。

    圖2 多核并行化劑量計(jì)算中速度與核數(shù)量的關(guān)系Fig.2 Relationship of speed and number of cores in multi-core based parallelized dose calculation

    從上述圖表可以看出,多核并行化的劑量計(jì)算只使用一個(gè)核進(jìn)行計(jì)算時(shí),與不采用多核并行的傳統(tǒng)計(jì)算方法耗時(shí)相當(dāng),無(wú)明顯差異;多核并行化的劑量計(jì)算速度與所使用核的數(shù)量接近正比,但比理想情況稍差,原因是多線程運(yùn)算時(shí),線程間切換和數(shù)據(jù)同步需要占用少量的時(shí)間。因此,多核并行化的劑量計(jì)算方法更進(jìn)一步提高了劑量計(jì)算速度,其速度提升與使用的核接近正比,當(dāng)采用8核進(jìn)行并行計(jì)算時(shí),其計(jì)算速度可提升7-8倍。

    5 小結(jié)

    本文提出一種基于錐形坐標(biāo)系和多核并行化的劑量計(jì)算加速算法?;阱F形坐標(biāo)系,該方法在保證計(jì)算精度的前提下可以極大提高劑量計(jì)算速度;同時(shí)利用多核并行化計(jì)算方法,還可進(jìn)一步提高劑量計(jì)算的速度,其速度提升與使用的核接近正比。從實(shí)驗(yàn)結(jié)果和分析看出,該方法可以將計(jì)算時(shí)間從傳統(tǒng)計(jì)算方法的584.3秒減少到30.3秒,節(jié)約95%的時(shí)間,即速度提升到傳統(tǒng)方法的20倍左右。在實(shí)際應(yīng)用中,還可以通過(guò)減小計(jì)算區(qū)域和降低計(jì)算網(wǎng)格精度等方法,得到劑量計(jì)算的快算結(jié)果,在可接受范圍內(nèi)降低計(jì)算精度,進(jìn)一步提高計(jì)算速度。

    [1]陳炳桓.立體定向放射[M].北京: 北京出版社,1994

    [2]嚴(yán)玉龍.立體定向放射療法的治療規(guī)劃——適形放療與立體定向放射外科治療規(guī)劃的理論與技術(shù)研究[D].東南大學(xué)博士論文,1997

    [3]Kawrakow I,Fippel M.Investigation of variance reduction techniques for Monte Carlo photon dose calculation using XVMC[J].Phys.Med.Biol,2000,45: 2163-2183.

    [4]Kawrakow I,Fippel M.VMC++,a MC algorithm optimized for electron and photon beam dose calculations for RTP[C].Proceedings of the 22th annual EMBS international conference,2000,1490-1493.

    [5]Miften M,Wiesmeyer M,Monthofer S,et al.Implement at ion of FFT convolution and multigrid superposition models in the FOCUS RTP system[J].Phys.Med.Biol.,2000,45(4): 817- 833.

    [6]Zhengdong Zhou,Wei Song.A Study of Dose Calculation Based on Collapsed Cone Convolution Algorithm[C].BMEI2010,2010,1331-1335.

    [7]羅浩,郭吉豐,劉涵.CT機(jī)器人伽瑪?shù)吨谢谌S區(qū)域生長(zhǎng)算法的劑量快速搜索計(jì)算研究[J].中國(guó)醫(yī)療器械雜志,2005,29(3):161-163.

    [8]于金明,殷蔚伯,李寶生主編.腫瘤精確放射治療學(xué)[M].山東科學(xué)技術(shù)出版社,2004

    [9]於文雪,李松毅,羅立民.鈷-60三維精確適形放射治療計(jì)劃系統(tǒng)的研制.中國(guó)醫(yī)療器械雜志[J],2007,31(6): 407-410.

    [10]周正東,羅立民,舒華忠,等.放射治療計(jì)劃的優(yōu)化方法[J].中國(guó)醫(yī)療器械雜志,2007,31(6): 391-394.

    [11]CAO Guogang,LUO Limin,RONG Chengcheng.Multicore-based parallelized differential evolution for medical image registration[C].MIPPR2009,2009,7497: 74972L1-74972L7.

    [12]曹國(guó)剛.計(jì)算機(jī)輔助手術(shù)導(dǎo)航關(guān)鍵技術(shù)研究[D].東南大學(xué)博士論文,2010.

    猜你喜歡
    計(jì)算速度錐形線程
    下頜管在下頜骨內(nèi)解剖結(jié)構(gòu)的錐形束CT測(cè)量
    淺談小學(xué)數(shù)學(xué)教學(xué)中學(xué)生計(jì)算能力的培養(yǎng)與提高
    學(xué)周刊(2019年3期)2019-01-11 11:33:48
    小學(xué)生數(shù)學(xué)思維能力培養(yǎng)的幾種策略
    淺談linux多線程協(xié)作
    錐形束CT結(jié)合顯微超聲技術(shù)診治老年鈣化根管的應(yīng)用
    宮頸錐形切除術(shù)后再次妊娠分娩方式的探討
    錐形流量計(jì)尾流流場(chǎng)分析
    美國(guó)將造超級(jí)計(jì)算機(jī)之王?速度超天河二號(hào)30倍
    探析小學(xué)數(shù)學(xué)教學(xué)中如何提升學(xué)生的計(jì)算能力
    Linux線程實(shí)現(xiàn)技術(shù)研究
    久久99热这里只频精品6学生| 性色av一级| 亚洲高清免费不卡视频| 91精品伊人久久大香线蕉| 99久久中文字幕三级久久日本| 亚洲美女视频黄频| 国产精品成人在线| 香蕉国产在线看| 亚洲av.av天堂| 一区二区av电影网| 色吧在线观看| 中文字幕最新亚洲高清| 男男h啪啪无遮挡| 人妻少妇偷人精品九色| 免费看av在线观看网站| 久久久国产欧美日韩av| 午夜影院在线不卡| 久久久国产欧美日韩av| 尾随美女入室| 欧美日韩综合久久久久久| 天天操日日干夜夜撸| 香蕉丝袜av| 美女大奶头黄色视频| 亚洲av成人精品一二三区| 极品人妻少妇av视频| 国产麻豆69| av片东京热男人的天堂| 一级片'在线观看视频| 国产亚洲欧美精品永久| 亚洲情色 制服丝袜| 国产亚洲欧美精品永久| 丰满迷人的少妇在线观看| 91精品伊人久久大香线蕉| 91精品国产国语对白视频| 男男h啪啪无遮挡| 国内精品宾馆在线| 男人爽女人下面视频在线观看| 三级国产精品片| 一级毛片黄色毛片免费观看视频| 久久久久国产网址| 精品一品国产午夜福利视频| 女人被躁到高潮嗷嗷叫费观| 99热这里只有是精品在线观看| 国产老妇伦熟女老妇高清| 国产免费福利视频在线观看| 满18在线观看网站| 亚洲av国产av综合av卡| 亚洲精品乱久久久久久| 激情五月婷婷亚洲| 精品少妇久久久久久888优播| 日韩成人av中文字幕在线观看| xxxhd国产人妻xxx| 黄色怎么调成土黄色| 国产精品一区二区在线不卡| 国国产精品蜜臀av免费| 男女边吃奶边做爰视频| 国产精品欧美亚洲77777| 最近最新中文字幕免费大全7| 亚洲av成人精品一二三区| av又黄又爽大尺度在线免费看| 在线观看国产h片| av卡一久久| 中文乱码字字幕精品一区二区三区| 久热久热在线精品观看| 国产男女超爽视频在线观看| 久久人人爽人人爽人人片va| 在线 av 中文字幕| 宅男免费午夜| 国产精品无大码| 亚洲内射少妇av| 国产视频首页在线观看| 国产高清不卡午夜福利| 久久av网站| 色94色欧美一区二区| 欧美xxxx性猛交bbbb| 人人澡人人妻人| av一本久久久久| 免费久久久久久久精品成人欧美视频 | 老司机影院成人| 中文字幕亚洲精品专区| 成年av动漫网址| 亚洲美女搞黄在线观看| 久久精品熟女亚洲av麻豆精品| 久久av网站| 国产av码专区亚洲av| av在线观看视频网站免费| 各种免费的搞黄视频| 国精品久久久久久国模美| a级毛色黄片| 国产精品人妻久久久久久| 欧美精品人与动牲交sv欧美| 我的女老师完整版在线观看| av女优亚洲男人天堂| 成人午夜精彩视频在线观看| 国产精品一国产av| 99九九在线精品视频| 涩涩av久久男人的天堂| 22中文网久久字幕| 精品少妇久久久久久888优播| 久久久久久久国产电影| 永久网站在线| 久久综合国产亚洲精品| 秋霞在线观看毛片| 精品国产一区二区三区四区第35| 国产成人精品在线电影| 午夜视频国产福利| 久久亚洲国产成人精品v| 人妻一区二区av| 成年美女黄网站色视频大全免费| 欧美xxⅹ黑人| 久久精品国产自在天天线| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品久久午夜乱码| 亚洲精品日本国产第一区| 乱人伦中国视频| 高清视频免费观看一区二区| 久久免费观看电影| 日韩成人av中文字幕在线观看| 精品99又大又爽又粗少妇毛片| 黑丝袜美女国产一区| 黑人高潮一二区| 97人妻天天添夜夜摸| 人妻 亚洲 视频| 边亲边吃奶的免费视频| 成人国产麻豆网| 亚洲精品色激情综合| 国产片内射在线| 精品久久国产蜜桃| 一级毛片 在线播放| 国产精品欧美亚洲77777| 中文字幕av电影在线播放| 欧美国产精品va在线观看不卡| 在线免费观看不下载黄p国产| 亚洲图色成人| 久久av网站| a级毛片在线看网站| 黑人高潮一二区| 一本久久精品| 蜜桃在线观看..| av免费观看日本| 日韩欧美精品免费久久| 国产高清三级在线| 亚洲经典国产精华液单| 国产免费一级a男人的天堂| 成人毛片60女人毛片免费| 99国产综合亚洲精品| 国产免费一级a男人的天堂| 26uuu在线亚洲综合色| 欧美精品国产亚洲| 亚洲精品久久午夜乱码| 亚洲精品国产av成人精品| 三级国产精品片| 欧美97在线视频| 香蕉丝袜av| www.色视频.com| 久久99热这里只频精品6学生| 麻豆乱淫一区二区| 精品视频人人做人人爽| 一级毛片黄色毛片免费观看视频| 日韩精品有码人妻一区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av电影在线观看一区二区三区| 在线 av 中文字幕| 99热6这里只有精品| 九九爱精品视频在线观看| 亚洲国产精品999| 熟女电影av网| 99香蕉大伊视频| 有码 亚洲区| 丝袜美足系列| 如何舔出高潮| 欧美激情 高清一区二区三区| 精品人妻偷拍中文字幕| 18禁裸乳无遮挡动漫免费视频| 国产成人a∨麻豆精品| 欧美日韩亚洲高清精品| a级片在线免费高清观看视频| 久久久久久久久久久久大奶| 精品一区二区三区四区五区乱码 | 高清视频免费观看一区二区| 国产成人精品无人区| 少妇人妻 视频| 欧美另类一区| 美女内射精品一级片tv| 亚洲精品国产av蜜桃| 国产毛片在线视频| 蜜臀久久99精品久久宅男| 亚洲av电影在线观看一区二区三区| 国产爽快片一区二区三区| 国产成人精品久久久久久| av天堂久久9| 高清毛片免费看| 母亲3免费完整高清在线观看 | 国产成人a∨麻豆精品| 亚洲三级黄色毛片| 汤姆久久久久久久影院中文字幕| 精品一区在线观看国产| 一级毛片我不卡| 丝袜脚勾引网站| a级毛色黄片| 成人综合一区亚洲| 捣出白浆h1v1| 成人黄色视频免费在线看| 亚洲色图综合在线观看| 黄色视频在线播放观看不卡| 99久久人妻综合| 看免费成人av毛片| 69精品国产乱码久久久| 亚洲色图综合在线观看| 亚洲国产成人一精品久久久| 女人久久www免费人成看片| 日本vs欧美在线观看视频| 美国免费a级毛片| 国产在线一区二区三区精| 美女中出高潮动态图| 亚洲,欧美精品.| 免费观看av网站的网址| 免费在线观看黄色视频的| 在线观看国产h片| 国产免费一区二区三区四区乱码| 久久久久精品性色| 一级爰片在线观看| 人人妻人人爽人人添夜夜欢视频| 全区人妻精品视频| 亚洲天堂av无毛| 日本av手机在线免费观看| 亚洲国产精品成人久久小说| 日韩视频在线欧美| 日韩中字成人| 国产日韩欧美视频二区| 久久国产精品大桥未久av| 狂野欧美激情性xxxx在线观看| 国产深夜福利视频在线观看| 一区在线观看完整版| 成年人午夜在线观看视频| 在线天堂中文资源库| 一边亲一边摸免费视频| 国产乱人偷精品视频| 18在线观看网站| 最近手机中文字幕大全| 免费大片黄手机在线观看| 色哟哟·www| 久久精品国产亚洲av天美| 满18在线观看网站| 黑人猛操日本美女一级片| 下体分泌物呈黄色| 久久久国产欧美日韩av| 午夜视频国产福利| 国产福利在线免费观看视频| 桃花免费在线播放| 热re99久久国产66热| 欧美成人精品欧美一级黄| 日本黄色日本黄色录像| 久久99精品国语久久久| 亚洲情色 制服丝袜| 欧美成人精品欧美一级黄| 国产免费现黄频在线看| 精品一区二区免费观看| 在线 av 中文字幕| 久久国产精品大桥未久av| 在线观看美女被高潮喷水网站| 久久久a久久爽久久v久久| 亚洲天堂av无毛| 一二三四在线观看免费中文在 | 黑人高潮一二区| 精品99又大又爽又粗少妇毛片| 2022亚洲国产成人精品| 大香蕉97超碰在线| 日韩一区二区视频免费看| 只有这里有精品99| 视频中文字幕在线观看| 亚洲精品第二区| 丝瓜视频免费看黄片| 高清视频免费观看一区二区| 这个男人来自地球电影免费观看 | 国产亚洲精品第一综合不卡 | 最后的刺客免费高清国语| 午夜影院在线不卡| 日韩人妻精品一区2区三区| 久久热在线av| 天天躁夜夜躁狠狠躁躁| 男女边摸边吃奶| 日韩免费高清中文字幕av| www.av在线官网国产| 在线观看国产h片| 又黄又爽又刺激的免费视频.| 国产熟女午夜一区二区三区| 少妇 在线观看| 亚洲欧美一区二区三区黑人 | av卡一久久| 男人添女人高潮全过程视频| 免费观看在线日韩| 国产国拍精品亚洲av在线观看| 日韩一区二区三区影片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲欧美日韩另类电影网站| 国产免费又黄又爽又色| 国产免费一级a男人的天堂| 国产片内射在线| 亚洲成人一二三区av| 国产成人一区二区在线| 亚洲色图 男人天堂 中文字幕 | 三上悠亚av全集在线观看| 九九在线视频观看精品| 国产精品免费大片| 啦啦啦啦在线视频资源| 亚洲人成77777在线视频| 国产精品三级大全| 晚上一个人看的免费电影| 日韩,欧美,国产一区二区三区| 精品一区二区三区四区五区乱码 | 国产 一区精品| 国产综合精华液| 18禁动态无遮挡网站| 久久久欧美国产精品| 亚洲一区二区三区欧美精品| 免费观看性生交大片5| 十八禁高潮呻吟视频| 天天躁夜夜躁狠狠久久av| 亚洲欧美日韩卡通动漫| 91成人精品电影| 久久99精品国语久久久| 国产精品久久久久久av不卡| 中文字幕av电影在线播放| 下体分泌物呈黄色| 欧美人与善性xxx| 精品久久久久久电影网| 2021少妇久久久久久久久久久| 久久婷婷青草| 精品久久久久久电影网| 国产精品一二三区在线看| 亚洲成人av在线免费| 中文精品一卡2卡3卡4更新| 九九爱精品视频在线观看| 免费高清在线观看视频在线观看| 大片免费播放器 马上看| 另类亚洲欧美激情| 亚洲国产精品专区欧美| 国产精品久久久久久精品电影小说| 90打野战视频偷拍视频| 亚洲欧洲日产国产| 伊人久久国产一区二区| av在线观看视频网站免费| 大片免费播放器 马上看| 另类亚洲欧美激情| 亚洲成国产人片在线观看| 国产精品一区www在线观看| 欧美日韩综合久久久久久| 天美传媒精品一区二区| 青春草视频在线免费观看| 这个男人来自地球电影免费观看 | 亚洲,欧美,日韩| 香蕉丝袜av| 亚洲国产精品专区欧美| 狂野欧美激情性bbbbbb| 成人免费观看视频高清| av天堂久久9| 国产精品久久久av美女十八| 国产国语露脸激情在线看| a级毛色黄片| 99国产综合亚洲精品| 9色porny在线观看| 18禁观看日本| av网站免费在线观看视频| 观看美女的网站| 纯流量卡能插随身wifi吗| 日韩一本色道免费dvd| 丰满少妇做爰视频| 国产精品不卡视频一区二区| 亚洲伊人久久精品综合| 亚洲欧美日韩另类电影网站| 99re6热这里在线精品视频| 日韩中文字幕视频在线看片| 蜜桃国产av成人99| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕最新亚洲高清| 97在线人人人人妻| 香蕉国产在线看| 美国免费a级毛片| 国产伦理片在线播放av一区| 日韩制服骚丝袜av| 日韩 亚洲 欧美在线| 一级毛片黄色毛片免费观看视频| 欧美少妇被猛烈插入视频| av黄色大香蕉| 最黄视频免费看| 国产精品国产三级国产专区5o| 最近最新中文字幕大全免费视频 | 另类精品久久| 99精国产麻豆久久婷婷| 欧美性感艳星| 又黄又爽又刺激的免费视频.| av在线观看视频网站免费| 精品国产国语对白av| 亚洲精华国产精华液的使用体验| 国产精品国产三级国产av玫瑰| 成年人午夜在线观看视频| 久久久国产精品麻豆| 欧美国产精品一级二级三级| 99久久人妻综合| 只有这里有精品99| 欧美国产精品va在线观看不卡| 岛国毛片在线播放| 男女高潮啪啪啪动态图| 国产亚洲av片在线观看秒播厂| 男人操女人黄网站| 下体分泌物呈黄色| 久久久久久久久久久免费av| 激情五月婷婷亚洲| 亚洲图色成人| 亚洲av国产av综合av卡| 欧美日韩国产mv在线观看视频| 亚洲成av片中文字幕在线观看 | 精品熟女少妇av免费看| 亚洲 欧美一区二区三区| 日韩一区二区视频免费看| 十分钟在线观看高清视频www| 欧美 日韩 精品 国产| 一级,二级,三级黄色视频| 国产精品久久久久成人av| 国产成人精品福利久久| 色视频在线一区二区三区| 午夜免费男女啪啪视频观看| 黑人高潮一二区| 黄色 视频免费看| 蜜臀久久99精品久久宅男| a 毛片基地| 日本-黄色视频高清免费观看| 国精品久久久久久国模美| 精品视频人人做人人爽| 一区在线观看完整版| 曰老女人黄片| av.在线天堂| 久久99热这里只频精品6学生| 国产午夜精品一二区理论片| 精品少妇久久久久久888优播| 久久久国产精品麻豆| 午夜福利影视在线免费观看| 天天躁夜夜躁狠狠久久av| 久久精品久久久久久久性| 免费在线观看黄色视频的| 涩涩av久久男人的天堂| 国国产精品蜜臀av免费| 国产探花极品一区二区| 国精品久久久久久国模美| 最近中文字幕2019免费版| 欧美国产精品va在线观看不卡| 日韩av不卡免费在线播放| 久久久精品区二区三区| 亚洲国产av影院在线观看| 亚洲第一av免费看| 久久久久人妻精品一区果冻| 女性生殖器流出的白浆| 欧美国产精品va在线观看不卡| 亚洲av中文av极速乱| 亚洲国产欧美日韩在线播放| 中国三级夫妇交换| 亚洲人成网站在线观看播放| 黄色 视频免费看| 最新的欧美精品一区二区| 亚洲精品456在线播放app| 少妇被粗大猛烈的视频| 国产色爽女视频免费观看| 免费av中文字幕在线| 91精品三级在线观看| 18禁在线无遮挡免费观看视频| 成年美女黄网站色视频大全免费| 人妻系列 视频| 26uuu在线亚洲综合色| videos熟女内射| 精品久久国产蜜桃| 久久鲁丝午夜福利片| 久久精品国产综合久久久 | 久久久精品区二区三区| www.色视频.com| 日韩电影二区| 亚洲精品久久午夜乱码| 中文字幕人妻丝袜制服| 99热6这里只有精品| 少妇被粗大猛烈的视频| 国产免费又黄又爽又色| 视频中文字幕在线观看| 黑丝袜美女国产一区| 黄色视频在线播放观看不卡| 中文字幕人妻熟女乱码| 夜夜爽夜夜爽视频| 午夜免费鲁丝| 你懂的网址亚洲精品在线观看| 97精品久久久久久久久久精品| 少妇被粗大猛烈的视频| 国产精品久久久av美女十八| 欧美性感艳星| 精品国产露脸久久av麻豆| 少妇 在线观看| 久久久久精品久久久久真实原创| 青春草视频在线免费观看| 综合色丁香网| 亚洲第一区二区三区不卡| 亚洲精品成人av观看孕妇| 免费播放大片免费观看视频在线观看| 国产色婷婷99| 看免费av毛片| 中文字幕人妻丝袜制服| 丰满饥渴人妻一区二区三| 日韩制服骚丝袜av| 赤兔流量卡办理| 国产免费视频播放在线视频| 欧美精品一区二区大全| 中文字幕制服av| 啦啦啦在线观看免费高清www| 免费在线观看完整版高清| 欧美日韩国产mv在线观看视频| 男女啪啪激烈高潮av片| 丝瓜视频免费看黄片| 夫妻午夜视频| 亚洲成色77777| 丝袜脚勾引网站| 观看av在线不卡| 国产精品久久久av美女十八| 色视频在线一区二区三区| 你懂的网址亚洲精品在线观看| 国产日韩欧美视频二区| 精品国产一区二区久久| 欧美xxⅹ黑人| 春色校园在线视频观看| 亚洲av日韩在线播放| 大话2 男鬼变身卡| 中国国产av一级| 午夜福利,免费看| √禁漫天堂资源中文www| 午夜影院在线不卡| 精品午夜福利在线看| 精品一区二区三卡| 日本猛色少妇xxxxx猛交久久| h视频一区二区三区| 又粗又硬又长又爽又黄的视频| 高清在线视频一区二区三区| 日韩精品有码人妻一区| 自线自在国产av| 亚洲国产av影院在线观看| 国产深夜福利视频在线观看| 精品一区二区免费观看| 精品亚洲乱码少妇综合久久| 国产成人91sexporn| 十八禁高潮呻吟视频| 桃花免费在线播放| 午夜视频国产福利| 嫩草影院入口| 蜜臀久久99精品久久宅男| 精品人妻一区二区三区麻豆| 久久热在线av| 99视频精品全部免费 在线| 国产亚洲一区二区精品| 18禁裸乳无遮挡动漫免费视频| a级片在线免费高清观看视频| 丰满少妇做爰视频| 超色免费av| 国产亚洲av片在线观看秒播厂| 青青草视频在线视频观看| 狠狠婷婷综合久久久久久88av| 妹子高潮喷水视频| 国产成人精品婷婷| 国产欧美亚洲国产| 欧美精品亚洲一区二区| 美女国产视频在线观看| 高清视频免费观看一区二区| 国产老妇伦熟女老妇高清| 又大又黄又爽视频免费| 赤兔流量卡办理| 一区二区三区乱码不卡18| 在线观看免费日韩欧美大片| 亚洲第一av免费看| 国精品久久久久久国模美| 一级片免费观看大全| 精品99又大又爽又粗少妇毛片| 丰满乱子伦码专区| 久久国产精品大桥未久av| 国产亚洲一区二区精品| 中文精品一卡2卡3卡4更新| 亚洲精品美女久久av网站| 熟女人妻精品中文字幕| 久久久久久伊人网av| 在线精品无人区一区二区三| 高清欧美精品videossex| 中文字幕人妻熟女乱码| 另类亚洲欧美激情| 高清欧美精品videossex| 国产一区二区在线观看日韩| 午夜免费男女啪啪视频观看| 中文字幕精品免费在线观看视频 | 人妻系列 视频| 亚洲精品色激情综合| 国产成人精品福利久久| 亚洲精品成人av观看孕妇| 一本色道久久久久久精品综合| 亚洲国产日韩一区二区| 欧美xxxx性猛交bbbb| 色吧在线观看| 看免费成人av毛片| 久久精品夜色国产| 色吧在线观看| 久久鲁丝午夜福利片| 久久精品夜色国产| 18禁裸乳无遮挡动漫免费视频| 久久久久国产网址| 少妇熟女欧美另类| 日韩中文字幕视频在线看片| 亚洲中文av在线| 人妻系列 视频| 久久人人爽人人爽人人片va| 免费高清在线观看日韩| 国产精品久久久久久精品古装| 国产精品三级大全| 日本免费在线观看一区| 三上悠亚av全集在线观看|