This invention discloses a method for
elastic modulus inversion based on
diffraction wave information. The method includes: acquiring a three-dimensional seismic data volume of
diffraction waves; obtaining near-, middle-, and far-field angularly stacked data of
diffraction waves based on the three-dimensional seismic data volume; obtaining a P-wave diffraction coefficient formula suitable for pre-stack inversion based on the P-wave
reflection coefficient formula; determining the diffraction coefficient based on the P-wave diffraction coefficient formula; obtaining near-, middle-, and far-field angularly stacked data of the whole wavefield based on layer interpretation; obtaining
wavelet data corresponding to near-, middle-, and far-field gathers based on the near-, middle-, and far-field angularly stacked data of the whole wavefield; converting the layer velocity in the whole-wavefield
seismic processing into layer velocity P-wave and S-wave modulus and density; using the layer velocity P-wave and S-wave modulus and density as background parameters of the inversion objective function; constructing and inverting the diffraction wave
elastic modulus inversion objective function based on the
wavelet data and diffraction coefficients; obtaining the
elastic modulus anomaly in the target layer caused by heterogeneous disturbance by solving the objective function; and using the elastic modulus anomaly to identify and predict heterogeneous bodies in the strata.