0105064森林面積資源的計算機仿真及其應用(doc)
森林面積資源的計算機仿真及其應用* 張冰峰** （北京林業(yè)大學信息學院，100083） 摘要 本文論述了森林面積資源的計算機仿真理論、方法及其應用。森林面積資源仿真是 進行森林面積資源現(xiàn)代化管理的理論依據和重要的系統(tǒng)方法之一，具有廣泛的應用 。 關鍵詞 森林面積資源 計算機仿真 應用 1 引言 仿真技術已成為研究各種復雜系統(tǒng)的重要工具，在生態(tài)環(huán)境和林業(yè)上的應用已成為研 究的熱點之一。森林面積資源的計算機仿真就是對描述森林面積資源的數學模型在計算 機上進行各種模擬、試驗和分析，以便描述森林面積資源的動態(tài)變化規(guī)律和運動軌跡。 目前國內外林業(yè)上進行的森林和森林資源的計算機仿真所使用的數學模型不是統(tǒng)計模型 就是經驗模型，并且多為常微分方程，這些模型都不是根據森林資源的動態(tài)變化規(guī)律和 生滅過程而推導出來的數理方程，因而它們不能精確的對森林資源系統(tǒng)進行描述。本文 所使用的模型是根據森林面積資源的動態(tài)變化規(guī)律和生滅過程而推導出來的數理方程， 是偏微分方程，它能對森林面積資源的動態(tài)變化規(guī)律和生滅過程做出精確的表達，能建 立科學的森林面積資源的仿真理論。它無論是在理論上，還是在實際應用中都有重要的 意義。 2 森林面積資源仿真的數學模型 2.1仿真的連續(xù)型數學模型 描述森林面積資源的方程一般有： 2.1.1 同齡純林的林齡面積分布結構變化方程[1]： P([pic],t)： 表示同齡純林面積的林齡分布密度函數 μ([pic],t)： 表示同齡純林的相對采消率函數 [pic] : 表示單位時間內林齡的變化率，即林齡與時間的換算量[pic] P([pic],0)： 表示方程的初始條件 φ(t) ： 表示 t 時刻，單位時間內造林更新并郁閉成林的面積 β(t) : 表示t時刻單位時間內的造林更新率 γ(t) : 表示t時刻單位時間內的郁閉成林率 [pic]： 為干擾項 2.1.2 面積按樹種混交系數·平均年齡分布結構變化方程[2]： A： 表示林分的樹種組成系數或比例 [pic]： 表示平均年齡 [pic]：表示t時刻樹種混交的隸屬度為[pic],平均年齡數為[pic]的分布密度函數 連續(xù)型模型在理論分析中是比較方便的，但不適合實際計算和統(tǒng)計分析。為了與現(xiàn)實 中采集的數據相結合和實現(xiàn)森林面積資源數學模型的計算機仿真，須對模型離散化。 2.2 離散型模型 為了便于討論，下面以同齡純林的林齡面積分布結構變化方程為例，給出該模型離散 化后以向量形式表示的代數方程和以矩陣形式表示的矩陣方程。 （1）式離散化后的向量代數方程形式為[1]： [pic] （3） S(k+1): 表示k+1 分期同齡純林面積向量 H(k): 表示面積的林齡轉移矩陣 S(k): 表示k分期同齡純林面積向量 M(k): 表示消耗率矩陣 G(k) : 表示同齡純林因突然大面積的火災、濫砍亂伐等干擾因素使森林減少的面積向量S(0 ) : 表示初始條件向量 （1）式離散化后的矩陣方程為[1]： [pic] =[pic].[pic]+ γ(k)β(k) [pic][pic]-[pic] （4） M(k) ： 為消耗率矩陣 H(k) ： 為林齡面積轉移矩陣 Si(k+1) ： 為k+1 分期第i齡級的面積 Si(k) ： 為k分期第i齡級的面積 β(k) ： 為k分期的造林更新率 γ(k) ： 為k分期的郁閉成林率 [pic] ： 為k分期第i齡級的同齡純林因突然的干擾因素使森林減少的面積 離散化后上述矩陣方程可視為是同齡純林面積資源仿真的仿真模型，便可根據分布參 數的仿真語言編制仿真軟件，從而在計算機上實現(xiàn)仿真計算。 3 森林面積資源仿真的計算設計 3．1仿真的輔助計算設計 要在計算機上實現(xiàn)森林面積資源的仿真計算，需要許多參數和數據，而這些參數和數 據必須預先從許多圖面資料、調查數據、統(tǒng)計資料經過提取、統(tǒng)計、分析或外推而得到 。有了這些準確的參數和數據才能進行準確的計算機仿真。而為了獲取這些參數和數據 所進行的計算設計稱為仿真的輔助計算設計，一般有： 3.1.1參數的提取 這些參數包括H(k)、M(k)、 γ(k)、β(k)等，要提取出這些參數可以通過林業(yè)上的連續(xù)清查（一類調查）、二類調查 (經理調查)、檔案材料以及林場自行設計的記錄檔案，也可通過遙感圖象的衛(wèi)片、航片 的判讀材料獲取。 3.1.2參數的外推 要利用上述模型進行森林面積資源的預測，必須進行參數的外推。因為在預測時我們 首先必須對這些參數函數的動態(tài)特性，未來的變化趨勢給出合理的假定和較準確的估計 ，使之符合未來實際發(fā)展情況。參數的外推，通常的做法是對過去的統(tǒng)計資料或調查資 料進行統(tǒng)計分析研究，采用平行外推、曲線擬合、參數辯識、時序模型等方法得出它的 函數方程，就可進行參數外推。 3.1.3指標參數計算設計 在計算機仿真時往往需要對森林的面積資源進行評估，則需要預先知道評估時所使用 的許多指標參數，而這只要確定了各種指標參數的計算公式和方法，在仿真時就很容易 計算出各種指標參數，對森林的面積資源進行適時而客觀的評價。 例如：計算某林業(yè)局某樹種的總面積： 已知某樹種各個林齡的面積數即分布結構，則可利用下面的公式計算出它的結構總面 積： 結構總面積：S(k)= [pic] Si(k)表示k分期林齡為i的面積數 ； 計算某林業(yè)局某樹種的面積占林業(yè)局總面積的比率[pic]: [pic]=[pic] [pic]=S(k)為k時某樹種分布結構的總面積 S總(k)為k時某林業(yè)局的總面積 同樣利用上式可計算出某樹種在林業(yè)局的森林覆被率。 同樣還可以建立許多指標參數的計算公式和方法，如結構齡級的齡組面積、林齡分布 率及其特征數等，以及各種經營指標參數如成過熟林系數、育林系數、可采系數、撫育 系數、擇伐指數等一些指標參數，都可直接給出計算的公式和方法。這些指標參數為森 林面積資源管理提供各種定量的評價指標。 3.2 仿真計算的設計 通過仿真計算設計，確定仿真的目標和要達到的目的。計算機仿真具有模擬、試驗、 協(xié)調、控制、優(yōu)化、決策等功能和作用。森林面積資源的計算機仿真計算其主要的內容 有： 3.2.1 模擬森林面積資源的動態(tài)變化和生滅過程的仿真計算 造林更新率β(k)反映了森林資源的“生”，消耗率矩陣M(k)反映了森林資源的“死”。從 方程（3）和（4）式都不難看出，當已知k分期的森林面積資源S(k)，將它和提取或預估 出的參數β(k)、M(k)、γ(k)的值代入方程，就可仿真計算出k+1分期的森林面積資源S(k +1)。對上述參數在不同的時間給出不同的狀態(tài)值，則可模擬任意時間的森林面積資源。 因此，該仿真模型能模擬任意時間森林面積資源的動態(tài)變化和生滅過程。利用這種計算 可進行森林面積資源的監(jiān)測和預測、森林面積資源的再現(xiàn)和現(xiàn)狀分析等。森林面積資源 的仿真計算為森林面積資源的科學管理、系統(tǒng)管理、信息管理、動態(tài)管理提供了理論依 據和技術方法。 3.2.2 各種經營措施、方案的優(yōu)化仿真計算 若已知k分期和k+1分期的森林面積資源S(k)和S(k+1)，就可以仿真計算出k分期的造 林更新率β(k) 、郁閉成林率γ(k) 、消耗率M(k)。而這些不同的參數就代表了不同的經營措施，因為不同的參數β(k)、M( k)反映了不同的造林更新措施和不同的采伐方式或措施。同樣，它也反映了森林面積資 源的不同經營方案，例如改變任一林齡的采伐面積，引起參數M(k)改變，而參數的改變 必然引起森林面積資源的狀態(tài)改變。從而可以對各種參數不同的取值進行模擬、試驗和 優(yōu)化，實現(xiàn)森林面積資源各種經營措施、方案的優(yōu)化仿真計算。 3.2.3 結構調整最優(yōu)控制的仿真計算 從方程（3）和（4）式都可以看出，若已知k和k+1時的狀態(tài)值S(k)和S(k+1)，k時的 狀態(tài)值為起始值，k+1時的狀態(tài)值為要求達到的值或理想的狀態(tài)值，如何使目前的狀態(tài)值 平穩(wěn)的過渡到理想的狀態(tài)值。我們可以通過調控參數μ(死)、β（生）、γ（保存率）的值 來實現(xiàn)，這就是森林面積資源的控制計算。實際生產中對同齡純林的面積資源變化控制 就是通過采伐（μ）和造林更新（β或γ） 的調整而實現(xiàn)的。同樣還可進行森林面積資源的結構調整的計算機仿真計算；森林面積 資源可持續(xù)發(fā)展定量描述的模擬、試驗、控制、優(yōu)化、決策的計算機仿真計算；森林面 積資源與生態(tài)環(huán)境和人類社會關系協(xié)調的試驗、模擬、控制、優(yōu)化、決策的計算機仿真 計算；還可以用森林面積資源的計算機仿真控制計算代替?zhèn)鹘y(tǒng)森林資源經營方案的編制 。 4 森林面積資源的計算機仿真在森林資源系統(tǒng)管理中的應用 通過森林面積資源的計算機仿真，可以模擬森林面積資源在任意時間內的動態(tài)變化規(guī) 律和生滅過程；可以在很短時間里或在隨意的時間里就能觀察到森林面積資源動態(tài)變化 和生滅過程的全過程；以及估計某些參數如造林更新率、采伐率和枯死率等對森林資源 變化的影響，為定量地了解科學地分析森林面積資源的動態(tài)變化規(guī)律及其生滅過程的特 點和森林環(huán)境產生的影響等提供科學而有效的技術方法，也為森林面積資源的戰(zhàn)略發(fā)展 、林業(yè)方針政策、各種經營措施、林業(yè)規(guī)劃方案等等的制定、決策等提供大量的參數和 決策依據。 通過森林面積資源的計算機仿真，可以對各種經營措施、規(guī)劃、設計、決策等方案和 策略進行各種實驗、比較、優(yōu)化，為其決策找出最佳方案，最好措施，最優(yōu)策略。可以 協(xié)調和調整森林面積資源的內部結構組成要素與要素、要素與整體（系統(tǒng)）的關系，維 持森林面積資源系統(tǒng)的動態(tài)平衡和穩(wěn)定，并促使系統(tǒng)走向組織化、有序化；可以協(xié)調、 調控森林面積資源與環(huán)境、森林面積資源與人（社會）的關系，從而使森林生態(tài)系統(tǒng)處 于協(xié)調、穩(wěn)定、平衡、發(fā)展的狀態(tài)。 通過森林面積資源的計算機仿真，為森林面積資源的各種調控提供依據和方法。森林 面積資源系統(tǒng)是一個影響因素比較多組成要素比較復雜的巨大控制系統(tǒng)，很難用物理模 型來進行仿真，如果用描述森林面積資源系統(tǒng)動態(tài)變化和生滅過程的數學模型進行計算 機數字仿真，則森林面積資源系統(tǒng)的復合控制、最優(yōu)控制、自適應控制、過程控制等等 新的調控方法將成為可能，從而為森林面積資源的預測和監(jiān)測、動態(tài)分析、結構調整、 指標參數體系等...
0105064森林面積資源的計算機仿真及其應用(doc)