基于多元統(tǒng)計(jì)分析絕緣紙工頻擊穿強(qiáng)度影響因素的研究

吳東樂(lè) 劉群華 孫勝然 徐凱麗 劉 文，*

（1.中國(guó)制漿造紙研究院有限公司，北京，100102；2.國(guó)家紙張質(zhì)量監(jiān)督檢驗(yàn)中心，北京，100102；3.制漿造紙國(guó)家工程實(shí)驗(yàn)室，北京，100102）

隨著特高壓輸變電工程的大規(guī)模建設(shè)，電工絕緣材料的需求呈現(xiàn)階段式增長(zhǎng)，作為一種優(yōu)良的絕緣材料——纖維素絕緣紙的性能及其油紙復(fù)合絕緣體系性能的提升面臨巨大挑戰(zhàn)［1］。此外，由于油紙絕緣體系的復(fù)雜性，能否獲得更準(zhǔn)確的影響機(jī)制和計(jì)算方法，以縮短電力裝備的研發(fā)周期并降低成本，是研究人員長(zhǎng)期面臨的重要問(wèn)題。工頻擊穿強(qiáng)度是絕緣紙的宏觀電氣性能，這與纖維形態(tài)和紙張結(jié)構(gòu)密切關(guān)聯(lián)，準(zhǔn)確分析絕緣紙工頻擊穿強(qiáng)度的影響因素，并對(duì)擊穿強(qiáng)度進(jìn)行預(yù)測(cè)和工藝優(yōu)化研究，有助于提升特高壓變壓器用絕緣紙產(chǎn)品的制造技術(shù)和質(zhì)量水平。

植物纖維形態(tài)包括纖維長(zhǎng)度、纖維寬度、纖維壁厚和細(xì)胞腔直徑等參數(shù)［2］，纖維形態(tài)對(duì)絕緣紙的多種性能影響顯著。黃建文等［3］研究針葉木原漿種類對(duì)絕緣紙擊穿性能的影響發(fā)現(xiàn)，不同種類針葉木漿的纖維形態(tài)不同，相應(yīng)制得的絕緣紙?jiān)谟椭械膿舸?qiáng)度也不同。Zhou等［4］采用多元統(tǒng)計(jì)分析方法研究纖維形態(tài)各項(xiàng)參數(shù)對(duì)直流電場(chǎng)下絕緣紙擊穿強(qiáng)度的影響；結(jié)果表明，纖維長(zhǎng)度和細(xì)小纖維含量對(duì)其擊穿強(qiáng)度影響顯著。

從絕緣紙的擊穿理論來(lái)說(shuō)，絕緣紙是由纖維相互交織、層疊而成、具有三維多孔結(jié)構(gòu)的材料，在纖維網(wǎng)絡(luò)內(nèi)部有大量的孔隙；因此，無(wú)論以油浸漬還是以原有狀態(tài)進(jìn)行擊穿強(qiáng)度實(shí)驗(yàn)，實(shí)際上都是對(duì)由固體纖維和液體或氣體等絕緣介質(zhì)組成的混合結(jié)構(gòu)進(jìn)行的實(shí)驗(yàn)［5］。因此，在研究絕緣紙擊穿強(qiáng)度性能的影響因素時(shí)，必須考慮這種混合結(jié)構(gòu)。紙張結(jié)構(gòu)可以用緊度、孔隙率、透氣度和孔徑分布等指標(biāo)進(jìn)行描述，這些指標(biāo)從不同的角度描述紙張結(jié)構(gòu)特性［6］。針對(duì)紙張的多孔結(jié)構(gòu)，Mohanty等［7-10］通過(guò)在薄膜上制造出一個(gè)圓形孔，將其放在絕緣紙上來(lái)模擬絕緣紙的孔隙局部放電，采用模糊邏輯方法、多層前饋網(wǎng)絡(luò)、徑向基函數(shù)網(wǎng)絡(luò)和人工神經(jīng)網(wǎng)絡(luò)等軟件計(jì)算技術(shù)，分別對(duì)幾類絕緣紙的擊穿電壓進(jìn)行建模，并針對(duì)直流和交流電壓條件提出了不同的預(yù)測(cè)模型；結(jié)果表明，建立的預(yù)測(cè)模型能夠較好預(yù)測(cè)絕緣紙的擊穿電壓。Ghosh等［11］則通過(guò)人工神經(jīng)網(wǎng)絡(luò)模型研究了絕緣材料的孔隙尺寸和厚度與其局部放電之間的關(guān)系。以上研究都是基于在薄膜上制造一個(gè)圓形孔，其大小取決于所制造孔的直徑和薄膜的厚度；但實(shí)際上，絕緣紙內(nèi)部的孔隙具有不同的幾何形狀，并且分布隨機(jī)。所以，以上研究與實(shí)際情況存在差異。Singh等［12］以厚度、相對(duì)介電常數(shù)、損耗角正切和體積電阻率為參數(shù)，利用人工神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)纖維素絕緣材料的擊穿強(qiáng)度，取得良好的預(yù)測(cè)結(jié)果，但研究存在神經(jīng)網(wǎng)絡(luò)模型運(yùn)用較為復(fù)雜的問(wèn)題。

