离散裂隙渗流方法与裂隙化渗透介质建模

流体渗流模拟的连续介质方法通常适用于多孔地质体，并不一定适用于裂隙岩体，由于裂隙分布及其特征与孔隙差异较大。若流体渗流主要受裂隙的控制，对于一定尺寸的裂隙岩体，多孔介质假设则较难刻划裂隙岩体的渗流特征。离散裂隙渗流方法不但可直接用于模拟裂隙岩体非均质性和各向异性等渗流特征，而且可用其确定所研究的裂隙岩体典型单元体及其水力传导（渗透）张量大小。主要讨论了以下问题：（1）饱和裂隙介质中一般的离散流体渗流模拟；（2）裂隙岩体中的REV（典型单元体）及其水力传导（渗透）张量的确定；（3）利用离散裂隙网络流体渗流模型研究裂隙方向几何参数对水力传导系数和REV的影响；（4）在二维和三维离散裂隙流体渗流模型中对区域大裂隙和局部小裂隙的处理方法。调查结果显示离散裂隙流体渗流数学模型可用来评价不同尺度上的裂隙岩体的水力特征，以及裂隙方向对裂隙化岩体的水力特征有着不可忽视的影响。同时，局部小裂隙、区域大裂隙应当区别对待，以便据其所起的作用及水力特征，建立裂隙化岩体相应的流体渗流模型。     

Groundwater-flow modeling in the Yucatan karstic aquifer, Mexico

The current conceptual model of the unconfined karstic aquifer in the Yucatan Peninsula, Mexico, is that a fresh-water lens floats above denser saline water that penetrates more than 40 km inland. The transmissivity of the aquifer is very high so the hydraulic gradient is very low, ranging from 7–10 mm/km through most of the northern part of the peninsula. The computer modeling program AQUIFER was used to investigate the regional groundwater flow in the aquifer. The karstified zone was modeled using the assumption that it acts hydraulically similar to a granular, porous medium. As part of the calibration, the following hypotheses were tested: (1) karstic features play an important role in the groundwater-flow system; (2) a ring or belt of sinkholes in the area is a manifestation of a zone of high transmissivity that facilitates the channeling of groundwater toward the Gulf of Mexico; and (3) the geologic features in the southern part of Yucatan influence the groundwater-flow system. The model shows that the Sierrita de Ticul fault, in the southwestern part of the study area, acts as a flow barrier and head values decline toward the northeast. The modeling also shows that the regional flow-system dynamics have not been altered despite the large number of pumping wells because the volume of water pumped is small compared with the volume of recharge, and the well-developed karst system of the region has a very high hydraulic conductivity. Electronic Publication     

Strike-slip fault bridge fluid pumping mechanism: insights from field-based palaeostress analysis and numerical modelling

We present a finite-element study of stress perturbation in evolving compressive and extensional strike-slip fault bridges. The results are compared with a fracture study of a compressive bridge at St Donats, UK. Horizontally interbedded calcareous mudstone and bioclastic calcilutite at St Donats have a distinct vertical permeability anisotropy. This sedimentary sequence behaves as a set of horizontal aquifers. The fluid flow in these aquifers is sensitive to mean stress gradients. Paleostress analysis of field fracture data, verified by finite-element modelling, indicates a rotation of σ₁ towards parallelism with boundary faults inside the growing compressive bridge. Boundary faults and bridge faults recorded numerous fluid flow events. The modelled mean stress pattern shows a regional maximum within the bridge and local maxima/minima pairs at boundary fault tips.Finite-element modelling of an extensional bridge indicates that σ₃ rotates towards parallelism with boundary faults. The mean stress pattern is similar to the pattern in compressive bridge but with maxima and minima locations interchanged. The stress patterns are reestablished by each stress build-up preceding the rupturation of the boundary faults throughout the development stages of strike-slip fault bridges. Mean stress gradients developed pre-failure control the fluid flow in fractures of the strike-slip fault system at and after the end of each stress build-up and the fluid flow in boundary faults post-failure. Fracture reactivation and new fracture generation within an evolving bridge is a process consisting of multiple successive events that retain the storage capacity of the bridge. Rupture and sealing of the main bounding-faults is a step-wise process that opens and closes fluid conduits between areas with different pressures.     

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…     

Computational modeling of cyclic mobility and post-liquefaction site response

Under seismic excitation, liquefied clean medium to dense cohesionless soils may regain a high level of shear resistance at large shear strain excursions. This pattern of response, known as a form of cyclic mobility, has been documented by a large body of laboratory sample tests and centrifuge experiments. A plasticity-based constitutive model is developed with emphasis on simulating the cyclic mobility response mechanism and associated pattern of shear strain accumulation. This constitutive model is incorporated into a two-phase (solid–fluid), fully coupled finite element code. Calibration of the constitutive model is described, based on a unique set of laboratory triaxial tests (monotonic and cyclic) and dynamic centrifuge experiments. In this experimental series, Nevada sand at a relative density of about 40% is employed. The calibration effort focused on reproducing the salient characteristics of dynamic site response as dictated by the cyclic mobility mechanism. Finally, using the calibrated model, a numerical simulation is conducted to highlight the effect of excitation frequency content on post-liquefaction ground deformations.     

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.     

Improved treatments for evaluating horizontal magnetic components through the 3-D FDM in E -polarization induction problems

To improve the accuracy of the numerical evaluation through the 3-D finite difference method, the surface boundary conditions are added to modify the old program. The author has tested the new program by making calculations for the model constructed by Wanamaker, et al (1984). The comparison between the numerical results obtained from this paper and those by Wannamaker, et al (1984) indicates that a pronounced improvement is realized in the evaluation of the horizontal magnetic components. Moreover, better calculations for the vertical magnetic components are also obtainable by using the new program.     

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

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

An oases dynamics model and its application in the Manas oasis

This study applies theories of systems and system dynamics to ecological economic systems in the oasis environment, with the intention of finding the basic characteristics of the oasis environment and its feedback structures. This study explores the inner mechanisms of economic development in the oasis environment and its relationship to policy variables. This model is applied to economic development in the Manas oasis of the Xinjiang Uygur Autonomous Region of China. System dynamics method starts from studying feedback structure within the system study and then searches for the best solution through identifying the right "policy point" and through model simulation in the computer. It embodies the cream of "systems cybernetics theory". This study is justified because the nature of the oasis system and the simulation results are scientifically sound and practically feasible.