Seismic modeling study of the Earth's deep crust

We use seismic modeling methods to validate the interpretation of deep-crust seismic exploration. An approximation of the stacked section is obtained with the nonreflecting wave equation and the exploding-reflector approach. Using this technique and ray-tracing algorithms, we obtain a geological model by comparing the synthetic section with the real stacked section. An isotropic constitutive equation is assumed in this phase. The exact synthetic stacked section is then obtained by applying the standard processing sequence to a set of synthetic common-shot profiles computed with the variable-density acoustic wave equation. We introduce elliptical P-wave anisotropy and the effects of small-scale inhomogeneities by using a von Kármán autocovariance probability function that simulates scattering Q effects. Verification of the geological model by poststack migration constitutes an additional test. The methodology, which is suitable for areas of complex geology, is applied to a seismic line acquired in the northern Apennines as part of the Italian deep-crust exploration project, CROP. This area is particularly difficult to interpret because of the presence of a complex tectonic setting.