Exploration
parnian javadisharif; Alireza Arab Amiri; Behzad Tokhmechi (غیرفعال); Fereydoun Sharifi
Abstract
The technique referred to as Complex Resistivity (CR) or Spectral Induced Polarization (SIP) possesses the capability to distinguish between various kinds of minerals or the sources of induced polarization by utilizing the physical characteristics of minerals or polarizable inclusions. The Generalized ...
Read More
The technique referred to as Complex Resistivity (CR) or Spectral Induced Polarization (SIP) possesses the capability to distinguish between various kinds of minerals or the sources of induced polarization by utilizing the physical characteristics of minerals or polarizable inclusions. The Generalized Effective Medium Theory of Induced Polarization (GEMTip) model is utilized to derive physical characteristics from SIP data. Different inversion methods are applied for this task, though they encounter difficulties such as computational costs, non-linearity, and the intricacy of the inverse issue. To tackle this, a new inversion approach based on Deep Learning (DL) via Convolutional Neural Network (CNN) is proposed for predicting the parameters of polarizable particles from SIP data. The CNN is trained on 20000 synthetic datasets produced using the GEMTip forward model. While DL networks address non-linearities, specific modifications are applied to synthetic datasets to evaluate the influence of non-linearity and correlation on the results. In the Kervian region southwest of Saqqez city, gold mineralization is linked to quartz and pyrite minerals, with two types of pyrite recognized - coarse-grained barren and fine-grained auriferous. The existence of sulfide mineral pyrite, along with variations in pyrite sizes, presents an attractive target for SIP exploration in the investigated area. The trained network is also validated on Gravian data and effectively retrieves parameters as evidenced by the data. The proposed methodology simplifies the inversion process by estimating parameters in one step, enabling a direct and efficient procedure.
Saeed Nazari; Alireza Arab Amiri; Abolghasem Kamkar Rouhani; Fereydoun Sharifi
Abstract
In this work, we simulate the frequency-domain helicopter-borne electromagnetic (HEM) data over the two-dimensional (2D) and three-dimensional (3D) earth models. In order to achieve this aim, the vector Helmholtz equation is used to avoid the convergence problems in Maxwell’s equations, and the ...
Read More
In this work, we simulate the frequency-domain helicopter-borne electromagnetic (HEM) data over the two-dimensional (2D) and three-dimensional (3D) earth models. In order to achieve this aim, the vector Helmholtz equation is used to avoid the convergence problems in Maxwell’s equations, and the corresponding fields are divided into primary and secondary components. We use the finite difference method on a staggered grid to discretize the equations, which can be performed in two ways including the conventional and improved finite difference methods. The former is very complex in terms of programming, which causes errors. Furthermore, it requires different programming loops over each point of the grid, which increases the program’s running time. The latter is the improved finite difference method (IFDM), in which pre-made derivative matrices can be used. These pre-made derivative matrices can be incorporated into the derivative equations and convert them directly from the derivative form to the matrix form. After having the matrix form system of linear equations, Ax = b is solved by the quasi-minimal residual (QMR). IFDM does not have the complexities of the conventional method, and requires much less execution time to form a stiffness or coefficient matrix. Moreover, its programing process is simple. Our code uses parallel computing, which gives us the ability to calculate the fields for all transmitter positions at the same time, and because we use sparse matrices thorough the code memory space, requires to store the files is less than 100 MB compared with normal matrices that require more than 15 GB space in the same grid size. We implement IFDM to simulate the earth’s responses. In order to validate, we compare our results with various models including the 3D and 2D models, and anisotropic conductivity. The results show a good fit in comparison with the FDM solution of Newman and the appropriate fit integral equations solution of Avdeev that is because of the different solution methods.
F. Jamali; A.R. Arab Amiri; A. Kamkar Rouhani; A. Bahrami
Abstract
In any geophysical exploration, the final goal is to achieve an accurate image of the relevant underground property. In order to achieve this, the geophysical operation using the electrical resistivity and induced polarization (IP) methods is conducted to explore the sub-surface sulfide mineralization. ...
Read More
In any geophysical exploration, the final goal is to achieve an accurate image of the relevant underground property. In order to achieve this, the geophysical operation using the electrical resistivity and induced polarization (IP) methods is conducted to explore the sub-surface sulfide mineralization. Considering the mineralization evidence in the Kaboudan area near the Bardeskan city, first, geophysical surveying of the polymetallic deposit is carried out using the electrical resistivity and IP methods by employing the rectangle array in order to detect the electrical anomalies in the area. Then for delineation of the identified anomalies and investigation of the mineralization in the area, the 2D resistivity and chargeability cross-sections are prepared and interpreted with the help of the geological information. This geophysical survey in the area has led to the identification of several potential areas for mineralization. Then in order to obtain a detailed picture of the sub-surface mineralization and an overview of the in-depth mineralization distribution, a 3D modeling of the acquired data is made, and the results of this modeling are shown in 3D forms. The mineralization zones are identified in the studied area from their high chargeability values as well as the low to medium electrical resistivity amounts. This can be attributed to the metal mineralization and the presence of sulfide minerals in the mineralization zones. Mineralization in many places of the studied area is determined with an approximate east-west trend as well as somewhat varying the intensities of the electrical resistivity and chargeability amounts. The geological and drilling information obtained from the area confirm the interpretations.
S. Najafi Ghoshebolagh; A. Kamkar Rouhani; A.R. Arab Amiri; H. Bizhani
Abstract
As many gold deposits are associated with sulfide zones, and the direct exploration of gold deposits using the geophysical methods is very difficult due to its low amount in the sub-surface, the direct exploration of sulfide zones by the geophysical electrical resistivity and induced polarization (IP) ...
Read More
As many gold deposits are associated with sulfide zones, and the direct exploration of gold deposits using the geophysical methods is very difficult due to its low amount in the sub-surface, the direct exploration of sulfide zones by the geophysical electrical resistivity and induced polarization (IP) methods may lead to the indirect exploration of gold deposits. The gold deposit in the Kervian area is located in the Kurdistan shear zone, and is directly related to the sulfide, silica, and carbonate alteration units. After acquiring the resistivity and IP data, 2D modeling of the data is made in order to indirectly identify the gold-bearing zones in the surveyed area. As some of the identified geophysical anomalies indicating the sulfide zones may not be associated with the economic amounts of gold, in order to obtain an exploration pattern for the gold deposit in the studied area, a combination of the geophysical data modeling and interpretation results with the geological information and other exploratory data is used to reduce the uncertainty in identifying the gold-bearing zones in the studied area. Thus, modeling and interpretation of the geophysical data lead to identify the sub-surface anomalies as the locations of possible gold mineralization in the area, and then the drilling points are suggested in the area. Considering the geological studies and chemical analysis of the samples taken from the drilled boreholes crossing some of the geophysical anomalies, we conclude that the geophysical anomalies occurring inside the phyllite and carbonate units in the area can contain an economic amount of gold, and thus are recommended as the top priority for further exploration.
F. Sharifi; A.R. Arab Amiri; A. Kamkar Rouhani
Abstract
The generalized effective-medium theory of induced polarization (GEMTIP) is a newly developed relaxation model that incorporates the petro-physical and structural characteristics of polarizable rocks in the grain/porous scale to model their complex resistivity/conductivity spectra. The inversion of the ...
Read More
The generalized effective-medium theory of induced polarization (GEMTIP) is a newly developed relaxation model that incorporates the petro-physical and structural characteristics of polarizable rocks in the grain/porous scale to model their complex resistivity/conductivity spectra. The inversion of the GEMTIP relaxation model parameter from spectral-induced polarization data is a challenging issue because of the highly non-linear dependency of the observed data on the model parameter and non-uniqueness of the problem. To solve these problems as well as scape the local minima of the highly complicated cost function, the genetic algorithm (GA) can be applied but it has proven to be time-intensive computationally. However, this drawback can be resolved by incorporating a faster algorithm, e.g. particle swarm optimization (PSO). The aim of this work is to investigate whether recovering the model parameter of the ellipsoidal GEMTIP model from SIP data using the combined GA and PSO algorithms is possible. To achieve this aim, we set the best calculated individuals using GA as the search space of PSO, and then the best location achieved by PSO in each iteration is assigned as the updated model parameters. The results of our research work reveal that the model parameters can effectively be recovered using the approach proposed in this paper but the time constant of a noisy data that arises from the adverse dependency of this parameter on the ellipticity of a polarizable grain. Moreover, the execution time of the ellipsoidal GEMTIP modeling of complex resistivity data can be significantly improved using the proposed algorithm.
Exploitation
M. Bavand Savadkoohi; B. Tokhmechi; E. Gloaguen; A.R. Arab-Amiri
Abstract
Computer graphics offer various gadgets to enhance the reconstruction of high-order statistics that are not correctly addressed by the two-point statistics approaches. Almost all the newly developed multiple-point geostatistics (MPS) algorithms, to some extent, adapt these techniques to increase the ...
Read More
Computer graphics offer various gadgets to enhance the reconstruction of high-order statistics that are not correctly addressed by the two-point statistics approaches. Almost all the newly developed multiple-point geostatistics (MPS) algorithms, to some extent, adapt these techniques to increase the simulation accuracy and efficiency. In this work, a scrutiny comparison between our recently developed MPS algorithm, the cross-correlation-wavelet simulation (CCWSIM), and a well-known MPS algorithm, FILTERSIM, is performed. The main motivation to benchmark these two algorithms is that both exploit some digital image processing filters for feature extraction. Indeed, both algorithms compute the similarity (or dissimilarity) between data events in simulation grid and training image in the feature space. In order to compare the accuracy of the algorithms, some statistics such as facies proportion, variogram, and connectivity function are computed. The results obtained reveal an excellent agreement of the CCWSIM realizations with the training image rather than FILTERSIM. Furthermore, on average, the required simulation runtime for CCWSIM is at least 10 times less than that for FILTERSIM.
Alireza Arab-Amiri; Fereydoun Sharifi; Abolghasem Kamkar-Rouhani
Abstract
The need for clean groundwater resources to have sustainable development in a country is undoubted. Due to the importance and high quality of karstic waters in supplying water in Iran especially in Shahrood city, it is attempted in this research work to recognize and explore karstic waters in southwest ...
Read More
The need for clean groundwater resources to have sustainable development in a country is undoubted. Due to the importance and high quality of karstic waters in supplying water in Iran especially in Shahrood city, it is attempted in this research work to recognize and explore karstic waters in southwest of Tepal area, Shahrood. For this purpose, integration of the results obtained from the methods of vertical electrical sounding (VES) and resistivity profiling has been used in this research work. The VES surveys have been performed in 10 sounding points using the Schlumberger array with electrode separations of a maximum 500 meters. The resistivity profiling surveys using dipole-dipole electrode array with 75m electrode spacing and dipole steps 1 to 8 have been carried out along four lines having a length of more than four kilometers in the study area. Then, one-dimensional (1-D) modeling and interpretation of the sounding results using master curves and IX1D software, and two-dimensional (2-D) modeling and interpretation of the profiling results using Res2DINV have been made. As a result of the interpretation and integration of the results, karstic water zones in the study area have been recognized, and based on that, suitable locations for drilling to access and extract karstic groundwater have been introduced.
