Potential-field geophysical data inversion for 3D modelling and reserve estimation (Example of the Hajjar mine, Guemassa massif, Morocco): magnetic and gravity data case

Geophysical data inversion is a tool, which can be used to recover the subsurface distribution of physical properties from field data. Each type of geophysical data can be inverted using one or more inversion algorithms. In this paper, a set of geophysical magnetic and gravity data of the Hajjar area in Morocco, covering an extent of $3.2\times 1.6~\mathrm{km}^{2}$, were used to make a 3D model of an orebody and to estimate the mineral reserve by potential-field geophysical data inversion and excess mass estimation. We thus promote the development and application of potential-field geophysical data inversion using the softwares Geosoft Oasis Montaj and Voxi Earth Modelling™ and the evaluation of its power compared to the excess mass estimation method. The process of inversion begins with data processing, then moves to analysis and interpretation, and ends with unconstrained Cartesian cut cell inversion. The results show a variation of $-$0.22 mGal to 1.59 mGal for the gravity residual anomaly map, leading to have density variations from $2.45~\mathrm{g}/\mathrm{cm}^{3}$ to $4.22~\mathrm{g}/\mathrm{cm}^{3}$, and a variation of $-$232 nT to 1018 nT for the reduced magnetic anomaly map.Moreover, data inversion allowed us to create a 3D model of the orebody and of the adjacent geological formation, and to estimate the different parameters that characterize the orebody derived from the inversion results, which have been confirmed from survey data: (depth $\approx 160$ m; maximum density $\approx 4.22~\mathrm{g}/\mathrm{cm}^{3}$; minimum density $\approx 3~\mathrm{g}/\mathrm{cm}^{3}$; mean density $\approx 3.61~\mathrm{g}/\mathrm{cm}^{3}$; thickness of the overburden $\approx$120 m; dip $\approx 45^{\circ }$; morphology ${\approx }$ lens; volume $\approx 4.8 \times 10^{6}~\mathrm{m}^{3}$).It was therefore possible to evaluate the reserve, and to validate the reliability of the inversion by having a root mean square error between the exploited reserve and the calculated reserve of 13.5%, i.e. an insignificant difference between the real and calculated magnetic and gravity orebody responses, which support the validity of the results

Potential-field geophysical data inversion for 3D modelling and reserve estimation (Example of the Hajjar mine, Guemassa massif, Morocco): magnetic and gravity data case

Geophysical data inversion is a tool, which can be used to recover the subsurface distribution of physical properties from field data. Each type of geophysical data can be inverted using one or more inversion algorithms. In this paper, a set of geophysical magnetic and gravity data of the Hajjar area in Morocco, covering an extent of $3.2\times 1.6~\mathrm{km}^{2}$, were used to make a 3D model of an orebody and to estimate the mineral reserve by potential-field geophysical data inversion and excess mass estimation. We thus promote the development and application of potential-field geophysical data inversion using the softwares Geosoft Oasis Montaj and Voxi Earth Modelling™ and the evaluation of its power compared to the excess mass estimation method. The process of inversion begins with data processing, then moves to analysis and interpretation, and ends with unconstrained Cartesian cut cell inversion. The results show a variation of $-$0.22 mGal to 1.59 mGal for the gravity residual anomaly map, leading to have density variations from $2.45~\mathrm{g}/\mathrm{cm}^{3}$ to $4.22~\mathrm{g}/\mathrm{cm}^{3}$, and a variation of $-$232 nT to 1018 nT for the reduced magnetic anomaly map.Moreover, data inversion allowed us to create a 3D model of the orebody and of the adjacent geological formation, and to estimate the different parameters that characterize the orebody derived from the inversion results, which have been confirmed from survey data: (depth $\approx 160$ m; maximum density $\approx 4.22~\mathrm{g}/\mathrm{cm}^{3}$; minimum density $\approx 3~\mathrm{g}/\mathrm{cm}^{3}$; mean density $\approx 3.61~\mathrm{g}/\mathrm{cm}^{3}$; thickness of the overburden $\approx$120 m; dip $\approx 45^{\circ }$; morphology ${\approx }$ lens; volume $\approx 4.8 \times 10^{6}~\mathrm{m}^{3}$).It was therefore possible to evaluate the reserve, and to validate the reliability of the inversion by having a root mean square error between the exploited reserve and the calculated reserve of 13.5%, i.e. an insignificant difference between the real and calculated magnetic and gravity orebody responses, which support the validity of the results

Airborne data processing and 3D modeling of magnetic structures between Nador region and the Chaffarine Islands: Geological significances

