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

    基于改進(jìn)Tennant法的河流生態(tài)基流量計算

    2017-03-22 03:07:43楊裕恒曹升樂
    中國農(nóng)村水利水電 2017年8期
    關(guān)鍵詞:沂河小清河補(bǔ)水

    楊裕恒,曹升樂,劉 陽

    (山東大學(xué)土建與水利學(xué)院,濟(jì)南 250061)

    0 引 言

    河流生態(tài)基流量定義為“維持河流中動植物等生物群落穩(wěn)定,河流生態(tài)環(huán)境與功能不受破壞所需要的基本水量[1,2]”,其計算值應(yīng)是一個長期觀測、較為穩(wěn)定的值。根據(jù)Tharme[3]對生態(tài)基流量計算方法的總結(jié),其計算方法大致可分為四類:水文學(xué)方法、水力學(xué)方法、生態(tài)流量分析法和生態(tài)環(huán)境模擬法。目前,水力學(xué)方法、生態(tài)流量分析法和生態(tài)環(huán)境模擬法在我國實際運(yùn)用中易遇到生態(tài)監(jiān)測資料不足、人力物力有限等問題[4],所以國內(nèi)使用較多的方法是水文學(xué)法。該方法核心是應(yīng)用相關(guān)的水文指標(biāo),通過長系列流量資料來推求相應(yīng)的河流生態(tài)需水量、生態(tài)基流量等指標(biāo)。在水文學(xué)方法中,Tennant法是應(yīng)用較廣泛的方法,它是田納特等專家研究了11條受人類活動影響很小的河流,觀測了河流在不同的流量與過流面積等條件下,對漁業(yè)的影響并匯總命名成Tennant法。Tennant法優(yōu)點在于只需獲得當(dāng)?shù)匾欢屋^長時間的流量資料,就能根據(jù)相應(yīng)系數(shù)快捷計算出當(dāng)?shù)睾恿魃鷳B(tài)基流量,且計算結(jié)果較為合理[5]。但它側(cè)重考慮流量年際變化,淡化了河流豐水期與枯水期流量的差距[6],對于不同功能定位河流沒有分類討論。我國北方地區(qū)河流呈現(xiàn)出明顯的季節(jié)性變化,汛期與非汛期流量差異顯著,河流生態(tài)基流也隨季節(jié)顯著變化[7]。受人類活動的影響,我國北方地區(qū)大部分天然河流、城市河流徑流量逐年減少[8],定義這種河流為未補(bǔ)水河流;但也有小部分河流由于調(diào)水、人工補(bǔ)水、河流改道等原因流量逐年增大,定義這種河流為補(bǔ)水河流。鑒于此,本文針對不同河流功能定位與徑流量的變化趨勢,提出了一種對徑流過程還原或還現(xiàn)的方法修正年均流量,并考慮河流年內(nèi)徑流變化以及生態(tài)過程對水文過程的影響,對Tennant法進(jìn)行改進(jìn)。

    1 研究方法

    1.1 Mann-Kendall法

    Mann-Kendall法簡稱M-K法,是一種適用于非正態(tài)分布[9]的無分布檢驗方法。該方法的優(yōu)點在于可以排除數(shù)據(jù)中極端值的干擾,也無需樣本具有一定的統(tǒng)計分布規(guī)律,并且計算簡便[10],在水文學(xué)上廣泛應(yīng)用于研究降雨、徑流等方面序列變異診斷[11]。

    對于某一個樣本容量為n的時間序列X(X1,X2,…,Xn),構(gòu)造它的秩序列:

    (5)

    式中:Dτ為第i時刻時所有大于j時刻數(shù)值的個數(shù)累計值;τ為i的取值范圍,τ=1,2,…,n;n為樣本容量;Ri為i時刻時依據(jù)Xi與Xj的大小取值0或1;Xi為i時刻的數(shù)值;Xj為j時刻的數(shù)值,j=1,2,…,i;UFτ為Dτ標(biāo)準(zhǔn)化計算后得到的值;E(Dτ)、V(Dτ)為Dτ的均值和方差。

    通過查找正態(tài)分布函數(shù)表,可以得出在一個顯著水平為α?xí)r的臨界值Uα,若公式(3)所計算的|UFτ|>Uα,則該序列存在明顯的上升或者下降的變化趨勢。同理,以逆序列{Xn,Xn-1,…,1}重復(fù)以上過程,令UBτ=-UFτ,其中τ=n,n-1,…,1,根據(jù)UFτ和UBτ值繪制UF和UB曲線圖,若UFτ和UBτ值小于0,則序列呈現(xiàn)下降趨勢,反之則為上升趨勢。若UFτ和UBτ曲線相交在臨界值Uα=±1.96之間(α=0.05),則交點對應(yīng)的時刻有可能為序列變異開始時刻。

    1.2 徑流還原與還現(xiàn)計算

    1.2.1 未補(bǔ)水河流的還原計算

    一般情況下受人類活動影響的未補(bǔ)水河流,河流徑流量呈現(xiàn)逐年減少的趨勢。部分北方河流甚至出現(xiàn)斷流,嚴(yán)重危害了河流健康。在此大背景下,導(dǎo)致徑流突變點前后的徑流量不能視為同一個整體。因此,為了保障北方地區(qū)未補(bǔ)水河流的水文情勢不受破壞,避免極端徑流序列資料對計算結(jié)果造成影響,計算生態(tài)基流量時有必要對徑流序列進(jìn)行還原[12]。

    本文采用降雨-徑流關(guān)系法進(jìn)行徑流量的還原計算。通過1.1節(jié)的方法找出河流徑流突變點,找出受人類活動影響較小的降雨-徑流資料,建立降雨徑流關(guān)系。通過建立的降雨-徑流關(guān)系和估算年份的降雨資料,計算出受人類影響較小的下墊面條件下的還原年徑流量,以此完成還原計算。

    1.2.2 補(bǔ)水河流的還現(xiàn)計算

    我國近年來對河流的開發(fā)利用程度提高,尤其是部分城市河流為滿足景觀、生活需求大量補(bǔ)水,徑流量往往較未補(bǔ)水前有明顯的增大,其水文過程已與未補(bǔ)水河流有較大差別。經(jīng)過補(bǔ)水的河流利用還原徑流序列進(jìn)行生態(tài)基流量計算必然會使計算結(jié)果偏小,造成誤差[13]。因此,在計算此類河流生態(tài)基流量之前,需篩選出突變點之后受到人為影響較大的徑流資料,進(jìn)行徑流量還現(xiàn)計算,增加徑流系列資料的可靠性與一致性,使其更符合實際情況。將河流還現(xiàn)計算歸納為,為滿足流域經(jīng)濟(jì)、生態(tài)要求所補(bǔ)水的河流,其核心在于維持河流當(dāng)前的水文情勢過程,因此在計算時應(yīng)把早期下墊面狀況下產(chǎn)生的徑流量修正為現(xiàn)狀(或近期)條件下的徑流量。

    對于徑流還現(xiàn)方法與1.2.1節(jié)還原計算相似,不再贅述。

    1.3 水文過程與生態(tài)過程的耦合的Tennant法

    1.3.1 典型年以及分期百分?jǐn)?shù)的選取改進(jìn)

    本文選用典型年的年內(nèi)流量過程來代替河流多年平均流量,典型年流量過程更能體現(xiàn)流量的年內(nèi)豐枯變化規(guī)律,即使用典型年的Tennant法求出的生態(tài)基流量更為合理[14]。對于補(bǔ)水河流,從突變點后選取典型年,對于未補(bǔ)水河流則相反。根據(jù)北方河流周期變化及生態(tài)需水的特點,10月至次年3月降水較少,河道水量較枯,水生動植物進(jìn)入冬眠期,生態(tài)需水量相對較小[15],將其定義為一般用水期;4月至9月用水量增加,魚類逐漸進(jìn)入產(chǎn)卵期、幼魚生長發(fā)育期,生態(tài)需水量相對較大[5],將其定義為魚類產(chǎn)育期。在一般用水期河流徑流量較小,Tennant法認(rèn)為多年平均流量的10%為河道流量的最低下限[16],故選用多年平均流量的10%作為生態(tài)基流的計算基準(zhǔn)。而在魚類產(chǎn)育期,河流需水量加大,10%的河道流量難以滿足生活、生態(tài)用水的最低要求。參考相關(guān)文獻(xiàn)的設(shè)定標(biāo)準(zhǔn)[17,18],選取4月至9月修正后流量的15%作為河流生態(tài)基流量能較好地滿足河流與生物的最低要求。具體生態(tài)基流量分期百分比的選取如表1所示。

    1.3.2 月均流量改進(jìn)

    魚類產(chǎn)育期按照生態(tài)基流量的定義,在滿足河流不受破壞的同時,還應(yīng)滿足魚類等水生生物能夠正常產(chǎn)卵、繁育等需求。

    表1 生態(tài)基流量分期百分比的選取 m3/s

    對于不同河流,首先需要明確重要保護(hù)魚類繁殖習(xí)性和生境條件,4月上中旬魚類開始產(chǎn)卵,需水量增大;4月下旬-5月大部分魚類產(chǎn)卵,需水量持續(xù)增大。當(dāng)河流水位下降或者河流流速達(dá)不到要求時,魚類會停止產(chǎn)卵行為[19];6-9月水溫較高,魚類索餌旺盛,為生長期,需要洪水將大量的樹葉草枝葉等有機(jī)物帶入河道。若河流沒有特定保護(hù)魚類,則可以參照四大家魚的習(xí)性。4月開始,魚類逐漸開始為產(chǎn)卵進(jìn)行準(zhǔn)備,根據(jù)相關(guān)研究,河流流速0.20~0.25 m/s是魚類產(chǎn)卵所需的最低流速[20,21],魚卵安全漂浮下限為0.25 m/s[19]。8月為魚類產(chǎn)育高峰期,從4月漲水期開始至8月流速最大期,河流流速至少增加0.1 m/s。因此,針對特定保護(hù)魚類產(chǎn)育習(xí)性,應(yīng)將魚類生長繁育規(guī)律與水文情勢變化的耦合在原Tennant法的基礎(chǔ)上增加一個流速條件,滿足魚類生存繁衍要求,計算流程如圖1所示。

    圖1 河流生態(tài)基流量計算流程圖Fig.1 Flow chart of river ecological basic flow

    2 實例應(yīng)用

    2.1 天然河流的生態(tài)基流量計算

    2.1.1 河流介紹

    沂河是山東境內(nèi)發(fā)源的第一大河。隨著沂河流域城市化速度的加快,流域經(jīng)濟(jì)發(fā)展對沂河水資源的依賴程度不斷提高,綜合用水量也不斷上升[22]。沂河流域內(nèi)由于地形與降雨等原因,徑流年內(nèi)波動劇烈,汛期徑流量較大,易造成難以利用的洪水資源;而枯水期徑流量偏小,部分河段甚至出現(xiàn)了斷流,影響了河流健康與流域經(jīng)濟(jì)的發(fā)展。因此,通過沂河徑流年際變化規(guī)律、年內(nèi)分配特征,有針對性的計算沂河生態(tài)基流量,對沂河流域水資源的合理開發(fā)與保護(hù)有著重要意義。

    2.1.2 徑流變異診斷分析與還原計算

    對沂河流域逐年流量序列通過M-K法進(jìn)行變異診斷分析,計算結(jié)果如圖2所示。圖2顯示,曲線的交點在1967年,經(jīng)過查證,60年代沂河修建大量水庫開始對河流流量造成擾動[21],故選取1967年為該階段變異點。從圖2,可以看出變異點之后流量呈現(xiàn)逐年減小的趨勢明顯,因此主要研究變異點之后的流量序列,對其進(jìn)行還原計算。

    圖2 沂河流域流深序列的Mann-Kendall法變異診斷結(jié)果Fig.2 Mann-Kendall variation diagnosis results of runoff series of Yihe River

    同一流域內(nèi),在降水量相同的條件下,下墊面不同,徑流產(chǎn)流量也不同。以圖3中1955-1966年降水-徑流散點為依據(jù),建立降水-徑流關(guān)系式,得到1955-1966年的降水-徑流關(guān)系式為y=0.150x-76.11,序列相關(guān)系數(shù)R2=0.889,表明變異點前沂河受人為活動影響較小,降雨徑流相關(guān)性良好。應(yīng)用此關(guān)系式對突變點后的徑流進(jìn)行模擬,結(jié)果如圖4所示。圖4顯示,突變點之后由于氣候變化、下墊面變化和其他人類活動等原因,實測流量明顯減小。因此,流量還原之后的生態(tài)基流分析,對恢復(fù)河流的天然水文情勢具有一定意義。

    圖3 沂河流域降水徑流線性相關(guān)分析Fig.3 Linear correlation analysis of rainfall and runoff in YiheRiver Basin

    圖4 沂河實測及還原流量序列圖Fig.4 The measured flow and flow reduction sequence diagramin Yihe River

    2.1.3 典型年的Tennant法計算生態(tài)基流量

    根據(jù)還原后的流量序列,計算得1955-2013年沂河平均流量為51.76 m3/s,比修正之前的流量提高了39.9%。通過選取

    突變點前典型年的方法,選取1961年實測逐月流量序列通過同倍比縮放至修正后的多年平均流量,再乘以不同時期所對應(yīng)的流量百分?jǐn)?shù),計算得經(jīng)過還原分析的沂河逐月生態(tài)基流量,見表2、圖5所示。從表2、圖5中可以看出原Tennant法未能滿足魚類產(chǎn)育要求,而選取典型年逐月流量考慮的魚類產(chǎn)育最低流速,更易滿足生態(tài)功能最小需水量的要求。綜合分析,利用流量還原修正年均流量并使用典型年改進(jìn)的Tennant方法能夠更好地反映沂河天然水文情勢,有助于保障河流流速、水深度、河流生物量,保障河流環(huán)境良性存在和發(fā)展,結(jié)果較為合理。

    2.2 補(bǔ)水河流的生態(tài)基流量計算

    2.2.1 河流介紹

    小清河位于山東省中部地區(qū),流經(jīng)省內(nèi)經(jīng)濟(jì)最發(fā)達(dá)的區(qū)域之一,是濟(jì)南市重要的一條獨流入海河流。小清河歷史上曾是一條具有防洪、排澇、灌溉、航運(yùn)、景觀等多功能的人工運(yùn)河。近年來,為滿足小清河污染治理及航運(yùn)的需求,政府通過臥虎山水庫、提取黃河水的方式補(bǔ)充小清河水量[23],小清河黃臺橋斷面結(jié)構(gòu)圖及補(bǔ)水前后年均水位關(guān)系如圖6所示。目前,小清河流量較補(bǔ)水之前有較大提升,并初步滿足通航要求及污染治理要求。對于小清河流域生態(tài)基流量的研究,需主要分析小清河補(bǔ)水之后年際變化規(guī)律、年內(nèi)分配特征。有針對性的計算小清河生態(tài)基流量,對小清河流域污染治理、流域經(jīng)濟(jì)社會發(fā)展有著重要意義。

    表2 修正前后的沂河逐月生態(tài)基流量計算表 m3/s

    圖5 修正前后的沂河逐月生態(tài)基流量曲線圖Tab.5 The monthly ecological base flow curves before and after correction in Yihe River

    圖6 小清河河流結(jié)構(gòu)及水位圖Tab.6 The structure and water level of the Xiaoqing River

    2.2.2 徑流變異診斷分析與還現(xiàn)計算

    采用式(1)-(5)對小清河流域徑流列進(jìn)行變異診斷分析,計算結(jié)果如圖7所示。圖7顯示,曲線在1997年出現(xiàn)第一個交點,在2007年出現(xiàn)第2個交點。經(jīng)查證,小清河于1996年開始進(jìn)行河道整治工程,因此突變點為1997年具有合理性。1997年為徑流的第一個突變點,且與2007年較為接近,故選取1997年為小清河徑流變異點。從圖7可以看出,變異點之后流量呈現(xiàn)逐年增大的趨勢,因此主要研究變異點之前的流量序列,對其進(jìn)行還現(xiàn)計算。

    圖7 小清河流域流量序列的Mann-Kendall法變異診斷結(jié)果Fig.7 Mann-Kendall variation diagnosis results of runoff series of Xiaoqing River

    以第一個突變點后1997-2015年的降水-徑流序列為基準(zhǔn),建立降水-徑流關(guān)系曲線如圖8所示。得到1997-2015年的降水-徑流關(guān)系式為y=0.113 9x0.722 5,序列相關(guān)系數(shù)R2=0.600 4,曲線相關(guān)系數(shù)較低,表明小清河近年受人工補(bǔ)水等人為活動影響較大。應(yīng)用此關(guān)系式對突變點前的徑流進(jìn)行模擬,結(jié)果如圖9所示。圖9顯示,突變點之前由于河流未經(jīng)補(bǔ)水,實測流量明顯減小。因此,流量還原之后的生態(tài)基流分析,對恢復(fù)河流的天然水文情勢具有一定意義。

    圖8 小清河流域降水徑流線性相關(guān)分析Fig.8 Linear correlation analysis of rainfall and runoff in Xiaoqing River Basin

    2.2.3 典型年的Tennant法計算生態(tài)基流量

    根據(jù)還現(xiàn)后的流量序列,計算得1962-2015年小清河平均流量為12.48 m3/s,比修正之前的流量提高了22.8%。通過選取突變點后典型年的方法,選取2007年實測逐月流量序列通過同倍比縮放至修正后的多年平均流量,再乘以對應(yīng)流量百分?jǐn)?shù),與對應(yīng)滿足魚類產(chǎn)育最低流量中的較大值,計算得經(jīng)過還現(xiàn)分析的小清河逐月生態(tài)基流量,見表2、圖10所示。從圖表中可以看出原Tennant法計算值扁平化了小清河實際來水過程,典型年流量過程則體現(xiàn)了補(bǔ)水河流的用水特點,更易滿足流域經(jīng)濟(jì)、社會、自然發(fā)展的需求。綜合分析,利用流量還現(xiàn)技術(shù)改進(jìn)的Tennant方法能夠更好地反映城市河流的用水特點,結(jié)果較為合理。

    圖9 小清河實測及還現(xiàn)流量序列圖Fig.9 TheMeasured flow and flow sequence diagram in Xiaoqing River

    圖10 修正前后的小清河逐月生態(tài)基流量曲線圖Tab.10 The monthly ecological base flow curves before and after correction in Xiaoqing River

    表3 修正前后的小清河逐月生態(tài)基流量計算表 m3/s

    3 結(jié) 語

    本文通過找出河流流量突變點,對不同類型的河流流量進(jìn)行修正,同時通過選取典型年以及考慮魚類產(chǎn)卵期最低流量對現(xiàn)行Tennant法進(jìn)行改進(jìn),主要得出以下結(jié)論:

    (1)以沂河為例,通過還原計算得出修正后的沂河平均流量為51.76 m3/s,比修正之前的流量提高了39.9%。通過修正的Tennant法解決了原始生態(tài)基流量多年平均流量值平坦化了沂河實際來水過程,典型年流量過程更好地體現(xiàn)了研究河段流量豐枯的季節(jié)性,更易滿足“生態(tài)功能最小需水量”這一要求。綜合分析,利用流量還原技術(shù)改進(jìn)的Tennant方法能夠更好地反映天然水文情勢,有助于保障河流流速及水面寬度、水深度等生態(tài)棲息環(huán)境良性存在和發(fā)展。

    (2)以小清河為例,通過還現(xiàn)計算得出修正后的小清河平均流量為12.48 m3/s,比修正之前的流量提高了22.8%,修正后的流量更符合小清河通航、排污的要求。從圖表中可以看出原Tennant法中多年平均流量值削弱了小清河豐水期實際來水過程,改進(jìn)方法能夠更好地體現(xiàn)了補(bǔ)水區(qū)河段的用水特點,更易流域經(jīng)濟(jì)、社會、自然發(fā)展的需求。綜合分析,利用流量還現(xiàn)技術(shù)改進(jìn)的Tennant方法能夠更好地反映城市河流的用水特點。

    [1] 楊志峰, 尹 民, 崔保山.城市生態(tài)環(huán)境需水量研究----理論與方法[J]. 生態(tài)學(xué)報, 2005,25(3):389-396.

    [2] 粟曉玲, 康紹忠. 生態(tài)需水的概念及其計算方法[J]. 水科學(xué)進(jìn)展, 2003,14(6):740-744.

    [3] Tharme R E. A global perspective on environmental flow assessment: emerging trends in the development and application of environmental flow methodologies for rivers[J]. River Research & Applications, 2003,19(5-6):397-441.

    [4] 陳 菡, 邵東國, 吳俊校, 等. 廣東省河道生態(tài)基流定量分析研究[J]. 南水北調(diào)與水利科技, 2011,9(1):92-95.

    [5] 王西琴, 劉昌明, 楊志峰. 生態(tài)及環(huán)境需水量研究進(jìn)展與前瞻[J]. 水科學(xué)進(jìn)展, 2002,13(4):507-514.

    [6] 韋雨婷, 逄 勇, 羅 縉. 清河干流生態(tài)基流計算分析[J]. 水資源與水工程學(xué)報, 2014,25(6):151-155.

    [7] 鄭紅星, 劉昌明, 豐華麗. 生態(tài)需水的理論內(nèi)涵探討[J]. 水科學(xué)進(jìn)展, 2004,15(5):626-633.

    [8] 鄭志宏, 張澤中, 黃 強(qiáng),等. 生態(tài)需水量計算Tennant法的改進(jìn)及應(yīng)用[J]. 四川大學(xué)學(xué)報(工程科學(xué)版), 2010,42(2):34-39.

    [9] 趙 陽,余新曉,鄭江坤,等. 氣候和土地利用變化對潮白河流域徑流變化的定量影響[J]. 農(nóng)業(yè)工程學(xué)報,2012,28(22):252-260.

    [10] 天 霄, 付 強(qiáng), 孟凡香, 等. 三江平原年降水量1959-2013年演變趨勢及突變分析[J]. 中國農(nóng)村水利水電, 2016,(9):201-204.

    [11] 許 浩,雷曉輝,宋萬禎,等. 水庫建設(shè)運(yùn)行對三岔河流域徑流與洪水特征變化影響分析[J]. 中國農(nóng)村水利水電,2016,(10):108-115.

    [12] 金新芽. 徑流還原實用方法研究[D]. 南京:河海大學(xué),2006.

    [13] 孫娟絨. 坪上水庫徑流還現(xiàn)計算分析[J]. 太原理工大學(xué)學(xué)報, 2005,36(5):589-592.

    [14] 黃 強(qiáng), 李 群, 張澤中, 等.計算黃河干流生態(tài)環(huán)境需水Tennant法的改進(jìn)及應(yīng)用[J]. 水動力學(xué)研究與進(jìn)展(A輯), 2007,22(6):774-781.

    [15] 楊裕恒, 曹升樂, 劉 陽, 等. 基于改進(jìn)Tennant法的小清河生態(tài)基流計算[J]. 水資源與水工程學(xué)報, 2016,27(5):97-101.

    [16] 徐宗學(xué), 武 瑋, 于松延. 生態(tài)基流研究:進(jìn)展與挑戰(zhàn)[J]. 水力發(fā)電學(xué)報, 2016,35(4):1-11.

    [17] 田景環(huán), 王 軼, 茹松楠, 等. 基于河流功能的Tennant法改進(jìn)及其應(yīng)用[J]. 人民黃河, 2011,33(11):100-102.

    [18] 李昌文. 基于改進(jìn)Tennant法和敏感生態(tài)需求的河流生態(tài)需水關(guān)鍵技術(shù)研究[D]. 武漢:華中科技大學(xué), 2015.

    [19] 林俊強(qiáng), 彭期冬, 黃真理. 河流魚類魚卵運(yùn)動的水力學(xué)研究進(jìn)展[J]. 水利學(xué)報, 2015,46(7):869-876.

    [20] 陳永柏, 廖文根, 彭期冬, 等.四大家魚產(chǎn)卵水文水動力特性研究綜述[J]. 水生態(tài)學(xué)雜志, 2009,30(2):130-133.

    [21] 柏海霞, 彭期冬, 李 翀, 等. 長江四大家魚產(chǎn)卵場地形及其自然繁殖水動力條件研究綜述[J]. 中國水利水電科學(xué)研究院學(xué)報, 2014,12(3):249-257.

    [22] 薛麗芳,譚海樵. 沂河流域水文特征變化及其驅(qū)動因素[J]. 地理科學(xué)進(jìn)展, 2011,30(11):1 354-1 360.

    [23] 譚永明.濟(jì)南市小清河水質(zhì)評價及環(huán)境需水量研究[D]. 濟(jì)南:山東大學(xué), 2009.

    猜你喜歡
    沂河小清河補(bǔ)水
    沂河濕地生態(tài)系統(tǒng)現(xiàn)狀及生態(tài)恢復(fù)對策
    園林科技(2021年1期)2022-01-19 03:14:02
    《肥皂》等
    董東管道穿越小清河對航道條件的影響分析
    愛運(yùn)動更要會補(bǔ)水
    風(fēng)吹過沂河淌(組詩)
    惱人的干燥
    時尚北京(2018年1期)2018-01-22 20:11:30
    降溫補(bǔ)水
    都市麗人(2017年6期)2017-07-25 20:40:35
    沂河臨沂站洪水預(yù)報影響因素分析
    山東小清河流域黃土的發(fā)現(xiàn)及勘探意義
    山東小清河分洪道治理方案探討
    中國水利(2015年4期)2015-02-28 15:12:22
    日韩一卡2卡3卡4卡2021年| 嫩草影院精品99| 国产熟女xx| 国产1区2区3区精品| 久久青草综合色| 又黄又粗又硬又大视频| 日韩精品中文字幕看吧| 精品熟女少妇八av免费久了| 国产极品粉嫩免费观看在线| 99久久精品国产亚洲精品| 日本一本二区三区精品| 亚洲中文字幕日韩| 天堂动漫精品| 亚洲成人久久爱视频| 黄片播放在线免费| 久久婷婷成人综合色麻豆| 日韩欧美一区二区三区在线观看| a在线观看视频网站| 午夜久久久在线观看| 一本大道久久a久久精品| 国产精品影院久久| 国产黄色小视频在线观看| 欧美日韩福利视频一区二区| 成人精品一区二区免费| 国产精品久久久久久亚洲av鲁大| 极品教师在线免费播放| 国产成人av教育| 国产激情偷乱视频一区二区| 成人国产一区最新在线观看| 亚洲人成77777在线视频| 亚洲人成电影免费在线| 啦啦啦韩国在线观看视频| 亚洲 欧美 日韩 在线 免费| 黄片小视频在线播放| 午夜福利在线观看吧| 午夜福利成人在线免费观看| 手机成人av网站| 国产真人三级小视频在线观看| 久久人妻福利社区极品人妻图片| 精品久久蜜臀av无| 一级a爱视频在线免费观看| avwww免费| 国产成人精品久久二区二区91| 免费看日本二区| 美国免费a级毛片| 国产野战对白在线观看| 久热爱精品视频在线9| 国产成人av激情在线播放| 欧美 亚洲 国产 日韩一| 妹子高潮喷水视频| 亚洲第一青青草原| 青草久久国产| 亚洲 欧美一区二区三区| 成人特级黄色片久久久久久久| 免费av毛片视频| 亚洲精华国产精华精| 99热这里只有精品一区 | 两个人视频免费观看高清| 欧美黄色淫秽网站| 曰老女人黄片| 在线免费观看的www视频| 俺也久久电影网| 亚洲欧洲精品一区二区精品久久久| 一级a爱视频在线免费观看| 两个人免费观看高清视频| 免费人成视频x8x8入口观看| 国产aⅴ精品一区二区三区波| 一级a爱视频在线免费观看| 国产精品久久久人人做人人爽| a级毛片在线看网站| 法律面前人人平等表现在哪些方面| 婷婷精品国产亚洲av| 日本 av在线| 夜夜爽天天搞| 黑人欧美特级aaaaaa片| 身体一侧抽搐| 91成年电影在线观看| 非洲黑人性xxxx精品又粗又长| 中文字幕高清在线视频| 午夜影院日韩av| 欧美性猛交╳xxx乱大交人| 国产亚洲欧美98| 精品卡一卡二卡四卡免费| 国产av在哪里看| 中文字幕另类日韩欧美亚洲嫩草| 国产激情久久老熟女| 白带黄色成豆腐渣| 亚洲国产欧洲综合997久久, | 久久国产乱子伦精品免费另类| 91av网站免费观看| 欧美国产日韩亚洲一区| 一边摸一边做爽爽视频免费| 国产亚洲欧美在线一区二区| 特大巨黑吊av在线直播 | 精品欧美国产一区二区三| 看免费av毛片| 国产精品久久久久久亚洲av鲁大| 成人午夜高清在线视频 | 老汉色av国产亚洲站长工具| aaaaa片日本免费| 久久精品亚洲精品国产色婷小说| 在线永久观看黄色视频| 可以在线观看的亚洲视频| 成人精品一区二区免费| 久久狼人影院| 在线观看午夜福利视频| 中文字幕av电影在线播放| 国产精品久久久久久亚洲av鲁大| 色播亚洲综合网| 国产精品永久免费网站| 天天躁夜夜躁狠狠躁躁| 亚洲欧美一区二区三区黑人| 久久热在线av| 婷婷精品国产亚洲av在线| 欧美激情高清一区二区三区| 亚洲欧洲精品一区二区精品久久久| 动漫黄色视频在线观看| 亚洲中文av在线| 桃色一区二区三区在线观看| 成人国产综合亚洲| 波多野结衣高清作品| 久久香蕉激情| 日韩国内少妇激情av| av欧美777| 日韩大码丰满熟妇| 中文字幕高清在线视频| 草草在线视频免费看| 十八禁人妻一区二区| 国内精品久久久久精免费| 97碰自拍视频| 亚洲国产精品成人综合色| 伦理电影免费视频| 少妇粗大呻吟视频| 天天添夜夜摸| 级片在线观看| 日本五十路高清| а√天堂www在线а√下载| 法律面前人人平等表现在哪些方面| 亚洲 欧美 日韩 在线 免费| 听说在线观看完整版免费高清| 狂野欧美激情性xxxx| 亚洲色图av天堂| 久久中文看片网| 国产麻豆成人av免费视频| 熟妇人妻久久中文字幕3abv| 久久天堂一区二区三区四区| 人人妻人人澡人人看| 欧美 亚洲 国产 日韩一| 免费av毛片视频| 中文字幕人妻熟女乱码| 岛国在线观看网站| 亚洲精品美女久久av网站| 精品无人区乱码1区二区| 欧美性长视频在线观看| 在线免费观看的www视频| 黄色片一级片一级黄色片| 老司机深夜福利视频在线观看| 国产av一区二区精品久久| 男人操女人黄网站| 熟妇人妻久久中文字幕3abv| 亚洲精品美女久久av网站| 成人国语在线视频| 悠悠久久av| 国产伦人伦偷精品视频| 人妻丰满熟妇av一区二区三区| 欧美又色又爽又黄视频| 丝袜人妻中文字幕| 亚洲专区字幕在线| www.999成人在线观看| 精品一区二区三区四区五区乱码| 午夜两性在线视频| 性色av乱码一区二区三区2| 日韩欧美国产在线观看| 无遮挡黄片免费观看| 久久精品91无色码中文字幕| 黑人操中国人逼视频| 日本 av在线| 一级作爱视频免费观看| 黄色成人免费大全| 免费看美女性在线毛片视频| 亚洲在线自拍视频| 一a级毛片在线观看| 制服诱惑二区| 白带黄色成豆腐渣| 国产黄a三级三级三级人| 精品久久蜜臀av无| 国产一区二区三区视频了| 国产三级黄色录像| 国产人伦9x9x在线观看| 国产在线观看jvid| 很黄的视频免费| 欧美精品啪啪一区二区三区| 中文字幕人成人乱码亚洲影| 国产成人精品久久二区二区91| 国产精品永久免费网站| 午夜福利18| 日韩高清综合在线| 欧美 亚洲 国产 日韩一| 在线免费观看的www视频| 桃色一区二区三区在线观看| 99久久久亚洲精品蜜臀av| 搡老妇女老女人老熟妇| 久久 成人 亚洲| 国产不卡一卡二| 后天国语完整版免费观看| 亚洲专区字幕在线| 色综合站精品国产| 午夜成年电影在线免费观看| 国产精品1区2区在线观看.| 丁香欧美五月| 精品一区二区三区四区五区乱码| 高潮久久久久久久久久久不卡| 丝袜在线中文字幕| 99国产综合亚洲精品| 首页视频小说图片口味搜索| 婷婷精品国产亚洲av| 日韩欧美三级三区| 国产成人影院久久av| 一本精品99久久精品77| 亚洲 欧美 日韩 在线 免费| videosex国产| 精品日产1卡2卡| 亚洲三区欧美一区| 成人亚洲精品一区在线观看| 亚洲成人久久爱视频| 国产伦在线观看视频一区| 少妇被粗大的猛进出69影院| 久9热在线精品视频| 91九色精品人成在线观看| av在线播放免费不卡| 亚洲精品久久国产高清桃花| 中文字幕久久专区| 久久久久久人人人人人| 久久久精品欧美日韩精品| 亚洲国产中文字幕在线视频| 一边摸一边做爽爽视频免费| 变态另类丝袜制服| 亚洲免费av在线视频| 精品福利观看| 国产精品一区二区三区四区久久 | 一本精品99久久精品77| 999精品在线视频| 丝袜美腿诱惑在线| 制服人妻中文乱码| 国产爱豆传媒在线观看 | 免费人成视频x8x8入口观看| xxx96com| 欧美亚洲日本最大视频资源| 听说在线观看完整版免费高清| 特大巨黑吊av在线直播 | 中文字幕人妻熟女乱码| 日韩国内少妇激情av| 在线观看一区二区三区| 国产成人影院久久av| 50天的宝宝边吃奶边哭怎么回事| 999久久久国产精品视频| 亚洲一码二码三码区别大吗| 欧美一级a爱片免费观看看 | 女同久久另类99精品国产91| www.999成人在线观看| 老司机午夜福利在线观看视频| 亚洲熟女毛片儿| 天天躁狠狠躁夜夜躁狠狠躁| 18美女黄网站色大片免费观看| 亚洲人成网站高清观看| 日韩欧美国产一区二区入口| 国产成人精品无人区| av片东京热男人的天堂| 一二三四在线观看免费中文在| 一进一出好大好爽视频| 操出白浆在线播放| 一级毛片高清免费大全| 亚洲三区欧美一区| 一个人免费在线观看的高清视频| 校园春色视频在线观看| 夜夜看夜夜爽夜夜摸| 黑丝袜美女国产一区| 成熟少妇高潮喷水视频| 国产午夜福利久久久久久| 女性生殖器流出的白浆| 成熟少妇高潮喷水视频| 婷婷精品国产亚洲av| 岛国视频午夜一区免费看| 黄色片一级片一级黄色片| 老司机深夜福利视频在线观看| 亚洲真实伦在线观看| 久久久水蜜桃国产精品网| 国产麻豆成人av免费视频| 中文字幕av电影在线播放| 欧美日韩黄片免| 亚洲久久久国产精品| 一级作爱视频免费观看| 一区二区三区国产精品乱码| 亚洲自拍偷在线| 久久午夜亚洲精品久久| 国产成人精品久久二区二区免费| 午夜福利18| 亚洲aⅴ乱码一区二区在线播放 | 两个人看的免费小视频| 久久久久久免费高清国产稀缺| 19禁男女啪啪无遮挡网站| 别揉我奶头~嗯~啊~动态视频| 男男h啪啪无遮挡| 亚洲精品一区av在线观看| 91国产中文字幕| 欧美激情久久久久久爽电影| 在线永久观看黄色视频| 两个人视频免费观看高清| 热99re8久久精品国产| 两个人看的免费小视频| 欧美精品亚洲一区二区| 欧美激情 高清一区二区三区| 一夜夜www| 一进一出抽搐动态| 国产在线精品亚洲第一网站| 亚洲国产精品合色在线| av欧美777| 亚洲精品久久国产高清桃花| 亚洲第一欧美日韩一区二区三区| 国产午夜福利久久久久久| 久久久久久大精品| 日韩欧美国产在线观看| 欧美日本亚洲视频在线播放| 欧美黑人巨大hd| 久久国产亚洲av麻豆专区| 日韩 欧美 亚洲 中文字幕| 一级作爱视频免费观看| 久久久久国产精品人妻aⅴ院| 国产成人系列免费观看| 18禁黄网站禁片午夜丰满| 变态另类成人亚洲欧美熟女| 免费高清视频大片| 国产亚洲欧美精品永久| 99久久久亚洲精品蜜臀av| 一夜夜www| 亚洲 欧美一区二区三区| 夜夜看夜夜爽夜夜摸| 国产精品爽爽va在线观看网站 | 国产精品一区二区精品视频观看| 国产精品自产拍在线观看55亚洲| 国产成人欧美在线观看| 国产精品免费视频内射| 一进一出好大好爽视频| 精品一区二区三区四区五区乱码| 侵犯人妻中文字幕一二三四区| 亚洲精品色激情综合| 狂野欧美激情性xxxx| 亚洲一区二区三区色噜噜| 黑丝袜美女国产一区| 精品国产乱子伦一区二区三区| 男人舔女人下体高潮全视频| 欧美人与性动交α欧美精品济南到| 亚洲自偷自拍图片 自拍| 亚洲,欧美精品.| 在线av久久热| 日本免费一区二区三区高清不卡| av在线天堂中文字幕| 正在播放国产对白刺激| www.自偷自拍.com| 麻豆av在线久日| 午夜免费鲁丝| 国产成人系列免费观看| bbb黄色大片| 欧美日韩精品网址| 免费高清在线观看日韩| 高清在线国产一区| 免费看a级黄色片| 久久性视频一级片| 欧美性长视频在线观看| 国产精品日韩av在线免费观看| svipshipincom国产片| 波多野结衣av一区二区av| 午夜精品在线福利| 又黄又爽又免费观看的视频| 久久久久久久久久黄片| 国产亚洲欧美98| 亚洲成人免费电影在线观看| 久久久国产精品麻豆| 校园春色视频在线观看| 伊人久久大香线蕉亚洲五| 满18在线观看网站| 久久久久久久久中文| 国产精品,欧美在线| tocl精华| 午夜免费鲁丝| 久久久国产欧美日韩av| 香蕉国产在线看| 人妻丰满熟妇av一区二区三区| 午夜视频精品福利| 成人亚洲精品一区在线观看| 精品一区二区三区四区五区乱码| netflix在线观看网站| 久久久久亚洲av毛片大全| 久久中文字幕一级| 中出人妻视频一区二区| 亚洲色图av天堂| 少妇被粗大的猛进出69影院| 高清毛片免费观看视频网站| 夜夜躁狠狠躁天天躁| 久久久久久人人人人人| a级毛片a级免费在线| 精品久久久久久久久久久久久 | 日日干狠狠操夜夜爽| 夜夜躁狠狠躁天天躁| 精品国产一区二区三区四区第35| 久久香蕉激情| 欧美乱色亚洲激情| 99在线人妻在线中文字幕| 国产精品久久电影中文字幕| 老熟妇仑乱视频hdxx| 亚洲精品中文字幕一二三四区| 国产激情偷乱视频一区二区| 丰满的人妻完整版| 亚洲av电影不卡..在线观看| av中文乱码字幕在线| 亚洲一卡2卡3卡4卡5卡精品中文| 手机成人av网站| 亚洲人成77777在线视频| 欧美最黄视频在线播放免费| 国产成年人精品一区二区| 亚洲第一电影网av| 一区福利在线观看| 狂野欧美激情性xxxx| 久久欧美精品欧美久久欧美| 男女之事视频高清在线观看| 在线观看舔阴道视频| 黄色视频不卡| 久久热在线av| 日韩精品免费视频一区二区三区| 国产成人欧美在线观看| 日韩免费av在线播放| 香蕉国产在线看| 午夜两性在线视频| 精品久久久久久久末码| 91九色精品人成在线观看| 长腿黑丝高跟| 久久草成人影院| 欧美亚洲日本最大视频资源| 满18在线观看网站| 国产成人精品久久二区二区免费| 精品国产美女av久久久久小说| 两人在一起打扑克的视频| 亚洲成人久久性| 午夜精品在线福利| 亚洲中文字幕日韩| 日韩av在线大香蕉| 91九色精品人成在线观看| 老汉色av国产亚洲站长工具| 国产成人啪精品午夜网站| 伊人久久大香线蕉亚洲五| 不卡av一区二区三区| 女生性感内裤真人,穿戴方法视频| 一级毛片高清免费大全| 最近最新中文字幕大全免费视频| 中文亚洲av片在线观看爽| 国产欧美日韩一区二区精品| 欧美日本亚洲视频在线播放| 日本三级黄在线观看| 久久99热这里只有精品18| 精品一区二区三区av网在线观看| 久久久国产精品麻豆| 欧美绝顶高潮抽搐喷水| 中国美女看黄片| 一区二区三区高清视频在线| www国产在线视频色| 一区二区日韩欧美中文字幕| 午夜免费鲁丝| 日韩欧美免费精品| 国产精品免费视频内射| 亚洲人成电影免费在线| 18禁黄网站禁片免费观看直播| svipshipincom国产片| 一二三四社区在线视频社区8| 亚洲av成人一区二区三| 久久久久久久精品吃奶| 亚洲精品国产区一区二| aaaaa片日本免费| av欧美777| 夜夜躁狠狠躁天天躁| 亚洲性夜色夜夜综合| 亚洲色图av天堂| 精品欧美国产一区二区三| 亚洲真实伦在线观看| 亚洲国产精品久久男人天堂| 国产黄a三级三级三级人| 欧美乱色亚洲激情| 99久久精品国产亚洲精品| 99re在线观看精品视频| 久久人妻av系列| 欧美性猛交黑人性爽| 国产免费男女视频| 一区二区三区精品91| 久久欧美精品欧美久久欧美| 日韩欧美在线二视频| 哪里可以看免费的av片| 99国产精品一区二区三区| 午夜免费激情av| 久久久久亚洲av毛片大全| www.精华液| www.999成人在线观看| 一本综合久久免费| 嫁个100分男人电影在线观看| 男人舔奶头视频| 亚洲国产中文字幕在线视频| 国产午夜精品久久久久久| 欧美精品啪啪一区二区三区| 亚洲色图 男人天堂 中文字幕| 人人妻人人澡欧美一区二区| 久久久久精品国产欧美久久久| 琪琪午夜伦伦电影理论片6080| 女性被躁到高潮视频| 在线观看免费日韩欧美大片| 久久人人精品亚洲av| 午夜久久久久精精品| 精品久久久久久久毛片微露脸| av免费在线观看网站| 可以在线观看的亚洲视频| 国产一区在线观看成人免费| 黄色a级毛片大全视频| 制服丝袜大香蕉在线| 亚洲 欧美 日韩 在线 免费| 久久久久久亚洲精品国产蜜桃av| 黄片小视频在线播放| 久久久久免费精品人妻一区二区 | 久久青草综合色| 亚洲欧美一区二区三区黑人| 麻豆久久精品国产亚洲av| 老司机午夜福利在线观看视频| 88av欧美| 日本熟妇午夜| 1024手机看黄色片| 美女免费视频网站| 欧美最黄视频在线播放免费| 69av精品久久久久久| 日韩欧美 国产精品| 午夜福利成人在线免费观看| 国产一区二区在线av高清观看| 婷婷亚洲欧美| bbb黄色大片| 婷婷六月久久综合丁香| 一级毛片高清免费大全| 午夜福利在线观看吧| 久久欧美精品欧美久久欧美| 午夜久久久久精精品| 欧美性猛交╳xxx乱大交人| 午夜两性在线视频| 国产精品亚洲av一区麻豆| 国产精品,欧美在线| 国产伦在线观看视频一区| 午夜福利视频1000在线观看| 亚洲国产精品sss在线观看| 精品久久久久久久久久免费视频| 男女之事视频高清在线观看| 老司机靠b影院| 校园春色视频在线观看| videosex国产| 久久香蕉精品热| 国产精品亚洲av一区麻豆| 男女床上黄色一级片免费看| 亚洲最大成人中文| 在线观看免费午夜福利视频| 国产av一区二区精品久久| 精品国产国语对白av| 十八禁人妻一区二区| 波多野结衣av一区二区av| 高潮久久久久久久久久久不卡| 亚洲一码二码三码区别大吗| 亚洲专区中文字幕在线| 日本三级黄在线观看| 亚洲人成77777在线视频| 亚洲三区欧美一区| 久久草成人影院| 亚洲精品粉嫩美女一区| 一二三四在线观看免费中文在| xxxwww97欧美| 十八禁网站免费在线| 亚洲午夜精品一区,二区,三区| 在线观看免费午夜福利视频| 在线永久观看黄色视频| 国产精品亚洲美女久久久| 欧美色视频一区免费| 国产主播在线观看一区二区| 精品国产国语对白av| 久久九九热精品免费| 美女免费视频网站| 99国产精品99久久久久| www.www免费av| 国产精品久久久av美女十八| 国产99久久九九免费精品| 亚洲精华国产精华精| 国产亚洲精品第一综合不卡| 欧美一级a爱片免费观看看 | 国产精品 欧美亚洲| av在线天堂中文字幕| 国产一区二区三区视频了| 欧美中文综合在线视频| 露出奶头的视频| 国产午夜福利久久久久久| 亚洲avbb在线观看| 国产一区二区三区在线臀色熟女| 又紧又爽又黄一区二区| 在线观看免费午夜福利视频| 亚洲国产精品999在线| 国产成人av激情在线播放| 啪啪无遮挡十八禁网站| 韩国精品一区二区三区| 成人三级做爰电影| 日本精品一区二区三区蜜桃| 国产成+人综合+亚洲专区| 国产真实乱freesex| 精品国产美女av久久久久小说| 不卡av一区二区三区| 热99re8久久精品国产| 美女高潮喷水抽搐中文字幕| 成年版毛片免费区|