System and method for generating a seismic attenuation model

A computer system, comprising:

• one or more processors;
• memory;
• and
• one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, the one or more programs including instructions that when executed by the one or more processors cause the system to:
• a. receive, at the one or more processors, a pre-migration seismic dataset D(s, r;
• t) representative of a subsurface volume of interest wherein s indicates source location, r indicates receiver location, and t is the recorded travel time;
• b. calculate a pre-migration attenuated travel time t*(s, r;
• t);
• c. compute a time derivative of the pre-migration attenuated travel time wherein 1/Q(s, r;
• t)=∂t*(s, r;
• t)/∂t;
• d. perform a first migration on D(s,r;t) to generate common image point (CIP) gathers G(x, h) wherein x is subsurface image point and h is angle or offset;
• e. perform a second migration on D(s, r;
• t)*1/Q(s, r;
• t) to generate weighted common image point (CIP) gathers G1/Q(x, h);
• and
• f. calculate a conditioned ratio of the weighted CIP gathers G1/Q(x, h) over the CIP gathers G(x, h) to get CIP gathers of 1/Q(x, h) wherein the conditioning prevents division by zero and Q is the seismic attenuation Q-factor.

A non-transitory computer readable storage medium storing one or more programs, the one or more programs comprising instructions, which when executed by an electronic device with one or more processors and memory, cause the device to:

• a. receive, at the one or more processors, a pre-migration seismic dataset D(s, r;
• t) representative of a subsurface volume of interest wherein s indicates source location, r indicates receiver location, and t is the recorded travel time;
• b. calculate a pre-migration attenuated travel time t*(s, r;
• t);
• c. compute a time derivative of the pre-migration attenuated travel time wherein 1/Q(s, r;
• t)=∂t*(s, r;
• t)/∂t;
• d. perform a first migration on D(s,r;t) to generate common image point (CIP) gathers G(x, h) wherein x is subsurface image point and h is angle or offset;
• e. perform a second migration on D(s, r;
• t)*1/Q(s, r;
• t) to generate weighted common image point (CIP) gathers G1/Q(x, h);
• and
• f. calculate a conditioned ratio of the weighted CIP gathers G1/Q(x, h) over the CIP gathers G(x, h) to get CIP gathers of 1/Q(x, h) wherein the conditioning prevents division by zero and Q is the seismic attenuation Q-factor.