IP Library › Granted Patent US 12,422,582
Granted Patent B2
US 12,422,582 · App. 18/263,319 · Granted Sep 23, 2025

Updating subsurface structural maps with well-measured orientation data while preserving local geological structures

Inventors: Yue Ma (Beijing, CN); Xu Ji (Beijing, CN); Nasher Muqbel Albinhassan (Beijing, CN)
Assignee: SAUDI ARABIAN OIL COMPANY
G01V1/302E21B49/00G01V1/307E21B41/00G01V2210/61
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,422,582
App. No.
18/263,319
Granted
Sep 23, 2025
Kind
B2
Abstract

The methods may include obtaining a seismic dataset regarding a subsurface region of interest and obtaining a well log for each of multiple wellbores penetrating the subsurface region of interest. The methods may also include determining a geological surface from the seismic dataset, wherein the geological surface includes seismic-estimated orientation data estimated at multiple points on the geological surface. The methods may further include determining an intersection point for each of the multiple wellbores with the geological surface, wherein the intersection point includes well-measured orientation data. The methods may still further include generating an updated geological surface by updating the seismic-estimated orientation data at the multiple points on the geological surface based, at least in part, on the well-measured orientation data.

Claims (63)

1. A method comprising:

obtaining, using a seismic acquisition system, a seismic dataset regarding a subsurface region of interest;

obtaining, using a logging tool, a well log for each of a plurality of wellbores penetrating the subsurface region of interest;

determining a sequence of surfaces from the seismic dataset, wherein each surface among the sequence of surfaces represents a geological surface within the subsurface region of interest, wherein the sequence of surfaces is ordered by increasing depth within the subsurface region of interest; and

for each surface in order:

obtaining, using the logging tool, a partial well log for a partially-drilled wellbore within the subsurface region of interest,

associating a plurality of points to each surface,

determining seismic-estimated orientation data at the plurality of points based on the seismic dataset,

determining a plurality of intersection points among the plurality of points, wherein each of the plurality of intersection points represents where each of the plurality of wellbores and the partially-drilled wellbore intersects with the geological surface,

determining well-measured orientation data at the plurality of intersection points based on the well log for the plurality of wellbores and the partial well log,

generating an updated surface by determining a region of influence around each of the plurality of intersection points and updating the seismic-estimated orientation data at the plurality of points within the region of influence based on the well-measured orientation data,

determining, using a wellbore planning system, an updated wellbore trajectory for an undrilled portion of the partially-drilled wellbore based on the updated surface, and

drilling, using a drilling system, a portion of the undrilled portion of the partially drilled wellbore along the updated wellbore trajectory.

2. The method of claim 1 , wherein the seismic-estimated orientation data comprises a dip angle and an azimuth.

3. The method of claim 1 , wherein generating the updated surface comprises performing 2D cubic spline interpolation.

4. The method of claim 1 , wherein the region of influence is of elliptical shape centered at each of the plurality of intersection points.

5. The method of claim 1 , wherein generating the updated surface further comprises:

obtaining structure orientation data of a geological structure in the subsurface region of interest;

determining a plurality of local points by applying a first coordinate transformation to the plurality of points in the region of influence based on the structure orientation data;

determining rotated orientation data at the plurality of local points based, at least in part, on the well-measured orientation data; and

updating the seismic-estimated orientation data at the plurality of points in the region of influence by applying a second coordinate transformation to the rotated orientation data.

6. A non-transitory computer-readable medium comprising computer-executable instructions stored thereon that, when executed on a processor, cause the processor to perform:

obtaining, from a seismic acquisition system, a seismic dataset regarding a subsurface region of interest;

obtaining, from a logging tool, a well log for each of a plurality of wellbores penetrating the subsurface region of interest,

determining a sequence of surfaces from the seismic dataset, wherein each surface among the sequence of surfaces represents a geological surface within the subsurface region of interest, wherein the sequence of surfaces is ordered by increasing depth within the subsurface region of interest; and

for each surface in order:

obtaining, using the logging tool, a partial well log for a partially-drilled wellbore within the subsurface region of interest,

