Search

Hydrothermal ore‐forming systems usually involve coupled processes among pore‐fluid flow, heat transfer, mass transport and chemical reactions in fluid‐saturated porous media. The finite element method (FEM) has provided a powerful tool to solve hydrothermal ore‐forming problems associated with pore‐fluid convective flow. Depending on different interests of research, a hydrothermal ore‐forming system may be simulated at two different length‐scale models, such as a regional (i.e., large length‐scale) model or a deposit (i.e., small length‐scale) model. However, there is an unsolved issue to guarantee the consistency between the applied boundary conditions to a deposit model and the predicted results, at exactly the same location as the boundary of the deposit model, from using the regional model. This paper proposes a FEM‐based dual length‐scale simulation procedure to address this issue. In the proposed procedure, a regional model of coarse mesh is first simulated by the FEM. Since a deposit model is located within the regional model, the applied boundary conditions to the deposit model of fine mesh can be determined from the simulation results of the regional model, so that the consistency can be guaranteed at the common boundary between the deposit model and the regional model. Main outcomes of this study have demonstrated that: (1) the proposed FEM‐based dual length‐scale simulation procedure is suitable and applicable for solving coupled problems involving pore‐fluid flow, heat transfer, mass transport and chemical reactions in fluid‐saturated porous media; (2) the proposed procedure is predictive because it is unnecessary to know the location of an ore deposit prior to the simulation of a hydrothermal ore‐forming system; (3) the boundary condition of a deposit model of fine mesh is consistent with the simulation results obtained from the regional model of coarse mesh, so that the reliability of the simulation results obtained from the deposit model can be ensured. [ABSTRACT FROM AUTHOR]

Skarn-hosted deposits are commonly recognized as the consequence of magma intrusion within the Earth's upper crust. The Ergu Pb-Zn deposit can be regarded as a typical skarn-hosted deposit in the hydrothermal ore-forming system within the central Lesser Xing'an Range (LXR), Heilongjiang, China. Although extensive studies were conducted to understand the ore-forming system associated with the Ergu Pb-Zn deposit through using the traditional geoscience methods, the ore-forming process involved in this deposit has not been justified in a strictly scientific manner to date. In this paper, the hydrothermal ore-forming process of the Ergu Pb-Zn deposit is computationally simulated through using the dual length-scale approach associated with the finite element method (FEM). The related computational simulation results have demonstrated that: (1) the pore-fluid convection provides continuous ore-forming fluid and material sources for the Ergu Pb-Zn deposit at the quartz-Pb-Zn sulfide stage; (2) the convective flow of the pore fluid is the main dynamic mechanism, which controls the temperature, chemical species and pore-fluid velocity distributions in the Ergu Pb-Zn deposit; (3) the localized structure plays a key role in controlling the localized pore-fluid flow pattern, which can further control the location and formation of the orebody grade in the Ergu Pb-Zn deposit; (4) the dual length-scale approach associated with the FEM is very useful in dealing with the computational simulation of the hydrothermal ore-forming mechanism involved in the Ergu Pb-Zn deposit. [ABSTRACT FROM AUTHOR]

[Display omitted] • To computationally simulate the ore-forming process in the Luziyuan deposit. • To determine the boundary conditions of the deposit model in a scientific way. • To identify the controlling ore-forming mechanism in the Luziyuan deposit. • To use computational simulations for identifying the location of an ore deposit. The giant Luziyuan Pb-Zn deposit is a typical hydrothermal ore-forming deposit in the southern Baoshan block of SE Tibet. Although extensive studies were conducted to understand the ore-forming system associated with the Luziyuan Pb-Zn deposit through using the traditional geoscience methods, the ore-forming process involved in this deposit has not been justified in a strictly scientific manner to date. In this paper, the hydrothermal ore-forming process of the Luziyuan Pb-Zn deposit is computationally simulated through using the dual length-scale approach associated with the finite element method (FEM). In the dual length-scale approach, the FEM is first used to simulate the regional model of a large length-scale in the geological region where the Luziyuan Pb-Zn deposit is located. Afterward, the FEM is then used to simulate the deposit model of the Luziyuan Pb-Zn deposit. Particularly, the boundary conditions of the deposit model are consistently determined from using the computational simulation results of the regional model. The related numerical results demonstrated that: (1) the pore-fluid convection provides a continuous source of ore-forming fluid and ore-forming material for the Luziyuan Pb-Zn deposit; (2) the convective pore-fluid flow is the main dynamic mechanism, which controls the temperature, pore-fluid velocity and chemical species distribution in the Luziyuan Pb-Zn deposit; (3) the localized structure is the main controlling factor for the localized pore-fluid flow pattern, which further controls the formation and location of the orebody in the Luziyuan Pb-Zn deposit; (4) the dual length-scale approach associated with the FEM is very useful for simulating the hydrothermal ore-forming process of the Luziyuan Pb-Zn deposit. [ABSTRACT FROM AUTHOR]

