CN109258041A - 一种幼年油茶林配方施肥方法 - Google Patents

Abstract

本发明公开了一种幼年油茶林配方施肥方法，属于配方施肥领域，是基于地形、气候和土壤多因子影响下的测土配方油茶林施肥方法；所述施肥方法是测定不同地形指标、气候指标和土壤指标以及相应的幼年油茶林生长量；采用主成分分析，德尔菲法和皮尔森相关分析法，筛选出与幼年油茶林生长量具有相关性的地形指标、气候指标和/或土壤指标构成最小数据集；根据所述最小数据集，建立养分效益模型，计算不同养分分配比例；通过数据包络分析养分利用效率；根据所述不同养分分配比例和养分利用效率，确定不同地形、气候和土壤条件下油茶林施肥量。本发明提出的幼年油茶林配方施肥方法，适用于不同地形、气候和土壤条件，从而使肥料投入合理化，提高肥料利用率和幼年油茶林的生长量。

Description

一种幼年油茶林配方施肥方法

技术领域

本发明涉及配方施肥领域，特别是指一种幼年油茶林配方施肥方法。

背景技术

油茶是世界四大木本油料之一，它生长于中国南方亚热带地区的高山及丘陵地带，是中国特有的一种纯天然高级油料。油茶作为高经济价值油料作物，其适应能力强，且不与农作物竞争耕地。

油茶幼林期是指从定植后到进入结果期前的阶段，一般为4-6年。油茶幼林期营养生长旺盛，促进春梢、夏梢的生长，根系和树冠迅速扩大，叶面积增加，培养良好的树体结构，促进树体养分积累，为开花结实和进入盛果期打下基础。油茶幼林期春梢、夏梢生长，根系和树冠扩大不仅受土壤肥力质量影响，同时受地形、气候和土壤共同影响。已有研究表明幼林期油茶春梢、夏梢的生长，根系和树冠迅速扩大与土壤肥力质量成正相关。

目前已有的施肥方案均基于一个试验地取样，只适用于实验开展地区，提高一个特定试验地的土壤肥力质量，而忽略了地形、气候和土壤多因子等多因素对油茶产量的影响，并不适用于不同地形和不同气候的油茶林。因此，现有施肥技术无法解决以湖南为例的地区，地形多样条件下幼年油茶抚育过程中养分投入效率低的现状。

发明内容

有鉴于此，本发明的目的在于提出一种幼年油茶林配方施肥方法，能够得出针对不同的地形、气候和土壤多因子影响下的配方施肥方法，提高幼年油茶林产量。

基于上述目的，本发明提供了一种幼年油茶林施肥方法，包括：测定不同地形指标、气候指标和土壤指标以及相应的幼年油茶林生长量；采用主成分分析，德尔菲法和皮尔森相关分析，筛选出与幼年油茶林生长量有相关性的地形指标、气候指标和/或土壤指标构成最小数据集；根据所述最小数据集，建立养分效益模型，计算不同养分分配比例；通过数据包络分析养分利用效率；根据所述不同养分分配比例和养分利用效率，确定不同地形、气候和土壤条件下幼年油茶林施肥量。

可选的，所述地形包括盆地、平原、丘陵或高山区。

可选的，所述幼年油茶林生长量包括幼年油茶春梢、夏梢生长量和树冠增长量；所述幼年油茶春梢、夏梢生长量包括春梢、夏梢平均长度；所述树冠增长量包括树冠面积增长量。

可选的，运用主成分分析法、德尔菲法和皮尔森相关分析法筛选出影响所述幼年油茶林生长量的影响因子的最小数据集。

可选的，所述最小数据集包括土壤容重、速效磷、有机质、日照量、土壤蓄水量和锌。

可选的，所述养分效益模型的计算公式如下：

幼年油茶林植株养分总需要量Y(kg/ha)＝优势株生物量增长量×养分系数；

幼年油茶林土壤养分有效性系数K(％)＝植株吸收量(kg/ha)/[土壤养分含量+养分供给量×养分利用率(％)]×100％；

土壤养分供应量：N(kg/ha)＝实测土壤养分含量×2.25×K(％)；

