劉彥花 葉國華

[摘 要]基于2008-2015年歷史數(shù)據(jù)，分析廣西北部灣經(jīng)濟(jì)區(qū)人口-資源-環(huán)境-發(fā)展（PRED）現(xiàn)狀，利用系統(tǒng)動力學(xué)原理，構(gòu)建廣西北部灣經(jīng)濟(jì)區(qū)PRED結(jié)構(gòu)模型，仿真模擬研究區(qū)2016-2025年原型發(fā)展模式、高速發(fā)展模式、穩(wěn)定發(fā)展模式三種不同情景下的PRED發(fā)展?fàn)顩r。結(jié)果表明：以工業(yè)發(fā)展為主導(dǎo)的原型發(fā)展模式和單方面追求經(jīng)濟(jì)增長的高速發(fā)展模式均不利于廣西北部灣經(jīng)濟(jì)區(qū)的可持續(xù)發(fā)展，新常態(tài)下的經(jīng)濟(jì)穩(wěn)步增長，不斷優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)，提高資源利用率，注重環(huán)境治理與保護(hù)的穩(wěn)定發(fā)展模式是當(dāng)前廣西北部灣經(jīng)濟(jì)區(qū)PRED協(xié)調(diào)發(fā)展的最佳方案。

[關(guān)鍵詞]系統(tǒng)動力學(xué)；PRED；情景模擬；廣西北部灣經(jīng)濟(jì)區(qū)

[中圖分類號]F127.9 [文獻(xiàn)標(biāo)識碼]A

系統(tǒng)動力學(xué)（System Dynamic，簡稱SD）是美國Forrester教授于1956年創(chuàng)建的一種研究系統(tǒng)復(fù)雜行為的方法。該方法基于系統(tǒng)論、信息論、控制論、計(jì)算機(jī)仿真模擬等技術(shù)，通過分析復(fù)雜系統(tǒng)的結(jié)構(gòu)、功能及動態(tài)行為間的相互關(guān)系，模擬復(fù)雜系統(tǒng)處于不同情景下的狀態(tài)和變化趨勢。從系統(tǒng)的角度，PRED（Population、Resources、Environment、Development）是人口、資源、環(huán)境與發(fā)展相互制約且相互促進(jìn)的復(fù)雜巨系統(tǒng)，PRED系統(tǒng)各要素間的協(xié)同狀態(tài)表征了一個(gè)區(qū)域的可持續(xù)發(fā)展水平。文獻(xiàn)分析表明：近年來，系統(tǒng)動力學(xué)作為地球科學(xué)重要的研究方法，被國外學(xué)者應(yīng)用于能源、土地、水資源、慢性疾病等領(lǐng)域的分析與模擬；國內(nèi)學(xué)者則主要應(yīng)用于模擬資源、環(huán)境及區(qū)域經(jīng)濟(jì)的可持續(xù)發(fā)展?fàn)顟B(tài)。

廣西北部灣經(jīng)濟(jì)區(qū)包括南寧、北海、欽州、防城港、玉林、崇左六市所轄行政區(qū)域，總面積達(dá)7.27萬km2，占廣西全區(qū)的30.7%。廣西北部灣經(jīng)濟(jì)區(qū)自2008年國務(wù)院批準(zhǔn)實(shí)施《廣西北部灣經(jīng)濟(jì)區(qū)發(fā)展規(guī)劃》以來，發(fā)展勢頭強(qiáng)勁。2008年末區(qū)域生產(chǎn)總值3031.82億元，占廣西全區(qū)生產(chǎn)總值的43.18%。到2015年末區(qū)域生產(chǎn)總值達(dá)7995.88億元，占廣西全區(qū)生產(chǎn)總值的47.59%。然而，隨著廣西北部灣經(jīng)濟(jì)區(qū)經(jīng)濟(jì)的快速發(fā)展，不可避免地伴生著人口增長、資源消耗和環(huán)境污染，如何保持區(qū)域經(jīng)濟(jì)發(fā)展的可持續(xù)狀態(tài)，實(shí)現(xiàn)區(qū)域人口、資源、環(huán)境與經(jīng)濟(jì)增長的協(xié)調(diào)發(fā)展，是目前廣西北部灣經(jīng)濟(jì)區(qū)亟需解決的問題。因此，科學(xué)地分析并動態(tài)識別廣西北部灣經(jīng)濟(jì)區(qū)PRED的協(xié)調(diào)關(guān)系，對于了解區(qū)域可持續(xù)發(fā)展?fàn)顟B(tài)具有重要的意義。本研究基于系統(tǒng)動力學(xué)方法，構(gòu)建廣西北部灣經(jīng)濟(jì)區(qū)六市人口、資源、環(huán)境與經(jīng)濟(jì)系統(tǒng)的協(xié)調(diào)發(fā)展動態(tài)模型，通過改變各子系統(tǒng)的模型參數(shù)，對研究區(qū)未來發(fā)展?fàn)顟B(tài)進(jìn)行仿真模擬預(yù)測，以期幫助人們更好地了解廣西北部灣經(jīng)濟(jì)區(qū)人口、資源、環(huán)境與經(jīng)濟(jì)的協(xié)調(diào)發(fā)展?fàn)顟B(tài)，為制定區(qū)域可持續(xù)發(fā)展措施提供依據(jù)。

