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

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

流域年均含沙量的PP回归预测

U'WTvd B!Fw0T4_0
李祚泳 邓新民侯宇光
(成都气象学院)(四川联合大学)

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

关键词流域 输沙量 投影寻踪回归 预测水利论文 @9s D/Ax!Zk&d b%H

1 引言水利论文X7Rsz"vt6vU

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

2@ b'P:_w)v0

2 PPR原理及算法简介

z TD _gev:s0

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

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

:A[ `)}1_ uF UP0

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

P Uihl"n%_0

990112e1.gif (1513 bytes)

$G @#?3\1V8j0

(1)

{.W#ZzfH(cr X(q$d0

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

Friedman和Stuetzle提出了实现PPR的SMART多重平滑回归技术,SMART模型具有如下形式水利论文8U:Yb4K!EX(y

990112e5.gif (1470 bytes)水利论文+c1iC7U~"Ys%N t9x

(2)

F,^Yjw-L?0

  它实际上是采用分层分组迭代交替优化方法对式(2)中的参数α,β,Mu和岭函数Gm寻优。实现步骤为水利论文ZJ"Sn;a,h WJ!V3K X

  ①给定一个初始模型;

7F8vP*jpD3T,E3P0

  ②把数据投影到一个低维空间上,找出数据与现有模型相差最大的投影,这就表明在这个投影中含有现有模型中没有反映的结构;水利论文%a1L hD|3T{ L

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

bUK"Ic+n vr4@:af0

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

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

990112e6.gif (1600 bytes)水利论文Y2?]mH4H4\+y

(3)水利论文X|#y_BLY+j6L*m'Aw

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

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

P? ^X.p [0

3 某流域的年均含沙量的PPR回归预测水利论文&eu7U`XlwS:N

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

     该流域的年均含沙量及有关的4个因子的历史监测数据见表1,数据取自文献[6]。水利论文@+N[ G$Y^

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

xK F1_4`B N*|0

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

` E&\W&K1@~7}_ K0

Yearly average sediment concentration and measured data of related factors

Z?+~Q&c!Q4M0

序号水利论文*Z%xIX7v&}_ ~,t

