劉克寶，陸忠軍，付 斌，辛 蕊，黃 楠，劉艷霞，李 楊

（1. 黑龍江省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)遙感與信息研究所，哈爾濱150086；2. 農(nóng)業(yè)農(nóng)村部農(nóng)業(yè)遙感重點(diǎn)實(shí)驗(yàn)室，北京100081）

0 引言

農(nóng)作物空間格局是一個(gè)地區(qū)或生產(chǎn)單位作物種植結(jié)構(gòu)、熟制與種植方式的空間表達(dá)［1］，是了解農(nóng)作物種類、結(jié)構(gòu)、分布特征的重要信息，也是進(jìn)行作物結(jié)構(gòu)調(diào)整和優(yōu)化的依據(jù)［2］，其時(shí)空分布受自然資源條件和市場經(jīng)濟(jì)行為影響，具有空間集聚性和時(shí)序動(dòng)態(tài)性［3］，其特征及時(shí)空動(dòng)態(tài)變化是農(nóng)業(yè)地理和可持續(xù)農(nóng)業(yè)發(fā)展的重要基礎(chǔ)，特別是在農(nóng)業(yè)生態(tài)系統(tǒng)的陸地碳排放、全球變化對區(qū)域農(nóng)業(yè)生產(chǎn)、農(nóng)業(yè)發(fā)展的影響以及作物模式的動(dòng)力學(xué)機(jī)理及模擬模型等方面［4］。

可持續(xù)的農(nóng)作物空間格局有利于糧食安全、農(nóng)業(yè)資源的利用和保護(hù)，為保障國家糧食安全，保護(hù)農(nóng)業(yè)資源，國家自2014年相繼出臺(tái)了相關(guān)的文件，一是政策方面：（1）種植結(jié)構(gòu)調(diào)整政策：2015年11月農(nóng)業(yè)部下發(fā)《鐮刀灣地區(qū)玉米結(jié)構(gòu)調(diào)整的意見》，調(diào)減玉米面積；2017年2月5日，國家公開發(fā)布《關(guān)于深入推進(jìn)農(nóng)業(yè)供給側(cè)結(jié)構(gòu)性改革加快培育農(nóng)業(yè)農(nóng)村發(fā)展新功能的若干意見》，統(tǒng)籌調(diào)整糧經(jīng)飼種植結(jié)構(gòu)，優(yōu)化農(nóng)業(yè)區(qū)域布局；（2）耕地保護(hù)政策：2015年，黑土地保護(hù)利用試點(diǎn)項(xiàng)目實(shí)施；2016年5月20日，《探索實(shí)行耕地輪作休耕制度試點(diǎn)方案》實(shí)施。二是補(bǔ)貼方面：2014年實(shí)施了大豆目標(biāo)價(jià)格補(bǔ)貼試點(diǎn)，自2017年起，將大豆目標(biāo)價(jià)格補(bǔ)貼調(diào)整為大豆生產(chǎn)者補(bǔ)貼，與玉米生產(chǎn)者補(bǔ)貼統(tǒng)籌銜接；水稻最低保護(hù)價(jià)政策從2017年開始降低。在國家供給側(cè)結(jié)構(gòu)改革和土地資源保護(hù)政策相繼出臺(tái)的背景下，對作物時(shí)空格局變化進(jìn)行分析，對國家農(nóng)業(yè)政策的制定和實(shí)施具有重要的理論和實(shí)際意義［5］。

