Numerical simulation of a stratigraphic model     

Véronique Gervais  Roland Masson 《Computational Geosciences》2008,12(2):163-179

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.  相似文献   

Structural evolution of the Bayanhongor region, west-central Mongolia     

Soichi Osozawa  Gonchig Tsolmon  Uondon Majigsuren  Jargalan Sereenen  Shigeaki Niitsuma  Naoyoshi Iwata  Terry Pavlis  Bor-ming Jahn 《Journal of Asian Earth Sciences》2008,33(5-6):337-352

Most of previous models suggest that the Central Asia Orogenic Belt grew southward in the Phanerozoic. However, in the Bayanhongor region in west-central Mongolia, volcanic arc, accretionary prism, ophiolite, and passive margin complexes accreted northeastward away from the Baydrag micro-continent, and hence the region constitutes the southwestern part of a crustal-scale syntaxis close to the west. The syntaxis should be original, because presumably reorientation due to strike-slip faulting can be ignored. It is reconfirmed that the Baydrag eventually collided with another micro-continent (the Hangai) to the northeast. A thick sedimentary basin developed along the southern passive margin of the Hangai micro-continent. This region is also characterized by an exhumed metamorphosed accretionary complex and a passive margin complex, which are both bounded by detachment faults as well as basal reverse faults which formed simultaneously as extrusion wedges. This part of the Central Asia Orogenic Belt lacks exhumed crystalline rocks as observed in the Himalayas and other major collisional orogenic belts. In addition, we identified two phases of deformation, which occurred at each phase of zonal accretion as D1 through Cambrian and Devonian, and a synchronous phase of final micro-continental collision of Devonian as D2. The pre-collisional ocean was wide enough to be characterized by a mid-ocean ridge and ocean islands. Two different structural trends of D1 and D2 are observed in accretionary complexes formed to the southwest of the late Cambrian mid-ocean ridge. That is, the relative plate motions on both sides of the mid-ocean ridge were different. Accretionary complexes and passive margin sediments to the northeast of the mid-ocean ridge also experienced two periods of deformation but show the same structural trend. Unmetamorphosed cover sediments on the accretionary prism and on the Hangai micro-continent experienced only the D2 event due to micro-continental collision. These unmetamorphosed sediments form the hanging walls of the detachment faults. Moreover, they were at least partly derived from an active volcanic arc formed at the margin of the Baydrag micro-continent.  相似文献   

Use of unit response functions for management of regional multilayered aquifers: application to the North Aquitaine Tertiary system (France)     

Francois Larroque  Wiktor Treichel  Alain Dupuy 《Hydrogeology Journal》2008,16(2):215-233

Although there are many numerical applications used in aquifer management, few methodologies seem to be sufficiently adapted to provide solutions for regional-scale problems. Their applications for multilayered aquifers are also too limited. A new groundwater management tool adapted to such large physical systems has been developed, using the unit response function (URF) approach. This tool is based on the link between a numerical groundwater flow model and a mathematical optimisation tool. Pre- and post processing are based on a geographical information system (GIS). The developed tool has been applied to a real management case—the regional Aquitaine multilayered aquifer (France). A representative hydrodynamic model was built using MODFLOW 2000 code. First tests show that the linearity condition needed for the application of URF methodology is respected despite the high complexity of the system. Two management scenarios were tested. Optimal solutions thus obtained allow for viable exploitation alternatives to be proposed, which integrate complex environmental constraints.

---

Résumé Malgré un nombre important d’applications existant dans le domaine de la gestion des aquifères, peu de méthodologies semblent adaptées pour proposer des solutions à des problèmes réels de dimension régionale. Leurs application pour des systèmes multicouches sont souvent trop limitées. Un outil de gestion adapté à ces systèmes de large taille a été développé, basé sur la théorie des Fonctions de Réponses Unitaires (FRU). Cet outil est basé sur le couplage d’un modèle numérique d’écoulement et d’un code d’optimisation mathématique. Les entrées et sorties du modèle sont réalisées à l’aide d’un Système d’Information Géographique (SIG). L’outil développé est appliqué à un cas de gestion réel—le système aquifère régional multicouche aquitain (France). Un modèle hydrodynamique représentatif a été construit, basé sur le code MODFLOW 2000. Les premiers tests montrent que les conditions de linéarité nécessaires à l’application de la méthode FRU sont respectées, malgré l’importante complexité du système. Deux scénarii de gestion sont testés. Les solutions optimales ainsi obtenues permettent de proposer des solutions d’exploitation alternatives viables, intégrant des contraintes environnementales complexes.

---

