随着社会经济的快速发展，人类对资源利用与环境保护之间的矛盾日趋尖锐。环境变化与社会经济之间的关系日益受到重视^{[1, 2, 3]}，国外学者经常通过对不同地区、不同污染物的实证研究来检验环境库兹涅茨曲线(EKC)是否存在，并将经济增长与生态环境之间是否存在EKC曲线作为论证工业化、经济增长与生态环境能否协调发展的依据^{[4, 5, 6]}。而流域水环境作为人类生存和发展的重要依托，自20世纪90年代以来，以流域为单元开展社会经济发展与水环境变化相互作用的研究方兴未艾^[7]。

国内学者近年来也开始关注水环境的研究。部分学者从战略的高度提出了我国流域水环境保护与社会经济协调发展等问题^{[8, 9]}；大量学者对于水环境的演变规律进行了分析^{[10, 11, 12, 13, 14, 15]}，部分学者对于水环境的演变成因进行了定性^{[10, 11, 12, 14]}或是定量^{[4, 7, 13, 16, 17, 18]}探讨。从国内研究现状来看，应加强研究方法的多元化以及研究的多视角化，从水环境的社会经济影响机制、发展趋势全方位进行探讨。

辽河流域是我国七大流域之一，位于我国重工业发达的东北地区，流域内城市群集中、人口密集、工农业发达，该区域曾为全国的经济发展作出过重大贡献，多年的工农业生产导致流域水环境严重受损，辽河被列为我国重点治理的“三湖三河”之一^[19]，其所面临的人口、资源、环境与发展问题比其它区域更加深刻而复杂。近年来，针对辽河流域水环境的研究日益增多。特别是在水环境演变趋势、水环境承载力、生态影响以及管理对策方面出现了较多成果^{[19, 20, 21, 22, 23, 24, 25, 26, 27, 28]}。从总体来看，这些研究工作对辽河流域水环境保护及水污染治理起到了重要作用，但未能将各类社会经济活动与流域的污染压力结合起来，没有考虑流域内社会经济分布与污染产生的关系。而这些都是流域水环境管理中不可忽视的视角。

基于以上原因，本文以辽河辽宁省段为例，尝试运用多元线性回归分析和环境库兹涅茨曲线模拟，定量研究分析辽河流域自2001—2010年10年间水环境演变与流域内社会经济发展之间的关系，并作出一定的发展判断和分析。

辽河流域是我国七大流域之一，包括辽河水系和大辽河水系。辽河全长1390公里，流域面积21.9万平方公里，由西辽河、东辽河以及发源于吉林、内蒙古的招苏台河、条子河等支流在辽宁省境内汇合而成，于盘锦市入海，河流干流主要在辽宁省境内。

本文研究范围为辽宁省境内的辽河流域，主要由辽河、太子河、浑河、大辽河4条干流和40多条一级支流组成，流经沈阳、鞍山、抚顺、本溪、盘锦、营口、辽阳、铁岭八地市，其干流由辽河水系和大辽河水系组成，分别由盘锦的双台子河和营口的大辽河入海。

辽河流域辽宁省境内是经济较发达的工业集聚区和都市密集区。流域内社会经济密度较大，社会经济相对活跃，辽河流域承载了超过全省半数以上的社会经济产值。以沈阳为核心的“五带十群”发展规划指出至 2015 年地区生产总值年均增长达14%以上，规模以上工业增加值年均增长16%，城镇化率要达到80%以上，要支撑如此大规模快速的地区发展，环境压力巨大^[22]。

2 001年至今，辽河流域辽宁省段内经济高速增长。2010年全流域地区生产总值达12284.41亿元，是2001年的3.93倍，人均国内生产总值是2001年的将近4倍。三次产业发展中，2010年生产总值分别是2001年的3.06、4.48、3.82倍(图 1)，由此可以看出，区域发展过程中，二产增长速度较快，尤以工业增长速度明显，为4.53倍。同时，2001—2010年辽河流域人口与GDP发展趋势类似，也呈现增长趋势，由2247.25万人增长到2302.15万人(图 2)。

