生態(tài)移民工程與生態(tài)系統(tǒng)可持續(xù)發(fā)展的系統(tǒng)動(dòng)力學(xué)研究

鐘水映　馮英杰

摘要 生態(tài)移民工程是近年來(lái)移民和生態(tài)環(huán)境研究方面的熱點(diǎn)問(wèn)題之一，它是在生態(tài)脆弱性地區(qū)面臨“生態(tài)”和“貧困”雙重短板下，基于人口、資源環(huán)境與經(jīng)濟(jì)社會(huì)的協(xié)調(diào)可持續(xù)發(fā)展原則而形成的一種由政府主導(dǎo)的系統(tǒng)的保護(hù)生態(tài)和扶貧的專項(xiàng)工程，著重實(shí)現(xiàn)生態(tài)脆弱性區(qū)域脫貧致富和生態(tài)保護(hù)，而實(shí)現(xiàn)生態(tài)移民的可持續(xù)安置和發(fā)展的基本前提是生態(tài)移民工程和生態(tài)系統(tǒng)可持續(xù)發(fā)展。通過(guò)構(gòu)建生態(tài)移民工程中的生態(tài)可持續(xù)系統(tǒng)動(dòng)力學(xué)模型，揭示了生境破碎化、生態(tài)移民工程和生態(tài)福利績(jī)效三個(gè)子系統(tǒng)之間的因果關(guān)系和動(dòng)態(tài)反饋機(jī)制，并以三江源地區(qū)生態(tài)移民為例對(duì)移民前后遷入地和遷出地的生態(tài)系統(tǒng)的可持續(xù)性進(jìn)行仿真模擬。對(duì)模型進(jìn)行的結(jié)構(gòu)適合性和行為一致性檢驗(yàn)表明，模型真實(shí)有效地反映了生態(tài)移民工程中生態(tài)系統(tǒng)的實(shí)際系統(tǒng)結(jié)構(gòu)和功能。最后，通過(guò)調(diào)整生境破碎度指數(shù)、生態(tài)移民環(huán)境保護(hù)投資、生態(tài)恢復(fù)因子和生態(tài)發(fā)展因子四個(gè)政策變量的取值，從而設(shè)定了六種模式以進(jìn)行生態(tài)移民生態(tài)可持續(xù)發(fā)展政策的模擬。政策模擬的結(jié)果顯示：政策模式1～6反映了生態(tài)移民環(huán)境保護(hù)投資的增加，有助于降低生態(tài)移民區(qū)域生境破碎指數(shù)和提高生態(tài)福利績(jī)效；政策模式3～6反映了生態(tài)移民環(huán)境保護(hù)投資應(yīng)側(cè)重于移民遷入地的生態(tài)發(fā)展投資，適當(dāng)降低移民遷出地的生態(tài)修復(fù)投資比例；政策模式5～6說(shuō)明了生態(tài)移民環(huán)境保護(hù)投資運(yùn)作效率，影響生態(tài)移民生態(tài)系統(tǒng)可持續(xù)發(fā)展。因此，加大生態(tài)移民環(huán)境保護(hù)投資力度，調(diào)整遷出地和遷入地的生態(tài)投資比重對(duì)于實(shí)現(xiàn)生態(tài)移民生態(tài)系統(tǒng)重建和可持續(xù)發(fā)展至關(guān)重要，在此基礎(chǔ)上提高生態(tài)投資運(yùn)作效率是促進(jìn)生態(tài)移民工程與生態(tài)系統(tǒng)可持續(xù)協(xié)調(diào)發(fā)展的最佳政策模式。

關(guān)鍵詞 系統(tǒng)動(dòng)力學(xué)；生境破碎化；生態(tài)移民工程；生態(tài)福利績(jī)效

中圖分類號(hào) F062.2 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 1002-2104（2018）11-0010-10 DOI：10.12062/cpre.20180516

隨著中國(guó)進(jìn)入資源消耗與財(cái)富集聚的加速更迭消長(zhǎng)階段，區(qū)域發(fā)展不平衡不充分問(wèn)題愈發(fā)突出，尤其是生態(tài)等因素長(zhǎng)期制約著區(qū)域間的平衡與協(xié)調(diào)發(fā)展。一方面，由于人為因素或者自然因素而導(dǎo)致完整的生態(tài)系統(tǒng)不斷被分割破碎或生態(tài)系統(tǒng)服務(wù)功能降低，生境破碎化已經(jīng)成為生態(tài)系統(tǒng)空間演化的重要特征，直接影響生態(tài)系統(tǒng)的可持續(xù)發(fā)展；另一方面，生態(tài)環(huán)境脆弱性往往與貧困現(xiàn)象陷入一種互為因果的惡性循環(huán)，生態(tài)移民是改善生態(tài)環(huán)境和消除區(qū)域性貧困而不得不實(shí)施的戰(zhàn)略工程，生態(tài)移民群體已經(jīng)成為除水電等工程移民之外的第二大移民群體。生態(tài)移民可以被理解為是一種調(diào)整生態(tài)系統(tǒng)服務(wù)供給與需求的空間不匹配的一種措施，通過(guò)改變?nèi)藗兝蒙鷳B(tài)系統(tǒng)服務(wù)的位置或變換生態(tài)系統(tǒng)服務(wù)供給的地點(diǎn)來(lái)實(shí)現(xiàn)人與自然關(guān)系的重新調(diào)整。生態(tài)移民過(guò)程中，面臨遷出地的自然生態(tài)環(huán)境問(wèn)題和遷入地的發(fā)展性生態(tài)環(huán)境問(wèn)題，如何實(shí)現(xiàn)生態(tài)移民涉及區(qū)域生態(tài)系統(tǒng)的重建與可持續(xù)發(fā)展，進(jìn)行生境破碎化研究與生態(tài)移民相關(guān)活動(dòng)的分析顯得十分必要。

1 文獻(xiàn)綜述

生境破碎化主要表現(xiàn)在環(huán)境因素自變和人為因素它變引起的形態(tài)破碎和功能破碎兩個(gè)方面，其所產(chǎn)生的面積、隔離、邊緣效應(yīng)通過(guò)作用于生物和非生物條件影響生態(tài)系統(tǒng)可持續(xù)發(fā)展，學(xué)者們分別從耕地和森林景觀等生境面積的喪失和城市等生境空間格局的改變展開了研究[1-5]。武正軍等[6]提出應(yīng)確定物種的適宜生境比例、物種分布和絕滅閾值，以研究異質(zhì)種群動(dòng)態(tài)、生境破碎化對(duì)物種絕滅的相對(duì)影響以及生態(tài)系統(tǒng)的生態(tài)過(guò)程。張?zhí)煨碌萚7]利用地塊數(shù)量與面積相關(guān)指數(shù)、通道密度，結(jié)合Google地圖數(shù)據(jù)，構(gòu)建了城市邊緣區(qū)空間破碎化指標(biāo)，論證了城市空間破碎化程度過(guò)高直接制約著城市空間的可持續(xù)發(fā)展。定量研究生境格局破碎化是景觀生態(tài)學(xué)的主要內(nèi)容，常用的景觀格局分析和模擬軟件為Fragstats和Apack[8]。在研究生境破碎對(duì)動(dòng)物物種群動(dòng)態(tài)的影響時(shí)，生態(tài)學(xué)家常用顯空間模型模擬動(dòng)物個(gè)體在景觀碎片間的活動(dòng)（如遷移等），并借助于遙感、全球定位系統(tǒng)、地理信息系統(tǒng)（3S）和GIS等技術(shù)研究生境破碎化對(duì)種群動(dòng)態(tài)活動(dòng)的影響[9-12]。

生態(tài)移民是指因環(huán)境脆弱而不具備生產(chǎn)生活所必需的條件或減少人類活動(dòng)對(duì)重要生態(tài)功能區(qū)的生態(tài)環(huán)境壓力而引發(fā)的人口遷移活動(dòng)[13-14]，同時(shí)，生態(tài)移民不是簡(jiǎn)單的人口生存空間轉(zhuǎn)移，更是對(duì)移民傳統(tǒng)的生產(chǎn)生活方式的重大轉(zhuǎn)型[15]。即，生態(tài)移民可持續(xù)安置與發(fā)展既

是人與生態(tài)系統(tǒng)關(guān)系調(diào)整的過(guò)程，又是包含區(qū)域經(jīng)濟(jì)系統(tǒng)與社會(huì)系統(tǒng)的整體可持續(xù)發(fā)展與重構(gòu)過(guò)程。而目前關(guān)于生態(tài)移民的研究多以經(jīng)濟(jì)社會(huì)問(wèn)題為核心訴求，體現(xiàn)出生存型移民的部分特征與功能，生態(tài)移民現(xiàn)象在生態(tài)與經(jīng)濟(jì)社會(huì)之間表現(xiàn)出一種倒置狀態(tài)[16]。許多學(xué)者從可持續(xù)發(fā)展的角度對(duì)生態(tài)移民問(wèn)題進(jìn)行了研究。李笑春等[17-18]

指出了保護(hù)西部地區(qū)生態(tài)環(huán)境脆弱性的緊迫性，從人口壓力與生態(tài)環(huán)境承載力角度論證了實(shí)施生態(tài)移民的必要性。侯東民等[19-20]提出借助國(guó)家財(cái)政能力以提高生態(tài)移民效益，把生態(tài)移民與城鎮(zhèn)化、“退耕還林”“退牧還草”等政策結(jié)合以來(lái)，可以更好地促進(jìn)生態(tài)恢復(fù)與經(jīng)濟(jì)持續(xù)發(fā)展。劉學(xué)敏等[21-22]認(rèn)為西北地區(qū)人地矛盾是導(dǎo)致生態(tài)環(huán)境惡化與經(jīng)濟(jì)落后的主要原因，生態(tài)移民的實(shí)施可以改變農(nóng)牧區(qū)人口分布的不合理，通過(guò)實(shí)現(xiàn)農(nóng)牧區(qū)移民的城鎮(zhèn)化，可以改善西部地區(qū)生態(tài)環(huán)境壓力和促使產(chǎn)業(yè)結(jié)構(gòu)合理化。劉慧等[23]

認(rèn)為生態(tài)移民生態(tài)扶貧包含原地扶貧和離地扶貧，前者依托當(dāng)?shù)氐纳鷳B(tài)資源環(huán)境發(fā)展特色綠色產(chǎn)業(yè)體系，創(chuàng)造生態(tài)管護(hù)等就業(yè)崗位，后者依托城鎮(zhèn)化等戰(zhàn)略，引導(dǎo)移民向城鎮(zhèn)和東部沿海地區(qū)轉(zhuǎn)移，實(shí)現(xiàn)異地就業(yè)，并提出了人力資本開發(fā)、綠色優(yōu)勢(shì)產(chǎn)業(yè)建設(shè)、制度創(chuàng)新建設(shè)的政策建議。張建軍等[24]

通過(guò)梳理塔里木河流域生態(tài)、生態(tài)移民狀況，認(rèn)為遷入地有限的環(huán)境資源不能完全滿足生態(tài)移民的發(fā)展需要，應(yīng)探索以人口素質(zhì)提高等為主要內(nèi)容的非農(nóng)化轉(zhuǎn)移方式，從而實(shí)現(xiàn)人口遷移與環(huán)境資源配置的可持續(xù)發(fā)展。定量分析中，層次分析法是通過(guò)權(quán)衡費(fèi)用和效益以評(píng)估生態(tài)移民綜合效益的主要方法；效率雙重差分法是分析生態(tài)移民前后移民收入變化、移民和非移民農(nóng)戶差異、地區(qū)人口遷移效率差別的主要方法，“壓力—狀態(tài)—相應(yīng)”（PSR）模型是評(píng)估資源可持續(xù)發(fā)展的主要方法[25]。

一些學(xué)者運(yùn)用系統(tǒng)動(dòng)力學(xué)方法對(duì)生態(tài)可持續(xù)發(fā)展（如資源環(huán)境耦合協(xié)調(diào)發(fā)展模式、生態(tài)承載力預(yù)警、生態(tài)服務(wù)功能價(jià)值仿真、生態(tài)安全評(píng)估）進(jìn)行了研究。李柏山[26]

基于距離指數(shù)法和GIS技術(shù)構(gòu)建了漢江流域水資源承載力系統(tǒng)動(dòng)力學(xué)模型，提出了提高水資源承載力的最優(yōu)方案，并分析了流域水資源承載力、水電梯級(jí)開發(fā)對(duì)漢江流域生態(tài)環(huán)境的影響及主要影響因子，從而對(duì)漢江流域的生態(tài)系統(tǒng)健康進(jìn)行評(píng)價(jià)。薛建春[27]

通過(guò)構(gòu)建生態(tài)足跡系統(tǒng)動(dòng)力學(xué)模型，分析了20年間平朔礦區(qū)生態(tài)系統(tǒng)的生態(tài)赤字狀況，并借助生態(tài)足跡指數(shù)、效率、多樣性等探索了生態(tài)系統(tǒng)的管理方式，提倡建立生態(tài)產(chǎn)業(yè)鏈?zhǔn)降陌l(fā)展模式。李華[28]

運(yùn)用系統(tǒng)動(dòng)力學(xué)方法對(duì)崇明生態(tài)系統(tǒng)的發(fā)展趨勢(shì)和生態(tài)風(fēng)險(xiǎn)進(jìn)行仿真，構(gòu)建并評(píng)價(jià)了河口島嶼的生態(tài)安全及空間分異性，在揭示生態(tài)系統(tǒng)作用機(jī)制的基礎(chǔ)上，提出生態(tài)管理和保護(hù)措施。張晶[29]利用產(chǎn)業(yè)生態(tài)系統(tǒng)的生態(tài)可承載、產(chǎn)業(yè)多樣性與主導(dǎo)型協(xié)調(diào)、生態(tài)效率高效三個(gè)定量化特征構(gòu)建了產(chǎn)業(yè)生態(tài)系統(tǒng)的健康評(píng)價(jià)系統(tǒng)動(dòng)力學(xué)模型，并探討了徐州市的環(huán)境壓力、產(chǎn)業(yè)經(jīng)濟(jì)、社會(huì)發(fā)展與可持續(xù)發(fā)展四個(gè)子系統(tǒng)的演化狀況。

綜上所述，當(dāng)前學(xué)者對(duì)于生態(tài)可持續(xù)發(fā)展的研究，涉及生境破碎化、生態(tài)移民中的生態(tài)修復(fù)與生態(tài)—經(jīng)濟(jì)—社會(huì)系統(tǒng)的協(xié)調(diào)發(fā)展、基于系統(tǒng)動(dòng)力學(xué)的對(duì)于生態(tài)系統(tǒng)健康安全評(píng)估和生態(tài)系統(tǒng)服務(wù)價(jià)值估算等多個(gè)方面，而缺少對(duì)生境破碎化、生態(tài)移民及生態(tài)福利績(jī)效三者相互影響機(jī)制的研究。本文試圖構(gòu)建生境破碎化、生態(tài)移民、生態(tài)福利績(jī)效的系統(tǒng)動(dòng)力學(xué)（SD）模型，并借助系統(tǒng)動(dòng)力學(xué)可視化建模工具 Vensim PLE 6.2版本軟件，對(duì)生態(tài)移民遷入地和遷出地的生態(tài)系統(tǒng)的可持續(xù)性進(jìn)行仿真模擬，深入系統(tǒng)地分析生態(tài)移民工程中生態(tài)可持續(xù)系統(tǒng)的功能、結(jié)構(gòu)與移民行為之間的關(guān)系，并系統(tǒng)仿真和政策模擬生態(tài)移民涉及區(qū)域的生態(tài)系統(tǒng)可持續(xù)發(fā)展的演化，以為生態(tài)移民政策的制定提供具有操作性的建議。

2 生態(tài)移民工程生態(tài)系統(tǒng)可持續(xù)系統(tǒng)動(dòng)力學(xué)模型構(gòu)建

2.1 系統(tǒng)邊界及子系統(tǒng)劃分

本文考慮的生態(tài)移民工程中生態(tài)可持續(xù)系統(tǒng)邊界包含由于實(shí)施生態(tài)移民工程所要實(shí)現(xiàn)的生態(tài)系統(tǒng)可持續(xù)的發(fā)展，即研究生態(tài)效果改善與生境破碎化、生態(tài)移民工程之間的關(guān)系，從而尋找促進(jìn)生態(tài)可持續(xù)發(fā)展的生態(tài)移民途徑。因此，生態(tài)可持續(xù)發(fā)展系統(tǒng)由生境破碎化子系統(tǒng)、生態(tài)移民子系統(tǒng)和生態(tài)福利績(jī)效子系統(tǒng)組成。需要說(shuō)明的是，經(jīng)濟(jì)社會(huì)系統(tǒng)必定會(huì)對(duì)生態(tài)可持續(xù)發(fā)展產(chǎn)生影響，但考慮到本研究中側(cè)重生態(tài)移民工程中生態(tài)系統(tǒng)的重建，因而沒有單獨(dú)設(shè)置經(jīng)濟(jì)和社會(huì)發(fā)展子系統(tǒng)，但是將受其影響的變量作為輸入變量，以表函數(shù)形式呈現(xiàn)，同時(shí)采用擬合方程的形式，使系統(tǒng)設(shè)定更貼近生態(tài)移民現(xiàn)實(shí)數(shù)據(jù)，能夠較好地反映生態(tài)移民的實(shí)際狀況。

生境破碎化子系統(tǒng)主要研究移民遷出地區(qū)人們以土地過(guò)度開墾、砍伐森林、污水排放等形式對(duì)生態(tài)環(huán)境的干擾，以及交通、水利、電力通信、教育、醫(yī)療衛(wèi)生條件的密度和差異性，并且研究了其與周邊產(chǎn)業(yè)基地（園區(qū)）、集（城）鎮(zhèn)、鄉(xiāng)村旅游景區(qū)和其他滿足水土資源平衡的生態(tài)盈余地區(qū)的距離，設(shè)定密度差異、干擾度和鄰近指數(shù)等變量；生態(tài)移民子系統(tǒng)主要研究生態(tài)移民工程如何通過(guò)人口密度重組、資源重配和結(jié)構(gòu)重調(diào)，實(shí)現(xiàn)移民與遷入地區(qū)的人口、資源和環(huán)境的可持續(xù)發(fā)展，表現(xiàn)在移民選址與安置、移民就業(yè)、移民社區(qū)各個(gè)方面，設(shè)定移民投資、移民數(shù)量、移民就業(yè)等變量；生態(tài)福利績(jī)效子系統(tǒng)主要研究遷出地區(qū)的生態(tài)修復(fù)與遷入地區(qū)發(fā)展生態(tài)產(chǎn)業(yè)如何通過(guò)調(diào)整人類發(fā)展指數(shù)和生態(tài)足跡以影響生態(tài)可持續(xù)發(fā)展，設(shè)定人類發(fā)展指數(shù)、生態(tài)足跡等變量。這三個(gè)子系統(tǒng)是相互影響與相互制約的，每個(gè)子系統(tǒng)的運(yùn)行既取決于子系統(tǒng)內(nèi)部的結(jié)構(gòu)，也取決于與其他子系統(tǒng)的關(guān)系，例如生境破碎化子系統(tǒng)直接影響生態(tài)能否可持續(xù)發(fā)展，但是生態(tài)移民工程的實(shí)施卻有助于改善生境破碎化，通過(guò)生態(tài)效果的逐步改善實(shí)現(xiàn)生態(tài)可持續(xù)發(fā)展。在相關(guān)文獻(xiàn)、調(diào)查研究和新聞報(bào)道的基礎(chǔ)上，以及前文關(guān)于關(guān)鍵要素的分析，構(gòu)建了生態(tài)移民工程中生態(tài)可持續(xù)發(fā)展的系統(tǒng)理論圖（見圖1）。

2.2 因果關(guān)系模型構(gòu)建

2.2.1 生境破碎化子系統(tǒng)因果關(guān)系分析

生態(tài)可持續(xù)的生境破碎化子系統(tǒng)主要是指生態(tài)移民遷出地的脆弱性生態(tài)系統(tǒng)。由于鄰近指數(shù)涉及生態(tài)移民遷出地與遷入地的實(shí)際距離和預(yù)期距離等不確定性因素，為了表征遷出地的脆弱性生態(tài)系統(tǒng)，僅考慮了生境破碎度量化中的關(guān)鍵因素，即以森林面積、草地面積、水土流失面積和耕地面積作為密度和差異指數(shù)，以人為砍伐面積和過(guò)度放牧（放牧牲畜數(shù)量）作為干擾度指數(shù)來(lái)表征生境破碎度。

生境破碎化子系統(tǒng)的主要反饋回路有：

森林面積→砍伐成耕地比例→砍伐成耕地面積→耕地面積→耕地收入→移民貧困→人為砍伐面積→森林面積；

草地面積→其他收入→移民貧困→過(guò)度放牧→草地面積；

水土流失面積→耕地面積→耕地收入→移民貧困→人為砍伐面積→森林面積→水土流失面積；

生境破碎化子系統(tǒng)的演化，始于森林面積、草地面積、水土流失面積和耕地面積等生境密度和差異，并受制于砍伐面積、過(guò)度放牧等人為因素的干擾。自然因素和人為因素的雙重作用，影響遷出地生態(tài)脆弱性的程度和區(qū)域面積，并影響了遷出地與遷入地的（具備優(yōu)越的自然條件）距離，從而影響了生態(tài)移民的選址與安置方式，移民的遷出反過(guò)來(lái)又緩解了生境破碎度的人地矛盾。

2.2.2 生態(tài)移民子系統(tǒng)因果關(guān)系分析

生態(tài)移民工程子系統(tǒng)是一個(gè)包含諸多因素的非線性系統(tǒng)，本文僅側(cè)重分析促使生態(tài)可持續(xù)的生態(tài)移民工程子系統(tǒng)，即生態(tài)移民工程中與生態(tài)有關(guān)并影響生態(tài)可持續(xù)的因素是分析的重點(diǎn)。因此，生態(tài)可持續(xù)的生態(tài)移民子系統(tǒng)主要是指生態(tài)移民工程中，通過(guò)實(shí)施遷出地生態(tài)恢復(fù)投資和遷入地生態(tài)發(fā)展投資，以解決遷出地的自然性環(huán)境問(wèn)題和遷入地的發(fā)展性環(huán)境問(wèn)題。在該系統(tǒng)中，重點(diǎn)考慮遷出地的水土治理投資、退耕還林（植樹造林）、退牧還草等生態(tài)恢復(fù)措施，遷入地移民發(fā)展中的生態(tài)休閑產(chǎn)業(yè)和現(xiàn)代服務(wù)業(yè)，以動(dòng)態(tài)反映生態(tài)移民工程中生態(tài)系統(tǒng)可持續(xù)的實(shí)現(xiàn)過(guò)程。

生態(tài)移民子系統(tǒng)的主要反饋回路有：

遷出地生態(tài)恢復(fù)投資→水土治理投資→遷出地的生態(tài)赤字→遷出地生態(tài)恢復(fù)投資；

遷出地生態(tài)恢復(fù)投資→植樹造林投資→遷出地的生態(tài)赤字→遷出地生態(tài)恢復(fù)投資；

遷出地生態(tài)恢復(fù)投資→退耕還林投資→遷出地的生態(tài)赤字→遷出地生態(tài)恢復(fù)投資；

遷出地生態(tài)恢復(fù)投資→退牧還草投資→遷出地的生態(tài)赤字→遷出地生態(tài)恢復(fù)投資；

遷入地生態(tài)發(fā)展投資→移民產(chǎn)業(yè)生態(tài)化效率→第一產(chǎn)業(yè)所占比重→產(chǎn)業(yè)結(jié)構(gòu)能耗系數(shù)→遷入地生態(tài)發(fā)展投資；

遷入地生態(tài)發(fā)展投資→移民產(chǎn)業(yè)生態(tài)化效率→第三產(chǎn)業(yè)所占比重→產(chǎn)業(yè)結(jié)構(gòu)能耗系數(shù)→遷入地生態(tài)發(fā)展投資。

2.2.3 生態(tài)福利績(jī)效子系統(tǒng)因果關(guān)系分析

生態(tài)可持續(xù)的生態(tài)福利績(jī)效子系統(tǒng)是影響生態(tài)可持續(xù)水平的重要子系統(tǒng)，生態(tài)福利績(jī)效是衡量自然消耗轉(zhuǎn)化為福利水平的效率的重要指標(biāo)，一方面，在生態(tài)環(huán)境的承載能力內(nèi)實(shí)現(xiàn)較高的福利水平是可持續(xù)發(fā)展的目標(biāo)[30]；另一方面，生態(tài)移民具有扶貧和生態(tài)修復(fù)的雙重目標(biāo)，以生態(tài)福利績(jī)效考核生態(tài)移民中的生態(tài)可持續(xù)發(fā)展，若生態(tài)福利績(jī)效提升，表明生態(tài)移民用較少的自然消耗（較多的生態(tài)修復(fù)）實(shí)現(xiàn)了更高比例的福利水平提升，生態(tài)正朝著可持續(xù)發(fā)展的目標(biāo)不斷靠近。因此，在生態(tài)福利績(jī)效子系統(tǒng)中，主要考慮了GDP、教育、預(yù)期壽命和生態(tài)足跡等因素。

生態(tài)福利績(jī)效子系統(tǒng)的主要反饋回路有：

GDP→財(cái)政支出→生態(tài)移民環(huán)境保護(hù)投資→GDP漏出→GDP；

GDP→財(cái)政支出→社區(qū)醫(yī)療投入→死亡率降低比例→預(yù)期壽命→移民收入→財(cái)政收入→GDP；

GDP→財(cái)政支出→教育投入→移民平均受教育年限→移民可就業(yè)比例→移民收入→財(cái)政收入→GDP；

教育投入→移民平均受教育年限→平均受教育年限對(duì)環(huán)保的支持度→環(huán)保態(tài)度→人均生態(tài)足跡→遷入地生態(tài)發(fā)展投資→生態(tài)移民環(huán)境保護(hù)投資→GDP漏出→教育投入；

預(yù)期壽命→移民收入→財(cái)政收入→GDP漏出→社區(qū)醫(yī)療投入→預(yù)期壽命；

人均生態(tài)足跡→生態(tài)盈余→環(huán)境質(zhì)量→生態(tài)移民環(huán)境保護(hù)投資→遷入地生態(tài)發(fā)展投資→人均生態(tài)足跡；

人均生態(tài)足跡→生態(tài)盈余→環(huán)境質(zhì)量→生態(tài)移民環(huán)境保護(hù)投資→遷出地生態(tài)恢復(fù)投資→人均生態(tài)足跡；

綜合生境破碎度子系統(tǒng)、生態(tài)移民子系統(tǒng)和生態(tài)福利績(jī)效子系統(tǒng)的因果關(guān)系分析，構(gòu)建出生態(tài)移民工程中生態(tài)可持續(xù)系統(tǒng)的因果關(guān)系模型。

3 系統(tǒng)動(dòng)力學(xué)總流圖模型構(gòu)建

通過(guò)對(duì)生態(tài)移民工程中生態(tài)可持續(xù)系統(tǒng)的因果關(guān)系圖的深入分析，根據(jù)各變量在系統(tǒng)中的不同作用，構(gòu)建出生態(tài)移民工程中生態(tài)可持續(xù)系統(tǒng)的流圖模型，如圖2所示。

3.1 系統(tǒng)變量性質(zhì)設(shè)置

在生態(tài)移民工程中生態(tài)可持續(xù)系統(tǒng)的流圖中，狀態(tài)變量包括森林面積等10個(gè)，速率變量包括森林面積增加量、森林面積減少量等17個(gè)，表函數(shù)為水土治理投資因子等17個(gè)，其他變量均為輔助變量，共74個(gè)，具體設(shè)置如下：

（1）生境破碎化子系統(tǒng)中，森林面積、草地面積、水土流失面積和耕地面積為狀態(tài)變量；森林面積增加量、森林面積減少量、草地面積增加量、草地面積減少量、水土流失面積增加量、水土流失面積減少量、耕地面積增加量、耕地面積減少量為速率變量；水土治理投資因子、植樹造林投資因子、退耕還林投資因子、退牧還草投資因子、人均生態(tài)承載能力、水土流失影響系數(shù)、森林面積影響系數(shù)、放牧因子、全國(guó)農(nóng)牧民純收入為表函數(shù)；其余為輔助變量。

（2）生態(tài)移民子系統(tǒng)中，遷出地生態(tài)恢復(fù)投資和遷入地生態(tài)發(fā)展投資為狀態(tài)變量；生態(tài)恢復(fù)投資增加量、生態(tài)恢復(fù)投資減少量、生態(tài)發(fā)展投資增加量和生態(tài)發(fā)展投資減少量為速率變量；生態(tài)修復(fù)因子、生態(tài)因子、地區(qū)人口總數(shù)、人均生態(tài)承載能力為表函數(shù)；其余為輔助變量。本文設(shè)置的移民投資主要是針對(duì)生態(tài)治理的投資，包括兩個(gè)方面：遷出地生態(tài)恢復(fù)投資和遷入地生態(tài)發(fā)展投資。移民選址和移民安置方式，主要取決于遷出地距離擬遷入地（生態(tài)盈余地區(qū)）的鄰近指數(shù)和遷入地的資源稟賦條件，鑒于遷入地與遷出地之間距離的不易測(cè)量性，本文不對(duì)此展開討論。

（3）生態(tài)福利績(jī)效子系統(tǒng)中，GDP、預(yù)期壽命、教育和生態(tài)足跡為狀態(tài)變量；GDP增加量、GDP減少量、預(yù)期壽命減少量、教育增加量、生態(tài)足跡減少量為速率變量；GDP增長(zhǎng)速度、生態(tài)移民數(shù)量、教育因子、產(chǎn)業(yè)生態(tài)化因子為表函數(shù)；其余為輔助變量。

3.2 系統(tǒng)變量間參數(shù)確定

三江源地區(qū)，是長(zhǎng)江、黃河和瀾滄江三大河流的發(fā)源地，屬于青海省南部范圍，行政區(qū)域上包括四個(gè)藏族自治州的16個(gè)縣和格爾木市的唐古拉山鎮(zhèn)[31-32]。選取三江

源生態(tài)移民為研究對(duì)象，以對(duì)生態(tài)移民生態(tài)可持續(xù)進(jìn)行模擬仿真。2005年《青海三江源自然保護(hù)區(qū)生態(tài)保護(hù)和建設(shè)總體規(guī)劃》的實(shí)行，標(biāo)志著三江源地區(qū)正式開始移民安置遷移工作，考慮到生態(tài)移民前后生態(tài)狀況變化的時(shí)間跨度性和相關(guān)數(shù)據(jù)獲取的難度，選取2005—2013年作為仿真模擬的歷史期間，2014—2020年為生態(tài)移民政策預(yù)測(cè)期間，仿真步長(zhǎng)為1年。

生態(tài)移民生態(tài)可持續(xù)系統(tǒng)變量間的參數(shù)估計(jì)方法及變量間方程如下：

（1）資料查閱法。直接通過(guò)查閱既有文獻(xiàn)、青海省統(tǒng)計(jì)年鑒、中國(guó)區(qū)域經(jīng)濟(jì)統(tǒng)計(jì)年鑒及相關(guān)的統(tǒng)計(jì)公報(bào)獲取系統(tǒng)中一些狀態(tài)變量數(shù)據(jù)，如森林面積、草地面積等歷年數(shù)據(jù)。通過(guò)尋找替代變量，考核系統(tǒng)中的特殊變量，如生產(chǎn)方式主要考察經(jīng)濟(jì)發(fā)展與資源消耗的關(guān)系，用萬(wàn)元生產(chǎn)總能消耗（t標(biāo)煤）、單位GDP二氧化硫排污強(qiáng)度（kg/萬(wàn)元）和工業(yè)固體廢物綜合利用率（%）來(lái)表示，數(shù)據(jù)來(lái)源為青海省年度環(huán)境統(tǒng)計(jì)公報(bào)。

（2）算術(shù)平均法。通過(guò)對(duì)既有數(shù)據(jù)進(jìn)行初步的分析與處理，確定系統(tǒng)中一些變量，如老年人口比例，移民社區(qū)疾病發(fā)生率，森林覆蓋率，砍伐成耕地比例等變量；通過(guò)考核三江源生態(tài)移民地區(qū)的環(huán)境相關(guān)數(shù)據(jù)，依據(jù)《生態(tài)環(huán)境狀況評(píng)價(jià)技術(shù)規(guī)范（試行）》評(píng)價(jià)研究區(qū)域的環(huán)境質(zhì)量指標(biāo)[33]；參照全國(guó)農(nóng)牧民人均純收入，定義三江源生態(tài)移民的移民貧困程度。在生境破碎指數(shù)的計(jì)算中，

結(jié)合人類發(fā)展指數(shù)（HDI）和生態(tài)足跡（EF），將生態(tài)福利績(jī)效（EP）定義為：EP=HDI/EF[30]。以PPSD表示森林面積、草地面積、耕地面積和水土流失面積的標(biāo)準(zhǔn)差，αi為四種類型的地塊面積，MPS為平均地塊面積，則生境破碎指數(shù)（HFI）的計(jì)算公式為：

HFI=PPSDMPS，PPSD=∑4i=1（αi-MPS）23

（3）回歸分析法。如系統(tǒng)中的砍伐因子，人為砍伐面積常量，二者與人為砍伐面積存在如下關(guān)系：移民醫(yī)療投入=財(cái)政支出×移民醫(yī)療投入因子+移民醫(yī)療投入常量，通過(guò)將移民醫(yī)療投入與財(cái)政支出進(jìn)行回歸，R2為0.861 5，移民醫(yī)療投入因子為0.096 5，移民醫(yī)療投入常量為5 351.32。對(duì)于難于估計(jì)的參數(shù)，采取可能值模擬測(cè)試的方法，直至模型行為無(wú)顯著變化時(shí)，確定相應(yīng)值為參數(shù)值[34]。

（4）灰色預(yù)測(cè)模型。根據(jù)仁善英的《三江源生態(tài)移民生產(chǎn)生活效益研究》和杜發(fā)春的《三江源生態(tài)移民研究》三江源地區(qū)生態(tài)移民研究成果，只能獲取2005—2013年間的三江源地區(qū)生態(tài)移民部分?jǐn)?shù)據(jù)（移民數(shù)量，移民收入等），而本系統(tǒng)的仿真年份是2005—2020年，因此缺失年間數(shù)據(jù)需要通過(guò)預(yù)測(cè)得到。根據(jù)鄧聚龍的灰色系統(tǒng)預(yù)測(cè)理論，即從2005—2013年間的已知的數(shù)據(jù)中尋找其內(nèi)在規(guī)律， 再通過(guò)建立灰色模型以預(yù)測(cè)未來(lái)數(shù)值，并將預(yù)測(cè)值輸入生態(tài)移民生態(tài)可持續(xù)SD[35]。

通過(guò)綜合運(yùn)用以上方法，構(gòu)建出生態(tài)移民生態(tài)可持續(xù)SD模型變量間的方程，主要方程如下：

（1）生態(tài)恢復(fù)投資增加量=生態(tài)移民環(huán)境保護(hù)投資×生態(tài)修復(fù)因子；

（2）遷出地生態(tài)恢復(fù)投資=INTEG（生態(tài)恢復(fù)投資增加量×0.322-遷出地的生態(tài)赤字×73.61+4 489.36，5 500）；

（3）生態(tài)發(fā)展投資增加量=生態(tài)移民環(huán)境保護(hù)投資×（生態(tài)因子-生態(tài)修復(fù)因子）；

遷入地生態(tài)發(fā)展投資=INTEG（生態(tài)發(fā)展投資增加量×0.777+產(chǎn)業(yè)結(jié)構(gòu)能耗系數(shù)0.093+2.206， 5 300）；

（4）生態(tài)移民環(huán)境保護(hù)投資=財(cái)政支出×0.315+36.89×DELAY3（生態(tài)盈余，1）-1 092.44；

（5）植樹造林面積=遷出地生態(tài)恢復(fù)投資×植樹造林投資因子；

（6）退耕還林面積=遷出地生態(tài)恢復(fù)投資×退耕還林投資因子；

（7）水土治理面積=遷出地生態(tài)恢復(fù)投資×水土治理投資因子；

（8）退牧還草面積=遷出地生態(tài)恢復(fù)投資×退牧還草投資因子；

（9）森林面積=INTEG（森林面積增加量-森林面積減少量，138.53）；

（10）草地面積=INTEG（草地恢復(fù)面積-草地退化面積，399.989）；

（11）水土流失面積=INTEG（水土流失面積增加量-水土流失治理面積，74.43）；

（12）耕地面積=INTEG（耕地面積增加量-耕地面積減少量，11.16）；

（13）遷出地生態(tài)赤字=生態(tài)移民數(shù)量×（人均生態(tài)承載能力-33.95+0.052×水土流失面積+0.058×草地面積+0.008 96×森林面積+0.260 3×耕地面積）；

（14）GDP=INTEG（GDP增長(zhǎng)，329 673）；

（15）移民教育投入=INTEG（財(cái)政支出×教育因子，24 588）；

（16）移民醫(yī)療投入=財(cái)政支出×0.096 5+5 351.32；

（17）預(yù)期壽命=INTEG（1 069.92×人口死亡率-65.613，66.03）；

（18）移民平均受教育年限=移民教育投入×0.000 024 4+6.944；

（19）人口死亡率=社區(qū)醫(yī)療投入×（-8.51×10-9）+0.062 3；

（20）移民后人均收入=預(yù)期壽命×295.346+12 808.1×移民可就業(yè)比例-21 500.1；

（21）產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)保的支持度=1-產(chǎn)業(yè)結(jié)構(gòu)能耗系數(shù)；

（22）移民產(chǎn)業(yè)生態(tài)化效率=遷入地生態(tài)發(fā)展投資×產(chǎn)業(yè)生態(tài)化因子；

（23）移民產(chǎn)業(yè)蘊(yùn)含工作崗位程度=移民產(chǎn)業(yè)生態(tài)化效率×0.208 9+0.581 7；

（24）平均受教育年限對(duì)環(huán)保的支持度=移民平均受教育年限/22；

（25）移民收入對(duì)環(huán)保的支持度=移民收入×0.000 117 8-0.353；

（26）移民環(huán)保態(tài)度=MIN（MIN（產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)保的支持度，平均受教育年限對(duì)環(huán)保的支持度），移民收入對(duì)環(huán)保的支持度）；

（27）產(chǎn)業(yè)結(jié)構(gòu)能耗系數(shù)=2.495 5-第三產(chǎn)業(yè)所占比重×3.683-第一產(chǎn)業(yè)所占比重×5.91；

（28）生態(tài)足跡減少量=環(huán)保態(tài)度×生態(tài)足跡因子產(chǎn)業(yè)結(jié)構(gòu)能耗系數(shù)；

（29）生態(tài)盈余=（人均生態(tài)承載能力-人均生態(tài)足跡）×地區(qū)人口總數(shù)。

3.3 模型有效性檢驗(yàn)

通過(guò)建立生態(tài)可持續(xù)系統(tǒng)動(dòng)力學(xué)總流圖和確定系統(tǒng)變量間參數(shù)，生態(tài)移民工程中生態(tài)系統(tǒng)重建的SD模型已經(jīng)基本構(gòu)建起來(lái)。為了驗(yàn)證模型的有效性，對(duì)所構(gòu)建的模型進(jìn)行結(jié)構(gòu)適合性檢驗(yàn)和行為一致性檢驗(yàn)，結(jié)合建模目的及三江源生態(tài)移民研究現(xiàn)有文獻(xiàn)，認(rèn)為生態(tài)移民生態(tài)可持續(xù)動(dòng)力學(xué)系統(tǒng)指標(biāo)設(shè)置和方程量綱反映了實(shí)際系統(tǒng)的結(jié)構(gòu)和功能，可以反映三江源地區(qū)生態(tài)移民過(guò)程中生態(tài)系統(tǒng)現(xiàn)狀，也可以用來(lái)預(yù)測(cè)生態(tài)系統(tǒng)發(fā)展的可持續(xù)性趨勢(shì)。

4 系統(tǒng)仿真與政策模擬

本節(jié)將運(yùn)用所構(gòu)建的生態(tài)移民生態(tài)可持續(xù)SD模型進(jìn)行系統(tǒng)仿真與政策模擬，以為生態(tài)移民工程實(shí)施中實(shí)現(xiàn)生態(tài)可持續(xù)提供決策依據(jù)。選取生境破碎度指數(shù)、生態(tài)移民環(huán)境保護(hù)投資、生態(tài)恢復(fù)因子和生態(tài)發(fā)展因子四個(gè)變量為政策變量，自然發(fā)展?fàn)顩r下，三江源地區(qū)生態(tài)移民生態(tài)系統(tǒng)可持續(xù)的仿真結(jié)果如表1所示。

由表1可以看出，生態(tài)福利績(jī)效呈現(xiàn)上升趨勢(shì)，且隨著生態(tài)因子增長(zhǎng)速度加快，生境破碎度明顯下降。主要得益于仿真期間政府對(duì)三江源地區(qū)生態(tài)移民遷出地生態(tài)恢復(fù)和遷入地生態(tài)發(fā)展投資力度的加大。2005年，國(guó)務(wù)院批準(zhǔn)了《青海三江源自然保護(hù)區(qū)生態(tài)保護(hù)和建設(shè)總體規(guī)劃》，投資75億元以開展三江源地區(qū)生態(tài)保護(hù)與建設(shè)一期工程，該工程于2005—2013年實(shí)施，主要包括生態(tài)保護(hù)與建設(shè)項(xiàng)目、農(nóng)牧民生產(chǎn)生活基礎(chǔ)設(shè)施建設(shè)項(xiàng)目和支撐項(xiàng)目三大類[36]。

根據(jù)中國(guó)科學(xué)院地理科學(xué)與資源研究所對(duì)三江源地區(qū)一期工程生態(tài)成效開展了連續(xù)9年的科學(xué)監(jiān)測(cè)與評(píng)估資料，結(jié)果顯示三江源地區(qū)的生態(tài)保護(hù)和改善效果明顯[37]。三江源生態(tài)保護(hù)和建設(shè)二期工程規(guī)劃于2013年12月開始實(shí)施，治理范圍從15.23萬(wàn)km2擴(kuò)大至39.5萬(wàn)km2，旨在通過(guò)生態(tài)系統(tǒng)修復(fù)與生態(tài)工程建設(shè)的有機(jī)結(jié)合，以實(shí)現(xiàn)生態(tài)系統(tǒng)保護(hù)和建設(shè)。

隨著對(duì)生態(tài)系統(tǒng)恢復(fù)和發(fā)展的治理，生境破碎度對(duì)生態(tài)系統(tǒng)可持續(xù)的制約作用減小，三江源地區(qū)居民人類發(fā)展指數(shù)逐步提升，從而提高了生態(tài)福利績(jī)效。

4.1 不同政策方案設(shè)置

在對(duì)自然發(fā)展?fàn)顩r下生態(tài)移民生態(tài)可持續(xù)SD模型進(jìn)行系統(tǒng)仿真的基礎(chǔ)上，發(fā)現(xiàn)生態(tài)移民生態(tài)系統(tǒng)的可持續(xù)發(fā)展存在提升的潛力。生境破碎度指數(shù)直接影響生態(tài)移民生態(tài)系統(tǒng)的可持續(xù)性，而生態(tài)因子（即生態(tài)移民環(huán)境保護(hù)投資）通過(guò)生態(tài)恢復(fù)因子與生態(tài)發(fā)展因子（生態(tài)發(fā)展因子=生態(tài)因子-生態(tài)恢復(fù)因子）二者的分配比例，間接影響生態(tài)移民生態(tài)系統(tǒng)的可持續(xù)性，這些是實(shí)現(xiàn)生態(tài)系統(tǒng)重建與提高生態(tài)系統(tǒng)福利績(jī)效的重要問(wèn)題。本文通過(guò)調(diào)整生境破碎度指數(shù)、生態(tài)移民環(huán)境保護(hù)投資、生態(tài)恢復(fù)因子和生態(tài)發(fā)展因子四個(gè)變量的不同取值，以形成不同的政策模式，從而通過(guò)觀察不同模式下生態(tài)移民生態(tài)可持續(xù)發(fā)展的運(yùn)行情況，為生態(tài)移民生態(tài)系統(tǒng)重建提供決策依據(jù)。

生境破碎度指數(shù)與森林面積、草地面積、耕地面積和水土流失面積直接相關(guān)，受生態(tài)移民環(huán)境保護(hù)投資和生態(tài)恢復(fù)因子的影響，這些變量與生態(tài)發(fā)展因子共同影響生態(tài)福利績(jī)效。通過(guò)改變這些政策變量參數(shù)的不同組合，可以進(jìn)行生態(tài)系統(tǒng)可持續(xù)的政策方針。在這里，主要討論以下六種政策方案，如表2所示。表中第一行為三江源地區(qū)生態(tài)移民生態(tài)系統(tǒng)的自然狀態(tài)，政策模式1～4的設(shè)置假定生態(tài)移民環(huán)境保護(hù)投資在原有的基礎(chǔ)上均增加50%，以此作為其他政策參數(shù)的調(diào)控基礎(chǔ)，每個(gè)方案中針對(duì)遷出地的生態(tài)修復(fù)投資和遷入地的生態(tài)發(fā)展投資的側(cè)重點(diǎn)不一樣。在政策模式1中，假定生態(tài)修復(fù)因子和生態(tài)發(fā)展因子協(xié)調(diào)發(fā)展，增加比重均為25%。為了探討生態(tài)移民生態(tài)福利績(jī)效的提升過(guò)程中，是否應(yīng)該平等對(duì)待遷出地的生態(tài)修復(fù)和遷入地的生態(tài)發(fā)展投資，分別調(diào)整二者的投資比例，以生態(tài)修復(fù)因子和生態(tài)發(fā)展因子增加比重為（30%，20%）、（20%，30%）和（10%和40%）為例，設(shè)置政策模式2、3和4。為了實(shí)現(xiàn)生境破碎指數(shù)最小化和生態(tài)福利績(jī)效最大化，在政策模式5中，設(shè)置生境破碎指數(shù)為0.5的目標(biāo)，遷出地生態(tài)修復(fù)投資和遷入地生態(tài)發(fā)展投資增加比例為（10%和40%）；在政策模式6中，設(shè)置生態(tài)福利績(jī)效為1的目標(biāo)，遷出地生態(tài)修復(fù)投資和遷入地的生態(tài)發(fā)展投資比例為（10%和40%）。

4.2 不同方案的對(duì)比分析

自然狀態(tài)和六種政策模式之下的生態(tài)福利績(jī)效變化對(duì)比情況如圖3所示。從七條仿真結(jié)果折線圖可以看到，政策模式6的生態(tài)福利績(jī)效效果最好，其次為政策模式5。由表2可知，這二種模式的生態(tài)修復(fù)因子占比小于生態(tài)發(fā)展因子，即生態(tài)移民遷出地的生態(tài)修復(fù)投資小于遷入地的生態(tài)發(fā)展投資。這表明通過(guò)改善生態(tài)移民遷入地的生態(tài)發(fā)展性問(wèn)題可以顯著提高生態(tài)移民生態(tài)福利績(jī)效，而單純性地增加生態(tài)移民環(huán)境保護(hù)投資，忽視生態(tài)福利績(jī)效最大化或者生境破碎指數(shù)最小化的政策導(dǎo)向目標(biāo)，會(huì)出現(xiàn)增加生態(tài)移民環(huán)境保護(hù)投資的效果反倒不如自然狀態(tài)下生態(tài)移民效果的情況。

自然狀態(tài)和六種政策模式之下的生境破碎指數(shù)對(duì)比情況如圖4所示。從七條仿真結(jié)果折線圖可以看到，政策模式5的生境破碎指數(shù)運(yùn)行效果最好，在2010年以前政策模式6為次優(yōu)政策選擇，2010年以后政策模式4為次優(yōu)政策選擇，體現(xiàn)了隨著三江源地區(qū)移民規(guī)模的擴(kuò)大，生態(tài)移民的生活安置、就業(yè)和生產(chǎn)發(fā)展等經(jīng)濟(jì)活動(dòng)強(qiáng)度與生境破碎指數(shù)之間的矛盾凸顯。由表2可知，政策模式4、5和6的生態(tài)修復(fù)因子占比小于生態(tài)發(fā)展因子的同時(shí)，分別具備生態(tài)移民環(huán)境保護(hù)投資向遷入地生態(tài)保護(hù)傾斜、以生境破碎指數(shù)最小化和生態(tài)福利績(jī)效最大化為目標(biāo)導(dǎo)向的特點(diǎn)，實(shí)現(xiàn)了生境破碎指數(shù)改善的最佳運(yùn)行效果，即當(dāng)生態(tài)移民環(huán)境保護(hù)投資有限的時(shí)候，只要協(xié)調(diào)遷出地和遷入地的生態(tài)投資，充分利用既有資源和資金條件，是可以實(shí)現(xiàn)生態(tài)移民生態(tài)系統(tǒng)可持續(xù)發(fā)展的。

對(duì)于政策模式1～4，雖然這幾種模式并沒有實(shí)現(xiàn)生態(tài)福利績(jī)效最佳，但是相比于自然狀態(tài)，這幾種模式均對(duì)生境破碎指數(shù)有所改善，主要原因在于生態(tài)移民環(huán)境保護(hù)的穩(wěn)步增加。而在這三種政策模式中，政策模式4的運(yùn)行效果最好，原因在于生態(tài)移民環(huán)境保護(hù)不變的情況下，遷出地的生態(tài)修復(fù)投資占比最小。

綜合對(duì)比自然狀態(tài)和六種政策模式的生態(tài)福利績(jī)效和生境破碎指數(shù)仿真結(jié)果，政策模式5為最優(yōu)選擇。即只注重提高生態(tài)移民環(huán)境保護(hù)投資，而忽視遷出地的生態(tài)修復(fù)投資和遷入地的生態(tài)發(fā)展投資比重分配，并不能有效提高生態(tài)移民的生態(tài)福利績(jī)效；同樣，只調(diào)整遷出地的生態(tài)修復(fù)投資和遷入地的生態(tài)發(fā)展投資比重，而不注重生態(tài)投資的運(yùn)作效率，也會(huì)影響生境破碎指數(shù)的改善。因此，為了實(shí)現(xiàn)生態(tài)移民生態(tài)系統(tǒng)的可持續(xù)發(fā)展，不僅要提高生態(tài)移民環(huán)境保護(hù)投資，協(xié)調(diào)生態(tài)移民遷出地和遷入地的投資分配比重，而且還要提高生態(tài)投資的運(yùn)作效率，即實(shí)現(xiàn)生境破碎指數(shù)最小化的運(yùn)作模式是促使生態(tài)移民生態(tài)可持續(xù)發(fā)展的最佳政策實(shí)施模式。

