家蠶血液型膿病流行規(guī)律及其與氣象因子的相關(guān)性

唐亮 蔣滿貴 石美寧 董桂清 王霞 胡文娟 黃深惠 陳小青 黃旭華 潘志新

摘要：【目的】摸清廣西家蠶血液型膿病流行規(guī)律，并研究其與氣象因子的相關(guān)性，建立病蟲害預(yù)測(cè)預(yù)報(bào)的動(dòng)態(tài)數(shù)學(xué)模型，為精準(zhǔn)防控蠶病提供科學(xué)依據(jù)?！痉椒ā恳詮V西2008—2018年的蠶病調(diào)查數(shù)據(jù)和氣象數(shù)據(jù)為基礎(chǔ)，將家蠶血液型膿病危害情況與主要?dú)庀笠蜃舆M(jìn)行相關(guān)性分析和通徑分析，構(gòu)建回歸預(yù)測(cè)模型。【結(jié)果】2008—2018年廣西家蠶血液型膿病呈波動(dòng)發(fā)生態(tài)勢(shì)，年均發(fā)生率為11.35%，危害等級(jí)為重度（Ⅲ級(jí)）和超重度（Ⅳ級(jí)）的出現(xiàn)頻率達(dá)55.56%;月均發(fā)病趨勢(shì)呈正態(tài)分布，7月發(fā)病率最高（18.04%），其次是6和8月，發(fā)病率分別為16.00%和15.36%。夏蠶期（6—8月）廣西蠶區(qū)家蠶血液型膿病危害等級(jí)大部分達(dá)Ⅳ級(jí)，占比為62.79%，各蠶區(qū)發(fā)病率由高至低排序?yàn)橘F港（18.15%）>南寧（17.89%）>河池（16.27%）>百色（15.82%）>柳州（14.16%）;春蠶期（3—5月）發(fā)病率絕大部分在中度（Ⅱ級(jí)）以內(nèi)，發(fā)病率由高至低排序?yàn)榱荩?0.13%）>河池（9.65%）>百色（7.43%）>貴港（7.00%）>南寧（6.68%）;秋蠶期（9—10月）各蠶區(qū)發(fā)病率均在Ⅱ級(jí)以內(nèi)，差異不顯著（P>0.05）。降雨量和氣溫是影響廣西蠶區(qū)家蠶血液型膿病發(fā)病率的最主要?dú)庀笠蜃?，?jù)此建立的預(yù)測(cè)預(yù)報(bào)模型為y=0.15116+0.0037371x1+0.00018527x2+0.000051826x3-0.0000030415x4-0.00000082702x5+0.00063798x6（y為家蠶血液型膿病發(fā)病率，x1為旬平均氣溫與23.75 ℃差值的二次方，x2為上一旬降水量，x3為上兩旬平均氣溫的二次方，x4為上兩旬日照時(shí)數(shù)的二次方，x5為前三旬降雨量之和的二次方，x6為近四旬降雨量之和），發(fā)病等級(jí)的建模擬合準(zhǔn)確率為78.89%，對(duì)2018年9批監(jiān)測(cè)數(shù)據(jù)的預(yù)測(cè)準(zhǔn)確率為66.67%。【結(jié)論】廣西家蠶血液型膿病危害較嚴(yán)重，夏蠶期（6—8月）是發(fā)病高峰期。建立的預(yù)測(cè)模型可用于廣西蠶桑生產(chǎn)中家蠶血液型膿病發(fā)生危害的預(yù)測(cè)預(yù)報(bào)。

關(guān)鍵詞： 家蠶;血液型膿病;氣象因子;回歸模型;相關(guān)性分析;廣西

中圖分類號(hào)： S884.5? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)志碼： A 文章編號(hào)：2095-1191（2020）05-1217-07

