Transient Radial Flow Toward a Well Aquifer Equation, based on assumptions becomes a 1D PDE for h(r,t) : Transient Radial Flow Toward a Well

Similar documents
11. HAFAT İş-Enerji Power of a force: Power in the ability of a force to do work

MEAN GRAVITY ALONG PLUMBLINE. University of New Brunswick, Department of Geodesy and Geomatics Engineering, Fredericton, N.B.

IrrItrol Products 2016 catalog

ELECTRIC & MAGNETIC FIELDS I (STATIC FIELDS) ELC 205A

Sections 3.1 and 3.4 Exponential Functions (Growth and Decay)

2. The units in which the rate of a chemical reaction in solution is measured are (could be); 4rate. sec L.sec

Part 2 KINEMATICS Motion in One and Two Dimensions Projectile Motion Circular Motion Kinematics Problems

Algebra 2A. Algebra 2A- Unit 5

WYSE Academic Challenge Sectional Mathematics 2006 Solution Set

Single Degree of Freedom System Forced Vibration

March 15. Induction and Inductance Chapter 31

Algebra-based Physics II

OBJECTIVE To investigate the parallel connection of R, L, and C. 1 Electricity & Electronics Constructor EEC470

Visco-elastic Layers

A) N B) 0.0 N C) N D) N E) N

The 37th International Physics Olympiad Singapore. Experimental Competition. Wednesday, 12 July, Sample Solution

Neutron Slowing Down Distances and Times in Hydrogenous Materials. Erin Boyd May 10, 2005

Chapter Finite Difference Method for Ordinary Differential Equations

Kinematics Review Outline

Lecture 17: Kinetics of Phase Growth in a Two-component System:

Example

A) (0.46 î ) N B) (0.17 î ) N

Electric Fields and Electric Forces

Lecture 3. Electrostatics

Microelectronics Circuit Analysis and Design. ac Equivalent Circuit for Common Emitter. Common Emitter with Time-Varying Input

Chapter 12: Kinematics of a Particle 12.8 CURVILINEAR MOTION: CYLINDRICAL COMPONENTS. u of the polar coordinate system are also shown in

The Gradient and Applications This unit is based on Sections 9.5 and 9.6, Chapter 9. All assigned readings and exercises are from the textbook

Notes on Inductance and Circuit Transients Joe Wolfe, Physics UNSW. Circuits with R and C. τ = RC = time constant

Maxwell Equations. Dr. Ray Kwok sjsu

( t) Steady Shear Flow Material Functions. Material function definitions. How do we predict material functions?

AP Physics 1 MC Practice Kinematics 1D

ME 3600 Control Systems Frequency Domain Analysis

Lecture 18: Kinetics of Phase Growth in a Two-component System: general kinetics analysis based on the dilute-solution approximation

β A Constant-G m Biasing

The Components of Vector B. The Components of Vector B. Vector Components. Component Method of Vector Addition. Vector Components

Prerna Tower, Road No 2, Contractors Area, Bistupur, Jamshedpur , Tel (0657) ,

5/20/2011. HITT An electron moves from point i to point f, in the direction of a uniform electric field. During this displacement:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

The Production of Polarization

5.1 Angles and Their Measure

Lecture 4 ( ) Some points of vertical motion: Here we assumed t 0 =0 and the y axis to be vertical.

A L A BA M A L A W R E V IE W

Chapter 7. Interference

156 There are 9 books stacked on a shelf. The thickness of each book is either 1 inch or 2

Module 4. Analysis of Statically Indeterminate Structures by the Direct Stiffness Method. Version 2 CE IIT, Kharagpur

, on the power of the transmitter P t fed to it, and on the distance R between the antenna and the observation point as. r r t

Control Volume Derivation

Convective Heat Transfer (6) Forced Convection (8) Martin Andersson

Graphs of Sine and Cosine Functions

1 2 U CV. K dq I dt J nqv d J V IR P VI

Lecture 4. Electrons and Holes in Semiconductors

Work, Energy, and Power. AP Physics C

AP Physics Electric Potential Energy

KINEMATICS OF RIGID BODIES

Unit-I (Feedback amplifiers) Features of feedback amplifiers. Presentation by: S.Karthie, Lecturer/ECE SSN College of Engineering

Today - Lecture 13. Today s lecture continue with rotations, torque, Note that chapters 11, 12, 13 all involve rotations

LEAP FROG TECHNIQUE. Operational Simulation of LC Ladder Filters ECEN 622 (ESS) TAMU-AMSC

Comparing Means: t-tests for One Sample & Two Related Samples

