999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

多普勒非對稱空間外差干涉儀調(diào)制度分析

2018-11-12 04:00:12施海亮李志偉羅海燕熊偉
航天返回與遙感 2018年5期
關(guān)鍵詞:制度

施海亮 李志偉 羅海燕 熊偉

?

多普勒非對稱空間外差干涉儀調(diào)制度分析

施海亮 李志偉 羅海燕 熊偉

(中國科學(xué)院安徽光學(xué)精密機(jī)械研究所,中國科學(xué)院通用光學(xué)定標(biāo)與表征技術(shù)重點(diǎn)實(shí)驗(yàn)室,合肥 230031)

多普勒非對稱空間外差光譜儀是近年得以迅速發(fā)展的新型中高層大氣風(fēng)場探測儀器,具有視場展寬、大光通量、高光譜分辨和多譜線同時探測等特點(diǎn)。該儀器的干涉調(diào)制度是影響風(fēng)場反演精度的關(guān)鍵指標(biāo)。文章基于非對稱空間外差光譜技術(shù)基本原理,推導(dǎo)了干涉儀視場展寬前后的干涉理論積分公式,分析了干涉調(diào)制度與儀器視場角、最優(yōu)光程差的關(guān)系,通過對干涉方程中像元采樣積分推導(dǎo)得到干涉調(diào)制度與采樣間隔的理論形式,仿真了干涉調(diào)制度與儀器視場角、最優(yōu)光程差和采樣間隔的定量關(guān)系曲線,并通過實(shí)驗(yàn)室搭建試驗(yàn)裝置對分析結(jié)果進(jìn)行驗(yàn)證。研究結(jié)果表明:隨著非對稱空間外差干涉儀視場角和光程差偏置量的增加,干涉調(diào)制度降低;調(diào)制度沿著波數(shù)增大的方向上升,在同一波數(shù)處,隨著采樣間隔的增加,干涉調(diào)制度降低。理論分析和試驗(yàn)測量結(jié)果調(diào)制度最大偏差0.02,分析與實(shí)測結(jié)果一致性較好。

多普勒外差干涉儀 調(diào)制度 視場角 最優(yōu)光程差 相位頻移 空間大氣探測

0 引言

中高層大氣對于研究全球氣候變化、大氣物理模型和航空航天活動具有重要的指導(dǎo)意義[1-4],其中風(fēng)場信息表征了大氣動態(tài)變化,是獲取大氣行為的關(guān)鍵因子。多普勒非對稱空間外差光譜技術(shù)(Doppler Asymmetric Spatial Heterodyne Spectroscopy,DASHS)作為新興的被動式風(fēng)場探測手段成為近些年的研究熱點(diǎn)。該技術(shù)2006年由J. M. Harlander等人提出,基于對稱空間外差光譜技術(shù)(Spatial Heterodyne Spectroscopy,SHS)[5]而改進(jìn),綜合空間干涉與光柵衍射于一體,不僅繼承了傳統(tǒng)邁克爾遜干涉光譜儀的優(yōu)點(diǎn),避免了運(yùn)動部件和窄帶濾波,而且能夠同時測量多個譜線,具備系統(tǒng)誤差同步探測的特點(diǎn),具有較高的風(fēng)場探測精度以及較強(qiáng)的空間環(huán)境適應(yīng)性[6]。

采用非對稱空間外差技術(shù)進(jìn)行風(fēng)場測量是依據(jù)多普勒以及傅里葉變換原理,通過對干涉圖進(jìn)行傅里葉變換獲取氣輝目標(biāo)光譜,然后逆變換獲取復(fù)數(shù)形式的干涉圖,由干涉圖的實(shí)虛部比求出光程差范圍內(nèi)每個點(diǎn)的相位,最終通過相位信息的變化量反演得到風(fēng)速[7-9]。干涉調(diào)制度是非對稱空間外差干涉儀進(jìn)行風(fēng)場探測過程中需要考慮的一個關(guān)鍵因素。調(diào)制度的大小直接影響了干涉條紋的品質(zhì),進(jìn)而經(jīng)過目標(biāo)光譜獲取和相位頻移信號提取后,導(dǎo)致風(fēng)速反演精度降低。影響調(diào)制度的因素是多方面的,具體包括干涉儀的視場角、最優(yōu)光程差和采樣間隔等。其中,由于擴(kuò)視場前后的非對稱空間外差干涉方程的差異,調(diào)制度與視場角的關(guān)系也并不相同。