Abstract：【Objective】In order to provide scientific basis for the accurate prevention and control of silkworm disease，the epidemic law of Bombyx mori nuclear polyhedrosis in Guangxi and the correlation between it and meteorological factors was studied，and a regional dynamic warning model was established. 【Method】Based on the survey data of silkworm diseases and meteorological data in Guangxi from 2008 to 2018， the correlation analysis and path analysis between the damage of B. mori nuclear polyhedrosis and the main meteorological factors were carried out，and a regression prediction model was built. 【Result】The results showed that the incidence of B. mori nuclear polyhedrosis in Guangxi was fluctua-ting from 2008 to 2018，with an average annual rate of 11.35%，and the frequency of serious（grade Ⅲ） and highly severe（grade Ⅳ） was 55.56%. The monthly incidence was normally distributed in a year， with the highest incidence rate of 18.04% in July， followed by 16.00% and 15.36% in June and August，respectively. The rate of incidence of B. mori nuclear polyhedrosis in Guangxi was mostly grade Ⅳin summer（June-August），accounting for 62.79%， and the incidence of nuclear polyhedrosis in different sericultural areas ranged from high to low was Guigang（18.15%）>Nanning（17.89%）>Hechi（16.27%）>Baise（15.82%）>Liuzhou（14.16%）. And the incidence grade of spring（March-May） was mostly in mo-derate level（grade II），and the incidence for different areas was Liuzhou（10.13%）>Hechi（9.65%）>Baise（7.43%）>Guigang（7.00%）>Nanning（6.68%）. While autumn silkworm rearing season（September-October） were below grade II and no significant difference（P>0.05）. The main factors that influenced the incidence rate of nuclear polyhedrosis of silkworm in Guangxi were rainfall and temperature. Based on this， the prediction model of the harm degree of B. mori nuclear polyhedrosis in Guangxi was y=0.15116+0.0037371x1+0.00018527x2+0.000051826x3-0.0000030415x4-0.00000082702x5+0.00063798x6（y was the incidence rate of B. mori nuclear polyhedrosis，x1 was the quadratic power of the difference between the average temperature of ten-day period and the 23.75 ℃，the x2 was the previous ten-day period precipitation，x3 was the square of the previous twenty-day period average temperature，x4 was the square of previous twenty-day period sunshine hours，x5 was the square of the rainfall of previous thirty-day period，and x6 was the sum of the rainfall in the past forty days），from which the fitting accuracy rate of the incidence grade of nuclear polyhedrosis in silkworm was 78.89% and the prediction accuracy of nine batches of monitoring data in 2018 was 66.67%. 【Conclusion】The B. mori nuclear polyhedrosis is serious in Guangxi，and its incidence peak period is in summer silkworn rearing season（June-August）. Finally，the prediction model can be used in the production of sericulture to predict the occurrence and damage of B. mori nuclear polyhedrosis in Guangxi.

Key words： Bombyx mori; nuclear polyhedrosis; meteorological factors; regression model;correlation analysis;Guangxi

Foundation item： Guangxi Key Research and Development Project（Guike AB16380102）;Construction Project of National Modern Agricultural Industry Technology System（CARS-18-SYZ14）

0 引言

【研究意義】蠶病是影響蠶繭產(chǎn)量、質(zhì)量及蠶農(nóng)經(jīng)濟(jì)效益的重要因素（祁廣軍，2013）。據(jù)不完全統(tǒng)計(jì)，廣西每年因蠶病導(dǎo)致蠶繭減產(chǎn)達(dá)15%～20%，給蠶農(nóng)造成巨大的經(jīng)濟(jì)損失，其中又以家蠶血液型膿病的危害最嚴(yán)重（閉立輝，2005）。因此，掌握家蠶血液型膿病發(fā)病規(guī)律，建立科學(xué)的防控策略和措施，對(duì)有效降低蠶病損失及促進(jìn)蠶業(yè)發(fā)展和農(nóng)民增收具有重要意義?！厩叭搜芯窟M(jìn)展】以廣西和云南等為代表的南方蠶區(qū)家蠶血液型膿病發(fā)病高峰期多出現(xiàn)在夏蠶、秋蠶期，以江蘇和安徽等為代表的長(zhǎng)江流域蠶區(qū)發(fā)病高峰期則主要出現(xiàn)在秋蠶、晚秋蠶期（白興榮等，2011;石美寧等，2015;王必華，2015;杭笑冬等，2017），即全國(guó)各大蠶區(qū)的家蠶血液型膿病發(fā)病規(guī)律因氣候條件不同而呈現(xiàn)一定差異。已有研究報(bào)道，春季冷害、氣溫升降波動(dòng)幅度大、晚秋持續(xù)低溫陰雨光照少的氣候條件，可影響大蠶齡期經(jīng)過異常，家蠶抗病力弱，導(dǎo)致細(xì)菌病、膿病等蠶病大范圍（大面積）流行或暴發(fā)（連武等，1996;韓正安和張曉紅，2014;肖建京等，2015）。牟志美等（1999）研究認(rèn)為，山東家蠶真菌病的發(fā)生與氣象因子及野外昆蟲的發(fā)生密切相關(guān)，高溫多濕是誘發(fā)家蠶多數(shù)種類真菌病的主要?dú)庀笠蛩?劉永裕等（2010）指出影響桂中地區(qū)原種蠶死蛹的主要因子是低溫冷害、高溫悶熱、高濕度和空氣干燥等;魯興萌（2012）指出溫度（5齡期或熟蠶期高溫）和桑葉（過度饑餓和桑葉偏嫩等）可直接影響家蠶對(duì)核型多角體病毒的抗性，而濕度等環(huán)境因子能從病原微生物的擴(kuò)散方面間接影響家蠶血液型膿病的流行與暴發(fā);關(guān)俊華等（2015）提出溫度是影響家蠶生長(zhǎng)發(fā)育的最重要?dú)庀笠蜃??！颈狙芯壳腥朦c(diǎn)】氣象因子對(duì)蠶病直接或間接的影響已有較多研究報(bào)道，其中家蠶微粒子病可實(shí)現(xiàn)利用氣象因子進(jìn)行預(yù)測(cè)預(yù)報(bào)（劉吉平等，2002;潘志新等，2017;胡文娟，2018），但利用氣象資料對(duì)家蠶血液型膿病進(jìn)行預(yù)測(cè)預(yù)報(bào)的研究鮮見報(bào)道。【擬解決的關(guān)鍵問題】以2008—2018年廣西典型蠶區(qū)血液型膿病發(fā)病調(diào)查數(shù)據(jù)和當(dāng)?shù)貧庀筚Y料為基礎(chǔ)，分析其相關(guān)性，構(gòu)建數(shù)學(xué)模型，為廣西防控家蠶血液型膿病提供測(cè)報(bào)依據(jù)。

