Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
details [2013/01/03 13:57]
espen [Seismic Modeling]
details [2013/01/04 13:36] (current)
barn
Line 2: Line 2:
 //This page is under development// //This page is under development//
  
-Our main research areas are seismic modeling, wave equation migration velocity analysis (WEMVA) and full waveform inversion (FWI). In this page we give short introductions to the different areas. To make the text easy to read, and not require too much background knowledge, the most technical details are leaved out. +Our main research areas are seismic modeling, wave equation migration velocity analysis (WEMVA) and full waveform inversion (FWI). Below we give short introductions to our main research areasseismic modeling, wave equation migration velocity analysis (WEMVA) and full waveform inversion (FWI).
-==== Seismic Modeling ==== +
-The purpose of a seismic modeling (or seismic stimulationis a process where computers are used to simulate how seismic waves propagate through the Earth. The results from the simulations are used to understand how seismic waves in a real experience behave and propagate. Other motivations for performing seismic modeling are for instance the fact that seismic modeling is a major part of several algorithms which purpose is to give detailed images of the subsurface based on real data collected from for instance a seismic survey.+
  
-The starting point for the modeling is to find a mathematical equation which describes the wave phenomenon in interest. Two standard equations are the acoustic wave equation (AWEand the elastic wave equation (EWE). The major difference between these two equations is that the latter includes [[http://en.wikipedia.org/wiki/S-wave|shear waves]] in addition to [[http://en.wikipedia.org/wiki/P-wave|pressure waves]] which the former equation describes. In that sensethe latter equation is physically more correct since it includes more true physics. On the other hand, the EWE is more complex than the AWE, and thus needs more computer resources.+==== Seismic Modelling ==== 
 +Seismic modeling is the process of modeling (simulatinghow seismic waves propagate through the earth 
 +using a computer 
 +The results from the simulations are used to understand how seismic waves behave and propagate and a typical practical application is designing and optimising new seismic surveysAlso seismic modelling is an integral part of seismic software used to obtain images and mechanical parameters of the subsurface. In order to actually perform seismic modellingsoftware based on a modeling algorithm is needed in addition to computer hardware.
  
-The next point is to transform the mathematical equation into a "language" which the computer understandThere exist several useful methods to do this, where the two most popular maybe are the [[http://en.wikipedia.org/wiki/Finite_difference_method|finite difference method]] and the [[http://en.wikipedia.org/wiki/Finite_element_method|finite element method]]. Once the equation is transformed into computer languageone is able to simulate a seismic wave inside a medium.+The starting point for a modeling algorithm is mathematical equations which describe seismic wavesTwo standard equations are the acoustic wave equation (AWE) and the elastic wave equation (EWE). The major difference between these two equations is that the latter includes [[http://en.wikipedia.org/wiki/S-wave|shear waves]] in addition to [[http://en.wikipedia.org/wiki/P-wave|pressure waves]] which the former equation describesIn that sense, the latter equation is physically more correct since it includes more true physics. On the other handthe EWE is more complex than the AWE, and thus needs more computer resources.
  
-The final step for the setup of seismic modeling is the construction of a model of the medium one is interested in. This model is used as input for the modeling software.+The mathematical equation must be translated into "language" which the computer understand. There exist several useful methods to do this, where the two most popular maybe are the [[http://en.wikipedia.org/wiki/Finite_difference_method|finite difference method]] and the [[http://en.wikipedia.org/wiki/Finite_element_method|finite element method]]. In addition to a wave equation understandable by a computer, a model of the earth subsurface is also needed in the form of  
 +mechanical parametersThese parameters are used as input data to the modelling software.
  
-An important observation of seismic modeling is that it is impossible to simulate wave propagation phenomena which are identical to the one observed in nature. This fact is due to the principle that when performing modeling several //approximations// are necessary; i.e the mathematical equations involved in the modeling are not able to describe //all// physical phenomena occurring in a real wave propagation. In addition, when transforming the equation to computer language several approximations are used. However, the most important wave phenomena are described by the equation and the software, and are not too far away from the one observed in nature.+It is impossible to simulate wave propagation phenomena which are completely identical to the real ones observed in nature. This fact is due to the principle that when performing modelling several //approximations// are necessary; i.e the mathematical equations involved in the modeling are not able to describe //all// physical phenomena occurring in a real wave propagation. In addition, when transforming the equation to computer language several approximations are used. However, the most important wave phenomena are described by the equation and the software, and are not too far away from the one observed in nature.
  
 === Example === === Example ===
 +Consider a two-dimensional medium with a box which is placed in the middle of the model, see the image below. The density, P- and S-wave velocities are used as physical model in the simulation. Below is an image of the P-wave velocities. The topology of the density and S-wave velocity are the same as the P-wave velocity.
 +{{ :vp.png?300 |The P-wave velocities used in the modeling}} 
  
 +Now, a seismic wave is starting to propagate from the left in the model. The wave propagation at different times after the start are given in the images below.
  
 +{{:sxx1.png?215|}}{{:sxx2.png?215|}}{{:sxx3.png?215|}}
 ==== Wave Equation Migration Velocity Analysis (WEMVA) ==== ==== Wave Equation Migration Velocity Analysis (WEMVA) ====
 +The ultimate objective of seismic surveys is to compute mechanical parameters like bulk modulus, shear modulus and density of the subsurface using seismic data observed at the surface. If we imagine that the subsurface is divided into small regular cells of size, say, 10 meters in all three directions, the seismic //inverse problem// then consists of estimating an average mechanical quantity for each cell. Usually seismic surveys covers an area of several thousand square kilometers to a depth of about 5-10 kilometer, implying that  a very large number of unknown parameters need to be estimated. 
 +
 +Wave equation migration velocity analysis is one of the approaches used to estimate primarily the bulk modulus, or equivalently the wave velocity for pressure-waves. The main idea for performing this is to
 +create an image of the subsurface by focusing the reflected seismic waves. The focusing is entirely artificial and performed by software on a computer system. It turns out that the sharpness of the image depends on the wave velocity. By systematically changing the wave velocity until a sharp image is found, one obtains an estimate of the correct wave velocity. Focusing is mainly sensitive to the average of the seismic velocities over many cells, which implies that the estimated velocities are a smoothed version of the true velocity as can be seen on the example below.
 +
 ==== Full Waveform Inversion (FWI) ==== ==== Full Waveform Inversion (FWI) ====
 === Example === === Example ===
Line 30: Line 40:
  
 {{ :fwi-optimvp.png?360 | Best velocity model using FWI with WEMVA as initial model}} {{ :fwi-optimvp.png?360 | Best velocity model using FWI with WEMVA as initial model}}
- 
 ~~DISCUSSION:off~~ ~~DISCUSSION:off~~
details.1357221438.txt.gz · Last modified: 2013/01/03 13:57 by espen
CC Attribution-Share Alike 4.0 International
www.chimeric.de Valid CSS Driven by DokuWiki do yourself a favour and use a real browser - get firefox!! Recent changes RSS feed Valid XHTML 1.0