图 1 辽河流域内2001—2010年各产业总产值 Fig. 1 Total production value of industries 2001—2010 in Liaohe River watershed

图选项

图 2 辽河流域内2001—2010年人口变化 Fig. 2 Population sorts of 1990—2000 in Liaohe River watershed

图选项

在最近的10a中，辽河流域水环境质量变化较大(图 3)。辽河水系辽宁省段内包括辽河、浑河、太子河和大辽河4条主要河流。2002年辽河铁岭、沈阳、盘锦3个河段均为劣Ⅴ类水质，污染严重。主要污染物是化学需氧量、生化需氧量和氨氮，分别超过国家Ⅴ类水质标准0.8倍、0.1倍和0.8倍。由于境外污水排入，对辽河铁岭段水质影响较大，浑河、太子河和大辽河的18个监测断面中，Ⅳ类和劣Ⅴ类水质分别占27.7％和66.7％。其中浑河抚顺段为Ⅳ类，沈阳段为劣Ⅴ类；太子河本溪段为Ⅴ类，辽阳段为Ⅳ类，鞍山段为劣Ⅴ类；大辽河盘锦段、营口段均为劣Ⅴ类水质。氨氮污染突出，最高断面超过国家Ⅴ类水质标准6.9倍。各断面中，沈阳七台子、于家房、本溪兴安和盘锦三岔河断面污染严重。大凌河水系大凌河是全省水质污染最重的河流。全河段化学需氧量、生化需氧量浓度均值分别超国家Ⅴ类水质标准3.6倍和7.4倍。

图 3 辽河水体2001—2010年平均水质结果 Fig. 3 Annual average on water quality of 2001—2010 in Liaohe river (mg／L)

图选项

2 010年，辽河流域水质相对好转，26个干流断面的化学需氧量平均浓度为18mg/L，辽河、浑河、太子河化学需氧量年均浓度在12—24 mg/L之间，为2001年来最低值。干流断面化学需氧量年均浓度全部符合Ⅴ类水质标准。但是氨氮污染仍然较重，57.7%的干流和59.5%的支流断面氨氮超Ⅴ类水质标准。

运用多元线性回归分析方法，建立辽河流域内2001—2010年经济发展与水质之间的相关关系模型。其中，社会经济指标为自变量，水质指标为因变量。选用的社会经济指标包括第一产业生产总值(FDP)、第二产业生产总值(SDP)、年末总人口(POP)，数据资源来源于辽宁省统计年鉴；考虑到辽河现有资料的实际情况，选择研究期间逐年的化学需氧量(COD_Mn)、氨氮(NH₃-N)等年平均值作为水质评价指标，数据来源于环境统计报表。

多元线性回归是研究多个变量之间关系的回归分析方法，根据多个自变量的最优组合建立方程来预测因变量。多元回归分析的模型可以表达为：

点击浏览公式

环境库兹涅茨曲线是通过人均收入与环境污染指标之间的演变模拟，说明经济发展对环境污染程度的影响，环境库兹涅茨曲线的含义是沿着一个国家的发展轨迹，尤其是在工业化的起飞阶段，不可避免地会出现一定程度的环境恶化；在人均收入达到一定水平后，经济发展会有利于环境质量的改善^[30]。

本研究对辽河流域的水环境与社会经济之间的关系进行EKC曲线模拟。

这里采用经济增长指标(x)：考虑指标的代表性，选取人均GDP作为模型中经济增长指标。

水环境指标(y)：根据辽河流域辽宁段区域目前发展阶段和生态环境问题的突出表现，分别选取COD_Mn、NH₃-N作为模型中环境污染指标。

将数据输入到模型中，得到COD_Mn多元回归模型结果表 1—表 3。

由表 1中的R²以及修正的R²值可以看出回归关系较好。表 2方差分析结果表明，当回归方程中包含自变量第一产业生产总值与第二产业生产总值时，其显著性概率值Sig．小于0.05，拒绝总体回归系数均为0的原假设，因此回归方程应该包含这两个自变量(年末总人口由于Sig．大于0.05，因此舍弃)。