基于變壓器體系特性，以及各影響因素與擊穿強(qiáng)度之間的復(fù)雜關(guān)系，可以采用灰色關(guān)聯(lián)分析（GRA）量化各影響因素與擊穿強(qiáng)度之間的關(guān)系。王雁凌等［13-14］采用GRA方法構(gòu)建了一種模型，模擬電力系統(tǒng)中電力負(fù)荷和輸電線路的損耗與各影響因素之間的關(guān)系。

本課題根據(jù)GRA理論，研究纖維長(zhǎng)度、纖維寬度、細(xì)小纖維含量、厚度、緊度、孔隙率和透氣度與絕緣紙工頻擊穿強(qiáng)度的灰色關(guān)聯(lián)度，構(gòu)建工頻擊穿強(qiáng)度影響因素量化分析模型，明確各因素對(duì)工頻擊穿強(qiáng)度的影響機(jī)制；其次，采用多元統(tǒng)計(jì)分析方法中的主成分分析和最佳子集選擇方法構(gòu)建工頻擊穿強(qiáng)度的預(yù)測(cè)模型，可以對(duì)產(chǎn)品性能進(jìn)行預(yù)測(cè)，預(yù)先做出風(fēng)險(xiǎn)評(píng)估，并指導(dǎo)絕緣紙的研究和生產(chǎn)。

1 實(shí)驗(yàn)

1.1 實(shí)驗(yàn)原料

纖維原材料采用加拿大Canfor公司的未漂白針葉木硫酸鹽漿，灰分0.43%；用于浸漬的變壓器油選用北京欣環(huán)誠(chéng)油脂分裝有限公司的25#變壓器油。

1.2 樣品的制備

本課題通過(guò)正交實(shí)驗(yàn)，對(duì)打漿度、打漿負(fù)荷、定量與干燥方式4個(gè)因素進(jìn)行設(shè)計(jì)并制備相應(yīng)的絕緣紙樣品，實(shí)驗(yàn)方案如表1所示。將絕緣紙樣品置于105℃烘箱內(nèi)烘干4 h，并將變壓器油預(yù)熱至70℃，然后將變壓器油倒入盛有絕緣紙樣品的容器中，在真空條件下浸漬48 h，浸油完成后除去絕緣紙樣品表面多余的變壓器油，將樣品放置于密封袋中保存。

1.3 纖維形態(tài)和紙張結(jié)構(gòu)參數(shù)測(cè)試

采用瑞典L&W公司的纖維分析儀測(cè)定纖維的長(zhǎng)度、寬度和細(xì)小纖維含量；按照GB/T 451.3—2002測(cè)定絕緣紙的厚度和緊度；采用全自動(dòng)真密度分析儀測(cè)定絕緣紙的孔隙率；按照GB/T 458—2002測(cè)定絕緣紙的透氣度。