基于以上考慮,本文從非對稱空間外差技術(shù)的基本原理出發(fā),分析并推導(dǎo)了與視場角和采樣間隔相關(guān)的積分干涉方程,得到干涉調(diào)制度公式;定量分析了儀器視場角、光程差偏置量和采樣間隔等系統(tǒng)參數(shù)對干涉調(diào)制度的影響;進(jìn)而用于指導(dǎo)多普勒非對稱空間外差干涉儀風(fēng)場探測儀器的方案設(shè)計和工程實(shí)現(xiàn)。

1 非對稱空間外差光譜技術(shù)原理

由于風(fēng)速引起的多普勒頻移量非常小,傳統(tǒng)光譜儀或干涉儀難以探測到如此微小的光譜頻移。在傳統(tǒng)的空間外差干涉儀的一臂上通過增加光程差偏置,可以增大由風(fēng)速引起的干涉相位變化量,探測此相位變化可反演獲得風(fēng)速。然而,考慮到干涉調(diào)制度會隨光程差的增加而減小,因此多普勒外差干涉儀存在一個最優(yōu)光程差opt[6-7]

單色光照明理想條件下,干涉圖調(diào)制度為1,調(diào)制度定義為

式中max為干涉條紋最大值;min為干涉條紋最小值。

2 儀器視場角與干涉調(diào)制度

為了獲得一定的輸入輻射強(qiáng)度,實(shí)際均采用有特定尺寸的擴(kuò)展光源,對儀器有一定的立體角。對于某一單色光p,擴(kuò)展光源將引起譜線展寬,使干涉條紋對比度降低,降低分辨率。為了滿足儀器的探測需求,干涉圖要有較高的調(diào)制度,DASHS干涉方程可以直接由空間外差干涉方程推導(dǎo)得到[11-14],這是因?yàn)楦缮娣匠淌菍αⅢw角積分,而與光程差無關(guān)。在非擴(kuò)視場的情況下,當(dāng)立體角為時,干涉條紋強(qiáng)度可以寫成如下形式

對方程(4)的立體角從0到s進(jìn)行積分得

式中為光柵有效截面長度。

式(5)積分號前面的s代表由于有限視場角而引起的條紋強(qiáng)度的增加值,非擴(kuò)視場的干涉調(diào)制函數(shù)可以寫為

在實(shí)際應(yīng)用中,為了增加信噪比,在非對稱空間外差干涉儀中通過擴(kuò)視場技術(shù)增大儀器視場。擴(kuò)視場后,上述干涉方程不再適用,根據(jù)光柵方程和擴(kuò)視場棱鏡出射光角度計算推導(dǎo)出空間外差干涉儀的擴(kuò)視場方程[15-17],進(jìn)而得到非對稱空間外差干涉儀的擴(kuò)視場方程

式中表示擴(kuò)視場棱鏡的折射率;表示光束經(jīng)過擴(kuò)視場棱鏡后的出射角度;表示擴(kuò)視場棱鏡的頂角。

對立體角從0到s進(jìn)行積分,根據(jù)方程(7)可得出

則擴(kuò)展視場DASHS干涉儀的調(diào)制度函數(shù)為

根據(jù)極限視場角的定義[5],可計算DASHS干涉儀的極限視場半角。

表1 DASHS干涉儀視場角參數(shù)

Tab.1 Angle of field of viewof the DASHS interferometer