5 結(jié)論和建議

生態(tài)移民的過(guò)程，是遷出地的自然生態(tài)環(huán)境問(wèn)題和遷入地的發(fā)展性生態(tài)環(huán)境問(wèn)題與“人”的因素交織在一起的過(guò)程，為了實(shí)現(xiàn)生態(tài)移民涉及區(qū)域生態(tài)系統(tǒng)的重建與可持續(xù)發(fā)展，本文構(gòu)建了生態(tài)移民生態(tài)可持續(xù)發(fā)展系統(tǒng)動(dòng)力學(xué)模型，并輔以三江源地區(qū)的自然資源和生態(tài)投資等數(shù)據(jù)進(jìn)行了系統(tǒng)仿真和相關(guān)政策模擬，得出如下結(jié)論：

第一，生態(tài)移民生態(tài)可持續(xù)發(fā)展與生境破碎指數(shù)和生態(tài)福利績(jī)效密切相關(guān)，生境破碎指數(shù)與生態(tài)移民環(huán)境投資密切相關(guān)。三江源生態(tài)移民的實(shí)現(xiàn)過(guò)程，伴隨著生態(tài)移民生產(chǎn)生活條件的改善，實(shí)現(xiàn)了生態(tài)環(huán)境保護(hù)的斐然成績(jī)，截止2015年實(shí)現(xiàn)生態(tài)移民受益人口約達(dá)234 771人，農(nóng)牧民人均純收入年均增長(zhǎng)10%以上，同時(shí)實(shí)現(xiàn)了草地均產(chǎn)量、森林新增覆蓋面積等五個(gè)增加。在建立生態(tài)移民可持續(xù)發(fā)展系統(tǒng)仿真時(shí)，發(fā)現(xiàn)生態(tài)移民環(huán)境保護(hù)投資影響遷出地的生態(tài)修復(fù)投資，并通過(guò)構(gòu)造植樹造林投資因子、退耕還林投資因子、水土治理投資因子和退牧還草投資因子描述了生態(tài)脆弱與貧困的惡性循環(huán)關(guān)系影響的同時(shí)，也在影響著生境破碎指數(shù)。

第二，生態(tài)移民環(huán)境保護(hù)投資的增加有助于降低生境破碎指數(shù)和提高生態(tài)福利績(jī)效，尤其是對(duì)遷入地的生態(tài)發(fā)展投資比重大于對(duì)遷出地的生態(tài)修復(fù)投資時(shí)。根據(jù)生態(tài)移民生態(tài)可持續(xù)系統(tǒng)仿真結(jié)果可知，在六種政策模式中，生態(tài)移民環(huán)境保護(hù)投資均有所增加，政策模式5和6的生態(tài)福利績(jī)效相比三江源自然狀態(tài)的生態(tài)移民效果均有所改善。政策模式1～5的生境破碎指數(shù)仿真模擬結(jié)果優(yōu)于自然狀態(tài)，原因在于生態(tài)移民環(huán)境保護(hù)投資一定時(shí)，遷入地的生態(tài)發(fā)展投資比重高于對(duì)遷出地的生態(tài)修復(fù)投資。

第三，生態(tài)移民生態(tài)可持續(xù)發(fā)展應(yīng)以生境破碎指數(shù)最小化為目標(biāo)，并注重生態(tài)移民環(huán)境保護(hù)投資資金的分配比重和運(yùn)作效率。在六種政策模式中，政策模式5運(yùn)行實(shí)現(xiàn)的生境破碎指數(shù)最低和生態(tài)福利績(jī)效最高，即運(yùn)行效果最佳。自然狀態(tài)下，三江源地區(qū)生態(tài)移民過(guò)程中的生境破碎指數(shù)為0.927 7、生態(tài)福利績(jī)效為0.592 6，在政策模式5之下，仿真期間所實(shí)現(xiàn)的生境破碎指數(shù)平均值為0.5，生態(tài)福利績(jī)效的平均值為1.187 5。雖然政策模式2～6的生境破碎指數(shù)和生態(tài)福利績(jī)效仿真模擬結(jié)果有所改善，但是政策模式5和6以生境破碎指數(shù)最小化和生態(tài)福利績(jī)效最大化為導(dǎo)向的運(yùn)行結(jié)果遠(yuǎn)優(yōu)于政策模式3和4，并且政策模式5的運(yùn)行效果最佳，說(shuō)明了生態(tài)移民生態(tài)可持續(xù)發(fā)展不但依賴于對(duì)遷出地和遷入地的生態(tài)環(huán)境保護(hù)投資的增加，更依賴于遷出地和遷入地的生態(tài)投資資金分配比重以及投資資金運(yùn)作效率。

基于上述結(jié)論，為實(shí)現(xiàn)三江源地區(qū)生態(tài)移民生態(tài)系統(tǒng)的可持續(xù)發(fā)展與重建，本文提出以下政策建議：一是增加三江源地區(qū)生態(tài)移民環(huán)境保護(hù)投資力度，保證遷出地生態(tài)修復(fù)與遷入地生態(tài)發(fā)展的資金供給，實(shí)現(xiàn)移民遷出地的生態(tài)修復(fù)問(wèn)題，并解決移民遷入地的生態(tài)發(fā)展問(wèn)題。二是嚴(yán)格監(jiān)測(cè)評(píng)估生態(tài)移民遷出地和遷入地的生態(tài)投資比例，適當(dāng)?shù)亟档鸵泼襁w出地區(qū)的生態(tài)修復(fù)投資，提高移民遷入地的生態(tài)發(fā)展投資。由于生態(tài)移民改變了人們利用生態(tài)系統(tǒng)服務(wù)的位置或變換了生態(tài)系統(tǒng)服務(wù)供給的地點(diǎn)，因此，應(yīng)協(xié)調(diào)好移民遷出地和遷入地的生態(tài)投資比重，并通過(guò)提高遷出地的森林面積、草地面積等來(lái)降低生境破碎指數(shù)，通過(guò)提高遷入地的GDP、發(fā)展教育和醫(yī)療、降低生態(tài)足跡等來(lái)提高生態(tài)福利績(jī)效。三是提高生態(tài)移民環(huán)境保護(hù)投資運(yùn)作效率，促進(jìn)生態(tài)移民生態(tài)系統(tǒng)可持續(xù)發(fā)展。研究發(fā)現(xiàn)，在含自然狀態(tài)的7中政策模式中，政策模式5的運(yùn)作效果最好。以生境破碎指數(shù)最小化為目標(biāo)，政策5在實(shí)現(xiàn)生境破碎指數(shù)最小化的同時(shí)，也實(shí)現(xiàn)了生態(tài)福利績(jī)效的最大化。然而，處于自然狀態(tài)的三江源地區(qū)的生態(tài)移民工作，還未實(shí)現(xiàn)生境破碎指數(shù)的最小化和生態(tài)福利績(jī)效的最大化。因此，本文建議未來(lái)三江源地區(qū)生態(tài)移民工作應(yīng)該審視生態(tài)移民遷出地生態(tài)修復(fù)投資與遷入地生態(tài)發(fā)展投資的辯證關(guān)系，以及生態(tài)移民環(huán)境保護(hù)投資的運(yùn)作效率問(wèn)題，將生境破碎指數(shù)最小化和生態(tài)福利績(jī)效最大化目標(biāo)導(dǎo)向納入生態(tài)移民生態(tài)系統(tǒng)重建與可持續(xù)發(fā)展戰(zhàn)略中。

（編輯：王愛萍）

參考文獻(xiàn)

[1]FAHRIG L. Effects of habitat fragmentation on biodiversity[J].Annual review of ecology， evolution， and systematics， 2003，34：487-515.

[2]DEBINSKI D M， HOLT R D. A survey and overview of habitat fragmentation experiments[J].Conservation biology， 2000，14（2）：342-355.

[3]ROBINSON G R，HOLT R D， GAINES M S， et al. Diverse and contrasting effects of habitat fragmentation[J].Science， 1992，257（5069）：524-526.

[4]KRUESS A. Effects of landscape structure and habitat type on a plantherbivoreparasitoid community[J].Ecography，2003， 26（3）：283-290.

[5]戚仁海.生境破碎化對(duì)城市化地區(qū)生物多樣性影響的研究——以蘇州為例[D].上海：華東師范大學(xué)，2008：6-7.

[6]武正軍，李義明.生境破碎化對(duì)動(dòng)物種群存活的影響[J].生態(tài)學(xué)報(bào)，2003（11）：2424-2435.

