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

    黃河流域典型三級區(qū)蒸發(fā)序列非一致性分析

    2021-09-05 08:29:06張彩玲魯帆方文堵傳聰
    人民黃河 2021年8期
    關(guān)鍵詞:黃河流域

    張彩玲 魯帆 方文 堵傳聰

    摘 要:為研究蒸發(fā)序列非一致性,增強水文頻率計算的穩(wěn)定性,以黃河流域典型三級區(qū)蒸發(fā)序列為研究對象,以GAMLSS模型為基礎(chǔ),選擇黃河流域5個典型水資源三級區(qū),引入時間t作為解釋變量進行參數(shù)擬合,分析各三級區(qū)年、汛期、非汛期蒸發(fā)序列的非一致性特征。采用Kendall法進行趨勢檢驗,Mann-Kendall法和5 a滑動T檢驗法進行突變檢驗。分析結(jié)果表明:①伊洛河、三門峽—小浪底區(qū)間、花園口以下干流區(qū)間的蒸發(fā)序列突變較為明顯,大夏河與洮河、湟水蒸發(fā)序列的突變存在一定的差異;②5個三級區(qū)不同時間尺度的蒸發(fā)序列以Generalized Gamma為最優(yōu)分布的最多,Logistics的最少;③5個三級區(qū)不同時間尺度的蒸發(fā)序列除部分序列外,均與時間t成負(fù)相關(guān),大部分蒸發(fā)數(shù)據(jù)均包含在置信水平為5%和95%的置信區(qū)間內(nèi),表明不同時間尺度的蒸發(fā)序列的GAMLSS模型擬合效果較好;④5個三級區(qū)蒸發(fā)序列非一致性特征較為明顯,以花園口以下干流區(qū)間年、非汛期蒸發(fā)序列的非一致性特征最為明顯。

    關(guān)鍵詞:GAMLSS;非一致性;蒸發(fā)序列;三級區(qū);黃河流域

    中圖分類號:P333;TV882.1 文獻標(biāo)志碼:A

    doi:10.3969/j.issn.1000-1379.2021.08.004

    引用格式:張彩玲,魯帆,方文,等.黃河流域典型三級區(qū)蒸發(fā)序列非一致性分析[J].人民黃河,2021,43(8):19-25.

    Abstract: In order to study the inconsistency of evaporation series and increase the certainty of hydrological frequency calculation, this paper took the evaporation series of typical third-level regions in the Yellow River Basin as the research object. Based on the GAMLSS model, it selected five typical water resources third-level regions in the Yellow River Basin, introduced the time t as an explanatory variable to fit the parameters, and analyzed the non-uniform characteristics of the evaporating sequences in the three-stage regions of the year, flood season and non-flood period. The Kendall method was used for trend test and the mutation test was performed by using Mann-kendall and 5-year sliding T method. The results show that: a) the evaporation sequence of the Yiluo River, Sanmenxia-Xiaolangdi reach and the main stream below the Huayuankou is more obvious and the mutations of evaporation sequence of Daxia, Taohe and Huangshui Rivers are different; b) the evaporation sequence at different time scales in the third-level region is the most optimal distribution with the most Generalized Gamma, and the Logistics is the least; c) the evaporation sequence on the different time scales of the five tertiary regions is negatively correlated with the time t except for some sequences. Most of the evaporation data are included in the confidence interval of 5% and 95% confidence level, indicating that the GAMLSS model of the evaporation sequence on different time scales has a better fitting effect; d) the five third-level evaporation sequences are not consistent. The sexual characteristics are more obvious. The non-uniform characteristics of the annual and non-flood season evaporation sequences in the main stream below the Huayuankou are the most obvious.

    Key words: GAMLSS; non-stationary; evaporation sequence; tertiary zone; Yellow River Basin

    水文頻率分析的一個重要前提是水文資料序列具有一致性,屬于同一分布。氣候變化和人類活動使黃河流域蒸發(fā)序列發(fā)生了不同程度的變化[1],致其不再滿足一致性的假設(shè),基于一致性假設(shè)理論的水文頻率計算方法不再適用。為進一步反映環(huán)境變化對水文頻率計算的影響,需要對水文序列進行非一致性分析,這對減小水文頻率分析結(jié)果的不確定性具有重要意義[2]。國內(nèi)外專家學(xué)者在對水文序列進行計算時會將非一致性假設(shè)條件考慮在內(nèi)[3],并對非一致性假設(shè)條件下水文頻率的計算模型[4]與方法進行研究[5]。魯帆等[6-7]利用Metropolis-Hastings算法,將廣義極值分布(GEV)模型和廣義帕雷托分布(GPD)模型應(yīng)用于大渡河流域。熊立華等[8-9]引入時變矩模型及基于P-Ⅲ型分布的時變矩模型,對非一致性進行歸因分析。但水文序列的計算模型對序列復(fù)雜程度和長度有一定要求,各具針對性。而GAMLSS模型包含多個分布函數(shù)及相應(yīng)的程序包,應(yīng)用范圍較廣。顧西輝等[10]將GAMLSS模型應(yīng)用于珠江流域,分析了洪水極值序列的平穩(wěn)性。張冬冬等[11]將GAMLSS模型應(yīng)用于大渡河流域,分析了降水頻率非一致性。熊斌等[12]應(yīng)用GAMLSS模型對渭河華縣站的枯水頻率進行了研究。江聰?shù)萚13]采用GAMLSS模型對宜昌站的年平均流量和年最小月流量序列的趨勢進行了研究。

    筆者選取黃河流域上游、中游、下游共計5個典型水資源三級區(qū)的面蒸發(fā)序列為研究數(shù)據(jù),采用Kendall法檢驗序列的趨勢和顯著性,以Mann-Kendall法和滑動T檢驗法(步長為5)檢驗序列的突變性[14],采用GAMLSS模型[15]檢驗序列的非一致性,比較不同統(tǒng)計建模方法的計算結(jié)果,以期為后續(xù)黃河流域典型水資源三級區(qū)蒸發(fā)序列的變化規(guī)律研究提供理論依據(jù)。

    1 GAMLSS模型

    1.1 模型定義

    1.2 評價準(zhǔn)則

    以最小GD(Global Deviance)值對應(yīng)的模型為最優(yōu)模型。AIC(Akaike Information Criterion)和SBC(Schwarz Bayesian Criterion)準(zhǔn)則[16]用來判別序列適用一致性還是非一致性模型,若非一致性模型的AIC和SBC值小于一致性模型擬合值,則表明GAMLSS模型更傾向采用非一致性模型對序列進行擬合,序列體現(xiàn)出非一致性特征。

    采用概率點據(jù)相關(guān)系數(shù)法,檢驗正態(tài)標(biāo)準(zhǔn)化的殘差序列ri是否服從標(biāo)準(zhǔn)正態(tài)分布,該方法是由Filliben[17]在1975年提出并廣泛應(yīng)用于統(tǒng)計學(xué)中。Filliben值越接近1,則待檢驗序列越接近標(biāo)準(zhǔn)正態(tài)分布。

    1.3 分布類型

    本文采用Gumbel、Gamma、Logistic、Generalized Gamma四種分布函數(shù),選擇一次函數(shù)、三次樣條函數(shù)(cs,自由度分別為1、2、3)、冪次分段多項式函數(shù)(bfp,自由度分別為1、2、3)三個函數(shù)[18],共7個表達式。以SBC值最小為準(zhǔn)則選擇最優(yōu)表達式。

    2 研究區(qū)域及研究數(shù)據(jù)

    2.1 典型區(qū)域選擇

    黃河流域有29個水資源三級區(qū),根據(jù)蒸發(fā)情況和氣候特征等因素選擇了大夏河與洮河、湟水、伊洛河、三門峽—小浪底區(qū)間、花園口以下干流區(qū)間5個三級區(qū)。

    大夏河與洮河中,大夏河海拔為2 500 m以上,屬寒冷溫潤氣候區(qū),洮河海拔為1 629~4 260 m,屬高原氣候區(qū);湟水海拔為1 727~4 520 m,屬于典型的大陸性季風(fēng)氣候區(qū),其氣候垂直變化明顯,越向上游蒸發(fā)量越小,年平均氣溫0.6~7.9 ℃;伊洛河是伊河和洛河的簡稱,伊河屬暖溫帶大陸性季風(fēng)氣候(半濕潤)區(qū),洛河海拔為169~2 449 m,屬暖溫帶季風(fēng)氣候區(qū);三門峽—小浪底區(qū)間屬亞熱帶和溫帶過渡地帶;花園口以下干流區(qū)間屬于北溫帶大陸性季風(fēng)氣候區(qū),年平均氣溫為14.0~14.3 ℃。

    2.2 研究數(shù)據(jù)分析

    本文所用的蒸發(fā)皿日蒸發(fā)量數(shù)據(jù)

    資料來源于中國氣象局氣候中心和黃河水利委員會。選用1963—2013年的蒸發(fā)皿蒸發(fā)數(shù)據(jù)作為基礎(chǔ)數(shù)據(jù),研究黃河流域蒸發(fā)序列的非一致性演變規(guī)律。

    基礎(chǔ)數(shù)據(jù)由Matlab處理,GAMLSS模型的計算是基于R軟件平臺的GAMLSS程序包,在RStudio軟件中完成的。

    3 計算結(jié)果

    3.1 趨勢及突變檢驗

    采用Kendall法檢驗大夏河與洮河、湟水、伊洛河、三門峽—小浪底區(qū)間和花園口以下干流區(qū)間5個三級區(qū)年、汛期(6—10月)、非汛期(11—5月)的面蒸發(fā)序列的變化趨勢。結(jié)果顯示,5個三級區(qū)不同時間尺度上的面蒸發(fā)量均呈顯著性變化,除大夏河與洮河的蒸發(fā)量呈顯著增加外,其他4個三級區(qū)的蒸發(fā)量均呈顯著減少狀態(tài)。

    采用Mann-Kendall(M-K)檢驗法和滑動T檢驗法(步長為5)檢驗5個三級區(qū)不同時間尺度上蒸發(fā)序列的突變情況,見表1。M-K檢驗法屬于定性的突變檢驗方法,統(tǒng)計量u0.05=1.96;滑動T檢驗法屬于定量的突變檢驗方法,統(tǒng)計量t0.05=2.31。兩種檢驗方法各有優(yōu)勢[19],所以本文綜合考慮來確定突變年份,兩種檢驗方法檢驗結(jié)果一致時以檢驗結(jié)果為準(zhǔn),檢驗結(jié)果不一致時,綜合考慮兩種檢驗方法及不同時間尺度蒸發(fā)序列的變化特征。檢驗結(jié)果顯示,伊洛河、三門峽—小浪底區(qū)間、花園口以下干流區(qū)間的蒸發(fā)序列突變較為明顯,突變多發(fā)生在2000年左右,而大夏河與洮河、湟水蒸發(fā)序列的突變存在一定的差異。大夏河與洮河年蒸發(fā)序列在1997年以前一直處于平緩變化狀態(tài),于1997年開始進入急劇升高狀態(tài);汛期蒸發(fā)序列在1965—1981年之間處于升高—平緩變化狀態(tài),1995年之后趨于平緩變化狀態(tài),于1986年發(fā)生低蒸發(fā)量區(qū)向相對高蒸發(fā)量區(qū)的突變;非汛期蒸發(fā)序列于1995年之前處于平緩變化狀態(tài),1995—1998年處于蒸發(fā)量減少的狀態(tài),以1998年為節(jié)點發(fā)生急劇增加的變化。湟水年蒸發(fā)序列于1973年發(fā)生高蒸發(fā)量區(qū)向相對低蒸發(fā)量區(qū)的突變;汛期蒸發(fā)序列突變不明顯;非汛期蒸發(fā)序列呈升高—急降—急升—下降的變化狀態(tài),于1976年發(fā)生突變。

    3.2 一致性分析

    對5個三級區(qū)在不同時間尺度上的蒸發(fā)序列采用一致性模型進行擬合,將概率密度函數(shù)中的參數(shù)設(shè)為固定值,根據(jù)選取的4種備用分布函數(shù)擬合5個三級區(qū)的蒸發(fā)序列,根據(jù)4種備用分布函數(shù)計算所得的全局?jǐn)M合偏差GD值優(yōu)選各三級區(qū)年、汛期、非汛期蒸發(fā)序列的分布函數(shù)類型,結(jié)果見表2。

    由表2可看出:對于年蒸發(fā)序列,上游的大夏河與洮河、湟水均以Logistics為最優(yōu)分布,下游的花園口以下干流區(qū)間以Generalized Gamma為最優(yōu)分布,中游的伊洛河以Gumbel為最優(yōu)分布、三門峽—小浪底區(qū)間以Generalized Gamma為最優(yōu)分布;對于汛期蒸發(fā)序列,大夏河與洮河、三門峽—小浪底區(qū)間均以Generalized Gamma為最優(yōu)分布,伊洛河、花園口以下干流區(qū)間均以Gumbel為最優(yōu)分布,湟水以Logistics為最優(yōu)分布;對于非汛期蒸發(fā)序列,除大夏河與洮河以Gumbel為最優(yōu)分布外,其他4個三級區(qū)均以Generalized Gamma為最優(yōu)分布。總體而言,5個三級區(qū)不同時間尺度上的蒸發(fā)序列以Generalized Gamma為最優(yōu)分布的最多,Gumbel次之,Logistics最少。

    3.3 非一致性分析

    對黃河流域的5個三級區(qū)不同時間尺度上的蒸發(fā)序列的擬合結(jié)果選擇適合序列的最優(yōu)分布函數(shù),采用以時間t為解釋變量的非一致性模型,擬合相關(guān)參數(shù)成果見表3,重現(xiàn)水平擬合成果如圖1所示,相關(guān)系數(shù)成果見表4。

    由表3可看出,位置參數(shù)與時間t的函數(shù)關(guān)系種類較多。對于年蒸發(fā)序列,除大夏河與洮河、花園口以下干流區(qū)間與時間t建立了自由度為3的冪次分段多項式函數(shù)外,其他3個三級區(qū)均與時間t建立自由度為1~3的三次樣條函數(shù);另外,除大夏河與洮河與時間t成正相關(guān)關(guān)系外,其他4個三級區(qū)均與時間t成負(fù)相關(guān)關(guān)系。對于汛期蒸發(fā)序列,5個三級區(qū)均與時間t成負(fù)相關(guān)關(guān)系;除花園口以下干流區(qū)間與時間建立冪次多項式函數(shù),其他4個三級區(qū)均與時間t建立了自由度為1~3的三次樣條函數(shù)關(guān)系。對于非汛期蒸發(fā)序列,其代表函數(shù)比較單一,均與時間t建立了自由度為1~3的三次樣條函數(shù),且5個三級區(qū)均與時間t成負(fù)相關(guān)關(guān)系。

    由圖1可看出,除伊洛河的年蒸發(fā)量點據(jù)有7個數(shù)據(jù)點位于5%和95%分位數(shù)線外,其余三級區(qū)蒸發(fā)量點據(jù)落入5%和95%分位數(shù)線外的不多于5個。對于年蒸發(fā)量點據(jù),除三門峽—小浪底區(qū)間、花園口以下干流區(qū)間2個三級區(qū)于50%分位數(shù)線周圍相對較散之外,其他3個三級區(qū)均相對較為集中于50%分位數(shù)線周圍。另外,花園口以下干流區(qū)間有5個數(shù)據(jù)點位于95%分位數(shù)線上。對于汛期蒸發(fā)序列,除三門峽—小浪底區(qū)間的蒸發(fā)數(shù)據(jù)點相對較散之外,其余三級區(qū)均相對較集中,重現(xiàn)水平的擬合效果較好。非汛期蒸發(fā)序列,5個三級區(qū)的擬合效果都相對較好,除大夏河與洮河外,其余三級區(qū)均有不多于5個數(shù)據(jù)點位于95%分位數(shù)線上。綜合來說,5個三級區(qū)大部分蒸發(fā)數(shù)據(jù)均包含在置信水平為5%和95%的置信區(qū)間內(nèi),僅有極少數(shù)數(shù)據(jù)位于置信區(qū)間之外,表明非一致性擬合效果較好。序列擬合圖的分位數(shù)線為曲線,表明蒸發(fā)序列跨度較大,非一致性特征更為明顯。

    分析蒸發(fā)序列的殘差序列的分布情況來檢驗非一致性模型是否合理,對5個三級區(qū)蒸發(fā)序列的GAMLSS模型進行殘差檢驗。對殘差序列進行標(biāo)準(zhǔn)化處理,采用PPCC法對處理后的序列進行檢驗,計算蒸發(fā)序列最優(yōu)分布?xì)埐畹那?階中心矩值及Filliben值,計算結(jié)果見表4,不同時間尺度上5個三級區(qū)Filliben值在0.926~0.996之間,接近于1,說明殘差序列服從標(biāo)準(zhǔn)正態(tài)分布。

    殘差序列正態(tài)分布圖如圖2所示,由最優(yōu)分布?xì)埐钚蛄袛M合圖可直觀看出殘差序列的分布情況,5個三級區(qū)的正態(tài)殘差序列均沿紅色曲線且分布在理論曲線的上、下限區(qū)間內(nèi)。大夏河與洮河年蒸發(fā)序列有4~5個正態(tài)殘差值點據(jù)位于理論曲線下限之下,非汛期有2個正態(tài)殘差值點據(jù)位于理論曲線下限之下;湟水年蒸發(fā)序列有1個正態(tài)殘差值點據(jù)位于理論曲線上限之上,汛期蒸發(fā)序列有2個點據(jù)位于理論曲線之外。其他各蒸發(fā)序列均位于理論曲線之內(nèi)。但綜合來看,5個三級區(qū)不同時間尺度上的蒸發(fā)序列的非一致性模型擬合效果較好。

    4 結(jié)果分析

    5個三級區(qū)不同時間尺度上的蒸發(fā)序列通過考慮統(tǒng)計參數(shù)的非一致性模型和不考慮統(tǒng)計參數(shù)的一致性模型的計算,得到兩組GD、AIC、SBC值,其擬合效果對比情況見表5。由表5可以看出,5個三級區(qū)不同時間尺度的蒸發(fā)序列的全局?jǐn)M合偏差GD、AIC和SBC值擬合效果均比不考慮統(tǒng)計參數(shù)非一致性特征的情況下有所改善。

    對于年蒸發(fā)序列,以SBC值為例,改善效果最明顯的三級區(qū)是花園口以下干流區(qū)間,為44.40,大夏河與洮河改善幅度最小,為7.66,其改善效果明顯程度由大到小排序:花園口以下干流區(qū)間、三門峽—小浪底區(qū)間、伊洛河、湟水、大夏河與洮河。以AIC值為例,改善效果最明顯的三級區(qū)是花園口以下干流區(qū)間,為46.33,大夏河與洮河改善幅度最小,為9.59,其改善效果明顯程度由大到小排序:花園口以下干流區(qū)間、三門峽—小浪底區(qū)間、伊洛河、湟水、大夏河與洮河。

    對于汛期蒸發(fā)序列,以SBC值為例,改善效果最明顯的三級區(qū)是伊洛河,為31.24,大夏河與洮河改善幅度最小,為1.43,其改善效果明顯程度由大到小排序:伊洛河、花園口以下干流區(qū)間、三門峽—小浪底區(qū)間、湟水、大夏河與洮河。以AIC值為例,改善效果最明顯的三級區(qū)是伊洛河,為37.04,大夏河與洮河改善幅度最小,為9.16,其改善效果明顯程度由大到小排序:伊洛河、三門峽—小浪底區(qū)間、花園口以下干流區(qū)間、湟水、大夏河與洮河。

    對于非汛期蒸發(fā)序列,以SBC值為例,改善效果最明顯的三級區(qū)是花園口以下干流區(qū)間,為43.68,伊洛河改善幅度最小,為5.88,其改善效果明顯程度由大到小排序:花園口以下干流區(qū)間、湟水、大夏河與洮河、三門峽—小浪底區(qū)間、伊洛河。以AIC值為例,改善效果最明顯的三級區(qū)是花園口以下干流區(qū)間,為47.55,伊洛河改善幅度最小,為13.91,其改善效果明顯程度由大到小排序:花園口以下干流區(qū)間、湟水、大夏河與洮河、三門峽—小浪底區(qū)間、伊洛河,其中伊洛河、三門峽—小浪底區(qū)間兩個三級區(qū)的SBC和AIC值非常接近。

    總體來說,在考慮統(tǒng)計參數(shù)非一致性特征情況下,5個三級區(qū)不同時間尺度的蒸發(fā)序列的擬合效果具有不同幅度的改善,說明這5個三級區(qū)的蒸發(fā)序列存在較為明顯的非一致性特征。以花園口以下干流區(qū)間的年、非汛期的蒸發(fā)序列的非一致性特征最為明顯,伊洛河汛期的蒸發(fā)序列次之,大夏河與洮河的年、汛期蒸發(fā)序列相對最差,湟水的汛期蒸發(fā)序列、伊洛河的非汛期蒸發(fā)序列、三門峽—小浪底區(qū)間的非汛期蒸發(fā)序列的擬合效果相對次之。

    5 結(jié) 論

    以黃河流域的5個水資源三級區(qū)為研究對象,采用Kendall秩次相關(guān)法、Mann-Kendall檢驗法和滑動T檢驗法,并引入GAMLSS模型,對5個三級區(qū)年、汛期、非汛期的蒸發(fā)序列進行非一致性研究。主要結(jié)論如下。

    (1)伊洛河、三門峽—小浪底區(qū)間、花園口以下干流區(qū)間的蒸發(fā)序列突變較為明顯,突變多發(fā)生在2000年左右;大夏河與洮河、湟水蒸發(fā)序列的突變存在一定的差異,大夏河與洮河突變發(fā)生在1990年左右,湟水蒸發(fā)序列突變發(fā)生在1975年左右。

    (2)5個三級區(qū)不同時間尺度的蒸發(fā)序列以Generalized Gamma為最優(yōu)分布的最多,Gumbel次之,Logistics最少。

    (3)5個三級區(qū)不同時間尺度上的蒸發(fā)序列除部分序列外,均與時間t建立了自由度為1~3的三次樣條函數(shù),且與時間t成負(fù)相關(guān);大部分蒸發(fā)數(shù)據(jù)均包含在置信水平為5%和95%的置信區(qū)間內(nèi),僅有極少數(shù)數(shù)據(jù)位于置信區(qū)間之外,表明不同時間尺度的蒸發(fā)序列的非一致性模型擬合效果較好。

    (4)在考慮統(tǒng)計參數(shù)非一致性特征情況下,5個三級區(qū)不同時間尺度的蒸發(fā)序列的擬合效果具有不同幅度的改善,說明這5個三級區(qū)的蒸發(fā)序列存在較為明顯的非一致性特征;以花園口以下干流區(qū)間的年、非汛期的蒸發(fā)序列的非一致性特征最為明顯。

    參考文獻:

    [1] 邱新法,劉昌明,曾燕.黃河流域近40年蒸發(fā)皿蒸發(fā)量的氣候變化特征[J].自然資源學(xué)報,2003,18(4):437-442.

    [2] 張冬冬.洪水頻率分析與預(yù)報中的不確定性問題研究[D].北京:中國水利水電科學(xué)研究院,2015:44-50.

    [3] XIONG L, DU T, XU C Y, et al. Non-Stationary Annual Maximum Flood Frequency Analysis Using the Norming Constants Method to Consider Non-Stationarity in the Annual Daily Flow Series[J]. Water Resources Management, 2015, 29(10):3615-3633.

    [4] 魯帆,肖偉華,嚴(yán)登華,等.非平穩(wěn)時間序列極值統(tǒng)計模型及其在氣候-水文變化研究中的應(yīng)用綜述[J].水利學(xué)報, 2017,48(4):379-389.

    [5] 嚴(yán)登華,袁喆,王浩,等.水文學(xué)確定性和不確定性方法及其集合研究進展[J].水利學(xué)報,2013,44(1):73-82.

    [6] 魯帆,嚴(yán)登華.基于廣義極值分布和Metropolis-Hastings抽樣算法的貝葉斯MCMC洪水頻率分析方法[J].水利學(xué)報,2013,44(8):942-949.

    [7] 魯帆,宋昕熠,朱奎,等.大渡河流域超閾值降雨樣本模擬及不確定性分析[J].人民長江,2016,47(7):23-27,33.

    [8] 熊立華,江聰,杜濤,等.變化環(huán)境下非一致性水文頻率分析研究綜述[J].水資源研究,2015,4(4):310-319.

    [9] 熊立華,江聰.考慮非一致性的渭河流域設(shè)計洪水過程線研究[J].水資源研究,2015,4(2):109-119.

    [10] 顧西輝,張強,王宗志.1951—2010年珠江流域洪水極值序列平穩(wěn)性特征研究[J].自然資源學(xué)報,2015,30(5):824-835.

    [11] 張冬冬,魯帆,周翔南,等.基于GAMLSS模型的大渡河流域極值降水非一致性分析[J].水利水電技術(shù),2016,47(5):12-15,20.

    [12] 熊斌,熊立華.基于基流退水過程的非一致性枯水頻率分析[J].水利學(xué)報,2016,47(7):873-883.

    [13] 江聰,熊立華.基于GAMLSS模型的宜昌站年徑流序列趨勢分析[J].地理學(xué)報,2012,67(11):1505-1514.

    [14] 魏鳳英.現(xiàn)代氣候統(tǒng)計診斷與預(yù)測技術(shù)[M].北京:氣象出版社,2007:63-71.

    [15] RIGBY B A, STASINOPOULOS D M. A Flexible Regression Approach Using GAMLSS in R[R]. London: London Metroditon University,2009:106-112.

    [16] MIKIS S, RIGBY R B, CALLIOPE A. Introductions on How to Use the GAMLSS Package in R[J].Chemical Communications, 2018(1): 67-69.

    [17] FILLIBEN J J. The Probability Plot Correlation Coefficient Test for Normality[J].Technometrics, 1975, 17(1):111-117.

    [18] 張彩玲.變化環(huán)境下黃河流域水文非一致性演變規(guī)律研究[D].鄭州:華北水利水電大學(xué),2018:16-20.

    [19] 賈文雄,何元慶,李宗省,等.近50年來河西走廊平原區(qū)氣候變化的區(qū)域特征及突變分析[J].地理科學(xué),2008,28(4):525-531.

    【責(zé)任編輯 張 帥】

    猜你喜歡
    黃河流域
    實施積極財政政策 推進黃河流域生態(tài)保護和高質(zhì)量發(fā)展
    山西財稅(2023年8期)2023-12-28 08:05:28
    生態(tài)環(huán)境部啟動新一年度黃河流域“清廢行動”
    黃河流域文學(xué)生態(tài)意識的地域透視
    “黃河流域生態(tài)保護和高度質(zhì)量發(fā)展”
    人民黃河(2022年9期)2022-09-08 13:04:38
    “黃河流域生態(tài)保護和高質(zhì)量發(fā)展”專題征稿
    人民黃河(2022年7期)2022-07-07 06:52:24
    “黃河流域生態(tài)保護和高質(zhì)量發(fā)展”專題征稿
    人民黃河(2022年4期)2022-04-07 09:04:04
    黃河流域甘肅段生態(tài)治理的法治保障分析
    在黃河流域生態(tài)保護和高質(zhì)量發(fā)展中展現(xiàn)陜西擔(dān)當(dāng)
    黃河流域燦爛的齊家文化——陶器
    收藏界(2019年2期)2019-10-12 08:26:10
    增強大局意識 提升黃河流域生態(tài)保護發(fā)展水平
    床上黄色一级片| 国产高清不卡午夜福利| 欧美潮喷喷水| 欧美潮喷喷水| 日本成人三级电影网站| 12—13女人毛片做爰片一| 人人妻人人澡欧美一区二区| 最新中文字幕久久久久| 男人舔女人下体高潮全视频| 婷婷色综合大香蕉| 搡女人真爽免费视频火全软件 | 少妇熟女aⅴ在线视频| 免费人成在线观看视频色| 国产探花在线观看一区二区| 亚洲在线自拍视频| 99热这里只有是精品50| 淫妇啪啪啪对白视频| 成人综合一区亚洲| 高清在线国产一区| 老司机午夜福利在线观看视频| 久久精品国产亚洲网站| 亚洲成a人片在线一区二区| 中文字幕精品亚洲无线码一区| 一级av片app| 十八禁国产超污无遮挡网站| 中文字幕久久专区| 国产探花在线观看一区二区| 看片在线看免费视频| 中文在线观看免费www的网站| 欧美日韩亚洲国产一区二区在线观看| 免费看日本二区| 国产精品久久久久久久久免| 一卡2卡三卡四卡精品乱码亚洲| 国产中年淑女户外野战色| 男插女下体视频免费在线播放| 在现免费观看毛片| 久久午夜亚洲精品久久| 精品欧美国产一区二区三| 国产麻豆成人av免费视频| 韩国av在线不卡| 国产淫片久久久久久久久| 欧美xxxx性猛交bbbb| 一级黄片播放器| 性欧美人与动物交配| 日本撒尿小便嘘嘘汇集6| 狠狠狠狠99中文字幕| 我的女老师完整版在线观看| 亚洲精品日韩av片在线观看| 国产精品福利在线免费观看| 美女免费视频网站| 人人妻人人澡欧美一区二区| 嫩草影视91久久| 国产精品综合久久久久久久免费| 身体一侧抽搐| 欧美又色又爽又黄视频| 久久国产乱子免费精品| 男女下面进入的视频免费午夜| 自拍偷自拍亚洲精品老妇| 少妇人妻一区二区三区视频| www日本黄色视频网| 他把我摸到了高潮在线观看| 99久久精品热视频| 亚洲av.av天堂| 99热网站在线观看| 草草在线视频免费看| 国产高清激情床上av| 深夜精品福利| 亚洲av免费在线观看| 国产精品人妻久久久影院| a级毛片a级免费在线| 国产精品乱码一区二三区的特点| 日韩在线高清观看一区二区三区 | 内射极品少妇av片p| 亚洲无线观看免费| 国产精品野战在线观看| 波多野结衣巨乳人妻| 18+在线观看网站| 一个人看视频在线观看www免费| 欧美绝顶高潮抽搐喷水| 天天一区二区日本电影三级| 日本熟妇午夜| 国产精品98久久久久久宅男小说| 亚洲专区国产一区二区| 久久香蕉精品热| av福利片在线观看| 小说图片视频综合网站| 女同久久另类99精品国产91| 女同久久另类99精品国产91| 国产精品一区二区性色av| 一区二区三区激情视频| 91久久精品电影网| 五月伊人婷婷丁香| 一级a爱片免费观看的视频| 男女边吃奶边做爰视频| 男女边吃奶边做爰视频| 性插视频无遮挡在线免费观看| 国产真实乱freesex| 午夜久久久久精精品| 特级一级黄色大片| 1000部很黄的大片| 国产精品久久久久久精品电影| 亚洲精品粉嫩美女一区| а√天堂www在线а√下载| 国产亚洲精品久久久com| 99在线视频只有这里精品首页| 国产麻豆成人av免费视频| 亚洲一区高清亚洲精品| 美女高潮的动态| 日韩大尺度精品在线看网址| 在线播放无遮挡| 又紧又爽又黄一区二区| 久久热精品热| 国产免费av片在线观看野外av| 国内揄拍国产精品人妻在线| 国产精品国产三级国产av玫瑰| 色吧在线观看| 久久亚洲真实| 欧美丝袜亚洲另类 | 国产黄片美女视频| 国产午夜福利久久久久久| 97人妻精品一区二区三区麻豆| 搡老岳熟女国产| 日韩欧美三级三区| 国产黄色小视频在线观看| 日韩高清综合在线| 一级黄片播放器| 1024手机看黄色片| 国产精华一区二区三区| 一区二区三区免费毛片| 国产美女午夜福利| 中文字幕精品亚洲无线码一区| 久久精品国产亚洲网站| 欧美中文日本在线观看视频| 久久久久精品国产欧美久久久| 国内久久婷婷六月综合欲色啪| 级片在线观看| 天堂影院成人在线观看| 成人精品一区二区免费| 可以在线观看毛片的网站| 欧美精品啪啪一区二区三区| 国产探花极品一区二区| 国产午夜福利久久久久久| 精品人妻视频免费看| 成人特级av手机在线观看| 一个人免费在线观看电影| 在线观看午夜福利视频| 成年人黄色毛片网站| 美女 人体艺术 gogo| 女人被狂操c到高潮| 亚洲不卡免费看| 最近最新中文字幕大全电影3| 免费av观看视频| av在线亚洲专区| 伦精品一区二区三区| 黄色配什么色好看| 日韩国内少妇激情av| 欧美日韩瑟瑟在线播放| 婷婷精品国产亚洲av| 婷婷亚洲欧美| h日本视频在线播放| 亚洲色图av天堂| 精品福利观看| 12—13女人毛片做爰片一| 精品人妻1区二区| 黄色女人牲交| 国产精品自产拍在线观看55亚洲| 国产伦人伦偷精品视频| 夜夜爽天天搞| 人妻少妇偷人精品九色| 国产精品久久电影中文字幕| 啦啦啦观看免费观看视频高清| 国产亚洲精品av在线| .国产精品久久| 黄色视频,在线免费观看| 99热网站在线观看| 长腿黑丝高跟| 一夜夜www| 桃红色精品国产亚洲av| 蜜桃亚洲精品一区二区三区| 日日摸夜夜添夜夜添小说| 99精品在免费线老司机午夜| 天天一区二区日本电影三级| 亚洲欧美日韩东京热| 麻豆成人av在线观看| 国内少妇人妻偷人精品xxx网站| 99热网站在线观看| 午夜福利在线观看吧| 两人在一起打扑克的视频| 91麻豆av在线| 国产精品电影一区二区三区| 99精品在免费线老司机午夜| 熟妇人妻久久中文字幕3abv| 女同久久另类99精品国产91| 久久久久久九九精品二区国产| 午夜福利高清视频| 婷婷色综合大香蕉| 乱人视频在线观看| 嫩草影院入口| 春色校园在线视频观看| 人妻少妇偷人精品九色| 美女 人体艺术 gogo| 婷婷丁香在线五月| 中文字幕人妻熟人妻熟丝袜美| 成人欧美大片| 国产精品乱码一区二三区的特点| 精品不卡国产一区二区三区| 国产高清视频在线播放一区| 国产精品福利在线免费观看| 婷婷六月久久综合丁香| 免费高清视频大片| or卡值多少钱| 嫩草影院入口| aaaaa片日本免费| 男人狂女人下面高潮的视频| 99热这里只有精品一区| 国产 一区 欧美 日韩| 人妻丰满熟妇av一区二区三区| 国产伦精品一区二区三区四那| 国产伦一二天堂av在线观看| 永久网站在线| 国产精品福利在线免费观看| 国产精品一区二区三区四区久久| 亚洲av熟女| 99精品久久久久人妻精品| 亚洲国产精品久久男人天堂| a级一级毛片免费在线观看| 欧美区成人在线视频| 久久久久久九九精品二区国产| 色在线成人网| 日本精品一区二区三区蜜桃| 欧美性猛交╳xxx乱大交人| 国产欧美日韩一区二区精品| 亚洲精品乱码久久久v下载方式| 人妻丰满熟妇av一区二区三区| 国内精品一区二区在线观看| 在线观看av片永久免费下载| 亚洲精品日韩av片在线观看| 国产高清视频在线观看网站| 欧美成人性av电影在线观看| 中文字幕精品亚洲无线码一区| 国内久久婷婷六月综合欲色啪| 精品久久久久久久久久免费视频| 免费一级毛片在线播放高清视频| 日本黄大片高清| 2021天堂中文幕一二区在线观| 精品一区二区免费观看| 成人特级av手机在线观看| 村上凉子中文字幕在线| 亚洲18禁久久av| 免费人成视频x8x8入口观看| 99热6这里只有精品| 欧洲精品卡2卡3卡4卡5卡区| a级一级毛片免费在线观看| 免费无遮挡裸体视频| 亚洲成人久久性| 色噜噜av男人的天堂激情| 一夜夜www| 给我免费播放毛片高清在线观看| www.www免费av| 极品教师在线视频| 美女被艹到高潮喷水动态| 精品久久久噜噜| 麻豆成人av在线观看| 五月伊人婷婷丁香| 91av网一区二区| 女生性感内裤真人,穿戴方法视频| 国模一区二区三区四区视频| 亚洲欧美精品综合久久99| 成人一区二区视频在线观看| 看十八女毛片水多多多| 丰满的人妻完整版| 色精品久久人妻99蜜桃| 国产 一区 欧美 日韩| 一个人看视频在线观看www免费| 亚洲成人久久性| 久久99热这里只有精品18| 欧美一区二区国产精品久久精品| 午夜福利在线观看吧| 成人二区视频| 中文字幕av在线有码专区| 亚洲中文字幕一区二区三区有码在线看| 欧美日韩精品成人综合77777| 国产一区二区在线av高清观看| www.www免费av| 欧美日韩亚洲国产一区二区在线观看| 在线观看免费视频日本深夜| 欧美不卡视频在线免费观看| 日本免费一区二区三区高清不卡| 国产主播在线观看一区二区| 欧美性猛交黑人性爽| 日韩高清综合在线| 中文字幕人妻熟人妻熟丝袜美| 一个人看的www免费观看视频| 久久香蕉精品热| 国产免费av片在线观看野外av| 成人午夜高清在线视频| 国内精品久久久久久久电影| 久久久精品大字幕| 午夜福利在线观看免费完整高清在 | 天天一区二区日本电影三级| 久久精品国产自在天天线| 99在线人妻在线中文字幕| 亚洲aⅴ乱码一区二区在线播放| 亚洲av免费高清在线观看| 国产大屁股一区二区在线视频| 亚洲国产精品合色在线| 精品人妻1区二区| 亚洲精品色激情综合| 欧美zozozo另类| 久久精品国产亚洲网站| 成人欧美大片| 在线观看66精品国产| 亚洲精品色激情综合| 亚洲欧美清纯卡通| 成人一区二区视频在线观看| 少妇人妻一区二区三区视频| 嫁个100分男人电影在线观看| 少妇的逼水好多| 人人妻人人澡欧美一区二区| 国产伦人伦偷精品视频| 亚洲av电影不卡..在线观看| a级一级毛片免费在线观看| 国语自产精品视频在线第100页| 91av网一区二区| 亚洲五月天丁香| 99久久成人亚洲精品观看| 日韩国内少妇激情av| 色吧在线观看| 男女视频在线观看网站免费| 中文字幕免费在线视频6| av福利片在线观看| 国产精品,欧美在线| 18禁黄网站禁片午夜丰满| 欧美一级a爱片免费观看看| 日韩精品中文字幕看吧| av国产免费在线观看| 琪琪午夜伦伦电影理论片6080| 男插女下体视频免费在线播放| 国产av一区在线观看免费| 免费无遮挡裸体视频| 久久这里只有精品中国| 毛片一级片免费看久久久久 | 国产在视频线在精品| 伦理电影大哥的女人| 国内揄拍国产精品人妻在线| 亚洲黑人精品在线| 欧美日韩亚洲国产一区二区在线观看| 亚洲欧美日韩东京热| av国产免费在线观看| 制服丝袜大香蕉在线| 国产视频内射| 亚洲精华国产精华液的使用体验 | 久久午夜亚洲精品久久| 在线播放国产精品三级| 国产白丝娇喘喷水9色精品| eeuss影院久久| 丝袜美腿在线中文| 一卡2卡三卡四卡精品乱码亚洲| 国产精品国产高清国产av| 在线观看一区二区三区| 免费人成视频x8x8入口观看| 免费看av在线观看网站| 91狼人影院| 美女大奶头视频| 嫩草影院新地址| 老司机福利观看| h日本视频在线播放| 亚洲综合色惰| 一边摸一边抽搐一进一小说| 蜜桃久久精品国产亚洲av| 精品福利观看| 色吧在线观看| 国产淫片久久久久久久久| 亚洲精华国产精华液的使用体验 | 久久99热6这里只有精品| 欧美人与善性xxx| 国产av在哪里看| 床上黄色一级片| 亚洲一区高清亚洲精品| 亚洲色图av天堂| 亚洲精品色激情综合| 简卡轻食公司| 在线观看舔阴道视频| 男女下面进入的视频免费午夜| 性色avwww在线观看| 一本精品99久久精品77| 成人av一区二区三区在线看| 国产一区二区亚洲精品在线观看| 97超级碰碰碰精品色视频在线观看| 一个人看视频在线观看www免费| 亚洲美女视频黄频| 午夜精品久久久久久毛片777| 99精品在免费线老司机午夜| 亚洲精品在线观看二区| 国产精品女同一区二区软件 | 91av网一区二区| 人妻少妇偷人精品九色| 白带黄色成豆腐渣| 淫秽高清视频在线观看| 欧美3d第一页| 欧美日韩综合久久久久久 | 久久午夜福利片| 国产大屁股一区二区在线视频| 成人毛片a级毛片在线播放| 亚洲精品亚洲一区二区| 久久久久久九九精品二区国产| 亚洲18禁久久av| 国产伦人伦偷精品视频| 97超级碰碰碰精品色视频在线观看| 久久精品国产亚洲av涩爱 | 99热6这里只有精品| 亚洲国产日韩欧美精品在线观看| 精品一区二区免费观看| bbb黄色大片| 久久精品国产亚洲av涩爱 | 日日摸夜夜添夜夜添小说| 全区人妻精品视频| 亚洲精品国产成人久久av| 亚洲图色成人| 免费看a级黄色片| 日韩 亚洲 欧美在线| 国产精品一区二区三区四区免费观看 | 久久精品国产亚洲av天美| 熟女电影av网| 色综合站精品国产| 亚洲第一电影网av| 男人舔奶头视频| 一本久久中文字幕| netflix在线观看网站| 久久精品国产亚洲av香蕉五月| 国产淫片久久久久久久久| 国产成人一区二区在线| 在线观看一区二区三区| 麻豆成人午夜福利视频| 欧美绝顶高潮抽搐喷水| av中文乱码字幕在线| 亚洲精华国产精华精| 淫秽高清视频在线观看| 可以在线观看的亚洲视频| 午夜视频国产福利| 欧美日韩精品成人综合77777| 日韩亚洲欧美综合| 国产精品亚洲美女久久久| 美女黄网站色视频| 麻豆av噜噜一区二区三区| 国产精品亚洲一级av第二区| 国产精品爽爽va在线观看网站| 内地一区二区视频在线| 91午夜精品亚洲一区二区三区 | 久久人妻av系列| 国产av在哪里看| 嫁个100分男人电影在线观看| 亚洲,欧美,日韩| 十八禁网站免费在线| 非洲黑人性xxxx精品又粗又长| 尾随美女入室| 亚洲一级一片aⅴ在线观看| 欧美色欧美亚洲另类二区| 久久香蕉精品热| av在线天堂中文字幕| av黄色大香蕉| 国产男靠女视频免费网站| 性色avwww在线观看| 欧美日本视频| 久久人人爽人人爽人人片va| 99久久成人亚洲精品观看| 亚洲精品在线观看二区| 性色avwww在线观看| 搡老熟女国产l中国老女人| 国产91精品成人一区二区三区| 在线看三级毛片| av在线蜜桃| 成人无遮挡网站| 看十八女毛片水多多多| 欧美日韩乱码在线| 午夜精品一区二区三区免费看| 99久久精品热视频| 国产单亲对白刺激| 午夜福利18| 精品无人区乱码1区二区| 99久久久亚洲精品蜜臀av| 69av精品久久久久久| 国产激情偷乱视频一区二区| avwww免费| 中文字幕精品亚洲无线码一区| 一级黄色大片毛片| 免费大片18禁| 欧美一级a爱片免费观看看| 搡老岳熟女国产| av女优亚洲男人天堂| 婷婷亚洲欧美| 九九热线精品视视频播放| 亚洲第一电影网av| 亚洲乱码一区二区免费版| 午夜福利视频1000在线观看| 一个人看视频在线观看www免费| 欧美精品国产亚洲| 亚洲男人的天堂狠狠| 97碰自拍视频| 国产精华一区二区三区| 看免费成人av毛片| 男女边吃奶边做爰视频| 天堂影院成人在线观看| 国产午夜精品久久久久久一区二区三区 | 亚洲不卡免费看| 97碰自拍视频| 亚洲精品久久国产高清桃花| 亚洲天堂国产精品一区在线| 久久久久久九九精品二区国产| 欧美日韩乱码在线| 久久草成人影院| 国内精品美女久久久久久| 天天躁日日操中文字幕| av在线亚洲专区| 床上黄色一级片| 99久久成人亚洲精品观看| 国产 一区精品| 美女 人体艺术 gogo| 日本黄色视频三级网站网址| 欧美激情在线99| 亚洲性久久影院| 99在线人妻在线中文字幕| 尤物成人国产欧美一区二区三区| 精品人妻视频免费看| 国产成人av教育| 亚洲一区二区三区色噜噜| 色精品久久人妻99蜜桃| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品日韩av在线免费观看| 国产色爽女视频免费观看| 日日撸夜夜添| 免费黄网站久久成人精品| 精品久久久久久,| 国产中年淑女户外野战色| 不卡一级毛片| 在线观看午夜福利视频| 老熟妇乱子伦视频在线观看| 最近中文字幕高清免费大全6 | 人妻夜夜爽99麻豆av| 18禁裸乳无遮挡免费网站照片| 国产免费男女视频| 免费看美女性在线毛片视频| 欧美日本视频| 亚洲av日韩精品久久久久久密| 99精品在免费线老司机午夜| 成人午夜高清在线视频| 美女被艹到高潮喷水动态| 亚洲国产欧洲综合997久久,| 国产综合懂色| 99久久精品热视频| 成熟少妇高潮喷水视频| 国内少妇人妻偷人精品xxx网站| 午夜精品在线福利| 国内精品美女久久久久久| 99精品久久久久人妻精品| 国产高潮美女av| 国产欧美日韩一区二区精品| 三级国产精品欧美在线观看| 美女cb高潮喷水在线观看| 我的老师免费观看完整版| 亚洲av免费高清在线观看| 国产乱人伦免费视频| 亚洲黑人精品在线| 国产乱人伦免费视频| 亚洲七黄色美女视频| 亚洲最大成人中文| 一区二区三区四区激情视频 | 国产一区二区三区在线臀色熟女| 可以在线观看毛片的网站| 国产av不卡久久| 成人美女网站在线观看视频| 免费不卡的大黄色大毛片视频在线观看 | 熟女电影av网| 日韩欧美精品免费久久| 网址你懂的国产日韩在线| 男女之事视频高清在线观看| 老司机深夜福利视频在线观看| 俄罗斯特黄特色一大片| 亚洲,欧美,日韩| 欧美成人性av电影在线观看| 国产精品一区二区三区四区免费观看 | 国产爱豆传媒在线观看| 久99久视频精品免费| 欧美又色又爽又黄视频| 成年人黄色毛片网站| 成年女人看的毛片在线观看| 久久久久九九精品影院| 女同久久另类99精品国产91| 男女视频在线观看网站免费| 看十八女毛片水多多多| 99国产精品一区二区蜜桃av| 麻豆一二三区av精品| 色吧在线观看| 99热精品在线国产| 国产一区二区激情短视频| 精品久久久久久久人妻蜜臀av| 国产精品自产拍在线观看55亚洲| 99精品在免费线老司机午夜| 国产精品一及| 久久久久久久亚洲中文字幕| 99热只有精品国产| 久久精品国产亚洲av天美| 色在线成人网| 蜜桃亚洲精品一区二区三区| 午夜福利欧美成人| 免费黄网站久久成人精品| 九九久久精品国产亚洲av麻豆| 黄色丝袜av网址大全| www日本黄色视频网| 99热这里只有精品一区| 黄色女人牲交| 亚洲人成网站高清观看| 天堂√8在线中文| 少妇的逼水好多| 18禁在线播放成人免费|