vkazei / fastFWI

Simple frequency domain full-waveform inversion (FWI) regularized by Sobolev space norm
MIT License
51 stars 19 forks source link

Wavefield Discrepancy between Finite Difference Modeling and Analytical Results #7

Open tjk9501 opened 6 months ago

tjk9501 commented 6 months ago

Dear Dr Vladimir Kazei

I use your published program fastFWI to run some simulations. And compare received data with 2D analytical results, modeling frequency is 90 kHz. The first model I used is a 1 meter by 1 meter homogeneous model, with transducers deployed along a circle and each transducers acted as both transmitters and receivers 1-1 The real part of received data at transducer locations obtained by finite difference modeling (FDM) and 2D analytical green functions are shown as follows 1-2 It can be seen that FDM data and analytical results matched very well.

However, when I increase the size of model to 2 meter (x direction) by 3 meter (y direction). And use parallel array setup, as shown in the following figure (grid spacing 2mm) 1-3

The real part of FDM wavefield and analytical wavefield (obtained by green functions) are shown as follows 1-4

It shows that with the increase of model size, the analytical results (homogeneous model) and FDM results begin to differ significantly. Such discrepancy influence the imaging results, leading to the algorithm failure at some circumstances. But I cannot figure out the possible reasons.

What is the possible reason that causes discrepancy between the FDM modeling wavefield and analytical results? Is such discrepancy caused by modeling grid size or some other reasons?

Hope to hear from your answers!