摘要分析了林業(yè)發(fā)展的重要性、“低碳經(jīng)濟”和“森林碳匯”的主要性能指標(biāo)及概念，指出林業(yè)可持續(xù)發(fā)展的重要性，介紹數(shù)學(xué)生態(tài)學(xué)的基礎(chǔ)知識，并建立林業(yè)資源的開發(fā)與管理的數(shù)學(xué)模型，在此基礎(chǔ)上，考慮到森林數(shù)量由于自然環(huán)境的惡化，出現(xiàn)因火災(zāi)而發(fā)生的數(shù)量的突變。根據(jù)森林的密度制約型及單種群生物，建立Logistic隨機模型。
　　關(guān)鍵詞現(xiàn)代林業(yè)；數(shù)學(xué)生態(tài)學(xué)；可持續(xù)發(fā)展；密度制約型；隨機模型
　　生物數(shù)學(xué)是一門年輕的邊緣學(xué)科，近幾十年來發(fā)展尤為迅速，而數(shù)學(xué)生態(tài)學(xué)作為其中最為基礎(chǔ)的分支，相對來說發(fā)展得比較早，也比較成熟。它用數(shù)學(xué)模型來描述生物生存與環(huán)境的關(guān)系，并利用數(shù)學(xué)方法進行研究，以使一些生態(tài)現(xiàn)象得到解釋和控制?，F(xiàn)在數(shù)學(xué)已經(jīng)滲透到生態(tài)學(xué)研究的各領(lǐng)域，在種群生態(tài)學(xué)（包括種群的數(shù)量動態(tài)、空間分布、抽樣理論等）和生態(tài)系統(tǒng)的物流、能流分析中，以及在農(nóng)業(yè)生態(tài)系統(tǒng)的研究中，運用了大量數(shù)學(xué)方法，準(zhǔn)確地刻畫各系統(tǒng)的內(nèi)在聯(lián)系，經(jīng)分析和決策，做出最優(yōu)管理和決策支持系統(tǒng)[1]。而在林業(yè)可持續(xù)發(fā)展的開發(fā)與管理中，既要保護環(huán)境，也要保持一定的森林覆蓋率，又要使林業(yè)具有最大的經(jīng)濟效益。因此，數(shù)學(xué)生態(tài)學(xué)在林業(yè)的可持續(xù)發(fā)展中有了用武之地[2]。
　　林業(yè)屬于可再生能源，并且具有一定的密度制約。因此，該文利用數(shù)學(xué)生態(tài)學(xué)的基礎(chǔ)知識建立了林業(yè)資源的開發(fā)與管理的數(shù)學(xué)模型，并在此基礎(chǔ)上，考慮到森林數(shù)量由于自然環(huán)境的惡化，出現(xiàn)因火災(zāi)而發(fā)生數(shù)量的突變。根據(jù)森林的密度制約型及單種群生物，建立Logistic隨機模型，為林業(yè)的發(fā)展提供一個可能的新思路。
　　1林業(yè)的可持續(xù)發(fā)展和數(shù)學(xué)生態(tài)學(xué)
　　1.1林業(yè)可持續(xù)發(fā)展的重要性
　　林業(yè)在國民經(jīng)濟的發(fā)展中起著重要的作用。2010年3月，北京召開的“兩會”中，林業(yè)建設(shè)成為國家發(fā)展大計中的亮點，生態(tài)理念已滲透到“兩會”的方方面面，許多省份將林業(yè)單獨列出來，成為與經(jīng)濟建設(shè)、民生保障等并列的議題。政府工作報告中也明確提出，2010年新增造林面積不低于592萬hm2，這既體現(xiàn)政府對林業(yè)建設(shè)的高度重視和極大關(guān)注，也為林業(yè)發(fā)展、生態(tài)建設(shè)指明前進的方向[3]。
　　高消耗、重污染、低效益的傳統(tǒng)發(fā)展模式，在實現(xiàn)GDP數(shù)字增長的同時，帶來嚴重的環(huán)境污染和生態(tài)破壞，威脅著廣大人民群眾的健康和生存，制約人的全面發(fā)展與小康社會建設(shè)。“低碳經(jīng)濟”首次成為“兩會”的十大熱點。2009年11月，中國公布到2020年的減排目標(biāo)：單位國內(nèi)生產(chǎn)總值二氧化碳排放比2005年下降40%~45%。發(fā)展低碳經(jīng)濟的目標(biāo)就在于將環(huán)境保護作為重要的民生問題，在發(fā)展中加以統(tǒng)籌解決。通過開展大規(guī)模的植樹造林，可有效保護生態(tài)環(huán)境，減少水土流失。據(jù)科學(xué)測定，森林每生長1 m3的木材，約可有效吸收二氧化碳1.83 t。由此可見，轉(zhuǎn)變中國經(jīng)濟發(fā)展方式、促進經(jīng)濟結(jié)構(gòu)調(diào)整正在成為全社會的共識。
　　根據(jù)國家林業(yè)局2009年發(fā)布的《應(yīng)對氣候變化林業(yè)行動計劃》，我國已充分認識并高度重視森林在減緩氣候變化中的獨特作用，將采取22項林業(yè)減緩和適應(yīng)氣候變化行動，計劃到2020年，全國年均造林育林面積500萬hm2以上，森林覆蓋率增加到23%，森林蓄積量達到140億 m3，森林碳匯能力進一步提高；到2050年，比2020年凈增森林面積4 700萬hm2，森林覆蓋率達到并穩(wěn)定在26%以上，森林碳匯能力保持相對穩(wěn)定。所謂“森林碳匯”，是指森林植物能夠吸收大氣中的二氧化碳，并將其固定在植被或土壤中，從而減少二氧化碳在空氣中的濃度。總之，林業(yè)是關(guān)系中華民族生存發(fā)展根基、全球生態(tài)安全的偉大事業(yè)。
　　1.2 數(shù)學(xué)生態(tài)學(xué)
　　20世紀是生態(tài)學(xué)成熟和發(fā)展的時期。數(shù)學(xué)生態(tài)學(xué)作為生態(tài)學(xué)中的分支學(xué)科，在20世紀20—40年代經(jīng)歷了成長期。60年代以后，有2個重要因素對生態(tài)模型的發(fā)展起至關(guān)重要的作用：一個是電子計算機技術(shù)的快速發(fā)展，另一個是工業(yè)化高速的發(fā)展。人們?nèi)找嬲J識到保護生態(tài)環(huán)境的重要性，對環(huán)境治理、資源合理開發(fā)、能源持續(xù)利用越來越關(guān)心。面對這些復(fù)雜生態(tài)系統(tǒng)的研究，只有借助于系統(tǒng)分析及計算機模擬，才能解決諸如預(yù)測系統(tǒng)的行為及提出治理的最佳方案等問題。在國際生物學(xué)計劃（IBP）的促進下，20多年來經(jīng)過Watt K E、Van Dyne G M、Holling C S、Odum H T、Patten B C等生態(tài)學(xué)家的創(chuàng)造性的工作，形成一門新興的生態(tài)學(xué)分支學(xué)科——系統(tǒng)生態(tài)學(xué)。這時系統(tǒng)模型的應(yīng)用已緊密聯(lián)系到實際的生態(tài)系統(tǒng)，例如植被演替模型、害蟲綜合管理模型、漁業(yè)開發(fā)管理模型、全球氣候變化模型、有毒物質(zhì)在生態(tài)系統(tǒng)中的遷移和累積的模型、人口模型和世界模型等[4]。
　　我國數(shù)學(xué)生態(tài)學(xué)的發(fā)展也緊跟國際的步伐。1962年馬世駿先生在中國科學(xué)院動物研究所建立我國第一個數(shù)學(xué)生態(tài)研究組，并開展這方面的研究工作，于1965年發(fā)表有關(guān)我國電子計算機應(yīng)用于生態(tài)學(xué)的第一篇論文，在國際數(shù)學(xué)生態(tài)學(xué)發(fā)展的舞臺上，引起強烈的反響[5]。我國數(shù)學(xué)生態(tài)學(xué)的發(fā)展與國際相比，更緊密聯(lián)系于國民經(jīng)濟的發(fā)展與實際應(yīng)用。中國數(shù)學(xué)生態(tài)學(xué)的發(fā)展源于昆蟲生態(tài)學(xué)，其后很快涉及到植物生態(tài)學(xué)、動物生態(tài)學(xué)、環(huán)境評價與生態(tài)規(guī)劃、可更新資源管理及保護生物學(xué)等領(lǐng)域。
　　2林業(yè)資源開發(fā)與管理的數(shù)學(xué)模型
　　為了能長期地利用各種各樣的生物資源，人類就必須科學(xué)地進行開發(fā)與管理，不僅要考慮到當(dāng)前的高產(chǎn)，而且應(yīng)考慮到保護生態(tài)平衡，以保證長期的高產(chǎn)；不僅考慮到產(chǎn)量的高低，還應(yīng)考慮到投入產(chǎn)出所獲得的經(jīng)濟利潤。生物經(jīng)濟學(xué)模型就是研究生物資源的開發(fā)與管理。在過去的幾十年中，這類模型引起生態(tài)學(xué)家和經(jīng)濟學(xué)專家們的濃厚興趣，比較系統(tǒng)地從事這方面研究的是Clark，他于1976年出版了這方面的專著[6]。此后的若干年，這個領(lǐng)域的研究有了長足的進步，提出了一系列數(shù)學(xué)模型，簡單介紹如下。
　　對于可以用單變量描述的資源存儲，用差分方程描述資源改變量：
　　Xt+1-Xt=F（Xt）（1）
　　式中，Xt為在t時刻資源的生物量，函數(shù)F（Xt）反映了資源的凈增長以及環(huán)境負載能力等多種因素的效應(yīng)。如果用Logistic方程來描述這種效應(yīng)，公式為：
　　Xt+1-Xt=γXt（1-Xt/K）（2）
　　在沒有人為收獲的情況下，種群將趨于平衡點K，趨于平衡點的方式取決于內(nèi)稟增長率γ的不同取值范圍，可分為漸近趨于K、衰減振蕩到K、無阻尼振蕩以及混沌的情況。就是說，對差分方程形式的Logistic資源動態(tài)方程，即使沒有人為地開發(fā)利用，種群也會呈現(xiàn)出許多負載的變化形式。為了使林業(yè)的發(fā)展具有經(jīng)濟利益，當(dāng)然就需要有人為收獲的影響。因此，可用收獲函數(shù)來描述：
　　Yt=H（Et，Xt）（3）
　　式中，Et為人為投入的努力，Xt為現(xiàn)存的生物量。
　　對林業(yè)來講，若每單位面積的投入與生產(chǎn)量成正比，則有：
　　Yt=qEtXt（4）
　　式中，q為砍伐系數(shù)。這時的資源方程可寫成：
　　Xt+1-Xt=F（Xt）-H（Et，Xt）（5）
　　當(dāng)砍伐者希望收獲區(qū)域內(nèi)的資源能長期維持下去，就是說資源改變不能有過度的大起大落，而應(yīng)該維持在某個平衡水平，即：
　　Xt=X*，Et=E*，Yt=Y*
對于Logistic增長模型，以及單位投入量與產(chǎn)量成正比的砍伐，有Gordon-schaefer模型[7]：
　　qE*X*=rX*（1-X*/K）（6）
　　把方程寫成X作為E的函數(shù)，則有：
　　Y=Y（E*）=qKE*（1-qE*/r）（7）
　　對開發(fā)可供公共享用的資源，以及此類問題的最優(yōu)管理策略，已經(jīng)有了較深入的研究，其中的部分結(jié)果已經(jīng)在漁業(yè)管理中得到應(yīng)用，但是在林業(yè)方面應(yīng)用的還不多。
　　3數(shù)學(xué)生態(tài)學(xué)隨機模型
　　在已有的數(shù)學(xué)模型中，通常把種群之間的影響及環(huán)境對種群的影響歸結(jié)到模型的分布參數(shù)中，使問題得以簡化為微分方程和差分方程。從局部來看，或者在一段相對較短的時期內(nèi)，用上述模型來描述生物種群的發(fā)展和變化是可行的。但從整體看來，或者說在一段相對較長的時期內(nèi)，用上述模型描述生物種群的變化就存在明顯的缺陷。例如某種昆蟲當(dāng)環(huán)境突然發(fā)生變化（例如溫度突然下降或突然上升），其數(shù)量也可能會發(fā)生突然的變化。對于林業(yè)來說，則是人為或天然火災(zāi)，會使森林數(shù)量急劇下降；加之近年來自然環(huán)境的惡化，出現(xiàn)火災(zāi)的頻率上升，所以此種情況必須考慮在內(nèi)，即在環(huán)境發(fā)生變化時考慮森林數(shù)量的突變。具有以下特性：一是這種突變發(fā)生的時刻是不確定的；二是突變的強度一般來說是不確定的。因此，根據(jù)森林的密度制約型及單種群生物，建立Logistic隨機模型如下[8]。
　　密度制約的Logistic模型由以下微分方程決定：
　　8）
　　其解曲線為：
　　N（＜1（13）
　　4結(jié)論
　　我國目前對生態(tài)林的資本投入相對較多，是因為過去缺乏環(huán)境意識，只為了暫時的經(jīng)濟利益而亂砍濫伐，導(dǎo)致森林被大量破壞，現(xiàn)在為了恢復(fù)林業(yè)資本投入較大。但這只是一個過渡時期，從長期均衡的角度來看，在生態(tài)林的資本投入上是減少，并可以獲得很大的經(jīng)濟利益。
　　當(dāng)前林業(yè)建設(shè)依然是按照工程類型分類建設(shè)的模式發(fā)展，在具體建設(shè)中，往往因缺乏系統(tǒng)規(guī)劃，林地的增長模式是以數(shù)量增長而非效益增長的模式。就可持續(xù)發(fā)展的目標(biāo)而言，當(dāng)前林業(yè)建設(shè)存在一定的問題，對林業(yè)的持續(xù)性發(fā)展帶來巨大的隱患。而深入解決當(dāng)前的問題，需要從系統(tǒng)的角度入手，需要建立良性、健康、具有生態(tài)安全格局的建設(shè)模式。該文為我國林業(yè)的可持續(xù)發(fā)展提供了一個可行的理論依據(jù)。
　　5參考文獻
　　[1] 劉利民，吳素文.生態(tài)學(xué)中的數(shù)學(xué)模型[J].沈陽農(nóng)業(yè)大學(xué)學(xué)報，2000， 31（3）：295-299.
　　[2] 尹玉顰，聶華.關(guān)于我國林業(yè)資源最優(yōu)配置的經(jīng)濟學(xué)思考[J].經(jīng)濟工作導(dǎo)刊，2003，4（8）：17-18.
　　[3] 田新程.低碳經(jīng)濟趨熱生態(tài)文明升溫——兩會林業(yè)特點透視[J].中國林業(yè)，2010（6）：4-5.
　　[4] 李典謨，馬祖飛.展望數(shù)學(xué)生態(tài)學(xué)與生態(tài)模型的未來[J].生態(tài)學(xué)報，2000，20（6）：1083-1089.
　　[5] 藍仲雄，李典謨.數(shù)學(xué)生態(tài)學(xué)進展[M].成都：成都科技大學(xué)出版社，1994.
　　[6] CLARK C W.Mathematical Bioeconomics[M].New York：Wiley-intersci-ence，1976：1-32.
　　[7] 馬世駿.現(xiàn)代生態(tài)學(xué)透視[M].北京：科學(xué)出版社，1990：202-211.
　　[8] 侯振挺，劉再明.數(shù)學(xué)生態(tài)學(xué)隨機模型[J].生物數(shù)學(xué)學(xué)報，2000，15（3）：301-307.