基于标准样地的省级耕地质量监测样地布设方法——以内蒙古自治区为例

doi:10.18402/resci.2016.11.03

资源科学 ›› 2016, Vol. 38 ›› Issue (11): 2037-2048.doi: 10.18402/resci.2016.11.03

基于标准样地的省级耕地质量监测样地布设方法——以内蒙古自治区为例

张玉臻^1,², 孔祥斌^1,²(), 刘炎³, 张蚌蚌^1,², 张青璞^1,², 王峰³, 李立强⁴, 魏利利⁴

1. 中国农业大学资源与环境学院,北京 100193
2. 国土资源部农用地质量与监测重点实验室,北京 100193
3. 国土资源部咨询研究中心,北京 100035
4. 内蒙古自治区土地整治中心,呼和浩特 010020

收稿日期:2016-06-14 修回日期:2016-09-29 出版日期:2016-11-16 发布日期:2016-11-16
作者简介:
作者简介：张玉臻,女,山东潍坊市人,硕士生,研究方向为持续土地利用与土地评价。E-mail：zhangyuzhen66@126.com
基金资助:
国土资源部公益性行业科研专项课题（201011006-3）

Layout method of monitoring samples of cultivated land at the provincial level based on standard plots in Inner Mongolia

ZHANG Yuzhen^1,², KONG Xiangbin^1,²(), LIU Yan³, ZHANG Bangbang^1,², ZHANG Qingpu^1,², WANG Feng³, LI Liqiang⁴, WEI Lili⁴

1.College of Resources and Environmental Sciences ,China Agricultural University,Bejing 100193,China
2. Key Laboratory for Farmland Quality,Monitoring and Control,National Ministry of Land and Resources,Beijing 100193,China
3. Consulting and Research Center Ministry of Land and Resources,Beijing 100035,China
4. Land Consolidation and Rehabilitation Center,The Inner Mongolia Autonomous Region,Huhehaote 010020,China

Received:2016-06-14 Revised:2016-09-29 Online:2016-11-16 Published:2016-11-16

摘要/Abstract

摘要：

耕地质量是国家粮食安全的基础,科学布设耕地质量监测样地是支撑实时掌握国家耕地质量的关键技术环节。当前对耕地质量监测的研究主要集中在国家级和县级尺度,缺乏省级尺度。本文以内蒙古自治区为例,提出在省级标准样地基础上进行耕地监测样地数量和质量代表性两方面的优化,进而形成省级耕地质量监测布点的技术方法。本文运用最邻近点指数和面积代表性指数对标准样地的空间分布和数量进行分析和验证,结果表明,原144个标准样地在空间上呈均匀分布,但在数量上不能满足监测布样的要求;运用地统计学方法确定理论监测样地数量为194个;以“耕作制度二级区-土壤类型-地貌-土地利用现状”为基础,结合《中国土地资源图集》中的中国土地潜力分区图,将内蒙古划分为14个监测控制分区;依据“等别-监测分区”的控制顺序和优先选择标准样地的原则修正监测样地数量和确定其空间位置,最终确定省级监测样地数量为219个,其中标准样地124个,新增监测样地95个。依据本文提出的方法进行省级监测样地布设,一方面可以建立起省级尺度符合耕地质量空间变异特征的耕地质量监测布设样地体系,另一方面可以促进国家和县级耕地质量监测工作之间的衔接。

关键词: 标准样地, 省级尺度, 耕地质量, 监测样地, 监测控制分区, 内蒙古自治区

Abstract:

Food security depends on the quality and quantity of arable land. The implementation of dynamic balance of total amount of arable land has made a significant impact to the quantity of farmland,but less regarding the quality of cultivated land protection. Farmland quality monitoring includes three levels：national,provincial and county levels. At present,domestic research mostly focusses at the county scale and national level,and provincial level research for Inner Mongolia is lacking. Here we propose optimization of the quantity and quality of standard samples,and form a technical method for monitoring sample layout at the provincial level. Inner Mongolia setting 144 provincial standard samples,we verify standard sample number and spatial distribution using the nearest neighbor index and area representative index. The results show that the spatial distribution of the original standard samples is uniform,but the distribution of number does not meet monitoring requirements. The theoretical numbers of monitoring samples are 194 using the geo-statistical method. Based on ‘secondary zones-soil types-landscapes-land use status’ combined with the distribution of land resource potential,we analyze limiting factors of land use and divide the area into 14 monitoring zones. According to the principles of one monitoring sample in different monitoring partitions in one farmland gradation at least and area index,the number of monitoring samples is revised to 219,made up of 124 original standard samples and 95 new samples. Layout of provincial monitoring samples based on our method minimizes the cost to meet the requirements of provincial arable land quality monitoring and can build arable land quality monitoring systems at the provincial scale in line with spatial variability of farmland quality. This approach provides a scientific basis for the national implementation of cultivated land quality monitoring.

