Hendra Grandis 1 and Prihadi Sumintadireja 2

Size: px
Start display at page:

Download "Hendra Grandis 1 and Prihadi Sumintadireja 2"

Transcription

1 A BRIEF REVIEW FOR THE PROPER APPLICATION OF MAGNETOTELLURIC (MT) AND CONTROLLED-SOURCE AUDIO-FREQUENCY MAGNETOTELLURIC (CSAMT) IN GEOTHERMAL EXPLORATION Hendra Grandis 1 and Prihadi Sumintadireja 2 ABSTRACT (1) Applied Geophysics Research Division Faculty of Mining and Petroleum Engineering Institut Teknologi Bandung grandis@geoph.itb.ac.id, grandis@earthling.net (2) Applied Geology Research Division Faculty of Earth Science and Technology Institut Teknologi Bandung prihadi@gc.itb.ac.id Magnetotelluric (MT) method is a frequency domain sounding technique that uses naturally existing electromagnetic (EM) fields as the source to infer the subsurface resistivity distribution. Controlled- Source Audio-frequency Magnetotellurics (CSAMT) is derived from MT and uses an artificial EM source (typically grounded electric dipole) in the range 0.1 Hz to 10 khz. Both MT and CSAMT are well recognized as geophysical methods in geothermal exploration since both results in resistivity structure that is more or less related to the temperature regime of the subsurface. However, from empirical point of view sometimes one method (e.g. CSAMT) delivers superior results than the other method (e.g. MT) in geothermal exploration, and vice versa. Successful applications of one technique were then used to penalize the other technique, especially in Indonesia. This situation is mainly due to lack of knowledge of the basic concept of both MT and CSAMT and also their processing and modeling aspects. Unfortunately, this point of view influences policy and decision making in adopting one method by disregarding the other in the exploration program. This paper describes both MT and CSAMT techniques from theoretical and practical perspectives in order to underline the possible difference of outcomes when they are applied in geothermal exploration. Synthetic data as the response of simple models as well as real (field) data are used to exemplify pitfalls that may generate misleading interpretations. We hope that in the future any exploration program employing geophysical methods, especially MT or CSAMT will be based on correct and well-founded scientific reasoning related to specific target of the survey. INTRODUCTION In geothermal exploration, electrical geophysical methods are of prime importance due to the close and direct relationships between temperature and resistivity of materials. Hydrothermal fluids and also their overlying cap rock have large resistivity contrasts relative to their host environment. Magnetotellurics (MT) and its controlled-source, high-frequency complement CSAMT are the main tools to delineate such geological features and the basic structural controls based on their resistivity expression. Case histories have shown that both MT (e.g. Newman et al., 2008; Cumming and Mackie, 2010) and CSAMT (e.g. Wannamaker, 1997) have been successful in establishing geothermal resources. There were less-successful results from MT and CSAMT as well but with very limited or even no publication at all. Successful applications of one technique were then used to penalize the other technique. Unfortunately, this point of view influences policy and decision making in adopting one method by disregarding the other in the geothermal exploration program. This situation is mainly due to lack of knowledge of the basic concept of both MT and CSAMT and also their processing and modeling aspects. This paper describes both MT and CSAMT techniques from theoretical and practical perspectives in order to underline the possible difference of outcomes when they are applied in

2 geothermal exploration. Synthetic data as the response of simple models as well as real (field) data are used to exemplify pitfalls that may generate misleading interpretations. MT AND CSAMT METHODS Both MT and CSAMT are frequency domain sounding techniques employing variations of electromagnetic (EM) fields as the primary source. Orthogonal components of electric (E x, E y ) and magnetic (H x, H y ) fields are recorded as time series data and then processed by using spectral analysis that results in complex frequencydependent impedance (Z). Impedance is the ratio of the electric and magnetic fields, from which the information on the subsurface can be inferred. For a complete set of orthogonal components, electric and magnetic fields are related by the impedance tensor (Z), E E x y Z = Z xx yx Z Z xy yy H H x y or E= ZH (1) MT and CSAMT sounding data, i.e. apparent resistivity (ρ a ) and phase (φ) as function of period (T) or frequency (f = T -1 ) can be calculated from each component of the impedance tensor (Z ij ), 1 ρ ; 2 ( a) ij = Zij ωµ 0 φ ij Im Z ij = tan 1 (2) Re Zij where ω = 2πf, µ 0 = 4π 10-7 (in SI unit) is freespace magnetic permeability. Data along a profile line can also be presented as pseudo-section, i.e. apparent resistivity and phase contour plot with distance on the horizontal scale and period or frequency on the vertical scale. The increasing period (or decreasing frequency) vertically downward represents increasing depth. It allows apparent resistivity and phase at many MT or CSAMT stations to be shown in a single figure and facilitates preliminary and qualitative interpretation. Quantitative interpretation can only be performed after modeling of the data using 1-D, 2-D or 3-D models to represent the subsurface resistivity distribution. Multi-dimensional MT and CSAMT modeling necessitates complete impedance tensor, while scalar impedance may result in 1-D model where resistivity varies only with depth. The investigation depth for both MT and CSAMT is usually expressed as the skin depth (in meter), δ 500 ρt or ρ δ 500 (3) f where ρ is resistivity of the homogeneous medium equivalent (see Figure 1). The skin depth serves as an approximation only and it tends to overestimate the investigation depth of an EM method. The lower frequency for CSAMT is in the range of Hz (or T between sec.). With longer periods reaching 100 sec. or more, MT investigation depth is significantly deeper than CSAMT. Typically, CSAMT is for the target depth around 1-2 km while MT is for 4-5 km. SPECIFIC CHARACTERISTICS OF CSAMT Near-Field Effect By using a transmitter (usually a grounded electric dipole, see Figure 2) to generate the source fields, CSAMT method overcomes problems associated with weak natural EM fields used in conventional audio-frequency MT surveys. However, in addition to survey logistic difficulties in providing high-powered transmitter, there are effects related to the proximity of the signal source. The assumption of uniform plane-waves as in MT implying that the wave s source is distant, is no longer valid in CSAMT. With finite transmitter-receiver distance (R), CSAMT data are likely to contain far-field and near-field data, separated by a transition zone. Generally, the near-field is defined as the zone where the separation between the transmitter and receiver is less than 3 times the skin depth for a uniform, isotropic material underlying the zone (Zonge and Hughes, 1991). The far-field is where the MT approximation is valid. Figures 3 and 4 show comparisons of CSAMT and MT sounding curves for two three-layered 1-D synthetic models with R = 2 km. In the farfield zone apparent resistivities for CSAMT and MT are identical up to a certain period. Beyond that zone, near-field zone prevails. The near-field zone for Model-1 starts at T = 0.1 sec., while in the Model-2 it starts at shorter period, i.e. approximately sec. These facts confirm that for a single transmitter-receiver separation in CSAMT sounding site, division of the data into near-field

