Page 109 - 2022年第53卷第8期
P. 109

表 1 4站日降水量统计值及模型参数
                                    实测日降水量序列统计值                                    模型参数
                  站点
                                                            r
                            x mean  C v     C s     R m1     1      π 0     π 1     λ 1     λ 11     β 11
                  51330     0.7    3.610   6.911    57.0   0.152   0.730   0.270   0.291   0.445    0.407
                  57278     2.3    3.469   8.188   293.9   0.174   0.690   0.310   0.367   0.503    0.268
                  57512     3.1    3.082   7.125   232.1   0.151   0.554   0.446   0.311   0.372    0.311
                  59849     6.6    2.861   6.227   373.5   0.380   0.462   0.538   0.421   0.687    0.278

              注:λ 1 代表 DAR(1)模型参数;λ 11 和 β 11 分别代表 DARMA(1,1)模型参数

                  所选择的 4个站点多年日平均降水量分别为 0.7、2.3、3.1和 6.6mm,且地理位置横跨东西南北。
              通过日降水量统计值可以发现,不同站点 C和 C值无明显差异,沿海地区的最大一日降水量大于内陆地
                                                         s
                                                     v
              区,大部分地区不降水概率大于降水概率,而降水概率越大,日降水量序列自相关性越大。
                  由式( 5)和式(13)分别确定 DARMA(1,1)和 DAR(1)模型的理论 ACF,图解法比较样本和理论
              ACF如图 2。总体而言,各站点的理论 ACF衰减缓慢,最终在第 8天附近趋于零,符合 DARMA(1,1)
              模型的特征,证实 DARMA(1,1)模型适用于中国气象站点的日降水序列模拟。DAR(1)模型的 ACF
              呈指数衰减,以站点 59849最为明显,在 ACF滞时为 2后出现明显衰减,导致与样本数据 ACF出现较
              大误差。理论 ACF检验结果显示,对于中国降水事件干湿序列的模拟,DARMA(1,1)更能保持较高
              滞后阶数的相关性。

































                                                  图 2 理论和实测 ACF对比图

                  在确定模型结构和阶数时,湿游程和干游程同样重要                          [29] 。因此,根据干湿游程样本和理论概率分
              布的平方和误差的最小值进行模型选择。通过式( 19)—(30)计算 DARMA(1,1)模型和 DAR(1)模型
              的理论干湿游程概率分布,与样本序列进行对比,如图 3所示。
                  从图 3可知:(1)DARMA(1,1)模型能够在 1到 10个连续雨日和连续干日内生成概率误差最小,
              优于 DAR(1)。对于 51330站,DARMA(1,1)模型单个干日概率为 0.2427,而实测数据计算概率为
              0.2368;对于 51330站,DARMA(1,1)模型连续 4个雨天的概率为 0.0903,而实测数据计算概率为
              0.0943。(2)DAR(1)模型连续干日概率分布的模拟效果较差。对于 57278站,DAR(1)模型单个干日
              概率为 0.1959,而实测数据计算概率为 0.2555;对于 59849站,DAR(1)模 型 连续 4个干 日 概率为

                     6
                —  9 9  —
   104   105   106   107   108   109   110   111   112   113   114