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

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

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

/@y}-Z!]mI'r0
李祚泳 邓新民侯宇光
(成都气象学院)(四川联合大学)

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

关键词流域 输沙量 投影寻踪回归 预测水利论文 q$Imv*AO)}

1 引言水利论文Q Of@ ];FP

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

D9l&Mi8`*`*v0

2 PPR原理及算法简介

L Q] ^1OO ^b u8r {0

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

a8bf~Q0

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

t1Zg]"z0_0

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

990112e1.gif (1513 bytes)水利论文)Su/K/MnTSb~ o|t

(1)

c3Sl5zz d8E&n4`0

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

Friedman和Stuetzle提出了实现PPR的SMART多重平滑回归技术,SMART模型具有如下形式水利论文5m$a(Ht&ed)O

990112e5.gif (1470 bytes)水利论文)JJ&[X2Gu

(2)

C(`#F6nY)o |M0

  它实际上是采用分层分组迭代交替优化方法对式(2)中的参数α,β,Mu和岭函数Gm寻优。实现步骤为水利论文C)w R;F/Vwb2W

  ①给定一个初始模型;水利论文W%kF.|U*U*_'w

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

*N9p eu f;xk0

     ③把上述投影中所包含的结构并到现有模型上,得到改进了的新模型;水利论文_7}R4\L C

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

     模型的关键是最终估计出式(1)的参数:岭函数最优化项数Mu,岭函数Gm,系数αjm,βm。其判别准则仍是:选择适当的参数组合,使式水利论文2[9g/L T{4c

990112e6.gif (1600 bytes)

iM(x+| T%\7^ x^0

(3)水利论文|7c-w?[

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

     用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]。水利论文4t+A1u"AfzB

3 某流域的年均含沙量的PPR回归预测水利论文5V&Z!Y6NccfgxiP

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

X&{R ] `%R&x6c:s0

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

     用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个因子对年均含沙量影响从大到小依次是年径流量>降雨量>采伐面积>采伐量。水利论文.L N$U&q3p)_7f EAo8D:@

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

xo?:whGbU0

Yearly average sediment concentration and measured data of related factors

iasD,c-Pc^R0

序号

7sE,k+W+qi0

年份水利论文@-L!N)y)?*zUb

采伐面积

xk.~1f0v*ll0

X1(M2

(c;|zJ|7J*l0D8g `1d0

采伐量水利论文Neh$u(r0Xt

X2(M3

(TeHxa9\7kY'g-gJ0

降雨量水利论文RT"s2Wi7I

X3(mm)水利论文,Zb$Q l R ~k

年平均流量

p'Jj?B}o!^M:o0

X4(m3/s)水利论文IM3o&_W;e k"TJ

年平均含沙量水利论文cU1C6OT.Q3t'N6V

Y(kg/m3水利论文 t3a#vHm$iSSe/iT


1

`'_#k+j-X#[4]0

1961

q nt#yQW0

15100

?{i)Z-Np0

29000水利论文{K(dyb \

508.0水利论文WjZH7z0Gor)v}|

180

k+f Zr\7l0

3.40

`/HT? ~l(Y)K2Z0

2

p^`!S:O| D#dt6?s7A*U0

1962水利论文2v GM+w'\_D#^{

20500

7EQ,c${!B~0

43350水利论文,YD#| q%U q4Sb7n5Z

453.7

:K N3f h ?2c K7oE`0

141水利论文^3ST ~R/Ph(Vr T

3.00

7N1Zy?iZB!]\8v(T0

3水利论文s2EV1SeceK#K&m%t p

1963

h,aFT9[/Jl7?0

80000

VnGW8nx]z~0

92610水利论文Y"H#K@Y-Z

487.9水利论文+Vw+U@!z0ph+]9A%r G9R

132

9X],Ct.iJ0

2.70

"rlc } YV0

4

G7r]3zQ0

1964

IP4U.bgZdV'X6g k0

91500水利论文 N;}%rg:?PCT&W(_N

142775

:c Y4dL_`0

572.3

y Bq-^bWG s0

182水利论文/^J+y7~&Y

3.37水利论文a OfqoG

5

){/R5x&@!\ F7H,R0

1965

1lx9V8N4N7`gG0

82500

bS+\O}2mZ:G C5R0

2123160

@Ed:YBu+Q0

455.7

9rf$P:S/nL0?X M2Fd0

113水利论文.Z#J5@/Bu

6.894水利论文)}+O:vQv7EW&X

6水利论文E5w+?$F yefg

1966

%@ QwO H9~0f_E0

20000水利论文6c*|@\J H6F

227800水利论文x N9X s#Zg

481.3

Nl H5?h0j0

170水利论文d7S$xf_q g

5.03

4gfs4LrKq-L+r x~|0

7水利论文~+z#~#H(F)^z

1967

DmDzrc s-b$K0

17800

$j@oCm ?v]0

140000水利论文X2mT0VM F"o.Q

541.3

9j;a!n:r#a}E Y0

199

$eC'Obx_G:fg0

3.55水利论文?wn4M,o$D%\Q

8水利论文}5B~B7Xw&N/cr]Y:j