1 材料與方法

1. 1 蠶病數(shù)據(jù)

廣西家蠶血液型膿病發(fā)病原始數(shù)據(jù)來源于廣西蠶桑病蟲害監(jiān)測(cè)網(wǎng)的常規(guī)蠶病監(jiān)測(cè)數(shù)據(jù)和相關(guān)科技項(xiàng)目蠶病調(diào)查數(shù)據(jù);數(shù)據(jù)匯總年份為2008—2018年，月份為每年3—10月。每批次全齡期蠶病調(diào)查數(shù)據(jù)為四齡眠期、五齡起蠶和五齡熟蠶期（上蔟第1 d）的血液型膿病發(fā)病率之和。

1. 2 氣象數(shù)據(jù)

氣象數(shù)據(jù)來源于中國(guó)天氣網(wǎng)，主要采集2008—2018年每年2—10月的旬平均氣溫、旬降雨量和旬日照時(shí)數(shù)等。

1. 3 數(shù)據(jù)整理

為更準(zhǔn)確地反映氣象因素對(duì)血液型膿病的影響，所引用數(shù)據(jù)不包括使用抗血液型膿病蠶品種情況下的蠶病采樣數(shù)據(jù)。蠶病數(shù)據(jù)與氣象因素時(shí)間匹配方案：五齡蠶期至上蔟第1 d所在某月某旬最長(zhǎng)的，即計(jì)入該旬，匹配對(duì)應(yīng)的氣象資料進(jìn)行分析。因家蠶血液型膿病具有極端氣候即在高溫和低溫發(fā)生概率較高的特點(diǎn)，故設(shè)置氣溫轉(zhuǎn)換算法，算法如下：

相關(guān)氣溫=（平均氣溫-某固定值）2

1. 4 蠶病程度等級(jí)劃分

為便于氣象數(shù)據(jù)預(yù)測(cè)預(yù)報(bào)，結(jié)合歷史調(diào)查數(shù)據(jù)的主要發(fā)病區(qū)間，將家蠶血液型膿病發(fā)病率劃分為4個(gè)等級(jí)，其劃分標(biāo)準(zhǔn)見表1。

1. 5 統(tǒng)計(jì)分析與建模

由于蠶病發(fā)病率為百分?jǐn)?shù)，對(duì)數(shù)據(jù)進(jìn)行平方根的反正弦值轉(zhuǎn)換，以便更符合正態(tài)分布規(guī)律（宋根鑫，2014）。使用SPSS 19.0的單因素方差分析（One-way ANOVA）進(jìn)行數(shù)據(jù)顯著性差異分析，通過Excel 2007對(duì)數(shù)據(jù)進(jìn)行整理和轉(zhuǎn)換，并匹配對(duì)應(yīng)的氣象因子，再利用DPS v7.05對(duì)數(shù)據(jù)進(jìn)行逐步回歸分析，篩選對(duì)蠶病影響顯著的氣象因子，生成回歸方程;分析氣象因子對(duì)蠶病的作用方式及影響程度，再對(duì)逐步回歸模型中的主要影響因子進(jìn)行通徑分析。采用2008—2017年發(fā)病數(shù)據(jù)和相關(guān)氣象數(shù)據(jù)進(jìn)行建模，并以2018年的氣象數(shù)據(jù)驗(yàn)證模型的可行性。

2 結(jié)果與分析

2. 1 不同年份家蠶血液型膿病發(fā)生規(guī)律分析