施用养分元素含量：C(kg/ha)＝(Y-N)/R。

可选的，所述养分利用效率通过数据包络分析，包括：综合技术效率、纯技术效率或规模效率。

可选的，所述综合技术效率＝纯技术效率×规模效率，所述综合技术效率由所述纯技术效率和所述规模效率共同决定。

可选的，所述纯技术效率和所述规模效率的值为1时，养分投入和所述幼年油茶林生长量效率达到最大化，所述规模效率值为1说明养分投入和所述幼年油茶林生长量达到最优状态。

从上面所述可以看出，本发明提供的一种幼年油茶林配方施肥方法与现有技术相比，通过建立针对不同地形，不同气候，不同油茶养分需求量的情况，引入养分吸收效率模型，根据评价的土壤指标与生长量的关系建立模型，合理投入肥料，提高肥料利用率和幼年油茶春梢、夏梢和树冠增长量。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例提供的幼年油茶林配方施肥方法的流程示意图；

图2为本发明实施例提供的幼年油茶林配方施肥方法的土壤样本位置示意图；

图3为本发明实施例提供的幼年油茶林配方施肥方法的不同施肥处理对树高增长量的影响结果示意图；

图4为本发明实施例提供的幼年油茶林配方施肥方法的不同施肥处理对冠幅生长量的影响结果示意图；

图5为本发明实施例提供的幼年油茶林配方施肥方法的不同施肥处理对春梢、夏梢生长量的影响结果示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明进一步详细说明。

本发明提供了一种幼年油茶林配方施肥方法的一种实施例，可适用于不同地形、气候和土壤条件。

参考图1所示，所述幼年油茶林配方施肥方法，包括：

步骤101：测定不同地形、气候和土壤指标以及相应的幼年油茶春梢、夏梢和树冠增长量；

步骤102：经相关性分析，得到不同地形、气候和土壤条件下与幼年油茶春梢、夏梢和树冠增长量的相关性，并筛选出最具有相关性的最小数据集；

步骤103：依据所筛选出的最小数据集，计算不同养分分配比例及施肥量；

步骤104：通过数据包络分析养分利用效率；

步骤105：确定不同地形、气候和土壤条件下油茶林施肥量。

从上述实施例可以看出，本发明提供的幼年油茶林配方施肥方法，通过筛选出与油茶产量相关的不同地形、气候的土壤指标最小数据集，并利用所述最小数据集计算相应参数，最后得到合适的施肥量，使得该施肥方法适用于不同地形、气候和土壤条件；该方法能充分考虑地形、气候和土壤多因子影响下有效提高幼年油茶春梢、夏梢和树冠增长量，减少化肥使用，有效防止林地土壤板结。

本发明还提供了所述幼年油茶林配方施肥方法的一种实施例。该幼年油茶林配方施肥方法包括如下步骤：

1.确定地形

参考图2所示，湖南省不同区域取盆地(东面罗霄山脉、西面雪峰山脉和南面南岭山脉包围的盆地区，以下称为region I)、平原(洞庭湖平原区，以下称为region II)，300米以下丘陵(南岭山脉区，以下为region III)、300米以上高山区(武陵山脉和雪峰山脉区，以下为region IV)四种不同的地形种植的生幼年油茶林。

2.取样

上述四种不同地形油茶林中每种地形选取试验样地30个，共计120个试验样地，所选取的试验样地在种植油茶林前均为生长了野草和灌木的废弃的次生林，采用机械翻耕土壤50cm。

试验选取的油茶林树龄为1-4年的幼林期，造林密度约为2m×2m。每个试验样地中选取3个能代表该样地特征的标准样地(20m×20m)。

选取的标准样地每个油茶林在土壤质地基本一致的林区，沿S型线路，设置6个采样点，均匀分布于样地中，同时避开路边、地角和堆积过肥料的地方，土钻取土深度为0～20cm，然后用四分法获得混合样。

土壤风干后全部过2mm筛、部分过0.149mm和1mm筛，用于测定土壤酶活性和理化性质的土壤样品指标。每个样地取环刀2～3个，设置一个土壤剖面，并进行土壤剖面调查。GPS采集每个标准样地海拔、坡度、坡向数据；年平均气温、年平均降雨量和年积温量由当地气象部门提供。

