slimgroup / JUDI.jl

Julia Devito inversion.
https://slimgroup.github.io/JUDI.jl
MIT License
94 stars 29 forks source link

Describe Wavefield Attenuation Not At Absorbing PML Regions #219

Closed tjk9501 closed 6 months ago

tjk9501 commented 6 months ago

Hello,

In simulations I input some laboratory sampled signals as source data. However, even if I carefully deploy the transducers to make sure they do not located at predefined PML boundary region. The amplitude of sampled signals at receivers seems significantly reduced compared to source data. Here I presented the waveform of source data, which is a 16 cycled sinusodial pulse. SourcePulse The received data, which propagated in a homogeneous medium several meters away from the source, lookes like this DevitoFWISimul If no attenuation is presented, I assume the pulse might correspond to the analytical expression of Hann windowed sinusodial pulse, and it looks like HannWindowedSimul

Why is attenuation occured? I wonder what causes the amplitude difference between finite difference modeling and analytical expression of pulse, and how to compensate this effect in the analytical expression?

Thanks!

mloubout commented 6 months ago

There is no attenuation here, this is the expected propagation effect which is in 2D the convolution of your source with a Hankel function. The energy is preserved in the isotropic acoustic case which means the energy of "25" of your source at a single point is distributed over all space after propagation.

You can see a more detailed overview of 2D propagation effects for a ricker wavelt here:

https://nbviewer.org/github/devitocodes/devito/blob/master/examples/seismic/acoustic/accuracy.ipynb

A simple rule of thumb is that the amplitude follow spherical spreading for a point source which means at time t the amplitude is reduced by 1/sqrt(t) in 2D and 1/t in 3D

tjk9501 commented 6 months ago

Thank you for your generous and detailed explanation!