Key words: standard farmland, provincial scale, cultivated land quality, monitoring sample, monitoring partion, Inner Mongolia

张玉臻, 孔祥斌, 刘炎, 张蚌蚌, 张青璞, 王峰, 李立强, 魏利利. 基于标准样地的省级耕地质量监测样地布设方法——以内蒙古自治区为例[J]. 资源科学, 2016, 38(11): 2037-2048.

ZHANG Yuzhen,KONG Xiangbin,LIU Yan,ZHANG Bangbang,ZHANG Qingpu,WANG Feng,LI Liqiang,WEI Lili. Layout method of monitoring samples of cultivated land at the provincial level based on standard plots in Inner Mongolia[J]. Resources Science, 2016, 38(11): 2037-2048.

图/表 12

图1

图2

表1

图3

Table 2

表3

监测样地数量优化"

	依据不同耕地质量等指数确定样本量
	自然等指数	利用等指数	经济等指数
总样本量/个	4 121	4 121	4 121
标准差	519	263	278
均值	1 303.7	610.9	572.5
$u a$	1.96	1.96	1.96
相对误差/%	5	5	5
$n 0$	243.53	284.80	362.34
传统简单随机抽样样本	229.94	266.39	333.05
MORAN’s I	0.71	0.59	0.42
显著性区间/%	>99	>99	>99
校正后监测样地数量 (取整)/个	67	110	194

表3

图4

表4

监测样地数量修正"

耕地等别	监测分区数量 $y p$ /个	耕地面积 /万hm²	面积比例 /%	合理数量 $h p$ /个	修正数量 $r p$ /个
7	2	0.15	0.02	0	2
8	5	0.56	0.08	0	5
9	6	2.76	0.37	1	6
10	11	11.50	1.55	3	11
11	11	25.20	3.40	7	11
12	14	49.70	6.71	13	14
13	14	77.30	10.44	20	20
14	14	136.00	18.37	35	35
15	14	437.00	59.04	115	115
总计	14	740.17	100.00	194	219

表4

图5

表5

表6

图6

参考文献 36

[1]	Kong X B.China must protect high-quality arable land[J]. Nature,2014,doi:10.1038/506007a.
[2]	肖建英,谭术魁. 中国粮食产量省级尺度下的空间分异规律[J]. 中国土地科学,2013,27(8):26-32.
	[Xiao J Y,Tan S K.Regularity of the spatial variation of grain yield at provincial scale in China[J]. China Land Sciences,2013,27(8):26-32.]
[3]	陈秧分,李先德. 中国粮食产量变化的时空格局与影响因素[J]. 农业工程学报,2013,29(20):1-10.
	[Chen Y F,Li X D.Spatial-temporal characteristics and influencing factors of grain yield change in China[J]. Transactions of the Chinese Society of Agricultural Engineering,2013,29(20):1-10.]
[4]	刘会玉,林振山,张明阳. 基于EMD的我国粮食产量波动及其成因多尺度分析[J]. 自然资源学报,2005,20(5):745-751.
	[Liu H Y,Lin Z S,Zhang M Y.Analysis on the fluctuation of grain output in China and its causes at multi-time scale based on empirical mode decomposition method[J]. Journal of Natural Resources,2005,20(5):745-751.]
[5]	中华人民共和国国家统计局. 2012年全国粮食生产再获丰收[N]. 中国信息报,2012-12-03(001).
	[National Bureau of Statistics of the People’s Republic of China. The National Grain Production Was a Very Good Again in2012[N]. China Infor-mation News,2012-12-03(001).]
