基于GIS的黃河北干流岸線土地利用變化分析

梁之豪 時(shí)芳欣 陳融旭 穆磊

摘 要：河岸帶是水陸交錯(cuò)帶，資源十分寶貴，其質(zhì)量關(guān)系到陸地及河流兩者的生態(tài)狀況，黃河流域岸線資源更為寶貴，為摸清其土地利用時(shí)空變化規(guī)律，以黃河北干流岸線為研究區(qū)域，采用GIS手段分析了黃河大北干流岸線利用需求較為迫切的河曲—準(zhǔn)格爾旗（河準(zhǔn)段）與府谷—保德（府保段）兩段河道和小北干流岸線土地利用時(shí)空變化情況，然后結(jié)合土地利用變化轉(zhuǎn)移矩陣，利用馬爾柯夫模型分別對(duì)其土地利用結(jié)構(gòu)進(jìn)行預(yù)測(cè)。分析可知，北干流岸線多年間土地利用在耕地、水域和建筑用地之間轉(zhuǎn)化相對(duì)劇烈，其中大部分區(qū)域耕地、建筑用地呈增長(zhǎng)態(tài)勢(shì)，水域呈減少態(tài)勢(shì)。原因可能是隨著北干流岸線區(qū)域經(jīng)濟(jì)發(fā)展，工業(yè)化、城鎮(zhèn)化進(jìn)程加快，建筑用地、耕地需求增加。利用馬爾柯夫模型對(duì)北干流岸線土地利用結(jié)構(gòu)的預(yù)測(cè)結(jié)果表明，隨著經(jīng)濟(jì)社會(huì)的持續(xù)發(fā)展，沿河人民拓寬生存空間的需求會(huì)持續(xù)增長(zhǎng)，北干流岸線將逐漸形成以耕地、建筑用地和水域?yàn)橹鞯牟季中问健?/p>

關(guān)鍵詞：北干流岸線;土地利用變化;GIS;轉(zhuǎn)移矩陣;馬爾柯夫預(yù)測(cè)

中圖分類號(hào)：TV212.4;TV882.1 文獻(xiàn)標(biāo)志碼：A

doi：10.3969/j.issn.1000-1379.2021.10.024

引用格式：梁之豪，時(shí)芳欣，陳融旭，等.基于GIS的黃河北干流岸線土地利用變化分析[J].人民黃河，2021，43（10）：123-127.

Abstract： Riparian zone is a waterway crisscross zone， and its resources are very precious. Its quality is related to the ecological conditions of land and rivers. The shoreline resources of the Yellow River Basin are more precious. In order to find out the temporal and spatial changes of land use， this paper takes the shoreline of the north main stream of the Yellow River as the research area， and uses GIS to analyze the temporal and spatial changes of land use of the shoreline of the great north main stream and the small north main stream of the Yellow River Basin， and then combines the transfer matrix of land use change and Markov model to predict the land quantitative structure. The analysis shows that the land use in the north main stream coastline has changed dramatically between cultivated land， water area and construction land for many years， and the cultivated land and construction land in most areas are increasing， while the water area is decreasing. This may be because with the economic development of the coastal area along the north main stream， the process of industrialization and urbanization is speeding up， and the demand for construction land and cultivated land is increasing. Using Markov model to predict the quantitative structure of land use along the north main stream coastline， it shows that with the sustainable development of economy and society， the demand of people along the river to broaden living space will continue to grow， and the coastline of the north main stream will gradually form a layout of mainly cultivated land， construction land and water area.

Key words： north main stream coastline; land use change; GIS; transfer matrix; markov forecast

河岸帶是水陸交錯(cuò)帶，資源十分寶貴，其質(zhì)量關(guān)系到陸地及河流兩者的生態(tài)狀況。黃河流域河道岸線是有限的寶貴資源，其獨(dú)特的地理位置使其具有豐富的自然功能和社會(huì)價(jià)值，已成為沿河經(jīng)濟(jì)社會(huì)建設(shè)的重要基礎(chǔ)和寶貴戰(zhàn)略資源。對(duì)其進(jìn)行合理開發(fā)利用和保護(hù)治理是黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展的重點(diǎn)工作之一。土地利用變化是反映人類活動(dòng)程度的重要因子，分析土地利用時(shí)空變化規(guī)律是揭示人類活動(dòng)程度的有效方式[1]。

