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

    光中子在現(xiàn)代放射治療領(lǐng)域的研究進(jìn)展*

    2017-01-09 09:12:34陳亞正許敬輝李廚榮廖雄飛
    中國醫(yī)學(xué)裝備 2016年12期
    關(guān)鍵詞:電離室中子加速器

    陳亞正 徐 珂 許敬輝 李廚榮 廖雄飛 黎 杰 王 培*

    光中子在現(xiàn)代放射治療領(lǐng)域的研究進(jìn)展*

    陳亞正①徐 珂①許敬輝①李廚榮①廖雄飛①黎 杰①王 培①*

    隨著各種醫(yī)用加速器新技術(shù)的提高和腫瘤患者數(shù)量的急劇增加,放射治療已經(jīng)成為治療腫瘤的主要手段之一。醫(yī)用加速器技術(shù)的不斷進(jìn)步,始終圍繞著照射腫瘤的同時最大限度地保護(hù)正常器官而發(fā)展,但仍無法解決由于高能射線在治療患者時所帶來的輻射污染。高能粒子產(chǎn)生的次級粒子如光子、光中子及質(zhì)子等對人體組織產(chǎn)生的生物效應(yīng)同樣不可忽略,并逐漸受到重視和廣泛的研究。由于加速器機頭內(nèi)部構(gòu)造的金屬元素為鉛,靶區(qū)主要為鎢,光子能量達(dá)到>7 MeV時,會超過光核反應(yīng)閾值從而產(chǎn)生中子,這一反應(yīng)稱為光核反應(yīng)。國外放療機構(gòu)對光中子產(chǎn)生的劑量開展了廣泛研究?;趨⑴c復(fù)雜相互作用的粒子特點,通過蒙特卡羅模擬和實驗測量兩個途徑來研究醫(yī)用加速器產(chǎn)生的光中子通量和能譜等物理量。由于醫(yī)用加速器設(shè)計不同(主要是機頭內(nèi)部),性能也有很大差別,因此光中子的通量、能譜和劑量也有不同的分布規(guī)律。在文獻(xiàn)評閱的基礎(chǔ)上,從蒙特卡羅模擬和實驗測量兩個方面闡述目前對醫(yī)用加速器產(chǎn)生光中子的研究現(xiàn)狀以及放射治療過程中的輻射防護(hù)和引發(fā)第二原發(fā)癌的物理因素。

    光中子;光核反應(yīng);劑量;通量;醫(yī)用直線加速器

    [First-author’s address]Radiotherapy Center, Sichuan Cancer Hospital, Chengdu 610041, China.

    放射治療是目前腫瘤治療中最重要的手段,而醫(yī)用直線加速器則是放射治療的主要工具。隨著計算機和機電一體化技術(shù)的進(jìn)步,醫(yī)用直線加速器近年有了極大的發(fā)展,隨之帶來的放射治療技術(shù)的革命更是日新月異。以調(diào)強放射治療(intensity modulated radiation therapy,IMRT),容積旋轉(zhuǎn)調(diào)強放射治療(volumetric modulated arc therapy,VMAT),螺旋斷層調(diào)強放射治療(helical tomotherapy,TOMO)為代表的現(xiàn)代放射治療技術(shù)巳經(jīng)在國內(nèi)外廣泛開展,每年數(shù)以百萬計的癌癥患者在接受這些設(shè)備和技術(shù)的治療。由于醫(yī)用直線加速器是一種用高能光子和電子來殺滅癌瘤細(xì)胞的設(shè)備。在高能光子和電子產(chǎn)生及輸運過程中要接觸不同靶材料、空氣、準(zhǔn)直器等。當(dāng)光子能量超過7MeV時會通過光核反應(yīng)(γ,n)或者電核反應(yīng)(e,e’n)產(chǎn)生中子,導(dǎo)致治療過程中存在一定程度的中子污染危害。電子與核子相互碰撞(電核反應(yīng))的散射截面比光子與核子碰撞(光核反應(yīng))的截面小2個數(shù)量級,因此電核反應(yīng)產(chǎn)生的中子可以不予考慮[1-2]。由于癌癥患者數(shù)量的急劇增加以及腫瘤新病種的出現(xiàn),醫(yī)院對于醫(yī)用加速器的需求不斷提高,治療機房內(nèi)產(chǎn)生的中子劑量貢獻(xiàn)也會越來越大,無論是從輻射防護(hù)的角度還是從次生粒子引起的新癌變風(fēng)險來考慮,都有必要研究其物理特性和風(fēng)險評估。通過對醫(yī)用加速器產(chǎn)生的光中子研究現(xiàn)狀進(jìn)行回顧分析,為中子防護(hù)輻射措施的規(guī)范化和第二原發(fā)癌風(fēng)險評估提供可參考性資料。

    1 醫(yī)用直線加速器產(chǎn)生的光中子

    中子是中性粒子,不帶電荷,國際輻射測量與測量委員會(International Commission on Radiological Protection,ICRP)最新報告指出中子的品質(zhì)因子約為20,即生物效應(yīng)遠(yuǎn)遠(yuǎn)大于光子的生物效應(yīng)[3]。當(dāng)中子達(dá)到一定能量時會對人體產(chǎn)生生物效應(yīng),對患者及放射治療工作人員產(chǎn)生潛在的輻射危險和健康隱患[4-5]。有研究表明,對于8~10 MV的電子束,出射線18 000 MU的跳數(shù)后在等中心處產(chǎn)生的光中子通量已經(jīng)超過了108cm2。隨著環(huán)境污染的加劇,飲食的不安全性,以及生活習(xí)慣的不規(guī)律,導(dǎo)致腫瘤患者的數(shù)量近年來不斷增加,國家衛(wèi)生計生委的統(tǒng)計資料表明,目前我國每年新生腫瘤患者總數(shù)約為250萬人,其中有>50%為惡性腫瘤新生患者,因癌癥死亡人數(shù)為140萬。我國每年新增腫瘤患者呈逐漸增長趨勢,年均增長率在10%左右。各省市及縣級單位醫(yī)院都相繼購進(jìn)醫(yī)療加速器用于放射治療,但是往往忽略了中子污染所帶來的影響。

    1.1 光中子污染

    研究電子直線加速器所產(chǎn)生的光中子污染已成為當(dāng)務(wù)之急。國家食品藥品監(jiān)督管理總局最新提示“關(guān)注醫(yī)用電子直線加速器風(fēng)險”,據(jù)2013年統(tǒng)計,自2010年1月1日至2013年6月20日,國家藥品不良反應(yīng)監(jiān)測中心共收到涉及醫(yī)用電子直線加速器的可疑不良事件報告223份[6];其中嚴(yán)重傷害報告99份,主要為放射治療引發(fā)的并發(fā)癥,表現(xiàn)為白細(xì)胞減少、血小板減少、脫發(fā)、放射性皮炎(皮膚潰瘍、紅腫及色素沉著等)、非照射部位炎癥(口腔黏膜炎、食管黏膜炎等)、惡心、嘔吐、厭食、腹脹和腹瀉等;其中白細(xì)胞減少53例,占53.5%;皮膚潰瘍20例,占20.2%;皮膚紅腫8例,占8.1%;惡心和嘔吐11例,占11.1%;包括了光中子污染帶來的危害。

    1.2 光中子污染與第二腫瘤

    光中子污染帶來的健康隱患之一是有可能引發(fā)第二腫瘤。惡性腫瘤采用放射治療后,腫瘤細(xì)胞DNA被X射線殺死,與此同時,人體正常的細(xì)胞也暴露于X光以及產(chǎn)生的光中子損傷下。這些射線足以導(dǎo)致細(xì)胞誘變和畸變,如果在較長的時間內(nèi)不能得到恢復(fù),第二惡性腫瘤就很有可能發(fā)生。腫瘤放射治療后誘發(fā)二次腫瘤的標(biāo)準(zhǔn)包括3個方面:①患者有過放射治療的經(jīng)歷;②誘發(fā)腫瘤需要有較長的潛伏期,一般在放射治療后4~5年;③需要病理組織學(xué)組織證實第二腫瘤與原發(fā)惡性腫瘤有差異,排除原發(fā)灶轉(zhuǎn)移或者復(fù)發(fā)的可能性。此外,放射治療技術(shù)的更新?lián)Q代也是誘發(fā)第二腫瘤的因素之一。如從三維適形放射治療(threedimensional conformal radiotherapy,3D-CRT)到IMRT的轉(zhuǎn)變,由于3D-CRT是從三個方向進(jìn)行固定野出射線,未被照到的正常組織能得到很好的保護(hù);而IMRT技術(shù)在出射線時會有多野光闌(multi leaf colimator,MLC)的跟隨變動,會出現(xiàn)更多角度的射野,同時機器要消耗更多的能量來增加跳數(shù)以滿足調(diào)強對腫瘤處方劑量的需要,更多的正常組織甚至距離腫瘤較遠(yuǎn)的器官也會受到照射。不同物理師對同一計劃的射野角度有不同的考慮,因此腫瘤以外的人體組織受到低劑量沉積的范圍更大。由于這些原因?qū)е律嫦聛淼姆派渲委熀蠡颊哒T發(fā)第二次腫瘤的概率比3D-CRT的情況高出0.5%[7]。

    相比非調(diào)強野,調(diào)強野為等中心點傳輸一個特定的劑量時,不管是動態(tài)IMRT還是靜態(tài)調(diào)強方法,均會要求機器激發(fā)更多的跳數(shù),因此落在患者身上的劑量必然增加。這種劑量累計有多種因素:即多葉光柵的設(shè)計及運動形式、治療計劃調(diào)強的方式以及調(diào)強過程中每個射野的MLC形成一系列的形狀。IMRT計劃會使跳數(shù)因子增加2~3,這種劑量泄露及散射輻射遠(yuǎn)遠(yuǎn)增加了腫瘤邊界以外的劑量。國外有研究[8]用電離室對模體進(jìn)行了一項測量,使用6 MV加速器對宮頸進(jìn)行4野照射70 Gy,乳腺受量為0.25 Gy。如果用20 MV的電子直線加速器進(jìn)行照射,相同的測量點受量是0.5 Gy的X射線和1 cGy的光中子。一般臨床治療患者不適用20 MV的能量檔。乳腺0.25 Gy的劑量換算為第二腫瘤危險概率是0.25%。IMRT技術(shù)所導(dǎo)致的第二腫瘤誘發(fā)率比3D-CRT高出0.5%[9]。光中子作為醫(yī)用加速器輻射源之一,對第二腫瘤的誘發(fā)有著極為關(guān)鍵的影響,因此有必要獲取大量詳細(xì)可靠的臨床實際數(shù)據(jù)來進(jìn)行研究。

    2 蒙特卡羅模擬光中子的產(chǎn)生

    加速器機頭產(chǎn)生光中子的來源包括靶、一級準(zhǔn)直器、二級準(zhǔn)直器、楔形板以及多葉光柵。鎢元素和鉛元素與光子發(fā)生光核反應(yīng)的散射截面最大,產(chǎn)生光中子的貢獻(xiàn)最大。鐵、銅和鋁雖然也是機頭的構(gòu)成元素,但是產(chǎn)生光中子的貢獻(xiàn)遠(yuǎn)遠(yuǎn)小于前者。鎢和鉛的光核反應(yīng)閾值分別是6.74 MeV和6.19 MeV,而銅和鐵的反應(yīng)閾值是9.91 MeV和7.65 MeV。從光子能量的角度而言,光中子相互作用概率隨著光子能量增加而急劇上升,并且在13~18 MeV的范圍內(nèi)會達(dá)到最大值[9]。同時,鎢和鉛對中子的散射截面很小,中子不被機頭吸收,有很大概率穿過機頭內(nèi)部的各種結(jié)構(gòu),到達(dá)患者身上和墻壁。

    近年有眾多研究使用不同劑量儀器測量的中子劑量,然而實驗方法無法追蹤中子產(chǎn)生的來源,解決這個問題的辦法就是借助蒙特卡羅程序。蒙特卡羅程序通過建立加速器機頭模型和模擬電子轟擊靶的相互作用過程,來追蹤電子、光子和中子從產(chǎn)生到消失的輸運過程,以及各種離子在不同位置所沉積的能量。有眾多研究報道,用蒙特卡羅模擬光中子的產(chǎn)生來分析加速器產(chǎn)生的光中子對人體以及機房室內(nèi)的劑量貢獻(xiàn)。Ghiasi等[10]用蒙卡程序模擬了考慮機頭所有部件及簡單部件兩者情況下產(chǎn)生的光中子產(chǎn)額,結(jié)果表明機頭內(nèi)部構(gòu)造如準(zhǔn)直器、多葉光柵對光子散射和中子的形成有重要影響。Chibani等[11]用蒙特卡羅中子光子(Monte Carlo neutron photon,MCNP)程序模擬了仿真人體在Varian 和Simens加速器15 MV、18 MV的光核反應(yīng)劑量,產(chǎn)生的粒子包括中子、質(zhì)子和α粒子。同時,研究了射野大小的影響,離軸曲線上的劑量分布,機頭各組成部分產(chǎn)生中子的貢獻(xiàn),以及仿真人體的不均勻性效應(yīng)。結(jié)果表明,對于18 MV的Siemens加速器(10 cm×10 cm)所有粒子的劑量當(dāng)量率(粒子最大劑量與光子最大劑量的比值)0.66 cSv/ Gy,中子當(dāng)量率隨著射野減小而減小,但是質(zhì)子和α粒子的當(dāng)量率不受射野大小的影響。Varian加速器(15 MV、18 MV)產(chǎn)生的各種粒子(中子、質(zhì)子和α粒子)多于Siemens產(chǎn)生的粒子,其主要原因是Varian的原初電子能量大于Siemens的電子能量。對于所有的射線束,產(chǎn)生的中子劑量占所有劑量的75%。對于15 MV和18 MV的Varian機器(10 cm×10 cm),中子當(dāng)量率分別為1.52 cSv/Gy和2.86 cSv/Gy。Ghiasi等[10]使用蒙特卡羅程序MCNP 2.4.0對加速器機頭進(jìn)行了完整模型和簡化模型的建立,計算了在10 cm×10 cm下中子通量的計算,認(rèn)為兩個模型計算出來的中子通量有4%~53%的差別。在完整模型里面,中子通量隨著射野面積增加而增加,然而在簡化模型里面,通量趨勢與前者完全相反。Zanini等[12]模擬了18 MV的加速器不同機頭產(chǎn)生的光中子在患者身上沉積的劑量,機頭構(gòu)造不一樣,中子通量也有所差別;距離等中心點越遠(yuǎn),中子劑量越小。Pena等[13]計算了Primus加速器15 MV的光子在10 cm×10 cm射野下機頭各組成部分產(chǎn)生的光中子,其中一級準(zhǔn)直器產(chǎn)生的中子占52%,二級準(zhǔn)直器產(chǎn)生的中子占21%,鎢靶占12%,多葉光柵占6.6%,機頭外殼占5%,均整器占0.41%。Mao等[14]對Varian 2100C/2300C加速器進(jìn)行了不同光子能量的模擬,得出的結(jié)果與前面類似,但是有稍許偏差。一級準(zhǔn)直器產(chǎn)生的中子最多,其次是鉛門、靶、均整器以及其他部件(見表1)。

    表1 蒙特卡羅模擬Varian加速器機頭各組成部分產(chǎn)生的光中子

    Becker等[15]做了類似的工作,在其蒙特卡羅模擬實驗結(jié)果中,準(zhǔn)直器、MLC和鉛門、靶(由鎢元素構(gòu)成)和均整器各產(chǎn)生的中子貢獻(xiàn)分別為54.85%、26.72%、10.08%和1.74%,與Pena等[13]的結(jié)果比較接近。Martínez-Ovalle等[16]使用MCNPX(V.2.5)計算了4臺不同的加速器(Varian Clinac2100C、Elekta Inor、Elekta SL25和Simens Mevatron KDS)在電子能量為15 MV、18 MV和20 MV時產(chǎn)生的中子劑量當(dāng)量和能譜,位置在射線軸延長線的模體表面和模體表面以下1 cm處以及射野外不同的深度。在相同額定能量下,Varian機器產(chǎn)生的中子劑量遠(yuǎn)遠(yuǎn)多于其他機器產(chǎn)生的中子劑量。Barquero等[17]使用蒙特卡羅程序模擬18 MV的醫(yī)用直線加速器產(chǎn)生光中子在仿真人體重要器官所沉積的能量,該研究從4個不同角度的射野(0o、90o、180o和270o)給出射線,計算射線方向距離皮膚最近的重要器官的中子劑量,射野大小為18 cm×18 cm。模擬結(jié)果顯示在0o方向,小腸所受光中子劑量為190 μSv/Gy,90o方向結(jié)腸所受光中子劑量為51 μSv;180o方向直腸所受光中子劑量為719 μSv/Gy;270o方向皮膚沉積的中子劑量為45 μSv。同時射野內(nèi)的中子劑量遠(yuǎn)遠(yuǎn)大于射野以外區(qū)域的中子劑量。

    大部分研究者使用的蒙特卡羅模擬程序為MCNP(版本有所不同),建立的機頭模型有所差別,使用的散射截面數(shù)據(jù)庫和程序參數(shù)也有所不同,因此計算結(jié)果不可能完全相同,但是相差并不大。

    3 加速器產(chǎn)生光中子的實驗測量

    測量光中子常用的儀器主要有Bonner球譜儀(由多個慢化球和中心探測器組成)、電離室、LiF熱釋光片(TLD-600、TLD-700)、CR-39徑跡探測器以及由生產(chǎn)商制造的中子探測儀。按測量工具分別闡述國內(nèi)外對醫(yī)用電子直線加速器產(chǎn)生光中子的研究。

    曹磊等[18]利用中子氣泡探測器測量了15 MV電子直線加速器不同射野的中子輻射場,并且用MCNP程序進(jìn)行了簡單的模擬。范悅鋒等[19]使用銀活化系統(tǒng)測了12 MV的X光機中子注量;2003年和2006年黃薇玲等[20-21]采用被動光中子劑量測量方法,分別使用Unidose劑量儀、電離室、固體行跡CR-39片測了不同能量下的中子劑量當(dāng)量,同時使用MCNP模擬了中子能譜及其通量。隨著醫(yī)療加速器技術(shù)的不斷發(fā)展,尤其是IMRT、VMAT、TOMO技術(shù)的出現(xiàn),機器的束流劑量率已高達(dá)2400 cGY/分,子野數(shù)達(dá)數(shù)千個。加速器機頭的多葉光柵設(shè)計越來越多元化,過去的加速器機頭只有一級、二級準(zhǔn)直器,而現(xiàn)在的加速器加入了跟隨系統(tǒng)多葉光柵,在治療技術(shù)和效果得到明顯改善的同時,治療室內(nèi)的光中子污染特點也有了新的變化,有必要進(jìn)行重新測量和分析。Exposito等[22]使用中子數(shù)字設(shè)備(neutron digital device)測量了瓦里安2100和西門子Primus加速器的不同電子能量檔(6 MV、9 MV、12 MV和15 MV)及光子能量檔(6 MV、15 MV)的中子劑量。Krmar等[23]使用Meridian model 5085型中子探測儀測量了準(zhǔn)直器和MLC在開閉及不同射野下15 MV加速器產(chǎn)生的中子在迷路出口處的劑量以及距離等中心處100 cm的中子通量,結(jié)果顯示準(zhǔn)直器完全關(guān)閉時迷路出口處測得的光中子劑量最大,準(zhǔn)直器和多葉光柵完全開起時,光中子劑量最低,表明光子與準(zhǔn)直器及MLC碰撞產(chǎn)生的中子數(shù)量十分明顯。

    電離室也可以測量中子劑量,但是存在一定的局限性。不同元素的原子與中子相互作用散射截面大小有較大差異,導(dǎo)致電離室器件對中子的敏感度也不同。根據(jù)粒子相互作用的這一特點,電離室空腔內(nèi)壁涂上不同的材料,可以測量不同能量段的光子或者中子。等效組織電離室(tissue equivalent chamber,TE)、鎂電離室(magnesium chamber,Mg)、硼涂層電離室(boron coated magnesium chamber,MgB)這3套電離室可以相互組合。TE電離室對所有的劑量都敏感,包括光子和所有能量段的中子,因此測得的劑量是混合場,無法區(qū)別光子及各能量段的中子劑量。Mg電離室主要對光子敏感,對快中子的敏感度稍差,對熱中子的敏感度最小,但是為了能精確的測量劑量,仍然不能忽略熱中子的貢獻(xiàn)。MgB電離室對熱中子非常敏感,這是由于熱中子與MgB電離室的10B元素散射截面很大;而對光子和快中子的敏感度與Mg電離室效果一樣。根據(jù)文獻(xiàn)[24]提供的三維方程組分別計算出光子、快中子、熱中子的劑量。Koivunoro等[25]使用雙電離室(鎂電離室,硼涂層電離室)測量了水中射線軸上不同深度的熱中子劑量。但是由于電離室的特點,只能測得熱中子劑量,其他能量段的中子劑量需要3套電離室組合以及方程解析來求得。Burmeister等[26]利用了相似的電離室系統(tǒng)測量了快中子放射治療過程中產(chǎn)生的熱中子通量;Howell等[27]比較了對中子有不同敏感度的熱釋光器件(thermoluminescence devices,TLDs)和電離室系統(tǒng)測量的中子劑量差異。需要說明的是電離室無法測量中子能譜、平均能量和通量。

    Bonner球中子譜儀(Bonner sphere spectrometer)也稱為中子多球譜儀,由多個直徑不同的聚乙烯慢化球殼組成[28]。球殼中心放置體積較小的熱中子探測器,由3He正比管或6Li探測器,有時也用BF3管和對熱中子敏感的活化片或轉(zhuǎn)換體(如10B、6Li和235U)等徑跡探測器組成,其中3He對熱中子的敏感度比其他探測器高出一個量級以上。球殼的外徑通常為5~40 cm,取決于中子的最大能量。直徑越大,所能測量的中子能量越高,最內(nèi)層球殼的內(nèi)經(jīng)通常為2 cm左右,中子探測器必須足夠小能置入其中。作為測量中子通量的儀器由于涵蓋的能量范圍很廣(從熱中子到GeV的中子),各向同性響應(yīng),操作簡單,受到大多數(shù)實驗室的使用。根據(jù)中心探測器的原理,分為主動式球譜儀(中心探測器為3He正比計數(shù)器)和被動式球譜儀(中心探測器為金箔)。意大利的Frascati國家實驗室[28]在中子發(fā)生器產(chǎn)生的中子場環(huán)境下比較了主動式多球譜儀和被動式多球譜儀測量中子通量和能譜的差異,在允許的誤差范圍內(nèi),兩種探測器的結(jié)果基本是一致的。但是在醫(yī)用電子直線加速器產(chǎn)生的光中子環(huán)境中,主動式球譜儀往往不能得到滿意的測量結(jié)果,這是因為主動式多球譜儀中心的探測器有3He正比計數(shù)器構(gòu)成,高能射線的直接照射和室內(nèi)的光子散射形成脈沖式的光子中子混合場,與He元素通過光核反應(yīng)產(chǎn)生質(zhì)子,影響探測器的數(shù)據(jù)讀取工作。另外,過高的輻射頻率(一般在103~104MHz)會擾動相關(guān)電子器件的正常運作。因此在這種環(huán)境下一般采用被動式球譜儀來測量中子通量和能譜。Bedogni等[29]使用Bonner球譜儀測量動態(tài)的18 MV加速器動態(tài)MLC產(chǎn)生的中子劑量,以輸出處方為45 Gy的條件下,比較了IMRT技術(shù)與傳統(tǒng)放射治療產(chǎn)生的中子劑量差異。儀器放在SSD=80 cm的治療床平面刪,距離等中心處21 cm的位置處,結(jié)果表明IMRT技術(shù)在這一位置產(chǎn)生中子劑量為300 mSv,而傳統(tǒng)放射治療產(chǎn)生的貢獻(xiàn)為105mSv。這是因為IMRT技術(shù)要滿足相同的處方必須輸出更多的射線,隨著射野面積減小,產(chǎn)生的中子劑量是增加的;這是因為射野越小,機頭內(nèi)部互相碰撞的粒子越多,導(dǎo)致產(chǎn)生更多的中子。Amgarou等[30]使用Bonner球譜儀測量了治療室內(nèi)電中子(由電子碰撞直接產(chǎn)生的中子)及光中子的通量譜。Howell等[31]使用Bonner球譜儀測量了瓦里安21EX(15 MV、18 MV和20 MV),西門子ONCOR(15 MV、18MV)以及醫(yī)科達(dá)Presice(15MV、18 MV)所產(chǎn)生的光中子的平均能量、中子通量及機房內(nèi)周邊的劑量,以及中子源(機頭)的強度。數(shù)據(jù)表明對于不同的機器和能量檔,中子通量譜及能量譜的曲線形狀趨勢是一樣的。對于同一能量檔的不同機器,測得瓦里安加速器的數(shù)據(jù)最高,文獻(xiàn)認(rèn)為與機頭結(jié)構(gòu)如多葉光柵有關(guān)。

    Takam等[32]對患者和模體同時測量了光子劑量和中子劑量。患者身上的光子劑量用LiF熱釋光和AL2 O3探測器,中子劑量用CR-39徑跡探測器,對于模體所受的中子劑量是使用25.4 cm的Bonner球譜儀來測量,球譜儀中心的探測元件是熱釋光TLD-600和TLD-700。這兩種對于熱中子的敏感度不同,TLD-600與熱中子相互作用的散射截面遠(yuǎn)遠(yuǎn)大于TLD-700(1380/1),后者用來測量光子劑量。根據(jù)嵌套的球?qū)訑?shù)量不同以及中心探測器件不同,所測中子能量范圍也有所差別,使用的Bonner球譜儀探測的中子能量范圍在100 KeV和7 MeV。根據(jù)實際情況可以對多譜球?qū)舆M(jìn)行添加或減少以滿足所測中子能量需要。

    4 展望

    醫(yī)療技術(shù)在飛速發(fā)展的同時也要盡量避免其帶來的負(fù)影響,光中子作為醫(yī)用加速器產(chǎn)生的高能射線中不需要的粒子逐步得到重視和研究。但是與此相關(guān)的輻射防護(hù)措施并沒有得到及時地處理。蒙特卡羅模擬和測量中子劑量各有優(yōu)勢,只要建立了機頭構(gòu)造的模型,蒙特卡羅程序就能夠模擬所有粒子在不同位置的運動學(xué)參數(shù),但是對于加速器機頭多葉光柵運動過程中的粒子通量和能量變化目前無法計算,只能通過實驗儀器來監(jiān)測。儀器無法追蹤粒子的運動軌跡,只能測量固定位置的粒子物理參數(shù)。無論是用于基礎(chǔ)研究的高能物理(能量在GeV甚至TeV以上)還是醫(yī)用直線加速器產(chǎn)生的殺死癌細(xì)胞的X射線,蒙特卡羅模擬與實驗測量是技術(shù)上互補,數(shù)據(jù)結(jié)果上互相驗證的兩個重要工具。光中子作為X射線治療腫瘤過程中出現(xiàn)的“副產(chǎn)品”對于患者和工作人員既有劑量貢獻(xiàn)又有輻射危害,也是引起第二原發(fā)癌產(chǎn)生的重要物理因素。對于超過光核反應(yīng)閾值的光子能量所產(chǎn)生的中子能貢獻(xiàn)多少劑量,以及對人體產(chǎn)生多少生物損害,是當(dāng)前腫瘤放射治療領(lǐng)域不能忽視的一個研究內(nèi)容。

    [1]Zohar E.Particle physics:Quantum simulation of fundamental physics[J].Nature,2016,534(7608):480-481.

    [2]Grigorenko AN.Classical physics of particles and fields[J].Journal of Geometry and Physics,2014,86:94-100.

    [3]ICRU,Internal Commission on Radiation Units and Measurements.Quantities and units in radiation protection dosimetry.[A].ICRU Report 90,2016.

    [4]于一,趙燕,董放.醫(yī)用電子直線加速器的可疑不良事件研究[J].中國醫(yī)療設(shè)備,2014,29(6):61-63.

    [5]趙新春,周艷,張會敏,等.廣西醫(yī)用電子加速器驗收檢測泄漏輻射水平分析[J].中國醫(yī)學(xué)裝備,2014,11(2):9-12.

    [6]Chen W,Zheng R,Baade PD,et al.Cancer statistics in China,2015[J].CA Cancer J Clin,2016,66(2):115-132.

    [7]Syam Kumar SA,Sukumar P,Sriram P,et al. A patient-specific quality assurance study on absolute dose verification using ionizationchambers of different volumes in RapidArc treatments[J].Medical Dosimetry,2012,37(4):3595.

    [8]Cardenas CE,Nitsch PL,Kudchadker RJ,et al. Out-of-field doses and neutron dose equivalents for electron beams from modern Varian and Elekta linear accelerators[J].J Appl Clin Med Phys,2016,17(4):6216.

    [9]Berrington de Gonzalez A,Wong J,Kleinerman R,et al.Risk of second cancers according to radiation therapy technique and modality in prostate cancer survivors[J].Int J Radiat Oncol Biol Phys,2015,91(2):295-302.

    [10]Ghiasi H,Mesbahi A.Monte Carlo characterization of photoneutrons in the radiation therapy with high energy photons:a Comparison between simplified and full Monte Carlo models[J]. Iranian Journal of Radiation Research,2011,8(3):187-193.

    [11]Chibani O,Ma CM.Photonuclear dose calculations for high-energy photon beams from Siemens and Varian linacs[J].Med Phys,2003,30(8):1990-2000.

    [12]Zanini A,Durisi E,F(xiàn)asolo F,et al.Monte Carlo simulation of the photoneutron field in linac radiotherapy treatments with different collimation systems[J].Phys Med Biol,2004,49(4):571-582.

    [13]Pena J,F(xiàn)ranco L,Gómez F,et al.Monte Carlo study of Siemens PRIMUS photoneutron production[J].Phys Med Biol,2005,50(24):5921-5933.

    [14]Mao XS,Kase KR,Liu JC,et al.Neutron sources in the Varian Clinac 2100C/2300C medical accelerator calculated by the EGS4 code[J]. Health Physics,1997,72(4):524-529.

    [15]Becker J,Brunckhorst E,Schmidt R.Investigation of the neutron contamination in IMRT deliveries with a paired magnesium and boron coated magnesium ionization chamber system[J]. Radiother Oncol,2008,86(2):182-186.

    [16]Martínez-Ovalle SA,Barquero R,Gómezros JM,et al.Neutron dose equivalent and neutron spectra in tissue for clinical linacs operating at 15,18 and 20MV[J].Radiat Prot Dosimetry,2011,147(4):498-511.

    [17]Barquero R,Edwards T M,I?iguez M P,et al. Monte Carlo simulation estimates of neutron doses to critical organs of a patient undergoing 18MV x-ray LINAC-based radiotherapy[J]. Medical Physics,2005,32(12):3579-3588.

    [18]曹磊,張貴英,鄧君,等.15MeV醫(yī)用電子直線加速器治療室內(nèi)光核中子輻射場的研究[J].輻射防護(hù),2010(3):131-134.

    [19]范銳鋒,馬景芳,艾杰,等.銀活化探測系統(tǒng)測量X光機中(γ,n)中子注量實驗研究[J].核電子學(xué)與探測技術(shù),2009,29(6):1288-1291.

    [20]黃蔚玲,李泉鳳.中高能醫(yī)用電子直線加速器中光中子問題的研究[C]//全國醫(yī)用加速器技術(shù)交流大會,2003.

    [21]黃蔚玲,李泉鳳,林郁正,等.無損檢測用15MeV電子直線加速器主束內(nèi)光中子劑量的測量與MC模擬計算[J].高能物理與核物理,2006,30(S1):63-65.

    [22]Expósito MR,Romero-Hermida MI,Terrón JA,et al.Neutron contamination in medical linear accelerators operating at electron mode[M]//World Congress on Medical Physics and Biomedical Engineering,IFMBE Proceedings,39.2013:1225-1228.

    [23]Krmar M,Nikoli? D,Kuzmanovi? A,et al.The effect of a paraffin screen on the neutron dose at the maze door of a 15 MV linear accelerator[J].Medical Physics,2013,40(8):083902-083902.

    [24]Takam R,Bezak E,Liu G,et al.The use of enriched 6Li and 7Li Lif:Mg,Cu,Pglassrod thermoluminescent dosemeters for linear,accelerator out-of-field radiation dose measurements[J].Radiat Protect Dosimetry,2012,150(1):22-33.

    [25]Koivunoro H,Hyv?nen H,Uusisimola J,et al. Effect of the calibration in water and the build-up cap on the Mg(Ar) ionization chamber measurements[J].Appl Radiat Isot,2011,69(12):1901-1903.

    [26]Burmeister J,Kota C,Yudelev M,et al. Paired Mg and Mg(B)ionization chambers for the measurement of boron neutron capture dose in neutron beams[J].Medical Physics,1999,26(11):2482-2487.

    [27]Howell RM,F(xiàn)erenci MS,Hertel NE,et al. Measurements of secondary neutron dose from 15 MV and 18 MV IMRT[J].Radiat Prot Dosim etry,2005,115(1-4):508-512.

    [28]Aza E,Dinar N,Manessi GP,et al.A B onner Sphere Spectrometer for pulsed fields[J]. Radiation Prot Dosimetry,2016,168(2):149-153.

    [29]Bedogni R,Esposito A,Gentile A,et al. Comparing active and passive Bonner Sphere Spectrometers in the 2.5MeV quasi monoenergetic neutron field of the ENEA Frascati Neutron Generator (FNG)[J].Radiation Measure ments,2011,46(12):1757-1760.

    [30]Amgarou K,Lacoste V,Martin A,et al.Nerutron Spectrometry With a Passive Bonner Sphere System Around a Medical LINAC and Evaluation of the Associated Unfolding Uncertainties[J].IEEE Transactions on Nuclear Science,2009,56(5):2885-2895.

    [31]Howell RM,F(xiàn)erenci MS,Hertel N E,et al. Investigation of secondary neutron dose for 18MV dynamic MLC IMRT delivery[J].Medical Physics,2005,32(3):786-793.

    [32]Takam R,Bezak E,Liu G,et al.The use of enriched 6Li and 7Li Lif:Mg,Cu,P glassrod thermoluminescent dosemeters for linear,accelerator out-of-field radiation dose measurements[J].Radiat Protect Dosimetry,2012,150(1):22-33.

    The research progress of photo-neutron in the modern radiotherapy field/

    CHEN Ya-zheng,XU Ke, XU Jing-hui, et al// China Medical Equipment,2016,13(12):154-159.

    With the rapid development of new medical accelerators and the rapid increase in the number of patients with cancer, radiotherapy has become one of the main means to cure cancer. The unceasing progressive of medical accelerator technology always considers two aspects which includes the irradiation of the tumor, at the same time, and the maximum protection for normal organs, but these progressives never solve the radiation contamination problem from high-energy radiation when treat patient. Besides, the biological effects on human tissue of secondary particles produced by high energy particles, such as photons, photoneutrons, protons and so on, also shouldn’t be neglected. These problems has been paid more and more attention and widely research. Because the internal structure of the head of accelerator is Pb elemental and the mainly target area is W (Tungsten) element, it will exceed the photon reaction threshold to produce neutrons when the photon energy reaches more than 7 MeV. This reaction is called photonuclear reaction. Foreign radiotherapy institutes have carried out a wide range studies on the dose of photo-neutrons. Based on the characteristics of particles involved in complex interactions, Monte Carlo simulation and experimental measurements are performed to study the photon flux and energy spectrum of medical accelerators. There are different distribution regulation of flux of light neutrons, energy spectrum and dose, due to different design of medical accelerator ( mainly is the interior of machine head) lead to different properties. This paper based on documents review, and the status of photon neutrons generated by medical accelerators, the radiation protection during radiotherapy and the exploration of the physical factors of the second primary cancer are expounded from Monte Carlo simulation and experimental measurement in this paper.

    Photo-neutron; Photonuclear raction; Dose; Flux; Medical linear accelerator

    10.3969/J.ISSN.1672-8270.2016.12.045

    1672-8270(2016)12-0154-06

    R815.2

    A

    2016-10-10

    四川省青年科技基金(2015JQ0053)“現(xiàn)代放療技術(shù)中引發(fā)第二原發(fā)癌風(fēng)險的輻射劑量因素研究”;四川省衛(wèi)生和計劃生育委員會課題(16PJ517)“醫(yī)用加速器高檔位能量(10MV)產(chǎn)生的中子輻射測試研究”

    ①四川省腫瘤醫(yī)院放療中心 四川 成都 610041

    *通訊作者:dengwangpei@163.com

    陳亞正,男,(1983- ),博士,助理研究員。四川省腫瘤醫(yī)院放療中心,研究方向:醫(yī)學(xué)物理及輻射防護(hù)。

    猜你喜歡
    電離室中子加速器
    輪滑加速器
    化學(xué)工業(yè)的“加速器”
    全民小康路上的“加速器”
    少先隊活動(2021年6期)2021-07-22 08:44:24
    PTW 729電離室矩陣不同驗證方法用于宮頸癌術(shù)后調(diào)強放療計劃驗證結(jié)果分析
    3D打印抗中子輻照鋼研究取得新進(jìn)展
    等待“加速器”
    基于PLC控制的中子束窗更換維護(hù)系統(tǒng)開發(fā)與研究
    DORT 程序進(jìn)行RPV 中子注量率計算的可靠性驗證
    中子深度定量分析的相對分析法
    計算物理(2014年2期)2014-03-11 17:01:27
    Precise直線加速器電離室的維護(hù)和故障排查
    国产一区二区激情短视频| 色94色欧美一区二区| 狠狠婷婷综合久久久久久88av| 91av网站免费观看| 日本一区二区免费在线视频| 成人18禁高潮啪啪吃奶动态图| 国产成+人综合+亚洲专区| 激情视频va一区二区三区| 精品一品国产午夜福利视频| 亚洲av欧美aⅴ国产| 国产成人av教育| 久久人妻av系列| 美女国产高潮福利片在线看| 美女午夜性视频免费| 国产午夜精品久久久久久| 这个男人来自地球电影免费观看| 国产激情久久老熟女| 欧美日韩国产mv在线观看视频| 老熟妇仑乱视频hdxx| a级片在线免费高清观看视频| 精品福利永久在线观看| 多毛熟女@视频| 黑人猛操日本美女一级片| 亚洲精品国产色婷婷电影| 国产一区二区在线观看av| 午夜福利视频在线观看免费| 天堂俺去俺来也www色官网| 大型黄色视频在线免费观看| 男女午夜视频在线观看| 亚洲自偷自拍图片 自拍| 午夜日韩欧美国产| 午夜福利影视在线免费观看| 国产高清视频在线播放一区| 色在线成人网| 亚洲国产毛片av蜜桃av| 十八禁高潮呻吟视频| 午夜两性在线视频| 久久人人97超碰香蕉20202| 亚洲欧美色中文字幕在线| 人妻 亚洲 视频| 午夜福利一区二区在线看| 考比视频在线观看| 精品久久久久久久毛片微露脸| 狠狠狠狠99中文字幕| 亚洲av欧美aⅴ国产| 99riav亚洲国产免费| 日本撒尿小便嘘嘘汇集6| 大码成人一级视频| 91麻豆av在线| 亚洲免费av在线视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲成国产人片在线观看| 亚洲国产中文字幕在线视频| 精品国产乱子伦一区二区三区| 黑人操中国人逼视频| 欧美 亚洲 国产 日韩一| 日本av手机在线免费观看| 国产成+人综合+亚洲专区| 真人做人爱边吃奶动态| 国产精品电影一区二区三区 | 日韩大片免费观看网站| 中文字幕色久视频| 久久午夜综合久久蜜桃| 成人影院久久| 脱女人内裤的视频| 午夜福利影视在线免费观看| 亚洲五月色婷婷综合| 亚洲美女黄片视频| 亚洲,欧美精品.| 日韩视频一区二区在线观看| 精品免费久久久久久久清纯 | 国产av又大| 国产97色在线日韩免费| 多毛熟女@视频| 在线观看免费高清a一片| 50天的宝宝边吃奶边哭怎么回事| 在线 av 中文字幕| 美女午夜性视频免费| 18禁美女被吸乳视频| av片东京热男人的天堂| 在线观看舔阴道视频| 嫩草影视91久久| av有码第一页| 亚洲综合色网址| 色老头精品视频在线观看| 久久精品熟女亚洲av麻豆精品| 亚洲欧美一区二区三区黑人| svipshipincom国产片| 一边摸一边抽搐一进一出视频| 国产不卡一卡二| 久久99一区二区三区| 十分钟在线观看高清视频www| 高清欧美精品videossex| 少妇精品久久久久久久| 亚洲精华国产精华精| 在线播放国产精品三级| 动漫黄色视频在线观看| 人人妻人人澡人人爽人人夜夜| 女警被强在线播放| 精品国产国语对白av| 在线观看人妻少妇| 母亲3免费完整高清在线观看| 午夜成年电影在线免费观看| 黄色成人免费大全| 岛国毛片在线播放| 18禁裸乳无遮挡动漫免费视频| 黑人操中国人逼视频| 久热爱精品视频在线9| 一个人免费看片子| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩视频精品一区| 在线观看舔阴道视频| 国产xxxxx性猛交| 亚洲av欧美aⅴ国产| 91大片在线观看| 999久久久精品免费观看国产| 99国产精品免费福利视频| 99热国产这里只有精品6| 我的亚洲天堂| 日本五十路高清| 国产有黄有色有爽视频| 国产精品欧美亚洲77777| 精品少妇一区二区三区视频日本电影| 欧美日本中文国产一区发布| 狂野欧美激情性xxxx| 老司机亚洲免费影院| 国产高清激情床上av| 夜夜爽天天搞| 日日夜夜操网爽| 伊人久久大香线蕉亚洲五| 国产男女超爽视频在线观看| 首页视频小说图片口味搜索| 欧美成人午夜精品| 精品久久久久久久毛片微露脸| 久久久久久亚洲精品国产蜜桃av| 亚洲欧美激情在线| 国产男女超爽视频在线观看| 在线播放国产精品三级| 女人被躁到高潮嗷嗷叫费观| 精品国产一区二区三区久久久樱花| 在线av久久热| 欧美成人午夜精品| 国产在线免费精品| 菩萨蛮人人尽说江南好唐韦庄| av不卡在线播放| 国产日韩欧美在线精品| 十八禁人妻一区二区| 99精品在免费线老司机午夜| 精品一品国产午夜福利视频| 国产熟女午夜一区二区三区| 欧美精品啪啪一区二区三区| 久久久水蜜桃国产精品网| 国产一区二区在线观看av| 热99国产精品久久久久久7| 人妻 亚洲 视频| 欧美成狂野欧美在线观看| 精品国产乱子伦一区二区三区| 99久久人妻综合| 天天操日日干夜夜撸| 亚洲精品国产色婷婷电影| 久久婷婷成人综合色麻豆| 精品久久久久久久毛片微露脸| 亚洲avbb在线观看| 亚洲三区欧美一区| 90打野战视频偷拍视频| 黑人猛操日本美女一级片| 欧美 日韩 精品 国产| 亚洲九九香蕉| 青青草视频在线视频观看| 亚洲精品在线观看二区| 欧美日韩中文字幕国产精品一区二区三区 | 成人手机av| 国产成人精品久久二区二区免费| 国产精品国产av在线观看| 中文字幕人妻熟女乱码| 五月天丁香电影| 1024香蕉在线观看| aaaaa片日本免费| 老司机在亚洲福利影院| 成人免费观看视频高清| 国产精品 欧美亚洲| 欧美 日韩 精品 国产| 叶爱在线成人免费视频播放| 我的亚洲天堂| 国产主播在线观看一区二区| 国产精品秋霞免费鲁丝片| 亚洲va日本ⅴa欧美va伊人久久| 国产伦理片在线播放av一区| 桃花免费在线播放| 精品国产超薄肉色丝袜足j| 亚洲性夜色夜夜综合| 午夜激情久久久久久久| 亚洲全国av大片| 亚洲七黄色美女视频| 9色porny在线观看| 丁香六月天网| 波多野结衣一区麻豆| 不卡av一区二区三区| 亚洲,欧美精品.| 老司机午夜十八禁免费视频| 超碰97精品在线观看| 亚洲九九香蕉| 免费在线观看日本一区| 国产日韩欧美在线精品| 一区在线观看完整版| 一本大道久久a久久精品| 久久毛片免费看一区二区三区| 国产成人av教育| 黄色视频在线播放观看不卡| 欧美日本中文国产一区发布| av天堂久久9| 精品亚洲乱码少妇综合久久| 成人手机av| 超碰成人久久| 国产精品免费一区二区三区在线 | 欧美另类亚洲清纯唯美| 嫩草影视91久久| 久久精品91无色码中文字幕| 国产一区二区三区综合在线观看| 大码成人一级视频| 久久精品国产a三级三级三级| 免费看十八禁软件| 超碰成人久久| 成在线人永久免费视频| 精品一区二区三区av网在线观看 | 视频区图区小说| 黄色视频在线播放观看不卡| 女性生殖器流出的白浆| 欧美国产精品一级二级三级| 亚洲全国av大片| 亚洲黑人精品在线| 免费在线观看视频国产中文字幕亚洲| 女人久久www免费人成看片| av网站在线播放免费| 久久久久久人人人人人| 亚洲一区二区三区欧美精品| 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品美女久久av网站| 自拍欧美九色日韩亚洲蝌蚪91| 激情视频va一区二区三区| 国产亚洲精品久久久久5区| 我的亚洲天堂| 国产免费av片在线观看野外av| 亚洲国产av新网站| 纯流量卡能插随身wifi吗| 在线看a的网站| 99精品在免费线老司机午夜| av欧美777| 亚洲国产av影院在线观看| 水蜜桃什么品种好| 精品国产乱码久久久久久男人| 日本av手机在线免费观看| 国产xxxxx性猛交| 欧美日韩一级在线毛片| 日韩中文字幕欧美一区二区| 久久99一区二区三区| 中文字幕av电影在线播放| av欧美777| aaaaa片日本免费| 亚洲三区欧美一区| 久久久久久久久久久久大奶| 国产成+人综合+亚洲专区| 人妻一区二区av| 可以免费在线观看a视频的电影网站| 妹子高潮喷水视频| 亚洲国产欧美一区二区综合| 欧美 亚洲 国产 日韩一| 在线天堂中文资源库| 丰满人妻熟妇乱又伦精品不卡| 老熟女久久久| 别揉我奶头~嗯~啊~动态视频| 亚洲精品一卡2卡三卡4卡5卡| 国产成人精品在线电影| 国精品久久久久久国模美| 两个人免费观看高清视频| 亚洲成a人片在线一区二区| 国产成人欧美| 久久性视频一级片| 变态另类成人亚洲欧美熟女 | 亚洲专区国产一区二区| 18禁美女被吸乳视频| 久久 成人 亚洲| 国产精品免费一区二区三区在线 | 大型av网站在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 丝瓜视频免费看黄片| 亚洲av日韩在线播放| 亚洲欧洲精品一区二区精品久久久| 韩国精品一区二区三区| 母亲3免费完整高清在线观看| 在线播放国产精品三级| 国产欧美日韩一区二区三区在线| 叶爱在线成人免费视频播放| 中文字幕最新亚洲高清| 国产深夜福利视频在线观看| 久久久久国内视频| 亚洲av成人一区二区三| 国产精品香港三级国产av潘金莲| 久久性视频一级片| 老熟妇乱子伦视频在线观看| 日韩 欧美 亚洲 中文字幕| 亚洲精品自拍成人| 一级毛片精品| 欧美精品高潮呻吟av久久| 久久免费观看电影| 99国产精品99久久久久| 国产又色又爽无遮挡免费看| 在线观看免费午夜福利视频| 精品少妇内射三级| 80岁老熟妇乱子伦牲交| 亚洲熟女精品中文字幕| 精品久久久精品久久久| 午夜福利乱码中文字幕| 欧美日韩福利视频一区二区| 超碰成人久久| 精品少妇黑人巨大在线播放| 成人国语在线视频| 国产精品亚洲一级av第二区| 捣出白浆h1v1| 精品一区二区三区四区五区乱码| 亚洲欧美一区二区三区久久| 日本欧美视频一区| 久久天躁狠狠躁夜夜2o2o| 高潮久久久久久久久久久不卡| 99国产极品粉嫩在线观看| 欧美在线一区亚洲| 国产日韩一区二区三区精品不卡| 午夜视频精品福利| 狠狠婷婷综合久久久久久88av| 国产高清videossex| 久久久久久人人人人人| 久久中文字幕一级| 成人18禁高潮啪啪吃奶动态图| 欧美日韩av久久| 国产xxxxx性猛交| 超色免费av| 国产激情久久老熟女| 黑人巨大精品欧美一区二区mp4| 免费观看av网站的网址| 美国免费a级毛片| 亚洲精品美女久久av网站| 国产aⅴ精品一区二区三区波| 99re在线观看精品视频| 大码成人一级视频| 成年动漫av网址| 一本色道久久久久久精品综合| 亚洲专区中文字幕在线| 亚洲人成电影免费在线| 嫩草影视91久久| 欧美日韩一级在线毛片| 韩国精品一区二区三区| 精品一区二区三区四区五区乱码| 99久久人妻综合| 久久久精品国产亚洲av高清涩受| 午夜精品国产一区二区电影| 亚洲人成电影免费在线| 乱人伦中国视频| 最黄视频免费看| 另类精品久久| 国产熟女午夜一区二区三区| 不卡一级毛片| 69精品国产乱码久久久| 亚洲全国av大片| 亚洲人成电影免费在线| 大型黄色视频在线免费观看| 亚洲人成伊人成综合网2020| 乱人伦中国视频| 一区二区av电影网| 免费看a级黄色片| 午夜老司机福利片| 国产在线一区二区三区精| 人人妻人人澡人人看| 少妇粗大呻吟视频| 成人黄色视频免费在线看| 亚洲色图av天堂| 国产欧美亚洲国产| 国产一区二区在线观看av| 啪啪无遮挡十八禁网站| 久久青草综合色| 女人精品久久久久毛片| 一区二区三区激情视频| 成人亚洲精品一区在线观看| 国产精品久久电影中文字幕 | aaaaa片日本免费| 亚洲精品在线观看二区| 可以免费在线观看a视频的电影网站| 亚洲成人免费av在线播放| 国产成人免费观看mmmm| 亚洲欧美色中文字幕在线| 国产单亲对白刺激| 国产又爽黄色视频| 男女高潮啪啪啪动态图| a级片在线免费高清观看视频| 久久午夜综合久久蜜桃| 亚洲熟妇熟女久久| 中文字幕高清在线视频| 日韩欧美免费精品| 老司机亚洲免费影院| 亚洲国产精品一区二区三区在线| 亚洲色图av天堂| av超薄肉色丝袜交足视频| 国产精品98久久久久久宅男小说| 国产精品成人在线| 午夜视频精品福利| 亚洲人成伊人成综合网2020| 国产日韩欧美在线精品| 国产亚洲午夜精品一区二区久久| 亚洲第一青青草原| 午夜免费鲁丝| 亚洲熟妇熟女久久| 在线观看舔阴道视频| 日韩视频一区二区在线观看| 在线 av 中文字幕| 亚洲av国产av综合av卡| 12—13女人毛片做爰片一| 99精品久久久久人妻精品| 天天躁夜夜躁狠狠躁躁| 十分钟在线观看高清视频www| 一个人免费看片子| 欧美国产精品一级二级三级| 亚洲男人天堂网一区| 天堂8中文在线网| 亚洲久久久国产精品| 看免费av毛片| videosex国产| 岛国在线观看网站| 久久久久国内视频| 日韩免费av在线播放| 亚洲国产av影院在线观看| 丁香六月天网| 丰满人妻熟妇乱又伦精品不卡| 久久天堂一区二区三区四区| 啦啦啦在线免费观看视频4| 国产精品.久久久| 咕卡用的链子| 久久99一区二区三区| 日韩免费高清中文字幕av| 亚洲av第一区精品v没综合| 国产精品久久久久久精品电影小说| 中亚洲国语对白在线视频| 国产视频一区二区在线看| 中文欧美无线码| 国产欧美日韩一区二区三区在线| 99热网站在线观看| 亚洲美女黄片视频| 首页视频小说图片口味搜索| 天天躁狠狠躁夜夜躁狠狠躁| 日本撒尿小便嘘嘘汇集6| videos熟女内射| 亚洲国产av新网站| 国产深夜福利视频在线观看| 国产一区二区在线观看av| 国产在线一区二区三区精| av片东京热男人的天堂| 国产精品麻豆人妻色哟哟久久| 久久久久国产一级毛片高清牌| 两个人免费观看高清视频| av不卡在线播放| 国产伦人伦偷精品视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲五月婷婷丁香| 亚洲专区中文字幕在线| 婷婷丁香在线五月| 亚洲精品一二三| 国精品久久久久久国模美| 日韩视频一区二区在线观看| 激情在线观看视频在线高清 | 精品第一国产精品| 国产精品.久久久| 无限看片的www在线观看| 最近最新免费中文字幕在线| 亚洲精品中文字幕在线视频| 亚洲人成伊人成综合网2020| 男女无遮挡免费网站观看| 在线播放国产精品三级| 日韩熟女老妇一区二区性免费视频| av电影中文网址| 母亲3免费完整高清在线观看| 亚洲精品在线美女| 老司机午夜十八禁免费视频| 精品人妻在线不人妻| 色视频在线一区二区三区| 露出奶头的视频| 国产老妇伦熟女老妇高清| a在线观看视频网站| 好男人电影高清在线观看| 可以免费在线观看a视频的电影网站| 国产欧美日韩一区二区三| 精品国内亚洲2022精品成人 | 久久精品成人免费网站| 纯流量卡能插随身wifi吗| 18禁裸乳无遮挡动漫免费视频| 久热这里只有精品99| 黄色怎么调成土黄色| 99国产综合亚洲精品| 色老头精品视频在线观看| 在线观看免费日韩欧美大片| 午夜免费成人在线视频| 亚洲av片天天在线观看| 免费观看av网站的网址| 在线观看www视频免费| 两人在一起打扑克的视频| 国产又色又爽无遮挡免费看| 久久青草综合色| 亚洲九九香蕉| 日韩一区二区三区影片| 亚洲人成伊人成综合网2020| av视频免费观看在线观看| 人人妻人人添人人爽欧美一区卜| 国产精品av久久久久免费| 黄色成人免费大全| 国产欧美日韩综合在线一区二区| 日韩视频在线欧美| 午夜免费鲁丝| 亚洲七黄色美女视频| 成人特级黄色片久久久久久久 | 国产一区二区三区视频了| 日本撒尿小便嘘嘘汇集6| 国产主播在线观看一区二区| 真人做人爱边吃奶动态| 交换朋友夫妻互换小说| 18禁观看日本| 国产精品免费大片| 性高湖久久久久久久久免费观看| 热99国产精品久久久久久7| 在线观看免费日韩欧美大片| 丝袜喷水一区| 国产在线精品亚洲第一网站| 欧美激情高清一区二区三区| 国产精品亚洲一级av第二区| 国产成人免费观看mmmm| 18禁国产床啪视频网站| 午夜福利影视在线免费观看| 欧美日韩成人在线一区二区| 免费黄频网站在线观看国产| 最近最新中文字幕大全电影3 | 日韩人妻精品一区2区三区| 欧美日韩亚洲高清精品| 老司机午夜福利在线观看视频 | 国产高清videossex| 欧美日韩中文字幕国产精品一区二区三区 | 久久精品国产a三级三级三级| 午夜激情久久久久久久| 大香蕉久久网| 下体分泌物呈黄色| 无人区码免费观看不卡 | 中文字幕av电影在线播放| 最近最新中文字幕大全免费视频| 999精品在线视频| 人妻久久中文字幕网| 一本一本久久a久久精品综合妖精| 又黄又粗又硬又大视频| 国产色视频综合| 国产日韩欧美在线精品| xxxhd国产人妻xxx| 国产成人精品久久二区二区91| 飞空精品影院首页| 自拍欧美九色日韩亚洲蝌蚪91| 99在线人妻在线中文字幕 | av欧美777| 黄色毛片三级朝国网站| 亚洲精华国产精华精| 丝瓜视频免费看黄片| 久久午夜综合久久蜜桃| 亚洲九九香蕉| 999精品在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 满18在线观看网站| 女人被躁到高潮嗷嗷叫费观| 国产精品 欧美亚洲| 亚洲性夜色夜夜综合| 精品国产超薄肉色丝袜足j| 人人妻人人澡人人看| 久久久精品免费免费高清| 亚洲中文日韩欧美视频| 亚洲一区二区三区欧美精品| 精品亚洲成国产av| 久久国产精品人妻蜜桃| 男女之事视频高清在线观看| 老司机靠b影院| 大片电影免费在线观看免费| 色94色欧美一区二区| 女同久久另类99精品国产91| 国产成人欧美在线观看 | 国产精品偷伦视频观看了| 波多野结衣av一区二区av| 亚洲色图综合在线观看| 国产精品 国内视频| 国产亚洲精品第一综合不卡| 亚洲情色 制服丝袜| 精品免费久久久久久久清纯 | 女人久久www免费人成看片| 一级a爱视频在线免费观看| 色婷婷av一区二区三区视频| 精品少妇黑人巨大在线播放| 一边摸一边抽搐一进一出视频| 国产麻豆69| 91麻豆精品激情在线观看国产 | 99国产精品99久久久久| 亚洲av国产av综合av卡| 久久久国产欧美日韩av| 日日摸夜夜添夜夜添小说| a级毛片在线看网站| 免费看a级黄色片| 一级毛片女人18水好多| 欧美精品人与动牲交sv欧美| 人人妻,人人澡人人爽秒播| 两个人看的免费小视频| 国产成人精品无人区| 丝瓜视频免费看黄片| 国产成人av激情在线播放| 狠狠狠狠99中文字幕| 99久久人妻综合| 天堂8中文在线网| 色精品久久人妻99蜜桃|