Method and device for enhanced seismic imaging based on one-way wave equation
A method for seismic processing includes steps of receiving seismic data comprising seismic waveforms acquired from a plurality of seismic receivers and a plurality of seismic sources in a subterranean area; generating a plurality of full wavefields at a plurality of corresponding propagation time steps based on a seismic velocity and a density model of the subterranean area; and applying a wavefield filter to the plurality of full wavefields to obtain a plurality of filtered wavefields. Each of the plurality of full wavefields propagates in all directions. On the other hand, each of the plurality of filtered wavefields is attenuated in or opposite to a target direction. The target direction is selected by the user.
1. A method for seismic processing, comprising:
conducting seismic survey using a plurality of seismic receivers and one or more seismic sources to acquire seismic data comprising seismic waveforms in a subterranean area;
generating, using the acquired seismic data, a plurality of full wavefields at a plurality of corresponding propagation time steps based on a seismic velocity and a density model of the subterranean area; and
applying a wavefield filter to the plurality of full wavefields to obtain a plurality of filtered wavefields,
wherein each of the plurality of filtered wavefields is attenuated in or opposite to a target direction,
wherein each of the plurality of full wavefields propagate in all directions,
wherein a time derivative of each of the plurality of filtered wavefields comprises half of a difference between a time derivative of the corresponding full wavefield and a product between a wave function in the target direction and a spatial gradient of the corresponding full wavefield,
wherein the time derivative of each of the plurality of filtered wavefields is expressed as
∂
p
n
→
∂
t
=
1
2
(
∂
p
∂
t
±
rv
∇
p
)
wherein ∥⋅∥ represents L2 norm, v represents acoustic velocity, and
wherein r is a continuous weighting factor defined as:
r
=
{
1.0
x
≥
x
1
sin
(
π
2
*
x
x
1
)
-
x
1
<
x
<
x
1
-
1.0
x
≤
-
x
1
,
wherein x={right arrow over (n)}·∇p, x 1 =α∥∇ p ∥, and α is a pre-set value between 0 to 1.
2. The method of claim 1 , wherein the plurality of seismic receivers are geophones.
3. The method of claim 1 , wherein the one or more seismic sources are chosen from a seismic truck, a wireline tool, an explosive source, or combinations thereof.