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

    改善八里湖水質(zhì)的最佳引水流量研究

    2022-10-27 09:42:16董文遜張艷軍周小國陳秀篁羅蘭
    中國農(nóng)村水利水電 2022年10期
    關(guān)鍵詞:九江市效益水體

    張 俊,董文遜,張艷軍,周小國,陳秀篁,羅蘭

    (1.長江生態(tài)環(huán)保集團(tuán)有限公司,湖北武漢 430062;2.中國長江三峽集團(tuán)有限公司,北京 100038;3.武漢大學(xué)水資源與水電工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,湖北武漢 430072)

    0 引言

    隨著城市化的進(jìn)程、人類社會的進(jìn)步,人與自然的矛盾進(jìn)一步凸顯,城市湖泊的水環(huán)境問題也越來越嚴(yán)重。為了治理城市湖泊,提高水體流速和更新速度,增強(qiáng)其生態(tài)系統(tǒng)穩(wěn)定性,除了削減入湖污染物以外,常見的工程措施是修建引水工程。通過引水工程來改善水環(huán)境的嘗試,可以追溯到20 世紀(jì)70 年代[1]。此后亦有諸多成功范例,例如通過對美國華盛頓州東部的Moses 湖修建引水工程,成功將Moses 湖的營養(yǎng)狀態(tài)由富營養(yǎng)化轉(zhuǎn)變?yōu)檩p度富營養(yǎng)化[2];針對水體缺乏活力等問題,為太湖修建了引水工程,取得了不錯(cuò)的治理成果和經(jīng)濟(jì)效益[3,4]引水工程具有見效快、改善明顯、能充分利用現(xiàn)有工程等特點(diǎn),在世界各地均取得了較好的成果。

    修建引水工程對湖泊的影響,歷來也多有研究。通過對鄭州鄭東新區(qū)的龍湖引水工程進(jìn)行研究,探索了引水工程對水體更新能力的影響[5,6];利用回歸神經(jīng)網(wǎng)絡(luò)模型的訓(xùn)練,成功預(yù)測了太湖公湖灣工程引水期間的水質(zhì)參數(shù)[7];對常熟尚湖進(jìn)行跟蹤研究,探明了引水對小型通江湖泊富營養(yǎng)化問題的改善效果等[8]。過去針對湖泊引水工程的研究,都在一定程度上肯定了引水工程對改善城市湖泊水環(huán)境的作用。為了制定合適的引水工程方案,學(xué)者們往往運(yùn)用水動(dòng)力、水質(zhì)模型等進(jìn)行計(jì)算。例如針對武漢東湖生態(tài)環(huán)境穩(wěn)定性較差的問題,建立了水動(dòng)力和水質(zhì)模型,確定了最優(yōu)引水方案[9];使用基于DEM 的水量水質(zhì)模型算法,針對湖北黃石磁湖污染,設(shè)計(jì)并模擬了多種水體連通方案,印證了引水工程對水環(huán)境的改善效果[10];使用MIKE21 軟件建立數(shù)學(xué)模型,探究了不同引水方案對青海省樂都主城段南側(cè)人工湖換水率的影響等[11]。

    八里湖作為九江市受人類活動(dòng)影響最大的城市湖泊之一,人與自然的矛盾突出,水質(zhì)問題也極為嚴(yán)重。八里湖長期保持V 類水甚至劣V 類水水平,水質(zhì)情況堪憂;部分水域?yàn)樗浪畢^(qū)、滯水區(qū),黑臭問題難以解決;近年來人類城市的高速發(fā)展新增了大量污染源,工業(yè)、生活廢水往往就近排入河道,缺乏監(jiān)管;2020年夏季在十里河入八里湖湖口部分,更是出現(xiàn)了大規(guī)模藍(lán)藻爆發(fā)問題。研究城市化對城市湖泊的影響,分析如何治理與修復(fù)城市湖泊,協(xié)調(diào)人與城市湖泊的關(guān)系,歷來是學(xué)者比較關(guān)心的問題[12-14],修建引水工程便是常見的工程治理措施。

    針對目前九江市面臨的一系列水環(huán)境水生態(tài)問題,尤其是八里湖的水質(zhì)問題,結(jié)合資料以及前期實(shí)地考察勘探工作,綜合考慮九江市獨(dú)特的“山-城-河-湖-江”復(fù)雜城市水系統(tǒng)特性,本文將基于DEM(Digital Elevation Model,數(shù)字高程模型)的水量水質(zhì)模型算法創(chuàng)新性地運(yùn)用于八里湖的水質(zhì)模擬,并針對八里湖氮磷超標(biāo)嚴(yán)重的問題,以總磷為主要模擬指標(biāo),設(shè)計(jì)多種工況的引水工程,“以清沖污”,分析改善八里湖水質(zhì)的最佳引水流量,評估其直接效益、間接效益和凈效益,從工程的角度,為九江市生態(tài)環(huán)境治理提供有效科學(xué)的建議。

    1 研究區(qū)域

    1.1 研究區(qū)域概況

    九江市位于江西省北部,長江中游,地理?xiàng)l件復(fù)雜,南靠廬山、東臨鄱陽湖、北接長江。全境南北寬140 km,東西長270 km,地理坐標(biāo)為113°57'E~116°53'E,28°47'N~30°06'N,市轄區(qū)面積1 352 km2。九江市的山地產(chǎn)流迅速匯流入各小河流,再流入城市內(nèi)湖,經(jīng)湖泊出口處順接的河道或泵閘匯入長江,形成了獨(dú)特的“山-城-河-湖-江”復(fù)雜城市水系統(tǒng)。長江九江中心城區(qū)段常年洪水位18~18.25 m(7-8 月)。每年1-2 月最低水位約為9.4 m左右,最大流量為8.17 萬m3/s,年均流量為2.43萬m3/s,水量充沛。八里湖、龍開河等可以通過自由出流匯入長江,但每年長江水位約有1~2個(gè)月的時(shí)間高于湖面。八里湖與長江相聯(lián),也能對長江起到一定的調(diào)蓄作用。九江市八里湖的具體地理位置如圖1所示。

    圖1 九江市地理位置Fig.1 Location of Jiujiang City

    八里湖地處九江市區(qū)西南部,東面與十里河以及九江經(jīng)濟(jì)技術(shù)開發(fā)區(qū)相接,西面在賽湖壩頭處與賽城湖相鄰,南面為廬山風(fēng)景區(qū)以及蛟灘河、截流河、沙河等水系,北面通過新開河、龍開河與長江相連。八里湖是長江南岸直入長江的湖泊,流域面積273 km2,現(xiàn)狀常水位為16.5 m,相應(yīng)湖區(qū)面積約18 km2,主要水源是廬山北部的來水,湖水主要經(jīng)新開河由八里湖閘和八里湖泵站排入長江。八里湖流域內(nèi)主要水系有十里河、蛟灘河、沙河等。八里湖多年平均徑流量2.5 億m3,五十年一遇防洪水位為18.97 m,相應(yīng)庫容0.68 億m3,主要依靠八賽隔堤、八里湖堤、向陽堤、十里河堤等堤防防御八里湖洪水。

    1.2 研究數(shù)據(jù)

    研究數(shù)據(jù)主要包括水下地形數(shù)據(jù)和水質(zhì)數(shù)據(jù)。水樣測點(diǎn)的分布如圖2 所示。整體上點(diǎn)位分布比較均勻,能輻射北半湖大部分區(qū)域;對八里湖各個(gè)出入口(如賽城湖、龍開河、十里河等)均有布置測點(diǎn),能兼顧重點(diǎn)區(qū)域。水質(zhì)監(jiān)測數(shù)據(jù)如表1所示。

    表1 八里湖實(shí)測水質(zhì)濃度 mg/LTab.1 Observed water quality concentration of Bali Lake

    圖2 八里湖水樣檢測分布示意圖Fig.2 Schematic diagram of testing distribution of Bali Lake water samples

    2 研究方法

    2.1 二維水動(dòng)力水質(zhì)模型

    2.1.1 水動(dòng)力模型

    基于笛卡爾坐標(biāo)系下的平面二維淺水控制方程組,構(gòu)建八里湖二維水動(dòng)力水質(zhì)模型的水動(dòng)力模擬部分。

    x方向上的動(dòng)量方程如下式所示:

    y方向上的動(dòng)量方程如下式所示:

    式中:h代表研究區(qū)域八里湖的水深;t為模型模擬的時(shí)間步長;u,v代表平面二維坐標(biāo)系橫向和縱向的流速;εx、εy分別為x、y方向的渦動(dòng)黏滯系數(shù);g為重力常數(shù);n為研究區(qū)域八里湖的糙率;q為八里湖的支流入流流量。

    2.1.2 污染物遷移轉(zhuǎn)化模型

    基于平面二維水質(zhì)對流擴(kuò)散方程(平面二維水質(zhì)遷移轉(zhuǎn)化基本方程),構(gòu)建八里湖二維水動(dòng)力水質(zhì)模型的水質(zhì)模擬。其表達(dá)式參照下式:

    式中:∑Si為湖泊水體污染物的源項(xiàng),包括內(nèi)部源漏項(xiàng)和外部源漏項(xiàng);Ex為x方向的分子擴(kuò)散系數(shù)、紊動(dòng)擴(kuò)散系數(shù)和離散系數(shù)之和;Ey為y方向的分子擴(kuò)散系數(shù)、紊動(dòng)擴(kuò)散系數(shù)和離散系數(shù)之和。

    2.2 成本和效益評估方法與水質(zhì)評價(jià)標(biāo)準(zhǔn)

    隨著城市建設(shè)高速發(fā)展,人類活動(dòng)對城市湖泊的影響越來越大。九江市的水體普遍存在水系溝通不良、水系萎縮、引排不暢、河床淤積等問題,除了控制污染源頭、減少排污量以外,還可選用引水工程改善環(huán)境。引水工程是通過對引水的時(shí)空調(diào)度和控制來綜合配置水資源,從而達(dá)到改善湖區(qū)水質(zhì)的目的[15],對于水資源調(diào)配、水質(zhì)改善、水災(zāi)防治等有著重要的意義,有助于提高區(qū)域經(jīng)濟(jì)發(fā)展水平。因此本文結(jié)合八里湖二維水動(dòng)力水質(zhì)模型,研究不同引水流量對湖泊水質(zhì)改善效果的影響,并基于成本和效益評估方法,分析八里湖最佳引水流量。

    引水工程的凈效益為直接效益減去成本,即:

    式中:E指工程的凈效益,元;Ed指直接效益,主要包括引水工程建設(shè)后減少的湖泊污水治理費(fèi)Ea以及增加的農(nóng)業(yè)效益En,元;Cp指泵站工作的運(yùn)行成本,元。

    2.2.1 直接效益

    引水工程實(shí)際上是將重度污染區(qū)的污染物通過水流,轉(zhuǎn)移到其他水體之中,污染物本質(zhì)上并未減少,但只需要控制好引水工程的流量,就可以保證引水工程的攜污廢水不會導(dǎo)致長江九江段(納污水體)超過其水質(zhì)目標(biāo)。因此為了研究方便,不考慮下游水體容納污染物帶來的經(jīng)濟(jì)損失。九江中心城區(qū)面積共554 km2,農(nóng)業(yè)用地136.56 km2,各類水利工程大部分取水用于農(nóng)業(yè)灌溉,因此本研究中,引水工程的直接效益主要來自于減少的污染物處理費(fèi)用以及農(nóng)業(yè)效益。

    (1)湖泊污水治理效益。湖泊污水治理效益可以通過下式計(jì)算,即:

    式中:M0指基于初始水質(zhì)條件通過八里湖二維水動(dòng)力水質(zhì)模型計(jì)算得到的污染物初始總量,g;M1指污染物的引水后總量,基于不同引水工況,由八里湖二維水動(dòng)力水質(zhì)模型模擬計(jì)算獲得,g;Pa指處理1 g污染物的成本,元/g,根據(jù)2014年國家發(fā)展和改革委員會發(fā)布的排污費(fèi)征收標(biāo)準(zhǔn)管理辦法[16],可知,污水排污費(fèi)應(yīng)該按照排污者排放污染物的種類和數(shù)量,以污染當(dāng)量計(jì)征,因此確定計(jì)算公式如下:

    式中:R指污染當(dāng)量的征收標(biāo)準(zhǔn),單位為元;K是污染物的污染當(dāng)量,kg;b1為換算系數(shù),這里取1 000。針對八里湖嚴(yán)重的水質(zhì)超標(biāo)問題,本研究中的典型污染物選定為總磷(TP)??偭讓?yīng)的征收標(biāo)準(zhǔn)為每污染當(dāng)量0.7元,對應(yīng)的污染當(dāng)量值為0.25 kg,因此Pa的值取0.002 8元/g。

    (2)農(nóng)業(yè)效益。農(nóng)業(yè)效益主要來自于八里湖引水工程運(yùn)行后水體劣V 類面積減小、污水灌溉比例減小帶來的種植業(yè)產(chǎn)量增加。農(nóng)業(yè)效益可按下式計(jì)算:

    式中,En為農(nóng)業(yè)效益,元;A0、A1分別為引水工程前、后劣V 類水體比例,%;B為增產(chǎn)幅度,%;C為九江市濂溪區(qū)和潯陽區(qū)的農(nóng)業(yè)產(chǎn)值,元。

    2.2.2 間接效益

    水環(huán)境污染除了會對水質(zhì)、農(nóng)業(yè)灌溉造成直接影響外,還會對周圍居民的健康狀況、水體觀賞價(jià)值、水體自凈效益等方面造成間接影響。因此本研究中,引水工程改善水質(zhì)的間接效益主要包括水質(zhì)改善帶來的社會效益(健康、觀賞、宣傳效益)和生態(tài)效益(水體自凈效益)兩個(gè)方面,可以起到輔助確定最佳引水流量和最佳工況的作用。

    (1)健康效益。水環(huán)境污染會造成人類健康損失,因此水質(zhì)改善的健康收益可以用患病造成的損失進(jìn)行計(jì)算。基于美國環(huán)保署提出的人體健康風(fēng)險(xiǎn)評價(jià)模型(HHRA)和國內(nèi)關(guān)于垃圾填埋場對居民健康風(fēng)險(xiǎn)評價(jià)的研究等[17],可以確定公式如下:

    式中:E1為健康收益,元;Eyl為人均醫(yī)療費(fèi),元/(人?a),可參考北京市癌癥平均住院費(fèi)用[18];Egz為年均工資,元/(人?a),可參考2020年江西九江平均工資為5 847 元/(人?月);Twg為誤工年數(shù),年;Rc為暴露人群經(jīng)飲水途徑的人均致癌年風(fēng)險(xiǎn);P為研究區(qū)域總?cè)藬?shù),以濂溪區(qū)和潯陽區(qū)總?cè)藬?shù)計(jì)算,共計(jì)58.59 萬人(2019年末)。

    Rc可以按照下式計(jì)算:

    式中:q(人)為致癌強(qiáng)度系數(shù),(kg?d)/mg;ADD為暴露劑量,mg/(kg?d),由于九江自來水廠主要從長江提水,因此這里不計(jì)算經(jīng)口暴露劑量,只計(jì)算皮膚吸收劑量。ADD可以按照下式計(jì)算:

    式中:CW為水中污染物濃度,mg/L;SA為皮膚接觸表面積,m2,可參考美國EPA 的手冊[19]等,取均值1.36 m2;PC為化學(xué)物質(zhì)皮膚滲透常數(shù),cm/h,可取均值0.003 2 cm/h;ET為暴露時(shí)間,以0.2 h/d 計(jì)算;EF為暴露頻率,以365 d/a 計(jì)算;ED為暴露持續(xù)時(shí)間;b2為換算系數(shù),這里取0.01;BW為人體質(zhì)量,以均值66.5 kg計(jì)算;AT為平均時(shí)間。

    (2)觀賞效益。水質(zhì)改善的觀賞效益可用下式計(jì)算:

    式中:E2為水體觀賞效益,元;A為八里湖水域面積,為18 km2;n為水體美學(xué)價(jià)值的當(dāng)量因子,參考關(guān)于生態(tài)系統(tǒng)服務(wù)價(jià)值化方法的研究[20]可取1.89;D為九江市單位面積凈產(chǎn)出損失,元/km2,可以通過查閱《九江統(tǒng)計(jì)年鑒——2020》等獲得;其他參數(shù)同上。

    (3)宣傳效益。宣傳效益來自于水質(zhì)提升的信息傳遞附加價(jià)值,可以近似使用互聯(lián)網(wǎng)上關(guān)于水污染信息的關(guān)注程度進(jìn)行折算,參考公式如下:

    式中:E3為宣傳效益,元;Pw為研究區(qū)域內(nèi)網(wǎng)民人數(shù),人;F1為網(wǎng)絡(luò)信息的經(jīng)濟(jì)轉(zhuǎn)化率,元/h,F(xiàn)2網(wǎng)民在新聞上花費(fèi)時(shí)間的百分比;F3為水污染相關(guān)信息的網(wǎng)絡(luò)關(guān)注度,可通過查詢百度指數(shù)獲得;Tw為上網(wǎng)時(shí)間,h,通過查閱《第47次中國互聯(lián)網(wǎng)發(fā)展?fàn)顩r統(tǒng)計(jì)報(bào)告》可知,截止2020 年12月,我國網(wǎng)民的周上網(wǎng)時(shí)間為26.2 h。

    網(wǎng)民人數(shù)可按照下式計(jì)算:

    式中:PC、PN分別為城鎮(zhèn)、農(nóng)村人口,可以通過查閱《九江統(tǒng)計(jì)年鑒——2020》等獲得;CC、CN分別為城鎮(zhèn)、農(nóng)村地區(qū)網(wǎng)民比例。

    通過查閱《第47次中國互聯(lián)網(wǎng)發(fā)展?fàn)顩r統(tǒng)計(jì)報(bào)告》得知,截止2020 年12月,我國城鎮(zhèn)、農(nóng)村地區(qū)互聯(lián)網(wǎng)普及率分別為79.8%和55.9%,可以此作為城鎮(zhèn)、農(nóng)村地區(qū)網(wǎng)民比例進(jìn)行計(jì)算。

    (4)水體自凈效益。八里湖流域因?yàn)橐挪粫?,容易?dǎo)致水體活力下降,形成死水,進(jìn)一步降低水體的自凈能力[21]。因此水體自凈效益即為八里湖污染物減少帶來的水體自凈能力增加值,可按下式計(jì)算:

    式中:Fj為自凈系數(shù)[22];Fin為入河污染物的太陽能值轉(zhuǎn)換率;EDR為能值貨幣比率,sej/元,可以用九江市全年能值投入總量除以當(dāng)年生產(chǎn)總值GDP獲得;其他參數(shù)意義同上。

    2.2.3 運(yùn)行成本

    實(shí)際引水工程從建成到運(yùn)行需要涉及多個(gè)方面的成本投入,如招標(biāo)施工費(fèi)用、維護(hù)運(yùn)行費(fèi)用等。為簡便起見,本研究中考慮的成本主要指泵站工作時(shí)產(chǎn)生的運(yùn)行成本,可以使用下式計(jì)算,即:

    式中:Pp指泵站單位引水量產(chǎn)生的運(yùn)行費(fèi)用,元/m3,已知八里湖泵站規(guī)模為51 m3/s,裝機(jī)為3 000 kW,根據(jù)江西省九江市電量分檔可知月均350 kWh以上的店家為0.90 元/kWh,因此估算得到泵抽水的價(jià)格Pp為0.014 7 元/m3;Q指引水工程的流量,m3/s;T是引水工程運(yùn)行持續(xù)的時(shí)間,s;其他符號意義同上。

    2.2.4 水質(zhì)評價(jià)標(biāo)準(zhǔn)

    根據(jù)《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)GB/T 3838-2002》(如表2)規(guī)定的標(biāo)準(zhǔn)對水質(zhì)進(jìn)行評價(jià)。其中,劣V 類水指,實(shí)測污染物濃度高于《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)GB/T 3838-2002》中V 類水的濃度標(biāo)準(zhǔn)。

    表2 地表水環(huán)境質(zhì)量基本項(xiàng)目標(biāo)準(zhǔn)限值 mg/LTab.2 Standard limits for basic items of surface water environmental quality

    3 工況情景設(shè)置與結(jié)果分析

    3.1 水體連通工況

    考慮到八里湖總磷背景濃度為劣V類,而長江九江段以及賽城湖的水質(zhì)較好,因此考慮通過八賽樞紐泵站,從賽城湖、新開河引清水進(jìn)入八里湖,“以清沖污”,最后通過改建龍開河閘,改暗河暗渠為明渠,增強(qiáng)過水能力,使污水加速排出,達(dá)到凈化湖泊水質(zhì)、保護(hù)水體功能的目的。具體的引水路線圖如圖3所示。

    圖3 八里湖引水工程路線圖Fig.3 Route map of Bali Lake water diversion project

    為充分利用現(xiàn)有工程,基于賽城湖排澇站(位于八里湖賽城湖聯(lián)通處)設(shè)計(jì)排澇流量為104 m3/s 的實(shí)際情況,設(shè)置了20、60、100 m3/s三個(gè)不同擋位的賽城湖引水流量;基于八里湖泵站(位于新開河上)設(shè)計(jì)流量為50 m3/s的實(shí)際情況,設(shè)置了10、30、50 m3/s三個(gè)不同擋位的新開河引水流量。兩處泵站的3種不同擋位引水流量兩兩結(jié)合,加上無引水工程時(shí)的工況,一共組合出以下10 種工況,按照總引水流量從小到大進(jìn)行排列,最小總引水流量為30 m3/s,最大總引水流量為150 m3/s,如表3所示。

    表3 計(jì)算工況表Tab.3 Calculation condition table

    3.2 二維水動(dòng)力水質(zhì)模型的率定與驗(yàn)證

    選取總磷和COD 對二維水動(dòng)力水質(zhì)模型模型進(jìn)行率定和驗(yàn)證,率定期時(shí)段設(shè)置為從2020年7月1日開始,到2020年7月8日結(jié)束。率定期的模擬結(jié)果如表4、5所示。

    從表4、5 中可以看出大部分點(diǎn)位的模擬誤差均在25%以內(nèi),誤差絕對值的均值也分別保持在17.97%和21.10%,僅有少量點(diǎn)位因環(huán)境波動(dòng)產(chǎn)生了偏差,因此該模型能夠較好地模擬八里湖的水質(zhì)情況。

    表4 八里湖水環(huán)境模型率定成果(7月8日,總磷)Tab.4 Calibration results of Bali Lake water environment model(July 8th,total phosphorus)

    表5 八里湖水環(huán)境模型率定成果(7月8日,COD)Tab.5 Calibration results of Bali Lake water environment model(July 8th,COD)

    3.3 最佳引水流量分析

    3.3.1 水動(dòng)力改善效果分析

    根據(jù)八里湖水動(dòng)力水質(zhì)模型,可以對引水工程建設(shè)后的流場進(jìn)行模擬,如圖4 所示。由流場分布圖可以直觀地分析八里湖引水工程投產(chǎn)使用后流場的空間分布特征:引水工程的建設(shè)改變了原始的水體流動(dòng)方向,原本蛟灘河與截流河的來水會向北流動(dòng),最終藉由新開河匯入長江,現(xiàn)在變?yōu)榕c賽城湖、新開河匯合,共同通過龍開河匯入長江。水動(dòng)力改善效果主要可以從平均流速、最大流速和滯水區(qū)面積三個(gè)指標(biāo)進(jìn)行評價(jià),其中滯水區(qū)是指流速小于0.000 6 m/s 的水體。由于九江市城市湖泊普遍存在水體動(dòng)力不足、河道淤塞、水體易黑易臭等問題,通過對以上3 個(gè)水動(dòng)力指標(biāo)進(jìn)行評價(jià),可以直觀反映引水工程對環(huán)境的改善效果。

    圖4 八里湖流場分布情況Fig.4 Distribution of flow field in Bali Lake

    3.3.2 水質(zhì)改善效果分析

    引水工程在8 天時(shí)能達(dá)到削減劣V 類水體面積的最大效果,因此繪制各種工況下引水8 天時(shí)的總磷濃度分布圖,如圖5所示。

    由圖5 可以看出,修建引水工程能對環(huán)境起到較大的改善作用,IV 類、V 類水體面積顯著提升;盡管引水工程對八里湖湖心的影響較小,但是八里湖北半湖的水質(zhì)質(zhì)量明顯好轉(zhuǎn),重點(diǎn)區(qū)域(如十里河入湖口附近)的水質(zhì)狀況好轉(zhuǎn);引水口附近的水質(zhì)改善效果最明顯,距離引水口越遠(yuǎn),水質(zhì)越差;修建引水工程后,各工況下的濃度整體分布差別不大,南半湖濃度分布幾乎完全相同,但引水流量越大,八里湖北半湖的非劣V 類水體占比越大;引水工程修建改變了原本的水體循環(huán)路徑,水流難以經(jīng)過的地方新增了小部分死角,以引水總流量相等的工況7、8為例,因?yàn)楣r8 兩個(gè)引水口的流量差別過大(新開河10 m3/s、賽城湖100 m3/s),在新開河西南方向出現(xiàn)了總磷濃度偏高的死角,所以在選擇工況時(shí)應(yīng)盡量保證引水口之間的濃度接近??傮w而言,八里湖引水工程能夠有效改善重點(diǎn)區(qū)域的水質(zhì)超標(biāo)問題,改善八里湖水環(huán)境,優(yōu)化八里湖水生態(tài),充分發(fā)揮九江水優(yōu)勢,對未來八里湖的治理起到了一定的參考價(jià)值。

    圖5 不同工況下總磷濃度分布圖(單位:mg/L)Fig.5 Distribution map of total phosphorus concentration under different working conditions

    3.3.3 最佳引水流量分析

    表6 給出了幾個(gè)關(guān)鍵的評價(jià)指標(biāo)??梢钥闯觯谝? 天后,即使在引水總量30 m3/s 的情況下,也能達(dá)到每次77.42 萬元/次的直接效益,說明引水工程確實(shí)有助于八里湖水質(zhì)改善;由表6 可知隨著引水流量的增加,整體上直接效益是遞增的。但是與此同時(shí),對應(yīng)的運(yùn)行成本也隨著引水時(shí)間的增加而增加,最小為工況2 的30.48 萬元/次,最大為工況10 的152.41 萬元/次,因此凈效益的增加幅度越來越小,甚至減小。根據(jù)凈效益的計(jì)算結(jié)果,可以確定最佳工況為工況4,對應(yīng)的直接效益為130.10 萬元/次,運(yùn)行成本為71.12 萬元/次,凈效益為58.98 萬元/次。

    表6 各工況引水沖污的效益與花費(fèi)表Tab.6 Table of benefits and costs of water diversion and flushing under various working conditions

    基于2.2.2 的計(jì)算方法,可以求解出8 天時(shí)各工況下的間接效益,如圖6 所示。由圖可知,隨著引水流量的增加,各類間接效益以及總間接效益的值逐漸增加,工況4 的總間接效益為70.89 萬元,和其他工況差距不大。在四種間接效益中,水體自凈效益占比較大,說明引水工程的生態(tài)效益比較顯著。同時(shí)宣傳效益的占比也較大,說明引水工程能夠增加居民對水污染治理工程的認(rèn)同感,提升城市品質(zhì)和影響力。

    圖6 各工況間接效益情況Fig.6 Indirect benefit of each working condition

    4 結(jié)論

    本文針對九江市八里湖區(qū)域嚴(yán)重的水質(zhì)超標(biāo)問題,基于實(shí)地考察獲得的水下地形資料和水質(zhì)資料,建立了八里湖二維水動(dòng)力水質(zhì)模型,并將該模型運(yùn)用于設(shè)定的引水工程之中,分析了引水工程的水動(dòng)力和水質(zhì)改善效果,基于成本和效益評估方法確定了八里湖最佳引水流量,提出了相應(yīng)的環(huán)境治理建議。

    對10 種不同工況的水動(dòng)力改善效果和水質(zhì)改善效果的分析表明,在選擇合適工況時(shí)應(yīng)該充分考慮引水口的流速和地理位置影響;濃度分布比例的變化表明,引水工程影響下的各類水體面積占比變化趨勢相同,且隨著引水時(shí)間的增加劣V 類水體面積逐漸減小,但減小幅度逐步下降,最終在第8天的時(shí)候劣V 類水體面積占比基本達(dá)到最小值,因此引水時(shí)間在8 天左右最為合適。

    在引水工程改善效果研究的基礎(chǔ)上,對比各工況下污染物總量和效益隨時(shí)間變化的關(guān)系,在少許下降或者波動(dòng)之后,直接效益曲線呈現(xiàn)先增長再穩(wěn)定的趨勢,凈效益曲線呈現(xiàn)先增長再降低的趨勢,污染物總量一般在第7~8 d 即可達(dá)到最小值。因此以引水時(shí)間8天為例,分析了不同引水流量組合下,可獲得的直接效益、間接效益和凈效益。最佳引水流量的選擇應(yīng)該考慮凈效益較大,同時(shí)兼顧水動(dòng)力和水質(zhì)改善效果以及間接效益情況。因此選擇工況4 進(jìn)行引水最為合適,相應(yīng)的最佳引水流量為70 m3/s(即新開河引水50 m3/s、賽城湖引水20 m3/s),引水的凈效益為58.98萬元/次。

    因此,引水工程可以通過提升水體動(dòng)力、溝通水體調(diào)配污染物兩個(gè)方面對八里湖環(huán)境進(jìn)行治理。

    猜你喜歡
    九江市效益水體
    草粉發(fā)酵 喂羊效益高
    蓮魚混養(yǎng) 效益提高一倍
    九江市審計(jì)局開展“七個(gè)一”系列活動(dòng)慶“七一”
    農(nóng)村黑臭水體治理和污水處理淺探
    生態(tài)修復(fù)理念在河道水體治理中的應(yīng)用
    我的“自畫像”
    冬棚養(yǎng)蝦效益顯著,看技術(shù)達(dá)人如何手到“錢”來
    不老的情緣(合唱)
    ——江西省九江市潯陽區(qū)老年大學(xué)校歌
    果園有了“鵝幫工” 一舉多得效益好
    廣元:治理黑臭水體 再還水清岸美
    久久精品国产a三级三级三级| 亚洲精品亚洲一区二区| 欧美人与性动交α欧美精品济南到 | 日韩中文字幕视频在线看片| 一级毛片电影观看| 国产精品嫩草影院av在线观看| 久久久久久人妻| 日韩在线高清观看一区二区三区| 欧美亚洲 丝袜 人妻 在线| 十分钟在线观看高清视频www| 卡戴珊不雅视频在线播放| 秋霞在线观看毛片| 欧美变态另类bdsm刘玥| 亚洲精品久久久久久婷婷小说| 少妇的逼好多水| 欧美 亚洲 国产 日韩一| 亚洲图色成人| 熟女电影av网| 大码成人一级视频| 老司机影院成人| 精品亚洲成a人片在线观看| www.av在线官网国产| 如日韩欧美国产精品一区二区三区 | 国产老妇伦熟女老妇高清| av在线app专区| 欧美+日韩+精品| 蜜臀久久99精品久久宅男| 多毛熟女@视频| 女性被躁到高潮视频| 精品亚洲成国产av| 成人亚洲欧美一区二区av| 欧美bdsm另类| 久久国产亚洲av麻豆专区| 超碰97精品在线观看| 午夜福利,免费看| 18在线观看网站| 男女高潮啪啪啪动态图| 亚洲五月色婷婷综合| 精品一品国产午夜福利视频| 91午夜精品亚洲一区二区三区| 国产男女内射视频| 99热全是精品| 色5月婷婷丁香| 国产极品天堂在线| 啦啦啦中文免费视频观看日本| 妹子高潮喷水视频| 精品少妇黑人巨大在线播放| 26uuu在线亚洲综合色| 亚洲欧美日韩卡通动漫| 国产亚洲精品第一综合不卡 | 亚洲国产av影院在线观看| 91精品国产国语对白视频| 国产精品久久久久久久久免| 国产精品一区二区在线观看99| 亚洲欧美成人综合另类久久久| 少妇猛男粗大的猛烈进出视频| 精品99又大又爽又粗少妇毛片| 成年人免费黄色播放视频| 日本av免费视频播放| 国产精品女同一区二区软件| 久热久热在线精品观看| 夫妻性生交免费视频一级片| 国产成人免费观看mmmm| 中文字幕精品免费在线观看视频 | 寂寞人妻少妇视频99o| 久久久久视频综合| 久久久久国产精品人妻一区二区| 91精品一卡2卡3卡4卡| 久久久午夜欧美精品| 久久午夜综合久久蜜桃| 狠狠精品人妻久久久久久综合| 黄色配什么色好看| 亚洲欧美成人综合另类久久久| 亚洲成人手机| av黄色大香蕉| 制服丝袜香蕉在线| 欧美国产精品一级二级三级| 91精品国产九色| 人妻制服诱惑在线中文字幕| 日本vs欧美在线观看视频| 亚洲国产最新在线播放| 久久婷婷青草| 亚洲色图 男人天堂 中文字幕 | 国产精品久久久久成人av| 午夜福利视频在线观看免费| 成人毛片a级毛片在线播放| 少妇的逼好多水| 狂野欧美激情性bbbbbb| 亚洲精品美女久久av网站| 欧美人与性动交α欧美精品济南到 | 成人亚洲欧美一区二区av| 最近最新中文字幕免费大全7| 国产精品久久久久成人av| 黄色欧美视频在线观看| 91国产中文字幕| 精品少妇久久久久久888优播| 如何舔出高潮| 久久久久久久久久人人人人人人| 国产国拍精品亚洲av在线观看| 国产精品一区www在线观看| 18+在线观看网站| 只有这里有精品99| 亚洲欧美日韩卡通动漫| 中文字幕久久专区| av播播在线观看一区| 亚洲色图综合在线观看| 一级爰片在线观看| 大香蕉久久成人网| 欧美日韩国产mv在线观看视频| 丝袜喷水一区| 免费人妻精品一区二区三区视频| videossex国产| 乱码一卡2卡4卡精品| 制服丝袜香蕉在线| 人人妻人人添人人爽欧美一区卜| 国产精品蜜桃在线观看| 亚洲第一区二区三区不卡| 18禁裸乳无遮挡动漫免费视频| 美女国产高潮福利片在线看| 人妻 亚洲 视频| 亚洲精品久久久久久婷婷小说| 色视频在线一区二区三区| 中文字幕亚洲精品专区| 国产av码专区亚洲av| 高清av免费在线| 久久精品夜色国产| 国产成人精品在线电影| 九色成人免费人妻av| 天堂俺去俺来也www色官网| 免费看光身美女| 国产伦精品一区二区三区视频9| 免费观看的影片在线观看| h视频一区二区三区| 大香蕉久久成人网| 国产成人免费观看mmmm| 亚洲欧美日韩另类电影网站| 日本黄色日本黄色录像| 亚洲av男天堂| 国产精品国产三级国产专区5o| 一个人看视频在线观看www免费| 国产欧美日韩综合在线一区二区| 国产精品久久久久久精品电影小说| 少妇的逼水好多| 街头女战士在线观看网站| 日日爽夜夜爽网站| 日韩精品免费视频一区二区三区 | 亚洲怡红院男人天堂| 美女内射精品一级片tv| 在线免费观看不下载黄p国产| 老女人水多毛片| 涩涩av久久男人的天堂| 精品酒店卫生间| 国产熟女午夜一区二区三区 | 亚洲精品456在线播放app| 晚上一个人看的免费电影| 国产亚洲av片在线观看秒播厂| 五月开心婷婷网| 精品国产露脸久久av麻豆| 美女主播在线视频| 国产成人91sexporn| 国产精品免费大片| 欧美日韩国产mv在线观看视频| 国产亚洲一区二区精品| 亚洲精品乱码久久久v下载方式| 日本欧美视频一区| 大香蕉久久网| 日韩av在线免费看完整版不卡| 18禁裸乳无遮挡动漫免费视频| 中文字幕av电影在线播放| 少妇 在线观看| 日韩视频在线欧美| 夫妻性生交免费视频一级片| 日本午夜av视频| 久久久精品94久久精品| 国内精品宾馆在线| 久久精品国产鲁丝片午夜精品| 亚洲久久久国产精品| 母亲3免费完整高清在线观看 | 如何舔出高潮| 中国美白少妇内射xxxbb| 最新中文字幕久久久久| 亚洲国产色片| 国产一区二区三区综合在线观看 | 亚洲欧美清纯卡通| 亚洲无线观看免费| 国产男人的电影天堂91| 纯流量卡能插随身wifi吗| 午夜福利在线观看免费完整高清在| 国产精品一区二区在线不卡| 国产深夜福利视频在线观看| 汤姆久久久久久久影院中文字幕| 我要看黄色一级片免费的| 久久久欧美国产精品| 久久精品久久久久久噜噜老黄| 五月天丁香电影| 高清视频免费观看一区二区| 一边亲一边摸免费视频| 中文字幕精品免费在线观看视频 | 日韩人妻高清精品专区| 成人国产麻豆网| 国国产精品蜜臀av免费| 国产无遮挡羞羞视频在线观看| 亚洲精品美女久久av网站| 人人妻人人澡人人看| 亚洲欧美日韩卡通动漫| 极品少妇高潮喷水抽搐| 精品国产露脸久久av麻豆| 免费观看av网站的网址| 欧美日韩成人在线一区二区| 丰满迷人的少妇在线观看| 久久精品夜色国产| 日日摸夜夜添夜夜添av毛片| 制服诱惑二区| 国产日韩一区二区三区精品不卡 | 女的被弄到高潮叫床怎么办| 日韩欧美一区视频在线观看| 亚洲av日韩在线播放| 一本—道久久a久久精品蜜桃钙片| 丝袜美足系列| 18禁观看日本| 看免费成人av毛片| 街头女战士在线观看网站| 视频中文字幕在线观看| 少妇 在线观看| 国产精品人妻久久久久久| 亚洲精品乱码久久久久久按摩| 精品一区在线观看国产| 麻豆精品久久久久久蜜桃| 欧美激情 高清一区二区三区| 久久久久久久久大av| 久久人人爽人人片av| 欧美日韩综合久久久久久| 亚洲av电影在线观看一区二区三区| 国产日韩欧美亚洲二区| 3wmmmm亚洲av在线观看| 国产乱来视频区| 插逼视频在线观看| 一级二级三级毛片免费看| 久久久国产欧美日韩av| 国产国语露脸激情在线看| 最近的中文字幕免费完整| 夫妻性生交免费视频一级片| 日日爽夜夜爽网站| 亚洲,欧美,日韩| 国产亚洲欧美精品永久| 国产日韩欧美亚洲二区| 免费播放大片免费观看视频在线观看| 久久久精品区二区三区| 亚洲精品乱码久久久久久按摩| 能在线免费看毛片的网站| 麻豆成人av视频| 五月玫瑰六月丁香| videos熟女内射| 欧美日韩在线观看h| 在线 av 中文字幕| 精品国产露脸久久av麻豆| 校园人妻丝袜中文字幕| 国产又色又爽无遮挡免| 国产国拍精品亚洲av在线观看| 欧美精品高潮呻吟av久久| av在线观看视频网站免费| 大香蕉97超碰在线| 波野结衣二区三区在线| 黄色毛片三级朝国网站| 亚洲欧美成人精品一区二区| 一级a做视频免费观看| 伊人久久国产一区二区| 又大又黄又爽视频免费| 欧美日韩亚洲高清精品| 国产精品国产三级专区第一集| 国产亚洲精品久久久com| 欧美日本中文国产一区发布| 成人18禁高潮啪啪吃奶动态图 | 午夜av观看不卡| 曰老女人黄片| 成年美女黄网站色视频大全免费 | 国产免费一级a男人的天堂| 18禁在线无遮挡免费观看视频| 欧美激情 高清一区二区三区| 国产男人的电影天堂91| 亚洲av成人精品一区久久| 久久精品国产自在天天线| 日韩一区二区三区影片| 我的老师免费观看完整版| av国产久精品久网站免费入址| 国产毛片在线视频| 久久午夜综合久久蜜桃| 蜜桃在线观看..| 国产男人的电影天堂91| 又黄又爽又刺激的免费视频.| 成人国产麻豆网| 国产无遮挡羞羞视频在线观看| 日韩人妻高清精品专区| 水蜜桃什么品种好| 91久久精品电影网| 日韩亚洲欧美综合| 大陆偷拍与自拍| 啦啦啦视频在线资源免费观看| 国产av一区二区精品久久| 日产精品乱码卡一卡2卡三| 精品少妇内射三级| 黄色配什么色好看| 九色成人免费人妻av| 美女视频免费永久观看网站| 欧美日韩在线观看h| 久久久久国产精品人妻一区二区| 日韩人妻高清精品专区| 狂野欧美激情性bbbbbb| 男女啪啪激烈高潮av片| 国产黄色视频一区二区在线观看| 高清黄色对白视频在线免费看| 丰满少妇做爰视频| 伦理电影大哥的女人| 日韩成人av中文字幕在线观看| 人妻 亚洲 视频| 777米奇影视久久| 亚洲少妇的诱惑av| 美女国产视频在线观看| 精品亚洲乱码少妇综合久久| 日韩制服骚丝袜av| 久久久久久久久久久丰满| 精品久久久久久久久亚洲| 人妻系列 视频| 久久久亚洲精品成人影院| av不卡在线播放| 99视频精品全部免费 在线| 日日爽夜夜爽网站| 女人久久www免费人成看片| 天美传媒精品一区二区| 精品少妇黑人巨大在线播放| xxxhd国产人妻xxx| 国产成人一区二区在线| 国产综合精华液| 国产精品不卡视频一区二区| 久久ye,这里只有精品| 一级毛片aaaaaa免费看小| 亚洲国产欧美在线一区| 国产淫语在线视频| 有码 亚洲区| 大又大粗又爽又黄少妇毛片口| 成人国语在线视频| 亚洲第一av免费看| 2018国产大陆天天弄谢| 精品视频人人做人人爽| 黑人猛操日本美女一级片| 国产色婷婷99| 国产亚洲午夜精品一区二区久久| 秋霞伦理黄片| 狂野欧美激情性xxxx在线观看| 精品人妻在线不人妻| 中文精品一卡2卡3卡4更新| 国产成人精品一,二区| 成人国产av品久久久| 免费观看a级毛片全部| 日本黄色片子视频| 国产午夜精品久久久久久一区二区三区| 人妻系列 视频| 汤姆久久久久久久影院中文字幕| 久久久精品区二区三区| 日本色播在线视频| 熟妇人妻不卡中文字幕| 国产精品秋霞免费鲁丝片| 自线自在国产av| 国产欧美日韩综合在线一区二区| 亚洲精品乱码久久久v下载方式| 久久久久精品久久久久真实原创| 亚洲不卡免费看| 久久精品国产自在天天线| 在线亚洲精品国产二区图片欧美 | av在线老鸭窝| 久久久国产一区二区| 日韩电影二区| 精品一区在线观看国产| 人人妻人人爽人人添夜夜欢视频| 日本爱情动作片www.在线观看| 十八禁网站网址无遮挡| av国产精品久久久久影院| 看免费成人av毛片| 大片免费播放器 马上看| 亚洲人成77777在线视频| 亚洲成人手机| av.在线天堂| 热99国产精品久久久久久7| 亚洲av日韩在线播放| 成年人午夜在线观看视频| 边亲边吃奶的免费视频| 韩国av在线不卡| 免费人妻精品一区二区三区视频| 超碰97精品在线观看| 亚洲天堂av无毛| 女的被弄到高潮叫床怎么办| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品一区蜜桃| 黑人高潮一二区| 精品久久久噜噜| 中文精品一卡2卡3卡4更新| 久久精品熟女亚洲av麻豆精品| 少妇的逼好多水| 考比视频在线观看| 亚洲情色 制服丝袜| 精品人妻熟女毛片av久久网站| 国产黄片视频在线免费观看| 波野结衣二区三区在线| 99久国产av精品国产电影| 我的老师免费观看完整版| 国产 精品1| 亚洲精品视频女| 久久99热6这里只有精品| 国产午夜精品久久久久久一区二区三区| 蜜桃国产av成人99| 亚洲欧洲国产日韩| 国产高清三级在线| 亚洲内射少妇av| 久久久欧美国产精品| 亚洲国产精品一区三区| 久久久久国产网址| 国语对白做爰xxxⅹ性视频网站| 午夜福利在线观看免费完整高清在| 精品一区二区三区视频在线| 超色免费av| 韩国高清视频一区二区三区| 看非洲黑人一级黄片| 99久久人妻综合| 亚洲精品乱久久久久久| 成人午夜精彩视频在线观看| 色视频在线一区二区三区| av有码第一页| 一区二区三区乱码不卡18| 大话2 男鬼变身卡| 一边摸一边做爽爽视频免费| 欧美日韩视频精品一区| 国产乱人偷精品视频| 在线观看三级黄色| 免费看不卡的av| 黄片播放在线免费| 性高湖久久久久久久久免费观看| 日韩免费高清中文字幕av| 亚洲综合精品二区| 美女脱内裤让男人舔精品视频| 热99国产精品久久久久久7| 久久久久久久久久久久大奶| 国产视频内射| 日日撸夜夜添| 久久婷婷青草| 国产精品女同一区二区软件| 色94色欧美一区二区| 高清欧美精品videossex| 美女大奶头黄色视频| 久久99热6这里只有精品| 久久精品国产自在天天线| 精品亚洲成国产av| 欧美另类一区| 人妻制服诱惑在线中文字幕| 夜夜骑夜夜射夜夜干| 91久久精品电影网| 999精品在线视频| 国产伦精品一区二区三区视频9| 亚洲av成人精品一区久久| 免费av中文字幕在线| 人妻一区二区av| 国产成人免费无遮挡视频| 大陆偷拍与自拍| a 毛片基地| 亚洲美女黄色视频免费看| 亚洲美女搞黄在线观看| 亚洲综合色网址| 久久毛片免费看一区二区三区| xxxhd国产人妻xxx| 欧美日韩一区二区视频在线观看视频在线| a 毛片基地| 国产一区二区三区av在线| av线在线观看网站| 日本猛色少妇xxxxx猛交久久| 色5月婷婷丁香| 嫩草影院入口| 国产成人精品无人区| 性色av一级| 成人国产麻豆网| 久久国产精品男人的天堂亚洲 | 久久精品国产亚洲网站| 精品国产乱码久久久久久小说| 午夜福利视频在线观看免费| 国产高清国产精品国产三级| 少妇高潮的动态图| 日本欧美国产在线视频| 免费人妻精品一区二区三区视频| 久久精品国产a三级三级三级| 大片电影免费在线观看免费| 高清视频免费观看一区二区| 国产在视频线精品| 欧美日韩成人在线一区二区| 伊人久久国产一区二区| 丰满乱子伦码专区| 久久久国产精品麻豆| 日韩亚洲欧美综合| 欧美日韩视频精品一区| 又黄又爽又刺激的免费视频.| 久久ye,这里只有精品| 成人免费观看视频高清| a级毛片黄视频| 亚洲色图 男人天堂 中文字幕 | 国产精品国产av在线观看| 黄色毛片三级朝国网站| 纵有疾风起免费观看全集完整版| 亚洲精品中文字幕在线视频| 日韩三级伦理在线观看| 亚洲,欧美,日韩| 亚洲综合精品二区| 亚洲精品,欧美精品| 十八禁高潮呻吟视频| 免费av中文字幕在线| av有码第一页| 中文乱码字字幕精品一区二区三区| 日本与韩国留学比较| 国产欧美日韩一区二区三区在线 | a级片在线免费高清观看视频| 高清不卡的av网站| 两个人的视频大全免费| 欧美日韩亚洲高清精品| 一边亲一边摸免费视频| 国产成人精品福利久久| 精品一区二区三区视频在线| 精品久久久久久久久亚洲| 黄色欧美视频在线观看| 啦啦啦中文免费视频观看日本| 欧美xxⅹ黑人| 飞空精品影院首页| 自拍欧美九色日韩亚洲蝌蚪91| 国产伦理片在线播放av一区| 一本久久精品| 成人手机av| av在线老鸭窝| 黄色怎么调成土黄色| 国产精品久久久久久久电影| 热re99久久精品国产66热6| 日日撸夜夜添| 亚洲无线观看免费| 欧美 亚洲 国产 日韩一| 如日韩欧美国产精品一区二区三区 | 色94色欧美一区二区| 美女cb高潮喷水在线观看| 只有这里有精品99| 国产色爽女视频免费观看| 国产综合精华液| 亚洲第一区二区三区不卡| 999精品在线视频| 亚洲av福利一区| 久久久久国产网址| 精品人妻熟女av久视频| 草草在线视频免费看| 午夜免费男女啪啪视频观看| 一区二区三区乱码不卡18| 日本与韩国留学比较| 欧美精品高潮呻吟av久久| 日本av手机在线免费观看| 永久网站在线| 欧美成人精品欧美一级黄| 九色成人免费人妻av| 女性生殖器流出的白浆| 少妇 在线观看| 日日撸夜夜添| 欧美激情国产日韩精品一区| 国产精品嫩草影院av在线观看| 制服诱惑二区| 国产亚洲av片在线观看秒播厂| 9色porny在线观看| 狠狠婷婷综合久久久久久88av| 热re99久久精品国产66热6| 亚洲国产欧美在线一区| 国产精品一区二区三区四区免费观看| 精品国产国语对白av| 2022亚洲国产成人精品| 久久久国产欧美日韩av| xxxhd国产人妻xxx| 伊人亚洲综合成人网| 飞空精品影院首页| 18禁裸乳无遮挡动漫免费视频| 精品人妻熟女毛片av久久网站| 亚洲国产色片| 亚洲怡红院男人天堂| 欧美激情 高清一区二区三区| 欧美 亚洲 国产 日韩一| 婷婷色综合大香蕉| h视频一区二区三区| 国产免费一区二区三区四区乱码| 日本与韩国留学比较| 中文字幕久久专区| 99久久综合免费| 国产精品国产三级专区第一集| 国产精品麻豆人妻色哟哟久久| 国精品久久久久久国模美| 亚洲精品乱码久久久久久按摩| 久久久国产欧美日韩av| 一区二区三区四区激情视频| 久久免费观看电影| 久久久国产欧美日韩av| 国产av国产精品国产| 亚洲激情五月婷婷啪啪| 91成人精品电影| 综合色丁香网| 国产永久视频网站| 高清在线视频一区二区三区| 综合色丁香网| 久久久亚洲精品成人影院| 日本91视频免费播放| 精品少妇久久久久久888优播| 天天躁夜夜躁狠狠久久av| 午夜福利视频在线观看免费| 国精品久久久久久国模美| 2022亚洲国产成人精品| 黄色欧美视频在线观看| 亚洲婷婷狠狠爱综合网| 制服诱惑二区| 国产午夜精品一二区理论片| 国产在视频线精品|