• <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毛片| a级一级毛片免费在线观看| 亚洲国产精品久久男人天堂| 国产单亲对白刺激| 性欧美人与动物交配| 久久99热6这里只有精品| 亚洲成人久久爱视频| x7x7x7水蜜桃| 成人av在线播放网站| 91麻豆精品激情在线观看国产| 久久精品国产清高在天天线| 精品午夜福利在线看| 男女下面进入的视频免费午夜| 国产一区二区三区av在线 | 毛片女人毛片| 欧美色视频一区免费| 亚洲欧美日韩高清专用| 国产人妻一区二区三区在| 热99re8久久精品国产| 美女 人体艺术 gogo| 国产单亲对白刺激| 联通29元200g的流量卡| 亚洲av中文字字幕乱码综合| 一个人看视频在线观看www免费| 欧美精品啪啪一区二区三区| 免费无遮挡裸体视频| 波野结衣二区三区在线| 午夜日韩欧美国产| 色综合婷婷激情| 国产精品亚洲美女久久久| 伊人久久精品亚洲午夜| 日韩人妻高清精品专区| 97人妻精品一区二区三区麻豆| 少妇丰满av| 五月玫瑰六月丁香| 69人妻影院| 天天躁日日操中文字幕| 久久久久久久精品吃奶| 久99久视频精品免费| 校园春色视频在线观看| 99热精品在线国产| 一个人看的www免费观看视频| 在线天堂最新版资源| 人妻夜夜爽99麻豆av| 人人妻,人人澡人人爽秒播| 少妇被粗大猛烈的视频| 久久久午夜欧美精品| 国产精品一区二区性色av| 男女啪啪激烈高潮av片| 黄片wwwwww| 久久中文看片网| a级毛片免费高清观看在线播放| 在线观看美女被高潮喷水网站| 亚洲五月天丁香| 久久久国产成人免费| 国产日本99.免费观看| 婷婷丁香在线五月| 精品国内亚洲2022精品成人| 日本熟妇午夜| 国产高潮美女av| 午夜影院日韩av| 九九久久精品国产亚洲av麻豆| 少妇的逼好多水| 国内揄拍国产精品人妻在线| 精品久久国产蜜桃| 日日干狠狠操夜夜爽| 国产精品一区二区三区四区久久| 国产精品三级大全| av视频在线观看入口| 亚洲精品久久国产高清桃花| 又爽又黄无遮挡网站| 亚洲美女搞黄在线观看 | 亚洲欧美清纯卡通| 精品一区二区三区视频在线| 在线免费观看不下载黄p国产 | 色哟哟哟哟哟哟| 国产老妇女一区| 日本免费一区二区三区高清不卡| 精品一区二区三区人妻视频| 国产免费一级a男人的天堂| 国产女主播在线喷水免费视频网站 | 国产高潮美女av| 久久欧美精品欧美久久欧美| 亚洲一级一片aⅴ在线观看| 欧美一区二区精品小视频在线| 亚洲成人中文字幕在线播放| 亚洲美女视频黄频| 国产精品久久久久久精品电影| 色噜噜av男人的天堂激情| 美女高潮喷水抽搐中文字幕| 色哟哟·www| 国产精品日韩av在线免费观看| 最近视频中文字幕2019在线8| 观看美女的网站| 久久久久久久精品吃奶| 麻豆成人午夜福利视频| 亚洲精品乱码久久久v下载方式| 国产在线精品亚洲第一网站| 久久久久久久久久成人| 亚洲成人免费电影在线观看| 亚洲 国产 在线| www日本黄色视频网| 精品一区二区免费观看| 久久精品影院6| 免费看美女性在线毛片视频| 亚洲精华国产精华精| 最近最新免费中文字幕在线| 久99久视频精品免费| 亚洲久久久久久中文字幕| 国产精品野战在线观看| 熟女电影av网| 18禁黄网站禁片免费观看直播| 国产精品国产三级国产av玫瑰| 欧美性猛交╳xxx乱大交人| 成人午夜高清在线视频| 99久久精品国产国产毛片| 在线看三级毛片| bbb黄色大片| 精品久久久久久久末码| 午夜激情福利司机影院| 一本精品99久久精品77| 久久久久免费精品人妻一区二区| 午夜福利高清视频| 欧美黑人欧美精品刺激| 一个人免费在线观看电影| 国产精品嫩草影院av在线观看 | 别揉我奶头 嗯啊视频| 在线观看午夜福利视频| 欧美另类亚洲清纯唯美| 午夜日韩欧美国产| 麻豆久久精品国产亚洲av| eeuss影院久久| 99在线人妻在线中文字幕| 国国产精品蜜臀av免费| 真人一进一出gif抽搐免费| 免费看av在线观看网站| 亚洲成av人片在线播放无| 12—13女人毛片做爰片一| eeuss影院久久| 观看美女的网站| 亚洲精品456在线播放app | 国产一区二区在线av高清观看| 性色avwww在线观看| 真实男女啪啪啪动态图| 亚洲电影在线观看av| 亚州av有码| 久久久色成人| 中文资源天堂在线| 男人舔女人下体高潮全视频| 欧美激情在线99| 欧洲精品卡2卡3卡4卡5卡区| 亚洲中文字幕一区二区三区有码在线看| 亚洲av成人av| 99热6这里只有精品| 免费一级毛片在线播放高清视频| 成年版毛片免费区| 99久久精品热视频| 桃红色精品国产亚洲av| 高清在线国产一区| 三级男女做爰猛烈吃奶摸视频| 成人三级黄色视频| 中文字幕熟女人妻在线| 亚洲七黄色美女视频| 成人av在线播放网站| 级片在线观看| 亚洲精华国产精华液的使用体验 | 国产91精品成人一区二区三区| 乱人视频在线观看| 国产成人影院久久av| 在线免费十八禁| 内地一区二区视频在线| 赤兔流量卡办理| 老熟妇乱子伦视频在线观看| 美女免费视频网站| 国内毛片毛片毛片毛片毛片| 日本黄色片子视频| 91狼人影院| 999久久久精品免费观看国产| avwww免费| 超碰av人人做人人爽久久| 午夜精品久久久久久毛片777| 国内精品久久久久久久电影| 国产探花在线观看一区二区| 一区二区三区高清视频在线| 91麻豆精品激情在线观看国产| 在线观看美女被高潮喷水网站| 日本 欧美在线| 精品久久久久久久末码| 一卡2卡三卡四卡精品乱码亚洲| 日日摸夜夜添夜夜添小说| 欧美日韩瑟瑟在线播放| 亚洲精品成人久久久久久| 中文字幕熟女人妻在线| 午夜激情欧美在线| 黄色欧美视频在线观看| 亚洲成人免费电影在线观看| 免费av毛片视频| 最近中文字幕高清免费大全6 | 午夜日韩欧美国产| 欧美日韩黄片免| www日本黄色视频网| 乱人视频在线观看| 美女高潮的动态| 成人精品一区二区免费| 男女啪啪激烈高潮av片| 亚洲最大成人av| 久久久色成人| 男女那种视频在线观看| 精品一区二区三区人妻视频| 成年女人永久免费观看视频| 在线看三级毛片| 日日摸夜夜添夜夜添av毛片 | 听说在线观看完整版免费高清| 非洲黑人性xxxx精品又粗又长| 亚洲,欧美,日韩| 一a级毛片在线观看| 最后的刺客免费高清国语| 国产大屁股一区二区在线视频| 亚洲av二区三区四区| 免费不卡的大黄色大毛片视频在线观看 | 给我免费播放毛片高清在线观看| 日本爱情动作片www.在线观看 | 亚洲国产精品成人综合色| 久久精品国产99精品国产亚洲性色| 亚洲三级黄色毛片| 亚洲精品亚洲一区二区| 亚洲最大成人av| 成人国产麻豆网| 日日干狠狠操夜夜爽| 免费看美女性在线毛片视频| 亚洲欧美日韩高清专用| 日本一本二区三区精品| 久久香蕉精品热| 老熟妇仑乱视频hdxx| 99在线人妻在线中文字幕| 禁无遮挡网站| 久久久久精品国产欧美久久久| 亚洲中文字幕一区二区三区有码在线看| 国产91精品成人一区二区三区| 嫩草影院入口| 久久久久精品国产欧美久久久| 人妻制服诱惑在线中文字幕| 日韩欧美精品v在线| 极品教师在线免费播放| 欧美性感艳星| 在线看三级毛片| 成人一区二区视频在线观看| 欧美又色又爽又黄视频| 琪琪午夜伦伦电影理论片6080| 国产精品美女特级片免费视频播放器| 成人欧美大片| 12—13女人毛片做爰片一| 女生性感内裤真人,穿戴方法视频| 91在线精品国自产拍蜜月| 在线免费观看不下载黄p国产 | 国产精品亚洲美女久久久| 欧美三级亚洲精品| 亚洲av成人av| 国产精品美女特级片免费视频播放器| 久久精品国产99精品国产亚洲性色| 午夜久久久久精精品| 久久99热6这里只有精品| 一卡2卡三卡四卡精品乱码亚洲| 两性午夜刺激爽爽歪歪视频在线观看| 日本五十路高清| 亚洲性夜色夜夜综合| 国产麻豆成人av免费视频| 欧美性猛交黑人性爽| 亚洲精品456在线播放app | videossex国产| 国产aⅴ精品一区二区三区波| 成人国产一区最新在线观看| 一边摸一边抽搐一进一小说| netflix在线观看网站| 精品一区二区免费观看| 国产69精品久久久久777片| 麻豆成人av在线观看| 欧美日韩中文字幕国产精品一区二区三区| 精品久久久久久成人av| 在线观看av片永久免费下载| 午夜免费激情av| 日韩欧美在线乱码| 国产免费av片在线观看野外av| 乱系列少妇在线播放| 国产乱人伦免费视频| 成人特级av手机在线观看| 亚洲图色成人| 村上凉子中文字幕在线| 午夜福利视频1000在线观看| 国内精品一区二区在线观看| 长腿黑丝高跟| 乱码一卡2卡4卡精品| 日韩欧美国产一区二区入口| 日韩亚洲欧美综合| 小说图片视频综合网站| 校园人妻丝袜中文字幕| 免费一级毛片在线播放高清视频| 99久久中文字幕三级久久日本| 两人在一起打扑克的视频| 亚洲人成网站在线播放欧美日韩| 国产伦人伦偷精品视频| 精品久久久久久久人妻蜜臀av| 99热网站在线观看| 成年女人毛片免费观看观看9| 国产极品精品免费视频能看的| 最近中文字幕高清免费大全6 | 亚洲av.av天堂| 一级黄片播放器| 熟女电影av网| 啪啪无遮挡十八禁网站| 成人永久免费在线观看视频| 亚州av有码| 午夜福利成人在线免费观看| 午夜免费成人在线视频| 国产欧美日韩一区二区精品| 一级黄色大片毛片| 亚洲七黄色美女视频| 看片在线看免费视频| 国产伦精品一区二区三区视频9| 热99re8久久精品国产| 中文亚洲av片在线观看爽| 午夜亚洲福利在线播放| 天堂av国产一区二区熟女人妻| 国产私拍福利视频在线观看| 免费黄网站久久成人精品| 国产三级中文精品| 国产欧美日韩一区二区精品| 亚洲专区国产一区二区| 狠狠狠狠99中文字幕| 99久久精品一区二区三区| 久久99热这里只有精品18| 桃色一区二区三区在线观看| 此物有八面人人有两片| 午夜免费激情av| 看免费成人av毛片| 91久久精品国产一区二区三区| 国产在线精品亚洲第一网站| 亚洲美女视频黄频| 欧美极品一区二区三区四区| 日韩一本色道免费dvd| 亚洲五月天丁香| 两性午夜刺激爽爽歪歪视频在线观看| 九色成人免费人妻av| 精品人妻一区二区三区麻豆 | 久久中文看片网| 欧美极品一区二区三区四区| 日本色播在线视频| 成人鲁丝片一二三区免费| 22中文网久久字幕| 亚洲av免费高清在线观看| 少妇高潮的动态图| 两个人视频免费观看高清| 色综合婷婷激情| 久久精品影院6| 中文字幕人妻熟人妻熟丝袜美| 国产真实伦视频高清在线观看 | 五月伊人婷婷丁香| 美女高潮的动态| 日韩人妻高清精品专区| 婷婷丁香在线五月| 精品久久久久久,| 一级黄片播放器| 国产免费一级a男人的天堂| 色精品久久人妻99蜜桃| 看免费成人av毛片| 亚洲欧美日韩卡通动漫| 看免费成人av毛片| 岛国在线免费视频观看| 搞女人的毛片| a级毛片a级免费在线| 精品久久久噜噜| www.色视频.com| 久久99热6这里只有精品| 国产一级毛片七仙女欲春2| 国产一区二区在线av高清观看| 日本三级黄在线观看| 午夜影院日韩av| 一进一出抽搐动态| 亚洲精品色激情综合| 日韩av在线大香蕉| 大又大粗又爽又黄少妇毛片口| 美女大奶头视频| 成年免费大片在线观看| 精品人妻1区二区| 别揉我奶头~嗯~啊~动态视频| 欧美区成人在线视频| 啪啪无遮挡十八禁网站| h日本视频在线播放| 岛国在线免费视频观看| 黄色一级大片看看| 岛国在线免费视频观看| 成人美女网站在线观看视频| 日韩欧美精品v在线| 18禁裸乳无遮挡免费网站照片| 色综合婷婷激情| 亚洲自偷自拍三级| 女的被弄到高潮叫床怎么办 | 亚洲精华国产精华精| 99热精品在线国产| 国产不卡一卡二| 久久午夜亚洲精品久久| 国产精品乱码一区二三区的特点| 十八禁网站免费在线| 欧美中文日本在线观看视频| 国产午夜福利久久久久久| 国产伦人伦偷精品视频| 日韩欧美国产在线观看| 精品人妻偷拍中文字幕| 亚洲av一区综合| 99热这里只有是精品50| 午夜影院日韩av| 免费观看的影片在线观看| 亚洲国产日韩欧美精品在线观看| 丰满的人妻完整版| 97超级碰碰碰精品色视频在线观看| av福利片在线观看| a级毛片免费高清观看在线播放| 欧洲精品卡2卡3卡4卡5卡区| 大型黄色视频在线免费观看| 夜夜爽天天搞| 午夜免费男女啪啪视频观看 | 免费高清视频大片| 男人狂女人下面高潮的视频| 日韩欧美 国产精品| 亚洲av熟女| a级一级毛片免费在线观看| 一a级毛片在线观看| 无人区码免费观看不卡| 国产女主播在线喷水免费视频网站 | 国产亚洲av嫩草精品影院| 国产白丝娇喘喷水9色精品| 搡老妇女老女人老熟妇| 欧美激情久久久久久爽电影| 亚洲国产精品成人综合色| 成人高潮视频无遮挡免费网站| or卡值多少钱| 亚洲中文日韩欧美视频| 亚洲性夜色夜夜综合| 免费人成在线观看视频色| 99久久中文字幕三级久久日本| 99热这里只有精品一区| 色在线成人网| 久久久国产成人精品二区| 91av网一区二区| 人妻久久中文字幕网| 欧美bdsm另类| 国产精品永久免费网站| 亚洲精品色激情综合| 成人美女网站在线观看视频| 色哟哟·www| 少妇的逼水好多| 色综合站精品国产| 中文资源天堂在线| 岛国在线免费视频观看| 免费av不卡在线播放| 狂野欧美激情性xxxx在线观看| 黄色丝袜av网址大全| 免费av毛片视频| 国产在线精品亚洲第一网站| 女生性感内裤真人,穿戴方法视频| 午夜免费男女啪啪视频观看 | 欧美丝袜亚洲另类 | 国产久久久一区二区三区| 精华霜和精华液先用哪个| 成熟少妇高潮喷水视频| 国产精品一及| 99精品在免费线老司机午夜| 成人特级av手机在线观看| 尾随美女入室| 久久6这里有精品| 一级黄色大片毛片| 久9热在线精品视频| 国产三级中文精品| 国产v大片淫在线免费观看| eeuss影院久久| 久久久久精品国产欧美久久久| 亚州av有码| 99久久精品一区二区三区| a级毛片a级免费在线| 欧美激情国产日韩精品一区| 亚洲av熟女| 在线观看免费视频日本深夜| 欧美色视频一区免费| 波多野结衣高清无吗| 久久草成人影院| 99久久精品热视频| 成人综合一区亚洲| 99九九线精品视频在线观看视频| 99久久九九国产精品国产免费| 18禁裸乳无遮挡免费网站照片| 免费在线观看影片大全网站| 在线免费观看不下载黄p国产 | 国产色爽女视频免费观看| 国产国拍精品亚洲av在线观看| 真人做人爱边吃奶动态| 日本爱情动作片www.在线观看 | 国产精品,欧美在线| 天堂网av新在线| 国产精品爽爽va在线观看网站| 不卡视频在线观看欧美| 中国美女看黄片| 国产女主播在线喷水免费视频网站 | 特大巨黑吊av在线直播| 99久久中文字幕三级久久日本| 窝窝影院91人妻| 免费大片18禁| 小说图片视频综合网站| 国产人妻一区二区三区在| 亚洲乱码一区二区免费版| 亚洲不卡免费看| 丰满人妻一区二区三区视频av| 国产激情偷乱视频一区二区| 久久久久免费精品人妻一区二区| 久久久久久久精品吃奶| 日本黄色视频三级网站网址| 国产精品美女特级片免费视频播放器| 日韩欧美精品v在线| 国内久久婷婷六月综合欲色啪| 欧美黑人欧美精品刺激| 国产亚洲精品av在线| 岛国在线免费视频观看| 成年女人看的毛片在线观看| 热99re8久久精品国产| 精华霜和精华液先用哪个| 一本精品99久久精品77| 欧美性猛交黑人性爽| 日韩精品中文字幕看吧| 久久午夜亚洲精品久久| 国产精品福利在线免费观看| 黄色丝袜av网址大全| 精品久久国产蜜桃| 成人高潮视频无遮挡免费网站| 一本久久中文字幕| 中文字幕av成人在线电影| 91麻豆av在线| 五月玫瑰六月丁香| 欧美日韩国产亚洲二区| 亚洲精品久久国产高清桃花| 久久精品人妻少妇| 特大巨黑吊av在线直播| 天堂动漫精品| 韩国av一区二区三区四区| 夜夜爽天天搞| 老熟妇仑乱视频hdxx| 中出人妻视频一区二区| 国产91精品成人一区二区三区| 中文字幕高清在线视频| 欧美一区二区亚洲| 欧美激情久久久久久爽电影| 日韩亚洲欧美综合| 亚洲精品日韩av片在线观看| 国产又黄又爽又无遮挡在线| 国产一区二区三区在线臀色熟女| 亚洲av二区三区四区| 国产精品一区www在线观看 | 黄色欧美视频在线观看| 免费观看人在逋| 一级黄色大片毛片| 在线免费观看不下载黄p国产 | 国产av麻豆久久久久久久| 一级毛片久久久久久久久女| 全区人妻精品视频| 韩国av在线不卡| 精品一区二区免费观看| 日韩欧美在线乱码| 国产精品久久视频播放| 1000部很黄的大片| 老司机福利观看| 久久人人爽人人爽人人片va| 国产探花极品一区二区| 春色校园在线视频观看| 人人妻,人人澡人人爽秒播| 高清在线国产一区| 狂野欧美激情性xxxx在线观看| 亚洲精品一卡2卡三卡4卡5卡| bbb黄色大片| 国产精品av视频在线免费观看| a级一级毛片免费在线观看| 色播亚洲综合网| 午夜福利高清视频| 九九爱精品视频在线观看| 久久久久久久午夜电影| 亚洲色图av天堂| bbb黄色大片| 黄色女人牲交| 亚洲中文日韩欧美视频| 欧美xxxx黑人xx丫x性爽| 欧美又色又爽又黄视频| 尤物成人国产欧美一区二区三区| 床上黄色一级片| 免费无遮挡裸体视频| bbb黄色大片| 国产蜜桃级精品一区二区三区| 日韩欧美在线乱码| 亚洲精品影视一区二区三区av| 国产免费一级a男人的天堂| 老司机福利观看| 九色成人免费人妻av| 亚洲色图av天堂| 美女 人体艺术 gogo| 搡女人真爽免费视频火全软件 | 少妇熟女aⅴ在线视频| 午夜福利视频1000在线观看| 亚洲精品在线观看二区| 级片在线观看| 日韩欧美免费精品| 日本-黄色视频高清免费观看| 国产单亲对白刺激| АⅤ资源中文在线天堂| 亚洲国产精品合色在线| 99精品久久久久人妻精品| 成人国产麻豆网| 亚洲午夜理论影院| 一个人观看的视频www高清免费观看|