余　超，晁麗君，張夢(mèng)婕，朱家昌，孫靖智

(1.河海大學(xué) 理學(xué)院， 南京　210098； 2.河海大學(xué) 文天學(xué)院， 安徽 馬鞍山　243031；

3.河海大學(xué) 水文水資源學(xué)院， 南京　210098； 4.天津市水務(wù)局,天津　300074；

5.遼寧省水利廳， 沈陽(yáng)　110003)

?

基于回歸模型的方法分析半干旱流域的降雨徑流關(guān)系

余超1，晁麗君2，張夢(mèng)婕3，朱家昌4，孫靖智5

(1.河海大學(xué) 理學(xué)院， 南京210098； 2.河海大學(xué) 文天學(xué)院， 安徽 馬鞍山243031；

3.河海大學(xué) 水文水資源學(xué)院， 南京210098； 4.天津市水務(wù)局,天津300074；

5.遼寧省水利廳， 沈陽(yáng)110003)

摘要：以半干旱流域大閣流域?yàn)槔?，通過回歸模型改進(jìn)半干旱流域降雨徑流關(guān)系，用3種方法建立回歸方程，然后對(duì)3種方法的結(jié)果進(jìn)行評(píng)定。比較3種方法的合格率、許可誤差和確定性系數(shù)。結(jié)果表明：采用方法3在大閣流域模擬降雨徑流達(dá)到了預(yù)報(bào)精度，很好地描述了降雨徑流關(guān)系，從而為改進(jìn)半干旱流域的降雨徑流關(guān)系提供新思路。

關(guān)鍵詞：回歸模型； 降雨徑流；相關(guān)性；半干旱流域；評(píng)定方法；多因子；可分離變量

基于回歸模型方法[1]分析降雨徑流關(guān)系[2]是將數(shù)據(jù)集在空間或時(shí)間上外推到其他情形的基本統(tǒng)計(jì)計(jì)算機(jī)技術(shù)。數(shù)據(jù)用于以因變量和一個(gè)或多個(gè)自變量或解釋變量的收集整理，統(tǒng)計(jì)計(jì)算用于推導(dǎo)建立一個(gè)因變量與自變量的相關(guān)關(guān)系的方程。

一次降雨的產(chǎn)流量和很多因素有關(guān)，也就是R=f(P,Pa,i,季節(jié)，…)[3]的函數(shù)關(guān)系。本方法以流域產(chǎn)流的物理機(jī)理為基礎(chǔ)，以主要影響因素為參變量，建立各因素與產(chǎn)流量R之間的關(guān)系。這里的各因素包括降雨量P、前期雨量指數(shù)Pa、降雨強(qiáng)度i、季節(jié)等等。

1回歸方程的描述

多元線性回歸模型[4]是描述因變量y依n項(xiàng)自變量x1,x2,x3,…,xn和誤差項(xiàng)ε取值的方程。

1.1線性回歸模型

設(shè)線性回歸模型的一般形式為

(1)

若有m組由因變量y和n項(xiàng)自變量x1,x2,x3,…,xn數(shù)據(jù)構(gòu)成的樣本，依式(1)可組成線性回歸模型，可用矩陣表示為

(2)

其中：

利用最小二乘法原理[5]求方程組(2)，從而建立線性回歸方程。

1.2非線性回歸模型

當(dāng)擬合曲線是非線性時(shí)，往往按下述方法求其回歸方程：

假設(shè)有n個(gè)影響因子：

(3)

首先假設(shè)y是關(guān)于各因子的可分離變量函數(shù)[6]，即

(4)

在自變量x1,x2,x3,…,xn中找出對(duì)倚變量y影響最大的主要因子。假設(shè)x1對(duì)y影響最大，建立數(shù)學(xué)模型：

(5)

以此確定主要因素與y的函數(shù)關(guān)系。當(dāng)取p=2時(shí)，即為最小二乘法。

