环状河道二维浅水流动的全隐数值求解方法

基于通用河网二维水流模式,采用全隐式矩阵追赶与回代技术,提出了环状河道单元二维数值模拟的通用模型。该模型仿照一维河网的建模方法,通过建立节点水位方程组,能方便地实现全流域二维河网或者一二维河网的隐式耦合求解。通过对长江下游南京八卦洲段分汊河道的数值模拟,表明该方法精度较高,切实有效,能满足实际工程的需要。     

基于未确知信息的环状河网水质数学模型研究

以圣维南方程组、对流扩散方程以及未确知信息理论为基础,建立了环状河网一维水动力与未确知水质数学模型。采用Preissman四点隐格式法对圣维南方程组进行离散和采用显式差分对对流扩散方程进行离散,对河网水动力与水质分别进行了合理的编码,并考虑了泵闸等控制工程、降雨径流的影响,采用河网的三级解法编制河网非恒定流的计算程序求解各断面的水位Z和流量Q,实现了水动力和水质两模块之间的数据连接,且对微分形式的水质模型运用未确知数学理论进行了含有两个未确知参数的水质计算,从而不仅能获得污染物浓度区间值,还能得到相应的可信度。实例研究表明,建立的河网水动力模型和未确知信息的水质模型是可靠的,可以用于河网水质治理工程的数值模拟研究和水环境评估。     

一维-二维耦合的河湖系统整体水动力模型

为对自然界中的复杂水流系统进行整体水动力数值模拟,建立了适用于河湖系统的一维-二维耦合模型.将河网计算中用于处理河段间耦合的汊点水位预测校正法应用到一维-二维耦合边界的处理中,即耦合边界水位预测校正法.该方法具有一维、二维模型计算完全独立,可有效利用现有模型的优点.与传统一维-二维耦合处理方法相比,耦合连接条件的满足程度较高且可通过对计算容差的调整进行有效控制.利用理想及实际例子进行了计算验证,结果表明该方法具有较好的实用性.     

基于Godunov格式的溃坝水流数学模型

为了更好地把握溃坝洪水风险,减小因溃坝洪水而造成的人员生命和财产损失,建立了基于Godunov格式的一维、二维溃坝水流耦合数学模型。一维溃坝水流模型采用HLL格式的有限体积法求解,二维溃坝水流模型采用基于非结构网格的Roe格式离散求解,在一维、二维模型的链接处采用重叠计算区域的方法实现一维模型和二维模型之间的水力要素信息交换。经弯道溃坝算例和断面突变溃坝算例验证,该耦合模型具有良好的可靠性和适用性,验证后的耦合模型为大尺度的溃坝水流数值模拟打下了基础。     

三峡水库调节作用对长江近河口段水文水动力特性影响

为了评估三峡水库蓄水对长江感潮河段水文水动力特性的影响,建立了具有预测功能的水动力模型。考虑径流与潮流之间的相互作用,模型在只有上边界流量过程的条件下可以给出下边界的潮位过程。在此基础上应用一维河网非恒定流数学模型,对三峡水库蓄水前后长江近河口段大通—徐六泾的水流运动进行了10年的数值模拟。统计分析结果表明,水库调蓄作用对感潮河段各站的潮位累积频率及水文水动力因素的年内过程变化影响甚大。     

基于广义垂线坐标系的三维非结构数学模型及其在珠江口的应用

基于广义垂线坐标变换,构建了非结构网格的三维水动力数学模型。模型采用半隐式有限体积法对控制方程进行数值离散,其中半隐式法用于水位梯度和垂向紊动扩散以及垂向对流项的离散,显式控制体积分法用于水平对流项等的离散。广义垂线坐标系使得模型能够灵活地对垂向网格进行布置,平面非结构网格使得模型能够适应河口海岸复杂的岸线,并可对局部进行网格加密。模型通过具有解析解的风生流和异重流对模型进行检验,应用模型模拟了珠江口的三维潮流过程,计算结果与实测水文数据的比较表明,模型能够较好地模拟珠江河口的盐水楔和三维分层水动力过程。     

闸控河网水文-水动力-水质耦合数学模型——I.理论

针对流域级闸控大型河网水环境日常管理的应用需求与实际特点,将资料适应性强的水文学方法与数据要求相对严苛的水动力-水质数值模型相融合,构建了一维与二维嵌套、分块组合的闸控大型河网水文-水动力-水质耦合数学模型DHQM（Hydrology, Hydrodynamics, and Water Quality Model for Impounded Rivers）。模型由河道径流模拟、闸坝调度过程模拟、河道水质模拟、区间入流及入河污染负荷估算和水质预警实时校正等5个模块组成。模型可服务于水环境实时预警和调度,也可为闸坝水文环境效应的量化提供基础工具。     

