Delft3D在鳌江口外平阳咀海域流场模拟中的应用

水文Ｖｏｌ．２７Ｎｏ．６Ｄｅｃ．，２００７

ＪＯＵＲＮＡＬＯＦＣＨＩＮＡＨＹＤＲＯＬＯＧＹ

Ｄｅｌｆｔ３Ｄ在鳌江口外平阳咀海域流场模拟中的应用

左书华

（交通部天津水运工程科学研究所工程泥沙交通行业重点实验室，天津３００４５６）

摘要：Ｄｅｌｆｔ３Ｄ是以ＦＬＯＷ水动力模块为主体，建立在Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的基础上，由一系列模块组成的模型软件，可

以对水流、波浪、泥沙、水质、地貌等建立二、三维模型。本文介绍了Ｄｅｌｆｔ３Ｄ数值模拟软件的程序模块及数值理论基本方程，对鳌江口外平阳咀海域的流场进行了模拟研究。研究表明：计算结果与实测结果吻合较好，能够较好地反映该海区真实的流场情况，可以为工程研究提供基本的流场背景。关键词：数值模拟；Ｄｅｌｆｔ３Ｄ；河口海岸工程；应用中图分类号：Ｐ３３８＋．９

文献标识码：Ａ

文章编号：１０００－０８５２（２００７）０６－００５５－０４

传播、非线性波浪之间的折射以及由于形条件下，波浪的产生、

波浪破碎和底摩擦引起的耗散。此模块里主要包含着两个模型，ＨＩＳＷＡ波浪模型和ＳＷＡＮ波浪模型，ＳＷＡＮ模型可以研究宽阔海域中的波浪，而ＨＩＳＷＡ模型主要研究近岸波浪，其中

随着计算机技术的发展，数值模拟技术也得到空前的发展，数值模拟是综合流体力学、计算机数学以及各种生产应用技术而发展成长起来的一门科学。它具有耗资省、速度快、修改灵便等优点，其应用的领域已日益广泛，在研究河口及河口工程的可行性方面也越来越普遍，并被工程界所认可和采纳。目前国际上应用较多的河口海洋数值模拟软件有美国普林斯顿大学的

ＨＩＳＷＡ模型是Ｄｅｌｆｔ３Ｄ提供的标准波浪模块，不过ＳＷＡＮ模型

也是可以用的。

（３）质点跟踪模块ＰＡＲＴ。它能通过及时的质点跟踪评估一个动态的浓度分布过程，用浓度等值线的形式来详细描述盐度、油、温度或其他极易衰减物质的瞬间或连续释放过程。主要应用在近岸污水排放、海上溢油扩散、河口赤潮的发生、以及电厂的热力排放等方面。

它是一个三维的水质模拟框架，它能解（４）水质模块ＷＡＱ。

决预先确定好的计算网格里的水平对流扩散问题，它的水流信息来自Ｆｌｏｗ模块。它可以用来模拟营养盐、有机物质、溶解氧、重金属等问题。ＢＯＤ、ＣＯＤ、

ＰＯＭ、ＥＣＯＭ模式，丹麦水力学研究所的ＭＩＫＥ２１软件，荷兰ＤＥＬＦＴ水力学研究所的Ｄｅｌｆｔ３Ｄ软件，还有美国Ｂｒｉｇｈａｍ大学环

境模型研究实验室开发的ＳＭＳ软件等等。

Ｄｅｌｆｔ３Ｄ软件是由荷兰ＤＥＬＦＴ水力学研究所开发的，集水

流、泥沙、环境于一体的程序软件包，既可以进行二维计算，也可以进行三维计算，可以进行潮流泥沙输移计算、台风风暴潮计算、温排水计算、水质计算、溢油扩散、质点跟踪模拟等。

１Ｄｅｌｆｔ３Ｄ包含的程序模块

Ｄｅｌｆｔ３Ｄ软件是由一系列模块组成，具有高度的可视化性，

每个模块都有自己单独的可视化界面，图１为Ｄｅｌｆｔ３Ｄ３．２６主菜单界面。Ｄｅｌｆｔ３Ｄ主要包括：水动力模块ＦＬＯＷ、波浪模块水质模块ＷＡＱ、泥沙输移模块ＳＥＤ、地形模块ＭＯＲ、质ＷＡＶＥ、

