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

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

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

'A:^rXZ6L*f!Z"p0
李祚泳 邓新民侯宇光
(成都气象学院)(四川联合大学)

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

关键词流域 输沙量 投影寻踪回归 预测水利论文]$_QJ-h

1 引言水利论文yx M8Pc-g'be2bx#R+`

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

s6Y'b Z7[:I0GR0

2 PPR原理及算法简介水利论文j1R {9Ly7G@r

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

H n'meY0v2`W5O/ltD0

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

0b8R-l+P'w)]y0

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

6t tt2V,N0

990112e1.gif (1513 bytes)

[Ya+X&XI$c0

(1)水利论文 Y c1p Q-J [b u8S)n

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

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

n a/\1F8hRTD0

990112e5.gif (1470 bytes)水利论文B0Q[&W+U/O^ P%W

(2)

d&S!o:_ yx.JE Uq0

  它实际上是采用分层分组迭代交替优化方法对式(2)中的参数α,β,Mu和岭函数Gm寻优。实现步骤为水利论文5V6~f/r,s ^

  ①给定一个初始模型;水利论文?;kTD!J0b

  ②把数据投影到一个低维空间上,找出数据与现有模型相差最大的投影,这就表明在这个投影中含有现有模型中没有反映的结构;水利论文/eJ+DF3E6~"^.b J

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

"B7K[:h3|%r,qJ0

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

)vg8Hbj#SG(M?:|Lr0

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

990112e6.gif (1600 bytes)

:~ri?n a[0

(3)水利论文4t nv e1U7YQ/|ir

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

uDj"G4C;wtH"N0

     用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]。水利论文9] \u3`}J7P

3 某流域的年均含沙量的PPR回归预测水利论文/^.O1?!B?#y%I0g:P

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

[I!i`U6Tp0

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

     用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个因子对年均含沙量影响从大到小依次是年径流量>降雨量>采伐面积>采伐量。水利论文?Gdrn^2e }o%l

表1 某流域年均含沙量及有关因子实测值水利论文_3PdNRmR2G

Yearly average sediment concentration and measured data of related factors

X bt1yWym3D0

序号

5mx&A6Ol1o'?l0

年份水利论文U;zj;k.x Y!I6m/q

采伐面积水利论文 m"fP8t6m8so

X1(M2水利论文0FV!Hs#Lqu&r

采伐量

!BII,U0Bg"[V0

X2(M3水利论文H1j7G"] U

降雨量水利论文'`;Y ?;`$` Z}-`8gn-W

X3(mm)

'JiN&s;Dd|0

年平均流量

T lMJ}j }9`0

X4(m3/s)水利论文Q6`BX"}Pw

年平均含沙量

U'n `kx5g G6[?0

Y(kg/m3

`n#u4S'A1s:sWU0

1水利论文1w5V(Bu%f\(k0\0]

1961水利论文3`U(eu a)}&[6}5Su

15100

%ZS tQIJ0

29000水利论文7|klb7E4X5CR?

508.0

l5{*O]C0

180

M-P1^)p#{7f*[K0

3.40

3XV/l"{Rh&Ov0

2水利论文{?[,EMco-q

1962水利论文9Y#o3uz/^RJ

20500水利论文5ray)p'vB\1a

43350

d tX7t6b iw5f*z7y0

453.7水利论文$J8z})J r p+ZJk5F

141水利论文8R+o(JER;z Z*L r

3.00水利论文\ yp$wk^:t9E~

3水利论文8z&P'q G l#G.i'Q5u

1963

'~!Qp#j O+|t(_){0

80000水利论文/b%c4x/y |

92610水利论文;hs X5V qu

487.9水利论文:?j.A:},rZh

132水利论文8s7UBM)KY.e

2.70水利论文7T$}o!u_#P1q

4水利论文 ~oQTn2pV-L e1~

1964水利论文$KA6e,N~n~`4`

91500水利论文\@&abh

142775水利论文h-M `6S/`(q/NG*K#_:hA

572.3水利论文$b8fYu1A't&g}V1}~

182

R+B+r q#r${ G0

3.37

rm4uU Y9i%h0

5

/q$^ZV:Yb#_0

1965

N.?"{)x0OVp/Q$f0

82500

U|7`NHr7E0

2123160水利论文W%a.xm*?3^

455.7

({G%~;flW{m0

113

\'RQ.JT;t.c5ZR1DI~0

6.894

Ce TxBF F6L%p o"Q0

6

Bb5] y;T f7Q9^}"A3R}0

1966

y7C%^N ZD0

20000水利论文vv s7Ts0e_!Vw~

227800水利论文R8I4UU~!q

481.3

&GH3n'k,j$^;X]0

170水利论文-PRiiw%BE*|s Q.L

5.03水利论文lT$}+C3S4k[KE

7

P'fX8{Pa.l0

1967水利论文H/d-syfG_ I

17800水利论文 Cu;gT4r^u4s

140000水利论文?)z np f U&Nq

541.3水利论文n5O'y3k_5ZG ~ol

199水利论文W1M9D d V

3.55

U8~F"SW0

8

H0vwFCg0

1968水利论文9{9FDJQ ~*S] z

3900

"qWPl6~1P z+M0Lz0

15980水利论文%r'k5Dmxi`

538.6水利论文u8\8I PKL

186水利论文^2cQ;}#zs

