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

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

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

2v@%ki~Dzk0
李祚泳 邓新民侯宇光
(成都气象学院)(四川联合大学)

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

!kvtYq6j%_i$b;K|0

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

Kez+t?(T0

1 引言

0RI|4`*jl a#AL0

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

2 PPR原理及算法简介

BRG/o#F4|G2K0

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

d(_4Dzq9^h6A!A0

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

oaZ:oE0

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

990112e1.gif (1513 bytes)水利论文5v3eUsY

(1)水利论文*r\D)|al~r&Y

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

Friedman和Stuetzle提出了实现PPR的SMART多重平滑回归技术,SMART模型具有如下形式水利论文9DT+PC-y.^_[

990112e5.gif (1470 bytes)水利论文4r a |!Mzz

(2)水利论文~!GK3Nr'L2GI

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

&uNFN o0

  ①给定一个初始模型;

$Iy LzM5t i p.X5l0

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

$k9|y9e ^ I K.O0

     ③把上述投影中所包含的结构并到现有模型上,得到改进了的新模型;水利论文;L9tr8uF1s,]H

     ④再从这个新模型出发,重复以上步骤,直到数据与模型在任何投影空间都没有明显的差别为止。水利论文.KWZ;E:|v c-lX(}

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

990112e6.gif (1600 bytes)水利论文x7Z)CQ]Jw)@|

(3)

ue7if#[0

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

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

%E v.v1n9I&d)OF0

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

UQ @}_ea;\0

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

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

t`K(d4p n5Ki0

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

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

!Ir5o;i_Xy B0

Yearly average sediment concentration and measured data of related factors

(A|SR&?}0

序号水利论文5kq/c7^k.J

年份

(ZNBK(X q5{0

采伐面积

l$P5r6{G,Q@0

X1(M2

WM#{v${r P0

采伐量

(KRZC+k0

X2(M3

I7[N I mj0CQ/u~0

降雨量水利论文 ]"]R9y4[B

X3(mm)水利论文Nd`]6frt

年平均流量

0s6qJy$B;L\J}0

X4(m3/s)水利论文]#tE }'b9l#Z'C$?4D

年平均含沙量水利论文Fg7F^vFk

Y(kg/m3水利论文1ykY1\ LV[#EN d6|U


1

Ma-cc0Oz*pf7y0

1961水利论文jTVh1eX)[jU

15100水利论文ui"@8eG8[,Y'u1c

29000

4D*L X4JU:b5q,O0

508.0水利论文#u fa+Q\P%W o/r~

180水利论文)K2t3psq [ d6N#G

3.40

/qZ2KoloV0

2水利论文bBA;[E,B"i&m

1962

0{&A8]rnu0

20500水利论文5E%B'm}`!n

43350

8l4z X J1i4Mk0

453.7水利论文)A:GGR`2gq@

141

Q;T$P iJ2x^0

3.00水利论文 Q![E(jY&k MOF9r

3水利论文F$Ia @Mv

1963

G;B7[,a%j J&r+{X/O0

80000水利论文'^[*{7qP0X@HBl

92610

2[g:cS K;xg`'y0

487.9水利论文-U'Y^~Y#i/@

132水利论文_@8Ky3~Y3c7z`

2.70水利论文%GYA!e6d)Bc1vt

4

-Vl#x)r,B8r0

1964

h0`!w)bz*Q@8S0

91500水利论文a)U"A"d H

142775水利论文4o]1?1u2ll"WoI

572.3水利论文E:g ~ H ?.NB&Lc9M b

182水利论文7j/t5R'[D

3.37

9_RE/t TO0

5

9NM Scz U[J7r0

1965水利论文F$@f!P;Pn7L

82500水利论文3vE4@'z ba+m&T

2123160

K)a-yBg0N0

455.7水利论文 J5YRS|x5H]O

113水利论文'e`P5kF/G

6.894水利论文,d4P4oTz3W

6水利论文+l z#ih%O6OS

1966水利论文W'p9AYZ~S1lh

20000水利论文2S6[T?z4F:Yt

227800

2s&Py dpU!jY0

481.3

+e"T"Wfp q k}0

170水利论文 E;M:Y(D#[;Qm;Y/f

5.03水利论文]&{{?zo

7水利论文[H l}k7i;L2c

1967水利论文)OV0L8N~9yQ`

17800水利论文t1SWJz|

140000

F4NfJ w ur(I [0

541.3水利论文H_$g'BxH}/N

199水利论文?#X Of0}6y5~6Z

3.55水利论文8cf)w0WpmY(W

8

FQ0h T ? [n2g0

1968

5H+Ln!Ck8u1q Z0

3900水利论文wjG5Ah!io6z},n'e8Q

15980

9I IAo1_f#Y0

538.6

n:C8@ FRy)L7a0

186

4VsKFW*uLi0

2.72水利论文 bsnAN

9水利论文$jW{J*{:pU&{yY o

1969水利论文lY;d9GxX'y

17300水利论文t rsd`6Y4@%\

223200

a`r4Aq2n0

460.6

#z A ?:n'k0

99.8水利论文#\,Y$]8q'Oxk,Jh

4.05水利论文r b5TkF+np

10

%[-`EmW H_:Hv0

1970水利论文 ffdP@oO1] |(w

25700水利论文zG&_8n*[_

229400水利论文Jv f Bk

393.1

N FA l ?\uhnB0

133水利论文Y4}`$vf0R9xP `o

3.22水利论文:Z3aK}a;~$H_4D

11

.^#_ \W9Nts'x0

1971

#{d%Il |:|b~d0

49400水利论文 z#L-m"@J ORZ

424500水利论文y&vQ[Itl

373.9水利论文+H`^'xX N;J'r,ck

106

v1Ca$l$t:cs+r0

2.65

cXFQ G!X0

12水利论文8Fn0sGK7t#k3Z {

1972水利论文L*h N_7U(f%~U

40700水利论文a5L:_)OE;pcJnf

561700

&B]z%L@1[8k4M0

428.8

$JTF3} q1aBLg-DN0

107

[ E+y~^2?0

1.91

9jIn h%GM.q0

13水利论文t3R"}r S-Z:x,f1k3n

1973水利论文b\&}.MiKV

77000水利论文Z]O|'v$d$R

563600

+~%ZjL-IVw0

482.1

T-hv Njl8I0

140水利论文X D;K#n|9G&l}c

3.00水利论文HYK"{ ?F

14水利论文(L$b1y _U(^^

1974

%uv/D ~)U1D2N&J_0

72900

]&aQ3BlI3?&ZADh0

557600

&I`(S(McE |-Q0

415.1

G5DX%u%x0

121

!jT'q[)a9o0

1.31水利论文-Lm\ x*WC]A D:@

15

aX)k@3]7GG?0

1975

5g y;oT@Y6?0

63300水利论文W ^2_&j8Oy(\

528300水利论文0WGv Y1mRY#[

536.7

,m1d~ b9a'hM0

144

L`a-t%Z:n-f)q*R0

2.33水利论文4fMx9ZR)F0q

16水利论文/\i?nM5Z Y9Hu

1976水利论文U @fi-Kf0m

51600水利论文 t%b:d.|Q xM

488940

wjLb7[H&P&z0

385.1

!N ?$i+i}*u4f0

154水利论文I } M'i3W7B6`

3.55水利论文)Dp1EeNj,o:WL

17

v'`;{)k P0

1977

UaE aq;\&^ @0

60000水利论文fQ5J`F/l

480500

,^:c!l O(e+f0

412.2水利论文WW,J_a;v8A }1}

111

JgM B+ET r0

3.35水利论文P9X6if o(j#[

18

'a)R:G,SC"QZuVN(] p0

1978

eAD,UIt0

70000

q'E.a|7N s,s},s1e0

530500

8N \7b(bFux0

567.1

,Xn/Pbn.So+J K0

139

Cd7i*}8M Y)d0

2.55

"z!N8JT C[%m0

表2 年平均含沙量PPR模型拟合和预留检验结果水利论文%OuwWS,n,P

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

'`1?i8y$qN\G1h0

实测值

W"o G2A,s/d0

计算值

;I#Mz-`'q0

绝对误差

V e#Y*j;M@h0

相对误差(%)

N ?+T1X P z;t-u KZ0

实测值水利论文5sn/k.LJg(I

计算值

7p vz4Wb'I\N0

绝对误差

rBS|\0

相对误差(%)水利论文X%wk#Z%kGZ6D


3.40

6|'g1J Et)Br0

3.131

[)s J+sPP6o{?0

-0.087

P6},[$cW R0

-2.5

PU*_Ni?$s0

2.72水利论文mp|;PQz0U

2.934

H9g+?2N@f0

0.214水利论文YX ?b9zUV A2Z(@

7.9水利论文#G;e g$r5bf

3.00水利论文*NP {(w~aP[

2.738水利论文.e$v;B1]\5r

-0.262水利论文{&j8CS\D

-8.7

N[ O0Ce/Mq0

4.05

4g1D [(ct)T2c h0

3.996

Jn#`'V/b/]+p0

-0.054

AU&CPe,F;w0

-1.3

RQE{K&Fqk0

2.70

y-t&BL7E5}^Y1b+c0

2.659水利论文3U p0wm ^MO s

-0.041

$gR'apc-I z0

-1.5水利论文c#K1J&y;D!| mT

3.22

E&@ z7@3r8xd0

3.214水利论文9p;S~G@b5{4b%v

-0.006

r%F*g9l7Mv0

-0.2

Gt~a Y&S6q-i'FJ0

3.37水利论文R;?Z/dZ4x?5q[b

3.400

knk2Ca\'n.~0

0.030

'Yp1h\u&~@0

0.9水利论文sz h4w$F

2.65水利论文/XW6DmXu2Q)x{

2.628水利论文R5my!{G!j

0.022

,YLCsS0

-0.8

,R!_,b*v;z*X"m0

6.894水利论文:Hizw@I3wA'h

6.926水利论文 MYQ)Y%d-r1_:~

0.032

5h7k6t|?{"T0

0.5

&Q7^0]G@ I0

1.91

4q%}m,e6g/S-^${V0

2.135

!A~g;lDv Nd0

0.225水利论文l[1D4~"}9PI

11.8

8A S(z"L f0cv%T0

5.03水利论文b/Mb.t7n|%_-M

4.729

"A4Kb `h$nZJ0

-0.301水利论文+Hn0|e!fb9f

-6.0

W-VlL(v7]6HXy2h0

3.00水利论文*[2k`^9Z3V0RoD

2.987

,QW5Uph!Z.RU3n0

-0.013

&Il7@.L`F's0

-0.4水利论文9D8U0l-u5F^

3.55水利论文M5w&PdoT_b

3.826

k2ZCq ?kanS I0

0.276

sW]TM2v }s~(m;N0

7.8水利论文?quXg7QAxnN

1.31水利论文a5_A9O5}e

1.319水利论文 r$Z0@-r*EV^%BO}w

0.009水利论文&`YZGx

0.7水利论文)f6@Xz0D{*r[ O

2.33

ilB.e;K&^X0

2.461

P7ZZp)@p0

0.131

A c$d-hh1K0

5.6水利论文^P4|#zx/Y

3.35水利论文1Wx[.{!I3cO'H

2.780水利论文5b/[v bl {C)p

-0.570水利论文S.s-[{Jr0tb h2x0G$]

-17.0水利论文D.X8G&phUl5Tq3W

3.55

6Qp1s+^QWv|wz0

2.869水利论文K)Uyq0zz;Z

-0.681

+h)J] S F t`1?0

-19.2水利论文*pU-x3Iu)F n

2.55水利论文 CN4lT#D8U \j

4.305

F eD DX|0

1.755

2E)F%ez5O5w0

68.8水利论文6Jkj8g@A9M


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

表3 各因子的相对权值

HN7hH"[J^Rd0

Relative weighted values of factors水利论文g0g2s0[o


权序水利论文,O(\2H ~F"T,x0|E

1

X{%|!Y@;R0

2

4m-V2~I6?iv0

3水利论文/pDW|H GO~

4

UUQ/oH0

因子

RJ)W:Ml+pQa[0

X4

V3{&h'{6Ok0

X3水利论文9}K p3R U9TE

X1

&dB)L2X2S7{0

X2水利论文 e(g%LSuW&fkv)wp

相对权值

W"e#a`!I I5U:_0

1.00000

]O|5I'?.\9Z0

0.85190水利论文6y1_o`3QqL

0.71617水利论文 LOJ)jOg'l0n

0.29765水利论文#^*@W(U_


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

M[^#XH0Y yT0

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

4 结论水利论文l"~%sOrb E

990112t1.gif (1658 bytes)水利论文9|/W{L1vs-pkU

图1 某流域产沙量的PP回归拟合及预留检验效果水利论文4Q8ErO(q|

Effects of fitting and verification of sediment yield

$B#l~'lSJg(Yr{0

in a watershed by using PPR model水利论文$J_#He ?/O'TI

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

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

     3.PP回归用于流域输沙量等预测建模,只需直接使用原始监测数据,不需要对数据作预处理,避免了人为干预,客观性好。水利论文5h*a+s9MB9X7}3f

     4.PP回归也可用于水土流失中其它指标的预测建模。水利论文'@9n:mI~ B

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

r3j(v){-XW ^;SI.tc0

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

参 考 文 献水利论文"I,N1xsSHBO"P

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

9U ?!gMb^0

2 张醒。泥沙研究。1991,(1):15-19.水利论文~:c R8yB^JU C-g

3 Friedman J.H.&. Stuetzle W.,J. Amer. Statis. Assoc. 1981, 76:817.水利论文"K5VE E6}Omb

4 Huber P.J..Ann. Statistics, 1985,13:435.水利论文5v4Hh)E6J1H6^

5 郑祖国。水文。1994,(4):6-10.水利论文*Ea4p E3?1fe

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

2d:~0m4S bct0

 水利论文kg+iCd9O^O

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

广告投放

广告投放