Amplitude preserving prestack migration method
A method is disclosed for amplitude preserving Kirchhoff time migration. The method includes calculating estimates of the geometrical spreading terms using a constant velocity approximation. Takeoff and emergence angles are accurately calculated based on a model of velocity which varies with respect to travel time. An image can be calculated at at least one travel time using the estimated weights and calculated takeoff and emergence angles.
1. A method for amplitude preserving Kirchhoff time migration, comprising:
calculating, for at least one image point, estimates of geometrical spreading terms using a constant velocity approximation;
calculating, for the at least one image point, takeoff and emergence angles using a model of velocity which varies with respect to at least one of migrated time and depth;
estimating, for the at least one image point, a value of a weight function from the estimated geometrical spreading terms and calculated takeoff and emergence angles; and
storing the estimated value of the weight function in a computer readable medium.
2. The method of claim 1 further comprising generating an image at at least one sample time using the estimated weight function values.
3. The method of claim 2 wherein the generating an image comprises calculating travel times to at least one image point for a plurality of seismic source and seismic receiver positions, and calculating a Kirchhoff boundary integral based on the travel times and the weight function values.
4. A computer program stored in a computer readable medium, the program having logic operable to cause a programmable computer to perform steps comprising:
calculating, for at least one image point, estimates of geometrical spreading terms for amplitude preserving Kirchhoff time migration using a constant velocity approximation;
calculating, for at least one image point takeoff and emergence angles using a model of velocity which varies with respect to at least one of migrated time and depth; and
estimating, for the at least one image point, a value of a weight function from the estimated geometrical spreading terms and calculated takeoff and emergence angles.
5. The program of claim 4 further comprising logic operable to cause the computer to perform generating an image at at least one sample time using the estimated weight function values.
6. The program of claim of claim 5 wherein the generating an image comprises calculating travel times to at least one image point for a plurality of seismic source and seismic receiver positions, and calculating a Kirchhoff boundary integral based on the travel times and the weight function values.
7. A method for amplitude preserving Kirchhoff time migration, comprising:
calculating, for at least one image point, estimates of geometrical spreading terms using a constant velocity approximation;
calculating, for the at least one image point, takeoff and emergence angles using a model of velocity which varies with respect to at least one of migrated time and depth;
estimating, for the at least one image point, a value of a weight function from the estimated geometrical spreading terms and calculated takeoff and emergence angles;
generating an image at at least one sample time using the estimated weight function values; and
displaying the image.
8. The method of claim 7 wherein the generating an image comprises calculating travel times to at least one image point for a plurality of seismic source and seismic receiver positions, and calculating a Kirchhoff boundary integral based on the travel times and the weight function values.
9. A method for amplitude preserving Kirchhoff time migration processing of seismic data, comprising:
accessing recorded signals, the signals corresponding to measured seismic amplitude with respect to time of actuation of a seismic energy source recorded at a plurality of geodetic source and receiver positions proximate a top of an area of the Earth's subsurface to be surveyed;
calculating, for at least one subsurface point in the area, estimates of geometrical spreading terms for each of the accessed signals using a constant velocity approximation;
calculating, for the at least one image point, takeoff and emergence angles using a model of velocity which varies with respect to at least one of migrated time and depth for each of the accessed signals;
estimating, for the at least one image point, values of a weight function from the estimated geometrical spreading terms and calculated takeoff and emergence angles.
10. The method of claim 9 further comprising generating an image at at least one sample time using the estimated weight function values.
11. The method of claim 10 wherein the generating an image comprises calculating a Kirchhoff boundary integral based on the weight function values.