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

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

流域年均含沙量的PP回归预测水利论文UBE9IR1{8`g

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

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

关键词流域 输沙量 投影寻踪回归 预测水利论文*u i;M$Y9~7c E5C

1 引言水利论文$^5Q2Bh/O9r

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

M(dLu x(x#}\3n0

2 PPR原理及算法简介水利论文hP/e9d5T

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

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

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

FLXF WBZa+qL-@0

990112e1.gif (1513 bytes)水利论文;M&Qyw0t.@

(1)水利论文&o8M9b&Xx D a

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

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

#w!n4f3s+F"pI0

990112e5.gif (1470 bytes)

e GRH xf7^)o6\0

(2)

1p&O/~kTq8j8~8]0

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

l7z&iD O;rR Z8kZ0

  ①给定一个初始模型;水利论文$d(z ?%a.Vk8yd&E

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

8]B~R_U;o#S0

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

-bB.hn$Blp7]0

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

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

4IO eWn]0

990112e6.gif (1600 bytes)

:xg&q^/R,[~0

(3)水利论文;] q4} ~ J+T;H,UV

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

     用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]。水利论文1G2y5b2K9G

3 某流域的年均含沙量的PPR回归预测水利论文b$_ q1pORl+WS!U

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

*U qIjMe0

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

C$Z5w,p"w'R0

     用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个因子对年均含沙量影响从大到小依次是年径流量>降雨量>采伐面积>采伐量。水利论文Yv'O1QM(\{

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

*gx4W{H*^#w'a0

Yearly average sediment concentration and measured data of related factors

"Y;s Ua}0

序号水利论文x7A v^@9v#\D

年份

` Y5jTEX0

采伐面积水利论文3k2s?:i#\ @ Ru#H

X1(M2水利论文l.uI R Xb ??#G

采伐量水利论文T1Q;@ q(X's

X2(M3水利论文@r S1]$K5a'L1w

降雨量水利论文HE@K~:T [D"j1N!^

X3(mm)水利论文H*~9lf9n&B6UI@

年平均流量

a$OX8k Q6C)zF0

X4(m3/s)

,gho)uH aj0

年平均含沙量水利论文6HTQ4oo+T5C(Dp

Y(kg/m3水利论文a%k(I!W*e-C4Y a


1

Pql8_ I$_0

1961水利论文$sn'qKWRD5U2I

15100水利论文!iG-~AR`7n!i:E

29000水利论文1d XYn a e&J!g s3W5^

508.0

7L@}\%y0

180水利论文#LB~:or!z

3.40水利论文eZTU9W;]YU.]

2

WdX-s\;O/I3[c&q0

1962

,|l8_tBb_0

20500水利论文2I.x5J,_)XH-Yf

43350水利论文!ytK(`2iD

453.7

9LH }%Wl$M_0

141水利论文 V?C^|b tyK!T~

3.00

}-M.r:_R+Oc`0

3水利论文&N-L%te5jy,[

1963

,K9^+fa*If0

80000

sv!Z i/s_ hXC0

92610

*pvd8~u,`&r0

487.9

RY A8_#FM6{T#J0

132

'_mS(?Zy0

2.70

:L`c }U_}0

4

du`aeU \? E0

1964水利论文 h8`g4y[V(W_ m

91500水利论文&dI0_(rc8?

142775水利论文4o8ZO1np z][#P

572.3水利论文T wwk,n"T

182水利论文6S2fQ _W$sE$?:R

3.37水利论文G O1x!MTc WrD ]

5水利论文G-BoO#h `l;X

1965水利论文o"Ou6r"X)HJG+i

82500

#hN-q3C$]k1i!R0

2123160

1g,R+x4Y6sK]L-p@5]0

455.7水利论文hv4L8v]0k/x

113