本文采用GIS手段對(duì)黃河流域大北干流岸線的河準(zhǔn)段（河曲—準(zhǔn)格爾旗）、府保段（府谷—保德）和小北干流岸線進(jìn)行土地利用結(jié)構(gòu)變化、土地利用轉(zhuǎn)移矩陣和土地利用空間變化情況分析，得到2010—2017年大北干流岸線和1990—2015年小北干流岸線土地利用時(shí)空變化情況，然后結(jié)合土地利用轉(zhuǎn)移矩陣，利用馬爾柯夫模型對(duì)大北干流岸線2017—2037年和小北干流岸線2015—2035年土地利用結(jié)構(gòu)進(jìn)行預(yù)測(cè)。

1 研究區(qū)概況

黃河干流的河口鎮(zhèn)至禹門口段也稱大北干流河段，流域面積11.16萬km2，河長(zhǎng)725.1 km，水力資源豐富;禹門口至潼關(guān)段稱為小北干流河段，該河段有汾河、渭河兩大支流匯入，河谷展寬，河長(zhǎng)約130 km。

大北干流區(qū)域內(nèi)諸多地帶并無岸線開發(fā)需求，本文選用人口數(shù)量相對(duì)較多、岸線利用需求較為迫切的河準(zhǔn)段和府保段對(duì)岸線資源進(jìn)行分析。河準(zhǔn)段河道屬山區(qū)峽谷型河道，是黃河著名的彎曲性河道之一，段內(nèi)農(nóng)業(yè)經(jīng)濟(jì)較發(fā)達(dá)，涉河建筑較多，岸線相對(duì)穩(wěn)定，岸線邊界采用歷史最高洪水位854.7 m的淹沒范圍邊界。府保段河道也屬山區(qū)峽谷型河道，河灘地主要是耕地、林地和少量園地，岸線相對(duì)穩(wěn)定，岸線邊界采用歷史最高洪水位820.2 m的淹沒范圍邊界。

小北干流河道左岸為山西省運(yùn)城市，右岸為陜西省渭南市，屬淤積性游蕩型河道，劇烈的河勢(shì)變化使得沿河群眾經(jīng)常受到洪水威脅，岸線邊界采用20 a一遇淹沒范圍邊界。

2 研究方法及數(shù)據(jù)來源

2.1 研究方法

（1）單一土地利用動(dòng)態(tài)度。單一土地利用動(dòng)態(tài)度簡(jiǎn)單來說就是指一定時(shí)間內(nèi)某種土地利用類型的變化率，它對(duì)比較土地利用變化的區(qū)域差異和預(yù)測(cè)未來土地利用變化趨勢(shì)都有積極的作用[2]，公式為

（2）土地利用轉(zhuǎn)移矩陣。采用土地利用轉(zhuǎn)移矩陣可以衡量土地利用類型的動(dòng)態(tài)變化過程，其來源于系統(tǒng)分析中對(duì)系統(tǒng)狀態(tài)與狀態(tài)轉(zhuǎn)移的定量描述，較為常用[3]的轉(zhuǎn)移矩陣為

式中：n為土地利用類型數(shù);Pij代表研究時(shí)段內(nèi)第i類轉(zhuǎn)換為第j類土地利用的面積或百分比。在實(shí)際應(yīng)用中，通常將矩陣用表格形式來表示。

（3）土地利用變化預(yù)測(cè)。馬爾柯夫預(yù)測(cè)法是由俄國(guó)數(shù)學(xué)家馬爾柯夫創(chuàng)造并命名的一種方法，利用轉(zhuǎn)移矩陣分析某一系統(tǒng)的現(xiàn)狀及其發(fā)展動(dòng)向，預(yù)測(cè)該系統(tǒng)未來狀況。利用這種方法可以揭示不同土地利用類型之間的轉(zhuǎn)移概率，更好地分析其時(shí)空演變過程[4-12]。從某一現(xiàn)狀出發(fā)，假定在其發(fā)展過程中有n個(gè)可能的狀態(tài)，即K1、K2、…、Kn。記Sij為從狀態(tài)Ki轉(zhuǎn)變?yōu)闋顟B(tài)Kj的狀態(tài)轉(zhuǎn)移概率，則矩陣為

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

