“…In the equations mentioned above, the subscripts i and j in brackets refer to spatial index of nodes in x and z directions, respectively, the superscript n represents the time discrete index, τ rr , τ θθ and τ zz refer to normal stress in x, θ and z directions, τ rz is shear stress, ∆t is the time step, L 1 R and L 1 Z are second-order difference operators of space, and the parameters ρ and µ represent the density and shear moduli of the medium. On the medium interface, the arithmetic mean value of ρ and the harmonic mean value of µ are employed; on the solid-liquid interface, µ is set to 0 [6,7]. σ i f (i,j) and σ j f (i,j) are arithmetic average operators and defined as:…”