# US 10,310,111 B2 · Wave-fields separation for seismic recorders distributed at non-flat recording surfaces > Agents: this data is available as JSON from the patentgrep API and MCP server — see https://patentgrep.com/llms.txt - Number: US 10,310,111 B2 - Document ID: US10310111B2 - Status: Granted - Kind: B2 - Application number: 13/670,900 - Filed: 2012-11-07 - Published: 2013-06-20 - Granted: 2019-06-04 - Earliest claimed: 2011-12-15 - Assignees: Cgg Services SAS - Applicants: Cggveritas Services SA - Inventors: Sergio Grion - Web page: https://patentgrep.com/patent/US10310111B2 ## CPC - G01V 1/28: Processing seismic data, e.g. for interpretation or for event detection (G01V1/48 takes precedence) - G01V 1/368: {Inverse filtering} - G01V 1/38: specially adapted for water-covered areas (G01V1/28 takes precedence) - G01V 1/3852: {to the seabed} - G01V 2210/56: De-ghosting; Reverberation compensation - G01V: GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS (means for indicating the location of accidentally buried, e.g. snow-buried, persons A63B29/02) ## Abstract Apparatus, computer instructions and method for separating up-going and down-going wave fields (U, D) from seismic data recorded within or beneath a body of water, or in general below the surface of the earth. The method includes a step of receiving seismic data (Po, Zo) recorded in the time-space domain with seismic recorders distributed on a first datum, wherein the first datum is non-flat; a step of establishing a mathematical relation between transformed seismic data (P, Z) and the up-going and down-going wave fields (U, D) on a second planar datum; and a step of solving with an inversion procedure, run on a processor, the mathematical relation to obtain the up-going and down-going wave fields (U, D) for the second datum. The second datum is different from the first datum. ## Claims ### Claim 1 A method for generating an image based on separating up-going and down-going wave fields (U, D) in seismic data related to a subsurface of a body of water, or to a subsurface of a body of rock, the method comprising: receiving seismic data (Po, Zo) recorded in the time-space domain with seismic recorders distributed on a first datum, wherein the first datum is non-flat; transforming the seismic data (P0, Z0) from the time-space domain to transformed seismic data for water pressure (P) and transformed seismic data for particle velocity (Z), in another domain; establishing, in the another domain, a mathematical relation that relates (1) the transformed seismic data for water pressure (P) to desired up-going and down-going wave fields (U, D), wherein the transformed seismic data for water pressure (P) is expressed as a first function of the desired up-going and down-going wave fields (U, D) on a second planar datum, and (2) the transformed seismic data for particle velocity (Z) to the desired up-going and down-going wave fields (U, D), wherein the transformed seismic data for particle velocity (Z) is expressed as a second function of the desired up-going and down-going wave fields (U, D) on the second planar datum; solving with an inversion procedure, run on a processor, the first and second functions to obtain the desired up-going and down-going wave fields (U, D) for the second datum; and generating the image based on at least one of the desired up-going and down-going wave fields (U, D) for the second datum, wherein the second datum is different from the first datum, wherein each of the first and second functions includes extrapolator operators for extrapolating the desired up-going and down-going wave fields (U, D), respectively, to the second datum. ### Claim 2 (depends on 1) The method of claim 1, wherein the first and second functions are sums over a horizontal wave-number. ### Claim 3 (depends on 1) The method of claim 1, wherein the step of establishing further comprises: using first and second extrapolation operators to multiply the up-going and down-going wave fields (U, D), respectively, in the another domain. ### Claim 4 (depends on 1) The method of claim 1, wherein the another domain is an f-k domain, wherein f is a frequency and k is a wave-number. ### Claim 5 (depends on 1) The method of claim 1, wherein the seismic data is recorded with seismic receivers placed on a bottom of the body of water. ### Claim 6 (depends on 1) The method of claim 1, wherein the seismic data is recorded with seismic receivers placed on streamers within the body of water having variable depths. ### Claim 7 (depends on 1) The method of claim 1, where the seismic data is recorded with seismic receivers placed beneath the earth's surface. ### Claim 8 (depends on 1) The method of claim 1, wherein the seismic data is recorded with multi-component seismic receivers. ### Claim 9 (depends on 1) The method of claim 1, wherein the second datum is flat. ### Claim 10 A computing device for generating an image based on separating up-going and down-going wave fields (U, D) in seismic data related to a subsurface of a body of water or to a subsurface of a body of rock, the computing device comprising: an interface configured to receive seismic data (Po, Zo) recorded in the time-space domain with seismic recorders distributed on a first datum, wherein the first datum is non-flat; and a processor connected to the interface and configured to, transform the seismic data (P0, Z0) from the time-space domain to transformed seismic data for water pressure (P) and transformed seismic data for particle velocity (Z), in another domain, receive a mathematical relation that relates (1) the transformed seismic data for water pressure (P) to desired up-going and down-going wave fields (U, D), wherein the transformed seismic data for water pressure (P) is expressed as a first function of the desired up-going and down-going wave fields (U, D) on a second planar datum, and (2) the transformed seismic data for particle velocity (Z) to the desired up-going and down-going wave fields (U, D), wherein the transformed seismic data for particle velocity (Z) is expressed as a second function of the desired up-going and down-going wave fields (U, D) on the second planar datum, solve with an inversion procedure the first and second functions to obtain the up-going and down-going wave fields (U, D) for the second datum, and generate the image based on at least one of the desired up-going and down-going wave fields (U, D) for the second datum, wherein the second datum is different from the first datum, and wherein each of the first and second functions includes extrapolator operators for extrapolating the desired up-going and down-going wave fields (U, D), respectively, to the second datum. ### Claim 11 (depends on 10) The device of claim 10, wherein the first and second functions are sums over a horizontal wave-number. ### Claim 12 (depends on 10) The device of claim 10, wherein the processor is further configured to: use first and second extrapolation operators to multiply the up-going and down-going wave fields (U, D), respectively, in the another domain. ### Claim 13 (depends on 10) The device of claim 10, wherein the another domain is an f-k domain, wherein f is a frequency and k is a wave-number. ### Claim 14 (depends on 10) The device of claim 10, wherein the seismic data is recorded with seismic receivers placed on a bottom of the body of water. ### Claim 15 (depends on 10) The device of claim 10, wherein the seismic data is recorded with seismic receivers placed on streamers within the body of water having variable depths. ### Claim 16 (depends on 10) The method of claim 10, where the seismic data is recorded with seismic receivers placed beneath the earth's surface. ### Claim 17 (depends on 10) The device of claim 10, wherein the seismic data is recorded with multi-component seismic receivers. ### Claim 18 (depends on 10) The device of claim 10, wherein the second datum is flat. ### Claim 19 (depends on 10) The device of claim 10, further comprising: a monitor connected to the processor and configured to display the image of the subsurface that is calculated based on the up-going and down-going wave-fields. ### Claim 20 A non-transitory computer readable medium including computer executable instructions, wherein the instructions, when executed by a processor, implement instructions for generating an image based on separating up-going and down-going wave fields (U, D) in seismic data related to a subsurface of a body of water or to a subsurface of a body of rock, the instructions comprising: receiving seismic data (Po, Zo) recorded in the time-space domain with seismic recorders distributed on a first datum, wherein the first datum is non-flat; transforming the seismic data (P0, Z0) from the time-space domain to transformed seismic data for water pressure (P) and transformed seismic data for particle velocity (Z), in another domain; establishing, in the another domain, a mathematical relation that relates (1) the transformed seismic data for water pressure (P) to desired up-going and down-going wave fields (U, D), wherein the transformed seismic data for water pressure (P) is expressed as a first function of the desired up-going and down-going wave fields (U, D) on a second planar datum, and (2) the transformed seismic data for particle velocity (Z) to the desired up-going and down-going wave fields (U, D), wherein the transformed seismic data for particle velocity (Z) is expressed as a second function of the desired up-going and down-going wave fields (U, D) on the second planar datum; solving with an inversion procedure, run on the processor, the first and second functions to obtain the up-going and down-going wave fields (U, D) for the second datum; and generating the image based on at least one of the desired up-going and down-going wave fields (U, D) for the second datum, wherein the second datum is different from the first datum, and wherein each of the first and second functions includes extrapolator operators for extrapolating the desired up-going and down-going wave fields (U, D), respectively, to the second datum.