研究采用的小北干流數(shù)據(jù)為1990年、1995年、2000年、2005年、2010年、2015年的1 km柵格數(shù)據(jù)，采用的土地利用數(shù)據(jù)由國(guó)家科技基礎(chǔ)條件平臺(tái)——國(guó)家地球系統(tǒng)科學(xué)數(shù)據(jù)中心（http：//www.geodata.cn）提供。大北干流研究區(qū)區(qū)域較小，數(shù)據(jù)有限，使用的數(shù)據(jù)為清華大學(xué)提供的2010年和2017年的30 m柵格數(shù)據(jù)。數(shù)據(jù)集均通過遙感影像人機(jī)交互解譯，結(jié)合外業(yè)考察和驗(yàn)證，通過柵格化處理獲得，數(shù)據(jù)真實(shí)可靠。

3 研究區(qū)土地利用變化分析

3.1 時(shí)間變化分析

2010年和2017年大北干流河準(zhǔn)段及府保段岸線土地利用情況見圖1～圖4，1990年、2015年小北干流岸線土地利用情況見圖5、圖6。

分別統(tǒng)計(jì)2010—2017年大北干流河準(zhǔn)段及府保段岸線和1990—2015年小北干流岸線土地利用類型變化情況，見表1。

分析表1可知：大北干流河準(zhǔn)段岸線從2010年到2017年土地利用類型變化明顯，其中建筑用地和水域變化顯著，建筑用地呈現(xiàn)上升的變化趨勢(shì)最為明顯，水域呈現(xiàn)下降的變化趨勢(shì)較明顯，耕地、林地、草地、未利用地相對(duì)變化較小。原因可能是隨著岸線區(qū)域工業(yè)化、城鎮(zhèn)化的快速發(fā)展，建筑用地需求增大，岸線側(cè)城市為盡力發(fā)展工業(yè)區(qū)、基礎(chǔ)建設(shè)等出現(xiàn)了城區(qū)改造及圍河造地現(xiàn)象。

府保段岸線從2010年到2017年土地利用類型變化也比較明顯，其中耕地、林地、水域下降趨勢(shì)明顯，建筑用地上升趨勢(shì)最為明顯，草地、未利用地變化相對(duì)較小。原因可能是隨著工業(yè)化、城鎮(zhèn)化快速發(fā)展，大量農(nóng)村人口尤其是青壯年勞動(dòng)力外出務(wù)工，岸線區(qū)域常住人口減少導(dǎo)致對(duì)耕地的需求逐漸萎縮。而同時(shí)，為了快速發(fā)展城市的基礎(chǔ)設(shè)施及工業(yè)，對(duì)建筑用地產(chǎn)生了巨大的需求，沿河居民將岸線區(qū)域的大量土地劃為建筑用地，甚至出現(xiàn)自發(fā)性圍河造地現(xiàn)象。

小北干流1990—2015年岸線土地利用類型變化較大，其中耕地面積呈現(xiàn)先升后降的趨勢(shì)，原因可能是在2010年之前為了滿足岸線經(jīng)濟(jì)發(fā)展及農(nóng)村居民對(duì)耕地需求的增長(zhǎng)，將岸線區(qū)域的大量土地劃為耕地，以及出現(xiàn)圍河造地現(xiàn)象，而在2010年以后隨著農(nóng)村人口尤其是青壯年勞動(dòng)力外出務(wù)工數(shù)量增多，耕地面積需求減小。水域面積呈現(xiàn)先降后升的變化規(guī)律，這與耕地面積變化規(guī)律正好相反，反映出耕地與水域面積變化之間的關(guān)聯(lián)性。建筑用地、未利用地面積整體變化幅度不大。

3.2 空間變化分析

采用ArcGIS空間分析工具，對(duì)2010—2017年大北干流河準(zhǔn)段及府保段、1990—2015年小北干流岸線土地利用空間分布變化進(jìn)行分析，結(jié)果見圖7～圖9。

由圖7可知：岸線土地利用類型在河準(zhǔn)段北部臨近河道處出現(xiàn)了小范圍變化，2010—2017年水域和建筑用地的空間轉(zhuǎn)化相對(duì)劇烈。其中北部出現(xiàn)了較多水域轉(zhuǎn)化為建筑用地的情況，原因是岸線區(qū)域工業(yè)化、城鎮(zhèn)化快速發(fā)展，建筑用地需求增加。

