王志明，張 亮，迮興業(yè)，趙思萌，王選遙，史麗雯，張 健

(1.長春科技學院，吉林 長春 130600；2.吉林省乾安縣林業(yè)局第三機械林場，吉林 松原 131400；3.吉林農業(yè)大學，吉林 長春 130118；4.長春師范大學，吉林 長春 130031；5.長春醫(yī)學高等?？茖W校，吉林 長春 130013；6.吉林省林業(yè)科學研究院，吉林 長春 130031)

花布燈蛾[Camptolomainteriorata(Walker)]是蒙古櫟(Quercusspp.)的重要害蟲，在中國南方和北方都曾爆發(fā)成災[1-13]。在中國東北區(qū)域，該蟲主要危害蒙古櫟(QuercusmongolicaFisch. ex Ledeb.)芽苞和葉片。在吉林省發(fā)生年份中，林分內蒙古櫟被害率可達90 %以上。盡管在花布燈蛾發(fā)生時，各級主管部門予以高度重視，調動一切可能資源給予撲滅，但所有措施均為災情過程的處理，實質為治蟲于臨時。近30年來經驗表明，花布燈蛾屬食葉性突發(fā)物種，對其缺少預警措施，是目前防控該蟲的最大短板。

害蟲的發(fā)生除與其內稟性相關外，尚與林分因子關系密切。林分因子對于林內害蟲既有提供生存條件，又有制約增殖的作用。林分因子包括樹種組成、林齡、胸徑、樹高、郁閉度等因子?？傮w來看，害蟲的爆發(fā)，其所居的林分因子應是發(fā)生的基礎因素以及成災的主導因子之一。

截至目前，有關花布燈蛾預報方面的研究極少報道，更無有關其發(fā)生危害與林分因子關系方面的研究報道。筆者從林分因子與花布燈蛾發(fā)生危害的關系為切入點，試圖分析出利于花布燈蛾的發(fā)生主要因子，為今后深入研究防治花布燈蛾技術提供支撐。

在本文內容研究時，筆者曾用傳統(tǒng)的生物統(tǒng)計分析方法分析數據，但囿于數據難于滿足生物統(tǒng)計的要求，未獲成功。后改用粗糙集理論分析，取得較好結果。

粗糙集理論是1982年由波蘭數學家Pawlak教授提出的一種研究不完整、不確定知識和數據的表達、學習、歸納的方法。其主要思想是在保持分類能力不變的前提下，通過知識簡約而得到問題的分類規(guī)則和決策規(guī)則[14]。由于其具有只利用數據本身而不需要任何先驗知識的特點，故在很多領域得到成功應用[15-17]。

1 研究方法

1.1數據來源

依據花布燈蛾在吉林省最近一次規(guī)模發(fā)生情況(2011年)檔案記載，在吉林省東豐縣、輝南縣、梅河口市、磐石市等4縣(市)區(qū)域內，以小班為單位，收集花布燈蛾發(fā)生程度和相應林分各項因子，包括樹種構成、郁閉度、平均胸徑、平均樹高、林齡和林分面積。

1.2 危害指標與指標評分

經對上述4縣(市)161個小班樣本收集，以《吉林省森林病蟲鼠調查技術規(guī)程》規(guī)定的花布燈蛾危害程度，按幼蟲蟲口密度(蟲苞/株)劃分：

0 <輕度發(fā)生(低)< 0.08

0.09 <中度發(fā)生(中)< 0.14

0.15 <重度發(fā)生(高)

最終將各林分因子指標含量及受害程度劃分為“高”“中”“低”三檔。各林分因子指標等級劃分、評分規(guī)則見表1。

表1 林分因子指標評分規(guī)則

2 結果與分析

2.1決策表構建

根據現場調查，以林分因子指標為條件屬性，以受害程度等級為決策屬性構建決策表，其中論域U={1,2,3……,120}，條件屬性集C={C1C2C3C4C5C6}，決策屬性集D=ik0wgia，其中C1為蒙古櫟所占組分，C2為林分郁閉度，C3為林分平均胸徑，C4為林分平均樹高，C5為林分平均林齡，C6林分面積，d林分受害程度等級。離散化處理林因分子指標數據，再將屬性值進行約簡，同時刪去屬性值相同的數據，得到簡化的決策表，見表2。

表2 屬性值約簡后林分因子評價決策

序號條件屬性C1C2C3C4C5C6決策屬性d37101037101033810107710733910103773240101033372411073333342107333324310771010101441073331014510333310146103777314710710107714810310101073491037771015073333103517333373523333310253733333254737773255103333725633737315773333715873337315910333371601033333161103101073362373333363773337364107333736510777373661033373267107337326810771071026910771077270103377327137333101723733331

