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

    基于LS-SVM的多系統(tǒng)融合GNSS-MR潮位反演

    2022-04-06 10:42:56游高沖羅孝文尹海博王朝陽
    海洋學(xué)研究 2022年1期
    關(guān)鍵詞:潮位乘法分辨率

    游高沖,郭 杭,羅孝文,尹海博,王朝陽

    (1.南昌大學(xué)信息工程學(xué)院,江西 南昌 330031; 2.自然資源部第二海洋研究所,浙江 杭州 310012)

    0 引言

    監(jiān)測海面高度變化對研究全球氣候變化、構(gòu)建海洋高程基準(zhǔn)、預(yù)防沿海臺風(fēng)等自然災(zāi)害具有重要意義。全球?qū)Ш蕉ㄎ幌到y(tǒng)多徑反射(GNSS Multipath Reflectometry,GNSS-MR)技術(shù)是基于全球定位系統(tǒng)多路經(jīng)效應(yīng)發(fā)展起來的一種新型地基遙感技術(shù),利用地表反射信號與導(dǎo)航衛(wèi)星直射信號的相干信號獲取地表環(huán)境信息[1]。相比于全球定位系統(tǒng)干涉(GNSS Reflectometry,GNSS-R)遙感技術(shù)[2],GNSS-MR不需要另外架設(shè)左旋圓極化天線[3],具有設(shè)備要求低、信號源穩(wěn)定、全天候、覆蓋范圍廣等優(yōu)勢,被廣泛應(yīng)用于地表土壤水分監(jiān)測[4]、積雪監(jiān)測[5]和海洋潮汐監(jiān)測[6-7]等領(lǐng)域中。

    LARSONKM et al[8-11]提出基于載波SNR(Signal-to-Noise Ratio)值的GNSS-MR技術(shù),利用該技術(shù)對土壤水分和雪深進(jìn)行了反演研究,并推廣到潮位反演。由于GNSS-MR僅使用低高度角信噪比數(shù)據(jù),在進(jìn)行潮位監(jiān)測時(shí)常存在時(shí)間分辨率和反演精度低的問題,為此,國內(nèi)外學(xué)者開展了一系列的改進(jìn)研究。針對GNSS-MR中低階多項(xiàng)式提取反射信號存在信號混雜的問題,王瑞芳[12]、ZHANG et al[13]等采用經(jīng)驗(yàn)?zāi)B(tài)分解提取反射信號的方法,提純了海面反射信號,提高了潮位反演精度。為提高GNSS-MR潮位監(jiān)測的時(shí)間分辨率,王杰 等[14]采用小波分析的方法,提取反射信號瞬時(shí)頻率進(jìn)行潮位反演,有效地增加了潮位反演點(diǎn)數(shù),但該方法提高潮位監(jiān)測時(shí)間分辨率有限,仍難以滿足高精度、高時(shí)間分辨率潮位監(jiān)測的要求。陳發(fā)德 等[15]利用GPS、BDS、Galileo多系統(tǒng)數(shù)據(jù)進(jìn)行融合潮位反演,有效提高了潮位反演的連續(xù)性和時(shí)間分辨率,但該方法將不同系統(tǒng)的反演結(jié)果無差別地組合在一起作為反演結(jié)果,聯(lián)合反演潮位的均方根誤差(Root Mean Squared Error,RMSE)為0.323 m,反演精度偏低。何秀鳳 等[16]在進(jìn)行風(fēng)暴潮反演時(shí),采用滑動窗口最小二乘法,對多系統(tǒng)反演數(shù)據(jù)進(jìn)行融合,在提高了風(fēng)暴潮監(jiān)測的時(shí)間分辨率的同時(shí),也提高了監(jiān)測的精度,但反演精度仍為分米級。

    為進(jìn)一步提高GNSS-MR潮位反演的時(shí)間分辨率和精度,本文提出基于LS-SVM的多系統(tǒng)數(shù)據(jù)融合潮位反演方法,利用多個(gè)系統(tǒng)的潮位反演數(shù)據(jù),對LS-SVM進(jìn)行訓(xùn)練,得到多系統(tǒng)數(shù)據(jù)融合的LS-SVM潮位反演模型,進(jìn)而估計(jì)出融合后的潮位值,實(shí)現(xiàn)高時(shí)間分辨率、高精度潮位監(jiān)測。

    1 數(shù)據(jù)與方法

    1.1 GNSS-MR潮位反演方法

    圖1為利用GNSS-MR技術(shù)反演潮位變化的示意圖。只考慮1次反射的情況下,載波多路徑相位延遲φ可以表示成:

    (1)

    式中:λ表示載波波長,h表示天線相位中心到海面的高度,θ表示衛(wèi)星高度角。受多路徑效應(yīng)的影響,海面高度信息以直、反射相干涉的方式保存在混合信號的信噪比中。海面多路徑效應(yīng)下SNR滿足如下等式:

    (2)

    式中:Ad是直射信號幅值,Ar是反射信號幅值,φ是多路徑效應(yīng)相位延遲量。SNR可視為直射信號作用的趨勢項(xiàng)與反射信號作用的周期項(xiàng)的疊加,圖2給出了GPS 01號衛(wèi)星的SNR序列和去趨勢后的殘差序列SNRm。由圖2a可知,SNR整體呈拋物線形式,通過二次擬合去趨勢,即可得到反射信號作用的周期項(xiàng)(圖2b),周期項(xiàng)SNRm滿足如下等式:

    SNRm=Acos(4πhλ-1sinθ)

    (3)

    式中:A是信號振幅。記f=2hλ-1,t=sinθ,則上式可進(jìn)一步整理為如下等式:

    SNRm=Acos(2πft)

    (4)

    周期信號SNRm的頻率包含著與潮位反演相關(guān)的h值,SNRm是以高度角正弦值為變量的不等間距信號,無法通過傅里葉變換獲取。因此,本文采用Lomb-Scargle(L-S)頻譜分析[17],提取SNRm的頻率f,再由h=f·λ/2、hs=H-h,求出海面高度值hs。圖3給出了2020年214年積日部分GPS衛(wèi)星反射信號SNRm及相應(yīng)的L-S頻譜圖,L-S頻譜圖(圖3b)中曲線峰值所對應(yīng)的橫坐標(biāo)為有效反射高h(yuǎn)。

    圖1 GNSS-MR監(jiān)測潮位變化示意圖

    圖2 GPS 01號衛(wèi)星SNR序列(a)和二次擬合去趨勢后的殘差SNRm序列(b)Fig.2 SNR sequence diagram (a) and residual SNRm sequence chart after two fitting detrended trend (b) of GPS 01 satellite

    圖3 HKQT站GPS部分衛(wèi)星海面反射信號SNRm(a)和相應(yīng)的Lomb-Scargle頻譜(b)Fig.3 Partial satellite sea surface reflection signal SNRm (a) and corresponding Lomb-Scargle spectrum (b) of GPS at HKQT station

    1.2 LS-SVM多系統(tǒng)數(shù)據(jù)融合方法

    根據(jù)SVM思想,在高維空間的回歸方程為

    f(X)=〈ω,φ(x)〉+b

    (5)

    式中:〈·,·〉表示點(diǎn)積;ω∈Rnk,是原權(quán)重空間中的權(quán)重向量;φ(·):Rn→Rnk,是將輸入數(shù)據(jù)映射到高維特征空間的非線性函數(shù);b是偏置量,b∈R。

    對于LS-SVM可將優(yōu)化問題轉(zhuǎn)化為[18]

    (6)

    式中:ek為樣本誤差;γ表示正則化參數(shù),在噪聲較小的情況下,較小的γ可以避免過度擬合。此時(shí)的約束條件為

    yk=〈ω,φ(xk)〉+b+ek,k=1,…,N

    (7)

    拉格朗日函數(shù)為

    (8)

    式中:αk為拉格朗日乘子,αk∈R。根據(jù)KKT條件消去ω和ek,得到如下線性方程式:

    (9)

    式中:y=[y1,y2…,yN]T;1v=[1,1,…,1]T;α=[α1,α2,…,αN]T;Ωk,l=〈φ(xk),φ(xl)〉=K(xk,xl),k,l=1,…,N。

    利用最小二乘法對公式進(jìn)行求解,即可得到LS-SVM 非線性回歸模型:

    (10)

    基于LS-SVM多系統(tǒng)數(shù)據(jù)融合方法的潮位反演性能與核函數(shù)K、核參數(shù)σ和正則化參數(shù)γ的選取密切相關(guān)[18],本文選用能很好反映模型復(fù)雜程度和普適性較好的徑向基函數(shù)為核函數(shù),通過網(wǎng)格搜索法優(yōu)選σ和γ參數(shù)?;贚S-SVM的多系統(tǒng)融合模型潮位反演流程如圖4所示。

    圖4 基于LS-SVM的多系統(tǒng)融合潮位反演流程圖Fig.4 Flow chart of multi system fusion tidal level inversion based on LS-SVM

    1.3 實(shí)驗(yàn)數(shù)據(jù)及處理

    為驗(yàn)證基于LS-SVM的多系統(tǒng)融合潮位監(jiān)測方法,本文選取2020年香港北岸HKQT站(22.29103351°N,114.21322098°E)214至244年積日的衛(wèi)星觀測數(shù)據(jù)進(jìn)行潮位反演實(shí)驗(yàn)。該測站配備有TrimbleNETR9型號接收機(jī),可以接收GPS、BDS、Galileo、GLONASS四個(gè)系統(tǒng)的衛(wèi)星信號。此外,該站同一位置建設(shè)有驗(yàn)潮站,提供1次/6 min的潮位數(shù)據(jù),可以用于對LS-SVM多系統(tǒng)融合潮位監(jiān)測方法的精度進(jìn)行評估。

    潮位數(shù)據(jù)由國際海平面實(shí)施網(wǎng)( http://www.ioc-sealevelmonitoring.org/)提供,該潮位數(shù)據(jù)以平均海平面(Mean Sea Level, MSL)為基準(zhǔn)。為使反演潮位值與驗(yàn)潮站的潮位數(shù)據(jù)具有統(tǒng)一的基準(zhǔn),本文依據(jù)該網(wǎng)站提供的驗(yàn)潮站修訂本地參考圖(Revised Local Reference, RLR)進(jìn)行相應(yīng)的高程轉(zhuǎn)換,可得到HKQT站天線相位中心到MSL的高。圖5中RLR的橢球高為-8.274±0.004 m,經(jīng)精密單點(diǎn)定位解算出天線相位中心的橢球高為5.683 m,將天線橢球高減去RLR橢球高,再減去7.01 m即可得到天線相位中心到MSL的高度值,經(jīng)解算高度值為6.407 m。

    潮位反演所需的載波波長、高度角、方向角、天線高等基本參數(shù)設(shè)定如表1所示。

    圖5 驗(yàn)潮站修訂本地參考圖Fig.5 Revising local reference diagram of tide gauge station

    表1 GNSS-MR潮位反演參數(shù)表Tab.1 Inversion parameters for GNSS-MR tide level

    2 結(jié)果與討論

    2.1 單系統(tǒng)潮位反演結(jié)果分析

    利用香港HKQT站提供的2020年214年積日至2020年244年積日連續(xù)30 d的GPS、BDS、Galileo、GLONASS衛(wèi)星數(shù)據(jù)對HKQT站附近潮位進(jìn)行反演,結(jié)果如圖6所示。圖6a、6c、6e和6g分別為GPS、BDS、Galileo和GLONASS反演結(jié)果及驗(yàn)潮站實(shí)測的潮位變化曲線,從中可以看出,HKQT站附近潮位在 ±1 m 區(qū)間范圍內(nèi)周期性波動,各系統(tǒng)反演點(diǎn)均落在曲線上或曲線附近,反演點(diǎn)的變化趨勢同驗(yàn)潮站實(shí)測潮位變化趨勢相同,反演潮位值和實(shí)測值具有明顯的一致性,表明4個(gè)系統(tǒng)均能有效實(shí)現(xiàn)潮位監(jiān)測。將反演潮位值減去同時(shí)刻驗(yàn)潮站實(shí)測潮位值,即可得到反演潮位的偏差值,圖6b、6d、6f和6h分別GPS、BDS、Galileo和GLONASS系統(tǒng)反演潮位的偏差序列,各系統(tǒng)反演偏差均在±0.5 m范圍內(nèi)變化。本文采用均方根誤差、相關(guān)系數(shù)和時(shí)間分辨率作為潮位反演的性能評定指標(biāo),評定結(jié)果如表2所示。

    表2 單系統(tǒng)潮位反演性能對比Tab.2 Comparison of inversion performance of single systems tidal level

    由表2可知,4個(gè)系統(tǒng)均能用于HKQT站附近的潮位監(jiān)測,且各系統(tǒng)反演潮位的精度相當(dāng),反演精度在0.120~0.139 m之間,相關(guān)性系數(shù)在0.953 6~0.963 2 之間,反演潮位值同驗(yàn)潮站實(shí)測值之間具有極強(qiáng)的相關(guān)性。各系統(tǒng)反演潮位的時(shí)間分辨率差異較大,其中GLONASS系統(tǒng)潮位監(jiān)測時(shí)間分辨率最高,平均1.23 h一個(gè)反演值;BDS系統(tǒng)潮位監(jiān)測時(shí)間分辨率最低,平均5.56 h一個(gè)反演值??傮w而言,單系統(tǒng)GNSS-MR潮位反演精度和時(shí)間分辨率均較差,難以滿足高精度、高時(shí)間分辨率潮位反演性能需求。

    2.2 LS-SVM融合潮位反演結(jié)果分析

    為提高潮位監(jiān)測的時(shí)間分辨率和精度,考慮到不同系統(tǒng)之間的互補(bǔ)性,本文采用LS-SVM對GPS、BDS、Galileo、GLONASS反演潮位數(shù)據(jù)進(jìn)行融合。首先,將4個(gè)系統(tǒng)的反演潮位值組合在一起,作為 LS-SVM 的訓(xùn)練集;然后對訓(xùn)練集進(jìn)行LS-SVM訓(xùn)練,并采用十倍交叉驗(yàn)證方法分析確定正則優(yōu)化參數(shù)γ和核參數(shù)σ,得到潮位LS-SVM模型;最后,由LS-SVM潮位反演模型估計(jì)出潮位值。

    圖6 GPS、BDS、Galileo、GLONASS反演潮位結(jié)果、反演偏差及驗(yàn)潮站實(shí)測潮位曲線Fig.6 GPS, BDS, Galileo and GLONASS respectively inversion tidal level results, inversion deviations and measuring tide level curves at tide gauge station

    圖7a是GPS、BDS、Galileo、GLONASS反演結(jié)果組合得到的訓(xùn)練集序列,4個(gè)系統(tǒng)組合后反演點(diǎn)數(shù)目明顯增多,反演點(diǎn)數(shù)目達(dá)到1 487個(gè),平均每0.5 h一個(gè)反演值,相比于GLONASS系統(tǒng),時(shí)間分辨率提高了59.3%以上,極大地增強(qiáng)了GNSS-MR潮位監(jiān)測的時(shí)間分辨率。圖7b是與訓(xùn)練集相對應(yīng)的反演偏差序列,反演偏差在±0.5 m區(qū)間內(nèi)變化,訓(xùn)練集的RMSE值為0.140 m,由此可見,直接將不同系統(tǒng)反演結(jié)果疊加在一起,并不能提高潮位反演精度。圖7c 是訓(xùn)練集經(jīng)LS-SVM訓(xùn)練后得到的LS-SVM潮位反演模型估計(jì)的潮位時(shí)間序列。由圖7d可知,LS-SVM潮位模型估計(jì)值的偏差幅值多集中在-0.2~0.2 m區(qū)間內(nèi),較訓(xùn)練集偏差明顯減小。經(jīng)分析,LS-SVM潮位反演模型估計(jì)潮位值的RMSE值為 0.053 m,較訓(xùn)練集RMSE值減小了62.1%,較單系統(tǒng)潮位反演RMSE最小的BDS系統(tǒng),RMSE值減小了55.8%,較單系統(tǒng)潮位反演時(shí)間分辨率最高的GLONASS系統(tǒng),時(shí)間分辨率提高了59.3%。由此可見,LS-SVM融合模型在提高潮位反演時(shí)間分辨率的同時(shí),也提高了潮位反演精度。

    圖7 4種系統(tǒng)組合的訓(xùn)練集和經(jīng)LS-SVM融合后的潮位反演結(jié)果及相應(yīng)的反演偏差序列Fig.7 The training set of the four system combination and the inversion result of tidal level after LS-SVM fusion and the corresponding inversion deviation sequence

    圖8是LS-SVM融合模型反演潮位值與實(shí)測潮位值的相關(guān)性散點(diǎn)圖,經(jīng)LS-SVM融合后,反演潮位值與實(shí)測潮位值更緊密地分布在直線附近,反演潮位值與實(shí)測值的相關(guān)系數(shù)為0.992 5,較訓(xùn)練集提高了4.2%;較單系統(tǒng)潮位反演相關(guān)性系數(shù)最大的GLONASS系統(tǒng),相關(guān)性系數(shù)提高了3.0%,有效提高了潮位監(jiān)測的精度。

    圖8 LS-SVM多系統(tǒng)融合反演潮位值與驗(yàn)潮站實(shí)測潮位值相關(guān)性圖Fig.8 Correlation diagram of LS-SVM multi system fusion inversion of tidal level and measured tidal level at tide gauge station

    2.3 多種潮位融合算法反演結(jié)果對比分析

    為了分析LS-SVM潮位融合算法的性能,本節(jié)將對比LS-SVM、SVR、滑動窗口最小二乘法3種潮位融合算法。基于SVR的多系統(tǒng)融合方法同LS-SVM方法類似,對多系統(tǒng)組合訓(xùn)練集進(jìn)行SVR訓(xùn)練,得到潮位反演模型,進(jìn)而估計(jì)出潮位值。本文在利用SVR進(jìn)行多系統(tǒng)數(shù)據(jù)融合時(shí),選用高斯徑向核函數(shù)RBF和網(wǎng)格尋優(yōu)法進(jìn)行參數(shù)優(yōu)選,不敏感損失系數(shù)ε設(shè)定為0.1,SVR算法詳見文獻(xiàn)[19]。

    圖9a和9b分別為基于SVR的數(shù)據(jù)融合潮位反演結(jié)果圖和反演偏差圖。經(jīng)SVR融合后,反演潮位的RMSE值為0.111 m,較訓(xùn)練集RMSE值減小了20.7%;較BDS潮位反演RMSE值減小了7.5%,反演精度有所提高。圖9c和9d分別為基于滑動窗口最小二乘法的多系統(tǒng)融合潮位反演結(jié)果和偏差圖,同訓(xùn)練集相比,偏差變化幅值明顯減小。經(jīng)滑動窗口最小二乘法融合后,RMSE值為0.090 m,較訓(xùn)練集RMSE減小了35.7%,較BDS系統(tǒng)RMSE值減小了25%。圖10為基于SVR算法和滑動窗口最小二乘法的多系統(tǒng)融合潮位反演結(jié)果與驗(yàn)潮站實(shí)測值的相關(guān)性圖,經(jīng)分析,SVR反演結(jié)果與驗(yàn)潮站實(shí)測值之間的相關(guān)系數(shù)為0.971 1,較訓(xùn)練集相關(guān)性系數(shù)提高了2.0%;滑動窗口最小二乘法反演結(jié)果與驗(yàn)潮站實(shí)測值的相關(guān)系數(shù)為0.980 4,較訓(xùn)練集提高了2.9%。由此可知,這兩種方法均能提高潮位反演精度。此外,基于滑動窗口最小二乘法的多系統(tǒng)融合潮位反演方法的反演精度為厘米級,優(yōu)于基于SVR的多系統(tǒng)融合潮位反演方法。

    圖9 基于SVR算法和滑動窗口最小二乘法的多系統(tǒng)數(shù)據(jù)融合潮位反演結(jié)果及偏差Fig.9 The inversion result and deviation of multi system data fusion tide level based on SVR algorithm and sliding window least square method

    圖10 基于SVR算法和滑動窗口最小二乘法反演潮位與實(shí)測潮位的相關(guān)性圖Fig.10 Correlation diagram between the inversion tide level based on the SVR algorithm and the sliding window least square method and the measured tide level

    由表3可知,3種多系統(tǒng)融合模型均能提高GNSS-MR反演潮位的精度和時(shí)間分辨率,且時(shí)間分辨率相同,均為0.5 h。LS-SVM模型較SVR模型,RMSE值減小了52.3%,相關(guān)系數(shù)提高2.2%;較滑動窗口最小二乘法,RMSE值減小了41.1%,相關(guān)系數(shù)提高了1.2%。LS-SVM融合模型優(yōu)于SVR融合模型和滑動窗口最小二乘法融合模型。

    表3 3種融合模型潮位反演性能對比Tab.3 Comparison of inversion performance of three kinds of fusion models

    3 結(jié)論

    本文針對單系統(tǒng)GNSS-MR潮位監(jiān)測中時(shí)間分辨率和反演精度低的問題,提出一種基于LS-SVM的多系統(tǒng)融合潮位反演方法,利用香港HKQT站提供的衛(wèi)星觀測數(shù)據(jù),對比了滑動窗口最小二乘法、SVR算法和LS-SVM算法在多系統(tǒng)融合潮位反演中的表現(xiàn),結(jié)果表明:

    (1)單系統(tǒng)GNSS-MR潮位反演的RMSE最小值為0.120 m、最大值為0.139 m;相關(guān)系數(shù)最大值為 0.963 2,最小值為0.953 6;GLONASS潮位反演的時(shí)間分辨率最高,BDS潮位反演的時(shí)間分辨率最低。單系統(tǒng)潮位監(jiān)測的時(shí)間分辨率和反演精度偏差無法滿足高精度、高時(shí)間分辨率潮位監(jiān)測的要求。

    (2)在多系統(tǒng)融合GNSS-MR潮位反演方面,基于LS-SVM的多系統(tǒng)融合潮位反演比滑動窗口最小二乘法、SVR算法具有更優(yōu)的潮位反演性能。相比于SVR模型,LS-SVM模型的RMSE值減小了52.3%,相關(guān)系數(shù)提高了2.2%;相比于滑動窗口最小二乘法,LS-SVM模型的RMSE值減小了41.1%,相關(guān)系數(shù)提高了1.2%。

    本文所用方法較單系統(tǒng)而言,時(shí)間分辨率提升到0.5 h,為進(jìn)一步提高時(shí)間分辨率,在未來可以采用小波變換、NTFT等視頻變換方法,提取瞬時(shí)反射信號的頻率值,達(dá)到提升時(shí)間分辨率的效果。

    猜你喜歡
    潮位乘法分辨率
    算乘法
    基于距離倒數(shù)加權(quán)的多站潮位改正方法可行性分析
    我們一起來學(xué)習(xí)“乘法的初步認(rèn)識”
    唐山市警戒潮位標(biāo)志物維護(hù)研究
    《整式的乘法與因式分解》鞏固練習(xí)
    把加法變成乘法
    EM算法的參數(shù)分辨率
    多潮位站海道地形測量潮位控制方法研究
    原生VS最大那些混淆視聽的“分辨率”概念
    基于改進(jìn)的OLS-RBF模型的感潮河段潮位預(yù)測研究
    欧美精品啪啪一区二区三区| 99久久成人亚洲精品观看| 观看美女的网站| 五月玫瑰六月丁香| 国产高清激情床上av| 精华霜和精华液先用哪个| 国产伦精品一区二区三区视频9| www.色视频.com| 国产av不卡久久| 国产一级毛片七仙女欲春2| 美女高潮的动态| 欧美日韩福利视频一区二区| 午夜免费激情av| 欧美xxxx黑人xx丫x性爽| 欧美激情在线99| 久久久久国内视频| 91在线精品国自产拍蜜月| 国产精品av视频在线免费观看| 午夜免费成人在线视频| 美女免费视频网站| 欧美激情久久久久久爽电影| 老司机深夜福利视频在线观看| 中文字幕免费在线视频6| 欧美一级a爱片免费观看看| 婷婷精品国产亚洲av| 五月伊人婷婷丁香| 日本a在线网址| 亚洲欧美日韩高清专用| 露出奶头的视频| 97人妻精品一区二区三区麻豆| 亚洲色图av天堂| 特级一级黄色大片| 天堂av国产一区二区熟女人妻| 日日夜夜操网爽| 亚洲av二区三区四区| 国产乱人伦免费视频| 俺也久久电影网| 精品熟女少妇八av免费久了| 欧美中文日本在线观看视频| 国产熟女xx| 久久精品人妻少妇| av在线天堂中文字幕| 好男人在线观看高清免费视频| 中文字幕人妻熟人妻熟丝袜美| 一个人看的www免费观看视频| 日韩中字成人| 97热精品久久久久久| 免费黄网站久久成人精品 | 国产真实伦视频高清在线观看 | or卡值多少钱| 日本免费a在线| 桃红色精品国产亚洲av| 久9热在线精品视频| 国产精品综合久久久久久久免费| 男人舔女人下体高潮全视频| 亚洲自偷自拍三级| 亚洲av第一区精品v没综合| 久久九九热精品免费| 91麻豆av在线| 欧美日本视频| 国产午夜福利久久久久久| 在线观看av片永久免费下载| www.色视频.com| 亚洲av熟女| 亚洲人成电影免费在线| 日本三级黄在线观看| 国产亚洲精品av在线| 1000部很黄的大片| 午夜福利在线在线| 国产精品一及| 18禁在线播放成人免费| 一区二区三区免费毛片| 少妇人妻精品综合一区二区 | 又粗又爽又猛毛片免费看| av在线蜜桃| 最新中文字幕久久久久| 成年女人永久免费观看视频| 国产av在哪里看| 成人美女网站在线观看视频| 国产伦精品一区二区三区四那| 欧美日本视频| 一区二区三区高清视频在线| 九九久久精品国产亚洲av麻豆| 最近视频中文字幕2019在线8| 亚洲五月婷婷丁香| 黄色丝袜av网址大全| 成人国产综合亚洲| 深夜a级毛片| 男女床上黄色一级片免费看| 狠狠狠狠99中文字幕| 亚洲人成电影免费在线| 久久久精品大字幕| 国产精品不卡视频一区二区 | 69人妻影院| 国产美女午夜福利| 丰满人妻熟妇乱又伦精品不卡| 国产久久久一区二区三区| 亚洲成人精品中文字幕电影| 精品一区二区三区人妻视频| 亚洲欧美精品综合久久99| 毛片女人毛片| 无人区码免费观看不卡| 久久久国产成人精品二区| 国产大屁股一区二区在线视频| 欧美日韩瑟瑟在线播放| 免费看a级黄色片| 国产成人福利小说| 亚洲五月天丁香| 亚洲久久久久久中文字幕| 国产伦精品一区二区三区四那| 在线观看66精品国产| 久久国产乱子伦精品免费另类| 美女大奶头视频| 少妇被粗大猛烈的视频| 91麻豆精品激情在线观看国产| 久久婷婷人人爽人人干人人爱| 99热只有精品国产| 91麻豆精品激情在线观看国产| 免费观看人在逋| 简卡轻食公司| 欧美高清性xxxxhd video| 欧美日韩乱码在线| 色综合站精品国产| 观看美女的网站| 麻豆国产av国片精品| 中文字幕精品亚洲无线码一区| 蜜桃久久精品国产亚洲av| 男女那种视频在线观看| 欧美激情国产日韩精品一区| 精品不卡国产一区二区三区| 午夜免费成人在线视频| 十八禁人妻一区二区| 非洲黑人性xxxx精品又粗又长| 欧美日韩亚洲国产一区二区在线观看| 一本精品99久久精品77| 欧美+日韩+精品| 日本黄色视频三级网站网址| 精品久久久久久久久亚洲 | 每晚都被弄得嗷嗷叫到高潮| 亚洲国产精品sss在线观看| 精品一区二区免费观看| 变态另类丝袜制服| 国产高清有码在线观看视频| 欧美日韩福利视频一区二区| 91麻豆av在线| 怎么达到女性高潮| 看十八女毛片水多多多| 精品人妻1区二区| 97超视频在线观看视频| 麻豆国产97在线/欧美| 琪琪午夜伦伦电影理论片6080| 性色av乱码一区二区三区2| 丝袜美腿在线中文| 99热6这里只有精品| 蜜桃久久精品国产亚洲av| 丁香六月欧美| 久久香蕉精品热| 亚洲av成人不卡在线观看播放网| 亚洲av美国av| 亚洲av日韩精品久久久久久密| 欧美在线一区亚洲| 亚洲欧美激情综合另类| 赤兔流量卡办理| 亚洲最大成人中文| 欧美三级亚洲精品| 日韩成人在线观看一区二区三区| 高清日韩中文字幕在线| 国产精品爽爽va在线观看网站| 久久亚洲真实| 女人被狂操c到高潮| 国产av麻豆久久久久久久| 国产精品98久久久久久宅男小说| 国产三级中文精品| 婷婷六月久久综合丁香| 1024手机看黄色片| 国产毛片a区久久久久| 国产私拍福利视频在线观看| 欧美黄色淫秽网站| 国产精品一区二区性色av| 露出奶头的视频| 亚洲美女搞黄在线观看 | 欧美日韩瑟瑟在线播放| 色综合站精品国产| 日本三级黄在线观看| 日本免费一区二区三区高清不卡| 男女之事视频高清在线观看| 日日夜夜操网爽| 精华霜和精华液先用哪个| 青草久久国产| 白带黄色成豆腐渣| 久久精品影院6| 欧美丝袜亚洲另类 | av国产免费在线观看| 国模一区二区三区四区视频| 18+在线观看网站| 最新中文字幕久久久久| 99久久精品国产亚洲精品| h日本视频在线播放| 一二三四社区在线视频社区8| 性欧美人与动物交配| 一级黄色大片毛片| 长腿黑丝高跟| 欧美高清性xxxxhd video| ponron亚洲| 精品免费久久久久久久清纯| a级一级毛片免费在线观看| 18禁黄网站禁片免费观看直播| 欧美成狂野欧美在线观看| 好看av亚洲va欧美ⅴa在| 色吧在线观看| 久久久久九九精品影院| 夜夜夜夜夜久久久久| 久久久久久久久久黄片| 黄片小视频在线播放| 欧美成人a在线观看| 99在线视频只有这里精品首页| 国产久久久一区二区三区| 国产亚洲精品av在线| 成人毛片a级毛片在线播放| 久久精品夜夜夜夜夜久久蜜豆| 丰满乱子伦码专区| 中国美女看黄片| 欧美成人免费av一区二区三区| 99国产精品一区二区三区| ponron亚洲| www.999成人在线观看| 69人妻影院| 在线免费观看不下载黄p国产 | 久久久久久久久中文| 男女做爰动态图高潮gif福利片| 久99久视频精品免费| 亚洲在线观看片| 国产爱豆传媒在线观看| 性插视频无遮挡在线免费观看| 欧美潮喷喷水| 色吧在线观看| 久久99热6这里只有精品| 小蜜桃在线观看免费完整版高清| 香蕉av资源在线| 毛片女人毛片| 小蜜桃在线观看免费完整版高清| 丰满人妻一区二区三区视频av| 国产69精品久久久久777片| 亚洲中文字幕一区二区三区有码在线看| 日本免费a在线| 亚洲人成网站高清观看| 精品国产三级普通话版| 男人狂女人下面高潮的视频| 午夜福利在线在线| 久久久国产成人免费| 国产黄片美女视频| 日日干狠狠操夜夜爽| 欧美最黄视频在线播放免费| 精品乱码久久久久久99久播| 欧美国产日韩亚洲一区| 美女免费视频网站| 麻豆av噜噜一区二区三区| 国产一区二区三区视频了| 亚洲aⅴ乱码一区二区在线播放| 99久久精品一区二区三区| 欧美在线黄色| 国产单亲对白刺激| 人人妻人人澡欧美一区二区| 日韩 亚洲 欧美在线| 三级男女做爰猛烈吃奶摸视频| 亚洲性夜色夜夜综合| 淫妇啪啪啪对白视频| eeuss影院久久| 国产成人av教育| 全区人妻精品视频| 毛片一级片免费看久久久久 | 国产乱人伦免费视频| 国产精品精品国产色婷婷| 国产成人a区在线观看| 欧美在线黄色| 国产成人欧美在线观看| 国产男靠女视频免费网站| 高清毛片免费观看视频网站| 脱女人内裤的视频| 不卡一级毛片| 免费电影在线观看免费观看| 18禁黄网站禁片免费观看直播| 波野结衣二区三区在线| 日韩成人在线观看一区二区三区| 观看美女的网站| aaaaa片日本免费| 丁香欧美五月| 久久亚洲真实| 丰满人妻熟妇乱又伦精品不卡| 亚洲专区国产一区二区| 亚洲,欧美精品.| 精品一区二区三区视频在线| 小蜜桃在线观看免费完整版高清| 亚洲国产精品sss在线观看| 夜夜爽天天搞| 老司机午夜福利在线观看视频| 精品一区二区三区av网在线观看| 久久香蕉精品热| 国产成人福利小说| 国内精品久久久久久久电影| 熟妇人妻久久中文字幕3abv| 国产高清激情床上av| 国产精华一区二区三区| 性欧美人与动物交配| 9191精品国产免费久久| 日韩欧美精品v在线| 国内毛片毛片毛片毛片毛片| 无人区码免费观看不卡| 看免费av毛片| 欧美又色又爽又黄视频| or卡值多少钱| 首页视频小说图片口味搜索| 国产精品99久久久久久久久| 亚洲精品在线观看二区| 极品教师在线免费播放| 老司机福利观看| 久久精品国产亚洲av涩爱 | 久久99热这里只有精品18| 日韩av在线大香蕉| 黄色丝袜av网址大全| 精品人妻偷拍中文字幕| 欧美成人一区二区免费高清观看| 午夜视频国产福利| 3wmmmm亚洲av在线观看| 中文在线观看免费www的网站| 亚洲av五月六月丁香网| 一二三四社区在线视频社区8| 国产精品伦人一区二区| 麻豆成人av在线观看| 免费电影在线观看免费观看| 99热这里只有是精品在线观看 | 国产主播在线观看一区二区| 亚洲黑人精品在线| a级毛片a级免费在线| 波多野结衣高清作品| 亚洲欧美日韩卡通动漫| 首页视频小说图片口味搜索| 伊人久久精品亚洲午夜| 91麻豆精品激情在线观看国产| 国产成人影院久久av| 亚洲av成人精品一区久久| 日本黄大片高清| 深爱激情五月婷婷| 中文字幕久久专区| 国产真实伦视频高清在线观看 | 免费大片18禁| 最近最新中文字幕大全电影3| 亚洲一区二区三区色噜噜| 欧美bdsm另类| 日韩亚洲欧美综合| 美女被艹到高潮喷水动态| 亚洲精品一区av在线观看| 亚洲精华国产精华精| 日日夜夜操网爽| 欧美又色又爽又黄视频| 欧美zozozo另类| 中文在线观看免费www的网站| 熟女电影av网| 久久中文看片网| 午夜激情欧美在线| 好看av亚洲va欧美ⅴa在| 精品国产三级普通话版| 真人一进一出gif抽搐免费| 麻豆一二三区av精品| 97碰自拍视频| 俄罗斯特黄特色一大片| 久久精品国产99精品国产亚洲性色| 51国产日韩欧美| 蜜桃久久精品国产亚洲av| av在线蜜桃| a级毛片a级免费在线| 色综合欧美亚洲国产小说| 亚洲一区二区三区色噜噜| 免费观看人在逋| 人妻丰满熟妇av一区二区三区| 亚洲国产精品999在线| 男插女下体视频免费在线播放| 国产成+人综合+亚洲专区| 精品久久久久久久久av| 国产精品野战在线观看| 久久久精品欧美日韩精品| 在线播放无遮挡| 欧美性感艳星| 伊人久久精品亚洲午夜| 乱人视频在线观看| 久久精品国产亚洲av涩爱 | 少妇高潮的动态图| 成人无遮挡网站| 午夜两性在线视频| 国产欧美日韩一区二区精品| 国产爱豆传媒在线观看| 国产日本99.免费观看| 国产三级在线视频| 午夜福利成人在线免费观看| 丰满的人妻完整版| 国产成人a区在线观看| www.999成人在线观看| 99热精品在线国产| 成人特级av手机在线观看| 亚洲av中文字字幕乱码综合| 一个人免费在线观看的高清视频| 熟妇人妻久久中文字幕3abv| 国产精品久久久久久人妻精品电影| 国产乱人伦免费视频| 一边摸一边抽搐一进一小说| 999久久久精品免费观看国产| 亚洲精华国产精华精| 国产精品亚洲一级av第二区| 亚洲欧美日韩高清专用| 日韩欧美国产一区二区入口| 久久草成人影院| eeuss影院久久| 日韩欧美在线二视频| 婷婷六月久久综合丁香| 国产成人a区在线观看| 变态另类丝袜制服| 日日摸夜夜添夜夜添av毛片 | 又爽又黄无遮挡网站| av国产免费在线观看| 99热这里只有精品一区| 一进一出好大好爽视频| 欧美一级a爱片免费观看看| 亚洲精品色激情综合| 精品国内亚洲2022精品成人| 我要看日韩黄色一级片| 欧美激情久久久久久爽电影| 十八禁人妻一区二区| 人妻夜夜爽99麻豆av| 狂野欧美白嫩少妇大欣赏| 国产三级中文精品| 国产精品乱码一区二三区的特点| 超碰av人人做人人爽久久| 熟妇人妻久久中文字幕3abv| 嫩草影院入口| 亚洲精品在线美女| 三级毛片av免费| 亚洲成人精品中文字幕电影| 欧美高清性xxxxhd video| 国产精品一及| a级一级毛片免费在线观看| 黄色视频,在线免费观看| 亚洲欧美日韩东京热| 一个人观看的视频www高清免费观看| 精品不卡国产一区二区三区| 99精品久久久久人妻精品| 亚洲av成人不卡在线观看播放网| 国产精品亚洲av一区麻豆| avwww免费| 中文字幕人妻熟人妻熟丝袜美| 国产成年人精品一区二区| 天堂av国产一区二区熟女人妻| 美女xxoo啪啪120秒动态图 | 别揉我奶头~嗯~啊~动态视频| 中文字幕av成人在线电影| 欧美丝袜亚洲另类 | 色播亚洲综合网| 亚洲中文字幕日韩| 亚洲成人久久爱视频| 亚洲乱码一区二区免费版| 午夜久久久久精精品| 国产精品av视频在线免费观看| 精品久久久久久久久久久久久| 久久久久久久久大av| 真实男女啪啪啪动态图| 国产精品日韩av在线免费观看| 亚洲无线观看免费| 久久天躁狠狠躁夜夜2o2o| 久久精品91蜜桃| 熟女人妻精品中文字幕| 免费看a级黄色片| 99热这里只有精品一区| 亚洲熟妇中文字幕五十中出| 好看av亚洲va欧美ⅴa在| 久久亚洲真实| 国产高潮美女av| 国产日本99.免费观看| 国产不卡一卡二| 国产高清三级在线| 乱人视频在线观看| av福利片在线观看| 日韩 亚洲 欧美在线| 成人三级黄色视频| 直男gayav资源| 九色成人免费人妻av| 中文资源天堂在线| 中国美女看黄片| 成人毛片a级毛片在线播放| 国产精品亚洲av一区麻豆| 日韩欧美一区二区三区在线观看| 国产精品野战在线观看| 两个人视频免费观看高清| 日本黄色视频三级网站网址| 亚洲欧美激情综合另类| 自拍偷自拍亚洲精品老妇| 亚洲国产精品久久男人天堂| 亚洲精品粉嫩美女一区| 午夜视频国产福利| 女生性感内裤真人,穿戴方法视频| 99久久精品一区二区三区| 精品久久国产蜜桃| 丁香六月欧美| 高清日韩中文字幕在线| 日韩欧美在线二视频| 男女做爰动态图高潮gif福利片| 欧美高清性xxxxhd video| 久久精品久久久久久噜噜老黄 | 毛片女人毛片| 免费无遮挡裸体视频| 搡老妇女老女人老熟妇| 成人精品一区二区免费| 亚洲最大成人av| 日日夜夜操网爽| 国产成人福利小说| 精品久久久久久久久久免费视频| 久久久成人免费电影| 51国产日韩欧美| 成人国产一区最新在线观看| 中文字幕精品亚洲无线码一区| 91午夜精品亚洲一区二区三区 | 高清在线国产一区| 精品午夜福利在线看| 日韩大尺度精品在线看网址| 亚州av有码| 高潮久久久久久久久久久不卡| 天堂动漫精品| 成人一区二区视频在线观看| www.999成人在线观看| 中文字幕久久专区| 床上黄色一级片| 啦啦啦韩国在线观看视频| 国产国拍精品亚洲av在线观看| 精品久久久久久久久久免费视频| 亚洲自偷自拍三级| 国产探花在线观看一区二区| 久久精品综合一区二区三区| 欧美潮喷喷水| 国产亚洲欧美在线一区二区| 国产精品日韩av在线免费观看| 中文在线观看免费www的网站| 久久精品国产亚洲av香蕉五月| 国产又黄又爽又无遮挡在线| 亚洲av.av天堂| 怎么达到女性高潮| 欧美日韩中文字幕国产精品一区二区三区| 欧美最黄视频在线播放免费| 成年女人永久免费观看视频| 美女被艹到高潮喷水动态| 国产一区二区激情短视频| 久久久精品大字幕| 国产精品,欧美在线| 亚洲人成网站高清观看| 一二三四社区在线视频社区8| 欧美潮喷喷水| 露出奶头的视频| 搞女人的毛片| 亚洲aⅴ乱码一区二区在线播放| 网址你懂的国产日韩在线| 午夜福利在线观看吧| 全区人妻精品视频| 好男人在线观看高清免费视频| 国产乱人视频| 国产日本99.免费观看| 国产精品一区二区性色av| 欧美黑人巨大hd| 熟女电影av网| 午夜精品久久久久久毛片777| 老熟妇仑乱视频hdxx| 午夜激情欧美在线| 亚洲中文字幕一区二区三区有码在线看| 丝袜美腿在线中文| 一个人看视频在线观看www免费| 免费观看人在逋| 在线观看舔阴道视频| 国产高清有码在线观看视频| 在线观看免费视频日本深夜| 99久久精品国产亚洲精品| 欧美黄色片欧美黄色片| 757午夜福利合集在线观看| 能在线免费观看的黄片| 男女做爰动态图高潮gif福利片| 欧美高清成人免费视频www| 午夜福利在线观看吧| 国产男靠女视频免费网站| 亚洲av成人精品一区久久| 岛国在线免费视频观看| 国产乱人视频| 国产综合懂色| 91在线精品国自产拍蜜月| 久久热精品热| 丰满人妻熟妇乱又伦精品不卡| 欧美性猛交╳xxx乱大交人| 永久网站在线| 蜜桃久久精品国产亚洲av| 亚洲av电影在线进入| 日本黄色视频三级网站网址| 亚洲欧美精品综合久久99| 国产一区二区在线av高清观看| 网址你懂的国产日韩在线| 97超级碰碰碰精品色视频在线观看| 亚州av有码| 伦理电影大哥的女人| 亚洲欧美激情综合另类| 国产高清视频在线观看网站| 欧美不卡视频在线免费观看| 免费在线观看亚洲国产| 在线免费观看的www视频| 伊人久久精品亚洲午夜| 精品人妻熟女av久视频| 亚洲av.av天堂| 伊人久久精品亚洲午夜| 一进一出抽搐动态| 我的女老师完整版在线观看| 亚洲欧美日韩东京热| 亚洲人成电影免费在线| 国产私拍福利视频在线观看|