J~^![/Rd0

6.894水利论文 o(t.[zs

6

EX3p4_ t0

1966水利论文/Y6^b!f6ep9\6L.N

20000

Pb,i8Jc0

227800水利论文L:VyC_"y

481.3

8k"j*i fSa0R0

170水利论文$lel3y w|8A

5.03水利论文4fK|-VU5?i5[

7

K'X0iw2gBa2{+U0

1967水利论文:`WI6I kEV([o3^-T

17800水利论文5V$J;IT`q G

140000

#G'i s/YR3T*d%F0

541.3

&g5x#Y'grM&_j]0

199水利论文(U'EZ[ A kX,} J

3.55水利论文6B\ Hy.Wk w4A C3[

8水利论文#u6F(}.@+{ \8A5Tu

1968

{h:b"a6y0

3900水利论文4CfGU2u2OV.~

15980

\(l i&n`P0

538.6水利论文/s'|*n"U9a-SB

186

0_] s4@*l0

2.72

%F]D1q[5T0

9

"GF)}#YN ANA(E%p$n-?0

1969

o1?x&]%LZUQ0

17300水利论文aT N$@6IJM"tD

223200

R|"B-X/{H3T ^"d0

460.6水利论文;dx1IL?!a

99.8

,nI VfV?Xc)b0

4.05

yU5Ky4X|0

10

H:\#Gqb:\ T:^(^0

1970水利论文S[U*^*hQo|T&A%L

25700

"]c.I0`C?3pF0

229400水利论文 Jq,?'cX Duq

393.1水利论文*BgPd8l1nK0P.B

133

%V(lbQH)\)c,HN0

3.22

/u)syw@0r2G0

11水利论文(G:P {%fG[#v

1971水利论文z[2Z/S9o hW

49400水利论文-r {Z7dZ6]9D

424500水利论文#["y:NXuV

373.9水利论文p o(XCPNw3T{g"J

106

'_"FO/P$Ak0

2.65

7X)Du6K8k*~;K0

12

Eav @4d}2b.o0

1972水利论文6G BCp$i

40700水利论文GaBvt-aObk(L-];J

561700水利论文6n"O_m(f+K

428.8

)aw&j:NN&bl0

107

)mE;Vz Mg0

1.91

qz2snr [(h#a*Q2W0

13

5Hb5^ BKOE0

1973水利论文GT`np!L,B e

77000水利论文 GO3W6r^5zm{Qt

563600

!P"BQ#EE#a _+N0

482.1

tK S/pSv5^M(H0

140水利论文Y#D)B!Gc1F#B3q,\?3a

3.00

w+RB/SY9m0

14水利论文 W3xI9Sp&B{&W*t.f7k

1974水利论文N,M A4e7u D?/a

72900

-y ku K[0

557600

*q%m;Z3z Y j5B3|:uS0

415.1水利论文3oLYqo6A

121

'k^ PD&l9p7n&}:V#nY0

1.31

:F r wb s't,DO.R0

15水利论文 fVs ?&QG(L

1975

`0J'rV g,X0

63300水利论文R"U;Pf/F(N;M

528300水利论文.D y2vTq L

536.7

*y:x!e e"Gl"dQK0

144

A#H4lf:vNXU-r9d0

2.33

XW| c2e5p0

16

4M)e0aBu8Bw$q0

1976

0q8s&e:K'zu7`$?W^G0

51600

$B?]4Q w?C,z0

488940

cr3A*B4M*Vq b9t-y0

385.1水利论文\P|!\#C4p.W

154

U c"ED k6U9D0

3.55

{Q T:@5J0w8|O0

17水利论文;M0T+~,ZP"n

1977水利论文,|,Jo b*Fg4@!}2i

60000

.FT{ gIlU Gr3L0

480500

n+k8l P%@1k5{0

412.2

4VZLN^R\0

111水利论文 Z(]S|3E%f`

3.35水利论文2bV7tw6mJ

18

eXd1?cL mt_)}0

1978

t4d"IzZ+lb0

70000

%Y{;r;_J-} a~{)?0

530500水利论文 J*A(U(xv c e;{

567.1水利论文@3H9])vb0B)D"W

139水利论文o cf,e,UTi

2.55

3H}o!Bn,u0

表2 年平均含沙量PPR模型拟合和预留检验结果水利论文T9fkDso4z-|;t ? I[

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

I L.M\4Uv,t?0

实测值

}$IQ v_.e0VDp'J0

计算值

P3yV-u P)v"FZ)F(M0

绝对误差水利论文 m r7N `7@GQV9z$OY

相对误差(%)

_fy hq0

实测值水利论文K&yZ[:P#c;Z{

计算值

Y5p,M;]:Vr0

绝对误差

n/mD H v&y l-z;d0

相对误差(%)

l&yGfIP _7\6m%s0

3.40水利论文3R3wU%i%~-f

3.131水利论文0t!fb%A%E c R|

-0.087

%I'O^hT?0

-2.5

vo5R xom0

2.72

]bIV1g.YrG0

2.934水利论文V\ S}`F7R T d

0.214水利论文/cF1ZC a)}

7.9水利论文S:g h0y3X e-xzV_ a

3.00

-dr wETR5hei2M0

2.738

E7V D)a[0

-0.262水利论文7U(KO9\3{ l'Rpm

-8.7

u7q:l ev{0

4.05

-F%f)W {#d"Xp n0

3.996水利论文;c!p8K0d%vO,d @|P

-0.054水利论文.\:y |1M6\*f'U6PQ_

-1.3

-D \jy$w#J7_%t0

2.70水利论文.M%Y1\|}g!lK$S

2.659

ny].x7~z:U[hW0

-0.041

r%c5|GRh7B0

-1.5

^jsV~?;T0

3.22水利论文%_iX1Z ?@

3.214

,A4I"Uk7V0

-0.006水利论文B~(Ww$X&d"Bs

-0.2

P G`^A,|(`0

3.37水利论文#q bs(^mcJ[5_

3.400水利论文3D8[i}$XW;QY

0.030水利论文POx(J8A,b8b2m

0.9

6O9t#^:^{0

2.65

rcnv6K0

2.628水利论文4^7P-m+oGB8A.s+L!l

0.022

2\RY zts's1p0

-0.8

__|YEU;p0

6.894

_7G!M7tE7K"m}$_0

6.926

/c'U!l-F$I[W0

0.032水利论文AMF;qC4F [G]wI1U

0.5

|2{]yjw0

1.91

$Gd+c+Y3Ne]0

2.135水利论文n;yvA%|FO

0.225水利论文;H)p({z%f;]&cp

11.8水利论文Pf9M}W-d:wm

5.03水利论文(A9@9c5s:k'\+i pq2gGi

4.729水利论文Y6l\,{.r)]s

-0.301水利论文HxA1]aq

-6.0

.sz*e}:Ne Sw/_lW0

3.00水利论文-Q4v1^ZwXF

2.987水利论文f'kH)ag

-0.013水利论文1],C0G?5V+JI,l

-0.4

.Rq2d7K&?q0

3.55

)a)qt:_8W1l0

3.826

,P)y {)H;|:h)h0

0.276水利论文^~;M)h(I Z i8_

7.8

[J-r-vc\0

1.31水利论文)X B3O'{6k6L

1.319水利论文5WDYRT"yds z

0.009

\*V1ndL5L V0

0.7水利论文t-^nv&l5HZ

2.33水利论文)T$Vp![5{}-n5{+[,N;O-L

2.461水利论文T\9E`~!u bkei

0.131

(Se kMag!y0

5.6

I/mgC0a7Fs0

3.35

k3Xq&M C"p2{F;tm H0

2.780水利论文[i{i#| j

-0.570水利论文 tF/p1O4V[C!D

-17.0

L9i ]/];a] deDi0

3.55水利论文f/yr\$X c!yV

2.869

1eZ _E&W.c+g5C0

-0.681

o:H1Ng5q$C8W0

-19.2

+M#X-E hErr0

2.55

@I q Pv3goU0

4.305

"p1o],mw$|yO0zA0

1.755

j*]4GdU l0

68.8水利论文5B.N7UO%pUu1f j


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

表3 各因子的相对权值水利论文$I$},C(L9[m's0s

Relative weighted values of factors水利论文8P5S:L6a8n*t;U


权序水利论文[wDf j ^L9p2D/o

1水利论文G5D W@X,aW

2

1k`$`2` xvu??0

3水利论文7H _3X|8V;wF{i

4水利论文zr%}9^(r5Ya/U


因子

V2Fqs jp%]0

X4水利论文:lX1F E e{

X3

Q9k/U R#u0

X1水利论文7|1_SyTM r$R

X2水利论文F j mozHv`*r7j/u)W

相对权值水利论文nQi!x A m7[2KB.j

1.00000

M6?T%w OXvx0M0

0.85190水利论文-H.hI3Av

0.71617

8O2h)M0z3O4_a0

0.29765水利论文{rHJe n_j


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

J!l$u!A"T$Sq0

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

4 结论水利论文vW3t7U3G

990112t1.gif (1658 bytes)

By6EK/L[;A&S+r0

图1 某流域产沙量的PP回归拟合及预留检验效果水利论文&yiXaV3j#~|M

Effects of fitting and verification of sediment yield

6H no)v;I0

in a watershed by using PPR model

pj4kBX-]evi c0

     1.PP回归采用了“审视数据—模拟—预测”探索性数据分析新途径,建立的模型稳健性和抗干扰性好,因而预测模型有较高的精度。水利论文'G!r.HG}E$ye ? gO+@

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

     3.PP回归用于流域输沙量等预测建模,只需直接使用原始监测数据,不需要对数据作预处理,避免了人为干预,客观性好。水利论文{\7J.g-R aj1gj L

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

^yb*X2x0

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

5AZ2IE@w0

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

ZIi2B_NNb0

参 考 文 献水利论文N'Q;D d"\:H

1 陈浩。水土保持学报。1992,6(2):17-23.水利论文l~Y.XfS0dE

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

;pc#{%| d?m&P+_ N0

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

iW U],EPv0

4 Huber P.J..Ann. Statistics, 1985,13:435.水利论文,u"HmiM K

5 郑祖国。水文。1994,(4):6-10.水利论文:@!@jr,K }7a

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

4Wp Abm5c$`0

 

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

广告投放

广告投放