基于因子分析和熵权法的赣江源流域水资源承载力研究

基于因子分析和熵权法的赣江源流域水资源承载力研究
刘慧，蔡定建，许宝泉，余长林
（江西理工大学材料与化学工程学院，江西赣州341000）
摘要以赣江源流域地区为例，运用因子分析法，从经济社会、生态环境、自然水资源和人口现状等16个指标中筛选出影响该区水资源
承载力变化的4个主因子C 1、C 2、C 3、C 4，并通过熵权法对4个因子的得分赋权重，综合评价出2001 2009年该源流域水资源承载力的状
况。

结果表明，近10年来，该区域水资源承载力总体走向一个相对健康的方向。

通过对赣江源流域水资源利用现状的分析，从水量分配、水权归属及产业结构优化等方面提出提高水资源承载能力的措施，以期为该区域经济社会发展规划和生态环境保护战略的制定，寻
求可持续发展模式提供水资源承载信息的参考和保障。

关键词水资源承载力；因子分析法；熵权法；赣江源流域中图分类号F323．213文献标识码A 文章编号0517－6611（2011）23－14264－04Study on Water Resource Carrying Capacity of the Headstream of Ganjiang River ’s Drainage area Based on Component Analysis and Entropy Weight Method
LIU Hui et al （School of Material and Chemistry Engineering ，Jiangxi University of Science and Technology ，Ganzhou ，Jiangxi 341000）Abstract Taking the headstream of Ganjiang River ’s drainage area as an example ，and by means of component analysis ，4extracted factors out of the 16factors ，which selected from the field of society ，economy ，population and natural water resource that assumed affecting the water resource carrying capacity of this area ，were identified as major components which include C 1，C 2，C 3and C 4．Then using the entropy value method ，weights were respectively bestowed to the scores of the 4components ，to give a comprehensive analysis representing the water re-source carrying capacity of the Ganjiang River ’s headstream from the year of 2001to 2009．As indicated by the annual scores ，the water re-source carrying capacity in this area had been developing towards a healthy direction in recent years．In order to improve the water resource carrying capacity and give reference to policy making on sustainable development in the future ，measures should be taken to govern the water right ，distribute the water in a reasonable way and adjust the industry structure．
Key words Water Resource Carrying Capacity ；Component Analysis ；Entropy Weight Method ；The Headstream of Ganjiang River ’s Drain-age Area
基金项目
国家自然科学基金项目（21067004）；江西省科技厅自然科
学基金项目（2007GZ0698）；江西省高校省级教改项目（JXJG-07-6-20）。

作者简介蔡定建（1959－），男，江西南昌人，副教授，硕士，硕士生导
师，从事循环经济与自然资源教学与科研等工作，
E-mail ：caidingjian@126．com 。

收稿日期2011-
05-18水资源承载力是可持续发展过程中各种自然资源承载能力的重要组成部分，是制定经济社会发展与生态环境保护
战略，寻求可持续发展模式的重要依据［1］。

从20世纪80年
代我国学者引入水资源承载力概念以来，
通过供需分析法、模糊综合评判法、主成分分析法和熵值法等一系列方法，在
流域、区域和城市研究等方面取的了丰富的研究成果
［2－5］。

笔者采用因子分析和熵权法
［6－10］
对赣江源流域的水资源承
载力进行研究，利用因子分析法的降维思想，充分揭示经济、社会、人口和生态环境等变量因子之间的规律性，并通过熵值法进行客观赋权，分析出近9年来水资源承载力的综合状况，
以期为该区经济社会的可持续发展和人民生活提供水资源优化配置方面的支撑和保障。

1
赣江源流域赣州地区的自然概况
古城赣州被誉为“千里赣江第一城”，山区为赣江正源贡水的发源地。

上犹江和章水汇成章江；梅江、琴江、绵江、湘水、濂水、平江和桃江等汇成贡江。

章贡两江汇集流经赣州18个县市的数百条河流，于章贡区八境台下交汇，始称赣江。