表1 實(shí)驗(yàn)方案Table 1 Experimental scheme

1.4 工頻擊穿強(qiáng)度測(cè)試

按照GB/T 3333—1999測(cè)定絕緣紙的工頻擊穿強(qiáng)度。具體為，采用不對(duì)稱黃銅電極，表面均進(jìn)行拋光處理，高壓電極直徑25 mm，低壓電極直徑75 mm；采用逐步均勻升壓的方式進(jìn)行測(cè)試，以測(cè)試數(shù)據(jù)的中位數(shù)作為最終實(shí)驗(yàn)結(jié)果。

1.5 GRA

GRA是一種基于灰色系統(tǒng)理論、定量描述和對(duì)比一個(gè)系統(tǒng)發(fā)展和變化的方法，是不確定系統(tǒng)研究領(lǐng)域的重要成果之一［15］。其基本思想是，將關(guān)注的研究對(duì)象及其影響因素看作一條曲線上的各點(diǎn)，參照曲線幾何相似性，即形狀越相似，兩個(gè)變量之間的關(guān)聯(lián)程度越高。通過(guò)計(jì)算多因素與同一參考序列之間的關(guān)聯(lián)度，確定主導(dǎo)因素。該方法適用于解決多因素之間復(fù)雜的相互關(guān)系，其優(yōu)點(diǎn)是不要求有大量的樣本數(shù)據(jù)，也不要求樣本數(shù)據(jù)有某種分布規(guī)律；具體過(guò)程如下［16］。

第一步，對(duì)獲得的原始數(shù)據(jù)樣本進(jìn)行數(shù)據(jù)轉(zhuǎn)換，使數(shù)據(jù)間可以進(jìn)行比較。這是因?yàn)椋饔绊懸蛩氐奈锢硪饬x不一樣，為了確保構(gòu)建模型的準(zhǔn)確性，采用區(qū)間值化轉(zhuǎn)換進(jìn)行數(shù)據(jù)轉(zhuǎn)換，分別對(duì)正、負(fù)序列采用下式處理。

第二步，確定絕對(duì)差序列：|x0（k）?xi（k）|。

第三步，根據(jù)式（3）計(jì)算灰色關(guān)聯(lián)系數(shù)（GRC）。

通過(guò)比較GRG的大小來(lái)量化影響因素對(duì)指標(biāo)的關(guān)聯(lián)程度；GRG越高，則相關(guān)程度越高，影響程度越大。

2 結(jié)果與討論

2.1 絕緣紙樣品的各項(xiàng)參數(shù)

絕緣紙樣品的纖維長(zhǎng)度、纖維寬度、細(xì)小纖維含量、厚度、緊度、孔隙率和透氣度7項(xiàng)參數(shù)測(cè)試結(jié)果如表2所示。為便于統(tǒng)計(jì)分析，將各參數(shù)映射到變量x1～x7。由表2可知，各項(xiàng)參數(shù)的變化幅度比較明顯。

表2 纖維形態(tài)和紙張結(jié)構(gòu)參數(shù)結(jié)果Table 2 Results of fiber morphology and paper structure parameters

絕緣紙的工頻擊穿強(qiáng)度數(shù)據(jù)如表3所示；其中，工頻擊穿強(qiáng)度（y）最大值為56.69 kV/mm，最小值為24.02 kV/mm，表明各絕緣紙樣品間工頻擊穿強(qiáng)度差異明顯。

表3 絕緣紙的工頻擊穿強(qiáng)度Table 3 AC breakdown strength of insulating paper samples

2.2 灰色關(guān)聯(lián)分析

