Utilizing of Geophysical Method for Geothermal Exploration In Aceh Besar (Indonesia)

Similar documents
2-D Resistivity Study: The Horizontal Resolution Improvement by Introducing the Enhancing Horizontal Resolution (EHR) Technique

The Preliminary Study of Meteorite Impact Crater at Bukit Bunuh, Lenggong

The Effect of 2-D Resistivity on Magnetic Field

Case Study: Shallow Subsurface Geology Mapping Using 2-D Resistivity Imaging with EHR Technique

Integration of Seismic Refraction and 2D Electrical Resistivity in Locating Geological Contact

Geoelectrical characterization for liquefaction at coastal zone in South Aceh

The Efficacy of Enhancing Horizontal Resolution (EHR) Technique in Shallow Subsurface Study Using 2D Resistivity Method at Bukit Bunuh

Interpretation of Subsurface Geological Structure of Massepe Geothermal Area Using Resistivity Data

Chapter 3. Geology & Tectonics

The Investigation of Groundwater Using Geophysical Techniques in Baitussalam Sub-district, Aceh Besar District, Aceh Province, Indonesia

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

GRAVITY AND MAGNETIC METHODS

Applying Seismic Refraction Method in Depicting Geological Contact at Bukit Bunuh, Lenggong, Perak, Malaysia

2D Resistivity Cross Section Interpretation of Shallow Hydrocarbon Reservoir in East Aceh, Indonesia

INTERGRATED GEOPHYSICAL METHODS USED TO SITE HIGH PRODUCER GEOTHERMAL WELLS

Vertical electrical sounding (VES) for subsurface geophysical investigation in Kanigiri area, Prakasam district, Andhra Pradesh, India

ELECTRICAL RESISTIVITY TOMOGRAPHY

The Role of Magnetotellurics in Geothermal Exploration

Overview of geophysical methods used in geophysical exploration

A New Combinational Terminology for Geothermal Systems

High Resolution Time-domain Induced Polarization Tomography with Merging Data Levels by Two Different Optimized Arrays for Slope Monitoring Study

Relevance of 2D Electrical Imaging in Subsurface Mapping: Case Study of National Animal Production Research Institute (NAPRI), Zaria.

Geophysics for Environmental and Geotechnical Applications

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

Visualizing Sungai Batu Ancient River, Lembah Bujang Archeology Site, Kedah Malaysia using 3-D Resistivity Imaging

Integrated Geophysical Interpretation On The Groundwater Aquifer (At The North Western Part of Sinai, Egypt)

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

Horizontal gradient and band-pass filter of aeromagnetic data image the subsurface structure; Example from Esh El Mellaha Area, Gulf of Suez, Egypt.

Determination of Geothermal Gradient in the Eastern Niger Delta Sedimentary Basin from Bottom Hole Temperatures

Ground Magnetic and Electrical Resistivity Mapping for Basement Structurs over Charnokitic Terrain in Ado-Ekiti Area, Southwestern Nigeria

Received 12 February 2010 Accepted 25 March Abstract. Keywords: sounding resistivity, gem-bearing gravel layer, Wenner, Schlumberger

Analysis of electrical resistivity data for the determination of aquifer depth at Sapele RD in Benin city

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

3D Electrical Resistivity Tomography (ERT) Survey of a Typical Basement Complex Terrain

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

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

APPLICATION OF ELECTRICAL RESISTIVITY TOMOGRAPHY FOR SAND UNDERWATER EXTRACTION

SPECTRAL ANALYSIS OF GROUND MAGNETIC DATA IN MAGADI AREA, SOUTHERN KENYA RIFT

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

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

1. are most likely to study the images sent back from Mars. A. Astronomers B. Geologists C. Doctors D. Engineers

Groundwater Assessment in Apapa Coast-Line Area of Lagos Using Electrical Resistivity Method

Coal Layer Identification using Electrical Resistivity Imaging Method in Sinjai Area South Sulawesi

Groundwater Sustainability at Wadi Al Bih Dam, Ras El Khaimah, United Arab Emirates (UAE) using Geophysical methods

Geophysics Course Introduction to DC Resistivity

Overview of Indonesian Geothermal System

A case study of crystalline limestone intrusion and fault zone identication using 2d eri technique in Ramco cements, pandalgudi mines, Tamilnadu