[6]	Kong X B,Lal R,Li B G,et al. Fertilizer intensification and its impacts in China’s HHH plains[J]. Advances in Agronomy,2014,125:135-169.
[7]	李宇轩. 中国化肥产业政策对粮食生产的影响研究[D]. 北京:中国农业大学,2014
	[Li Y X.The Research on Effect of Fertilizer Industry Policy on Grain Production in China[D]. Beijing:China Agricultural University,2014.]
[8]	Lobell D B,Schlenker W,Costa-Roberts J.Climate trends and global crop production since 1980[J]. Science,2011,333(6042):616-620.
[9]	孙宏岭,韩国山. 我国粮食安全的主要指标与保障途径研究[J]. 粮食科技与经济,2015,40(1):5-7.
	[Sun H L,Han G S.The main indicators and guarantee of food security in China[J]. Food Echnology and Economy,2015,40(1):5-7.]
[10]	付国珍,摆万奇. 耕地质量评价研究进展及发展趋势[J]. 资源科学,2015,37(2):226-236.
	[Fu G Z,Bai W Q.Advances and prospects of evaluating cultivated land quality[J]. Resources Science,2015,37(2):226-236.]
[11]	欧胜彬,苏雪华. 粮食安全:耕地质量保护是关键[J]. 农业经济,2012,(8):54-56.
	[Ou S B,Su X H.Food security:The quality of cultivated land protection is the key[J]. Agricultural Economy,2012,(8):54-56.]
[12]	雷勋平,吴杨,叶松,等. 基于熵权可拓决策模型的区域粮食安全预警[J]. 农业工程学报,2012,28(6):233-239.
	[Lei X P,Wu Y,Ye S,et al. Regional grain security pre-warning based on entropy weight extension decision model[J]. Transactions of the Chinese Society of Agricultural Engineering,2012,28(6):233-239.]
[13]	王洪波,郧文聚,吴次芳,等. 基于农用地分等的耕地产能监测体系研究[J]. 农业工程学报,2008,24(4):122-126.
	[Wang H B,Yun W J,Wu C F,et al. Farmland productivity monitoring system based on farmland gradation[J]. Transactions of the Chinese Society of Agricultural Engineering,2008,24(4):122-126.]
[14]	祝锦霞,徐保根,章琳云. 基于半方差函数与等别的耕地质量监测样点优化布设方法[J]. 农业工程学报,2015,31(19):254-261.
	[Zhu J X,Xu B G,Zhang L Y.Optimization layout method of monitoring sample points of cultivated land quality based on semi-variance analysis and grade combination[J]. Transactions of the Chinese Society of Agricultural Engineering,2015,31(19):254-261.]
[15]	Nusser S M,Breidt E J,Fuller W A.Design and estimation for investigating the dynamics of natural resources[J]. Ecological Applications,1998,8(2):234-245.
[16]	Nusser S M,Goebel J J.The national resources inventory:A long-term multi-resource monitoring programme[J]. Journal of Agricul-tural Biological and Environmental Statistics,1997,4(3):181-204.
[17]	Mol G,Vriend S P,Gaans P F M V. Future trends,detectable by soil monitoring networks?[J]. Journal of Geochemical Exploration,1998,62(1-3):61-66.
[18]	Maiti T,Miller C P,Mukhopadhyay P K.Neural network imputation:An experience with the national resources inventory survey[J]. Journal of Agricultural Biological and Environmental Statistics,2008,13(3):255-269.
[19]	张蚌蚌,孔祥斌,相慧,等. 中国西部区国家级标准样地分布特征及代表性[J]. 农业工程学报,2014,30(7):227-235.
	[Zhang B B,Kong X B,Xiang H,et al. Distribution characteristics and representative test of national standard farmland in Western China[J]. Transactions of the Chinese Society of Agricultural Engi-neering,2014,30(7):227-235.]
[20]	何江华,周明中. 基于标准样地的耕地质量监测研究[J]. 河北农业科学,2010,14(9):125-127.
	[He J H,Zhou M Z.Research on monitoring of cultivated land quality based on standard plot[J]. Journal of Hebei Agricultural Sciences,2010,14(9):125-127.]
[21]	伍育鹏,郧文聚,李武艳. 用标准样地进行耕地质量动态监测与预警探讨[J]. 中国土地科学,2006,20(4):40-45.
	[Wu Y P,Yun W J,Li W Y.Research on standard-plot based monitoring and early-warning of arable land quality[J]. Resources Science,2006,20(4):40-45.]