確定y1=f1(x1)及殘差Δy1=y-y1。由于y是x1,x2,x3,…,xn的函數(shù)，自y中分離x1的函數(shù)y1=f1(x1)后，殘差Δy1應(yīng)近似為x2,x3,…,xn的函數(shù)。再分析Δy1與x2,x3,…,xn的相關(guān)關(guān)系。比較各因素與殘差的相關(guān)性，往往會(huì)應(yīng)用3種方法。

(6)

方法2利用所建立的數(shù)學(xué)模型(5)分別求出Δy1與x2,x3,…,xn的函數(shù)，即Δy1=f(x2),Δy1=f(x3),…,Δy1=f(xn)，然后建立超定方程：

(7)

用最小二乘法求出a2,a3,…,an，由此得出

(8)

方法3利用所建立的數(shù)學(xué)模型(5)求出Δy1與x2,x3,…,xn的函數(shù)：

(9)

2相關(guān)分析

在研究?jī)蓚€(gè)水文序列之間的關(guān)系時(shí)，可采用相關(guān)分析法[7]。例如對(duì)于一次洪水的徑流深和流域一次洪水的降雨量，可以通過互相關(guān)來表示同一次洪水的徑流量與降雨量之間的關(guān)系。

水文序列Xt和Yt的相關(guān)系數(shù)ρk(X,Y)定義為

(10)

式中C(X,Y)為序列Xt和Yt的互協(xié)方差。通常應(yīng)用式(11)由樣本序列(Xt,Yt;t=1,2,…,n)估計(jì)互相關(guān)系數(shù)。

(11)

3降雨徑流關(guān)系中各要素的計(jì)算原理

3.1流域平均次降雨量計(jì)算

流域平均次降雨量又稱面雨量，是指一次降雨過程中整個(gè)流域面上的平均累計(jì)降雨量。推求流域平均次降雨量最常用的方法有算術(shù)平均法和垂直平分法(又叫做泰森多邊形法)，也有用繪制等雨量線圖來推求的。

(12)

其中：ωi為流域內(nèi)各雨量站權(quán)重(無量綱)；Pi為流域內(nèi)各雨量站的時(shí)段降雨量(mm)；P為時(shí)段流域平均降雨量(mm)；n為流域內(nèi)雨量站的個(gè)數(shù)(無量綱)。

3.2前期影響雨量的計(jì)算

前期影響雨量Pa為衡量流域濕潤(rùn)程度的指標(biāo)[8]，反映流域蓄水量的大小。Pa越小表示流域的干旱程度越大；Pa越大表示流域的濕潤(rùn)程度越大。流域前期影響雨量(Pa值)采用遞推公式計(jì)算，并以流域土壤最大損失量(Im值)進(jìn)行控制。Im為流域最大損失量，范圍一般在80～120 mm。本次研究中Im統(tǒng)一取為120 mm。單站前期影響雨量Pa的計(jì)算式為[9]

(13)

其中：Pa,t，Pa,t+1分別為第t天和第t+1天開始時(shí)刻的前期影響雨量(mm)，Pa,t+1≤Im；Pt為第t天該站的日降雨量(mm)；Ka為流域土壤含水量消退系數(shù)(無量綱)，取值范圍一般為0.8～0.95。

3.3雨強(qiáng)的計(jì)算

根據(jù)前面計(jì)算的流域平均次雨量P可以推出雨強(qiáng)i：

(14)

其中：P為時(shí)段流域平均降雨量；t為本次洪水的歷時(shí)；i為雨強(qiáng)。

3.4地表徑流量的計(jì)算

產(chǎn)流是流域內(nèi)各種徑流成分的生成過程，也就是流域下墊面對(duì)降雨的再分配過程。當(dāng)降雨強(qiáng)度小于下滲強(qiáng)度時(shí)，降落在地面的雨水將全部滲入土壤；當(dāng)降雨強(qiáng)度大于下滲能力時(shí)，雨水除按下滲能力入滲外，超出下滲能力的部分便形成地面徑流。滲入土壤中的水量除被土壤蒸發(fā)和植物吸收散發(fā)而損耗掉外，余下的補(bǔ)充了土壤含水量。當(dāng)土壤包氣帶的水量超過田間持水量時(shí)，則產(chǎn)生壤中徑流，同時(shí)通過淺層地下徑流和深層地下徑流補(bǔ)給河流。因此，一次降水的產(chǎn)流過程是由不同水源組成的徑流過程。

