A Poisson s Ratio Distribution From Wadati Diagram as Indicator of Fracturing of Lahendong Geothermal Field, North Sulawesi, Indonesia

Similar documents
DETAILED IMAGE OF FRACTURES ACTIVATED BY A FLUID INJECTION IN A PRODUCING INDONESIAN GEOTHERMAL FIELD

Hijiori HDR Reservoir Evaluation by Micro-Earthquake Observation

Iwan Yandika Sihotang, Tommy Hendriansyah, Nanang Dwi Ardi

Microearthquake (MEQ) Investigation Reveals the Sumatran Fault System in Hululais Geothermal Field, Bengkulu, Indonesia

The Coso Geothermal Area: A Laboratory for Advanced MEQ Studies for Geothermal Monitoring

Earthquakes.

DOWN-HOLE SEISMIC SURVEY AND VERTICAL ELECTRIC SOUNDINGS RABASKA PROJECT, LÉVIS, QUÉBEC. Presented to :

Earthquakes and Earth s Interior

Tu P8 08 Modified Anisotropic Walton Model for Consolidated Siliciclastic Rocks: Case Study of Velocity Anisotropy Modelling in a Barents Sea Well

NUMERICAL MODELING STUDY OF SIBAYAK GEOTHERMAL RESERVOIR, NORTH SUMATRA, INDONESIA

Characterization of Induced Seismicity in a Petroleum Reservoir: A Case Study

Geothermometer, Geoindicator and Isotope Monitoring in Lahendong Wells during

What is an Earthquake?

Module 7: Plate Tectonics and Earth's Structure Topic 4 Content : Earthquakes Presentation Notes. Earthquakes

ERTH2020 Introduction to Geophysics The Seismic Method. 1. Basic Concepts in Seismology. 1.1 Seismic Wave Types

Chapter 6. Conclusions. 6.1 Conclusions and perspectives

they help tell Earth s story! Regents Earth Science With Ms. Connery

Global geophysics and wave propagation

GeothermEx, Inc. GEOTHERMAL RESERVOIR ASSESSMENT METHODOLOGY FOR THE SCIENTIFIC OBSERVATION HOLE PROGRAM, KILAUEA EAST RIFT ZONE, HAWAII TASK 1 REPORT

Earthquakes. Forces Within Eartth. Faults form when the forces acting on rock exceed the rock s strength.

Dangerous tsunami threat off U.S. West Coast

Determining the Earthquake Epicenter: Japan

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

Microearthquake Focal Mechanisms

Internal Layers of the Earth

Overview of Microseismic in Sustainable Monitoring of Geothermal Reservoirs in Indonesia*

OVERVIEW OF THE WAIRAKEI-TAUHARA SUBSIDENCE INVESTIGATION PROGRAM

of other regional earthquakes (e.g. Zoback and Zoback, 1980). I also want to find out

Numerical Simulation Study of the Mori Geothermal Field, Japan

NOTICE CONCERNING COPYRIGHT RESTRICTIONS

Geophysical Methods in Geothermal Exploration

Empirical Green s Function Analysis of the Wells, Nevada, Earthquake Source

Using Repeat Microgravity Measurements to Track Reinjection in Liquid-Dominated Fields

Extending the magnitude range of seismic reservoir monitoring by Utilizing Hybrid Surface Downhole Seismic Networks

TIME-DEPENDENT SEISMIC TOMOGRAPHY OF GEOTHERMAL SYSTEMS

Earthquakes. Pt Reyes Station 1906

Topic 5: The Dynamic Crust (workbook p ) Evidence that Earth s crust has shifted and changed in both the past and the present is shown by:

UNIT - 7 EARTHQUAKES

Synthetic Seismogram A Tool to Calibrate PP & PS Seismic Data

Section 19.1: Forces Within Earth Section 19.2: Seismic Waves and Earth s Interior Section 19.3: Measuring and Locating.

5. What is an earthquake 6. Indicate the approximate radius of the earth, inner core, and outer core.

GEOPHYSICAL PROSPECTING: DYNAMIC RESERVOIR CHARACTERIZATION AND TIME-LAPSE MULTICOMPONENT SEISMOLOGY FOR RESERVOIR MONITORING UNESCO EOLSS

The Role of Magnetotellurics in Geothermal Exploration

Chapt pt 15 er EARTHQUAKES! BFRB P 215 ages -226