3 and far-field depends on the skin depth. This in turn is function of the period and resistivity of the medium. The combination of the first two layers in Model-1 (1000 Ohm.m and 10 Ohm.m) are significantly lower than those of the Model-2 (100 Ohm.m and 1000 Ohm). Therefore, the far-field (R >> δ) extends to longer period in Model-1 than in the Model-2. From Figures 3 and 4 it is also obvious that in the near-field zone, the CSAMT apparent resistivity curves do not make an asymtote to the resistivity of lowest layer as in MT. Therefore, qualitative interpretation of apparent resistivity pseudosection and even quantitative 1-D modeling using MT approximation are definitely misleading. Yamashita et al. (1985) and Bartels and Jacobson (1987) advocate the correction of the near-field data such that CSAMT data can be interpreted using MT modeling tools. Grandis (2000) presented in detail the implementation of such correction. Apparent resistivity for the near-field is calculated by using a formula taking into account the effects of array geometry (Zonge and Hughes, 1991), i.e. R ( ρ a) ij = Zij or ( ρ a) ij = kn R Zij (4) 2 where R is the trasmitter-receiver distance and k n is the generalized multiplication factor for the impedance to obtain the apparent resistivity. In the transition zone from the near-field side, k n is interpolated from coefficients determined from a generalization of the homogeneous medium case. Similar approach is employed in the transition zone from the far-field side. Grandis (2000) proposed the normalized frequency for easier identification of the transition zone. Figure 5 shows the correction of the synthetic data presented in Figures 3 and 4. Corrected CSAMT data are fairly similar with the synthetic MT data, especially when the near-field effect is not too severe as in Model-1 (see Figure 5). Application of the correction technique for the field CSAMT data from Kamojang geothermal field is presented in Figure 6. Qualitative interpretation of the uncorrected pseudo-section would under-estimate the low resistivity layer. With relatively dense sounding points (100 m apart), the lateral boundary at the profile's right-end would also be missed. Source Effects Zonge and Hughes (1991) and Mitsuhata (2000) discuss the effects of the artificial source in CSAMT by using 2-D numerical model studies. The source overprint effect is illustrated in Figure 7 where geological environment near the transmitter affects CSAMT data. In addition to the subsurface resistivity distribution under the sounding point, CSAMT data also reflect the one under the transmitter. The shadow effect is another source related problems in CSAMT. Figure 8 shows the apparent resistivity pseudo-section of a line perpendicular to the transmitter as a result of 2-D numerical modeling by Zonge and Hughes (1991). Soundings at the line parallel to and on the broadside of the transmitter (for example at 7 km from the transmitter, at the right-most of Figure 8) would reflect the effect of the conductive anomaly located between trasmitter and receiver. Scalar Mode Measurement Complete tensor CSAMT measurement is possible by employing two perpendicular transmitters to provide two independent source polarisations. However, survey logistic difficulties prohibit tensor CSAMT mode for routine operations, except for research and studies for very specific purposes. In most CSAMT surveys, scalar mode measurement is always implemented. Typically, a transmitter is positioned along the x- axis and the survey grid is arranged at the broadside of the transmitter, along a line parallel to the x-axis. A multi-channel receiver records electric field simultaneously at seven dipoles (E x ) and one perpendicular magnetic field (H y ). Due to restriction in cable transmission, the distance between CSAMT soundings is limited to around 100 m (see Figure 2). In the processing stage, each electric field from a single set-up is paired with associated magnetic field from the same set-up. The results are scalar impedances Z xy = E x /H y from which apparent resistivity and phase can be calculated. Owing to redundant and spatially high density data, apparent resistivity pseudo-section may appear as a high-resolution image of the subsurface. Even so, modeling of scalar impedance data from each sounding point can only be performed by using 1-D approximation of the subsurface.

4 DISCUSSION CSAMT technique was devised to overcome disadvantages of the natural field mode, i.e. MT method. Both MT and CSAMT methods result in resistivity estimation of the subsurface. There is no reason to prefer one method over the other based only on empirical evidences related to success or failure of one of them in determining geothermal resources. CSAMT method is the most appropriate if all the following conditions are true: (i) target is relatively shallow, (ii) the study is intended as a preliminary or reconnaissance survey, (iii) the survey area suffers from both cultural and natural noise, and (iv) accessibility of survey area does not prohibit the use of high-powered electrical generator for the transmitter. For deeper target and for advanced stage of the geothermal exploration where detailed subsurface resistivity distribution shall be achieved, then MT method is more suitable. MT data quality can be improved by employing Remote Referenced MT (RRMT) measurement (e.g. Krings et al., 2007). Combination of MT and AMT (Audio-frequency MT) can still be considered to achieve an efficient compromise between delineation and detailed purposes. This is without sacrificing the completeness of the data by measuring complete tensor impedance even for AMT frequency band. In such case, dimensionality analysis can be performed on the tensor impedance data. When CSAMT is chosen, the transmitter should be placed as far as possible from the sounding points to obtain only far-field data. Far from the transmitter, the electromagnetic fields are planar, and CSAMT data can be interpreted using conventional approaches developed for MT data. Contrary to the previous requirement, there is a need to minimize transmitter-receiver separation for maximum signal. Therefore, most CSAMT data sets will contain some near-field data. However, if most of the data are far-field, then we can discard a small number of data in the nearfield and consider far-field data only in all subsequent CSAMT data analysis. Cautious CSAMT data analysis and interpretation must be performed, when they contain mostly near-field data. Although the validity of near-field CSAMT data correction is still debatable, at least it can produce better and more "correct" pseudosections. Rather than interpreting directly apparent resistivity pseudo-sections from uncorrected CSAMT data, the use of existing MT 1-D modeling to near-field corrected data is advisable. Nowadays, CSAMT 1-D inversion is available employing full-solution of the EM fields due to artificial electrical line source (e.g. Routh and Oldenburg, 1999; Wang et al., 2012). In such case, there is no need for near-field correction of CSAMT data. There are also methods for CSAMT 2-D inversion (e.g. Mitsuhata, 2000; Wang et al., 2009). However, such 2-D inversion necessitates CSAMT tensor or at least vector impedance data, both resulting from measurement using two perpendicular sources. For more practical purposes, we can consider that the profile line is oriented across a presumed geologic strike, then the scalar impedance Z xy = E x /H y may be regarded as TM-mode data. Such "incomplete" data can be inverted using available codes or softwares for MT 2-D inversion (e.g. Uchida, 1993). Figure 9 presents the result of 2-D inversion of scalar CSAMT (far-field data only) from Lahendong geothermal field. The geological interpretation of the model is beyond the scope of this paper. However, it can be observed that the 2-D subsurface resistivity distribution is very informative. Combined with results from other profiles, the model may generate meaningful interpretation related to the local geology of the survey area. CONCLUSION There is no doubt that both MT and CSAMT play significant role as main geophysical methods in geothermal exploration. Both techniques are capable to produce resistivity structure that is more or less related to the temperature regime of the subsurface. As in any other geophysical methods, data analysis and interpretation shall be performed with great cautions. In MT, the low signal to noise ratio demands data acquisition and processing that can minimize the influence of noise. The term Robust Remote Referenced MT (RRRMT) has become a standard practice in MT. The complete tensor impedance data shall be exploited as far as possible to extract