根據(jù)表1計算的極限視場半角,在DASHS干涉儀擴(kuò)展視場前,視場角小于極限視場半角時的調(diào)制度曲線,如圖2所示。隨著視場半角的增加,干涉調(diào)制度逐漸降低,當(dāng)?shù)竭_(dá)極限視場角0.14°時,干涉調(diào)制度約0.64。擴(kuò)展視場后,極限視場增大,可以看出,視場角為0.91°時,理論上的干涉調(diào)制度才降至0.64,而在0.14°的地方,干涉調(diào)制度為0.999 79。分析結(jié)果表明,干涉儀擴(kuò)展視場后,調(diào)制度得到了極大的改善。

3 光程差偏置量與干涉調(diào)制度

根據(jù)調(diào)制度方程,干涉調(diào)制度不僅是視場角的函數(shù),也是光程差偏置量的函數(shù)。下面將對不同光程差偏置量off下的調(diào)制度(干涉儀擴(kuò)視場條件下)對比,如圖3所示。

由圖3可知,對于擴(kuò)展視場的DASHS干涉儀,其干涉調(diào)制度不僅隨著視場半角的增加而減小,而且隨著光程差偏置量off的增加而減小,因此必須選擇合適的光程差才能保證滿足調(diào)制度的需要。

4 數(shù)據(jù)采樣間隔與干涉調(diào)制度

方程(1)給出的干涉圖方程是在滿足采樣定理的條件下,將探測器上每個像元尺寸近似為無限小的條件下成立。事實(shí)上,探測器的像元實(shí)際都具有一定的尺寸,從而在一定程度影響了干涉調(diào)制度[18-20]。設(shè)探測器上單像元的尺寸為,在光譜維兩個相鄰的像元坐標(biāo)記為1,2,1像元中心點(diǎn)到2像元中心點(diǎn)之間某一位置的長度為m,則有

對積分項(xiàng)展開并計算可得

式中為探測器半尺寸。

因此,與像元尺寸有關(guān)的DASHS干涉儀調(diào)制度可以寫為

已知CCD探測器像元尺寸為13 μm,Littrow波長為764.78 nm,Littrow角為13.2638°,帶入得到調(diào)制度M與波數(shù)的關(guān)系曲線,見圖4。

由圖4看出,調(diào)制度沿著波數(shù)增大的方向上升,即在長波方向干涉調(diào)制度較低,在短波方向干涉調(diào)制度相對較高。另外,在同一波長處,隨著像元尺寸的增加,干涉調(diào)制度降低,圖中所給的像元尺寸最大為15μm。因此,滿足采樣定理前提下,為了減小探測器像元尺寸對調(diào)制度的影響,像元尺寸越小越好。

通過以上分析可知,影響DASHS干涉儀干涉調(diào)制度的因素是多方面的。由于干涉儀具有一定的視場角,因此通過對干涉方程中的立體角積分,可以得出干涉調(diào)制度隨著視場角和光程差的增加而降低;通過對干涉方程中的探測器面上像元積分,可以得出干涉調(diào)制度隨著像元尺寸的增加而降低,并且在短波方向干涉調(diào)制度較高。

5 試驗(yàn)結(jié)果與分析

中科院安徽光機(jī)所已成功搭建了基于DASHS原理的風(fēng)場探測試驗(yàn)裝置,針對760nm氧氣氣輝譜線進(jìn)行風(fēng)速探測[10],根據(jù)儀器視場角和干涉調(diào)制度的分析結(jié)果,為了保證其中干涉儀采用擴(kuò)視場一體化膠合設(shè)計,其中偏置光柵未進(jìn)行最后膠合,通過移動光柵位置采集不同偏置量下的干涉數(shù)據(jù)。試驗(yàn)裝置主要器件參數(shù)及系統(tǒng)的主要性能參數(shù)如表2所示。

