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

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

流域年均含沙量的PP回归预测水利论文Y*{~DM6{]/`!X

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

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

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

Ur#r*o+@9RD^/A0

1 引言

~,?-f-?X P4K0

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

2 PPR原理及算法简介

W4yH~H }0

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

7c2h!n{]3G-g!a0

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

G Z7Uzw!c0

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

990112e1.gif (1513 bytes)

MV[)Y8u4~GM9[0

(1)水利论文R ](og2cH[2j"i

式中Gm(Z)表示第m个岭函数,Z=(990112e2.gif (928 bytes))为岭函数的自变量,它是向量990112e3.gif (888 bytes)990112e4.gif (879 bytes)方向上的投影,990112e4.gif (879 bytes)也为某方向的P维向量,M为岭函数的个数。水利论文HA0]l"H;l0} Q

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

U}.m!v@c)W/h0

990112e5.gif (1470 bytes)

(r(VBW&e`!fr~0

(2)

L/y)x i1F7f0

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

;PjmPi p0

  ①给定一个初始模型;

zA W+U`e"VS0

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

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

Nh!X w a J"o zY0

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

d XN2?G+\0

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

'Y/e|h,T Wpv0

990112e6.gif (1600 bytes)水利论文4Ec F/{J

(3)

4H%|C5d+k0

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

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

]D;ik3m0

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

Yd xS&P8`4l0

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

y*\J/`-V&~u S0

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

2tp%i_7DMX'P \{0

     用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+E4A%lN4u o

表1 某流域年均含沙量及有关因子实测值水利论文`-{r,k!A3T l

Yearly average sediment concentration and measured data of related factors

er+Nj8lv+j)^0

序号水利论文teE PI+s6of

年份水利论文dS1fV0U/|'?

采伐面积

_6\'`"ww\@i0{#L0

X1(M2

+k'n8uG1X9o0

采伐量

VfmFBXb%C9i0

X2(M3水利论文"ne$F }^+LF~

降雨量水利论文_yfr@5b

X3(mm)

\DBg bK0

年平均流量水利论文vG/CY ?'IC@8~

X4(m3/s)

B5l;\*GX7VQ0

年平均含沙量水利论文*Z*?&HZ8lY&{'o

Y(kg/m3水利论文1rv%y:wsc \V5Eh%]/~(M


1

D[2dI p Q3L,v0

1961水利论文 hQ@"dDU?)M6N

15100水利论文/f!s6S,}a$zG B aB

29000

)dh&E-YaW'K _0

508.0水利论文pA)TEvI

180

k w}P%`B0

3.40

(_c`vNy:U-\S0

2

FgYI%N;QDn*j0

1962

SN J4B } K0

20500水利论文#V%}z5|1LO8e^

43350

YXie#K)l&~ Ka-w0

453.7水利论文)zf'UyQC&d!X!^

141

3LfTk,F-gcl0

3.00水利论文J)SL7hrZb2l

3

0w&~U$X8s$u0

1963水利论文!WG,tDh5l KH[

80000

7WP iCS{px K0

92610

mj3b!C,\$iA"N9g*C9px2}0

487.9水利论文K lL m6Vf^4J7Bp3G

132

:gKXSC){?;H0

2.70

^2B%^P @~U0

4水利论文4E*o+^|)oX w4t~9L`

1964水利论文8}F2J$B/L'ZJI)wVtM?w

91500水利论文,D5Ey0y9Vy KR

142775

~H*f'YLEQxE0

572.3

hLOKB N:?1@0

182水利论文6g0Wd#m0P0p

3.37

'O8dr)T'`rc0

5水利论文4N$DFpU|e2h7S

1965

lY^2z m0

82500水利论文[CD$l!Ca`1AAS

2123160水利论文9[+uu^1h%Eu

455.7

j;UO-Y3S D+VLf4L0

113水利论文s2H S1W%@c)n5D6o5Z

6.894水利论文y H6P@,[ Z:U,E

6

%S1O6z8BLN H0

1966水利论文6{'Dh4h~

20000水利论文b;Y Qp,Ab4C

227800

rz uo n5m'Oe0

481.3

Z3B gGU0

170

#Nx v-RR1N]g0

5.03水利论文H4p i6z0Ta%U*|?

7

._1g1Y?gqOH0

1967水利论文)@}hB?JB~*u,Q

17800

(x,o.d Z[ _DJ0

140000水利论文.q%fJm`Q

541.3

W|7NP%M gh [0

199水利论文/Tzo/W ^:ye"MU}

3.55

t:K;C)R_e0c h/F1J.G0

8水利论文Mgg].c/B _&y*_0d

1968

{/g;G0N/?I0

3900

.}wa~O,C4_2x0

15980

P lRV:I U7q#[0

538.6

i$KeK&Y n8hY#M _E'j0

186水利论文-s4k&I%G.T`Sf9Y^

2.72

Z'j7S?Y \&]0

9

z v4G!@'O$c/{^7f;\0

1969

/Q8p Re(q`W)f0

