3D electromagnetic modeling of graphitic faults in the Athabasca Basin using a finite-volume time-domain approach with unstructured grids

Geophysics ◽  
2021 ◽  
pp. 1-73
Author(s):  
Xushan Lu ◽  
Colin Farquharson ◽  
Jean-Marc Miehé ◽  
Grant Harrison

Uranium exploration in the Athabasca Basin, Canada, relies heavily on ground-based transient electromagnetic (TEM) surveys to target thin, steeply dipping graphitic conductors that are often closely related to the uranium ore deposits. The interpretation of TEM data is important in identifying the locations and trends of conductors in order to guide subsequent drilling campaigns. We present a trial-and-error modeling approach and its application to the interpretation of a data set acquired at Close Lake in the Athabasca Basin. The modeling process has two key tasks: building geo-electric models and computing their TEM responses. The modeling process is repeated with the geo-electric model being iteratively refined based on the match between three-component calculated and measured data from early to late times. To create geo-electric models, we first build a realistic geological model and discretize it using an unstructured tetrahedral mesh, with each mesh cell populated with appropriate resistivities. To calculate the TEM responses of the geo-electric model, we use a 3D finite-volume time-domain (FVTD) algorithm. We construct our initial model based on existing geologic information and drilling data. We show that this modeling process is flexible and can easily handle thin, steeply dipping conductive graphitic fault models with variable resistivities in the fault and background, and with topography. Our interpretation of the Close Lake data matches well with the trend and location of the main conductor as revealed by drilling data, and also confirms the existence of a smaller conductor which only caused noticeable anomalous responses in early-time horizontal-component data. The smaller conductor was suggested by previous electromagnetic data but was missed in a recent interpretation based on the modeling of only late-time vertical component data with plate-based approximate modeling methods.

Geophysics ◽  
2013 ◽  
Vol 78 (4) ◽  
pp. E149-E159 ◽  
Author(s):  
Joel E. Podgorski ◽  
Esben Auken ◽  
Cyril Schamper ◽  
Anders Vest Christiansen ◽  
Thomas Kalscheuer ◽  
...  

Helicopter time-domain electromagnetic (HTEM) surveying has historically been used for mineral exploration, but over the past decade it has started to be used in environmental assessments and geologic and hydrologic mapping. Such surveying is a cost-effective means of rapidly acquiring densely spaced data over large regions. At the same time, the quality of HTEM data can suffer from various inaccuracies. We developed an effective strategy for processing and inverting a commercial HTEM data set affected by uncertainties and systematic errors. The delivered data included early time gates contaminated by transmitter currents, noise in late time gates, and amplitude shifts between adjacent flights that appeared as artificial lineations in maps of the data and horizontal slices extracted from inversion models. Multiple processing steps were required to address these issues. Contaminated early time gates and noisy late time gates were semiautomatically identified and eliminated on a record-by-record basis. Timing errors between the transmitter and receiver electronics and inaccuracies in absolute amplitudes were corrected after calibrating selected HTEM data against data simulated from accurate ground-based TEM measurements. After editing and calibration, application of a quasi-3D spatially constrained inversion scheme significantly reduced the artificial lineations. Residual lineations were effectively eliminated after incorporating the transmitter and receiver altitudes and line-to-line amplitude factors in the inversion process. The final inverted model was very different from that generated from the original data provided by the contractor. For example, the average resistivity of the thick surface layer decreased from [Formula: see text] to [Formula: see text], the depths to the layer boundaries were reduced by 15%–23%, and the artificial lineations were practically eliminated. Our processing and inversion strategy is entirely general, such that with minor system-specific modifications it could be applied to any HTEM data set, including those recorded many years ago.


2021 ◽  
Vol 34 (04) ◽  
pp. 1180-1199
Author(s):  
Hossein Shahnazari- Aval ◽  
Mirsattar Meshinchi-Asl ◽  
Mahmoud Mehramuz

In this study, we have implemented an edge-based finite element method for the numerical modeling of the transient electromagnetic method. We took the Helmholtz equation of the electric field as the governing equation for the edge-based finite element analysis. The modeling domain was discretized using linear tetrahedral mesh supported by Whitney-type vector basis functions. We inferred the equations by applying the Galerkin method. The system of equation was solved using a corrected version of the Bi Conjugate Gradient Stabilized Method (BiCGStab) algorithm to reduce the computational time. We obtained numerical solution for electric field in the Laplace domain; then the field was transformed into the time domain using the Gaver-Stehfest algorithm. Following this, the impulse response of the magnetic field was obtained through the Faraday law of electromagnetic induction as it is considerably more stable and computationally more efficient than inversion using the Fourier Transform. 3D geoelectric models were used to investigate the convergence of the edge-based finite element method with the analytic solution. The results are in good agreement with the analytical solution value for two resistivity contrasts in the 3D geoelectric brick model. We also compared the results of tetrahedral elements with the brick element in the 3D horizontal sheet and 3D conductive brick model. The results indicated that these two elements show very similar errors, but tetrahedral reflects fewer relative errors. For the low resistivity geoelectric model, numerical checks against the analytical solution, integral-equation method, and finite-difference time-domain solutions showed that the solutions would provide accurate results.


2020 ◽  
Vol 0 (0) ◽  
Author(s):  
Tawfik Yahya ◽  
Nur Azah Hamzaid ◽  
Sadeeq Ali ◽  
Farahiyah Jasni ◽  
Hanie Nadia Shasmin

