鐘水映　馮英杰

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

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

中圖分類號 F062.2 文獻標識碼 A 文章編號 1002-2104（2018）11-0010-10 DOI：10.12062/cpre.20180516

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

1 文獻綜述

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

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

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

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

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

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

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

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

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

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

綜上所述，當前學者對于生態(tài)可持續(xù)發(fā)展的研究，涉及生境破碎化、生態(tài)移民中的生態(tài)修復與生態(tài)—經(jīng)濟—社會系統(tǒng)的協(xié)調(diào)發(fā)展、基于系統(tǒng)動力學的對于生態(tài)系統(tǒng)健康安全評估和生態(tài)系統(tǒng)服務(wù)價值估算等多個方面，而缺少對生境破碎化、生態(tài)移民及生態(tài)福利績效三者相互影響機制的研究。本文試圖構(gòu)建生境破碎化、生態(tài)移民、生態(tài)福利績效的系統(tǒng)動力學（SD）模型，并借助系統(tǒng)動力學可視化建模工具 Vensim PLE 6.2版本軟件，對生態(tài)移民遷入地和遷出地的生態(tài)系統(tǒng)的可持續(xù)性進行仿真模擬，深入系統(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)動力學模型構(gòu)建

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

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

生境破碎化子系統(tǒng)主要研究移民遷出地區(qū)人們以土地過度開墾、砍伐森林、污水排放等形式對生態(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)移民工程如何通過人口密度重組、資源重配和結(jié)構(gòu)重調(diào)，實現(xiàn)移民與遷入地區(qū)的人口、資源和環(huán)境的可持續(xù)發(fā)展，表現(xiàn)在移民選址與安置、移民就業(yè)、移民社區(qū)各個方面，設(shè)定移民投資、移民數(shù)量、移民就業(yè)等變量；生態(tài)福利績效子系統(tǒng)主要研究遷出地區(qū)的生態(tài)修復與遷入地區(qū)發(fā)展生態(tài)產(chǎn)業(yè)如何通過調(diào)整人類發(fā)展指數(shù)和生態(tài)足跡以影響生態(tài)可持續(xù)發(fā)展，設(shè)定人類發(fā)展指數(shù)、生態(tài)足跡等變量。這三個子系統(tǒng)是相互影響與相互制約的，每個子系統(tǒng)的運行既取決于子系統(tǒng)內(nèi)部的結(jié)構(gòu)，也取決于與其他子系統(tǒng)的關(guān)系，例如生境破碎化子系統(tǒng)直接影響生態(tài)能否可持續(xù)發(fā)展，但是生態(tài)移民工程的實施卻有助于改善生境破碎化，通過生態(tài)效果的逐步改善實現(xiàn)生態(tài)可持續(xù)發(fā)展。在相關(guān)文獻、調(dià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)移民遷出地與遷入地的實際距離和預期距離等不確定性因素，為了表征遷出地的脆弱性生態(tài)系統(tǒng)，僅考慮了生境破碎度量化中的關(guān)鍵因素，即以森林面積、草地面積、水土流失面積和耕地面積作為密度和差異指數(shù)，以人為砍伐面積和過度放牧（放牧牲畜數(shù)量）作為干擾度指數(shù)來表征生境破碎度。

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

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

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

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

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

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

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

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

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

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

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

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

遷入地生態(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)福利績效子系統(tǒng)因果關(guān)系分析

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

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

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

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

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

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

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

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

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

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

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

通過對生態(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個，速率變量包括森林面積增加量、森林面積減少量等17個，表函數(shù)為水土治理投資因子等17個，其他變量均為輔助變量，共74個，具體設(shè)置如下：

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

遷入地生態(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)境保護投資=財政支出×0.315+36.89×DELAY3（生態(tài)盈余，1）-1 092.44；

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

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

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

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

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

（10）草地面積=INTEG（草地恢復面積-草地退化面積，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增長，329 673）；

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

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

（17）預期壽命=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）移民后人均收入=預期壽命×295.346+12 808.1×移民可就業(yè)比例-21 500.1；

（21）產(chǎn)業(yè)結(jié)構(gòu)對環(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è)蘊含工作崗位程度=移民產(chǎn)業(yè)生態(tài)化效率×0.208 9+0.581 7；

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

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

（26）移民環(huán)保態(tài)度=MIN（MIN（產(chǎn)業(yè)結(jié)構(gòu)對環(huán)保的支持度，平均受教育年限對環(huán)保的支持度），移民收入對環(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 模型有效性檢驗

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

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

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

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

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

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

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

在對自然發(fā)展狀況下生態(tài)移民生態(tài)可持續(xù)SD模型進行系統(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)境保護投資）通過生態(tài)恢復因子與生態(tài)發(fā)展因子（生態(tài)發(fā)展因子=生態(tài)因子-生態(tài)恢復因子）二者的分配比例，間接影響生態(tài)移民生態(tài)系統(tǒng)的可持續(xù)性，這些是實現(xiàn)生態(tài)系統(tǒng)重建與提高生態(tài)系統(tǒng)福利績效的重要問題。本文通過調(diào)整生境破碎度指數(shù)、生態(tài)移民環(huán)境保護投資、生態(tài)恢復因子和生態(tài)發(fā)展因子四個變量的不同取值，以形成不同的政策模式，從而通過觀察不同模式下生態(tài)移民生態(tài)可持續(xù)發(fā)展的運行情況，為生態(tài)移民生態(tài)系統(tǒng)重建提供決策依據(jù)。

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

4.2 不同方案的對比分析

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

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

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

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

5 結(jié)論和建議

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

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

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

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

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

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

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

[8]楊芳，賀達漢.生境破碎化對生物多樣性的影響[J].生態(tài)科學，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]姜廣順，張明海，馬建章.黑龍江省完達山地區(qū)馬鹿生境破碎化及其影響因子[J].生態(tài)學報，2005（7）：1691-1698.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

[35]鄧聚龍.灰預測與灰決策[M].武漢：華中科技大學出版社，2002：45-130.

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

[37]中國科學院地理科學與資源研究所.三江源生態(tài)保護和建設(shè)一期工程生態(tài)成效評估報告[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