1 Widrow-Hoff Algorithm

Outline. Steady Heat Transfer with Conduction and Convection. Review Steady, 1-D, Review Heat Generation. Review Heat Generation II

ECE Spring Prof. David R. Jackson ECE Dept. Notes 5

Ch 13 Universal Gravitation

Integral Control via Bias Estimation

MATH 31B: MIDTERM 2 REVIEW. x 2 e x2 2x dx = 1. ue u du 2. x 2 e x2 e x2] + C 2. dx = x ln(x) 2 2. ln x dx = x ln x x + C. 2, or dx = 2u du.

Lecture #2 : Impedance matching for narrowband block

Example 11: The man shown in Figure (a) pulls on the cord with a force of 70

Experiment I Voltage Variation and Control

The Global Trade and Environment Model: GTEM

CHAPTER GAUSS'S LAW

The shortest path between two truths in the real domain passes through the complex domain. J. Hadamard

Announcements Candidates Visiting Next Monday 11 12:20 Class 4pm Research Talk Opportunity to learn a little about what physicists do

Propagation of Torsional Surface Waves. in Heterogeneous Half-Space. with Irregular Free Surface

Physics 201 Lecture 18

MECHANICS OF MATERIALS Poisson s Ratio

Lecture 4. Electrons and Holes in Semiconductors

Electrostatics. 1. Show does the force between two point charges change if the dielectric constant of the medium in which they are kept increase?

Electric Field, Potential Energy, & Voltage

dm dt = 1 V The number of moles in any volume is M = CV, where C = concentration in M/L V = liters. dcv v

An Open cycle and Closed cycle Gas Turbine Engines. Methods to improve the performance of simple gas turbine plants

Maximum Flow and Minimum Cut

Variance and Covariance Processes

AP Calculus AB Exam Review Sheet B - Session 1

CHAPTER 3. Section 1. Modeling Population Growth

Application of Fractional Calculus in Food Rheology E. Vozáry, Gy. Csima, L. Csapó, F. Mohos


Kinematics in 2-D (II)

Convective Heat Transfer (6) Forced Convection (8) Martin Andersson

Small Combustion Chamber. Combustion chamber area ratio

a. (1) Assume T = 20 ºC = 293 K. Apply Equation 2.22 to find the resistivity of the brass in the disk with

Assignment 6. Tyler Shendruk December 6, 2010

T h e C S E T I P r o j e c t

x i v x t a dx dt t x

Brace-Gatarek-Musiela model

GEOID-QUASIGEOID CORRECTION IN FORMULATION OF THE FUNDAMENTAL FORMULA OF PHYSICAL GEODESY. Robert Tenzer 1 Petr Vaníček 2

CHAPTER 24 GAUSS LAW

Physics 2212 GH Quiz #2 Solutions Spring 2016

Consider the simple circuit of Figure 1 in which a load impedance of r is connected to a voltage source. The no load voltage of r

Derivatives of Inverse Trig Functions

Physics Courseware Physics I Constant Acceleration

EM-2. 1 Coulomb s law, electric field, potential field, superposition q. Electric field of a point charge (1)

Transcription:

ansien Radial Flw wad a Well Aqife Eqain, based n assmpins becmes a D PDE f h(,) : -ansien flw in a hmgenes, ispic aqife -flly peneaing pmping well & infinie, hiznal, cnfined aqife f nifm hickness, hs essenially hiznal gndwae flw -flw symmey: adially symmeic flw Ss h h K Bnday cndiins - nd de in space PDE f h(,), need w BC s in space, a w s, say and ; in geneal we ll pick w, -One Diichle BC a infiniy and ne Neman a well adis, w, whee we assme we knw he pmping ae, w h h + x y S h h h S h + Diffsin eqains ansien Radial Flw wad a Well Ding a pmping es, we ypically mease dawdwn (as ppsed head), s le s se p he eqain in ems f dawdwn. Assme iniial penimeic sface is hiznal eveywhee. ( ) h h s s h h 0 0 h s, h s s s S s +

ansien Flw Well es Analysis slve, we need bnday cndiins and iniial cndiins I.C.: dawdwn is ze, head is nifm iniially. s 0 @ 0 B.C.s : nd de eqain, we need w BCs Use 0, Use cniniy a he pmping well, eaed as a line sink ( 0) A 0, h s KA KA A b s b K s ansien Flw Well es Analysis s b K s s ( ) s % lim& 0' ) $ (eqies cnsan ) s 0 as (f all )