FINAL REPORT GEOPHYSICAL INVESTIGATION VILLAGE ALHAMBRA RETENTION POND SITE THE VILLAGES, FLORIDA

LIST OF FIGURES APPENDICES

Assessing our untapped energy resources. Derek Reay Geological Survey of Northern Ireland

Chapter 13. Groundwater

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

CONTENTS 1. INTRODUCTION. 2. THE D.C. RESISTIVITY METHOD 2.1 Equipment 2.2 Survey Procedure 2.3 Data Reduction

2-D Resistivity Imaging in Taiping, Perak, Malaysia

Airborne Geophysics to Map Groundwater. Bill Brown

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

Geology : Embedded Inquiry

Geothermal Exploration in Eritrea

Determination of Thickness of Aquifer with Vertical Electrical Sounding

Geophysical Surveys of The Geothermal System of The Lakes District Rift, Ethiopia

Geophysical Investigation of Ground Water Using Vertical Electrical Sounding and Seismic Refraction Methods

Wind Mountain Project Summary Memo Feeder Program

The Use of Vertical Electrical Sounding (VES) for Subsurface Geophysical Investigation around Bomo Area, Kaduna State, Nigeria

Rock Identification. invisible rhyolite andesite basalt komatiite. visible granite diorite gabbro peridotite

Scholars Research Library

Gravity Interpretation of the RUW (Rendingan-Ulubelu-Waypanas) Geothermal System in Tanggamus Regency, Lampung, Indonesia

MAGNETIC MODELING OF THE SUBSURFACE STRUCTURE OF SHIP ROCK, NEW MEXICO

Geochemical Study of Ampallas Geothermal Area, Mamuju District, West Sulawesi Province

RESISTIVITY IMAGING AND BOREHOLE INVESTIGATION OF THE BANTING AREA AQUIFER, SELANGOR, MALAYSIA. A.N. Ibrahim Z.Z.T. Harith M.N.M.

Engineering Geology. Igneous rocks. Hussien Al - deeky

RECHARGE ELEVATION OF HOT SPRING STUDY IN THE MT

Practice Test Rocks and Minerals. Name. Page 1

Conceptual model for non-volcanic geothermal resources - examples from Tohoku Japan

Geophysical technique applied to gemexploration in Chanthaburi, Eastern Thailand

Use of a Mise-a-la-Masse Survey to Determine New Production Targets in Sibayak Field, Indonesia

B7 Applications of DC resistivity exploration

RESISTIVITY OF ROCKS

DETECTION OF GROUNDWATER POLLUTION USING RESISTIVITY IMAGING AT SERI PETALING LANDFILL, MALAYSIA

Geophysical Investigation of the Precambrian Marble Occurrence in Itobe Area, Central Nigeria

3D Geometry of the Xade Complex inferred from Gravity and Magnetic Data

Complex geophysical and geological study of the Balaton Highlands: inferences from magnetic and geoelectric surveys

HIGH TEMPERATURE HYDROTHERMAL ALTERATION IN ACTIVE GEOTHERMAL SYSTEMS A CASE STUDY OF OLKARIA DOMES

transform boundary Photograph by Robert E. Wallace, USGS.

Subsurface Characterization using Electrical Resistivity(Dipole-Dipole) method at Lagos State University (LASU) Foundation School, Badagry

Soil Mechanics and Foundation Engineering, Vol. 53, No. 1, March, 2016 (Russian Original No. 1, January-February, 2016)

Japan Engineering Consultants, Inc., Energy and Industrial Technology Development Organization,Tokyo, Japan

1. What is the most important agent of chemical weathering on Earth? a. oxygen b. salt c. carbon dioxide d. carbonic acid e. water

Groundwater. (x 1000 km 3 /y) Oceans Cover >70% of Surface. Groundwater and the. Hydrologic Cycle

A Case Study on Ground Subsidence Using Ground Penetrating Radar

Geo-scientific Data Integration to Evaluate Geothermal Potential Using GIS (A Case for Korosi-Chepchuk Geothermal Prospects, Kenya)

2-D potential field modeling across the Hawtmi and Wiahatya faults in search of geothermal resources within the Umatilla Indian Reservation

DATA ACQUISITION METHODS FOR GROUNDWATER INVESTIGATION AND THE SITING OF WATER SUPPLY WELLS

APPLICATION OF GEOPHYSICS TO GEOTHERMAL ENERGY EXPLORATION AND MONITORING OF ITS EXPLOITATION

