Effects of geological inhomogeneity on high Rayleigh number steady state heat and mass transfer in fluid-saturated porous media heated from below

A parametric study is carried out to investigate how geological inhomogeneity affects the pore-fluid convective flow field, the temperature distribution, and the mass concentration distribution in a fluid-saturated porous medium. The related numerical results have demonstrated that (1) the effects of both medium permeability inhomogeneity and medium thermal conductivity inhomogeneity are significant on the pore-fluid convective flow and the species concentration distribution in the porous medium; (2) the effect of medium thermal conductivity inhomogeneity is dramatic on the temperature distribution in the porous medium, but the effect of medium permeability inhomogeneity on the temperature distribution may be considerable, depending on the Rayleigh number involved in the analysis; (3) if the coupling effect between pore-fluid flow and mass transport is weak, the effect of the Lewis number is negligible on the pore-fluid convective flow and temperature distribution, hut it is significant on the species concentration distribution in the medium.

Instantaneous outbursts in underground coal mines: An overview and association with coal type

Instantaneous outbursts in underground coal mines have occurred in at least 16 countries, involving both methane (CH4) and carbon dioxide (CO2). The precise mechanisms of an instantaneous outburst are still unresolved but must consider the effects of stress, gas content and physico-mechanical properties of the coal. Other factors such as mining methods (e.g., development heading into the coal seam) and geological features (e.g., coal seam disruptions from faulting) can combine to exacerbate the problem. Prediction techniques continue to be unreliable and unexpected outburst incidents resulting in fatalities are a major concern for underground coal operations. Gas content thresholds of 9 m(3)/t for CH4 and 6 m(3)/t for CO2 are used in the Sydney Basin, to indicate outburst-prone conditions, but are reviewed on an individual mine basis and in mixed as situations. Data on the sorption behaviour of Bowen Basin coals from Australia have provided an explanation for the conflicting results obtained by coal face desorption indices used for outburst-proneness assessment. A key factor appears to be different desorption rates displayed by banded coals, which is supported by both laboratory and mine-site investigations. Dull coal bands with high fusinite and semifusinite contents tend to display rapid desorption from solid coal, for a given pressure drop. The opposite is true for bright coal bands with high vitrinite contents and dull coal bands with high inertodetrinite contents. Consequently, when face samples of dull, fusinite-or semifusinite-rich coal of small particle size are taken for desorption testing, much gas has already escaped and low readings result. The converse applies for samples taken from coal bands with high vitrinite and/or inertodetrinite contents. In terms of outburst potential, it is the bright, vitrinite-rich and the dull, inertodetrinite-rich sections of a coal seam that appear to be more outburst-prone. This is due to the ability of the solid coal to retain gas, even after pressure reduction, creating a gas content gradient across the coal face sufficient to initiate an outburst. Once the particle size of the coal is reduced, rapid gas desorption can then take place. (C) 1998 Elsevier Science.

Type 1 and Type 2 cytokines: from basic science to fungal infections

At the present time, it is clear that Th1 responses afford protection against the fungi; however, the development, maintenance and function of the protective immune responses are complex mechanisms and are influenced by multiple factors. The route of infection has been shown to affect initial cytokine production and, consequently, the induction of protective Th1 responses. The ability of different isolates of the same fungal agent to induce and sustain a protective response has also been emphasized. Protective immune responses have been shown to vary in genetically different mouse strains after infection. In addition, these protective responses, such as cellular influx and cytokine production, also vary within the same animal depending on the tissue infected. The functional dominance of certain cytokines over others in influencing development and maintenance of protective responses has been discussed. Certain cytokines may act differently in hosts lacking important components of their innate or immune repertoire. It is evident from these presentations that a more comprehensive understanding of the protective mechanisms against different fungal agents is emerging. However, there is still much to learn before cytokine modulatory therapy can be used effectively without risk in the human host.

Factors controlling the origin of gas in Australian Bowen Basin coals