將工頻擊穿強(qiáng)度定義為參考序列，各參數(shù)定義為比較序列，按照式（1）～式（4）計(jì)算得到各影響因素的GRG，結(jié)果如表4所示。由表4可知，各參數(shù)對(duì)絕緣紙工頻擊穿強(qiáng)度的GRG在0.524～0.797之間；按GRG的大小排序?yàn)椋豪w維長(zhǎng)度＞細(xì)小纖維含量＞纖維寬度＞緊度＞透氣度＞孔隙率＞厚度。

表4 灰色關(guān)聯(lián)度Table 4 Grey relational grades

圖1為絕緣紙擊穿的過(guò)程機(jī)理圖，反映了擊穿與電荷的關(guān)系［17］。由圖1可知，絕緣紙的擊穿往往是由內(nèi)部纖維間孔隙的局部放電引起的，更小的孔隙可以承受更高的電場(chǎng)。較短的纖維可以改善纖維間的排列，減小孔隙；細(xì)小纖維也會(huì)填補(bǔ)在纖維間的孔隙中，影響孔徑分布，從而降低絕緣紙的孔隙率；纖維寬度越小，緊度越大，紙張內(nèi)部的纖維排列越緊密，孔隙越小。因此，纖維越短、細(xì)小纖維含量越高、纖維寬度越小、緊度越大，均有助于提高絕緣紙的工頻擊穿強(qiáng)度。此外，與添加納米材料增強(qiáng)絕緣紙工頻擊穿強(qiáng)度的原理類似，細(xì)小纖維可以改善絕緣紙的表面均勻性，使絕緣紙表面電場(chǎng)均勻分布；同時(shí)增加絕緣紙內(nèi)部的陷阱密度，可以限制載流子的運(yùn)輸，使得絕緣紙內(nèi)部電子的有效平均自由行程變短，電子傳輸時(shí)不容易在電場(chǎng)中移動(dòng)，并轉(zhuǎn)變?yōu)樗俣容^慢的電子，從而減輕對(duì)纖維素分子鏈的轟擊，降低電離的概率，使得放電通道難以形成，在一定程度上抑制絕緣紙內(nèi)部放電，從而提高絕緣紙的工頻擊穿強(qiáng)度［18-21］。

圖1 絕緣紙擊穿機(jī)理Fig.1 Breakdown mechanism of insulating paper

2.3 預(yù)測(cè)模型的建立

2.3.1 小樣本和多重共線性問(wèn)題

多元統(tǒng)計(jì)分析對(duì)自變量和測(cè)試樣本的數(shù)量有要求，一般來(lái)說(shuō)，構(gòu)建多元回歸模型所需的樣本數(shù)量通常為自變量數(shù)量的2～5倍［4］。本課題研究中自變量有7個(gè)，而用來(lái)建立模型的模擬樣本數(shù)量只有12個(gè)。因此，需要解決樣本數(shù)量少的問(wèn)題。

