Im given an initial distribution of the wall of water, shall we say, and its velocity, the normal thing. This equation is simply the elastodynamic wave equation fourier transformed over space and time. The seismic wave equation x 1 x 2 x 3 t x 1 t x 1 dx 1 dx 2 dx 3 figure 3. The algorithms based on the elastic wave equation section l. Elastic wave equation migration for isotropic and hti media richard a. Bostock and jeanphilippe mercier department of earth and ocean sciences, university of british columbia, 6339 stores road, vancouver, bc, canada v6t 1z4. Teleseismic waveform modelling with a oneway wave equation. After wavefield separation we extrapolate p and s wave scalar potentials using a shot profile wave equation migration operator with an angle gather imaging condition.
Introduction to the theory of seismic wave propagation indico. First, a linear relationship between the traveltime residual 1tdtobs tsyn and the relative velocity perturbation c. Wave equation datuming applied to marine obs data and to land. Elastic wave equation university of calgary in alberta. At a sufficiently great depth, the upcoming wave u is taken to be 0. The wave equation is a partial di erential equation that relates second time and spatial derivatives of propagating wave disturbances in a simple way. Since then, surface waves have been used to treat many specific problems. In this paper, we propose a waveequationbased traveltime seismic tomography method with a detailed description of its stepbystep process. The full wave finitedifference method is used to calculate the waves from the source and receiver to the subsurface. Velocity tends to increase with depth through earths crust and mantle, but drops sharply going from the mantle to the outer core. Matrix transformation and solutions of wave equation of. A robust and efficient iterative method for the numerical solution of the helmholtz equation phd thesis, delft, 2005. Seismology and the earths deep interior the elastic wave equation solutions to the wave equation solutions to the wave equation ggeneraleneral let us consider a region without sources.
For p wave seismology it is common to approximate wave propagation using the scalar wave equation. This approximation is often sufficient to produce reasonable p wave structural images. The threedimensional current and surface wave equations. Wave equation based seismic illumination analysis xiaobi xie1, shengwen jin2, and rushan wu1 abstract we present a wave equation based method for seismic illumination analysis. The propagation velocity of seismic waves depends on density and elasticity of the medium as well as the type of wave. On the other hand, model b has a pwave velocity of v p 4522 m s. We studied two models a and b of an elastic homogeneous and isotropic medium. The standard form for seismic wave equation in homogeneous media is. Displacements occurring from a harmonic plane p wave top and s wave bottom traveling horizontally across the page. Waveequationbased traveltime seismic tomography part. Claerbout 1970 and methods that factorize the wave equation e. In principle, what we need is a formulation of the seismic source, equations to describe elastic wave propagation once motion has started somewhere, and a theory for coupling the source description to the solution for the equations of motion. The acoustic, sourcefree, wave equation is given by the helmholtz equation and is written as neglecting density. So the solution is 12 of a delta function thats traveling.
This type of wave forces the point in earth from where it passes to vibrate back and forth in the direction in which the wave is travelling. Basic principles of the seismic method tu delft ocw. And those waves are 12 of a delta function each way. The velocities of these waves are given in terms of the bulk and shear modulii, k, and density. A oneway wave equation based, generalized screen propagator is used to extrapolate the wave. Mulder pdf, bibtex a robust iterative solver for the twoway wave equation based on a complex shiftedlaplace operator. Seismic wave propagation in stratified media has been studied for a long time aki and. Raypaths are lines that show the direction that the seismic. The christoffel equation takes the form of a simple eigenvalueeigenvector problem, as follows. Seismology and the earths deep interior elasticity and seismic waves linearlinear and nonlinear stress and strainand nonlinear stress and strain stress vs.
Elastic wave propagation and generation in seismology. The seismic wave equation in isotropic media if velocity is only a function of depth, then the material can be modeled as a series of homogeneous layers. For a nondispersive system where all frequencies of excitation. Introduction the very modern technology nonlinear realtime expert seismology is the main and best tool which can be used by the hydrocarbon industry to map petroleum and gas deposits in the earth. Waveequationbased seismic illumination analysis xiaobi xie1, shengwen jin2, and rushan wu1 abstract we present a waveequationbased method for seismic illumination analysis.
A factor defining an exponential decrease with frequency f and propagation time t of a seismic body wave. Seismic waves are waves of energy that travel through the earth, and are a result of an earthquake, explosion, or a volcano. For a nondispersive system where all frequencies of excitation propagate at the same velocity, the formula for sinusoidal or harmonic waves. In the prestack domain, wed is applied in two steps. Solving the eigensystem in equation is straightforward. Examples of this would include many applications of ray theory which do. Seismic wave fields are recorded by a seismometer, hydrophone in water, or accelerometer. Crustal attenuation in the region of the maltese islands. Regional pn wave teleseismic p wave and s wave propagation surface wave propagation at different scales psvrayleigh wave propagation plane sh love wave propagation plane technologies for imaging the moho. Animations are provided to view wave propagation in a 3dimensional solid for compressional p, shear s, rayleigh r and love l waves. Pdf in this paper, we propose a waveequationbased traveltime seismic tomography method with a detailed description of its stepbystep. I want to solve the wave equation, study it further for other initial conditions.
Direct wave reflected wave refracted wave geophones laid out in a line to record arrivals from a shot. Critical refraction is also recorded at distance beyond which angle of. Seismic phases and 3d seismic waves main seismic phases. The elastic wave equation equation of motion the force balance equation can be written as. It specifies the propagation velocity and particlemotion also called polarization direction for each plane wave component in the fourier domain. Equation of motion p and s waves plane wave solutions. The wave equation based approach can properly handle forward multiplescattering phenomena, including focusingdefocusing, diffraction, and interference effects. Direct wave attenuation, intrinsic attenuation, scattering attenuation, medium heterogeneity, coda normalization method. Examples of such signals bold lines and corresponding hilbert. In reflection seismology, the anelastic attenuation factor, often expressed as seismic quality factor or q which is inversely proportional to attenuation factor, quantifies the effects of anelastic attenuation on the seismic wavelet caused by fluid movement and grain boundary friction. Waveequationbased traveltime seismic tomography part 1. We simplify it to the standard form by modeling the material as series of homogeneous layers.
Pdf waveequationbased traveltime seismic tomography part. As a seismic wave propagates through a medium, the elastic energy associated with the wave is gradually. Modeling based on oneway wave equations does not include multiples, while modeling based on twoway wave equations includes. In the earth, seismic waves can propagate as longitudinal p or as shear s waves. We discussed two types of waves pwavescompressional swavesshear finally, if we assume no shearing then we reduced it to an acoustic wave equation. Pdf we present a waveequationbased method for seismic illumination analysis. Chapter 2 the seismic wave equation the classical physics behind most of seismology starts with newtons second law of motion f ma 2. Numerical simulation of seismic wave propagation produced. Rochester2 2 1 york university, earth sciences, downsview, ontario m3j 1p3 canada department of physics, memorial university of newfoundland, st.
Berryhill, 1979, berryhill, 1984 can provide a basis computation to deal with irregular surfaces wiggins, 1984 and variable velocities. For example, for a moderate number of threecomponent seismograms, 0. The p waves are the fastest moving waves and are simply sound waves. It is integrated upward with 1 112 where the product of the reflection coefficient c and the downgoing wave d act as a source for the upgoing wave. Separation of the elastodynamic equation in a homogeneous isotropic.
Engineering, helmholtz differential equation, wave equation, real time expert system. In these cases, more complex methods, such as wave equation datuming wed, are adopted. A display in the y, t plane is called a seismic section. The kirchhoff integral solution to the scalar wave equation using both nearfield and farfield terms. A display in the f, t plane is called a common midpoint gather or a common reflection point gather, which. Among the many types of seismic waves, one can make a broad distinction between body waves, which travel through the earth, and surface waves, which travel at the earths surface 4850. Seismic wave animations are used to illustrate different types of wave propagation through elastic materials. The seismic wave equation rick aster february 15, 2011 waves in one dimension. To generate images m pp and m ps simultaneously, we can write. Sauter we now consider the wave equation, where, for simplicity, we restrict ourselves to a homogeneous righthand side.
So for the wave equation, what comes out of a delta function in 1d. Numerical simulation of seismic wave propagation produced by. Wave equation datuming applied to marine obs data and to. Mar 24, 2018 the nature of general solutions to the wave equation are discussed along with the significance of a constant that emerges between the second spatial and temporal derivatives.
Recording at each geophone is a waveform called a seismogram. In the last section, we transform the the generalized wave equation and. And of course, initial conditions plural because the wave equation is second order in time. The derivations of oneway or parabolic wave equations can be split into two categories. To contrast the developments in this paper with conventional logic and to simplify discussion, i will temporarily address deep water kh k1 propagating waves such that hx, t, c 5 ax, t cosc, where c5 kx 2 st.
Seismology is a datadriven science and its most important discoveries usually result from analysis of new data. Pdf waveequationbased seismic illumination analysis. A oneway waveequationbased, generalized screen propagator is used to extrapolate the wave. Nonlinear seismic wave motion in elastodynamics with. A oneway wave equation based, generalized screen propagator is used to extrapolate the wave fields from sources and. Furthermore, the slow wind, tide, and densitydriven horizontal velocities are often stipulated to be independent of z a priori. Yang university of california summary we propose a full wave equation based method for seismic illumination analysis. The p wave is a longitudinal wave made up of a series of compressions and rarefactions. A oneway waveequationbased, generalized screen propagator is used. Indeed, when the medium is smoothly heterogeneous and does not contain scales smaller than the minimum wavelength of the wave.
The wave equation to motivate our discussion, consider the onedimensional wave equation. The relative amplitudes of the transmitted and reflected waves depend on. A fullwave equation based seismic illumination analysis. We will see that the fundamental theory of wave propagation is primarily based on two equations. Matrix transformation and solutions of wave equation of free. Seismic scalar wave equation modeling by a convolutional differentiator article pdf available in bulletin of the seismological society of america 821.
Abstract prestack wave equation migration of isotropic or anisotropic elastic seismic data is described as vector wavefield extrapolation, plus an imaging condition for combinations of shot and receiver wave modes. For free space, the onedimensional wave equation is derived. Within each layer, there are no gradients in the lames parameters and so these terms go to zero. Generalized to a continuous medium1, this equation becomes. The wave phenomena occurring at a boundary between two layers are discussed, such as snells. Elastic waveequation migration for isotropic and hti media. Elastic wave equation has been widely used to describe wave propagation in an elastic medium, such as seismic waves in earth and ultrasonic waves in human body. Direct signal from shot travels along top of first layer. The seismic wave equation the classical physics behind most of seismology starts with newtons second law of motion f ma 2. S wave propagation is pure shear with no volume change, whereas pwaves involve both a volume change and shearing change in shape in the material. Well, a wave goes to the right, and a wave goes to the left. For wave propagation problems assuming linear elasticity is usually sufficient. The wave equation is a partial differential equation that relates second. Seismic reflection imaging acoustic pressure source is set off near the surface sound waves propagate in all directions from the source 0.
Basic principles of the seismic method in this chapter we introduce the basic notion of seismic waves. Seismic rays are used instead of the wave front to describe the wave propagation. I see that let me write down the other half thats traveling the other way. Least squares migration of converted wave seismic data. Every point on the wave front is a source of a new wave that travels out of it in the form of spherical shells. Multiplying f0 1 by f 00 1 results in a function f that is identical to the solution of the generalized wave equation. Least squares migration of converted wave seismic data cseg. A fullwave equation based seismic illumination analysis method. A similar equation should be written for each atom in the lattice, resulting in n coupled differential equations, which should be solved simultaneously n total. The derivative matrix defines the direction of plane wave propagation. We present a wave equation based method for seismic illumination analysis. Seismology and the earths deep interior elasticity and seismic waves stressstrain regimesstressstrain regimes linear elasticity teleseismic waves rupture b, reaknig stable slip aseismic. In both cases, factorization involves choosing a preferred axis or direction of propagation followed by splitting the. A full wave equation based seismic illumination analysis method x.
797 1058 139 544 1430 888 1049 591 459 1360 374 93 930 1085 972 206 561 322 1282 899 256 1276 1421 1241 372 362 397 116 709 1085 712 940 50 1354 806 32 224 438 193 940 731 1182 359