本發(fā)明屬于電氣工程技術(shù)領(lǐng)域，尤其涉及一種ENSO時(shí)期覆冰響應(yīng)程度計(jì)算和分析方法。

背景技術(shù)：

隨著能源互聯(lián)網(wǎng)的推進(jìn)和跨區(qū)電網(wǎng)的建設(shè)，輸電線路走廊受到的諸多環(huán)境因素影響越來(lái)越受到人們的重視。覆冰嚴(yán)重威脅輸電線路安全穩(wěn)定運(yùn)行。2008年的南方地區(qū)冰災(zāi)造成了輸電線路倒塔超過(guò)1300基，財(cái)產(chǎn)損失超過(guò)120億元，京廣電氣化鐵路停運(yùn)，對(duì)社會(huì)穩(wěn)定和人民生產(chǎn)與生活造成嚴(yán)重威脅。2015年遼寧“11.7”冰災(zāi)造成38萬(wàn)用戶(hù)停電，跳閘超過(guò)300條次，累計(jì)電力損失64萬(wàn)千瓦時(shí)，暴露出北方地區(qū)輸電線路在抵御覆冰災(zāi)害方面存在不足。因此，開(kāi)展冬季覆冰預(yù)測(cè)和覆冰特征分析對(duì)于合理調(diào)配資源、降低覆冰災(zāi)害損失的必要性不言而喻。而ENSO事件(包括厄爾尼諾事件和拉尼娜事件)通過(guò)影響我國(guó)氣候特征，在一定程度上控制了冬季覆冰發(fā)生條件。目前的電網(wǎng)覆冰程度預(yù)測(cè)方法雖然已經(jīng)將ENSO影響考慮在內(nèi)，但主要針對(duì)一段時(shí)間的數(shù)據(jù)進(jìn)行無(wú)差別統(tǒng)計(jì)，忽視了ENSO時(shí)期的特殊性，因此在分析ENSO時(shí)期覆冰響應(yīng)的一般性特征時(shí)仍然存在針對(duì)性不夠明顯的情況。為總結(jié)覆冰發(fā)生和增長(zhǎng)規(guī)律，進(jìn)行ENSO時(shí)期覆冰響應(yīng)程度的研究和分析已經(jīng)勢(shì)在必行，以期開(kāi)展針對(duì)性的覆冰預(yù)測(cè)預(yù)警，盡早制定電網(wǎng)防災(zāi)措施，減少電網(wǎng)災(zāi)害可能帶來(lái)的損失，維護(hù)電網(wǎng)安全穩(wěn)定運(yùn)行。

技術(shù)實(shí)現(xiàn)要素：

針對(duì)ENSO時(shí)期輸電線路覆冰一般性特征及其對(duì)ENSO事件的響應(yīng)情況缺乏分析的現(xiàn)狀，本發(fā)明提供一種ENSO時(shí)期覆冰響應(yīng)程度分析方法，該方法思路新穎、流程清晰、操作簡(jiǎn)單。

為實(shí)現(xiàn)上述目的，本發(fā)明的技術(shù)方案如下：

一種ENSO時(shí)期覆冰響應(yīng)程度分析方法，包括如下步驟：

(1)、獲取待分析地區(qū)過(guò)去一段時(shí)間對(duì)應(yīng)的覆冰數(shù)據(jù)和覆冰期每個(gè)月的ENSO指數(shù)數(shù)據(jù)，覆冰數(shù)據(jù)包括覆冰期每個(gè)月待分析地區(qū)的最大覆冰日數(shù)和最大覆冰厚度；

(2.1)、依據(jù)ENSO事件的性質(zhì)，將ENSO指數(shù)數(shù)據(jù)按照厄爾尼諾時(shí)期數(shù)據(jù)、拉尼娜時(shí)期數(shù)據(jù)、非ENSO時(shí)期數(shù)據(jù)劃分為三類(lèi)，每一類(lèi)中的數(shù)據(jù)仍然按照時(shí)間順序排列；

