生態(tài)系統(tǒng)碳匯能力仿真模型結(jié)構(gòu)分析與模擬

1、1生態(tài)系統(tǒng)碳匯能力仿真模型結(jié)構(gòu)分析與模擬生態(tài)系統(tǒng)碳匯能力仿真模型結(jié)構(gòu)分析與模擬以以獼猴桃獼猴桃果園為例果園為例盧花（?？诮?jīng)濟(jì)學(xué)院財(cái)務(wù)會(huì)計(jì)學(xué)院，海南 ?？?71127） 摘摘 要要: : 論文運(yùn)用系統(tǒng)動(dòng)力學(xué)方法及原理，配合相關(guān)數(shù)據(jù)處理。依據(jù)系統(tǒng)動(dòng)力學(xué)模型結(jié)構(gòu)、因果關(guān)系分析的基礎(chǔ)上,建立了由七個(gè)子模型、52 個(gè)方程式組成的獼猴桃果園生態(tài)系統(tǒng)的碳匯能力的動(dòng)態(tài)模型。在計(jì)算機(jī)上模擬出 20112021 年獼猴桃果園生態(tài)系統(tǒng)的碳匯能力變化情況。模型通過(guò)不斷的修改、調(diào)試、檢驗(yàn)證明該模型結(jié)構(gòu)是可行的，而且也能夠用來(lái)支持系統(tǒng)仿真模擬。其研究結(jié)果表明：2011 年獼猴桃果園生態(tài)系統(tǒng)的凈碳匯實(shí)際值是4 262.41

2、 kgC/mu，系統(tǒng)仿真值是 4 260.45 kgC/mu，誤差僅有 1.96 kgC/mu。同時(shí)，利用系統(tǒng)仿真模型還預(yù)測(cè)出獼猴桃果園生態(tài)系統(tǒng)在未來(lái) 10 年碳匯能力的大小及變化趨勢(shì)。關(guān)鍵詞：關(guān)鍵詞：獼猴桃果園生態(tài)系統(tǒng)；模型結(jié)構(gòu)；因果分析；系統(tǒng)仿真模型 中圖分類(lèi)號(hào)：中圖分類(lèi)號(hào)：F30F30 文獻(xiàn)標(biāo)識(shí)碼：文獻(xiàn)標(biāo)識(shí)碼：A A 文章編號(hào)：文章編號(hào)：20172017（0101）- -Structure of Carbon Sink Capacity Simulation Model of Kiwifruit Orchard EcosystemLU Hua(School of Finance and

3、Accounting, Haikou College of Economics, Haikou, Hainan 571127)Abstract: The methods and principles of system dynamics are used to cooperate with data processing. Based on the analysis of systems dynamics model structure and causality, the dynamic model concerning of carbon sink capacity of kiwifrui

4、t orchard ecosystem composed of seven sub-models and 52 equations is established. The changes of carbon sink capacity of kiwifruit orchard ecosystem from 2011 to 2021 are simulated on a computer. The model is proved to be feasible by continuously rectifying, adjusting, and examining, and can be used

5、 to support the simulation of the system. The results show that in the year of 2011, the net carbon sink value of kiwifruit orchard ecosystem was 4262.41 kgC/mu, the simulation value was 4,260.45kgC/mu, with error only 1.96 kgC/mu. Furthermore, the system simulation model is used to predict the size

6、 and change trends of carbon sink capacity of kiwifruit orchard ecosystem in the next 10 years. Keywords: Ecosystem of kiwifruit orchard; model structure; causality analysis; system simulation model 一、引言引言 全球氣候變暖問(wèn)題，已成為當(dāng)今國(guó)內(nèi)外科學(xué)界關(guān)注的熱點(diǎn)話題之一。大部分的氣象專(zhuān)家認(rèn)為：溫室氣體，尤其是CO2的濃度劇增引起了氣候變暖。同時(shí)氣候變暖還可能引起一些極端惡劣天氣的發(fā)生，隨之而來(lái)將是

