（环境科学专业论文）区域生态安全承载力的系统动力学研究——以重庆市云阳县为例.pdf

关系进行了具体的分析，找出系统内部之间复杂的反馈关系，建立了系统动力 学模型的反馈回路图。 第六章在生态安全承载力各子系统反馈机制的基础上，建立了系统流程图 和各子系统的系统动力学方程。 第七章在设定模型基本运行参数后，用系统动力学专用v e n s i m 软件运行模 型，并对模型的结构、历史行为和灵敏度进行检验分析，判断模型的有效性； 进而对系统未来趋势进行仿真运行及分析。 第八章选取四个方面的政策指标对云阳县生态安全承载力系统进行政策实 验研究。 第九章在政策实验基础上，调整运行参数，同时考虑在适当的时候调入粮 食，从而实现对云阳县生态安全承载力模型的优化。 第十章得出结论。云阳县生态安全承载力系统动力学模型能清晰地表征云 阳县生态安全承载力系统各要素间的相互关系，可利用该模型对云阳县的生态 安全承载力发展趋势进行了预测，其结果可作为云阳县制订未来发展规划的科 学依据和参考。 关键词：生态安全生态安全承载力系统动力学模型云阳县 s y s t e md y n a m i c ss t u d yo ne c o s e c u r i t y c a p a c i t yi ny u n y a n gc o u n t y ，c h o n g q i n g m a j o r ：e n v i r o n m e n t a ls c i e n c e g r a d u a t es t u d e n t ：x i u j i ey a n g a d v i s e r ：w e n f e n gl u o i nr e c e n ty e a r s ，t h ee c o l o g i c a lp r o b l e m sa r eb e c o m i n gm o r ea n dm o r es e r i o u s d u et op e o p l e sd e s t r o y , w h i c hh a v eg i v e nd i r e c t l yt h r e a to nt h ee x i s t e n c eo fh u m a n b e i n ga n dt h e n a t i o n a ld e v e l o p m e n t s ot h ee c o l o g i c a lp r o b l e m sa r ee v o l v i n g e c o s e c u r i t yo n e ，w h i c ha l s ob e c o m et h ef o c u sb e t w e e ng l o b a la c a d e m ya n d g o v e r n m e n to fa l lt h ec o u n t r i e s e c o s e c u r i t yc a r r y i n gc a p a c i t y ( e c c 、i sak i n do fn e wm e t h o dt h a te v a l u a t e s e e o s e c u r i t y ，w h i c hh a db e e np u tf o r w a r db yac h i n e s es c h o l a r ，j i n g p i n gy a n g b a s e do nt h ea p p l i e dp r o d u c t i o no fc a r r y i n gc a p a c i t yi ne n v i r o n m e n t a ls c i e n c e ，t h i s p a p e rr e d e f i n e se c o - s e c u r i t yc a r r y i n gc a p a c i t y a tt h es a m et i m e ，u s i n gt h em e t h o d o fs y s t e md y n a m i c s ( s d ) a n dt a k i n gy u n y a n gc o u n t ya sa ne x a m p l e ，t h ep a p e r m a k e sat r i a lo fr e s e a r c ho ne c o s e c u d t yc a r r y i n gc a p a c i t yo fd i s t r i c t t h ea r t i c l ei sd i v i d e di n t ot e nc h a p t e r s t h ew h o l es t r u c t u r ea sf o l l o w s ： t h eb a c k g r o u n d ，s i g n i f i c a n c eo fe c o - s e c u r i t yc a r r y i n gc a p a c i t ys t u d ya n dt h e w o r ko ft h i sa r t i c l ea r ei n t r o d u c e di nc h a p t e r1 t h et h e o r yo fe c o s e c u r i t yc a r r y i n g c a p a c i t yi sd i s c u s s e di nc h a p t e r2 i t c o n s i s t e do ft w op a r t st h a ta r ee c o s e c u r i t ya n d e c o s e c u r i t yc a r r y i n gc a p a c i t y t h e o r i g i n 。c o n c e p t i o n ，c h a r a c t e r i s t i c s ，c o n n o t a t i o n ，c u r r e n ts i t u a t i o na n dr e l a t i o n s b e t w e e nt w op a r t sa r es e p a r a t e l yi n t r o d u c e d i ti si n t r o d u c e dt h em e t h o do fs y s t e md y n a m i c si nc h a p t e r3 a na p p r o a c ht o t h ep r o b l e m so fs y s t e md y n a m i c si sd i v i d e di n t of i v es t e p si n c l u d ei d e n t i f y i n g s y s t e m ，a n a l y s i sc o n f i g u r a t i o n ，e s t a b l i s h i n gm o d e l a n a l y s i sm o d e la n de v a l u a t i o n m o d e l i nc h a p t e r4 ，t h es y s t e mo fe c o s e c u r i t yc a r r y i n gc a p a c i t yi ny u n y a n gc o u n t y i si d e n t i f i e d ，w h i c hc o n s i s t e do ft h ep h y s i c a ls u b s y s t e mt h a ti sc o m p o s e do fn a t u r a l g e o g r a p h y , n a t u r a lr e s o u r c e s ，e c o l o g i c a le n v i r o n m e n t ，t h es o c i a ls u b s y s t e m t h a ti n e l u d e st h ep o p u l a t i o nc o n d i t i o n s ，s c i e n c ea n dt e c h n o l o g y ，t h ee c o n o m i cs u b s y s t e m w h i c hc o m p r i s e sn a t i o n a le c o n o m ya n da g r i c u l t u r e ，i n d u s 时，b u s i n e s s ，e t c o nt h eb a s i so fa n a l y z i n gt h ep o p u l a t i o n ，r e s o n r c e s ，e n v i r o n m e n ta n ds o c i a l e c o n o m i cc o n d i t i o n so fy u n y a n gc o u n t y , t h ec o n s t r u c t i o no fe c o s e c u r i t yc a r r y i n g c a p a c i t ys y s t e mi ny u n y a n gc o u n t yi se s t a b l i s h e d ，t h es y s t e md y n a m i c sm o d e l a n df e e d b a c kl o o p sa r eb u i l ti nc h a p t e r5 b a s e do nt h ef e e d b a c km e c h a n i s mo fs u b s y s t e ma n dt h ep r i n c i p l eo fc a t c h i n g t h em a i np r o b l e m s ，t h ef l o wc h a r to fe c o s e c u r i t yc a r r y i n gc a p a c i t ys y s t e ma n dt h e e q u a t i o n so fs y s t e md y n a m i c sf o rs u b s y s t e ma r es e tu pi nc h a p t e r6 i nc h a p t e r7 ，u s i n gv e n s i ms o f t w a r et h es y s t e md y n a m i c sm o d e lw i l lb e c i r c u l a t e da f t e rt h ep a r a m e t e r so fm o d e li ss e tu p t h r o u g hv e r i f y i n gt h ec o n s t r u c t i o n ， h i s t o r i c a lb e h a v i o r st h a ti fo rn o tf i th i s t o r i c a ld a t aa n dt h ea n a l y s i so fs e n s i t i v i t y ，t h e v a l i d i t ya n dp r a c t i c a b i l i t yo fm o d e li se s t i m a t e d f o l l o w i n g ，e c o s e c u r i t yc a r r y i n g c a p a c i t yi ny u n y a n g i sp r e d i c t e db yt h em o d e l t h ep o l i c ye x p e r i m e n t so fe c o s e c u r i t yc a r r y i n gc a p a c i t ym o d e la r ed i s c u s s e d b yt h ep o l i c ys y n t h e s i si n d e x e so ff o u ra s p e c t si uc h a p t e r8 a c c o r d i n gt ot h ep o l i c ye x p e r i m e n t s ，i nc h a p t e r9 ，a f t e ra d j u s t i n gt h ep o l i c y p a r a m e t e r s ，a tt h es a m et i m e ，t a k i n gp u r c h a s eg r a i ni n t oa c c o u n ta ts u i t a b l et i m e ， e c o s e c u r i t yc a n t i n gc a p a c i t ym o d e li ny u n y a n g i so p t i m i z e d f i n a l l y , t h i sa r t i c l eg i v e so u tt h ec o n c l u s i o n s i tc l e a r l yi n d i c a t e st h a tt h i sm o d e l c a ne x p l a i nt h ei n t e r r e l a t i o n so fe a c hf a c t o ro fe c o s e c u r i t yc a r r y i n gc a p a c i t y s y s t e ma b o u ty u n y a n g ，s oi tc a np r e d i c te c o s e c u r i t yc a r r y i n gc a p a c i t y ，a n dt h e r e s u l t sc a nb ea p p l i e da sas c i e n t i f i cf o u n d a t i o na n dr e f e r e n c ef o rg o v e r n m e n t a l d e v e l o p m e n ta n dp l a n k e y w o r d s ：e c o s e c u r i t y ；e c o - s e c u r i t yc a r r y i n gc a p a c i t y s y s t e md y n a m i c sm o d e l ；y u n y a n gc o u n t y ； 四川大学硕_ i 学位论文 1 导论 1 1 选题背景与意义 1 1 1 选题背景 近年来，人类对生态的破坏越来越严重已直接威胁到人类的生存和国家 发展，生态问题演化为生态安全问题，并引起国际学术界和各国政府的广泛关 注。生态安全研究成为当前地学、资源学和环境科学研究的前沿任务和重要领 域。 曲格平指出，环境污染导致的生态环境恶化，毁林、毁牧引起的土地退化， 严重的水土流失使大片土地变为荒漠，这意味着大片国土将失去对国民经 济的承载能力；生态环境的破坏，会造成工农业生产能力的降低和人民生活水 平的下降，进而可能引起经济衰退、暴力冲突和政治动荡口】。因此，生态安全 是一个国家发展、进步的基础安全问题。 生态安全是实现可持续发展的前提和保障。维护生态安全的目的，就是要 使人类的活动在不突破资源和环境承载力的条件下，提高经济活动效率，改善 人民健康状况和提高生活质量，避免因自然资源衰竭、资源生产率下降、环境 污染和退化给社会生活和生产造成的短期灾害和长期不利影吼，以实现经济社 会的可持续发膨”。 西部大开发是我国2 1 世纪发展中的重大战略之一。脆弱的西部生态环境的 保护和建设是西部大开发的根本，是西部实现可持续发展的关键。也直接关系 到东部地区乃至全国的生态安全。故而，西部地区的生态安全状况，必将对我 国二十一世纪的可持续发展产生深远影响川1 5 j 。 基于以上原因，本文拟用系统动力学( s y s t e m d y n a m i c s ，s d ) 模型对地处 西部地区的重庆市云阳县的生态安全承载力进行实证研究，以期对生态安全的 研究提供新的思路，并对区域经济和社会发展规划、决策的制订提供理论参考。 1 1 2 研究意义 从系统功能层面而言，本研究可以动态仿真和模拟生态安全承载力系统， 模拟多个指标的演进。进行生态安全承载力系统的政策实验和分析。 从理论研究层面而言，基于系统动力学模型研究生态安全承载力系统有一 1 导论 定的新意：而基于系统动力学方法进行生态安全承载力研究，也是生态安全研 究方法的一种新的尝试。本研究在生态安全承载力系统的结构、功能、指标量 化、变量关系分析等方面进行了一定的理论探讨和研究。 实证研究层面而言，基于系统动力学方法研究某个具体区域的生态安全承 载力系统并应用于生态安全管理具有一定现实意义。 1 2 本文的工作 本文的主要工作有以下三个方面： 第一，云阳县历史与现状调查及发展战略研究 详细了解云阳县近几年来经济社会形势以及当前与将来要采取的发展战略 和政策措施，同时了解经济发展与资源间的关系以及对生态安全方面的影响。 在对云阳县人口、资源、环境和社会经济系统辩识的基础上，形成对生态 安全承载力系统的总体认识，为定量研究提供更为详实的前提准备。 第二，数据收集、整理与分析以及系统动力学模型的建立 根据云阳县生态安全承载力系统实际的结构、建模构思收集有关云阳县解 放以来的国民经济、环境保护等方面的资料，并对其进行整理和分析。 探讨获得模型数据的可能性，使数据和模型在研究目的的框架下统一起来。 通过确定系统研究的因果关系，建立起云阳县生态安全承载力的系统动力 学模型。 第三，政策实验和系统优化 运用模型从投资倾斜、调整结构等方面进行系列的政策实验，进而对云 阳县生态安全承载力系统进行优化，为云阳县中、长期的发展战略和政策措施 的制定提供理论参考。 网川大学硕：学位论文 2 生态安全承载力综述 2 1 生态安全概述 2 1 1 生态安全概念 安全是针对于风险来说的。近年来，人类对生态的破坏越来越严重，水土 流失、干旱洪涝、沙尘暴、泥石流、水污染、大气污染等都在直接或间接地威 胁着人类的健康和社会经济的发展，部分地区由于人口激增，人类对自然资源的 开发消耗过度，资源的不当利用活动同益频繁，使得生态环境日趋恶化，已直 接威胁人类的生存。生态风险的凸现使人们对生态问题予以高度重视并提出了 生态安全这种新的安全观【6 】。 早在1 9 7 7 年，美国的著名学者l e s t e rr b r o w n 在“r e d e f i n i n gn a t i o n a l s e c u r i t y ”报告中f _ ”，首次将环境引入安全概念，并力图将环境问题纳入国家安 全，然而未能引起美国政府和国际学术界的重视。直到1 9 8 9 年，国际应用系统 分析研究所才对生态安全的含义作了阐述：生态安全是指在人的生活、健康、 安乐、基本权利、生活保障来源、必要的资源、社会秩序和人类适应环境变化 的能力等方面不受威胁，其中包括自然生态安全、经济生态安全和政治生态安 全 1 】。1 9 9 3 年，n o r m a nm y e r s y 在“u l t i m a t es e c u r i t y _ 1 h ee n v i r o n m e n t a lb a s i s o f p o l i t i c a ls t a b i l i t y ”一书中指出，生态安全涉及由地区间的资源战争和全球的生 态威胁而引起的环境退化，这些问题继而波及经济和政治的不安全；并认为生 态安全是人类的终极安全【i 】。k a t r i n as r o g e r s 认为，生态安全是指人类社会创 造一种条件，使得自然环境在满足人类需求的同时使自然资源存量不致减少【8 】。 d e n n i sp h a g e s 则认为，生态安全是建立在维持人类与自然环境四类平衡人 类需求与环境系统可持续承受力的平衡、人口内部的平衡、人口与其它物种的 平衡以及人口与致病微生物的平衡基础之上的安全【9 】。 尹希成认为，生态安全是指人类赖以生存的生物圈处于自然平衡状态，在 这种状态下，人能够与自然界共生、共荣，协同进化【l 0 】。 郭中伟从生态系统的角度定义生态安全。他认为，生态安全是指一个生态 系统的结构是否受到破坏，其生态功能是否受到损害。生态安全的显性特征是 生态系统所提供的服务的质量或( 和) 数量的状态。当一个生态系统所提供的服 务的质量或数量出现异常时，则表明该系统的生态安全受到了威胁，即处于生 2 生态安全承载力综述 态不安全状态。因此，生态安全包含两重含义：一是生态系统自身是否安全， 即自身结构是否受到破坏；二是生态系统对于人类是否安全、即生态系统所提 供的服务是否满足人类的生存需要f 川。 全国人大环境与资源保护委员会主任曲格平从两个方面解释生态安全：一 是防止生态环境的退化对经济基础构成威胁，主要指环境质量状况和自然资源 的减少、退化，削弱了经济可持续发展的支撑能力；二是防止环境问题引发公 众的不满，特别是导致环境难民的大量产生，从而导致国家的动荡射。 陈国阶认为，生态安全是指人类赖以生存的生态与环境，包括聚落、聚区、 区域、国家乃至全球，不受生态条件、状态极其变化的胁迫、威胁、危害、损 害乃至毁灭，能处于正常的生存和发展的状态【l ”。 肖笃宁将生态安全与保障程度相联系，把生态安全定义为：人类在生产、 生活和健康等方面不受生态破坏与环境污染等影响的保障程度，包括饮用水与 食物安全、空气质量与绿色环境等基本要素【l 引。 显然，从不同角度可以对生态安全做出不同的理解和诠释。但就生态安全 的本质而言，生态安全是围绕人类社会可持续发展的目的，促进经济、社会和 生态三者之间和谐统一，由生物安全、环境安全和生态系统安全组成的安全体 系l l4 。 2 1 - 2 生态安全研究现状 2 1 2 1 国外研究现状 在宏观上国外研究主要围绕生态安全的概念及生态安全与国家安全、民 族问题、军事战略、可持续发展和全球化的相互关系而展开。1 9 9 6 年地球公 约的面对全球生态安全的市民条约中，约有1 0 0 多个国家的2 0 0 多万人 签字。该条约建立在生态安全、可持续发展和生态责任的基础之上，要求各成 员国和各团体组织互相协调利益，履行责任和义务【j l 。 从微观角度看，目前国外关于生态安全的研究主要集中在两个方面：一是 基因工程生物的生态( 环境) 风险与生态( 环境) 安全 二是化学品的施用对农业生 态系统健康及生态( 环境) 安全的影响【1 5 1 。 2 1 2 2 国内研究现状 1 9 9 8 年1 1 月，史培军“人地系统动力学与生态安全建设”一文在科技曰 4 四川大学硕士学位论文 报的发表【l6 1 ，拉开了我国学术界生态安全研究的序幕。广大学者除了对生态安 全的概念进行深入探讨之外，还对生态安全的各个方面进行了广泛研究。张雷 和孟旭光对我国的资源安全问题进行了比较系统而全面地论述【1 7 - 1 9 o 也有些学 者对生物安全、粮食安全、食物安全、水安全的现状与发展作了阐述并提出对 策【2 0 _ 2 4 1 。国家环境保护总局在全国范围开展了生物安全调查，于2 0 0 0 年制定了 国家生物安全框架。结合朱镕基总理在政府工作报告中的讲话，曲格平概 括了生态安全对国家安全的影响，认为生态安全和经济安全、政治安全一样， 也是国家安全的重要基石，指出我国的生态安全已经在国土、水、人民健康和 生物4 个方面突现出来f 2 列。曲格平建议，为应对可能出现的生态环境灾害和社 会事件，需要建立必要的预警、预防和应急措施【2 6 1 。俞孔坚等将生态安全概 念引入生物保护的景观格局研究中【2 7 2 酊。 2 1 ，3 生态安全研究方法“4 1 目前，用于生态安全研究的方法主要有如下几种 2 1 3 1 综合指数法 综合指数法评价的方法过程是： ( 1 ) 分析研究评价的生态安全因子的程度与变化规律； ( 2 ) 建立表征各生态安全因子特性的指标体系； ( 3 ) 确定评价标准： ( 4 ) 建立评价函数曲线，将评价的环境因子的现状值和预测值转化为统 一的无量纲的生态安全与质量指标( 用1 o 表示) ； ( 5 ) 确定生态安全的标准值，计算人类活动对生态安全及环境影响变化 值； ( 6 ) 根据各评价因子的相对重要性赋予权重； ( 7 ) 将各因子的变化值综合，得出综合影响评价值，即 a e = b 。一e 。扛w ； ( 2 1 ) i = l 式中：ae 表示人类活动前后生态安全及质量变化值；凰表示人类活动后i 因子的质量指标；表示人类活动前i 因子的质量标准；w i 表示i 因子的权重。 5 2 生态安全承载力综述 2 1 3 ，2 景观生态学方法 景观生态学方法评价的基本思路是： ( 1 ) 计算景观多样性指数h ，公式为： 抒= 瓴x l n p ，) ( 2 - 2 ) i = 1 式中：p i 为某类型景观所占百分比面积：m 为景观类型数。 ( 2 ) 优势度指数d 的计算，公式为： d = + r ，) 2 + l ，k l o o ( 2 _ 3 ) 式中：r d ( 密度) - ( 拼块i 的数目，拼块总数) 1 0 0 r f ( 频度) = ( 拼块i 出现的样方数，总样方数) 1 0 0 ； l p ( 景观比例) = ( 拼块i 的面积，样地总面积) 1 0 0 ( 3 ) 生态环境质量( 功能与稳定性) e q 的计算： e q = a i n ( 2 4 ) 式中：a i = 当地生态适宜性( 域值0 1 0 0 ) i a 2 = 植被覆盖度( 以土地的实际覆盖度为权值，域值按实际覆盖度降以 1 0 0 ) ； a 3 = 抗退化能力赋值( 群落抗退化能力强时赋值1 0 0 ，较强者赋值6 0 ，一 般水平赋值4 0 。以下赋值o ) ； a 4 = 恢复能力赋值( 群落恢复能力强时赋值8 0 ，较强赋值6 0 一般赋值 4 0 ，以下赋值o ) 。 生态环境质量与生态环境安全级别对照表如表2 - 1 所示。 塞! 二! 皇查墅堕堕曼皇塞查堡丝室全堡型型堕塞 望q 堡! 螋= ! ! ! = ! ! 塑= 垫! ! = ! ! = ! 生态环境安全的级别i i ii i i i v v 2 1 3 3 景现生态安全格局方法 景观生态安全格局的分析识别步骤为： 6 四川大学硕士学位论文 ( 1 ) 源的确定 在大多数情况下，景观生态安全规划的保护对象是多个物种和群体，能充 分反映保护地的多种生境特点，要在区系成分调查基础上，确定作为主要保护 对象的物种和相应的栖息地( 源) 。 ( 2 ) 建立阻力面 物种对景观与生境的利用可以看作是对空间竞争性控制和覆盖过程，而这 种控制和覆盖必须通过克服阻力来实现。因而，阻力面反映了物种空间运动的 趋势。 ( 3 ) 根据阻力面来判别安全格局 阻力面是用等阻力线来表示的一种矢量图，是反映物种运动的时空连续体， 类似地形表面。阻力面在源区下陷，在最不易到达的地区阻力面呈峰突起，而 两陷之间有低阻力的谷线相连，两峰之间有高阻力的脊线相连，每个谷线和脊 线上都各有一鞍，它们是谷线或脊线上的极值。根据阻力面进行空间分析，可 以判别缓冲区、源区联接、辐射道和战略点。 对上述各组分及各种存在的和潜在的景观结构组分叠加、组合，就形成了 某一安全水平上的生物安全保护格局。 2 1 3 4 层次分析法 层次分析法是一种对复杂现象的决策思维过程进行系统化、模型化、数量 化的方法。其具体分析步骤如下： ( 1 ) 明确问题 确定评价范围和评价g t 的、对象，进行生态因子相关分析，明确各因子之 间的相互关系。 ( 2 ) 建立层次结构 将被评价关系按其组成层次构筑成一个树状层次结构，一般地分成为3 个 层次：目标层、指标层、策略层。 ( 3 ) 构造判断矩阵 在每一层次上，按照上一层次的对应准则要求，对该层次的元素( 指标) 进行逐对比较，并用标度1 、3 、5 、7 、9 和2 、4 、6 、8 以及倒数来比较相对重 要性( 表2 2 ) ，并写成矩阵形式，即为判断矩阵。 7 2 生态安全承载虚堡整 ( 4 ) 层次排序计算和一致性检验 即是计算判断矩阵的最大特征值及相应的特征向量。另外，对层次分析所 得结果是否基本合理。还需要对判断矩阵进行一致性检验。其具体计算为： 建立最大特征值方程 a n n w 咄m x w ( 2 5 ) 式中：a 。表示n x n 降互反性判断矩阵；k 。为最大特征值；w 为最大特 征值所对应的特征向量。 计算特征向量w 第一步取归一化的初值向量 断。 w o = ( w o l ，w 0 2 ，w o t i ) ( 2 6 ) 第二步计算a 。w 1 1 = w 直至收敛，则w = ( w l ，w 2 i ，w n i ) a x = a u + ，川 ( i = 1 2 n ) ( 2 - 7 ) j = l 注意：求得的w 便是要求的排序权重，xm a x 可用于矩阵的一磁性检验判 矩阵的一致性判断 判断方程为： c r = ( h n a x n ) ( n - 1 ) r i ( 2 8 ) 网川大学硕上学位论文 当c r 值 o ， r c s lr c cf 时，超载 当0 = o ，lr c si = ir c ci 时，超载 当9 ： lr m c d a k = l t m c d a j + d p ( r m c z r j k r m c j r j k ) 年人工牧草地增量( r c z r ) ： rr c z r k l = g r z r k + d e l a y l ( n x k m a 2 k + n t g m a k ，s t i m e 2 ) 年人工牧草地退化量( r c j r ) ：r r c j r k l = m c t v 2 k + r m c d a k 四川大学硕十学位论文 宜牧荒草地面积( y m c a ) ：ly m c a k - - y m c a j + d t ( n x k m a 2 j k ) 年新开发牧草地( n x k m a 2 ) ： rn x k m a 2 k l = m z xy - k 8 m c kv k y m c a k 天然牧草产量( t c c l ) ： at c c l k = 1 m c v k + t m c a k 人工牧草产量( r c c l ) ： ar c c l k = r m c v k * r m c a k 草地可载畜量( c z x l ) ： ac z x l k = ( t c c l k x c v l + r c c l k ， r x c v 2 ) 羊存栏量( y c l ) ：l y c l k = y c l j + d t * ( n y z r j k n y c r j k ) 羊年增量( n y z r ) ： rn y z r k l = n y z v 4 z x l k * y c l k 羊年出栏( n y c r ) ： rn y c r k l = n y c v k * y c l k 羊肉产量( y r c l ) ：ay r c l k = y r c r * n y c r k 牧草承载力系数( c c z r ) ： ac c z r k = c z ) ( l k z x l k 饲粮承载猪量( l z z l ) ：ac z z l k = s l l k z x l v k 饲粮承载力系数( l c z r ) ： al c z r k ；c z 盈i g z z l k 肉总产量( r z c l ) ： ar z c l k = z r cl k + y r cl k + n r c l k 式中：n x k m a ：年新开牧草地： n t g m a ：年退耕还牧