(2.2)、對(duì)照ENSO指數(shù)所表征的厄爾尼諾時(shí)期、拉尼娜時(shí)期、非ENSO時(shí)期，將待分析地區(qū)覆冰數(shù)據(jù)分別劃分為厄爾尼諾時(shí)期數(shù)據(jù)、拉尼娜時(shí)期數(shù)據(jù)、非ENSO時(shí)期數(shù)據(jù)三類(lèi)，每一類(lèi)中的數(shù)據(jù)仍然按照時(shí)間順序排列；

(3.1)、選取步驟(2.2)中得到的待分析地區(qū)厄爾尼諾時(shí)期覆冰數(shù)據(jù)，按照公式(1)，計(jì)算得到其厄爾尼諾時(shí)期最大覆冰日數(shù)的平均值和最大覆冰厚度的平均值再選取待分析地區(qū)拉尼娜時(shí)期覆冰數(shù)據(jù)，按照公式(1)，計(jì)算得到其拉尼娜時(shí)期最大覆冰日數(shù)的平均值和最大覆冰厚度的平均值選取待分析地區(qū)非ENSO時(shí)期覆冰數(shù)據(jù)，按照公式(1)，計(jì)算得到其非ENSO時(shí)期最大覆冰日數(shù)的平均值和最大覆冰厚度的平均值

$\stackrel{\‾}{x}=\frac{\underset{i=1}{\overset{n}{\Σ}}{x}_i}{n}---\left(1\right)$

式中，為非ENSO時(shí)期、厄爾尼諾時(shí)期或拉尼娜時(shí)期中某一時(shí)期最大覆冰日數(shù)或最大覆冰厚度的平均值，x_i為這一時(shí)期最大覆冰日數(shù)或最大覆冰厚度的數(shù)據(jù)，n為該時(shí)期數(shù)據(jù)總量；

(3.2)、以ENSO指數(shù)數(shù)據(jù)代替覆冰數(shù)據(jù)，重復(fù)步驟(3.1)，得到不同時(shí)期ENSO指數(shù)平均值；

(4.1)、將步驟(3.1)中得到的待分析地區(qū)厄爾尼諾時(shí)期最大覆冰日數(shù)平均值與待分析地區(qū)非ENSO時(shí)期最大覆冰日數(shù)平均值相減，得到的結(jié)果作為待分析地區(qū)厄爾尼諾時(shí)期最大覆冰日數(shù)距平值

(4.2)、將步驟(3.1)中得到的待分析地區(qū)厄爾尼諾時(shí)期最大覆冰厚度平均值與待分析地區(qū)非ENSO時(shí)期最大覆冰厚度平均值相減，得到的結(jié)果作為待分析地區(qū)厄爾尼諾時(shí)期最大覆冰厚度距平值

(4.3)、分別以待分析地區(qū)拉尼娜時(shí)期最大覆冰日數(shù)平均值和最大覆冰厚度平均值代替和重復(fù)步驟(4.1)～(4.2)，得到待分析地區(qū)拉尼娜時(shí)期最大覆冰日數(shù)距平值A(chǔ)_a2和最大覆冰厚度距平值A(chǔ)_b2；

(5.1)、選取步驟(2.1)～(4.3)中得到的厄爾尼諾時(shí)期ENSO指數(shù)數(shù)據(jù)與待分析地區(qū)厄爾尼諾時(shí)期最大覆冰日數(shù)，按照公式(2)，計(jì)算得到待分析地區(qū)最大覆冰日數(shù)與厄爾尼諾強(qiáng)度的相關(guān)系數(shù)r_a1：

$r=\frac{\underset{i=1}{\overset{n}{\Σ}}(x_i-\stackrel{\‾}{x})(y_i-\stackrel{\‾}{y})}{\sqrt{\underset{i=1}{\overset{n}{\Σ}}{(x_i-\stackrel{\‾}{x})}^2\underset{i=1}{\overset{n}{\Σ}}{(y_i-\stackrel{\‾}{y})}^2}}---\left(2\right)$

式中，r為待分析地區(qū)最大覆冰日數(shù)或最大覆冰厚度與厄爾尼諾強(qiáng)度或拉尼娜強(qiáng)度的相關(guān)系數(shù)，x_i為厄爾尼諾時(shí)期或拉尼娜時(shí)期中某一時(shí)期的最大覆冰日數(shù)或最大覆冰厚度，為這一時(shí)期該種覆冰數(shù)據(jù)的平均值，y_i為這一時(shí)期ENSO指數(shù)的數(shù)據(jù)，為這一時(shí)期ENSO指數(shù)的平均值，n為這一時(shí)期數(shù)據(jù)總量；

