WO2012072923A2 - Procede de simulation de formation geologique d'une zone fluviale - Google Patents

Procede de simulation de formation geologique d'une zone fluviale Download PDF

Info

Publication number
WO2012072923A2
WO2012072923A2 PCT/FR2011/052769 FR2011052769W WO2012072923A2 WO 2012072923 A2 WO2012072923 A2 WO 2012072923A2 FR 2011052769 W FR2011052769 W FR 2011052769W WO 2012072923 A2 WO2012072923 A2 WO 2012072923A2
Authority
WO
WIPO (PCT)
Prior art keywords
zone
particles
term
simulating
occurrence
Prior art date
Application number
PCT/FR2011/052769
Other languages
English (en)
Inventor
Gérard MASSONNAT
Francis Morandini
Original Assignee
Total Sa
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 Total Sa filed Critical Total Sa
Priority to RU2013130263/12A priority Critical patent/RU2013130263A/ru
Priority to BR112013013491A priority patent/BR112013013491A2/pt
Priority to US13/991,360 priority patent/US10007741B2/en
Priority to NO20130806A priority patent/NO345728B1/no
Publication of WO2012072923A2 publication Critical patent/WO2012072923A2/fr

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V20/00—Geomodelling in general
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V99/00—Subject matter not provided for in other groups of this subclass
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F9/00—Arrangements for program control, e.g. control units
    • G06F9/06—Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/44—Arrangements for executing specific programs
    • G06F9/455—Emulation; Interpretation; Software simulation, e.g. virtualisation or emulation of application or operating system execution engines
    • G—PHYSICS
    • G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
    • G09B—EDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
    • G09B23/00—Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes
    • G09B23/40—Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes for geology
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/70—Other details related to processing