序號條件屬性C1C2C3C4C5C6決策屬性d7377710771747733310175773773176777107317777333317810777371791077733180107777718110771073182103377718331031010738471077732851033373186310373718771037310188101033331891033377190101033332917103333192337107739333101010101943371010101953771010739637710103397377107103983771077399377101010310037771010310177710101031023771010721033777101011043777107110531071010731063107101010310731010101073108310101010103

序號條件屬性C1C2C3C4C5C6決策屬性d1093107710103110310771073111710710101031123101010101021133107710101114310771071115710771010111671010101010311771010101073118710101010102119101010101010212031010101071

2.2 等價類劃分

定義1：設關系R為論域U上的等價關系，則稱U/IND(R)為U上的劃分，其中，

IND(R)={(x,y)∈U×U|?a∈C∪D,f(x,a)=f(y,a)}

對于表2中的數據論域分別按照和進行等價類劃分得：

|U/IND(C)|=104

|U/IND(D)|=3

依次減去C中的條件屬性Ci，i=1,2,3,4,5,6,再將表2中的數據論域進行等價類劃分得：

|U/IND(C-C1)|=69

|U/IND(C-C2)|=70

|U/IND(C-C3)|=95

|U/IND(C-C4)|=86

|U/IND(C-C5)|=87

|U/IND(C-C6)|=62

2.3 權重確定

定義2：令C′為條件屬性集C的一個子集，則決策屬性集D的確定域，PosC′(D)={x∈U,[x]C′?D}是指論域U中那些在現有的條件屬性集C′下能夠確定地歸入D的元素的集合。

定義3：令C′為條件屬性集C的一個子集，則決策屬性集D對C′的依賴度γC′(D)定義如下：

定義4：令C′為條件屬性集C的一個子集，則C′關于決策屬性集D的重要度SigD(C′)可定義如下：

SigD(C′)=γC(D)-γC′(D)

C′的權重?(C′)定義為：

按照給定定義2、定義3和定義4可得各條件屬性下的決策屬性的確定域：

|PosC(D)|=88

|PosC-C1(D)|=53

|PosC-C2(D)|=60

|PosC-C3(D)|=76

|PosC-C4(D)|=68

|PosC-C5(D)|=69

|PosC-C6(D)|=49

各條件屬性關于決策屬性的依賴度：

各條件屬性關于決策屬性的重要度：

SigD(C1) =γC(D)-γC-C1(D)

SigD(C2)=γC(D)-γC-C2(D)

SigD(C3)=γC(D)-γC-C3(D)

SigD(C4)=γC(D)-γC-C4(D)

SigD(C5)=γC(D)-γC-C5(D)

SigD(C6)=γC(D)-γC-C6(D)

各條件屬性關于決策屬性的權重：

?(C1)=0.228

?(C2)=0.183

?(C3)=0.078

?(C4)=0.131

?(C5)=0.124

?(C6)=0.255

各評價指標權重及排序結果見表3。

表3 各林分因子的權重及排序

從表3中影響花布燈蛾爆發(fā)的林分因子評價指標權重排序上看，林分面積因子和蒙古櫟百分比構成因子對花布燈蛾發(fā)生具有重要影響，其次為林分郁閉度因子、平均樹高因子和林齡因子對于花布燈蛾發(fā)生的影響大致相同。

3 結論與討論

對花布燈蛾發(fā)生危害程度與相關林分因子的關系進行分析，需要考慮多個不直接相關指標屬性，但為達到這個目標，運用傳統(tǒng)、常用的生物統(tǒng)計方法很難進行。為此，筆者幾經探討，最終選定利用粗糙集理論和方法在消除冗余信息和處理不確定性信息等方面的優(yōu)勢，將林分因子指標權重確定轉化為粗糙集中的屬性重要性評價，從而確定花布燈蛾危害的林分因子評價指標權，有效地減少其計算主觀性，做到模糊性和精準性的統(tǒng)一。

應用粗糙集理論判定林分因子與花布燈蛾發(fā)生相關性，首先遵循原則是盡可能提高精準及確定因子是否對提供的信息有用。在諸多林分因子中，對花布燈蛾危害成災起重要作用的當屬林分面積因子和林分所占蒙古櫟百組分因子。在長白林區(qū)，蒙古櫟多以大面積次生闊葉林存在。其林分內，蒙古櫟所占比例通常高達50 %以上。事實表明，花布燈蛾爆發(fā)危害均是大面積突發(fā)，沒有在1個或幾個小班的小面積成災。究其原因，林分環(huán)境對花布燈蛾影響有緊密關系。大面積林分可以保持森林小氣候遠優(yōu)于小面積林分。森林小氣候特點為光照低、風速小、濕度大、最高氣溫低和最低氣溫高，林中空間和林地溫度日變化和年變化都比較小，對林內溫度、濕度、光照、水分、風和空氣的穩(wěn)定和優(yōu)化具有促進作用，特別是連片林分，利于花布燈蛾發(fā)育和繁殖。在較大面積林分內，既穩(wěn)定和強化了環(huán)境因子對花布燈蛾的物理生存條件，也為花布燈蛾種群基數積累提供幫助，更為花布燈蛾居群的林間基因交換和信息交流提供保證。