(5.2)、以待分析地區(qū)厄爾尼諾時(shí)期最大覆冰厚度數(shù)據(jù)和最大覆冰厚度平均值分別代替厄爾尼諾時(shí)期最大覆冰日數(shù)數(shù)據(jù)和最大覆冰日數(shù)平均值，代入公式(2)，計(jì)算得到待分析地區(qū)最大覆冰厚度與厄爾尼諾強(qiáng)度的相關(guān)系數(shù)r_b1；

(5.3)、分別以待分析地區(qū)拉尼娜時(shí)期最大覆冰日數(shù)和覆冰厚度代替厄爾尼諾時(shí)期最大覆冰日數(shù)和覆冰厚度，重復(fù)步驟(5.1)～(5.2)，得到待分析地區(qū)最大覆冰日數(shù)與拉尼娜強(qiáng)度的相關(guān)系數(shù)r_a2和最大覆冰厚度與拉尼娜強(qiáng)度的相關(guān)系數(shù)r_b2；

(6.1)、將步驟(4.1)～(4.2)計(jì)算得到的厄爾尼諾時(shí)期距平值進(jìn)行分析，若待分析地區(qū)厄爾尼諾時(shí)期最大覆冰日數(shù)距平值A(chǔ)_a1或最大覆冰厚度距平值A(chǔ)_b1為正，則認(rèn)為厄爾尼諾容易造成最大覆冰日數(shù)或最大覆冰厚度的增加；否則認(rèn)為厄爾尼諾容易造成最大覆冰日數(shù)或最大覆冰厚度的減少；

再將上述A_a1和A_b1代入公式(3)，計(jì)算得到厄爾尼諾影響系數(shù)：

$E=5\×\left(\right|\frac{A_a}{10}|+|\frac{A_b}{30}\left|\right)---\left(3\right)$

式中，E為ENSO影響系數(shù)，其代入厄爾尼諾時(shí)期數(shù)據(jù)后得到厄爾尼諾影響系數(shù)，代入拉尼娜時(shí)期數(shù)據(jù)后得到拉尼娜影響系數(shù)，A_a為對(duì)應(yīng)厄爾尼諾時(shí)期或拉尼娜時(shí)期最大覆冰日數(shù)距平值，A_b為這一時(shí)期最大覆冰厚度距平值；

基于厄爾尼諾影響系數(shù)的覆冰距平響應(yīng)程度判斷方式如表(1)所示：

表(1)覆冰距平響應(yīng)程度表

(6.2)、對(duì)步驟(5.1)～(5.2)得到的待分析地區(qū)覆冰與厄爾尼諾強(qiáng)度相關(guān)系數(shù)進(jìn)行顯著性檢驗(yàn)，查相關(guān)系數(shù)顯著性檢驗(yàn)表，若待分析地區(qū)最大覆冰日數(shù)與厄爾尼諾強(qiáng)度的相關(guān)系數(shù)r_a1的絕對(duì)值或最大覆冰厚度與厄爾尼諾強(qiáng)度的相關(guān)系數(shù)r_b1的絕對(duì)值大于表中對(duì)應(yīng)自由度下的顯著相關(guān)閾值r₀，則認(rèn)為該相關(guān)系數(shù)通過(guò)顯著性檢驗(yàn)，厄爾尼諾事件對(duì)待分析地區(qū)覆冰持續(xù)時(shí)間或覆冰厚度的影響是顯著的；否則認(rèn)為該相關(guān)系數(shù)沒(méi)有通過(guò)顯著性檢驗(yàn)，厄爾尼諾事件對(duì)待分析地區(qū)覆冰持續(xù)時(shí)間或覆冰厚度的影響可能存在，但并不顯著的；

