流域年均含沙量的PP回归预测(李祚泳 邓新民 侯宇光)

热度196票 浏览64次 【共0条评论】【我要评论 时间:1999年1月01日 14:33

流域年均含沙量的PP回归预测水利论文5A*_wqB"y!P/S

李祚泳 邓新民侯宇光
(成都气象学院)(四川联合大学)

摘 要应用投影寻踪回归技术,建立了流域年均含沙量的预测模型。用降雨量和年平均径流等4个因子建立的某流域平均含沙量的PPR预测结果的拟合合格率达100%,预留检验样本报准率为75%,表明PPR用于泥沙输移规律的预测研究是可行的。水利论文8E9G$v*w b5tp5S

关键词流域 输沙量 投影寻踪回归 预测

q+O Qo2@0u"q3p)d0

1 引言

3n ^p$Y5V0

  我国是一个水土流失严重的国家。严重的水土流失给工农业生产和国民经济建设造成巨大危害。产沙量是反映水土流失的一个重要指标。而气象要素、地形、土质状况、植被系统及人类活动均对产沙量有重要的影响。国内外不少学者针对不同的地域特征,对流域的产沙机理,泥沙输移规律及其防治对策等进行了大量的研究[1,2],但对诸多因子与产沙量之间的定量研究进行得较少。考虑到引起泥沙流失的诸多因子与产沙量之间的关系具有高维和非线性的特点,而传统的统计预报方法是采用“从某些假定出发,按照一定准则,找出最优拟合”这样一条途径,难以适应千变万化的客观世界,也就无法真正找出数据的内在规律。这种传统的预报方案往往是还原拟合较好,但预留实况检验的精度很差。近20年来,在统计学中提出了一条“审视数据,模拟,预报”称为探索性数据分析(EDA)新途径。本文正是采用基于这种新思路,应用投影寻踪回归技术(PPR),建立流域产沙量的多因子预报模型。

2\rPbf].m,Ol0

2 PPR原理及算法简介

(_D)Nk${9ig3`O0

     投影寻踪是国际统计界70年代兴起的高新技术,是应用数学,统计学和计算机技术的交叉学科,属前沿领域。

6n)iL#MlPcJ0

     PP是用来分析和处理高维数据,尤其是来自非正态总体高维数据的一类新兴统计方法[3]。其基本思想是:利用计算机技术把高维数据通过某种组合投影到低维子空间上,寻找出能反映原高维数据结构或特征的投影,在低维上对数据结构进行分析,以达到分析研究高维数据的目的。水利论文^)o I/Xgu0lJ.uj2A

     传统的诸多线性模型通常局限于正态分布,但多数实际问题却不呈线性,因此勉强用线性手段进行辩识和预报,很难取得好的效果。而PP与其它非参数法一样,它可用来解决某些非线性问题。它虽然是以数据线性投影为基础,但它寻找的是线性投影中的非线性结构。因此,它可用来解决一定程度的非线性问题。PPR模型如下:设X=(X1\:\:XP)是一P维随机向量,Y=f(X)是一维随机变量,为了避免线性回归不能反映实际非线性情况的矛盾,PPR采用一系列岭函数的和来逼近回归函数的方法,即水利论文R$`I"aT3A d0C

990112e1.gif (1513 bytes)

'}s Xew].jmx0

(1)

AciYYpNK/v0

式中Gm(Z)表示第m个岭函数,Z=(990112e2.gif (928 bytes))为岭函数的自变量,它是向量990112e3.gif (888 bytes)990112e4.gif (879 bytes)方向上的投影,990112e4.gif (879 bytes)也为某方向的P维向量,M为岭函数的个数。

bdF#{C!tj0

Friedman和Stuetzle提出了实现PPR的SMART多重平滑回归技术,SMART模型具有如下形式

r4JG X Xf'T0

990112e5.gif (1470 bytes)

/}TzT jY0

(2)