[7]張?zhí)煨?，劉鳳.城市邊緣區(qū)空間破碎化的指標(biāo)體系及其應(yīng)用研究——以北京市海淀區(qū)上地、清河街道為例[J].城市發(fā)展研究，2014（5）：24-29.

[8]楊芳，賀達(dá)漢.生境破碎化對(duì)生物多樣性的影響[J].生態(tài)科學(xué)，2006（12）：564-567.

[9]HEILMAN G E， STRITTHOLT J R， SLOSSER N C， et al. Forest fragmentation of the conterminous United States：assessing forestintactness through road density and spatial characteristics[J].Bioscience，2002，52： 411-422.

[10]CHEN L D，LIU X H，F(xiàn)U B J.Evaluation on giant panda habitat fragmentation in Wolong Nature Reserve[J].Acta ecologica sinica，1999，19（ 3）： 291-297.

[11]BEISSINGER S R，WESTPHAL M I. On the use of demographic models of population viability in endang ered species management[J].Journal of wildlife management， 1998， 62： 821-841.

[12]姜廣順，張明海，馬建章.黑龍江省完達(dá)山地區(qū)馬鹿生境破碎化及其影響因子[J].生態(tài)學(xué)報(bào)，2005（7）：1691-1698.

[13]劉學(xué)敏.西北地區(qū)生態(tài)移民的效果與問(wèn)題探討[J].中國(guó)農(nóng)村經(jīng)濟(jì)，2002（4）：47-52.

[14]包智明.關(guān)于生態(tài)移民的定義、分類及若干問(wèn)題[J].中央民族大學(xué)學(xué)報(bào)（哲學(xué)社會(huì)科學(xué)版），2006（1）：27-31.

[15]張麗君，王菲.中國(guó)西部牧區(qū)生態(tài)移民后續(xù)發(fā)展對(duì)策探析[J].中央民族大學(xué)學(xué)報(bào)（哲學(xué)社會(huì)科學(xué)版），2011（4）：31-36.

[16]色音，張繼焦.生態(tài)移民的環(huán)境社會(huì)學(xué)研究[M].北京：民族出版社，2009：236-239.

[17]李笑春，陳志.對(duì)生態(tài)移民的理性思考——以渾善達(dá)克沙地為例[J].內(nèi)蒙古大學(xué)學(xué)報(bào)（人文社會(huì)科學(xué)版），2004（5）：34-38.

[18]侯東民.草原人口生態(tài)壓力持續(xù)增長(zhǎng)態(tài)勢(shì)與解決方法——經(jīng)濟(jì)誘導(dǎo)式生態(tài)移民工程的可行性分析[J].中國(guó)人口科學(xué)，2002（5）：32-37.

[19]侯東民.移民經(jīng)濟(jì)政策市場(chǎng)化：內(nèi)地“買地移民”的調(diào)查[J].人口與發(fā)展，2008（4）：57-62.

[20]孟琳琳，包智明.生態(tài)移民研究綜述[J].中央民族大學(xué)學(xué)報(bào)（哲學(xué)社會(huì)科學(xué)版），2004（6）：48-52.

[21]劉學(xué)敏，陳靜.生態(tài)移民、城鎮(zhèn)化與產(chǎn)業(yè)發(fā)展——對(duì)西北地區(qū)城鎮(zhèn)化的調(diào)查與思考[J].中國(guó)特色社會(huì)主義研究，2002（2）：61-63.

[22]葛根高娃，烏云巴圖.內(nèi)蒙古牧區(qū)生態(tài)移民的概念、問(wèn)題與對(duì)策[J].內(nèi)蒙古社會(huì)科學(xué)，2003（2）：118-122.

[23]劉慧，葉兒肯·吾扎提.中國(guó)西部地區(qū)生態(tài)扶貧策略研究[J].中國(guó)人口·資源與環(huán)境，2013（10）：52-58.

[24]張建軍，生態(tài)移民與環(huán)境治理——基于塔里木河流域的經(jīng)驗(yàn)考察[J].西南民族大學(xué)學(xué)報(bào)（人文社會(huì)科學(xué)版），2015（9）：56-62.

[25]賈耀鋒，中國(guó)生態(tài)移民效益評(píng)估研究綜述[J].資源科學(xué)，2016（8）：1550-1560.

[26]李柏山，水資源開發(fā)利用對(duì)漢江流域水生態(tài)環(huán)境影響及生態(tài)系統(tǒng)健康評(píng)價(jià)研究[D].武漢：武漢大學(xué)，2013：48-51.

[27]薛建春，基于生態(tài)足跡模型的礦區(qū)復(fù)合生態(tài)系統(tǒng)分析及動(dòng)態(tài)預(yù)測(cè)[D].北京：中國(guó)地質(zhì)大學(xué)，2010：53-55.

[28]李華，基于系統(tǒng)仿真和情景模擬的崇明生態(tài)安全評(píng)估[D].上海：華東師范大學(xué)，2010：40-43.

[29]張晶，產(chǎn)業(yè)生態(tài)系統(tǒng)的定量解析與評(píng)價(jià)及仿真[D].徐州：中國(guó)礦業(yè)大學(xué)，2012：20-23.

[30]諸大建，張帥.生態(tài)福利績(jī)效與深化可持續(xù)發(fā)展的研究[J].同濟(jì)大學(xué)大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版），2014（5）：106-115.

[31]仁善英.三江源生態(tài)移民生產(chǎn)生活效益研究[M].天津：天津大學(xué)出版社，2016：19-20.

[32]杜發(fā)春.三江源生態(tài)移民研究[M].北京：中國(guó)社會(huì)科學(xué)出版社，2016：19-20.

[33]楊皓然.三江源生態(tài)經(jīng)濟(jì)系統(tǒng)耦合機(jī)制研究[M].北京：中國(guó)社會(huì)科學(xué)出版社，2016：46-47.

[34]曾麗君，隋映輝，申玉三.科技產(chǎn)業(yè)與資源型城市可持續(xù)協(xié)同發(fā)展的系統(tǒng)動(dòng)力學(xué)研究[J].中國(guó)人口·資源與環(huán)境，2014（10）：85-93.

[35]鄧聚龍.灰預(yù)測(cè)與灰決策[M].武漢：華中科技大學(xué)出版社，2002：45-130.

[36]芳旭.三江源生態(tài)保護(hù)和建設(shè)十年成績(jī)斐然[N/OL].青海：西寧晚報(bào)報(bào)，2016-09-16[2018-01-13]. http：//www.forestry.gov.cn/main/441/content-905683.html.

[37]中國(guó)科學(xué)院地理科學(xué)與資源研究所.三江源生態(tài)保護(hù)和建設(shè)一期工程生態(tài)成效評(píng)估報(bào)告[R].2015.

Abstract Ecological migration project is one of the hot issues in the study of migration and ecological environment in recent years， and it is a governmentled systematic project for ecological protection and poverty alleviation will focus on poverty alleviation and ecological protection in ecologically vulnerable areas， based on the principle of coordinated and sustainable development of population， resources， environment， economy and society.The basic premise of the sustainable resettlement and development of ecological immigrants is the sustainable development of ecological migration projects and ecosystems.Through constructing the ecological sustainable dynamic system model of ecological migration project， reveals the dynamic feedback mechanism among habitat fragmentation， ecological migration engineering and Ecological Welfare Performance，makes the case of ecological immigrants in Threeriversource region to simulate the immigration sustainable ecological system before and after the immigration， to provide the basis for decisionmaking in the implementation of ecological migration project ecological sustainable development. The structure suitability and behavior consistency test of the model show that the model is true and effective in reflecting the structure and function of the actual ecological immigration system. By adjusting the four policy variables of habitat fragmentation index， ecological migration environmental protection investment， ecological restoration factor and ecological development factor， set up six models to simulate the ecological immigrants sustainable development policy.Policy simulation shows：policy model 1-6 reflect the increasing of environmental protection investment helps to reduce the habitat fragmentation index and improve the ecological welfare performance in ecological migration areas； policy model 3-6 reflect environmental protection investment should focus on the development of ecological immigration immigration，appropriate to reduce the proportion of ecological emigration investment；policy model 5-6 illustrates the operation efficiency of ecological environmental protection investment affects the ecological migration ecosystem sustainable development. Therefore， it is important to increase environmental protection investment and adjust the proportion of ecological investment in the areas of relocation for achieving the reconstruction and sustainable development of ecological migration ecosystems. On this basis， improving the efficiency of ecological investment is the best policy model to promote the sustainable development of ecological migration projects and ecosystems.

Key words system dynamics； habitat fragmentation； ecological migration project； ecological welfare performance