主权项 |
1. A computer-implemented method of attenuating interbed multiples in multiply-reflected seismic waves, the method comprising:
at a computer system having a processor and memory storing one or more program modules to be executed by the processor:
providing (i) seismic data comprising a plurality of beams and (ii) an earth model related to a geological volume, wherein the earth model includes an upper interface and a lower interface in the geological volume;selecting a beam in the plurality of beams as an input beam, wherein the input beam is associated with a source and a detector located near a top surface of the geological volume, and wherein the source and the detector are paired for collecting the seismic data;determining at least one of (i) a source-side stationary pegleg arrival between the upper interface and the lower interface and a corresponding detector-side primary beam and (ii) a detector-side stationary pegleg arrival between the upper interface and the lower interface and a corresponding source-side primary beam, wherein determining the source-side stationary pegleg arrival between the upper and lower interfaces and the corresponding detector-side primary beam further includes:
performing a raytracing between the source, a pegleg reflection point on the lower interface, an interbed reflection point on the upper interface, and a predefined location on the top surface using a source ray parameter;identifying a difference between (i) the raytracing property at the predefined location and (ii) a predefined raytracing property derived from the seismic data associated with the predefined location;if the difference satisfies a predefined threshold:determining the source-side stationary pegleg arrival between the upper and lower interfaces in accordance with the raytracing result; andidentifying, among the beams, a beam between the predefined location on the top surface and the detector as the corresponding detector-side primary beam;if the difference does not satisfy the predefined threshold:repeating the raytracing and identifying steps using an updated source ray parameter;predicting an interbed multiples beam between the upper interface and the lower interface using at least one of (i) the detector-side primary beam delayed by the source-side stationary pegleg arrival and (ii) the source-side primary beam delayed by the detector-side stationary pegleg arrival; anddeconvolving the predicted interbed multiples beam with the input beam thereby attenuating interbed multiples present in the input beam. |