分析圖8可知：岸線土地利用類型在府保段臨近河道處出現(xiàn)了小范圍變化，2010—2017年耕地、水域、建筑用地的空間轉(zhuǎn)化相對(duì)劇烈。其中在府保段北部、府谷和保德兩縣縣城交界的區(qū)域出現(xiàn)了大量水域、耕地轉(zhuǎn)化為建筑用地的情況，說明青壯年勞動(dòng)力外出，對(duì)耕地的需求減小。而同時(shí)，城市建設(shè)和發(fā)展對(duì)建筑土地產(chǎn)生了巨大的需求，由此導(dǎo)致了侵占水域面積的現(xiàn)象。

分析圖9可知：1990—2015年小北干流岸線土地利用變化主要集中在北部及中部，南部變化較小，其中耕地與水域的變化最為明顯。該河段為山區(qū)峽谷型河段，灘地本身難以固定，但是隨著經(jīng)濟(jì)社會(huì)發(fā)展、人口的增加及上游大型水利樞紐的投入使用，導(dǎo)致水情發(fā)生巨大變化，為了擴(kuò)寬生存空間及改善自身生活環(huán)境，部分峽谷區(qū)的群眾搬遷至河灘上，當(dāng)?shù)卣疄楸Ｗo(hù)群眾生命財(cái)產(chǎn)安全修建了灘岸工程，隨著群眾圍墾土地逐年增多，部分灘地已變成了基本農(nóng)田。

3.3 土地利用轉(zhuǎn)移變化

根據(jù)式（2）計(jì)算的大北干流河準(zhǔn)段及府保段、小北干流岸線土地利用轉(zhuǎn)移矩陣見表2～表4。

分析可知：2000—2017年河準(zhǔn)段土地類型變化較大的為耕地、水域和建筑用地，其余土地類型變化較小。其中水域有3.14 km2演變成耕地、1.37 km2演化為建筑用地，同時(shí)期內(nèi)有1.43 km2耕地演化為建筑用地，有1.05 km2耕地演化為水域。因此，在水域內(nèi)占用河道進(jìn)行耕種和圍河造地進(jìn)行城市建設(shè)可能是這段時(shí)間水域面積大幅度減小的首要原因。

府保段土地類型變化較大的為耕地、水域和建筑用地，其余土地類型變化較小。其中，耕地有5.50 km2和2.32 km2演變?yōu)榻ㄖ玫睾筒莸兀瑫r(shí)期內(nèi)只有0.53 km2建筑用地演化為耕地，農(nóng)村勞動(dòng)力外流與占用耕種面積進(jìn)行城市建設(shè)可能是主要原因;水域面積減小，有2.89 km2演變?yōu)榻ㄖ玫亍?.61 km2演變成耕地，同時(shí)期建筑用地僅有0.09 km2演變?yōu)樗?，由此可見，在該期間存在人為侵占岸線資源的情況。

大北干流河段部分岸線土地利用類型在耕地、水域與建設(shè)用地之間反復(fù)轉(zhuǎn)換，這是近年來隨著大北干流沿岸地區(qū)的經(jīng)濟(jì)發(fā)展，人口遷移變化和經(jīng)濟(jì)發(fā)展造成城鎮(zhèn)化進(jìn)程加快等原因造成的。

分析表4可知：1990—2015年小北干流土地類型變化較大的為耕地、草地和水域，其余土地類型變化較小。其中水域有144 km2演變成耕地，同時(shí)期內(nèi)只有30 km2耕地演化為水域。因此，該時(shí)期存在較為突出的與水爭(zhēng)地現(xiàn)象，有46 km2水域演變?yōu)椴莸兀莸赜?9 km2和42 km2演變?yōu)樗蚝透?，耕地僅有6 km2演變?yōu)椴莸?。小北干流區(qū)域?yàn)榈湫偷暮拥烙问巺^(qū)，上游水利樞紐的泥沙調(diào)控使黃河的部分水流沖出原有河道，導(dǎo)致部分耕地和草地被水域覆蓋，同時(shí)也使部分河道變?yōu)闉┑兀芙?jīng)濟(jì)發(fā)展和人口增長(zhǎng)等因素的驅(qū)動(dòng)，這些區(qū)域存在“占灘為耕”的現(xiàn)象，導(dǎo)致部分水域和草地面積減小，轉(zhuǎn)為耕地。

3.4 土地利用變化預(yù)測(cè)