其次，與花布燈蛾成災有關的林分因子為林分郁閉度。根據FAO(聯合國糧農組織)規(guī)定，郁閉度≥0.70的林分為密林，郁閉度0.20～0.69的林分為中度郁閉林，郁閉度<0.20以下為疏林。林分在不同郁閉度下，由于小氣候等物理因素影響，林內生物結構將發(fā)生相應變化?；ú紵舳暝谥袊|南以北區(qū)域，以3齡幼蟲在樹干基部的枯枝落葉中的蟲苞內越冬，蛹期也是在樹根基附近的枯枝落葉層中[3,10,11]。幼蟲期和蛹期的隱蔽度和生存條件極為重要。調查結果已知，花布燈蛾發(fā)生的林分大多在中度郁閉林分以上，且集中在密林。這些因素的單個或綜合作用直接影響到花布燈蛾對不同郁閉度條件下微環(huán)境適應性，從而產生對不同郁閉度林生存規(guī)律的差異。調查表明，在適當范圍內，花布燈蛾種群數量會隨著林分郁閉度增加而上升，這或許與郁閉度帶來的對于林分條件改變有關。更深層次的研究有待今后進一步進行。

上述3項屬性指標占整個6項屬性指標2/3，客觀上基本反映出林分主要因子對花布燈蛾影響程度。其余平均樹高因子和林齡因子對花布燈蛾爆發(fā)的影響大致相同，而平均胸徑因子對花布燈蛾爆發(fā)影響最弱。

從結果上來看，應用粗糙集理論判定林分因子與花布燈蛾發(fā)生相關性符合實際情況，其研究結果豐富了傳統(tǒng)的分析害蟲相關因子方法，且計算簡單、方便，結果合理。

[1]陳天璘.花布燈蛾(CamptolomainteriorataWalker)生活習性及防治的研究[J].林業(yè)科學,1963，8(1)：74-77.

[2]席惠蘭.花布燈蛾的初步觀察[J].沈陽農業(yè)大學學報，1986，17(4)：81-83.

[3]侯彬.花布燈蛾觀察初報[J].吉林林業(yè)科技, 1988,(5)：30-32.

[4]徐璞，韓英材，徐少風.花布燈蛾防治對策[J].吉林林業(yè)科技, 1989,(4)：30.

[5]徐和光,李東軍,秦緒兵,等.花布燈蛾生物學特性與防治技術報告[J].山東林業(yè)科技，1998，(5)：32-34.

[6]席忠誠. 花布燈蛾生物學特性及綜合防治技術[J].甘肅林業(yè)科技，1999，(2)：35-37.

[7]殷吉龍，史金富，戈文達，等.花布燈蛾的生物學特性及人工防治技術研究[J].經濟技術協(xié)作信息,2006，(13)：63.

[8]鄭文云，高德三.柞樹害蟲花布燈蛾生物學特性的研究[J].林業(yè)科技，2001，26(4)：22-25.

[9]王衛(wèi)平，李貴波，宋書利，等.寬甸縣柞蠶場花布燈蛾危害及防治[J].遼寧絲綢，2009，(2)：38-39.

[10]婁杰，張鐵利，鄭柏華，等. 花布燈蛾生物學特性及防治[J].遼寧林業(yè)科技，2010，(4)：46-47.

[11]房立新，陳連正，張聯合，等. 花布燈蛾生物學特性及綜合防治技術[J]. 吉林林業(yè)科技， 2012，41(4)：53-54.

[12]趙瑞興,邢禮國,李素梅. 遼寧省花布燈蛾發(fā)生的原因分析和防治對策[J].遼寧林業(yè)科技2011，(6)：32-33.

[13]賈中鵬，張曉軍，張健.花布燈蛾生物學特性及綜合防治技術[J].吉林農業(yè)，2013，(3)：64-64.

[14]Pawlak Z.RoughSets:TheoreticalAspectsofReasoningaboutData[M].Dordrecht：Kluwer Academic Publishsers,1991.

[15]張文修，吳偉志，梁吉業(yè)，等.粗糙集理論與方法[M].北京：科學出版社，2001.

[16]王禹萌，鄒春紅，呂晗,等. 基于粗糙集理論的長白山火山泥石流影響因子權重分析[J].防災減災學報，2013,29(2)：28-32.

[17]葉回春，張世文，黃元仿，等.粗糙集理論在土壤肥力評價指標權重確定中的應用[J].中國農業(yè)科學,2014，47(4)：710-717.