李思雯
(黑龍江省林業(yè)監(jiān)測(cè)規(guī)劃院,黑龍江哈爾濱150080)
風(fēng)倒是森林中的常見(jiàn)現(xiàn)象,研究其機(jī)理具有理論和應(yīng)用兩方面的意義。在理論研究上,風(fēng)是一個(gè)生態(tài)因子,風(fēng)倒是林隙的主要成因,林隙對(duì)森林動(dòng)態(tài)有重要的影響[1-4],研究風(fēng)倒機(jī)理有助于揭示林隙動(dòng)態(tài)規(guī)律及其生態(tài)學(xué)影響[5]。在實(shí)際應(yīng)用中,風(fēng)倒是森林中多發(fā)的災(zāi)害,會(huì)造成木材產(chǎn)量和質(zhì)量的下降。英、美等國(guó)每次災(zāi)害性風(fēng)倒造成的木材損失都達(dá)數(shù)百萬(wàn)立方米,每年非災(zāi)害性風(fēng)倒造成的木材減產(chǎn)均在15%以上[6]。目前我國(guó)全國(guó)風(fēng)倒造成的損失也很?chē)?yán)重[7-9],僅1996年“9615”號(hào)臺(tái)風(fēng)給湛江市林業(yè)帶來(lái)的直接損失就達(dá)21.1億元[10]。只有弄清風(fēng)倒機(jī)理,建立相關(guān)模型,建立科學(xué)的風(fēng)倒風(fēng)險(xiǎn)評(píng)估體系[11],才能在森林經(jīng)營(yíng)中采取科學(xué)的措施來(lái)減少風(fēng)倒帶來(lái)的損失。該文對(duì)風(fēng)倒機(jī)理的研究現(xiàn)狀進(jìn)行綜述,并對(duì)未來(lái)的研究提出建議。
風(fēng)倒包括干折(stem breakage)和根倒(uprooting)2種,主要受林木生物力學(xué)特性、土壤支持力(root anchorage)和風(fēng)的擾動(dòng)3個(gè)因素影響[12]。其中,林木的生物力學(xué)特性與立木年齡、樹(shù)干強(qiáng)度、干型、樹(shù)冠形狀和大小及樹(shù)木健康狀況等因素有關(guān);土壤支持力受土壤剪切強(qiáng)度、含水率、樹(shù)木根系深度、寬度和重量等影響,二者決定了樹(shù)木的總抗風(fēng)性。風(fēng)的擾動(dòng)主要受風(fēng)速、風(fēng)壓、頻率、不同林分密度和復(fù)雜地形下的風(fēng)場(chǎng)分布影響。因此,風(fēng)倒機(jī)理的研究涉及力學(xué)中的固體力學(xué)、流體力學(xué)、土力學(xué)和林學(xué)中的森林生態(tài)學(xué)、森林土壤學(xué)、森林培育學(xué)、森林氣象學(xué)、森林病理學(xué)等多個(gè)學(xué)科,需要這些學(xué)科的交叉與合作。
風(fēng)倒機(jī)理的研究在國(guó)外開(kāi)展較早。英國(guó)、加拿大、美國(guó)等國(guó)家對(duì)樟子松(Pinus sylvestris)、西加云杉(Picea sitchensis)、挪威云杉(Picea abies)、西鐵杉(Tsuga heterophylla)等樹(shù)種的風(fēng)倒機(jī)理進(jìn)行了深入、廣泛的研究。通過(guò)動(dòng)力學(xué)理論推導(dǎo),結(jié)合野外試驗(yàn)中林木和土壤力學(xué)參數(shù)的測(cè)定、地形數(shù)字高程模型和實(shí)物模型風(fēng)洞試驗(yàn),獲得了復(fù)雜地形下和不同林分密度時(shí)的風(fēng)場(chǎng)分布,建立了一系列力學(xué)方程來(lái)描述林木的抗風(fēng)性[13-15],對(duì)風(fēng)倒的力學(xué)機(jī)理進(jìn)行了闡述,形成若干林木風(fēng)倒機(jī)理模型,如 MC2[16]、HWIND[17]、FORESTGALES[18]和WINDARC[19]等,它們代表了國(guó)際風(fēng)倒機(jī)理研究的先進(jìn)水平。其中,典型的是Peltola[13]將樟子松簡(jiǎn)化成一個(gè)地上的簡(jiǎn)支梁,分析了樹(shù)木在持續(xù)穩(wěn)定風(fēng)載作用下的穩(wěn)定性,然后計(jì)算能夠?qū)?shù)吹倒的風(fēng)速。Gardiner[18]研究了針葉樹(shù)的樹(shù)干模型,把它等效成一個(gè)逐漸變細(xì)的懸臂梁,在樹(shù)高的70%處加一個(gè)質(zhì)量體,對(duì)系統(tǒng)進(jìn)行了動(dòng)力學(xué)分析,獲得系統(tǒng)的自然頻率、平均和波動(dòng)位移。
我國(guó)目前還沒(méi)有系統(tǒng)的研究風(fēng)倒學(xué)機(jī)理。在國(guó)內(nèi)涉及風(fēng)倒的研究中,李國(guó)旗等[20]根據(jù)現(xiàn)有公式對(duì)不同風(fēng)壓和高度處樹(shù)木應(yīng)力的分布進(jìn)行了研究。陳少雄等[21-23]采用野外調(diào)查和統(tǒng)計(jì)方法分別研究了泡桐、桉樹(shù)、木麻黃等樹(shù)種的抗風(fēng)性與個(gè)體特征、林分特征之間的關(guān)系。其他基本是從生態(tài)學(xué)角度出發(fā)研究林隙的[24-26]。上述研究涉及風(fēng)倒機(jī)理較少,在2000年以前國(guó)內(nèi)關(guān)于風(fēng)倒機(jī)理的專項(xiàng)研究基本屬于空白。之后,開(kāi)展了較多的研究,主要是利用國(guó)外模型,對(duì)國(guó)內(nèi)一些樹(shù)種,主要是對(duì)東北的云杉進(jìn)行了研究[27-31]。其中宋曉鶴等[28-30]以云杉為研究對(duì)象,從不同角度研究了樹(shù)木風(fēng)倒的力學(xué)機(jī)理。他主要研究了云杉的風(fēng)倒靜力學(xué)特性,測(cè)定了相應(yīng)的力學(xué)性能參數(shù),建立了該樹(shù)種的風(fēng)倒靜力學(xué)模型。王琳[29]將云杉樹(shù)干看作橫截面隨高度呈冪函數(shù)的彈性桿,一端固定,一端自由并連有質(zhì)量團(tuán),分析了該彈性桿樹(shù)干的振動(dòng)特性。賴秋明[30]基于經(jīng)典梁振動(dòng)理論和Ansys有限元法對(duì)云杉風(fēng)倒動(dòng)力學(xué)問(wèn)題進(jìn)行研究,得到了云杉模型的固有頻率和各階模態(tài),分析了影響云杉模型固有頻率的幾個(gè)因子,建立了該樹(shù)種的風(fēng)倒動(dòng)力學(xué)模型。賈杰等[31]在此基礎(chǔ)上,將林木進(jìn)一步簡(jiǎn)化成一根無(wú)質(zhì)量彈性桿連結(jié),給出林木的動(dòng)態(tài)微分方程以及模型邊界條件,考慮了樹(shù)干底部彎矩和風(fēng)力的關(guān)系,建立了無(wú)阻尼和有阻尼情況下林木風(fēng)倒的簡(jiǎn)化動(dòng)態(tài)模型。此外,賈杰等[27]利用Von-Karman風(fēng)譜模型,在頻率域內(nèi)分析了懸鈴木風(fēng)倒的機(jī)理,通過(guò)對(duì)力譜、基礎(chǔ)彎矩譜的理論推導(dǎo),給出陣風(fēng)頻率的能量分配問(wèn)題,從而獲得相關(guān)的參數(shù)變量,尤其是振動(dòng)彎矩極值,并指出當(dāng)這個(gè)極值超過(guò)莖和根系所能承受的彎矩時(shí),風(fēng)倒便會(huì)發(fā)生。這些研究對(duì)于森林防風(fēng)具有重要意義。
現(xiàn)有的模型正確描述了風(fēng)倒機(jī)理的許多方面,但受研究條件等限制在系統(tǒng)設(shè)計(jì)時(shí)對(duì)關(guān)鍵環(huán)節(jié)進(jìn)行了簡(jiǎn)化,造成模型與實(shí)際情況之間的較大誤差[14]。其中最主要的問(wèn)題是:活立木具有生長(zhǎng)適應(yīng)能力,其生物力學(xué)特性決定了它對(duì)風(fēng)的響應(yīng)。這種生物力學(xué)響應(yīng)既有彈性形變,又不完全等同。在低于臨界風(fēng)速(critical windspeed)的持續(xù)風(fēng)的作用下樹(shù)干(乃至樹(shù)冠)會(huì)發(fā)生形變積累,而不是如現(xiàn)有模型中所描述的重新回到直立平衡位置。活立木在風(fēng)倒時(shí)不僅發(fā)生脆性斷裂,還會(huì)出現(xiàn)塑性形變,而且這2種過(guò)程的比例和作用隨立木年齡變化而變化。以前的一些模型把樹(shù)干作為非生物彈性體、將風(fēng)倒過(guò)程僅按脆性斷裂處理而未考慮活立木的生物力學(xué)特性。因此,它們只能預(yù)測(cè)超過(guò)固定臨界風(fēng)速時(shí)發(fā)生的風(fēng)倒,而不能揭示活立木生物特性對(duì)風(fēng)倒過(guò)程的影響,也無(wú)法解釋樹(shù)木在較低風(fēng)速時(shí)發(fā)生的風(fēng)倒和形變積累對(duì)臨界風(fēng)速的影響。這是現(xiàn)有模型與實(shí)際情況出入最大的地方。國(guó)內(nèi)研究[28-30]對(duì)這方面進(jìn)行了改進(jìn),將樹(shù)干看作是彈性桿,但其模型過(guò)于簡(jiǎn)化,沒(méi)有考慮樹(shù)冠的形狀對(duì)風(fēng)倒的影響。
其他存在的問(wèn)題有:①模型將樹(shù)冠看成具有連續(xù)介質(zhì)的傘蓋式形體,把復(fù)雜的根系簡(jiǎn)化為規(guī)則幾何體,未考慮其孔性。②模型在計(jì)算土壤對(duì)林木地上部分的支持時(shí),僅通過(guò)林木地上與地下的質(zhì)量比來(lái)估測(cè),并未充分考慮土壤的力學(xué)特征。③模型雖然考慮了林木個(gè)體的固有頻率、阻尼系數(shù)等動(dòng)力學(xué)量,但未考慮林木間共振對(duì)風(fēng)倒在林分中的傳播效應(yīng)。④模型沒(méi)有考慮動(dòng)力學(xué)以及靜力學(xué)量在風(fēng)(流場(chǎng))、林木和土壤(固場(chǎng))等多物理場(chǎng)中的耦合效應(yīng)。此外,由于風(fēng)倒的突發(fā)性和危險(xiǎn)性,其過(guò)程難以記錄。許多研究是在風(fēng)倒后才進(jìn)行,往往采集不到所需的數(shù)據(jù),這極大限制了對(duì)風(fēng)倒機(jī)理的研究并制約了風(fēng)倒機(jī)理模型在生態(tài)學(xué)研究和森林經(jīng)營(yíng)中的應(yīng)用。
在方法上,建議進(jìn)一步開(kāi)展多學(xué)科的交叉,將動(dòng)力學(xué)原理和林木生物學(xué)、生態(tài)學(xué)特性結(jié)合,采用非線性有限元法,完善樹(shù)冠、樹(shù)根的多孔彈性介質(zhì)三維模型。合理簡(jiǎn)化樹(shù)干模型,可以把樹(shù)干作為生物彈性體處理,模擬自然風(fēng)對(duì)樹(shù)干造成的彎曲變形累積、確定不同彎曲下的臨界風(fēng)速和年齡對(duì)風(fēng)倒中塑變和脆裂過(guò)程的影響。還應(yīng)該加強(qiáng)與試驗(yàn)資料的結(jié)合,利用土力學(xué)計(jì)算根系中應(yīng)力與阻力彎矩間的定量關(guān)系,通過(guò)野外現(xiàn)場(chǎng)力學(xué)試驗(yàn)和風(fēng)洞試驗(yàn)對(duì)復(fù)雜地形下的風(fēng)場(chǎng)分布和力學(xué)方程參數(shù)化,最終揭示多物理場(chǎng)耦合的風(fēng)倒的非線性生物動(dòng)力學(xué)機(jī)理。這種思路是一種生物動(dòng)力學(xué)機(jī)理思想,可以在現(xiàn)有單純動(dòng)力學(xué)基礎(chǔ)上更好地揭示風(fēng)倒的機(jī)理。
在研究對(duì)象上,建議今后應(yīng)研究不同樹(shù)種的風(fēng)倒機(jī)理。目前國(guó)際上對(duì)針葉樹(shù)的風(fēng)倒研究較多,這與一些地區(qū)針葉樹(shù)風(fēng)倒嚴(yán)重有關(guān)。但在一些地方,闊葉樹(shù)風(fēng)倒也時(shí)有發(fā)生,今后應(yīng)加強(qiáng)闊葉樹(shù)風(fēng)倒機(jī)理的研究。
[1]于振良,趙士洞.林隙(Gap)模型研究進(jìn)展[J].生態(tài)學(xué)雜志,1996,16(2):42-46.
[2]臧潤(rùn)國(guó),徐化成.林隙干擾研究進(jìn)展[J].林業(yè)科學(xué),1998,34(1):90-98.
[3]臧潤(rùn)國(guó).林隙更新動(dòng)態(tài)研究進(jìn)展[J].生態(tài)學(xué)雜志,1998,17(2):50-58.
[4]SCHAETZL R J,BURN SF,JOHNSON D L,et al.Tree uprooting:Review of impacts on forest ecology[J].Vegetatio,1989,74:132-145.
[5]ULANOVA N G.The effects of windthrow on forests at different spatial scales:A review[J].Forest ccology and management,2000,135:155-167.
[6]MOORE J,QUINE CP.A comparison of the relative risk of wind damage to planted forests in Border Forest Park,Great Britain,and Central North Island,New Zealand[J].Forest ecology and management,2000,135:345-353.
[7]陳少雄,楊民勝,王理平.尾葉桉造林密度與蓄積量、抗風(fēng)和材性關(guān)系研究[J].林業(yè)科學(xué)研究,1998,11(4):435-438.
[8]王玉岱,朱愛(ài)萍,陳建華,等.臺(tái)風(fēng)對(duì)泰安市綠化樹(shù)木危害情況調(diào)查[J].山東林業(yè)科技,1998(3):45-47.
[9]劉震,劉巖.大興安嶺林區(qū)伐區(qū)風(fēng)倒風(fēng)折原因探討[J].林業(yè)勘察設(shè)計(jì),2007(3):38-40.
[10]陳綬柱,岑奮,吳澤鵬.沿海沙岸防風(fēng)固沙木麻黃試驗(yàn)示范林抗御臺(tái)風(fēng)分析[J].廣東林業(yè)科技,1999,15(1):26-30.
[11]GARDINER BA,QUINE CP.Management of forests to reduce the risk of abiotic damage-a review with particular reference to the effects of strong winds[J].Forest Ecology and management,2000,135:261-277.
[12]SCHAETZL R J,JOHNSON D L,BURNSSF,et al.Tree uprooting:Review of terminology,process,and environmental implications[J].Canadian journal of forest research,1989,19:1-11.
[13] PELTOLA H,KELLOMAKI S.A mechanistic model fore calculating windthrow and stem breakage at stand edge[J].Silva Fenn,1993,27(2):99-111.
[14]RUEL J C,PIN D,SPACE L,et al.The estimation of wind exposure for windthrow hazard rating:Comparison between Strongblow,MC2,Topex and a wind tunnel study[J].Forestry,1997,70(3):254-266.
[15]ENGLANDA H,BAKERDJ,SAUNDERSONST.A dynamic analysis of windthrow of trees[J].Forestry,2000,73(3):225-237.
[16]BENOIT R,BESGAGNé M,PELLERIN P,et al.The Canadian MC2:A semi-Lagrangian,semi-implicit wide-band atmospheric model suited for fine-scale process studies and simulation[J].Mon Wea Rev,1997,125:2382-2415.
[17]PELTOLA H,KELLOMAKIH,IKONENV P.A mechanistic model for assessing the risk of wind and snow damage to single trees and stands of Scots pine,Norway spruce,and birch[J].Canadian journal of forest research,1999,29:647-661.
[18]GARDINER B A,PELTOLA H,KELLOMAKI S.The development and testing of models for predicting the critical wind speed to damage trees[J].Ecological modeling,1999,119:110-123.
[19]LEKESV,DANDUL I.Using airflow modeling and spatial analysis for defining winddamage risk classification(WINDARC)[J].Forest ecology and management,2000,135:331-344.
[20]李國(guó)旗,安樹(shù)青,張紀(jì)林,等.海岸帶防護(hù)林4種樹(shù)木的風(fēng)壓應(yīng)力分析[J].南京林業(yè)大學(xué)學(xué)報(bào),1999,23(4):76-80.
[21]陳少雄,王觀明,羅建中.桉樹(shù)幼苗不同株行距配置抗臺(tái)風(fēng)效果的研究[J].林業(yè)科學(xué)研究,1995,8(5):582-585.
[22]施士爭(zhēng),倪善慶,周友仁,等.泡桐樹(shù)冠結(jié)構(gòu)與抗風(fēng)性研究[J].江蘇林業(yè)科技,1996,23(2):6-9,18.
[23]陳士銀.木麻黃不同株行距配置抗臺(tái)風(fēng)效果[J].廣東林業(yè)科技,1997,13(2):44-46.
[24]國(guó)慶喜,葛劍平,馬承慧,等.長(zhǎng)白山紅松混交林林隙狀況與更新研究[J].東北林業(yè)大學(xué)學(xué)報(bào),1998,26(1):4-7.
[25]臧潤(rùn)國(guó),徐化成,高文韜.紅松闊葉林主要樹(shù)種對(duì)林隙大小及其發(fā)育階段更新反應(yīng)規(guī)律的研究[J].林業(yè)科學(xué),1999,35(3):2-9.
[26]代力民,徐振邦,楊麗韞.紅松闊葉林倒木貯量動(dòng)態(tài)的研究[J].應(yīng)用生態(tài)學(xué)報(bào),1999(10):513-517.
[27]賈杰,李洪峰,李靜輝.基于Von-Karman譜的獨(dú)立林木風(fēng)倒機(jī)理的頻率域分析[J].東北林業(yè)大學(xué)學(xué)報(bào),2013,41(8):165-169.
[28]宋曉鶴.云杉風(fēng)倒靜力學(xué)模型的建立及其分析[D].哈爾濱:哈爾濱工業(yè)大學(xué),2006.
[29]王琳.云杉風(fēng)倒動(dòng)力學(xué)模型的建立與分析[D].哈爾濱:哈爾濱工業(yè)大學(xué),2006.
[30]賴秋明.云杉風(fēng)倒動(dòng)力學(xué)問(wèn)題的研究[D].哈爾濱:哈爾濱工業(yè)大學(xué),2008.
[31]賈杰,李靜輝.林木風(fēng)倒動(dòng)態(tài)模型的建立與分析[J].森林工程,2013 29(4):63-67.