每个标准样地中选取平均株、优势株测定其生物量。同时采集植株叶片样本，分别在每株树冠的东、南、西和北面，植株内膛及顶部6个部位采集植树新梢叶和功能叶，每株取100片，共采集50株油茶树，混合叶片样本。采集的叶片样本保持新鲜，去除灰尘等脏物，于105℃杀青30min后，80℃烘干至恒重，将烘干处理后的叶片样本粉碎后保存，用于测定叶片样品营养元素含量。

3.测定

测定所述土壤样本指标及所述叶片样品营养元素含量。

(1)土壤容重、总孔隙度和蓄水量的测定：依据LY/T 1225-1999标准，采用环刀法；

所述LY/T 1225-1999为森林土壤颗粒组成(机械组成)的测定标准，其中LY/T1225-1999为标准号。

(2)土壤酸碱度(pH)的测定：依据LY/T1239-1999标准，采用电位法。

(3)土壤有机质的测定：依据LY/T1237-1999标准，采用重铬酸钾氧化—外加热法测定。

(4)土壤硝态氮(N)的测定：依据LY/T 1230-1999标准，采用酚二磺酸比色法。

(5)土壤全氮的测定：依据LY/T1228-1999标准，采用半微量凯氏法测定。

(6)森林土壤全钾(K)的测定：依据LY/T1254-1999标准，利用氢氟酸—高氯酸溶液消煮土壤样品，使土壤样品中的钾矿物分解成钙、镁、钾、钠、锰、铝等形成高氯酸盐类，用盐酸溶解残渣，成为可溶性的氯化物，制成钾的待测液，之后应用火焰光度计法测定土壤全钾量。

(7)森林土壤全磷(P)的测定：采用碱熔，在间断化学分析仪测定Smartchem 200(WestCo Scientific Instruments，Brookfield，CT，USA)。样品在银坩埚中用氢氧化钠高温熔融分解后，其中的不溶性磷酸盐转变成可溶性磷酸盐，待测液供测定全磷量。

(8)森林土壤有效硫(S)的测定：依据LY/T 1255-1999标准，采用燃烧碘量法和EDTA间接滴定法。

(9)森林土壤速效磷、速效钾、有效钙、有效镁、有效铁、有效锰、有效锌依据有效磷Mehlich 3方法提取，速效磷在间断化学分析仪测定Smartchem 200(WestCo ScientificInstruments，Brookfield，CT，USA)。速效钾在火焰光度计下测定；有效钙、有效镁、有效铁、有效锰、有效硼和有效锌在原子吸收分光光度计下测定；

上述Mehlich3试剂(简称M3)是适用于多类土壤、多种大量和微量养分元素测试的通用浸提剂。

(10)叶片中全氮测定：依据LY/T 1269-1999标准，采用凯氏法测定。

(11)叶片中全钾、磷、硫、锰、铜、锌测定：依据LY/T 1270-1999标准，采用硫酸(硫酸)--高氯酸消煮法；钙、镁、硼的测定采用ICP-AES法；锰、铁、锌和铜的测定采用原子吸收光谱测定；

所述ICP-AES为电感耦合等离子体原子发射光谱法，是以电感耦合等离子矩为激发光源的光谱分析方法，用于元素的测定。

如上所述，测定每个标准样地中取两年平均每年油茶鲜果重量和优势株生物量。

4.数据处理与结果分析

选取的湖南省四种不同地形中26个土壤、地形和气候指标方差分析表明，不同地形之间指标差异显著(p<0.01，见表1)。参考表1所述，盆地(Region I)中土壤指标铜、铁、钙、镁、速效磷、有效钾、有机质、硝态氮、全氮、总孔隙度、土壤蓄水量显著高于其他地区(p<0.01)；同时锰、容重显著低于其他地形(p<0.01)；平原区(Region II)硼、硫、全磷显著高于其他地形(p<0.01)；丘陵区(Region III)锌含量最高且具有显著性差异(p<0.01)；铁、有机质、全氮含量最低(p<0.01)；高山区(Region IV)钙、全磷、容重最高，而铜、硫、锌、镁、pH、速效磷、有效钾、硝态氮，全钾、总孔隙度、土层厚度、土壤蓄水量含量较低。