在考慮多個(gè)自變量進(jìn)行研究時(shí)，當(dāng)各變量之間存在較強(qiáng)的關(guān)系時(shí)，進(jìn)行多元回歸分析時(shí)經(jīng)常產(chǎn)生多重共線性問(wèn)題，導(dǎo)致不符合多元線性回歸數(shù)學(xué)模型的基本假設(shè)。多重共線性會(huì)使模型的回歸系數(shù)失去實(shí)際意義，其正負(fù)號(hào)也可能出現(xiàn)倒置。皮爾遜相關(guān)系數(shù)適用于表示兩變量之間的相關(guān)性，其可量化不同變量之間的關(guān)聯(lián)程度。一般認(rèn)為，|r|≥0.8時(shí)，兩變量高度相關(guān)；0.5≤|r|＜0.8時(shí)，兩變量中度相關(guān)；0.3≤|r|＜0.5，兩變量低度相關(guān)；|r|＜0.3，兩變量不相關(guān)［22］。本課題中，自變量與因變量（共8個(gè)變量）的皮爾遜相關(guān)系數(shù)和雙尾顯著性概率以及變量之間的散點(diǎn)圖和置信橢圓，如圖2所示。散點(diǎn)圖和置信橢圓可以直觀地顯示變量之間的關(guān)聯(lián)程度。從圖2可以看出，纖維寬度、細(xì)小纖維含量和緊度與絕緣紙工頻擊穿強(qiáng)度呈正相關(guān)關(guān)系；纖維長(zhǎng)度、孔隙率和透氣度與絕緣紙工頻擊穿強(qiáng)度呈負(fù)相關(guān)關(guān)系。各參數(shù)與絕緣紙工頻擊穿強(qiáng)度的皮爾遜相關(guān)系數(shù)基本大于0.5，說(shuō)明它們之間具有很強(qiáng)的線性相關(guān)性；同樣，部分變量之間存在很強(qiáng)的相關(guān)性，如細(xì)小纖維含量與纖維寬度、纖維長(zhǎng)度和透氣度之間。由此推斷，可以使用多元線性回歸進(jìn)行絕緣紙工頻擊穿強(qiáng)度影響因素的分析，但各變量之間存在嚴(yán)重的多重共線性問(wèn)題，將所有參數(shù)作為自變量直接進(jìn)行回歸，建立的回歸方程模型是不合適的。

針對(duì)小樣本和多重共線性問(wèn)題，一方面可以增加樣本數(shù)量，但增加樣本數(shù)量有時(shí)是不現(xiàn)實(shí)的，如很多自變量不能控制，或者由于實(shí)際因素?zé)o法再獲得一些樣本數(shù)據(jù)；在一些情況下，當(dāng)自變量的個(gè)數(shù)較多時(shí)，在選擇增加樣本數(shù)量的情況下，增加什么樣的數(shù)據(jù)才能避免多重共線性問(wèn)題是無(wú)法確定的。另一方面可以降低自變量的數(shù)量，可以剔除解釋意義不強(qiáng)的自變量，直接選擇解釋意義非常強(qiáng)的自變量；還可以根據(jù)變量間的相關(guān)性對(duì)自變量進(jìn)行重新組合，提取較少的變量。

考慮到增加測(cè)試樣本的數(shù)量受到許多實(shí)際條件的限制，本課題選擇減少變量的數(shù)量，即從7個(gè)參數(shù)中選擇2～4個(gè)變量。主成分分析（PCA）適用于變量之間存在嚴(yán)重共線性的情況，在確保數(shù)據(jù)信息丟失最少的情況下，把多個(gè)變量轉(zhuǎn)化為更少的幾個(gè)新變量，是一種對(duì)多個(gè)變量進(jìn)行降維的統(tǒng)計(jì)處理方法［22］。提取的新成分是原始多個(gè)變量數(shù)據(jù)的線性組合，它們之間無(wú)相關(guān)性，可以有效地代表原變量，從而簡(jiǎn)化復(fù)雜問(wèn)題。最佳子集選擇法是基于枚舉的思想，對(duì)于給定數(shù)量的自變量遍歷所有可能的組合，并列出變量子集中具有最高擬合度的子集；然后根據(jù)變量組合的R2和誤差平方和（Mallows＇Cp）值選定參數(shù)的個(gè)數(shù)，可以確定適于多元線性回歸建模的參數(shù)集［4］。本課題選擇PCA和最優(yōu)子集選擇方法進(jìn)行變量的選擇，然后進(jìn)行多元線性建模。

2.3.2 PCA的多元線性回歸模型

在進(jìn)行計(jì)算前，需對(duì)所有變量進(jìn)行標(biāo)準(zhǔn)化變化，然后利用標(biāo)準(zhǔn)化變化的結(jié)果求解主成分會(huì)更加簡(jiǎn)便，標(biāo)準(zhǔn)化公式如下［22］。

圖2 變量的皮爾遜相關(guān)系數(shù)Fig.2 Pearson correlation coefficients among variables

