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

    定常風(fēng)對(duì)鄱陽湖水動(dòng)力的影響*

    2016-02-23 10:21:08國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃973項(xiàng)目2012CB417003江西省重大生態(tài)安全問題監(jiān)控協(xié)同創(chuàng)新中心項(xiàng)目JXSEW00和國家自然科學(xué)基金項(xiàng)目4130102341371062聯(lián)合資助20141205收稿20150418
    湖泊科學(xué) 2016年1期
    關(guān)鍵詞:鄱陽湖

    國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃“973”項(xiàng)目(2012CB417003)、江西省重大生態(tài)安全問題監(jiān)控協(xié)同創(chuàng)新中心項(xiàng)目(JXS-EW-00)和國家自然科學(xué)基金項(xiàng)目(41301023,41371062)聯(lián)合資助.2014-12-05收稿;2015-04-18

    收修改稿.姚靜(1980~),女,助理研究員; E-mail: jyao@niglas.ac.cn。

    ?

    定常風(fēng)對(duì)鄱陽湖水動(dòng)力的影響*

    *國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃“973”項(xiàng)目(2012CB417003)、江西省重大生態(tài)安全問題監(jiān)控協(xié)同創(chuàng)新中心項(xiàng)目(JXS-EW-00)和國家自然科學(xué)基金項(xiàng)目(41301023,41371062)聯(lián)合資助.2014-12-05收稿;2015-04-18

    收修改稿.姚靜(1980~),女,助理研究員; E-mail: jyao@niglas.ac.cn。

    (1:中國科學(xué)院南京地理與湖泊研究所,中國科學(xué)院流域地理學(xué)重點(diǎn)實(shí)驗(yàn)室,南京 210008)

    (2:江西師范大學(xué),鄱陽湖濕地與流域研究教育部重點(diǎn)實(shí)驗(yàn)室,南昌 330022)

    (3:中國科學(xué)院大學(xué),北京 100049)

    摘要:鄱陽湖屬大風(fēng)區(qū),風(fēng)場(chǎng)作為僅次于流域“五河”傾瀉和長江頂托作用的另一重要驅(qū)動(dòng)力,或在某些時(shí)刻影響局部區(qū)域的水流結(jié)構(gòu),進(jìn)而影響局部水體中泥沙、污染物、營養(yǎng)鹽等物質(zhì)的輸移和擴(kuò)散.基于鄱陽湖二維水動(dòng)力數(shù)學(xué)模型,模擬定常風(fēng)場(chǎng)條件下的鄱陽湖流場(chǎng)分布及環(huán)流形式,并與無風(fēng)條件下的水流時(shí)空結(jié)構(gòu)進(jìn)行對(duì)比.結(jié)果表明:3.03m/s的NE向和SSW向定常風(fēng)對(duì)湖泊水位影響微弱;對(duì)流速的影響主要集中在7月中旬至9月底的“湖相”期;其影響區(qū)域主要分布在湖區(qū)中部大湖面偏西岸及東部湖灣,約占湖泊最大水面積的16%;上述區(qū)域出現(xiàn)明顯環(huán)流,環(huán)流結(jié)構(gòu)具有時(shí)空異質(zhì)性特點(diǎn),環(huán)流區(qū)流速普遍增至無風(fēng)時(shí)的兩倍以上;NE向和SSW向風(fēng)場(chǎng)產(chǎn)生的環(huán)流位置相近,方向相反.相比于以往鄱陽湖水動(dòng)力研究中對(duì)風(fēng)場(chǎng)的忽略,本次研究揭示了定常風(fēng)場(chǎng)對(duì)鄱陽湖的重點(diǎn)影響區(qū)域、影響程度及影響形式,可為泥沙及污染物輸移模擬中對(duì)風(fēng)場(chǎng)條件的處理及可能帶來的誤差與誤差的空間分布提供重要依據(jù)。

    關(guān)鍵詞:水動(dòng)力模擬;定常風(fēng);風(fēng)生環(huán)流;鄱陽湖

    湖泊水動(dòng)力過程對(duì)湖泊中的物質(zhì)、能量的輸移擴(kuò)散起著決定性作用.湖泊水流的數(shù)值模擬是泥沙、污染物、生態(tài)等各要素模擬的基礎(chǔ).對(duì)大型寬淺且相對(duì)封閉的湖泊而言,風(fēng)為水流的主要驅(qū)動(dòng)力.風(fēng)場(chǎng)的存在,改變了水流的流速和運(yùn)動(dòng)軌跡,影響泥沙、污染物、營養(yǎng)物質(zhì)等輸移、擴(kuò)散及再懸浮,進(jìn)而影響水質(zhì)及水生生態(tài)系統(tǒng)[1-4].這使得準(zhǔn)確模擬并評(píng)估風(fēng)場(chǎng)對(duì)水動(dòng)力場(chǎng)的影響作用顯得尤為重要。

    在我國,太湖即為典型的風(fēng)生流湖泊,太湖風(fēng)生流的數(shù)值模擬研究,開展得較早且深入[5-10].與太湖湖盆地形相對(duì)平坦、水位變幅較小的水情地貌特點(diǎn)相比,鄱陽湖情況截然不同.鄱陽湖作為我國最大的通江湖泊,受贛江、撫河、信江、饒河、修水“五河”來水及長江來水的雙重影響,屬吞吐型、季節(jié)性湖泊.鄱陽湖地形變化極具空間異質(zhì)性特點(diǎn),湖盆自東向西、由南向北傾斜,南北高程落差達(dá)10m,南北最大長度為173km,最寬處70km,入江通道最窄處約2.8km[11].“五河”來水自上游注入后,在重力作用下,沿河道及洲灘,順地勢(shì)向下游推進(jìn),至湖口匯入長江.受“五河”和長江來水的季節(jié)性變化影響,鄱陽湖水位變幅巨大(可達(dá)10m以上)[11],高水湖相、低水河相,湖流類型以重力型吞吐流為主,風(fēng)生流僅為次要類型[12]。

    以往有關(guān)鄱陽湖水動(dòng)力模擬方面的研究,水流驅(qū)動(dòng)力主要考慮“五河”來水傾瀉和長江頂托作用,通過建立鄱陽湖水動(dòng)力數(shù)學(xué)模型,模擬鄱陽湖的水位、濕地面積等變化評(píng)價(jià)環(huán)境影響,取得了良好的模擬效果[13-17].這些研究均未考慮風(fēng)作為次要因素的影響,但并不影響整體水位模擬效果.但需要強(qiáng)調(diào)的是,鄱陽湖屬大風(fēng)區(qū),據(jù)文獻(xiàn)[18]報(bào)道,星子湖區(qū)6級(jí)以上大風(fēng)年內(nèi)最多達(dá)86d,年均45d,棠蔭站曾實(shí)測(cè)到風(fēng)速31m/s.加之鄱陽湖幅員遼闊,湖泊地形復(fù)雜,河道、島嶼、碟形洼地與灘地并存,大湖面時(shí)期,在局部地形和風(fēng)場(chǎng)的共同作用下,或可能出現(xiàn)風(fēng)生流占主導(dǎo)的情形.目前為止,鄱陽湖盛行風(fēng)在哪些時(shí)段有著頗為顯著的影響作用?在空間尺度上其影響范圍究竟有多大?其影響形式和程度是怎樣的?不同湖區(qū)對(duì)風(fēng)場(chǎng)有著怎樣的響應(yīng)程度?這些問題至今仍缺乏系統(tǒng)認(rèn)識(shí)和深入研究。

    圖1 鄱陽湖地形及水文站點(diǎn)分布Fig.1 Bathymetry of Lake Poyang and locations of hydrological gauging stations in the lake

    鑒于以上幾點(diǎn),本文在鄱陽湖二維水動(dòng)力數(shù)學(xué)模型基礎(chǔ)上,設(shè)置定常風(fēng)場(chǎng),模擬定常風(fēng)場(chǎng)影響下鄱陽湖水流空間結(jié)構(gòu)及變化過程,分析風(fēng)場(chǎng)對(duì)鄱陽湖的主要影響區(qū)域、影響量級(jí)及流態(tài)響應(yīng)形式.本研究結(jié)果不僅有助于闡釋空間湖流形態(tài)對(duì)風(fēng)場(chǎng)的響應(yīng)機(jī)制,還可為泥沙及污染物輸移模擬中對(duì)風(fēng)場(chǎng)條件的處理及可能帶來的誤差與誤差的空間分布提供重要依據(jù),更重要的是該成果可用于進(jìn)一步解釋鄱陽湖水環(huán)境的空間異質(zhì)特性及其影響因素。

    1 研究方法與方案設(shè)計(jì)

    1.1 水動(dòng)力學(xué)模型

    鄱陽湖屬寬淺型湖泊,水體垂向混合較好[19],適用于二維水動(dòng)力數(shù)學(xué)模型.針對(duì)鄱陽湖岸線曲折復(fù)雜、地形起伏多變、水位變幅巨大的特點(diǎn),選擇基于無結(jié)構(gòu)網(wǎng)格的MIKE 21模型.一方面,無結(jié)構(gòu)網(wǎng)格可以更好地?cái)M合地形和岸線,而且通過靈活的加密技術(shù),對(duì)曲折的湖中河道進(jìn)行局部加密,可達(dá)到對(duì)地形的準(zhǔn)確刻畫;另一方面,模型中的干濕判別方法,通過設(shè)置最小干、濕水深,判斷網(wǎng)格是否參與計(jì)算,可準(zhǔn)確模擬鄱陽湖大面積頻繁露灘、淹沒過程。

    基于MIKE 21的鄱陽湖二維水動(dòng)力數(shù)學(xué)模型計(jì)算范圍及岸線邊界根據(jù)湖泊歷史洪水淹沒范圍確定(圖1,基面為85國家高程).流域五河的入湖徑流過程作為水動(dòng)力模型上游開邊界條件,湖泊出口湖口處水位過程作為下游開邊界條件.模型采用三角形網(wǎng)格,為了更好地刻畫湖中河道地形變化,對(duì)主河道區(qū)域進(jìn)行局部網(wǎng)格加密處理,最小網(wǎng)格為70m,而遠(yuǎn)離主河道的洲灘區(qū)域網(wǎng)格分辨率相對(duì)較粗,最大網(wǎng)格尺寸為1500m,共計(jì)剖分網(wǎng)格20450個(gè).該模型作為鄱陽湖水文水動(dòng)力聯(lián)合模擬系統(tǒng)中的水動(dòng)力模塊,已成功再現(xiàn)鄱陽湖水動(dòng)力時(shí)空變化過程[15-16].模型參數(shù)設(shè)置和率定過程詳見文獻(xiàn)[15-16].反映風(fēng)場(chǎng)的表面風(fēng)應(yīng)力計(jì)算公式為:

    (1)

    1.2 鄱陽湖水情

    鄱陽湖屬過水型湖泊,湖流以吞吐流為主[19],平均換水周期約為21d[14].為配合風(fēng)場(chǎng)資料時(shí)間序列,選取2010年作為代表年份進(jìn)行水動(dòng)力數(shù)值模擬.據(jù)2010年實(shí)測(cè)資料統(tǒng)計(jì),“五河”流量年內(nèi)平均為5075.2 m3/s,受暴雨影響,6月21日達(dá)年內(nèi)峰值44719 m3/s.湖泊水位以星子站為代表,最高水位20.28m,最低水位7.75m,水位最大變幅12.53m;平均水位13.84m,比該站1960-2010年多年平均水位高0.49m.基于水位-面積及水位-庫容曲線[22]可得,2010年湖盆區(qū)最大水面面積3265.5km2,對(duì)應(yīng)庫容275.8×108m3,最小水面面積550km2,對(duì)應(yīng)庫容13.1×108m3。

    圖2 2010年鄱陽湖水位驗(yàn)證Fig.2 Validation of water levels at four hydrological gauging stations in 2010

    1.3 模型驗(yàn)證

    圖3 2010年湖口站流量驗(yàn)證Fig.3 Validation of discharge at Hukou hydrological gauging station in 2010

    表1 2010年模型驗(yàn)證誤差

    該模型已對(duì)2000-2008年的水位、水面積、流量、流場(chǎng)等水動(dòng)力指標(biāo)進(jìn)行了全面的率定和驗(yàn)證[15-16],表明模型具有關(guān)鍵過程模擬方面的能力和優(yōu)勢(shì),故此處僅對(duì)2010年的水動(dòng)力過程進(jìn)行進(jìn)一步驗(yàn)證.對(duì)湖區(qū)南北分布的4個(gè)主要水文站點(diǎn)(圖1)的水位過程(圖2)和湖口站的流量過程(圖3)進(jìn)行驗(yàn)證,驗(yàn)證誤差見表1.結(jié)果顯示,水位和流量過程吻合較好,除湖口流量驗(yàn)證誤差稍大外,其余各水位驗(yàn)證相對(duì)誤差均在±3%之內(nèi),確定性系數(shù)均在0.98以上,Nash-Sutcliffe效率系數(shù)也均達(dá)0.95以上.流速驗(yàn)證方面,1)湖口為流量驗(yàn)證,其余4站為水位驗(yàn)證.從以往鄱陽湖水動(dòng)力模擬研究來看,缺乏流速資料驗(yàn)證是一直以來共同面臨的問題.如何在缺乏資料的情況下,有效地驗(yàn)證流速的真實(shí)性和準(zhǔn)確性?自2013年底開始,于星子和湖口附近,進(jìn)行多次ADCP斷面流速監(jiān)測(cè)(圖1),積累了部分流速斷面資料.本文選取2014年5月23日星子斷面和2013年12月4日湖口斷面的流速監(jiān)測(cè)資料,分別代表洪、枯季(高、低水位)的流速情況.由于缺乏2010年之后的水文數(shù)據(jù),無法模擬2013、2014年的同時(shí)刻流速斷面并與實(shí)測(cè)值進(jìn)行比較.為使不同年份的流速具有可比性,選擇本模型水情與2010年相近的時(shí)期進(jìn)行流速對(duì)比.2014年5月23日測(cè)得星子附近斷面水位約為18m,流量為11280.1 m3/s,根據(jù)2010年的模擬情況,選擇5月26日水位17.8m、流量12200.6 m3/s條件下的星子斷面流速與之比較;2013年12月4日測(cè)得湖口斷面水位7.4m、流量1386.9 m3/s,選取計(jì)算的2010年12月8日水位8.1m、流量1207.1 m3/s條件下的湖口斷面流速與之對(duì)比.斷面以河道左岸為起點(diǎn),沿垂直河道方向,至右岸結(jié)束,枯水時(shí)由于水量減少、河寬收縮,實(shí)際監(jiān)測(cè)距離小于模型中的左右岸邊界距離.從斷面流速對(duì)比可知(圖4),在相似水情條件下,流速斷面從量級(jí)到分布趨勢(shì)都比較接近,但畢竟不同年份水情不可能完全相同,流速量值存在一定的偏差,且由于不同年份斷面局部地形變化,存在小范圍內(nèi)流速變化不一致的情形(圖4a).總體而言,這種流速對(duì)比方法在缺乏驗(yàn)證資料的情況下也不失為一種可行的替代方式.綜上所述,模型在2010年水位、流量和流速方面模擬較好,能真實(shí)反映水動(dòng)力變化過程。

    圖4 2010年星子斷面和湖口斷面流速驗(yàn)證Fig.4 Validation of velocity at Xingzi and Hukou cross sections in 2010

    1.4 風(fēng)場(chǎng)分析及模擬方案

    風(fēng)場(chǎng)數(shù)據(jù)采用位于星子縣的鄱陽湖湖泊濕地觀測(cè)研究站標(biāo)準(zhǔn)氣象場(chǎng)實(shí)測(cè)的2010年1-12月每日3次(8:00、14:00、20:00)的風(fēng)速、風(fēng)向資料,其中11、12月只有風(fēng)向資料,風(fēng)速資料缺測(cè).鄱陽湖水面遼闊,各處風(fēng)場(chǎng)并不均勻,但由于缺乏多站點(diǎn)實(shí)測(cè)風(fēng)場(chǎng)資料,因此假定全湖風(fēng)場(chǎng)空間均一,以星子站實(shí)測(cè)風(fēng)場(chǎng)為代表進(jìn)行計(jì)算.鄱陽湖水情復(fù)雜,受“五河”和長江來水雙重影響,不同季節(jié)“河相”、“湖相”相互轉(zhuǎn)換.模擬不同水情疊加變化的風(fēng)場(chǎng)組合,難以區(qū)分風(fēng)場(chǎng)影響.本文側(cè)重關(guān)注鄱陽湖盛行風(fēng)風(fēng)向?qū)λ畡?dòng)力場(chǎng)的影響,因此從最基本的定常風(fēng)模擬入手,統(tǒng)計(jì)出現(xiàn)頻次最多的風(fēng)向及其對(duì)應(yīng)的平均風(fēng)速,模擬該風(fēng)向、風(fēng)速對(duì)不同季節(jié)、不同水情的湖泊水流影響方式及差異,而不同風(fēng)速大小的影響差異需要在此基礎(chǔ)上做進(jìn)一步測(cè)試研究,本文中不予以考慮。

    以往觀測(cè)及研究結(jié)果[12,18,23]表明,除6-8月盛行偏南風(fēng)外,鄱陽湖全年大多時(shí)間為偏北風(fēng).從星子站2010年全年及6-8月的風(fēng)玫瑰圖(圖5)可知,星子站全年的常風(fēng)向?yàn)镹E向,6-8月的常風(fēng)向?yàn)镾SW向,與以往結(jié)論相吻合。

    圖5 星子站風(fēng)玫瑰圖Fig.5 Rose diagram of wind direction at Xingzi station

    分別模擬2010年整年無風(fēng)、NE向定常風(fēng)、SSW向定常風(fēng)3種情景的水動(dòng)力變化過程.根據(jù)實(shí)測(cè)資料統(tǒng)計(jì)的NE向風(fēng)對(duì)應(yīng)的平均風(fēng)速為3.03m/s,將該值設(shè)為NE向定常風(fēng)風(fēng)速,為方便比較,SSW向定常風(fēng)風(fēng)速也取為3.03m/s。

    2 結(jié)果與分析

    2.1 風(fēng)場(chǎng)對(duì)主河道的影響

    為分析風(fēng)場(chǎng)對(duì)鄱陽湖各水動(dòng)力要素的影響,首先提取4個(gè)水文站點(diǎn)有風(fēng)、無風(fēng)條件下的水位、流速、流向過程曲線(圖6).由圖6可知,在3.03m/s的風(fēng)速條件下,無論是NE向還是SSW向風(fēng),各站點(diǎn)水位變化過程與無風(fēng)條件下的基本吻合,即該風(fēng)場(chǎng)條件對(duì)水位的影響很微弱.流速、流向方面,7月中旬至9月底的“湖相”期,棠蔭、康山站受風(fēng)的影響相對(duì)明顯,主要是因?yàn)椤昂唷逼谑荛L江頂托作用,重力流變緩,水面開闊,有利于風(fēng)生流的形成,而棠蔭、康山站正處于開闊大湖面。

    2.2 風(fēng)場(chǎng)對(duì)流速場(chǎng)的影響

    為進(jìn)一步分析風(fēng)場(chǎng)對(duì)空間流速場(chǎng)的影響,針對(duì)影響顯著時(shí)段,分別計(jì)算NE向和SSW向風(fēng)場(chǎng)作用下7月20日至9月30日期間全場(chǎng)平均流速分布,并與無風(fēng)條件下的該時(shí)段全場(chǎng)平均流速分布相比.NE向風(fēng)場(chǎng)作用下,湖區(qū)大部分淺灘區(qū)域流速減小,量級(jí)小于0.5cm/s,相反,河道流速存在不同程度的增大;SSW向風(fēng)場(chǎng)作用下,與NE向相反,大部分淺灘區(qū)域流速增大,增幅也在0.5cm/s以內(nèi),而河道流速存在不同程度的減小(圖7).主要原因?yàn)椋簻\灘相對(duì)河道而言,水深較淺,更易受風(fēng)的影響,NE向風(fēng)與主要流向相逆,對(duì)淺灘水流產(chǎn)生阻礙作用,而SSW向風(fēng)順應(yīng)主流向,對(duì)淺灘水流有一定的推動(dòng)作用,主河道因水深較深,流速較大,在重力流作用下并不直接受風(fēng)的影響,而是在周圍淺灘水體流速變化帶動(dòng)下,產(chǎn)生反方向的流速補(bǔ)償.無論哪種風(fēng)向,在湖區(qū)中部湖灣及中部西岸,均存在明顯的流速增大區(qū),流速普遍增至兩倍以上,最大增幅為2.0~2.5cm/s.對(duì)照流場(chǎng)圖可以發(fā)現(xiàn),上述區(qū)域流速增大主要是由環(huán)流引起的,由于地處開闊的大湖面,流速相對(duì)緩慢,而水深變化較大,極易形成風(fēng)生環(huán)流.明顯環(huán)流區(qū)面積可達(dá)600km2,約占湖泊最大水面積的16%。

    2.3 典型環(huán)流區(qū)分析

    鑒于上節(jié)流速變化的空間分析結(jié)果,結(jié)合流場(chǎng)隨時(shí)間的變化過程,劃分出風(fēng)生環(huán)流比較明顯的3個(gè)局部區(qū)域(圖8),以便進(jìn)行局部流場(chǎng)的詳細(xì)比對(duì).圖中A區(qū)位于東部湖灣,B區(qū)位于東南部湖灣,C區(qū)位于湖區(qū)中部偏西.重點(diǎn)分析受風(fēng)場(chǎng)影響最為明顯的時(shí)段,即7月中旬至9月底。

    圖6 有風(fēng)、無風(fēng)條件湖區(qū)水位、流速與流向變化過程曲線Fig.6 Comparisons of water level, current speed and current direction of four hydrological gauging stations with no wind, NE wind and SSW wind

    圖7 有風(fēng)與無風(fēng)條件相比全場(chǎng)平均流速變幅分布(7月20日-9月30日,正值代表流速增大,負(fù)值代表流速減小)Fig.7 Spatial distribution of average current speed difference between wind conditions and no wind condition(July 20th-September 30th, positive for increasing current speed and negative for decreasing current speed)

    圖8 典型環(huán)流區(qū)劃分Fig.8 Typical circulation regions

    從流場(chǎng)變化來看,7月中旬至9月底,A區(qū)、B區(qū)在兩種風(fēng)場(chǎng)作用下均存在穩(wěn)定的環(huán)流,故選擇8月8日的流場(chǎng)作為典型流場(chǎng)進(jìn)行分析。

    A區(qū)(圖9)地形南高北低,北部近岸區(qū)常年有水,但流速極小,近乎停滯.在NE向風(fēng)場(chǎng)作用下,東南部的湖濱淺水區(qū)水體最易受影響,順風(fēng)向自東北流向西南,脅迫東南入湖口處水流也向西南向偏移,拉動(dòng)北部深水湖區(qū)水體形成補(bǔ)償流,依地形走勢(shì)呈順時(shí)針方向旋轉(zhuǎn)形成最大直徑18km的環(huán)流渦旋;西北部小湖灣,在主河道和束窄的湖岸共同影響下,也形成了小尺度順時(shí)針環(huán)流.而在SSW風(fēng)場(chǎng)作用下,東南部淺水區(qū)仍受風(fēng)場(chǎng)直接影響,流向與風(fēng)向一致,指向東北,導(dǎo)致東部整體流態(tài)呈S型,在與NE向風(fēng)生環(huán)流相同位置,形成了逆時(shí)針環(huán)流,同時(shí)在島嶼的東南部,受微地形和周圍流態(tài)影響,也存在逆時(shí)針環(huán)流渦旋;西北部湖灣雖未形成封閉環(huán)流,但流向發(fā)生了偏逆時(shí)針方向的變形.A區(qū)不同風(fēng)場(chǎng)引起的主要環(huán)流方向符合以往研究規(guī)律,即淺水區(qū)風(fēng)生流向與風(fēng)向一致,深水區(qū)則與風(fēng)向相反[24-25]。

    圖9 A區(qū)流場(chǎng)及環(huán)流示意圖Fig.9 Flow field and circulation patterns in region A

    B區(qū)(圖10)在NE向風(fēng)場(chǎng)作用下,產(chǎn)生順時(shí)針環(huán)流;在SSW風(fēng)場(chǎng)作用下,產(chǎn)生逆時(shí)針環(huán)流,均與A區(qū)的主要環(huán)流方向一致。

    圖10 B區(qū)流場(chǎng)及環(huán)流示意圖Fig.10 Flow field and circulation patterns in region B

    簡而言之,東部湖灣A、B區(qū),在NE向風(fēng)場(chǎng)作用下,會(huì)形成穩(wěn)定的順時(shí)針環(huán)流,而在SSW風(fēng)場(chǎng)作用下,相同區(qū)域形成逆時(shí)針環(huán)流.風(fēng)生環(huán)流的存在使得東部湖灣主河道外的大部分區(qū)域流速普遍增大,但流速量級(jí)總體仍較小。

    C區(qū)流場(chǎng)變化最為復(fù)雜,環(huán)流并不穩(wěn)定,其產(chǎn)生位置和尺度隨時(shí)空變化顯著.對(duì)照水位過程曲線(圖2)可知,7月底至8月底,水位急劇下降;8月底至9月底,水位變化處于相對(duì)滯緩期.因此,分別選取8月8日和9月15日為兩時(shí)期典型代表進(jìn)行流場(chǎng)對(duì)比分析(圖11).無風(fēng)狀態(tài),9月15日流速與8月8日相比整體偏小,這與兩時(shí)段的水位變化相呼應(yīng).在NE向風(fēng)場(chǎng)作用下,9月15日西北部出現(xiàn)多處尺度不一的環(huán)流,以逆時(shí)針為主,因該區(qū)域處于大湖面,水深較淺,水下地形復(fù)雜,加之受島嶼掩蔽作用,極易生成環(huán)流;該區(qū)南部,大湖面沿岸帶,水流明顯發(fā)生順時(shí)針偏轉(zhuǎn),從湖盆地形(圖1)可知,該處地形變化較大,對(duì)風(fēng)場(chǎng)作用較為敏感;而8月8日水位急劇下降,流速較大,相對(duì)而言受風(fēng)影響較小,因而整個(gè)西岸都沒有形成環(huán)流,只在局部有微弱的流向變化.SSW向風(fēng)場(chǎng)作用下,9月15日,島嶼南側(cè)水流偏轉(zhuǎn)呈“幾”字形,并伴有順時(shí)針環(huán)流;西岸在NE向出現(xiàn)明顯偏轉(zhuǎn)流的相近位置,形成逆時(shí)針環(huán)流;而8月8日,雖然也有環(huán)流形成,但整體流態(tài)的變化程度不及前者。

    圖11 C區(qū)流場(chǎng)及環(huán)流示意圖Fig.11 Flow field and circulation patterns in region C

    總體而言,在NE向和SSW向風(fēng)場(chǎng)作用下,C區(qū)在西岸及湖中島的南側(cè)都易生成環(huán)流,環(huán)流位置大致相同,但方向相反;水位停滯期整體流速偏小,受風(fēng)場(chǎng)影響更明顯,環(huán)流尺度更大。

    3 討論

    本文利用有限的風(fēng)場(chǎng)資料,采用定常風(fēng)場(chǎng),模擬鄱陽湖風(fēng)生流的產(chǎn)生時(shí)間、位置、環(huán)流形態(tài)等.真實(shí)的風(fēng)場(chǎng)均為非定常風(fēng),本文對(duì)于定常風(fēng)的影響模擬研究,并不是為了反演某真實(shí)風(fēng)場(chǎng)過程,而是針對(duì)真實(shí)風(fēng)場(chǎng)中出現(xiàn)頻率較高的風(fēng)向,研究其影響形式和區(qū)域.定常風(fēng)風(fēng)場(chǎng)作用下風(fēng)向和風(fēng)速固定不變,故環(huán)流形式相對(duì)穩(wěn)定,持續(xù)時(shí)間也較長;而真實(shí)風(fēng)場(chǎng)作用下,由于風(fēng)向、風(fēng)速的頻繁變化,導(dǎo)致環(huán)流的大小、方向及結(jié)構(gòu)形式隨之變化,并非定常風(fēng)下的一成不變,環(huán)流持續(xù)時(shí)間也與某風(fēng)向持續(xù)時(shí)間相關(guān).此外需要說明的是,本文定常風(fēng)風(fēng)速采用的是資料統(tǒng)計(jì)的平均值,對(duì)極端風(fēng)速可能引起的短暫情景未予以考慮。

    定常風(fēng)風(fēng)速根據(jù)2010年1-10月的風(fēng)速序列做平均,缺少11-12月大風(fēng)月資料,故得到的平均風(fēng)速3.03m/s比文獻(xiàn)[12]里的略小.但本研究主要針對(duì)鄱陽湖的兩種盛行風(fēng),側(cè)重風(fēng)向的影響,風(fēng)速大小的影響并非研究重點(diǎn).類似的風(fēng)向?qū)椿蚝乘畡?dòng)力影響研究已有不少,如Razmi等[26]通過設(shè)置多種風(fēng)向、相同風(fēng)速(3m/s)的定常風(fēng),測(cè)試大型湖泊中開敞式湖灣湖流結(jié)構(gòu)對(duì)不同風(fēng)場(chǎng)的響應(yīng),結(jié)果表明:在大型湖泊中,即使是很小的湖灣,也有可能產(chǎn)生復(fù)雜流態(tài);風(fēng)向發(fā)生小角度改變,湖流結(jié)構(gòu)即發(fā)生明顯變化.Anderson等[27]通過8種主風(fēng)向和2種風(fēng)速組合的數(shù)值模擬發(fā)現(xiàn),圣克萊爾湖及其三角洲對(duì)某些風(fēng)向較為敏感,響應(yīng)程度也存在較大的空間異質(zhì)性特點(diǎn),某些風(fēng)場(chǎng)組合條件下,甚至?xí)霈F(xiàn)倒灌流.Guo等[2]模擬了同一風(fēng)應(yīng)力的4種定常風(fēng)(NE、NW、SW、SE)對(duì)切薩皮克海灣流場(chǎng)和鹽度場(chǎng)的時(shí)空結(jié)構(gòu)影響,結(jié)果顯示,上、下游海灣對(duì)不同風(fēng)向的響應(yīng)并不相同.本研究與上述研究結(jié)論相似,即不同的風(fēng)向?qū)α鲌?chǎng)的影響形式不同,不同區(qū)域?qū)ν伙L(fēng)向的響應(yīng)程度也存在較大差異;區(qū)別在于本文沒有對(duì)風(fēng)向影響做廣泛測(cè)試,只是集中模擬分析了鄱陽湖兩種盛行風(fēng)的影響。

    以往風(fēng)生流的研究結(jié)果[28-30]表明:相似風(fēng)向條件下,風(fēng)速增大或減小主要影響流速大小,對(duì)環(huán)流空間分布結(jié)構(gòu)不會(huì)產(chǎn)生顯著影響.由此認(rèn)為,本文的風(fēng)速模擬結(jié)果基本能代表NE、SSW風(fēng)向下的流態(tài)結(jié)構(gòu),可據(jù)此大致了解風(fēng)場(chǎng)對(duì)鄱陽湖水流的作用方式.同時(shí),風(fēng)生流強(qiáng)弱也受吞吐流量影響[24],隨著風(fēng)速量級(jí)增大或者吞吐流量的變化,是否會(huì)在其他區(qū)域產(chǎn)生新的環(huán)流,環(huán)流形式或形成時(shí)間會(huì)不會(huì)有所差異,需要在以后的研究中做進(jìn)一步的數(shù)值模擬測(cè)試。

    有風(fēng)、無風(fēng)條件下的湖泊水流空間結(jié)構(gòu)的差異體現(xiàn)了風(fēng)場(chǎng)對(duì)鄱陽湖的作用,而無風(fēng)條件下的水動(dòng)力場(chǎng)則反映了“五河”和長江的作用,因此本文結(jié)果也可闡釋“五河”、長江共同作用及風(fēng)場(chǎng)作用對(duì)鄱陽湖水動(dòng)力場(chǎng)的貢獻(xiàn)差異.即:對(duì)大部分湖區(qū)尤其是主河道而言,水流主要驅(qū)動(dòng)力為“五河”和長江,只在湖相期,湖區(qū)中部開闊區(qū)部分淺灘和河道,風(fēng)場(chǎng)的作用可與“五河”和長江作用相當(dāng);對(duì)流速極小的局部湖灣如東部湖灣,風(fēng)場(chǎng)貢獻(xiàn)甚至可超過“五河”和長江.由于本文只選用一年情景,并未考慮“五河”和長江的多種水情及與不同風(fēng)場(chǎng)的組合,因此對(duì)兩種驅(qū)動(dòng)力的貢獻(xiàn)無法給出詳細(xì)的對(duì)比,這將在下一步工作中繼續(xù)完善。

    風(fēng)生環(huán)流的潛在重要作用使得營養(yǎng)鹽在局部區(qū)域富集,影響浮游植物的分布,有可能導(dǎo)致藻類水華.Ishikawa等[31]通過研究琵琶湖中的藍(lán)藻暴發(fā),指出環(huán)流渦旋對(duì)營養(yǎng)物質(zhì)和藻類運(yùn)輸產(chǎn)生重要作用,進(jìn)而對(duì)藍(lán)藻分布產(chǎn)生重要影響.鄱陽湖水質(zhì)總體較好,但是近年的監(jiān)測(cè)資料對(duì)比發(fā)現(xiàn),鄱陽湖水質(zhì)呈下降趨勢(shì),某些時(shí)段局部湖區(qū)已經(jīng)發(fā)生富營養(yǎng)化,并且整個(gè)湖區(qū)正緩慢向富營養(yǎng)化趨勢(shì)發(fā)展[32-34].Wu等[35]給出鄱陽湖夏季葉綠素a濃度空間分布(2009-2012年),其中高濃度區(qū)主要分布在東部湖灣,與本研究中的A、B區(qū)域高度吻合,同時(shí)指出,該區(qū)域夏季曾觀測(cè)到藍(lán)藻水華.這恰好印證了本文環(huán)流區(qū)分布的合理性,同時(shí)本研究成果也可從水動(dòng)力角度對(duì)這一現(xiàn)象的解釋給予理論支撐。

    鄱陽湖濕地生態(tài)系統(tǒng)的可持續(xù)發(fā)展與洲灘濕地的發(fā)育和演化密不可分[36].在淺水湖泊中,風(fēng)成浪對(duì)泥沙再懸浮起著重要作用[37],而風(fēng)生流則影響泥沙的遷移、輸運(yùn)、堆積,進(jìn)而影響泥沙沖淤、洲灘發(fā)育.從本研究結(jié)果來看,鄱陽湖局部湖區(qū)對(duì)風(fēng)場(chǎng)作用響應(yīng)明顯,說明風(fēng)場(chǎng)對(duì)局部湖區(qū)洲灘發(fā)育存在潛在影響。

    4 結(jié)論

    本文基于鄱陽湖二維水動(dòng)力數(shù)學(xué)模型,通過無風(fēng)、NE向和SSW向定常風(fēng)(3.03m/s)3種情景模擬,對(duì)考慮及不考慮風(fēng)場(chǎng)對(duì)鄱陽湖水動(dòng)力場(chǎng)的影響,以及兩種定常風(fēng)場(chǎng)的影響差異開展了對(duì)比分析.主要結(jié)論為:

    1) 該量級(jí)風(fēng)場(chǎng)對(duì)湖泊水位影響微弱;對(duì)流速影響主要集中在7月中旬至9月底的“湖相”期,尤其水位變化越緩,流速越小,影響越明顯。

    2) 風(fēng)場(chǎng)影響區(qū)域主要分布在湖區(qū)中部大湖面偏西岸及東部湖灣;上述區(qū)域出現(xiàn)明顯環(huán)流,環(huán)流區(qū)面積可達(dá)600km2,約占湖泊最大水面積的16%;同一風(fēng)場(chǎng)作用下不同區(qū)域環(huán)流結(jié)構(gòu)及持續(xù)時(shí)間也不盡相同,具有顯著的時(shí)空異質(zhì)性特點(diǎn);無論哪種風(fēng)向作用,環(huán)流區(qū)流速平均增至無風(fēng)狀態(tài)時(shí)兩倍以上,最大增幅2~2.5cm/s;其它大部分湖區(qū)流速變幅均在±0.5cm/s以內(nèi)。

    3) NE向和SSW向風(fēng)場(chǎng)產(chǎn)生的環(huán)流區(qū)位置相近,環(huán)流尺度相似,但方向相反。

    4) 涉及“湖相”時(shí)期湖區(qū)中部西岸及東部湖灣等局部湖區(qū)的水動(dòng)力、物質(zhì)輸移、泥沙沖淤、洲灘發(fā)育等研究時(shí),風(fēng)的影響不容忽視。

    致謝:感謝鄱陽湖湖泊濕地觀測(cè)研究站提供風(fēng)場(chǎng)資料。

    參考文獻(xiàn)5

    [1]Bloesch J. Mechanisms, measurement and importance of sediment resuspension in lakes.MarineandFreshwaterResearch, 1995, 46(1): 295-304。

    [2]Guo XY, Valle-Levinson A. Wind effects on the lateral structure of density-driven circulation in Chesapeake Bay.ContinentalShelfResearch, 2008, 28: 2450-2471。

    [3]秦伯強(qiáng), 胡維平, 陳偉民等. 太湖梅梁灣水動(dòng)力及相關(guān)過程的研究. 湖泊科學(xué), 2000, 12(4): 327-334. DOI 10. 18307/2000. 0406。

    [4]池俏俏, 朱廣偉, 張戰(zhàn)平等. 風(fēng)浪擾動(dòng)對(duì)太湖水體懸浮物重金屬含量的影響. 湖泊科學(xué), 2006, 18(5): 495-498. DOI 10. 18307/2006. 0509。

    [5]梁瑞駒, 仲金華. 太湖風(fēng)生流的三維數(shù)值模擬. 湖泊科學(xué), 1994, 6(4): 289-297. DOI 10. 18307/1994. 0401。

    [6]逄勇, 濮培民, 高光等. 非均勻風(fēng)場(chǎng)作用下太湖風(fēng)成流風(fēng)涌水的數(shù)值模擬及驗(yàn)證. 海洋湖沼通報(bào), 1994, (4): 9-15。

    [7]逄勇, 濮培民. 太湖風(fēng)生流三維數(shù)值模擬試驗(yàn). 地理學(xué)報(bào), 1996, 51(4): 322-328。

    [8]胡維平, 濮培民, 秦伯強(qiáng). 太湖水動(dòng)力學(xué)三維數(shù)值試驗(yàn)研究——1. 風(fēng)生流和風(fēng)涌增減水的三維數(shù)值模擬. 湖泊科學(xué), 1998, 10(4): 17-25. DOI 10. 18307/1998. 0403。

    [9]胡維平, 濮培民, 秦伯強(qiáng). 太湖水動(dòng)力學(xué)三維數(shù)值試驗(yàn)研究——2. 典型風(fēng)場(chǎng)風(fēng)生流的數(shù)值計(jì)算. 湖泊科學(xué), 1998, 10(4): 26-34. DOI 10. 18307/1998. 0404。

    [10]朱永春, 蔡啟銘. 太湖梅梁灣三維水動(dòng)力學(xué)的研究——I. 模型的建立及結(jié)果分析. 海洋與湖沼, 1998, 29(1): 79-85。

    [11]譚國良, 郭生練, 王俊等. 鄱陽湖生態(tài)經(jīng)濟(jì)區(qū)水文水資源演變規(guī)律研究. 北京: 中國水利水電出版社, 2013。

    [12]熊道光. 鄱陽湖湖流特性分析與研究. 海洋與湖沼, 1991, 22(3): 200-207。

    [13]汪迎春, 賴錫軍, 姜加虎等. 三峽水庫調(diào)節(jié)典型時(shí)段對(duì)鄱陽湖濕地水情特征的影響. 湖泊科學(xué), 2011, 23(2): 191-195. DOI 10. 18307/2011. 0205。

    [14]賴錫軍, 姜家虎, 黃群等. 鄱陽湖二維水動(dòng)力和水質(zhì)耦合數(shù)值模擬. 湖泊科學(xué), 2011, 23(6): 893-902. DOI 10. 18307/2011. 0611。

    [15]李云良, 張奇, 姚靜等. 鄱陽湖湖泊流域系統(tǒng)水文水動(dòng)力聯(lián)合模擬. 湖泊科學(xué), 2013, 25(2): 227-235. DOI 10. 18307/2013. 0208。

    [16]Li YL, Zhang Q, Yao Jetal. Hydrodynamic and hydrological modeling of the Poyang Lake catchment system in China.JournalofHydrologicEngineering, 2014, 19(3): 607-616。

    [17]王鵬, 賴格英, 黃小蘭. 鄱陽湖水利樞紐工程對(duì)湖泊水位變化影響的模擬. 湖泊科學(xué), 2014, 26(1): 29-36. DOI 10. 18307/2014. 0104。

    [18]尹宗賢, 張俊才. 鄱陽湖水文特征(Ⅱ). 海洋與湖沼, 1987, 18(2): 208-214。

    [19]程時(shí)長, 盧兵. 鄱陽湖湖流特征. 江西水利科技, 2003, 29(2): 105-108。

    [20]Wu J. The sea surface is aerodynamically rough even under light winds.BoundaryLayerMeteorology, 1994, 69: 149-158。

    [21]Wu J. Wind-stress coefficients over sea surface and near neutral conditions-a revisit.JournalofPhysicalOceanography, 1980, 10: 727-740。

    [22]譚國良, 李國文, 喻中文等. 鄱陽湖動(dòng)態(tài)水位-面積-容積關(guān)系研究. 中國水利學(xué)會(huì)2013學(xué)術(shù)年會(huì)論文集—S2湖泊治理開發(fā)與保護(hù). 北京: 中國水利水電出版社, 2013: 736-742。

    [23]程時(shí)長, 李良文. 鄱陽湖的風(fēng)情及風(fēng)浪特性. 海洋湖沼通報(bào), 1993, 4: 1-12。

    [24]韓紅娟, 胡維平, 晉義泉. 風(fēng)速變化對(duì)竹湖流場(chǎng)結(jié)構(gòu)影響的數(shù)值試驗(yàn). 海洋與湖沼, 2008, 39(6): 567-576。

    [25]Schoen JH, Stretch DD, Tirok K. Wind-driven circulation patterns in a shallow estuarine lake: St Lucia, Estuarine, South Africa.CoastalandShelfScience, 2014, 146: 49-59。

    [26]Razmi AM, Barry DA, Bakhtyar Retal. Current variability in a wide and open lacustrine embayment in Lake Geneva(Switzerland).JournalofGreatLakesResearch, 2013, 39: 455-465。

    [27]Anderson EJ, Schwab DJ. Relationships between wind-driven and hydraulic flow in Lake St. Clair and the St. Clair River Delta.JournalofGreatLakesResearch, 2011, 37: 147-158。

    [28]姜加虎. 貢湖及其相關(guān)水域風(fēng)生流模擬研究. 海洋湖沼通報(bào), 1997, (4): 1-7。

    [29]曹慧江, 楊蕓, 朱建榮. 滴水湖夏季環(huán)流的三維數(shù)值模擬. 華東師范大學(xué)學(xué)報(bào): 自然科學(xué)版, 2006, 4: 78-83, 90。

    [30]韓龍喜, 陸東燕, 李洪晶等. 高鹽度湖泊艾比湖風(fēng)生流三維數(shù)值模擬. 水科學(xué)進(jìn)展, 2011, 22(1): 97-103。

    [31]Ishikawa K, Kumagai M, Vincent WFetal. Transport and accumulation of bloom-forming cyanobacteria in a large, mid-latitude lake: the gyreMicrocystishypothesis.Limnology, 2002, 3(2): 87-96。

    [32]呂蘭軍. 鄱陽湖富營養(yǎng)化調(diào)查與評(píng)價(jià). 湖泊科學(xué), 1996, 8(3): 241-247. DOI 10. 18307/1996. 0308。

    [33]王毛蘭, 周文斌, 胡春華. 鄱陽湖區(qū)水體氮、磷污染狀況分析. 湖泊科學(xué), 2008, 20(3): 334-338. DOI 10. 18307/2008. 0311。

    [34]胡春華, 周文斌, 王毛蘭等. 鄱陽湖氮磷營養(yǎng)鹽變化特征及潛在性富營養(yǎng)化評(píng)價(jià). 湖泊科學(xué), 2010, 22(5): 723-728. DOI 10. 18307/2010. 0514。

    [35]Wu ZS, He H, Cai YJetal. Spatial distribution of chlorophyll a and its relationship with the environment during summer in Lake Poyang: a Yangtze-connected lake.Hydrobiologia, 2014, 732: 61-70。

    [36]葛剛, 紀(jì)偉濤, 劉成林等. 鄱陽湖水利樞紐工程與濕地生態(tài)保護(hù). 長江流域資源與環(huán)境, 2010, 19(6): 606-613。

    [37]Qian J, Zheng SS, Wang PFetal. Experimental study on sediment resuspension in Taihu Lake under different hydrodynamic disturbances.JournalofHydrodynamics, 2011, 23(6): 826-833。

    J.LakeSci.(湖泊科學(xué)), 2016, 28(1): 225-236

    ?2016 byJournalofLakeSciences

    The influence of uniform winds on hydrodynamics of Lake Poyang

    YAO Jing1, ZHANG Qi1,2**, LI Yunliang1& LI Mengfan1,3

    (1:KeyLaboratoryofWatershedGeographicSciences,NanjingInstituteofGeographyandLimnology,ChineseAcademyofSciences,Nanjing210008,P.R.China)

    (2:KeyLaboratoryofPoyangLakeWetlandandWatershedResearch,MinistryofEducation,JiangxiNormalUniversity,Nanchang330022,P.R.China)

    (3:UniversityofChineseAcademyofSciences,Beijing100049,P.R.China)

    Abstract:Lake Poyang is a strongly wind-affected region. The wind stress, as the secondary driver after the catchment forcing and effect of the Yangtze River, may influence the flow pattern of local lake body at certain times, which affects the transport and diffusion processes of sediment, pollutants and nutrient substances. Based on two-dimensional hydrodynamic model, the flow distributions and circulation patterns were simulated under uniform wind conditions. Comparisons were made between cases of no wind and uniform wind in order to investigate the spatiotemporal patterns of the flow. The results showed that, for NE and SSW steady winds at 3.03m/s speed, the impact of wind on water level was negligible, and the effect on flow was mainly restricted to the west coast and eastern bay between the middle of July and the end of September. In the affected regions(about 16% of the total lake area), obvious circulations were noted, which showed spatial and temporal variability, and where average flow velocity was approximately twice as much as that for no wind condition. NE and SSW wind-driven circulations occurred in similar locations but corresponding flow directions were opposite. In comparison to previous hydrodynamic simulations of Lake Poyang that ignored wind effect, this study enhanced the understanding of the location, magnitude and pattern of uniform wind forcing effect on flow. It hints the potential uncertainty in sediment and pollutant transport modelling when dealing with wind conditions。

    Keywords:Hydrodynamic simulation; uniform wind; wind-driven circulation; Lake Poyang

    通信作者姚靜1,張奇1,2*;E-mail: qzhang@niglas.ac.cn.,李云良1,李夢(mèng)凡1,3

    DOI10.18307/2016.0126

    猜你喜歡
    鄱陽湖
    鄱陽湖水系之潦河
    鄱陽湖
    快樂語文(2021年8期)2021-05-06 06:09:08
    《鄱陽湖生態(tài)系列插畫》
    鄱陽湖好風(fēng)光
    老友(2017年4期)2017-02-09 00:26:04
    鄱陽湖鳥語
    中國攝影家(2014年6期)2014-04-29 14:54:47
    鄱陽湖國家級(jí)自然保護(hù)區(qū)濕地植被的干旱響應(yīng)及影響因素
    淺析太平軍鄱陽湖大捷
    軍事歷史(1985年4期)1985-08-20 07:26:34
    久久 成人 亚洲| 最黄视频免费看| 国产精品一区二区在线观看99| 男女无遮挡免费网站观看| 人人妻人人看人人澡| 国产黄片美女视频| 日韩,欧美,国产一区二区三区| 伦理电影大哥的女人| 国产精品一区二区三区四区免费观看| 亚洲国产高清在线一区二区三| 免费av不卡在线播放| 亚洲丝袜综合中文字幕| 成人高潮视频无遮挡免费网站| 亚洲久久久国产精品| 欧美日韩综合久久久久久| av网站免费在线观看视频| 1000部很黄的大片| 最黄视频免费看| 午夜激情久久久久久久| 国产男女内射视频| 91久久精品国产一区二区三区| 超碰av人人做人人爽久久| 久久韩国三级中文字幕| 热99国产精品久久久久久7| 日韩 亚洲 欧美在线| 久久精品久久久久久久性| 伊人久久国产一区二区| 男人添女人高潮全过程视频| 伦精品一区二区三区| 人妻夜夜爽99麻豆av| 一本色道久久久久久精品综合| 精品久久久久久久末码| 亚洲av成人精品一区久久| 欧美日韩综合久久久久久| 久久99热这里只有精品18| 少妇裸体淫交视频免费看高清| 久久久久网色| 欧美97在线视频| 日本免费在线观看一区| 久久久国产一区二区| av不卡在线播放| 观看免费一级毛片| 两个人的视频大全免费| 成人影院久久| 如何舔出高潮| 亚洲成人av在线免费| 菩萨蛮人人尽说江南好唐韦庄| 美女内射精品一级片tv| av在线老鸭窝| 日本免费在线观看一区| 亚洲欧洲日产国产| 熟女人妻精品中文字幕| 国产在线男女| 日本av免费视频播放| 少妇人妻久久综合中文| 久久 成人 亚洲| 秋霞在线观看毛片| 一个人免费看片子| 三级国产精品欧美在线观看| 黑人猛操日本美女一级片| 狂野欧美激情性bbbbbb| 99re6热这里在线精品视频| 欧美日韩在线观看h| 国产精品女同一区二区软件| 熟妇人妻不卡中文字幕| 全区人妻精品视频| 色哟哟·www| 国产美女午夜福利| 欧美人与善性xxx| 国产真实伦视频高清在线观看| 久久人人爽av亚洲精品天堂 | 亚洲,欧美,日韩| 国产精品偷伦视频观看了| 一级二级三级毛片免费看| 超碰97精品在线观看| 高清视频免费观看一区二区| 欧美日韩在线观看h| 午夜免费鲁丝| 狠狠精品人妻久久久久久综合| 五月开心婷婷网| 男女下面进入的视频免费午夜| 欧美成人a在线观看| 高清在线视频一区二区三区| 22中文网久久字幕| 久久ye,这里只有精品| 久久久久性生活片| 亚洲色图av天堂| 亚洲伊人久久精品综合| 99久久精品国产国产毛片| 高清毛片免费看| av天堂中文字幕网| 一级a做视频免费观看| 少妇 在线观看| 人妻一区二区av| freevideosex欧美| 亚洲婷婷狠狠爱综合网| tube8黄色片| 国产免费一级a男人的天堂| 九色成人免费人妻av| 国产美女午夜福利| 欧美+日韩+精品| 欧美日韩视频精品一区| 噜噜噜噜噜久久久久久91| 亚洲精品国产成人久久av| 久久人妻熟女aⅴ| 寂寞人妻少妇视频99o| 亚洲精品乱码久久久久久按摩| 26uuu在线亚洲综合色| 有码 亚洲区| 久久午夜福利片| 国产视频内射| 国产成人精品一,二区| 黄色配什么色好看| 日韩一本色道免费dvd| 99热网站在线观看| 亚洲色图综合在线观看| 水蜜桃什么品种好| 高清在线视频一区二区三区| 国产欧美亚洲国产| 国产男人的电影天堂91| 性色av一级| 九草在线视频观看| 国产成人免费观看mmmm| 国产视频内射| 一级毛片黄色毛片免费观看视频| 自拍欧美九色日韩亚洲蝌蚪91 | 成年免费大片在线观看| 美女脱内裤让男人舔精品视频| 国产男人的电影天堂91| 香蕉精品网在线| 国产白丝娇喘喷水9色精品| 美女国产视频在线观看| a 毛片基地| 国产成人一区二区在线| 国产一区亚洲一区在线观看| 热re99久久精品国产66热6| 五月开心婷婷网| 爱豆传媒免费全集在线观看| 国产美女午夜福利| 日韩,欧美,国产一区二区三区| 日日撸夜夜添| 久久久久国产精品人妻一区二区| 国产精品一区二区在线观看99| 777米奇影视久久| 亚洲精品第二区| 国产成人freesex在线| 日韩制服骚丝袜av| 国产淫片久久久久久久久| 熟女av电影| 精品国产三级普通话版| 日韩 亚洲 欧美在线| 国产一区二区三区综合在线观看 | 成年女人在线观看亚洲视频| 久久精品国产自在天天线| 能在线免费看毛片的网站| 亚洲丝袜综合中文字幕| 搡女人真爽免费视频火全软件| 日本欧美国产在线视频| 亚洲国产日韩一区二区| 亚洲国产欧美人成| 亚洲精品一区蜜桃| 日日啪夜夜撸| 日韩免费高清中文字幕av| 国产欧美日韩精品一区二区| 精品久久久久久久久亚洲| 国语对白做爰xxxⅹ性视频网站| 在线观看免费视频网站a站| 自拍偷自拍亚洲精品老妇| 美女福利国产在线 | 久久久精品94久久精品| 直男gayav资源| 成人毛片a级毛片在线播放| 深爱激情五月婷婷| 精品一区二区三卡| 亚洲美女视频黄频| 麻豆精品久久久久久蜜桃| 性高湖久久久久久久久免费观看| 免费久久久久久久精品成人欧美视频 | 久久久久久久亚洲中文字幕| 成人黄色视频免费在线看| 精品久久久精品久久久| 高清视频免费观看一区二区| 国产男女内射视频| 欧美 日韩 精品 国产| 女人久久www免费人成看片| 亚洲av福利一区| 少妇熟女欧美另类| av福利片在线观看| 精品少妇久久久久久888优播| 久久久a久久爽久久v久久| 亚洲精品自拍成人| 亚洲精品国产av成人精品| 久久亚洲国产成人精品v| 亚洲伊人久久精品综合| 成人二区视频| av网站免费在线观看视频| 亚洲欧美日韩卡通动漫| 国产午夜精品一二区理论片| 97超视频在线观看视频| 18禁动态无遮挡网站| 大话2 男鬼变身卡| 91精品伊人久久大香线蕉| 精品视频人人做人人爽| 久久99热这里只有精品18| 亚洲人与动物交配视频| 国产av码专区亚洲av| 搡老乐熟女国产| 久久久久久久久久久丰满| 中文字幕av成人在线电影| 久久 成人 亚洲| 日韩一区二区视频免费看| 26uuu在线亚洲综合色| 在线免费十八禁| 看免费成人av毛片| 成人特级av手机在线观看| av线在线观看网站| 91aial.com中文字幕在线观看| 插阴视频在线观看视频| 久久鲁丝午夜福利片| 久久人人爽人人爽人人片va| 久久久久久人妻| 亚洲一区二区三区欧美精品| 能在线免费看毛片的网站| 免费看不卡的av| 色吧在线观看| 中文字幕免费在线视频6| 日韩人妻高清精品专区| 亚洲av免费高清在线观看| 99国产精品免费福利视频| 国产欧美另类精品又又久久亚洲欧美| 中文字幕制服av| 国产乱人视频| 午夜老司机福利剧场| 中国美白少妇内射xxxbb| 十分钟在线观看高清视频www | 极品少妇高潮喷水抽搐| 日日啪夜夜撸| 国产精品99久久99久久久不卡 | 国产欧美亚洲国产| 99九九线精品视频在线观看视频| 国产欧美日韩精品一区二区| 爱豆传媒免费全集在线观看| 精品亚洲乱码少妇综合久久| 成人亚洲欧美一区二区av| 亚洲国产精品专区欧美| 国产成人精品婷婷| 国产精品精品国产色婷婷| 亚洲精品日韩av片在线观看| 免费久久久久久久精品成人欧美视频 | 国产爱豆传媒在线观看| 丰满乱子伦码专区| 久久精品国产亚洲网站| 免费看不卡的av| 日本午夜av视频| 91精品国产国语对白视频| 久久久久久久久久久丰满| 在线天堂最新版资源| 久久久国产一区二区| 少妇人妻精品综合一区二区| 91狼人影院| 精品久久久久久电影网| 18禁动态无遮挡网站| 日韩一本色道免费dvd| 18禁裸乳无遮挡免费网站照片| 亚洲在久久综合| 一级a做视频免费观看| 精品人妻一区二区三区麻豆| 男女下面进入的视频免费午夜| 少妇熟女欧美另类| 91午夜精品亚洲一区二区三区| 一区二区三区乱码不卡18| 久久鲁丝午夜福利片| 伦理电影免费视频| 免费人成在线观看视频色| 国产伦精品一区二区三区视频9| 久久99热6这里只有精品| 在线看a的网站| 欧美xxxx黑人xx丫x性爽| 中国美白少妇内射xxxbb| 国产国拍精品亚洲av在线观看| 熟妇人妻不卡中文字幕| 欧美bdsm另类| 国产精品久久久久久久久免| 一区二区三区精品91| 免费观看的影片在线观看| 国产精品99久久久久久久久| 丰满迷人的少妇在线观看| 欧美一级a爱片免费观看看| 又大又黄又爽视频免费| 国产v大片淫在线免费观看| 亚洲精品亚洲一区二区| 最近2019中文字幕mv第一页| 交换朋友夫妻互换小说| 女性被躁到高潮视频| 久久久久久久久久人人人人人人| 国产成人a∨麻豆精品| 欧美激情国产日韩精品一区| 偷拍熟女少妇极品色| 亚洲丝袜综合中文字幕| 九色成人免费人妻av| 最近最新中文字幕大全电影3| 日韩伦理黄色片| 妹子高潮喷水视频| 亚洲精品久久久久久婷婷小说| 在线免费十八禁| 国产熟女欧美一区二区| 国语对白做爰xxxⅹ性视频网站| 亚洲中文av在线| 一级二级三级毛片免费看| 亚洲av中文字字幕乱码综合| 欧美另类一区| 赤兔流量卡办理| 国产精品一及| 欧美xxxx性猛交bbbb| 国产在线视频一区二区| 午夜福利在线在线| 我要看日韩黄色一级片| 99热这里只有是精品在线观看| kizo精华| 99久久精品热视频| 成人特级av手机在线观看| 国产亚洲一区二区精品| 久久久久久久久久久免费av| 久久人人爽人人爽人人片va| 国产一级毛片在线| 99视频精品全部免费 在线| 多毛熟女@视频| 熟妇人妻不卡中文字幕| 日韩欧美 国产精品| 色视频在线一区二区三区| 极品少妇高潮喷水抽搐| 汤姆久久久久久久影院中文字幕| 国产成人a∨麻豆精品| 精品人妻一区二区三区麻豆| 91aial.com中文字幕在线观看| 久久久久性生活片| 亚洲欧美精品专区久久| 欧美3d第一页| 国产 一区精品| 99热这里只有是精品在线观看| 亚洲一区二区三区欧美精品| 国产乱来视频区| 国产乱人视频| 久久精品国产亚洲网站| 日韩一本色道免费dvd| 精品久久久久久久末码| kizo精华| 美女视频免费永久观看网站| 日本欧美视频一区| 国产精品一区二区性色av| 高清av免费在线| 插阴视频在线观看视频| 一级二级三级毛片免费看| 2021少妇久久久久久久久久久| 少妇人妻精品综合一区二区| 国产成人精品福利久久| 少妇人妻久久综合中文| 国产中年淑女户外野战色| 欧美最新免费一区二区三区| 成人特级av手机在线观看| 亚洲美女搞黄在线观看| 欧美精品人与动牲交sv欧美| 国产精品久久久久久av不卡| 久久久欧美国产精品| 国产伦精品一区二区三区视频9| 亚洲欧美中文字幕日韩二区| 久久久久久久久久成人| 日韩不卡一区二区三区视频在线| 日本一二三区视频观看| 国产成人午夜福利电影在线观看| 国产一区二区三区av在线| 中文天堂在线官网| 啦啦啦在线观看免费高清www| 国产精品一二三区在线看| 国产av码专区亚洲av| 超碰av人人做人人爽久久| 国产成人a区在线观看| 欧美bdsm另类| 少妇人妻精品综合一区二区| 熟女电影av网| 久久青草综合色| 国产午夜精品一二区理论片| 亚洲av.av天堂| 一级av片app| 久久热精品热| 最近手机中文字幕大全| 国产成人一区二区在线| 日韩中文字幕视频在线看片 | 边亲边吃奶的免费视频| av卡一久久| 色5月婷婷丁香| 一级片'在线观看视频| 97在线人人人人妻| 免费不卡的大黄色大毛片视频在线观看| tube8黄色片| 亚洲av福利一区| 日韩一区二区视频免费看| 国产国拍精品亚洲av在线观看| 欧美高清性xxxxhd video| 国产亚洲5aaaaa淫片| 欧美zozozo另类| 99久久中文字幕三级久久日本| 五月天丁香电影| 高清视频免费观看一区二区| 在线观看三级黄色| 在线亚洲精品国产二区图片欧美 | 毛片女人毛片| 伊人久久精品亚洲午夜| 国产精品免费大片| 国产高清三级在线| 免费播放大片免费观看视频在线观看| 久久精品国产亚洲网站| 久久久久久久久久久丰满| 在线观看人妻少妇| 国产极品天堂在线| 免费观看的影片在线观看| 十八禁网站网址无遮挡 | 成人国产麻豆网| 建设人人有责人人尽责人人享有的 | 国产免费福利视频在线观看| 国产av国产精品国产| 一区二区三区四区激情视频| 日韩在线高清观看一区二区三区| 最新中文字幕久久久久| 精品人妻熟女av久视频| 波野结衣二区三区在线| 亚洲国产精品一区三区| 成人美女网站在线观看视频| 国产伦精品一区二区三区视频9| 国产69精品久久久久777片| 亚洲av欧美aⅴ国产| 日本欧美国产在线视频| 午夜日本视频在线| 久久久久久久大尺度免费视频| 久久99热这里只有精品18| 国产高清三级在线| a级毛片免费高清观看在线播放| 人妻一区二区av| 女人十人毛片免费观看3o分钟| 色网站视频免费| 男人和女人高潮做爰伦理| 中文天堂在线官网| 国产69精品久久久久777片| 少妇丰满av| 国产大屁股一区二区在线视频| 高清av免费在线| 欧美老熟妇乱子伦牲交| 毛片一级片免费看久久久久| av.在线天堂| 国产美女午夜福利| 日韩中字成人| 久久久a久久爽久久v久久| 精品久久久久久久久av| 欧美人与善性xxx| 欧美成人a在线观看| 永久网站在线| 青春草国产在线视频| 26uuu在线亚洲综合色| 在线观看人妻少妇| 伦理电影免费视频| 性高湖久久久久久久久免费观看| 国产精品熟女久久久久浪| 少妇人妻精品综合一区二区| 国产精品福利在线免费观看| 直男gayav资源| av天堂中文字幕网| 日本一二三区视频观看| 国产极品天堂在线| 亚洲精品成人av观看孕妇| 久久99精品国语久久久| 99热这里只有是精品在线观看| 久久99精品国语久久久| 一区二区三区免费毛片| 亚洲av综合色区一区| 亚洲精品视频女| 51国产日韩欧美| 男人舔奶头视频| 狂野欧美白嫩少妇大欣赏| 精品人妻偷拍中文字幕| 日韩 亚洲 欧美在线| 国产女主播在线喷水免费视频网站| 日本一二三区视频观看| 永久网站在线| 国产淫片久久久久久久久| 日韩一本色道免费dvd| 国产视频内射| 日本午夜av视频| 国产精品蜜桃在线观看| 又大又黄又爽视频免费| 亚洲美女搞黄在线观看| 久久久久精品久久久久真实原创| 成年av动漫网址| 久久精品国产自在天天线| 精品99又大又爽又粗少妇毛片| 亚洲欧美日韩无卡精品| 妹子高潮喷水视频| 国产精品无大码| 丝袜脚勾引网站| 亚洲精品久久久久久婷婷小说| 亚洲精品第二区| 国产极品天堂在线| 97精品久久久久久久久久精品| 国产女主播在线喷水免费视频网站| 最近中文字幕高清免费大全6| 国产精品国产av在线观看| 亚洲国产精品专区欧美| 日日啪夜夜爽| 国产真实伦视频高清在线观看| 国产精品一二三区在线看| 亚洲激情五月婷婷啪啪| 人妻少妇偷人精品九色| 亚洲电影在线观看av| 97热精品久久久久久| av福利片在线观看| 久久午夜福利片| 女人久久www免费人成看片| 天美传媒精品一区二区| 久久久欧美国产精品| 久久人人爽人人爽人人片va| 国产色爽女视频免费观看| 99久久精品热视频| 久久久久久久久久久丰满| 日本vs欧美在线观看视频 | 国产成人午夜福利电影在线观看| 欧美成人a在线观看| 中文字幕av成人在线电影| 亚洲人成网站在线播| 国产精品三级大全| 成人免费观看视频高清| 乱系列少妇在线播放| 日韩精品有码人妻一区| h视频一区二区三区| 一二三四中文在线观看免费高清| 寂寞人妻少妇视频99o| 亚洲av综合色区一区| 18禁裸乳无遮挡动漫免费视频| 91精品一卡2卡3卡4卡| 毛片一级片免费看久久久久| 22中文网久久字幕| 成年av动漫网址| 麻豆国产97在线/欧美| 人妻制服诱惑在线中文字幕| 国产老妇伦熟女老妇高清| 3wmmmm亚洲av在线观看| 男女无遮挡免费网站观看| 啦啦啦视频在线资源免费观看| 亚洲国产日韩一区二区| 嫩草影院新地址| 大话2 男鬼变身卡| 18禁动态无遮挡网站| 亚洲内射少妇av| 看十八女毛片水多多多| 乱系列少妇在线播放| 草草在线视频免费看| 80岁老熟妇乱子伦牲交| 18禁在线无遮挡免费观看视频| 国产黄色免费在线视频| 精品亚洲成a人片在线观看 | 国产熟女欧美一区二区| 国产69精品久久久久777片| 五月玫瑰六月丁香| www.av在线官网国产| 国产成人aa在线观看| 少妇人妻久久综合中文| 色5月婷婷丁香| 热99国产精品久久久久久7| 黄片wwwwww| 另类亚洲欧美激情| 97热精品久久久久久| 黄片无遮挡物在线观看| 中文字幕久久专区| 91精品国产国语对白视频| 日韩av在线免费看完整版不卡| 国产乱人偷精品视频| 在线观看国产h片| 九九在线视频观看精品| 男人狂女人下面高潮的视频| 麻豆国产97在线/欧美| 国产男女超爽视频在线观看| 欧美3d第一页| 九色成人免费人妻av| 久久女婷五月综合色啪小说| 天堂俺去俺来也www色官网| xxx大片免费视频| 亚洲色图av天堂| 麻豆精品久久久久久蜜桃| 精品午夜福利在线看| 女性生殖器流出的白浆| 三级国产精品片| 在线观看三级黄色| 三级经典国产精品| 99热全是精品| xxx大片免费视频| 国产黄片美女视频| 性色av一级| 亚洲内射少妇av| 国产精品一区二区在线观看99| 性色av一级| 一级a做视频免费观看| 最近中文字幕2019免费版| 嘟嘟电影网在线观看| 成人特级av手机在线观看| 亚洲电影在线观看av| 国产91av在线免费观看| 激情五月婷婷亚洲| 你懂的网址亚洲精品在线观看| 91精品国产九色| 国产亚洲最大av| 久久人妻熟女aⅴ| 高清毛片免费看| 青青草视频在线视频观看| 国产爱豆传媒在线观看| 国产 一区 欧美 日韩| a级毛片免费高清观看在线播放| 亚洲成色77777|