亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        外差激光干涉儀周期非線性誤差形成機(jī)理與補(bǔ)償方法

        2020-06-11 06:43:56胡鵬程付海金楊宏興楊睿韜譚久彬
        關(guān)鍵詞:外差干涉儀光束

        王 越,胡鵬程,付海金,楊宏興,楊睿韜,譚久彬

        (1.哈爾濱工業(yè)大學(xué) 超精密光電儀器工程研究所,哈爾濱 150080;2.精密儀器技術(shù)及智能化工業(yè)和信息化部重點(diǎn)實(shí)驗(yàn)室(哈爾濱工業(yè)大學(xué)),哈爾濱 150080)

        外差激光干涉測(cè)量?jī)x具有測(cè)量精度高、測(cè)量速度快、非接觸、可直接溯源至長(zhǎng)度基準(zhǔn)等優(yōu)點(diǎn),是高精度位移測(cè)量的代表性精密儀器[1-4].隨著超精密加工制造、精密計(jì)量、引力波探測(cè)等科學(xué)研究與工業(yè)技術(shù)領(lǐng)域的發(fā)展,對(duì)位移測(cè)量精度的需求逐漸由納米量級(jí)過渡為亞納米甚至皮米量級(jí).在以光刻機(jī)為代表的超精密加工制造領(lǐng)域,為實(shí)現(xiàn)集成電路特征線寬達(dá)到5 nm甚至1 nm,外差激光干涉儀對(duì)硅片臺(tái)和掩膜臺(tái)位置的測(cè)量精度需要達(dá)到0.1 nm甚至數(shù)十皮米量級(jí)[5-6].在質(zhì)量量子化基準(zhǔn)裝置中,外差激光干涉儀負(fù)責(zé)對(duì)能量/功率線圈的位移的精密測(cè)量,若要實(shí)現(xiàn)10-8量級(jí)以上的溯源精度,外差激光干涉儀的測(cè)量精度需要突破0.1 nm[7-11].在引力波探測(cè)裝置LISA中,為了能夠在數(shù)千至百萬千米的等效臂長(zhǎng)下實(shí)現(xiàn)不同頻段引力波的探測(cè),外差激光干涉儀需要精確識(shí)別出皮米量級(jí)的位移[12-15].

        面臨著上述亞納米/皮米量級(jí)超精密位移測(cè)量精度需求的挑戰(zhàn),周期非線性誤差的問題逐漸成為外差激光干涉儀測(cè)量精度進(jìn)一步提升的障礙.周期非線性誤差是外差激光干涉儀的一種原理性誤差,其隨著位移的變化表現(xiàn)為周期性的變化規(guī)律,當(dāng)目標(biāo)每運(yùn)動(dòng)半個(gè)波長(zhǎng)時(shí),該誤差至少出現(xiàn)一次數(shù)納米至數(shù)十納米的峰值.外差激光干涉測(cè)量技術(shù)誕生后不久,Quenelle就發(fā)現(xiàn)了周期非線性誤差的存在[16].隨后,Sutton通過實(shí)驗(yàn)進(jìn)一步驗(yàn)證了非線性誤差隨測(cè)量位移周期性變化的規(guī)律[17].此后,國(guó)內(nèi)外研究人員對(duì)周期非線性誤差的理論模型和抑制方法進(jìn)行了廣泛的理論與實(shí)驗(yàn)研究[18-24].針對(duì)外差激光干涉儀中的周期非線性誤差問題,本文重點(diǎn)分析外差激光干涉儀中的兩類周期非線性誤差的形成機(jī)理,對(duì)具有代表性的周期非線性誤差優(yōu)化補(bǔ)償方法進(jìn)行了研究.本文旨在為抑制外差激光干涉測(cè)量技術(shù)中非線性誤差,進(jìn)一步提高外差激光干涉儀測(cè)量精度提供有益的參考.

        1 外差激光干涉儀中的周期非線性誤差形成機(jī)理

        在外差激光干涉儀中,根據(jù)誤差來源不同,周期非線性誤差可以分為兩類:第1類周期非線性誤差來源于共光路的測(cè)量光、參考光無法完全分離引起的“雙頻激光交叉混疊”,第2類周期非線性誤差來源于主要測(cè)量光束引起的部分虛反射光束的“多階多普勒頻移虛反射光束自混疊”.下面詳細(xì)分析兩類周期非線性誤差的形成機(jī)理.

        1.1 外差激光干涉儀中的第1類周期非線性誤差

        在理想外差激光干涉儀中,偏振分光鏡PBS將共軸傳輸?shù)膬蓚€(gè)頻率光f1、f2完全分離.但是實(shí)際中,如圖1所示,由于雙頻激光器存在橢偏化、安裝誤差、偏振分光鏡不理想等原因[25-30],導(dǎo)致在測(cè)量光路或者參考光路中存在著另一頻率的泄漏光,這種現(xiàn)象就是雙頻激光交叉混疊.

        圖1 外差激光干涉儀中雙頻激光交叉混疊示意圖

        由于雙頻激光交叉混疊,實(shí)際外差激光干涉儀中的測(cè)量與參考光束表示為

        (1)

        式中:A、B分別為主要測(cè)量和參考光束的幅值,α、β分別為參考和測(cè)量光路中泄漏光的幅值;Γ1=α/A,Γ2=β/B稱為泄漏系數(shù);Δφ=2πfd為測(cè)量位相,fd為目標(biāo)運(yùn)動(dòng)所引入的多普勒頻移;ω1=2πf1,ω2=2πf2.如圖2(a)所示,含有泄漏光的測(cè)量光束與參考光束干涉后形成測(cè)量信號(hào)為

        Im∝ABcos(Δωt+Δφ)+(Aβ+Bα)cos(Δωt)+αβcos(Δωt-Δφ).

        (2)

        式中:Δω=ω1-ω2;第1項(xiàng)為主要測(cè)量信號(hào)(main measurement signal, MMS),后兩項(xiàng)分別為0階和-1階寄生干涉信號(hào)(parasite interference signal,PIS),對(duì)其進(jìn)行仿真可以得到測(cè)量信號(hào)頻譜如圖2(b)所示,其中設(shè)置雙頻激光頻差Δf=5 MHz,多普勒頻移fd=1 MHz,泄漏系數(shù)Γ1=Γ2=0.1.可見受到雙頻激光交叉混疊影響的信號(hào)中,除了主要測(cè)量信號(hào)外,還包含有一階和二階寄生干涉信號(hào).

        圖2 雙頻激光交叉混疊誤差形成機(jī)理和測(cè)量信號(hào)頻譜分布

        對(duì)式(2)進(jìn)行變形,可得

        Im∝ABcos(Δωt+Δφ+Δφnonlin).

        (3)

        式中Δφnonlin1為非線性誤差項(xiàng),可以表示為

        (4)

        對(duì)于一般的外差激光干涉儀,Γ1?1,Γ2?1,則式(4)可近似為

        (5)

        根據(jù)外差激光干涉儀測(cè)量位移與位相的關(guān)系,可計(jì)算得外差激光干涉儀第1類周期非線性誤差為

        (6)

        由式(6)可見,第1類周期非線性誤差中包含有兩項(xiàng),其中第1項(xiàng)與第2項(xiàng)分別為一階和二階周期非線性誤差.設(shè)定主要參考和測(cè)量光束幅值A(chǔ)=B=1,對(duì)式(6)進(jìn)行仿真如圖3所示.圖3(a)為不同泄漏系數(shù)下的非線性誤差曲線,從圖中可以看出,隨著泄漏系數(shù)的增大,非線性誤差也隨之增大.圖3(b)為泄漏系數(shù)為0.1的情況下,一階和二階非線性誤差曲線,從圖中可以看出,一階非線性誤差要遠(yuǎn)大于二階非線性誤差,這也說明了非線性誤差大小主要取決于一階非線性誤差的大小,同時(shí)一階和二階寄生干涉信號(hào)引入的非線性誤差都具有明顯的周期性,且周期分別為2π和π.

        圖3 雙頻激光交叉混疊導(dǎo)致的非線性誤差曲線

        1.2 外差激光干涉儀中的第2類周期非線性誤差

        伴隨著位移測(cè)量精度要求逐漸過渡到亞納米至皮米量級(jí),哈爾濱工業(yè)大學(xué)胡鵬程等[23-24]發(fā)現(xiàn)外差激光干涉儀中存在著一種來源與表現(xiàn)形式都不同于第1類周期非線性誤差的新型周期非線性誤差——第2類周期非線性誤差.研究表明第2類周期非線性誤差來源于外差激光干涉儀中測(cè)量光束在光學(xué)界面處所產(chǎn)生的部分虛反射光引入的多階多普勒頻移虛反射光束自混疊.外差激光干涉儀的干涉光路由多個(gè)光學(xué)元件構(gòu)成,存在著大量的光學(xué)界面.理論上主要測(cè)量光經(jīng)過的每個(gè)光學(xué)界面上都會(huì)產(chǎn)生虛反射光束.圖4為多階多普勒頻移(Doppler frequency shift,DFS)光束形成的示意圖.圖4(a)為主要測(cè)量光束的傳播路徑,光束在目標(biāo)反射鏡上反射1次,因此具有一階多普勒頻移;圖4(b)中,主要測(cè)量光束穿過PBS,在1/4波片(QWP)外表面發(fā)生反射,產(chǎn)生零階多普勒頻移虛反射光束;圖4(c)、4(d)分別表示具有二階和三階多普勒頻移的虛反射光束產(chǎn)生過程.

        圖4 多階多普勒頻移形成過程示意圖[23]

        在外差激光干涉儀中,主要測(cè)量光束表示為

        Em=Acos(ω1t+Δφ).

        (7)

        根據(jù)圖4分析,包含多階多普勒頻移分量的虛反射光束可表示為

        Egh=γ0cos(ω1t)+γ1cos(ω1t+Δφ)+

        γ2cos(ω1t+2Δφ)+γkcos(ω1t+kΔφ).

        (8)

        式中γk=AR|k-1|/2為各階虛反射光束幅值,R為虛反射率.

        參考鏡固定不動(dòng),在參考光路中沒有多階多普勒頻移產(chǎn)生,因此,參考光束可表示為

        Er=Bcos(ω2t).

        (9)

        如圖5(a)所示,虛反射光束、測(cè)量光束與參考光束發(fā)生干涉后,得到測(cè)量信號(hào)為

        Im∝(Em+Er+Egh)·(Em+Er+Egh)*.

        (10)

        式(10)經(jīng)過化簡(jiǎn),得到

        Im∝Γ0cos(Δωt)+Γ1cos(Δωt+Δφ)+

        Γ2cos(Δωt+2Δφ)+…+Γkcos(Δωt+kΔφ).

        (11)

        式中:Γ0-Γk為各干涉分量的幅值,Γk=Bγk=ABR|k-1|/2;Γ1cos(Δωt+Δφ)為主要測(cè)量信號(hào);Γkcos(Δωt+kΔφ)為k階寄生干涉信號(hào).

        設(shè)置雙頻激光頻差Δf=5 MHz,多普勒頻移fd=1 MHz,虛反射率R=1%,對(duì)式(11)進(jìn)行仿真可以得到測(cè)量信號(hào)頻譜如圖5(b)所示.相比于雙頻激光交叉混疊模型,在多階多普勒頻移虛反射光束自混疊模型中,沒有-1階寄生干涉信號(hào),并且信號(hào)中含有大量高階寄生干涉信號(hào).

        對(duì)各階寄生干涉信號(hào)的影響進(jìn)行獨(dú)立分析,k階寄生干涉信號(hào)引入的非線性誤差可表示為

        (12)

        由于Γk?Γ1,對(duì)式(12)進(jìn)行一階泰勒展開后,可得由k階寄生干涉信號(hào)引入的周期非線性誤差為

        (13)

        根據(jù)外差激光干涉儀測(cè)量位移與位相的關(guān)系,可計(jì)算得外差激光干涉儀第2類周期非線性誤差為

        (14)

        由式(14)可以看出,各階非線性誤差隨虛反射率的增大而增大,且階數(shù)越大,誤差越小.對(duì)式(14)進(jìn)行仿真得到不同階次非線性誤差曲線如圖6所示,可見多階多普勒頻移虛反射引入的非線性誤差隨著測(cè)量位相呈現(xiàn)周期性變化,且多普勒頻移的階次越高,變化周期越??;非線性誤差的大小隨著多普勒頻移階次的增大而減小,隨著虛反射率的增大而增大.

        圖5 多階多普勒頻移虛反射光束自混疊模型非線性誤差形成機(jī)理與測(cè)量信號(hào)頻譜分布

        綜上分析,第1類與第2類周期非線性誤差形成機(jī)理和表現(xiàn)形式均有所不同.第1類周期非線性誤差來源于雙頻激光交叉混疊,非線性誤差中僅包含有一階和二階非線性誤差,誤差大小主要與外差激光干涉儀中的光學(xué)泄漏系數(shù)相關(guān).第2類周期非線性誤差來源于多階多普勒頻移虛反射光束自混疊,誤差中包含有高階非線性誤差,誤差大小主要與外差激光干涉儀中光學(xué)界面的虛反射率相關(guān).一般來說,在外差激光干涉儀中,第1類周期非線性誤差大于第2類周期非線性誤差,并且二者同時(shí)存在.

        2 周期非線性誤差補(bǔ)償方法

        非線性誤差補(bǔ)償方法是抑制周期非線性誤差影響,提高外差激光干涉儀測(cè)量精度的重要措施.根據(jù)第1節(jié)的分析,第1類和第2類周期非線性誤差的來源與形成原理互不相同,因此針對(duì)兩者的非線性誤差補(bǔ)償措施也有所不同,下面分別介紹.

        2.1 第1類周期非線性誤差補(bǔ)償方法

        針對(duì)第1類周期非線性誤差的補(bǔ)償,研究人員已經(jīng)進(jìn)行了廣泛而深入的研究,并提出了多種補(bǔ)償方法.這些非線性誤差補(bǔ)償方法大體上可以分為3種:第1種補(bǔ)償方法主要通過調(diào)整光路來減小非線性誤差[31-35],德國(guó)聯(lián)邦物理技術(shù)研究院提出了基于正交測(cè)量的外差激光干涉儀周期非線性誤差補(bǔ)償方法[33],如圖7(a)所示,該方法利用PBS將測(cè)量信號(hào)分為兩個(gè)偏振方向正交的信號(hào)同時(shí)進(jìn)行探測(cè),得到兩個(gè)測(cè)量結(jié)果并相加,測(cè)量的周期非線性誤差項(xiàng)由于符號(hào)相反相消而減小,從而達(dá)到補(bǔ)償周期非線性誤差的目的,這種方法能將外差激光干涉儀周期非線性誤差補(bǔ)償至幾納米.第2種補(bǔ)償方法一般通過信號(hào)處理算法的改進(jìn)來對(duì)包含有非線性誤差的信號(hào)進(jìn)行補(bǔ)償[36-40].韓國(guó)標(biāo)準(zhǔn)與科學(xué)研究院提出了基于橢圓擬合的外差激光干涉儀周期非線性誤差補(bǔ)償方法[36],如圖7(b)所示,該方法通過橢圓擬合法提取橢圓參數(shù)對(duì)非線性誤差進(jìn)行補(bǔ)償從而去除非線性誤差,這種方法實(shí)際中可以將非線性誤差補(bǔ)償至亞納米量級(jí).可見現(xiàn)有的周期非線性誤差補(bǔ)償方法普遍比較復(fù)雜,而且難以將周期非線性誤差抑制到0.1 nm以下.第3種非線性誤差補(bǔ)償方法為空間分離式外差激光干涉結(jié)構(gòu)[41],如圖8所示,該方法采用獨(dú)立產(chǎn)生的不共光路雙頻光源,其測(cè)量/參考光束在進(jìn)行外差干涉之前從不共光路,避免了雙頻激光無法完全分離的問題.該方法可以從源頭上消除雙頻激光交叉混疊引入的第1類周期非線性誤差.為此,美國(guó)NIST[42-43]、美國(guó)Rochester大學(xué)[44]、德國(guó)PTB[45]、荷蘭Delft理工大學(xué)[46]、我國(guó)清華大學(xué)[47]和哈爾濱工業(yè)大學(xué)[48]等設(shè)計(jì)了多種空間分離式外差激光干涉儀.這些空間分離式外差激光干涉儀原理上可完全消除第1類周期非線性誤差,實(shí)驗(yàn)表明剩余非線性誤差為數(shù)十皮米到一納米,并且剩余非線性誤差均為多階多普勒頻移虛反射光束自混疊引起的第2類周期非線性誤差[23-24].

        圖7 第1類周期非線性誤差補(bǔ)償方法

        圖8 空間分離式外差激光干涉結(jié)構(gòu)[23]

        2.2 第2類周期非線性誤差補(bǔ)償方法

        根據(jù)第2.1節(jié)的分析,空間分離式外差激光干涉結(jié)構(gòu)可以從原理上去除第1類周期非線性誤差.在此基礎(chǔ)上,研究人員提出了多種第2類周期非線性誤差的補(bǔ)償方法.根據(jù)第1.2節(jié)的分析,第2類周期非線性誤差與虛反射率的大小成正相關(guān)的關(guān)系.因此,鍍?cè)鐾改ひ越档凸鈱W(xué)元件表面反射率來減弱虛反射光束的強(qiáng)度的方法被提出來[49].該方法基于薄膜干涉原理,在光學(xué)元件表面鍍適當(dāng)厚度的薄膜,使在薄膜兩個(gè)面反射光的光程差恰好等于半個(gè)波長(zhǎng),從而相互抵消.在此基礎(chǔ)上,為了進(jìn)一步減低反射率,又提出了多層增透膜設(shè)計(jì)[50].普通寬帶增透膜反射率為0.25%,而針對(duì)單一波長(zhǎng)的窄帶增透膜反射率可達(dá)0.1%.根據(jù)第1.2節(jié)中式(14)可以計(jì)算出不同反射率下周期非線性誤差大小,見表1,可見僅僅通過降低反射率無法使周期非線性誤差降至皮米量級(jí).

        此外,根據(jù)虛反射光束方向與產(chǎn)生該光束的光學(xué)界面角度相關(guān)的特點(diǎn),可以通過調(diào)節(jié)虛反射面角度與空間濾波的方法來減小虛反射光束對(duì)外差激光干涉測(cè)量的影響[51-52].該方法通過調(diào)節(jié)虛反射面角度使虛反射光束與主要測(cè)量光束之間產(chǎn)生失配角,同時(shí)在接收器前加裝透鏡,利用透鏡的聚焦作用濾除與主要測(cè)量光束失配角過大的虛反射光束,從而消除虛反射光束引起的第2類周期非線性誤差.如圖9所示,空間濾波透鏡焦距為100 mm,光纖接收器直徑為10 μm,所有器件均鍍有增透膜,反射率約為0.1%.當(dāng)失配角大于200 μrad時(shí),該虛反射面產(chǎn)生的虛反射光束引入的周期非線性誤差基本被去除.德國(guó)聯(lián)邦物理研究院的Weichert-type空間分離式外差激光干涉結(jié)構(gòu)中采用了這些措施,該干涉結(jié)構(gòu)的周期非線性誤差能夠達(dá)到±10 pm[45].然而,由于外差激光干涉儀中若干光學(xué)界面之間虛反射行為復(fù)雜,存在著單次、多次虛反射情況,簡(jiǎn)單調(diào)節(jié)部分虛反射面角度不能保證所有虛反射光束完全逸出主要測(cè)量光路.因此,亟需研發(fā)能夠精確調(diào)整虛反射光束傳輸路徑的優(yōu)化和分析技術(shù)來消除虛反射光束引入的周期非線性誤差.

        表1 不同反射率下各階多普勒頻移虛反射光束引入的周期非線性誤差

        Tab.1 Nonlinear errors induced by DFS ghost beam of each order with different reflectivity

        類型界面反射率/%各階虛反射引入的非線性誤差/nm1階2階3階未鍍膜4.0010.1002.0000.400寬帶增透膜0.252.5000.1300.006窄帶增透膜0.101.6000.0500.001

        圖9 調(diào)節(jié)虛反射角度與空間濾波去除第2類周期非線性誤差原理與誤差曲線

        綜上所述,第1類周期非線性誤差與第2類周期非線性誤差的補(bǔ)償方法各不相同.針對(duì)第1類周期非線性誤差,現(xiàn)有的補(bǔ)償方法可將誤差抑制至0.1 nm量級(jí),近年來的空間分離式外差激光干涉結(jié)構(gòu)可以從原理上消除第1類周期非線性誤差.針對(duì)第2類周期非線性誤差,通過降低虛反射率以及調(diào)節(jié)虛反射面角度和空間濾波的方法,可以將非線性誤差減小至數(shù)十皮米.然而由于外差激光干涉儀中虛反射情況復(fù)雜,目前尚無針對(duì)性的優(yōu)化和分析技術(shù),不能完全滿足皮米測(cè)量的需求.

        3 結(jié) 論

        1)為提高外差激光干涉儀測(cè)量精度,對(duì)外差激光干涉儀中周期非線性誤差的形成機(jī)理和補(bǔ)償方法進(jìn)行研究,對(duì)比分析了外差激光干涉儀中的兩類周期非線性誤差,即雙頻激光交叉混疊引起的第1類周期非線性誤差和多階多普勒頻移虛反射光束自混疊引起的第2類周期非線性誤差.前者周期非線性誤差中只包含一階和二階非線性誤差,后者除一階和二階非線性誤差外,還包含有高階非線性誤差.此外,前者幅值可達(dá)數(shù)納米到數(shù)十納米,后者幅值為數(shù)納米到數(shù)皮米.

        2)現(xiàn)有第1類周期非線性誤差補(bǔ)償方法,可將誤差抑制至0.1 nm量級(jí),特別是空間分離式外差激光干涉儀,從原理上完全消除了這一類誤差;對(duì)于第2類周期非線性誤差,現(xiàn)有方法可以將其減小到數(shù)十皮米.

        3)隨著外差激光干涉測(cè)量精度需求逐漸由納米量級(jí)到皮米量級(jí),現(xiàn)有非線性誤差補(bǔ)償方法尚不能完全滿足皮米測(cè)量的需求,亟需研發(fā)新的第2類周期非線性誤差優(yōu)化補(bǔ)償技術(shù).

        猜你喜歡
        外差干涉儀光束
        2維Airy光束陣列強(qiáng)度的調(diào)控技術(shù)研究
        激光外差探測(cè)陣列互相關(guān)相干累積法
        基于改進(jìn)的邁克爾遜干涉儀對(duì)熱變形特性的研究
        基于結(jié)構(gòu)光三維重建系統(tǒng)的改進(jìn)相位研究*
        詭異的UFO光束
        奧秘(2021年3期)2021-04-12 15:10:26
        用于原子干涉儀的光學(xué)鎖相環(huán)系統(tǒng)
        非對(duì)稱干涉儀技術(shù)及工程實(shí)現(xiàn)
        基于外差-分?jǐn)?shù)階傅里葉變換的線性調(diào)頻連續(xù)波主動(dòng)聲吶處理
        基于最優(yōu)模糊的均勻圓陣干涉儀測(cè)向算法
        激光共焦顯微光束的偏轉(zhuǎn)掃描
        色偷偷亚洲第一成人综合网址| 亚洲av无码一区二区三区鸳鸯影院| av在线免费观看你懂的| 女人扒开下面无遮挡| 国产无套乱子伦精彩是白视频 | 亚洲国产精品二区三区| 日韩区在线| 久久精品中文字幕免费| 91福利国产在线观看一区二区| 日韩毛片免费无码无毒视频观看| 激情精品一区二区三区| 综合人妻久久一区二区精品| 亞洲綜合一區二區三區無碼| 一性一交一口添一摸视频| 亚洲线精品一区二区三区| 精品国产3p一区二区三区| 无码中文av有码中文av| 中出人妻中文字幕无码| 少妇下面好爽好紧好湿一区二区| 国产免费一区二区三区最新不卡| 国产精品久久久久久久免费看 | 亚洲精品中文字幕一二三 | a毛片全部免费播放| 国产无遮挡裸体免费视频| 爆乳熟妇一区二区三区霸乳 | 亚洲美女主播一区二区| 亚洲国产成人精品激情资源9| 乌克兰少妇xxxx做受6| 日本人妻免费在线播放| 99热婷婷一区二区三区| 亚洲欧洲中文日韩久久av乱码| 国产一区二区三区在线电影| 免费观看a级毛片| 国产人妖网站在线视频| 成人国产高清av一区二区三区| av网址不卡免费在线观看| 人妻无码人妻有码不卡| 国产综合自拍| 可以免费观看的毛片| 亚洲爱婷婷色婷婷五月| 白又丰满大屁股bbbbb|