[22]	郭力娜,张凤荣,马仁会,等. 基于标准样地的国家级农用地等别质量监测点设置方法探讨-以冀豫鄂三省为例[J]. 资源科学,2009,31(11):1957-1966.
	[Guo L N,Zhang F R,Ma R H,et al. Methods for setting up farmland grade monitoring points from standard plots in Hebei,Henan and Hubei Provinces[J]. Resources Science,2009,31(11):1957-1966.]
[23]	吴克宁,焦雪瑾,梁思源,等. 基于标准样地国家级汇总的耕地质量动态监测点构架研究[J]. 农业工程学报,2008,24(10):74-79.
	[Wu K N,Jiao X J,Liang S Y,et al. Framework of the arable-land quality dynamic monitoring sites through national gathering of standard farmland[J]. Transactions of the Chinese Society of Agricultural Engineering,2008,24(10):74-79.]
[24]	郭力娜,马仁会,徐东瑞,等. 农用地等别质量监测带布设方法探讨-以冀豫鄂三省为例[J]. 资源科学,2008,30(8):1199-1205.
	[ Guo L N,Ma R H,Xu D R,et al. Methods of setting up farmland grade monitoring belts in Hebei,Henan and Hubei[J]. Resources Science,2008,30(8):1199-1205.]
[25]	余述琼,张蚌蚌,相慧,等. 基于因素组合的耕地质量等级监测样点布控方法[J]. 农业工程学报,2014,30(24):288-297.
	[Yu S Q,Zhang B B Xiang H,et al. Layout method for monitoring sample point of arable land quality level based on combination of factors[J]. Transactions of the Chinese Society of Agricultural Engineering,2014,30(24):288-297.]
[26]	王佳宁,孔祥斌,张青璞,等. 基于农用地分等的县级耕地质量等级监测指标选取方法-以内蒙古达拉特旗为例[J]. 中国农业大学学报,2015,20(3):175-187.
	[Wang J N,Kong X B,Zhang Q P,et al. Factors selection method for monitoring the cultivated land level based on farmland grading:In Dalate County[J]. Journal of China Agricultural University,2015,20(3):175-187.]
[27]	孙亚彬,吴克宁,胡晓涛,等. 基于潜力指数组合的耕地质量等级监测布点方法[J]. 农业工程学报,2013,29(4):245-254.
	[Sun Y B,Wu K N,Hu X T,et al. Layout method for monitoring quality level of arable land based on combination of potential index[J]. Transactions of the Chinese Society of Agricultural Engi-neering,2013,29(4):245-254.]
[28]	相慧,孔祥斌,陈培雄,等. 县域耕地质量等别监测样点布控研究-以内蒙古自治区达拉特旗为例[J]. 资源科学,2014,36(6):1203-1210.
	[Xiang H,Kong X B,Chen P X,et al. Establishing monitoring samples of arable land at the county level:A case study of Dalate District in Inner Mongolia[J]. Resources Science,2014,36(6):1203-1210.]
[29]	刘霈珈,吴克宁,赵华甫. 基于农用地分等与土地质量地球化学评估的耕地质量监测类型研究[J]. 资源科学,2015,37(1):37-44.
	[Liu P J,Wu K N,Zhao H F.Cultivated land quality monitoring types base on agricultural land classification and land quality geochemical assessment[J]. Resources Science,2015,37(1):37-44.]
