Innovative imaging of iron deposits using cross‐gradient joint inversion of potential field data with petrophysical correlation

This study demonstrates the application of the cross‐gradient joint inversion method to investigate iron mineralization zones within a volcano‐sedimentary environment. The presence of minerals with intense contrasts in density or magnetic susceptibility, such as hematite or magnetite, facilitates modelling the distribution of ore bodies with depth. Our approach involves establishing a unified interpretation of reconstructed density and susceptibility models through both independent and joint inversion with sparsity regularization in conjunction with a petrophysical model resulting from core data. This approach provides an ideal strategy to uncover the realistic geologic setting of iron ore deposits. We initially simulated a synthetic model closely resembling real‐case scenarios to assess the efficacy of the cross‐gradient joint inversion algorithm in comparison to independent inversion. Subsequently, the inversion algorithms were implemented on gravity and magnetic data, collected over an area of 500 × 600 m2 in Shavaz iron‐bearing deposits located in the central Iranian block. The primary iron oxide–apatite type mineralization in the study area is associated with the Nain–Dehshir–Baft fault as a NW–SE trending strike‐slip fault. Although both inversion methods yield satisfactory models, incorporating the cross‐gradient constraint in joint inversion resulted in a more constrained delineation of iron–oxide ore deposits in the fault system. This improvement facilitates the differentiation between hematite and a small percentage of magnetite, providing a more accurate estimation of ore depth. Inversion results suggest that the magnetite mineralization is coated with extensive hematite mineralization and both are positioned relatively within the same depth interval, covered by approximately a 15–25 m sequence of sediments.