試驗(yàn)采用鉀元素空心陰極燈作為光源,在波段內(nèi)具有兩條波長分別為766.49nm和769.90nm的特征發(fā)射線。光源發(fā)出的光經(jīng)過準(zhǔn)直鏡頭后進(jìn)入干涉儀,形成Fizeau條紋。在干涉儀出射端,干涉條紋被成像鏡頭按照1.6∶1的縮放比進(jìn)行縮放,最終成像于CCD探測器面,移動偏置臂光柵距離分別為10、30、50、60、70mm,測量實(shí)測不同偏置量下的干涉圖,如圖5所示。

對原始干涉圖扣除暗電流、去除探測器非均勻性校正和非線性校正影響、并且進(jìn)行干涉圖像相位校正后,計算干涉條紋的調(diào)制度,獲得不同偏置量下實(shí)際干涉數(shù)據(jù)調(diào)制度,與相同偏置量下理論干涉圖調(diào)制度對比,結(jié)果如表3所示。根據(jù)上文分析,獲取視場半角0.6°對應(yīng)不同偏置量的調(diào)制度因子、探測器像元13μm對應(yīng)的調(diào)制度曲線均值及成像鏡頭MTF,計算不同偏置量理論干涉圖調(diào)制度。

表2 多普勒非對稱空間外差光譜儀主要器件參數(shù)

Tab.2 The main parameters of DASHS breadboard

表3 不同光柵偏置量調(diào)制度理論與實(shí)測調(diào)制度

Tab.3 The theoretical modulation efficiency and the measured with different grating offsets

由表3對比分析可知,理論獲得的干涉圖調(diào)制度與實(shí)測干涉圖調(diào)制度一致性較好。需要指出的是,表3中理論調(diào)制度計算時探測器像元大小和成像鏡頭MTF的影響,均是以濾光片在整個波段內(nèi)的均值,理論調(diào)制度與實(shí)測干涉圖調(diào)制度的微小偏差是由于探測器的噪聲、計算模型的簡化、成像鏡頭的MTF誤差等因素造成。

6 結(jié)論

本文在空間外差干涉方程的基礎(chǔ)上推導(dǎo)了擴(kuò)視場前后的非對稱空間外差干涉方程,得出干涉調(diào)制度與儀器視場角和最優(yōu)光程差的定量關(guān)系,并對有限尺寸的探測器像元推導(dǎo)了像元尺寸與調(diào)制度的關(guān)系,首次從理論上分析了多普勒非對稱空間外差干涉儀調(diào)制度的影響因素,在實(shí)驗(yàn)室搭建試驗(yàn)裝置采集不同光柵偏置量下的干涉數(shù)據(jù),對實(shí)測干涉圖調(diào)制度和理論調(diào)制度進(jìn)行了對比分析,得到以下結(jié)論:

1)隨著視場角的增加,干涉調(diào)制度降低;

2)對比擴(kuò)展視場前后的調(diào)制度曲線得出,視場展寬有效提高了干涉調(diào)制度;

3)對比同一視場下的干涉調(diào)制度得出,隨著光程差偏置量的增加,干涉調(diào)制度降低;

4)隨著像元尺寸的增加,調(diào)制度降低,且調(diào)制度沿著波數(shù)增大的方向上升。

通過定量分析儀器視場角、光程差偏置量與探測器像元尺寸對多普勒空間外差干涉儀調(diào)制度的影響,可保證多普勒非對稱空間外差干涉通過優(yōu)選參數(shù)提升調(diào)制度,獲取更有效的干涉頻移信號。該結(jié)論可為儀器系統(tǒng)設(shè)計參數(shù)選擇提供重要理論依據(jù),對儀器方案優(yōu)化及工程研制具有重要的指導(dǎo)意義。

[1] SHEPHERD G G, THUILLIER G, GAULT W A, et al. WINDII, the Wind Imaging Interferometer on the Upper Atmosphere Research Satellite[J]. Journal of Geophysical Research: Atmospheres (1984-2012), 1993, 98(D6): 10725-10750.

[2] HAYS P B, ABREU V J, DOBBS M E, et al. The High-resolution Doppler Imager on the Upper Atmosphere Research Satellite[J]. Journal of Geophysical Research, 1993, 98(D6): 10713-10723.