17300

'Yj+v)lk1{e0

223200

v4P4h b)v0

460.6

z W[K s'|"_B9yf0

99.8水利论文N5D Z2|2SDqV

4.05水利论文.`/H H Y)bWT#FP

10水利论文z3VYo0XA a

1970水利论文,r{y5D)u

25700

^Q6@9l$p3rl0

229400

P H*U q/n?J*t0

393.1水利论文3T;o,A'H-D`/ba

133

f.j%D~2m*c0

3.22

*H;\X7f4?0

11

5bsC*^ v'aO:xD$m0

1971

Fk2QP-s'G u}0

49400

5B1C+t%zy z9`/s _0

424500

YSz,lDyLR0

373.9

?j(u&K0q2fp{_0

106

t0PL9u$~.tFz4Y0

2.65

Wwz*z tQ0

12

'l5Pv }S8wc/d0

1972

T(Xr r9bW0

40700

Jl:`5JV']qoVq"r0

561700水利论文DX2|#[y\(R|iH

428.8

NI8jw.C/[2L)I0

107

|;Vkx q:M2`]7Z1A]u0

1.91水利论文 X.?$O8y0FH:V

13

_~%H)g-U;X-p0

1973水利论文V wHj&m#D_O} a

77000水利论文3y%A8xu q@8~

563600

g8u P-se!e2k,OVZ0

482.1水利论文{%e P)l0W t

140

5WX"@ |B @[z0

3.00

:_ r _F)z~O-A0

14

}$kVN*^{G/T0

1974

H{k"u{5S0

72900水利论文;XX(y,x`)G6x

557600

*aB g,`&a5MI,Zpw0

415.1

&`"{_!G!G_0

121水利论文6^\J%v0]0Fs3m

1.31水利论文7@+_&| m(T;N G0lE3E L

15

FH:xou#p\1zz t0

1975水利论文G9G!R MP8m-eKo

63300

)B4C1s"b4VJQ)Y+}0

528300

*L7W$[M6J1S6J0

536.7

T*jgy&t(E0l0

144水利论文{6h"A.@u,m

2.33水利论文zCQ#[{/G*se9Gr

16

$WteZ[7~NRg0

1976

.\,Z!s|g,Tb0

51600

Wy%F"Qv#a1oZ-lmH0

488940

I9^Y Cu2l]l0

385.1

![8b OX hF9S0

154水利论文3d Dk%ZF@B

3.55水利论文J I*?4vXg

17

D'n to zV\D.Z0

1977

%Nv^sH!h5lB0

60000水利论文y+{#j`1Q

480500水利论文2d"[2`~2Fu+c`.`Q i-td

412.2

/j/IWGF6`0

111

LlJ"B)m*P!pt0

3.35水利论文(gTT [6c

18水利论文m%X.K Fn!n6jH+ic

1978

!gtx$J9uWLK0

70000水利论文,N;i.g1}5s\&y"@ w

530500

oBpp+{w h&k+Pw2^"\0

567.1水利论文D ?)J] X5bXj:i

139

{8]0i;s ]e0

2.55

e B`,fG_| A*YI0

表2 年平均含沙量PPR模型拟合和预留检验结果水利论文6E'}@+X$v-fR5@(X3i

Results of fitting verification of yearly average concentration by using PPR model水利论文 @,a:w\+m!t?2J


实测值

*X JuQ2H}Qq Vz0

计算值水利论文U? vZw,zzw&a%P

绝对误差水利论文*fmcqLu T

相对误差(%)

[.EF0e7PZ0

实测值

kQyP\*rw0

计算值

3i9_B3^3j H r Yh!T0

绝对误差水利论文_f.n7p x7Il-d'T

相对误差(%)水利论文9j:],q}]zJ\ ^N#t.B


3.40水利论文2F6HcUO9J(y

3.131

H/x#wlg0_%Cl4yQ0

-0.087水利论文F hr7Ax

-2.5

2kc#lw{y,u B;g0

2.72

(^'p ].j7mqQ4o0

2.934

\'b/N-H bKl%n*w"z0

0.214

#f5x+q$L EGr0

7.9水利论文,},utn+ipy0sf o

3.00

C4sx'Z QwD!ih.k7]0

2.738水利论文(zOUuEz

-0.262水利论文7_'[#l!Jc gUO

-8.7

Z8[`$|RMe0

4.05水利论文)UW [ nIC f

3.996水利论文$O"@YQFZY

-0.054

w$gVBv+Y0

-1.3

l!aV$d[f^"_"f"A0

2.70

iDfGD2~$xX|T0

2.659水利论文T;Dr ly.{!Y

-0.041水利论文|/t2L){%S#v

-1.5

%`0A:J!\ E6V0

3.22

5@'s7Q&iO2f7`0

3.214

{7b5b*^q%U+{5c0

-0.006

SW n+de0

-0.2水利论文qA rbwB Qd)~

3.37

)}&_-bKS$P)i8~0

3.400

?7L~ Gmvw KGI0

0.030

GCq)ANHZ0

0.9

3Jwds6X0

2.65

1kAc-llXK0

2.628水利论文f*at-PmdxaR~$L

0.022

/{wIX*VUR%QA0

-0.8水利论文6\9\_8sN&Z

6.894

p C)T#])r'W0

6.926水利论文%igG)u#oMa T

0.032

5a/X0p(mP.aW&Tc0

0.5

3[&BY M$} x'BVy0

1.91

"h0E;LM nw0

2.135水利论文 ]r%C]a

