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     

Thermal and rheological structures of the Xisha Trough, South China Sea

The Xisha Trough, located in the northwest of the South China Sea (SCS) mainly rifted 30 Ma ago, has been a failed rift since the cessation of the seafloor spreading of the NW subbasin. Based on the velocity–depth model along Profile OBH-4 across the Xisha Trough, a seven-layer density–depth model is used to estimate density structure for the profile. The relationship between seismic velocity and radiogenic heat production is used to estimate the vertical distribution of heat sources in the lower crust. The 2-D temperature field is calculated by applying a 2-D numerical solution of the heat conduction equation and the thermal lithosphere thickness is obtained from the basalt dry solidus (BDS). The rheology of the profile is estimated on the basis of frictional failure in the brittle regime and power-law steady-state creep in the ductile regime. Rheological model is constructed for a three-layer model involving a granitic upper crust, a quartz diorite lower crust and an olivine upper mantle. Gravity modeling supports basically the velocity–depth model. The Moho along Profile OBH-4 is of relatively high heat flow ranging from 46 to 60 mW/m² and the Moho heat flow is higher in the trough than on the flanks. The depth of the “thermal” lithospheric lower boundary is about 54 km in the center, deepens toward two sides, and is about 75 km at the northern slope area and about 70 km at the southern Xisha–Zhongsha Block. Rheological calculation indicates that the two thinnest ductile layers in the crust and the thickest brittle layer in the uppermost mantle lie in the central region, showing that the Xisha Trough has been rheologically strengthened, which are mainly due to later thermal relaxation. In addition, the strengthening in rheology during rifting was not the main factor in hampering the breakup of the Xisha Trough.     

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.     

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.     

用ANSYS模拟结构-地基动力相互作用振动台试验的建模方法

本文以结构－地基动力相互作用振动台模型试验为基础，结合通用有限元软件ANSYS，对结构－地基动力相互作用体系进行有限元计算建模的一些问题作了研究，主要包括柔性土容器的模拟、土体动力本构模型的选用、土体与结构交界面上的状态非线性模型、网格划分、重力的考虑、结构中钢筋的处理以及对称性的应用等。文中给出了利用上述建模方法对结构－地基动力相互作用体系进行计算的一些加速度时程结果，并与试验结果相对照，吻合较好。通过计算分析，验证了简化处理方法的合理性和计算模型的可行性。     

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.     

基于OpenGL的多层三维空间实体叠加建模

多层三维空间实体普遍存在于现实生活中。借助于目前国际最广泛应用的三维图形接口OpenGL将其真实地显示出来，具有重要的现实意义。其中多层三维空间实体建模技术是显示过程的关键步鼻。在实现单层空间实体三维建模的基础上，利用OpenGL实现了多层三维空间实体叠加建模，叠加实现的过程及其属性数据的表示方法。     

不同电阻率测深方法对旁侧不均匀体的反映

利用三维有限元数值模拟,重点研究了旁侧三维不均匀体所产生的视电阻率的影响问题,展示了4种常规电阻率测深法在旁侧不均匀体上的异常分布特点;探讨了测深点正下方目标体异常与旁侧不均匀体异常的区别问题,指出旁侧不均匀体的影响是产生电阻率法假异常的主要原因之一.     

生态地理建模中的多尺度问题

本文在分析生态地理建模内涵的基础上,讨论了生态地理建模中的尺度转换问题、跨尺度相互作用问题、空间尺度与时间尺度的关联问题和多尺度数据处理问题.由于生态地理问题的非线性、生态环境的异质性和随机事件,简单的线性尺度转换方法远不能满足生态地理建模的要求.为了从根本上解决生态地理建模中的时空尺度问题,除需要运用微分几何学和等级理论等经典方法外,还需要引入格点生成法和网格计算等现代理论和技术手段.