3.4.1流量過程的分割

由于地面徑流和地下徑流有不同的匯流特性[10]，根據(jù)實(shí)測(cè)洪水流量過程線，采用斜線分割法對(duì)地面徑流和地下徑流進(jìn)行劃分。首先，從整個(gè)流量過程線中分割出前、后期徑流和深層基流，并獲得本次降雨所形成的流量過程，然后進(jìn)行地表、地下流量過程的分割，這就是所謂的基流分割[11]。該方法的基本思路：先尋找洪水過程的直接徑流終止點(diǎn)，然后用斜線連接漲點(diǎn)與終止點(diǎn)，則斜線上部為地面徑流，下部為地下徑流。

對(duì)于半干旱、干旱地區(qū)，地面徑流終止點(diǎn)可用目估法確定。對(duì)于濕潤(rùn)、半濕潤(rùn)地區(qū)，地面徑流終止點(diǎn)可通過經(jīng)驗(yàn)法確定。該方法是通過對(duì)實(shí)測(cè)資料的分析確定洪峰時(shí)間到地面徑流終止點(diǎn)的時(shí)距N[12]，并認(rèn)為同一流域的N值是常數(shù)。求取N的經(jīng)驗(yàn)公式為

(15)

其中：F為流域面積(km2)；N為洪峰時(shí)間到地面徑流終止點(diǎn)的時(shí)距(h)。

3.4.2地表徑流量的計(jì)算

在分割出一次降雨所形成的地面徑流過程線后，即可根據(jù)起漲點(diǎn)、終止點(diǎn)及洪水過程線包圍的面積計(jì)算次洪徑流深。計(jì)算公式為

(16)

其中：Rs為流域地面徑流深(mm)；Qi為流域時(shí)段流量(m3/s)；Δti為計(jì)算時(shí)段(h)，常取1 h；F為流域面積(km2)；3.6為單位換算系數(shù)(無量綱)。

4精度評(píng)定

徑流深預(yù)報(bào)以實(shí)測(cè)值的20%作為許可誤差，當(dāng)一次預(yù)報(bào)的誤差小于許可誤差時(shí)為合格預(yù)報(bào)。

5流域概況

大閣流域位于密云水庫(kù)流域。密云水庫(kù)建于1960年，蓄水以來，為京、津、冀部分地區(qū)工農(nóng)業(yè)與城市生活用水發(fā)揮了巨大作用。潮白河密云水庫(kù)以上流域位于東經(jīng)115°30′～117°30′，北緯 40°20′～41°30′。東部和東北部與灤河流域?yàn)榻纾泵媾c內(nèi)陸河閃電河為鄰，西與永定河支流洋河流域相鄰。密云水庫(kù)流域面積15 788km2，占潮白河流域面積18 000km2的88%。本流域除水庫(kù)庫(kù)區(qū)附近地勢(shì)稍低，其高程一般在海拔500m以下外，其他地方地勢(shì)均較高，具有山高、坡陡、溝深、流急的特點(diǎn)，河道比降較大。東支潮河發(fā)源于河北省承德地區(qū)，河源至水庫(kù)全長(zhǎng)220km，河道平均比降為1.87%。西支白河發(fā)源于河北省張家口地區(qū)，河源至水庫(kù)全長(zhǎng)248km，河道平均比降為4.87%。流域內(nèi)基本屬于土石山區(qū)，一般土層較薄，植被較好，但流域內(nèi)各地差異較大，不少地方巖石裸露，裂隙較大。