MICRO LEVEL GEO-RESISTIVITY SURVEY THROUGH V.E.S. TEST FOR GROUNDWATER FEASIBILITY STUDY AND SELECTION OF BORE WELL SITES IN PIPILI BLOCK OF PURI


Mapping of Manganese Ore Deposits by Using Geomagnetic Method in Aceh Jaya District, Nangro Aceh Darussalam Province, Indonesia

2D Resistivity Imaging Investigation of Solid Waste Landfill Sites in Ikhueniro Municipality, Ikpoba Okha Local Government Area,Edo State, Nigeria.

Does the SDCP need inputs from geology?

Integrated interpretation of multimodal geophysical data for exploration of geothermal resources Case study: Yamagawa geothermal field in Japan

Case Study: University of Connecticut (UConn) Landfill

Transcription:

Utilizing of Geophysical Method for Geothermal Exploration In Aceh Besar (Indonesia) Z.M Taqiuddin 1, M.M Nordiana 2, S Rosli 3 1Postgraduate Student, Geophysics Section, School of Physics, 11800 Universiti Sains Malaysia, Penang, Malaysia 2Senior Lecturer, Geophysics Section, School of Physics, 11800 Universiti Sains Malaysia, Penang, Malaysia 3 Assoc. Prof, Geophysics Section, School of Physics, 11800 Universiti Sains Malaysia, Penang, Malaysia ---------------------------------------------------------------------***--------------------------------------------------------------------- Abstract - Magnetic and 2-D resistivity imaging method were widely used in geothermal exploration. In this study geophysical method were conducted in two different locality situated in Ie-Seu 'Um and Iejue area. The purpose of this study to identify the geothermal environment and subsurface characteristics. Magnetic method was performed with proton precession magnetometer device and Global Positioning System (GPS) navigation for real-time measurements while 2- D resistivity imaging method were conducted with ABEM SAS4000 system with pole-dipole array. Two survey lines were conducted for each locality with 2 m electrode spacing for Ie- Seu 'Um area while 10 m electrode spacing for Iejue area. The magnetic data were processed using Surfer10 software for contouring to produce magnetic residual map while resistivity data were processed using Res2Dinv software. The magnetic residual map shows high magnetic value in geothermal area which is about >60 nt while resistivity give low values <7 Ωm. High magnetic intensity for geothermal area were caused due to increasing of iron oxide while low in resistivity are causes due to increasing the temperature as the depth increase. Low in resistivity <7 Ωm were suspected as hot water with composition of sulphur and clays. Key Words: Magnetic, 2-D resistivity imaging, Seulimeum fault 1.INTRODUCTION Geothermal described as heat contained within the earth. Geothermal system can be located in the region with normal or slightly above normal geothermal gradient and the region around plate margin where geothermal gradient higher than the average values. A volcanoes, hotsprings, and other thermal geological phenomena is a part from the geothermal system. The characteristic of geothermal system consist of the heat sources, reservoirs, a fluid which carries and transfer the heat, and recharge area [1]. The heat source can be either a very high temperature (>600 C) magmatic intrusion that reached shallow depths (5-10 km) or in certain low temperature as it increase with depth. The reservoir refer as a volume of hot permeable rocks from which the circulating fluids extract a heat and normally consists of hot fluids, vapors and gases depending on its temperature and pressure. The reservoir is surrounded by a cover of impermeable rocks and connected to recharge area where water (fluids) flows into the reservoir due the buoyancy forces. Geophysical technique such as resistivity and magnetic are widely successfully used in identifying and delineating the geothermal system. In resistivity exploration the parameter of interest is electrical resistivity of the rocks which reflects the properties of geothermal system. Resistivity is highly sensitive to temperature and geothermal alteration process and controlled by several important geothermal parameter such as temperature, fluid type and salinity, porosity, composition of the rocks, and also the presence of alteration minerals [2]. In a geothermal environment, the resistivity value decrease as a result increasing in temperature and salt concentration in geothermal reservoir [3] and almost volcanic rocks a magnetic due to presence of small amount of primary magnetic minerals (magnetite and titanomagnetite). Magnetic exploration in geothermal are important to identify the potential area of reducing magnetization due thermal activity [4]. Induced magnetization depends on the magnetic susceptibility of rocks and the magnitude of earth magnetizing field. The present of ferrimagnetic minerals such as magnetite and titanomagnetite in significant amount contributes a high values of magnetic susceptibility of volcanic rocks. The magnetic susceptibility in geothermal environment depends on the temperature and type of rocks in the environment it is in. Normally, in geothermal environment the heated rocks increase the bulk susceptibility due to growth of iron oxides ([5];[6]). Figure 1 shows the idealized model of geothermal system. Fig -1: Idealized model of geothermal system [7]. 2016, IRJET Impact Factor value: 4.45 ISO 9001:2008 Certified Journal Page 2234

