塔里木河流域胡杨生长模型研究

doi:10.13466/j.cnki.lyzygl.2020.05.015

摘要/Abstract

摘要：

以塔里木河流域78组胡杨树木为研究对象,进行胡杨树高、胸径、材积的生长模型研究。通过将建模和检验数据带入16个经验方程以及理论方程中,得到各模型参数值和决定系数(R²),依据平均绝对误差(MAE)、均方根误差(RMSE)、平均误差(ME)和平均预测误差百分比(MPSE),选出其中拟合效果最优的胡杨生长模型。最适宜胡杨树高的生长模型是柯列尔模型,其决定系数为0.771;最适宜胡杨胸径的生长模型也是柯列尔模型,其决定系数为0.969;最适宜胡杨材积的生长模型是逻辑斯蒂模型,其决定系数为0.966。通过设定的检验、评价标准,对多个拟合模型进行对比研究得到胡杨生长模型,以期为塔里木河流域胡杨的造林设计提供参考依据。

关键词: 塔里木河流域, 胡杨, 生长模型, 适应性检验

Abstract:

This paper takes 78 groups of Populus euphratica trees in the Tarim River Basin as the research object,aiming to study the growth model of P euphratica tree height,diameter at breast height and volume.The study introduces the modeling and testing data into 16 empirical equations and theoretical equations,and obtains each model parameter value and coefficient of determination (R ²).Based on the average absolute error (MAE),root mean square error (RMSE) and average error (ME ) and the average prediction error percentage (MPSE),the study selects the P euphratica growth model of the best fitting effect.It finds that the most suitable growth model for P euphratica tree height is the Collier model with a determination coefficient of 0.771,and the most suitable growth model for P euphratica at breast height is also the Kolier model with a determination coefficient of 0.969. The most suitable growth model for P euphratica volume is the logistic model and its coefficient of determination is 0.966.The growth model of P euphratica is obtained through comparative study of multiple fitting models through the set inspection and evaluation standards,in order to provide a reference for the afforestation design of P euphratica in the Tarim River Basin.

Key words: Tarim River Basin, Populus euphratica, growth model, fitness test

中图分类号:

S758
S792.11

王庆杰, 李国东, 刘晓曼, 安博文, 计怀峰. 塔里木河流域胡杨生长模型研究[J]. 林业资源管理, 2020,(5): 100-107.

WANG Qingjie, LI Guodong, LIU Xiaoman, AN Bowen, JI Huaifeng. Study on Populus euphratica Growth Model in Tarim River Basin[J]. FOREST RESOURCES WANAGEMENT, 2020,(5): 100-107.

图/表 12

表1

表2

经验方程和理论方程

编号	方程名称	表达式	编号	方程名称	表达式	编号	方程名称	表达式
(1)	柯列尔	y=at^be^-^ct	(7)	对数	y=a+blg(t)	(13)	坎派兹	y=a $e - b e - ct$ (c>0)
(2)	豪斯费尔德	y= $a (1 + b t - c)$	(8)	理查德	y=a(1-e^-^bx)^c	(14)	考尔夫	y=a $e - b t - c$
(3)	修正威布尔	y=a $1 - e - b t c$	(9)	双曲线	y=a- $b t + c$	(15)	莱瓦科威尔	y= $a (1 + b t - d) c$
(4)	吉田正男	y= $a 1 + b t - c$ +d	(10)	混合型	y= $(a + b t) - c$	(16)	逻辑斯蒂	y=a/(1+be^-^ct)
(5)	斯洛波达	y=a $e - b e - c t d$	(11)	舒马切尔	y=a $e - b t$
(6)	幂函数	y=at^b	(12)	单分子式	y=A(1-e^-^ct)(c>0)

表2

表3

表4

表5

表6

图1

图2

图3

图4

图5

图6

参考文献 21