而后曲折北上，经由万安等地到南昌，分4支注入鄱阳湖，属长江流域的一级水系。

赣江全长为751km ，以赣州和新干为区划点，分为上游、中游和下游3段，总流域面积达8．28万km 2，约占江西省总面积的50%左右，在古城赣州境内的流域
面积达3．64万km 2
，约占总流域面积的44%左右，赣州全区
总体处在赣江源流域的位置
［11］。

赣州地处中亚热带南缘，属亚热带湿润季风气候区，具有春夏降水集中、气候温和和雨量充沛等气候特征。

该区年
均降水量为1601mm ，
年均水资源总量约335．7亿m 3
，人均占有量约为3900m 3，比全国人均2300m 3
高出近70%，水
资源量丰富，
属富水区［12］。

但是，水量在年内及年际间变化
较大，呈现出明显的季节性和区域性，分配上出现极不均匀的情况，降水主要集中在4 6月，占年降水量的46．8%左右；7 9月，降水少，气温高，工农业用水正值高峰，江河却处在少水期
［13］。

近年来，随着该区经济、人口的不断增长，取
用水消耗量，废污水排放量也急剧增加，水质型和管理型缺水的问题逐渐凸显出来。

根据2009年
《江西省水资源公报》显示，绵江、贡水、章江、平江和桃江等赣江源头处的支流均出现不同程度的污染，且非汛期比汛期的污染更加严重。

与此同时，沿岸水土流失的加剧，致使赣江水中泥沙含量急增，出现河道淤积，河床抬高等威胁沿岸居民生产和生活的问题，亟待规划出长远的水资源承载方案。

2研究方法2．1
因子分析法
因子分析的目的在于揭示观测变量之
间的内在关联性，在尽可能保留原有信息量的前提下，用较少的维度去表示原来的数据结构，简化数据（降维），以形成少数几组潜在起支配作用的公共因子与特殊因子的线性组合，从而揭示整体变量的规律或本质，降低在多元分析处理过程中由于变量多，关系混乱而导致局面错综复杂的问题
［14－15］。

因子分析的计算步骤如下：①原始数据标准化，消
除变量间在数量级和量纲上的不同；②求标准化数据的相关矩阵；③求相关矩阵的特征值和特征向量；④计算方差贡献率与累积方差贡献率；⑤确定因子：前m 个因子包含的数据信息总量（即其累积贡献率）不低于85%时，可取前m 个公因子来反映原有的评价指标；⑥进行适当的因子旋转，
为公责任编辑姜丽责任校对卢瑶
安徽农业科学，Journal of Anhui Agri．Sci．2011，39（23）：14264－14267，14277
因子寻求最佳的实际解释方式；⑦采用回归分析方法求各因子得分。

该研究采用因子分析法产生的主成分得分新变量，再通过熵值法对主因子得分进行客观赋权，计算出赣江源流域水资源承载力的综合得分。

2．2熵权法熵是信息论中对不确定性的一种度量，信息量越大，不确定性就越小，熵也就越小；反之亦然［9］。

利用熵权法计算各指标的权重系数，根据评价指标的离散程度客观的计算权重，因子的离散程度越大，对水资源承载力变化的影响就越大。

熵权法的计算步骤如下。

（1）数据指标同趋势处理：非负化处理各成分因子的得分数据，为避免求熵值时对数的无意义，需要将数据进行平移。

处理如下：
x ij =［X
ij
－min（X
ij
）］/［max（X
ij
）－min（X
ij
）］+1（i=1，2，…，n；j=1，2，…，m）（1）
式中，X
ij
为某一个地区i时段水资源承载力第j个影响因子的值。