Definitions

  • a method of simulating the geological formation of a fluvial zone is provided.
  • the invention relates to the field of simulation of geological processes for the study of the subsoil. We are particularly interested in the formation phenomena of fluvial zones.
  • the sands are areas of interest for oil exploration. We therefore seek to predict their presence and their characteristics.
  • a method for simulating the geological formation of a fluvial zone comprising:
  • the term stochastic may include the superposition of a meandering term and a random disturbance. In doing so, channel modeling is made more realistic.
  • meandering can include a superposition of at least one trigonometric function. Such a representation is realistic for a meandering term, and easily parameterizable.
  • the observation data may comprise at least one of the following types of data:
  • the probability of the occurrence of a separation of a particle into two independent particles can be estimated and, in case of occurrence, proceed to such separation;
  • the probability of the occurrence of a disappearance of a particle can be estimated based on the observation data.
  • a channeling belt defining a width and a thickness in which the particles are allowed to move.
  • a computer program product for geological formation simulation of a fluvial zone intended to be stored in a memory of a central unit, and / or stored on a memory medium intended to cooperate with a reader of said central unit and / or downloaded via a telecommunication network, characterized in that it comprises instructions for implementing such a method when it is executed on a programmable machine.
  • FIG. 1 is a diagrammatic view of a channel
  • FIG. 2a is a schematic view from above of a channel model
  • FIG. 2b is a diagrammatic sectional view of the channel
  • FIG. 3 is an illustrative diagram of the zoning in a channel
  • FIG. 4 is an overall diagram illustrating particle trajectory modeling in a channel part
  • FIG. 5 is an illustrative diagram similar to that of FIG. 3 describing a separation process
  • FIG. 6 is a view similar to FIG. 3 describing a melting process
  • FIG. 7 is a view similar to FIG. 3 describing a disappearing process
  • Figure 8 is a schematic top view of a model in which is symbolized a displacement field imposed by an external source.
  • FIG. 1 schematically depicts a fluvial zone 1 extending between an upstream zone 2a and a downstream zone 2b.
  • the fluvial zone to be modeled is considered over a certain width 1, a certain thickness e and a certain length L. Thickness and width are for example variable along this length.
  • a river zone 1 generally comprises a channel
  • This channel was formed by the flow of particles over time from the upstream zone to the downstream zone.
  • a channel may have at least one of the following configurations; here being listed here in the order in which they are generally found between the upstream zone and the downstream zone:
  • a meandering system 4c in which the channels have a high sinuosity
  • a deltaic system 4e in which a main channel splits into a number of secondary channels opening into the sea, and
  • turbid system 4f related to the interaction with another flow, such as a maritime flow.
  • a spatial model of the fluvial zone can be realized in the form of a surface or volume mesh in which each of the elements represents a location of the fluvial zone.
  • observation data can for example include:
  • a gradient of flow velocity in the fluvial zone which can for example be defined on the basis of a hydraulic gradient in this zone, and / or from the altitude data of the upstream and downstream regions,
  • geometric parameters of the fluvial zone such as, for example, the current geometry of the surface of the zone, the discontinuities observed, etc.
  • geological data resulting from imagery for example seismic imagery carried out in the zone
  • FIG. 2a thus represents from above a spatial model of a region of the fluvial zone, presenting a set of mesh elements, for example parallelepipedal (squares in plan view), in which a channel model 5 has already been generated by simulation.
  • FIG. 2b represents a sectional view along the line B-B of FIG. 2a, thus showing a slice of the calculated channel model 5.
  • this channel has three of the distinct zones mentioned above in relation to FIG. 1, which are identified in FIG. 3 by references 6a, 6b and 6c.
  • Each of these zones can be the subject of a separate modeling, while respecting the constraints of connectivity at the interface between two zones. The modeling of one of these zones will be described hereinafter with reference to FIG.
  • this is zone 6a.
  • the channel 8 defines a width in which the particles are allowed to flow, this width possibly being variable along the average flow curve (the width is here represented constant).
  • the channeling belt may extend not only on the surface, but also to a certain depth below the representation plane of the window 7a.
  • the area shown comprises an upstream hydraulic load 10a and a downstream hydraulic load 10b.
  • the displacement of particles in the spatial model is simulated by superimposing a deterministic term and a stochastic term. For example, we write:
  • ⁇ represents the displacement of the particle
  • AD represents the deterministic term
  • ⁇ represents the stochastic term.
  • the deterministic term is defined from the flow field shown on window 7a of FIG.
  • the term stochastic is added to this term to define a perturbation of the deterministic term.
  • the stochastic term can itself be realized as the superposition of two terms.
  • the first of these terms is a meandering term that can easily be realized, from the mathematical point of view, by one or more trigonometric terms (sinusoidal) as represented on the window 7b of FIG. 4.
  • the curve 11 a period and a half of a sine extending between points a and g, and of a given amplitude.
  • the points a and g correspond to the upstream and downstream zones of the modeled zone, and the number of periods as well as the amplitude of the sinus can be parameterized independently for this zone.
  • the meandriform term is superimposed on a random perturbation represented on the window 7b by the arrows 12.
  • This random term is defined as a disturbance of the sinusoidal curve 11.
  • the orientation and / or the amplitude of the perturbation can be determined randomly on a plurality of points 13 of the sinusoidal curve 11, for example equidistant from each other, or of equidistant abscissa, as shown.
  • the window 7c represents a simulated trajectory 13 of particles, within the channeling belt 8, obtained from a superposition of the type mentioned above.
  • the parameterization (period, amplitude) of the sinusoidal function is determined according to the zone modeled for the fluvial zone.
  • the parameterization (frequency, intensity) of the random disturbance can be determined according to the zone considered. This parameterization can also be performed taking into account other observation data, among the observation data mentioned above.
  • Such divisions are visible at 15 in Figure 5.
  • the probability of having a division may depend on the area. In the example of FIG. 5, it is possible, for example, to provide that no division is possible in zone 6a. In zone 6b, it can be foreseen that divisions are possible, and that, in the case of division, a particle divides into a maximum number of daughter particles (here the maximum number of daughter particles is two). Furthermore, it is also possible to set the maximum number of divisions allowed for a particle in a given zone (at most two divisions in zone 6b shown in FIG. 5 in the example shown).
  • FIG. 6 Another phenomenon that can be modeled is shown in FIG. 6. If two particles are located at a location at one and the same time, it is a matter of junction or fusion of these particles. Such a merger is visible in FIG. 6 as represented by the references 16. It will be noted that the merger process can be implemented in parallel with the division process represented above, with divisions occurring on the example of Figure 6, in point 15.
  • the possibility of implementing a merger phenomenon is probabilized. This probability may for example depend on the area considered. For example, it can be expected that one proceeds systematically to a merger in the zone 6a. In zone 6b, as shown, it is possible to never implement fusion.
  • FIG. 7 Another phenomenon that can be modeled is shown in FIG. 7. It is, in this case, a phenomenon of disappearance of a particle. For example, it is possible to predict that a particle will disappear after having traveled a certain distance, or to probabilize the possibility of seeing a disappearance as a function of the path traveled by the particle. Such disappearances are visible, referenced by reference 17 in FIG. 7. Separation, junction and disappearance can be probabilized according to observation data.
  • Another phenomenon that can be modeled is the superposition of another energy source to the flow. in the river zone.
  • this other energy source is a littoral drift as represented by the arrow 18 in FIG. 8.
  • This term may for example be added to the deterministic term directly determined by the hydraulic gradient in the area when calculating the particle displacement.
  • the displacement of particles in the fluvial zone is simulated according to one and / or the other of the equations and phenomena presented below.
  • This simulation is for example implemented on a programmable machine, such as a computer, having stored in memory a computer program stored in memory locally or accessible via a network, and programmed to implement the above steps. It is possible to implement an end-of-simulation determination criterion, for example based on the number of simulated particles, on the simulation time, and / or on the similarity of certain results of the simulation with observation data.
  • the simulation can evaluate one or both of the following parameters:
  • the realistic modeling of channel formation, obtained above, and the above parameters may be useful in the context of the oil production of certain channel locations. For example, one can predict, on a channel sometimes several hundred kilometers long, areas likely to contain hydrocarbon deposits. If an area is predicted to be of particular interest for such an operation, provision may be made to place an operating facility and, where appropriate, to produce hydrocarbons from the results of the simulation process.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pure & Applied Mathematics (AREA)
  • Educational Technology (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Geology (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Business, Economics & Management (AREA)
  • Educational Administration (AREA)
  • Paleontology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geophysics (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)
  • Geophysics And Detection Of Objects (AREA)

Description

Procédé de simulation de formation géologique d' une zone fluviale .
L'invention se rapporte au domaine de la simulation de processus géologiques pour l'étude du sous-sol. On s'intéresse particulièrement aux phénomènes de formation des zones fluviales.
Les sables constituent des régions d'intérêt pour la prospection pétrolière. On cherche par conséquent à prédire leur présence et leurs caractéristiques.
On peut en particulier tenter de simuler le processus géologique de leur formation. Toutefois, un défi réside encore dans l'élaboration d'un modèle qui sache, tout en s 'appuyant sur des données d'observation, prendre en compte les spécificités des zones fluviales pour fournir un résultat exploitable pour la prospection.
Selon un aspect, il est prévu un procédé de simulation de la formation géologique d'une zone fluviale comprenant :
- on fournit des données d'observation de la zone fluviale,
on définit un modèle spatial de la zone fluviale comprenant au moins une région amont et une région aval, à partir des données d'observation,
- on simule le déplacement de particules dans le modèle spatial en superposant :
un terme déterministe défini à partir des données d'observation, et
un terme stochastique paramétré au moins par les données d'observation.
Grâce à ces dispositions, on peut prendre en compte à la fois l'écoulement fluide des particules dans la zone, et introduire une perturbation probabiliste.
Quand on parle de superposition de deux termes, on comprendra que le déplacement simulé se compose comme la somme du terme déterministe et du terme stochastique.
Dans des modes de réalisation préférés de l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions ci-après.
Le terme stochastique peut comprendre la superposition d'un terme méandriforme et d'une perturbation aléatoire. Ce faisant, la modélisation du chenal est rendue plus réaliste.
Le terme méandriforme peut comprendre une superposition d'au moins une fonction trigonométrique . Une telle représentation est réaliste pour un terme méandriforme, et facilement paramétrable.
Les données d'observation peuvent comprendre au moins l'un des types de données suivant :
- un gradient de vitesse d'écoulement dans la zone,
des paramètres géométriques de la zone,
des données issues de forage,
des données issues d'imagerie, notamment d'imagerie sismique.
Selon certains aspects, on peut représenter quelques phénomènes physiques propres aux chenaux, lorsqu'on simule le déplacement de particules:
on peut estimer la probabilité de l'occurrence d'une séparation d'une particule en deux particules indépendantes et, en cas d'occurrence, procéder à une telle séparation ;
lorsque deux particules indépendantes se trouvent à un même emplacement, on peut estimer la probabilité de l'occurrence d'une jonction des deux particules indépendantes en une et, en cas d'occurrence, on procède à une telle jonction ;
on peut estimer la probabilité de l'occurrence d'une disparition d'une particule et, en cas d'occurrence, procéder à une telle disparition. Par exemple, la probabilité d'occurrence de tels événements peut être estimée en fonction des données d ' observation .
En variante, lorsqu'on simule le déplacement de particules, on peut superposer également un champ de déplacement imposé par une source extérieure. On peut ainsi prendre en compte des sources extérieures telles que, par exemple, une dérive littorale.
Selon un mode de réalisation, on peut définir une pluralité de zones du modèle spatial entre la région amont et la région aval, et faire différer un paramétrage du terme stochastique (en particulier, la fonction trigonométrique ) dans une première zone d'un paramétrage du terme stochastique dans au moins une deuxième zone. On peut ainsi prendre en compte le fait qu'une zone fluviale peut présenter diverses régions très différentes.
Par exemple, lorsqu'on définit un modèle spatial, on peut définir une ceinture chenalisante définissant une largeur et une épaisseur dans laquelle les particules sont autorisées à se déplacer.
Selon un autre aspect, on prévoit un procédé de fabrication d'une installation d'extraction d'hydrocarbure comprenant la mise en œuvre d'un tel procédé de simulation.
Selon un autre aspect, on prévoit un produit programme d'ordinateur de simulation de formation géologique d'une zone fluviale, destiné à être stocké dans une mémoire d'une unité centrale, et/ou stocké sur un support mémoire destiné à coopérer avec un lecteur de ladite unité centrale et/ou téléchargé via un réseau de télécommunication, caractérisé en ce qu'il comprend des instructions pour mettre en œuvre un tel procédé lorsqu'il est exécuté sur une machine programmable.
D'autres caractéristiques et avantages de l'invention apparaîtront au cours de la description suivante d'une de ses formes de réalisation, donnée à titre d'exemple non limitatif, en regard des dessins joints.
Sur les dessins :
la figure 1 est une vue schématique d'un chenal, la figure 2a est une vue schématique de dessus d'un modèle de chenal,
la figure 2b est une vue schématique en coupe de chenal ,
la figure 3 est un schéma illustratif du zonage dans un chenal,
- la figure 4 est un schéma d'ensemble illustrant la modélisation de trajectoire de particules dans une partie de chenal,
la figure 5 est schéma illustratif similaire à celui de la figure 3 décrivant un procédé de séparation,
- la figure 6 est une vue similaire à la figure 3 décrivant un procédé de fusion,
- la figure 7 est vue similaire à la figure 3 décrivant un procédé de disparition,
la figure 8 est une vue schématique de dessus d'un modèle dans lequel on symbolise un champ de déplacement imposé par une source extérieure.
Sur les différentes figures, les mêmes références désignent des éléments identiques ou similaires.
La figure 1 décrit schématiquement une zone fluviale 1 s 'étendant entre une zone amont 2a et une zone aval 2b. La zone fluviale à modéliser est considérée sur une certaine largeur 1, une certaine épaisseur e et une certaine longueur L. Epaisseur et largeur sont par exemple variables le long de cette longueur.
Une zone fluviale 1 comporte généralement un chenal
3 s 'étendant entre la zone amont et la zone aval. Ce chenal a été formé par l'écoulement de particules au cours du temps de la zone amont à la zone aval.
Par observation, il est connu qu'un chenal peut présenter au moins une des configurations suivantes, celle- ci étant listée ici dans l'ordre où on les trouve généralement entre la zone amont et la zone aval :
- un système enchevêtré 4a, dans lequel de nombreux petits canaux individuels se croisent et se séparent en étant interposé par de grandes plages de terre,
- un système à faible sinuosité 4b, où le canal s'étend sensiblement rectiligne,
- un système à méandres 4c, dans lequel les canaux présentent une forte sinuosité,
- un système à anastomose 4d, présentant de nombreuses connections entre canaux,
- un système deltaïque 4e, dans lequel un canal principal se divise en un nombre de canaux secondaires débouchant dans la mer, et
- un système turbide 4f lié à l'interaction avec un autre écoulement, tel qu'un écoulement maritime.
Comme on peut le voir, en particulier sur la figue 1, ces différents systèmes présentent des caractéristiques géométriques très différentes, et partant, les propriétés géologiques du chenal dans ces différentes zones peuvent être très diverses.
On cherche en particulier à modéliser la formation géologique d'une zone fluviale telle que celle représentée sur la figure 1.
Pour cela, on définit un modèle spatial de la zone fluviale. Un tel modèle spatial peut être réalisé sous la forme d'un maillage surfacique ou volumique dans lequel chacun des éléments représente un emplacement de la zone fluviale .
On dispose également de données d'observations de la zone fluviale dont on veut modéliser la formation géologique. Ces données d'observation peuvent par exemple comprendre :
- la définition de l'emplacement d'une région amont et d'une région aval définies sur la base de la géométrie actuelle de la zone modélisée,
- un gradient de vitesse d'écoulement dans la zone fluviale, qui peut par exemple être défini sur la base d'un gradient hydraulique dans cette zone, et/ou à partir des données d'altitude des régions amont et aval,
- des paramètres géométriques de la zone fluviale, telle que, par exemple, la géométrie actuelle de la surface de la zone, des discontinuités observées, ...
- des données géologiques issues de forages réalisés dans la zone,
- des données géologiques issues d'imagerie, par exemple d'imagerie sismiques réalisées dans la zone,
La figure 2a représente ainsi de dessus un modèle spatial d'une région de la zone fluviale, présentant un ensemble d'éléments de maillage, par exemple parallélépipédiques (des carrés en vue de dessus), dans lequel un modèle de chenal 5 a déjà été généré par simulation. La figure 2b représente une vue en coupe selon la ligne B-B de la figure 2a représentant ainsi une tranche du modèle de chenal calculé 5.
Comme cela est représenté sur la figure 3, en partant des données d'observations, on peut représenter une trajectoire du chenal correspondant à la direction moyenne d'écoulement, s 'étendant de la zone amont 2a à la zone aval 2b. Dans l'exemple simplifié représenté sur la figure 3, ce chenal présente trois des zones distinctes mentionnées ci- dessus par relation avec la figure 1, qui sont identifiées sur la figure 3 par les références 6a, 6b et 6c. Chacune de ces zones peut faire l'objet d'une modélisation séparée, en respectant des contraintes de connexité à l'interface entre deux zones. La modélisation d'une de ces zones va être décrite ci-après par référence à la figure 4.
On supposera par exemple qu'il s'agit ici de la zone 6a. Comme cela est visible sur la fenêtre 7a de la figure 4, on définit une ceinture chenalisante 8 qui entoure la courbe moyenne d'écoulement 9 dans cette zone. La ceinture chenalisante 8 définit une largeur dans laquelle les particules sont autorisées à s'écouler, cette largeur étant éventuellement variable le long de la courbe moyenne d'écoulement (la largeur est ici représentée constante) . De plus, la ceinture chenalisante peut s'étendre non seulement sur la surface, mais également sur une certaine profondeur en dessous du plan de représentation de la fenêtre 7a. La zone représentée comporte une charge hydraulique amont 10a et une charge hydraulique aval 10b.
Le déplacement des particules dans le modèle spatial est simulé en superposant un terme déterministe et un terme stochastique. Par exemple, on écrit :
ΔΧ = AD + ΔΨ,
où :
ΔΧ représente le déplacement de la particule, AD représente le terme déterministe, et
ΔΨ représente le terme stochastique.
Le terme déterministe est défini à partir du champ d'écoulement représenté sur la fenêtre 7a de la figure 4.
Le terme stochastique est ajouté à ce terme pour définir une perturbation du terme déterministe. En particulier, le terme stochastique peut être lui-même réalisé comme la superposition de deux termes. Le premier de ces termes est un terme méandriforme qui peut facilement être réalisé, du point de vue mathématique, par un ou plusieurs termes trigonométriques (sinusoïdaux) comme représenté sur la fenêtre 7b de la figure 4. Sur cette fenêtre, on a représenté par la courbe 11 une période et demi d'un sinus s 'étendant entre les points a et g, et d'une amplitude donnée. Les points a et g correspondent aux zones amont et aval de la zone modélisée, et le nombre de périodes ainsi que l'amplitude du sinus peuvent être paramétrés indépendamment pour cette zone.
On superpose au terme méandriforme une perturbation aléatoire représentée sur la fenêtre 7b par les flèches 12. Ce terme aléatoire est défini comme une perturbation de la courbe sinusoïdale 11. L'orientation et/ou l'amplitude de la perturbation peuvent être déterminées de manière aléatoire sur une pluralité de points 13 de la courbe sinusoïdale 11, par exemple équidistant entre eux, ou d'abscisse équidistante, comme représenté.
La fenêtre 7c représente une trajectoire simulée 13 de particules, au sein de la ceinture chenalisante 8, obtenue à partir d'une superposition du type mentionné ci- dessus. On peut en particulier prévoir que le paramétrage (période, amplitude) de la fonction sinusoïdale soit déterminé en fonction de la zone modélisée pour la zone fluviale. Le paramétrage (fréquence, intensité) de la perturbation aléatoire peut être déterminé en fonction de la zone considérée. Ce paramétrage peut également être effectué en prenant en compte d'autres données d'observation, parmi les données d'observation mentionnées ci-dessus .
Selon un aspect particulier, on peut en outre modéliser certains phénomènes survenant aux particules au cours de leur écoulement.
Un de ces phénomènes représenté sur la figure 5, est la séparation ou mitose. Au cours de chaque calcul de déplacement, on peut déterminer une probabilité que la particule, en l'emplacement donné, se divise.
De telles divisions sont visibles en 15 sur la figure 5. La probabilité de voir subvenir une division peut dépendre de la zone considérée. Sur l'exemple de la figure 5, on peut par exemple prévoir qu'aucune division ne soit possible dans la zone 6a. Dans la zone 6b, on peut prévoir que des divisions soient possibles, et que, en cas de division, une particule se divise en un nombre maximal de particules filles (ici le nombre maximal de particules filles est de deux) . Par ailleurs, on peut également paramétrer le nombre maximal de divisions autorisées pour une particule dans une zone donnée (au maximum deux divisions dans la zone 6b représentées sur la figure 5 dans l'exemple représenté).
Un autre phénomène pouvant être modélisé est représenté sur la figure 6. Il s'agit, si deux particules sont situées à un emplacement à un même instant, de procéder à une jonction ou fusion de ces particules. Une telle fusion est visible, sur la figure 6, telle que représentée par les références 16. On notera que le processus de fusion pourra être mis en œuvre en parallèle du processus de division représenté ci-dessus, des divisions survenant, sur l'exemple de la figure 6, au point 15.
La possibilité de mettre en œuvre un phénomène de fusion est probabilisée . Cette probabilité peut par exemple dépendre de la zone considérée. Par exemple, on peut prévoir que l'on procède systématiquement à une fusion dans la zone 6a. Dans la zone 6b, telle que représentée, on peut prévoir ne jamais mettre en œuvre de fusion.
Un autre phénomène pouvant être modélisé est représenté sur la figure 7. Il s'agit, en l'occurrence, d'un phénomène de disparition d'une particule. On peut par exemple prévoir qu'une particule va disparaître après avoir parcouru une certaine distance, ou probabiliser la possibilité de voir survenir une disparition en fonction du trajet parcouru par la particule. De telles disparitions sont visibles, référencées par la référence 17 sur la figure 7. Séparation, jonction et disparition peuvent être probabilisées en fonction de données d'observation.
Un autre phénomène pouvant être modélisé est la superposition d'une autre source d'énergie à l'écoulement dans la zone fluviale. A titre d'exemple, on peut par exemple prévoir que cette autre source d'énergie soit une dérive littorale telle que représentée par la flèche 18 sur la figure 8. Ce terme peut par exemple venir s'ajouter au terme déterministe directement déterminé par le gradient hydraulique dans la zone lors du calcul du déplacement des particules .
Lorsqu'on procède à la simulation, on simule le déplacement de particules dans la zone fluviale selon l'une et/ou l'autre des équations et phénomènes présentés ci- dessous. Cette simulation est par exemple mise en œuvre sur une machine programmable, telle qu'un ordinateur, ayant stocké en mémoire un programme informatique stocké en mémoire localement ou accessible via un réseau, et programmé pour mettre en œuvre les étapes ci-dessus. On peut mettre en œuvre un critère de détermination de fin de simulation, par exemple basé sur le nombre de particules simulées, sur le temps de simulation, et/ou sur la ressemblance de certains résultats de la simulation avec des données d'observation.
La simulation peut permettre d'évaluer l'un et/ou l'autre des paramètres suivants :
la trajectoire du chenal,
l'azimut des points qui composent sa trajectoire (en particulier l'orientation par rapport au nord) ,
- le rayon de courbure des trajectoires,
- leur terminaison,
pour chaque point de la trajectoire, la distance par rapport à la base (fond) et au bord du chenal,
- pour chaque maille du chenal, une propriété indiquant une distance par rapport au bord le plus proche, au centre, à la base ou au toit (surface du chenal) ,
- la zonation, La modélisation réaliste de la formation du chenal, obtenue ci-dessus, et les paramètres ci-dessus peuvent être utiles dans le cadre de l'exploitation pétrolifère de certains emplacements du chenal. On peut par exemple prédire, sur un chenal faisant parfois plusieurs centaines de kilomètres de long, les zones susceptibles de contenir des gisements d'hydrocarbures. Si une zone est prédite comme étant particulièrement intéressante pour une telle exploitation, on peut prévoir d'y placer une installation d'exploitation et, le cas échéant, de produire des hydrocarbures à partir des résultats du procédé de simulation.

Claims

REVENDICATIONS
1. Procédé de simulation de la formation géologique d'une zone fluviale comprenant :
- on fournit des données d'observation de la zone fluviale ,
on définit un modèle spatial de la zone fluviale comprenant au moins une région amont (2a) et une région aval (2b), à partir des données d'observation,
- on simule le déplacement de particules dans le modèle spatial en superposant :
un terme déterministe défini à partir des données d'observation, et
un terme stochastique paramétré au moins par les données d'observation.
2. Procédé selon la revendication 1, dans lequel le terme stochastique comprend la superposition d'un terme méandriforme et d'une perturbation aléatoire.
3. Procédé selon la revendication 2 dans lequel le terme méandriforme comprend une superposition d'au moins une fonction trigonométrique .
4. Procédé selon l'une des revendications 1 à 3 dans lequel les données d'observation comprennent au moins l'un des types de données suivant :
- un gradient de vitesse d'écoulement dans la zone,
des paramètres géométriques de la zone,
des données issues de forage,
des données issues d'imagerie, notamment d'imagerie sismique.
5. Procédé selon l'une quelconque des revendications 1 à 4, dans lequel, lorsqu'on simule le déplacement de particules :
on estime la probabilité de l'occurrence d'une séparation (15) d'une particule en deux particules indépendantes et, en cas d'occurrence, on procède à une séparation .
6. Procédé selon l'une quelconque des revendications 1 à 5, dans lequel, lorsqu'on simule le déplacement de particules, lorsque deux particules indépendantes se trouvent à un même emplacement :
on estime la probabilité de l'occurrence d'une jonction (16) des deux particules indépendantes en une et, en cas d'occurrence, on procède à une jonction.
7. Procédé selon l'une quelconque des revendications 1 à 6, dans lequel, lorsqu'on simule le déplacement de particules :
on estime la probabilité de l'occurrence d'une disparition (17) d'une particule et, en cas d'occurrence, on procède à une disparition.
8. Procédé selon l'une des revendications 5 à 7, dans lequel ladite probabilité est estimée en fonction des données d'observation.
9. Procédé selon l'une quelconque des revendications 1 à 8, dans lequel, lorsqu'on simule le déplacement de particules, on superpose également un champ de déplacement (18) imposé par une source extérieure.
10. Procédé selon l'une quelconque des revendications 1 à 9, dans lequel on définit une pluralité de zones (6a, 6b, 6c) du modèle spatial entre la région amont et la région aval, et dans lequel un paramétrage du terme stochastique dans une première zone diffère d'un paramétrage du terme stochastique dans au moins une deuxième zone.
11. Procédé selon la revendication 3 et la revendication 10, dans lequel la fonction trigonométrique est choisie en fonction de la zone.
12. Procédé selon l'une quelconque des revendications 1 à 11, dans lequel, lorsqu'on définit un modèle spatial, on définit une ceinture chenalisante (8) définissant une largeur et une épaisseur dans laquelle les particules sont autorisées à se déplacer.
13. Procédé de fabrication d'une installation d'extraction d'hydrocarbure comprenant la mise en œuvre d'un procédé de simulation selon l'une quelconque des revendications 1 à 12.
14. Produit programme d'ordinateur de simulation de formation géologique d'une zone fluviale, destiné à être stocké dans une mémoire d'une unité centrale, et/ou stocké sur un support mémoire destiné à coopérer avec un lecteur de ladite unité centrale et/ou téléchargé via un réseau de télécommunication, caractérisé en ce qu'il comprend des instructions pour mettre en œuvre le procédé selon l'une quelconque des revendications 1 à 12 lorsqu'il est exécuté sur une machine programmable.
PCT/FR2011/052769 2010-12-03 2011-11-25 Procede de simulation de formation geologique d'une zone fluviale WO2012072923A2 (fr)

Priority Applications (4)

Application Number Priority Date Filing Date Title
RU2013130263/12A RU2013130263A (ru) 2010-12-03 2011-11-25 Способ моделирования формирования флювиальной зоны
BR112013013491A BR112013013491A2 (pt) 2010-12-03 2011-11-25 procedimento de simulação da formação geológica de uma zona fluvial, procedimento de fabricação de uma instalação de extração de hidrocarboneto e produto de programa de computador
US13/991,360 US10007741B2 (en) 2010-12-03 2011-11-25 Method for simulating the geological formation of a river zone
NO20130806A NO345728B1 (no) 2010-12-03 2011-11-25 Fremgangsmåte for simulering av geologisk formasjon i en elvesone

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1060053 2010-12-03
FR1060053A FR2968443A1 (fr) 2010-12-03 2010-12-03 Procede de simulation de formation geologique d'une zone fluviale

Publications (1)

Publication Number Publication Date
WO2012072923A2 true WO2012072923A2 (fr) 2012-06-07

Family

ID=44201540

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/FR2011/052769 WO2012072923A2 (fr) 2010-12-03 2011-11-25 Procede de simulation de formation geologique d'une zone fluviale

Country Status (7)

Country Link
US (1) US10007741B2 (fr)
AR (1) AR084068A1 (fr)
BR (1) BR112013013491A2 (fr)
FR (1) FR2968443A1 (fr)
NO (1) NO345728B1 (fr)
RU (1) RU2013130263A (fr)
WO (1) WO2012072923A2 (fr)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2994315A1 (fr) * 2012-08-06 2014-02-07 Total Sa Procede de determination de trajectoires de chenaux
WO2014053423A3 (fr) * 2012-10-05 2014-07-24 Total Sa Procédé pour déterminer une région karstique
WO2019243858A1 (fr) 2018-06-20 2019-12-26 Total Sa Procédé de détermination d'une composition de sous-sol réel
WO2019243857A1 (fr) 2018-06-20 2019-12-26 Total Sa Procédé de détermination de composition de sous-sol
WO2019243865A1 (fr) 2018-06-20 2019-12-26 Total Sa Procédé de détermination d'une composition de sous-sol réel
EP3819682A1 (fr) 2019-11-08 2021-05-12 Total Se Procédé de détermination de formation géologique de sous-sol réel
EP3819678A1 (fr) 2019-11-08 2021-05-12 Total Se Procédé de détermination de formation géologique de sous-sol réelle
EP3819679A1 (fr) 2019-11-08 2021-05-12 Total Se Procédé de détermination de formation géologique de sous-sol réelle
EP3819680A1 (fr) 2019-11-08 2021-05-12 Total Se Procédé de détermination de formation géologique de sous-sol réelle
EP3819681A1 (fr) 2019-11-08 2021-05-12 Total Se Procédé de détermination de formation géologique de sous-sol réelle

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109098143A (zh) * 2018-10-08 2018-12-28 浙江省水利河口研究院 一种增强涌潮的人工河道设计方法
US20220236448A1 (en) * 2019-05-10 2022-07-28 Totalenergies Se Method for modelling a water current induced by a river in a geological gridded model
WO2020229863A1 (fr) * 2019-05-10 2020-11-19 Total Se Procédé de modélisation de la formation d'une zone sédimentaire par simulation de transport de particules induit par du courant
US11232242B2 (en) * 2019-07-22 2022-01-25 Red Hat, Inc. Sensory data generator
CN110672803B (zh) * 2019-09-29 2021-08-31 浙江水利水电学院 一种实验水槽隔断结构及其安装方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2744224B1 (fr) * 1996-01-26 1998-04-17 Inst Francais Du Petrole Methode pour simuler le remplissage d'un bassin sedimentaire
US6480790B1 (en) * 1999-10-29 2002-11-12 Exxonmobil Upstream Research Company Process for constructing three-dimensional geologic models having adjustable geologic interfaces
FR2849211B1 (fr) * 2002-12-20 2005-03-11 Inst Francais Du Petrole Methode de modelisation pour constituer un modele simulant le remplissage multilithologique d'un bassin sedimentaire
EP1618414A2 (fr) * 2003-03-13 2006-01-25 ExxonMobil Upstream Research Company Procede de prediction de la granulometrie a partir de l'epaisseur d'un gisement
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
US8190414B2 (en) * 2008-03-26 2012-05-29 Exxonmobil Upstream Research Company Modeling of hydrocarbon reservoirs containing subsurface features
WO2012078217A1 (fr) * 2010-12-08 2012-06-14 Exxonmobil Upstream Research Company Construction de modèles géologiques à partir de concepts géologiques

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2994315A1 (fr) * 2012-08-06 2014-02-07 Total Sa Procede de determination de trajectoires de chenaux
WO2014023908A1 (fr) * 2012-08-06 2014-02-13 Total Sa Procede de determination de trajectoires de chenaux
US10180516B2 (en) 2012-08-06 2019-01-15 Total Sa Method of determining channelway trajectories
WO2014053423A3 (fr) * 2012-10-05 2014-07-24 Total Sa Procédé pour déterminer une région karstique
US10101497B2 (en) 2012-10-05 2018-10-16 Total Sa Method for determining a karstic region
WO2019243858A1 (fr) 2018-06-20 2019-12-26 Total Sa Procédé de détermination d'une composition de sous-sol réel
WO2019243857A1 (fr) 2018-06-20 2019-12-26 Total Sa Procédé de détermination de composition de sous-sol
WO2019243865A1 (fr) 2018-06-20 2019-12-26 Total Sa Procédé de détermination d'une composition de sous-sol réel
EP3819682A1 (fr) 2019-11-08 2021-05-12 Total Se Procédé de détermination de formation géologique de sous-sol réel
EP3819678A1 (fr) 2019-11-08 2021-05-12 Total Se Procédé de détermination de formation géologique de sous-sol réelle
EP3819679A1 (fr) 2019-11-08 2021-05-12 Total Se Procédé de détermination de formation géologique de sous-sol réelle
EP3819680A1 (fr) 2019-11-08 2021-05-12 Total Se Procédé de détermination de formation géologique de sous-sol réelle
EP3819681A1 (fr) 2019-11-08 2021-05-12 Total Se Procédé de détermination de formation géologique de sous-sol réelle
US11487034B2 (en) 2019-11-08 2022-11-01 Total Se Method for determination of real subsoil geological formation
US11614556B2 (en) 2019-11-08 2023-03-28 Total Se Method for determination of real subsoil geological formation
US11630235B2 (en) 2019-11-08 2023-04-18 Total Se Method for determination of real subsoil geological formation
US11656380B2 (en) 2019-11-08 2023-05-23 Total Se Method for determination of real subsoil geological formation
US11995382B2 (en) 2019-11-08 2024-05-28 Totalenergies Onetech Method for determination of real subsoil geological formation

Also Published As

Publication number Publication date
NO20130806A1 (no) 2013-06-10
RU2013130263A (ru) 2015-01-10
US10007741B2 (en) 2018-06-26
BR112013013491A2 (pt) 2016-09-06
US20130262063A1 (en) 2013-10-03
FR2968443A1 (fr) 2012-06-08
NO345728B1 (no) 2021-07-05
AR084068A1 (es) 2013-04-17

Similar Documents

Publication Publication Date Title
WO2012072923A2 (fr) Procede de simulation de formation geologique d'une zone fluviale
EP3144468B1 (fr) Procédé pour caracteriser le réseau de fractures d'un gisement fracture et procédé pour l'exploiter
Allen et al. Orogenic plateau growth: Expansion of the Turkish‐Iranian Plateau across the Zagros fold‐and‐thrust belt
Ielpi et al. Planform architecture, stratigraphic signature and morphodynamics of an exhumed Jurassic meander plain (Scalby Formation, Yorkshire, UK)
Quigley et al. Recurrent liquefaction in Christchurch, New Zealand, during the Canterbury earthquake sequence
Beaubouef Deep-water leveed-channel complexes of the Cerro Toro Formation, Upper Cretaceous, southern Chile
Liang et al. Quantifying the patterns and dynamics of river deltas under conditions of steady forcing and relative sea level rise
MacInnes et al. Tsunami geomorphology: Erosion and deposition from the 15 November 2006 Kuril Island tsunami
EP2530493B1 (fr) Méthode pour construire un maillage d'un réseau de fractures à partir de diagramme de voronoi
Van der Merwe et al. Depositional architecture of sand-attached and sand-detached channel-lobe transition zones on an exhumed stepped slope mapped over a 2500 km2 area
EP2356493B1 (fr) Procede de modelisation geologique de donnees sismiques par correlation de traces
Darke et al. Beach–dune sediment budgets and dune morphodynamics following coastal dune restoration, Wickaninnish Dunes, Canada
EP2622385B1 (fr) Simulation de karstification
Kereszturi et al. LiDAR-based quantification of lava flow susceptibility in the City of Auckland (New Zealand)
Kostic et al. Role of upper-flow-regime bedforms emplaced by sediment gravity flows in the evolution of deltas
Reston et al. Extreme crustal thinning in the south Porcupine Basin and the nature of the Porcupine Median High: implications for the formation of non-volcanic rifted margins
EP2650471B1 (fr) Procédé de sélection des positions de puits a forer pour l'exploitation d'un gisement pétrolier
Reston et al. Lithospheric extension from rifting to continental breakup at magma-poor margins: rheology, serpentinisation and symmetry
Selles et al. The eastern flank of the Merapi volcano (Central Java, Indonesia): Architecture and implications of volcaniclastic deposits
Authemayou et al. Unraveling the roles of asymmetric uplift, normal faulting and groundwater flow to drainage rearrangement in an emerging karstic landscape
Åberg et al. GIS-based 3D sedimentary model for visualizing complex glacial deposition in Kersilö, Finnish Lapland
Xu et al. Detecting the spatial differentiation in settlement change rates during rapid urbanization in the Nanjing metropolitan region, China
CA2374896C (fr) Methode pour deformer graduellement une repartition initiale d'objets dans un milieu heterogene pour l'adapter a des contraintes physiques imposees
Bourillot et al. Structure and evolution of a Messinian mixed carbonate‐siliciclastic platform: the role of evaporites (Sorbas Basin, South‐east Spain)
Almeida Junior et al. River‐dominated and tide‐influenced shelf‐edge delta systems: Coarse‐grained deltas straddling the Early–Middle Jurassic shelf–slope break and transforming downslope, Lajas–Los Molles formations, Neuquén Basin, Argentina

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 11802517

Country of ref document: EP

Kind code of ref document: A2

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 13991360

Country of ref document: US

ENP Entry into the national phase

Ref document number: 2013130263

Country of ref document: RU

Kind code of ref document: A

122 Ep: pct application non-entry in european phase

Ref document number: 11802517

Country of ref document: EP

Kind code of ref document: A1

REG Reference to national code

Ref country code: BR

Ref legal event code: B01A

Ref document number: 112013013491

Country of ref document: BR

ENP Entry into the national phase

Ref document number: 112013013491

Country of ref document: BR

Kind code of ref document: A2

Effective date: 20130531