0.225水利论文pS6jCz Q

11.8

,mgc#\+|%| M5r N0

5.03水利论文&d+E0h:Yz9l9U7U

4.729

)l+F9Dlmx,v3P0

-0.301水利论文 L3hO%}SuP9t

-6.0水利论文9yH[m X J&P

3.00

3X.Dp&XQQh"G0

2.987

)y9h)o`0I f)z0

-0.013

%^x.Hd)p#_&uXU)]0

-0.4水利论文7y"Lq T3Ty0M G#Y

3.55水利论文sM7e&^9\E i$Kp3Y

3.826水利论文+p/e~lx6J{(v

0.276水利论文O/RO ^[&g

7.8

u(?O|:tV0

1.31水利论文 zc G%m(pF#{4es0j

1.319水利论文_#~{%[ jH

0.009

+z5}+t-yPg5AQU0

0.7水利论文7c+G;w,o.E/XB

2.33

P9F.LZU0

2.461

$rnD2S&U&]$N!L0

0.131水利论文8DnJ4X/I?j Nie

5.6水利论文I)fgD r5x]_;K

3.35水利论文K vO`F6A(b

2.780

M-Ja d3C_ Q0

-0.570水利论文)yEX3N7_E YVN

-17.0

av8nGVU0

3.55水利论文.n)|pa3bx1_q

2.869

,kqZ8Ix6?:M.s'E,Em0

-0.681水利论文)Qt'[,kB,NMND

-19.2

_@"K%ALH6e0

2.55

IvY7P I(RcM0

4.305水利论文&vKX6O3W^ i

1.755水利论文(bnHzdt

68.8

Jj_Ey-LW m0R:E0

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

表3 各因子的相对权值水利论文7f7K _%S/^;d\{:Ri~:i

Relative weighted values of factors水利论文3}O MF/r


权序

`r*G1Za j L0

1

Xp)w/C$o NM)~D0

2水利论文bo t;o'lue;\

3水利论文tAgAu/uk

4

[3{y-q)b(s\:[0

因子水利论文(r/cj j5C9B

X4水利论文@ aI7C:BK

X3水利论文w'p1V9x6`!dQ*V

X1水利论文;O s|I ?$C;cu&H

X2

2v5pw K u8rdO0

相对权值

S,]^.X P D0

1.00000水利论文6pv,J*@0S.c

0.85190水利论文UB&_gEqBE&UdTk

0.71617

1N5u4sF,P5L"~F5J0

0.29765水利论文$U~ Z,V4}v ^!ol


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

GyO,S*l5o,r/z0

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

uT!og!e8t x3l0

4 结论水利论文$m1U+G/t,a

990112t1.gif (1658 bytes)

vrIO$J ?0

图1 某流域产沙量的PP回归拟合及预留检验效果水利论文mQvn&m

Effects of fitting and verification of sediment yield水利论文;IM1_;zLFF&P

in a watershed by using PPR model

m o;i?!|G0

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

e}9ycx~N:I0

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

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

     4.PP回归也可用于水土流失中其它指标的预测建模。水利论文w0K1uq8w)Q%E

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

R M#qyUe5i1@0

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

4JkW"PAu*P:p'OC V0

参 考 文 献水利论文,J7BN]_$B ?6wiW

1 陈浩。水土保持学报。1992,6(2):17-23.水利论文.gq$p7B&G;~'S!HrO

2 张醒。泥沙研究。1991,(1):15-19.水利论文{+f0\ e*H

3 Friedman J.H.&. Stuetzle W.,J. Amer. Statis. Assoc. 1981, 76:817.水利论文'p0F,k:F#E

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

.bLYNvA0

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

`N[]/?7m-kb$e0

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

Y%OU,y C"Y ].@ j9l0

 

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

广告投放

广告投放