2.72

+k thh&C0

9水利论文9T5ejr`+^4y0{

1969水利论文6Sn(cfTv4[:] U

17300

c7a2D{fkek0

223200

Nw-J*B+GKTg`0

460.6水利论文0J[(d;m5Vxw

99.8

}^-e.wjt~x |I0

4.05

\ X\7[ ^.Fy*v;Q$w0

10

W:_0wD!I\\0

1970

O ~"_1H p0y B5_5I0

25700

I%{K8hjhtL/~0

229400水利论文0?&E%R8^Gl~~

393.1

D C.glS `#W!n0

133水利论文qjJJCQ3z,W

3.22水利论文F6P ]"B"Z7t

11

O J9om4W mM0

1971水利论文s8w*pR {,L |+r._

49400

R-VaP G Hh1v M0

424500

%X Nv$hd0

373.9水利论文+Zn5]9rZ

106

P`#W:A8Cb*B@!X2|h0

2.65

ZXq$nYO0

12

MJp.vZI-E+zL b*t0

1972水利论文SXw.^#U"WE&f

40700水利论文$g5Z.u9Yj1v1l o

561700

\8JjH%]3JI vz0

428.8

(q*Vl |{:Tbv0

107水利论文XV4j4em+l}a~

1.91

%}!i Mvj0

13水利论文D+_WUi8dr6i

1973水利论文,Dv D+[s9? P

77000

y(ov'mlS0

563600水利论文'x%qm6SJKf7`5S

482.1水利论文 g_ |Q]C

140

E7z\&w;A5aiU Wy0

3.00水利论文*MT0Y/ATH!F EP

14

i RZKit0

1974

Wu%dz4~ {9oe7@,Vj9?0

72900

$I~ |,s"T0

557600水利论文#jo*on&Y3p'C+b6u

415.1

x"Y&\ Oi:uIo0

121水利论文}}.SG Q8}Fn

1.31水利论文 U!M6r-n@ S$L6i

15水利论文L ~^:E1tWk r

1975

q_#S+}'T0

63300

.ihWm4P/g6_e0

528300

\'E:rh Oy`0

536.7

,^{&] VM+cl6n`$n0

144

~u}~aW2ut0

2.33

A ?3HS9@[:[L0

16

qwplkg0

1976水利论文4~ n'r{(o O-C7Nf7y

51600水利论文5e&r#l0u5\ Z

488940水利论文xi+Q,J;D%@p-j.i3O)Lby

385.1

IKa9WPjOI)]6n0

154水利论文&z6L^YL Hz J

3.55

"Rg hc.}N r(^0

17水利论文xgrhsht p&W!N+k

1977水利论文pa%Mh2W-y

60000水利论文(A'X|+vcV&NKi[

480500水利论文,[h0jRhe-}

412.2水利论文 sfuJ4wlb

111

$KF0QJ i+[p^v0

3.35水利论文\? Qtf q8mc}Wq

18水利论文.rK0a,{)D*L-N

1978

M9{R.{F)m\*O0

70000水利论文/@1qtN;ke k1sh T

530500

R)u%x0s B0

567.1

:otY#Q!i L0d\!\0

139

'gt#R0uAd|eANe5p0q%j0

2.55

[*a.f3UT&gE0

表2 年平均含沙量PPR模型拟合和预留检验结果水利论文YV ewQT g

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

i1g(hv#t0

实测值水利论文s;{2Mh7PWx1C

计算值

:W+RrRZJ0

绝对误差水利论文e*y c9{B q

相对误差(%)

5}7|)vrv$D0

实测值

A,y7^;J*K x&tV/b%{0

计算值水利论文Y$~5p s6^ }

绝对误差水利论文;yVXg!{&Lio(n

相对误差(%)

2mC4_7J9V;l:Z-ER.]S0

3.40水利论文w:oe1jl8zXzp7M

3.131水利论文)l)~KK;l8j4M;|#z

-0.087水利论文Vy%yVj

-2.5

-u t J#z^Hw X'? b0

2.72水利论文9~*I#n:ZnLm)x([$M'r

2.934

2y|9E mQH5_])EQ0

0.214水利论文t}O7MQK Fr+?

7.9水利论文2@0L Z{T

3.00水利论文`I9t1z+N)d]jq t

2.738水利论文rf{,_9p@B?-Eu

-0.262

1Ko@Z\G0

-8.7水利论文/oBc4\8hH.B

4.05

L&]6f7d W_'X0

3.996水利论文+[ oSL!jdQ

-0.054水利论文9f$F[1d|!z X,|.}

-1.3水利论文&S [:Uj/z.Z

2.70

3^ N4yv { X0

2.659水利论文ze}/M3kp3E"h

-0.041水利论文,]8F} L7~B b'dX

-1.5

%gag{;WB|.I+}0

