湖泊生态承载力研究

　　根据MIEP，在约束条件(2)下信息熵SJ最大化，即利用Lagrange乘算法使式(1)在式(2)下取极值，由泰勒展开式得：

　　式中，系数σ为由Lagrange优化决定的参数.由此可推出动力学演化方程：

　　式中，Fk和Sk反映指标或组元间的相互作用，aki表示组元xi间的连接权值，ξk(同ξ)是xi的宏观组合模式的结构参量，ξ值代表生态承载力水平，若指标数据xi都是正向的，则ξ值越大，生态承载力越高.

　　该模型提供了一种描述世界万象的新方法，这样的“有生于无”的过程通过方程(4)这个动力系统来实现.一个特征值λk对应着一个可能的宏观结构模式ξk，λk>0的情况对应着活跃的组合模式，将在竞争的过程中获胜进而决定生态承载力的宏观结构，λk<0的情况对应着不活跃的组合模式，在竞争过程中被淘汰.由此揭示生态承载力结构模式的变化特征，这样我们就建立了基于信息熵的生态承载力评价模型.

　　自组织特征映射神经网络(SOM)是一个与生态承载力等价的动力系统，每个神经元相当于一个元胞，神经元之间存在复杂的相互作用.将生态承载力的指标组元映射到SOM网络中，激发某个神经元产生了权值aki，式(1)中的信息会发生变化，通过竞争继续激发新的神经元，以此迭代可以呈现出生态承载力的动态演化过程，同时完成式(4)的求解.在这样的复杂系统中需要力所能及地考虑所有联系，因此，需要利用MATLAB得到其演化的动力学过程和生态承载力的表征值ξ.

　　3 洞庭湖生态经济区生态承载力评价

　　3.1 数据来源与计算

　　本文的研究区域包括洞庭湖生态经济区的常德、益阳、岳阳3个中心城市及长沙望城区共25个县(市、区)，对其2008—2012年生态承载力进行分析，其中，社会经济与资源环境类指标的数据来源于2008—2012年各地的《统计年鉴》、《国民经济与发展统计公报》、《环境状况公报》和《国民经济和社会发展统计公报》等，土壤与土地类指标数据由其他论文和网络数据综合得到.

　　评价指标中涉及的经济数据包括人均GDP、万元GDP能耗和万元GDP水耗.考虑到经济发展中价格变化的影响，以现价人均GDP、居民消费水平绝对额和保险福利费用不能直接对比，故先将上述经济数值折算成2008年不变价，然后利用不变价的经济数值进行比较.并且由于原数据受单位和数量级的影响，指标耦合的差异性较大，故需要对原始数据进行无量纲化，从而保持原始数据的一致性，在消除量纲、数量级影响的同时，保留各指标的变异程度.对生态城市建设有负影响的逆向指标先取倒数使其正向化，然后采取式(6)所示极值法进行归一化处理后，在冯綮一(2010)开发的MFP模拟程序包的基础上进行程序的编写，借助MATLAB进行计算.

　　式中，i=1，2，…，m;Ci，max为该指标最大值.

　　3.2 结果与分析 3.2.1 洞庭湖生态经济区生态承载力综合评价

　　将指标层的数据标准化，采用3×4网格数进行400步演化模拟，计算洞庭湖生态经济区2008—2012年的ξ值(ξ1~ξ5)，结果如图 2所示.在演化过程中，指标组元相互关联、耦合，不断进行竞争和自组织作用，从而产生反映其生态承载力水平的获胜模式，获得了每一组给定x1，x2，…，xn数据下的ξ值.可以看出，前150步模拟ξ值波动较大，之后逐渐趋于稳定，达到相对稳定的模式结构，此时的ξ值可代表生态承载力水平.

　　图 2 2008—2012年洞庭湖生态经济区生态承载力动态演化过程