The Cap of Trois Fourches-Nador-Chaffarine Islands zone (CTF-N-CI) located in the Moroccan Oriental Rif is recognized by Neogene-Quaternary volcanic outcrops. Those later formed volcanic cones and are considered the main source of the huge Nador magnetic anomaly. This study aims to enhance the magnetic context using geophysical techniques and identify the geometry of magnetic sources by performing 3D modeling using Talwani method. We have achieved the following findings: (1) The elaborated residual magnetic map (Rs) that shows three aligned magnetic anomalies. They are oriented W-E and register a maximum intensity of 380 nT. Applying reduction to the pole to (RTP) the Rs map makes the anomalies on the top of their body sources. Thus, we have found that those anomalies have resulted essentially from the Gourougou volcano and Chaffarine Islands volcanoes outcropping in the study area. Except for the middle Pm anomaly, which is related to volcanic rocks of the Kariate Arekmane drillhole. (2) Performing horizontal derivative (HD) technique leads to establishing the edges of the source bodies. Hence, the inversed map of the magnetization distribution proves that the three bodies are homogeneous and probably splintered from a unique mass. (3) The downward up continuation (DC) to different heights (from 500 to 2000 m) indicates the source's depth. Thus, we have demonstrated that the common source could reach 0.6 km of altitude. (4) Finally, we have constructed a 3D Talwani model of the volcanic edifices. The released 3D model demonstrates that the anomalies are related to a very extensive mass of 80 km in length, 20 km in width, 1.3 km in thickness. This model supports that the Neogene-Quaternary volcanism of the study area was fissural and related to the W-E fracture zone of up to 80 km in length. This volcanism was also affected by the sinistral Nekor fault, which probably made the Trois Fourches Cap body isolated from the other volcanic edifices.Previously mentioned elements probably prove the simultaneous establishment of this large magnetic structure. Moreover, they emphasize the crucial role of the magmatic phenomena which occurred in Northeastern Morocco at the beginning of the Messinian

Semi-Automatic Image Processing System of Aeromagnetic Data for Structural and Mining Investigations (Case of Bou Azzer Inlier, Central Anti-Atlas, Morocco)

Numerical analysis of geophysical data to uncover Precambrian belts and probably to enclose mineral deposits is becoming once more communal in mining activity. The method is founded on typifying zones branded to comprehend deposits and looking for analogous areas. The proposed work outlines a semi-automatic image processing system for the structural and mining investigation of the Bou Azzer inlier, which varies from preceding approaches as it is centered only on aeromagnetic data. The aeromagnetic signature of what seem to be geologically expressive features are pursued within the aeromagnetic items. Cobalt and associated mineralizations in the Bou Azzer inlier are recognized to arise nearby main crustal discontinuities revealing as significant shear zones, which turn act as drains for mineralizing fluids. Mineralization occurs in sectors of structural complexity beside the shear zones. Developing towards the semi-automatic uncovering of such regions, the furthermost prospective extents are those everywhere inferred structural complexity occurs next to the regions of magnetic discontinuity. The proposed method is mainly based on the approach developed by the center for exploration targeting. The study was led by means of aeromagnetic data from the Bou Azzer inlier, which is considered one of the most productive and prospective regions for minerals and base metal mineralization in Morocco. The combined results obtained from geological and geophysical data prove that prospective areas have a dominant trend of NNE-SSW, NW-SE, NNW-SSE, E-W, and NE-SW directions. The CET Grid and Porphyry Analyses show that the probable porphyry mineral deposit locations mainly concentrated in the center of inlier, the Foum Zguid dyke, and northern and eastern part, which correspond to the Bou Azzer ophiolitic complex and platform deposits of the Lower Neoproterozoic Tachdamt-Bleïda

Mapping of Hydrothermal Alteration Zones in the Kelâat M’Gouna Region Using Airborne Gamma-Ray Spectrometry and Remote Sensing Data: Mining Implications (Eastern Anti-Atlas, Morocco)

The mapping of hydrothermal alteration zones associated with mineralization is of paramount importance in searching for metal deposits. For this purpose, targeting alteration zones by analyzing airborne geophysical and satellite imagery provides accurate and reliable results. In the Kelâat M’Gouna inlier, located in the Saghro Massif of the Moroccan Anti Atlas, natural gamma-ray spectrometry and ASTER satellite data were used to map hydrothermal alteration zones. Natural gamma-ray spectrometry data were processed to produce maps of Potassium (K in %), Uranium (eU in ppm), Thorium (eTh in ppm) and ratios of K/eTh and K/eU. In addition, four-band ratios were computed, on ASTER data, to map the distribution of clay minerals, phyllitic minerals, propylitic minerals, and iron oxides. The combined results obtained from geophysical and satellite data were further exploited by fuzzy logic modelling using the Geographic Information System (GIS) to generate a mineral prospectivity map. Seven hydrothermal alteration zones likely to be favorable for mineralization have been identified. They show a spatial correlation with (i) known surface prospects and mineral occurrences, (ii) the granite-encasing contact zone, and (iii) the fault zones (Sidi Flah and Tagmout faults). This research therefore provides important information on the prospecting of mineral potential in the study area