地形因子显示盆地多为缓坡和向阳面，平原地区海拔最低，高山地区海拔最高，同时坡度最大，向阴面面积也多于其他地区；气候因子显示平原地区降雨量和光照亮最丰富，而高山地区最少。

湖南省幼年油茶生长量冠幅增长量范围为32.45-52.32m²，。盆地地形幼年油茶林油茶产量最高(57.98m²)，随后是平原区(45.68m²)，低山丘陵区(39.32m²)，高山区(33.24m²)。偏相关分析产量与各因子之间的关系表明铜、锌、速效磷、有效钾、有机质、硝态氮、全氮、总孔隙度、土层厚度、土壤蓄水量与油茶生长量呈显著正相关，而钙、容重、海拔、坡向、温度和坡度与油茶生长量呈显著负相关。

表1土壤、地形和气候指标值及偏相关关系

**表示显著性小于0.01；

*表示显著性小于0.05；

小写字母表示Fisher's LSD多重比较不同样本之间显著性差异(p<0.05)，所述Fisher's LSD为最小显著性差异法；

n表示样本重复数；

方位:1＝南面；2＝东南面；3＝西南面；4＝东面；5＝西面；6＝东北面；7＝西北面；8＝北面。

对26个土壤、地形和气候相关因子进行主成分分析，运用spss软件算出6个主成分因子解释了大于82.31％累计率。运用德尔菲法(Delphi)进行30个专家评分及皮尔森相关分析确定最小数据集。各因子载荷在方差极大旋转后显示主成分1中权重前10％的指标有容重和速效磷，而Delphi法经过两轮问卷调查，其和谐系数为0.26，有显著统计意义(p＝0.01)，表明容重和速效磷均为油茶产量关键影响因子。主成分2显示有机质权重最高；主成分3显示日照权重最高；主成分4表明土壤蓄水量权重最高；主成分5表明锌权重最高。

因此，选取土壤容重、速效磷、有机质、日照量、土壤蓄水量和锌为衡量湖南省不同地形幼年油茶林生长量的最小数据集。

作为一种可选的实施方式，最小数据集中各指标得分采用隶属度函数(方程1-6)计算，对各指标进行标准化。其中，所述土壤有机质、有效磷和锌采用方程1计算，土壤容重采用方程2计算，日照量采用方程3，土壤蓄水量采用方程4。最小数据集中各指标隶属度函数式如下：

i)土壤有机质、速效磷和锌隶属度函数:

ii)土壤容重隶属度函数:

f(x)＝1/[1+17.18(x_x1-x)] 方程2

iii)日照量隶属度函数:

iv)土壤蓄水量隶属度函数:

其中，所述f(x)为土壤指标隶属度函数；x表示变量；x₀表示变量的最小值；x₁表示变量最大值。

从上述实施例可以看出，所有最小数据集中的变量通过层次分析确定权重。其中有机质权重最大，表明有机质对幼年油茶林生长量影响最为显著，随后排列为容重、速效磷、日照量、土壤蓄水量和锌。

依据方程5，计算不同地形油茶林质量指数。

其中，所述PQI表示油茶林质量指数，W_i表示指标的权重；Si表示指标分数；n表示最小数据集的个数。

结果显示不同地形油茶幼年林质量指数PQI值在0.54到0.82之间,其中盆地(Region I)的幼年林地质量优于其他地形。幼年油茶林质量等级排序如下：盆地(RegionI)，0.82±0.04；平原(Region II)，0.73±0.02；丘陵(Region III)，0.69±0.05；高山区(Region IV)，0.54±0.06。采用分段回归法建立油茶林质量指数与油茶生长量的关系，林地质量指数与生长量呈显著相关，方程描述如下:

Y₁＝96×10^-3x+0.2786(n＝360,r²＝0.758,p<0.05)

其中，所述Y₁表示冠幅增长量，x表示质量指数。

Y₂＝74×10^-3x+0.2064(n＝360,r²＝0.812,p<0.05)

