• <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é)校歌
    果園有了“鵝幫工” 一舉多得效益好
    廣元:治理黑臭水體 再還水清岸美
    久热爱精品视频在线9| 日本精品一区二区三区蜜桃| 国产精品香港三级国产av潘金莲| 免费少妇av软件| 久久人人97超碰香蕉20202| 久久婷婷成人综合色麻豆| 国产精品免费视频内射| 日本黄色视频三级网站网址| 9热在线视频观看99| 国产一区二区激情短视频| 国语自产精品视频在线第100页| 亚洲国产日韩欧美精品在线观看 | 亚洲精品久久国产高清桃花| 真人一进一出gif抽搐免费| 免费看十八禁软件| 99在线视频只有这里精品首页| 国产单亲对白刺激| 在线十欧美十亚洲十日本专区| 丝袜在线中文字幕| 欧美日本中文国产一区发布| 中文字幕人妻熟女乱码| 免费一级毛片在线播放高清视频 | 99久久精品国产亚洲精品| 一级毛片女人18水好多| 国产亚洲av高清不卡| 欧美av亚洲av综合av国产av| 亚洲性夜色夜夜综合| 禁无遮挡网站| 看免费av毛片| 女人被躁到高潮嗷嗷叫费观| 国产欧美日韩精品亚洲av| 嫁个100分男人电影在线观看| 精品久久久久久成人av| 亚洲最大成人中文| 久久精品国产亚洲av香蕉五月| 午夜亚洲福利在线播放| 熟女少妇亚洲综合色aaa.| 欧美乱色亚洲激情| 99re在线观看精品视频| 久久亚洲真实| 午夜福利成人在线免费观看| 欧美国产精品va在线观看不卡| 免费看a级黄色片| 成人免费观看视频高清| 免费不卡黄色视频| 又黄又粗又硬又大视频| 黑人巨大精品欧美一区二区mp4| 精品国产亚洲在线| 色哟哟哟哟哟哟| 亚洲国产欧美日韩在线播放| 女人被躁到高潮嗷嗷叫费观| www.自偷自拍.com| 日韩精品青青久久久久久| 少妇被粗大的猛进出69影院| 成在线人永久免费视频| 国产精品久久视频播放| 1024视频免费在线观看| 免费看十八禁软件| 久久久久九九精品影院| 99在线视频只有这里精品首页| 国产欧美日韩综合在线一区二区| 欧美一级毛片孕妇| xxx96com| 欧美绝顶高潮抽搐喷水| 亚洲熟女毛片儿| 亚洲美女黄片视频| 国产成人精品久久二区二区91| 亚洲成国产人片在线观看| 国产亚洲精品一区二区www| 神马国产精品三级电影在线观看 | 欧美一区二区精品小视频在线| 国产精品一区二区三区四区久久 | x7x7x7水蜜桃| 欧美绝顶高潮抽搐喷水| 国产精品久久久久久精品电影 | 日韩欧美在线二视频| 9热在线视频观看99| 国产高清激情床上av| 亚洲欧美一区二区三区黑人| 亚洲精品国产色婷婷电影| 变态另类成人亚洲欧美熟女 | 九色国产91popny在线| 一边摸一边做爽爽视频免费| 久久精品影院6| 大香蕉久久成人网| 91精品国产国语对白视频| 国产成人欧美| 亚洲av电影在线进入| 国内精品久久久久久久电影| 少妇熟女aⅴ在线视频| 97人妻精品一区二区三区麻豆 | 国产成人影院久久av| 日韩欧美国产在线观看| 女人被躁到高潮嗷嗷叫费观| 亚洲熟妇中文字幕五十中出| 国产亚洲欧美98| 一本大道久久a久久精品| 久久国产精品男人的天堂亚洲| 搡老妇女老女人老熟妇| 99久久综合精品五月天人人| 久久精品亚洲熟妇少妇任你| 日韩精品免费视频一区二区三区| 精品欧美国产一区二区三| 一级毛片女人18水好多| 日本免费一区二区三区高清不卡 | 亚洲国产精品合色在线| 18禁国产床啪视频网站| 亚洲色图 男人天堂 中文字幕| 一级毛片女人18水好多| 国产野战对白在线观看| 69av精品久久久久久| 日韩欧美在线二视频| 女同久久另类99精品国产91| 久久久水蜜桃国产精品网| 人人妻人人澡人人看| 欧美乱码精品一区二区三区| 两个人免费观看高清视频| 免费看十八禁软件| 丝袜美足系列| 一a级毛片在线观看| 大陆偷拍与自拍| 久久久国产成人精品二区| 不卡av一区二区三区| 国产不卡一卡二| 色精品久久人妻99蜜桃| 午夜亚洲福利在线播放| 欧美中文综合在线视频| 99久久久亚洲精品蜜臀av| 琪琪午夜伦伦电影理论片6080| 精品一区二区三区四区五区乱码| 777久久人妻少妇嫩草av网站| 久久天堂一区二区三区四区| 在线观看免费午夜福利视频| 国产成人精品久久二区二区免费| 国产xxxxx性猛交| 一级,二级,三级黄色视频| 欧美久久黑人一区二区| 免费看a级黄色片| 国产成人精品无人区| 日本三级黄在线观看| 日韩中文字幕欧美一区二区| 天堂√8在线中文| 久久性视频一级片| 久久久久精品国产欧美久久久| 色综合欧美亚洲国产小说| 国产三级黄色录像| 好男人在线观看高清免费视频 | 大香蕉久久成人网| 法律面前人人平等表现在哪些方面| 99精品在免费线老司机午夜| 一边摸一边做爽爽视频免费| 国产精品久久久人人做人人爽| 久久精品成人免费网站| 级片在线观看| 日韩精品中文字幕看吧| 国产精品久久久av美女十八| 日韩欧美三级三区| 久久人妻福利社区极品人妻图片| 国产精品久久视频播放| 久久草成人影院| 午夜日韩欧美国产| 成人特级黄色片久久久久久久| 成人三级黄色视频| 91成年电影在线观看| 精品欧美国产一区二区三| 黄色丝袜av网址大全| 超碰成人久久| 久久久久久人人人人人| 91麻豆精品激情在线观看国产| cao死你这个sao货| 97超级碰碰碰精品色视频在线观看| 高潮久久久久久久久久久不卡| 免费无遮挡裸体视频| 少妇 在线观看| 亚洲人成电影免费在线| 首页视频小说图片口味搜索| 别揉我奶头~嗯~啊~动态视频| 欧美日韩亚洲国产一区二区在线观看| 国产精品永久免费网站| 亚洲va日本ⅴa欧美va伊人久久| 国产av在哪里看| 女生性感内裤真人,穿戴方法视频| 国产精品 国内视频| 18美女黄网站色大片免费观看| 操出白浆在线播放| 成熟少妇高潮喷水视频| 黄色丝袜av网址大全| 亚洲av五月六月丁香网| 成人国产一区最新在线观看| 变态另类丝袜制服| 男人的好看免费观看在线视频 | 国产亚洲精品第一综合不卡| 欧美成狂野欧美在线观看| 91在线观看av| 欧美黑人欧美精品刺激| 欧美日韩亚洲国产一区二区在线观看| 国产精品久久久久久亚洲av鲁大| 久久久久九九精品影院| 搡老熟女国产l中国老女人| 少妇被粗大的猛进出69影院| 成年女人毛片免费观看观看9| 中文亚洲av片在线观看爽| 美女高潮到喷水免费观看| 色av中文字幕| 少妇粗大呻吟视频| 亚洲国产精品成人综合色| 亚洲欧美日韩高清在线视频| avwww免费| 一区在线观看完整版| 国产精品乱码一区二三区的特点 | 国产熟女xx| 婷婷六月久久综合丁香| 操出白浆在线播放| 国产精品亚洲一级av第二区| 俄罗斯特黄特色一大片| 精品第一国产精品| 午夜影院日韩av| 给我免费播放毛片高清在线观看| 亚洲一区中文字幕在线| 嫩草影院精品99| av天堂久久9| av天堂在线播放| 丝袜人妻中文字幕| 欧美不卡视频在线免费观看 | 亚洲激情在线av| e午夜精品久久久久久久| 中文亚洲av片在线观看爽| 日韩欧美免费精品| 国产精品电影一区二区三区| 亚洲国产看品久久| 国产亚洲欧美98| 亚洲中文字幕一区二区三区有码在线看 | videosex国产| 久久久久久久久免费视频了| 国产主播在线观看一区二区| 一进一出抽搐动态| 亚洲人成网站在线播放欧美日韩| 久久中文字幕人妻熟女| 91精品国产国语对白视频| av视频免费观看在线观看| 老司机午夜十八禁免费视频| 精品久久久久久久人妻蜜臀av | 丝袜人妻中文字幕| 搞女人的毛片| 久久香蕉国产精品| 精品一区二区三区av网在线观看| 一区二区日韩欧美中文字幕| 黄色片一级片一级黄色片| 国产精品1区2区在线观看.| 视频区欧美日本亚洲| 在线播放国产精品三级| av在线播放免费不卡| 18禁黄网站禁片午夜丰满| 久久国产精品影院| 国产成人精品无人区| 麻豆一二三区av精品| 在线永久观看黄色视频| 久久久国产成人精品二区| 每晚都被弄得嗷嗷叫到高潮| 国产精华一区二区三区| 国产精品98久久久久久宅男小说| 免费高清在线观看日韩| 国产激情欧美一区二区| 91麻豆精品激情在线观看国产| 成人国语在线视频| 97碰自拍视频| 欧美人与性动交α欧美精品济南到| 亚洲av五月六月丁香网| 日本a在线网址| 久久人人爽av亚洲精品天堂| 亚洲全国av大片| 亚洲性夜色夜夜综合| 国产麻豆成人av免费视频| 欧美日本视频| 久久国产精品男人的天堂亚洲| 久久久久久大精品| 非洲黑人性xxxx精品又粗又长| 欧美黄色淫秽网站| 不卡av一区二区三区| 久久国产乱子伦精品免费另类| 亚洲精华国产精华精| 久久影院123| 午夜精品国产一区二区电影| 天堂动漫精品| 亚洲精品久久国产高清桃花| 国产成人啪精品午夜网站| 女性被躁到高潮视频| 97人妻精品一区二区三区麻豆 | 日本 av在线| 一区二区三区精品91| 久久 成人 亚洲| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲aⅴ乱码一区二区在线播放 | 变态另类丝袜制服| 欧美色欧美亚洲另类二区 | 亚洲精品中文字幕一二三四区| √禁漫天堂资源中文www| 日本黄色视频三级网站网址| 久久久久久久午夜电影| 首页视频小说图片口味搜索| 在线观看66精品国产| 免费看美女性在线毛片视频| 国产亚洲精品av在线| 国产欧美日韩精品亚洲av| 日韩欧美国产在线观看| 国产精品香港三级国产av潘金莲| e午夜精品久久久久久久| 欧美一区二区精品小视频在线| 黄片大片在线免费观看| 亚洲欧美日韩另类电影网站| 1024视频免费在线观看| 在线观看66精品国产| 国产成人精品久久二区二区91| 搡老熟女国产l中国老女人| av有码第一页| 丁香六月欧美| 国产午夜精品久久久久久| 最新美女视频免费是黄的| 久久精品aⅴ一区二区三区四区| 久久久国产成人精品二区| 黑人巨大精品欧美一区二区蜜桃| 制服诱惑二区| 亚洲成人久久性| 嫁个100分男人电影在线观看| 成人国产综合亚洲| 亚洲aⅴ乱码一区二区在线播放 | 99国产精品99久久久久| 国产精品日韩av在线免费观看 | 女生性感内裤真人,穿戴方法视频| 两个人看的免费小视频| 欧美老熟妇乱子伦牲交| 最近最新中文字幕大全电影3 | 久久精品影院6| 久久人妻av系列| 在线永久观看黄色视频| 中文字幕最新亚洲高清| 国产精品免费一区二区三区在线| 午夜久久久久精精品| 日韩欧美一区视频在线观看| 国产激情久久老熟女| 在线视频色国产色| 欧美中文日本在线观看视频| 亚洲国产精品sss在线观看| 在线观看日韩欧美| 亚洲 国产 在线| 中文字幕人妻熟女乱码| 亚洲专区中文字幕在线| 国产成人一区二区三区免费视频网站| 国产精品影院久久| 精品无人区乱码1区二区| ponron亚洲| 欧美色视频一区免费| 亚洲熟女毛片儿| 丝袜人妻中文字幕| 久久精品91蜜桃| 大香蕉久久成人网| 999精品在线视频| 国产色视频综合| 亚洲五月色婷婷综合| 亚洲五月天丁香| 美女国产高潮福利片在线看| 国产精品 国内视频| 老司机在亚洲福利影院| 18美女黄网站色大片免费观看| 午夜激情av网站| 亚洲aⅴ乱码一区二区在线播放 | 免费在线观看视频国产中文字幕亚洲| 俄罗斯特黄特色一大片| 欧美绝顶高潮抽搐喷水| 亚洲一码二码三码区别大吗| 亚洲av五月六月丁香网| 久久精品成人免费网站| 精品国产超薄肉色丝袜足j| a在线观看视频网站| 在线观看66精品国产| 麻豆av在线久日| 国产精品免费视频内射| 在线观看免费日韩欧美大片| 亚洲国产欧美日韩在线播放| 热99re8久久精品国产| 欧美+亚洲+日韩+国产| 无人区码免费观看不卡| 免费看a级黄色片| 亚洲免费av在线视频| 午夜福利18| 少妇的丰满在线观看| 两个人免费观看高清视频| 欧美国产日韩亚洲一区| 色综合亚洲欧美另类图片| 色尼玛亚洲综合影院| 国产精品久久久久久精品电影 | 亚洲中文字幕日韩| 后天国语完整版免费观看| 国产一区二区三区视频了| 99在线视频只有这里精品首页| 少妇粗大呻吟视频| x7x7x7水蜜桃| 美女大奶头视频| 精品日产1卡2卡| 国产真人三级小视频在线观看| 人妻久久中文字幕网| 国产精品香港三级国产av潘金莲| 亚洲欧洲精品一区二区精品久久久| av视频免费观看在线观看| 国产高清videossex| 国产精品久久久久久精品电影 | 久久久久久免费高清国产稀缺| 大香蕉久久成人网| 老熟妇仑乱视频hdxx| 99精品久久久久人妻精品| 一进一出抽搐gif免费好疼| 中文字幕人妻熟女乱码| 日本五十路高清| 女人被狂操c到高潮| 无遮挡黄片免费观看| 欧美一级毛片孕妇| av视频免费观看在线观看| 久久午夜综合久久蜜桃| 19禁男女啪啪无遮挡网站| 国产成人精品久久二区二区免费| 免费女性裸体啪啪无遮挡网站| 岛国在线观看网站| 99riav亚洲国产免费| 在线观看www视频免费| 好男人在线观看高清免费视频 | 欧美另类亚洲清纯唯美| 首页视频小说图片口味搜索| 亚洲色图综合在线观看| 国产一区二区三区视频了| 亚洲av美国av| 久久性视频一级片| 欧美一级a爱片免费观看看 | АⅤ资源中文在线天堂| 可以在线观看的亚洲视频| 亚洲国产欧美网| 欧美在线黄色| 一级,二级,三级黄色视频| 欧美久久黑人一区二区| 91字幕亚洲| 免费看美女性在线毛片视频| 法律面前人人平等表现在哪些方面| 黄网站色视频无遮挡免费观看| 免费高清视频大片| 美女 人体艺术 gogo| 国产精品国产高清国产av| 叶爱在线成人免费视频播放| 成人精品一区二区免费| 亚洲精品美女久久久久99蜜臀| 久久久久精品国产欧美久久久| 亚洲国产看品久久| 精品第一国产精品| 视频区欧美日本亚洲| 欧美绝顶高潮抽搐喷水| 亚洲精品久久国产高清桃花| 精品电影一区二区在线| e午夜精品久久久久久久| 一区福利在线观看| 黄片播放在线免费| 成在线人永久免费视频| 欧美精品亚洲一区二区| 国产欧美日韩一区二区精品| 日本五十路高清| 欧美另类亚洲清纯唯美| 91精品三级在线观看| 亚洲av成人一区二区三| 国产精品久久久久久人妻精品电影| 亚洲三区欧美一区| 国产麻豆成人av免费视频| 美女大奶头视频| 正在播放国产对白刺激| 亚洲av日韩精品久久久久久密| 熟妇人妻久久中文字幕3abv| 国产成人影院久久av| 国产精品综合久久久久久久免费 | 午夜福利高清视频| 久久精品国产99精品国产亚洲性色 | 黄频高清免费视频| 一边摸一边做爽爽视频免费| 午夜老司机福利片| 亚洲精品久久成人aⅴ小说| 亚洲在线自拍视频| 亚洲九九香蕉| 老司机在亚洲福利影院| 欧美日韩一级在线毛片| 欧美在线黄色| 怎么达到女性高潮| 免费看美女性在线毛片视频| 国产亚洲av高清不卡| 好男人在线观看高清免费视频 | 男女之事视频高清在线观看| 亚洲欧美日韩另类电影网站| 三级毛片av免费| 在线观看一区二区三区| 一级毛片女人18水好多| 亚洲欧美一区二区三区黑人| 国产一级毛片七仙女欲春2 | 欧美精品啪啪一区二区三区| 亚洲成av人片免费观看| 久久久国产欧美日韩av| 精品第一国产精品| 久久中文字幕一级| 麻豆成人av在线观看| 99香蕉大伊视频| 亚洲国产毛片av蜜桃av| 欧美久久黑人一区二区| 国内毛片毛片毛片毛片毛片| 精品不卡国产一区二区三区| 国产精品秋霞免费鲁丝片| 在线观看免费视频网站a站| 淫妇啪啪啪对白视频| 久久中文字幕一级| 亚洲国产高清在线一区二区三 | 亚洲第一电影网av| 久久久久久久午夜电影| 国语自产精品视频在线第100页| 午夜精品久久久久久毛片777| 天天躁狠狠躁夜夜躁狠狠躁| 国产在线精品亚洲第一网站| 亚洲九九香蕉| 激情视频va一区二区三区| 夜夜夜夜夜久久久久| 国产一卡二卡三卡精品| 91在线观看av| 午夜免费激情av| 神马国产精品三级电影在线观看 | 可以免费在线观看a视频的电影网站| 欧美av亚洲av综合av国产av| 成人亚洲精品av一区二区| 日本 av在线| 国产人伦9x9x在线观看| 国产精品一区二区三区四区久久 | 亚洲avbb在线观看| av视频免费观看在线观看| 91大片在线观看| 91精品三级在线观看| 日日爽夜夜爽网站| 天天添夜夜摸| 久久精品亚洲精品国产色婷小说| 狠狠狠狠99中文字幕| 激情在线观看视频在线高清| 午夜福利在线观看吧| 久久香蕉精品热| 午夜福利视频1000在线观看 | 欧美成狂野欧美在线观看| 悠悠久久av| 日韩有码中文字幕| 亚洲精品在线观看二区| 久久草成人影院| 亚洲片人在线观看| 亚洲人成电影免费在线| 久久欧美精品欧美久久欧美| 夜夜夜夜夜久久久久| 搡老岳熟女国产| 亚洲国产中文字幕在线视频| 久久中文看片网| 脱女人内裤的视频| 亚洲人成77777在线视频| 中文字幕精品免费在线观看视频| 亚洲一码二码三码区别大吗| 黑人巨大精品欧美一区二区mp4| 给我免费播放毛片高清在线观看| 午夜精品国产一区二区电影| 99久久精品国产亚洲精品| 美女免费视频网站| 久久伊人香网站| 黄片小视频在线播放| 欧美最黄视频在线播放免费| 天天添夜夜摸| 一边摸一边抽搐一进一出视频| 在线观看午夜福利视频| 国产一级毛片七仙女欲春2 | 日本撒尿小便嘘嘘汇集6| 免费av毛片视频| 精品卡一卡二卡四卡免费| 18禁裸乳无遮挡免费网站照片 | 1024香蕉在线观看| 一本综合久久免费| 如日韩欧美国产精品一区二区三区| 一级a爱片免费观看的视频| 可以在线观看毛片的网站| 国产精品久久久久久人妻精品电影| 伊人久久大香线蕉亚洲五| 久久久久精品国产欧美久久久| 丝袜美腿诱惑在线| 国产成人欧美| 午夜成年电影在线免费观看| 国产精品自产拍在线观看55亚洲| 午夜两性在线视频| 99久久综合精品五月天人人| 国产亚洲欧美在线一区二区| 亚洲欧美日韩另类电影网站| 亚洲成国产人片在线观看| 啦啦啦 在线观看视频| 久久久久国产精品人妻aⅴ院| 18禁黄网站禁片午夜丰满| 亚洲激情在线av| 日韩国内少妇激情av| 可以在线观看的亚洲视频| 亚洲激情在线av| 99国产精品一区二区三区| 操出白浆在线播放| 国产私拍福利视频在线观看| 欧美性长视频在线观看| 在线观看午夜福利视频| 好男人电影高清在线观看| 国产亚洲欧美精品永久| 欧美激情 高清一区二区三区| 级片在线观看| 国产精品九九99| 色婷婷久久久亚洲欧美| 国产成人影院久久av| 黄网站色视频无遮挡免费观看| 国产欧美日韩一区二区精品| 精品人妻1区二区| 久久精品国产99精品国产亚洲性色 |