5 the dimensionality of the structures. It is wellknown that most geothermal fields related to volcanic area have very complex structures. In this regard, the use of the tensor decomposition technique is therefore essential (see Weckmann et al., 2003 for a relatively recent review). Multidimensional MT modelings (2-D or even 3-D) are also of interest in order to take into account the complicated subsurface structure of a geothermal prospect. The advantage of CSAMT, in addition to high signal to noise ratio, is the redundancy of the data along a profile line. The CSAMT appears to provide laterally high-resolution image of the subsurface. However, CSAMT data and therefore its pseduo-section suffer from erroneous effects due to finite transmitter-receiver separation. Therefore, direct qualitative interpretation from uncorrected CSAMT data is not advisable. Existing full-solution for CSAMT modelings in 1-D and 2-D are advantageous for CSAMT interpretation. These exonerate the need for nearfield correction. However, the fact that most CSAMT data are scalar prohibit us to fully take advantage from such modelling tools. For most practical use of CSAMT, it would be sufficient to consider the existing tools at hands. First, perform CSAMT only if all requirements described in the discussion above are fulfilled. Second, select only far field data or apply nearfield correction, then qualitatively analyse the resulting pseudo-sections. Third, perform MT 1-D and/or 2-D (with TM-mode data only) inversion modeling to the far-field data or near-field corrected CSAMT data. REFERENCES Bartel, L. C., and R.D. Jacobson, 1987, Results of a controlled-source audio frequency magnetotelluric survey at the Puhimau thermal area, Kilauea Volcano, Hawaii, Geophysics, 52, Cumming, W., and R. Mackie, 2010, Resistivity imaging of geothermal resources using 1D, 2D and 3D MT inversion and TDEM static shift correction illustrated by a Glass Mountain case history, Proceedings World Geothermal Congress 2010 Bali, Indonesia. Grandis, H., 2000, Koreksi efek sumber pada data Controlled-Source Audio-Magnetotellurics (CSAMT), Jurnal Teknologi Mineral, VII/1, Krings, T., O. Ritter, G. Munoz and U. Weckmann, 2007, MT robust remote reference processing revisited, 22 Kolloquium Elektromagnetische Tiefenforschung, Czech Republic. Mitsuhata, Y., 2000, 2-D electromagnetic modeling by finite-element method with a dipole source and topography, Geophysics, 65, Newman, G.A., E. Gasperikova, G.M. Hoversten and P.E. Wannamaker, 2008, Three-dimensional magnetotelluric characterization of the Coso geothermal field, Geothermics, 37-4, Routh, S.P., and D.W. Oldenburg, 1999, Inversion of controlled-source audio frequency magnetotelluric data for a horizontally layered earth, Geophysics, 64, Uchida, T., 1993, Smooth 2-D inversion of magnetotelluric data based on statistical criterion ABIC, Journal of Geomagnetism & Geoelectricity, 45, Wang, R., M. Wang, Q.Y. Di and G.J. Wang, 2009, 2- D numerical study on the effect of conductor between the transmitter and survey area in CSEM exploration, Applied geophysics, 6, Wang, R., C. Yin, M. Wang, and G. Wang, 2012, Simulated annealing for controlled-source audiofrequency magnetotelluric data inversion, Geophysics, 77, E127-E133. Wannamaker, P., 1997, Tensor CSAMT survey over the Sulphur Springs thermal area, Valles Calder, New Mexico, USA, Part 1: Implications for structure of the western Caldera, Geophysics, 62, Weckmann, U., O. Ritter and V. Haak, 2003, Images of the magnetotelluric apparent resistivity tensor, Geophysical Journal International, 155, Yamashita, M., P.G. Hallof and W.H. Pelton, 1985, CSAMT case histories with a multi-channel CSAMT system and discussion of near-field data correction, 55th SEG Annual Meeting Expanded Abstracts. Zonge, K.L., and L.J. Hughes, 1991, Controlled-source audio-frequency magnetotellurics, in Nabighian, M.N., Ed., Electromagnetic Methods in Applied Geophysics, Vol. 2, Society of Exploration Geophysicists.

6 Figure 1. Skin depth for MT and CSAMT methods as function of the period for different homogeneous medium equivalent. Figure 3. MT and CSAMT apparent resistivity sounding curves (top) and phase (bottom) for the synthetic Model-1. Figure 2. Typical lay-out of transmitter (Tx) and receiver (R) in a CSAMT survey. The survey grid is located at the broad-side of the transmitter. The scalar mode measurement at each sounding point is composed by a single component of electric and magnetic fields. Figure 4. MT and CSAMT apparent resistivity sounding curves (top) and phase (bottom) for the synthetic Model-2.

7 Figure 5. Comparison of corrected CSAMT and MT apparent resistivity sounding curves for Model-1 (left) and Model-2 (right). Figure 6. Apparent resistivity pseudo-section of CSAMT data from one line at Kamojang geothermal prospect, uncorrected (top) and corrected (bottom) using the method described by Grandis (2000). Figure 7. Apparent resistivity sounding curves for different geological environment at the transmitter location illustrating the source overprint effect modeled by Zonge and Hughes (1991).