associating a plurality of points to each surface,

determining seismic-estimated orientation data at the plurality of points based on the seismic dataset, determining a plurality of intersection points among the plurality of points, wherein each of the plurality of intersection points represents where each of the plurality of wellbores and the partially-drilled wellbore intersects with the geological surface,

determining well-measured orientation data at the plurality of intersection points based on the well log for the plurality of wellbores and the partial well log,

generating an updated surface by determining a region of influence around each of the plurality of intersection points and updating the seismic-estimated orientation data at the plurality of points within the region of influence based on the well-measured orientation data, and

determining an updated wellbore trajectory for an undrilled portion of the partially-drilled wellbore based on the updated surface, wherein a drilling system is configured to drill a portion of the undrilled portion of the partially-drilled wellbore along the updated wellbore trajectory.

7. The non-transitory computer-readable medium of claim 6 , wherein the computer-executable instructions further cause the processor to perform:

generating the updated surface by performing 2D cubic spline interpolation.

8. The non-transitory computer-readable medium of claim 6 , wherein the computer-executable instructions further cause the processor to perform:

obtaining structure orientation data of a geological structure in the subsurface region of interest;

determining a plurality of local points by applying a first coordinate transformation to the plurality of points in the region of influence based on the structure orientation data;

determining rotated orientation data at the plurality of local points based, at least in part, on the well-measured orientation data; and

updating the seismic-estimated orientation data at the plurality of points in the region of influence by applying a second coordinate transformation to the rotated orientation data.

9. A system comprising:

a seismic acquisition system configured to record a seismic dataset regarding a subsurface region of interest;

a logging tool configured to record a well log for each of a plurality of wellbores penetrating the subsurface region of interest;

a seismic processor configured to:

receive, from the seismic acquisition system, the seismic dataset,

receive, from the logging tool, the well log for the plurality of wellbores, and

determine a sequence of surfaces from the seismic dataset, wherein each surface among the sequence of surfaces represents a geological surface within the subsurface region of interest, wherein the sequence of surfaces is ordered by increasing depth within the subsurface region of interest,

for each surface in order:

receive, from the logging tool, a partial well log for a partially-drilled wellbore within the subsurface region on interest;

associate a plurality of points to each surface;

determine seismic-estimated orientation data at the plurality of points based on the seismic dataset;

determine a plurality of intersection points among the plurality of points, wherein each of the plurality of intersection points represents where each of the plurality of wellbores and the partially-drilled wellbore intersects with the geological surface;

determine well-measured orientation data at the plurality of intersection points based on the well log for the plurality of wellbores and the partial well log; and

generate an updated surface by determining a region of influence around each of the plurality of intersection points and updating the seismic-estimated orientation data at the plurality of points within the region of influence based on the well-measured orientation data;

a wellbore planning system configured to determine an updated wellbore trajectory for an undrilled portion of the partially-drilled wellbore based on the updated surface; and

a drilling system configured to drill a portion of the undrilled portion of the partially-drilled wellbore along the updated wellbore trajectory.

10. The system of claim 9 , wherein the seismic-estimated orientation data comprises a dip angle and an azimuth.

11. The system of claim 9 , wherein the seismic processor is further configured to generate the updated surface by performing 2D cubic spline interpolation.

12. The system of claim 9 , wherein the region of influence is of elliptical shape centered at each of the plurality of intersection points.

13. The system of claim 9 , wherein the seismic processor is further configured to generate the updated surface by:

obtaining structure orientation data of a geological structure in the subsurface region of interest;

determining a plurality of local points by applying a first coordinate transformation to the plurality of points in the region of influence based on the structure orientation data;

determining rotated orientation data at the plurality of local points based, at least in part, on the well-measured orientation data; and

