的特点。为此,郑邦民、齐鄂荣17]提出了用水尺站的实测值与计算值的误差函数和控制方程离散后非线性代数方程组的余量函数共同组成目标函数,并用无约束求极值的方法求解此目标函数。这一方法有助于天然河道糙率的深入研究。Р70年代以来,计算机技术的发展,促进了明渠~维非恒定流计算,出现了~批通用程序,其中规模较大功能较完善的软件包有DWOPER【81、PROGStg]等。DWOPER 是美国气象局的河流预报系统的一部分。这些软件系统能够迅速为特定河流建立数学模型,可联机修改和率定,但其适用性还不够普遍,差分格式少,边界条件的组合少。Р谭维炎【1哪针对上述软件包的不足,开发了明渠一维非恒定流计算程序包 MYBC,它采用完全形式的圣维南方程组,提供了上下游9种边界条件组合,可选择显式差分、隐式差分和特征格式共三类九种典型算法。但未考虑断面流速分布不均匀系数及其变化、局部水头损失和漫滩的影响。庚维德等【【I】采用Preissmarm隐式解,扩充了MYBC,开发出明渠一维非恒定流计算程序包MYBC/YI,计入了断面突然变化引起的局部水头损失和漫滩影响。Р近年来一个更为先进和功能较为完整的水力分析软件~Hee-Ras--是由美国陆军工程兵团开发的‘12】,系统主要用于一维恒定流及非恒定流的水面线计算,同时包括泥沙等沉积物的输移和一些水力设计的计算。Hec.Ras模型将丁坝考虑成一个堰的形式,该形式不仅考虑了丁坝的挡水作用,而且考虑了丁坝的壅水作用。该软件一个突出的优点是地形数据的输入方便快捷,计算结果输出生动直观。Р尽管上述一维恒定流及非恒定流的水面线计算方法在实际工程中得到广泛应用,但对河床地形复杂或象丁坝这样引起断面形状突变的情况,计算误差仍是较大的。另外,其中关键性的糙率和局部损失系数一般都通过实测资料来估计,这就大大削弱了该计算方法对未建成工程的预测功能。因此,直到现在,在航道整治中的