Source-receiver position estimation using direct arrival modeling and inversion

Patent metadata

Field	Value
Publication number	US-11442186-B2
Application number	US-201816633787-A
Country	US
Kind code	B2
Filing date	Aug 29, 2018
Priority date	Aug 31, 2017
Publication date	Sep 13, 2022
Grant date	Sep 13, 2022

Abstract

A method for correcting physical positions of seismic sensors and/or seismic sources for a seismic data acquisition system. The method includes estimating a respective energy generated by each source element, which belongs to a source array; calculating a respective energy recorded by each individual seismic sensor, which belongs to a composite receiver; summing, for each individual seismic sensor, all the generated energies from the all the source elements; estimating a model of direct arrival waves that propagate from the source elements to the individual seismic sensors; calculating positions of the individual seismic sensors based on the model of direct arrival waves; comparing calculated positions of the individual seismic sensors with observed positions of the individual seismic sensors; selecting a best calculated position for each of the individual seismic sensors based on an objective function; and correcting the observed positions of the individual seismic sensors with corresponding best calculated positions.

First claim

What is claimed is:

1. A method for correcting observed positions of seismic sensors and/or seismic sources for a seismic data acquisition system, the method comprising:
estimating a source element energy generated by each source element, which belongs to a source array;
estimating a model of direct arrival waves that propagate straight from the source elements to the individual seismic sensors, respectively, the model being based on source modelling for each of the source elements, propagation modelling and characteristics of the seismic data acquisition system related to each of individual seismic sensors of a composite receiver;
for at least two sets of positions among which at least one set of positions is obtained by perturbing at least a position of one of the seismic sensors and/or of the source elements from another of the at least two sets of positions,
calculating direct arrival energies of the direct arrival waves expected to be recorded by each individual seismic sensor and due to each of the source elements, respectively, a direct arrival energy being calculated using
the model,
the estimated source element energy of a respective source element, and
positions of the respective source element within the source array and of a respective individual seismic sensor within the composite receiver according to one of the at least two sets of positions, and

summing, for each individual seismic sensor, the direct arrival energies expected to be recorded from the source elements,
calculating values of an objective function for the at least two sets of positions, each of the values indicating a mismatch between
the direct arrival waves computed using the model for a set of positions among the at least two sets of positions, and
direct arrival waves according to observed seismic data;

selecting a best set of positions among the at least two sets of positions, based on the values of the objective function, the best set of positions corresponding to a lowest value of the objective function; and
correcting one or more of the observed positions of the individual seismic sensors and/or of the source elements according to the best set of positions.

2. The method of claim 1 , wherein the characteristics of the acquisition system include clipping a voltage recorded by an individual seismic sensor.

3. The method of claim 1 , wherein the characteristics of the acquisition system include filtering a signal recorded by an individual seismic sensor.

4. The method of claim 1 , wherein the characteristics of the acquisition system include a sampling of the individual seismic sensor.

5. The method of claim 1 , wherein the characteristics of the acquisition system include a calibration of the individual seismic sensor.

6. The method of claim 1 ,
wherein the model estimates a direct arrival wave due to a source element and reaching a seismic sensor based on
a signature of the source element,
a ghost associated with propagation from the source element to the seismic sensor, and
the estimated source element energy,

to determine an amplitude and a phase of the direct arrival wave.

7. The method of claim 1 , wherein for a first one of the at least two sets of positions, the direct arrival energy reaching each individual seismic sensor is calculated by taking into account the position of the individual seismic sensor based on global positioning system, GPS, information recorded by the acquisition system for a head and a tail of a streamer on which the individual seismic sensor is located.

8. The method of claim 1 , further comprising:
towing with a streamer vessel the individual seismic sensors along sail lines, and
towing the source array so that the source array is located directly above the individual seismic sensors, along a gravity direction.

9. A computing device for correcting observed positions of seismic sensors and/or seismic sources for a seismic data acquisition system, the computing device comprising:
an interface configured to receive a position of each source element within the source array and a position of each individual seismic sensor within a composite receiver; and
a processor connected to the interface and configured
to estimate a source element energy generated by each source element, which belongs to the source array,
to estimate a model of direct arrival waves that propagate straight from the source elements to the seismic sensors, respectively, the model being based on source modelling, propagation modelling and characteristics of the seismic data acquisition system related to each of individual seismic sensors of a composite receiver,
for at least two sets of positions among which at least one set of positions is obtained by perturbing at least a position of one of the seismic sensors and/or of the source elements from another of the at least two sets of positions,
to calculate direct arrival energies expected to be recorded by each individual seismic sensor and due to each of the source elements of the source array, respectively, a direct arrival energy being calculated using
the model,
the estimated source element energy of a respective source element, and
positions of the respective source element within the source array and of a respective individual seismic sensor within the composite receiver according to one of the at least two sets of positions, and

to sum, for each individual seismic sensor, the direct arrival energies from the source elements,
to calculate values of an objective function for the at least two sets of positions, respectively, each of the values indicating a mismatch between
the direct arrival waves computed using the model of direct arrival waves for a set of positions among the at least two sets of positions, and
direct arrival waves according to observed seismic data,

to select a best set of positions among the at least two sets of positions based on the values of the objective function, the best set of positions corresponding to a lowest value of the objective function, and
to correct one or more of the observed positions of the individual seismic sensors and/or of the source elements according to the best set of positions.

10. A method for correcting observed positions of seismic sensors and/or seismic sources for a seismic data acquisition system, the method comprising:
estimating a source element energy generated by each source element, which belongs to a source array;
estimating a model of direct arrival waves that propagate straight from the source elements to the individual seismic sensors, respectively, the model being based on source modeling, propagation modeling and characteristics of the seismic data acquisition system related to each of individual seismic sensors of a composite receiver;
calculating direct arrival energies of the direct arrival waves expected to be recorded by each individual seismic sensor and due to each source element of the source array, using the model, the estimated source element energy of a source element, an observed position of the source element within the source array and an observed position of a respective individual seismic sensor within the composite receiver;
summing, for each individual seismic sensor, the direct arrival energies from the source elements;
extracting a parameter representing mismatches between the direct arrival waves according to the model and direct arrival waves observed in seismic data recorded by the individual seismic sensors;
applying an inversion process to the parameter for determining one or more positioning residuals for one or more of the individual seismic sensors; and
correcting one or more of the observed positions of the individual seismic sensors wit

Assignees

• Cgg Services Sas
• Sercel Sas

Inventors

• Deschizeaux Bernard
• Vinje Vetle
• Nilsen Carl-Inge

Classifications

• G01V2210/67

  Wave propagation modeling · CPC title

• G01V1/3835 Primary

  measuring position, e.g. by GPS or acoustically · CPC title

• G01V1/282 Primary

  Application of seismic models, synthetic seismograms · CPC title