工程水文及水利计算-第8讲

工程水文及水利计算工程水文及水利计算门宝辉门宝辉可再生能源学院可再生能源学院提提 纲纲第一章第一章 绪论绪论第二章第二章 水循环及径流形成水循环及径流形成第三章第三章 水文资料的观测、收集与处理水文资料的观测、收集与处理 第四章第四章 水文统计基本知识水文统计基本知识第五章第五章 设计年径流及径流随机模拟设计年径流及径流随机模拟 （分（分3 3讲）讲）第六章第六章 由流量资料推求洪水由流量资料推求洪水 （分（分3 3讲）讲）第七章第七章 流域产、汇流计算流域产、汇流计算 （分（分4 4讲）讲）第八章第八章 由暴雨资料推求设计洪水由暴雨资料推求设计洪水 （分（分3 3讲）讲）第九章第九章 水文预报水文预报第十章第十章 水库兴利调节计算水库兴利调节计算 （分（分2 2讲）讲）第十一章第十一章 水电站的水能计算水电站的水能计算 （分（分2 2讲）讲）第十二章第十二章 水库防洪计算水库防洪计算 （分（分2 2讲）讲）第五章第五章 目录目录第五章第五章 设计年径流及径流随机模拟设计年径流及径流随机模拟5-1 5-1 概述概述5-2 5-2 影响年径流的因素影响年径流的因素5-3 5-3 具有长期实测径流资料时设计年径流量及其年内分配具有长期实测径流资料时设计年径流量及其年内分配的分析计算的分析计算 5-4 5-4 具有短期实测径流资料时设计年径流量及其年内分配具有短期实测径流资料时设计年径流量及其年内分配的分析计算的分析计算5-5 5-5 缺乏实测径流资料时设计年径流量及其年内分配的分缺乏实测径流资料时设计年径流量及其年内分配的分析计算析计算5-6 5-6 流量历时曲线流量历时曲线5-7 5-7 设计枯水流量的分析计算设计枯水流量的分析计算5-8 5-8 径流随机模拟径流随机模拟5-4 5-4 具具有有短短期期实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内分配的分析计算内分配的分析计算资料条件：资料条件：当实测资料年数当实测资料年数n30，为具有长期实测资；，为具有长期实测资；当当30n12，为具有短期实测资料；，为具有短期实测资料；当当n12，为缺乏实测资料。，为缺乏实测资料。关键：关键：展延年径流系列的长度。展延年径流系列的长度。实质：实质：寻求与设计断面径流有密切关系并有较长观测系列寻求与设计断面径流有密切关系并有较长观测系列的参证变量，通过设计断面年径流与其参证变量的相关关的参证变量，通过设计断面年径流与其参证变量的相关关系，将设计断面年径流系列适当地加以延长至规范要求的系，将设计断面年径流系列适当地加以延长至规范要求的长度。长度。当年径流系列适当延长以后，其频率分析方法与本书第四当年径流系列适当延长以后，其频率分析方法与本书第四章所述完全一样。章所述完全一样。1 1、参证变量的选择、参证变量的选择 最常采用的参证变量有最常采用的参证变量有：设计断面的：设计断面的水位水位、上下游测站或、上下游测站或邻近河流测站的邻近河流测站的径流量径流量、流域的、流域的降水量降水量。参证变量应具备下列条件参证变量应具备下列条件：参证变量与设计断面径流量在成因上有密切关系。参证变量与设计断面径流量在成因上有密切关系。参证变量与设计断面径流量有较多的同步观测资料。参证变量与设计断面径流量有较多的同步观测资料。参证变量的系列较长，并有较好的代表性。参证变量的系列较长，并有较好的代表性。5-4 5-4 具具有有短短期期实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内分配的分析计算内分配的分析计算 利用本站的水位资料延长年径流系列利用本站的水位资料延长年径流系列有些测站开始只观测水位，后有些测站开始只观测水位，后来增加了流量测验。可根据其来增加了流量测验。可根据其水位流量关系，将水位资料水位流量关系，将水位资料转化成径流资料。转化成径流资料。2 2、年径流系列的展延、年径流系列的展延 5-4 5-4 具具有有短短期期实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内分配的分析计算内分配的分析计算 利用上下游站或邻近河流测站径流资料，延长设计断面的径流系列利用上下游站或邻近河流测站径流资料，延长设计断面的径流系列同一河流上下游的水量存在着有机同一河流上下游的水量存在着有机联系，因此，当设计断面上下游不联系，因此，当设计断面上下游不太远处有实测径流资料时，常是很太远处有实测径流资料时，常是很好的参证变量。好的参证变量。同一水文气候区内的邻近河流，当同一水文气候区内的邻近河流，当流域面积与设计流域面积相差不太流域面积与设计流域面积相差不太悬殊时，其径流资料也可试选为参悬殊时，其径流资料也可试选为参证变量。证变量。2 2、年径流系列的展延、年径流系列的展延 5-4 5-4 具具有有短短期期实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内分配的分析计算内分配的分析计算 利用年降水资料延长设计断面的年径流系列利用年降水资料延长设计断面的年径流系列径流是降水的产物。流域的年径流量与流域的年降水量往往有良好的相径流是降水的产物。流域的年径流量与流域的年降水量往往有良好的相关关系。又因降水观测系列在许多情况下较径流观测系列长，因此降水关关系。又因降水观测系列在许多情况下较径流观测系列长，因此降水系列常被用来作为延长径流系列的参证变量。系列常被用来作为延长径流系列的参证变量。2 2、年径流系列的展延、年径流系列的展延 5-4 5-4 具具有有短短期期实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内分配的分析计算内分配的分析计算 注意：注意：（1 1）找出有物理成因的真实关系；）找出有物理成因的真实关系；（2 2）不要通过展转相关建立相关关系；）不要通过展转相关建立相关关系；（3 3）不要过分延长相关线展延资料。）不要过分延长相关线展延资料。2 2、年径流系列的展延、年径流系列的展延 5-4 5-4 具具有有短短期期实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内分配的分析计算内分配的分析计算5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算在部分中小设计流域内，有时只有零星的径流观测资料，且无在部分中小设计流域内，有时只有零星的径流观测资料，且无法延长其系列，甚至完全没有径流观测资料，则只能利用一些法延长其系列，甚至完全没有径流观测资料，则只能利用一些间接的方法，对其设计径流量进行估算。间接的方法，对其设计径流量进行估算。前提：前提：设计流域所在的区域内，有水文特征值的地区综合分析设计流域所在的区域内，有水文特征值的地区综合分析成果，或在水文相似区内有径流系列较长的参证站可资利用。成果，或在水文相似区内有径流系列较长的参证站可资利用。方法方法：采用地理综合法间接推求年径流量的统计参数均值、：采用地理综合法间接推求年径流量的统计参数均值、Cv和和Cs，然后推求满足设计保证率的年径流量，最后采用同倍比，然后推求满足设计保证率的年径流量，最后采用同倍比放大法推求设计年径流过程。放大法推求设计年径流过程。常用方法：常用方法：水文比拟法水文比拟法、参数、参数等值线图法等值线图法和和经验公式法经验公式法。5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算1 1、水文比拟法、水文比拟法 水文比拟法：水文比拟法：在水文相似区内，无资料流域移置在水文相似区内，无资料流域移置(经过修正经过修正)有资料流域实测水文特征值的一种方法，特别适用于年径流有资料流域实测水文特征值的一种方法，特别适用于年径流的分析估算。的分析估算。当设计断面（设计站）缺乏实测径流资料，但其上下游或水当设计断面（设计站）缺乏实测径流资料，但其上下游或水文相似区内有实测水文资料可以选作参证站时，可采用本法文相似区内有实测水文资料可以选作参证站时，可采用本法估算设计年径流。估算设计年径流。要点：要点：将参证站的径流特征值，经过适当的修正后移用于设将参证站的径流特征值，经过适当的修正后移用于设计断面。计断面。5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算修正参变量：修正参变量：流域面积、多年平均降水量。流域面积、多年平均降水量。多年平均年径流量的计算多年平均年径流量的计算均值均值流域面积为主要参变量，二者应比较接近，通常以不超过流域面积为主要参变量，二者应比较接近，通常以不超过1515为宜；如径流的相似性较好，也可以适当放宽上述限制。为宜；如径流的相似性较好，也可以适当放宽上述限制。当设计流域无降水资料时，可以不采用降水为参变量。当设计流域无降水资料时，可以不采用降水为参变量。1 1、水文比拟法、水文比拟法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算 多年平均年径流量的计算多年平均年径流量的计算1 1、水文比拟法、水文比拟法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算 多年平均年径流量的计算多年平均年径流量的计算Cv值值1 1、水文比拟法、水文比拟法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算Cs值值Cv一般可以直接采用，不需要进行修正。一般可以直接采用，不需要进行修正。Cs一般不需要进行计算，取用一般不需要进行计算，取用Cs=23Cv 设计成果的移植设计成果的移植如果参证站已有年径流分析成果，也可以用下列公式，将参如果参证站已有年径流分析成果，也可以用下列公式，将参证站的设计年径流直接移用于设计流域。证站的设计年径流直接移用于设计流域。下标下标P代表频率，其他符号的意义同前。代表频率，其他符号的意义同前。1 1、水文比拟法、水文比拟法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算 设计成果的移植设计成果的移植水文比拟法成果的精度水文比拟法成果的精度，取决于设计流域和参证流域的相似，取决于设计流域和参证流域的相似程度，特别是流域下垫面的情况要比较接近。程度，特别是流域下垫面的情况要比较接近。年内分配：年内分配：直接移植典型年或地区综合的年内分配过程。直接移植典型年或地区综合的年内分配过程。同倍比放大法：同倍比放大法：推求设计年径流分配过程。推求设计年径流分配过程。特点：特点：对于较小流域，由于缺乏实测水文资料，等值线图又对于较小流域，由于缺乏实测水文资料，等值线图又难以应用，水文比拟法具有较大的实用价值。难以应用，水文比拟法具有较大的实用价值。1 1、水文比拟法、水文比拟法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算 原理原理等值线图可以反映水文特征值的地理分布规律，如等值线图可以反映水文特征值的地理分布规律，如中国年降中国年降雨量等值线图雨量等值线图反映的是多年平均年降雨量随地理坐标而发生反映的是多年平均年降雨量随地理坐标而发生的变化。的变化。影响水文特征值的因素分类：影响水文特征值的因素分类：当影响水文特征值（如年径流当影响水文特征值（如年径流深、年降雨量等）的因素主要是深、年降雨量等）的因素主要是区域性因素区域性因素（如气候因素）（如气候因素）时，则该特征值随地理坐标不同而发生连续的变化，利用这时，则该特征值随地理坐标不同而发生连续的变化，利用这种特性就可以在地图上作出它的等值线图；种特性就可以在地图上作出它的等值线图；2 2、等值线图法、等值线图法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算影响水文特征值的因素分类：影响水文特征值的因素分类：当影响水文特征值（如洪峰流量、特征水位等）的因素主要当影响水文特征值（如洪峰流量、特征水位等）的因素主要是是非区域性因素非区域性因素（如流域面积、河床下切深度等）时，则该（如流域面积、河床下切深度等）时，则该特征值不随地理坐标而连续变化，就无法作出它的等值线图；特征值不随地理坐标而连续变化，就无法作出它的等值线图；对于同时受对于同时受区域性和非区域性区域性和非区域性两种因素影响的特征值（如年两种因素影响的特征值（如年径流总量等），在消除非区域性的影响后，才能得出该特征径流总量等），在消除非区域性的影响后，才能得出该特征值的地理分布规律。值的地理分布规律。原理原理2 2、等值线图法、等值线图法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算闭合流域多年水量平衡方程为：闭合流域多年水量平衡方程为：R0=P0-E0从上式可以看出，影响闭合流域多年平均年径流量从上式可以看出，影响闭合流域多年平均年径流量的因素主要是气候因素，即多年平均降雨量和蒸发的因素主要是气候因素，即多年平均降雨量和蒸发量。由于降雨量和蒸发量具有地理分布规律，所以量。由于降雨量和蒸发量具有地理分布规律，所以多年平均年径流量也具有这一规律，因而可以绘制多年平均年径流量也具有这一规律，因而可以绘制多年平均年径流深等值线图。多年平均年径流深等值线图。原理原理2 2、等值线图法、等值线图法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算 多年平均年径流深等值线图多年平均年径流深等值线图绘制绘制为了消除流域面积这一非区域性因素的影响，等值线图总是为了消除流域面积这一非区域性因素的影响，等值线图总是以径流深（毫米）或径流模数（升以径流深（毫米）或径流模数（升秒秒-1-1 公里公里-2-2）来表示。）来表示。绘制气候特征值（如降雨量、蒸发量）等值线图时，是把各绘制气候特征值（如降雨量、蒸发量）等值线图时，是把各观测值点注在地图中各对应观测值点注在地图中各对应观测点观测点的位置上。的位置上。2 2、等值线图法、等值线图法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算 多年平均年径流深等值线图多年平均年径流深等值线图绘制绘制绘点位置：绘点位置：把多年平均年径流深点注在流域形心处。把多年平均年径流深点注在流域形心处。形心：形心：沿不同方向将流域面积大致等分为两等分的几条直线沿不同方向将流域面积大致等分为两等分的几条直线的交点。的交点。方法：方法：将一些有实测资料流域的多年平均年径流深点注在各将一些有实测资料流域的多年平均年径流深点注在各流域形心处后，再考虑降雨及地形特征勾绘等值线，最后用流域形心处后，再考虑降雨及地形特征勾绘等值线，最后用大、中流域的资料加以校核调整，并和多年平均年降雨量等大、中流域的资料加以校核调整，并和多年平均年降雨量等值线图对照，消除不合理现象。值线图对照，消除不合理现象。2 2、等值线图法、等值线图法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算 多年平均年径流深等值线图多年平均年径流深等值线图应用应用通过形心插值：通过形心插值：在小流域中，流域内通过的等值线很少，可在小流域中，流域内通过的等值线很少，可按通过流域形心的直线距离比例内插法，计算流域平均径流按通过流域形心的直线距离比例内插法，计算流域平均径流深。深。2 2、等值线图法、等值线图法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算 多年平均年径流深等值线图多年平均年径流深等值线图应用应用首先绘出通过流域长短轴的首先绘出通过流域长短轴的两条垂线，其交点两条垂线，其交点O即近似即近似地位于流域的形心位置。通地位于流域的形心位置。通过过O点作大致垂直于点作大致垂直于700、800mm两条等值线的直线，两条等值线的直线，交点分别为交点分别为A、B，则流域，则流域的平均径流深为的平均径流深为2 2、等值线图法、等值线图法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算 多年平均年径流深等值线图多年平均年径流深等值线图应用应用面积加权法：面积加权法：如果设计流域内通过多如果设计流域内通过多条年径流深等值线，可条年径流深等值线，可以用面积加权法推求流以用面积加权法推求流域的平均径流深。域的平均径流深。2 2、等值线图法、等值线图法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算 多年平均年径流深等值线图多年平均年径流深等值线图精度精度大流域：大流域：一般具有长期的实测资料，很少用等值线图；一般具有长期的实测资料，很少用等值线图；中流域：中流域：精度一般较高，具有很大的实用价值；精度一般较高，具有很大的实用价值；小流域：小流域：误差很大，且由等值线图查得的径流深有可能偏误差很大，且由等值线图查得的径流深有可能偏大。因为绘制等值线图时主要依据的是中等流域的资料，大。因为绘制等值线图时主要依据的是中等流域的资料，由等值线图查得的径流深应该是小流域的径流量全部流出由等值线图查得的径流深应该是小流域的径流量全部流出时的径流量，但小流域不闭合，河槽下切不深，不能全部时的径流量，但小流域不闭合，河槽下切不深，不能全部汇集本流域所形成的地下径流，即实际资料要小于这个数。汇集本流域所形成的地下径流，即实际资料要小于这个数。2 2、等值线图法、等值线图法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算 年径流年径流Cv值的估算值的估算年径流深的年径流深的Cv值，也有等值线图可供查算，方法与年径流均值，也有等值线图可供查算，方法与年径流均值估算方法类似，但可更简单一点，即按比例内插出流域重值估算方法类似，但可更简单一点，即按比例内插出流域重心的心的Cv值就可以了。值就可以了。大流域：大流域：一般具有长期的实测资料，很少用等值线图；一般具有长期的实测资料，很少用等值线图；中流域：中流域：精度一般较高，具有很大的实用价值；精度一般较高，具有很大的实用价值；小流域：小流域：误差很大，且由等值线图查得的误差很大，且由等值线图查得的Cv值有可能偏小。值有可能偏小。因为绘制等值线图时主要依据的是中等流域的资料，一般来因为绘制等值线图时主要依据的是中等流域的资料，一般来说流域越大，径流变化越均匀，说流域越大，径流变化越均匀，Cv值越小。值越小。2 2、等值线图法、等值线图法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算 年径流年径流Cs值的估算值的估算年径流的年径流的Cs值，一般采用值，一般采用Cv的倍比。按照规范规定，一般可的倍比。按照规范规定，一般可采用采用Cs2 23 3Cv。2 2、等值线图法、等值线图法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算年径流的地区综合年径流的地区综合，常以经验公式表示。这类公式主要是与，常以经验公式表示。这类公式主要是与年径流影响因素建立关系。例如，年径流均值的经验公式有年径流影响因素建立关系。例如，年径流均值的经验公式有如下类型：如下类型：b1，b2、n1、n2、m参数，通过实测资料分析确定，或按参数，通过实测资料分析确定，或按已有分析成果采用。已有分析成果采用。Cv、Cs也有地区综合的成果。也有地区综合的成果。3 3、经验公式法、经验公式法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算不同设计频率的年平均流量不同设计频率的年平均流量Qp，也可以建立类似的关系，只，也可以建立类似的关系，只是其参数的定量亦各有不同。是其参数的定量亦各有不同。这类方法的精度，一般较这类方法的精度，一般较等值线图法等值线图法低。但在进行低。但在进行流域初步流域初步规划规划，需要快速估算流域的地表水资源量及水力蕴藏量时，需要快速估算流域的地表水资源量及水力蕴藏量时，有实用价值有实用价值。3 3、经验公式法、经验公式法 5-5 5-5 缺缺乏乏实实测测径径流流资资料料时时设设计计年年径径流流量量及及其其年年内内分分配的分析计算配的分析计算径流的分配过程除用流量过程表示之外，还可以用径流的分配过程除用流量过程表示之外，还可以用流量历时流量历时曲线曲线来表示。来表示。这种曲线是按其时段内这种曲线是按其时段内所出现的流量数值及其所出现的流量数值及其历时历时(或相对历时或相对历时)而绘而绘成的，成的，说明径流分配的说明径流分配的一种特性曲线。一种特性曲线。5-6 5-6 流量历时曲线流量历时曲线多年综合日流量历时曲线是根据所有各年份的实测日平均流多年综合日流量历时曲线是根据所有各年份的实测日平均流量资料绘成的，它能反映流量在多年期间的历时情况量资料绘成的，它能反映流量在多年期间的历时情况。在工程设计中，有时要求绘制在工程设计中，有时要求绘制丰水年丰水年(或枯水年或枯水年)的综合日流的综合日流量历时曲线，它是根据各丰水年量历时曲线，它是根据各丰水年(或枯水年或枯水年)的实测日平均流的实测日平均流量资料绘成的。量资料绘成的。1 1、多年综合日流量历时曲线多年综合日流量历时曲线 所谓所谓丰水期丰水期(枯长期、灌溉期枯长期、灌溉期)的综合日流量历时曲线，它是的综合日流量历时曲线，它是根据所有各年丰水期根据所有各年丰水期(枯水期，灌溉期枯水期，灌溉期)的实测日平均流量资的实测日平均流量资料绘成的。料绘成的。5-6 5-6 流量历时曲线流量历时曲线代表年日流量历时曲线是根据某一年份的实测日平均流量资代表年日流量历时曲线是根据某一年份的实测日平均流量资料绘成的，曲线的纵坐标为日平均流量或其相对值料绘成的，曲线的纵坐标为日平均流量或其相对值(模比系数模比系数)，横坐标则为历时日数或相对历时，横坐标则为历时日数或相对历时(占全年的百分数占全年的百分数)。2 2、代表年日流量历时曲线代表年日流量历时曲线 平均历时曲线是以各年同历时的日平均流量的平均值为纵坐平均历时曲线是以各年同历时的日平均流量的平均值为纵坐标，其相应历时为横坐标点绘的曲线。标，其相应历时为横坐标点绘的曲线。3 3、平均日流量历时曲线、平均日流量历时曲线 平均历时曲线是一种平均历时曲线是一种虚拟的曲线虚拟的曲线。与综合历时曲线相比，它。与综合历时曲线相比，它的上部较低而下部较高，中间则大致与综合曲线重合。的上部较低而下部较高，中间则大致与综合曲线重合。5-6 5-6 流量历时曲线流量历时曲线5-7 5-7 设计枯水流量的分析计算设计枯水流量的分析计算枯水流量枯水流量亦称亦称最小流量最小流量，是河川径流的一种特殊形态。枯，是河川径流的一种特殊形态。枯水流量往往制约着城市的发展规模、灌溉面积、通航的容水流量往往制约着城市的发展规模、灌溉面积、通航的容量和时间，同时，也是决定量和时间，同时，也是决定水电站保证出力水电站保证出力的重要因素。的重要因素。按设计时段的长短按设计时段的长短，枯水流量又可分为瞬时、日、旬、月、，枯水流量又可分为瞬时、日、旬、月、最小流量，其中又以日、旬、月最小流量对水资源利最小流量，其中又以日、旬、月最小流量对水资源利用工程的规划设计关系最大。用工程的规划设计关系最大。时段枯水流量与时段径流量在分析方法上没有本质区别，时段枯水流量与时段径流量在分析方法上没有本质区别，主要在选样方法上有所不同。时段径流在时序上往往是主要在选样方法上有所不同。时段径流在时序上往往是固固定定的，而枯水流量则在一年中选其最小值，在时序上是的，而枯水流量则在一年中选其最小值，在时序上是变变动动的。此外，在一些具体环节上也有一些差异。的。此外，在一些具体环节上也有一些差异。当设计代表站有长系列实测径流资料时，可按年最小选样当设计代表站有长系列实测径流资料时，可按年最小选样原则选取一年中最小的时段径流量，组成样本系列。原则选取一年中最小的时段径流量，组成样本系列。枯水流量常采用不足概率枯水流量常采用不足概率q，即以，即以小于和等于该径流的概率来表示，小于和等于该径流的概率来表示，它和年最大选样的概率它和年最大选样的概率P有有ql一一P的关系。因此在系列排队时按由小的关系。因此在系列排队时按由小到大排列。除此之外，年枯水流量到大排列。除此之外，年枯水流量频率曲线的绘制与时段径流频率曲频率曲线的绘制与时段径流频率曲线的绘制基本相同，也常采用线的绘制基本相同，也常采用P-型频率曲线适线。右图为某水型频率曲线适线。右图为某水文站不同天数的枯水流量频率曲线文站不同天数的枯水流量频率曲线的示例。的示例。1 1、有实测水文资料时的枯水流量计算有实测水文资料时的枯水流量计算5-7 5-7 设计枯水流量的分析计算设计枯水流量的分析计算年枯水流量频率曲线，在某些河流上，特别是在干旱半干年枯水流量频率曲线，在某些河流上，特别是在干旱半干旱地区的中小河流上，还会出现旱地区的中小河流上，还会出现时段径流量为零时段径流量为零的现象，的现象，可用下面一种简易的实用方法来可用下面一种简易的实用方法来处理处理。设系列的全部项数为设系列的全部项数为n，其中非零项数为，其中非零项数为k，零值项数为，零值项数为nk。首先把。首先把k项非零资料视作一个独立系列，按一般方法求项非零资料视作一个独立系列，按一般方法求出其频率曲线。然后通过下列转换，即可求得全部系列的出其频率曲线。然后通过下列转换，即可求得全部系列的频率曲线。其转换关系为：频率曲线。其转换关系为：1 1、有实测水文资料时的枯水流量计算有实测水文资料时的枯水流量计算5-7 5-7 设计枯水流量的分析计算设计枯水流量的分析计算在枯水流量频率曲线上，往往会出现在两端接近在枯水流量频率曲线上，往往会出现在两端接近P20和和P9090处曲线转折现象。在处曲线转折现象。在P=20%以下的部分是河网及潜以下的部分是河网及潜水逐渐枯竭，径流主要靠深层地下水补给。在水逐渐枯竭，径流主要靠深层地下水补给。在P=90%，可，可能是某些年份有地表水补给，枯水流量偏大所致。能是某些年份有地表水补给，枯水流量偏大所致。1 1、有实测水文资料时的枯水流量计算有实测水文资料时的枯水流量计算5-7 5-7 设计枯水流量的分析计算设计枯水流量的分析计算当设计断面短缺径流资料时，设计枯水流量主要借助于参当设计断面短缺径流资料时，设计枯水流量主要借助于参证站延长系列或成果移置，与本章第五节所述方法基本相证站延长系列或成果移置，与本章第五节所述方法基本相同。同。但枯水流量较之固定时段的径流，其时程变化更为稳定。但枯水流量较之固定时段的径流，其时程变化更为稳定。因此，在与参证站建立径流相关时，效果会好一些。或者因此，在与参证站建立径流相关时，效果会好一些。或者说，条件可以适当放宽。说，条件可以适当放宽。例如，当设计站只有少数几年资料，与参证站的相似性较例如，当设计站只有少数几年资料，与参证站的相似性较好时，也可建立较好的枯水流量相关关系。在这种情况下，好时，也可建立较好的枯水流量相关关系。在这种情况下，甚至可以不进行设计站的径流系列延长和频率分析，而直甚至可以不进行设计站的径流系列延长和频率分析，而直接移用参证站的频率分析成果，经上述相关关系，转化为接移用参证站的频率分析成果，经上述相关关系，转化为本站的相应频率的设计枯水流量。本站的相应频率的设计枯水流量。2 2、短缺水文资料时的枯水流量估算短缺水文资料时的枯水流量估算5-7 5-7 设计枯水流量的分析计算设计枯水流量的分析计算在设计站完全没有径流资料的情况下还可以临时进行资料在设计站完全没有径流资料的情况下还可以临时进行资料的补充收集工作，以应需要。如果能施测一个枯水季的流的补充收集工作，以应需要。如果能施测一个枯水季的流量过程，对于建立量过程，对于建立30天以下时段的枯水流量关系，有很大天以下时段的枯水流量关系，有很大用处；如只研究日最小流量，那么在枯水期只施测几次流用处；如只研究日最小流量，那么在枯水期只施测几次流量量(如如10次流量次流量)，就可与参证站径流建立相关关系。，就可与参证站径流建立相关关系。2 2、短缺水文资料时的枯水流量估算短缺水文资料时的枯水流量估算5-7 5-7 设计枯水流量的分析计算设计枯水流量的分析计算5-8 5-8 径流随机模拟径流随机模拟由于实际水文资料往往比较短，难于满足实际水文计算和由于实际水文资料往往比较短，难于满足实际水文计算和水利计算的要求，故需要对径流过程作随机模拟。这种随水利计算的要求，故需要对径流过程作随机模拟。这种随机模拟的目的之一在于机模拟的目的之一在于充分利用现有水文实测资料的信息充分利用现有水文实测资料的信息，而不是而不是用来延长资料长度用来延长资料长度。1 1、随机过程基本知识随机过程基本知识 随机过程和时间序列的定义随机过程和时间序列的定义 随机过程：随机过程：在实际问题中，常涉及试验过程随某个参变量在实际问题中，常涉及试验过程随某个参变量的变化而变化的随机变量，如河流某处的流量、水位是随的变化而变化的随机变量，如河流某处的流量、水位是随时间变化的随机变量，气温是随时间、高度而变的随机变时间变化的随机变量，气温是随时间、高度而变的随机变量等。通常称这种随机变量为随机函数，并称以时间为参量等。通常称这种随机变量为随机函数，并称以时间为参数的数的随机函数随机函数为为随机过程随机过程，记为，记为（t），），tT，T是是t变化变化的范围。的范围。1 1、随机过程基本知识随机过程基本知识 随机过程和时间序列的定义随机过程和时间序列的定义 时间序列：时间序列：若时间参变量若时间参变量T是连续时刻的集合，则称这种随是连续时刻的集合，则称这种随机过程为连续参数随机过程，如水位过程、流量过程等。机过程为连续参数随机过程，如水位过程、流量过程等。若时间参变量若时间参变量T T是离散时刻的集合，则称这种随机过离散参是离散时刻的集合，则称这种随机过离散参数随机过程，也称为数随机过程，也称为随机序列随机序列或或时间序列时间序列。如年、月径流。如年、月径流过程，年最大流量过程都是时间序列，也称水文时间序列。过程，年最大流量过程都是时间序列，也称水文时间序列。5-8 5-8 径流随机模拟径流随机模拟 随机过程的数字特征随机过程的数字特征 数学期望数学期望 方差方差 自相关函数自相关函数 随机过程的基本分类随机过程的基本分类 按统计性质的稳定性分类按统计性质的稳定性分类：若随机过程统计数字特征不随时间若随机过程统计数字特征不随时间的平移而变化，则称为的平移而变化，则称为平稳过程平稳过程，否则为，否则为非平稳过程非平稳过程。按不同时刻状态间的关系分类按不同时刻状态间的关系分类：可分成独立过程和马尔柯夫过可分成独立过程和马尔柯夫过程。若过程各状态相互独立，则称为程。若过程各状态相互独立，则称为独立随机过程独立随机过程。在非独立。在非独立随机过程中，最重要的一类是随机过程中，最重要的一类是马尔柯夫过程马尔柯夫过程，其特点是，其特点是tn时刻时刻状态只与状态只与tn-1时刻有关，而与时刻有关，而与tn-1以前各时刻无关。以前各时刻无关。1 1、随机过程基本知识随机过程基本知识5-8 5-8 径流随机模拟径流随机模拟时间序列组成分析；时间序列组成分析；模型的建立；模型的建立；序列的生成；序列的生成；模型及生成系列的检验。模型及生成系列的检验。2 2、径流、径流随机模拟的一般步骤随机模拟的一般步骤5-8 5-8 径流随机模拟径流随机模拟水文序列水文序列Qt一般可以表示成下式：一般可以表示成下式：Tt，Ct，Pt，St分别为趋势项、跳跃项、周期项和随机项。分别为趋势项、跳跃项、周期项和随机项。3 3、水文时间序列的组成分析、水文时间序列的组成分析5-8 5-8 径流随机模拟径流随机模拟对年月径流序列建立随机模型一般是对原始年月径流序列排除对年月径流序列建立随机模型一般是对原始年月径流序列排除趋势、跳跃等确定性成分后的随机成分而建立的模型。趋势、跳跃等确定性成分后的随机成分而建立的模型。设用于建模的年径流序列为：设用于建模的年径流序列为：Q1、Q2、Qn目前常见的随机模型有：目前常见的随机模型有：线性平稳模型线性平稳模型、非线性平稳模型非线性平稳模型，马马尔柯夫模型尔柯夫模型及及非马尔柯夫模型等非马尔柯夫模型等。对年径流序列而言，常采用对年径流序列而言，常采用线性平稳模型线性平稳模型或或马尔柯夫模型马尔柯夫模型。建立随机模型的一般步骤为：建立随机模型的一般步骤为：选择模型；选择模型；确定阶数；确定阶数；估估计模型参数。计模型参数。4 4、单站年、月径流随机模型的建立、单站年、月径流随机模型的建立5-8 5-8 径流随机模拟径流随机模拟模型的选择：模型的选择：通常采用线性自回归模型。通常采用线性自回归模型。线性自回归模型的一船形式线性自回归模型的一船形式：令令 单站年径流随机模型的建立单站年径流随机模型的建立 4 4、单站年、月径流随机模型的建立、单站年、月径流随机模型的建立5-8 5-8 径流随机模拟径流随机模拟模型参数的估计：模型参数的估计：单站年径流随机模型的建立单站年径流随机模型的建立 均值均值方差方差 自回归系数自回归系数 4 4、单站年、月径流随机模型的建立、单站年、月径流随机模型的建立5-8 5-8 径流随机模拟径流随机模拟模型参数的估计：模型参数的估计：单站年径流随机模型的建立单站年径流随机模型的建立 在在n较大，较大，K较小时较小时 ，K阶样本自相关系数阶样本自相关系数 4 4、单站年、月径流随机模型的建立、单站年、月径流随机模型的建立5-8 5-8 径流随机模拟径流随机模拟常见常见AR（1 1）及）及AR（2 2）模型参数估计公式）模型参数估计公式 单站年径流随机模型的建立单站年径流随机模型的建立 AR（1 1）模型形式为：）模型形式为：模型参数估计公式为模型参数估计公式为：4 4、单站年、月径流随机模型的建立、单站年、月径流随机模型的建立5-8 5-8 径流随机模拟径流随机模拟常见常见AR（1 1）及）及AR（2 2）模型参数估计公式）模型参数估计公式 单站年径流随机模型的建立单站年径流随机模型的建立 AR（2 2）模型形式为：）模型形式为：模型参数估计公式为模型参数估计公式为：4 4、单站年、月径流随机模型的建立、单站年、月径流随机模型的建立5-8 5-8 径流随机模拟径流随机模拟模型阶数模型阶数P P的确定的确定 单站年径流随机模型的建立单站年径流随机模型的建立 4 4、单站年、月径流随机模型的建立、单站年、月径流随机模型的建立5-8 5-8 径流随机模拟径流随机模拟统计推断方法是：统计推断方法是：取显著水平取显著水平 ，若若 ，则可接受，则可接受 异于异于0 0的的假设。假设。对已有对已有n年实测月径流资料的单站通常有两条途径建立其月径流年实测月径流资料的单站通常有两条途径建立其月径流随机模型。一是先建立年径流模型，再通过建立随机模型。一是先建立年径流模型，再通过建立解集模型解集模型把年径把年径流分解成各月月径流；二是直接建立月径流随机模型，通常采用流分解成各月月径流；二是直接建立月径流随机模型，通常采用季节性一阶自回归模型，即假定可用季节性一阶自回归模型，即假定可用1212个一阶自回归模型来描述个一阶自回归模型来描述各月月径流流量及相关系数。各月月径流流量及相关系数。解集模型解集模型 单站年径流随机模型的建立单站年径流随机模型的建立 4 4、单站年、月径流随机模型的建立、单站年、月径流随机模型的建立5-8 5-8 径流随机模拟径流随机模拟各月月径流模型各月月径流模型 单站年径流随机模型的建立单站年径流随机模型的建立 4 4、单站年、月径流随机模型的建立、单站年、月径流随机模型的建立5-8 5-8 径流随机模拟径流随机模拟 纯随机变量随机数的生成纯随机变量随机数的生成 年、月径流序列模拟涉及到纯随机变量年、月径流序列模拟涉及到纯随机变量t的生成及由的生成及由t生成年、生成年、月径流序列的问题月径流序列的问题。纯随机变量纯随机变量t t的分布可以是正态，也可以是偏态。它的生成，的分布可以是正态，也可以是偏态。它的生成，一般先生成一般先生成00，11均匀分布随机数均匀分布随机数u u，再通过变换生成指定分，再通过变换生成指定分布的随机数。布的随机数。5 5、年、月径流序列的生成（模拟）、年、月径流序列的生成（模拟）5-8 5-8 径流随机模拟径流随机模拟 纯随机变量随机数的生成纯随机变量随机数的生成 生成方法有生成方法有随机数表法随机数表法、物理方法物理方法及及数学方法数学方法。由于前两种方。由于前两种方法存在严重缺陷，故常用法存在严重缺陷，故常用数学方法数学方法生成，其中应用最广的是生成，其中应用最广的是乘乘同余法同余法。均匀分布随机数的生成均匀分布随机数的生成 乘同余法生成随机数递推公式是乘同余法生成随机数递推公式是:5 5、年、月径流序列的生成（模拟）、年、月径流序列的生成（模拟）5-8 5-8 径流随机模拟径流随机模拟这种方法生成的随机数存在着循环周期，因此，这种方法生成的随机数存在着循环周期，因此，un+1 不是真正不是真正意义上的随机数，俗称意义上的随机数，俗称“伪随机数伪随机数”。但由于。但由于M往往取值很大，往往取值很大，周期也很长，目前微机上周期可达周期也很长，目前微机上周期可达109以上，实用上完全能满以上，实用上完全能满足需要。正因如此，实际大都使用该法生成足需要。正因如此，实际大都使用该法生成0，1均匀分布随均匀分布随机数。不过使用前要对生成的随机数作均匀性、独立性等检验。机数。不过使用前要对生成的随机数作均匀性、独立性等检验。纯随机变量随机数的生成纯随机变量随机数的生成 均匀分布随机数的生成均匀分布随机数的生成 5 5、年、月径流序列的生成（模拟）、年、月径流序列的生成（模拟）5-8 5-8 径流随机模拟径流随机模拟通常用通常用BoxMullerBoxMuller变换生成变换生成 纯随机变量随机数的生成纯随机变量随机数的生成 正态分布随机数的生成正态分布随机数的生成 5 5、年、月径流序列的生成（模拟）、年、月径流序列的生成（模拟）5-8 5-8 径流随机模拟径流随机模拟利用舍选法生成利用舍选法生成P型分布随机数型分布随机数 纯随机变量随机数的生成纯随机变量随机数的生成 PP型分布随机数的生成型分布随机数的生成 5 5、年、月径流序列的生成（模拟）、年、月径流序列的生成（模拟）5-8 5-8 径流随机模拟径流随机模拟设所建立的模型为设所建立的模型为AR(1)(1)，且参数已估计出。其模型为，且参数已估计出。其模型为：考虑年径流为正态分布考虑年径流为正态分布 年、月径流的生成年、月径流的生成 5 5、年、月径流序列的生成（模拟）、年、月径流序列的生成（模拟）5-8 5-8 径流随机模拟径流随机模拟考虑年径流为正态分布考虑年径流为正态分布 年、月径流的生成年、月径流的生成 5 5、年、月径流序列的生成（模拟）、年、月径流序列的生成（模拟）5-8 5-8 径流随机模拟径流随机模拟考虑年径流为偏态分布考虑年径流为偏态分布 年、月径流的生成年、月径流的生成 5 5、年、月径流序列的生成（模拟）、年、月径流序列的生成（模拟）5-8 5-8 径流随机模拟径流随机模拟残差独立性检验残差独立性检验 模型检验是指所选用的模型残差模型检验是指所选用的模型残差t为独立随机变量是否成立，为独立随机变量是否成立，t分布是否为假定分布。生成序列检验是指所生成年月径流序分布是否为假定分布。生成序列检验是指所生成年月径流序列是否能反映实测径流的统计特性。列是否能反映实测径流的统计特性。生成系列检验生成系列检验 一般要求生成序列与实测序列统计特征相近。一般要求生成序列与实测序列统计特征相近。6 6、模型及生成系列的检验、模型及生成系列的检验5-8 5-8 径流随机模拟径流随机模拟小小 结结1 1、具有短期实测径流资料时设计年径流量及其年内、具有短期实测径流资料时设计年径流量及其年内分配的分析计算分配的分析计算2 2、缺乏实测径流资料时设计年径流量及其年内分配、缺乏实测径流资料时设计年径流量及其年内分配的分析计算的分析计算3 3、流量历时曲线、流量历时曲线4 4、设计枯水流量的分析计算、设计枯水流量的分析计算5 5、径流随机模拟、径流随机模拟