1 廣西北部灣經(jīng)濟(jì)區(qū)PRED演變過程

1.1 人口

2008～2015年間，廣西北部灣經(jīng)濟(jì)區(qū)戶籍總?cè)丝谠鲩L8.76%，常住總?cè)丝谧兓幻黠@（圖1），八年間，戶籍總?cè)丝谶h(yuǎn)大于常住人口，說明研究區(qū)內(nèi)人口的流動表現(xiàn)為流出高于流入人口，凈流出人口數(shù)占戶籍人口比例由2008年的6.28%增加到2015年的13.6%，研究區(qū)的經(jīng)濟(jì)發(fā)展情況還不至于留住更多的本人人口，也不足以吸引更多的外地人口流入。八年間區(qū)域內(nèi)城鎮(zhèn)人口增長19.86%。城鎮(zhèn)人口的增長速度明顯高于總?cè)丝诘脑鲩L速度，說明研究區(qū)內(nèi)快速的城鎮(zhèn)化進(jìn)程。而農(nóng)村常住人口在不斷減少，表明大量的農(nóng)村勞動力遷入城市服務(wù)于城鎮(zhèn)建設(shè)或流出外地工作。通過對研究區(qū)內(nèi)各市轄區(qū)三大產(chǎn)業(yè)從業(yè)人口統(tǒng)計(jì)，八年間，各市轄區(qū)第一產(chǎn)業(yè)從業(yè)人數(shù)減少率為33.17%，第二產(chǎn)業(yè)從業(yè)人數(shù)增長率為107.29%，第三產(chǎn)業(yè)從業(yè)人數(shù)增長率為40.83%。

1.2 資源

2008～2015年間，廣西北部灣經(jīng)濟(jì)區(qū)常用耕地面積有小幅增長，增長率約16.47%，但其增幅低于農(nóng)村人口的增幅18.27%；研究區(qū)內(nèi)森林面積增長率約9.33%，相比于人口的增長，八年間人均森林面積增長率為0.53%；研究區(qū)內(nèi)水資源總量小幅波動且整體呈現(xiàn)減少趨勢，八年間水資源減少率為18.95%，在人口增長的同時(shí)，其人均水資源量減少率達(dá)25.48%。

1.3 環(huán)境

2008～2015年間，廣西北部灣經(jīng)濟(jì)區(qū)億元GDP能源消耗量減少率為32.7%，生活垃圾清運(yùn)量增長率為57.4%，城鎮(zhèn)人口人均生活垃圾排放量增長率為31.4%，農(nóng)田化肥使用量增長率為16.9%，單位面積耕地化肥使用量增長率為4.1%，污水排放量增長率為42.9%，人均污水排放量增長率為19.2%，單位產(chǎn)值污水排放量減少率為48.5%，污水處理量增長率為61.4%，污水處理廠集中處理能力增長率為62.9%，污水處理率由2008年77.9%提高到2015年88%。

1.4 經(jīng)濟(jì)發(fā)展