2008—2018年廣西家蠶血液型膿病平均發(fā)病率和危害等級(jí)如表2和圖1所示。從圖1可看出，廣西家蠶血液型膿病多年發(fā)病規(guī)律呈波動(dòng)態(tài)勢(shì)，其中2008年平均發(fā)病率最高，達(dá)13.24%，2012年平均發(fā)病率最低，為10.01%;2008—2018年年均發(fā)病率為11.35%，年發(fā)病率大于平均值的年份有6年（2008年、2011年、2014—2017年），占54.55%，家蠶血液型膿病發(fā)生率屬偏高。

以旬為單元，統(tǒng)計(jì)分析各危害等級(jí)的出現(xiàn)頻率，結(jié)果發(fā)現(xiàn)其比例約為Ⅰ級(jí)∶Ⅱ級(jí)∶Ⅲ級(jí)∶Ⅳ級(jí)=1∶4.5∶3.5∶3.4，出現(xiàn)頻率最高的是Ⅱ級(jí)，占36.36%，其次是Ⅲ級(jí)，占28.28%，第三是Ⅳ級(jí)，占27.27%;發(fā)病率大于10.00%的出現(xiàn)頻率達(dá)55.56%，說明廣西家蠶血液型膿病發(fā)生及危害較嚴(yán)重。

2. 2 不同月份家蠶血液型膿病發(fā)生規(guī)律分析

以月為單元，對(duì)家蠶血液型膿病發(fā)病率（表2）進(jìn)行分析，從圖2可看出，廣西家蠶血液型膿病發(fā)病率呈正態(tài)分布，即兩頭低中間高。經(jīng)單因素方差分析顯示，僅6和8月的發(fā)病率間差異不顯著（P>0.05，下同），其他月份差異均達(dá)顯著水平（P<0.05，下同）。病害發(fā)生最嚴(yán)重的月份是7月，平均發(fā)病率為18.04%;其次是6和8月，平均發(fā)病率分別為16.00%和15.36%。危害程度Ⅳ級(jí)的批次全部來自以上3個(gè)月份，且Ⅳ級(jí)出現(xiàn)頻率占6—8月調(diào)查總批次的62.79%，說明6—8月是廣西家蠶血液型膿病發(fā)病高峰期，危害嚴(yán)重。而3—5月和9—10月發(fā)病率總體處于較低水平，危害程度Ⅱ級(jí)以內(nèi)約占70.07%。

2. 3 不同蠶區(qū)不同季節(jié)家蠶血液型膿病發(fā)生規(guī)律分析

整理南寧、河池、柳州、貴港和百色五大蠶區(qū)2008—2018年春、夏、秋不同季節(jié)家蠶血液型膿病發(fā)生數(shù)據(jù)，分析其發(fā)病規(guī)律。從圖3可看出，春蠶期（3—5月）各蠶區(qū)發(fā)病率絕大部分在Ⅱ級(jí)以內(nèi)，由高至低依次為柳州（10.13%）>河池（9.65%）>百色（7.43%）>貴港（7.00%）>南寧（6.68%），經(jīng)單因素方差分析顯示，柳州和河池蠶區(qū)與南寧、百色、貴港蠶區(qū)差異達(dá)顯著水平。桂中、桂西北（柳州和河池）蠶區(qū)的發(fā)病率較其他蠶區(qū)偏高，可能與桂中、桂西北在春季低溫陰雨天氣偏多，桑葉質(zhì)量差，同時(shí)該季節(jié)常出現(xiàn)倒春寒天氣，氣溫劇烈波動(dòng)從而影響蠶生命力、抵抗力等因素有關(guān)，春季養(yǎng)蠶應(yīng)加強(qiáng)升溫防寒。

夏蠶期（6—8月）各蠶區(qū)發(fā)病率絕大部分達(dá)Ⅳ級(jí)，由高至低依次為貴港（18.15%）>南寧（17.89%）>河池（16.27%）>百色（15.82%）>柳州（14.16%）。經(jīng)單因素方差分析發(fā)現(xiàn)，貴港蠶區(qū)與柳州、百色、河池3個(gè)蠶區(qū)及南寧蠶區(qū)與柳州、百色2個(gè)蠶區(qū)的差異達(dá)顯著水平。地區(qū)發(fā)病率呈現(xiàn)出桂東南（南寧和貴港）高于桂中、桂西北（柳州、河池和百色），可能是桂東南地區(qū)夏季的高溫多濕氣候?qū)πQ生理活動(dòng)產(chǎn)生負(fù)面影響，且有利于病原微生物的繁殖和傳播，這一時(shí)段養(yǎng)蠶應(yīng)加強(qiáng)消毒防病及飼養(yǎng)抗病抗逆的蠶品種。

秋蠶期（9—10月）各蠶區(qū)發(fā)病率均在Ⅱ級(jí)以內(nèi)，由高至低依次為柳州（9.12%）>河池（8.46%）>貴港（8.44%）>南寧（8.14%）>百色（7.87%）。經(jīng)單因素方差分析顯示，各蠶區(qū)差異不顯著，說明廣西秋季氣候較適合養(yǎng)蠶，可選用優(yōu)質(zhì)高產(chǎn)的蠶品種，以提高蠶繭質(zhì)量和產(chǎn)量。

