INVESTIGATION OF BURSA, ESKIKARAAGAC USING VERTICAL ELECTRICAL SOUNDING METHOD

Similar documents
P62 In-mine Geoelectric Methods for Detection of Tectonic Disturbances of Coal Seams

A027 Uncertainties in Local Anisotropy Estimation from Multi-offset VSP Data

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS

Soil Investigation for Determining Aquifer Position at Pasir Impun Area, Bandung, West Java Using Resistivity Method Wenner Configuration

Predict Global Earth Temperature using Linier Regression

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

Review of Calculus, cont d

Math& 152 Section Integration by Parts

Motion of Electrons in Electric and Magnetic Fields & Measurement of the Charge to Mass Ratio of Electrons

5.7 Improper Integrals

Applications of Bernoulli s theorem. Lecture - 7

Math 1B, lecture 4: Error bounds for numerical methods

Measuring Electron Work Function in Metal

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

CHAPTER 4a. ROOTS OF EQUATIONS

POLYPHASE CIRCUITS. Introduction:

Week 10: Line Integrals

UNIT 1 FUNCTIONS AND THEIR INVERSES Lesson 1.4: Logarithmic Functions as Inverses Instruction

The Velocity Factor of an Insulated Two-Wire Transmission Line

4.4 Areas, Integrals and Antiderivatives

New Expansion and Infinite Series

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

APPROXIMATE INTEGRATION

Foundation Pit Leakage Detection Method Based on Electrical Resistivity Imaging Technology

SESSION 2 Exponential and Logarithmic Functions. Math 30-1 R 3. (Revisit, Review and Revive)

KINEMATICS OF RIGID BODIES

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

The Properties of Stars

ragsdale (zdr82) HW2 ditmire (58335) 1

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

Estimation of the particle concentration in hydraulic liquid by the in-line automatic particle counter based on the CMOS image sensor

MATHS NOTES. SUBJECT: Maths LEVEL: Higher TEACHER: Aidan Roantree. The Institute of Education Topics Covered: Powers and Logs

CHM Physical Chemistry I Chapter 1 - Supplementary Material

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

Chapter 0. What is the Lebesgue integral about?

Space Curves. Recall the parametric equations of a curve in xy-plane and compare them with parametric equations of a curve in space.

New data structures to reduce data size and search time

14.4. Lengths of curves and surfaces of revolution. Introduction. Prerequisites. Learning Outcomes

CONTRIBUTION TO THE EXTENDED DYNAMIC PLANE SOURCE METHOD

