Method and systems for computational efficiency 3D prestack Kirchhoff depth migration
Methods and systems for forming a three-dimensional (“3D”) seismic image of a subterranean region of interest is disclosed. The method includes obtaining a seismic dataset a seismic trace for each of a plurality of pairs of one source and one receiver location and obtaining a 3D travel-time cube for each source location and each receiver location. The method further includes dividing the seismic dataset into a plurality of seismic subsets composed of set of source locations, set of receiver locations a seismic trace for each pair of source and receiver location and the 3D travel-time cube for each source for each receiver location. The method still further includes transmitting, to a random-access memory block of a computer processing unit the seismic subset, and forming a seismic partial image based on the seismic subset, and determining the 3D seismic image based on a combination of the seismic partial images.
1. A method of forming a three-dimensional (“3D”) seismic image of a subterranean region of interest, comprising:
obtaining a seismic dataset for the subterranean region of interest, wherein the seismic dataset comprises a plurality of source locations, a plurality of receiver locations, and a seismic trace for each pair of one source location and one receiver location;
obtaining a 3D travel-time cube for each source location and each receiver location;
dividing the seismic dataset into a plurality of seismic subsets, each subset comprising:
a first set of source locations;
a second set of receiver locations;
the seismic trace for each pair of one source location from the first set and one receiver location from the second set; and
the 3D travel-time cube for each source in the first set and 3D travel-time cube for each source in the second set,
wherein a first size of the first set and a second size of the second set differ by less than fifty percent of a sum of the first size and the second size;
for each seismic subset;
transmitting, to a random-access memory block of a computer processing unit (“CPU”), the seismic subset; and
forming, using the CPU, a seismic partial image based, at least in part, on the seismic subset, and
determining the 3D seismic image based, at least in part, on a combination of the seismic partial image for each seismic subset.
2. The method of claim 1 , further comprising:
identifying a portion of the subterranean region of interest with a likelihood of containing hydrocarbons based, at least in part, on the 3D seismic image of the subterranean region of interest;
determining a well path through the subterranean region of interest based, at least in part, on the identified portion of the subterranean region of interest; and
performing the well path using a drilling system.
3. The method of claim 1 , wherein the first size is equal to the second size.
4. The method of claim 1 , wherein the CPU comprises:
the CPU; and
at least one graphical processing unit.
5. The method of claim 1 , wherein the 3D seismic image comprises a seismic partial image generated by a Kirchhoff migration.
6. The method of claim 5 , wherein the seismic partial image comprises a common image gather.
7. The method of claim 1 , wherein the seismic dataset is acquired by utilizing a cross-spread acquisition geometry.
8. A non-transitory computer readable medium storing instructions executable by a computer processor, the instructions comprising functionality for:
receiving a seismic dataset for a subterranean region of interest, wherein the seismic dataset comprises a plurality of source locations, a plurality of receiver locations, and a seismic trace for each pair of one source location and one receiver location;
receiving a 3D travel-time cube for each source location and each receiver location;
dividing the seismic dataset into a plurality of seismic subsets, each subset comprising:
a first set of source locations;
a second set of receiver locations;
the seismic trace for each pair of one source location from the first set and one receiver location from the second set; and
the 3D travel-time cube for each source in the first set and 3D travel-time cube for each source in the second set,
wherein a first size of the first set and a second size of the second set differ by less than fifty percent of an average of the first size and the second size;
for each seismic subset;
transmitting, to a random-access memory block of a computer processing unit (“CPU”), the seismic subset; and
forming, using the CPU, a seismic partial image based, at least in part, on the seismic subset,
determining a 3D seismic image based, at least in part, on a combination of the seismic partial image for each seismic subset.
9. The non-transitory computer readable medium of claim 8 , the instructions further comprising:
identifying a portion of the subterranean region of interest with a high likelihood of containing hydrocarbons based, at least in part, on the 3D seismic image of the subterranean region of interest; and
determining a well path through the subterranean region of interest based, at least in part, on the identified portion of the subterranean region of interest.
10. The non-transitory computer readable medium of claim 8 , wherein the first size is equal to the second size.
11. The non-transitory computer readable medium of claim 8 , wherein the CPU comprises:
the CPU; and
at least one graphical processing unit.
12. The non-transitory computer readable medium of claim 8 , wherein the 3D seismic image comprises a seismic partial image generated by a Kirchhoff migration.
13. The non-transitory computer readable medium of claim 12 , wherein the seismic partial image comprises a common image gather.
14. The non-transitory computer readable medium of claim 8 , wherein the seismic dataset is acquired by utilizing a cross-spread acquisition geometry.
15. A system for forming three-dimensional (“3D”) seismic image of a subterranean region of interest, comprising:
a seismic acquisition system; and
a seismic processor configured to:
receive a seismic dataset from the seismic acquisition system, wherein the seismic dataset comprises a plurality of source locations, a plurality of receiver locations, and a seismic trace for each pair of one source location and one receiver location;
receive a 3D travel-time cube for each source location and each receiver location;
divide the seismic dataset into a plurality of seismic subsets, each subset comprising:
a first set of source locations;
a second set of receiver locations;
the seismic trace for each pair of one source location from the first set and one receiver location from the second set; and
the 3D travel-time cube for each source in the first set and 3D travel-time cube for each source in the second set,
wherein a first size of the first set and a second size of the second set differ by less than fifty percent of an average of the first size and the second size;
for each seismic subset;
transmit, to a random-access memory block of a computer processing unit (“CPU”), the seismic subset;
form, using the CPU, a seismic partial image based, at least in part, on the seismic subset; and
determine the 3D seismic image based, at least in part, on a combination of the seismic partial image for each seismic subset.
16. The system of claim 15 , wherein forming an objective function comprises:
identify a portion of the subterranean region of interest with a high likelihood of containing hydrocarbons based, at least in part, on the 3D seismic image of the subterranean region of interest; and
determine a well path through the subterranean region of interest based, at least in part, on the identified portion of the subterranean region of interest.
17. The system of claim 15 , wherein forming the plurality of trace pairs comprises:
the first size is equal to the second size.
18. The system of claim 15 , wherein the CPU comprises:
the CPU comprises the CPU and at least one graphical processing unit.
19. The system of claim 15 , wherein the 3D seismic image comprises a seismic partial image generated by a Kirchhoff migration.
20. The system of claim 15 , wherein the seismic partial image comprises a common image gather.