(6.3)、分別以拉尼娜時(shí)期距平值和待分析地區(qū)最大覆冰日數(shù)/最大覆冰厚度與拉尼娜強(qiáng)度相關(guān)系數(shù)分別代替厄爾尼諾時(shí)期距平值和待分析地區(qū)最大覆冰日數(shù)/最大覆冰厚度與厄爾尼諾強(qiáng)度相關(guān)系數(shù)，重復(fù)步驟(6.1)～(6.2)，進(jìn)行拉尼娜時(shí)期距平值和相關(guān)性的響應(yīng)分析，步驟(6.2)中以步驟(5.3)代替步驟(5.1)～(5.2)。

本發(fā)明的有益效果是：

1、可以較為全面的認(rèn)識(shí)ENSO時(shí)期輸電線路覆冰一般性特征及其對(duì)ENSO事件的響應(yīng)情況；

2、可操作性強(qiáng)；

3、提高了輸電線路覆冰防治工作的針對(duì)性。根據(jù)分析結(jié)果，可及時(shí)做好相應(yīng)的應(yīng)急處置措施，對(duì)于覆冰對(duì)ENSO事件響應(yīng)比較明顯的地區(qū)，在ENSO時(shí)期的冬季做好提前應(yīng)對(duì)工作，減少電網(wǎng)損失。

附圖說(shuō)明

圖1為本發(fā)明實(shí)施例分析方法的流程圖。

具體實(shí)施方式

下面結(jié)合附圖及實(shí)例，對(duì)本發(fā)明做進(jìn)一步說(shuō)明。

本發(fā)明實(shí)施例的思路是：一種ENSO時(shí)期覆冰響應(yīng)程度分析方法，系考慮各省輸電線路覆冰在ENSO時(shí)期的一般性特征以及ENSO事件強(qiáng)度對(duì)輸電線路覆冰的影響，建立基于ENSO時(shí)期覆冰距平值和覆冰與ENSO強(qiáng)度相關(guān)性的計(jì)算方法，得到ENSO時(shí)期覆冰響應(yīng)程序的分析結(jié)論。

如圖1所示，本發(fā)明實(shí)施例分析方法包括下述步驟。

(1)、數(shù)據(jù)獲取。獲取時(shí)間對(duì)應(yīng)的覆冰數(shù)據(jù)和ENSO指數(shù)數(shù)據(jù)，獲取各省過(guò)去20年的覆冰數(shù)據(jù)，包括覆冰期每個(gè)月各省的最大覆冰日數(shù)和最大覆冰厚度。獲取過(guò)去20年覆冰期每個(gè)月的ENSO指數(shù)數(shù)據(jù)；

(2.1)、ENSO指數(shù)數(shù)據(jù)劃分。依據(jù)ENSO事件的性質(zhì)(ENSO事件包括厄爾尼諾事件和拉尼娜事件)，將ENSO指數(shù)數(shù)據(jù)按照厄爾尼諾時(shí)期數(shù)據(jù)、拉尼娜時(shí)期數(shù)據(jù)、非ENSO時(shí)期數(shù)據(jù)劃分為三類(lèi)，每一類(lèi)中的數(shù)據(jù)仍然按照時(shí)間順序排列；

(2.2)、覆冰數(shù)據(jù)劃分。對(duì)照ENSO指數(shù)所表征的厄爾尼諾時(shí)期、拉尼娜時(shí)期、非ENSO時(shí)期，將各省覆冰數(shù)據(jù)分別劃分為厄爾尼諾時(shí)期數(shù)據(jù)、拉尼娜時(shí)期數(shù)據(jù)、非ENSO時(shí)期數(shù)據(jù)三類(lèi)，每一類(lèi)中的數(shù)據(jù)仍然按照時(shí)間順序排列；

(3.1)、平均值計(jì)算。選取步驟(2.2)中得到的某省厄爾尼諾時(shí)期覆冰數(shù)據(jù)，按照公式(1)，計(jì)算得到其厄爾尼諾時(shí)期最大覆冰日數(shù)的平均值和最大覆冰厚度的平均值再選取該省拉尼娜時(shí)期覆冰數(shù)據(jù)，按照公式(1)，計(jì)算得到其拉尼娜時(shí)期最大覆冰日數(shù)的平均值和最大覆冰厚度的平均值選取該省非ENSO時(shí)期覆冰數(shù)據(jù)，按照公式(1)，計(jì)算得到其非ENSO時(shí)期最大覆冰日數(shù)的平均值和最大覆冰厚度的平均值

