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. In all cases, this desired migration was obtained using the phaseshift method with appropriate parameters and velocities. 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 accurate estimation of gpr migration velocity and. 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. 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. Matlab tutorial this tutorial series introduces you to the matlabr2006b suite of applications from mathworks, inc. On the other hand, vermeer 1995 has advocated a fine. The kirchhoff migration is the most commonly used migration for depth and time imaging. Kirchhoff prestack depth and time migration schlumberger. The algorithm will be described along with its assumptions and limitations. Kirchhoff prestack time migration pstm can accurately image steep or overturned events, provided that turning rays and anisotropy are taken into account. The magic of seismic migration tutorial for nongeophysicists. The algorithms are developed using matlab simulations over different.
Click download or read online button to get seismic migration book now. A kirchhoff elliptic vortex is a 2d elliptical region or patch of uniform vorticity embedded in an inviscid, incompressible and irrotational fluid. This article presents a method to run 3d kirchhoff prestack depth migration on. Request pdf prestack time migration of 2d seismic records by the kirchhoff. Choose a web site to get translated content where available and see local events and offers. Writing kirchhoff migrationmodelling in a matrix form. Matlab provides the features of a traditional programming language, as well as layout tools for designing custom graphical user interfaces. Fresnel kirchoff integral in paraxial approximation. 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. In terms of engineering calculations, it is not the most important topic that will be covered. 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. An opencl approach of prestack kirchhoff time migration. Prestack time migration of 2d seismic records by the kirchhoff. Pdf ground penetrating radar gpr based land mine detection has a main challenge.
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. 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. 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. In previous chapters, you have learned how to enter commands from the matlab command prompt. 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 massively parallel nature of graphics processing units has made them an attractive platform for some computationally tintensive algorithms. Threedimensional prestack kirchhoff migration of deep seismic. 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. Opencl is an open standard for portable, parallel programming across heterogeneous platforms. Migration of seismic data ieno gazdag and piero sguazzero. By exploiting the spatial continuity of the impulse response, an efficient implementation of kirchhoff. Diffraction stack migration then becomes a crosscorrelation of the data with an estimate of the diffraction. One is by kirchhoffs migration using geometrical approach and other one is fk migration algorithms with.
In this paper, we presented how to implement and optimize prestack kirchhoff time migration algorithm. Radar imaging of concealed targets international journal of. This interactive document is called a live script, and you can share it with anyone, even if they dont have matlab. 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. It directly provides rotary motion and, coupled with. This site is like a library, use search box in the widget to get ebook that you want. High resolution seismic imaging using least squares migration.
An example of one such system was modeled, namely the diffraction and. Writing kirchhoff migrationmodelling in a matrix form crewes. According to matchedfilter theory, the amplitudes on the estimated diffraction must. 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. As in much of processing the choice of which migration algorithm to apply is. Herman jaramillos personal website, powered by workfolio. The resulting migration procedure, known as kirchhoff migration, compares favor ably with other methods. An example of synthetic data with extreme levels of noise illustrates that. This example is derived from gerard schusters matlab example and book seismic interferometry addpath faultmodeldata addpath migration. Kirchhoff migration the kirchhoff migration technique finds its origin in the field of seismics. How to solve this big system of odes kirchhoffs law 2. Matlab symbol definitions for plotting symbol replaces an x in the plot function.
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 software in matlab is written which is capable of migrating common shot records. In kirchhoff migration, this is done by using the kirchhoff integral representation of. It is particularly attractive due to its speed and targetoriented capability, which enable efficient velocity model building and. In order to keep this paper short, many aspects of kirchhoff migration and modelling are ignored. Lesson 17 seismic processing iris earthquake science. Readers need a little background knowledge on kirchhoff migration and matlab programming language to follow the procedure mentioned in this paper easier.
For example, since the derivative withrespecttox of x2 is 2x, we can say that an indefinite integral of 2x is x2. The idea of using kirchhoff modelingmigration operators in the explicit forms is explained. Davide teixeira, alex yeh and sampath gajawada, tgs. 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. 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. The impulse response of kirchhoff migration, although it may be very long, is continuos in space. 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 programming language to follow the procedure mentioned in this paper. Migration of seismic data moves dipping events to their correct positions, collapses. Kirchhoff algorithms kirchhoff migration schneider, 1978 is applied to the data as a convolution with a space and time variant operator. Professional individual with highlevel knowledge of mathematics, physics, geophysics, and computer science. Cmp collection with handpicking, semblance analysis and crosscorrelation analysis. Matlab tutorial sorting es 111 16 sorting sorting lists of numbers is a commonly discussed skill covered in introductory computer programming courses. Seismic migration download ebook pdf, epub, tuebl, mobi. System modeling key matlab commands used in this tutorial are. Rtm has a proven track record for generating final images and is increasingly used for refining structural. Amplitude preserving vz prestack kirchhoff migration. Mentioned kirchhoff migration using dsr equation is an example.
Kirchhoff migrations fast, accurate subsurface imaging and efficient velocity model building. Implementation of kirchhoff prestack depth migration on gpu. A visual relationship between kirchhoff migration and. Although the kirchhoff migration has been developed for the backpropagation of scalar pressure wavefields, it is often applied. Readers need a little background knowledge on kirchhoff migration and. Older versions% of matlab can copy and paste entirebloc. An introduction to migration without a single equation by claudio. As the name suggests, matlab is especially designed for matrix computations. If youre sharing your work with other matlab users, live controls are a great tool for drawing attention to important parameters. Introduction kirchhoff migration has been the workhorse of prestack seismic imaging for over a decade. Since the elements of the matrices are very small, the determinants are becoming very close to zero. Pdf kirchhoff and fk migration to focus ground penetrating radar. 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. For example if only the beige area of the previous figure were recorded then.
The m files matlab allows writing two kinds of program files. However, sorting is used by engineers just enough to make it. Bohlen universitat kiel wrote the matlab tool for the. Kirchhoff migration and inversion crewes research report s volume 14 2002 7 in the migration problem, the signal to be detected is the diffraction.
Arbitrary weights can be applied to the various input trace. I am happy finally to see this simple idea in print. This paper aims at getting better resolution by applying two migration algorithms. 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. Migration of seismic data university of california, davis. So this boils down to you being able to solve a system of four linear equations with four unknowns in matlab.
It can be published as a pdf, word, html, or latex document. When loops are reordered, we think differently and opportunities arise for speed improvements. A method of seismic migration that uses the integral form kirchhoff equation of the wave equation. Control tutorials for matlab and simulink motor position. Kirchhoff migration schlumberger oilfield glossary.
547 1503 1353 694 1239 1441 672 883 1597 1127 1357 1361 633 1424 374 1386 1388 1449 1480 170 484 1305 713 979 1019 1008 573 1314 1495 355 579 763 1030 86 594 1314 1315 1099 1043 297 416 406 776 268 739 1122 288