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

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

流域年均含沙量的PP回归预测水利论文5Vzlq__\nS

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

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

关键词流域 输沙量 投影寻踪回归 预测水利论文mQ| TdsD Q~

1 引言水利论文Q%Y(WR } sO

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

2 PPR原理及算法简介

w'y \0K5]0Q"I}0

     投影寻踪是国际统计界70年代兴起的高新技术,是应用数学,统计学和计算机技术的交叉学科,属前沿领域。水利论文T&@/B\\$P;\;| a$l

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

oS n'N0i-|'\-vBbD0

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

990112e1.gif (1513 bytes)

?Z'Vw5~D5h4B0

(1)水利论文;Df&X6R.M@+n$IKb

式中Gm(Z)表示第m个岭函数,Z=(990112e2.gif (928 bytes))为岭函数的自变量,它是向量990112e3.gif (888 bytes)990112e4.gif (879 bytes)方向上的投影,990112e4.gif (879 bytes)也为某方向的P维向量,M为岭函数的个数。水利论文E%o%qrt:T9{ DItHv

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

'`-\t5|@0

990112e5.gif (1470 bytes)

R}([/EI#V0

(2)

WWM^%@0jA4o0

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

  ①给定一个初始模型;

-l0J*B P%~TN0

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

a#A/nLp/hX ^0

     ③把上述投影中所包含的结构并到现有模型上,得到改进了的新模型;

J d-o u?v6b2X0

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

2T"{w [ Ms0

     模型的关键是最终估计出式(1)的参数:岭函数最优化项数Mu,岭函数Gm,系数αjm,βm。其判别准则仍是:选择适当的参数组合,使式水利论文 @u ?r"l

990112e6.gif (1600 bytes)水利论文0b_P o\%c

(3)

9z#\ p3yJ/o6`0

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

6L\&qp n'Oy0~6K0

     用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]。

W)[NE0Fk7N0

3 某流域的年均含沙量的PPR回归预测水利论文9H#Y)t7C)F@3y u |0H

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

     该流域的年均含沙量及有关的4个因子的历史监测数据见表1,数据取自文献[6]。水利论文+m@X:N`s0~#qE

     用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个因子对年均含沙量影响从大到小依次是年径流量>降雨量>采伐面积>采伐量。

0s0s/k:k.K5Q:p0

表1 某流域年均含沙量及有关因子实测值

_4A)S?U6i-Y0

Yearly average sediment concentration and measured data of related factors

T&F`A$U,Vc0

序号

`Vj{,`w4sN4rc0

年份

7c,i xS~z,Qu0~m0

采伐面积

nI%jmsV0

X1(M2水利论文} uHn!c.RCW^

采伐量

*_k!k!H xr0

X2(M3

a]2d'yb'k"q(W/l0

降雨量

^x.v&z$W.K(J0

X3(mm)水利论文j&k/Dm.VFKY-g*c

年平均流量水利论文YR'vsz$qUO

X4(m3/s)水利论文 f9L"cel5D

年平均含沙量

,o,s.eeW5a&@0

Y(kg/m3水利论文-PE}5\1DbA:t'`"T


1

Q u6g8Vt)V8W$zP0

1961

f Db/zU(o8AgEz6o0

15100水利论文3c,lT?Xb/f)X$B

29000

n(^ P9B s D:VdyH0

508.0水利论文z6g_/|~U dD\g

180水利论文,VlWIv

3.40

,s'S0g7e u+lcv#Br0

2

yU#Y!R$Y-o0

1962

t2M@0Gu.e4dv0

20500水利论文9@KF:d@uA8WRHt

43350

7f/AC{4Z*c~Xq7N[0

453.7

+@.C!R TqV!VC0

141

^Z:N6L~ i'A@0

3.00

u*r`"p5IZ0

3

\Zt D BZLB"@0

1963水利论文}y2A2wm K%dx

80000水利论文~C K\k7J m8f u

92610水利论文7@N9c*a-xR

487.9

'_$K4c*a\)x8`0

132

'kh,VU1H2E6A0

2.70水利论文gkV&`}

4

R9M`@ZH8\ E+d0

1964

2z6T)Y |x$n0

91500水利论文4~+n:fa7^B

142775水利论文F8]ZH8Sp:G

572.3水利论文1Ng^W j R8zx

182

GVt+p$iL"E0

3.37水利论文[b4ZmX-q6Z

5

s |d8lC d0

1965水利论文O Lr)bu cby

82500水利论文z8BL`+z/dH;l

2123160水利论文 _3u,O4iqTm,f4u

455.7

l,Kc+{@yL2S0

113

cp6huh:Zs0

6.894

9s0|_|6n4G5]0

6

+\f Fu M4eo1q'R0

1966

U z)_7U3I n0

20000

-U$C/c }5L+`K L*ux0

227800

BO @6i$l I0f0

481.3水利论文)~2z.wZje D4_

170水利论文,MPd V,e%p H

5.03水利论文nq4u [.t.[

7

T{ j#WTPM5n g&Q0

1967

:e]8[4T@:zkwL0

17800水利论文}x"Rl!eo^6pw

140000水利论文"yWZ,v)?~ Q%I g

541.3水利论文$A8eU+mi:M*k

199

zM:Zu]~R r\0

3.55

#[ pV+Wwi5RL6}1f1I0

8水利论文h7Xi sn-c7Z

1968水利论文Pk/ko7v;c!r

3900水利论文 f9?;g9o J:N/m+q,O

15980水利论文"^v/d0K w

538.6

^1D?w%iE0

186

7J&A_9CO Z0{F0

2.72水利论文mIgKo%o*vz5r

9水利论文RV t.B-F

1969水利论文^%^;b0j6ve8u

17300水利论文B7S a^3h0b A

223200水利论文Y#q$Z3Ws1E7~g{

460.6

0Cw3f Ml/G brx0

99.8水利论文 E$LJ:kA:D

4.05水利论文&W.s5dU E2]&c

10水利论文d)a0dTB#C

1970

-| YGH:P0

25700水利论文 g|3O#@X7T\

229400

?d'K/U3h0

393.1水利论文O*?M)_.p/JM

