2. 兰州大学生命科学学院;
3. 兰州大学大气科学学院;
4. 陕西省榆林市气象局
2. School of Life Scie nce, Lanzhou University;
3. School of Atmospheric Sciences, Lanzhou University;
4. Yulin Meteorological Office, Shaanxi Province
近年来,全球气候持续变暖,各类极端天气事件日益频繁,天气气候异常变化对我国国民经济、社会发展和人民生产生活的影响日益增大。气象风险具有自然的和社会的双重属性,是一种自然-社会现象。由于天气气候预测中的不确定性是不可避免的,尤其作为气象风险源的高影响天气事件具有难以度量和不确定性强等特性,因此对其进行风险评估仍是一个难度较大的课题[1]。罗慧等人在相关前期研究的基础上,以西安市为例,挑选在2004年1月至2007年10月与公众生活关系密切的5类高影响天气事件作为气象风险源,即高温、雷暴、大雨、中雨、小雨等,结合这5类风险因子发生时所对应的当天12121的日拨打次数,利用基于模糊数学的风险评估方法计算得到了风险超越概率[2],那么如何进一步定量判断不同的气象风险源等级,以及在不同气象风险等级上的社会公众关注人数,是值得进一步研究的问题。本文试图通过多项式拟合方法,通过求出不同风险源的拐点,定量划分风险关注等级程度,进而计算对应等级的关注人数,这有利于政府决策者及社会公众在响应相关风险时有较为定量的判断。
1 风险评估的具体研究思路 1.1 基于模糊数学的超越概率计算方法高影响天气事件属于较小样本事件,应用建立在大数定理基础之上的传统概率统计方法给出的概率统计结果有时很不可靠,为了弥补信息不足而考虑优先利用样本的模糊信息,将单值样本变成集值样本,最成熟的模型是正态扩散模型[3-4]。结合西安市的实际情况,选择影响公众生活的天气事件作为产生气象相关风险的风险源,按照气象部门的定义挑选、选取了2004年1月至2007年10月的高温、雷暴、大到暴雨、中雨、小雨五种天气为研究对象。
1.2 整理气象风险源对应的社会公众关注样本量收集2004年1月至2007年10月逐日12121气象电话拨打次数,可以看出期间西安市12121气象咨询电话日拨打量最大峰值不超过9万次。通过整理期间的气象和12121拨打量数据,共得到108个高温样本量;37个雷暴样本量;16个大到暴雨样本量;38个中雨样本量;255个小雨样本量。
对应高温、雷暴、大雨、中雨和小雨五种高影响天气事件,根据实际拨打量的具体情况,从而确定连续论域μ为[0, 90000]。把一维实数空间上的集合[0, 90000]作为的x1论域,将连续论域[0, 90000]按等间距5000取点,转变为离散论域,根据文献中[1-3]的公式计算,可得五类风险源的超越概率估计值,代表不确定意义下的气象风险估计。
1.3 风险关注度计算举例以高温天气为例,2004年1月至2007年10月共有高温天气事件108个,持续出现2天70个、持续出现3天及以上48个,获得详细高温风险源所对应的12121气象电话拨打次数。设y j=(j=1, 2, …,108)为历次高温事件。考虑到计算精度的要求,取18个控制点(n=18),构成离散论域:U={u1u2,…,un}={0,5000,10000,…,90000},此时样本个数m=108,样本最大值b=89027,最小值a=9845,可算得扩散系数h= 1051.42。继而可以得出高温天气事件发生时的超越概率风险估计值。由于关注度越高、电话拨打次数越大,风险水平高,危险概率大。为了便于理解,转为计算(1-超越概率值)作为关注风险的累积概率值,来表征不同风险水平下的风险关注度,其余风险源计算方法类似,详见表1。
![]() |
表 1 不同气象风险源的风险关注度表 |
根据4次多项式拟合计算其启动点和拐点,可以定量对社会公众关注程度进行分类,继而得到不同风险等级的关注人数。
2 气象风险社会公众关注度等级分类以高温天气风险源为例,首先计算其拟合公式,由于次数低拟合程度不理想,次数高计算复杂,所以去除多余的“0.0000”和“0.0097”,均保留一个,风险水平为10000到65000,数据量为12个,通过excel添加趋势线获取四次多项式拟合公式(详见图 1)。即可求得“启动点”、“加速点”和“减速点”。“启动点”为四次拟合公式(拟合函数)的第一个有意义的极值点。由风险关注度计算方法及四次多项式性质可知,其第一个有意义的极值点必为极小值点。因此,这个点可用来表征社会公众关注度刚刚启动。“加速点”和“减速点”分布为四次多项式的二阶导数为零的两个点,即其凹凸性发生改变的拐点,可用其分别表征中、高等级社会公众关注度的分界点。中等社会公众关注度为从“启动点”到第一个拐点的区间,其间社会关注人数(拨打量)保持在中等水平,并快速增长;高等社会公众关注度为从第一个拐点到第二个拐点的区间,其间关注人数(拨打量)维持在较高水平,并持续增长,但增长速度明显减慢。
![]() |
图 1 高温风险源四次多项式拟合曲线图 |
高温多项式拟合公式为:P=0.0004x4- 0.0097x3+0.0725x2 -0.0263x-0.0616,其中R2=0.9920。其中X是拟合过程中的一个数据,B 为风险水平,计算得到关注度概率P。
启动点:X1=0.18834;B1=10000+(X1-1)×5000=5942;P1 =0.0641
拐点:X2=3.34024;B2=10000+(X2-1)×5000=21701;P2= 0.342931
X3=9.63293;B3=10000+(X3-1)×5000=53165;P3=0.936225
其余风险源的计算过程以此类推,得到不同气象风险源的启动点和拐点分布图(详见表2)。需要指出的是:由于中雨和大到暴雨两种气象风险源的拟合多项式的第一个极值点处于其定义域之外,故其“启动点”为其拟合多项式的第二个极值点。这一拟合函数的定义域内就不存在第二个拐点。出现这种情况的原因可能是中雨和大到暴雨这两种气象风险源的特殊性造成的,即所有人对其都非常地关注,一旦出现这两种天气状况,人们的关注程度将不断上升,直到扩大到所有公众。
![]() |
表 2 不同气象风险源的启动点和拐点分布图 |
根据《2006年陕西统计年鉴》的数据,西安市常住人口533.21万人,全市拥有固定电话用户324.04万,每百人固定电话普及率60.77%,依据2007年度陕西省气象局开展的相关社会调查结果显示,约有20.45%的社会公众通过电话获得气象信息服务。根据实际情况和多项式求解的数学意义,定义启动点以上为一般等级社会关注度,拐点1对应中等级社会关注度,拐点2对应高等级社会关注度。依据表2中三级关注度对应的不同风险概率值,高温时社会公众关注人数计算如下:
一般等级关注度以上的关注人数:533.21×20.45%×0.0641=6.9896万人;
中等级关注度以上的关注人数:533.21×20.45%×0.342931=37.3937万人;
高等级关注度以上的关注人数:533.21×20.45%×0.936225=102.0873万人;
以此类推,得到西安社会公众关注高影响天气事件在不同风险等级的分布人数,如表3所示。图 2给出了更为直观的比较说明。研究发现,在研究期内,西安市533多万常住人口中,一旦发生高温天气约有6.99万人就开始关注,对于小雨、中雨和大到暴雨天气事件发生过程中,分别至少约有44万、60万和48万人比较关注,超过107万西安公众很关注雷阵雨天气的发生演变。
![]() |
表 3 西安气象风险社会公众在不同等级气象关注度的关注人数(单位:万人) |
![]() |
图 2 西安气象风险社会公众不同等级气象关注人数堆积图 |
以高影响天气事件作为重要气象风险源,综合应用气象热线电话拨打量信息,通过计算出超越概率和累积概率值,来表征不同风险水平下的风险关注度。通过4次多项式的拟合和求解,定量计算了风险关注度等级。根据实际情况和多项式求解的数学意义,通过求一阶导得到其“启动点”,表征一般等级社会公众关注度刚刚启动的点;然后,通过求二阶导计算其“加速点”和“减速点”,对应着4次多项式的两个拐点,分别表征中、高等级社会公众关注度的分界点。
综上分析,定量计算了不同等级气象风险的关注人数。研究发现,在研究期内,一旦高温、小雨、中雨和大到暴雨和雷阵雨天气事件开始发生,西安社会公众分别有不同人数开始关注,当事件持续演变到高风险关注程度时,均有超过102万西安公众十分关注各类高影响天气事件的进一步发展演变,关心其是否成灾对生产生活所造成的影响等。通过公共气象服务平台发布定量计算的风险关注度和风险水平等级,有利于提醒政府决策者和社会公众关注高影响乃至灾害性天气气候事件,关注其对日常生活、生产、交通运输和旅游的影响,提高对各种气象风险的防灾减灾意识并及时采取防范应对措施,从而达到预防和减轻其所造成负面影响的目的。
[1] |
叶笃正, 严中伟, 戴新刚, 等. 未来的天气气候预测体系[J]. 气象, 2006, 32(4): 3-8. |
[2] |
罗慧, 张雅斌, 刘璐, 等. 高影响天气事件公众关注度的风险评估[J]. 气象, 2007, 33(10): 15-22. |
[3] |
黄崇福. 自然灾害风险评价理论与实践[M]. (第一版). 北京: 科学出版社, 2006: 2-45.
|
[4] |
Huang Chongfu. Principle of information diffusion[J]. Fuzzy Sets and Systems, 1997, 91: 69-90. DOI:10.1016/S0165-0114(96)00257-6 |