Near-surface geomorphological characterization based on remote sensing data

• US 8,615,362 B2
• Filed: 09/28/2009
• Issued: 12/24/2013
• Est. Priority Date: 10/10/2008
• Status: Expired due to Fees

• Current Assignee: WesternGeco LLC
• Original Assignee: WesternGeco LLC
• Inventors: Laake, Andreas
• Primary Examiner(s): Hellner, Mark
• Assistant Examiner(s): Diacou, Ari M
• Application Number:
  US12/568,322
• Publication Number:
  US 20100091611A1
• Time in Patent Office: 1,548 Days
• Field of Search: 702/5
• US Class Current: 702/5
• CPC Class Codes: G01S 13/865 Combination of radar systems with lidar systems
• CPC Class Codes: G01V 1/30 Analysis G01V1/50 takes precedence
• CPC Class Codes: G01V 8/00 Prospecting or detecting by optical means
• CPC Class Codes: G06T 17/05 Geographic models

Abstract

A technique includes receiving first data belonging to a first type of remote sensing data for a region of interest and receiving second data belonging to a different second type of remote sensing data for the region of interest. The technique includes determining at least one geomorphological feature of the region of interest based at least in part on the first and second data.

Claims

1. A method comprising:

receiving first data belonging to a first type of remote sensing data for a region of interest;

receiving second data belonging to a different second type of remote sensing data for the region of interest; and

processing the first and second data in a processor-based machine to determine a risk map for the region of interest based at least in part on the first and second data, wherein the risk map represents a risk of scattering for a geologic survey conducted in the region of interest.

2. The method of claim 1, wherein:

the first data comprises data representing one of the members of a set consisting essentially of an elevation model of the region of interest, satellite-acquired spectroscopic imagery of the region of interest and satellite-acquired radar imagery of the region of interest, and

the second data comprises data representing another one of the members of the set.

3. The method of claim 1, wherein:

the first data comprises data representing an elevation model of the region of interest, and

the second data comprises data representing a satellite-acquired spectroscopic image of the region of interest.

4. The method of claim 1, wherein the risk map further represents logistical risks for the survey.

5. The method of claim 1, wherein processing the first and second data to determine the risk map comprises:

determining topography features of the region of interest based on the first data; and

determining lithology features of the region of interest based on the second data.

6. The method of claim 5, wherein the act of determining the topography features comprises:

determining a histogram of the elevations of the region of interest;

determining a hypsographic curve based on the histogram; and

determining a classification map of the topography of the region of interest based on the hypsographic curve.

7. The method of claim 5, wherein the act of determining the topography features comprises:

determining locations of escarpments based on a multi-azimuth edge detection algorithm.

8. The method of claim 5, wherein the act of determining the lithology features comprises:

analyzing color bands of the satellite-acquired spectroscopic imagery.

9. The method of claim 8, wherein the act of analyzing the color bands comprises:

for each mineral of a predetermined group of minerals, determining a band ratio mask and applying the band ratio mask to generate a map of the mineral in the region of interest.

10. The method of claim 5, wherein

the act of determining the topography features comprises determining an escarpment map,

the act of determining the lithology features comprises determining a hard rock mask and determining a clay mask, and

processing the first and second data to determine the risk map comprises generating a lithology risk map based on the escarpment map, the hard rock mask and the clay mask.

11. The method of claim 5, wherein

the act of determining the topography features comprises determining an escarpment map,

the act of determining the lithology features comprises determining a lithology boundary map, and

processing the first and second data to determine the risk map comprises generating a data quality scatter risk map based on the escarpment map and the lithology boundary.

12. The method of claim 5, wherein

the act of determining the topography features comprises determining an escarpment map,

the act of determining the lithology features comprises determining a hard rock map, a clay mask map and a soft rock mask map, and

processing the first and second data to determine the risk map comprises generating a data quality scatter risk map based on the escarpment map, the hard rock map, the clay mask map and the soft rock mask map.

13. The method of claim 1, wherein

the scattering comprises lithologic scattering or topographic scattering.

14. The method of claim 1, wherein

the risk map further represents a risk of vibrator point loading in connection with the survey.

15. The method of claim 1, wherein

the risk map further indicates a risk for vehicle access in the region of interest.

16. The method of claim 1, further comprising

determining static corrections based at least in part on the first and second data.

17. The method of claim 16, wherein the act of determining the static corrections comprises determining a geologic near surface model for the region of interest based at least in part on the first and second data.

18. The method of claim 17, wherein the act of determining the static corrections further comprises determining an elastic near surface model for the region of interest based at least in part on the geologic near surface model.

19. The method of claim 18, wherein the act of determining the static corrections further comprises calibrating the elastic near surface model based on refraction static estimates based on near surface seismic data acquired in a seismic survey.

20. The method of claim 1, wherein

the first data comprises data representing an elevation model of the region of interest,

the second data comprises data representing a satellite-acquired spectroscopic image of the region of interest, and

processing the first and second data comprises determining a pre-erosion map of the region of interest.

21. The method of claim 20 wherein the act of determining the pre-erosion map comprises:

