Both sides previous revision
Previous revision
Next revision
|
Previous revision
|
details [2013/01/03 14:41] espen [Seismic Modeling] |
details [2013/01/04 13:36] (current) barn |
//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 areas; seismic modeling, wave equation migration velocity analysis (WEMVA) and full waveform inversion (FWI). |
==== Seismic Modeling ==== | |
The purpose of a seismic modeling (or seismic stimulation) is 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 (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 describes. In that sense, the 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 (simulating) how 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 surveys. Also 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 modelling, software 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 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]]. Once the equation is transformed into computer language, one is able to simulate a seismic wave inside a medium. | The starting point for a modeling algorithm is mathematical equations which describe seismic waves. Two 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 describes. In that sense, the 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. |
| |
The final step for the setup of a 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 a "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 parameters. These 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 === |
{{:sxx1.png?215|}}{{:sxx2.png?215|}}{{:sxx3.png?215|}} | {{: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 === |
| |
{{ :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~~ |