IP Library Granted Patent US 12,704,652
Granted Patent B2
US 12,704,652 · App. 17/951,250 · Granted Aug 11, 2026

Faulted seismic horizon mapping

Inventors: Daniel James Possee (Southampton, GB); Graham Baines (Abingdon, GB)
Assignee: Landmark Graphics Corporation
G01V1/345E21B49/00G01V1/30G01V2210/642
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,704,652
App. No.
17/951,250
Filed
Sep 23, 2022
Granted
Aug 11, 2026
Kind
B2
Examiner
HAGOS, EYOB
Art Unit
2857
USPC
702/16
Abstract

Disclosed herein are embodiments of a method, a non-transitory computer readable medium, and an apparatus for faulted seismic horizon mapping. In one example, a method comprises: obtaining seismic data for a seismic volume that corresponds to a subsurface formation; generating a map of at least one horizon in the subsurface formation based on the seismic volume; identifying at least one fault intersecting the at least one horizon; determining a throw of the at least one fault; and updating the map of the at least one horizon to incorporate the at least one fault based on the throw of the at least one fault.

Claims (77)

1 . A method comprising:

obtaining, via one or more seismic sensors, seismic data for a seismic volume that corresponds to a subsurface formation;

generating a model of at least one horizon in the subsurface formation based on the seismic data for the seismic volume;

identifying at least one fault intersecting the at least one horizon;

determining an estimated throw of the at least one fault;

iteratively updating the model of the at least one horizon to incorporate the at least one fault as a part of the model by substituting equations solving for points on a surface of the horizon that lie on either side of the fault, based on the estimated throw of the at least one fault and an update of the estimated throw of the at least one fault, until a convergence criterion is satisfied; and

performing at least one downhole operation for recovery of hydrocarbons based on the updated model of the at least one horizon, wherein the at least one operation for the recovery of hydrocarbons is one of completion of a wellbore, updating drilling operations, perforating, fracking, logging operations, and additional sampling of the subsurface formation, and wherein updating drilling operations includes changing a drilling operation, a perforating operation, or a fracking operation.

2 . The method of claim 1 , wherein generating the model of the at least one horizon comprises:

generating at least one constraint condition for generating the model;

generating a number of grid points along the at least one horizon, based on at the at least one constraint, and wherein each of the number of grid points is derived from a seismic trace;

determine a weighting factor for the at least one horizon on local grid edges of the number of grid points; and

determining a connecting line between adjacent grid points of the number of grid points indicating a seismic attribute.

3 . The method of claim 2 , wherein the seismic attribute comprises a seismic dip.

4 . The method of claim 2 , wherein updating the model of the at least one horizon comprises:

removing any grid point on the at least one fault; and

removing any connecting lines that connects any grid point that has been removed.

5 . The method of claim 2 , wherein updating the model of the at least one horizon comprises removing any connecting line that crosses the at least one fault.

6 . The method of claim 1 , wherein updating the model of the at least one horizon comprises determining, independently, the model of the at least one horizon on each side of the at least one fault.

7 . The method of claim 2 , wherein updating the model of the at least one horizon comprises determining a connecting line between two grid points of the number of grid points on opposite sides of the fault relative to the seismic attribute.

8 . The method of claim 7 , wherein the two grid points includes grid points nearest the at least one fault on each side of the at least one fault.

9 . The method of claim 2 ,

wherein the at least one horizon includes a first horizon partition and a second horizon partition, and

wherein updating the model of the at least one horizon comprises,

adding a first new grid point positioned on the at least one fault at a location wherein the first horizon partition is predicted to intersect with the at least one fault;

adding a second new grid point positioned on the at least one fault at a location wherein the second horizon partition is predicted to intersect with the at least one fault;

determining a connecting line between the first new grid point and a grid point of the number of grid points that is nearest the at least one fault in the first horizon partition;

determining a connecting line between the second new grid point and a grid point of the number of grid points that is nearest the at least one fault in the second horizon partition; and

determining a connecting line between the first new grid point and the second new grid point.