applying a spatial gradient filter to extract terrain edges and escarpments from the elevation model to produce a filtered image; and

generating the pre-erosion map based on the filtered image and the multispectral imagery.

22. The method of claim 21, further comprising:

identifying a graben fault in the pre-erosion map; and

based on the pre-erosion map determining a throw of the graben fault.

23. A system comprising:

an interface to receive first data belonging to a first type of remote sensing data for a region of interest and receive second data belonging to a different second type of remote sensing data for the region of interest; and

a processor to determine a risk map for the region of interest based at least in part on the first and second data, the risk map representing a risk of scattering for a geologic survey conducted in the region of interest.

24. The system of claim 23, wherein

the first data comprises data representing one of the members of a set consisting essentially of an elevation model of the region of interest, satellite-acquired spectroscopic imagery of the region of interest and satellite-acquired radar imagery of the region of interest, and

the second data comprises data indicative of another one of the members of the set.

25. The system of claim 23, wherein

the first data comprises data representing an elevation model of the region of interest, and

the second data comprises data indicative of a satellite-acquired spectroscopic image of the region of interest.

26. The system of claim 23, wherein the risk map further represents logistical risks for the survey.

27. The system of claim 23, wherein the processor is adapted to:

determine topography features of the region of interest based on the first data;

determine lithology features of the region of interest based on the second data; and

combine the determined topography and lithology features to generate the risk map.

28. The system of claim 23, wherein the processor is adapted to:

determine a histogram of the elevations of the region of interest;

determine a hypsographic curve based on the histogram; and

determine a classification map of the topography of the region of interest based on the hypsographic curve.

29. The system of claim 23, wherein the processor is adapted to:

determine locations of escarpments based on a multi-azimuth edge detection algorithm.

30. The system of claim 23, wherein at least one of the first and second data represents

color bands of the satellite-acquired spectroscopic imagery, and the processor is adapted to analyze color bands of the satellite-acquired spectroscopic imagery.

31. The system of claim 30, wherein the processor is adapted to use the analysis of the color bands to, for each mineral of a predetermined group of minerals, determine a band ratio mask and apply the band ratio mask to generate a map of the mineral in the region of interest.

32. The system of claim 23, wherein the processor is adapted to:

determine an escarpment map;

determine a hard rock mask;

determine a clay mask; and

generate a lithology risk map based on the escarpment map, the hard rock mask and the clay mask.

33. The system of claim 23, wherein the processor is adapted to:

determine an escarpment map;

determine a lithology boundary map; and

generate a data quality scatter risk map based on the escarpment map and the lithology boundary map.

34. The system of claim 23, wherein the processor is adapted to:

determine an escarpment map,

determine a hard rock map;

determine a clay mask map;

determine a soft rock mask map, and

generate a data quality velocity risk map based on the escarpment map, the hard rock map, the clay mask map and the soft rock mask map.

35. The system of claim 23, wherein:

the risk map represents of a risk of lithologic scattering or-topographic scattering.

36. The system of claim 23, wherein:

the risk map further represents a risk of vibrator point loading in connection with the survey.

37. The system of claim 23, wherein:

the risk map further represents a risk of noise from trapped surface modes in connection with a seismic survey conducted in the region of interest.

38. The system of claim 23, wherein:

the risk map further represents a risk for vehicle access in the region of interest.

39. The system of claim 23, wherein:

the first data comprises data representing an elevation model of the region of interest,

the second data comprises data representing satellite-acquired multispectral imagery of the region of interest, and

the processor is adapted to determine static corrections based at least in part on the first and second data.

40. The system of claim 23, wherein the processor is adapted to determine a geologic near surface model for the region of interest based at least in part on the first and second data.

41. The system of claim 40, wherein the processor is adapted to determine an elastic near surface model for the region of interest based at least in part on the geologic near surface model.

42. The system of claim 41, wherein the processor is adapted to calibrate the elastic near surface model based on refraction static estimates based on near surface seismic data acquired in a seismic survey.

43. The system of claim 23, wherein:

the first data comprises data representing an elevation model of the region of interest,

the second data comprises data representing a satellite-acquired spectroscopic image of the region of interest, and

the processor is adapted to determine a pre-erosion map of the region of interest.

44. The system of claim 43, wherein the processor is adapted to:

apply a spatial gradient filter to extract terrain edges and escarpments from the elevation model to produce a filtered image; and

combine the filtered image with the satellite-acquired multispectral imagery to generate the pre-erosion map.

45. The system of claim 44, wherein the processor is adapted to identify a graben fault in the pre-erosion map and based on the pre-erosion map, determine a throw of the graben fault.

46. A method comprising:

receiving first data representing different bands of satellite acquired multispectral imagery;

processing the first data in a processor-based machine to generate second data representing a band ratio image;

identifying wadis based at the band ratio image; and

applying a pattern tracking algorithm to identify at least one path formed by the wadis to identify an active fault.

47. The method of claim 46, wherein the act of processing the first data comprises:

processing data indicative of pan-chromatic and thermal infrared bands such that the band ratio image is formed from a ratio of the pan-chromatic and thermal infrared bands.