　　将模拟得到的表征生态承载力水平的每一组给定x1，x2，…，xn数据下的ξ值整理后作折线图，从图 3以看出，洞庭湖生态经济区2008—2012年期间生态承载力呈总体上升的趋势，其中，2008年和2009年两年生态承载力水平基本一致，2009—2010年呈现较明显的上升，而2010—2011年呈现下降趋势，2011—2012年生态承载力水平大幅度提高.说明近年来对洞庭湖生态经济区的生态承载力建设效果显著，该区域生态系统组元间的耦合关系越来越紧密.但2011年生态承载力水平的异常情况需要进一步模拟分析，从而为该地区下一步的规划建设提出针对性的建议.

　　图 3 洞庭湖生态经济区生态承载力综合评价结果

　　3.2.2 生态承载力准则层评价研究

　　为了深入了解洞庭湖生态经济区生态承载力发展状况，根据其各方面的发展状态提出针对性的对策建议，提高生态承载力水平，可以基于准则层进行分析.将标准化的数据依次代入SOM网络中进行分类分析模拟，借助MATLAB，分别演化出能够反映该地区生态弹性(B1)、资源环境承载力(B2)、环境压力(B3)和人类潜力(B4)的指标特征值，从而得到洞庭湖生态经济区在这4个方面的发展趋势(图 4).

　　图 4 洞庭湖生态经济区2008—2012年准则层指标发展趋势

　　洞庭湖生态经济区生态弹性(B1)的发展状态分两个阶段：2008—2009年呈现下降趋势，2009—2012年呈稳步上升趋势.可见前期对生态系统的利用和改造忽略了生态弹性方面的发展，导致生态系

　　统自我调节能力下降.2009年以后政府开始认识到生态弹性对整个生态系统承载力的重要性，使其发展水平慢慢上升.资源环境承载力(B2)的发展很不稳定，特别是2011年急剧下降，可见洞庭湖生态经济区的资源状况易受影响，需要政府在这方面加强监管，采取一定的措施促进其平稳发展，以免对生态承载力综合水平造成较大的影响.环境压力(B3)在这5年间变化不大，说明该方面各指标的发展有上升有下降，在一定程度上造成了整体水平平缓的现象.人类潜力(B4)的发展逐年上升，这与经济发展和科技进步密不可分，还有政府在环境保护方面的投资和政策支持，人类环保意识的加强，都不同程度地促进了生态承载力的发展.

　　3.2.3 ξ值各分量雷达图表征的生态承载力的发展演化

　　MIEP模型将区域生态系统看作是一个充满耦合的复杂网络，连接权值表示网络中各指标与周围环境交换广义流能力的大小，对系统发展所做贡献能力的大小.为了分析洞庭湖生态经济区生态承载力发展演化状况及各指标对系统所做贡献的多少，用aki×xi描述ξ的形态，据其作雷达图，将i轴的值用直线连接起来，得到的平面图形即可表示各年生态承载力状态特征值ξ的结构形态特征(图 5).

　　图 5 洞庭湖生态经济区生态承载力的演化过程

　　整体来看，从2008年至2012年，雷达图覆盖的面积越来越增大，说明洞庭湖生态经济区生态承载力发展过程的信息交流越来越强烈，系统结构逐渐完善.2008年到2009年覆盖面积相似，说明这两年间各指标的结构状态没有发生较大变化.相比于前两年，2010年到2012年万元GDP能耗(C16)、农药施用强度(C17)、环保支出占财政比重(C18)和人均国内生产总值(C19)的较大发展使其对生态承载力的贡献变大，指标均衡度提高，促进生态承载力发展结构的完善.2010年和2012年洞庭湖生态经济区的地表水资源量(C9)、地下水资源量(C10)和人均水资源量(C11)三方面指标的较高水平使该地区生态承载力的系统结构有了明显的改善，各指标发展更均衡，生态承载力水平大幅提高.这与当地的实际情况相符，2010年和2012年洞庭湖生态经济区年人均水资源总量分别为2430.914 m3和2480.617 m3，远高于其他年份，促使这两年洞庭湖生态经济区的生态承载力达到两个较协调稳定的发展状态，从而造成2011年生态承载力水平异常发展的现象.

　　表 2 模型的预测值与实际值的对比

　　3.3 生态承载力方面的预测