自2008年以來，廣西北部灣進(jìn)入了新的發(fā)展時(shí)期，GDP呈現(xiàn)逐年增長的趨勢。在產(chǎn)業(yè)結(jié)構(gòu)調(diào)整上，政府著力推動產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化升級，加大對戰(zhàn)略性新興產(chǎn)業(yè)、現(xiàn)代服務(wù)業(yè)和現(xiàn)代農(nóng)業(yè)等領(lǐng)域的支持力度，深入實(shí)施“壯二提三強(qiáng)一”的戰(zhàn)略，構(gòu)建現(xiàn)代產(chǎn)業(yè)體系。通過不斷優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)，完善產(chǎn)業(yè)布局，廣西北部灣第一、二、三產(chǎn)業(yè)比重由2008年的21.40：36.64：41.96調(diào)整為2015年的15.31：43.03：41.66，各產(chǎn)業(yè)的增長軌跡如圖2。由圖可知，2008年以來，廣西北部灣經(jīng)濟(jì)區(qū)各產(chǎn)業(yè)均為增速狀態(tài)，其中第二產(chǎn)業(yè)和第三產(chǎn)業(yè)總體上比第一產(chǎn)業(yè)增長快速，第二產(chǎn)業(yè)的增速主要體現(xiàn)在工業(yè)產(chǎn)值的快速增長，說明廣西北部灣經(jīng)濟(jì)區(qū)工業(yè)和服務(wù)業(yè)發(fā)展勢頭強(qiáng)勁。且在經(jīng)濟(jì)區(qū)開發(fā)的前幾年，經(jīng)濟(jì)增速過快，隨著經(jīng)濟(jì)發(fā)展對資源與環(huán)境的影響，經(jīng)濟(jì)發(fā)展速度逐漸放緩，經(jīng)濟(jì)區(qū)的經(jīng)濟(jì)發(fā)展呈現(xiàn)不斷調(diào)整與平衡的趨勢。

2 區(qū)域PRED情景模擬

2.1 SD模型數(shù)據(jù)及參數(shù)設(shè)定

基于2008-2016年《廣西統(tǒng)計(jì)年鑒》《廣西北部灣統(tǒng)計(jì)年鑒》、《廣西年鑒》等文字報(bào)告及數(shù)據(jù)資料，通過提取廣西北部灣經(jīng)濟(jì)區(qū)人口、資源、環(huán)境及經(jīng)濟(jì)發(fā)展相關(guān)信息，分析廣西北部灣經(jīng)濟(jì)區(qū)PRED演變過程，建立廣西北部灣經(jīng)濟(jì)區(qū)PRED系統(tǒng)SD流圖（圖3）。系統(tǒng)模型的時(shí)間界限為2008-2025年，基底年份為2008年，模擬時(shí)間步長為1年，2008-2015年為模型模擬與檢驗(yàn)階段，即運(yùn)用該階段的歷史數(shù)據(jù)設(shè)定模型參數(shù)、調(diào)整并檢驗(yàn)?zāi)Ｐ瓦\(yùn)行結(jié)果；2016-2025為情景模擬階段，即基于三種不同的設(shè)定方案改變模型參數(shù)，開展不同情景下的PRED系統(tǒng)仿真模擬。

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

選取研究區(qū)域2008-2015年數(shù)據(jù)，將模型中各變量的模擬值與實(shí)際值進(jìn)行對比分析，結(jié)果顯示真實(shí)值與模擬值的相對誤差基本處于-10%～10%之間，說明系統(tǒng)模型對于模擬研究區(qū)PRED狀態(tài)具有較好的可靠性。本論文選取總?cè)丝?、總用水量、有效灌溉面積、能源消耗、地區(qū)生產(chǎn)總值五個(gè)指標(biāo)變量展示模型行為與實(shí)際系統(tǒng)的一致性檢驗(yàn)結(jié)果（表1）。

2.3 情景模擬

根據(jù)模型檢驗(yàn)結(jié)果，說明模型的結(jié)構(gòu)合理，可用于反映廣西北部灣經(jīng)濟(jì)區(qū)PRED系統(tǒng)的狀態(tài)特征，從而可用于模擬未來不同情景下的系統(tǒng)動態(tài)發(fā)展過程。為了解未來10年研究區(qū)人口、資源、環(huán)境與經(jīng)濟(jì)發(fā)展的協(xié)調(diào)途徑，根據(jù)研究區(qū)歷史數(shù)據(jù)和區(qū)域發(fā)展規(guī)劃，利用SD模型，仿真模擬三種不同假設(shè)條件下的發(fā)展情景。

