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

    綠色屋頂降雨徑流削減效果監(jiān)測與過程模擬

    2020-04-24 06:46:14申紅彬徐宗學(xué)張書函殷瑞雪
    農(nóng)業(yè)工程學(xué)報 2020年5期
    關(guān)鍵詞:徑流系數(shù)場次產(chǎn)流

    申紅彬,徐宗學(xué),張書函,殷瑞雪

    綠色屋頂降雨徑流削減效果監(jiān)測與過程模擬

    申紅彬1,2,徐宗學(xué)2,張書函3,殷瑞雪3

    (1. 華北水利水電大學(xué),河南省水利水運(yùn)工程技術(shù)研究中心,鄭州 450045;2. 北京師范大學(xué)水科學(xué)研究院,城市水循環(huán)與海綿城市技術(shù)北京市重點(diǎn)實(shí)驗室,北京 100875;3. 北京市水科學(xué)技術(shù)研究院,北京市非常規(guī)水資源開發(fā)利用與節(jié)水工程技術(shù)研究中心,北京 100048)

    選擇平面、斜面、綠色、滯蓄4種不同形式屋頂,采用“液位計+三角堰”的徑流測量方法,開展降雨徑流對比監(jiān)測。結(jié)果表明,相比平面與斜面屋頂平均徑流系數(shù)(分別為0.68和0.84),綠色與滯蓄屋頂平均徑流系數(shù)分別為0.41和0.43,均明顯有所減小,并且數(shù)值相差不大,反映出相近的徑流總量削減效應(yīng)。比較分析綠色與滯蓄屋頂降雨徑流過程,2類屋頂徑流過程接近,產(chǎn)流后水流運(yùn)動機(jī)制相似,均受上部積水水頭作用影響,并以孔流(孔隙、孔口)形式出流,可以采用相似的水庫概念模型進(jìn)行簡化模擬。根據(jù)土壤水飽和下滲理論,考慮到綠色屋頂基質(zhì)層厚度有限,推導(dǎo)建立了描述綠色屋頂飽和產(chǎn)流后降雨徑流過程的基本方程,并應(yīng)用于模擬綠色屋頂降雨徑流過程。結(jié)果表明,徑流過程計算值與實(shí)測值間決定系數(shù)值為0.85、Nash-Sutcliffe效率系數(shù)值為0.84,初步驗證了模型的合理性。研究成果有助于深入認(rèn)識綠色屋頂對降雨徑流的削減效應(yīng)與出流機(jī)制。

    降雨;徑流;模擬;綠色屋頂;監(jiān)測

    0 引 言

    綠色基礎(chǔ)設(shè)施(green infrastructure,GI)是實(shí)現(xiàn)社會與經(jīng)濟(jì)(城市、農(nóng)村)可持續(xù)發(fā)展的一種重要理念。綠色屋頂(green roof)作為綠色基礎(chǔ)設(shè)施重要組成部分之一,因能有效緩解熱島效應(yīng)、減少建筑能耗、削減暴雨徑流,從而引起人們的廣泛關(guān)注。綠色屋頂也稱屋頂綠化,是以傳統(tǒng)不透水屋頂為載體,以植物為配置主體的一種綠化形式,結(jié)構(gòu)由上至下依次為植被層、基質(zhì)層、過濾層、(蓄)排水層和防水層[1],坡度一般小于4%[2]。其中,基質(zhì)層對雨水的滲透吸持作用是決定雨水蓄存能力的主要因素[3-4],植被層通過蒸散發(fā)耗水以維持基質(zhì)層的雨水滲透吸持能力[5-6]。另外,植物莖葉也對雨水具有截留作用[7],并隨株高、覆蓋度和生物量的增大而逐漸增大[8-10]。

    已有大量研究證實(shí)了綠色屋頂對降雨徑流削減作用(包括徑流總量、峰值的削減、峰現(xiàn)時間的延遲等)的有效性[11-13],并有不少學(xué)者關(guān)注綠色屋頂結(jié)構(gòu)特征(包括植被種類、基質(zhì)厚度、屋頂坡度、運(yùn)行時間等)對降雨徑流削減效果差異性的影響,以用于確定不同的徑流模型參數(shù)[14-18]。根據(jù)質(zhì)量守恒原理,綠色屋頂徑流量為降水與4個組成(植物截留、基質(zhì)滲持、植物蒸散發(fā)、基質(zhì)蒸發(fā)量)之和的差值。綠色屋頂基質(zhì)層內(nèi)土壤滲流可以采用Richards方程進(jìn)行描述,相關(guān)模型有Hydrus-1D、SWMM等[19-21]。但是,因綠色屋頂基質(zhì)層厚度有限,故雨水蓄存能力有限,在基質(zhì)層內(nèi)土壤達(dá)到飽和含水量之前不會產(chǎn)生徑流。對于綠色屋頂,在基質(zhì)層內(nèi)土壤達(dá)到飽和含水量之后的徑流過程和規(guī)律,尚少有結(jié)合基質(zhì)屋厚度有限性條件進(jìn)行研究的簡化模擬方法。

    本文選擇平面屋頂、斜面屋頂、綠色屋頂、滯蓄屋頂4種不同形式屋頂為研究對象,采用“液位計+三角堰”的徑流測量方法[22],分別系統(tǒng)開展4種屋頂降雨徑流的現(xiàn)場監(jiān)測,通過比較4種不同形式屋頂?shù)慕涤陱搅魈匦?,結(jié)合綠色屋頂基質(zhì)屋厚度有限性的特點(diǎn),探討建立綠色屋頂降雨徑流過程的簡化模擬方法,以期能深入理解綠色屋頂對降雨徑流的削減效應(yīng)與出流機(jī)制。

    1 屋頂類型與監(jiān)測方案

    1.1 屋頂布置概況

    4種不同形式屋頂(平面屋頂、斜面屋頂、綠色屋頂、滯蓄屋頂)位于北京市水科學(xué)技術(shù)研究院屋頂平臺,地理位置與平面布置示意如圖1a所示。普通平面屋頂(坡度:0)面積約96 m2,屋頂鋪有防水油毛氈;斜面屋頂(坡度:5%)面積約85 m2,屋頂材料為鋼板;綠色屋頂(坡度:0;綠化比例:100%)面積約65 m2;滯蓄屋頂(坡度:0)面積約24 m2,屋頂鋪有防水油毛氈。每個屋頂對應(yīng)1個排水管(管徑110 mm)。其中,平面屋頂、斜面屋頂為原有屋頂,綠色屋頂與滯蓄屋頂系在原平面屋頂基礎(chǔ)上改建而成。綠色屋頂有關(guān)構(gòu)造如下:1)植被層草種選擇佛甲草,具有抗旱節(jié)水、隔熱降溫、易于管理等特點(diǎn);2)基質(zhì)層為采用草炭土、蛭石和砂土混合而成的填料,配比為4:2:1,基質(zhì)層厚度為6 cm,具有質(zhì)量輕、透水性好、持水性好、性能穩(wěn)定、養(yǎng)護(hù)方便等特點(diǎn);3)過濾層位于基質(zhì)層底部,材料為過濾布(厚度1~2 cm),防止介質(zhì)流失;4)排水層厚5 cm,用輕質(zhì)塑料制成,均勻布置碗狀結(jié)構(gòu)以承載徑流,并有排水出口。見圖1b。滯蓄屋頂主要在屋面排水管入口附近加裝攔污格柵與穿孔擋板,攔污格柵底部抬高對雨水具有一定的蓄存作用,穿孔擋板對雨水徑流具有一定的調(diào)控作用,穿孔擋板結(jié)構(gòu)形式參見圖1a。

    圖1 不同形式屋頂布置示意圖

    1.2 降雨徑流監(jiān)測方案

    降雨量監(jiān)測:屋頂安裝雨量計,可獲得1和5 min雨量值。從中國氣象數(shù)據(jù)共享服務(wù)網(wǎng)獲取北京氣象站(54511)逐日降雨數(shù)據(jù),以用于比較補(bǔ)充。

    徑流過程監(jiān)測:采用“液位計+三角堰”,可以對水位及流量進(jìn)行連續(xù)監(jiān)測。液位計選用WL-16壓力液位計(美國Global Water公司生產(chǎn)),該產(chǎn)品測量范圍包括0~0.9、0~4.6、0~9 m等多種選擇,恒溫條件下測量精度為±0.1%滿量程,存儲容量為81 759個記錄,用2塊9 V鋰電池供電。試驗過程中,分別在4類屋頂排水管出口下方各配置1套三角堰測流槽,在槽內(nèi)安裝液位計,見圖1a。三角堰頂角角度均為30°,平坡、斜坡屋頂對應(yīng)三角堰堰底面積為0.48 m2,綠色、滯蓄屋頂對應(yīng)三角堰堰底面積為0.4 m2。

    流量計算方法:三角堰流量計算在溢流前、后有所區(qū)別,溢流前計算時段內(nèi)屋頂徑流量為時段始末水位差值與堰底面積的乘積,溢流后過堰流量采用Kindsvater-Shen公式[23]:

    根據(jù)質(zhì)量守恒原理,三角堰溢流后水流連續(xù)方程為

    結(jié)合式(1)和式(2),在堰內(nèi)水位與過堰流量測算基礎(chǔ)上,通過逐步反向遞推,可以得到當(dāng)三角堰溢流后屋頂徑流過程。

    徑流削減效應(yīng)計算方法:徑流削減效應(yīng)可以分別采用徑流系數(shù)或降雨截留率衡量。徑流系數(shù)越小或降雨截留率越大,則徑流削減效應(yīng)越好。徑流系數(shù)定義為場次降雨在匯水面積上產(chǎn)生的徑流深和降雨量的比值:

    式中為徑流系數(shù);為降雨量,mm;為徑流深,mm,等于徑流量與匯水面積的比值。

    據(jù)質(zhì)量守恒原理,因降雨天氣下土壤植被蒸發(fā)蒸騰量相對較小可暫予忽略,則場次降雨截留深與截留率計算公式分別為

    =?(4)

    D=/=1?(5)

    式中為場次降雨截留深;D為場次降雨截留率。

    2 結(jié)果與分析

    2.1 不同屋頂徑流總量削減效應(yīng)

    2015年6—9月期間,對4種類型屋頂(平面屋頂、斜面屋頂、綠色屋頂、滯蓄屋頂)降雨徑流過程進(jìn)行連續(xù)監(jiān)測,共計監(jiān)測到11場有效降雨徑流數(shù)據(jù),統(tǒng)計各場次降雨徑流系數(shù)如表1所示。通過比較不同類型屋頂徑流系數(shù)的大小,可用以分析判別綠色屋頂?shù)膹搅骺偭肯鳒p效應(yīng)。有研究表明對于防水油毛氈,當(dāng)降雨量大于1.0 mm時就會產(chǎn)生徑流[10]。從表1可以看出,監(jiān)測各場次降雨量位于1.2~61.4 mm之間,平面、斜面、綠色、滯蓄屋頂徑流系數(shù)變化區(qū)間分別為:0.11~0.78、0.17~0.96、0.01~0.58、0.01~0.55,結(jié)合降雨量加權(quán)平均值分別為:0.68、0.84、0.41、0.43,相比斜面與平面屋頂,綠色與滯蓄屋頂徑流系數(shù)均有所減小,反映出明顯的徑流削減效應(yīng),且綠色與滯蓄屋頂徑流系數(shù)減小程度較為接近。分析徑流系數(shù)變化規(guī)律,4類屋頂徑流系數(shù)總體均呈現(xiàn)出隨降雨量增大而增大的趨勢,同時還受前期降雨量、場次降雨間隔時間等因素的影響。綠色屋頂因蒸散發(fā)時間較長,徑流系數(shù)受前期降雨量、場次降雨間隔時間因素影響尤為明顯。如對于2015年7月28日07:15-08:45降雨,因上場降雨2015年7月27日20:00-7月28日00:40雨量較大,累計降雨量達(dá)到54.8 mm,2場降雨間隔時間較短,間隔時長小于7 h,導(dǎo)致綠色屋頂徑流系數(shù)達(dá)到0.25,相比雨量較為接近的2015年6月30日06:05-10:00與2015年8月2日01:15-08:15降雨,增大效應(yīng)明顯;對于2015年6月26日20:30-21:35降雨,相比上場降雨2015年6月26日01:50-06:35雨量相差不大,但2場降雨間隔時間較短,間隔時長約為14 h左右,導(dǎo)致徑流系數(shù)增大明顯,達(dá)到0.50。

    表1 4類屋頂降雨徑流監(jiān)測時期監(jiān)測結(jié)果

    為進(jìn)一步比較綠色與滯蓄屋頂徑流削減效果,基于表1按降雨量從小到大依次排序,分別計算綠色、滯蓄屋頂?shù)膱龃谓涤杲亓羯钆c截留率,結(jié)果如圖2所示。

    圖2 綠色、滯蓄屋頂場次降雨截留深與截留率

    在圖2中,綠色屋頂場次降雨截留率變化范圍為42.0%~99.0%,平均值為58.9%,滯蓄屋頂場次降雨截留率變化范圍為45.0%~99.0%,平均值為57.2%,兩者平均截留效果相差不大。分序列分析場次降雨截留深與截留率,對于降雨序列7~11,降雨量相對較大,綠色、滯蓄屋頂各場次降雨截留深除序列8降雨(2015年6月26日20:30-21:35)因距上場降雨(2015年6月26日01:50-06:35)間隔時間較短而有所減小外,其余各場變化不大,平均截留深分別為25.7、24.5 mm,相差不大,各場降雨截留率均遠(yuǎn)小于100%,為部分截留。對于降雨序列4~6,降雨量相對不大,滯蓄屋頂場次降雨截留深近似等于降雨量、截留率達(dá)到98%,基本實(shí)現(xiàn)完全截留、剩余2%監(jiān)測徑流可能是屋頂攔污格柵與排水管入口所圍面積產(chǎn)生的徑流;綠色屋頂除序列6降雨(2015年8月1日21:50-22:55)外場次降雨截留深近似等于降雨量、截留率達(dá)到97%以上,基本實(shí)現(xiàn)完全截留,而序列6降雨(2015年8月1日21:50-22:55)截留深略小于降雨量,可能是由于前期降雨(2015年7月27日20:00-7月28日00:40,雨量54.8 mm)雨量較大且間隔時間不長、屋頂植被土壤通過蒸散發(fā)恢復(fù)調(diào)蓄能力過程較慢所致。對于降雨序列1~3,降雨量相對較小,滯蓄屋頂場次降雨截留深近似等于降雨量,除序列2降雨(2015年6月30日06:05-10:00)外截留率可以達(dá)到98%以上,基本實(shí)現(xiàn)完全截留,序列2降雨(2015年6月30日06:05-10:00)截留率略小可能因距上場降雨(2015年6月29日13:15-15:35)間隔時間較短所致;綠色屋頂場次降雨截留深近似等于降雨量,除序列1降雨(2015年7月28日07:15-08:45)外截留率可以達(dá)到93%以上,基本實(shí)現(xiàn)完全截留,序列1降雨(2015年7月28日07:15-08:45)截留率略小主要是因上場降雨(2015年7月27日20:00-7月28日00:40,雨量54.8 mm)雨量較大且間隔時間較短所致。

    2.2 徑流過程比較分析

    從2.1節(jié)可以看出,綠色與滯蓄屋頂多場次降雨平均截留深與截留率相差不大。進(jìn)一步深入比較2類屋頂?shù)慕涤陱搅鬟^程,宜選擇降雨量較大、只能部分截留的降雨序列進(jìn)行比較。以2015年7月27日20:00-7月28日00:40和2015年8月7日19:00-20:40降雨為例,2場降雨量均較大并且數(shù)值相近(分別為54.8、52.0 mm),而且2場降雨的上場降雨量均相對較?。ǚ謩e為6.7、1.8 mm)、間隔時間較長均為5 d左右,分別繪制2場降雨條件下綠色與滯蓄屋頂逐時段截留深、徑流深的變化過程,如圖3a~圖3d所示??梢钥闯觯G色與滯蓄屋頂降雨徑流均經(jīng)歷了蓄集雨水(無徑流產(chǎn)生)與產(chǎn)生徑流2個階段,2類屋頂截留深與徑流深變化過程均較為相近,表明2類屋頂具有相似的雨水徑流調(diào)控過程。其中,不同場次降雨累計截留深在增大至一定程度后基本保持穩(wěn)定,分別對應(yīng)為該場降雨的最大截留量。最大截留量出現(xiàn)時間相比產(chǎn)流時間略有延后,這是由于基質(zhì)層土壤濕潤鋒運(yùn)移至底部產(chǎn)生徑流后,土壤并未完全飽和。分析比較綠色與滯蓄屋頂雨水徑流調(diào)控機(jī)理,綠色屋頂蓄集雨水主要依靠植被屋植被截留、基質(zhì)層土壤恢復(fù)飽和等,并且土水勢作為基質(zhì)層內(nèi)土壤水運(yùn)動的驅(qū)動因子,包括重力勢、壓力勢與基模勢,并依土壤是否飽和而組成不同。根據(jù)土壤水飽和下滲理論,忽略最大截留量出現(xiàn)時間與產(chǎn)流時間的延后差異,簡化認(rèn)為基質(zhì)層內(nèi)土壤水運(yùn)動以濕潤鋒為界,上部土壤含水率達(dá)到飽和,下部土壤保持初始含水率,呈活塞狀向下運(yùn)動。則綠色屋頂在基質(zhì)層濕潤鋒運(yùn)移至底部、土壤恢復(fù)飽和前不會產(chǎn)生徑流,只有在基質(zhì)層內(nèi)濕潤鋒運(yùn)移至底部、土壤達(dá)到飽和后才會產(chǎn)生徑流。滯蓄屋頂雨水蓄集主要是攔污格柵底部抬高形成的水庫蓄水作用(圖3e和圖3f),當(dāng)蓄水深度超過擋水高度后產(chǎn)生徑流。綠色屋頂產(chǎn)生徑流后,基質(zhì)層內(nèi)土壤水運(yùn)動主要受到上部積水水頭的作用,經(jīng)土壤孔隙流至排水層排放。相比滯蓄屋頂產(chǎn)流后受上部積蓄水頭影響的水庫多孔出流過程,2類屋頂水流運(yùn)動動力機(jī)制十分相似,產(chǎn)流后徑流過程可以采用相似的水庫概念模型進(jìn)行簡化模擬。

    圖3 綠色與滯蓄屋頂降雨徑流過程比較

    3 徑流過程簡化模擬模型及驗證

    模型方程:根據(jù)土壤水飽和下滲理論,下滲土柱上的作用合力包括表面積水深、下滲土柱長度、濕潤鋒處毛管上升高度以及空氣剩余壓力。因綠色屋頂基質(zhì)層厚度的有限性,當(dāng)基質(zhì)層內(nèi)土壤全部恢復(fù)飽和產(chǎn)生徑流后,下滲土柱長度將恒為基質(zhì)屋厚度,成為定值,如忽略濕潤鋒處毛管上升高度及空氣剩余壓力項,則下滲速率可以根據(jù)Darcy定律表示為式(6)形式,并因質(zhì)量守恒應(yīng)與基質(zhì)層底部出流速率相等:

    式中f(?0)為基質(zhì)層產(chǎn)流后土壤下滲速率,mm/min;(?0)為基質(zhì)層產(chǎn)流后底部出流速率,轉(zhuǎn)換為徑流深單位,mm/min;h(?0)為下滲土柱表面積水深,mm;為基質(zhì)層厚度,mm;(?0)為作用水頭,mm;K為飽和水力傳導(dǎo)度,mm/min;K為系數(shù),min-1;為計算時間(以降雨開始為起始時間),min;0為產(chǎn)流時間,min。

    對于作用水頭(?0),根據(jù)質(zhì)量守恒原理存在如下平衡方程:

    式中(?0)為基質(zhì)層產(chǎn)流后降雨強(qiáng)度,mm/min。

    聯(lián)合求解式(6)和式(7),得綠色屋頂產(chǎn)流后徑流過程的連續(xù)卷積形式:

    其中

    表示為離散形式:

    其中

    式中(0,)為瞬時單位線;(Δ,)為時段單位線;()為(0,)積分函數(shù)曲線;為計算時段數(shù);0為產(chǎn)流開始時段序號;為計算時段序號;為0~之間的時間變量,min。

    模型評價:選決定系數(shù)(2)和Nash-Sutcliffe效率系數(shù)(Nash-Sutcliffe efficiency,NSE)作為評價指標(biāo)[24-26]。利用Excel軟件對流量計算值與實(shí)測值進(jìn)行線性回歸,可以直接得到回歸方程的2值,2值越趨近于1,說明流量模型擬合效果越好。NSE主要用于水文模型效率評價,是判定殘差與實(shí)測值數(shù)據(jù)方差相對量的標(biāo)準(zhǔn)化統(tǒng)計值。NSE值一般在-∞~1之間,當(dāng)NSE=1時說明流量計算值與實(shí)測值完全吻合,當(dāng)0

    模型應(yīng)用:應(yīng)用模型方程離散形式式(10),分別對綠色屋頂2015年6月26日01:50—06:35、2015年6月26日20:30-21:35、2015年7月27日20:00—7月28日00:40及2015年8月7日19:00—20:40降雨徑流過程進(jìn)行模擬,該4場降雨雨量相對較大,均為部分截留,計算結(jié)果(轉(zhuǎn)換為時段徑流深)如圖4所示。其中,系數(shù)K經(jīng)過率定,取值為0.2,對應(yīng)飽和水力傳導(dǎo)度K值為12 mm/min,數(shù)值較大可能是因基質(zhì)層土壤內(nèi)混合含有砂土有關(guān);對于不同場次降雨綠色屋頂產(chǎn)流時間0,因基質(zhì)層土壤初始含水率不同而有所區(qū)別,可以根據(jù)不同場次降雨實(shí)測最大截留量與累計降雨量比較確定。比較4場降雨徑流計算值與實(shí)測值,結(jié)果見表2。結(jié)合圖4和表2可以看出,綠色屋頂降雨徑流過程計算值與實(shí)測值變化趨勢基本符合,綜合所有降雨徑流的2值及NSE值分別為0.85、0.84,初步驗證了所建模型方法的合理性。

    圖4 綠色屋頂降雨徑流過程模擬結(jié)果

    表2 綠色屋頂降雨徑流模擬效果評價

    4 結(jié) 論

    1)平面、斜面、綠色、滯蓄屋頂降雨徑流監(jiān)測結(jié)果表明,降雨量位于1.2~61.4 mm之間,4類屋頂徑流系數(shù)變化區(qū)間分別為0.11~0.78、0.17~0.96、0.01~0.58、0.01~0.55,加權(quán)平均值分別為0.68、0.84、0.41、0.43。相比平面與斜面屋頂,綠色與滯蓄屋頂徑流系數(shù)均明顯減小,并且數(shù)值相差不大,反映出相近的徑流總量削減效果。

    2)比較綠色屋頂與滯蓄屋頂?shù)慕涤陱搅鬟^程,2類屋頂徑流過程較為接近,均經(jīng)歷了蓄集雨水與產(chǎn)生徑流2個階段,產(chǎn)流后的水流運(yùn)動機(jī)制十分相似,均主要受到上部積水水頭作用的影響,并以孔流(孔隙、孔口)形式出流,可以采用相似的水庫概念模型進(jìn)行簡化模擬。

    3)根據(jù)土壤水飽和下滲理論,考慮綠色屋頂基質(zhì)層厚度的有限性,推導(dǎo)建立了描述綠色屋頂基質(zhì)層土壤飽和產(chǎn)流后降雨徑流過程的模型方程,并對監(jiān)測綠色屋頂降雨徑流過程進(jìn)行了模擬應(yīng)用。結(jié)果表明,徑流過程計算值與實(shí)測值變化趨勢基本符合,2值為0.85、Nash-Sutcliffe效率系數(shù)為0.84,初步驗證了模型的合理性。

    本文根據(jù)綠色、滯蓄屋頂徑流過程的相似性,結(jié)合理論分析,建立了用于描述綠色屋頂徑流過程的簡化水庫模型,尚未涉及綠色屋頂不同綠化形式對降雨徑流的影響。綠色屋頂不同綠化形式對降雨徑流的影響主要表現(xiàn)在不同植被雨水截流能力的差異以及土壤基質(zhì)含水率恢復(fù)速率的差異等,這些尚有待于后期進(jìn)一步的深入研究。除此之外,除氣溫、濕度等影響因子外,綠色屋頂不同綠化形式及土壤基質(zhì)對本文簡化模型參數(shù)(飽和水力傳導(dǎo)度和產(chǎn)流時間)的影響仍需進(jìn)一步深入研究。

    [1] Berndtsson J C. Green roof performance towards management of runoff water quantity and quality: A review[J]. Ecological Engineering, 2010, 36(4): 351-360.

    [2] 唐穎. SUSTAIN支持下的城市降雨徑流最佳管理BMP規(guī)劃研究[D]. 北京:清華大學(xué),2010.

    Tang Ying. Study of Urban Stormwater Runoff BMPs Planning with Support of SUSTAIN[D]. Beijing: Qsinghua University, 2010. (in Chinese with English abstract)

    [3] Monterusso M A, Rowe D B, Rugh C L. Establishment and persistence of Sedum spp. and native taxa for green roof applications[J]. Hortscience, 2005, 40(2): 391-396.

    [4] Dunnett N, Kingsbury N. Planting Green Roofs and Living Walls[M]. Portland: Timber Press, 2008.

    [5] Nagase A, Dunnett N. Amount of water runoff from different vegetation types on extensive green roofs: Effects of plant species, diversity and plant structure[J]. Landscape and Urban Planning, 2012, 104(3/4): 356-363.

    [6] 曹傳生,劉慧民,王南. 屋頂花園雨水利用系統(tǒng)設(shè)計與實(shí)踐[J]. 農(nóng)業(yè)工程學(xué)報,2013,29(9):76-85.

    Cao Chuansheng, Liu Huimin, Wang Nan. Design and practice of rainwater utilization system for green garden[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(9): 76-85. (in Chinese with English abstract)

    [7] Zhang S H, Guo Y P. Analytical probabilistic model for evaluating the hydrologic performance of green roofs[J]. Journal of Hydrologic Engineering, 2013, 18(1): 19-28.

    [8] Wolf D, Lundholm J T. Water uptake in green roof microcosms: effects of plant species and water availability[J]. Ecological Engineering, 2008, 33(2): 179-186.

    [9] Teemusk A, Mander U. Rainwater runoff quantity and quality performance from a green roof: The effects of short-term events[J]. Ecological Engineering, 2007, 30(3): 271-277.

    [10] 葛德,張守紅. 不同降雨條件下植被對綠色屋頂徑流調(diào)控效益影響[J]. 環(huán)境科學(xué),2018,39(11):5015-5023.

    Ge De, Zhang Shouhong. Impacts of vegetation on hydrological performances of green roofs under different rainfall conditions[J]. Environmental Science, 2018, 39(11): 5015-5023. (in Chinese with English abstract)

    [11] Mentens J, Raes D, Hermy M. Green roofs as a tool for solving the rainwater runoff problem in the urbanized 21stcentury[J]. Landscape Urban Plan, 2006, 77(3): 217-226.

    [12] Bliss D J, Neufeld R D, Ries R J. Storm water runoff mitigation using a green roof[J]. Environmental Engineering Science, 2009, 26(2): 407-418.

    [13] Bemdtsson J C. Green roof performance towards management of runoff water quantity and quality: A review[J]. Ecological Engineering, 2010, 36(4): 351-360.

    [14] Berndtsson J C, Bengtsson L, Jinno K. Runoff water quality from intensive and extensive vegetated roofs[J]. Ecological Engineering, 2009, 35(3): 369-380.

    [15] Beck D A, Johnson G R, Spolek G A. Amending greenroof soil with biochar to affect runoff water quantity and quality[J]. Environmental Pollution, 2011, 159(8/9): 2111-2118.

    [16] 盧浩,康成,譚松明,等. 屋頂綠化基質(zhì)對雨水徑流水質(zhì)的影響[J]. 中國給水排水,2016,32(7):135-138.

    Lu Hao, Kang Cheng, Tan Songming, et al. Effect of green roof substance on water quality of stromwater runoff[J]. China Water & Wastewater, 2016, 32(7): 135-138. (in Chinese with English abstract)

    [17] 秦華鵬,唐女,唐巧玲. 蓄水層對綠色屋頂徑流削減能力的影響分析[J]. 中國給水排水,2016,32(13):132-135.

    Qin Huapeng, Tang Nü, Tang Qiaoling. Impact of water storage layer on runoff reduction of green roofs[J]. China Water & Wastewater, 2016, 32(13): 132-135. (in Chinese with English abstract)

    [18] Getter K L, Rowe D B, Andresen J A. Quantifying the effect of slope on extensive green roof stormwater retention[J]. Ecological Engineering, 2007, 31(4): 225-231.

    [19] Hilten R N, Lawrence T M, Tollner E W. Modeling stormwater runoff from green roofs with Hydrus-1D[J]. Journal of Hydrology, 2008, 358(3/4): 288-293.

    [20] 唐莉華,倪廣恒,劉茂峰,等. 綠化屋頂?shù)漠a(chǎn)流規(guī)律及雨水滯蓄效果模擬研究[J]. 水文,2011,31(4):18-22.

    Tang Lihua, Ni Guangheng, Liu Maofeng, et al. Study on runoff and rainwater retention capacity of green roof by experiment and model simulation[J]. Journal of China Hydrology, 2011, 31(4): 18-22. (in Chinese with English abstract)

    [21] Palla A, Gnecco I, Lanza L G. Compared performance of a conceptual and a mechanistic hydrologic models of a green roof[J]. Hydrological Processes, 2012, 26(1): 73-84.

    [22] 廖愛民,劉九夫,張建云,等. 實(shí)驗小流域尺度下高精度水位計的比測分析[J]. 水科學(xué)進(jìn)展,2019,30(3):337-347.

    Liao Aimin, Liu Jiufu, Zhang Jianyun, et al. Inter comparison of high-accuracy water level gauges in the scale of small experiment catchment[J]. Advance in Water Science, 2019, 30(3): 337-347. (in Chinese with English abstract)

    [23] 吳長文,王禮先. 小集水區(qū)試驗中三角堰量水問題的研究[J]. 北京林業(yè)大學(xué)學(xué)報,1993,15(3):66-74.

    Wu Changwen, Wang Lixian. Study on the discharge characteristics of triangular-notch thin-plate weirs[J]. Journal of Beijing Forestry University, 1993, 15(3): 66-74. (in Chinese with English abstract)

    [24] 胡勝,曹明明,邱海軍,等. CFSR氣象數(shù)據(jù)在流域水文模擬中的適用性評價-以灞河流域為例[J]. 地理學(xué)報,2016,71(9):1571-1586.

    Hu Sheng, Cao Mingming, Qiu Haijun, et al. Applicability evaluation of CFSR climate data for hydrologic simulation: A case study in the Bahe River Basin[J]. ACTA Geographica Sinica, 2016, 71(9): 1571-1586. (in Chinese with English abstract)

    [25] Nash J E, Sutcliffe J V. River flow forecasting through conceptual models (Part 1): A discussion of principles[J]. Journal of Hydrology, 1970, 10(3): 282-290.

    [26] Gupta H V, Sorooshian S, Yapo P O. Status of automatic calibration for hydrologic models: Comparison with multilevel expert calibration[J]. Journal of Hydrologic Engineering, 1999, 4(2): 135-143.

    Monitoring and simulation of rainfall-runoff reduction effects and process on green roofs

    Shen Hongbin1,2, Xu Zongxue2, Zhang Shuhan3, Yin Ruixue3

    (1.,,450045,; 2.,,,100875,; 3.,,100048,)

    Green roof is an important facility to mitigate urban flood resulted by increasing urbanization and impermeable surface. A thick substrate material layer with plants, grass, and/or trees, and extra structural support is one of main components of green roof. Due to the complexity of rainwater storage mechanism in the substrate material layer, the simulations results on rainfall-runoff on green roof are mostly unsatisfactory. In this study, the runoff reduction effect of green roof was simulated by considering the effect of substrate material layer on runoff. The green roof was compared with slope roof, plane roof and retention roof in order to understand the effect of substrate material layer on runoff. The experiment was carried out in Beijing Water Science and Technology Institute, China. The precipitation was measured by using a self-recording tipping bucket rain gauge. Runoff was monitored by using an equipment including water level gauges and triangular weir. The runoff coefficient and rainwater retention rate were calculated. A total of 11 rainfall-runoff events were monitored in 2015. The results showed that the change of the runoff coefficient was highly complex and it was influenced by many factors, e.g., precipitation, average/maximum rainfall intensity, and interval time between rainfall events. For all the monitored rainfall events, the average runoff coefficient values for plane roof, slope roof, green roof and retention roof were 0.68, 0.84, 0.41 and 0.43, respectively. The average rainwater retention rate values of green roof and retention roof were 58.9% and 57.2%, respectively. It indicated that the runoff reduction effects of green roof and retention roof are similar but both were better than the other two types of roofs. The runoff process of green roof was also similar with that of retention roof under the same rainfall conditions. Both the outflow processes from green roof and retention roof could be regarded as an unsteady orifice flow under the effect of accumulated water head on the roof-top. Thus, a unit hydrograph model to simulate runoff from green roof was proposed based on the saturated infiltration theory that included the thickness limitation of the substrate material layer and the influence of water accumulated on the surface. By using the proposed unit hydrograph model, the rainfall-runoff processes on green roof were simulated. The simulation results showed that the trends of the simulated and measured runoff were consistent well. The coefficient of determination and Nash-Sutcliffe efficiency were 0.85 and 0.84, respectively. This study will be helpful for further understanding the runoff reduction effect and the outflow mechanism of green roofs.

    precipitation; runoff; simulation; green roof; monitoring

    2019-08-25

    2019-12-10

    城市水循環(huán)與海綿城市技術(shù)北京市重點(diǎn)實(shí)驗室開放基金(HYD2019OF02);國家科技重大科技專項課題(2013ZX07304-001)

    申紅彬,博士,講師,主要從事水文學(xué)、河流動力學(xué)方面研究工作。Email:hongbinshen 2012@163.com

    10.11975/j.issn.1002-6819.2020.05.020

    TV121+.1

    A

    1002-6819(2020)-05-0175-07

    申紅彬,徐宗學(xué),張書函,殷瑞雪. 綠色屋頂降雨徑流削減效果監(jiān)測與過程模擬[J]. 農(nóng)業(yè)工程學(xué)報,2020,36(5):175-181. doi:10.11975/j.issn.1002-6819.2020.05.020 http://www.tcsae.org

    Shen Hongbin, Xu Zongxue, Zhang Shuhan, Yin Ruixue. Monitoring and simulation of rainfall-runoff reduction effects and process on green roofs[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(5): 175-181. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2020.05.020 http://www.tcsae.org

    王笑亞 9787512343733 254 王笑亞 DXBL202002 122

    猜你喜歡
    徑流系數(shù)場次產(chǎn)流
    產(chǎn)流及其研究進(jìn)展
    降雨特征對半透水道路徑流系數(shù)的影響
    長江上游高洪水期泥沙輸移特性
    不同坡面單元人工降雨產(chǎn)流試驗與分析
    北京山區(qū)側(cè)柏林地坡面初始產(chǎn)流時間影響因素
    基于運(yùn)行場次用時誤差的載人設(shè)備故障預(yù)警可視化研究
    排考場次分配方法及其SQL實(shí)現(xiàn)
    地表粗糙度對黃土坡面產(chǎn)流機(jī)制的影響
    無資料山丘小流域徑流系數(shù)的計算研究
    地鐵觀影指南
    電影故事(2015年33期)2015-09-06 01:05:30
    婷婷精品国产亚洲av在线| 日韩亚洲欧美综合| 国产不卡一卡二| 免费观看在线日韩| 亚洲无线观看免费| 亚洲在线自拍视频| 性欧美人与动物交配| 日本一二三区视频观看| 国产精品一及| 国产真实伦视频高清在线观看 | 我的老师免费观看完整版| 日日摸夜夜添夜夜添小说| 99热这里只有是精品在线观看| 亚洲av熟女| 变态另类丝袜制服| 国产精品免费一区二区三区在线| 久久久午夜欧美精品| 在线天堂最新版资源| 18+在线观看网站| 可以在线观看毛片的网站| 麻豆国产av国片精品| а√天堂www在线а√下载| 国产成年人精品一区二区| 91久久精品电影网| 不卡视频在线观看欧美| 久久精品国产鲁丝片午夜精品 | 亚洲精品一卡2卡三卡4卡5卡| 国产精品不卡视频一区二区| 啦啦啦观看免费观看视频高清| 色av中文字幕| xxxwww97欧美| 亚洲欧美清纯卡通| 日韩av在线大香蕉| 亚洲七黄色美女视频| a级毛片免费高清观看在线播放| 99久久中文字幕三级久久日本| 欧美日韩精品成人综合77777| 午夜福利18| 婷婷精品国产亚洲av在线| 亚洲av中文av极速乱 | 观看美女的网站| 久久久久久国产a免费观看| 天堂动漫精品| 久久精品国产自在天天线| 亚洲真实伦在线观看| 身体一侧抽搐| av天堂中文字幕网| 国产精品不卡视频一区二区| 可以在线观看毛片的网站| 免费观看精品视频网站| 国产av一区在线观看免费| 亚洲av不卡在线观看| 啦啦啦啦在线视频资源| 久99久视频精品免费| 久久天躁狠狠躁夜夜2o2o| 毛片一级片免费看久久久久 | 美女 人体艺术 gogo| 熟女人妻精品中文字幕| 可以在线观看的亚洲视频| av女优亚洲男人天堂| 美女cb高潮喷水在线观看| 亚洲精品影视一区二区三区av| 岛国在线免费视频观看| 制服丝袜大香蕉在线| 欧美日韩中文字幕国产精品一区二区三区| 亚洲精品色激情综合| 久久久成人免费电影| 久久亚洲精品不卡| 波多野结衣高清无吗| 三级国产精品欧美在线观看| 国产精品一区二区性色av| 国产精品98久久久久久宅男小说| 国产午夜精品论理片| 成人欧美大片| 内地一区二区视频在线| 校园春色视频在线观看| 中出人妻视频一区二区| 久久99热6这里只有精品| 乱系列少妇在线播放| 国产真实乱freesex| 亚洲四区av| 欧美成人一区二区免费高清观看| 欧美日韩乱码在线| 欧美日韩乱码在线| 99国产精品一区二区蜜桃av| 舔av片在线| 欧美性猛交╳xxx乱大交人| 听说在线观看完整版免费高清| 欧美性猛交╳xxx乱大交人| 免费大片18禁| 亚洲国产精品成人综合色| 亚洲内射少妇av| 舔av片在线| 国产伦人伦偷精品视频| 欧美激情在线99| 日本免费一区二区三区高清不卡| 日韩 亚洲 欧美在线| 亚洲精品久久国产高清桃花| 又紧又爽又黄一区二区| 国产伦精品一区二区三区四那| 国产高清三级在线| 99热精品在线国产| 看黄色毛片网站| 亚洲成av人片在线播放无| 日韩大尺度精品在线看网址| 成人综合一区亚洲| 欧美+日韩+精品| 赤兔流量卡办理| 可以在线观看毛片的网站| 国产激情偷乱视频一区二区| 在线国产一区二区在线| 亚洲午夜理论影院| 欧美+日韩+精品| 老司机深夜福利视频在线观看| 欧美+日韩+精品| 少妇高潮的动态图| 国产白丝娇喘喷水9色精品| 永久网站在线| 亚洲乱码一区二区免费版| 亚洲真实伦在线观看| 亚洲av中文av极速乱 | 久久亚洲真实| 国产成人影院久久av| 给我免费播放毛片高清在线观看| 简卡轻食公司| 中文字幕av成人在线电影| 日韩一区二区视频免费看| 国产精品野战在线观看| 人妻丰满熟妇av一区二区三区| 99热这里只有是精品50| 国产激情偷乱视频一区二区| 波多野结衣巨乳人妻| 啦啦啦韩国在线观看视频| 国产 一区 欧美 日韩| 大又大粗又爽又黄少妇毛片口| 亚洲午夜理论影院| 亚洲国产精品合色在线| 伦理电影大哥的女人| 国产精品久久久久久久电影| 麻豆久久精品国产亚洲av| 大型黄色视频在线免费观看| 午夜福利18| 草草在线视频免费看| 国产精品,欧美在线| 久久久久久久久久黄片| 国产高清视频在线观看网站| 久久久久久久午夜电影| 久久久精品欧美日韩精品| 日本免费a在线| 一级毛片久久久久久久久女| 国产亚洲91精品色在线| 成人欧美大片| 又爽又黄a免费视频| 99riav亚洲国产免费| 一本一本综合久久| 深夜精品福利| 日本 av在线| 欧美另类亚洲清纯唯美| 午夜老司机福利剧场| 麻豆成人午夜福利视频| 美女高潮的动态| 精品久久久久久成人av| 欧美极品一区二区三区四区| 日本三级黄在线观看| 久久久精品大字幕| 波多野结衣巨乳人妻| 亚洲午夜理论影院| 亚洲经典国产精华液单| 动漫黄色视频在线观看| 午夜福利在线观看免费完整高清在 | 大又大粗又爽又黄少妇毛片口| 欧美日本视频| 精品久久久久久久久久免费视频| 欧美高清成人免费视频www| 亚洲色图av天堂| 日本黄大片高清| 欧美+亚洲+日韩+国产| 久久精品国产自在天天线| 久久精品夜夜夜夜夜久久蜜豆| 一进一出好大好爽视频| 一区福利在线观看| 久久久久久久久久黄片| 嫩草影院精品99| 啪啪无遮挡十八禁网站| 99久久精品国产国产毛片| 午夜激情欧美在线| 国产精品伦人一区二区| 国产大屁股一区二区在线视频| 免费在线观看日本一区| 18+在线观看网站| 国产精品98久久久久久宅男小说| 国产精品一区二区性色av| 免费看a级黄色片| 久久精品夜夜夜夜夜久久蜜豆| 精品久久国产蜜桃| 免费看美女性在线毛片视频| 午夜免费激情av| 麻豆av噜噜一区二区三区| 日本一本二区三区精品| 国产亚洲精品综合一区在线观看| av天堂中文字幕网| 国产精品久久久久久亚洲av鲁大| 丰满人妻一区二区三区视频av| 亚洲一级一片aⅴ在线观看| 国产一区二区三区在线臀色熟女| 夜夜爽天天搞| 热99re8久久精品国产| 天堂影院成人在线观看| 久久久久久久亚洲中文字幕| 午夜免费男女啪啪视频观看 | 国产在线精品亚洲第一网站| 精品久久久久久久久久免费视频| 欧美绝顶高潮抽搐喷水| 国产精品久久久久久久电影| 免费看a级黄色片| 长腿黑丝高跟| 精品免费久久久久久久清纯| 一进一出抽搐动态| av黄色大香蕉| 欧美+日韩+精品| 1000部很黄的大片| 五月伊人婷婷丁香| 国产精品综合久久久久久久免费| 欧美区成人在线视频| 老熟妇仑乱视频hdxx| 3wmmmm亚洲av在线观看| 国产高清视频在线观看网站| 亚洲av一区综合| 欧美在线一区亚洲| 亚洲国产色片| 中文字幕免费在线视频6| 午夜免费男女啪啪视频观看 | 伊人久久精品亚洲午夜| 欧美最新免费一区二区三区| 亚洲成人免费电影在线观看| 国产真实伦视频高清在线观看 | 欧美极品一区二区三区四区| 免费一级毛片在线播放高清视频| 成年版毛片免费区| 欧美又色又爽又黄视频| 99久久精品热视频| 深爱激情五月婷婷| 男人的好看免费观看在线视频| 亚洲av电影不卡..在线观看| 国内精品久久久久精免费| 亚洲av五月六月丁香网| 午夜福利高清视频| 免费av毛片视频| 嫩草影院精品99| 国产淫片久久久久久久久| 中文字幕免费在线视频6| 国产成人一区二区在线| 悠悠久久av| 日日撸夜夜添| 一夜夜www| 国产成人av教育| 一本精品99久久精品77| 亚洲,欧美,日韩| 中文字幕av成人在线电影| 色综合亚洲欧美另类图片| 国产一区二区亚洲精品在线观看| 在线观看午夜福利视频| 国产真实乱freesex| 精品人妻一区二区三区麻豆 | 99精品久久久久人妻精品| 级片在线观看| 亚洲色图av天堂| 日韩欧美精品v在线| 亚洲精品国产成人久久av| 精品一区二区三区视频在线| 中文在线观看免费www的网站| 亚洲欧美日韩东京热| 三级国产精品欧美在线观看| 国产主播在线观看一区二区| 美女xxoo啪啪120秒动态图| 国内精品宾馆在线| 亚洲第一电影网av| 欧美日韩国产亚洲二区| 国产高潮美女av| 99热网站在线观看| 黄色视频,在线免费观看| 国模一区二区三区四区视频| 美女高潮的动态| 精品久久久久久成人av| 九色国产91popny在线| 搞女人的毛片| 亚洲电影在线观看av| 精品一区二区三区人妻视频| 五月伊人婷婷丁香| 国产精品一及| 长腿黑丝高跟| 久久午夜福利片| 国产精品久久电影中文字幕| 亚洲精华国产精华液的使用体验 | 琪琪午夜伦伦电影理论片6080| 99久久精品国产国产毛片| 国产精品永久免费网站| 欧美精品啪啪一区二区三区| 深夜a级毛片| 深夜精品福利| 99久久精品国产国产毛片| 午夜久久久久精精品| 白带黄色成豆腐渣| 国产不卡一卡二| 欧美激情在线99| 色av中文字幕| 欧美日韩国产亚洲二区| 日本五十路高清| 舔av片在线| 国产精品嫩草影院av在线观看 | 少妇高潮的动态图| 亚洲成人中文字幕在线播放| 一区二区三区免费毛片| 成熟少妇高潮喷水视频| 中国美白少妇内射xxxbb| 日日啪夜夜撸| 俄罗斯特黄特色一大片| 国内精品久久久久久久电影| 午夜视频国产福利| 欧美又色又爽又黄视频| 淫秽高清视频在线观看| 麻豆成人午夜福利视频| 国产成人a区在线观看| 亚洲性久久影院| 91午夜精品亚洲一区二区三区 | 成熟少妇高潮喷水视频| 成年版毛片免费区| 搡女人真爽免费视频火全软件 | 91久久精品电影网| 欧美精品国产亚洲| 麻豆国产av国片精品| 啦啦啦韩国在线观看视频| 99久久久亚洲精品蜜臀av| 国产亚洲91精品色在线| 日韩 亚洲 欧美在线| 麻豆成人av在线观看| 99久久精品热视频| 亚洲精品久久国产高清桃花| 直男gayav资源| 亚洲国产精品合色在线| 少妇被粗大猛烈的视频| 嫩草影院入口| 黄色丝袜av网址大全| 日本黄色视频三级网站网址| 丰满的人妻完整版| 舔av片在线| 夜夜夜夜夜久久久久| 人妻制服诱惑在线中文字幕| 日本 av在线| 国产精品久久视频播放| 99久久成人亚洲精品观看| 亚洲第一区二区三区不卡| 精品久久久久久久久av| 久久久色成人| 99热这里只有是精品50| 1000部很黄的大片| 色综合站精品国产| 成人av在线播放网站| 亚洲精华国产精华液的使用体验 | 乱系列少妇在线播放| 午夜视频国产福利| 亚洲美女黄片视频| 99国产精品一区二区蜜桃av| 精品午夜福利视频在线观看一区| 欧美日韩瑟瑟在线播放| 精品人妻1区二区| 欧美成人性av电影在线观看| 久久久国产成人精品二区| .国产精品久久| 亚洲av免费在线观看| 18禁黄网站禁片午夜丰满| 欧美人与善性xxx| 草草在线视频免费看| 毛片女人毛片| 窝窝影院91人妻| 国内少妇人妻偷人精品xxx网站| 国产精品自产拍在线观看55亚洲| 亚洲人成伊人成综合网2020| 天天一区二区日本电影三级| 91午夜精品亚洲一区二区三区 | 亚洲欧美清纯卡通| 非洲黑人性xxxx精品又粗又长| 久久久国产成人精品二区| 国产aⅴ精品一区二区三区波| 国产精品人妻久久久久久| 日韩精品中文字幕看吧| 亚洲成人久久性| avwww免费| 听说在线观看完整版免费高清| 成人无遮挡网站| 亚洲欧美精品综合久久99| 日韩一区二区视频免费看| 中文亚洲av片在线观看爽| 美女免费视频网站| 久久精品综合一区二区三区| 日韩欧美精品免费久久| 人妻久久中文字幕网| 欧美不卡视频在线免费观看| 欧美精品国产亚洲| 国产视频内射| 变态另类成人亚洲欧美熟女| 亚洲18禁久久av| 午夜爱爱视频在线播放| 国产亚洲精品av在线| 国产精品久久久久久精品电影| 毛片女人毛片| 欧美日本视频| 国内揄拍国产精品人妻在线| 欧美色视频一区免费| 国产探花极品一区二区| 精品一区二区免费观看| 能在线免费观看的黄片| 成年版毛片免费区| 亚洲av成人av| 麻豆av噜噜一区二区三区| 小说图片视频综合网站| 亚洲av日韩精品久久久久久密| 91麻豆av在线| 成年女人看的毛片在线观看| 极品教师在线免费播放| 亚洲人成网站在线播| 成人av一区二区三区在线看| 一个人免费在线观看电影| 亚洲av成人精品一区久久| 人人妻人人澡欧美一区二区| 在线观看午夜福利视频| 午夜日韩欧美国产| 婷婷六月久久综合丁香| 亚洲欧美日韩高清专用| 日韩欧美在线乱码| 九九久久精品国产亚洲av麻豆| 国产成人av教育| 小说图片视频综合网站| or卡值多少钱| 日本一本二区三区精品| 午夜免费成人在线视频| 国产精品久久视频播放| 亚洲精华国产精华液的使用体验 | 成人二区视频| 成年女人看的毛片在线观看| 日本免费a在线| 国产色婷婷99| 亚洲av.av天堂| 国产精品美女特级片免费视频播放器| 久久精品国产亚洲av涩爱 | av黄色大香蕉| 一进一出抽搐gif免费好疼| 国产精品女同一区二区软件 | 99久久精品热视频| 国产视频内射| 国产精品久久久久久久电影| 深爱激情五月婷婷| 99热这里只有是精品50| 97人妻精品一区二区三区麻豆| 久9热在线精品视频| 久久久久免费精品人妻一区二区| 男女下面进入的视频免费午夜| 成人性生交大片免费视频hd| 麻豆成人av在线观看| 极品教师在线免费播放| 精品久久久噜噜| 国产主播在线观看一区二区| 国产成人福利小说| 日韩强制内射视频| 日本在线视频免费播放| 久久热精品热| 免费人成视频x8x8入口观看| av专区在线播放| 色综合亚洲欧美另类图片| 亚洲一区高清亚洲精品| 欧美三级亚洲精品| 丝袜美腿在线中文| 国产伦精品一区二区三区四那| 精品国产三级普通话版| 美女 人体艺术 gogo| 小说图片视频综合网站| 免费大片18禁| 亚洲在线观看片| 精品人妻1区二区| 久久午夜亚洲精品久久| 欧美绝顶高潮抽搐喷水| 美女高潮喷水抽搐中文字幕| 欧美高清性xxxxhd video| 久久久久久九九精品二区国产| 亚洲av熟女| 日日摸夜夜添夜夜添小说| 精品国产三级普通话版| 国产69精品久久久久777片| 色哟哟哟哟哟哟| 亚洲精品影视一区二区三区av| 中文字幕av在线有码专区| 搡女人真爽免费视频火全软件 | 久久久久久久午夜电影| 真实男女啪啪啪动态图| 男人舔女人下体高潮全视频| x7x7x7水蜜桃| 身体一侧抽搐| 一边摸一边抽搐一进一小说| 亚洲va日本ⅴa欧美va伊人久久| 亚洲一级一片aⅴ在线观看| 亚洲欧美精品综合久久99| 国产女主播在线喷水免费视频网站 | 在现免费观看毛片| 又黄又爽又免费观看的视频| 国产成人福利小说| 尤物成人国产欧美一区二区三区| 国产一区二区在线观看日韩| aaaaa片日本免费| 999久久久精品免费观看国产| 亚洲五月天丁香| 久久久久国产精品人妻aⅴ院| 俄罗斯特黄特色一大片| 中文字幕av在线有码专区| 国产亚洲91精品色在线| 精品人妻一区二区三区麻豆 | 免费看光身美女| 狂野欧美激情性xxxx在线观看| 大又大粗又爽又黄少妇毛片口| 国产私拍福利视频在线观看| 亚洲久久久久久中文字幕| 18禁黄网站禁片免费观看直播| 很黄的视频免费| 精品人妻偷拍中文字幕| 国产亚洲av嫩草精品影院| 又紧又爽又黄一区二区| 嫩草影院入口| 色综合色国产| 欧美zozozo另类| 色5月婷婷丁香| 狂野欧美激情性xxxx在线观看| 日日夜夜操网爽| 欧美黑人欧美精品刺激| 久久久久精品国产欧美久久久| 精品久久久久久久久亚洲 | 日韩欧美 国产精品| 久久天躁狠狠躁夜夜2o2o| 成人高潮视频无遮挡免费网站| 免费大片18禁| 免费av观看视频| 成人特级av手机在线观看| 日韩一区二区视频免费看| 亚洲国产欧美人成| 国产激情偷乱视频一区二区| 亚洲av成人av| 久久久精品大字幕| 高清日韩中文字幕在线| 一本精品99久久精品77| 久久精品夜夜夜夜夜久久蜜豆| 久久草成人影院| 午夜福利在线在线| 麻豆成人午夜福利视频| 亚洲成人免费电影在线观看| 久久天躁狠狠躁夜夜2o2o| 日韩精品中文字幕看吧| 免费看美女性在线毛片视频| 免费在线观看日本一区| 日韩欧美一区二区三区在线观看| 精品午夜福利视频在线观看一区| 男人舔奶头视频| 热99re8久久精品国产| 一区二区三区高清视频在线| 国产亚洲精品久久久com| 欧美日韩综合久久久久久 | 国产中年淑女户外野战色| 12—13女人毛片做爰片一| 免费观看的影片在线观看| 国语自产精品视频在线第100页| 国产高清有码在线观看视频| 日韩欧美在线乱码| 国国产精品蜜臀av免费| 国产欧美日韩精品一区二区| 免费高清视频大片| 真实男女啪啪啪动态图| 国产私拍福利视频在线观看| 亚洲精品乱码久久久v下载方式| 久久精品久久久久久噜噜老黄 | 我要看日韩黄色一级片| 一卡2卡三卡四卡精品乱码亚洲| 亚洲美女黄片视频| 啦啦啦啦在线视频资源| 国产男靠女视频免费网站| 人妻少妇偷人精品九色| .国产精品久久| 国产一区二区三区在线臀色熟女| 麻豆av噜噜一区二区三区| 中文字幕av在线有码专区| 精品人妻视频免费看| 制服丝袜大香蕉在线| 精品99又大又爽又粗少妇毛片 | 国产精品三级大全| 欧美一级a爱片免费观看看| 亚洲av美国av| 中文字幕久久专区| 琪琪午夜伦伦电影理论片6080| 少妇的逼水好多| or卡值多少钱| 国产免费男女视频| 国产高清有码在线观看视频| 成年版毛片免费区| 午夜免费激情av| 三级国产精品欧美在线观看| 熟女人妻精品中文字幕| 午夜福利高清视频| bbb黄色大片| 免费看a级黄色片| xxxwww97欧美| 91av网一区二区| 蜜桃久久精品国产亚洲av| 丰满的人妻完整版| 久久国产精品人妻蜜桃| 哪里可以看免费的av片| 国产精品电影一区二区三区| 亚洲无线在线观看| 99久国产av精品| 亚洲精品一区av在线观看|