（2）计算第i时段第j项因子占所有时段因子之和的比例：
y ij =x
ij
/（∑n
i=1
x
ij
）（i=1，2，…，n；j=1，2，…，m）（2）
（3）求各因子的权重：
w
j =d
j
/（∑m
j=1
d
j
）（1≤j≤m）（3）
式中，d
j =（1－e
j
）/（m－∑m
j=1
e
j
）（0≤d
j
≤1，∑m
j=1
d
j
=1），为第j项
指标的信息效用值。

第j项因子的熵值：
e
j
=－∑n
i=1
y
ij
ln y
ij
/（ln n）（e
j
＞0）
（4）计算不同时段水资源承载力的综合得分：
s
i
=∑m
j=1
w
j
p
ij
（i=1，2，…，n）（4）
3实例分析
3．1驱动因子的选取流域源区水资源承载力的研究涉及
经济社会和资源环境等诸多因素，影响承载力的驱动因子有
很多，因此在选取时应遵循实用性、可行性和代表性等原则，
尽可能包括影响系统的主要因素。

以赣江源流域赣州地区
为例，从该区的经济、社会、自然和人文等领域出发，进行驱
动因子成分的选取。

根据2001 2009年的《赣州年鉴》、《江西统计年鉴》、
《江西省水资源公报》和《革命老区铸伟业，战略对接大跨
越———辉煌六十年崛起新赣州》等统计资料，分析自然资源、
经济社会以及人口增长等方面的指标，选取以下具有代表性
的16个驱动力影响因子：水资源总量（X
1
，亿m3），总供水量
（X
2
，亿m3），农业灌溉用水量（X
3
，万m3），工业用水量（X
4
，
万m3），居民用水量（X5，万m3），生态环境用水量（X6，万
m3），废污水排放量（X
7
，万t/a），全市年总降水量（X
8
，mm），
总人口（X
9
，万人），全市年生产总值（X
10
，亿元），粮食产量
（X
11
，万t），人均水资源量（X
12
，m3/人），耕地灌溉率
（X
13
，%），城镇化率（X
14
，%），全社会固定资产投资（X
15
，亿
元），全市财政收入（X
16
，亿元）（表1）。

表12001 2009年赣州地区的自然、经济、社会统计指标
Table1Statistics of economic，social and natural resources in Ganzhou region
年份Year
X
1
亿m3
X
2
亿m3
X
3
万m3
X
4
万m3
X
5
万m3
X
6
万m3
X
7
万t/a
X
8
mm
X
9
万人
X
10
亿元
X
11
万t
X
12
m3/人
X
13
%
X
14
%
X
15
亿元
X
16
亿元
2001405．6028．7722900031200317004000345001694．3807．73290．34243．09502179．0917．4372．9223．35 2002417．2025．1516180039500315003800352001967．8818．49321．82225．30509786．6918．19109．2126．04 2003260．5024．1117010047000269003000380001078．5831．20366．39199．07313489．3819．16155．4931．40 2004181．8032．9024100040100308001800321001236．8845．65426．23257．39215088．9919．88158．0038．35 2005320．3324．6318130035200278002000284001636．4845．69500．11259．33378888．3920．03186．0046．71 2006400．0423．6817930029700225002100284001823．6861．15582．73260．53464588．1520．33223．0256．65 2007304．1729．2122860026200309002200338001423．4877．04701．97263．11346885．5320．64300．3274．52 2008315．6529．1721190027400317002200355001558．6888．95834．85269．04355186．2720．85406．00100．13 2009207．0932．5223850027400366002300386001329．2896．99940．02277．71230983．3221．02594．70111．29
3．2计算与分析运用SPSS13．0软件对原始数据进行因子分析处理，得到赣江源流域水资源承载力变化驱动因素的相关系数矩阵（表2）、特征值和贡献率（表3）。

从表3可以看出，前4个因子的累积贡献率达到96．395%，故取前4个因子为主成分进行分析计算。

为了更清楚地解释各变量在主成分因子上的载荷意义，对因子载荷作方差最大化正交旋转，旋转后的载荷矩阵见表4。

载荷矩阵旋转后得到4个主成分因子。