7、全球生態(tài)問(wèn)題，對(duì)人類(lèi)生產(chǎn)生活造成極大的威脅。因此，當(dāng)前碳排放、碳儲(chǔ)匯等相關(guān)問(wèn)題成為各國(guó)學(xué)術(shù)界的研究重點(diǎn)。與此同時(shí)，眾多估算碳匯的研究方法也層出不窮。以下我是目前國(guó)內(nèi)外對(duì)森林碳匯的一系列估算方法。 例如：樣地清查法，是指通過(guò)設(shè)立典型樣地， 準(zhǔn)確測(cè)定 森林生態(tài)系統(tǒng)中的植被、枯落物或土壤等碳庫(kù)的碳 儲(chǔ)量，并可通過(guò)連續(xù)觀測(cè)來(lái)獲知一定時(shí)期內(nèi)的儲(chǔ)量 變化情況的推算方法。歸納起來(lái)主要分為三種 方法，即平均生物量法，平均換算因子法和換算因 子連續(xù)函數(shù)法，這三種方法具有相同的數(shù)學(xué)推理方 法基礎(chǔ)，即都是在推算出生物量的基礎(chǔ)上再乘作者簡(jiǎn)介：作者簡(jiǎn)介：盧花（1982），女，江西九江人，?？诮?jīng)濟(jì)學(xué)財(cái)務(wù)會(huì)計(jì)學(xué)院講師，研究

8、方向：農(nóng)業(yè)經(jīng)濟(jì)。2以一 個(gè)換算系數(shù)求得碳儲(chǔ)量的方法，換算系數(shù)通常在0.440.55之間。獼猴桃果園是一個(gè)復(fù)雜的生態(tài)系統(tǒng)，具有復(fù)雜性、因果反饋性、動(dòng)態(tài)性、非線性等特點(diǎn)。其碳匯能力的大小不僅與獼猴桃果園生態(tài)系統(tǒng)本身密切相關(guān)，還與人口、管理投入、獼猴桃果種植面積以及農(nóng)民的生活水平等各種因素密切相關(guān)。根據(jù)獼猴桃果園生態(tài)系統(tǒng)的匯能力，針對(duì)其內(nèi)部結(jié)構(gòu)和運(yùn)行機(jī)制，建立獼猴桃果園生態(tài)系統(tǒng)碳匯能力預(yù)測(cè)仿真模型。該模型通過(guò)定性與定量，宏觀與微觀相結(jié)合，全面的對(duì)獼猴桃果園生態(tài)系統(tǒng)進(jìn)行的碳匯能力進(jìn)行估算，并通過(guò)系統(tǒng)仿真模型，宏觀、定量、動(dòng)態(tài)的預(yù)測(cè)獼猴桃果園生態(tài)系統(tǒng)未來(lái)10年碳匯能力大小。 二、二、模型結(jié)構(gòu)模型結(jié)構(gòu) 獼

9、猴桃果園生態(tài)系統(tǒng)碳匯能力模型設(shè)計(jì)的目的是進(jìn)行碳儲(chǔ)量的估算，影響排碳量的因子構(gòu)成了該系統(tǒng)的各要素。通過(guò)對(duì)農(nóng)業(yè)生態(tài)系統(tǒng)排碳、匯碳現(xiàn)狀的透視可知，人口因素模塊、管理投入模塊、獼猴桃收入模塊、獼猴桃產(chǎn)量模塊、獼猴桃能源消費(fèi)模塊、減排模塊、以及碳匯模塊。七個(gè)模塊構(gòu)成總模型結(jié)構(gòu)，具體見(jiàn)圖1.1。 圖 1.1 獼猴桃果園生態(tài)系統(tǒng)碳匯能力估算仿真模型結(jié)構(gòu)圖Fig 1.1 Structure of Estimating Carbon Sink Systematic Model of Kiwifruit Agro-ecosystems （一）模型子模塊的介紹 前面主要針對(duì)總體模型框架的敘述，下面是對(duì)模型中子模塊的