樣本的充分性Kaiser-Meyer-Olkin（KMO）檢測(cè)是檢驗(yàn)變量之間偏相關(guān)的大小，它的取值范圍在0～1之間，數(shù)值越大越適合進(jìn)行PCA，小于0.5則不適合進(jìn)行PCA；巴特利特球形度檢驗(yàn)用于驗(yàn)證相關(guān)系數(shù)矩陣是不是單位矩陣，表明是不是適合進(jìn)行PCA［23-24］。KMO檢測(cè)和巴特利特球形度檢驗(yàn)的結(jié)果如表5所示。由表5可知，本課題樣本的KMO值為0.718，巴特利特球形度檢驗(yàn)的顯著性概率是0，說(shuō)明變量間存在相關(guān)關(guān)系，可以進(jìn)行PCA。

PCA方法的目的是用最少的主成分來(lái)表示原始變量最多的信息，可以根據(jù)主成分的累積貢獻(xiàn)率和特征根來(lái)確定主成分個(gè)數(shù)。一般來(lái)說(shuō)，累積貢獻(xiàn)率達(dá)到80%以上就為滿意；特征根是各個(gè)主成分影響大小的參數(shù)，一般保留特征根大于1所對(duì)應(yīng)的主成分［22］。本課題研究中，7個(gè)自變量的主成分的特征根和累積貢獻(xiàn)率如表6所示。由表6可知，第一、第二和第三主成分的特征根都大于1，同時(shí)，這3個(gè)主成分的累積貢獻(xiàn)率為95.76%，所以，可選擇提取這3個(gè)變量進(jìn)行PCA。

表5 KMO檢驗(yàn)和巴特利特球形度檢驗(yàn)Table 5 KMO and Bartlett’s tests

表6 主成分的特征根和方差貢獻(xiàn)率Table 6 Eigenvalues and its total variance explained

主成分的因子載荷反映變量在主成分上的重要性，3個(gè)主成分的因子載荷如表7所示。由表7可知，纖維長(zhǎng)度、纖維寬度、細(xì)小纖維含量、緊度和透氣度在第一主成分占較高的載荷，孔隙率在第二主成分占較高的載荷，厚度在第三主成分占有較高的載荷。如此，原來(lái)的7項(xiàng)參數(shù)轉(zhuǎn)變?yōu)?個(gè)新變量，而這3個(gè)新變量能夠解釋原參數(shù)95.76%的信息。

表7 主成分的因子載荷Table 7 Component matrix of principal components

進(jìn)行多元線性回歸分析，建立標(biāo)準(zhǔn)化的因變量絕緣紙工頻擊穿強(qiáng)度（zy）與主成分F1、F2、F3的回歸方程，如式（8）所示。

此模型的擬合度R2為0.917，說(shuō)明模型具有較好的模擬性。為了進(jìn)一步說(shuō)明模型的合理性，圖3給出了PCA回歸模型標(biāo)準(zhǔn)化殘差的結(jié)果。由圖3可知，PCA回歸模型標(biāo)準(zhǔn)化殘差的正態(tài)P-P圖呈近似線性趨勢(shì)，表明殘差服從正態(tài)分布，說(shuō)明回歸方程有效。對(duì)上述所求得的回歸模型進(jìn)行顯著性檢驗(yàn)，即方差分析。方差分析是對(duì)模型中多個(gè)自變量與因變量之間的線性關(guān)系在總體上是不是顯著作出判斷。方差分析結(jié)果表明，F(xiàn)=29.434，P=0＜0.05，說(shuō)明通過(guò)檢驗(yàn)，多元回歸模型有意義。

圖3 PCA回歸模型標(biāo)準(zhǔn)化殘差的正態(tài)P-P圖Fig.3 The normal P-P plot of standardized residuals of the PCA regression model

通過(guò)轉(zhuǎn)化，最終可得各參數(shù)與絕緣紙工頻擊穿強(qiáng)度的多元回歸模型，如式（9）所示。

