《论文（设计）基于粮食安全的中国粮食需求预测与耕地资源配置研究30529.doc》由会员分享，可在线阅读，更多相关《论文（设计）基于粮食安全的中国粮食需求预测与耕地资源配置研究30529.doc（8页珍藏版）》请在三一办公上搜索。
1、基于粮食安全的中国粮食需求预测与耕地资源配置研究马永欢 牛文元摘 要：本文以系统动力学原理为基础，仿真模拟了我国2010、2015和2020年的粮食需求。结果表明，随着人口的增加和经济的发展，口粮和种子粮的需求相对稳定，而饲料粮和工业用粮的需求增长明显，粮食总需求呈现刚性增长趋势。在95%和100%两种粮食自给标准下，通过对我国粮食需求与耕地配置的模拟发现，粮食需求与耕地供给呈现出一种紧平衡。确保有效的耕地面积和粮食作物播种面积，实现粮食单产的稳步提高和严格控制人口增长是实现粮食安全的“三重”底线。关键词：粮食安全，粮食需求，耕地配置，调控我国耕地资源的稀缺性日益突出，己成为制约农业生产和国民
2、经济可持续发展的瓶颈因素。随着人口的持续增长和工业化、城镇化的快速发展，我国粮食需求呈现刚性增长，而耕地供给将会表现出不断减少的态势。又由于气候变化、旱涝灾害和水资源短缺等约束因素对粮食生产的影响日益突出，保障粮食安全面临着严峻挑战。对于一个占世界22%的人口大国来说，立足国内为主解决粮食问题是我国一项长期的战略性任务。粮食需求预测是未来耕地配置的重要依据。虽然学术界采用线性方程、指数平滑法和自回归移动平均（ARIMA）模型对口粮、饲料粮、种子粮和工业用粮的时间序列变化进行了分析，对粮食需求预测具有一定的指导作用。然而，粮食生产的自然系统与需求的社会经济系统存在内在联系，只有考虑以工业化为标志
3、的经济发展水平和以城市化为标志的城乡结构对粮食需求的影响，才能从动态角度对粮食需求做出客观判断。对此，本文从社会经济系统相互影响的角度，以系统动力学模型为主模型，以ARIMA模型、计量经济模型和回归分析模型等为辅助模型，对其进行了对接和耦合，建立了粮食需求与耕地优化配置的动态仿真模型。一、研究方法（一）粮食需求的预测方法理顺粮食消费结构是预测粮食需求的先决条件。一般认为，粮食消费需求主要由生活用粮和生产用粮两部分构成。前者包括城乡居民口粮，而后者主要包括饲料粮、工业用粮和种子粮等。虽然国家统计局公布了城镇居民和农村居民平均每天的主要食品消费量，但是没有包括外出消费，对此，城乡居民实际的口粮消费
4、量需要在统计消费量的基础上进行修正。从系统论角度来看，口粮需求与人口数量、城乡结构和城乡消费结构密切相关，而人口的城乡结构与工业化水平一脉相承。本文在此基础上构建了口粮需求的系统动力学模型。城乡居民对肉、蛋、奶和水产品的消费量是间接由粮食转换而来，因此，对饲料粮的预测也采用了系统动力学的原理和方法。工业用粮是粮食需求的重要组成部分。研究表明，工业用粮与工业化发展水平密切相关，因此，对工业用粮的预测采用了回归分析法。而种子粮主要取决于单位耕地面积播种量和播种面积，因此，对种子粮的需求也采用了系统动力学方法。由口粮需求、饲料粮需求、工业用粮需求和种子粮需求，共同构成了粮食总需求模型： 上式中，G表
5、示粮食总需求，G1、G2、G3、G4分别代表口粮、饲料粮、种子粮和工业用粮；pi代表城乡人口，ei代表城乡人均直接口粮消费量；hi代表肉、蛋、奶和水产品消费量，fi代表对应的饲料转换效率，q表示粮食对饲料的转换系数；su表示粮食作物播种面积结构，gu表示对应的单位面积种子投入量；m代表工业消耗的粮食。（二）耕地配置的建模方法从粮食安全的角度来看，一定时期的耕地需求量不但取决于区域的粮食自给率、人均粮食需求量和人口数量，还受耕地的复种指数、粮食作物播种面积和粮食单产的影响。对此，构建了最低耕地需求函数：式中，S为最小耕地需求面积，为粮食自给率，Gp为人均粮食需求量，E为复种指数，R为粮食作物播种
6、面积占农作物总播种面积的比重，U为粮食单产，P为全国的总人口。为便于动态研究我国最低耕地需求，需要预测未来不同年份的粮食单产。通过1983-2007年粮食单产的时间序列数据分析（图1），发现其变化具有一定的周期性波动，符合ARIMA模型的建模特点。通过计算机平台上的Eviews程序进行调试和优化，建立了粮食单产的ARIMA（1，1，1）模型。这样便以粮食需求为主线，通过Vensim_PLE程序建立了粮食需求与耕地资源配置的系统仿真模型（图2）。二、模型的动态模拟与调控方案的确定（一）模型调试与调控变量的确定模型调试就是通过模型的模拟运行，改进模型中不合理的反馈结构，优化模型的控制变量。一般认为
7、，模型的结构是决定模型行为的主要因素。在调试过程中，如果不改变模型的结构，当控制变量在合理的范围内取值时，发现系统的行为模式相对稳定，系统的变化表现为输出结果的变化，就认定模型具有较好的稳定性，可以用来模拟粮食需求与耕地资源配置。在综合分析近20年来我国人口数量、粮食单产、耕地变化和经济发展水平等数据的基础上，将对系统的目标函数起决定影响的人口自然增长率、城乡人均口粮直接消费量、饲料转换效率、工业化率、粮食单产、粮食作物播种比重和粮食自给标准等变量设定为调控变量（表1），其背景值的确定在参考全国土地利用总体规划纲要（2006-2020年）、国家粮食安全中长期规划纲要、中国统计年鉴、中国农业统计
8、年鉴和已有研究成果的基础上，还对预测年份的粮食单产、工业化率等相关变量采用了自回归移动平均和回归分析等处理手段，从而满足模型动态仿真模拟的需要。经常反复调试，得到了比较完整的仿真模型。经过检验，仿真模拟的结果与实际数据基本吻合，相对误差小于5%，说明模型较好地反映了系统的本质特征，有足够的真实性和有效性，能够满足研究需要，是一个具有较低数值敏感度的模型。（二）调控方案的比选粮食自给水平的高低关系到粮食安全与社会稳定。虽然学术界对我国粮食自给率标准上存在争议，但近两年来全球粮食市场供需吃紧，国际粮价连创新高，在亚非拉的数个国家引发了社会骚乱、经济动荡和政治危机。因此，我国粮食自给标准的确定需要紧
9、密结合我国人多地少的复杂国情。在此基础上，本文确定了95%和100%两种粮食自给标准的调控方案。三、结果与分析（一）粮食需求预测人口数量直接决定粮食需求。以2006年的人口数量为基数，以该年5.28的人口自然增长率为控制变量，可以预测2010、2015和2020年的总人口分别为13.42、13.78和14.15亿。从城乡结构来看（表2），城市化水平在不断提高，即农村人口的比重在减少，而城镇人口的比重在增加。又由于农村人均口粮消费量高于城镇，所以，口粮的需求量相对稳定，但在总需求中的比重相对减少，由2010年的50.8%减少到2020年的47.9%。城市化水平提高的另一个显著效应就是增加了对肉、
10、蛋、奶和水产品的需求，这就决定了饲料粮的需求量在不断增加。仿真模拟结果表明（表2），2020年的饲料粮需求量比2010年增加1257.5万吨，增长了7.5%。种子粮的需求变化比较平稳，基本稳定在1210万吨的水平上。根据预测工业化率的ARIMA（1，1，4）模型，工业化水平将由2010年的47%提高到2020年的52.5%。这在很大程度上带动了白酒、啤酒、乙醇和味精等粮食工业的快速发展，对粮食的消耗量不断增加，年均增长2.5%。在预测口粮、饲料粮、种子粮和工业用粮的基础上，便可计算粮食的总需求（表2）。从需求结构的时序变化来看，粮食总需求的增长主要是由饲料粮和工业用粮需求增加引起的，这在很大程
11、度上说明新增粮食需求量取决于饲料粮和工业用粮。（二）耕地资源配置方案 实现我国的粮食安全就需要保障耕地安全和优化配置耕地资源。尽管2007年末我国耕地资源总量为18.26亿亩，距18亿亩耕地的红线还有一定距离，但是建设用地挤占农田的现象在短期内还不可避免，这就表明耕地减少的惯性依然存在。从单位耕地面积的产出来看，粮食作物优于经济作物。如果粮食作物播种比重过高，虽有利于保障粮食安全，但不利于提高农业收益；如果粮食作物播种面积过低，虽能提高农业收益，但不利于粮食安全。因此，粮食作物的播种面积需要占适当的比例。在这样的背景下，本文模型模拟了100%和95%两种自给标准下的耕地需求与耕地配置。1.10
12、0%粮食自给标准下的耕地配置 粮食单产的变化是影响粮食综合生产能力的重要因素，根据其变动趋势的ARIMA（1，1，1）模型，可以预测不同年份的耕地需求。据预测，2010年的粮食单产为320公斤亩，根据系统动力学模型优化的粮食作物播种比重（表3），在2006年的复种指数条件下，耕地需求为181920.3万亩，如果耕地存量为182000万亩，尚有79.7万亩耕地盈余。2015年的粮食单产为336公斤亩，耕地需求为180768.2万亩，如果耕地保有量为181000万亩，尚有231.8万亩盈余。2020年的粮食单产将达到350公斤亩，在粮食作物播种比重为67.7%的条件下，耕地需求为179829.7万
13、亩，如果耕地保有量为180000万亩，尚有170.3万亩盈余。盈余的肥沃耕地可以发展优质高效的经济作物，贫瘠的耕地可以适度退耕还林（草）。2.95%粮食自给标准下的耕地配置 根据预测粮食单产的ARIMA（1，1，1）模型和系统优化的粮食作物播种比重（表3），参照当前条件下的耕地复种指数，可以量化预测年份的耕地需求。据预测，2010年的耕地需求为181781.5万亩，如果耕地保有量为182000万亩，尚有218.5万亩盈余。2015年的耕地需求为180754.2万亩，比2010年的需求量低0.56%。如果耕地保有量为181000万亩，尚有245.8万亩盈余。2020年的耕地需求为179787.7
14、万亩，比2015年的需求低966.5万亩。如果坚守18亿亩的耕地保有量，将会有212. 3万亩的盈余（表3）。通过100%和95%两种粮食自给率模拟，说明我国坚守18亿亩耕地“红线”的科学性、紧迫性和必要性。四、结论（一）粮食需求增长的总趋势不变粮食是社会稳定和经济发展的基础。虽然口粮的需求相对减少，但由于城市化水平的不断提高，对肉、蛋、奶和水产品需求不断增加，由此产生了饲料粮需求的增长。今后20年，我国粮食作物的播种面积相对稳定，因此，种子粮的需求也相对稳定。相反，工业的稳步发展带动了粮食加工业的快速发展，工业用粮将保持年均2.5%的增长速度。纵观粮食需求变化，饲料粮和工业用粮的增加带动了粮
15、食需求总量的增加，粮食需求总量将由2010年的5.16亿吨增加到2020年的5.48亿吨。（二）必须恪守粮食安全的“三重”底线 保障粮食安全是一项长期的战略任务，对此，必须恪守粮食安全的“三重”底线。一是，全国耕地面积底线。有效的耕地面积是确保粮食安全的基础，只有耕地供给得到保障，才能实现粮食安全。据此，2010年的耕地面积不低于18.2亿亩，粮食作物播种比重不低于65.5%；2015年的耕地面积不低于18.1亿亩，粮食作物播种比重不低于64.7%；而2020年的耕地面积不低于18亿亩，粮食作物播种比重不低于64.3%。二是，单产底线。粮食单产的提升是提高粮食综合生产能力的必要条件。据此，在2
16、010年的单产水平要达到320公斤亩，2015年的单产要达到326公斤亩，到2020年，粮食单产要达到350公斤亩。三是，人口底线。人口的增长是粮食需求增加的直接原因，必须实施谨慎的人口生育政策，严格控制人口增长。据此，要把今后10年的人口自然增长率控制在5.3以内，确保2010年的总量不突破13.42亿，2015年的总量不突破13.78亿，把2020年的人口规模控制在14.15亿以内。（三）严格管控，确保粮食安全 粮食安全关系国家安全。实施最严格的土地管制制度，进一步强化支农惠农政策，是实现粮食生产稳步增长的先决条件。在此基础上，创新粮食储备机制，实行“动静”相结合的储备机制，严禁使用粮食提
17、取生物燃料，工业用粮只能消耗储藏两年以上的储备粮，严防消耗当年新粮。此外，实行适度从严的市场监管。只有严防走私、囤积居奇等扰乱市场经济秩序等现象的发生，才能保障我国的粮食安全，才能在风云变幻的国际局势中始终处于主动地位。（四）设立国家粮食生产主体功能区，构筑生存安全保障线 鉴于耕地资源相对数量的减少和粮食需求总量的刚性增长，有必要依照耕地资源的空间分布和粮食生产潜力，在东北平原、黄淮海平原和长江中下游平原和若干生产潜力大的地区设立国家粮食生产主体功能区，为国家的粮食需求提供基础安全保障。该区包括8亿亩基本农田，通过加强农田基础设施建设，完善排涝灌溉系统，健全病虫害预警防治体系，提高耕地的抗灾能
18、力和应对特大自然灾害的能力，为粮食的高产、稳产提供基础条件。与此同时，在国家粮食生产主体功能区内实施合理的区域分工，发挥地区粮食生产的比较优势和专业化优势，进行专业化、集约化、现代化的粮食生产，将有利于提高我国粮食的综合生产能力和保障供给能力，这不仅是应对粮食“海啸”的战略抉择，还是中华民族实现百年生存安全的保障线。参考文献： 石玉林，封志明.开展农业资源高效利用研究J.自然资源学报，1997，12（4）：293298. 刘晓俊，李春萍，侯 聪.我国粮食需求分析与预测J.金融教学与研究，2006，（3）：34，49. 王西萍。青海省粮食需求趋势预测J.青海科技，1998，5（3）：1217.
19、何忠伟.中国粮食供求模型及其预测研究J.北京电子科技学院学报，2005，13（1）：1922. 韩 俊.当前我国粮食供求形势分析OL.中国农垦推广，2004，（2）：1112. 徐建华.现代地理学中的数学方法（第二版）M.北京：高等教育出版社，2004.370390. 张利庠.中国饲料产业发展报告M.中国农业出版社，2006. 刘晓俊.我国粮食需求分析与预测J.金融教学与研究，2006，（3）：34，49. （澳）郜若素，马国南.中国粮食研究报告M.北京农业大学出版社，1993. 孙宗风.系统动力学在水资源管理中的应用J.水利水电技术，2005，36（6）：1421. 陈百明，周小萍.中国粮食
20、自给率与耕地资源安全底线的探讨J.经济地理，2005，25（2）：145148. 张效军，欧名豪，李景刚，等.中国区域耕地赤字盈余预测J.经济学家，2006，（3）：41-48. 刘振伟.我国粮食安全的几个问题J.农业经济问题，2004，（12）：813. 中国粮食经济学会.国家粮食安全新战略研究和政策建议J.中国粮食经济，2005，（3）：813. 马永欢，牛文元，汪云林，等.我国粮食生产的空间差异与安全战略J.中国软科学，2008，（9）：1-9.作者：中国科学院 科技政策与管理科学研究所 马永欢 中国科学院 自然与社会交又科学研究中心 牛文元 来源：中国软科学2009年第3期Editor
21、s note: Judson Jones is a meteorologist, journalist and photographer. He has freelanced with CNN for four years, covering severe weather from tornadoes to typhoons. Follow him on Twitter: jnjonesjr (CNN) - I will always wonder what it was like to huddle around a shortwave radio and through the crack
22、ling static from space hear the faint beeps of the worlds first satellite - Sputnik. I also missed watching Neil Armstrong step foot on the moon and the first space shuttle take off for the stars. Those events were way before my time.As a kid, I was fascinated with what goes on in the sky, and when
23、NASA pulled the plug on the shuttle program I was heartbroken. Yet the privatized space race has renewed my childhood dreams to reach for the stars.As a meteorologist, Ive still seen many important weather and space events, but right now, if you were sitting next to me, youd hear my foot tapping rap
24、idly under my desk. Im anxious for the next one: a space capsule hanging from a crane in the New Mexico desert.待添加的隐藏文字内容2Its like the set for a George Lucas movie floating to the edge of space.You and I will have the chance to watch a man take a leap into an unimaginable free fall from the edge of
25、space - live.The (lack of) air up there Watch man jump from 96,000 feet Tuesday, I sat at work glued to the live stream of the Red Bull Stratos Mission. I watched the balloons positioned at different altitudes in the sky to test the winds, knowing that if they would just line up in a vertical straig
26、ht line we would be go for launch.I feel this mission was created for me because I am also a journalist and a photographer, but above all I live for taking a leap of faith - the feeling of pushing the envelope into uncharted territory.The guy who is going to do this, Felix Baumgartner, must have tha
27、t same feeling, at a level I will never reach. However, it did not stop me from feeling his pain when a gust of swirling wind kicked up and twisted the partially filled balloon that would take him to the upper end of our atmosphere. As soon as the 40-acre balloon, with skin no thicker than a dry cle
28、aning bag, scraped the ground I knew it was over.How claustrophobia almost grounded supersonic skydiverWith each twist, you could see the wrinkles of disappointment on the face of the current record holder and capcom (capsule communications), Col. Joe Kittinger. He hung his head low in mission contr
29、ol as he told Baumgartner the disappointing news: Mission aborted.The supersonic descent could happen as early as Sunday.The weather plays an important role in this mission. Starting at the ground, conditions have to be very calm - winds less than 2 mph, with no precipitation or humidity and limited
30、 cloud cover. The balloon, with capsule attached, will move through the lower level of the atmosphere (the troposphere) where our day-to-day weather lives. It will climb higher than the tip of Mount Everest (5.5 miles/8.85 kilometers), drifting even higher than the cruising altitude of commercial ai
31、rliners (5.6 miles/9.17 kilometers) and into the stratosphere. As he crosses the boundary layer (called the tropopause), he can expect a lot of turbulence.The balloon will slowly drift to the edge of space at 120,000 feet (22.7 miles/36.53 kilometers). Here, Fearless Felix will unclip. He will roll
32、back the door.Then, I would assume, he will slowly step out onto something resembling an Olympic diving platform.Below, the Earth becomes the concrete bottom of a swimming pool that he wants to land on, but not too hard. Still, hell be traveling fast, so despite the distance, it will not be like div
33、ing into the deep end of a pool. It will be like he is diving into the shallow end.Skydiver preps for the big jumpWhen he jumps, he is expected to reach the speed of sound - 690 mph (1,110 kph) - in less than 40 seconds. Like hitting the top of the water, he will begin to slow as he approaches the m
34、ore dense air closer to Earth. But this will not be enough to stop him completely.If he goes too fast or spins out of control, he has a stabilization parachute that can be deployed to slow him down. His team hopes its not needed. Instead, he plans to deploy his 270-square-foot (25-square-meter) main
35、 chute at an altitude of around 5,000 feet (1,524 meters).In order to deploy this chute successfully, he will have to slow to 172 mph (277 kph). He will have a reserve parachute that will open automatically if he loses consciousness at mach speeds.Even if everything goes as planned, it wont. Baumgar
36、tner still will free fall at a speed that would cause you and me to pass out, and no parachute is guaranteed to work higher than 25,000 feet (7,620 meters).It might not be the moon, but Kittinger free fell from 102,800 feet in 1960 - at the dawn of an infamous space race that captured the hearts of many. Baumgartner will attempt to break that record, a feat that boggles the mind. This is one of those monumental moments I will always remember, because there is no way Id miss this.
