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

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

流域年均含沙量的PP回归预测水利论文Z$EB+kXl K

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

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

9a"w V{ r6n S#[%L0

关键词流域 输沙量 投影寻踪回归 预测水利论文ZnC-y8s%I|CU0k

1 引言

$Gw6`Aa:o0

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

-Knn1GZ\+qa0

2 PPR原理及算法简介

(_/v;@tA0

     投影寻踪是国际统计界70年代兴起的高新技术,是应用数学,统计学和计算机技术的交叉学科,属前沿领域。水利论文7Va'}3h*R/R oW

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

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

990112e1.gif (1513 bytes)水利论文Z$q+O-b~3\ m~

(1)

y!O m1v*T7\@0

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

z2H,Pd#Cy&P b0

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

TC)u2wVx#Z0

990112e5.gif (1470 bytes)

U\,y y? {0

(2)

9^~*Q BNmtoa$|0

  它实际上是采用分层分组迭代交替优化方法对式(2)中的参数α,β,Mu和岭函数Gm寻优。实现步骤为水利论文9jQz0~x"z _'W

  ①给定一个初始模型;水利论文6@7bX%Pz@:|L

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

%|M4wEj9{ Ae/o0

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

}uM+p2k"F0

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

Ew4P kr0

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

990112e6.gif (1600 bytes)

qeWd)X0

(3)

[m$w?8g;X0

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

*V h9tvD0

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

*G }5E9U'~w&h0

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

9iH#|.],J:_/CG~0

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

     该流域的年均含沙量及有关的4个因子的历史监测数据见表1,数据取自文献[6]。水利论文Yx(T:Y%PU

     用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个因子对年均含沙量影响从大到小依次是年径流量>降雨量>采伐面积>采伐量。水利论文+Fhb'B$I/au

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

Kw*HC!^9d0n|:}x0

Yearly average sediment concentration and measured data of related factors

g?~0d2L.g$}1KZ0

序号水利论文IW O4~%ks'q

年份

-[1bKP!X%T2j0

采伐面积

_ujN"l6K"d0

X1(M2

_,j?HK:R_Q-h0

采伐量水利论文'wh u tg0e_

X2(M3

A~6\~0PH1m;B0

降雨量水利论文 ])l j6Mj y!To t

X3(mm)

Pf0aAShAw.L0

年平均流量

5n.[&muT0

X4(m3/s)

Wk2V?m1Q'R0

年平均含沙量水利论文$Li0Li:e}5T(]N[{

Y(kg/m3水利论文7pQ/Lm | O2t


1水利论文!t3EO5^UmR.P [?

1961水利论文$}Eiz{j7WS

15100水利论文-X3z&wU p o(LG

29000水利论文,Y rva&s

508.0

)R8?U)d'W+qYM1O0

180水利论文`#j4rN9}^o

3.40

#w"IPN2TB%Ou0

2

Ag2K+e({i0

1962

/Z6?]B`c b0

20500

`5T ].`4a.J p D1r0

43350

3BP!p'Z9NK,r0

453.7

WL7Y4m$??%O{X0

141

&u6|O L}xHE0

3.00

q7{I*u+|%HM0

3

Rlp%UO,h:DY0

1963水利论文4Y4~[N&i0c;}K(o

80000水利论文 V,n3R~}6gk

92610水利论文,VAr_!V-`@8X-m

487.9水利论文O gIbU3JX c)T ~

132

&y.vMsi.`6y F0

2.70水利论文^ Q8H^.od&O

4水利论文2U,wm3TlSb6q

1964水利论文%ba2S"},b"L%^

91500水利论文a1d#x_/M7B'?k

142775

xh"r)t [Lk6N'L0

572.3

f#^"E1\ D6@g0

182

"QWe s v XX9W0

3.37水利论文,q)F$~8Zt1jx"R

5水利论文 ~l R5Z&`c5r'TB

1965

W;E+C-wC0

82500

1p-j-V%t3P-Z7d\-b3Z0

2123160

1S2D*B1y7S:]0

455.7水利论文p oWF0cA

113

X Wqt9X/E.W7?8I2\1L0

6.894

]y9Y9x-Ah,n4z1D0

6

0q"{-[%\jqf~ T0