8 PROCEEDINGS The 12TH ANNUAL INDONESIAN GEOTHERMAL ASSOCIATION MEETING & CONFERENCE pseudo section of a line perpendicular to the transmitter as a result of 2-D Figure 8. Apparent resistivity pseudo-section numerical modeling by Zonge and Hughes (1991). Soundings at the line parallel to and on the broad broad-side of the transmitter would reflect the effect of the conductive conductive anomaly located between trasmitter and receiver. Figure 9. Result of 2-D MT (TM-mode) (TM inversion modeling of scalar CSAMT far-field far data only from Lahendong geothermal field.

Magnetotelluric (MT) Method

Magnetotelluric (MT) Method Magnetotelluric (MT) Method Dr. Hendra Grandis Graduate Program in Applied Geophysics Faculty of Mining and Petroleum Engineering ITB Geophysical Methods Techniques applying physical laws (or theory) to

More information

MT Prospecting. Map Resistivity. Determine Formations. Determine Structure. Targeted Drilling

MT Prospecting. Map Resistivity. Determine Formations. Determine Structure. Targeted Drilling MT Prospecting Map Resistivity Determine Formations Determine Structure Targeted Drilling Cross-sectional interpretation before and after an MT survey of a mineral exploration prospect containing volcanic

More information

A Brief Introduction to Magnetotellurics and Controlled Source Electromagnetic Methods

A Brief Introduction to Magnetotellurics and Controlled Source Electromagnetic Methods A Brief Introduction to Magnetotellurics and Controlled Source Electromagnetic Methods Frank Morrison U.C. Berkeley With the help of: David Alumbaugh Erika Gasperikova Mike Hoversten Andrea Zirilli A few

More information

Geophysics for Geothermal Exploration

Geophysics for Geothermal Exploration Geophysics for Geothermal Exploration Dr. Hendra Grandis Geofisika - ITB Agenda Basic concept of Geophysics: : Review Geophysical signatures of a geothermal system Geophysical methods for geotermal exploration

More information

RESISTIVITY IMAGING IN EASTERN NEVADA USING THE AUDIOMAGNETOTELLURIC METHOD FOR HYDROGEOLOGIC FRAMEWORK STUDIES. Abstract.

RESISTIVITY IMAGING IN EASTERN NEVADA USING THE AUDIOMAGNETOTELLURIC METHOD FOR HYDROGEOLOGIC FRAMEWORK STUDIES. Abstract. RESISTIVITY IMAGING IN EASTERN NEVADA USING THE AUDIOMAGNETOTELLURIC METHOD FOR HYDROGEOLOGIC FRAMEWORK STUDIES Darcy K. McPhee, U.S. Geological Survey, Menlo Park, CA Louise Pellerin, Green Engineering,

More information

APPENDIX A: Magnetotelluric Data in Relation to San Pedro Mesa Structural. The San Pedro Mesa structural high (discussed in main text of paper) was

APPENDIX A: Magnetotelluric Data in Relation to San Pedro Mesa Structural. The San Pedro Mesa structural high (discussed in main text of paper) was Page of DR for GSA Special Paper 9, Chapter, Geophysical constraints APPENDIX A: Magnetotelluric Data in Relation to San Pedro Mesa Structural High The San Pedro Mesa structural high (discussed in main

More information

3D MAGNETOTELLURIC CHARACTERIZATION OF THE COSO GEOTHERMAL FIELD

3D MAGNETOTELLURIC CHARACTERIZATION OF THE COSO GEOTHERMAL FIELD PROCEEDINGS, Thirtieth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 31-February 2, 2005 SGP-TR-176 3D MAGNETOTELLURIC CHARACTERIZATION OF THE COSO GEOTHERMAL

More information

Integrated Geophysical Model for Suswa Geothermal Prospect using Resistivity, Seismics and Gravity Survey Data in Kenya

Integrated Geophysical Model for Suswa Geothermal Prospect using Resistivity, Seismics and Gravity Survey Data in Kenya Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Integrated Geophysical Model for Suswa Geothermal Prospect using Resistivity, Seismics and Gravity Survey Data in Kenya

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION Supplementary online material for Bai et al., (2). EHS3D MT data collection Broadband magnetotelluric (MT) data were recorded on profiles P, P2 and P4 in the frequency band -.5

More information

1D and 2D Inversion of the Magnetotelluric Data for Brine Bearing Structures Investigation

1D and 2D Inversion of the Magnetotelluric Data for Brine Bearing Structures Investigation 1D and 2D Inversion of the Magnetotelluric Data for Brine Bearing Structures Investigation Behrooz Oskooi *, Isa Mansoori Kermanshahi * * Institute of Geophysics, University of Tehran, Tehran, Iran. boskooi@ut.ac.ir,

More information

Magnetic and Gravity Methods for Geothermal Exploration

Magnetic and Gravity Methods for Geothermal Exploration Magnetic and Gravity Methods for Geothermal Exploration Dr. Hendra Grandis Geophysics - ITB method and survey procedure Aero- or ground magnetic (covers a large area) Schlumberger resistivity mapping and

More information

A Note on the Use of the Second Vertical Derivative (SVD) of Gravity Data with Reference to Indonesian Cases

A Note on the Use of the Second Vertical Derivative (SVD) of Gravity Data with Reference to Indonesian Cases J. Eng. Technol. Sci., Vol. 50, No. 1, 018, 17-139 17 A Note on the Use of the Second Vertical Derivative (SVD) of Gravity Data with Reference to Indonesian Cases Prihadi Sumintadireja 1, Darharta Dahrin

More information

Inversion of controlled source audio-frequency magnetotellurics data for a horizontally layered earth

Inversion of controlled source audio-frequency magnetotellurics data for a horizontally layered earth GEOPHYSICS, VOL 64, NO 6 NOVEMBER DECEMBER 1999; P 1689 1697, 8 FIGS Inversion of controlled source audio-frequency magnetotellurics data for a horizontally layered earth Partha S Routh and Douglas W Oldenburg

More information

Comparison of 1-D, 2-D and 3-D Inversion Approaches of Interpreting Electromagnetic Data of Silali Geothermal Area

Comparison of 1-D, 2-D and 3-D Inversion Approaches of Interpreting Electromagnetic Data of Silali Geothermal Area Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Comparison of 1-D, 2-D and 3-D Inversion Approaches of Interpreting Electromagnetic Data of Silali Geothermal Area Charles

More information

THREE DIMENSIONAL INVERSIONS OF MT RESISTIVITY DATA TO IMAGE GEOTHERMAL SYSTEMS: CASE STUDY, KOROSI GEOTHERMAL PROSPECT.

THREE DIMENSIONAL INVERSIONS OF MT RESISTIVITY DATA TO IMAGE GEOTHERMAL SYSTEMS: CASE STUDY, KOROSI GEOTHERMAL PROSPECT. Proceedings, 6 th African Rift Geothermal Conference Addis Ababa, Ethiopia, 2 nd 4 th November 2016 THREE DIMENSIONAL INVERSIONS OF MT RESISTIVITY DATA TO IMAGE GEOTHERMAL SYSTEMS: CASE STUDY, KOROSI GEOTHERMAL

More information

Resistivity structure of Sumikawa geothermal field, northeastern Japan, obtained from magnetotelluric data. Toshihiro Uchida

Resistivity structure of Sumikawa geothermal field, northeastern Japan, obtained from magnetotelluric data. Toshihiro Uchida Resistivity structure of Sumikawa geothermal field, northeastern Japan, obtained from magnetotelluric data Toshihiro Geological Survey of Japan 1-1-3 Higashi, Tsukuba, Ibaraki 35, Japan ABSTRACT: Resistivity

More information

A resolution comparison of horizontal and vertical magnetic transfer functions

A resolution comparison of horizontal and vertical magnetic transfer functions Journal of the Earth and Space Physics, Vol. 40, No. 3, 2014, P. 47-53 A resolution comparison of horizontal and vertical magnetic transfer functions Habibian Dehkordi, B. 1 and Oskooi, B. 2 1 Assitant

More information

Data Repository Comeau et al., (2015)

Data Repository Comeau et al., (2015) Data Repository 2015087 Comeau et al., (2015) 1. Magnetotelluric data Owing to the remote location, and large distance from sources of electromagnetic noise, most stations were characterized by high quality

More information

Numerical simulation of magnetotelluric fields at Stromboli

Numerical simulation of magnetotelluric fields at Stromboli Numerical simulation of magnetotelluric fields at Stromboli A. Franke, S. Kütter, R.-U. Börner and K. Spitzer TU Bergakademie Freiberg, Germany 1 Summary Stromboli is a small volcanic island in the Mediterranean

More information

Three Dimensional Inversions of MT Resistivity Data to Image Geothermal Systems: Case Study, Korosi Geothermal Prospect

Three Dimensional Inversions of MT Resistivity Data to Image Geothermal Systems: Case Study, Korosi Geothermal Prospect PROCEEDINGS, 42nd Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 13-15, 2017 SGP-TR-212 Three Dimensional Inversions of MT Resistivity Data to Image Geothermal

More information

Exploration of Geothermal High Enthalpy Resources using Magnetotellurics an Example from Chile

Exploration of Geothermal High Enthalpy Resources using Magnetotellurics an Example from Chile Exploration of Geothermal High Enthalpy Resources using Magnetotellurics an Example from Chile Ulrich Kalberkamp, Federal Institute for Geosciences and Natural Resources (BGR), Stilleweg 2, 30655 Hannover,

More information

Keywords: geophysics field camp, low to moderate, temperature, geothermal system, Mt Pancar, Indonesia

Keywords: geophysics field camp, low to moderate, temperature, geothermal system, Mt Pancar, Indonesia Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Geophysics Field Camp (GFC): A Student Project Aimed to Investigate the Low to Moderate Temperature Geothermal System in

More information

Subsurface Conductivity Obtained from DC Railway Signal Propagation with a Dipole Model

Subsurface Conductivity Obtained from DC Railway Signal Propagation with a Dipole Model Subsurface Conductivity Obtained from DC Railway Signal Propagation with a Dipole Model Anne Neska Institute of Geophysics PAS, Ul. Ks. Janusza 64, 01-452 Warszawa, Poland, email: anne@igf.edu.pl Abstract

More information

An Experimental Magnetotelluric Short Sounding During High Sun Spot Activity Cycle in 2013

An Experimental Magnetotelluric Short Sounding During High Sun Spot Activity Cycle in 2013 Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 An Experimental Magnetotelluric Short Sounding During High Sun Spot Activity Cycle in 2013 Lendriadi AGUNG and Imam B.

More information

Magnetotelluric and Controlled-Source Electromagnetic Study of Aquistore CO 2 Sequestration Site, near Estevan, Saskatchewan

Magnetotelluric and Controlled-Source Electromagnetic Study of Aquistore CO 2 Sequestration Site, near Estevan, Saskatchewan Magnetotelluric and Controlled-Source Electromagnetic Study of Aquistore CO 2 Sequestration Site, near Estevan, Saskatchewan Joe, B, McLeod 1, Ian, J, Ferguson 1, Jim Craven 2, Brian Roberts 2, Eric Roots

More information

Determination of EM Coupling Effect in Dipole-Dipole Method : Complex Image Approach

Determination of EM Coupling Effect in Dipole-Dipole Method : Complex Image Approach 6 IJCSNS International Journal of Computer Science and Network Security, VOL. No.3, March Determination of EM Coupling Effect in Dipole-Dipole Method : Complex Image Approach Mester Sitepu Universitas

More information

Peter R. Milligan. Stephan Thiel. Geoscience Australia. The University of Adelaide. MT acquisition & processing

Peter R. Milligan. Stephan Thiel. Geoscience Australia. The University of Adelaide. MT acquisition & processing Magnetotelluric acquisition & processing, with examples from the Gawler Craton, Curnamona Province and Curnamona- Gawler Link transects in South Australia Peter R. Milligan Geoscience Australia Stephan

More information

Vertical and horizontal resolution considerations for a joint 3D CSEM and MT inversion

Vertical and horizontal resolution considerations for a joint 3D CSEM and MT inversion Antony PRICE*, Total E&P and Don WATTS, WesternGeco Electromagnetics Summary To further explore the potential data content and inherent limitations of a detailed 3D Controlled Source ElectroMagnetic and

More information

Verification and Validation Calculations Using the STAR Geophysical Postprocessor Suite

Verification and Validation Calculations Using the STAR Geophysical Postprocessor Suite SAIC-03/1040 Final Report Verification and Validation Calculations Using the STAR Geophysical Postprocessor Suite Prepared by: J. W. Pritchett December 2003 Verification and Validation Calculations Using

More information

Geophysics Course Introduction to DC Resistivity

Geophysics Course Introduction to DC Resistivity NORAD supported project in MRRD covering Capacity Building and Institutional Cooperation in the field of Hydrogeology for Faryab Province Afghanistan Geophysics Course Introduction to DC Resistivity By

More information

SEAM2: 3D NON-SEISMIC MODELING OF A COMPLEX MIDDLE EAST O&G PROSPECT

SEAM2: 3D NON-SEISMIC MODELING OF A COMPLEX MIDDLE EAST O&G PROSPECT 6th International Symposium on Three-Dimensional Electromagnetics Berkeley, California, USA, March 28-3, 217 SEAM2: 3D NON-SEISMIC MODELING OF A COMPLEX MIDDLE EAST O&G PROSPECT M. D. Watts 1, R. L. Mackie

More information

VLF -MT Survey around Nakadake crater at Aso Volcano

VLF -MT Survey around Nakadake crater at Aso Volcano IOP Conference Series: Earth and Environmental Science OPEN ACCESS VLF -MT Survey around Nakadake crater at Aso Volcano To cite this article: A Harja et al 2016 IOP Conf. Ser.: Earth Environ. Sci. 29 012016

More information

Early Analogue Modeling Experiments and Related Studies to Today s Problems of Geo-electromagnetic Exploration

Early Analogue Modeling Experiments and Related Studies to Today s Problems of Geo-electromagnetic Exploration Journal of Earth Science, Vol. 20, No. 3, p. 618 625, June 2009 ISSN 1674-487X Printed in China DOI: 10.1007/s12583-009-0051-y Early Analogue Modeling Experiments and Related Studies to Today s Problems

More information

Two-dimensional Inversion of Resistivity and IP Data with Topography

Two-dimensional Inversion of Resistivity and IP Data with Topography Two-dimensional Inversion of Resistivity and IP Data with Topography a) Smooth-Model Resistivity (ohm-m) Elevation (m) b) Calculated Apparent Resistivity (ohm-m) 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9

