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

    ALOS PALSAR雷達(dá)影像InSAR數(shù)據(jù)處理中的基線和地形誤差分析*

    2012-11-14 13:47:46喬學(xué)軍金銀龍汪漢勝江利明
    關(guān)鍵詞:基線差分高程

    沈 強(qiáng) 喬學(xué)軍 金銀龍 汪漢勝 江利明

    (1)中國科學(xué)院測量與地球物理研究所大地測量與地球動(dòng)力學(xué)國家重點(diǎn)實(shí)驗(yàn)室,武漢 430077 2)中國地震局地震研究所,武漢 430071 3)武漢大學(xué)水資源與水電工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,武漢430072)

    ALOS PALSAR雷達(dá)影像InSAR數(shù)據(jù)處理中的基線和地形誤差分析*

    沈 強(qiáng)1)喬學(xué)軍2)金銀龍3)汪漢勝1)江利明1)

    (1)中國科學(xué)院測量與地球物理研究所大地測量與地球動(dòng)力學(xué)國家重點(diǎn)實(shí)驗(yàn)室,武漢 430077 2)中國地震局地震研究所,武漢 430071 3)武漢大學(xué)水資源與水電工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,武漢430072)

    從理論上分析了InSAR處理中的兩大誤差源——基線和地形對形變測量精度的影響,并建立了評價(jià)模型。在此基礎(chǔ)上,研究了上述兩種誤差源對ALOS PALSAR的影響,獲取了兩種誤差對形變影響的空間分布特征。結(jié)果表明,ALOS PALSAR的軌道精度為3~5 m,在現(xiàn)有的軌道精度情況下,形變測量精度較差,需對基線誤差進(jìn)行糾正;地形誤差的影響很小,主要誤差源于軌道誤差。在現(xiàn)有DEM精度條件下(SRTM DEM),利用ALOS PALSAR,要獲得優(yōu)于1cm的形變測量精度,基線長度最好小于400 m。

    ALOS PALSAR;InSAR;基線誤差;地形誤差;基線長

    1 引言

    InSAR技術(shù)成像特點(diǎn)使得該技術(shù)在地殼形變(如地震,火山運(yùn)動(dòng))監(jiān)測中具有獨(dú)特優(yōu)勢,其在衛(wèi)星視線向(Line Of Sight:LOS)的形變測量精度可高達(dá)3~10 mm[1,2]。但在形變測量應(yīng)用中,其結(jié)果的精度仍受到諸多因素的影響。這些因素主要包括衛(wèi)星軌道精度、地形精度、時(shí)空去相關(guān)以及其他因素,如大氣延遲,電離層等,而軌道精度是影響形變測量精度的主要因素。對于兩通(two-pass)差分測量,外部數(shù)字高程模型(DEM)的精度也是影響形變測量精度的另一重要因素。此前已有學(xué)者以The European Remote Sensing Satellite(ERS-1)為例,對軌道、地形和相位誤差對形變測量的影響進(jìn)行過研究[3];Zebker等人[4]研究說明基線誤差呈弱二次曲面影響。但這些研究都僅在理論情況下得出的,其結(jié)果僅具指導(dǎo)性,且不能獲得這些因素影響的空間分布特征。另外不同的衛(wèi)星傳感器,由于雷達(dá)頻率和軌道等參數(shù)的不同,這些誤差源對形變結(jié)果的影響是否存在一致性也是值得認(rèn)真研究的問題。ALOS PALSAR因其L波段波長較長,受地物影響較小等,在地殼運(yùn)動(dòng)、火山、滑坡等研究方面有著廣泛的應(yīng)用,特別在地震的同震形變測量方面應(yīng)用較多[5-7]。然而,由于該衛(wèi)星軌道精度較低,定量評估上述因素對該衛(wèi)星的InSAR測量精度的影響是評價(jià)ALOS PALSAR形變測量精度的關(guān)鍵。因此本文擬分別從理論和實(shí)踐兩方面對影響形變測量精度的因素(基線和地形)進(jìn)行深入探討,以獲取基線和地形誤差對ALOS PALSAR形變測量精度的影響情況。首先從理論上分析基線和地形誤差對形變測量的影響;其次利用汶川震區(qū)的ALOS PALSAR雷達(dá)影像,在建立評價(jià)模型的基礎(chǔ)上,對以上兩種誤差引起的InSAR形變測量誤差進(jìn)行估計(jì),確定無精密軌道支持下的ALOS PALSAR基線誤差和地形誤差對InSAR形變測量精度影響的空間分布特征,定量化評估兩種誤差對形變測量精度的影響。這些研究將是評價(jià)ALOS PALSAR形變測量精度的基礎(chǔ)。

    2 基線和地形誤差理論分析

    2.1 基線誤差

    由于衛(wèi)星軌道誤差是系統(tǒng)誤差[4],而基線是由兩不同衛(wèi)星軌道獲得的,因此基線誤差也具系統(tǒng)性,這種系統(tǒng)誤差可以通過一定的處理方法加以消除或者減弱。我們更加感興趣的是基線誤差引起的不同地面點(diǎn)形變量的誤差的變率,簡稱相對誤差[3],尤其是近距點(diǎn)和遠(yuǎn)距點(diǎn)誤差之差。圖1給出了重復(fù)軌道InSAR測量原理示意圖[8]。A1和A2分別代表兩不同衛(wèi)星位置,其中,B、Bx,By、B⊥、B∥分別是基線、水平基線、垂直基線、有效垂直基線、有效平行基線; ρ、δρ分別代表衛(wèi)星在A1點(diǎn)到地面測點(diǎn)的距離及不同的兩個(gè)位置到同一地面點(diǎn)的距離之差;θ為衛(wèi)星視角,α為軌道傾角。h地面高程,H為衛(wèi)星高度。

    圖1 重復(fù)軌道InSAR測量幾何關(guān)系Fig.1 Simplified interferometric SAR geometry

    由重復(fù)干涉測量的相位-距離轉(zhuǎn)換公式,則可得到相位差

    利用幾何關(guān)系,可得到地面高程h為:

    如果重復(fù)飛行期間發(fā)生地表形變,并設(shè)定地表形變量在衛(wèi)星視線向投影為Δρ(圖1),簡稱形變量,則重復(fù)測量的相位差φ應(yīng)該為:

    指示 Bx、By微分表達(dá)式??紤]到 Bx/ρ?1 000/ 850 000=0.001 2,略去微小項(xiàng)有:

    可以看出,基線誤差引起的地表形變測量誤差呈系統(tǒng)性,并且與衛(wèi)星視角相關(guān)。由于其系統(tǒng)性,評估其影響的空間分布特征更具現(xiàn)實(shí)意義。Bx、By誤差引起地表形變量相對誤差分別為:

    圖2 基線誤差對形變量Δρ的影響Fig.2 Influence of the baseline errors on surface deformation measurment

    為了模擬不同視角條件下基線誤差對形變結(jié)果的影響,先設(shè)ALOS PALSAR幾何參數(shù),距離ρ,波長λ,衛(wèi)星視角θ0和平臺(tái)高度H分別為ρ=844 000 m、λ=0.236 m、θ0=34.3°、H=698 000 m。視角的不同則對應(yīng)的測區(qū)范圍的變化,因此本文取θ分別為33°、34°和36°,以近距為基準(zhǔn),大約對應(yīng)于測區(qū)寬度L分別為10 km、25 km和50 km進(jìn)行模擬計(jì)算。圖2給出了基線誤差在不同測區(qū)范圍對形變量測量精度的影響。其中L1=10 km、L2=25 km、L3=50 km,可以看出Bx對Δρ影響是By影響的1.5倍。測區(qū)范圍越大則基線誤差對形變測量精度的影響越大。當(dāng)區(qū)域范圍為50km,水平和垂直基線誤差均為10 m時(shí),對Δρ的影響分別達(dá)到60 cm和40 cm。因此對于ALOS PALSAR要獲取1 cm的形變精度,則需要優(yōu)于10 cm的基線精度,也說明基線精度是制約ALOS PALSAR地殼形變監(jiān)測精度重要因素。

    2.2 地形精度

    對公式(4)、(6)分別求h的微分,有

    解方程得到

    可以看出,DEM高程誤差對形變測量精度的影響,正比于有效垂直基線B⊥,反比于距離ρ和衛(wèi)星視角θ。有效垂直基線越長,誤差影響越大;距離和角度越大,誤差影響越小。說明DEM高程誤差的影響也存在系統(tǒng)性。式(15)給出了DEM誤差引起形變誤差的相對誤差。

    其中,θN、θR分別代表近距和遠(yuǎn)距時(shí)的衛(wèi)星視角(θN=32°,θR=36°),B⊥為有效垂直基線。

    圖3 DEM高程及高差誤差對形變測量精度的影響Fig.3 Influence of height and height difference errors on surface deformation measurment

    圖3分別給出了不同有效基線條件下,外部DEM精度引起形變量的絕對誤差和相對誤差。其中設(shè)定有效垂直基線B⊥分別為B1=100 m,B2= 500 m,B3=1 000 m,ρ=844 000 m,測區(qū)寬度L=50 km。可以看出,有效垂直基線B⊥越長,地形誤差對形變測量精度影響越大,其中高程誤差對形變量測量精度的影響接近10倍于高程相對誤差對觀測精度的影響。因此在差分測量中需要選擇有效基線較短的像對。圖3給出的是可能出現(xiàn)的各種極限情況,在實(shí)際處理中采用的SRTM DEM作為兩通差分干涉測量的外部DEM,其絕對精度一般優(yōu)于16 m。可以得出在以上三種有效基線條件下,形變測量誤差分別為0.34 cm、1.68 cm、3.36 cm。高程相對誤差對形變測量精度的影響分別為0.035 cm、0.18 cm、0.36 cm。因此也進(jìn)一步說明了在當(dāng)前DEM精度情況下,要獲得1 cm的形變監(jiān)測精度,需要選擇有效垂直基線的長度最好小于300 m。

    3 基線和DEM誤差對形變測量的影響的空間分布

    根據(jù)理論上分別對基線誤差和DEM誤差引起形變測量精度分析可以得出,要獲得高精度的形變監(jiān)測結(jié)果,要求高精度基線,選擇合理的干涉像對(如基線長度,時(shí)間間隔等),在進(jìn)行兩通差分時(shí),還需選擇精度較高的外部DEM。以上獲得的結(jié)果只具有指導(dǎo)性。定量化評價(jià)兩種誤差對結(jié)果的影響,特別是這些誤差對結(jié)果影響的空間分布特征是ALOS PALSAR形變測量精度評定的關(guān)鍵。因此需要研究一套切實(shí)可行的方法來分析和評價(jià)兩種誤差影響的空間分布特征。

    由于ALOS PALSAR raw1.0格式的雷達(dá)影像未提供精密的軌道,因此在實(shí)際處理過程中,需要采用基線重估計(jì)處理方法來消除或減弱軌道誤差的影響[9-12]。本文也將利用該技術(shù)來估計(jì)基線誤差和DEM誤差對形變結(jié)果的影響。如果假設(shè)形變測量誤差δ△LOS可由式(16)表示,主要包括基線誤差δB、DEM誤差δDEM、相位測量誤差δφ以及其他誤差(如電離層、大氣延遲等)。

    分別獲取基線估計(jì)前后的經(jīng)解纏的形變影像,設(shè)定兩次解纏的其他參數(shù)一致,因此可以假設(shè)其他3個(gè)因素的影響在兩次解纏過程中的影響一致,再對基線估計(jì)前后的解纏影像相差處理,可消除其他3個(gè)因素的影響,差分圖像則為基線誤差引起的誤差。在進(jìn)行DEM誤差評定時(shí),則利用兩種不同高程的DEM,分別進(jìn)行兩次解纏處理,誤差項(xiàng)中除DEM影像,其他項(xiàng)在兩次解纏過程中的影響相同,因此對兩解纏影像相差處理,可獲得DEM誤差對形變量的影響。

    主要技術(shù)流程(圖4)是首先采用兩通差分干涉測量方法對像對進(jìn)行處理[5,8],然后利用外部DEM和解纏后的總相位對基線進(jìn)行重新估計(jì),再進(jìn)行兩通差分處理。分別獲得了基線重新估計(jì)前后形變影像。另外為了評估不同精度的DEM對結(jié)果的影響,本文利用SRTM DEM提取了兩種高程基準(zhǔn)的數(shù)字高程模型,一種的高程基準(zhǔn)是似大地水準(zhǔn)面,而另一種是參考橢球面。利用這兩種DEM分別進(jìn)行兩次差分處理,獲得基于不同DEM的形變分布圖。利用這些形變分布圖可分析和研究ALOS PALSAR的基線精度和DEM精度對形變測量精度的影響。

    圖4 兩通差分干涉測量流程及策略Fig.4 Flow chart and processing strategies of two-pass interferometry

    為了充分考慮不同基線、不同軌跡對形變結(jié)果的影響程度,本文基于以下4個(gè)標(biāo)準(zhǔn)選擇干涉雷達(dá)像對:1)選擇不同基線;2)影像對在重復(fù)測量期間無明顯的地表形變;3)選擇丘陵、平原地區(qū)的影像,減少地形及植被對相位測量精度的影響;4)為了充分評價(jià)不同軌道條件下,上述參數(shù)對結(jié)果的影響,選擇了不同軌道的3幅干涉像對進(jìn)行分析。干涉雷達(dá)像對基線參數(shù)見表1。

    表1 干涉像對的參數(shù)(單位:m)Tab.1 Parameters of interferometry(unit:m)

    4 結(jié)果分析

    由于選擇的干涉像對在成像期間地表未發(fā)生明顯形變,因此在進(jìn)行兩通差分后,獲得的形變量分布圖被認(rèn)為是誤差分布圖,主要包括相位測量誤差、基線誤差和DEM誤差及其他誤差(如電離層、大氣延遲等)。進(jìn)行基線重新估計(jì)后再進(jìn)行的差分干涉測量結(jié)果則主要包含的是相位測量誤差和DEM誤差及其他誤差。對同一外部DEM的兩次差分處理結(jié)果再進(jìn)行相差處理,則可獲得基線誤差的空間分布。對不同DEM條件下的解纏相位圖進(jìn)行相差處理,即可獲得DEM誤差對形變測量精度的影響。

    圖5分別給出了像對2(表1)干涉測量及基線估計(jì)后的形變相位分布圖,其中圖(a)、(b)、(c)分別為基線估計(jì)前的相位分布、基線估計(jì)后相位的分布和以上兩相位分布的差分結(jié)果。(e)、(f)、(g)分別對應(yīng)于(a)、(b)、(c)圖中表示的剖面線。單位是弧度,乘以λ/4π,即可換算成視線距離(LOS)。從圖中可以看出,基線估計(jì)前相位殘差分布呈現(xiàn)系統(tǒng)變化,由近距到遠(yuǎn)距分別從2.5~-6.0(圖(a)、(e))?;€糾正后,殘差分布呈現(xiàn)隨機(jī)分布特征。大小在-1~1之間(圖(b)、(f))。獲得的軌道誤差從-5~2變化。對圖(c)的統(tǒng)計(jì)發(fā)現(xiàn)最大值為-5.5,轉(zhuǎn)換到距離單位則為-10.3 cm。說明基線產(chǎn)生的最大誤差為10.3 cm。根據(jù)式(11)、(12)可得到最大基線誤差為2.57 m。同樣對另外兩像對進(jìn)行統(tǒng)計(jì)發(fā)現(xiàn),基線最大誤差分別為13.9 cm與17.8 cm,對應(yīng)最大基線誤差為3.47 m、4.25 m。通過對3個(gè)不同軌跡的ALOS PALSAR的基線精度的評價(jià)可以得出,ALOS PALSAR的軌道精度為3~5 m。也說明在無精密軌道支持下,ALOS PALSAR兩通InSAR測量結(jié)果的精度較差,因此必須對基線進(jìn)行重新估計(jì),消除或減弱軌道誤差對形變測量精度的影響。

    圖5 基線估計(jì)前后殘差分布Fig.5 Residual distribution before and after baseline estimation

    圖6 DEM誤差引起的形變測量誤差Fig.6 Residual distribution of terrain errors

    圖6給出了3組像對DEM誤差對形變測量的評價(jià)結(jié)果,其中圖6(e)、(f)、(g)分別對應(yīng)于圖6 (a)、(b)、(c)中線條所示位置的剖面線。圖6(e)、(f)、(g)左縱軸代表是DEM誤差引起的形變誤差。右縱軸代表對應(yīng)位置兩DEM的高差。DEM對像對1、2、3中影響分別為-0.28(-0.5 cm)、1.6(3.0 cm)、2.5(4.7 cm)。3幅影像覆蓋范圍的兩DEM高差變化為40~45 m。從圖中可以明顯看出,DEM誤差的影響差異較大,對照表1,可以看出有效垂直基線越長,則DEM誤差對結(jié)果的影響越大,因此要獲得1 cm的形變精度,則要求有效垂直基線小于150 m。對于現(xiàn)有的SRTM DEM,其高程精度優(yōu)于16 m,在基線小于400m的情況下,DEM引起的誤差小于1 cm。這也說明了差分測量中的DEM誤差對形變測量精度的敏感性較低。

    5 結(jié)語

    分別從理論和實(shí)際對影響InSAR形變測量因素的基線精度和 DEM精度進(jìn)行分析探討,得出ALOS PALSAR的軌道精度為3~5 m,要使InSAR形變測量結(jié)果具可量測性,必須對基線進(jìn)行糾正,消除或減少基線誤差對結(jié)果精度的影響。DEM誤差對形變量的影響根據(jù)有效基線的不同差異較大。對于現(xiàn)有的SRTM DEM情況下,有效垂直基線小于400m的像對時(shí),DEM誤差對形變測量結(jié)果的影響可以忽略。

    致謝 感謝Jaxa提供的ALOS PALSAR雷達(dá)影像!

    1 Massonnet D,et al.The displacement field of the Landers earthquake mapped by radar interferometry[J].Nature,1993,364:138-142.

    2 Klees R and Massonnet D.Deformation measurement using SAR interferometry:Potential and limitations[J].Geologie en Mijnbouw,1998,77(2):161-176.

    3 李振洪,劉經(jīng)南,許才軍.InSAR數(shù)據(jù)處理中的誤差分析[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2004,29(1):72-76.(Li Zhenhong,Liu Jingnan and Xu Caijun.Error analysis in InSAR data processing[J].Geomatics and Information Science of Wuhan University,2004,29(1):72-76)

    4 Zebker H A,et al.Accuracy of topographic maps derived from ERS-1 interferometric radar[J].IEEE Transactions on Geoscience and Remote Sensing,1994,32:823-836.

    5 Shen Zhengkang,et al.Slip maxima at fault junctions and rupturing of barriers during the 2008 Wenchuan earthquake[J].Nature Geoscience,2009,2:718-724.

    6 沈強(qiáng),譚凱.日本Mw9.0地震多傳感器InSAR形變初步分析[J].大地測量學(xué)與地球動(dòng)力學(xué),2011,(5):21-25.(Shen Qiang and Tan Kai.Characteristic analysis of coseismic deformation from InSAR images[J].Journal of Geodesy and Geodynamics,2011,(5):21-25)

    7 Manabu Hashimotoa,Mari Enomotoa and Yo Fukushima.Coseismic deformation from the 2008 Wenchuan,China,Earthquake Derived from ALOS/PALSAR Images[J].Tectonophysics,2010,491(1-4):59-71.

    8 舒寧.雷達(dá)影像干涉測量原理[M].武漢大學(xué)出版社,2003.(Shu Ning.Principles of rader interferometry[M].Wuhan:Press of Wuhan University,2003)

    9 A Monti Guarnieri,et al.Accurate and robust baseline estimation[A].Proceedings of ESA Fringe’99[C].Available on the WEB at the URL:http://home.dei.polimi.it/monti/papers/fringes99_mg.pdf.

    10 Kuldip S,et al.Baseline estimation in interferometric SAR[J].Geoscience and Remote Sensing,1997,1:454-456.

    11 Nico G.A baseline estimation algorithm for SAR interferometry[A].Proceedings of ESA Fringe’97,1997[C].A-vailable on the WEB at the URL:http://earth.esa.int/ workshops/ers97/papers/singh/.

    12 鄭芳,鄔鈞霆,馬德寶.基于軌道誤差的基線估計(jì)及其對高程精度的影響[J].信息工程大學(xué)學(xué)報(bào),2005,6(1):77 -79.(Zheng Fang,Wu Junting and Ma Debao.Baseline estimation based on the track error and the effect on height decision[J].Journal of Information Engineering University,2005,6(1):77-79)

    ERROR ANALYSIS OF BASELINE AND TERRAIN IN InSAR DATA PROCESSING USING ALOS PALSAR

    Shen Qiang1),Qiao Xuejun2),Jin Yinlong3),Wang Hansheng1)and Jiang Liming1)

    (1)State Key Laboratory of Geodesy and Earth’s Dynamics,Institute of Geodesy and Geophysics,CAS,Wuhan 430077 2)Institute of Seismology of China Earthquake Administration,Wuhan 430071 3)State Key Laboratory of Water Resources and Hydropower Engineering Science,Wuhan University,Wuhan 430072)

    Seeing that there are still many factors governing the accuracy of InSAR although they have appeared to good application in various field as earthquake,landslide,etc.this paper focus on the two sources of errors controlling the accuracy of InSAR.In the first place,theoretical analysis has been first made for the evaluation of influence of baseline and terrains errors on surface deformation.On the basis of above theoretical analysis,we estimated the influence of baseline errors and terrain errors on InSAR for ALOS PALSAR radar images without accurate orbit information.The spatial patterns of additional surface deformation due to the two error sources have also been retrieved.The results show that the orbit accuracy of ALOS PALSAR is 3-5 m,in which the displacements must be rectified for the low orbit accuracy.The simulation and experimentation all show that baseline is main error sources in the use of ALOS PALSAR and terrain errors is negligible.For the accuracy of better than 1 cm in deformation measurements,the length of baseline must be less than 400 m under current SRTM DEM if two-pass interferometry was used.

    ALOS PALSAR;InSAR;baseline errors;terrain errors;length of baseline

    1671-5942(2012)02-0001-06

    2011-12-21

    水資源與水電工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室開放基金(2009B054);國家自然科學(xué)基金(41004008,40825012);冰凍圈國家重點(diǎn)實(shí)驗(yàn)室開放基金(SKLCS 09-03)

    沈強(qiáng),1980年生,博士,副研究員,現(xiàn)從事衛(wèi)星遙感地殼形變研究.E-mail:cl980606@hotmail.com

    P227

    A

    猜你喜歡
    基線差分高程
    數(shù)列與差分
    適用于MAUV的變基線定位系統(tǒng)
    8848.86m珠峰新高程
    航天技術(shù)與甚長基線陣的結(jié)合探索
    科學(xué)(2020年5期)2020-11-26 08:19:14
    GPS控制網(wǎng)的高程異常擬合與應(yīng)用
    一種改進(jìn)的干涉儀測向基線設(shè)計(jì)方法
    基于差分隱私的大數(shù)據(jù)隱私保護(hù)
    SDCORS高程代替等級水準(zhǔn)測量的研究
    相對差分單項(xiàng)測距△DOR
    太空探索(2014年1期)2014-07-10 13:41:50
    回歸支持向量機(jī)在區(qū)域高程異常擬合中的應(yīng)用
    一本—道久久a久久精品蜜桃钙片| 久久精品熟女亚洲av麻豆精品| 午夜老司机福利剧场| 欧美人与性动交α欧美精品济南到 | 亚洲av成人精品一二三区| 高清黄色对白视频在线免费看| 青春草国产在线视频| 亚洲国产av新网站| av在线老鸭窝| 久久久久国产精品人妻一区二区| 日本爱情动作片www.在线观看| 亚洲欧洲国产日韩| 婷婷色av中文字幕| 国产 精品1| 成年女人毛片免费观看观看9 | 性少妇av在线| 亚洲精品aⅴ在线观看| 日韩电影二区| 婷婷色av中文字幕| av女优亚洲男人天堂| 在线观看免费日韩欧美大片| 国产一区二区三区综合在线观看| 久久毛片免费看一区二区三区| 欧美在线黄色| 久久久久国产精品人妻一区二区| 天天躁夜夜躁狠狠久久av| 午夜影院在线不卡| 亚洲成人一二三区av| 日本午夜av视频| 老汉色av国产亚洲站长工具| 韩国精品一区二区三区| 亚洲一区中文字幕在线| 国产精品女同一区二区软件| 久久精品国产自在天天线| 最近中文字幕2019免费版| 伊人久久大香线蕉亚洲五| 男女边吃奶边做爰视频| 国产免费福利视频在线观看| 久久亚洲国产成人精品v| 天天躁夜夜躁狠狠躁躁| 亚洲av在线观看美女高潮| 国产日韩一区二区三区精品不卡| 亚洲av成人精品一二三区| 一级,二级,三级黄色视频| 人人妻人人澡人人爽人人夜夜| 亚洲国产看品久久| 久久久久久人人人人人| 美女国产高潮福利片在线看| 9色porny在线观看| 七月丁香在线播放| 日本av免费视频播放| av.在线天堂| 久久精品国产亚洲av涩爱| 香蕉精品网在线| 精品少妇久久久久久888优播| 精品第一国产精品| 观看美女的网站| 黑人欧美特级aaaaaa片| 母亲3免费完整高清在线观看 | 国产成人欧美| 亚洲av免费高清在线观看| 精品国产乱码久久久久久小说| 一区二区日韩欧美中文字幕| 久久久久久人人人人人| 亚洲一区中文字幕在线| 久久久久久久精品精品| 亚洲欧美成人精品一区二区| 国语对白做爰xxxⅹ性视频网站| 国产精品偷伦视频观看了| 久久精品国产亚洲av高清一级| 日韩在线高清观看一区二区三区| 天美传媒精品一区二区| 日本色播在线视频| 美女大奶头黄色视频| 三上悠亚av全集在线观看| 777久久人妻少妇嫩草av网站| 国产在视频线精品| 视频在线观看一区二区三区| 午夜免费男女啪啪视频观看| 男女午夜视频在线观看| 日韩免费高清中文字幕av| 欧美精品av麻豆av| 在线观看免费日韩欧美大片| 母亲3免费完整高清在线观看 | 亚洲精品国产色婷婷电影| 国产女主播在线喷水免费视频网站| 亚洲三级黄色毛片| 最近中文字幕高清免费大全6| 91午夜精品亚洲一区二区三区| 国产深夜福利视频在线观看| 日日爽夜夜爽网站| 18禁裸乳无遮挡动漫免费视频| 亚洲三级黄色毛片| 久久午夜综合久久蜜桃| 另类亚洲欧美激情| 婷婷成人精品国产| 极品人妻少妇av视频| 亚洲精品久久成人aⅴ小说| 日本色播在线视频| 综合色丁香网| 观看美女的网站| 成人国语在线视频| 免费高清在线观看日韩| 久久99一区二区三区| 美女大奶头黄色视频| 婷婷色麻豆天堂久久| 亚洲欧美一区二区三区黑人 | 爱豆传媒免费全集在线观看| 91在线精品国自产拍蜜月| 国产片特级美女逼逼视频| 亚洲国产看品久久| 精品国产乱码久久久久久男人| 欧美日韩视频精品一区| 啦啦啦在线观看免费高清www| 亚洲国产精品一区三区| 免费黄色在线免费观看| 国产亚洲av片在线观看秒播厂| 大香蕉久久网| 看非洲黑人一级黄片| 中文字幕最新亚洲高清| 国产精品久久久久久精品电影小说| 久久久久久人妻| 少妇被粗大的猛进出69影院| 国产精品.久久久| 国产精品成人在线| 久久久久人妻精品一区果冻| 建设人人有责人人尽责人人享有的| 韩国av在线不卡| 精品亚洲成国产av| 国产熟女欧美一区二区| 少妇人妻久久综合中文| 精品少妇内射三级| 亚洲成人一二三区av| 老汉色av国产亚洲站长工具| 欧美日韩成人在线一区二区| 两个人看的免费小视频| 麻豆乱淫一区二区| 亚洲视频免费观看视频| 久久99精品国语久久久| 国产人伦9x9x在线观看 | 国产福利在线免费观看视频| 一本—道久久a久久精品蜜桃钙片| 黄色毛片三级朝国网站| 三上悠亚av全集在线观看| 国产毛片在线视频| 亚洲av国产av综合av卡| 国产视频首页在线观看| 麻豆精品久久久久久蜜桃| 日韩视频在线欧美| 午夜福利网站1000一区二区三区| 亚洲一码二码三码区别大吗| 韩国精品一区二区三区| 欧美xxⅹ黑人| 国产亚洲精品第一综合不卡| 国产精品三级大全| 亚洲国产欧美在线一区| 国产精品秋霞免费鲁丝片| 久久99蜜桃精品久久| 国产精品一国产av| 美女中出高潮动态图| 欧美在线黄色| 国产麻豆69| 国产色婷婷99| 成年美女黄网站色视频大全免费| 乱人伦中国视频| 中国国产av一级| 自拍欧美九色日韩亚洲蝌蚪91| 26uuu在线亚洲综合色| 欧美亚洲 丝袜 人妻 在线| 成人影院久久| 精品亚洲乱码少妇综合久久| 成年人午夜在线观看视频| 免费久久久久久久精品成人欧美视频| 人人妻人人澡人人看| 亚洲人成电影观看| 欧美精品一区二区大全| 久久国产精品大桥未久av| 黄片小视频在线播放| 国产成人午夜福利电影在线观看| 成年女人在线观看亚洲视频| 亚洲图色成人| 黑丝袜美女国产一区| 在线观看免费日韩欧美大片| 美女xxoo啪啪120秒动态图| 香蕉精品网在线| 亚洲男人天堂网一区| 国产精品偷伦视频观看了| 成人影院久久| 国产欧美日韩一区二区三区在线| 国产免费现黄频在线看| 国产亚洲午夜精品一区二区久久| 国产老妇伦熟女老妇高清| 国产免费又黄又爽又色| 精品福利永久在线观看| 国产国语露脸激情在线看| 精品少妇久久久久久888优播| 妹子高潮喷水视频| 亚洲av电影在线进入| 日韩不卡一区二区三区视频在线| 免费在线观看视频国产中文字幕亚洲 | 国产激情久久老熟女| 9热在线视频观看99| 精品亚洲成a人片在线观看| 亚洲精品aⅴ在线观看| 中文精品一卡2卡3卡4更新| 免费在线观看视频国产中文字幕亚洲 | 久久亚洲国产成人精品v| 婷婷色麻豆天堂久久| 成人漫画全彩无遮挡| 激情视频va一区二区三区| www.av在线官网国产| 99香蕉大伊视频| 午夜福利视频精品| 亚洲国产看品久久| 亚洲天堂av无毛| 春色校园在线视频观看| 欧美在线黄色| 久久99精品国语久久久| 99九九在线精品视频| 久久久久久人人人人人| 国产女主播在线喷水免费视频网站| 亚洲欧洲日产国产| 超碰97精品在线观看| 国产一区二区三区综合在线观看| 日韩大片免费观看网站| 久久久久久人人人人人| 亚洲第一区二区三区不卡| 午夜激情久久久久久久| 亚洲精品第二区| 99re6热这里在线精品视频| 天天影视国产精品| 精品人妻偷拍中文字幕| 亚洲 欧美一区二区三区| 欧美精品亚洲一区二区| 色网站视频免费| 成年女人毛片免费观看观看9 | 五月开心婷婷网| 日韩制服骚丝袜av| 九色亚洲精品在线播放| av.在线天堂| 97在线视频观看| 热re99久久精品国产66热6| 看十八女毛片水多多多| 免费av中文字幕在线| 国产一级毛片在线| 亚洲av日韩在线播放| 欧美另类一区| 国产亚洲精品第一综合不卡| 亚洲一区二区三区欧美精品| 久久午夜福利片| 久久久久久久亚洲中文字幕| 最新的欧美精品一区二区| 一区二区三区四区激情视频| 日本欧美视频一区| 两性夫妻黄色片| 咕卡用的链子| 免费在线观看完整版高清| 亚洲欧洲国产日韩| 黄色视频在线播放观看不卡| 精品国产一区二区三区久久久樱花| 日韩av免费高清视频| 麻豆乱淫一区二区| 老司机影院毛片| 国产日韩一区二区三区精品不卡| 免费黄网站久久成人精品| 国产熟女欧美一区二区| 在线观看国产h片| 久久久久国产一级毛片高清牌| 建设人人有责人人尽责人人享有的| 国产欧美亚洲国产| 2018国产大陆天天弄谢| 久久久久久人人人人人| 国产熟女欧美一区二区| 中文字幕av电影在线播放| 新久久久久国产一级毛片| 久久99一区二区三区| 美女xxoo啪啪120秒动态图| 另类精品久久| 欧美日韩成人在线一区二区| 国产在视频线精品| 免费在线观看完整版高清| 欧美日韩精品成人综合77777| 久久久久国产精品人妻一区二区| 日韩,欧美,国产一区二区三区| 各种免费的搞黄视频| 少妇人妻精品综合一区二区| 免费高清在线观看日韩| 精品少妇内射三级| 久久精品aⅴ一区二区三区四区 | 男的添女的下面高潮视频| 老汉色av国产亚洲站长工具| 亚洲欧美日韩另类电影网站| 卡戴珊不雅视频在线播放| 日本猛色少妇xxxxx猛交久久| 丝袜人妻中文字幕| 国产在线视频一区二区| 久热这里只有精品99| 亚洲三级黄色毛片| 国产亚洲av片在线观看秒播厂| 国产精品99久久99久久久不卡 | 中文字幕精品免费在线观看视频| www.精华液| 欧美日韩视频精品一区| 国产精品 国内视频| 高清黄色对白视频在线免费看| 国产白丝娇喘喷水9色精品| 伊人久久大香线蕉亚洲五| 欧美 日韩 精品 国产| 亚洲国产精品国产精品| 久久国产精品男人的天堂亚洲| 丝袜脚勾引网站| 夜夜骑夜夜射夜夜干| 超碰成人久久| 亚洲成人手机| 久久 成人 亚洲| 久久 成人 亚洲| 午夜激情久久久久久久| 青春草亚洲视频在线观看| 国产精品 欧美亚洲| 夫妻午夜视频| 亚洲综合色网址| 欧美少妇被猛烈插入视频| 欧美精品国产亚洲| 亚洲精品国产一区二区精华液| 99热网站在线观看| 如日韩欧美国产精品一区二区三区| 1024视频免费在线观看| 中文字幕人妻丝袜制服| 宅男免费午夜| 亚洲国产精品999| 岛国毛片在线播放| 免费观看av网站的网址| 欧美日韩一区二区视频在线观看视频在线| 欧美激情高清一区二区三区 | 精品国产一区二区久久| 久久久久久久久久久久大奶| 一区福利在线观看| 男人舔女人的私密视频| 女的被弄到高潮叫床怎么办| 亚洲国产精品国产精品| 久久毛片免费看一区二区三区| 久久久国产精品麻豆| 婷婷成人精品国产| 老女人水多毛片| 国产精品不卡视频一区二区| 99久久综合免费| 日韩制服丝袜自拍偷拍| 成年女人在线观看亚洲视频| 18禁国产床啪视频网站| 国产精品香港三级国产av潘金莲 | 有码 亚洲区| 如日韩欧美国产精品一区二区三区| 日韩中字成人| 一级毛片 在线播放| 国产精品免费大片| 青春草亚洲视频在线观看| 久久人人爽av亚洲精品天堂| 卡戴珊不雅视频在线播放| 欧美日韩一级在线毛片| 日韩中字成人| 精品一区二区三区四区五区乱码 | av线在线观看网站| 成年av动漫网址| 亚洲精品中文字幕在线视频| 乱人伦中国视频| 韩国av在线不卡| 日本91视频免费播放| 欧美激情 高清一区二区三区| 看十八女毛片水多多多| 亚洲人成77777在线视频| 2022亚洲国产成人精品| 看非洲黑人一级黄片| 99热网站在线观看| 亚洲国产精品一区三区| 国产av国产精品国产| 女性被躁到高潮视频| 婷婷成人精品国产| 国产乱人偷精品视频| 精品人妻在线不人妻| 国产一区二区三区综合在线观看| 免费高清在线观看日韩| 亚洲国产精品一区二区三区在线| 永久网站在线| 成人手机av| 免费不卡的大黄色大毛片视频在线观看| 成人国语在线视频| 国产精品 欧美亚洲| 国产高清国产精品国产三级| 青青草视频在线视频观看| 国精品久久久久久国模美| 久久国产亚洲av麻豆专区| 午夜福利视频精品| 热99国产精品久久久久久7| 亚洲精品中文字幕在线视频| 91成人精品电影| 日本黄色日本黄色录像| 人体艺术视频欧美日本| 婷婷色麻豆天堂久久| 久久久久视频综合| 伦理电影免费视频| videossex国产| 亚洲精品一区蜜桃| 欧美激情高清一区二区三区 | 少妇人妻精品综合一区二区| 国产精品麻豆人妻色哟哟久久| 国产av精品麻豆| 成人18禁高潮啪啪吃奶动态图| 久久精品国产鲁丝片午夜精品| 国产精品 国内视频| av卡一久久| 丝袜喷水一区| 最新的欧美精品一区二区| 一二三四中文在线观看免费高清| 天天影视国产精品| 只有这里有精品99| 天天躁日日躁夜夜躁夜夜| 亚洲精华国产精华液的使用体验| 99香蕉大伊视频| 中文乱码字字幕精品一区二区三区| 国产午夜精品一二区理论片| 欧美日韩一级在线毛片| 欧美人与性动交α欧美软件| 精品视频人人做人人爽| 日韩av在线免费看完整版不卡| 国产精品 欧美亚洲| 亚洲av综合色区一区| 亚洲国产精品999| 又大又黄又爽视频免费| 亚洲国产欧美网| 丰满饥渴人妻一区二区三| 亚洲精品一二三| 亚洲精品日韩在线中文字幕| 久久ye,这里只有精品| a 毛片基地| 人妻 亚洲 视频| 国产精品成人在线| 黄色毛片三级朝国网站| 久久精品国产亚洲av天美| 黑丝袜美女国产一区| 亚洲综合色惰| 99久国产av精品国产电影| 涩涩av久久男人的天堂| 母亲3免费完整高清在线观看 | 亚洲av男天堂| 十分钟在线观看高清视频www| 男人爽女人下面视频在线观看| 欧美最新免费一区二区三区| 在线观看国产h片| 91成人精品电影| 国产精品无大码| 啦啦啦视频在线资源免费观看| 伊人久久大香线蕉亚洲五| 黄频高清免费视频| 亚洲精品乱久久久久久| 777久久人妻少妇嫩草av网站| 亚洲人成网站在线观看播放| 免费高清在线观看视频在线观看| 十八禁网站网址无遮挡| 中文字幕人妻熟女乱码| 久久国产亚洲av麻豆专区| 国产精品麻豆人妻色哟哟久久| 天天操日日干夜夜撸| 校园人妻丝袜中文字幕| a级片在线免费高清观看视频| 亚洲综合色网址| 香蕉精品网在线| 久久人妻熟女aⅴ| 伦理电影免费视频| 乱人伦中国视频| 亚洲精品aⅴ在线观看| 久久影院123| av女优亚洲男人天堂| 看十八女毛片水多多多| 免费久久久久久久精品成人欧美视频| 18禁动态无遮挡网站| 精品福利永久在线观看| 亚洲,一卡二卡三卡| 国产一区有黄有色的免费视频| 免费人妻精品一区二区三区视频| 亚洲婷婷狠狠爱综合网| 国产在线视频一区二区| av国产精品久久久久影院| 捣出白浆h1v1| 国产av国产精品国产| 九草在线视频观看| 日韩成人av中文字幕在线观看| www.av在线官网国产| 91精品国产国语对白视频| 久久久久国产网址| 久久韩国三级中文字幕| 成人免费观看视频高清| 一二三四中文在线观看免费高清| 在线看a的网站| 亚洲国产色片| 精品人妻偷拍中文字幕| 侵犯人妻中文字幕一二三四区| 黄频高清免费视频| 亚洲欧美一区二区三区国产| 国产极品天堂在线| 99热全是精品| 亚洲精品自拍成人| 国产乱来视频区| 免费高清在线观看视频在线观看| 免费看av在线观看网站| av福利片在线| 一级爰片在线观看| 免费久久久久久久精品成人欧美视频| av视频免费观看在线观看| www.熟女人妻精品国产| 久久久精品免费免费高清| 成人漫画全彩无遮挡| 亚洲精品美女久久久久99蜜臀 | 在线看a的网站| 水蜜桃什么品种好| 一级片免费观看大全| 97人妻天天添夜夜摸| 亚洲av中文av极速乱| 欧美精品亚洲一区二区| 亚洲一码二码三码区别大吗| 老司机影院成人| 国产福利在线免费观看视频| 成年人午夜在线观看视频| 波野结衣二区三区在线| 欧美精品国产亚洲| 在线观看免费高清a一片| 中国三级夫妇交换| 麻豆乱淫一区二区| 美女中出高潮动态图| 最新的欧美精品一区二区| 欧美国产精品va在线观看不卡| 美女高潮到喷水免费观看| 国产片内射在线| 国产日韩欧美在线精品| 91午夜精品亚洲一区二区三区| 国产成人免费无遮挡视频| 午夜免费男女啪啪视频观看| 99久久中文字幕三级久久日本| 999精品在线视频| 男女免费视频国产| 毛片一级片免费看久久久久| 国产午夜精品一二区理论片| 中文欧美无线码| 人妻少妇偷人精品九色| 色94色欧美一区二区| 天堂8中文在线网| 热99久久久久精品小说推荐| 十分钟在线观看高清视频www| 午夜精品国产一区二区电影| av在线观看视频网站免费| 国产片特级美女逼逼视频| 亚洲人成77777在线视频| xxxhd国产人妻xxx| 18+在线观看网站| 一本—道久久a久久精品蜜桃钙片| 国产精品久久久久久久久免| 久久久久久久久久久久大奶| 9色porny在线观看| 美女午夜性视频免费| 大片电影免费在线观看免费| 久热这里只有精品99| 美女xxoo啪啪120秒动态图| 亚洲久久久国产精品| 咕卡用的链子| 亚洲美女黄色视频免费看| 国产一区有黄有色的免费视频| 亚洲精品久久午夜乱码| 精品福利永久在线观看| 亚洲图色成人| 国产一区二区 视频在线| 日本-黄色视频高清免费观看| 妹子高潮喷水视频| av网站在线播放免费| 国产精品.久久久| 久久久久国产精品人妻一区二区| 亚洲美女视频黄频| 久久亚洲国产成人精品v| 最近手机中文字幕大全| av免费观看日本| 精品人妻在线不人妻| 欧美日韩国产mv在线观看视频| 青草久久国产| 午夜影院在线不卡| 亚洲精品久久午夜乱码| 电影成人av| 精品亚洲成a人片在线观看| 色视频在线一区二区三区| 精品福利永久在线观看| 国产乱人偷精品视频| av天堂久久9| 成年美女黄网站色视频大全免费| 如何舔出高潮| 97在线视频观看| 欧美日韩综合久久久久久| 日韩中文字幕欧美一区二区 | 91精品三级在线观看| 亚洲一区中文字幕在线| 免费播放大片免费观看视频在线观看| 国产xxxxx性猛交| 久久午夜综合久久蜜桃| 欧美日本中文国产一区发布| 精品少妇内射三级| 香蕉丝袜av| 在线免费观看不下载黄p国产| 色婷婷av一区二区三区视频| 99国产精品免费福利视频| 97人妻天天添夜夜摸| 十分钟在线观看高清视频www| 久久鲁丝午夜福利片| 水蜜桃什么品种好| 日韩大片免费观看网站| 成年人免费黄色播放视频| 亚洲综合色网址| 精品酒店卫生间| 欧美人与性动交α欧美软件| 日韩精品有码人妻一区|