1966水利论文&].o:Z6P-B${dIM$x$c

20000水利论文2Xm1]4x{Vm

227800

Z:T*z/Xzp&s0

481.3

TY m|t OTk Z0

170水利论文cZ|.G([Q

5.03水利论文)}Spfz

7水利论文z6ep6e2?#Ky"B,n

1967

H@'_RW0

17800水利论文9u bE6| qc.T

140000

iPTk VI0

541.3

7g[-y(J4]0

199水利论文] V^ h2r"q_

3.55水利论文.XZ ](MtZ6ZO

8水利论文%C]K(jF

1968水利论文 nLZ[8g@m`

3900水利论文)F0Y!l0{}&K&?p*Zy

15980

Z} ue(R#f8B?L-N0

538.6

}$fZ$R6S:]$R0

186

s'A%M;K|i T0

2.72水利论文qWo)j5oZ3Z5l3p,N,F

9

oStc S,j(v|l9i0

1969

#R L!xLx c$Z@u0

17300

,Cty(l(o3R7~p0

223200

(k EH D4o3h:L ]0

460.6

*~!cz1~CEg-d+k0

99.8

(Xf!V#c5V?0

4.05水利论文jc&E$e/C?K

10

o {@ \ Y0

1970水利论文*I:Br \0qZ-fJ

25700水利论文6ur&W-um

229400

DMI"}Vo(`9i0

393.1

0rqqPYO Y0

133

)a Ag2ng9? B9yn:{0

3.22水利论文4iv-O f9DQ

11

)|Q1} UvA0l%G0

1971水利论文s*h/f w6P`9M

49400

/Bee Y"`L.wb0

424500水利论文0sZ6n-v/[9Z;I&?I

373.9水利论文DEr$p4Q5U0GR

106

/fb9dE5T.Ym-g0

2.65水利论文Z&I0}:p!Pu7Xw

12

,K.O3Y ?lttF0

1972水利论文 kyW"{8D @:`

40700

\9nZ1O\+QA0

561700

*cC.qdX(tY%y%y0

428.8水利论文6N*zy3Fq!\8g$K

107

J#aM6Mp!g M0

1.91

#}.O+I8Y-C a"^@Ni0

13水利论文G3RP[4`%G"Q u*b

1973水利论文*AZU/F8o

77000

u|iKY%t7K2kD_0

563600

,E-`eh7V0

482.1水利论文#ZZ1\x*[c+_:J8I

140水利论文TU)n:I0~:d!Z3K6p

3.00水利论文},|9|(}NW*~$N

14

[8Q!XU0}&a2y0

1974水利论文}'R,f0jcdZ

72900

5|-b*xltXQ0

557600

g LE2`'[0

415.1

r)b]gE7c2V$[%br9~0

121水利论文'J2[y-S3e%L@o

1.31

Z^-@m&Q0

15

P.Y:? z~?Nv0

1975

4A!]9Ko2gY8u TW0

63300水利论文*SSeZ5[6r

528300

Pq.LR6Z/es0g0

536.7水利论文 tF pS c.N'e&M

144水利论文*ck(~nb@L9h*kx

2.33水利论文 a_mOak i

16水利论文rBcfi,Q

1976水利论文1_w#W-R6e5T:M

51600水利论文5[D6~aV n

488940水利论文%H-KD:c;U%I

385.1水利论文 wv-iR(w@}J"RQV

154

;JZs7Syw6Ct/Y0

3.55

MydE)[M1L5D0

17水利论文{l5t K/V;s_3X\p

1977水利论文HS%Dx$l%d

60000水利论文 ~%H.N'^;I P n

480500

.sq0Xz$iVut;?0

412.2

5w2W2biM;nADc;I;U0

111水利论文grrK*AC0j4A @

3.35

"S(LH3Q6? D?0

18水利论文,q/ElU1vf4@x

1978水利论文Ozt h$w ^Ga q e

70000

-W`s%{(oS[0

530500

?iD3G2m3w-?r&L0

567.1

Y m@8h W~2dLEJ0

139水利论文[,{)ij;C~C$c8O)u g

2.55

9C)~E} w`0

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

!lU/A:}VoO|0

Results of fitting verification of yearly average concentration by using PPR model水利论文(ew0IwV6g;@u-I


实测值水利论文 G+f(yaZPq+z

计算值

(C`a4Q ]]dRO0

绝对误差

/kzG0Yn0

相对误差(%)水利论文YO Z:`Q:B,R.l

实测值

d?&aVbUbg0

计算值

;g(icmKysdNg6[0

绝对误差

G9C3\4AK/\*Xn5E0

相对误差(%)

JEm*v2A;N\0F%W)vr0

3.40水利论文-FV*s4w2u/_6u d

3.131

+mF`w/lU0

-0.087水利论文 tCe:P6j

