径流随机模拟

1、径流随机模拟由于实际水文资料往往比较短，难于满足实际水文程作随机模拟。这种随机模拟的目的之一在于充分利用是用来延长资料长度。当所建模型及参数准确时，这种年月径流随机模拟方法，对于多站及更深入的随机模拟一、随机过程基本知识 (一)随机过程和时间序列的定义 在实际问题中，常涉及试验过程随某个参变量的变的流量、水位是随时间变化的随机变量，气温是随时间这种随机变量为随机函数，并称以时间为参数的随机函数为随机过程，记为，T是t变化的范围。 随机过程在一次试验或观测中所得结果，称为随机过程的一个实现。 若时间参变量T是连续时刻的集合，则称这种随创机过程为连续参数随机过程，如水位过程、流量过程等。若时间参变

2、量T是程为离散时刻的集合，则称这种随机过离散参数随机过程，也称为随机序列或时间序列。如年、月径流程，年最大流量过程都是时间序列，也称水文时间序列。 (二)随机过程的数字特征随机过程在任一固定时刻的状态是随机变量，因此可按与前述随机变量同样的方法定义随机过程的数学期望和方差。定义如下数学期望 (811)方差 (812)为了规划随机变量两个不同时刻状态间关系的密切程度，可定义随机变量的自相关函数为 （8-13） (三)随权过程基本分类 l。按统计性质的稳定性分类 按随机过程的统计性质是否随时间而变化，可分成平稳和非平稳过程。若随机过程统计数字特征不随时间的平移而变化，则称为平稳过程，否则为非平稳过

3、程。 2按不同时刻状态间的关系分类 可分成独立过程和马尔柯夫过程。若过程各状态相互独立，则称为独立随机过程。在非独立随机过程中，最重要的一类是马尔柯夫过程，其特点是时刻状态只与时刻有关，而与以前各时刻无关。以上各个概念的严格定义，可参见有关文献。二、径流随机模拟一船步骤图87给出了单站年径流随机模拟的一般步骤：时间序列组成分析；模型的建立；序列的生成；模型及生成系列的检验。三、水文时间序列的组成分析水文序列一殷可按下式表示： （8-14）Tt，Ct，Pt，St分别为趋势项、跳跃项、周期项和随机项。 当水文序列中不含Tt，Ct，Pt等确定性成分时，如St，即仅包括随机成分的序列。对年径流序列而言

4、，这种情况是出较常见的。但月径流序列因存在明显的年周期，所以不是仅包括随机成分的序列。 趋势项指的是水文变量的统计参数(主要是均值)的长期系统性升降，如图88即存在明显趋势(增长)。一般是由于气候因子或下垫面因子逐步改变而引起的缓慢变化。对实测水文序列，可用假设检验或滑动平均的方法查明是否存在趋势。若存在趋势，呈线性变化时，常用线性方程拟合，然后从序列中将趋势滤掉。跳跃项是指水文序列急剧变化的一种形式， 当水文序列从一种状态过渡到另一种状态时表现出来。见图8-9。跳跃是人为或自然原因造成的。如建库后面积增大，蒸发量等损失增加，有可能出现跳跃，并反映在年径流序列的均值等参数。突变可看作跳跃的一种

5、特殊情况，如由于地震塌方，拦截江河，形成水库以后溃坝，这样引起流量的突变，随着临时水坝的冲毁，又恢复到原来状态。 跳跃是否存在于序列中，多用分割样本的方法检验。若存在较显著的跳跃成分，应从序列中排除掉，使得剩余序列具有原始状态或一致条件。 周期项(含近似周期)是由于天体运动的周期性影响造成的，如地球公转、自转引起以年和日为周期的变化，以及太阳黑子活动引起的旱涝多年变化。通常可用谐波分析的方法析出，再从序列中滤掉。 随机成分是由于不规则及随机振荡引起的，一般由相依成分和纯随机成分组成。严格地讲，几乎所有水文变量应是非平稳过程。不过除了人为影响及自然灾变外，水文环境的变化在数十年或几百年期间都相当

6、小，因此，从实用观点，常把水文序列中的随机成分看作平稳过程。图88序列趋势变化示意图图89跳跃成分及包含跳跃成分的合成序列四、单站年月径流随机模型的建立 对年月径流序列建立随机模型一般是对原始年月径流序列排除趋势跳跃等确定性成分后的随机成分而建立的模型。设用于建模的年径流序列为，。 径流序列随机模型是实测时间序列一种概化的数学式子。目前常见的随机模型有：线性平稳模型、非线性平稳模型，马尔柯夫模型及非马尔柯夫模型等。对年径流序列而言，常采用线性平稳模型或马尔柯夫模型。 建立随机模型的一般步骤为：选择模型；确定阶数；估计模型参数 (一)单站年径流随机模型的建立 通常采用线性自回归模型，即马尔柯夫模

