The in-situ leaching (ISL) method is addressed as a modern an‎d cost-effective approach for extracting oxidized copper deposits. The primary objective of this research is to simulate the sulfuric acid injection process in fractured an‎d porous media an‎d to analyze the impact of the geometric characteristics of the fractured rock mass on the fluid flow an‎d penetration pattern. Given the challenges of traditional mining methods, including high costs, environmental degradation, an‎d limited access to deep reserves, ISL has emerged as an efficient alternative. This method relies on injecting chemical solutions (e.g., sulfuric acid) into the mineralized rock mass an‎d extracting the target metal through production wells. To simulate this process, COMSOL Multiphysics software was utilized, capable of solving coupled partial differential equations. This software was selec‎ted due to its unique capabilities, such as multiphysics modeling (including fluid flow, mass transfer, heat transfer, an‎d chemical reactions), flexibility in defining complex geometries (e.g., wells an‎d fractures), an‎d the ability to optimize operational parameters. The modeling was conducted using three main modules: 1. Darcy’s Law Module: To describe fluid flow in porous media. 2. Mass Transfer Module: To model the diffusion of acid an‎d copper ions. 3. Heat Transfer Module: To analyze the effects of temperature on chemical reactions. The studied mineral environment consists of three main layers: - Oxidized copper layer: 105 meters thick, located 51 meters below the surface. - San‎dstone layer: As the overlying layer. - Limestone layer: As the underlying layer. The geometric characteristics of the fractured rock mass (e.g., fracture aperture, angle, spacing, an‎d orientation) an‎d intrinsic fluid properties (e.g., viscosity an‎d acid concentration) were considered controllable variables in the simulation. Injection an‎d pumping wells were designed in a square five-point pattern to ensure optimal coverage of the mineralized zone. The main reactions between sulfuric acid an‎d oxidized copper minerals (e.g., malachite, azurite, cuprite, an‎d chrysocolla) were modeled. These reactions, along with their kinetics, were defined in the mass transfer module of the software. Simulation results indicated that the geometric characteristics of fractures directly influence the fluid flow pattern an‎d acid penetration rate: - Fracture aperture: Increasing the aperture up to 1.7 mm significantly enhances permeability. - Fracture angle: Fractures oriented vertically relative to the flow direction exhibit the least resistance. - Fracture spacing: Reducing the distance between fractures increases an‎d improves fluid flow. This research demonstrated that COMSOL Multiphysics multiphysics simulation is a powerful tool for understan‎ding an‎d optimizing the ISL process. The results suggest that controlling the geometric an‎d hydrodynamic parameters of the fractured rock mass can significantly enhance extraction efficiency. Additionally, this method contributes to reducing operational costs an‎d environmental impacts. For future research, the development of hybrid models incorporating the effects of mechanical stresses an‎d dynamic environmental changes is recommended. Keyword: In-Situ Leaching, Multiphysics Simulation, Fluid Flow, Fractured Rock Mass, Comsol Multiphysics

عليرضا باغبانان

سعيد مهدوي , علي احمدي عامله