Amplitude extraction

[fjsouzageo]

Hello! I've been working for a few days with the wedge model code (1412_Tuning_and_AVO/tuning_wedge_v2.py). Thanks for sharing it. I'm not a programmer and it is really helping me.

I would like to share a question and hear your opinion about it: in the algorithm, the graph "thickness vs amplitude" is constructed by extracting the amplitude over the upper time of the modeled wedge layer, that is, a time-slice extraction. However, as the algorithm itself shows, sometimes the seismic reflector mismatches the modeled interface.

So, my question is: wouldn't it be more appropriate to make te code extract the maximum amplitude of the synthetic reflector? By doing this we would be simulating how the amplitude would actually be if the reflector was mapped.

To illustrate, I made this change in the code and generated two figures, comparing this graph to the scenarios:

First, as the code is now:

And with the suggested amplitude extraction to the maximum:

Thanks!