其中，所述Y₂表示树高增长量，x表示质量指数。

Y₃＝14.6×10^-3x+0.1187(n＝360,r²＝0.792,p<0.05)

其中，所述Y₃表示春梢、夏梢增长量，x表示质量指数。

作为另一种可选的实施方式，采用灰色关联分析法分析26个土壤、地形和气候指标，其排序依次为：土壤有机质、速效磷、坡向、光照时间、土壤蓄水量、容重、锌、温度、土层厚度、、坡度、海拔、镁、降雨量、pH、值总孔隙度、全钾、硝态氮、全磷、钙、硫、铁、锰、硼、全氮和铜。其中与油茶生长量显著相关的前三个因子分别为有机质、速效磷、坡向。灰色关联系数显示幼年油茶林质量排序为盆地，平原，丘陵和高山区。

其中，所述灰色关联分析法是根据因素之间发展趋势的相似或相异程度，亦即“灰色关联度”，作为衡量因素间关联程度的一种方法。

从上述实施例可以看出，所述土壤因素对幼年油茶生长量增长的影响显著高于气候及地形指标对油茶生长量增长的影响。因此，施肥抚育是幼年油茶考虑的主要因素。

5.建立不同幼年油茶林质量养分管理模型

本发明的一个优选实施例为建立不同立地条件养分管理模型。

其中所述立地条件为与油茶生长发育有关的自然环境因子(如地形、土壤、水分等)，统称为它的立地条件。

本实施例中，所述养分管理模型，公式如下：

油茶林植株养分总需要量Y(kg/ha)＝优势株产量×养分系数；

油茶林土壤养分有效性系数K(％)＝植株吸收量(kg/ha)/[土壤养分含量+养分供给量×养分利用率(％)]×100％

土壤养分供应量：N(kg/ha)＝实测土壤养分含量×2.25×K(％)；

施用养分元素含量：C(kg/ha)＝(Y-N)/R。

(1)确定不同地形油茶林土壤有效养分校正系数

设立四种不同地形油茶林土壤硝态氮、有效磷和速效钾为自变量X_NI-IV、X_PI-IV和X_KI-IV，对应的土壤有效养分校正系数为因变量Y_NI-IV、Y_PI-IV和Y_KI-IV，通过回归分析筛选出的最佳数学模型均为双曲线模型。

模型1:

模型2:

模型3:

以上12个数学模型均达到极显著水平，对土壤有效养分校正系数的拟合程度较高，可用来计算土壤大量养分含量所对应的有效养分校正系数。

(2)确定实际施用养分元素含量

作为一种优选实施方式，分别以氮、磷、钾的实际施肥纯量(Y)为因变量Y_N、Y_P和Y_K，取目标产量(X₁)、土壤养分测试值(X₂)为自变量，建立二元一次回归方程：

盆地RegionI施肥模型如下

施N量，Y_NI＝-23.120+0.1123X₁-0.4865X₂(R²＝0.799)

施P量，Y_PI＝＝-36.24+0.0876X₁-0.2545X₂(R²＝0.813)

施K量，Y_KI＝-62.35+0.0721X₁-0.0224X₂(R²＝0.845)

平地RegionII施肥模型如下

施N量，Y_NII＝-37.652+0.0423X₁-0.1123X₂(R²＝0.822)

施P量，Y_PII＝＝-65.32+0.0645X₁-0.0845X₂(R²＝0.834)

施K量，Y_KII＝-33.65+0.0775X₁-0.0454X₂(R²＝0.889)

丘陵RegionIII施肥模型如下

施N量，Y_NIII＝-45.25+0.0985X₁-0.1982X₂(R²＝0.782)

施P量，Y_PIII＝-57.35+0.05456X₁-0.0654X₂(R²＝0.784)

施K量，Y_KIII＝-36.35+0.05214X₁-0.0325X₂(R²＝0.823)

高山区RegionIV施肥模型如下

施N量，Y_NIV＝-72.35+0.04582X₁-0.0546X₂(R²＝0.795)

施P量，Y_PIV＝-62.79+0.1951X₁-0.0627X₂(R²＝0.846)

