A method of seismic migration that uses the integral form kirchhoff equation of the wave equation. A software in matlab is written which is capable of migrating common shot records. Kirchhoff migration the kirchhoff migration technique finds its origin in the field of seismics. Migration of seismic data ieno gazdag and piero sguazzero. Readers need a little background knowledge on kirchhoff migration and. As in much of processing the choice of which migration algorithm to apply is. 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. 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.
Migration of seismic data university of california, davis. Although the kirchhoff migration has been developed for the backpropagation of scalar pressure wavefields, it is often applied. In order to keep this paper short, many aspects of kirchhoff migration and modelling are ignored. System modeling key matlab commands used in this tutorial are. Cmp collection with handpicking, semblance analysis and crosscorrelation analysis. In terms of engineering calculations, it is not the most important topic that will be covered. This article presents a method to run 3d kirchhoff prestack depth migration on. Opencl is an open standard for portable, parallel programming across heterogeneous platforms. Matlab tutorial sorting es 111 16 sorting sorting lists of numbers is a commonly discussed skill covered in introductory computer programming courses.
Amplitude preserving vz prestack kirchhoff migration. One is by kirchhoffs migration using geometrical approach and other one is fk migration algorithms with. It directly provides rotary motion and, coupled with. Choose a web site to get translated content where available and see local events and offers. So this boils down to you being able to solve a system of four linear equations with four unknowns in matlab. Introduction kirchhoff prestack time migration is widely used in 3d seismic data processing while kirchhoff demigration is a useful tool to obtain a zerooffset data volume after kirchhoff prestack time migration and stack. The accurate estimation of gpr migration velocity and. In this paper, we presented how to implement and optimize prestack kirchhoff time migration algorithm.
When loops are reordered, we think differently and opportunities arise for speed improvements. Prestack time migration of 2d seismic records by the kirchhoff. Kirchhoff prestack depth and time migration schlumberger. However, sorting is used by engineers just enough to make it. On the other hand, vermeer 1995 has advocated a fine. 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.
Control tutorials for matlab and simulink motor position. Pdf kirchhoff and fk migration to focus ground penetrating radar. Rtm has a proven track record for generating final images and is increasingly used for refining structural. Matlab also allows you to write series of commands into a file and execute the file as complete unit, like writing a function and calling it. This interactive document is called a live script, and you can share it with anyone, even if they dont have matlab. Bohlen universitat kiel wrote the matlab tool for the. The kirchhoff migration is the most commonly used migration for depth and time imaging.
I am happy finally to see this simple idea in print. The magic of seismic migration tutorial for nongeophysicists. An example of synthetic data with extreme levels of noise illustrates that. Matlab symbol definitions for plotting symbol replaces an x in the plot function. Kirchhoff prestack time migration pstm can accurately image steep or overturned events, provided that turning rays and anisotropy are taken into account. Fresnel kirchoff integral in paraxial approximation. 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. Davide teixeira, alex yeh and sampath gajawada, tgs. Herman jaramillos personal website, powered by workfolio.
As the name suggests, matlab is especially designed for matrix computations. In addition, from the view of image display, the overall profile obtained by kirchhoff migration is more blurred than that obtained by using the other two algorithms. Addon toolboxes extend the matlab environment to solve problems in a range of applications, including signal processing and communications, video and image processing, control design, and computational. Diffraction stack migration then becomes a crosscorrelation of the data with an estimate of the diffraction. 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.
Introduction kirchhoff migration has been the workhorse of prestack seismic imaging for over a decade. It can be published as a pdf, word, html, or latex document. In kirchhoff migration, this is done by using the kirchhoff integral representation of. Finding indefinite integral using matlab by definition, if the derivative of a function fx is fx, then we say that an indefinite integral of fx with respect to x is fx. The algorithms are developed using matlab simulations over different. 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. Migration of seismic data moves dipping events to their correct positions, collapses. Older versions% of matlab can copy and paste entirebloc. The massively parallel nature of graphics processing units has made them an attractive platform for some computationally tintensive algorithms.
This paper aims at getting better resolution by applying two migration algorithms. An introduction to migration without a single equation by claudio. An opencl approach of prestack kirchhoff time migration. It is particularly attractive due to its speed and targetoriented capability, which enable efficient velocity model building and. Pdf ground penetrating radar gpr based land mine detection has a main challenge. In previous chapters, you have learned how to enter commands from the matlab command prompt. The idea of using kirchhoff modelingmigration operators in the explicit forms is explained.
If youre sharing your work with other matlab users, live controls are a great tool for drawing attention to important parameters. Seismic migration download ebook pdf, epub, tuebl, mobi. Arbitrary weights can be applied to the various input trace. Click download or read online button to get seismic migration book now.
According to matchedfilter theory, the amplitudes on the estimated diffraction must. Writing kirchhoff migrationmodelling in a matrix form crewes. This example is derived from gerard schusters matlab example and book seismic interferometry addpath faultmodeldata addpath migration. Implementation of kirchhoff prestack depth migration on gpu. Kirchhoff algorithms kirchhoff migration schneider, 1978 is applied to the data as a convolution with a space and time variant operator. An example of one such system was modeled, namely the diffraction and. Kirchhoff migrations fast, accurate subsurface imaging and efficient velocity model building. Lesson 17 seismic processing iris earthquake science. Matlab provides the features of a traditional programming language, as well as layout tools for designing custom graphical user interfaces. 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. Kirchhoff migration and inversion crewes research report s volume 14 2002 7 in the migration problem, the signal to be detected is the diffraction.
Kirchhoff migration schlumberger oilfield glossary. By exploiting the spatial continuity of the impulse response, an efficient implementation of kirchhoff. It allows time and depth migration methods to be incorporated within a single basic program, facilitates targetoriented migration, and enables straightforward migration velocity analysis. A visual relationship between kirchhoff migration and. Since the elements of the matrices are very small, the determinants are becoming very close to zero.
Request pdf prestack time migration of 2d seismic records by the kirchhoff. How to solve this big system of odes kirchhoffs law 2. For example if only the beige area of the previous figure were recorded then. 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. The impulse response of kirchhoff migration, although it may be very long, is continuos in space. The m files matlab allows writing two kinds of program files. This site is like a library, use search box in the widget to get ebook that you want. For example, since the derivative withrespecttox of x2 is 2x, we can say that an indefinite integral of 2x is x2. Radar imaging of concealed targets international journal of. Professional individual with highlevel knowledge of mathematics, physics, geophysics, and computer science. Matlab programming language to follow the procedure mentioned in this paper.
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. The algorithm will be described along with its assumptions and limitations. Mentioned kirchhoff migration using dsr equation is an example. Readers need a little background knowledge on kirchhoff migration and matlab programming language to follow the procedure mentioned in this paper easier. The resulting migration procedure, known as kirchhoff migration, compares favor ably with other methods. Kirchhoff migration of seismic data compressed by matching pursuit decomposition, by bin wang and keh pann, shows how to perform kirchhoff summation migration directly in the wavelet domain. Writing kirchhoff migrationmodelling in a matrix form. Matlab tutorial plotting es 111 34 to create a smoother curve for the two functions, use more x values to plug into the sine and cosine functions. In all cases, this desired migration was obtained using the phaseshift method with appropriate parameters and velocities. A kirchhoff elliptic vortex is a 2d elliptical region or patch of uniform vorticity embedded in an inviscid, incompressible and irrotational fluid. High resolution seismic imaging using least squares migration. 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. Matlab tutorial this tutorial series introduces you to the matlabr2006b suite of applications from mathworks, inc. In some cases for r and l, they are even zero, because matlab rounds numbers which are matlab tutorial matlab matrix laboratory is an interactive software system for numerical computations and graphics.
1068 626 1040 1032 324 1418 590 1346 325 625 1477 246 1117 1381 1185 968 1119 557 391 188 1007 1087 114 1477 1107 198 27 297 334 851 656 1336 1013 598 1425 873 1417 127 579 227