"b0X%n^3Z0H Y!E'~0

  它实际上是采用分层分组迭代交替优化方法对式(2)中的参数α,β,Mu和岭函数Gm寻优。实现步骤为

6|7S v7iIMo0

  ①给定一个初始模型;水利论文6UW tM*E g

  ②把数据投影到一个低维空间上,找出数据与现有模型相差最大的投影,这就表明在这个投影中含有现有模型中没有反映的结构;水利论文 [y ybwL|&UV |

     ③把上述投影中所包含的结构并到现有模型上,得到改进了的新模型;水利论文,P.A*XEY{']{

     ④再从这个新模型出发,重复以上步骤,直到数据与模型在任何投影空间都没有明显的差别为止。

d[)G~6@EF0

     模型的关键是最终估计出式(1)的参数:岭函数最优化项数Mu,岭函数Gm,系数αjm,βm。其判别准则仍是:选择适当的参数组合,使式水利论文'{p:g/@7?QD| E M

990112e6.gif (1600 bytes)

$`t2b6r`3ur9yw5Vd Qr0

(3)

p e&fh&MQ0

具体作法是:把全体参数分成几组,除其中一组外,都给定一初值,然后对留下的一组参数寻优。得到结果后,把这一组参数的极值点作初值,另选一组参数在这一初值下寻优。多次重复直到参数收敛为止,即将αjm,j=1,2…………P,βm及岭函数Gm划入一组,m=1,2,………M,共有M组。固定其中的M-1组,而对这一组ajm,βm,Gm优化求解。此时,又将其分成三个子组,分别固定其中的两个子组,对第3子组优化。然后重复这一过程,直到收敛为止,即L2不再减小为止。

mS4J |&r0? o0

     用SMART计算软件进行PPR计算时,实际操作十分简单,因为模型只有岭函数的光滑系数S,岭函数个数的上限M及最优个数Mu3个参数需要调整。光滑系数S确定着数据平滑窗口的大小,其取值范围S∈[0.1,0.9],考虑到在满足一定精度条件下,尽可能缩短模型运行时间,要求Mu≤M≤9,通常取下面这些组合:“M,Mu”为“5,3”;“6,5”;“6,4”和“9,6”等,其最佳值最终由计算结果分析确定。此外,样本容量值N可在不超过500范围内自行指定,屏显值与指定值N之差由程序自动视为预留检验样本数。PPR的理论分析及算法详见文献[4,5]。水利论文OAU R0]4V-k2f1WJ"z

3 某流域的年均含沙量的PPR回归预测

AV)H5?:ct/Jz0

     一般说来,流域的年均含沙量受多种因素影响,但对具体流域必须具体分析。其因子选择原则应尽可能选择那些与年均含沙量相关性好,而又比较稳定的有明确物理成因的因子作为PP回归建模因子。经分析,引起该流域含沙量增加的突发性偶然因素发生的可能性较小,年均含沙量比较稳定。该流域含沙量与气象、水文要素和下垫面植被状况关系密切。因此,选择了采伐面积(X1),采伐量(X2),降雨量(X3)和年平均径流(X4)等4个因子作该流域产沙量的预报因子。水利论文?u8E-I@a#l(B

     该流域的年均含沙量及有关的4个因子的历史监测数据见表1,数据取自文献[6]。水利论文/TA~dWJR

     用PPR的SMART软件建立该流域年均含沙量的预测模型时,用表1中的前14组数据建模,后4组数据预留检验。将表1中全部18组数据输入PPR计算程序,在固定因子数P=4和建模样本数N=14情况下,反复调试模型中3个参数S.M和Mu的不同组合,使模型计算输出的拟合检验效果达最佳。由于式(2)中的岭函数Gm(Z)不是解析函数,而是数值函数表,因此,PP回归最终只给出达效果最佳时的模型参数组合值为S=0.1,M=6,Mu=4,此时年均含沙量的拟合和预留样本检验效果见表2。表2的后两行为4个预留样本检验结果。由于用PP回归计算的拟合值是多重平滑后的结果,因此,拟合值不是时序值,它与时序值有一个差异。不过这种拟合差异如用相对误差绝对值的平均值表示,不超过4%。用SMART计算软件实现该流域年均含沙量预测建模过程中,还同时给出如表3所示的各因子对年均含沙量的贡献大小的相对权值。这些相对权值就反映了这些因子对流域年均含沙量的影响程度。从表3可见,4个因子对年均含沙量影响从大到小依次是年径流量>降雨量>采伐面积>采伐量。水利论文1b n ?_*?X

