A HORIZONTAL 2-D HYDRAULIC NUMERICAL MODEL AND IT’S APPLICATIONS TO FLOOD FORECAST

1 INTRODUCTION There were many flood disasters in China in recent years. When the water level in a river is very high, weak parts of its dike may be destroyed resulted in the submersion of the protected land and properties. It is of significance for decision-makers to exactly predict the processes of flood propagation during flood control. There are many modes of dike bursting, such as seepage destroying by overflow on top of dike caused by dike body sinking induced by piping and soil fl…     

Hybrid vibration experiments with a bridge foundation system model

In order to improve seismic design technology of bridges, it is necessary to evaluate the vibration characteristics of a bridge–soil system that consists of soil, foundation structure, pier and superstructure. However, there have been few experimental studies on seismic behavior of bridge–soil system. In this paper, we conducted the hybrid vibration experiment on seismic behavior of bridge–soil system, and examined the applicability of hybrid vibration experiment to study seismic response of bridge–soil system. Based on the experiment results, seismic response of bridge was quantitatively studied.     

Efficient methods for solving water flow in variably saturated soils under prescribed flux infiltration

The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.     

地形变场与现代地壳垂直运动及地震活动关系初探(Ⅰ)——地形变场的计算方法和物理意义

地震与地壳形变有着密切的关系 ,但台站地倾斜资料空间连续性不好 .为此在场分析的基础上 ,利用逐步订正法将地倾斜观测资料内插到 1°× 1°的经纬度网格点 ,引入曲面总曲率的概念计算格点上的Gauss曲率 ,并推广到格点化的有限区域上得到表征曲面形态和变形强度的地表形变场 .结果表明 :该方法可以合理地弥补资料不连续的缺陷 ,并在地表形变场上捕捉到一些关键的地表形变特征 ,如地凸区、地凹区、鞍形场等 .为地震分析和预报中如何拓展和利用大范围时空信息提供了一种有益的思路 .     

零价铁去除Cr（Ⅵ）的批实验研究

零价铁去除Cr(Ⅵ)的一系列批实验结果表明：水溶液中发生的氧化还原反应符合准一级反应；对于相同粒度的铁屑(比表面积2．89m^2／g)，固液比由05g／100mL增加到2．0g／100mL，Cr(Ⅵ)质量浓度达到排放标准所需反应时间从26．7min降为3．6min；在固液比相同、原水中Cr(Ⅵ)质量浓度不超过30mg／L时，原水质量浓度增大，达到饮用水标准的反应时间增加；其他条件相同，升高温度可提高零价铁与Cr(Ⅵ)的反应速率。     

Land subsidence caused by groundwater exploitation in Suzhou City,China

Suzhou City, located at the lower reaches of the Yangtze River in southeastern Jiangsu Province, is one of the few cities in China which suffer from severe ground settlement. A research project was carried out to investigate this problem. Geological and hydrogeological studies show that there is a multi-layered aquifer system with three distinct, soft mud layers of marine and lagoonal origins. An examination of historical records of groundwater extraction, water levels, and ground settlement shows that the ground subsidence is associated with the continuously increasing groundwater extraction in the deep, confined aquifer. It is believed that the consolidation of the soft mud layers, especially the third layer which is thick and close to the main pumped aquifer, contributes to the ground settlement. A three-dimensional finite difference numerical model representing the multi-layered aquifer system was developed to study the ground settlement in response to groundwater extraction. By calibrating the model with both the measured groundwater level and ground settlement, the aquifer parameters were estimated. The model outputs fit reasonably well with the observed results, which indicates that the numerical model can reproduce the dynamic processes of both groundwater flow and soil consolidation. The hydraulic conductivity of the third mud layer near the center of the ground settlement has been reduced by over 30% in the last 14 years. The gradual deterioration in the hydraulic conductivity of the mud may have significant adverse effect on the sustainable groundwater resource of the deep confined aquifer, since the recharge from the shallow aquifers through the mud layer is the only source of water to the deep aquifer. An analysis of the spatial distributions of groundwater drawdown and ground settlement shows that the area with maximum drawdown is not necessarily the area with maximum ground settlement due to the occurrence of the soft mud layer. A simple reallocation in pumping rates on the basis of the spatial distribution of the thick mud layer could significantly reduce the ground settlement. Electronic Publication     

Pull-down basin in the central part of Japan due to subduction-induced mantle flow

The driving force for the basin subsiding against isostatic balance in and around Lake Biwa in the Kinki district, Japan is discussed. The lake region is characterized by strong negative Bouguer anomalies, especially by a steep horizontal gradient zone of gravity anomaly running along the western margin of the lake. The large negative anomaly (>50 mgal) cannot be explained by low-density sediments beneath it. A down-warping structure extending to the Moho depth should be taken into account. This conjecture has been strongly supported by a short-period receiver function imaging, which shows a clear offset of about 8 km for the Moho discontinuity under the steep gravity gradient zone.A question arises as to what is the driving force to create such a large down-warping structure. We consider that the subduction of the shallow-dipping slab under the region (Philippine Sea Slab) may cause crustal deformation by dragging the viscous mantle downward. In order to verify this model, we simulated the induced mantle flow due to the subduction of the Philippine Sea Slab and the pressure distribution on the crust–mantle boundary. This numerical experiment showed that the induced flow makes a strong negative pressure zone under the lake region if the slab has a vertical offset along the direction of subduction. This offset of the slab is consistent with plate models deduced from hypocentral distributions and Sp phases of the deep-focus earthquakes.     

南海各月月平均海流场的数值模拟

本文用一个区域三维、自由表面、斜压海洋模式对南海全年各月的三维平均流场进行了数值模拟，模拟得到的环境特征与己有的实测资料十分相近，这表明该模拟结果是可信的，该数值模拟也是成功的。     

湖南省2001年夏秋季对流云降水潜力数值模拟结果

应用对流云数值模式模拟了2001年7～9月湖南省三个探空站0000和1200(世界时)462个对流云降水算例.模拟结果表明(1)其中53天的198个算例属于有利人工增雨的天气形势,133个算例的可播度大于零,54个算例的增雨率大于零;(2)7～9月天气形势有利于人工增雨日的算例,对流云平均含水量为575万吨,平均降水效率11.4%;(3)8月降水潜力大于7月,9月降水潜力最小.这表明即使大旱的2001年湖南省夏秋季对流云仍然有一定的增雨潜力.     

一种在异构系统中实现负载平衡的方法

提出了在异构系统实现负载平衡的区域分解算法和实现负载平衡的计算方法，利用它的负反馈性质解决了异构系统处理机计算速度测量误差造成的负载测量不准问题，并对处理机速度变化，速度测量误差、处理机数量、网格点计算量的分布等因素的影响进行了计算，结果表明本方法具有很强的平衡负载能力和较强的适应性；根据计算结果提出了解决模式网格点计算量不易测量问题的解决方案，并用扩散方程和模拟物理过程进行试验，试验表明这种方法是可行的，平衡负载的效果十分显著。