径流模拟方法探究——P-Ⅲ频率曲线与copula函数

方法简介

关于三峡梯级水库群径流资料模拟的方法，尝试使用copula函数构建不同水库区间来水的相关函数，进行径流模拟。因此需要求得各区间来水的边际分布函数，套用P-Ⅲ频率分布曲线，进行适线法拟合得出各水库各年各旬来水的边际分布函数，而后准备进一步进行copula函数拟合。

copula函数

Copula 是定义域为 [0,1] 均匀分布的多维联合分布 函数, 它可以将多个随机变量的边缘分布连接起来得到它们的联合分布

Sklar 定理: 令 H 为一个 n- 维分布函数, 其边缘分 布为 F1, F2, …Fn。则存在一个 n- Copula 函数 C, 使得 对任意 x∈Rn
H( x1,x2,…xn) =C( F1( x1) , F2, ( x2) , …, Fn( xn) ) ——（1）
如果 F1, F2, …, Fn 是连续的, 则 C 是唯一的, 相反地, 如果 C 是一个
n- Copula, F1, F2, …, Fn 为分布函数,
则式( 1) 中所定义的函数 H 是一个 n- 维分布函数, 其 边缘分布为 F1, F2, …, Fn

使用工具

EXCEL、P-Ⅲ线型适线法软件（VB）

具体实施

根据已提供的74年三库区间逐旬的来水资料，以旬为组，每库划分36组资料，每组为74年单库单旬数据，对每组数据进行频率适线，得到各库单旬的P-Ⅲ型频率分布曲线，通过查图可得其不同保证率下的区间来水量，以P=10%，50%，90%作为控制条件，对应丰、平、枯三种典型年状况，分别查的各库某旬的来水量大小，每库有3组，共9组区间来水数据。






已有成果
1、通过适线法采用模拟典型年方法得到的径流模拟系列，与老师提供的以总水量划分丰平枯年组选取代表年方法数据对比，来水曲线更加平缓没有大量锯齿，综合考虑了74年数据的影响，使得趋势更加有代表性。已经可以作为来水模拟的一种方法。
2、得到各库逐旬的P-Ⅲ频率分布曲线的各项参数，得到其编辑分布函数，以此可以进一步求解得，以每两库逐旬来水的边际分布函数所构成的联合分布函数，即copula联合分布，进一步归纳各库来水逐旬相关性。

存在问题
1、采用的这个适线法求频率构造典型年的方法，人为影响比较大，枯季和汛期不同也会对其适线有影响，所得数据缺乏定量校核。
2、在后续构造copula函数的理论中，目前想选用GH型函数，目前这种函数大多用于构造不同种变量间的相关关系函数，如洪峰和洪量，洪量与洪水历时等，对于上下游区间来水同种变量关系的构造尚不多见，但我认为他们确实有很强的时空关联性，可以尝试。

后续工作
1、继续进行来水变量相关关系的构造，做出各库之间的逐旬的来水相关性函数，通过单库数据可以直接模拟后续区间来水数据。
2、进行多目标联合优化调度增益分配机制研究，基于Nash讨价还价模型展开。