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

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

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

%e5q9\rDCT n4T}0
李祚泳 邓新民侯宇光
(成都气象学院)(四川联合大学)

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

F#@6[%c~9_$S%l+|l2v0

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

`b?N? B:gD$I0

1 引言水利论文`HX*X2?!sw

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

2 PPR原理及算法简介水利论文P)D:aX;N#lj.],ns

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

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

Wp'V3Q N1c7M g0

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

~8j sy ye;rG0

990112e1.gif (1513 bytes)水利论文!k8KW"_i.R

(1)水利论文:ZnTQ}n.|

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

Friedman和Stuetzle提出了实现PPR的SMART多重平滑回归技术,SMART模型具有如下形式水利论文-c9aG v!b5PsH8G|

990112e5.gif (1470 bytes)水利论文 nfG5l+u"_ ~@

(2)

r"l k#a#l0

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

e2{'c%I!{1ie5`b0

  ①给定一个初始模型;水利论文^:r#_q9V c0ZN

  ②把数据投影到一个低维空间上,找出数据与现有模型相差最大的投影,这就表明在这个投影中含有现有模型中没有反映的结构;水利论文'n@-`4X4Pn&@

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

-^ X3c*vW'_bt"jO*c0

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

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

v$GV/ts[~0

990112e6.gif (1600 bytes)水利论文w9L5t_oF{

(3)水利论文9k8{"f\ {iJ+z

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

q'G?#]e%my'`/u0

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

n*[!oz;js,T0

3 某流域的年均含沙量的PPR回归预测水利论文 yf3huz]4s

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

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

x1y1| _.j0

     用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个因子对年均含沙量影响从大到小依次是年径流量>降雨量>采伐面积>采伐量。水利论文/^0[ y(|#S9J4oY

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

"M%Np#RcXB0

Yearly average sediment concentration and measured data of related factors水利论文/W3@Po#{g6]0j


序号

2e)O,U{6c8h3d0

年份水利论文O^aUh d

采伐面积水利论文nP(G%sZA q\/o L

X1(M2

9[Z:} C qXaCm0

采伐量水利论文!wt+z(h+q\

X2(M3水利论文/A7E0l$zm)E#}

降雨量水利论文/A1UP Ru@Q+V

X3(mm)

O1gCT&S9^+}6kfnt.r0

年平均流量水利论文(V&P:H8T)g&zM)]

X4(m3/s)水利论文nsVD GW

年平均含沙量

;a8q8g`F.\0

Y(kg/m3水利论文2_,U)n mg


1

.T8a6U ]2i6{0

1961水利论文 V.Kp bKEiB }V+r

15100水利论文,|~&} UQ

29000

I3^'A;Xt"h0

508.0

VO0DVN#y9B?0

180水利论文#b v ~&q2v^g-fc

3.40

7h0DNt)odr0

2水利论文cnqqc|;E j

1962

^a;RIP,B"A1C0

20500水利论文+hw s+zU:v3?/{-y

43350水利论文zH%c }8M q4{

453.7水利论文E:Ss#quV#}

141

8X s j @0{,@*E/?@0

3.00水利论文 M'Eh{i!r,^;W }/L

3水利论文Ga-pMF4F$uCa

1963

v/_BB)G!Z_0

80000水利论文oz1BpPR2d

92610水利论文#hC%i1?*U~*i8Z

487.9

jOeE2@3s$K1M*U0

132水利论文 L.XM*R)D8Bv } b0vP

2.70水利论文*m.z4L7j%T6f

4

1F)z6_NK.oP0

1964水利论文6Po.m%iP y

91500水利论文 IVo n;}u

142775

*YOq,HQs D9^%hE!G0

572.3水利论文 n*|"|STbUVw

