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

    錢塘江涌潮觀測及其動(dòng)力學(xué)特性研究

    2015-08-22 08:25:49劉文虎朱小華張鐘哲張傳正馬云龍
    關(guān)鍵詞:錢塘江波速水深

    劉文虎,朱小華,張鐘哲,張傳正,馬云龍,

    龍鈺2,樊孝鵬2,朱澤南2、3,趙瑞祥2,任晨悅2

    (1.大連海洋大學(xué)海洋科技與環(huán)境學(xué)院,遼寧大連116023;2.國家海洋局第二海洋研究所 衛(wèi)星海洋環(huán)境動(dòng)力學(xué)國家重點(diǎn)實(shí)驗(yàn)室,浙江杭州310012;3.浙江大學(xué)海洋學(xué)院,浙江杭州310058)

    錢塘江涌潮觀測及其動(dòng)力學(xué)特性研究

    劉文虎1、2,朱小華2、1、3,張鐘哲1,張傳正2,馬云龍1、2,

    龍鈺2,樊孝鵬2,朱澤南2、3,趙瑞祥2,任晨悅2

    (1.大連海洋大學(xué)海洋科技與環(huán)境學(xué)院,遼寧大連116023;2.國家海洋局第二海洋研究所 衛(wèi)星海洋環(huán)境動(dòng)力學(xué)國家重點(diǎn)實(shí)驗(yàn)室,浙江杭州310012;3.浙江大學(xué)海洋學(xué)院,浙江杭州310058)

    為了研究錢塘江涌潮的動(dòng)力特征,于2013年11月5—6日在錢塘江杭州段采用沿海聲層析 (Coastal acoustic tomography,CAT)和溫鹽深儀 (Conductivity temperature depth,CTD)進(jìn)行了涌潮觀測試驗(yàn)。兩個(gè)觀測站位分別設(shè)置在沿江兩岸,站間距為3050 m,試驗(yàn)歷時(shí)26 h,成功觀測到兩次涌潮的全過程。結(jié)果表明:兩次涌潮到來前后,由雙向傳播時(shí)間差計(jì)算得到的斷面平均流速的變化范圍分別為0.39~-1.08、0.31~-0.85 m/s(東向流為正),水位升高的最大幅度分別為1.38、1.14 m;兩次涌潮通過觀測站位的波速分別為-8.69、-8.43 m/s;觀測斷面流量的變化范圍為-7753.00~3504.49 m3/s,整個(gè)試驗(yàn)期間的平均流量為644.14 m3/s;涌潮的弗勞德數(shù)介于1.19~1.29之間,表明觀測涌潮為波狀涌潮;涌潮到達(dá)前后能量的損耗率小于0.6%,符合波狀涌潮的能量守恒關(guān)系。

    錢塘江;涌潮;沿海聲層析;動(dòng)力特性

    錢塘江是浙江省第一大河,發(fā)源于安徽南部,流經(jīng)皖、浙兩省,最后經(jīng)杭州灣注入東海,全長600多千米,流域面積約為55 558 km2。涌潮是入??诤涌诙蔚囊环N潮水暴漲現(xiàn)象[1]。當(dāng)海水漲潮時(shí),潮水涌入錢塘江漏斗狀河口,潮水堆積會(huì)形成涌潮。每逢大潮期間,錢塘江都會(huì)形成雄偉壯觀的涌潮。涌潮在傳播過程中,由于河床地形和河道平面形態(tài)變化而呈現(xiàn)不同的形態(tài)潮景,最具代表性的有 “一線潮”、 “交叉潮” 和 “回頭潮等[1]。 涌潮蘊(yùn)藏著巨大的能量,同時(shí)也是一種自然災(zāi)害。它對(duì)河道、岸堤和涉水建筑沖擊嚴(yán)重,不但影響航運(yùn)交通和沿江兩岸的開發(fā),而且對(duì)兩岸人民的生命和財(cái)產(chǎn)也構(gòu)成潛在危脅。因此,研究涌潮對(duì)開發(fā)涌潮和防御潮災(zāi)都具有重要意義。

    目前國內(nèi)常用數(shù)值模擬的方法研究涌潮的宏觀性質(zhì)和微觀性質(zhì)。在涌潮宏觀性質(zhì)方面,蘇銘德等[2-3]采用激波捕捉法初步模擬了錢塘江涌潮的形成和發(fā)展過程;李紹武等[4]建立了準(zhǔn)三維涌潮數(shù)學(xué)模型,模擬河口區(qū)涌潮流場;潘存鴻等[5]應(yīng)用基于KFVS格式的二維數(shù)值模型模擬了錢塘江涌潮形成、發(fā)展和衰減的全過程。在涌潮微觀性質(zhì)方面,曾劍等[6]用徑向基函數(shù) (RBF)建立模型預(yù)測了鹽官站涌潮潮頭的高度;楊火其等[7]通過水槽試驗(yàn)建立了涌潮流速與潮前水深、涌潮高度和涌潮傳播速度的經(jīng)驗(yàn)關(guān)系;潘存鴻等[5]結(jié)合實(shí)測數(shù)據(jù)分析了涌潮高度、流速、形態(tài)和壓力等涌潮特性;謝東風(fēng)等[8]建立了涌潮潮頭坡度的解析式,復(fù)演了涌潮到達(dá)時(shí)刻潮頭急劇增加的過程;林炳堯[1]綜合分析了錢塘江涌潮的特性。在國外,Wolanskia等[9]基于實(shí)測數(shù)據(jù)分析了澳大利亞戴利河波狀涌潮的動(dòng)力特性;Bonneton等[10]首次在法國加倫河上觀測了大振幅涌潮,研究了涌潮的流態(tài)和斷面特征;Chanson[11]結(jié)合觀測和數(shù)值模擬對(duì)涌潮進(jìn)行了詳細(xì)的綜述。

    由于涌潮能量巨大,并且錢塘江涌潮段水淺流急,使用常規(guī)儀器觀測流速具有一定困難,故有關(guān)涌潮流速的資料相對(duì)匱乏[1]。目前,沿海聲層析(Coastal acoustic tomography,CAT)為涌潮流速觀測提供了新的方法,這是一種利用雙向聲傳播方法觀測流速的儀器,該儀器已經(jīng)在中國舟山、臺(tái)州三門灣、膠州灣、瓊州海峽和錢塘江等地區(qū)成功完成了測流試驗(yàn)[12-16]。張傳正等[15]首次利用CAT方法對(duì)錢塘江局部流速進(jìn)行了觀測,通過傳播時(shí)間差計(jì)算的流速與多普勒流速剖面儀 (Acoustic doppler profiler,ADP)走航觀測的結(jié)果基本一致,驗(yàn)證了CAT觀測涌潮河段流速的可行性。隨后,Zhu等[16]根據(jù)傳播時(shí)間差與走航ADP觀測流量建立了經(jīng)驗(yàn)關(guān)系,得到了試驗(yàn)期間斷面流量的時(shí)間序列,成功捕捉到了因錢塘江涌潮而引起的流量變化。本研究中,利用CAT和溫鹽深儀 (Conductivity temperature depth,CTD)對(duì)錢塘江杭州段涌潮進(jìn)行了觀測試驗(yàn),對(duì)涌潮動(dòng)力學(xué)特性進(jìn)行了分析。

    1 觀測地點(diǎn)與方法

    涌潮觀測試驗(yàn)在距離杭州灣大約90 km的錢塘江杭州段進(jìn)行。試驗(yàn)從2013年11月5日11:00開始到次日13:00結(jié)束,歷時(shí)26 h。觀測站位E站和W站的布置如圖1所示,E站位于復(fù)興大橋的橋墩上面,W站位于靠近錢塘江大橋附近的閘口水文站,兩站相距3050 m,站位間的平均水深約為6 m。兩站各設(shè)置一套CAT系統(tǒng),儀器連接GPS和收發(fā)兩用的換能器,換能器置于水中約3 m處。試驗(yàn)使用聲信號(hào)為5 kHz的10階M序列信號(hào),發(fā)聲間隔為15 s,換能器同時(shí)接收對(duì)方發(fā)射的聲信號(hào),數(shù)據(jù)自動(dòng)儲(chǔ)存在內(nèi)部SD卡中。

    圖1 試驗(yàn)站位設(shè)置示意圖Fig.1 Schematic view of experimental set-up

    通過對(duì)數(shù)據(jù)的相關(guān)處理,可以獲得對(duì)方站位發(fā)射的聲信號(hào)到達(dá)的時(shí)間。假設(shè)沿著聲傳播路徑的平均流速為Vm,聲音在水中傳播的平均速度為Cm,則從W站到E站的傳播時(shí)間τ1和從E站到W站的傳播時(shí)間τ2可以表示為[17]

    其中:L為站位之間的距離。將式 (1)和 (2)進(jìn)行簡單推導(dǎo)可得:

    其中:ˉτ=(τ1+τ2)/2, Δτ=τ2-τ1。 假設(shè) ˉτ≈τ1≈τ2,對(duì)式 (3)和式 (4)進(jìn)行微分和變換,得到Vm的相對(duì)誤差為

    由于GPS時(shí)鐘精確同步,ˉτ?δτ,則Δτ?δ(Δˉτ), 故式 (5)中第二項(xiàng)和第三項(xiàng)可以忽略,式 (5)可簡化為

    在本次試驗(yàn)中,兩站位間距L=3 km,GPS誤差為10 m,流速相對(duì)誤差δVm/Vm=0.003 3,即當(dāng)平均流速為Vm=1 m/s時(shí),對(duì)應(yīng)的誤差值為0.003 3 m/s。這樣,在保證高精度傳播時(shí)間的觀測條件下,由距離誤差導(dǎo)致平均流速的觀測誤差極小。

    由于兩個(gè)站位離岸距離相對(duì)于整個(gè)江面的寬度(1200 m)都較小且站位靠近岸邊的水深較淺,本研究中假定觀測站位距岸邊的流量忽略不計(jì)。根據(jù)張傳正[18]提出的觀測區(qū)域雙向聲信號(hào)傳播時(shí)間差Δτ與流量Q的經(jīng)驗(yàn)關(guān)系式對(duì)斷面流量進(jìn)行推算:

    為了觀測涌潮期間水位、溫度和鹽度的變化,本研究中在兩個(gè)觀測站位處各放置了一臺(tái)錨系CTD(RBR,XR-420)。CTD距河床高度為0.35 m,采樣間隔為2 s。

    如圖2所示,涌潮沿錢塘江自下游向上游傳播,通常把潮波傳播的速度稱為波速。這里將涌潮波速和涌潮高度的計(jì)算近似地視為一維問題。涌潮波速可用下面的理論公式[9]進(jìn)行估算:

    其中:h1和h2分別為涌潮到來前后的水深 (m); g為重力加速度 (m/s2)。根據(jù)波速的定義,實(shí)際觀測波速的表達(dá)式如下:

    其中:R為涌潮從E站到W站的傳播距離;t1和t2分別為涌潮到達(dá)E站和W站的時(shí)間。

    圖2 涌潮示意圖Fig.2 Sketch of the tidal bores

    涌潮高度是衡量涌潮大小的重要參數(shù)之一。用Δη=h2-h1表示涌潮前后高度,Δv=v2-v1表示斷面平均流速差,則涌潮前后流速差與涌潮高度的理論關(guān)系為[17]

    其中,Δv可通過傳播時(shí)間差計(jì)算得到,涌潮前后的水深可通過CTD觀測得到。假設(shè)涌潮前后能量守恒,則可以通過涌潮前后的流速差來推測當(dāng)時(shí)理論上的涌潮高度。

    2 結(jié)果與分析

    2.1 聲信號(hào)

    圖3為W站和E站接收到信號(hào)信噪比 (Signal to noise ratio,SNR)的時(shí)間序列。試驗(yàn)觀測顯示, SNR在2.056~2.066 s之間有明顯的峰值,SNR的峰值介于100~300之間,SNR最大峰值對(duì)應(yīng)的時(shí)間認(rèn)為是聲信號(hào)到達(dá)的時(shí)間。從圖3可知,聲信號(hào)實(shí)際到達(dá)時(shí)間與預(yù)期到達(dá)時(shí)間基本一致 (3050/ 1480=2.06 s,3050 m為兩站距離,1480 m/s為根據(jù)現(xiàn)場觀測的溫度、鹽度和深度計(jì)算所得到的聲波在水中的傳播速度)。由于試驗(yàn)期間W站在11月5日16:00至17:30之間GPS天線發(fā)生故障,故出現(xiàn)少量數(shù)據(jù)缺失。

    圖3 信噪比時(shí)間序列Fig.3 Tim e series of signal-to-noise ratio

    2.2 水深、鹽度、溫度、聲傳播時(shí)間、流速和流量的觀測

    圖4-A為試驗(yàn)期間在E站和W站用CTD觀測得到的水深時(shí)間序列。第一次涌潮到來時(shí),E站和W站水位升高的最大幅度分別為1.31、1.38 m;第二次涌潮到來時(shí),E站和W站水位升高的最大幅度分別為0.94、1.14 m。根據(jù)涌潮到來前后的水深,由式 (8)計(jì)算得到兩次涌潮到來時(shí)E站和W站的理論波速的平均值分別為8.62、8.54 m/s。另外,根據(jù)站位水深的變化可確定涌潮到達(dá)的時(shí)刻,由式 (9)可得到實(shí)際波速分別為8.69、8.43 m/s。理論波速與實(shí)際波速基本一致。

    圖4-B為在E站和W站用CTD觀測得到的鹽度時(shí)間序列。第一次涌潮開始約2 h,E站鹽度由0.4逐漸增大至1.34,W站鹽度由0.4逐漸增大至0.71。表明第一次涌潮期間,有海水從下游入侵到上游。第二次涌潮過后,觀測區(qū)域的鹽度沒有出現(xiàn)明顯的變化。

    圖4-C為在E站和W站用CTD觀測得到的水溫時(shí)間序列。在整個(gè)試驗(yàn)觀測期間,水溫變化不大,其變化范圍為19.4~20.2℃。涌潮過后,水溫整體有降低的趨勢。

    圖4-D為根據(jù)CAT觀測得到的雙向聲信號(hào)平均傳播時(shí)間的時(shí)間序列。整個(gè)試驗(yàn)期間,雙向聲信號(hào)的平均傳播時(shí)間為2.058~2.062 s。兩次涌潮到來前后,雙向聲信號(hào)平均傳播時(shí)間均發(fā)生較大變化,與之相對(duì)應(yīng)的是聲速減小。如果傳播時(shí)間的增大是由溫度變化引起,則溫度變化必須降低大約0.52℃,而從圖中可以看出,涌潮過后水溫變化不大,而且僅第一次涌潮過后溫度稍微升高,因此,雙向聲信號(hào)的平均傳播時(shí)間增大的主要原因不是溫度的變化,而是涌潮的到來讓換能器的相對(duì)位置發(fā)生了一定變化。

    圖4 水深、鹽度、溫度、雙向聲信號(hào)平均傳播時(shí)間和計(jì)算得到的流速和流量時(shí)間序列Fig.4 Time series of water depth,salinity,tem perature,mean reciprocal travel time,velocity and calculated discharge

    根據(jù)式 (3)和雙向聲信號(hào)傳播時(shí)間差求得斷面平均流速的時(shí)間序列如圖4-E所示。兩次涌潮到來之前,流速基本穩(wěn)定在0.39、0.31 m/s。涌潮到來之后,涌潮流 (西向流)與上游淡水流(東向流)的流速相互抵消,流向發(fā)生逆轉(zhuǎn),流速在20 min內(nèi)達(dá)到最大。兩次涌潮中觀測到的最大流速分別為-1.08、-0.85 m/s。流向逆轉(zhuǎn)為西向流之后大約1.5 h,漲潮流轉(zhuǎn)變成落潮流,流速逐漸減小,直到流向轉(zhuǎn)向東向流之后流速緩慢增大,并逐漸趨于平穩(wěn)。從涌潮開始到結(jié)束的整個(gè)作用過程大約持續(xù)了5 h 40 min。

    根據(jù)式 (7)和雙向傳播時(shí)間差求得斷面流量的時(shí)間序列如圖4-F所示。試驗(yàn)觀測期間,最大東向流 (最大落潮流量)出現(xiàn)在無涌潮影響期間,最大西向流 (最大漲潮流量)出現(xiàn)在第一次涌潮到來之后。對(duì)缺失數(shù)據(jù)進(jìn)行插值,得到觀測期間最大東向流量和西向流量分別為3504.49、-7753.00 m3/s,平均流量為644.14 m3/s。兩次涌潮帶來的西向輸運(yùn)量分別為9.09×107、6.83×107m3,相當(dāng)于5~6個(gè)西湖的總蓄水量 (1.43×107m3)[19]。

    3 討論

    為了判別涌潮的類型,本研究中引入了弗勞德數(shù) (Froude number,Fr),其為水躍強(qiáng)度的動(dòng)力參數(shù),可用潮流速與波速間的關(guān)系表示[9]:

    當(dāng)Fr<1時(shí),潮波推進(jìn)的波列未形成涌潮;當(dāng)1≤Fr<1.7時(shí),涌潮為波狀涌潮;當(dāng)Fr≥1.7時(shí),流態(tài)發(fā)生根本變化,不再具有波動(dòng)性,變成破碎涌潮。根據(jù)觀測數(shù)據(jù)求得兩次涌潮的 Fr=1.19~1.29,表明觀測期間的涌潮類型屬于波狀涌潮。

    當(dāng)波速和流速相等時(shí)的水流為臨界流,臨界流對(duì)應(yīng)的水深即為臨界水深,臨界水流對(duì)應(yīng)的Fr為1,臨界水深 (dc)可由潮流速(v1)、水深(h1)和波速(C)表示[11]:

    當(dāng)實(shí)際水深大于臨界水深時(shí),對(duì)應(yīng)的水流為亞臨界流,流速小于臨界流速,則Fr<1;當(dāng)實(shí)際水深小于臨界水深時(shí),對(duì)應(yīng)的水流為超臨界流,流速大于臨界流速,則Fr>1。在觀測的流場中,計(jì)算得出兩次涌潮中E站臨界水深分別為4.98、4.80 m, W站臨界水深分別為6.63、6.48 m,E站和W站的臨界水深均介于涌潮前后觀測的水深之間,說明涌潮導(dǎo)致流向逆轉(zhuǎn),本次觀測到的兩次涌潮流場中均存在亞臨界流和超臨界流轉(zhuǎn)換的流動(dòng)。

    涌潮的種類不同,其能量的耗散也不同。根據(jù)涌潮前后的水深關(guān)系,可得出涌潮前后單位水體的能量差,且涌潮的能量耗散為[20]

    式 (13)表明,涌潮的能量耗散與Fr有關(guān)。假如涌潮發(fā)生前后的能量耗散忽略不計(jì),可認(rèn)為能量近似守恒。根據(jù)能量和動(dòng)量守恒得到如下關(guān)系[11]:

    其中:M為動(dòng)量;E為單位質(zhì)量的能量 (比能)。

    圖5-A為根據(jù)式 (13)計(jì)算得到的能量耗散和Fr的關(guān)系。當(dāng) Fr<1.2時(shí),能量耗散率小于0.22%;當(dāng)Fr=1.2~2.0時(shí),能量損耗率會(huì)增加到9.08%,增幅較快;當(dāng)Fr>2.0時(shí),能量損耗增加趨勢變緩。對(duì)錢塘江涌潮的觀測表明,涌潮能量損失小于0.6%,可以忽略不計(jì),涌潮到來前后的能量基本守恒。圖5-B為根據(jù)式 (14)、(15)計(jì)算得到本次觀測涌潮能量與動(dòng)量的關(guān)系圖,可以看出,試驗(yàn)所觀測的數(shù)據(jù)完全處于參數(shù)化的M-E曲線上,表明觀測數(shù)據(jù)滿足能量近似守恒關(guān)系。

    圖6為根據(jù)式 (10)計(jì)算得到的涌潮前后流速差與涌潮高度的關(guān)系,可以看出,流速差越大對(duì)應(yīng)的潮頭高度就越高,試驗(yàn)觀測值與理論值基本吻合。

    4 結(jié)論

    2013年11月5—6日在錢塘江杭州段進(jìn)行了一次涌潮觀測試驗(yàn),成功地觀測到兩次涌潮期間水位、流速的變化,得到主要結(jié)論如下:

    (1)兩次涌潮發(fā)生時(shí),涌潮前后觀測斷面流速的變化范圍分別為 0.39~-1.08 m/s和0.31~-0.85 m/s(東向流為正),水位升高的最大幅度分別為1.38、1.14 m,兩次涌潮通過觀測站位的波速分別為-8.69、-8.43 m/s。

    圖5 涌潮能量耗散率與弗勞德數(shù)及涌潮能量與動(dòng)量的關(guān)系Fig.5 General relationship between energy dissipation rates and Froude num ber,and the relationship between energy and momentum of the tidal bores

    圖6 流速差與涌潮高度關(guān)系Fig.6 Relationship between velocity difference and wave height

    (2)觀測期間,最大東向流 (最大落潮流量)出現(xiàn)在無涌潮影響期間,最大西向流 (最大漲潮流量)出現(xiàn)在涌潮導(dǎo)致江水逆流階段,最大東向流和西向流流量分別為3504.49、-7753.00 m3/s,平均流量為644.14 m3/s。兩次涌潮帶來的西向輸運(yùn)量分別為9.09×107、6.83×107m3。

    (3)涌潮的弗勞德數(shù)為1.19~1.29,表明觀測的涌潮是波狀涌潮。臨界水深均介于涌潮前后的水深之間,涌潮流場中存在亞臨界流和超臨界流轉(zhuǎn)換的流動(dòng)。

    (4)涌潮到達(dá)前后能量的損耗率小于0.6%,表明波狀涌潮的能量損耗可基本忽略不計(jì),觀測數(shù)據(jù)都處于參數(shù)化的M-E曲線上,涌潮的能量基本保持守恒。

    [1] 林炳堯.錢塘江涌潮特性[M].北京:海洋出版社,2008:1-52.

    [2] 蘇銘德,徐聽,朱錦林,等.數(shù)值模擬在錢塘江涌潮分析中的應(yīng)用:Ⅰ.數(shù)值計(jì)算方法[J].力學(xué)學(xué)報(bào),1999,31(5):521-533.

    [3] 蘇銘德,徐聽,朱錦林,等.數(shù)值模擬在錢塘江涌潮分析中的應(yīng)用:Ⅱ.計(jì)算結(jié)果和分析[J].力學(xué)學(xué)報(bào),1999,31(6):700-716.

    [4] 李紹武,盧麗鋒,時(shí)鐘.河口準(zhǔn)三維涌潮數(shù)學(xué)模型研究[J].水動(dòng)力學(xué)研究與進(jìn)展A輯,2004,19(4):407-415.

    [5] 潘存鴻,魯海燕,曾劍.錢塘江涌潮特性及其數(shù)值模擬[J].水利水運(yùn)工程學(xué)報(bào),2008,6(2):1-9.

    [6] 曾劍,孫志林,熊紹隆,等.基于徑向基函數(shù)網(wǎng)絡(luò)的錢塘江涌潮模型及其應(yīng)用[J].浙江大學(xué)學(xué)報(bào):工學(xué)版,2007,41(9):1457-1461.

    [7] 楊火其,潘存鴻,周建炯,等.涌潮水動(dòng)力學(xué)特性試驗(yàn)研究[J].水電能源科學(xué),2008,26(4):136-138.

    [8] 謝東風(fēng),潘存鴻,陸波,等.基于實(shí)測數(shù)據(jù)的錢塘江水動(dòng)力學(xué)特性研究[J].水動(dòng)力學(xué)研究與進(jìn)展A輯,2012,27(5):501-508.

    [9] Wolanski E,Williams D,Spagnol S,et al.Undular tidal bore dynamics in the Daly Estuary,Northern Australia[J].Estuarine, Coastal and Shelf Science,2004,60(4):629-636.

    [10] Bonneton P,Van de Loock J,Parisot JP,et al.On the occurrence of tidal bores-the Garonne River case[J].Journal of Coastal Research,2011,SI64:1462-1466.

    [11] Chanson H.Tidal bores,Aegir,Eagre,Mascaret,Pororoca:Theory and Observations[M].Singapore:World Scientific Publishing Company,2012:9-20

    [12] Zhu X H,Kaneko A,Wu Q,etal.Mapping tidal current structures in Zhitouyang Bay,China,using coastal acoustic tomography[J]. IEEE Journal of Oceanic Engineering,2013,38(2):285-296.

    [13] 朱澤南,朱小華,張傳正,等.三門灣沿海聲層析潮流觀測實(shí)驗(yàn)[J].地球物理學(xué)報(bào),2015,58(4):1-12.

    [14] Zhu X H,Zhu Z,Guo X,et al.Measurement of tidal and residual currents and volume transport through the Qiongzhou Strait using coastal acoustic tomography[J].Continental Shelf Research, 2015,108:65-75.

    [15] 張傳正,朱小華,吳清松.錢塘江雙向聲傳播測流實(shí)驗(yàn)研究[J].地球物理學(xué)報(bào),2010,53(7):1741-1748.

    [16] Zhu X H,Zhang C,Wu Q,et al.Measuring discharge in a river with tidal bores by use of the coastal acoustic tomography system [J].Estuarine,Coastal and Shelf Science,2012,104/105:54-55.

    [17] Munk W,Worcester P,Wunsch C.Ocean Acoustic Tomography [M].New York:Cambridge University Press,1995.

    [18] 張傳正.沿海聲層析技術(shù)在我國近海的應(yīng)用研究[D].杭州:國家海洋局第二海洋研究所,2011.

    [19] 杭州市地方志編纂委員會(huì).杭州年鑒[M].北京:方志出版社,2004.

    [20] 張瑋,徐金環(huán),李國臣,等.涌潮分析與波速計(jì)算[J].水利水運(yùn)科學(xué)研究,1999,6(2):158-164.

    Observation and dynam ic characteristics of tidal bore in Qiantang River,China

    LIUWen-hu1,2,ZHU Xiao-hua2,1,3,ZHANG Zhong-zhe1,ZHANG Chuan-zheng2,MA Yun-long1,2, LONG Yu2,FAN Xiao-peng2,ZHU Ze-nan2,3,ZHAO Rui-xiang2,REN Chen-yue2
    (1.College of Marine Science and Environment,Dalian Ocean University,Dalian 116023,China;2.State Key Laboratory of Satellite Ocean Environment Dynamics,Second Institute of Oceanography,State Oceanic Administration,Hangzhou310012,China;3.Ocean College,Zhejiang University,Hangzhou 310058,China)

    A 26-hour observation of tidal boreswas carried out by a coastal acoustic tomography(CAT)system and conductivity temperature depth(CTD)in two stations at both sides with a distance of 3050 m diagonally across Qiantang River from 5 to 6 November,2013 to valuate the dynamic characteristics of tidal bore in Qiantang River, China.The range-averaged current velocities calculated from the reciprocal sound transmission data were changed from 0.39 m/s to-1.08 m/s before arrival of tidal bore with themaximalwater level changes of 1.38 m and from 0.31 m/s to-0.85 m/s after arrival of tidal borewith themaximalwater level changes of1.14 m(the positive velocity towards the east).Themean phase speeds of two tidal bores were estimated to be-8.69 m/s and-8.43 m/s at the two stations,respectively.During the whole experiment,the volume transport was varied from -7753.00 m3/s to 3504.49 m3/s through the section,with amean of644.14 m3/s.The Froude numbers of tidal bores were estimated to be from 1.19 to 1.29,indicating that the observed tidal boreswere undular bore,with tidal bore energy losses of less than 0.6%,conforming to energy conservation relations for undular bores.

    Qiantang River;tidal bore;coastal acoustic tomography system;dynamic characteristics

    P731.23

    A

    2015-01-12

    國家自然科學(xué)基金資助項(xiàng)目 (41276095,41476020,41321004,41176021,41576001);國家海洋局第二海洋研究所基本科研業(yè)務(wù)費(fèi)專項(xiàng) (JT1402);衛(wèi)星海洋環(huán)境動(dòng)力學(xué)國家重點(diǎn)實(shí)驗(yàn)室自主項(xiàng)目 (SOEDZZ1403,SOEDZZ1501);國家海洋局第二海洋研究所SOED開放課題基金資助項(xiàng)目 (SOED1511);遼寧省教育廳科研項(xiàng)目 (L2012268)

    劉文虎 (1989—),男,碩士研究生。E-mail:liuwh026@163.com

    朱小華 (1963—),男,博士,研究員。E-mail:xhzhu@sio.org.cn

    張鐘哲 (1965—),男,博士,教授。E-mail:zzz@dlou.edu.cn

    10.16535/j.cnki.dlhyxb.2015.05.022

    2095-1388(2015)05-0567-06

    猜你喜歡
    錢塘江波速水深
    為什么錢塘江的浪潮格外壯觀
    書法靜水深流
    河北水利(2022年10期)2022-12-29 11:48:12
    基于水深分段選擇因子的多光譜影像反演水深
    我在錢塘江邊長大
    基于實(shí)測波速探討地震反射波法超前預(yù)報(bào)解譯標(biāo)志
    錢塘江觀潮
    小讀者(2021年2期)2021-03-29 05:03:18
    吉林地區(qū)波速比分布特征及構(gòu)造意義
    GPS RTK技術(shù)在水深測量中的應(yīng)用
    錢塘江海塘及涌潮景觀申遺研究
    中國水利(2015年4期)2015-02-28 15:12:25
    浸入式水深監(jiān)測儀器的設(shè)計(jì)
    亚洲,欧美精品.| 精品久久久久久久久久久久久| 一本久久中文字幕| 日本一二三区视频观看| 在线观看免费日韩欧美大片| 亚洲成av人片免费观看| 99在线视频只有这里精品首页| 国产av又大| 热99re8久久精品国产| 精品一区二区三区av网在线观看| 欧美日韩精品网址| 他把我摸到了高潮在线观看| 男女之事视频高清在线观看| 久久久国产成人精品二区| 午夜日韩欧美国产| 国产精品久久久久久亚洲av鲁大| 亚洲狠狠婷婷综合久久图片| 国内精品久久久久久久电影| 啦啦啦韩国在线观看视频| 国语自产精品视频在线第100页| 午夜福利18| 十八禁网站免费在线| 久久精品成人免费网站| 亚洲精品国产一区二区精华液| 香蕉久久夜色| 久久国产乱子伦精品免费另类| 色综合站精品国产| 国产精品美女特级片免费视频播放器 | 亚洲七黄色美女视频| 国内精品久久久久久久电影| 国产精品亚洲美女久久久| 身体一侧抽搐| 国产精品九九99| 亚洲色图 男人天堂 中文字幕| 91大片在线观看| 777久久人妻少妇嫩草av网站| www日本黄色视频网| 亚洲无线在线观看| 黄色视频不卡| 欧美性长视频在线观看| 亚洲最大成人中文| 亚洲va日本ⅴa欧美va伊人久久| av免费在线观看网站| 又黄又爽又免费观看的视频| 国产成人精品无人区| 少妇裸体淫交视频免费看高清 | 特级一级黄色大片| 国产真实乱freesex| 午夜免费激情av| 成年免费大片在线观看| 在线永久观看黄色视频| 不卡一级毛片| 日本黄大片高清| 亚洲精品中文字幕在线视频| 精品人妻1区二区| 97碰自拍视频| 国产真人三级小视频在线观看| 男女午夜视频在线观看| aaaaa片日本免费| 国产成人aa在线观看| 欧美成狂野欧美在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜免费激情av| 亚洲人与动物交配视频| 婷婷丁香在线五月| 亚洲在线自拍视频| 香蕉久久夜色| 一边摸一边做爽爽视频免费| 男人舔女人的私密视频| 一a级毛片在线观看| 一二三四社区在线视频社区8| 日韩欧美免费精品| 精品久久久久久久末码| 一进一出抽搐动态| 可以在线观看的亚洲视频| 欧美一级a爱片免费观看看 | 老鸭窝网址在线观看| 欧美日韩亚洲国产一区二区在线观看| 小说图片视频综合网站| 日韩三级视频一区二区三区| 亚洲成a人片在线一区二区| 日韩欧美在线乱码| 男女午夜视频在线观看| 男女视频在线观看网站免费 | 免费在线观看完整版高清| 在线观看免费日韩欧美大片| 在线国产一区二区在线| 一级a爱片免费观看的视频| 亚洲av熟女| 蜜桃久久精品国产亚洲av| 白带黄色成豆腐渣| 99热只有精品国产| 麻豆国产av国片精品| 国产精品电影一区二区三区| 91字幕亚洲| 国产精品自产拍在线观看55亚洲| 国产v大片淫在线免费观看| 亚洲国产精品合色在线| 国产蜜桃级精品一区二区三区| www.www免费av| 可以在线观看的亚洲视频| 国内精品久久久久久久电影| 亚洲成人国产一区在线观看| 在线免费观看的www视频| 变态另类丝袜制服| 久久久精品欧美日韩精品| 国产精品爽爽va在线观看网站| 欧美成人一区二区免费高清观看 | 18美女黄网站色大片免费观看| 非洲黑人性xxxx精品又粗又长| 伦理电影免费视频| 亚洲一区二区三区不卡视频| 久久精品91蜜桃| 欧美av亚洲av综合av国产av| 一二三四社区在线视频社区8| 亚洲黑人精品在线| 成年免费大片在线观看| 看免费av毛片| 国产精品综合久久久久久久免费| 免费看美女性在线毛片视频| 美女扒开内裤让男人捅视频| 国产99白浆流出| 久久香蕉国产精品| 特大巨黑吊av在线直播| 亚洲 欧美 日韩 在线 免费| 五月玫瑰六月丁香| 亚洲免费av在线视频| 不卡av一区二区三区| 国内揄拍国产精品人妻在线| 一级作爱视频免费观看| 搡老熟女国产l中国老女人| 五月伊人婷婷丁香| 久久久久国内视频| 午夜日韩欧美国产| 在线观看一区二区三区| 精品欧美国产一区二区三| 亚洲av中文字字幕乱码综合| 美女午夜性视频免费| 久久久久国内视频| 欧美日本视频| 国产av一区二区精品久久| 国产精华一区二区三区| 天堂av国产一区二区熟女人妻 | tocl精华| 又黄又粗又硬又大视频| 99热只有精品国产| e午夜精品久久久久久久| 免费在线观看视频国产中文字幕亚洲| 国产精品免费一区二区三区在线| 欧美丝袜亚洲另类 | 美女大奶头视频| 一个人观看的视频www高清免费观看 | 国产69精品久久久久777片 | 宅男免费午夜| 国产精品美女特级片免费视频播放器 | 日本熟妇午夜| 日本精品一区二区三区蜜桃| 国产单亲对白刺激| 久久精品91蜜桃| 亚洲av日韩精品久久久久久密| 可以免费在线观看a视频的电影网站| 亚洲色图av天堂| 欧美黑人欧美精品刺激| 很黄的视频免费| 757午夜福利合集在线观看| 久久中文看片网| 狂野欧美白嫩少妇大欣赏| 超碰成人久久| 变态另类丝袜制服| 狂野欧美激情性xxxx| 亚洲精品一卡2卡三卡4卡5卡| 国产欧美日韩精品亚洲av| av片东京热男人的天堂| 日本黄色视频三级网站网址| 免费观看人在逋| 级片在线观看| 亚洲精品在线观看二区| a在线观看视频网站| 一级毛片精品| 在线永久观看黄色视频| 欧美3d第一页| 色哟哟哟哟哟哟| 成在线人永久免费视频| 久久国产精品人妻蜜桃| 国产精品久久久av美女十八| 日本一二三区视频观看| 成人高潮视频无遮挡免费网站| 欧美黄色片欧美黄色片| 一卡2卡三卡四卡精品乱码亚洲| 夜夜看夜夜爽夜夜摸| 国产一区二区三区在线臀色熟女| 欧洲精品卡2卡3卡4卡5卡区| 精品电影一区二区在线| 麻豆av在线久日| 亚洲精品av麻豆狂野| 一级a爱片免费观看的视频| 午夜福利视频1000在线观看| www日本在线高清视频| 国产精品免费一区二区三区在线| 婷婷丁香在线五月| 日韩欧美国产一区二区入口| 亚洲国产精品成人综合色| 黄色a级毛片大全视频| 人人妻人人澡欧美一区二区| 久久久久九九精品影院| 91在线观看av| 亚洲熟妇熟女久久| 日本在线视频免费播放| 久久精品人妻少妇| 午夜福利视频1000在线观看| 老熟妇乱子伦视频在线观看| 国产亚洲欧美在线一区二区| 一本精品99久久精品77| 国产男靠女视频免费网站| 国模一区二区三区四区视频 | 999久久久国产精品视频| 最近视频中文字幕2019在线8| 很黄的视频免费| 又黄又粗又硬又大视频| 男女床上黄色一级片免费看| 黄片大片在线免费观看| 日日爽夜夜爽网站| 久久中文看片网| 操出白浆在线播放| 香蕉久久夜色| 日韩成人在线观看一区二区三区| 亚洲人成77777在线视频| 国产视频内射| 又黄又粗又硬又大视频| 日本 av在线| 亚洲七黄色美女视频| 成人亚洲精品av一区二区| 婷婷亚洲欧美| 香蕉av资源在线| 国产成人精品久久二区二区91| 99精品欧美一区二区三区四区| 久久国产精品影院| 伊人久久大香线蕉亚洲五| 亚洲中文av在线| 首页视频小说图片口味搜索| 9191精品国产免费久久| 国产午夜福利久久久久久| 亚洲国产欧洲综合997久久,| 非洲黑人性xxxx精品又粗又长| 国内精品久久久久精免费| 黄色 视频免费看| 亚洲人成77777在线视频| 国产精品亚洲av一区麻豆| 国产亚洲精品一区二区www| 在线视频色国产色| 亚洲无线在线观看| 国内精品久久久久久久电影| 亚洲一区二区三区不卡视频| 丁香欧美五月| 51午夜福利影视在线观看| 精品久久久久久成人av| 免费无遮挡裸体视频| 男女做爰动态图高潮gif福利片| 欧美成人一区二区免费高清观看 | 在线看三级毛片| 国产亚洲精品一区二区www| 一级a爱片免费观看的视频| 一个人观看的视频www高清免费观看 | 夜夜爽天天搞| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美一级毛片孕妇| 精品人妻1区二区| 国产av又大| 婷婷精品国产亚洲av| 舔av片在线| 夜夜躁狠狠躁天天躁| 国产一区在线观看成人免费| 欧美中文日本在线观看视频| 精品人妻1区二区| 后天国语完整版免费观看| 777久久人妻少妇嫩草av网站| 叶爱在线成人免费视频播放| 精品久久久久久久人妻蜜臀av| 宅男免费午夜| 久久久久精品国产欧美久久久| 2021天堂中文幕一二区在线观| 好男人在线观看高清免费视频| 国产精品一区二区精品视频观看| 69av精品久久久久久| 色综合婷婷激情| 精品乱码久久久久久99久播| 国产av一区在线观看免费| 欧美日本亚洲视频在线播放| 老司机午夜十八禁免费视频| 国产爱豆传媒在线观看 | 又粗又爽又猛毛片免费看| 在线观看午夜福利视频| 一二三四在线观看免费中文在| videosex国产| 日日摸夜夜添夜夜添小说| 国产91精品成人一区二区三区| 首页视频小说图片口味搜索| 午夜影院日韩av| 久久性视频一级片| 亚洲午夜理论影院| 美女高潮喷水抽搐中文字幕| 熟女少妇亚洲综合色aaa.| 国产av在哪里看| 亚洲成人久久爱视频| 国产99白浆流出| 欧美成人一区二区免费高清观看 | 午夜a级毛片| 成人国语在线视频| 亚洲美女黄片视频| xxxwww97欧美| 人人妻,人人澡人人爽秒播| 国产欧美日韩一区二区精品| 免费在线观看影片大全网站| 真人一进一出gif抽搐免费| 国产亚洲av高清不卡| 免费人成视频x8x8入口观看| 白带黄色成豆腐渣| 大型av网站在线播放| 成人永久免费在线观看视频| 制服丝袜大香蕉在线| 一区二区三区国产精品乱码| 欧美一区二区国产精品久久精品 | 99久久精品热视频| 亚洲国产精品久久男人天堂| 国产av一区二区精品久久| 特大巨黑吊av在线直播| 日韩 欧美 亚洲 中文字幕| 成人国语在线视频| 亚洲国产精品999在线| 老司机午夜十八禁免费视频| 男女午夜视频在线观看| 成人国语在线视频| 成人三级黄色视频| 久久久久国产一级毛片高清牌| 国产野战对白在线观看| 国产69精品久久久久777片 | 久久久国产精品麻豆| 国产精品一区二区三区四区久久| 欧美成人一区二区免费高清观看 | 日本在线视频免费播放| 一二三四在线观看免费中文在| 亚洲精品美女久久久久99蜜臀| 日本撒尿小便嘘嘘汇集6| 国产午夜精品论理片| 欧美一级a爱片免费观看看 | 日本免费a在线| √禁漫天堂资源中文www| 757午夜福利合集在线观看| 色综合亚洲欧美另类图片| 一边摸一边抽搐一进一小说| 精品国产乱子伦一区二区三区| 一级a爱片免费观看的视频| 精品无人区乱码1区二区| 国产乱人伦免费视频| 免费在线观看完整版高清| 免费看日本二区| 日日摸夜夜添夜夜添小说| 成人三级黄色视频| a级毛片在线看网站| 亚洲第一欧美日韩一区二区三区| 18美女黄网站色大片免费观看| 亚洲av成人精品一区久久| 国产精品影院久久| 欧美性猛交╳xxx乱大交人| 成人18禁高潮啪啪吃奶动态图| 国产精品综合久久久久久久免费| 日日干狠狠操夜夜爽| 蜜桃久久精品国产亚洲av| 日韩中文字幕欧美一区二区| 亚洲av日韩精品久久久久久密| 亚洲一区高清亚洲精品| 手机成人av网站| 亚洲国产欧洲综合997久久,| 人妻夜夜爽99麻豆av| 一个人免费在线观看电影 | 校园春色视频在线观看| 两个人的视频大全免费| 国产成人精品无人区| 麻豆一二三区av精品| 久久99热这里只有精品18| 精品无人区乱码1区二区| 又紧又爽又黄一区二区| 一个人免费在线观看的高清视频| 成人18禁高潮啪啪吃奶动态图| 亚洲国产精品999在线| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品一区二区精品视频观看| 国产一区二区三区在线臀色熟女| 国产精品永久免费网站| 天天躁夜夜躁狠狠躁躁| 国产亚洲av嫩草精品影院| 欧美日韩黄片免| 真人一进一出gif抽搐免费| 两个人的视频大全免费| 99热6这里只有精品| 日韩欧美国产一区二区入口| 国产单亲对白刺激| 草草在线视频免费看| √禁漫天堂资源中文www| 黄色女人牲交| 国产在线精品亚洲第一网站| 亚洲国产欧美人成| 欧美日韩国产亚洲二区| 国产高清有码在线观看视频 | 国产精品av久久久久免费| 91九色精品人成在线观看| 香蕉丝袜av| 成人高潮视频无遮挡免费网站| 免费人成视频x8x8入口观看| 亚洲中文字幕日韩| 亚洲人成电影免费在线| 人人妻,人人澡人人爽秒播| 首页视频小说图片口味搜索| 亚洲精品久久成人aⅴ小说| 伦理电影免费视频| 久久精品影院6| 精品久久久久久久人妻蜜臀av| 国产一区二区三区在线臀色熟女| 日本熟妇午夜| 亚洲精品久久国产高清桃花| 麻豆久久精品国产亚洲av| 女警被强在线播放| 免费看a级黄色片| 亚洲欧美精品综合久久99| 国产精品香港三级国产av潘金莲| 一本精品99久久精品77| 黄色视频不卡| 国产成+人综合+亚洲专区| 大型av网站在线播放| 三级男女做爰猛烈吃奶摸视频| 久久久久九九精品影院| 天天躁夜夜躁狠狠躁躁| 又爽又黄无遮挡网站| 男女下面进入的视频免费午夜| 又爽又黄无遮挡网站| 麻豆久久精品国产亚洲av| 亚洲aⅴ乱码一区二区在线播放 | 欧美又色又爽又黄视频| 午夜精品在线福利| 白带黄色成豆腐渣| av有码第一页| 黑人欧美特级aaaaaa片| 黑人巨大精品欧美一区二区mp4| 香蕉丝袜av| 精品一区二区三区av网在线观看| 亚洲国产精品999在线| 国产成人aa在线观看| 麻豆成人av在线观看| 丝袜美腿诱惑在线| 成人一区二区视频在线观看| 色综合亚洲欧美另类图片| 特大巨黑吊av在线直播| 色在线成人网| 99久久精品国产亚洲精品| 我的老师免费观看完整版| 桃红色精品国产亚洲av| 国产乱人伦免费视频| 国产黄色小视频在线观看| 9191精品国产免费久久| 国产一区二区在线av高清观看| 久久国产精品人妻蜜桃| 久久精品国产亚洲av香蕉五月| 精品欧美一区二区三区在线| 神马国产精品三级电影在线观看 | 身体一侧抽搐| 又大又爽又粗| 三级国产精品欧美在线观看 | 欧美高清成人免费视频www| 日本a在线网址| 国产亚洲av高清不卡| 精品熟女少妇八av免费久了| 国产在线观看jvid| 国产人伦9x9x在线观看| 欧美日韩中文字幕国产精品一区二区三区| 欧美精品啪啪一区二区三区| 青草久久国产| 久久欧美精品欧美久久欧美| 精品日产1卡2卡| 欧美午夜高清在线| 一区二区三区国产精品乱码| 欧美zozozo另类| 欧美三级亚洲精品| 欧美色欧美亚洲另类二区| 国产真人三级小视频在线观看| 国产又黄又爽又无遮挡在线| 欧美一区二区国产精品久久精品 | 精品久久久久久久毛片微露脸| 日本熟妇午夜| 给我免费播放毛片高清在线观看| 亚洲真实伦在线观看| 蜜桃久久精品国产亚洲av| www.999成人在线观看| 国产熟女xx| av超薄肉色丝袜交足视频| 妹子高潮喷水视频| 久久久久久亚洲精品国产蜜桃av| 久久99热这里只有精品18| 亚洲 欧美一区二区三区| 少妇裸体淫交视频免费看高清 | 午夜亚洲福利在线播放| 全区人妻精品视频| 美女 人体艺术 gogo| 国产真人三级小视频在线观看| 黑人欧美特级aaaaaa片| av在线天堂中文字幕| 国产精品1区2区在线观看.| 最近在线观看免费完整版| 男女午夜视频在线观看| 中文字幕最新亚洲高清| 国产一区二区三区在线臀色熟女| 亚洲精品一卡2卡三卡4卡5卡| 在线观看日韩欧美| 巨乳人妻的诱惑在线观看| 亚洲人成网站高清观看| 在线观看美女被高潮喷水网站 | 性色av乱码一区二区三区2| 一边摸一边做爽爽视频免费| 国产日本99.免费观看| 少妇被粗大的猛进出69影院| 免费一级毛片在线播放高清视频| 国产成人精品无人区| 麻豆久久精品国产亚洲av| 少妇的丰满在线观看| 岛国视频午夜一区免费看| 香蕉久久夜色| 久久性视频一级片| 日韩欧美精品v在线| 一本精品99久久精品77| 黄色视频,在线免费观看| 欧美成人一区二区免费高清观看 | 免费av毛片视频| 中文字幕高清在线视频| 啪啪无遮挡十八禁网站| 亚洲专区中文字幕在线| √禁漫天堂资源中文www| 久久久久久国产a免费观看| 手机成人av网站| 亚洲人成电影免费在线| 国产主播在线观看一区二区| 国产1区2区3区精品| 免费搜索国产男女视频| 大型黄色视频在线免费观看| 国产精品 欧美亚洲| videosex国产| 麻豆成人av在线观看| 久久久久性生活片| 18禁国产床啪视频网站| 久久精品国产综合久久久| www.www免费av| 人人妻人人澡欧美一区二区| 久久热在线av| 午夜亚洲福利在线播放| 我的老师免费观看完整版| 成人三级做爰电影| 天堂√8在线中文| 国产精品一区二区三区四区免费观看 | 欧美日本亚洲视频在线播放| 男女那种视频在线观看| 久久中文字幕一级| 欧美三级亚洲精品| 亚洲第一欧美日韩一区二区三区| 哪里可以看免费的av片| 亚洲国产精品合色在线| 一进一出抽搐动态| 一进一出好大好爽视频| 久久精品亚洲精品国产色婷小说| tocl精华| 男女下面进入的视频免费午夜| 看片在线看免费视频| 免费看十八禁软件| 色在线成人网| 国产探花在线观看一区二区| 欧美在线一区亚洲| 日韩精品青青久久久久久| 一本精品99久久精品77| 制服丝袜大香蕉在线| 国产99白浆流出| 成人三级做爰电影| videosex国产| 夜夜看夜夜爽夜夜摸| 亚洲第一欧美日韩一区二区三区| 黄色视频不卡| 在线a可以看的网站| 国产69精品久久久久777片 | 麻豆成人午夜福利视频| 国产主播在线观看一区二区| 欧美日本视频| 一本一本综合久久| 日韩 欧美 亚洲 中文字幕| 男男h啪啪无遮挡| 午夜福利在线观看吧| av中文乱码字幕在线| 久久热在线av| av视频在线观看入口| 日韩免费av在线播放| 国产人伦9x9x在线观看| 性色av乱码一区二区三区2| 老熟妇仑乱视频hdxx| 国产高清videossex| 国产黄色小视频在线观看| ponron亚洲| 中文在线观看免费www的网站 | 成人特级黄色片久久久久久久| 精品欧美一区二区三区在线| 欧美一区二区国产精品久久精品 | av视频在线观看入口| 国产三级中文精品| 国产aⅴ精品一区二区三区波| 国产高清有码在线观看视频 | 午夜精品在线福利| 久久亚洲精品不卡| 一本综合久久免费| 欧美日韩一级在线毛片| 他把我摸到了高潮在线观看| 99国产精品一区二区三区|