Lithology prediction and fluid discrimination in Block A6 offshore Myanmar

Modeling pressure response into a fractured zone of Precambrian basement to understand deep induced-earthquake hypocenters from shallow injection

Earthquakes Chapter 19

WAMUNYU EDWARD MUREITHI I13/2358/2007

EVALUATING HEAT FLOW AS A TOOL FOR ASSESSING GEOTHERMAL RESOURCES

Renewability Assessment of the Reykjanes Geothermal System, SW-Iceland. Gudni Axelsson et al. (see next slide) Iceland GeoSurvey (ÍSOR)

Section Forces Within Earth. 8 th Grade Earth & Space Science - Class Notes

Earthquakes and Earthquake Hazards Earth - Chapter 11 Stan Hatfield Southwestern Illinois College

Seismicity and the SWD-C4A well: An ongoing UIC case study in the Denver Basin, Colorado

Earthquake patterns in the Flinders Ranges - Temporary network , preliminary results

Earthquakes. Building Earth s Surface, Part 2. Science 330 Summer What is an earthquake?

Chapter 4 Earthquakes and Tsunamis

Chapter 4 Earthquakes and Tsunamis. Geology of the Hawaiian Islands. Any Questions? Class March Mean = 78.

A look into Gassmann s Equation

NOTICE CONCERNING COPYRIGHT RESTRICTIONS

C) 10:20:40 A) the difference between the arrival times of the P -wave and the S -wave

A VOLCANOSEISMIC APPROACH TO GEOTHERMAL EXPLORATION AND RESERVOIR MONITORING: OLKARIA, KENYA AND CASA DIABLO, USA

Dallas-Fort Worth and Cleburne, TX

SEISMIC MONITORING OF EGS STIMULATION TESTS AT THE COSO GEOTHERMAL FIELD, CALIFORNIA, USING MICROEARTHQUAKE LOCATIONS AND MOMENT TENSORS

LAB 6: Earthquakes & Faults

Dynamic Triggering Semi-Volcanic Tremor in Japanese Volcanic Region by The 2016 Mw 7.0 Kumamoto Earthquake

Earthquake Lab Introduction: Purpose: Method: Equipment needed: Fundamental Concepts and Key Terms:

THREE-DIMENSIONAL COMPRESSIONAL- AND SHEAR-WAVE SEISMIC VELOCITY MODELS FOR THE SOUTHEAST GEYSERS GEOTHERMAL FIELD, CALIFORNIA

ES Ch 19 Earthquakes 1

National Geophysical Survey Science Opportunities

Geothermal Model of the Lahendong Geothermal Field, Indonesia

Location uncertainty for a microearhquake cluster

Kinematic inversion of pre-existing faults by wastewater injection-related induced seismicity: the Val d Agri oil field case study (Italy)

SHALE GAS AND HYDRAULIC FRACTURING

Geophysical model response in a shale gas

Activity Submitted by Tim Schroeder, Bennington College,

Hypocenter Estimation of Detected Event near Venlo on September 3rd 2018

Some aspects of seismic tomography

Probabilistic Inversion Technique for Seismic Data

GEOTHERMAL POTENTIAL OF ST. KITTS AND NEVIS ISLANDS

The Initial-State Geochemistry as a Baseline for Geochemical Monitoring at Ulubelu Geothermal Field, Indonesia

An empirical method for estimation of anisotropic parameters in clastic rocks

Calpine Corporation. NW Geysers Enhanced Geothermal System Demonstration Project. Community Update Calpine Geysers Visitor Center Middletown, CA

Heat (& Mass) Transfer. conceptual models of heat transfer. large scale controls on fluid movement. distribution of vapor-saturated conditions

Geophysical Site Investigation (Seismic methods) Amit Prashant Indian Institute of Technology Gandhinagar

Estimation of fracture porosity and permeability using radon as a tracer

SEG/New Orleans 2006 Annual Meeting

Locating Events using Borehole Microseismic Monitoring by Inclusion of Particle Motion Analysis

Earthquakes Earth, 9th edition, Chapter 11 Key Concepts What is an earthquake? Earthquake focus and epicenter What is an earthquake?

Seismicity Associated with Geothermal Systems

The Earthquake of Padang, Sumatra of 30 September 2009 scientific information and update

Lecture 9 faults, folds and mountain building

RESERVOIR MODELING OF GEOTHERMAL ENERGY PRODUCTION FROM STRATIGRAPHIC RESERVOIRS IN THE GREAT BASIN