　　由于区域生态系统是由信息作用耦合在一起形成的稳定结构模式，由其结构参量控制，如果掌握组元与ξ之间作用关系的规律后，MIEP模型也可预测生态系统的发展演化趋势.调控方法是在xi-ξ关联图上作相应延展或插值，以此再现组元相互作用的动力学过程和结构模式的演化.

　　为找出与洞庭湖生态经济区生态承载力发展变化关联性较高的评价指标，将根据MIEP模型的输入与输出做趋势回归分析.通过xi-ξ间的拟合关系，得到xi与ξ值的决定系数R2，发现趋势回归关联性最好的3个指标分别为地下水资源量(C10)、人均水资源占有量(C11)和教育投资占GDP比重(C22).将地下水资源量2013年的实际值带入拟合方程求得ξ为16.0059，然后根据其求得其他指标的理论值并与实际值对比后，结果如表 3所示.结果表明，指标的实际值与理论预测值之间的误差均比较小，在一定程度上证明了MIEP模型预测复杂系统的潜力.在实际工作中，通过该模型可对未来发展进行预测分析，为管理决策人员提供参考.

　　3.4 MIEP模型与主成分分析法的比较

　　MIEP模型给出的是ξ值，而主成分分析法(丁扬，2014)给出的一个指数值，两种评价方法虽然不能进行数值比较，但它们对洞庭湖生态经济区2008—2012年生态承载力发展评价结果的趋势基本相同，201年生态承载力水平剧烈下降，2011—2012年间逐渐回升(图 6).二者的主要区别是本文使用的MIEP模型显示，2012年生态承载力水平超过了2010年，而主成分分析法则相反.主成分分析法依靠每个主成分的贡献率作为权重构造评价函数，但通常各主成分特征值的相差较大，使其方差贡献差别较大，无法使所有指标完全综合，主成分权重的设定带有强烈主观性，易使个别指标成为影响评价结果的决定性因素，耦合性差，因此，造成生态承载力水平受到水资源因素影响震荡较大.

　　图 6 MIEP模型与主成分分析法的比较

　　本文所使用的基于MIEP的系统动力学模型，不必预先确定各个指标的权重，在模拟过程中能自动调整各指标数据之间的耦合关系，不断进行自组织作用，克服了其他方法权重确定的主观性的缺陷;而且具有鲁棒性，更适用于复杂系统的模拟分析，可以呈现出各种复杂系统的动态演化过程.具体参见 污水处理技术资料或污水技术资料更多相关技术文档。

　　4 结论

　　本文从复杂系统组元相互作用的角度建立了评价生态承载力发展水平的MIEP模型，并运用该模型对洞庭湖生态经济区2008—2012年的生态承载力发展情况进行了评价.结果表明，在此期间生态承载力呈总体上升的趋势，其中，2008年和2009年两年生态承载力水平基本一致，2010年呈现较明显的上升，而2011年的生态承载力剧烈下降后在2012年呈现大幅回升.ξ值各分量的雷达图揭示了洞庭湖生态经济区复杂系统的演化结构，根据其变化判断各指标贡献程度的变化，通过对指标层雷达图及准则层的分析，发现大幅波动的原因是2010年和2012年水资源方面的巨大优势使这两年资源环境承载力水平较高，从而引起生态承载力综合水平的大幅提升.根据MIEP的工作机理，对该地人均水资源占有量(C11)和教育投资占GDP比重(C22)在2013年的发展情况进行预测，发现误差较小，显示了其预测生态城市建设发展的较大潜力.最后将本文的评价结果与主成分分析法结果进行了比较，发现MIEP模型更适用于复杂系统的分析研究，更能表现其动态演化的本质.

　　5 建议

　　由于水资源指标更易影响洞庭湖生态经济区生态承载力综合水平，因此，当地政府应着重提高其水资源承载力，如改变传统的灌溉方式，加强节水的宣传教育和推行节水工艺等.传统工农业对环境影响较大，易造成水和大气污染，因此，转变经济发展方式有助于生态承载力的提高.同时，还要加强环保监管力度，政策落实到位才能从源头上遏制污染，提高环境质量，促进生态承载力的发展.