Open system pyrolysis (heating rate 10 degrees C/min) of coal maturity (vitrinite reflectance, VR) sequence (0.5%, 0.8% and 1.4% VR) demonstrates that there are two stages of thermogenic methane generation from Bowen Basin coals. The first and major stage shows a steady increase in methane generation maximising at 570 degrees C, corresponding to a VR of 2-2.5%. This is followed by a less intense methane generation which has not as yet maximised by 800 degrees C (equivalent to VR of 5%). Heavier (C2+) hydrocarbons are generated up to 570 degrees C after which only the C-1 (CH4, CO and CO2) gases are produced. The main phase of heavy hydrocarbon generation occurs between 420 and 510 degrees C. Over this temperature range,methane generation accounts for only a minor component, whereas the wet gases (C-2-C-5) are either in equal abundance or are more abundant by a factor of two than the liquid hydrocarbons. The yields of non-hydrocarbon gases CO2 and CO are greater then methane during the early stages of gas generation from an immature coal, subordinate to methane during the main phase of methane generation after which they are again dominant. Compositional data for desorbed and produced coal seam gases from the Bowen show that CO2 and wet gases are a minor component. This discrepancy between the proportion of wet gas components produced during open system pyrolysis and that observed in naturally matured coals may be the result of preferential migration of wet gas components, by dilution of methane generated during secondary cracking of bitumen, or kinetic effects associated with different activations for production of individual hydrocarbon gases. Extrapolation of results of artificial pyrolysis of the main organic components in coal to geological significant heating rates suggests that isotopically light methane to delta(13)C of -50 parts per thousand can be generated. Carbon isotope depletions in C-13 are further enhanced, however, as a result of trapping of gases over selected rank levels (instantaneous generation) which is a probable explanation for the range of delta(13)C values we have recorded in methane desorbed from Bowen Basin coals (-51 +/- 9 parts per thousand). Pervasive carbonate-rich veins in Bowen Basin coals are the product of magmatism-related hydrothermal activity. Furthermore, the pyrolysis results suggest an additional organic carbon source front CO2 released at any stage during the maturation history could mix in varying proportions with CO2 from the other sources. This interpretation is supported by C and O isotopic ratios, of carbonates that indicate mixing between magmatic and meteoric fluids. Also, the steep slope of the C and O isotope correlation trend suggests that the carbonates were deposited over a very narrow temperature interval basin-wide, or at relatively high temperatures (i.e., greater than 150 degrees C) where mineral-fluid oxygen isotope fractionations are small. These temperatures are high enough for catagenic production of methane and higher hydrocarbons from the coal and coal-derived bitumen. The results suggests that a combination of thermogenic generation of methane and thermodynamic processes associated with CH4/CO2 equilibria are the two most important factors that control the primary isotope and molecular composition of coal seam gases in the Bowen Basin. Biological process are regionally subordinate but may be locally significant. (C) 1998 Published by Elsevier Science Ltd. All rights reserved.

Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method. Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately. Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.

Theoretical and numerical analyses of convective instability in porous media with upward throughflow

Exact analytical solutions have been obtained for a hydrothermal system consisting of a horizontal porous layer with upward throughflow. The boundary conditions considered are constant temperature, constant pressure at the top, and constant vertical temperature gradient, constant Darcy velocity at the bottom of the layer. After deriving the exact analytical solutions, we examine the stability of the solutions using linear stability theory and the Galerkin method. It has been found that the exact solutions for such a hydrothermal system become unstable when the Rayleigh number of the system is equal to or greater than the corresponding critical Rayleigh number. For small and moderate Peclet numbers (Pe less than or equal to 6), an increase in upward throughflow destabilizes the convective flow in the horizontal layer. To confirm these findings, the finite element method with the progressive asymptotic approach procedure is used to compute the convective cells in such a hydrothermal system. Copyright (C) 1999 John Wiley & Sons, Ltd.

Numerical and experimental study of seepage in unconfined aquifers with a periodic boundary condition

The assessment of groundwater conditions within an unconfined aquifer with a periodic boundary condition is of interest in many hydrological and environmental problems. A two-dimensional numerical model for density dependent variably saturated groundwater flow, SUTRA (Voss, C.I., 1984. SUTRA: a finite element simulation model for saturated-unsaturated, fluid-density dependent ground-water flow with energy transport or chemically reactive single species solute transport. US Geological Survey, National Center, Reston, VA) is modified in order to be able to simulate the groundwater flow in unconfined aquifers affected by a periodic boundary condition. The basic flow equation is changed from pressure-form to mixed-form. The model is also adjusted to handle a seepage-face boundary condition. Experiments are conducted to provide data for the groundwater response to the periodic boundary condition for aquifers with both vertical and sloping faces. The performance of the numerical model is assessed using those data. The results of pressure- and mixed-form approximations are compared and the improvement achieved through the mixed-form of the equation is demonstrated. The ability of the numerical model to simulate the water table and seepage-face is tested by modelling some published experimental data. Finally the numerical model is successfully verified against present experimental results to confirm its ability to simulate complex boundary conditions like the periodic head and the seepage-face boundary condition on the sloping face. (C) 1999 Elsevier Science B.V. All rights reserved.

