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