第1主成分在全市生产总值、全市财政收入、全市总人口、全社会固定资产投资、城镇化水平和粮食年产量上具有较大的载荷，综合性较强，充分体现出赣江源流域经济社会的发展水平。

第2主成分在人均水资源量、水资源总量和全市年降水总量上具有较大的载荷，代表着总体的水资源承载量情况。

第3主成分在总供用水量、农业灌溉用水量和居民用水量上具有较大的载荷，主要代表了该区生产生活中水资源量的消耗情况，是技术水平及配置设施在水资源利用方面的表现；第3主成分与耕地灌溉率呈负相关，主要是该区的农田灌溉多数为大水漫灌，在渠系和节水方面的有效利用率较低。

第4主成分在废污水排放量和生态环境用水方面具有较高的载荷矩阵，主要代表着该区生态环境与水环境纳污自净能力。

通过分析可知，4个主成分比较全面地涵盖了赣江源流域水资源承载力的驱动因素，客观、清楚地反映了该区的水资源状况。

所以将赣江源流域水资源承载力的驱动因子归纳为经济社会发展、水资源自然状况、生产生活用水状况及生态水环境纳污自净方面的4个综合指标。

由此，可以通过这4个主成分的得分值F
1
、F
2
、F
3
和F
4
动态地来反映赣江源流域水资源承载力的状况。

图1为2001 2009年该区水资源承载力4个主成分的得分情况。

56241
39卷23期刘慧等基于因子分析和熵权法的赣江源流域水资源承载力研究
表2水资源承载力变化驱动因素相关系数矩阵
Table2Correlation coefficients matrix of driving force factors
指标Index X
1
X
2
X
3
X
4
X
5
X
6
X
7
X
8
X
9
X
10
X
11
X
12
X
13
X
14
X
15
X
16
X
1
1．000
X
2
－0．6401．000
X
3
－0．5270．9321．000
X
4
－0．160－0．351－0．5041．000
X
5
－0．3720．7690．619－0．2841．000
X
6
0．589－0．259－0．2770．2630．1711．000
X
7
－0．3200．3570．2040．1140．6590．4021．000
X
8
0．876－0．422－0．400－0．297－0．1890．372－0．4461．000
X
9
－0．4630．4210．404－0．6470．268－0．7210．120－0．2711．000
X
10
－0．3880．4370．422－0．7160．377－0．5950．205－0．1980．9791．000
X
11
－0．2100．5560．621－0．8370．314－0．617－0．2900．1170．7340．7491．000
X
12
0．993－0．636－0．531－0．078－0．3550．657－0．2960．857－0．558－0．480－0．2831．000
X
13
－0．287－0．377－0．5050．526－0．561－0．543－0．330－0．2330．092－0．073－0．234－0．3051．000
X
14
－0．5340．2910．278－0．4760．045－0．889－0．088－0．3370．9470．8750．670－0．6260．3641．000
X
15
－0．4730．4910．429－0．6050．507－0．4820．382－0．2790．9220．9660．652－0．546－0．1330．7921．000
X
16
－0．3840．4620．436－0．7110．428－0．5460．264－0．2000．9650．9960．731－0．473－0．1150．8410．9711．000
表3相关矩阵的特征根、贡献率和累积贡献率
Table3Eigen values and squared loadings of the principal components
主成分Main component
旋转前Before rotation
特征根
Characteristic
root
方差贡献率∥%
Variance
contributes
累积贡献率∥%
Accumulated
contribution rate
旋转后After rotation
特征根
Characteristic
root
方差贡献率∥%
Variance
contributes
累积贡献率∥%
Accumulated
contribution rate
17．95849．73749．7376．07837．98737．987 23．09319．33169．0693．94124．62962．615 33．00618．78987．8583．26220．38883．004 41．3668．53796．3952．14313．39196．395
表4旋转后的因子载荷矩阵
Table4Rotated component matrix
变量Variables 第1主成分
First
component
第2主成分
Second
component
第3主成分
Third
component
第4主成分
Fourth
component
X
10
0．971－0．1360．1750．079
X
16
0．959－0．1280．2010．139
X
9
0．958－0．2490．103－0．041
X
15
0．908－0．2100．1980．261
X
14
0．880－0．388－0．068－0．259
X
11
0．7160．0590．543－0．410
X
4
－0．716－0．426－0．4940．144
X
1
－0．2270．944－0．206－0．050
X
12
－0．3300．922－0．186－0．022
X
8
－0．0410．911－0．057－0．205
X
3
0．226－0．3270．899－0．037
X
2
0．233－0．4410．8320．108
X
13
0．027－0．482－0．733－0．384
X
5
0．208－0．1450．6770．572
X
7
0．073－0．2550．1350．946
X
6
－0．5980．5240．0540．603 3．3熵权法的应用及分析赣江源流域水资源承载力是受赣州地区自然和人文环境中诸多因子制约的一个有机体，通过对该区水资源系统、自然生态环境系统以及社会经济活动系统之间联系和制约规律的解析，具体地阐述出赣江源流域水资源承载力这个有机体的运行和变化。