10、功能進(jìn)行介紹。 1.人口因素模塊 人口因素是決定著獼猴桃的需求量，進(jìn)而影響產(chǎn)獼猴桃所需能源消耗的傳統(tǒng)因素。較高的人口增長(zhǎng)會(huì)增加獼猴桃的需求量，從而增加能源的消費(fèi)需求量。因此而構(gòu)建人口因素模塊。 2.獼猴桃產(chǎn)量狀況模塊 這一模塊主要考慮在不同的獼猴桃果園生態(tài)系統(tǒng)中，獼猴桃產(chǎn)量由于采用不同的種植、管理、其能源消耗以及農(nóng)用材料的使用狀況也會(huì)不同，從而影響排碳和匯碳都會(huì)有所差別。 3.獼猴桃收入模塊 這一部分主要是考慮到收入這塊，獼猴桃收入增加，用于獼猴桃的投資就會(huì)增加，相應(yīng)的獼猴桃產(chǎn)量就增加。因此，相關(guān)的管理投入以及用于獼猴桃林固定資產(chǎn)投入等因素都會(huì)相應(yīng)發(fā)生變化，從而影響到獼猴桃果園生態(tài)系統(tǒng)的碳源/

11、匯發(fā)生相應(yīng)的變化。 4.獼猴桃能源消費(fèi)模塊3 在兩個(gè)不同的獼猴桃果園生態(tài)系統(tǒng)中，能源消耗狀況必定會(huì)有所不同的。綜合考慮到不同獼猴桃果園生態(tài)系統(tǒng)的獼猴桃收入和單產(chǎn)耗能，得出總能源消費(fèi)量。然后再考慮其他相關(guān)因素的影響。比如人口因素、管理投入、不同類(lèi)型能源、以及不同的排碳因子等因素，進(jìn)行這一模塊的構(gòu)建。 5.獼猴桃減排模塊 我國(guó)明確提出：到 2020 年我國(guó) CO2排放量要比 2005 年下降 4045。從這我們不難看出，不管農(nóng)業(yè)系統(tǒng)是碳匯還是碳源，我們都要做好增匯減排的工作。因此，建立這一模塊主要考慮人口因素、能源減排政策、能源減排比例、能源減排成本等因素的影響。 6.獼猴桃管理投入模塊 隨著科學(xué)

12、管理的進(jìn)步，對(duì)于能源的消耗會(huì)有一定程度的減少，新的管理投入會(huì)提高能源的使用效率，進(jìn)而影響到 CO2的排放。換句話說(shuō)可能會(huì)在原來(lái)的排放基礎(chǔ)上減少 CO2的排放。 7.獼猴桃碳匯模塊 本研究的核心就是碳匯，而增匯所涉及到的相關(guān)因子有，獼猴桃的種植面積、管理投入以及獼猴桃產(chǎn)業(yè)耗能和農(nóng)用材料的使用狀況等因素。為此，建立此模塊。 三、因果關(guān)系分析三、因果關(guān)系分析 在 SD 理論中進(jìn)行因果分析通常是由“果”進(jìn)行逆推，推到所定系統(tǒng)邊界，即所要研究的因素為止。圖 2.1 為獼猴桃果園生態(tài)系統(tǒng)碳匯能力模型簡(jiǎn)化因果關(guān)系圖。而獼猴桃果園生態(tài)系統(tǒng)的碳匯能力間接的是由能源及農(nóng)用材料消耗量決定，消耗量越大，排碳量也越大，

13、排碳量越大相應(yīng)的碳匯要抵消排碳，從而凈碳匯就越小。圖中的“-”號(hào)也正是表現(xiàn)了這一層意思。而決定獼猴桃果園生態(tài)系統(tǒng)能源及農(nóng)用材料消耗量的又是由獼猴桃種植面積、獼猴桃固定資產(chǎn)投資、人口數(shù)量等因素決定。 為便于讀者更好的理解因果關(guān)系中的反饋回路，在這里首先要介紹一下正（負(fù)）反饋回路的概念。 正反饋回路（正循環(huán)）：正反饋回路產(chǎn)生增長(zhǎng)、放大偏移并且加強(qiáng)變化。 負(fù)反饋回路（負(fù)循環(huán)）：負(fù)反饋回路尋求平衡、均衡和停滯。圖 2.1 獼猴桃果園生態(tài)系統(tǒng)碳匯能力模型簡(jiǎn)化因果關(guān)系圖4Fig 2.1 Simplified Cause and Effect Relationship Figure of Carbon Si

