Paleo-stratigraphic permeability anisotropy controls supergene mimetic martite goethite deposits

The Hamersley Basin in Western Australia is one of the world's largest iron ore-producing regions, hosting two types of ore in banded iron formations: the high-grade martite-microplaty haematite and the supergene martite-goethite ores. With the high-grade ores almost entirely mined in the last decade, the supergene ores have more recently become the dominant resource of interest. Consequently, understanding the genesis of these martite-goethite deposits is a critical step for exploration. Yet, although various models exist, there is still no consensus on how these mineral resources formed, complicating the prediction of resource volume and location. Here, we show that the paleo-stratigraphic permeability anisotropy (with higher permeability along strata than across) controls the supergene mimetic enrichment transport process and, subsequently, the mineralisation distribution. We introduce a flow model that implicitly represents strata with a potential function that orients the permeability tensor accurately. The numerical solver uses automatic mesh adaptivity to deliver robust solutions. By accurately reproducing the mineralisation patterns in specific deposits, we identify and quantify the paleo-water table level and permeability anisotropy ratio as the two main controlling parameters for the mineralisation distribution. These insights provide new timing constraints for the mineralisation and the physical process of iron enrichment, suggesting much more potential mineralisation volume in the paleo-reconstructed zones than previously anticipated. These flow models allow us to draw geological conclusions with few a priori assumptions required for the genetic model in which the transport component is dominant. The predictive power of this methodology will allow targeted drilling to narrow down the prospective areas and lower exploration costs. Furthermore, the methodology's generality applies to other commodities in sedimentary basins involving supergene processes and will improve our understanding of various genetic models.     

The Continuum Voronoi Block Model for simulation of fracture process in hard rocks

In this paper, a numerical model is presented to represent the fracture process in hard rocks based on a pseudo-discontinuum method called the Continuum Voronoi Block Model (CVBM). To validate this tool, numerical models for one Brazilian test, one unconfined compression test, and multiple triaxial compression tests with different confining stress were calibrated to match laboratory test results for Creighton granite. The model proved robust and matched the following macro-properties: crack initiation (CI) stress, (CD) stress, peak strength, tensile strength, Young's Modulus, and Poisson's ratio. The calibrated model served as a basis for a sensitivity study to analyze how micro-properties influence the rock's macroscopic responses. From the sensitivity study, a calibration methodology was proposed, which shall facilitate the use of the CVBM in future works.