[3] KILLEEN T L, SKINNER W R, JOHNSON R M, et al. TIMED Doppler Interferometer (TIDI)[J]. Optical Spectroscopic Techniques and Instrumentation for Atmospheric and Space Research III, 1999, 3756: 289-301.

[4] MCLANDRESS C, SHEPHERD G G, SOLHEIM B H, et al. Combined Mesosphere/ Thermosphere Winds Using WINDII and HRDI Data from the Upper Atmosphere Research Satellite[J]. Journal of Geophysical Research, 1996, 101: 10441-10453.

[5] HARLANDER J M. Spatial Heterodyne Spectroscopy: Interferometric Performance at Any Wavelength Without Scanning[D]. Wisconain, USA: The University of Wisconain, 1991.

[6] ENGLERT C R, HARLANDER J M, BABCOCK D D, et al. Doppler Asymmetric Spatial Heterodyne Spectroscopy (DASH: An Innovative Concept for Measuring Winds in Planetary Atmospheres[J]. SPIE Atmospheric Optical Modeling, Measurement, and Simulation II Optics+ Photonics, 2006, 6303: 290-296.

[7] ENGLERT C R, BABCOCK D D, HARLANDER J M. Doppler Asymmetric Spatial Heterodyne Spectroscopy (DASH): Concept and Experimental Demonstration[J]. Applied optics, 2007, 46(29): 7297-7307.

[8] ENGLERT C R, HARLANDER J M, EMMERT J T, et al. Initial Ground-based Thermospheric Wind Measurements Using Doppler Asymmetric Spatial Heterodyne Spectroscopy (DASH)[J]. Optics Express, 2010, 18(26): 27416-27430.

[9] HARLANDER J M, ENGLERT C R, BABCOCK D D, et al. Design and Laboratory Tests of a Doppler Asymmetric Spatial Heterodyne (DASH) Interferometer for Upper Atmospheric Wind and Temperature Observations[J]. Optics Express, 2010, 18(25): 26430-26440.

[10] 沈靜, 熊偉, 施海亮, 等. 用于風(fēng)場探測的多普勒非對稱空間外差干涉數(shù)據(jù)處理方法研究[J]. 光譜學(xué)與光譜分析, 2016, 36(9): 3014-3019. SHEN Jing, XIONG Wei, SHI Hailiang, et al. Doppler Asymmetric Spatial Heterodyne Interference Data Processing Methods Used in Wind Field Detection Research[J]. Spectroscopy and Spectral Analysis, 2016, 36 (9): 3014-3019. (in Chinese)

[11] 沈靜. 中高層大氣風(fēng)場探測多普勒非對稱空間外差技術(shù)研究[D]. 合肥: 中國科學(xué)技術(shù)大學(xué), 2016. SHEN Jing. Doppler's Asymmetric Spatial Heterodyne Technology for Mid and Upper Level Wind Field Detection[D]. Hefei: University of Science & Technology China, 2016. (in Chinese)

[12] SHEN J, XIONG W, SHI H L. et al. Data Processing Error Analysis Based on Doppler Asymmetric Spatial Heterodyne (DASH) measurement[J]. Applied Optics, 2017, 56 (12): 3531-3537.

[13] 李志偉, 熊偉, 施海亮, 等. 非對稱空間外差光譜技術(shù)研究[J]. 光譜學(xué)與光譜分析, 2016, 36(7): 2291-2295. LI Zhiwei, XIONG Wei, SHI Hailiang, et al. Research on Asymmetric Spatial Heterodyne Spectroscopy[J]. Spectroscopy and Spectral Analysis, 2016, 36(7): 2291-2295. (in Chinese)

[14] SHEN Jing, XIONG Wei, SHI Hailiang, et al. Phase Detection and Drift Correction for Doppler Asymmetric Spatial Heterodyne Interferometer[J]. Acta Photomica Sinica, 2017, 46(9): 114-120.

[15] 葉松. 空間外差光譜技術(shù)研究[D]. 北京: 中國科學(xué)院研究生院, 2006. YE Song. Spatial Heterodyne Spectroscopy Technology[D]. Beijing: Graduate University of Chinese Academy of Sciences, 2006. (in Chinese)

[16] 施海亮.空間外差光譜儀定標(biāo)技術(shù)研究[D]. 北京: 中國科學(xué)院研究生院, 2012. SHI Hailiang. Spatial Heterodyne Spectrometer Calibration Technology[D]. Beijing: Graduate University of Chinese Academy of Sciences, 2012. (in Chinese)

[17] 羅海燕. 空間外差干涉成像光譜技術(shù)研究[D]. 北京: 中國科學(xué)院大學(xué), 2015. LUO Haiyan. Spatial Heterodyne Interferometric Imaging Spectroscopy Technology[D]. Beijing: University of Chinese Academy of Sciences, 2015. (in Chinese)

[18] 陶然, 相里斌, 杜述松, 等. 高通量傅里葉變換成像光譜儀調(diào)制度分析[J]. 光子學(xué)報, 2006, 35(5): 675-679. TAO Ran, XIANG Libin, DU Shusong, et al. High Throughput Fourier Transform Imaging Spectrometer Modulation Analysis[J]. Acta Photorica Sinica, 2006, 35(5): 675-679. (in Chinese)

[19] 楊建峰, 相里斌, 高瞻, 等. 干涉成象光譜儀中前置光學(xué)系統(tǒng)遠(yuǎn)心角對干涉調(diào)制度的影響[J]. 光子學(xué)報, 1997, 26(4): 368-372. YANG Jianfeng, XIANG Libin, GAO Zhan, et al. Influence of the Far Optical Angle of the Front Optical System on Interferometric Modulation in Interferometric Imaging Spectrometer[J]. Acta Photorica Sinica, 1997, 26(4): 368-372. (in Chinese)

[20] 張淳民, 相里斌, 楊建峰, 等. 穩(wěn)態(tài)大視場偏振干涉成像光譜儀原理及調(diào)制度分析[J]. 西安交通大學(xué)學(xué)報, 2000(6): 93-97. ZHANG Chunmin, XIANG Libin, YANG Jianfeng, et al. The Principle and Modulation Analysis of the Steady-state Large Field Polarization Interference Imaging Spectrometer[J]. Journal of Xi'an Jiaotong University, 2000(6): 93-97. (in Chinese)

The Modulation Efficiency Analysis of Doppler Asymmetric Spatial Heterodyne Spectrometer

SHI Hailiang LI Zhiwei LUO Haiyan XIONG Wei

(Anhui Institute of Optics and Fine Mechanics, Chinese Academy of Sciences, Key Laboratory of Optical Calibration and Characterization, Hefei 230031, China)

As a new type of the detection instrument for upper atmospheric wind field, Doppler asymmetric spatial heterodyne spectrometer has been rapidly developed in recent years, with many simultaneous characteristics such as broadening view, large flux, high spectral resolution and multispectral detection ability. The interference modulation efficiency is the key parameter affecting the wind field retrieval precision. Based on the basic principle of asymmetric spatial heterodyne spectroscopy, the integral formula in the interference theory is derived before and after enlarging the view, and the relationship among the interference modulation efficiency, angle of field of view and the optimal optical path difference is analyzed. The theoretical expressions of the interference modulation efficiency and the sampling interval are derived by the pixel sampling integration in the interference equation, and the quantitative relationship curve among the interference modulation efficiency, the angle of field of view and the optimal optical path difference is simulated, and then validated by the laboratory breadboard. The results show that the interference modulation efficiency decreases while increasing the viewing angle and the offset of the optical path difference, and increases with the rise of the wave number or adding the sampling interval at same wave number. The derived modulation has good agreement with the theoretical calculation results, with the maximum error less than 0.02.

Doppler heterodyne interferometer; modulation efficiency; field of view; optimal path difference; phase shift; space atmosphere detection

O433.4

A

1009-8518(2018)05-0057-09

10.3969/j.issn.1009-8518.2018.05.008

施海亮,男,1983年生,2012年獲中國科學(xué)院研究生院光學(xué)專業(yè)博士學(xué)位,副研究員。研究方向?yàn)楦吖庾V遙感探測技術(shù)。E-mail:hlshi@aiofm.ac.cn。

