Numerical modeling for analyzing thermal surface anomalies induced by underground coal fires

Coal seams burning underneath the surface are recognized all over the world and have drawn increasing public attention in the past years. Frequently, such fires are analyzed by detecting anomalies like increased exhaust gas concentrations and soil temperatures at the surface. A proper analysis presumes the understanding of involved processes, which determine the spatial distribution and dynamic behavior of the anomalies.In this paper, we explain the relevance of mechanical and energy transport processes with respect to the occurrence of temperature anomalies at the surface. Two approaches are presented, aiming to obtain insight into the underground coal fire situation: In-situ temperature mapping and numerical simulation. In 2000 to 2005, annual temperature mapping in the Wuda (Inner Mongolia, PR China) coal fire area showed that most thermal anomalies on the surface are closely related to fractures, where hot exhaust gases from the coal fire are released. Those fractures develop due to rock mechanical failure after volume reduction in the seams. The measured signals at the surface are therefore strongly affected by mechanical processes.More insight into causes and effects of involved energy transport processes is obtained by numerical simulation of the dynamic behavior of coal fires. Simulations show the inter-relation between release and transport of thermal energy in and around underground coal fires. Our simulation results show a time delay between the coal fire propagation and the observed appearance of the surface temperature signal. Additionally, the overall energy flux away from the burning coal seam into the surrounding bedrock is about 30-times higher than the flux through the surface. This is of particular importance for an estimation of the energy released based on surface temperature measurements. Finally, the simulation results also prove that a fire propagation rate estimated from the interpretation of surface anomalies can differ from the actual rate in the seam.     

Three-dimensional discontinuous deformation analysis (3-D DDA) using a new contact resolution algorithm

This paper presents a new contact calculating algorithm for contacts between two polyhedra with planar boundaries in the three-dimensional discontinuous deformation analysis (3-D DDA). In this algorithm, all six type contacts in 3-D (vertex-to-face, vertex-to-edge, vertex-to-vertex, face-to-face, edge-to-edge, and edge-to-face) are simply transformed into the form of point-to-face contacts. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. In this paper, formulations of normal contact, shear contact and frictional force submatrices based on the new method are derived and the algorithm has been programmed in VC⁺⁺. Examples are provided to demonstrate the new contact rule between two blocks.     

Application of MODFLOW and geographic information system to groundwater flow simulation in North China Plain, China

MODFLOW is a groundwater modeling program. It can be compiled and remedied according to the practical applications. Because of its structure and fixed data format, MODFLOW can be integrated with Geographic Information Systems (GIS) technology for water resource management. The North China Plain (NCP), which is the politic, economic and cultural center of China, is facing with water resources shortage and water pollution. Groundwater is the main water resource for industrial, agricultural and domestic usage. It is necessary to evaluate the groundwater resources of the NCP as an entire aquifer system. With the development of computer and internet information technology it is also necessary to integrate the groundwater model with the GIS technology. Because the geological and hydrogeological data in the NCP was mainly in MAPGIS format, the powerful function of GIS of disposing of and analyzing spatial data and computer languages such as Visual C and Visual Basic were used to define the relationship between the original data and model data. After analyzing the geological and hydrogeological conditions of the NCP, the groundwater flow numerical simulation modeling was constructed with MODFLOW. On the basis of GIS, a dynamic evaluation system for groundwater resources under the internet circumstance was completed. During the process of constructing the groundwater model, a water budget was analyzed, which showed a negative budget in the NCP. The simulation period was from 1 January 2002 to 31 December 2003. During this period, the total recharge of the groundwater system was 49,374 × 10⁶ m³ and the total discharge was 56,530 × 10⁶ m³ the budget deficit was −7,156 × 10⁶ m³. In this integrated system, the original data including graphs and attribution data could be stored in the database. When the process of evaluating and predicting groundwater flow was started, these data were transformed into files that the core program of MODFLOW could read. The calculated water level and drawdown could be displayed and reviewed online.     

On CO<Subscript>2</Subscript> fluid flow and heat transfer behavior in the subsurface,following leakage from a geologic storage reservoir

Geologic storage of CO₂ is expected to produce plumes of large areal extent, and some leakage may occur along fractures, fault zones, or improperly plugged pre-existing wellbores. A review of physical and chemical processes accompanying leakage suggests a potential for self-enhancement. The numerical simulations presented here confirm this expectation, but reveal self-limiting features as well. It seems unlikely that CO₂ leakage could trigger a high-energy run-away discharge, a so-called “pneumatic eruption,” but present understanding is insufficient to rule out this possibility. The most promising avenue for increasing understanding of CO₂ leakage behavior is the study of natural analogues.     

Determination of historical emission of heavy metals into the atmosphere: Celje case study