Resumen Aunque hay muchas aplicaciones numéricas usadas en la gestión de acuíferos, unas pocas metodologías parecen estar adaptadas suficientemente para dar las soluciones a los problemas de escala regional. También son limitadas sus aplicaciones para los acuíferos multicapa. Se ha desarrollado una nueva herramienta para el manejo del agua subterránea, adaptada a tales sistemas físicos grandes, usando el avance de la Función de Respuesta Unitaria (URF). Esta herramienta se basa en la combinación entre un modelo numérico de flujo de agua subterránea y una herramienta de optimización matemática. El pre-proceso y el post-proceso se basan en un Sistema de Información Geográfico (GIS). La herramienta desarrollada se ha aplicado a un caso real de manejo en el acuífero regional multicapa de Aquitania (Francia). Se construyó un modelo hidrodinámico representativo, que usó el código MODFLOW 2000. Las primeras pruebas muestran que la condición de linearidad, necesaria para la aplicación de la metodología URF, se respeta a pesar de la complejidad alta del sistema. Se probaron dos escenarios de manejo. Las soluciones óptimas así obtenidas permiten la proyección de alternativas viables de explotación, las cuales integran restricciones medioambientales complejas.

  相似文献   

Evaluation of distributed recharge in an upland semi-arid karst system: the West Bank Mountain Aquifer,Middle East   总被引：2，自引：0，他引：2  

A. G. Hughes  M. M. Mansour  N. S. Robins 《Hydrogeology Journal》2008,16(5):845-854

Assessment of recharge in a structurally complex upland karst limestone aquifer situated in a semi-arid environment is difficult. Resort to surrogate indicators such as measurement of spring outflow and borehole discharge, is a common alternative, and attempts to apply conventional soil moisture deficit analysis may not adequately account for the intermittent spate conditions that arise in such environments. A modelling approach has been made using the West Bank Mountain Aquifer system in the Middle East as a trial. The model uses object oriented software which allows various objects to be switched on and off. Each of the main recharge processes identified in the West Bank is incorporated. The model allows either conventional soil moisture deficit analysis calculations or wetting threshold calculations to be made as appropriate, and accommodates both direct recharge and secondary recharge. Daily time steps enable recharge and runoff routing to be calculated for each node. Model runs have enabled a series of simulations for each of the three aquifer basins in the West Bank and for the whole of the West Bank. These provide recharge estimates comparable to those prepared by earlier workers by conventional means. The model is adaptable and has been successfully used in other environments.  相似文献   

Estimating hydraulic conductivity using grain-size analyses, aquifer tests, and numerical modeling in a riverside alluvial system in South Korea   总被引：2，自引：0，他引：2  

Jae-Yeol Cheong  Se-Yeong Hamm  Hyoung-Soo Kim  Eun-Joung Ko  Kyounghee Yang  Jeong-Hwan Lee 《Hydrogeology Journal》2008,16(6):1129-1143

Hydraulic conductivity (K) for an alluvial system in a riverbank filtration area in Changwon City, South Korea, has been studied using grain-size distribution, pumping and slug tests, and numerical modeling. The alluvial system is composed of layers: upper fine sand, medium sand, lower fine sand, and a highly conductive sand/gravel layer at the base. The geometric mean of K for the sand/gravel layer (9.89?×?10^?4 m s^?1), as determined by grain-size analyses, was 3.33 times greater than the geometric mean obtained from pumping tests (2.97?×?10^?4 m s^?1). The geometric mean of K estimates obtained from slug tests (3.08?×?10^?6 m s^?1) was one to two orders of magnitude lower than that from pumping tests and grain-size analyses. K estimates derived from a numerical model were compared to those derived from the grain-size methods, slug tests and pumping tests in order to determine the degree of deviation from the numerical model. It is considered that the K estimates determined by the slug tests resemble the uppermost part of the alluvial deposit, whereas the K estimates obtained by grain-size analyses and pumping tests are similar to those from the numerical model for the sand/gravel layer of the riverside alluvial system.  相似文献   

Pumping decisions for sustainable development of groundwater resources in areas of grassland degradation: a case study in Lanqi Banner, Inner Mongolia, China     

Bagen Chaolun  Bin He  Ruizhong Gao  Liang Wang 《Hydrogeology Journal》2008,16(6):1101-1111

The best planting alternatives for satisfying high water use demands of forage and fodder crops in a region of Inner Mongolia, China, were determined by a multiobjective distributed-parameter groundwater management model. These alternatives took account of different cropping patterns and pumping decisions associated with both temporal and spatial aspects of water allocation. The model was developed for phreatic, homogenous, and isotropic aquifers using the response matrix technique of quadratic programming theory and, in this case, using the alternative direction implicit (ADI) scheme. Model solutions using effective rainfall with a probability of 50%, show that average water table drawdown in the planning period (2006–2017) is 0.22 m and the groundwater fluctuation in each pumping well is very low. In order to evaluate the pumping decisions under an effective rainfall with a probability of 75%, a sensitive analysis was also conducted. Analysis shows that it is useful to apply the results from the proposed model to control the landscape degradation due to overgrazing and overpumping activities.  相似文献   

The Glueckstadt Graben of the North-German Basin: new insights into the structure from 3D and 2D gravity analyses   总被引：1，自引：0，他引：1  

Tamara Yegorova  Yuriy Maystrenko  Ulf Bayer  Magdalena Scheck-Wenderoth 《International Journal of Earth Sciences》2008,97(5):915-930