$\stackrel{\‾}{x}=\frac{\underset{i=1}{\overset{n}{\Σ}}{x}_i}{n}---\left(1\right)$

式中，為某一時(shí)期(非ENSO時(shí)期、厄爾尼諾時(shí)期或拉尼娜時(shí)期)最大覆冰日數(shù)或最大覆冰厚度的平均值，x_i為這一時(shí)期(非ENSO時(shí)期、厄爾尼諾時(shí)期或拉尼娜時(shí)期)最大覆冰日數(shù)或最大覆冰厚度的數(shù)據(jù)，n為該時(shí)期數(shù)據(jù)總量。

(3.2)、以ENSO指數(shù)數(shù)據(jù)代替覆冰數(shù)據(jù)，重復(fù)步驟(3.1)，得到不同時(shí)期ENSO指數(shù)平均值；

(4.1)、最大覆冰日數(shù)距平值計(jì)算。將步驟(3)中得到的某省厄爾尼諾時(shí)期最大覆冰日數(shù)平均值與對(duì)應(yīng)省份非ENSO時(shí)期最大覆冰日數(shù)平均值相減得到的結(jié)果作為該省厄爾尼諾時(shí)期最大覆冰日數(shù)距平值A(chǔ)_a1；

(4.2)、最大覆冰厚度距平值計(jì)算。將步驟(3)中得到的某省厄爾尼諾時(shí)期最大覆冰厚度平均值與對(duì)應(yīng)省份非ENSO時(shí)期最大覆冰厚度平均值相減得到的結(jié)果作為該省厄爾尼諾時(shí)期最大覆冰厚度距平值A(chǔ)_b1。

(4.3)、分別以所選省份拉尼娜時(shí)期最大覆冰日數(shù)平均值和最大覆冰厚度平均值代替和重復(fù)步驟(4.1)～(4.2)，得到該省拉尼娜時(shí)期最大覆冰日數(shù)距平值A(chǔ)_a2和最大覆冰厚度距平值A(chǔ)_b2；

(5.1)、最大覆冰日數(shù)與ENSO強(qiáng)度相關(guān)性計(jì)算。選取步驟(2.1)～(4.3)中得到的厄爾尼諾時(shí)期ENSO指數(shù)數(shù)據(jù)與某省厄爾尼諾時(shí)期最大覆冰日數(shù)，按照公式(2)，計(jì)算得到該省最大覆冰日數(shù)與厄爾尼諾強(qiáng)度的相關(guān)系數(shù)r_a1：

$r=\frac{\underset{i=1}{\overset{n}{\Σ}}(x_i-\stackrel{\‾}{x})(y_i-\stackrel{\‾}{y})}{\sqrt{\underset{i=1}{\overset{n}{\Σ}}{(x_i-\stackrel{\‾}{x})}^2\underset{i=1}{\overset{n}{\Σ}}{(y_i-\stackrel{\‾}{y})}^2}}---\left(2\right)$

式中，r為待分析地區(qū)最大覆冰日數(shù)或最大覆冰厚度與厄爾尼諾強(qiáng)度或拉尼娜強(qiáng)度的相關(guān)系數(shù)，x_i為某一時(shí)期(厄爾尼諾時(shí)期或拉尼娜時(shí)期)的某種覆冰數(shù)據(jù)(最大覆冰日數(shù)或最大覆冰厚度)，為這一時(shí)期該種覆冰數(shù)據(jù)的平均值，y_i為這一時(shí)期ENSO指數(shù)的數(shù)據(jù)，為這一時(shí)期ENSO指數(shù)的平均值，n為這一時(shí)期數(shù)據(jù)總量；

(5.2)、最大覆冰厚度與ENSO強(qiáng)度相關(guān)性計(jì)算。以某省厄爾尼諾時(shí)期最大覆冰厚度數(shù)據(jù)和最大覆冰厚度的平均值分別代替厄爾尼諾時(shí)期最大覆冰日數(shù)數(shù)據(jù)和最大覆冰日數(shù)平均值，代入公式(2)，計(jì)算得到該省最大覆冰厚度與厄爾尼諾強(qiáng)度的相關(guān)系數(shù)r_b1。