情景A——原型發(fā)展模式：保持模型現(xiàn)有參數(shù)不變，按目前經(jīng)濟(jì)發(fā)展模式進(jìn)行仿真模擬。

情景B——經(jīng)濟(jì)高速發(fā)展模式：保持資源環(huán)境方面的參數(shù)不變，提高各產(chǎn)業(yè)增率、城鎮(zhèn)化增率與耕地面積增率等。

情景C——穩(wěn)定發(fā)展模式：參考《廣西北部灣經(jīng)濟(jì)區(qū)“十三五”發(fā)展規(guī)劃》，基于經(jīng)濟(jì)新常態(tài)的背景，穩(wěn)步發(fā)展經(jīng)濟(jì)的同時(shí)加大資源利用與環(huán)境保護(hù)力度，適當(dāng)調(diào)整三大產(chǎn)業(yè)增率；提高人口出生率、城鎮(zhèn)化增率、污水集中處理率、第二產(chǎn)業(yè)單位產(chǎn)值耗水量減率、每萬元地區(qū)生產(chǎn)總值消耗能源減率等指標(biāo)；降低單位面積耕地化肥使用量、城鎮(zhèn)人均生活垃圾排放增率等。

2.4 SD系統(tǒng)模擬運(yùn)行結(jié)果

選取2018年、2020年、2022年和2015年作為系統(tǒng)PRED仿真模擬的代表年份，對系統(tǒng)總?cè)丝凇⒌貐^(qū)生產(chǎn)總值、用水總量、能源消耗量、第一產(chǎn)業(yè)GDP、第二產(chǎn)業(yè)GDP、第三產(chǎn)業(yè)GDP等變量進(jìn)行模擬預(yù)測（表2）。

3 討論

3.1 原型發(fā)展模式

原型模式是基于研究區(qū)2008-2015歷史數(shù)據(jù)計(jì)算出來的系統(tǒng)參數(shù)，這一時(shí)期的經(jīng)濟(jì)處于快速增長的階段，各產(chǎn)業(yè)的GDP增速均保持在10%以上。按該模式發(fā)展，研究區(qū)2020年三大產(chǎn)業(yè)結(jié)構(gòu)比例為11.76：48.41：39.83，到2025年三大產(chǎn)業(yè)結(jié)構(gòu)的比例為8.94：53.04：38.02。從產(chǎn)業(yè)結(jié)構(gòu)的比例可以看出，研究區(qū)的發(fā)展模式以大力發(fā)展第二產(chǎn)業(yè)為主導(dǎo)，這與廣西北部灣開放開發(fā)初期打造的重工業(yè)體系是相一致的，如總裝機(jī)120萬千瓦的北?；痣姀S、240萬千瓦的欽州火電廠、240萬千瓦的防城港火電廠、中石油1000萬噸煉油廠、防城港千萬噸級鋼鐵項(xiàng)目、總裝機(jī)600萬千瓦的防城港核電項(xiàng)目、年產(chǎn)300萬噸重油瀝青項(xiàng)目等一大批重工業(yè)項(xiàng)目均在這一階段投產(chǎn)并為第二產(chǎn)業(yè)GDP增長作出很大的貢獻(xiàn)。然而，這些重工業(yè)的投產(chǎn)在增長GDP的同時(shí)也帶來了資源過度利用和環(huán)境污染的負(fù)面效應(yīng)，當(dāng)研究區(qū)經(jīng)濟(jì)發(fā)展達(dá)到一定程度，再過度注重第二產(chǎn)業(yè)的GDP增長勢必導(dǎo)致研究區(qū)PRED的不協(xié)調(diào)發(fā)展和產(chǎn)業(yè)結(jié)構(gòu)的不平衡。因此，這一發(fā)展模式不能順應(yīng)“十三五”以后的廣西北部灣發(fā)展趨勢。

3.2 高速發(fā)展模式

