Implementation of the Marchenko multiple elimination algorithm

The Marchenko multiple elimination (MME) and transmission compensation schemes retrieve primary reflections in the two-way traveltime domain without model information or using adaptive subtraction. Both schemes are derived from projected Marchenko equations and are similar to each other, but they use different time-domain truncation operators. The MME scheme retrieves a new data set without internal multiple reflections. The transmission-compensated Marchenko multiple elimination scheme does the same and additionally compensates for transmission losses in the primary reflections. Both schemes can be solved with an iterative algorithm based on a Neumann series. At each iteration, a convolution or correlation between the projected focusing function and the measured reflection response is performed, and, after each convolution or correlation, a truncation in the time domain is applied. After convergence, the resulting projected focusing function is used for retrieving the transmission-compensated primary reflections and the projected Green's function is used for the physical primary reflections. We have determined that internal multiples are removed by using time-windowed input data that only contain primary reflections. We evaluate both schemes in detail and develop an iterative implementation that reproduces the presented numerical examples. The software is part of our open-source suite of programs and fits into the Seismic Unix software suite of the Colorado School of Mines.