7、型 1线性自回归模型的一船形式 （8-15） 第t年的年径流量，t=1，2，常称式（8-15）的为自回归系列。序列的平均值；，自回归系数或偏相关系数，反映在时间上相依性大小； 模型残差项，纯随机成分， 与，无关且是独立随机变量，其均值为0，方差为。 由于与的方差有确定关系，因此，一般自回归模型中参数有：、和，共p十2个参数。 该模型说明第t年年径流量仅依赖于第t一1年、第t一2年、第t一p年的年径流量和一个纯随机变量。 若令，则式(815)变为： (816)式(816)是中心化变量表示的自回归模型。2模型参数的估计 （8-19）其中K阶样本自相关系数在n较大，K较小时，计算公式为 （8-20）

8、据推导 （8-21）在数学上序列，一般假定为正态分布，故不需计算其偏态系数。但对于具有偏态的水文系列，一般把当作P-型分布，因此还必须估计的偏态系数。 （8-22）（t=P+1，P+2，n）是根据估计出的以上P+2个参数及观测序列，利用公式（8-15）反推得到的。平均值是根据反推序列利用矩法估计的。3常见AR（1）及AR（2）模型参数估计公式AR（1）模型形式为： （8-23）模型参数估计公式为： （8-26）AR（2）模型形式为： （8-31）4模型阶数P的确定对于AR（P）序列，可以证明：它的自相关系数随滞时增大而减小，呈拖尾状，而偏相关系数则呈截尾状，在时出现一个截止点，即在时，当时，。

9、因此从理论上讲，可以通过计算不同的进行模式识别。例如，当从样本序列估计在时具有明显的截尾现象，那么可以推断该水文序列即适合于AR（3）模型。但是由于实际水文样本容量较小，故统计量抽样误差较大，即使是AR（P）序列，当时，可能并不为零，这样就难于做直观判断，必须进行统计推断。统计推断方法是：取显著水平，若，则可接受异于0的假设。例如，某河流年径流偏相关系数如图8-10所示，该图表明，只有超过95%容许极限，即，故该模型阶数应为1，即AR（1）模型。（二） 单站月径流随机模型的建立对已有n年实测月径流资料的单站通常有两条途径建立其月径流随机模型。一是先建立年径流模型，再通过建立解集模型把年径流分解

10、成各月月径流。解集模型公式为： （8-32）公式中 -各月月径流流量，；-年径流量；-模型参数，反映各月月径流流量平均分配水平；-12×12的参数矩阵，反映各月之间的相关关系程度；-模型残差项，相互独立，可以是正态或偏态分布。以上参数n由年实测资料估算。本模型结构简单，概念清晰，但因参数多，故所需实测资料较长。二是直接建立月径流随机模型，通常采用季节性一阶自回归模型，即假定可用12个一阶自回归模型来描述各月月径流流量及相关系数。各月月径流模型如下： （8-33）公式中-年份，；-月份，；-第年第月的月径流量；，-第月的月径流均值和均方差，；-第月和第-1月月径流之间的相关系数，表示第

11、一月和上一年第十二月月径流量相关系数；-第月纯随机变量，是模型残差想，可以是标准正态分布或标准P-分布，各月之间相互独立，且与相互独立。以上各有关参数可由n年实测月径流资料用矩法估算。若采用标准正态分布，则月径流量也是正态分布。若采用标准P-分布，则月径流量为近似P-分布，生成P-分布时，还需估算的偏态系数 （8-34）公式中-第月径流偏态系数，时，表示的是第12月的偏态系数，该参数可用实测n年第月月径流序列估算。五、年月径流序列的生成（模拟）年月径流序列模拟涉及到纯随机变量的生成及由生成年月径流序列的问题。（一） 纯随机变量随机数的产生纯随机变量的分布可以是正态，也可以是偏态。它的生成一般先

12、生成0，1均匀分布随机数u，再通过变换生成指定分布的随机数。 1均匀分布随机数的生成 生成方法有随机数表法、物理方法及数学方法。由于前两种方法存在严重缺陷，故常用数学方法生成，其中应用最广的是乘同余法。 乘同余法生成随机数递推公式是： (835) (836) 式(835)、式(836)中为初值，为乘子，为模，它们均为非负整数，而且M。是被M整除后的余数，于是M，故即为0，1上的随机数。 这种方法生成的随机数存在着循环周期，因此，不是真正意义上的随机数，俗称“伪随机数”。但由于M往往取值很大，周期也很长，目前微机上周期可达109以上，实用上完全能满足需要。正因如此，实际大都使用该法生成0，1均匀

13、分布随机数。不过使用前要对生成的随机数作均匀性、独立性等检验。 2正态分布随机数的生成 通常用BoxMuller变换生成，即 (837)中(837)中，为0，1上的随机数，为相互独立标准化正态随机数。对于任意正态分布，随机数 生成公式，为标准化正态分布随机数。 3P型分布随机数生利用舍选法生成P型分布随机数，详见图8一11。 其中（）为为0，1均匀分布随机数，z为所需生成P型分布随机数。该分布3个参数为：均值，离势系数及偏态系数。在，已知的情况下，框图中n、p、a、的计算公式如下： (838) (839)或 （8-40） （8-41） (二)年月径流的生成下面仅介绍年径流生成的方法，月径流生成

