Reverse Time Migration (RTM) MatLab Code

Abstract RTM is a depth migration method that implements a solution of the full wave equation. In this work, a methodology for the application of RTM to synthetic and real data is developed. We discuss the components of RTM: Acoustic wave modeling, boundary conditions and imaging condition. A Finite Difference Method (FDM) based on a staggered grid scheme was implemented for the wave modeling. A Hibryd ABC method that combines two of the most common solutions to the problem of boundary conditions: predictive boundary condition and absorbing boundary condition. RTM images were obtained using a Cross-correlation Imaging Condition with Wavefield Decomposition. The algorithm was applied to synthetic data obtained by modeling a complex geological structure, with a known velocity model. On real data, an image of the complex structure was obtained, which allowed us to know the real location of the reflection events. The obtained seismic image better resolved the seismic events compared with conventional methods of migration. Keywords: Migration, Imaging, RTM, Finite Difference Method, Full Wave Equation, Boundary Condition, Imaging Condition.

is distributed under the terms of the Apache 2.0 license.

Copyright 2020 Emiro Chica Quiñones - echica@terralica.com

RTM MatLab Code is distributed under the Apache License, Version 2.0 (the "License"); you may not use this file except in compliance with the License.

You may obtain a copy of the License at

http://www.apache.org/licenses/LICENSE-2.0

Unless required by applicable law or agreed to in writing, software distributed under the License is distributed on an "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. See the License for the specific language governing permissions and limitations under the License.