高速發(fā)展模式是在原型模式的基礎(chǔ)上繼續(xù)提高經(jīng)濟(jì)發(fā)展速度，并進(jìn)行產(chǎn)業(yè)結(jié)構(gòu)調(diào)整與優(yōu)化，保持第二產(chǎn)業(yè)發(fā)展優(yōu)勢，重點(diǎn)發(fā)展第三產(chǎn)業(yè)。按該模式發(fā)展，研究區(qū)2020年三大產(chǎn)業(yè)結(jié)構(gòu)比例為9.67：44.46：45.87，到2025年三大產(chǎn)業(yè)結(jié)構(gòu)的比例為6.76：46.91：46.33。地區(qū)生產(chǎn)總值將于2021年在原型模式同期地區(qū)生產(chǎn)總值的基礎(chǔ)上翻一番。在該模式下，經(jīng)濟(jì)得到高速發(fā)展，但資源利用率沒有相應(yīng)提高，也沒有采取更有力的措施應(yīng)對經(jīng)濟(jì)發(fā)展所帶來的環(huán)境污染，從而導(dǎo)致水資源等不可再生資源過度消耗，不符合廣西北部灣經(jīng)濟(jì)區(qū)“十三五”規(guī)劃的生態(tài)優(yōu)先，綠色發(fā)展的理念。

3.3 穩(wěn)定發(fā)展模式

依據(jù)《廣西北部灣經(jīng)濟(jì)區(qū)“十三五”規(guī)劃》發(fā)展目標(biāo)，到2020年，廣西北部灣經(jīng)濟(jì)區(qū)率先實(shí)現(xiàn)服務(wù)業(yè)為主導(dǎo)產(chǎn)業(yè)，地區(qū)生產(chǎn)總值達(dá)13300萬元，三次產(chǎn)業(yè)結(jié)構(gòu)比為11：42：47，常住人口城鎮(zhèn)化率達(dá)56%。情景C模擬結(jié)果顯示，到2020年，研究區(qū)地區(qū)生產(chǎn)總值為12961.7億元，三大產(chǎn)業(yè)結(jié)構(gòu)比例為11.38：42.07：46.55，常住人口城鎮(zhèn)化率達(dá)55.84%，基本符合規(guī)劃發(fā)展目標(biāo)。按該模式發(fā)展，到2025年，研究區(qū)地區(qū)生產(chǎn)總值達(dá)24200萬元，三次產(chǎn)業(yè)結(jié)構(gòu)比為8.58：43.74：47.68，常住人口城鎮(zhèn)化率達(dá)65%。同時(shí)，在發(fā)展經(jīng)濟(jì)的同時(shí)，資源利用和能源消耗得到了有效控制。由此可見，在三種情景模擬中，情景C的穩(wěn)定發(fā)展模式是保障廣西北部灣經(jīng)濟(jì)區(qū)PRED協(xié)調(diào)發(fā)展的最佳方案，有利于促進(jìn)研究區(qū)社會經(jīng)濟(jì)的可持續(xù)發(fā)展。

4 結(jié)論

本研究基于廣西北部灣2008-2015歷史數(shù)據(jù)，建立研究區(qū)PRED系統(tǒng)SD模型，預(yù)測三種不同情景下的PRED發(fā)展?fàn)顟B(tài)。通過分析比較三種情景下的各變量模擬仿真值，結(jié)合《廣西北部灣“十三五”發(fā)展規(guī)劃》目標(biāo)，選取情景C——穩(wěn)定發(fā)展模式作為廣西北部灣經(jīng)濟(jì)區(qū)PRED可持續(xù)發(fā)展最佳方案。

根據(jù)情景C的發(fā)展模式，廣西北部灣經(jīng)濟(jì)區(qū)在今后一段時(shí)間內(nèi)可采取以下發(fā)展對策：

（1）穩(wěn)步控制人口的增長速度，通過政策優(yōu)惠吸引區(qū)內(nèi)外更多高層次人才流入本區(qū)域，同時(shí)通過發(fā)展鄉(xiāng)鎮(zhèn)企業(yè)和服務(wù)業(yè)，盡可能減少本地人口的流出量，從而逐步降低凈流出人口數(shù)量。

（2）在提高城鎮(zhèn)化率的同時(shí)，確保耕地的質(zhì)量和數(shù)據(jù)，控制森林和綠地資源的開發(fā)力度，提高各類土地資源的生態(tài)效益。

（3）廣西北部灣經(jīng)濟(jì)區(qū)水資源豐富，人們的節(jié)水意識淡薄，從而造成了很多水資源的浪費(fèi)。今后需加強(qiáng)資源利用與環(huán)境保護(hù)宣傳力度，提高人們的資源利用與環(huán)境保護(hù)意識，將資源利用與環(huán)境保護(hù)貫穿到人們的日常生活和工作當(dāng)中。