14、的方法类同。设所建立的模型为AR(1)，且参数已估计出。其模型为 （8-42）式(842)中分布参数：均值为0，均方差为、，偏态系数为。下面分两种情况介绍年径流生成步骤：1. 考虑年径流为正态分布这种情况下，为正态分布，即可认为等于0。生成步骤：(1)以或(tl、2、n)为；(2)生成一个符合的正态随机数；(3)以及代入式(842)生成一个年径流Q1；(4)同步骤(2)，生成一个；(5)以Q1和代人式(842)，计算出Q2；(6)重复上述步骤，可得到一个很长生成序列，如容量为的序列；(7)考虑到前50项可能受初值影响，应舍去，故剩下年生成的年径流系列。2考虑年径流为偏态分布这种情况一般考虑为P

15、型分布，三个参数是：均值（均值为0），方差，偏态系数。生成年径流序列的方法与考虑年径流为正态分布时几乎一样，唯一不同的是上述第二步改用了P型分布随机数生成。这样生成的序列可近似认为是P型分布。六、模型及生成系列的检验 模型检验是指所选用的模型残差为独立随机变量是否成立，分布是否为假定分布。生成序列检验是指所生成年月径流序列是否能反映实测径流的统计特性。 1残差独立性检验在模型及参数确定后，根据实测样本，用式（8-15）可推出残差序列（t=P+1，n），由序列可计算其各阶自相关系数，再对作独立性假设检验。当检验通过，即满足独立性，说明建模时对独立性假定是成立的，否则要分析产生的原因。若序列存在密

16、切自相关，应考虑使用其他模型。 2生成系列检验 一般要求生成序列与实测序列统计特征相近。如差异很大，要分析原因，确实是模型结构问题，应考虑改变模型。七、实例红水河龙滩站年径流序列生成红水河龙滩水库为一多年调节水库，坝址处有自19461979年共34年的资料。为了分析工作需要，要求生成年径流系列。1径流组成分析经过分析，末发现有趋势、突变、周期等确定成分，故可直接对实测序列建立平稳模型。2模型选择和参数估计为了便于选择模型，已估算了年径流序列统计参数、及自相关系数、，见表83及图812(a)。从图812(a)可看出，自相关系数呈指数衰减趋势，故选用常用的AR(P)模型作为年径流模型。为了确定模型

17、阶数，计算了（），互1，2，15)见图812(b)。从图中可知，K3时，超过置信限范围，而K3时，均在置信限范围内，说明在K3久可以认为是0，因此，应选定阶数为3。 模型参数利用前面介绍的矩法估算，结果为：一0021，0143，0333，0931。此外，利用实测系列仅推出、计算出的，结果是接近于0。 3年径流序列生成 由于很小，故把当作正态分布，即用正态分布生成，利用三阶自回归模型生成1000F年径流系列。 4模型及生成系列的检验 (1)对于反推序列、，求自相会系数，经检验可认为是独立的，因此，模型的假定是成立的。 (2)生成系列的检验。对生成1000F年径流系列计算、及自相关系数、，见表83

18、。通过对比发现两序列各项统计参数相对误差较小。 以上检验结果表明，所建模型是可接受的。 总之，径流随机模型是随着数学中随机过程理论和电子计算机技术在水文学中应用而逐渐发展起来的，所有各种随机模型都是建立在水文现象的统计特性基础上，模型中参须基于实测资料所提供信息加以估计，因此，只有深入了解水文现象的特性以及获得尽能多的可靠信息，才能更有效地使随机模型应用于工程水文中。图812龙滩站自、偏相关函数示意图第七节 设计枯水流量分析计算（未吃透）枯水流量亦称最小流量，是河川径流的一种特殊形态。枯水流量往往制约着城市的发展规模、灌溉面积、通航的容量和时间，同时，也是决定水电站保证出力的重要因素。按设计时

19、段的长短，枯水流量又可分为瞬时、日、旬、最小流量。其中又以日、旬、月最小流量对水资源利用工程的规划设计关系最大。一、有实测水文资料时的枯水流量计算 当设计代表站有长系列实测径流资料时，可按年最小选样原则选取一年中最小的时段径流量，组成样本系列。 枯水流量常采用不足概率q，即以小于和等于该径流的概率来表示，它和年最大选择的概率P有ql一P的关系。因此在系列排对时按右小到大排列。除此之外，年枯水流量频率曲线的绘制与时段径流频率曲线的绘制基本相同，也常采用P-型频率曲线适线。图8-13为某水文站不同天数的枯水流量频率曲线的事例。 年枯水流量频率由线，在某些河流上，特别是在干旱半干旱地区的中小河流上，还会出现