Synoptic Meteorology I: Finite Differences September Partial Derivatives (or, Why Do We Care About Finite Differences?

Operations with Polynomials

Prof. Anchordoqui. Problems set # 4 Physics 169 March 3, 2015

Name Solutions to Test 3 November 8, 2017

Physics 202H - Introductory Quantum Physics I Homework #08 - Solutions Fall 2004 Due 5:01 PM, Monday 2004/11/15

7.2 The Definite Integral

Design Data 1M. Highway Live Loads on Concrete Pipe

Math 8 Winter 2015 Applications of Integration

Definition of Continuity: The function f(x) is continuous at x = a if f(a) exists and lim

Precalculus Spring 2017

Estimation of Global Solar Radiation at Onitsha with Regression Analysis and Artificial Neural Network Models

Available online Journal of Scientific and Engineering Research, 2017, 4(7): Research Article

7.6 The Use of Definite Integrals in Physics and Engineering

Recitation 3: More Applications of the Derivative

Riemann is the Mann! (But Lebesgue may besgue to differ.)

D i k B(n;D)= X LX ;d= l= A (il) d n d d+ (A ) (kl)( + )B(n;D); where c B(:::;n c;:::) = B(:::;n c ;:::), in prticulr c n = n c. Using the reltions [n

Algebra Readiness PLACEMENT 1 Fraction Basics 2 Percent Basics 3. Algebra Basics 9. CRS Algebra 1

The steps of the hypothesis test

Finite Element Determination of Critical Zones in Composite Structures

approaches as n becomes larger and larger. Since e > 1, the graph of the natural exponential function is as below

Designing Information Devices and Systems I Spring 2018 Homework 7

Lecture 20: Numerical Integration III

Abstract inner product spaces

Non-Linear & Logistic Regression

NOT TO SCALE. We can make use of the small angle approximations: if θ á 1 (and is expressed in RADIANS), then

Vorticity. curvature: shear: fluid elements moving in a straight line but at different speeds. t 1 t 2. ATM60, Shu-Hua Chen

PART 1 MULTIPLE CHOICE Circle the appropriate response to each of the questions below. Each question has a value of 1 point.

Intro to Nuclear and Particle Physics (5110)

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.

3.1 Review of Sine, Cosine and Tangent for Right Angles

Jackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

C3D. Conductivity In Less Than Three Dimensions

pivot F 2 F 3 F 1 AP Physics 1 Practice Exam #3 (2/11/16)

19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007 A NON-CONTACT SYSTEM FOR TRANSPORTING OBJECTS USING ULTRASONIC LEVITATION

WELL SCHEDULE. Longitude -84 0' 24.10'' Ground Elevation Casing Elevation OWNER WELL CONSTRUCTION DESCRIPTION. Name of Aquifer: TYPE OF DRILLING

Purpose of the experiment

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

Chapter 6 Electrostatic Boundary Value Problems. Dr. Talal Skaik

Heat flux and total heat

7.1 Integral as Net Change and 7.2 Areas in the Plane Calculus

MAT187H1F Lec0101 Burbulla

fractions Let s Learn to

Section 4.8. D v(t j 1 ) t. (4.8.1) j=1

Math 113 Exam 2 Practice

(See Notes on Spontaneous Emission)

4 7x =250; 5 3x =500; Read section 3.3, 3.4 Announcements: Bell Ringer: Use your calculator to solve

Numerical Integration

Tangent Line and Tangent Plane Approximations of Definite Integral

a < a+ x < a+2 x < < a+n x = b, n A i n f(x i ) x. i=1 i=1

5.2 Volumes: Disks and Washers

Supplementary Information for Directional Reflective Surface Formed via Gradient- Impeding Acoustic Meta-surfaces

Chapter 14. Matrix Representations of Linear Transformations

Pre-Session Review. Part 1: Basic Algebra; Linear Functions and Graphs

n f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1

Practical exercise 7. Surge tank

Physics 1402: Lecture 7 Today s Agenda

Minimum Energy State of Plasmas with an Internal Transport Barrier

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

SUPPLEMENTARY INFORMATION

A-Level Mathematics Transition Task (compulsory for all maths students and all further maths student)

Transcription:

INVESTIGATION OF BURSA, ESKIKARAAGAC USING VERTICAL ELECTRICAL SOUNDING METHOD Gökçen ERYILMAZ TÜRKKAN, Serdr KORKMAZ Uludg University, Civil Engineering Deprtment, Burs, Turkey geryilmz@uludg.edu.tr, skorkmz06@gmil.com Abstrct: The im of this reserch is to study the ground resistivity of region locted t Burs, Eskikrğç villge by using verticl electricl resistivity (VES) method which is one of the geoelectricl methods. For tht purpose, two Wenner rrys were pplied on the ground. The mesured VES dt is nlyzed vi softwre nmed RES2DINV. The progrm gives lyer thicknesses nd rel resistivity vlues on two dimensionl underground resistivity sections. Interprettions on the possible loctions of wter cn lso be mde. These interprettions cn lter be verified by drilling groundwter wells. Keywords: Verticl Electricl Sounding, RES2DINV, Wenner rry Introduction Tody, climte chnge leds to shortge of wter while the demnd from limited resources increses with popultion. Drilling is unknowingly used for the purpose of serching nd finding wter. Wheres, geoelectricl methods provides some informtion bout quifer lyers nd wter existence before drilling. Geoelectricl methods re the most common methods in explortion geophysics. The method, which mkes use of the resistivity feture of erth, is used in serch for underground ssets such s wter, geotherml energy, minerls nd petrol. In this study, geophysics resistivity mesurements hve been conducted on previously chosen field over the quifer t Eskikrğç villge. Resistivity mesurement ws tken using Wenner rrys. Two pieces of Wenner rrys were mde in totl. A progrm clled RES2DINV (Loke, 2000) ws used in the evlution of these mesurements. Consequently, rel resistivity vlues nd thicknesses of geologicl lyers were obtined nd mpped in two dimensions. By considering rel resistivity vlues nd the two-dimensionl mp, inference bout wter existence ws mde. Geoelectricl Methods Bsiclly, geoelectricl methods re implemented in three different wys. These re Self Potentil (SP), Induced Polriztion (IP) nd Verticl Electricl Sounding (VES). Verticl electricl sounding method is bsed on the principle of mesuring the potentil difference formed by electric current delivered to the ground by two current electrodes. The im of this method is to determine the depth nd resistivity of underground lyers by using surfce potentil mesurements. VES method hs different geometricl rrys depending on the problem on the ground. The most common used rrys re Schlumberger, Wenner nd Dipol-Dipol. Bsic differences between rrys re investigtion depth nd horizontl solution ccurcy. Compring Wenner nd Schlumberger; both electrode rrys hve the sme current trnsmission horizontlly, however, investigtion depth is more in Schlumberger method. For the sme opening width, ccurcy of mesurement is much better in Wenner rry. According to these comprisons; generlly, Wenner is preferred for shllow studies. As to Schlumberger, it cn be used both shllow nd deep reserch (Bşokur, 2004). In Figure 1, A nd B re current, M nd N re potentil electrodes in Wenner rry. In prctice, electrode current is emitted from two points like A nd B nd then potentil difference is mesured from seprte two points like M nd N. www.tojst.net Copyright The Online Journl of Science nd Technology 69

A M N B K=2 π Figure 1. Wenner rry After clculting the geometric fctor (K), pprent resistivity is clculted with Eq. 1. ρ = V I K (1) ρ= Apprent resistivity (ohm-m) V= Voltge (mv) I= Current (ma) K= Geometric fctor Essentilly, erth consists of complex structures rther thn isotropic semi-infinite medi. Becuse of this reson, resistivity reltionships do not give rel ground resistivity vlues. Resistivity vlue obtined with current given into the ground nd obtined using severl equtions is nmed s pprent resistivity. Apprent resistivity is equivlent to rel resistivity in cse of isotropic semi-infinite medi (Bşokur, 2004). To obtin rel resistivity of lyers nd thickness of lyers grphicl methods re used. Generlly, on logrithmic pper, pprent resistivity vlue is plced on verticl xis nd electrode spcing on horizontl xis nd severl methods re used to obtin vlue of rel resistivity nd thickness of lyers. These methods re given below;. Successive pproximtion methods b. Asymptotic estimte c. Rough estimte methods d. Estimte with model curves d.1. Estimte with double lyer model curves d.2. Estimte with three lyer model curves d.3. Estimte with ssistnce point crds e. Direct comment methods While Successive pproximtion methods hve been the lest used methods, Asymptotic estimte were the most used methods until computers improved (Coşkun, 2005). Progrms, such s IPI2WIN (Bobchow, 2002), RES2DINV (Loke, 1997) re widely used tody. www.tojst.net Copyright The Online Journl of Science nd Technology 70

Field Study Wenner 1 Wenner 2 Figure 2. Mesuring points Wenner rrys (Multiple-electrode resistivity method) is implemented in two profiles (Figures 2-3-4). 39 electrodes were used, electrode spcing ws selected s 10 meters. Totlly, 380 meters of opening is mde for Wenner 1. 42 electrodes were used, electrode spcing ws selected s 10 meters. Totlly, 410 meters of opening is mde for Wenner 2. Figure 3. Wenner 1 www.tojst.net Copyright The Online Journl of Science nd Technology 71

Figure 4. Wenner 2 RES2DINV s Results After records tken from multichnnel resistivity device were trnsferred to computer in DAT formt, they were evluted by RES2DINV progrm for the two Wenner rrys (Figure 5 nd Figure 6). In the evlution, rel resistivity, lyer thicknesses nd lso possible well loction were obtined. Formtions re indicted in Figure 5 nd Figure 6. Alluvil formtions re vilble in different regions in Figure 5, Figure 6. For this resons possible well loction is chosen t x=85 m long Wenner 1 rry. Alluvil formtion Metmorphic rocks Pebble Sndstone Figure 5. RES2DINV s results for Wenner 1 rry Pebble Sndstone Metmorphic rocks Alluvil formtion Figure 6. RES2DINV s results for Wenner 2 rry www.tojst.net Copyright The Online Journl of Science nd Technology 72

Conclusion Considering the two-dimensionl mps in Figure 5 nd Figure 6, when resistivity vlues re smll, wter contining formtions re thought to be vilble. If the loction is indicted with blue color, it my be lluvil formtion. If the loction is indicted with yellow or green colors, it my be pebble, sndstone. Drk red unit is considered to be the metmorphic rocks. By referring to the two-dimensionl mps, comments cn be mde bout the loction of the well to be drilled. While determining the loction of well, it is necessry to consider drilling method nd formtion. The suggested well loction in this study lies long Wenner 1 line t round x=85 m which is locted wy from lluvil formtions. Mud rotry drilling method is preferred for suggested well. The reson for this method is tht it is more economicl thn others. References Bşokur, A.T. (2004), Düşey Elektrik Sondjı Verilerinin Yorumu. Ankr Üniversitesi Mühendislik Fkültesi, 225 s. Bobchow, A. (2002), "Ipi2win user s guide." Moscow Stte University. Coşkun, A. (2005), Yerltısu Kynklrının Terrmeter Ss 300C Cihzı ile Tespiti. Yüksek Lisns Tezi, Fırt Üniversitesi, Ypı Eğitimi Anbilim Dlı, Elzığ. Loke, M. H. (1997), "Res2DINV softwre user s mnul." University Sins Mlysi, Penng. Loke, M. H. (2000), Electricl Imging Surveys for Environmentl nd Engineering Studies, Prcticl Guide to 2- D nd 3-D Surveys: Mnul of RES2DINV nd RES3DINV softwre. www.tojst.net Copyright The Online Journl of Science nd Technology 73