ENSO综合指标及其径流预测方法与流程

本发明涉及一种ENSO综合指标及其径流预测方法，属于环保技术领域。

背景技术：

多元ENSO(MEI，Wolter and Timlin,1993,1998)指标是根据6种热带海洋观测要素经主成分分析得到的。这6种要素分别是海平面气压、纬向风、经向风、海表面温度、海面上空的气温、天空总云量。目前多元ENSO指标只是一个统计数据，它是海洋上6种气象要素的主成分分析结果。尽管它表达了每个月份海洋要素的综合状况，但是，却不能给出每次厄尔尼诺(拉尼娜)事件的真实评价，这与厄尔尼诺(拉尼娜)事件具有持续性、区域扩展性和强度特性等综合特征有密切关系。

技术实现要素：

目的：为了克服现有技术中存在的不足，本发明提供一种ENSO综合指标及其径流预测方法。

技术方案：为解决上述技术问题，本发明采用的技术方案为：

一种ENSO综合指标，构建ENSO综合指标，用于全面阐述一个ENSO事件的强度及等级：

首先，定义特征值用来描述每一个ENSO事件，用持续时间(MEI_month)反映厄尔尼诺/拉尼娜事件的持续性；用平均强度(MEI_mean)和极端强度(MEI_max)表示厄尔尼诺/拉尼娜事件的强度特性；用赤道太平洋海温增温≥0.5℃或≤-0.5℃范围代表空间扩展度(MEI_spa)；用温度变化达到或超过±2℃的区域代表极端空间扩展度(MEI_exspa)，利用这几个参数，构建ENSO综合指标(MEI_Z)；

以±0.5为临界值，提取每个连续5次达到±0.5的厄尔尼诺/拉尼娜事件过程，计算极端最大MEI值MEI_max、平均MEI值MEI_mean、达到标准后持续时间MEI_month，按月计算；提取赤道东太平洋海表温度≥0.5℃(或≤-0.5)的空间格点数，提取达到或超过±2℃的空间格点数，形成月文件，并将海表温度≥0.5℃(或≤-0.5)的空间格点数和达到或超过±2℃的空间格点数进行标准化，形成MEI_0.5和MEI_2.0序列；在达到厄尔尼诺/拉尼娜事件标准的时段内，计算MEI_spa和MEI_exspa；

将平均强度MEI_mean、极端强度MEI_max、持续时间MEI_month、空间扩展度MEI_spa和极端空间扩展度MEI_exspa按0.3：0.2：0.1：0.2：0.2的权重衡量一个ENSO事件，对历次ENSO事件进行综合指标计算，得到综合MEI值，即ENSO综合指标MEIz；

MEI_Z＝0.3MEI_mean+0.2MEI_max+0.1MEI_month+0.2MEI_spa+0.2MEI_exspa

计算拉尼娜La的综合值时，对MEI的数值先取绝对值，然后再进行权重计算的，因此，在MEIz划分的等级中，拉尼娜La与厄尔尼诺El等级标准一致，只是符号相反；MEIz的等级划分如下：

ENSO事件对流域流量的预测方法为：

有益效果：本发明提供的ENSO综合指标及其径流预测方法，具有以下优点：该方法提出了用持续时间(MEI_month)反映厄尔尼诺/拉尼娜事件的持续性；用平均强度(MEI_mean)和极端强度(MEI_max)表示厄尔尼诺/拉尼娜事件的强度特性；用赤道太平洋海温增温≥0.5℃或≤-0.5℃范围代表空间扩展度(MEI_spa)；用温度变化达到或超过±2℃的区域代表极端空间扩展度(MEI_exspa)，利用这五个参数，构建ENSO综合指标(MEI_Z)；

这是一种全新的构建方法，综合考虑了每个厄尔尼诺/拉尼娜事件的多种特性，是一种全新的描述厄尔尼诺/拉尼娜事件的方法。该方法全面考察了厄尔尼诺/拉尼娜事件对大气的影响，其作用不仅科学地辨别了每一次事件，而且客观地评价了厄尔尼诺/拉尼娜事件对大气的影响程度。这是目前任何一个指标都不具备的。

附图说明

图1是极强拉尼娜事件与红水河流域流量的关系图：(a)1954/1957事件，(b)1973/1976事件；

图2是极强厄尔尼诺事件与红水河流域流量的关系图：(a)1982/1983事件，(b)1991/1993事件，(c)1997/1998事件；

图3是极强拉尼娜事件与乌江流域流量的关系图：(a)1954/1957事件，(b)1973/1976事件；

图4是强拉尼娜事件与乌江流域流量的关系图：(a)1970/1972事件，(b)1998/2000事件，(c)2010/2011事件；

图5是强厄尔尼诺事件与乌江流域流量的关系图：(a)1957/1958事件，(b)1986/1988事件；

图6是极强拉尼娜事件与金沙江流域流量的关系图：(a)1954/1957事件，(b)1973/1976事件；

图7是强厄尔尼诺事件与金沙江流域流量的关系：(a)1957/1958事件，(b)1986/1988事件；

图8是极强拉尼娜事件与东江流域流量的关系：(a)1954/1957事件，(b)1973/1976事件；

图9是强拉尼娜事件与东江流域流量的关系：(a)1970/1972事件，(b)1998/2000事件，(c)2010/2011事件；

图10是极强厄尔尼诺事件与东江流域流量的关系：(a)1982/1983事件，(b)1991/1993事件，(c)1997/1998事件；

图11是强厄尔尼诺事件与东江流域流量的关系：(a)1957/1958事件，(b)1986/1988事件。

具体实施方式

下面结合具体实施例对本发明作更进一步的说明。