More information

X004 3D CSEM Inversion Strategy - An Example Offshore West of Greenland

X004 3D CSEM Inversion Strategy - An Example Offshore West of Greenland X004 3D CSEM Inversion Strategy - An Example Offshore West of Greenland A. Lovatini* (WesternGeco EM - Geosystem), M.D. Watts (WesternGeco EM - Geosystem), K. Umbach (EnCana Corp.) & A. Ferster (EnCana

More information

Introducing the Stratagem EH4. Electrical Conductivity Imaging System Hybrid-Source Magnetotellurics

Introducing the Stratagem EH4. Electrical Conductivity Imaging System Hybrid-Source Magnetotellurics Introducing the Stratagem EH4 Electrical Conductivity Imaging System Hybrid-Source Magnetotellurics Stratagem EH4 Hybrid-Source Magnetotellurics Frequency range of 10 Hz to 90k Hz Approx. depth of investigation

More information

THREE DIMENSIONAL MODELING OF GEOELECTRICAL STRUCTURE BASED ON MT AND TDEM DATA IN MORI GEOTHERMAL FIELD, HOKKAIDO, JAPAN

THREE DIMENSIONAL MODELING OF GEOELECTRICAL STRUCTURE BASED ON MT AND TDEM DATA IN MORI GEOTHERMAL FIELD, HOKKAIDO, JAPAN THREE DIMENSIONAL MODELING OF GEOELECTRICAL STRUCTURE BASED ON MT AND TDEM DATA IN MORI GEOTHERMAL FIELD, HOKKAIDO, JAPAN Tatsuya Kajiwara 1, Tohru Mogi 2, Elena Fomenko 3 and Sachio Ehara 4 1 JMC Geothermal

More information

CSAMT exploration at Sellafield: Characterization of a potential radioactive waste disposal site

CSAMT exploration at Sellafield: Characterization of a potential radioactive waste disposal site GEOPHYSICS, VOL. 65, NO. 4 (JULY AUGUST 2000); P. 1070 1079, 10 FIGS. CSAMT exploration at Sellafield: Characterization of a potential radioactive waste disposal site Martyn J. Unsworth, Xinyou Lu, and

More information

C5 Magnetic exploration methods data analysis techniques

C5 Magnetic exploration methods data analysis techniques C5 Magnetic exploration methods data analysis techniques C5.1 Data processing and corrections After magnetic field data have been collected a number of corrections are applied to simplify the interpretation.

More information

Geoelectricity. ieso 2010

Geoelectricity. ieso 2010 Geoelectricity ieso 2010 1 RESISTIVITY SURVEY AT VENETO VILLA GRITTI AT THE TOWN OF TREVISO (VENETO REGION) The survey was carried out to verify the underground presence of the fondations of a rustic building.

More information

Downloaded 08/29/13 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 08/29/13 to Redistribution subject to SEG license or copyright; see Terms of Use at New approach to 3D inversion of MCSEM and MMT data using multinary model transform Alexander V. Gribenko and Michael S. Zhdanov, University of Utah and TechnoImaging SUMMARY Marine controlled-source electromagnetic

More information

Hamed Aber 1 : Islamic Azad University, Science and Research branch, Tehran, Iran. Mir Sattar Meshin chi asl 2 :

Hamed Aber 1 : Islamic Azad University, Science and Research branch, Tehran, Iran. Mir Sattar Meshin chi asl 2 : Present a Proper Pattern for Choose Best Electrode Array Based on Geological Structure Investigating in Geoelectrical Tomography, in order to Get the Highest Resolution Image of the Subsurface Hamed Aber

More information

Three-dimensional geophysical modelling of the Alalobeda geothermal field

Three-dimensional geophysical modelling of the Alalobeda geothermal field Three-dimensional geophysical modelling of the Alalobeda geothermal field Daniele Rizzello 1, Egidio Armadillo 2, Claudio Pasqua 2,3, Massimo Verdoya 2, Solomon Kebede 4, Andarge Mengiste 4, Getenesh H/Giorgis

More information

POTASH DRAGON CHILE GEOPHYSICAL SURVEY TRANSIENT ELECTROMAGNETIC (TEM) METHOD. LLAMARA and SOLIDA PROJECTS SALAR DE LLAMARA, IQUIQUE, REGION I, CHILE

POTASH DRAGON CHILE GEOPHYSICAL SURVEY TRANSIENT ELECTROMAGNETIC (TEM) METHOD. LLAMARA and SOLIDA PROJECTS SALAR DE LLAMARA, IQUIQUE, REGION I, CHILE POTASH DRAGON CHILE GEOPHYSICAL SURVEY TRANSIENT ELECTROMAGNETIC (TEM) METHOD LLAMARA and SOLIDA PROJECTS SALAR DE LLAMARA, IQUIQUE, REGION I, CHILE OCTOBER 2012 CONTENT Page I INTRODUCTION 1 II FIELD

More information

CHARACTERIZATION OF SOIL PROFILE OF DHAKA CITY USING ELECTRICAL RESISTIVITY TOMOGRAPHY (ERT)

CHARACTERIZATION OF SOIL PROFILE OF DHAKA CITY USING ELECTRICAL RESISTIVITY TOMOGRAPHY (ERT) CHARACTERIZATION OF SOIL PROFILE OF DHAKA CITY USING ELECTRICAL RESISTIVITY TOMOGRAPHY (ERT) Mehedi Ahmed ANSARY 1, B.S. Pushpendue BISWAS 2 and Abul KHAIR 3 1 Professor, Department of Civil Engineering

More information

Anisotropic 2.5D Inversion of Towed Streamer EM Data from Three North Sea Fields Using Parallel Adaptive Finite Elements

Anisotropic 2.5D Inversion of Towed Streamer EM Data from Three North Sea Fields Using Parallel Adaptive Finite Elements Anisotropic 2.5D Inversion of Towed Streamer EM Data from Three North Sea Fields Using Parallel Adaptive Finite Elements K. Key (Scripps Institution of Oceanography), Z. Du* (PGS), J. Mattsson (PGS), A.

More information

Titiribi Gold: Sleeping Giant in the Colombian Cordilleras (a Geophysical Study)

Titiribi Gold: Sleeping Giant in the Colombian Cordilleras (a Geophysical Study) Titiribi Gold: Sleeping Giant in the Colombian Cordilleras (a Geophysical Study) Emmett Van Reed, P.Gp. Zonge Engineering & Research Amit Tripathi, Ph.D. Terra Resources Development Technologies Jennifer

More information

ARRAYS FOR TENSOR MEASUREMENTS OF THE ELECTRIC FIELD. Abstract

ARRAYS FOR TENSOR MEASUREMENTS OF THE ELECTRIC FIELD. Abstract ARRAYS FOR TENSOR MEASUREMENTS OF THE ELECTRIC FIELD Aleksandr Mousatov, Mexican Petroleum Institute, Mexico E-mail: amousat@imp.mx Evgueni Pervago, Mexican Petroleum Institute, Mexico E-mail: epervago@imp.mx

More information

LIST OF FIGURES APPENDICES

LIST OF FIGURES APPENDICES RESISTIVITY / INDUCED POLARIZATION SURVEY EL PORVENIR PROJECT MUNICIPALITY OF REMEDIOS, ANTIOQUIA, COLOMBIA LOGISTICS REPORT M-17427 APRIL 2017 TABLE OF CONTENTS Abstract... 1 1. The Mandate... 2 2. El

More information

. Slide 1. Geological Survey of Ethiopia, P.O.Box 2302, Addis Ababa, Ethiopia

. Slide 1. Geological Survey of Ethiopia, P.O.Box 2302, Addis Ababa, Ethiopia The electrical resistivity structure of the Southern Main Ethiopian Rift (around the Lake Abaya geothermal prospect) revealed by magneto-telluric telluric imaging Mohammednur Desissa and Yohanes Lema Geological

More information

Electrical prospecting involves detection of surface effects produced by electrical current flow in the ground.

Electrical prospecting involves detection of surface effects produced by electrical current flow in the ground. Electrical Surveys in Geophysics Electrical prospecting involves detection of surface effects produced by electrical current flow in the ground. Electrical resistivity method Induced polarization (IP)

More information

RESISTIVITY SURVEYING AND ELECTROMAGNETIC METHODS

RESISTIVITY SURVEYING AND ELECTROMAGNETIC METHODS Presented at Short Course VII on Surface Exploration for Geothermal Resources, organized by UNU-GTP and LaGeo, in Santa Tecla and Ahuachapán, El Salvador, March 14-22, 2015. GEOTHERMAL TRAINING PROGRAMME

More information

Towed Streamer EM data from Barents Sea, Norway

Towed Streamer EM data from Barents Sea, Norway Towed Streamer EM data from Barents Sea, Norway Anwar Bhuiyan*, Eivind Vesterås and Allan McKay, PGS Summary The measured Towed Streamer EM data from a survey in the Barents Sea, undertaken in the Norwegian

More information

Western EARS Geothermal Geophysics Context Mar-2016 William Cumming

Western EARS Geothermal Geophysics Context Mar-2016 William Cumming Western EARS Geothermal Geophysics Context Mar-2016 William Cumming Cumming Geoscience, Santa Rosa CA wcumming@wcumming.com Office: +1-707-546-1245 Mobile: +1-707-483-7959 Skype: wcumming.com Kibiro Geophysics

More information

Pseudo-seismic wavelet transformation of transient electromagnetic response in engineering geology exploration

Pseudo-seismic wavelet transformation of transient electromagnetic response in engineering geology exploration GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L645, doi:.29/27gl36, 27 Pseudo-seismic wavelet transformation of transient electromagnetic response in engineering geology exploration G. Q. Xue, Y. J. Yan, 2 and

More information

Distortion Effects on Magnetotelluric Sounding Data Investigated by 3D Modeling of High-Resolution Topography

Distortion Effects on Magnetotelluric Sounding Data Investigated by 3D Modeling of High-Resolution Topography GRC Transactions, Vol. 37, 2013 Distortion Effects on Magnetotelluric Sounding Data Investigated by 3D Modeling of High-Resolution Topography Mitchel A. Stark 1, Wolfgang Soyer 2, Stephen Hallinan 2, and

More information

Modeling of 3D MCSEM and Sensitivity Analysis

Modeling of 3D MCSEM and Sensitivity Analysis PIERS ONLINE, VOL. 3, NO. 5, 2007 641 Modeling of 3D MCSEM and Sensitivity Analysis Zhanxiang He 1, 2, Zhigang Wang 2, Gang Yu 3, Kurt Strack 3, and Haiying Liu 2 1 Department of Geosciences, University

More information

GLE 594: An introduction to applied geophysics

GLE 594: An introduction to applied geophysics GLE 594: An introduction to applied geophysics Electrical Resistivity Methods Fall 2004 Earth Properties and Basic Theory Reading Today : 207-218 Next Lecture : 218-228 1 Introduction Link resistivity

More information

Magnetotelluric tensor decomposition: Part II, Examples of a basic procedure

Magnetotelluric tensor decomposition: Part II, Examples of a basic procedure GEOPHYSICS, VOL. 63, NO. 6 (NOVEMBER-DECEMBER 1998); P. 1898 1907, 5 FIGS. Magnetotelluric tensor decomposition: Part II, Examples of a basic procedure F. E. M. (Ted) Lilley ABSTRACT The decomposition

More information

First Technical Report Geophysical experiments near Kajiado town

First Technical Report Geophysical experiments near Kajiado town First Technical Report Geophysical experiments near Kajiado town Synthesis Report of First case study within the ISGEAG VIA Water project By: Michel Groen *), Harry Rolf **), and Ammon Muiti ***) *) Acacia