点跟踪模块ＰＡＲＴ、生态模块ＥＣＯ、化学模块ＣＨＥＭ、工具模块

Ｕｔｉｌｉｔｉｅｓ，其中在３．２６版本里ＳＥＤ、ＭＯＲ模块包含在水动力模块

生态模块ＥＣＯ包含在ＷＡＱ里。ＦＬＯＷ里，化学模块ＣＨＥＭ、

（１）水动力模块ＦＬＯＷ。Ｆｌｏｗ模块是一个三维的水动力和输初运模块，也可以进行二维计算，在此模块里，给定开边界条件、始条件以及有关的参数，如水位、流速、流量和糙率等，关于泥沙和地形的计算也包含在此模块里。与此模块相关联的模块有

ＷＡＱ、ＷＡＶＥ、ＰＡＲＴ。

地（２）波浪模块ＷＡＶＥ。此模块主要是用来研究一定水流、

收稿日期：２００７－０３－２７

作者简介：左书华（１９７９－），男，河北邯郸人，硕士，主要从事河口海岸泥沙及海岸工程研究。

图１

Ｄｅｌｆｔ３Ｄ３．２６主菜单界面

５６

水文第２７卷

物质传输方程

传输方程在平面上采用正交曲线坐标系，垂向上采用&坐

另外，软件还包含着一些前期、后期处理的程序，如网络生产、处理程序ＲＧＦＧＲＩＤ和ＱＵＩＣＫＩＮ；后期图像显示与处理程序

２．３

ＧＰＰ和ＱＵＩＣＫＰＬＯＴ；潮汐调和分析程序ＴＩＤＥ。标系，方程如下：

２Ｄｅｌｆｔ３Ｄ软件数值模型理论基础［１］

Ｄｅｌｆｔ３Ｄ是以Ｆｌｏｗ水动力模块为主体的多模块模型，水动

!（ｄ＋\"）ｃ＋１

!ｔ!Ｇξξ!Ｇ$$!%!Ｇ&$$（ｄ＋\"）ｕｃ

$!\"$（!Ｇξξ（ｄ＋\"）ｖｃ

＋

Ｈ

$$ξξ

力模块建立在Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的基础上，求解基于有限差分法—ＡＤＩ（ＡｌｔｅｒｎａｔｉｎｇＤｉｒｅｃｔｉｏｎＩｍｐｌｉｃｉｔ）法［２］。

!ξＡＤＩ法实质上是把时间步长Δｔ分成２个半步长，前半步在

方向用隐格式，η方向用显格式，后半步在η方向用隐格式，在ξ

方向用显格式，这样可以把较大的系数矩阵化为２个三角形ξ

系数矩阵，可用求解。“追赶法”

，η），在垂直Ｄｅｌｆｔ３Ｄ基本方程是建立在正交曲线坐标上（ξ方向上采用&坐标，表示如下：

ｄ＋\"!Ｇξξ!Ｇ$$

!Ｄ!Ｇ&%!Ｇ!ξ＋!%ｃ＝’!&!$!ｃＤ!Ｇ!ｃ$＋!%$’!$!$!Ｇ!$Ｈ

ξξ$$

!ｃ

＋１!ｃＤＶ－’ｄ（ｄ＋\"）ｃ＋Ｓ

!&ｄ＋\"!&%$（６）

式中：!Ｇξ!Ｇ$$为曲线坐标系至直角坐标系的转换系数；ξ、、、Ｕ、Ｖ分别为ξ$方向的平均流速；ｕ、ｖ和%分别为ξ$和&方向的流速；ｆ为科氏力系数；Ｐξ和Ｐ$是ξ、$方向压强梯度；Ｆξ和

（１）

、、$方向的紊动动量通量；Ｍξ和Ｍ$分别是ξ$方向上动Ｆ$是ξ

量的源或汇；ｃ为保守物质浓度；Ｓ为源汇项，’ｄ为第一项衰减。

&＝ｚ－\"＝ｚ－\"

ｄ＋\"Ｈ式中：ｚ为在垂直方向的坐标，在模型参照平面上取为０，往下到底部为ｄ；ｄ为相对模型参照平面的水深；\"为相对于模型参照平面的水位；Ｈ为全水深，Ｈ＝ｄ＋\"。

３Ｄｅｌｆｔ３Ｄ在鳌江口外平阳咀海域流场模拟中的应用

在河口海岸工程的研究中，数学模型已经成为重要的研究

