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

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

流域年均含沙量的PP回归预测水利论文|1fO `0yup R B

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

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

关键词流域 输沙量 投影寻踪回归 预测水利论文9` oq-PXu~

1 引言水利论文V:F%c _4q

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

K:g;O7{E)Q|VU0

2 PPR原理及算法简介

VH ].]2SH;U y Q?0

     投影寻踪是国际统计界70年代兴起的高新技术,是应用数学,统计学和计算机技术的交叉学科,属前沿领域。水利论文$l+P8}_3|7z-~${wq

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

)C(] }*|L4Ql5`0

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

990112e1.gif (1513 bytes)水利论文i[}e%rvM4X

(1)

-P/ln3bL0

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

/Tm/GYU^ O0

Friedman和Stuetzle提出了实现PPR的SMART多重平滑回归技术,SMART模型具有如下形式水利论文"q*_@V$U

990112e5.gif (1470 bytes)水利论文bD!O7H E S

(2)水利论文R)B;_6Emp9Msi&q

  它实际上是采用分层分组迭代交替优化方法对式(2)中的参数α,β,Mu和岭函数Gm寻优。实现步骤为水利论文#t}'s$y9KY x ?

  ①给定一个初始模型;水利论文5Ii+^4P|+v

  ②把数据投影到一个低维空间上,找出数据与现有模型相差最大的投影,这就表明在这个投影中含有现有模型中没有反映的结构;水利论文DN}4y3O+wtA JL0U7l.D

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

-ZX%TQ N0

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

+W&C2s/cK7@0

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

990112e6.gif (1600 bytes)水利论文g })fa D1_ v

(3)水利论文0V'as"A4F(H

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

2Cwu_*?2k{3J0

     用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]。水利论文#rl nPFu'xuc|

3 某流域的年均含沙量的PPR回归预测水利论文wI{iYZ

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

Yt*O_!?:[2I!y0

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

     用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个因子对年均含沙量影响从大到小依次是年径流量>降雨量>采伐面积>采伐量。水利论文,W K-H'}7a8_ c'HF

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

kky+V3F0p)` F0

Yearly average sediment concentration and measured data of related factors水利论文7TY/ax5ck


序号水利论文u7Kg!Qhs;f,I

年份

9r-IdD*B2BkF0

采伐面积水利论文Dwx8C n%oV;}#K+{.w

X1(M2

Tk IwN0

采伐量水利论文]5@? J V0tO

X2(M3

3A#C cP QD.~0

降雨量

DKHI&k0

X3(mm)水利论文/n)Vw {;Iffst2M V

年平均流量水利论文,F4r-? R~%kF

X4(m3/s)水利论文*b*n5}h f8N

年平均含沙量

?8a7L.s.qq5Q0

Y(kg/m3

Umvcq0

1水利论文uk&gl`|

1961

w B t J-m6~0

15100水利论文"Me-P0f;v:^

29000水利论文N{4jI0g!]'Z

508.0

5u,_r~b0

180水利论文j'{?jjA

3.40

/]_X7Y5ITx,D0

2

{Lfz jgD}0

1962水利论文!U tq3o?Z S1b

20500水利论文U'{+k`Y

43350

%ZY'~F+oY }(N0

453.7

;\+L[3WfX L(|mM0

141水利论文M PH1crXwO

3.00水利论文8}~8}c'R

3

;N?Fvh:u)w0

1963水利论文A%P1xG#VOL0V(@_

80000水利论文+w V8L%_0u

92610

0o*oWV O0

487.9

+],oHy C)z0

