论文部分内容阅读
摘要：立足于系统论和协同学理论，应用计量经济分析方法将耕地资源视作由经济、社会、生态环境3个子系统相互关联、相互影响、相互制约所构成的复合系统，从三维视角研究耕地资源系统在一定时空范围内，经济-社会-生态环境子系统的有序度、相互胁迫关系与变化规律。研究表明，耕地资源系统有序度是3个子系统的综合反映，20年间经济与生态环境子系统有序度变化幅度较大，各子系统之间存在交叉相对变化的趋势，相对来说，社会子系统有序度变化不明显。与此同时，沧州市耕地各子系统间确实存在一定的交互胁迫关系，在1994—2013年间，除泊头市外，经济-社会子系统表现出明显的协同发展关系，经济社会与生态环境子系统之间的胁迫关系较明显。本研究为不同学科间的交叉研究提供了新的理念，并为耕地资源可持续发展研究提出了新思路。
关键词：耕地资源系统；系统有序度；交互胁迫；沧州市；可持续发展
中图分类号： F301.24文献标志码： A
文章编号：1002-1302（2017）05-0256-09
耕地作为宝贵资源，是由自然、经济、社会、生态、环境等各个子系统相互关联、相互制约所构成的复合系统的综合体。在某一耕地，经济、社会、生态环境子系统之间能够相互配合、相互促进、协调发展，耕地资源系统就会实现健康、有序、可持续发展的良性循环；反之，耕地资源系统就会陷入生态环境恶化、社会经济紊乱的恶性循环。当前，随着城镇化步伐的加快，建設用地不断占用耕地，尽管我国耕地保护制度较为严格，但耕地数量锐减、耕地质量下降的问题依然存在[1]；与此同时，面对人地关系紧张与粮食安全问题的严峻形势，为满足粮食需求和获得更大的产出效益，耕地利用集约化显著，大量的化肥、农药投入在一定程度上产生了一系列耕地生态安全问题[2-4]。
面对耕地生态保护与社会经济发展之间的博弈问题，将耕地資源作为经济-社会-生态环境的复合系统，探讨各子系统之间的相互胁迫关系与作用机制，探寻耕地资源系统经济-社会-生态环境的协调发展途径，对于实现研究区乃至其他地区耕地资源系统的可持续发展具有重要意义；同时，耕地资源系统的相关研究也将作为土地利用修编的重要命题之一，为政府部门制定和改进相关政策建议提供一定的理论依据。
1研究区概况及数据来源
1.1研究区概况
沧州市位于河北省东南部，环渤海湾经济区，地处 116°10′～117°44′E，37°53′～38°51′N，水文地质条件较差，中低产田和盐碱地分布较广，是河北省粮、棉、油集中产区之一。截至2013年末，全市常住人口为730.95万人，耕地面积为 78.69×104 hm2，沧州市耕地质量总体水平不高，且部分地区生态环境脆弱，土壤盐碱化问题突出。
1.2基础数据来源
本研究数据资料包括耕地数据和社会经济数据2个方面，耕地数据主要来自1994—2013年《河北省土地调查统计年鉴》，社会经济数据来源于1994—2013年《河北农村统计年鉴》及沧州市各县（市）国民经济统计年鉴。
2耕地资源系统有序度的内涵界定
关于耕地的定义有很多种，但立足于多学科理论，耕地资源同时隶属于多个系统的子系统——土地资源、社会经济、生态环境，而其本身又可视作由自然、社会、经济、生态环境等要素等构成，所有要素之间相互关联、相互影响而形成的复杂系统[5-9]。因此本研究认为，耕地资源系统是指能够种植农作物，在一定时空范围内，以实现耕地资源保护为目标，伴随人类的开发利用与影响活动而不断进行物质、能量、信息的交流与交换过程，由自然、经济、社会、生态等要素所构成的复合系统的综合体。
将系统发展状态定量化引入有序度的概念。有序度是度量系统各个要素之间的和谐与配合程度，是对系统发展状态的量化表征。系统的有序度高，说明系统各个要素之间协调、有序，系统的结构和功能得以维持；系统的有序度低，则表明系统自组织混乱、无序，系统的结构和功能遭到破坏，最终导致系统“代谢”的恶性循环。
3耕地资源系统有序度评价模型构建及评价
[BT2 *2]3.1耕地资源系统有序度评价模型构建
耕地资源系统是综合了自然、社会经济、生态环境等要素的复合系统，这些要素较多，难以一一描述。为此，引入序参量的概念。序参量是协同学理论中的重要概念，是指主导系统有序结构的几个变化缓慢的少数变量，同时又具有支配系统发展和演变的重要作用，也称为慢弛豫参量[10]。因此，科学合理地选取耕地资源系统的关键序参量以及对参序量的建模分析，成为系统有序性测度的关键内容。
承接前面的分析，以及对大量文献的总结与归纳，本研究借鉴孟庆松等的有序度模型[11]，对其中部分公式进行适当修正，构建以下耕地资源系统有序性测度模型。
3.2沧州市耕地资源系统有序度评价
3.2.1耕地资源系统序参量的选取
为了评价耕地资源系统的有序度，选取和辨识系统的关键序参量成为有序度评价模型构建的重要内容。耕地资源是以人为核心，由社会-经济-生态环境子系统构成的复合系统，因此，耕地资源系统序参量的选取可分别从社会、经济、生态环境3个子系统来把握。遵循整体性、合理性、可操作性的基本原则，结合研究区实际情况，立足现有数据资料进行序参量的拟定。本研究从社会稳定、经济发展、生态环境安全的目标出发，构建了沧州市耕地资源系统有序度评价的序参量体系（表1）。
针对表1选取的序参量，作几点说明：人口密度、农村人口文盲率是对社会子系统的有序度起负效应的序参量；农业在生产总值占GDP比重越大，经济发展越滞后，因此该指标对经济子系统起负效应；化肥负荷指数、农药负荷指数、农膜负荷指数、盐碱化率是对生态环境子系统起负效应的序参量，而年均降水量对耕地生态环境子系统的效应随其值的变化而 3.2.3沧州市耕地资源系统有序度评价分析
通过模型计算，可分别得到1994—2013年沧州市耕地资源系统及3大子系统的有序度值，现分别从各子系统及综合系统加以评述与说明。
3.2.3.1耕地资源系统有序度整体分析
如图1所示，耕地资源系统有序度与社会-经济-生态环境子系统有序度存在密切的联系。耕地资源系统有序度是3大子系统的综合反映，经济、生态环境子系统有序度变化幅度较大，各子系统之间存在交叉相对变化的趋势，相对来说社会子系统有序度变化不明显。在1994—2013年20年间，沧州市耕地资源系统有序度值由0.38上升至0.65，年均上升1.35%，耕地资源系统有序度等级由最初的基本有序发展到今天的较有序，发展趋势总体良好。大体可分为2个阶段，从1994—2007年，系统整体有序度变化不明显，系统有序度等级一直维持在基本有序状态，究其原因发现，虽然经济子系统有序度一直稳步上升，但由于耕地生态环境子系统有序度处于下降的趋势，生态环境恶化明显，拉低了系统有序度综合水平；从2008—2013年，系统整体发展状态呈现出较明显地上升趋势，有序度值从0.47增长到0.65，5年平均年增长3.6%，系统等级由基本有序到较有序的转变，此期间耕地资源生态环境子系统的有序度也呈现上升趋势。
3.2.3.2耕地资源社会子系统有序度分析
沧州市地处环渤海和京津冀都市圈重要区位，是国务院确定的经济开发区重点城市，近几年，社会发展水平显著提高，集中表现在城镇化发展水平明显改观，已由1994年的13.23%提高到2013年的45.18%；人民生活質量显著提高，已基本实现小康水平；农村人口文化素质也基本实现了无文盲状态，社会科学技术水平显著进步，基础设施建设也在进一步加强。由图1可知，耕地资源社会子系统有序度值由0.35上升到0.70，年均增长1.75%，系统等级由基本有序上升到较有序状态，说明社会子系统有序度对维持耕地资源系统有序度作出了重要的贡献。
3.2.3.2耕地资源经济子系统有序度分析
沧州市20年耕地资源经济子系统保持着平稳快速的增长趋势，总体发展态势良好，有序度值已由最初的0.00上升到2013年的 0.89，其有序度等级呈现由无序到基本有序再到较有序最后到高度有序的变化历程，2012年达到经济有序度的顶峰0.95，2013年表现出向下降的趋势，原因可能是2013年农业科服人员数量显著减少致使经济子系统有序度水平降低。由此说明，应该加大对农业科技的重视程度，增加科服技术人员的投入水平尤为重要。
3.2.3.3耕地资源生态环境子系统有序度分析
随着近几年耕地集约化水平的不断提高，耕地投入水平与利用强度不断加大，但由于农业科技相对落后，农民耕地生态环保意识相对淡薄等，造成耕地资源出现土壤污染、盐碱化、水土流失等恶化现象。由图1可以看出，耕地资源生态环境子系统变化情况较为复杂，子系统有序度值呈现出上升—下降—上升—平缓不变—下降—上升的态势，其系统等级由高度有序—较有序—基本有序—无序—基本有序的趋势变化，说明了生态环境子系统在20年间发生了螺旋上升的“S”形发展态势。1994—2002年，生态环境子系统有序度发生了曲折变化，但变化幅度较小；2003—2007年，生态环境子系统有序度值明显下降，究其原因发现，4年间耕地集约利用水平显著提高，为了追求更高的产出水平，化肥、农药、农膜投入量越来越大，利用程度也越来越高，农民只知种地，而不注重养地，导致耕地生态环境明显恶化；2008年之后，在生态文明的影响下，人们逐渐意识到耕地生态问题的严重性，从而增强了耕地生态保护意识，林地覆盖率进一步加大，进而增强了耕地资源生态环境系统的生物多样性和稳定性，加大了土地整理复垦工作，盐碱地治理水平进一步提高。
4耕地资源系统交互胁迫关系分析与验证
4.1耕地资源经济-社会-生态环境交互胁迫关系分析
耕地资源社会-经济-生态环境3个子系统之间存在复杂的交互胁迫关系[12]。
4.1.1社会子系统与生态环境子系统的交互胁迫关系
不同的社会发展阶段，2个子系统呈现出不同的胁迫关系。在社会发展的早期，人口的快速增长导致就业和生活需求压力增大，刺激了经济子系统的工业化进程和农业科技的进步，但随着社会城镇化的推进和人口基数的增长，就业压力进一步增大，导致劳动力剩余，某种程度上造成了社会的不稳定并阻碍正常的经济活动；其次，过多的人口增加了国民经济的生活消费负担，进一步制约了经济产业结构的资金、技术密集型产业的高级化发展，严重影响了经济的发展水平[13]。反过来，经济子系统的发展也促使了社会子系统的进步，具体表现在人口素质的提高、生活质量的改善、科学技术的进步，但不可避免产生的社会老龄化问题将再次阻碍经济的发展。
4.1.2社会子系统与生态环境子系统的交互胁迫关系
随着社会的发展和人口压力的增大，生活消费水平日益提高，社会子系统对资源环境的负担也越来越重。就耕地资源来说，在一定时期和一定区域内，耕地是不可再生资源，数量有限，与此同时，随着社会的发展，城镇化进程的加速，建设用地不断占用耕地，导致耕地数量锐减。在外部化利益的诱惑下，人们所实施的增减挂钩政策也总是以次充好补充耕地，导致耕地资源的生态环境质量退化。若没有及时保护与建设耕地，耕地资源生态子系统将不能维持自身结构和功能的稳定，从而导致耕地资源系统的无序和混乱状态；与此相反，若社会保护耕地资源的意识加大并予以落实，积极改进耕作生态科技水平及不断促进生态管理水平，则能够改善耕地生态环境，实现社会的生态文明与精神文明的高层次目标，从而达到社会进步与生态环境和谐发展的良性循环。
4.1.3经济子系统与生态环境子系统的交互胁迫关系
随着经济水平的提高和科技水平的进步，为了追求更大的经济效益，耕地利用越来越集约化，化肥、农药等投入强度越来越大，造成了耕地资源土壤污染、土地板结、耕地生态环境进一步恶化，而弱化的耕地生态环境子系统将反过来制约经济子系统的发展，于是形成惡性循环。由此可见，耕地资源生态环境子系统是保证经济子系统良好发展的前提和基础，而经济子系统也可为耕地生态环境保护提供良好的物质保障，若进行合理的耕地生态保护与建设，将促进耕地资源系统的良性循环。 由此可见，在耕地资源复合系统中，各子系统之间既有恶性制约也有良性循环，并且处于动态变化的过程之中。在某一个发展阶段，系统之间胁迫作用强烈，整体表现出恶性循环；当某一个子系统得到积极的调控措施发挥作用时，2个系统之间的胁迫作用得以减缓，新一轮的恶性、良性循环又将开始。系统在相互制约、相互作用下，表现出“S”形演变规律[14]。
4.2耕地资源系统相互胁迫关系验证
众多专家学者已经证实并得出社会发展水平-经济发展水平-生态环境水平两两之间存在的交互作用曲线关系[15-18]。本研究在耕地资源社会-经济-生态环境子系统两两之间发生交互胁迫关系时，引入上述的3种关系模型进行3个子系统两两拟合分析，通过曲线拟合，不仅能够将一些离散分布的点找出一定的规律性，还能进行预测分析。但须要作如下定义：耕地资源生态环境子系统的无序度指的是耕地生态环境的恶化指数。为此，需要对耕地资源生态环境子系统有序度评价的序参量进行变向处理，即对生态环境子系统的有序性产生正效应的序参量对耕地资源子系统的无序性产生负效应，反之亦然。依次可计算耕地资源生态环境子系统的无序度值。具体的3个曲线模型公式如下。
耕地资源社会-经济子系统存在对数曲线关系，表示为：
[JZ（]y=algx b。[JZ）][JY]（6）
式中：y为耕地资源社会子系统的有序度值；x为耕地资源经济子系统的有序度值；a、b为待定系数。
耕地资源经济-生态环境子系统存在倒“U”形曲线关系，表示为：
可知，拟合优度可以反映曲线拟合的水平，由SPSS处理后，沧州市各县（市）经济-社会对数曲线拟合度处于中等水平，相比较而言，黄骅市、肃宁县拟合度水平较高，沧县、南皮县、任丘市等拟合水平较差，河间市、青县拟合水平处于中等水平。由图2可以看出，除泊头市外，沧州市各县（市）耕地资源经济-社会子系统表现出明显的协同发展的关系，由曲线的拟合形状来看，随着经济的发展，肃宁县、黄骅市社会子系统的变化趋势较快，而吴桥县、献县社会子系统发展水平较缓慢。
4.2.2耕地资源社会子系统与生态环境子系统的交互胁迫关系验证
根据公式z=m-n[10（y-b）/a-p]2倒“U”形曲线，应用Matlab 7.0软件对1994—2013年沧州市及各县（市）社会-生态子系统发展水平进行双指数曲线拟合，拟合函数关系见表4、图3。
由表4、图3可看出，除青县、孟村回族自治县、泊头市、河间市外，其他各县市20年来耕地社会-生态环境子系统的双指数拟合水平较理想，其中，沧县、海兴县、肃宁县、东光县等表现出明显的倒“U”形的特征。
由表4双指数函数的几个参系数可知，m值代表曲线达到拐点时生态环境子系统的恶化程度。由图3可以看出，沧县、东光县、海兴县等生态环境子系统的恶化程度已达到拐点，生态环境发展态势已出现好转的迹象，然而，河间市的耕地生态环境子系统仍有继续恶化的趋势，所以须要对当地的耕地生态加以保护，防止耕地生态环境继续恶化下去。n值表示耕地生态环境子系统有序度发展变化的速率大小，值越大，说明耕地生态环境子系统有序度变化速率越快。由n值来看，南皮县>沧县>河间市>……>青县>任丘市，表明随着社会子系统的变化，南皮县耕地生态环境的恶化速度较快。b值表示曲线出现最高点的位置，其值越大，表明出现最高点时社会子系统的有序度发展水平越高。由此可知，孟村回族自治县、任丘市生态环境恶化程度出现最高点比较晚，而相对应的耕地社会子系统发展也较快；而盐山县、黄骅市b值最小，表明其生态环境子系统恶化程度最高點出现比较早，生态环境从此得到了有效控制。
4.2.3耕地资源经济子系统与生态环境子系统的交互胁迫关系验证
根据方程z=m-n（x-p）2倒“U”形曲线，应用SPSS 17.0软件对1994—2013年沧州市及各县（市）经济-生态子系统发展程度进行二项式曲线拟合（表5、图4）。
可看出，除海兴县、沧县等外，沧州市其他各县市1994—2013年耕地资源经济-生态环境子系统的二项式曲线拟合水平较理想，其中，青县、盐山县、南皮县、任丘市等表现出明显的倒“U”形相互胁迫关系。m值反映了耕地生态环境子系统的承载力阈值，由m值来看，沧县>海兴县>献县>孟村县>……>泊头市>河间市，由此可知，沧县的耕地生态环境子系统承载力水平较高，而河间市较差。另外，曲线存在一个最高点，由此被分为2个部分，在最高点以左，随着耕地资源经济子系统有序度值的提高，耕地生态环境恶化程度加重；在最高点以右，生态环境子系统无序度值随着经济子系统的发展而降低，即耕地生态环境出现好转的趋势。黄骅市经过20年的发展，耕地生态环境子系统刚达到拐点，因此预测今后几年黄骅市生态环境子系统有序度提高并且经济保持稳步增长水平。而从东光县、海兴县、肃宁县、河间市等几个县市的拟合曲线可以看出，从1994—2013年，随着经济子系统有序度水平的提高，生态环境子系统的有序度有进一步恶化的趋势，这就要求耕地利用应因地制宜，采取合理的耕地生态保护措施，遏制耕地生态环境继续恶化。
5结论与讨论
5.1结论
本研究对沧州市1994—2013年耕地资源系统有序度进行评价，研究表明，耕地资源系统有序度是3个子系统的综合反映，经济子系统与生态环境子系统有序度变化幅度较大，各子系统之间存在交叉相对变化的趋势，相对来说社会子系统有序度变化不明显。具体来说，在1994—2013年的20年间，沧州市耕地资源系统有序度值由0.38上升至0.65，年均上升1.35%，耕地资源系统有序度等级由最初的基本有序发展到至今的较有序，发展趋势总体良好。
通過对耕地经济-社会子系统间交互胁迫关系与模型验证分析，研究表明，在1994—2013年间，除泊头市外，沧州市各[CM（25]县（市）耕地资源经济-社会子系统表现出明显的协同发 展关系，随着经济的发展，肃宁县、黄骅市社会子系统的变化趋势较快，而吴桥县、献县社会子系统发展水平较缓慢。通过耕地社会-生态环境子系统的交互胁迫关系验证表明，社会的发展与生态环境子系统的胁迫关系较明显，其中南皮县耕地生态环境的恶化速度较快，孟村县、任丘市生态环境恶化程度出现最高点比较晚，而相对应的耕地社会子系统发展也较快；而盐山县、黄骅市的生态环境子系统恶化程度最高点出现比较早，生态环境从此得到了有效控制。通过耕地经济-生态环境子系统的交互胁迫关系验证表明，沧县的耕地生态环境子系统承载力水平较高，而河间市较差，而东光县、海兴县、肃宁县、河间市等耕地生态环境子系统有进一步恶化的趋势，这就要求耕地利用应因地制宜，采取合理的耕地生态保护措施，遏制耕地生态环境继续恶化。
5.2讨论
本研究立足于系统论和协同学理论，应用计量经济分析方法将耕地资源看作由经济、社会、生态环境3大子系统相互关联、相互影响、相互制约所构成的复合系统，从三维视角研究耕地资源系统在一定时空范围内，经济-社会-生态环境子系统间的协调发展机制、相互胁迫关系与变化规律。由于数据资料获取比较困难，对于耕地资源系统有序度评价序参量的选择存在一定的局限性；但该研究涉及生态学、环境学、经济学、社会学等多个学科，为自然和社会学科间的交叉研究提供了新理念，并为耕地资源可持续发展研究提出了新的思路。
参考文献：
[1]王培东，赵丽，杨皓，等. 基于PSR模型的沧州市农用地集约利用评价研究[J]. 中国农业资源与区划，2013，34（4）：57-63.
[2]Alcamo J，Endqan M B，Kaspar F，et al. The glass model：a strategy for quantifying global environmental security[J]. Environmental Science and Policy，2001，4（1）：1-12.
[3]杨京平. 生态安全的系统分析[M]. 北京：化学工业出版社，2002.
[4]陈印军，王晋臣，肖碧林，等. 我国耕地质量变化态势分析[J]. 中国农业资源与区划，2011，32（2）：1-5.
[5]唐建. 耕地价值评价实证研究[M]. 北京：中国农业出版社，2011.
[6]陈宏民. 系统工程导论[M]. 北京：高等教育出版社，2006.
[7]宇振荣. 土地利用系统分析方法及实践[M]. 北京：中国农业科技出版社，1998.
[8]朱曉华，蔡运龙，蒋为國，等. 自组织的土地利用系统[J]. 中国土地科学，2005，19（2）：21-28.
[9]黄成毅. 成都平原耕地资源系统协同性分析与调控研究[D]. 雅安：四川农业大学，2011.[ZK）]
[10]龙妍. 基于物质流、能量流与信息流协同的大系统研究[D]. 武汉：华中科技大学，2009.
[11]孟庆松，韩文秀. 复合系统整体协调度模型研究[J]. 河北师范大学学报（自然科学版），1999，23（2）：177-179，187.
[12]高乐华. 我国海洋生态经济系统协调发展测度与优化机制研究[D]. 中国海洋大学，2012.
[13]左学金. 人口增长对经济发展的影响[J]. 国际经济评论，2010（6）：127-135.
[14]方创琳. 河西走廊：绿洲支撑着城市化——与仲伟志先生商榷[J]. 中国沙漠，2003，23（3）：128-130.
[15]周一星. 城市地理学[M]. 北京：商务印书馆，1995.
[16]许学强，周一星，宁越敏. 城市地理学[M]. 北京：高等教育出版社，1997.
[17]Grossman G M，Krueger A B. Economic growth and the environment[J]. Quarterly Journal of Economics，1995，110（2）：353-377.
[18]乔标，方创琳，黄金川. 干旱区城市化与生态环境交互耦合的规律性及其验证[J]. 生态学报，2006，26（7）：2183-2190.
关键词：耕地资源系统；系统有序度；交互胁迫；沧州市；可持续发展
中图分类号： F301.24文献标志码： A
文章编号：1002-1302（2017）05-0256-09
耕地作为宝贵资源，是由自然、经济、社会、生态、环境等各个子系统相互关联、相互制约所构成的复合系统的综合体。在某一耕地，经济、社会、生态环境子系统之间能够相互配合、相互促进、协调发展，耕地资源系统就会实现健康、有序、可持续发展的良性循环；反之，耕地资源系统就会陷入生态环境恶化、社会经济紊乱的恶性循环。当前，随着城镇化步伐的加快，建設用地不断占用耕地，尽管我国耕地保护制度较为严格，但耕地数量锐减、耕地质量下降的问题依然存在[1]；与此同时，面对人地关系紧张与粮食安全问题的严峻形势，为满足粮食需求和获得更大的产出效益，耕地利用集约化显著，大量的化肥、农药投入在一定程度上产生了一系列耕地生态安全问题[2-4]。
面对耕地生态保护与社会经济发展之间的博弈问题，将耕地資源作为经济-社会-生态环境的复合系统，探讨各子系统之间的相互胁迫关系与作用机制，探寻耕地资源系统经济-社会-生态环境的协调发展途径，对于实现研究区乃至其他地区耕地资源系统的可持续发展具有重要意义；同时，耕地资源系统的相关研究也将作为土地利用修编的重要命题之一，为政府部门制定和改进相关政策建议提供一定的理论依据。
1研究区概况及数据来源
1.1研究区概况
沧州市位于河北省东南部，环渤海湾经济区，地处 116°10′～117°44′E，37°53′～38°51′N，水文地质条件较差，中低产田和盐碱地分布较广，是河北省粮、棉、油集中产区之一。截至2013年末，全市常住人口为730.95万人，耕地面积为 78.69×104 hm2，沧州市耕地质量总体水平不高，且部分地区生态环境脆弱，土壤盐碱化问题突出。
1.2基础数据来源
本研究数据资料包括耕地数据和社会经济数据2个方面，耕地数据主要来自1994—2013年《河北省土地调查统计年鉴》，社会经济数据来源于1994—2013年《河北农村统计年鉴》及沧州市各县（市）国民经济统计年鉴。
2耕地资源系统有序度的内涵界定
关于耕地的定义有很多种，但立足于多学科理论，耕地资源同时隶属于多个系统的子系统——土地资源、社会经济、生态环境，而其本身又可视作由自然、社会、经济、生态环境等要素等构成，所有要素之间相互关联、相互影响而形成的复杂系统[5-9]。因此本研究认为，耕地资源系统是指能够种植农作物，在一定时空范围内，以实现耕地资源保护为目标，伴随人类的开发利用与影响活动而不断进行物质、能量、信息的交流与交换过程，由自然、经济、社会、生态等要素所构成的复合系统的综合体。
将系统发展状态定量化引入有序度的概念。有序度是度量系统各个要素之间的和谐与配合程度，是对系统发展状态的量化表征。系统的有序度高，说明系统各个要素之间协调、有序，系统的结构和功能得以维持；系统的有序度低，则表明系统自组织混乱、无序，系统的结构和功能遭到破坏，最终导致系统“代谢”的恶性循环。
3耕地资源系统有序度评价模型构建及评价
[BT2 *2]3.1耕地资源系统有序度评价模型构建
耕地资源系统是综合了自然、社会经济、生态环境等要素的复合系统，这些要素较多，难以一一描述。为此，引入序参量的概念。序参量是协同学理论中的重要概念，是指主导系统有序结构的几个变化缓慢的少数变量，同时又具有支配系统发展和演变的重要作用，也称为慢弛豫参量[10]。因此，科学合理地选取耕地资源系统的关键序参量以及对参序量的建模分析，成为系统有序性测度的关键内容。
承接前面的分析，以及对大量文献的总结与归纳，本研究借鉴孟庆松等的有序度模型[11]，对其中部分公式进行适当修正，构建以下耕地资源系统有序性测度模型。
3.2沧州市耕地资源系统有序度评价
3.2.1耕地资源系统序参量的选取
为了评价耕地资源系统的有序度，选取和辨识系统的关键序参量成为有序度评价模型构建的重要内容。耕地资源是以人为核心，由社会-经济-生态环境子系统构成的复合系统，因此，耕地资源系统序参量的选取可分别从社会、经济、生态环境3个子系统来把握。遵循整体性、合理性、可操作性的基本原则，结合研究区实际情况，立足现有数据资料进行序参量的拟定。本研究从社会稳定、经济发展、生态环境安全的目标出发，构建了沧州市耕地资源系统有序度评价的序参量体系（表1）。
针对表1选取的序参量，作几点说明：人口密度、农村人口文盲率是对社会子系统的有序度起负效应的序参量；农业在生产总值占GDP比重越大，经济发展越滞后，因此该指标对经济子系统起负效应；化肥负荷指数、农药负荷指数、农膜负荷指数、盐碱化率是对生态环境子系统起负效应的序参量，而年均降水量对耕地生态环境子系统的效应随其值的变化而 3.2.3沧州市耕地资源系统有序度评价分析
通过模型计算，可分别得到1994—2013年沧州市耕地资源系统及3大子系统的有序度值，现分别从各子系统及综合系统加以评述与说明。
3.2.3.1耕地资源系统有序度整体分析
如图1所示，耕地资源系统有序度与社会-经济-生态环境子系统有序度存在密切的联系。耕地资源系统有序度是3大子系统的综合反映，经济、生态环境子系统有序度变化幅度较大，各子系统之间存在交叉相对变化的趋势，相对来说社会子系统有序度变化不明显。在1994—2013年20年间，沧州市耕地资源系统有序度值由0.38上升至0.65，年均上升1.35%，耕地资源系统有序度等级由最初的基本有序发展到今天的较有序，发展趋势总体良好。大体可分为2个阶段，从1994—2007年，系统整体有序度变化不明显，系统有序度等级一直维持在基本有序状态，究其原因发现，虽然经济子系统有序度一直稳步上升，但由于耕地生态环境子系统有序度处于下降的趋势，生态环境恶化明显，拉低了系统有序度综合水平；从2008—2013年，系统整体发展状态呈现出较明显地上升趋势，有序度值从0.47增长到0.65，5年平均年增长3.6%，系统等级由基本有序到较有序的转变，此期间耕地资源生态环境子系统的有序度也呈现上升趋势。
3.2.3.2耕地资源社会子系统有序度分析
沧州市地处环渤海和京津冀都市圈重要区位，是国务院确定的经济开发区重点城市，近几年，社会发展水平显著提高，集中表现在城镇化发展水平明显改观，已由1994年的13.23%提高到2013年的45.18%；人民生活質量显著提高，已基本实现小康水平；农村人口文化素质也基本实现了无文盲状态，社会科学技术水平显著进步，基础设施建设也在进一步加强。由图1可知，耕地资源社会子系统有序度值由0.35上升到0.70，年均增长1.75%，系统等级由基本有序上升到较有序状态，说明社会子系统有序度对维持耕地资源系统有序度作出了重要的贡献。
3.2.3.2耕地资源经济子系统有序度分析
沧州市20年耕地资源经济子系统保持着平稳快速的增长趋势，总体发展态势良好，有序度值已由最初的0.00上升到2013年的 0.89，其有序度等级呈现由无序到基本有序再到较有序最后到高度有序的变化历程，2012年达到经济有序度的顶峰0.95，2013年表现出向下降的趋势，原因可能是2013年农业科服人员数量显著减少致使经济子系统有序度水平降低。由此说明，应该加大对农业科技的重视程度，增加科服技术人员的投入水平尤为重要。
3.2.3.3耕地资源生态环境子系统有序度分析
随着近几年耕地集约化水平的不断提高，耕地投入水平与利用强度不断加大，但由于农业科技相对落后，农民耕地生态环保意识相对淡薄等，造成耕地资源出现土壤污染、盐碱化、水土流失等恶化现象。由图1可以看出，耕地资源生态环境子系统变化情况较为复杂，子系统有序度值呈现出上升—下降—上升—平缓不变—下降—上升的态势，其系统等级由高度有序—较有序—基本有序—无序—基本有序的趋势变化，说明了生态环境子系统在20年间发生了螺旋上升的“S”形发展态势。1994—2002年，生态环境子系统有序度发生了曲折变化，但变化幅度较小；2003—2007年，生态环境子系统有序度值明显下降，究其原因发现，4年间耕地集约利用水平显著提高，为了追求更高的产出水平，化肥、农药、农膜投入量越来越大，利用程度也越来越高，农民只知种地，而不注重养地，导致耕地生态环境明显恶化；2008年之后，在生态文明的影响下，人们逐渐意识到耕地生态问题的严重性，从而增强了耕地生态保护意识，林地覆盖率进一步加大，进而增强了耕地资源生态环境系统的生物多样性和稳定性，加大了土地整理复垦工作，盐碱地治理水平进一步提高。
4耕地资源系统交互胁迫关系分析与验证
4.1耕地资源经济-社会-生态环境交互胁迫关系分析
耕地资源社会-经济-生态环境3个子系统之间存在复杂的交互胁迫关系[12]。
4.1.1社会子系统与生态环境子系统的交互胁迫关系
不同的社会发展阶段，2个子系统呈现出不同的胁迫关系。在社会发展的早期，人口的快速增长导致就业和生活需求压力增大，刺激了经济子系统的工业化进程和农业科技的进步，但随着社会城镇化的推进和人口基数的增长，就业压力进一步增大，导致劳动力剩余，某种程度上造成了社会的不稳定并阻碍正常的经济活动；其次，过多的人口增加了国民经济的生活消费负担，进一步制约了经济产业结构的资金、技术密集型产业的高级化发展，严重影响了经济的发展水平[13]。反过来，经济子系统的发展也促使了社会子系统的进步，具体表现在人口素质的提高、生活质量的改善、科学技术的进步，但不可避免产生的社会老龄化问题将再次阻碍经济的发展。
4.1.2社会子系统与生态环境子系统的交互胁迫关系
随着社会的发展和人口压力的增大，生活消费水平日益提高，社会子系统对资源环境的负担也越来越重。就耕地资源来说，在一定时期和一定区域内，耕地是不可再生资源，数量有限，与此同时，随着社会的发展，城镇化进程的加速，建设用地不断占用耕地，导致耕地数量锐减。在外部化利益的诱惑下，人们所实施的增减挂钩政策也总是以次充好补充耕地，导致耕地资源的生态环境质量退化。若没有及时保护与建设耕地，耕地资源生态子系统将不能维持自身结构和功能的稳定，从而导致耕地资源系统的无序和混乱状态；与此相反，若社会保护耕地资源的意识加大并予以落实，积极改进耕作生态科技水平及不断促进生态管理水平，则能够改善耕地生态环境，实现社会的生态文明与精神文明的高层次目标，从而达到社会进步与生态环境和谐发展的良性循环。
4.1.3经济子系统与生态环境子系统的交互胁迫关系
随着经济水平的提高和科技水平的进步，为了追求更大的经济效益，耕地利用越来越集约化，化肥、农药等投入强度越来越大，造成了耕地资源土壤污染、土地板结、耕地生态环境进一步恶化，而弱化的耕地生态环境子系统将反过来制约经济子系统的发展，于是形成惡性循环。由此可见，耕地资源生态环境子系统是保证经济子系统良好发展的前提和基础，而经济子系统也可为耕地生态环境保护提供良好的物质保障，若进行合理的耕地生态保护与建设，将促进耕地资源系统的良性循环。 由此可见，在耕地资源复合系统中，各子系统之间既有恶性制约也有良性循环，并且处于动态变化的过程之中。在某一个发展阶段，系统之间胁迫作用强烈，整体表现出恶性循环；当某一个子系统得到积极的调控措施发挥作用时，2个系统之间的胁迫作用得以减缓，新一轮的恶性、良性循环又将开始。系统在相互制约、相互作用下，表现出“S”形演变规律[14]。
4.2耕地资源系统相互胁迫关系验证
众多专家学者已经证实并得出社会发展水平-经济发展水平-生态环境水平两两之间存在的交互作用曲线关系[15-18]。本研究在耕地资源社会-经济-生态环境子系统两两之间发生交互胁迫关系时，引入上述的3种关系模型进行3个子系统两两拟合分析，通过曲线拟合，不仅能够将一些离散分布的点找出一定的规律性，还能进行预测分析。但须要作如下定义：耕地资源生态环境子系统的无序度指的是耕地生态环境的恶化指数。为此，需要对耕地资源生态环境子系统有序度评价的序参量进行变向处理，即对生态环境子系统的有序性产生正效应的序参量对耕地资源子系统的无序性产生负效应，反之亦然。依次可计算耕地资源生态环境子系统的无序度值。具体的3个曲线模型公式如下。
耕地资源社会-经济子系统存在对数曲线关系，表示为：
[JZ（]y=algx b。[JZ）][JY]（6）
式中：y为耕地资源社会子系统的有序度值；x为耕地资源经济子系统的有序度值；a、b为待定系数。
耕地资源经济-生态环境子系统存在倒“U”形曲线关系，表示为：
可知，拟合优度可以反映曲线拟合的水平，由SPSS处理后，沧州市各县（市）经济-社会对数曲线拟合度处于中等水平，相比较而言，黄骅市、肃宁县拟合度水平较高，沧县、南皮县、任丘市等拟合水平较差，河间市、青县拟合水平处于中等水平。由图2可以看出，除泊头市外，沧州市各县（市）耕地资源经济-社会子系统表现出明显的协同发展的关系，由曲线的拟合形状来看，随着经济的发展，肃宁县、黄骅市社会子系统的变化趋势较快，而吴桥县、献县社会子系统发展水平较缓慢。
4.2.2耕地资源社会子系统与生态环境子系统的交互胁迫关系验证
根据公式z=m-n[10（y-b）/a-p]2倒“U”形曲线，应用Matlab 7.0软件对1994—2013年沧州市及各县（市）社会-生态子系统发展水平进行双指数曲线拟合，拟合函数关系见表4、图3。
由表4、图3可看出，除青县、孟村回族自治县、泊头市、河间市外，其他各县市20年来耕地社会-生态环境子系统的双指数拟合水平较理想，其中，沧县、海兴县、肃宁县、东光县等表现出明显的倒“U”形的特征。
由表4双指数函数的几个参系数可知，m值代表曲线达到拐点时生态环境子系统的恶化程度。由图3可以看出，沧县、东光县、海兴县等生态环境子系统的恶化程度已达到拐点，生态环境发展态势已出现好转的迹象，然而，河间市的耕地生态环境子系统仍有继续恶化的趋势，所以须要对当地的耕地生态加以保护，防止耕地生态环境继续恶化下去。n值表示耕地生态环境子系统有序度发展变化的速率大小，值越大，说明耕地生态环境子系统有序度变化速率越快。由n值来看，南皮县>沧县>河间市>……>青县>任丘市，表明随着社会子系统的变化，南皮县耕地生态环境的恶化速度较快。b值表示曲线出现最高点的位置，其值越大，表明出现最高点时社会子系统的有序度发展水平越高。由此可知，孟村回族自治县、任丘市生态环境恶化程度出现最高点比较晚，而相对应的耕地社会子系统发展也较快；而盐山县、黄骅市b值最小，表明其生态环境子系统恶化程度最高點出现比较早，生态环境从此得到了有效控制。
4.2.3耕地资源经济子系统与生态环境子系统的交互胁迫关系验证
根据方程z=m-n（x-p）2倒“U”形曲线，应用SPSS 17.0软件对1994—2013年沧州市及各县（市）经济-生态子系统发展程度进行二项式曲线拟合（表5、图4）。
可看出，除海兴县、沧县等外，沧州市其他各县市1994—2013年耕地资源经济-生态环境子系统的二项式曲线拟合水平较理想，其中，青县、盐山县、南皮县、任丘市等表现出明显的倒“U”形相互胁迫关系。m值反映了耕地生态环境子系统的承载力阈值，由m值来看，沧县>海兴县>献县>孟村县>……>泊头市>河间市，由此可知，沧县的耕地生态环境子系统承载力水平较高，而河间市较差。另外，曲线存在一个最高点，由此被分为2个部分，在最高点以左，随着耕地资源经济子系统有序度值的提高，耕地生态环境恶化程度加重；在最高点以右，生态环境子系统无序度值随着经济子系统的发展而降低，即耕地生态环境出现好转的趋势。黄骅市经过20年的发展，耕地生态环境子系统刚达到拐点，因此预测今后几年黄骅市生态环境子系统有序度提高并且经济保持稳步增长水平。而从东光县、海兴县、肃宁县、河间市等几个县市的拟合曲线可以看出，从1994—2013年，随着经济子系统有序度水平的提高，生态环境子系统的有序度有进一步恶化的趋势，这就要求耕地利用应因地制宜，采取合理的耕地生态保护措施，遏制耕地生态环境继续恶化。
5结论与讨论
5.1结论
本研究对沧州市1994—2013年耕地资源系统有序度进行评价，研究表明，耕地资源系统有序度是3个子系统的综合反映，经济子系统与生态环境子系统有序度变化幅度较大，各子系统之间存在交叉相对变化的趋势，相对来说社会子系统有序度变化不明显。具体来说，在1994—2013年的20年间，沧州市耕地资源系统有序度值由0.38上升至0.65，年均上升1.35%，耕地资源系统有序度等级由最初的基本有序发展到至今的较有序，发展趋势总体良好。
通過对耕地经济-社会子系统间交互胁迫关系与模型验证分析，研究表明，在1994—2013年间，除泊头市外，沧州市各[CM（25]县（市）耕地资源经济-社会子系统表现出明显的协同发 展关系，随着经济的发展，肃宁县、黄骅市社会子系统的变化趋势较快，而吴桥县、献县社会子系统发展水平较缓慢。通过耕地社会-生态环境子系统的交互胁迫关系验证表明，社会的发展与生态环境子系统的胁迫关系较明显，其中南皮县耕地生态环境的恶化速度较快，孟村县、任丘市生态环境恶化程度出现最高点比较晚，而相对应的耕地社会子系统发展也较快；而盐山县、黄骅市的生态环境子系统恶化程度最高点出现比较早，生态环境从此得到了有效控制。通过耕地经济-生态环境子系统的交互胁迫关系验证表明，沧县的耕地生态环境子系统承载力水平较高，而河间市较差，而东光县、海兴县、肃宁县、河间市等耕地生态环境子系统有进一步恶化的趋势，这就要求耕地利用应因地制宜，采取合理的耕地生态保护措施，遏制耕地生态环境继续恶化。
5.2讨论
本研究立足于系统论和协同学理论，应用计量经济分析方法将耕地资源看作由经济、社会、生态环境3大子系统相互关联、相互影响、相互制约所构成的复合系统，从三维视角研究耕地资源系统在一定时空范围内，经济-社会-生态环境子系统间的协调发展机制、相互胁迫关系与变化规律。由于数据资料获取比较困难，对于耕地资源系统有序度评价序参量的选择存在一定的局限性；但该研究涉及生态学、环境学、经济学、社会学等多个学科，为自然和社会学科间的交叉研究提供了新理念，并为耕地资源可持续发展研究提出了新的思路。
参考文献：
[1]王培东，赵丽，杨皓，等. 基于PSR模型的沧州市农用地集约利用评价研究[J]. 中国农业资源与区划，2013，34（4）：57-63.
[2]Alcamo J，Endqan M B，Kaspar F，et al. The glass model：a strategy for quantifying global environmental security[J]. Environmental Science and Policy，2001，4（1）：1-12.
[3]杨京平. 生态安全的系统分析[M]. 北京：化学工业出版社，2002.
[4]陈印军，王晋臣，肖碧林，等. 我国耕地质量变化态势分析[J]. 中国农业资源与区划，2011，32（2）：1-5.
[5]唐建. 耕地价值评价实证研究[M]. 北京：中国农业出版社，2011.
[6]陈宏民. 系统工程导论[M]. 北京：高等教育出版社，2006.
[7]宇振荣. 土地利用系统分析方法及实践[M]. 北京：中国农业科技出版社，1998.
[8]朱曉华，蔡运龙，蒋为國，等. 自组织的土地利用系统[J]. 中国土地科学，2005，19（2）：21-28.
[9]黄成毅. 成都平原耕地资源系统协同性分析与调控研究[D]. 雅安：四川农业大学，2011.[ZK）]
[10]龙妍. 基于物质流、能量流与信息流协同的大系统研究[D]. 武汉：华中科技大学，2009.
[11]孟庆松，韩文秀. 复合系统整体协调度模型研究[J]. 河北师范大学学报（自然科学版），1999，23（2）：177-179，187.
[12]高乐华. 我国海洋生态经济系统协调发展测度与优化机制研究[D]. 中国海洋大学，2012.
[13]左学金. 人口增长对经济发展的影响[J]. 国际经济评论，2010（6）：127-135.
[14]方创琳. 河西走廊：绿洲支撑着城市化——与仲伟志先生商榷[J]. 中国沙漠，2003，23（3）：128-130.
[15]周一星. 城市地理学[M]. 北京：商务印书馆，1995.
[16]许学强，周一星，宁越敏. 城市地理学[M]. 北京：高等教育出版社，1997.
[17]Grossman G M，Krueger A B. Economic growth and the environment[J]. Quarterly Journal of Economics，1995，110（2）：353-377.
[18]乔标，方创琳，黄金川. 干旱区城市化与生态环境交互耦合的规律性及其验证[J]. 生态学报，2006，26（7）：2183-2190.