[1]	Ayjamal K, Ümüt Halik, Maierdang K, et al. Gap dynamics of natural Populus euphratica floodplain forests affected by hydrological alteration along the Tarim River:implications for restoration of the riparian forests[J]. Forest Ecology and Management, 2019,438:103-113. doi: 10.1016/j.foreco.2019.02.009
[2]	Wang Dandan, Yu Zhitong, Peng Gang, et al. Water use strategies of Populus euphratica seedlings under groundwater fluctuation in the Tarim River Basin of Central Asia[J]. Catena, 2018,166:89-97. doi: 10.1016/j.catena.2018.03.020
[3]	张胜利, 代金玲, 姜海燕. 胡杨分子生物学研究进展[J]. 世界林业研究, 2019,32(1):22-27.
[4]	朱雅丽, 张绘芳, 张景路, 等. 基于无人机遥感的胡杨地上生物量估测模型构建[J]. 林业资源管理, 2019(2):80-87.
[5]	高生峰, 叶茂. 塔里木河上、中游极端干旱区淹灌对胡杨(Populus euphratica)长势的影响[J]. 生态科学, 2020,39(2):50-57.
[6]	周朝彬, 宋于洋, 王炳举, 等. 干旱胁迫对胡杨光合和叶绿素荧光参数的影响[J]. 西北林学院学报, 2009,24(4):5-9.
[7]	李熙萌, 杨琼, 李征珍, 等. 胡杨(Populus euphratica)叶片呼吸作用对地下水埋深的响应[J]. 生态科学, 2016,35(3):29-36.
[8]	陈元宏, 路静. 胡杨温室扦插育苗技术[J]. 辽宁农业科学, 2019(1):91-92.
[9]	Cailin Li, Shaowei Qin, Lianghong Bao, et al. Identification and functional prediction of circRNAs in Populus Euphratica Oliv.heteromorphic leaves[J]. Genomics, 2020,112(1):92-98. doi: 10.1016/j.ygeno.2019.01.013 pmid: 30707937
[10]	陈玉萍. 发扬“胡杨精神” 开创新疆高等教育发展新局面[J]. 新疆财经大学学报, 2019(1):5-8.
[11]	许业洲, 吴楚, 杜超群. 罗田垂枝杉单木生长规律及其生物量测定[J]. 东北林业大学学报, 2020,48(5):1-6.
[12]	孟得干, 于浩然, 刘剑, 等. 长白山区3种天然阔叶林分生长规律研究[J]. 林业资源管理, 2019(3):147-150.
[13]	曾伟生, 陈新云, 杨学云. 我国人工杨树生物量建模和生产力分析[J]. 林业科学, 2019,55(11):1-8.
[14]	韩金, 王新杰, 张鹏, 等. 不同立地条件下杉木径生长过程研究[J]. 山东农业大学学报:自然科学版, 2019,50(5):768-773.
[15]	薛春泉, 徐期瑚, 林丽平, 等. 基于异速生长和理论生长方程的广东省木荷生物量动态预测[J]. 林业科学, 2019,55(7):86-94.
[16]	付爱红, 陈亚宁, 李卫红. 塔里木河流域生态系统健康评价[J]. 生态学报, 2009,29(5):2418-2426.
[17]	胡春宏, 王延贵, 郭庆超, 等. 塔里木河干流河道演变与整治[M]. 北京: 科学出版社, 2011: 2-9.
[18]	孟宪宇. 测树学[M]. 北京: 中国林业出版社, 2011: 2-9.
[19]	袁慧, 杜超群, 李斌成, 等. 湖北省杉木生长模型研究[J]. 中国农学通报, 2019,35(34):50-58.
[20]	王佳慧, 董利虎, 李凤日. 帽儿山阔叶混交林天然更新幼苗幼树地径-树高模型[J]. 应用生态学报, 2019,30(11):3811-3823.
[21]	陈纪龙, 刘真真. 基于三维结构的胡杨茎生长模型的研究与可视化[J]. 塔里木大学学报, 2014,26(3):96-100.

变量	建模数据(n=41)				检验数据(n=37)
变量	最大值	最小值	均值	标准差	最大值	最小值	均值	标准差
树龄/a	102.00	5.00	30.93	23.46	45.00	8.00	21.59	8.270
树高/m	12.62	3.15	7.09	2.03	8.68	2.45	6.34	1.26
胸径/m	0.350	0.035	0.130	0.070	0.140	0.039	0.100	0.025
材积/m³	1.1715	0.0030	0.1602	0.2534	0.1344	0.003	0.0577	0.0325