Lecture 4: Earthquakes and Seismic Waves

Magnetic and Gravity Methods for Geothermal Exploration

Possible Fluid Driven Open Crack Events Observed in TCDP Borehole Seismometers

Data Repository: Seismic and Geodetic Evidence For Extensive, Long-Lived Fault Damage Zones

2, from which K = # " 2 % 4 3 $ ) ( )

Faults, folds and mountain building

Earthquakes. Photo credit: USGS

Transcription:

Proceedings World Geothermal Congress Antalya, Turkey, - April A Poisson s Ratio Distribution From Wadati Diagram as Indicator of Fracturing of Lahendong Geothermal Field, North Sulawesi, Indonesia Silitonga T. H., Siahaan E. E., Suroso Mailing address, PT. Geodipa Energi, Jl. Karawaitan No., Bandung 6, Indonesia E-mail address, thsilitonga@geodipa.co.id Keywords: Poisson s ratio, Wadati diagram, fracturing, permeability. ABSTRACT Poisson s ratio from Wadati Diagram using micro earthquake monitoring is relatively simple and useful for direct indicator of fracture. A short term micro earthquake monitoring in the Lahendong geothermal field shows two distinct shallow and deep hypocenters distributions. This paper describes the Poisson s ratio s distribution from Wadati method in the Lahendong geothermal field. The higher Poisson s ratio associated with shallow hypocenters indicated by the fractured rocks, and the lower Poisson s ratio associated with deeper hypocenters and un-fractured rocks with silicified rocks.. INTRODUCTION The Lahendong geothermal field situated in the North Sulawesi Province which is approximately km from the capital city of Manado. One of the geophysical method used to detect reservoir permeability is micro earth quake. Each event of micro earthquake can be associated with permeable zone. As a result, this geophysical method is very helpful to develop a geothermal field. Poisson s ratio is one of the important parameters that can be used as a fracture indicator through simple calculation of Wadati diagram. Whereas, well data like cores, minerals and loss circulation can provide reservoir permeability are very expensive and limited. This paper describes the calculation and interpretation of Poisson s ratio derived from the short monitoring of micro earthquake at the beginning of production in the Lahendong geothermal field. More over, one layer modeling of earthquake distribution and Poisson s ratio calculation by using Wadati diagram indicate different fracture intensity between shallow and deep reservoir.. MICRO-EARTHQUAKE DISTRIBUTION The Lahendong geothermal field is confined to an area about km as shown Figure. Five MEQ s stations with spacing less than km were utilized to cover interest area. Two seismometers L-C (-components) and three seismometers SS- ranges (-component) with natural frequency Hz were connected to DATAMARK LS- 8SH. By using sample interval of msec, gains times, pre trigger data, and low pass filter Hz were applied to detect micro-earthquake with magnitude < Richter s scale focal depths < km. The Lahendong geothermal field was produce steam about 6 tonnes/hr and about 7 tonnes/hr brine were injected back into the earth through an injection well LHD- at WHP up to -.8 bar. It is believed that local changes in the stress field are obviously caused by production and reinjection; however; the induced seismicity is not confined to the immediate regions of production and re-injection, but also occurs at greater depths (Figure ). DATAMARK LS-8SH is a high precision data logger capable of recording induced or tectonic earthquakes data. Since it is synchronize with UTC by GPS s time calibration function, each logger can record the data independently. It has been claimed by the manufacture that the time accuracy is ± µsec. In the Lahendong geothermal field, the typical micro earthquakes, see Figure, are as follow: duration < second, the depths < km, and the time differences between P and S arrival (S-P) < second. Therefore, a simple assumption that model is a horizontal layered velocity structure, assuming constant, isotropic velocity within one layer, is reasonable enough to get a quick interpretation as preliminary result. A small error in picking of the arrival time and the inaccuracies of velocity model can produce a large error in determining the hypocenter and epicenter location. Therefore, both factors must be controlled, carefully. Firstly, as we know, it is easier to pick P wave than S wave, so P wave data from stations is much better than P and S wave data from stations. Using this procedure, which controlled by data from the shortest station usually associated with good event, it is expected to minimize the error. Secondly, sensitivity analysis was used to obtain the representative velocity of the area. Since there is no direct velocity measurement, and no velocity data from petrophysics, the whole micro-earthquakes data were calculated using a simple linier inversion program (Gede Swantika) with various velocities ranges from. to km /sec. Most of data converged quickly to provide satisfactory hypocenter locations using a compressional wave velocity of. km/sec. Since full production began, as well as during commissioning, four production wells in the Southern block were used to produce MWe. The brine from a production separator were injected into the ground by using one injection well, also in the Eastern block. Although the mechanism which causes induced seismicity is not understood yet, but it has been studied that injection of fluid into the ground already deformed by regional stress can cause man made earthquakes. Another way, detailed study of the Geyser reservoir showed that seismicity can also increase due to production (Denlinger et al, 8)