（4）加強(qiáng)對工礦等排污企業(yè)的監(jiān)管，在關(guān)注其創(chuàng)造經(jīng)濟(jì)效益的同時(shí)，注重其環(huán)境保護(hù)與治理力度，以及資源開發(fā)利用效率，強(qiáng)化各企業(yè)創(chuàng)造的生態(tài)效益評價(jià)。

（5）加大土地整治力度和步伐，為農(nóng)業(yè)規(guī)模化、現(xiàn)代化、機(jī)械化創(chuàng)造更便利的條件；注重引進(jìn)高新技術(shù)加大工業(yè)企業(yè)創(chuàng)新力度，逐步削減傳統(tǒng)生產(chǎn)方式的工業(yè)企業(yè)；充分發(fā)揮廣西北部灣經(jīng)濟(jì)區(qū)濱海和沿邊旅游資源，大力發(fā)展旅游服務(wù)業(yè)，適當(dāng)加大第三產(chǎn)業(yè)投資力度，不斷調(diào)整與優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)。

本研究的SD模型是在假設(shè)保持當(dāng)前發(fā)展?fàn)顟B(tài)且參考研究區(qū)“十三五”規(guī)劃的基礎(chǔ)上進(jìn)行短期的仿真預(yù)測，在技術(shù)發(fā)展迅猛的今天，社會經(jīng)濟(jì)的發(fā)展模式也會隨著高新技術(shù)的引進(jìn)和改革的不斷深入而進(jìn)行不斷的調(diào)整和優(yōu)化，在“十三五”過后的一段時(shí)期，政府可能根據(jù)“十三五”的改革成果和實(shí)際情況重新調(diào)整經(jīng)濟(jì)發(fā)展導(dǎo)向，該模型也需要進(jìn)行適當(dāng)?shù)恼{(diào)整。

[參考文獻(xiàn)]

[1] 許聯(lián)芳，張建新，陳坤，等.基于SD模型的湖南省土地利用變化情景模擬[J].熱帶地理，2014（06）.

[2] 滕宇思.基于系統(tǒng)動力學(xué)的西安市土地綜合承載力研究[J].未來與發(fā)展，2016（01）.

[3] 王耕，魏遼生.基于系統(tǒng)動力學(xué)的大連市可持續(xù)發(fā)展模型模擬分析[J].海洋開發(fā)與管理，2015（02）.

[4] 鄭德鳳，臧正，張雨，等.基于新型城鎮(zhèn)化視角的區(qū)域PRED系統(tǒng)綜合評價(jià)——以大連為例[J].地理科學(xué)進(jìn)展，2014（03）.

[5] Alireaz Aslani， Petri Helo， Marja Naaranoja. Role of renewable enerty policies in energy dependency in Finland： System dynamics approach[J]. Applied Energy， 2013，113：758-765.

[6] Nastaran Ansari，Abbas Seifi. A system dynamics model for analyzing energy consumption and CO2 emission in Iranian cement industry under various production and export scenarios [J]. Energy Pollicy， 2013，58：75-89.

[7] Maria Niedertscheider， Tobias Kuemmerle， Daniel Muller， et al. Exploring the effects of drastic institutional and socio-economic changes on land system dynamics in Germany between 1883 and 2007[J].Global Environmental Change， 2014，28：98-108.

[8] Xi Xi， Kim Leng Poh. Using system dynamics for sustainable water resources management in Singapore[J]. Procedia Computer Science， 2013，16：157-166.

[9] Heather Karina Loyo， Cynthia Batcher， RN， et al. From Model to Action： Using a System Dynamics Model of Chronic Disease Risks to Align Community Action [J]. Health Promotion Practice， 2013，14（1）： 53–61.

[10] 劉承良，顏琪，羅靜.武漢城市圈經(jīng)濟(jì)資源環(huán)境耦合的系統(tǒng)動力學(xué)模擬[J].地理研究，2013（05）.

[11] 張子珩，濮勵(lì)杰，周秀慧.烏海市可持續(xù)發(fā)展的系統(tǒng)動力學(xué)模型仿真[J].干旱區(qū)資源與環(huán)境，2010（12）.

[12] 祝秀芝，李憲文，賈克敬，等.上海市土地綜合承載力的系統(tǒng)動力學(xué)研究[J].中國土地科學(xué)，2014（02）.