A kirchhoff elliptic vortex is a 2d elliptical region or patch of uniform vorticity embedded in an inviscid, incompressible and irrotational fluid. A library of matlab code for analysing post migration prestack seismic data. Array imaging using kirchhoff migration and beamforming, resolution theory for broad and narrow band array imaging in homogeneous media, topics in highfrequency, variable background imaging. Kirchhoff prestack time migration pstm can accurately image steep or overturned events, provided that turning rays and anisotropy are taken into account. A matlab version of equivalent offset migration eom john c. Finite difference modeling of acoustic waves in matlab. Either c or fortran subroutine similar to matlab functions. Equivalent offset migration is based on kirchhoff prestack migration but forms prestack migration gathers that are referred to as common scatterpoint csp gathers or. Good that some stuff is available for beginners, they r not required to rub their heads to write from scratch. In order to keep this paper short, many aspects of kirchhoff migration and modelling are ignored. Osa phaseless computational imaging with a radiating. In kirchhoff migration in practice, finitedifference migration in practice, frequencyspace migration in practice, and frequencywavenumber migration in practice, migration results of different algorithms using various parameters are compared with a desired migration.
Todays seismic processing centers need a broad portfolio. Kirchhoff and fk migration to focus ground penetrating radar. In this paper, a matlab software is presented that synthesizes 2d shot records with up to four flat reflectors with different velocities and then migrate them using. How to solve this big system of odes kirchhoffs law 2 with a trapezoidal implicit solver. There is currently no complete application that is built entirely with python i.
Kirchhoff prestack migration and modelling are linear operators. Kirchhoff migrations fast, accurate subsurface imaging and efficient velocity model building. Introduction the objective of this project is to develop a simple 2d kirchhoff migration program in java that can be used as a tutorial for others to see how a simple summationbased migration algorithm can be coded. Figure 3 illustrated that one output sample could be migrated. Writing kirchhoff migrationmodelling in a matrix form crewes. This software package is based on the wavefront construction wfc technique that represents the same robustness. Kirchhoff migration offers adaptability to irregular acquisition. Kirchhoff vortex contour dynamics simulation file exchange. Kirchhoff migration integral has two inverse integrals in an. Common offset migration using kirchhoff depth github. Although there are many methods offered by different scientists, there is not any complete migrationfocusing method that works perfectly for all scenarios. Analytical and numerical methods for solving ordinary differential equations arising in engineering applications.
Finite difference modeling of acoustic waves in matlab carrie f. The impulse response of kirchhoff migration, although it may be very long, is continuos in space. Readers need a little background knowledge on kirchhoff migration and matlab programming language to follow the procedure mentioned in this paper easier. In this thesis, the norsar3d software has been employed to trace rays.
It is particularly attractive due to its speed and targetoriented capability, which enable efficient velocity model building and updating. Alternative processing software such as promax or vista is required to read the. Reverse time migration is an ultrahighend prestack twoway waveequation migration for accurate imaging in and below areas with both great structural and velocity complexities, such as sedimentary areas with steep salt inclusions. Kirchhoff migration methods are based on the diffraction summation technique, which sums the. Downloads are available as either unix gzipped tar files or windows zip files. On your matlab prompt, type run kirmiggom to see a movie of migration as a function of maximum propagation angle amax. Travel time approximations, including higherorder terms, to extend the applicability of kirchhoff time migration. This paper aims at getting better resolution by applying two migration algorithms.
It uses central finite difference schemes to approximate derivatives to the scalar wave equation. In addition, the kirchhoff migration offers adaptability to irregular acquisition geometries, and flexibility in handling anisotropic and converted wave velocity models and surface topography. Parallel programming in fortran 90 fortran 90 is a new fortran standard that has been implemented on a few massively parallel computers and a few vector computers. Kirchhoff migration produces the section shown in figure 4. The following matlab project contains the source code and matlab examples used for kirchhoff vortex contour dynamics simulation. Simulation for 2d reverse time migration rtm, with 100 shots and 100 sensors on the free surface. Ground penetrating radar gpr based land mine detection has a main challenge of having an accurate image analysis method that is capable of reducing false alarms. The main advance of fortran 90, with respect to fortran 77, is that arrays are treated as firstclass objects.
Prestack time migration of 2d seismic records by the. The impulse response indicates that kirchhoff migration can accurately handle dips up to 90 degrees. General forms for excel is a very good software, study, work and research of large amounts of data may be stored in excel tables. Thanks for contributing an answer to stack overflow. Piratheeben kulathilagan master thesis in geosciences discipline. However this image analysis depends on having sufficient spatial resolution in the backscattered signal. When you watch the movie, notice how the migrated image changes. Note that migration can be limited to smaller dips. Range migration algorithm file exchange matlab central.
This will allow us to use the parallel framework and libraries to quickly implement new algorithms. The seg advanced modeling program seam is a partnership between industry and seg designed to advance geophysical science and technology through the construction of subsurface models and generation of synthetic data sets. This process is necessary to overcome the limitations of geophysical methods imposed by areas of. Image in and below areas with structural and velocity complexities. Feb 04, 2016 this short video explains to a nonseismologist the concept of seismic migration. This paper introduces how to optimize a practical prestack kirchhoff time migration program by the compute unified device architecture cuda on a general purpose gpu gpgpu. Prestack time migration of 2d seismic records by the kirchhoff method. Asking for help, clarification, or responding to other answers. Prestack kirchhoff time migration is applied via software written with matlab. Excel in matrix computation, and optimization algorithm of data fitting function is insufficient, combining excel and matlab are effective methods for dealing with comple. Bohlen universitat kiel wrote the matlab tool for the true amplitude. Threedimensional prestack kirchhoff migration of deep seismic. By exploiting the spatial continuity of the impulse response, an efficient implementation of kirchhoff.
We implemented the contour dynamics algorithm in matlab in order to reexamine the evolution of the kirchhoff vortex, with an emphasis on the modes of the system. Terwilliger1 abstract in this computer exercise you will modify the kirchho. A gui graphical user interface was designed to work with codes for convenience of operator. Two fitting routines are included which decompose the solutions into constituent linear eigenmodes. The controversial topic of 3d spatial sampling and its impact on the final seismic image has generated several articles in the last few years advocating coarsely sampled acquisition geometries e. Apr 15, 2016 ground penetrating radar gpr based land mine detection has a main challenge of having an accurate image analysis method that is capable of reducing false alarms. Kirchhoff and fk migration to focus ground penetrating. This short video explains to a nonseismologist the concept of seismic migration. To investigate further the quality of the migrated cube, one representative. Writing kirchhoff migrationmodelling in a matrix form. This page presents a list of programs for seismic processing, imaging and analysis. Travel time approximations, including higherorder terms. Modeling by demigration a feasibility study muhammad junaid yaqoob. Although there are many methods offered by different scientists, there is not any complete migration focusing method that works perfectly for all scenarios.
Margrave abstract a matlab toolkit, called the afd package, has been written to model waves using acoustic finite differences. However, the algorithms and the programs must still be designed carefully if good performance is to be obtained. Kirchhoff wavefront migration for efficient and accurate imaging. However, it may of course be adapted to other uses, within the limits of the license. Computing prestack kirchhoff time migration on general. Even though ground penetrating radar has been well studied and applied by many researchers for the last couple of decades, the focusing problem in the measured gpr images is still a challenging task.
The migration code is being developed in both matlab and c with careful attention to modular design. One of the main advantages of kirchhoff migration is its ability to use arbitrary input or output geometry. This article presents a method to run 3d kirchhoff prestack depth migration on. Here are the guidelines for geophysical software and algorithms papers. Search kirchhoff migration in matlab, 300 results found. The idea of using kirchhoff modeling migration operators in the explicit forms is explained. On the other hand, vermeer 1995 has advocated a fine. Subroutine kirchslow below is the best tutorial kirchhoff migrationmodeling program i could devise. P51 2d kirchhoff prestack time migration of seismic synthetic. Geophysics department of geosciences faculty of mathematics and. In this article, a new scalable 3d kirchhoff depth migration algorithm is presented on state of the art multicore cpu based cluster. A matlab version of equivalent offset migration eom.
Bohlen universitat kiel wrote the matlab tool for the trueamplitude. Parallelization of 3d kirchhoff depth migration is challenging due to its high demand of compute time, memory, storage and io along with the need of their effective management. Kirchhoff time migration seismic processing software. Im not sure of what ways youre familiar with doing this, but i was taught to do it using matrices, and ive stuck with it. The matlab program for crosscorrelating two length. First you will limit the propagation angles in the subroutine, and then you will modify the routine to handle vz. Bg group plc or any of its respective subsidiaries, affiliates and associated companies or by any of their respective officers, employees or agents makes no representation or warranty, express or implied, in respect to the quality, accuracy or usefulness of this repository. Kirchhoff algorithms kirchhoff migration schneider, 1978 is applied to the data as a convolution with a space and time variant operator. Seismic migration is the process by which seismic events are geometrically relocated in either space or time to the location the event occurred in the subsurface rather than the location that it was recorded at the surface, thereby creating a more accurate image of the subsurface. Computational imaging modalities support a simplification of the active architectures required in an imaging system and these approaches have been validated across the electromagnetic spectrum. The kirchhoff migration is the most commonly used migration for depth and time imaging. The algorithm will be described along with its assumptions and limitations.
Kirchhoff vortex contour dynamics simulation in matlab. Kirchhoff prestack depth and time migration schlumberger. This routine was an assignment for a university class and of course is highly specific. In all cases, this desired migration was obtained using the phaseshift method with appropriate parameters. A library of matlab code for analysing postmigration pre. Jan 11, 2015 simulation for 2d reverse time migration rtm, with 100 shots and 100 sensors on the free surface. Common midpoint stack, which was state of the art until the early 1970s, was able to produce seismic images. The layout shown is the one used for the kirchhoff migration algorithm described in the. Also known as kirchhoff migration, a method of seismic migration that uses the integral form kirchhoff equation of the wave equation. Solution of initial and boundary value problems, series solutions, laplace transforms, and nonlinear equations. Geomages kpstm is a full kirchhoff diffraction summation algorithm that includes advanced antialiasing filters and is capable of imaging extremely large. Migration operator aliasing and its effect on interpolation. Reverse time migration is an ultrahighend prestack twoway waveequation migration for accurate imaging in and below areas with both great structural and velocity complexities, such as. All methods of seismic migration involve the backpropagation or continuation of the seismic wavefield from the region where it was measured earths surface or along a borehole into the region to be imaged.
Kirchhoff wavefront migration for efficient and accurate imaging 3d prestack depth migrations support multiple imaging. Parallel 3d prestack depth migration using recursive. So this boils down to you being able to solve a system of four linear equations with four unknowns in matlab. An overview of all taskcentric pages can be found in the main taskcentric program list the most frequent question encountered from a new user is. Davide teixeira, alex yeh and sampath gajawada, tgs. Kirchhoff prestack time migration kpstm is an excellent tool that accurately images geological settings characterized by complex structure with relatively gentle velocity gradients. The massively parallel nature of graphics processing units has made them an attractive platform for some computationally tintensive algorithms. Migration of seismic data moves dipping events to their correct positions. A nice feature of this program is that it works ok while the. How to solve this big system of odes kirchhoffs law 2. A gui graphical user interface was designed to work with codes for convenience of.
General forms for excel is a very good software, study, work and research of large amounts of. The amplitude weighting of prestack kirchhoff migration is a very complex issue and is typically the topic of various solutions to the inversion problem. This repository includes a seismic migration routine for common offset gathers using a kirchhoff depth migration. The dip on a migration impulse response is measured as the angle.
884 760 339 1180 1383 1255 115 1003 335 60 1434 1186 739 105 1036 867 1332 1335 919 1226 893 118 1477 629 271 1017 882 767 1031 65 1207 285 934 840 1212 8 1431 301 640 568 1093