132水利论文 tYL3Q5b](Mz

2.70

5W#z{+Wv0

4

fv-`AAX8in0St0

1964水利论文u t'b g7n8fK3m

91500

g'olK9q.L8^-?7R0AM0

142775水利论文g,Cv'o#w#h

572.3水利论文da7l%YT/m4iO!YU:o

182水利论文N;n1eOa[7`d5UO

3.37水利论文A `(EF(SUq0[j

5水利论文.yi;U#E E

1965水利论文5d;DS4E"ZZ%_

82500

9T+~N-Rp([.u0

2123160

:O z:^J8r0

455.7水利论文_+?sz*KlU8p"J{

113

h-X s#NA&D0

6.894水利论文 @f"]6O0R;B

6水利论文j4i#W&L{tB

1966水利论文w x*AcD$AV0c} S

20000水利论文2O{C@!Yn5b8Ta)F H

227800水利论文yr3{Q0M0L

481.3水利论文r'd5XI8BC7U4S$@8DL

170

Dz7AEH X0

5.03

s7r3v8kwx0

7

?NlbR6wy/M'rp kh0

1967

({*LV3Ng/G0

17800水利论文8_$T$Mw3l8O

140000水利论文!g#@y6i%YL7b

541.3水利论文#@'l[]!mI

199水利论文/X4S uCgk$hll

3.55

E"m L6Pn%^8~ y m0

8

:_n s[}3?o0

1968

s%u4n \,D:]o Ywk0

3900水利论文B:Z+bhv%u(dT,C

15980水利论文^HwG4m

538.6水利论文)qo F6E8@#~ @

186水利论文/f } W\hU

2.72水利论文&Q q4vO2SK/} ] osi

9

0M5hts2UN0

1969水利论文R'qv(V9k

17300

}Thv TY0

223200水利论文i"r.U2fs/\,_0Uu6y

460.6水利论文u7a%wN:{C)i

99.8水利论文oNJ(ZG;^F X

4.05

dt:J#Zh)y-CX0

10水利论文 m6cX&f{!D N0[j

1970水利论文2I9b/Z7f:\8u,M

25700

9k6y b&^)J/u6u$f8O{0

229400水利论文+j.k^R,F9|9I ^;pS)m

393.1水利论文d6FLk XN4CB?

133

}3c rf\0@eA|T0

3.22水利论文,] FKc;T;z/ov

11水利论文}3Nj d.Eq*A N

1971水利论文VXd3HevO^c8R

49400

#}^wO#HzK-s0

424500水利论文_G }'vY2C&f

373.9

m9\FxP$s7Y$y0

106水利论文 Nd ^g+hrS

2.65

^.FP0J&s8a E0

12水利论文"f"Vv#Rs.j

1972

sTd9kL8m8F5T pq0

40700

4f"e!a:M:op&C0

561700水利论文v ny7_c!SH#fn

428.8

1N@B1F-l5W+\0

107

J#[I wQ0]L0

1.91水利论文}&y5|lL

13水利论文:X9O&pG8{/Ub

1973

N(T x-g"Kdd0

77000水利论文d Lx4u5pY

563600水利论文9yi+~ VyB)m

482.1水利论文/h8K'pk.tj lr$E

140水利论文a+NLGx

3.00

;Gm(o~1{vr0

14水利论文'M }1S J/x]

1974水利论文 k+d1n]Ife K'o

72900水利论文WGAo4K3H's1q&A

557600

3h1`#qx(fZ0

415.1水利论文Hd6[:s/k

121水利论文%\2n Gg(},Pn)[S

1.31水利论文4m4F,SZ-G#N

15

8?#OJ)B_$R0

1975水利论文y+k)E1V.]]%T,N y

63300水利论文$bF5W7c:QK{]/{

528300

Y:k/N$fI.vKl0

536.7

:ycd:?p}'|X0

144

1x&| S whv0

2.33水利论文 ABOUbw4L

16水利论文$r }ro/} z$W

1976水利论文-c w)tNxM%g?

51600水利论文xV-VDW"T!\(tI

488940水利论文x#sk;DI:E6l

385.1水利论文xj`]2K _ K

154水利论文s'PK2R }^"|/f

3.55水利论文0R)JR8\.]

17水利论文 Y@}}&||E

1977水利论文:v qwN2^~I$E(e

60000

{GY2q&Sm EY c5~0

480500

/M6d;D8aGU*Q0

412.2

7SCB6nB@4x)A,u2H0

111水利论文,|7U%J5W#{c#@*d`

3.35水利论文u;J&^~%m%ih

18

FS X~ pZ"i q0

1978水利论文fB}%\2?R

70000

cK#A'?r.ZR0

530500水利论文"R5elo D?

567.1

)F8HZ/G)k"\? a0

139水利论文p?XIP]c^

2.55水利论文 r| Y5X?B%eV}


表2 年平均含沙量PPR模型拟合和预留检验结果

jq'[?2ya#K0o0

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

Yt1I;b+rHe `.WXp0

实测值

$b m FT3cL7VvIja0

计算值

5kb8A`zw2o+{0

绝对误差

-oqUIR1T Yq0

相对误差(%)

s3Ut9V!]M0

实测值

!?f3{mEH I8O0E0

计算值

x@!]h]\8^f0

绝对误差

^&H6q `&a(W!T r0

相对误差(%)水利论文1w l1lV(F g n.Y v2@


3.40水利论文W dl;Qe7k1V

3.131

6A*?[W)](wWfY^0

-0.087

E'@ FzjQkgf2nuP)~0

-2.5

:Yc-blY sAD1B0

2.72

s\9j,y$A(f/S2i(C0

2.934水利论文/x$c\G4Z

0.214

0qwZ-{_h0

7.9

h n8` A*DT.Fx:W,F%_,o0