updating the seismic-estimated orientation data at the plurality of points in the region of influence by applying a second coordinate transformation to the rotated orientation data.

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 Sep 5, 2023
From: MA, YUE
To: ARAMCO FAR EAST (BEIJING) BUSINESS SERVICES CO., LTD.
Reel/Frame 064794/0113 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Sep 5, 2023
From: JI, XU; ALBINHASSAN, NASHER MUQBEL
To: SAUDI ARABIAN OIL COMPANY
Reel/Frame 064794/0142 →
Continuity (1)
Related Publication 20240427041A1 · Dec 26, 2024
References Cited (37)
US 7557581B2 · Ostermeier · 2009 [cited by examiner]
US 8095318B2 · Heliot · 2012 [cited by examiner]
US 8931580B2 · Cheng et al. · 2015 [cited by applicant]
US 9581710B2 · Leiceaga · 2017 [cited by examiner]
US 10877171B2 · Ramsay · 2020 [cited by examiner]
US 20080236270A1 · Denichou et al. · 2008 [cited by applicant]
US 20090157361A1 · Toghi · 2009 [cited by examiner]
US 20100149917A1 · Imhof · 2010 [cited by examiner]
US 20130085676A1 · Sonneland et al. · 2013 [cited by applicant]
US 20160320512A1 · Zhao et al. · 2016 [cited by applicant]
US 20160364508A1 · Glazkova · 2016 [cited by examiner]
US 20200011167A1 · Zhao et al. · 2020 [cited by applicant]
US 20200033501A1 · Nyrnes et al. · 2020 [cited by applicant]
US 20200300064A1 · Gee et al. · 2020 [cited by applicant]
US 20220129788A1 · Zhang et al. · 2022 [cited by applicant]
CA 2804354A1 · 2014 [cited by applicant]
CN 102253411A · 2011 [cited by applicant]
CN 104520734A · 2015 [cited by applicant]
CN 106934858A · 2017 [cited by applicant]
CN 109267996A · 2019 [cited by applicant]
CN 112987086A · 2021 [cited by applicant]
CN 114562244A · 2022 [cited by applicant]
EP 0796442B1 · 2001 [cited by applicant]
WO 2009126375A1 · 2009 [cited by applicant]
WO 2013164685A2 · 2013 [cited by applicant]
WO 2018156354A1 · 2018 [cited by applicant]
WO 2020222050A1 · 2020 [cited by applicant]
WO 2022050967A1 · 2022 [cited by applicant]
Y. Ma, et al., “Automatic integration of 3D reflection seismic and well-measured orientation in sedimentary basins”, First International Meeting for Applied Geoscience & Energy, 2021, pp. 1131-1135 (5 pages). [cited by applicant]
F. Yong et al., “Estimation of lateral correlation length from deep seismic reflection profile based on stochastic model”, Acta Geophysica, 2021, vol. 69, pp. 1297-1312 (16 pages). [cited by applicant]
International Search Report and Written Opinion of the International Searching Authority issued in corresponding International Application No. PCT/CN2022/107763, mailed Apr. 26, 2023 (9 pages). [cited by applicant]
International Search Report and Written Opinion of the International Searching Authority issued in corresponding International Application No. PCT/CN2023/074354, mailed May 23, 2023 (9 pages). [cited by applicant]
Stewart, S.A., 2018. Generalization and multiscale structure of subsurface structural maps. Interpretation 6, T1045-T1054 (10 pages). [cited by applicant]
Stewart, S.A., 2020, Scale dependence of strike and dip in sedimentary basins: Implications for field measurements and integrating subsurface datasets, Journal of Structural Geology, 131, 103943 (8 pages). [cited by applicant]
Thore, P., Shtuka, A., Lecour, M., Ait-Ettajer, T., Cognot, R., 2002. Structural uncertainties: Determination, management, and applications. Geophysics 67, 840-852 (13 pages). [cited by applicant]
V. Tschannen et al; “Extracting Horizon Surfaces from 3D Seismic Data using Deep Learning”; Geophysics, vol. 85, No. 3, pp. N17-N26; May-Jun. 2020 (10 pages). [cited by applicant]
Non-Final Office Action issued by U.S. Patent Office for corresponding U.S. Appl. No. 18/008,081, mailed Dec. 6, 2024 (21 pages). [cited by applicant]