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

    自動(dòng)站數(shù)據(jù)質(zhì)控及其在實(shí)時(shí)業(yè)務(wù)中的應(yīng)用*

    2021-03-09 00:27:34許駿寧楊志勇田功平葉棟水
    海峽科學(xué) 2021年1期
    關(guān)鍵詞:游程落點(diǎn)觀測(cè)

    袁 偉 許駿寧 楊志勇 田功平 葉棟水

    (1.福建省氣象信息中心,福建 福州 350001;2.福建省減災(zāi)中心,福建 福州 350001;

    3.海峽氣象開放實(shí)驗(yàn)室,福建 廈門 361000)

    1 概述

    區(qū)域自動(dòng)站觀測(cè)資料往往因站點(diǎn)建設(shè)位置、數(shù)量多不易精細(xì)化管理等原因,數(shù)據(jù)質(zhì)量一定程度上難以保證,同時(shí)也鑒于目前技術(shù)手段的投入與限制等問題,一旦數(shù)據(jù)進(jìn)入相關(guān)業(yè)務(wù)系統(tǒng),往往不易察覺,并帶來系列不良影響。本文圍繞區(qū)域自動(dòng)站觀測(cè)數(shù)據(jù),研究實(shí)時(shí)業(yè)務(wù)中的質(zhì)量控制問題,期望在關(guān)鍵時(shí)候能夠阻止顯性錯(cuò)誤,業(yè)務(wù)應(yīng)用前對(duì)疑似錯(cuò)誤提示告警,必要時(shí)對(duì)大概率異常值先行剔除。

    此外,在結(jié)合業(yè)務(wù)需求的研究過程中,我們將隨機(jī)取一次包含降雨的天氣過程的區(qū)域自動(dòng)站、雷達(dá)資料,用于對(duì)文中算法的學(xué)習(xí)與驗(yàn)證,其中區(qū)域自動(dòng)站資料5分鐘觀測(cè)頻率,資料窗口包含天氣過程頭尾共7天的數(shù)據(jù)量進(jìn)行演算。

    2 技術(shù)內(nèi)容

    2.1 誤差的來源分類

    通常來說,觀測(cè)值與真值之間的誤差有系統(tǒng)誤差、測(cè)量誤差、環(huán)境誤差等[1],表達(dá)如下:

    X(t)=X′(t)+E(t)

    其中,X′(t)為真值,X(t)為觀測(cè)值,E(t)為觀測(cè)誤差,t為觀測(cè)時(shí)次。

    立足區(qū)域自動(dòng)站的觀測(cè)數(shù)據(jù),我們使用誤差來表述觀測(cè)值正常與否,觀測(cè)誤差的基本組成部分理解如下:系統(tǒng)誤差(記為:S(t))可理解為設(shè)備故障、接觸不良、通信中斷等情況代入的數(shù)據(jù)誤差;測(cè)量誤差(記為:R(t))可理解為測(cè)量設(shè)備因可能的老化、校正等情況引起的測(cè)量變動(dòng);環(huán)境誤差(記為:V(t))可理解為測(cè)試、人為性質(zhì)代入的數(shù)據(jù)誤差。為簡化分析,我們僅研究E(t)=S(t)+R(t)+V(t)誤差模型,忽略其他誤差源。對(duì)于不同的誤差來源,需采用相應(yīng)的一些技術(shù)手段,本文將通過重構(gòu)或改造或利用一些分析方法,提出誤差抑制的相關(guān)性分析。值得關(guān)注的是,在較長序列的數(shù)據(jù)中,上述誤差可能會(huì)相互轉(zhuǎn)換,因此需要考慮不同計(jì)算法的混合應(yīng)用。

    根據(jù)誤差模型,應(yīng)用歷史資料研究誤差分布,進(jìn)而用于實(shí)時(shí)觀測(cè)的預(yù)測(cè),也是一種研究思路。

    此時(shí)根據(jù)原始觀測(cè)值與經(jīng)過人工訂正后的歷史資料,做相應(yīng)研究,如差分、樣條函數(shù)、最小二乘擬合等相關(guān)性研究,挖掘誤差分布情況。

    2.2 數(shù)據(jù)質(zhì)量控制技術(shù)

    對(duì)于氣象觀測(cè)數(shù)據(jù)的質(zhì)量控制,結(jié)合觀測(cè)設(shè)備的類型和觀測(cè)對(duì)象,國內(nèi)外做了大量的算法研究[2-4],如氣候?qū)W界限值檢查、氣候極值檢查、空間一致性檢驗(yàn)、時(shí)間一致性檢驗(yàn)、統(tǒng)計(jì)分析等。本文的自動(dòng)站數(shù)據(jù)質(zhì)量控制,綜合考慮實(shí)時(shí)業(yè)務(wù)中的時(shí)效性、準(zhǔn)確性,除了結(jié)合這些成熟的技術(shù)外,還將根據(jù)具體的實(shí)時(shí)應(yīng)用,改造或重構(gòu)計(jì)算法,來滿足現(xiàn)有業(yè)務(wù)的部分急需應(yīng)用。

    此外,考慮到業(yè)務(wù)實(shí)時(shí)應(yīng)用,將采用更多的輕量級(jí)的質(zhì)控技術(shù),遴選若干種后,最后做綜合評(píng)分,如果質(zhì)控評(píng)估分值超大,則認(rèn)為數(shù)據(jù)顯性異常,如果評(píng)估分值較大,則認(rèn)為數(shù)據(jù)存疑。其他情況對(duì)數(shù)據(jù)放行。

    根據(jù)觀測(cè)數(shù)據(jù)的統(tǒng)計(jì)分析,本文復(fù)合應(yīng)用IQR、Z-Score、DBSCAN、LSM等方法進(jìn)行平滑、濾波、外推研究,并根據(jù)其計(jì)算結(jié)果作為質(zhì)量控制的依據(jù)。

    2.3 差分游程技術(shù)

    相鄰時(shí)次的觀測(cè)值之差,可以有效削弱系統(tǒng)誤差所帶來的影響,尤其是觀測(cè)對(duì)象本身變動(dòng)有規(guī)律的情況下,差分結(jié)果跳變較大時(shí),比較容易發(fā)現(xiàn)異常值,如氣溫、氣壓等要素。與此同時(shí),當(dāng)測(cè)量值出現(xiàn)不變或變動(dòng)范圍很小時(shí),需要改造差分算法,引入長序列穩(wěn)態(tài)判斷,即需要對(duì)長序列數(shù)據(jù)進(jìn)行監(jiān)控,我們通過實(shí)時(shí)計(jì)算差分游程,來探測(cè)數(shù)據(jù)誤差的影響。

    設(shè)觀測(cè)時(shí)次t1,t2,t3,...,tn,n為當(dāng)前觀測(cè),觀測(cè)值差分:DX(ti)=X(ti)-X(ti-1) ,則差分游程計(jì)算式如下:

    其中α是游程標(biāo)尺閾值,當(dāng)α取0時(shí),為嚴(yán)格的游程過程,容易描述長時(shí)序數(shù)據(jù)的穩(wěn)定性,當(dāng)取α非零值,不同觀測(cè)要素下,適當(dāng)取值,容易描述長時(shí)序數(shù)據(jù)的微變性。統(tǒng)計(jì)連續(xù)0態(tài)或連續(xù)1態(tài)的資料觀測(cè)時(shí)次的次數(shù),稱為穩(wěn)定周期。表1為實(shí)際數(shù)據(jù)在不同標(biāo)尺下的游程穩(wěn)定性估測(cè)。

    表1 游程穩(wěn)定周期評(píng)估

    從表1可見,隨著α增大,穩(wěn)定周期由小逐漸變大,并開始趨于相對(duì)固定。在實(shí)際業(yè)務(wù)中,我們選取穩(wěn)定周期相對(duì)固定的最小α值作為差分游程技術(shù)的游程標(biāo)尺閾值。

    2.4 鄰近判識(shí)

    通過對(duì)區(qū)域自動(dòng)站進(jìn)行地理分析,從地理相似條件、距離、海拔等位置出發(fā),嘗試匹配出每個(gè)站的若干鄰近站,用于對(duì)當(dāng)前站的觀測(cè)數(shù)據(jù)質(zhì)量進(jìn)行研判。與此同時(shí),本文嘗試引入雷達(dá)資料進(jìn)行格點(diǎn)化,用于對(duì)某區(qū)域的自動(dòng)站進(jìn)行降水量輔助質(zhì)量控制判斷。綜合這些方式方法,構(gòu)建區(qū)域站鄰近判識(shí)技術(shù)。為找到最佳鄰近站,首先從地理相似條件出發(fā)進(jìn)行篩選:

    d={D[Fk(B,L,H),F(B,L,H)]}

    其中F(B,L,H)為本站的GCCS地理坐標(biāo),F(xiàn)k(B,L,H)為除本站外的其他站GCCS地理坐標(biāo),D為判斷函數(shù),最簡單的一種方式是取歐氏距離,較復(fù)雜的方式可以考慮等緯度線、等溫線、等壓線等條件。d為篩選集合,在實(shí)際使用時(shí),可取TOP5(例如,按歐氏距離由小到大排序后,取前5個(gè))。通過計(jì)算,則每個(gè)站均能得到相應(yīng)的鄰近站。

    其次,針對(duì)每個(gè)站的鄰近站集,再通過歷史數(shù)據(jù)進(jìn)行對(duì)比,找出其觀測(cè)要素相關(guān)性更好的站。要素相關(guān)性采用均方差計(jì)算(考慮到天氣過程的尺度問題,不同的觀測(cè)要素有不同的相關(guān)性屬性,需要設(shè)計(jì)更復(fù)雜的計(jì)算法,本文假設(shè)在同一尺度下進(jìn)行各要素的相關(guān)性研究),取均方差最小的作為觀測(cè)一致性符合程度的判識(shí)依據(jù):

    代入實(shí)際區(qū)域站相關(guān)信息運(yùn)算后,并不能保證每個(gè)站都有理想的鄰近站。同時(shí)在實(shí)際運(yùn)算過程中,該算法過程可以設(shè)置為動(dòng)態(tài)學(xué)習(xí)過程,確保找到的鄰近站更能反映本站的一些特性,從而為本站數(shù)據(jù)質(zhì)控提供參考依據(jù)。

    天氣過程的尺度選擇依據(jù),我們從某次天氣過程中的雷達(dá)回波演變過程,來粗略劃定一個(gè)區(qū)域,得到一個(gè)回波過程區(qū)域:經(jīng)度范圍[118.00,118.45],緯度范圍[26.75,27.25](按3小時(shí)演變,統(tǒng)計(jì)24小時(shí)的回波覆蓋區(qū)域)。

    在這個(gè)區(qū)域內(nèi),再考慮測(cè)站的地理特性,并對(duì)選定的觀測(cè)站點(diǎn)進(jìn)行遍歷計(jì)算,得到每一個(gè)參考站所屬的鄰近集。以F9742觀測(cè)站為例,得到相應(yīng)的鄰近集如表2所示。

    表2 TOP5鄰近集

    針對(duì)每個(gè)參考站及其鄰近集,分析觀測(cè)要素的相關(guān)性(這里采用要素的均方差):

    進(jìn)一步鎖定一個(gè)時(shí)間段內(nèi)的最佳鄰近集,以此確定每個(gè)觀測(cè)站在實(shí)時(shí)業(yè)務(wù)中相當(dāng)可靠的鄰近集。表3是部分鄰近判識(shí)的測(cè)算過程和最佳鄰近集選擇參考,統(tǒng)計(jì)表明歐式距離遠(yuǎn)近能夠反映出部分相關(guān)性,距離越近相關(guān)性要強(qiáng)些,以F9742為例,對(duì)應(yīng)最佳鄰近集的相關(guān)性計(jì)算結(jié)果如表3所示。

    表3 最佳鄰近集

    根據(jù)每個(gè)參考站的最佳鄰近集,可以采用大數(shù)判決、最小方差等方法,進(jìn)行觀測(cè)預(yù)測(cè),預(yù)測(cè)值可作為參考站的當(dāng)前觀測(cè)值的質(zhì)量控制依據(jù)。整個(gè)計(jì)算過程平滑滾動(dòng),其結(jié)果將隨學(xué)習(xí)過程而作動(dòng)態(tài)調(diào)整,便于實(shí)時(shí)業(yè)務(wù)的應(yīng)用展開。

    2.5 概率落點(diǎn)技術(shù)

    針對(duì)各種誤差所帶入的觀測(cè)要素,相關(guān)質(zhì)量控制技術(shù)在國內(nèi)外也做過大量的研究如空間一致性、時(shí)間一致性等均有不同程度的質(zhì)量控制效果。本文針對(duì)相關(guān)地區(qū)的區(qū)域站在特殊環(huán)境,設(shè)計(jì)較短時(shí)序的概率落點(diǎn)技術(shù),即統(tǒng)計(jì)分析較短時(shí)序下的觀測(cè)要素的相鄰觀測(cè)值相關(guān)性,構(gòu)建落點(diǎn)統(tǒng)計(jì)矩陣,獲取相應(yīng)落點(diǎn)概率,以此作為判斷短時(shí)效內(nèi)觀測(cè)要素的變動(dòng)范圍及質(zhì)控依據(jù)。

    記ti時(shí)次的觀測(cè)值對(duì)應(yīng)ti+1時(shí)次的觀測(cè)值為落點(diǎn)C[X(tk)],相同觀測(cè)值的相同落點(diǎn)計(jì)數(shù)加1,累積落點(diǎn)值∑X(tki),得到落點(diǎn)統(tǒng)計(jì)矩陣:

    針對(duì)落地統(tǒng)計(jì)矩陣,到落點(diǎn)概率分布:

    其中,p(ki)=∑X(tki)/(∑X(tk1)+…+∑X(tkn))。

    當(dāng)C[X(tk)]是ti時(shí)次的觀測(cè)值對(duì)應(yīng)ti+2時(shí)次的觀測(cè)值時(shí),稱為二階落點(diǎn),并對(duì)應(yīng)相應(yīng)的二階落點(diǎn)矩陣和二階落點(diǎn)概率分布,依此類推,本文暫不贅述。

    經(jīng)過一定時(shí)間的算法訓(xùn)練后,落點(diǎn)概率分布將趨于穩(wěn)態(tài),學(xué)習(xí)時(shí)間越長,分布越穩(wěn)定,但同時(shí)落點(diǎn)也變得更加離散,擴(kuò)大判斷范圍,不利于質(zhì)量控制應(yīng)用。本文采用實(shí)際數(shù)據(jù)作為算法學(xué)習(xí)訓(xùn)練依據(jù),不同的觀測(cè)要素存在相應(yīng)的最佳時(shí)間窗口。

    從上述雷達(dá)回波過程中選取一塊區(qū)域進(jìn)行研究,以F9742對(duì)應(yīng)的降水量要素落點(diǎn)進(jìn)行研究,過程描述如圖1所示。

    (a)8日降水落點(diǎn)

    計(jì)算表明,降水在8~9日之間有切變,與實(shí)際天氣過程較吻合;氣溫變化穩(wěn)定,落點(diǎn)概率平穩(wěn),表明預(yù)測(cè)參考前一值(或前一天的當(dāng)前時(shí)刻值)有一定意義;氣壓落點(diǎn)概率固定,后經(jīng)數(shù)據(jù)驗(yàn)證為氣壓計(jì)輸出值長時(shí)段內(nèi)不改變。

    3 結(jié)果分析

    融合上述幾種計(jì)算法,我們針對(duì)每個(gè)站,利用一定時(shí)序的數(shù)據(jù),預(yù)測(cè)下一個(gè)時(shí)次的數(shù)值,并結(jié)合實(shí)際值進(jìn)行對(duì)比,觀測(cè)與預(yù)測(cè)對(duì)比如圖2所示。

    圖2 觀測(cè)與預(yù)測(cè)對(duì)比過程(實(shí)際數(shù)據(jù)分析)

    觀測(cè)值與動(dòng)態(tài)預(yù)測(cè)值的相關(guān)系數(shù)為:

    觀測(cè)值與動(dòng)態(tài)預(yù)測(cè)值的均方差為:

    圖3 觀測(cè)與預(yù)測(cè)對(duì)比過程(模擬異常值)

    此時(shí),觀測(cè)值與動(dòng)態(tài)預(yù)測(cè)值的相關(guān)性系數(shù)為:Correl(X,Y)=0.293 ,觀測(cè)值與動(dòng)態(tài)預(yù)測(cè)值的均方差為Stdev(X,Y)=7.01。這反映出相關(guān)性變低,兩者之間表現(xiàn)出背離。從圖3對(duì)比結(jié)果來看,預(yù)測(cè)值能夠較好符合實(shí)時(shí)觀測(cè)值,當(dāng)出現(xiàn)顯性異常時(shí),相關(guān)性背離加快,這在實(shí)時(shí)業(yè)務(wù)應(yīng)用中具備了風(fēng)險(xiǎn)研判和數(shù)據(jù)質(zhì)量訂正依據(jù)。

    4 結(jié)論

    自動(dòng)站數(shù)據(jù)在災(zāi)害性天氣監(jiān)測(cè)預(yù)警中應(yīng)用最廣泛、最直接,其數(shù)據(jù)質(zhì)量對(duì)于預(yù)報(bào)分析和服務(wù)效果有直接影響。目前也有很多數(shù)據(jù)質(zhì)控平臺(tái)專門對(duì)數(shù)據(jù)進(jìn)行研究與應(yīng)用,且也已經(jīng)實(shí)現(xiàn)了快速質(zhì)控、人工疑誤處理等多種處理方式。然而在實(shí)際業(yè)務(wù)中,仍然存在各式各樣的要素冒“大數(shù)”的現(xiàn)象,如雨量、風(fēng)速、溫度等。本文主要從實(shí)際業(yè)務(wù)出發(fā),業(yè)務(wù)數(shù)據(jù)入口之前通過一定的技術(shù)研究增加數(shù)據(jù)質(zhì)量防火墻或預(yù)警墻,并取得了一定的預(yù)期效果。本文所提的計(jì)算法還有不斷完善和改進(jìn)的空間,需要繼續(xù)探索國內(nèi)外技術(shù),以期繼續(xù)豐富數(shù)據(jù)的質(zhì)量控制技術(shù)。

    猜你喜歡
    游程落點(diǎn)觀測(cè)
    觀測(cè)到恒星死亡瞬間
    軍事文摘(2023年18期)2023-11-03 09:45:42
    基于劃分組參考數(shù)的差值編碼壓縮方法
    中國羽毛球組合鄭思維/黃雅瓊連續(xù)得失分規(guī)律研究
    基于空間分層組合設(shè)計(jì)的火箭落點(diǎn)實(shí)時(shí)計(jì)算模型
    改進(jìn)型相對(duì)游程長度編碼方法
    天測(cè)與測(cè)地VLBI 測(cè)地站周圍地形觀測(cè)遮掩的討論
    美火星軌道器拍到歐洲著陸器落點(diǎn)圖像
    太空探索(2016年12期)2016-07-18 11:13:43
    拼搶第二落點(diǎn)新聞打好新聞競爭戰(zhàn)
    新聞傳播(2016年4期)2016-07-18 10:59:21
    可觀測(cè)宇宙
    太空探索(2016年7期)2016-07-10 12:10:15
    探訪江蘇地方立法:百姓關(guān)切成立法落點(diǎn)
    欧美 日韩 精品 国产| 人妻人人澡人人爽人人| 青青草视频在线视频观看| 亚洲国产精品一区三区| 亚洲伊人久久精品综合| 国产一区二区在线观看av| 婷婷色综合大香蕉| 久久这里只有精品19| 男女边摸边吃奶| 天天操日日干夜夜撸| 亚洲精品久久久久久婷婷小说| 91国产中文字幕| 亚洲精品,欧美精品| 一级,二级,三级黄色视频| av免费观看日本| 超色免费av| 日本与韩国留学比较| 国产极品粉嫩免费观看在线| 亚洲欧美成人综合另类久久久| 日韩三级伦理在线观看| 男女高潮啪啪啪动态图| 久久99热这里只频精品6学生| 一二三四中文在线观看免费高清| 久久久久国产精品人妻一区二区| 青春草亚洲视频在线观看| 成人漫画全彩无遮挡| 国产免费一区二区三区四区乱码| 一区二区三区乱码不卡18| 久久久久网色| 欧美3d第一页| 老司机亚洲免费影院| 精品少妇内射三级| 中文字幕精品免费在线观看视频 | 少妇猛男粗大的猛烈进出视频| 啦啦啦中文免费视频观看日本| 免费观看无遮挡的男女| 成人影院久久| 日本爱情动作片www.在线观看| 91精品伊人久久大香线蕉| av卡一久久| 精品久久久精品久久久| 欧美日韩av久久| 免费黄色在线免费观看| 欧美成人午夜精品| 中文字幕制服av| 国产亚洲精品第一综合不卡 | 少妇人妻 视频| 国产精品国产三级国产av玫瑰| 99热网站在线观看| 有码 亚洲区| 久久午夜综合久久蜜桃| 国精品久久久久久国模美| 99九九在线精品视频| 一级毛片我不卡| 日本欧美国产在线视频| 久久精品久久久久久久性| 日韩制服丝袜自拍偷拍| 成人综合一区亚洲| 亚洲精华国产精华液的使用体验| 日韩视频在线欧美| 男人添女人高潮全过程视频| 在线观看一区二区三区激情| 人妻 亚洲 视频| 18禁国产床啪视频网站| 欧美激情国产日韩精品一区| 老司机影院毛片| 精品福利永久在线观看| 国产精品无大码| 人人妻人人澡人人看| 黑人高潮一二区| 两个人看的免费小视频| 少妇猛男粗大的猛烈进出视频| 亚洲综合精品二区| 日韩一区二区三区影片| 国产欧美日韩一区二区三区在线| 久久婷婷青草| 久久精品国产亚洲av涩爱| 色94色欧美一区二区| 中文字幕人妻丝袜制服| 亚洲av国产av综合av卡| 五月伊人婷婷丁香| 亚洲精品乱久久久久久| 成人18禁高潮啪啪吃奶动态图| 在线观看免费高清a一片| 九色亚洲精品在线播放| 97在线人人人人妻| 亚洲国产精品999| 中文字幕另类日韩欧美亚洲嫩草| 大香蕉久久成人网| 最黄视频免费看| 五月玫瑰六月丁香| 欧美日韩视频高清一区二区三区二| 少妇猛男粗大的猛烈进出视频| 日韩欧美精品免费久久| 一边摸一边做爽爽视频免费| 日本91视频免费播放| 亚洲欧美色中文字幕在线| 免费在线观看完整版高清| 亚洲av日韩在线播放| 中文字幕人妻熟女乱码| 桃花免费在线播放| 日日摸夜夜添夜夜爱| av女优亚洲男人天堂| 十八禁高潮呻吟视频| 亚洲,欧美精品.| 在线天堂最新版资源| 一区在线观看完整版| 免费久久久久久久精品成人欧美视频 | 国产爽快片一区二区三区| 成年动漫av网址| 亚洲第一区二区三区不卡| 男人操女人黄网站| 欧美另类一区| 亚洲综合精品二区| 欧美成人精品欧美一级黄| 国产日韩欧美在线精品| 国产毛片在线视频| 免费女性裸体啪啪无遮挡网站| 亚洲天堂av无毛| 少妇的逼好多水| 午夜福利视频精品| 最近的中文字幕免费完整| 国产亚洲av片在线观看秒播厂| 久久这里有精品视频免费| 国产精品久久久久久久电影| av播播在线观看一区| 国产精品国产三级国产av玫瑰| 欧美日韩视频高清一区二区三区二| 国产精品久久久久成人av| 性色av一级| 一个人免费看片子| 亚洲性久久影院| 久久久久久久久久人人人人人人| 精品福利永久在线观看| 日本wwww免费看| 久久久久人妻精品一区果冻| 成人无遮挡网站| 国产精品一二三区在线看| 巨乳人妻的诱惑在线观看| 亚洲国产日韩一区二区| 欧美精品人与动牲交sv欧美| 又大又黄又爽视频免费| 黑人猛操日本美女一级片| 精品国产国语对白av| 免费黄色在线免费观看| 国产免费又黄又爽又色| 波野结衣二区三区在线| 18在线观看网站| 久久精品久久久久久噜噜老黄| 欧美精品av麻豆av| 99久国产av精品国产电影| 亚洲国产欧美日韩在线播放| 亚洲国产毛片av蜜桃av| 久久女婷五月综合色啪小说| 2022亚洲国产成人精品| 日韩三级伦理在线观看| 91精品国产国语对白视频| 少妇 在线观看| 精品99又大又爽又粗少妇毛片| 欧美+日韩+精品| 亚洲av日韩在线播放| xxx大片免费视频| 丝袜美足系列| 啦啦啦啦在线视频资源| 精品视频人人做人人爽| 你懂的网址亚洲精品在线观看| 亚洲熟女精品中文字幕| 免费av中文字幕在线| 一本久久精品| 九九爱精品视频在线观看| 免费在线观看完整版高清| 男女啪啪激烈高潮av片| 亚洲成国产人片在线观看| 久久韩国三级中文字幕| 性色av一级| 一区二区日韩欧美中文字幕 | 菩萨蛮人人尽说江南好唐韦庄| 人人妻人人添人人爽欧美一区卜| 亚洲欧美日韩另类电影网站| 一级爰片在线观看| 美女脱内裤让男人舔精品视频| 成人黄色视频免费在线看| 搡老乐熟女国产| 久久久精品区二区三区| 午夜福利,免费看| 欧美日韩成人在线一区二区| 色哟哟·www| 色婷婷久久久亚洲欧美| 久久人人97超碰香蕉20202| 成人免费观看视频高清| 日本午夜av视频| 免费大片18禁| 免费看av在线观看网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日韩精品免费视频一区二区三区 | 一区二区av电影网| 黄色怎么调成土黄色| 久久人人爽人人片av| 80岁老熟妇乱子伦牲交| 满18在线观看网站| 亚洲综合色惰| 久久国产亚洲av麻豆专区| 9热在线视频观看99| 免费黄频网站在线观看国产| 成人免费观看视频高清| 久久精品国产自在天天线| 9色porny在线观看| 欧美日韩av久久| tube8黄色片| 全区人妻精品视频| 成人国产麻豆网| 免费观看性生交大片5| 欧美 日韩 精品 国产| 国产伦理片在线播放av一区| 亚洲av国产av综合av卡| 男女国产视频网站| 亚洲国产精品成人久久小说| 国产一区有黄有色的免费视频| 色94色欧美一区二区| 国产国语露脸激情在线看| 午夜免费鲁丝| 国产精品久久久久久久久免| 久久精品国产亚洲av涩爱| 国产 一区精品| 黄色怎么调成土黄色| 国产精品熟女久久久久浪| 一个人免费看片子| 成年美女黄网站色视频大全免费| 免费久久久久久久精品成人欧美视频 | 黄色毛片三级朝国网站| 黄色 视频免费看| 蜜桃在线观看..| 男女边摸边吃奶| www.熟女人妻精品国产 | 日产精品乱码卡一卡2卡三| 欧美少妇被猛烈插入视频| 最近2019中文字幕mv第一页| 日本欧美视频一区| 天堂中文最新版在线下载| 亚洲av福利一区| 国产精品一区www在线观看| 亚洲精华国产精华液的使用体验| 超色免费av| 亚洲国产精品一区三区| 亚洲成人一二三区av| 欧美日本中文国产一区发布| 91aial.com中文字幕在线观看| 亚洲,一卡二卡三卡| 久久久久久人人人人人| 国产精品不卡视频一区二区| 亚洲国产欧美在线一区| 日日摸夜夜添夜夜爱| 国产精品久久久久久久久免| 精品少妇黑人巨大在线播放| 成人手机av| 天堂8中文在线网| 飞空精品影院首页| 黄色配什么色好看| 日本与韩国留学比较| 久久 成人 亚洲| 黑人巨大精品欧美一区二区蜜桃 | 中文精品一卡2卡3卡4更新| 国产老妇伦熟女老妇高清| 男女无遮挡免费网站观看| 男女边摸边吃奶| 国产av国产精品国产| 国产爽快片一区二区三区| 你懂的网址亚洲精品在线观看| 2018国产大陆天天弄谢| 久久久久精品性色| 国产成人午夜福利电影在线观看| 日本av免费视频播放| 久久久久久久精品精品| 九九爱精品视频在线观看| 久久久久久久久久人人人人人人| 亚洲av福利一区| 韩国高清视频一区二区三区| 亚洲在久久综合| 美女主播在线视频| 9191精品国产免费久久| 巨乳人妻的诱惑在线观看| 中文天堂在线官网| 久久ye,这里只有精品| 香蕉精品网在线| 国产av码专区亚洲av| 高清av免费在线| 曰老女人黄片| 日韩在线高清观看一区二区三区| 国产精品成人在线| 亚洲精品久久午夜乱码| 自线自在国产av| 成人毛片60女人毛片免费| 成人18禁高潮啪啪吃奶动态图| 婷婷成人精品国产| 老司机影院毛片| 国产国拍精品亚洲av在线观看| 青春草国产在线视频| 国产免费现黄频在线看| 爱豆传媒免费全集在线观看| 侵犯人妻中文字幕一二三四区| 亚洲国产精品999| 26uuu在线亚洲综合色| videossex国产| 在线观看一区二区三区激情| 一级片免费观看大全| 天美传媒精品一区二区| 欧美成人午夜免费资源| 夫妻性生交免费视频一级片| 国产成人免费观看mmmm| 久久久久精品性色| 色5月婷婷丁香| 国产欧美日韩综合在线一区二区| 中国国产av一级| 天堂俺去俺来也www色官网| www.av在线官网国产| 成人国产麻豆网| 中文乱码字字幕精品一区二区三区| 日韩在线高清观看一区二区三区| 只有这里有精品99| 校园人妻丝袜中文字幕| 国产亚洲av片在线观看秒播厂| 欧美精品人与动牲交sv欧美| 观看美女的网站| 国产高清三级在线| av在线观看视频网站免费| 一区二区三区四区激情视频| 免费少妇av软件| 亚洲美女黄色视频免费看| 亚洲综合精品二区| 另类亚洲欧美激情| 日本av免费视频播放| 亚洲欧美清纯卡通| 亚洲精品久久成人aⅴ小说| 亚洲精品,欧美精品| 王馨瑶露胸无遮挡在线观看| 午夜福利乱码中文字幕| 中文欧美无线码| 18禁在线无遮挡免费观看视频| 成人综合一区亚洲| 日本色播在线视频| 高清欧美精品videossex| 国产不卡av网站在线观看| 精品国产一区二区久久| 欧美性感艳星| 国产爽快片一区二区三区| 少妇高潮的动态图| 又粗又硬又长又爽又黄的视频| 另类精品久久| 在线天堂最新版资源| 男女国产视频网站| 三级国产精品片| 人体艺术视频欧美日本| 制服人妻中文乱码| av在线老鸭窝| 中文欧美无线码| 欧美日韩av久久| 亚洲经典国产精华液单| www.熟女人妻精品国产 | 亚洲美女黄色视频免费看| 夜夜骑夜夜射夜夜干| 精品少妇黑人巨大在线播放| 国产精品人妻久久久影院| 婷婷成人精品国产| 免费观看在线日韩| 啦啦啦中文免费视频观看日本| 男人添女人高潮全过程视频| 高清黄色对白视频在线免费看| 一区二区三区四区激情视频| 高清在线视频一区二区三区| 在线天堂最新版资源| 国产日韩一区二区三区精品不卡| 97在线人人人人妻| 国产极品天堂在线| 一本大道久久a久久精品| 亚洲成av片中文字幕在线观看 | a级毛片在线看网站| 精品少妇内射三级| 亚洲人成网站在线观看播放| 精品一区二区免费观看| a级毛色黄片| 午夜福利网站1000一区二区三区| 王馨瑶露胸无遮挡在线观看| 熟女av电影| 欧美精品国产亚洲| 有码 亚洲区| 黄网站色视频无遮挡免费观看| 免费观看在线日韩| 日日摸夜夜添夜夜爱| 午夜视频国产福利| 日本wwww免费看| 欧美3d第一页| 狠狠婷婷综合久久久久久88av| 国产激情久久老熟女| 亚洲国产欧美日韩在线播放| 亚洲av欧美aⅴ国产| 王馨瑶露胸无遮挡在线观看| 熟女av电影| 亚洲av.av天堂| 亚洲av综合色区一区| 大话2 男鬼变身卡| 男男h啪啪无遮挡| 成人午夜精彩视频在线观看| 久久久久精品久久久久真实原创| 亚洲精品第二区| 肉色欧美久久久久久久蜜桃| 国产日韩一区二区三区精品不卡| av黄色大香蕉| 1024视频免费在线观看| 亚洲精品国产色婷婷电影| 高清毛片免费看| 亚洲精品日韩在线中文字幕| 一边亲一边摸免费视频| 香蕉精品网在线| 一级,二级,三级黄色视频| 国产精品嫩草影院av在线观看| 国产av精品麻豆| 国产1区2区3区精品| 亚洲精品一二三| 2018国产大陆天天弄谢| 97在线人人人人妻| 欧美人与性动交α欧美精品济南到 | 成年美女黄网站色视频大全免费| 免费观看无遮挡的男女| 亚洲,欧美精品.| 九色成人免费人妻av| av线在线观看网站| av播播在线观看一区| 国产亚洲精品久久久com| 少妇高潮的动态图| 男女高潮啪啪啪动态图| av在线app专区| 少妇被粗大猛烈的视频| 九色亚洲精品在线播放| 全区人妻精品视频| av国产精品久久久久影院| 午夜视频国产福利| 各种免费的搞黄视频| 久久久久久久大尺度免费视频| 天天操日日干夜夜撸| 久久国产亚洲av麻豆专区| 日韩在线高清观看一区二区三区| 啦啦啦在线观看免费高清www| 纯流量卡能插随身wifi吗| 日产精品乱码卡一卡2卡三| 亚洲国产色片| 99精国产麻豆久久婷婷| 亚洲欧美日韩卡通动漫| 日韩在线高清观看一区二区三区| 国产国拍精品亚洲av在线观看| 国产日韩欧美亚洲二区| 日韩中字成人| 性色av一级| 天天躁夜夜躁狠狠久久av| 国内精品宾馆在线| 国产男女超爽视频在线观看| 秋霞伦理黄片| 久久久久久伊人网av| 中文乱码字字幕精品一区二区三区| 国产午夜精品一二区理论片| 日本wwww免费看| 久久久久久久大尺度免费视频| 这个男人来自地球电影免费观看 | 宅男免费午夜| 精品一区二区免费观看| 亚洲精品中文字幕在线视频| 天天躁夜夜躁狠狠久久av| 性色avwww在线观看| 女性生殖器流出的白浆| 久久综合国产亚洲精品| 亚洲丝袜综合中文字幕| 亚洲国产精品专区欧美| 欧美成人精品欧美一级黄| 1024视频免费在线观看| 男男h啪啪无遮挡| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 免费看光身美女| 搡老乐熟女国产| 亚洲精品久久成人aⅴ小说| 亚洲av.av天堂| 如何舔出高潮| 中文欧美无线码| 亚洲欧美日韩另类电影网站| 国产成人精品无人区| 一本大道久久a久久精品| 日韩熟女老妇一区二区性免费视频| 久久人人爽av亚洲精品天堂| 考比视频在线观看| 大片电影免费在线观看免费| 亚洲美女搞黄在线观看| 欧美日韩视频精品一区| 免费黄色在线免费观看| 母亲3免费完整高清在线观看 | 欧美国产精品一级二级三级| 国产一区有黄有色的免费视频| 草草在线视频免费看| 国产爽快片一区二区三区| 日韩中字成人| 黄色 视频免费看| 哪个播放器可以免费观看大片| 青春草国产在线视频| 18在线观看网站| 丝袜脚勾引网站| 永久网站在线| 老司机亚洲免费影院| 午夜福利在线观看免费完整高清在| 精品国产国语对白av| 免费高清在线观看日韩| 亚洲,一卡二卡三卡| 国产成人精品一,二区| 这个男人来自地球电影免费观看 | 十八禁网站网址无遮挡| 人人澡人人妻人| 久久这里有精品视频免费| 亚洲内射少妇av| 成年动漫av网址| 国产亚洲欧美精品永久| 亚洲欧美色中文字幕在线| 亚洲精品乱久久久久久| 男女午夜视频在线观看 | 久久久久久久大尺度免费视频| 日本av免费视频播放| 青春草国产在线视频| 亚洲激情五月婷婷啪啪| 久久人人爽人人片av| 丝袜人妻中文字幕| 熟女人妻精品中文字幕| a级毛片在线看网站| 97人妻天天添夜夜摸| 亚洲av欧美aⅴ国产| 九色亚洲精品在线播放| 国产福利在线免费观看视频| 欧美国产精品va在线观看不卡| 国产熟女欧美一区二区| 国产免费一区二区三区四区乱码| 久久久久久久久久久免费av| av卡一久久| 国产av精品麻豆| 欧美 亚洲 国产 日韩一| 亚洲av日韩在线播放| 亚洲国产最新在线播放| 大陆偷拍与自拍| 性色avwww在线观看| 99视频精品全部免费 在线| 春色校园在线视频观看| 国产熟女欧美一区二区| 高清欧美精品videossex| 永久网站在线| 亚洲精品日本国产第一区| 亚洲av日韩在线播放| 丝袜脚勾引网站| 边亲边吃奶的免费视频| 九色亚洲精品在线播放| 1024视频免费在线观看| 99热这里只有是精品在线观看| 2022亚洲国产成人精品| 青春草国产在线视频| 久久99蜜桃精品久久| 国产成人欧美| 亚洲欧洲国产日韩| 国产爽快片一区二区三区| 少妇人妻久久综合中文| 精品亚洲乱码少妇综合久久| 七月丁香在线播放| 日本欧美视频一区| 亚洲精品乱久久久久久| 国产色婷婷99| 黑人猛操日本美女一级片| 国产成人午夜福利电影在线观看| 狂野欧美激情性xxxx在线观看| 久久av网站| 国产1区2区3区精品| 精品国产一区二区三区久久久樱花| 午夜91福利影院| 日本免费在线观看一区| av女优亚洲男人天堂| 男女啪啪激烈高潮av片| 青春草国产在线视频| 久久99蜜桃精品久久| 国产成人欧美| 精品国产一区二区久久| 欧美成人午夜免费资源| 国产精品人妻久久久影院| 在线观看免费视频网站a站| 久久ye,这里只有精品| 欧美激情国产日韩精品一区| 国产成人免费无遮挡视频| 国产av码专区亚洲av| 亚洲av免费高清在线观看| 日韩一本色道免费dvd| 啦啦啦啦在线视频资源| 婷婷色综合大香蕉| 亚洲精品一区蜜桃| 日韩三级伦理在线观看| 亚洲国产看品久久| 91久久精品国产一区二区三区| 欧美人与善性xxx| 99久久综合免费| 国产成人精品无人区| 青春草国产在线视频| 日本午夜av视频| 男人舔女人的私密视频| 十分钟在线观看高清视频www| 久久狼人影院| 久久人人爽av亚洲精品天堂| 国产精品.久久久| 咕卡用的链子| 永久网站在线| 韩国av在线不卡| 久久久久精品久久久久真实原创| 又黄又爽又刺激的免费视频.| 一级片'在线观看视频| 人妻一区二区av| 99久久综合免费| 日韩一本色道免费dvd| 国产成人欧美|