133水利论文 k1Zq0|,Ju XIp

3.22

3{l,Xl'?'vs*K0

11水利论文 V|E"dkl|q e

1971

dh%K V X8cNT~0

49400

Oa!uu{"W8KP-s G,O3`0

424500

j-gQ0h@-ne:f*F0

373.9

%W6Q.c1yX0

106

PT1`'jZX2C(i)]0

2.65

Ts sxr0Xc0

12

"RU/oPi,a ^H+L0

1972水利论文'b s1`(L_4zqQ

40700水利论文l`6MO(X7e ]

561700水利论文2F RyY+]#VV#C

428.8水利论文%y{ bbL!UX*a

107水利论文#AQc FfBk

1.91

Uzj%~ M4U0

13水利论文cb"S)S'D%{|

1973水利论文YiH }yAgzq V

77000

;b{ N-YC6J9d cs0

563600

.Mf'Qo.jn8Um0

482.1水利论文'W"O_B)u"l6|P

140水利论文2gZ5K7fP!\Y

3.00

*K/J-e'@&Y T0

14水利论文rX*DG/g&ht

1974

Fy5E;t eRu0

72900

)f OTI"W.qi0

557600水利论文 `0fp0b^0O

415.1水利论文;dE'~2[6GjP

121

Nt&L$t#W K#l L_%}.{'[0

1.31

n+\*U`j@0

15

@ke/S$w Pp0

1975

+o {4U4R-C0

63300

ep ` JC#`m ld hl*^0

528300水利论文/k#ZZA%Ac \*JY5J

536.7

,l&J1E/x5Nx ~;bZ0

144水利论文Db.qf M)V'I

2.33水利论文 sb9_a^wr(q

16

7p p$~5N f]H0

1976

+C'o/p_Ug eO*A0

51600水利论文CPV0P5}A+i{5rZ

488940

e+yu^fQ`j0

385.1

{#}J4tS{df v0

154

WF U+X'eC/z&^,_Y0

3.55水利论文5_9kJfQ~

17水利论文$eoCwFi_)N

1977

{"J%y-YH*]4^4|0

60000

C.FMv7N#SO P0p7U0

480500

f3n7R*Z Vwd.q0

412.2水利论文6zj!l&U"Q)q+^

111

Z[G0xpl5q0

3.35

EP/OPN!L0

18

eD2Ws A:ESv0

1978水利论文8_B L;}1V)v&U%S

70000水利论文9Q5d7v$};z/Y7L)\

530500水利论文"W`m,|C)y&V+I5aS

567.1

!IW {Y1Pyr[0

139

3~eJq;SN0~8R0U D+|^0

2.55水利论文5an)Q|jF


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

^+x"j(i%D!f4~0

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

3][q0S[{*u0

实测值

#G S-Bh+m8r0

计算值水利论文8p Tw6}8X;H%T/t

绝对误差水利论文|nU{ [,q\

相对误差(%)

:pl/?zi g0al0

实测值

|5op?#dC od0

计算值

1_s7H{ML0

绝对误差水利论文8t$E?'\$SW

相对误差(%)

^4q2C(y A/_jm?0

3.40水利论文vR6i&q&aIf$k3}!`

3.131

(p` z$K6q!sz0

-0.087水利论文kA;V#CS5T1bPL

-2.5

w{9Ew B$U#Ap'QV ]0

2.72水利论文&K_;B H)en0{/a?s!f B

2.934水利论文6cx,\,s+of

0.214

_:ieZ,{ DOQ0

7.9水利论文MhvM'F:Ud5A9s

3.00水利论文F(V `r&^-@@)x\

2.738水利论文*r'M;kS:}Kj9w

-0.262水利论文S/dFDW+u|H L;a

-8.7

p P(H&U-p0r0

4.05水利论文G&Bg'|NNp`ibc

3.996水利论文9@+TT9a f

-0.054

M `T*x6ONLs0

-1.3

j+w2_+Nq H0

2.70水利论文Ec/X9vz:w@9c

2.659水利论文J2\.i2y n

-0.041

;F9M|s {1r }/n0

-1.5水利论文Z ?vF&O

3.22水利论文.As}aF&X

3.214水利论文tn)i;M2{}(D%|%p

-0.006水利论文6G;gML3b LM"A

-0.2水利论文q L6bBm

3.37水利论文4V$O9U5O@J6A

3.400

+d%s]7I.rx(@#Y]0

0.030

VKJ`x\#?z0

0.9水利论文)K#dE qeK?)x

2.65

6l+_]${xt&^)I0

2.628

m]5iZ-x}P0

0.022水利论文.fyoH^!\:Uu.[.i,R8I

-0.8水利论文m&^5`'I#BV_w vr

6.894水利论文P]-Qt IiT&x

6.926

r2R{,a'v? q0

0.032

xG%oJ/Q{ KVd0

0.5

&E.B}WHUBS!~0

1.91水利论文Ur?!K2` ~y,vK

2.135水利论文E xk6nhOh4q?

0.225

fz6\)q.Z:O0

11.8

$f&D7L6|w0

5.03

J{7j4ob `4P nL{0

4.729水利论文&\.Jiy&q^n

-0.301

^'_K0@)B1g-q7j ]!X0

-6.0

(z-{9e9X gV0

3.00

XMG,a8n0

2.987

(NJQ Y"c0t0

-0.013

!@,R;w w8Ke l%B0

-0.4

s\ I5O4S!sz-I X0

3.55水利论文-QiJ ?I(y F kr

3.826

5n!S)D|}2D~!H4i#ma0

0.276水利论文+R3gt-C T`zjm

7.8

UU _/l[!X0

1.31

$FHi'e"hM6y$T}0

1.319水利论文3Z u5H]XE

0.009

ja,b5}v\1Nxj0

0.7水利论文x3H}|-sW P

2.33

OS xj_u {\+T0

2.461

t"R[oY&ca0

0.131水利论文:s4~)E{-iM,Vt

5.6

9|%\/bon#F$RG0

3.35水利论文:plSka^H

2.780

)c U%t u)zA5U'?e0

-0.570

,CZ9v3d;H9f(_V#W]1zC0

-17.0水利论文 I ^#_"[_7XB

3.55

:M4ST]&b!B6Ia0

2.869水利论文?Z9c$q4l

-0.681

6MM!\*S-k3G0

-19.2

0g)@5rrg0

2.55

w6N&Zo:Q2]m0

4.305水利论文 R t4[%L1M7K}6p``

1.755水利论文#}(_Yso8H+P.~aP

68.8

|RK ^ r!Q7^_!y q0

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

表3 各因子的相对权值水利论文:j,AyC9K

Relative weighted values of factors

9Q)[+`(a*Hi-?0

权序水利论文_-{J`C G

1

N bygr/z/d0

2

ahpJkY/o%H0

3

0ph*V.hK5g&K|#`-F0

4

NP7R @ u,L&z0

因子

kx'v~7Y0

X4水利论文D7V;f5i7m{

X3水利论文.J o7VMS,TEE

X1

b T$^4b$]0

X2

m9{7T*Ju.Cx |0

相对权值水利论文v@'x kXg d!b

1.00000

$DArxO%G%B&D0

0.85190水利论文 N:H5?1{e?$v

0.71617

H |"v%B$i(Q0

0.29765水利论文+h j9s0s)C:g g.w


     若以相对误差δ<20%算合格,从表2可见,年均含沙量PPR预测模型拟合样本数为全部合格,其拟合合格率为100%,而预留4个检验样本合格率为75%,其拟合和预留检验效果如图1所示。   

["SDVI/H0

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

#F(EGh[Lei_D0

4 结论

0c?v Y }O wX M0

990112t1.gif (1658 bytes)

_CMu;eM0

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

)j%e6C1?jc X(x0

Effects of fitting and verification of sediment yield水利论文 \ Cn\g$Y n \

in a watershed by using PPR model

$G^k^ _l0

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

J$\Puo0

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

8p6z4nP5QV,]0

     3.PP回归用于流域输沙量等预测建模,只需直接使用原始监测数据,不需要对数据作预处理,避免了人为干预,客观性好。水利论文:d7Le'h*}$R

     4.PP回归也可用于水土流失中其它指标的预测建模。水利论文-b)F G1{|hO

   5.本文用PP回归技术建立了流域年均含沙量的预测模型,求得年均含沙量S后,只要将年径流量QT×年均含沙量S就可得到流域的年输沙量QsT。水利论文,P&y MC)v"v zYujm

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

9S,X9n!H:_$z q0

参 考 文 献

)Qxj"_u iN0

1 陈浩。水土保持学报。1992,6(2):17-23.水利论文 D gn4e/B#|,x

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

Yy%p+T1o"RRB` alS0

3 Friedman J.H.&. Stuetzle W.,J. Amer. Statis. Assoc. 1981, 76:817.水利论文4|,P-Q Nm A_

4 Huber P.J..Ann. Statistics, 1985,13:435.

;k [8^/` A IE:W0

5 郑祖国。水文。1994,(4):6-10.

@U%U4@:lB0

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

5Whv?8}]0i,v {Y0

 水利论文,S9o Of0aO-u&C

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

广告投放

广告投放