1968

RS!x$I3Lo?sD0

3900

D/@c7|-g"t0

15980水利论文kFd+q}It

538.6水利论文F3A;OtEG

186水利论文-N[(ZD9^

2.72

)vCzqU6Xe-AcQS7C0

9水利论文v5~"J8yp7s l

1969水利论文"ei NM#P8W

17300水利论文T U ~T"x7a

223200

Yp \z"D8?wg0

460.6

B7gen-Me0

99.8

NS ?T4v%Ls3\0

4.05水利论文0s"i8RM#Q2b:K;s

10水利论文!I"Ep4S[4U

1970

RGBA1r |;Q0

25700水利论文-S3D:D`~I5^

229400

6d.}#s'AO)tc gV0

393.1水利论文X@Qoph

133水利论文 X%zwKWZ,Vny

3.22水利论文`;I6iZ O+R-{rlL

11水利论文Vk'{2|d @

1971

6qHU9G ATxO-S6}}0

49400水利论文IP ^9}5Sh R5n

424500水利论文 e!`@|bk

373.9

;J3U R0f&l8V&W0

106水利论文4~r&MfL Ei+VB

2.65

2^d],d-B)F u0

12水利论文(\5X)t&p2\1w"fP"t

1972水利论文-E0G7^1Wv\

40700水利论文 j(M;k C/};~%UjO(QS

561700水利论文@!N}5F8z'|3i

428.8

M@`xjCg zc0

107