[30]	中华人民共和国国土资源部. 农用地分等规程(TD/T 1004-2003)[S]. 北京:中国标准出版社,2003.
	[ Ministry of Land and Resources of the People 's Republic of China. Regulation for Gradation on Agricultural Land Quality(TD/T 1004-2003)[S].Beijing:China Standard Press,2003.]
[31]	王劲峰,廖一兰,刘鑫. 空间数据分析教程[M]. 北京:科学出版社,2010.
	[Wang J F,Liao Y L,Liu X.A Course in Spatial Data Analysis[M]. Beijing:Science Press,2010.]
[32]	姜成晟,王劲峰,曹志冬. 地理空间抽样理论研究综述[J]. 地理学报,2009,64(3):368-380.
	[Jiang C S,Wang J F,Cao Z D.A review of geo-spatial sampling theory[J]. Acta Geogra Sinica,2009,64(3):368-380.]
[33]	王迪,陈仲新,周清波,等. 冬小麦种植面积空间抽样样本布局的优化设计[J]. 中国农业科学,2014,47(18):3545-3556.
	[Wang D,Chen Z X,Zhou Q B,et al. Optimization of samples layout in spatial sampling schemes for estimating winter wheat planting acreage[J]. Scientia Agricultura Sinica,2014,47(18):3545-3556.]
[34]	李连发,王劲峰,姜成晟. 空间抽样与统计[M]. 北京:科学出版社,2009.
	[Li L F,Wang J F,Jiang C S.Spatial Sampling and Statistics[M]. Beijing:Science Press,2009.]
[35]	魏晓峰,吴健平. 基于ArcGIS的空间自相关分析模块的开发与应用[J]. 测绘与空间地理信息,2005,28(6):77-80
	[Wei X F,Wu J P.The development and application of spatial auto-correlation based on ArcGIS[J]. Geomatics & Spatial Information Technology,2005,28(6):77-80.]
[36]	石玉林,岳燕珍,石竹筠.中国土地资源图集[M]. 北京:中国大地出版社,2006.
	[Shi Y L,Yue Y Z,Shi Z Y.China Land Re-sources Atlas[M]. Beijing:China Land Press,2006.]

二级区名称	二级区代称	最邻近点指数D	空间结构类型	面积/万hm²	标准样地数量/个
阿拉善高原区	ALS	6.21	均匀型	2.73	4
大小兴安岭山地区	DXXAL	4.64	均匀型	110.59	14
河套银川平原区	HTYC	2.38	均匀型	90.53	23
后山坝上高原区	HSBS	4.76	均匀型	144.93	40
辽吉西蒙东南冀北山地区	LJXM	2.99	均匀型	257.91	45
内蒙古草原区	NMG	4.70	均匀型	109.74	17
松嫩平原区	SN	-	-	23.74	1
全省	-	3.56	均匀型	740.37	144

二级区	7等		8等		9等		10等		11等		12等		13等		14等		15等		合计
二级区	面积	数量	面积	数量	面积	数量	面积	数量	面积	数量	面积	数量	面积	数量	面积	数量	面积	数量	样本总数
阿拉善高原区	-	-	-	-	0.31	-	1.06	2	0.96	2	0.40	-	-	-	-	-	-	-	4
大小兴安岭山地区	-	-	-	-	-	-	-	-	-	-	0.76	-	5.52	2	28.5	4	75.81	8	14
河套银川平原区	0.13	-	0.37	2	2.16	2	10.20	3	22.15	4	22.18	6	10.04	2	8.89	2	14.41	2	23
后山坝上高原区	0.02	2	0.04	-	0.06	2	0.16	2	0.89	4	7.04	1	11.16	5	27.86	12	97.7	12	40
辽吉西蒙东南冀北山地区	-	-	-	-	0.21	2	0.13	-	1.17	1	17.40	3	22.30	13	47.70	14	169	12	45
内蒙古草原区	-	-	0.15	-	0.02	-	0	-	0.03	-	1.94	3	28.30	2	15.60	3	63.7	9	17
松嫩平原区	-	-	-	-	-	-	-	-	-	-	-	-	0.01	-	7.31	-	16.62	1	1
总计	0.15	2	0.56	2	2.76	6	11.50	7	25.20	11	49.70	13	77.30	24	136	35	437	44	144

监测控制区编号	名称	面积比例 /%	实际样本量 /个	保留标准样地 /个	新增监测样地 /个
J-01	阿拉善高原-多限制因素区	0.12	4	4	0
J-02	东北及松嫩平原-较少限制区	7.35	15	4	11
J-03	东北-坡度限制区	8.02	23	10	13
J-04	河套银川平原-坡度与侵蚀限制区	5.50	12	9	2
J-05	河套银川平原-盐碱限制区	7.53	13	10	3
J-06	后山坝上高原-坡度侵蚀限制区	10.80	23	11	12
J-07	后山坝上高原-水分限制区	10.35	10	9	1
J-08	后山坝上高原-土质限制区	1.49	13	10	3
J-09	辽吉西蒙东南冀北山地-坡度与土质限制区	24.70	54	26	29
J-10	辽吉西蒙东南冀北山地-轻微土质限制区	13.10	18	10	8
J-11	内蒙古草原-水分土质限制区	8.72	18	11	7
J-12	内蒙古草原-土质限制区	0.86	8	5	3
J-13	内蒙古草原-水分限制区	0.37	5	3	2
J-14	内蒙古草原-多种限制因素区	1.09	3	2	1
总计	-	100	219	124	95

监测分区	类别	监测样地序号	标准样地序号	利用等别
河套银川平原	标准样地	JCYD-89	SJ-123	14
坡度与侵蚀		JCYD-90	SJ-124	12
限制区		JCYD-91	GJJ-125	8
		JCYD-96	GJJ-128	13
		JCYD-99	GJJ-136	12
		JCYD-100	SJ-137	11
		JCYD-101	SJ-138	12
		JCYD-104	GJJ-84	10
		JCYD-108	SJ-129	13
	新增监测	JCYD-103	-	14
	样地	JCYD-146	-	15
		JCYD-147	-	15

基于标准样地的省级耕地质量监测样地布设方法——以内蒙古自治区为例

Layout method of monitoring samples of cultivated land at the provincial level based on standard plots in Inner Mongolia

RichHTML

PDF (PC)

赞

可视化

摘要/Abstract

引用本文

使用本文

图/表 12

参考文献 36

相关文章 15

Metrics

本文评价

推荐阅读 0

[1]	张亚丽, 白云丽, 辛良杰. 耕地质量与土地流转行为关系研究[J]. 资源科学, 2019, 41(6): 1102-1110.
[2]	杨江燕, 殷守强, 张利, 门明新, 陈影. 基于空间聚类分层抽样的黄骅市县域耕地质量等别监测样点布设[J]. 资源科学, 2019, 41(2): 257-267.
[3]	辛芸娜, 范树印, 孔祥斌, 张蚌蚌, 温良友, 张青璞. 四重质量维度下的县域耕地质量评价方法研究[J]. 资源科学, 2018, 40(4): 737-747.
[4]	杜挺, 朱道林. 中国土地流转价格时空演化与宏观机制研究[J]. 资源科学, 2018, 40(11): 2202-2212.
[5]	慕正隆, 石培基, 刘春芳, 郑海松. 基于形态优化和边界识别的基本农田保护区研究[J]. 资源科学, 2018, 40(1): 22-31.
[6]	龙云, 任力. 农地流转制度对农户耕地质量保护行为的影响——基于湖南省田野调查的实证研究[J]. 资源科学, 2017, 39(11): 2094-2103.
[7]	程子良, 孔祥斌, 张蕾娜, 史婧然, 程佳. 基于质量、数量和经济发展水平的耕地占用税征收标准研究[J]. 资源科学, 2017, 39(1): 61-73.
[8]	殷守强, 王鑫, 贺文龙, 门明新, 张利. 耕地质量监测中不同抽样方法的精度比较——以河北省黄骅市为例[J]. 资源科学, 2016, 38(11): 2049-2057.
[9]	陈藜藜, 宋戈, 邹朝晖, 周浩, 张旭. 基于改进局部空间自相关分析的萝北县耕地保护分区[J]. 资源科学, 2016, 38(10): 1871-1882.
[10]	魏洪斌, 吴克宁, 赵华甫, 夏敏峰, 刘巧芹. 中国中部粮食主产区耕地等别空间分布特征[J]. 资源科学, 2015, 37(8): 1552-1560.
[11]	侯现慧, 王占岐, 杨俊. 富硒区耕地质量评价及利用分区研究——以福建省三元区为例[J]. 资源科学, 2015, 37(7): 1367-1375.
[12]	白雪红, 闫慧敏, 黄河清, 胡云峰, 匡文慧. 1991-2010年内蒙古耕地转出时空格局及分布规律[J]. 资源科学, 2015, 37(2): 360-369.
[13]	刘霈珈, 吴克宁, 赵华甫. 基于农用地分等与土地质量地球化学评估的耕地质量监测类型研究[J]. 资源科学, 2015, 37(1): 37-44.
[14]	胡晓涛, 马建辉, 王倩, 吴克宁, 赵华甫. 北京市大兴区耕地质量等级监测控制点布设[J]. , 2012, 34(10): 1891-1897.
[15]	金晓斌, 宋佳楠, 周寅康. 基于多层线性模型的耕地集约利用对粮食生产力贡献度分析[J]. , 2010, 32(6): 1161-1168.