Prediction of rock alteration patterns: A potential tool in mineral exploration

Hydrothermal alteration of a quartz‐K‐feldspar rock is simulated numerically by coupling fluid flow and chemical reactions. Introduction of CO₂ gas generates an acidic fluid and produces secondary quartz, muscovite and/or pyrophyllite at constant temperature and pressure of 300°C and 200 MPa. The precipitation and/or dissolution of the secondary minerals is controlled by either mass‐action relations or rate laws. In our simulations the mass of the primary elements are conserved and the mass‐balance equations are solved sequentially using an implicit scheme in a finite‐element code. The pore‐fluid velocity is assumed to be constant. The change of rock volume due to the dissolution or precipitation of the minerals, which is directly related to their molar volume, is taken into account. Feedback into the rock porosity and the reaction rates is included in the model. The model produces zones of pyrophyllite quartz and muscovite due to the dissolution of K‐feldspar. Our model simulates, in a simplified way, the acid‐induced alteration assemblages observed in various guises in many significant mineral deposits. The particular aluminosilicate minerals produced in these experiments are associated with the gold deposits of the Witwatersrand Basin.     

Strain-dependent Fluid Flow Defined Through Rock Mass Classification Schemes

Summary ?Strain-dependent hydraulic conductivities are uniquely defined by an environmental factor, representing applied normal and shear strains, combined with intrinsic material parameters representing mass and component deformation moduli, initial conductivities, and mass structure. The components representing mass moduli and structure are defined in terms of RQD (rock quality designation) and RMR (rock mass rating) to represent the response of a whole spectrum of rock masses, varying from highly fractured (crushed) rock to intact rock. These two empirical parameters determine the hydraulic response of a fractured medium to the induced-deformations. The constitutive relations are verified against available published data and applied to study one-dimensional, strain-dependent fluid flow. Analytical results indicate that both normal and shear strains exert a significant influence on the processes of fluid flow and that the magnitude of this influence is regulated by the values of RQD and RMR.     

Shear Band Formation in Numerical Simulations Applying a Continuum Damage Rheology Model

In seismically active regions, faults nucleate, propagate, and form networks that evolve over time. To simulate crustal faulting processes, including the evolution of fault-zone properties, a rheological model must incorporate concepts such as damage rheology that describe the various stages of the seismic cycle (strain localization, subcritical crack growth and macroscopic failure) while accounting for material degradation and healing and off-fault deformation. Here we study the fundamental patterns of strain-localisation within the framework of a continuum damage rheology by performing a shear band analysis (linear instability analysis) and comparing predictions of shear band orientations with numerical results of the nonlinear problem. We find (analytically and numerically) that the angle between the shear band and the less compressive (transverse) direction is between 47° in compression tests with a strain ratio of 0.25 (highly confined compression test), and 60° for a strain ratio higher than 1.4 (axial compression and transverse extension). In addition we find that shear bands exhibit local dilation (I ₁ > 0) in a wide range of strain ratios excluding only simulations with highly confined compression (which yield compacting shear bands or non-localized deformation). Finally, we discuss the applicability of the damage model for simulating deformation in the seismogenic, brittle crust.     

Towards realistic simulations of lava dome growth using the level set method

The level set method has been implemented in a computational volcanology context. New techniques are presented to solve the advection equation and the reinitialisation equation. These techniques are based upon an algorithm developed in the finite difference context, but are modified to take advantage of the robustness of the finite element method. The resulting algorithm is tested on a well documented Rayleigh–Taylor instability benchmark [19], and on an axisymmetric problem where the analytical solution is known. Finally, the algorithm is applied to a basic study of lava dome growth.     

PRACTICAL AND PREDICTIVE MODELLING OF ORE DEPOSITS IN HYDROTHERMAL SYSTEMS

Over the past five years,we have been making efforts to develop a practical and predictive tool to explore for giant ore deposits in hydrothermal systems.Towards this goal,a significant progress has been made towards a better understanding of the basic physical and chemical processes behind ore body formation and mineralization in hydrothermal systems.