２．１连续方程

沿水深积分平均的二维潮流连续性方程：

手段，逐步替代了过去必须是物理模型试验才能解决的问题，而且应用范围越来越广，几乎涉及河口海岸工程的所有分支。

$!%（ｄ＋\"）Ｕ!Ｇξ!\"＋１ξ

＋

!!ｔ!Ｇξξξ!Ｇ$$$!%（ｄ＋\"）Ｖ!Ｇξ１ξ

＝Ｑ

!$ＧＧ!ξξ!$$

纬度；Ｒ为地球半径。

（２）

Ｄｅｌｆｔ３Ｄ数值模拟软件在河口海岸工程中的应用也越来越广，储

天文潮鏖［３］利用Ｄｅｌｆｔ３Ｄ以长江口地区为例建立了台风风暴潮、耦合数值预报模型，对台风风暴潮、天文潮两潮耦合预报模式进行探研和分析；刘成等［４］利用Ｄｅｌｆｔ３Ｄ对长江口水动力条件，上海市现有及拟建排放口污水排放的稀释扩散场进行了模拟；栗苏文等

［５］

式中：Ｑ为每个单元上的源项；!Ｇξξ＝ＲＣＯＳ)，!Ｇ$$＝Ｒ；*为

在污染源调查和污染负荷估算的基础上，利用

［６］

２．２动量方程

，$方向的动量方程为：ξ

Ｄｅｌｆｔ３Ｄ对大鹏湾现状和未来条件下的水质进行了数值模拟，并

采用模型试算法估算出大鹏湾的水环境容量；黄坚等江口流场的影响进行研究。

应用

!ｕ＋ｕ!ｕ＋ｖ!ｕ＋%!ｕ

ξｄ＋\"!&!ｔ!Ｇξ!!$!Ｇ$$ξ

ｕｖｖ２!!Ｇξ!!Ｇ$$－ｆｖξ

＋－

!$!ξ!Ｇξ!Ｇξξ!Ｇ$$ξ!Ｇ$$

１＝－Ｐξ＋Ｆξ＋１２!（ｖｖ!ｕ）＋Ｍξ

（ｄ＋\"）!&!&+０!Ｇξξ!ｖ＋ｕ!ｖ＋ｖ!ｖ＋%!ｖξｄ＋\"!&!ｔ!Ｇξ!!$!Ｇ$$ξ

ｕｖｕ２!!Ｇ$$－!!Ｇξξ

＋－ｆｕ

!ξ!$ＧＧＧ!Ｇξ!!!ξξξ$$$$

１＝－Ｐ$＋Ｆ$＋１２!（ｖｖ!ｖ）＋Ｍ$

（ｄ＋\"）!&!&+０!Ｇ$$

垂向速度%在&坐标系中由下式计算得出：

Ｄｅｌｆｔ３Ｄ通过数值模拟深水航道工程及横沙东滩促淤工程对长３．１３．１．１

（３）

二维潮流数学模型的建立边界条件

图３。最大网格间距５００ｍ，最模型计算区域及网格如图２、

小网格间距８０ｍ，南、东、北三侧为开边界，有潮位过程控制，潮位由中国海域潮汐模型［７］计算给出；由于鳌江口流量较小，给一多年平均流量值６５ｍ３／ｓ。

３．１．２初始条件和参数选取

初始条件为上节方程初始时刻的潮位和流速，计算时取流

（４）

速为０；初始潮位取各边界的平均值。根据海域水深不同，海底糙率取０．０１８～０．０２２。水平涡粘系数和扩散系数均取１００ｍ２／ｓ。

３．１．３模型建立过程

本文应用Ｄｅｌｆｔ３Ｄ对鳌江口外平阳咀水域进行了二维潮流

!\"（ｄ＋\"）ｕ!Ｇ$$１!\"＋

!ｔ!Ｇξ!ξξ!Ｇ$$

#（５）

数值模拟，过程主要包括：①根据陆地边界在ＲＧＦＧＲＩＤ里生产网格文件；②根据实际需要在ＱＵＩＣＫＩＮ里对网格文件进行处理，如水深文件的生产、糙率的设定、观测点的设置等；③然后

$!%!\"（ｄ＋\"）ｖ!Ｇξ１ξ＋＋＝Ｈ（ｑｉｎ－ｑｏｕｔ）