密云水庫(kù)流域面積較大，本次研究采用潮河的大閣流域。大閣位于東經(jīng)116°41′，北緯41°11′，流域面積為1 850km2。該流域?qū)儆谥芯暥却箨懶约撅L(fēng)氣候，據(jù)1951—1981年資料統(tǒng)計(jì)，流域平均降雨量為525mm，暴雨分布以庫(kù)區(qū)為中心，分別沿潮白河向東北和西北方向漸減。潮河上游多年平均降雨量約為500mm。降雨量主要集中在6～9 月份，尤其集中在7～8月份。

6算例

選取大閣流域?yàn)槔?，倚變量為流域地面徑流深Rs，自變量分別是流域平均次降雨量P、前期影響雨量Pa、雨強(qiáng)i，共44組資料。本次算例采用泰森多邊形法計(jì)算流域平均次降雨量P[13-18]。大閣流域雨量站及權(quán)重和位置如表1所示。

表1　大閣流域雨量站及其權(quán)重、位置

6.1計(jì)算相關(guān)性

采用式(10)計(jì)算出倚變量流域地面徑流深Rs與自變量流域平均次降雨量P、前期影響雨量Pa、雨強(qiáng)i各自的相關(guān)系數(shù)，結(jié)果如表2所示。

表2　相關(guān)系數(shù)

通過比較可以看出：倚變量流域地面徑流深Rs與流域平均次降雨量P相關(guān)系數(shù)最大，其次是雨強(qiáng)i，最小的是前期影響雨量Pa。這里采用的是Pearson 相關(guān)性和顯著性雙側(cè)檢驗(yàn)。

6.2計(jì)算結(jié)果

本算例分別采用方法1～3進(jìn)行計(jì)算。選取大閣流域44場(chǎng)洪水中的前30場(chǎng)洪水用于模型的參數(shù)率定，后14場(chǎng)用于預(yù)報(bào)模型驗(yàn)證。由相關(guān)系數(shù)可知：倚變量流域地面徑流深Rs與流域平均次降雨量P的相關(guān)系數(shù)最大。利用式(5)擬定Rs=f(P)，計(jì)算結(jié)果如表3所示，式(5)中取p=2時(shí)即為最小二乘法。

表3　Rs= f(P)計(jì)算結(jié)果

序號(hào)Rs實(shí)測(cè)PPaif(P)Δy1162.4217.0974.470.280.002.42170.677.2245.280.120.000.67183.7336.8825.660.613.80-0.07196.2140.1276.080.396.140.07200.9613.5726.770.130.000.96210.9314.5127.150.160.000.93221.7817.0352.410.200.001.78233.2531.6982.680.380.053.20242.1439.2576.600.382.19-0.05250.7430.4641.580.300.070.67261.2313.0945.490.210.001.23270.7515.6047.120.220.000.75280.9512.9351.840.180.000.95297.1066.0079.150.637.12-0.02303.5249.0544.080.653.57-0.04

采用方法1計(jì)算出Δy1后用式(5)擬定Δy1=f(i)，而后計(jì)算Δy2，即Δy2=f(P)+f(i)-Rs實(shí)測(cè)。用式(5)擬定Δy2=f(Pa)，最后得到Rs預(yù)報(bào)值=f(P)+f(i)+f(Pa)，計(jì)算結(jié)果如表4所示。

表4　 方法1計(jì)算結(jié)果

采用方法3計(jì)算出Δy1后用式(5)擬定Δy1=g(i,Pa)，得到Rs預(yù)報(bào)值=f(P)+g(i,Pa)，方法2、3的計(jì)算結(jié)果如表5所示。

表5　方法2、3計(jì)算結(jié)果

6.3誤差評(píng)定