14、nk in kiwifruit Agro-ecosystems （一）因果反饋回路介紹 從圖中可知，有以下主要反饋回路：1、獼猴桃種植面積 獼猴桃產(chǎn)量 獼猴桃收入 獼猴桃固定資產(chǎn)投資 獼猴桃產(chǎn)量（正循環(huán)）2、獼猴桃種植面積 獼猴桃產(chǎn)量 獼猴桃耗能 獼猴桃收入 獼猴桃固定資產(chǎn)投資 獼猴桃產(chǎn)量（負(fù)循環(huán)）3、獼猴桃種植面積 獼猴桃產(chǎn)量 相關(guān)農(nóng)用材料使用量 獼猴桃收入 獼猴桃固定資產(chǎn)投資 獼猴桃產(chǎn)量（負(fù)循環(huán)）4、獼猴桃收入 管理投入 單產(chǎn)耗能 獼猴桃耗能 獼猴桃收入（正循環(huán)）5、獼猴桃收入 管理投入 農(nóng)用材料使用率 農(nóng)用材料排碳量 總排碳量（負(fù)循環(huán)） 6、獼猴桃收入 管理投入 單產(chǎn)耗能 獼猴桃耗能 耗

15、能排碳量 總排碳量 減排比例 減排量 減排成本 獼猴桃收入（正循環(huán)）7、獼猴桃產(chǎn)量 相關(guān)農(nóng)用材料使用量 農(nóng)用材料排碳量 總排碳量（正循環(huán)）8、獼猴桃產(chǎn)量 獼猴桃耗能 耗能排碳量 總排碳量 減排比例 減排量 減排成本 獼猴桃收入 獼猴桃固定資產(chǎn)投資 獼猴桃產(chǎn)量 （負(fù)循環(huán)）9、獼猴桃種植面積 獼猴桃林碳匯量 總排碳量 減排比例 減排量 減排成本 獼猴桃收入 獼猴桃固定資產(chǎn)投資 獼猴桃產(chǎn)量（正循環(huán)）10、獼猴桃種植面積 相關(guān)農(nóng)用廢棄物 廢棄物排碳量 總排碳量 減排比例 減排量 減排成本 獼猴桃收入 獼猴桃固定資產(chǎn)投資獼猴桃產(chǎn)量（負(fù)循環(huán)） 圖中還有幾條類(lèi)似反饋回路在這里就不一一敘述了，通過(guò)對(duì)系統(tǒng)各要素

16、的因果關(guān)系分析，可以明確了解在獼猴桃生產(chǎn)過(guò)程中的排碳和匯碳機(jī)制，還需指出的是系統(tǒng)中的各參數(shù)因子在不同時(shí)期會(huì)具有不同的地位和作用，這主要通過(guò)因子在整個(gè)系統(tǒng)中對(duì)其影響對(duì)象的強(qiáng)度進(jìn)行體現(xiàn)。 （二）因果反饋回路分析 在第一個(gè)正循環(huán)中，由于獼猴桃種植面積增加，使獼猴桃產(chǎn)量的增加，從而促進(jìn)了獼猴桃收入的增加，獼猴桃收入增加后，可用于獼猴桃固定資產(chǎn)投資的費(fèi)用亦相應(yīng)增加，而這一行為又會(huì)返回來(lái)刺激獼猴桃產(chǎn)量的增長(zhǎng)；那么獼猴桃產(chǎn)量是否就會(huì)這樣無(wú)上限的增加呢？我們?cè)賮?lái)看看第二個(gè)反饋回路；由于獼猴桃種植面積增加，獼猴桃產(chǎn)量也相應(yīng)增加，獼猴桃耗能也會(huì)增加，隨之就是生產(chǎn)費(fèi)用及相關(guān)費(fèi)用的增加，而這個(gè)增長(zhǎng)并不會(huì)增加獼猴桃的收