More information

Available online Journal of Scientific and Engineering Research, 2016, 3(2):1-7. Research Article

Available online   Journal of Scientific and Engineering Research, 2016, 3(2):1-7. Research Article Available online www.jsaer.com, 2016, 3(2):1-7 Research Article ISSN: 2394-2630 CODEN(USA): JSERBR Assessment of the Reliability of Magnetic Method to Delineate Geologic Features in a Basement Complex:

More information

Effects of Dipping Layers on TEM 1D-Inversion

Effects of Dipping Layers on TEM 1D-Inversion Effects of Dipping Layers on TEM 1D-Inversion S. Hölz, T. Hiller & H. Burkhardt Angew. Geowissenschaften, FR Geophysik TU Berlin Introduction In the interpretation of central-loop or coincident-loop TEM-measurements,

More information

DEEP GROUND-WATER EXPLORATION IN THE FLAGSTAFF, ARIZONA AREA USING CSAMT

DEEP GROUND-WATER EXPLORATION IN THE FLAGSTAFF, ARIZONA AREA USING CSAMT DEEP GROUND-WATER EXPLORATION IN THE FLAGSTAFF, ARIZONA AREA USING CSAMT Norman R. Carlson Zonge Engineering & Research Organization, Inc. 3322 E Fort Lowell Rd Tucson, Arizona 85716 Tel: (520)327-5501

More information

National Geophysical Survey Science Opportunities

National Geophysical Survey Science Opportunities National Geophysical Survey Science Opportunities Anton Ziolkowski Professor of Petroleum Geoscience University of Edinburgh BGS Keyworth 4 April 2016 Outline Brief biography Science opportunities High

More information

