Acoustic multi-parameter full waveform inversion based on the wavelet method

Figures & data

Figure 1. The snapshot of wavefield at propagation time $1.1\mathrm{s}$ by the wavelet method (a) without absorbing boundary conditions and (b) with PML absorbing boundary conditions.

Figure 2. The sketch map of data acquisition system. The blue points denote the shots and the red crosses denote the receivers. The shots and receivers are arranged under the surface $24\mathrm{m}$. The thickness of PML is $420\mathrm{m}$.

Figure 3. The shot gather for the Marmousi model computed by the wavelet method. The shot located at two different positions near the surface, i.e. (a) the left side of the model, x = 0; (b) the middle of the model, x = 3 km. Each shot has 200 receivers. The horizontal axis is the receiver points. The vertical axis is the time points.

Figure 4. The shot gather for the Marmousi model computed by the staggered-grid method. The shot located at two different positions near the surface, i.e. (a) the left side of the model, x = 0; (b) the middle of the model, $x=3\mathrm{km}$. Each shot has 200 receivers. The horizontal axis is the receiver points. The vertical axis is the time points.

Figure 5. The Marmousi exact model. (a) velocity, (b) density.

Figure 6. The initial model for FWI. (a) velocity, (b) density.

Figure 7. The inversion result with the 0–5 Hz data. (a) velocity, (b) density.

Figure 8. The inversion result with the 0–15 Hz data. (a) velocity, (b) density.

Figure 9. The inversion result with the 0–25 Hz data. (a) velocity, (b) density.

Figure 10. The inversion result with the whole frequency band data. (a) velocity, (b) density.

Figure 11. The single parameter inversion for velocity with the data in four different frequency bands. (a) 0–2.5 Hz, (b) 0–5 Hz, (c) 0–15 Hz, (d) all frequencies.