System and method for multiple prediction with angular dependent reflectivity

A computing system for removing multiples m from seismic data d associated with a subsurface, the computing system comprising:

• an interface configured to receive the seismic data d associated with the subsurface, wherein the seismic data d includes primaries and multiples m recorded by a receiver;
• and
• a processor connected to the interface and configured to,
• forward propagate the seismic data d into the subsurface to form forward propagated data dΦ;
• receive angular dependent reflectivities r74 , associated with an angular range Δθ, in the subsurface, for a given point;
• generate an angle dependent reflecting wavefield DΦθrθ based on the forward propagated data dΦ and the angular dependent reflectivities rθ;
• calculate a multiple model ΦDΦθrθ by forward propagating the angle dependent reflecting wavefield DΦθrθ to the receiver;
• attenuate multiplies m associated with the multiple model, from the seismic data d, by subtraction of the multiple model to calculate demultiple data dd;
• and
• generate an image of the subsurface indicative of geophysical features associated with a natural resource, carbon capture monitoring or wind turbine placement, based on the demultiple data dd.