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

    城市人工淺水湖泊群水系連通方案研究

    2020-06-09 10:09:06志,馮
    人民長江 2020年4期
    關鍵詞:換水水系湖泊

    楊 志,馮 民 權

    (1.天津社會科學院 資源環(huán)境與生態(tài)研究所,天津 300191; 2.西安理工大學 省部共建西北旱區(qū)生態(tài)水利國家重點實驗室,陜西 西安 710048)

    城市大規(guī)模小型人工湖泊濕地群的建設,作為區(qū)域用水的新型載體,在生態(tài)效益、景觀效益、人文效益和社會效益等方面發(fā)揮重要作用,河湖水系連通作為一個新形勢下的治水方略已經(jīng)受到高度重視,逐漸成為被國內外學者關注的熱點問題。相關研究聚焦于水系連通的概念、指標、度量和應用,并廣泛借鑒了景觀生態(tài)學[1]、水文學[2]、地形學和經(jīng)濟學等學科領域的研究方法[3]。河湖水系連通的基本概念、分類體系[4]、理論范疇、演化制約[5],以及面臨的問題和挑戰(zhàn)等已被學者們大量探討。而河湖水系連通作為河網(wǎng)連通中的特殊單元,其典型特征在于水動力特性和污染物遷移擴散規(guī)律在河網(wǎng)一維特性到平面二維特性的突變?,F(xiàn)有研究多從湖泊水動力條件的改善出發(fā),探尋其影響因素[6],確定調水時機和調度方案[7],量化水系連通效果[8-9],評估水系連通性能[10],預判連通工程風險[11],或基于水系連通理論研究流域水資源合理配置[12]與動態(tài)管理[13]。而一維河網(wǎng)與二維湖泊的耦合仿真是上述問題的研究基礎,國內外學者從徑流形成的水文過程著手,開展河湖水系的水動力條件模擬研究[14-15]。在此基礎上,考慮水系連通對河湖之間污染物質的遷移擴散作用,蘇苑君等[16]評估了大東湖水系建成前后外沙湖和水果湖水質的變化情況,楊衛(wèi)等[17]從水動力學、水質、社會經(jīng)濟學等方面建立了水系聯(lián)通方案的綜合評價體系。

    上述研究在水系連通性和水量水質調控領域取得了長足的進展,然而仍存在以下兩方面問題值得深入。一是水質研究受不同水域污染物質成分和水生態(tài)轉化過程的共同作用,因地制宜建立相適應的生態(tài)動力學模型尤為重要。二是在遷移、擴散、轉化體系下,將水生境結構相對簡單的城市人工湖泊置于多閘壩河網(wǎng)之中,其人工河湖水系的水體交換方式與水質控制措施,仍有待進一步探討。基于此,本文以淮河流域沙潁河水系污染最為嚴重的清潩河流域為例,構建流域水生態(tài)模型,分析不同河湖水體交換方式下,城市人工湖泊的流場及多種污染物質濃度場的分布特征及變化規(guī)律;并根據(jù)污染物質分布及削減量,優(yōu)化清潩河流域河湖水系連通方案。

    1 研究區(qū)域概況及數(shù)據(jù)來源

    清潩河位于河南省境內,屬淮河流域沙潁河水系,全長149 km,流域面積2 362 km2。清潩河流域許昌段有人工湖泊及濕地8個,其中,北海為該段規(guī)模最大、邊界最復雜的人工湖泊。以北海為例開展城市人工湖的流場及水質濃度場的變化規(guī)律分析,再以北海研究結果為依據(jù),擬定清潩河流域湖泊群的整體水系連通方案,優(yōu)化城市人工河湖水系連通性調控措施,盤活清潩河流域水系,提升流域水環(huán)境質量。湖泊位置分布見圖1,基本情況見表1。

    表1 清潩河流域湖泊基本參數(shù)

    水系圖通過來源于地理空間數(shù)據(jù)中心的DEM數(shù)據(jù)自動提取,并根據(jù)河南省地圖院提供的水系圖進行校正,河網(wǎng)斷面和湖泊地形數(shù)據(jù)通過委托監(jiān)測獲得。氣象數(shù)據(jù)采用許昌氣象局提供和rp5.ru天氣網(wǎng)下載的許昌站(東經(jīng)113°52’、北緯34°02’)監(jiān)測資料。水文和水質數(shù)據(jù)來源于水文年鑒、重點污染源在線監(jiān)測、環(huán)境統(tǒng)計年鑒和排污口調研。

    2 研究方法

    2.1 流域水動力-水質-水生態(tài)模型

    流域水動力-水質-水生態(tài)模型將一維河道、河網(wǎng)汊點和二維湖泊嵌套,以表達強人工干擾流域的河湖連通關系;并將簡單的水質方程和二次開發(fā)的復雜生態(tài)動力學模型相耦合,模擬污染物質的輸移、擴散、降解和轉化。

    一維河網(wǎng)水動力水質模型采用一維非恒定流Saint-Venant方程組和水質方程構建,河網(wǎng)汊點單元水流運動遵循質量守恒原理和能量守恒原理[18],采用Abbott-Ionescu六點隱式有限差分格式求解方程組,進一步采用追趕法求解離散后的線形方程組。汊點單元通過一條無斷面的短河流將支流與干流連接,僅有3個計算點,水位點-流量點-水位點,并采用相同的河底高程連接水位點-水位點,以保證模型計算的穩(wěn)定性??刂品匠倘缦?

    (1)

    (2)

    式中,A為過水斷面面積,m2;t為時間,s;Q為過流流量,m3/s;x為空間坐標,m;q為旁側入流流量,m3/s;α為動量校正系數(shù);h為水位,m;g為重力加速度,m/s2;C為謝才系數(shù);R為水力半徑;C0為濃度,mg/L;D為擴散系數(shù);K為降解系數(shù),1/s;C2是源、匯項濃度,mg/L。

    二維湖泊水動力模型基于三角形非結構網(wǎng)格,由水流連續(xù)性方程、水流運動方程和水質方程構建,并采用單元中心有限體積法求解。一、二維模型可根據(jù)“水位-流量”銜接關系進行標準連接耦合[19],鑒于清潩河流域河網(wǎng)復雜、相鄰斷面高程差、河道型湖泊體量小、河湖連接方式簡單等原因,本次研究簡化耦合過程,將湖泊源、匯概化為出流與入流邊界嵌套于一維河網(wǎng)模型參與計算。

    水流連續(xù)方程:

    (3)

    水流運動方程:

    (4)

    水質方程:

    (5)

    水生態(tài)模型包括7個狀態(tài)變量、30個常數(shù)、6個作用力、13個輔助變量等。7個狀態(tài)變量為溶解氧、生化需氧量、氨氮、硝酸鹽氮、亞硝酸鹽氮、總氮和總磷,6個作用力包括風場、溫度、鹽度、水深、當前流速和分層數(shù),除溫度和鹽度要從外界進行輸入之外,其他均從水動力模型中獲取。

    (1)溶解氧。

    phtsyn·F(N,P)-respT-sod

    (6)

    reaera=K2(CS-DO)

    (7)

    (8)

    phtsyn=

    (9)

    (10)

    (11)

    (12)

    上述式中,

    reaera代表大氣復氧;phtsyn指光合作用過程中的實際產氧量,g/(m2·d);respT代表生物呼吸效率,g/(m2·d);BODd為BOD降解需氧量,nitriDO為硝化作用需氧量,sod為底泥需氧量;Y1為復氧系數(shù);F(N,P)為光和作用營養(yǎng)鹽限制函數(shù);K2為大氣復氧速率,1/d;Cs為水體中飽和狀態(tài)下的溶解氧的濃度,mg/L;K4為20℃時硝化率,1/d;NH3為氨氮的濃度,mg/L;θ4為硝化作用的溫度系數(shù);DO為實際溶解氧濃度,mg/L;HS-nitr為硝化作用的半飽和濃度,mg/L;Pmax為中午最大產氧量,g/(m2·d);F1(H)為光衰減函數(shù);τ為中午的實際時間;α為實際相對日長;tup,tdown分別為日出和日落時間;R1為自養(yǎng)生物在20℃下的光合作用呼吸效率,g/(m2·d);θ1為光合作用和呼吸作用的溫度系數(shù);R2為動物和細菌等異養(yǎng)生物的呼吸效率,g/(m2·d);θ2為異養(yǎng)呼吸的溫度系數(shù);K3為20℃時有機物的降解系數(shù),1/d;θ3為阿列紐斯溫度系數(shù);BOD為生化需氧量濃度,mg/L;HS-BOD為BOD的半飽和氧濃度,mg/L;HS-SOD為SOD的半飽和氧濃度,mg/L;θ7為SOD的溫度系數(shù);IN為無機氮總和,mg N/L;KSN為限制植物和藻類光合作用的氮的半飽和濃度,mg N/L;PO4為磷酸根濃度,mg P/L;KSP為限制植物和藻類光合作用的磷的半飽和濃度,mg P/L;S為鹽度,ppt;T為溫度,℃;V為水深的平均流速,m/s;H為水深,m;WV為風速,m/s;k為消光系數(shù)。

    (2)生化需氧量。

    (13)

    (3)氨氮。

    (14)

    (15)

    (16)

    (17)

    (18)

    (19)

    式中,BODn為降解BOD釋放氨氮的過程;plantn為植物吸收氨氮的過程;bactn為細菌分解氨氮的過程;nitrif為氨氮轉化為亞硝酸鹽的過程;hetersn為異養(yǎng)呼吸釋放氨氮的過程;YBOD為BOD降解釋放的氨氮,mg NH3-N/mg BOD;UNP為被植物吸收的氨氮,mg N/mg O2;UNb為被細菌吸收的氨氮,mg N/mg BOD;HS-NH3為細菌吸收氨氮半飽和濃度,mg N/L;K4為20℃時的硝化率,1/d;θ4為硝化作用的溫度系數(shù)。

    (4)亞硝酸鹽氮。

    (20)

    (21)

    (22)

    式中,nitri為亞硝酸鹽轉化為硝酸鹽的過程;HS-nitr為硝化作用的半飽和濃度,mg O2/L;NO2為亞硝酸鹽的濃度,mg/L;K5為20℃時NO2轉化為NO3效率,1/d;θ5為NO2轉化為NO3的溫度系數(shù)。

    (5)硝酸鹽氮。

    (23)

    (24)

    式中,deni為反硝化過程;K6為反硝化速率,1/d;θ6為阿列紐斯溫度系數(shù)。

    (6)磷酸鹽。

    (26)

    (26)

    (27)

    (28)

    (29)

    式中,BODp為生化需氧量衰減釋放磷的過程;plantp為植物吸收磷的過程;bactp為細菌吸收磷的過程;hetersp為異養(yǎng)呼吸釋放磷的過程;Y2為溶解態(tài)BOD中磷的含量,mg P/mg BOD;UPp為植物吸收的磷,g P/m3·d-1;UPb為細菌吸收的磷,g P/m3·d-1;HS-PO4為細菌吸收的磷酸鹽的半飽和濃度,mgO2/L。

    (7)總氮。

    (30)

    2.2 邊界條件

    采用非結構網(wǎng)格,將北海模擬區(qū)域剖分為3 634個網(wǎng)格(見圖2),并通過調整邊界節(jié)點,控制湖泊邊界處和湖心島處的網(wǎng)格密度和大小,根據(jù)實測湖泊高程點資料,插值湖泊地形如圖3所示。由于清潩河流域湖泊群普遍具有狹長型城市人工淺水湖泊的特征,且均位于河道之上或與河道連通,具有河道型湖泊的特點;又因流域內湖泊建成通水時間較短,研究過程中未能獲取完備的數(shù)據(jù)資料;因此后續(xù)選取數(shù)據(jù)基礎更為充實的河網(wǎng)進行水生態(tài)模型參數(shù)敏感性分析和水動力水質水生態(tài)模型驗證。水動力水質水生態(tài)模型中,由于水質模型在研究過程中,對COD表現(xiàn)出較好的表達效果,因此COD僅考慮溶質輸移和簡單降解作用,其他指標采用水生態(tài)模型表達。

    圖2 北海網(wǎng)格劃分

    圖3 北海地形

    河網(wǎng)入流邊界:天然徑流采用水文模擬得到的流量和非點源污染物質濃度[20];污水處理廠采用逐日實測出水流量和濃度過程;工業(yè)點源采用逐日實測出水流量過程和水質濃度過程,以及環(huán)境統(tǒng)計報表數(shù)據(jù)計算年平均值計入模型;上述數(shù)據(jù)中缺失的水質組份數(shù)據(jù)采用許昌市流域斷面監(jiān)測逐月數(shù)據(jù)計入模型。河網(wǎng)出流邊界:采用斷面水位流量關系,并采用國控、省控斷面的實測水質數(shù)據(jù)進行校正。湖泊入流邊界:水量邊界采用擬定的換水和循環(huán)流量方案,水質邊界采用河道供水斷面的水質計算結果。湖泊出流邊界:水量邊界采用恒定水位控制,水質邊界采用零梯度邊界。閘壩節(jié)點根據(jù)閘門過水形式,主要包括閘下過水和閘上過水兩種,分為泄流閘和越流閘兩類,分別控制其閘門開啟度和壩頂高程。

    2.3 參數(shù)率定與模型驗證

    根據(jù)清潩河流域2014年(率定期)和2015年(驗證期)的實測流量和COD濃度率定和驗證水動力水質模型的參數(shù)值。清潩河流域糙率和污染物質的擴散和降解系數(shù)見表2。由于清潩河流域無實測糙率資料,其取值通過數(shù)值實驗率定求取,率定值較類似河流偏大,主要是由于實際情況掌握有限,模擬中未能盡現(xiàn)諸多閘壩的水流阻礙作用,一定程度上由糙率得到彌補與平衡。清潩河流域平均流速0.15 m/s,擴散系數(shù)隨流速動態(tài)變化的關系式為D=aVb,率定可知當a=125,b=1,最小值和最大值分別限制在5~20之間,模擬誤差最小。采用高村橋斷面2014年的模擬結果平均值,對水生態(tài)模型的30個參數(shù)進行敏感性分析,每次僅改變一個參數(shù),變化范圍為基準值上下調整50%,進一步對高敏感度的參數(shù)進行率定,結果見表3。

    表2 清潩河流域水動力-水質模型參數(shù)率定結果

    高村橋斷面率定期和驗證期的流量模擬誤差分別為9.40%和8.07%,相關系數(shù)分別為0.92和0.78,見圖4;COD濃度各斷面誤差平均值為16.05%,見表4。其他指標個別月份實測值與年平均水平相差極大,去除奇異值,整體吻合程度較好。BOD整體上吻合程度好,平均相對誤差為11.71%,一致呈現(xiàn)下降趨勢,在3月份的誤差最大;溶解氧在2月、4月和5月份的實測值分別低于年平均水平49.78%,46.19%和58.74%,除去這3個月份之外,整體相對誤差為15.01%;總氮在6月份的實測值低于年平均水平71.58%,除6月份外平均相對誤差為21.59%;總磷在整體上吻合較好,平均相對誤差為8.54%,4月份誤差較大;氨氮年總負荷誤差為2.48%。因此,所建模型在極值的捕捉上有所偏差,但整體上能夠較好模擬清潩河流域的水生態(tài)變化規(guī)律。

    表3 水生態(tài)模型參數(shù)敏感性分析及率定結果

    圖4 高村橋流量模型驗證結果

    表4 COD降解系數(shù)誤差分析

    3 結果與分析

    對于水資源量缺乏的城市人工湖系統(tǒng),以換水方式改善景觀蓄水水體水質,較為經(jīng)濟節(jié)水。但周期性換水雖然一定程度上改善了小型人工湖泊封閉式的水體狀態(tài),卻無法從根本上建立起穩(wěn)定循環(huán)的湖泊系統(tǒng)。鑒于此,本研究一方面針對清潩河流域現(xiàn)有湖泊換水方案進行優(yōu)化,另一方面從恒定流量和波動流量兩種類型探索湖泊與河網(wǎng)的穩(wěn)定、循環(huán)的水體交換方案。

    3.1 不同頻率換水方案

    清潩河流域人工河湖水系規(guī)劃現(xiàn)狀(保證率50%年份)換水次數(shù)為每年5次,分別于1月、4月、7月、9月和11月進行換水,換水流量為2 m3/s,一次換水大約4 d。其他年份每年換水4次,特殊枯水年份減少至2次。根據(jù)上述實際情況,本節(jié)擬定每年于3月、7月和11月?lián)Q水3次,和分別于1,4,7,9月和11月?lián)Q水5次兩種工況,分析不同換水頻率對人工湖水質的影響。本節(jié)針對化學需氧量、氨氮、總磷和總氮4項指標分析其調控效果。

    表5 2種工況下出水水質年均值和污染物質削減水平對比

    圖6 2種工況下的COD、氨氮和TP年均分布

    換水3次和5次時的北海流速平均值分別為3.2 mm/s和3.4 mm/s,湖內分布如圖5所示,可見,易形成環(huán)流的湖泊上游分布規(guī)律較為相似,差別主要體現(xiàn)在較為順直的湖泊下游。不同換水頻率下的湖泊內污染物質濃度分布情況如圖6所示。COD平均年濃度分別為13.03 mg/L和14.48 mg/L,氨氮平均年濃度分別為0.32 mg/L和0.42 mg/L,換水5次均高于換水3次。TP平均年濃度分別為0.359 mg/L和0.357 mg/L,換水5次略低于換水3次。換水頻率的增加導致湖泊水動力狀況略有改善,但并不一定可以改善水體水質狀況。從湖泊出水平均濃度看來,換水3次的出水水質優(yōu)于換水5次,說明換水3次過程中進入河網(wǎng)的水體水質較好(見表5)。從換水對污染物質總量削減情況看來,換水3次削減量大于換水5次(見表5)。綜上所述,換水3次更有利于污染負荷和水系水質控制。

    3.2 恒定循環(huán)流量方案

    鑒于周期性換水沖擊穩(wěn)定湖泊系統(tǒng)的弊端,本節(jié)提出分別以0.1 m3/s(工況1)、0.2 m3/s(工況2)、0.3 m3/s(工況3)、0.4 m3/s(工況4)和0.5 m3/s(工況5)的恒定小流量持續(xù)循環(huán)進出湖泊的水系連通方式。根據(jù)5種工況下的流速平均值分布,流速整體分布隨入湖流量增加而變大,且主要集中在環(huán)流區(qū)和收縮段。表6為5種工況下的流速平均值統(tǒng)計結果,入湖流量在0.1~0.5 m3/s,流速集中于0.001 ~0.1m/s之間,其中0.005m/s以下占比最大。綜上可見,控制入湖流量是改變流速值分布的有效手段,增加入湖流量可提高湖泊整體流速水平。

    圖5 2種工況下流速年均分布

    表6 5種工況下的流速平均值統(tǒng)計結果

    圖7為進湖流量分別為0.1,0.3 ,0.5 m3/s時的北海局部流場分布。在北海中選取兩個局部區(qū)域進行重點分析,A區(qū)為環(huán)流區(qū),B區(qū)為近似河道式湖泊,包含收縮段和湖心島。A區(qū)的環(huán)流結構隨入湖流量增加,逐漸從一個中心環(huán)流加西北部環(huán)流跡象(a1)發(fā)展為南北兩個環(huán)流(a3),說明入湖流量可改變環(huán)流中心位置與環(huán)流結構。B區(qū)隨入湖流量增加,破壞了前段的局部環(huán)流結構,隨后的收縮段流速明顯增大,湖心島北部流速隨南部主流增加而降低,說明入湖流量可明顯改變河道型湖泊的流場分布。根據(jù)湖泊上6個觀測點的瞬時流速值,t1在環(huán)流中心,流速呈增加趨勢;t2和t3處由于局部環(huán)流結構的改變,流速隨機性變化;t4在環(huán)流區(qū)邊界處,流速明顯大于其他區(qū)域;t5為收縮段,流速明顯增加;t6為湖心島,環(huán)流強度隨主流增加而降低。

    圖7 進湖流量分別為0.1,0.3 ,0.5 m3/s時的北海局部流場分布

    綜上所述,流速整體分布隨入湖流量增加而變大,且主要集中在環(huán)流區(qū)和收縮段,入湖流量可改變環(huán)流區(qū)的環(huán)流中心位置與環(huán)流結構,可明顯改變河道型湖泊的流場分布。隨入湖流量的增加,環(huán)流中心流速和收縮段流速明顯增加,湖心島環(huán)流強度隨之降低。控制入湖流量是改變流速值分布的有效手段,增加入湖流量可提高湖泊整體流速水平。

    根據(jù)湖泊內COD,DO,TP和氨氮的年平均濃度(見圖8),可見進出湖流量從0.1 m3/s 增至0.2 m3/s,湖內污染物質濃度變幅較大,整體上變幅隨進出湖流量增加而逐漸減小。5種工況下整個區(qū)域在0.1~0.5 m3/s下的DO平均水平分別為4.76,4.30,4.20,4.23,4.29 mg/L,進出湖流量為0.1 m3/s時湖內DO水平最高,湖泊DO水平與入湖水體DO水平相關。整體上COD和氨氮污染負荷隨進出湖流量增加而加重,5種工況下COD年平均濃度分別在16.00~18.77 mg/L之間,氨氮年平均濃度在0.56 ~0.96 mg/L之間。說明恒定小流量循環(huán)受進湖水體污染物質水平影響較大,換水方式下由于換水期流量遠遠大于點源污染排水流量,且周期較短,入湖污染物質負荷較低,對湖水水質影響較小。恒定小流量循環(huán)使得調水水體與工業(yè)點源水體充分混合進入湖泊,導致湖泊水體受到一定影響。整體上TP污染負荷隨進出湖流量增加而得到改善,5種工況下TP年平均濃度分別為0.46,0.43,0.41,0.40 mg/L和0.39 mg/L。

    圖8 5種工況下湖泊內的污染物質濃度年均值

    表7為5種工況下出湖水體水質平均水平,COD和氨氮隨循環(huán)流量增加而增加,TN和TP隨之減少。圖9為5種循環(huán)流流量下COD和氨氮污染物質總量的削減情況,可見0.1 m3/s工況下COD削減能力明顯低于其他幾種工況,0.2~0.5 m3/s 工況下COD削減量在26.45~26.91 t/a,差距并不明顯,但0.3 m3/s工況對COD削減效果略好。5種工況下的氨氮年削減量分別為5.10,6.89,3.63,1.79,0.82 t/a,其中0.2m3/s工況下的氨氮削減能力顯著于其他工況。

    表7 湖泊出水污染物濃度年均值

    圖9 5種工況下污染物質削減水平對比

    綜合上述分析,湖泊水動力水平隨進出湖流量增加而增強;湖內及出湖水體的COD和氨氮濃度隨循環(huán)流量增加而增加,TN和TP隨之減小,DO變化趨勢不明確;5種工況對COD和氨氮污染物質總量的削減水平對比中,循環(huán)流量為0.2 m3/s工況最佳;結合進出湖水量的來源和經(jīng)濟等因素,可考慮采用0.2 m3/s的流量進出湖泊,將之與流域水系連通循環(huán)。

    3.3 波動循環(huán)流量方案

    根據(jù)3.2計算結果,采用0.2 m3/s的流量進出北海,將之與流域水系連通循環(huán)較為合理。本節(jié)設計相同水量消耗下,采用波動流量進出湖泊的3種工況,分別為以0.4 m3/s單數(shù)周供為(工況1)、單數(shù)月供水(工況2),以及1,2,5,6,9,10月份供水(工況3)。統(tǒng)計模擬期內,湖泊內水動力和水質的平均水平,并將之與進出水為恒定流量0.2 m3/s工況下的統(tǒng)計結果進行對比,結果見表8,可見波動進水工況下該湖泊水動力和水質條件有所改善,但效果并不顯著。整體看來,工況3對湖泊水動力和水質改善效果最佳。

    表8 湖內污染物濃度年均值對比

    根據(jù)出湖水體COD和氨氮濃度變化過程(見圖10)可知,波動進水使湖內污染物質濃度呈現(xiàn)周期性變化規(guī)律,且一個周期內的最大值和最小值的差值隨周期長度增加而增加,說明波動進水在一定程度上破壞了污染物質在湖內的降解過程。3種工況下出湖水體的COD濃度年平均值分別為16.00,16.22 mg/L和15.94 mg/L,氨氮濃度年平均值分別為0.55,0.54 mg/L和0.54 mg/L,工況3效果最佳,但3種工況均大于恒定流量0.2 m3/s工況下的污染物濃度平均值。

    根據(jù)污染物年削減量統(tǒng)計結果(見表9),3種工況均沒有0.2 m3/s削減效果好,但工況1比工況2和工況3削減量大。綜合上述分析,相同水量消耗下,恒定進湖流量和波動進湖流量對湖內及出水的水動力和水質條件改善效果相差不大,但對污染物的削減效果,恒定進水效果更佳,隔周進水的工況次之??紤]到進水流量的調節(jié)所消耗的人員和經(jīng)濟成本,恒定流量進出湖的循環(huán)模式更具優(yōu)勢。

    圖10 出水污染物濃度年均值

    表9 污染物年削減量均值

    4 清潩河流域湖泊群調控方案討論

    根據(jù)前述恒定流量循環(huán)模式與換水模式中的最優(yōu)工況,本節(jié)擬定循環(huán)+換水的湖泊群調控措施,即按照0.2 m3/s的恒定小流量持續(xù)循環(huán)進出湖泊,并于每年的3,7月和11月?lián)Q水3次的水系連通方式。該工況下的水動力和水質平均水平的分布情況如圖11所示,湖內流速年平均值為3.62 mm/s,比0.2 m3/s進湖工況和3次換水工況分別增加了10.78%和14.97%。湖內和出水各項水質指標年均值如表10所示,表明湖內的COD和氨氮平均濃度分別比0.2 m3/s進湖工況增加了1.21%和3.10%,但TP減少了11.26%。出水COD和氨氮平均濃度比0.2 m3/s進湖工況略有增加,但TP大幅降低。COD的年削減量為31.20 t,比0.2 m3/s進湖工況增加了17.48%;氨氮的年削減量為6.65 t,比0.2 m3/s進湖工況降低了3.46%;總磷的年增加量為0.67 t,比0.2 m3/s進湖工況降低了20.81%。

    圖11 北海調控后各水質指標年均分布

    表10 湖內和出水各項水質指標年均值

    綜合上述分析,循環(huán)+換水模式對湖內和出水水體中不同污染物的控制水平上有所差異,整體改善水平與0.2 m3/s進湖工況相差不大,但就污染物質負荷的削減量來看,循環(huán)+換水模式效果顯著。根據(jù)該結果,結合清潩河上8個人工湖泊的蓄水量差別,分別擬定各湖泊的進水方式如表11所示,清潩河流域湖泊群的循環(huán)流量在0.03~0.20 m3/s之間,并于每年的3,7月和11月統(tǒng)一換水3次。清潩河流域水資源短缺,生態(tài)用水需充分考慮經(jīng)濟因素,因此上述方案在理論上相對保守。在實際河湖水系連通過程中,可根據(jù)不同河段的供退水能力差異加以調整,更大限度的保障水生態(tài)需求。

    表11 清潩河流域湖泊進水方式優(yōu)化

    5 結 論

    將水生態(tài)模型應用于河湖水系連通及其水量水質調控研究,分析了不同的換水頻率、不同恒定流量和不同波動流量等換水方式下人工湖泊的流場及水質濃度場,根據(jù)污染物質濃度分布及削減量提出了循環(huán)+換水的湖泊群調控模式,得到主要結論如下。

    (1)在參數(shù)敏感性分析的基礎上,率定了水生態(tài)模型中9個重要參數(shù)值并根據(jù)實測資料驗證了模型的準確性。

    (2)換水頻率的增加導致湖泊水動力狀況略有改善,但并不一定可以改善水體水質狀況,北海換水3次更有利于污染負荷和水系水質控制。湖泊水動力水平隨進出湖流量增加而增強,湖內及出湖水體的COD和氨氮濃度隨循環(huán)流量增加而增加,TN和TP隨之減小,DO變化趨勢不明確,河湖水系連通循環(huán)流量為0.2 m3/s時,對北海COD和氨氮污染物質總量的削減水平最佳。相同水量消耗下,恒定進湖流量和波動進湖流量對湖內及出水的水動力和水質條件改善效果相差不大,但對污染物的削減效果,恒定進水效果更佳,隔周進水的工況次之,考慮到調節(jié)方式所消耗的人員和經(jīng)濟成本,恒定流量進出湖的循環(huán)模式更具優(yōu)勢。

    (3)確定了循環(huán)+換水的湖泊群調控模式,該模式對湖內和出水水體中不同污染物的控制水平上有所差異,整體改善水平與0.2 m3/s進湖工況相差不大,但就污染物質負荷的削減效果顯著。量化并優(yōu)化了清潩河流域各湖泊的循環(huán)流量在0.03~0.20 m3/s之間,換水次數(shù)為每年的3,7月和11月統(tǒng)一換水3次。

    猜你喜歡
    換水水系湖泊
    杭州亞運泳池不換水的“黑科技”
    軍事文摘(2023年24期)2023-12-19 06:49:50
    鄱陽湖水系之潦河
    水培植物霎換水
    自動換水
    環(huán)水系旅游方案打造探析——以臨沂市開發(fā)區(qū)水系為例
    你相信嗎?湖泊也可以“生死輪回”
    水系魔法之止水術
    “害羞”的湖泊
    奇異的湖泊
    我給小烏龜換水
    免费久久久久久久精品成人欧美视频| 成人手机av| 精品国产一区二区三区久久久樱花| 深夜精品福利| h视频一区二区三区| 国产熟女欧美一区二区| 七月丁香在线播放| 人人澡人人妻人| 国产精品免费视频内射| 免费看不卡的av| 大型av网站在线播放| 操美女的视频在线观看| 一本大道久久a久久精品| 欧美日韩av久久| 啦啦啦中文免费视频观看日本| 在线观看国产h片| 亚洲午夜精品一区,二区,三区| 亚洲伊人色综图| www.999成人在线观看| 亚洲一区二区三区欧美精品| 在线观看国产h片| 亚洲专区国产一区二区| 精品一品国产午夜福利视频| 久久精品国产亚洲av高清一级| 十分钟在线观看高清视频www| 久久国产亚洲av麻豆专区| 丰满饥渴人妻一区二区三| 人妻 亚洲 视频| 亚洲av日韩在线播放| 女人精品久久久久毛片| av电影中文网址| 久久精品国产a三级三级三级| 国产三级黄色录像| 午夜老司机福利片| netflix在线观看网站| 亚洲精品国产av蜜桃| 男女床上黄色一级片免费看| 十八禁人妻一区二区| 激情五月婷婷亚洲| 国产精品一二三区在线看| 午夜福利一区二区在线看| 大话2 男鬼变身卡| 一边亲一边摸免费视频| 国产精品一区二区精品视频观看| kizo精华| 一区在线观看完整版| 久久综合国产亚洲精品| 美女视频免费永久观看网站| 国产亚洲欧美在线一区二区| 欧美日韩一级在线毛片| 天天躁日日躁夜夜躁夜夜| 每晚都被弄得嗷嗷叫到高潮| 亚洲伊人色综图| 男人舔女人的私密视频| 日本五十路高清| 亚洲av日韩精品久久久久久密 | 久久久久国产一级毛片高清牌| 精品久久久精品久久久| 精品人妻一区二区三区麻豆| 大话2 男鬼变身卡| 国产精品三级大全| 又紧又爽又黄一区二区| 欧美激情极品国产一区二区三区| 别揉我奶头~嗯~啊~动态视频 | 日本欧美国产在线视频| 国产野战对白在线观看| 亚洲欧美精品自产自拍| av又黄又爽大尺度在线免费看| 天天躁日日躁夜夜躁夜夜| 精品人妻1区二区| 日韩熟女老妇一区二区性免费视频| 国产一区二区三区av在线| 男女边吃奶边做爰视频| 超碰97精品在线观看| 深夜精品福利| 久久综合国产亚洲精品| 女人精品久久久久毛片| av又黄又爽大尺度在线免费看| 成人国语在线视频| 一边亲一边摸免费视频| 97人妻天天添夜夜摸| 国产精品香港三级国产av潘金莲 | 亚洲国产精品成人久久小说| 另类亚洲欧美激情| 日本av免费视频播放| 国产精品一区二区在线不卡| 97人妻天天添夜夜摸| 人妻 亚洲 视频| 1024视频免费在线观看| 我的亚洲天堂| 自线自在国产av| 看免费成人av毛片| 日本wwww免费看| 中文字幕高清在线视频| videos熟女内射| 成在线人永久免费视频| 免费不卡黄色视频| 欧美日韩国产mv在线观看视频| 男女午夜视频在线观看| 国产福利在线免费观看视频| xxx大片免费视频| 亚洲自偷自拍图片 自拍| 国产亚洲av片在线观看秒播厂| 香蕉丝袜av| 成人国产av品久久久| 汤姆久久久久久久影院中文字幕| 国产日韩一区二区三区精品不卡| 亚洲av综合色区一区| 久久精品国产亚洲av涩爱| 中文字幕人妻熟女乱码| 婷婷色综合www| 巨乳人妻的诱惑在线观看| 久久中文字幕一级| 一级片'在线观看视频| 午夜激情av网站| 亚洲成人国产一区在线观看 | 精品人妻1区二区| 男女之事视频高清在线观看 | 久久久国产欧美日韩av| 天天躁夜夜躁狠狠久久av| 国产人伦9x9x在线观看| 中文字幕制服av| 免费高清在线观看日韩| 777米奇影视久久| 五月天丁香电影| 国产av一区二区精品久久| 免费高清在线观看视频在线观看| 91精品国产国语对白视频| 亚洲av在线观看美女高潮| 久久精品国产亚洲av高清一级| 丰满人妻熟妇乱又伦精品不卡| 婷婷色av中文字幕| 亚洲人成77777在线视频| 成年女人毛片免费观看观看9 | 国产成人精品久久二区二区91| √禁漫天堂资源中文www| 一边摸一边做爽爽视频免费| 欧美激情极品国产一区二区三区| 免费av中文字幕在线| 狂野欧美激情性bbbbbb| h视频一区二区三区| 国产福利在线免费观看视频| 免费看不卡的av| 久久人人爽人人片av| 久久久精品区二区三区| 一个人免费看片子| 国产成人一区二区在线| 黄频高清免费视频| 婷婷色综合大香蕉| 国产高清videossex| 超碰成人久久| 成年av动漫网址| 少妇被粗大的猛进出69影院| 日本黄色日本黄色录像| 久久人妻福利社区极品人妻图片 | 精品少妇内射三级| 日本欧美视频一区| 国产色视频综合| 成人影院久久| 午夜久久久在线观看| 国产黄色免费在线视频| 日本一区二区免费在线视频| 男女高潮啪啪啪动态图| 老熟女久久久| 丝袜人妻中文字幕| 欧美 亚洲 国产 日韩一| 成人午夜精彩视频在线观看| 亚洲国产av新网站| 下体分泌物呈黄色| 亚洲精品美女久久av网站| 操美女的视频在线观看| 免费看av在线观看网站| 天天操日日干夜夜撸| 黄片播放在线免费| 美女午夜性视频免费| 欧美精品一区二区免费开放| 9191精品国产免费久久| a级片在线免费高清观看视频| 欧美日韩福利视频一区二区| 亚洲欧美中文字幕日韩二区| 午夜福利在线免费观看网站| 国产高清国产精品国产三级| 日本91视频免费播放| 日韩av在线免费看完整版不卡| 国产成人av激情在线播放| 国产在线一区二区三区精| 大话2 男鬼变身卡| 爱豆传媒免费全集在线观看| 激情五月婷婷亚洲| 午夜福利影视在线免费观看| 999精品在线视频| 视频区欧美日本亚洲| 欧美日韩亚洲国产一区二区在线观看 | a级毛片黄视频| 欧美黑人欧美精品刺激| 午夜激情久久久久久久| 日韩制服丝袜自拍偷拍| 亚洲av欧美aⅴ国产| 久久久久久久精品精品| 在线观看一区二区三区激情| 欧美日韩福利视频一区二区| 99久久精品国产亚洲精品| 777久久人妻少妇嫩草av网站| 成年人免费黄色播放视频| 成人三级做爰电影| 黄网站色视频无遮挡免费观看| 最近中文字幕2019免费版| 日韩大片免费观看网站| 天天躁日日躁夜夜躁夜夜| 国产日韩欧美视频二区| 巨乳人妻的诱惑在线观看| 性色av乱码一区二区三区2| 精品久久久久久电影网| 老司机在亚洲福利影院| 香蕉丝袜av| 人人澡人人妻人| 国产在线一区二区三区精| 麻豆av在线久日| 国产成人a∨麻豆精品| h视频一区二区三区| 精品国产乱码久久久久久小说| 一边摸一边抽搐一进一出视频| 久久精品国产亚洲av涩爱| 观看av在线不卡| 亚洲国产中文字幕在线视频| 黄片播放在线免费| 国产日韩欧美视频二区| 国产成人精品无人区| 亚洲精品久久午夜乱码| 免费看不卡的av| 最新在线观看一区二区三区 | 人体艺术视频欧美日本| 日日夜夜操网爽| 欧美黑人精品巨大| 首页视频小说图片口味搜索 | 日韩大片免费观看网站| 久久人人97超碰香蕉20202| 男的添女的下面高潮视频| 国产1区2区3区精品| 国产精品国产av在线观看| 操美女的视频在线观看| 嫩草影视91久久| 一区二区三区精品91| av电影中文网址| 十八禁网站网址无遮挡| 日本猛色少妇xxxxx猛交久久| 嫩草影视91久久| 十八禁人妻一区二区| 免费日韩欧美在线观看| 天天添夜夜摸| 国精品久久久久久国模美| 欧美变态另类bdsm刘玥| 亚洲国产欧美在线一区| 久久久久久久大尺度免费视频| 男女边摸边吃奶| 欧美日韩国产mv在线观看视频| 亚洲成av片中文字幕在线观看| 国产亚洲欧美精品永久| 91成人精品电影| 制服诱惑二区| 久久久精品区二区三区| 免费看十八禁软件| a级毛片黄视频| 黑丝袜美女国产一区| 50天的宝宝边吃奶边哭怎么回事| 少妇被粗大的猛进出69影院| 少妇的丰满在线观看| 每晚都被弄得嗷嗷叫到高潮| 中文欧美无线码| 老司机影院成人| 性高湖久久久久久久久免费观看| xxxhd国产人妻xxx| 热99久久久久精品小说推荐| 久久久久国产一级毛片高清牌| 伊人久久大香线蕉亚洲五| 丝袜人妻中文字幕| 又紧又爽又黄一区二区| a级毛片在线看网站| av电影中文网址| 黄色视频不卡| 黄色a级毛片大全视频| 男女无遮挡免费网站观看| 亚洲av男天堂| 韩国精品一区二区三区| 亚洲欧洲国产日韩| 国产高清视频在线播放一区 | 亚洲国产精品一区三区| 国产av一区二区精品久久| 国产一区亚洲一区在线观看| 丝袜喷水一区| 日韩伦理黄色片| 日韩中文字幕欧美一区二区 | 老司机靠b影院| 国产成人影院久久av| 欧美日韩av久久| 性色av乱码一区二区三区2| 99久久综合免费| 亚洲精品美女久久av网站| 久久ye,这里只有精品| 久久狼人影院| 啦啦啦在线观看免费高清www| 男女午夜视频在线观看| 七月丁香在线播放| 新久久久久国产一级毛片| 18禁黄网站禁片午夜丰满| 久久久精品区二区三区| 久久九九热精品免费| 少妇 在线观看| 国产免费福利视频在线观看| 欧美 亚洲 国产 日韩一| 国产高清不卡午夜福利| 国产成人精品久久二区二区免费| 欧美精品高潮呻吟av久久| 51午夜福利影视在线观看| 国产一区二区在线观看av| 在线 av 中文字幕| 国产av精品麻豆| 高清不卡的av网站| 国产一区二区激情短视频 | 亚洲熟女毛片儿| 国产精品 国内视频| 亚洲精品日本国产第一区| 色网站视频免费| 亚洲,一卡二卡三卡| 国产伦理片在线播放av一区| 黄片播放在线免费| 一二三四社区在线视频社区8| 男女国产视频网站| 国产91精品成人一区二区三区 | 少妇的丰满在线观看| 女人高潮潮喷娇喘18禁视频| 极品人妻少妇av视频| 亚洲 国产 在线| 精品熟女少妇八av免费久了| 久久久久久久国产电影| 99热国产这里只有精品6| 妹子高潮喷水视频| 国产熟女午夜一区二区三区| 亚洲国产欧美日韩在线播放| 后天国语完整版免费观看| 麻豆乱淫一区二区| 十八禁人妻一区二区| 97精品久久久久久久久久精品| 亚洲一码二码三码区别大吗| 国产老妇伦熟女老妇高清| 免费观看人在逋| 亚洲成人手机| 欧美日韩亚洲综合一区二区三区_| 美女福利国产在线| 久久精品人人爽人人爽视色| 9191精品国产免费久久| 肉色欧美久久久久久久蜜桃| 五月天丁香电影| 最近最新中文字幕大全免费视频 | 欧美日韩成人在线一区二区| 热re99久久精品国产66热6| 嫁个100分男人电影在线观看 | 日本vs欧美在线观看视频| 老熟女久久久| 国产精品亚洲av一区麻豆| videos熟女内射| 两人在一起打扑克的视频| 欧美+亚洲+日韩+国产| 亚洲av成人精品一二三区| 中文字幕亚洲精品专区| 日本wwww免费看| 老熟女久久久| 爱豆传媒免费全集在线观看| 欧美在线黄色| 日本猛色少妇xxxxx猛交久久| 亚洲精品久久成人aⅴ小说| 视频区欧美日本亚洲| 亚洲中文日韩欧美视频| 久久女婷五月综合色啪小说| 欧美日韩国产mv在线观看视频| 九草在线视频观看| 搡老乐熟女国产| 亚洲精品国产av蜜桃| 制服诱惑二区| 秋霞在线观看毛片| 最近最新中文字幕大全免费视频 | 久久久久久久久久久久大奶| 一二三四在线观看免费中文在| 纵有疾风起免费观看全集完整版| 国产精品久久久久成人av| 久久精品国产亚洲av高清一级| 国产老妇伦熟女老妇高清| 纵有疾风起免费观看全集完整版| 精品少妇一区二区三区视频日本电影| a级毛片在线看网站| 亚洲欧美日韩另类电影网站| 后天国语完整版免费观看| 午夜视频精品福利| 精品国产一区二区久久| 久久久久久久国产电影| 亚洲激情五月婷婷啪啪| 99久久综合免费| 久久毛片免费看一区二区三区| 久久这里只有精品19| 午夜激情av网站| 亚洲国产精品一区二区三区在线| 亚洲欧美中文字幕日韩二区| 免费高清在线观看视频在线观看| 午夜视频精品福利| 七月丁香在线播放| h视频一区二区三区| 亚洲精品国产色婷婷电影| 免费看不卡的av| 久久综合国产亚洲精品| 成人免费观看视频高清| 亚洲精品国产av成人精品| 满18在线观看网站| 国产日韩欧美视频二区| 最黄视频免费看| 国产深夜福利视频在线观看| 久久久久国产一级毛片高清牌| 日韩精品免费视频一区二区三区| 国产日韩一区二区三区精品不卡| 国产在线免费精品| 天天躁日日躁夜夜躁夜夜| 久久人妻福利社区极品人妻图片 | 午夜福利视频在线观看免费| 欧美精品高潮呻吟av久久| 国产精品成人在线| 美国免费a级毛片| 99re6热这里在线精品视频| 18禁裸乳无遮挡动漫免费视频| 999精品在线视频| a级毛片在线看网站| 一级a爱视频在线免费观看| 在线观看免费日韩欧美大片| 老司机靠b影院| 国产一区二区在线观看av| 欧美精品高潮呻吟av久久| 国产亚洲午夜精品一区二区久久| 久久亚洲国产成人精品v| 久久99热这里只频精品6学生| 亚洲久久久国产精品| 欧美国产精品一级二级三级| 在线观看人妻少妇| 国产色视频综合| 美国免费a级毛片| 国产一区二区激情短视频 | 亚洲欧洲国产日韩| 久久精品久久久久久久性| 视频区欧美日本亚洲| 亚洲免费av在线视频| 欧美成人精品欧美一级黄| 18在线观看网站| 国产视频一区二区在线看| 香蕉国产在线看| 国产欧美日韩一区二区三区在线| 两人在一起打扑克的视频| 精品一区二区三卡| 丰满迷人的少妇在线观看| 人人妻,人人澡人人爽秒播 | 国产无遮挡羞羞视频在线观看| 婷婷色综合大香蕉| 天天躁狠狠躁夜夜躁狠狠躁| 免费久久久久久久精品成人欧美视频| 欧美人与善性xxx| 国产成人免费无遮挡视频| 国产有黄有色有爽视频| 欧美性长视频在线观看| 人体艺术视频欧美日本| 女人爽到高潮嗷嗷叫在线视频| 国产成人精品无人区| 九色亚洲精品在线播放| 交换朋友夫妻互换小说| 手机成人av网站| 国产精品国产三级专区第一集| 午夜精品国产一区二区电影| 18禁裸乳无遮挡动漫免费视频| 国产av一区二区精品久久| 久久99热这里只频精品6学生| 国产精品av久久久久免费| 成年人午夜在线观看视频| 18在线观看网站| 91字幕亚洲| 各种免费的搞黄视频| 蜜桃在线观看..| 国产av国产精品国产| 99国产综合亚洲精品| 曰老女人黄片| 免费人妻精品一区二区三区视频| 91字幕亚洲| 啦啦啦在线观看免费高清www| 欧美老熟妇乱子伦牲交| 午夜福利视频精品| 国产成人精品无人区| 免费看av在线观看网站| 国产精品亚洲av一区麻豆| 国产不卡av网站在线观看| 中文字幕高清在线视频| 纵有疾风起免费观看全集完整版| 亚洲国产欧美日韩在线播放| 国产精品人妻久久久影院| 热99久久久久精品小说推荐| 国产成人精品久久二区二区免费| 亚洲国产欧美一区二区综合| 飞空精品影院首页| 色综合欧美亚洲国产小说| 中文字幕人妻丝袜一区二区| 我要看黄色一级片免费的| 视频在线观看一区二区三区| 国产爽快片一区二区三区| 久久久久国产一级毛片高清牌| 欧美黄色片欧美黄色片| 一区二区三区乱码不卡18| 麻豆av在线久日| 欧美亚洲日本最大视频资源| 亚洲av日韩在线播放| 国产精品 国内视频| 手机成人av网站| a级毛片在线看网站| 777久久人妻少妇嫩草av网站| 亚洲七黄色美女视频| 久久精品亚洲熟妇少妇任你| 丝袜人妻中文字幕| 人人妻,人人澡人人爽秒播 | 性色av一级| 日韩大片免费观看网站| 人妻一区二区av| 中文字幕精品免费在线观看视频| 亚洲欧美一区二区三区黑人| 五月天丁香电影| 亚洲五月色婷婷综合| 男女下面插进去视频免费观看| 欧美性长视频在线观看| 日本a在线网址| 久久久久网色| 精品免费久久久久久久清纯 | 51午夜福利影视在线观看| 免费看av在线观看网站| 91麻豆av在线| 国产在线视频一区二区| 亚洲欧美中文字幕日韩二区| 日韩伦理黄色片| 少妇精品久久久久久久| 成人18禁高潮啪啪吃奶动态图| 最近中文字幕2019免费版| 亚洲三区欧美一区| 日韩大片免费观看网站| 18禁国产床啪视频网站| 一本大道久久a久久精品| 欧美成狂野欧美在线观看| 午夜影院在线不卡| 亚洲欧美日韩高清在线视频 | 亚洲成色77777| 下体分泌物呈黄色| 日本wwww免费看| 91麻豆精品激情在线观看国产 | 丝袜喷水一区| 精品少妇久久久久久888优播| 亚洲国产精品一区二区三区在线| 午夜视频精品福利| 黑人猛操日本美女一级片| 国产高清不卡午夜福利| 在线观看www视频免费| 国产精品久久久久久精品古装| 国产麻豆69| 18禁国产床啪视频网站| 免费在线观看视频国产中文字幕亚洲 | 欧美日韩福利视频一区二区| 久热这里只有精品99| netflix在线观看网站| 在线观看免费日韩欧美大片| 亚洲第一av免费看| 国产1区2区3区精品| 免费看十八禁软件| 丝袜喷水一区| av国产精品久久久久影院| 永久免费av网站大全| 午夜精品国产一区二区电影| 久久久久精品国产欧美久久久 | 欧美精品人与动牲交sv欧美| 乱人伦中国视频| 精品一区在线观看国产| 久久久国产精品麻豆| 一级黄色大片毛片| 国产成人欧美在线观看 | 亚洲午夜精品一区,二区,三区| 亚洲av成人精品一二三区| 亚洲成人国产一区在线观看 | 亚洲美女黄色视频免费看| 搡老岳熟女国产| av在线播放精品| 亚洲精品自拍成人| 亚洲国产毛片av蜜桃av| 亚洲 欧美一区二区三区| 麻豆av在线久日| 免费黄频网站在线观看国产| av国产精品久久久久影院| 一级,二级,三级黄色视频| 丝袜喷水一区| 免费一级毛片在线播放高清视频 | 国产成人精品久久二区二区91| 好男人电影高清在线观看| 亚洲av电影在线观看一区二区三区| 香蕉丝袜av| 亚洲精品国产av蜜桃| 久久午夜综合久久蜜桃| 最近最新中文字幕大全免费视频 | svipshipincom国产片| 美女脱内裤让男人舔精品视频| 啦啦啦啦在线视频资源| 日韩av免费高清视频| 99国产综合亚洲精品| 飞空精品影院首页| 伦理电影免费视频| 亚洲国产精品一区二区三区在线| 国产深夜福利视频在线观看| 黄色a级毛片大全视频| 日韩人妻精品一区2区三区| 国产日韩欧美亚洲二区|