(5.3)、分別以所選省份拉尼娜時(shí)期最大覆冰日數(shù)和覆冰厚度代替厄爾尼諾時(shí)期最大覆冰日數(shù)和覆冰厚度，重復(fù)步驟(5.1)～(5.2)，得到該省最大覆冰日數(shù)與拉尼娜強(qiáng)度的相關(guān)系數(shù)r_a2和最大覆冰厚度與拉尼娜強(qiáng)度的相關(guān)系數(shù)r_b2；

(6.1)、ENSO時(shí)期距平值響應(yīng)情況。將步驟(4.1)～(4.2)計(jì)算得到的厄爾尼諾時(shí)期距平值進(jìn)行分析，若某省厄爾尼諾時(shí)期最大覆冰日數(shù)距平值A(chǔ)_a1或最大覆冰厚度距平值A(chǔ)_b1為正，則認(rèn)為厄爾尼諾容易造成最大覆冰日數(shù)或最大覆冰厚度的增加。否則認(rèn)為厄爾尼諾容易造成最大覆冰日數(shù)或最大覆冰厚度的減少。

再將上述A_a1和A_b1代入公式(3)，計(jì)算得到厄爾尼諾影響系數(shù)：

$E=5\×\left(\right|\frac{A_a}{10}|+|\frac{A_b}{30}\left|\right)---\left(3\right)$

式中，E為ENSO影響系數(shù)(代入厄爾尼諾時(shí)期數(shù)據(jù)后得到厄爾尼諾影響系數(shù)，代入拉尼娜時(shí)期數(shù)據(jù)后得到拉尼娜影響系數(shù))，A_a為對(duì)應(yīng)時(shí)期(厄爾尼諾時(shí)期或拉尼娜時(shí)期)最大覆冰日數(shù)距平值，A_b為這一時(shí)期最大覆冰厚度距平值。

基于厄爾尼諾影響系數(shù)的覆冰距平響應(yīng)程度判斷方式如表(1)所示：

表(1)覆冰距平響應(yīng)程度表

(6.2)、ENSO時(shí)期相關(guān)性響應(yīng)情況。對(duì)步驟(5.1)～(5.2)得到的各省覆冰與厄爾尼諾強(qiáng)度相關(guān)系數(shù)進(jìn)行顯著性檢驗(yàn)，查相關(guān)系數(shù)顯著性檢驗(yàn)表，若某省最大覆冰日數(shù)與厄爾尼諾強(qiáng)度的相關(guān)系數(shù)r_a1的絕對(duì)值或最大覆冰厚度與厄爾尼諾強(qiáng)度的相關(guān)系數(shù)r_b1的絕對(duì)值大于表中對(duì)應(yīng)自由度下的顯著相關(guān)閾值r₀，則認(rèn)為該相關(guān)系數(shù)通過(guò)顯著性檢驗(yàn)，厄爾尼諾事件對(duì)該省覆冰持續(xù)時(shí)間或覆冰厚度的影響是顯著的。否則認(rèn)為該相關(guān)系數(shù)沒(méi)有通過(guò)顯著性檢驗(yàn)，厄爾尼諾事件對(duì)該省覆冰持續(xù)時(shí)間或覆冰厚度的影響可能存在，但并不顯著的；

(6.3)分別以拉尼娜時(shí)期距平值和各省最大覆冰日數(shù)(最大覆冰厚度)與拉尼娜強(qiáng)度相關(guān)系數(shù)分別代替厄爾尼諾時(shí)期距平值和各省最大覆冰日數(shù)(最大覆冰厚度)與厄爾尼諾強(qiáng)度相關(guān)系數(shù)，重復(fù)步驟(6.1)～(6.2)，進(jìn)行拉尼娜時(shí)期距平值和相關(guān)性的響應(yīng)分析；此時(shí)步驟(6.2)中需以步驟(5.3)代替步驟(5.1)～(5.2)。

還可重復(fù)步驟(3.1)～(6.2)，完成所有省份的響應(yīng)分析。