2. General Geology Aceh province is located on the western Sumatra island with consist of four active major volcanoes mountains namely as Burni Telong mountain, Peuet Sagoe mountain, Jaboi mountain, and Seulawah Agam mountain which are more active compare from others. Generally, high temperature of 150-225 C dominated for this geothermal environment in Indonesia region [8]. Pre-Tertiary basements rocks outcrop mainly along the central spine of the Barisan Mountains, which extend the length of the island parallel to the southwest cost. The area from northeast and southwest is overlain by Tertiary sedimentary and volcanic rocks. Rock unit of all ages are transacted by Sumatra fault which follow the NW-SE trend [9]. The geology of Banda Aceh Quadrangle has been mapped (figure 2) [11]. The lithology of the study area, Ieju is dominated by Lam Tuba volcanic which is composed of andesitic to dacitic volcanic, pumiceous breccia, tuffs, agglomerate and ash flows with composition of tuffaceous and calcareous sandstones, conglomerates and mudstones [10]. The geological formation formed a topographic depression, occupied with alluvial flat and low flat-topped hills within Barisan Range ([11]; [12]). Ie Seu 'Um The apparent resistivity, ρ a are calculated based on equation 1.1. ρ a= k I v where; (1.1) V: voltage (V) I : current (A) or (ma) k : geometric factor depends on electrode arrangement In resistivity survey, the value k plays an important role, the depth of penetration of resistivity sounding depends on spacing of electrode which is refer to value of k. For this survey pole-dipole electrode array was chosen based on the deepest depth penetration beside provide a good resolution. Pole-dipole array is an asymmetrical arrangement which influence the anomalies by giving asymmetrical apparent resistivity in 2-D inversion model. To eliminate this asymmetry affect, the electrode are arrange in reverse manner as shown in figure 3 below. Forward measurement C 1 P 1 P 2 na a Iejue Reverse measurement P 1 P 2 a na C 1 Fig -3: Pole-dipole electrode array. Fig -2: Geological Map of Banda Aceh Quadrangle, Sumatra [11]. The resistivity method basically measures the resistivity distribution of the subsurface material. Rock types such as igneous and metamorphic typically have high resistivity values. Variations in electrical resistivity may indicate changes in composition, layer or contaminant levels [14]. The resistivity of these rocks is mainly dependent on the degree of fracturing. High fracturing of rocks will result a lower resistivity values. Table 1 shows the resistivity values of rocks and soil types [15]. 3. Theory 2-D resistivity imaging technique utilize direct current to investigate the electrical properties of the subsurface. The resistivity method basically measures the resistivity distribution of the overburden materials. Electrical resistivity in the earth depends on a combination of ohmic and dielectric effects related to lithology of subsurface [13]. 2016, IRJET Impact Factor value: 4.45 ISO 9001:2008 Certified Journal Page 2235

