The Elastic Wave Equation Equationsequations Ofof Motionmotion

Total Page:16

File Type:pdf, Size:1020Kb

The Elastic Wave Equation Equationsequations Ofof Motionmotion TheThe ElasticElastic WaveWave EquationEquation • Elastic waves in infinite homogeneous isotropic media Numerical simulations for simple sources • Plane wave propagation in infinite media Frequency, wavenumber, wavelength • Conditions at material discontinuities Snell’s Law Reflection coefficients Free surface • Reflection Seismology: an example from the Gulf of Mexico Seismology and the Earth’s Deep Interior The elastic wave equation EquationsEquations ofof motionmotion 2 ρ∂t ui = fi + ∂ jσ ij What are the solutions to this equation? At first we look at infinite homogeneous isotropic media, then: σ ij = λθδ ij + 2µε ij σ ij = λ∂ kukδ ij + µ(∂iu j + ∂ jui ) ρ∂ 2u = f + ∂ λ∂ u δ + µ(∂ u + ∂ u ) t i i j ( k k ij i j j i ) ρ∂ 2u = f + λ∂ ∂ u + µ∂ ∂ u + µ∂ 2 u t i i i k k i j j j i Seismology and the Earth’s Deep Interior The elastic wave equation EquationsEquations ofof motionmotion –– hhomogeneousomogeneous mediamedia ρ∂ 2u = f + λ∂ ∂ u + µ∂ ∂ u + µ∂ 2u t i i i k k i j j j i We can now simplify this equation using the curl and div operators 2 2 2 2 ∇ •u = ∂iui ∇ = ∆ = ∂ x + ∂ y + ∂ z and ∆u= ∇∇ •u-∇×∇×u ρ∂ 2u= f+ (λ + 2µ)∇∇ •u- µ∇×∇×u t … this holds in any coordinate system … This equation can be further simplified, separating the wavefield into curl free and div free parts Seismology and the Earth’s Deep Interior The elastic wave equation EquationsEquations ofof motionmotion –– PP waveswaves ρ∂ 2u= (λ + 2µ)∇∇ •u- µ∇×∇×u t Let us apply the div operator to this equation, we obtain ρ∂ 2θ = (λ + 2µ)∆θ t where Acoustic wave equation θ = ∇u= ∂iui P wave velocity or 2 1 λ + 2µ ∂ θ = ∆θ α = t α 2 ρ Seismology and the Earth’s Deep Interior The elastic wave equation EquationsEquations ofof motionmotion –– sshearhear waveswaves ρ∂ 2u= (λ + 2µ)∇∇ •u- µ∇×∇×u t Let us apply the curl operator to this equation, we obtain ρ∂ 2∇×u = (λ + µ)∇×∇θ + µ∆(∇×u ) t i i we now make use of ∇×∇θ = 0 and define ∇×ui = ϕi to obtain Shear wave Wave equation for velocity shear waves 2 1 µ ∂ ϕ = ∆ϕ β = t β 2 ρ Seismology and the Earth’s Deep Interior The elastic wave equation ElastodynamicElastodynamic PotentialsPotentials Any vector may be separated into scalar and vector potentials u = ∇Φ + ∇× Ψ where P is the potential for Φ waves and Ψ the potential for shear waves ⇒ θ = ∆Φ ⇒ ϕ = ∇×u = ∇×∇× Ψ = −∆Ψ P-waves have no rotation Shear waves have no change in volume ∂ 2θ = α 2∆θ ∂ 2ϕ = β 2∆ϕ t t Seismology and the Earth’s Deep Interior The elastic wave equation SeismicSeismic VelocitiesVelocities Material and Source P-wave velocity (m/s) shear wave velocity (m/s) Water 1500 0 Loose sand 1800 500 Clay 1100-2500 Sandstone 1400-4300 Anhydrite, Gulf Coast 4100 Conglomerate 2400 Limestone 6030 3030 Granite 5640 2870 Granodiorite 4780 3100 Diorite 5780 3060 Basalt 6400 3200 Dunite 8000 4370 Gabbro 6450 3420 Seismology and the Earth’s Deep Interior The elastic wave equation SolutionsSolutions toto thethe wavewave equationequation -- ggeneraleneral Let us consider a region without sources ∂ 2η = c2∆η t Where n could be either dilatation or the vector potential and c is either P- or shear-wave velocity. The general solution to this equation is: η(xi ,t) = G(a j x j ± ct) Let us take a look at a 1-D example Seismology and the Earth’s Deep Interior The elastic wave equation SolutionsSolutions toto thethe wavewave equationequation -- hharmonicarmonic Let us consider a region without sources ∂ 2η = c2∆η t The most appropriate choice for G is of course the use of harmonic functions: ui (xi ,t) = Ai exp[ik(a j x j − ct)] Seismology and the Earth’s Deep Interior The elastic wave equation SolutionsSolutions toto thethe wavewave equationequation -- hharmonicarmonic … taking only the real part and considering only 1D we obtain u(x,t) = Acos[k(x − ct)] 2π 2π 2π k(x − ct) = kx − kct = x −ωt = x − t λ λ T c wave speed kwavenumber λ wavelength Tperiod ω frequency A amplitude Seismology and the Earth’s Deep Interior The elastic wave equation SphericalSpherical WavesWaves 2 2 Let us assume that η is a function ∂ η = c ∆η of the distance from the source t 2 1 r ∆η = ∂ 2η + ∂ η = ∂t η r r r c2 2 where we used the definition of the Laplace operator in spherical coordinates let us define η η = to obtain r 2 2 ∂ η = c ∆η with the known solution η = f (r −αt) t Seismology and the Earth’s Deep Interior The elastic wave equation GeometricalGeometrical spreadingspreading so a disturbance propagating away with spherical wavefronts decays like 1 1 η = f (r −αt) η ≈ r r r ... this is the geometrical spreading for spherical waves, the amplitude decays proportional to 1/r. If we had looked at cylindrical waves the result would have been that the waves decay as (e.g. surface waves) 1 η ≈ r Seismology and the Earth’s Deep Interior The elastic wave equation PlanePlane waveswaves ... what can we say about the direction of displacement, the polarization of seismic waves? u = ∇Φ + ∇× Ψ ⇒ u = P + S P = ∇Φ S = ∇× Ψ ... we now assume that the potentials have the well known form of plane harmonic waves Φ = Aexpi(k• x−ωt) Ψ = B expi(k• x−ωt) P = Akexpi(k• x−ωt) S = k× B expi(k• x−ωt) P waves are longitudinal as P is shear waves are transverse parallel to k because S is normal to the wave vector k Seismology and the Earth’s Deep Interior The elastic wave equation HeterogeneitiesHeterogeneities .. What happens if we have heterogeneities ? Depending on the kind of reflection part or all of the signal is reflected or transmitted. • What happens at a free surface? • Can a P wave be converted in an S wave or vice versa? • How big are the amplitudes of the reflected waves? Seismology and the Earth’s Deep Interior The elastic wave equation BoundaryBoundary ConditionsConditions ... what happens when the material parameters change? ρ1 v1 welded interface ρ v 2 2 At a material interface we require continuity of displacement and traction A special case is the free surface condition, where the surface tractions are zero. Seismology and the Earth’s Deep Interior The elastic wave equation ReflectionReflection andand TransmissionTransmission –– SSnell’snell’s LawLaw What happens at a (flat) material discontinuity? sin i1 v1 Medium 1: v1 i1 = sin i2 v2 i2 Medium 2: v2 But how much is reflected, how much transmitted? Seismology and the Earth’s Deep Interior The elastic wave equation ReflectionReflection andand TransmissionTransmission coefficientscoefficients Let’s take the most simple example: P-waves with normal incidence on a material interface R A R ρ α − ρ α = 2 2 1 1 Medium 1: r1,v1 A ρ2α2 + ρ1α1 Medium 2: r2,v2 T 2ρ α = 1 1 T A ρ2α2 + ρ1α1 At oblique angles conversions from S-P, P-S have to be considered. Seismology and the Earth’s Deep Interior The elastic wave equation ReflectionReflection andand TransmissionTransmission –– AAnsatznsatz How can we calculate the amount of energy that is transmitted or reflected at a material discontinuity? We know that in homogeneous media the displacement can be described by the corresponding potentials u = ∇Φ + ∇× Ψ in 2-D this yields ux = ∂ xΦ − ∂ y Ψz u y = ∂ z Ψx − ∂ xΨz uz = ∂ zΦ + ∂ xΨy an incoming P wave has the form ω Φ = A expi (a x −αt) 0 α j j Seismology and the Earth’s Deep Interior The elastic wave equation ReflectionReflection andand TransmissionTransmission –– AAnsatznsatz ... here ai are the components of the vector normal to the wavefront : ai=(cos e, 0, -sin e), where e is the angle between surface and ray direction, so that for the free surface Φ = A0 expik(x1 − x3 tan e − ct) + Aexpik(x1 + x3 tan e − ct) Ψ = B expik'(x1 + x3 tan f − c't) α β e f where c = c'= cose cos f Pr ω ω ω k = cose = k'= cos f α c β P SVr what we know is that σ xz = 0 σ zz = 0 Seismology and the Earth’s Deep Interior The elastic wave equation ReflectionReflection andand TransmissionTransmission –– CCoefficientsoefficients ... putting the equations for the potentials (displacements) into these equations leads to a relation between incident and reflected (transmitted) amplitudes A 4 tan e tan f − (1− tan 2 f )2 RPP = = 2 2 A0 4 tan e tan f + (1− tan f ) A 4 tan e − (1− tan 2 f ) RPS = = 2 2 A0 4 tan e tan f + (1− tan f ) These are the reflection coefficients for a plane P wave incident on a free surface, and reflected P and SV waves. Seismology and the Earth’s Deep Interior The elastic wave equation CaseCase 1:1: ReflectionsReflections atat aa freefree surfacesurface A P wave is incident at the free surface ... i j P P SV The reflected amplitudes can be described by the scattering matrix S ⎛ Pu Pd Su Pd ⎞ S = ⎜ ⎟ ⎝ Pu Sd Su Sd ⎠ Seismology and the Earth’s Deep Interior The elastic wave equation CaseCase 2:2: SHSH waveswaves For layered media SH waves are completely decoupled from P and SV waves SH There is no conversion only SH waves are reflected or transmitted ⎛ Su Sd Su Sd ⎞ S = ⎜ ⎟ ⎝ Su Sd Su Sd ⎠ Seismology and the Earth’s Deep Interior The elastic wave equation CaseCase 3:3: Solid-solidSolid-solid interfaceinterface SVr P P r SVt Pt To account for all possible reflections and transmissions we need 16 coefficients, described by a 4x4 scattering matrix.
Recommended publications
  • Glossary Physics (I-Introduction)
    1 Glossary Physics (I-introduction) - Efficiency: The percent of the work put into a machine that is converted into useful work output; = work done / energy used [-]. = eta In machines: The work output of any machine cannot exceed the work input (<=100%); in an ideal machine, where no energy is transformed into heat: work(input) = work(output), =100%. Energy: The property of a system that enables it to do work. Conservation o. E.: Energy cannot be created or destroyed; it may be transformed from one form into another, but the total amount of energy never changes. Equilibrium: The state of an object when not acted upon by a net force or net torque; an object in equilibrium may be at rest or moving at uniform velocity - not accelerating. Mechanical E.: The state of an object or system of objects for which any impressed forces cancels to zero and no acceleration occurs. Dynamic E.: Object is moving without experiencing acceleration. Static E.: Object is at rest.F Force: The influence that can cause an object to be accelerated or retarded; is always in the direction of the net force, hence a vector quantity; the four elementary forces are: Electromagnetic F.: Is an attraction or repulsion G, gravit. const.6.672E-11[Nm2/kg2] between electric charges: d, distance [m] 2 2 2 2 F = 1/(40) (q1q2/d ) [(CC/m )(Nm /C )] = [N] m,M, mass [kg] Gravitational F.: Is a mutual attraction between all masses: q, charge [As] [C] 2 2 2 2 F = GmM/d [Nm /kg kg 1/m ] = [N] 0, dielectric constant Strong F.: (nuclear force) Acts within the nuclei of atoms: 8.854E-12 [C2/Nm2] [F/m] 2 2 2 2 2 F = 1/(40) (e /d ) [(CC/m )(Nm /C )] = [N] , 3.14 [-] Weak F.: Manifests itself in special reactions among elementary e, 1.60210 E-19 [As] [C] particles, such as the reaction that occur in radioactive decay.
    [Show full text]
  • Effect of Fabric Anisotropy on the Dynamic Mechanical Behavior Of
    EFFECT OF FABRIC ANISOTROPY ON THE DYNAMIC MECHANICAL BEHAVIOR OF GRANULAR MATERIALS by Bo Li Submitted in partial fulfillment of the requirements For the degree of Doctor of Philosophy Dissertation Advisor: Prof. Xiangwu Zeng Department of Civil Engineering Case Western Reserve University January, 2011 CASE WESTERN RESERVE UNIVERSITY SCHOOL OF GRADUATE STUDIES We hereby approve the thesis/dissertation of Bo Li Candidate for Ph. D. degree*. (signed) Xiangwu Zeng Chair of Commitee Adel Saada Wojbor A. Woyczynski Xiong Yu (date) 11/15/2010 *We also certify that written approval has been obtained for any proprietary material contained therein Table of Contents Table of Contents I List of Tables IV List of Figures VI Acknowledgement XI Abstract XII Notation XIV Chapter 1 Introduction ........................................................................................ 1 1.1 Effect of Anisotropy on Granular Materials ........................................................................... 1 1.2 Research Objectives ............................................................................................................... 1 1.3 Outline of the Dissertation ...................................................................................................... 2 1.4 Review of the Effect of Fabric Anisotropy on Clay ............................................................... 3 1.4.1 Experimental Study on Effect of Fabric Anisotropy on Clay ......................................... 3 1.4.2 Analytical Method on Effect of Fabric Anisotropy
    [Show full text]
  • Confidential Manuscript Submitted to Engineering Geology
    Confidential manuscript submitted to Engineering Geology 1 Landslide monitoring using seismic refraction tomography – The 2 importance of incorporating topographic variations 3 J S Whiteley1,2, J E Chambers1, S Uhlemann1,3, J Boyd1,4, M O Cimpoiasu1,5, J L 4 Holmes1,6, C M Inauen1, A Watlet1, L R Hawley-Sibbett1,5, C Sujitapan2, R T Swift1,7 5 and J M Kendall2 6 1 British Geological Survey, Environmental Science Centre, Nicker Hill, Keyworth, 7 Nottingham, NG12 5GG, United Kingdom. 2 School of Earth Sciences, University 8 of Bristol, Wills Memorial Building, Queens Road, Bristol, BS8 1RJ, United 9 Kingdom. 3 Lawrence Berkeley National Laboratory (LBNL), Earth and 10 Environmental Sciences Area, 1 Cyclotron Road, Berkeley, CA 94720, United 11 States of America. 4 Lancaster Environment Center (LEC), Lancaster University, 12 Lancaster, LA1 4YQ, United Kingdom 5 Division of Agriculture and Environmental 13 Science, School of Bioscience, University of Nottingham, Sutton Bonington, 14 Leicestershire, LE12 5RD, United Kingdom 6 Queen’s University Belfast, School of 15 Natural and Built Environment, Stranmillis Road, Belfast, BT9 5AG, United 16 Kingdom 7 University of Liege, Applied Geophysics, Department ArGEnCo, 17 Engineering Faculty, B52, 4000 Liege, Belgium 18 19 Corresponding author: Jim Whiteley ([email protected]) 20 21 Copyright British Geological Survey © UKRI 2020/ University of Bristol 2020 22 1 Confidential manuscript submitted to Engineering Geology 23 Abstract 24 Seismic refraction tomography provides images of the elastic properties of 25 subsurface materials in landslide settings. Seismic velocities are sensitive to 26 changes in moisture content, which is a triggering factor in the initiation of many 27 landslides.
    [Show full text]
  • On the Wave Bottom Shear Stress in Shallow Depths: the Role of Wave Period and Bed Roughness
    water Article On the Wave Bottom Shear Stress in Shallow Depths: The Role of Wave Period and Bed Roughness Sara Pascolo *, Marco Petti and Silvia Bosa Dipartimento Politecnico di Ingegneria e Architettura, University of Udine, 33100 Udine, Italy; [email protected] (M.P.); [email protected] (S.B.) * Correspondence: [email protected]; Tel.: +39-0432-558-713 Received: 6 September 2018; Accepted: 25 September 2018; Published: 28 September 2018 Abstract: Lagoons and coastal semi-enclosed basins morphologically evolve depending on local waves, currents, and tidal conditions. In very shallow water depths, typical of tidal flats and mudflats, the bed shear stress due to the wind waves is a key factor governing sediment resuspension. A current line of research focuses on the distribution of wave shear stress with depth, this being a very important aspect related to the dynamic equilibrium of transitional areas. In this work a relevant contribution to this study is provided, by means of the comparison between experimental growth curves which predict the finite depth wave characteristics and the numerical results obtained by means a spectral model. In particular, the dominant role of the bottom friction dissipation is underlined, especially in the presence of irregular and heterogeneous sea beds. The effects of this energy loss on the wave field is investigated, highlighting that both the variability of the wave period and the relative bottom roughness can change the bed shear stress trend substantially. Keywords: finite depth wind waves; bottom friction dissipation; friction factor; tidal flats; bottom shear stress 1. Introduction The morphological evolution of estuarine and coastal lagoon environments is a widely discussed topic in literature, because of its complexity and the great variability of morphological patterns.
    [Show full text]
  • Evaluation of Liquefaction Potential of Soil by Down Borehole Method
    Missouri University of Science and Technology Scholars' Mine International Conferences on Recent Advances 2010 - Fifth International Conference on Recent in Geotechnical Earthquake Engineering and Advances in Geotechnical Earthquake Soil Dynamics Engineering and Soil Dynamics 26 May 2010, 4:45 pm - 6:45 pm Evaluation of Liquefaction Potential of Soil by Down Borehole Method Sumit Ghose Bengal Engineering and Science University, India Ambarish Ghosh Bengal Engineering and Science University, India Follow this and additional works at: https://scholarsmine.mst.edu/icrageesd Part of the Geotechnical Engineering Commons Recommended Citation Ghose, Sumit and Ghosh, Ambarish, "Evaluation of Liquefaction Potential of Soil by Down Borehole Method" (2010). International Conferences on Recent Advances in Geotechnical Earthquake Engineering and Soil Dynamics. 15. https://scholarsmine.mst.edu/icrageesd/05icrageesd/session01/15 This work is licensed under a Creative Commons Attribution-Noncommercial-No Derivative Works 4.0 License. This Article - Conference proceedings is brought to you for free and open access by Scholars' Mine. It has been accepted for inclusion in International Conferences on Recent Advances in Geotechnical Earthquake Engineering and Soil Dynamics by an authorized administrator of Scholars' Mine. This work is protected by U. S. Copyright Law. Unauthorized use including reproduction for redistribution requires the permission of the copyright holder. For more information, please contact [email protected]. EVALUATION OF LIQUEFACTION
    [Show full text]
  • Defining P-Wave and S-Wave Stratal Surfaces with Nine-Component Vsps
    INTERPRETER’S CORNER Coordinated by Rebecca B. Latimer Defining P-wave and S-wave stratal surfaces with nine-component VSPs BOB A. HARDAGE, MICHAEL DEANGELO, and PAUL MURRAY, Bureau of Economic Geology, Austin, Texas, U.S. Nine-component vertical seismic profile (9C VSP) data were acquired across a three-state area (Texas, Kansas, Colorado) to evaluate the relative merits of imaging Morrow and post-Morrow stratigraphy with compressional (P-wave) seismic data and shear (S-wave) seismic data. 9C VSP data generated using three orthogonal vector sources were used in this study rather than 3C data generated by only a ver- tical-displacement source because the important SH shear mode would not have been available if only the latter had been acquired. The popular SV converted mode (or C wave) utilized in 3C seismic technology is created by P-to-SV mode conversions at nonvertical angles of incidence at subsurface interfaces rather than propagating directly from an SV source as in 9C data acquisition. In contrast, 9C acquisition generates a P-wave mode and both fundamental S-wave modes (SH and SV) directly at the source station and allows these three basic compo- nents of the elastic seismic wavefield to be isolated from each other for imaging purposes. The purpose of these field tests was to compare SH, SV, and P reflectivity at targeted inter- faces. Because C waves were not an imaging objective, we used zero-offset acquisition geometry, which caused C-wave reflectivity to be quite small. We consider zero offset a source- to-receiver offset that is less than one-tenth of the depth to the shallowest receiver station.
    [Show full text]
  • Seismic Site Classifications for the St. Louis Urban Area
    Missouri University of Science and Technology Scholars' Mine Geosciences and Geological and Petroleum Geosciences and Geological and Petroleum Engineering Faculty Research & Creative Works Engineering 01 Jun 2012 Seismic Site Classifications for the St. Louis Urban Area Jaewon Chung J. David Rogers Missouri University of Science and Technology, [email protected] Follow this and additional works at: https://scholarsmine.mst.edu/geosci_geo_peteng_facwork Part of the Geology Commons Recommended Citation J. Chung and J. D. Rogers, "Seismic Site Classifications for the St. Louis Urban Area," Bulletin of the Seismological Society of America, vol. 102, no. 3, pp. 980-990, Seismological Society of America, Jun 2012. The definitive version is available at https://doi.org/10.1785/0120110275 This Article - Journal is brought to you for free and open access by Scholars' Mine. It has been accepted for inclusion in Geosciences and Geological and Petroleum Engineering Faculty Research & Creative Works by an authorized administrator of Scholars' Mine. This work is protected by U. S. Copyright Law. Unauthorized use including reproduction for redistribution requires the permission of the copyright holder. For more information, please contact [email protected]. Bulletin of the Seismological Society of America, Vol. 102, No. 3, pp. 980–990, June 2012, doi: 10.1785/0120110275 Seismic Site Classifications for the St. Louis Urban Area by Jae-won Chung and J. David Rogers Abstract Regional National Earthquake Hazards Reduction Program (NEHRP) soil class maps have become important input parameters for seismic site character- ization and hazard studies. The broad range of shallow shear-wave velocity (VS30, the average shear-wave velocity in the upper 30 m) measurements in the St.
    [Show full text]
  • Imaging Through Gas Clouds: a Case History in the Gulf of Mexico S
    Imaging Through Gas Clouds: A Case History In The Gulf Of Mexico S. Knapp1, N. Payne1, and T. Johns2, Seitel Data1, Houston, Texas: WesternGeco2, Houston, Texas Summary Results from the worlds largest 3D four component OBC seismic survey will be presented. Located in the West Cameron area, offshore Gulf of Mexico, the survey operation totaled over 1000 square kilometers and covered more than 46 OCS blocks. The area contains numerous gas invaded zones and shallow gas anomalies that disturb the image on conventional 3D seismic, which only records compressional data. Converted shear wave data allows images to be obtained that are unobstructed by the gas and/or fluids. This reduces the risk for interpretation and subsequent appraisal and development drilling in complex areas which are clearly petroleum rich. In addition, rock properties can be uniquely determined from the compressional and shear data, allowing for improved reservoir characterization and lithologic prediction. Data acquisition and processing of multicomponent 3D datasets involves both similarities and differences when compared to conventional techniques. Processing for the compressional data is the same as for a conventional OBC survey, however, asymmetric raypaths for the converted waves and the resultant effects on fold-offset-azimuth distribution, binning and velocity determination require radically different processing methodologies. These issues together with methods for determining the optimal parameters of the Vp/Vs ratio (γ) will be discussed. Finally, images from both the compressional (P and Z components summed) and converted shear waves will be shown to illuminate details that were not present on previous 3D datasets. Introduction Multicomponent 3D surveys have become one of the leading areas where state of the art technology is being applied to areas where the exploitation potential has not been fulfilled due to difficulties in obtaining an interpretable seismic image (MacLeod et al 1999, Rognoe et al 1999).
    [Show full text]
  • Shear Wave Velocity and Its Effect on Seismic Design Forces and Liquefaction Assessment
    Missouri University of Science and Technology Scholars' Mine International Conference on Case Histories in (2004) - Fifth International Conference on Case Geotechnical Engineering Histories in Geotechnical Engineering 17 Apr 2004, 10:30am - 12:30pm Shear Wave Velocity and Its Effect on Seismic Design Forces and Liquefaction Assessment Thomas G. Thomann URS Corporation, New York, New York Khaled Chowdhury URS Corporation, New York, New York Follow this and additional works at: https://scholarsmine.mst.edu/icchge Part of the Geotechnical Engineering Commons Recommended Citation Thomann, Thomas G. and Chowdhury, Khaled, "Shear Wave Velocity and Its Effect on Seismic Design Forces and Liquefaction Assessment" (2004). International Conference on Case Histories in Geotechnical Engineering. 2. https://scholarsmine.mst.edu/icchge/5icchge/session03/2 This work is licensed under a Creative Commons Attribution-Noncommercial-No Derivative Works 4.0 License. This Article - Conference proceedings is brought to you for free and open access by Scholars' Mine. It has been accepted for inclusion in International Conference on Case Histories in Geotechnical Engineering by an authorized administrator of Scholars' Mine. This work is protected by U. S. Copyright Law. Unauthorized use including reproduction for redistribution requires the permission of the copyright holder. For more information, please contact [email protected]. SHEAR WAVE VELOCITY AND ITS EFFECT ON SEISMIC DESIGN FORCES AND LIQUEFACTION ASSESSMENT Thomas G. Thomann Khaled Chowdhury URS Corporation URS Corporation New York, NY (USA) New York, NY (USA) ABSTRACT The shear wave velocity of soil and rock is one of the key components in establishing the design response spectra, and therefore the seismic design forces, for a building, bridge, or other structure.
    [Show full text]
  • Linear Elastodynamics and Waves
    Linear Elastodynamics and Waves M. Destradea, G. Saccomandib aSchool of Electrical, Electronic, and Mechanical Engineering, University College Dublin, Belfield, Dublin 4, Ireland; bDipartimento di Ingegneria Industriale, Universit`adegli Studi di Perugia, 06125 Perugia, Italy. Contents 1 Introduction 3 2 Bulk waves 6 2.1 Homogeneous waves in isotropic solids . 8 2.2 Homogeneous waves in anisotropic solids . 10 2.3 Slowness surface and wavefronts . 12 2.4 Inhomogeneous waves in isotropic solids . 13 3 Surface waves 15 3.1 Shear horizontal homogeneous surface waves . 15 3.2 Rayleigh waves . 17 3.3 Love waves . 22 3.4 Other surface waves . 25 3.5 Surface waves in anisotropic media . 25 4 Interface waves 27 4.1 Stoneley waves . 27 4.2 Slip waves . 29 4.3 Scholte waves . 32 4.4 Interface waves in anisotropic solids . 32 5 Concluding remarks 33 1 Abstract We provide a simple introduction to wave propagation in the frame- work of linear elastodynamics. We discuss bulk waves in isotropic and anisotropic linear elastic materials and we survey several families of surface and interface waves. We conclude by suggesting a list of books for a more detailed study of the topic. Keywords: linear elastodynamics, anisotropy, plane homogeneous waves, bulk waves, interface waves, surface waves. 1 Introduction In elastostatics we study the equilibria of elastic solids; when its equilib- rium is disturbed, a solid is set into motion, which constitutes the subject of elastodynamics. Early efforts in the study of elastodynamics were mainly aimed at modeling seismic wave propagation. With the advent of electron- ics, many applications have been found in the industrial world.
    [Show full text]
  • 1 Fundamental Solutions to the Wave Equation 2 the Pulsating Sphere
    1 Fundamental Solutions to the Wave Equation Physical insight in the sound generation mechanism can be gained by considering simple analytical solutions to the wave equation. One example is to consider acoustic radiation with spherical symmetry about a point ~y = fyig, which without loss of generality can be taken as the origin of coordinates. If t stands for time and ~x = fxig represent the observation point, such solutions of the wave equation, @2 ( − c2r2)φ = 0; (1) @t2 o will depend only on the r = j~x − ~yj. It is readily shown that in this case (1) can be cast in the form of a one-dimensional wave equation @2 @2 ( − c2 )(rφ) = 0: (2) @t2 o @r2 The general solution to (2) can be written as f(t − r ) g(t + r ) φ = co + co : (3) r r The functions f and g are arbitrary functions of the single variables τ = t± r , respectively. ± co They determine the pattern or the phase variation of the wave, while the factor 1=r affects only the wave magnitude and represents the spreading of the wave energy over larger surface as it propagates away from the source. The function f(t − r ) represents an outwardly co going wave propagating with the speed c . The function g(t + r ) represents an inwardly o co propagating wave propagating with the speed co. 2 The Pulsating Sphere Consider a sphere centered at the origin and having a small pulsating motion so that the equation of its surface is r = a(t) = a0 + a1(t); (4) where ja1(t)j << a0.
    [Show full text]
  • Seismic Reflection and Refraction Methods
    Seismic Reflection and Refraction Methods A. K. Chaubey National Institute of Oceanography, Dona Paula, Goa-403 004. [email protected] Introduction and radar systems. Whereas, in seismic refraction Seismic reflection and refraction is the principal method, principal portion of the wave-path is along the seismic method by which the petroleum industry explores interface between the two layers and hence hydrocarbon-trapping structures in sedimentary basins. approximately horizontal. The travel times of the Its extension to deep crustal studies began in the 1960s, refracted wave paths are then interpreted in terms of the and since the late 1970s these methods have become depths to subsurface interfaces and the speeds at which the principal techniques for detailed studies of the deep wave travels through the subsurface within each layer. crust. These methods are by far the most important For both types of paths the travel time depends upon the geophysical methods and the predominance of these physical property, called elastic parameters, of the methods over other geophysical methods is due to layered Earth and the attitudes of the beds. The various factors such as the higher accuracy, higher objective of the seismic exploration is to deduce resolution and greater penetration. Further, the information about such beds especially about their importance of the seismic methods lies in the fact that attitudes from the observed arrival times and from the data, if properly handled, yield an almost unique and variations in amplitude, frequency and wave form. unambiguous interpretation. These methods utilize the principle of elastic waves travelling with different Both the seismic techniques have specific velocities in different layer formations of the Earth.
    [Show full text]