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

    臺(tái)風(fēng)三維風(fēng)場(chǎng)數(shù)據(jù)處理方法研究

    2024-12-06 00:00:00練江峰
    中國新技術(shù)新產(chǎn)品 2024年15期
    關(guān)鍵詞:插值法數(shù)據(jù)處理臺(tái)風(fēng)

    摘 要:合理地處理脈動(dòng)風(fēng)場(chǎng)風(fēng)速實(shí)測(cè)數(shù)據(jù)是獲取準(zhǔn)確脈動(dòng)風(fēng)場(chǎng)參數(shù)的基礎(chǔ)。為了提升三維脈動(dòng)風(fēng)場(chǎng)實(shí)測(cè)數(shù)據(jù)分析結(jié)果的準(zhǔn)確性和合理性,本文基于超聲風(fēng)速儀無效風(fēng)速數(shù)據(jù)自動(dòng)識(shí)別提出了一種基于多倍截?cái)喾讲罘ǖ拿}動(dòng)風(fēng)速實(shí)測(cè)數(shù)據(jù)處理方法,開展了基于實(shí)測(cè)風(fēng)場(chǎng)特征分析的數(shù)據(jù)修正工作,并利用臺(tái)風(fēng)“瑪莉亞”實(shí)測(cè)脈動(dòng)風(fēng)速數(shù)據(jù)驗(yàn)證了該方法的有效性。該方法可以為風(fēng)場(chǎng)實(shí)測(cè)脈動(dòng)風(fēng)數(shù)據(jù)處理以及風(fēng)障等構(gòu)造研發(fā)中風(fēng)場(chǎng)精細(xì)化模擬提供參考。

    關(guān)鍵詞:臺(tái)風(fēng);數(shù)據(jù)處理;脈動(dòng)風(fēng)速;插值法

    中圖分類號(hào):P 44 文獻(xiàn)標(biāo)志碼:A

    臺(tái)風(fēng)是風(fēng)災(zāi)害的主要形式之一,屬于極端風(fēng)氣候,其造成的受災(zāi)范圍廣,發(fā)生頻率高。隨著我國柔性橋梁結(jié)構(gòu)的建設(shè)日益增多,例如在本公司承接的濱海西大道提升改造工程中,后田人行天橋采用雙肢外斜變寬鋼斜塔斜拉橋,海滄鰲冠大道工程采用雙層橋梁,上層為現(xiàn)澆連續(xù)箱梁、下層為懸掛式鋼箱梁漂浮體系結(jié)構(gòu)。為了保障極端風(fēng)氣候下橋面行車及人行安全,亟需加強(qiáng)橋梁風(fēng)障等構(gòu)造阻風(fēng)性能研究[1],而獲取準(zhǔn)確臺(tái)風(fēng)風(fēng)場(chǎng)參數(shù)是實(shí)現(xiàn)風(fēng)障等構(gòu)造風(fēng)場(chǎng)精細(xì)化模擬的基礎(chǔ)性工作。作為獲取復(fù)雜天氣三維脈動(dòng)風(fēng)速的有效手段,風(fēng)場(chǎng)實(shí)測(cè)被國內(nèi)外學(xué)者廣泛應(yīng)用于風(fēng)工程研究中[2-3]。由于實(shí)測(cè)三維風(fēng)速是隨時(shí)間而變化的動(dòng)態(tài)數(shù)據(jù),時(shí)序數(shù)據(jù)往往受儀器設(shè)計(jì)缺陷、暴雨、雷擊、風(fēng)致飛射物沖擊、電壓不穩(wěn)定以及大氣湍流隨機(jī)性等影響,存在數(shù)據(jù)壞點(diǎn)或不合理點(diǎn),從而影響實(shí)測(cè)數(shù)據(jù)的準(zhǔn)確性。因此對(duì)觀測(cè)風(fēng)場(chǎng)數(shù)據(jù)進(jìn)行處理、判定觀測(cè)數(shù)據(jù)的可靠性是準(zhǔn)確描述臺(tái)風(fēng)場(chǎng)特性的重要基礎(chǔ)。本文基于超聲風(fēng)速儀無效風(fēng)速數(shù)據(jù)自動(dòng)識(shí)別,探討了臺(tái)風(fēng)風(fēng)場(chǎng)數(shù)據(jù)預(yù)處理方法,開展了基于實(shí)測(cè)風(fēng)場(chǎng)特征分析的數(shù)據(jù)修正工作,并利用實(shí)測(cè)的臺(tái)風(fēng)“瑪莉亞”實(shí)測(cè)脈動(dòng)風(fēng)速數(shù)據(jù)驗(yàn)證了此方法的有效性。該方法可為風(fēng)場(chǎng)實(shí)測(cè)脈動(dòng)風(fēng)數(shù)據(jù)處理以及風(fēng)障等構(gòu)造風(fēng)場(chǎng)精細(xì)化模擬提供參考。

    1 實(shí)測(cè)風(fēng)場(chǎng)數(shù)據(jù)來源

    1.1 觀測(cè)站

    觀測(cè)站位于福建省平潭縣嶼頭島,觀測(cè)站主要包括觀測(cè)系統(tǒng)及數(shù)據(jù)采集系統(tǒng)。觀測(cè)系統(tǒng)主要采用英國Gill公司生產(chǎn)的高性能三維超聲風(fēng)速儀Windmaster Pro作為風(fēng)數(shù)據(jù)采樣儀器,分別布置 10m、80m、100 m三個(gè)觀測(cè)層,該風(fēng)速儀風(fēng)速量程為0m/s~65m/s,分辨率為0.01m/s;風(fēng)向量程為0°~359°,分辨率為0.1°;支持32Hz高頻采樣,然而,過高的采樣頻率并不一定能有較好的還原度,反而會(huì)使接近電壓頻率的信號(hào)折疊引入數(shù)據(jù)低頻端,導(dǎo)致數(shù)據(jù)量過大,存在一定的誤差,因而現(xiàn)場(chǎng)選用的樣本輸出頻率為10Hz,即每0.1s記錄一組u(x)、u(y)、u(z)三維瞬時(shí)風(fēng)速。

    數(shù)據(jù)采集系統(tǒng)采用CR3000作為數(shù)據(jù)采集器,該設(shè)備性能可靠,運(yùn)行穩(wěn)定,達(dá)到歐盟CE、EMC要求,能有效防止過壓下強(qiáng)電流對(duì)儀器的損傷,CR3000具有CSI/O和RS-232接口,支持SDM外圍設(shè)備,可以直接與計(jì)算機(jī)或PDA連接,也能夠利用Internet、無線電等通信方式進(jìn)行數(shù)據(jù)無線傳輸,以實(shí)現(xiàn)臺(tái)風(fēng)期間連續(xù)性觀測(cè)及數(shù)據(jù)實(shí)時(shí)傳輸?shù)墓δ堋?/p>

    1.2 實(shí)測(cè)臺(tái)風(fēng)

    2018年第8號(hào)臺(tái)風(fēng)“瑪莉亞”于7月11日05:00位于福建霞浦東偏南方大約120km的海面上,北緯26.3°、東經(jīng)121.0°,中心附近最大風(fēng)力有十五級(jí)(48m/s),中心最低氣壓為94500Pa,七級(jí)風(fēng)圈半徑280km~400km,十級(jí)風(fēng)圈半徑150km~200km,十二級(jí)風(fēng)圈半徑40km~50km。并于7月11日09時(shí)10分在福建兩江黃岐半島登陸,登陸時(shí)附近最大風(fēng)力十四級(jí)(42m/s),中心最低氣壓96000Pa。臺(tái)風(fēng)“瑪莉亞”的登陸路徑與觀測(cè)站最短距離為18.1km。

    1.3 臺(tái)風(fēng)數(shù)據(jù)初步分析

    獲取可靠的、具有代表性的風(fēng)觀測(cè)數(shù)據(jù)是風(fēng)場(chǎng)研究的重要工作。對(duì)結(jié)構(gòu)風(fēng)工程來說,風(fēng)觀測(cè)數(shù)據(jù)的精度要求往往比一般的區(qū)域氣象預(yù)報(bào)要高,因此對(duì)一手的實(shí)測(cè)風(fēng)數(shù)據(jù)進(jìn)行初步分析預(yù)處理是科學(xué)開展結(jié)構(gòu)風(fēng)工程抗風(fēng)研究的必要過程,但影響風(fēng)場(chǎng)數(shù)據(jù)精度因素較多,為了完成針對(duì)性降低誤差的工作,利用臺(tái)風(fēng)“瑪莉亞”期間觀測(cè)站10m、80m高度層風(fēng)場(chǎng)三維風(fēng)速數(shù)據(jù),對(duì)基于超聲風(fēng)速儀無效風(fēng)速數(shù)據(jù)自動(dòng)識(shí)別的數(shù)據(jù)進(jìn)行初步處理,得到“瑪莉亞”期間10m、80m處ux、uy、uz風(fēng)速序列如圖1所示。

    由圖1可知,臺(tái)風(fēng)“瑪莉亞”期間,ux、uy、uz部分最大值接近100m/s,這主要是臺(tái)風(fēng)期間受環(huán)境等因素的影響,導(dǎo)致部分?jǐn)?shù)據(jù)失真,產(chǎn)生較多野點(diǎn),發(fā)生“毛刺”現(xiàn)象。這將降低實(shí)測(cè)數(shù)據(jù)的有效性,在某種程度上,也會(huì)使臺(tái)風(fēng)期間實(shí)測(cè)平均風(fēng)速數(shù)值出現(xiàn)不真實(shí)的“虛大”情況,如果將不真實(shí)的實(shí)測(cè)結(jié)果應(yīng)用于橋梁等結(jié)構(gòu)的設(shè)計(jì)基準(zhǔn)風(fēng)速等參數(shù)評(píng)估上,勢(shì)必會(huì)影響結(jié)構(gòu)設(shè)計(jì)參數(shù)取值的合理性和科學(xué)性,為了避免發(fā)生上述情況,有必要在上述超聲風(fēng)速儀無效風(fēng)速數(shù)據(jù)自動(dòng)識(shí)別數(shù)據(jù)特征分析的基礎(chǔ)上,進(jìn)一步修正臺(tái)風(fēng)期間實(shí)測(cè)風(fēng)速。

    2 數(shù)據(jù)處理方法

    為了提升三維脈動(dòng)風(fēng)場(chǎng)實(shí)測(cè)數(shù)據(jù)分析結(jié)果的準(zhǔn)確性和合理性,降低實(shí)測(cè)數(shù)據(jù)“毛刺”現(xiàn)象對(duì)結(jié)果產(chǎn)生的不利影響,基于多倍截?cái)喾讲罘▽?duì)上述特征分析中“壞點(diǎn)”數(shù)據(jù)或者不合理數(shù)據(jù)進(jìn)行逐個(gè)時(shí)間序列的原始數(shù)據(jù)平滑估計(jì),通過檢測(cè)數(shù)據(jù)突變位置,確定該點(diǎn)是否為平滑估計(jì)范圍內(nèi)的為正常點(diǎn),其處理過程如公式(1)~公式(5)所示[4-6]。

    du(t)=u(t+2)-u(t) (1)

    式中:du(t)為構(gòu)建的時(shí)間序列;t為樣本序列號(hào);u(t)為樣本系列中第t個(gè)樣本風(fēng)速;u(t+2)為樣本系列中第t+2個(gè)樣本風(fēng)速。

    (2)

    式中:為du(t)時(shí)間序列的平均值;n為時(shí)間序列的樣本數(shù)。

    (3)

    式中:為du(t)2時(shí)間序列的平均值;n為時(shí)間序列的樣本數(shù)。

    σ=-2 (4)

    式中:σ為截?cái)喾讲睢?/p>

    Δ=c·σ0.5 (5)

    式中:Δ為數(shù)據(jù)壞點(diǎn)判斷標(biāo)準(zhǔn)參數(shù);c為常系數(shù)。

    通過系數(shù)c可以設(shè)置正常數(shù)據(jù)偏離平滑估計(jì)的范圍,這里取c=4,當(dāng)滿足|du(t)|gt;Δ或|du(t+2)|gt;Δ時(shí),即該樣本點(diǎn)與總體樣本均值的差的絕對(duì)值大于4倍標(biāo)準(zhǔn)差時(shí),則判斷u(t+2)是不合理的壞點(diǎn),需對(duì)該壞點(diǎn)進(jìn)行處理。需要對(duì)數(shù)據(jù)壞點(diǎn)進(jìn)行插值替換,創(chuàng)建一個(gè)更平滑的數(shù)據(jù)序列作為參考[7],

    對(duì)壞點(diǎn)u(t+2)的具體處理過程如下。1)u(1)為從樣本點(diǎn)

    u(t)至u(t+4)共5個(gè)數(shù)據(jù)的中位數(shù)。2)u(2)為從樣本點(diǎn)u(1)t+1至u(1)t+3共3個(gè)數(shù)據(jù)的中位數(shù)。3)取u(2)t+1、u(2)t+2、

    u(2)t+2、u(2)t+3的算術(shù)平均值,即得u(3)。4)u(3)替代壞點(diǎn)u(t+2),完成該樣本壞點(diǎn)數(shù)據(jù)訂正。

    當(dāng)完成第一個(gè)矩形窗的數(shù)據(jù)質(zhì)量控制后,以1s為單位進(jìn)行滑移,繼續(xù)計(jì)算下一個(gè)矩形窗內(nèi)的數(shù)據(jù)質(zhì)量控制參數(shù),并完成壞點(diǎn)數(shù)據(jù)訂正任務(wù),并一直不斷重復(fù)上述過程,直到經(jīng)過4倍標(biāo)準(zhǔn)差過濾篩選后,沒有壞點(diǎn)被檢測(cè)出來為止。

    采用上述處理方法對(duì)臺(tái)風(fēng)“瑪莉亞”期間測(cè)風(fēng)塔10m和80m高度上,超聲風(fēng)速儀在三維方向ux、uy、uz記錄的0.1s時(shí)距風(fēng)速序列進(jìn)行判定、插值,得到處理后的風(fēng)場(chǎng)數(shù)據(jù)如圖2所示。

    由圖2可知,經(jīng)過質(zhì)量控制后,大量不合理的壞點(diǎn)被剔除,在圖上反映為大量失真的數(shù)據(jù)“毛刺”消失,大量不合理的數(shù)據(jù)被剔除,例如10m高度處z方向未經(jīng)處理的數(shù)據(jù)風(fēng)速達(dá)到100m/s,顯然為不合理數(shù)據(jù),經(jīng)過上述數(shù)據(jù)質(zhì)量控制處理后,風(fēng)速區(qū)間位于10m/s之內(nèi),提高了數(shù)據(jù)的合理性和可靠性。

    4 數(shù)據(jù)對(duì)比分析

    為了直觀分析數(shù)據(jù)修正前后風(fēng)速序列對(duì)平均風(fēng)速產(chǎn)生影響,對(duì)基礎(chǔ)風(fēng)速樣本進(jìn)行處理。由于Windmaster Pro三維超聲風(fēng)速儀同步輸出的是3個(gè)正交方向上的風(fēng)速分量ux、uy、uz,分別對(duì)應(yīng)正北向(x軸)、正西向(y軸)和豎直方向(z軸),儀器探頭夾角、儀器自身坐標(biāo)系以及地理坐標(biāo)系之間均存在轉(zhuǎn)換關(guān)系,因此有必要通過坐標(biāo)轉(zhuǎn)換進(jìn)行基礎(chǔ)風(fēng)速樣本的處理。對(duì)10min水平平均風(fēng)速U、水平風(fēng)向角?和豎向平均風(fēng)速W(豎直方向與儀器坐標(biāo)一致)[8-9]進(jìn)行計(jì)算,如公式(6)~公式(8)所示。

    (6)

    式中:U為平均風(fēng)速;、分別為風(fēng)速儀在x、y方向上分量的10min平均值。

    (7)

    式中:?為風(fēng)向角。

    W= (8)

    式中:W為豎向平均風(fēng)速;為風(fēng)速儀在z方向上分量的10min平均值。

    三維平均風(fēng)速如公式(9)所示。

    (9)

    式中:N為樣本數(shù),N=6000。

    將儀器坐標(biāo)旋轉(zhuǎn)?角,得到自然坐標(biāo),x、y、z軸分別代表主風(fēng)、側(cè)風(fēng)和垂直風(fēng)向,則u(t)在x軸上的投影u'(t)為縱風(fēng)向脈動(dòng)風(fēng)速,v(t)在y軸上的投影v'(t)為橫風(fēng)向脈動(dòng)風(fēng)速,w'(t)豎向脈動(dòng)風(fēng)速方向與儀器坐標(biāo)相同,如公式(10)~公式(12)所示。

    u'(t)=u(t)cos?+v(t)sin?-U (10)

    v'(t)=-u(t)sin?+v(t)cos? (11)

    w'(t)=w(t)-W (12)

    式中:u'(t)、v'(t)、w'(t)分別為3個(gè)方向脈動(dòng)風(fēng)速分量;u(t)、v(t)、w(t)分別為3個(gè)方向?qū)崪y(cè)脈動(dòng)風(fēng)速分量。

    完成上述計(jì)算即可得到3個(gè)方向上的脈動(dòng)風(fēng)速,為脈動(dòng)風(fēng)特性的數(shù)據(jù)誤差分析奠定基礎(chǔ)。以臺(tái)風(fēng)“瑪莉亞”期間觀測(cè)站10m高度層數(shù)據(jù)處理為例,對(duì)現(xiàn)場(chǎng)實(shí)測(cè)的風(fēng)場(chǎng)數(shù)據(jù)進(jìn)行處理比較,10m、80m高度處質(zhì)量控制前后平均風(fēng)速對(duì)比結(jié)果如圖3所示。

    由圖3可知,在10m及80m高度處,均在第70~80個(gè)樣本質(zhì)量控制前后誤差較大。其主要原因?yàn)?0~80樣本時(shí)間為臺(tái)風(fēng)登陸期間,期間風(fēng)速較大,氣候較復(fù)雜,引起的數(shù)據(jù)誤差較大,失真數(shù)據(jù)較多,從而導(dǎo)致質(zhì)量控制前后誤差較大。經(jīng)過數(shù)據(jù)質(zhì)量控制處理,有效地對(duì)不合理點(diǎn)進(jìn)行剔除和插補(bǔ),能夠有效減少復(fù)雜環(huán)境帶來的誤差,為后續(xù)脈動(dòng)風(fēng)場(chǎng)特性的分析以及實(shí)現(xiàn)臺(tái)風(fēng)氣候下風(fēng)障等構(gòu)造風(fēng)場(chǎng)有限元精細(xì)化模擬奠定基礎(chǔ)。

    5 結(jié)語

    臺(tái)風(fēng)風(fēng)場(chǎng)觀測(cè)通常在極端氣候下進(jìn)行,因此極容易受到復(fù)雜環(huán)境(例如雨水、儀器振動(dòng)等條件)的影響,同時(shí)三維風(fēng)速儀采樣頻率的選取也影響風(fēng)場(chǎng)數(shù)據(jù)的準(zhǔn)確性。本文基于超聲風(fēng)速儀無效風(fēng)速數(shù)據(jù)自動(dòng)識(shí)別提出了一種基于多倍截?cái)喾讲罘ǖ拿}動(dòng)風(fēng)速實(shí)測(cè)數(shù)據(jù)預(yù)處理方法,并利用臺(tái)風(fēng)“瑪莉亞”實(shí)測(cè)脈動(dòng)風(fēng)速數(shù)據(jù)驗(yàn)證了此方法的有效性。結(jié)果表明,經(jīng)儀器自動(dòng)判別生成的風(fēng)場(chǎng)數(shù)據(jù),采用四倍截?cái)喾軌蛴行У淖R(shí)別數(shù)據(jù)中壞點(diǎn),從而采取插值法進(jìn)行補(bǔ)充,得到完整可靠的風(fēng)場(chǎng)數(shù)據(jù)。通過對(duì)比質(zhì)量控制前后的實(shí)測(cè)數(shù)據(jù)、平均風(fēng)數(shù)據(jù)及脈動(dòng)風(fēng)數(shù)據(jù),驗(yàn)證了方法的有效性,可以為后續(xù)風(fēng)場(chǎng)參數(shù)分析,特別是海邊柔性橋梁結(jié)構(gòu)以及風(fēng)障等構(gòu)造研發(fā)中風(fēng)場(chǎng)精細(xì)化模擬提供可靠的數(shù)據(jù)基礎(chǔ)。

    參考文獻(xiàn)

    [1]林曉波,侯海濤,練江峰,等.雙幅式梁橋擋風(fēng)障阻風(fēng)性能影響參數(shù)分析[J].福州大學(xué)學(xué)報(bào)(自然科學(xué)版),2021,49(1):87-94.

    [2]LIN L,CHEN K,XIA D,et al.Analysis on the Wind Characteristics

    under Typhoon Climate at the Southeast Coast of China[J].Journal of Wind Engineering and Industrial Aerodynamics,2018(182):37-48.

    [3]趙林,潘晶晶,梁旭東,等.臺(tái)風(fēng)邊緣/中心區(qū)域經(jīng)歷平坦地貌時(shí)平均風(fēng)剖面特性[J].土木工程學(xué)報(bào),2016,49(8):45-50.

    [4]宋麗莉,陳雯超,黃浩輝.工程抗臺(tái)風(fēng)研究中風(fēng)觀測(cè)數(shù)據(jù)的可靠性和代表性判別[J].氣象科技進(jìn)展,2011,1(1):35-39.

    [5]卞林根,陸龍驊,程彥杰,等.青藏高原南部昌都地區(qū)近地層湍流輸送的觀測(cè)研究[J].應(yīng)用氣象學(xué)報(bào),2001,12(1):1-13.

    [6]陳紅巖,胡非,曾慶存.處理時(shí)間序列提高計(jì)算湍流通量的精度[J].氣候與環(huán)境研究,2000,5(3):304-311.

    [7]CASTELAO G P. A Flexible System for Automatic Quality Control

    of Oceanographic Data[J].Physics, 2016(11):1-15.

    [8]XU Y L, ZHAN S. Field Measurements of Di Wang Tower during

    Typhoon York[J].Journal of Wind Engineering and Industrial Aerodynamics,

    2001,89(1): 73-93.

    [9]李家亮.基于近地觀測(cè)的登陸臺(tái)風(fēng)湍流積分尺度研究[D].哈爾濱:哈爾濱工業(yè)大學(xué),2012.

    猜你喜歡
    插值法數(shù)據(jù)處理臺(tái)風(fēng)
    臺(tái)風(fēng)過韓
    認(rèn)知診斷缺失數(shù)據(jù)處理方法的比較:零替換、多重插補(bǔ)與極大似然估計(jì)法*
    ILWT-EEMD數(shù)據(jù)處理的ELM滾動(dòng)軸承故障診斷
    臺(tái)風(fēng)來了
    小讀者(2020年4期)2020-06-16 03:33:46
    《計(jì)算方法》關(guān)于插值法的教學(xué)方法研討
    臺(tái)風(fēng)愛搗亂
    臺(tái)風(fēng)來時(shí)怎樣應(yīng)對(duì)
    基于希爾伯特- 黃變換的去噪法在外測(cè)數(shù)據(jù)處理中的應(yīng)用
    基于二次插值法的布谷鳥搜索算法研究
    Newton插值法在光伏發(fā)電最大功率跟蹤中的應(yīng)用
    丰满乱子伦码专区| 精品人妻一区二区三区麻豆 | 又爽又黄无遮挡网站| 免费看日本二区| 国产成人欧美在线观看| 亚洲男人的天堂狠狠| 尤物成人国产欧美一区二区三区| 色尼玛亚洲综合影院| 久久久久精品国产欧美久久久| 久久国产精品人妻蜜桃| 十八禁国产超污无遮挡网站| 国内精品久久久久精免费| 丁香六月欧美| 欧美日韩国产亚洲二区| 在线观看av片永久免费下载| 综合色av麻豆| 国产成人欧美在线观看| 少妇高潮的动态图| 国产高清激情床上av| а√天堂www在线а√下载| 精品无人区乱码1区二区| 国产探花在线观看一区二区| 一区二区三区高清视频在线| 亚洲18禁久久av| 午夜老司机福利剧场| 欧美日韩亚洲国产一区二区在线观看| 成年版毛片免费区| 我的女老师完整版在线观看| 在线观看午夜福利视频| 午夜精品一区二区三区免费看| 久久香蕉精品热| 黄色女人牲交| 亚洲黑人精品在线| 国产av一区在线观看免费| 日韩精品中文字幕看吧| 热99在线观看视频| 亚洲五月天丁香| 亚洲av电影不卡..在线观看| 亚洲av二区三区四区| 婷婷色综合大香蕉| 亚洲五月天丁香| 日本在线视频免费播放| 永久网站在线| 久久精品综合一区二区三区| 两人在一起打扑克的视频| 免费av观看视频| 又粗又爽又猛毛片免费看| 国产精品三级大全| 亚洲人与动物交配视频| 搡老熟女国产l中国老女人| 午夜免费成人在线视频| 岛国在线免费视频观看| 久久精品国产亚洲av天美| 国产白丝娇喘喷水9色精品| 午夜福利欧美成人| 日韩高清综合在线| 久久久久精品国产欧美久久久| 麻豆成人av在线观看| 一区福利在线观看| 男人和女人高潮做爰伦理| 日本一本二区三区精品| 全区人妻精品视频| 在线观看一区二区三区| 动漫黄色视频在线观看| 精品人妻视频免费看| 国产精品一区二区性色av| 男人舔女人下体高潮全视频| 五月玫瑰六月丁香| 一夜夜www| 在线观看美女被高潮喷水网站 | 国产亚洲精品久久久com| 麻豆成人午夜福利视频| 国产v大片淫在线免费观看| 夜夜躁狠狠躁天天躁| 日日摸夜夜添夜夜添av毛片 | 国产精品综合久久久久久久免费| 亚洲美女黄片视频| 深夜精品福利| 久久久久久国产a免费观看| 老鸭窝网址在线观看| 首页视频小说图片口味搜索| 国产熟女xx| 国产精品久久久久久久电影| 哪里可以看免费的av片| 最新在线观看一区二区三区| 国产精品国产高清国产av| 国产伦精品一区二区三区视频9| 啪啪无遮挡十八禁网站| 淫妇啪啪啪对白视频| 免费在线观看影片大全网站| 亚洲性夜色夜夜综合| aaaaa片日本免费| 搡老熟女国产l中国老女人| 少妇人妻一区二区三区视频| 久久午夜福利片| 最近最新中文字幕大全电影3| 一区二区三区激情视频| 一个人免费在线观看的高清视频| 欧美成人免费av一区二区三区| 99久久九九国产精品国产免费| 亚洲国产精品sss在线观看| 亚洲精品成人久久久久久| 99精品在免费线老司机午夜| 亚洲午夜理论影院| 亚洲,欧美,日韩| 91av网一区二区| 亚洲七黄色美女视频| 亚洲男人的天堂狠狠| 99视频精品全部免费 在线| www.色视频.com| 久久久久久国产a免费观看| 久久久色成人| 国产 一区 欧美 日韩| 又黄又爽又免费观看的视频| 精品午夜福利在线看| 久99久视频精品免费| 在线播放无遮挡| 丁香欧美五月| 毛片女人毛片| av专区在线播放| 国产精品影院久久| 久久午夜福利片| 午夜精品一区二区三区免费看| 国产爱豆传媒在线观看| 一区二区三区四区激情视频 | eeuss影院久久| 嫩草影院入口| 99精品在免费线老司机午夜| 五月伊人婷婷丁香| 精品一区二区三区视频在线| 搞女人的毛片| 美女高潮喷水抽搐中文字幕| 色视频www国产| 激情在线观看视频在线高清| 丰满人妻一区二区三区视频av| 好男人在线观看高清免费视频| 国产亚洲精品综合一区在线观看| 国产欧美日韩一区二区三| 成年人黄色毛片网站| 成人一区二区视频在线观看| 五月伊人婷婷丁香| 国产人妻一区二区三区在| 91九色精品人成在线观看| 午夜福利欧美成人| 18禁黄网站禁片午夜丰满| 国产伦一二天堂av在线观看| 亚洲美女视频黄频| 亚洲人与动物交配视频| 久久国产精品影院| 亚洲av成人av| 极品教师在线视频| 亚洲第一区二区三区不卡| 欧美性猛交黑人性爽| 特大巨黑吊av在线直播| 日韩欧美 国产精品| 国产一区二区亚洲精品在线观看| 免费黄网站久久成人精品 | 97超级碰碰碰精品色视频在线观看| 亚洲最大成人中文| 亚洲熟妇熟女久久| 精品午夜福利视频在线观看一区| 国产一区二区三区视频了| 国产精品自产拍在线观看55亚洲| 精品人妻一区二区三区麻豆 | 99久久99久久久精品蜜桃| 色噜噜av男人的天堂激情| 亚洲欧美清纯卡通| 国产免费男女视频| 国产91精品成人一区二区三区| 欧美日韩福利视频一区二区| 婷婷精品国产亚洲av| 国产成人福利小说| 亚洲美女搞黄在线观看 | 中文亚洲av片在线观看爽| 精品久久国产蜜桃| 精品人妻一区二区三区麻豆 | 国产伦在线观看视频一区| 日本 av在线| 久久久色成人| 一个人看的www免费观看视频| 久久精品影院6| 最新在线观看一区二区三区| 亚洲成人久久性| 欧美日韩乱码在线| 天堂av国产一区二区熟女人妻| 少妇人妻一区二区三区视频| 亚洲成人久久爱视频| 国产毛片a区久久久久| 国产精品一区二区三区四区久久| 人妻久久中文字幕网| 久久精品91蜜桃| 一二三四社区在线视频社区8| 国产日本99.免费观看| 窝窝影院91人妻| 久久热精品热| 色哟哟·www| 久久久久久久精品吃奶| 3wmmmm亚洲av在线观看| 色综合站精品国产| 精品熟女少妇八av免费久了| 九色国产91popny在线| 国产一区二区三区视频了| 欧美黑人欧美精品刺激| 日本 欧美在线| 国产精品伦人一区二区| 久久久久久久久中文| 日韩免费av在线播放| 天堂√8在线中文| 亚洲五月婷婷丁香| 国产成人欧美在线观看| 国产精品亚洲美女久久久| 神马国产精品三级电影在线观看| 欧美精品国产亚洲| 午夜老司机福利剧场| 又爽又黄a免费视频| 国产白丝娇喘喷水9色精品| 少妇被粗大猛烈的视频| 搡老妇女老女人老熟妇| 免费看日本二区| 老司机深夜福利视频在线观看| 无遮挡黄片免费观看| 观看美女的网站| 欧美中文日本在线观看视频| 亚洲va日本ⅴa欧美va伊人久久| 久久国产乱子伦精品免费另类| 国产精品女同一区二区软件 | 国产亚洲精品久久久com| 精品国内亚洲2022精品成人| 日韩欧美一区二区三区在线观看| 又爽又黄无遮挡网站| 精品一区二区三区人妻视频| 久久午夜福利片| 91九色精品人成在线观看| 久久这里只有精品中国| 18美女黄网站色大片免费观看| 亚洲天堂国产精品一区在线| www.熟女人妻精品国产| 精品久久久久久成人av| 搞女人的毛片| 最好的美女福利视频网| 国产精品自产拍在线观看55亚洲| 亚洲欧美激情综合另类| 日韩欧美在线乱码| av在线观看视频网站免费| 丰满乱子伦码专区| 国产亚洲欧美98| 国产伦一二天堂av在线观看| 国产一区二区在线观看日韩| 国产精品久久久久久人妻精品电影| 成熟少妇高潮喷水视频| 欧美不卡视频在线免费观看| eeuss影院久久| 一区二区三区免费毛片| 欧美绝顶高潮抽搐喷水| 淫秽高清视频在线观看| 成人美女网站在线观看视频| 午夜a级毛片| 少妇人妻精品综合一区二区 | 亚洲18禁久久av| 婷婷丁香在线五月| 亚洲精品一区av在线观看| 男人和女人高潮做爰伦理| 免费无遮挡裸体视频| 国产精品精品国产色婷婷| 九色国产91popny在线| eeuss影院久久| 一级黄片播放器| 亚洲成人免费电影在线观看| 搡老妇女老女人老熟妇| 欧美乱色亚洲激情| 嫩草影院精品99| 精品久久久久久久久久免费视频| www.色视频.com| 免费看日本二区| 成人国产一区最新在线观看| 成人性生交大片免费视频hd| 日本黄大片高清| 性插视频无遮挡在线免费观看| 国产熟女xx| 变态另类丝袜制服| 免费观看精品视频网站| 波多野结衣高清作品| 高清在线国产一区| 国产精品av视频在线免费观看| 久久久久性生活片| 好男人电影高清在线观看| 日韩欧美精品免费久久 | 日韩亚洲欧美综合| 亚洲中文字幕日韩| 十八禁人妻一区二区| 黄色配什么色好看| 午夜亚洲福利在线播放| 欧美+日韩+精品| 久久久久久国产a免费观看| 九色国产91popny在线| 久久亚洲精品不卡| 婷婷丁香在线五月| 日韩中字成人| 色综合欧美亚洲国产小说| 日韩精品青青久久久久久| 国产精品嫩草影院av在线观看 | 亚洲专区中文字幕在线| 日韩有码中文字幕| 91久久精品电影网| 1024手机看黄色片| 色哟哟哟哟哟哟| 一本久久中文字幕| 三级男女做爰猛烈吃奶摸视频| 国产黄色小视频在线观看| 91在线观看av| 亚洲狠狠婷婷综合久久图片| 色综合婷婷激情| 久久精品影院6| 男女那种视频在线观看| 色吧在线观看| 国产高清三级在线| 男女床上黄色一级片免费看| av在线老鸭窝| 日韩欧美三级三区| 久久久久久久亚洲中文字幕 | 亚洲av中文字字幕乱码综合| 国产一区二区亚洲精品在线观看| 真人做人爱边吃奶动态| 亚洲18禁久久av| 观看免费一级毛片| 免费电影在线观看免费观看| 国产aⅴ精品一区二区三区波| 国产伦精品一区二区三区视频9| 特大巨黑吊av在线直播| avwww免费| 国产单亲对白刺激| 少妇人妻一区二区三区视频| 91麻豆av在线| 五月伊人婷婷丁香| 亚洲久久久久久中文字幕| 久久国产乱子免费精品| 婷婷丁香在线五月| 国内揄拍国产精品人妻在线| 日本在线视频免费播放| 高清日韩中文字幕在线| 亚洲av成人精品一区久久| 久久中文看片网| 在线播放国产精品三级| 国产私拍福利视频在线观看| 99热这里只有是精品50| 天堂动漫精品| 色哟哟·www| 中文字幕人妻熟人妻熟丝袜美| 一区福利在线观看| 精品久久久久久久久久久久久| 久久亚洲精品不卡| 尤物成人国产欧美一区二区三区| 久久中文看片网| 国产又黄又爽又无遮挡在线| 亚洲18禁久久av| 淫妇啪啪啪对白视频| 黄片小视频在线播放| 国产大屁股一区二区在线视频| 婷婷六月久久综合丁香| 中文字幕人妻熟人妻熟丝袜美| 国产精品久久久久久精品电影| 直男gayav资源| 日韩欧美一区二区三区在线观看| 中文字幕人妻熟人妻熟丝袜美| 久久久久久久久大av| 久久精品国产99精品国产亚洲性色| 日本五十路高清| 97超级碰碰碰精品色视频在线观看| 日本五十路高清| 桃色一区二区三区在线观看| 久久国产乱子免费精品| 精品久久久久久成人av| 九九久久精品国产亚洲av麻豆| 亚洲美女搞黄在线观看 | 无遮挡黄片免费观看| 国产精品亚洲av一区麻豆| 精品人妻视频免费看| 成人鲁丝片一二三区免费| 日韩成人在线观看一区二区三区| 校园春色视频在线观看| 亚洲av.av天堂| 一本精品99久久精品77| 国产成年人精品一区二区| 日本一本二区三区精品| 免费高清视频大片| 免费在线观看日本一区| 午夜福利免费观看在线| 国产熟女xx| 国产伦人伦偷精品视频| 亚洲av五月六月丁香网| 中亚洲国语对白在线视频| 亚洲经典国产精华液单 | 国产三级中文精品| 成年免费大片在线观看| 中文字幕人成人乱码亚洲影| 国产精品久久电影中文字幕| 欧美午夜高清在线| 精品一区二区三区av网在线观看| 国产欧美日韩精品亚洲av| 国产欧美日韩精品一区二区| 听说在线观看完整版免费高清| 亚洲成av人片在线播放无| 国产精品女同一区二区软件 | 国内揄拍国产精品人妻在线| 国产大屁股一区二区在线视频| 亚洲自偷自拍三级| 免费一级毛片在线播放高清视频| 偷拍熟女少妇极品色| 韩国av一区二区三区四区| 亚洲熟妇中文字幕五十中出| 日本在线视频免费播放| 免费看a级黄色片| 97热精品久久久久久| 一本综合久久免费| 免费高清视频大片| av黄色大香蕉| eeuss影院久久| 桃色一区二区三区在线观看| 两个人的视频大全免费| 国模一区二区三区四区视频| 内射极品少妇av片p| 变态另类丝袜制服| 乱码一卡2卡4卡精品| 久久精品国产自在天天线| 99热这里只有是精品在线观看 | 国产美女午夜福利| 国产极品精品免费视频能看的| 91久久精品国产一区二区成人| 中文资源天堂在线| 国产黄a三级三级三级人| 级片在线观看| avwww免费| 国产单亲对白刺激| 亚洲精品亚洲一区二区| 精品一区二区三区视频在线| 色哟哟哟哟哟哟| 无人区码免费观看不卡| 亚洲国产精品999在线| 校园春色视频在线观看| 亚州av有码| 麻豆一二三区av精品| 免费高清视频大片| 极品教师在线视频| 精品人妻一区二区三区麻豆 | 日本免费a在线| 99精品在免费线老司机午夜| 国产av不卡久久| 国产av一区在线观看免费| 亚洲成人精品中文字幕电影| 欧美成狂野欧美在线观看| 久久九九热精品免费| 永久网站在线| aaaaa片日本免费| 91av网一区二区| 国产av一区在线观看免费| 十八禁网站免费在线| 成人精品一区二区免费| 精品免费久久久久久久清纯| 午夜两性在线视频| 中文字幕精品亚洲无线码一区| 美女免费视频网站| 最近最新中文字幕大全电影3| 精品久久国产蜜桃| 亚洲最大成人中文| 毛片女人毛片| 亚洲精品一区av在线观看| 午夜福利成人在线免费观看| 丰满人妻一区二区三区视频av| 日日摸夜夜添夜夜添小说| 免费观看人在逋| 一级黄色大片毛片| 精品欧美国产一区二区三| 精品人妻1区二区| 人妻久久中文字幕网| 久久中文看片网| 小蜜桃在线观看免费完整版高清| 91在线精品国自产拍蜜月| 九九在线视频观看精品| 非洲黑人性xxxx精品又粗又长| 精品人妻视频免费看| 精品久久久久久久久久久久久| 精品久久国产蜜桃| 又黄又爽又免费观看的视频| 一本一本综合久久| 日本与韩国留学比较| 国产精品久久电影中文字幕| 99在线人妻在线中文字幕| 久久久久久久久久成人| 午夜免费激情av| 久久欧美精品欧美久久欧美| 91麻豆av在线| 欧美精品国产亚洲| 日韩国内少妇激情av| 亚洲五月婷婷丁香| 最好的美女福利视频网| 欧美国产日韩亚洲一区| 琪琪午夜伦伦电影理论片6080| 亚洲最大成人中文| 欧美色视频一区免费| 精品久久久久久,| 色尼玛亚洲综合影院| 女生性感内裤真人,穿戴方法视频| 日本熟妇午夜| 亚洲欧美日韩高清专用| 精品久久久久久,| 最近最新中文字幕大全电影3| or卡值多少钱| 亚洲精品一卡2卡三卡4卡5卡| 久99久视频精品免费| 国产一区二区在线观看日韩| 99热精品在线国产| 岛国在线免费视频观看| 悠悠久久av| 亚洲一区高清亚洲精品| 美女高潮喷水抽搐中文字幕| 丁香六月欧美| 国产精品日韩av在线免费观看| 国产成人啪精品午夜网站| 中文资源天堂在线| 听说在线观看完整版免费高清| 色5月婷婷丁香| 在线观看美女被高潮喷水网站 | 在线天堂最新版资源| 天天躁日日操中文字幕| 国产毛片a区久久久久| 色播亚洲综合网| 女同久久另类99精品国产91| 深夜a级毛片| 国产单亲对白刺激| 久久精品国产自在天天线| 亚洲不卡免费看| 亚洲综合色惰| 99热这里只有是精品50| 亚洲内射少妇av| 国产在线男女| 久久久精品欧美日韩精品| 国产亚洲精品av在线| 深爱激情五月婷婷| 久久精品国产亚洲av涩爱 | 日韩有码中文字幕| 97碰自拍视频| 国产日本99.免费观看| 亚洲美女黄片视频| 极品教师在线视频| 最新中文字幕久久久久| 欧美国产日韩亚洲一区| 我的女老师完整版在线观看| 国产精品一区二区免费欧美| 国产国拍精品亚洲av在线观看| 美女被艹到高潮喷水动态| 亚洲精品粉嫩美女一区| 在线国产一区二区在线| 国产精品自产拍在线观看55亚洲| 欧美乱色亚洲激情| 99久久无色码亚洲精品果冻| 在线观看66精品国产| 久久久久九九精品影院| 我的老师免费观看完整版| 国产人妻一区二区三区在| 午夜福利在线在线| 中出人妻视频一区二区| 日本 av在线| 中文字幕人妻熟人妻熟丝袜美| 国产精品一区二区三区四区久久| 99久久精品一区二区三区| 欧美性感艳星| 在线十欧美十亚洲十日本专区| 嫩草影院入口| 午夜福利在线观看免费完整高清在 | 亚洲综合色惰| 好看av亚洲va欧美ⅴa在| 欧美乱妇无乱码| 一a级毛片在线观看| 午夜免费男女啪啪视频观看 | 国产高清有码在线观看视频| 成人av一区二区三区在线看| 精品久久久久久成人av| 久久久精品欧美日韩精品| 欧美三级亚洲精品| 校园春色视频在线观看| 亚洲最大成人av| 亚洲精品在线美女| 校园春色视频在线观看| 亚洲最大成人av| 精品一区二区三区视频在线观看免费| 亚洲av不卡在线观看| 九九热线精品视视频播放| 亚洲国产高清在线一区二区三| 极品教师在线视频| 午夜福利免费观看在线| 成年女人看的毛片在线观看| 热99re8久久精品国产| 97人妻精品一区二区三区麻豆| 桃红色精品国产亚洲av| 高清在线国产一区| 国产91精品成人一区二区三区| 日本一本二区三区精品| 色哟哟·www| 欧美日韩国产亚洲二区| 啪啪无遮挡十八禁网站| 国产一级毛片七仙女欲春2| 欧美xxxx性猛交bbbb| 无遮挡黄片免费观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲综合色惰| 窝窝影院91人妻| 亚洲七黄色美女视频| 亚洲专区国产一区二区| 国产欧美日韩一区二区精品| 亚洲成人精品中文字幕电影| 级片在线观看| 欧美一区二区国产精品久久精品| 黄色配什么色好看| 一本一本综合久久| 深夜a级毛片|