施K量，Y_KIV＝-62.47+0.0247X₁-0.0728X₂(R²＝0.841)

将不同试验地的平均产量以及硝态氮、有效磷和速效钾平均值代入模型，计算出氮、磷、钾施肥量。

6.建立幼年油茶林养分效益模型计算不同养分分配比例及施肥量，通过数据包络分析(date envelopment analysis，DEA)分析养分利用效率。

构建肥料效率投入、油茶果产量指标，运用Deap 2.1软件，对四个不同地形的养分利用率进行测算，包括：综合技术效率、纯技术效率和规模效率。

本实施例中，模型中综合技术效率＝纯技术效率×规模效率，说明综合技术效率由纯技术效率和规模效率共同决定。根据表2所述，RegionI和RegionII两个地形纯技术效率和规模效率的值为1，说明这2个地形养分投入和幼年油茶生长量出效率达到最大化，而规模效率值为1说明养分投入和油茶生长量产出已达到最优状态，规模收益处于不变的状态。RegionIII和RegionIV两个地形纯技术效率值和规模效率值都分别为0.461和0.449，说明这两个地形的养分投入和油茶果生长量结构不合理，没有实现养分的优化配置，规模效率值不为1说明养分投入和油茶果产量没有达到最优状态，养分投入水平还需要调整；RegionIII规模报酬递增，表明适当增加投入量并合理的利用，会带来更高比例的产出。而其余RegionIV规模报酬递减，表明增加投入量已经不可能带来更大的产出，只会造成更严重的资源浪费。

油茶养分DEA效率为1的地区有RegionI和RegionII，这两个地形综合效率是处于DEA有效的状态；其余RegionIII和RegionIV的DEA值小于1，即处于DEA非有效的状态。(所述DEA分析法效率值为1则表示DEA有效，效率值不为1则表示DEA非有效)。

表2湖南省四种不同地形油茶林养分利用相对效率值

“—”表示规模报酬不变；irs表示规模收益递增；drs表示规模收益递减，且综合技术效率＝纯技术效率×规模效率。

根据表3所述，四个不同地形养分效率投入、产出指标松弛变量的测算结果可以看出，四个不同地形中RegionI和RegionII综合DEA效率为1，表示地形综合效率处于有效状态，不存在投入、产出松弛；RegionIII和RegionIV DEA值小于1，表示地形综合效率处于非有效状态，存在不同程度的投入冗余或是产出不足情况。其中，RegionII规模效率呈现递增的趋势，从输入指标看，氮存在投入冗余，冗余量分别为：12m²。虽然各个投入指标都有一定量的投入冗余，但是表现出一定的产出不足。说明RegionIII养分投入要素没有得到合理的分配，导致投入、产出不符，需要通过调整养分，从而提高效率。RegionIV同样存在投入冗余和产出不足同时存在的现状，但是却表现出规模效率递减的趋势，过多的投入不能带来更高比例的产出，此时要减少非有效的投入，将有效投入量合理利用转化为有效的产出。

表3湖南省幼年油茶林养分投入、产出指标的松弛变量均值

7.验证油茶成林专用肥施肥量

本实施例验证油茶成林专用肥施肥量的效果，即不同施肥后对幼年油茶林生长的影响。

从图3和表4中可以看出，油茶专用肥施肥后增长最大，较不施肥树高平均值增长了79.39％，较施用常规施肥树高平均值增长了13.89％。方差分析结果显示：不同施肥处理对树高影响差异显著(F＝27.672，p＝0.000)。

表4不同施肥处理对树高的影响

从图4和表5中可以看出油茶专用肥对油茶冠幅增长大于常规施肥和空白，其中专用肥冠幅增长较常规施肥提高27.47％，较空白提高51.09％，体现了不同肥料的特点。常规施肥为复混肥，特点是速溶养分比例固定，油茶并不能全部利用其中的养分，建议施用常规复混肥料需要采用坑施，少量多次施用，提高肥料利用率。油茶专用肥按照油茶生长需要的养分进行配比，适当水平的速效肥对冠幅增长有促进作用，在生长期施用可以明显提高冠幅的增长。不同施肥处理对冠幅影响差异显著(F＝93.621，p＝0.000)。

