Review of recent results on excursion set models

Many images consist of two or more 'phases', where a phase is a collection of homogeneous zones. For example, the phases may represent the presence of different sulphides in an ore sample. Frequently, these phases exhibit very little structure, though all connected components of a given phase may be similar in some sense. As a consequence, random set models are commonly used to model such images. The Boolean model and models derived from the Boolean model are often chosen. An alternative approach to modelling such images is to use the excursion sets of random fields to model each phase. In this paper, the properties of excursion sets will be firstly discussed in terms of modelling binary images. Ways of extending these models to multi-phase images will then be explored. A desirable feature of any model is to be able to fit it to data reasonably well. Different methods for fitting random set models based on excursion sets will be presented and some of the difficulties with these methods will be discussed.

Finite element modelling of temperature gradient driven rock alteration and mineralization in porous rock masses

We present finite element simulations of temperature gradient driven rock alteration and mineralization in fluid saturated porous rock masses. In particular, we explore the significance of production/annihilation terms in the mass balance equations and the dependence of the spatial patterns of rock alteration upon the ratio of the roll over time of large scale convection cells to the relaxation time of the chemical reactions. Special concepts such as the gradient reaction criterion or rock alteration index (RAI) are discussed in light of the present, more general theory. In order to validate the finite element simulation, we derive an analytical solution for the rock alteration index of a benchmark problem on a two-dimensional rectangular domain. Since the geometry and boundary conditions of the benchmark problem can be easily and exactly modelled, the analytical solution is also useful for validating other numerical methods, such as the finite difference method and the boundary element method, when they are used to dear with this kind of problem. Finally, the potential of the theory is illustrated by means of finite element studies related to coupled flow problems in materially homogeneous and inhomogeneous porous rock masses. (C) 1998 Elsevier Science S.A. All rights reserved.

A numerical study of pore-fluid, thermal and mass flow iu fluid-saturated porous rock basins

We present a numerical methodology for the study of convective pore-fluid, thermal and mass flow in fluid-saturated porous rock basins. lit particular, we investigate the occurrence and distribution pattern of temperature gradient driven convective pore-fluid flow and hydrocarbon transport in the Australian North West Shelf basin. The related numerical results have demonstrated that: (1) The finite element method combined with the progressive asymptotic approach procedure is a useful tool for dealing with temperature gradient driven pore-fluid flow and mass transport in fluid-saturated hydrothermal basins; (2) Convective pore-fluid flow generally becomes focused in more permeable layers, especially when the layers are thick enough to accommodate the appropriate convective cells; (3) Large dislocation of strata has a significant influence off the distribution patterns of convective pore;fluid flow, thermal flow and hydrocarbon transport in the North West Shelf basin; (4) As a direct consequence of the formation of convective pore-fluid cells, the hydrocarbon concentration is highly localized in the range bounded by two major faults in the basin.

Sedimentologic, petrographic, and sulfur isotope constraints on fine-grained pyrite formation at Mount Isa Mine and environs, Northwest Queensland, Australia

Fine-grained pyrite is the earliest generation of pyrite and the most abundant sulfide within the Urquhart Shale at Mount Isa, northwest Queensland. The pyrite is intimately interbanded with ore-grade Pb-Zn miner alization at the Mount Isa mine but is also abundant north and south of the mine at several stratigraphic horizons within the Urquhart Shale. Detailed sedimentologic, petrographic, and sulfur isotope studies of the Urquhart Shale, mostly north of the mine, reveal that the fine-grained pyrite (delta(34)S = -3.3 to +26.3 parts per thousand) formed by thermochemical sulfate reduction during diagenesis. The sulfate source was local sulfate evaporites, pseudo morphs of which are present throughout the Urquhart Shale (i.e., gypsum, anhydrite, and barite). Deep-burial diagenetic replacement of these evaporites resulted in sulfate-bearing ground waters which migrated parallel to bedding. Fine-grained pyrite formed where these fluids infiltrated and then interacted with carbon-rich laminated siltstones. Comparison of the sulfur isotope systematics of fine-grained pyrite and spatially associated base metal sulfides from the Mount Isa Pb-Zn and Cu orebodies indicates a common sulfur source of ultimately marine origin for all sulfide types. Different sulfur isotope ratio distributions for the various sulfides are the result of contrasting formation mechanisms and/or depositional conditions rather than differing sulfur sources. The sulfur isotope systematics of the base metal and associated iron sulfide generations are consistent with mineralization by reduced hydrothermal fluids, perhaps generated by bulk reduction of evaporite-sourced sulfate-bearing waters generated deeper in the Mount Isa Group, the sedimentary sequence which contains the Urquhart Shale. The available sulfur isotope data from the Mount Isa orebodies are consistent with either a chemically and thermally zoned, evolving Cu-Pb-Zn system, or discrete Cu and Pb-Zn mineralizing events linked by a common sulfur source.

A review of the Arhythmacanthidae (Acanthocephala) with a description of Heterosentis hirsutus n. sp. from Cnidoglanis macrocephala (Plotosidae) in Australia

Heterosentis hirsutus n. sp. is described from Cnidoglanis macrocephalo (Siluriformes: Plotosidae) from the Swan Estuary, Western Australia. It is distinguished by having 14 longitudinal rows of 6-7 hooks per row on the proboscis, a trunk armed anteriorly and posteriorly (=genital spines) with minute spines and lemnisci that may extend to the poster;or margin of the proboscis receptacle The new species also has prominent fragmented nuclei in its trunk well. New information is given for Heterosentis plotosi Yamoguti, 1935 from Plotosus lineatus (Siluriformes: Plotosidae) and H. poraplagusiarum (Nickol, 1972) Amin, 1985 from Paraplogusia guttata (Pleuronectiformes: Cynoglossidoe), both from Queensland. A key to the species of Heterosentis Van Cleave, 1931 is provided. The Arhythmacanthidae subfamilies are reviewed: there is little utility in the recognition of these taxa because of the small number of genera involved and the validity/ of the characters on which they ore based is in doubt, particularly whether trunk spines are present or absent. Only Acanthocephaloides Meyer, 1932, Breizocanthus Golvon, 1969, Euzetocanthus Golvan & Houin, 1964, Heterosentis, Hypoechinorhynchus Yamaguti, 1939 and Paracanthocepholoides Golvan, 1969 of the Arhythmacanthidae are considered valid. A key to these genera is provided. The monotypic genus Neocanthocepholoides Cable & Quick, 1954 is considered a new synonym of Acanthocephaloides thus creating Acanthocephaloides spinicaudatus (Cable & Quick, 1954) n. comb. Arhythmocanthus Yamaguti, 1935 is maintained as a synonym of Heterosentis because the distinction between two and three hook types is made equivocal when the transition between the opical and subapical hooks is gradual.

Computer simulations of coupled problems in geological and geochemical systems

The finite element method is used to simulate coupled problems, which describe the related physical and chemical processes of ore body formation and mineralization, in geological and geochemical systems. The main purpose of this paper is to illustrate some simulation results for different types of modelling problems in pore-fluid saturated rock masses. The aims of the simulation results presented in this paper are: (1) getting a better understanding of the processes and mechanisms of ore body formation and mineralization in the upper crust of the Earth; (2) demonstrating the usefulness and applicability of the finite element method in dealing with a wide range of coupled problems in geological and geochemical systems; (3) qualitatively establishing a set of showcase problems, against which any numerical method and computer package can be reasonably validated. (C) 2002 Published by Elsevier Science B.V.

Finite element modelling of reactive fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins

We present the finite element simulations of reactive mineral carrying fluids mixing and mineralization in pore-fluid saturated hydrothermal/sedimentary basins. In particular we explore the mixing of reactive sulfide and sulfate fluids and the relevant patterns of mineralization for Load, zinc and iron minerals in the regime of temperature-gradient-driven convective flow. Since the mineralization and ore body formation may last quite a long period of time in a hydrothermal basin, it is commonly assumed that, in the geochemistry, the solutions of minerals are in an equilibrium state or near an equilibrium state. Therefore, the mineralization rate of a particular kind of mineral can be expressed as the product of the pore-fluid velocity and the equilibrium concentration of this particular kind of mineral Using the present mineralization rate of a mineral, the potential of the modern mineralization theory is illustrated by means of finite element studies related to reactive mineral-carrying fluids mixing problems in materially homogeneous and inhomogeneous porous rock basins.

Effects of hot intrusions on pore-fluid flow and heat transfer in fluid-saturated rocks

We use the finite element method to solve coupled problems between pore-fluid flow and heat transfer in fluid-saturated porous rocks. In particular, we investigate the effects of both the hot pluton intrusion and topographically driven horizontal flow on the distributions of the pore-flow velocity and temperature in large-scale hydrothermal systems. Since general mineralization patterns are strongly dependent on distributions of both the pore-fluid velocity and temperature fields, the modern mineralization theory has been used to predict the general mineralization patterns in several realistic hydrothermal systems. The related numerical results have demonstrated that: (1) The existence of a hot intrusion can cause an increase in the maximum value of the pore-fluid velocity in the hydrothermal system. (2) The permeability of an intruded pluton is one of the sensitive parameters to control the pore-fluid flow, heat transfer and ore body formation in hydrothermal systems. (3) The maximum value of the pore-fluid velocity increases when the bottom temperature of the hydrothermal system is increased. (4) The topographically driven flow has significant effects on the pore-fluid flow, temperature distribution and precipitation pattern of minerals in hydrothermal systems. (5) The size of the computational domain may have some effects on the pore-fluid flow and heat transfer, indicating that the size of a hydrothermal system may affect the pore-fluid flow and heat transfer within the system. (C) 2003 Elsevier Science B.V. All rights reserved.

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.