將2017年大北干流河準(zhǔn)段、府保段各土地類型的面積設(shè)置為初始狀態(tài)向量，每經(jīng)過1 a為1步，利用轉(zhuǎn)移概率矩陣，預(yù)測(cè)2037年土地利用結(jié)構(gòu)，結(jié)果見表5、表6。

大北干流河準(zhǔn)段、府保段未來土地利用變化趨勢(shì)形成的主要原因可能有：①截至2017年大北干流河準(zhǔn)段共修建河道工程14處，府保段共修建河道工程6處，隨著這些水利樞紐工程的使用，加之上游用水量增加等，使得河道干流水量減少，從而形成大量的沿河灘地;②隨著區(qū)域經(jīng)濟(jì)社會(huì)的持續(xù)發(fā)展，人們對(duì)生存空間拓展的需求激增，部分群眾為改善生活條件搬遷至河岸地區(qū)，將河灘陸續(xù)改造成耕地、建筑用地等;③河段沿岸礦產(chǎn)、旅游等資源豐富，如磧口、壺口瀑布等成為“網(wǎng)紅打卡地”，各地政府因勢(shì)利導(dǎo)，大力投資，在當(dāng)前發(fā)展趨勢(shì)下，如果沒有有效的政策加以調(diào)控，圍河造地現(xiàn)象會(huì)日趨嚴(yán)重。

將2015年小北干流各土地類型的面積設(shè)置為初始狀態(tài)向量，每經(jīng)過1 a為1步，預(yù)測(cè)2035年土地利用結(jié)構(gòu)，結(jié)果見表7。

分析可知：未來小北干流的土地利用結(jié)構(gòu)變化最劇烈的是耕地和水域，這也說明了未來隨著經(jīng)濟(jì)社會(huì)的不斷發(fā)展，與水爭(zhēng)地現(xiàn)象會(huì)日趨嚴(yán)重。小北干流這種發(fā)展趨勢(shì)的形成可能有以下原因：①隨著城鎮(zhèn)化進(jìn)程的加快，城市人口急劇增多，加之沿河旅游公路等交通條件改善，群眾對(duì)生存空間拓展的愿望日趨強(qiáng)烈;②當(dāng)?shù)卣疄榇_保沿河居民生命財(cái)產(chǎn)安全，陸續(xù)修建了一些護(hù)岸、護(hù)村工程，灘地位置逐漸穩(wěn)定，一些灘地被群眾自發(fā)性開墾成基本農(nóng)田;③與大北干流情況類似，為提升河岸帶居民生活幸福感，一些旅游開發(fā)項(xiàng)目、城市景觀等工程陸續(xù)修建。

為了確保黃河流域水生態(tài)高質(zhì)量發(fā)展，黃河河岸帶健康發(fā)展是今后的重點(diǎn)，因此有以下建議：①沿河城鎮(zhèn)應(yīng)實(shí)施最嚴(yán)格的河道管理制度和水資源管理制度;②定時(shí)進(jìn)村、進(jìn)社區(qū)開展法制宣傳教育活動(dòng)，將保護(hù)黃河的理念植入民心;③將沿河居民對(duì)改善生活環(huán)境的愿景與河岸治理規(guī)劃結(jié)合起來，努力做到“雙贏”。

4 結(jié) 論

本文以黃河北干流岸線為研究區(qū)域，利用GIS等技術(shù)手段，對(duì)2010—2017年大北干流河準(zhǔn)段及府保段岸線、1990—2015年小北干流岸線土地利用變化進(jìn)行了分析，得出以下結(jié)論。

（1）黃河北干流岸線土地利用總體結(jié)構(gòu)上呈現(xiàn)建筑用地增長(zhǎng)、水域面積減少的態(tài)勢(shì)，可能原因是隨著北干流岸線區(qū)域經(jīng)濟(jì)發(fā)展，工業(yè)化、城鎮(zhèn)化進(jìn)程加快，沿黃岸線城市出于拓寬城鎮(zhèn)居民生存空間的考量，對(duì)位于岸線內(nèi)的主城區(qū)基礎(chǔ)設(shè)施、建筑用地等進(jìn)行了改造利用，并且出現(xiàn)了圍河造地的現(xiàn)象。