農(nóng)作物空間格局及其時(shí)空動(dòng)態(tài)變化研究已經(jīng)成為地理學(xué)和生態(tài)學(xué)的前沿和熱點(diǎn)研究問題，其信息獲取主要有3種方法：基于統(tǒng)計(jì)數(shù)據(jù)的方法、基于遙感信息的方法以及基于空間模型的方法。劉珍環(huán)等［6］利用縣級(jí)農(nóng)業(yè)統(tǒng)計(jì)數(shù)據(jù)，綜合運(yùn)用時(shí)序變化趨勢、空間集聚分析等方法，從種植結(jié)構(gòu)類型和種植比例變化趨勢分析了1980年以來中國縣域種植結(jié)構(gòu)的時(shí)空特征。潘曉卉等［7］基于2005—2015年東北地區(qū)40個(gè)市/盟主要糧食作物播種面積數(shù)據(jù)（統(tǒng)計(jì)數(shù)據(jù)），分析東北地區(qū)主要糧食作物種植結(jié)構(gòu)時(shí)空變化。余強(qiáng)毅等［8］在農(nóng)作物空間格局變化模擬模型（an agent-based model for simulating crop pattern dynamics，CroPaDy）概念化設(shè)計(jì)的基礎(chǔ)上，借助MATLAB平臺(tái)開放性、矩陣運(yùn)算能力強(qiáng)等特點(diǎn)，實(shí)現(xiàn)CroPaDy模型的數(shù)值模擬，并以黑龍江省賓縣調(diào)查數(shù)據(jù)為依據(jù)，完成模型的區(qū)域?qū)嵶C；夏天等［9］研究并提出農(nóng)作物空間格局動(dòng)態(tài)變化模型（Crop Pattern Simulator，CROPS）的可行架構(gòu)，并選取我國東北三省進(jìn)行試驗(yàn)，耕地空間格局模擬精度達(dá)到90%，農(nóng)作物空間格局模擬精度達(dá)到75%；譚杰揚(yáng)等［10］、劉珍環(huán)等［11］利用作物空間分配模型（SPAM—China）模擬東北三省作物時(shí)空分布特征。劉航等［12］利用Landsat TM/OLI遙感影像，以遙感、地理信息系統(tǒng)技術(shù)為支撐，利用隨機(jī)森林方法提取大豆空間分布信息，并結(jié)合景觀格局指數(shù)分析松嫩平原大豆種植結(jié)構(gòu)時(shí)空演變特征。陸忠軍等［13］利用人工目視解譯中高分辨率衛(wèi)星數(shù)據(jù)（TM5、SPOT2/4），結(jié)合自然生態(tài)背景（坡度、積溫、地表徑流），分時(shí)段分析1993—2011年黑龍江省水稻時(shí)空變化。

以上基于統(tǒng)計(jì)數(shù)據(jù)和空間模型的作物空間格局變化分析的研究中，均未描述作物之間的面積相互轉(zhuǎn)換以及種植方式（輪作、連作）；而基于遙感手段的作物種植結(jié)構(gòu)時(shí)空變化分析多以長時(shí)段年間節(jié)點(diǎn)主要作物或單一作物時(shí)空變化分析為主，忽略了農(nóng)業(yè)政策背景下年際間作物空間格局的劇烈變化，同樣無法獲得作物種植方式（輪作、連作）。基于此，文章選擇松嫩平原中部的海倫市為研究區(qū)域，以2014—2018年為研究時(shí)段，以中高分辨率衛(wèi)星數(shù)據(jù)提取的農(nóng)作物空間分布數(shù)據(jù)為基礎(chǔ)，系統(tǒng)深入分析該時(shí)段內(nèi)年際間作物種植結(jié)構(gòu)時(shí)空變化特征及種植方式變化，為黑龍江省種植結(jié)構(gòu)調(diào)整與優(yōu)化布局提供科學(xué)依據(jù)和對策以及土地可持續(xù)利用建議。

1 材料與方法

1.1 研究區(qū)域

松嫩平原是國家重要的商品糧基地，是國家種植結(jié)構(gòu)調(diào)整的重點(diǎn)區(qū)域，農(nóng)民種植作物主要受大田作物間比較效益、國家政策的影響，但目前針對該區(qū)域作物空間格局的研究較少。

圖1 2014—2018年研究區(qū)位置和作物空間分布圖：a.研究區(qū)位置；b.作物空間分布Fig.1 Location of study area and crops spatial distribution map from 2014 to 2018

為兼顧研究的典型性和代表性，該研究選擇海倫市為研究區(qū)域，海倫市隸屬于黑龍江省綏化市，地處松嫩平原和小興安嶺余脈的過度地帶，松嫩平原中部，是世界三大黑土地之一，屬寒溫帶大陸性氣候，四季分明，降水集中，有效積溫2 200～2 400℃，年降水量500～600 mm。主要農(nóng)作物有大豆、玉米、水稻等糧食作物和亞麻、烤煙等經(jīng)濟(jì)作物，享有中國優(yōu)質(zhì)大豆之鄉(xiāng)、中國高淀粉玉米之鄉(xiāng)、中國甜菇娘之鄉(xiāng)、中國甜菜之鄉(xiāng)、中國籽鵝之鄉(xiāng)“六鄉(xiāng)”美譽(yù)，是國家重要的商品糧基地縣，是大豆、玉米種植的交錯(cuò)地帶，是2015年東北4省（區(qū)）17個(gè)縣（市）開展黑土地保護(hù)利用試點(diǎn)項(xiàng)目縣之一。

1.2 數(shù)據(jù)

行政邊界線采用2012年全國1∶100萬分縣行政邊界線；作物空間分布數(shù)據(jù)為遙感手段獲得的2014、2015、2016、2017和2018年共5年的遙感數(shù)據(jù)，該數(shù)據(jù)的耕地和水稻采用人工目視解譯獲得，旱田作物（大豆、玉米和其他作物）采用計(jì)算機(jī)自動(dòng)分類輔助人工目視解譯，其中2014年和2015年作物空間分布使用的衛(wèi)星數(shù)據(jù)源為Landsat8 OLI（30 m），2016、2017和2018年3年使用Sentinel-2A（10 m）為主要數(shù)據(jù)源，Landsat8 OLI（30 m）為補(bǔ)充數(shù)據(jù)源。

1.3 分析方法

該文綜合運(yùn)用數(shù)理統(tǒng)計(jì)和GIS空間分析等方法，分析年際間主要農(nóng)作物類型變化量及幅度、輪作方式等方面的時(shí)空變化規(guī)律。具體如下。

（1）作物種植比例S，反映耕地內(nèi)部種植結(jié)構(gòu)特征。

式（1）中，Ai是作物類型i的總面積，其中A1+A2+…+Ai=100。

（2）作物類型變化幅度ΔS。表示某種作物在一段時(shí)期內(nèi)的面積變化比例。

式（2）中，S（i，t1）和S（i，t2）分別為研究初期t1和末期t2某一作物類型i的面積。

（3）作物變化方向時(shí)空變化分析。作物變化方向時(shí)空變化分析為分別疊加鄰近兩年的作物空間分布數(shù)據(jù)，分析各作物（水稻、玉米、大豆和其他作物）不變區(qū)域、增加區(qū)域和減少區(qū)域。

（4）旱田作物輪作方式時(shí)空變化分析。疊加年際間作物空間分布數(shù)據(jù)，將公共區(qū)域內(nèi)的作物輪作方式劃分為：玉米—大豆、玉米—其他作物、大豆—玉米、大豆—其他作物、其他作物—玉米和其他作物—大豆6種輪作方式，并進(jìn)行年際間面積統(tǒng)計(jì)分析。

（5）旱田作物輪作次數(shù)分析。疊加2014、2015、2016、2017和2018年作物空間分布結(jié)果，選取5年旱田公共區(qū)域，將旱田作物輪作次數(shù)分為0、1、2、3和4次，并對面積進(jìn)行統(tǒng)計(jì)，分析各輪作次數(shù)中各輪作時(shí)段所占的面積。

2 結(jié)果及分析

2.1 農(nóng)作物空間分布精度評價(jià)

農(nóng)作物空間分布精度評價(jià)使用野外調(diào)查的驗(yàn)證點(diǎn)作為評價(jià)標(biāo)準(zhǔn)，將驗(yàn)證點(diǎn)與耕地矢量解譯圖斑進(jìn)行空間疊加，計(jì)算落在各作物地塊內(nèi)的點(diǎn)數(shù)，進(jìn)行精度評價(jià)。2018年共采集地面驗(yàn)證點(diǎn)991個(gè)，其中水稻219個(gè)、玉米283個(gè)、大豆465個(gè)、其他作物24個(gè)，結(jié)果顯示水稻、玉米、大豆和其他作物定性精度分別為97.26%、95.41%、96.99%和95.83%，總體精度為96.57%。

2.2 農(nóng)作物種植比例及面積時(shí)空變化分析

圖2 2014—2018年農(nóng)作物面積及比例Fig.2 Map of crop area and planting scale from 2014 to 2018

表1 2014—2018年作物面積統(tǒng)計(jì)Table 1 Statistical of crop area from 2014 to 2018 hm2

據(jù)圖2，表1可知，2014年作物種植比例（水稻∶玉米∶大豆∶其他作物）為14∶40∶45∶1，旱田占耕地面積的86.64%，水稻占13.36%；2015年作物種植比例為15∶54∶30∶1，旱田占耕地面積的85.08%，水稻占14.92%；2016年作物種植比例為17∶38∶43∶2，旱田占耕地面積的82.48%，水稻占17.52%；2017年作物種植比例為20∶21∶54∶5，旱田占耕地面積的80.34%，水稻占19.66%；2018年作物種植比例為19∶38∶38∶5，旱田占耕地面積的80.75%，水稻占19.25%。

表2 2014—2018年作物年際變化Table 2 Interannual variation of crops from 2014 to 2018

由表2可知，4個(gè)時(shí)段中，水稻在前3個(gè)時(shí)段的種植面積持續(xù)增加，2017—2018年減少，4個(gè)時(shí)段面積變化幅度分別為11.66%、17.89%、12.13%和-3.31%；玉米2014—2015年種植面積增加，2015—2017年面積持續(xù)減少，2017—2018年面積增加，4個(gè)時(shí)段面積變化幅度分別為34.21%、-29.65%、-44.00%和78.42%；大豆2014—2015年種植面積減少，2015—2017年面積持續(xù)增加，2017—2018年面積減少，4個(gè)時(shí)段面積變化幅度分別為-34.02%、44.95%、25.74%和-31.88%；其他作物種植面積在4個(gè)時(shí)段持續(xù)增加，且面積變化幅度分別為7.58%、19.71%、205.12%和4.38%。

2.3 作物變化方向時(shí)空特征分析

由圖3a可知，2014—2015、2015—2016、2016—2017和2017—2018年4時(shí)段水稻不變區(qū)域的面積分別為45 011.92、50 541.44、59 446.14、63 805.81 hm2，分別占各時(shí)段初期水稻面積的98.90%、99.46%、99.22%、94.98；增加區(qū)域的面積分別為5 805.82、9 369.26、7 734.54、1 269.67 hm2；減少區(qū)域的面積分別為498.54、276.30、464.57、3 374.87 hm2。

圖3 2014—2018年作物4個(gè)時(shí)段面積變化：a.水稻；b.玉米；c.大豆；d.其他作物Fig.3 Crop area variation in 4 periods from 2014 to 2018

由圖3b可知，4時(shí)段玉米不變區(qū)域的面積分別為85 100.10、84 221.16、39 059.24、33 271.83 hm2，分別占各時(shí)段初期玉米面積的62.24%、45.89%、30.25%、46.02%；增加區(qū)域的面積分別為98 415.50、44 882.12、33 240.39、95 728.38 hm2；減少區(qū)域的面積分別為 51 636.89、99 294.44、90 044.04、39 027.80 hm2。

由圖3c可知，4時(shí)段大豆不變區(qū)域的面積分別為54 187.09、54 747.05、103 527.78、84 931.86 hm2，分別占各時(shí)段初期大豆面積的35.14%、53.81%、70.20%、45.80%；增加區(qū)域的面積分別為47 553.20、92 726.55、81 907.17、41 391.72 hm2；減少區(qū)域的面積分別為 100 013.14、46 993.24、43 945.83、100 503.08 hm2。

由圖3d可知，4時(shí)段其他作物不變區(qū)域的面積分別為231.02、433.60、791.81、3 137.08 hm2，分別占各時(shí)段初期其他作物面積的5.38%、9.39%、14.32%、18.59%；增加區(qū)域的面積分別為4 388.60、5 096.41、16 081.43、14 475.45 hm2；減少區(qū)域的面積分別為 4 063.15、4 186.04、4 738.19、13 736.16 hm2。

2.4 旱田輪作方式時(shí)空變化分析

表3為作物輪作面積統(tǒng)計(jì)結(jié)果。

表3 年際間作物輪作面積統(tǒng)計(jì)Table 3 Statistical of interannual crop rotation area hm2

結(jié)果顯示，2014—2015年，旱田公共區(qū)面積為288 982.23 hm2，輪作區(qū)面積為149 464.02 hm2，占年際間旱田公共區(qū)面積的51.72%。其中旱田輪作方式玉米—大豆面積為45 956.50 hm2，占輪作區(qū)域面積的30.75%；玉米—其他作物面積為1 609.58 hm2，占輪作區(qū)域面積的1.08%；大豆—玉米面積為95 169.88 hm2，占輪作區(qū)域面積的63.67%；大豆—其他作物面積為2 775.36 hm2，占輪作區(qū)域面積的1.85%；其他作物—玉米面積為2 463.48 hm2，占輪作區(qū)域面積的1.65%；其他作物—大豆面積為1 489.22 hm2，占輪作區(qū)域面積的1.00%。

2015—2016年，旱田公共區(qū)面積為280 843.19 hm2，輪作區(qū)面積為141 441.41 hm2，占年際間旱田公共區(qū)面積的50.36%。其中旱田輪作方式玉米—大豆面積為90 372.91 hm2，占輪作區(qū)域面積的63.89%；玉米—其他作物面積為3 118.42 hm2，占輪作區(qū)域面積的2.21%；大豆—玉米面積為41 934.65 hm2，占輪作區(qū)域面積的29.65%；大豆—其他作物面積為1 897.56 hm2，占輪作區(qū)域面積的1.34%；其他作物—玉米面積為2 108.85 hm2，占輪作區(qū)域面積的1.49%；其他作物—大豆面積為2 009.01 hm2，占輪作區(qū)域面積的1.42%。

2016—2017年，旱田公共區(qū)面積為273 429.48 hm2，輪作區(qū)面積為130 075.51 hm2，占年際間旱田公共區(qū)面積的47.57%。其中旱田輪作方式玉米—大豆面積為78 218.54 hm2，占輪作區(qū)域面積的60.13%；玉米—其他作物面積為7 648.36 hm2，占輪作區(qū)域面積的5.88%；大豆—玉米面積為31 411.76 hm2，占輪作區(qū)域面積的24.15%；大豆—其他作物面積為8 178.62 hm2，占輪作區(qū)域面積的6.29%；其他作物—玉米面積為1 440.00 hm2，占輪作區(qū)域面積的1.11%；其他作物—大豆面積為3 178.23 hm2，占輪作區(qū)域面積的2.44%。

2017—2018年，旱田公共區(qū)面積為271 578.64 hm2，輪作區(qū)面積為150 241.42 hm2，占年際間旱田公共區(qū)面積的55.32%。其中旱田輪作方式玉米—大豆面積為33 711.87 hm2，占輪作區(qū)域面積的22.44%；玉米—其他作物面積為4 531.63 hm2，占輪作區(qū)域面積的3.02%；大豆—玉米面積為89 128.55 hm2，占輪作區(qū)域面積的59.32%；大豆—其他作物面積為9 716.00 hm2，占輪作區(qū)域面積的6.47%；其他作物—玉米面積為6 144.19 hm2，占輪作區(qū)域面積的4.09%；其他作物—大豆面積為7 009.18 hm2，占輪作區(qū)域面積的4.66%。

2.5 旱田作物輪作次數(shù)分析

圖4 2014—2018年輪作次數(shù)空間分布圖Fig.4 Map of crops rotation number from 2014 to 2018

表4 2014—2018年各輪作情況面積統(tǒng)計(jì)Table 4 Statistical data of area for each crop rotation period from 2014 to 2018

續(xù)表4

由表4可知4個(gè)時(shí)段旱田公共區(qū)作物輪作次數(shù)分別為0次、1次、2次、3次、4次的旱田面積分別為28 248.60、48 695.32、96 812.60、69 538.50、25 731.23 hm2，分別占旱田公共區(qū)面積的10.50%、18.10%、35.99%、25.85%、9.56%。

其中輪作0次和4次的只有1種輪作模式。輪作1次的旱田區(qū)域中，包含4種輪作模式，即僅2014—2015年輪作的面積為15 174.05 hm2，占1次輪作面積的31.16%；僅2015—2016年時(shí)段輪作的面積為10 690.97 hm2，占1次輪作面積的21.95；僅2016—2017年輪作的面積為8 056.74 hm2，占1次輪作面積的16.55%；僅2017—2018年輪作的面積為14 773.56 hm2，占1次輪作面積的30.34%。

輪作2次的旱田區(qū)域中，包含6種輪作模式，即僅2014—2015年和2015—2016年輪作的面積為23 412.25 hm2，占2次輪作旱田面積的24.18%；僅2014—2015年和2016—2017年輪作的面積為10 717.45 hm2，占11.07%；僅2014—2015年和2017—2018年輪作面積為13 175.41 hm2，占13.61%；僅2015—2016年和201—2017年輪作面積為11 992.86 hm2，占12.39%；僅2015—2016年和2017—2018年輪作的面積為14 528.80 hm2，占15.01%；僅2016—2017年和2017—2018年輪作面積為22 985.83 hm2，占23.74%。