310 m 2.37 km International Research Journal of Engineering and Technology (IRJET) e-issn: 2395-0056 Table -1: Resistivity values of rocks and soil [15]. Material Resistivity (Ωm) Alluvium 10-800 Sand 60-1000 Clay 1-100 Groundwater (fresh) 10-100 Sandstone 8-4x10 3 Shale 20-2x10 3 Limestone 50-4x10 3 Granite 5x10 3-1x10 6 conducted near the hotspring area with ABEM SAS 4000 system. Pole-dipole array were selected for this survey with electrode spacing of 2 m. In Iejue area three magnetic survey lines were conducted with distance of 50 m distance interval for each magnetic station while two resistivity survey line are conducted with electrode spacing of 10 m (Figure 4). The magnetic were processed using Surfer10 software to produced full magnetic residual map. Resistivity data processing involve standard processing and model resistivity using Res2Dinv software. The data were then outputted into Surfer10 software for contouring and final presentation. Magnetic survey measure the magnitude and the orientation of the Earth's magnetic field intensity. The purpose of magnetic survey is to identify the subsurface geology on the basis of the anomalies in the earth's magnetic field due to magnetic properties of underlying rock [16]. When a magnet undergo magnetization it will have a remanent that refer to magnetization left behind after an external magnetic field is removed. It also refer as magnetic memory in magnetic storage and source of information on the past Earth's field in paleomagnetism [17]. Induced magnetization depends on the magnetic susceptibility of the material and magnitude and direction of ambient magnetic field while remanent magnetization reflects the previous past of magnetic history of the material [18]. Table 2 shows the general magnetic susceptibility value in common rocks and ores. The magnetic properties of highly magnetic rocks tend to be extremely variable and their magnetization is not directly proportional to the applied field [19]. A) B) 380 m 2.15 km Table -2: Magnetic susceptibility in common rocks and ores. Rock types Susceptibility (k) Altered ultra basics 10-4 to 10-2 Basalt 10-4 Gabbro 10-4 to 10-3 Granite 10-5 to 10-3 Andesite 10-4 Rhyolite 10-5 to 10-4 Metamorphic rocks 10-4 to 10-6 Most sedimentary rocks 10-6 to 10-5 Limestone and chert 10-6 Shale 10-5 to 10-4 4. METHODOLOGY Two survey area which is Ie-Seu 'Um and Iejue locality are choose for this study. Ie-Seu 'Um locality are famous established hotspring area used for eco-tourism purpose while Iejue area are located about 10 km from Ie-Seu 'Um hotspring. A magnetic survey was conducted surrounding the Ie-Seu 'Um hotspring area. Data was measured randomly by proton magnetometer with 20 m interval spacing for each station (rover) while 2-D resistivity survey line was Legend Magnetic cover area Resistivity line Fig -4: Study area; Ie-Seu 'Um (A) and Iejue (B). 2016, IRJET Impact Factor value: 4.45 ISO 9001:2008 Certified Journal Page 2236

5. RESULTS AND DISCUSSIONS Figure 5 shows the magnetic residual of Ie-Seu 'Um area. The magnetic value cover from -240 nt to 220 nt. High magnetic value from range of 0 nt to 220 nt assembles at western part while low magnetic value from range -180 nt to 0 nt dominated the eastern part of the map. The contrast zone of this magnetic map suspected a presence of fault zone. The black rectangle area shows the hotspring point with magnetic values covers from 60-120 nt. Line 2 b) Geothermal fluid Resistivity values (Ωm) Fig -6: 2-D inversion model resistivity line 1(a)and line 2(b). Hotspring area Figure 7 shows the magnetic residual map of Iejue area. The magnetic values covers from -150 nt to 550 nt. From this map a low magnetic values <-50 nt assembles at western part while high magnetic values >100 nt assemble at eastern part of the map. The highly contrast of magnetic zone may indicate the presence of fault zone. High magnetic value (>100 nt) are indicated as the presence of geothermal fluid which flow from northwest to southeast direction. Fault line Magnetic values (nt) Fig -5 : Magnetic residual map of Ie-Seu 'Um locality. Figure 6 (a and b) shows the 2-D inversion model of resistivity line 1 and 2 in Ie-Seu 'Um area. The resistivity values covers from 0-1800 Ωm with exploration depth up to 30 m. High resistivity values >200 Ωm may indicate the presence of volcanic tuff and rocks. A resistivity <50 Ωm indicate the saturated zone for this area. Low resistivity values <0-3.5 Ωm indicates the geothermal fluid consist of hot water with composition of sulphur and clays. Line 1 a) Geothermal fluid Magnetic values (nt) Fig -7: Magnetic residual map of Iejue area. Figure 8 shows 2-D inversion model resistivity of line 1 and line 2 in Iejue area. Low in resistivity value at range of 0-3.5 Ωm interpreted as geothermal fluid which consist of hot water and hot mud with composition of sulphur and clays. Low in resistivity value due to increase in temperature as it is increase with depth. From this result, two main zone are identify which are a top layer with thickness of 120 m at resistivity of >50 Ωm. The second layer are identify at resistivity value of 5-50 Ωm that interpreted as a cap layer for a geothermal system. 2016, IRJET Impact Factor value: 4.45 ISO 9001:2008 Certified Journal Page 2237