借助熵权法来判断各主因子的离散程度，得出评价指标对水资源承载力变化的综合影响程度［10］。

图1赣江源流域赣州地区水资源承载力历年得分Fig．1Water resources carrying capacity score of Ganzhou in past years
将水资源承载力的4个主因子得分值带入式（1） （4），得到各主因子的权重分别为：W
1
=0．263207，W
2
=
0．219067，W
3
=0．206428，W
4
=0．311298。

2001 2009年赣江源流域水资源承载力的综合得分变化趋势图（图2）。

熵权法的综合得分越高，水资源承载力状况将越好。

从图2可以看出，近年来水资源承载力相对稳定，略有增势，2003年有弱化趋势，2004年则达最低限度。

主要原因是2004年大旱，水资源总量为181．10亿m3，较多年平均减少46．05%，年降水量为1236．8mm，较多年平均减少22．75%［11］，使得该年水资源承载力总体承受较大压力而出现降低情况。

2004年之后呈现缓增的趋势，表明该区域水资源向着不断优化的方向发展。

与2009年相比，2001年水资
66241安徽农业科学2011年
源承载力综合得分状况有明显的起色，说明随着社会经济结构的调整及可持续发展战略的指导，该区水资源承载力步入一个良性的发展趋势。

图22001 2009年水资源承载力综合得力
Fig．2Score of water resource carrying capacity from the year of2001to2009
在水资源严重制约地区发展的今天，赣江源流域有着独特的水资源优势，挖掘潜力较大。

从2005年开始，江西省历时4年之久，于2008年底率先在全国范围内完成水量分配任务，从用水总量控制、额定制度管理到明确取用水权和优化配置等方面，摒弃以需定供，粗放管理及“天赐之水”取用不竭等局限性的观念，贯彻一种前瞻性的水资源战略理念［12］。

通过水量分配方案，对赣江流域进行统筹规划，在满足上游地区经济社会和生态环境可持续发展的同时，兼顾中下游地区水量、水质和水环境的协同发展，以推进整个流域水资源承载力的健康发展。

3．4赣江源流域水资源承载力的综合分析图1为影响赣江源流域水资源承载力的主因子得分趋势，正负值反映出与多年平均水平的离差情况。

从图1可以看出，第1主成分代表了经济社会的发展水平，得分值呈明显的上升趋势，展现出该区近年来经济社会发展取得的骄人业绩。

第1主成分在成为该区水资源承载力的主要驱动因子的同时，也对发展提出了挑战。

因为经济的迅猛发展必然会导致水资源消耗量的规模性增加，在现行的技术水平条件下，施加给水资源承载的压力也将不断增大，与之呈负相关的工业用水量的大增，也将进一步使该区水资源承载力面临严峻的考验。

第2主成分代表了赣江源流域水资源的总体状况。