相對(duì)誤差是指誤差在真實(shí)結(jié)果中所占的比例，這對(duì)于比較在各種情況下所得結(jié)果的準(zhǔn)確性更為方便。PCA回歸模型預(yù)測(cè)值與實(shí)測(cè)值（見(jiàn)表3）的對(duì)比結(jié)果如圖4所示。由圖4可知，除了樣本3號(hào)和7號(hào)的相對(duì)偏差為13.65%和11.84%，其余樣本的相對(duì)偏差均在10%以內(nèi)，說(shuō)明該模型具有較好的模擬性。

為了對(duì)所得回歸模型進(jìn)行驗(yàn)證，需要制備驗(yàn)證樣本以檢驗(yàn)所構(gòu)建模型的準(zhǔn)確性。驗(yàn)證樣本的各項(xiàng)參數(shù)如表8所示。將2組驗(yàn)證樣本（13和14）數(shù)據(jù)代入上述回歸模型中可得，兩樣本的工頻擊穿強(qiáng)度的預(yù)測(cè)值分別為44.83和30.65 kV/mm，與實(shí)測(cè)值（44.91和32.31 kV/mm）對(duì)比，計(jì)算所得偏差分別為?0.18%和?5.16%。可見(jiàn)，本課題建立的回歸模型合理可信，具有良好的預(yù)測(cè)精度。

圖4 PCA回歸模型預(yù)測(cè)值與實(shí)測(cè)值的對(duì)比Fig.4 Contrast between fitting values of PCA regression model and actual values

表8 驗(yàn)證樣本參數(shù)Table 8 Parameters of verification samples

2.3.3 最佳子集選擇法的多元線性回歸模型

由于本課題所選擇變量的期望數(shù)量為2～4，因此在進(jìn)行最佳子集選擇時(shí)，給定變量的范圍為1～6。絕緣紙工頻擊穿強(qiáng)度的最佳子集選擇結(jié)果如表9所示。由表9可知，通過(guò)比較不同變量數(shù)量的最佳子集可以看出，隨著變量數(shù)量的增加，擬合度R2增大，說(shuō)明回歸模型的擬合效果得到改善；但是，當(dāng)自變量數(shù)量達(dá)到一定值時(shí)，R2的增大非常有限。變量數(shù)量從1增加到3，R2增大了0.134；從3增加到4時(shí)，R2只增大了0.012。因此，當(dāng)變量數(shù)量達(dá)到3時(shí)，繼續(xù)增加變量數(shù)量并不會(huì)使R2值顯著增大。

Mallows＇Cp描述了模型的偏差和準(zhǔn)確性，可以幫助在多個(gè)回歸模型中進(jìn)行選擇，只有當(dāng)Mallows＇Cp值接近預(yù)測(cè)變量數(shù)量加上常量，模型預(yù)測(cè)精度才更準(zhǔn)確［4］。由表9可知，當(dāng)變量數(shù)量為3時(shí)，Mallows＇Cp值為6.1，接近于6，符合要求。綜合考慮擬合度R2值和Mallows＇Cp值可知，適合建立絕緣紙工頻擊穿強(qiáng)度的最佳子集選擇多元線性回歸模型（以下簡(jiǎn)稱最佳子集選擇回歸模型）的參數(shù)是3個(gè)變量，分別為細(xì)小纖維含量、厚度和孔隙率。

表9 最佳子集選擇結(jié)果Table 9 Results for the best subset selection

對(duì)絕緣紙工頻擊穿強(qiáng)度與細(xì)小纖維含量、厚度和孔隙率進(jìn)行多元線性回歸，得到模型如式（10）所示。