2018-03-12

遙感科學(xué)國家重點(diǎn)實(shí)驗(yàn)室開放基金(Y7Y00100KZ);國家自然科學(xué)基金(41301373)

(編輯:王麗霞)

猜你喜歡
制度
淺探遼代捺缽制度及其形成與層次
重大誤解制度“重大”之認(rèn)定
刑事申訴制度的三重檢討與完善進(jìn)路
法大研究生(2018年2期)2018-09-23 02:20:40
完善我國人大制度的幾點(diǎn)思考
簽約制度怎么落到實(shí)處
構(gòu)建好制度 織牢保障網(wǎng)
一項(xiàng)完善中的制度
論讓與擔(dān)保制度在我國的立法選擇
翻譯人員參與刑事訴訟制度的構(gòu)建
應(yīng)訴管轄制度適用之探討
主站蜘蛛池模板: 午夜国产精品视频| 日韩精品无码一级毛片免费| 国产一区成人| 精品国产乱码久久久久久一区二区 | 国产精品私拍在线爆乳| 成年A级毛片| 午夜啪啪网| 色综合日本| 国产男女免费完整版视频| jizz在线免费播放| 九九久久99精品| 亚洲第一黄片大全| 国产成人亚洲综合a∨婷婷| 久久久受www免费人成| 国产麻豆福利av在线播放| a级毛片在线免费| 乱系列中文字幕在线视频| 日韩免费毛片| 国产精品亚欧美一区二区三区| 亚洲人成网站色7799在线播放| 亚洲天堂精品视频| 亚洲人成亚洲精品| 亚洲欧洲AV一区二区三区| 秘书高跟黑色丝袜国产91在线| 日韩在线中文| 婷婷六月综合| 久久精品无码一区二区国产区| 精品国产女同疯狂摩擦2| 亚洲二三区| 日本道中文字幕久久一区| 婷婷六月天激情| 精品视频一区在线观看| 美女无遮挡免费网站| 国产男人的天堂| 国产精品30p| 国产在线第二页| 久久久久久国产精品mv| 欧美 国产 人人视频| 在线看片国产| 99国产精品免费观看视频| 国产在线91在线电影| 亚洲视频影院| 亚洲AⅤ无码日韩AV无码网站| 国产在线精彩视频论坛| 免费高清a毛片| 午夜国产理论| 国产二级毛片| 久久久久人妻一区精品| 无码中文字幕精品推荐| 黄色三级网站免费| 色偷偷一区| 久久午夜夜伦鲁鲁片无码免费| 免费播放毛片| 97se亚洲综合不卡| 国产精品青青| 国产精品网址在线观看你懂的| 亚洲国产中文综合专区在| 国产精品理论片| 国内精品小视频福利网址| 国产AV毛片| 日韩一区二区三免费高清| 日韩在线第三页| 精品视频91| 爽爽影院十八禁在线观看| 91探花国产综合在线精品| 精品一區二區久久久久久久網站| 日本午夜网站| 国产SUV精品一区二区6| 久久久久国产精品免费免费不卡| 5388国产亚洲欧美在线观看| 91精品国产自产在线老师啪l| 亚洲综合九九| 日本在线视频免费| 欧洲av毛片| 国产精品lululu在线观看| 亚洲精品第1页| 国产成+人+综合+亚洲欧美| AV片亚洲国产男人的天堂| 亚洲视频免| 欧美成人看片一区二区三区| 人人爱天天做夜夜爽| 国产精品欧美激情|