模型	参数				模型检验指标
模型	a	b	c	d	R²	RMSE	MAE
柯列尔	2.737	0.272	-0.002		0.771	0.959	0.696
豪斯费尔德	13435164.603	5799842.965	0.342		0.766	0.969	0.687
莱瓦科威克	4322.05658	3608.18774	0.91951	0.37281	0.747	0.969	0.687
修正威布尔	342.735	0.007	0.346		0.766	0.970	0.687
吉田正男	12154069.084	8079928.468	0.413	1.245	0.767	0.968	0.692
斯洛波达	1539.39322	7.72629	0.2043	0.17359	0.745	0.967	0.695
幂函数	0.347	2.257			0.736	0.972	0.700
对数	2.380	-0.469			0.746	1.011	0.684
双曲线	18.34322	92.57468	1350.15059		0.742	0.993	0.713
混合型	-0.1069	0.46291	0.18909		0.757	0.964	0.697
单分子式	9.37788		0.06403		0.599	1.255	0.837
坎派兹	13.13569	1.11695	0.02083		0.732	1.073	0.752
考尔夫	43126.06768	0.03945	9.90784		0.753	0.973	0.684
舒马切尔	9.991	7.081			0.645	1.196	0.790
理查德	89.72211	2.31296E-5	0.34259		0.753	0.970	0.940
逻辑斯蒂	12.56733	1.79302	0.029		0.725	1.024	0.776

模型	参数				模型检验指标
模型	a	b	c	d	R²	RMSE	MAE
柯列尔	0.024	0.449	-0.006		0.969	0.012	0.008
豪斯费尔德	7655523.413	656158279.64	0.721		0.957	0.014	0.010
莱瓦科威克	8765.62717	29248.40462	1.31571	0.54798	0.953	0.014	0.011
修正威布尔	44.967	0.000	-0.722		0.957	0.148	0.010
吉田正男	1.243	77.309	0.714	-0.018	0.937	0.017	0.012
斯洛波达	315.97725	9.6861	0.05751	0.39088	0.964	0.012	0.010
幂函数	0.666	0.014			0.942	0.016	0.011
对数	0.087	-0.145			0.829	0.029	0.020
双曲线	5.58583	1778.00234	9863.0299		0.962	0.013	0.010
混合型	-5.8752	45.81237	0.35849		0.969	0.012	0.009
单分子式	0.46715		0.0114		0.926	0.019	0.016
坎派兹	0.55248	2.3396	0.01551		0.956	0.014	0.012
考尔夫	2.47778E7	0.03817	21.65695		0.952	0.015	0.012
舒马切尔	0.334	23.629			0.792	0.010	0.032
理查德	146.71653	2.05359E-6	0.72083		0.955	0.014	0.011
逻辑斯蒂	0.42517	6.12994	0.03134		0.949	0.015	0.012

模型	参数				模型检验指标
模型	a	b	c	d	R²	RMSE	MAE
柯列尔	0.00419	0.72777	-0.02245		0.965	0.045	0.026
豪斯费尔德	131490.7952	2.65646E9	2.16987		0.954	0.052	0.036
莱瓦科威克	9.08549E8	321.43643	5.30102	0.41634	0.953	0.05	0.036
修正威布尔	1	1	1		-0.484	—	—
吉田正男	71033.69166	160.52015	160.52015	0.31013	0.158	0.221	0.141
斯洛波达	759.04193	10.78983	0.01014	0.8484	0.964	0.045	0.028
幂函数	4.95382E-5	2.16969			0.955	0.087	0.041
对数	-0.628	0.248			0.721	0.173	0.122
双曲线	38.87394	3863.0469	150743.338		0.825	0.102	0.079
混合型	-6.82096	88.33189	0.52962		0.943	0.058	0.038
单分子式	11859.06412		5.79522E-7		0.709	0.133	0.104
坎派兹	46.2551	7.73845	0.00733		0.965	0.045	0.029
考尔夫	87426.6377	0.17601	25.4439		0.948	0.055	0.039
舒马切尔	4.3275	143.77232			0.914	0.073	0.057
理查德	141551.82202	4.42627E-5	2.17274		0.954	0.052	0.036
逻辑斯蒂	2.4785	97.85629	0.0439		0.966	0.045	0.030