1vv4aZ[B {G0

1.91

M7s6|$POW+BV0

13

R#\E Ri|7f:i%g0

1973

k#?'s%K4ux'C0

77000

1rU/B(G*k%x1f6u6j0

563600

.mzh(C+R0

482.1

P4OL h!{9OrI [0

140

$mOp4U3{M|f*C0

3.00

G_gK,yQ*?0

14

rEb"NvU)A0

1974水利论文)u3M/H)X"BI

72900水利论文AVM[+oVj"l

557600

[ ~8dImFW!y(R:Dv0

415.1水利论文 H:HjkPXn"D

121水利论文7i"R W C9\?}

1.31水利论文'Wo'@X7ODB)U8v

15水利论文hV y$y/I,n @?

1975

A4rB7eVf0

63300

+@wb7e6Y6Tc tk0

528300水利论文+AUZDc7hJ BQ h

536.7水利论文,b[+hM-pn$v

144

@6\v0j qhD0

2.33

`Q'm uae0

16水利论文ARg0eD^(v&q

1976水利论文FZFh z-w9N c;x0Wu

51600水利论文?R,U5EX#h2Lo!A"V

488940

,Q!\r4@,P[K/{0

385.1水利论文`y2E4z;u p NP)j

154

5Hab$u\0

3.55

U,f+XSgz0

17

WIfu*v:x0

1977水利论文 b)}?5w*Z5L3n{

60000水利论文9G |R~2o

480500

h6d I j9sYRZSW0

412.2

r.by!vc$d[0

111

[$u.Ga0FHe0

3.35水利论文 [6LQm8D9aAs

18

`R-DNY0

1978水利论文,pz5kX R_.c

70000

w*qt7KM0

530500

1~'D6G}(} Y!Y$]W \-q0

567.1

b+^T lK0

139

sPEnz Inn8B0

2.55

SG:d"H7QX0

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

|;bYuO(F0

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

&B.^G3F[/D0

实测值水利论文V U/?q3oq\q

计算值

A"P1r;]8a/c)WP|%N0

绝对误差

1Z5WlFp$p f0

相对误差(%)

^8RrS M s0

实测值

f$t1l.Xcny5t0

计算值水利论文G@g%t1t n+f#\2n1N

绝对误差水利论文FM T k?

相对误差(%)

@"B1}"MMC0

3.40水利论文\O8G6bF)rP3h l

3.131水利论文 b I AHFEI:~3v

-0.087

$uo3o`Z$]0

-2.5水利论文 ^8}] y"\0~8]#N)^E3c

2.72水利论文?;bJ.z*`djYwT

2.934水利论文 zm)w0h d1N2K'c M

0.214水利论文[D:`"v-l!sc&UW

7.9水利论文\q3E:]7h

3.00水利论文H;n8N~#q5fh5j9T4_

2.738

t2pm,qLj4gb0

-0.262

u9]:o9S:C;}#Q i"@0

-8.7

iWQ#Vr0

4.05

lH}@ r(g0

3.996

&Tu/y-|TE&D8tj:KD0

-0.054

:he(R2nUy9Y0

-1.3水利论文2^j1A3E!V

2.70

8oLso Ln q4v%W.n0

2.659

qgO/a \"m s2|0

-0.041水利论文t:t&A5|,vI-PI#i

-1.5水利论文@b;[|Yrbk"g

3.22

!XDJ;_}0I S lJ0

3.214

UO"K(}Z7}8z8}@f0

-0.006水利论文%`k#|#`|0M D~%UG

-0.2水利论文'a&j2V/e8Ex

3.37水利论文1o?7M\(|8?3Y\:n(X

3.400水利论文[T6J"H!i7p.J3SY nG0D

0.030

Mxrdg)X0

0.9水利论文 @q y Q.`!b%q$KM&W

2.65水利论文[R!u@;F:B@

2.628水利论文.~:TL p j)_ EdO8rrW

0.022水利论文 q9{ zME"J

-0.8

U'TFZ6^ A&\0

6.894

(Ms*FA5vt{y0

6.926

#]2Y@zw9C0

0.032水利论文 b`\ yH*N3J

0.5

G9QuVC+M*n0

1.91水利论文{5ok"vz-a^z Z

2.135

kHm$t.H0

0.225水利论文3{4X'n|&u

11.8

H{]&n Z7|P$TN0

5.03水利论文\I2dKH{u

4.729

3z fQZ'dW0

-0.301

;pnF0]0Tv J8v]3G0

-6.0水利论文o2Z}:IxwS#O

3.00

f!g*t9b_~/WL0

2.987水利论文My*\$YN W

-0.013

eWT0`UN0

-0.4水利论文_2M8h0T"l&V?"m

3.55

z r$m#P~S(f0

3.826

c,ODzYQ;aD[/Y Z0

0.276水利论文o.O/S-?_.A

7.8水利论文5{Xc3BD9V*^4m r

1.31

jW+Ye _sXs0

1.319

'J~9h4R9n:Y:S#^0

0.009

W+@9{{tP5^#}kh0

0.7

Pg$gh s(rz0

2.33水利论文A$Y8vm~ j|

2.461水利论文6jZw_g"\1Q

0.131

T4l[G!sp)_%G M0

5.6水利论文v/O.\z"Ge

3.35

j5j9i[n?0

2.780

Q2X3s,Bk:NK fS0

-0.570

s+Rk4wj&m0

-17.0水利论文{:vWi1h,x

3.55

B"`4lZN T0

2.869水利论文 F i-G)t.E Y

-0.681

-t M] AF8Hx~&Z0

-19.2

aL2^"~0J;mh6Iv0

2.55

$Xd9~nX H'C7l0

4.305水利论文q5AG P(C%Oa

1.755

J @{(wv9C7G0

68.8

K:?([[X S]N?A,}0

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

表3 各因子的相对权值水利论文e@]5`-O\

Relative weighted values of factors水利论文z[ @3fDdz


权序水利论文5c^&f;MQjBK z

1

'z6o2DJ5AK*jx0

2水利论文sz(Ij6CS(e,w"w

3水利论文FP jTJRv-h

4水利论文P'u t;w2@h)CA3k


因子

mZ1Y!K k@&a0

X4

@5D~[U4N0

X3水利论文pVMI1ea*h

X1水利论文C n*a.f(k

X2水利论文 _!{s z:]'N|9J9N S

相对权值水利论文)[G-Jm PF a Q

1.00000

JTb lQ dj`5o!N0

0.85190

J N pvs`/Y0

0.71617水利论文\4q4JG ^"Z6U"EW"U

0.29765水利论文5\ JbX?#A8Tjr


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

6P:g5X_\!l4d8z0

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

;TnTE1zB:ey D0

4 结论水利论文'qn7Gd0^?ZCO \

990112t1.gif (1658 bytes)

F r[.n8Q4t%V"u5K0

图1 某流域产沙量的PP回归拟合及预留检验效果

Y_T-U"u5s |y7u0

Effects of fitting and verification of sediment yield水利论文8e/u5\&vp/a/\ S

in a watershed by using PPR model

zg!] j2k x0

     1.PP回归采用了“审视数据—模拟—预测”探索性数据分析新途径,建立的模型稳健性和抗干扰性好,因而预测模型有较高的精度。水利论文0Bx@`u|2h

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

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

     4.PP回归也可用于水土流失中其它指标的预测建模。水利论文 R ~/BZ?D

   5.本文用PP回归技术建立了流域年均含沙量的预测模型,求得年均含沙量S后,只要将年径流量QT×年均含沙量S就可得到流域的年输沙量QsT。水利论文L*Z*i0u V/B)JnL

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

L-\,F8?z5~R{0

参 考 文 献水利论文.B/uFU/~)jB

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

[;}3nTTYJT0

2 张醒。泥沙研究。1991,(1):15-19.水利论文DE&EU/F^

3 Friedman J.H.&. Stuetzle W.,J. Amer. Statis. Assoc. 1981, 76:817.水利论文6E ]"b~H)r@)UI8r

4 Huber P.J..Ann. Statistics, 1985,13:435.水利论文8S X8h y{p2KE u p

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

6 吴建明编著。环境统计学。北京:中国环境科学出版社,1991:504.水利论文0AJ7U!j,hZ

 

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

广告投放

广告投放