17、入，相反使收入降低了，而獼猴桃收入的減少后，使可用于獼猴桃固定資產(chǎn)投資的資金必然減少，最終導(dǎo)致獼猴桃的產(chǎn)量不會(huì)隨著再增加。同理第三個(gè)回路對(duì)獼猴桃產(chǎn)量起到了同樣的抑制作用。那么獼猴桃產(chǎn)量是增、是減還是保持不變呢?這就由這三個(gè)回路共同作用決定了。其它循環(huán)都是同理，在這里就不一一介紹。 而人口數(shù)量是這個(gè)小系統(tǒng)的外部輸入條件。隨著人口數(shù)量增加，對(duì)獼猴桃的需求量增5加。而人口數(shù)量的下降，又會(huì)導(dǎo)致獼猴桃的需求量降低。當(dāng)然，影響需求的因素不止這一個(gè)，但在不考慮其它因素條件的情況下這種關(guān)系是成立的。 從系統(tǒng)因果分析中我們可以看出，能源消耗以及獼猴桃果園的一些基礎(chǔ)設(shè)施和農(nóng)用輔助工具所需的農(nóng)用材料消耗量因子，在系

18、統(tǒng)中是決定排碳量大小的。它的數(shù)量變化關(guān)系到估算內(nèi)部對(duì)象排碳量的變化(即能源消耗及農(nóng)用材料消費(fèi)量增加，排碳量相應(yīng)上升，反之下降)，它是一個(gè)“累積性數(shù)量”。而影響能源消耗及農(nóng)用材料使用量的主要因子，是獼猴桃產(chǎn)量受價(jià)格、成本、管理、獼猴桃固定資產(chǎn)投資等信息流的制約與調(diào)節(jié)。該系統(tǒng)中的四條負(fù)反饋回路具有平抑、穩(wěn)定排碳量激烈變化的主要作用。當(dāng)我們引入人口因素及減排比例影響因子的時(shí)候，就會(huì)使整個(gè)排碳系統(tǒng)在這些反饋回路的作用下，相互影響、相互藕合、并且使其處于動(dòng)態(tài)平衡中。 在實(shí)際系統(tǒng)中，影響能源消耗和農(nóng)用材料消耗量的因子還有其消耗強(qiáng)度。而消耗強(qiáng)度又受管理水平所制約，其能源消耗量又與獼猴桃產(chǎn)量有直接聯(lián)系。在估算

19、對(duì)象內(nèi)部，雖然獼猴桃產(chǎn)量從理論上講應(yīng)主要表現(xiàn)為供應(yīng)與需求之間的關(guān)系變動(dòng)，但這一結(jié)果在實(shí)際中并不由需求直接決定。影響需求的因素除上述人口因素外。還包括經(jīng)濟(jì)增長(zhǎng)速度、獼猴桃價(jià)格等一系列因素。而這些因素的背后各自又有著許多相關(guān)因子，并且各因子間存在著錯(cuò)綜復(fù)雜的聯(lián)系和相互作用的關(guān)系。 四、四、獼猴桃獼猴桃果園生態(tài)系統(tǒng)碳匯能力仿真模型的應(yīng)用果園生態(tài)系統(tǒng)碳匯能力仿真模型的應(yīng)用 （一）模型條件假設(shè) 通過(guò)對(duì)獼猴桃果園生態(tài)系統(tǒng)碳匯能力的分析，鑒于條件的約束，本研究假設(shè)能源有效燃燒所含碳量最終都以 CO2氣體形式排放。獼猴桃果園能正常存活，不遭受?chē)?yán)重的病蟲(chóng)害以及自然災(zāi)害所成的獼猴桃樹(shù)不正常死亡或過(guò)早死亡的假設(shè)條件