表1 某流域年均含沙量及有关因子实测值水利论文 dBo1}7p

Yearly average sediment concentration and measured data of related factors

w@"e s1\I+Xsv0

序号

,\^$?[ b{0

年份

9k+~(r4n6r$[0

采伐面积

'YOr'Y,d1aayG1|0

X1(M2水利论文$W$y8jo.N#\{e

采伐量水利论文qX n)@4G

X2(M3

@3i1d.c w?+b0

降雨量水利论文#wh y_6za,sY}

X3(mm)水利论文[ h]aJG5HI1{

年平均流量

YS^@N6w"Fw`0

X4(m3/s)

a6}K'e7dt6s0

年平均含沙量

$X5R^"t5D0

Y(kg/m3

H7GV-g'w j0

1水利论文)}uz nqb8N

1961水利论文,k$Wv%p wA,SU-o*dj L)q

15100水利论文,Mu^B0zQ6c

29000

Rd/Dmul0

508.0水利论文;U4`+C7fo@

180

{2m|$S7?,{ ON I0

3.40水利论文2Z l;q.Y(V#|8r+`Q

2水利论文 {yn#]4IT {

1962水利论文3D6for`EK'MHmC

20500水利论文lm,qibl:FZ

43350

g!qs@5od9Y'I0

453.7

Z^IR2Ck'x0

141

4o}Y;}9u+|.Ut0

3.00

:}+o%[/[-pP7`I0

3

j+xV;v`9h?4m|0

1963

"D:q)b E$v$aR:vR0

80000水利论文 N'Pg3C H/E#ez4^

92610

;D-W]2Rq!p0

487.9水利论文AD*u8KmM(WL

132水利论文$h*@%@}]&mL

2.70

~S*N(z/b5w4d${Ln$H0

4

/yX SJ1c0

1964

1}.y#|8\5Yj_|9zW0

91500

jY4_(Gn3~0

142775

-J,mU,?^r0

572.3

9G%FCi!{0S0

182

LE(L1t:Jg*]?&T8]0

3.37

d?;QU J;O9l0

5水利论文U)q2gxsO c1X'b7[}L

1965水利论文#c7x3^ ca4\%PGZ

82500水利论文sgGh!t,fZ

2123160水利论文 o)oHB&[n'V

455.7

*C,Ev-sF.?9`t0

113

3OT"N6t1W.nbS'Y0

6.894

DdL(FC9\$Yw#i0

6

,w'z$e pq+b1r0

1966水利论文t?u1nd'X

20000水利论文E0|?H)|Z_l/B

227800水利论文G zM9dPq]

481.3水利论文 ?+h2p.TIw(BF

170

S} c)o _*eZ.Lz0

5.03水利论文~h;Ar,|#N ^P^Ex

7水利论文$MT5g.`&D:O n?I

1967

5b9xR+]R@#Y.p\#[0

17800

/\6uUy!kw T0

140000水利论文W4@,? p\D(S7b-V1Q*C

541.3

@ `"C h `+~0

199

)r6C9|v+rg3@%c;R0

3.55水利论文JQ#_ ^fpN4j

8水利论文1w*m8BXk'Z'f;jzY We

1968

f}pp$M7L%W8j)E3`A0

3900

*?a6I@5O0

15980

~4ojn-jEi0

538.6

[.v0T"Z I"rN+j? V0

186水利论文5C n#rxjh2~ @ g:BE a A

2.72水利论文~i2\*[#I8g

9

Xc"h5L n(U0

1969水利论文5ZUR h0kC3C8` L!{b

17300

6g hq\~"f*a"F0

223200水利论文g4`t#U f1~p

460.6水利论文 p~h Q7j{%fF

99.8水利论文#Ey S1`Fr+HS,S4l

4.05

C;K(tDx ^9F0

10水利论文8n4C6?2qg'I$h9r*?

1970

dAf { lr ?t0

25700

5yMXn!M{!sUO0

229400水利论文|!b1`)YR*G)Z

393.1水利论文Eke)X b-P4?(E

133水利论文5c O9z? k@I$_|

3.22

Fx-B{ T0

11水利论文v b5HkfJ\6u.O

1971

-])N#h"f$K~e-U0

49400

&_8BM;b6}H7L \(G0

424500

!A7O;ZEa:h0

373.9

+Q4zmX9?,xz;oi0

106水利论文)nr+}mjVo

2.65

Zm2M8D}h%N:n1~V0

12水利论文5{In |%t0Y

1972水利论文h'h?;yn^'O1b

40700

BW/V[L,QPm0

561700

"_CxS/L8S/[0

428.8

*K tBi7zEg0

107水利论文 ];fnB c'V.S0t

1.91

d3Kb,T NG^0

13

?"?^ ?/l!q[],Hy0

1973水利论文U\K6EL5D/O.XR f

77000

4Nw.].Al!i0

563600

~| V4{5oLT0

482.1水利论文A)F!R;m3g K

140水利论文/tNx1C0s0L-V7j4B,xD

3.00

6d0n5_mE @0

14

PML5j6A:H;w.K0

1974水利论文2O'@H%e0x8BJN

72900

$bsL npS Njh"M#Q0

557600水利论文@&r@ ?"E.oV"_

415.1水利论文:r9EqQS k1JS u%\$sS*v

121水利论文lt@9^'Q

1.31

F UqSe1\-E0

15

DT7s7`Akm C2n0

1975水利论文D6b6H6rn_8Cu(z[

63300水利论文9E l"s k3IC B3|

528300水利论文4R~C6Y{[T z

536.7

H)YG,l f*[w0

144

3u])|8LPDl0

2.33水利论文{@~ o{u"\

16水利论文s7ir*_;T V:Qdj

1976

[xW4Fj0

51600水利论文})YB G P"v6V

488940

SMO W|6e0

385.1水利论文mck`5R!C

154

X P!So5IB@7X%s ~0

3.55水利论文Ie&c]Vd'b{*O

17

M(W[rK0mR:o0

1977水利论文 ]D,Wezx?

60000

%ZqbaP0So!c&~0

480500

'\ yM%o.Oe:n0

412.2水利论文1n/mG:U/]Ex

111水利论文!V}jK8@7V)?*j|

3.35

6Z-Lz"Uh#t0

18

]2o/nR7C2E+z%sY?0

1978

g&Ly*l(};A PZva0

70000水利论文-FFP-sQa"u f

530500水利论文g1}hhV\'Z1n Vr'H

567.1水利论文 lw/n j1B,`E2Z

139水利论文 w4T1xx)?

2.55水利论文Z jc9`#l'f


表2 年平均含沙量PPR模型拟合和预留检验结果

6X.c"Kol0

Results of fitting verification of yearly average concentration by using PPR model

"H#M;]dv*Sf0

实测值

/wbC+^#E0

计算值水利论文TU$WQ B#L&d/j

绝对误差

TD{}7b&] ^0

相对误差(%)

nUt$o W+E }0

实测值水利论文:UO6DNZ j1d {

计算值水利论文 _V"b4T(s/i

绝对误差水利论文KjUt,zX6I

相对误差(%)水利论文bBj8m!a gTp,O0g


3.40水利论文sH kE&@`

3.131水利论文ng*v;| [M

-0.087水利论文B,?2ip ftB"W W

-2.5水利论文kn*w8Y&c.e4ub

2.72水利论文/bhJE,Z"pQ M

2.934水利论文,]&JIg/p7e,r

0.214水利论文 CBf C"C7Yi [

7.9水利论文ir%i uBb\3]"t

3.00

IESc` S Voj0

2.738

B PALB aw,z0

-0.262水利论文M9]?k,s+]n

-8.7

3u{'a G5`6G)U0

4.05水利论文)x!l E\&P F&J)h

3.996

jV UB:tp*Q0

-0.054

.HnCANT8n z m0

-1.3

GcoXzwE0

2.70水利论文9d7k,w0HR R

2.659水利论文4n@ d G5CHt

-0.041

V E*E e`b0

-1.5

6mZi6~5Pr0

3.22水利论文\2l Cv O

3.214水利论文)Rl:N/Zh Mm(@

-0.006

7}%D|(b |5_ N"~0

-0.2水利论文x"C{ B}/k

3.37水利论文 _7I,};@ S;n1lT

3.400水利论文Z"h3EY ?@

0.030

E8O R!G-WX;T"s0

0.9

,V ^ ^K-XD/Byj~0

2.65

gU1{ ]mk kj:b0

2.628水利论文o)NNr5vgrJ

0.022

z&Z,Qi&o ZM-k0

-0.8

2b1D.|f@'_:ja0

6.894

6J6|"o)T5B}%l%]0

6.926水利论文x `[,za9T;jG

0.032水利论文*MU+t:u)AzS*gi@

0.5水利论文-X.n3X] `

1.91水利论文/s7kEd7O\

2.135水利论文jc#r+{Y` D

0.225水利论文&zn dO"{

11.8

uFC1j K.EJ0

5.03水利论文;L{ Rd(|-`0kD

4.729水利论文giZ~ ]/T.r5k

-0.301

,N$Iv VQ:m5}0

-6.0水利论文y s]9j(X!K'y#?+u5N

3.00水利论文+m-z9JlVp6IrnN

2.987

^ gndUDj\d!|0

-0.013水利论文`/T1B;JF:\:NE)vO9k

-0.4

)DRn)kH:IG+`0

3.55水利论文"R l7W1k]k

3.826

C1iZy5a8X0

0.276

&WBE h/w0

7.8

B1N J3Q5Q&c/U^O@0

1.31水利论文$_xap+q ~Rr

1.319

&Y_c)t%@ d N$[0

0.009水利论文6}ONUX&|V0|

0.7水利论文3sz Z XTP/j

2.33水利论文bg$l,Wv!u

2.461水利论文f+o8c h*_

0.131水利论文K;w];{u)txk;{

5.6

8W G {%@ Lo'Pjz0

3.35水利论文 C?g8hR(J

2.780水利论文Y7J5mQ2B?/Gw;W

-0.570水利论文l_p!loe!w;w:j

-17.0

p7P~+VWI6T ~0

3.55

QT*u:h%?g"p0

2.869

_F I;gXU0

-0.681

Aw%[n8Q\ tM~0

-19.2水利论文f1{Sk\#]T$i

2.55水利论文E/I bqz#w j

4.305

cv9uC4i"HQG0

1.755水利论文 iw+Ka'd0k"Ij S$@)G Y

68.8

/hc0Ba"wz0

拟合合格项数:14 拟合率:100% 预留检验合格项数:3 合格率:75%

表3 各因子的相对权值

H2By[m|aR0

Relative weighted values of factors水利论文Y&} R2X,BUG"X)U|


权序水利论文7srn1G/~n9} c8K\

1水利论文.Kc rD5R}3x"\e,f

2水利论文4B em\0Q/X [

3

$P!w z%FD(vXF0

4

7z;ty,Q7q!h)g+v'd0

因子水利论文#vI,I&Hxn

X4

,Zz3?7_4_-gK*`0

