• <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ì)館
    欧美日韩一级在线毛片| 久久99一区二区三区| 在线精品无人区一区二区三| 亚洲成人av在线免费| 国产淫语在线视频| 欧美激情 高清一区二区三区| 在线精品无人区一区二区三| 大话2 男鬼变身卡| 欧美人与性动交α欧美精品济南到 | 亚洲婷婷狠狠爱综合网| 成人18禁高潮啪啪吃奶动态图| √禁漫天堂资源中文www| 成人18禁高潮啪啪吃奶动态图| 久久韩国三级中文字幕| 日韩 亚洲 欧美在线| 天天影视国产精品| 人成视频在线观看免费观看| 色哟哟·www| 欧美 亚洲 国产 日韩一| 久久精品国产亚洲av天美| 欧美日本中文国产一区发布| 欧美人与性动交α欧美精品济南到 | 亚洲成人手机| 日日啪夜夜爽| 狠狠婷婷综合久久久久久88av| 亚洲欧美成人综合另类久久久| 国产成人av激情在线播放| 少妇人妻精品综合一区二区| 一区福利在线观看| 亚洲国产精品一区二区三区在线| 十分钟在线观看高清视频www| 少妇猛男粗大的猛烈进出视频| 999精品在线视频| 久久久国产欧美日韩av| 成年人免费黄色播放视频| 国产男女超爽视频在线观看| 国产激情久久老熟女| 免费久久久久久久精品成人欧美视频| 国产男女内射视频| 最近中文字幕高清免费大全6| 精品少妇内射三级| 欧美成人午夜精品| 欧美在线黄色| 亚洲一区中文字幕在线| 精品卡一卡二卡四卡免费| 中国国产av一级| 日韩av在线免费看完整版不卡| 国产精品二区激情视频| 97人妻天天添夜夜摸| 热re99久久国产66热| 国产极品天堂在线| 久久精品国产亚洲av涩爱| videos熟女内射| 在线 av 中文字幕| 中文乱码字字幕精品一区二区三区| 成年人午夜在线观看视频| 一边摸一边做爽爽视频免费| av线在线观看网站| 久久久a久久爽久久v久久| 国产成人精品一,二区| 久久国产精品男人的天堂亚洲| 亚洲第一av免费看| 18+在线观看网站| 一级爰片在线观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲四区av| 久久精品久久久久久噜噜老黄| 日韩成人av中文字幕在线观看| 国产精品av久久久久免费| 日本91视频免费播放| 日韩 亚洲 欧美在线| av天堂久久9| 日韩不卡一区二区三区视频在线| 午夜福利在线观看免费完整高清在| 国产精品麻豆人妻色哟哟久久| 狂野欧美激情性bbbbbb| 日日摸夜夜添夜夜爱| 成人漫画全彩无遮挡| 亚洲精品av麻豆狂野| 老司机影院成人| 国产日韩欧美视频二区| 男女午夜视频在线观看| 亚洲第一青青草原| 男女啪啪激烈高潮av片| 国产精品亚洲av一区麻豆 | 欧美变态另类bdsm刘玥| 亚洲久久久国产精品| 99国产综合亚洲精品| 亚洲人成77777在线视频| 亚洲av日韩在线播放| 捣出白浆h1v1| 老鸭窝网址在线观看| 69精品国产乱码久久久| 午夜久久久在线观看| 色视频在线一区二区三区| 超碰成人久久| av线在线观看网站| 伦理电影大哥的女人| 国产成人精品无人区| 日韩不卡一区二区三区视频在线| 一级毛片电影观看| 婷婷色综合www| 国产黄色视频一区二区在线观看| 精品人妻一区二区三区麻豆| 看免费成人av毛片| 中文字幕人妻丝袜制服| 欧美亚洲 丝袜 人妻 在线| 国产精品偷伦视频观看了| 爱豆传媒免费全集在线观看| 黄色视频在线播放观看不卡| 爱豆传媒免费全集在线观看| 亚洲av男天堂| av国产久精品久网站免费入址| 精品人妻一区二区三区麻豆| 肉色欧美久久久久久久蜜桃| 免费久久久久久久精品成人欧美视频| 国产精品一二三区在线看| 日本黄色日本黄色录像| 亚洲伊人久久精品综合| 亚洲精品国产av成人精品| tube8黄色片| 亚洲国产精品一区二区三区在线| 伦理电影免费视频| 99久久人妻综合| 最近中文字幕2019免费版| 美女福利国产在线| 欧美精品高潮呻吟av久久| 国产有黄有色有爽视频| www.熟女人妻精品国产| 黑人巨大精品欧美一区二区蜜桃| 母亲3免费完整高清在线观看 | 亚洲三区欧美一区| 岛国毛片在线播放| xxx大片免费视频| 在线免费观看不下载黄p国产| 国产成人一区二区在线| 七月丁香在线播放| 亚洲精品成人av观看孕妇| 人妻一区二区av| 欧美日韩一级在线毛片| 亚洲一区中文字幕在线| 久久久久国产网址| 亚洲综合色网址| 精品少妇一区二区三区视频日本电影 | 久久久久久久国产电影| 国产白丝娇喘喷水9色精品| 国产成人91sexporn| 色播在线永久视频| 久久精品aⅴ一区二区三区四区 | 老司机影院成人| 国产成人午夜福利电影在线观看| 少妇人妻 视频| 99久国产av精品国产电影| 久久精品aⅴ一区二区三区四区 | 啦啦啦啦在线视频资源| 亚洲成人av在线免费| 亚洲av福利一区| 久久久久国产精品人妻一区二区| 大陆偷拍与自拍| 欧美中文综合在线视频| 男人舔女人的私密视频| 天天操日日干夜夜撸| 欧美国产精品一级二级三级| 人妻人人澡人人爽人人| 香蕉精品网在线| 午夜久久久在线观看| 9191精品国产免费久久| 91成人精品电影| 亚洲成人av在线免费| 久久精品久久久久久噜噜老黄| 久久国产亚洲av麻豆专区| 蜜桃国产av成人99| 免费在线观看完整版高清| 丝袜人妻中文字幕| 午夜久久久在线观看| 欧美精品亚洲一区二区| 一级爰片在线观看| 精品国产一区二区三区久久久樱花| 婷婷色综合大香蕉| 国产视频首页在线观看| 少妇猛男粗大的猛烈进出视频| 欧美精品亚洲一区二区| videos熟女内射| 99re6热这里在线精品视频| 久久久久久久久久久久大奶| 男人操女人黄网站| √禁漫天堂资源中文www| 春色校园在线视频观看| 韩国精品一区二区三区| 国产成人一区二区在线| 男女啪啪激烈高潮av片| 一级片'在线观看视频| 嫩草影院入口| 午夜精品国产一区二区电影| 久久亚洲国产成人精品v| 亚洲精品aⅴ在线观看| 成人漫画全彩无遮挡| 亚洲国产精品国产精品| 亚洲人成电影观看| 女性被躁到高潮视频| 亚洲精品美女久久久久99蜜臀 | 国产伦理片在线播放av一区| 可以免费在线观看a视频的电影网站 | 精品99又大又爽又粗少妇毛片| 欧美日韩一区二区视频在线观看视频在线| av天堂久久9| 精品人妻一区二区三区麻豆| 街头女战士在线观看网站| 日韩免费高清中文字幕av| 国产精品国产av在线观看| 免费在线观看黄色视频的| 9热在线视频观看99| 精品酒店卫生间| 最黄视频免费看| 日本av手机在线免费观看| 久久精品国产鲁丝片午夜精品| 麻豆av在线久日| 丝袜美足系列| 97精品久久久久久久久久精品| 丝袜喷水一区| 涩涩av久久男人的天堂| 黄片无遮挡物在线观看| 99久久中文字幕三级久久日本| 精品第一国产精品| 久久久久久免费高清国产稀缺| 99国产精品免费福利视频| 久久99精品国语久久久| 久久午夜综合久久蜜桃| 亚洲精品中文字幕在线视频| 伦理电影大哥的女人| 亚洲欧美成人综合另类久久久| 精品久久久精品久久久| 精品人妻一区二区三区麻豆| 免费女性裸体啪啪无遮挡网站| 啦啦啦在线观看免费高清www| 久久久久网色| 亚洲国产看品久久| 成人毛片a级毛片在线播放| 男女国产视频网站| 妹子高潮喷水视频| 在线精品无人区一区二区三| 婷婷色综合大香蕉| 大片免费播放器 马上看| 午夜影院在线不卡| 边亲边吃奶的免费视频| 日本猛色少妇xxxxx猛交久久| 下体分泌物呈黄色| 国产成人午夜福利电影在线观看| 大香蕉久久成人网| av在线播放精品| 亚洲国产精品国产精品| 免费观看在线日韩| 中国三级夫妇交换| 欧美日韩精品网址| 免费在线观看完整版高清| 欧美日韩精品成人综合77777| 亚洲国产av新网站| 咕卡用的链子| 日韩精品免费视频一区二区三区| 国产熟女欧美一区二区| 一区二区av电影网| 亚洲精品av麻豆狂野| 午夜福利视频在线观看免费| 在线看a的网站| 制服人妻中文乱码| 亚洲欧美日韩另类电影网站| 少妇被粗大的猛进出69影院| 成人手机av| 国产一区二区在线观看av| www日本在线高清视频| 宅男免费午夜| 亚洲综合色惰| 97人妻天天添夜夜摸| 黄片小视频在线播放| 欧美日韩视频精品一区| 一级毛片我不卡| h视频一区二区三区| 国产精品偷伦视频观看了| 午夜老司机福利剧场| 久久国产精品男人的天堂亚洲| 男女啪啪激烈高潮av片| 少妇猛男粗大的猛烈进出视频| 亚洲国产毛片av蜜桃av| 午夜福利网站1000一区二区三区| 国产激情久久老熟女| 一级毛片我不卡| 久久ye,这里只有精品| 飞空精品影院首页| 日本-黄色视频高清免费观看| 黑人猛操日本美女一级片| 久久人人97超碰香蕉20202| 国产女主播在线喷水免费视频网站| 制服人妻中文乱码| tube8黄色片| 亚洲色图 男人天堂 中文字幕| 男男h啪啪无遮挡| 亚洲精品久久成人aⅴ小说| 69精品国产乱码久久久| 亚洲综合精品二区| 男女无遮挡免费网站观看| 一区二区三区精品91| av在线播放精品| 美女中出高潮动态图| 如日韩欧美国产精品一区二区三区| 一本久久精品| 午夜福利,免费看| 精品一区二区三区四区五区乱码 | a级片在线免费高清观看视频| 69精品国产乱码久久久| 亚洲精品日本国产第一区| 宅男免费午夜| 99久久中文字幕三级久久日本| 丝袜脚勾引网站| 国产在线免费精品| 国产亚洲欧美精品永久| 蜜桃国产av成人99| 亚洲av日韩在线播放| 亚洲,一卡二卡三卡| av在线老鸭窝| 男女边摸边吃奶| 亚洲色图综合在线观看| 一级爰片在线观看| 男女边摸边吃奶| 久久青草综合色| av线在线观看网站| 欧美日韩综合久久久久久| 可以免费在线观看a视频的电影网站 | 亚洲av福利一区| 热re99久久精品国产66热6| 亚洲第一区二区三区不卡| 日韩成人av中文字幕在线观看| 亚洲国产精品国产精品| 亚洲一码二码三码区别大吗| 亚洲欧洲日产国产| 激情视频va一区二区三区| 久久精品人人爽人人爽视色| 国产一区亚洲一区在线观看| 老司机影院成人| 亚洲国产精品999| 免费播放大片免费观看视频在线观看| av在线播放精品| 精品国产一区二区三区四区第35| 桃花免费在线播放| tube8黄色片| 国产成人免费观看mmmm| 97在线人人人人妻| 中国三级夫妇交换| 精品亚洲成a人片在线观看| 欧美bdsm另类| 亚洲欧美日韩另类电影网站| 叶爱在线成人免费视频播放| 黄色毛片三级朝国网站| 成年人免费黄色播放视频| 日韩中文字幕视频在线看片| 精品99又大又爽又粗少妇毛片| 五月伊人婷婷丁香| 另类精品久久| 欧美日韩视频精品一区| 在线亚洲精品国产二区图片欧美| av免费观看日本| 水蜜桃什么品种好| 一级,二级,三级黄色视频| 人妻少妇偷人精品九色| 亚洲激情五月婷婷啪啪| 亚洲精品成人av观看孕妇| 日韩av不卡免费在线播放| 美女国产高潮福利片在线看| 欧美精品一区二区大全| 在线免费观看不下载黄p国产| 大话2 男鬼变身卡| 不卡av一区二区三区| 香蕉丝袜av| 午夜老司机福利剧场| 最近中文字幕2019免费版| 日韩欧美精品免费久久| 老司机影院毛片| 亚洲,一卡二卡三卡| 九草在线视频观看| 久久韩国三级中文字幕| 叶爱在线成人免费视频播放| 热99国产精品久久久久久7| 免费观看av网站的网址| 亚洲av免费高清在线观看| 九草在线视频观看| av国产久精品久网站免费入址| 国产高清不卡午夜福利| 国产在视频线精品| 亚洲一码二码三码区别大吗| 国产精品免费视频内射| 视频在线观看一区二区三区| 建设人人有责人人尽责人人享有的| 在线观看免费日韩欧美大片| 青春草国产在线视频| 黄片无遮挡物在线观看| 狠狠精品人妻久久久久久综合| 欧美日韩精品成人综合77777| 999精品在线视频| 婷婷色综合大香蕉| 久久精品久久久久久久性| 成人黄色视频免费在线看| 香蕉精品网在线| 日韩伦理黄色片| 国产欧美日韩一区二区三区在线| 久久精品国产自在天天线| 韩国精品一区二区三区| 国产成人免费观看mmmm| 美女中出高潮动态图| 国产精品 欧美亚洲| 蜜桃在线观看..| 国产片内射在线| 人妻一区二区av| 久久青草综合色| 女性生殖器流出的白浆| 最近手机中文字幕大全| 欧美精品人与动牲交sv欧美| 韩国高清视频一区二区三区| 国产精品二区激情视频| √禁漫天堂资源中文www| 国产成人午夜福利电影在线观看| xxxhd国产人妻xxx| 日本黄色日本黄色录像| 亚洲国产最新在线播放| 日本-黄色视频高清免费观看| 亚洲国产色片| 麻豆av在线久日| 免费黄色在线免费观看| 欧美精品一区二区免费开放| 亚洲美女搞黄在线观看| 中国三级夫妇交换| 久久久精品免费免费高清| 在线观看国产h片| 深夜精品福利| 卡戴珊不雅视频在线播放| 国产精品av久久久久免费| www日本在线高清视频| a级毛片在线看网站| 99久久精品国产国产毛片| 亚洲国产精品一区二区三区在线| 不卡av一区二区三区| 美女大奶头黄色视频| 三级国产精品片| 日本av免费视频播放| 国产爽快片一区二区三区| 一级毛片 在线播放| 天天影视国产精品| 中国国产av一级| 成年人免费黄色播放视频| 99九九在线精品视频| 免费女性裸体啪啪无遮挡网站| 亚洲欧美成人综合另类久久久| 成年av动漫网址| 搡女人真爽免费视频火全软件| 国产日韩欧美亚洲二区| 高清av免费在线| 国产又色又爽无遮挡免| 黄色毛片三级朝国网站| 免费观看av网站的网址| 少妇熟女欧美另类| 视频区图区小说| 欧美精品高潮呻吟av久久| 免费av中文字幕在线| 91精品三级在线观看| 精品卡一卡二卡四卡免费| 在线看a的网站| 国产日韩欧美视频二区| 人人妻人人澡人人看| 国产精品久久久av美女十八| √禁漫天堂资源中文www| 国产精品二区激情视频| 午夜老司机福利剧场| 成人毛片a级毛片在线播放| 国产黄色视频一区二区在线观看| 国产熟女午夜一区二区三区| 丝袜喷水一区| 亚洲av电影在线观看一区二区三区| 黄色一级大片看看| 又黄又粗又硬又大视频| 男女边摸边吃奶| 晚上一个人看的免费电影| av免费在线看不卡| 18+在线观看网站| 中文字幕av电影在线播放| 99国产综合亚洲精品| 97人妻天天添夜夜摸| 哪个播放器可以免费观看大片| av在线app专区| 亚洲综合色网址| 国产不卡av网站在线观看| 日本爱情动作片www.在线观看| 久久女婷五月综合色啪小说| 国产 精品1| av免费观看日本| 免费黄频网站在线观看国产| 久久久精品94久久精品| 久久久久国产网址| 在现免费观看毛片| 亚洲色图 男人天堂 中文字幕| 极品少妇高潮喷水抽搐| 亚洲精品国产色婷婷电影| 精品酒店卫生间| www日本在线高清视频| 久久99蜜桃精品久久| 韩国精品一区二区三区| www.精华液| 国产欧美亚洲国产| 中文乱码字字幕精品一区二区三区| 91国产中文字幕| 丝瓜视频免费看黄片| 中文精品一卡2卡3卡4更新| 亚洲精品一二三| 国产成人免费无遮挡视频| 天堂中文最新版在线下载| 丰满饥渴人妻一区二区三| 99热网站在线观看| 蜜桃在线观看..| 精品卡一卡二卡四卡免费| 母亲3免费完整高清在线观看 | 男女边摸边吃奶| 成年av动漫网址| 免费久久久久久久精品成人欧美视频| 黄色配什么色好看| 国产精品久久久久久精品电影小说| av免费观看日本| 国产免费又黄又爽又色| 久久青草综合色| 亚洲av日韩在线播放| 亚洲精品中文字幕在线视频| 国产精品人妻久久久影院| 黄片小视频在线播放| 精品国产国语对白av| 电影成人av| 精品99又大又爽又粗少妇毛片| 2018国产大陆天天弄谢| 大香蕉久久网| 水蜜桃什么品种好| 亚洲国产av新网站| 亚洲视频免费观看视频| 天天躁日日躁夜夜躁夜夜| 麻豆av在线久日| 国产在线一区二区三区精| 伦理电影免费视频| 国产精品一区二区在线观看99| 久久久久精品性色| 少妇被粗大猛烈的视频| 亚洲国产看品久久| 一二三四中文在线观看免费高清| 亚洲精品aⅴ在线观看| 中文字幕制服av| 欧美中文综合在线视频| 久久久久久久久久久免费av| 最近的中文字幕免费完整| 精品一区二区三区四区五区乱码 | 麻豆乱淫一区二区| 美女国产视频在线观看| 亚洲国产精品999| 中文字幕精品免费在线观看视频| 国产一区二区三区av在线| 国产免费又黄又爽又色| 啦啦啦视频在线资源免费观看| 老鸭窝网址在线观看| 黑人巨大精品欧美一区二区蜜桃| 国产高清国产精品国产三级| 亚洲成人av在线免费| 免费黄色在线免费观看| 亚洲av男天堂| 亚洲久久久国产精品| 中文乱码字字幕精品一区二区三区| 午夜av观看不卡| 巨乳人妻的诱惑在线观看| av网站在线播放免费| 五月开心婷婷网| 九色亚洲精品在线播放| 久久精品国产自在天天线| 80岁老熟妇乱子伦牲交| 性色avwww在线观看| 欧美亚洲 丝袜 人妻 在线| 新久久久久国产一级毛片| 亚洲欧洲精品一区二区精品久久久 | 欧美变态另类bdsm刘玥| 一级毛片我不卡| 午夜av观看不卡| 久久婷婷青草| 国产精品久久久久久av不卡| 1024视频免费在线观看| 久久久a久久爽久久v久久| 久久青草综合色| a级毛片在线看网站| 精品亚洲成a人片在线观看| 黄色怎么调成土黄色| 日本-黄色视频高清免费观看| 久久久精品区二区三区| 久久久国产一区二区| 波野结衣二区三区在线| 精品国产超薄肉色丝袜足j| av网站在线播放免费| 亚洲成色77777| 久久人人爽人人片av| 国产国语露脸激情在线看| 国产精品国产av在线观看| 欧美人与性动交α欧美软件| 26uuu在线亚洲综合色| 1024香蕉在线观看| 高清欧美精品videossex| 欧美bdsm另类| 宅男免费午夜| 中文乱码字字幕精品一区二区三区| 久久这里只有精品19| 欧美人与善性xxx| 国产成人aa在线观看| 日本vs欧美在线观看视频| 精品人妻熟女毛片av久久网站| 男人爽女人下面视频在线观看| 搡女人真爽免费视频火全软件| 亚洲男人天堂网一区| 97人妻天天添夜夜摸|