US 12,436,305 B2
Automatic tying structure maps of subsurface horizons to well-derived orientation information
Yue Ma, Beijing (CN); Xu Ji, Dhahran (SA); Yubing Li, Beijing (CN); and Yi Luo, Dhahran (SA)
Assigned to SAUDI ARABIAN OIL COMPANY, Dhahran (SA)
Appl. No. 18/008,081
Filed by SAUDI ARABIAN OIL COMPANY, Dhahran (SA)
PCT Filed Jul. 26, 2022, PCT No. PCT/CN2022/107763
§ 371(c)(1), (2) Date Dec. 2, 2022,
PCT Pub. No. WO2024/020763, PCT Pub. Date Feb. 1, 2024.
Prior Publication US 2024/0264324 A1, Aug. 8, 2024
Int. Cl. G01V 1/30 (2006.01); G01V 11/00 (2006.01)
CPC G01V 1/302 (2013.01) [G01V 11/00 (2013.01); G01V 2210/643 (2013.01)] 12 Claims
OG exemplary drawing
 
1. A method, comprising:
obtaining, using a seismic acquisition system, a seismic image volume for a subsurface region of interest;
obtaining, using a logging system, a well log for each of a plurality of wellbores penetrating the subsurface region of interest;
determining, using a computer system, a seismic map of a geological surface from the seismic image volume, wherein the seismic map comprises an estimated depth and an estimated vector normal to the seismic map at a plurality of horizontal locations;
determining, using the computer system, an intersection point for each of the plurality of wellbores with the geological surface, wherein the intersection point comprises a depth, a horizontal location, and a strike and a dip of the geological surface;
forming, using the computer system, a cost function based, at least in part, on the seismic map and the intersection points of the plurality of wellbores; and
constructing, using the computer system, a subsurface map by solving a constrained optimization problem based on the cost function;
determining, using the computer system and based, at least in part, on the constructed subsurface map, a wellbore path configured to penetrate a hydrocarbon reservoir; and
drilling, using a drilling system, along the wellbore path,
wherein forming the cost function further comprises:
determining a warped map of the geological surface based, at least in part, on the seismic map and the depth of the intersection point for each of wellbore with the geological surface; and
forming the cost function based, at least in part, on a square of a difference between the estimated depth of the seismic map and a depth of the warped map summed over a plurality of horizontal location, and on the vector normal to the seismic map and the strike and dip of the geological surface at the intersection point for each of the plurality of wellbores.