年份水利论文2G6by&[|*^

采伐面积水利论文d a%?3h @ K

X1(M2

c(Rn"C$g0

采伐量

^p?B}Uo:q5@0

X2(M3

S0F1^$[+L!a0

降雨量水利论文`'] O/m}8|%nI

X3(mm)水利论文v4UN7CCfo {

年平均流量

~{ Wv [4^/Z'S5?%d0

X4(m3/s)水利论文)Y+heZ2EAg

年平均含沙量水利论文`lC6y.scDZ9gkl%v

Y(kg/m3水利论文 cMae-K{ Nx*r


1

5t*jm(d6@ {Y;F0

1961

/Pm/L6c3S0

15100水利论文 G,`:b!\#v `6_%Q z.x

29000

'y ~5i9c`F:j:|0

508.0

J3_YCco,?JT0

180水利论文*I{,v4KW(}

3.40

1F]2Kh3g:[@o'o0

2水利论文:S(m_U,['{&H

1962水利论文$U0^n.J$?I)T M

20500

+hY9V)j$I/b0

43350

WW,`*Ug5r-mU0

453.7水利论文Ei,T*sv O

141

(X-vd%nxz2~ p0

3.00水利论文0N ?MZX+zC'Gn

3水利论文)s@4i!b4]"N,TCr

1963水利论文l]o&^/E1R

80000水利论文fH0Q&C;M*LuY

92610

'E,jv b)k0

487.9水利论文5L L\I1e!C^ Jk

132

3B5F$zm0Y0

2.70水利论文H`u lvt

4

4A9o'Bw-wB0

1964

V E'ys.y+GLI0

91500

2rJD3ug,g[A C0

142775水利论文6cHv d[ty%R|

572.3

0U2r D cj]r0

182

XA(x1Y8w`,v'{w0

3.37

G6[%yH(e0

5水利论文 ^e { l-?GJX8Y|"t

1965水利论文}$E/i8L"M z5~h

82500

k4N&V^uf{0

2123160水利论文 Q6bwe9`zU x&K

455.7水利论文nB4E!Go@

113水利论文W h*W'[i$G1\l"m

6.894

"sVN`.\(\[@h0

6

o7A:{AP1f%Sq_SB0

1966水利论文*C0O'u Cx

20000水利论文:uH6TmDM

227800

*yGQuq#]CDS#Y;p0

481.3

i4sxb#Y[^0

170

(_!\Y4ld@0

5.03水利论文"ng L+]/zulE!w5?

7水利论文3L(y&z p"^s1V gX

1967

2ib1R T3n$Y0

17800水利论文"]'JX#VdKX3L

140000水利论文zK3bQ?JaaA0J4s;n)X

541.3水利论文7Fec;rY

199

G\\`:[(\-BZ0

3.55水利论文2}5u7Mp%^7z"c:MQ1D

8

?NG]r,E0

1968

f7Tv.p!Qv8Ag0

3900水利论文c-i(nh}P&a,L!l_

15980水利论文8_,O_RWp5Jk

538.6水利论文)~;T-aX7Tl

186水利论文g3crV"e a O3W

2.72

s g(q2_\jq4{!k0

9水利论文:^:c.Vq3^4m:GK

1969

nJD ?B/_0

17300水利论文,L5Vw z \{6Y)IbU

223200

!\:r"l{{(?-V+|0

460.6水利论文ik1?XA$PH

99.8水利论文_&]-^U |)V

4.05

4k/Zm%Ht?;M0

10水利论文 S?8\']g8U Y

1970

k-sw'D8o3RZ0

25700水利论文#Vb"e2lpl(w

229400水利论文H'K C-Ukt)U"[S

393.1水利论文 I t|)T6v u2l9~;EYW H

133

!yt2JxK|-H7?3s0

3.22水利论文"S T W6jpH

11

/X V kr'E5tXM~!x0

1971水利论文HG I$q0tZH'MHm

49400

8M8h)R;u{3j8ZAN e3z{0

424500水利论文|;Iy*]8P

373.9

X;H(kX!C$s h'X0

106水利论文Ix,G bZ)lu

2.65水利论文!w!E!{N Y

12水利论文)n#D u q^}!N1n h

1972水利论文*L#_c1l lFt"K_ j

40700

.x#U%}5y/P0[D:Kd}*t0

561700水利论文TQdkz)Eq

428.8

A&F!K"P2tkd0

107水利论文w'I-kZ.q

1.91水利论文(tIMA4Jq0Yx

13水利论文+~"]5jk2Y#|7w&{ ez-RC

1973

-d4M1W+p)b|#^0

77000水利论文x*f"E&K4xQF

563600

C-tSs1h5D4sG0

482.1

?9_.i&Hd\/`$C0

140水利论文m n#B t%Gsrpr+O

3.00

0DG!NA)_4?$iJ0

14水利论文8}-_%AOU7_

1974

2E0mL C6u:m4_*A5nyqF0

72900

&H1H a O7r Wp0

557600水利论文`,\8JK(j e"I+Xk| P

415.1

6sc Tk;_:t0

121

z#[j\DUJs0

1.31

h!~(h hL9y+Wk1|I0

15

/K,g G|ENh0

1975水利论文e!B/qK5q A|

63300

#~ p$p0`8iF_0

528300水利论文D4^ eub^m|

536.7

"H-dZ$OE0

144

8F ?$Dt x b0^ r%? di.q0

2.33水利论文8M x"z*y'tqh

16水利论文:zJ4IM ]

1976水利论文P&| [2X;g8r;Y8c5or g

51600水利论文wb p(I k

488940

[bKd WV'lR0

385.1

2j [ [$g+{}lCe0

154

k*hCE/|*S+^(E)T0

3.55水利论文d2M9f;Bhyw%]$m

17

~oVbm MdG,` L0

1977水利论文U\;\w1R;Xl8K

60000水利论文 s)~M E-o

480500水利论文/VxY3q!rCP y-d*g

412.2

$r;c+^C Y0

111

3d M7vK+T??'_0

3.35

/L s1X2m'StNL0

18水利论文Ph2@+G0J%] s2?({7G

1978

6B AI&x9T,vIpgeU0

70000

%M Vq UJ;Fb0

530500水利论文ofA BE5}

567.1水利论文!S~#IY:i!RLIE,q

139水利论文,RK8a&uP Q9HOP |Q

2.55

c"iVfz0

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

b(H#ou0eJ#[-[s3X0

Results of fitting verification of yearly average concentration by using PPR model水利论文9FuC _&x4v


实测值水利论文am(jF.a jU0C

计算值水利论文V |2qK O

绝对误差

Q%[K0s D8Fw4V0

相对误差(%)水利论文i%{*z1n$Z9X6n

实测值

\ Z+B*DdP]+a8V0

计算值

`d w R-~l4vh0

绝对误差

"G{L;D9h~|}0

相对误差(%)水利论文%i)gR4q#g1^


3.40水利论文e*fp'H~ D-aWU

3.131

3U+lP6k.H W0

-0.087水利论文5C8ta5EE0l

-2.5水利论文.J2lyt}*nNO3A;h

2.72

N [ GU l h0

2.934水利论文zJW!e7Ut

0.214

WJ:X&b8mY}2U0

7.9水利论文 W.R5a:kd:@@D

3.00水利论文w+_.A I?(Y5Q+E

2.738

5?Q.R"zI.u1e0

-0.262

v/c k(@"Y1a0h[F0

-8.7

QF[!G"Ij a0

4.05水利论文 D`7gQ!? F,z6L

3.996水利论文c6XG[Im h^ Z

-0.054

v*X k'dN-i-R0

-1.3

!Cz*c1d#WvW"T0

2.70水利论文YI_eW2F3^1N8O

2.659

t4iSn2_#p0

-0.041水利论文J V*j+X^8p,C6t3j

-1.5水利论文yt3D P!waB

3.22水利论文O'D#K4Oo4gP y

3.214水利论文:q;u D}Y$g8C/I,i

-0.006水利论文_F.L }5mf

-0.2

o c(W| d0F0

3.37

Z g'o6DuyC:Z0

3.400

y$DqppCl0G0

0.030

C+H"]HGN)|Y0

0.9水利论文C'}] _7L m

2.65水利论文HV8h e;Ih4N.d#}R

2.628

lH8q*~OW0

0.022水利论文a$}Rp4fkc

-0.8

K-rw%TC,o;p2W0

6.894

Eg? aY6J9iQ)I0

6.926

_%|Fo9C(J(V$e0

0.032

m[4a x7r)v/j0

0.5水利论文+@\_D ?j h'G e

1.91

N/ebr?Kf~+X0

2.135

7V-VV(uh$J0

0.225水利论文OtEn B

11.8

9@+Z&fP%r kED0

5.03水利论文iz T+I|,rI's{ }#y

4.729水利论文.[o#E5fIt

-0.301水利论文xSsl5s

-6.0水利论文1] RQp_E)`m

3.00水利论文3u0W9E Yf

2.987

D1x8W)Ka0

-0.013

WNy5Kmo0

-0.4

{+p IW1LIk0

3.55水利论文)i+?i'k#cJ

3.826

(EG$G*w?R0

0.276

@y~ B4UV0

7.8水利论文 D)KJE5f'I7i Z;S

1.31水利论文 V/?`6~R6Z7B!}#Z Y

1.319

}X;Q g N \K9d? V%{0

0.009

1G7n4Zh%o cr8O0

0.7

[)\B`c.h&L0

2.33水利论文 o6z P7w'{g R d\(W9eb

2.461

1j4a0C"Mu8Fv@1{0

0.131水利论文 Yl1wP7f;W5t

5.6水利论文 zZ2b| Ra.W

3.35水利论文Tv?m&{&]v

2.780

!z/s0?9E'[7ma0

-0.570

ID?td0X{ l)f0

-17.0水利论文F-i8[|k H

3.55

Te7K-?;@$L"}9E0

2.869水利论文"A'N8Xlv f

-0.681水利论文M'Y.zsZ&u6](my

-19.2水利论文s5Bb b(L2n3v l

2.55

d&c@"EV\!x |o2H0

4.305

P8P!m { e3Ph$i-g2v0

1.755水利论文 ?#v5y3a"g

68.8

;j!`)HU'W+fPd6o$M0

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

表3 各因子的相对权值水利论文kn$]{*m~4G%^G

Relative weighted values of factors

V$W'`Q&E a?4P"Z0

权序水利论文[nOwt-I

1水利论文}7g-^*U f

2

yaN;s(OXx0

3水利论文5V L?rU GN([

4

L(zm HNM0

因子水利论文0F2E/?|e RR

X4

#p(yF1]-qm0

X3

dnJT]0

X1水利论文HR}:Uj |;S-e

X2

'^v1]3\a}V0

相对权值

0_W,JOJK7jTi8r0

1.00000

j!t'p3]2_*gu0

0.85190

.fom1t]}"L/] x UD0

0.71617水利论文%X7vUq$aEy

0.29765

7gu1| r,D%Y9s0

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

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

4R$IS q-v@#R0

4 结论水利论文4Mw(S-_ R/`

990112t1.gif (1658 bytes)

E8_E6pX0

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

Q7tyyaiDB0

Effects of fitting and verification of sediment yield

z%rb&dwk0

in a watershed by using PPR model水利论文/BQ$s%uRJU0p;}

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

;wZq0pUQ@/kq0

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

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

lO*@lf8x#F5I0

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

[.W&\B&og g ]9m0

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

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

参 考 文 献水利论文C\+Ji%Y3G \"x h

1 陈浩。水土保持学报。1992,6(2):17-23.水利论文Yf$n7c4vy.c

2 张醒。泥沙研究。1991,(1):15-19.水利论文W j+Dm*da6h+m

3 Friedman J.H.&. Stuetzle W.,J. Amer. Statis. Assoc. 1981, 76:817.水利论文 w(S%Fwk

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

{"c GTWLr2k$a0

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

`8q1E N~%gF0

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

&i| x&l'cr| `-P'd#V0

 水利论文t{8\ n-r:O4k

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

广告投放

广告投放