3.22

Hg(dk+J2oE)_ Z-v0

3.214

j"om)q8SW u O7M0

-0.006水利论文4H(gc4}K r,A

-0.2

*g0g [+edS iJ0

3.37

!zY2s5] NZ{K/R0

3.400水利论文-y!{&vU8](WC*G|

0.030

F9p1L-qR$s0

0.9水利论文 a'F#}?FKxW

2.65水利论文8EQsrL_4^

2.628

7[5S!?7^~g2|dE0

0.022

*Y7i@ma~K0

-0.8水利论文_[*]|^!Z;Ij? B

6.894

b&t'X$_~$l!n4_6G n-fQ0

6.926

/\f N @0HN0

0.032水利论文 k!~ Y;kAE!BT(C0I n(W

0.5

IDN&[$T(SH?sr0

1.91

2TD1x'qKsme\0t,z0

2.135

.@&b A9{ TD0

0.225

p%q5\(A'n"q UM!G5i;F0

11.8

@2xF6pZ;m!l%I&}0

5.03

'b{?7s:`#j\ I!D0

4.729

4G+cX:a"Rr0

-0.301

)EIz Da+Q0

-6.0水利论文O*MhtbX'\{ `"j

3.00

9~ l'_"E2m)^rl0

2.987水利论文"P0u,W|P@v;f&W

-0.013

L6?GBR~L2M0

-0.4

|*`{,~Jz)o FC0

3.55水利论文%b:]f+\4\y}

3.826

.GTP D[7v:ri0

0.276水利论文}^1y3Y;{w!Ng

7.8水利论文|b@H"M[^

1.31

]$[ jX9UC U N0P0

1.319水利论文,y%D8va1i:p-o O

0.009

@,ueP4{4C p"wVX0

0.7

l.z$d:L9i0

2.33水利论文_)o1\.hr:s [

2.461

"C1a;q:Q:roz0

0.131水利论文 {)v?3C"ie5I2c;b

5.6水利论文1T _i2tVz{n

3.35水利论文@Z P^u@/|G

2.780

.lTZ9kqS:r0

-0.570水利论文(P wPR1_|)Q

-17.0

d X b(a)T@8R~ O$ER0

3.55水利论文F} n4M0i)d

2.869

*}uk!|6qn!x0

-0.681水利论文C4?Rt)hA9?9|*c

-19.2

U#X(F.M*m:\#lsY A0

2.55水利论文9OrG8uc|#pu

4.305水利论文9D'k^+{$j3B1I3NuE2^

1.755水利论文rrjR`

68.8水利论文O"X j3Q`:B.FlZE'd


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

表3 各因子的相对权值水利论文myM#h `Y!x

Relative weighted values of factors

:fy.I@r-tS[0\']0

权序

6A}-MlR j(?E4cH_0^l0

1

KO8ze2_Sb0

2水利论文MFGhWq

3

&JUy4P[ x9hZ0u0

4

[4v x+Rg0

因子水利论文3|[xy},pZ

X4水利论文3\K$`_L2K

X3水利论文9t3{Y?g5h

X1水利论文0|)I8K%B4k[Z1S8qC

X2水利论文"zK%l_? C

相对权值水利论文i4pk.`FW

1.00000水利论文;pR._X9b:H#x [c {,NB

0.85190水利论文9K2x*CY0u%p

0.71617

?&J`d,tyms)G0

0.29765水利论文!e5R!_,T@k I@F&v9NAM


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

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

4 结论

iY'M9C:X ma0

990112t1.gif (1658 bytes)

`|k4K0X Uk0

图1 某流域产沙量的PP回归拟合及预留检验效果水利论文K@ l&Yb~9hg

Effects of fitting and verification of sediment yield水利论文7E0s2P$\}*lG)xW

in a watershed by using PPR model

N4lQ+oD5g"F0

     1.PP回归采用了“审视数据—模拟—预测”探索性数据分析新途径,建立的模型稳健性和抗干扰性好,因而预测模型有较高的精度。水利论文9r wte#Y:?a _

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

#`7m7U H6q+^-s0

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

4T4B+iO3o,z;W0

     4.PP回归也可用于水土流失中其它指标的预测建模。水利论文%J/i%yJ@

   5.本文用PP回归技术建立了流域年均含沙量的预测模型,求得年均含沙量S后,只要将年径流量QT×年均含沙量S就可得到流域的年输沙量QsT。水利论文#n*Y Xx+O)Z#qf N

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

参 考 文 献

D6[g'cI9z0

1 陈浩。水土保持学报。1992,6(2):17-23.水利论文BR)V6q[)lY7t

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

V5mARSjaN*V0

3 Friedman J.H.&. Stuetzle W.,J. Amer. Statis. Assoc. 1981, 76:817.水利论文!qU6x"_#sp

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

r.O2L,p2Flc}#k0

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

-E2yzH)?k"F/Gox5A0

6 吴建明编著。环境统计学。北京:中国环境科学出版社,1991:504.水利论文mN jt9Qn$KJCB&|

 

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

广告投放

广告投放