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

    EEMD-Holt-Winters的電離層TEC預(yù)報(bào)模型

    2021-05-13 03:56:58周強(qiáng)波
    測(cè)繪工程 2021年3期
    關(guān)鍵詞:低緯度電離層殘差

    周強(qiáng)波

    (核工業(yè)二三〇研究所, 湖南 長(zhǎng)沙 410007)

    電離層是距地面60~2 000 km的電離區(qū)域,其包含的大量自由電子使穿越其中的無(wú)線信號(hào)發(fā)生折射,是GNSS測(cè)量、導(dǎo)航、定位的主要誤差源之一。電離層總電子含量(Total Electric Content, TEC)是表征電離層結(jié)構(gòu)的一個(gè)重要參數(shù),精確TEC模型對(duì)衛(wèi)星導(dǎo)航定位、無(wú)線電傳播通訊、自然現(xiàn)象的認(rèn)知反演具有重要意義。電離層模型可分為理論模型和經(jīng)驗(yàn)?zāi)P蛢深?,理論模型可以定性地描述電離層的結(jié)構(gòu)特征和變化規(guī)律,通常較為復(fù)雜且精度較低。經(jīng)驗(yàn)?zāi)P褪歉鶕?jù)長(zhǎng)期的觀測(cè)數(shù)據(jù)利用一系列數(shù)學(xué)物理方法擬合建立,使用這類模型在進(jìn)行某一地區(qū)的一定時(shí)間段內(nèi)電離層TEC預(yù)報(bào)可獲得較高的精度。

    目前常見的電離層TEC預(yù)報(bào)模型有灰色系統(tǒng)模型[1-2]、時(shí)間序列模型[3-4]、神經(jīng)網(wǎng)絡(luò)模型[5-6]等。其中,灰色系統(tǒng)模型在小數(shù)據(jù)量預(yù)報(bào)中具有優(yōu)勢(shì),在面對(duì)大量TEC數(shù)據(jù)時(shí)預(yù)報(bào)精度較差。時(shí)間序列模型在短期預(yù)報(bào)領(lǐng)域精度較高,隨著預(yù)報(bào)時(shí)間的增加精度隨之下降。神經(jīng)網(wǎng)絡(luò)模型需要大量的數(shù)據(jù)進(jìn)行訓(xùn)練才能擬合出精確的模型,且參數(shù)設(shè)置復(fù)雜,預(yù)報(bào)結(jié)果不穩(wěn)定。針對(duì)TEC序列非線性、不平穩(wěn)、復(fù)雜無(wú)序的特點(diǎn),部分學(xué)者基于先分解再預(yù)測(cè)的思路建立預(yù)報(bào)模型[7-8],在TEC預(yù)報(bào)領(lǐng)域取得了較好的效果。文獻(xiàn)[9-10]針對(duì)時(shí)空?qǐng)龅臅r(shí)間序列預(yù)報(bào),利用經(jīng)驗(yàn)正交函數(shù)(Empirical Orthogonal Function, EOF)和方差貢獻(xiàn)率提取主時(shí)間成分剔除冗余信息,更有效利用時(shí)間信息取得了更好的效果。文獻(xiàn)[11-12]利用經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition, EMD)首先做數(shù)據(jù)預(yù)處理,然后利用預(yù)報(bào)模型對(duì)各分量進(jìn)行預(yù)測(cè)重構(gòu),但EMD分解會(huì)出現(xiàn)模態(tài)混疊的現(xiàn)象,進(jìn)而影響預(yù)報(bào)精度。

    為進(jìn)一步提高預(yù)報(bào)精度,本文利用EMD的改進(jìn)方法集合經(jīng)驗(yàn)?zāi)B(tài)分解(Ensemble Empirical Mode Decomposition, EEMD)與Holt-Winters模型結(jié)合對(duì)電離層TEC序列進(jìn)行預(yù)報(bào)研究,建立EEMD-Holt-Winters模型。

    1 模型基本原理

    1.1 EEMD方法

    EEMD算法是時(shí)頻領(lǐng)域常用的分解方法,在分解過(guò)程中無(wú)需設(shè)置參數(shù),具有極強(qiáng)的自適應(yīng)性,其主要作用是將非線性、不平穩(wěn)信號(hào)分解為不同頻率的特征模態(tài)分量(Intrinsic Mode Function, IMF),針對(duì)EMD分解結(jié)果中不同時(shí)間尺度出現(xiàn)在同一IMF分量中的模態(tài)混疊問(wèn)題,提出添加輔助白噪聲的方法進(jìn)行分析。通過(guò)在原始信號(hào)中多次添加高斯白噪聲,使信號(hào)中不同特征的時(shí)間尺度以白噪聲為參考自動(dòng)投射到合適的尺度上,利用白噪聲均值為0的特性能夠有效抑制EMD算法在分解時(shí)產(chǎn)生的模態(tài)混疊問(wèn)題[13]。EEMD算法的具體步驟如下:

    1)在原始信號(hào)中加入一定幅值的高斯白噪聲得到新的待分解信號(hào),表示為:

    xn(t)=x(t)+n(t).

    (1)

    式中:x(t)為原始信號(hào);n(t)為高斯白噪聲;xn(t)為加噪后的待處理信號(hào);

    2)使用EMD算法將xn(t)分解為若干特征模態(tài)分量imfi(t),i=1,2,…,n和一個(gè)殘余項(xiàng)rn(t);

    3)每次加入均方根值相等的不同高斯白噪聲重復(fù)步驟1)、2),得到k組不同的imf分量和殘余項(xiàng);

    4)將分解得到的k組imf分量和殘余項(xiàng)取平均值,即可得到最終的分解結(jié)果。

    (2)

    1.2 Holt-Winters模型

    Holt-Winters模型是一種基于數(shù)理統(tǒng)計(jì)原理的時(shí)間序列預(yù)報(bào)模型,該模型的基本思想是將數(shù)據(jù)分解研究,與指數(shù)平滑法結(jié)合對(duì)時(shí)間序列的趨勢(shì)、季節(jié)性、波動(dòng)變化進(jìn)行估計(jì),適用于預(yù)測(cè)具有明顯趨勢(shì)和季節(jié)性的時(shí)間序列,對(duì)周日重復(fù)的數(shù)據(jù)預(yù)報(bào)效果較好。其主要包括無(wú)季節(jié)模型、Holt-Winters加法模型和Holt-Winters乘法模型3種[14]。在電離層TEC序列預(yù)報(bào)領(lǐng)域,文獻(xiàn)[15]驗(yàn)證了Holt-Winters加法模型對(duì)此類序列的預(yù)報(bào)效果最優(yōu)。其基本公式為:

    St=α(Xt-It-L)+(1-α)(St-1-bt-1),

    (3)

    It=β(Xt-St)+(1-β)It-L,

    (4)

    bt=γ(St-St-1)+(1-γ)bt-1,

    (5)

    Ft+m=St+mbt+It-L+m.

    (6)

    式中:Xt為觀測(cè)值;St,It,bt分別為穩(wěn)定成分、季節(jié)成分、趨勢(shì)成分;L為季節(jié)長(zhǎng)度;m為預(yù)測(cè)期數(shù);Ft+m為第m期預(yù)測(cè)值;α,β,γ為平滑參數(shù)。

    本文建立EEMD-Holt-Winters模型首先利用EEMD算法將電離層TEC序列分解得到多個(gè)imf分量和殘余項(xiàng),考慮電離層TEC序列特點(diǎn)選用Holt-Winters加法模型對(duì)imf分量和殘余項(xiàng)進(jìn)行預(yù)報(bào),最后將預(yù)測(cè)得到的各個(gè)分量重構(gòu)得到最終預(yù)測(cè)值。

    2 實(shí)驗(yàn)分析

    2.1 數(shù)據(jù)準(zhǔn)備

    本文利用IGS組織提供的2018年全年數(shù)據(jù),選取年積日152~171和249~268兩個(gè)時(shí)段低緯度(120°E5°N)和中緯度(120°E35°N)時(shí)間分辨率1 h的TEC數(shù)據(jù)進(jìn)行分析。分別利用前15 d數(shù)據(jù)作為訓(xùn)練樣本,使用本文所建模型預(yù)報(bào)后5 d的TEC變化趨勢(shì)。首先對(duì)選取的數(shù)據(jù)進(jìn)行EEMD分解,然后利用Holt-Winters加法模型預(yù)測(cè)各分量數(shù)據(jù),最后將所有分量預(yù)測(cè)值相加得到后5 d預(yù)報(bào)值,將預(yù)報(bào)結(jié)果與IGS實(shí)際值比對(duì),選用均方根誤差(RMSE)、平均殘差(Δ)和平均相對(duì)精度(P)來(lái)進(jìn)行精度評(píng)定。具體公式為:

    (7)

    (8)

    (9)

    地磁活動(dòng)和太陽(yáng)黑子是影響電離層TEC值的主要因素,一般使用Kp指數(shù)表示地磁活動(dòng)的強(qiáng)弱,當(dāng)Kp≤3時(shí)認(rèn)為無(wú)磁暴現(xiàn)象,電離層不受地磁活動(dòng)影響。本文選取兩個(gè)時(shí)段的Kp指數(shù)和F10.7指數(shù),如圖1所示,年積日152~171時(shí)段僅有1 d的Kp指數(shù)超過(guò)4,17 d的Kp指數(shù)在2以下,平均Kp指數(shù)為1.33,表明此時(shí)段地磁活動(dòng)平靜。年積日249~268時(shí)段Kp指數(shù)在年積日254和265超過(guò)4,有弱磁暴,其余Kp指數(shù)也大多在2以下,平均值為1.93,此時(shí)段地磁活動(dòng)較活躍。2018年是太陽(yáng)活動(dòng)低峰年,無(wú)黑子日達(dá)224 d,從圖1可以看出F10.7指數(shù)基本在70左右,對(duì)電離層TEC影響較小。

    圖1 Kp指數(shù)和F10.7指數(shù)

    2.2 實(shí)驗(yàn)結(jié)果與分析

    利用本文建立模型對(duì)所選數(shù)據(jù)進(jìn)行預(yù)報(bào)分析,選取低緯度地區(qū)平靜日TEC序列進(jìn)行分解,低緯度(120°E5°N)地區(qū)平靜日TEC序列經(jīng)EMD和EEMD分解后的結(jié)果如圖2和圖3所示。TEC序列經(jīng)EMD分解后得到由高頻到低頻排列的6個(gè)imf分量和1個(gè)殘余項(xiàng),可以看出模態(tài)混疊現(xiàn)象嚴(yán)重,多個(gè)分量雜糅了不同特征時(shí)間尺度的信息。經(jīng)EEMD分解后得到7個(gè)imf分量和1個(gè)殘余項(xiàng),分解更精細(xì),從第3和第4個(gè)分量可以看出EEMD與EMD相比有效抑制了模態(tài)混疊現(xiàn)象。

    圖2 低緯度地區(qū)TEC序列EMD分解結(jié)果

    圖3 低緯度地區(qū)TEC序列EEMD分解結(jié)果

    地磁平靜時(shí)段低緯度和中緯度地區(qū)兩種模型5 d預(yù)報(bào)結(jié)果和實(shí)際值的對(duì)比和殘差如圖4和圖5所示,可以看出TEC序列在第3日發(fā)生擾動(dòng),低緯度地區(qū)TEC值接近30 TECu,中緯度地區(qū)達(dá)23 TECu,較其余4 d高出8 TECu。單一Holt-Winters模型和EEMD-Holt-Winters模型均能較好地預(yù)報(bào)出TEC序列的變化趨勢(shì),TEC最小值預(yù)報(bào)效果較為準(zhǔn)確,與實(shí)際變化趨勢(shì)相吻合。從殘差分布圖可以看出,EEMD-Holt-Winters模型的預(yù)報(bào)殘差與單一Holt-Winters模型分布基本相同但數(shù)值更小,證明與單一模型相比組合模型的預(yù)報(bào)結(jié)果與實(shí)際值更接近,具有更好的預(yù)報(bào)精度。

    地磁平靜時(shí)段兩種模型預(yù)報(bào)精度對(duì)比見表1??梢钥闯鲈诘卮牌届o時(shí)段低緯度地區(qū)EEMD-Holt-Winters模型和單一Holt-Winters的RMSE分別為1.67 TECu和1.90 TECu,平均殘差分別為1.27 TECu和1.51 TECu,平均相對(duì)精度分別為85.16%和80.90%。在中緯度地區(qū),EEMD-Holt-Winters模型和單一Holt-Winters的RMSE分別為1.54 TECu和1.62 TECu,平均殘差分別為0.99 TECu和1.08 TECu,平均相對(duì)精度分別為91.71%和90.65%。從具體數(shù)值來(lái)看,中緯度地區(qū)的預(yù)報(bào)效果優(yōu)于低緯度地區(qū),組合模型的預(yù)報(bào)誤差更小,相對(duì)精度更優(yōu)。

    圖4 地磁平靜時(shí)段兩種模型預(yù)報(bào)結(jié)果示意圖

    圖5 地磁平靜時(shí)段兩種模型預(yù)報(bào)殘差示意圖

    表1 地磁平靜時(shí)段兩種模型預(yù)報(bào)精度對(duì)比

    地磁平靜時(shí)段低緯度地區(qū)5 d預(yù)報(bào)結(jié)果的殘差分布見表2,從中可以看出EEMD-Holt-Winters模型有53.33%的預(yù)報(bào)殘差在1 TECu以內(nèi),30%的預(yù)報(bào)殘差在1~2 TECu之間,16.67%的預(yù)報(bào)殘差大于3 TECu。單一Holt-Winters模型預(yù)報(bào)殘差在1 TECu以內(nèi)的僅有36.67%,在1~2 TECu之間的預(yù)報(bào)殘差占比40.83%,殘差大于3 TECu的預(yù)報(bào)值占比22.50%。整體來(lái)看EEMD-Holt-Winters組合模型的預(yù)報(bào)殘差較單一模型更優(yōu),基本保持在2 TECu以內(nèi),僅有少部分超過(guò)2 TECu,而單一模型的預(yù)報(bào)誤差更大,年積日170有43.48%的預(yù)報(bào)殘差大于2 TECu,高誤差的預(yù)報(bào)結(jié)果更多。

    為進(jìn)一步評(píng)定本文模型預(yù)報(bào)效果,利用本文模型預(yù)測(cè)地磁活躍時(shí)段的TEC變化趨勢(shì)。年積日264~268時(shí)段兩種模型的預(yù)報(bào)結(jié)果如圖6所示,可以看出此時(shí)段的TEC序列波動(dòng)劇烈,低緯度地區(qū)TEC最大值基本穩(wěn)定在30 TECu左右,最高達(dá)33 TECu,中緯度地區(qū)波動(dòng)明顯,第2日伴隨弱磁暴TEC最大值高達(dá)22 TECu,可見磁暴對(duì)電離層TEC有一定影響,其余基本保持在13 TECu左右。兩種模型在地磁活躍時(shí)段預(yù)報(bào)殘差的分布結(jié)果如圖7所示,可以看出在地磁活躍時(shí)段兩種模型的預(yù)報(bào)殘差分布大致相同,具體數(shù)值較地磁平靜期更大,EEMD-Holt-Winters組合模型的預(yù)報(bào)殘差略優(yōu)于單一Holt-Winters模型。

    表2 地磁平靜時(shí)段低緯度地區(qū)兩種模型殘差統(tǒng)計(jì)

    圖6 地磁活躍時(shí)段兩種模型預(yù)報(bào)結(jié)果示意圖

    圖7 地磁活躍時(shí)段兩種模型預(yù)報(bào)殘差示意圖

    不同模型在地磁活躍時(shí)段5 d預(yù)報(bào)結(jié)果的精度評(píng)定指標(biāo)對(duì)比見表3,從中可以看出在年積日265擾動(dòng)期間預(yù)報(bào)誤差明顯,中緯度地區(qū)單一Holt-Winters模型的預(yù)報(bào)殘差高達(dá)4.37 TECu,低緯度地區(qū)達(dá)2.59 TECu。從整體上看單日的預(yù)報(bào)結(jié)果兩種模型大致相同,組合模型的預(yù)報(bào)效果較單一模型好,5 d的平均RMSE組合模型低、中緯度分別為2.32 TECu和2.18 TECu,平均相對(duì)精度分別為83.33%和86.83%,單一Holt-Winters模型5 d的平均RMSE組合模型低、中緯度分別為2.40 TECu和2.24 TECu,平均相對(duì)精度分別為82.23%和85.55%,預(yù)報(bào)效果與EEMD-Holt-Winters組合模型相比略差。

    地磁活躍時(shí)段低、高緯度地區(qū)兩種模型的殘差統(tǒng)計(jì)見表4,從中可以看出低緯度地區(qū)單一Holt-Winters模型|Δ|<1.0 TECu的預(yù)報(bào)值占比為25%,EEMD-Holt-Winters組合模型為28.83%,單一模型|Δ|≥3.0 TECu的預(yù)報(bào)值占比為21.67%,組合模型為18.33%,兩種模型1.0≤|Δ|<2.0 TECu和2.0≤|Δ|<3.0 TECu的預(yù)報(bào)值占比相同,均為30.83%和22.50%,這表明在低緯度地區(qū)組合模型高誤差預(yù)報(bào)值少于單一模型,大多數(shù)預(yù)報(bào)誤差在3 TECu以內(nèi)。在中緯度地區(qū),兩種模型|Δ|<1.0 TECu的預(yù)報(bào)值占比分別為47.50%和55.00%,1.0≤|Δ|<2.0 TECu的預(yù)報(bào)值占比分別為24.17%和20.83%,2.0≤|Δ|<3.0 TECu的預(yù)報(bào)值占比分別為15.00%和11.67%,|Δ|≥3.0 TECu的預(yù)報(bào)值占比分別為13.33%和12.50%,進(jìn)一步驗(yàn)證了EEMD-Holt-Winters組合模型誤差較小的預(yù)報(bào)值數(shù)量?jī)?yōu)于單一Holt-Winters模型,預(yù)報(bào)值較大的預(yù)報(bào)值數(shù)量少于單一Holt-Winters模型。

    表3 不同模型年積日264~268精度評(píng)定對(duì)比

    表4 地磁活躍時(shí)段兩種模型殘差統(tǒng)計(jì)

    3 結(jié)束語(yǔ)

    本文利用IGS提供的電離層TEC數(shù)據(jù),分別利用單一Holt-Winters模型和EEMD-Holt-Winters組合模型對(duì)地磁平靜時(shí)段和地磁活躍時(shí)段的TEC序列進(jìn)行建模預(yù)報(bào)。實(shí)驗(yàn)結(jié)果表明,在地磁活動(dòng)平靜時(shí)段,組合模型預(yù)報(bào)效果顯著優(yōu)于單一模型,在地磁活動(dòng)活躍時(shí)段,組合模型與單一模型預(yù)報(bào)結(jié)果略優(yōu)于單一模型,驗(yàn)證了EEMD算法與Holt-Winters組合模型的可行性和優(yōu)勢(shì)。然而,組合模型在地磁活躍期的預(yù)報(bào)效果穩(wěn)定性較差,精度并無(wú)明顯提高,仍需進(jìn)行下一步優(yōu)化改進(jìn)。

    猜你喜歡
    低緯度電離層殘差
    中國(guó)低緯度地區(qū)的夏季太陽(yáng)總輻照度研究
    基于雙向GRU與殘差擬合的車輛跟馳建模
    一種電離層TEC格點(diǎn)預(yù)測(cè)模型
    Kalman濾波估算電離層延遲的一種優(yōu)化方法
    基于殘差學(xué)習(xí)的自適應(yīng)無(wú)人機(jī)目標(biāo)跟蹤算法
    基于FY-3C/MWHTS觀測(cè)資料反演中低緯度海面氣壓
    高海拔低緯度條件下的8個(gè)家蠶品種比較試驗(yàn)
    基于遞歸殘差網(wǎng)絡(luò)的圖像超分辨率重建
    電離層對(duì)中高軌SAR影響機(jī)理研究
    Linux Shell語(yǔ)言在電離層解算中的應(yīng)用
    国产精品麻豆人妻色哟哟久久| 久久久国产精品麻豆| 国产有黄有色有爽视频| 成人国产麻豆网| 久久影院123| 午夜福利影视在线免费观看| 亚洲成人一二三区av| 日韩欧美一区视频在线观看| 最近最新中文字幕大全免费视频 | 中文字幕人妻丝袜一区二区 | 亚洲欧洲日产国产| 老司机影院成人| 丁香六月天网| 欧美日韩一区二区视频在线观看视频在线| 国产97色在线日韩免费| 国产成人a∨麻豆精品| 超碰97精品在线观看| 久久久欧美国产精品| 成人毛片a级毛片在线播放| 99久久中文字幕三级久久日本| 国产免费视频播放在线视频| 啦啦啦啦在线视频资源| 久久久精品免费免费高清| 国产福利在线免费观看视频| 国产男女内射视频| 日韩 亚洲 欧美在线| 国产免费又黄又爽又色| 男女下面插进去视频免费观看| 成人国语在线视频| 日本欧美国产在线视频| 久久国产精品男人的天堂亚洲| 人体艺术视频欧美日本| 亚洲国产av影院在线观看| 一二三四中文在线观看免费高清| 精品亚洲乱码少妇综合久久| 中文字幕另类日韩欧美亚洲嫩草| 自拍欧美九色日韩亚洲蝌蚪91| 日韩欧美精品免费久久| 亚洲av在线观看美女高潮| 国产一区二区激情短视频 | 国产日韩欧美亚洲二区| 欧美最新免费一区二区三区| 中国三级夫妇交换| 咕卡用的链子| 免费少妇av软件| 精品国产一区二区三区久久久樱花| 少妇人妻 视频| 免费观看无遮挡的男女| 成年av动漫网址| 美女高潮到喷水免费观看| 日本午夜av视频| 欧美日韩视频高清一区二区三区二| 国产精品国产三级国产专区5o| 日韩在线高清观看一区二区三区| 亚洲熟女精品中文字幕| 少妇猛男粗大的猛烈进出视频| 在线亚洲精品国产二区图片欧美| 90打野战视频偷拍视频| 色吧在线观看| 亚洲精品日本国产第一区| 巨乳人妻的诱惑在线观看| 免费观看在线日韩| 亚洲内射少妇av| 黄色视频在线播放观看不卡| 熟女av电影| 这个男人来自地球电影免费观看 | 亚洲欧美一区二区三区国产| 午夜福利,免费看| 中文欧美无线码| 国产一区二区 视频在线| 丰满饥渴人妻一区二区三| 成人午夜精彩视频在线观看| 国产黄频视频在线观看| 免费观看性生交大片5| 国产精品久久久久久精品古装| 欧美人与善性xxx| 国产精品成人在线| 久久久久久久久免费视频了| 晚上一个人看的免费电影| 欧美日本中文国产一区发布| 日韩一卡2卡3卡4卡2021年| 精品卡一卡二卡四卡免费| 少妇人妻久久综合中文| 黄片播放在线免费| 亚洲国产欧美日韩在线播放| 美女脱内裤让男人舔精品视频| 女人精品久久久久毛片| 在线亚洲精品国产二区图片欧美| 国产欧美亚洲国产| 日韩制服骚丝袜av| 亚洲四区av| 青春草视频在线免费观看| 午夜福利在线观看免费完整高清在| 成年女人毛片免费观看观看9 | 极品人妻少妇av视频| 久久午夜福利片| av在线app专区| 欧美日韩一区二区视频在线观看视频在线| 久久久国产一区二区| 欧美精品亚洲一区二区| 亚洲视频免费观看视频| 欧美bdsm另类| 亚洲欧美一区二区三区久久| h视频一区二区三区| 亚洲av国产av综合av卡| 中文字幕av电影在线播放| 欧美成人午夜免费资源| 精品福利永久在线观看| 亚洲情色 制服丝袜| 国产成人欧美| 中文字幕另类日韩欧美亚洲嫩草| 久久久久人妻精品一区果冻| 精品99又大又爽又粗少妇毛片| 老司机影院毛片| av免费在线看不卡| av女优亚洲男人天堂| 国产不卡av网站在线观看| 如日韩欧美国产精品一区二区三区| 久久久久久久精品精品| 精品人妻一区二区三区麻豆| 亚洲av日韩在线播放| 亚洲欧洲日产国产| 久久久久久久国产电影| 久久精品国产亚洲av高清一级| 国产激情久久老熟女| 欧美成人午夜精品| 成人亚洲精品一区在线观看| 国产精品三级大全| 黄网站色视频无遮挡免费观看| 国产一区二区三区综合在线观看| 人人妻人人澡人人看| 欧美人与性动交α欧美软件| 宅男免费午夜| av免费在线看不卡| 性色av一级| 色94色欧美一区二区| 精品久久蜜臀av无| 黄网站色视频无遮挡免费观看| 久久久久国产网址| av又黄又爽大尺度在线免费看| 一二三四在线观看免费中文在| 精品人妻熟女毛片av久久网站| 日韩一卡2卡3卡4卡2021年| www日本在线高清视频| 少妇的逼水好多| 如日韩欧美国产精品一区二区三区| 男男h啪啪无遮挡| 国产亚洲av片在线观看秒播厂| 国产成人精品无人区| 国产人伦9x9x在线观看 | 亚洲av福利一区| 最近的中文字幕免费完整| 伦精品一区二区三区| 99热国产这里只有精品6| 日韩免费高清中文字幕av| 久久久精品免费免费高清| 午夜av观看不卡| 久久 成人 亚洲| 免费黄频网站在线观看国产| 亚洲美女黄色视频免费看| 99久久综合免费| 精品国产一区二区三区四区第35| 日本午夜av视频| 又黄又粗又硬又大视频| 午夜福利,免费看| 精品国产一区二区久久| 一级a爱视频在线免费观看| 在线精品无人区一区二区三| 久久精品国产自在天天线| 亚洲欧美清纯卡通| 国产又色又爽无遮挡免| 亚洲一级一片aⅴ在线观看| 亚洲五月色婷婷综合| 日本黄色日本黄色录像| 免费久久久久久久精品成人欧美视频| 亚洲精品中文字幕在线视频| 成人毛片60女人毛片免费| 观看美女的网站| 极品少妇高潮喷水抽搐| 精品人妻偷拍中文字幕| 好男人视频免费观看在线| 成年女人在线观看亚洲视频| 精品国产一区二区久久| 免费黄色在线免费观看| 99精国产麻豆久久婷婷| 人人妻人人添人人爽欧美一区卜| 一区福利在线观看| 国产 一区精品| 免费在线观看视频国产中文字幕亚洲 | 久久av网站| 亚洲美女视频黄频| 国产精品国产三级国产专区5o| 夫妻性生交免费视频一级片| 有码 亚洲区| 晚上一个人看的免费电影| 亚洲第一区二区三区不卡| 成人影院久久| 亚洲三级黄色毛片| 日本wwww免费看| av女优亚洲男人天堂| 成年人免费黄色播放视频| 日韩三级伦理在线观看| 免费看av在线观看网站| 永久免费av网站大全| 日日摸夜夜添夜夜爱| 久久婷婷青草| 欧美bdsm另类| 久久精品久久久久久久性| 中文字幕另类日韩欧美亚洲嫩草| 天美传媒精品一区二区| h视频一区二区三区| 视频区图区小说| 久久久久精品性色| 色网站视频免费| av天堂久久9| 新久久久久国产一级毛片| 男女啪啪激烈高潮av片| 亚洲国产最新在线播放| 十八禁网站网址无遮挡| 日韩中文字幕视频在线看片| 免费黄网站久久成人精品| 久久久亚洲精品成人影院| 精品少妇内射三级| 18禁国产床啪视频网站| 欧美精品国产亚洲| 国产深夜福利视频在线观看| 亚洲美女黄色视频免费看| 老司机亚洲免费影院| 99国产精品免费福利视频| 少妇的丰满在线观看| 婷婷成人精品国产| 麻豆av在线久日| 亚洲婷婷狠狠爱综合网| 99re6热这里在线精品视频| 中文字幕色久视频| 18禁动态无遮挡网站| a级毛片在线看网站| 啦啦啦视频在线资源免费观看| 精品一区二区三卡| 最黄视频免费看| 97人妻天天添夜夜摸| 国产免费现黄频在线看| 国产综合精华液| 美女中出高潮动态图| 五月伊人婷婷丁香| 在线观看www视频免费| 狠狠婷婷综合久久久久久88av| 观看美女的网站| 久久狼人影院| 欧美精品一区二区免费开放| 最近的中文字幕免费完整| 99国产综合亚洲精品| 黄片播放在线免费| kizo精华| 又粗又硬又长又爽又黄的视频| 午夜精品国产一区二区电影| 午夜av观看不卡| 欧美日韩精品网址| 国产免费福利视频在线观看| 国产精品久久久久久av不卡| 人人妻人人爽人人添夜夜欢视频| 色哟哟·www| 久久人人97超碰香蕉20202| 国产精品国产三级专区第一集| 免费在线观看视频国产中文字幕亚洲 | www.精华液| av线在线观看网站| 只有这里有精品99| 亚洲精品日本国产第一区| 亚洲av电影在线进入| 亚洲av电影在线观看一区二区三区| 天天影视国产精品| av免费在线看不卡| 高清在线视频一区二区三区| 国产国语露脸激情在线看| 久久久久国产精品人妻一区二区| 精品福利永久在线观看| www日本在线高清视频| 亚洲av电影在线进入| 日韩成人av中文字幕在线观看| 久热这里只有精品99| 国产精品亚洲av一区麻豆 | 久久女婷五月综合色啪小说| 免费观看性生交大片5| 另类亚洲欧美激情| 18禁裸乳无遮挡动漫免费视频| 纯流量卡能插随身wifi吗| freevideosex欧美| 欧美精品亚洲一区二区| 成人毛片a级毛片在线播放| 国语对白做爰xxxⅹ性视频网站| 欧美老熟妇乱子伦牲交| 狠狠婷婷综合久久久久久88av| 久久人人97超碰香蕉20202| 99国产精品免费福利视频| 欧美日韩一级在线毛片| 国产成人a∨麻豆精品| 丝袜喷水一区| 纵有疾风起免费观看全集完整版| 国产综合精华液| 久久久a久久爽久久v久久| 老鸭窝网址在线观看| 精品一区二区免费观看| 超碰成人久久| 在线天堂中文资源库| 亚洲精品乱久久久久久| 成人18禁高潮啪啪吃奶动态图| 久久国内精品自在自线图片| 日韩制服丝袜自拍偷拍| 久久久久精品人妻al黑| 亚洲av成人精品一二三区| 最近手机中文字幕大全| 亚洲精品久久成人aⅴ小说| 一区在线观看完整版| 咕卡用的链子| 国产av一区二区精品久久| 国产熟女午夜一区二区三区| 春色校园在线视频观看| 1024香蕉在线观看| 国产欧美日韩综合在线一区二区| 99久久精品国产国产毛片| 国产不卡av网站在线观看| 久久这里只有精品19| 曰老女人黄片| 欧美黄色片欧美黄色片| 日韩不卡一区二区三区视频在线| 另类精品久久| 国产熟女午夜一区二区三区| 女的被弄到高潮叫床怎么办| 久久精品久久久久久噜噜老黄| 婷婷色av中文字幕| 日本午夜av视频| 久久久久久免费高清国产稀缺| freevideosex欧美| 日韩中文字幕欧美一区二区 | 久久亚洲国产成人精品v| 国产激情久久老熟女| 涩涩av久久男人的天堂| 亚洲精品自拍成人| 99久久综合免费| 久久狼人影院| 一级,二级,三级黄色视频| 有码 亚洲区| 热re99久久精品国产66热6| www.av在线官网国产| 国产精品人妻久久久影院| 超碰97精品在线观看| 欧美日韩一区二区视频在线观看视频在线| 成人漫画全彩无遮挡| 精品人妻一区二区三区麻豆| 2022亚洲国产成人精品| 91精品国产国语对白视频| 日韩电影二区| 一区在线观看完整版| 热99国产精品久久久久久7| 在线精品无人区一区二区三| 精品久久蜜臀av无| 蜜桃国产av成人99| 最近最新中文字幕大全免费视频 | 青春草亚洲视频在线观看| 18+在线观看网站| 亚洲欧洲精品一区二区精品久久久 | 成年av动漫网址| 精品国产乱码久久久久久小说| 欧美日韩亚洲高清精品| 99热国产这里只有精品6| 国产欧美日韩一区二区三区在线| 久久精品国产自在天天线| 纯流量卡能插随身wifi吗| 中文字幕人妻熟女乱码| 99久久人妻综合| 菩萨蛮人人尽说江南好唐韦庄| 老熟女久久久| 亚洲精品乱久久久久久| 国产精品香港三级国产av潘金莲 | 最近最新中文字幕免费大全7| 丁香六月天网| 久久久国产一区二区| 成年动漫av网址| 自拍欧美九色日韩亚洲蝌蚪91| 日韩视频在线欧美| 色哟哟·www| 亚洲人成网站在线观看播放| 国产精品女同一区二区软件| 久久久久国产网址| 午夜91福利影院| 在线天堂最新版资源| 欧美 亚洲 国产 日韩一| 肉色欧美久久久久久久蜜桃| 国产片特级美女逼逼视频| 成人毛片a级毛片在线播放| 欧美精品亚洲一区二区| 久久精品亚洲av国产电影网| 国产成人欧美| 男人操女人黄网站| 亚洲国产成人一精品久久久| 日韩制服骚丝袜av| 一级爰片在线观看| 欧美国产精品一级二级三级| 日韩在线高清观看一区二区三区| 久久青草综合色| 在线天堂中文资源库| 国产免费视频播放在线视频| 国产熟女欧美一区二区| 久久久久精品性色| 亚洲精品,欧美精品| 韩国av在线不卡| 国产精品无大码| 如何舔出高潮| 色婷婷av一区二区三区视频| 久久国产精品男人的天堂亚洲| 免费久久久久久久精品成人欧美视频| 国产精品不卡视频一区二区| 日日爽夜夜爽网站| 在线天堂最新版资源| 亚洲精品日本国产第一区| 国产精品 国内视频| 国产日韩欧美视频二区| 哪个播放器可以免费观看大片| 亚洲国产成人一精品久久久| 一级毛片 在线播放| 丝袜喷水一区| 熟女少妇亚洲综合色aaa.| 另类亚洲欧美激情| 久久青草综合色| 国产精品一区二区在线观看99| 日韩免费高清中文字幕av| 成人毛片60女人毛片免费| 侵犯人妻中文字幕一二三四区| 国产精品国产三级专区第一集| 999精品在线视频| 日韩不卡一区二区三区视频在线| 一级爰片在线观看| 国产人伦9x9x在线观看 | 亚洲欧美清纯卡通| 男女啪啪激烈高潮av片| 香蕉精品网在线| 视频区图区小说| 丝袜喷水一区| 亚洲欧洲精品一区二区精品久久久 | 美女中出高潮动态图| 亚洲精品国产色婷婷电影| 欧美+日韩+精品| 美女福利国产在线| 欧美xxⅹ黑人| 国产乱来视频区| 看非洲黑人一级黄片| 国产在线一区二区三区精| 成人国产av品久久久| 久久久国产精品麻豆| 又黄又粗又硬又大视频| 亚洲精品视频女| 久久婷婷青草| 亚洲精品自拍成人| 性色avwww在线观看| 亚洲,欧美精品.| 可以免费在线观看a视频的电影网站 | 丰满迷人的少妇在线观看| 久久精品aⅴ一区二区三区四区 | 中文乱码字字幕精品一区二区三区| 99精国产麻豆久久婷婷| 大片免费播放器 马上看| 老汉色av国产亚洲站长工具| 黑人欧美特级aaaaaa片| 免费看不卡的av| 午夜久久久在线观看| 国产又爽黄色视频| 色吧在线观看| 国产乱人偷精品视频| 又大又黄又爽视频免费| 欧美 日韩 精品 国产| 两性夫妻黄色片| 天堂俺去俺来也www色官网| freevideosex欧美| 777久久人妻少妇嫩草av网站| 亚洲精品aⅴ在线观看| 建设人人有责人人尽责人人享有的| 免费在线观看黄色视频的| 久久久久久久亚洲中文字幕| 国产一区二区 视频在线| 久久精品人人爽人人爽视色| 下体分泌物呈黄色| 日本午夜av视频| 天天躁狠狠躁夜夜躁狠狠躁| 男女国产视频网站| 久久久久久久久久久免费av| 久久久久精品性色| 欧美日韩精品成人综合77777| 波多野结衣av一区二区av| 亚洲精品国产色婷婷电影| 大香蕉久久网| 在线 av 中文字幕| 亚洲色图 男人天堂 中文字幕| 国产在线一区二区三区精| 欧美中文综合在线视频| 捣出白浆h1v1| 老司机影院成人| 91aial.com中文字幕在线观看| 女人高潮潮喷娇喘18禁视频| 欧美少妇被猛烈插入视频| 不卡视频在线观看欧美| 国产精品熟女久久久久浪| 国产黄频视频在线观看| 在线观看人妻少妇| 两个人看的免费小视频| 久久久久久久久久人人人人人人| 97精品久久久久久久久久精品| 最新中文字幕久久久久| 男女无遮挡免费网站观看| 最近中文字幕高清免费大全6| 狠狠婷婷综合久久久久久88av| 少妇熟女欧美另类| 日韩制服骚丝袜av| 18+在线观看网站| 午夜日韩欧美国产| 午夜日本视频在线| 十八禁高潮呻吟视频| 一级毛片电影观看| 亚洲人成77777在线视频| 亚洲国产色片| 观看av在线不卡| 日韩视频在线欧美| 亚洲色图综合在线观看| 国产成人精品无人区| 欧美 日韩 精品 国产| 国产男人的电影天堂91| 高清视频免费观看一区二区| 性色avwww在线观看| 欧美亚洲 丝袜 人妻 在线| 97在线人人人人妻| 亚洲av电影在线观看一区二区三区| 成人亚洲欧美一区二区av| 九色亚洲精品在线播放| 一区在线观看完整版| 日韩视频在线欧美| 黄色毛片三级朝国网站| 日韩人妻精品一区2区三区| 亚洲美女搞黄在线观看| 伦理电影免费视频| 99国产综合亚洲精品| 天堂俺去俺来也www色官网| 日韩大片免费观看网站| 亚洲久久久国产精品| 春色校园在线视频观看| 日本av免费视频播放| 天天影视国产精品| 亚洲,一卡二卡三卡| 亚洲av欧美aⅴ国产| 伊人久久大香线蕉亚洲五| 色哟哟·www| 90打野战视频偷拍视频| 99热全是精品| 18禁裸乳无遮挡动漫免费视频| 蜜桃在线观看..| 美国免费a级毛片| 91aial.com中文字幕在线观看| 日本欧美视频一区| 啦啦啦在线免费观看视频4| 母亲3免费完整高清在线观看 | 国产 一区精品| 国产免费又黄又爽又色| 黄色毛片三级朝国网站| 麻豆av在线久日| 国产亚洲精品第一综合不卡| 色视频在线一区二区三区| 亚洲一码二码三码区别大吗| 免费黄网站久久成人精品| 国产人伦9x9x在线观看 | 青春草国产在线视频| 成人国语在线视频| 国产在线视频一区二区| 下体分泌物呈黄色| 极品人妻少妇av视频| 午夜福利,免费看| 老司机亚洲免费影院| 婷婷色综合大香蕉| 一区福利在线观看| 好男人视频免费观看在线| 大香蕉久久成人网| 欧美精品高潮呻吟av久久| 日日爽夜夜爽网站| freevideosex欧美| 国产精品久久久av美女十八| 亚洲精品av麻豆狂野| 午夜日本视频在线| 少妇精品久久久久久久| 久久99热这里只频精品6学生| 永久免费av网站大全| 久久精品国产a三级三级三级| 黑人欧美特级aaaaaa片| 纵有疾风起免费观看全集完整版| 老熟女久久久| 国产成人av激情在线播放| 人成视频在线观看免费观看| 精品亚洲乱码少妇综合久久| 色视频在线一区二区三区| 精品福利永久在线观看| 精品人妻一区二区三区麻豆| 校园人妻丝袜中文字幕| 久久久亚洲精品成人影院| 久久久a久久爽久久v久久| 国精品久久久久久国模美| 这个男人来自地球电影免费观看 | 亚洲av欧美aⅴ国产| 中文字幕亚洲精品专区| 在线观看三级黄色| xxxhd国产人妻xxx| 久久久久久久久久人人人人人人| a级毛片黄视频| 人妻 亚洲 视频| 1024视频免费在线观看| 国产精品免费视频内射| 欧美人与性动交α欧美精品济南到 | 少妇猛男粗大的猛烈进出视频| 丝袜美腿诱惑在线|