Electroseismics for CO 2 storage and hydrocarbon reservoirs. Fabio I. Zyserman
|
|
- Bennett Newton
- 5 years ago
- Views:
Transcription
1 Electroseismics for CO 2 storage and hydrocarbon reservoirs Fabio I. Zyserman EOST, March 12th, 213
2 Electroseismics: field experiments Thompson et al, 27 Electroseismic response in gas and oil reservoirs from about 15 m depth. Signal between two and six orders of magnitude less than ambient noise. It is necessary to optimize the source power ( megawatt), the injected signal (few thousands A), and the detection equipment (digital accelerometers). EOST, March 12th, 213
3 Electroseismic Modeling I Glover and Jackson, TLE, 21 When an applied electric field acts on an electrolyte saturated porous material, besides driving σe, it acts as a body force on the excess charge, giving rise to a net fluid filtration; this is called electro-osmosis electroseismic phenomena. Reciprocally, an applied pressure gradient generates an electric current; this is called electro-filtration seismoelectric phenomena. EOST, March 12th, 213
4 Electroseismic modeling II Assuming an e +iωt time dependence, Pride (1994) proposed (σ + iɛω)e H + L(ω)ηκ 1 [ iωu f ] L(ω)E = J ext e, E + iωµh = J ext m, ω 2 ρ b u s ω 2 ρ f u f τ(u) = F (s), ω 2 ρ f u s + ηκ 1 [ iωu f ] L(ω)E + p f = F (f ), τ lm (u) = 2G ε lm (u s ) + δ lm (λ c u s + αk av u f ), p f (u) = αk av u s K av u f. φ porosity, ρ s, ρ f solid and fluid densities, ρ b = (1 φ)ρ s + φρ f, η fluid viscosity, κ(ω) dinamic permeability In the constitutive equations λ c = K c 2/3G and K c = K m + α 2 K av, [ α = 1 Km α φ, K av = + φ ] 1. K s K s K f K s, K m and K f : bulk moduli of the solid grains, the dry matrix and the fluid. EOST, March 12th, 213
5 Electroseismic modeling III 2D Sources and Modes Infinite solenoid: Jm ext generates electromagnetic fields (E x (x, z), E z (x, z)), and H y (x, z), coupled with solid displacements (ux(x, s z), uz(x, s z)) and fluid displacements (ux(x, f z), uz(x, f z)). This is the so-called PSVTM-mode, in which compressional and vertically polarized shear seismic waves (PSV-waves) are present. Infinite current line: Je ext generates electromagnetic fields (H x (x, z), H z (x, z)) and E y (x, z), coupled with solid displacements uy(x, s z) and fluid displacements uy(x, f z). This is known as the SHTE-mode, where only horizontally polarized seismic waves (SH-waves) are present. EOST, March 12th, 213
6 Electroseismic modeling IV Some assumptions We work in the seismic frequency range, then Re(η/κ(ω)) η/κ and 1 ω Im(η/κ(ω)) g = 1.5 ρ f T φ, T being the tortuosity factor so that the low-frequency Biot s equations are recovered. The electroseismic coupling coefficient L is assumed to be frequency independent, L = φ T F (s) = F (f ) =, and ωε/σ 1 ε κ f ζ (1 2 d η Λ ), We consider lossy media using Liu s model. Electro-filtration feedback negligible; this decouples the EM fields from the poroviscoelastic response. (This makes calculations easier) EOST, March 12th, 213
7 Scheme of the Finite Element Procedure Create a partition of the domain (elements). Transform original equations into a "weak form". Choose appropriate polynomial functions to approximate the solution in each element, (dofs). Transform the weak form into a linear sistem, and solve it. ( unknowns) EOST, March 12th, 213
8 Single horizontal layer 7 m 4 m 4 m Medium 1 Medium 2 (layer) σ (S/m).1.1 φ ( ).2.33 Ks (Pa) vp (m/s) vs (m/s) ρs (kg/m 3 ) k (m 2 ) L Q ( ) 9 9 ρf (kg/m 3 ) 1 1 η (kg/(m s)).1.1 Kf (Pa) Sf ( ) 1 1 EOST, March 12th, 213
9 Single horizontal layer SHTE-mode a b t=.3 s b t=.4 s c a b t=.1 s c d t=.2 s EOST, March 12th, 213
10 Single horizontal layer Different layer widths, SHTE-mode 4e-1 3e-1 2e-1 4 m Acceleration [m/s^2] 5e-1 4e-1 3e-1 2e-1 1 m 1e-1-1e-1-2e-1-3e-1 1e-1-1e-1-2e-1-3e-1-4e-1-5e-1-4e-1-6e e-1 5 m 3e-1 2 m 4e-1 2e-1 2e-1 1e-1-2e-1-1e-1-4e-1-2e-1-6e-1-3e-1-8e-1-4e t[s] t[s] EOST, March 12th, 213
11 Single horizontal layer PSTVM-mode 1e-12 8e-13 Acceleration traces - X component [m/s 2 ] -7 m m 7 m 1e-11 8e-12 Acceleration traces - Z component [m/s 2 ] -7 m m 7 m 6e-13 6e-12 4e-13 4e-12 2e-13 2e-12-2e-12-2e-13-4e-12-4e-13-6e-12-6e-13-8e-12-8e t [s] -1e t [s] EOST, March 12th, 213
12 Single horizontal layer PSTVM-mode 2e e-14 Div and Curl of an acceleration trace SV Div Curl 1e-14 5e-15-5e-15-1e e-14-2e e-14 P t [s] EOST, March 12th, 213
13 Wedge First Layer Second layer Third layer Wedge σ c (S/m) φ ( ) Ks (Pa) m vp (m/s) vs (m/s) ρs (kg/m 3 ) m 2 m k (m 2 ) L m 65 m 4 m Q ( ) ρf (kg/m 3 ) η (kg/(m s)) Kf (Pa) Sf ( ) EOST, March 12th, 213
14 Surface gather x-component acceleration Traces 2 4 Traces a a time(s).2 b time(s).2 b.3 c.3 d.4.4 EOST, March 12th, 213
15 Surface gather z-component acceleration Traces 2 4 Traces time(s).2 a time(s).2 a b b.3 c EOST, March 12th, 213
16 Methane hydrates (GH)......form stable ice-like crystals in permafrost regions and beneath the ocean floor along continental margins. Ellis, 28...are considered as a potentially huge energy resource....have the highest energy density of any naturally occurring form of methane (about 16 times that of methane gas)...decrease the electrical conductivity of the medium. EOST, March 12th, 213
17 Composite Media We use an extended Biot theory for composite matrix rocks with non uniform porosity distributions. The solid matrix can be formed by mixtures of different mineral grains. A fraction of the GH (ice) is assumed to form a second matrix occupying the pore space, and the rest of it is assumed to cement the mineral grains. LettingV = (V c gh + V nc gh ) + V mg + V f, V p = V V mg, φ a = V p /V ; and C gh = V c gh /V gh,cementation coefficient, S gh = V gh /V p GH saturation, all (Biot) model parameters are obtained in terms of C gh, S gh, φ a, K j mg and µ j mg forming the solid matrix; and K gh, µ gh. EOST, March 12th, 213
18 The Model x z 48 m Permafrost 31 m 25 m 6 m Sandstone 1 Sandstone 2 21 m GH Saturated Slab Permafrost Sandstone Slab (Sgh =.1) Slab (Sgh =.8) φ Vp m/s Vs m/s σ S/m L A/(Pa m) EOST, March 12th, 213
19 SHTE-mode Acceleration surface gather Surface detectors locations (m) Surface detectors locations (m) A.4 A Time (s) Time (s).6.6 B.8.8 C Without hydrate With hydrate EOST, March 12th, 213
20 SHTE-mode Acceleration well gather 8 8 F I 6 C D 6 J E C H Depth (m) A 4 Depth (m) 4 A B B G Time (s) Without hydrate Time (s) With hydrate EOST, March 12th, 213
21 SHTE-mode Acceleration surface trace for S gh =.1 and S gh =.8 5e-9 4e-9 3e-9 2e-9 1e-9-1e-9-2e-9-3e-9-4e-9 Solid Acceleration (m/s 2 ) S gh =.1 S gh =.8-5e t (s) 8e-1 6e-1 4e-1 2e-1-2e-1-4e-1-6e-1 Solid Acceleration (m/s 2 ) S gh =.1 S gh =.8-8e t (s) EOST, March 12th, 213
22 PSVTM-mode Acceleration (x-component) well gather 8 8 E F 6 C D 6 C D G Depth (m) 4 A Depth (m) 4 A H 2 B 2 B Time (s) Without hydrate Time (s) With hydrate EOST, March 12th, 213
23 PSVTM-mode Acceleration well traces for for S gh =.1 and S gh =.8 X-component solid acceleration (m/s 2 ) 6e-13 Well Trace #61 4e-13 2e-13-2e-13-4e-13-6e-13 S gh =. S gh =.1 S gh =.8-8e t (s) Z-component solid acceleration (m/s 2 ) 2.5e-12 S gh =. 2e-12 S gh =.1 S gh =.8 1.5e-12 1e-12 5e-13-5e-13-1e e-12 Well Trace #61-2e t (s) EOST, March 12th, 213
24 CO 2 storage monitoring The Model 1 mts (a) z Geophones (b) Air x Electromagnetic source 7mts 1 mts 4 mts 4 mts Seal layer Reservoir Homogeneous Soil EOST, March 12th, 213
25 CO 2 storage monitoring PSVTM mode, acceleration traces 4e-15 3e-15 x-component of acceleration traces [m/s 2 ] % 15 % 55 % 95 % 1.5e-14 1e-14 z-component of acceleration traces [m/s 2 ] % 15 % 55 % 95 % 2e-15 1e-15 5e-15-1e-15-5e-15-2e-15-3e-15-1e-14-4e t [s] -1.5e t [s] EOST, March 12th, 213
26 CO 2 storage monitoring PSVTM mode, acceleration traces 4e-15 3e-15 x-component of acceleration traces [m/s 2 ] (a) (b) 1.2e-14 1e-14 8e-15 z-component of acceleration traces [m/s 2 ] (a) (b) 2e-15 6e-15 4e-15 1e-15 2e-15-2e-15-1e-15-4e-15-2e-15-6e-15-8e-15-3e t [s] -1e t [s] EOST, March 12th, 213
27 Summary We have developed a numerical tool to simulate electroseismic (and seismoelectric) phenomena. It was observed that the response is sensitive to changes in fluid conductivities (mixtures of gas an brine). We have shown that methane hydrates can be detected by means of electroseismics on land in permafrost regions. We have observed that the seismic response is sensitive to methane hydrate concentration Preliminar results indicate that it could be very interesting to consider electroseismics/ seismoelectrics as a monitoring tools for CO 2 storage sites. Heaps of work ahead! EOST, March 12th, 213
Coupled seismoelectric wave propagation in porous media. Mehran Gharibi Robert R. Stewart Laurence R. Bentley
Coupled seismoelectric wave propagation in porous media Mehran Gharibi Robert R. Stewart Laurence R. Bentley a Introduction Seismic waves induce electric and magnetic fields. Conversion of acoustic energy
More informationTheoretical and Experimental Studies of Seismoelectric Conversions in Boreholes
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 3, No. 1, pp. 109-120 Commun. Comput. Phys. January 2008 Theoretical and Experimental Studies of Seismoelectric Conversions in Boreholes Zhenya Zhu 1,, Shihong
More informationAcoustics in Porous Media.
Acoustics in Porous Media. JUAN E. SANTOS work in collaboration with J. M. Carcione, S. Picotti, P. M. Gauzellino and R. Martinez Corredor. Department of Mathematics, Purdue University, W. Lafayette, Indiana,
More informationHydrogeophysics - Seismics
Hydrogeophysics - Seismics Matthias Zillmer EOST-ULP p. 1 Table of contents SH polarized shear waves: Seismic source Case study: porosity of an aquifer Seismic velocities for porous media: The Frenkel-Biot-Gassmann
More informationSite Characterization & Hydrogeophysics
Site Characterization & Hydrogeophysics (Source: Matthew Becker, California State University) Site Characterization Definition: quantitative description of the hydraulic, geologic, and chemical properties
More informationJuan E. Santos a,b,c, Gabriela B. Savioli a and Robiel Martínez Corredor c a
Juan E. Santos a,b,c, Gabriela B. Savioli a and Robiel Martínez Corredor c a Universidad de Buenos Aires, Fac. Ing., IGPUBA, ARGENTINA b Department of Mathematics, Purdue University, USA c Universidad
More informationA look into Gassmann s Equation
A look into Gassmann s Equation Nawras Al-Khateb, CHORUS Heavy Oil Consortium, Department of Geoscience, University of Calgary nawras.alkhateb@ucalgary.ca Summary By describing the influence of the pore
More informationPropagation of seismic-induced electromagnetic waves in a semi-infinite porous medium: A Fourier transform approach
Seismoelectric Simulations Propagation of seismic-induced electromagnetic waves in a semi-infinite porous medium: A Fourier transform approach R.Arief Budiman, Aqsha Aqsha, Mehran Gharibi and Robert R.
More informationWe Simultaneous Joint Inversion of Electromagnetic and Seismic Full-waveform Data - A Sensitivity Analysis to Biot Parameter
We-09-04 Simultaneous Joint Inversion of Electromagnetic and Seismic Full-waveform Data - A Sensitivity Analysis to Biot Parameter J. Giraud* (WesternGeco Geosolutions), M. De Stefano (WesternGeco Geosolutions)
More information3D simulations of an injection test done into an unsaturated porous and fractured limestone
3D simulations of an injection test done into an unsaturated porous and fractured limestone A. Thoraval *, Y. Guglielmi, F. Cappa INERIS, Ecole des Mines de Nancy, FRANCE *Corresponding author: Ecole des
More informationEstimating Permeability from Acoustic Velocity and Formation Resistivity Factor
5th Conference & Exposition on Petroleum Geophysics, Hyderabad-2004, India PP 582-587 and Formation Resistivity Factor Majid Nabi-Bidhendi Institute of Geophysics, University of Tehran, P.O. Box 14155-6466,
More informationScale-up in poroelastic systems and applications to reservoirs
Scale-up in poroelastic systems and applications to reservoirs James G. Berryman 1 ABSTRACT A fundamental problem of heterogeneous systems is that the macroscale behavior is not necessarily well-described
More informationReservoir Rock Properties COPYRIGHT. Sources and Seals Porosity and Permeability. This section will cover the following learning objectives:
Learning Objectives Reservoir Rock Properties Core Sources and Seals Porosity and Permeability This section will cover the following learning objectives: Explain why petroleum fluids are found in underground
More informationFluid Flow Fluid Flow and Permeability
and Permeability 215 Viscosity describes the shear stresses that develop in a flowing fluid. V z Stationary Fluid Velocity Profile x Shear stress in the fluid is proportional to the fluid velocity gradient.
More informationCold production footprints of heavy oil on time-lapse seismology: Lloydminster field, Alberta
Cold production footprints of heavy oil on time-lapse seismology: Lloydminster field, Alberta Sandy Chen, Laurence R. Lines, Joan Embleton, P.F. Daley, and Larry Mayo * ABSTRACT The simultaneous extraction
More informationScienceDirect. Model-based assessment of seismic monitoring of CO 2 in a CCS project in Alberta, Canada, including a poroelastic approach
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 63 (2014 ) 4305 4312 GHGT-12 Model-based assessment of seismic monitoring of CO 2 in a CCS project in Alberta, Canada, including
More informationWave Propagation in Fractured Poroelastic Media
Wave Propagation in Fractured Poroelastic Media WCCM, MS170: Advanced Computational Techniques in Geophysical Sciences, Barcelona, Spain, July 2014 Juan E. Santos Instituto del Gas y del Petróleo (IGPUBA),
More information13th International Symposium on Nondestructive Characterization of Materials (NDCM-XIII), May 2013, Le Mans, France
3th International Symposium on Nondestructive Characterization of Materials (NDCM-XIII), 2-24 May 23, Le Mans, France www.ndt.net/?id=5532 Biot waves in porous ceramic plates : influence of boundary conditions
More informationLawrence Berkeley National Laboratory
Lawrence Berkeley National Laboratory Peer Reviewed Title: Fracture permeability and seismic wave scattering--poroelastic linear-slip interface model for heterogeneous fractures Author: Nakagawa, S. Publication
More informationPartial Saturation Fluid Substitution with Partial Saturation
Fluid Substitution with 261 5 4.5 ρ fluid S w ρ w + S o ρ o + S g ρ g Vp (km/s) 4 3.5 K fluid S w K w + S o K o + S g K g Patchy Saturation Drainage 3 2.5 2 Fine-scale mixing 1 = S w + S o + S g K fluid
More informationGG 454 March 19, EFFECTIVE STRESS AND MOHR-COULOMB FAILURE (26)
GG 454 March 19, 2002 1 EFFECTIVE STRESS AND MOHR-COULOMB FAILURE (26) I Main Topics A Driving and resisting stresses at the base of an inclined block B Factor of safety C Effective stress D Mohr-Coulomb
More informationFour-D seismic monitoring: Blackfoot reservoir feasibility
Four-D seismic monitoring Four-D seismic monitoring: Blackfoot reservoir feasibility Laurence R. Bentley, John Zhang and Han-xing Lu ABSTRACT The Blackfoot reservoir has been analysed to determine the
More informationA Thesis. Presented to. the Faculty of the Department of Earth and Atmospheric Sciences. University of Houston. In Partial Fulfillment
METHODS OF FITTING COMPRESSIONAL AND SHEAR WAVE VELOCITIES VERSUS SATURATION CURVES AND THE INTERPRETATION OF LABORATORY VELOCITY MEASUREMENTS IN PARTIALLY SATURATED ROCKS ------------------------------------------------------------
More informationwater L v i Chapter 4 Saturation
4. Resistivity The presence of hydrocarbons is identified by the electrical resistance of the formation. These electrical properties of rocks depend on the pore geometry and fluid distribution. That is,
More informationUncertainties in rock pore compressibility and effects on time lapse seismic modeling An application to Norne field
Uncertainties in rock pore compressibility and effects on time lapse seismic modeling An application to Norne field Amit Suman and Tapan Mukerji Department of Energy Resources Engineering Stanford University
More informationSeismic Velocity Dispersion and the Petrophysical Properties of Porous Media
Seismic Velocity Dispersion and the Petrophysical Properties of Porous Media L. Flora Sun* University of Toronto, Toronto, ON lsun@physics.utoronto.ca and B. Milkereit University of Toronto, Toronto, ON,
More informationSRC software. Rock physics modelling tools for analyzing and predicting geophysical reservoir properties
SRC software Rock physics modelling tools for analyzing and predicting geophysical reservoir properties Outline About SRC software. Introduction to rock modelling. Rock modelling program structure. Examples
More informationResearch Article Study on p-wave Attenuation in Hydrate-Bearing Sediments Based on BISQ Model
Geological Research Volume 23, Article ID 76579, 8 pages http://dx.doi.org/.55/23/76579 Research Article Study on p-wave Attenuation in Hydrate-Bearing Sediments Based on BISQ Model Chuanhui Li, Kai Feng,
More informationWELLBORE STABILITY ANALYSIS IN CHEMICALLY ACTIVE SHALE FORMATIONS
S911 WELLBORE SABILIY ANALYSIS IN CHEMICALLY ACIVE SHALE FORMAIONS by Xiang-Chao SHI *, Xu YANG, Ying-Feng MENG, and Gao LI State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest
More informationscaling parameters of laboratory modelling of
Influence of low and high IFT fluid systems on scaling parameters of laboratory modelling of CO 2 injection into saline aquifers The 6 th Trondheim Conference on CO 2 Capture, Transport, and Storage 14-16
More informationEvaporation-driven soil salinization
Evaporation-driven soil salinization Vishal Jambhekar 1 Karen Schmid 1, Rainer Helmig 1,Sorin Pop 2 and Nima Shokri 3 1 Department of Hydromechanics and Hydrosystemmodeling, University of Stuttgart 2 Department
More informationShear Wave Velocity Estimation Utilizing Wireline Logs for a Carbonate Reservoir, South-West Iran
Iranian Int. J. Sci. 4(2), 2003, p. 209-221 Shear Wave Velocity Estimation Utilizing Wireline Logs for a Carbonate Reservoir, South-West Iran Eskandari, H. 1, Rezaee, M.R., 2 Javaherian, A., 3 and Mohammadnia,
More informationMethane hydrate rock physics models for the Blake Outer Ridge
Stanford Exploration Project, Report 80, May 15, 2001, pages 1 307 Methane hydrate rock physics models for the Blake Outer Ridge Christine Ecker 1 ABSTRACT Seismic analyses of methane hydrate data from
More informationGeophysics 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 informationSENSITIVITY ANALYSIS OF THE PETROPHYSICAL PROPERTIES VARIATIONS ON THE SEISMIC RESPONSE OF A CO2 STORAGE SITE. Juan E. Santos
SENSITIVITY ANALYSIS OF THE PETROPHYSICAL PROPERTIES VARIATIONS ON THE SEISMIC RESPONSE OF A CO2 STORAGE SITE Juan E. Santos Instituto del Gas y del Petróleo, Facultad de Ingeniería UBA and Department
More informationComparative review of theoretical models for elastic wave attenuation and dispersion in partially saturated rocks
1 Comparative review of theoretical models for elastic wave attenuation and dispersion in partially saturated rocks J. Toms a*, T. Müller b, R. Ciz c and B. Gurevich a,c a Dept. Exp. Geophys., Curtin University
More informationIDENTIFYING PATCHY SATURATION FROM WELL LOGS Short Note. / K s. + K f., G Dry. = G / ρ, (2)
IDENTIFYING PATCHY SATURATION FROM WELL LOGS Short Note JACK DVORKIN, DAN MOOS, JAMES PACKWOOD, AND AMOS NUR DEPARTMENT OF GEOPHYSICS, STANFORD UNIVERSITY January 5, 2001 INTRODUCTION Gassmann's (1951)
More informationROCK PHYSICS DIAGNOSTICS OF NORTH SEA SANDS: LINK BETWEEN MICROSTRUCTURE AND SEISMIC PROPERTIES ABSTRACT
ROCK PHYSICS DIAGNOSTICS OF NORTH SEA SANDS: LINK BETWEEN MICROSTRUCTURE AND SEISMIC PROPERTIES PER AVSETH, JACK DVORKIN, AND GARY MAVKO Department of Geophysics, Stanford University, CA 94305-2215, USA
More informationABSTRACT 1. INTRODUCTION
ABSTRACT In this paper the three-dimensional transient wave propagation is investigated due to a point force applied at the interface of a fluid and a poroelastic solid. Using the total response, it is
More informationImproved Exploration, Appraisal and Production Monitoring with Multi-Transient EM Solutions
Improved Exploration, Appraisal and Production Monitoring with Multi-Transient EM Solutions Folke Engelmark* PGS Multi-Transient EM, Asia-Pacific, Singapore folke.engelmark@pgs.com Summary Successful as
More information6. THE BOREHOLE ENVIRONMENT. 6.1 Introduction. 6.2 Overburden Pressures
6. THE BOREHOLE ENVIRONMENT 6.1 Introduction Wireline logging has a single clearly defined purpose: to give accurate and representative data on the physical properties of the rock formations and fluids
More informationPart II: Self Potential Method and Induced Polarization (IP)
Part II: Self Potential Method and Induced Polarization (IP) Self-potential method (passive) Self-potential mechanism Measurement of self potentials and interpretation Induced polarization method (active)
More informationSUMMARY NUMERICAL SIMULATION OF BRINE-CO 2 FLOW AND WAVE PROPAGATION
Sensitivity analysis of the petrophysical properties variations on the seismic response of a CO 2 storage site Juan E. Santos, Universidad de Buenos Aires, Facultad Ingeniería, Instituto del Gas y del
More informationProceedings, International Snow Science Workshop, Banff, 2014
SNOW SLAB FAILURE DUE TO BIOT-TYPE ACOUSTIC WAVE PROPAGATION Rolf Sidler 2 Simon Fraser University, 8888 University Drive, Burnaby V5A1S6, Canada ABSTRACT: Even though seismic methods are among the most
More informationROCK PHYSICS MODELING FOR LITHOLOGY PREDICTION USING HERTZ- MINDLIN THEORY
ROCK PHYSICS MODELING FOR LITHOLOGY PREDICTION USING HERTZ- MINDLIN THEORY Ida Ayu PURNAMASARI*, Hilfan KHAIRY, Abdelazis Lotfy ABDELDAYEM Geoscience and Petroleum Engineering Department Universiti Teknologi
More informationUvA-DARE (Digital Academic Repository) Electrokinetics in porous media Luong, D.T. Link to publication
UvA-DARE (Digital Academic Repository) Electrokinetics in porous media Luong, D.T. Link to publication Citation for published version (APA): Luong, D. T. (2014). Electrokinetics in porous media General
More informationSeismic behaviour of CO 2 saturated Fontainebleau sandstone under in situ conditions
Seismic behaviour of CO 2 urated Fontainebleau sandstone under in situ conditions Md Mizanul Huq Chowdhury*, University of Alberta, Edmonton, AB, Canada, mhchowdh@ualberta.ca and Douglas R. Schmitt, University
More informationSome consideration about fluid substitution without shear wave velocity Fuyong Yan*, De-Hua Han, Rock Physics Lab, University of Houston
ain enu Some consideration about fluid substitution without shear wave velocity Fuyong Yan*, De-Hua Han, Rock Physics Lab, University of Houston Summary When S-wave velocity is absent, approximate Gassmann
More informationA thermo-hydro-mechanically coupled analysis of clay using a thermo-elasto-viscoplastic model
JHUWS05 A thermo-hydro-mechanically coupled analysis of clay using a thermo-elasto-viscoplastic model by F. Oka, S. Kimoto, Y.-S. Kim, N. Takada Department of Civil & Earth Resources Engineering, Kyoto
More informationCorrelation Between Resistivity Index, Capillary Pressure and Relative Permeability
Proceedings World Geothermal Congress 2010 Bali, Indonesia, 25-29 April 2010 Correlation Between Resistivity Index, Capillary Pressure and Kewen Li Stanford Geothermal Program, Stanford University, Stanford,
More informationinter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE
Copyright SFA - InterNoise 2000 1 inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering 27-30 August 2000, Nice, FRANCE I-INCE Classification: 7.5 IMPEDANCE CONSIDERATION
More informationAcoustic landmine detection: a 3D poroelastic model
Acoustic landmine detection: a 3D poroelastic model Y. Q. Zeng a and Q. H. Liu a a Department of Electrical and Computer Engineering, Duke University, Durham, North Carolina 2778, U.S.A. ABSTRACT Acoustic
More informationEffects of Fracture Parameters in an Anisotropy Model on P-Wave Azimuthal Amplitude Responses
PROC. ITB Eng. Science Vol. 38 B, No. 2, 2006, 159-170 159 Effects of Fracture Parameters in an Anisotropy Model on P-Wave Azimuthal Amplitude Responses Fatkhan Program Studi Teknik Geofisika FIKTM-ITB
More informationCompeting Effect of Pore Fluid and Texture -- Case Study
Competing Effect of Pore Fluid and Texture -- Case Study Depth (m) Sw Sxo. m Poisson's Ratio.. GOC.1 5 7 8 9 P-Impedance OWC 15 GR.. RHOB.5 1 Saturation...5. 1. 1.5 Vs (km/s).. Poisson's Ratio 5 7 P-Impedance
More informationAppendix: Nomenclature
Appendix: Nomenclature 209 Appendix: Nomenclature Chapter 2: Roman letters A amplitude of surface wave m A j (complex) amplitude of tidal component [various] a right ascension rad or a magnitude of tidal
More informationRock Physics Laboratory at the PIRC (rooms 1-121, B-185)
Rock Physics Laboratory at the PIRC (rooms 1-121, B-185) Description: The Rock Physics Lab at The PIRC examines the properties of rocks throughout experimental, imagining and numerical techniques. It focuses
More informationSPE These in turn can be used to estimate mechanical properties.
SPE 96112 Pressure Effects on Porosity-Log Responses Using Rock Physics Modeling: Implications on Geophysical and Engineering Models as Reservoir Pressure Decreases Michael Holmes, SPE, Digital Formation,
More informationSEG/New Orleans 2006 Annual Meeting
On the applicability of Gassmann model in carbonates Ravi Sharma*, Manika Prasad and Ganpat Surve (Indian Institute of Technology, Bombay), G C Katiyar (Third Eye Centre, Oil and Natural Gas Corporation
More informationKey Laboratory of Geo-detection (China University of Geosciences, Beijing), Ministry of Education, Beijing , China
Pet.Sci.(20118:43-48 DOI 10.1007/s12182-011-0113-5 43 Experimental study of the relationship between fluid density and saturation and sonic wave velocity of rock samples from the WXS Depression, South
More informationAgain we will consider the following one dimensional slab of porous material:
page 1 of 7 REVIEW OF BASIC STEPS IN DERIVATION OF FLOW EQUATIONS Generally speaking, flow equations for flow in porous materials are based on a set of mass, momentum and energy conservation equations,
More informationNumerical simulation of the Biot slow wave in water-saturated Nivelsteiner Sandstone
GEOPHYSICS, VOL. 66, NO. 3 (MAY-JUNE 2001; P. 890 896, 4 FIGS., 2 TABLES. Short Note Numerical simulation of the Biot slow wave in water-saturated Nivelsteiner Sstone Børge Arntsen José M. Carcione INTRODUCTION
More informationHydrological geophysical relationships
International PhD Course in HYDROGEOPHYSICS Hydrological geophysical relationships Andrew Binley Lancaster University Overview In the course we will concentrate on electrical, electromagnetic and radar
More informationEffective Medium Theories Effective Medium Models
Effective Medium Models 3 Bounding Methods for Estimating Effective Elastic Moduli For many reasons we would like to be able to model or estimate the effective elastic moduli of rocks in terms of the properties
More information42. POROSITY AND VELOCITY VS. DEPTH AND EFFECTIVE STRESS IN CARBONATE SEDIMENTS 1
Duncan, R. A., Backman, J., Peterson, L. C, et al., 1990 Proceedings of the Ocean Drilling Program, Scientific Results, Vol. 115 42. POROSITY AND VELOCITY VS. DEPTH AND EFFECTIVE STRESS IN ONATE SEDIMENTS
More informationP- and S-Wave Velocity Measurements and Pressure Sensitivity Analysis of AVA Response
P- and S-Wave Velocity Measurements and Pressure Sensitivity Analysis of AVA Response Tiewei He* University of Alberta, Edmonton, Alberta, Canada tieweihe@phys.ualberta.ca and Douglas Schmitt University
More informationAN EXPERIMENTAL INVESTIGATION OF BOILING HEAT CONVECTION WITH RADIAL FLOW IN A FRACTURE
PROCEEDINGS, Twenty-Fourth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 25-27, 1999 SGP-TR-162 AN EXPERIMENTAL INVESTIGATION OF BOILING HEAT CONVECTION
More informationReservoir properties inversion from AVO attributes
Reservoir properties inversion from AVO attributes Xin-gang Chi* and De-hua Han, University of Houston Summary A new rock physics model based inversion method is put forward where the shaly-sand mixture
More informationVelocity Measurements of Pore Fluids at Pressure and Temperature: Application to bitumen
Velocity Measurements of Pore Fluids at Pressure and Temperature: Application to bitumen Arif Rabbani 1*, Douglas R Schmitt 1, Jason Nycz 2, and Ken Gray 3 1 Institute for Geophysical Research, Department
More informationCalibration of the petro-elastic model (PEM) for 4D seismic studies in multi-mineral rocks Amini, Hamed; Alvarez, Erick Raciel
Heriot-Watt University Heriot-Watt University Research Gateway Calibration of the petro-elastic model (PEM) for 4D seismic studies in multi-mineral rocks Amini, Hamed; Alvarez, Erick Raciel DOI: 10.3997/2214-4609.20132136
More informationModeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures
Modeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures ENRU LIU, SERAFEIM VLASTOS, and XIANG-YANG LI, Edinburgh
More informationNumeric Simulation for the Seabed Deformation in the Process of Gas Hydrate Dissociated by Depressurization
Numeric Simulation for the Seabed Deformation in the Process of Gas Hydrate Dissociated by Depressurization Zhenwei Zhao 1,3 and Xinchun Shang 2 1 Department of Civil Engineering, University of Science
More informationJihoon Kim, George J. Moridis, John Edmiston, Evan S. Um, Ernest Majer. Earth Sciences Division, Lawrence Berkeley National Laboratory 24 Mar.
TOUGH+ROCMECH for the Analysis of coupled Flow, Thermal, Geomechanical and Geophysical Processes Code Description and Applications to Tight/Shale Gas Problems Jihoon Kim, George J. Moridis, John Edmiston,
More informationPermeability Estimates & Saturation Height Functions: A talk of two halves. Dr Joanne Tudge LPS Petrophysics 101 Seminar 17 th March 2016
Permeability Estimates & Saturation Height Functions: A talk of two halves Dr Joanne Tudge LPS Petrophysics 101 Seminar 17 th March 2016 Permeability: What is it? How do we measure it? Why do we need it?
More informationTemperature Dependence of Acoustic Velocities in Gas-Saturated Sandstones
Temperature Dependence of Acoustic Velocities in Gas-Saturated Sandstones Andreas Bauer, Christian Lehr, Frans Korndorffer, Arjan van der Linden Shell Global Solutions International, Rijswijk. Netherlands
More informationElements of Rock Mechanics
Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider
More informationRock physics of a gas hydrate reservoir. Gas hydrates are solids composed of a hydrogen-bonded ROUND TABLE
ROUND TABLE Rock physics of a gas hydrate reservoir JACK DVORKIN and AMOS NUR, Stanford University, California, U.S. RICHARD UDEN and TURHAN TANER, Rock Solid Images, Houston, Texas, U.S. Gas hydrates
More information16 Rainfall on a Slope
Rainfall on a Slope 16-1 16 Rainfall on a Slope 16.1 Problem Statement In this example, the stability of a generic slope is analyzed for two successive rainfall events of increasing intensity and decreasing
More informationGeology 228/378 Applied & Environmental Geophysics Lecture 8. Induced Polarization (IP) and Nuclear Magnetic Resonance (NMR)
Geology 228/378 Applied & Environmental Geophysics Lecture 8 Induced Polarization (IP) and Nuclear Magnetic Resonance (NMR) Induced Polarization (IP) and Nuclear Magnetic Resonance (NMR) 1. Time domain
More informationFlow and Transport. c(s, t)s ds,
Flow and Transport 1. The Transport Equation We shall describe the transport of a dissolved chemical by water that is traveling with uniform velocity ν through a long thin tube G with uniform cross section
More informationInvestigation of Compositional Grading in Petroleum Reservoirs
Investigation of Compositional Grading in Petroleum Reservoirs Zhangxing Chen University of Calgary Outline Importance of the Research Factors Leading to Compositional Variations Compositional Grading
More informationAPPLICATION OF 1D HYDROMECHANICAL COUPLING IN TOUGH2 TO A DEEP GEOLOGICAL REPOSITORY GLACIATION SCENARIO
PROCEEDINGS, TOUGH Symposium 2015 Lawrence Berkeley National Laboratory, Berkeley, California, September 28-30, 2015 APPLICATION OF 1D HYDROMECHANICAL COUPLING IN TOUGH2 TO A DEEP GEOLOGICAL REPOSITORY
More information6298 Stress induced azimuthally anisotropic reservoir - AVO modeling
6298 Stress induced azimuthally anisotropic reservoir - AVO modeling M. Brajanovski* (Curtin University of Technology), B. Gurevich (Curtin University of Technology), D. Nadri (CSIRO) & M. Urosevic (Curtin
More informationModeling two-phase flow in strongly heterogeneous porous media
Presented at the COMSOL Conference 2010 China COMSOL 2010 I«^rc Modeling two-phase flow in strongly heterogeneous porous media Zhaoqin Huang Research Center for Oil & Gas Flow in Reservoir, School of Petroleum
More informationUse of Seismic and EM Data for Exploration, Appraisal and Reservoir Characterization
Use of Seismic and EM Data for Exploration, Appraisal and Reservoir Characterization Anton Ziolkowski and Folke Engelmark Petroleum Geo-Services CSEG, Calgary, 6 May 2009 Outline Exploration, appraisal,
More informationAn empirical method for estimation of anisotropic parameters in clastic rocks
An empirical method for estimation of anisotropic parameters in clastic rocks YONGYI LI, Paradigm Geophysical, Calgary, Alberta, Canada Clastic sediments, particularly shale, exhibit transverse isotropic
More informationRATE OF FLUID FLOW THROUGH POROUS MEDIA
RATE OF FLUID FLOW THROUGH POROUS MEDIA Submitted by Xu Ming Xin Kiong Min Yi Kimberly Yip Juen Chen Nicole A project presented to the Singapore Mathematical Society Essay Competition 2013 1 Abstract Fluid
More informationThe role of capillary pressure curves in reservoir simulation studies.
The role of capillary pressure curves in reservoir simulation studies. M. salarieh, A. Doroudi, G.A. Sobhi and G.R. Bashiri Research Inistitute of petroleum Industry. Key words: Capillary pressure curve,
More informationDynamic permeability in poroelasticity
Stanford Exploration Project, Report 113, July 8, 2003, pages 1 454 Dynamic permeability in poroelasticity James G. Berryman 1 ABSTRACT The concept of dynamic permeability is reviewed. Modeling of seismic
More informationCalculation types: drained, undrained and fully coupled material behavior. Dr Francesca Ceccato
Calculation types: drained, undrained and fully coupled material behavior Dr Francesca Ceccato Summary Introduction Applications: Piezocone penetration (CPTU) Submerged slope Conclusions Introduction Porous
More informationLab 5 - Aquifer Elasticity and Specific Storage. Due October 12 14, 2010
Lab 5 - Aquifer Elasticity and Specific Storage Due October 12 14, 2010 The goal of this experiment is to measure the specific storage S s of a balloon, which simulates aquifer elasticity. The experiment
More informationFirst laboratory measurements of seismo-magnetic conversions in fluid-filled Fontainebleau sand
GEOPHYSICAL RESEARCH LETTERS, VOL. 33,, doi:10.1029/2005gl024582, 2006 First laboratory measurements of seismo-magnetic conversions in fluid-filled Fontainebleau sand C. Bordes, 1 L. Jouniaux, 2 M. Dietrich,
More informationFull waveform inversion of shot gathers in terms of poro-elastic parameters
Full waveform inversion of shot gathers in terms of poro-elastic parameters Louis De Barros, M. Dietrich To cite this version: Louis De Barros, M. Dietrich. Full waveform inversion of shot gathers in terms
More informationShort Note. A simple derivation of the effective stress coefficient for seismic velocities in porous rocks. Boris Gurevich. n = 1 K 0 /K s, (4)
GEOPHYSICS, VOL. 69, NO. 2 (MARCH-APRIL 2004); P. 393 397, 1 FIG. 10.1190/1.1707058 Short Note A simple derivation of the effective stress coefficient for seismic velocities in porous rocks Boris Gurevich
More informationJoint inversion of borehole electromagnetic and sonic measurements G. Gao, A. Abubakar, T. M. Habashy, Schlumberger-Doll Research
Joint inversion of borehole electromagnetic and sonic measurements G. Gao, A. Abubakar, T. M. Habashy, Schlumberger-Doll Research SUMMARY New-generation electromagnetic (EM) and sonic logging tools are
More informationOn pore fluid pressure and effective solid stress in the mixture theory of porous media
On pore fluid pressure and effective solid stress in the mixture theory of porous media I-Shih Liu Abstract In this paper we briefly review a typical example of a mixture of elastic materials, in particular,
More informationSUMMARY INTRODUCTION EXPERIMENTAL PROCEDURE
Frequency dependent elastic properties and attenuation in heavy-oil sands: comparison between measured and modeled data Agnibha Das, and Michael Batzle, Colorado School of Mines SUMMARY We have measured
More informationMicro-mechanics in Geotechnical Engineering
Micro-mechanics in Geotechnical Engineering Chung R. Song Department of Civil Engineering The University of Mississippi University, MS 38677 Fundamental Concepts Macro-behavior of a material is the average
More informationWavefields in porous media: lithology, pressure and saturation effects
Nam Hoai Pham Wavefields in porous media: lithology, pressure and saturation effects A thesis presented for the degree of Dr. Ing. (Doctor of Engineering) of the Norwegian University of science and Technology
More informationUnjacketed bulk compressibility of sandstone in laboratory experiments. R. M. Makhnenko 1 and J. F. Labuz 1
481 Unjacketed bulk compressibility of sandstone in laboratory experiments R. M. Makhnenko 1 and J. F. Labuz 1 1 Department of Civil Engineering, University of Minnesota, Minneapolis, MN 55455; PH (612)
More informationElectrical 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