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

    感潮河網(wǎng)降雨徑流污染空間分析與模擬

    2021-04-30 03:47:16張鳳山尚明珠趙朋曉程開宇唐穎棟
    中國環(huán)境科學 2021年4期
    關(guān)鍵詞:茅洲排口河網(wǎng)

    張鳳山,尚明珠,趙朋曉,程開宇,唐穎棟,魏 俊

    感潮河網(wǎng)降雨徑流污染空間分析與模擬

    張鳳山,尚明珠,趙朋曉,程開宇,唐穎棟,魏 俊*

    (中國電建集團華東勘測設(shè)計研究院有限公司,浙江 杭州 311122)

    為了探究茅洲河流域感潮河網(wǎng)面源污染空間分布特征和降雨徑流污染規(guī)律,基于空間分析、統(tǒng)計分析與流域水動力-水質(zhì)耦合模擬方法,對典型降雨情景下河網(wǎng)水質(zhì)情況進行模擬分析,提出基于水質(zhì)改善目標的生態(tài)補水點位空間布局優(yōu)化策略.研究表明,層次聚類凝聚算法和K-均值法迭代組合可以較好地實現(xiàn)面源污染分級與分類;茅洲河各支流中,石巖渠、松崗河中上游等河道(段)由于面源污染負荷相對較高且缺乏生態(tài)補水,雨后水質(zhì)恢復(fù)緩慢;基于補水總量不變原則,對生態(tài)補水方案進行局部優(yōu)化,優(yōu)化結(jié)果可使雨后受污染重點河道(段)水質(zhì)恢復(fù)速度加快一倍以上,提高了流域水質(zhì)的整體穩(wěn)定性.研究結(jié)論可為進一步認識茅洲河流域水污染特征、實現(xiàn)流域水環(huán)境精細化管理提供支撐.

    感潮河網(wǎng);面源污染;數(shù)值模擬;茅洲河

    降雨導(dǎo)致的面源污染由于其不確定性高,時空變異性大,是影響城市水環(huán)境的主要因素之一[1-4],目前的研究涵蓋流域[5-6]、城市[7-9]、小區(qū)或?qū)W校[10-11]等不同尺度,重點關(guān)注面源污染對下游受納水體水質(zhì)的沖擊[6,8-9].車蕊等[6]對東江流域近38a極端降雨事件分析發(fā)現(xiàn),氨氮在降雨過程中呈前期高,中后期低的特征.許淑敏[9]對引水工程緩解海河流域面源污染程度進行分析,發(fā)現(xiàn)雨量較小時調(diào)水引流效果更佳.對于感潮河網(wǎng)地區(qū),污染物遷移轉(zhuǎn)化受到徑流和潮流共同作用,不利于擴散和降解[12],更容易受到降雨徑流污染的沖擊.近年來,珠江三角洲各地開展了一系列感潮河網(wǎng)區(qū)水動力水質(zhì)模擬研究,論證了閘泵優(yōu)化調(diào)度[13]、引水增流[14-15]等措施對河網(wǎng)水質(zhì)恢復(fù)的效果.研究[15]表明,珠江三角洲地區(qū)受到污染-咸潮雙重影響,利用閘泵-河庫聯(lián)合調(diào)度是改善珠江三角洲感潮河網(wǎng)水環(huán)境狀況的有效措施之一.

    茅洲河流域?qū)僦榻侵薷谐焙泳W(wǎng),地處深圳市西部,于伶仃洋交椅灣入海,流域內(nèi)為高密度建成區(qū),老城區(qū)和老工業(yè)區(qū)眾多,產(chǎn)污量大,河流徑污比嚴重偏低,面源污染負荷高,面源污染中COD和NH3-N分別占總排放污染負荷的19%和11.2%[16],其中老城區(qū)特別是城中村面源污染風險最為突出.由于徑流是驅(qū)動污染物向下游擴散的主要動力[17],而伶仃洋為弱潮河口,潮差較小,潮動力不足[18-19],不利于污染物遷移擴散.同時,受氣候變化和城市化發(fā)展等因素影響,深圳西部城區(qū)近40a間年降水量、汛期降水量和極端降水指標均呈現(xiàn)增加趨勢[20],下墊面改變和降雨特征變化影響了茅洲河流域面源污染的時空分布規(guī)律.特定的自然環(huán)境條件和社會發(fā)展背景導(dǎo)致茅洲河流域面源污染風險居高不下.茅洲河經(jīng)過系統(tǒng)治理,截止2019年底,流域水質(zhì)明顯改善,旱季水質(zhì)較穩(wěn)定.但由于茅洲河較高的面源污染負荷和不利的河口潮動力條件,雨后污染風險仍然較高.

    本研究以茅洲河感潮河網(wǎng)為研究區(qū)域,基于GIS空間分析和聚類分析理論,對沿河雨水排口污染等級進行分級分類,提高了大量排口數(shù)據(jù)統(tǒng)計分析的效率;在此基礎(chǔ)上通過流域水動力-水質(zhì)耦合模型開展典型降雨情景面源污染規(guī)律研究,提出改善雨后河道水質(zhì)恢復(fù)規(guī)律的工程措施,為揭示茅洲河流域水污染特征、保障河網(wǎng)水質(zhì)穩(wěn)定提供參考.

    1 研究區(qū)域與方法

    1.1 研究區(qū)域概況

    圖1 茅洲河流域示意

    茅洲河是深圳第一大河,發(fā)源于羊臺山北麓,地跨深圳、東莞兩市,在沙井民主村匯入伶仃洋.流域面積344.23km2,干流全長30.69km.其中,寶安區(qū)境內(nèi)流域面積122.65km2,干流河長19.71km,感潮河段長約13km,下游河口段11.4km為深圳市與東莞市界河.茅洲河流域?qū)毎矃^(qū)境內(nèi)共有干、支流19條,河道總長度96.56km.目前,茅洲河流域(寶安片區(qū))已經(jīng)形成以松崗污水廠和沙井污水廠再生水為主要水源的生態(tài)補水系統(tǒng),如圖1所示,補水規(guī)模80萬m3/d,實測再生水水質(zhì)主要指標優(yōu)于地表水IV類,是改善茅洲河干支流水質(zhì)的重要手段之一.

    1.2 研究方法

    1.2.1 聚類分析法 聚類分析廣泛應(yīng)用于水環(huán)境數(shù)據(jù)分析和水污染綜合評價[21-23].層次聚類方法是一種常用的聚類分析算法,可分為凝聚和分裂兩種方法[24].K-均值算法是聚類算法中最基礎(chǔ)也最重要的無監(jiān)督聚類算法,使類內(nèi)具有較高的相似度,而類間的相似度較低,適用于數(shù)值型數(shù)據(jù)且易于實現(xiàn),時間復(fù)雜度低,算法的可解釋度較強[25].

    本文運用SPSS分析工具進行雨水排口聚類分析,首先應(yīng)用層次聚類凝聚算法得到結(jié)果類的數(shù)目,在此基礎(chǔ)上應(yīng)用K-均值法改進聚類結(jié)果.

    1.2.2 面源污染計算 降雨徑流污染強度由污染物累積過程和沖刷過程共同決定,采用飽和函數(shù)式(1)和指數(shù)函數(shù)式(2)分別計算污染物累積和沖刷過程[26].

    式中:1為最大增長可能, kg/hm2;2為半飽和常數(shù), (達到最大增長一半時的天數(shù)), d.

    式中:為污染物沖刷量, kg/h;1為沖刷系數(shù);2為沖刷指數(shù);為單位面積的徑流速率, mm/h;為污染物增長質(zhì)量, kg/hm2.

    本文基于典型降雨情景及排口聚類結(jié)果,采用雨水管理模型(SWMM)模型模擬面源污染情況.

    1.2.3 河網(wǎng)水動力水質(zhì)模擬 水動力模塊控制方程為圣維南方程[27](式3、式4):

    式中:為河道過水面積, m2;為流量, m3/s;為側(cè)向來流在河道方向的流速, m/s;為時間, s;為沿水流方向的水平坐標, m;為河道的側(cè)向來流量, m3/s;為動量修正系數(shù);為重力加速度,m/s2;為水位, m;S為摩阻坡降.

    水質(zhì)模塊控制方程為對流擴散方程[27](式5):

    式中:為模擬物質(zhì)的濃度, mg/L;為河流斷面平均流速, m/s;E為對流擴散系數(shù), m2/s;為模擬物質(zhì)的一級衰減系數(shù), mg/(L·s);為空間坐標, m;為時間, s.

    本文采用MIKE11模型對茅洲河流域河網(wǎng)水動力、水質(zhì)過程進行模擬.

    2 結(jié)果與討論

    2.1 降雨特征及典型降雨分析

    近40a來茅洲河流域所在的深圳市西部城區(qū)年降水量整體呈現(xiàn)增加趨勢[20],表現(xiàn)出一定的年際變化規(guī)律.為了研究近期茅洲河流域降雨特征,本文基于近5a茅洲河流域日降雨過程進行統(tǒng)計分析,結(jié)果表明近5a流域年平均降雨日數(shù)102d,主要集中在4~9月份,日降雨以小雨(小于10.0mm)為主,占總降雨日數(shù)的64.7%,中雨(10.0~24.9mm)占20.3%,大雨及以上(大于25.0mm)占15%.統(tǒng)計表明超過76%的降雨間隔時間在3d以內(nèi),僅約10%的降雨時間間隔在7d以上,從圖2各量級降雨事件的降雨間隔分布圖可看出,茅洲河流域降雨時間間隔普遍較短,連續(xù)性降雨頻發(fā).

    圖2 茅洲河流域降雨特征統(tǒng)計

    圖3 茅洲河流域典型降雨過程

    在降雨量相當時,雨峰偏后的降雨污染負荷大于雨峰偏前情況[28];當其他條件一定時,降雨量累積值、最大降雨強度和平均降雨強度越大,污染負荷越大[29].結(jié)合前人研究經(jīng)驗與茅洲河流域降雨規(guī)律分析結(jié)果,按照同一量級降雨事件中降雨量較大、前期干燥天數(shù)較長、雨型集中、雨峰偏后的原則選取共3場典型降雨(小雨、中雨、大雨各1場),作為河網(wǎng)水環(huán)境模擬的典型降雨事件,日降雨量分別為8.8mm(小雨)、17.6mm(中雨)以及43.4mm(大雨),3場典型降雨逐小時降雨過程如圖3所示.

    2.2 面源污染特征空間分析

    面源污染負荷強度與下墊面密切相關(guān).依據(jù)《深圳市面源污染整治管控技術(shù)路線及技術(shù)指南(試行)》[30]中對深圳市下墊面分類的方法,按面源污染負荷從低到高將茅洲河寶安片區(qū)用地分為A、B、C、D 4類,各等級代表性下墊面見表1.

    表1 面源污染等級標準

    按以上地塊劃分原則對茅洲河流域?qū)毎财瑓^(qū)下墊面面源污染等級進行空間分析(圖4),由于茅洲河流域老城區(qū)和老工業(yè)區(qū)眾多,C類和D類等高污染負荷用地是該區(qū)域主要的用地類型,而A類和B類的用地明顯較少,反映了茅洲河寶安片區(qū)段干支流沿岸開發(fā)密度高、面源污染風險大的特點.

    圖4 下墊面類型

    經(jīng)統(tǒng)計茅洲河寶安片區(qū)沿河共分布428個雨水排口,各排口服務(wù)范圍面積與用地性質(zhì)差異顯著.通過GIS空間分析,以各沿河雨水排口服務(wù)范圍為統(tǒng)計單元,按排口污染物濃度由低到高將面源污染等級分為I~V級,如圖5所示.當降雨事件一定時,降雨徑流污染過程主要與排口服務(wù)面積及下墊面類型密切相關(guān),特異性明顯;而排口服務(wù)范圍內(nèi)用地類型相似且匯流特性接近的區(qū)域污染物累積與沖刷過程又呈現(xiàn)相似特征.因此,為了統(tǒng)計分析茅洲河流域雨水排口特性,同時提高建模效率,有必要對雨水排口進行分類統(tǒng)計.本文以排口服務(wù)面積與污染負荷等級為聚類因子,采用先驗策略,應(yīng)用層次聚類凝聚算法和K-均值法進行迭代優(yōu)化,最終將雨水排口分為18個類別,排口分類結(jié)果見表2.

    圖5 排口污染等級

    表2 排口分類結(jié)果

    圖6 排口聚類效果檢驗

    圖7 不同污染負荷等級代表性排口分布示意

    為驗證聚類效果,如圖6所示,排口呈現(xiàn)出組內(nèi)相似、組間差異明顯的特點,表明分類效果良好.其中,第2、12類排口污染負荷高且集水面積較大,負荷總量高.同時,聚類分析結(jié)果可以表征不同類型排口之間的污染等級及其空間分布,如圖7所示.其中,第3類排口為高污染排口,主要分布在沙井河及其支流、排澇河沿岸;第4類排口為低污染排口,主要分布在老虎坑、龜嶺東沿岸;第11類排口污染負荷中等,在各支流沿線廣泛分布.

    2.3 茅洲河流域水環(huán)境模型建立與驗證

    模型范圍為茅洲河流域?qū)毎财瑓^(qū)干、支流,共概化河道(段)25條、節(jié)點452個、河道斷面322個、水閘14座.模型選取2019年3月21日~3月22日26h全潮期為率定期,2019年11月為驗證期,以共和村斷面為參證斷面,通過水動力、水質(zhì)同步監(jiān)測數(shù)據(jù)進行模型參數(shù)率定和驗證,采用納什效率系數(shù)(NSE)對模擬精度進行評價.結(jié)合相關(guān)研究成果[31]與流域資料情況,河道糙率取值為0.028~0.032,擴散系數(shù)取10m2/s.模型率定和驗證結(jié)果如圖8所示.其中圖8(a)、圖8(b)為率定期共和村斷面流量、水位模擬精度檢驗,NSE分別達到0.993、0.988,表明模型具有可靠的水動力模擬精度;圖8(c)、圖8(d)分別為率定期和驗證期共和村斷面氨氮濃度模擬精度檢驗,NSE分別為0.715和0.841,水質(zhì)模擬精度較高.率定和驗證結(jié)果表明,茅洲河流域水環(huán)境模型可以較好地反應(yīng)流域水動力和水污染特征,模擬結(jié)果較為可靠.

    圖8 模型率定和驗證

    2.4 雨后河道水質(zhì)變化規(guī)律及生態(tài)補水優(yōu)化分析

    目前茅洲河流域已建成較為完善的補水系統(tǒng),但尚未建立針對雨后河道水質(zhì)改善的補水調(diào)度策略[32].在現(xiàn)狀工程條件下,基于已建立的茅洲河流域模型,在2.1節(jié)選定的典型降雨情景下分析茅洲河干流、支流雨后水質(zhì)變化規(guī)律.茅洲河流域降雨期間和雨后水質(zhì)情況如圖9所示,由于污染風險和工程背景的差異,降雨對干支流各河道水質(zhì)呈現(xiàn)出不同影響,其中,潭頭河、沙井河、潭頭渠由于沿線高污染風險排口較為密集,降雨期間水質(zhì)惡化明顯,但由于現(xiàn)狀補水系統(tǒng)完善,水質(zhì)恢復(fù)較快;石巖渠、松崗河中上游、七支渠上游、萬豐河上游等河道(段)由于污染負荷相對較高且現(xiàn)狀缺少生態(tài)補水,水環(huán)境容量不足,雨后水質(zhì)恢復(fù)緩慢.

    生態(tài)補水是改善河道水質(zhì)的有效措施,受到補水水質(zhì)、補水量、補水位置等因素的影響,不同補水調(diào)度方式的效果差異顯著[33].為了加快雨后河道水質(zhì)改善過程,基于補水總量不變原則,對現(xiàn)狀生態(tài)補水方案進行局部優(yōu)化,調(diào)整部分河道補水量和補水位置,優(yōu)化策略見表3.以中雨為例,補水方案優(yōu)化后七支渠、萬豐河、松崗河、石巖渠雨后水質(zhì)恢復(fù)至V類水的速度加快一倍以上;沙井河、上寮河由于現(xiàn)狀補水量較大,適當縮減補水量對其水質(zhì)恢復(fù)影響較小.雨后流域水質(zhì)情況如圖10所示,與圖9結(jié)果相比,在保持補水總量不變的前提下,優(yōu)化方案可顯著提高流域水質(zhì)的整體穩(wěn)定性.

    表3 補水方案優(yōu)化策略

    圖10 補水方案優(yōu)化后茅洲河干支流雨后水質(zhì)情況(NH3-N濃度)

    3 結(jié)論

    3.1 茅洲河流域?qū)毎财瑓^(qū)開發(fā)密度高,下墊面類型以高污染負荷地塊為主,面源污染風險大.基于GIS空間分析和聚類分析理論,以研究區(qū)域內(nèi)沿河污染排口服務(wù)范圍為統(tǒng)計單元,排口服務(wù)范圍和污染等級為聚類因子,采用先驗策略,應(yīng)用層次聚類算法和K-均值算法進行迭代優(yōu)化,將沿河雨水排口分為18類,經(jīng)檢驗聚類效果良好.

    3.2 流域水環(huán)境模型選取26h全潮期為率定期, 30d為驗證期,以共和村斷面為參證站,采用同步水動力、水質(zhì)實測數(shù)據(jù)進行參數(shù)率定和驗證.率定期和驗證期的流量、水位模擬結(jié)果NSE均達到0.99左右;率定期和驗證期水質(zhì)模擬結(jié)果NSE在0.7以上,模型參數(shù)較為可靠.

    3.3 在現(xiàn)狀工程條件下分析茅洲河流域降雨期間和雨后水質(zhì)變化情況,干、支流各河道水質(zhì)變化呈現(xiàn)不同規(guī)律,部分河段沿線污染風險較高,雨后水質(zhì)明顯惡化,但現(xiàn)狀較為完善的生態(tài)補水系統(tǒng)可使其水質(zhì)較快恢復(fù);但仍部分游河道(段)由于污染負荷相對較高且生態(tài)補水不足,雨后水質(zhì)不易恢復(fù).基于生態(tài)補水總量不變的原則優(yōu)化現(xiàn)有補水策略,結(jié)果表明優(yōu)化后重點污染河道(段)雨后水質(zhì)恢復(fù)速率明顯加快,同時對其他河道無明顯不利影響,流域水質(zhì)整體穩(wěn)定性提高.

    [1] 丁程程,劉 健.中國城市面源污染現(xiàn)狀及其影響因素 [J]. 中國人口.資源與環(huán)境, 2011,127(S1):86-89. Ding C C, Liu J. Discussion on urban non-point source pollution and control technologies in China [J]. China population, resources and environment, 2011,127(S1):86-89.

    [2] 耿潤哲,梁璇靜,殷培紅,等.面源污染最佳管理措施多目標協(xié)同優(yōu)化配置研究進展 [J]. 生態(tài)學報, 2019,39(8):2667-2675. Geng R Z, Liang X J, Yin P H, et al. A review: multi-objective collaborative optimization of best management practices for non-point sources pollution contro1 [J]. Aeta Ecologica Sinica, 2019,39(8): 2667-2675.

    [3] Charters F J, Cochrane T A, O’Sullivan A D. Predicting event-based sediment and heavy metal loads in untreated urban runoff from impermeable surfaces [J]. Water, 2020,12:969.

    [4] 朱聞博,王 健,薛 菲,等.從海綿城市到多維海綿——系統(tǒng)解決城市水問題 [M]. 南京:江蘇鳳凰科學出版社, 2018:37. Zhu W B, Wang J, Xue F, et al. Integrated solutions for urban water issues [M]. Nanjing: Phoenix Science Press, 2018:37.

    [5] 彭兆弟,李勝生,劉 莊,等.太湖流域跨界區(qū)農(nóng)業(yè)面源污染特征 [J]. 生態(tài)與農(nóng)村環(huán)境學報, 2016,32(3):458-465. Peng Z D, Li S S, Liu Z, et al. Characteristics of transboundary non-point source agricultural pollution in the Taihu valley [J]. Journal of Ecology and Rural Environment, 2016,32(3):458-465.

    [6] 車 蕊,林 澍,范中亞,等.連續(xù)極端降雨對東江流域水質(zhì)影響分析 [J]. 環(huán)境科學, 2019,40(10):140-149. Che R, Lin S, Fan Z Y, et al. Effects of continuous extreme rainfall on water quality of the Dongjiang river basin [J]. Environmental Science, 2019,40(10):140-149.

    [7] 何夢男,張 勁,陳 誠,等.上海市淀北片降雨徑流過程污染時空特性分析 [J]. 環(huán)境科學學報, 2018,38(2):103-112. He M N, Zhang J, Chen C, et al. Analysis of the temporal and spatial characteristics of rainfall-runoff pollution in Dianbei basin of Shanghai [J]. Acta Scientiae Circumstantiae, 2018,38(2):103-112.

    [8] Li C, Zheng X, Zhao F, et al. Effects of urban non-point source pollution from Baoding city on Baiyangdian lake [J]. Water, 2017,9(4): 249-266.

    [9] 許淑敏.天津市降雨徑流對海河干流水質(zhì)的影響研究 [D]. 天津:天津大學, 2018. Xu S M. Study on the effects of rainfall-runoff on the water quality in the Haihe mainstream in Tianjin [D]. Tianjin: Tianjin University, 2018.

    [10] Gong Y W, Liang X, Li X, et al. Influence of rainfall characteristics on total suspended solids in urban runoff: a case study in Beijing [J]. Water, 2016,8:278-300.

    [11] Qiu S, Yin H, Deng J, et al. Cost-effectiveness analysis of green–gray stormwater control measures for non-point source pollution [J]. International Journal of Environmental Research and Public Health, 2020,17:998-1010.

    [12] 杜 建,陳曉宏,陳志和,等.珠江三角洲感潮河網(wǎng)區(qū)水環(huán)境引水調(diào)控研究 [J]. 水文, 2012,32(4):18-23,98. Du J, Chen X H, Chen Z H, et al. Water diversion and control for water environment in tidal network of Pearl River Delta [J]. Journal of China Hydrology, 2012,32(4):18-23,98.

    [13] 江 濤,朱淑蘭,張 強,等.潮汐河網(wǎng)閘泵聯(lián)合調(diào)度的水環(huán)境效應(yīng)數(shù)值模擬 [J]. 水利學報, 2011,42(4):388-395. Jiang T, Zhu S L, Zhang Q, et al. Numerical simulation on effects of gate-pump joint operation on water environment in tidal river network [J]. Journal of Hydraulic Engineering, 2011,42(4):388-395.

    [14] 何民輝,舒 展.感潮河網(wǎng)地區(qū)引水增流改善水環(huán)境實踐探討 [J]. 南昌工程學院學報, 2013,32(3):64-67. He M H, Shu Z. Improving water environment by water diversion and flow strengthening in tidal river networks [J]. Journal of Nanchang Institute of Technology, 2013,32(3):64-67.

    [15] 賀新春,黃芬芬,汝向文,等.珠江三角洲典型河網(wǎng)區(qū)水資源調(diào)度策略與技術(shù)研究 [J]. 華北水利水電大學學報:自然科學版, 2016,6:55- 60. He X C, Huang F F, Ru X W, et al. Study on water resources dispatching strategy and technology in typical river network area of the Pearl River [J]. Journal of North China University of Water Resources and Electric Power, 2016,6:55-60.

    [16] 樓少華,呂權(quán)偉,任珂君,等.從深圳治水歷程研究高密度建成區(qū)排水系統(tǒng)的選擇與改造 [J]. 中國給水排水, 2018,34(18):28-31. Lou S H, Lyu Q W, Ren K J, et al. Study on the selection and reconstruction of urban drainage system in high density construction area from the course of water control in Shenzhen [J]. China Water & Wastewater, 2018,34(18):28-31.

    [17] 何振強,方詩標,陳永明,等.錢塘江感潮河段污染物遷移擴散數(shù)值分析 [J]. 環(huán)境科學學報, 2017,37(5):1668-1673. He Z Q, Fang S B, Chen Y M, et al. Numerical analysis of pollutant transport and diffusion in the tidal reach of Qiantang River [J]. Acta Scientiae Circumstantiae, 2017,37(5):1668-1673.

    [18] 季榮耀,陸永軍,詹小磊,等.伶仃洋茅洲河口動力地貌演變過程 [J]. 水科學進展, 2019,30(6):781-788. Ji R Y, Lu Y J, Zhan X L, et al. Study on the morphodynamic evolution processes in the Maozhou Estuary of the Lingding Bay [J]. Advances in Water Science, 2019,30(6):781-788.

    [19] 徐 群,莫思平,季榮耀,等.港珠澳大橋?qū)α尕暄蠛涌诔绷鳝h(huán)境的影響 [J]. 水利水運工程學報, 2012,2:79-83. Xu Q, Mo S P, Ji R Y, et al. Impacts of the Hong Kong-Zhuhai-Macao Bridge on flow environment of Lingdingyang estuary [J]. Hydro- Science and Engineering, 2012,2:79-83.

    [20] 趙彥軍,徐宗學,譚秋陽,等.深圳市降水特征時空演變規(guī)律分析 [J]. 北京師范大學學報(自然科學版), 2019,5:564-571. Zhao Y J, Xu Z X, Tan Q Y, et al. Spatiotemporal variations in precipitation in Shenzhen [J]. Journal of Beijing Normal University (Natural Science), 2019,5:564-571.

    [21] 馮存岸,于金珍,馮依蕾,等.陸海統(tǒng)籌下的陸源排放口聚類分析[J]. 水動力學研究與進展(A輯), 2020,35(3):42-48. Feng C A, Yu J Z, Feng Y L, et al. Cluster analysis of land-sourced sewage outlets based on Jaccard similarity [J]. Chinese Journal of Hydrodynamics, 2020,35(3):42-48.

    [22] 鄒永福,田 檢,黃曉航.因子分析和聚類分析在長江水污染綜合評價中的應(yīng)用 [J]. 科技信息, 2009,19:728-729. Zou Y F, Tian J, Huang X H. Application of factor analysis and ouster analysis to water quality comprehensive evaluation on Yangtze River [J]. Science & technology information, 2009,19:728- 729.

    [23] 沈迅偉,袁春偉.基于模糊聚類分析的水污染評價 [J]. 南京化工大學學報, 1998,20:30-34. Shen X W, Yuan C W. Water pollution assessment based on fuzzy cluster analysis [J]. Journal of Nanjing University of chemical technology, 1998,20:30-34.

    [24] 馬曉艷,唐 雁.層次聚類算法研究 [A]//2008年計算機應(yīng)用技術(shù)交流會論文集 [C]. 重慶:重慶計算機學會, 2008. Ma X Y, Tang Y. Research on hierarchical clustering algorithm [A]//Proceedings of 2008computer application technology exchange conference [C]. Chongqing: Chongqing Computer federation, 2008.

    [25] 尉景輝,何丕廉,孫越恒.基于K-Means的文本層次聚類算法研究 [J]. 計算機應(yīng)用, 2005,10:111-112. Yu J H, He P L, Sun Y H. Research on text hierarchical clustering algorithm based on K-Means [J]. Journal of computer applications, 2005,10:111-112.

    [26] Rossman L A, Storm water management model user’s manual version 5.0 [R].USA: United States Environmental Protection Agency, 2009.

    [27] Danish Hydraulic Institute. Water and environment. MIKE11Reference Manual [Z]. Denmark: DHI, 2014.

    [28] 黃國如,聶鐵鋒.廣州城區(qū)雨水徑流非點源污染特性及污染負荷 [J]. 華南理工大學學報:自然科學版, 2012,40(2):146-152. Huang G R, Nie T F. Characteristics and load of non-point source pollution of urban rainfall runoff in Guangzhou, China [J]. Journal of South China University of Technology (Natural Science Edition), 2012,40(2):146-152.

    [29] 羅鴻兵,羅 麟,黃 鵠,等.城市入河徑流排放口總污染特征研究 [J]. 環(huán)境科學, 2009,30(11):109-117. Luo H B, Luo L, Huang G, et al. Total pollution features of urban runoff outlet for urban river [J]. Environmental Science, 2009,30(11): 109-117.

    [30] 深圳市水務(wù)局.深圳市面源污染整治管控技術(shù)路線及技術(shù)指南(試行) [Z]. 深圳:深圳市水務(wù)局, 2018. Shenzhen Water Affairs Bureau: Shenzhen non point source pollution control technology route and technical guidelines [Z]. Shenzhen: Shenzhen Water Affairs Bureau, 2018.

    [31] 林希晨,倪紅珍,王 琳.東南沿海平原河網(wǎng)區(qū)域水質(zhì)水量模型研究 [J]. 水利水電技術(shù), 2019,50(6):150-157. Lin X C, Ni H Z, Wang L. Application of hydrodynamic and water quality model in the river network of southeast coastal plain area [J]. Water Resources and Hydropower Engineering, 2019,50(6):150-157.

    [32] 熊鴻斌,陳 雪,張斯思.基于MIKE11模型提高污染河流水質(zhì)改善效果的方法 [J]. 環(huán)境科學, 2017,38(12):5063-5073. Xiong H B, Chen X, Zhang S S. Method of Improving the Water Quality of Polluted Rivers Based on the MIKE11Model [J]. Environmental Science, 2017,38(12):5063-5073.

    [33] 樓少華,唐穎棟,陶 明,等.深圳市茅洲河流域水環(huán)境綜合治理方法與實踐 [J]. 中國給水排水, 2020,36(10):1-6. Lou S H, Tang Y D, Tao M, et al. Methods and practice of comprehensive improvement of Maozhou River water environment in Shenzhen [J]. China Water & Wastewater, 2020,36(10):1-6.

    致謝:感謝深圳市寶安區(qū)水務(wù)局提供的部分數(shù)據(jù).

    Spatial analysis and simulation study of rainfall runoff pollution for a tidal river network.

    ZHANG Feng-shan, SHANG Ming-zhu, ZHAO Peng-xiao, CHENG Kai-yu, TANG Ying-dong, WEI Jun*

    (Power China Huadong Engineering Corporation Limited, Hangzhou 311122, China)., 2021,41(4):1834~1841

    To explore the spatial distribution characteristics and pollution pattern of non-point source pollution of the tidal river network in the Maozhouhe River basin, the river network water quality under typical rainfall scenarios was simulated and analysed based on spatial analysis, statistical analysis and hydrodynamic-water quality coupling simulation method. An ecological water supply optimization strategy based on the target of water quality improvement was proposed. The results showed that the combination of hierarchical clustering aggregation algorithm and K-means can preferably distinguish the level and class of non-point source pollution. The water quality of the Shiyanqu River and the middle and upper reaches of the Songganghe River recovered slowly after rain fall due to the high non-point source pollution load and the lack of ecological water supply. A local optimized ecological water supply scheme was proposed based on the principle of constant amount of water replenishment. The recovery speed of water quality in key polluted rivers after rainfall was doubled by the optimized results, and the overall stability of river water quality in the basin was improved. The research conclusions provide support to the further understanding on the water pollution characteristics of the Maozhouhe River basin and the delicacy management of watershed water environment.

    tidal river networks;non-point source pollution;numerical modelling;Maozhouhe River

    X522

    A

    1000-6923(2021)04-1834-08

    張鳳山(1992-),男,內(nèi)蒙古豐鎮(zhèn)人,工程師,碩士,主要從事水環(huán)境數(shù)值模擬研究.發(fā)表論文8篇.

    2020-08-23

    廣東省重點領(lǐng)域研發(fā)計劃“污染防治與修復(fù)”重點專項(2019B110205005)

    * 責任作者, 正高級工程師, wei_j@ecidi.com

    猜你喜歡
    茅洲排口河網(wǎng)
    沿著茅洲河散步(組詩)
    連云港文學(2023年6期)2024-01-26 21:20:40
    基于小世界網(wǎng)絡(luò)的海河流域河網(wǎng)結(jié)構(gòu)及功能響應(yīng)
    水科學進展(2023年4期)2023-10-07 11:23:44
    多點網(wǎng)格采樣在燃煤機組總排口的應(yīng)用
    能源工程(2022年2期)2022-05-23 13:51:54
    Targeting cannabidiol to specific areas of the brain: an ultrasound-based strategy
    基于TMDL理念的流域排口污染物削減研究
    基于PSR模型的上海地區(qū)河網(wǎng)脆弱性探討
    天脊集團:總排口廢水深度處理裝置正式運行
    深圳:全流域治理共同系統(tǒng)治污
    環(huán)境(2019年2期)2019-03-29 01:16:00
    美了美了!茅洲河上賽龍舟
    南風窗(2018年12期)2018-06-12 04:25:24
    不同引水水源對平原河網(wǎng)影響分析
    国产熟女午夜一区二区三区| 看片在线看免费视频| 午夜激情av网站| 久久久久亚洲av毛片大全| АⅤ资源中文在线天堂| 两性夫妻黄色片| 黄色丝袜av网址大全| 国产亚洲精品av在线| 久久久久久免费高清国产稀缺| 91成人精品电影| 亚洲中文av在线| 99热这里只有精品一区 | 精品免费久久久久久久清纯| aaaaa片日本免费| 午夜久久久久精精品| 又黄又粗又硬又大视频| 国产精品久久久av美女十八| 亚洲性夜色夜夜综合| 日本一区二区免费在线视频| 淫妇啪啪啪对白视频| 麻豆成人av在线观看| 精华霜和精华液先用哪个| 成人亚洲精品一区在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 两个人看的免费小视频| 亚洲欧美精品综合一区二区三区| 国产爱豆传媒在线观看 | 成年人黄色毛片网站| 亚洲天堂国产精品一区在线| 国语自产精品视频在线第100页| 亚洲国产精品999在线| 黄色视频,在线免费观看| 亚洲成人久久爱视频| 黄色a级毛片大全视频| 成人特级黄色片久久久久久久| 久热这里只有精品99| 国产成人精品无人区| 国产单亲对白刺激| 成人国产综合亚洲| 国产免费av片在线观看野外av| 人成视频在线观看免费观看| 亚洲精品色激情综合| a级毛片a级免费在线| 两性夫妻黄色片| 国产精品一区二区免费欧美| 亚洲人成网站在线播放欧美日韩| 亚洲片人在线观看| 精品国产乱码久久久久久男人| 99热6这里只有精品| 日日摸夜夜添夜夜添小说| 麻豆成人av在线观看| 精品国产美女av久久久久小说| 久久天躁狠狠躁夜夜2o2o| 亚洲国产中文字幕在线视频| 成人国产综合亚洲| 欧美激情 高清一区二区三区| 久久中文字幕一级| 国产亚洲欧美精品永久| 女同久久另类99精品国产91| 色综合亚洲欧美另类图片| 在线观看午夜福利视频| 看免费av毛片| 日韩精品青青久久久久久| 国产午夜精品久久久久久| 国产熟女午夜一区二区三区| 窝窝影院91人妻| 日韩免费av在线播放| 亚洲av日韩精品久久久久久密| 久久久国产成人免费| 高潮久久久久久久久久久不卡| 欧美日韩黄片免| 日本精品一区二区三区蜜桃| 欧美人与性动交α欧美精品济南到| 啦啦啦观看免费观看视频高清| 妹子高潮喷水视频| 99在线人妻在线中文字幕| 久久久国产成人免费| 看片在线看免费视频| 成人国产一区最新在线观看| 国产精品爽爽va在线观看网站 | 亚洲中文字幕日韩| 亚洲精品美女久久久久99蜜臀| 久久午夜亚洲精品久久| 91在线观看av| 91成人精品电影| 亚洲七黄色美女视频| 特大巨黑吊av在线直播 | 91国产中文字幕| 亚洲人成网站高清观看| 国产精品久久久久久亚洲av鲁大| 亚洲av电影在线进入| 亚洲成a人片在线一区二区| 久久久久免费精品人妻一区二区 | bbb黄色大片| 中文字幕高清在线视频| av在线播放免费不卡| 啪啪无遮挡十八禁网站| 精品国产乱子伦一区二区三区| 黄色 视频免费看| 啦啦啦 在线观看视频| 可以在线观看的亚洲视频| 久久 成人 亚洲| 此物有八面人人有两片| 久久精品国产亚洲av高清一级| 一卡2卡三卡四卡精品乱码亚洲| 变态另类成人亚洲欧美熟女| 国产一卡二卡三卡精品| 亚洲色图 男人天堂 中文字幕| 亚洲人成网站高清观看| 午夜影院日韩av| 国产黄色小视频在线观看| 色综合欧美亚洲国产小说| 一二三四在线观看免费中文在| 国产私拍福利视频在线观看| 一本一本综合久久| 女人爽到高潮嗷嗷叫在线视频| 国产成人系列免费观看| 久久精品91无色码中文字幕| 成人18禁高潮啪啪吃奶动态图| 国产一级毛片七仙女欲春2 | 亚洲午夜精品一区,二区,三区| 一区二区日韩欧美中文字幕| 18禁裸乳无遮挡免费网站照片 | 日韩欧美一区二区三区在线观看| 国产av又大| 一级毛片高清免费大全| av欧美777| 久久久久九九精品影院| 亚洲国产欧美网| 午夜福利18| 两个人视频免费观看高清| 日本三级黄在线观看| 国产亚洲精品久久久久久毛片| 久久久水蜜桃国产精品网| 久久久精品国产亚洲av高清涩受| 中文字幕精品免费在线观看视频| avwww免费| 一边摸一边抽搐一进一小说| 国产又黄又爽又无遮挡在线| 老司机靠b影院| 国产精品影院久久| 亚洲精品国产精品久久久不卡| 免费在线观看成人毛片| 麻豆成人av在线观看| 变态另类丝袜制服| 女人被狂操c到高潮| 曰老女人黄片| 久久久水蜜桃国产精品网| e午夜精品久久久久久久| 亚洲国产精品合色在线| 日韩 欧美 亚洲 中文字幕| 亚洲五月天丁香| 欧美一级a爱片免费观看看 | 亚洲av电影在线进入| 婷婷丁香在线五月| 亚洲成人免费电影在线观看| 亚洲精品粉嫩美女一区| 亚洲一区高清亚洲精品| 日韩中文字幕欧美一区二区| 亚洲av日韩精品久久久久久密| 成年女人毛片免费观看观看9| 777久久人妻少妇嫩草av网站| 国产精品久久久久久精品电影 | 91成人精品电影| 中亚洲国语对白在线视频| 91老司机精品| 国产在线观看jvid| 男人操女人黄网站| 欧美日韩瑟瑟在线播放| 19禁男女啪啪无遮挡网站| 可以免费在线观看a视频的电影网站| 黄色成人免费大全| 嫩草影院精品99| 亚洲国产精品久久男人天堂| 国内少妇人妻偷人精品xxx网站 | 亚洲人成77777在线视频| 精品欧美一区二区三区在线| 国产亚洲精品久久久久久毛片| 成人亚洲精品一区在线观看| 成熟少妇高潮喷水视频| 国内少妇人妻偷人精品xxx网站 | 叶爱在线成人免费视频播放| 亚洲国产精品久久男人天堂| 大型黄色视频在线免费观看| 日日干狠狠操夜夜爽| 最新美女视频免费是黄的| 日韩欧美国产一区二区入口| 老鸭窝网址在线观看| 可以在线观看毛片的网站| 午夜福利18| 午夜免费激情av| 男男h啪啪无遮挡| 午夜福利欧美成人| 国产麻豆成人av免费视频| 手机成人av网站| 亚洲国产欧洲综合997久久, | 日韩欧美一区二区三区在线观看| 亚洲成人久久爱视频| 亚洲午夜精品一区,二区,三区| 老汉色av国产亚洲站长工具| 狂野欧美激情性xxxx| 动漫黄色视频在线观看| 免费搜索国产男女视频| 久久精品国产综合久久久| 国产主播在线观看一区二区| 别揉我奶头~嗯~啊~动态视频| 亚洲专区字幕在线| 在线免费观看的www视频| 欧美色欧美亚洲另类二区| av免费在线观看网站| 欧美激情极品国产一区二区三区| 国产成人欧美在线观看| 欧美一区二区精品小视频在线| 免费女性裸体啪啪无遮挡网站| 国产av一区在线观看免费| 国内少妇人妻偷人精品xxx网站 | 中文字幕另类日韩欧美亚洲嫩草| 国产又黄又爽又无遮挡在线| www.自偷自拍.com| 久久久久九九精品影院| 少妇熟女aⅴ在线视频| 免费高清视频大片| 成人18禁在线播放| 久久精品人妻少妇| 中亚洲国语对白在线视频| 精品久久久久久久人妻蜜臀av| 2021天堂中文幕一二区在线观 | 人人妻人人澡欧美一区二区| 欧美乱妇无乱码| 超碰成人久久| 欧美性猛交黑人性爽| 精品久久久久久久久久久久久 | 精品欧美国产一区二区三| 超碰成人久久| 国产精品自产拍在线观看55亚洲| 天堂影院成人在线观看| 一级a爱片免费观看的视频| 精品欧美国产一区二区三| 久久久久国产精品人妻aⅴ院| 久久香蕉国产精品| 欧美久久黑人一区二区| 亚洲片人在线观看| 久久精品国产99精品国产亚洲性色| 精品国产亚洲在线| 亚洲黑人精品在线| 国产精品亚洲美女久久久| www.999成人在线观看| 精品午夜福利视频在线观看一区| 国产精品综合久久久久久久免费| av超薄肉色丝袜交足视频| 男男h啪啪无遮挡| 亚洲成人国产一区在线观看| 国产精品亚洲av一区麻豆| 身体一侧抽搐| 在线观看午夜福利视频| 午夜激情av网站| 黄色视频,在线免费观看| 淫妇啪啪啪对白视频| 男人舔女人下体高潮全视频| 国产视频内射| 波多野结衣高清无吗| 18美女黄网站色大片免费观看| 91在线观看av| 国产精品爽爽va在线观看网站 | 男人舔女人下体高潮全视频| 日韩一卡2卡3卡4卡2021年| 亚洲中文字幕日韩| 国产精品 国内视频| 久久久久久久久中文| 美女高潮喷水抽搐中文字幕| 亚洲国产中文字幕在线视频| 久久人人精品亚洲av| 男人操女人黄网站| 精品国产亚洲在线| 黄片大片在线免费观看| 久久国产精品影院| 久久中文字幕一级| 91九色精品人成在线观看| 欧美黑人欧美精品刺激| 91大片在线观看| 美女免费视频网站| 麻豆成人午夜福利视频| 国产亚洲av高清不卡| 久久精品国产亚洲av高清一级| 人人妻人人看人人澡| 欧美黑人欧美精品刺激| 黑人欧美特级aaaaaa片| 婷婷六月久久综合丁香| 国产精品亚洲av一区麻豆| 大型黄色视频在线免费观看| 一区二区三区精品91| 国产精品99久久99久久久不卡| 亚洲av成人不卡在线观看播放网| 国产成人影院久久av| 日本熟妇午夜| 国产黄色小视频在线观看| 国产激情偷乱视频一区二区| 欧美激情 高清一区二区三区| 精品欧美一区二区三区在线| 国产精品影院久久| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲av五月六月丁香网| 成人亚洲精品一区在线观看| 欧美久久黑人一区二区| 啦啦啦免费观看视频1| 俺也久久电影网| 亚洲国产精品合色在线| 成人午夜高清在线视频 | 亚洲精华国产精华精| 日韩国内少妇激情av| 精品久久久久久久毛片微露脸| 欧美一级a爱片免费观看看 | 真人一进一出gif抽搐免费| av免费在线观看网站| 国产高清有码在线观看视频 | 国产精品免费视频内射| 欧美黑人巨大hd| 老司机深夜福利视频在线观看| 亚洲av日韩精品久久久久久密| 怎么达到女性高潮| 欧美黑人欧美精品刺激| 一本精品99久久精品77| 亚洲片人在线观看| 久久伊人香网站| 日韩欧美国产一区二区入口| 99精品久久久久人妻精品| 夜夜看夜夜爽夜夜摸| 午夜福利在线观看吧| 一区二区三区国产精品乱码| 国产一区二区在线av高清观看| 日韩欧美国产在线观看| 一边摸一边做爽爽视频免费| 亚洲aⅴ乱码一区二区在线播放 | 午夜福利欧美成人| 国产精品久久电影中文字幕| 一a级毛片在线观看| 久久草成人影院| 国内少妇人妻偷人精品xxx网站 | 视频在线观看一区二区三区| 黄色成人免费大全| 国产成人影院久久av| 侵犯人妻中文字幕一二三四区| 亚洲av熟女| 免费在线观看视频国产中文字幕亚洲| 日韩视频一区二区在线观看| 精品高清国产在线一区| 91成年电影在线观看| 国产激情欧美一区二区| 亚洲中文av在线| 久久热在线av| 亚洲,欧美精品.| 久久久久精品国产欧美久久久| 亚洲第一av免费看| 亚洲男人的天堂狠狠| 亚洲国产看品久久| 国产久久久一区二区三区| 亚洲国产看品久久| 99热这里只有精品一区 | 国产激情欧美一区二区| 欧美激情极品国产一区二区三区| 日日摸夜夜添夜夜添小说| 在线天堂中文资源库| 欧美一级a爱片免费观看看 | 一进一出抽搐动态| 日本五十路高清| 热99re8久久精品国产| 国产高清视频在线播放一区| 亚洲成人久久爱视频| 伊人久久大香线蕉亚洲五| 99久久国产精品久久久| 神马国产精品三级电影在线观看 | 亚洲国产欧美网| 黑丝袜美女国产一区| netflix在线观看网站| 欧美激情极品国产一区二区三区| 精品久久久久久久人妻蜜臀av| 成人免费观看视频高清| 精品卡一卡二卡四卡免费| 亚洲av日韩精品久久久久久密| 一区二区三区国产精品乱码| 成在线人永久免费视频| 亚洲最大成人中文| 大型av网站在线播放| 美女国产高潮福利片在线看| 满18在线观看网站| 亚洲专区国产一区二区| 欧美乱妇无乱码| 12—13女人毛片做爰片一| 青草久久国产| 18禁观看日本| 亚洲一区高清亚洲精品| 午夜a级毛片| 夜夜夜夜夜久久久久| 一边摸一边做爽爽视频免费| 欧美黑人巨大hd| 久久国产精品影院| 精品欧美国产一区二区三| 一进一出抽搐gif免费好疼| tocl精华| 中文字幕精品免费在线观看视频| 女生性感内裤真人,穿戴方法视频| 亚洲 欧美一区二区三区| 午夜a级毛片| 亚洲人成伊人成综合网2020| av天堂在线播放| 99久久精品国产亚洲精品| 亚洲成人精品中文字幕电影| 欧美另类亚洲清纯唯美| 在线永久观看黄色视频| 怎么达到女性高潮| 嫩草影院精品99| 欧美乱色亚洲激情| av超薄肉色丝袜交足视频| 在线观看免费视频日本深夜| 成年人黄色毛片网站| 亚洲精品中文字幕在线视频| 国产麻豆成人av免费视频| 日韩欧美国产一区二区入口| 桃色一区二区三区在线观看| 美女大奶头视频| 女人被狂操c到高潮| 日韩 欧美 亚洲 中文字幕| 色尼玛亚洲综合影院| 欧美激情 高清一区二区三区| 夜夜躁狠狠躁天天躁| 亚洲av成人av| 色综合婷婷激情| 制服丝袜大香蕉在线| 亚洲午夜精品一区,二区,三区| 18美女黄网站色大片免费观看| 日本黄色视频三级网站网址| 麻豆成人午夜福利视频| 国产亚洲精品第一综合不卡| 久久久久国产一级毛片高清牌| 国产一区在线观看成人免费| 欧美又色又爽又黄视频| 18禁黄网站禁片午夜丰满| 2021天堂中文幕一二区在线观 | 欧美黑人欧美精品刺激| 国产成人精品久久二区二区免费| 男女之事视频高清在线观看| 侵犯人妻中文字幕一二三四区| 成年人黄色毛片网站| 久久精品人妻少妇| 婷婷丁香在线五月| 亚洲第一青青草原| 一级毛片精品| 51午夜福利影视在线观看| 亚洲午夜精品一区,二区,三区| 日韩大码丰满熟妇| 国内少妇人妻偷人精品xxx网站 | 中文字幕精品亚洲无线码一区 | 精品一区二区三区四区五区乱码| 美女免费视频网站| 精品久久久久久久久久久久久 | 天堂影院成人在线观看| 女生性感内裤真人,穿戴方法视频| 久久久久免费精品人妻一区二区 | 好看av亚洲va欧美ⅴa在| 成人午夜高清在线视频 | 国产日本99.免费观看| 亚洲国产精品999在线| 欧美国产日韩亚洲一区| 不卡av一区二区三区| 久久亚洲真实| 亚洲成人久久性| 精品无人区乱码1区二区| 97人妻精品一区二区三区麻豆 | 国产精品久久久av美女十八| 中文字幕最新亚洲高清| 黄色视频不卡| 欧美国产精品va在线观看不卡| 日韩高清综合在线| 国产爱豆传媒在线观看 | 757午夜福利合集在线观看| www.www免费av| ponron亚洲| 性欧美人与动物交配| 两个人免费观看高清视频| 黄色a级毛片大全视频| 久久久精品欧美日韩精品| 亚洲va日本ⅴa欧美va伊人久久| 色精品久久人妻99蜜桃| 日本免费一区二区三区高清不卡| www国产在线视频色| www日本黄色视频网| 观看免费一级毛片| 一本精品99久久精品77| 午夜福利免费观看在线| 女人爽到高潮嗷嗷叫在线视频| 男女视频在线观看网站免费 | 久久香蕉国产精品| 久久久久国产精品人妻aⅴ院| 国产亚洲av嫩草精品影院| 久久久久亚洲av毛片大全| 精品久久久久久成人av| 一进一出好大好爽视频| 久久 成人 亚洲| av视频在线观看入口| 亚洲 欧美 日韩 在线 免费| 热99re8久久精品国产| av天堂在线播放| 国产成人啪精品午夜网站| 欧美黄色片欧美黄色片| av电影中文网址| 最近最新免费中文字幕在线| 亚洲全国av大片| 91麻豆av在线| 亚洲av美国av| 国产亚洲精品av在线| 久久天躁狠狠躁夜夜2o2o| 亚洲专区中文字幕在线| www.精华液| 国产一区二区激情短视频| 欧美激情久久久久久爽电影| 好男人在线观看高清免费视频 | 日本三级黄在线观看| 国产精品一区二区三区四区久久 | 哪里可以看免费的av片| 日韩国内少妇激情av| 国产色视频综合| 99国产极品粉嫩在线观看| 丰满人妻熟妇乱又伦精品不卡| 琪琪午夜伦伦电影理论片6080| 亚洲国产精品合色在线| 久久香蕉激情| 日韩欧美国产一区二区入口| 人妻久久中文字幕网| 女人高潮潮喷娇喘18禁视频| 欧美另类亚洲清纯唯美| 欧美黑人欧美精品刺激| 日韩国内少妇激情av| 久久午夜亚洲精品久久| 波多野结衣av一区二区av| 国产av一区在线观看免费| 亚洲中文av在线| 成年版毛片免费区| 亚洲一卡2卡3卡4卡5卡精品中文| 一a级毛片在线观看| 中文字幕av电影在线播放| 最近最新中文字幕大全免费视频| 亚洲中文字幕一区二区三区有码在线看 | 日本黄色视频三级网站网址| 欧美丝袜亚洲另类 | 精品国产一区二区三区四区第35| 50天的宝宝边吃奶边哭怎么回事| 亚洲av第一区精品v没综合| 国产99白浆流出| 精品午夜福利视频在线观看一区| 啦啦啦 在线观看视频| 亚洲午夜精品一区,二区,三区| 免费在线观看黄色视频的| 草草在线视频免费看| 99国产精品99久久久久| 中国美女看黄片| 神马国产精品三级电影在线观看 | 777久久人妻少妇嫩草av网站| 99久久国产精品久久久| 香蕉av资源在线| 国产亚洲精品一区二区www| 十八禁人妻一区二区| 久久国产乱子伦精品免费另类| 国产精品一区二区免费欧美| 久久精品国产亚洲av高清一级| 黄色 视频免费看| 老司机福利观看| 亚洲中文av在线| 老熟妇乱子伦视频在线观看| 久久这里只有精品19| 国产欧美日韩一区二区三| 国产成+人综合+亚洲专区| 真人一进一出gif抽搐免费| 国产成人系列免费观看| 国产av一区在线观看免费| 亚洲av成人av| 天天一区二区日本电影三级| 可以在线观看的亚洲视频| 国产精品久久久久久精品电影 | 欧美性猛交黑人性爽| 国产精品乱码一区二三区的特点| 中国美女看黄片| 午夜福利高清视频| 精品欧美一区二区三区在线| 国产精品,欧美在线| 国产激情偷乱视频一区二区| 男女之事视频高清在线观看| 精品福利观看| 后天国语完整版免费观看| av天堂在线播放| 欧美日韩一级在线毛片| 99国产极品粉嫩在线观看| 久久草成人影院| 两个人免费观看高清视频| 久久午夜综合久久蜜桃| 亚洲精品美女久久av网站| 亚洲av片天天在线观看| 满18在线观看网站| 亚洲中文字幕日韩| 国产国语露脸激情在线看| 日韩成人在线观看一区二区三区| 日本在线视频免费播放| 国产欧美日韩一区二区三| 成人一区二区视频在线观看| 免费高清在线观看日韩| 亚洲va日本ⅴa欧美va伊人久久| 免费看十八禁软件| 日韩欧美 国产精品| 欧美丝袜亚洲另类 | 免费在线观看成人毛片| 精品国产一区二区三区四区第35| 亚洲第一电影网av| 变态另类丝袜制服| 欧美绝顶高潮抽搐喷水| 最近最新免费中文字幕在线| 一区二区三区激情视频|