主槽一维和滩地二维侧向耦合洪水演进模型

为实现对复式河道大尺度漫滩洪水的快速模拟,结合洪水在滩槽行洪特点和河道一维、二维水动力学模型研究成果,在主槽与滩地分别划分为矩形和三角形网格,以左右滩唇为分界实现两网格系统嵌套;建立主槽一维水动力学数学模型,内插大断面间网格处的水力因子,在滩地采用简化的二维水动力学模型模拟滩槽间、滩地网格间水量交换,实现主槽一维与滩地二维模型侧向耦合。利用黄河下游“96·8”大洪水实测资料,对花园口至夹河滩河段洪水演进过程验证,结果表明该方法模拟精度较高,能够利用该模型模拟复式河道大漫滩洪水过程。     

煤层气羽状水平井数值模拟模型的交替隐式算法

煤层气开采过程中通常采用数值模拟方法进行生产指标预测,其中的数值模型普遍采用三维全隐式差分法求解,模型涉及的参数变量多,求解难度大,计算复杂。根据煤层气的储气特征、吸附特征和解吸机理建立了反映煤层气解吸、扩散及渗流过程的气、水两相流动煤层气羽状水平井三维数学模型,模型考虑了煤层非均质各向异性的影响,运用块中心差分格式,对气、水相偏微分方程进行空间差分和时间差分,采用交替隐式(ADI)法求解,实现了煤层气羽状水平井开发数值模拟计算。交替方向格式可以把二维问题转化成一维问题,对x,y两个方向的迭代矩阵均为三对角矩阵,结构相同,易于编程计算。所建立的差分格式具有计算量小,稳定性好等优点,数值试验的结果表明效果良好。     

闸控河网水文-水动力-水质耦合数学模型——Ⅱ.应用

针对中国多闸坝重污染河流的典型代表——淮河中游河网,为提高其流域水环境预警及应急响应能力,依据水文-水动力-水质耦合数学模型DHQM(Hydrology, Hydrodynamics, and Water Quality Model )的理论框架,构建了由淮河干流、沙颍河和涡河3个一维模块和鲁台子、蚌埠闸2个二维模块组成的淮河中游河网水文-水动力-水质耦合数学模型.采用经过验证的模型针对一个虚拟的水污染事故情景,分析了不同调度方案与淮河干流来水组合下污染团下泄对下游水质的影响;并以保障淮河干流蚌埠市饮用水源地的取水水质安全和减轻事故对下游水质的影响程度为目标,优选出分期小流量慢速下泄污染水体为最佳的水量水质联合应急调度方案.     

浅谈地下水流模拟中河流的处理方法

本文分析了河水与地下水不同作用形式下的补排关系,给出了河流在地下水流模拟中的两种处理方法.根据河流与地下水的相互作用规律将其概化为四种边界条件,并给出了在不同边界条件中河水与地下水交换量的计算公式,初步建立了河水与地下水的解析耦合模型.     

Theis不稳定潜水井流模型的改进——具入渗补给

地下水的补排主要包括垂向的地面入渗补给、蒸发排泄（蒸发可视为入渗的负值）及侧向的地表水补给、排泄。水文地质学最基本的问题之一——地下水可持续开釆量的评价准则,涉及补给的增量与排泄的减量,因此地下水开采的预测模型必须包含上述两类的补给、排泄因素,否则不能满足要求。然而,经典的Theis不稳定井流模型,即使在傍河抽水,也只有侧边界的补给、排泄作用,而不涉及上边界的地面入渗补给。这样一来,这个解析模型基本上不能够用于预测,而只能在旱季用于井流试验求取含水系统的参数。为此,文章的目标是发展具地面入渗补给的Theis不稳定潜水井流模型。对于潜水流问题,不能再用承压水流的以水头为应变量的方程来建立,应采用第二类线性化方法的势函数来建立潜水流问题。对于既有降雨入渗补给,又有抽水井作用的复杂的水文地质问题所概化数学模型的求解,采取的方法是把它分解成若干个简单的子模型问题求解,然后将其合成为原来复杂数学模型的解。基于质量守恒原理,假定渗流服从Darcy定律并满足Dupuit徦定建立了水流基本微分方程。然后对于两平行河流及一河流平行一隔水边界形成的两类条形区域,具地面均匀稳定入渗补给的井流问题,获得通用水位方程和几类常见的特定条件水位方程及其流量方程。此外,提出并采用“边界对边界的反映法”用以求解一河流平行一隔水边界条形区域的同一问题,减少了许多推导过程。最后,作为上述理论成果的初步应用,也是一个重要的应用,即在河水水质不能满足要求的河流附近,设有一口抽水井,计算该抽水井在不汲取河水的前提下的临界流量方程,获得具重要意义的结构简洁的关系式。该方程也可以用于滨海区的抽水井,在不发生海水入侵前提下的临界抽水流量计算。给出了上述条件不稳定井流过程某时刻的地下水流网图,其流网与文献中常见的傍河井流的流网相比,具显明的特征。     