20、下。 （二）模型檢驗(yàn) 1.表達(dá)準(zhǔn)確性 系統(tǒng)是結(jié)構(gòu)和功能的統(tǒng)一體。在建模過(guò)程中，不僅要考慮系統(tǒng)的功能和行為，而且還要考慮系統(tǒng)的結(jié)構(gòu)是否真實(shí)的反映實(shí)際情況。本文通過(guò)對(duì)人口、經(jīng)濟(jì)、能源以及不同獼猴桃果園生態(tài)系統(tǒng)，環(huán)境狀況的全面分析和準(zhǔn)確把握。同時(shí)借鑒相關(guān)模型，在應(yīng)用中具體情況具體分析，最后通過(guò)不斷的調(diào)試、修改和運(yùn)行過(guò)程中建立此模型結(jié)構(gòu)。 在 Stella 軟件中運(yùn)行已表達(dá)的數(shù)學(xué)模型，經(jīng)過(guò) Stella 軟件的編譯檢錯(cuò)以及跟蹤功能檢驗(yàn)，驗(yàn)證了模型表達(dá)正確性和量綱的一致性。同時(shí)研究中的模型界限通過(guò)檢驗(yàn)，證明該模型結(jié)構(gòu)是可行的。 2.模型有效性 將調(diào)研的歷史統(tǒng)計(jì)數(shù)據(jù)與模型模擬的數(shù)據(jù)進(jìn)行比較，就可以檢驗(yàn)?zāi)Ｐ?/p>

21、的有效性。本模型主要從獼猴桃果園生態(tài)系統(tǒng)的模型仿真凈碳匯值與實(shí)際生活中獼猴桃果園生態(tài)系統(tǒng)的碳匯值，來(lái)進(jìn)行檢驗(yàn)擬合程度。如 4.1 表：表 4.1 凈碳匯總量的實(shí)際值與仿真值的對(duì)比Table 4.1 The comparison of the net carbon aggregate amount of the actual value and simulation value年份實(shí)際值（kgC/667m2)仿真值（kgC/667 m2)誤差20071,665.671,664.641.0320082,190.922,193.00-2.0820092,826.162,827.36-1.202010

22、3,569.293,568.181.1120114,262.414,260.451.96 3.靈敏性分析 6 靈敏度分析主要是結(jié)構(gòu)靈敏度和參數(shù)靈敏度這兩個(gè)方面，通過(guò)獼猴桃果園生態(tài)系統(tǒng)碳匯能力進(jìn)行檢驗(yàn)，結(jié)果表明： 就參數(shù)靈敏度來(lái)說(shuō)，模型中大部分參數(shù)值的改變，只會(huì)讓運(yùn)行結(jié)果的趨勢(shì)發(fā)生一定的變化，而不會(huì)導(dǎo)致系統(tǒng)行為發(fā)生根本性的變化。這些說(shuō)明了模型對(duì)大多參數(shù)敏感性不是很強(qiáng)烈，對(duì)參數(shù)的變化是不敏感的，同時(shí)也是有效的。而對(duì)于模型的總體結(jié)構(gòu)，而言在測(cè)試的過(guò)程中發(fā)現(xiàn)，模型的正反饋結(jié)構(gòu)靈敏度是高于負(fù)反饋靈敏度的。根據(jù)系統(tǒng)動(dòng)力學(xué)理論知識(shí)表明出現(xiàn)這種結(jié)果，說(shuō)明模型的結(jié)構(gòu)是合理可行的。 綜上所述表明了模型經(jīng)過(guò)了表達(dá)正

23、確性、模型有效性以及靈敏度分析三方面的驗(yàn)證，因此確定了本研究的模型是可以用來(lái)反映真實(shí)系統(tǒng)的，而且也能夠用來(lái)支持系統(tǒng)仿真模擬。 （三）傳統(tǒng)散戶(hù)種植模式碳匯能力仿真模型 1.系統(tǒng)仿真流程圖 見(jiàn)附錄（1） 2.模型方程式 見(jiàn)附錄（2） 3.參數(shù)說(shuō)明見(jiàn)附錄（3） 4.仿真結(jié)果 在上述假設(shè)及約束下建立模型并運(yùn)行，得出獼猴桃果園生態(tài)系統(tǒng)碳匯能力的估算，結(jié)果如 4.2 表：表 4.2 獼猴桃果園生態(tài)系統(tǒng)碳匯能力仿真估算表Table 4.2 kiwifruit orchard ecosystem carbon sink capacity simulation estimates in table年份凈碳匯（k