X3水利论文0D1jh c2]`{3~

X1水利论文e)H/{#g#P yQ,W

X2水利论文8Y B3ow"_@f

相对权值水利论文u&od|W

1.00000水利论文 p/f3V8iSQ

0.85190水利论文,}-x4xD7q}j _?0n

0.71617水利论文B4AK&| Fq8u?

0.29765

bm PYX5dQ.r)x}0

     若以相对误差δ<20%算合格,从表2可见,年均含沙量PPR预测模型拟合样本数为全部合格,其拟合合格率为100%,而预留4个检验样本合格率为75%,其拟合和预留检验效果如图1所示。   水利论文/}K s*b"?jn7x"W _1O

  应用PP回归技术建立该流域的年均含沙量预测模型的拟合精度较高,但预留的4个检验效果不太理想。这是因为流域年均含沙量除了受上述几个比较稳定的关系密切的因子影响外,也会受到某些突发性偶然因素比如特大暴雨,滑坡等灾害事故及人为因素的影响。而用PP回归建模时,无法考虑这些偶然因素作用。第18个预测值与实况值相差较大,有可能就是某些突发性因素影响的结果。

6k7n)VR.]1I;x2n%w0

4 结论水利论文Gr u Z ?7g

990112t1.gif (1658 bytes)

{lOh,D.v.^{c*}0

图1 某流域产沙量的PP回归拟合及预留检验效果

;A"a[[|1zy0E0

Effects of fitting and verification of sediment yield水利论文5B`[ p!jw

in a watershed by using PPR model

EA%g9V6w/f%i'U0

     1.PP回归采用了“审视数据—模拟—预测”探索性数据分析新途径,建立的模型稳健性和抗干扰性好,因而预测模型有较高的精度。水利论文On!NSK:`h

     2.PP回归无论对于正态或偏态分布,线性或非线性,独立或非独立分布的数据都能有效地处理。因此,它适用于分析和处理流域含沙量等非正态和非线性等有关问题。

tl zRj0

     3.PP回归用于流域输沙量等预测建模,只需直接使用原始监测数据,不需要对数据作预处理,避免了人为干预,客观性好。水利论文B#rST| t

     4.PP回归也可用于水土流失中其它指标的预测建模。

0{1n(f8X9|)Z,t5{o)^gD0

   5.本文用PP回归技术建立了流域年均含沙量的预测模型,求得年均含沙量S后,只要将年径流量QT×年均含沙量S就可得到流域的年输沙量QsT。

Y['_'F%` F0

     6.本文用PP回归建立的流域年均含沙量预测模型未考虑对含沙量有影响的突发性因子在内,只适用于对含沙量有稳定影响因素的预测建模,如何把这些影响考虑在内,还有待进一步研究。水利论文1z E0e:Zbrc

参 考 文 献水利论文M.@ gRq0HQ@4Iu

1 陈浩。水土保持学报。1992,6(2):17-23.

}1gI7F:RuE0

2 张醒。泥沙研究。1991,(1):15-19.

+|U+kk},JK%}c#ad0

3 Friedman J.H.&. Stuetzle W.,J. Amer. Statis. Assoc. 1981, 76:817.

l5G/s#Iluv0

4 Huber P.J..Ann. Statistics, 1985,13:435.水利论文b`.`kQg$E.k

5 郑祖国。水文。1994,(4):6-10.水利论文9O1o(r4]ZJ0c)TEV

6 吴建明编著。环境统计学。北京:中国环境科学出版社,1991:504.

cB+C7oL#[0

 

A vS K9{/B&u:kW0
TAG: 含沙量 邓新民 侯宇光 李祚泳 流域
顶:28 踩:24
【已经有134人表态】
16票
极差
14票
很差
9票
较差
23票
稍差
26票
稍好
20票
较好
11票
很好
15票
极好
下一篇:海滦河流域水系分形(冯金良 张稳)
上一篇:垂线平均含沙量两种表述方法之比较(黄才安 梅小文)
查看全部回复【已有0位网友发表了看法】

广告投放

广告投放