Petrographic and isotope constraints on the origin of authigenic carbonate minerals and the associated fluid evolution in Late Permian coal measures, Bowen Basin (Queensland), Australia

Authigenic carbonate minerals are ubiquitous throughout the Late Permian coal measures of the Bowen Basin, Queensland, Australia. In the northern Bowen Basin, carbonates include the following assemblages: siderite I (delta O-18(SMOW) = +11.4 to + 17%, delta C-13(PDB) = - 5.3 to + 120), Fe-Mg calcite-ankerite-siderite II mineral association (delta O-18(SMOW) = +7.2 to + 10.20, delta C-13(PDB) = 10.9 to - 1.80 for ankerite) and a later calcite (delta O-18(SMOW) = +5.9 to + 14.60, delta C-13(PDB) = -11.4 to + 4.40). In the southern Bowen Basin, the carbonate phase consists only of calcite (delta O-18(SMOW) = +12.5 to + 14.80, delta C-13(PDB) = -19.4 to + 0.80), where it occurs extensively throughout all stratigraphic levels. Siderite I occurs in mudrocks and sandstones and predates all other carbonate minerals. This carbonate phase is interpreted to have formed as an early diagenetic mineral from meteoric waters under cold climate and reducing conditions. Fe-Mg calcite-ankerite-siderite Il occur in sandstones as replacement of volcanic rock fragments. Clay minerals (illite-smectite, chlorite and kaolinite) postdate Ca-Fe-Mg carbonates, and precipitation of the later calcite is associated with clay mineral formation. The Ca-Fe-Mg carbonates and later calcite of the northern Bowen Basin are regarded as having formed as a result of hydrothermal activity during the latest Triassic extensional tectonic event which affected this part of the basin, rather than deep burial diagenesis during the Middle to Late Triassic as previously reported. This hypothesis is based on the timing relationships of the authigenic mineral phases and the low delta O-18 values of ankerite and calcite, together with radiometric dating of illitic clays and recently published regional geological evidence. Following the precipitation of the Ca-Fe-Mg carbonates from strongly O-18-depleted meteoric-hydrothermal fluids, continuing fluid circulation and water-rock interaction resulted in dissolution of these carbonate phases as well as labile fragments of volcaniclastic rocks. Subsequently, the later calcite and day minerals precipitated from relatively evolved (O-18-enriched) fluids. The nearly uniform delta O-18 values of the southern Bowen Basin calcite have been attributed to very low water/rock ratio in the system, where the fluid isotropic composition was buffered by the delta O-18 values of rocks. (C) 2000 Elsevier Science B.V. All rights reserved.

Wave-induced seepage flux into anisotropic seabeds

Considerable effort has been devoted to quantifying the wave-induced soil response in a porous seabed in the last few decades. Most previous investigations have focused on the analysis of pore pressure and effective stresses within isotropic sediments, despite strong evidence of anisotropic soil behaviour reported in the literature. Furthermore, the seepage flux, which is important in the context of contaminant transport, has not been examined. In this paper, we focus on water wave-driven seepage in anisotropic marine sediments of finite thickness. The numerical results predict that the effects of hydraulic anisotropy and anisotropic soil behaviour on the wave-driven seepage in marine sediment are significant. Copyright (C) 2001 John Wiley & Sons, Ltd.

A numerical study of flexural buckling of foliated rock slopes

The occurrence of foliated rock masses is common in mining environment. Methods employing continuum approximation in describing the deformation of such rock masses possess a clear advantage over methods where each rock layer and each inter-layer interface (joint) is explicitly modelled. In devising such a continuum model it is imperative that moment (couple) stresses and internal rotations associated with the bending of the rock layers be properly incorporated in the model formulation. Such an approach will lead to a Cosserat-type theory. In the present model, the behaviour of the intact rock layer is assumed to be linearly elastic and the joints are assumed to be elastic-perfectly plastic. Condition of slip at the interfaces are determined by a Mohr-Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformations. The model is incorporated into the finite element program AFENA and validated against an analytical solution of elementary buckling problems of a layered medium under gravity loading. A design chart suitable for assessing the stability of slopes in foliated rock masses against flexural buckling failure has been developed. The design chart is easy to use and provides a quick estimate of critical loading factors for slopes in foliated rock masses. It is shown that the model based on Euler's buckling theory as proposed by Cavers (Rock Mechanics and Rock Engineering 1981; 14:87-104) substantially overestimates the critical heights for a vertical slope and underestimates the same for sub-vertical slopes. Copyright (C) 2001 John Wiley & Sons, Ltd.

