EP3811123B1 - Procédé de détermination d'une composition réelle de sous-sol - Google Patents

Procédé de détermination d'une composition réelle de sous-sol Download PDF

Info

Publication number
EP3811123B1
EP3811123B1 EP18759377.7A EP18759377A EP3811123B1 EP 3811123 B1 EP3811123 B1 EP 3811123B1 EP 18759377 A EP18759377 A EP 18759377A EP 3811123 B1 EP3811123 B1 EP 3811123B1
Authority
EP
European Patent Office
Prior art keywords
constraint
geological formation
geological
satisfied
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
EP18759377.7A
Other languages
German (de)
English (en)
Other versions
EP3811123A1 (fr
Inventor
Gérard MASSONNAT
David LEDEZ
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
TotalEnergies Onetech SAS
Original Assignee
TotalEnergies Onetech SAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by TotalEnergies Onetech SAS filed Critical TotalEnergies Onetech SAS
Publication of EP3811123A1 publication Critical patent/EP3811123A1/fr
Application granted granted Critical
Publication of EP3811123B1 publication Critical patent/EP3811123B1/fr
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V20/00—Geomodelling in general

Definitions

  • the present invention relates to determination of subsoil compositions/property, especially for the use in the hydrocarbon industries.
  • the accurate determination of the subsoil composition and structure is a key feature for determining hydrocarbon reservoirs and enabling a proper industrial extraction of hydrocarbons.
  • log data may be transformed with blocking methods to adapt the precision of the log data to the dimension of the cells.
  • US2010332205 discloses a model for modeling and determining real subsoil composition.
  • the invention relates to a method, a non-transitory computer readable storage medium and a device for determination of real subsoil composition according to claims 1, 5 and 6.
  • a last step the method is to determine real subsoil composition based on the meshed formation and output said result for future use in geophysical tools (e.g. prevision of hydrocarbon production, determination of a correct location of a well to be drilled, estimation of the reservoir capacity).
  • geophysical tools e.g. prevision of hydrocarbon production, determination of a correct location of a well to be drilled, estimation of the reservoir capacity.
  • the determination of the at least one parametric surface may be based on NURBS curve or NURBS surface or NURBS volume.
  • the determination of the fluvial trajectory may be based on a stochastic process.
  • the meshing of the fluvial zone may be based on a method in a group comprising Quadtree meshing method, the Octree meshing method, Front method, Delaunay method, Prograding grid, divergent grid, and aggrading grid.
  • the method may further comprise:
  • the method may further comprise:
  • the distortion of the parametric surface comprises a prior step of determining a point in said surface minimizing a distance between a distance of said point and a position of said constraint.
  • the geological constraints may be based on a constraint in a group comprising a data log constraint and a seismic constraint.
  • the constraint may be satisfied when a position of said constraint is inside the geological formation.
  • Another aspect relates to a computer program product comprising a computer-readable medium, having thereon a computer program comprising program instructions.
  • the computer program is loadable into a data-processing unit and adapted to cause the data-processing unit to carry out the method(s) described above when the computer program is run by the data-processing unit.
  • fluvial geological formations are described to exemplify the invention but it applies to any possible geological formations (e.g. lobes, turbiditic systems, etc.)
  • Figure 1 is a chart describing a possible process of the invention.
  • the manipulated model is possibly a parametric model of the subsoil.
  • a parametric model M p (u,v,t) of the subsoil is a transformation of a 3D model M(x,y,z) of the subsoil.
  • a 3D model M(x,y,z) of a real subsoil describes the subsoil according to its real geographical coordinates (x,y,z) (i.e. at the present time).
  • a parametric model M p (u,v,t) of said subsoil describes the state of the subsoil at a geological time t: each layer represents the state of the subsoil at the time t where the sedimentation occurs.
  • the parametric model M p (u,v,t) restore the horizontal layer for a given sedimentation/geological time t.
  • Said geological formation determination of a fluvial zone may comprise, as described in WO 2012/072923 , the displacement of particles (step 101) in the manipulated model by superimposing:
  • the stochastic term can comprise the superposition of a meandriform term and of a random perturbation. So doing, the modeling of the channel is rendered more realistic.
  • the meandriform term can comprise a superposition of at least one trigonometric function. Such a representation is realistic for a meandriform term, and easily parametrizable.
  • the observation data can comprise at least one of the following types of data:
  • model does not need to be meshed up to this point.
  • a parametric description of the trajectory may be easily determined in the model.
  • NURBS or non-uniform Rational B-Splines.
  • Non-uniform rational basis spline is a mathematical model used in computer graphics for generating and representing curves and surfaces. It offers great flexibility and precision for handling both analytic (surfaces defined by common mathematical formulae) and modeled shapes.
  • trajectory 201 is perpendicular to said view.
  • trajectory 201 is represented as a single point.
  • point 201 one means the point corresponding to the intersection between the view plan of figure and the trajectory 201.
  • NURBS 310, 311 and 312 To represent the extend of the determined geological formation, it is possible to represent the shape of bed by a plurality of NURBS 310, 311 and 312, ensuring that the contacts of these NURBS satisfy some parametric continuity (for instance C 0 , C 1 and/or C 2 ) and/or some geometric continuity (for instance G 0 , G 1 and/or G 2 ).
  • point 303 and 309 may be defined such that 303, 201 and 309 are aligned, the respective distance between 303-201 and 201-309 may be function of the local curvature of the trajectory 201 (see Figure 2b ).
  • point 303 represents a point on the convex side of the curvature of the trajectory 201
  • point 309 represents a point on the concave side of the curvature of the trajectory 201.
  • Point 306 may be determined such that the line (306;201) is perpendicular to the line (303;309) and such that the distance 306-201 is either a predetermined value or function of the local curvature or function of the distance 303-309.
  • points 306, 303 and 309 are determined, it is possible to determine a plurality of set of points 304, 305 307 and 308. These points may have a location function of the positions of points 306, 303 and 309 and/or the curvature of trajectory 201.
  • Points 303, 304, 305, 306 may define a first NURBS 311.
  • Points 309, 308, 307, 306 may define a second NURBS 310.
  • a third NURBS 312 may be defined to close the shape of the bed.
  • NURBS curves may be determined thanks to the Cox-de Boor's Algorithm.
  • NURBS surfaces may be defined thanks to the NURBS curves (or thanks to points that have been used to define the NURBS curves).
  • NURBS surfaces are very effective as it is very simple to "substract" volumes of determined geological formations. For instance, referring to Figure 3b , if two geological formations 320 and 321 are determined and if their respective volumes intersect each other, available NURBS intersection algorithms provide accurate and effective method to subtract the volume of formation 321 to formation 320 (representing the replacement of the formation 320 by formation 321 due to hydrodynamic/fluvial erosion) to create the new volume 322 ( Figure 3c ). Algorithms based on cells / meshes (i.e. identifying cells that are in both formations) are not as accurate and effective as the NURBS algorithms.
  • the model has some constraints (e.g. a seismic or geological indication that a formation is present at a given location), it is possible to distort the closest NURBS 320 to ensure that this constraint 401 (see Figure 4a , which is a top view of the model) is satisfied.
  • the distortion methods for NURBS surfaces are well available. Therefore it is possible to satisfy all constraints of the model.
  • a plurality of meshing is possible as described in Figure 5 .
  • meshes 501, 502, 503 or 104 are possible.
  • Well known methods such as the Quadtree meshing method, the Octree meshing method, Front method, Delaunay method, Prograding grid, divergent grid, aggrading grid, etc. are possible.
  • shape of other geological formations may be determined based on a distance to a previously determined formation or/and to a probability of existence of a surface of said other geological formations (function, for instance, of a distance (e.g. radial or lateral) to previously determined formations).
  • This association may be based on:
  • the quality index distributions are function of the distance to the vertical axis 201-306 of the shape of the bed and of the distance to the base 306 of the bed, but other geological formations may have quality index distributions function of other parameters set by the operator (i.e. the person setting the modeling).
  • a cell quality index QI c is computed as being ( i being the current driver, N being the total number of drivers).
  • the following process may be used for associating a facies to a cell in a formation (having a plurality of cells):
  • This association may be based on property quality index distributions QI i ( p ) for each driver i function of the considered (absolute or relative) position p (element 713 and 714 being respectively the property quality index distribution function of the distance to the vertical axis 201-306 of the shape of the bed and the property quality index distribution function of the distance to the base 306 of the bed) and on property distribution (curve 711).
  • the property quality index distributions are function of the distance to the vertical axis 201-306 of the shape of the bed and of the distance to the base 306 of the bed but other geological formations may have property quality index distributions function of other parameters set by the operator (i.e. the person setting the modeling).
  • a cell quality index QI c is computed as being ( i being the current driver, N being the total number of drivers).
  • the following process may be used for associating a property to a cell in a formation (having a plurality of cells):
  • step 105 of Figure 1 For part(s) of the model that is/are not determined geological formations (i.e. background zone), it is possible to mesh (step 105 of Figure 1 ) it/them with common method(s) and to associate its/their cells with property/facies (step 106 of Figure 1 ) according to the relative position of the cells in regard of the geological formation determined (e.g. facies probabilities may be set by the operator according to a radial/lateral distance of a determined geological formation).
  • facies probabilities may be set by the operator according to a radial/lateral distance of a determined geological formation.
  • Figure 8 is a possible embodiment for a device that enables the present invention.
  • the device 800 comprise a computer, this computer comprising a memory 805 to store program instructions loadable into a circuit and adapted to cause circuit 804 to carry out the steps of the present invention when the program instructions are run by the circuit 804.
  • the memory 805 may also store data and useful information for carrying the steps of the present invention as described above.
  • the circuit 804 may be for instance:
  • This computer comprises an input interface 803 for the reception of data/model/input used for the above method according to the invention and an output interface 806 for providing a complete model.
  • a screen 801 and a keyboard 802 may be provided and connected to the computer circuit 804.

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Claims (6)

  1. Procédé de détermination d'une composition ou d'une structure de sous-sol réel comprenant
    - la réception d'un modèle représentant le sous-sol réel, ledit modèle comprenant au moins une surface paramétrique décrivant une formation géologique dans ledit modèle ; caractérisé en ce que le procédé comprend :
    - la réception d'une contrainte géologique devant être satisfaite par ladite formation géologique ;
    - la détermination si la contrainte est satisfaite par ladite formation géologique ;
    - si la contrainte n'est pas satisfaite par ladite formation géologique, la déformation de la surface paramétrique de sorte que la contrainte soit satisfaite par ladite formation géologique ; dans lequel la déformation de la surface paramétrique comprend une étape préalable de détermination d'un point dans ladite surface minimisant une distance entre une distance dudit point et une position de ladite contrainte ;
    - la détermination d'une composition de sous-sol réel sur la base de la surface paramétrique déformée ; et
    - la production d'une composition de sous-sol réel pour une utilisation future dans des outils de géophysique.
  2. Procédé selon la revendication 1, dans lequel la surface paramétrique est basée sur une courbe NURBS ou une surface NURBS ou un volume NURBS.
  3. Procédé selon l'une des revendications précédentes, dans lequel la contrainte géologique est basée sur une contrainte dans un groupe comprenant une contrainte de registre de données et une contrainte sismique.
  4. Procédé selon l'une des revendications précédentes, dans lequel la contrainte est satisfaite lorsqu'une position de ladite contrainte se trouve à l'intérieur de la formation géologique.
  5. Support de stockage non transitoire lisible par ordinateur (805), dans lequel est stocké un programme informatique comprenant des instructions de programme, le programme informatique pouvant être chargé dans une unité de traitement de données (804) et adapté pour amener l'unité de traitement de données à effectuer les étapes de l'une quelconque des revendications 1 à 4 lorsque le programme informatique est exécuté par l'unité de traitement de données.
  6. Dispositif de détermination d'une composition ou d'une structure de sous-sol réel comprenant
    - une interface (803) pour la réception d'un modèle représentant le sous-sol réel, ledit modèle comprenant au moins une surface paramétrique décrivant une formation géologique dans ledit modèle ; caractérisé en ce que le dispositif comprend un circuit (804) pour :
    - recevoir une contrainte géologique devant être satisfaite par ladite formation géologique ;
    - déterminer si la contrainte est satisfaite par ladite formation géologique ;
    - si la contrainte n'est pas satisfaite par ladite formation géologique, déformer la surface paramétrique de sorte que la contrainte soit satisfaite par ladite formation géologique, dans lequel la déformation de la surface paramétrique comprend une étape préalable de détermination d'un point dans ladite surface minimisant une distance entre une distance dudit point et une position de ladite contrainte ;
    - déterminer une composition de sous-sol réel sur la base de la surface paramétrique déformée ; et
    - une interface de sortie (806) pour produire ladite composition de sous-sol réel pour une utilisation future dans des outils de géophysique.
EP18759377.7A 2018-06-20 2018-06-20 Procédé de détermination d'une composition réelle de sous-sol Active EP3811123B1 (fr)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/IB2018/000876 WO2019243858A1 (fr) 2018-06-20 2018-06-20 Procédé de détermination d'une composition de sous-sol réel

Publications (2)

Publication Number Publication Date
EP3811123A1 EP3811123A1 (fr) 2021-04-28
EP3811123B1 true EP3811123B1 (fr) 2023-11-29

Family

ID=63364096

Family Applications (1)

Application Number Title Priority Date Filing Date
EP18759377.7A Active EP3811123B1 (fr) 2018-06-20 2018-06-20 Procédé de détermination d'une composition réelle de sous-sol

Country Status (2)

Country Link
EP (1) EP3811123B1 (fr)
WO (1) WO2019243858A1 (fr)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2868848B1 (fr) * 2004-04-07 2006-06-23 Earth Decision Sciences Sa Procede, dispositif et produit-programme de simulation stochastique de paleo chenaux ou de paleo structures geologiques
FR2947345B1 (fr) * 2009-06-26 2011-07-15 Inst Francais Du Petrole Methode pour modifier des proportions de facies lors du calage d'historique d'un modele geologique
FR2968443A1 (fr) 2010-12-03 2012-06-08 Total Sa Procede de simulation de formation geologique d'une zone fluviale
WO2016070073A1 (fr) * 2014-10-31 2016-05-06 Exxonmobil Upstream Research Company Gestion des discontinuités dans des modèles géologiques

Also Published As

Publication number Publication date
EP3811123A1 (fr) 2021-04-28
WO2019243858A1 (fr) 2019-12-26

Similar Documents

Publication Publication Date Title
US10598819B2 (en) Systems and methods to build sedimentary attributes
US10795053B2 (en) Systems and methods of multi-scale meshing for geologic time modeling
US9053570B1 (en) Systems and processes for building multiple equiprobable coherent geometrical models of the subsurface
EP3293552B1 (fr) Système et procédé pour l'édition de modèles géologiques par commutation entre des modèles basés sur le volume et des modèles structuraux basés sur la surface augmentée par des faisceaux de fibres stratigraphiques
US11555944B2 (en) Method for determination of subsoil composition
AU2016223281B2 (en) Method for parameterizing a 3D domain with discontinuities
EP3811123B1 (fr) Procédé de détermination d'une composition réelle de sous-sol
US11656380B2 (en) Method for determination of real subsoil geological formation
US11614556B2 (en) Method for determination of real subsoil geological formation
US20210270997A1 (en) A method for determination of real subsoil composition
US10310116B2 (en) Method for determining a stacked model describing architectural elements
US11630235B2 (en) Method for determination of real subsoil geological formation
US11995382B2 (en) Method for determination of real subsoil geological formation
US11487034B2 (en) Method for determination of real subsoil geological formation

Legal Events

Date Code Title Description
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: UNKNOWN

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE

PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20210113

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

DAV Request for validation of the european patent (deleted)
DAX Request for extension of the european patent (deleted)
RAP1 Party data changed (applicant data changed or rights of an application transferred)

Owner name: TOTALENERGIES ONETECH

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

17Q First examination report despatched

Effective date: 20220818

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: GRANT OF PATENT IS INTENDED

INTG Intention to grant announced

Effective date: 20230821

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE PATENT HAS BEEN GRANTED

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602018061841

Country of ref document: DE

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: NO

Ref legal event code: T2

Effective date: 20231129

REG Reference to a national code

Ref country code: NL

Ref legal event code: FP

P01 Opt-out of the competence of the unified patent court (upc) registered

Effective date: 20231214

REG Reference to a national code

Ref country code: LT

Ref legal event code: MG9D

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20240301

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20240329

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20231129

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: ES

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20231129

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20231129

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20240329

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20240301

Ref country code: ES

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20231129

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20240229

REG Reference to a national code

Ref country code: AT

Ref legal event code: MK05

Ref document number: 1636703

Country of ref document: AT

Kind code of ref document: T

Effective date: 20231129