The structure of the Glueckstadt Graben has been investigated by use of 3D gravity backstripping technique and by 2D gravity and magnetic modelling. Subtracting the gravity effects of the Meso-Cenozoic sediments together with Permian salt reveals a positive residual anomaly within the Glueckstadt Graben. This anomaly includes two local maxima over the Westholstein and Eastholstein Troughs. The 2D gravity models point to the presence of a high-density body within the lower crust of the Glueckstadt Graben. In addition, the results of 2D magnetic modelling indicate that the central part of the high-density body is overlain by an area with high susceptibility. Most probable, the formation of this high-density body is a result of complex poly-phase tectonic history of the study area. Finally, the results of gravity modelling indicate that Permian salt is not homogeneous. 3D gravity analysis and, especially, 2D gravity modelling have distinguished the differences in degree of salt saturation in salt-rich bodies, and elucidate the proportion of Rotliegend salt.  相似文献   

A virtual voyage through 3D structures generated by chaotic mixing of magmas and numerical simulations: a new approach for understanding spatial and temporal complexity of magma dynamics     

Diego Perugini  Maurizio Petrelli  Giampiero Poli 《Visual Geosciences》2008,13(1):1-24

In this contribution, we present a virtual voyage through 3D structures generated by chaotic mixing of magmas and numerical simulations with the aim to highlight the power of 3D representations in the understanding of this geological phenomenon. In particular, samples of mixed juveniles from Salina island (Southern Italy) are reconstructed in 3D by serial lapping and digital montage and numerical simulations are performed by using a 3D chaotic dynamical system. Natural and simulated magma mixing structures are visualized by using several multimedia tools including animations and “virtual reality” models. It is shown that magma interaction processes can generate large spatial and temporal compositional heterogeneities in magmatic systems. The same topological structures are observed in both 3D reconstructed rock samples and chaotic numerical simulations, indicating that the mixing of magmas is governed by chaotic dynamics. The use of 3D multimedia models gives the opportunity to penetrate into magma mixing structures and to understand their significance in the context of magma dynamics. Such an approach is very powerful since multimedia tools can strongly capture the attention of the reader bringing him/her into an interactive and memorable geological experience. Electronic supplementary material  enclosed:  相似文献   

Comparison of numerical models of two debris flows in the Cortina d’ Ampezzo area,Dolomites, Italy   总被引：2，自引：0，他引：2  

Maria Concetta Armento  Rinaldo Genevois  Pia Rosella Tecca 《Landslides》2008,5(1):143-150

The accurate prediction of runout distances, velocities and the knowledge of flow rheology can reduce the casualties and property damage produced by debris flows, providing a means to delineate hazard areas, to estimate hazard intensities for input into risk studies and to provide parameters for the design of protective measures. The application of most of models that describe the propagation and deposition of debris flow requires detailed topography, rheological and hydrological data that are not always available for the debris-flow hazard delineation and estimation. In the Cortina d’Ampezzo area, Eastern Dolomites, Italy, most of the slope instabilities are represented by debris flows; 325 debris-flow prone watersheds have been mapped in the geomorphological hazard map of this area. We compared the results of simulations of two well-documented debris flows in the Cortina d’Ampezzo area, carried on with two different single-phase, non-Newtonian models, the one-dimensional DAN-W and the two-dimensional FLO-2D, to test the possibility to simulate the dynamic behaviour of a debris flow with a model using a limited range of input parameters. FLO-2D model creates a more accurate representation of the hazard area in terms of flooded area, but the results in terms of runout distances and deposits thickness are similar to DAN-W results. Using DAN-W, the most appropriate rheology to describe the debris-flow behaviour is the Voellmy model. When detailed topographical, rheological and hydrological data are not available, DAN-W, which requires less detailed data, is a valuable tool to predict debris-flow hazard. Parameters obtained through back-analysis with both models can be applied to predict hazard in other areas characterized by similar geology, morphology and climate.  相似文献   

Groundwater flow modeling for effective implementation of landslide stability enhancement measures     

H. K. Shrestha  R. Yatabe  N. P. Bhandary 《Landslides》2008,5(3):281-290

This paper deals with groundwater hydrology at a prominent fracture zone landslide slope (Nuta–Yone landslides) in Japan with an objective to explore an efficient method for the application of landslide stability enhancement measures. The correlation analyses between the hydrological parameters and ground surface movement data at this landslide resulted in low correlation values indicating that the geological formation of the area is extremely complex. For the purpose of understanding the groundwater flow behavior in the landslide area, a three-dimensional transient groundwater flow model was prepared for a part of the landslide slope, where the levels of effectiveness of applied landslide stability enhancement measures (in the form of multilayered deep horizontal drains) are different, and was calibrated against the measured water surface elevations at different piezometer locations. The parameter distributions in the calibrated model and the general directions of the groundwater flow in terms of flow vectors and the results of particle tracking at the model site were interpreted to understand the reasons for variations in effectiveness of existing landslide stability enhancement measures and to find potentially better locations for the implementation of future landslide stability enhancement measures. From the modeling results, it was also understood that groundwater flow model can be effectively used in better planning and locating the landslide stability enhancement measures.  相似文献