扩散波方程数值解及其在洪水演算中的应用
|
摘要
扩散波洪水演算法是一种既具有足够精度又便于求解的方法。本文对扩散波洪水演算进行了较为系统的研究，着重探讨了其数值差分解。
基于运动波数值扩散在一定条件下可以模拟扩散波的物理扩散，探讨了马斯京根—康吉洪水演算法和一种具有预见期的差分解法；分析了方法的相容性、稳定性和收敛性，讨论了采用不同方式处理边界的相容性、稳定性和收敛性条件；介绍了运动波差分方程模拟扩散波方程近似解的精度准则，给出了二阶、三阶和四阶精度解的最优条件。利用直接差分格式的多样性，给出了扩散波方程的两种显式直接差分解、一种隐式直接差分解和三种混合直接差分解，讨论了以上6种差分方程及其边界处的相容性、稳定性和收敛性条件。
对汉江襄阳至皇庄河段和西江梧州至高要河段，应用扩散波解析解法、马斯京根—康吉洪水演算法、具有预见期的差分解法和6种直接差分法对其进行洪水演算；采用遗传算法根据一定的准则对训练样本的各场次洪水的模型参数进行率定，进而利用信息扩散技术由训练样本的模型参数获得校核样本的模型参数；并给出通过不同方法获得的模型参数的校核样本的预报结果。最后，比较了不同方法的扩散波洪水演算结果，结果表明，不同方法对两条河段的演算精度均较高，通过优化和信息扩散方法得到的样本模型参数好于直接采用物理方法获得模型参数的计算结果。
It is known that the diffusion-wave flood routing method is a good one with accuracy and simplicity. In this thesis, this method is systematically studied, especially the numerical difference solution to the equations.
Based on principle that numerical dispersion of the linear kinematic-wave can simulate physical dispersion of the diffusion wave on certain condition, Muskingum-Cunge method and a flood routing method with forecasting period are discussed; the compatibility, stability and convergence of the two methods are analyzed, and the conditions to handle the compatibility, stability and convergence at the boundary in different methods; accurate principles of the approximation of that simulation are introduced, and optimal conditions necessary to obtain second-, third-and fourth-order accurate solutions are also given. By making use of the diversity of the finite difference scheme, two explicit direct difference, one implicit direct difference and three mixed direct difference schemes are given, and the above six kinds of difference equations together with their compatibility, stability and convergence conditions at the boundary are discussed.
In the case of Xiangyang-Huangzhuang channel (part of Hanjiang river of Changjiang) and Wuzhou-Gaoyao channel (part of Xijiang river of Zhujiang), methods of the diffusion wave analysis and Muskingum-Cunge, flood routing method with forecasting period, and six finite difference solutions are applied to the calculation of the diffusion wave; according to certain rules, real coding based accelerating genetic algorithm is used to optimize the model parameters of studied samples, from which the model parameters of the verified samples are obtained by the way of information diffusive technique; in addition, the forecasting outcome the model parameters of the verified samples through different methods are given.
After comparing the results of different methods, the outcome shows that all methods have satisfactory precision, and model parameters of samples gained by optimization and information diffusive technique are better than those by physical means.
基于运动波数值扩散在一定条件下可以模拟扩散波的物理扩散，探讨了马斯京根—康吉洪水演算法和一种具有预见期的差分解法；分析了方法的相容性、稳定性和收敛性，讨论了采用不同方式处理边界的相容性、稳定性和收敛性条件；介绍了运动波差分方程模拟扩散波方程近似解的精度准则，给出了二阶、三阶和四阶精度解的最优条件。利用直接差分格式的多样性，给出了扩散波方程的两种显式直接差分解、一种隐式直接差分解和三种混合直接差分解，讨论了以上6种差分方程及其边界处的相容性、稳定性和收敛性条件。
对汉江襄阳至皇庄河段和西江梧州至高要河段，应用扩散波解析解法、马斯京根—康吉洪水演算法、具有预见期的差分解法和6种直接差分法对其进行洪水演算；采用遗传算法根据一定的准则对训练样本的各场次洪水的模型参数进行率定，进而利用信息扩散技术由训练样本的模型参数获得校核样本的模型参数；并给出通过不同方法获得的模型参数的校核样本的预报结果。最后，比较了不同方法的扩散波洪水演算结果，结果表明，不同方法对两条河段的演算精度均较高，通过优化和信息扩散方法得到的样本模型参数好于直接采用物理方法获得模型参数的计算结果。
It is known that the diffusion-wave flood routing method is a good one with accuracy and simplicity. In this thesis, this method is systematically studied, especially the numerical difference solution to the equations.
Based on principle that numerical dispersion of the linear kinematic-wave can simulate physical dispersion of the diffusion wave on certain condition, Muskingum-Cunge method and a flood routing method with forecasting period are discussed; the compatibility, stability and convergence of the two methods are analyzed, and the conditions to handle the compatibility, stability and convergence at the boundary in different methods; accurate principles of the approximation of that simulation are introduced, and optimal conditions necessary to obtain second-, third-and fourth-order accurate solutions are also given. By making use of the diversity of the finite difference scheme, two explicit direct difference, one implicit direct difference and three mixed direct difference schemes are given, and the above six kinds of difference equations together with their compatibility, stability and convergence conditions at the boundary are discussed.
In the case of Xiangyang-Huangzhuang channel (part of Hanjiang river of Changjiang) and Wuzhou-Gaoyao channel (part of Xijiang river of Zhujiang), methods of the diffusion wave analysis and Muskingum-Cunge, flood routing method with forecasting period, and six finite difference solutions are applied to the calculation of the diffusion wave; according to certain rules, real coding based accelerating genetic algorithm is used to optimize the model parameters of studied samples, from which the model parameters of the verified samples are obtained by the way of information diffusive technique; in addition, the forecasting outcome the model parameters of the verified samples through different methods are given.
After comparing the results of different methods, the outcome shows that all methods have satisfactory precision, and model parameters of samples gained by optimization and information diffusive technique are better than those by physical means.
引文
[1] 刘国纬，沈国昌．中国年最大致洪暴雨落区研究[J]．水科学进展，2006，17(2)：151-159．
[2] 魏立和．北运河洪水预报模型研究[D]．天津：天津大学，2004：1-15．
[3] 周轶．水文学与水力学方法相结合的洪水预报模型[D]．南京：河海大学，2004：1-20．
[4] 韩程起，张儒生，刘默．河道洪水演算方法应用研究[J]．黑龙江水专学报，1999，26(3)：36-38．
[5] 芮孝芳．产汇流理论[M]．北京：水利电力出版社，1995：35-95．
[6] 芮孝芳．水文学原理[M]．北京：中国水利水电出版社，2004：212-245．
[7] 左广巍．河道洪水演算方法的研究与应用[M]．华中科技大学．2004：1-25．
[8] 芮孝芳，李琼芳，王伶俐．线性扩散波洪水演算模型研究[J]．水文，1999，(6)：3-7．
[9] 黄国如．扩散波方程解的研究及应用[D]．南京：河海大学，1999：1-60．
[10] 芮孝芳，姜广斌．洪水演算理论与计算方法的若干进展与评论[J]．水科学进展，1998，9(4)：299-395．
[11] 黄国如，芮孝芳．扩散波洪水演算研究进展[J]．水利水电技术，2000，31(8)：31-34．
[12] 丁伯良，黄国如，芮孝芳．线性扩散波洪水演算方法的应用研究[J]．水电能源科学，2001，19 (3)：43-46．
[13] Hayami S. On the propagation of the flood waves[J]. Disaster prevention research. 1951,(1):55-80.
[14] Gonwa, W. S., Kawas, M. L, A modified diffusion equation for flood propagation in rapezoidal channels[J], Journal of Hydrology, 1986, (83): 119-136.
[15] Keefer, T. N., Mcquivey, R S., Multiple linearization flow routing modes[J], Journal of hydraulics division, ASCE, 100 (7): 103-104.
[16] Suthcrland, A. J., Barnett, A. G., Diffusion solution to flows with upstream control[J], Journal of hydraulics division, ASCE 1972, 98 (11): 1969-1982.
[17] Dooge J. A. Practical aspects of computation river hydraulics[M]. Pitman Publishing Limited. London.1979:1-52
[18] Hayami S. On the propagation of the flood waves. Disaster prevention research. Inst. Bull. 1, 1951, (1):55-80.
[19] Tawatchai Tingsanchali. Analytical diffusion model for flood routing. Journal of hydraulic Engineering. 1985, 111 (1): 435-454.
[20] Gonwa, W. S., Kawas, M. L, A modified diffusion equation for flood propagation in rapezoidal channels[J], Journal of Hydrology, 1986, (83): 119-136.
[21] Bernard Cappelaere. Accurate diffusive wave routing. J. of hydr. engineering, 1997, 23 (3): 174-181.
[22] 芮孝芳．运动波数值扩散与洪水演算方法[J]．水利学报，1987，(2)：37-43．
[23] 芮孝芳，王伶俐．具有预见期的洪水演算方法研究[J]．水科学进展，2000，11(3)：291-295．
[24] Cunge J. A.. On the subject of a flood propagation computation method (Muskingum method). J. Hydrol. Res., IAHR, 1969, 7: 205-230.
[25] Ponce V. M. Simplified Muskingum routing equation. J. Hydrol. Div, ASCE, 1979, 105 (HY1): 85-90.
[26] Koussis A. D., Unified theory for flood and pollution routing. J. Hydrol. Eng., ASCE, 1983, 109 (12): 1552-1664.
[27] Bajracharya K., Barry D. A. Accuracy criteria for linearised diffusion wave flood routing. Journal of Hydrology, 1997, (195): 200-217.
[28] 芮孝芳．Nuskingum法及其分段连续演算的若干理论探讨[J]．水科学进展，2002，12(6)：682-688．
[29] 芮孝芳，王伶俐．具有预见期的洪水演算方法研究[J]．水科学进展，2000，11(3)：291-295．
[30] 黄国如，芮孝芳．基于运动波数值扩散的洪水演算方法[J]．河海大学学报，2001，29(2)：10-113．
[31] 王船海，李光炽．实用河网水流计算，河海大学水文系讲义，1995：1-50．
[32] 李庆扬，王能超，易大义．数值分析[M]．北京：清华大学出版社，2001：1-102．
[33] 王文洽．求解扩散方程的一类交替分组显式方法[J]．山东大学学报，2002，37(3)：194-199．
[34] 黄素珍．求解对流扩散方程的一类交替分组显式方法[J]．盐城工学院学报，2004，17(1)：12-16．．
[35] 郑兴华．二维对流扩散方程的分步交替分组显式格式[J]．华侨大学学报，2002，23(2)： 122-128.
[36] Holland J. H. Adaptation in Natural and Artificial Systems[M]. (2nd ed) Cambridge MA: MIT Press, 1992:1-65.
[37] Goldberg D. E. Genetic Algorithms in Search, Optimization, and Machine Learning[M]. New York: Addison-wesley Publishing Company, INC, 1989:1-180.
[38] 周明，孙树栋．遗传算法及其应用[M]．．北京：国防工业出版社，1999：1-60．
[39] 金菊良，丁晶．遗传算法及其在水科学中的应用[M]．成都：四川大学出版社，2000：1-120．
[40] 潘正君，康立山，陈毓屏．演化计算[M]．北京：清华大学出版社，1998：1-60．
[41] 云庆夏，黄光球，王战权、遗传算法和遗传规划-一种寻优搜索技术[M]．北京：冶金工业出版社，1997：15-85．
[42] 周玉良．水安全问题中的智能非参数方法[D]．合肥：合肥工业大学，2006：12-83．
[43] 刘值勤．人工生命理论及其应用[M]．北京：冶金工业出版社，1997：1-100．
[44] 黄崇福，王家鼎．模糊信息优化处理技术及其应用[M]．北京：北京航空航天大学出版社，1995：1-74．
[45] 王新洲．基于信息扩散原理的估计理论、方法及其抗差性[J]．武汉测绘科技大学学报，1999，24(3)：240-244．
[46] 黄崇福，张俊香，刘静．模糊信息优化处理技术应用简介[J]．信息与控制，2004，33(1)：61-66
[47] Huang C F, Shi Y. Towards Efficient Fuzzy Information Progressing - Using the Principle of Information Diffusion [M]. Hei delberg: Physica Verlag (Springer), 2002:1-63.
[48] Huang C F. Deriving samples from incompleted data [M]. Anchorage alaska: ALTEC, 1998:645-650.
[49] 陈志芬，黄崇福，张俊香．基于扩散函数的内集-外集模型[J]．模糊系统与数学，2006，20(1)：42-48．
[50] 张有富．模糊集方法在水安全问题中的应用研究[D]．合肥：合肥工业大学，2006：1-75．
[2] 魏立和．北运河洪水预报模型研究[D]．天津：天津大学，2004：1-15．
[3] 周轶．水文学与水力学方法相结合的洪水预报模型[D]．南京：河海大学，2004：1-20．
[4] 韩程起，张儒生，刘默．河道洪水演算方法应用研究[J]．黑龙江水专学报，1999，26(3)：36-38．
[5] 芮孝芳．产汇流理论[M]．北京：水利电力出版社，1995：35-95．
[6] 芮孝芳．水文学原理[M]．北京：中国水利水电出版社，2004：212-245．
[7] 左广巍．河道洪水演算方法的研究与应用[M]．华中科技大学．2004：1-25．
[8] 芮孝芳，李琼芳，王伶俐．线性扩散波洪水演算模型研究[J]．水文，1999，(6)：3-7．
[9] 黄国如．扩散波方程解的研究及应用[D]．南京：河海大学，1999：1-60．
[10] 芮孝芳，姜广斌．洪水演算理论与计算方法的若干进展与评论[J]．水科学进展，1998，9(4)：299-395．
[11] 黄国如，芮孝芳．扩散波洪水演算研究进展[J]．水利水电技术，2000，31(8)：31-34．
[12] 丁伯良，黄国如，芮孝芳．线性扩散波洪水演算方法的应用研究[J]．水电能源科学，2001，19 (3)：43-46．
[13] Hayami S. On the propagation of the flood waves[J]. Disaster prevention research. 1951,(1):55-80.
[14] Gonwa, W. S., Kawas, M. L, A modified diffusion equation for flood propagation in rapezoidal channels[J], Journal of Hydrology, 1986, (83): 119-136.
[15] Keefer, T. N., Mcquivey, R S., Multiple linearization flow routing modes[J], Journal of hydraulics division, ASCE, 100 (7): 103-104.
[16] Suthcrland, A. J., Barnett, A. G., Diffusion solution to flows with upstream control[J], Journal of hydraulics division, ASCE 1972, 98 (11): 1969-1982.
[17] Dooge J. A. Practical aspects of computation river hydraulics[M]. Pitman Publishing Limited. London.1979:1-52
[18] Hayami S. On the propagation of the flood waves. Disaster prevention research. Inst. Bull. 1, 1951, (1):55-80.
[19] Tawatchai Tingsanchali. Analytical diffusion model for flood routing. Journal of hydraulic Engineering. 1985, 111 (1): 435-454.
[20] Gonwa, W. S., Kawas, M. L, A modified diffusion equation for flood propagation in rapezoidal channels[J], Journal of Hydrology, 1986, (83): 119-136.
[21] Bernard Cappelaere. Accurate diffusive wave routing. J. of hydr. engineering, 1997, 23 (3): 174-181.
[22] 芮孝芳．运动波数值扩散与洪水演算方法[J]．水利学报，1987，(2)：37-43．
[23] 芮孝芳，王伶俐．具有预见期的洪水演算方法研究[J]．水科学进展，2000，11(3)：291-295．
[24] Cunge J. A.. On the subject of a flood propagation computation method (Muskingum method). J. Hydrol. Res., IAHR, 1969, 7: 205-230.
[25] Ponce V. M. Simplified Muskingum routing equation. J. Hydrol. Div, ASCE, 1979, 105 (HY1): 85-90.
[26] Koussis A. D., Unified theory for flood and pollution routing. J. Hydrol. Eng., ASCE, 1983, 109 (12): 1552-1664.
[27] Bajracharya K., Barry D. A. Accuracy criteria for linearised diffusion wave flood routing. Journal of Hydrology, 1997, (195): 200-217.
[28] 芮孝芳．Nuskingum法及其分段连续演算的若干理论探讨[J]．水科学进展，2002，12(6)：682-688．
[29] 芮孝芳，王伶俐．具有预见期的洪水演算方法研究[J]．水科学进展，2000，11(3)：291-295．
[30] 黄国如，芮孝芳．基于运动波数值扩散的洪水演算方法[J]．河海大学学报，2001，29(2)：10-113．
[31] 王船海，李光炽．实用河网水流计算，河海大学水文系讲义，1995：1-50．
[32] 李庆扬，王能超，易大义．数值分析[M]．北京：清华大学出版社，2001：1-102．
[33] 王文洽．求解扩散方程的一类交替分组显式方法[J]．山东大学学报，2002，37(3)：194-199．
[34] 黄素珍．求解对流扩散方程的一类交替分组显式方法[J]．盐城工学院学报，2004，17(1)：12-16．．
[35] 郑兴华．二维对流扩散方程的分步交替分组显式格式[J]．华侨大学学报，2002，23(2)： 122-128.
[36] Holland J. H. Adaptation in Natural and Artificial Systems[M]. (2nd ed) Cambridge MA: MIT Press, 1992:1-65.
[37] Goldberg D. E. Genetic Algorithms in Search, Optimization, and Machine Learning[M]. New York: Addison-wesley Publishing Company, INC, 1989:1-180.
[38] 周明，孙树栋．遗传算法及其应用[M]．．北京：国防工业出版社，1999：1-60．
[39] 金菊良，丁晶．遗传算法及其在水科学中的应用[M]．成都：四川大学出版社，2000：1-120．
[40] 潘正君，康立山，陈毓屏．演化计算[M]．北京：清华大学出版社，1998：1-60．
[41] 云庆夏，黄光球，王战权、遗传算法和遗传规划-一种寻优搜索技术[M]．北京：冶金工业出版社，1997：15-85．
[42] 周玉良．水安全问题中的智能非参数方法[D]．合肥：合肥工业大学，2006：12-83．
[43] 刘值勤．人工生命理论及其应用[M]．北京：冶金工业出版社，1997：1-100．
[44] 黄崇福，王家鼎．模糊信息优化处理技术及其应用[M]．北京：北京航空航天大学出版社，1995：1-74．
[45] 王新洲．基于信息扩散原理的估计理论、方法及其抗差性[J]．武汉测绘科技大学学报，1999，24(3)：240-244．
[46] 黄崇福，张俊香，刘静．模糊信息优化处理技术应用简介[J]．信息与控制，2004，33(1)：61-66
[47] Huang C F, Shi Y. Towards Efficient Fuzzy Information Progressing - Using the Principle of Information Diffusion [M]. Hei delberg: Physica Verlag (Springer), 2002:1-63.
[48] Huang C F. Deriving samples from incompleted data [M]. Anchorage alaska: ALTEC, 1998:645-650.
[49] 陈志芬，黄崇福，张俊香．基于扩散函数的内集-外集模型[J]．模糊系统与数学，2006，20(1)：42-48．
[50] 张有富．模糊集方法在水安全问题中的应用研究[D]．合肥：合肥工业大学，2006：1-75．