根據(jù)徑流深預(yù)報(bào)以實(shí)測(cè)值的20%作為許可誤差,且許可誤差小于20 mm時(shí)點(diǎn)據(jù)評(píng)定為合格。對(duì)大閣流域的14場(chǎng)洪水進(jìn)行方案評(píng)定，結(jié)果如下：方法1得到的合格率QR為71.4%，相對(duì)平均誤差為18.9%，確定性系數(shù)DC為0.997；方法2得到的合格率QR為42.9%，相對(duì)平均誤差為26%，確定性系數(shù)DC為0.995；方法3得到的合格率QR為78.6%，相對(duì)平均誤差為18.3%，確定性系數(shù)DC為0.995。誤差計(jì)算結(jié)果如表6所示。

從評(píng)定結(jié)果可以看出：3種方法的確定性系數(shù)都很高，達(dá)到評(píng)定要求，但是只有方法3合格率達(dá)到乙級(jí)精度要求。分析數(shù)據(jù)可以看出：產(chǎn)生這種現(xiàn)象的原因是大閣流域的徑流深相對(duì)較小，但是其中會(huì)出現(xiàn)一些大數(shù)，這就導(dǎo)致確定性系數(shù)普遍很高。

根據(jù)誤差評(píng)定出的結(jié)果，大閣流域應(yīng)選取方法3進(jìn)行徑流深的預(yù)報(bào)。通過與文獻(xiàn)[9]中應(yīng)用的降雨徑流相關(guān)圖法進(jìn)行比較，運(yùn)用降雨徑流相關(guān)圖法的合格率QR只有34%，相對(duì)平均誤差為36%；運(yùn)用回歸模型的合格率QR為78.6%，相對(duì)平均誤差為18.3%，大大提高了精度。

表6　誤差計(jì)算結(jié)果

7結(jié)論

1) 基于回歸模型的分析方法可用于估算模型參數(shù)或無觀察資料地區(qū)的徑流深。

2) 通過與文獻(xiàn)[9]在大閣流域的比較可以看出：預(yù)報(bào)精度得到明顯提高，說明基于回歸模型的分析方法可以彌補(bǔ)降雨徑流相關(guān)圖在半干旱地區(qū)應(yīng)用的不足。

3) 模擬的關(guān)鍵點(diǎn)是在最終確定方程時(shí)如何選取自變量。

4) 回歸的預(yù)報(bào)值取決于模型與原始數(shù)據(jù)的擬合情況。

5) 所有預(yù)報(bào)成果都應(yīng)當(dāng)計(jì)算誤差聯(lián)系起來，這種估算標(biāo)準(zhǔn)誤差可表述預(yù)報(bào)結(jié)果的不確定性。

6) 預(yù)報(bào)成果所依據(jù)數(shù)據(jù)的代表性好，其估算的標(biāo)準(zhǔn)誤差就??；代表性差，標(biāo)準(zhǔn)誤差就大。

參考文獻(xiàn)：

[1]袁慰平.計(jì)算方法與實(shí)習(xí)[M].南京:東南大學(xué)出版社,1991.

[2]貝文.降雨—徑流模擬[M].馬俊，譯.北京:中國(guó)水利出版社,2006.

[3]芮孝芳.產(chǎn)流模式的發(fā)現(xiàn)與發(fā)展[J].水利水電科技進(jìn)展,2013,33(1):2-5.

[4]張潤(rùn)楚.多元統(tǒng)計(jì)分析[M].北京:科學(xué)出版社,2006.

[5]李慶揚(yáng).科學(xué)計(jì)算方法基礎(chǔ)[M].北京:清華大學(xué)出版社,2006.

[6]任英.變量可分離函數(shù)的重要性質(zhì)和應(yīng)用[J].大連輕工學(xué)院學(xué)報(bào),1992,12(30)：66-71.

[7]丁晶,高榮松,鄧育仁.隨機(jī)水文學(xué)[J].四川水力發(fā)電,1985(2)：128-129.

[8]芮孝芳.流域水文模型研究中的若干問題[J].水科學(xué)進(jìn)展,1997,8(1):94-98.

[9]李致家，于莎莎，李巧玲，等.降雨-徑流關(guān)系的區(qū)域規(guī)律[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版)， 2012,40(6):598-603.

[10]包為民.水文預(yù)報(bào)[M].北京:中國(guó)水利水電出版社，2009:25-28,328-329.

