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

    X射線熒光CT成像中熒光產(chǎn)額、退激時(shí)間、散射、偏振等關(guān)鍵物理問題計(jì)算與分析*

    2021-11-01 06:10:26張芝振李亮
    物理學(xué)報(bào) 2021年19期
    關(guān)鍵詞:康普頓散射截面偏振度

    張芝振 李亮?

    1) (清華大學(xué)工程物理系,北京 100084)

    2) (清華大學(xué),粒子技術(shù)與輻射成像教育部重點(diǎn)實(shí)驗(yàn)室,北京 100084)

    X 射線熒光CT(X-ray fluorescence computed tomography,XFCT)是一種使用X 射線熒光(X-ray fluorescence,XRF)實(shí)現(xiàn)功能性成像的新技術(shù),在生物醫(yī)學(xué)成像中表現(xiàn)出較大潛力.但是,X 射線穿過生物體的同時(shí)還會(huì)產(chǎn)生大量康普頓散射光子,對(duì)XRF 信號(hào)的采集形成很強(qiáng)的背景噪聲;因此,如何有效消除康普頓散射噪聲對(duì)于提高XFCT 成像質(zhì)量至關(guān)重要.本文研究總結(jié)了XFCT 成像過程中涉及的物理過程,包括:熒光的產(chǎn)額、退激發(fā)時(shí)間、熒光發(fā)射角分布、熒光偏振態(tài)、康普頓散射角分布與散射光偏振態(tài),并通過研究熒光與散射光物理性質(zhì)的差異尋找去除康普頓散射噪聲的方法.經(jīng)過物理過程推導(dǎo)和分析計(jì)算,發(fā)現(xiàn):1) 高原子序數(shù)元素的K 層熒光退激發(fā)時(shí)間極短,在現(xiàn)有探測(cè)器的時(shí)間分辨率條件下,無法分辨散射光與熒光;2) K 層發(fā)射熒光的角分布各向同性,康普頓散射角分布在與入射光偏振方向附近取得最小值,而且入射光線偏振度越高,散射光的微分截面越小,偏振光源將有利于減少康普頓散射噪聲;3) K 層熒光線偏振度為零,而康普頓散射光子在一些散射方向上具有一定線偏振度,因此偏振態(tài)的差異可能用于降低康普頓散射噪聲.

    1 引言

    X 射線熒光CT(X-ray fluorescence computed tomography,XFCT)是一種針對(duì)高原子序數(shù)元素的功能成像技術(shù).XFCT 通過探測(cè)X 射線與物體內(nèi)高原子序數(shù)元素發(fā)生光電效應(yīng)后產(chǎn)生的XRF 得到物質(zhì)內(nèi)該元素的濃度分布圖像,所以兼具透射型CT 結(jié)構(gòu)成像和發(fā)射型CT 功能成像的優(yōu)點(diǎn).Boisseau 和Grodzins[1]利用美國國家同步輻射光源首次實(shí)現(xiàn)XFCT 實(shí)驗(yàn).同步輻射光源具有準(zhǔn)單能和高線偏振度的特性,是XFCT 成像的理想光源,但是因?yàn)橥捷椛湔加脠龅睾统杀靖甙旱膯栴},所以基于同步輻射的XFCT 系統(tǒng)很難應(yīng)用于實(shí)際醫(yī)學(xué)應(yīng)用中.X 光機(jī)相比于同步輻射具有巨大的靈活性和經(jīng)濟(jì)性優(yōu)勢(shì),近年來有很多基于X 光機(jī)的XFCT 系統(tǒng)設(shè)計(jì)的研究[2?9].因?yàn)閄 光機(jī)光源采用電子軔致輻射的方式產(chǎn)生X 射線,所以其能譜是分布范圍較大的連續(xù)譜并且不具備線偏振特性.相比于同步輻射光源,X 光機(jī)光源會(huì)產(chǎn)生更大的散射背景,嚴(yán)重影響了XFCT 的成像質(zhì)量.

    去除散射背景的方法主要有兩種.第一種方法采用略高于和略低于熒光物質(zhì)的K-edge 的單能入射光分別采取兩組數(shù)據(jù),兩組數(shù)據(jù)之差可近似為熒光信號(hào)[10].第一種方法對(duì)光源單能性要求很高,一般采用同步輻射源.第二種方法利用散射光子能量連續(xù)分布的特點(diǎn),通過三次多項(xiàng)式擬合的方式估計(jì)散射,進(jìn)而在數(shù)據(jù)中減去散射的貢獻(xiàn)[2],該方法對(duì)探測(cè)器的能量分辨率要求高,一般使用單像素X 射線譜儀采集熒光信號(hào),導(dǎo)致掃描時(shí)間很長.

    逆康普頓散射源相比于同步輻射源不僅具有占地面積小和價(jià)格適中的優(yōu)勢(shì),而且具有同步輻射源準(zhǔn)單能和線偏振的特性,適合應(yīng)用于XFCT 成像系統(tǒng)設(shè)計(jì)中.逆康普頓散射源與同步輻射源降低XFCT 散射背景的原理都是應(yīng)用了線偏振X 光康普頓散射的角分布特性.Chi 等[11]蒙卡模擬研究了逆康普頓光源線偏振特性對(duì)XFCT 成像質(zhì)量的影響,發(fā)現(xiàn)逆康普頓光源能抑制XFCT 成像的康普頓散射背景約1.6 倍.

    此外,熒光與散射光的物理性質(zhì)的差別可能用于降低散射背景,比如是否有可能利用熒光與散射光飛行時(shí)間、角分布和偏振態(tài)的差別減少散射背景或?qū)ふ疑⑸浔尘白钚〉奶綔y(cè)方位.本文分析了XFCT中熒光與散射光相關(guān)物理過程的性質(zhì):熒光與散射光的角分布,熒光退激發(fā)時(shí)間與熒光產(chǎn)額,熒光與散射光的偏振態(tài).這些物理過程決定了XFCT 熒光成像的極限,通過這些物理過程性質(zhì)能夠指導(dǎo)XFCT 的系統(tǒng)設(shè)計(jì),減少XFCT 的康普頓散射背景和提高XRF 信號(hào)強(qiáng)度.

    2 XFCT 成像方法和物理原理

    在XFCT 所涉及的能量范圍內(nèi),X 射線與物質(zhì)相互作用方式有光電效應(yīng)、相干散射和康普頓散射.在X 射線能量為幾十keV 時(shí),相干散射截面與康普頓散射截面相當(dāng),所以在處理XFCT 散射問題時(shí)相干散射和康普頓散射一般都要考慮,其中相干散射主要沿著小角方向散射,對(duì)光子的傳播會(huì)有影響,但通過后面的討論可以看出,在大角度方向康普頓散射是主要的.

    XFCT 成像系統(tǒng)幾何設(shè)計(jì)有很多種方式[2?12],其中基于小孔成像、扇束X 射線源和XRF 探測(cè)器陣列的成像系統(tǒng)設(shè)計(jì)能大幅提高掃描速度[5],其XFCT 實(shí)驗(yàn)布置示意圖如圖1 所示.小孔準(zhǔn)直器作用是使XRF 探測(cè)器具有空間分辨能力.示意圖中被掃描物體是一個(gè)包含四個(gè)小圓柱的大圓柱模型,四個(gè)小圓柱里含有不同濃度的某種高原子序數(shù)元素的溶液.圖1 中描繪了某一條光路以及該路徑上可能產(chǎn)生的XRF 光子和散射光子,入射光與溶液中的重金屬元素原子發(fā)生光電效應(yīng),重金屬原子發(fā)射熒光(特征X 射線),熒光經(jīng)過小孔后,被XRF探測(cè)器陣列探測(cè),同時(shí)散射光子也會(huì)進(jìn)入XRF 探測(cè)器,形成散射背景.XRF 探測(cè)器可通過設(shè)置能窗閾值得到入射到每個(gè)探測(cè)器單元的相應(yīng)能窗的光子計(jì)數(shù).X 射線探測(cè)器陣列用于傳統(tǒng)的透射型CT成像,其得到的衰減系數(shù)分布圖可用于XFCT 成像的衰減校正.

    圖1 基于小孔成像方式的XFCT 實(shí)驗(yàn)設(shè)置Fig.1.Experimental setup based on pinhole imaging.

    為描述入射和散射X 射線的偏振態(tài),引入Stokes矢量[13]:

    其中I表示光強(qiáng)度;P1和P2表示在與傳播方向垂直的平面內(nèi)成 4 5?夾角的兩對(duì)正交方向下分別測(cè)量得到的平面偏振度;P3表示圓偏振度.光束的線偏振度和圓偏振度定義分別為[14]

    其中PL表示線偏振度;PC表示圓偏振度;Imax和Imin分別表示經(jīng)過理想線偏振片后光強(qiáng)的最大與最小值;IR和IL分別表示光束中右旋光和左旋光的強(qiáng)度.定義歸一化Stokes 矢量為

    Stokes 矢量與正交方向選擇有關(guān),選定ξ1?1 對(duì)應(yīng)的線偏振方向?yàn)閤軸,則ξ1+1 對(duì)應(yīng)的線偏振方向?yàn)閥軸,傳播方向?yàn)閦軸正向,并規(guī)定:迎著光子傳播方向(z軸正向)看時(shí),繞z軸順時(shí)針旋轉(zhuǎn)角度為正,繞z軸逆時(shí)針旋轉(zhuǎn)角度為負(fù),則在旋轉(zhuǎn)ψ角度后的新正交方向下的Stokes 矢量 [ 1,ζ]T與原正交方向下的Stokes 矢量 [ 1,ξ]T由變換矩陣M聯(lián)系,即 [ 1,ζ]TM[1,ξ]T.M的非零元為

    一束光經(jīng)過某個(gè)相互作用后,其Stokes 矢量會(huì)改變,通過一個(gè) 4×4 變換矩陣T聯(lián)系相互作用前后的Stokes 矢量:

    3 K 層熒光物理性質(zhì)

    3.1 熒光產(chǎn)額

    某個(gè)殼層或子殼層的熒光產(chǎn)額定義為該殼層產(chǎn)生一個(gè)空穴后,通過發(fā)射熒光退激發(fā)的概率.K 殼層產(chǎn)生空穴后,不同躍遷方式產(chǎn)生的熒光的能量和強(qiáng)度不同.一般地,殼層越高的電子通過退激發(fā)填充K 層空穴,產(chǎn)生的熒光的能量越高,某種熒光的強(qiáng)度正比于該熒光對(duì)應(yīng)躍遷的分支比.熒光產(chǎn)額的半經(jīng)驗(yàn)公式[15]:

    式中,ωK表示K 層熒光產(chǎn)額;Z表示原子序數(shù);Ci表示擬合參數(shù).(6)式中參數(shù)Ci可以由實(shí)驗(yàn)得到的不同元素?zé)晒猱a(chǎn)額數(shù)據(jù)擬合得到.Bambynek 等[15]和Hubbell 等[16]結(jié)合理論計(jì)算與實(shí)驗(yàn)給出擬合參數(shù)和元素?zé)晒猱a(chǎn)額數(shù)據(jù)表.

    K 層熒光產(chǎn)額ωK是原子序數(shù)Z的增函數(shù),對(duì)Z >60 的元素,ωK>90%,所以高原子序數(shù)原子的K 層空穴態(tài)通過發(fā)射熒光退激的概率最大.因此XFCT 應(yīng)該選擇高原子序數(shù)元素作為熒光元素,在高Z元素范圍內(nèi)熒光產(chǎn)額基本飽和,提高Z對(duì)增加熒光產(chǎn)額的收益很低,但提高Z可增大熒光能量,能量高的光子更容易穿透物體被探測(cè)器采集到,從而降低統(tǒng)計(jì)噪聲.

    K 殼層空穴態(tài)不同躍遷方式產(chǎn)生不同能量的熒光,熒光按末態(tài)空穴所在層分成兩組,即Kα代表來自L 層電子的躍遷產(chǎn)生的熒光,Kβ表示來自M 層及更高層電子躍遷產(chǎn)生的熒光,同一組內(nèi)的熒光能量相近.Ertu?ral 等[17]給出了原子序數(shù)16≤Z≤92 范圍內(nèi)59 種元素的Kβ/Kα熒光強(qiáng)度比測(cè)量值,Scofield[18]給出了10≤Z≤98 原子序數(shù)范圍內(nèi)Kβ/Kα熒光強(qiáng)度比理論計(jì)算值,從實(shí)驗(yàn)和理論計(jì)算結(jié)果可以看出,熒光強(qiáng)度比Kβ/Kα隨原子序數(shù)增大而增大,對(duì)于Z>60 的元素,Kβ/Kα熒光強(qiáng)度比在0.3 左右.通過熒光強(qiáng)度比和K 層熒光產(chǎn)額ωK可以得到K 層不同能量熒光的產(chǎn)額.

    入射X 射線也會(huì)與重金屬原子的L,M 和N等更高電子殼層的電子發(fā)生光電效應(yīng),且發(fā)射的熒光能量隨層數(shù)增高而減小,其中熒光能量最高的L 層熒光的能量一般為十幾keV,屬于低能區(qū)的光子,幾乎無法穿透模體,所以L,M 和N 等更高電子殼層的熒光在XFCT 熒光成像中可以忽略.

    3.2 熒光退激發(fā)時(shí)間

    如果熒光退激發(fā)半衰期很長,且長于探測(cè)器的時(shí)間分辨率,那么可以采用極短時(shí)間的X 光脈沖作為光源,在散射光子到達(dá)探測(cè)器后,延遲一段時(shí)間采集熒光光子,這樣就可以在時(shí)間上分開散射和熒光從而去除康普頓散射背景;如果探測(cè)器的時(shí)間分辨率大于熒光退激發(fā)半衰期,那么熒光和散射光子到達(dá)探測(cè)器的平均時(shí)間差小于探測(cè)器時(shí)間分辨率,從而無法分辨熒光和散射光子.為了證明這種方法是否可行,需要計(jì)算熒光退激發(fā)時(shí)間.原子發(fā)生光電效應(yīng)后處于激發(fā)態(tài),激發(fā)態(tài)是不穩(wěn)定的,假設(shè)t=0 時(shí)刻原子處于激發(fā)態(tài),則在以后的t時(shí)刻體系仍然停留在這個(gè)態(tài)的概率是

    式中,τ是激發(fā)態(tài)平均存在時(shí)間,亦稱激發(fā)態(tài)壽命,其與激發(fā)態(tài)能級(jí)寬度Γ關(guān)系為

    K層熒光分寬度ΓR與K層熒光產(chǎn)額的關(guān)系是[15]

    對(duì)于Z>40 的元素,K層空穴態(tài)能級(jí)寬度半經(jīng)驗(yàn)公式為[15]

    通過(8)式和(9)式得到K層熒光壽命τK為

    可見,K 層熒光壽命隨原子序數(shù)增加快速下降,對(duì)已知的XFCT 應(yīng)用的最低原子序數(shù)的元素釓(Z64),K 層熒光壽命為0.03 ps.當(dāng)前對(duì)伽馬光子探測(cè)器時(shí)間分辨率要求最高的PET 設(shè)備的符合時(shí)間分辨率在百皮秒量級(jí)[19],遠(yuǎn)大于中高原子序數(shù)元素的K 層熒光退激發(fā)時(shí)間,因此,無法通過時(shí)間差別來分辨熒光與散射光子進(jìn)而減少散射背景.

    3.3 熒光發(fā)射角分布與線偏振度

    原子內(nèi)殼層熒光發(fā)射角分布和偏振態(tài)在理論和實(shí)驗(yàn)上已經(jīng)有很多研究[20?26]:若空穴態(tài)的總角動(dòng)量量子數(shù)J1/2,則熒光發(fā)射角分布各向同性且線偏振度為零[20];若J>1/2,則熒光發(fā)射一般各向異性并且部分偏振[22].因?yàn)镵 層空穴態(tài)J1/2[22],所以K 層空穴態(tài)發(fā)射熒光的角分布各向同性且線偏振度為零.L,M 等不同子層具有不盡相同的總角動(dòng)量量子數(shù),而且某些子層的空穴態(tài)會(huì)通過Coster-Kronig 躍遷將空穴轉(zhuǎn)移到其他子層,因此L,M 等殼層的空穴態(tài)發(fā)射熒光更加復(fù)雜多樣,剔除Coster-Kronig 效應(yīng)的影響后,不同子層空穴態(tài)發(fā)射熒光角分布與線偏振度也滿足上述角動(dòng)量規(guī)則[24].

    4 XFCT 中的散射

    4.1 散射光角分布

    實(shí)驗(yàn)室參考系下單個(gè)自由電子與線偏振光的康普頓散射截面Klein-Nishina 公式為[25]

    式中,r0是經(jīng)典電子半徑;k和k0分別是出射和入射光子能量,以電子靜質(zhì)量能為單位;e'和e分別是出射光偏振方向和入射光偏振方向; d?是立體角元.假設(shè)在出射方向n放置偏振器,該偏振器透光軸沿e'方向,一個(gè)偏振態(tài)為e的入射光子經(jīng)過康普頓散射到出射方向n,則該光子通過偏振器的概率正比于n方向上的 dσKN/d?[27].散射方向和偏振方向如圖2 所示,入射光沿Z軸正向入射,用矢量k0(|k0|k0)表示;散射光沿方向出射,用矢量k(|k|k)表示.e是描述入射光偏振方向的單位矢量,θ∠SOZ是散射角,φ∠AOX是散射的方位角.y和y'表示垂直于散射平面SOX的線偏振方向,x和x'表示平行于平面SOX的線偏振方向,位于平面x′Sy′內(nèi)的e'表示散射光的偏振方向,e'與y'夾角為β,顯然,x',y'和e'垂直kf.從圖2描述的幾何關(guān)系得到 (e·e′)2cos2β(1?sin2θcos2φ).

    圖2 偏振康普頓散射示意圖Fig.2.Schematic diagram of polarized Compton scattering.

    對(duì)出射光相互垂直的兩個(gè)偏振方向x'和y'的分截面求和得到入射光為線偏振光且不測(cè)量散射光偏振的實(shí)驗(yàn)中的截面[25]為

    一般地,如果入射光是線偏振度為PL(0

    式中,pp 表示部分偏振(partially polarized),相應(yīng)的圖2 中X軸方向取為入射光透過理想偏振片后的透射強(qiáng)度最大方向.注意到公式中k/k01/[1+k0(1?cosθ)] 與 方位角φ無關(guān),所以公式對(duì)φ積分后得到截面沿散射角θ的分布與線偏振度PL無關(guān),進(jìn)而總截面與PL無關(guān).

    下面討論(14)式的最小值,做變換xk/k0,則

    式中,x∈[1/(1+2k0),1].記(15)式右邊方括號(hào)內(nèi)的項(xiàng)為函數(shù)g(x,φ),則

    ?g/?x在x1/(1+2k0) 處的值為

    令f(k0,a)0,得到k0關(guān)于a的唯一解,記為k0K(a).易知,當(dāng)f(k0,a)≤0 時(shí),方程h(x,a)0 在定義域內(nèi)存在唯一解,記為xX(a).討論函數(shù)(15)的最小值得到如下結(jié)論:若k0

    由于K(1+PL)≥K(1)≈1.74,熒光CT 所用X 射線能量一般小于0.2 個(gè)電子靜能(100 keV),遠(yuǎn)小于K(1),所以最終得到截面最小值點(diǎn)對(duì)應(yīng)的θmin:

    不同入射光能量下的截面最小點(diǎn)θmin隨偏振度P的變化如圖3 所示,可見熒光CT 所用的X光能量范圍內(nèi),微分截面最小值對(duì)應(yīng)散射角在90°附近.偏振度越高且入射能量越低,則θmin越接近90°.由于部分偏振光可以統(tǒng)計(jì)描述為完全偏振光(PL1 )和一個(gè)完全非偏振光(PL0)的疊加,所以最小微分截面方向接近入射光完全偏振成分的電矢量方向.

    圖3 不同入射能量時(shí),θmin 隨線偏 振度變化Fig.3.θmin varying with polarization for different incident energy.

    圖4 描述了不同入射能量下,微分截面最小值隨偏振度的變化.可見在熒光CT 成像的能量范圍內(nèi),隨入射光偏振度的增大,最小微分截面值減小.這為熒光CT 去散射提供了思路,即入射光偏振度越高,在θθmin,φ0 or π 處的康普頓散射本底越低.

    圖4 不同入射能量時(shí)最小微分截面隨偏振度的變化Fig.4.Minimum differential cross section varying with polarization for different incident energy.

    前面討論的情形是單個(gè)自由電子與入射光的康普頓散射截面,這種情形很理想,方便分析.實(shí)際上物質(zhì)中的電子都不是自由的,此時(shí),一個(gè)原子的總散射截面與其各個(gè)電子在自由狀態(tài)下散射截面的總和一般不相等,引入非相干散射函數(shù)S(q,Z),則一個(gè)原子的非相干散射微分截面為[29]

    式中,Z是原子序數(shù);xsin(θ/2)/λ0,

    λ0是入射光子波長(單位為?).非相干散射函數(shù)的理論計(jì)算方法和數(shù)值表可在相關(guān)文獻(xiàn)中查尋[30?33].非相干散射函數(shù)滿足 0 ≤S(x,Z)/Z≤1,對(duì)于低Z原子的散射,只有在x較小(軟X 射線或小角度散射,x<1 )時(shí),S(x,Z)/Z才顯著小于1.例如100 keV 入射光子,波長約為0.124 ?,入射到氧原子(Z8 ),散射角滿足θ<4?才使S(x,Z)/Z <0.5,實(shí)際實(shí)驗(yàn)中,入射光束總有發(fā)散角和寬度,在小散射角處測(cè)量熒光信號(hào)存在很強(qiáng)的入射光背景,所以實(shí)驗(yàn)設(shè)置應(yīng)該在較大角度處測(cè)量熒光信號(hào),因此在這個(gè)測(cè)量角度范圍內(nèi)單個(gè)自由電子微分截面最小值的分析依然是有效的,即考慮非相干散射函數(shù)后,入射光線偏振優(yōu)勢(shì)方向是依然是近似的散射強(qiáng)度最小點(diǎn).

    單個(gè)自由電子的Klein-Nishina 散射微分截面公式在kk0→0 時(shí)得到單個(gè)自由電子的Thomson散射截面公式:

    若入射光線偏振度為PL,則相應(yīng)的散射截面公式為

    由(22)式易知,最小微分截面方向?yàn)棣圈?2 且φ0,π,所以最小微分截面方向?yàn)槿肷涔馔耆癯煞值碾娛噶糠较?從經(jīng)典電動(dòng)力學(xué)角度出發(fā)也可以得到相同的散射截面公式,即電子在入射光電場驅(qū)動(dòng)下振動(dòng)輻射電磁波.

    一個(gè)原子的相干散射截面為[25]

    其中F(x,Z) 是原子形狀因子,可通過文獻(xiàn)[30?32]查詢?cè)有螤钜蜃雍屠碚撚?jì)算方法.原子形狀因子隨x的變化趨勢(shì)與非相干散射函數(shù)相反,x越大F(x,Z) 越接近零,相反地,x趨于零時(shí)F(x,Z) 趨于Z.

    根據(jù)(20)式和(23)式,一個(gè)原子的康普頓散射與相干散射微分截面的總和為

    (20)式、(23)式和(24)式忽略光的偏振對(duì)原子形狀因子和散射因子的影響[26],所以根據(jù)(14)式和(22)式的最小值的討論,一個(gè)原子的總微分截面的最小值點(diǎn)一定滿足φ0,π.

    一個(gè)分子的散射截面角分布為[34]

    其中m(x) 表示分子干涉函數(shù);s(x) 和f(x) 表示采用獨(dú)立原子模型(independent atomic modelling,IAM)計(jì)算得到的非相干散射函數(shù)和分子形狀因子[33],即

    其中Fi(x,Zi) 表示分子中第i種元素原子的形狀因子;Si(x,Zi) 表示第i種元素的非相干散射函數(shù);ni表示分子中第i種元素的原子個(gè)數(shù);m(x) 表示分子干涉函數(shù)[34].

    計(jì)算得到水分子在方位角φ0,π 處的微分散射截面隨散射角θ的分布如圖5 所示.從圖5(a)可以看出,相干散射微分截面隨散射角增大快速減小,在大角度處,微分散射截面由主要由非相干散射貢獻(xiàn).由圖5(b)中總微分截面角分布隨線偏振度的變化可知,線偏振度越高,垂直方向散射的微分截面越小.PL0.5 時(shí)不同入射光能量下的總微分截面如圖5(c)所示,可見在XFCT 的光源能量范圍內(nèi),不同能量光子的水分子散射截面最小值均在入射光線偏振優(yōu)勢(shì)方向附近,且該方向上微分截面相近.

    圖5 φ=0,π 處微分截面隨散射角變化 (a) PL=1 時(shí)相干、非相干和總微分截面;(b) 不同線偏振度時(shí)的總微分截面;(c) PL=0.5 時(shí)不同入射能量下的總微分截面Fig.5.Differential cross section varying with θ at φ=0,π :(a) Incoherent,coherent and total differential cross section at PL=1 ;(b) total differential cross section varying with θ for different PL ;(c) total differential cross section varying with θ for different incident energy at PL=0.5.

    4.2 散射光的偏振

    聯(lián)系康普頓散射前后光束Stokes 矢量的變換矩陣T為[25]

    其中k0k0n0,kkn,S是電子初始自旋方向,n0是入射光動(dòng)量方向,n是散射光動(dòng)量方向.在各向同性材料中,電子初始自旋方向的平均值為0,所以矩陣T中含S的矩陣元為零.在變換矩陣T成立的坐標(biāo)系中要求垂直散射平面方向完全線偏振光的歸一化Stokes 矢量為即y和y'軸垂直于散射平面,由圖2 易知,坐標(biāo)系xOy可 由XOY繞OZ軸旋轉(zhuǎn)ψ?φ得 到,則 由(4a)式、(4b)式和(4c)式可知,在XOY坐標(biāo)系內(nèi)Stokes 矢量為的入射光轉(zhuǎn)換到坐標(biāo)系xOy后歸一化Stokes 矢量變?yōu)?/p>

    在各向同性材料某點(diǎn)發(fā)生康普頓散射后,散射光Stokes 矢量變?yōu)?[If,ζ′]T,即:

    由(29)式亦可得到散射光偏振的主方向與散射平面法向夾角β,符號(hào)規(guī)定為:迎著光子動(dòng)量方向,在表示散射光偏振的局部坐標(biāo)系x′Sy′內(nèi),y'軸以最小角度旋轉(zhuǎn)到e',若其旋轉(zhuǎn)方向?yàn)轫槙r(shí)針,則β為正,否則為負(fù).規(guī)定β的符號(hào)后,在局部坐標(biāo)x′Sy′旋轉(zhuǎn)β角得到的新坐標(biāo)系中,散射光Stokes矢量的第二個(gè)分量應(yīng)達(dá)到正的最大值,據(jù)此條件可求得

    其中函數(shù) a tan2(y,x) 為雙參數(shù)反正切函數(shù).

    在φ0 or π 散射平面上散射光的歸一化Stokes矢量、線偏振度和偏振主方向角為β分別為

    特別地,在θπ/2,φ0,π 散射方向上散射光的歸一化Stokes 矢量和線偏振度為

    圖6 展示了在不同入射光能量下φ0,π 散射平面內(nèi)隨散射角θ的變化.圖7 展示了不同能量下θπ/2,φ0,π 散射方向上散射光線偏振度隨入射光線偏振度的變化.若入射光為完全線偏振光(1 ),則θπ/2,φ0,π 散射方向上散射光的線偏振度為零,這與(12)式是相符的.從圖6 可以看出,在XFCT 所涉及的能量范圍內(nèi),相同散射角處不同能量下散射光線偏振度差距較小,而在θπ/2 附近,部分偏振的入射光,其散射光也具有較高偏振度.從圖7 可知,只有當(dāng)入射光線偏振度接近完全線偏振時(shí),θπ/2 的散射光的偏振度才顯著降低.

    圖6 φ=0 or π 處 隨散 射角變化Fig.6. varying with θ at φ=0, π.

    圖7 不同能量下 θ =π/2,φ=0, π 處散射光線偏振度隨入射光線偏振度變化Fig.7.The linear polarization of scattering photons varying with at θ =π/2,φ=0, π for different incident energy.

    從(30)式、(31)式和(32a)式可知:完全非偏振入射光的散射光的偏振主方向始終垂直于散射平面,部分偏振入射光的散射光偏振;部分偏振入射光在φ0,π 平面內(nèi)的散射光偏振度是散射角的函數(shù),且存在兩個(gè)零點(diǎn),左零點(diǎn)左側(cè)和右零點(diǎn)的右側(cè)偏振優(yōu)勢(shì)方向平行于散射平面,在兩個(gè)零點(diǎn)之間,偏振優(yōu)勢(shì)方向垂直于散射平面,在兩個(gè)零點(diǎn)處散射光完全非偏振;在φ0,π 散射平面內(nèi),完全線偏振入射光的散射光在θπ/2 時(shí)的偏振優(yōu)勢(shì)方向均平行于散射平面,在θπ/2 時(shí)散射光完全非偏振.

    根據(jù)以上關(guān)于散射光的線偏振度及熒光線偏振度的分析,并假設(shè)存在理想的X 光偏振片,因?yàn)镵 層熒光的線偏振度為零,K 層熒光經(jīng)過理想偏振片后減少50%,而散射光具有一定線偏振度,比如當(dāng)入射光完全非偏振時(shí),在單個(gè)自由電子這個(gè)理想的情況下,垂直于入射光方向的散射光的線偏振度可達(dá)到90%,經(jīng)過理想偏振片后散射光可減少95%,因此,在經(jīng)過偏振片后,雖然熒光強(qiáng)度減半,但是熒光強(qiáng)度與散射強(qiáng)度比值即信噪比在理想情況下可提高十倍,所以在原理上,散射光的線偏振度和K 層熒光線偏振度的差異可用于提高熒光CT 成像的信噪比.

    5 結(jié)論

    本文分析了XFCT 的物理過程及其性質(zhì).計(jì)算熒光退激發(fā)時(shí)間和熒光產(chǎn)額發(fā)現(xiàn),中高Z元素?zé)晒猱a(chǎn)額基本飽和,而熒光退激發(fā)時(shí)間遠(yuǎn)遠(yuǎn)小于現(xiàn)有探測(cè)器時(shí)間分辨率,因此,熒光退激發(fā)時(shí)間無法用于分辨散射光與熒光,提高原子序數(shù)可以增加熒光的能量,熒光能量越高,其穿透性越強(qiáng),有助于提高探測(cè)到的熒光光子數(shù),減小統(tǒng)計(jì)噪聲,提高信噪比.偏振光源有助于減少康普頓散射背景,入射光偏振度不影響康普頓散射總截面,但對(duì)截面的角分布有影響,部分偏振入射光在偏振方向附近具有最小的散射截面,在此方向上探測(cè)熒光光子有助于降低康普頓散射背景,并且對(duì)于完全非偏振或部分偏振的入射光,在大部分方向上的散射光仍然具有一定的線偏振度,而K 層熒光線偏振度為零,因此熒光和散射光線偏振度差異可以用于進(jìn)一步降低熒光探測(cè)中的康普頓散射背景,由于實(shí)際成像過程的復(fù)雜性,通過偏振性質(zhì)來減少散射背景需要進(jìn)一步模擬和實(shí)驗(yàn)驗(yàn)證.

    猜你喜歡
    康普頓散射截面偏振度
    基于Geant4模擬的康普頓散射研究
    LHCb =8 TeV的Drell-Yan-Z→e+e-數(shù)據(jù)對(duì)部分子分布函數(shù)的影響
    基于偏振度的目標(biāo)微觀表面偏振特性分析
    基于微波倍頻源太赫茲頻段雷達(dá)散射截面測(cè)量
    基于可調(diào)偏振度源驗(yàn)證偏振光譜強(qiáng)度調(diào)制系統(tǒng)
    115In中子非彈性散射截面的實(shí)驗(yàn)測(cè)量及蒙特卡羅修正
    核技術(shù)(2016年4期)2016-08-22 09:05:22
    史上著名的物理實(shí)驗(yàn)——康普頓效應(yīng)的發(fā)現(xiàn)
    物理之友(2015年9期)2015-02-25 03:11:07
    平板形目標(biāo)的量子雷達(dá)散射截面計(jì)算
    基于康普頓散射本底扣除的航空γ能譜測(cè)量譜線比大氣氡校正方法
    自然光在兩種介質(zhì)分界面上反射和折射時(shí)的偏振特性討論
    成年版毛片免费区| 欧美国产精品va在线观看不卡| 人人妻,人人澡人人爽秒播| 久久午夜亚洲精品久久| 每晚都被弄得嗷嗷叫到高潮| 老熟妇乱子伦视频在线观看| 国产单亲对白刺激| 搡老妇女老女人老熟妇| 中文字幕人成人乱码亚洲影| 黑人巨大精品欧美一区二区mp4| 正在播放国产对白刺激| 搡老妇女老女人老熟妇| 好男人电影高清在线观看| 正在播放国产对白刺激| av视频在线观看入口| 少妇被粗大的猛进出69影院| 成人18禁高潮啪啪吃奶动态图| 亚洲性夜色夜夜综合| 亚洲精品国产区一区二| 99久久精品国产亚洲精品| 国产一卡二卡三卡精品| 手机成人av网站| 18禁裸乳无遮挡免费网站照片 | 日韩精品青青久久久久久| 亚洲第一电影网av| 一级作爱视频免费观看| 乱人伦中国视频| 国产精品野战在线观看| 两人在一起打扑克的视频| 亚洲人成77777在线视频| 夜夜看夜夜爽夜夜摸| 成年女人毛片免费观看观看9| 久久久久久久久久久久大奶| 亚洲久久久国产精品| av在线播放免费不卡| 久久欧美精品欧美久久欧美| 999久久久国产精品视频| 精品福利观看| 精品国产乱子伦一区二区三区| 午夜精品久久久久久毛片777| 国产免费av片在线观看野外av| 动漫黄色视频在线观看| 熟女少妇亚洲综合色aaa.| 高清毛片免费观看视频网站| 亚洲七黄色美女视频| 亚洲人成电影观看| 狂野欧美激情性xxxx| 国产精品一区二区在线不卡| 亚洲国产欧美网| 日韩精品中文字幕看吧| 12—13女人毛片做爰片一| 久久国产精品男人的天堂亚洲| 天堂√8在线中文| a级毛片在线看网站| 免费在线观看黄色视频的| 婷婷丁香在线五月| 97超级碰碰碰精品色视频在线观看| 叶爱在线成人免费视频播放| 日韩三级视频一区二区三区| 久久久水蜜桃国产精品网| 亚洲男人天堂网一区| 18禁国产床啪视频网站| 91成人精品电影| 国语自产精品视频在线第100页| 黑人操中国人逼视频| 伦理电影免费视频| 免费不卡黄色视频| 亚洲人成77777在线视频| 一区二区三区激情视频| 一区二区三区国产精品乱码| 正在播放国产对白刺激| 午夜免费激情av| 国产精品二区激情视频| 亚洲视频免费观看视频| 一级毛片高清免费大全| 欧美久久黑人一区二区| 精品国产一区二区久久| 日韩大码丰满熟妇| 两个人视频免费观看高清| 国产一区二区激情短视频| 亚洲,欧美精品.| 日本五十路高清| 欧美日本视频| 国产高清videossex| 别揉我奶头~嗯~啊~动态视频| 国产欧美日韩精品亚洲av| 精品乱码久久久久久99久播| 亚洲片人在线观看| 身体一侧抽搐| 一本综合久久免费| 人人澡人人妻人| 欧美丝袜亚洲另类 | 精品国产国语对白av| 夜夜夜夜夜久久久久| 99国产综合亚洲精品| 国产欧美日韩一区二区三| 国产蜜桃级精品一区二区三区| 欧美一区二区精品小视频在线| 精品国产乱子伦一区二区三区| 午夜福利成人在线免费观看| 好男人在线观看高清免费视频 | 久久久精品欧美日韩精品| 男女午夜视频在线观看| 大型黄色视频在线免费观看| 日韩中文字幕欧美一区二区| 在线观看舔阴道视频| 精品久久久久久成人av| 精品第一国产精品| 亚洲精品中文字幕一二三四区| 嫩草影院精品99| 亚洲黑人精品在线| 99国产综合亚洲精品| 亚洲av电影在线进入| 久久香蕉国产精品| 国产精品二区激情视频| 久久狼人影院| 亚洲国产看品久久| 自线自在国产av| 亚洲欧美激情在线| 日本黄色视频三级网站网址| 精品不卡国产一区二区三区| 97人妻天天添夜夜摸| 在线国产一区二区在线| 精品国产乱子伦一区二区三区| 久久人人爽av亚洲精品天堂| 免费观看精品视频网站| 国内久久婷婷六月综合欲色啪| 国产成人啪精品午夜网站| 色综合欧美亚洲国产小说| 亚洲人成伊人成综合网2020| 桃色一区二区三区在线观看| www国产在线视频色| 自拍欧美九色日韩亚洲蝌蚪91| av片东京热男人的天堂| 久久久久久久久久久久大奶| 国产精品美女特级片免费视频播放器 | 久久久久久久精品吃奶| 51午夜福利影视在线观看| 午夜福利免费观看在线| 免费人成视频x8x8入口观看| 久久人妻福利社区极品人妻图片| 国产三级黄色录像| АⅤ资源中文在线天堂| 88av欧美| 国产主播在线观看一区二区| av中文乱码字幕在线| avwww免费| 国产精品影院久久| 日本欧美视频一区| av欧美777| 久久人人97超碰香蕉20202| 国产精品一区二区精品视频观看| 搞女人的毛片| 搡老熟女国产l中国老女人| 亚洲 欧美一区二区三区| 最新美女视频免费是黄的| 久久国产精品男人的天堂亚洲| ponron亚洲| 黄网站色视频无遮挡免费观看| 色婷婷久久久亚洲欧美| 国产精品一区二区在线不卡| 精品国产乱子伦一区二区三区| 黄色丝袜av网址大全| 50天的宝宝边吃奶边哭怎么回事| 最新在线观看一区二区三区| 大码成人一级视频| 久久草成人影院| 啦啦啦观看免费观看视频高清 | 久久精品aⅴ一区二区三区四区| 日韩精品免费视频一区二区三区| 国产三级在线视频| 欧美午夜高清在线| 久久中文看片网| 长腿黑丝高跟| 国产精品久久久久久精品电影 | 久热爱精品视频在线9| 亚洲一区二区三区不卡视频| 一级a爱视频在线免费观看| 欧美久久黑人一区二区| 我的亚洲天堂| 99久久国产精品久久久| 69精品国产乱码久久久| 久久天躁狠狠躁夜夜2o2o| 国产一区二区激情短视频| 美女大奶头视频| 免费在线观看影片大全网站| 久久中文看片网| 精品久久久久久久人妻蜜臀av | svipshipincom国产片| 伦理电影免费视频| 久久久久久久久免费视频了| 精品久久久久久久久久免费视频| 一级毛片精品| 国内久久婷婷六月综合欲色啪| 中文字幕高清在线视频| 香蕉久久夜色| av电影中文网址| 老熟妇乱子伦视频在线观看| 女人爽到高潮嗷嗷叫在线视频| 成人精品一区二区免费| 国产精品亚洲美女久久久| 欧美日韩瑟瑟在线播放| 国产精品98久久久久久宅男小说| 老熟妇乱子伦视频在线观看| 亚洲色图综合在线观看| 可以在线观看毛片的网站| 国产成人免费无遮挡视频| 一级毛片女人18水好多| 免费在线观看视频国产中文字幕亚洲| 老司机深夜福利视频在线观看| 成人精品一区二区免费| 黑人欧美特级aaaaaa片| 嫁个100分男人电影在线观看| 午夜亚洲福利在线播放| 国产蜜桃级精品一区二区三区| 久久婷婷成人综合色麻豆| 亚洲精品在线美女| 韩国av一区二区三区四区| 亚洲色图av天堂| 欧美精品亚洲一区二区| 成人av一区二区三区在线看| 日韩中文字幕欧美一区二区| 免费在线观看影片大全网站| 亚洲专区中文字幕在线| 国产精品久久久久久亚洲av鲁大| 国产精品久久久av美女十八| 搡老岳熟女国产| 久久婷婷人人爽人人干人人爱 | 国产精品 欧美亚洲| 久久久久久久久免费视频了| 精品一区二区三区av网在线观看| 成年女人毛片免费观看观看9| 最新美女视频免费是黄的| 国产激情欧美一区二区| 国产精品综合久久久久久久免费 | 色综合亚洲欧美另类图片| 一个人免费在线观看的高清视频| 国产精品 国内视频| 男人操女人黄网站| 欧美黄色淫秽网站| √禁漫天堂资源中文www| 人妻丰满熟妇av一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 女同久久另类99精品国产91| 啦啦啦韩国在线观看视频| 国产高清videossex| 欧美一级a爱片免费观看看 | 男女床上黄色一级片免费看| 久久中文看片网| 嫩草影院精品99| 欧美成人性av电影在线观看| 波多野结衣一区麻豆| 黄片大片在线免费观看| 亚洲精华国产精华精| 日韩高清综合在线| 黄片大片在线免费观看| 99久久国产精品久久久| svipshipincom国产片| 两个人看的免费小视频| 精品国产乱子伦一区二区三区| 国产精品98久久久久久宅男小说| 天天躁夜夜躁狠狠躁躁| 亚洲av美国av| 久久精品国产亚洲av高清一级| 一区在线观看完整版| 桃色一区二区三区在线观看| 国产成人av激情在线播放| 美国免费a级毛片| 国产精品乱码一区二三区的特点 | 制服丝袜大香蕉在线| 久久久久久久精品吃奶| 欧美日韩一级在线毛片| 麻豆一二三区av精品| 一级黄色大片毛片| 午夜免费成人在线视频| 黄色视频不卡| 夜夜爽天天搞| 亚洲国产中文字幕在线视频| 国产1区2区3区精品| 一边摸一边抽搐一进一出视频| 亚洲av成人不卡在线观看播放网| 国产不卡一卡二| 18禁美女被吸乳视频| 久99久视频精品免费| 国产1区2区3区精品| 亚洲欧美日韩另类电影网站| 久久久久久免费高清国产稀缺| 美女免费视频网站| 亚洲精品粉嫩美女一区| 天天添夜夜摸| a级毛片在线看网站| 精品一区二区三区四区五区乱码| 国产精品一区二区免费欧美| 精品国内亚洲2022精品成人| 欧美绝顶高潮抽搐喷水| 亚洲精品国产区一区二| 桃色一区二区三区在线观看| 大香蕉久久成人网| 国内精品久久久久久久电影| 午夜福利在线观看吧| 视频在线观看一区二区三区| 免费在线观看黄色视频的| 亚洲avbb在线观看| 欧美精品啪啪一区二区三区| 国产午夜福利久久久久久| 搡老岳熟女国产| 成人免费观看视频高清| 国产熟女午夜一区二区三区| 多毛熟女@视频| 99精品在免费线老司机午夜| 精品久久久久久久人妻蜜臀av | 不卡一级毛片| 99国产精品免费福利视频| 中文字幕另类日韩欧美亚洲嫩草| 可以在线观看的亚洲视频| 亚洲成国产人片在线观看| 91老司机精品| 女性被躁到高潮视频| ponron亚洲| 精品人妻在线不人妻| 99精品欧美一区二区三区四区| 日本免费a在线| e午夜精品久久久久久久| 亚洲情色 制服丝袜| 制服人妻中文乱码| 好看av亚洲va欧美ⅴa在| 欧美乱色亚洲激情| 久久香蕉精品热| 一级a爱片免费观看的视频| 久久久久久国产a免费观看| 久久性视频一级片| 亚洲avbb在线观看| 一级毛片高清免费大全| 亚洲一区中文字幕在线| 亚洲激情在线av| 嫩草影院精品99| 久久这里只有精品19| 91大片在线观看| 韩国av一区二区三区四区| 999久久久精品免费观看国产| 国产真人三级小视频在线观看| 91av网站免费观看| 在线观看日韩欧美| 夜夜躁狠狠躁天天躁| 欧美+亚洲+日韩+国产| 亚洲av日韩精品久久久久久密| 18禁裸乳无遮挡免费网站照片 | 大陆偷拍与自拍| 人人妻人人澡人人看| 免费在线观看日本一区| 亚洲五月色婷婷综合| 91麻豆精品激情在线观看国产| 久久久久国内视频| svipshipincom国产片| 成人国产综合亚洲| 少妇被粗大的猛进出69影院| 欧美在线一区亚洲| 狂野欧美激情性xxxx| 精品福利观看| 一进一出好大好爽视频| 亚洲av第一区精品v没综合| www.熟女人妻精品国产| 亚洲人成网站在线播放欧美日韩| av在线播放免费不卡| 久久久久久久久久久久大奶| 国产日韩一区二区三区精品不卡| 久久久久久免费高清国产稀缺| 99热只有精品国产| 欧美黄色淫秽网站| 天天一区二区日本电影三级 | 别揉我奶头~嗯~啊~动态视频| 亚洲欧美日韩高清在线视频| 一区二区三区激情视频| 一个人免费在线观看的高清视频| 国产97色在线日韩免费| 成人精品一区二区免费| 亚洲精品一区av在线观看| 黄色视频,在线免费观看| 欧美丝袜亚洲另类 | 免费观看人在逋| av片东京热男人的天堂| 又黄又爽又免费观看的视频| 19禁男女啪啪无遮挡网站| 成年女人毛片免费观看观看9| 中文字幕精品免费在线观看视频| 婷婷精品国产亚洲av在线| 成人免费观看视频高清| e午夜精品久久久久久久| 麻豆一二三区av精品| 国产精品影院久久| av在线天堂中文字幕| 欧美另类亚洲清纯唯美| 黑人操中国人逼视频| 纯流量卡能插随身wifi吗| 久久久久久大精品| 亚洲欧美激情综合另类| 无限看片的www在线观看| 国产精品,欧美在线| 如日韩欧美国产精品一区二区三区| 欧美久久黑人一区二区| 亚洲电影在线观看av| or卡值多少钱| 亚洲精品国产一区二区精华液| 黄网站色视频无遮挡免费观看| 一二三四社区在线视频社区8| 亚洲人成电影免费在线| 亚洲最大成人中文| 免费av毛片视频| 一区二区三区精品91| 国产成+人综合+亚洲专区| 欧美激情 高清一区二区三区| 制服诱惑二区| 日日摸夜夜添夜夜添小说| 日韩免费av在线播放| 国产精品二区激情视频| 后天国语完整版免费观看| 欧美最黄视频在线播放免费| 身体一侧抽搐| 天堂√8在线中文| 一卡2卡三卡四卡精品乱码亚洲| av中文乱码字幕在线| 99香蕉大伊视频| 日韩中文字幕欧美一区二区| 欧美在线一区亚洲| 亚洲av成人一区二区三| 在线免费观看的www视频| 午夜免费激情av| 久久久久久久久免费视频了| 亚洲精品中文字幕在线视频| 亚洲无线在线观看| 91大片在线观看| 日日摸夜夜添夜夜添小说| 午夜福利18| 国产午夜精品久久久久久| 亚洲欧美激情综合另类| 一级a爱视频在线免费观看| 亚洲精品久久成人aⅴ小说| 日本 欧美在线| 久久青草综合色| 亚洲第一av免费看| 如日韩欧美国产精品一区二区三区| 脱女人内裤的视频| 国产一区在线观看成人免费| 国产区一区二久久| 日本黄色视频三级网站网址| 天堂√8在线中文| 日本撒尿小便嘘嘘汇集6| 国产精品综合久久久久久久免费 | 女人被躁到高潮嗷嗷叫费观| 日本 欧美在线| 久久精品国产综合久久久| 丁香欧美五月| 久久九九热精品免费| 免费无遮挡裸体视频| 国产三级黄色录像| 少妇的丰满在线观看| 国语自产精品视频在线第100页| 在线免费观看的www视频| 国产欧美日韩综合在线一区二区| 精品少妇一区二区三区视频日本电影| 一二三四社区在线视频社区8| 操美女的视频在线观看| 国产精品1区2区在线观看.| 成人免费观看视频高清| 国产高清videossex| 国产三级在线视频| 18禁裸乳无遮挡免费网站照片 | 精品欧美国产一区二区三| 国产伦一二天堂av在线观看| 亚洲,欧美精品.| 久久久精品国产亚洲av高清涩受| 日韩国内少妇激情av| 亚洲精品久久成人aⅴ小说| 99在线视频只有这里精品首页| 国产成+人综合+亚洲专区| 国产免费av片在线观看野外av| 一区在线观看完整版| 无限看片的www在线观看| 一区二区日韩欧美中文字幕| 给我免费播放毛片高清在线观看| 亚洲国产欧美网| 国内精品久久久久久久电影| xxx96com| 中文字幕人成人乱码亚洲影| 好男人在线观看高清免费视频 | 日韩大尺度精品在线看网址 | 国产精品美女特级片免费视频播放器 | 欧美激情极品国产一区二区三区| 亚洲精品国产一区二区精华液| 亚洲专区中文字幕在线| www.www免费av| 亚洲va日本ⅴa欧美va伊人久久| 深夜精品福利| 亚洲第一电影网av| 亚洲视频免费观看视频| 精品日产1卡2卡| 悠悠久久av| 一区在线观看完整版| 国语自产精品视频在线第100页| 露出奶头的视频| 色综合婷婷激情| 午夜福利,免费看| 熟女少妇亚洲综合色aaa.| 美女免费视频网站| 在线观看日韩欧美| 免费在线观看视频国产中文字幕亚洲| 美女高潮喷水抽搐中文字幕| 久久精品国产99精品国产亚洲性色 | 色婷婷久久久亚洲欧美| 19禁男女啪啪无遮挡网站| 午夜福利成人在线免费观看| 纯流量卡能插随身wifi吗| av超薄肉色丝袜交足视频| 亚洲va日本ⅴa欧美va伊人久久| 黄色视频,在线免费观看| 黑人欧美特级aaaaaa片| 日韩 欧美 亚洲 中文字幕| 可以免费在线观看a视频的电影网站| 色哟哟哟哟哟哟| 黄片大片在线免费观看| 在线观看66精品国产| 午夜亚洲福利在线播放| 色播在线永久视频| 国产av精品麻豆| 国产在线精品亚洲第一网站| 国产精品,欧美在线| √禁漫天堂资源中文www| 搡老岳熟女国产| 亚洲人成伊人成综合网2020| 免费在线观看影片大全网站| 国产av又大| 国产精品一区二区精品视频观看| or卡值多少钱| 夜夜躁狠狠躁天天躁| 天堂动漫精品| 亚洲国产中文字幕在线视频| 国产精品久久久av美女十八| 国产麻豆69| 欧美日韩黄片免| 两性夫妻黄色片| 午夜激情av网站| 91av网站免费观看| 一级毛片高清免费大全| 久久中文看片网| 法律面前人人平等表现在哪些方面| 此物有八面人人有两片| a在线观看视频网站| 一级黄色大片毛片| 91av网站免费观看| 丝袜美足系列| 88av欧美| 亚洲欧洲精品一区二区精品久久久| 亚洲美女黄片视频| 欧美日韩福利视频一区二区| 黑人巨大精品欧美一区二区mp4| 亚洲性夜色夜夜综合| 黄色a级毛片大全视频| 亚洲成人国产一区在线观看| 成熟少妇高潮喷水视频| 亚洲欧美日韩无卡精品| 精品一区二区三区四区五区乱码| 18禁观看日本| 好男人在线观看高清免费视频 | 精品国产国语对白av| 男女下面进入的视频免费午夜 | 亚洲片人在线观看| 两人在一起打扑克的视频| 久热爱精品视频在线9| 他把我摸到了高潮在线观看| 成熟少妇高潮喷水视频| ponron亚洲| 啦啦啦 在线观看视频| 黑人欧美特级aaaaaa片| 亚洲久久久国产精品| 精品一区二区三区av网在线观看| 亚洲第一青青草原| 校园春色视频在线观看| 国产又爽黄色视频| 国产亚洲欧美98| 午夜影院日韩av| 老司机午夜十八禁免费视频| 男人舔女人的私密视频| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看免费视频日本深夜| 啪啪无遮挡十八禁网站| 亚洲国产精品合色在线| 国内毛片毛片毛片毛片毛片| 欧美日本视频| 少妇裸体淫交视频免费看高清 | 亚洲精品中文字幕在线视频| 国产精华一区二区三区| 精品不卡国产一区二区三区| 在线观看午夜福利视频| 淫妇啪啪啪对白视频| 国产免费av片在线观看野外av| 亚洲第一电影网av| 18禁国产床啪视频网站| 99国产精品99久久久久| 国产一区二区三区在线臀色熟女| 咕卡用的链子| 国内毛片毛片毛片毛片毛片| 香蕉久久夜色| 免费看美女性在线毛片视频| 亚洲一区二区三区不卡视频| 亚洲av第一区精品v没综合| 亚洲三区欧美一区| 成年人黄色毛片网站| 丝袜在线中文字幕| 久久久国产成人免费| 免费在线观看影片大全网站| 日韩一卡2卡3卡4卡2021年| 欧美在线黄色| 丰满人妻熟妇乱又伦精品不卡| 91麻豆精品激情在线观看国产| 丰满人妻熟妇乱又伦精品不卡|