一种ENSO综合指标及其径流预测方法，对多元ENSO(MEI)进行了重新研究和划分。以±0.5为标准值，出现MEI≥0.5的正值频率大约为30.94％，出现MEI≤-0.5的负值频率为29.01％，其余40.05％为正常年的频率。

任何一个ENSO事件达到的强度、持续的月数差异很大，我们定义了五个特征值，包括极端值、平均值和持续月数来描述每一个事件。以±0.5为临界值，提取每个连续5次达到±0.5的El/La过程，计算极端最大MEI值(MEI_max)、平均MEI值(MEI_mean)、达到标准后持续的月数(MEI_month)、空间扩展度(MEI_spa)和极端异常区扩展度(MEI_exspa)。以不同的权重将这五个特征参数进行组合，可以得到综合MEI值(MEI_Z)。

ENSO事件对降水、气温和径流等现象的影响可通过事件极端强度、平均强度和持续时间等指标来体现。一个短暂(例如持续时间仅有半年)而很强(峰值很大)的ENSO事件和一个持续时间很长(例如持续时间仅有1年半以上)和一般强度(峰值较小)的ENSO事件对全球气候的影响是不同的。因此，一次ENSO(包括El/La)是极端强度、平均强度、持续过程、空间扩展度和极端异常区扩展度等多个指标的综合函数，仅用某个指标不能全面描述一个ENSO事件，因此，需要构建综合指标，全面阐述一个ENSO事件的强度及等级。

对于一个ENSO事件，极端强度、平均强度和持续时间可以较全面地描述它的特征，并使它可以被划分为不同等级，以此，探讨不同ENSO事件对气候的影响程度。对于一个ENSO事件，极端值虽然很重要，但是，由于极值持续时间很短，通常只有1-2个月，而平均强度可以很好地表征一个ENSO事件，因此，对于ENSO事件的综合描述，平均强度将是一个重要方面，而极端强度也应有一定代表性。于是，我们将平均强度、极端强度和持续时间按0.3：0.2：0.1：0.2：0.2的权重衡量一个ENSO事件，对历次ENSO事件进行综合指标计算，得到综合MEI值，即MEIz指标。表1是根据MEIz指标划分的ENSO事件的等级。

计算La的综合值时，对MEI的数值先取绝对值，然后再进行权重计算的，因此，在MEIz划分的等级中，La与El等级标准一致，只是符号相反。

表1 MEIz的等级划分

按照表1的标准，可以将1950-2015年期间的历次El/La事件挑选出来。表2和表3分别为最近60多年来El/La事件表。可以看到，用MEIz指标挑选出19次El事件(表2)和14次La事件(表3)。比较表2和表3可以发现，El事件的次数多于La事件，但是，每个El事件持续的时间都短于La在表2中，最长的El事件是1991/93，持续了32个月，其余事件的持续时间都小于20个月。1991/93的El事件在其他人的研究被分成了两个事件，一个是1991/92，另一个是1993事件，中间间隔了7个月；而MEI指标却一直持续正值，且都在El的标准之上，没有中断现象。在表3中，有两个持续时间最长的La事件，一个是1954/57，持续了34个月，该事件与其他人的结果是一致的，唯一不同是MEI指标给出的持续时间更长，其他人描述的该事件基本是20多个月，即持续期为2年多。另一个是1973/76，持续了35个月，该事件在其他人的研究中也被划分成两次事件，即1973/74和1974/76，两次事件中间隔了4个月，在MEI指标中，该次事件自1973年6月形成后，一直以较强的负值持续到1976年4月。

表2 El事件表(截止到2015年11月MEI指标统计)

表3 La事件表

综合MEI指标对流量的影响也显示出独特的效果。图1-图2是ENSO事件对红水河流域流量的影响。在极强拉尼娜事件的当年年底(11-12月，见图1)，MEI综合值达到2.0以上时，转年夏季，红水河流域的天一、光照、右江、龙滩和长洲流量将增大80％到150％。在转年秋季，这些站的流量将增大100％到200％。在极强厄尔尼诺事件夏季(特别是8月份，见图2)，MEI综合值达到1.8以上时，当年年底的流量将增大50％-100％，转年春季流量仍增大50％-200％。

在乌江流域，极强拉尼娜事件的强度在MEI_Z只有1.8左右时(见图3)，乌江流域转年夏季流量即可出现100-200％的增大现象。而强拉尼娜事件MEI_Z只有1.5左右时，乌江流域转年夏季流量即可出现100-150％的增大现象(见图4)。在强厄尔尼诺事件的夏季或当年年底，MEI_Z在1.3左右时(见图5)，乌江流域转年夏季流量即可出现100-200％的增大现象。

在金沙江流域，极强拉尼娜事件可引起转年夏季流量增大50％左右(见图6)，而强厄尔尼诺可引起转年流量减少40％左右(见图7)。

在东江流域，极强和强拉尼娜事件均可引起东江流域转年夏季流域减少50-80％(见图8和图9)，而极强厄尔尼诺可引起流量转年2-3月份增大200-700％，强厄尔尼诺事件可引起流量转年2-3月份增大50-100％(见图10和图11)。

对于ENSO事件的影响和预测，可归纳为表4的结果。综合而言，利用极强和强拉尼娜事件，可预测红水河、和乌江和金沙江流域流量增大，但是，这三个流域流量增大的幅度逐渐减小。利用拉尼娜事件，可预测东江的流量减小。利用极强和强厄尔尼诺事件，可预测红水河、乌江和东江流域流量增大，其中东江流域增大幅度最大，红水河次之，乌江最小。可预测金沙江流域的流量减小。

表4 ENSO事件对流域流量的预测

以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。