表5不同施肥处理对冠幅的影响

从图5和表6中可以看出油茶专用肥对油茶春梢、夏梢增长大于常规施肥和空白，其中专用肥冠幅增长较常规施肥提高25.10％，较空白提高42.47％，不同施肥处理对春梢、夏梢影响差异显著(F＝14.31，p＝0.000)。幼林主要以枝叶生长为主，成年后的油茶树主要依靠春梢、夏梢结果，春梢、夏梢生长越旺，发育就越良好，这是油茶增产的重要途径。

表6不同施肥处理对春梢、夏梢的影响

需要说明的是，本实施例中的常规施肥为复混肥，特点是速溶养分比例固定，油茶并不能全部利用其中的养分，建议施用常规复混肥料需要采用坑施，少量多次施用，提高肥料利用率。

所属领域的普通技术人员应当理解：以上任何实施例的讨论仅为示例性的，并非旨在暗示本公开的范围(包括权利要求)被限于这些例子；在本发明的思路下，以上实施例或者不同实施例中的技术特征之间也可以进行组合，步骤可以以任意顺序实现，并存在如上所述的本发明的不同方面的许多其它变化，为了简明它们没有在细节中提供。

本发明的实施例旨在涵盖落入所附权利要求的宽泛范围之内的所有这样的替换、修改和变型。因此，凡在本发明的精神和原则之内，所做的任何省略、修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种幼年油茶林配方施肥方法，其特征在于，包括：

测定不同地形指标、气候指标和土壤指标以及相应的幼年油茶林生长量；

采用主成分分析，德尔菲法和皮尔森相关分析法，筛选出与幼年油茶林生长量具有相关性的地形指标、气候指标和/或土壤指标构成最小数据集；

根据所述最小数据集，建立养分效益模型，计算不同养分分配比例；

通过数据包络分析养分利用效率；

根据所述不同养分分配比例和养分利用效率，确定不同地形、气候和土壤条件下油茶林施肥量。

2.根据权利要求1所述的一种幼年油茶林配方施肥方法，其特征在于，所述地形包括盆地、平原、丘陵或高山区。

3.根据权利要求1所述的一种有幼年油茶林配方施肥方法，其特征在于，所述幼年油茶林生长量包括幼年油茶春梢、夏梢生长量和树冠增长量；所述幼年油茶春梢、夏梢生长量包括春梢、夏梢平均长度；所述树冠增长量包括树冠面积增长量。

4.根据权利要求1所述的一种幼年油茶林配方施肥方法，其特征在于，运用主成分分析法、德尔菲法和皮尔森相关分析法筛选出影响所述油幼年油茶林生长量的影响因子的最小数据集。

5.根据权利要求1所述的一种幼年油茶林配方施肥方法，其特征在于，所述最小数据集包括土壤容重、速效磷、有机质、日照量、土壤蓄水量和锌。

6.根据权利要求1所述的一种幼年油茶林配方施肥方法，其特征在于，所述养分效益模型的计算公式如下：

幼年油茶林植株养分总需要量Y(kg/ha)＝优势株生物量增长量×养分系数；

幼年油茶林土壤养分有效性系数K(％)＝植株吸收量(kg/ha)/[土壤养分含量+养分供给量×养分利用率(％)]×100％；

土壤养分供应量：N(kg/ha)＝实测土壤养分含量×2.25×K(％)；

施用养分元素含量：C(kg/ha)＝(Y-N)/R。

7.根据权利要求1所述的一种幼年油茶林配方施肥方法，其特征在于，所述养分利用效率通过数据包络分析，包括：综合技术效率、纯技术效率或规模效率。

8.根据权利要求7所述的一种幼年油茶林配方施肥方法，其特征在于，所述综合技术效率＝纯技术效率×规模效率，所述综合技术效率由所述纯技术效率和所述规模效率共同决定。

9.根据权利要求8所述的一种幼年油茶林配方施肥方法，其特征在于，所述纯技术效率和所述规模效率的值为1时，养分投入和所述幼年油茶林生长量效率达到最大化，所述规模效率值为1说明养分投入和所述幼年油茶林生长量达到最优状态。