基于主成分和聚类分析的永泰县土壤养分地球化学综合评价_江晓龙.pdf

基于主成分和聚类分析的永泰县土壤养分地球化学综合评价江晓龙(福建省煤田地质勘查院，福建福州 350000)摘要 目的 运用统计学方法对永泰县土壤养分地球化学综合等级进行评价，最大限度地减少人为主观因素的影响。方法 对工作区内土壤进行系统采样，送检测试土壤养分大量元素(有机质、碱解氮、有效磷、速效钾、全氮、全磷、全钾)含量及 pH，通过采样数据的一般性描述，初步判断工作区内土壤养分元素富集分等情况，随后利用 SPSS 软件对数据进行标准化处理，建立相关系数矩阵，判断其进行统计学分析的可行性。通过主成分分析计算筛选出评价因子主成分，计算综合评价得分，并利用 K均值聚类分析法将综合得分划分为缺乏、较缺乏、中等、较富集、富集 5 类。结果 工作区内土壤养分元素分布不均匀，变异性高，土壤总体呈酸性，富氮缺磷钾。土壤地球化学养分综合状况较差，综合得分属较缺乏、缺乏等级占比达 5468%。结论 该研究为土壤养分评价提供科学依据。关键词主成分分析;聚类分析;土壤养分;地球化学;综合评价;永泰县中图分类号S1519文献标识码A文章编号05176611(2023)01006804doi:103969/jissn05176611202301015开放科学(资源服务)标识码(OSID):Comprehensive Evaluation of Soil Nutrients Geochemistry in Yongtai County Based on Principal Component Analysis and ClusterAnalysisJIANG Xiao-long(Fujian Coalfield Geological Prospecting Institute，Fuzhou，Fujian 350000)Abstract ObjectiveTo evaluate the comprehensive geochemical grade of soil nutrients in Yongtai County by using statistical methods，soas to minimize the impact of human subjective factors MethodSystematically sample the soil in the working area，the content of soil nutrientmacroelements(organic matter，alkali hydrolyzable nitrogen，available phosphorus，available potassium，total nitrogen，total phosphorus，to-tal potassium)and pH were submitted for testing，through the general description of the sampling data，the enrichment and grading of soil nu-trient elements in the work area were preliminarily judgedThen，SPSS software was used to standardize the data and establish a correlation co-efficient matrix to judge the feasibility of statistical analysisThe principal components of the evaluation factors were selected through principalcomponent analysis，and the comprehensive evaluation scores were calculated The comprehensive scores were divided into five categories:lack，relatively lack，medium，relatively enriched and enriched by using K-means cluster analysis esult The soil nutrients elements in theworking area was uneven distribution，high variability，and the soil was generally acidic，nitrogen enrichment，lack of potassium and phosphor-us The comprehensive nutrient status was poor，the lack and relatively lack of comprehensive level accounted for 5468%ConclusionThisstudy provides scientific basis for soil nutrient evaluationKey wordsPrincipal component analysis;Cluster analysis;Soil nutrients;Geochemistry;Comprehensive evaluation;Yongtai County基金项目福建省财政项目(2018001)。作者简介江晓龙(1992)，男，福建龙岩人，工程师，硕士，从事农业地质调查工作。收稿日期20220227土壤是物质生产最基本的载体，是作物吸收营养的主要渠道。农业生产中土壤的贡献率在 50%60%，土壤养分的丰缺是影响作物产量、品质的基本因素12，是地力评价的主要指标，因此，查明土壤养分状况对农业生产的合理化调配显得尤其重要。通过对工作区内的土壤养分地球化学情况进行综合评价，可为土壤资源开发管理提供一定的数据参考。土壤养分综合评价的方法主要有聚类分析法3、模糊综合评判法45、BP 神经网络法6、主成分分析法78、层次分析法9 等。近年来，主成分分析法因其可通过降维方式找出具有代表性、综合性的影响因子而避免指标重叠的特性，得到学者的广泛应用，成为土壤定量研究中应用最为广泛的数理统计方法1011。该研究以永泰县土壤为研究目标，共采集土壤样品 1 955 件，测试氮、磷、钾、有机质等大量养分指标，通过统计学方法，运用主成分分析、聚类分析筛选出影响土壤肥力的主成分及其影响权重大小，对土壤养分地球化学等级进行较为客观的综合评价，尽可能地减少人为主观因素干扰。1材料与方法11研究区概况永泰县隶属福州市管辖，地理位置为11807081191230E、253900260421N，土地面积共2 22986 km2。全县常年平均气温 146260，年无霜期近300 d，日照时数 1 4452 193 h，年降水量 1 4002 000 mm。凭借其得天独厚的自然环境，孕育出芙蓉李、绿茶、白云槟榔芋等一系列名特优农产品，获得国家地理标志证明商标、农产品地理标志产品等。成土母质以沉积岩、中酸性侵入岩、火山岩为主，在多种因素的综合作用下，形成了诸多土壤类型。根据第二次土壤普查工作，查明县域地表土壤有 8 个土类，17 个亚类，包括红壤、黄壤、黄棕壤、草甸土、紫色土、石灰(岩)土、潮土、水稻土。12评价方法121土壤样品采集。以 福建省农业地质调查评价工作手册12 为依据，对永泰县域内耕地现状图斑进行评价单元划分，在评价单元中面积较大的地块内通过 GPS 定位采集土壤样品，根据评价单元面积大小采集 35 个子样品，再以等量混合的形式组合成 1 件样品，共采集耕层土壤 1 955 件。经室内风干、过筛后按四分法缩分留样备用。122样品测试与方法。样品的检测根据 土地质量地球化学评价规范(DZ/T 02952016)、福建省农业地质调查评价工作手册 等相关规范规程由福建省 121 地质大队化验测试中心完成。土壤氮测定采用凯氏定氮法;土壤磷、钾测定采用 X荧光光谱法(XF);有机质测定采用化学容量法安徽农业科学，JAnhui AgricSci 2023，51(1):6871(VOL)，碱解氮测定采用碱解扩散法;速效磷测定采用氟化铵盐酸提取(碳酸氢钠溶液浸提)钼锑抗比色法;速效钾测定采用乙酸铵提取等离子体发射光谱法。13数据处理利用 SPSS 统计软件对样品测试结果进行统计学分析，在标准化处理后构建相关系数矩阵并进行KMO、Bartlett 球形度检验，以判断进行因子分析的可行性。累计贡献率大于 80%作为选择主成分的前提条件，随后建立主成分的因子荷载矩阵，通过换算得出主成分得分系数，以计算主成分分值，进一步通过综合得分公式求取各评价单元土壤养分综合得分值(F)，随即根据综合得分值对土壤样品进行 K均值聚类统计分析，评价工作区内土壤养分地球化学综合等级。主成分得分计算公式7:Fni=zj=1(ijZnj)(1)式中，Fni代表第 n 个采样点第 i 项主成分的得分;ij代表第 i项主成分第 j 项原始指标的载荷;Znj代表第 n 个采样点第 j项原始指标的标准化数据。综合得分计算公式13:F=3i=1(iFni)(2)式中，F 表示综合得分，i表示第 i 项主成分贡献率。2结果与分析21土壤养分指标一般性描述根据参数统计结果(表 1)可见，工作区内土壤养分含量分布不均匀，从偏度、峰度可看出，养分含量的分布形态多呈高尖峰右偏态分布，即养分含量总体处于相对低值区，但存在少量的异常值，其中有效磷变化幅度最大，极差达 1 01847 mg/kg。从样品变异系数可看出，有机质、碱解氮、全氮、全钾变异系数为 34%39%，属于中等变异;有效磷、速效钾、全磷变异系数不小于 69%，呈高变异性，异常高值的存在一方面可能是受人为活动影响，另一方面可能与成土母质不同有关。pH 中位数为 502，众数为 479，均值 504，变异系数仅为 8%，属于弱变异，表明工作区内土壤酸碱度差异性较小，区内土壤总体呈酸性。表 1土壤养分指标描述统计Table 1The descriptive statistics of nutrients content index项目Item有机质OMg/kg碱解氮Alkelin-Nmg/kg有效磷Olsen-Pmg/kg速效钾AKmg/kg全氮TNg/kg全磷TPg/kg全钾TKg/kgpH极小值 Minimum246465005737009006110345极大值 Maximum9931331981 01852511603304132609754均值 Mean244611161537779321210501211504众数 Mode1662102910052757108006110479变 异 系 数 CV%3936166693775348偏度 Skewness141088404242094214022045峰度 Kurtosis457174248294312283503416622土壤养分指标等级根据全国第二次土壤普查养分分级标准，从表 2 可以看出，有机质含量 3 等(中等)以上评价单元比例达 6450%;约 6829%的评价单元内碱解氮含量水平在 3 等(中等)以上;全氮含量 3 等(中等)以上比例为63.68%;有 效 磷 含 量 在 33 50%的 评 价 单 元 中 超 过 了40 mg/kg，表现为 1 等(丰富);全磷含量水平较低，7208%的评价单元全磷等级在 4 等(较缺乏)以下;速效钾、全钾等级较低，以 4 等(较缺乏)为主。表 2土壤养分含量等级分布比例Table 2Distribution proportion of soil nutrient content grades等级Grade有机质 OM含量Contentg/kg比例atio%