水信息技术课程设计.docx

水信息技术课程设计 水文预报课程设计 学院：水利土木工程学院 专业：水文与水资源管理 班级：* 学号：* 姓名：* 指导老师：程* 时间： 2022年12月20日 第一章设计要求 1、根据所给的资料，分析稳定的水位流量关系的特点，并用单一线法定线（单一线的拟合采用正交函数法通过编程计算完成。正交函数法参考课本P269）； 2、对所定的水位流量关系曲线进行符号检验、适线检验和偏离数值检验； 3、推求6月5日至6月11日逐时流量，并计算日平均流量。 第二章设计资料 1、测站1976年实测流量成果表见表1（说明：测流断面为基本测流断面， 测流方法为流速仪测流）； 2、测站1976年6月5日至6月12日部分洪水水文要素摘录表见表2； 3、测站水文概况 测验河段顺直，河床整齐多卵石，土质坚实，两岸无显著冲淤变化，水位超过10米时，两岸发生漫滩。左岸漫滩约700米，右岸漫滩约1200米。测站上游约五公里处为两个支流的汇合口。测站下游800米处急滩上有一水坝，下游15 公里处有一较大支流汇入。 第三章水位流量关系特点 2.1 水位流量关系简介 水位流量关系stage-discharge relation 河渠中某断面的流量与其水位之间的对应状态。这项关系要根据当地的不同水流情况的多次实测流量和其相应的水位来确定，通常用经验曲线、经验方程或表格等形式表达。它是由各种水力因素（如水面宽、断面面积、水力比降、糙率等）决定的。 2.2 稳定水位流量关系 由上面表1数据作图: 从上图可知，在较长时期内，该断面的实测流量与相应水位的点据呈密集带状分布，可用一条单一曲线来表示。这是稳定的水位流量关系。水位流量关系维持稳定，必须具备下列条件之一：断面面积、水力比降和糙率等水力因素在同一水位时,维持不变;在同一水位时，上述各因素虽有变动，但其变动对水位流量关系的影响可以互相补偿。在这种条件下，同一个水位，就只有一个相应的流量，水位流量关系就成为单一的曲线。水文测站断面总是尽可能选择在那些具备基本稳定条件的地方，必要时，用人工铺设固定河床或在断面下游附近建造拦河 低堰等方法来创造这种稳定条件。应用稳定的水位流量关系对于编制逐日平均流量、洪水水文要素摘录表和因工程设计或防汛需要推求各种水位或流量都十分方便。 第四章单一曲线定线 4.1 单一曲线定线原理 应根据定线的原则，用目估的方法，通过点群中心，分别绘出水位流量、水位面积、水位流速单一的光滑曲线，测点应均匀地分面于关系曲线的两旁，不允许测点有系统偏离现象。对初步制定的关系曲线，应检查调正。检查方法是：在三关系曲线的各级水位上，分别查出流量、面积、流速；进行面积乘流速是否等于流量的计算，若误差不超过±2±3%时，即认为合理，否则，应调正曲线，直到符合为止。 4.2 单一曲线定线方法 1、图解法-人工定线 2、解析法-人工或计算机选配方程定线 4.3 单一曲线定线步骤 1、在同一张图纸上依次点绘Z-Q 、Z-A、Z-V关系曲线。 (1)三条曲线比例尺选择，应使它们与横轴的夹角近似为45?、60?、60?，且互不相交。 (2)测点右侧注明测次号 (3)3条曲线左右顺序为Z-Q 、Z-A 、Z-V。 2、分析Z-Q 、Z-A 、Z-V的曲线特征。 (1)符合单一曲线的标准 (2)有无特殊情况出现。 (3)突出点分析。 A、定义：明显偏离所定平均曲线的反常关系点。 B、原因：测验和数据处理各环节中出现差错；因特殊水情造成关系点的突然变化。 C、判断标准：ISO标准用二号标准差。 D 、处理：系水力因素变化所致，则应作可靠资料看待；系测算错误造成，则应该设法予以改正；无法改正者，暂时不用，留待核实之后再定是否安全。 3、通过点群中心定一条单一线。 (1)误差最小原则。 (2)光滑无反曲 4、进行曲线检验 (1)符号检验； (2)适线检验； (3)偏离数值检验； (4)两条曲线需要合并时，还要进行t 检验。 4.4 正交函数法 正交多项式拟合法： 定义关于变量的一组相互正交函数： （ 正交多项式的阶数） 采用如下的递推方法求解正交函数： 当m=0时,P 0,j =1， 当m=1时，P 1，j =X j - 1， m m Z Z b Z Z b Z Z b b Q )()()(0202210-+-+-+= 0Z Z x -=Q y =m m x b x b x b b y += 22100 0=+k P P k i i ，m i ， 32,10=n j ， 32,1=（实测水位流量关系点总数） 当i =2,3m 时， 最后得到用于拟合多项式方程式的正交多项式： 4.5正交多项式拟合法程序设计 ' 正交多项式拟合水位流量关系 Private Sub Form_Load() Dim z(50) As Single, z0 As Single, Q(50) As Single, Qc(50) As Single, Z1(50) As Single, Q1(50) As Single Dim x(50) As Single, y(50) As Single, n1 As Single Dim a(5) As Single, alf(5) As Single, bet(5) As Single, bt(30) As Single Dim P(5, 50) As Single, Sum_p(5) As Single ' 输入实测数据 Open "C:Users郗鸿峰Desktop课件水信息技术-PPT1.txt" For Input As #1 n x p x n j j j j n j j = = 1 1 2,011 j i i j i i j j i p p x p ,21,1,)(- -= =-=-?= n j j i n j j i j i p p x 1 211 21， =-= n j j i n j j i i p p 1211 2，j m m j j j j p a p a p a p a y ,22,11,00+= =?= n j j i n j j j i i p y p a 1 21 ， ， Input #1, n For j = 1 To n Input #1, z(j) Next j For i = 1 To n Input #1, Q(i) Next i Input #1, n1 For i = 1 To n1 Input #1, Z1(i) Next i Close #1 ' 拟合多项式 z0 = 1.5 For j = 1 To n x(j) = Log(z(j) - z0) y(j) = Log(Q(j) Next j For j = 1 To n P(0, j) = 1 Next j Sum_p(0) = n s = 0 For j = 1 To n s = s + x(j) Next j alf(1) = s / Sum_p(0) For j = 1 To n P(1, j) = x(j) - alf(1) Next j m = 3 For i = 1 To m Sum_p(i) = 0 For j = 1 To n Sum_p(i) = Sum_p(i) + P(i, j) 2 Next j bet(i) = Sum_p(i) / Sum_p(i - 1) s = 0 For j = 1 To n s = s + x(j) * P(i, j) 2 Next j alf(i + 1) = s / Sum_p(i) For j = 1 To n P(i + 1, j) = (x(j) - alf(i + 1) * P(i, j) - bet(i) * P(i - 1, j) Next j Next i Print alf(1), alf(2), alf(3), bet(1), bet(2), bet(3) Sum_p(3) = 0 For j = 1 To n Sum_p(3) = Sum_p(3) + P(3, j) 2 Next j For i = 0 To m s = 0 For j = 1 To n s = s + P(i, j) * y(j) Next j a(i) = s / Sum_p(i) Next i Print a(0), a(1), a(2), a(3) For j = 1 To n Qc(j) = Exp(a(0) + a(1) * P(1, j) + a(2) * P(2, j) + a(3) * P(3, j) Next j '计算曲线拟合标准差Se s = 0 For i = 1 To n s = s + (Q(i) - Qc(i) * (Q(i) - Qc(i) Next i Se = Sqr(s / (n - 1) Open "C:Users郗鸿峰Desktop课件水信息技术-PPT2.txt" For Output As #2 Print #2, Print #2, Tab(8); "正交多项式拟合水位流量关系", Print #2, Tab(8); "_" Print #2, Print #2, Tab(3); "时刻" Tab(9); "水位" Tab(17); "流量" Tab(24); " 计算流量" Tab(35); " 相对误差%", Print #2, For j = 1 To n Print #2, Tab(3); Format(j, "0"); Tab(11); Format(z(j), "0.00"); Print #2, Tab(21); Format(Q(j), "0.0"); Tab(34); Format(Qc(j), "0.0"); Tab(45); Format(Qc(j) - Q(j) / Q(j) * 100, "0.00"), Next j Print #2, Print #2, Print #2, Tab(8); "标准差=" Format(Se, "0.00"); ' 符号检验 k = 0 For j = 1 To n If Qc(j) <= Q(j) Then k = k + 1 End If Next j u = (Abs(k - 0.5 * n) - 0.5) / (0.5 * Sqr(n) u_1 = 1.15