随机水文学-第2章

随机水文学随机水文学-第第2 2章章2.1 2.1 随机过程的概念随机过程的概念注：(1)随机过程 X(t)或时间序列 Xt 中的t通常表示时间，但也可以表示空间、长度等其他非时间变量。(2)随机过程和时间序列在许多方面存在相互平等的理论，但两者不完全相同，在实际工作中有时没有加以严格区分。(3)不管是随机过程还是时间序列，它的一个显著特征是预测值前后的相依性，因此，数据的时间顺序十分重要。(4)一个随机变量 X 与它的一个样本值 x，意义不同。2.2 2.2 随机过程的分布函数随机过程的分布函数维数越大，随机过程的统计特性的描述也越趋完善。一般认为，分布函数族(F1，F2，)完全地确定了随机过程的全部统计特性。2.3 2.3 随机过程的数字特征随机过程的数字特征5.变差系数描述随机过程对于均值的相对偏离程序。6.偏态系数描述随机过程对于正态分布的相对偏离程序。均值、方差、变差系数、偏态系数描述了随机过程 X(t)在各个孤立时刻 t 的统计特性，称为随机过程的统计参数。统计参数的值反映了随机过程的数字特征。统计参数是随时间变化的。如表2-12.4 2.4 随机过程的分类随机过程的分类按是否相依：相关随机过程、独立随机过程按变量多少：多变量随机过程、单变量随机过程按参数是否随机时间变化：平稳随机过程、非平稳随机过程按随机过程的分布函数的不同性质：独立随机过程、平稳随机过程、独立增量随机过程、马尔柯夫随机过程1.独立随机过程随机过程任意时刻的截口与其他时刻的截口之间互不影响。相关系数为零2.独立增量随机过程在任意时间间隔上过程状态的改变并不影响未来任一时间间隔上过程状态的改变。2.5 2.5 平稳随机过程平稳随机过程平稳随机过程在水文水资源计算中最平稳随机过程在水文水资源计算中最为常用。为常用。一、平稳随机过程的概念 一个随机过程 X(t)，若对任何 n 与 k，X(t)的 n 维分布函数满足：则X(t)被称为被称为平稳随机过程平稳随机过程，否则，否则被称为被称为非平稳随机过程非平稳随机过程。平稳随机过程的 n 维分布函数不因所选开始时刻的改变而不同，即平稳随机过程的统计特性与所选取的时间起点无关。也就是说，平稳随机过程的统计特性不随时间 t 的变化而改变。例如：利用相当长的年径流序列，计算的任一年径流量n维分布函数是相等的。若产生随机过程的主要物理条件在时间进程中没有变化，则该随机过程的统计特性也不会随机随时间而变化。如果产生年径流的气候条件与下垫面条件都没有重大变化，则年径流的统计特性也不会随时间而变化，因而不同开始时刻的年径流n维分布函数也不会有变化。平稳随机过程具有一系列简单的特性，使问题的分析计算大为简化，实际应用广泛。二、平稳随机过程的数字特征1、均值平稳根据平稳过程的定义，当n=1时，对任意有当=-t时时，即平稳过程 X(t)的一维分布函数和一维概率密度与时间t无关。均值与时间t无关，为常数。平稳随机过程的均值平稳，又称一阶平稳。2、方差平稳、方差平稳方差函数平稳随机过程 X(t)的方差函数2与时间t无关，为常数，称为方差平稳，又称二阶平稳。同理标准差函数(x)也是平稳的。3、偏态系数平稳、偏态系数平稳平稳随机过程 X(t)的偏态系数与时间t无关，为常数，称为偏态系数偏态系数平稳。4、协方差平稳、协方差平稳根据平稳过程的定义，当n=2时，对任意 k 有令k=-t1，t2-t1=时时，即平稳过程 X(t)的二维分布函数与时间t无关，只与时间间隔有关有关。自协方差函数与时间t无关，只与时间间隔有关有关。称自协方差平稳。5、自相关函数平稳、自相关函数平稳平稳随机过程 X(t)的自相关函数与时间位置无关，只与时间间隔有关，有关，称自相关函数平稳。三、平稳随机过程的分类平稳随机过程可分为两类：u严平稳随机过程：即满足定义的平稳随机过程，又称狭义平稳过程或高阶平稳过程。现实中不存在u宽平稳随机过程：即均值和协方即均值和协方差平稳的过程，也称广义平稳过程差平稳的过程，也称广义平稳过程或二阶平稳过程。或二阶平稳过程。一般平稳过程如一般平稳过程如没加特别说明都是指宽平稳过程。没加特别说明都是指宽平稳过程。水文水资源系统中，当影响它的主水文水资源系统中，当影响它的主要因素（气候、下垫面及人类活动要因素（气候、下垫面及人类活动等）相对稳定时，以年为时间尺度等）相对稳定时，以年为时间尺度的水文序列可近似作为平稳随机序的水文序列可近似作为平稳随机序列。列。年降水量序列、年径流序列、年降水量序列、年径流序列、年蒸发量序列年蒸发量序列等。等。四、平稳随机过程的各态历经性通过大量的样本函数计算的各截口数字特征，能真实反映随机过程的统计特性。在水文学中，经常难以获取大量的样本函数，实际上仅仅有其中一个样本函数。在这种情况下，能否用一个样本函数来分析随机过程的统计特性？1、平稳随机过程的各态历经性【定义】在一定条件下，平稳随机过程的一个相当长的样本资料可以用来分析计算平稳随机过程的统计特性。这样的随机过程称为具备各态历经性或遍历性。该随机过程称为各态历经过程。在样本容量很大的情况下，每个样本函数能够代表过程的所有可能样本函数。因而任何一个样本函数都可以代表平稳过程的统计特性，可以由任何一个样本函数估计平稳过程的统计特征。非平稳随机过程不具备各态历经性，平稳随机过程也不全具备各态历经性。水文水资源系统，假定平稳随机过程具有各态历经性。2、基于时间域的数字特征计算 设平稳随机过程X(t)的任意一个样本函数x(t)(0tT)当 T 或 n 足够长时，平稳过程计算的数字特征，时间平均等于统计平均。2.6 2.6 马尔柯夫过程马尔柯夫过程一、马尔柯夫过程的定义及特征一、马尔柯夫过程的定义及特征 若随机过程 X(t)满足：则X(t)被称为被称为马尔柯夫过程马尔柯夫过程。称为称为马尔柯夫过程从时刻马尔柯夫过程从时刻 tn 状态状态 Xn，转移到时刻，转移到时刻 tn+k 状态状态 Xn+k 的概率，的概率，简称转移概率。简称转移概率。在在 tn 时刻所处的状态已知的条件下，时刻所处的状态已知的条件下，马尔柯夫过程在时刻马尔柯夫过程在时刻 tn+k 所处的状态只所处的状态只与其在与其在 tn 时刻所处的状态有关，而与时刻所处的状态有关，而与其在其在tn 时刻以前所处的状态无关。这种时刻以前所处的状态无关。这种特性称为马尔柯夫过程的无后效性。特性称为马尔柯夫过程的无后效性。过程过程“现在现在”的状态已知，其的状态已知，其“将将来来”的状态与的状态与“过去过去”的状态无关。的状态无关。马尔柯夫过程的统计特性完全由它马尔柯夫过程的统计特性完全由它的的初始分布和转移概率初始分布和转移概率确定，因此，确定，因此，研究马尔柯夫过程，只需确定初始分研究马尔柯夫过程，只需确定初始分布和转移概率。布和转移概率。马尔柯夫过程可以分为三类：马尔柯夫过程可以分为三类：n时间和状态都连续的马尔柯夫过程时间和状态都连续的马尔柯夫过程维纳过程维纳过程n时间连续、状态离散的马尔柯夫过程时间连续、状态离散的马尔柯夫过程散粒噪声过程散粒噪声过程n时间和状态都离散的马尔柯夫过程时间和状态都离散的马尔柯夫过程马尔柯夫链马尔柯夫链马尔柯夫链是最简单的马氏过程，马尔柯夫链是最简单的马氏过程，在水文学中广泛应用在水文学中广泛应用二、马尔柯夫链设马尔柯夫链有设马尔柯夫链有m个状态个状态（如径流的特丰、丰、中、枯、特枯）（如径流的特丰、丰、中、枯、特枯），记转移时刻为。某一时，记转移时刻为。某一时刻的状态为刻的状态为m个状态之一。个状态之一。为过为过程从时刻程从时刻 tn状态状态 ai 经过经过 k 步转移到状步转移到状态态 aj 的概率。的概率。一般情况下，一般情况下，pij(n,k)与与 i,j,k和和 n 有有关，当关，当pij(n,k)与与n(初始时刻初始时刻)无关时，无关时，称为齐次马尔柯夫链。称为齐次马尔柯夫链。实际工作中，一般考虑齐次马尔柯实际工作中，一般考虑齐次马尔柯夫链。夫链。取取k=1，则，则 称为一步转移概率称为一步转移概率。由一步转移概率可构成一步转移概率由一步转移概率可构成一步转移概率矩阵矩阵 式中：式中：当当k2时变成多步转移概率矩阵。时变成多步转移概率矩阵。一步转移概率矩阵与多步转移概率矩一步转移概率矩阵与多步转移概率矩阵关系：阵关系：令时刻令时刻t的无条件概率分布或边际概的无条件概率分布或边际概率分布为率分布为 Pt=pt(1),pt(2),pt(m)若时刻若时刻t已发生，则已发生，则Pt 已知，已知，t+1时刻时刻的条件分布为的条件分布为 年份年份年份年份195219521953195319541954195519551956195619571957195819581959195919601960196119611962196219631963年径流年径流年径流年径流540540478478466466273273378378422422251251508508307307465465375375190190年份年份年份年份196419641965196519661966196719671968196819691969197019701971197119721972197319731974197419751975年径流年径流年径流年径流404404279279336336351351570570280280528528374374329329515515356356432432年份年份年份年份197619761977197719781978197919791980198019811981198219821983198319841984198519851986198619871987年径流年径流年径流年径流466466499499386386395395386386445445434434480480314314335335303303382382年份年份年份年份198819881989198919901990199119911992199219931993199419941995199519961996199719971998199819991999年径流年径流年径流年径流301301282282352352260260418418568568633633405405455455500500518518411411例：桂江流域中游控制站平乐站例：桂江流域中游控制站平乐站48年年(1952-1999年年)径流径流量资料。量资料。经计算，年径流样本均值经计算，年径流样本均值402 m3/s，样，样本标准差本标准差96.2 m3/s。将年径流划分为将年径流划分为5个状态：枯、偏枯、个状态：枯、偏枯、平、偏丰、丰。分别用平、偏丰、丰。分别用1，2，3，4，5表示。表示。枯枯枯枯 0,305.8 0,305.8偏枯偏枯偏枯偏枯(305.8,353.9(305.8,353.9平平平平(353.9,450.1(353.9,450.1偏丰偏丰偏丰偏丰(450.1,498.2(450.1,498.2丰丰丰丰(498.2,+)(498.2,+)解题过程：解题过程：1、列出平均流量状态表。列出平均流量状态表。年份年份年份年份195219521953195319541954195519551956195619571957195819581959195919601960196119611962196219631963年径流年径流年径流年径流540540478478466466273273378378422422251251508508307307465465375375190190状态状态状态状态5 54 44 41 13 33 31 15 52 24 43 31 1年份年份年份年份196419641965196519661966196719671968196819691969197019701971197119721972197319731974197419751975年径流年径流年径流年径流404404279279336336351351570570280280528528374374329329515515356356432432状态状态状态状态3 31 12 22 25 51 15 53 32 25 53 33 3年份年份年份年份197619761977197719781978197919791980198019811981198219821983198319841984198519851986198619871987年径流年径流年径流年径流466466499499386386395395386386445445434434480480314314335335303303382382状态状态状态状态4 44 43 33 33 33 33 34 42 22 21 13 3年份年份年份年份198819881989198919901990199119911992199219931993199419941995199519961996199719971998199819991999年径流年径流年径流年径流301301282282352352260260418418568568633633405405455455500500518518411411状态状态状态状态1 11 12 21 13 35 55 53 34 44 45 53 32、统计一步转移频率矩阵。、统计一步转移频率矩阵。3、一步转移概率矩阵。一步转移概率矩阵。4、估算概率。估算概率。1999年为径流平水年年为径流平水年,则无条件概率分布为则无条件概率分布为 P1999=0,0,1,0,0则则2000年径流概率分布为年径流概率分布为年份年份年份年份196919691970197019711971197219721973197319741974197519751976197619771977197819781979197919801980平均流量平均流量平均流量平均流量662662656656542542576576446446719719644644650650635635544544701701776776年份年份年份年份198119811982198219831983198419841985198519861986198719871988198819891989199019901991199119921992平均流量平均流量平均流量平均流量677677737737533533500500773773480480681681582582777777779779672672583583年份年份年份年份199319931994199419951995199619961997199719981998199919992000200020012001200220022003200320042004平均流量平均流量平均流量平均流量876876471471569569662662533533760760774774880880713713453453830830737737 某河年平均流量资料见下表，当某河年平均流量资料见下表，当Q 535m3/s时，为枯水年；当时，为枯水年；当535m3/sQ 775m3/s时，为中水年；当时，为中水年；当Q 775m3/s时，为丰水年时，为丰水年。试。试估算估算2005年径流为丰、中、枯状态的概率。年径流为丰、中、枯状态的概率。年份年份年份年份196919691970197019711971197219721973197319741974197519751976197619771977197819781979197919801980平均流平均流平均流平均流量量量量662662656656542542576576446446719719644644650650635635544544701701776776状态状态状态状态2 22 22 22 23 32 22 22 22 22 22 21 1年份年份年份年份198119811982198219831983198419841985198519861986198719871988198819891989199019901991199119921992平均流平均流平均流平均流量量量量677677737737533533500500773773480480681681582582777777779779672672583583状态状态状态状态2 22 23 33 32 23 32 22 21 11 12 22 2年份年份年份年份199319931994199419951995199619961997199719981998199919992000200020012001200220022003200320042004平均流平均流平均流平均流量量量量876876471471569569662662533533760760774774880880713713453453830830737737状态状态状态状态1 13 32 22 23 32 22 21 12 23 31 12 2 解题过程：解题过程：1、列出平均流量状态表。、列出平均流量状态表。2、统计一步转移频率矩阵。、统计一步转移频率矩阵。3、写出一步转移矩阵。、写出一步转移矩阵。4、估算概率。、估算概率。2004年为径流中水年年为径流中水年,则无条件概率分布为则无条件概率分布为 P2004=0,1,0则则2005年径流概率分布为年径流概率分布为 某地区年降水量资料见下表，当某地区年降水量资料见下表，当Q 350mm时，为干旱年；当时，为干旱年；当350mmQ 500mm时，为半干旱年；当时，为半干旱年；当 500mmQ 700mm时，为中时，为中水年；当水年；当 700mmQ 800mm时，为半湿润年；当时，为半湿润年；当Q800mm时，为湿时，为湿润年润年。试估算。试估算2001年降水量为干旱、半干旱、中、半湿润、湿润状态年降水量为干旱、半干旱、中、半湿润、湿润状态的概率。的概率。年份年份年份年份195519551956195619571957195819581959195919601960196119611962196219631963196419641965196519661966降水量降水量降水量降水量661661744.4744.4369369471.1471.1744.8744.8552.1552.11113.51113.5482.7482.7914.1914.1967.2967.2343.2343.2819.3819.3年份年份年份年份196719671968196819691969197019701971197119721972197319731974197419751975197619761977197719781978降水量降水量降水量降水量615.5615.5365.3365.3476.3476.3682.1682.1663.9663.9541541736.3736.3596.4596.4440.6440.6787.5787.5653.8653.8496.7496.7年份年份年份年份197919791980198019811981198219821983198319841984198519851986198619871987198819881989198919901990降水量降水量降水量降水量728.4728.4405.3405.3415.7415.7404.6404.6450.4450.4626.3626.3752.6752.6418.6418.6736.7736.7418.4418.4459.8459.8716.4716.4年份年份年份年份19911991199219921993199319941994199519951996199619971997199819981999199920002000降水量降水量降水量降水量552.6552.6532532405.1405.1837.8837.8579.6579.6593.2593.2504.8504.8686.7686.7314.4314.4532.7532.7作业：作业：