根据得分值所示，水资源在年际间分布不均匀，波动性较大，对承载力的表现值造成不同程度的影响。

该区水资源的主要来源为地表径流，受降水状况的影响较大，尤以2003和2004年为重，异常的大旱使水资源承载力经受了较大的威胁。

从提高防汛抗旱的水利设施和配置方面出发，对水资源进行宏观调控以维护流域自身的环境稳态，避免上游流域挤占、抢夺下游用水及工业挤占生活、生态和农业用水等情况，从而维持该区水资源承载力的发展稳态。

第3主成分代表了农业和生活水量的耗用情况。

虽然2001 2003年得分值波动较大，但总趋势是在接近一个均值的水平。

随着连年来人口的增加，耕地灌溉和生活用水量也持续增加，尤其是2004年大旱期间，水资源总量较少，农业和生活等方面的大耗水量凸显出来，成为削弱水资源承载力的重要因素。

因此，在调整社会产业结构方面，须下大力气，合理安排三产业的比重，符合社会进步发展的要求，以减轻水资源承载的压力。

第4主成分代表着生态水环境系统的纳污自净能力，是河道水系健康承载的重要指标。

近年来，由于河岸植被破坏严重、森林资源过载砍伐及生产生活废污水的大量排放，使得赣江上游水土流失加剧，泥沙大量下泄，掩埋农田，淤塞江河湖库等，对中下游河道及鄱阳湖生态区等地构成严重的威胁，并且出现赣江流域两头重中间轻的污染状况，诸多支流河段的氨氮、总磷含量超标［16］。

因此，生态水环境系统也是流域水资源承载力综合评价的重要指标，只有在人与自然和谐相处的前提下，人类社会的可持续发展才会有坚定的基石。

4结语
（1）运用因子分析法，得出影响赣江源流域水资源承载力的4个主成分因子，综合了在经济社会、人口现状、自然水资源及生态环境方面的多项指标；通过熵权法，取得近年来水资源承载力的一个良性趋势，以图表的形式将潜在的状况直观地展现出来，以期为该区经济社会和水资源的优化配置提供有效参考。

（2）针对目前该区的发展态势，提出一些提高水资源承载力的切实对策。

首先，从水资源自身即承载主体出发，可通过雨水集存，改善水利工程等“开源”措施，增大水资源的可利用量。

其次，从水资源承载的客体即广义的与水资源相关的经济、社会、生态和环境系统的指标出发，通过控制人口增长速度，调整产业结构规模，普及推广节水和污水资源化技术以及出台相应的政策法规等措施，降低水资源的使用量即“节流”，降低经济社会发展形成的需求压力［1］。

在江西省推出水量分配的实际情况下，从节流方面采取措施将会是更加主动，合理的途径。

（3）赣江源流域水资源承载能力是该区可持续发展过程中的重要组成部分。

总体看来，随着社会经济的发展，技术水平的不断进步，人们对水资源忧患意识的提高以及水资源优化政策的出台，该区水资源承载力将逐渐趋于健康稳定状态，在推进该区综合实力提升方面将发挥不可估量的作用。

