The invention relates to a
time domain acoustic wave equation explicit finite difference seismic response
simulation method. The method includes the following steps that: 1) the initial value and allowable maximum
approximation error xi of
spatial difference operator length 2M+1 are set; 2) on the basis of space-wave number domain
function approximation-based spatial
finite difference coefficient solving method, the difference coefficient and wave number domain
error function value of an explicit finite difference formula which is built based on K-space operator compensation numerical value
simulation error are solved according to the
acoustic wave speed of a target region; 3) if the wave number domain
error function value is larger than the allowable maximum
approximation error xi, 1 is added to the value of M, and the method returns to step 2), if the wave number domain
error function value is smaller than the allowable maximum
approximation error xi, and M is larger than 1, 1 is subtracted from the value of M, and the method returns to step 2), otherwise, a minimum
spatial difference operator length 2M+1 which satisfies the allowable maximum approximation error xi and a difference coefficient corresponding to the minimum
spatial difference operator length 2M+1 are determined, a next step is continued; and 4) seismic response numerical value
simulation is performed on the target region, so that the simulation seismic records or seismic
waves of the target region can be obtained.