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

    城市建設(shè)用地海綿化程度對地表徑流水質(zhì)的影響
    ——以浙江省嘉興市為例

    2019-10-11 04:38:22陳前虎李玉蓮黃初冬王賢萍
    水土保持通報 2019年4期
    關(guān)鍵詞:不透水徑流綠地

    陳前虎, 李玉蓮, 黃初冬, 王賢萍

    (1.浙江工業(yè)大學(xué) 建筑工程學(xué)院, 浙江 杭州 310023; 2.嘉興市規(guī)劃設(shè)計研究院, 浙江 嘉興 314000)

    為解決城市化進程中因地表硬鋪率[1](impervious surface ratio,ISR)增加帶來的徑流水質(zhì)惡化問題,我國提出了基于“滲、滯、蓄、凈、用、排”為導(dǎo)向的海綿城市建設(shè)策略,出臺了一系列有關(guān)的標(biāo)準(zhǔn)規(guī)范與技術(shù)導(dǎo)則[2-6]。這些文件均強調(diào)利用海綿設(shè)施對雨水徑流的截留作用來控制徑流污染負(fù)荷,并將年徑流總量控制率作為控制徑流污染物的重要指標(biāo),而年徑流總量控制與地塊的海綿化程度——海綿設(shè)施的類型、數(shù)量與布局等息息相關(guān)[7-10]。為此,在海綿城市的規(guī)劃設(shè)計實踐中,需要探究的兩大科學(xué)問題主要是哪些因素對徑流污染削減起關(guān)鍵性作用;為達(dá)到一定的水質(zhì)控制目標(biāo),該如何處理這些因素間的量化結(jié)構(gòu)關(guān)系。在各地加快實施海綿城市建設(shè)的背景下,有關(guān)建設(shè)用地海綿化程度與地表徑流水質(zhì)響應(yīng)關(guān)系的研究迫在眉睫。圍繞上述問題,國內(nèi)外學(xué)者從時間和空間兩維視角出發(fā),采用多元統(tǒng)計方法或借助雨洪管理模型,分析和模擬不同用地類型和降雨強度下典型海綿設(shè)施對徑流污染的影響[11-16]。劉文等[17-18]研究提出,海綿設(shè)施對污染物濃度的控制具有長期顯著性,不同類型的設(shè)施對污染物濃度的削減程度差異明顯;胡愛兵等[19-20]研究表明,優(yōu)化海綿設(shè)施的布局結(jié)構(gòu)有利于改善用地徑流水質(zhì);Roy等[21-22]發(fā)現(xiàn),低影響開發(fā)主要通過改變總不透水下墊面(total impervious area, TIA)與排水收集系統(tǒng)的連通性來控制對雨水徑流及污染負(fù)荷的攔截程度,且TIA分為有效不透水下墊面(effective impervious area, EIA,指與排水管渠直接連接的不透水下墊面)和無效不透水下墊面(ineffective impervious area, IIA,指與排水管渠斷接的不透水下墊面),其中EIA是導(dǎo)致雨水中污染負(fù)荷增加的主要因素;Sohn等[23-25]研究認(rèn)為,EIA是評價LID設(shè)施截留雨水徑流能力的重要參考指標(biāo)。這些研究表明海綿化程度對徑流水質(zhì)有重要影響,為從土地利用視角控制徑流污染提供了理論基礎(chǔ)和新的視角與方法。但如何量化、揭示并管控海綿化程度與徑流水質(zhì)的響應(yīng)關(guān)系,在理論和實踐中均具有重大意義,值得進一步探索?;诖?,本文以首批海綿城市試點——嘉興市的20個海綿化改造不同程度的項目為研究單元,依托定點測量的數(shù)據(jù),建立海綿化與地表徑流水質(zhì)兩套指標(biāo);采用冗余分析、偏最小二乘法及Origin擬合方程,研究海綿化程度對徑流水質(zhì)的影響,明確對徑流水質(zhì)顯著作用的海綿化指標(biāo),并嘗試探究在水質(zhì)達(dá)到地表水Ⅳ類及污水排放二級標(biāo)準(zhǔn)時,改造程度不同的用地需要管控的指標(biāo)閾值,為海綿城市建設(shè)的控制性詳細(xì)規(guī)劃提供技術(shù)指導(dǎo)依據(jù)。

    1 研究區(qū)概況

    嘉興市地處長三角沖積平原地帶,地勢平坦且水系豐富,是典型的水質(zhì)型缺水城市。隨著嘉興市污水治理力度的推進,點源污染問題逐步解決,以徑流污染為主的非點源污染對河道水質(zhì)的影響日益增大。研究選取嘉興市10個已海綿化改造半年以上(多位于海綿城市改造示范區(qū)內(nèi))和10個未海綿化改造(多位于示范區(qū)外)共20個單一性質(zhì)的建設(shè)用地為樣本(圖1)。為保證各樣本排水管網(wǎng)系統(tǒng)的獨立性以及用地面積的均等性,將府南三期小區(qū)以戚家村浜為界分為東區(qū)和西區(qū),中央公園和植物園分別以由拳路、戚家浜為界分為南區(qū)和北區(qū)。樣本用地類型有居住小區(qū)(10個)、公共建筑(6個)和公園(4個)3類,各用地容積率均在1.5及以下。

    圖1 研究區(qū)建設(shè)用地樣本點分布示意圖

    2 數(shù)據(jù)獲取與研究方法

    2.1 海綿化指標(biāo)的獲取

    結(jié)合用地改造現(xiàn)狀,所有研究區(qū)地面(含屋面)的類型可用EIA、IIA、綠地(含各種LID設(shè)施)以及弱透水面(透水路面及水體)4種類型來涵蓋。依據(jù)下墊面的透水性強度、不透水下墊面與排水管網(wǎng)的斷接或連接程度,選取EIA比例、IIA比例、綠地比例以及弱透水比例作為衡量海綿化程度的指標(biāo)。利用Excel軟件及Sutherland經(jīng)驗公式[26]計算出各樣本用地具體數(shù)值,其中Sutherland經(jīng)驗公式主要用于計算EIA比例,該公式解決了美國地質(zhì)勘探局提出的僅適用于小流域(8~28 hm2)內(nèi)EIA值計算方法的局限性問題,具體方式參考表1。結(jié)合各研究區(qū)海綿設(shè)施的建設(shè)施工現(xiàn)狀及表1可知,研究樣本中已改造居住小區(qū)和公共建筑的不透水區(qū)域與排水收集系統(tǒng)部分?jǐn)嘟?,公園用地的不透水區(qū)域與排水收集系統(tǒng)高度斷接;未改造用地的不透水區(qū)域與排水收集系統(tǒng)完全連接。

    表1 依據(jù)Sutherland方程式求得的EIA[23,26]

    2.2 水質(zhì)數(shù)據(jù)的采集

    水質(zhì)數(shù)據(jù)采用人工時間間隔(每隔10 min)采樣法,選取2017年3—8月10場平均降雨強度等級為中雨的有效降雨,并在每塊建設(shè)用地選取2~3個連接河道的地下管網(wǎng)雨水排出口作為水質(zhì)采樣點(共45個),由此獲取1391份徑流水質(zhì)數(shù)據(jù)。研究選擇懸浮物(SS)、生化需氧量(BOD5)、化學(xué)需氧量(COD)、氨氮(NH3-N)、總磷(TP)為水質(zhì)指標(biāo),并利用Grubbs檢驗法[27]篩選數(shù)據(jù)并剔除離群值,最終求得多場降雨徑流中各指標(biāo)的平均濃度[28]。

    2.3 研究方法

    采用冗余分析定性判斷海綿化程度與徑流水質(zhì)的相關(guān)性以及海綿化指標(biāo)對水質(zhì)指標(biāo)的控制作用方向,運用偏最小二乘法定量分析海綿化指標(biāo)對徑流水質(zhì)的影響權(quán)重、篩選出影響水質(zhì)的關(guān)鍵指標(biāo),利用Origin擬合方程預(yù)測改造程度不同的用地徑流水質(zhì)達(dá)到目標(biāo)要求時,海綿化指標(biāo)的響應(yīng)閾值。

    2.3.1 冗余分析 借助Canoco 5.0軟件對各樣本點不同水質(zhì)指標(biāo)進行降趨勢對應(yīng)分析可知,水質(zhì)指標(biāo)數(shù)據(jù)第一軸的梯度值為0.171,遠(yuǎn)小于3,故選擇冗余分析[29]。冗余分析是一種基于排序技術(shù)的線性分析方法,通過樣本點在象限的分布形態(tài)等方式來揭示研究區(qū)域的特點。

    2.3.2 偏最小二乘法 偏最小二乘法是一種多因變量對多自變量的回歸建模方法,通過提取自變量和因變量的主成分,篩選出對兩變量解釋性最強的綜合變量,并解決了變量間的多重共線性問題[30]。

    2.3.3 Origin擬合方程 Origin軟件[31]自帶多種非線性擬合函數(shù),能方便、準(zhǔn)確地擬合出兩組變量數(shù)據(jù)間的最佳關(guān)系曲線,并可直觀地了解變量的變化趨勢及偏差。

    3 結(jié)果與分析

    3.1 海綿化指標(biāo)特征

    研究表明,樣本用地的下墊面類型以EIA為主(公園用地除外),已改造建設(shè)用地的下墊面類型較未改造建設(shè)用地豐富(圖2)。已改造樣本R01—R06的下墊面類型以EIA和綠地為主,R07—R10的下墊面類型以綠地為主;未改造樣本U01—U10的下墊面類型主要為EIA,綠地次之。不同用地類型的綠地比例和透水比例(綠地和弱透水比例之和)差異明顯,其中綠地比例由高到低依次為:公園(超過70%)>居住小區(qū)(超過30%)>公共建筑(除市政府U05外均低于25%);透水比例由高到低依次為:公園(超過95%)>居住小區(qū)(超過52%)>公共建筑(低于45%)。公共建筑中汽車站的兩項指標(biāo)值最低(均為11%)。當(dāng)綠地比例相同時,已改造建設(shè)用地EIA比例(低于45%)均低于未改造建設(shè)用地(超過60%)。

    3.2 水質(zhì)指標(biāo)特征

    由圖3可知,樣本中部分已改造用地各水質(zhì)指標(biāo)的平均濃度未滿足地表水Ⅳ類及污水二級排放標(biāo)準(zhǔn),所有未改造用地各水質(zhì)指標(biāo)的平均濃度均未達(dá)到水質(zhì)要求。三種用地類型中公園的水質(zhì)最好,居住小區(qū)次之,公共建筑的水質(zhì)明顯劣于前兩者。公園的綠地比例和透水比例較高、EIA比例較低,海綿設(shè)施布局較好;而公共建筑的綠地和透水比例較低、EIA比例較高,海綿設(shè)施布局較差。20個樣本中,嘉興汽車北站(U10)的水質(zhì)最差,SS,BOD5等五項指標(biāo)的平均濃度明顯高于其他樣本。對比樣本R01—R06與U01—U06并結(jié)合圖2發(fā)現(xiàn),在相同的綠地比例下,已改造建設(shè)用地的EIA比例低于未改造建設(shè)用地,而水質(zhì)則前者更優(yōu);隨著綠地比例的上升,建設(shè)用地各項水質(zhì)指標(biāo)的平均濃度逐漸下降(R06與U06除外)。

    圖2 研究區(qū)各樣本用地不同海綿化指標(biāo)的比例

    3.3 海綿化程度與徑流水質(zhì)響應(yīng)分析

    3.3.1 海綿化指標(biāo)與水質(zhì)指標(biāo)的關(guān)系 借助Canoc 5.0軟件對樣本數(shù)據(jù)進行冗余分析,得出海綿化程度與水質(zhì)的關(guān)系圖(圖4)。圖中,箭頭長度及箭頭間的夾角表示水質(zhì)指標(biāo)對各海綿化指標(biāo)的響應(yīng)程度,長度越長、夾角(銳角表示海綿化指標(biāo)與水質(zhì)指標(biāo)呈正相關(guān),鈍角表示負(fù)相關(guān))的余弦絕對值越大,說明二者相關(guān)性越大,反之越小。

    由圖4可知,影響各樣本徑流水質(zhì)的海綿化指標(biāo)類型差異顯著。已改造公共建筑類樣本(R01—R06)主要受IIA比例影響,公園類(R08,R09)和(R07,R10)分別受綠地和弱透水比例影響明顯;未改造樣本(U01—U10)則主要受綠地和EIA比例影響。EIA比例與各水質(zhì)指標(biāo)呈顯著正相關(guān),表明建設(shè)用地內(nèi)EIA是導(dǎo)致徑流污染負(fù)荷增加的主要因素;綠地、弱透水和IIA三者比例與各項水質(zhì)指標(biāo)呈負(fù)相關(guān)關(guān)系,且綠地比例與水質(zhì)指標(biāo)的相關(guān)性最大,表明綠地對徑流污染負(fù)荷的截留程度高于另外兩者。SS,BOD5,COD和NH3-N這4項水質(zhì)指標(biāo)受海綿化程度影響明顯,TP受海綿化程度影響不顯著。對比EIA比例、IIA比例與水質(zhì)指標(biāo)的相關(guān)性后發(fā)現(xiàn),二者對各水質(zhì)指標(biāo)的作用方向互為相反,且EIA比例與各水質(zhì)指標(biāo)的相關(guān)性較IIA比例高。

    圖3 研究區(qū)不同建設(shè)用地水質(zhì)指標(biāo)空間變化

    圖4 海綿化指標(biāo)與水質(zhì)指標(biāo)的冗余分析

    3.3.2 海綿化與水質(zhì)關(guān)鍵指標(biāo)的選取 本文運用SIMCA-P軟件構(gòu)建偏最小二乘回歸模型,并在Analysis中點擊Autofit功能對模型進行自動擬合。利用軟件的交叉有效性原則提取出2個PLS成分,模型對X的信息利用率為R2X(cum)=0.942,對Y的解釋能力為R2Y(cum)=0.828,且求得的交叉有效性值Q2(cum)=0.781 9(大于0.5),表示模型的精度較高。最終,得到經(jīng)過標(biāo)準(zhǔn)化的偏最小二乘回歸方程為:

    (1)

    式中:f(x)——水質(zhì)指標(biāo)濃度(mg/L);αρ——第p種海綿化指標(biāo)對水質(zhì)指標(biāo)的影響常數(shù)項;μρ——第p種海綿化指標(biāo)的影響系數(shù);βρ——第p種海綿化指標(biāo)(%)。

    不同海綿化指標(biāo)對水質(zhì)指標(biāo)的影響系數(shù)如表2所示,影響大小順序為:綠地比例>EIA比例>IIA比例>弱透水比例。其中,SS受四項海綿化指標(biāo)影響明顯;BOD5和COD主要受綠地和EIA比例影響;海綿化指標(biāo)對NH3-N的影響不如其他水質(zhì)指標(biāo)顯著,說明其不是導(dǎo)致NH3-N平均濃度變化的最重要因素。綜合來看,EIA比例越大,各污染物的平均濃度越高,水質(zhì)越差;綠地、IIA和弱透水比例越大,各污染物的平均濃度越低,水質(zhì)越好。降雨時,綠地及透水路面等透水區(qū)域的雨水經(jīng)過滲透進入地下,不易形成地表徑流,污染物負(fù)荷較難通過排水管網(wǎng)進入河道;IIA雖為不透水區(qū)域,但與排水管道斷接,其地表徑流主要流向附近的海綿設(shè)施,未能直接進入排水管道。因此,提高建設(shè)用地內(nèi)綠地及透水路面比例或?qū)⒉煌杆访娴膹搅饕刖G地等海綿設(shè)施、降低EIA比例,都將有利于控制用地內(nèi)徑流污染從而改善河道水質(zhì)。

    表2 海綿化指標(biāo)對水質(zhì)指標(biāo)的影響系數(shù)

    海綿化指標(biāo)對水質(zhì)指標(biāo)影響的重要程度,可通過變量投影重要性指標(biāo)[32](variable importance in the projection, VIP)表示(圖5)。根據(jù)VIP排序圖可知,綠地和EIA比例的VIP值大于1,IIA比例的VIP值介于0.8,1之間,而弱透水比例的VIP值小于0.8,表明綠地和EIA比例對水質(zhì)指標(biāo)具有顯著解釋意義,二者為影響水質(zhì)的關(guān)鍵性指標(biāo)。

    圖5 影響水質(zhì)的海綿化指標(biāo)VIP值排序

    為了驗證模型建立的優(yōu)劣度和可靠性,需對模型進行預(yù)測分析。利用偏最小二乘回歸方程求得各水質(zhì)指標(biāo)的預(yù)測值,結(jié)合實測值繪制的散點趨勢圖分析可知,擬合度較好(R2>0.8)的水質(zhì)指標(biāo)為SS,BOD5和COD,樣本點多數(shù)均勻分布在對角線附近,表明SS,BOD5和COD是受海綿化程度影響的主要指標(biāo)。

    3.3.3 海綿化指標(biāo)的閾值探討 綜上分析,以EIA比例為x軸,以SS,BOD5及COD三項指標(biāo)的平均濃度為y軸,根據(jù)樣本中不透水區(qū)域與排水管網(wǎng)系統(tǒng)斷接或連接程度的不同,利用Origin 2017軟件建立9個擬合方程;再以綠地比例為x軸,與其他三項水質(zhì)指標(biāo)建立3個擬合方程。根據(jù)研究獲取的樣本數(shù)據(jù),可預(yù)測中度降水下,不同改造程度的用地徑流水質(zhì)達(dá)到地表水Ⅳ類及污水排放二級標(biāo)準(zhǔn)時,需要管控的海綿化指標(biāo)的范圍。

    (1) 當(dāng)EIA≈0時,即針對不透水區(qū)域與排水收集系統(tǒng)高度或完全斷接的建設(shè)用地,隨著綠地比例的上升,各水質(zhì)指標(biāo)的平均濃度均相應(yīng)下降。由表3方程可知,若使建設(shè)用地的水質(zhì)滿足標(biāo)準(zhǔn),綠地比例應(yīng)不低于31%。

    (2) 當(dāng)EIA=0.04(TIA)1.7時,即針對不透水區(qū)域與排水收集系統(tǒng)部分?jǐn)嘟拥慕ㄔO(shè)用地,各項水質(zhì)指標(biāo)的平均濃度隨著綠地比例的上升而降低,且斜率逐漸減小,其中COD對綠地比例變化的敏感性最高。由圖6可知,綠地比例超過38.5%時,建設(shè)用地的水質(zhì)可滿足標(biāo)準(zhǔn)。由此可見,目前嘉興市已改造居住小區(qū)和公共建筑的綠地比例均有待提高。

    (3) 當(dāng)EIA=TIA時,即針對未經(jīng)過海綿化改造的建設(shè)用地,隨著綠地比例的上升,各水質(zhì)指標(biāo)的平均濃度均下降,且SS的下降斜率逐漸增大,說明SS對綠地比例變化的敏感性更高。由圖7可知,若使建設(shè)用地的徑流水質(zhì)滿足要求,綠地比例應(yīng)不低于47%。目前看來,嘉興市未改造居住小區(qū)和公共建筑的綠地比例均未超過47%,徑流污染較為嚴(yán)重。

    (4) 在常規(guī)綠地比例建設(shè)標(biāo)準(zhǔn)(綠地比例≈35%)下,隨著EIA比例的上升,各水質(zhì)指標(biāo)的平均濃度均上升,且斜率由低變高。由圖8可知,EIA比例低于9%時,水質(zhì)可滿足目標(biāo)要求。然而,目前嘉興市除公園用地外其他類型建設(shè)用地的EIA比例均超過9%。

    圖8 常規(guī)綠地比例下,SS,BOD5和COD與EIA比例的回歸預(yù)測

    由分析結(jié)果可知,在特定的海綿化程度下,當(dāng)改變綠地和EIA比例使COD滿足水質(zhì)要求時,其余水質(zhì)指標(biāo)均在標(biāo)準(zhǔn)范圍之內(nèi)。因此,COD是評價研究樣本用地徑流水質(zhì)的關(guān)鍵性指標(biāo)。

    4 結(jié) 論

    (1) 基于冗余分析,定性地揭示了海綿化程度與徑流水質(zhì)的響應(yīng)關(guān)系。影響各樣本用地徑流水質(zhì)的海綿化指標(biāo)不同;綠地、IIA、弱透水比例與水質(zhì)指標(biāo)呈負(fù)相關(guān),EIA比例與水質(zhì)指標(biāo)呈正相關(guān);SS,BOD5,COD和NH3-N與海綿化指標(biāo)相關(guān)性較大。

    (2) 基于偏最小二乘回歸分析,定量地探究出不同海綿化指標(biāo)對各項水質(zhì)指標(biāo)的影響權(quán)重,篩選得到影響水質(zhì)的關(guān)鍵指標(biāo)為綠地比例和EIA比例,受海綿化程度影響的主要水質(zhì)指標(biāo)為SS,BOD5和COD。

    (3) 基于關(guān)鍵指標(biāo)的Origin擬合分析,得出用地徑流水質(zhì)達(dá)到Ⅳ類及污水排放二級標(biāo)準(zhǔn)時綠地或EIA比例的響應(yīng)閾值: ①當(dāng)EIA≈0時,綠地比例的閾值為31%; ②當(dāng)EIA=0.04(TIA)1.7時,綠地比例的閾值為38.5%; ③當(dāng)EIA=TIA時,綠地比例的閾值為47%; ④在常規(guī)綠地比例(約為35%)建設(shè)標(biāo)準(zhǔn)下,EIA比例的閾值為9%。

    當(dāng)前,《城市綠化規(guī)劃建設(shè)指標(biāo)的規(guī)定》中規(guī)定城市建設(shè)用地綠地比例應(yīng)不低于35%,其中,居住小區(qū)和公共建筑的綠地比例標(biāo)準(zhǔn)分別為40%(新小區(qū))、35%(舊小區(qū)和公共建筑)。由此可知,當(dāng)EIA≈0時,則無需改變城市建設(shè)用地內(nèi)綠地比例標(biāo)準(zhǔn);當(dāng)EIA=0.04(TIA)1.7時,建議已改造的居住小區(qū)和公共建筑的綠地比例標(biāo)準(zhǔn)均提高3.5%;當(dāng)EIA=TIA時,建議未改造建設(shè)用地中居住小區(qū)和公共建筑的綠地比例分別提高7%(新小區(qū))和12%(舊小區(qū)和公共建筑);在常規(guī)綠地比例(約為35%)建設(shè)標(biāo)準(zhǔn)下,可通過增加透水路面面積或?qū)⑵胀ňG地改為下凹式綠地、雨水花園及雨落管等方式降低EIA比例。

    致謝:嘉興市規(guī)劃設(shè)計研究院為本研究提供了詳細(xì)的海綿城市建設(shè)施工圖和降雨數(shù)據(jù)資料,特此致謝。

    猜你喜歡
    不透水徑流綠地
    基于無人機可見光影像與OBIA-RF算法的城市不透水面提取
    可進入式綠地的養(yǎng)護與管理的研究
    北京園林(2021年2期)2022-01-19 02:17:34
    老城區(qū)綠地更新策略——以北京西城區(qū)綠地提升為例
    北京園林(2020年1期)2020-01-18 05:19:30
    Landsat8不透水面遙感信息提取方法對比
    杭州市區(qū)春季綠地對PM2.5消減作用的研究
    走進綠地
    歲月(2017年9期)2017-09-08 14:45:50
    Topmodel在布哈河流域徑流模擬中的應(yīng)用
    探秘“大徑流”
    攻克“大徑流”
    江埡水庫降雨徑流相關(guān)圖的建立
    91字幕亚洲| 一本久久精品| 在线观看免费高清a一片| 男男h啪啪无遮挡| 亚洲中文日韩欧美视频| 日韩免费高清中文字幕av| 欧美人与性动交α欧美软件| 91精品国产国语对白视频| 脱女人内裤的视频| 9191精品国产免费久久| 久久99热这里只频精品6学生| 黑人操中国人逼视频| 91九色精品人成在线观看| 热99国产精品久久久久久7| 精品少妇一区二区三区视频日本电影| 欧美日韩黄片免| 亚洲全国av大片| videos熟女内射| 国产淫语在线视频| 在线天堂中文资源库| 中文字幕另类日韩欧美亚洲嫩草| 满18在线观看网站| 两性夫妻黄色片| 欧美日韩一级在线毛片| av片东京热男人的天堂| 免费久久久久久久精品成人欧美视频| 色视频在线一区二区三区| 一区二区三区国产精品乱码| 久久精品亚洲熟妇少妇任你| 国产亚洲一区二区精品| 亚洲专区中文字幕在线| 少妇精品久久久久久久| 国产主播在线观看一区二区| 国产精品九九99| 亚洲精品中文字幕一二三四区 | 精品午夜福利视频在线观看一区 | 国产精品亚洲一级av第二区| 欧美日韩黄片免| 咕卡用的链子| 国产亚洲欧美精品永久| 欧美人与性动交α欧美精品济南到| 黄频高清免费视频| 国产一卡二卡三卡精品| 国产高清激情床上av| 亚洲精品中文字幕在线视频| 成年版毛片免费区| 国产在视频线精品| 亚洲中文日韩欧美视频| 91老司机精品| 国产97色在线日韩免费| 高清av免费在线| 黄色a级毛片大全视频| 丝袜美足系列| 亚洲午夜精品一区,二区,三区| 老汉色∧v一级毛片| 亚洲第一av免费看| 国产成人系列免费观看| 成年动漫av网址| 亚洲精品一二三| 视频区欧美日本亚洲| 天堂中文最新版在线下载| 国产福利在线免费观看视频| 波多野结衣一区麻豆| 18禁裸乳无遮挡动漫免费视频| 中文字幕色久视频| 自拍欧美九色日韩亚洲蝌蚪91| 搡老乐熟女国产| 高清av免费在线| 亚洲成国产人片在线观看| 国产精品偷伦视频观看了| 免费在线观看黄色视频的| 国产亚洲欧美精品永久| 91大片在线观看| 9191精品国产免费久久| 亚洲成人免费电影在线观看| videos熟女内射| 99国产精品免费福利视频| 亚洲精品中文字幕一二三四区 | 我要看黄色一级片免费的| 欧美日韩精品网址| 成人国产一区最新在线观看| 久久亚洲精品不卡| h视频一区二区三区| 欧美黑人欧美精品刺激| 亚洲成av片中文字幕在线观看| 国内毛片毛片毛片毛片毛片| 亚洲精品自拍成人| 99热网站在线观看| 黄色a级毛片大全视频| 老熟妇乱子伦视频在线观看| 999久久久国产精品视频| 国产精品香港三级国产av潘金莲| 91成年电影在线观看| 汤姆久久久久久久影院中文字幕| 亚洲人成伊人成综合网2020| 高清av免费在线| 高潮久久久久久久久久久不卡| 亚洲国产欧美网| 久久精品亚洲熟妇少妇任你| 一级黄色大片毛片| 999精品在线视频| 色老头精品视频在线观看| 成人影院久久| 久久国产精品人妻蜜桃| 亚洲精品乱久久久久久| av又黄又爽大尺度在线免费看| 久久久久视频综合| 亚洲精品自拍成人| 精品一区二区三卡| 成人av一区二区三区在线看| 一边摸一边抽搐一进一出视频| netflix在线观看网站| 美女福利国产在线| 免费在线观看日本一区| 黄片播放在线免费| 免费不卡黄色视频| 天堂中文最新版在线下载| 黄频高清免费视频| 久久午夜亚洲精品久久| 午夜老司机福利片| 亚洲精品国产一区二区精华液| 国产成人精品久久二区二区免费| 亚洲av片天天在线观看| 欧美 亚洲 国产 日韩一| 亚洲国产欧美网| 成年版毛片免费区| 免费在线观看影片大全网站| 日本五十路高清| 国产在线视频一区二区| 满18在线观看网站| 三级毛片av免费| 高潮久久久久久久久久久不卡| 国产精品99久久99久久久不卡| 午夜福利视频在线观看免费| 国产一区二区三区在线臀色熟女 | 亚洲精品国产区一区二| 成年动漫av网址| 欧美日韩亚洲综合一区二区三区_| 后天国语完整版免费观看| 久久人妻av系列| 色婷婷久久久亚洲欧美| 人人妻,人人澡人人爽秒播| 十八禁网站网址无遮挡| 一边摸一边抽搐一进一小说 | 久久久久久久国产电影| 久久av网站| 久久亚洲真实| 自拍欧美九色日韩亚洲蝌蚪91| 国产成人精品久久二区二区免费| 亚洲国产成人一精品久久久| 99精品欧美一区二区三区四区| a级毛片在线看网站| av又黄又爽大尺度在线免费看| 午夜91福利影院| 丝瓜视频免费看黄片| 亚洲第一青青草原| 日韩人妻精品一区2区三区| 精品久久蜜臀av无| av电影中文网址| 纯流量卡能插随身wifi吗| 18禁国产床啪视频网站| 中亚洲国语对白在线视频| 精品熟女少妇八av免费久了| 国产成人啪精品午夜网站| 五月天丁香电影| 美女视频免费永久观看网站| 一区福利在线观看| 久久精品aⅴ一区二区三区四区| 亚洲精品一二三| 日韩人妻精品一区2区三区| e午夜精品久久久久久久| 一级毛片精品| 黑人巨大精品欧美一区二区mp4| 成年动漫av网址| 亚洲欧美一区二区三区黑人| 日韩大片免费观看网站| 精品少妇内射三级| 少妇粗大呻吟视频| www.999成人在线观看| 欧美成狂野欧美在线观看| 中文字幕最新亚洲高清| 麻豆成人av在线观看| av在线播放免费不卡| 免费女性裸体啪啪无遮挡网站| 又大又爽又粗| 91麻豆av在线| www日本在线高清视频| 国产在线观看jvid| 男女之事视频高清在线观看| 久久精品国产99精品国产亚洲性色 | 欧美成人午夜精品| 涩涩av久久男人的天堂| 精品国产一区二区三区久久久樱花| aaaaa片日本免费| 国产成人精品久久二区二区免费| 亚洲专区中文字幕在线| 美女高潮到喷水免费观看| 国产熟女午夜一区二区三区| 多毛熟女@视频| 亚洲一区中文字幕在线| av在线播放免费不卡| 精品福利永久在线观看| 黑人猛操日本美女一级片| 老熟妇乱子伦视频在线观看| 一区二区日韩欧美中文字幕| 欧美在线一区亚洲| 欧美中文综合在线视频| 一边摸一边抽搐一进一出视频| 亚洲欧美色中文字幕在线| 久久国产精品大桥未久av| 国产色视频综合| 我的亚洲天堂| 激情视频va一区二区三区| 高清视频免费观看一区二区| 国产精品偷伦视频观看了| 久久人妻av系列| 另类精品久久| 天堂动漫精品| 美女主播在线视频| 丁香六月欧美| 乱人伦中国视频| 大香蕉久久网| 色精品久久人妻99蜜桃| aaaaa片日本免费| 免费一级毛片在线播放高清视频 | 老司机深夜福利视频在线观看| 一边摸一边做爽爽视频免费| 91字幕亚洲| 国产亚洲精品第一综合不卡| 国产免费av片在线观看野外av| 亚洲熟妇熟女久久| 满18在线观看网站| 国产成人免费观看mmmm| 一区二区三区国产精品乱码| 久久婷婷成人综合色麻豆| 欧美精品高潮呻吟av久久| 久久久久精品人妻al黑| 男女免费视频国产| 中文欧美无线码| 久久久国产一区二区| 天天躁日日躁夜夜躁夜夜| 国产日韩欧美亚洲二区| 国产免费福利视频在线观看| 欧美另类亚洲清纯唯美| 侵犯人妻中文字幕一二三四区| 日本欧美视频一区| 色94色欧美一区二区| 亚洲免费av在线视频| 超碰成人久久| 另类亚洲欧美激情| 99riav亚洲国产免费| 亚洲精品久久午夜乱码| 丰满饥渴人妻一区二区三| 久久久久精品人妻al黑| 国产成人av教育| 成人永久免费在线观看视频 | 亚洲黑人精品在线| 大片电影免费在线观看免费| 少妇 在线观看| 久久久久精品人妻al黑| 国产成人av教育| 亚洲va日本ⅴa欧美va伊人久久| 精品视频人人做人人爽| 亚洲av国产av综合av卡| 淫妇啪啪啪对白视频| 纯流量卡能插随身wifi吗| 丝袜美腿诱惑在线| 欧美黄色淫秽网站| 一边摸一边抽搐一进一出视频| 18禁国产床啪视频网站| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩国产mv在线观看视频| 巨乳人妻的诱惑在线观看| av超薄肉色丝袜交足视频| 精品少妇久久久久久888优播| 久久久久久人人人人人| 人人妻人人澡人人爽人人夜夜| 欧美黄色淫秽网站| 国产精品99久久99久久久不卡| 中文字幕人妻丝袜一区二区| 欧美日韩一级在线毛片| 9191精品国产免费久久| 日韩 欧美 亚洲 中文字幕| 18在线观看网站| 久久国产亚洲av麻豆专区| 国产日韩欧美亚洲二区| 国产精品久久久久久人妻精品电影 | 激情在线观看视频在线高清 | 国产一区二区三区综合在线观看| 一区二区av电影网| 下体分泌物呈黄色| 一本一本久久a久久精品综合妖精| 免费在线观看黄色视频的| 91九色精品人成在线观看| 免费在线观看影片大全网站| 国产男女超爽视频在线观看| 亚洲成国产人片在线观看| 咕卡用的链子| 人人妻,人人澡人人爽秒播| 51午夜福利影视在线观看| 久久精品国产99精品国产亚洲性色 | 久久久久久久精品吃奶| svipshipincom国产片| 飞空精品影院首页| 黄色视频,在线免费观看| 亚洲avbb在线观看| 欧美日韩黄片免| 国产精品av久久久久免费| 一级片'在线观看视频| 手机成人av网站| 99re6热这里在线精品视频| 一区二区av电影网| 国产精品一区二区精品视频观看| 欧美黑人欧美精品刺激| 亚洲成a人片在线一区二区| 91九色精品人成在线观看| 人妻 亚洲 视频| 欧美日韩成人在线一区二区| 成人三级做爰电影| 在线观看免费视频日本深夜| 久久精品国产亚洲av香蕉五月 | 亚洲精品自拍成人| 美女福利国产在线| 人妻 亚洲 视频| 国产精品美女特级片免费视频播放器 | 国产成人一区二区三区免费视频网站| 不卡av一区二区三区| 国产主播在线观看一区二区| 亚洲人成电影免费在线| 女性被躁到高潮视频| 国产精品二区激情视频| 成人国产av品久久久| 亚洲熟妇熟女久久| 免费久久久久久久精品成人欧美视频| 亚洲精品中文字幕在线视频| 国产免费福利视频在线观看| 99精品在免费线老司机午夜| 一个人免费看片子| 高清毛片免费观看视频网站 | 日本精品一区二区三区蜜桃| 18禁裸乳无遮挡动漫免费视频| 亚洲三区欧美一区| 国产一区二区三区在线臀色熟女 | 叶爱在线成人免费视频播放| 国产精品久久电影中文字幕 | 久久久久网色| 老汉色av国产亚洲站长工具| 国产成人欧美在线观看 | 精品少妇内射三级| 波多野结衣av一区二区av| 男女之事视频高清在线观看| 高清在线国产一区| 成人特级黄色片久久久久久久 | 极品教师在线免费播放| 丰满人妻熟妇乱又伦精品不卡| 亚洲成人免费电影在线观看| 啦啦啦 在线观看视频| 午夜老司机福利片| 精品国产国语对白av| 日韩精品免费视频一区二区三区| 桃红色精品国产亚洲av| 久久 成人 亚洲| 欧美 日韩 精品 国产| av网站免费在线观看视频| 国产成人精品久久二区二区免费| 在线永久观看黄色视频| 亚洲人成伊人成综合网2020| 亚洲成人免费电影在线观看| 国产亚洲欧美在线一区二区| 精品午夜福利视频在线观看一区 | 亚洲精品国产精品久久久不卡| 亚洲avbb在线观看| 国产深夜福利视频在线观看| 这个男人来自地球电影免费观看| 女人爽到高潮嗷嗷叫在线视频| 天堂动漫精品| 大片电影免费在线观看免费| 在线亚洲精品国产二区图片欧美| 高清黄色对白视频在线免费看| 99精品在免费线老司机午夜| 精品一品国产午夜福利视频| 精品国产一区二区久久| 变态另类成人亚洲欧美熟女 | 69精品国产乱码久久久| 曰老女人黄片| 色94色欧美一区二区| 成人18禁高潮啪啪吃奶动态图| 无遮挡黄片免费观看| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美精品综合一区二区三区| 亚洲人成77777在线视频| 欧美久久黑人一区二区| 精品国产一区二区三区四区第35| 视频区图区小说| 午夜两性在线视频| 日本a在线网址| 国产1区2区3区精品| 久久久精品免费免费高清| 精品一区二区三区四区五区乱码| 在线观看www视频免费| 男女下面插进去视频免费观看| 午夜日韩欧美国产| 欧美黄色片欧美黄色片| 欧美日韩亚洲国产一区二区在线观看 | 欧美一级毛片孕妇| 成人国产一区最新在线观看| 午夜福利在线观看吧| 午夜91福利影院| 黄色视频不卡| 欧美大码av| 一区福利在线观看| 一本—道久久a久久精品蜜桃钙片| 蜜桃国产av成人99| av网站免费在线观看视频| 免费在线观看视频国产中文字幕亚洲| 老汉色av国产亚洲站长工具| av一本久久久久| 日韩 欧美 亚洲 中文字幕| 日本vs欧美在线观看视频| 美女视频免费永久观看网站| 久久免费观看电影| 99国产精品一区二区三区| 90打野战视频偷拍视频| 精品少妇久久久久久888优播| 亚洲视频免费观看视频| 一级片'在线观看视频| 在线观看66精品国产| 在线观看免费视频日本深夜| 国产免费福利视频在线观看| 真人做人爱边吃奶动态| 一本一本久久a久久精品综合妖精| 国产成人啪精品午夜网站| 日韩成人在线观看一区二区三区| 别揉我奶头~嗯~啊~动态视频| 一区二区三区激情视频| 少妇裸体淫交视频免费看高清 | 一区二区三区国产精品乱码| 黄片小视频在线播放| 国产精品免费大片| 天天影视国产精品| 国产成人免费观看mmmm| 在线亚洲精品国产二区图片欧美| 欧美亚洲 丝袜 人妻 在线| 老熟妇乱子伦视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 99国产精品99久久久久| 啦啦啦 在线观看视频| 50天的宝宝边吃奶边哭怎么回事| 国产av一区二区精品久久| 大片电影免费在线观看免费| 丁香六月天网| 视频区图区小说| 老汉色∧v一级毛片| 国产精品成人在线| 日韩一区二区三区影片| 电影成人av| 亚洲一区中文字幕在线| 香蕉丝袜av| netflix在线观看网站| 成人免费观看视频高清| 黄色怎么调成土黄色| 亚洲成人国产一区在线观看| 午夜成年电影在线免费观看| 免费少妇av软件| 免费不卡黄色视频| 国产在线一区二区三区精| 99久久国产精品久久久| 欧美 日韩 精品 国产| 亚洲色图 男人天堂 中文字幕| 一区二区三区激情视频| 汤姆久久久久久久影院中文字幕| 国产精品99久久99久久久不卡| 99国产精品一区二区蜜桃av | 国产精品亚洲一级av第二区| 亚洲一区二区三区欧美精品| 曰老女人黄片| 亚洲国产精品一区二区三区在线| 欧美人与性动交α欧美软件| 俄罗斯特黄特色一大片| 丝袜美足系列| 亚洲av日韩精品久久久久久密| 一级毛片精品| 久久久久精品人妻al黑| 国产成人av激情在线播放| 人妻一区二区av| 我的亚洲天堂| 十八禁网站网址无遮挡| 看免费av毛片| 亚洲精品粉嫩美女一区| 久久99热这里只频精品6学生| 亚洲成av片中文字幕在线观看| 丰满迷人的少妇在线观看| 欧美老熟妇乱子伦牲交| 高清在线国产一区| 丝瓜视频免费看黄片| 欧美日韩视频精品一区| 777久久人妻少妇嫩草av网站| 大香蕉久久网| 久久人妻av系列| 天天影视国产精品| 国产无遮挡羞羞视频在线观看| 亚洲性夜色夜夜综合| 国产激情久久老熟女| 成人国产一区最新在线观看| 香蕉国产在线看| 亚洲三区欧美一区| 精品国产乱码久久久久久男人| 中文字幕高清在线视频| 国产免费av片在线观看野外av| 夜夜夜夜夜久久久久| 欧美 日韩 精品 国产| 青草久久国产| 欧美日韩精品网址| 亚洲av第一区精品v没综合| 天天躁日日躁夜夜躁夜夜| 丝袜美腿诱惑在线| 老司机深夜福利视频在线观看| 高清黄色对白视频在线免费看| 9热在线视频观看99| 国产精品 国内视频| 亚洲精品成人av观看孕妇| 1024视频免费在线观看| av免费在线观看网站| 国产成人精品无人区| a级毛片黄视频| 人成视频在线观看免费观看| 97人妻天天添夜夜摸| 久久天躁狠狠躁夜夜2o2o| 欧美成人午夜精品| 日韩大码丰满熟妇| 亚洲精品在线美女| 欧美日韩视频精品一区| 久久久久久久国产电影| 日本vs欧美在线观看视频| 亚洲欧美色中文字幕在线| 欧美另类亚洲清纯唯美| 午夜福利一区二区在线看| 搡老乐熟女国产| 91精品三级在线观看| 啦啦啦中文免费视频观看日本| 国产精品.久久久| 一区二区三区乱码不卡18| 在线看a的网站| 欧美精品亚洲一区二区| a在线观看视频网站| 国产一区二区三区视频了| av福利片在线| 国产在线一区二区三区精| 久久久久国产一级毛片高清牌| 午夜福利乱码中文字幕| 国产精品一区二区在线不卡| 人人妻人人爽人人添夜夜欢视频| 亚洲精品久久午夜乱码| 五月开心婷婷网| 一区二区三区乱码不卡18| 满18在线观看网站| 少妇 在线观看| 男人操女人黄网站| 国产在线一区二区三区精| 久久久久视频综合| 在线观看免费视频网站a站| 国产aⅴ精品一区二区三区波| 午夜福利,免费看| 中文字幕av电影在线播放| 国产有黄有色有爽视频| 欧美国产精品va在线观看不卡| 人成视频在线观看免费观看| 在线观看免费视频日本深夜| 久久性视频一级片| 大型黄色视频在线免费观看| e午夜精品久久久久久久| 美女午夜性视频免费| www.熟女人妻精品国产| 91精品三级在线观看| 国产成人精品久久二区二区免费| 亚洲七黄色美女视频| 在线观看66精品国产| 欧美乱妇无乱码| 亚洲伊人色综图| 成年女人毛片免费观看观看9 | 亚洲精品一卡2卡三卡4卡5卡| 黄片大片在线免费观看| 老熟妇仑乱视频hdxx| 国产精品亚洲av一区麻豆| 久久久精品免费免费高清| 精品国产乱子伦一区二区三区| 亚洲av电影在线进入| 国产精品熟女久久久久浪| 中国美女看黄片| 狠狠精品人妻久久久久久综合| 18禁裸乳无遮挡动漫免费视频| 少妇被粗大的猛进出69影院| 欧美久久黑人一区二区| 国产男女超爽视频在线观看| 老熟妇乱子伦视频在线观看| 欧美成狂野欧美在线观看| 精品亚洲乱码少妇综合久久| 一级片'在线观看视频| 99精品欧美一区二区三区四区| 热re99久久国产66热| 国产成人精品久久二区二区91| 国产日韩欧美视频二区| 视频区图区小说| 欧美在线黄色| 精品一区二区三区av网在线观看 | 色婷婷av一区二区三区视频| 国产又色又爽无遮挡免费看| 精品第一国产精品| 久久精品人人爽人人爽视色| 下体分泌物呈黄色| 一区在线观看完整版| 亚洲成人免费av在线播放| 国产高清videossex| 国产精品.久久久| 成年人午夜在线观看视频|