CN110009009B - 一种基于主成分计算的绿化树种选择方法 - Google Patents
一种基于主成分计算的绿化树种选择方法 Download PDFInfo
- Publication number
- CN110009009B CN110009009B CN201910208015.XA CN201910208015A CN110009009B CN 110009009 B CN110009009 B CN 110009009B CN 201910208015 A CN201910208015 A CN 201910208015A CN 110009009 B CN110009009 B CN 110009009B
- Authority
- CN
- China
- Prior art keywords
- value
- tree species
- ecological
- unit
- forest
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/10—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
- Y02A40/22—Improving land use; Improving water use or availability; Controlling erosion
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P60/00—Technologies relating to agriculture, livestock or agroalimentary industries
- Y02P60/40—Afforestation or reforestation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cultivation Of Plants (AREA)
Abstract
本发明涉及基于主成分计算的绿化树种选择方法，可有效解决树种的选择，保证绿化质量的问题，方法是，对不同树种的生态适应性与生态功能进行测定，以这些树种生态适应性总指标值，气体调节、土壤保持、水调节、营养积累、生产木材的价值为依据，采用主成份分析法，对各树种以上指标进行主成分测定，计算树种综合指标值，对造林树种进行量化选择，本发明方法新颖独特，科学合理，易操作，可有效用于绿化造林树种的选择，具有很强的实际应用价值，为实现绿化造林提供了技术支撑，有显著的经济和社会效益。
Description
技术领域
本发明涉及绿化树种，特别是用于土壤气候条件相近的平原区绿化树种的选择，也适用于立地条件相近的丘陵山区绿化树种选择的一种基于主成分计算的绿化树种选择方法。
背景技术
树种选择关系到造林工作的成败，如何筛选绿化树种值得深入探讨。国外通过制定树种选择手册或树种特性表，对树种进行选择。国内树种选择主要依据“适地适树”原则，比较宏观，在实际造林工作中难以贯彻执行，树种选择的基础研究较少，成果不系统，对树种选择的支撑能力较差。查阅国内外树种选择的相关文献，尚未建立系统的树种量化选择方法。
森林是陆地生态系统的主体，在全球范围内发挥着重要的作用。森林生态功能及其价值是衡量森林生态系统质量的一个重要指标。选择生态适应性强、生态服务价值较高的树种造林，可以提高造林成活率及其综合效益。树种的生态适应性主要是指其抗寒、抗旱、抗风、抗病虫害情况，生态适应性较强的树种一般生长迅速，具有较高的生态价值。但是不同树种的生物学特性差别较大，在适应环境之后，生态价值差异也较大。森林中乔木层是影响生态系统功能的最重要层次，其生态价值主要表现在气体调节、土壤保持、水调节、营养积累等方面。
本方法以平原区主要绿化树种为研究对象，对不同树种的生态适应性与生态服务功能进行测定，以这些树种生态适应性总指标值，气体调节、土壤保持、水调节、营养积累、生产木材的价值为依据，采用主成份分析法，计算各树种综合指标值，对造林树种进行量化选择，以保证绿化质量，但至今未见有公开报导。
发明内容
针对上述情况，为克服现有技术之缺陷，本发明之目的就是提供一种基于主成分计算的绿化树种选择方法，可有效解决树种的选择，保证绿化质量的问题。
本发明解决的技术方案是，一种基于主成分计算的绿化树种选择方法是，对不同树种的生态适应性与生态功能进行测定，以这些树种生态适应性总指标值，气体调节、土壤保持、水调节、营养积累、生产木材的价值为依据，采用主成份分析法，对各树种以上指标进行主成分测定，计算树种综合指标值，对造林树种进行量化选择，具体包括以下步骤：
(1)建立树种选择指标体系：根据树种的生态适应性，对树种进行筛选，在树种适应环境条件下，才能够提供较高的生态价值，在构建树种选择指标体系时，选择了生态适应性与生态价值指标；
(2)设置标准地：在土壤气候条件相近的苗圃设为标准地，作为树种生态适应性与生态功能测定的场地；
(3)生态适应性指标值测定与计算：在标准地中，观测记录不同树种抗寒、抗旱、抗病虫害、抗风情况，对以上4个指标进行测定分级量化，然后合计得出不同树种生态适应性总指标值；
(4)气体调节价值测定与计算：在标准地中，用收获法测算林木地上与地下部分生物量，根据每个树种树龄及其生物量，计算出单位面积林木的年平均总生长量；采集不同树种地上与地下部分样品，测定样品的含碳量；根据不同树种生长量、含碳量，计算单位面积林木的总固碳量，采用瑞典的碳税率计算不同树种的固碳价值；根据林木年生长量，计算不同树种释氧量，然后根据释氧量与当年医用氧气的市场价格，计算林木释氧价值；
(5)土壤保持价值测定与计算：在有林地与无林地标准地中，随机布设小样方并取样，测定样方内的土壤侵蚀模数，林木保土量为无林地土壤侵蚀模数与有林地土壤侵蚀模数的差值，水土流失使氮、磷、钾和有机质大量流失，从而增加土壤的化肥施用量，因此，把森林保持这些养分元素的价值换算成减少化肥施用的费用，计算林木土壤保持的价值；
(6)水调节价值测定与计算：根据水库工程的蓄水成本确定林木涵养水源的价值，水库库容造价采用水库蓄洪工程投资费用，林木净化水质的单价采用当年城市居民用水平均价格；
(7)营养积累价值测定与计算：分别在各标准地采集各树种枝、干与根样品，测定植物样品中氮、磷、钾含量，根据林木生长量与氮、磷、钾含量，计算不同林木积累氮、磷、钾量，把树木积累的氮、磷折合化肥，分别按照当年化肥市场价格，计算出林木营养积累的价值；
(8)综合指标值计算：用Excel 2007与SPSS 19.0软件对不同树种的生态适应性指标值与生态服务价值进行主成分分析，采用隶属函数法对备选树种的以上指标进行综合评价，计算各项指标的隶属度指数，求出各公因子负荷量与权重，算出各树种的综合指标值；
(9)树种选择：根据不同树种的综合指标值排序，综合指标值高的树种适应环境的能力强，生态服务价值高，采用计算机及其数学模型计算机软件，计算不同树种的综合指标值，确定造林绿化选择的树种。
本发明方法新颖独特，科学合理，易操作，可有效用于绿化造林树种的选择，具有很强的实际应用价值，为实现绿化造林提供了技术支撑，有显著的经济和社会效益。
附图说明
图1为本发明的工艺流程图。
图2为本发明树种选择指标体系框示图。
图3为本发明的计算机软件操作流程图。
具体实施方式
以下结合附图和实施例对本发明的具体实施方式作详细说明。
以平原区绿化树种选择为例，说明本发明的具体实施方式。
由图1给出，一种基于主成分计算的绿化树种选择方法，包括以下步骤：
(1)树种选择指标体系构建：由图2所示，根据树种选择需要，结合平原区生态环境条件，选择了生态适应性(1个指标类别，4个指标)与生态服务价值(5个指标类别，10个指标)共6个指标类别14个指标，构建了树种选择指标体系；
(2)测定林地选择：在平原区选择立地条件相近的纯林地，树龄10—15年生，避免由于立地条件的差异而产生的误差；
(3)标准地设置：每个树种设置3个标准地，标准地规格为20×20米，以20×20米的无林地样地作为对照，作为树种生态适应性与生态服务功能测定的场地；
(4)生态适应性测定：在设置的标准地中，测定记录不同树种抗寒、抗旱、抗风、抗病虫害情况；
(5)生态适应性指标值计算：对不同树种抗寒、抗旱、抗风、抗病虫害情况进行分级量化，具体是：
A、不同树种的生态适应性指标值
| 树种 | 抗寒性 | 抗旱性 | 抗风性 | 抗病性 | 总指标值 |
| 悬铃木 | 4 | 4 | 4 | 5 | 17 |
| 石楠 | 3 | 4 | 4 | 4 | 15 |
| 广玉兰 | 3 | 4 | 3 | 4 | 14 |
| 107杨 | 5 | 4 | 3 | 4 | 16 |
| 枇杷 | 3 | 4 | 3 | 4 | 14 |
| 雪松 | 5 | 4 | 4 | 5 | 18 |
| 紫叶李 | 4 | 5 | 4 | 3 | 16 |
| 黄山栾 | 5 | 4 | 4 | 4 | 17 |
| 白蜡 | 5 | 5 | 4 | 5 | 19 |
| 旱柳 | 5 | 4 | 3 | 3 | 15 |
| 垂柳 | 5 | 4 | 3 | 3 | 15 |
| 女贞 | 2 | 4 | 4 | 4 | 14 |
| 侧柏 | 5 | 5 | 5 | 5 | 20 |
| 泡桐 | 4 | 5 | 3 | 3 | 15 |
| 毛白杨 | 4 | 5 | 3 | 4 | 16 |
| 圆柏 | 5 | 4 | 5 | 5 | 19 |
| 国槐 | 5 | 5 | 4 | 5 | 19 |
| 银杏 | 5 | 4 | 4 | 4 | 17 |
B、不同树种的生态服务功能
C、不同树种的生态服务价值(元/(公顷·年)
| 类别 | 气体调节 | 水调节 | 土壤保持 | 营养积累 | 木材价值 |
| 悬铃木 | 62168 | 14928 | 27194 | 3406 | 17644 |
| 石楠 | 36125 | 25380 | 27670 | 1648 | 10226 |
| 广玉兰 | 30509 | 14274 | 27712 | 2056 | 8642 |
| 107杨 | 50639 | 26153 | 35317 | 3915 | 19900 |
| 枇杷 | 35614 | 25237 | 33475 | 1703 | 8662 |
| 雪松 | 32106 | 34393 | 28398 | 1827 | 6264 |
| 紫叶李 | 9146 | 28166 | 29078 | 616 | 1453 |
| 黄山栾 | 27281 | 21522 | 28959 | 1913 | 7743 |
| 白蜡 | 28330 | 29880 | 31865 | 1801 | 25241 |
| 旱柳 | 33621 | 16612 | 28194 | 2566 | 9247 |
| 垂柳 | 29540 | 15789 | 27965 | 2221 | 8125 |
| 女贞 | 29662 | 16862 | 28564 | 1634 | 8397 |
| 侧柏 | 9261 | 33567 | 26939 | 493 | 3918 |
| 泡桐 | 34679 | 16137 | 34122 | 1845 | 17055 |
| 毛白杨 | 29260 | 15821 | 26006 | 1556 | 8305 |
| 圆柏 | 9171 | 34599 | 28342 | 488 | 3879 |
| 国槐 | 40114 | 16116 | 32994 | 2132 | 16430 |
| 银杏 | 48325 | 17498 | 29127 | 2570 | 14626 |
D、不同树种的隶属函数值
分别按照式(1)计算其指标值：
S＝(1X1+2X2+3X3+4X4)/∑X式(1)
式中，S为抗寒、抗旱、抗风、抗病虫害指标值；X1、X2、X3、X4为各树种各等级株数，抗性越强前面的常数越大；∑X为观测总株数；
按照式(2)计算出不同树种的生态适应性总指标值：
S总＝S抗寒+S抗旱+S抗风+S抗病虫害式(2)
式中，S总为各树种生态适应性总指标值，S抗寒为各树种抗寒性指标值，S抗旱为各树种抗旱性指标值，S抗风为各树种抗风指标值，S抗病虫害为各树种抗病虫害指标值；
从表1可以看出，根据生态适应性测定与计算结果，选择的18种主要乔木树种生态适应性总指标值变化范围为14—20：
表1不同树种的生态适应性指标值
(6)气体调节价值测定与计算：在标准地中，用收获法测定林木地上与地下部分生物量，根据每个树种树龄及其生物量，计算出单位面积林木地上与地下部分的年平均生长量。采集不同树种地上与地下部分样品，并测定样品的含碳量，分别计算地上与地下部分的含碳量，合计得到林木的总含碳量。按照式(3)计算林木固碳释氧价值：
U气体调节＝(C碳R碳+1.19C氧)×B年式(3)
式中，U气体调节为林木固碳释氧价值，单位：元/(公顷·年)；C碳为固碳价格，单位：元/吨；C氧为氧气价格，单位：元/吨；R碳为林木总含碳量，单位：‰；B年为林木年总生长量，单位：吨/(公顷·年)；
(7)水调节价值测定与计算：在标准地测定林地年径流量，根据水库工程的蓄水成本确定林木涵养水源的价值，水库库容造价采用水库蓄洪工程投资费用来代替水库工程单位库容造价；林木净化水质费用采用河南省当年城市居民用水平均价格，按照式(4)计算不同生长涵养水源的价值：
U涵养水源＝(10C库+10K)(P－E－C) 式(4)
式中，U涵养水源为水库建设单位库容投资，单位：元/(公顷·年)；K为水的净化费用，单位：元/吨；P为降水量，单位：毫米/年；E为林分蒸散量，单位：毫米/年；C为地表径流量，单位：毫米/年；C库为水库建设单位库容投资，单位：元/立方米；
(8)土壤保持价值测定与计算：在标准地中，测定土壤侵蚀模数，在每个标准地内随机布设3个小样方(0.5×0.5米)，在0—60厘米土层取土3次，测定样方内的土壤有机质、氮、磷、钾含量，按照式(5)计算林木保育土壤的价值：
U固土保肥＝C土(X2－X1)/ρb+(X2－X1)(NC1/R1+PC1/R2+KC2/R3+MC3)式(5)
式中，U固土保肥为固土保肥价值，单位：元/(公顷·年)；X1为林地土壤侵蚀模数，单位：吨/(公顷·年)；X2为无林地土壤侵蚀模数，单位：吨/(公顷·年)；C土为挖取和运输单位体积土方所需费用，单位：元/立方米；ρb为土壤容重，单位：吨/立方米；N为林分土壤平均含氮量，单位：‰；P为林分土壤平均含磷量，单位：‰；K为林分土壤平均含钾量，单位：‰；M为林分土壤平均有机质含量，单位：‰；R1为磷酸二铵化肥含氮量，单位：‰；R2为磷酸二铵化肥含磷量，单位：‰；R3为氯化钾化肥含钾量，单位：‰；C1为磷酸二铵化肥价格，单位：元/吨；C2为氯化钾化肥价格，单位：元/吨；C3为有机质价格，单位：元/吨；
(9)营养积累价值测定与计算：分别在各标准地采集各树种地上与地下部分样品，测定植物样品中氮、磷、钾含量，按照地上与地下部分生长量占总生长量百分比，计算各树种氮、磷、钾含量，按照式(6)计算林木积累营养价值：
U积累营养＝B年(N营养C1/R1+P营养C1/R2+K营养C2/R3)式(6)
式中，U积累营养为林木积累营养价值，单位：元/(公顷·年)；N营养为林木平均含氮量，单位：‰；P营养为林木平均含磷量，单位：‰；K营养为林木平均含钾量，单位：‰；R1为磷酸二铵化肥含氮量，单位：‰；R2为磷酸二铵化肥含磷量，单位：‰；R3为氯化钾化肥含钾量，单位：‰；C1为磷酸二铵化肥价格，单位：元/吨；C2为氯化钾化肥价格，单位：元/吨；B年为林木年总生长量，吨/(公顷·年)；
(10)木材价值测定与计算：测定各树种树干的年生长量与木材密度，计算其单位面积木材材积，按照不同树种的材积及其当年的平均市场价格，扣除木材生产成本(植树与管理费用)，按照式(7)计算各树种生产的价值：
V木＝G年T/ρ(J—C成本)式(7)
式中，V木为木材生产的价值，单位：元/(公顷·年)；G年为林干年生长量，单位：吨/(公顷·年)；T为平均出材率，单位：％；ρ为木材密度，单位：千克/立方米；J为木材销售价格，单位：元/立方米；C成本为平均木材生产成本，单位：元/立方米；
(11)不同树种的生态服务价值：根据测定的不同树种生态服务功能，根据式(3)—(7)计算其生态服务价值；
根据测定结果，在表2中列出了不同树种的生态服务功能，这些基础数据是计算林木生态服务价值的依据；
表2不同树种的生态服务功能
不同树种生态适应性与生态服务价值测定与计算结果列于表3中，由表3可知，不同树种气体调节价值，悬铃木、107杨、银杏、国槐、石楠与枇杷较高，白蜡、黄山栾、圆柏、侧柏与紫叶李较低，其它树种中等；不同树种水调节价值，圆柏、雪松、侧柏、白蜡、紫叶李与107杨较高，泡桐、国槐、毛白杨、垂柳、悬铃木与广玉兰较低，其它树种中等；不同树种土壤保持价值，107杨、泡桐、枇杷、国槐、白蜡与银杏较高，垂柳、广玉兰、石楠、悬铃木、侧柏与毛白杨低，其它树种中等；不同树种营养积累价值，107杨、悬铃木、银杏、旱柳、垂柳与国槐较高，石楠、女贞、毛白杨、紫叶李、侧柏与圆柏较低，其它树种中等；
表3不同树种的生态服务价值(元/(公顷·年)
| 类别 | 气体调节 | 水调节 | 土壤保持 | 营养积累 | 木材价值 |
| 悬铃木 | 62168 | 14928 | 27194 | 3406 | 17644 |
| 石楠 | 36125 | 25380 | 27670 | 1648 | 10226 |
| 广玉兰 | 30509 | 14274 | 27712 | 2056 | 8642 |
| 107杨 | 50639 | 26153 | 35317 | 3915 | 19900 |
| 枇杷 | 35614 | 25237 | 33475 | 1703 | 8662 |
| 雪松 | 32106 | 34393 | 28398 | 1827 | 6264 |
| 紫叶李 | 9146 | 28166 | 29078 | 616 | 1453 |
| 黄山栾 | 27281 | 21522 | 28959 | 1913 | 7743 |
| 白蜡 | 28330 | 29880 | 31865 | 1801 | 25241 |
| 旱柳 | 33621 | 16612 | 28194 | 2566 | 9247 |
| 垂柳 | 29540 | 15789 | 27965 | 2221 | 8125 |
| 女贞 | 29662 | 16862 | 28564 | 1634 | 8397 |
| 侧柏 | 9261 | 33567 | 26939 | 493 | 3918 |
| 泡桐 | 34679 | 16137 | 34122 | 1845 | 17055 |
| 毛白杨 | 29260 | 15821 | 26006 | 1556 | 8305 |
| 圆柏 | 9171 | 34599 | 28342 | 488 | 3879 |
| 国槐 | 40114 | 16116 | 32994 | 2132 | 16430 |
| 银杏 | 48325 | 17498 | 29127 | 2570 | 14626 |
(12)主成分分析：选用1个生态适应性总指标、5个生态服务价值指标类别，共6个指标类别(图2)，筛选乔木树种(表1、2)，由图3所示，用SPSS 19.0软件对不同树种的生态适应性总指标值与生态服务价值进行主成分测定，再用式(8)对树种进行隶属函数值计算；
隶属函数值计算
U(xi)＝(xi-xmin)/(xmax-xmin)i＝1,2,3,…,n 式(8)
式中，U(xi)为隶属度函数值，xi为指标值；xmin和xmax为某一指标的最小值和最大值；
按照公式(9)，根据表2中数据，用Excel 2007软件，计算各因子的隶属函数值(表4)：
表4不同树种的隶属函数值
表5、6给出了不同树种生态适应性指标值与生态服务价值的成分矩阵、解释的总方差，从因子主成分的贡献率来看，第一主成分的贡献率最大，在第一主成分中，气体调节的贡献率均为0.924，营养积累的贡献率为0.922，木材价值的贡献率为0.784，具有较大的荷载；在第二主成份中，生态适应性总指标值的贡献率是0.739，水调节的贡献率为0.618，土壤保持的贡献率均为0.568，具有较大的荷载(表5)：
表5成份矩阵
公因子贡献率计算(Pi)：根据SPSS19.0软件分析结果，得到公因子贡献率，其结果见表6：
表6公因子贡献率
公因子权重计算：根据表6列出了不同因子的主成分分析结果(表7)，由表5、7可知，第一主成分的贡献率不能完全代表原变量的信息，所以计算中用了第一、第二主成分的值。用公式(9)，求出第一、第二主成分的因子权重，然后用Excel软件求出公因子权重(表8)：
式中，Wi表示第i个公因子在所有公因子中的权重，i为各树种第i个公因子，Pi表示了各树种第i个公因子的贡献率；
表7不同树种生态适应性指标值与生态服务价值因子贡献率
| 主成分 | 1 | 2 | 3 | 4 | 5 | 6 |
| 贡献率 | 51.613 | 24.946 | 12.950 | 6.213 | 2.798 | 1.480 |
| 累计贡献率 | 51.613 | 76.559 | 89.509 | 95.723 | 98.520 | 100 |
表8因子负荷与权重
(13)综合指标值计算：根据隶属度函数值与公因子权重，根据式(10)，用Excel软件计算各树种的综合指标值(表10)：
式中，D值为各树种生态适应性与生态服务价值综合指标值；U(xi)为隶属度函数值；Wi表示第i个公因子在所有公因子中的权重；
单一的因子评价的重点集中于某些主要因子的变化上，不能反映树木生态适应性指标值与生态服务价值的总体情况，所以树木生态适应性指标值与生态服务价值总体变化难以精确的比较，综合指标值能够全面反映各因子的总体情况，其值越大，生态适应性与生态服务价值指标值越高；
(14)树种选择：由表10可知，不同乔木树种综合指标值：107杨＞悬铃木＞国槐＞白蜡＞泡桐＞银杏＞枇杷＞旱柳＞黄山栾＞垂柳＞石楠＞广玉兰与雪松＞女贞＞毛白杨＞圆柏＞侧柏＞紫叶李，综合指标值反映了不同树种生态适应性与生态服务价值的总体差异性，综合指标值可作为树种选择的依据，在造林绿化工作中优先选择综合指标值较高，即综合表现较好的树种。
表10不同树种的综合指标值
上述量化指标的计算，均采用计算机及其计算机数学模型软件，用Excel 2007软件计算公因子权重，根据隶属度指标值与公因子权重，用Excel 2007软件计算不同树种的综合指标值，确定树种选择，具体操作方法见图3所示。
由上述可以看出，本发明是以平原区主要绿化树种为研究对象，对不同树种的生态适应性与生态服务功能进行测定，以这些树种生态适应性总指标值，气体调节、土壤保持、水调节、营养积累、生产木材的价值为依据，采用主成份分析法，对各树种以上指标进行主成分分析，计算各树种综合指标值，对造林树种进行量化选择，以保证绿化质量，根据树种生态的自然规律，采用科学的技术手段，以数学模型，采用计算机及其软件，最终实现绿化树种的选择，方法易操作，符合当地选择树种绿化造林的实际，并经反复多次实验，均取得了相同和相近似的结果，并与当地实地情况验证，非常一致，表明方法稳定可靠，具有实际的应用价值，可为绿化造林提供有效的技术手段，实现科学绿化造林的种植和管理，有显著的经济和社会效益。
Claims (2)
1.一种基于主成分计算的绿化树种选择方法，其特征在于，对不同树种的生态适应性与生态功能进行测定，以这些树种生态适应性总指标值，气体调节、土壤保持、水调节、营养积累、生产木材的价值为依据，采用主成分分析法，计算树种综合指标值，对造林树种进行量化选择，具体包括以下步骤：
(1)建立树种选择指标体系：根据树种的生态适应性，对树种进行筛选，在树种适应环境条件下，才能够提供较高的生态价值，在构建树种选择指标体系时，选择了生态适应性与生态价值指标；
(2)设置标准地：在土壤气候条件相近的苗圃设为标准地，作为树种生态适应性与生态功能测定的场地；
(3)生态适应性指标值测定与计算：在标准地中，观测记录不同树种抗寒、抗旱、抗病虫害、抗风情况，对以上4个指标进行测定分级量化，然后合计得出不同树种生态适应性总指标值；
(4)气体调节价值测定与计算：在标准地中，用收获法测算林木地上与地下部分生物量，根据每个树种树龄及其生物量，计算出单位面积林木的年平均总生长量；采集不同树种地上与地下部分样品，测定样品的含碳量；根据不同树种生长量、含碳量，计算单位面积林木的总固碳量，采用瑞典的碳税率计算不同树种的固碳价值；根据林木年生长量，计算不同树种释氧量，然后根据释氧量与当年医用氧气的市场价格，计算林木释氧价值；
(5)土壤保持价值测定与计算：在有林地与无林地标准地中，随机布设小样方并取样，测定样方内的土壤侵蚀模数，林木保土量为无林地土壤侵蚀模数与有林地土壤侵蚀模数的差值，水土流失使氮、磷、钾和有机质大量流失，从而增加土壤的化肥施用量，因此，把森林保持这些养分元素的价值换算成减少化肥施用的费用，计算林木土壤保持的价值；
(6)水调节价值测定与计算：根据水库工程的蓄水成本确定林木涵养水源的价值，水库库容造价采用水库蓄洪工程投资费用，林木净化水质的单价采用当年城市居民用水平均价格；
(7)营养积累价值测定与计算：分别在各标准地采集各树种枝、干与根样品，测定植物样品中氮、磷、钾含量，根据林木生长量与氮、磷、钾含量，计算不同林木积累氮、磷、钾量，把树木积累的氮、磷折合化肥，分别按照当年化肥市场价格，计算出林木营养积累的价值；
(8)综合指标值计算：用Excel 2007与SPSS 19.0软件对不同树种的生态适应性指标值与生态服务价值进行主成分分析，采用隶属函数法对备选树种的以上指标进行综合评价，计算各项指标的隶属度指数，求出各公因子负荷量与权重，算出各树种的综合指标值；
(9)树种选择：根据不同树种的综合指标值排序，综合指标值高的树种适应环境的能力强，生态服务价值高，采用计算机及其数学模型计算机软件，计算不同树种的综合指标值，确定造林绿化选择的树种。
2.根据权利要求1所述的基于主成分计算的绿化树种选择方法，其特征在于，包括以下步骤：
(1)树种选择指标体系构建：根据树种选择需要，结合平原区生态环境条件，选择生态适应性指标抗寒性、抗旱性、抗风性、抗病性和生态服务价值的气体调节、水调节、土壤保持、营养积累、木材价值，构建树种选择指标体系；
(2)测定林地选择：在平原区选择立地条件相近的纯林地，树龄10—15年生，避免由于立地条件的差异而产生的误差；
(3)标准地设置：每个树种设置3个标准地，标准地规格为20×20米，以20×20米的无林地样地作为对照，作为树种生态适应性与生态服务功能测定的场地；
(4)生态适应性测定：在设置的标准地中，测定记录不同树种抗寒、抗旱、抗风、抗病虫害情况；
(5)生态适应性指标值计算：对不同树种抗寒、抗旱、抗风、抗病虫害情况进行分级量化，具体是：
A、不同树种的生态适应性指标值
B、不同树种的生态服务功能
C、不同树种的生态服务价值(元/(公顷·年)
D、不同树种的隶属函数值
分别按照式(1)计算其指标值：
S＝(1X1+2X2+3X3+4X4)/∑X 式(1)
式中，S为抗寒、抗旱、抗风、抗病虫害指标值；X1、X2、X3、X4为各树种各等级株数，抗性越强前面的常数越大；∑X为观测总株数；
按照式(2)计算出不同树种的生态适应性总指标值：
S总＝S抗寒+S抗旱+S抗风+S抗病虫害 式(2)
式中，S总为各树种生态适应性总指标值，S抗寒为各树种抗寒性指标值，S抗旱为各树种抗旱性指标值，S抗风为各树种抗风指标值，S抗病虫害为各树种抗病虫害指标值；
根据生态适应性测定与计算结果，选择的18种主要乔木树种生态适应性总指标值变化范围为14—20：
(6)气体调节价值测定与计算：在标准地中，用收获法测定林木地上与地下部分生物量，根据每个树种树龄及其生物量，计算出单位面积林木地上与地下部分的年平均生长量；采集不同树种地上与地下部分样品，并测定样品的含碳量，分别计算地上与地下部分的含碳量，合计得到林木的总含碳量，按照式(3)计算林木固碳释氧价值：
U气体调节＝(C碳R碳+1.19C氧)×B年 式(3)
式中，U气体调节为林木固碳释氧价值，单位：元/(公顷·年)；C碳为固碳价格，单位：元/吨；C氧为氧气价格，单位：元/吨；R碳为林木总含碳量，单位：‰；B年为林木年总生长量，单位：吨/(公顷·年)；
(7)水调节价值测定与计算：在标准地测定林地年径流量，根据水库工程的蓄水成本确定林木涵养水源的价值，水库库容造价采用水库蓄洪工程投资费用来代替水库工程单位库容造价；林木净化水质费用采用当年城市居民用水平均价格，按照式(4)计算不同生长涵养水源的价值：
U涵养水源＝(10C库+10K)(P－E－C) 式(4)
式中，U涵养水源为水库建设单位库容投资，单位：元/(公顷·年)；K为水的净化费用，单位：元/吨；P为降水量，单位：毫米/年；E为林分蒸散量，单位：毫米/年；C为地表径流量，单位：毫米/年；C库为水库建设单位库容投资，单位：元/立方米；
(8)土壤保持价值测定与计算：在标准地中，测定土壤侵蚀模数，在每个标准地内随机布设3个0.5×0.5米的样方，在0—60厘米土层取土3次，测定样方内的土壤有机质、氮、磷、钾含量，按照式(5)计算林木保育土壤的价值：
U固土保肥＝C土(X2－X1)/ρb+(X2－X1)(NC1/R1+PC1/R2+KC2/R3+MC3) 式(5)
式中，U固土保肥为固土保肥价值，单位：元/(公顷·年)；X1为林地土壤侵蚀模数，单位：吨/(公顷·年)；X2为无林地土壤侵蚀模数，单位：吨/(公顷·年)；C土为挖取和运输单位体积土方所需费用，单位：元/立方米；ρb为土壤容重，单位：吨/立方米；N为林分土壤平均含氮量，单位：‰；P为林分土壤平均含磷量，单位：‰；K为林分土壤平均含钾量，单位：‰；M为林分土壤平均有机质含量，单位：‰；R1为磷酸二铵化肥含氮量，单位：‰；R2为磷酸二铵化肥含磷量，单位：‰；R3为氯化钾化肥含钾量，单位：‰；C1为磷酸二铵化肥价格，单位：元/吨；C2为氯化钾化肥价格，单位：元/吨；C3为有机质价格，单位：元/吨；
(9)营养积累价值测定与计算：分别在各标准地采集各树种地上与地下部分样品，测定植物样品中氮、磷、钾含量，按照地上与地下部分生长量占总生长量百分比，计算各树种氮、磷、钾含量，按照式(6)计算林木积累营养价值：
U积累营养＝B年(N营养C1/R1+P营养C1/R2+K营养C2/R3) 式(6)
式中，U积累营养为林木积累营养价值，单位：元/(公顷·年)；N营养为林木平均含氮量，单位：‰；P营养为林木平均含磷量，单位：‰；K营养为林木平均含钾量，单位：‰；R1为磷酸二铵化肥含氮量，单位：‰；R2为磷酸二铵化肥含磷量，单位：‰；R3为氯化钾化肥含钾量，单位：‰；C1为磷酸二铵化肥价格，单位：元/吨；C2为氯化钾化肥价格，单位：元/吨；B年为林木年总生长量，吨/(公顷·年)；
(10)木材价值测定与计算：测定各树种树干的年生长量与木材密度，计算其单位面积木材材积，按照不同树种的材积及其当年的平均市场价格，扣除木材生产成本，按照式(7)计算各树种生产的价值：
V木＝G年T/ρ(J—C成本) 式(7)
式中，V木为木材生产的价值，单位：元/(公顷·年)；G年为林干年生长量，单位：吨/(公顷·年)；T为平均出材率，单位：％；ρ为木材密度，单位：千克/立方米；J为木材销售价格，单位：元/立方米；C成本为平均木材生产成本，单位：元/立方米；
(11)不同树种的生态服务价值：根据测定的不同树种生态服务功能，根据式(3)—(7)计算其生态服务价值；
根据测定结果，列出不同树种的生态服务数据，作为计算林木生态服务价值的依据；
不同树种的生态服务数据如下：
悬铃木、107杨、银杏、国槐、石楠与枇杷高，白蜡、黄山栾、圆柏、侧柏与紫叶李低；不同树种水调节价值，圆柏、雪松、侧柏、白蜡、紫叶李与107杨高，泡桐、国槐、毛白杨、垂柳、悬铃木与广玉兰低；不同树种土壤保持价值，107杨、泡桐、枇杷、国槐、白蜡与银杏高，垂柳、广玉兰、石楠、悬铃木、侧柏与毛白杨低；不同树种营养积累价值，107杨、悬铃木、银杏、旱柳、垂柳与国槐高，石楠、女贞、毛白杨、紫叶李、侧柏与圆柏低；
(12)主成分分析：选用1个生态适应性总指标、5个生态服务价值指标类别，共6个指标类别，分别是生态适应性、气体调节、水调节、土壤保持、营养积累、木材生产，筛选乔木树种，用SPSS 19.0软件对不同树种的生态适应性总指标值与生态服务价值进行主成分测定，再用式(8)对树种进行隶属函数值计算；
U(xi)＝(xi-xmin)/(xmax-xmin)i＝1,2,3,…,n 式(8)
式中，U(xi)为隶属度函数值，xi为指标值；xmin和xmax为某一指标的最小值和最大值；
按照式(9)，用Excel 2007软件，计算各因子的隶属函数值；
气体调节的贡献率均为0.924，营养积累的贡献率为0.922，木材价值的贡献率为0.784，生态适应性总指标值的贡献率是0.739，水调节的贡献率为0.618，土壤保持的贡献率均为0.568；
公因子贡献率计算(Pi)：根据SPSS 19.0软件计算，得到公因子贡献率，分别是：
公因子权重计算：根据不同因子成分计算结果，由表5、7可知，第一主成分的贡献率不能完全代表原变量的信息，所以计算中用了第一、第二主成分的值，用式(9)，求出第一、第二主成分的因子权重，然后用Excel软件求出公因子权重：
式中，Wi表示第i个公因子在所有公因子中的权重，i为各树种第i个公因子，Pi表示了各树种第i个公因子的贡献率；
(13)综合指标值计算：根据隶属度函数值与公因子权重，根据式(10)，用Excel软件计算各树种的综合指标值：
式中，D值为各树种生态适应性与生态服务价值综合指标值；U(xi)为隶属度函数值；Wi表示第i个公因子在所有公因子中的权重；
单一的因子评价的重点集中于某些主要因子的变化上，不能反映树木生态适应性指标值与生态服务价值的总体情况，所以树木生态适应性指标值与生态服务价值总体变化难以精确的比较，综合指标值能够全面反映各因子的总体情况，其值越大，生态适应性与生态服务价值指标值越高；
(14)树种选择：根据数学模型、计算机软件，由计算机计算出不同乔木树种综合指标值：107杨＞悬铃木＞国槐＞白蜡＞泡桐＞银杏＞枇杷＞旱柳＞黄山栾＞垂柳＞石楠＞广玉兰与雪松＞女贞＞毛白杨＞圆柏＞侧柏＞紫叶李，综合指标值反映了不同树种生态适应性与生态服务价值的总体差异性，以综合指标值作为树种选择的依据，选出造林绿化所需的树种。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910208015.XA CN110009009B (zh) | 2019-03-19 | 2019-03-19 | 一种基于主成分计算的绿化树种选择方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910208015.XA CN110009009B (zh) | 2019-03-19 | 2019-03-19 | 一种基于主成分计算的绿化树种选择方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110009009A CN110009009A (zh) | 2019-07-12 |
| CN110009009B true CN110009009B (zh) | 2022-11-22 |
Family
ID=67167379
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910208015.XA Active CN110009009B (zh) | 2019-03-19 | 2019-03-19 | 一种基于主成分计算的绿化树种选择方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110009009B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3865764B1 (ja) * | 2006-08-15 | 2007-01-10 | アルスマエヤ株式会社 | 森林資源調査方法および森林資源調査装置 |
| RU2489846C2 (ru) * | 2011-09-16 | 2013-08-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Брянская государственная инженерно-технологическая академия" | Способ комплексной оценки состояния лесных экосистем в районах техногенного воздействия промышленных объектов |
| CN104778626A (zh) * | 2015-03-25 | 2015-07-15 | 山西大学 | 一种造林树种生态适应性等级评价方法 |
| CN106547948A (zh) * | 2016-10-10 | 2017-03-29 | 天津农学院 | 用于引种的园林绿化小乔木生长适应性的综合评价方法 |
-
2019
- 2019-03-19 CN CN201910208015.XA patent/CN110009009B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3865764B1 (ja) * | 2006-08-15 | 2007-01-10 | アルスマエヤ株式会社 | 森林資源調査方法および森林資源調査装置 |
| RU2489846C2 (ru) * | 2011-09-16 | 2013-08-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Брянская государственная инженерно-технологическая академия" | Способ комплексной оценки состояния лесных экосистем в районах техногенного воздействия промышленных объектов |
| CN104778626A (zh) * | 2015-03-25 | 2015-07-15 | 山西大学 | 一种造林树种生态适应性等级评价方法 |
| CN106547948A (zh) * | 2016-10-10 | 2017-03-29 | 天津农学院 | 用于引种的园林绿化小乔木生长适应性的综合评价方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110009009A (zh) | 2019-07-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Tian et al. | Climate and land use controls over terrestrial water use efficiency in monsoon Asia | |
| Morales et al. | Comparing and evaluating process‐based ecosystem model predictions of carbon and water fluxes in major European forest biomes | |
| Lindroth et al. | Environmental controls on the CO2 exchange in north European mires | |
| Herbert et al. | Forest productivity and efficiency of resource use across a chronosequence of tropical montane soils | |
| Targhetta et al. | Effects of the hydroedaphic gradient on tree species composition and aboveground wood biomass of oligotrophic forest ecosystems in the central Amazon basin | |
| Malmer et al. | Accumulation and release of organic matter in ombrotrophic bog hummocks‐processes and regional variation | |
| CN109816207B (zh) | 一种基于大数据的灌区节水潜力测定方法及装置 | |
| Vacek et al. | Structural diversity and production of alder stands on former agricultural land at high altitudes | |
| Di Vittorio et al. | Development and optimization of an Agro-BGC ecosystem model for C4 perennial grasses | |
| CN111340644A (zh) | 基于碳储量分布状态分析小兴安岭针阔混交林碳汇结构的方法 | |
| CN111340643A (zh) | 基于碳储量分析小兴安岭针阔混交林类型及群落学特征的方法 | |
| CN113706012A (zh) | 一种耕地资源资产生态价值评估方法及模型 | |
| Han et al. | Modeling grassland net primary productivity and water-use efficiency along an elevational gradient of the Northern Tianshan Mountains | |
| Trouillier et al. | A unifying concept for growth trends of trees and forests–the “potential natural forest” | |
| CN114088916A (zh) | 一种用于碳汇农业的增加土壤碳储备的方法 | |
| Farooq et al. | Spatial distribution of carbon dynamics and nutrient enrichment capacity in different layers and tree tissues of Castanopsis eyeri natural forest ecosystem | |
| Kennedy et al. | Amounts, forms, and management of nitrogen and phosphorus export from agricultural peatlands | |
| Wang et al. | Dynamics of carbon and nitrogen storage in two typical plantation ecosystems of different stand ages on the Loess Plateau of China | |
| Yang et al. | Alpine wetland degradation reduces carbon sequestration in the Zoige Plateau, China | |
| Huang et al. | Soil organic carbon content and storage of raised field wetlands in different functional zones of a typical shallow freshwater lake, China | |
| Li et al. | Seasonal variation characteristics of c, n, and p stoichiometry and water use efficiency of Messerschmidia sibirica and its relationship with soil nutrients | |
| CN110009009B (zh) | 一种基于主成分计算的绿化树种选择方法 | |
| Schlutow et al. | Rule-based classification and mapping of ecosystem services with data on the integrity of forest ecosystems | |
| McLaughlin et al. | Growth Responses of 53 Open‐Pollinated Loblolly Pine Families to Ozone and Acid Rain | |
| Lélé | Sustainable use of biomass resources: A note on definitions, criteria, and practical applications |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |