李帥旗，牟鳳云，王 力，劉 曦

(重慶交通大學(xué)河海學(xué)院，重慶 400074)

商丘市土地利用變化及其對糧食安全的影響分析

李帥旗，牟鳳云，王 力，劉 曦

(重慶交通大學(xué)河海學(xué)院，重慶 400074)

基于1990年和2004年的TM遙感影像數(shù)據(jù)，借助遙感圖像處理軟件ERDAS、ENVI和地理信息系統(tǒng)軟件ArcGIS進(jìn)行分析，得到商丘市兩個時相間的土地利用轉(zhuǎn)移矩陣;采用了土地利用動態(tài)度、土地利用開發(fā)度、土地利用耗減度、土地利用程度綜合指數(shù)等模型對土地利用狀況進(jìn)行定量分析評價，并通過Markov鏈對商丘市2018年、2032年的土地利用結(jié)構(gòu)進(jìn)行了預(yù)測;研究結(jié)果表明，14年間商丘市主城區(qū)土地利用格局發(fā)生了巨大變化，尤以耕地和建設(shè)用地變化顯著，期間耕地面積減少了208.38 km2，而建設(shè)用地增加了141.89 km2，最后，基于耕地壓力指數(shù)分析了耕地數(shù)量的變化對于該農(nóng)業(yè)大市的糧食安全的影響，商丘市的糧食安全問題日益嚴(yán)重，迫切需要保護(hù)耕地的數(shù)量和質(zhì)量。

土地利用;轉(zhuǎn)移矩陣;馬爾可夫;糧食安全

土地利用/覆被類型是自然因素與人為因素共同作用結(jié)果。對土地利用/覆被變化(Land use/land cover change，LUCC)過程的數(shù)量、結(jié)構(gòu)和景觀特性的動態(tài)描述研究日趨成熟，定量化成為主流［1］。將遙感(RS)和地理信息系統(tǒng)(GIS)技術(shù)引入LUCC研究領(lǐng)域，可以快速、準(zhǔn)確、及時地獲取大面積的土地利用/土地覆蓋方面的實(shí)時信息，對土地利用變化進(jìn)行動態(tài)監(jiān)測，為解決復(fù)雜的空間問題提供條件，為可持續(xù)發(fā)展研究提供技術(shù)支撐［2］。2012年11月，國務(wù)院批復(fù)《中原經(jīng)濟(jì)區(qū)規(guī)劃指導(dǎo)意見》將商丘市定位為國家“三化”(新型工業(yè)化、新型城鎮(zhèn)化、新型農(nóng)業(yè)現(xiàn)代化)協(xié)調(diào)發(fā)展先行區(qū)、國家現(xiàn)代農(nóng)業(yè)示范區(qū)、中原經(jīng)濟(jì)區(qū)承接產(chǎn)業(yè)轉(zhuǎn)移示范市及內(nèi)陸開放高地［3］。在此政策的引導(dǎo)下，商丘市工業(yè)化和城鎮(zhèn)化步伐必將加快，人口增長與快速城市化過程導(dǎo)致的土地資源減少之間的矛盾將日益突出。商丘市是農(nóng)業(yè)大市，是國家糧食主產(chǎn)核心示范區(qū)、豫東糧倉，因而對土地利用變化和糧食安全的研究具有很強(qiáng)的現(xiàn)實(shí)意義。

1 研究區(qū)概況

商丘市位于河南省最東部，地理坐標(biāo)為114°49'-116°39'E，33°43'-34°52'N，豫、魯、蘇、皖4省結(jié)合處。商丘市總面積10 704 km2，約占河南省總面積的6.4%。截至2011年底全市總?cè)丝?90萬，城鎮(zhèn)化水平31. 5%。商丘市是農(nóng)業(yè)大市，是國家糧食主產(chǎn)核心示范區(qū)、豫東糧倉，2011年糧食總產(chǎn)量608.19萬噸，在河南省18個地級市中排名第3位。商丘市轄2個區(qū)、6個縣、1個縣級市。隨著工業(yè)化和城鎮(zhèn)化的逐步推進(jìn)，商丘市越來越成為土地利用變化最頻繁的熱點(diǎn)區(qū)域之一。論文的研究區(qū)域主要指商丘的兩個市轄區(qū)即梁園區(qū)和睢陽區(qū)，這兩個區(qū)域是商丘市土地利用變化較快的區(qū)域，具有代表意義。

2 數(shù)據(jù)來源與處理方法

2.1 數(shù)據(jù)來源

(1)遙感圖像:研究采用的遙感圖像通過國際科學(xué)數(shù)據(jù)服務(wù)平臺和中國科學(xué)院對地觀測與數(shù)字地球科學(xué)中心下載獲得，圖幅編號是123-36，成像時間分別是1990年9月2日和2004年7月22日。兩景圖像云量少，質(zhì)量較好，能滿足應(yīng)用需求。

(2)非遙感影像數(shù)據(jù):主要包括1∶50 000商丘市行政區(qū)劃圖、商丘市地形圖、商丘市土地利用規(guī)劃圖等，以及1990年、2004年和其他相關(guān)年份的商丘市人口、耕地面積、糧食產(chǎn)量等統(tǒng)計(jì)資料。

2.2 數(shù)據(jù)處理

基于地形圖，對兩期遙感數(shù)據(jù)進(jìn)行了幾何校正和配準(zhǔn)，幾何校正的精度控制在1個像元內(nèi)。結(jié)合國家土地利用類型標(biāo)準(zhǔn)和商丘市土地資源實(shí)際利用情況，將商丘市土地利用類型劃分為4種，即林地、水體、耕地和建設(shè)用地。在遙感圖像處理軟件ERDAS中利用最大似然法(Maximum Likelihood)進(jìn)行監(jiān)督分類，通過聚類統(tǒng)計(jì)、去除分析、分類重編碼方法相結(jié)合對監(jiān)督分類后的圖像進(jìn)行處理，以提高分類的精度，得到兩個時期的遙感分類圖像。

3 土地利用變化狀況分析

3.1 土地利用的數(shù)量變化

對土地利用現(xiàn)狀進(jìn)行分析，可以知道各類土地的特點(diǎn)和土地的利用結(jié)構(gòu)，從而掌握土地的動態(tài)變化及可持續(xù)利用的趨勢，為土地合理開發(fā)利用、整治保護(hù)，為國家政策制定提供依據(jù)［4］。利用ArcGIS9.3的統(tǒng)計(jì)功能，得到商丘市1990年和2004年的土地利用狀況(表1)，進(jìn)而得到兩個時期各類型土地所占比例。

表1 1990-2004年商丘市土地利用情況 km2

可以看出，商丘市土地利用以耕地為主，但耕地面積變化較大，從1990年到2004年間耕地面積減少了208.38 km2，很大程度上是由于期間商丘市大力推動工業(yè)化和城鎮(zhèn)化建設(shè)所引起;其次面積較大的是建設(shè)用地，14年間建設(shè)用地大幅增加，增幅達(dá)73.99%;林地、水體的面積所占比例較小，雖然變化幅度很大，但是面積并沒有太大的增加。

3.2 土地利用轉(zhuǎn)換分析

土地利用轉(zhuǎn)移矩陣能夠定量說明類型間的相互轉(zhuǎn)化情況，揭示他們之間相互轉(zhuǎn)化的比例和速率等特征。轉(zhuǎn)移矩陣不僅可以反映研究初、末期土地利用類型結(jié)構(gòu)，而且還可以反映研究期內(nèi)各土地利用轉(zhuǎn)移變化情況，便于了解研究初期各土地利用類型的轉(zhuǎn)出方向和研究末期各土地利用類型的轉(zhuǎn)入來源與數(shù)量特征［5］。在ArcGIS9.3空間疊置分析和Envi4.5中的Change Detection Analysis(變化監(jiān)測分析)功能支持下，統(tǒng)計(jì)了兩個時期間4種利用類型的土地相互轉(zhuǎn)換情況，得到1990-2004年間土地利用轉(zhuǎn)移矩陣(表2)，從變化數(shù)量和類型上反映了土地利用變化的總體特征。

表2 1990-2004年商丘市各土地利用類型面積轉(zhuǎn)移矩陣 km2

(1)林地轉(zhuǎn)換分析。從表2中可以看出，林地的變化包括輸入和輸出兩部分。在輸出部分，14年間，共有16.74 km2的林地轉(zhuǎn)換成其他類型的用地，其中有0.27 km2轉(zhuǎn)換成水體，9.86 km2轉(zhuǎn)換成耕地，6.61 km2轉(zhuǎn)換成建設(shè)用地。在輸入部分，14年間共有50.76 km2的其他類型用地轉(zhuǎn)換成林地，即水體0.27 km2，耕地46.07 km2，建設(shè)用地4.42 km2。林地的輸入部分比輸出部分要高34.02 km2，即在14年間林地面積增加34.02 km2，入出顯示出不平衡的態(tài)勢，主要是由耕地轉(zhuǎn)換而來。

(2)水體轉(zhuǎn)換分析。表2中顯示，14年間水體的輸入面積明顯高于輸出面積，在此期間有大量的其他用地轉(zhuǎn)換成水體，這包括0.27 km2的林地，28.22 km2的耕地，7.47 km2的建設(shè)用地，共計(jì)35.96 km2。而水體轉(zhuǎn)換成其他類型共計(jì)只有3.72 km2，其中轉(zhuǎn)換成林地0.27 km2，轉(zhuǎn)換成耕地1.42 km2，轉(zhuǎn)換成建設(shè)用地2.03 km2。在這些增加的面積中，也主要是由耕地面積轉(zhuǎn)換而來。

(3)耕地轉(zhuǎn)換分析?？梢钥闯觯?4年間，耕地絕大部分用于輸出轉(zhuǎn)換為其他類型的用地，收支明顯不平衡。具體來說，輸出方面，耕地轉(zhuǎn)換為林地、水體、建設(shè)用地的面積分別為:46.07 km2，28.22 km2，212.75 km2，輸入方面，林地、水體、建設(shè)用地分別有9.86 km2，1.42 km2，67.61 km2轉(zhuǎn)換為耕地。輸出比輸入面積多208.15 km2(輸出面積287.04 km2，輸入面積78.89 km2)，約占1990年耕地面積的13.56%，耕地輸出主要轉(zhuǎn)換為建設(shè)用地。

(4)建設(shè)用地轉(zhuǎn)換分析。14年間共有221.39 km2其他利用類型的土地轉(zhuǎn)換為建設(shè)用地，建設(shè)用地轉(zhuǎn)換成其他土地面積79.5 km2，呈現(xiàn)出明顯的畸形轉(zhuǎn)換趨勢，建設(shè)用地主要轉(zhuǎn)換成耕地，約占全部輸出面積的85. 04%，其來源也主要是耕地，約占全部轉(zhuǎn)入的96.10%，但是兩種轉(zhuǎn)換性質(zhì)的量明顯不同?？傮w說來，1990-2004年商丘市土地利用類型相互轉(zhuǎn)化明顯。隨著城市迅速擴(kuò)張，城市化水平的不斷加快，商丘市區(qū)發(fā)展迅速，城市建設(shè)日新月異，各種建設(shè)項(xiàng)目占用大量耕地，導(dǎo)致耕地面積銳減，建設(shè)用地越來越多。其次，商丘市政府積極調(diào)動農(nóng)民開展平原綠化防沙治沙的積極性，堅(jiān)持植樹造林，使水土流失得到有效控制，植被覆蓋率不斷提高［6］。

4 土地利用動態(tài)變化分析

土地利用數(shù)量變化特征可以用土地利用變化動態(tài)度指數(shù)來描述，它對比較土地利用變化的區(qū)域差異和預(yù)測未來土地利用變化趨勢都具有積極的作用［7］。

4.1 土地利用動態(tài)度

包括單一類型動態(tài)度和綜合變化動態(tài)度。土地利用變化動態(tài)度方法能夠定量描述土地利用變化的幅度和速度，對預(yù)測未來土地利用變化趨勢有積極的作用［8］。單一土地利用變化動態(tài)度反映的是區(qū)域一定時間間隔內(nèi)某種土地利用類型的數(shù)量變化情況，可以通過公式(1)計(jì)算:

綜合變化動態(tài)度反映的是一個地區(qū)土地利用數(shù)量變化的綜合情況，可用式(2)計(jì)算:

式(1)(2)中，Ki和K分別代表某一單一土地利用類型變化的動態(tài)度和區(qū)域綜合土地利用變化的動態(tài)度;Li，T1為研究初期i類型土地利用類型的面積;ΔLi，T1為研究期內(nèi)i類型土地利用類型未發(fā)生變化的面積;T1，T2分別代表研究初期和末期［9］。將相關(guān)數(shù)據(jù)帶入上述公式，可得到不同土地利用類型的單一動態(tài)度與綜合土地利用動態(tài)度(表3)。

表3 1990-2004年商丘市土地利用動態(tài)變化指數(shù) %

由表3可以看出，1990-2004年間，水體的變化速率最大，耕地變化速率最小，各土地利用類型變化的速度依次為水體＞林地＞建設(shè)用地＞耕地。

4.2 土地利用開發(fā)度

土地利用開發(fā)度是指單位時間內(nèi)某類型土地利用實(shí)際新開發(fā)的程度，以單位時間內(nèi)新開發(fā)的該類型土地面積占初期該類型總面積的百分比來表示［4］:

式(3)中，Dab為從a時刻到b時刻新開發(fā)某類型土地利用的面積，即由其他土地利用類型轉(zhuǎn)變?yōu)樵擃愋偷拿娣e總和;Ua表示a時刻該土地利用類型的面積;T為研究時長，以年為單位。

根據(jù)公式(3)計(jì)算出1990-2004年商丘市各種土地利用開發(fā)度，其中林地19．20%，水體23．85%，耕地0．37%，建設(shè)用地8．25%。顯然，耕地的開發(fā)度最小，這是因?yàn)樵诖似陂g耕地面積大量轉(zhuǎn)換為其他類型的土地，并沒有太多其他類型的土地新開發(fā)為耕地。

4.3 土地利用耗減度

為了測量單位時間內(nèi)某類型土地利用被實(shí)際消耗的程度，以單位時間內(nèi)已耗減的該類型土地面積占初期該類型總面積的百分比來表示［10］:式(4)中，Cab為從a時刻到b時刻某種土地利用類型被消耗的面積，即該種土地利用類型轉(zhuǎn)變?yōu)槠渌恋乩妙愋偷拿娣e總和;Ua表示a時刻該土地利用類型的面積;T為研究長，以年為單位。

計(jì)算可得，林地、水體、耕地、建設(shè)用地4種土地利用類型的耗減度分別為6．33%，2．47%，1．34%，2．96%。可以看出，14年間，林地每年轉(zhuǎn)換為其他類型的用地的速率最快，耗減度最大，這與大量的城市建設(shè)密切相關(guān)。

4.4 土地利用程度綜合指數(shù)和土地利用程度變化

土地利用程度及其變化可以定量表達(dá)該地區(qū)土地利用的綜合水平和變化趨勢。將土地利用的理想狀態(tài)定為4種土地利用級，并對4種土地利用級賦予其本身類別的值，則可得到4種土地利用程度的分級指數(shù)［11］(表4)。

表4 區(qū)域土地利用程度分級表

表4中的4種土地利用級僅僅是4種理想型，在實(shí)際狀態(tài)下，這4種類型都是混合存在于同一區(qū)域，各自占據(jù)不同的面積比例，并對當(dāng)?shù)赝恋乩贸潭?，按其自身的?quán)重做出自己的貢獻(xiàn)。據(jù)此，土地利用程度的綜合量化指標(biāo)必須在此基礎(chǔ)上進(jìn)行數(shù)學(xué)綜合，形成一個在1～4之間連續(xù)分布的綜合指數(shù)，其值的大小則綜合反映了某一區(qū)域土地利用程度。

土地利用程度綜合指數(shù)為

式中La為土地利用程度綜合指數(shù)，Ai為第i級的土地利用程度分級指數(shù)，Ci為第i級土地利用程度分級面積百分比，n為土地利用程度分級數(shù)。土地利用程度綜合量化指標(biāo)是一個從100～400之間連續(xù)變化的指標(biāo)，綜合指數(shù)的大小反映了土地利用程度的高低。

土地利用程度變化指數(shù)

若R為正值，則該區(qū)域土地利用處于發(fā)展期，否則處于衰退期。

代入公式(6)計(jì)算可得，1990年商丘市土地利用程度綜合指數(shù)為309．23;相應(yīng)地，2004年土地利用程度綜合指數(shù)為313．53，顯然R為正值，這表明商丘市土地利用正處于發(fā)展期，且開發(fā)利用強(qiáng)度越來越大。

5 未來土地利用預(yù)測

馬爾可夫(Markov)預(yù)測法，是一種預(yù)測事件發(fā)生概率的方法，它是基于馬爾可夫鏈，根據(jù)事件的目前狀況預(yù)測其將來各個時刻變動狀況的一種預(yù)測方法［12］。馬爾可夫過程是一種無后效性的隨機(jī)過程，而馬爾可夫鏈預(yù)測是根據(jù)變量目前狀態(tài)預(yù)測未來如何變化的分析方法，所以它不需要連續(xù)的歷史資料，而僅需要最近或現(xiàn)在的動態(tài)資料，便可預(yù)測未來。由于土地利用類型間的轉(zhuǎn)化滿足馬爾可夫鏈預(yù)測理論的特點(diǎn)，故馬爾可夫鏈理論可用于土地利用動態(tài)過程模擬與預(yù)測［13］。

5.1 馬爾可夫預(yù)測模型

根據(jù)馬爾可夫過程的無后效性及Bayes條件概率公式，有P(n)=P(n-1)P0或P(n)=P(0)Pn，P(n)為經(jīng)過n-1次狀態(tài)轉(zhuǎn)移到達(dá)n次轉(zhuǎn)移后的狀態(tài)概率向量;P(n-1)為經(jīng)過n-2次狀態(tài)轉(zhuǎn)移到達(dá)n-1次轉(zhuǎn)移后的狀態(tài)概率向量;P或Pij為轉(zhuǎn)移概率矩陣;P(0)為它的初始狀態(tài)概率向量。

5.2 土地利用轉(zhuǎn)移概率矩陣

根據(jù)前述研究結(jié)果，帶入公式(7)可得到1990-2004年的土地利用轉(zhuǎn)移概率矩陣(表5)。

表5 土地利用轉(zhuǎn)移概率矩陣

5.3 結(jié)果預(yù)測

借助Matlab7.0進(jìn)行編程，將1990年各種土地利用類型的面積和轉(zhuǎn)移概率矩陣作為初始值進(jìn)行運(yùn)算，可以得到對2018年和2032年土地利用結(jié)構(gòu)的預(yù)測結(jié)果(Y2018，Y2032分別代表2018年和2032年土地利用結(jié)果預(yù)測):

Y2018=［54．1，70．1，1222，410．8］，即2018年林地、水體、耕地、建設(shè)用地面積預(yù)測值分別為54．1 km2，70．1 km2，1 222 km2，410．8 km2。

Y2032=［53．5，86．9，1170．8，445．7］，即2032年林地、水體、耕地、建設(shè)用地面積預(yù)測值分別為53．5 km2，86．9 km2，1 170．8 km2，445．7 km2。

預(yù)測結(jié)果表明，在未來一段時間內(nèi)耕地數(shù)量將進(jìn)一步減少，但減少幅度將逐漸下降;建設(shè)用地進(jìn)一步增加，增加的速度將放緩。

6 土地利用變化對糧食安全影響分析

近年來，隨著經(jīng)濟(jì)技術(shù)的進(jìn)步，城鎮(zhèn)化步伐的加快，未利用土地大量開發(fā)，耕地面積迅速減少，嚴(yán)重威脅到維持人們生存的農(nóng)業(yè)安全［14］。耕地資源是農(nóng)民和農(nóng)業(yè)賴以生存和發(fā)展的基礎(chǔ)條件，耕地數(shù)量和質(zhì)量的變化直接影響糧食產(chǎn)量，從而影響到糧食有效供給和糧食安全水平，所以研究耕地變化對糧食產(chǎn)糧的影響至關(guān)緊要［15］。糧食是人類生存的最基本要素，土地是糧食生產(chǎn)的最終源泉。保證一定的耕地和基本的糧食供給，不僅關(guān)系個人的生存問題，而且直接關(guān)系國家或集體的生存安全。商丘市作為典型的農(nóng)業(yè)大市，是國家糧食主產(chǎn)核心示范區(qū)、豫東糧倉，對地區(qū)甚至全省、全國的糧食安全都起著舉足輕重的作用。

從耕地總量來看，商丘市耕地數(shù)量呈減少趨勢，耕地面積占總面積比例不斷減少。從耕地結(jié)構(gòu)來看，由于近幾年農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整，耕地內(nèi)部結(jié)構(gòu)也發(fā)生了巨大的變化。從人均來看，商丘市人均耕地面積呈快速減少趨勢。

耕地壓力指數(shù)可以衡量一個地區(qū)耕地資源的緊張程度，給出耕地保護(hù)的閾值，可作為耕地保護(hù)的調(diào)控基準(zhǔn)［13］:

式中，K為耕地壓力指數(shù)，Smin為最小人均耕地面積(hm2/人);Sa為實(shí)際人均耕地面積(是區(qū)域耕地總面積與人口數(shù)量的函數(shù))。K=1是耕地保護(hù)的警戒線，此時必須防止耕地流失，在提高耕地的物質(zhì)投入水平和生產(chǎn)能力的情況下合理引導(dǎo)耕地的用途轉(zhuǎn)移。當(dāng)K＜1時，表明耕地資源有所富余，在這種情況下，可以適度轉(zhuǎn)換耕地用途以保證城市化、工業(yè)化對土地的需求，也可以大力調(diào)整農(nóng)業(yè)種植結(jié)構(gòu)，退耕還林還草以改善生態(tài)環(huán)境，甚至可以讓耕地休閑，以保持耕地的綜合生產(chǎn)能力。當(dāng)K＞1時，耕地資源承受著巨大的壓力。商丘市隸屬于我國中部地區(qū)，研究表明，我國中部地區(qū)最小耕地面積為0．117 8hm2/人［16］，依據(jù)近幾年的統(tǒng)計(jì)數(shù)據(jù)，計(jì)算可得近幾年商丘市的耕地壓力指數(shù)(圖1)。

由圖1可以看出，自進(jìn)入2000年以來，商丘市耕地壓力指數(shù)一直大于1，且呈現(xiàn)逐年上升的趨勢。這表明耕地資源承受著巨大的壓力，危及糧食安全的資源基礎(chǔ)。此時，不僅要限制耕地資源向非農(nóng)利用轉(zhuǎn)移，減少耕地被占用，更要增加物質(zhì)投入和科技投入來提高土地生產(chǎn)力，以確保地區(qū)耕地資源可持續(xù)利用和糧食安全。

圖1 2002-2008年商丘市耕地壓力指數(shù)

7 結(jié)論與討論

在RS和GIS技術(shù)的支持下，對商丘市1990年和2004年TM影像進(jìn)行處理和分析，得到其土地利用變化狀況并對2018年和2032年的土地利用結(jié)構(gòu)進(jìn)行了預(yù)測，得出以下結(jié)論:

(1)14年間商丘市主城區(qū)土地利用格局發(fā)生了巨大變化，尤以耕地和建設(shè)用地變化顯著，期間耕地面積減少了208.38 km2，而建設(shè)用地增加了141.89 km2，土地利用強(qiáng)度逐漸增大。

(2)運(yùn)用馬爾可夫轉(zhuǎn)移鏈對未來商丘市土地利用結(jié)構(gòu)進(jìn)行預(yù)測，結(jié)果表明商丘市耕地面積在未來一定階段仍將繼續(xù)減少，建設(shè)用地將占據(jù)大量耕地和其他類型的土地。

(3)耕地壓力指數(shù)表明，商丘市的糧食安全問題日益嚴(yán)重，保護(hù)耕地數(shù)量和質(zhì)量任務(wù)艱巨，迫在眉睫。

對1990-2004年間商丘市土地利用變化進(jìn)行了較為透徹的分析，總結(jié)出該市近年來土地利用的時空變化特點(diǎn)并對未來的土地利用進(jìn)行了預(yù)測，希望能為該市土地資源合理利用與保護(hù)提供理論基礎(chǔ)，為有關(guān)管理部門制定土地資源開發(fā)的政策法規(guī)提供較可靠的科學(xué)依據(jù)。但是，也存在許多不足之處，如論文僅對商丘市土地利用的變化進(jìn)行分析，沒有涉及與之緊密聯(lián)系的驅(qū)動力研究，還需進(jìn)一步開展工作。此外，論文所采用的TM影像的空間分辨率不夠、精度受限，提取的用地類型偏少;解譯遙感圖像采用的是監(jiān)督分類，不能消除誤差，這些都是在今后的學(xué)習(xí)和科研中更加努力重視和改進(jìn)的問題。

參考文獻(xiàn):

［1］陳穎彪，李雁，千慶蘭，等.基于TM數(shù)據(jù)的廣州市番禺區(qū)土地覆被格局分析［J］.地理科學(xué)，2012，32(4):458-463

［2］周楊，周廷剛，馮沖.基于RS和GIS的LUCC研究——以重慶市萬州區(qū)為例［J］.長江科學(xué)院院報，2012，29(5):67-71

［3］國務(wù)院關(guān)于支持河南省加快建設(shè)中原經(jīng)濟(jì)區(qū)的指導(dǎo)意見［OL］.(2011-10-07)［2013-04-09］.http://www.gov.cn/zwgk/2011-10/07/content_1963574.htm

［4］MENG J J，WU X Q，LI Z G.Land Use/Land Cover Changes in Zhangye Oasis of Hexi Corridor［J］.Journal of Geographical Sciences，2003，13(1):71-75

［5］朱會義，李秀彬.關(guān)于區(qū)域土地利用變化指數(shù)模型方法的討論［J］.地理學(xué)報，2003，58(5):643-650

［6］賀振.基于TM影像的商丘市區(qū)土地利用變化動態(tài)監(jiān)測［J］.資源與產(chǎn)業(yè)，2010，12(5):91-95

［7］王秀蘭，包玉海.土地利用動態(tài)變化研究方法探討［J］.地理科學(xué)進(jìn)展，1999，18(1):81-87

［8］李忠鋒，王彥麗.定邊縣土地利用動態(tài)變化研究［J］.干旱區(qū)地理，2004，27(4):520-524

［9］朱會義，李秀彬，何書金，等.環(huán)渤海地區(qū)土地利用的時空變化分析［J］.地理學(xué)報，2001，56(3):253-259

［10］張新長，梁金成.城市土地利用動態(tài)變化及預(yù)測模型研究［J］.中山大學(xué)學(xué)報:自然科學(xué)版，2004，43(2):121-125

［11］莊大方，劉紀(jì)遠(yuǎn).中國土地利用程度的區(qū)域分異模型研究［J］.自然資源學(xué)報，1997，12(2):105-111

［12］徐建華.計(jì)量地理學(xué)［M］.北京:高等教育出版社，2006

［13］何瑞珍，張敬東，胡艷紅，等.基于RS和GIS的鄭州市土地利用動態(tài)變化及預(yù)測［J］.甘肅農(nóng)業(yè)大學(xué)學(xué)報，2009(4): 126-131

［14］范玉枝.山東省土地利用變化對農(nóng)業(yè)安全的影響研究［J］.小城鎮(zhèn)建設(shè)，2004(3):62-63

［15］魯奇.中國耕地資源開發(fā)、保護(hù)與糧食安全保障問題［J］.資源科學(xué)，1999，14(4):5-8

［16］蔡運(yùn)龍，傅澤強(qiáng)，戴爾阜.區(qū)域最小人均耕地面積與耕地資源調(diào)控［J］.地理學(xué)報，2002，57(2):127-134

Analysis of the Change of Land Exploitation of Shanqiu City and Its Influence on Grain Safety

LI Shuai-qi，MOU Feng-yun，WANG Li，LIU Xi
(School of River and Sea，Chongqing Jiaotong University，Chongqing 400074，China)

Based on TM remote sensing image data in 1990 and in 2004，by the analysis based on remote sensing image processing software ERDAS，ENVI and geographic information system software ArcGIS，the land exploitation transfer matrix of Shanqiu City at two phases is obtained.By using the models such as dynamic degree of land exploitation，developing degree of land exploitation，depletion degree of land exploitation，comprehensive index of land exploitation degree，and so on，quantitative analysis and evaluation are conducted on the status quo of the land exploitation.By using Markov chain，the land exploitation structure of Shanqiu City is predicted in 2018 and in 2032.The results show that the land exploitation pattern of main urban area of Shanqiu City has greatly taken place in the past 14 years，especially the change of arable land and construction land is significant，in this period，arable land area reduced 208.38 km2while construction land increased 141.89 km2.Finally，the influence of the change of the number of arable land on grain safety of this big agricultural city is analyzed based on arable land pressure index and the result demonstrates that the grain problem in Shanqiu City is more and more serious and so the number of arable land and the quality arable land are urgently needing to be protected.

land exploitation;transfer matrix;Markov;grain safety

F301.24

A

1672-058X(2014)01-0074-08

責(zé)任編輯:代小紅

校 對:李翠薇

2013-06-23;

2013-10-10.

李帥旗(1988-)，男，河南商丘人，碩士研究生，從事區(qū)域分析與規(guī)劃方面的研究.