Dynamic routing modeling for flash flood forecast in river system

A real-time flood-forecasting method coupled with the one-dimensional unsteady flow model was developed for the Danshuei River system in northern Taiwan. Based on the flow at current time, the flow at new time is calculated to provide the water stage forecasting during typhoons. Data, from two typhoons in 2000: Bilis and Nari, were used to validate and evaluate the model capability. First, the developed model was applied to validate and evaluate with and without discharge corrections at the Hsin-Hai Bridge in Tahan Stream, Chung-Cheng Bridge in Hsintien Stream, and Sir-Ho Bridge in the Keelung River. The results indicate that the calculated water stage profiles approach the observed data. Moreover, the water stage forecasting hydrograph with discharge correction is close to the observed water stage hydrograph and yields a better prediction than that without discharge correction. The model was then used to quantify the difference in prediction between different methods of real-time water stage correction. The model results reveal that water stages using the 1–6 h forecast with real-time stage correction exhibits the best lead times. The accuracy for 1–3 h lead time is higher than that for 4–6 h lead time, suggesting that the flash flood forecast in the river system is reasonably accurate for 1–3 h lead time only. The method developed is effective for flash flood forecasting and can be adopted for flood forecasting in complicated river systems.     

钱塘江河口段盐水入侵的时空变化及预测模型

杭州城市供水85%取自钱塘江河口段,取水水质在枯水大潮期都不同程度地受到盐水入侵的威胁,分析钱塘江河口盐水入侵时空变化及研制二维数值预测模型对保障城市供水安全十分必要。根据钱塘江河口段实测水文氯度资料,分析了强潮作用下盐水入侵的时空变化特征;据此构建考虑斜压作用的二维水流、盐度输移的耦合数学模型,计算格式采用守恒性较好的有限体积法;在模型验证的基础上,数值分析了径流和潮汐对钱塘江河口段盐水入侵的影响,结果表明河口段的盐水入侵明显地受径流和潮汐的影响,据此可通过增大上游新安江水库的下泄流量抑制盐水入侵上溯以减小取水口氯度及超标时间,确保用水安全。     

An automatic quadrilateral mesh generation algorithm applied to 2-D overland flow simulations

This paper presents an automatic quadrilateral mesh generation algorithm to provide 2-D meshes for simulating overland flow and transport problems in watershed systems. The conformal mapping method is studied deeply and specific treatment strategies are established according to the width of rivers and the area of other water storage zones. A conformity treatment method is proposed to deal with the incompatibility problems during discretizing the complex geometry with a plurality of sub-domains. This method ensures grid conformity through removing certain overlapped nodes and creating new appropriate nodes on common boundaries of the multi-domain mesh. Aiming at the undesired deviation of mesh boundaries from geometric contours, a relative position-percent interpolation method is proposed and a supplementary match method related to dead ends is presented. And, the double effects of projection and smooth for boundary nodes are achieved. An objective function method is proposed to locally optimize the degenerated quadrilaterals located on concave or large-curvature curves. In order to accommodate the need for overland simulations, the 1-D/2-D correspondence on river reaches, junctions with and without storage, ponds, lakes, dead ends and control structures is established. Finally, practical applications are provided to demonstrate the accuracy and reliability of the quadrilateral meshing algorithms proposed in this paper.     

A modeling study of the Satilla River estuary,Georgia. II: suspended sediment