!$!&ＧＧ!ξξ!$$

第６期左书华：Ｄｅｌｆｔ３Ｄ在鳌江口外平阳咀海域流场模拟中的应用

５７

琵琶门

计算实测

在水动力模块Ｆｌｏｗ里输入网格文件、水深文件等，设置时间、初始条件、给定开边界和各种参数，生成ｍｄｆ文件，进行运行调试。

５．００

０．００

１２：００１６：００２０：０００：００－５．００５．００

０．００

１２：００１６：００２０：０００：００－５．００

图４

４：００平阳嘴

８：００１２：００１６：００２０：００

计算实测

４：００８：００１２：００１６：００２０：００

潮位站潮位验证（大潮）

计算实测

３．００２．００１．０００．００１４：００３．００２．００１．０００．００１４：００

１８：００

２２：００

１８：００

２２：００

１＃

６：００２＃

１０：００１４：００

计算实测

２：００３＃

６：００１０：００１４：００

计算实测

３．００２．００１．０００．００１４：００

１８：００

２２：００图５

２：００６：００１０：００１４：００

测站潮流流速验证

４００３００２００１０００１４：００

３６０３００２４０１８０１２０６００

１＃

计算实测

１８：００２２：００２：００２＃

６：００１０：００１４：００

计算实测

１４：００１８：００２２：００２：００３＃

６：００１０：００１４：００

计算实测

３．２模型验证和结果分析

模型验证采用２００５年３月２６日至２７日该海域实测大潮

潮位、水流资料进行验证。

各潮位站潮位以及流速过程验证结果见图４￣图６。潮位基面为８５基面；各潮流测点的流速过程图中表示为实测值，“■”实线代表计算结果。验证结果表明：从潮位验证的结果看，计算潮位的相位与实测潮位的相位吻合较好，相位误差一般不超过

４００

３００２００１０００１４：００

１８：００２２：００图６

２：００６：００１０：００１４：００

测站潮流流向验证

０．２５ｈ，峰谷值误差一般不大于０．１５ｍ；流速流向率定结果表明，

流速相位一般也不超过０．２５ｈ，流速峰谷值的绝对误差不大于涨落潮最大流速值绝对值之和的２０％。因此来说，计算结果能够较好地反映该海区真实的流场情况，可以为工程研究提供基本的流场背景。

４结语

本文对Ｄｅｌｆｔ３Ｄ软件的组成模块、理论基础（方程）作了介

绍，并对其具体应用作了实例分析。Ｄｅｌｆｔ３Ｄ是以ＦＬＯＷ水动力模块为主体，建立在Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的基础上，由一系列模块组成的模型，以采用正交曲线网格、计算稳定、精度高等特点

５８

水文第２７卷

在海岸带、河口、风暴潮的研究中得到了较好地应用，为我国的海岸、河口的研究提供了一种比较理想的模型。

不过，模型对网格的正交性要求相当的高，Ｄｅｌｆｔ３Ｄ是基于正交曲线网格的，对于岸线曲折复杂的海域来说，建立满足要求的曲线正交网格难度较大，也降低了工作效率，许多研究者多采用矩形网格进行计算，无疑是降低了岸线的拟和度。因此有待于集成一套更好的网格生成方法。参考文献：

［１］ＷＬ．ＤｅｌｆｔＨｙｄｒａｕｌｉｃｓ，Ｄｅｌｆｔ３Ｄ－ＦｌｏｗＵｓｅｒｍａｎｕａｌ［Ｒ］．２００５，１２．［２］Ｌｅｅｎｄｅｒｔｓｅ，Ｊ．Ｊ．ａｎｄＥ．Ｃ．Ｇｒｉｔｔｏｎ．Ａｗａｔｅｒｑｕａｌｉｔｙｓｉｍｕｌａｔｉｏｎｍｏｄｅｌ

ｆｏｒｗｅｌｌｍｉｘｅｄｅｓｔｕａｒｉｅｓａｎｄｃｏａｓｔａｌｓｅａｓ：Ｖｏｌ．ＩＩ，Ｃｏｍｐｕｔａｔｉｏｎ

Ｐｒｏｃｅｄｕｒｅｓ，ＲｅｐｏｒｔＲ－７０８－ＮＹＣ，ＴｈｅＲａｎｄＣｏｒｐｏｒａｔｉｏｎ（ＳａｎｔａＭｏｎｉｃａ）［Ｍ］，１９７１．

