实验一径流的随机分析法.docx

实验一：径流的随机分析法（古尔德法）（一）、古尔德法的原理古尔德法又称直接求总库容的概率演算法，它以实测年、月径流资料为基础，假设入库年径流过程为独立随机序列，以年初蓄水状态为条件，对各年进行年内各月的调节演算求年末的蓄水状态，然后再根据实际演算结果，统计出蓄水状态的年转移概率，并沿用马氏过程的原理分析水库蓄水状态的年转移概率矩阵，并推求出水库稳定的蓄水概率分布古尔德法的特点是它不仅考虑了水库对径流的年际调节，而且同时考虑了径流的年内调节二）、古尔德法步骤1、将水文年径流资料转为水利年径流资料：结果如表“径流资料”2、状态离散：将水库的兴利库容Vc从库空（0）到库满（Vn）,其间等分为（m—2）份，（m由精度等而定，本次取为10）,每份为DV=Vn/（m—2）,各状态序号i所对应的状态值A（i）及状态的数值域x（i）按下式计算：i=l时A（i）=0x（1）W01ViVm时A（i）=DVX[i—（3/2）]（i-2）DVvx（i）v（i-1）DVi=m时A（i）=Vnx（m）三Vn3、推求转移概率矩阵和条件破坏概率以一个指定的年初蓄水状态A（i）为出发点，根据某年实测径流资料和供水资料，采用简化运行策略进行年内调节计算，求出水库各月的蓄水过程及其年末的蓄水状态。

记K表示年份序号，L表示月份的序号，S（K，L）表示水库的第K年L月份的月初蓄水量，X（K，L）表示第K年L月份的月径流，d（L）表示第L月份的调节流量，则初始条件为S（K,1）=A（i）逐月径流调节计算时，实际上是执行下式的水量平衡演算，即S（K，L+1）=S（K，L）+X（K，L）—d（L）演算中规定：水库的蓄水量应满足下式的约束，即0WS（K，L）WVn当S（K，L+1）>Vn时，应取S（K，L+1）=Vn，其多余水量为S（K，L）+x（K，L）—d（L）—Vn,作为弃水泄放当S（K，L+1）V0时，应取S（K，L+1）=0,水库只能按S（K，L）+x（K，L）供水，由于S（K，L）+x（K,L）Vd（L），因此这个月将出现正常供水的破坏，规定若一年中任何一个月出现正常供水的不足则该年计为不能保证正常供水的年份即这一年供水破坏年末蓄水状态S（K，13）视其所在的状态数值域判断其状态序号jo按上述演算方法，可对实测各水利年每年都以年初S（K，L）=A（i）为初始蓄水状态进行各年的演算，求出各年的水库蓄水过程及其相应的年末蓄水状态S（K，13）=A（j）o然后统计出各年初蓄水状态序号i下各年末蓄水状态序号j的出现年数，并据此计算出蓄水状态的年转移概率q”，以及以年初蓄水状态为条件的供水破坏概率FY、FM。

蓄水状态的年转移概率可按下式计算qij年初状态序号为i出现年末状态序号为的年数总计算年数年初蓄水状态序号为i时出现供水破坏年数条件供水破坏概率按下式计算：总年数年初蓄水状态序号为i时出现供水破坏月数总月数对于年初蓄水状态X(i),i=l,2,3,...10,逐一进行相同的分析，最后可求得蓄水状态的年转移概率矩阵Q={qij}i,j=l，2,…，10i,j及各状态点的条件供水破坏概率：FY=[fy(1)......fy(10)]TFM=[fm(1)fm(10)]T4、推求水库稳定的蓄水概率分布(用迭代法进行演算)假定第一年初的蓄水状态的概率为：n1=[1000000000]利用蓄水概率分布演算式nt+1=ntq进行逐年迭代演算，求出第一年末的蓄水概率分布n2=n&，以第一年的年末的蓄水状态的概率分布n2作为第二年的初蓄，求第二年末的蓄水概率分布n3,直到此概率分布收效于一个稳定的蓄水概率分布n，此时n=nt+1^nt5、水库正常供水保证率计算前面在推求水库蓄水状态的年转移概率的同时，也求出了条件的洪水破坏概率FY、FM，这些条件的供水破坏概率只与年初蓄水状态有关而水库的稳定的蓄水状态的概率分布代表水库正常运行条件下的年初蓄水概率分布，根据概率论原理可知，年初蓄水出现状态X(i)与该年出现正常供水破坏这二个事件同时出现的概率等于二者各自的出现概率的乘积,故可按全概率公式求水库长期运行的正常供水破坏概率。