本發(fā)明實(shí)施例的有益效果是：

1、可以較為全面的認(rèn)識(shí)ENSO時(shí)期輸電線路覆冰一般性特征及其對(duì)ENSO事件的響應(yīng)情況；

2、可操作性強(qiáng)；

3、提高了輸電線路覆冰防治工作的針對(duì)性。根據(jù)分析結(jié)果，可及時(shí)做好相應(yīng)的應(yīng)急處置措施，對(duì)于覆冰對(duì)ENSO事件響應(yīng)比較明顯的地區(qū)，在ENSO時(shí)期的冬季做好提前應(yīng)對(duì)工作，減少電網(wǎng)損失。

下面以遼寧省為例對(duì)本發(fā)明方法進(jìn)行說(shuō)明。該方法包括如下步驟：

(1)、數(shù)據(jù)獲取。獲取時(shí)間對(duì)應(yīng)的覆冰數(shù)據(jù)和ENSO指數(shù)數(shù)據(jù)：獲取遼寧省過(guò)去20年的覆冰數(shù)據(jù)，包括覆冰期每個(gè)月的最大覆冰日數(shù)和最大覆冰厚度；獲取過(guò)去20年覆冰期每個(gè)月的ONI(海洋NINO指數(shù)，一種ENSO指數(shù))數(shù)據(jù)；

(2.1)、ENSO指數(shù)數(shù)據(jù)劃分。依據(jù)ENSO事件的性質(zhì)(ENSO事件包括厄爾尼諾事件和拉尼娜事件)，將ONI數(shù)據(jù)按照厄爾尼諾時(shí)期數(shù)據(jù)、拉尼娜時(shí)期數(shù)據(jù)、非ENSO時(shí)期數(shù)據(jù)劃分為三類(lèi)，每一類(lèi)中的數(shù)據(jù)仍然按照時(shí)間順序排列；

(2.2)、覆冰數(shù)據(jù)劃分。選取遼寧省的覆冰數(shù)據(jù)，對(duì)照ONI所表征的厄爾尼諾時(shí)期、拉尼娜時(shí)期、非ENSO時(shí)期，覆冰數(shù)據(jù)分別劃分為厄爾尼諾時(shí)期數(shù)據(jù)、拉尼娜時(shí)期數(shù)據(jù)、非ENSO時(shí)期數(shù)據(jù)三類(lèi)，每一類(lèi)中的數(shù)據(jù)仍然按照時(shí)間順序排列；

(3.1)、平均值計(jì)算。基于步驟(2.2)中得到的數(shù)據(jù)，選取非ENSO時(shí)期覆冰數(shù)據(jù)，按照公式(1)，計(jì)算得到其非ENSO時(shí)期最大覆冰日數(shù)和最大覆冰厚度的平均值。選取厄爾尼諾時(shí)期覆冰數(shù)據(jù)，按照公式(1)，計(jì)算得到其厄爾尼諾時(shí)期最大覆冰日數(shù)和最大覆冰厚度的平均值。

(3.2)、再選取拉尼娜時(shí)期覆冰數(shù)據(jù)，按照公式(1)，計(jì)算得到其拉尼娜時(shí)期最大覆冰日數(shù)和最大覆冰厚度的平均值。以O(shè)NI數(shù)據(jù)代替覆冰數(shù)據(jù)，重復(fù)步驟(3.1)，得到不同時(shí)期的ONI平均值；

(4.1)、最大覆冰日數(shù)距平值計(jì)算。將步驟(3)中得到的厄爾尼諾時(shí)期最大覆冰日數(shù)平均值與對(duì)應(yīng)的非ENSO時(shí)期最大覆冰日數(shù)平均值相減，得到厄爾尼諾時(shí)期最大覆冰日數(shù)距平值1.3天；

(4.2)、最大覆冰厚度距平值計(jì)算。將步驟(3)中得到的厄爾尼諾時(shí)期最大覆冰厚度平均值與對(duì)應(yīng)的非ENSO時(shí)期最大覆冰厚度平均值相減，得到厄爾尼諾時(shí)期最大覆冰厚度距平值2.1毫米；