［３］储鏖．Ｄｅｌｆｔ３Ｄ在天文潮与风暴潮耦合数值模拟中的应用［Ｊ］．海洋预报，２００４，２１（３）：２９￣３６．

［４］刘成，李行伟，韦鹤平，等．长江口水动力及污水稀释扩散模拟［Ｊ］．海

洋与湖沼，２００３，３４（５）：４７４￣４８２．

［５］栗苏文，李红艳，夏建新．基于Ｄｅｌｆｔ３Ｄ模型的大鹏湾水环境容量分

析［Ｊ］．环境科学研究，２００５，１８（５）：９１￣９５．

［６］黄坚，何青，桑永尧，等．长江口南北槽及横沙东滩工程流场数值模拟

［Ｊ］．水运工程，２００３，３５１（４）：５￣９．

［７］李孟国．中国海域潮汐预报软件Ｃｈｉｎａｔｉｄｅ的应用［Ｊ］．水道港口，

２００７，２８（１）：６５￣６８．

ＩｎｔｒｏｄｕｃｔｉｏｎｏｆＮｕｍｅｒｉｃａｌＳｉｍｕｌａｔｉｏｎＳｏｆｔｗａｒｅＤｌｅｆｔ３Ｄ＆ＩｔｓＡｐｐｌｉｃａｔｉｏｎｏｎ

ＯｆｆｓｈｏｒｅＡｒｅａｏｆＡｏｊｉａｎｇＥｓｔｕａｒｉｎｅ

ＺＵＯＳｈｕＨｕａ

（ＫｅｙｌａｂｏｒａｔｏｒｙｏｆＥｎｇｉｎｅｅｒｉｎｇＳｅｄｉｍｅｎｔｏｆＭｉｎｉｓｔｒｙｏｆＣｏｍｍｕｎｉｃａｔｉｏｎｓ，ＴｉａｎｊｉｎＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅｏｆＷａｔｅｒＴｒａｎｓｐｏｒｔＥｎｇｉｎｅｅｒｉｎｇ，Ｔｉａｎｊｉｎ

３００４５６，Ｃｈｉｎａ）

Ａｂｓｔｒａｃｔ：ＴｈｅＤｅｌｆｔ３Ｄｓｕｉｔｅｈａｓｄｅｖｅｌｏｐｅｄａｕｎｉｑｕｅ，ｆｕｌｌｙｉｎｔｅｇｒａｔｅｄｃｏｍｐｕｔｅｒｓｏｆｔｗａｒｅｓｕｉｔｅｆｏｒａｍｕｌｔｉ－ｄｉｓｃｉｐｌｉｎａｒｙａｐｐｒｏａｃｈａｎｄ３Ｄｃｏｍｐｕｔａｔｉｏｎｓｆｏｒｃｏａｓｔａｌ，ｒｉｖｅｒａｎｄｅｓｔｕａｒｉｎｅａｒｅａｓ．Ｉｔｃａｎｃａｒｒｙｏｕｔｓｉｍｕｌａｔｉｏｎｓｏｆｆｌｏｗｓ，ｓｅｄｉｍｅｎｔｔｒａｎｓｐｏｒｔｓ，ｍｏｒｐｈｏｌｏｇｉｃａｌｄｅｖｅｌｏｐｍｅｎｔｓ，ｗａｖｅｓ，ｗａｔｅｒｑｕａｌｉｔｙ，ａｎｄｅｃｏｌｏｇｙ．Ｉｔｈａｓｂｅｅｎｄｅｓｉｇｎｅｄｆｏｒｅｘｐｅｒｔｓａｎｄｎｏｎ－ｅｘｐｅｒｔｓａｌｉｋｅ．ＴｈｅＤｅｌｆｔ３ＤｓｏｆｔｗａｒｅｉｓｃｏｍｐｏｓｅｄｏｆｓｅｖｅｒａｌｍｏｄｕｌｅｓｂａｓｉｎｇｏｎＦＬＯＷｈｙｄｒｏｄｙｎａｍｉｃｍｏｄｕｌｅ，ｇｒｏｕｐｅｄａｒｏｕｎｄａｍｕｔｕａｌｉｎｔｅｒｆａｃｅ，ｗｈｉｌｅｂｅｉｎｇｃａｐａｂｌｅｔｏｉｎｔｅｒａｃｔｗｉｔｈｏｎｅａｎｏｔｈｅｒ．ＴｈｅａｕｔｈｏｒｄｅｓｃｒｉｂｅｓｔｈｅｔｈｅｏｒｅｔｉｃａｌｐｈｙｓｉｃｓｍｏｄｅｌｔｏＤｅｌｆｔ３Ｄ－ＦＬＯＷ，ｗｈｉｃｈｉｓｂａｓｅｄｏｎｔｈｅＮａｖｉｅｒ－ＳｔｏｋｅｓｅｑｕａｔｉｏｎｓａｎｄｉｎｔｒｏｄｕｃｅｓａｐｐｌｉｃａｔｉｏｎｏｆＤｅｌｆｔ３ＤｏｎｏｆｆｓｈｏｒｅａｒｅａｏｆＡｏｊｉａｎｇｅｓｔｕａｒｉｎｅｉｎｔｈｉｓｐａｐｅｒ．Ｋｅｙｗｏｒｄｓ：ＮｕｍｅｒｉｃａｌＳｉｍｕｌａｔｉｏｎ；Ｄｌｅｆｔ３Ｄ；ＥｓｔｕａｒｉｎｅａｎｄＣｏａｓｔａｌＥｎｇｉｎｅｅｒｉｎｇ；Ａｐｐｌｉｃａｔｉｏｎ

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!（上接第５页）