Abstract Convective seepage flow plays a key role not only in naturally forming mineral deposits and oil reservoirs, but also in the carbon-dioxide sequestration in fluid-saturated rocks. This paper presents a steady-state numerical solver, which is based on the finite element method (FEM) and progressive asymptotic approach procedure (PAAP), to solve steady-state convective seepage flow problems in fluid-saturated heterogeneous rocks. Although this kind of flow problem is commonly solved by using transient-state numerical solvers, the use of the proposed steady-state numerical solver, which is the main characteristic of this study, can have certain advantages. After the proposed steady-state numerical solver is verified by a benchmark testing problem, for which analytical solutions are available for comparison, it is further applied to simulate steady-state convective seepage flow in the Australia North West Shelf basin, which is composed of naturally-formed heterogeneous porous rocks. Through the simulated numerical results, it can be found that: (1) compared with the analytical solutions of the benchmark testing problem, the use of the steady-state numerical solver can produce correct and accurate numerical simulation results for dealing with convective seepage flow problems in fluid-saturated heterogeneous rocks. (2) The main advantage of using the steady-state numerical solver is to avoid the need of considering initial conditions of the problem, which are commonly difficult to be determined because the convective seepage flow within the upper crust is a past event in geology. (3) In the Australia North West Shelf basin, the convective seepage flow and temperature distribution patterns depend strongly on both the permeability contrast of the faults and the basin geometry (including the layer structures and fault locations) in the computational model. Highlights • To present and verify the numerical method for modeling large-scale convective seepage flow in heterogeneous rocks. • To simulate large-scale convective seepage flow in the Australia North West Shelf basin. • Rock permeability heterogeneity can affect large-scale convective seepage flow and temperature patterns. • Basin geometry can also affect large-scale convective seepage flow and temperature patterns. [ABSTRACT FROM AUTHOR]

Mountain topography can cause both pressure-gradient driven advective flow and temperature-gradient driven convective flow. Compared with the pressure-gradient driven advective flow, the temperature-gradient driven convective flow can last long enough to enable the hydrothermal ore-forming system to reach a steady state. This paper is focused on the numerical modeling of a hydrothermal ore-forming system, in which the temperature-gradient driven convective flow is the main driving force of orebody formation and mineralization. A generic model, which can be viewed as a generalized geological model, is used in the numerical modeling. To simulate the chemical aspects of the hydrothermal ore forming system, a simplified chemical process associated with the Broken Hill type Pb-Zn deposit is considered in this study. The related numerical results demonstrated that: (1) The mountain topography plays an important role in controlling Pb and Zn mineralization distribution patterns in a hydrothermal ore-forming system; (2) The mixing of the reactive sulfide and sulfate fluids can be enhanced by the temperature-gradient driven convective flow; (3) If the duration of mineralization time is known, then the computed mineral precipitation rate can be used to predict the average (weight) grade of a mineral in the hydrothermal ore-forming system. Similarly, if the average (weight) grade of a mineral is known from a mining site, then the computed mineral precipitation rate can be used to estimate the duration of mineralization time in the hydrothermal ore-forming system. [ABSTRACT FROM AUTHOR]

Abstract: For a three-dimensional vertically-oriented fault zone, we consider the coupled effects of fluid flow, heat transfer and reactive mass transport, to investigate the patterns of fluid flow, temperature distribution, mineral alteration and chemically induced porosity changes. We show, analytically and numerically, that finger-like convection patterns can arise in a vertically-oriented fault zone. The onset and patterns of convective fluid flow are controlled by the Rayleigh number which is a function of the thermal properties of the fluid and the rock, the vertical temperature gradient, and the height and the permeability of the fault zone. Vigorous fluid flow causes low temperature gradients over a large region of the fault zone. In such a case, flow across lithological interfaces becomes the most important mechanism for the formation of sharp chemical reaction fronts. The degree of rock buffering, the extent and intensity of alteration, the alteration mineralogy and in some cases the formation of ore deposits are controlled by the magnitude of the flow velocity across these compositional interfaces in the rock. This indicates that alteration patterns along compositional boundaries in the rock may provide some insights into the convection pattern. The advective mass and heat exchanges between the fault zone and the wallrock depend on the permeability contrast between the fault zone and the wallrock. A high permeability contrast promotes focussed convective flow within the fault zone and diffusive exchange of heat and chemical reactants between the fault zone and the wallrock. However, a more gradual permeability change may lead to a regional-scale convective flow system where the flow pattern in the fault affects large-scale fluid flow, mass transport and chemical alteration in the wallrocks. [Copyright &y& Elsevier]

Abstract: Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. [Copyright &y& Elsevier]

We use the finite element method to model three-dimensional convective pore-fluid flow in fluid-saturated porous media when they are heated from below. In particular, we employ the particle-tracking technique to mimic the trajectories of particles in three-dimensional fluid flow problems. The related numerical results demonstrated that: (1) The progressive asymptotic approach procedure, which was previously developed for the finite element modelling of two-dimensional convective pore-fluid flow problems, is equally applicable to the finite element modelling of three-dimensional convective pore-fluid flow in fluid-saturated porous media heated from below. (2) The perturbation of gravity at different planes has a significant effect on the pattern of three-dimensional convective pore-fluid flow and therefore, may influence the pattern of orebody formation and mineralization in three-dimensional hydrothermal systems. Copyright © 2001 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]