182水利论文-Y ic8[Qt.EU

3.37水利论文+Xb aVG7hNhS

5

#g{2Xe l M0c0

1965

,t*_Xt+Nci0

82500水利论文H e0L4i5t{6d\ t

2123160水利论文3AI5G4W vyZ ?m5G R a

455.7水利论文YX;mQHG.aw

113水利论文C'Kq'N7Op;dKU

6.894水利论文8?"n5[ri

6

!Q V fG-f r6We%H0

1966水利论文*Gj!`q{@[J {

20000

].wyy T0

227800

o?d9GvE2|0

481.3

!_zA9vv#Pt0

170

q?/K-|`2mi0

5.03

6On-q,P,H;u0

7水利论文$Z!G*y7h1D#L3g

1967水利论文XpRG^o ~ t8J

17800

o+B _ w{h^5x0

140000水利论文`p4U5L7}NN

541.3水利论文 U|4AQ N+? I(e#B

199

i b/o Y5D:K\X a-Y;g0

3.55

JD#q b't/u G*Y%[;O0

8水利论文0i'i1r]M@QlZ

1968

LNUT$Ghr*Mb!}0

3900水利论文6cFRb8_FFK

15980

@9dW|b?:o.}\l0

538.6水利论文"VoW'o R#c~]:^!{

186

rL\ W Ky[0

2.72

{W;C\]:?#fo0

9水利论文DXE$BevW*}4A

1969水利论文O;V"f1PL

17300

R a*]'r9Us,b-X0

223200

4N U6a8C3rG Q1H \~-|0

460.6水利论文1pi*ER1J*X#e#FE!i#~

99.8

T7cTH4T:~0M0

4.05

6bLPpgz0

10

g(vJ;B]@0

1970

`s6l&r?)a0

25700水利论文}8Bd/Z(R7gtn-?

229400水利论文I4sB4c9iR;z#|

393.1水利论文t/{)T j ht6G/n5q

133

E&Z$h;wm~?0

3.22

`L}-\9u0

11水利论文 cF0A ?5f!\

1971水利论文3OQk%s&U zK4n9mY

49400

4n)Y c [ez^0

424500

T O1D.@f0

373.9水利论文1wCt,|!L_

106

Ixr'Uv5lVjw0

2.65水利论文gv&r'h9]b

12水利论文4Z.V)`Py5`

1972水利论文7]'\U;T#UEagA5X

40700

r)k+x*`l(e0

561700

1q$^[:@t7~p0

428.8

h1h A3C&^ al"^0

107

'dZYcQ#~0

1.91

4\#T)VK9U(p2Re0

13水利论文b,Zxw,y`$r+w5W

1973

F^$|(?/OM,}0

77000水利论文H w o-x"s1y}

563600水利论文gT1o']4S{] R3f

482.1水利论文 j7F7Q6g/h

140水利论文/o3I&kyi2O9S6?

3.00水利论文 mCH3X'uQlUV

14

-Tghs\g,r O6R0

1974

+?+n$h*WDN0

72900

{:_f"l_Uj0

557600水利论文;?j5]:eg.@4g

415.1水利论文-I~k0|-J4\

121水利论文:yqrl(g7Do9Q

1.31水利论文#h7X%MZ1IU,zI9w8fL

15

uW9FnNW2])G0

1975水利论文/dZn0}?*_Z

63300水利论文g NT8Vf'@%h

528300

/I4^YG$Kf)rd\!X[0

536.7水利论文E:}Oy2M$[%|]+[

144

/MC2V.aq Z~Y0

2.33水利论文1B'OTX*ud

16水利论文MH-H8UdPm

1976水利论文.\ y6cI2I*y@]-ZD

51600水利论文$~ L:@"wj:}K5v;g3^

488940水利论文*v_$L;W.SY

385.1水利论文.ij }'K(GrfZWC"G

154水利论文y0g8W3NI)V's"C

3.55水利论文H iS(oN2~t7} mr

17

)k!M6xs N&K0

1977

] [`4E'Hf0

60000水利论文'ZF _:j$@8P{ sHm P

480500

G$Tmn/nQ(y3I{0

412.2水利论文"k,J&p @/XN`C ?8Z

111水利论文p^Avtv \sI

3.35

j8G i:c:N D%t} J0

18

~&Up lFs0

1978

^e;ko G})S0

70000水利论文-P6C AZ2nDO

530500

\vPmF$E0

567.1

f\e-U(LUT%\1O0

139

0U^'C I_0

2.55

M+P~be#pl0

表2 年平均含沙量PPR模型拟合和预留检验结果水利论文q"s Ne;@p

Results of fitting verification of yearly average concentration by using PPR model水利论文/oV]Ra|_


实测值水利论文y`^#l7Mk/G-z

计算值

v@7|$Be } v0

绝对误差水利论文-wm bm"F)[%t

相对误差(%)

cC"Mo&{@E9k0

实测值水利论文b+] c1N:S*r ]

计算值

7p#X$g%Q6TGX0

绝对误差水利论文,D5t{|9Ku

相对误差(%)水利论文;N"g G}^:~NT


3.40水利论文N5[3\q bN b

3.131

2G0la&R3|e ?J0

-0.087水利论文 e F&~6`qP6d3}@

-2.5水利论文+xRJ n R)s8J:h7v U7o1B

2.72水利论文2@yy7Hn }4l

2.934

sq%?Z4_0

0.214水利论文"`A4H,[r*M&G!K

7.9水利论文7GMm.BH:f nK

3.00

RT;b7|a'X%tg8m;i0