ansien Flw Well es Analysis Slin: e s (, ) d whee $ S e d In mahemaics: he expnenial inegal In hydgelgy: he well fncin e W ( ) d s W ( ) heis Eqain ansien Flw Well es Analysis S : disance fm pmping well bsevain well : ime since pmping saed i $ W ( ) $ % $ ln $ $ Ele s cnsan 0.5775 i i i Vales f W() f vais vales f ae als ablaed in nmes efeences 3

ansien Flw Well es Analysis s W ( ) S O gal in aqife esing is find and S; slving he dawdwn eqain shwn abve f is pblemaic appeas wice. One way slve f is cmpae a pl f s vess a nmal gaphic slin (when pling daa fm nmes bsevain wells, we can ypically nmalize s by pling i as a fncin f / insead f ) ansien Flw Well es Analysis ype cve mehd f slving he heis Eqain Pl s vess ( s vess / ) n lg-lg pape On anhe shee f he same pape, pl W() vess / Why /? O pl has in he nmea, b in he eqain defining, appeas in he denmina S Mach daa (keep gaph axes paallel) Pick mach pin easies pick a mach pin wih simple nmbes he mach pin des n have fall n he cve. Slve eqain

ansien Flw Well es Analysis (Schwaz and Zhang, 003) ansien Flw Well es Analysis 500 m 3 /d 300 m (Schwaz and Zhang, 003) 5

ansien Flw Well es Analysis W ( ) ; 0.; min; s 0.78 m s W ( 3 m 500 m ) d 5 (0.78 m) d S m d 5 & min $ & $ d & % % 0 min $ % ( 300 m) ( 0.) 3.6 x 0-6 B S is had mease ypically shld nly be eped ne significan digi id es 3 x 0-6 ansien Flw Well es Analysis Assmpins: Hmgenes, ispic aqife (we sed a simplified fm f he cniniy eqain insead f inclding he ens f K) Cnsan wae ppeies Infinie aqife (sed develp bnday cndiin) Aqife is cnfined; cnfining nis ae impemeable All flw is adial and hiznal Wae is wihdawn insananesly wih decline in head (n ypically he case f leakage fm cnfining nis lw k lenses) Well diamee is infiniesimal (wae can be sed in he well be, leading delayed dawdwn) Cnsan dischage fm well 6

Head f NM Hydlgy Pgam C.E. Jacb (965-970) Hansh s men, well hydalics, hey f leaky aqifes 3 facly, inclding Fank is and William Bsae meged wih Gelgy and Gephysics fm Gesciences Depamen died f hea aack in 970; Hansh ened biefly ansien Radial Flw wad a Well Simplified appach ansien analysis: he Jacb Appximain A simplifying assmpin ha makes slving he heis eqain easie C. E. Jacb (90) Jacb appximain Jacb-Cpe eqain Cpe-Jacb mehd Cpe-Jacb saigh-line mehd i $ W ( ) $ % $ ln $ i i i 3 W ( ) 0.5776 ln + + +... 3 3 7

8 ansien Flw Well es Analysis ) ( 0.0, F W i i i i << $ % S S we need be lage f be small ( ) ln 0.577 s 0.0, F ( ) $ % & ' ' ( ) * * +, - - S ln ln.78 s. ansien Flw Well es Analysis ( ) $ % & ' ' ( ) * * +, - - S ln ln.78 s. ( ) ln ln ( ) lg ln.3 ( ) $ % & ' ( ) * +, + S ln 0.565 ln s - ( ) $ % & S 0.565 ln s '

ansien Flw Well es Analysis s ' & ln$ % ( 0.565) S Jacb-Cpe simplificain ( ) lg ln.3 ansien Flw Well es Analysis Sppse we mease dawdwn a w imes:, & -.5 * -.5 * s ' s $ ln+ ( + ln ' ln+ ( ' ln. %, S ), S ) s s [ ln ] s ln s ln 9

ansien Flw Well es Analysis s s ln Pick 0 ; head change pe lg cycle %s lc s - s s s 0 ln.3 ansien Flw Well es Analysis &.5 A, s 0 ln$ ' % S 0 ( ) ' &.5 ln$ % S &.5 0 ln$ % S e 0 &.5 ln$ % S e 0

ansien Flw Well es Analysis e 0 &.5 ln$ % S e.5 S ansien Flw Well es Analysis Semi-lg pl (dawdwn vess lg ime) 500 m 3 /d (Schwaz and Zhang, 003)

