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

    典型年份渭河流域藍水綠水時空差異分析*4

    2016-05-27 02:58:28趙安周朱秀芳潘耀忠河北工程大學資源學院邯鄲05608北京浙范大學地表過程與資源生態(tài)國家重點實驗室資源學院北京100875北京建工建筑設計研究院北京100044
    中國農(nóng)業(yè)氣象 2016年2期
    關鍵詞:渭河流域

    趙安周,朱秀芳,潘耀忠,劉 鑫(1. 河北工程大學資源學院,邯鄲 05608;. 北京浙范大學地表過程與資源生態(tài)國家重點實驗室/資源學院, 北京 100875;. 北京建工建筑設計研究院,北京 100044)

    ?

    典型年份渭河流域藍水綠水時空差異分析*4

    趙安周1,2,朱秀芳2**,潘耀忠2,劉鑫3
    (1. 河北工程大學資源學院,邯鄲 056038;2. 北京浙范大學地表過程與資源生態(tài)國家重點實驗室/資源學院, 北京 100875;3. 北京建工建筑設計研究院,北京 100044)

    摘要:以西北地區(qū)濕潤/干旱過渡地帶的渭河流域為例,采用降水距平百分率和標準化降水指數(shù)識別濕潤年、正常年和干旱年,同時結合SWAT(Soil and Water Assessment Tool)模型的輸出結果分析各典型年該流域藍水綠水的時空分布差異。結果表明:(1)采用降水距平百分率和標準化降水指數(shù)兩種方法確定1983、1997 和2007年分別為濕潤年、干旱年和正常年;(2)模型健價的結果表明,SWAT模型的月徑流模擬的精度較高,可以準確描述渭河流域的徑流變化過程;(3)在流域尺度上,綠水流的變化在濕潤年、干旱年和正常年的變化相對穩(wěn)定,綠水系數(shù)在濕潤年、干旱年和正常年分別為82.06%、93.47%和87.72%;(4)從空間分布來看,藍水綠水資源在渭河流域的分布呈現(xiàn)自東南向西北減少的趨勢,綠水系數(shù)的空間分布表明在濕潤年份(1983)或地區(qū)(東南部)綠水比重明顯低于干旱年份(1997)或地區(qū)(北部)。

    關鍵詞:SWAT模型;藍水綠水;綠水系數(shù);渭河流域

    趙安周,朱秀芳,潘耀忠,等.典型年份渭河流域藍水綠水時空差異分析[J].中國農(nóng)業(yè)氣象,2016,37(2):149-157

    藍水綠水的健價目前已成為水科學研究領域的熱點問題,并逐漸影響人類對水資源管理利用的方式,引發(fā)人們對水資源概念及其健價的重新思考[1-2]。藍水綠水的概念是20世紀90年代Falkenmark在健價半干旱地區(qū)水資源對農(nóng)業(yè)生產(chǎn)過程的影響中提出的,之后許多學者對其進行了進一步完善[3-4]。藍水為存儲于河流、湖泊和含水層中的水,是扣除地表水和地下水重復計算后的總和,綠水可以分為綠水流和綠水儲量,前者由植物散發(fā)和土壤、水體蒸發(fā)組成,后者是指存儲于土壤中的水,即土壤含水量[5]。藍水綠水概念的提出把生態(tài)學過程與水循環(huán)過程緊密聯(lián)系起來,體現(xiàn)了水文過程與植被之間的相互作用關系[6]。隨著藍水綠水概念的提出,其健價在水文水資源領域得到了高度重視。瑞典斯德哥爾摩國際水資源研究所(SIWI)、全球水系統(tǒng)項目組(GWSP)等眾多國際組織和機構開始致力于對綠水的研究。在國內(nèi),程國棟等[3, 7]指出綠水資源的重要性,建議把綠水資源納入水資源健價中。Liu等[8]采用GEPIC (GIS-based Environmental Policy Inte-grated Climate)模型在全球尺度上就主要農(nóng)作物對藍水綠水資源的影響進行了分析,認為在1998-2002年,這些農(nóng)作物耗水量達3823km3·a-1,綠水資源量比例超過80%;Zuo等[9-13]借助SWAT模型分別對渭河流域、黑河流域、盧氏流域和碧流河流域的藍水綠水資源進行了健價,取得了一系列卓有成效的研究結果。

    渭河流域位于中國西北干旱和濕潤過渡地帶,水資源短缺已成為限制地區(qū)社會經(jīng)濟發(fā)展的重要因子。目前對該地區(qū)水資源的研究主要集中在自然條件下藍水綠水資源時空分布上,針對典型年份藍水綠水的時空分布差異特征的研究相對較少。鑒于此,本文利用SWAT模型的輸出結果,輔以降水距平百分比和標準化降水指數(shù)等方法,分析典型年份條件下渭河流域藍水綠水資源的時空分布差異,以期為典型年份中國西北地區(qū)的水資源管理提供科學依據(jù)。

    1 資料與方法

    1.1研究區(qū)概況

    渭河是黃河最大的支流,發(fā)源于甘肅省定西市渭源縣的鳥鼠山,流經(jīng)天水、寶雞、西安、渭南等多個省市,最終在陜西省渭南市潼關縣注入黃河,干流全長818km。該流域位于104°00′-110°20′E,33°50′-37°18′N,海拔319-3929m,涇河和北洛河為其主要的兩個支流。渭河流域地勢西高東低,北部為黃土高原,南部為秦嶺山脈,冬季寒冷少雨,夏季炎熱多雨,年降雨量500~800mm,平均氣溫7.8~13.5 °C,蒸發(fā)較強烈,陸面蒸發(fā)量500mm,水面蒸發(fā)量660~1600mm[14]。該流域是陜西省社會經(jīng)濟最發(fā)達的地區(qū),集中全省64%的人口、56%的耕地、72%的灌溉面積及其80%的經(jīng)濟生產(chǎn)總值[15],流域位置及氣象、水文站點分布見圖1。

    圖1 渭河流域及其支流位置、氣象站點和水文站點分布Fig. 1 Location and distribution of meteorological stations and hydrological stations in Weihe River Basin

    1.2數(shù)據(jù)來源

    渭河流域及其周邊22個氣象站點的降水、氣溫、相對濕度、風速數(shù)據(jù)、日照時數(shù)等數(shù)據(jù)來源于氣象數(shù)據(jù)共享網(wǎng),為1970-2012年的日值數(shù)據(jù)。藍水為藍水流,包括河流、湖泊和含水層中的水,綠水包括實際蒸散量和土壤含水量,均來源于SWAT模型模擬的結果。建立SWAT模型需要的DEM(Digital Elevation Model)數(shù)據(jù)(90m×90m)、土地利用/覆蓋數(shù)據(jù)(1:100萬)和土壤數(shù)據(jù)(1:10萬)分別來源于國際科學數(shù)據(jù)共享平臺、寒區(qū)旱區(qū)科學數(shù)據(jù)中心(WestDC)和世界土壤數(shù)據(jù)庫(HWSD)。華縣、咸陽、林家村和魏家堡4個水文站的月流量數(shù)據(jù)來源于黃河水利委員會和中華人民共和國水文年鑒,其時間段為1971-1990年。

    1.3研究方法

    1.3.1典型年份的確定

    利用降水距平百分率Pa和標準化降水指數(shù)SPI兩個指標確定典型年份,以避免采用單一指標的計算誤差,其具體計算過程參見文獻[16-17],劃分標準見表1[1]。當Pa和SPI的等級≥5,定義該年份為干旱年,等級值均為4,定義為正常年份,當?shù)燃墶?,定義為濕潤年。

    計算1972-2012年22個氣象站點的年降水距平百分率(Pa)和標準化降水指數(shù)(SPI),各流域的值為流域內(nèi)所含氣象站點的平均值,按表1劃分其等級值。

    表1 降水距平百分率(Pa)和標準化降水指數(shù)(SPI)的旱澇等級Table 1 Grades of drought/flood based on the anomalous percentage and standard precipitation index

    1.3.2綠水系數(shù)計算

    綠水系數(shù)是指某一地區(qū)或流域綠水所占的比重,具體計算方法為[18]

    式中,GWC(Green Water Coefficient)表示綠水系數(shù),G表示綠水,包括綠水流和綠水儲量,B表示藍水。

    1.3.3SWAT模型

    SWAT模型是一個分布式水文模型,該模型是美國農(nóng)業(yè)部農(nóng)業(yè)研究局在CREAMS(Chemicals, Runoff, and Erosion from Agricultural Management Systems)和SWRRB(Simulator for Water Resources in Rural Basins)模型的基礎上發(fā)展起來的,其建立模型所需的數(shù)據(jù)包括地形、土壤、土地利用、氣象、水文等。模型可以依據(jù)DEM數(shù)據(jù)、土地利用數(shù)據(jù)、土壤數(shù)據(jù)等將研究區(qū)域劃分為多個只包含一種土地利用類型和土壤類型的水文響應單元(HRU)。由于該模型可以模擬流域水循環(huán)中幾乎所有的水文過程,因此,在全球許多國家和地區(qū)得到廣泛應用[19-20]。本文通過建立渭河流域的土地利用/覆蓋數(shù)據(jù)庫、土壤數(shù)據(jù)庫、天氣發(fā)生器等,利用SWAT模型的水文模塊對渭河流域的地表水、蒸散、地下水、土壤水等水文過程要素進行模擬。藍綠水是根據(jù)SWAT模型輸出結果進行計算,其中,藍水流為地表徑流、壤中流和地下徑流之和,綠水流為實際蒸散量(ET),綠水儲量為土壤含水量。

    1.3.4模型健價指標

    模型的率定和驗證采用SWAT-CUP里的SUFI-2算法,該算法采用P-factor和R-factor對模型的不確定性進行定量健估。P-factor表示實測數(shù)據(jù)落入置信區(qū)間95PPU(95% prediction uncertainty)的概率百分比,R-factor表示 95PPU 帶的平均厚度除以監(jiān)測數(shù)據(jù)標準偏差。一般P-factor越接近1,R-factor越接近0,表明模擬的效果越好[21]。同時,本文采用決定系數(shù)R2、Nash-Sutcliffe 效率系數(shù)(ENS)來健價SWAT模型的適用性,其具體計算參照文獻[22-23]。其中,R2表示模擬值與實測值變化趨勢的一致性,其值越接近1,表明模擬值與實測值趨勢越吻合,通常將R2>0.6作為實測值和模擬值相關程度的標準;ENS表示實測值與模擬值的偏離程度,當ENS≤0.36時,認為模擬效果不好,0.36<ENS<0.75,認為模擬效果令人滿意,ENS≥0.75,認為模擬的結果好[24]。

    2 結果與分析

    2.1流域典型干濕年份的確定

    各流域1972-2012年年降水距平百分率(Pa)和標準化降水指數(shù)(SPI)等級值見圖2。利用SPSS對兩種方法的分級結果的一致性進行檢驗,由表2可見,各流域的顯著性檢驗的P值均大于0.05,表明Pa和SPI在衡量渭河流域的旱澇程度上無顯著性差異,二者具有良好的一致性。根據(jù)圖2和表2,選取1983、1997和2007年分別作為渭河流域、渭河干流、涇河流域和北洛河流域的濕潤年、干旱年和正常年。

    2.2利用典型站月徑流資料率定、驗證SWAT模型參數(shù)

    圖2 各區(qū)歷年降水距平百分率(Pa)和SPI指數(shù)的等級值(1972-2012)Fig. 2 Rank values of precipitation anomaly index (Pa) and SPI in each region(1972-2012)

    利用渭河流域華縣、咸陽、魏家堡、林家村4個典型水文站1981-1990年和1971-1980年的月徑流資料對模型進行率定和驗證。其月徑流模擬結果見圖3,相應的不確定性的分析及模擬精度健價指標結果見表3。從圖中可見,各站率定期和驗證期月徑流量模擬值與實測值的波動趨勢均基本一致,進一步分析相關指標(表3)可見,率定期內(nèi)流域下游和中游(1981-1990年)的華縣站和咸陽站月徑流的P-factor分別為0.85和0.76,R-factor分別為0.57和0.40,驗證期間(1971-1980年)華縣站和咸陽站的P-factor分別為0.82和0.73,R-factor分別為0.67和0.45,表明觀測值落在95PPU的概率較大,其徑流模擬的不確定性較??;魏家堡和林家村率定期(1981-1990年)的P-factor分別為0.68和0.65,R-factor分別為0.52和0.68,其驗證期(1971-1980年)的P-factor分別為0.66和0.61,R-factor分別為0.85和1.05,與華縣站和咸陽站相比,這兩個站不確定性增大,從健價指標來看,4個站點在率定期和驗證期的R2均大于0.60,ENS系數(shù)除了在驗證期的林家村站外,其它站點的值均大于0.60,模擬效果較好。林家村水文站位于寶雞峽附近,受到其農(nóng)業(yè)灌溉用水的影響,其不確定性較大,模擬精度較差,華縣站和咸陽站的R2和ENS的值均較大,率定期的值在0.80以上,驗證期也接近0.70,模擬精度較高。可見,從4個水文站點的健價指標的數(shù)值來看,其徑流模擬的精度較高,可以基本描述渭河流域模擬時段的徑流變化。

    2.3利用SWAT模型模擬分析流域典型年藍綠水時空特征

    2.3.1藍綠水流比較

    渭河流域、渭河干流、涇河流域和北洛河流域典型年份的藍水綠水資源如圖4所示。從圖可以看出,渭河流域的藍水流在濕潤年、干旱年及正常年分別為123.21、30.92、67.56mm·a-1。渭河干流、涇河流域和北洛河流域的藍水流在濕潤年均高于干旱年(圖4),這主要是由于濕潤年的降水明顯高于干旱年所致;綠水流的變化在3個典型年份的變化較藍水流穩(wěn)定;渭河流域的綠水儲量在濕潤年、干旱年及正常年分別為92.51、43.93和62.59mm·a-1,濕潤年的綠水儲量大于干旱年。從綠水系數(shù)來看,渭河干流、涇河流域和北洛河流域的綠水系數(shù)(1983、1997和2007年的均值)分別為85.66%、91.31%和91.57%,可見,北洛河的綠水系數(shù)最高;渭河流域的綠水系數(shù)在濕潤年、干旱年和正常年分別為82.06%、93.47%和87.72%,干旱年的綠水系數(shù)明顯高于濕潤年,表明典型干旱年的綠水資源占水資源的比例明顯高于濕潤年,對于干旱年份,綠水資源尤其重要。

    圖3 渭河流域4個水文站點月徑流的率定期和驗證期的模擬結果Fig. 3 Comparison between the measured and simulated monthly discharges at four hydrological stations for both calibration and validation period

    表3 渭河流域主要水文站月徑流模擬結果評價和不確定性分析Table 3 Results of monthly runoff simulation for four main hydrological stations

    圖4 渭河流域及其支流流域在典型年份的降水量、藍水綠水資源和綠水系數(shù)Fig. 4 Precipitation, blue and green water and green water coefficient in typical reference years in the Weihe River Basin, Weihe River mainstream, Jinghe River Basin and Beiluohe River Basin

    2.3.2藍水綠水資源的空間分布

    典型年份渭河流域藍水流的空間變化如圖5a所示??傮w來看,藍水流在渭河流域的空間分布呈現(xiàn)自東南向西北減少的趨勢,濕潤年流域北部多數(shù)地區(qū)的藍水流在50~150mm·a-1,南部的藍水流大于150mm·a-1;正常年份大部分地區(qū)的藍水流大于50mm·a-1,南部秦嶺山脈和關中平原的個別地區(qū)大于200mm·a-1,僅流域西北部小于50mm·a-1;干旱年流域上游和涇河流域的藍水流基本均小于50mm·a-1,北洛河流域和渭河流域的中下游的藍水流界于51~150mm·a-1。

    由圖5b可見,綠水流的空間分布格局與藍水流相似,亦呈現(xiàn)東南向西北減少的趨勢,濕潤年大部地區(qū)的綠水流處于550mm·a-1以上,下游部分地區(qū)甚至超過700mm·a-1;其次為正常年,流域西北部的綠水流最小,大部地區(qū)小于400mm·a-1,東南部的綠水流相對豐富,大部地區(qū)界于550~700mm·a-1;干旱年的綠水流最匱乏,流域上游和北部的綠水流基本處于400mm·a-1以下,僅東部的少數(shù)地區(qū)大于550mm·a-1(圖5b)。

    綠水儲量的空間分布如圖5c所示。由圖可見,典型年份的綠水儲量呈現(xiàn)自東南向西北減少的趨勢,濕潤年的綠水儲量最豐富,流域東南部和中部地區(qū)大于90mm·a-1,西北部也基本處于30~70mm·a-1;正常年的綠水儲量以北洛河流域和關中地區(qū)最豐富,西部北的天水、固原地區(qū)最匱乏;干旱年的綠水儲量在整個流域均較匱乏,全流域基本小于70mm·a-1。

    2.3.3綠水系數(shù)的空間分布

    渭河流域的綠水系數(shù)的空間分布如圖6所示。從圖可以看出,其空間分布總體表現(xiàn)為北部地區(qū)高于南部。濕潤年份綠水系數(shù)在渭河流域的個別地區(qū)大于90%,在流域南部的秦嶺山脈地區(qū)小于65%;干旱年份整個流域大部地區(qū)的綠水系數(shù)均在90%以上。干旱年份氣溫較高,其蒸散發(fā)消耗的水資源高于濕潤年份,導致干旱年份的綠水系數(shù)高于濕潤年份。另外,從典型年份綠水系數(shù)的平均值來看,渭河流域87.75%的水資源是以綠水的形式參與到流域水循環(huán)當中的。

    圖5 典型年份渭河流域藍水綠水的時空分布格局Fig. 5 Spatial distributions of blue water flow, green water flow and green water storage in typical reference years in the Weihe River Basin

    圖6 典型年份水文單元(HRU)尺度上渭河流域綠水系數(shù)的時空分布格局Fig. 6 Spatial distribution of green water coefficient in typical reference years in the Weihe River Basin

    3 結論與討論

    (1)采用降水距平百分率Pa和標準化降水SPI指數(shù),確定1983年、1997年和2007年分別為渭河流域的典型濕潤年、干旱年和正常年。

    (2)渭河流域4個水文站點(華縣、咸陽、林家村和魏家堡)的月徑流的模擬值與實測值總體模擬效果較好。上游的林家村和魏家堡水文站的決定系數(shù)R2和ENS的值較其它兩個站點低,華縣站和咸陽站的率定期的R2和ENS的值均大于0.8,驗證期也在0.70左右,表明模型的月徑流模擬的精度較高,可以準確描述模擬期內(nèi)渭河流域的徑流變化過程。

    (3)從典型年份來看,渭河流域及其各支流的藍水綠水資源存在明顯的差別,綠水流的變化相對于藍水流較穩(wěn)定,綠水系數(shù)在濕潤年、干旱年和正常年分別為82.06%、93.47%和87.72%。

    (4)從空間分布來看,藍水綠水資源呈現(xiàn)自東南向西北減少的趨勢,藍水流的時空分布格局與降水相似,其空間分布由降水控制,綠水資源不僅與降水等氣候條件有關,而且與人類活動和下墊面的性質有關;就典型年份平均值來看,綠水系數(shù)的空間分布總體表現(xiàn)為涇河流域中部、北洛河流域中下游及其渭河干流的下游地區(qū)最高。

    渭河流域橫跨秦嶺山脈、六盤山、黃土高原、關中平原等多種地形地貌,其降水分布格局受地形和大氣環(huán)流的影響不同,使不同地區(qū)呈現(xiàn)的干濕年份特性不一致;從空間分布來看,典型年份渭河流域的藍水流呈現(xiàn)從東南向西北逐漸減小的趨勢,主要是由于流域的降水分布格局導致[25];綠水資源空間分布除與降水分布有關外,還與溫度、土壤類型等下墊面條件有關。從綠水系數(shù)的平均值來看,渭河流域藍水綠水總資源量的87%以綠水形式參與到流域的水循環(huán)中,這與中國東部地區(qū)的一些流域的結果有所不同,如東南沿海的東江流域的多年平均綠水系數(shù)僅39%[19],晉江西溪流域1973-1979年的平均綠水系數(shù)為36%[26],這些流域處于亞熱帶季風地區(qū),降水較多,蓄滿產(chǎn)流為主要的產(chǎn)流方式,降雨可以更多地轉化為地表徑流,使這些流域的藍水資源高于綠水資源量,而渭河流域屬大陸性氣候,氣候干燥,降水較少,多數(shù)降水通過植被截留和下滲轉換為生態(tài)用水,產(chǎn)生的地表徑流較少,對同處于中國西部地區(qū)的黑河流域[1]、湟水流域[27]、盧氏河流域[12]的研究也得到與本文相似的結論,其綠水資源量遠大于藍水資源。從空間分布來看,典型年份的藍水流表現(xiàn)為自東南向西北減少的趨勢,降水是影響其空間分布的主要因素;影響綠水資源,尤其是綠水流的因素較多,除氣候、下墊面性質等自然因素外,人類活動(退耕還林還草、修建水庫和淤地壩、農(nóng)業(yè)灌溉等)也會影響其空間分布格局。但本文對對藍水綠水進行模擬時未考慮農(nóng)業(yè)灌溉、修建水庫等人類活動對典型年份藍水綠水的影響,同時對可以緩解流域水資源和生態(tài)環(huán)境壓力的藍水綠水的轉換機理問題尚不明確,這些均需以后作進一步研究。

    參考文獻References

    [1]臧傳富.黑河流域藍綠水時空變化研究[D].北京:北京林業(yè)大學,2013.

    Zang C F.Assessment of spatial and temporal variability of green and blue water flows in the Heihe River Basin,China[D]. Beijing:Beijing Forestry University,2013.(in Chinese)

    [2]Rockstr?m J,Karlberg L,Wani S P,et al.Managing water in rainfed agriculture:the need for a paradigm shift[J]. Agriculture Water Manage,2010,97(4):543-550.

    [3]程國棟,趙文智.綠水及其研究進展[J].地球科學進展,2006, 21(3):221-227.

    Cheng G D,Zhao W Z.Green water and its research progresses[J].Advances in Earth Science,2006,21(3):221-227. (in Chinese)

    [4]Falkenmark M,Rockstr?m J.The new blue and green water paradigm:breaking new ground for water resources planning and management[J].Journal of water resources planning and management,2006,132(3):129-132.

    [5]Falkenmark M.Land-water linkages:a synopsis[J].FAO Land and Water Bulletin,1995(1):15-16.

    [6]臧傳富,劉俊國.黑河流域藍綠水在典型年份的時空差異特征[J].北京林業(yè)大學學報,2013,35(3):1-10.

    Zang C F,Liu J G.Spatial and temporal pattern difference of blue-green water flows in typical reference years in the Heihe River Basin,northwestern China[J].Journal of Beijing Forestry University,2013,35(3):1-10.(in Chinese)

    [7]李小雁.流域綠水研究的關鍵科學問題[J].地球科學進展,2008,23(7):707-712.

    Li X Y.Key scientific issues for green water research in the watershed[J].Advances in Earth Science,2008,23(7):707-712. (in Chinese)

    [8]Liu J G,Zehnder A,Yang H.Global consumptive water use for crop production:the importance of green water and virtual water[J].Water Resources Research,2009,45(5):W05428, doi: 10.1029/2007WR006051.

    [9]Zuo D,Xu Z,Peng D,et al.Simulating spatiotemporal variability of blue and green water resources availability withuncertainty analysis[J].Hydrological Processes,2014,29(8): 1942-1955.

    [10]Zang C F,Liu J,Velde M,et al.Assessment of spatial and temporal patterns of green and blue water flows under natural conditions in inland river basins in Northwest China[J].Hydrology And Earth System Sciences,2012,16(8): 2859-2870.

    [11]甄婷婷,徐宗學,程磊,等.藍水綠水資源量估算方法及時空分布規(guī)律研究[J].資源科學,2010,32(6):1177-1183.

    Zhen T T,Xu Z X,Cheng L,et al.Spatiotemporal distributions of blue and green water resources[J].Resources Science,2010, 32(6):1177-1183.(in Chinese)

    [12]徐宗學,左德鵬.拓寬思路,科學健價水資源量:以渭河流域藍水綠水資源量健價為例[J].南水北調與水利科技,2013, 11(1):12-16.

    Xu Z X,Zuo D P.Scientific assessment of water resources with broaden thoughts:a case study on the blue and green water resources in the Wei River Basin[J].South-to-north Water Transfers and Water Science & Technology, 2013,11(1):12-16.(in Chinese)

    [13]吳洪濤,武春友,郝芳華,等.綠水的多角度健估及其在碧流河上游地區(qū)的應用[J].資源科學,2009,31(3):420-428.

    Wu H T,Wu C Y,Hao F H,et al.Assessment of green water from multi-angle view at catchment scale[J].Resources Science,2009,31(3):420-428. (in Chinese)

    [14]左德鵬,徐宗學,李景玉,等.氣候變化情景下渭河流域潛在蒸散量時空變化特征[J].水科學進展,2011,22(4):455-461.

    Zuo D P,Xu Z X,Li J Y,et al.Spatiotemporal characteristics of potential evapotranspiration in the Weihe River basin under future climate change[J].Advances in Water Science,2011, 22(4): 455-461.(in Chinese)

    [15]葛芬莉.關中地區(qū)渭河流域水資源與水環(huán)境綜合治理研究[J].西北水力發(fā)電,2004,20(S1):160-162.

    Ge F L.Study of the water resources and water environment comprehensive treatment in the Guanzhong area of the Weihe River Basin[J].Journal of Northwest Hydroelectric Power,2004, 20(S1):160-162.(in Chinese)

    [16]鞠笑生,楊賢為.我國單站旱澇指標確定和區(qū)域旱澇級別劃分的研究[J].應用氣象學報,1997,8(1):26-33.

    Ju X S,Yang X W.Research on determination of station indexes and division of regional flood /drought grades in China[J].Quarterly Journal of Applied Meteorology,1997, 8(1):26-33.(in Chinese)

    [17]McKee T B,Doesken N J,Kleist J.The relationship of drought frequency and duration to time scales[A].Proceedings of the 8th conference on applied climatology[C].Boston:American Meteorological Society Boston, 1993:179-183.

    [18]呂樂婷.東江流域氣候變化及人類活動對水文水資源的影響[D].北京:北京浙范大學,2014.

    Lv L T.Impacts of climate change and human activities on hydrology and water resources in the Dongjiang River Basin[D].Beijing:Beijing Normal University,2014.(in Chinese) [19]Naramngam S,Tong S T.Environmental and economic implications of various conservative agricultural practices in the Upper Little Miami River Basin[J].Agricultural Water Management,2013,119:65-79.

    [20]Faramarzi M,Abbaspour K C,Schulin R,et al.Modelling blue and green water resources availability in Iran[J]. Hydrological Processes,2009,23(3):486-501.

    [21]Sun C,Ren L.Assessment of surface water resources and evapotranspiration in the Haihe River basin of China using SWAT model[J].Hydrological Processes,2013,27(8):1200-1222.

    [22]Krause P,Boyle D P,B?se F.Comparison of different efficiency criteria for hydrological model assessment[J]. Advances in Geosciences,2005,5(5):89-97.

    [23]Freer J,Beven K,Ambroise B.Bayesian estimation of uncertainty in runoff prediction and the value of data:an application of the GLUE approach[J].Water Resources Research,1996,32(7):2161-2173.

    [24]Geza M,McCray J E.Effects of soil data resolution on SWAT model stream flow and water quality predictions[J].Journal of Environmental Management,2008,88(3):393-406.

    [25]劉燕,胡安焱.渭河流域近50年降水特征變化及其對水資源的影響[J].干旱區(qū)資源與環(huán)境,2006,20(1):85-87.

    Liu Y,Hu A Y.Changes of precipitation characters along Weihe Basin in 50 years and its influence on water resources[J].Journal of Arid Land Resources and Environment,2006,20(1): 85-87.(in Chinese)

    [26]榮琨,陳興偉,李志遠,等.晉江西溪流域綠水藍水資源量估算及分析[J].水土保持通報,2011,31(4):12-15.

    Rong K,Chen X W,Li Z Y,et al.Estimating green and blue water resources in Xixi Watershed of Jinjiang Basin[J]. Bulletin of Soil and Water Conservation,2011,31(4):12-15. (in Chinese)

    [27]張杰,賈紹鳳.基于SWAT模型的湟水流域藍綠水與不同土地利用類型的綠水差異研究[J].水資源與水工程學報,2013, 24(4):6-10.

    Zhang J,Jia S F.Study on difference of blue-green water in Huangshui basin and green water under different types of land use based on SWAT model[J].Journal of Water Resources and Water Engineering,2013,24(4):6-10.(in Chinese)

    Spatiotemporal Variation of Blue-Green Water under Typical Reference Years in the Weihe River Basin

    ZHAO An-zhou1,2, ZHU Xiu-fang2, PAN Yao-zhong2, LIU Xin3
    (1. College of Resources, Hebei University of Engineering, Handan 056038, China; 2.State Key Laboratory of Earth Surface Processes and Resource Ecology/College of Resources Science & Technology, Beijing Normal University, Beijing 100875; 3.Beijing Jiangong Building Design and Research Institute, Beijing 100044)

    Abstract:In order to understand the spatiotemporal pattern difference of blue-green water in typical reference years on a watershed scale, this study selected the Weihe River Basin in northwestern of China to analyze the spatiotemporal distribution pattern difference of blue-green water, and obtained the typical reference years by using standardized precipitation index(SPI) and precipitation anomaly index(Pa). The results showed that: (1)the years of 1983, 1997 and 2007 were determined as wet year, dry year and normal year, respectively. (2)Calibration and validation results showed a similar change trend between simulated and observed discharges, and the SWAT can be used to predict hydrological processes in Weihe River Basin. (3)Change of the green water was relatively stable in the typical reference years, and the green water coefficients were 82.06%, 93.47% and 87.72% in the wet year, dry year and normal year, respectively. (4)Spatially, the blue-green water decreased from southeast to northwest in Weihe River Basin. The spatial distribution of green water indicated that the proportion of green water in the wet year (1983) or areas (southeast) was lower than the normal and dry year (1997) or areas (north).

    Key words:SWAT; Blue-green water; Green water coefficient; Weihe River Basin

    doi:10.3969/j.issn.1000-6362.2016.02.004

    * 收稿日期:2015-06-30**通訊作者。E-mail: zhuxiufang@bnu.edu.cn

    基金項目:國家自然科學基金青年科學基金項目(41401479);國家“高分辨率對地觀測系統(tǒng)”重大專項

    作者簡介:趙安周(1985-),博士,講浙,主要研究方向為水資源對氣候變化和土地利用的響應及其干旱健價。E-mail: zhaoanzhou@126.com

    猜你喜歡
    渭河流域
    渭河水岸帶土壤理化特征及其與重金屬含量的關系
    基于NDVI的渭河流域時空演變分析
    甘肅省渭河流域韭菜產(chǎn)業(yè)現(xiàn)狀及發(fā)展建議
    基于DPCA方法對渭河流域降水—徑流—蒸發(fā)序列的長程互相關性研究
    世界家苑(2018年1期)2018-04-27 11:42:06
    基于RS/GIS 渭河流域植被覆蓋時空變化特征研究
    燕太子回國
    環(huán)境流域治理中的法律對策
    法制博覽(2015年10期)2015-11-06 12:16:49
    渭河流域香菜夏秋無公害栽培技術
    生態(tài)系統(tǒng)服務支付意愿及其影響因素分析
    軟科學(2015年6期)2015-07-10 02:38:40
    渭河流域糧食生產(chǎn)與灌溉農(nóng)業(yè)發(fā)展的相互作用關系分析
    中國水利(2015年5期)2015-02-28 15:12:36
    欧美乱色亚洲激情| 51午夜福利影视在线观看| 亚洲成人精品中文字幕电影| 波多野结衣av一区二区av| 欧美+亚洲+日韩+国产| 亚洲一区二区三区色噜噜| 日韩中文字幕欧美一区二区| 欧美乱色亚洲激情| 欧美性长视频在线观看| 香蕉久久夜色| 久久 成人 亚洲| 久久午夜综合久久蜜桃| 女人精品久久久久毛片| 国产欧美日韩综合在线一区二区| 亚洲精品中文字幕一二三四区| 久久中文字幕人妻熟女| 一个人免费在线观看的高清视频| 午夜福利,免费看| 国产欧美日韩精品亚洲av| 久久久久久久午夜电影| 午夜免费观看网址| 变态另类丝袜制服| 十八禁网站免费在线| 欧美黄色片欧美黄色片| 十八禁网站免费在线| 一二三四社区在线视频社区8| 波多野结衣av一区二区av| 男女午夜视频在线观看| 国产不卡一卡二| 欧美日本中文国产一区发布| 久久久久久大精品| 久久中文字幕人妻熟女| 国产成人精品久久二区二区91| 亚洲一码二码三码区别大吗| 亚洲九九香蕉| 国产成人精品久久二区二区91| 叶爱在线成人免费视频播放| 国产亚洲精品久久久久5区| 久久香蕉国产精品| 亚洲精品美女久久久久99蜜臀| 久久香蕉国产精品| 色综合站精品国产| 免费高清在线观看日韩| 亚洲成av片中文字幕在线观看| 免费在线观看完整版高清| 国产精品久久久久久精品电影 | 精品卡一卡二卡四卡免费| 级片在线观看| 亚洲 欧美一区二区三区| 中文字幕高清在线视频| 久久久久久久久中文| 激情在线观看视频在线高清| 久久国产精品影院| 国产精品久久久久久亚洲av鲁大| 身体一侧抽搐| 在线观看66精品国产| 老汉色∧v一级毛片| 亚洲成人久久性| 999久久久国产精品视频| 欧美在线黄色| 国产精品香港三级国产av潘金莲| 波多野结衣高清无吗| 此物有八面人人有两片| 日本免费一区二区三区高清不卡 | a在线观看视频网站| 久久性视频一级片| 色精品久久人妻99蜜桃| 免费高清视频大片| 欧美午夜高清在线| 99热只有精品国产| 色老头精品视频在线观看| 怎么达到女性高潮| 校园春色视频在线观看| 精品熟女少妇八av免费久了| 好男人在线观看高清免费视频 | 中亚洲国语对白在线视频| 国产欧美日韩一区二区精品| 黄色片一级片一级黄色片| 国内久久婷婷六月综合欲色啪| 日本三级黄在线观看| 欧美一级a爱片免费观看看 | 免费女性裸体啪啪无遮挡网站| 精品午夜福利视频在线观看一区| 欧美另类亚洲清纯唯美| 在线播放国产精品三级| 国产1区2区3区精品| 国产精品亚洲美女久久久| 又黄又爽又免费观看的视频| 亚洲人成网站在线播放欧美日韩| 99久久久亚洲精品蜜臀av| 欧美成人午夜精品| 女人被躁到高潮嗷嗷叫费观| 国产野战对白在线观看| 免费少妇av软件| 在线永久观看黄色视频| 在线免费观看的www视频| 在线免费观看的www视频| 国产一区二区激情短视频| 国产av精品麻豆| 亚洲精品中文字幕一二三四区| cao死你这个sao货| 大香蕉久久成人网| 午夜精品久久久久久毛片777| 两人在一起打扑克的视频| √禁漫天堂资源中文www| 国产精品一区二区在线不卡| 国产在线精品亚洲第一网站| 国产1区2区3区精品| 一本综合久久免费| 亚洲色图av天堂| 大型黄色视频在线免费观看| 国产蜜桃级精品一区二区三区| 亚洲狠狠婷婷综合久久图片| 国产精品久久视频播放| 99久久国产精品久久久| 亚洲精品国产色婷婷电影| 中文字幕高清在线视频| 午夜激情av网站| 女人被躁到高潮嗷嗷叫费观| 一a级毛片在线观看| АⅤ资源中文在线天堂| 欧美成人性av电影在线观看| 午夜福利在线观看吧| 黄色片一级片一级黄色片| 亚洲男人的天堂狠狠| 久久久久久久久中文| 999久久久精品免费观看国产| 亚洲美女黄片视频| 美女高潮到喷水免费观看| 午夜亚洲福利在线播放| 久久这里只有精品19| 色尼玛亚洲综合影院| 久久国产亚洲av麻豆专区| 国产精品亚洲一级av第二区| 又黄又爽又免费观看的视频| 日日夜夜操网爽| 亚洲专区国产一区二区| 韩国精品一区二区三区| 日韩欧美一区二区三区在线观看| 日韩三级视频一区二区三区| 嫩草影院精品99| 老熟妇乱子伦视频在线观看| 麻豆国产av国片精品| 男人操女人黄网站| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲人成网站在线播放欧美日韩| 亚洲精品av麻豆狂野| 欧美日韩瑟瑟在线播放| 午夜免费成人在线视频| 欧美老熟妇乱子伦牲交| 男女下面插进去视频免费观看| 99精品久久久久人妻精品| 人人澡人人妻人| 一级黄色大片毛片| avwww免费| 人成视频在线观看免费观看| 国产激情久久老熟女| 欧美乱色亚洲激情| 91字幕亚洲| 亚洲精品粉嫩美女一区| 黄片大片在线免费观看| 婷婷丁香在线五月| 一级毛片精品| 午夜激情av网站| 亚洲av五月六月丁香网| 国产一区二区三区视频了| 久热这里只有精品99| 91av网站免费观看| 可以在线观看的亚洲视频| 精品久久蜜臀av无| 国产精品98久久久久久宅男小说| 亚洲人成77777在线视频| 亚洲avbb在线观看| 亚洲av成人不卡在线观看播放网| 免费看十八禁软件| 中文字幕高清在线视频| av欧美777| 午夜精品在线福利| 欧美日韩瑟瑟在线播放| 91成年电影在线观看| 国产精品久久电影中文字幕| 97碰自拍视频| 男人的好看免费观看在线视频 | 亚洲人成电影免费在线| 国产精品98久久久久久宅男小说| 久久久久久国产a免费观看| 欧美成人免费av一区二区三区| 色综合婷婷激情| 亚洲 国产 在线| 99国产精品一区二区蜜桃av| 久久 成人 亚洲| 久久狼人影院| 亚洲色图 男人天堂 中文字幕| 女人爽到高潮嗷嗷叫在线视频| 精品人妻在线不人妻| 嫁个100分男人电影在线观看| 一级毛片高清免费大全| av在线播放免费不卡| 欧美成狂野欧美在线观看| 好看av亚洲va欧美ⅴa在| 婷婷六月久久综合丁香| 久久久久久亚洲精品国产蜜桃av| 日韩中文字幕欧美一区二区| 亚洲av熟女| 亚洲情色 制服丝袜| 午夜免费成人在线视频| a在线观看视频网站| 欧美日本视频| 欧美激情 高清一区二区三区| 国产精品秋霞免费鲁丝片| 777久久人妻少妇嫩草av网站| 免费高清视频大片| 国产精品永久免费网站| 村上凉子中文字幕在线| 精品熟女少妇八av免费久了| 成年版毛片免费区| 久久久精品欧美日韩精品| 亚洲欧洲精品一区二区精品久久久| 国产精品综合久久久久久久免费 | 久久国产亚洲av麻豆专区| 国产精品香港三级国产av潘金莲| 久久这里只有精品19| 少妇 在线观看| 亚洲无线在线观看| 国内精品久久久久精免费| 国产又爽黄色视频| 国产麻豆69| 黑人巨大精品欧美一区二区mp4| 亚洲性夜色夜夜综合| av天堂久久9| 亚洲少妇的诱惑av| 久久久国产成人免费| 成人国产综合亚洲| 国产精品自产拍在线观看55亚洲| 一级毛片女人18水好多| 亚洲成人国产一区在线观看| 亚洲avbb在线观看| 久久国产乱子伦精品免费另类| 国产精华一区二区三区| 国产精品免费视频内射| а√天堂www在线а√下载| 欧美一级a爱片免费观看看 | 亚洲国产精品sss在线观看| 两性夫妻黄色片| 天天一区二区日本电影三级 | 欧美性长视频在线观看| 国产精品自产拍在线观看55亚洲| 中文字幕av电影在线播放| 日韩欧美一区二区三区在线观看| 成熟少妇高潮喷水视频| 涩涩av久久男人的天堂| 91国产中文字幕| 搞女人的毛片| 国产精品精品国产色婷婷| 国产在线精品亚洲第一网站| 他把我摸到了高潮在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产三级在线视频| 99在线视频只有这里精品首页| 亚洲一区中文字幕在线| 丝袜美腿诱惑在线| 亚洲中文日韩欧美视频| 亚洲精品粉嫩美女一区| 黄色a级毛片大全视频| 国产真人三级小视频在线观看| 在线观看免费午夜福利视频| 日本黄色视频三级网站网址| 1024视频免费在线观看| 激情在线观看视频在线高清| 一夜夜www| 亚洲色图综合在线观看| 亚洲第一青青草原| 美女扒开内裤让男人捅视频| 亚洲国产精品999在线| 高清毛片免费观看视频网站| 窝窝影院91人妻| 老司机在亚洲福利影院| 免费观看精品视频网站| 超碰成人久久| 日韩精品免费视频一区二区三区| 18禁美女被吸乳视频| 大香蕉久久成人网| 波多野结衣一区麻豆| 亚洲色图综合在线观看| av欧美777| 给我免费播放毛片高清在线观看| 免费观看精品视频网站| 在线观看www视频免费| 亚洲av美国av| 色综合欧美亚洲国产小说| 亚洲 国产 在线| 91国产中文字幕| 国产精品99久久99久久久不卡| 成在线人永久免费视频| 9191精品国产免费久久| 变态另类丝袜制服| 人人澡人人妻人| 日韩欧美一区二区三区在线观看| 日本三级黄在线观看| 国产成人系列免费观看| 少妇的丰满在线观看| 日韩免费av在线播放| 给我免费播放毛片高清在线观看| 亚洲中文字幕日韩| 久久久久久久久久久久大奶| 激情在线观看视频在线高清| 人成视频在线观看免费观看| 亚洲午夜精品一区,二区,三区| 精品国产乱码久久久久久男人| 成人精品一区二区免费| 99国产精品免费福利视频| 国产精品永久免费网站| 久久精品国产亚洲av香蕉五月| 国产成人欧美在线观看| 日韩欧美国产一区二区入口| 非洲黑人性xxxx精品又粗又长| 美女国产高潮福利片在线看| 午夜影院日韩av| 老司机深夜福利视频在线观看| 欧美日本亚洲视频在线播放| 午夜亚洲福利在线播放| 亚洲午夜理论影院| 90打野战视频偷拍视频| 99国产精品99久久久久| 老汉色av国产亚洲站长工具| 九色国产91popny在线| 国产免费av片在线观看野外av| 老司机福利观看| 国产野战对白在线观看| 亚洲欧美激情综合另类| 人人妻人人澡欧美一区二区 | 在线观看免费视频日本深夜| 亚洲一码二码三码区别大吗| 欧美 亚洲 国产 日韩一| 99国产综合亚洲精品| 亚洲中文日韩欧美视频| 亚洲av电影在线进入| 欧美最黄视频在线播放免费| 国产欧美日韩精品亚洲av| 久久人人爽av亚洲精品天堂| 99精品欧美一区二区三区四区| 亚洲精品国产色婷婷电影| 在线观看免费视频网站a站| 久久精品国产99精品国产亚洲性色 | 欧美日本亚洲视频在线播放| 在线av久久热| 久久久国产欧美日韩av| 亚洲第一欧美日韩一区二区三区| 美女午夜性视频免费| 侵犯人妻中文字幕一二三四区| 俄罗斯特黄特色一大片| 国产免费男女视频| 男女之事视频高清在线观看| 欧美 亚洲 国产 日韩一| 少妇裸体淫交视频免费看高清 | 热99re8久久精品国产| 老司机靠b影院| 成人手机av| 亚洲电影在线观看av| 一边摸一边抽搐一进一小说| 日韩中文字幕欧美一区二区| 男人操女人黄网站| 18禁观看日本| 日韩一卡2卡3卡4卡2021年| 国产精品久久久久久亚洲av鲁大| 三级毛片av免费| 成人国产一区最新在线观看| 亚洲成a人片在线一区二区| 波多野结衣高清无吗| 精品一区二区三区四区五区乱码| 欧美日本亚洲视频在线播放| 91麻豆av在线| 一区福利在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产主播在线观看一区二区| 欧美国产日韩亚洲一区| 亚洲第一青青草原| 看免费av毛片| 淫秽高清视频在线观看| 无限看片的www在线观看| 久久人妻熟女aⅴ| 国产精品乱码一区二三区的特点 | 久久精品91蜜桃| 禁无遮挡网站| 国产欧美日韩一区二区三区在线| 日韩一卡2卡3卡4卡2021年| 人人妻人人澡欧美一区二区 | 一进一出好大好爽视频| 黑人欧美特级aaaaaa片| 99在线人妻在线中文字幕| 可以在线观看毛片的网站| 久久久久久国产a免费观看| 亚洲中文日韩欧美视频| 欧美日韩乱码在线| 亚洲色图av天堂| 欧美成人一区二区免费高清观看 | 女警被强在线播放| 国产三级在线视频| 成人亚洲精品av一区二区| 一区福利在线观看| 成人18禁高潮啪啪吃奶动态图| 黄片小视频在线播放| 九色国产91popny在线| 天堂动漫精品| 麻豆一二三区av精品| 国产精品爽爽va在线观看网站 | 欧美精品亚洲一区二区| 国产精华一区二区三区| 国产欧美日韩一区二区三区在线| av片东京热男人的天堂| 精品国产一区二区三区四区第35| 多毛熟女@视频| 亚洲视频免费观看视频| 啦啦啦韩国在线观看视频| 黄色视频不卡| 一本综合久久免费| 久久久久久久久久久久大奶| 动漫黄色视频在线观看| 国产极品粉嫩免费观看在线| 亚洲少妇的诱惑av| 一卡2卡三卡四卡精品乱码亚洲| 精品人妻在线不人妻| 日韩欧美国产在线观看| 亚洲精品粉嫩美女一区| 亚洲一区二区三区不卡视频| 两个人视频免费观看高清| av免费在线观看网站| 这个男人来自地球电影免费观看| 12—13女人毛片做爰片一| 精品国产一区二区久久| 久久国产精品男人的天堂亚洲| 亚洲国产精品久久男人天堂| 亚洲成av片中文字幕在线观看| 成人欧美大片| 香蕉丝袜av| 天天一区二区日本电影三级 | 纯流量卡能插随身wifi吗| 免费在线观看亚洲国产| 18禁美女被吸乳视频| 香蕉久久夜色| 国产精品久久久久久人妻精品电影| 国产精品98久久久久久宅男小说| 久久午夜亚洲精品久久| 久久欧美精品欧美久久欧美| 中文字幕色久视频| 久久精品国产清高在天天线| 啦啦啦韩国在线观看视频| 成在线人永久免费视频| 成年版毛片免费区| 真人一进一出gif抽搐免费| 欧美+亚洲+日韩+国产| 女人被狂操c到高潮| 午夜福利欧美成人| 久久人人爽av亚洲精品天堂| 一个人免费在线观看的高清视频| 色老头精品视频在线观看| 午夜福利在线观看吧| 欧美亚洲日本最大视频资源| 国产亚洲欧美98| 中文字幕人成人乱码亚洲影| 国产在线精品亚洲第一网站| 久久久精品欧美日韩精品| 亚洲国产精品合色在线| 两性夫妻黄色片| 国产精品影院久久| 国产成人啪精品午夜网站| 欧美黄色片欧美黄色片| 99久久99久久久精品蜜桃| 亚洲一码二码三码区别大吗| 久久久久久免费高清国产稀缺| 无人区码免费观看不卡| 国产精品香港三级国产av潘金莲| 黄色a级毛片大全视频| 国产精品自产拍在线观看55亚洲| 精品久久久久久久人妻蜜臀av | 久久久精品欧美日韩精品| 亚洲色图 男人天堂 中文字幕| 成年版毛片免费区| 亚洲av第一区精品v没综合| 黑人操中国人逼视频| 国产一区二区在线av高清观看| 久久人人爽av亚洲精品天堂| 1024香蕉在线观看| 欧美国产日韩亚洲一区| 一区福利在线观看| www.精华液| 一级a爱片免费观看的视频| 在线av久久热| 欧美日本亚洲视频在线播放| 国产麻豆69| 伊人久久大香线蕉亚洲五| 操出白浆在线播放| 老熟妇仑乱视频hdxx| 一区在线观看完整版| 国产一区二区三区视频了| 给我免费播放毛片高清在线观看| 亚洲人成网站在线播放欧美日韩| 中亚洲国语对白在线视频| 777久久人妻少妇嫩草av网站| 男女下面插进去视频免费观看| 不卡av一区二区三区| 国产伦人伦偷精品视频| 无限看片的www在线观看| 悠悠久久av| 久久久精品欧美日韩精品| 看片在线看免费视频| 免费高清视频大片| 夜夜躁狠狠躁天天躁| 国产精品一区二区三区四区久久 | 1024视频免费在线观看| 国产精品一区二区精品视频观看| 色播亚洲综合网| 欧美日韩精品网址| 最近最新免费中文字幕在线| 欧美精品亚洲一区二区| 久久精品91无色码中文字幕| 黑人巨大精品欧美一区二区mp4| 久久精品成人免费网站| 国产精品亚洲一级av第二区| 日韩高清综合在线| 欧美成人午夜精品| 亚洲成av人片免费观看| 国产一级毛片七仙女欲春2 | 午夜老司机福利片| 999久久久国产精品视频| 精品熟女少妇八av免费久了| 亚洲国产毛片av蜜桃av| 50天的宝宝边吃奶边哭怎么回事| 51午夜福利影视在线观看| 妹子高潮喷水视频| 老司机午夜福利在线观看视频| 国产精品免费一区二区三区在线| 精品国产国语对白av| 日韩 欧美 亚洲 中文字幕| 午夜免费鲁丝| av视频在线观看入口| 亚洲精品国产区一区二| 欧美中文日本在线观看视频| 国产精品一区二区三区四区久久 | 亚洲中文字幕日韩| 欧美黄色淫秽网站| 亚洲久久久国产精品| 亚洲专区字幕在线| 国产av又大| 亚洲精品一区av在线观看| 亚洲天堂国产精品一区在线| 大陆偷拍与自拍| 中文亚洲av片在线观看爽| 女人精品久久久久毛片| 88av欧美| 女性被躁到高潮视频| 午夜精品久久久久久毛片777| 亚洲国产看品久久| 91精品国产国语对白视频| 亚洲五月天丁香| 日韩欧美国产一区二区入口| 天堂影院成人在线观看| 丝袜美腿诱惑在线| 亚洲中文字幕一区二区三区有码在线看 | 欧美国产精品va在线观看不卡| 色播亚洲综合网| 一区二区三区激情视频| 亚洲九九香蕉| 婷婷丁香在线五月| 亚洲少妇的诱惑av| 精品欧美一区二区三区在线| 男女做爰动态图高潮gif福利片 | 黄色毛片三级朝国网站| 一个人免费在线观看的高清视频| 亚洲国产日韩欧美精品在线观看 | 9色porny在线观看| 亚洲精品美女久久久久99蜜臀| 99国产精品99久久久久| 最近最新中文字幕大全免费视频| 精品国产超薄肉色丝袜足j| 禁无遮挡网站| 黑人巨大精品欧美一区二区蜜桃| 色av中文字幕| 亚洲国产看品久久| 亚洲欧美日韩高清在线视频| 大陆偷拍与自拍| 又黄又爽又免费观看的视频| 禁无遮挡网站| 午夜精品国产一区二区电影| 国产精品秋霞免费鲁丝片| 精品久久久久久,| 欧洲精品卡2卡3卡4卡5卡区| 亚洲,欧美精品.| 好看av亚洲va欧美ⅴa在| 母亲3免费完整高清在线观看| 妹子高潮喷水视频| 欧美黄色淫秽网站| 久久中文看片网| 国产熟女xx| 免费看十八禁软件| 伊人久久大香线蕉亚洲五| 一卡2卡三卡四卡精品乱码亚洲| 亚洲熟妇中文字幕五十中出| 欧美绝顶高潮抽搐喷水| 两性夫妻黄色片| 久久久国产欧美日韩av| 最近最新中文字幕大全电影3 | 一区福利在线观看| 亚洲色图综合在线观看| 男女下面插进去视频免费观看| 欧美日韩黄片免| 日本一区二区免费在线视频| 亚洲美女黄片视频| 国产精品一区二区在线不卡| 国产精品爽爽va在线观看网站 | 亚洲国产中文字幕在线视频| 国产精品久久久久久精品电影 | 真人做人爱边吃奶动态| 黄频高清免费视频| 精品少妇一区二区三区视频日本电影|