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

    入太湖污染物通量削減目標(biāo)下水工程調(diào)控優(yōu)化探究

    2024-06-12 11:56:55秦灝柳子豪姚杰夫唐仁吳小靖
    人民長江 2024年5期
    關(guān)鍵詞:湖西孟河入湖

    秦灝 柳子豪 姚杰夫 唐仁 吳小靖

    摘要:入湖污染物通量削減,一直是太湖水環(huán)境治理工作的重點(diǎn)。新孟河作為現(xiàn)狀唯一一條可以直接連接長江與太湖流域上游區(qū)域腹部的流域性重要河道,前期工作發(fā)現(xiàn),在新孟河工程實(shí)施后,隨著引長江清水進(jìn)入太湖,湖西區(qū)主要入湖河道的水質(zhì)得到明顯改善,但總體入太湖的污染物通量也進(jìn)一步增加。在新孟河工程即將常態(tài)化運(yùn)行的背景下,針對工程運(yùn)行后湖西區(qū)入湖污染物通量仍可能進(jìn)一步增加的難點(diǎn)問題,通過靈活調(diào)整工程調(diào)度,合理部署區(qū)域產(chǎn)水出路,計(jì)劃增加沿江口門排水機(jī)會,利用新孟河向長江排水,以拉動腹部地區(qū)水質(zhì)差的本地產(chǎn)水由東入太湖轉(zhuǎn)而北排長江,從而有效削減入太湖污染物通量。結(jié)果表明:通過優(yōu)化新孟河工程現(xiàn)狀調(diào)度,典型平水年(p=50%)時(shí)可以通過新孟河排江水量新增4.68億m3,湖西區(qū)入湖氨氮通量較現(xiàn)狀調(diào)度方案削減184 t/a、總磷入湖通量削減35.5 t/a,說明優(yōu)化后的水工程調(diào)度方案是利用水利工程調(diào)度實(shí)現(xiàn)入湖污染物通量削減的可行性措施。相關(guān)研究成果可為太湖水環(huán)境治理工作提供新的思路。

    關(guān) 鍵 詞:污染物通量;水利工程調(diào)度;新孟河;太湖

    中圖法分類號:TV213.4

    文獻(xiàn)標(biāo)志碼:A

    DOI:10.16232/j.cnki.1001-4179.2024.05.010

    0 引 言

    入湖污染物通量居高不下,一直是太湖水環(huán)境治理的重點(diǎn)和難點(diǎn)[1-2。以總磷為例,2010年以來入太湖磷通量年均2 224 t,出湖磷通量年均638 t[3;而太湖總磷容量僅為(420±100)t,磷滯留量高達(dá)1 265 t/a[4。污染物是隨水流匯入湖泊5-6,因此通過限制入湖水量來控制污染物輸入,是進(jìn)行湖泊水環(huán)境治理的有效辦法7-8。在太湖水環(huán)境治理工作中,通過嚴(yán)格管控太湖北部武澄錫虞區(qū)沿湖口門,限制區(qū)域內(nèi)水質(zhì)差的本地產(chǎn)水進(jìn)入太湖,該分區(qū)入湖污染物通量得到極大削減,太湖北部的梅梁湖、貢湖等湖區(qū)水環(huán)境也相應(yīng)改善9。

    長期以來,湖西區(qū)一直是入太湖污染物通量的主要輸入?yún)^(qū)域。以總磷為例,2008年以來,湖西區(qū)年均入湖總磷通量達(dá)1 730 t/a,占入太湖總磷通量的78.1%。在《太湖流域水環(huán)境綜合治理總體方案(2021~2035年)》(以下簡稱《太湖水環(huán)境治理總體方案》)中,將湖西區(qū)所在的江蘇省上游地區(qū)作為入湖污染物防控的重點(diǎn)區(qū)域[10。但由于湖西區(qū)主要入湖河道均未實(shí)施工程控制,因此通過控制湖西區(qū)入湖水量來削減入湖污染物通量就變得十分困難。

    新孟河延伸拓浚工程是《太湖水環(huán)境治理總體方案》中部署的“提高水環(huán)境容量(納污能力)引排工程”重點(diǎn)實(shí)施項(xiàng)目,工程開拓了原新孟河過京杭運(yùn)河,延伸至洮、滆湖地區(qū)[11,是湖西區(qū)現(xiàn)狀唯一一條可以直接連接長江與太湖流域上游區(qū)域腹部的流域性重要河道。本文在新孟河工程即將常態(tài)化運(yùn)行背景下,針對湖西區(qū)入湖污染物通量居高不下的流域水環(huán)境治理難點(diǎn)問題,靈活調(diào)整新孟河工程調(diào)度,通過增加新孟河沿江口門排水機(jī)會,拉動腹部洮、滆地區(qū)水質(zhì)差的本地產(chǎn)水由東入太湖轉(zhuǎn)而北排長江,從而有效削減湖西區(qū)入太湖污染物通量。相關(guān)研究成果可以為太湖水環(huán)境治理工作提供新的思路,對于進(jìn)一步發(fā)揮水利工程綜合效益、挖掘區(qū)域乃至流域水環(huán)境治理潛力具有重要的理論研究意義。

    1 研究區(qū)概況

    湖西區(qū)位于太湖流域西北部,屬于太湖上游地區(qū),總面積7 671 km2。區(qū)域內(nèi)水系發(fā)達(dá),根據(jù)地形特點(diǎn)分為北部運(yùn)河水系、中部洮滆水系和南部南河水系。北部運(yùn)河水系分布在運(yùn)河以北地區(qū),向北直接匯入長江;中部洮滆水系和南部南河水系分布在運(yùn)河以南地區(qū),向東匯入太湖。其中,中部洮滆水系位于湖西腹部地區(qū),主要承接西部茅山及金壇、丹陽一帶丘陵、高地來水,在匯入洮湖后經(jīng)由北干河、湟里河、夏溪河等橫向連接河道進(jìn)入滆湖,最后匯入太湖(圖1)。沿線區(qū)域工、農(nóng)業(yè)發(fā)達(dá),現(xiàn)狀河道水質(zhì)相對較差。

    新孟河北起長江,向南過京杭運(yùn)河后,延伸至北干河,全長116.47 km。新孟河連通長江處建有界牌水利樞紐,配置80 m寬閘門和設(shè)計(jì)流量為300 m3/s的雙向引排泵站,具備動力引排水能力。新孟河延伸拓浚工程建設(shè)河段除拓浚延伸新孟河之外,同時(shí)包括拓浚北干河,以連接長蕩湖和滆湖;拓浚太滆運(yùn)河、漕橋河,暢通工程引水入太湖的通道。2023年7月,工程通過投入使用驗(yàn)收。

    新孟河延伸拓浚工程具有改善水環(huán)境、提高防洪排澇和水資源配置能力等綜合功能和效益。其中,工程引水任務(wù)為平水年(2000年型)增加湖西區(qū)引江入湖水量25.2億m3;排水規(guī)劃任務(wù)為遇100 a一遇1991年型洪水,造峰期北排長江洪水7.5億~7.9億m3。根據(jù)工程初步設(shè)計(jì),當(dāng)太湖水位位于防洪調(diào)度區(qū),滆湖水位高于4.2 m時(shí),界牌水利樞紐開閘排水;當(dāng)滆湖水位超過4.6 m時(shí),界牌水利樞紐開泵排水。當(dāng)太湖水位處于適時(shí)調(diào)度區(qū),滆湖水位高于4.2 m時(shí),界牌水利樞紐開閘排水;而當(dāng)滆湖水位低于3.7 m時(shí),界牌水利樞紐開閘引水。目前,工程尚未進(jìn)入常態(tài)化運(yùn)行階段,2022年10~12月,江蘇省水利廳組織實(shí)施“新孟河抗旱調(diào)水試運(yùn)行”試驗(yàn),調(diào)水試驗(yàn)期間,通過引長江清水有效增加了湖西區(qū)特別是洮、滆兩湖水資源量,緩解了地區(qū)旱情,太湖水資源量也得到補(bǔ)充。

    在前期工作中,對新孟河及其一級支流沿線片區(qū)污染源情況進(jìn)行搜集整理和調(diào)查,重點(diǎn)關(guān)注新孟河沿線流經(jīng)的11個(gè)鄉(xiāng)鎮(zhèn),涉及常州市新北區(qū)、武進(jìn)區(qū)和金壇市[12-13,作為本文構(gòu)建區(qū)域水質(zhì)模型的重要支撐資料。其中,點(diǎn)源污染上,包括直排工業(yè)、污水廠等;面源污染上,包括生活直排、農(nóng)業(yè)種植和城市建成區(qū)地表徑流等。結(jié)果表明,新孟河及其一級支流沿線COD、氨氮、總氮、總磷污染物入河量分別為3 043.6,285.1,674.5,41.4 t,其中面源污染占比要大于點(diǎn)源污染。在不設(shè)控的運(yùn)河以南地區(qū),主要污染物入河通量明顯高于運(yùn)河以北地區(qū)(表1)。

    2 水量水質(zhì)模型構(gòu)建

    為定量化分析新孟河工程實(shí)施前、后湖西區(qū)入湖污染物通量及主要入湖河道水質(zhì)變化情況,選用MIKE 11模型進(jìn)行湖西區(qū)水量水質(zhì)模型構(gòu)建工作。MIKE 11模型在水環(huán)境治理工作中應(yīng)用廣泛,是適用于太湖流域平原河網(wǎng)地區(qū)一維非恒定流穩(wěn)定計(jì)算的有效手段[14-15。本次研究應(yīng)用降雨徑流模型(NAM)、水動力模塊(HD)與水質(zhì)模塊(AD),結(jié)合長期實(shí)測資料,對湖西區(qū)入太湖水量、污染物通量等參數(shù)進(jìn)行模擬計(jì)算分析。

    2.1 典型年選擇

    在《太湖流域水資源綜合規(guī)劃》[16等相關(guān)規(guī)劃中,明確在遭遇流域平水年(2000年型)時(shí),通過新孟河工程實(shí)施增加湖西區(qū)入太湖水量25.2億m3的引水任務(wù)要求。在本文中,對應(yīng)選擇2000年型作為典型平水年型開展計(jì)算分析。

    在2000年,太湖流域全年面雨量為1 074.2 mm,保證率為52%;湖西區(qū)全年降雨保證率為55%。流域、湖西區(qū)保證率基本一致,都接近于50%,可以作為P=50%的典型平水年雨型進(jìn)行模型計(jì)算。

    2.2 水系概化與模型邊界選擇

    2.2.1 水系概化

    根據(jù)湖西區(qū)水系分布、匯水范圍及河道特征,對區(qū)域中骨干河道進(jìn)行河網(wǎng)概化(圖2)。

    2.2.2 模型邊界選擇

    (1)產(chǎn)匯流計(jì)算。根據(jù)湖西區(qū)32個(gè)雨量站點(diǎn)分布劃分泰森多邊形,劃分出的每個(gè)多邊形作為流域平均降雨量計(jì)算單元。采用氣象站點(diǎn)典型平水年(2000年)降雨蒸發(fā)數(shù)據(jù)進(jìn)行產(chǎn)匯流計(jì)算,計(jì)算結(jié)果作為河網(wǎng)中的側(cè)向入流。

    (2)水動力邊界。模型共設(shè)置23處水動力外部邊界。其中,沿長江潮位邊界選擇連江骨干河道典型平水年(2000年)實(shí)測潮位資料;主要入太湖河道、與武澄錫虞區(qū)內(nèi)河連接河道選擇典型平水年(2000年)實(shí)測水位或流量資料;湖西區(qū)沿江主要水利工程均設(shè)置閘控運(yùn)行條件。

    (3)水質(zhì)邊界。水質(zhì)邊界采用2020年水動力外部邊界對應(yīng)點(diǎn)位的國、省考斷面月際尺度實(shí)測水質(zhì)(氨氮、總磷)資料。

    (4)區(qū)域污染負(fù)荷輸入。文中所使用的湖西區(qū)污染負(fù)荷數(shù)據(jù),與太湖上游鎮(zhèn)江、常州、無錫三市“十四五”太湖綜合治理規(guī)劃中統(tǒng)計(jì)數(shù)據(jù)一致。其中,點(diǎn)源包括直排工業(yè)源、城市生活源與農(nóng)村生活源,面源包括未接管生活源、城市面源、畜禽養(yǎng)殖源、水產(chǎn)養(yǎng)殖源與農(nóng)業(yè)種植源,各污染源占比見表2。同時(shí),將湖西區(qū)劃分為28個(gè)控制單元,在模型中點(diǎn)面源污染負(fù)荷均按照各匯水范圍內(nèi)河道側(cè)向入流流量均分,并以側(cè)向入流污染物濃度的形式進(jìn)入河道。

    2.3 模型率定與驗(yàn)證

    2.3.1 水量模型率定與驗(yàn)證

    (1)水量模型參數(shù)率定。選擇2000年1~6月金壇、溧陽站、王母觀日水位數(shù)據(jù),進(jìn)行HD模型參數(shù)率定。主要率定參數(shù)為河床糙率,糙率率定值為0.033。NAM模型參數(shù)依據(jù)HD模型模擬結(jié)果以及區(qū)域經(jīng)驗(yàn)值選取。

    (2)水量模型驗(yàn)證誤差。選擇2000年7~12月金壇、王母觀、溧陽站日水位數(shù)據(jù)進(jìn)行模型驗(yàn)證。結(jié)果表明,各站水位模擬與實(shí)測值之間平均相對誤差均在±5%以內(nèi)(圖3),說明模擬擬合在可接受區(qū)間范圍,本次建立的區(qū)域水動力模型具有較好的計(jì)算精度。

    2.3.2 水質(zhì)率定與驗(yàn)證

    選擇2020年1~6月東潘橋、辛豐鎮(zhèn)、靜堂大橋斷面水質(zhì)數(shù)據(jù)進(jìn)行水質(zhì)參數(shù)率定。污染物衰減系數(shù)根據(jù)文獻(xiàn)[17]取固定經(jīng)驗(yàn)值,氨氮衰減系數(shù)率定值為0.13/d,總磷衰減系數(shù)率定值為0.008/d。利用2020年7~12月對應(yīng)斷面水質(zhì)數(shù)據(jù)進(jìn)行模型驗(yàn)證(圖4~5)。結(jié)果表明,東潘橋斷面氨氮平均絕對相對誤差為-14.6%,總磷平均絕對相對誤差為14.7%;辛豐鎮(zhèn)斷面氨氮平均絕對相對誤差為21.3%,總磷平均絕對相對誤差為11.1%;靜堂大橋斷面氨氮平均絕對相對誤差為12.1%,總磷平均絕對相對誤差為10.6%。說明構(gòu)建的模型具有較好的計(jì)算精度,可以滿足本次工作計(jì)算要求。

    3 現(xiàn)狀入湖污染物通量特征分析

    3.1 實(shí)測資料分析

    實(shí)測資料數(shù)據(jù)來源于江蘇省水文水資源勘測局,資料年限為1986~2020年。

    3.1.1 湖西區(qū)入湖水量

    湖西區(qū)多年平均入湖水量為51.4億m3,占環(huán)太湖入湖總水量的56.4%。以太湖水環(huán)境綜合治理啟動的2007年為時(shí)間節(jié)點(diǎn),在2007~2020年期間,湖西區(qū)年均入湖水量為71.8億m3,較1986~2006年期間增加34.0億m3,增幅達(dá)到90.2%(圖6)。湖西區(qū)在入湖水量中占比也相應(yīng)增加,在2007~2020年期間,湖西區(qū)占總?cè)牒康?3.0%,較1986~2006年期間增加16.4%,是環(huán)太湖入湖水量增加的主要來源。

    3.1.2 湖西區(qū)入湖水質(zhì)

    自2008年太湖啟動水環(huán)境治理工作以來,湖西區(qū)主要入湖河道水質(zhì)改善明顯。其中,新孟河工程拓浚的太滆運(yùn)河、漕橋河,在2008~2020年期間主要水質(zhì)指標(biāo)(高錳酸鹽指數(shù)、氨氮、總氮、總磷)均出現(xiàn)降低趨勢。以總磷為例,太滆運(yùn)河總磷濃度從0.339 mg/L(2008年)降低至0.290 mg/L(2020年),漕橋河總磷濃度從0.340 mg/L(2008年)降低至0.208 mg/L(2020年)(圖7)。其中,在2012年以后,主要水質(zhì)指標(biāo)的變化趨于穩(wěn)定。

    3.1.3 湖西區(qū)入湖污染物通量

    從湖西區(qū)入湖污染物通量變化來看,與水質(zhì)相對應(yīng),入湖污染物通量也表現(xiàn)出整體減少的趨勢。以總磷為例,從2008年的1 779 t/a,降低至2020年的1 537 t/a(圖8)。

    3.2 工程實(shí)施前、后入湖污染物通量變化分析

    參照新孟河工程初步設(shè)計(jì)中工程調(diào)度原則,對典型平水年(2000年型)時(shí)新孟河工程實(shí)施前、后湖西區(qū)入湖水量、主要入湖河道水質(zhì)及入湖污染物通量等參數(shù)進(jìn)行模擬計(jì)算,重點(diǎn)關(guān)注氨氮、總磷水質(zhì)指標(biāo)。一方面分析工程實(shí)施對流域水環(huán)境可能產(chǎn)生的影響,另一方面也為下一步進(jìn)行工程調(diào)度的優(yōu)化調(diào)整明確基礎(chǔ)背景。

    3.2.1 新孟河工程實(shí)施前

    (1)水量。在典型平水年時(shí),新孟河工程實(shí)施前,現(xiàn)狀工程調(diào)度下,湖西區(qū)入太湖水量為55.26億m3。其中洮滆水系入湖水量為26.86億m3,其主要入湖河道為太滆南運(yùn)河、太滆運(yùn)河、漕橋河,入湖水量分別為9.95億m3和5.56億m3(太滆運(yùn)河與漕橋河匯流)。該時(shí)期,通過新孟河引江水量為2.87億m3,通過新孟河排江水量僅為0.18億m3

    (2)水質(zhì)。在典型平水年時(shí),新孟河工程實(shí)施前,從主要入湖河道的水質(zhì)年均值變化來看,入湖河道平均氨氮濃度為0.427 mg/L;平均總磷濃度為0.140 mg/L。具體入湖河道來看,在氨氮指標(biāo)上,太滆運(yùn)河、漕橋河匯合后的分水橋斷面氨氮濃度為0.550 mg/L,太滆南運(yùn)河為0.312 mg/L,東湛瀆港為0.611 mg/L,南河為0.414 mg/L。在總磷指標(biāo)上,分水橋斷面總磷濃度為0.173 mg/L,太滆南運(yùn)河為0.230 mg/L,東湛瀆港為0.114 mg/L,南河為0.136 mg/L。湖西區(qū)主要入湖河道水質(zhì)指標(biāo)情況見圖9。

    (3)入湖污染物通量。在典型平水年時(shí),新孟河工程實(shí)施前,湖西區(qū)氨氮入湖通量為2 135 t/a,總磷入湖通量為840 t/a,計(jì)算結(jié)果與相關(guān)研究中對該年型下湖西區(qū)入湖總磷通量的估算值(893 t/a)基本一致[18。具體從入湖河道來看,在氨氮指標(biāo)上,太滆運(yùn)河、漕橋河匯合后的分水橋斷面氨氮入湖通量為285 t/a,太滆南運(yùn)河為297 t/a,東湛瀆港為166 t/a,南河為1 083 t/a。在總磷指標(biāo)上,分水橋斷面總磷入湖通量為98 t/a,太滆南運(yùn)河為230 t/a,東湛瀆港為32 t/a,南河為378 t/a。湖西區(qū)主要入湖河道污染物入湖通量情況見圖10。

    3.2.2 新孟河工程實(shí)施后

    (1)水量。在典型平水年時(shí),新孟河工程實(shí)施后,湖西區(qū)入太湖水量為88.11億m3,較工程實(shí)施前湖西區(qū)新增入湖水量32.85億m3。其中,洮滆水系入湖水量為56.86億m3,主要入湖河道為太滆運(yùn)河、漕橋河、太滆南運(yùn)河。新孟河工程拓浚的河線太滆運(yùn)河、漕橋河匯流后入太湖水量為35.32億m3,較工程實(shí)施前增加29.76億m3,是湖西區(qū)新增入湖水量的主要通道。該時(shí)期,通過新孟河引江水量為45.28億m3,通過新孟河排江水量僅為0.02億m3,工程運(yùn)行以引水為主。

    (2)水質(zhì)。在典型平水年時(shí),新孟河工程實(shí)施后,主要入湖河道的水質(zhì)均出現(xiàn)明顯改善,從年均值變化來看,入湖河道平均氨氮濃度為0.326 mg/L;平均總磷濃度為0.101 mg/L。具體從入湖河道來看,作為新孟河引水入太湖的重要通道,太滆運(yùn)河、漕橋河水質(zhì)改善最為明顯。在太滆運(yùn)河、漕橋河匯合后的分水橋斷面上,氨氮、總磷濃度分別從工程實(shí)施前的0.550 mg/L和0.173 mg/L,降低至工程實(shí)施后的0.351 mg/L和0.125 mg/L。湖西區(qū)主要入湖河道水質(zhì)指標(biāo)情況見圖11。

    (3)入湖污染物通量。相應(yīng)于水質(zhì)的明顯改善,在典型平水年時(shí),新孟河工程實(shí)施后,隨著入湖水量的大幅增加,主要入湖河道的污染物入湖通量也出現(xiàn)增加。其中,湖西區(qū)氨氮入湖通量增加至2 343 t/a,總磷入湖通量增加至1 042 t/a。具體從入湖河道來看,由于湖西區(qū)新增入湖水量主要是通過太滆運(yùn)河、漕橋河進(jìn)入太湖,這兩條河道的污染物入湖通量增幅明顯。在太滆運(yùn)河、漕橋河匯合后的分水橋斷面上,氨氮、總磷入湖通量分別從工程實(shí)施前的285 t/a和98 t/a,增加至工程實(shí)施后的1 025 t/a和427 t/a。湖西區(qū)主要入湖河道入太湖污染物通量情況見圖12。

    4 以入湖污染物通量削減為目標(biāo)的工程調(diào)度優(yōu)化

    4.1 調(diào)度優(yōu)化思路及方案

    4.1.1 優(yōu)化思路

    入湖污染物通量的變化,與入湖河道水質(zhì)及入湖水量有直接聯(lián)系。實(shí)測資料表明,在湖西區(qū)入湖水量整體增加背景下,入湖河道水質(zhì)也趨于穩(wěn)定,因此從進(jìn)一步改善水質(zhì)的角度來實(shí)現(xiàn)入湖污染物通量較大幅度的削減就十分困難。從新孟河工程的實(shí)施效果來看,工程實(shí)施后,湖西區(qū)入湖水量顯著增加,主要入湖河道水質(zhì)也得到明顯改善,但同時(shí)湖西區(qū)入湖污染物通量也出現(xiàn)了進(jìn)一步增加的情況,這可能會對太湖水環(huán)境造成負(fù)面影響。結(jié)合新孟河沿江口門的引排情況進(jìn)行分析,在平水年,界牌水利樞紐以引水為主,排江水量十分有限。這是由于在工程現(xiàn)狀調(diào)度中,界牌水利樞紐進(jìn)行排水調(diào)度的啟用水位相對較高,難以通過增加排水來減少該時(shí)期湖西區(qū)入太湖水量。

    在本次研究中,針對界牌水利樞紐現(xiàn)狀調(diào)度原則,考慮增加新孟河沿江口門的排水機(jī)會,利用新孟河拉動湖西區(qū)腹部水質(zhì)差的本地產(chǎn)水由東入太湖轉(zhuǎn)而北排長江,從而達(dá)到有效削減入湖污染物通量的目的。

    4.1.2 方案設(shè)計(jì)

    在本次方案設(shè)計(jì)中,以滿足新孟河工程引水任務(wù)為前提,同時(shí)考慮太湖、滆湖等湖泊最低水位保障要求。在工程運(yùn)行調(diào)度水位的確定上,一方面考慮合理增加新孟河沿江口門排水機(jī)會;另一方面考慮在排水調(diào)度時(shí)盡可能地使用水閘進(jìn)行排水,以節(jié)省工程運(yùn)行成本。經(jīng)模型模擬比選后,最終確定坊前站(滆湖)水位達(dá)到3.5 m時(shí)作為界牌水利樞紐調(diào)整后的排水調(diào)度水位。具體調(diào)度方式見表3。

    4.2 效果分析

    (1)水量。在界牌水利樞紐調(diào)度方案調(diào)整后,典型平水年(2000年型)時(shí),通過新孟河引江水量為32.55億m3,湖西區(qū)入太湖水量為80.84億m3,相較工程實(shí)施前新增入太湖水量達(dá)25.58億m3,滿足流域水資源規(guī)劃中對工程引水的任務(wù)要求[16。該時(shí)期,太湖年內(nèi)平均水位3.27 m,最低水位3.01 m,滿足太湖最低生態(tài)水位要求(2.65 m)[19;滆湖年內(nèi)平均水位3.44 m,最低水位3.15 m,滿足滆湖最低生態(tài)水位要求(2.67 m)[20。

    調(diào)度方案調(diào)整后,在年內(nèi)1月份和6月份,界牌水利樞紐執(zhí)行排水調(diào)度(圖13),在排水運(yùn)行時(shí),新孟河運(yùn)南段東西兩側(cè)支流(鶴溪河、夏溪河、湟里河、北干河)水流方向均轉(zhuǎn)為向新孟河匯入,新孟河流向也轉(zhuǎn)為自南向北進(jìn)入長江,有效拉動了洮滆腹部地區(qū)的本地產(chǎn)水向北排入長江。年內(nèi)湖西區(qū)通過新孟河北排長江水量為4.79億m3,較工程實(shí)施前新增排江水量4.61億m3;同時(shí)相較于現(xiàn)狀調(diào)度方案,減少湖西區(qū)入太湖水量7.27億m3,降幅約8.3%。

    (2)水質(zhì)。在優(yōu)化調(diào)度方案實(shí)行下,湖西區(qū)入湖河道水質(zhì)明顯改善。典型平水年時(shí),新孟河工程實(shí)施后,在年均值變化上,入湖河道平均氨氮濃度為0.309 mg/L,平均總磷濃度為0.120 mg/L。相較于現(xiàn)狀調(diào)度方案,優(yōu)化方案下入湖河道氨氮濃度進(jìn)一步降低,總磷濃度則略微升高。湖西區(qū)主要入湖河道水質(zhì)指標(biāo)情況見圖14。

    (3)入湖污染物通量。典型平水年時(shí),相較于現(xiàn)狀調(diào)度方案,隨著腹部地區(qū)水質(zhì)差的本地產(chǎn)水通過新孟河向北增加排江,湖西區(qū)入湖污染物通量也明顯減少。在優(yōu)化調(diào)度方案實(shí)行下,湖西區(qū)氨氮入湖通量減少至2 159 t/a,較現(xiàn)狀方案削減184 t/a;總磷入湖通量減少至1 006 t/a,較現(xiàn)狀方案削減35.5 t/a。

    從太滆運(yùn)河、漕橋河入湖污染物通量變化來看,優(yōu)化調(diào)度方案下,在太滆運(yùn)河、漕橋河匯合后的分水橋斷面上,氨氮、總磷入湖通量分別從現(xiàn)狀調(diào)度實(shí)行下的1 025 t/a和427 t/a,削減至907 t/a和367 t/a。湖西區(qū)主要入湖河道入太湖污染物通量情況見圖15。

    5 結(jié)論與建議

    5.1 結(jié) 論

    (1)在新孟河工程實(shí)施后,典型平水年(2000年)時(shí),通過新孟河引水,湖西區(qū)入湖水量顯著增加,主要入湖河道水質(zhì)明顯改善,但同時(shí)湖西區(qū)入湖污染物通量也存在一定增幅。

    (2) 通過對界牌水利樞紐調(diào)度規(guī)則進(jìn)行優(yōu)化調(diào)整,降低工程排水啟用條件,利用新孟河拉動湖西區(qū)腹部水質(zhì)差的本地產(chǎn)水由東入太湖轉(zhuǎn)而北排長江,可以在典型平水年時(shí)增加工程北排長江水量4.61億m3,湖西區(qū)入湖氨氮通量較現(xiàn)狀調(diào)度方案削減184 t/a,總磷通量削減35.5 t/a,該方案是通過靈活調(diào)度水利工程實(shí)現(xiàn)入湖污染物通量削減的有效措施。

    5.2 建 議

    通過水利工程優(yōu)化調(diào)度來促使本地水質(zhì)較差的產(chǎn)水由入湖轉(zhuǎn)為排江,可以在一定程度上削減入湖污染負(fù)荷,進(jìn)而有利于改善湖泊水環(huán)境。但作為以“控源截污”為主線的太湖流域水環(huán)境綜合治理工作中,水利工程的優(yōu)化調(diào)控僅可以作為一種輔助手段。對于流域上游地區(qū)的水環(huán)境治理,建議持續(xù)推進(jìn)城鎮(zhèn)生活污水和雨洪排口的規(guī)范化整治,同時(shí)多模式開展農(nóng)村生活污水治理設(shè)施建設(shè),通過源頭充分減量、過程精準(zhǔn)控制、末端截蓄治并舉的系統(tǒng)治理思路[21-22,消除旱天非雨出流,控制雨天溢流污染,進(jìn)一步削減點(diǎn)面源污染,提高入湖河道水環(huán)境質(zhì)量,實(shí)現(xiàn)從源頭上有效削減入太湖污染負(fù)荷。

    參考文獻(xiàn):

    [1]陳凱麒,祁昌軍,陶潔,等.關(guān)于太湖水環(huán)境治理的思考[J].環(huán)境保護(hù),2016,44(18):20-23.

    [2]周雨,王殿常,余甜雪,等.典型河湖治理經(jīng)驗(yàn)及對長江流域生態(tài)保護(hù)的啟示[J].人民長江,2024,55(1):45-50,59.

    [3]水利部太湖流域管理局.太湖健康狀況報(bào)告[R].上海:水利部太湖流域管理局,2010-2020.

    [4]朱偉,胡思遠(yuǎn),馮甘雨,等.特大洪水對淺水湖泊磷的影響:以2016年太湖為例[J].湖泊科學(xué),2020,32(2):325-336.

    [5]LI L Q,SHAN B Q,YIN C Q.Stormwater runoff pollution loads from an urban catchment with rainy climate in China[J].Frontiers of Environmental Sciences & Engineering,2012,6(5):672-677.

    [6]LUO L C,QIN B Q,YANG L Y,et al.Total inputs of phosphorus and nitrogen by wet deposition into Lake Taihu,China[J].Hydrobiologia,2007,581(1):63-70.

    [7]秦伯強(qiáng).淺水湖泊湖沼學(xué)與太湖富營養(yǎng)化控制研究[J].湖泊科學(xué),2020,32(5):1229-1243.

    [8]程浩淼,王安安,張健,等.調(diào)水工程對過水性湖泊中污染物遷移轉(zhuǎn)化的影響[J].人民長江,2023,54(3):49-55,71.

    [9]申琳燕.太湖環(huán)境治理研究[D].南京:南京師范大學(xué),2013.

    [10]國家發(fā)展與改革委員會.太湖流域水環(huán)境綜合治理總體方案(2021-2035年)[R].北京:國家發(fā)展與改革委員會,2022.

    [11]國家發(fā)展與改革委員會.太湖流域水環(huán)境綜合治理總體方案[R].北京:國家發(fā)展與改革委員會,2013.

    [12]常州市新北區(qū)人民政府.常州市新北區(qū)“十四五”太湖流域綜合治理規(guī)劃[R].常州:常州市新北區(qū)人民政府,2021.

    [13]江蘇省生態(tài)環(huán)境廳.重點(diǎn)控制斷面及上溯主要入湖河流水環(huán)境綜合狀況分析及支流支浜分期分類管理研究[R].南京:江蘇省生態(tài)環(huán)境廳,2020.

    [14]王海森,馮嵐,李百煉,等.基于MIKE 11的河網(wǎng)區(qū)域小型閘壩調(diào)控的水質(zhì)響應(yīng)模擬研究[J].環(huán)境科學(xué)學(xué)報(bào),2022,42(8):314-324.

    [15]WANG Q G,WANG Y P,LU X C,et al.Impact assessments of water allocation on water environment of river network:method and application[J].Physics & Chemistry of the Earth Parts A/b/c,2017,103:101-106.

    [16]水利部太湖流域管理局.太湖流域水資源綜合規(guī)劃[R].上海:水利部太湖流域管理局,2010.

    [17]馮帥,李敘勇,鄧建才.太湖流域上游平原河網(wǎng)污染物綜合衰減系數(shù)的測定[J].環(huán)境科學(xué)學(xué)報(bào),2017,37(3):878-887.

    [18]陸昊,楊柳燕,楊明月,等.太湖流域上游降水量對入湖總氮和總磷的影響[J].水資源保護(hù),2022,38(4):174-181.

    [19]水利部太湖流域管理局.太湖生態(tài)水位保障實(shí)施方案[R].上海:水利部太湖流域管理局,2020.

    [20]江蘇省人民政府.滆湖保護(hù)規(guī)劃[R].南京:江蘇省人民政府,2022.

    [21]秦伯強(qiáng),吳海斌.長江流域湖泊富營養(yǎng)化發(fā)展趨勢與展望[J].人民長江,2023,54(10):18-23.

    [22]張事,馬方凱,徐望朋,等.淺水湖泊草-藻型穩(wěn)態(tài)轉(zhuǎn)換理論研究進(jìn)展[J].水利水電快報(bào),2022,43(12):121-127.

    (編輯:劉 媛)

    Optimization of water engineering regulation aiming to reduce pollutant flux into Taihu Lake

    QIN Hao1,LIU Zihao1,YAO Jiefu2,TANG Ren1,WU Xiaojing1

    (1.Jiangsu Taihu Water Conservancy Planning and Design Institute Co.,Ld.,Suzhou 215006,China;2.Earth Sciences and Engineering,Nanjing University,Nanjing,210023,China)

    Abstract:The reduction of pollutants flux into lake is the focus of water environment improvement in Taihu Lake.The Xinmeng River Project is the only important river channel that can directly connect the Changjiang River and the inner area of the upper reaches of Taihu Lake Basin.Previous work has found that since the implementation of the project,with the flow diverted from the Changjiang River into Taihu Lake,the water quality of the rivers into the Huxi hydraulic region has been significantly improved,but the overall pollutant flux into Taihu Lake has also further increased.Aiming at this problem after the normal operation of Xinmeng River project,we are about to adjust the projects′ dispatch flexibly and arrange water flowing route reasonably,in which we will increase the drainage outlets along the Changjiang River.In detail,we utilize the Xinmeng River to drain the lake water to the Changjiang River so as to promote the hinterland poor quality water to be drained to the Changjaing River,so the pollutant flux into the Taihu Lake will be effectively reduced.The results show that by optimizing the current scheduling of the Xinmeng River project,the amount of water discharged from the Xinmeng River can be increased by 468 million m3in the typical normal year (p=50%).Compared with the current scheduling scheme,the ammonia nitrogen flux into the Huxi area is reduced by 184 t/a,and the total phosphorus flux into the lake is reduced by 35.5 t/a.It shows that the optimized water project scheduling scheme is a feasible measure to reduce the pollutant flux into the lake by using the water project scheduling.The relevant research results can provide new ideas for the water environment treatment of Taihu Lake.

    Key words:pollution flux;water conservancy project regulation;Xinmeng River;Taihu Lake

    猜你喜歡
    湖西孟河入湖
    降雨事件下洱海北部主要入湖河流污染物特征分析*
    新時(shí)代湖西紅色文化傳承與發(fā)展創(chuàng)新研究*
    從學(xué)術(shù)傳承角度探討孟河費(fèi)氏流派對濕證的認(rèn)識和發(fā)展
    環(huán)洱海主要入湖河流水質(zhì)特征及入湖污染負(fù)荷估算
    人民長江(2022年1期)2022-04-01 02:44:35
    傳承孟河醫(yī)派 弘揚(yáng)國醫(yī)精髓
    基于中醫(yī)傳承輔助平臺探討孟河醫(yī)派治療水氣病的組方用藥規(guī)律研究
    龍湖·首開湖西星辰美學(xué)館
    孟河醫(yī)派三大家及傳人檔案追蹤
    入湖河口濕地恢復(fù)與重建規(guī)劃設(shè)計(jì)初探——以資興市興寧河入湖河口濕地為例
    滇池入湖河流磷負(fù)荷時(shí)空變化及形態(tài)組成貢獻(xiàn)*
    国产伦精品一区二区三区视频9 | 久久草成人影院| 97超级碰碰碰精品色视频在线观看| 久久久久国产一级毛片高清牌| 亚洲精品456在线播放app | 国产成+人综合+亚洲专区| av福利片在线观看| 老汉色∧v一级毛片| 中文字幕人成人乱码亚洲影| 国产av不卡久久| 性色av乱码一区二区三区2| 成人无遮挡网站| 久久久久国产一级毛片高清牌| 久久久精品大字幕| 91在线观看av| 国产极品精品免费视频能看的| 美女大奶头视频| 国产欧美日韩精品亚洲av| 日韩av在线大香蕉| 国产麻豆成人av免费视频| 国产精品永久免费网站| 国产伦精品一区二区三区四那| 一个人看的www免费观看视频| 欧美+亚洲+日韩+国产| 丰满的人妻完整版| 亚洲国产欧美人成| 性色avwww在线观看| 亚洲精品456在线播放app | 午夜福利在线在线| 97人妻精品一区二区三区麻豆| 国产主播在线观看一区二区| 九色成人免费人妻av| 久久中文字幕人妻熟女| 午夜福利在线观看吧| 天天躁日日操中文字幕| 国产美女午夜福利| 一本精品99久久精品77| 超碰成人久久| 好男人电影高清在线观看| 国产欧美日韩一区二区三| 偷拍熟女少妇极品色| 亚洲欧美日韩高清专用| 欧美国产日韩亚洲一区| 久久精品亚洲精品国产色婷小说| 日韩免费av在线播放| 国产精品久久久久久精品电影| 叶爱在线成人免费视频播放| 麻豆国产av国片精品| 国产精品久久电影中文字幕| 国产伦一二天堂av在线观看| 好男人电影高清在线观看| 免费电影在线观看免费观看| 两个人视频免费观看高清| 国产精品美女特级片免费视频播放器 | 欧美一级a爱片免费观看看| 曰老女人黄片| 色在线成人网| 天堂√8在线中文| 久久天躁狠狠躁夜夜2o2o| 久久这里只有精品19| 亚洲av成人精品一区久久| 激情在线观看视频在线高清| 午夜影院日韩av| 国产午夜精品久久久久久| 中国美女看黄片| 中出人妻视频一区二区| 欧美日本视频| 午夜福利高清视频| 亚洲精品色激情综合| 欧美黑人巨大hd| 欧美丝袜亚洲另类 | 精品人妻1区二区| 精品日产1卡2卡| 亚洲人成网站高清观看| 曰老女人黄片| 成人三级黄色视频| 久久久精品欧美日韩精品| 久久久国产成人免费| 日本三级黄在线观看| 九色国产91popny在线| 欧美大码av| 国产探花在线观看一区二区| 国产伦一二天堂av在线观看| 精品人妻1区二区| 久久精品综合一区二区三区| 亚洲色图 男人天堂 中文字幕| 日本在线视频免费播放| 欧美激情久久久久久爽电影| 女人高潮潮喷娇喘18禁视频| 国产午夜精品久久久久久| 午夜精品在线福利| 九色国产91popny在线| 国产精品99久久99久久久不卡| 亚洲成人中文字幕在线播放| 国产精品精品国产色婷婷| 精品熟女少妇八av免费久了| 久久久久久久久免费视频了| 国产精品1区2区在线观看.| 美女免费视频网站| 久久香蕉国产精品| 给我免费播放毛片高清在线观看| 中文在线观看免费www的网站| 精品福利观看| 国产不卡一卡二| 国产黄色小视频在线观看| 99久久精品国产亚洲精品| 一区二区三区国产精品乱码| 12—13女人毛片做爰片一| 亚洲精品色激情综合| 国产成人影院久久av| 国产精品久久电影中文字幕| 久久久久国产精品人妻aⅴ院| 午夜福利18| 久久中文字幕一级| 99热6这里只有精品| 美女被艹到高潮喷水动态| 亚洲精品在线美女| 亚洲,欧美精品.| 18美女黄网站色大片免费观看| 成熟少妇高潮喷水视频| 熟妇人妻久久中文字幕3abv| 午夜福利18| 免费看美女性在线毛片视频| 美女高潮的动态| 韩国av一区二区三区四区| 人人妻,人人澡人人爽秒播| 久久久久久国产a免费观看| 精品熟女少妇八av免费久了| 99热精品在线国产| 人妻夜夜爽99麻豆av| 深夜精品福利| 日本一二三区视频观看| 真实男女啪啪啪动态图| 精品不卡国产一区二区三区| 欧美午夜高清在线| 欧美成人一区二区免费高清观看 | 1024香蕉在线观看| 久久天躁狠狠躁夜夜2o2o| 丰满人妻熟妇乱又伦精品不卡| 黄色成人免费大全| 99久久精品热视频| 中国美女看黄片| 黄色丝袜av网址大全| 在线观看免费视频日本深夜| 国产精品久久久久久精品电影| 麻豆av在线久日| 中文字幕最新亚洲高清| 黑人巨大精品欧美一区二区mp4| 久久99热这里只有精品18| 中亚洲国语对白在线视频| 国内久久婷婷六月综合欲色啪| 国产精品一区二区精品视频观看| av欧美777| 天天躁日日操中文字幕| 在线播放国产精品三级| 国语自产精品视频在线第100页| 在线观看一区二区三区| 中亚洲国语对白在线视频| 国产av一区在线观看免费| 婷婷精品国产亚洲av| 国产主播在线观看一区二区| 久久亚洲真实| 国产又黄又爽又无遮挡在线| 九九在线视频观看精品| 亚洲,欧美精品.| 麻豆国产av国片精品| 美女被艹到高潮喷水动态| 午夜视频精品福利| 在线永久观看黄色视频| 久久婷婷人人爽人人干人人爱| 成年女人永久免费观看视频| 欧美乱妇无乱码| 国产亚洲欧美98| 亚洲五月天丁香| 日韩有码中文字幕| 两个人看的免费小视频| 成人永久免费在线观看视频| 欧美又色又爽又黄视频| 色综合站精品国产| 精品无人区乱码1区二区| 韩国av一区二区三区四区| 亚洲国产精品999在线| 99国产精品一区二区蜜桃av| 亚洲精品乱码久久久v下载方式 | 成人国产一区最新在线观看| 亚洲性夜色夜夜综合| 成年女人毛片免费观看观看9| 老司机深夜福利视频在线观看| 欧美日本亚洲视频在线播放| 岛国在线免费视频观看| www.精华液| 中文亚洲av片在线观看爽| 国产av不卡久久| 亚洲人成伊人成综合网2020| 亚洲自偷自拍图片 自拍| 午夜a级毛片| 国产三级在线视频| 精品国产三级普通话版| 欧美日韩乱码在线| 久久草成人影院| 国产精品一区二区精品视频观看| 97超视频在线观看视频| 99精品久久久久人妻精品| 真人做人爱边吃奶动态| 2021天堂中文幕一二区在线观| 搡老妇女老女人老熟妇| 欧美中文综合在线视频| 一区二区三区激情视频| 十八禁网站免费在线| 女警被强在线播放| 好看av亚洲va欧美ⅴa在| 午夜福利欧美成人| 丝袜人妻中文字幕| 亚洲 欧美一区二区三区| 大型黄色视频在线免费观看| 日韩欧美免费精品| 国产激情久久老熟女| 三级男女做爰猛烈吃奶摸视频| 亚洲av免费在线观看| 欧美日韩一级在线毛片| 久久精品国产综合久久久| 久久人人精品亚洲av| 两个人视频免费观看高清| 国产av一区在线观看免费| 身体一侧抽搐| 听说在线观看完整版免费高清| 成人亚洲精品av一区二区| 国产精品久久久久久久电影 | 51午夜福利影视在线观看| 蜜桃久久精品国产亚洲av| 全区人妻精品视频| 亚洲国产欧美人成| 99久久综合精品五月天人人| 久久香蕉精品热| 久久久久久久午夜电影| 狂野欧美激情性xxxx| 三级男女做爰猛烈吃奶摸视频| 免费高清视频大片| 97超级碰碰碰精品色视频在线观看| 精品福利观看| 哪里可以看免费的av片| 99久久精品国产亚洲精品| 欧美激情久久久久久爽电影| 国产精品免费一区二区三区在线| 一进一出抽搐动态| 亚洲自偷自拍图片 自拍| 国产精品久久久久久精品电影| 别揉我奶头~嗯~啊~动态视频| 欧美日韩国产亚洲二区| 他把我摸到了高潮在线观看| 亚洲成人久久爱视频| 日本熟妇午夜| 日韩欧美三级三区| 在线看三级毛片| 国产精品99久久99久久久不卡| 亚洲av成人av| 亚洲精品456在线播放app | 舔av片在线| 国产激情久久老熟女| 精品国产超薄肉色丝袜足j| 欧美一区二区国产精品久久精品| 亚洲第一欧美日韩一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 变态另类丝袜制服| 少妇的丰满在线观看| 亚洲熟女毛片儿| 露出奶头的视频| 日韩欧美在线二视频| 1000部很黄的大片| а√天堂www在线а√下载| 欧美大码av| 午夜福利高清视频| 日韩欧美免费精品| 真实男女啪啪啪动态图| 亚洲乱码一区二区免费版| 天天躁日日操中文字幕| 国产精品野战在线观看| 美女午夜性视频免费| 成年人黄色毛片网站| 在线国产一区二区在线| 国产精品av视频在线免费观看| 久久久色成人| 欧美绝顶高潮抽搐喷水| 宅男免费午夜| 亚洲一区高清亚洲精品| 全区人妻精品视频| 少妇裸体淫交视频免费看高清| 一本精品99久久精品77| 亚洲精品一区av在线观看| 这个男人来自地球电影免费观看| 波多野结衣巨乳人妻| 午夜a级毛片| bbb黄色大片| 欧美乱色亚洲激情| 99精品欧美一区二区三区四区| 国产精品亚洲av一区麻豆| 成人特级黄色片久久久久久久| 亚洲狠狠婷婷综合久久图片| 久久午夜综合久久蜜桃| 亚洲精品一区av在线观看| 丝袜人妻中文字幕| 国产午夜精品论理片| 91麻豆精品激情在线观看国产| 99久久国产精品久久久| 变态另类成人亚洲欧美熟女| 成在线人永久免费视频| 九色国产91popny在线| 亚洲成人精品中文字幕电影| tocl精华| 亚洲国产精品合色在线| 国产午夜精品久久久久久| 精品99又大又爽又粗少妇毛片 | 制服丝袜大香蕉在线| 高清在线国产一区| 神马国产精品三级电影在线观看| 久久天躁狠狠躁夜夜2o2o| 亚洲国产欧美网| 一进一出好大好爽视频| 欧美丝袜亚洲另类 | 午夜福利免费观看在线| 级片在线观看| 2021天堂中文幕一二区在线观| 国产亚洲av嫩草精品影院| 不卡av一区二区三区| 国产黄a三级三级三级人| 国产成人啪精品午夜网站| 亚洲国产欧美网| 国产黄片美女视频| 久久中文看片网| 在线免费观看不下载黄p国产 | 国产精品电影一区二区三区| 午夜福利18| 日本一二三区视频观看| 一区二区三区国产精品乱码| 国产精品香港三级国产av潘金莲| 很黄的视频免费| 网址你懂的国产日韩在线| 亚洲中文av在线| 日韩欧美国产在线观看| 99国产极品粉嫩在线观看| 91av网站免费观看| 免费在线观看日本一区| 国产一区二区在线观看日韩 | 黄色女人牲交| 无人区码免费观看不卡| 国内久久婷婷六月综合欲色啪| 成人av一区二区三区在线看| 精品久久久久久久末码| 色综合欧美亚洲国产小说| 成年女人看的毛片在线观看| 久久中文字幕人妻熟女| 免费观看精品视频网站| 在线永久观看黄色视频| 99久久综合精品五月天人人| 精品人妻1区二区| 午夜福利在线观看吧| 欧美乱色亚洲激情| 国产一区二区在线av高清观看| av黄色大香蕉| 久久亚洲精品不卡| 欧美一区二区国产精品久久精品| 国产在线精品亚洲第一网站| 国产精品香港三级国产av潘金莲| 欧美色欧美亚洲另类二区| 免费在线观看影片大全网站| 亚洲欧美日韩卡通动漫| 人妻丰满熟妇av一区二区三区| avwww免费| 国产成人av教育| 黑人巨大精品欧美一区二区mp4| 国产亚洲欧美在线一区二区| 日本三级黄在线观看| 国产乱人视频| 18美女黄网站色大片免费观看| 欧美日本视频| 国产精品免费一区二区三区在线| 亚洲精品456在线播放app | 欧美性猛交╳xxx乱大交人| 午夜精品一区二区三区免费看| 精品国产三级普通话版| 国产成人系列免费观看| 日日夜夜操网爽| 草草在线视频免费看| 国产69精品久久久久777片 | 国产精品久久久人人做人人爽| 天堂动漫精品| 人妻久久中文字幕网| 少妇人妻一区二区三区视频| 99在线视频只有这里精品首页| 亚洲av五月六月丁香网| 黄色视频,在线免费观看| 天堂影院成人在线观看| 亚洲av成人不卡在线观看播放网| 日韩高清综合在线| 黑人操中国人逼视频| 久久久成人免费电影| 性欧美人与动物交配| 午夜福利在线观看免费完整高清在 | 在线观看一区二区三区| 免费人成视频x8x8入口观看| 国内精品久久久久精免费| 在线视频色国产色| 老鸭窝网址在线观看| 少妇熟女aⅴ在线视频| 精品久久久久久久毛片微露脸| 国产久久久一区二区三区| 男女午夜视频在线观看| 午夜激情欧美在线| 亚洲 欧美 日韩 在线 免费| 欧美在线一区亚洲| 亚洲国产日韩欧美精品在线观看 | 黄色视频,在线免费观看| 亚洲在线自拍视频| 91麻豆精品激情在线观看国产| 成年女人毛片免费观看观看9| 国产午夜福利久久久久久| 搡老妇女老女人老熟妇| 一级毛片精品| 国产91精品成人一区二区三区| 色综合欧美亚洲国产小说| 国产精品久久久久久亚洲av鲁大| 好男人在线观看高清免费视频| 一二三四在线观看免费中文在| 日本免费a在线| 看片在线看免费视频| 长腿黑丝高跟| 一本精品99久久精品77| 男人舔女人的私密视频| 国产精品女同一区二区软件 | 亚洲欧美日韩无卡精品| 精品久久久久久成人av| 欧美黑人欧美精品刺激| 国产亚洲精品综合一区在线观看| 欧美zozozo另类| 在线观看66精品国产| 两性夫妻黄色片| 一进一出抽搐gif免费好疼| 国产日本99.免费观看| 91av网一区二区| 日本在线视频免费播放| 久久久水蜜桃国产精品网| 少妇的逼水好多| 国内毛片毛片毛片毛片毛片| 一二三四在线观看免费中文在| 欧美xxxx黑人xx丫x性爽| 国产极品精品免费视频能看的| 18禁黄网站禁片午夜丰满| 国产免费av片在线观看野外av| 亚洲男人的天堂狠狠| 国产精品一及| 成人精品一区二区免费| 男女之事视频高清在线观看| 99视频精品全部免费 在线 | 18禁黄网站禁片免费观看直播| 国产激情久久老熟女| 97超视频在线观看视频| 99久久99久久久精品蜜桃| 两个人视频免费观看高清| 国产亚洲av高清不卡| 在线观看日韩欧美| 丁香欧美五月| 最新美女视频免费是黄的| 成年女人永久免费观看视频| 中文资源天堂在线| 麻豆av在线久日| 成人av一区二区三区在线看| 欧美中文日本在线观看视频| 精品国产美女av久久久久小说| 色尼玛亚洲综合影院| 午夜免费观看网址| 丁香六月欧美| 最近最新中文字幕大全免费视频| 熟女电影av网| 免费在线观看成人毛片| 草草在线视频免费看| 中文字幕久久专区| 性色av乱码一区二区三区2| 国产精品99久久99久久久不卡| 欧美色视频一区免费| 最好的美女福利视频网| 久久久久精品国产欧美久久久| 九色成人免费人妻av| 无限看片的www在线观看| 欧美日韩乱码在线| 一区二区三区激情视频| 亚洲成a人片在线一区二区| 国产精品av视频在线免费观看| 国产探花在线观看一区二区| 搡老妇女老女人老熟妇| 精品国产超薄肉色丝袜足j| 久久久久性生活片| 免费av毛片视频| 成人精品一区二区免费| 男女视频在线观看网站免费| 老熟妇仑乱视频hdxx| 99久国产av精品| 亚洲,欧美精品.| 亚洲成av人片在线播放无| 男女视频在线观看网站免费| 精品电影一区二区在线| 精品久久久久久成人av| 亚洲精品国产精品久久久不卡| 在线十欧美十亚洲十日本专区| a在线观看视频网站| 国产精品影院久久| 久久久久精品国产欧美久久久| 国产亚洲欧美在线一区二区| 免费无遮挡裸体视频| 别揉我奶头~嗯~啊~动态视频| 日韩免费av在线播放| 可以在线观看的亚洲视频| 又粗又爽又猛毛片免费看| 色综合欧美亚洲国产小说| 国产成年人精品一区二区| 香蕉久久夜色| 久久久久九九精品影院| 久久久久久九九精品二区国产| 成人av在线播放网站| 999久久久国产精品视频| 欧美中文日本在线观看视频| 美女大奶头视频| 别揉我奶头~嗯~啊~动态视频| 欧洲精品卡2卡3卡4卡5卡区| 欧美高清成人免费视频www| 波多野结衣高清作品| 国产激情偷乱视频一区二区| 亚洲国产欧洲综合997久久,| 高清在线国产一区| 亚洲国产看品久久| 国产精品亚洲av一区麻豆| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产中文字幕在线视频| 国产精品久久久人人做人人爽| 欧美黄色淫秽网站| 国产精品久久久久久亚洲av鲁大| 99久久精品国产亚洲精品| 国产男靠女视频免费网站| 97超级碰碰碰精品色视频在线观看| 亚洲专区中文字幕在线| 一级作爱视频免费观看| 好男人电影高清在线观看| 精品无人区乱码1区二区| 亚洲乱码一区二区免费版| 老司机午夜福利在线观看视频| 色噜噜av男人的天堂激情| 黄色丝袜av网址大全| 丰满人妻一区二区三区视频av | 精品久久久久久,| 午夜激情欧美在线| 久久欧美精品欧美久久欧美| 亚洲成a人片在线一区二区| 亚洲成人精品中文字幕电影| 国产亚洲精品久久久com| 精品国产超薄肉色丝袜足j| 成年人黄色毛片网站| 人妻夜夜爽99麻豆av| 亚洲人成伊人成综合网2020| 美女午夜性视频免费| 欧美大码av| 床上黄色一级片| 看免费av毛片| 综合色av麻豆| 久久久久免费精品人妻一区二区| 亚洲国产精品成人综合色| 欧美日韩瑟瑟在线播放| 嫩草影院入口| 欧美一区二区精品小视频在线| 欧美3d第一页| 国产精品久久久久久精品电影| 美女被艹到高潮喷水动态| 我要搜黄色片| 久久国产精品影院| 在线观看日韩欧美| 丁香欧美五月| 久久久久久久久免费视频了| 成在线人永久免费视频| 国产乱人伦免费视频| 国产淫片久久久久久久久 | 免费在线观看成人毛片| 国产精品久久久av美女十八| 国产av一区在线观看免费| АⅤ资源中文在线天堂| 免费一级毛片在线播放高清视频| 午夜亚洲福利在线播放| 欧美日韩综合久久久久久 | 欧美成人免费av一区二区三区| 亚洲激情在线av| 午夜精品一区二区三区免费看| 99国产综合亚洲精品| 久久国产精品影院| 国产麻豆成人av免费视频| 听说在线观看完整版免费高清| 亚洲五月婷婷丁香| 此物有八面人人有两片| 成人三级做爰电影| 宅男免费午夜| 真人做人爱边吃奶动态| 69av精品久久久久久| 国产一级毛片七仙女欲春2| 日本黄色片子视频| 亚洲激情在线av| 日韩欧美在线乱码| 午夜福利在线观看吧| 18禁观看日本| 久久精品国产99精品国产亚洲性色| 国产精品久久久人人做人人爽| 国产精品野战在线观看| 欧美在线一区亚洲| 亚洲成人久久性| 亚洲欧美精品综合久久99| 久久久久久久久久黄片| 中文字幕高清在线视频| 国产高清有码在线观看视频| 精品国产超薄肉色丝袜足j| 久久久国产成人精品二区|