Line 1 a) Line 2 b) Resistivity values (Ωm) Fig -8: 2-D inversion model resistivity of line 1(a) and line 2(b) in Iejue area. 6. CONCLUSIONS Magnetic and 2-D resistivity imaging technique capable in delineating the geothermal system. High magnetic value >100 nt and low in resistivity value of 0-3.5 Ωm indicate the presence of geothermal fluid. High magnetic values in geothermal system are causes due to growth in iron oxide while low in resistivity value causes by increasing temperature as it increase in depth. ACKNOWLEDGEMENT The authors wish to thank to lecturers, technical staffs and all member of Geophysics group, Universiti Sains Malaysia (USM), Penang, Malaysia, all students and staffs of Faculty of Science, Universitas Syiah Kuala (UNISYIAH), Banda Aceh (Indonesia) for their assistance during the geophysical field data acquisition. REFERENCES Geothermal fluid Geothermal fluid [1] Zira, A. M. (2013). Challenges and Prospects of Geophysical Exploration of Geothermal System for National Development. Academic Journal of Interdisciplinary Studies, 2(12), 137 143. [2] Hersir G. P., Flóvenz, Ó. G., (2013). Resistivity surveying and electromagnetic methods. IGA Academy Report, 1 10. [3] Neawsuparp, K., Soisa, T., and Charusiri, P. (2010). Physical Characteristics of Pong Kum Hot Spring, Chiang Mai, Thailand, Using Ground Geophysical Investigation, (April), 25-29. [4] Georgsson, L.S. (2009). Geophysical methods in geothermal exploration. Short Course on Surface Exploration for Geothermal Resources, UNU-GTP and La Geo, El Salvador. [5] Dunlop, DJ. Thermal enhancement of magnetic susceptibilty, J. Geophys., 40, 439-451, 1974 [6] Thompson, R. and Oldfield, F. Environmental Magnetism, 227 pp., Allen & Unwin, Boston, 1986. [7] Berktold, A. (1983) Electromagnetic studies in geothermal regions: Geophys. Surv., 6, 173-200. [8] Syukri, M. (2014). The Investigation of Hot Spring Flow Using Resistivity Method at Geothermal Field Ie-Seu Um, Aceh-Indonesia. Electronic Journal of Geotechnical Engineering, Vol.19 (Bund.K),n2419 2427. [9] Barber, a. J., and Crow, M. J. (2003). An Evaluation of Plate Tectonic Models for the Development of Sumatra. Gondwana Research, 6(1), 1 28. [10] Bennett, J. D., McC, D., Bridge, Cameron, N. R., Djunuddin, A., Ghazali, S. A., Jeffrey, D.H., Kartawa, W., Keats, W., Rock, N. M. S., Thomson, S. J. and Whandoyo, R. (1981) Geologic map of the Banda Aceh Quadrangle, Sumatra. [11] Katili, J., and Hehuwet, F. (1967), On the occurrence of large transcurrent faults on Sumatra Indonesia: Osaka, Japan, Osaka City University Journal of Geoscience, Vol. 10, pp 1-17. [12] Page, B. G. N., Bennett, J. D., Cameron, N.R., Bridge, D., McC Jeffery D. II., Keats, W. and Thaib, J. (1979), A review of the main structural and magmatic features of northern Sumatra, Journal of Geological Society of London, Vol. 1, 36, pp 569-579. [13] Edwin, S. Robinson, Cahit, Coruh. (1988). Basic Exploration Geophysics. [14] Loke, M. H. (2000). Electrical imaging surveys for environmental and engineering studies. A practical guide to 2-D and 3-D surveys, 6574525, 59. [15] Keller G.V. and Frischknecht F.C. (1966). Electrical methods in geophysical prospecting, Pergamon Press Inc., Oxford. [16] Keller, G. V. (1986). An Introduction to Geophysical Exploration. Eos, Transactions American Geophysical Union, Vol. 67 [17] Nurgalieva, N. D., & Yassonov, P. G. (2013). Magnetic Susceptibility, Remenent Magnetization and Coercivity Varitions Along Soil Profile, 8(5), 357 361. [18] William J. Hinze, Ralph R. B. von Frese and Afif H. Saad (2013), Gravity and MagneticnExploration (Principle, Practice, and Applications), Cambridge University Press, New York. [19] Milsom. J., (2003). Field Geophysics, 3rd Edition, University College London. 2016, IRJET Impact Factor value: 4.45 ISO 9001:2008 Certified Journal Page 2238