Marine surveys conducted with multiple source arrays

1. In a process for generating an image of a subterranean formation using marine seismic techniques in which multiple sources are activated in a body of water above the subterranean formation to create overlapping source wavefields that are reflected from structures and layers within the subterranean formation and recorded as blended seismic data by receivers located in the body of water, the specific improvement comprising:

transforming the recorded blended seismic data from a space-time domain to a space-frequency domain to obtain blended seismic data matrices, each blended seismic data matrix associated with a different frequency in the space-frequency domain;

for each blended seismic data matrix,

calculating a pseudo-deblended seismic data matrix based on the blended seismic data matrix and a source blending matrix that models physical interactions of the overlapping source wavefields,

attenuating noise in the pseudo-deblended seismic data matrix, and

calculating an estimated deblended seismic data matrix based at least in part on the pseudo-deblended seismic data matrix;

forming deblended seismic data from the estimated deblended seismic data matrices; and

using the deblended seismic data to generate the image of the subterranean formation, the image revealing structures and layers of the subterranean formation.

2. The process of claim 1 , wherein the multiple source array comprises a vertical source array in which the sources are substantially linearly arranged in a vertical direction.

3. The process of claim 1 , wherein the multiple source array comprises a bent source array in which the sources are arranged along diagonal directions with respect to a vertical direction.

4. The process of claim 1 , wherein calculating the pseudo-deblended seismic data matrix further comprises:

forming the source blending matrix based on times when each source of the multiple source array is activated; and

calculating a pseudo-deblended seismic data matrix as a product of the blended seismic data matrix and conjugate transpose of the source blending matrix.

5. The process of claim 1 , wherein calculating the estimated deblended seismic data matrix based at least in part on the source blending matrix further comprises:

sorting the blended seismic data into common-receiver-station gathers, each common-receiver-station gather having one or more traces associated with one or more sources of the multiple source array and one or more traces associated with one or more source ghosts;

for each common-receiver-station gather,

transforming the common-receiver-station gather from a space-time domain to a wavenumber-frequency domain,

muting seismic data with wavenumbers corresponding to up-going wavefields to generate a source-deghosted common-receiver-station gather without reflection events associated with ghosts of the multiple source array, and

iteratively calculating an estimated deblended common-receiver-station gather from the source-deghosted common-receiver-station gather, and

collecting the estimated deblended common-receiver-station gathers to form the estimated deblended seismic data matrix.

6. The process of claim 5 , wherein the common-receiver-station gathers comprise two-dimensional common-receiver-station gathers.

7. The process of claim 5 , wherein sorting the blended seismic data into the common-receiver-station gathers comprises filtering noise in a frequency-wavenumber domain of two or more wavenumber dimensions.

8. A computer system for generating an image of a subterranean formation from blended seismic data recorded in a marine survey of the subterranean formation computer system comprising:

one or more processors;

one or more data-storage devices; and

machine-readable instructions stored in the one or more data-storage devices that when executed using the one or more processors controls the system to carry out

transforming the recorded blended seismic data from a space-time domain to a space-frequency domain to obtain blended seismic data matrices, each blended seismic data matrix associated with a different frequency in the space-frequency domain;

for each blended seismic data matrix,

calculating a pseudo-deblended seismic data matrix based on the blended seismic data matrix and a source blending matrix that models physical interactions of overlapping source wavefields generated by a multiple source array,

attenuating noise in the pseudo-deblended seismic data matrix, and

calculating an estimated deblended seismic data matrix based at least in part on the pseudo-deblended seismic data matrix;

forming deblended seismic data from the estimated deblended seismic data matrices; and

using the deblended seismic data to generate the image of the subterranean formation, the image revealing structures and layers of the subterranean formation.

9. The system of claim 8 , wherein calculating the pseudo-deblended seismic data matrix further comprises:

forming the source blending matrix based on times when each source of the multiple source array is activated; and

calculating the pseudo-deblended seismic data matrix as a product of the blended seismic data matrix and conjugate transpose of the source blending matrix.

10. The system of claim 8 , wherein calculating the estimated deblended seismic data matrix based at least in part on the source blending matrix comprises:

sorting the blended seismic data into common-receiver-station gathers, each common-receiver-station gather having one or more traces associated with one or more sources of the multiple source array and one or more traces associated with one or more source ghosts;

for each common-receiver-station gather,

transforming the common-receiver-station gather from a space-time domain to a wavenumber-frequency domain,

muting seismic data with wavenumbers corresponding to up-going wavefields to generate a source-deghosted common-receiver-station gather without reflection events associated with ghosts of the multiple source array, and

iteratively calculating an estimated deblended common-receiver-station gather from the source-deghosted common-receiver-station gather; and

collecting the estimated deblended common-receiver-station gathers to form the estimated deblended seismic data matrix.

11. The system of claim 10 , wherein the common-receiver-station gathers comprise two-dimensional common-receiver-station gathers.

12. The system of claim 10 , wherein sorting the blended seismic data into the common-receiver-station gathers comprises filtering noise in a frequency-wavenumber domain of two or more wavenumber dimensions.

13. A non-transitory computer-readable medium having machine-readable instructions encoded thereon for enabling one or more processors of a computer system to perform the operations of

transforming blended seismic data from a space-time domain to a space-frequency domain to obtain blended seismic data matrices, each blended seismic data matrix associated with a different frequency in the space-frequency domain;

for each blended seismic data matrix,

calculating a pseudo-deblended seismic data matrix based on the blended seismic data matrix and a source blending matrix that models physical interactions of overlapping source wavefields generated by a multiple source array,

attenuating noise in the pseudo-deblended seismic data matrix, and

calculating an estimated deblended seismic data matrix based at least in part on the pseudo-deblended seismic data matrix;

forming deblended seismic data from the estimated deblended seismic data matrices; and

using the deblended seismic data to generate the image of the subterranean formation, the image revealing structures and layers of the subterranean formation.

14. The medium of claim 13 , wherein calculating the pseudo-deblended seismic data matrix further comprises:

forming the source blending matrix based on times when each source of the multiple source array is activated; and

calculating the pseudo-deblended seismic data matrix as a product of the blended seismic data matrix and conjugate transpose of the source blending matrix.

15. The medium of claim 13 , wherein calculating the estimated deblended seismic data matrix based at least in part on the source blending matrix further comprises:

sorting seismic data into common-receiver-station gathers, each common-receiver-station gather having one or more traces associated with one or more sources of the multiple source array and one or more traces associated with one or more source ghosts;

for each common-receiver-station gather,

transforming the common-receiver-station gather from a space-time domain to a wavenumber-frequency domain,

muting seismic data with wavenumbers corresponding to up-going wavefields to generate a source-deghosted common-receiver-station gather without reflection events associated with ghosts of the multiple source array, and

iteratively calculating an estimated deblended common-receiver-station gather from the source-deghosted common-receiver-station gather, and

collecting the estimated deblended common-receiver-station gathers to form the estimated deblended seismic data matrix.

16. The medium of claim 15 , wherein the common-receiver-station gathers comprise two-dimensional common-receiver-station gathers.

17. The medium of claim 15 , wherein sorting the blended seismic data into the common-receiver-station gathers comprises filtering noise in a frequency-wavenumber domain of two or more wavenumber dimensions.

18. A method of manufacturing a geophysical data product, the method comprising:

accessing blended seismic data acquired in a marine survey of a subterranean formation using a multiple source array;

transforming the blended seismic data from a space-time domain to a space-frequency domain to obtain blended seismic data matrices, each blended seismic data matrix associated with a different frequency in the space-frequency domain;

for each blended seismic data matrix,

calculating a pseudo-deblended seismic data matrix from the blended seismic data matrix and a source blending matrix that models physical interactions of overlapping source wavefields generated by the multiple source array,

attenuating noise in the pseudo-deblended seismic data matrix, and

calculating an estimated deblended seismic data matrix based at least in part on the pseudo-deblended seismic data matrix;

forming deblended seismic data from the estimated deblended seismic data matrices;

using the deblended seismic data to generate an image of the subterranean formation, the image revealing structures and layers of the subterranean formation; and

storing the image on a non-transitory computer-readable medium.