3.00水利论文w0L(J3j$gP%R Gc

2.738水利论文|'w'UoJ%_U

-0.262水利论文9V(\@9r(TC/l

-8.7

^wpR-k(^ f0

4.05

.JBV#tms0

3.996水利论文!] E3o.D-]+v:[{ g*a

-0.054水利论文"^8|{mU5w.E{+Ql ?;R,V

-1.3水利论文 M J^G7X r6z8w

2.70

:`m$_M7e P| Q b0

2.659

u"T fC_W-@0

-0.041

wtvHac%_0

-1.5水利论文 KtH,Z'chf

3.22

$U DA~ b;oE*][l0

3.214水利论文 `A)qFL ]8sF

-0.006水利论文R[9X%R(p;i,v M'X

-0.2

Lk%s6Q,AIf0

3.37水利论文G}:]Xx

3.400

ch*~-j*|[&A(Yt0

0.030

%c'i-l+r h%P0

0.9水利论文b?+X#|'l$t

2.65

r}Due0

2.628水利论文:Ihq0lJ9X lm

0.022水利论文!k5^5K%qW`1G Pg

-0.8水利论文1Or*q)Z1C

6.894

vUfY)kvj0

6.926

6O'YZV&y0

0.032水利论文~ i^xgq

0.5水利论文#yO,Q[IF0jQu$]

1.91水利论文J)N%A:R~'e

2.135水利论文UL1YNO!YD

0.225

)a.BA,I t0

11.8水利论文0^l6@m f&t"[$@9g

5.03水利论文i/f7]C4J OI j/A o v

4.729水利论文y3o$nxl5^ p*~

-0.301

P'Mqt,S0V {?5ya;d+P0

-6.0水利论文8}-q EE6l:|.?

3.00水利论文!q-br6e `8d7`

2.987

r_J,}1nf0

-0.013水利论文-Kt)q ^l"Zz

-0.4水利论文n sM:U}&K;W

3.55

$q^{:e*L)a2WS0

3.826

WF&G {(A3A m0

0.276水利论文H1A;u|&VU!q

7.8

,I-iLK |%@#}*u0

1.31

"CR`(lC(K4f\0

1.319

@*T`v OR;} i,}0

0.009水利论文.X~x&O-j(W*s ~|.W)G~6G

0.7水利论文t+}*C)?/`~_F1S

2.33

:x,P iy6m4O$i0

2.461

R&oI^1W,M8u7FE s0

0.131水利论文(T C:L.Kw4{3Z `

5.6

#g] FhUs M!]0

3.35水利论文&T p k[hA

2.780水利论文-G.H(K\ `,ZN%[

-0.570水利论文@1Fh:q3v)M

-17.0

W#oNF6f_H2c0

3.55水利论文5QPlfG7FJ[2t

2.869

KB_gY0

-0.681

;r5Lo,Vn^,Z0

-19.2

i9a7g.Wj-L9J0

2.55水利论文1oQ a}7E

4.305

)n1XlC.eh Z!Kw?g0

1.755

7|Yo;p9L0

68.8水利论文0p"l&atCd&g


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

表3 各因子的相对权值水利论文v)b'Ea{

Relative weighted values of factors

/Cgl5pJ Ys0

权序水利论文!ax J p&t+R)T2k

1水利论文q4j;T5~-j Nc7`RsF

2

(n#_ ~ W*B.v0

3

!y1h"R]7?MZ^ K0

4

b0``_n7q0

因子

2Rj1e2B&O AW I0

X4水利论文t'e|,f}%e

X3

RL+j&Z1^pb `H0

X1水利论文j8e6H-g&{J1tr6`z

X2

+?+eH*H0|0

相对权值水利论文1f8E9_?Gw@3\r

1.00000

d(iv*bkhS`8~8q ?-q0

0.85190

3f't @bG0

0.71617

['TF.|/M\1VK0

0.29765水利论文 t;Vj;BR%g


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

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

8N@"d;Q)G0

4 结论水利论文kT7jB#w w,h

990112t1.gif (1658 bytes)水利论文.Kt+`1g(w7@%K#V

图1 某流域产沙量的PP回归拟合及预留检验效果水利论文O8WS'S$k;X

Effects of fitting and verification of sediment yield

6WA"xXiIS0

in a watershed by using PPR model水利论文`6D!N(o%f)~%Oe&g

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

\$mKoi1h5Cw0

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

B$QNS+A'R1vk0

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

8|+@6ut^n0

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

g2O,bN4h.`"~P.ZjZ0

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

+k[sx&^0

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

r+q&Us0gg9By%j0

参 考 文 献

6Z?XDL0

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

bHa(o9}+kG0

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

)@BWxqC3\],X0

3 Friedman J.H.&. Stuetzle W.,J. Amer. Statis. Assoc. 1981, 76:817.水利论文6\g~^%u w y%Z

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

c(avO c*s-w"}0

5 郑祖国。水文。1994,(4):6-10.水利论文~Om w{m

6 吴建明编著。环境统计学。北京:中国环境科学出版社,1991:504.水利论文\?.j7_4`;Wp2ef!P

 水利论文 n#m.KE |5r_

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

广告投放

广告投放