-2.5水利论文c%v` s3ye+q

2.72

H^+TNp4]-P2}0

2.934水利论文mMd#WX

0.214

+|x C,[+T:r]}}0

7.9

lp:@,h4Sk^ _0

3.00

tuf-u zH@2O$l$dy9r0

2.738

M5j Ix%k9[q0

-0.262

;o"qdAk{6@]0

-8.7

l!] G'D%P7j+}1P0d0

4.05

w w+p4A~+{ {T0

3.996水利论文u?n6K,uZ0y

-0.054水利论文~Qq ~{Z"j#OE

-1.3水利论文!Tg2R dD$zk0?T

2.70水利论文J_-YjQp

2.659水利论文#i&n4i?+Bb-W

-0.041水利论文3Km#h&S(w YwkL

-1.5

F ]b;]o0

3.22水利论文 y#]%^'EV't&^2t

3.214水利论文$W1thE]"@

-0.006水利论文.[\1i5F5Q#b7V

-0.2

(@#?(FDNIdBWG0

3.37水利论文L!C3\m fL

3.400

&aEa'\,X_4Ex0

0.030

_!v3u8|D4cT5n}$w0

0.9水利论文t\7|eu ~{t.q8F*B

2.65

ON T }-\C2@9aOT0

2.628水利论文3S iI8r.KN ~

0.022

3rue*nM)epF0

-0.8

H_:UO4Q)Fwz8g0

6.894水利论文 oLy3pNdFS|

6.926水利论文{6|mP1Y#BN&o+|

0.032水利论文*yz.RV'`uO-G

0.5

2E3oPk`;`a3v u0

1.91水利论文QAqq^9^

2.135

Z.\,Qb3hO8c0

0.225

4A;QKk+mi~&W0

11.8水利论文a4| Y m(xB Q

5.03水利论文0}$IH*Y,H]![(u

4.729

C(AI\-zQ)A0

-0.301水利论文8?.BB A#|e

-6.0水利论文.MkDO$[Sv^"L^

3.00

KN8jR#i v0

2.987水利论文tSE%M4o*_9s

-0.013水利论文3u:u3us1ZG-{(z

-0.4水利论文-d,Zi j7e/X

3.55

:{+s's4U&_$u(w0

3.826

;_w9x T+dv0

0.276

Xu.B2y`0

7.8水利论文Z [X{(o%e9^ h

1.31水利论文%A}(|+]t5]-u

1.319

}j!\Q.P*u6J1vG9Y0

0.009

i gw#o3{A&j X0

0.7水利论文 g:PB2?^I

2.33

hMn mu3rW^0

2.461

^ui9U4Y2RJ(KUE0

0.131

&s5p5kU Y*ISU'A:l!~0

5.6

0Ba$b.a5\0

3.35水利论文a!v-l'?)R!hd-B

2.780

K#[[%w-x/^a0

-0.570

vF[4Pc%Lm6W0

-17.0

b}{/\Zq'h:]YK0

3.55

U,W/`dE?%f0

2.869水利论文8XiHm \HN]

-0.681水利论文zwm"Bd6{

-19.2

K'ZG&G'Y j'Fd0

2.55水利论文*K s!f4SO.E3E[/i

4.305水利论文&S!J$h9Tg4} ?)m!Of]

1.755

D y~o2l#^UL;|3Vf0

68.8

.U cpx9^%@0

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

表3 各因子的相对权值水利论文3dO$hqm9rC7b!mg%Uo

Relative weighted values of factors

@um y5Qtd6O0

权序

h!a9vb-]$fw0

1

[/eyekcD ~ R7l,u \0

2水利论文C8Y z-z6o1B L

3

b`8t,~QdN3j0

4

zu3@7zcN4{0

因子

k ^ C$R6J,BR~p0

X4

P.MS5t(~B0

X3水利论文.F6o(X)u"T],s

X1水利论文 ^4P2mpV6N|M

X2水利论文SDlF9B!Z,l0}

相对权值

G @6x x:Y,B7_2^%}1A0

1.00000

.I9lf'ck1A0

0.85190水利论文O D9s?$R z+k^:fb+K*{

0.71617

RS `[yY0d:h&x5h0

0.29765

(Rj{Q1aw!_b!b0

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

A+ng{E'|mk0

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

!ER5yu \3ZA,H0

4 结论

+n|V$]q"m+K&u3i1X0

990112t1.gif (1658 bytes)水利论文8{ N5V9H&T.};N

图1 某流域产沙量的PP回归拟合及预留检验效果水利论文J8^ pYnd)Y|?

Effects of fitting and verification of sediment yield

/f!jUXDvVO0

in a watershed by using PPR model

%C)l@4Z2h M lD0

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

?1g)jb4h-y j0

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

MM.]x6BY%I&G[4RH0

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

     4.PP回归也可用于水土流失中其它指标的预测建模。水利论文%Rm%e c?0E

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

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

参 考 文 献

R:VS f1z0

1 陈浩。水土保持学报。1992,6(2):17-23.水利论文n.UY8GS.rAi5Q3M!q/pz

2 张醒。泥沙研究。1991,(1):15-19.水利论文5S]JcTt%U#p4`

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

KTA Mj`!O0

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

5F a^:[K/}0

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

_%A4hq&[ d9t0

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

;l'c{w4ll0

 水利论文#h'M~Q@^&G5]

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

广告投放

广告投放