(4.3)、分別以拉尼娜時(shí)期最大覆冰日數(shù)和覆冰厚度代替厄爾尼諾時(shí)期最大覆冰日數(shù)和覆冰厚度，重復(fù)步驟(4.1)～(4.2)，得到拉尼娜時(shí)期最大覆冰日數(shù)距平值1.5天和最大覆冰厚度距平值1.1毫米；

(5.1)、最大覆冰日數(shù)與ENSO強(qiáng)度相關(guān)性計(jì)算。選取步驟(2.1)～(4.3)中得到的厄爾尼諾時(shí)期ONI數(shù)據(jù)與厄爾尼諾時(shí)期最大覆冰日數(shù)，按照公式(2)，計(jì)算得到最大覆冰日數(shù)與厄爾尼諾強(qiáng)度的相關(guān)系數(shù)為0.26；

(5.2)、最大覆冰厚度與ENSO強(qiáng)度相關(guān)性計(jì)算。選取步驟(2.1)～(4.2)中得到的厄爾尼諾時(shí)期ONI數(shù)據(jù)與厄爾尼諾時(shí)期最大覆冰厚度，按照公式(2)，計(jì)算得到最大覆冰厚度與厄爾尼諾強(qiáng)度的相關(guān)系數(shù)為0.19；

(5.3)、分別以拉尼娜時(shí)期最大覆冰日數(shù)和覆冰厚度代替厄爾尼諾時(shí)期最大覆冰日數(shù)和覆冰厚度，重復(fù)步驟(5.1)～(5.2)，得到最大覆冰日數(shù)與拉尼娜強(qiáng)度的相關(guān)系數(shù)0.21和最大覆冰厚度與拉尼娜強(qiáng)度的相關(guān)系數(shù)0.28；

(6.1)、ENSO時(shí)期距平值響應(yīng)情況。將步驟(4.1)～(4.2)計(jì)算得到的厄爾尼諾時(shí)期距平值進(jìn)行分析，由于遼寧省厄爾尼諾時(shí)期最大覆冰日數(shù)(最大覆冰厚度)距平值為正，認(rèn)為厄爾尼諾容易造成最大覆冰日數(shù)(最大覆冰厚度)的增加。通過(guò)公式(3)的計(jì)算得到厄爾尼諾影響系數(shù)為1.3，認(rèn)為覆冰距平對(duì)厄爾尼諾有輕微的響應(yīng)；

(6.2)、ENSO時(shí)期相關(guān)性響應(yīng)情況。對(duì)步驟(5.1)～(5.2)得到的遼寧省覆冰與厄爾尼諾強(qiáng)度相關(guān)系數(shù)進(jìn)行顯著性檢驗(yàn)，查相關(guān)系數(shù)顯著性檢驗(yàn)表，由于最大覆冰日數(shù)(最大覆冰厚度)與厄爾尼諾強(qiáng)度的相關(guān)系數(shù)絕對(duì)值小于表中對(duì)應(yīng)自由度下、95％顯著性水平下的顯著相關(guān)閾值0.29，則認(rèn)為厄爾尼諾事件對(duì)遼寧省覆冰持續(xù)時(shí)間(覆冰強(qiáng)度)的影響并不顯著；

(6.3)分別以拉尼娜時(shí)期距平值和遼寧省最大覆冰日數(shù)(最大覆冰厚度)與拉尼娜強(qiáng)度相關(guān)系數(shù)分別代替厄爾尼諾時(shí)期距平值和遼寧省最大覆冰日數(shù)(最大覆冰厚度)與厄爾尼諾強(qiáng)度相關(guān)系數(shù)，重復(fù)步驟(6.1)～(6.2)，進(jìn)行拉尼娜時(shí)期距平值和相關(guān)性響應(yīng)分析，此時(shí)步驟(6.2)中需以步驟(5.3)代替步驟(5.1)～(5.2)。同樣發(fā)現(xiàn)拉尼娜容易造成遼寧省最大覆冰日數(shù)(最大覆冰厚度)的增加，遼寧省覆冰對(duì)拉尼娜存在輕微響應(yīng)，同時(shí)拉尼娜事件對(duì)遼寧省覆冰持續(xù)時(shí)間(覆冰強(qiáng)度)的影響并不顯著。