2. 4 家蠶血液型膿病規(guī)律與氣象因子相關(guān)性分析

以表2中2008—2017年家蠶血液型膿病發(fā)病率為基礎(chǔ)，經(jīng)反正弦平方根轉(zhuǎn)化后，與對(duì)應(yīng)時(shí)段的旬平均氣溫、旬降雨量、旬日照時(shí)數(shù)等氣象因子及其不同組合的氣象數(shù)據(jù)等，運(yùn)用數(shù)理統(tǒng)計(jì)方法進(jìn)行相關(guān)性分析、逐步回歸分析和通徑分析，建立廣西家蠶血液型膿病發(fā)病規(guī)律與氣象因子的逐步回歸模型：y=0.15116+0.0037371x1+0.00018527x2+0.000051826x3-0.0000030415x4-0.00000082702x5+0.00063798x6。式中，x1為旬平均氣溫與23.75 ℃差值的二次方，x2為上一旬降水量，x3為上兩旬平均氣溫的二次方，x4為上兩旬日照時(shí)數(shù)的二次方，x5為前三旬降雨量之和的二次方，x6為近四旬降雨量之和?；貧w模型決定系數(shù)R2=0.8844，達(dá)極顯著水平（P=0.001），已滿足逐步回歸篩選顯著因子的要求（P<0.01），且各單因素均與發(fā)病率顯著相關(guān)。

采用通徑分析獲得各因素對(duì)發(fā)病率的貢獻(xiàn)值（表3），6個(gè)構(gòu)成因素直接通徑系數(shù)分別為py1=0.4445、py2=0.1010、py3=0.0996、py4=-0.0713、py5=-0.3531和py6=0.8254，即單因素影響前兩位的分別是近四旬降雨量之和和旬平均氣溫與23.75 ℃差值的二次方，同時(shí)這兩個(gè)因素與發(fā)病率的偏相關(guān)分析也達(dá)極顯著水平（P<0.01）。參考各因素對(duì)發(fā)病率的影響系數(shù)總量，可知降雨量和氣溫是影響廣西蠶區(qū)家蠶血液型膿病發(fā)病率的最主要?dú)庀笠蜃樱罢邔?duì)發(fā)病率的直接影響程度略高于后者，日照時(shí)數(shù)與家蠶血液型膿病發(fā)病率呈弱負(fù)相關(guān)。此外，通徑分析的決定系數(shù)R2=0.8844，誤差e對(duì)因變量的通徑系數(shù)為0.3400，即仍有34.00%的未知因素未能有效模擬，說明家蠶血液型膿病的發(fā)生除了與氣象因子相關(guān)外，可能還與當(dāng)?shù)仞B(yǎng)蠶環(huán)境、技術(shù)水平及病原積累等因素有關(guān)，但具體原因有待進(jìn)一步探究。

該回歸模型對(duì)2008—2017年廣西蠶病發(fā)病率轉(zhuǎn)換等級(jí)的擬合精度為78.89%（圖4）。為驗(yàn)證該預(yù)測(cè)模型的可靠性，利用該模型對(duì)2018年數(shù)據(jù)進(jìn)行回檢，將預(yù)測(cè)發(fā)病率轉(zhuǎn)換為發(fā)病等級(jí)，結(jié)果（表4）顯示，可準(zhǔn)確預(yù)測(cè)2018年9批次蠶病發(fā)生等級(jí)中的6批次，擬合精度為66.67%，另外3批次可預(yù)測(cè)到相鄰等級(jí)，說明該蠶病預(yù)測(cè)模型達(dá)實(shí)用水平。對(duì)照預(yù)測(cè)有偏差的數(shù)據(jù)，當(dāng)觀測(cè)值接近等級(jí)分界邊緣時(shí)，容易出現(xiàn)防控等級(jí)預(yù)測(cè)偏差。

為進(jìn)一步驗(yàn)證該模型的適應(yīng)性，對(duì)目前廣西蠶繭產(chǎn)量排名前三的河池、南寧和百色蠶區(qū)分別進(jìn)行預(yù)測(cè)分析，結(jié)果（表4）顯示，2018年家蠶血液型膿病危害等級(jí)的擬合精度均達(dá)77.78%。表明建立的回歸模型對(duì)具體蠶區(qū)預(yù)報(bào)也有較高的準(zhǔn)確性。

3 討論