同样的方法可以得到氨氮的方差分析和回归分析表。但由于Sig．大于0.05，表明氨氮与各指标间的关系不显著，因此略掉。

根据以上COD_Mn的多元回归分析结果，建立的线性相关关系如下：

COD_Mn=116.08-0.13FDP+0.002SDP

模型的输出结果表明：①2001—2010年期间第一产业总产值、第二产业总产值对水质污染因子COD_Mn的影响效果是显著的。可以推断流域内社会经济所带来的污染压力可能还会逐年增加，应加强管理与治理；② 氨氮与社会经济因子间建立多元回归方程时R²较小，表明其与各因子间的关系相对较弱，这可能与影响氨氮的因子众多有关，这也是辽河流域近年来氨氮指标一直居高不下的重要原因。

根据上文分析，利用2001—2010年的人均GDP以及COD_Mn、NH₃-N数据，可以得到COD_Mn、NH₃-N与人均GDP之间拟合方程，各个方程所描述的经济增长与水环境环境污染之间关系曲线图分别如图 4和图 5所示。

图 4 辽河流域CODMn浓度与人均GDP关系 Fig. 4 Relationship between per capita GDP and thickness of CODMn in Liaohe River watershed

图选项

图 5 辽河流域NH3-N浓度与人均GDP关系 Fig. 5 Relationship between per capita GDP and thickness of NH3-N in Liaohe River watershed

图选项

模拟结果表明：

(1)辽河辽宁段的COD_Mn浓度与流域人均GDP存在N型函数关系，曲线经历下降、上升、再下降的多个阶段；但趋势线呈现出一定的库兹涅茨曲线倒U型特征，在人均GDP达到27716元时出现拐点。总的来说，辽河流域在COD排放控制上已经开始逐步越过这个拐点位置，目前正处于由上升向下降转化的阶段。这表明伴随着经济的不断增长，环保技术的应用、产业结构的变化，水环境质量开始呈现出局部改善的趋势。

(2)辽河辽宁段的NH₃-N浓度与流域人均GDP之间，没有显示存在典型EKC曲线倒U型的特征。目前，随着人均GDP的增长，NH₃-N浓度呈现出逐步下降的趋势，但下降趋势不大，结合辽宁省近年来发布的水环境公报，综合分析其原因，主要是由于造成面源污染的因素众多，较难控制。

以流域为单元开展社会经济发展与环境演变相互作用及调控研究，是实现可持续发展的有效途径^[7]。辽河流域水环境问题的发生与发展同流域内经济发展阶段及产业结构的调整有很大的关系。通过辽河流域辽宁段的社会经济指标与水质指标的多元线性回归分析，可以发现，辽河的水环境质量的演变，同流域内以工业为主的产业结构和流域人口的增加有关。

环境库兹涅茨曲线模拟研究表明：在辽河流域经济增长与水环境环境质量变迁关系的曲线中，随着区域经济水平的提高，COD_Mn、NH₃-N污染在综合治理下开始出现一定程度的改善，但不完全符合典型库兹涅茨曲线特征，特别是氨氮变化趋势不明显。

由2001—2010年流域内经济发展与水环境演变的关系，并结合前人的研究成果可以推断，未来辽河流域的污染物来源中，农业及农村生活污水所占的比重可能将持续增加，并成为主要的污染来源^[13]。因此，在开展针对工业点源污染而采取的达标排放的同时，还必须加强农业面源污染的治理力度，“点源与面源相结合”，以利于辽河水环境的改善。

本文运用多元回归分析、环境库兹涅茨曲线模拟分析了辽河水环境演变与流域内社会经济发展的相关关系，为辽河流域水环境污染控制与治理提供一定的科学依据。流域水环境及其治理问题，受多种因素综合作用，在本研究中仅仅考虑了流域内社会经济的发展状况对辽河水质的影响，对气候变化、流域土地利用类型改变、政府行为及环境保护政策等方面的作用没有探讨，这可能导致研究结果的分析不够全面。此外，环境库兹涅茨曲线仅反映了多种环境-收入理论关系的一种形态，且更适用于流量污染物和短期的情况，其多形态的环境-收入关系的理论基础还需要深入研究和探讨。