The research shows that in the Celje area (Slovenia), the historical anthropogenical emissions are 1,712 tons of Zn and 9.1 tons of Cd. For Zn, this value represents approximately 0.3% of the total Zn production in that area. Close to the former zinc smelting plant, the “Zn precipitation” has been estimated to be up to 0.036 mm. The 100-year Zn production left behind a heavily contaminated area with maximum concentrations of Zn of up to 5.6% in attic dust and 0.85% in the soil, and 456 mg/kg of Cd in attic dust and 59.1 mg/kg in the soil. The calculation of historical emissions is based on the data of heavy metals concentration in the attic dust at 98 sampling points and on the data from 19 measurement sites of the weight of total monthly air deposit. The main idea behind determining past emissions is that when the weight of the deposited dust on a small area is multiplied by the concentration of the element in that area, the mass of the polluter which has been transported to the place of interest by air can be calculated. If we sum up all the weight over the whole geochemical anomaly, we get the quantity of historical emissions.     

An operational multiscale system for hazards prediction, mapping, and response

By definition, a crisis is a situation that requires assistance to be managed. Hence, response to a crisis involves the merging of local and non-local emergency response personnel. In this situation, it is critical that each participant: (1) know the roles and responsibilities of each of the other participants; (2) know the capabilities of each of the participants; and (3) have a common basis for action. For many types of natural disasters, this entails having a common operational picture of the unfolding events, including detailed information on the weather, both current and forecasted, that may impact on either the emergency itself or on response activities. The Consequences Assessment Tool Set (CATS) is a comprehensive package of hazard prediction models and casualty and damage assessment tools that provides a linkage between a modeled or observed effect and the attendant consequences for populations, infrastructure, and resources, and, hence, provides the common operational picture for emergency response. The Operational Multiscale Environment model with Grid Adaptivity (OMEGA) is an atmospheric simulation system that links the latest methods in computational fluid dynamics and high-resolution gridding technologies with numerical weather prediction to provide specific weather analysis and forecast capability that can be merged into the geographic information system framework of CATS. This paper documents the problem of emergency response as an end-to-end system and presents the integrated CATS–OMEGA system as a prototype of such a system that has been used successfully in a number of different situations.     

Flood hazard delineation combining geomorphological and hydrological methods: an example in the Northern Iberian Peninsula

Flood mapping requires the combination and integration of geomorphological and hydrological-hydraulic methods; however, despite this, there is very little scientific literature that compares and validates both methods. Two types of analysis are addressed in the present article. On the one hand, maps of flood plains have been elaborated using geomorphological evidence and historical flood data in the mountainous area of northwestern Spain, covering an area of more then 232 km² of floodplains. On the other hand, a hydrometeorological model has been developed (Clark semidistributed unit hydrograph) in the Sarria River basin (155 km², NW Spain). This basin is not gauged, hence the model was subjected to a goodness-of-fit test of its parameter (curve number) by means of Monte Carlo simulation. The peak flows obtained by means of the hydrological model were used for hydraulic modeling (one-phase, one-dimensional and steady flow) in a 4 km² urban stretch of the river bed. The delineation of surface areas affected by floods since 1918, as well as those analyzed subsequent to the geomorphological study, reveals a high degree of reliability in the delineation of the flooded areas with frequent recurrence intervals (<50 years). If we compare these flooded surface areas with the estimate obtained by the hydrological-hydraulic method we can see that the latter method overestimates the extent of the surface water by 144% for very frequent recurrence intervals (>10 years) and underestimates it as the recurrence interval increases, by up to 80% less floodplain for exceptional events (>500 years). Finally, a management map is put forth combining the most reliable results available by integrating both methods. Originally presented at the Sixth International Conference on Geomorphology.     

Numerical modeling of tide-surge interaction along Orissa coast of India

The Orissa coast of India is one of the most vulnerable regions of extreme sea levels associated with severe tropical cyclones. There was extensive loss of life and property due to the October 1999 super cyclone, which devastated large part of the Orissa coast. The shallow nature of the head bay, presence of a large number of deltas formed by major rivers of Orissa such as Mahanadi and Dhamra, and high tidal range are responsible for storm surge flooding in the region. Specifically, rising and falling tidal phases influence the height, duration, and arrival time of peak surge along the coast. The objective of the present study is to evaluate the tide-surge interaction during the 1999 Orissa cyclone by using nonlinear vertically integrated numerical models. The pure tidal solution for the head bay region of the Bay of Bengal provides the initial condition for the fine resolution nested grid Orissa model. However, the feedback from the Orissa model does not affect the head bay model as the study provides a one-way interaction. Numerical experiments are performed to study the tide-surge interaction by considering various relative phases of the tidal waves with the surge-wave produced by 1999 Orissa cyclone. The comparison, although utilizing only the limited estimates of tidal data, appears adequate to assert that the principal features are reproduced correctly.     

Ground motion estimation during the Kashmir earthquake of 8th October 2005

In this article, analytical methods have been used to estimate ground motion during the 8 October 2005, Kashmir earthquake. Peak ground acceleration (PGA) values at several stations in the epicentral region have been estimated by empirical analytical source mechanism models. As an alternate analysis, PGA estimates have also been obtained using the stochastic finite fault seismological model. The estimated PGAs are compared with that obtained from damage values. A PGA contour map in the near-source region is provided. It is found that very near to the epicenter, PGA would have reached more than 1 g. It is demonstrated that empirical analytical models can be effectively used to estimate ground motion due to rupture of active faults.