2008—2018年廣西家蠶血液型膿病的年均發(fā)病率呈波動(dòng)變化態(tài)勢(shì)，年均發(fā)生率為11.35%，危害等級(jí)在重度（Ⅲ級(jí)）和超重度（Ⅳ級(jí)）的出現(xiàn)頻率達(dá)55.56%，說明廣西家蠶血液型膿病發(fā)生和危害較嚴(yán)重，應(yīng)引起相關(guān)部門和蠶農(nóng)的重視，采取相應(yīng)的防控措施。廣西家蠶血液型膿病發(fā)生Ⅳ級(jí)危害幾乎都在夏蠶期（6—8月），其中該病發(fā)病率最嚴(yán)重的月份是7月，其次是6和8月，發(fā)病高峰期與廣西高溫多濕的季節(jié)基本重合，與閉立輝（2005）、關(guān)俊華等（2015）的報(bào)道基本一致。這個(gè)季節(jié)養(yǎng)蠶應(yīng)是全年防控的重點(diǎn)時(shí)期，除了加強(qiáng)消毒防病、溫濕度調(diào)節(jié)等技術(shù)措施外，還應(yīng)選擇適宜當(dāng)?shù)貧夂蛱攸c(diǎn)的蠶品種或抗性品種進(jìn)行飼養(yǎng)，以降低家蠶血液型膿病的損失。而3—5月（春蠶期）和9—10月（秋蠶期）發(fā)病率總體處于較低水平，蠶農(nóng)可在這兩個(gè)時(shí)期飼養(yǎng)優(yōu)質(zhì)高產(chǎn)蠶品種以提高蠶繭產(chǎn)量和質(zhì)量，從而獲得更大收益。

通過對(duì)廣西家蠶血液型膿病發(fā)生數(shù)據(jù)和氣象因子進(jìn)行相關(guān)性分析和通徑分析，明確了降雨量和氣溫是影響該病發(fā)生的主要?dú)庀笠蛩?，且前者的影響略高于后者，與關(guān)俊華等（2015）的研究結(jié)果略有差異。其原因除了研究對(duì)象的地域差異外，蠶的適溫范圍（20～28 ℃）較廣也不容忽視，即氣溫在蠶適宜溫度范圍內(nèi)變化對(duì)蠶病的影響小于其在適宜溫度范圍之外的影響，呈非線性關(guān)系。以廣西夏蠶期降雨量相近的河池和柳州（672.9和602.0 mm）為例，河池的平均氣溫27.96 ℃明顯高于柳州的25.81 ℃，家蠶血液型膿病發(fā)病率趨勢(shì)與溫度呈顯著正相關(guān);春蠶期柳州和貴港降雨量相近（436.1和401.4 mm），但柳州的平均氣溫18.50 ℃明顯低于貴港的22.95 ℃，家蠶血液型膿病發(fā)病率反而顯著高于貴港;秋蠶期降雨量相近的河池和貴港（247.7和239.2 mm），前者平均氣溫24.34 ℃明顯低于后者的26.21 ℃，但兩者的發(fā)病率無顯著差異。針對(duì)上述情況，本研究引入了相關(guān)氣溫的算法輔助分析氣溫對(duì)家蠶血液型膿病的影響，較好地解決了逐步回歸中血液型膿病發(fā)病率與氣溫的線性關(guān)系。

在預(yù)測(cè)模型中可知，2008—2010年的觀測(cè)值常大于測(cè)報(bào)模型擬合值，說明當(dāng)前氣象因子并不能完全解釋家蠶血液型膿病發(fā)病的原因，還存在其他未知因素，如病原致病力差異（黃深惠等，2012;唐芬芬等，2018）、蠶品種抗病性差異（陸瑞好等，2007）、蠶品種與當(dāng)?shù)貧夂虿贿m應(yīng)及養(yǎng)蠶消毒不徹底病原積累（閉立輝，2005）等。根據(jù)通徑分析的結(jié)果，近四旬降雨量之和和旬平均氣溫與23.75 ℃差值的二次方（相關(guān)氣溫）是血液型膿病防控中應(yīng)重點(diǎn)關(guān)注的氣象因素。在實(shí)際應(yīng)用過程中，可利用本研究的預(yù)測(cè)預(yù)報(bào)模型，以近1個(gè)月的氣象資料和中長(zhǎng)期天氣預(yù)報(bào)信息為基礎(chǔ)，能提供未來一段時(shí)間家蠶血液型膿病的危害等級(jí)預(yù)測(cè)，為該病精準(zhǔn)防控提供指導(dǎo)依據(jù)。

4 結(jié)論

