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

    硯瓦川流域河川基流變化規(guī)律及其驅(qū)動(dòng)因素

    2021-11-26 10:21:00畢如田宋孝玉呂春娟馬耘秀李懷有
    生態(tài)學(xué)報(bào) 2021年21期
    關(guān)鍵詞:方法研究

    夏 露,畢如田,宋孝玉,呂春娟,馬耘秀,李懷有

    1 山西農(nóng)業(yè)大學(xué)資源環(huán)境學(xué)院,太谷 030801 2 西安理工大學(xué)省部共建西北旱區(qū)生態(tài)水利國(guó)家重點(diǎn)實(shí)驗(yàn)室,西安 710048 3 黃河水利委員會(huì)西峰水土保持科學(xué)試驗(yàn)站,西峰 745000

    近年來(lái),在全球氣候變化和區(qū)域人類活動(dòng)的雙重影響下,黃河流域各支流徑流量大都呈現(xiàn)減少的趨勢(shì)[1- 3]。黃河徑流量的持續(xù)減少進(jìn)一步加劇了流域水資源短缺的矛盾,直接影響到了社會(huì)經(jīng)濟(jì)的穩(wěn)定發(fā)展,因此,黃河流域徑流變化及其驅(qū)動(dòng)機(jī)制研究已成為當(dāng)前水文水資源領(lǐng)域的研究熱點(diǎn)之一?;髯鳛閬?lái)源于地下水和其他延遲部分的流量, 它是黃河流域河川徑流中相對(duì)穩(wěn)定的組成部分,是枯水期河道水流的重要補(bǔ)給來(lái)源,在維持河川流量和河流生態(tài)系統(tǒng)健康方面發(fā)揮著重要作用,同時(shí)也是流域匯流計(jì)算和水文模擬中的重要研究?jī)?nèi)容[4],但是有關(guān)變化環(huán)境下的基流變化規(guī)律研究卻還比較少見(jiàn)。

    河川徑流量可以直接測(cè)定獲得,但基流量難以直接測(cè)定,一般會(huì)采用一定的方法進(jìn)行估算,因此基流的分割一直以來(lái)就是基流研究中的重點(diǎn)和難點(diǎn)[5]。目前,國(guó)內(nèi)外學(xué)者就基流的分割已提出不少方法,大致可以分為4大類:①圖解法[6],該方法作為傳統(tǒng)的基流分割法,雖然過(guò)程簡(jiǎn)單,但在處理長(zhǎng)時(shí)間序列數(shù)據(jù)時(shí)工作量很大,且?guī)в袕?qiáng)烈的主觀性;②水文模型法[7],該方法物理意義相對(duì)明確,但在模型構(gòu)建、參數(shù)率定等過(guò)程上較為復(fù)雜繁瑣,且具有一定的不確定性;③環(huán)境同位素和水化學(xué)方法[8],該方法通過(guò)分析降雨與地下水及河流中穩(wěn)定同位素的組成來(lái)分割基流,其結(jié)果較為可靠,但是需要投入較多的人力物力,難以廣泛應(yīng)用;④數(shù)值模擬法[9- 10],該方法可以通過(guò)計(jì)算機(jī)自動(dòng)實(shí)現(xiàn)基流分割,克服了人工方法的主觀性和隨意性,操作簡(jiǎn)便,高效實(shí)用,因此目前已成為最為廣泛應(yīng)用的基流分割方法。對(duì)于同一流域,不同方法所得出的基流分割結(jié)果有所不同,有時(shí)差異可能較大,而先前較多研究只采用了一種方法進(jìn)行基流分割從而可能造成了一定的不確定性,因此有必要對(duì)不同基流分割方法的適用性進(jìn)行明確和篩選,而目前相關(guān)的方法對(duì)比研究卻比較缺乏。

    目前,有關(guān)氣候變化和人類活動(dòng)對(duì)基流的影響研究也日益增多。Zomlot等[11]通過(guò)分析比利時(shí)Flanders地區(qū)基流變化特征及其驅(qū)動(dòng)機(jī)制,發(fā)現(xiàn)植被變化是研究區(qū)基流變化的主導(dǎo)因素;Rumsey等[12]利用主成分分析法研究北美科羅拉多河上游流域基流變化的影響因素,結(jié)果表明基流量與降水量、積雪、草地及自然荒地面積比等因素呈正相關(guān),而與溫度、潛在蒸散發(fā)量、農(nóng)地及灌木林地面積比等因素呈負(fù)相關(guān)。氣候變化對(duì)基流的正負(fù)效應(yīng)一般較為明確,而人類活動(dòng)對(duì)基流變化的影響效應(yīng)可正可負(fù),且不同學(xué)者的研究結(jié)論也存在著一定的差異,例如在土地利用/覆被變化引起的基流變化研究中,Wu等[13]認(rèn)為黃河中游大面積的退耕還林還草措施增強(qiáng)了流域?qū)嶋H蒸散發(fā)能力,造成了土壤干燥化,從而最終減少了基流量,而Zhang等[14]研究發(fā)現(xiàn)渭河流域水土保持措施引起的土地利用變化是流域河川基流增加的主要原因,其內(nèi)在原因在于土地利用變化導(dǎo)致降水入滲量增加,從而增加了基流量。由此可見(jiàn),人類活動(dòng)特別是土地利用/覆被變化對(duì)基流變化的影響機(jī)制尚不明確。此外,有關(guān)基流變化的定量影響評(píng)價(jià)相對(duì)較少,還有待進(jìn)一步研究。

    因此,本文以黃河中游硯瓦川流域?yàn)檠芯繉?duì)象,基于流域1981—2016年的水文、氣象及植被資料,應(yīng)用9種基流分割方法估算流域河川基流量及基流指數(shù)并對(duì)比分析不同方法的可靠性,確定出適用于研究區(qū)的最優(yōu)基流分割方法;在此基礎(chǔ)上,揭示流域基流時(shí)間變化規(guī)律,并定性和定量的評(píng)價(jià)氣候變化和人類活動(dòng)對(duì)流域基流變化的影響,以期為黃河流域水資源開(kāi)發(fā)利用及生態(tài)環(huán)境保護(hù)提供科學(xué)依據(jù)。

    1 研究區(qū)概況

    圖1 硯瓦川流域地理位置及測(cè)站分布情況 Fig.1 Geographical location of the Yanwachuan watershed and hydrological stations in the watershed

    硯瓦川流域地處甘肅省慶陽(yáng)市西峰區(qū)及寧縣境內(nèi),位于東經(jīng)107°37′—107°55′,北緯35°31′—35°44′,系涇河支流馬蓮河右岸的一條支溝。流域總面積為385.63 km2,海拔高度為947—1432 m(圖1),主溝長(zhǎng)為35.0 km,溝壑密度為1.59 km/km2。該區(qū)屬半濕潤(rùn)季風(fēng)氣候區(qū),多年平均降水量為535.8 mm,5—9月降水量占全年降水量的78.9%,年平均氣溫為8.1℃。流域地貌主要包括塬面、梁峁坡和溝谷這三種類型,分別占流域總面積的53.7%、17.7%和28.6%,具有典型的黃土高塬溝壑區(qū)地貌特征。硯瓦川流域地質(zhì)構(gòu)造比較單一,地表主要被第四紀(jì)黃土所覆蓋,厚度可達(dá)200 m左右,土壤侵蝕非常嚴(yán)重。為了有效地控制水土流失,黃委會(huì)西峰水保站于1975年選定硯瓦川流域作為黃土高塬溝壑區(qū)的典型中尺度流域,進(jìn)行流域綜合治理試驗(yàn)研究。

    流域水土保持治理總體上可分為以下3個(gè)階段[15]:

    第一階段為1975—1980年:1975年當(dāng)?shù)卣闪⒘顺幫叽饔蛑卫碇笓]部,開(kāi)始對(duì)硯瓦川流域集中開(kāi)展治理工作,到1979年,已修成梯田3886.4 hm2,人工造林1768.6 hm2,人工種草988.7 hm2。

    第二階段為1981—1993年:進(jìn)入80年代,隨著農(nóng)村實(shí)行包產(chǎn)到戶,土地使用權(quán)歸農(nóng)民所有,這一時(shí)期國(guó)家投資力度處于下降趨勢(shì),流域治理面積變化不大,有的地方治理面積有所減少。

    第三階段為1994—2016年:1994年,黃土高原世界銀行貸款馬蓮河流域治理項(xiàng)目開(kāi)始實(shí)施,這是當(dāng)?shù)卣状我M(jìn)外資進(jìn)行流域綜合治理,使得流域治理步伐開(kāi)始加快。至此之后,涉及本研究流域的的綜合治理項(xiàng)目主要還有:1999年開(kāi)始實(shí)施的退耕還林工程、2001—2005年實(shí)施的黃河水土保持生態(tài)工程齊家川示范區(qū)項(xiàng)目以及2008—2010年實(shí)施的黃河水土保持生態(tài)工程硯瓦川示范區(qū)項(xiàng)目。

    2 數(shù)據(jù)與方法

    2.1 數(shù)據(jù)來(lái)源與處理

    本研究所用觀測(cè)數(shù)據(jù)主要包括流域水文、氣象和植被資料。

    (1)水文資料:主要包括降水資料和徑流泥沙資料,均來(lái)源于黃河水利委員會(huì)西峰水土保持科學(xué)試驗(yàn)站,數(shù)據(jù)年限為1981—2016年。其中,降水資料選用了流域內(nèi)系列連續(xù)且分布均勻的10個(gè)雨量站(圖1)共36年的降水?dāng)?shù)據(jù),流域面降水量采用泰森多邊形法求得;徑流泥沙資料采用流域出口控制站—硯瓦川水文站的多年觀測(cè)數(shù)據(jù),該水文站控制流域面積為341.54 km2。

    (2)氣象資料:采用緊鄰硯瓦川流域邊界的西峰國(guó)家氣象觀測(cè)站1981—2016年的逐日數(shù)據(jù)資料,包括氣溫(最高氣溫、最低氣溫、平均氣溫)、相對(duì)濕度、太陽(yáng)輻射、日照時(shí)數(shù)和風(fēng)速等?;谠撜径嗄曛鹑諝庀髷?shù)據(jù),采用世界糧農(nóng)組織推薦的Penman-Monteith公式[16]計(jì)算流域多年潛在蒸散發(fā)量。

    (3)植被資料:本文用NDVI數(shù)據(jù)表征植被變化對(duì)基流的影響,采用的是AVHRR GIMMS NDVI(數(shù)據(jù)來(lái)源http://www.geodata.cn;數(shù)據(jù)年限為1981—2006年)和SPOT/VEGETATION NDVI(數(shù)據(jù)來(lái)源http://www.resdc.cn;數(shù)據(jù)年限為1998年至今)數(shù)據(jù)集。由于這兩套數(shù)據(jù)均不能完全覆蓋本研究期限,需要進(jìn)行插補(bǔ)獲得全時(shí)段NDVI,具體操作如下:首先對(duì)兩套數(shù)據(jù)重合時(shí)段(1998—2006年)的月NDVI值進(jìn)行了相關(guān)分析,結(jié)果表明兩者呈顯著性相關(guān),相關(guān)系數(shù)為0.96,因此可以利用這兩套數(shù)據(jù)進(jìn)行插補(bǔ)延長(zhǎng),然后建立起兩套數(shù)據(jù)1998—2006年的月NDVI值線性回歸方程,最終將GIMMS 1981—1997的NDVI數(shù)據(jù)代入回歸方程,從而得到插補(bǔ)延長(zhǎng)后全時(shí)段的SPOT/VEGETATION NDVI數(shù)據(jù)。

    2.2 研究方法

    2.2.1基流的分割方法

    本研究采用數(shù)值模擬法中的數(shù)字濾波法、平滑最小值法、時(shí)間步長(zhǎng)法共3類9種方法對(duì)流域基流進(jìn)行分割,并利用實(shí)際基流量對(duì)各種基流分割方法進(jìn)行參數(shù)率定。

    (1)數(shù)字濾波法

    數(shù)字濾波法源于信號(hào)分析和處理技術(shù),主要功能是將信號(hào)分為高頻信號(hào)和低頻信號(hào)。由于直接徑流的快速響應(yīng)特征與高頻信號(hào)類似,而基流的慢速響應(yīng)特征與低頻信號(hào)類似,因此該方法的原理將日徑流看作是直接徑流和基流的疊加,通過(guò)數(shù)字濾波器分解出高頻和低頻信號(hào),從而相應(yīng)的將直接徑流和基流從徑流中分割出來(lái)[17]。目前應(yīng)用最為廣泛的數(shù)字濾波法主要包括以下4種:

    ①Lyne-Hollick濾波法

    Lyne-Hollick濾波法是由Lyne和Hollick在1979年首次提出,并由Nathan和Mcmahon于1990年首次引入到水文中進(jìn)行基流分割,其濾波方程為[18]:

    Qd(i)=f1Qd(i-1)+f2(1+f1)[Q(i)-Q(i-1)]Qb(i)=Q(i)-Qd(i)

    (1)

    式中:Qd(i)和Qd(i-1)分別為第i和第i-1時(shí)刻的直接徑流(m3/s);Q(i)和Q(i-1)分別為第i和第i-1時(shí)刻的河川流量(m3/s);Qb(i)為第i時(shí)刻的基流(m3/s);f1和f2為濾波參數(shù),推薦取值分別為0.925和0.5,另外濾波次數(shù)(N)一般取1或3。本研究通過(guò)對(duì)f1(調(diào)參范圍在0.90—0.98之間)、f2(調(diào)參范圍在0.4—0.6之間)和N(調(diào)參范圍為1、3和5)進(jìn)行率定,發(fā)現(xiàn)當(dāng)f1為0.975、f2為0.5、N為1時(shí),基流分割結(jié)果精度較高。

    ②Chapman濾波法

    Chapman通過(guò)對(duì)Lyne-Hollick濾波方程中的問(wèn)題進(jìn)行系統(tǒng)的分析,在1991年對(duì)該濾波法進(jìn)行了改進(jìn),提出了Chapman濾波法,其濾波方程為[19]:

    (2)

    式中:f1為濾波參數(shù),推薦取值為0. 95,另外濾波次數(shù)(N)一般取1或3,其他參數(shù)含義同上,再用式(1)中的基流公式計(jì)算基流。本研究通過(guò)對(duì)f1(調(diào)參范圍在0.90—0.98之間)和N(調(diào)參范圍為1、3和5)進(jìn)行率定,發(fā)現(xiàn)當(dāng)f1為0.925、N為1時(shí),基流分割結(jié)果較優(yōu)。

    ③Chapman-Maxwell濾波法

    Chapman和Maxwell在1991年對(duì)以上濾波法進(jìn)行進(jìn)一步改進(jìn),他們假定某時(shí)刻的基流為該時(shí)刻的直接徑流和前一時(shí)刻基流的加權(quán)平均,其提出的濾波方程為[20]:

    (3)

    式中:f1為退水參數(shù),推薦取值為0. 95,另外濾波次數(shù)(N)一般取1或3,其他參數(shù)含義同上。本研究通過(guò)對(duì)f1(調(diào)參范圍在0.90—0.98之間)和N(調(diào)參范圍為1、3和5)進(jìn)行率定,發(fā)現(xiàn)當(dāng)f1為0.925、N為1時(shí),基流分割結(jié)果較優(yōu)。

    ④Boughton-Chapman濾波法

    Boughton于1993年提出了Boughton-Chapman濾波法,其濾波方程為[21]:

    (4)

    式中:f1和f2為方程參數(shù),推薦取值分別為0. 95和0.15,另外濾波次數(shù)(N)一般取1或3,其他參數(shù)含義同上。本研究通過(guò)對(duì)f1(調(diào)參范圍在0.90—0.98之間)、f2(調(diào)參范圍在0.05—0.35之間)和N(調(diào)參范圍為1、3和5)進(jìn)行率定,發(fā)現(xiàn)當(dāng)f1為0.95、f2為0.1、N為1時(shí),基流分割結(jié)果精度較高。

    (2)平滑最小值法

    平滑最小值法是英國(guó)水文研究所在1980年提出的一種基流計(jì)算方法,主要有標(biāo)準(zhǔn)平滑最小值(f)法和改進(jìn)平滑最小值(k)法,其方法原理為[22]:將連續(xù)的日徑流系列按照時(shí)間間隔N劃分成365/N個(gè)時(shí)間段,然后確定每個(gè)時(shí)間段的最小流量值,形成最小流量序列Q1、Q2、…、QN,若某個(gè)i時(shí)段內(nèi)的流量最小值與拐點(diǎn)檢驗(yàn)因子f或k的乘積小于等于左右相鄰兩個(gè)時(shí)段內(nèi)的流量最小值,即fQi≤min(Qi-1,Qi+1)[kQi≤min(Qi-1,Qi+1)],則將中間點(diǎn)確定為拐點(diǎn)。重復(fù)此過(guò)程,在流量過(guò)程線上確定出所有拐點(diǎn),將所有拐點(diǎn)用直線連接即可得到基流過(guò)程線。

    該方法需要確定2個(gè)參數(shù):時(shí)間間隔N和拐點(diǎn)檢驗(yàn)因子f或k。N一般取值為3—5 d,而f和k推薦值為0.9和0.979。本研究通過(guò)對(duì)f(調(diào)參范圍在0.85—0.95之間)、k(調(diào)參范圍在0.90—0.99之間)和N(調(diào)參范圍為3、4和5)進(jìn)行率定,發(fā)現(xiàn)當(dāng)f為0.9或k為0.98、N為3時(shí),基流分割結(jié)果較優(yōu)。

    (3)時(shí)間步長(zhǎng)法

    時(shí)間步長(zhǎng)法也稱為HYSEP法,是由Pettyjohn和Henning首先于1979年提出[23],并由美國(guó)地質(zhì)調(diào)查局開(kāi)發(fā)計(jì)算程序推薦使用。HYSEP法共有固定步長(zhǎng)法(Fixed interval,FI)、滑動(dòng)步長(zhǎng)法(Sliding interval,SI)和局部最小值法(Local minimum,LM)3 種不同的分割方法,這3種方法都首先需要利用經(jīng)驗(yàn)公式計(jì)算直接徑流的持續(xù)時(shí)間t:

    t=(2.95A)0.2

    (5)

    式中:t為直接徑流的持續(xù)時(shí)間(d);A為流域面積(km2)。

    通過(guò)計(jì)算直接徑流的持續(xù)時(shí)間,選擇與2t最為接近且介于3到11之間的奇數(shù)作為時(shí)間間隔N。

    硯瓦川水文站控制流域面積為341.54 km2,計(jì)算得到t為3.99,故時(shí)間間隔N為7。之后,這3種方法計(jì)算基流的原理有所不同,其中,固定步長(zhǎng)法(FI)是將該時(shí)間間隔內(nèi)的最小流量作為該時(shí)段內(nèi)任意一天的基流;滑動(dòng)步長(zhǎng)法(SI)是將某天前后(2N-1)/2d內(nèi)的最小流量作為該天的基流,由此計(jì)算出每一天的基流;局部最小值法(LM)是選擇時(shí)間步長(zhǎng)內(nèi)中心點(diǎn)前后(2N-1)/2d內(nèi)的最小流量作為相鄰時(shí)間步長(zhǎng)內(nèi)中心點(diǎn)的基流值,然后通過(guò)線性內(nèi)插得到步長(zhǎng)中心點(diǎn)之外時(shí)段的基流。3種方法均以本次計(jì)算的終點(diǎn)作為下次時(shí)間的起點(diǎn),重復(fù)以上過(guò)程便可得到基流過(guò)程線。

    2.2.2基流分割結(jié)果的評(píng)價(jià)

    通過(guò)確定硯瓦川流域歷年的實(shí)際基流量,利用統(tǒng)計(jì)指標(biāo)評(píng)價(jià)不同基流分割方法的精度,從而選取最優(yōu)基流計(jì)算方法的結(jié)果開(kāi)展進(jìn)一步的基流分析。

    (1)實(shí)際基流量的確定

    直接觀測(cè)河川基流量是一個(gè)非常困難的過(guò)程,因此目前國(guó)際上主要采用枯水指數(shù)法[24]來(lái)近似估計(jì)實(shí)際基流量,該方法的計(jì)算步驟為:利用歷年的逐日流量數(shù)據(jù)繪制出歷年的日流量歷時(shí)曲線,然后推求歷年的Q90和Q50值,這兩個(gè)值分別代表某一年的低流量和中流量,那么本研究用枯水指數(shù)(Q90/Q50)與年均流量相乘,即推求出年基流量的實(shí)際值。

    (2)基流分割方法的評(píng)價(jià)

    將實(shí)際年基流量值與基流分割估算結(jié)果進(jìn)行對(duì)比,利用納什效率系數(shù)(NSE)和平均絕對(duì)相對(duì)誤差(Re)來(lái)率定不同基流分割方法的參數(shù),并評(píng)價(jià)其優(yōu)劣性。其中,NSE的計(jì)算公式如下[25]:

    (6)

    在施工過(guò)程中錘頭脫落往往不是單一原因造成的,其主要因素是鋼絲繩的磨耗斷裂,卡扣失靈掉落,錘頭橫向軸承過(guò)勞斷裂等。

    Re可由下式計(jì)算得出[26]:

    (7)

    式中,Re值越小,估算效果就越好。一般認(rèn)為,當(dāng)Re值小于20%時(shí),估算結(jié)果就具有一定的精度。

    2.2.3年基流變化趨勢(shì)及突變分析

    本研究采用Mann-Kendall檢驗(yàn)法進(jìn)行年基流的變化趨勢(shì)及突變分析,并采用滑動(dòng)t檢驗(yàn)法對(duì)Mann-Kendall突變檢驗(yàn)的結(jié)果進(jìn)行驗(yàn)證。由于這兩種方法在水文研究中已被廣泛使用,因此有關(guān)方法原理在此并不詳述,具體方法詳見(jiàn)文獻(xiàn)[27-28]。

    2.2.4基流對(duì)氣候變化和人類活動(dòng)的響應(yīng)分析

    本文通過(guò)參考以往的基流驅(qū)動(dòng)因素研究,最終選取降水量和潛在蒸散發(fā)量這兩個(gè)指標(biāo)反映氣候變化的影響,同時(shí)考慮到研究區(qū)水土保持治理以植被措施為主,大型水利工程較少,流域內(nèi)人類取用水量也非常少,人類活動(dòng)主要通過(guò)改變土地利用/覆蓋導(dǎo)致基流的變化,因此,選取NDVI這一指標(biāo)反映人類活動(dòng)(植被變化)對(duì)基流的影響。

    選取出主要影響因素后,對(duì)流域年降水量、年潛在蒸散發(fā)量和年NDVI值進(jìn)行變化趨勢(shì)及其與基流的相關(guān)性分析,探討基流對(duì)氣候變化和植被變化的響應(yīng)規(guī)律,并定量評(píng)價(jià)氣候變化和植被變化對(duì)基流變化的影響程度。本研究采用多元線性回歸法來(lái)計(jì)算基流變化的貢獻(xiàn)率,其計(jì)算步驟如下[24]:以標(biāo)準(zhǔn)化年基流量為因變量,以主要影響因子(即標(biāo)準(zhǔn)化年降水量、年潛在蒸散發(fā)量和年NDVI)為自變量,建立標(biāo)準(zhǔn)化年基流量的多元線性回歸方程,其中,每個(gè)變量的標(biāo)準(zhǔn)化值由各變量系列值減去變量均值再除以該變量的標(biāo)準(zhǔn)差求得。那么,每個(gè)影響因子對(duì)基流量變化的貢獻(xiàn)率大小就可以由各自的回歸系數(shù)來(lái)估計(jì),計(jì)算公式如下:

    (8)

    3 結(jié)果與分析

    3.1 基流分割結(jié)果的對(duì)比與評(píng)價(jià)

    表1列出了9種不同基流分割方法的計(jì)算結(jié)果及評(píng)價(jià)指標(biāo),總體來(lái)看,除了Chapman濾波法和Chapman-Maxwell濾波法計(jì)算出來(lái)的BFI值低于0.45之外,其他7種方法的BFI計(jì)算結(jié)果都處于0.55—0.59之間,與枯水指數(shù)法計(jì)算出來(lái)的多年實(shí)際BFI平均值0.58較為接近,且都明顯高于Chapman濾波法和Chapman-Maxwell濾波法的計(jì)算結(jié)果,這可能是由于這兩種方法的計(jì)算原理有所不同,導(dǎo)致基流分割過(guò)程在低流量時(shí)更為平滑,基流量及基流指數(shù)低于其他方法的結(jié)果。不同基流分割方法的評(píng)價(jià)指標(biāo)結(jié)果顯示,Chapman濾波法、Chapman-Maxwell濾波法、平滑最小值法和時(shí)間步長(zhǎng)法的NSE都低于0.4,說(shuō)明分割效果不夠理想,而Lyne-Hollick濾波法的NSE和Re分別為0.49和17%,Boughton-Chapman濾波法的NSE和Re分別為0.64和16%,說(shuō)明這兩種方法具有較高的精度,其中最優(yōu)的還屬Boughton-Chapman濾波法。

    表1 不同基流分割方法的基流指數(shù)值及評(píng)價(jià)結(jié)果

    圖2 不同基流分割方法的1981年日流量和日基流量過(guò)程線Fig.2 Baseflow hydrographs estimated with the 9 baseflow separating methods

    為了進(jìn)一步驗(yàn)證不同基流分割方法的適用性,本研究以1981年(頻率為10%的豐水年)為例,對(duì)比分析不同日基流分割過(guò)程線的特征,結(jié)果見(jiàn)圖2(為了使日流量過(guò)程線展示清晰,已刪除3個(gè)日流量大于4 m3/s的數(shù)據(jù)點(diǎn))。在4種數(shù)字濾波法中,Boughton-Chapman濾波法日基流過(guò)程線起伏最大,不太符合實(shí)際的日基流漲落過(guò)程,Chapman濾波法和Chapman-Maxwell法波動(dòng)幅度次之,而Lyne-Hollick法的日基流過(guò)程線較為平滑,同時(shí)對(duì)汛期降水有一定的響應(yīng),在非汛期BFI值較大而在汛期BFI值小于其他3種濾波法,能夠較好地反映日基流的實(shí)際變化規(guī)律。兩種平滑最小值法的日基流過(guò)程相對(duì)于其他方法波動(dòng)幅度最小,但是過(guò)于平滑,沒(méi)能體現(xiàn)日基流對(duì)汛期降水的響應(yīng)過(guò)程。在3種時(shí)間步長(zhǎng)法中,固定步長(zhǎng)法和滑動(dòng)步長(zhǎng)法的日基流變化過(guò)程與日徑流變化過(guò)程非常一致,未能體現(xiàn)基流的滯后性這一特征,這是由于這2種方法的分割原理所決定的,而局部最小值法整體波動(dòng)非常小,日基流過(guò)程線過(guò)于平滑??傮w來(lái)看,Lyne-Hollick法的基流分割結(jié)果比較符合實(shí)際日基流變化規(guī)律。

    綜合年基流指數(shù)、年基流量誤差分析與日基流過(guò)程線結(jié)果,本研究認(rèn)為L(zhǎng)yne-Hollick法在所有基流分割方法中具有更好的合理性和可靠性??紤]到Lyne-Hollick法在年尺度上的基流估算還存在著一定的誤差,所以在后續(xù)的年際尺度基流研究中直接利用枯水指數(shù)法計(jì)算的實(shí)際基流量成果,而在年內(nèi)尺度基流研究中則采用Lyne-Hollick法的分割結(jié)果。

    3.2 基流年內(nèi)變化規(guī)律

    圖3 硯瓦川流域多年平均月徑流和月基流變化規(guī)律 Fig.3 Variation of annual average monthly streamflow and monthly baseflow in the Yanwachuan watershedBFI:基流指數(shù) Baseflow index

    根據(jù)Lyne-Hollick法的逐日基流分割結(jié)果,得到硯瓦川流域1981—2016年平均月基流量及月BFI變化過(guò)程(圖3)。硯瓦川流域1—3月份河川流量逐漸增加,同時(shí)月基流量也呈現(xiàn)上升的趨勢(shì),并在3月份達(dá)到第一個(gè)峰值,這和積雪融化補(bǔ)充徑流及基流關(guān)系密切;4—6月份河川流量與3月份相比有所下降,這一方面是由于這個(gè)時(shí)期降水量較小,另一方面是因?yàn)橹参锏纳L(zhǎng)使降水入滲量增加,從而減少了徑流,同時(shí)期月基流量也呈現(xiàn)下降的趨勢(shì);7月份河川流量達(dá)到全年最高峰值,這主要與7月份的降水量達(dá)到峰值有關(guān),但是該月份的基流量卻未達(dá)到最大,另外7—10月是徑流主要的集中期,河川流量隨著降水量的減少而逐漸降低,但是基流量在7—10月期間卻穩(wěn)步增加,在10月份達(dá)到了最大值,這體現(xiàn)了基流作為慢速流的滯后性;11—12月份河川流量徑流量和基流量均呈下降趨勢(shì)??傮w來(lái)看,基流在年內(nèi)的變化幅度小于徑流,且對(duì)降水/融雪補(bǔ)給的響應(yīng)時(shí)間也較為緩慢,這都與基流的補(bǔ)給機(jī)制密切相關(guān)。另外,BFI在年內(nèi)總體呈現(xiàn)先減少后增加的變化規(guī)律,汛期(5—9月份)的月BFI值明顯低于非汛期,非汛期各月BFI值達(dá)到0.56—0.81之高,這也驗(yàn)證了在枯水期基流量是河川徑流的主要組成部分這一結(jié)論。

    3.3 基流年際變化規(guī)律

    根據(jù)年徑流及基流統(tǒng)計(jì)分析結(jié)果(圖4),流域多年平均河川流量和基流量分別為0.264 m3/s和0.152 m3/s,其中年流量和年基流量最大值分別為0.564 m3/s和0.322 m3/s,均出現(xiàn)在1984年,而年流量和年基流量最小值分別為0.180 m3/s和0.129 m3/s,分別發(fā)生在1999年和2007年。在整個(gè)研究期間,流域年基流量與年流量均呈現(xiàn)出下降的趨勢(shì),且兩者的變化過(guò)程較為一致,基流在年際尺度上并未體現(xiàn)出一定的滯后性。流域多年平均BFI值為0.58,年BFI值在0.38—0.76之間浮動(dòng),整體表現(xiàn)出下降的趨勢(shì),特別是在2006年以后,流域BFI值維持在一個(gè)較低的水平,這說(shuō)明從2006年之后降水更多的形成了快速?gòu)搅?使慢速流在徑流中的占比較低。

    圖4 硯瓦川流域1981—2016年年徑流和年基流變化規(guī)律Fig.4 Variation of annual streamflow and annual baseflow in the Yanwachuan watershed from 1981 to 2016

    年基流量和年BFI值的M-K檢驗(yàn)和滑動(dòng)t檢驗(yàn)結(jié)果如圖5所示。年基流量和年BFI值的UF值分別為-3.787和-2.642,超過(guò)了P=0.01顯著性水平,這說(shuō)明硯瓦川流域年河川基流量和年BFI值均呈現(xiàn)極顯著的減少趨勢(shì)。M-K突變檢驗(yàn)結(jié)果顯示,年基流量UF和UB線相交于1993年,年基流指數(shù)UF和UB線相交于2006年,而t檢驗(yàn)顯示年基流量于1992年附近發(fā)生了突變,年基流指數(shù)在2006年和2007年附近發(fā)生了突變,兩種突變檢驗(yàn)的結(jié)果非常一致。因此,可以認(rèn)為硯瓦川流域年河川基流量和年基流指數(shù)在研究期間發(fā)生了極顯著的減少趨勢(shì)性變化,且分別于1993年和2006年附近發(fā)生了突變。

    圖5 硯瓦川流域年河川基流量和年基流指數(shù)M-K檢驗(yàn)和滑動(dòng)t檢驗(yàn)結(jié)果Fig.5 Results of M-K test and t test for annual baseflow and annual baseflow index in the Yanwachuan watershedUF: Mann-Kendall統(tǒng)計(jì)值組成的曲線; UB: Mann-Kendall反序列統(tǒng)計(jì)值組成的曲線; t: 滑動(dòng)t檢驗(yàn)的統(tǒng)計(jì)值

    3.4 基流變化驅(qū)動(dòng)因素分析

    通過(guò)對(duì)硯瓦川流域1981—2016年降水量(P)、潛在蒸散發(fā)量(ET0)和NDVI數(shù)據(jù)進(jìn)行趨勢(shì)分析及突變檢驗(yàn)(圖6和表2),結(jié)果表明:流域多年平均年降水量和年潛在蒸散發(fā)量分別為535.8 mm和1172.1 mm,在整個(gè)研究期內(nèi)年降水量呈現(xiàn)不顯著的增加趨勢(shì)(P>0.05),以平均每年1.351 mm的速度在增加,且未發(fā)生突變,而潛在蒸散發(fā)在年際上呈現(xiàn)極顯著的增加趨勢(shì)(P<0.01),以平均每年4.121 mm的速度在增加,且于1993年附近發(fā)生了突變;流域多年平均年NDVI值為0.535,年NDVI呈現(xiàn)極顯著的增加趨勢(shì)(P<0.01),以平均每年0.0056的速度在增加,且于2004年附近發(fā)生了突變。綜合以上因素的趨勢(shì)分析及突變檢驗(yàn)結(jié)果,可以發(fā)現(xiàn)年潛在蒸散發(fā)的突變時(shí)間與年基流量的突變時(shí)間均出現(xiàn)在1993年附近,且自1994年開(kāi)始的黃土高原世界銀行貸款使得流域治理步伐開(kāi)始加快,流域植被條件也因此開(kāi)始發(fā)生了較大的變化,這說(shuō)明硯瓦川流域潛在蒸散發(fā)和NDVI的增加可能是河川基流量減小的主要原因;另外,年NDVI的突變時(shí)間與年基流指數(shù)的突變時(shí)間較為接近,這表明NDVI的變化可能對(duì)于流域基流指數(shù)的減小趨勢(shì)起到了一定的作用,還有待進(jìn)一步驗(yàn)證。

    圖6 硯瓦川流域1981—2016年降水量、潛在蒸散發(fā)量和NDVI變化規(guī)律Fig.6 Variation of precipitation, potential evapotranspiration and NDVI in the Yanwachuan watershed during 1981—2016NDVI: 歸一化植被指數(shù)Normalized difference vegetation index

    表2 硯瓦川流域年降水量、潛在蒸散發(fā)量和NDVI趨勢(shì)分析及突變檢驗(yàn)結(jié)果

    為了進(jìn)一步揭示這三者對(duì)基流的影響規(guī)律,本研究利用Spearman秩相關(guān)檢驗(yàn)對(duì)三者與基流量和基流指數(shù)進(jìn)行相關(guān)性分析,結(jié)果見(jiàn)表3。根據(jù)Spearman相關(guān)分析結(jié)果,基流量與三者均顯著相關(guān),其中,潛在蒸散發(fā)量相關(guān)性最大(R= -0.510),這說(shuō)明潛在蒸散發(fā)應(yīng)該是影響基流量最重要的因素,這也從基流量的突變年份與潛在蒸散發(fā)的突變年份一致這一點(diǎn)上得到了驗(yàn)證;基流量與降水量和NDVI均顯著相關(guān)(P<0.05),不過(guò)其與降水量呈正相關(guān),而與NDVI卻呈負(fù)相關(guān),這表明了流域降水量的上升趨勢(shì)確實(shí)會(huì)使基流量有所增加,而NDVI的上升趨勢(shì)導(dǎo)致基流量有所減少。另外,在3個(gè)影響因素中,NDVI與基流指數(shù)的相關(guān)性最高,降水量與基流指數(shù)的相關(guān)性次之,但是呈負(fù)相關(guān)關(guān)系,而潛在蒸散發(fā)與基流指數(shù)的相關(guān)性很弱,相關(guān)系數(shù)僅為0.079。

    表3 硯瓦川流域基流與主要影響因素間的Spearman相關(guān)分析結(jié)果

    為了更深入的分析驅(qū)動(dòng)因素對(duì)基流變化的定量影響,采用多元線性回歸法評(píng)價(jià)基流變化的貢獻(xiàn)率。表4列出了流域基流量和基流指數(shù)的線性回歸系數(shù)及貢獻(xiàn)率結(jié)果,兩個(gè)線性回歸方程均通過(guò)了0.05顯著性水平上的F檢驗(yàn),但基流量線性回歸方程的擬合度(R2=0.446)要優(yōu)于基流指數(shù)線性回歸方程(R2=0.241)。根據(jù)基流變化貢獻(xiàn)率評(píng)價(jià)結(jié)果,流域降水量變化、潛在蒸散發(fā)量變化及NDVI變化對(duì)基流量變化的貢獻(xiàn)率分別為-99.1%、113.3%和85.8%,即降水量的增加使流域基流量增加了99.1%,而潛在蒸散發(fā)的增加以及NDVI的變化分別使基流量減少了113.3%和85.8%,可見(jiàn)潛在蒸散發(fā)和NDVI的變化是引起河川基流量減小的主要原因,這與前面得出的結(jié)論一致;流域降水量變化、潛在蒸散發(fā)量變化及NDVI變化對(duì)基流指數(shù)變化的貢獻(xiàn)率分別為41.3%、-27.7%和86.5%,這說(shuō)明流域NDVI的增加確實(shí)對(duì)于基流指數(shù)的降低起到了決定性的作用。

    表4 硯瓦川流域基流與基流指數(shù)的貢獻(xiàn)率評(píng)價(jià)結(jié)果

    4 討論

    不同基流分割方法的計(jì)算結(jié)果有所差異,在本文采用的9種數(shù)值模擬方法中,Chapman濾波法和Chapman-Maxwell濾波法計(jì)算出來(lái)的基流量結(jié)果明顯低于其他7種方法,且與枯水指數(shù)法計(jì)算出來(lái)的實(shí)際BFI值相差較遠(yuǎn),這與吳珍妮等[29]在北洛河上游基流分割方法適用性研究中得出的結(jié)論一致,這表明Chapman濾波法和Chapman-Maxwell濾波法可能不適用于黃河中游地區(qū)的基流分割。另外,本研究通過(guò)綜合分析不同方法的年基流誤差與日基流過(guò)程線結(jié)果,發(fā)現(xiàn)Lyne-Hollick法在所有基流分割方法中具有更好的可靠性。吳珍妮等[29]、申戀綿等[30]、亢小語(yǔ)等[31]和周嘉欣等[5]曾分別對(duì)北洛河上游、窟野河流域、昕水河流域和疏勒河上游不同基流分割方法的適用性進(jìn)行了評(píng)價(jià),結(jié)果均發(fā)現(xiàn)Lyne-Hollick法基流估算結(jié)果穩(wěn)定可靠且精度高,這說(shuō)明Lyne-Hollick濾波法確實(shí)是黃河中上游比較適宜的基流分割方法。

    硯瓦川流域河川基流量和基流指數(shù)在年際尺度上呈極顯著下降趨勢(shì)(P<0.01),氣候變化和植被變化是影響研究區(qū)河川基流的主要驅(qū)動(dòng)因素。在氣候要素中,降水量和潛在蒸散發(fā)量是最關(guān)鍵的兩個(gè)因子。降水是基流的補(bǔ)給來(lái)源,它直接決定著基流量的大小,因此在一般情況下降水量與基流量之間呈正相關(guān)關(guān)系,這一點(diǎn)也被眾多研究所證實(shí)。本研究還分析了年降水量和年基流指數(shù)之間的相關(guān)關(guān)系,發(fā)現(xiàn)兩者呈顯著負(fù)相關(guān),即降水量的增加會(huì)顯著減小基流指數(shù)。Xu等[32]曾采用彈性系數(shù)法分析降水和潛在蒸散發(fā)對(duì)基流的定量影響,他們也發(fā)現(xiàn)了基流的降水彈性系數(shù)為負(fù)值這一現(xiàn)象,但并未作相應(yīng)的解釋,本研究認(rèn)為這可能是由于流域降水量一般集中在汛期,年降水量的增加意味著汛期降水量的增加,而汛期降水量更容易產(chǎn)生快速流,從而減小了基流占徑流的比例即減小了基流指數(shù)。

    潛在蒸散發(fā)反映了大氣的蒸發(fā)能力,潛在蒸散發(fā)能力越強(qiáng),降水入滲量會(huì)更多的轉(zhuǎn)化為蒸散發(fā)量,基流量也會(huì)相應(yīng)地減少,因此兩者呈負(fù)相關(guān)關(guān)系(表3)。在氣象要素中,基流量與潛在蒸散發(fā)量相關(guān)性最大,而與降水量次之,這與傳統(tǒng)的氣候變化對(duì)徑流的影響結(jié)論不太一致,例如楊大文等[33]基于流域水熱耦合平衡模型對(duì)黃河流域38個(gè)典型流域1961—2010年間徑流變化進(jìn)行分析,結(jié)果表明在所有流域中徑流變化對(duì)降水變化比對(duì)潛在蒸散發(fā)變化更為敏感。這種不一致的原因可能在于潛在蒸散發(fā)能更大程度上決定降水入滲量的消耗途徑:它到底是會(huì)更多的形成蒸散發(fā)量還是會(huì)補(bǔ)給基流?這比降水量對(duì)基流的影響更加直接和密切。

    有關(guān)植被變化對(duì)徑流的影響研究非常多見(jiàn),但是針對(duì)植被變化對(duì)基流的影響研究卻相對(duì)匱乏,且其定性影響結(jié)論尚有一定的差異。例如黃河中游大規(guī)模的生態(tài)恢復(fù)特別是人工造林對(duì)基流的影響仍存在爭(zhēng)議:一些學(xué)者[4,8,14,34]認(rèn)為造林或植被恢復(fù)能增加土壤入滲,從而增加河川基流量;另一些學(xué)者[7,13,35]認(rèn)為造林或植被恢復(fù)增加了植被蒸騰和截留蒸發(fā),使得基流量有所減少。造成這種不一致的結(jié)果可能在于影響基流的因素錯(cuò)綜復(fù)雜,一般的研究很難將植被變化的影響從眾多因素中剝離出來(lái),容易受到其他因素的干擾,從而得到了不同的研究結(jié)論,因此在以后的研究中亟需利用對(duì)比流域法來(lái)直接評(píng)判這一效應(yīng)。本文的結(jié)論更傾向于后者,即在黃河中游地區(qū)植被覆蓋度的增加會(huì)減小基流量,NDVI與基流量的負(fù)相關(guān)關(guān)系主要是因?yàn)辄S河中游干旱半干旱地區(qū)作為水分限制區(qū),植被恢復(fù)增加的降水入滲不能抵消植被恢復(fù)引起的蒸散發(fā)(植被蒸騰和截留蒸發(fā))增加,最終導(dǎo)致可以補(bǔ)給基流的降水入滲量有所減少。硯瓦川流域人類活動(dòng)對(duì)基流的影響主要體現(xiàn)在植被變化上,因此該結(jié)果也有一定的參考性。另外,在3個(gè)影響因素中,基流指數(shù)與NDVI的相關(guān)性最高,這可能是因?yàn)镹DVI能更大程度上決定徑流的組成成分,穆興民等[36]指出黃土高原植被恢復(fù)深刻地改變了黃土高原土壤水文性質(zhì),使土壤入滲能力增強(qiáng)、流域蓄水容量增大,雖然這意味著直接徑流量的減少,但是由于流域NDVI的極顯著增加導(dǎo)致蒸散發(fā)能力的增強(qiáng),使補(bǔ)給基流的水量大大減小,最終導(dǎo)致基流的減少速度要快于直接徑流的減少速度。

    硯瓦川流域基流變化的定量評(píng)價(jià)結(jié)果表明,潛在蒸散發(fā)和NDVI的變化是引起流域河川基流量減少的主要原因,而流域NDVI的增加對(duì)于硯瓦川流域基流指數(shù)的降低起到了決定性的作用。Wu等[13]對(duì)黃河中游11個(gè)典型流域的基流變化進(jìn)行歸因分析,發(fā)現(xiàn)人類活動(dòng)是造成大部分流域河川基流量減少的原因,而基流指數(shù)的歸因結(jié)果在不同流域分異較大。不過(guò)這些流域的人類活動(dòng)不僅僅局限于植被變化,而如對(duì)地下水的開(kāi)采、煤炭開(kāi)采對(duì)地下含水層的破壞等等均是引起基流量減小的重要原因。

    5 結(jié)論

    本文以硯瓦川流域?yàn)檠芯繀^(qū)域,探討了不同基流分割方法在研究區(qū)的適用性,分析了流域1981—2016年期間河川基流的年際及年內(nèi)變化規(guī)律,定性和定量的評(píng)價(jià)了變化環(huán)境對(duì)基流變化的影響,主要結(jié)論如下:

    (1)Lyne-Hollick濾波法計(jì)算出來(lái)的流域多年平均BFI值為0.55,與枯水指數(shù)法計(jì)算出來(lái)的多年實(shí)際BFI平均值0.58較為接近,另外該方法的NSE和Re分別為0.49和17%,具有一定的估算精度,且其基流分割結(jié)果比較符合實(shí)際日基流變化規(guī)律。綜合年基流指數(shù)、年基流量誤差分析與日基流過(guò)程線結(jié)果,本研究認(rèn)為L(zhǎng)yne-Hollick法在所有基流分割方法中具有更好的合理性和可靠性。

    (2)基流量在年內(nèi)的變化幅度小于徑流,且對(duì)降水/融雪補(bǔ)給的響應(yīng)時(shí)間也較為緩慢,體現(xiàn)了基流補(bǔ)給的滯后性,而B(niǎo)FI在年內(nèi)總體呈現(xiàn)先減少后增加的變化規(guī)律,汛期(5—9月份)的月BFI值明顯低于非汛期。流域多年平均基流量為0.152 m3/s,基流量和BFI在年際上均呈現(xiàn)極顯著的減少趨勢(shì)(P<0.01),且分別于1993年和2006年附近發(fā)生了突變。

    (3)年潛在蒸散發(fā)的突變時(shí)間(1993年)以及流域治理步伐開(kāi)始加快的時(shí)間(1994年)均與年基流量的突變時(shí)間(1993年)較為接近,這說(shuō)明硯瓦川流域潛在蒸散發(fā)和NDVI的增加可能是河川基流量減小的主要原因;另外,年NDVI的突變時(shí)間與年基流指數(shù)的突變時(shí)間較為接近,這表明NDVI的變化對(duì)于河川基流指數(shù)的減小應(yīng)該起到了一定的作用。另外,根據(jù)Spearman相關(guān)分析結(jié)果,河川基流量與潛在蒸散發(fā)量相關(guān)性最大(R= -0.510),而基流指數(shù)與NDVI的相關(guān)性最高(R= -0.432)。

    (4)流域降水量、潛在蒸散發(fā)量及NDVI的變化對(duì)基流量變化的貢獻(xiàn)率分別為-99.1%、113.3%和85.8%,對(duì)BFI變化的貢獻(xiàn)率分別為41.3%、-27.7%和86.5%,可見(jiàn)潛在蒸散發(fā)和NDVI的變化是引起硯瓦川流域河川基流量變化的主要原因,而流域NDVI的增加確實(shí)對(duì)于河川基流指數(shù)的降低起到了決定性的作用。

    猜你喜歡
    方法研究
    FMS與YBT相關(guān)性的實(shí)證研究
    2020年國(guó)內(nèi)翻譯研究述評(píng)
    遼代千人邑研究述論
    視錯(cuò)覺(jué)在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    學(xué)習(xí)方法
    EMA伺服控制系統(tǒng)研究
    新版C-NCAP側(cè)面碰撞假人損傷研究
    可能是方法不對(duì)
    用對(duì)方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    亚洲av在线观看美女高潮| 在线播放无遮挡| 成人亚洲欧美一区二区av| av播播在线观看一区| 欧美另类一区| 51国产日韩欧美| videos熟女内射| 亚洲色图av天堂| 免费观看a级毛片全部| 有码 亚洲区| 日本一二三区视频观看| 亚洲人成网站高清观看| 国产成人精品久久久久久| 国产精品久久久久久精品电影| 观看美女的网站| 少妇熟女欧美另类| 国产免费又黄又爽又色| 国产高潮美女av| 国产欧美另类精品又又久久亚洲欧美| 久久久久九九精品影院| 搞女人的毛片| 2022亚洲国产成人精品| 激情五月婷婷亚洲| 国产精品一二三区在线看| 精品久久久久久久人妻蜜臀av| 身体一侧抽搐| 久久久久久伊人网av| av一本久久久久| 精品欧美国产一区二区三| 亚洲天堂国产精品一区在线| 日本-黄色视频高清免费观看| 国产成人91sexporn| 一个人观看的视频www高清免费观看| 超碰97精品在线观看| 午夜免费观看性视频| 精品一区二区三区人妻视频| 精品久久久久久久久av| 秋霞在线观看毛片| 国产伦理片在线播放av一区| 哪个播放器可以免费观看大片| 91久久精品国产一区二区三区| 一本一本综合久久| 99久久精品一区二区三区| 国产女主播在线喷水免费视频网站 | 草草在线视频免费看| 91在线精品国自产拍蜜月| 色视频www国产| 国产精品不卡视频一区二区| 少妇人妻一区二区三区视频| 丝袜美腿在线中文| 狠狠精品人妻久久久久久综合| 亚洲在久久综合| av播播在线观看一区| 免费大片18禁| 91精品国产九色| 日韩国内少妇激情av| 2021天堂中文幕一二区在线观| 噜噜噜噜噜久久久久久91| 午夜激情欧美在线| 不卡视频在线观看欧美| 草草在线视频免费看| 天天躁日日操中文字幕| videos熟女内射| 高清午夜精品一区二区三区| 2018国产大陆天天弄谢| 99热全是精品| 免费大片黄手机在线观看| eeuss影院久久| 女人十人毛片免费观看3o分钟| 成年免费大片在线观看| 99久久九九国产精品国产免费| 五月伊人婷婷丁香| 日日摸夜夜添夜夜爱| 亚洲精品国产av蜜桃| 黄色一级大片看看| 国产精品国产三级国产av玫瑰| 午夜免费激情av| 99热这里只有是精品50| 亚洲精品自拍成人| 啦啦啦中文免费视频观看日本| 婷婷色av中文字幕| 亚洲经典国产精华液单| 亚洲18禁久久av| 久久久久精品性色| 亚洲,欧美,日韩| 97人妻精品一区二区三区麻豆| 久久99热这里只频精品6学生| 如何舔出高潮| 夜夜看夜夜爽夜夜摸| 亚州av有码| 五月天丁香电影| 午夜免费观看性视频| 神马国产精品三级电影在线观看| 久久99热这里只频精品6学生| 亚洲精品国产av蜜桃| av福利片在线观看| 精品久久国产蜜桃| 欧美 日韩 精品 国产| 三级经典国产精品| 亚洲av福利一区| 国产精品久久久久久精品电影小说 | 午夜精品在线福利| 国产精品不卡视频一区二区| 日韩不卡一区二区三区视频在线| 五月天丁香电影| 赤兔流量卡办理| 2021天堂中文幕一二区在线观| 啦啦啦中文免费视频观看日本| 欧美高清性xxxxhd video| 国产白丝娇喘喷水9色精品| 国产精品久久久久久精品电影| 亚洲四区av| 成人一区二区视频在线观看| 人妻夜夜爽99麻豆av| 国产伦精品一区二区三区视频9| 欧美激情国产日韩精品一区| 亚洲成人av在线免费| 黄色一级大片看看| 亚洲欧美一区二区三区黑人 | 真实男女啪啪啪动态图| 亚洲av中文字字幕乱码综合| 91精品一卡2卡3卡4卡| 大香蕉久久网| 精品熟女少妇av免费看| 黄片无遮挡物在线观看| 爱豆传媒免费全集在线观看| 国产精品爽爽va在线观看网站| 亚洲欧美日韩卡通动漫| 最近的中文字幕免费完整| 午夜免费激情av| 亚洲国产成人一精品久久久| 欧美日韩综合久久久久久| av又黄又爽大尺度在线免费看| 能在线免费观看的黄片| 国产成人精品福利久久| 一边亲一边摸免费视频| 中文字幕免费在线视频6| 国产男女超爽视频在线观看| 91久久精品国产一区二区三区| 国产69精品久久久久777片| 国产精品一区二区性色av| 国产黄a三级三级三级人| 国产视频内射| 一个人看的www免费观看视频| 国产精品一及| 三级国产精品片| 欧美激情国产日韩精品一区| 菩萨蛮人人尽说江南好唐韦庄| 国产精品熟女久久久久浪| 午夜视频国产福利| 亚洲国产精品专区欧美| 99热6这里只有精品| 免费看光身美女| 美女国产视频在线观看| a级毛色黄片| 三级经典国产精品| 亚洲国产欧美人成| 亚洲图色成人| ponron亚洲| 一个人免费在线观看电影| 99热这里只有精品一区| 国产爱豆传媒在线观看| 两个人视频免费观看高清| 久久久a久久爽久久v久久| 国产免费福利视频在线观看| 国产精品女同一区二区软件| 午夜久久久久精精品| 久久久久久九九精品二区国产| 嫩草影院精品99| 成年av动漫网址| 三级国产精品欧美在线观看| 国产av国产精品国产| 国产色爽女视频免费观看| 日韩av在线免费看完整版不卡| 一个人免费在线观看电影| 老女人水多毛片| 久久精品人妻少妇| 欧美日韩国产mv在线观看视频 | 老师上课跳d突然被开到最大视频| 久久久久久九九精品二区国产| 亚洲精品久久午夜乱码| 欧美潮喷喷水| 精品一区二区免费观看| 天堂俺去俺来也www色官网 | 嫩草影院入口| 中国美白少妇内射xxxbb| 国产伦精品一区二区三区视频9| 亚洲欧美中文字幕日韩二区| 免费av不卡在线播放| 国产精品三级大全| 日日啪夜夜爽| 日韩欧美精品v在线| 亚洲欧美清纯卡通| 啦啦啦韩国在线观看视频| 只有这里有精品99| 好男人在线观看高清免费视频| 伦精品一区二区三区| 国产成人91sexporn| 99热6这里只有精品| 亚洲av二区三区四区| 色哟哟·www| 欧美成人精品欧美一级黄| 永久网站在线| 男女视频在线观看网站免费| 日韩成人伦理影院| av免费在线看不卡| 中文字幕亚洲精品专区| 2022亚洲国产成人精品| 最后的刺客免费高清国语| 人人妻人人看人人澡| 亚洲怡红院男人天堂| 三级毛片av免费| 国产精品1区2区在线观看.| 国产精品三级大全| 国产一区有黄有色的免费视频 | 水蜜桃什么品种好| 美女脱内裤让男人舔精品视频| 日韩强制内射视频| 乱码一卡2卡4卡精品| 亚洲真实伦在线观看| 精品久久久久久久末码| 欧美激情国产日韩精品一区| av一本久久久久| 日韩欧美国产在线观看| 美女主播在线视频| 久久久久性生活片| ponron亚洲| 精品久久久久久久末码| 亚洲欧美日韩东京热| 国产国拍精品亚洲av在线观看| 蜜桃亚洲精品一区二区三区| 亚洲精品日韩在线中文字幕| 男插女下体视频免费在线播放| 麻豆成人av视频| 亚洲熟妇中文字幕五十中出| 亚洲av在线观看美女高潮| 2022亚洲国产成人精品| 国产精品久久久久久精品电影小说 | 干丝袜人妻中文字幕| 亚洲欧美一区二区三区黑人 | 纵有疾风起免费观看全集完整版 | 伦精品一区二区三区| 亚洲18禁久久av| 久久久久久久国产电影| 2021天堂中文幕一二区在线观| 毛片女人毛片| 免费av不卡在线播放| 99九九线精品视频在线观看视频| 成人午夜精彩视频在线观看| 日韩三级伦理在线观看| 99久久人妻综合| 日产精品乱码卡一卡2卡三| 国语对白做爰xxxⅹ性视频网站| 中国美白少妇内射xxxbb| 国产精品一区二区三区四区久久| 午夜久久久久精精品| 哪个播放器可以免费观看大片| 久99久视频精品免费| 国产男女超爽视频在线观看| 色播亚洲综合网| 亚洲国产日韩欧美精品在线观看| 精品久久国产蜜桃| 日韩视频在线欧美| 亚洲一级一片aⅴ在线观看| 极品教师在线视频| 久久99蜜桃精品久久| 好男人在线观看高清免费视频| 一个人看的www免费观看视频| 一级二级三级毛片免费看| 欧美一级a爱片免费观看看| 免费播放大片免费观看视频在线观看| 男人狂女人下面高潮的视频| 久久99精品国语久久久| 97精品久久久久久久久久精品| av.在线天堂| 国产高清有码在线观看视频| av在线亚洲专区| 高清av免费在线| 高清在线视频一区二区三区| 日本wwww免费看| 国产精品一二三区在线看| 国产女主播在线喷水免费视频网站 | 亚洲欧美成人精品一区二区| 97在线视频观看| 国产色爽女视频免费观看| 又大又黄又爽视频免费| 久热久热在线精品观看| 少妇人妻一区二区三区视频| 亚洲一区高清亚洲精品| 国产极品天堂在线| 久99久视频精品免费| 在线a可以看的网站| 18禁在线播放成人免费| 日韩成人伦理影院| av一本久久久久| 夜夜爽夜夜爽视频| 777米奇影视久久| 又爽又黄a免费视频| 黄片wwwwww| 麻豆精品久久久久久蜜桃| 91久久精品电影网| 亚洲精品日本国产第一区| 网址你懂的国产日韩在线| 欧美一区二区亚洲| 国产亚洲精品av在线| 有码 亚洲区| 亚洲人成网站高清观看| 亚洲怡红院男人天堂| 久久综合国产亚洲精品| 亚洲精品国产av蜜桃| 国产精品无大码| 非洲黑人性xxxx精品又粗又长| 我的老师免费观看完整版| 性色avwww在线观看| 亚洲色图av天堂| 一夜夜www| 亚洲精华国产精华液的使用体验| 久久久精品94久久精品| 人人妻人人看人人澡| 少妇的逼好多水| 国产一区有黄有色的免费视频 | 国产精品久久久久久久久免| 亚洲熟女精品中文字幕| 超碰97精品在线观看| 亚洲精品久久久久久婷婷小说| 成年人午夜在线观看视频 | 69av精品久久久久久| 久久精品国产亚洲av涩爱| 国产亚洲最大av| 一级毛片aaaaaa免费看小| 亚洲一级一片aⅴ在线观看| 别揉我奶头 嗯啊视频| 久久精品国产自在天天线| 国产精品国产三级国产专区5o| 久久热精品热| 三级毛片av免费| 男的添女的下面高潮视频| 日韩亚洲欧美综合| 亚洲av日韩在线播放| 国语对白做爰xxxⅹ性视频网站| 国产精品一区二区在线观看99 | 伦精品一区二区三区| 亚洲av电影在线观看一区二区三区 | 国产成人一区二区在线| 欧美性感艳星| 欧美人与善性xxx| av在线天堂中文字幕| 国产午夜精品一二区理论片| 男女国产视频网站| 日日摸夜夜添夜夜添av毛片| 深爱激情五月婷婷| 国产成人freesex在线| 青青草视频在线视频观看| 国产亚洲精品av在线| 搡老乐熟女国产| 小蜜桃在线观看免费完整版高清| 舔av片在线| 久久久欧美国产精品| 一级毛片电影观看| 久久亚洲国产成人精品v| 国产高清国产精品国产三级 | 青春草视频在线免费观看| 亚洲自偷自拍三级| 国产欧美日韩精品一区二区| 久久久精品94久久精品| 亚洲国产日韩欧美精品在线观看| 综合色丁香网| 亚洲最大成人中文| 天天躁日日操中文字幕| 高清午夜精品一区二区三区| 男人和女人高潮做爰伦理| 亚洲av日韩在线播放| 老司机影院成人| 亚洲av日韩在线播放| 精品欧美国产一区二区三| 人妻夜夜爽99麻豆av| 亚洲国产色片| 亚洲精品乱码久久久久久按摩| 能在线免费观看的黄片| 2022亚洲国产成人精品| 韩国高清视频一区二区三区| 深夜a级毛片| 国产精品久久久久久久久免| 草草在线视频免费看| 精品亚洲乱码少妇综合久久| 黑人高潮一二区| 国产一区有黄有色的免费视频 | 尤物成人国产欧美一区二区三区| 亚洲怡红院男人天堂| 直男gayav资源| 搡女人真爽免费视频火全软件| av线在线观看网站| 亚洲美女视频黄频| 又爽又黄无遮挡网站| 两个人的视频大全免费| 久久精品夜色国产| 一区二区三区高清视频在线| 欧美xxxx黑人xx丫x性爽| 黄色配什么色好看| 大片免费播放器 马上看| 国产男女超爽视频在线观看| 2021少妇久久久久久久久久久| 三级国产精品欧美在线观看| 午夜激情福利司机影院| 亚洲精品日韩av片在线观看| 亚洲不卡免费看| 色综合亚洲欧美另类图片| 极品少妇高潮喷水抽搐| 国产av不卡久久| 美女被艹到高潮喷水动态| 一级毛片电影观看| 搞女人的毛片| 好男人视频免费观看在线| 成人亚洲欧美一区二区av| 亚洲欧美清纯卡通| 久久热精品热| 婷婷色av中文字幕| 亚洲国产av新网站| 尤物成人国产欧美一区二区三区| 中文精品一卡2卡3卡4更新| 国产成人精品一,二区| av又黄又爽大尺度在线免费看| 可以在线观看毛片的网站| a级毛片免费高清观看在线播放| 日日摸夜夜添夜夜添av毛片| 色播亚洲综合网| 精品国产露脸久久av麻豆 | 三级男女做爰猛烈吃奶摸视频| 午夜福利在线观看吧| 亚洲电影在线观看av| 男人爽女人下面视频在线观看| 亚洲激情五月婷婷啪啪| 男人狂女人下面高潮的视频| www.色视频.com| 免费少妇av软件| 激情五月婷婷亚洲| 伊人久久精品亚洲午夜| av国产久精品久网站免费入址| 欧美日韩一区二区视频在线观看视频在线 | 淫秽高清视频在线观看| 99re6热这里在线精品视频| 在线观看免费高清a一片| 欧美区成人在线视频| 麻豆成人午夜福利视频| 只有这里有精品99| 高清毛片免费看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美高清成人免费视频www| 一级毛片久久久久久久久女| 国产精品一区二区在线观看99 | 国产精品一区二区在线观看99 | 在线a可以看的网站| 国产精品熟女久久久久浪| 毛片一级片免费看久久久久| 搞女人的毛片| 18禁裸乳无遮挡免费网站照片| av在线蜜桃| 一级黄片播放器| 97人妻精品一区二区三区麻豆| 久久精品国产鲁丝片午夜精品| 亚洲熟妇中文字幕五十中出| 中文资源天堂在线| 国产精品久久久久久av不卡| 国产一区有黄有色的免费视频 | 国产亚洲av片在线观看秒播厂 | 又黄又爽又刺激的免费视频.| 夫妻午夜视频| 免费看光身美女| 18禁裸乳无遮挡免费网站照片| ponron亚洲| 成人二区视频| 免费电影在线观看免费观看| 久久久久久久久中文| 久久97久久精品| 肉色欧美久久久久久久蜜桃 | kizo精华| 啦啦啦中文免费视频观看日本| 熟女电影av网| 99热全是精品| 日韩制服骚丝袜av| 日韩一区二区视频免费看| 亚洲av一区综合| 老司机影院毛片| 乱系列少妇在线播放| 精品人妻视频免费看| 国产在线一区二区三区精| 一边亲一边摸免费视频| 日本熟妇午夜| 久久国内精品自在自线图片| 午夜福利在线在线| 亚洲图色成人| 久久精品综合一区二区三区| 免费大片黄手机在线观看| 大香蕉97超碰在线| 国产亚洲5aaaaa淫片| 精品久久久久久成人av| 国产av不卡久久| 一个人看视频在线观看www免费| 18禁裸乳无遮挡免费网站照片| 精品人妻偷拍中文字幕| 国产精品蜜桃在线观看| 爱豆传媒免费全集在线观看| 国产伦在线观看视频一区| 午夜老司机福利剧场| 男人舔奶头视频| 亚洲人成网站在线播| 欧美精品国产亚洲| 伦精品一区二区三区| 一本久久精品| 汤姆久久久久久久影院中文字幕 | 美女大奶头视频| 国产久久久一区二区三区| 精品久久久久久电影网| 免费黄频网站在线观看国产| 亚洲丝袜综合中文字幕| 日本wwww免费看| 免费黄频网站在线观看国产| freevideosex欧美| 天天躁日日操中文字幕| av天堂中文字幕网| 午夜福利视频精品| 成年女人在线观看亚洲视频 | 亚洲人成网站在线观看播放| 国产精品人妻久久久影院| 成年女人在线观看亚洲视频 | 最近中文字幕高清免费大全6| 伦理电影大哥的女人| 免费观看性生交大片5| 少妇熟女aⅴ在线视频| 国产精品久久久久久av不卡| 99视频精品全部免费 在线| 亚洲最大成人手机在线| 人妻夜夜爽99麻豆av| 亚洲欧洲日产国产| 真实男女啪啪啪动态图| 爱豆传媒免费全集在线观看| 青春草视频在线免费观看| 国产亚洲午夜精品一区二区久久 | 亚洲精品一二三| 久久综合国产亚洲精品| 国产成人freesex在线| 91精品伊人久久大香线蕉| 美女脱内裤让男人舔精品视频| 99久国产av精品国产电影| 永久免费av网站大全| 午夜激情欧美在线| 街头女战士在线观看网站| 亚洲精品456在线播放app| 精品国产三级普通话版| 久久久久久久国产电影| 免费av毛片视频| 日韩视频在线欧美| 亚洲精品国产成人久久av| 高清视频免费观看一区二区 | 精品酒店卫生间| 亚洲国产最新在线播放| 美女脱内裤让男人舔精品视频| 国产精品熟女久久久久浪| 国内少妇人妻偷人精品xxx网站| 久久久久性生活片| 美女主播在线视频| 国产乱人偷精品视频| 日韩欧美三级三区| 欧美潮喷喷水| 国产高清国产精品国产三级 | 不卡视频在线观看欧美| 成人av在线播放网站| 成人国产麻豆网| 免费av观看视频| 最后的刺客免费高清国语| 91aial.com中文字幕在线观看| 夜夜爽夜夜爽视频| 在线观看免费高清a一片| 亚洲av中文av极速乱| 老师上课跳d突然被开到最大视频| 色哟哟·www| 亚洲在久久综合| 午夜福利在线观看吧| 国产一级毛片七仙女欲春2| 日本欧美国产在线视频| 午夜福利高清视频| 韩国高清视频一区二区三区| 亚洲自拍偷在线| 成年av动漫网址| 日韩欧美三级三区| 小蜜桃在线观看免费完整版高清| 一级毛片久久久久久久久女| 国产白丝娇喘喷水9色精品| 国产亚洲精品久久久com| 久久精品人妻少妇| 国产v大片淫在线免费观看| 日日撸夜夜添| 日韩一区二区视频免费看| 欧美日韩视频高清一区二区三区二| 国产精品一区二区三区四区免费观看| 一区二区三区乱码不卡18| 精品人妻偷拍中文字幕| 久热久热在线精品观看| 色综合亚洲欧美另类图片| 欧美潮喷喷水| 一个人免费在线观看电影| 久久久久久久久久人人人人人人| 精品欧美国产一区二区三| 只有这里有精品99| av线在线观看网站| 少妇的逼好多水| 内射极品少妇av片p| 午夜老司机福利剧场| 天天躁夜夜躁狠狠久久av| 国产精品蜜桃在线观看| 久久久久久久久久黄片| 亚洲欧美成人精品一区二区| 国产精品嫩草影院av在线观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品一二三|