Media parameter-modified method for realizing an adaptive expression of an arbitrary discontinuous surface

A media parameter-modified method for realizing an adaptive expression of an arbitrary discontinuous surface, comprising following steps:

• S1:
• importing an initial forward model, and the initial forward model includes a velocity and a density of vertical and horizontal waves, setting a seismic wavelet type and a main frequency parameter, and setting a space step and a time step according to initial forward model parameters;
• S2:
• selecting the space step and the time step according to the initial forward model parameters;
• and setting a differential order-number and a seismic wavelet;
• S3:
• starting a stepped discretization of a free surface of the initial forward model when a surface of the initial forward model fluctuates, and turning an irregular undulating ground into a regular stepped grid;
• S4:
• using a media Parameter-modified method to correct a constitutive relationship and a density of non-horizontal points on the free surface;
• S5:
• using a corrected constitutive relationship to replace an original constitutive relationship on the free surface, and an influence of the free surface introduced without correcting an original finite difference method code;
• and
• S6:
• for generalized anisotropic media, including orthorhombic anisotropic media and triclinic anisotropic media, bringing the corrected constitutive relationship into a finite difference method code to introduce the influence of the free surface.