ansien Flw Well es Analysis.3.3 s lc 3 m ( 500 ) d 5 m d (.8m).5 S m.5 5 d S 300 m ( ) 3.min d ( 0min ) 3 x 0 6 ( ) ansien Flw Well es Analysis Remembe ha he Jacb-Cpe mehd is pedicaed n he vale f being less han ab 0.0 S 6 ( 300 m) ( 3 x 0 ) m d ( 5 )( 00 min ) d 0min 0.0 Accding y exbk, 0.0 is clse 0.0, and i is he maximm vale f ve lg cycle f daa, s he Jacb-Cpe mehd is accepable.

ansien Flw Well es Analysis Insead f lking a dawdwn in ne well wih ime, we can als lk a dawdwn meased a he same ime in w me wells: Disance-dawdwn mehd s ( ) s ( ) Meased a he same ime (Schwaz and Zhang, 003) ansien Flw Well es Analysis s ' &.5 ln$ % S '.5 $ '.5 $ ln % ln ln % i S & & S i &,.5 ),.5 ) s - s $ ln* ' - ln - ln* ' + ln. % + S ( + S ( s s ln + ln [ ] s s ln 3

s s ln ansien Flw Well es Analysis s ( s ' $ % ln & ' $ % ln & ( s Disance-dawdwn fmla.5 S ansien Flw Well es Analysis calclae S, we se a fmla simila ha sed in he ime-dawdwn fmla. We sed becase i ld s he ime a which dawdwn was ze; f a disance-dawdwn fmla, we wan knw he disance a which s 0 (we ll call i ) he es f he deivain is he same as f he disance-dawdwn mehd, s we ll jmp igh he fmla:.5 S Disance-dawdwn fmla Remembe ha hee ae w fmlae f applying he Jacb-Cpe mehd ime-dawdwn and disance-dawdwn

ansien Flw Well es Analysis 0 gpm 9. f 3 /min,350 f 3 /d Dawdwns meased a 0 min (Schwaz and Zhang, 003) ansien Flw Well es Analysis.3 s lc.3 3 f ( 350 ) d 987 f d ( 5.7 f).5 S S S ( ) 0min f.5 987 d 987 f ( ) 9 x 0 5 ( 900 f) ( ) d [( )( 0min) ] 0.5 ( 900 f) 9 x 0 5 ( d ) 0min 5

ansien Flw Well es Analysis ' $ % ln & s ' 0 $ ln % & ( h ( lc Disance-dawdwn fmla hiem eqain ansien Flw Well es Analysis ' $ % ln & s ' 0 $ ln % & ( h ( lc Disance-dawdwn fmla hiem eqain Fmlain f he DDF is idenical ha f he hiem eqain his means ha when Jacb s appximain applies, he diffeence in dawdwn beween any w pins sabilizes he penimeic sface is being dawn dwn nifmly wih ime 6

ansien Flw Well es Analysis Why desn his apply in ealy ime? Sabiliy is achieved when he slpe f he penimeic sface has he cec gadien spply he well We ae cnsanly lweing he penimeic sface we ms d his elease wae fm sage. In ealy ime, wae is nly eleased fm sage nea he well, b as ime ges lage, wae is spplied fm a lage aea. As he vlme f aqife spplying inceases, less change in head is needed yield he same amn f wae, s he cne f depessin inceases a a deceasing ae. 7

Hw d cnes f depessin vay fm nmal wih changes in aqife ppeies, and hw des his affec he dawdwn hydgaph f a mniing well? Hw migh dawdwn hydgaphs vay fm nmal in siains whee pmping es assmpins ae n me? Vaiains in dawdwn Cne f depessin as a fncin f ime min 0 min 0 min 000 gpd/f; S 0 - ; 300 gpm (Hall, 996) 8

Vaiains in dawdwn Vaiains wih pmping ae 00 gpm 300 gpm 600 gpm 000 gpd/f; S 0-; 00 f (Hall, 996) Vaiains in dawdwn Vaiains wih saiviy S 0. S 0.00000 S 0.000 S 0. 0. 0. 0 00 00 0. 0. S 0-; 300 gpm; 00 f 0 00 00 (Hall, 996) 9

Vaiains in dawdwn Vaiains wih ansmissiviy 500 gpd/f, 000 gpd/f,500 gpd/f S 000 gpd/f; 300 gpm; 00 f (Hall, 996) Sensiiviy and S smmaized Lw S: Cne shape emains simila, he cne is js bigge High S: Feeze and Chey (979) Lw : igh, deep cne; High : wide, shallw cne Lw S, ceaes geae dawdwn pdce a given vlme f wae; cne f depessin shape says he same Lw : had ansmi wae, ms efficien pdcin is fm clse in, ceaes seepe gadiens f s pdce wae fm clse in. 0