AbstractA transfemoral prosthesis is required to assist amputees to perform the activity of daily living (ADL). The passive prosthesis has some drawbacks such as utilization of high metabolic energy. In contrast, the active prosthesis consumes less metabolic energy and offers better performance. However, the recent active prosthesis uses surface electromyography as its sensory system which has weak signals with microvolt-level intensity and requires a lot of computation to extract features. This paper focuses on recognizing different phases of sitting and standing of a transfemoral amputee using in-socket piezoelectric-based sensors. 15 piezoelectric film sensors were embedded in the inner socket wall adjacent to the most active regions of the agonist and antagonist knee extensor and flexor muscles, i. e. region with the highest level of muscle contractions of the quadriceps and hamstring. A male transfemoral amputee wore the instrumented socket and was instructed to perform several sitting and standing phases using an armless chair. Data was collected from the 15 embedded sensors and went through signal conditioning circuits. The overlapping analysis window technique was used to segment the data using different window lengths. Fifteen time-domain and frequency-domain features were extracted and new feature sets were obtained based on the feature performance. Eight of the common pattern recognition multiclass classifiers were evaluated and compared. Regression analysis was used to investigate the impact of the number of features and the window lengths on the classifiers’ accuracies, and Analysis of Variance (ANOVA) was used to test significant differences in the classifiers’ performances. The classification accuracy was calculated using k-fold cross-validation method, and 20% of the data set was held out for testing the optimal classifier. The results showed that the feature set (FS-5) consisting of the root mean square (RMS) and the number of peaks (NP) achieved the highest classification accuracy in five classifiers. Support vector machine (SVM) with cubic kernel proved to be the optimal classifier, and it achieved a classification accuracy of 98.33 % using the test data set. Obtaining high classification accuracy using only two time-domain features would significantly reduce the processing time of controlling a prosthesis and eliminate substantial delay. The proposed in-socket sensors used to detect sit-to-stand and stand-to-sit movements could be further integrated with an active knee joint actuation system to produce powered assistance during energy-demanding activities such as sit-to-stand and stair climbing. In future, the system could also be used to accurately predict the intended movement based on their residual limb’s muscle and mechanical behaviour as detected by the in-socket sensory system.


Geophysics ◽  
1986 ◽  
Vol 51 (7) ◽  
pp. 1462-1471 ◽  
Author(s):  
Brian R. Spies ◽  
Dwight E. Eggers

Problems and misunderstandings arise with the concept of apparent resistivity when the analogy between an apparent resistivity computed from geophysical observations and the true resistivity structure of the subsurface is drawn too tightly. Several definitions of apparent resistivity are available for use in electromagnetic methods; however, those most commonly used do not always exhibit the best behavior. Many of the features of the apparent resistivity curve which have been interpreted as physically significant with one definition disappear when alternative definitions are used. It is misleading to compare the detection or resolution capabilities of different field systems or configurations solely on the basis of the apparent resistivity curve. For the in‐loop transient electromagnetic (TEM) method, apparent resistivity computed from the magnetic field response displays much better behavior than that computed from the induced voltage response. A comparison of “exact” and “asymptotic” formulas for the TEM method reveals that automated schemes for distinguishing early‐time and late‐time branches are at best tenuous, and those schemes are doomed to failure for a certain class of resistivity structures (e.g., the loop size is large compared to the layer thickness). For the magnetotelluric (MT) method, apparent resistivity curves defined from the real part of the impedance exhibit much better behavior than curves based on the conventional definition that uses the magnitude of the impedance. Results of using this new definition have characteristics similar to apparent resistivity obtained from time‐domain processing.


2012 ◽  
Vol 183 (8) ◽  
pp. 1634-1640 ◽  
Author(s):  
Adnane Hamiaz ◽  
Rudy Klein ◽  
Xavier Ferrieres ◽  
Olivier Pascal ◽  
Jean-Pierre Boeuf ◽  
...  

2013 ◽  
Vol 57 (03) ◽  
pp. 125-140
Author(s):  
Daniel A. Liut ◽  
Kenneth M. Weems ◽  
Tin-Guen Yen

A quasi-three-dimensional hydrodynamic model is presented to simulate shallow water phenomena. The method is based on a finite-volume approach designed to solve shallow water equations in the time domain. The nonlinearities of the governing equations are considered. The methodology can be used to compute green water effects on a variety of platforms with six-degrees-of-freedom motions. Different boundary and initial conditions can be applied for multiple types of moving platforms, like a ship's deck, tanks, etc. Comparisons with experimental data are discussed. The shallow water model has been integrated with the Large Amplitude Motions Program to compute the effects of green water flow over decks within a time-domain simulation of ship motions in waves. Results associated to this implementation are presented.


2021 ◽  
Author(s):  
Temirlan Zhekenov ◽  
Artem Nechaev ◽  
Kamilla Chettykbayeva ◽  
Alexey Zinovyev ◽  
German Sardarov ◽  
...  

SUMMARY Researchers base their analysis on basic drilling parameters obtained during mud logging and demonstrate impressive results. However, due to limitations imposed by data quality often present during drilling, those solutions often tend to lose their stability and high levels of predictivity. In this work, the concept of hybrid modeling was introduced which allows to integrate the analytical correlations with algorithms of machine learning for obtaining stable solutions consistent from one data set to another.


Sign in / Sign up

Export Citation Format

Share Document