Marchenko imaging can be seen as an internal multiple elimination scheme in imaging domain. By solving coupled Marchenko equations, the measured seismic data can be correctly redatumed to the subsurface with source positioned at subsurface and receivers positioned at acquisition surface. Then, the artefact-free image of the source point can be retrieved by multi-dimensionally deconvolving the redatumed up- and downgoing seismic data. The Marchenko Multiple Elimination scheme, derived from coupled Marchenko equations, is a data domain multiple elimination scheme. It can successfully predict and remove all orders of internal multiple reflections without model information or adaptive subtraction. The data domain multiple elimination scheme surgically removes internal multiple reflections without touching primary reflections. In this paper, we give a detailed comparison of both imaging domain and data domain Marchenko multiple elimination schemes. One synthetic model is used to illustrate the performance of both schemes and, the advantages and disadvantages of both schemes are discussed in detail.

Seismoelectric coupling in an electric isotropic and elastic anisotropic medium is developed using a primary–secondary formulation. The anisotropy is of vertical transverse isotropic type and concerns only the poroelastic parameters. Based on our finite difference time domain algorithm, we solve the seismoelectric response to an explosive source. The seismic wavefields are computed as the primary field. The electric field is then obtained as a secondary field by solving the Poisson equation for the electric potential. To test our numerical algorithm, we compared our seismoelectric numerical results with analytical results obtained from Pride's equation. The comparison shows that the numerical solution gives a good approximation to the analytical solution. We then simulate the seismoelectric wavefields in different models. Simulated results show that four types of seismic waves are generated in anisotropic poroelastic medium. These are the fast and slow longitudinal waves and two separable transverse waves. All of these seismic waves generate coseismic electric fields in a homogenous anisotropic poroelastic medium. The tortuosity has an effect on the propagation of the slow longitudinal wave. The snapshot of the slow longitudinal wave has an oval shape when the tortuosity is anisotropic, whereas it has a circular shape when the tortuosity is isotropic. In terms of the Thomsen parameters, the radiation anisotropy of the fast longitudinal wave is more sensitive to the value of ε, while the radiation anisotropy of the transverse wave is more sensitive to the value of δ.