（2）大北干流府保段耕地面積出現(xiàn)下降趨勢(shì)，小北干流的耕地面積以2010年為轉(zhuǎn)折點(diǎn)呈現(xiàn)先增長(zhǎng)后減少的態(tài)勢(shì)，原因是隨著城鎮(zhèn)化進(jìn)程的加快，大量農(nóng)村青壯年勞動(dòng)力外出務(wù)工，在一些村莊形成了“空巢村”的現(xiàn)象，農(nóng)村人口對(duì)耕地需求量減少，部分耕地隨著城鎮(zhèn)化被改造成了建筑用地，這與2008年后國(guó)內(nèi)城鎮(zhèn)化進(jìn)程速度加快相吻合。

（3）結(jié)合土地利用變化轉(zhuǎn)移矩陣，利用馬爾柯夫模型對(duì)大北干流岸線2037年和小北干流岸線2035年土地利用結(jié)構(gòu)進(jìn)行預(yù)測(cè)。結(jié)果發(fā)現(xiàn)：隨著經(jīng)濟(jì)社會(huì)的持續(xù)發(fā)展，建筑用地的需求量會(huì)持續(xù)增大，水域面積會(huì)持續(xù)減小;大北干流河準(zhǔn)段和小北干流岸線耕地資源需求會(huì)持續(xù)增大，大北干流府保段耕地資源需求會(huì)減小。由此可以推斷：在土地政策沒有突然變化的情況下，位于岸線內(nèi)的城市化建設(shè)會(huì)導(dǎo)致對(duì)建筑用地的需求繼續(xù)增大，為了拓寬發(fā)展空間，還會(huì)出現(xiàn)圍河造地現(xiàn)象，大北干流府保段農(nóng)村人口向城鎮(zhèn)轉(zhuǎn)移會(huì)日趨加劇，岸線內(nèi)居民對(duì)耕地的需求會(huì)逐漸減小;北干流岸線土地資源將逐漸形成以耕地、建筑用地和水域?yàn)橹鞯牟季中问健?/p>

參考文獻(xiàn)：

[1] 吳琳娜，楊勝天，劉曉燕，等.1976年以來北洛河流域土地利用變化對(duì)人類活動(dòng)程度的響應(yīng)[J].地理學(xué)報(bào)，2014，69（1）：54-63.

[2] 梁龍禛，岳彩榮.基于GIS的石家莊市土地利用變化分析及預(yù)測(cè)[J].林業(yè)調(diào)查規(guī)劃，2016，41（1）：1-5，10.

[3] 朱會(huì)義，李秀彬.關(guān)于區(qū)域土地利用變化指數(shù)模型方法的討論[J].地理學(xué)報(bào)，2003，58（5）：643-650.

[4] 田傳召，于陽(yáng)，肖虹雁，等.基于TM影像的昆明市區(qū)土地利用遙感動(dòng)態(tài)監(jiān)測(cè)[J].林業(yè)資源管理，2014（4）：103-108.

[5] 林秀芝，董晨燕，蘇林山，等.黃河小北干流沖淤與水沙響應(yīng)關(guān)系分析[J].人民黃河，2019，41（5）：5-8.

[6] 門昭宇，侯志軍，侯嬌建，等.黃河大北干流河段采砂量分析[J].人民黃河，2019，41（7）：57-60.

[7] 梁如心，李娟，張維江，等.基于GIS的寧夏六盤山地區(qū)干旱風(fēng)險(xiǎn)分析[J].人民黃河，2020，42（1）：48-52.

[8] 郭彥君，郭文炯.晉西黃土丘陵區(qū)土地利用變化及生態(tài)效應(yīng)分析[J].人民黃河，2021，43（2）：106-111.

[9] 馮永玖，韓震.基于遙感的黃浦江沿岸土地利用時(shí)空演化特征分析[J].國(guó)土資源遙感，2010（2）：91-96.

[10] 賀勝柏.馬爾柯夫預(yù)測(cè)模型及其實(shí)證分析[J].中南財(cái)經(jīng)政法大學(xué)學(xué)報(bào)，2002（2）：132-135.

[11] 張前進(jìn)，李笑笑.土地利用結(jié)構(gòu)優(yōu)化配置方法研究綜述[J].天津城市建設(shè)學(xué)院學(xué)報(bào)，2013，19（1）：10-15，33.

[12] 劉榮霞，薛安，韓鵬，等.土地利用結(jié)構(gòu)優(yōu)化方法述評(píng)[J].北京大學(xué)學(xué)報(bào)（自然科學(xué)版），2005，41（4）：655-662.

【責(zé)任編輯 張華巖】