A resistivity cross-section of Usu volcano, Hokkaido, Japan, by audiomagnetotelluric soundings

A resistivity cross-section of Usu volcano, Hokkaido, Japan, by audiomagnetotelluric soundings Earth Planets Space, 50, 339 346, 1998 A resistivity cross-section of Usu volcano, Hokkaido, Japan, by audiomagnetotelluric soundings Y. Ogawa 1, N. Matsushima 1, H. Oshima 2, S. Takakura 1, M. Utsugi

More information

Lawrence Berkeley Laboratory, MS 74R316C, One Cyclotron Road, Berkeley CA 94720, USA.

Lawrence Berkeley Laboratory, MS 74R316C, One Cyclotron Road, Berkeley CA 94720, USA. Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 The Importance of Full Impedance Tensor Analysis for 3D Magnetotelluric Imaging the Roots of High Temperature Geothermal

More information

Multi-source inversion of TEM data: with field applications to Mt. Milligan

Multi-source inversion of TEM data: with field applications to Mt. Milligan Multi-source inversion of TEM data: with field applications to Mt. Milligan D. Oldenburg, E. Haber, D. Yang Geophysical Inversion Facility, University of British Columbia, Vancouver, BC, Canada Summary

More information

Summary. Introduction

Summary. Introduction Noel Black*, Glenn A. Wilson, TechnoImaging, Alexander V. Gribenko and Michael S. Zhdanov, TechnoImaging and The University of Utah Summary Recent studies have inferred the feasibility of time-lapse controlled-source

More information

SEG Houston 2009 International Exposition and Annual Meeting SUMMARY

SEG Houston 2009 International Exposition and Annual Meeting SUMMARY Simultaneous joint inversion of MT and CSEM data using a multiplicative cost function Aria Abubakar, Maokun Li, Jianguo Liu and Tarek M. Habashy, Schlumberger-Doll Research, Cambridge, MA, USA SUMMARY

More information

15. MATHEMATICAL ALGORITHMS

15. MATHEMATICAL ALGORITHMS 15. MATHEMATICAL ALGORITHMS 15.1 STANDARD ERROR OF THE MEAN (SEM)... 2 15.2 3-POINT DECOUPLING... 3 15.3 A NOTE ON RESISTIVITY CALCULATIONS... 4 15.4 RESISTIVITY CALCULATIONS... 4 DIPOLE - DIPOLE...4 POLE

More information

Interpretation of out of loop data in large fixed-loop TEM surveys Les P. Beard, Zonge International, Tucson, Arizona, USA

Interpretation of out of loop data in large fixed-loop TEM surveys Les P. Beard, Zonge International, Tucson, Arizona, USA Interpretation of out of loop data in large fixed-loop TEM surveys Les P. Beard, Zonge International, Tucson, Arizona, USA Summary Interpretation of out-of-loop data from fixed-loop transient EM surveys

More information

3D MAGNETOTELLURIC SURVEY AT THE YANAIZU-NISHIYAMA GEOTHERMAL FIELD, NORTHERN JAPAN

3D MAGNETOTELLURIC SURVEY AT THE YANAIZU-NISHIYAMA GEOTHERMAL FIELD, NORTHERN JAPAN 3D MAGNETOTELLURIC SURVEY AT THE YANAIZU-NISHIYAMA GEOTHERMAL FIELD, NORTHERN JAPAN Toshihiro Uchida 1, Shinichi Takakura 1, Takumi Ueda 1, Masaho Adachi 2, Hitoshi Ozeki 2, Kunikazu Kamada 3, Tatsuya

More information

Full Tensor Gradient of Simulated Gravity Data for Prospect Scale Delineation

Full Tensor Gradient of Simulated Gravity Data for Prospect Scale Delineation J. Math. Fund. Sci., Vol. 46, No., 014, 107-14 107 Full Tensor Gradient of Simulated Gravity Data for Prospect Scale Delineation Hendra Grandis* & Darharta Dahrin Applied Geophysics Research Group, Faculty

More information

OUTLINE COURSE DESCRIPTION

OUTLINE COURSE DESCRIPTION GEOL 482.3 OUTLINE ELECTRICAL METHODS GEOL 482.3 email: jim.merriam@usask.ca ELECTRICAL METHODS IN GEOPHYSICAL PROSPECTING COURSE DESCRIPTION This course reviews the background theory for the electromagnetic

More information

The Role of Magnetotellurics in Geothermal Exploration

The Role of Magnetotellurics in Geothermal Exploration The Role of Magnetotellurics in Geothermal Exploration Adele Manzella CNR - Via Moruzzi 1 56124 PISA, Italy manzella@igg.cnr.it Foreword MT is one of the most used geophysical methods for geothermal exploration.

More information

Main Menu. Douglas Oldenburg University of British Columbia Vancouver, BC, Canada

Main Menu. Douglas Oldenburg University of British Columbia Vancouver, BC, Canada Using ERA low frequency E-field profiling and UBC 3D frequency-domain inversion to delineate and discover a mineralized zone in Porcupine district, Ontario, Canada. Vlad Kaminski* University of British

More information

Tu Efficient 3D MT Inversion Using Finite-difference Time-domain Modelling

Tu Efficient 3D MT Inversion Using Finite-difference Time-domain Modelling Tu 11 16 Efficient 3D MT Inversion Using Finite-difference Time-domain Modelling S. de la Kethulle de Ryhove* (Electromagnetic Geoservices ASA), J.P. Morten (Electromagnetic Geoservices ASA) & K. Kumar

More information

We N Geophysical Near-surface Characterization for Static Corrections: Multi-physics Survey in Reggane Field, Algeria

We N Geophysical Near-surface Characterization for Static Corrections: Multi-physics Survey in Reggane Field, Algeria We N114 01 Geophysical Near-surface Characterization for Static Corrections: Multi-physics Survey in Reggane Field, Algeria A. Pineda* (Repsol), S. Gallo (CGG) & H. Harkas (GRN Sonatrach) SUMMARY We are

More information

Geophysical Exploration in Water Resources Assessment. John Mundell, P.E., L.P.G., P.G. Ryan Brumbaugh, L.P.G. Mundell & Associates, Inc.

Geophysical Exploration in Water Resources Assessment. John Mundell, P.E., L.P.G., P.G. Ryan Brumbaugh, L.P.G. Mundell & Associates, Inc. Geophysical Exploration in Water Resources Assessment John Mundell, P.E., L.P.G., P.G. Ryan Brumbaugh, L.P.G. Mundell & Associates, Inc. Presentation Objective Introduce the use of geophysical survey methods

More information

3D INVERSION OF MT DATA IN GEOTHERMAL EXPLORATION: A WORKFLOW AND APPLICATION TO HENGILL, ICELAND

3D INVERSION OF MT DATA IN GEOTHERMAL EXPLORATION: A WORKFLOW AND APPLICATION TO HENGILL, ICELAND PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30 - February 1, 2012 SGP-TR-194 3D INVERSION OF MT DATA IN GEOTHERMAL EXPLORATION:

More information

Imaging Reservoir Structure of Mt. Pancar Geothermal Prospect Using Audio-Frequency Magnetotelluric (AMT) and Gravity Technology

Imaging Reservoir Structure of Mt. Pancar Geothermal Prospect Using Audio-Frequency Magnetotelluric (AMT) and Gravity Technology PROCEEDINGS, 42nd Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 13-15, 2017 SGP-TR-212 Imaging Reservoir Structure of Mt. Pancar Geothermal Prospect Using

More information

Geophysical Investigation of the Old Gaborone Dumpsite, Botswana SHEMANG, E M; MOLWALEFHE, L; CHAOKA, TR; MOSWEU E; NONDO, M

Geophysical Investigation of the Old Gaborone Dumpsite, Botswana SHEMANG, E M; MOLWALEFHE, L; CHAOKA, TR; MOSWEU E; NONDO, M JASEM ISSN 1119-8362 All rights reserved Full-text Available Online at www.bioline.org.br/ja J. Appl. Sci. Environ. Mgt. September, 2006 Vol. 10 (3) 87-92 Geophysical Investigation of the Old Gaborone