10 . The method of claim 1 , wherein changing the drilling operation, the perforating operation, or the fracking operation comprises initiating or stopping the drilling operation, the perforating operation, or the fracking operation.

11 . A non-transitory, machine-readable medium having instructions stored thereon that are executable by a processor to perform seismic data processing operations comprising:

obtaining, via one or more seismic sensors, seismic data for a seismic volume that corresponds to a subsurface formation;

generating a model of at least one horizon in the subsurface formation based on the seismic data for the seismic volume;

identifying at least one fault intersecting the at least one horizon;

determining an estimated throw of the at least one fault;

iteratively updating the model of the at least one horizon to incorporate the at least one fault as part of the model by substituting equations solving for points on a surface of the horizon that lie on either side of the fault, based on the estimated throw of the at least one fault and an update of the estimated throw of the at least one fault, until a convergence criterion is satisfied; and

performing an operation for recovery of hydrocarbons based on the updated model of the at least one horizon, wherein the operation for the recovery of hydrocarbons is one of completion of a wellbore, updating drilling operations, perforating, fracking, logging operations, and additional sampling of the subsurface formation, and wherein updating drilling operations includes changing a drilling operation, a perforating operation, or a fracking operation.

12 . The non-transitory, machine-readable medium of claim 11 , wherein generating the model of the at least one horizon comprises:

generating at least one constraint condition for generating the model;

generating a number of grid points along the at least one horizon, based on the at least one constraint, and wherein each of the number of grid points is derived from a seismic trace;

determine a weighting factor for the at least one horizon on local grid edges of the number of grid points; and

determining a connecting line between adjacent grid points of the number of grid points indicating a seismic attribute, wherein the seismic attribute comprises a seismic dip.

13 . The non-transitory, machine-readable medium of claim 12 , wherein updating the model of the at least one horizon comprises:

removing any grid point on the at least one fault; and

removing any connecting lines that connects any grid point that has been removed.

14 . The non-transitory, machine-readable medium of claim 12 , wherein updating the model of the at least one horizon comprises removing any connecting line that crosses the at least one fault.

15 . The non-transitory, machine-readable medium of claim 12 , wherein updating the model of the at least one horizon comprises determining a connecting line between two grid points of the number of grid points on opposite sides of the fault relative to the seismic attribute, wherein the two grid points includes grid points nearest the at least one fault on each side of the at least one fault.

16 . The non-transitory, machine-readable medium of claim 12 ,

wherein the at least one horizon is separated into a first horizon partition and a second horizon partition by the at least one fault, and

wherein updating the model of the at least one horizon comprises,

adding a first new grid point positioned on the at least one fault at a location wherein the first horizon partition is predicted to intersect with the at least one fault;

adding a second new grid point positioned on the at least one fault at a location wherein the second horizon partition is predicted to intersect with the at least one fault;

determining a connecting line between the first new grid point and a grid point of the number of grid points that is nearest the at least one fault in the first horizon partition;

determining a connecting line between the second new grid point and a grid point of the number of grid points that is nearest the at least one fault in the second horizon partition; and

determining a connecting line between the first new grid point and the second new grid point.

17 . An apparatus for performing seismic data processing comprising:

a processor; and

a non-transitory computer-readable medium having instructions stored thereon that are executable by the processor to cause the processor to:

obtain, via one or more seismic sensors, seismic data for a seismic volume that corresponds to a subsurface formation;

generate a model of at least one horizon in the subsurface formation based on the seismic data for the seismic volume;

identify at least one fault intersecting the at least one horizon;

determine an estimated throw of the at least one fault;

iteratively update the model of the at least one horizon to incorporate the at least one fault as part of the model by substituting equations solving for points on a surface of the horizon that lie on either side of the fault, based on the estimated throw of the at least one fault and an update of the estimated throw of the at least one fault, until a convergence criterion is satisfied; and

direct a drilling tool to perform an operation for recovery of hydrocarbons based on the updated model of the at least one horizon, wherein the operation for the recovery of hydrocarbons is one of completion of a wellbore, updating drilling operations, perforating, fracking, logging operations, and additional sampling of the subsurface formation, and wherein updating drilling operations includes changing a drilling operation, a perforating operation, or a fracking operation.