廣西家蠶血液型膿病危害較嚴(yán)重，6—8月是發(fā)病高峰期。夏蠶期廣西家蠶血液型膿病危害等級(jí)大部分達(dá)Ⅳ級(jí)，呈現(xiàn)桂東南蠶區(qū)比桂中桂西北蠶區(qū)高的趨勢(shì)，建議這一時(shí)期選用抗病抗逆蠶品種飼養(yǎng);春蠶期和秋蠶期的危害等級(jí)大部分在Ⅱ級(jí)以內(nèi)，可選用優(yōu)質(zhì)高產(chǎn)的蠶品種，以提高質(zhì)量和產(chǎn)量。降雨量和氣溫是影響廣西家蠶血液型膿病發(fā)病率的最主要?dú)庀笠蜃印Ｔ趶V西蠶桑生產(chǎn)中可應(yīng)用建立的預(yù)測(cè)模型對(duì)家蠶血液型膿病發(fā)生危害進(jìn)行預(yù)測(cè)預(yù)報(bào)。

參考文獻(xiàn)：

白興榮，江亞，黃平. 2011. 云南省家蠶血液型膿病的危害與流行分析[J]. 云南農(nóng)業(yè)科技，40（4）：31-32. [Bai X R，Jiang Y，Huang P. 2011. Harm and epidemic analysis of Bombyx mori nuclear polyhedrosis in silkworm in Yunnan Province[J]. Yunnan Agricultural Science and Technology，40（4）：31-32.]

閉立輝. 2005. 廣西當(dāng)前主要蠶病與防治[J]. 廣西蠶業(yè)，42（4）：10-14. [Bi L H. 2005. Sillkworm diseases and control in Guangxi[J]. Guangxi Sericulture，42（4）：10-14.]

關(guān)俊華，高俊喜，齊淑彩，周文東. 2015. 南宮市蠶桑生產(chǎn)主要?dú)庀鬄?zāi)害分析[J]. 中國(guó)蠶業(yè)，36（4）：76-79. [Guan J H， Gao J X，Qi S C，Zhou W D. 2015. Analysis of main meteorological disasters of sericulture production in Nangong City[J]. China Sericulture，36（4）：76-79.]

韓正安，張曉紅. 2014. 春蠶后期低溫冷害的影響及防范措施[J]. 北方蠶業(yè)，35（2）：58-59. [Han Z A，Zhang X H. 2014. Effect of low temperature and cold injury in late spring silkworm rearing season and its preventive measures[J]. North Sericulture，35（2）：58-59.]

杭笑冬，仲偉榮，陳俊驥. 2017. 家蠶血液型膿病的發(fā)生原因與防治對(duì)策[J]. 江蘇蠶業(yè)，39（2）：12-13. [Hang X D，Zhong W R，Chen J J. 2017. Cause and control strategy of Bombyx mori nuclear polyhedrosis in silkworm[J]. Jiang-su Sericulture，39（2）：12-13.]

胡文娟. 2018. 廣西家蠶微粒子病流行病學(xué)研究[D]. 鎮(zhèn)江：江蘇科技大學(xué). [Hu W J. 2018. Epidemiologic study of silkworm pebrine disease in Guangxi Province[D]. Zhenjiang：Jiangsu University of Science and Technology.]

黃深惠，石美寧，唐亮，黃旭華. 2012. 廣西不同蠶區(qū)家蠶核型多角體病毒對(duì)家蠶致病力的研究[J]. 廣西蠶業(yè)，49（2）：17-19. [Huang S H，Shi M N，Tang L，Huang X H. 2012. Study on virulence of Bombyx mori nuclear polyhedrovirus to silkworm in different areas of Guangxi[J]. Guangxi Sericulture，49（2）：17-19.]

連武，陳謀宇，武有初. 1996. 氣象環(huán)境條件與養(yǎng)蠶效果的關(guān)系[J]. 四川蠶業(yè)，（2）：36-38. [Lian W，Chen M Y，Wu Y C. 1996. Relationship between meteorological environment and sericulture effect[J]. Sichuan Sericulture，（2）：36-38.]

劉吉平，張國(guó)權(quán)，徐興耀. 2002. 家蠶微粒子病流行發(fā)生與氣象因子的關(guān)系分析研究[J]. 農(nóng)業(yè)系統(tǒng)科學(xué)與綜合研究，18（3）：164-168. [Liu J P，Zhang G Q，Xu X Y. 2002. Ana-lytical study on the relation with pebrine disease prevalence to climatic factors[J]. System Sciences and Comprehensive Studies in Agriculture，18（3）：164-168.]

劉永裕，劉梅，唐國(guó)敏，吳炫柯，陳家福，李家文. 2010. 桂中地區(qū)原種蠶死籠率與氣象條件分析[J]. 安徽農(nóng)業(yè)科學(xué)，38（15）：8088-8090. [Liu Y Y，Liu M，Tang G M，Wu X K，Chen J F，Li J W. 2010. Analysis on the percentage of dead worm cocoon of original silkworm and weather conditions in central Guangxi[J]. Journal of Anhui Agricultural Sciences，38（15）：8088-8090.]

