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

    長(zhǎng)江河口北支潮位與潮差的時(shí)空變化和機(jī)理

    2011-12-20 00:56:28宋永港朱建榮
    關(guān)鍵詞:潮差潮位青龍

    宋永港, 朱建榮, 吳 輝

    (華東師范大學(xué) 河口海岸學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海 200062)

    長(zhǎng)江河口北支潮位與潮差的時(shí)空變化和機(jī)理

    宋永港, 朱建榮, 吳 輝

    (華東師范大學(xué) 河口海岸學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海 200062)

    考慮長(zhǎng)江河口徑流、潮汐和風(fēng)場(chǎng)共同作用,數(shù)值模擬和定量分析北支潮位和潮差時(shí)空變化和動(dòng)力機(jī)制.北支月平均潮位呈現(xiàn)出從1月到7月逐漸增大,從8月到12月逐漸減小的變化趨勢(shì),主要決定于徑流量產(chǎn)生的余水位.潮差具有季節(jié)變化,一年中出現(xiàn)兩次極大值和兩次極小值.兩次極大值出現(xiàn)在3月(農(nóng)歷二月)和9月(農(nóng)歷八月),兩次極小值出現(xiàn)在6月(農(nóng)歷五月)和12月(農(nóng)歷十一月).青龍港最高潮位和最低潮位主要由潮汐強(qiáng)度決定,也與潮汐季節(jié)性變化和徑流量有關(guān).北支平均潮位、最高潮位和最低潮位縱向上基本呈現(xiàn)從下游往上游逐漸增加的變化規(guī)律,北支的潮差變化呈現(xiàn)出從口門向中段逐漸增大,再?gòu)闹卸蜗蛏隙沃饾u減小的規(guī)律.冬季偏北風(fēng)使青龍港潮位增大,夏季偏南風(fēng)使青龍港潮位和潮差略微減小,風(fēng)對(duì)北支潮差幾乎沒有影響.徑流產(chǎn)生的余水位增加潮位,徑流對(duì)潮汐具有抑制作用,使潮差減小.

    長(zhǎng)江河口; 北支; 潮位; 潮差; 影響因子

    0 引 言

    長(zhǎng)江河口三級(jí)分汊,第一級(jí)分為南支和北支,由崇明島分隔(見圖1).北支為喇叭形河道,從南北支分汊口崇頭至入海口連興港全長(zhǎng)約86 km.近60多年來(lái),由于大量圍墾促淤工程的實(shí)施,北支上段河道變淺,河寬縮窄,導(dǎo)致水道萎縮,北支洪、枯季的凈分流比逐漸減小至目前的5%以下[1].河勢(shì)的演變?cè)斐杀敝搅髯饔脺p弱、潮汐作用加強(qiáng),逐漸變成以潮流為主的潮汐通道,并進(jìn)一步影響到北支河道的沖淤和河床的演變[2,3].北支另外一個(gè)重要的現(xiàn)象就是出現(xiàn)了冬季北支鹽水倒灌進(jìn)入南支,潮差越大,倒灌越強(qiáng),倒灌強(qiáng)度隨潮差指數(shù)增長(zhǎng)[4].北支倒灌進(jìn)入南支的鹽水在徑流的作用下向下游移動(dòng),影響陳行水庫(kù)、寶鋼水庫(kù)和青草沙水庫(kù)取水,因此研究北支潮差的時(shí)空變化對(duì)上海供水安全具有重要的應(yīng)用意義.

    以往對(duì)長(zhǎng)江河口北支潮差的研究主要是基于實(shí)測(cè)資料從定性的角度分析北支潮汐和潮差的變化[5],沒有詳細(xì)研究時(shí)空變化,更沒有研究徑流和風(fēng)對(duì)它們的影響.本文應(yīng)用改進(jìn)后的ECOM-si數(shù)值模式模擬北支潮位和潮差,從定量上給出北支潮位潮差隨時(shí)間和空間的變化,并分析徑流和風(fēng)對(duì)潮位和潮差的影響.北支潮汐和潮差是基本的水文要素,揭示其時(shí)空變化規(guī)律和動(dòng)力機(jī)理,具有重要理論意義.

    圖1 長(zhǎng)江河口地形圖Fig.1 Map of the Changjiang Estuary

    1 數(shù)值模式設(shè)置和驗(yàn)證

    本文數(shù)值模式應(yīng)用改進(jìn)后的ECOM-si模式[6,7],該模式長(zhǎng)期應(yīng)用于長(zhǎng)江河口水動(dòng)力過(guò)程和鹽水入侵的模擬,已取得了眾多的成果[4-9].數(shù)值模式的計(jì)算區(qū)域包括整個(gè)長(zhǎng)江河口、杭州灣和鄰近海區(qū),上邊界在大通,外海邊界到東經(jīng)125°,南到27.5°N,北至33.5°N.外海開邊界以潮位驅(qū)動(dòng),考慮16個(gè)主要分潮 M2、S2、N2、K2、K1、O1、P1、Q1、U2、V2、T2、L2、2N2、J1、M1和OO1,由各分潮調(diào)和常數(shù)合成給出.上游大通開邊界條件以實(shí)測(cè)徑流量給出,海平面考慮風(fēng)應(yīng)力.初始條件水位和流速均取為零,初始溫度、鹽度取自實(shí)測(cè)資料結(jié)合海洋圖冊(cè)給出.長(zhǎng)江河口水深資料由各區(qū)最新測(cè)得的數(shù)據(jù)組合而成,北支上段、南支上段和南匯邊灘由2008年觀測(cè)水深資料給出,口門內(nèi)其它區(qū)域大部分由2007年觀測(cè)資料給出,口外大部分區(qū)域水深由海軍航保部海圖數(shù)值化給出.

    采用2006年12月份的蘆潮港實(shí)測(cè)潮位資料,以及2010年3月15日至4月4日北支青龍港的潮位資料(實(shí)測(cè)資料基面均已調(diào)整至平均海平面)驗(yàn)證數(shù)值模式.模型驗(yàn)證期間采用實(shí)測(cè)的風(fēng)場(chǎng)和大通站徑流量資料.從驗(yàn)證結(jié)果看,總體上模式計(jì)算值與實(shí)測(cè)值吻合良好(見圖2和圖3).蘆潮港站實(shí)測(cè)潮位與模型計(jì)算潮位的皮爾遜相關(guān)系數(shù)為0.888,青龍港站的皮爾遜相關(guān)系數(shù)為0.963;這說(shuō)明模型計(jì)算值與測(cè)量值的相關(guān)性好.另外,為了說(shuō)明潮周期中各個(gè)時(shí)段計(jì)算值與測(cè)量值的相關(guān)性,在每個(gè)時(shí)間點(diǎn)處前后各取12 h內(nèi)的數(shù)據(jù)進(jìn)行24 h內(nèi)的皮爾遜相關(guān)性分析,得到相關(guān)系數(shù)的時(shí)間序列(見圖2和圖3中的下圖).從圖2中可以看出蘆潮港站,模型計(jì)算的潮位在小潮時(shí)與實(shí)測(cè)值的相關(guān)性要高于大潮;青龍港站模型計(jì)算值與實(shí)測(cè)值的相關(guān)性是大潮好于小潮.這樣的誤差主要由地形、底摩擦、開邊界條件誤差和模式本身存在計(jì)算誤差等因素造成的.總的說(shuō)來(lái),ECOM-si模式能較好模擬長(zhǎng)江河口和杭州灣潮汐變化.

    圖2 2006年12月蘆潮港模式計(jì)算潮位(實(shí)線)與實(shí)測(cè)潮位(黑點(diǎn))(上圖)、皮爾遜相關(guān)系數(shù)(下圖)隨時(shí)間變化Fig.2 Comparison of the simulated tidal level(solid line)with the measured data(dots)(upper)and the Pearson correlation coefficient(lower)at the Luchaogang station in December 2006

    2 北支潮位和潮差時(shí)空變化

    本節(jié)考慮徑流、潮汐和風(fēng)場(chǎng)共同作用下北支潮位和潮差時(shí)空變化,計(jì)算方案稱為數(shù)值實(shí)驗(yàn)A.在時(shí)間上分析各月潮位和潮差的變化,在空間上分析北支上段、中段和下段3個(gè)代表性測(cè)站青龍港、三條港和連興港潮位和潮差的變化.

    各月模式上游邊界條件以大通站1950~2005年各月平均徑流量給出[10](見表1),徑流量季節(jié)變化明顯,12、1、2、3月徑流量分別為14 504、10 904、11 707、15 886 m3/s,為枯季徑流量;6、7、8、9月徑流量分別為40 638、50 759、44 596、40 171 m3/s,為洪季徑流量.風(fēng)場(chǎng)采用NOAA的NECP數(shù)據(jù),以多年半月平均的形式給出,平均的時(shí)間段從1999年到2008年,計(jì)算區(qū)域的半月平均風(fēng)場(chǎng)隨時(shí)間變化見圖4.圖中顯示9月至3月為偏北風(fēng),4月至8月為偏南風(fēng),季風(fēng)特征明顯.

    圖3 2010年3月15日至4月5日青龍港模式計(jì)算潮位(實(shí)線)與實(shí)測(cè)潮位(黑點(diǎn))(上圖)、皮爾遜相關(guān)系數(shù)(下圖)隨時(shí)間變化Fig.3 Comparison of the simulated tidal level(solid line)with the measured data(dots)(upper)and the Pearson correlation coefficient(lower)at the Qinglonggang station from March 15 to April 5 in 2010

    表1 大通多年月平均徑流量Tab.1 The monthly mean river discharge at the Datong station m3·s-1

    圖4 模式計(jì)算區(qū)域多年半月平均風(fēng)矢量Fig.4 The semimonthly mean wind vector in the model domain

    2.1 潮位和潮差時(shí)間變化

    本文分析青龍港潮位和潮差隨時(shí)間變化,對(duì)實(shí)驗(yàn)A的計(jì)算結(jié)果作統(tǒng)計(jì),得到青龍港12個(gè)月各月的最高潮位、最低潮位和平均潮位(見表2).可見青龍港月平均潮位呈現(xiàn)出從1月到7月逐漸增大,從8月到12月逐漸減小的變化趨勢(shì);平均潮位最大值出現(xiàn)在7、8月份,約0.74 m,平均潮位最小值出現(xiàn)在2月份,為0.36 m.平均潮位體現(xiàn)的是余水位,主要由徑流量決定.青龍港月平均潮位與長(zhǎng)江月平均徑流量呈現(xiàn)出明顯正相關(guān)性,即長(zhǎng)江徑流量增大時(shí),青龍港的平均潮位增大;徑流量減小時(shí),平均潮位減小.平均潮位與徑流量呈現(xiàn)出這種正相關(guān)性說(shuō)明徑流量是引起平均潮位變化的原因之一.

    表2 青龍港每月潮位和潮差統(tǒng)計(jì)(單位:m)Tab.2 The monthly statistics of the tidal level and tidal range at Qinglonggang

    長(zhǎng)江河口為半日潮,對(duì)青龍港每日相鄰的最高潮位和最低潮位相減,得到兩個(gè)潮差,取大的潮差作為當(dāng)日的潮汐強(qiáng)度,一年逐日潮差隨時(shí)間變化見圖5.徑流和風(fēng)的作用會(huì)產(chǎn)生余水位,同時(shí)體現(xiàn)在最高潮位和最低潮位中,兩者相減,基本消除了徑流和風(fēng)的作用,故潮差最能體現(xiàn)潮汐強(qiáng)度.從圖5可見,潮差具有季節(jié)和顯著的半月變化.半月變化即大小潮變化,一月中有兩次大潮和兩次小潮,大潮期間潮差大,小潮期間潮差小,但存在月不等現(xiàn)象.上半年1~6月,前半月最大潮差比后半月最大潮差大;下半年7~12月,后半月最大潮差比前半月最大潮差大.季節(jié)變化體現(xiàn)在一年中潮差出現(xiàn)兩次極大值和兩次極小值.兩次極大值出現(xiàn)在3月(農(nóng)歷二月)和9月(農(nóng)歷八月),量值分別為4.29 m和4.02 m(見表2);兩次極小值出現(xiàn)在6月(農(nóng)歷五月)和12月(農(nóng)歷十一月),量值分別為3.52 m和3.78 m.潮差季節(jié)變化呈現(xiàn)出這種極大潮差出現(xiàn)在3月和9月的原因在于3月和9月地球分別處在春分點(diǎn)和秋分點(diǎn)附近,此時(shí)的月球、地球和太陽(yáng)三者的位置最接近直線,太陽(yáng)的引潮力和月球的引潮力的合力達(dá)到一年中的極大值,形成一年中的大潮,潮差在一年中達(dá)到極大值.各月每日潮差中取最大值,作為該月潮汐的強(qiáng)度,量值見表2.

    潮差的這種季節(jié)變化對(duì)長(zhǎng)江河口具有重要的指導(dǎo)意義.潮差對(duì)長(zhǎng)江河口的鹽水入侵有著巨大的影響,潮差越大,倒灌越強(qiáng),倒灌強(qiáng)度隨潮差指數(shù)增長(zhǎng)[4].3月潮汐最強(qiáng),且2月和1月潮差在一年中強(qiáng)度為第二和第三,量值分別為4.22 m和4.29 m(見表2),加上1~3月長(zhǎng)江徑流量低,這種低徑流量和強(qiáng)潮差的組合,導(dǎo)致長(zhǎng)江河口冬季發(fā)生嚴(yán)重的鹽水入侵.9月盡管潮差出現(xiàn)一個(gè)極大值,但平均徑流量達(dá)到40 000 m3/s,故一般不會(huì)發(fā)生北支鹽水倒灌南支現(xiàn)象.三峽工程建成后,三峽水庫(kù)季節(jié)性調(diào)水,若9月大幅蓄水導(dǎo)致徑流量大減,因9月出現(xiàn)第二個(gè)潮差極大值,會(huì)導(dǎo)致提前和加劇長(zhǎng)江河口鹽水入侵,需引起注意.

    圖5 青龍港站模式計(jì)算逐日潮差隨時(shí)間變化Fig.5 The simulated daily tidal range at the Qinglonggang station

    最高潮位與平均潮位略有不同,最高潮位1~3月增大,再逐漸減小至6月,然后又逐漸增大直到8、9月,再逐漸減小,呈M形變化(見表2).最低潮位與平均潮位的變化趨勢(shì)類似,1月最低,量值-1.3 m,隨后增大,至8月達(dá)到最高值,量值-0.67 m,之后逐漸減小,至12月量值為-1.15 m.最高潮位和最低潮位主要由潮汐強(qiáng)度決定,與季節(jié)性潮汐變化有關(guān),也與余水位有關(guān),是由兩者共同決定的.

    2.2 潮位和潮差空間變化

    通過(guò)對(duì)比北支縱向上青龍港、三條港和連興港3個(gè)站點(diǎn)的潮位和潮差的變化,分析潮位和潮差的空間變化.表3給出了這3個(gè)站點(diǎn)各月最高潮位、最低潮位和平均潮位,可見除1~3月青龍港最高潮位比三條港低外,青龍港平均潮位、最高潮位和最低潮位均比三條港的高,而三條港的平均潮位、最高潮位和最低潮位又比連興港的高.因此,北支潮位縱向上基本呈現(xiàn)從下游往上游逐漸增加的變化規(guī)律.

    圖6給出了這3個(gè)站點(diǎn)潮差隨時(shí)間的變化,可以看到青龍港、三條港和連興港這3個(gè)站點(diǎn)中潮差最大的地方是三條港,其最大潮差達(dá)5.03 m,出現(xiàn)在3月,年平均潮差3.21 m;連興港的潮差其次,最大潮差4.83 m,年平均潮差3.05 m,比三條港略?。磺帻埜墼谶@3個(gè)站點(diǎn)中潮差最小,最大值4.29 m,年平均潮差只有2.64 m.這說(shuō)明北支的潮差變化呈現(xiàn)出從口門向中段逐漸增大,再?gòu)闹卸蜗蛏隙沃饾u減小的規(guī)律.這與沈煥庭等所描述的現(xiàn)象一樣[8],原因在于北支地形呈喇叭口,潮波進(jìn)入北支口門后,過(guò)水?dāng)嗝嫜杆偈湛s,潮波單寬能量增加,潮差增大;但隨著潮流繼續(xù)向內(nèi)推進(jìn),能量損耗越來(lái)越大,當(dāng)單寬能量增加量小于因損耗的能量時(shí),潮差會(huì)減小;北支上段由于徑流頂托作用和摩擦阻力作用導(dǎo)致的能量損耗超過(guò)能量聚集,潮差逐漸減?。?].

    3 風(fēng)和徑流量對(duì)北支潮位和潮差的影響

    本節(jié)設(shè)計(jì)兩個(gè)數(shù)值實(shí)驗(yàn)B和C,數(shù)值實(shí)驗(yàn)B考慮徑流和潮流作用,與數(shù)值實(shí)驗(yàn)A比較,分析風(fēng)對(duì)北支潮位和潮差的影響;數(shù)值實(shí)驗(yàn)C僅考慮潮流作用,與數(shù)值實(shí)驗(yàn)B比較,分析徑流對(duì)北支潮位和潮差的影響.

    表3 每月北支潮位縱向比較Tab.3 Longitudinal comparison of the monthly tidal level in the North Branch m

    圖6 北支三條港(實(shí)線)、連興港(虛線)、青龍港(黑粗線)每日潮差隨時(shí)間變化Fig.6 Temporal variation of the daily tidal range at the Santiaogang station(solid line),the Lianxinggang station(dash line)and the Qinglonggang station(solid bold line)in the North Branch

    3.1 風(fēng)對(duì)北支潮位和潮差的影響

    將數(shù)值實(shí)驗(yàn)A計(jì)算出的青龍港每日最大潮位和潮差減去數(shù)值實(shí)驗(yàn)B計(jì)算出的每日最大潮位和潮差,得到每日最大潮位差值和潮差差值隨時(shí)間變化,兩者差值可體現(xiàn)風(fēng)的作用.

    長(zhǎng)江河口冬季盛行偏北風(fēng),夏季盛行偏南風(fēng)(見圖4).總體上,冬季偏北風(fēng)使青龍港潮位增大,夏季偏南風(fēng)使青龍港潮位和潮差減小,但影響微?。ㄒ妶D7).北風(fēng)對(duì)潮位有增高作用是因?yàn)楸憋L(fēng)產(chǎn)生Ekman輸運(yùn),在北支口附近形成向西輸運(yùn)的風(fēng)生流[10],引起北支內(nèi)水位的抬升;而夏季南風(fēng)引起的Ekman輸運(yùn)方向指向口外,有降低水位的作用.另外,每日最大潮位差值存在著半月變化.原因在于北支水深淺,而平均潮差約在3 m,大潮期間水深的變幅遠(yuǎn)大于小潮期間水深的變幅,相同的風(fēng)場(chǎng)在大潮和小潮期間對(duì)每日最大潮位作用不同,即風(fēng)和潮汐之間存在著非線性相互作用.與圖5中青龍港潮差隨時(shí)間變化比較,大潮期間偏北風(fēng)對(duì)潮位增加不明顯,甚至降低;小潮期間偏北風(fēng)對(duì)潮位增加相對(duì)明顯.

    總體上風(fēng)對(duì)北支潮差幾乎沒有影響,但存在著更顯著的半月變化,且與潮差存在負(fù)相關(guān)(見圖8),原因同樣在于風(fēng)與潮汐的非線性相互作用.大潮期間偏北風(fēng)減弱潮差,小潮期間偏北風(fēng)增大潮差,量值在0.02~0.04 m之間.

    圖7 青龍港每日最大潮位差值(數(shù)值實(shí)驗(yàn)A-數(shù)值實(shí)驗(yàn)B)隨時(shí)間變化Fig.7 The difference of the daily highest tidal level between the numerical experiment A and B(A-B)at the Qinglonggang staion

    圖8 青龍港每日潮差差值(數(shù)值實(shí)驗(yàn)A-數(shù)值實(shí)驗(yàn)B)隨時(shí)間變化Fig.8 The difference of the daily tidal range between the numerical experiment A and B(A-B)at the Qinglonggang station

    3.2 徑流對(duì)北支潮位和潮差的影響

    將數(shù)值實(shí)驗(yàn)B計(jì)算出的青龍港每日最大潮位和潮差減去數(shù)值實(shí)驗(yàn)C計(jì)算出的每日最大潮位和潮差,得到最大潮位差值和潮差差值隨時(shí)間變化,兩者差值可體現(xiàn)徑流的作用.

    從物理概念上看,減去潮汐作用,最大潮位差值基本上為徑流引起的余水位,量值與徑流量有關(guān).徑流總會(huì)產(chǎn)生余水位,故潮位差值為正,1~7月逐漸增大,8月后又逐漸減小,變化趨勢(shì)與徑流的變化趨勢(shì)一致(見圖9).但潮位差值存在著極為顯著的半月變化,原因在于徑流與潮汐存在非線性相互作用,與圖5比較,可見大潮期間徑流引起的青龍港潮位增加大,小潮期間小,徑流與潮差相互作用對(duì)潮位為正相關(guān)關(guān)系.

    徑流對(duì)潮汐具有抑制作用,使潮差減小.徑流量越大,減弱潮差的作用越明顯(見圖10).同樣存在半月變化,與圖5比較,大潮期間徑流對(duì)潮差的抑制作用弱,小潮期間強(qiáng).

    圖9 青龍港每日最大潮位差值(數(shù)值實(shí)驗(yàn)B-數(shù)值實(shí)驗(yàn)C)隨時(shí)間變化Fig.9 The difference of the daily highest tidal level between the numerical experiment B and C(B-C)at the Qinglonggang staion

    圖10 青龍港潮差差值(數(shù)值實(shí)驗(yàn)B-數(shù)值實(shí)驗(yàn)C)隨時(shí)間變化Fig.10 The difference of the daily highest tidal level between the numerical experiment B and C(B-C)at the Qinglonggang station

    4 結(jié) 論

    本文考慮長(zhǎng)江河口徑流、潮汐和風(fēng)場(chǎng)共同作用,數(shù)值模擬和定量分析北支潮位和潮差時(shí)空變化和動(dòng)力機(jī)制,主要結(jié)論如下.

    (1)北支青龍港站的月平均潮位呈現(xiàn)出1~7月逐漸增大,8~12月逐漸減小的變化規(guī)律,這種變化規(guī)律主要取決于北支徑流引起的余水位的變化.潮差最能體現(xiàn)潮汐強(qiáng)度,存在月不等現(xiàn)象,具有季節(jié)和顯著的半月變化;潮差季節(jié)變化體現(xiàn)在一年中潮差出現(xiàn)兩次極大值和兩次極小值;兩次極大值出現(xiàn)在3月(農(nóng)歷二月)和9月(農(nóng)歷八月),兩次極小值出現(xiàn)在6月(農(nóng)歷五月)和12月(農(nóng)歷十一月);青龍港最高潮位和最低潮位主要由潮汐強(qiáng)度決定,也與潮汐季節(jié)性變化和徑流量有關(guān).

    (2)北支潮位和潮差隨空間變化體現(xiàn)在北支平均潮位、最高潮位和最低潮位縱向上基本呈現(xiàn)從下游往上游逐漸增加的變化規(guī)律,北支的潮差變化呈現(xiàn)出從口門向中段逐漸增大,再?gòu)闹卸蜗蛏隙沃饾u減小的規(guī)律.

    (3)風(fēng)對(duì)北支潮位和潮差的影響體現(xiàn)在冬季偏北風(fēng)使青龍港潮位增大,夏季偏南風(fēng)使青龍港潮位和潮差減小,但影響微小.風(fēng)和潮汐之間存在著非線性相互作用,大潮期間偏北風(fēng)對(duì)潮位增加不明顯,小潮期間偏北風(fēng)對(duì)潮位增加相對(duì)明顯;總體上風(fēng)對(duì)北支潮差幾乎沒有影響.

    (4)徑流對(duì)北支潮位和潮差的影響體現(xiàn)在徑流產(chǎn)生的余水位增加潮位,其變化趨勢(shì)與徑流的變化趨勢(shì)一致;大潮期間徑流引起的青龍港潮位增加大,小潮期間?。粡搅鲗?duì)潮汐具有抑制作用,使潮差減小;大潮期間徑流對(duì)潮差的抑制作用弱,小潮期間強(qiáng).

    [1] 李伯昌.1984年以來(lái)長(zhǎng)江口北支演變分析[J].水利水運(yùn)工程學(xué)報(bào),2006(9):9-17.

    LI B C.Channel evolution in North Branch of Changjiang River estuary since 1984[J].Hydro-Science and Engineering,2006(9):9-17.

    [2] 陳沈良,陳吉余,谷國(guó)傳.長(zhǎng)江口北支的涌潮及其對(duì)河口的影響[J].華東師范大學(xué)學(xué)報(bào):自然科學(xué)版,2003(2):74-80.

    CHEN S L,CHRN J Y,GU G C.The tidal bore on the North Branch of Changjiang Estuary and its effects on the estuary[J].Journal of Eastchina Normal University:Natural Science,2003(2):74-80.

    [3] 張靜怡,胡震云,黃志良.近年長(zhǎng)江口北支涌潮變化及其成因分析[J].水科學(xué)進(jìn)展,2007(9):724-729.

    ZHANG J Y,HU Z Y,HUANG Z L.Analysis of the changes and formation causes of tidal bore at the north bayou of the Yangtze River estuary in recent years[J].Advances in Water Science,2007(9):724-729.

    [4] WU H,ZHU J R,CHEN B R,et al.Quantitative relationship of runoff and tide to saltwater spilling over from the North Branch in the Changjiang Estuary:A numerical study[J].Estuarine,Coastal and Shelf Science,2006,69:125-132.

    [5] 陳吉余,沈煥庭,惲才興.長(zhǎng)江河口動(dòng)力地貌和地貌演變[M].上海:上??萍汲霭嫔?,1988:76-79.

    CHEN J Y,SHEN H T,YUN C X.Process of Dynamic and Geomorphology of the Changjiang Estuary[M].Shanghai:Shanghai Scientific and Technical Publishers,1988:76-79.

    [6] 朱建榮.海洋數(shù)值計(jì)算方法和數(shù)值模式[M].北京:海洋出版社,2003.

    ZHU J R.Ocean Numerical Calculation Methods and Numerical Model[M].Beijing:China Ocean Press,2003.

    [7] 朱建榮 朱首賢.ECOM模式的改進(jìn)及在長(zhǎng)江河口、杭州灣及鄰近海區(qū)的應(yīng)用[J].海洋與湖沼,2003(7):364-374.

    ZHU J R,ZHU S X.Improvement of the ecom with application to the changjiang river estuary,hangzhou bay and adjacent waters[J].Oceanologia et Limnologia Sinica,2003(7):364-374.

    [8] 沈煥庭,茅志昌,朱建榮.長(zhǎng)江河口鹽水入侵[M].北京:海洋出版社,2003.

    SHEN H T,MAO Z C,ZHU J R.Saltwater Intrusion in the Changjiang Estuary[M].Beijing:China Ocean Press,2003.

    [9] LI L,ZHU J R,WU H,et al.A numerical study on the water diversion ratio of the Changjiang Estuary during the dry season[J].Chinese Journal of Oceanology and Limnology,2001,28(3):700-712.

    [10] 陳義忠.黃海東海環(huán)流和長(zhǎng)江沖淡水季節(jié)連續(xù)變化的數(shù)值模擬[D].上海:華東師范大學(xué),2007:61.

    CHEN Y Z.Numerical simulation of seasonal changes of the circulation and changjiang diluted water in the Yellow Sea and East China Sea[D].Shanghai:East China Normal University,2007:61.

    Spatial and temporal variations and mechanism of the tidal level and range in the North Branch of the Changjiang Estuary

    SONG Yong-gang, ZHU Jian-rong, WU Hui

    (State Key Laboratory of Estuarine and Coastal Research,East China Normal University,Shanghai 200062,China)

    With the combined effects of the river discharge,tide and wind stress,the spatial/temporal variations as well as the dynamical mechanism of the tidal level and tidal range in the North Branch of the Changjiang Estuary were numerically simulated and quantitatively analyzed in this paper.The monthly mean tidal level in the North Branch is increased gradually from January to July and decreased gradually from October to December.It mainly depends on the residual water level,which is determined by the amount of the river discharge.The tidal range experiences a seasonal variation,with two maximum and two minimum values in a year.These two maximum values appear in March (the second lunar month)and September(the eighth lunar month),and the two minimum ones appear in June(the fifth lunar month)and December(the eleventh lunar month),respectively.The highest and lowest tidal levels at the Qinglonggang sta-tion are mainly determined by the intensity of tidal forcing,and associate with the seasonal variations of the tide and river discharge.The averaged,highest and lowest tidal levels represent a landward increment feature along the North Branch.The tidal range is increased gradually from the mouth to the middle reaches of the North Branch,and decreased from its middle reaches to its upper reaches.The northerly wind in winter tends to produce an increment of the tidal level at the Qinglonggang,and the southerly wind in summer tends to produce a slight decrement of the tidal level and range there.The wind stress has almost no effect on the tidal range in the North Branch.The river discharge increases the tidal level,restrains the tide,and decreased the tidal range.

    the Changjiang Estuary; the North Branch; tidal level; tidal range; impact factors

    P731.2

    A

    10.3969/j.issn.1000-5641.2011.06.002

    1000-5641(2011)06-0010-10

    2010-09

    國(guó)家自然科學(xué)基金(40976056,40806034);國(guó)家水體污染控制與治理科技重大專項(xiàng)(2008ZX07421-001)

    宋永港,男,碩士研究生.E-mail:syg1530216@163.com.

    朱建榮,男,教授,博士生導(dǎo)師,從事河口海洋學(xué)研究.E-mail:jrzhu@sklec.ecnu.edu.cn.

    猜你喜歡
    潮差潮位青龍
    基于距離倒數(shù)加權(quán)的多站潮位改正方法可行性分析
    唐山市警戒潮位標(biāo)志物維護(hù)研究
    長(zhǎng)江口潮差中長(zhǎng)期變化對(duì)河口生態(tài)環(huán)境的影響
    多潮位站海道地形測(cè)量潮位控制方法研究
    少林功夫拳(三)
    少林與太極(2018年9期)2018-09-28 08:37:34
    基于Copula函數(shù)的甬江流域設(shè)計(jì)潮位過(guò)程研究
    小青龍說(shuō)“角”
    基于改進(jìn)的OLS-RBF模型的感潮河段潮位預(yù)測(cè)研究
    青龍現(xiàn)身記
    奧秘(2016年8期)2016-09-06 16:58:26
    論地質(zhì)鉆探技術(shù)在海水域工程的方法措施
    地球(2015年8期)2015-10-21 19:30:08
    色噜噜av男人的天堂激情| 久久这里只有精品中国| a级毛片免费高清观看在线播放| 久久天躁狠狠躁夜夜2o2o| 夜夜夜夜夜久久久久| 99久久精品一区二区三区| 国产精品亚洲一级av第二区| 久久精品久久久久久噜噜老黄 | 亚洲成人久久性| 无遮挡黄片免费观看| 久久久成人免费电影| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲七黄色美女视频| АⅤ资源中文在线天堂| 中国美女看黄片| 99热全是精品| 熟妇人妻久久中文字幕3abv| 日本五十路高清| 男女边吃奶边做爰视频| 国产中年淑女户外野战色| 青春草视频在线免费观看| 插逼视频在线观看| 五月玫瑰六月丁香| 国产精品久久视频播放| 日本一本二区三区精品| 一区二区三区高清视频在线| 久久久精品94久久精品| 一级黄片播放器| 一级毛片aaaaaa免费看小| 精品一区二区三区视频在线观看免费| 国产精品国产三级国产av玫瑰| 国产成人a∨麻豆精品| 特大巨黑吊av在线直播| 网址你懂的国产日韩在线| 久久国产乱子免费精品| 精品免费久久久久久久清纯| 国产91av在线免费观看| 小说图片视频综合网站| 国产成人91sexporn| 一个人观看的视频www高清免费观看| 国产视频内射| av黄色大香蕉| 国产单亲对白刺激| 亚洲一区高清亚洲精品| 国产精品伦人一区二区| 国产精品免费一区二区三区在线| 国产高清视频在线观看网站| 色综合站精品国产| 一夜夜www| 男人的好看免费观看在线视频| 久99久视频精品免费| 成年女人毛片免费观看观看9| 亚洲国产精品国产精品| 国产在视频线在精品| 国产精品伦人一区二区| 成熟少妇高潮喷水视频| 中文字幕免费在线视频6| 国产又黄又爽又无遮挡在线| 国产亚洲精品av在线| 亚洲无线观看免费| 久久精品国产亚洲av天美| 亚洲欧美精品综合久久99| 女人被狂操c到高潮| 老师上课跳d突然被开到最大视频| 18+在线观看网站| 欧美+日韩+精品| 长腿黑丝高跟| 亚洲国产精品成人综合色| 成人性生交大片免费视频hd| 尤物成人国产欧美一区二区三区| 赤兔流量卡办理| 国语自产精品视频在线第100页| 日韩av不卡免费在线播放| 一a级毛片在线观看| av在线天堂中文字幕| 亚洲最大成人手机在线| 国内久久婷婷六月综合欲色啪| 婷婷亚洲欧美| 自拍偷自拍亚洲精品老妇| АⅤ资源中文在线天堂| 成人无遮挡网站| 免费看av在线观看网站| 亚洲天堂国产精品一区在线| 精品久久久久久久久久免费视频| 中国美白少妇内射xxxbb| 联通29元200g的流量卡| 午夜视频国产福利| 久久久久国产网址| 欧美又色又爽又黄视频| 色播亚洲综合网| 亚洲久久久久久中文字幕| 国语自产精品视频在线第100页| www.色视频.com| 久久久久九九精品影院| 赤兔流量卡办理| 国产男靠女视频免费网站| 听说在线观看完整版免费高清| 成年免费大片在线观看| 波野结衣二区三区在线| 亚洲国产欧洲综合997久久,| 成年免费大片在线观看| 国产精品一区二区三区四区免费观看 | 国产不卡一卡二| 综合色丁香网| 老女人水多毛片| 永久网站在线| 嫩草影院精品99| 国产视频内射| 小说图片视频综合网站| 亚洲无线在线观看| 99热6这里只有精品| 国产伦一二天堂av在线观看| 日本 av在线| 男女视频在线观看网站免费| 精品国内亚洲2022精品成人| 亚洲欧美日韩无卡精品| 18禁裸乳无遮挡免费网站照片| 国产乱人偷精品视频| 中文亚洲av片在线观看爽| 国产成人91sexporn| 日本爱情动作片www.在线观看 | 97在线视频观看| 观看美女的网站| 国产真实乱freesex| 国产精品久久电影中文字幕| av在线亚洲专区| 成年女人看的毛片在线观看| 少妇丰满av| 超碰av人人做人人爽久久| 一本一本综合久久| 级片在线观看| 国产免费男女视频| 亚洲一区二区三区色噜噜| 亚洲激情五月婷婷啪啪| 亚洲成人精品中文字幕电影| 97人妻精品一区二区三区麻豆| 国产精品一及| or卡值多少钱| 插阴视频在线观看视频| 国产亚洲精品综合一区在线观看| 一级a爱片免费观看的视频| 亚洲精品影视一区二区三区av| 国产 一区 欧美 日韩| 一进一出抽搐gif免费好疼| 我要看日韩黄色一级片| 亚洲欧美精品综合久久99| 国产高清不卡午夜福利| 日韩强制内射视频| 老女人水多毛片| 91午夜精品亚洲一区二区三区| 校园人妻丝袜中文字幕| 国产 一区 欧美 日韩| 22中文网久久字幕| 欧美成人a在线观看| 亚洲人成网站高清观看| 精品人妻视频免费看| 亚洲高清免费不卡视频| 搡老熟女国产l中国老女人| 少妇被粗大猛烈的视频| 狠狠狠狠99中文字幕| 乱系列少妇在线播放| 中文字幕人妻熟人妻熟丝袜美| 国产精品人妻久久久久久| 日本黄色视频三级网站网址| 久久精品国产自在天天线| 国产精品女同一区二区软件| 欧美3d第一页| 欧美极品一区二区三区四区| a级一级毛片免费在线观看| 亚洲国产精品合色在线| av在线蜜桃| 久久久久久久久大av| 在线免费十八禁| 又黄又爽又刺激的免费视频.| 久久精品国产亚洲av涩爱 | 亚洲国产日韩欧美精品在线观看| 人人妻人人澡人人爽人人夜夜 | 99九九线精品视频在线观看视频| 精品午夜福利视频在线观看一区| 成人特级黄色片久久久久久久| 国产一区二区三区在线臀色熟女| 免费黄网站久久成人精品| 国产一级毛片七仙女欲春2| 狂野欧美激情性xxxx在线观看| 欧美成人免费av一区二区三区| 午夜福利在线观看免费完整高清在 | 露出奶头的视频| 成人亚洲精品av一区二区| 夜夜看夜夜爽夜夜摸| 中文字幕熟女人妻在线| 国产aⅴ精品一区二区三区波| 久久欧美精品欧美久久欧美| 午夜免费激情av| 午夜免费男女啪啪视频观看 | 精品久久久噜噜| 日日啪夜夜撸| 日韩亚洲欧美综合| 亚洲精品一卡2卡三卡4卡5卡| 久久精品影院6| 夜夜夜夜夜久久久久| 在线观看免费视频日本深夜| 亚洲av一区综合| 久久久a久久爽久久v久久| 亚洲精品一区av在线观看| 欧美xxxx黑人xx丫x性爽| 看片在线看免费视频| 一a级毛片在线观看| 在线免费十八禁| 久久久久国产精品人妻aⅴ院| 国产高清视频在线观看网站| 免费黄网站久久成人精品| 日本免费一区二区三区高清不卡| 卡戴珊不雅视频在线播放| 欧美最黄视频在线播放免费| av国产免费在线观看| 欧美色视频一区免费| 日韩,欧美,国产一区二区三区 | 成人国产麻豆网| 亚洲国产精品成人久久小说 | 18禁黄网站禁片免费观看直播| 99久国产av精品国产电影| 国产欧美日韩精品一区二区| www.色视频.com| a级毛色黄片| 国产熟女欧美一区二区| 亚洲av一区综合| 欧美成人免费av一区二区三区| 久久久久久久久久成人| 波野结衣二区三区在线| 国产蜜桃级精品一区二区三区| 一区二区三区免费毛片| 精品人妻一区二区三区麻豆 | 亚洲精品一区av在线观看| 偷拍熟女少妇极品色| 免费看av在线观看网站| 日韩成人伦理影院| 日本免费一区二区三区高清不卡| 国产一区二区亚洲精品在线观看| 亚洲欧美日韩高清专用| 国产av麻豆久久久久久久| 男女啪啪激烈高潮av片| 成人漫画全彩无遮挡| 日本撒尿小便嘘嘘汇集6| 欧美日本视频| 久久人人精品亚洲av| 一个人免费在线观看电影| 极品教师在线视频| 看十八女毛片水多多多| 国产高清激情床上av| av卡一久久| 国产在视频线在精品| 你懂的网址亚洲精品在线观看 | 国产aⅴ精品一区二区三区波| 国产视频一区二区在线看| 黑人高潮一二区| 久久久国产成人免费| 在线播放国产精品三级| 国产老妇女一区| 香蕉av资源在线| 成人无遮挡网站| 人妻丰满熟妇av一区二区三区| 我的老师免费观看完整版| 亚洲国产欧美人成| 能在线免费观看的黄片| 日韩一区二区视频免费看| 99国产极品粉嫩在线观看| 日本a在线网址| 国产久久久一区二区三区| 99视频精品全部免费 在线| 天堂动漫精品| 国产黄色小视频在线观看| 成年av动漫网址| 熟女电影av网| 色哟哟·www| 99热全是精品| 国产伦在线观看视频一区| 最新中文字幕久久久久| 99riav亚洲国产免费| 啦啦啦观看免费观看视频高清| 男女啪啪激烈高潮av片| 日本免费a在线| 国内久久婷婷六月综合欲色啪| 日本黄大片高清| 丰满乱子伦码专区| 午夜激情欧美在线| 久久久久免费精品人妻一区二区| 亚洲精品成人久久久久久| 69av精品久久久久久| 国产在视频线在精品| 人妻夜夜爽99麻豆av| 亚洲精品国产成人久久av| 最近的中文字幕免费完整| 成人一区二区视频在线观看| 久久精品人妻少妇| а√天堂www在线а√下载| 国产精品人妻久久久影院| 成人性生交大片免费视频hd| 国产亚洲精品av在线| 少妇被粗大猛烈的视频| 午夜免费激情av| 看黄色毛片网站| 少妇丰满av| 男人和女人高潮做爰伦理| 18+在线观看网站| 国内揄拍国产精品人妻在线| 男人和女人高潮做爰伦理| 国产精品人妻久久久久久| 91av网一区二区| 午夜免费激情av| 免费av观看视频| eeuss影院久久| 女人十人毛片免费观看3o分钟| av在线老鸭窝| 日日啪夜夜撸| 性插视频无遮挡在线免费观看| 亚洲在线观看片| 国产男靠女视频免费网站| 亚洲精品影视一区二区三区av| 青春草视频在线免费观看| 欧美日韩在线观看h| 亚洲精品影视一区二区三区av| 天堂网av新在线| 村上凉子中文字幕在线| 国模一区二区三区四区视频| 成人综合一区亚洲| 国产男靠女视频免费网站| 亚洲高清免费不卡视频| 亚洲丝袜综合中文字幕| 97人妻精品一区二区三区麻豆| 高清毛片免费观看视频网站| 看黄色毛片网站| 日本 av在线| 成人av在线播放网站| 91狼人影院| 免费av毛片视频| 久久人人爽人人片av| 日韩精品中文字幕看吧| 欧美精品国产亚洲| 成年免费大片在线观看| 欧美最新免费一区二区三区| 国产精品三级大全| 精品日产1卡2卡| 国产精品一区二区性色av| АⅤ资源中文在线天堂| 久久精品国产亚洲网站| 人人妻人人澡人人爽人人夜夜 | 18禁黄网站禁片免费观看直播| av在线播放精品| 美女被艹到高潮喷水动态| 国产午夜精品久久久久久一区二区三区 | 国产精品爽爽va在线观看网站| 女人被狂操c到高潮| 三级国产精品欧美在线观看| 欧美性感艳星| 观看免费一级毛片| 成人三级黄色视频| 精品午夜福利视频在线观看一区| 男女下面进入的视频免费午夜| 国产精品久久久久久久久免| 婷婷色综合大香蕉| 国产91av在线免费观看| 精品一区二区免费观看| 十八禁国产超污无遮挡网站| 免费高清视频大片| 国产欧美日韩精品一区二区| 91久久精品国产一区二区三区| 午夜视频国产福利| av在线天堂中文字幕| 久久亚洲国产成人精品v| 最近视频中文字幕2019在线8| 午夜视频国产福利| 我要看日韩黄色一级片| 欧美+日韩+精品| 日韩欧美免费精品| 色5月婷婷丁香| 99久久中文字幕三级久久日本| 最近在线观看免费完整版| 精品少妇黑人巨大在线播放 | 亚洲av二区三区四区| 女生性感内裤真人,穿戴方法视频| 亚洲国产精品久久男人天堂| 身体一侧抽搐| 国产精品国产三级国产av玫瑰| 国产精品国产三级国产av玫瑰| 少妇熟女欧美另类| 亚洲中文字幕日韩| 午夜老司机福利剧场| 国产探花极品一区二区| 免费av毛片视频| 能在线免费观看的黄片| 不卡视频在线观看欧美| АⅤ资源中文在线天堂| 老司机福利观看| 久久久午夜欧美精品| 禁无遮挡网站| 精品少妇黑人巨大在线播放 | 中文字幕久久专区| 欧美bdsm另类| 国产成人aa在线观看| 在线观看免费视频日本深夜| 91久久精品国产一区二区三区| 淫妇啪啪啪对白视频| 草草在线视频免费看| 一进一出抽搐动态| 你懂的网址亚洲精品在线观看 | 国产精品无大码| 亚洲最大成人手机在线| 欧美不卡视频在线免费观看| 亚洲熟妇熟女久久| 国产在视频线在精品| 国产伦精品一区二区三区视频9| 亚洲图色成人| 国产精品国产三级国产av玫瑰| av在线老鸭窝| 一级毛片电影观看 | 俄罗斯特黄特色一大片| 黄色一级大片看看| 日韩欧美精品免费久久| 12—13女人毛片做爰片一| 精品人妻熟女av久视频| 日本免费一区二区三区高清不卡| 久久人人爽人人片av| 国产高清视频在线播放一区| 少妇熟女欧美另类| 黄片wwwwww| 尤物成人国产欧美一区二区三区| 国产一区亚洲一区在线观看| 99热精品在线国产| 国产高清视频在线观看网站| 日本一本二区三区精品| 久久精品影院6| 日韩欧美在线乱码| 岛国在线免费视频观看| 国产精品久久久久久久电影| 精品人妻熟女av久视频| 国产真实伦视频高清在线观看| 成年av动漫网址| 国产精品,欧美在线| 99精品在免费线老司机午夜| 最近视频中文字幕2019在线8| 18+在线观看网站| 成年女人永久免费观看视频| 尾随美女入室| 欧美性感艳星| av福利片在线观看| 一区二区三区免费毛片| 成人特级黄色片久久久久久久| 精品一区二区三区视频在线| 大又大粗又爽又黄少妇毛片口| 国产精品久久电影中文字幕| 人人妻人人澡欧美一区二区| 男人的好看免费观看在线视频| 白带黄色成豆腐渣| or卡值多少钱| 国产在线男女| 长腿黑丝高跟| 亚洲av成人精品一区久久| 日韩欧美三级三区| 久久久久久久久久久丰满| 成人三级黄色视频| 欧美3d第一页| 蜜臀久久99精品久久宅男| av女优亚洲男人天堂| 国产老妇女一区| 亚洲七黄色美女视频| 成人永久免费在线观看视频| 亚洲无线在线观看| 天堂网av新在线| 久久久久久久久大av| 又爽又黄无遮挡网站| 22中文网久久字幕| 日本欧美国产在线视频| 成人亚洲欧美一区二区av| 色播亚洲综合网| 国内精品一区二区在线观看| 国产一区二区三区在线臀色熟女| 国产av麻豆久久久久久久| 亚洲av美国av| 国产伦精品一区二区三区四那| 国产一区二区三区在线臀色熟女| 三级经典国产精品| 成人av在线播放网站| 少妇人妻一区二区三区视频| 午夜久久久久精精品| 国产黄a三级三级三级人| 色5月婷婷丁香| 亚洲五月天丁香| 我的老师免费观看完整版| 美女高潮的动态| 日韩欧美一区二区三区在线观看| 国产精品电影一区二区三区| а√天堂www在线а√下载| 国产熟女欧美一区二区| 在线看三级毛片| 久久婷婷人人爽人人干人人爱| 99热精品在线国产| 日日干狠狠操夜夜爽| 小说图片视频综合网站| 非洲黑人性xxxx精品又粗又长| 3wmmmm亚洲av在线观看| 国内精品久久久久精免费| 嫩草影院精品99| 人人妻人人澡人人爽人人夜夜 | 亚洲在线自拍视频| 我要搜黄色片| av.在线天堂| 在线看三级毛片| 精品一区二区免费观看| 久久久久久久久中文| 夜夜爽天天搞| 日韩一本色道免费dvd| 午夜福利成人在线免费观看| 免费高清视频大片| 成人高潮视频无遮挡免费网站| 午夜福利视频1000在线观看| 亚洲真实伦在线观看| 国产乱人偷精品视频| 亚洲精品国产成人久久av| 一个人看的www免费观看视频| 国产精品一区www在线观看| 亚洲国产精品sss在线观看| 国产 一区 欧美 日韩| 国产精品精品国产色婷婷| 国产老妇女一区| 免费在线观看成人毛片| 亚洲av.av天堂| 日韩欧美在线乱码| 国模一区二区三区四区视频| 波多野结衣高清无吗| 免费黄网站久久成人精品| 国产一区二区在线av高清观看| 又黄又爽又免费观看的视频| 别揉我奶头~嗯~啊~动态视频| 日韩中字成人| 人妻丰满熟妇av一区二区三区| 日本爱情动作片www.在线观看 | 美女高潮的动态| 精品乱码久久久久久99久播| 特级一级黄色大片| 亚洲欧美日韩高清专用| 黄片wwwwww| 国产欧美日韩精品一区二区| 毛片一级片免费看久久久久| 日韩人妻高清精品专区| 久久精品国产亚洲av天美| 成人特级av手机在线观看| 欧美不卡视频在线免费观看| 亚洲人与动物交配视频| 色综合亚洲欧美另类图片| 欧美高清成人免费视频www| 欧美成人一区二区免费高清观看| 国产乱人视频| 欧美国产日韩亚洲一区| 淫秽高清视频在线观看| 九九在线视频观看精品| 免费看av在线观看网站| 久久久久久久久大av| 亚洲无线观看免费| 亚洲四区av| 亚洲精品日韩av片在线观看| 精品久久久噜噜| 搡老妇女老女人老熟妇| 亚洲国产高清在线一区二区三| 97超级碰碰碰精品色视频在线观看| 国产精品久久久久久精品电影| 在线播放无遮挡| 99久久精品国产国产毛片| 2021天堂中文幕一二区在线观| 国产久久久一区二区三区| 国产亚洲91精品色在线| 人妻制服诱惑在线中文字幕| 最近在线观看免费完整版| 最好的美女福利视频网| 精品久久久噜噜| 色播亚洲综合网| 此物有八面人人有两片| 久久久久精品国产欧美久久久| 日韩三级伦理在线观看| 麻豆一二三区av精品| 亚洲av一区综合| 欧美日韩国产亚洲二区| 日韩一本色道免费dvd| 欧美3d第一页| 最新中文字幕久久久久| 国产色爽女视频免费观看| 床上黄色一级片| 国产精品亚洲美女久久久| 色5月婷婷丁香| a级毛片a级免费在线| 看免费成人av毛片| 国产亚洲精品久久久com| 国产男靠女视频免费网站| 亚洲人成网站在线播| 免费不卡的大黄色大毛片视频在线观看 | 成人性生交大片免费视频hd| 免费看a级黄色片| 精品久久久久久成人av| 亚洲精品色激情综合| 真人做人爱边吃奶动态| 久久国内精品自在自线图片| 亚洲成av人片在线播放无| 淫秽高清视频在线观看| 一个人看视频在线观看www免费| 观看美女的网站| 国产高清不卡午夜福利| 狂野欧美白嫩少妇大欣赏| avwww免费| 99精品在免费线老司机午夜| 成人午夜高清在线视频| 成人特级黄色片久久久久久久| 99热全是精品| 自拍偷自拍亚洲精品老妇| 午夜精品国产一区二区电影 | 人妻久久中文字幕网| 国产私拍福利视频在线观看|