輪作3次的旱田區(qū)域中，包含4種輪作模式，即僅2014—2015年、2015—2016年和2016—2017年輪作的面積為11 770.54 hm2，占3次輪作面積的16.93；僅2014—2015年時(shí)段、2015—2016年和2017—2018年輪作的面積為21 348.45 hm2，占30.70%；僅2014—2015年、2016—2017年和2017—2018時(shí)段年輪作的面積為19 864.33 hm2，占28.57%；僅2015—2016年、2016—2017年和2017—2018年輪作的面積為16 555.19 hm2，占23.81%。

3 結(jié)論

利用中高分辨率衛(wèi)星數(shù)據(jù)（Landsat 5 TM、Landsat 8 OLI、Sentinel-2A）獲得的耕地內(nèi)部作物分布數(shù)據(jù)，運(yùn)用數(shù)理統(tǒng)計(jì)方法，分析海倫市2014—2018年5年間作物空間格局時(shí)空變化特征，主要結(jié)論如下。

（1）5年間海倫市作物空間格局發(fā)生了較大的變化，水稻總體呈增長趨勢，其種植面積和比例2017年之前持續(xù)增加，年增長率超過10%，2018年略有降低；玉米、大豆種植面積和比例總體呈下降趨勢，年際間波動(dòng)劇烈，且變化趨勢相反；其他作物種植面積和比例持續(xù)增加，種植比例提高近5倍。

（2）2014—2015年、2015—2016年、2016—2017年、2017—2018年4個(gè)時(shí)段中旱田輪作區(qū)域面積高于45%，2017—2018年最高，達(dá)到55.32%；玉米和大豆之間的輪作比例總體呈下降趨勢，但玉米和大豆之間的輪作為主要的輪作方式，高于輪作面積的80%；4個(gè)時(shí)段均輪作的區(qū)域面積不足旱田公共區(qū)面積的10%，玉米、大豆的重茬和迎茬仍較為普遍。

4 討論

（1）國家農(nóng)業(yè)供給側(cè)改革政策的實(shí)施，對松嫩地區(qū)的作物種植結(jié)構(gòu)產(chǎn)生了較大的影響，對種植結(jié)構(gòu)的調(diào)整起到了積極的推進(jìn)作用：旱田中玉米、大豆的種植比例降低，其他作物比例升高，水稻面積回落，旱田輪作比例呈增高趨勢。同時(shí)，伴隨著農(nóng)村土地30年不變政策和土地流轉(zhuǎn)政策的實(shí)施，城鎮(zhèn)化的不斷推進(jìn)，以及農(nóng)村機(jī)械化和規(guī)?；潭炔粩嗵岣?，種植模式正一步一步從小農(nóng)戶經(jīng)營向土地規(guī)模經(jīng)營的方式轉(zhuǎn)變，同時(shí)，我國大豆、玉米收益受進(jìn)出口價(jià)格、政策性補(bǔ)貼等因素的影響較大，農(nóng)民種植大都根據(jù)前一年的種植收益來決定下一年的種植作物，導(dǎo)致黑龍江省的玉米、大豆面積年際間波動(dòng)巨大，農(nóng)民收入易受沖擊，制定穩(wěn)定、長期、綜合的農(nóng)業(yè)政策，合理規(guī)劃作物空間格局，對于保障國家糧食安全、合理引導(dǎo)農(nóng)民種植，保護(hù)農(nóng)民利益、提高農(nóng)民收入、保護(hù)耕地可持續(xù)利用起到積極的作用。

（2）黑龍江省的輪作主要是玉米和大豆輪作，不利于耕地的可持續(xù)利用，目前，小麥?zhǔn)呛邶埥〉牡谌蠛堤镒魑?，但是面積僅為10.2萬hm2［14］且主要種植于嫩江、北安、五大連池以北，收益低，農(nóng)民不愿意種植。建議根據(jù)地形、積溫、降水等自然條件，選擇適宜的作物，與玉米、大豆建立科學(xué)合理的“三三制”輪作制度。