A three-dimensional (3-D) suspended sediment model was coupled with a 3-D hydrodynamic numerical model and used to examine the spatial and temporal distribution of suspended sediments in the Satilla River estuary of Georgia. The hydrodynamic model was a modified ECOM-si model with inclusion of the flooding-drying cycle over intertidal salt marshes. The suspended sediment model consisted of a simple passive tracer equation with inclusion of sinking, resuspension, and sedimentation processes. The coupled model was driven by tidal forcing at the open boundary over the inner shelf of the South Atlantic Bight and real-time river discharge at the upstream end of the estuary, with a uniform initial distribution of total suspended sediment (TSS). The initial conditions for salinity were specified using observations taken along the estuary. The coupled model provided a reasonable simulation of both the spatial and temporal distributions of observed TSS concentration. Model-predicted TSS concentrations varied over a tidal cycle; they were highest at maximum flood and ebb tidal phases and lowest at slack tides. Model-guided process studies suggest that the spatial distribution of TSS concentration in the Satilla River estuary is controlled by a complex nonlinear physical process associated with the convergence and divergence of residual flow, a non-uniform along-estuary distribution of bottom stress, and the inertial effects of a curved shoreline.     

A 2.5-D coupled FE–BE model for the dynamic interaction between saturated soil and longitudinally invariant structures

A two-and-a-half-dimensional (2.5-D) coupled finite element–boundary element (FE–BE) model is presented to simulate the three-dimensional dynamic interaction between saturated soils and structures with longitudinally invariant geometries. A regularized 2.5-D boundary integral equation for saturated porous media is derived that avoids the evaluation of singular traction integrals. The 2.5-D coupled FE–BE model is established by using the continuity conditions on the soil–structure interface. The developed model is verified through comparison with an existing semi-analytical method. Two case studies of a tunnel embedded in a poroelastic half-space and the efficiency of a vibration isolating screen are presented.     

Development of Tide–Surge Interaction Model for the Coastal Region of Bangladesh

In this paper, a two-dimensional, vertically integrated hydrodynamic model is developed taking into account entrained air bubbles during storm surges as well as incorporating inverted barometer, and river and land dynamics. The model is specifically designed for the coastal region of Bangladesh. A nested scheme method with a fine mesh scheme (FMS), capable of incorporating the complex coastline and all major offshore islands accurately, nested into a coarse mesh scheme (CMS) covering up to 15° N latitude in the Bay of Bengal is used. To incorporate the small and big offshore islands in the Meghna estuarine region with its complex coastline accurately, a very fine mesh scheme (VFMS) is again nested into the FMS. Along the northeast corner of the VFMS, the Meghna river discharge is taken into account. The coastal and island boundaries are approximated through proper stair steps. The model equations are solved by a semi-implicit finite difference technique using a staggered C-grid. A stable appropriate tidal condition over the model domain is generated by applying tidal forcing with the four major tidal constituents M ₂, S ₂, K ₁, and O ₁ along the southern open boundary of the CMS. This tidal regime is introduced as the initial state of the sea for nonlinear interaction of tide and surge. The model is applied to simulate water levels due to the interaction of tide and surge associated with the cyclones April 1991 and Aila at different coastal and island locations along the coast of Bangladesh. The results are found to be quite satisfactory with root mean square error of ～0.50 m as calculated for both the storm events. Tests of sensitivities on water levels are carried out for air bubbles, offshore islands, river discharge, inverse barometer, and grid resolution. The presence of air bubbles increases simulated water levels a little bit in our model, and the contribution of air bubbles in increasing water level is found around 2 %. Further, water levels are found to be influenced by offshore islands, river discharge, inverse barometer as well as grid resolution.     

水质大涡模拟数学模型研究

建立了基于弱可压缩流基础上的水质大涡模拟数学模型及污染物质扩散系数与平均运动的关系。数值计算采用有限体积法和预测校正法,固壁边界条件采用"部分滑移条件"。并应用该模型对重庆嘉陵江、长江两江交汇段的流场,污染物质浓度分布进行了计算,计算结果与实测值比较吻合。研究表明,横向扩散与流动状态有关,当支流流量相对较小时,横向扩散较弱,形成岸边污染带,当支流流量相对较大时,横向扩散较强,污染带向江心移动。     

干旱区间歇性河流河道两侧地下水运动模型及其求解方法

本文根据间歇性河流河道充水和停水过程周期性重复出现的水文特征,建立了河道两侧区域地下水运动的一维非稳定流模型,提出了把流量边界条件与水位界条件相互转换来求解模型的一种方法.最后应用上述模型分析西北某内陆河下游间歇性输水条件下河道两侧地区地下水位恢复状况,为输水生态效应的定量评价提供了理论基础.