New Document

If you want to machinetranslate the articles into other languages, press her and then choose language on top of page

What happens when oil runs out? The technology must be used for anything else! Applications will be different with maths the same! !

Master-degree (- Ny versjon 2008) 


Under utarbeidelse.
Forum - Send gjerne rettelser
Epost
Kap.1: Innledning
Kap.2: Historie
Kap.3: Bølgeforplantning
i viskoelastiske media
Kap.4.Generering av syntetiske seismogram
Kap.5.Observasjon fra marine sedimenter
Kap.6.Dempning på syntetiske seismogram
Kap.7.Viskoelastiske dempningsmodeller
Kap.8.Dempningens innflytelse på seismiske pulser
Kap.8.b. Tidsvariant spiking filter og invers Q-filtrering
Krever passord
Kap.9.Generering av syntetiske seismogram med dempning. (krever passord)
Vannlagsfilteret
Kap.8.c. Tidsvariant spiking filter og invers Q-filtrering på Riccatiligningen
Litteraturliste

Andre artikler om seismikk
GEO 4280 Matlab excercises
Wavelet tuning - Mathematica
Wavelet tuning - and Riccatiequation
Viscoelastic damping in the Riccatiequation
Differentiation and integration in seismic
Real seismic data
Time-shifting theorem
Simple inverse problem
Damping of seismic wave
Z-transformation
Q-filter
Deconvolution
Time-varying Deconvolution
Time-varying Wiener - outline
Inverse Q-filter algorithm
Predictive Deconvolution
Predictive Deconvolution 2
Hilbert-transformation
Digitale filtre i geofysikk
frekvensanalyse av seismisk trase
Kap.10.Studium av reelle marin seismiske data
Legges ikke på nett

GEO 4280 Matlab excercises (word)
Wavelet tuning - Mathematica
 



Seismic invers Q-filtering

Spring term 2010, I will apply the theory of seismic data on an ordinary soundtrace and see the analogies. As an example: In Audacity we can apply a low pass filter on an audio recording, and we can have the same dampingeffect on a seismogram. In Audacity we must divide the trace into several windows, while in a seismic program we get a continuous solution. In a seismic inversion we try to restore the frequencies that are damped out. In Audacity we do this by applying inverse filters in each window, while we in our seismic program can do this with the entire trace at once.


Inverse Q-filtering applied on Riccatiequation

A direct application of my profession, I have in the topic inverse Q-filtering.In the fall I will hold a seminar on the subject at the University of Oslo, together with others. (The GEO-subject) I would then among other things, go through Yanguah Wang's latest book, 'Seismic inverse Q-filtering' which was published in May 2008. Inverse Q-filtering is a topic that has become more important in the seismic data, and Wang's book will undoubtedly increase interest. Damping caused by internal friction can be compensated for by deconvolution by applying a Wiener filter. But the method has its limitations - you have to share the seismic section into several windows to get good 'spiking' in that you get a non-continuous time-variation application. One should also have a signature in the minimum-phase. By inverse Q-filtering do not have such limitations while making a better 'spiking'.

Another book that I will try to teach is the book Basic Theory of Exploration Seismology By John K. Costain, Cahit Çoruh.It is now a Google book (click on the picture) so it's easy to keep track of the contents on the Internet. It has many Mathematica program that is good to work with. I will review some of these. The book is well up to the book of "Seismic inverse Q-filtering". Costain et al. has based his presentation on the Wawelet-shaping and Q-filtering with emphasis on including Hilbert-transformations and predictive Deconvolution.




Figure shows the frequency response of a seismic layer iterated with the Riccati-equation.
Figuren ved siden av til høyre viser en ferdig prosessert seksjon av sedimenter i havbunnen fra et felt i USA. Man bruker en eller flere såkalt tetthetslogger og hastighetslogger som utgangspunkt, tolker dem og lager et syntetisk seismogram som brukes til å få frem et bilde slik som til høyre.

Syntetisk seismogram

Bildet til høyre viser en seismisk seksjon som er blitt laget ved såkalt seismisk prosessering. Først er et seismogram blitt syntetisert ut fra en (eller to) logger som viser lydhastighet (eller også tetthet). Vi vet at dersom tetthet og/eller lydhastighet endrer seg får vi refleksjon, og dette kan vi lage på seismogrammet. Når vi syntetiserer kjenner vi disse data om det seismiske mediet slik at vi kan lage et foreløpig bilde av hvordan det ser ut. Da får vi med effekter som støy, multiple refleksjoner osv. Det er bl.a. naturlig at signalene også blir svakere jo dypere vi måler pga. transmisjonstap, sfærisk spredning og indre friksjon, slik at vi får et svært uklart bilde av den seismiske seksjonen etterhvert. Det er viktig å holde rede på alle de metoder som vi har brukt under syntetiseringsprosessen og kjenne dem ut og inn.



Research project
Inversion techniques for seismic data


There is a considerable interest in inversion techniques for seismic reflection data. The mathematical problem of how to obtain the properties of a reflecting medium from the properties of the reflected wave field has been studied in many papers. .

Read more here 

Inversjon

Når vi inverterer, forsøker vi å finne hvilke fysiske karakteristikker av bergarter og væsker som har produsert seismogrammet vi sitter med. Vi kan da f.eks. forsøke å få en del kunnskap ved å kartlegge nærliggende undersøkelsesområder og eventuelt samkjøre data fra andre kilder for å vite mer, men vi kjenner ikke lydhastighet og tetthet (såkalt impedans) som har skapt seismogrammet. For siden vårt foreløpige seismogram er laget fra kilder som inneholder støy og seismiske effekter som multippel-refleksjoner og dempning vet vi ikke om det gir et korrekt bilde av det reelle mediet. Likevel kan vi bruke dette til å finne impedansen ved "å gjette litt". Da må vi fjerne de effektene som vi antar har virket på seismogrammet som f.eks. støy, dempning, multiple refleksjoner osv. Og dersom vi ikke er fornøyd med resultatet har vi mulighet til å sammenligne det med nærliggende undersøkelser som vi er mer sikre på. Å invertere et seismogram er dermed det motsatte av å syntetisere et seismogram. Invertering er stort sett det samme som å gjøre et bilde skarpt (dekonvolvering). Når all prossessering er ferdig sitter vi igjen med en seismisk seksjon som på bilde til høyre.