Analysis of steady-state heat transfer through mid-crustal vertical cracks with upward throughflow in hydrothermal systems

We conduct a theoretical analysis of steady-state heat transfer problems through mid-crustal vertical cracks with upward throughflow in hydrothermal systems. In particular, we derive analytical solutions for both the far field and near field of the system. In order to investigate the contribution of the forced advection to the total temperature of the system, two concepts, namely the critical Peclet number and the critical permeability of the system, have been presented and discussed in this paper. The analytical solution for the far field of the system indicates that if the pore-fluid pressure gradient in the crust is lithostatic, the critical permeability of the system can be used to determine whether or not the contribution of the forced advection to the total temperature of the system is negligible. Otherwise, the critical Peclet number should be used. For a crust of moderate thickness, the critical permeability is of the order of magnitude of 10(-20) m(2), under which heat conduction is the overwhelming mechanism to transfer heat energy, even though the pore-fluid pressure gradient in the crust is lithostatic. Furthermore, the lower bound analytical solution for the near field of the system demonstrates that the permeable vertical cracks in the middle crust can efficiently transfer heat energy from the lower crust to the upper crust of the Earth. Copyright (C) 2002 John Wiley Sons, Ltd.

An equivalent algorithm for simulating thermal effects of magma intrusion problems in porous rocks

An equivalent algorithm is proposed to simulate thermal effects of the magma intrusion in geological systems, which are composed of porous rocks. Based on the physical and mathematical equivalence, the original magma solidification problem with a moving boundary between the rock and intruded magma is transformed into a new problem without the moving boundary but with a physically equivalent heat source. From the analysis of an ideal solidification model, the physically equivalent heat source has been determined in this paper. The major advantage in using the proposed equivalent algorithm is that the fixed finite element mesh with a variable integration time step can be employed to simulate the thermal effect of the intruded magma solidification using the conventional finite element method. The related numerical results have demonstrated the correctness and usefulness of the proposed equivalent algorithm for simulating the thermal effect of the intruded magma solidification in geological systems. (C) 2003 Elsevier B.V. All rights reserved.

Convective instability of 3-D fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the convective instability of 3-D fluid-saturated geological fault zones when they are heated uniformly from below. In particular, we have derived exact analytical solutions for the critical Rayleigh numbers of different convective flow structures. Using these critical Rayleigh numbers, three interesting convective flow structures have been identified in a geological fault zone system. It has been recognized that the critical Rayleigh numbers of the system have a minimum value only for the fault zone of infinite length, in which the corresponding convective flow structure is a 2-D slender-circle flow. However, if the length of the fault zone is finite, the convective flow in the system must be 3-D. Even if the length of the fault zone is infinite, since the minimum critical Rayleigh number for the 2-D slender-circle flow structure is so close to that for the 3-D convective flow structure, the system may have almost the same chance to pick up the 3-D convective flow structures. Also, because the convection modes are so close for the 3-D convective flow structures, the convective flow may evolve into the 3-D finger-like structures, especially for the case of the fault thickness to height ratio approaching zero. This understanding demonstrates the beautiful aspects of the present analytical solution for the convective instability of 3-D geological fault zones, because the present analytical solution is valid for any value of the ratio of the fault height to thickness. Using the present analytical solution, the conditions, under which different convective flow structures may take place, can be easily determined.

Mediatizing educational policy: The journalistic field, science policy, and cross-field effects

This paper is concerned to demonstrate the usefulness of the theory of Bourdieu, including the concepts of field, logics of practice and habitus, to understanding relationships between media and policy, what Fairclough has called the 'mediatization' of policy. Specifically, the paper draws upon Bourdieu's accessible account of the journalistic field as outlined in On television and journalism. The usefulness of this work is illustrated through a case study of a recent Australian science policy, The chance to change. As this policy went through various iterations and media representations, its naming and structure became more aphoristic. This is the mediatization of contemporary policy, which often results in policy as sound bite. The case study also shows the cross-field effects of this policy in education, illustrating how today educational policy can be spawned from developments in other public policy fields.