18 . The apparatus of claim 17 , wherein the instructions executable by the processor to cause the processor to generate the model of the at least one horizon comprises instructions executable by the processor to cause the processor to:

generating at least one constraint condition for generating the model;

generate a number of grid points along the at least one horizon, based on the at least one constraint, and wherein each of the number of grid points is derived from a seismic trace;

determine a weighting factor for the at least one horizon on local grid edges of the number of grid points; and

determine a connecting line between adjacent grid points of the number of grid points indicating a seismic attribute.

19 . The apparatus of claim 18 , wherein the instructions executable by the processor to cause the processor to update the model of the at least one horizon comprises instructions executable by the processor to cause the processor to determine a connecting line between two grid points of the number of grid points on opposite sides of the fault relative to the seismic attribute, wherein the two grid points includes grid points nearest the at least one fault on each side of the at least one fault.

20 . The apparatus of claim 18 ,

wherein the at least one horizon includes a first horizon partition and a second horizon partition, and

wherein the instructions executable by the processor to cause the processor to update the model of the at least one horizon comprises instructions executable by the processor to cause the processor to:

add a first new grid point positioned on the at least one fault at a location wherein the first horizon partition is predicted to intersect with the at least one fault;

add a second new grid point positioned on the at least one fault at a location wherein the second horizon partition is predicted to intersect with the at least one fault;

determine a connecting line between the first new grid point and a grid point of the number of grid points that is nearest the at least one fault in the first horizon partition;

determine a connecting line between the second new grid point and a grid point of the number of grid points that is nearest the at least one fault in the second horizon partition; and

determine a connecting line between the first new grid point and the second new grid point.

