陳 聰,徐菀鋆
(福建農(nóng)林大學(xué),福建 福州 350002)
據(jù)統(tǒng)計(jì),我國(guó)2012-2021 年近十年期間發(fā)生的火災(zāi)超過335 萬起,因火災(zāi)造成的死亡人數(shù)和受傷人數(shù)分別高達(dá)15 748 人和11 547 人,直接經(jīng)濟(jì)損失逼近420 億元[1]。影響火災(zāi)發(fā)生的因素有很多,如社會(huì)經(jīng)濟(jì)因素、自然因素等,氣象因素是主要的自然因素之一。研究氣象因素與火災(zāi)次數(shù)之間的相關(guān)關(guān)系,有助于在了解氣象相關(guān)數(shù)據(jù)時(shí),對(duì)于火災(zāi)的發(fā)生進(jìn)行預(yù)判,從而降低火災(zāi)發(fā)生的概率,有效地進(jìn)行火災(zāi)防治。
目前對(duì)城市火災(zāi)次數(shù)與氣象因素之間關(guān)系的研究主要集中在我國(guó)北方城市,南方的較少。黃韜等運(yùn)用自回歸模型,研究了北京市月度火災(zāi)起數(shù)與氣象因素之間的關(guān)系[2];譚常春等運(yùn)用逐步回歸模型和Adaptive-Lasso 分析方法研究了天津市月度火災(zāi)發(fā)生次數(shù)之間與氣象因子間的關(guān)系[3-4]。南北方氣候差異顯著,對(duì)火災(zāi)次數(shù)的影響也不一樣。本文將運(yùn)用最小二乘回歸和分位數(shù)回歸對(duì)福州市月度火災(zāi)次數(shù)與氣象因素之間的相關(guān)性進(jìn)行研究。
在傳統(tǒng)的多元線性回歸模型中,使用殘差平方和最小準(zhǔn)則即最小二乘法估計(jì)模型的回歸參數(shù)。
隨機(jī)抽取容量為n 的樣本觀察值{(x1i,x2i,…,yi):i:1,2,…,n},多元線性回歸模型可表示為:
根據(jù)最小二乘原理,參數(shù)估計(jì)值應(yīng)使
達(dá)到最小,求解待估參數(shù)的正規(guī)方程組就能得到
分位數(shù)回歸模型[5]可表示為:
其中,Y=(y1,y2,…,yn)' 為被解釋變量,X=(1,x1,x2,…,xk)'為解釋變量,βτ=(β0,β1,…,βk)為參數(shù)向量,ετ=(ε1τ,ε2τ,…,εnτ)'為對(duì)應(yīng)的誤差向量,τ(0<τ<1)表示特定的分位數(shù)。
分位數(shù)回歸通過給定被解釋變量Y 的分位數(shù)τ,然后利用加權(quán)殘差和最小求參數(shù)估計(jì),具體表示如下:
其中,ρτ=μ(τ-I(μ))是損失函數(shù),I(μ)為示性函數(shù),即
從而得到Y(jié) 的條件分位數(shù)估計(jì)為
本文火災(zāi)次數(shù)數(shù)據(jù)來源于福州市人民政府(http://www.fuzhou.gov.cn/),氣象因素?cái)?shù)據(jù)來源于中國(guó)氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn/)的地面資料中數(shù)據(jù)與產(chǎn)品中的中國(guó)地面國(guó)際交換站氣候資料月值數(shù)據(jù)集。
從圖1 可以看出,福州市2017-2019 年間的月度火災(zāi)次數(shù)有離群點(diǎn),為避免其影響模型擬合結(jié)果,將3個(gè)異常數(shù)據(jù)剔除。對(duì)所有數(shù)據(jù)取自然對(duì)數(shù),得到月火災(zāi)對(duì)數(shù)次數(shù)lnN、月平均對(duì)數(shù)風(fēng)速lnM、月平均對(duì)數(shù)氣溫lnW、月平均對(duì)數(shù)降水量lnP 及月平均對(duì)數(shù)相對(duì)濕度lnH。
圖1 福州市2017-2019 年月火災(zāi)次數(shù)箱線圖
考慮到月火災(zāi)次數(shù)可能與某些氣象因素的前期值有關(guān),引入滯后變量[2]。對(duì)滯后0、1、2 期的變量進(jìn)行Spearman 相關(guān)系數(shù)檢驗(yàn),根據(jù)表1 結(jié)果選擇對(duì)月平均相對(duì)濕度進(jìn)行滯后一期處理。最后建立回歸模型如下:
表1 滯后變量的Spearman 相關(guān)系數(shù)檢驗(yàn)
對(duì)模型(8)進(jìn)行最小二乘回歸,回歸結(jié)果和擬合效果分別見表2 和圖2。
表2 最小二乘回歸結(jié)果
圖2 最小二乘回歸擬合圖
回歸模型F 檢驗(yàn)的p 值為0.000 7,說明模型的回歸效果是顯著的。從圖2 可以看出,模型的擬合效果較好,可用于解釋福州市月火災(zāi)次數(shù)與氣象因素之間的相關(guān)關(guān)系。從表2 可以看出,氣溫、降水量、相對(duì)濕度這三個(gè)因素的系數(shù)都通過了顯著性檢驗(yàn),說明它們對(duì)福州市月火災(zāi)次數(shù)的影響是顯著的,且月火災(zāi)次數(shù)與風(fēng)速、溫度和降水量呈負(fù)相關(guān)關(guān)系;與相對(duì)濕度呈現(xiàn)正相關(guān)關(guān)系。從回歸系數(shù)的估計(jì)值可以看出,在各解釋變量中最為敏感的是相對(duì)濕度,當(dāng)上月相對(duì)濕度減少兩個(gè)單位時(shí),本月的月火災(zāi)次數(shù)平均減少1起。模型中系數(shù)都較小,說明火災(zāi)次數(shù)還受到其他因素的影響,但火災(zāi)次數(shù)和氣象因素之間存在密不可分的關(guān)系。
選擇0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8、0.9 九個(gè)分位點(diǎn)對(duì)模型(8)進(jìn)行分位數(shù)回歸,結(jié)果見表3。
表3 分位數(shù)回歸結(jié)果
從表3 可以看出,風(fēng)速對(duì)福州市月火災(zāi)次數(shù)的影響是不顯著的。但可以看到,風(fēng)速在不同分位點(diǎn)處對(duì)月火災(zāi)次數(shù)的影響還是有差異的,在高分位點(diǎn)對(duì)于福州市月火災(zāi)次數(shù)的影響從絕對(duì)值方面看要高于在低分位點(diǎn)上的影響,且整體與福州市月火災(zāi)次數(shù)呈現(xiàn)負(fù)相關(guān)關(guān)系。
氣溫與福州市月火災(zāi)次數(shù)呈負(fù)相關(guān)關(guān)系,在0.3~0.5 分位點(diǎn)處對(duì)福州市月火災(zāi)次數(shù)的影響顯著,說明月平均氣溫在溫度較高即春季即將進(jìn)入夏季,夏季即將進(jìn)入秋季時(shí),對(duì)于福州市月火災(zāi)次數(shù)有顯著的影響。在這兩個(gè)較高溫季節(jié),是溫度、濕度都在轉(zhuǎn)變的季節(jié),易燃物所處環(huán)境容易受到突然的變化。在春季即將進(jìn)入夏季時(shí),雖然溫度升高,但濕度也升高,所以對(duì)于火災(zāi)次數(shù)有一定的抑制作用。在夏季進(jìn)入秋季時(shí),溫度保持一定的“高度”,濕度也在增加,所以對(duì)于火災(zāi)次數(shù)也存在抑制作用。
降水量與福州市月火災(zāi)次數(shù)呈負(fù)相關(guān)關(guān)系,當(dāng)降水量越大時(shí),火災(zāi)所處的環(huán)境便會(huì)較為濕潤(rùn),故火災(zāi)次數(shù)較少。相較于其他影響顯著的因素,降水量對(duì)福州市月火災(zāi)次數(shù)的影響稍小,僅在0.3 分位點(diǎn)影響顯著,這可能與福州的地理位置所處的氣候帶常年雨水充足有關(guān)。
相對(duì)濕度在0.3~0.8 分位點(diǎn)處對(duì)福州市月火災(zāi)次數(shù)均有顯著影響,在中位數(shù)處的影響最大。不管在火災(zāi)高發(fā)期還是低發(fā)期,相對(duì)濕度都對(duì)于火災(zāi)次數(shù)有顯著的正相關(guān)影響。福州容易發(fā)生“回南天”的現(xiàn)象,出現(xiàn)“回南天”時(shí),相對(duì)濕度變高,容易發(fā)生電路短路火災(zāi)和生產(chǎn)經(jīng)營(yíng)性火災(zāi),使得火災(zāi)的發(fā)生概率提高。模型中的濕度因素為上月平均相對(duì)濕度,當(dāng)上月平均相對(duì)濕度低時(shí),容易引起相關(guān)部門的注意,從而對(duì)于易燃物及其所處的環(huán)境更加重視,進(jìn)而使本月的月火災(zāi)次數(shù)得到降低。
截距可以認(rèn)為是當(dāng)各個(gè)氣象因素值均為零時(shí),福州市月火災(zāi)次數(shù)所受到的影響。經(jīng)計(jì)算,截距的95%置信區(qū)間均為正值,說明福州市月火災(zāi)次數(shù)不僅受到本文提到的四種氣象因素的影響,還受到其他氣象因素或其他因素的影響。
本文分別運(yùn)用最小二乘回歸和分位數(shù)回歸研究了月平均風(fēng)速、月平均氣溫、月平均降水量、月平均相對(duì)濕度等氣象因素對(duì)福州市月火災(zāi)次數(shù)的影響。兩種方法比較發(fā)現(xiàn),分位數(shù)回歸不僅能夠反映出最小二乘回歸所得到的各氣象因素與火災(zāi)次數(shù)的相關(guān)性,還可以反映出不同分位點(diǎn)上各因素對(duì)火災(zāi)次數(shù)的影響。
結(jié)果表明,氣象因素中月平均氣溫、月平均降水量、月平均相對(duì)濕度對(duì)火災(zāi)次數(shù)影響顯著,而月平均風(fēng)速對(duì)火災(zāi)次數(shù)的影響不大。