参考文献
［1］朱启林，张海涛，游进军，等．提高水资源承载能力途径的研究［J］．南水北调与水利科技，2009，7（1）：31－35．
［2］赵军凯，张爱社．水资源承载力的研究进展与趋势展望［J］．水文，2006，26（6）：47－50，54．
［3］周琳，金辉．主成分分析法在江门市水资源承载力研究中的运用［J］．人民珠江，2007（5）：39－42．
［4］程美家，韩美，杜剑．基于主成分分析法的山东省水资源承载力评价［J］．资源开发与市场，2009，25（5）：410－412．
［5］陈腊娇，冯利华，毛小军．主成分分析法在水资源承载力影响因子评价中的应用［J］．水利科技与经济，2006，12（6）：362－364．
［6］张勇，刘娟，吕建树．山东省水资源承载力的主成分分析和熵值评价［J］．山东水利，2009（11）：56－59．
［7］杨娜，李慧明．基于因子分析与熵值法的水资源承载力研究———以天津市为例［J］．软科学，2010，24（6）：66－70．
［8］王好芳，郭乐，窦实．基于主成分分析的东江流域水资源承载能力评价［J］．水文，2008，28（4）：16－19．
［9］周亮广，梁虹．基于主成分分析和熵的喀斯特地区水资源承载力动态变化———以贵阳市为例［J］．自然资源学报，2006，21（5）：827－833．［10］邵磊，周孝德，杨方廷，等．基于主成分分析和熵权法的水资源承载能力及其演变趋势评价方法［J］．西安理工大学学报，2010，26（2）：170－176．
（下转第14277页）
76241
39卷23期刘慧等基于因子分析和熵权法的赣江源流域水资源承载力研究
＞Ni＞Mn＞Zn＞Pb。

推测该结果与各地工业类型、发达程度有关。

表6婺江表土沉积物重金属地质积累指数I
geo
Table6I
geo
in each sampling site of Wujiang River
采样点Sampling point
I
geo
Cr Mn Ni Cu Zn Pb
平均值
Mean
武义江2．261．442．253．910．790．921．93 Wuyijiang
义乌江3．601．363．518．121．141．243．16 Yiwujiang
环五百滩4．911．253．727．521．170．903．25 Huanwubeitan
金华江3．911．323．827．691．251．143．19 Jinhuajiang
工业过程中广泛使用重金属元素，工矿企业将未经严格处理的废水直接排放，使其周围土壤易富集高含量的有毒重金属。

企业排放的烟尘、废气中也含有重金属，并最终通过自然沉降和雨淋进入土壤。

矿业和工业固体废弃物在堆放或处理过程中，由于日晒、雨淋、水洗等作用，使重金属极易移动，以辐射状、漏斗状向周围土壤扩散，固体废弃物也可以通过风的传播而使污染范围扩大。

综上，重金属Cu和Pb在义乌江侧污染较为严重，而污染元素Pb、Cu、Zn等来源主要与汽车尾气排放及汽车轮胎磨损产生的大量含重金属的有害气体和粉尘的沉降有关。

义乌江上接东阳江，东阳市3家污染大户（化工厂、化肥厂和农药厂）、沿江几十家企业以及义乌江上游分布10余家电镀、印染企业，除此，义乌市高密度车流量与人流量均对义乌江水质产生重要影响；Mn污染侧重于武义江，其上游的永康、武义是全国最大的五金产品生产基地和集散中心。

除此还需考虑Mn、Fe元素本身土壤中含量就较高［13］；婺江流域干流金华江的重金属一部分来自于义乌江和武义江，也有一部分来自当地，该流域的主要工业有通用设备制造业、金属制品业、纺织业、有色金属冶炼及压延业、纺织服装鞋帽制造业，其中，重工业比重多于轻工业，普遍存在企业规模小、产业延伸能力不强、产业链短、高耗能、高污染等主要问题，另外还有多家建材企业和多种设备工程公司，如华海工业设备安装有限公司、华泰工业设备工程有限公司。

3结论
（1）婺江流域义乌江和武义江磁性矿物类型单一，主要由亚铁磁性矿物贡献，而不完整反铁磁性矿物含量低；干流金华江段存在多相磁性矿物，亚铁磁性矿物主导。

样品磁学参数指示金华江与武义江的晶粒特征主要以稳定单畴颗粒（SD）为主，义乌江的磁性矿物主要由多畴（MD）、超顺磁晶粒（SP）等贡献，婺江流域段多属于由粗多畴颗粒和稳定单畴颗粒组成的混合成分磁性物质。

（2）通过对婺江表层沉积物的部分重金属含量与磁性参数之间的相关性分析，发现Cr、Ni、Cu、Zn与X呈较强的正相关关系，其中，Cu与磁化率的相关性最强，相关系数达0．63。