2.738水利论文6[?reC2Y"eAy

-0.262水利论文R^ D(i:z `%Im*y

-8.7

4D3~2gdWP'\UH8H0

4.05

J7yK6yW9Z*hH0

3.996

c `$`+b'B0

-0.054

bl^8@dkn0

-1.3

wPzgw;d7m7n0

2.70水利论文/`%I h"\5G

2.659水利论文#Y"@sZN8HkL

-0.041

1x/Tw,QM'm9`#st0

-1.5水利论文!Lp k OE {;\GLWz&q

3.22水利论文l Ar:Q'~ ltyc

3.214

8e {4wz%bb b0

-0.006水利论文+ohG#H;[NA

-0.2

/O-f6@"d*cAU~0

3.37

+[cI7[0q0

3.400

!x3t2D1Ip jpF0

0.030水利论文 t(SR x nKS(qi

0.9

(~f0] k#H?;V0

2.65水利论文.d"~$DSADgLwQ

2.628水利论文 G$G,L1VOM9V4|

0.022

"iXE&U)\#jpqg0

-0.8

5pE.So`,F]3I0

6.894水利论文'p#jk7YXw4W

6.926水利论文ze)N9EBX r {[

0.032水利论文7K K-UM)g8|

0.5

c)dj;X^%ChkI5A0

1.91水利论文 ?0SF)_8l"]?o

2.135

6A Rj0k2xA{ p~8X0

0.225

#p^*D%i'W-]0

11.8水利论文i6tJSY&g9rB

5.03水利论文2dC}(N_+w0mJ

4.729

)cN |P ]-q i#G7c0

-0.301水利论文I l5E S6T,P V`

-6.0水利论文#Wy*C.o4bR8q

3.00

|:z M x,Bkp~`8PX0

2.987水利论文/`DMf-v"e/C R l

-0.013

7}d0z#r X+d0

-0.4水利论文C*N*s5fH:I3m)@|

3.55水利论文 K^ H2Ra2k

3.826

y#x-K~ jaG%s0

0.276水利论文N;{~$R NL3nhE

7.8水利论文O/^%iJp+y|a(j&~ d!F

1.31水利论文1e T0I}o1Lv5S3{

1.319

8EA'_p%L,VO,i0c&U0

0.009水利论文ZOK*S,v a@H3i ?

0.7水利论文y3XAI9Q/e

2.33

i0F he!BcLlBm0

2.461

A~2fr]2n-T$l0

0.131

/C$va"_Q5Ce5O0

5.6

m tt6IrL0

3.35

2F*{*y!_ } ]6Q6~0

2.780水利论文.j%W.D4n$S3h

-0.570水利论文h&W;d G2jX-Z

-17.0

k8no$fZ0

3.55水利论文HN$UmA w0O sb

2.869水利论文iPw-M#U'q:x

-0.681水利论文m/jQ Qo.O

-19.2

:bKLaw:G0

2.55水利论文.tYY%e y3O^F

4.305水利论文+G#z1B ]V(G/H

1.755

?[$VM;H/Z%iB0

68.8水利论文:}t6{~!T


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

表3 各因子的相对权值

W;w @)D5\}[0

Relative weighted values of factors水利论文6}6V^{m_0ZB


权序

g} eTfz&l0

1水利论文!ac&@l%_MNI4r][

2水利论文_P{$b8VZU_o

3水利论文$l})J7t3D.K_m

4

s2i&W(Gp/Ffr0

因子水利论文Bsb#@7r_#M2KV(\;Sv LE

X4

+K[u6Qkzk0m0z%Z0

X3

K5P]u1U0

X1水利论文|*H{AY'n9BV&f'Z'^

X2水利论文7~/^$Z/`u s4d F

相对权值

#nrEC6e/{ G4xN0

1.00000

p[A8L C9|WY0

0.85190

|7K0TTjgy4P&?0

0.71617水利论文e?aK5je[-~b?n

0.29765水利论文~'e^7k#pq T3C


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

I%Q8f*P~#n \0

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

4 结论水利论文6AY @-E \1IMG

990112t1.gif (1658 bytes)水利论文'GNOgE\:Y6A!R

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

%JE qQ\*KjJ*R}0

Effects of fitting and verification of sediment yield

mt(}Z u0

in a watershed by using PPR model水利论文ie6G h+UU'b:rUv:|

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

IN.CK!rtP6Z&p$F0

     2.PP回归无论对于正态或偏态分布,线性或非线性,独立或非独立分布的数据都能有效地处理。因此,它适用于分析和处理流域含沙量等非正态和非线性等有关问题。水利论文)[3b}4~(m:j&i/ub

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

     4.PP回归也可用于水土流失中其它指标的预测建模。水利论文4KD.\^J]rK}

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

'VR'G"a Eq`VKV B0

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

参 考 文 献水利论文eB$dKA!pOZ

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

#P:X9J-F;I5g#I'n7O0

2 张醒。泥沙研究。1991,(1):15-19.水利论文fo:V*K[$S+jL

3 Friedman J.H.&. Stuetzle W.,J. Amer. Statis. Assoc. 1981, 76:817.水利论文u4x.T;bdeA

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

n'I| P,r!?"G0

5 郑祖国。水文。1994,(4):6-10.水利论文 nq%N{{GO

6 吴建明编著。环境统计学。北京:中国环境科学出版社,1991:504.水利论文_;?-X*W5x _z!D

 

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

广告投放

广告投放