［４２］ＪａｙａｗａｒｄｅｎａＡＷ，ＧｕｒｕｎｇＡＢ．Ｎｏｉｓｅｒｅｄｕｃｔｉｏｎａｎｄｐｒｅｄｉｃｔｉｏｎｏｆｈｙｄｒｏｍｅ－

ｔｅｏｒｏｌｏｇｉｃａｌｔｉｍｅｓｅｒｉｅｓ：ｄｙｎａｍｉｃａｌｓｙｓｔｅｍｓａｐｐｒｏａｃｈｖｓ．ｓｔｏｃｈａｓｔｉｃａｐ－ｐｒｏａｃｈ［Ｊ］．ＪｏｕｒｎａｌｏｆＨｙｄｒｏｌｏｇｙ，２０００，２２８（３－４）：２４２－２６４．

［４３］ＥｌｓｈｏｒｂａｇｙＡ，ＳｉｍｏｎｏｖｉｃＳＰ，ＰａｎｕＵＳ．Ｎｏｉｓｅｒｅｄｕｃｔｉｏｎｉｎｈｙｄｒｏｌｏｇｉｃｔｉｍｅ

ｓｅｒｉｅｓ：ｆａｃｔｓａｎｄｄｏｕｂｔｓ［Ｊ］．ＪｏｕｒｎａｌｏｆＨｙｄｒｏｌｏｇｙ，２００２，２５６（３－４）：１４７－１６５．［４４］ＳｍｉｔｈＬＡ．Ｉｎｔｒｉｎｓｉｃｌｉｍｉｔｓｏｎｄｉｍｅｎｓｉｏｎｃａｌｃｕｌａｔｉｏｎｓ［Ｊ］．ＰｈｙｓＬｅｔｔＡ，

１９８８，１３３（６）：２８３－２８８．

［４５］ＮｅｒｅｎｂｅｒｇＭＡＨ，ＥｓｓｅｘＣ．Ｃｏｒｒｅｌａｔｉｏｎｄｉｍｅｎｓｉｏｎａｎｄｓｙｓｔｅｍａｔｉｃｇｅｏｍｅｔ－

ｒｉｃｅｆｆｅｃｔｓ［Ｊ］．ＰｈｙｓｉｃａｌＲｅｖｉｅｗＡ，１９９０，４２（１２）：７０６５－７０７４．

［４６］ＰａｓｔｅｒｎａｃｋＧＢ．Ｄｏｅｓｔｈｅｒｉｖｅｒｒｕｎｗｉｌｄ？Ａｓｓｅｓｓｉｎｇｃｈａｏｓｉｎｈｙｄｒｏｌｏｇｉｃａｌ

ｓｙｓｔｅｍｓ［Ｊ］．ＡｄｖａｎｃｅｓｉｎＷａｔｅｒＲｅｓｏｕｒｃｅｓ，１９９９，２３（３）：２５３－２６０．