24、gC/667m2)20114,262.4520125,102.6720135,870.0820146,749.6120157,553.6520168,475.5520179,381.56201810,302.47201911,352.04202012,324.41202113,433.38 五、結(jié)論五、結(jié)論 本文對(duì)獼猴桃果園生態(tài)系統(tǒng)碳匯能力大小進(jìn)行了研究，采用系統(tǒng)動(dòng)力學(xué)的方法構(gòu)建模型，并得出仿真結(jié)果。這種估算方法比較準(zhǔn)確的反應(yīng)了獼猴桃果園生態(tài)系統(tǒng)碳匯能力的大小，具有一定的實(shí)用價(jià)值。同時(shí)預(yù)測(cè)出，獼猴桃果園生態(tài)系統(tǒng)未來(lái) 10 年的碳匯能力。結(jié)果表明： 1、利用系統(tǒng)動(dòng)力學(xué)方法及原理，構(gòu)建了一個(gè)由七個(gè)

25、子模塊：人口因素模塊、獼猴桃產(chǎn)量狀況模塊、獼猴桃收入模塊、獼猴桃能源消費(fèi)模塊、獼猴桃減排模塊、獼猴桃管理投入模塊、獼猴桃碳匯模塊組成的仿真模型，模擬未來(lái)碳匯能力。通過(guò)不斷的修改、調(diào)試、檢7驗(yàn)證明該模型結(jié)構(gòu)是可行的，而且也能夠用來(lái)支持系統(tǒng)仿真模擬。 2、仿真模擬數(shù)據(jù)表明，在未來(lái) 10 年獼猴桃果園生態(tài)系統(tǒng)的碳匯增長(zhǎng)潛力都是呈漸步上升趨勢(shì)。這是系統(tǒng)仿真數(shù)據(jù)，前提是在模型假設(shè)條件都成立的情況下。 3、對(duì)未來(lái)農(nóng)業(yè)里園生態(tài)系統(tǒng)中，相關(guān)的農(nóng)業(yè)碳審計(jì)工作提供強(qiáng)有力的理論依據(jù)。另一方面為制定碳稅政策、擴(kuò)大碳匯范圍、建立碳交易市場(chǎng)及與其他國(guó)家進(jìn)行碳貿(mào)易也具有重要意義。 六、附錄六、附錄 （一）系統(tǒng)仿真流程圖nc

26、slcslrkslnswlswlmhtzzmjmhtclmhthnny jxhnny y shnmhthnpttsy lctshnmny ptctshnptmny jpzcmny zcjpljplzthzptny jxhnblny y shnblny sy lmhtxzlxgny fqwlmhtzty jfxgny fqwptny sy lptdwmjmhtlztlngzzlngztzlgdzcztnzzlmhtslmhtslnzzlngltrbhgltrngltrlmny jgbhmny jgbhzmny jgzlmny jgtjcctshnbljthsxsnzjlcsnzjhncsnptc

27、snzjhnblsnhntjny jxhntjny y shntjctshntj （二）模式模型方程式gdzczt(t)=gdzczt(t - dt) + (ngzzl) * dtINIT gdzczt=6276.44*0.44836INFLOWS:ngzzl=mhtsl*ngztzlgltr(t)=gltr(t - dt) + (ngltrbh) * dtINIT gltr=0.7INFLOWS:ngltrbh=gltr*ngltrl/100mhtsl(t)=mhtsl(t - dt) + (nzzl) * dtINIT mhtsl=32508INFLOWS:nzzl=mhtsl*mhtsl

28、nzzl/100mnyjgbh(t)=mnyjgbh(t - dt) + (mnyjgbhz) * dtINIT mnyjgbh=1INFLOWS:mnyjgbhz=if(mnyjgbh0)then(mnyjgzl+mnyjgtjc)/100 else(0)rksl(t)=rksl(t - dt) + (ncsl - nswl) * dtINIT rksl=2082INFLOWS:ncsl=rksl*csl/1000OUTFLOWS:nswl=rksl*swl/1000csnpt=csnzjhn*2.493csnzjhn=csnzjhnbl*sxsnzjlcsnzjhnbl=gdzczt*0.

