Analysis of the factors affecting trunk density and wood density of three native tree species in Guangdong Province of southern China
-
摘要:目的分析不同因子对树干密度和木材密度的影响,为林木选育、碳汇计量提供数据支撑。方法基于广东省樟树、木荷、枫香3个乡土阔叶树种树干密度和木材密度的实测数据,利用含协变量和无交互作用的多因子方差分析法,从5大类30个因子(11个定性因子、19个定量因子)中,筛选出与树干密度和木材密度相关的因子,进而用增强回归树(BRT)来分析不同因子对3个树种树干密度和木材密度影响程度的大小。结果(1)植被类型、枝下高、胸径、植被总覆盖度、冠幅东西向宽度是影响樟树树干密度的主要因子,地市、植被类型是影响木荷树干密度的主要因子,坡向、海拔、平均高度为影响枫香树干密度的主要因子;3个树种树干密度的主要影响因子不同,无共同主要影响因子。(2)影响樟树木材密度的主要因子有枝下高、植被类型、海拔、植被总覆盖度、平均高度、灌木盖度、年龄、胸径、林种、土层厚度,影响木荷木材密度的主要因子有年龄、草本盖度、枝下高、平均胸径、土层厚度、植被类型,影响枫香木材密度的主要因子为坡向、海拔、平均高度、枝下高;3个树种木材密度具有共同主要影响因子枝下高,其相对贡献率相近,均在10%左右。(3)林分因子和单木因子同为影响樟树树干密度和木材密度的主导因子,其相对贡献率之和分别为87.04%和76.92%。林分因子、单木因子和地域因子是影响木荷树干密度的主导因子,其相对贡献率之和为79.96%;影响木荷木材密度的主导因子为林分因子、单木因子和土壤因子,其相对贡献率之和为83.04%。地形因子、林分因子和单木因子同是影响枫香树干密度和木材密度的主导因子,其相对贡献率之和分别为83.98%和92.70%。结论本文通过多因子方差分析和增强回归树对不同因子进行分析,得出林分因子和单木因子同是影响樟树、木荷、枫香树干密度和木材密度的主导因子。
-
关键词:
- 树干密度 /
- 木材密度 /
- 乡土树种 /
- 影响因子 /
- 增强回归树(BRT)
Abstract:Objective The effects of different factors on trunk density and wood density were analyzed to provide data support for tree breeding and carbon sequestration measurement.Method Based on the measured data of trunk density and wood density of Cinnamomum camphora, Schima superba and Liquidambar formosana, which are three native tree species in Guangdong Province of southern China, using multifactor analysis of variance with covariate and no interaction to screen out factors affecting trunk density and wood density from 30 factors (11 qualitative factors and 19 quantitative factors) of five categories, and then the boosted regression trees(BRT) was used to analyze the influence of different factors on trunk density and wood density of three species.Result (1) Vegetation type, height under branch, DBH, vegetation coverage and the crown width from east to west are the main factors affecting trunk density of Cinnamomum camphora. City and vegetation type are the main factors affecting trunk density of Schima superba. Slope aspect, altitude and average height are the main factors affecting trunk density of Liquidambar formosana. The main factors affecting trunk density of three tree species are different and there was no common main factors in them. (2) The main factors affecting wood density of Cinnamomum camphora are height under branch, vegetation type, altitude, vegetation coverage, average height, shrub coverage, age, DBH, forest category and soil thickness. The main factors affecting wood density of Schima superba are age, herb coverage, height under branch, average DBH, soil thickness and vegetation type. The main factors affecting wood density of Liquidambar formosana are slope aspect, altitude, average height, height under branch. The height under the branch is the common main influencing factor, affecting wood density of the three tree species, and the relative influence upon three tree species was similar, all of which was about 10%. (3) Stand factor and single tree factor are the dominant factors affecting trunk density and wood density of Cinnamomum camphora, and their total relative influence rates were 87.04% and 76.92%, respectively. Stand factor, single tree factor and regional factor are the dominant factors affecting trunk density of Schima superba, and the total relative influence rate was 79.96%. The dominant factors affecting wood density of Schima superba are stand factor, single tree factor and soil factor, and the total relative influence rate was 83.04%. Terrain factor, stand factors and single tree factor are the main factors affecting trunk density and wood density of Liquidambar formosana, and their total relative influence rates were 83.98% and 92.70%, respectively.Conclusion This paper analyzes different factors by multifactor analysis of variance and BRT, and it is concluded that stand factor and single tree factor are the dominant factors affecting trunk density and wood density of Cinnamomum camphora, Schima superba and Liquidambar formosana. -
在众多人造板材产品中,中密度纤维板(MDF)具有质地细腻、性能稳定及良好加工性能等特点,故被广泛地应用于地板、家具、建筑、室内装饰等领域[1-2]。对于MDF而言,磨削多用于去除热压成型后板材表面的预固化层以及表面装饰处理前降低表面粗糙度[3]。
在前期针对MDF磨削过程的实验研究中,多采用宏观层面的砂带磨削形式展开,重点考察磨削工艺参数对磨削力以及磨削表面质量的影响规律等[4-6]。当获得最佳磨削工艺参数以后,进一步从参数优化角度提高磨削过程的材料去除效率和加工表面质量就会变得困难。因此,从节能提效角度出发,考察磨具本身的切削特性具有重要的理论和实际价值。
磨削过程较一般切削过程更为复杂,在木质材料磨削加工中多采用砂带磨削形式,其本质上体现为众多磨粒与被加工材料之间的综合作用[7]。砂带表面众多形态各异的磨粒均可被近似视作刀尖圆弧半径为几微米到几十微米的一把微型切刀[8],当其进行磨削加工(多为负前角)时,在单磨粒的挤压、耕犁、切削作用和多磨粒协同干涉作用下木质材料会发生一系列的弹性变形—塑性变形—去除行为,最终形成一定量的材料去除以及新的加工表面。
一段时间以来,国内外一批学者针对金属、陶瓷等塑性材料和脆性材料积极开展了以单磨粒磨削为主要形式的研究[9-10]。针对木质材料的单磨粒研究甚少且进展缓慢[11],但对于具有自身材性特点的木质材料(如延展性差、非均质性、多孔性)而言,从单磨粒磨削角度对材料去除、表面创成等磨削理论进行完善尤显迫切。在金属、陶瓷等均质材料的磨削研究中,为了简化分析过程,在进行仿真和试验时多将单磨粒形态设定为圆锥形[12-13],但含明显切削刃特征的磨具对于具有纤维结构的木质材料来说具有更重要的去除作用。
采用球头圆锥磨粒和五棱锥磨粒,分别对MDF试件进行单磨粒磨削试验,并对其磨削特性进行比较分析,具体研究指标包括磨削表面形貌、成屑过程、磨削力、磨削力比以及能量消耗,以探究磨具切削刃在磨削MDF时的作用机制。本研究为实现MDF高效磨削提供理论依据,并对完善木质材料磨削理论具有重要意义。
1. 材料和方法
1.1 材 料
中密度纤维板(MDF),密度0.72 g/cm3,试件尺寸为130 mm (长) × 30 mm (宽) × 30 mm (厚);所有试件表面平整,无明显缺陷。在实验开始之前,为了准确测定磨削表面形貌参数,中密度纤维板在被锯制成规格试件前用150目砂带进行砂光预处理,所得表面粗糙度为2.843 μm。
1.2 试验装置
单磨粒磨削试验装置如图1所示。铝制圆盘(直径100 mm,厚度10 mm)通过锥形卡套和锁紧螺母装夹在高速精密电主轴(额定功率1.5 kW,设定转速2 400 r/min)上。该试验平台可沿X、Y、Z方向移动,调节精度均为0.01 mm。本试验采用楔形式划擦法,为了在试件表面产生具有渐变磨削深度的连续可辨划痕,将试件一端用铜箔垫起0.6 mm,并调节直线滑台使得MDF试件进给速度为9 m/min。磨削力测量系统主要包括压电晶体式测力仪(3257A,KISTLER,瑞士)、电荷放大器(5806,KISTLER,瑞士)以及信号分析仪(Omniace Ⅱ RA2300,NEC,日本)。在本研究中,测力仪的采样周期设定为50 μs,即采样频率为20 kHz。根据信号分析仪上的力学信号波形,调节圆盘Z方向位置以确定单磨粒顶端与试件前端的临界接触,实现精确对刀之后使圆盘在Z方向下降一高度,即为初始设定的磨粒切深。磨削动态过程通过高速摄像机(i-SPEED 3,OLYMPUS,日本)记录,本试验中设定拍摄速度为每秒15 000帧。单磨粒选取球头圆锥磨粒(尖端圆弧半径为0.02 mm,顶锥角为80°)和五棱锥磨粒(顶角为80°)两种形态,其材质均为304不锈钢,密度为7.93 g/cm3,其中五棱锥磨粒的切削形态(切削方向)见图1,两种磨粒均通过螺纹在铝盘的圆周面上进行刚性固结。
1.3 表面形貌表征
使用日本基恩士公司(KEYENCE)的VR-5000系列3D轮廓仪对MDF试件磨削表面进行三维扫描,为了避免对刀时磨粒对试件表面的刻划影响,试件前端形成的3道磨削痕迹不予测定。高度测量的重复精度为0.4 μm,单次视野成像(12 ×)范围为24 mm × 18 mm。随后使用配套分析软件VR Series(version 3.2.0.277)进行表面轮廓测定,获得磨削深度、沟壑体积等参数。
1.4 磨削力分解
磨削力是磨削过程中的一个重要参量,反映试件与磨粒间的接触干涉—相互作用状态。法向磨削力Fn主要影响磨粒与工件之间的接触变形、加工表面粗糙度以及工件表层的磨削质量;切向磨削力Ft则主要与磨削过程的磨粒与工件之间接触变形、功率消耗和磨粒磨损状态有关;一般情况下,轴向磨削力较小,可以不计[14]。在本试验中,通过压电晶体式测力仪测得的磨削力为水平方向上的磨削力Fh以及竖直方向上的磨削力Fv,如图2所示。根据力的分解,法向磨削力Fn和切向磨削力Ft的计算如下:
Ft=Fvsinα+Fhcosα (1) Fn=Fvcosα−Fhsinα (2) 式中:Ft为切向磨削力(N),Fn为法向磨削力(N),Fv为竖直方向上的磨削力(N),Fh为水平方向上的磨削力(N),
α 为磨粒回转轴线与竖直方向的夹角(°)。本研究中的磨粒切深非常小,理论磨削深度为0 ~ 0.3 mm,旋转角
α 的绝对值在3°左右。因此,近似认为单磨粒磨削过程中的法向磨削力Fn等于竖直方向磨削力Fv,切向磨削力Ft等于水平方向磨削力Fh[15]。采用力学–电压标定法[16],电压值U(V)与磨削力F(N)的转换见公式(3)。标定试验确定Fn的转换系数k为137.35,Ft的转换系数k为53.05。在本文中,磨削力选取的是最大磨削深度附近的最大磨削力,此时的旋转角α 等于0,Ft等于Fh,Fn等于Fv。定义旋转角α 等于0时,Ft与Fn的比值为磨削力比,用来表示切削阶段的占比,数值越大,表明切削能力越好[17]。F=kU (3) 1.5 比磨削能计算
比磨削能一般是指去除单位体积工件材料所消耗的能量,反映了磨粒与工件的干涉机理和干涉程度,直接关系到磨削效率和磨削质量,是评价磨削效果的重要指标[18-19]。比磨削能的大小与切向磨削力有密切关系,在本研究中,其计算公式如下:
S=EV (4) 式中:S 为比磨削能(J/cm3);V是单道划痕形成的沟壑体积(cm3),通过3D轮廓仪测得;E为材料去除过程中消耗的能量(J),表现为切向磨削力Ft做功,可通过OriginPro 9.1软件的积分功能实现(如图3所示),计算公式如下:
E=∫t2t1Ft(t)vcdt (5) 式中:
vc 为切削速度(m/s),可根据电主轴转速进行换算,本试验中相应的切削速度为13.82 m/s;t1和t2分别为磨粒切入和切出试件的时间点,其差值为Δt ,单位为s。于是,式(4)可以改写为:S=EV=kvc∫t2t1U(t)dtV (6) 2. 结果与分析
2.1 磨削表面形貌
图4为MDF试件的磨削表面三维形貌图,结合图5所示的磨削表面二维轮廓,发现在大致相同的磨粒切入深度下,五棱锥磨粒相比球头圆锥磨粒划擦时在沟壑两侧形成更多的材料隆起。由于五棱锥磨粒具有锋利的切削刃,在切削时呈现1个主切削刃、2个侧切削刃以及2个侧前刀面,而球头圆锥磨粒不具有明显的切削刃,这就造成了沟壑双侧边沿的形态差异。具体来说,五棱锥磨粒形成的沟壑两侧存在较多的片状材料去除(多为切削刃引起的材料崩裂),进而破坏了初始表面的完整性;球头圆锥磨粒形成的沟壑两侧除去局部材料隆起后,受相对光滑前刀面的连续挤压作用,临近沟壑边沿的初始表面完整性较好。分析二维轮廓,发现同金属、陶瓷等均质且延展性好的材料相比,MDF试件在磨削中形成材料隆起的连续性和显著性差,可见构成MDF的纤维−胶粒单元的延展流动特性弱。
2.2 磨屑生成特征
高速相机拍摄的动态磨削过程如图6、图7所示。由于磨粒的切削形态不同,选择具有相近沟壑体积而非相近磨削深度的两道磨削痕迹,来比较两种磨粒的成屑特性。图6中t1为磨粒切入时间点,t2为磨粒切出时间点。观察球头圆锥磨粒划擦MDF试件时的动态成屑过程,发现整体上产生的磨屑相对较少,磨屑流相对不明显,磨屑流多沿着切削方向的锥形前刀面流出,少量磨屑滑移至锥形后刀面附近,并且磨屑多集中产生在磨削过程的中后期。观察五棱锥磨粒划擦MDF试件时的动态成屑过程,发现整体上产生明显的磨屑(流),肉眼可见的磨屑生成的相对时间点相较球头圆锥磨粒更早,在磨粒切入试件初期便产生磨屑,磨屑流集中沿着2个平面形前刀面流出。
图 6 球头圆锥磨粒划擦MDF试件时的动态成屑过程t1为磨粒切入时间点,t2为磨粒切出时间点。此道划痕的最大磨削深度为0.180 mm,沟壑体积为0.593 mm3。t1 is the time point when the grit penetrates into workpiece and t2 is the time point when the grit penetrates out workpiece. Maximum cutting depth of this scratch is 0.180 mm, and the groove volume is 0.593 mm3.Figure 6. Dynamic chip formation process of MDF workpiece scratching by single spherical cone grit图 7 五棱锥磨粒划擦MDF试件时的动态成屑过程t1为磨粒切入时间点,t2为磨粒切出时间点。此道划痕的最大磨削深度为0.240 mm,沟壑体积为0.574 mm3。t1 is the time point when the grit penetrates into workpiece and t2 is the time point when the grit penetrates out workpiece. Maximum cutting depth of this scratch is 0.240 mm, and the groove volume is 0.574 mm3.Figure 7. Dynamic chip formation process of MDF workpiece scratching by single pentagonal pyramid grit为了进一步分析磨屑生成特征,对作为重要过程量并且对材料去除起重要作用的切向磨削力进行考察(图8)。切向磨削力做功的过程包含多种能量间的相互转化,集中反映在材料的变形、破坏直至去除。从图中可以看出,在约为1 ms的试件–磨粒接触干涉过程中(对应图8中竖直虚线),即从磨粒切入试件到切出试件,Ft曲线表现为先增大至一峰值,而后开始下降(两种情况均下降至负值),不难发现,五棱锥磨粒的Ft达到最大值的相对时点要早于球头圆锥磨粒。在此,不考虑磨粒与试件材料之间复杂的摩擦关系,且假定Ft做功均被用于材料的变形、破坏和去除。可以把单磨粒切削过程看成是一个能量“蓄积—释放”的循环过程,即便是“连续成屑”过程也是有微小的时间间隔的。分散看,伴随磨屑生成存在着许多小型能量“蓄积—释放”过程,在形成单道划痕时,五棱锥磨粒产生磨屑的相对时点明显早于球头圆锥磨粒,这表明五棱锥磨粒的切削去除性能优于球头圆锥磨粒,而球头圆锥磨粒在磨削初期Ft做功多转换为材料的塑性变形(如压密、隆起)而非材料分离;整体看,当接近磨削痕迹中间位置及理论最大磨削深度处时,是峰值磨削力所在位置,五棱锥磨粒的Ft达到最大值的相对时间点要早于球头圆锥磨粒,这进一步说明五棱锥磨粒的成屑(切削去除)能力优于球头圆锥磨粒。
2.3 磨削力与磨削力比
相关研究表明,磨削深度在诸多磨削工艺参数中对磨削力的影响显著。在本研究中,采用单因素试验方法,主要考察各评价指标随磨削深度的变化规律。图9和图10分别表示了磨削力与磨削深度之间的线性回归和二次多项式回归。由于单磨粒在划擦MDF试件过程中,呈现较大的负前角,所以整体上法向磨削力Fn大于切向磨削力Ft。回归分析结果显示,二次多项式回归的决定系数R2值大于线性回归,对因变量的变化具有更好的解释能力。根据二次多项式回归分析,发现对于球头圆锥磨粒和五棱锥磨粒而言无论是法向磨削力Fn还是切向磨削力Ft,与磨削深度之间均是二次幂增加趋势,这与之前的宏观砂带磨削实验结果相符[16]。对于球头圆锥磨粒来说,较大磨削深度时,磨削力随磨削深度增大而增大的趋势增强;对于五棱锥磨粒来说,较大磨削深度时,磨削力随磨削深度增大而增大的趋势减弱。根据能量“蓄积—释放”分析模型,推断当磨削深度较大时,五棱锥磨粒的切削能力更优,而圆锥磨粒由于没有显著的切削刃作用,Ft做功除去用于成屑外更多地转换为材料的塑性变形。
磨削力比随磨削深度的变化如图11所示,对于球头圆锥磨粒和五棱锥磨粒而言,并无明显的变化规律。球头圆锥磨粒切削时的磨削力比多集中在0.35 ~ 0.45之间,而五棱锥磨粒切削时的磨削力比多集中在0.5 ~ 0.8之间,这进一步说明了切削刃的存在更易于对纤维类材料进行磨削去除。
2.4 能量消耗
假定所产生的沟壑体积与被去除的材料体积相等,如图12所示,为两种磨粒切削形态下比磨削能随沟壑体积的变化情况。实际上,被去除的材料体积明显小于沟壑体积(存在着材料的塑性非断裂变形和断裂材料未被完整去除),故计算所得的比磨削能相较真实值偏小,但可以用来做比较分析。从图中可以得出,整体上两种磨粒切削时比磨削能与沟壑体积之间为二次幂关系,且呈现先降低后增大趋势,其中存在一个拐点,可称之为临界沟壑体积(对应一个临界磨削深度),对应磨削沟壑整体形态来看,此点位于沟壑最深处略偏向于磨粒划出位置。未达到临界沟壑体积时,整体比磨削能较大,因为此时磨削力沿竖直方向的平均分力向下,磨粒不断深入工件内部,消耗大量能量将材料推挤、分割和压缩,比磨削能的变化规律也符合磨削力的尺寸效应,即磨粒切深(或平均磨削面积)越小时比磨削能越大,切除单位体积材料需要更多的能量[18]。超越临界沟壑体积后,磨削力沿竖直方向的分力向上,磨粒沿磨削弧方向的分力会越来越小,而堆积在磨粒前刀面的材料也会随之被带出,此阶段磨粒的主要作用是去除材料。在沟壑体积基本相同的情况下,计算结果显示球头圆锥磨粒和五棱锥磨粒的比磨削能基本相当,但根据上面的分析知道五棱锥磨粒切削时的成屑能力更好,对于近似相同的沟壑体积产生更多的实际材料去除(磨屑),所以球头圆锥磨粒的真实比磨削能相较于五棱锥磨粒来说更大。据此推断,球头圆锥磨粒划擦试件时将更多的能量耗散在材料的挤压、耕犁等塑性变形上。
3. 结 论
采用球头圆锥磨粒和五棱锥磨粒对MDF试件分别开展磨削试验,并对相关磨削性能指标展开比较分析,得到如下结论:
(1)从磨削表面质量来看,五棱锥磨粒在磨削沟壑两侧产生较少的材料隆起,球头圆锥磨粒在相对光滑前刀面的连续挤压作用下使沟壑两侧边沿附近区域具有更好的初始表面完整性。
(2)从磨屑生成特征来看,当产生的沟壑体积相近时,五棱锥磨粒切削MDF试件生成的磨屑更多,磨屑流更为明显地沿着前刀面流出;匹配切向磨削力和动态磨削过程,一致得出五棱锥磨粒完成单道磨削痕迹时产生磨屑的相对时间点更早,从能量“蓄积—释放”的角度认为五棱锥磨粒的切削特性优于球头圆锥磨粒。
(3)回归分析显示磨削力与磨削深度之间呈二次幂增加关系;对于五棱锥磨粒而言,当磨削深度较大时,磨削力随磨削深度增大而增大的趋势减弱,而球头圆锥磨粒在磨削深度较大时磨削力随磨削深度增大而增大的趋势增强;五棱锥磨粒体现更大的磨削力比,进一步说明其切削性能更优。
(4)从能量消耗角度来看,当产生的沟壑体积相近时,五棱锥磨粒的比磨削能小于球头圆锥磨粒,据此推断球头圆锥磨粒划擦试件时将更多的能量耗散在材料的挤压、耕犁等塑性变形上,切削能力弱。
-
表 1 树干密度、木材密度统计量
Table 1 Statistics of trunk density and wood density
树种
Tree species项目
Item平均数
Mean/(g·cm− 3)标准差
Standard deviation/
(g·cm− 3)最小值
Min./(g·cm− 3)最大值
Max./(g·cm− 3)变异系数
Coeffient of
variation/%樟树
Cinnamomum camphora树干密度
Trunk density0.449 1 0.083 0 0.229 2 0.810 0 18.48 木材密度
Wood density0.447 7 0.084 7 0.230 4 0.723 0 18.92 木荷
Schima superba树干密度
Trunk density0.501 8 0.103 1 0.257 6 0.798 1 20.54 木材密度
Wood density0.503 9 0.109 4 0.257 1 0.781 5 21.70 枫香
Liquidambar formosana树干密度
Trunk density0.479 6 0.106 2 0.104 5 0.905 7 22.15 木材密度
Wood density0.469 8 0.109 8 0.097 9 0.871 8 23.38 表 2 定性因子统计表
Table 2 Statistical table of qualitative factors
因子类别 Factor category 因子 Factor 水平数 Level 樟树 Cinnamomum camphora 木荷 Schima superba 枫香 Liquidambar formosana 地形因子 Terrain factor 坡向 Slope aspect 9 9 8 坡位 Slope position 4 4 3 土壤因子 Soil factor 土壤类型 Soil type 4 3 3 林分因子 Stand factor 植被类型 Vegetation type 8 6 6 起源 Origin 3 2 2 优势树种 Dominant species 2 2 2 林种 Forest category 5 4 3 龄组 Age group 5 5 4 单木因子 Single tree factor 径阶 Diameter class 10 10 10 地域因子 Region factor 地市 City 17 20 12 气候带 Climate zone 3 2 2 表 3 树干密度多因子方差分析模型R2与调整R2
Table 3 R2 and adjusted R2 of multifactor analysis of variance for trunk density
树种 Tree species R2 调整R2 Adjusted R2 樟树 Cinnamomum camphora 0.518 0.273 木荷 Schima superba 0.659 0.502 枫香 Liquidambar formosana 0.632 0.445 表 4 树干密度多因子方差分析F值与P值统计
Table 4 F value and P value statistics of multifactor variance analysis for trunk density
因子类别
Factor category源
Source樟树
Cinnamomum camphora木荷
Schima superba枫香
Liquidambar formosanaF P > Pr F P > Pr F P > Pr 地形因子
Terrain factor海拔 Altitude − − 9.23 0.003 5** 43.64 < 0.000 1** 坡度 Slope − − 10.63 0.001 8** 3.03 0.087 1 坡位 Slope position 2.47 0.070 4 − − 2.75 0.072 4 坡向 Slope aspect 3.82 0.055 3 4.43 0.039 5* 1.27 0.280 8 土壤因子
Soil factor腐殖层厚度 Thickness of humic layer − − − − 9.30 0.003 4** 枯落物厚度 Litter depth 4.66 0.034 9* 4.49 0.038 1* 6.34 0.014 5* 土层厚度 Soil thickness − − 2.56 0.115 0 2.32 0.133 2 土壤类型 Soil type − − − − 5.94 0.004 4** 林分因子
Stand factor草本盖度 Herb coverage 1.23 0.271 1 3.01 0.087 6 8.23 0.005 7** 灌木盖度 Shrub coverage 1.33 0.254 0 − − 9.10 0.003 7** 林种 Forest category 2.80 0.033 8* 1.61 0.196 5 − − 龄组 Age group 1.67 0.183 1 5.87 0.000 4** 2.68 0.055 0 平均高度 Average height 2.36 0.130 0 6.11 0.016 2* 23.59 < 0.000 1** 起源 Origin − − 4.52 0.037 6* − − 郁闭度 Canopy density − − 1.23 0.272 6 2.98 0.089 7 植被类型 Vegetation type 2.40 0.038 3* 2.98 0.018 1* − − 植被总覆盖度 Total vegetation coverage 1.46 0.232 4 − − − − 单木因子
Single tree factor地径 Ground diameter 11.05 0.001 5** 3.26 0.076 0 − − 冠幅东西向 Crown width from east to west 1.86 0.177 6 3.82 0.055 3 − − 冠幅南北向 Crown width from north to south − − 9.51 0.003 0** − − 径阶 Diameter class − − − − 10.01 0.000 2** 年龄 Age − − 25.42 < 0.000 1** 6.99 0.010 5* 树高 Tree height − − 4.84 0.031 6* 2.96 0.090 6 胸径 DBH 8.28 0.005 6** − − 21.37 < 0.000 1** 枝下高 Height under branch 1.22 0.274 8 − − 2.38 0.127 9 地域因子
Region factor地市 City 7.84 0.006 9** 37.48 < 0.000 1** − − 气候带 Climate zone 1.54 0.222 4 6.70 0.012 1* 5.81 0.019 1* 注:−表示无影响的因子,*表示0.05水平显著因子,**表示0.01水平极显著因子。表6同此。Notes: − means non-influential factors,* means significant factor at 0.05 level,** means extremely significant factor at 0.01 level. Same as Tab. 6. 表 5 木材密度多因子方差分析模型R2与调整R2
Table 5 R2 and adjusted R2 of multifactor analysis of variance for wood density
树种 Tree species R2 调整 R2 Adjusted R2 樟树 Cinnamomum camphora 0.674 0.472 木荷 Schima superba 0.587 0.398 枫香 Liquidambar formosana 0.696 0.550 表 6 木材密度多因子方差分析F值与P值统计
Table 6 F value and P value statistics of multifactor variance analysis for tree wood density
因子类别
Factor category源
Source樟树
Cinnamomum camphora木荷
Schima superba枫香
Liquidambar formosanaF P > Pr F P > Pr F P > Pr 地形因子
Terrain factor海拔 Altitude 1.62 0.209 0 1.67 0.200 7 9.48 0.003 1** 坡度 Slope − − − − − − 坡位 Slope position − − 1.25 0.299 2 4.46 0.015 6* 坡向 Slope aspect − − − − 2.57 0.022 0** 土壤因子
Soil factor腐殖层厚度 Thickness of humic layer 3.62 0.062 2 − − 12.84 0.000 7** 枯落物厚度 Litter depth 2.21 0.143 1 − − 24.07 < 0.000 1** 土层厚度 Soil thickness 13.87 0.000 5** 1.78 0.187 6 − − 土壤类型 Soil type − − 5.25 0.007 9** − − 林分因子
Stand factor草本盖度 Herb coverage − − 7.28 0.009 0** 19.30 < 0.000 1** 灌木盖度 Shrub coverage 3.13 0.082 6 1.63 0.206 6 7.03 0.010 2* 林种 Forest category 5.17 0.001 3** 1.03 0.387 4 − − 龄组 Age group 4.24 0.009 1** 3.43 0.013 5* − − 平均高度 Average height 4.56 0.037 2* − − 8.98 0.003 9** 平均年龄 Average age 2.47 0.121 8 − − − − 平均胸径 Average DBH − − 2.37 0.129 1 − − 起源 Origin 5.20 0.026 4* 3.17 0.080 0 − − 郁闭度 Canopy density − − − − 2.74 0.103 2 优势树种 Dominant species 2.65 0.109 5 − − 6.40 0.014 1* 植被类型 Vegetation type 3.97 0.003 8** 1.57 0.182 1 2.78 0.025 1* 植被总覆盖度 Total vegetation coverage 7.13 0.001 0** − − 7.29 0.009 0** 单木因子
Single tree factor地径 Ground diameter 13.74 0.000 5** − − − − 冠幅东西向 Crown width from east to west − − 5.23 0.025 7* 18.40 < 0.000 1** 冠幅南北向 Crown width from north to south − − 3.06 0.085 2 − − 径阶 Diameter class 7.84 0.001 0** − − 22.07 < 0.000 1** 年龄 Age 6.12 0.016 4* 17.98 < 0.000 1** 11.72 0.001 1** 树高 Tree height − − − − 6.46 0.013 6* 胸径 DBH 8.70 0.004 7** − − − − 枝下高 Height under branch 6.05 0.017 1* 7.83 0.006 9** 2.39 0.127 2 地域因子
Region factor地市 City 12.34 < 0.000 1** 5.24 0.025 5* − − 气候带 Climate zone 4.72 0.012 8* − − − − 表 7 5大类因子对树干密度和木材密度影响的相对贡献率
Table 7 Relative contribution rate of five factors to trunk density and wood density
因子类别
Factor category树干密度 Trunk density 木材密度 Wood density 樟树
Cinnamomum camphora木荷
Schima
superba枫香
Liquidambar formosana樟树
Cinnamomum camphora木荷
Schima
superba枫香
Liquidambar formosana地形因子 Terrain factor 6.08 12.44 36.15 7.95 12.03 32.23 土壤因子 Soil factor 3.66 7.60 15.12 13.02 13.63 7.30 林分因子 Stand factor 53.43 32.75 23.62 49.62 36.92 34.62 单木因子 Single tree factor 33.61 32.02 24.20 27.30 32.79 25.85 地域因子 Region factor 3.21 15.19 0.90 2.11 4.62 0.00 -
[1] Henry M, Besnard A, Asante W A, et al. Wood density, phytomass variations within and among trees, and allometric equations in a tropical rainforest of Africa[J]. Forest Ecology and Management, 2010, 260(8): 1375−1388. doi: 10.1016/j.foreco.2010.07.040
[2] 董利虎, 李凤日, 贾炜玮. 黑龙江省红松人工林立木生物量估算模型的研建[J]. 北京林业大学学报, 2012, 34(6):16−22. Dong L H, Li F R, Jia W W. Effects of tree competition on biomass and biomass models of Pinus koraiensis plantation[J]. Journal of Beijing Forestry University, 2012, 34(6): 16−22.
[3] 李海奎, 赵鹏祥, 雷渊才, 等. 基于森林清查资料的乔木林生物量估算方法的比较[J]. 林业科学, 2012, 48(5):44−52. doi: 10.11707/j.1001-7488.20120507 Li H K, Zhao P X, Lei Y C, et al. Comparison on estimation of wood biomass using forest inventory data[J]. Scientia Silvae Sinicae, 2012, 48(5): 44−52. doi: 10.11707/j.1001-7488.20120507
[4] 符利勇, 雷渊才, 曾伟生. 几种相容性生物量模型及估计方法的比较[J]. 林业科学, 2014, 50(6):42−54. Fu L Y, Lei Y C, Zeng W S. Comparison of several compatible biomass models and estimation approaches[J]. Scientia Silvae Sinicae, 2014, 50(6): 42−54.
[5] 曾伟生. 基于木材密度的34个树种组一元立木生物量模型建立[J]. 林业资源管理, 2017(6):41−46. Zeng W S. Developing one-variable individual tree biomass models based on wood density for 34 tree species in China[J]. Forest Resources Management, 2017(6): 41−46.
[6] 张会儒, 唐守正, 胥辉. 树干平均密度的估计方法研究[J]. 林业科学研究, 1998, 11(1):58−62. Zhang H R, Tang S Z, Xu H. Studies on the estimate method for average stem density[J]. Forest Research, 1998, 11(1): 58−62.
[7] He D, Deane D C. The relationship between trunk- and twigwood density shifts with tree size and species stature[J]. Forest Ecology and Management, 2016, 372: 137−142. doi: 10.1016/j.foreco.2016.04.015
[8] Dadzie P K, Amoah M, Boampong E, et al. Effect of density and moisture content on biological durability of stem and branch wood of Entandrophragma cylindricum, (sapele)[J]. Journal of the Indian Academy of Wood Science, 2015, 12(1): 44−53. doi: 10.1007/s13196-015-0143-5
[9] Wu R. Microstructural study of sanded and polished wood by replication[J]. Wood Science and Technology, 1998, 32(4): 247−260. doi: 10.1007/BF00702893
[10] Machado J S, Louzada J L, Santos A J A, et al. Variation of wood density and mechanical properties of blackwood (Acacia melanoxylon R. Br.)[J]. Materials & Design, 2014, 56(4): 975−980.
[11] Molteberg D, Hoibo O. Modelling of wood density and fibre dimensions in mature Norway spruce[J]. Canadian Journal of Forest Research, 2007, 37(37): 1373−1389.
[12] Poorter L, Wright S J, Paz H, et al. Are functional traits good predictors of demographic rates? Evidence from five neotropical forests[J]. Ecology, 2008, 89(7): 1908−1920. doi: 10.1890/07-0207.1
[13] 尹惠妍. 大区域森林生物量估算方法研究[D]. 北京: 中国林业科学研究院, 2014. Yin H Y. Study on estimation method to regional forest biomass[D]. Beijing: Chinese Academy of Forestry, 2014.
[14] 王仲锋. 森林生物量建模与精度分析[D]. 北京: 北京林业大学, 2006. Wang Z F. On the forest biomass's modeling and precision analysis[D]. Beijing: Beijing Forestry University, 2006.
[15] 罗云建. 华北落叶松人工林生物量碳计量参数研究[D]. 北京: 中国林业科学研究院, 2007. Luo Y J. Study on biomass carbon accounting factors of Larix principis-rupprechtii plantation[D]. Beijing: Chinese Academy of Forestry, 2007.
[16] Patino S, Lloyd J, Paiva R, et al. Branch xylem density variations across the Amazon Basin[J]. Biogeosciences, 2009, 6(4): 545−568. doi: 10.5194/bg-6-545-2009
[17] 徐金梅, 吕建雄, 鲍甫成, 等. 祁连山青海云杉木材密度对气候变化的响应[J]. 北京林业大学学报, 2011, 33(5):115−121. Xu J M, Lü J X, Bao F C, et al. Response of wood density of Picea crassifolia to climate change in Qilian Mountains of northwestern China[J]. Journal of Beijing Forestry University, 2011, 33(5): 115−121.
[18] 徐明锋, 柯娴氡, 张毅, 等. 粤东6种阔叶树木材密度及其影响因子研究[J]. 华南农业大学学报, 2016, 37(3):100−106. doi: 10.7671/j.issn.1001-411X.2016.03.016 Xu M F, Ke X D, Zhang Y, et al. Wood densities of six hardwood tree species in Eastern Guangdong and influencing factors[J]. Journal of South China Agricultural University, 2016, 37(3): 100−106. doi: 10.7671/j.issn.1001-411X.2016.03.016
[19] Koprowski M, Duncker P. Tree ring width and wood density as the indicators of climatic factors and insect outbreaks affecting spruce growth[J]. Ecological Indicators, 2012, 23(4): 332−337.
[20] 朱丽梅, 胥辉. 思茅松树干平均密度估计方法研究[J]. 西南林业大学学报, 2008, 28(2):17−20. doi: 10.3969/j.issn.2095-1914.2008.02.005 Zhu L M, Xu H. Studies on estimation methods for average stem density of Pinus kesiya var. langbianensis[J]. Journal of Southwest Forestry University, 2008, 28(2): 17−20. doi: 10.3969/j.issn.2095-1914.2008.02.005
[21] Nguyen H, Firn J, Lamb D, et al. Wood density: a tool to find complementary species for the design of mixed species plantations[J]. Forest Ecology and Management, 2014, 334: 106−113. doi: 10.1016/j.foreco.2014.08.022
[22] 王秀花, 陈柳英, 马丽珍, 等. 7年生木荷生长和木材基本密度地理遗传变异及种源选择[J]. 林业科学研究, 2011, 24(3):307−313. Wang X H, Chen L Y, Ma L Z, et al. Geographical provenance variation of growth and wood basic density of 7-year-old Schima superba and its provenance selection[J]. Forest Research, 2011, 24(3): 307−313.
[23] 黄继丰, 冯学周, 莫罗坚, 等. 粤北10种乡土阔叶树种生长及根际养分比较[J]. 林业与环境科学, 2006, 22(3):17−21. doi: 10.3969/j.issn.1006-4427.2006.03.004 Huang J F, Feng X Z, Mo L J, et al. Growth and rhizospheic fertility of 10 native broad-leaved trees in north Guangdong[J]. Forestry and Environmental Science, 2006, 22(3): 17−21. doi: 10.3969/j.issn.1006-4427.2006.03.004
[24] Franklin J, McCullough P, Gray C. Terrain variables used for predictive mapping of vegetation communities in Southern California[M]//Wilson J P, Gallant J C. Terrain analysis: principles and applications. New York: John Wiley and Sons, 2000: 331−353.
[25] De'ath G. Boosted trees for ecological modeling and prediction[J]. Ecology, 2007, 88(1): 243−251. doi: 10.1890/0012-9658(2007)88[243:BTFEMA]2.0.CO;2
[26] Elith J, Leathwick J R, Hastie T. A working guide to boosted regression trees[J]. Journal of Animal Ecology, 2008, 77(4): 802−813. doi: 10.1111/j.1365-2656.2008.01390.x
[27] Weisberg P J, Shandra O, Becker M E. Landscape influences on recent timberline shifts in the Carpathian Mountains: abiotic influences modulate effects of land-use change[J]. Arctic, Antarctic, and Alpine Research, 2013, 45(3): 404−414. doi: 10.1657/1938-4246-45.3.404
[28] Pittman S J, Costa B M, Battista T A. Using lidar bathymetry and boosted regression trees to predict the diversity and abundance of fish and corals[J]. Journal of Coastal Research, 2009, 25(6): 27−38.
[29] 焦琳琳, 常禹, 申丹, 等. 利用增强回归树分析中国野火空间分布格局的影响因素[J]. 生态学杂志, 2015, 34(8):2288−2296. Jiao L L, Chang Y, Shen D, et al. Using boosted regression trees to analyze the factors affecting the spatial distribution pattern of wildfire in China[J]. Chinese Journal of Ecology, 2015, 34(8): 2288−2296.
[30] 李春林, 刘淼, 胡远满, 等. 基于增强回归树和Logistic回归的城市扩展驱动力分析[J]. 生态学报, 2014, 34(3):727−737. Li C L, Liu M, Hu Y M, et al. Driving forces analysis of urban expansion based on boosted regression trees and Logistic regression[J]. Acta Ecologica Sinica, 2014, 34(3): 727−737.
[31] 刘一星, 赵广杰. 木材学[M]. 北京: 中国林业出版社, 2012. Liu Y X, Zhao G J. Wood science[M]. Beijing: China Forestry Publishing House, 2012.
[32] 张大鹏. 不同海拔高度对马尾松茎木质部结构的影响[J]. 生态学杂志, 1994, 13(2):12−18. Zhang D P. Xylem structure of Pinus massonions stem at different altitudes[J]. Chinese Journal of Ecology, 1994, 13(2): 12−18.
[33] 刘彦春, 张远东, 刘世荣, 等. 川西亚高山针阔混交林乔木层生物量、生产力随海拔梯度的变化[J]. 生态学报, 2010, 30(21):5810−5820. Liu Y C, Zhang Y D, Liu S R, et al. Changes of tree layer aboveground biomass, ANPP to altitudinal gradient in the subalpine secondary mixed forest of Western Sichuan, China[J]. Acta Ecologica Sinica, 2010, 30(21): 5810−5820.
[34] 国家林业局. 立木生物量模型及碳计量参数: 木荷: LY/T 2660—2016[S]. 北京: 中国标准出版社, 2017. State Forestry Administration. Tree biomass models and related parameters to carbon accounting for Schima superba: LY/T 2660—2016[S]. Beijing: Standards Press of China, 2017.
[35] 国家林业局. 立木生物量模型及碳计量参数: 枫香: LY/T 2661—2016[S]. 北京: 中国标准出版社, 2017. State Forestry Administration. Tree biomass models and related parameters to carbon accounting for Liquidambar formosana: LY/T 2661—2016[S]. Beijing: Standards Press of China, 2017.
[36] 闫晶, 罗云建, 郑德福, 等. 桉树生物量估算差异的源解析[J]. 林业科学, 2014, 50(2):92−98. Yan J, Luo Y J, Zheng D F, et al. Source appointment of differences in biomass estimates of eucalypt plantation[J]. Scientia Silvae Sinicae, 2014, 50(2): 92−98.
-
期刊类型引用(1)
1. 杨冲,张子木,朱麒轩,杜瑶,李黎,罗斌. 基于中密度纤维板磨削的定组态砂带寿命研究. 林业机械与木工设备. 2023(03): 62-66+75 . 百度学术
其他类型引用(4)