［４７］陈云浩，史培军，李晓兵．不同热力背景对城市降雨（暴雨）的影响（ｌ）一

降雨时序的混饨分析［Ｊ］．自然灾害学报，２００１，１０（３）：２０－２５．

［４８］杨思全，陈亚宁．河冰湖突发洪水的分形和混沌特征研究［Ｊ］．干旱区

地理，１９９９，２２（２）：７７－８２．

［４９］丁晶，王文圣，赵永龙．反演水文动力模型的探讨［Ｊ］．水力发电学报，

２００２，（７８）：４７－１１．

ＣｕｒｒｅｎｔＳｔａｔｕｓａｎｄＰｒｏｓｐｅｃｔｏｆＣｈａｏｓＴｈｅｏｒｙＲｅｓｅａｒｃｈｉｎＨｙｄｒｏｌｏｇｙ

ＬＩＨｏｎｇ－ｘｉａ，ＸＵＳｈｉ－ｇｕｏ，ＸＵＸｉａｎｇ－ｚｈｏｕ，ＷＡＮＧＦｕ－ｑｉａｎｇ

（ＳｃｈｏｏｌｏｆＣｉｖｉｌａｎｄＨｙｄｒａｕｌｉｃＥｎｇｉｎｅｅｒｉｎｇ，ＤａｌｉａｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ，Ｄａｌｉａｎ１１６０２４，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｔｈｅｏｂｊｅｃｔｉｖｅｏｆｔｈｉｓｐａｐｅｒｉｓｔｏｐｒｅｓｅｎｔｔｈｅａｃｈｉｅｖｅｍｅｎｔｓａｎｄｉｍｐｏｒｔａｎｔｉｓｓｕｅｓｉｎｔｈｅａｐｐｌｉｃａｔｉｏｎｏｆｃｈａｏｓｔｈｅｏｒｙｉｎｈｙｄｒｏｌｏｇｙ，ａｎｄｔｏｄｉｓｓｅｍｉｎａｔｅｔｈｅｐｏｔｅｎｔｉａｌｆｏｒｆｕｔｕｒｅｉｍｐｒｏｖｅｍｅｎｔｓ．Ｔｈｉｓｐａｐｅｒｒｅｖｉｅｗｅｄｔｈｅａｐｐｌｉｃａｔｉｏｎｏｆｃｈａｏｔｉｃｉｄｅｎｔｉｆｉｃａｔｉｏｎ，ｐｒｅｄｉｃｔｉｏｎ，ｐａｒａｍｅｔｅｒｓｆｏｒｐｈａｓｅｓｐａｃｅｒｅｃｏｎｓｔｒｕｃｔｉｏｎ，ｎｏｉｓｅｐｒｏｃｅｓｓｉｎｇａｎｄｄａｔａｓｉｚｅｉｎｈｙｄｒｏｌｏｇｙ．

Ｔｈｉｓｐａｐｅｒａｌｓｏｐｏｉｎｔｅｄｏｕｔｔｈａｔｃｈａｏｔｉｃａｐｐｒｏａｃｈｓｈｏｕｌｄｂｅｃｏｍｐａｒｅｄａｎｄｃｏｍｂｉｎｅｄ

ｗｉｔｈｏｔｈｅｒｔｈｅｏｒｉｅｓａｎｄｔｈｅｒｅｓｅａｒｃｈｏｆｍｕｌｔｉｖａｒｉａｔｅｐｒｅｄｉｃｔｉｏｎａｎｄｒｅｃｏｎｓｔｒｕｃｔｉｏｎｏｆｓｙｓｔｅｍｄｙｎａｍｉｃａｌｅｑｕａｔｉｏｎｓｍａｙｒｅｖｅａｌｔｈｅｅｖｏｌｕｔｉｏｎａｎｄｍｅｃｈａｎｉｓｍｏｆｃｈａｏｔｉｃｈｙｄｒｏｌｏｇｙ．

Ｋｅｙｗｏｒｄｓ：ｈｙｄｒｏｌｏｇｉｃｔｉｍｅｓｅｒｉｅｓ；ｃｈａｏｓｉｄｅｎｔｉｆｉｃａｔｉｏｎａｎｄｐｒｅｄｉｃｔｉｏｎ；ｉｍｐｏｒｔａｎｔｉｓｓｕｅｓ；ｒｅｓｅａｒｃｈｐｒｏｃｅｓｓ