29、23+gltr*0.5+snhntjctshn=tsyl*ctshnblctshnbl=gltr*0.5+gdzczt*0.45+mnyjgbh+ctshntjctshnpt=ctshn*2.493ctshntj=TIMEdwmjmhtlztl=(mhtztyjf*0.5*667)+(1*3.446*65+7143.53)*100jpl=mnyzcjpljth=(zth-zpt)/100mhtcl=mhtzzmj*2500mhthn=(mhtcl*nyjxhn*2.63)+(mhtcl*nyyshn*0.785)mhthnpt=mhthn*100mhtxzl=mhtzzmj*1000mhtzt

30、yjf=mhtxzl*0.5mhtzzmj=rksl*100mnyjgtjc=TIMEmnyjgzl=mhtsl/10000mnypt=mhthnpt+ctshnpt*0.97mnyzcjpl=mnypt*(1-mnyjpzc)nyjxhn=5000*0.785*nyjxhnblnyjxhnbl=gdzczt*0.043+gltr*0.5+nyjxhntjnyjxhntj=TIMEnysyl=mhtzzmj*46.5nysylpt=nysyl*2.493nyyshn=nyyshnbl*0.785nyyshnbl=gdzczt*0.043+gltr*0.5+nyyshntjnyyshntj=TI

31、MEsnhntj=TIMEsxsnzjl=mhtzzmj*75tsyl=mhtzzmj*1565.25xgnyfqwl=mhtcl*44.389xgnyfqwpt=xgnyfqwl*2.493zpt=csnpt+mnypt+nysylpt+xgnyfqwpt-jplzth=dwmjmhtlztl+jpl*0.9378csl=GRAPH(TIME)(2011, 0.083), (2012, 0.09), (2013, 0.081), (2014, 0.09), (2015, 0.078), (2016, 0.086), (2017, 0.0715), (2018, 0.079), (2019,

32、0.087), (2020, 0.081), (2021, 0.081)mhtslnzzl=GRAPH(TIME)(2011, 1.45), (2012, 1.50), (2013, 1.35), (2014, 1.75), (2015, 1.55), (2016, 1.65), (2017, 2.10), (2018, 1.70), (2019, 2.35), (2020, 1.75), (2021, 2.40)mnyjpzc=GRAPH(TIME)(2011, 0.1), (2012, 0.1), (2013, 0.1), (2014, 0.1), (2015, 0.1), (2016,

33、0.1), (2017, 0.1), (2018, 0.1), (2019, 0.1), (2020, 0.1), (2021, 0.1)ngltrl=GRAPH(mhtsl/10000)(0.00, 0.535), (10.0, 0.505), (20.0, 0.54), (30.0, 0.505), (40.0, 0.545), (50.0, 0.505), (60.0, 0.55), (70.0, 0.5), (80.0, 0.53), (90.0, 0.475), (100, 0.525)ngztzl=GRAPH(TIME)(2011, 3.00), (2012, 2.70), (20

34、13, 3.05), (2014, 2.75), (2015, 3.10), (2016, 3.00), (2017, 3.00), (2018, 3.35), (2019, 3.05), (2020, 3.40), (2021, 3.05)swl=GRAPH(TIME)(2011, 0.053), (2012, 0.057), (2013, 0.062), (2014, 0.059), (2015, 0.054), (2016, 0.068), (2017, 0.058), (2018, 0.075), (2019, 0.052), (2020, 0.059), (2021, 0.054) （三）模型參數(shù)說(shuō)明1Rksl 人口數(shù)量2Csl 出生率3Ncsl 年人口出生率4Swl 死亡率5Nswl 年人口死亡率6Mhtsl 獼猴桃收入7Mhtslnzzl 獼猴桃收入年增長(zhǎng)率8Nzzl 獼猴桃收入年增長(zhǎng)量9Gdzczt 固定資產(chǎn)總投資10Ngztzl 年固定資產(chǎn)投資率11Ngzzl 年固定資產(chǎn)增加值12Gltr 管理投入13Ngltrl 年管理投入率14Ngltrbh