6 8 Legend : : <M< : <M< : <M< : Deeper :Shallower LHD- : Proposed development area : Interpreted Natural recharge : Sta. Numb. & Seis. -component : Sta. Numb. & Seis. -component : Production Well : Injection Well : Shutin Well] 8 LHD- THE LAHENDONG PROSPECT AREA 8 Southern block 7 barrier? LHD- E Eastern block 6 LHD-7 LHD- 7 8 6 6 Recharge F8 7 7 8 7 6 6 8 LHD- 6 LHD- 8 LHD-6 North 6 SW km Figure The map of the Lahendong geothermal field showing the prospect area, production and injection wells, microearthquakes distributions in the period Desember - Februari (Modified from Sudarman S.) Z (km) - - SW Mt.Lengkoan Legend : : Production Well : Injection Well : Microearthquake & No. LHD- LHD- A=km - 6 NE W LHD- 6 Shallow Resrvoir 7 6 Propilitik 6 Deeper Reservoir LHD-7 7 Mt.Tampuso < C 8 6 6 Natural recharge (Interpreted) E 8 km Figure The cross section of the Lahendong geothermal field showing vertical distribution of micro-earthquake and injection and production wells (Modified from Sudarman S.).

.6.. Fig. b ST-NS -. -. -.6 :: ::7 :: ::7 :: ::7.6.. Fig. a ST-UD -. -. -.6 :: ::7 :: ::7 :: ::7.6.. Fig. c ST-EW -. -. -.6 :: ::7 :: ::7 :: ::7 West Indonesian Time (Hr:mnt:sec) Figure A Typical of micro-earthquake record from a -component seismometer in the Lahendong geothermal field, during monitoring (Fig. a, record from a vertical component, Fig. b and c record from horizontal components). WADATI DIAGRAM Wadati diagram is important in processing microearthquake data for two reasons. The first is that we can control the correct picks of travel time of P wave and S wave (ts and tp) from record data. The second reason is that we can compare the origin time of earth-quake (To, if the time of S and P are equal, it means the rock are fractured) with To from a computer program e.g. linier inversion. In many cases, it is hard to pick S wave accurately even though we use -component seismometers. Particularly, in the Lahendong geothermal field we only use seismometers -components. Therefore, the use of Wadati diagram is very important to control the correct picks. Table The corrects picks of micro-earthquacke data on st January Sta. To (Wadati) To (Inversion) Error Poisson s ratio Tp..7.7.7 Ts.....7.7.6.6 Ts Tp..8.77.6 After picking all of the ts and tp from all micro-earthquake stations as shown in Tabel-, we can construct a Wadati diagram by plotting all Ts Tp vs Tp for all the receiver stations, as shown in Figure.

