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

    調(diào)水期風(fēng)場(chǎng)對(duì)淮河入江水道淺水湖泊水動(dòng)力水質(zhì)的影響研究

    2023-06-05 13:11:26饒貴康王玲玲胡禮涵徐鵬程唐洪武
    關(guān)鍵詞:兩湖高郵湖區(qū)

    饒貴康,王玲玲,徐 津,胡禮涵,徐鵬程,唐洪武

    (1. 河海大學(xué)水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210098; 2. 河海大學(xué)水利水電學(xué)院,江蘇 南京 210098)

    我國平原河網(wǎng)地區(qū)地勢(shì)低洼,淺水湖泊分布眾多[1-2],大多兼有防洪供水、生態(tài)修復(fù)、灌溉航運(yùn)等功能,同時(shí)還具有平原水庫的調(diào)蓄功效。風(fēng)場(chǎng)是影響淺水湖泊湖流運(yùn)動(dòng)、泥沙沖淤、污染物輸移的主要?jiǎng)恿χ籟3-5],對(duì)塑造大型淺水湖泊水動(dòng)力場(chǎng)的形態(tài)和結(jié)構(gòu)有重要作用,因此研究淺水湖泊的風(fēng)生流問題對(duì)于湖泊資源的開發(fā)利用與保護(hù)具有重要意義。

    目前湖泊風(fēng)生流研究主要有現(xiàn)場(chǎng)觀測(cè)、實(shí)驗(yàn)室物理模型以及數(shù)值模擬等方法。數(shù)值模擬由于可以全方位、高效率地模擬湖泊水動(dòng)力在不同水文、氣象和地形要素驅(qū)動(dòng)下的時(shí)空演化特征,現(xiàn)已成為風(fēng)生流研究的重要手段[3]。許旭峰等[6]采用平面二維數(shù)學(xué)模型,研究了太湖夏季盛行風(fēng)對(duì)湖區(qū)水位、流速以及流態(tài)的影響;黨晚婷[7]采用平面二維水動(dòng)力水質(zhì)模型,研究了風(fēng)力驅(qū)動(dòng)下不同形態(tài)的封閉湖灣與湖區(qū)水體置換能力的響應(yīng)關(guān)系;陳志琦等[8]為準(zhǔn)確分析羅源灣水體的交換周期以及納潮量的大小,數(shù)值計(jì)算中施加了當(dāng)?shù)仫L(fēng)荷載對(duì)湖灣流場(chǎng)的影響;烏景秀等[9]通過對(duì)比風(fēng)向、湖泊地形、補(bǔ)水量等因素對(duì)靜水湖泊水流流態(tài)的影響,提出了改善湖泊動(dòng)力環(huán)境的調(diào)控方案;Liu等[10]基于二維水流-污染物耦合數(shù)值模型,分析了風(fēng)場(chǎng)驅(qū)動(dòng)下太湖環(huán)流及污染物的輸移過程;Jalil等[11]通過模擬風(fēng)場(chǎng)引起的水動(dòng)力變化對(duì)淺水湖泊水流結(jié)構(gòu)及沉積物再懸浮的影響,揭示了風(fēng)荷載在沉積物再懸浮中的作用機(jī)制;舒葉華等[12]利用水質(zhì)模型,模擬了風(fēng)生流驅(qū)動(dòng)下太湖湖區(qū)污染物的輸移規(guī)律,指出環(huán)流結(jié)構(gòu)是導(dǎo)致太湖北部湖灣污染物長(zhǎng)期聚集的主要原因。

    前人的研究成果主要集中在探討風(fēng)場(chǎng)對(duì)湖流結(jié)構(gòu)及污染物輸移特性的影響方面,而風(fēng)場(chǎng)對(duì)調(diào)蓄湖泊沿線引調(diào)水工程設(shè)計(jì)運(yùn)行影響方面的研究成果少有報(bào)道。高郵湖、邵伯湖位于淮河入江水道中段[13],是南水北調(diào)后續(xù)工程兩個(gè)重要的調(diào)蓄湖泊。按照規(guī)劃的要求,南水北調(diào)后續(xù)工程運(yùn)行時(shí),兩湖的湖流流向?qū)⒂稍瓉淼挠杀毕蚰细臑橛赡舷虮?兩湖的淮河水質(zhì)將被長(zhǎng)江水質(zhì)置換。由于高郵湖、邵伯湖是典型的淺水湖泊,風(fēng)生流對(duì)湖區(qū)流場(chǎng)以及水體置換作用不容忽視,風(fēng)涌增減水作用還可能影響兩湖之間抽水泵站提水功能的發(fā)揮。鑒于此,本文以高郵湖、邵伯湖為研究對(duì)象,建立二維水動(dòng)力水質(zhì)模型,分析風(fēng)場(chǎng)與調(diào)水動(dòng)力共同作用下兩湖的流場(chǎng)特征以及物質(zhì)輸移規(guī)律,研究成果可為工程方案的設(shè)計(jì)優(yōu)化、湖區(qū)水質(zhì)水生態(tài)的保護(hù)提供參考。

    1 研究區(qū)域概況

    高郵湖、邵伯湖位于蘇皖兩省交界處,是淮河入江水道上重要的蓄洪湖泊和行洪通道。兩湖屬于過水型湖泊,水域總面積分別約760.67km2和98km2,平均水深約1.44m和1.1m,周圍地形以平原為主,常年主導(dǎo)風(fēng)向?yàn)闁|南風(fēng)[14-16],冬季會(huì)刮西北風(fēng)和東北風(fēng),平均風(fēng)速約3.6m/s。湖區(qū)內(nèi)灘地密布,水利工程眾多,其中新民灘的7個(gè)漫水閘與高程較高的灘地組成了新民灘控制線,為兩湖的分界線(圖1)。南水北調(diào)后續(xù)工程規(guī)劃的全湖方案將在湖區(qū)開展抽槽工程,抽槽底寬邵伯湖段為300m,高郵湖段為200m。兩湖有多個(gè)濕地及生態(tài)保護(hù)區(qū),輸水槽道設(shè)計(jì)線位避開了這些生態(tài)敏感區(qū)域。擬建的高郵泵站緊鄰新民灘控制線的西岸堤防,共安裝8臺(tái)機(jī)組,設(shè)計(jì)揚(yáng)程1.3m。調(diào)水期間,高郵泵站在保證邵伯湖生態(tài)水位的條件下,抽引370m3/s流量進(jìn)入高郵湖。

    圖1 研究區(qū)域范圍

    2 模型建立

    2.1 模型原理

    MIKE21模型有穩(wěn)定、快速、可靠的運(yùn)行特點(diǎn),其計(jì)算原理是利用以單元為中心的有限體積法對(duì)水流連續(xù)方程和動(dòng)量方程進(jìn)行離散求解,可模擬各種作用力下二維自由表面流的水位及流速變化,對(duì)淺水湖泊風(fēng)生流的模擬效果較好,在太湖、巢湖、洪澤湖等流域已有諸多應(yīng)用。綜合考慮模型適用性、現(xiàn)有資料情況,采用該模型建立高郵湖、邵伯湖二維水動(dòng)力水質(zhì)模型。

    2.1.1 水動(dòng)力模型基本原理

    水深平均二維水動(dòng)力控制方程由連續(xù)性方程(式(1))和動(dòng)量方程(式(2))組成。

    (1)

    (2)

    (3)

    式中:t為時(shí)間;x、y、z為笛卡爾坐標(biāo);g為重力加速度;h為總水深;η為自由水面高程;d為靜止水深;u、v為x、y方向上的流速分量;f為科氏力參數(shù);ρ為水的密度;pa為當(dāng)?shù)卮髿鈮簭?qiáng);ρ0為水的相對(duì)密度;τsx、τsy為風(fēng)場(chǎng)摩擦力在x、y方向上的分量;τbx、τby為河床摩擦力在x、y方向的分量;sxx、sxy、syy、syx為輻射應(yīng)力分量;Txx、Txy、Tyy、Tyx為水平黏滯應(yīng)力項(xiàng);S為點(diǎn)源流量大小;us、vs為源匯項(xiàng)水流流速。

    2.1.2 水質(zhì)模型基本原理

    水質(zhì)指標(biāo)控制方程為對(duì)流擴(kuò)散方程,考慮源項(xiàng)與降解過程:

    (4)

    式中:C為垂線平均濃度;Dx、Dy為x、y方向的擴(kuò)散系數(shù);Kc為綜合降解系數(shù);Sa為源匯項(xiàng)。

    2.2 模型范圍及網(wǎng)格劃分

    模型采用兩湖最新的實(shí)測(cè)地形資料,利用AutoCAD提取高程點(diǎn)并生成用于模擬的地形文件。模型范圍包括高郵湖區(qū)、邵伯湖區(qū)、新民灘灘地與邵伯湖群灘。輸水槽道采用四邊形網(wǎng)格模擬,其余區(qū)域采用三角形網(wǎng)格,在湖區(qū)出入口、輸水槽道以及灘地區(qū)域進(jìn)行局部加密,空間步長(zhǎng)10~500 m。計(jì)算節(jié)點(diǎn)共7 715個(gè),計(jì)算單元11 158個(gè),網(wǎng)格剖分示意圖如圖2所示。

    圖2 二維模型網(wǎng)格劃分

    2.3 參數(shù)選取與邊界設(shè)置

    為滿足模型穩(wěn)定與計(jì)算精度要求,同時(shí)受限于柯朗條件[17],模型計(jì)算的時(shí)間步長(zhǎng)取120s;干濕邊界的水深系數(shù)分別取模型中的默認(rèn)參數(shù):干水深hdry= 0.005m,淹沒水深hflood= 0.05m,濕水深hwet= 0.1m;風(fēng)阻力系數(shù)依據(jù)風(fēng)速值的大小由模型自動(dòng)插值計(jì)算;渦黏系數(shù)取0.28。水質(zhì)模塊的參數(shù)主要包括綜合降解系數(shù)與擴(kuò)散系數(shù),根據(jù)前人的相關(guān)研究[18-20],高郵湖、邵伯湖總磷的降解系數(shù)可取0.04d-1,擴(kuò)散系數(shù)[21]取1。

    模型的邊界條件分別設(shè)置于邵伯湖南端和高郵湖的出湖口(施尖)。調(diào)水時(shí)邵伯湖南端為模型的入流邊界,給定調(diào)水流量為370m3/s;高郵湖的出湖口(施尖)為模型的出流邊界,給定設(shè)計(jì)水位5.13m;兩湖間高郵泵站的調(diào)水作用概化為模型的內(nèi)邊界,給定設(shè)計(jì)調(diào)水流量370m3/s。上述水位流量數(shù)據(jù)均來源于工程設(shè)計(jì)單位提供的資料。

    3 模型率定

    2007年淮河流域自北向南大部分地區(qū)先后降暴雨到大暴雨,降雨過程歷時(shí)長(zhǎng)、范圍廣,為緩解淮河中游的防洪壓力,入江水道沿線水閘全開泄水。由于該次洪水來流量大,洪峰期水位變化明顯,適合作為數(shù)學(xué)模型的率定資料。選取2007年7月9日至8月30日入江水道實(shí)測(cè)的洪水過程對(duì)模型進(jìn)行率定,最終確定高郵湖、邵伯湖的糙率分別為0.021、0.027。湖區(qū)內(nèi)兩個(gè)水位站點(diǎn)高郵(高)與六閘的水位驗(yàn)證結(jié)果如圖3所示,由圖可知模型計(jì)算的水位過程與實(shí)測(cè)過程吻合較好,納什效率系數(shù)分別達(dá)到了0.971、0.988,峰值水位誤差小于3%,表明所建模型精度較高,糙率取值合理,可用于后續(xù)調(diào)水工程的研究。

    4 計(jì)算工況及結(jié)果分析

    高郵湖、邵伯湖具有水深淺、湖面寬等特點(diǎn),受外界因素影響水動(dòng)力過程復(fù)雜多變。風(fēng)是影響淺水湖泊湖流運(yùn)動(dòng)的主要驅(qū)動(dòng)力之一。為研究盛行風(fēng)在調(diào)水期對(duì)高郵湖、邵伯湖輸水流場(chǎng)的影響,選用東南風(fēng)、西北風(fēng)、東北風(fēng)3個(gè)盛行方向的定常風(fēng)進(jìn)行數(shù)值試驗(yàn),風(fēng)速值取為該地區(qū)多年平均風(fēng)速3.6m/s,同時(shí)考慮風(fēng)速增高至7.2m/s時(shí)的影響,工況設(shè)置見表1,計(jì)算時(shí)初始條件分別給定為高郵湖、邵伯湖的設(shè)計(jì)輸水水位5.13m與4.13m。

    表1 計(jì)算工況設(shè)置

    4.1 風(fēng)場(chǎng)作用對(duì)湖區(qū)流場(chǎng)的影響

    圖4是7種風(fēng)況下高郵湖、邵伯湖的風(fēng)生流場(chǎng)圖,圖中黑色直線表示新民灘控制線。由圖4可知,無風(fēng)時(shí)高郵湖、邵伯湖僅受泵站調(diào)水動(dòng)力的影響,水流從邵伯湖南端入湖后一路向北流動(dòng),整個(gè)過程并未產(chǎn)生環(huán)流或漩渦等水動(dòng)力現(xiàn)象,湖區(qū)流態(tài)穩(wěn)定。

    圖4 不同風(fēng)況下兩湖區(qū)輸水流場(chǎng)

    當(dāng)受風(fēng)場(chǎng)影響時(shí),不同風(fēng)況下湖區(qū)北部會(huì)形成不同結(jié)構(gòu)的風(fēng)生環(huán)流,環(huán)流的流向、形狀與風(fēng)場(chǎng)風(fēng)向相關(guān)。隨著風(fēng)速增大,兩湖的水動(dòng)力強(qiáng)度提升,流速增幅為0.01~0.06m/s,但流場(chǎng)的總體形態(tài)改變不大,僅西北風(fēng)作用時(shí)高郵湖出現(xiàn)環(huán)流增多的現(xiàn)象。對(duì)比風(fēng)場(chǎng)作用前后兩湖區(qū)的流場(chǎng),可以發(fā)現(xiàn)風(fēng)荷載對(duì)高郵湖的影響程度明顯強(qiáng)于邵伯湖。由于邵伯湖為典型的河道型湖泊,狹長(zhǎng)的湖泊地形削弱了風(fēng)對(duì)湖面的剪切拖曳作用;而高郵湖具有平面尺度大、湖區(qū)水深淺等特點(diǎn),水動(dòng)力較弱的水域極易受風(fēng)荷載的影響,因此湖區(qū)北部流場(chǎng)復(fù)雜多變,在調(diào)水期可能會(huì)對(duì)湖區(qū)的物質(zhì)輸運(yùn)與水質(zhì)更新造成影響。

    4.2 風(fēng)場(chǎng)作用對(duì)高郵泵站樞紐運(yùn)行水位的影響

    風(fēng)涌增減水現(xiàn)象是指風(fēng)生湖流向前運(yùn)動(dòng)時(shí),若遇障礙物阻擋會(huì)在迎風(fēng)岸引起水位抬高、背風(fēng)岸引起水位下降的現(xiàn)象。由于高郵站連同新民灘控制線將高郵湖和邵伯湖隔離成上下兩個(gè)梯級(jí)湖泊,阻擋了風(fēng)生流的前行方向,導(dǎo)致高郵站的上下游易發(fā)風(fēng)涌增減水現(xiàn)象,影響泵站工程的正常運(yùn)行。

    為研究調(diào)水期風(fēng)荷載對(duì)高郵湖、邵伯湖水位的影響,繪制工況2~7與工況1的水位差等值線如圖5所示。由圖5可知,風(fēng)場(chǎng)作用會(huì)引起兩湖區(qū)的水位變化,沿吹程方向水位差呈梯級(jí)分布。東南風(fēng)時(shí),由于新民灘控制線阻擋了風(fēng)生流的前行方向,水流會(huì)在邵伯湖一側(cè)形成水位壅高,高郵湖一側(cè)形成水位下降,且風(fēng)速越高水位變幅越大;西北風(fēng)時(shí)兩湖區(qū)的水位差分布同東南風(fēng)相似,但趨勢(shì)相反,即靠近新民灘一側(cè)高郵湖水位壅高,而邵伯湖一側(cè)水位則下降;東北風(fēng)作用時(shí)整個(gè)湖面向西南側(cè)傾斜,兩湖區(qū)的迎風(fēng)岸水位壅高,而背風(fēng)岸水位下降。

    圖5 兩湖區(qū)水位差等值線(單位:m)

    6個(gè)工況的計(jì)算結(jié)果表明,新民灘控制線的附近湖面是風(fēng)涌水現(xiàn)象最為顯著的區(qū)域,風(fēng)速越高該區(qū)域的水位變幅越大,上下級(jí)湖面的水位差最大可超過0.2m。倘若高郵泵站在設(shè)計(jì)之初沒有考慮風(fēng)涌水的影響,可能導(dǎo)致設(shè)計(jì)揚(yáng)程偏低,調(diào)水動(dòng)力不足。此外,湖區(qū)水位過低也可能影響泵站取水口的取水效果。因此,建于湖泊沿岸的泵站工程、取水工程等水利設(shè)施,在規(guī)劃建設(shè)時(shí)有必要考慮風(fēng)涌水的影響,避免湖區(qū)水位變化過大而影響工程的正常運(yùn)行。

    4.3 風(fēng)場(chǎng)作用對(duì)兩湖區(qū)水體更新的影響

    污染物輸移擴(kuò)散與水體的更新置換是南水北調(diào)后續(xù)工程水質(zhì)水生態(tài)研究的一項(xiàng)重要內(nèi)容。長(zhǎng)時(shí)間的輸水過程中,湖泊水體將與長(zhǎng)江水體發(fā)生融合、置換,從而引起湖區(qū)水質(zhì)的變化。

    以高郵湖、邵伯湖的污染物總氮為研究對(duì)象,分析調(diào)水期不同風(fēng)場(chǎng)對(duì)兩湖區(qū)水質(zhì)分布與水體交換能力的影響。根據(jù)2019年中國科學(xué)院南京地理與湖泊研究所的現(xiàn)場(chǎng)實(shí)測(cè),高郵湖總氮質(zhì)量濃度在0.85~2.09mg/L之間,高濃度水域集中在湖區(qū)的西南區(qū)域;邵伯湖總氮質(zhì)量濃度在1.15~1.48mg/L之間;長(zhǎng)江干流的三江營(yíng)站點(diǎn)總氮質(zhì)量濃度為2mg/L左右,高于兩湖。因此從總氮來看,兩湖水質(zhì)優(yōu)于長(zhǎng)江水質(zhì)。

    根據(jù)實(shí)測(cè)資料,設(shè)定模型計(jì)算時(shí)高郵湖與邵伯湖的總氮初始質(zhì)量濃度如圖6所示,圖中A、B、C、D是4個(gè)水質(zhì)監(jiān)測(cè)點(diǎn)。模擬時(shí)邵伯湖南部邊界的總氮入流質(zhì)量濃度設(shè)定為2mg/L。經(jīng)計(jì)算,當(dāng)在無風(fēng)條件下(工況1)調(diào)水約90d,4個(gè)監(jiān)測(cè)點(diǎn)的總氮質(zhì)量濃度基本穩(wěn)定并保持不變,如圖7所示。

    圖6 高郵湖、邵伯湖總氮初始質(zhì)量濃度

    圖7 無風(fēng)條件下4個(gè)水質(zhì)監(jiān)測(cè)點(diǎn)的總氮質(zhì)量濃度變化

    為了研究風(fēng)荷載對(duì)兩湖區(qū)水體更新置換的影響,圖8比較了不同風(fēng)況下兩湖區(qū)調(diào)水90d后總氮質(zhì)量濃度分布結(jié)果。由圖8可知,無風(fēng)條件下,兩湖大部分區(qū)域的總氮質(zhì)量濃度已基本升至調(diào)水濃度,其中高郵湖、邵伯湖生態(tài)保護(hù)區(qū)均被長(zhǎng)江水體置換,僅高郵湖北端及湖區(qū)西南側(cè)水動(dòng)力較弱的區(qū)域總氮質(zhì)量濃度還維持初始較低值。當(dāng)風(fēng)向?yàn)闁|南風(fēng)時(shí)(圖8(b)(e)),兩湖水體的置換程度明顯提高,高郵湖北部原本難以更新的水域在東南風(fēng)的驅(qū)動(dòng)下也基本置換為長(zhǎng)江水體。相比于無風(fēng)工況,東南風(fēng)的作用能夠顯著提升湖區(qū)內(nèi)污染物的輸移能力。而當(dāng)湖區(qū)受冬季西北風(fēng)影響時(shí)(圖8(c)(f)),由于風(fēng)向與調(diào)水流向相反,逆風(fēng)作用使調(diào)水水體難以運(yùn)動(dòng)至湖區(qū)北部,加之湖區(qū)北部形成的風(fēng)生環(huán)流在很大程度上抑制了兩種水體間的融合交換,因此輸水過程中該水域的總氮質(zhì)量濃度始終較低,水體更新程度弱。東北風(fēng)也屬于逆風(fēng)類工況(圖8(d)(g)),調(diào)水期高郵湖的水體更新同樣受阻,但相比于西北風(fēng)其影響程度相對(duì)較低,水體更新有所增強(qiáng)。

    圖8 調(diào)水第90d時(shí)兩湖總氮質(zhì)量濃度分布

    上述分析表明,調(diào)水期風(fēng)場(chǎng)是影響淮河入江水道高郵湖水體更新與水質(zhì)分布的重要因素,風(fēng)場(chǎng)條件不同,湖區(qū)內(nèi)總氮等污染物的輸移擴(kuò)散規(guī)律也不同。此外,模擬結(jié)果還表明:高郵湖、邵伯湖生態(tài)保護(hù)區(qū)在調(diào)水時(shí)大部分區(qū)域會(huì)被長(zhǎng)江水體置換,水質(zhì)的改變可能影響原有的水生態(tài)環(huán)境,不利于保護(hù)區(qū)內(nèi)水生生物的生長(zhǎng),有必要采取相關(guān)的應(yīng)對(duì)保護(hù)措施,如加強(qiáng)保護(hù)區(qū)內(nèi)的水質(zhì)監(jiān)測(cè)和管理、設(shè)置有效隔離帶、建立人工濕地、利用水生植物的自然凈化作用降低來流水體的污染物含量等。

    5 結(jié) 論

    a.高郵湖是風(fēng)場(chǎng)作用的敏感區(qū)域,在風(fēng)力驅(qū)動(dòng)下湖區(qū)將形成不同結(jié)構(gòu)的風(fēng)生環(huán)流,環(huán)流的流向、形狀由風(fēng)場(chǎng)風(fēng)向決定;邵伯湖由于平面尺度較小,風(fēng)荷載不會(huì)對(duì)湖區(qū)的流場(chǎng)結(jié)構(gòu)產(chǎn)生過大的影響。

    b.風(fēng)場(chǎng)作用會(huì)引起湖區(qū)迎風(fēng)岸水位抬高、背風(fēng)岸水位下降,新民灘控制線的上下級(jí)湖面是風(fēng)涌增減水最為顯著的區(qū)域,本文研究工況下水位差最大可超過0.2 m。由于風(fēng)荷載對(duì)淺水湖泊的水位影響較大,因此沿湖岸修建的泵站工程、取水工程等水利設(shè)施,在規(guī)劃建設(shè)時(shí)有必要考慮當(dāng)?shù)仫L(fēng)場(chǎng)的影響,以免湖區(qū)水位變化過大而影響涉水工程的正常運(yùn)行。

    c.風(fēng)場(chǎng)是影響淮河入江水道高郵湖水體更新和水質(zhì)分布的重要因素:當(dāng)風(fēng)場(chǎng)風(fēng)向與調(diào)水方向同向時(shí),高郵湖水體的置換程度會(huì)顯著提高;而當(dāng)兩者逆向時(shí),湖區(qū)水體的置換程度將顯著下降。

    d.由于輸水過程中高郵湖、邵伯湖的生態(tài)保護(hù)區(qū)大部分區(qū)域水體將被長(zhǎng)江水體置換,水質(zhì)的改變可能影響原有的水生態(tài)環(huán)境,不利于保護(hù)區(qū)水生生物的生長(zhǎng),應(yīng)采取適當(dāng)?shù)谋Wo(hù)措施。

    猜你喜歡
    兩湖高郵湖區(qū)
    大通湖區(qū)河蟹產(chǎn)業(yè)發(fā)展綜述
    “兩湖”公園空氣負(fù)離子濃度時(shí)空變化特征
    好歌好酒好高郵
    民族音樂(2019年6期)2020-01-15 06:50:18
    生活在湖區(qū)
    海峽旅游(2018年4期)2018-06-01 11:20:00
    湖區(qū)航道風(fēng)速預(yù)警監(jiān)測(cè)點(diǎn)布設(shè)研究
    江西建材(2018年4期)2018-04-10 12:37:24
    看千帆競(jìng)渡,數(shù)高郵風(fēng)流
    ——高郵博物館精華演繹
    大眾考古(2016年5期)2017-01-16 09:32:16
    高郵風(fēng)光
    兩湖文學(xué):作為整體感覺文化區(qū)
    江蘇高郵鴨集團(tuán)
    鎮(zhèn)遠(yuǎn)·兩湖會(huì)館
    国产av国产精品国产| 国产在线视频一区二区| 香蕉久久夜色| 久9热在线精品视频| 久久狼人影院| 两性午夜刺激爽爽歪歪视频在线观看 | 男人舔女人的私密视频| 老司机深夜福利视频在线观看| 国产男女内射视频| 满18在线观看网站| 成人国产一区最新在线观看| 美女福利国产在线| 亚洲成国产人片在线观看| 性高湖久久久久久久久免费观看| 国产视频一区二区在线看| 久久久久精品国产欧美久久久| 亚洲伊人久久精品综合| 天天操日日干夜夜撸| 91成人精品电影| a级毛片黄视频| 天堂8中文在线网| 别揉我奶头~嗯~啊~动态视频| 亚洲中文av在线| avwww免费| 变态另类成人亚洲欧美熟女 | 亚洲av电影在线进入| 搡老岳熟女国产| 国产亚洲一区二区精品| 黑人猛操日本美女一级片| 欧美日韩一级在线毛片| tocl精华| 十八禁人妻一区二区| 18在线观看网站| 久久久久久久精品吃奶| 国产精品.久久久| 99国产精品免费福利视频| 亚洲精品在线观看二区| 妹子高潮喷水视频| 女人爽到高潮嗷嗷叫在线视频| 美女扒开内裤让男人捅视频| 久9热在线精品视频| av网站免费在线观看视频| 欧美日韩国产mv在线观看视频| 亚洲一码二码三码区别大吗| 人人妻人人澡人人爽人人夜夜| 国产午夜精品久久久久久| 亚洲精品久久成人aⅴ小说| 亚洲国产成人一精品久久久| 丝袜喷水一区| 捣出白浆h1v1| 日韩视频一区二区在线观看| 老司机影院毛片| 亚洲一区二区三区欧美精品| 久久这里只有精品19| 亚洲精品国产色婷婷电影| 精品少妇内射三级| 侵犯人妻中文字幕一二三四区| 91精品国产国语对白视频| 免费av中文字幕在线| 在线av久久热| 日本一区二区免费在线视频| 如日韩欧美国产精品一区二区三区| 满18在线观看网站| 免费观看人在逋| 久久天堂一区二区三区四区| 18禁国产床啪视频网站| 80岁老熟妇乱子伦牲交| 欧美在线黄色| 国产在线观看jvid| 午夜精品国产一区二区电影| 丰满饥渴人妻一区二区三| 色播在线永久视频| 丰满人妻熟妇乱又伦精品不卡| 国产精品秋霞免费鲁丝片| a在线观看视频网站| 欧美亚洲 丝袜 人妻 在线| 91国产中文字幕| 啪啪无遮挡十八禁网站| 波多野结衣一区麻豆| 日韩大片免费观看网站| www.熟女人妻精品国产| 亚洲第一欧美日韩一区二区三区 | 欧美久久黑人一区二区| 欧美成狂野欧美在线观看| 成人黄色视频免费在线看| 欧美精品一区二区免费开放| 超色免费av| 欧美精品人与动牲交sv欧美| 国产男女超爽视频在线观看| 国产精品九九99| 五月开心婷婷网| 免费女性裸体啪啪无遮挡网站| 黄色视频,在线免费观看| 国产精品一区二区在线观看99| 欧美日本中文国产一区发布| 91精品三级在线观看| 亚洲成人国产一区在线观看| 亚洲成人国产一区在线观看| 国产精品二区激情视频| 热99久久久久精品小说推荐| 少妇被粗大的猛进出69影院| 丝袜喷水一区| 久久精品aⅴ一区二区三区四区| 午夜福利一区二区在线看| h视频一区二区三区| 国产成人系列免费观看| 国产真人三级小视频在线观看| 国产xxxxx性猛交| 极品少妇高潮喷水抽搐| 免费看a级黄色片| 免费在线观看影片大全网站| 国产日韩欧美视频二区| 一边摸一边做爽爽视频免费| 午夜两性在线视频| 成人特级黄色片久久久久久久 | 妹子高潮喷水视频| 久久婷婷成人综合色麻豆| 12—13女人毛片做爰片一| 老熟妇仑乱视频hdxx| 黄色视频,在线免费观看| 考比视频在线观看| 国产又色又爽无遮挡免费看| 精品人妻1区二区| 日韩 欧美 亚洲 中文字幕| 热99久久久久精品小说推荐| 国产精品九九99| 日韩欧美一区二区三区在线观看 | 美女国产高潮福利片在线看| av电影中文网址| 欧美国产精品va在线观看不卡| 一边摸一边做爽爽视频免费| 成人av一区二区三区在线看| 女人被躁到高潮嗷嗷叫费观| 久久精品国产a三级三级三级| 久久久国产欧美日韩av| 国产日韩欧美在线精品| av视频免费观看在线观看| 亚洲国产欧美日韩在线播放| 变态另类成人亚洲欧美熟女 | 热re99久久精品国产66热6| 婷婷成人精品国产| 丝瓜视频免费看黄片| 久久久久久久大尺度免费视频| 美女福利国产在线| 亚洲欧美日韩高清在线视频 | 免费在线观看影片大全网站| 激情视频va一区二区三区| 老汉色∧v一级毛片| 自线自在国产av| 不卡一级毛片| 丰满饥渴人妻一区二区三| 悠悠久久av| 999久久久国产精品视频| 一级a爱视频在线免费观看| 久热爱精品视频在线9| 国产不卡一卡二| 亚洲av欧美aⅴ国产| 高清在线国产一区| 国产av国产精品国产| 999久久久国产精品视频| 色94色欧美一区二区| 国产麻豆69| 精品高清国产在线一区| 男女午夜视频在线观看| 国产一区二区激情短视频| 日本一区二区免费在线视频| 王馨瑶露胸无遮挡在线观看| 韩国精品一区二区三区| 在线观看66精品国产| 亚洲伊人色综图| 国产精品自产拍在线观看55亚洲 | 亚洲精品美女久久av网站| 亚洲精品国产精品久久久不卡| 欧美精品一区二区大全| 午夜福利在线观看吧| 精品久久蜜臀av无| 啦啦啦视频在线资源免费观看| 757午夜福利合集在线观看| 又紧又爽又黄一区二区| 美女午夜性视频免费| 日韩欧美免费精品| 最新在线观看一区二区三区| 国产极品粉嫩免费观看在线| 久久影院123| 一区二区三区激情视频| 深夜精品福利| 亚洲av欧美aⅴ国产| 黄色丝袜av网址大全| 色老头精品视频在线观看| 久久久精品区二区三区| 欧美激情 高清一区二区三区| 纵有疾风起免费观看全集完整版| 亚洲国产av新网站| 91大片在线观看| 日日夜夜操网爽| 99九九在线精品视频| 国产深夜福利视频在线观看| 久久久水蜜桃国产精品网| 黄色视频,在线免费观看| 亚洲黑人精品在线| 18禁裸乳无遮挡动漫免费视频| 法律面前人人平等表现在哪些方面| 色播在线永久视频| 精品国产亚洲在线| 精品一品国产午夜福利视频| 国内毛片毛片毛片毛片毛片| 99香蕉大伊视频| 高潮久久久久久久久久久不卡| 不卡一级毛片| 露出奶头的视频| 久久午夜亚洲精品久久| 99九九在线精品视频| 国产精品香港三级国产av潘金莲| 久久午夜亚洲精品久久| 欧美乱码精品一区二区三区| 日韩视频在线欧美| 国产福利在线免费观看视频| 久久久精品国产亚洲av高清涩受| 亚洲午夜理论影院| 丁香六月欧美| 亚洲欧美一区二区三区久久| 欧美精品亚洲一区二区| av线在线观看网站| 美女午夜性视频免费| 另类亚洲欧美激情| 精品少妇黑人巨大在线播放| 亚洲成av片中文字幕在线观看| 黄色视频,在线免费观看| 一区二区三区精品91| 别揉我奶头~嗯~啊~动态视频| 在线av久久热| 午夜免费鲁丝| 欧美人与性动交α欧美软件| 国产成人欧美| 日韩制服丝袜自拍偷拍| 日韩精品免费视频一区二区三区| 国产在线视频一区二区| 91av网站免费观看| 1024视频免费在线观看| 青草久久国产| 一本—道久久a久久精品蜜桃钙片| 亚洲自偷自拍图片 自拍| 老司机靠b影院| 精品人妻1区二区| 午夜福利,免费看| 精品福利永久在线观看| av网站在线播放免费| 在线观看免费高清a一片| 999久久久精品免费观看国产| 十八禁高潮呻吟视频| 精品人妻在线不人妻| 99精国产麻豆久久婷婷| 十八禁人妻一区二区| 亚洲国产欧美在线一区| 国产精品熟女久久久久浪| 国产av一区二区精品久久| 男人操女人黄网站| 精品少妇久久久久久888优播| 精品人妻1区二区| 日本a在线网址| 久久久久精品国产欧美久久久| 亚洲精品一卡2卡三卡4卡5卡| 天天躁夜夜躁狠狠躁躁| av在线播放免费不卡| 91成年电影在线观看| 十八禁人妻一区二区| 久久精品国产99精品国产亚洲性色 | 最近最新中文字幕大全电影3 | 999久久久精品免费观看国产| 男女午夜视频在线观看| 欧美黄色片欧美黄色片| 久久人妻福利社区极品人妻图片| 高清毛片免费观看视频网站 | 高潮久久久久久久久久久不卡| 两性午夜刺激爽爽歪歪视频在线观看 | 最近最新中文字幕大全电影3 | 国产在线免费精品| 一级片免费观看大全| 国产成人欧美在线观看 | 国产人伦9x9x在线观看| 老司机深夜福利视频在线观看| 亚洲国产欧美日韩在线播放| 亚洲视频免费观看视频| 叶爱在线成人免费视频播放| 亚洲一卡2卡3卡4卡5卡精品中文| 精品第一国产精品| 777久久人妻少妇嫩草av网站| 国产亚洲精品第一综合不卡| 亚洲精品美女久久av网站| 丝袜人妻中文字幕| 人妻久久中文字幕网| 亚洲专区字幕在线| 韩国精品一区二区三区| 亚洲精品美女久久久久99蜜臀| 高清毛片免费观看视频网站 | 午夜激情av网站| 国产日韩欧美在线精品| 国产精品一区二区在线观看99| 成人黄色视频免费在线看| 国产高清国产精品国产三级| 欧美人与性动交α欧美软件| 亚洲av成人不卡在线观看播放网| 亚洲男人天堂网一区| 久久亚洲真实| av又黄又爽大尺度在线免费看| 精品卡一卡二卡四卡免费| 久久婷婷成人综合色麻豆| 热99久久久久精品小说推荐| 男女下面插进去视频免费观看| 操美女的视频在线观看| 啦啦啦中文免费视频观看日本| 国产日韩欧美视频二区| 国产精品久久久人人做人人爽| 母亲3免费完整高清在线观看| 18禁裸乳无遮挡动漫免费视频| 王馨瑶露胸无遮挡在线观看| 欧美日韩亚洲高清精品| 欧美老熟妇乱子伦牲交| 亚洲av欧美aⅴ国产| 一级片'在线观看视频| 高清欧美精品videossex| 日韩欧美国产一区二区入口| 天堂动漫精品| 中文欧美无线码| 黄色丝袜av网址大全| 亚洲专区国产一区二区| av欧美777| 99国产精品一区二区蜜桃av | 制服诱惑二区| 欧美日韩亚洲国产一区二区在线观看 | 国产在线免费精品| 欧美日韩福利视频一区二区| 亚洲情色 制服丝袜| 成年版毛片免费区| 国产精品香港三级国产av潘金莲| 亚洲七黄色美女视频| 大码成人一级视频| 老司机深夜福利视频在线观看| 国产亚洲精品第一综合不卡| 老司机影院毛片| 日本wwww免费看| 国产1区2区3区精品| 丝瓜视频免费看黄片| 国产男女超爽视频在线观看| 国产麻豆69| 9191精品国产免费久久| 亚洲 国产 在线| 成人免费观看视频高清| 69精品国产乱码久久久| 99精品在免费线老司机午夜| 女人被躁到高潮嗷嗷叫费观| 欧美久久黑人一区二区| 欧美日韩av久久| 久久亚洲精品不卡| 满18在线观看网站| 国产男女内射视频| 亚洲专区国产一区二区| 热99久久久久精品小说推荐| 亚洲人成电影观看| 两个人看的免费小视频| 国产无遮挡羞羞视频在线观看| 亚洲成av片中文字幕在线观看| 欧美精品一区二区免费开放| 久久人妻av系列| 亚洲精品在线美女| 亚洲av美国av| 久9热在线精品视频| 日日夜夜操网爽| 国产深夜福利视频在线观看| 国产成人精品在线电影| 亚洲成a人片在线一区二区| 久久精品亚洲精品国产色婷小说| 国产精品影院久久| 免费女性裸体啪啪无遮挡网站| 亚洲情色 制服丝袜| 国产亚洲精品第一综合不卡| 大型av网站在线播放| www.熟女人妻精品国产| 欧美激情久久久久久爽电影 | 久久久国产欧美日韩av| 黄频高清免费视频| 可以免费在线观看a视频的电影网站| 夜夜骑夜夜射夜夜干| 日韩成人在线观看一区二区三区| 亚洲性夜色夜夜综合| 久久久久久亚洲精品国产蜜桃av| 一区二区三区乱码不卡18| 国产在线观看jvid| 最近最新中文字幕大全免费视频| 亚洲精品国产一区二区精华液| 亚洲一区中文字幕在线| av欧美777| 又大又爽又粗| 黄色怎么调成土黄色| 飞空精品影院首页| 成人亚洲精品一区在线观看| svipshipincom国产片| 亚洲精品在线观看二区| 精品国产一区二区三区久久久樱花| 99久久精品国产亚洲精品| 久久精品亚洲精品国产色婷小说| 国产成人系列免费观看| 国产成人精品久久二区二区91| 欧美另类亚洲清纯唯美| 少妇 在线观看| 亚洲男人天堂网一区| 午夜精品久久久久久毛片777| 色综合婷婷激情| 精品视频人人做人人爽| 国产欧美日韩综合在线一区二区| 不卡av一区二区三区| 久久久欧美国产精品| 美女主播在线视频| 我要看黄色一级片免费的| 日韩 欧美 亚洲 中文字幕| a在线观看视频网站| 国精品久久久久久国模美| 好男人电影高清在线观看| 99re6热这里在线精品视频| 精品午夜福利视频在线观看一区 | 久久久欧美国产精品| 女人精品久久久久毛片| 999久久久精品免费观看国产| 中文字幕人妻丝袜制服| xxxhd国产人妻xxx| 日韩欧美三级三区| 久久久精品94久久精品| 美女午夜性视频免费| 国产av国产精品国产| 国产一区二区 视频在线| 日韩成人在线观看一区二区三区| 欧美人与性动交α欧美精品济南到| www日本在线高清视频| 丝瓜视频免费看黄片| 国产免费现黄频在线看| 欧美成人免费av一区二区三区 | 国产欧美日韩一区二区精品| 免费在线观看视频国产中文字幕亚洲| 精品乱码久久久久久99久播| 日韩有码中文字幕| 夜夜骑夜夜射夜夜干| 久久毛片免费看一区二区三区| 交换朋友夫妻互换小说| 久久人妻av系列| 国产深夜福利视频在线观看| 亚洲一码二码三码区别大吗| 欧美久久黑人一区二区| 欧美+亚洲+日韩+国产| 久久久久久久精品吃奶| 成人18禁高潮啪啪吃奶动态图| 天天影视国产精品| 无人区码免费观看不卡 | 香蕉国产在线看| 精品福利观看| 欧美激情高清一区二区三区| 两个人看的免费小视频| 欧美在线黄色| 亚洲免费av在线视频| 免费在线观看影片大全网站| 欧美激情高清一区二区三区| 黄色怎么调成土黄色| 国产伦理片在线播放av一区| 国产成人一区二区三区免费视频网站| 亚洲色图av天堂| 五月开心婷婷网| 久久久久精品国产欧美久久久| 亚洲精品在线观看二区| 国产成人影院久久av| 无限看片的www在线观看| 精品少妇黑人巨大在线播放| 亚洲avbb在线观看| 老熟女久久久| 一级,二级,三级黄色视频| 日日摸夜夜添夜夜添小说| 精品少妇黑人巨大在线播放| 国产欧美日韩一区二区精品| 欧美人与性动交α欧美软件| www.999成人在线观看| 人妻久久中文字幕网| 国产一区有黄有色的免费视频| 亚洲伊人久久精品综合| 少妇粗大呻吟视频| 国产精品秋霞免费鲁丝片| 精品一区二区三区视频在线观看免费 | 波多野结衣av一区二区av| 精品国产一区二区三区久久久樱花| 可以免费在线观看a视频的电影网站| 国产精品久久久久久精品古装| 男女下面插进去视频免费观看| 真人做人爱边吃奶动态| 宅男免费午夜| 免费久久久久久久精品成人欧美视频| 日本五十路高清| 人成视频在线观看免费观看| 午夜91福利影院| 91麻豆精品激情在线观看国产 | 亚洲精品国产区一区二| 久久精品国产亚洲av高清一级| 怎么达到女性高潮| 精品一品国产午夜福利视频| 黄色片一级片一级黄色片| 亚洲欧美日韩另类电影网站| 国产精品av久久久久免费| www.精华液| 国产日韩欧美亚洲二区| 久久午夜亚洲精品久久| 美女国产高潮福利片在线看| 日韩熟女老妇一区二区性免费视频| 美女扒开内裤让男人捅视频| 韩国精品一区二区三区| 高清毛片免费观看视频网站 | 欧美av亚洲av综合av国产av| 别揉我奶头~嗯~啊~动态视频| 水蜜桃什么品种好| 午夜福利在线免费观看网站| 国产在线精品亚洲第一网站| 日本欧美视频一区| 国产麻豆69| 欧美 日韩 精品 国产| 一边摸一边抽搐一进一小说 | 亚洲精品在线美女| 亚洲国产欧美一区二区综合| 国产99久久九九免费精品| 欧美黑人精品巨大| 精品久久蜜臀av无| 国产在线一区二区三区精| 中文欧美无线码| 国产精品99久久99久久久不卡| 国产男女超爽视频在线观看| 黄色a级毛片大全视频| av视频免费观看在线观看| 精品久久蜜臀av无| 黑人巨大精品欧美一区二区蜜桃| 国产黄频视频在线观看| 亚洲精品在线美女| 国产黄频视频在线观看| 麻豆乱淫一区二区| 久久亚洲真实| 亚洲一区中文字幕在线| 久久久久久久精品吃奶| 欧美精品啪啪一区二区三区| 亚洲中文av在线| 日本av手机在线免费观看| 一区二区日韩欧美中文字幕| av又黄又爽大尺度在线免费看| 国产亚洲精品第一综合不卡| 免费不卡黄色视频| 大型黄色视频在线免费观看| 国产欧美日韩一区二区精品| 一区二区三区激情视频| 久久天躁狠狠躁夜夜2o2o| 国产精品久久久久久精品电影小说| 黄频高清免费视频| 香蕉久久夜色| 黄色视频,在线免费观看| 丰满少妇做爰视频| 午夜福利在线免费观看网站| 中文字幕av电影在线播放| 亚洲av美国av| 亚洲欧洲精品一区二区精品久久久| 午夜福利乱码中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 一区福利在线观看| 熟女少妇亚洲综合色aaa.| 国产又爽黄色视频| 精品一区二区三区四区五区乱码| 国产精品98久久久久久宅男小说| 男女下面插进去视频免费观看| 无限看片的www在线观看| 免费高清在线观看日韩| 美国免费a级毛片| 久久国产精品大桥未久av| 国产成人精品久久二区二区免费| 成人亚洲精品一区在线观看| 国产精品电影一区二区三区 | 考比视频在线观看| 国产亚洲欧美精品永久| av欧美777| 正在播放国产对白刺激| 欧美午夜高清在线| 91大片在线观看| 91麻豆精品激情在线观看国产 | 在线永久观看黄色视频| 久9热在线精品视频| 自线自在国产av| 国产精品美女特级片免费视频播放器 | 人妻 亚洲 视频| 午夜免费鲁丝| 国产三级黄色录像| 色播在线永久视频| 黄频高清免费视频| 亚洲人成伊人成综合网2020| 亚洲伊人色综图| 性少妇av在线| 涩涩av久久男人的天堂| 在线av久久热| 亚洲欧美色中文字幕在线| 大码成人一级视频| 日韩 欧美 亚洲 中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 菩萨蛮人人尽说江南好唐韦庄| 国产亚洲av高清不卡| 久久午夜亚洲精品久久| 男人舔女人的私密视频| 国产成人免费观看mmmm| 亚洲国产毛片av蜜桃av| 久久久国产一区二区| 国产不卡av网站在线观看| 免费在线观看日本一区| 亚洲av电影在线进入| 成人永久免费在线观看视频 | 精品高清国产在线一区| 黄频高清免费视频| 在线 av 中文字幕| 久久精品亚洲av国产电影网|