表10是最佳子集選擇回歸模型的基本信息。由表10可知，回歸模型擬合度R2達(dá)到0.943，說(shuō)明模型具有較好的模擬性。回歸系數(shù)顯著性檢驗(yàn)可反映自變量與因變量的線性關(guān)系是否顯著，x3、x4和x6的P值均小于0.05，說(shuō)明它們與工頻擊穿強(qiáng)度存在顯著的線性關(guān)系；x3、x4和x6的VIF非常接近于1，表明變量之間不存在多重共線性問(wèn)題，不會(huì)對(duì)回歸模型產(chǎn)生嚴(yán)重影響。

表10 最佳子集選擇回歸模型Table 10 Regression model for the best subset selection

圖5給出了最佳子集選擇回歸模型標(biāo)準(zhǔn)化殘差的結(jié)果。由圖5可知，最佳子集回歸模型標(biāo)準(zhǔn)化殘差的正態(tài)P-P圖近似呈線性趨勢(shì)，表明殘差服從正態(tài)分布，說(shuō)明建立的回歸方程有效。方差分析結(jié)果表明，F(xiàn)=43.883，P=0＜0.05，說(shuō)明通過(guò)檢驗(yàn)，回歸方程有意義。

最佳子集選擇回歸模型預(yù)測(cè)值與實(shí)測(cè)值（見(jiàn)表3）的對(duì)比結(jié)果如圖6所示。由圖6可知，預(yù)測(cè)值和實(shí)測(cè)值的最大偏差為9.55%，二者非常接近，說(shuō)明該模型具有較好的模擬性。將兩組驗(yàn)證樣本數(shù)據(jù)代入模型中，其預(yù)測(cè)值分別為44.62和36.57 kV/mm，與實(shí)測(cè)值對(duì)比，計(jì)算所得偏差分別為?0.64%和13.19%。可見(jiàn)，本課題建立的最佳子集選擇回歸模型合理可信，具有良好的預(yù)測(cè)精度。

圖5 最佳子集選擇回歸模型標(biāo)準(zhǔn)化殘差的正態(tài)P-P圖Fig.5 The normal P-P plot of standardized residuals of the best subset selection method regression model

圖6 實(shí)測(cè)值與最佳子集選擇回歸模型預(yù)測(cè)值對(duì)比Fig.6 Contrast between actual values and fitting values of regression model based on best subset selection method

綜上可知，PCA解決了建模過(guò)程中多重共線性的問(wèn)題，既不丟失原始數(shù)據(jù)的主要信息，又容易抓住主要矛盾，若將所有變量納入模型，模型變得較為復(fù)雜。最佳子集選擇方法簡(jiǎn)單直觀，所得模型有較高的擬合度，降低了模型的復(fù)雜程度。

3 結(jié)論

本課題提出一種基于灰色關(guān)聯(lián)分析的絕緣紙工頻擊穿強(qiáng)度影響因素的量化分析模型，分別采用主成分分析和最佳子集選擇方法進(jìn)行多參數(shù)優(yōu)化，構(gòu)建絕緣紙工頻擊穿強(qiáng)度的多元線性回歸模型，主要結(jié)論如下。

3.1 絕緣紙工頻擊穿強(qiáng)度影響因素的灰色關(guān)聯(lián)排序?yàn)椋豪w維長(zhǎng)度＞細(xì)小纖維含量＞纖維寬度＞緊度＞透氣度＞孔隙率＞厚度。

3.2 由相關(guān)性分析得出，纖維寬度、細(xì)小纖維含量和緊度與絕緣紙工頻擊穿強(qiáng)度呈正相關(guān)；纖維長(zhǎng)度、孔隙率和透氣度與絕緣紙工頻擊穿強(qiáng)度呈負(fù)相關(guān)。

3.3 利用主成分分析和最佳子集選擇方法構(gòu)建絕緣紙工頻擊穿強(qiáng)度的多元線型回歸模型，模型的R2分別為0.917和0.943，兩種模型的擬合度均較高；兩種模型中，模擬樣本和驗(yàn)證樣本預(yù)測(cè)結(jié)果的相對(duì)偏差基本在10%以內(nèi)，說(shuō)明兩種模型均具有良好的預(yù)測(cè)能力。