表土沉积物磁化率与重金属Cr、Ni、Cu、Zn均达到0．01极显著的正相关水平，与Mn也达到了0．01极显著负相关水平，与Pb的相关性不显著。

X
fd
与重金属Mn和Pb达到0．01极显著正相关，与Cr达到了0．05显著负相关水平。

表明磁化率对婺江表层沉积物重金属Cr、Ni、Cu、Zn、Mn等具有较好的指示作用。

（3）通过重金属的分布与对重金属污染评价，推测对其污染来源主要与河流两岸工业企业污水排放、汽车尾气排放、汽车轮胎摩损、大气烟尘沉降以及居民生活污染等因素有关，而前3项是重金属污染的主要来源。

参考文献
［1］俞立中．环境磁学在城市污染研究中的应用［J］．上海环境科学，1999，18（4）：175－178．
［2］R·汤普森，F·奥尔德费尔德．环境磁学［M］．严尧基，等，译．北京：地质出版社，1995：36－49．
［3］俞立中，张卫国．利用磁信息研究潮滩重金属污染的探讨［J］．环境污染治理技术与设备，1993，5（1）：37－44．
［4］张卫国，俞立中．沉积物磁性测量对铁还原的指示及其在重金属污染研究中的应用［J］．科学通报，1998，43（19）：37－43．
［5］张卫国，俞立中，HUTCHINSON S M．长江口南岸边滩沉积物重金属污染记录的磁诊断方法［J］．海洋与湖泊，2000，31（6）：616－623．
［6］THOMPSON R．Modeling magnetization data using SIMPLEX［J］．Phys Planet Int，1986，42：113－127．
［7］陈满荣．长江口潮滩沉积物重金属污染研究［M］．上海：华东师范大学出版社，2009：141－143．
［8］吕达．崇明岛沉积物重金属污染特征及磁学研究［D］．上海：华东师范大学，2007：7－8．
［9］CANBAY M，AYDIN A，KURTULUS C．Magnetic susceptibility and heavy-metal contamination in topsoils along the Izmit Gulf coastal area and IZAY-TAS（Turkey）［J］．Journal of Applied Geophysics，2010，70：46－57．［10］董岩翔，郑文，周建华，等．浙江省土壤地球化学背景值［M］．北京：地质出版社，2007：141－142．
［11］单红丹，卢升高．火电厂粉煤灰的矿物磁性及其环境意义［J］．矿物学报，2005，25（2）：141－146．
［12］MULLER G．Index of geo accumulation in sediments of Rhine River［J］．Geo J，1969，2：108－118．
［13］王虹艳．丽水市城市土壤重金属污染空间分异及磁学响应研究［D］．杭州：浙江大学，2008：
檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪檪
29－32．
（上接第14267页）
［10］邵磊，周孝德，杨方廷，等．基于主成分分析和熵权法的水资源承载能力及其演变趋势评价方法［J］．西安理工大学学报，2010，26（2）：170－176．
［11］吴速英，许洪胤，钟洪鸣，等．赣江流域水资源承载能力初步分析［J］．中国资源综合利用，2009，27（2）：28－30．
［12］龚向民，李昆，刘筱琴，等．赣江流域水土流失现状与发展态势研究［J］．人民长江，2006，37（8）：48－50．［13］李昆，谭振江，周家丽，等．人类活动对赣江流域水环境的影响探讨［J］．水资源研究，2007，28（3）：9－11，28．
［14］章文波，陈红艳．实用数据统计分析及SPSS12．0应用［M］．北京：人民邮电出版社，2006．
［15］倪青山，刘小丹．SPSS应用实验教程［M］．长沙：湖南大学出版社，2007．
［16］傅春，姜哲．赣江水源地保护的若干问题及对策［J］．水资源与水工程学报，2008，19（5）：58－60，64．
77241
39卷23期方慧青等婺江表层沉积物重金属元素的磁学响应。