Calculation method, storage medium and device for seabed reflection coefficient of point source elastic wave

A calculation method for a seabed reflection coefficient of point source elastic wave, comprising the following steps:

• step 1:
• initializing a calculation accuracy and a calculation range of the seabed reflection coefficient;
• step 2:
• discretizing a parameter space and obtaining the seabed reflection coefficient of point source elastic wave;
• step 3:
• combining an equivalent equation with a traditional calculation equation for the seabed reflection coefficient of point source elastic wave;
• wherein the equivalent equation and traditional calculation equation for the seabed reflection coefficient of point source elastic wave are combined by the following equation:
• f(W,m)=r←r′ (1),
• in equation (1), f (w,m)=r represents the equivalent equation, W is an undetermined coefficient of the equivalent equation, r is the seabed reflection coefficient, ← represents an assignment to a variable, and r′ is the seabed reflection coefficient obtained by the traditional calculation method for the seabed reflection coefficient of point source elastic wave;
• the combination is performed within a parameter space range to obtain the following:
• f(W,M)=R (2),
• in equation (2), R is the seabed reflection coefficient of point source elastic wave;
• step 4:
• solving the undetermined coefficient of the equivalent equation in step 3 with the following equation:
• W=G−1R (3),
• in equation (3), G is a symmetric matrix Gij=φ(mi,mj) composed of Gaussian kernel function, i,j∈{k|k>0&k∈N}, R is the seabed reflection coefficient of the point source elastic wave;
• step 5:
• obtaining a concise expression for the seabed reflection coefficient of point source under an accuracy in step 1 with the following equation:
• r=g(m)·W (4),
• in equation (4), g(m)=[φ(m1,m)φ(m2,m) . . . φ(mi,m)], W is the undetermined coefficient of the equivalent equation;
• step 6:
• calculating the seabed reflection coefficient of point source elastic wave within a given calculation accuracy range using an obtained expression in step 5.