# US 10,310,115 B2 · Prediction and subtraction of multiple diffractions > Agents: this data is available as JSON from the patentgrep API and MCP server — see https://patentgrep.com/llms.txt - Number: US 10,310,115 B2 - Document ID: US10310115B2 - Status: Granted - Kind: B2 - Application number: 15/191,887 - Filed: 2016-06-24 - Published: 2017-01-26 - Granted: 2019-06-04 - Earliest claimed: 2015-07-22 - Assignees: PGS Geophysical AS - Applicants: PGS Geophysical AS - Inventors: Robertus Franciscus Hegge; Paolo Terenghi; Andreas Klaedtke; Amr Ahmed Mahmoud Ibrahim - Web page: https://patentgrep.com/patent/US10310115B2 ## CPC - G01V 1/30: Analysis (G01V1/50 takes precedence) - G01V 1/325: {Transforming one representation into another} - G01V 1/364: {Seismic filtering (G01V1/37 takes precedence)} - G01V 1/38: specially adapted for water-covered areas (G01V1/28 takes precedence) - G01V 2210/121: Active source - G01V 2210/1423: Sea - G01V 2210/46: Radon transform - G01V 2210/48: Other transforms - G01V 2210/56: De-ghosting; Reverberation compensation - G01V 2210/675: Wave equation; Green's functions - 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 Prediction and subtraction of multiple diffractions may include transforming previously acquired seismic data from a time-space domain to a transformed domain using a dictionary of compressive basis functions and separating, within the transformed previously acquired seismic data, a first portion and a second portion of the transformed previously acquired seismic data. Prediction and subtraction of multiple diffractions may also include predicting a plurality of multiple diffractions based on the separated first and second portions and adaptively subtracting the predicted multiple diffractions from the previously acquired seismic data. Prediction and subtraction of multiple diffractions may also include inverse transforming a particular seismic data set from the transformed domain to the time-space domain. ## Claims ### Claim 1 A system, comprising: an analysis hardware engine configured to transform previously acquired seismic data acquired at a plurality of locations from a time-space domain to a transformed domain using a dictionary of asymptote- and apex-shifted hyperbolic radon transform (AASHRT) basis functions and a dictionary of Green's functions; a separation hardware engine configured to separate, within the transformed previously acquired seismic data, a first portion and a second portion of the transformed previously acquired seismic data based on an asymptote shift parameter associated with the dictionary of AASHRT basis functions and a firing delay parameter associated with the dictionary of Green's functions; a prediction hardware engine configured to predict, a plurality of multiple diffractions based on the separated first and second portions; a subtraction hardware engine configured to adaptively subtract the predicted multiple diffractions from the previously acquired seismic data; and a synthesis hardware engine configured to inverse transform a particular seismic data set from the transformed domain to the time-space domain such that the separated diffractive portion is synthesized at denser locations than the previously acquired seismic data, wherein the denser locations are more compact with respect to one another as compared to the plurality locations during acquisition of the previously acquired seismic data. ### Claim 2 (depends on 1) The system of claim 1, wherein the first portion is a predominantly reflective portion and the second portion is a predominantly diffractive portion. ### Claim 3 (depends on 2) The system of claim 2, wherein the prediction engine is configured to predict, within the transformed domain, the plurality of multiple diffractions based on the separated predominantly diffractive portion, the synthesis engine is configured to inverse transform the plurality of multiple diffractions, and the subtraction engine is configured to subtract, in the time-space domain, the inverse transformed multiple diffractions from the previously acquired seismic data. ### Claim 4 (depends on 2) The system of claim 2, wherein the synthesis engine is configured to inverse transform the separated predominantly diffractive portion, the prediction engine is configured to predict, within the time-space domain, the plurality of multiple diffractions based on the separated inverse-transformed diffractive portion, and the adaptive subtraction engine is configured to subtract, in the time-space domain, the inverse transformed predicted multiple diffractions from previously acquired seismic data. ### Claim 5 (depends on 2) The system of claim 2, wherein the prediction engine is configured to predict, within the transformed domain, the plurality of multiple diffractions based on the separated predominantly diffractive portion, the adaptive subtraction engine is configured to subtract, in the transformed domain, the predicted multiple diffractions from the transformed previously acquired seismic data, and the synthesis engine is configured to inverse transform the result to the time-space domain. ### Claim 6 (depends on 1) The system of claim 1, wherein the first portion is associated with multiple diffractions and the second portion is a remainder of the transformed previously acquired seismic data, and the synthesis engine is configured to inverse transform the remainder of the transformed previously acquired seismic data to the time-space domain. ### Claim 7 A method, comprising: transforming previously acquired seismic data acquired at a plurality of locations from a time-space domain to a transformed domain using a dictionary of asymptote- and apex-shifted hyperbolic radon transform (AASHRT) basis functions and a dictionary of Green's functions; separating, within the transformed previously acquired seismic data, a reflective portion and a diffractive portion based on an asymptote shift parameter associated with the dictionary of AASHRT basis functions and a firing delay parameter associated with the dictionary of Green's functions; predicting a plurality of multiple diffractions based on the separated diffractive portion; inverse transforming a first particular seismic data set from the transformed domain to the time-space domain such that the separated diffractive portion is synthesized at denser locations than the previously acquired seismic data, wherein the denser locations are more compact with respect to one another as compared to the plurality locations during acquisition of the previously acquired seismic data; adaptively subtracting a second particular seismic data set from the previously acquired seismic data; and recording a result of the adaptive subtraction on a non-transitory machine-readable medium. ### Claim 8 (depends on 7) The method of claim 7, wherein predicting the plurality of multiple diffractions comprises predicting, within the transformed domain, the plurality of multiple diffractions. ### Claim 9 (depends on 7) The method of claim 7, wherein predicting the plurality of multiple diffractions comprises predicting, within the time-space domain, the plurality of multiple diffractions. ### Claim 10 (depends on 9) The method of claim 9, wherein predicting the plurality of multiple diffractions within the time-space domain comprises predicting the plurality of multiple diffractions using a surface-related multiple elimination prediction based on a combination of the separated diffractive portion and the previously acquired seismic data. ### Claim 11 (depends on 7) The method of claim 7, wherein inverse transforming the first particular seismic data set comprises inverse transforming the predicted plurality of multiple diffractions. ### Claim 12 (depends on 7) The method of claim 7, wherein inverse transforming the first particular seismic data set comprises inverse transforming the separated diffractive portion. ### Claim 13 (depends on 12) The method of claim 12, wherein inverse transforming the separated diffractive portion such that the separated diffractive portion is synthesized at denser locations than the previously acquired seismic data comprises synthesizing the separated diffractive portion at a more compact location as compared to a synthesis location of the previously acquired seismic data. ### Claim 14 (depends on 7) The method of claim 7, wherein adaptively subtracting the second particular seismic data set comprises adaptively subtracting, within the transformed domain, the plurality of predicted multiple diffractions from the previously acquired seismic data. ### Claim 15 (depends on 7) The method of claim 7, wherein adaptively subtracting the second particular seismic data set comprises adaptively subtracting, within the time-space domain, the inverse transformed plurality of predicted multiple diffractions from the previously acquired seismic data.