Assignments (1)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Sep 23, 2022
From: POSSEE, DANIEL JAMES; BAINES, GRAHAM
To: LANDMARK GRAPHICS CORPORATION
Reel/Frame 061192/0365 →
Continuity (1)
Related Publication 20240111067A1 · Apr 4, 2024
References Cited (58)
US 6138076A · Graf · 2000 [cited by examiner]
US 6757615B2 · Zauderer et al. · 2004 [cited by applicant]
US 7283911B2 · Fitzsimmons et al. · 2007 [cited by applicant]
US 7523024B2 · Endres et al. · 2009 [cited by applicant]
US 8213261B2 · Imhof et al. · 2012 [cited by applicant]
US 8566069B2 · Pauget et al. · 2013 [cited by applicant]
US 9164192B2 · Maerten · 2015 [cited by examiner]
US 9759826B2 · Mallet · 2017 [cited by applicant]
US 9817142B2 · Bornhurst et al. · 2017 [cited by applicant]
US 9869783B2 · Jin et al. · 2018 [cited by applicant]
US 10215868B2 · Nguyen et al. · 2019 [cited by applicant]
US 10339709B1 · Mallet et al. · 2019 [cited by applicant]
US 11209561B2 · Klinger · 2021 [cited by examiner]
US 11320551B2 · Liu · 2022 [cited by examiner]
US 11520077B2 · Denli · 2022 [cited by examiner]
US 20030158669A1 · Davidson · 2003 [cited by applicant]
US 20060253759A1 · Wei · 2006 [cited by applicant]
US 20080015784A1 · Dorn et al. · 2008 [cited by applicant]
US 20100214870A1 · Pepper et al. · 2010 [cited by applicant]
US 20110118985A1 · Aarre · 2011 [cited by applicant]
US 20110307178A1 · Hoekstra · 2011 [cited by applicant]
US 20130235697A1 · Szydlik et al. · 2013 [cited by applicant]
US 20150066460A1 · Klinger · 2015 [cited by examiner]
US 20150103623A1 · Gersztenkorn · 2015 [cited by applicant]
US 20150293260A1 · Ghayour et al. · 2015 [cited by applicant]
US 20160124117A1 · Huang · 2016 [cited by examiner]
US 20170242141A1 · Cotton et al. · 2017 [cited by applicant]
US 20170337302A1 · Mezghani · 2017 [cited by examiner]
US 20180031721A1 · Etiene Queiroz · 2018 [cited by examiner]
US 20190243017A1 · Klinger · 2019 [cited by examiner]
US 20190302309A1 · Li · 2019 [cited by examiner]
US 20200166664A1 · Bin Gubair et al. · 2020 [cited by applicant]
US 20220091292A1 · Nguyen · 2022 [cited by applicant]
US 20220221604A1 · Skjæveland · 2022 [cited by examiner]
US 20220291418A1 · Noufal · 2022 [cited by examiner]
US 20230085023A1 · Possee et al. · 2023 [cited by applicant]
US 20230161061A1 · Denli et al. · 2023 [cited by applicant]
GB 2592203 · 2021 [cited by applicant]
KR 20200075727 · 2020 [cited by applicant]
WO 2018065684 · 2018 [cited by applicant]
WO WO2022256039A1 · 2022 [cited by examiner]
WO 2023042476 · 2023 [cited by applicant]
“U.S. Appl. No. 17/447,604, Final Office Action”, Dec. 21, 2022, 13 pages. [cited by applicant]
“PCT Application No. PCT/US2023/065971, International Search Report and Written Opinion”, Jul. 27, 2023, 10 pages. [cited by applicant]
“PCT Application No. PCT/US2022/077032, International Search Report and Written Opinion”, Jun. 9, 2023, 10 pages. [cited by applicant]
Wang, et al., “3D Seismic Horizon Extraction with Horizon Patch Constraints”, Paper presented at the 2015 SEG Annual Meeting, New Orleans, Louisiana, Oct. 2015, Oct. 18, 2015, 5 pages. [cited by applicant]
“PCT Application No. PCT/US2021/071473, International Search Report and Written Opinion”, Jun. 14, 2022, 8 pages. [cited by applicant]
“U.S. Appl. No. 17/447,604, Non-Final Office Action”, Sep. 6, 2022, 11 pages. [cited by applicant]
Branch, et al., “A Subspace, Interior, and Conjugate Gradient Method for Large-Scale Bound-Constrained Minimization Problems”, SIAM Journal on Scientific Computing, vol. 21, 23 pages, 1999. [cited by applicant]
Ferreira, et al., “Seismic Graph Analysis to Aid Seismic Interpretation”, Interpretation, vol. 7, Issue 3, May 28, 2019, 12 pages. [cited by applicant]
Hale, “Dynamic Warping of Seismic Images”, Geophysics, vol. 78, Issue 2, Mar. 1, 2013, pp. 233-244. [cited by applicant]
Hale, “Methods to Compute Fault Images, Extract Fault Surfaces, and Estimate Fault Throws From 3D Seismic Images”, Geophysics, vol. 78, Issue 2, Mar.-Apr. 2013, pp. O33-O43. [cited by applicant]
Höllt, et al., “Interactive Seismic Interpretation with Piecewise Global Energy Minimization”, 2011 IEEE Pacific Visualization Symposium, Hong Kong, Mar. 4, 2011, 8 pages. [cited by applicant]
Kidd, “Fundamentals of 3D Seismic Volume Visualization”, Offshore Technology Conference, Houston, TX, May 3-6, 1999, 13 pages. [cited by applicant]
Pauget, et al., “A Global Approach in Seismic Interpretation Based on Cost Function Minimization”, SEG Technical Program Expanded Abstracts, 2009, pp. 2592-2596. [cited by applicant]
Wu, et al., “Horizon vols. with interpreted constraints”, Geophysics, vol. 80, No. 2 (Mar.-Apr. 2015), 13 pages. [cited by applicant]
Wu, et al., “Least-squares horizons with local slopes and multi-grid correlations”, Geophysics, vol. 83; Jun. 2018, 15 pages. [cited by applicant]
“U.S. Appl. No. 17/950,995 Non-Final Office Action”, Feb. 3, 2025, 20 pages. [cited by applicant]