Selecting aperture for RTM

[kerim371]

Hi,

I get RTM like in the picture below.
I suspect I need to set an aperture for migration.
Is it possible? Or the problem isn't connected with aperture?

# Setup operators
Pr = judiProjection(d_obs.geometry)
F = judiModeling(model0; options=jopt)
Ps = judiProjection(q.geometry)
J = judiJacobian(Pr*F*adjoint(Ps), q)

# RTM
rtm = adjoint(J)*d_obs

[mloubout]

What are your options jopt ? Looks like usual FWI artifacts that can be removed with isic or a Laplacian filter

[kerim371]

RTM options:

# JUDI options
jopt = JUDI.Options(
    save_data_to_disk=true, 
    file_path=shots_out_dir, 
    file_name=shots_out_name,
    limit_m = true,
    buffer_size = 100f0,
    optimal_checkpointing=false)

FWI options:

jopt = JUDI.Options(
    save_data_to_disk=true, 
    file_path=shots_out_dir, 
    file_name=shots_out_name,
    limit_m = true,
    buffer_size = 100f0,
    optimal_checkpointing=false)

Buffersize is such small due to RAM limitations.

Model computed with FWI

[mloubout]

You may wanna try isic=true for RTM

[kerim371]

Thank you!
Now it looks much better.

What is possible the reason of random noise (especially in the upper part)?

[mloubout]

Sweet