More information

Optimizing Geophysical Inversions for Archean Orogenic Gold Settings

Optimizing Geophysical Inversions for Archean Orogenic Gold Settings Geophysical Inversion and Modeling Optimizing Geophysical Inversions for Archean Orogenic Gold Settings 1. University of British Columbia Paper 103 Mitchinson, D. E. [1], Phillips, N. D. [1] ABSTRACT Geophysical

More information

Effects of galvanic distortions on magnetotelluric data: Interpretation and its correction using deep electrical data

Effects of galvanic distortions on magnetotelluric data: Interpretation and its correction using deep electrical data Effects of galvanic distortions on magnetotelluric data: Interpretation and its correction using deep electrical data Jimmy Stephen 1, S G Gokarn 2, C Manoj 1 and S B Singh 1 1 National Geophysical Research

More information

BULLS-EYE! - SIMPLE RESISTIVITY IMAGING TO RELIABLY LOCATE THE GEOTHERMAL RESERVOIR

BULLS-EYE! - SIMPLE RESISTIVITY IMAGING TO RELIABLY LOCATE THE GEOTHERMAL RESERVOIR BULLS-EYE! - SIMPLE RESISTIVITY IMAGING TO RELIABLY LOCATE THE GEOTHERMAL RESERVOIR Errol Anderson, Daryl Crosby and Greg Ussher PB Power GENZL Division, P O Box 3935, Auckland, New Zealand Key Words:

More information

Numerical modeling of Z-TEM (airborne AFMAG) responses to guide exploration strategies Bob Lo*, Geotech Ltd., Michael Zang, Exploration Syndicate Inc.

Numerical modeling of Z-TEM (airborne AFMAG) responses to guide exploration strategies Bob Lo*, Geotech Ltd., Michael Zang, Exploration Syndicate Inc. Bob Lo*, Geotech Ltd., Michael Zang, Exploration Syndicate Inc. Summary Mineral exploration using the Z-TEM or Z-Axis Tipper Electromagnetics, an airborne AFMAG system, has been guided by numerical modeling

More information

INTRODUCTION TO APPLIED GEOPHYSICS

INTRODUCTION TO APPLIED GEOPHYSICS INTRODUCTION TO APPLIED GEOPHYSICS EXPLORING THE SHALL0W SUBSURFACE H. Robert Burger Anne F. Sheehan Craig H.Jones VERSITY OF COLORADO VERSITY OF COLORADO W. W. NORTON & COMPANY NEW YORK LONDON Contents

More information

Inductive source induced polarization David Marchant, Eldad Haber and Douglas W. Oldenburg, University of British Columbia

Inductive source induced polarization David Marchant, Eldad Haber and Douglas W. Oldenburg, University of British Columbia David Marchant, Eldad Haber and Douglas W. Oldenburg, University of British Columbia Downloaded /9/3 to 37.8.7.. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/

More information

RESISTIVITY STRUCTURE OF PAKA GEOTHERMAL PROSPECT IN KENYA

RESISTIVITY STRUCTURE OF PAKA GEOTHERMAL PROSPECT IN KENYA GEOTHERMAL TRAINING PROGRAMME Reports 2011 Orkustofnun, Grensasvegur 9, Number 26 IS-108 Reykjavik, Iceland RESISTIVITY STRUCTURE OF PAKA GEOTHERMAL PROSPECT IN KENYA Raymond M. Mwakirani Geothermal Development

More information

Physics-Based Characterization of UXO from Multi-Component TEM Data. Abstract

Physics-Based Characterization of UXO from Multi-Component TEM Data. Abstract Physics-Based Characterization of UXO from Multi-Component TEM Data Scott C. MacInnes, Donald D. Snyder, and Kenneth L. Zonge Zonge Engineering & Research Organization 3322 E Fort Lowell Rd Tucson, AZ,

More information

Optimization of galvanic TEM-system geometry to decrease induced polarization effect on transient electromagnetic data

Optimization of galvanic TEM-system geometry to decrease induced polarization effect on transient electromagnetic data IAGA WG. on Electromagnetic Induction in the Earth Extended Abstract 9 th Workshop S._E0 Optimization of galvanic TEM-system geometry to decrease induced polarization effect on transient electromagnetic

More information

Z i = a X i + b Y i + r i (1) 22. Kolloquium Elektromagnetische Tiefenforschung, Hotel Maxičky, Děčín, Czech Republic, October 1-5,

Z i = a X i + b Y i + r i (1) 22. Kolloquium Elektromagnetische Tiefenforschung, Hotel Maxičky, Děčín, Czech Republic, October 1-5, MT Robust Remote Reference Processing revisited Thomas Krings 1,2, Oliver Ritter 2, Gerard Muñoz 2, Ute Weckmann 2 1 Westfälische Wilhelms-Universität Münster 2 GeoForschungsZentrum Potsdam contact: t

More information

Using seismic guided EM inversion to explore a complex geological area: An application to the Kraken and Bressay heavy oil discoveries, North Sea

Using seismic guided EM inversion to explore a complex geological area: An application to the Kraken and Bressay heavy oil discoveries, North Sea Using seismic guided EM inversion to explore a complex geological area: An application to the Kraken and Bressay heavy oil discoveries, North Sea Zhijun Du*, PGS Summary The integrated analysis of controlled

More information

The relative influence of different types of magnetotelluric data on joint inversions

The relative influence of different types of magnetotelluric data on joint inversions Earth Planets Space, 55, 243 248, 23 The relative influence of different types of magnetotelluric data on joint inversions Anna Gabàs and Alex Marcuello GRC Geodinàmica i Anàlisi de Conques, Departament

More information

Geophysics for Environmental and Geotechnical Applications

Geophysics for Environmental and Geotechnical Applications Geophysics for Environmental and Geotechnical Applications Dr. Katherine Grote University of Wisconsin Eau Claire Why Use Geophysics? Improve the quality of site characterization (higher resolution and

More information

INTERGRATED GEOPHYSICAL METHODS USED TO SITE HIGH PRODUCER GEOTHERMAL WELLS

INTERGRATED GEOPHYSICAL METHODS USED TO SITE HIGH PRODUCER GEOTHERMAL WELLS Presented at Short Course VII on Exploration for Geothermal Resources, organized by UNU-GTP, GDC and KenGen, at Lake Bogoria and Lake Naivasha, Kenya, Oct. 27 Nov. 18, 2012. GEOTHERMAL TRAINING PROGRAMME

More information

Magnetotelluric Survey in an Extremely Noisy Environment at the Pohang Low-Enthalpy Geothermal Area, Korea

Magnetotelluric Survey in an Extremely Noisy Environment at the Pohang Low-Enthalpy Geothermal Area, Korea Proceedings World Geothermal Congress 2005 Antalya, Turkey, 24-29 April 2005 Magnetotelluric Survey in an Extremely Noisy Environment at the Pohang Low-Enthalpy Geothermal Area, Korea Toshihiro Uchida

More information

Achieving depth resolution with gradient array survey data through transient electromagnetic inversion

Achieving depth resolution with gradient array survey data through transient electromagnetic inversion Achieving depth resolution with gradient array survey data through transient electromagnetic inversion Downloaded /1/17 to 128.189.118.. Redistribution subject to SEG license or copyright; see Terms of

More information

High Resolution Geophysics: A Better View of the Subsurface. By John Jansen, P.G., Ph.D., Aquifer Science and Technology

High Resolution Geophysics: A Better View of the Subsurface. By John Jansen, P.G., Ph.D., Aquifer Science and Technology High Resolution Geophysics: A Better View of the Subsurface By John Jansen, P.G., Ph.D., Aquifer Science and Technology Geologist Use Only Part of the Information Available To Them Most Geologist rely

More information