按相对破坏年数表示，供水破坏概率的计算式为：HY=FY二工兀(i)f(i)yi按相对破坏月数表示，供水破坏概率：HM=□•FM二工兀(i)f(i)mi分别求出相对年数表示的正常供水保证率PYPY=1－HY和相对月数表示的正常供水保证率PMPM=1－HM三)、简化运行策略水库的运行策略是指水库在已知来水的情况下，如何放水的问题，而简化运行策略则是运行策略的一种概化策略，它是一种将水库的入流量作为一个确定性的径流过程的简化的计算方法，适用于水资源系统规划阶段的水库模拟运行(或长系列操作)，它与径流调节中的等流量调节时历法相似图3.1为水库简化运行策略的示意图图中abed为水库供水指示线，Vn为有效库容，M为图3.l水库简化运行策略示意图模拟实验方法将时段入流作为已知，死库容取0设时段K入库水量为X(K),时段初水库蓄水量为S(K),时段末水库蓄水量为S(K+1),水库的实际时段供水量为W(K)以如下规则按简化运行策略进行操作演算第一条规则：S(K)+X(K)VM,W(K)=S(K)+X(K),S(K+1)=0；第二条规则：MWS(K)+X(K)WM+Vn,W(K)=M,S(K+1)=S(K)+X(K)—M；第三条规则：S(K)十X(K)>M+Vn,W(K)=S(K)+X(K)-Vn,S(K+1)=Vn；第一条规则与指示线ab段相对应，即当时段初蓄水量S(K)与该时段入库水量X(K)之和小于该时段要求的额定供水量M时，水库只能按时段初蓄水量与时段入库水量之和进行供水，控制时段末蓄水量不致于消落至死水位(死库容)以下，因此ab段反映水库限制供水执行,执行结果是该时段不能满足正常供水要求(即正常供水遭到破坏)。

第二规则与指示线bc段相对应，即当时段初蓄水量S(K)与该时段入库水量X(K)之和大于该时段要求的额定供水量M时，水库在满足正常供水量的同时，尚有多余水量可充蓄水库第三条规则与指示线cd段相对应，即当时段初蓄水量S(K)与该时段入库水量X(K)之和除满足正常供水量M和充蓄有效库容Vn之外，仍有余水，可加大供水，但控制水库供水量不超过水库引水设备的最大过水能力，超过则通过泄洪建筑物泄放(弃水)，使水库水位维持在正常蓄水位上四)、程序/*x(i,j)--水库入流(年，月径流量)A(i)--蓄水状态变量S(i,j)--月初蓄水量S1(i)--年末蓄水量Q(i,j)--年转移概率矩阵IND(j)--年末出现小于等于某一状态的总次数PI(i)--水库稳定蓄水概率分布FY(i)--条件供水破坏概率(相对破坏年数)FM(i)--条件供水破坏概率(相对破坏月数)QP--调节流量(已转化为水量)Vn--有效库容M--状态总数K1--计算总年教PY--年供水保证率PM--月供水保证率*/#include#include#defineM60#defineK149#definem18voidmain(){/*状态离散*/inti,j,k,l,e;floatZ[M],V[M],A[m];floatZ0=162.0,Zn=207.5;floatV0,Vn,DV;ifstreamZV("H:\\水位库容关系.txt");for(i=0;i>Z[i]>>V[i];ZV.close();for(i=0;i=Z[i]&&Z0=Z[i]&&Zn>x[i][j];Qin.close();for(i=0;i=Vn)S[k][l+1]=Vn;}if(INF>0){IY=IY+1;IM=IM+INF;}S1[k]=S[k][12];}FY[i]=IY/float(K1);FM[i]=IM/float(K1*12);/*推求水库蓄水状态的年转移概率*/for(j=0;j=0.001)PAI[j]=PI[j];}/*计算水库正常供水保证率*/floatHY=0.0,HM=0.0;floatPY,PM;for(i=0;i