IP Library › Granted Patent US 12,631,778
Granted Patent B2
US 12,631,778 · App. 18/254,764 · Granted May 19, 2026

Method and system for determining a location of hydrocarbon reservoir within subterranean region

Inventors: Yuhan Sui (Beijing, CN); Yue Ma (Beijing, CN); Hongwei Liu (Dhahran, SA); Dongliang Zhang (Khobar, SA)
Assignee: SAUDI ARABIAN OIL COMPANY
G01V1/345G01V1/282G01V1/301E21B7/04E21B2200/20G01V2210/645
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,631,778
App. No.
18/254,764
Granted
May 19, 2026
Kind
B2
Abstract

A system and methods are disclosed. The method includes obtaining a seismic dataset including a plurality of recorded multiple events, generating a predicted multiple model using a multiple prediction method and the seismic dataset, and estimating a set of initial matching filters using a matching method, to match the plurality of estimated and recorded multiple events. The method further includes generating a tensor field based on the predicted multiple model, determining a set of structure-oriented matching filters based on the set of initial matching filters and the tensor field, generating a filtered multiple model based on the predicted multiple model and the set of structure-oriented matching filters, and generating a multiple-attenuated seismic dataset based on the filtered multiple model and the seismic dataset, forming a seismic image based, at least in part, on the multiple-attenuated seismic dataset, and determining a location of a hydrocarbon reservoir based on the seismic image.

Claims (51)

1 . A method, comprising:

obtaining a seismic dataset for a subterranean region, wherein the seismic dataset comprises a plurality of recorded multiple events;

generating a predicted multiple model using a multiple prediction method and the seismic dataset, wherein the predicted multiple model comprises a plurality of estimated multiple events;

estimating a set of initial matching filters using a matching method, wherein the set of initial matching filters matches the plurality of estimated multiple events to the plurality of recorded multiple events;

generating a tensor field based, at least in part, on the predicted multiple model;

determining a set of structure-oriented matching filters based, at least in part, on the set of initial matching filters and the tensor field;

generating a filtered multiple model based, at least in part, on the predicted multiple model and the set of structure-oriented matching filters;

generating a multiple-attenuated seismic dataset based, at least in part, on the filtered multiple model and the seismic dataset;

forming a seismic image based, at least in part, on the multiple-attenuated seismic dataset; and

determining a location of a hydrocarbon reservoir within the subterranean region based, at least in part, on the seismic image.

2 . The method of claim 1 , further comprising planning a well to penetrate the hydrocarbon reservoir based on the location, wherein the planned well comprises a planned wellbore path.

3 . The method of claim 2 , further comprising drilling a wellbore guided by the planned wellbore path.

4 . The method of claim 1 , wherein the predicted multiple model comprises interbed multiple energy.

5 . The method of claim 1 , wherein the matching method comprises a least-square matching method.

6 . The method of claim 1 , wherein the multiple prediction method comprises a Marchenko-based multiple prediction method.

7 . The method of claim 1 , wherein determining the set of structure-oriented matching filters comprises solving an anisotropic diffusion equation.

8 . A non-transitory computer-readable medium having computer-executable instructions stored thereon that, when executed by a processor, perform steps comprising:

receiving a seismic dataset for a subterranean region, wherein the seismic dataset comprises a plurality of recorded multiple events;

generating a predicted multiple model using a multiple prediction method and the seismic dataset, wherein the predicted multiple model comprises a plurality of estimated multiple events;

estimating a set of initial matching filters using a matching method, wherein the set of initial matching filters matches the plurality of estimated multiple events to the plurality of recorded multiple events;

generating a tensor field based, at least in part, on the predicted multiple model;

determining a set of structure-oriented matching filters based, at least in part, on the set of initial matching filters and the tensor field;

generating a filtered multiple model based, at least in part, on the predicted multiple model and the set of structure-oriented matching filters;

generating a multiple-attenuated seismic dataset based, at least in part, on the filtered multiple model and the seismic dataset;

forming a seismic image based, at least in part, on the multiple-attenuated seismic dataset; and

determining a location of a hydrocarbon reservoir within the subterranean region based, at least in part, on the seismic image.

9 . The non-transitory computer-readable medium of claim 8 , further comprising planning a well to penetrate the hydrocarbon reservoir based on the location, wherein the planned well comprises a planned wellbore path.

10 . The non-transitory computer-readable medium of claim 8 , wherein the predicted multiple model comprises interbed multiple energy.

11 . The non-transitory computer-readable medium of claim 8 , wherein the matching method comprises a least-square matching method.

12 . The non-transitory computer-readable medium of claim 8 , wherein the multiple prediction method comprises a Marchenko-based multiple prediction method.

13 . The non-transitory computer-readable medium of claim 8 , wherein determining the set of structure-oriented matching filters comprises solving an anisotropic diffusion equation.

14 . A system, comprising:

a seismic acquisition system configured to:

record a seismic dataset for a subterranean region, wherein the seismic dataset comprises a plurality of recorded multiple events;

a seismic processing system configured to:

receive the seismic dataset,

generate a predicted multiple model using a multiple prediction method and the seismic dataset, wherein the predicted multiple model comprises a plurality of estimated multiple events,

estimate a set of initial matching filters using a matching method, wherein the set of initial matching filters matches the plurality of estimated multiple events to the plurality of recorded multiple events,

generate a tensor field based, at least in part, on the predicted multiple model,

determine a set of structure-oriented matching filters based, at least in part, on the set of initial matching filters and the tensor field,

generate a filtered multiple model based, at least in part, on the predicted multiple model and the set of structure-oriented matching filters,

generate a multiple-attenuated seismic dataset based, at least in part, on the filtered multiple model and the seismic dataset, and

form a seismic image based, at least in part, on the multiple-attenuated seismic dataset; and

a seismic interpretation workstation configured to:

determine a location of a hydrocarbon reservoir within the subterranean region based, at least in part, on the seismic image.

15 . The system of claim 14 , further comprising a well planning system configured to plan a well to penetrate the hydrocarbon reservoir based on the location, wherein the planned well comprises a planned wellbore path.

16 . The system of claim 15 , further comprising a drilling system configured to drill a wellbore guided by the planned wellbore path.

17 . The system of claim 14 , wherein the predicted multiple model comprises interbed multiple energy.

18 . The system of claim 14 , wherein the matching method comprises a least-square matching method.

19 . The system of claim 14 , wherein the multiple prediction method comprises a Marchenko-based multiple prediction method.

20 . The system of claim 14 , wherein determining the set of structure-oriented matching filters comprises solving an anisotropic diffusion equation.

Assignments (3)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Oct 13, 2023
From: ARAMCO FAR EAST (BEIJING) BUSINESS SERVICES CO., LTD.
To: SAUDI ARABIAN OIL COMPANY
Reel/Frame 065238/0035 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Aug 25, 2023
From: SUI, YUHAN; MA, YUE
To: ARAMCO FAR EAST (BEIJING) BUSINESS SERVICES CO., LTD.
Reel/Frame 064707/0365 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Aug 25, 2023
From: LIU, HONGWEI; ZHANG, DONGLIANG
To: SAUDI ARABIAN OIL COMPANY
Reel/Frame 064707/0381 →
Continuity (1)
Related Publication 20240353584A1 · Oct 24, 2024
References Cited (21)
US 9971051B2 · Wu · 2018 [cited by applicant]
US 11391856B2 · Al Ibrahim · 2022 [cited by examiner]
US 11693140B2 · Tawil · 2023 [cited by examiner]
US 12282128B2 · Widatalla · 2025 [cited by examiner]
US 20130194893A1 · Nagarajappa · 2013 [cited by applicant]
US 20170123090A1 · Khalil et al. · 2017 [cited by applicant]
US 20230037886A1 · He et al. · 2023 [cited by applicant]
CN 105474048A · 2016 [cited by applicant]
CN 106932824A · 2017 [cited by applicant]
CN 113219535A · 2021 [cited by applicant]
CN 113687417A · 2021 [cited by applicant]
WO 2016154404A1 · 2016 [cited by applicant]
International Search Report issued in PCT/CN2023/074355 on Apr. 7, 2023 (5 pages). [cited by applicant]
Written Opinion of the International Searching Authority issued in PCT/CN2023/074355 on Apr. 7, 2023 (3 pages). [cited by applicant]
Y. Wang, “Multiple subtraction using an expanded multichannel matching filter,” Geophysics, vol. 68, No. 1, pp. 346-354, 2003 (9 pages). [cited by applicant]
S. Fomel, “Adaptive multiple subtraction using regularized nonstationary regression,” Geophysics, vol. 74, No. 1, pp. V25-V33, 2009 (9 pages). [cited by applicant]
D. Hale, “Structure-oriented smoothing and semblance.” CWP report, vol. 635, No. 635, 2009 (10 pages). [cited by applicant]
D.Hale, “Structure-oriented bilateral filtering of seismic images” SEG Technical Program Expanded Abstracts, 2011 (5 pages). [cited by applicant]
D. Donno, “Improving multiple removal using least-squares dip filters and independent component analysis” Geophysics, vol. 76, No. 5, pp. V91-V104, 2011 (14 pages). [cited by applicant]
Joost van der Neut et al. “Adaptive overburden elimination with the multidimensional Marchenko equation” Geophysics, vol. 81, No. 5, P. T265-284, Sep.-Oct. 2016 (14 pages). [cited by applicant]
Jakubowicz H. “Wave equation prediction and removal of interbed multiple” 68th Annual International Meeting, SEG, Expanded Abstracts, 1527-1530 (1998) (4 pages). [cited by applicant]