陸瑞好，石美寧，閉立輝，顧家棟，黃君霆. 2007. 廣西地區(qū)家蠶品種資源對(duì)BmNPV抵抗性的初步調(diào)查[J]. 蠶業(yè)科學(xué)，33（1）：117-120. [Lu R H，Shi M N，Bi L H，Gu J D，Huang J T. 2007. A preliminary investigation on the resistance of the silkworm variety resources in Guangxi to BmNPV[J]. Science of Sericulture，33（1）：117-120.]

魯興萌. 2012. 家蠶對(duì)血液型膿病的抗性與防治策略[J]. 中國(guó)蠶業(yè)，33（3）：4-7. [Lu X M. 2012. Resistance and control strategy of Bombyx mori nuclear polyhedrosis in silkworm[J]. China Sericulture，33（3）：4-7.]

牟志美，王彥文，崔為正，李衛(wèi)國(guó)，殷培偉. 1999. 家蠶真菌病發(fā)生及流行規(guī)律的研究[J]. 山東農(nóng)業(yè)大學(xué)學(xué)報(bào)，30（3）：199-207. [Mu Z M，Wang Y W，Cui W Z，Li W G，Yin P W. 1999. Study on the outbreak law of mycosis of silkworm Bombyx mori L.[J]. Journal of Shandong Agricultural University，30（3）：199-207.]

潘志新，黃旭華，湯慶坤，羅梅蘭，韋廷秀，蔣滿貴，何松濤，李標(biāo). 2017. 廣西家蠶微粒子病流行規(guī)律分析[J]. 廣西蠶業(yè)，54（1）：6-10. [Pan Z X，Huang X H，Tang Q K，Luo M L，Wei T X，Jiang M G，He S T，Li B. 2017. Analysis on the epidemic law of the pebrine disease of the silkwom in Guangxi[J]. Guangxi Sericulture，54（1）：6-10.]

祁廣軍. 2013. “東桑西移”背景下廣西蠶絲業(yè)發(fā)展實(shí)證研究[D]. 北京：中國(guó)農(nóng)業(yè)大學(xué). [Qi G J. 2013. An enpirical study on the development of Guangxi silk industry in the context of ‘sericulture moving from east to west[D]. Beijing：China Agricultural University.]

石美寧，唐亮，黃紅燕，韋廷秀，潘志新，湯慶坤，祁廣軍，浦月霞，莫云霞，賈雪峰. 2015. 桑蠶新品種“桂蠶N2”的推廣應(yīng)用初報(bào)[J]. 廣西蠶業(yè)，52（3）：27-31. [Shi M N，Tang L，Huang H Y，Wei T X，Pan Z X，Tang Q K，Qi G J，Pu Y X，Mo Y X，Jia X F. 2015. A preliminary report on the popularization and application of the new silkworm varie-ty “Guican N2”[J]. Guangxi Sericulture，52（3）：27-31.]

宋根鑫. 2014. 縣域農(nóng)田養(yǎng)分動(dòng)態(tài)監(jiān)測(cè)的方法優(yōu)化和應(yīng)用研究[D]. 杭州：浙江大學(xué). [Song G X. 2014. County farmland nutrient dynamic monitoring method optimization and application research[D]. Hangzhou：Zhejiang University.]

唐芬芬，張永紅，邵榆嵐，朱峰，白興榮. 2018. 云南家蠶核型多角體病毒分離株的毒力測(cè)定及系統(tǒng)進(jìn)化分析[J]. 昆蟲學(xué)報(bào)，61（10）：1160-1169. [Tang F F，Zhang Y H，Shao Y L，Zhu F，Bai X R. 2018. Virulence and phylogenetic analysis of Bombyx mori nucleopolyhedrovirus isolates from Yunnan，southwestern China[J]. Acta Entomologica Sinica，61（10）： 1160-1169.]

王必華. 2015. 家蠶血液型膿病的病因與防治對(duì)策[J]. 現(xiàn)代農(nóng)業(yè)科技，44（11）：284. [Wang B H. 2015. Etiology and prevention of Bombyx mori nuclear polyhedrosis in silkworm[J]. Modern Agricultural Science and Technology，44（11）：284.]

肖建京，余榮峰，豐洪平，項(xiàng)瓊斐，陳田飛，章朝凱，何曉峰，何荷根. 2015. 淳安縣2014年晚秋蠶病大面積暴發(fā)原因的調(diào)查分析[J]. 中國(guó)蠶業(yè)，36（2）：39-43. [Xiao J J，Yu R F，F(xiàn)eng H P，Xiang Q F，Chen T F，Zhang C K，He X F，He H G. 2015. Investigation and analysis on the causes of large-scale outbreak of silkworm disease in late autumn silkworm rearing season of 2014 in Chunan County[J]. China Sericulture，36（2）：39-43.]

（責(zé)任編輯 羅 麗）