Micro-earthquake date: st Jan time: :8:.7 The graph of Poisson's ratio versus depth Shallow reservoir Ts-Tp (Second) y =,8x -,87 R =, depth (Km) - - - deeper reservoir,,,,, Tp (Second) -,,,,, Poisson's Ratio Figure-: Wadati diagram of micro-earthquakes on st January. POISSON S RATIO CONCEPT Basically Poisson s ratio is defined by the ratio of the strain perpendicular to either type of deforming force to that in the direction of the force itself. In other way, by using compressional velocity and transversal velocity or the comparison of both velocity from Wadati diagram the Poisson s ratio is defined as follow : σ = ( Vp ( Vp Vs) Vs ) where σ, Vp dan Vs are Poisson s ratio, the compressional wave velocity (longitudinal) and the shear wave velocity (transversal) respectively. If we assume that the ray-paths are confined to homogeneous geology, experimental and theoretical results indicate that additional fracturing of a fluid-filled rock will cause an increase in the Poisson s ratio. The P-wave velocity decreases slightly and the S-wave velocityis significantly reduced. This phenomena is observed in a survey on Hawaii near the successful HGO well (Butler David). The study of refraction survey indicates the anomaly of the velocity structure in the high temperature; the velocity of the compression wave (Vp) reduced by the fluid composition and temperature; velocity of transversal wave also influenced by temperature but not affected by fluid composition. (Ito vide Soengkono ). In the steam dominated system, Vp decreases faster than Vs. The two parameters provide poisson ratio lower than.. In the water dominated system at the higher porosity pressure, the Poisson s ratio is bigger than.. The Vs decreases faster than Vp. Figure : The graph of Poisson s ratio value of Lahendong geothermal field versus depth relative to sea level. consolidated Soft, poorly hard, rigid ρ..... Lahendong shallow reservoir Lahendong deeper reservoir (Andesit) Geyser Wairakei (greywacke) Most crustal rocks Coso, California Fracturing Ciero Prieto,Mexico Steam dominated Hot water dominated Vp/Vs Figure 6 Poisson s ratio of some geothermal fields (from Soengkono, S.; Dobrin) over the world wide versus Vp/Vs. Figure 6 shows the significant differences between Poisson s ratio in the shallower zone and deeper zone. The best explanation of the Poisson s ratio anomaly is where the rocks type is different from in the anomalous area.. POISSON S RATIO OF THE LAHENDONG GEOTHERMAL FIELD. The Poisson s ratio of the Lahendong geothermal field indicate a significant differences between shallower and deeper zone. In the shallower zone, the Poisson s ratio varies between. up to.. In the deeper zone, the Poisson s ratio varies from.7 to..

6. CORRELATION OF POISSON RATIO WITH WELL BORE DATA The value of Poisson s ratio with depth in the Lahendong geothermal field can be divided into two zones which correspond to shallower and deeper zone. The shallower zone with interval of Poisson s ratio varies from. to. has depth from surface to meter depth and the deeper zone with interval of Poisson s ratio varies from.7 to. has depth is deeper than meter. The correlation between cores and Poisson s ratio data provide the two reservoir zones. The shallower zone is distinguished by the intensively fractured rock, the loss circulation, partial loss circulation and drilling break. All the steam in Lahendong field derived from this zone which is the main target to develop another X MWe. In the deeper zone, the rock is silicified without loss circulation zones, see Figure 7. depth: (m) 8 POISSON'S RATIO - MEQ Shallow Reservoir Well LHD- Well LHD- Well LHD- (Eastern side) (Western side) (Southern side) TD= m TD= 7m TD= 87m steam kick 6 lpm lpm drilling break drilling break. Higher Poisson s ratio is associated with deeper hypocenters characterized by un-fractured rocks and the silicified alteration. ACKNOWLEDGMENTS The authors would like to thank Pertamina Geothermal Management for permission to publish this paper. REFERENCES Hochstein M. P : Geothermal Exploration-66.6, Geothemal Institute, The University of Auckland,. Butler David : Notes on Seismic Exploration for Geothermal Energy, MicroGeophysics Corporation, W. th Aveneu, Wheat Ridge, Colorado 8, 8. Dobrin M. : Introduction To Geophysical Prospecting,. Denlinger, R. P.; Isherwood, W. F.; Kovach, R. L : A Geodetic Analysis of reservoir Depletion at The Geysers Steam Field (Denlinger et al, 8) - - - - - Deep Reservoir - lpm - Silicification - No loss circulation - lpm - Tight permeability Note: = Total Loss Circulation = Partial Loss Circulation Gede Swastika, Pembuatan Program Inversi Linier Untuk Menentukan Hypocenter Gempa, a Report in Master of Science Program, Bandung Institute of Technology,. Siahaan E.E. : Hydrothermal Mineralogy and veining system in well LHD-, LHD-, LHD-, -.......... poissons vs fracture Siahaan E.E. : Interpretation of Land sat and Aerial Photo of Lahendong (internal report Geothermal Pertaminaunpublished),. Figure 7 Poisson s ratio versus fracturing system (Total Loss Circulation) and lithology (Modified from Siahaan, E. E.) 7. CONCLUSIONS. Lower Poisson s ratio is associated with shallower hypocenters, indicative of fractured rock s. Soengkono, S:.Geothermal Geophysics Passive Seismic Method (66.6) Geothemal Institute, The University of Auckland,. Sudarman S. : Review of Lahendong Model and Development Strategy, (Internal Report, Geothermal Pertamina-unpublished)