[11]倪雅茜,張文華,郭生練.流量過程線分割方法的分析探討[J].水文,2005,25(3):10-14.

[12]RAY K,LINSLEY J,MAX A,et al.Hydrology for Engineers[M].3rd ed.New York:McGraw-Hill,Inc，1982：120-140.

[13]龍胤慧,廖梓龍.基于泰森多邊形法的慶陽(yáng)市面雨量計(jì)算[J].河北工程大學(xué)學(xué)報(bào)(自然科學(xué)版),2012(3):64-67.

[14]張德偉,崔永生.在計(jì)算機(jī)上應(yīng)用泰森多邊形法計(jì)算流域平均面雨量[J].水文,1991(1):52-53.

[15]王玉德.基于ArcGIS的泰森多邊形法計(jì)算區(qū)域平均雨量[J].吉林水利,2014(6):58-60.

[16]趙良民.泰森多邊形法計(jì)算機(jī)程序的研制[J].水利科技,2005(1):55-56.

[17]魯南,菅瑞卿,劉洪波.流域變雨量站泰森多邊形法的微機(jī)處理原理[J].黃河水利職業(yè)技術(shù)學(xué)院學(xué)報(bào),2001(1):1-2.

[18]楊瑋,韋宏鵠.利用泰森多邊形法分析汾河流域中段降水量[J].山西水利,2006(2):72-73.

(責(zé)任編輯劉舸)

Methods Analysis on Rainfall-Runoff of Semi-Arid Basin Relationship Based on Regression Model

YU Chao1, CHAO Li-jun2, ZHANG Meng-jie3, ZHU Jia-chang4,SUN Jing-zhi5

(1.College of Science, Hohai University, Nanjing 210098,China;2.Wentian College, Hohai University, Maanshan 243031, China;3.College of Hydrology and Water Resources, Hohai University, Nanjing 210098, China;4.Tianjin Water Authority Water Conservancy Department, Tianjin 300074,China;5. Liaoning Provincial Department of Water Resources, Shenyang 110003,China)

Abstract:The relationship between rainfall and runoff of the semi-arid region was revised by using the regression model. Taking the Dage basin that is a semi-arid basin for example, three methods were used to build the regression equation. The percent of pass, the permissible error and the deterministic coefficient of the three methods were calculated and compared to evaluate the performance of them. The results indicate that the third method shows the best precision in Dage basin and has a good description of the relationship between rainfall and runoff in semiarid region. It is an effective approach to revise the relationship between rainfall and runoff in semiarid region.

Key words:regression model; rainfall-runoff; correlation; semi-arid basin; evaluation method; multi-factor; separable variable

文章編號(hào)：1674-8425(2016)02-0058-08

中圖分類號(hào)：P338

文獻(xiàn)標(biāo)識(shí)碼：A

doi:10.3969/j.issn.1674-8425(z).2016.02.011

作者簡(jiǎn)介：余超(1989—)，男，江西人，碩士研究生，主要從事水文預(yù)報(bào)和流域水文模型研究；通訊作者 晁麗君(1987—)，女，陜西人，主要從事水文預(yù)報(bào)和流域水文模型研究。

基金項(xiàng)目：國(guó)家自然科學(xué)基金資助項(xiàng)目(41130639，51179045, 41101017，41201028)；院級(jí)科研項(xiàng)目(WT14001ZD)

收稿日期：2015-03-27

引用格式：余超，晁麗君，張夢(mèng)婕，等.基于回歸模型的方法分析半干旱流域的降雨徑流關(guān)系[J].重慶理工大學(xué)學(xué)報(bào)(自然科學(xué)版)，2016(2):58-65.

Citation format：YU Chao, CHAO Li-jun, ZHANG Meng-jie,et al.Methods Analysis on Rainfall-Runoff of Semi-Arid Basin Relationship Based on Regression Model[J].Journal of Chongqing University of Technology(Natural Science)，2016(2):58-65.