WO2014023916A1 - Procede d'assistance a la modelisation geologique par regroupements de mailles - Google Patents

Procede d'assistance a la modelisation geologique par regroupements de mailles Download PDF

Info

Publication number
WO2014023916A1
WO2014023916A1 PCT/FR2013/051900 FR2013051900W WO2014023916A1 WO 2014023916 A1 WO2014023916 A1 WO 2014023916A1 FR 2013051900 W FR2013051900 W FR 2013051900W WO 2014023916 A1 WO2014023916 A1 WO 2014023916A1
Authority
WO
WIPO (PCT)
Prior art keywords
column
meshes
determined
value
decompositions
Prior art date
Application number
PCT/FR2013/051900
Other languages
English (en)
Inventor
Gérard MASSONNAT
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 GB1501860.9A priority Critical patent/GB2519037B/en
Priority to NO20150173A priority patent/NO346799B1/en
Priority to US14/420,153 priority patent/US10274641B2/en
Publication of WO2014023916A1 publication Critical patent/WO2014023916A1/fr

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V20/00Geomodelling in general
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05Geographic models

Definitions

  • the present invention relates to the field of computerized simulation of data in a mesh model, particularly the field of computerized simulation of data representing geological information.
  • the number of meshes of such models can be close to several million when creating fine models. With current technologies, it is not always possible to use these fine models in dynamic simulations, as flow simulators can hardly handle models exceeding 500,000 meshes.
  • geologists and reservoir engineers simplify these models to reduce the number of meshes while preserving at best the information they possess: rough models are then obtained.
  • the present invention improves the situation.
  • the present invention proposes to combine the constraints of geologists and reservoir engineers and to propose several optimized solutions of mesh groupings in order to reduce the final number of meshes of the model.
  • the present invention thus aims at a method of grouping meshes in a one to three-dimensional meshed geological model, the geological model comprising a plurality of meshes, and each cell of said model being associated with at least one numerical parameter representing a geological property.
  • the method comprises the steps: for at least one column of at least n stitches of the model:
  • a2 / determining a dispersion value of said decomposition as a function, at least, of the new numerical parameters associated with said subsets determined in step a1 /;
  • n is a strictly positive integer.
  • geological property refers to any physical property of a soil related to its composition, structure, evolution, etc. By way of illustration, this property may correspond to the type of majority facies of the soil corresponding to said mesh. In addition, this property may correspond to the porosity of the soil, whether or not there is fracture, the direction of flow, etc. soil corresponding to said mesh.
  • a column can be constituted in the z direction of a set of cubic meshes "stacked" on top of each other, vertically. It may be advantageous for the position of this column in the model to correspond to the position of a wellbore.
  • a "subset of adjacent meshes" is a set of meshes each having at least one face in common with one of the other meshes.
  • the decomposition determined does not comprise a subset of disjoint meshes and may allow the subsequent creation of a new model having coarse meshes corresponding to these subsets.
  • the "new numerical parameter of a subset" can be determined, for example, by averaging, weighted averaging, thresholding, etc. different numeric parameters associated with the meshes of this subset.
  • the "dispersion value of a decomposition" can be determined, for example, by calculating the variance of the new numerical parameters of the subsets of this decomposition.
  • the geologist and the reservoir engineer try to optimize a certain number of sometimes antagonistic parameters and the provision of these decompositions as well as their dispersion values can allow them to perform such arbitration.
  • This arbitration can be, by way of illustration, an arbitration between the simulation times related to the number of meshes of the model and the accuracy of the simulation related to the proximity of the model to the geological reality (which can depend on the dispersion value) .
  • the method may further comprise:
  • step b ordering said decompositions provided in step b / according to the dispersion value determined for each of the decompositions;
  • Classification of the decompositions according to the determined dispersion value may allow the final operator to simplify his arbitration by allowing him to select the target decomposition more quickly according to his objectives (for example, the speed of simulation and / or the accuracy) .
  • each decomposition determined for one of the plurality of columns can have a homologous decomposition in each of the other columns.
  • the method may further include:
  • a model may include several columns representing these wells.
  • decomposition homologous to a decomposition of a column a representation decomposition identical to a decomposition in another column.
  • a decomposition in a first column can group the first three meshes of this column into a first subset, then the next ten meshes into a second subset, and so on.
  • the decomposition homologous to this decomposition then comprises the same groupings, but on a second distinct column: the first three meshes of this second column in a first subset, then the next ten meshes of the second column in a second subset, etc.
  • This correction can make it possible to better account for the geological reality of the soil and thus allow better arbitrations to be made between different coarse meshes. Indeed, it may be possible to choose a decomposition (and therefore a new mesh of a coarse model) among the decompositions having the best dispersion value and this for a majority of columns among the columns of the model.
  • This correction can consist in assigning the same dispersion value to all the columns for the homologous decompositions.
  • This same value may, for example, be the average of the dispersion values of the homologous decompositions of the different columns or be one of the dispersion values of the columns.
  • each column having an importance parameter the dispersion value of each decomposition can be a function of said importance parameter of the column.
  • An "importance parameter” is a parameter that can be used to weight the importance of certain wells and thus the columns of meshes associated with these wellbores.
  • this significance parameter may be a function of well flow, well capacity prediction, well size (main well, secondary wells, etc.).
  • all possible decompositions of the column can be determined in step a /.
  • each cell of the column may be associated with a value (for example numeric) completion value representing a degree of completion of the mesh.
  • the plurality of decompositions can be determined in step a / so that the completion value of each cell of each of the m subsets can be different from the completion values of the other cells of the same subset of the subset. plus a determined threshold.
  • the decomposition has two groupings (the first subset corresponding to the first ten meshes of the column and the second subset corresponding to the other meshes of the column), the difference in maximum completion values between two meshes of the column. first subset does not exceed this threshold.
  • the determined threshold may be less than half the difference between the highest completion value and the smallest completion value associated with the cells of the column.
  • the threshold is less than 5.
  • the decomposition can correspond to the grouping of the first three meshes (because the difference of the values 1, 3, 2 is at most 2, ie a difference of less than 5), following two (because the difference of the values 9, 8 is at the maximum of 1, ie a difference of less than 5) and of the following three (because the difference of the values 1, 5, 3 is at most 4, ie a difference less than 5).
  • the decomposition can not correspond to the grouping of the first 4 meshes (because the difference of the values 1, 3, 2, 9 is maximum of 8, ie a difference greater than 5), and of the following four meshes (because the difference between the values 8, 1, 5, 3 is at most 7, ie a difference greater than 5).
  • the geological model may be a two or three-dimensional mesh model, and each cell of each column may be associated with a completion value representing a degree of completion of the mesh.
  • each cell of the first column can be associated with a completion value associated with a mesh of the second column, called virtual completion value.
  • the plurality of decompositions for the first column can be determined in step a / so that the virtual completion value of each cell of each of the m subsets can be different from the virtual completion values of the other cells of the same. subset of not more than one determined threshold.
  • the use of constraints outside the column can make it possible to ensure that the decompositions determined for one column of the model are compatible with the decompositions determined for the other columns of the model.
  • the columns have at least one homologous decomposition with the other columns.
  • one of the thresholds determined above may be equal to 0.
  • the mesh determination values can be identical.
  • the determined threshold may be less than half the difference between the largest virtual completion value and the smallest virtual completion value associated with the meshes of the first column.
  • the dispersion value of a decomposition can be determined as a function of the interclass variance of the determined subsets.
  • the variance of the set constituted by the new numerical parameters associated with the subsets of a decomposition calculated in step a1 / is called "interclass variance of subsets".
  • each cell of the column can be associated with a value (for example, digital) of heterogeneity representing a geological parameter of the soil
  • the plurality of decompositions can be determined in step a / so that the value of heterogeneity each cell of each of the m subsets may be different from the heterogeneity values of the other meshes of the same subset by at most a determined threshold.
  • heterogeneity value is a numerical value representing a geological parameter whose heterogeneity is to be preserved during decompositions.
  • a device for grouping meshes in a mesh geological model can be advantageous in itself.
  • the present invention also provides a device for grouping meshes in a one to three-dimensional meshed geological model, the geological model comprising a plurality of meshes, each mesh of said model being associated with at least one numeric parameter representing a geological property, the device comprising:
  • n is a strictly positive integer.
  • a computer program, implementing all or part of the method described above, installed on a pre-existing equipment, is in itself advantageous, since it allows the grouping of meshes in a meshed geological model.
  • the present invention also relates to a computer program comprising instructions for implementing the method described above, when this program is executed by a processor.
  • FIGS. 1a and 1b illustrate a remeshing of a three-dimensional geological model
  • FIG. 3 illustrates a determination of a new subset numerical parameter and of dispersion value determination in one embodiment of the invention
  • FIG. 4 illustrates a flowchart of a possible method in one embodiment of the invention
  • FIG. 5 illustrates groupings constrained by mesh completion values in one embodiment of the invention
  • FIG. 6 illustrates a device that can implement the method of the invention in a particular embodiment.
  • Figures 1a and 1b illustrate a remeshing of a three-dimensional geological model.
  • the geological model of Figure 1a (called fine model) is represented by a 100 mesh cube of 216 meshes (i.e. 6x6x6 meshes). Each mesh, so-called fine mesh (i.e. 101, 102, 103, 104, etc.) has coordinates in the space formed by the orthonormal coordinate system (0, x, y, z).
  • this large number of meshes can be a handicap during mathematical calculations pushed on the geological model, like simulations of flows of fluids or geological deformations. It is thus advantageous to group some meshes, for example, in groups of 8 (e.g. 2x2x2) or in groups of 4 (e.g., 4x1 x1).
  • Figure 1b illustrates such a grouping.
  • the geological model (called coarse model) of Figure 1b is represented by a cube 100b mesh of 27 meshes (ie 3x3x3 meshes).
  • Each mesh, called coarse mesh (ie 106, etc.) has coordinates in the space formed by the orthonormal coordinate system (0, x, y, z) and can be matched with 8 fine meshes of the fine model.
  • the grouping of the 8 fine meshes corresponding to the cube 106 of Figure 1b is highlighted by the cube 105 of Figure 1a.
  • FIG. 2 illustrates a two-dimensional geological model.
  • This geological model has 32 meshes (i.e. 4x8 meshes) having coordinates in a formed space, the orthonormal coordinate system ( ⁇ , ⁇ , ⁇ ).
  • These meshes (201, 202, etc.) each include a geological parameter (Goo, G01, G-10, Gn, G20, etc.).
  • This geological parameter can notably correspond to:
  • This geological parameter (Goo, G01, G10, Gn, G20, etc.) is a numeric parameter, for example decimal.
  • FIG. 3 illustrates a determination of a new subset geological parameter in a one-dimensional model and according to an embodiment of the invention.
  • the set 351 of the five meshes constitutes an initial column and corresponds to the one-dimensional model before any decomposition into subsets.
  • Each mesh of this initial column is associated with a geological parameter (Go is associated with the mesh 300, G 4 is associated with the mesh 304).
  • the set 352 corresponds to a possible decomposition of the column 351 into four subsets (305 to 308).
  • the subassemblies 306, 307 and 308 respectively include the meshes 302, 303 and 304, while the subassembly 305 includes the two meshes 300 and 301. For each of these subsets (305 to 308), it is possible to determine a new geologic numerical parameter associated with this subset. For example :
  • the new geological parameter associated with the subset 305 is a function of the parameters Go and d (denoted by f (Go, Gi));
  • the new geological parameter associated with the subset 306 is a function of the parameter G2 (denoted f (G2));
  • the new geological parameter associated with the subset 307 is a function of the parameter Gs (denoted f (Gs));
  • the new geological parameter associated with the subset 308 is a function of the parameter G 4 (denoted f (G 4 )).
  • the set 353 corresponds to another decomposition of the column 351, in three subsets (309 to 31 1).
  • the subassemblies 309 and 31 1 respectively comprise the meshes 300 and 304, whereas the subassembly 310 includes the three meshes 301, 302 and 303.
  • For each of these subassemblies (309 to 31 1) it is possible to to determine a new geologic numerical parameter associated with this subset. For example :
  • the new geological parameter associated with the subset 309 is a function of the parameter Go (denoted f (GB));
  • the new geological parameter associated with the subset 310 is a function of the parameters d, G2 and Gs (denoted f (d, G2, G3));
  • the new geological parameter associated with the subset 31 1 is a function of the parameter G 4 (denoted f (G 4 )).
  • the function f allowing the determination of the new parameter can be, for example,
  • the function f can also be an expressed weighted mean function of the
  • geological 30i mesh as the permeability. It is also possible to note G ; this new parameter associated with the index subset j (conferring, if necessary, on each subset of a decomposition, an index). Moreover, a dispersion value can be determined for each of the decompositions 352 and 353. The dispersion value makes it possible in particular to quantify the loss of heterogeneity of the column during a decomposition.
  • this scatter value may correspond to the value of the interclass variance within the column. So, such a variance can be calculated from the
  • the variance ⁇ 353 is the maximum possible variance to be obtained during a decomposition with these assumptions.
  • Figure 4 illustrates a flowchart of a possible method in one embodiment of the invention.
  • the pattern received includes, for example, meshes, each mesh being cubic and each side of this cube being common with another cube.
  • this model may also indicate a set of mesh columns stacked on top of each other, adjacent and representing geological boreholes.
  • the object of such decompositions is to reduce the number of meshes of the geological model according to the vertical axis (ie for successive layers of the model). If all the decompositions are determined then their number is with i the initial number of meshes in the column, and f the number
  • subsets in the decomposition i.e. the number of meshes in the target model.
  • a decomposition for which no dispersion value has been determined is selected. So for each of the subsets of this decomposition, a a new geological parameter associated with the subset is calculated (step 403), for example, by averaging the geological parameters of the meshes contained in this subset.
  • a dispersion value associated with the selected decomposition is then calculated (step 404), for example by calculating the variance of the averages previously calculated for the subsets of this decomposition.
  • this ordered list may be provided (message 406) to a user or other processing system for display or recalculation.
  • the user for example, a geologist or a reservoir engineer
  • the decomposition that corresponds most closely to the geological reality found among the determined decompositions having the highest dispersion values.
  • Figure 5 illustrates groupings constrained by mesh completion values in one embodiment of the invention.
  • the total number of decompositions of a column can be relatively high and it can be tedious to analyze them all.
  • Completion constraints means any constraint related to the operation of the wellbore, as For example, the opening in the casing (or “casing” in English) of openings allowing the exploitation of layers of oil, gas, etc.
  • grouping R3 is not possible because the cells 504 and 505 do not have an identical completion value
  • each cell of the wells P1 and P2 is associated with a numerical value of completion representing a degree of completion of the cell.
  • this completion value is not necessarily 1 or 0. It can also be equal to 0.9 or 0.2 or any other numerical value.
  • the well columns P1 and P2 have meshes (respectively 500 to 507 and 510 to 517) each associated with a completion value. It is possible to match the meshes of these two wells according to their "depth" (i.e. their position in their respective columns).
  • the mesh 500 can be matched with the mesh 510
  • the mesh 501 can be matched with the mesh 51 1 ...
  • the mesh 507 can be matched with the mesh 517.
  • the process described above can then be performed by transposing the completion values with the virtual completion values.
  • the threshold ⁇ can be fixed in advance in a configuration file or dynamically according to the results of a first decomposition process.
  • An intermediate solution may consist in setting the ⁇ threshold to a low initial numerical value (for example, 0) and then an incremental upward change in its value, if the number of determined decompositions is not high enough.
  • FIG. 6 represents an example of device 600 of groupings of meshes in a meshed geological model in one embodiment of the invention.
  • the device 600 comprises a computer, comprising a memory 604 for storing instructions for implementing the method, received model data, and temporary data for performing the various steps of the method as described above. .
  • the computer further comprises a mesh grouping circuit 601.
  • This circuit can be, for example:
  • processor capable of interpreting instructions in the form of a computer program
  • a programmable electronic chip such as an FPGA (for "Field Programmable Gate Array”).
  • This computer has an input interface 602 for receiving mesh geological model data, and an output interface 603 for providing decomposition data. Finally, the computer comprises, to allow easy interaction with a user, a screen 605 and a keyboard 606.
  • FIG. 4 is a typical example of a program, some of whose instructions can be carried out with the device described. As such, FIG. 4 may correspond to the flowchart of the general algorithm of a computer program within the meaning of the invention.
  • the present invention is not limited to the embodiments described above as examples; it extends to other variants.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Geometry (AREA)
  • Databases & Information Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Remote Sensing (AREA)
  • Computer Graphics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

La présente invention concerne un procédé d'assistance à la modélisation d'un sous-sol géologique par regroupements de mailles dans un modèle géologique comportant une pluralité de mailles,chaque maille dudit modèle étant associée à au moins un paramètre numérique représentant une propriété géologique. Pour au moins une colonne d'au moins n mailles du modèle, il est possible de déterminer une pluralité de décompositions de ladite colonne en m sous-ensembles de mailles adjacentes. Puis, pour chaque décomposition déterminée, on détermine, pour chacun des m sous-ensembles, un nouveau paramètre numérique associé audit sous-ensemble en fonction des paramètres numériques associés aux mailles dudit sous-ensemble. Enfin, on peut déterminer une valeur de dispersion de ladite décomposition en fonction,au moins, des nouveaux paramètres numériques associés aux dits sous-ensembles.

Description

PROCEDE D'ASSISTANCE A LA MODELISATION GEOLOGIQUE PAR
REGROUPEMENTS DE MAILLES
La présente invention concerne le domaine de la simulation informatisée de données dans un modèle maillé, notamment le domaine de la simulation informatisée de données représentant des informations géologiques.
Dans le cadre d'exploitation pétrolière, l'amélioration du taux de récupération des hydrocarbures dans un gisement nécessite de connaître le plus précisément possible la structure géologique des sols contenant ce gisement. De nombreuses informations sont recueillies sur le terrain au cours de la vie du champ pétrolifère, puis sont analysées et enfin combinées pour la création de modèles géologiques numériques.
Le nombre de mailles de tels modèles peut avoisiner plusieurs millions lors de la création de modèles fins. Avec les technologies actuelles, il n'est pas toujours possible d'utiliser ces modèles fins lors de simulations dynamiques, les simulateurs d'écoulement ne parvenant que difficilement à manipuler des modèles excédant 500.000 mailles.
Ainsi, les géologues et ingénieurs-réservoir simplifient ces modèles pour diminuer le nombre de mailles tout en préservant au mieux les informations qu'ils possèdent : des modèles grossiers sont alors obtenus.
Cependant, de telles simplifications ne sont pas exemptes de défauts.
Cette étape de simplification est difficile car certaines simplifications utiles pour les géologues compromettent d'autres simplifications souhaitées par les ingénieurs- réservoir, et réciproquement. En effet, les géologues proposent des modèles en regroupant certaines mailles selon leur propre jugement ou analyse. Ces modèles sont utilisés pour les évaluations de la capacité de production puis validés (ou invalidés) par les ingénieurs-réservoir. Cette simplification est donc très fastidieuse et nécessite de nombreux échanges entre les géologues et les ingénieurs-réservoir afin d'obtenir le meilleur compromis envisageable. Il y a ainsi un besoin pour simplifier de manière efficace les modèles utilisés pour de telles simulations.
La présente invention vient améliorer la situation.
A cet effet, la présente invention propose de combiner les contraintes des géologues et des ingénieurs-réservoir et de proposer plusieurs solutions optimisées de regroupements de mailles afin de diminuer le nombre de mailles final du modèle.
La présente invention vise alors un procédé de regroupements de mailles dans un modèle géologique maillé de une à trois dimensions, le modèle géologique comportant une pluralité de mailles, et chaque maille dudit modèle étant associée à au moins un paramètre numérique représentant une propriété géologique.
Le procédé comporte les étapes : pour au moins une colonne d'au moins n mailles du modèle:
a/ déterminer une pluralité de décompositions de la colonne en m sous- ensembles de mailles adjacentes, m étant un entier positif ou nul inférieur ou égal à n ;
pour chaque décomposition déterminée :
a1/ pour chacun des m sous-ensembles, déterminer un nouveau paramètre numérique associé audit sous-ensemble en fonction des paramètres numériques associés aux mailles dudit sous- ensemble ;
a2/ déterminer une valeur de dispersion de ladite décomposition en fonction, au moins, des nouveaux paramètres numériques associés auxdits sous-ensembles déterminés à l'étape a1/ ;
b/ fournir lesdites décompositions associées à leur valeur de dispersion. n est un entier strictement positif.
On appelle « propriété géologique » toute propriété physique d'un sol ayant trait soit à sa composition, soit à sa structure, soit à son évolution, etc. A titre d'illustration, cette propriété peut correspondre au type de faciès majoritaire du sol correspondant à ladite maille. En outre, cette propriété peut correspondre à la porosité du sol, à la présence ou non de fracture, à la direction d'écoulement, etc. du sol correspondant à ladite maille.
On appelle « colonne » tout ensemble de mailles consécutives du modèle le long d'une direction de ce même modèle. Ainsi, si le modèle est maillé de manière cubique et de manière régulière, le modèle étant orienté selon un repère orthonormé
(0, x,y,z ), une colonne peut être constituée selon la direction z d'un ensemble de mailles cubiques « empilées» les unes sur les autres, verticalement. Il peut être avantageux que la position de cette colonne dans le modèle corresponde à la position d'un puits de forage.
On appelle « un sous-ensemble de mailles adjacentes » un ensemble de mailles disposant chacun d'au moins une face commune avec une des autres mailles.
Ainsi, la décomposition déterminée ne comporte pas de sous-ensemble de mailles disjointes et peut permettre la création ultérieure d'un nouveau modèle ayant des mailles grossières correspondant à ces sous-ensembles.
Le « nouveau paramètre numérique d'un sous-ensemble » peut être déterminé, par exemple, en réalisant une moyenne, une moyenne pondérée, un seuillage, etc. des différents paramètres numériques associés aux mailles de ce sous-ensemble.
La « valeur de dispersion d'une décomposition » peut être déterminée, par exemple, en calculant la variance des nouveaux paramètres numériques des sous- ensembles de cette décomposition.
Une fois les décompositions fournies en association avec leurs valeurs de dispersion, il est possible de réaliser un grand nombre d'opérations supplémentaires. En effet, le géologue et l'ingénieur-réservoir essayent d'optimiser un certain nombre de paramètres parfois antagonistes et la mise à disposition de ces décompositions ainsi que de leurs valeurs de dispersion peut leur permettre de réaliser un tel arbitrage. Cet arbitrage peut être, à titre d'illustration, un arbitrage entre les temps de simulation liés aux nombres de mailles du modèle et la précision de la simulation liée à la proximité du modèle avec la réalité géologique (pouvant dépendre de la valeur de dispersion). De plus, le procédé peut comprendre en outre :
c/ ordonner lesdites décompositions fournies à l'étape b/ selon la valeur de dispersion déterminée pour chacune des décompositions ;
d/ fournir une liste ordonnée desdites décompositions.
Le classement des décompositions selon la valeur de dispersion déterminée peut permettre à l'opérateur final de simplifier son arbitrage en lui permettant de sélectionner plus rapidement la décomposition cible en fonction de ses objectifs (par exemple, la rapidité de simulation et/ou la précision).
Avantageusement, les étapes pouvant être exécutées pour une pluralité de colonnes, les colonnes pouvant être distinctes et parallèles deux à deux, chaque décomposition déterminée pour une colonne de la pluralité de colonnes peut disposer d'une décomposition homologue dans chacune des autres colonnes. Alors, le procédé peut comprendre en outre :
el pour chaque décomposition des colonnes de la pluralité de colonnes, corriger les valeurs de dispersion déterminées en fonction de la valeur de dispersion de la décomposition et des valeurs de dispersion des décompositions homologues des autres colonnes.
En effet, il n'est pas rare que de nombreux puits de forage (d'exploitation ou de prospection) soient réalisés dans une zone géographique limitée. Ainsi, un modèle peut comprendre plusieurs colonnes représentant ces puits de forages.
On appelle « décomposition homologue à une décomposition d'une colonne » une décomposition de représentation identique à une décomposition dans une autre colonne. Par exemple, une décomposition dans une première colonne peut regrouper les trois premières mailles de cette colonne dans un premier sous- ensemble, puis les dix mailles suivantes dans un deuxième sous-ensemble, etc. La décomposition homologue à cette décomposition comprend alors les mêmes regroupements, mais sur une deuxième colonne distincte : les trois premières mailles de cette deuxième colonne dans un premier sous-ensemble, puis les dix mailles suivantes de la deuxième colonne dans un deuxième sous-ensemble, etc. Cette correction peut permettre de mieux rendre compte de la réalité géologique du sol et ainsi permettre de réaliser de meilleurs arbitrages entre différents maillages grossiers. En effet, il peut être possible de choisir une décomposition (et donc un nouveau maillage d'un modèle grossier) parmi les décompositions ayant la meilleure valeur de dispersion et ceci pour une majorité de colonnes parmi les colonnes du modèle.
Cette correction peut consister dans le fait d'attribuer une même valeur de dispersion à toutes les colonnes pour les décompositions homologues. Cette même valeur peut, par exemple, être la moyenne des valeurs de dispersion des décompositions homologues des différentes colonnes ou être une des valeurs de dispersion des colonnes.
En outre, chaque colonne comportant un paramètre d'importance, la valeur de dispersion de chaque décomposition peut être fonction dudit paramètre d'importance de la colonne.
On appelle « paramètre d'importance » un paramètre pouvant permettre de pondérer l'importance de certains puits et donc des colonnes de mailles associées à ces puits de forage. Par exemple, ce paramètre d'importance peut être fonction du débit du puits, de la prévision de capacité de production du puits, de la taille du puits (puits principal, puits secondaires, etc.).
En effet, il est possible de favoriser la modélisation de certains puits en attribuant aux colonnes correspondantes un paramètre d'importance (paramètre numérique, par exemple décimal) de valeurs différentes. Alors, la valeur de dispersion préalablement calculée peut être pondérée par ce paramètre d'importance.
Dans un mode de réalisation de l'invention, toutes les décompositions possibles de la colonne peuvent être déterminées à l'étape a/.
Ainsi, il est possible d'explorer de manière exhaustive, sans a priori ou postulat, l'ensemble des décompositions possibles de la colonne. Avantageusement, chaque maille de la colonne peut être associée à une valeur (par exemple numérique) de complétion représentant un degré de complétion de la maille.
Ainsi, la pluralité de décompositions peut être déterminée à l'étape a/ de sorte que la valeur de complétion de chaque maille de chacun des m sous-ensembles peut être différente des valeurs de complétion des autres mailles du même sous-ensemble d'au plus un seuil déterminé.
Par exemple, si la décomposition comporte deux regroupements (le premier sous- ensemble correspondant aux dix premières mailles de la colonne et le deuxième sous-ensemble correspondant aux autres mailles de la colonne), l'écart de valeurs de complétion maximal entre deux mailles du premier sous-ensemble n'excède pas ce seuil.
En outre, le seuil déterminé peut être inférieur à la moitié de la différence entre la plus grande valeur de complétion et la plus petite valeur de complétion associées aux mailles de la colonne.
Ainsi, si les valeurs de complétion des mailles de la colonne sont comprises entre 0 et 10, le seuil est inférieur à 5. Ainsi, si les mailles ont comme valeurs de complétion respectives (dans l'ordre de la colonne) 1 , 3, 2, 9, 8, 1 , 5, 3, alors la décomposition peut correspondre au regroupement des trois premières mailles (car l'écart des valeurs 1 , 3, 2 est au maximum de 2, soit un écart inférieur à 5), des deux suivantes (car l'écart des valeurs 9, 8 est au maximum de 1 , soit un écart inférieur à 5) et des trois suivantes (car l'écart des valeurs 1 , 5, 3 est au maximum de 4, soit un écart inférieur à 5).
En revanche, la décomposition ne peut pas correspondre au regroupement des 4 premières mailles (car l'écart des valeurs 1 , 3, 2, 9 est au maximum de 8, soit un écart supérieur à 5), et des quatre mailles suivantes (car l'écart des valeurs 8, 1 , 5, 3 est au maximum de 7, soit un écart supérieur à 5). Avantageusement, le modèle géologique peut être un modèle maillé à deux ou trois dimensions, et chaque maille de chaque colonne peut être associée à une valeur de complétion représentant un degré de complétion de la maille. En outre, chaque maille de la première colonne peut être associée à une valeur de complétion associée à une maille de la deuxième colonne, dite valeur de complétion virtuelle.
Ainsi, la pluralité de décompositions pour la première colonne peut être déterminée à l'étape a/ de sorte que la valeur de complétion virtuelle de chaque maille de chacun des m sous-ensembles peut être différente des valeurs de complétion virtuelle des autres mailles du même sous-ensemble d'au plus un seuil déterminé.
Ainsi, il est possible d'associer à une maille de la première colonne la valeur de complétion d'une maille correspondante dans une deuxième colonne. Il est possible de prendre une maille de la deuxième colonne ayant une même coordonnée avec la maille de la première colonne. Par exemple, dans un modèle à deux dimensions dont l'espace est défini par le repère orthonormé (0, x,z ), et si les colonnes sont orientées selon l'axe z , il est possible d'associer une maille de la première colonne avec une maille de la deuxième colonne si ces deux mailles possèdent une même coordonnée selon l'axe z (i.e. si ces deux mailles se trouvent à la « même profondeur » dans la colonne).
Avantageusement, l'utilisation de contraintes extérieures à la colonne peut permettre de s'assurer que les décompositions déterminées pour une colonne du modèle sont compatibles avec les décompositions déterminées pour les autres colonnes du modèle. Ainsi, il est possible de s'assurer que les colonnes possèdent au moins une décomposition homologue avec les autres colonnes.
Dans un mode de réalisation de l'invention, un des seuils déterminés ci-dessus peut être égal à 0.
Ainsi, dans un sous-ensemble, les valeurs de détermination des mailles peuvent être identiques. En outre, le seuil déterminé peut être inférieur à la moitié de la différence entre la plus grande valeur de complétion virtuelle et la plus petite valeur de complétion virtuelle associées aux mailles de la première colonne.
Dans un mode de réalisation, la valeur de dispersion d'une décomposition peut être déterminée en fonction de la variance interclasse des sous-ensembles déterminés.
On appelle « variance interclasse des sous-ensembles » la variance de l'ensemble constitué par les nouveaux paramètres numériques associés aux sous- ensembles d'une décomposition calculés à l'étape a1/.
Lorsque la variance interclasse est maximale, cela peut signifier que l'hétérogénéité des valeurs de mailles a été respectée lors des regroupements.
Avantageusement, chaque maille de la colonne pouvant être associée à une valeur (par exemple, numérique) d'hétérogénéité représentant un paramètre géologique du sol, la pluralité de décompositions peut être déterminée à l'étape a/ de sorte que la valeur d'hétérogénéité de chaque maille de chacun des m sous- ensembles peut être différente des valeurs d'hétérogénéité des autres mailles du même sous-ensemble d'au plus un seuil déterminé.
On appelle « valeur d'hétérogénéité », une valeur numérique représentant un paramètre géologique dont on souhaite préserver l'hétérogénéité lors des décompositions.
Un dispositif destiné aux regroupements de mailles dans un modèle géologique maillé peut être avantageux, en lui-même.
Ainsi, la présente invention vise également un dispositif de regroupements de mailles dans un modèle géologique maillé de une à trois dimensions, le modèle géologique comportant une pluralité de mailles, chaque maille dudit modèle étant associée à au moins un paramètre numérique représentant une propriété géologique, le dispositif comportant :
- un circuit électronique apte à, pour au moins une colonne d'au moins n mailles du modèle :
- déterminer une pluralité de décompositions de la colonne en m sous- ensembles de mailles adjacentes, m étant un entier positif ou nul inférieur ou égal à n ;
pour chaque décomposition déterminée :
a1/ pour chacun des m sous-ensembles, déterminer un nouveau paramètre numérique associé audit sous-ensemble en fonction des paramètres numériques associés aux mailles dudit sous-ensemble ;
a2/ déterminer une valeur de dispersion de ladite décomposition en fonction, au moins, des paramètres numériques associés auxdits sous-ensembles ;
- une interface de sortie pour fournir lesdites décompositions associées à leur valeur de dispersion, et dans lequel n est un entier strictement positif.
Un programme informatique, mettant en œuvre tout ou partie du procédé décrit ci- avant, installé sur un équipement préexistant, est en lui-même avantageux, dès lors qu'il permet le regroupement de mailles dans un modèle géologique maillé.
Ainsi, la présente invention vise également un programme informatique comportant des instructions pour la mise en œuvre du procédé précédemment décrit, lorsque ce programme est exécuté par un processeur.
La figure 4 décrite en détails ci-après, peut former l'organigramme de l'algorithme général d'un tel programme informatique. D'autres caractéristiques et avantages de l'invention apparaîtront encore à la lecture de la description qui va suivre. Celle-ci est purement illustrative et doit être lue en regard des dessins annexés sur lesquels :
- les figures 1 a et 1 b illustrent un remaillage d'un modèle géologique à trois dimensions ;
- la figure 2 illustre un modèle géologique à deux dimensions ;
- la figure 3 illustre une détermination d'un nouveau paramètre numérique de sous-ensemble et de détermination de valeur de dispersion dans un mode de réalisation de l'invention ;
- la figure 4 illustre un organigramme d'un procédé possible dans un mode de réalisation de l'invention ;
- la figure 5 illustre des regroupements contraints par des valeurs de complétion de mailles dans un mode de réalisation de l'invention ;
- la figure 6 illustre un dispositif pouvant mettre en œuvre le procédé de l'invention dans un mode de réalisation particulière.
Les figures 1 a et 1 b illustrent un remaillage d'un modèle géologique à trois dimensions.
Le modèle géologique de la figure 1 a (dit modèle fin) est représenté par un cube 100a maillé de 216 mailles (i.e. 6x6x6 mailles). Chaque maille, dite maille fine (i.e. 101 , 102, 103, 104, etc.) possède des coordonnées dans l'espace formé par le repère orthonormé (0, x,y,z ).
Néanmoins, ce grand nombre de mailles peut constituer un handicap lors de calculs mathématiques poussés sur le modèle géologique, comme des simulations d'écoulements de fluides ou de déformations géologiques. Il est ainsi avantageux de regrouper certaines mailles, par exemple, par groupe de 8 (e.g. 2x2x2) ou encore par groupe de 4 (e.g. 4x1 x1 ).
La figure 1 b illustre un tel regroupement. Le modèle géologique (dit modèle grossier) de la figure 1 b est représenté par un cube 100b maillé de 27 mailles (i.e. 3x3x3 mailles). Chaque maille, dite maille grossière (i.e. 106, etc.) possède des coordonnées dans l'espace formé par le repère orthonormé (0, x,y,z ) et peut être mise en correspondance avec 8 mailles fines du modèle fin. Le regroupement des 8 mailles fines correspondant au cube 106 de la figure 1 b est mis en évidence par le cube 105 de la figure 1 a.
La figure 2 illustre un modèle géologique à deux dimensions. Ce modèle géologique comporte 32 mailles (i.e. 4x8 mailles) ayant des coordonnées dans un espace formé, le repère orthonormé (Ο, χ,ζ ). Ces mailles (201 ,202, etc.) comportent chacune un paramètre géologique (Goo, G01 , G-io, Gn , G20, etc.). Ce paramètre géologique peut notamment correspondre :
- au type de faciès de la roche majoritaire dans la maille considérée,
- à la perméabilité de la maille,
- à la présence de failles dans la maille,
- à un paramètre d'écoulement de la maille,
- à un paramètre de porosité de la maille,
- à la concentration en calcite de la maille, etc.
Ce paramètre géologique (Goo, G01, G10, Gn, G20, etc.) est un paramètre numérique, par exemple décimal.
La figure 3 illustre une détermination d'un nouveau paramètre géologique de sous- ensemble dans un modèle à une dimension et selon un mode de réalisation de l'invention.
L'ensemble 351 des cinq mailles (300 à 304) constitue une colonne initiale et correspond au modèle à une dimension avant toute décomposition en sous- ensembles. Chaque maille de cette colonne initiale est associée à un paramètre géologique (Go est associé à la maille 300, G4 est associé à la maille 304).
L'ensemble 352 correspond à une décomposition possible de la colonne 351 en quatre sous-ensembles (305 à 308). Les sous-ensembles 306, 307 et 308 englobent respectivement les mailles 302, 303 et 304, tandis que le sous-ensemble 305 englobe les deux mailles 300 et 301 . Pour chacun de ces sous-ensembles (305 à 308), il est possible de déterminer un nouveau paramètre numérique géologique associé à ce sous-ensemble. Par exemple :
- le nouveau paramètre géologique associé au sous-ensemble 305 est fonction des paramètres Go et d (noté f(Go, Gi )) ;
- le nouveau paramètre géologique associé au sous-ensemble 306 est fonction du paramètre G2 (noté f(G2)) ;
- le nouveau paramètre géologique associé au sous-ensemble 307 est fonction du paramètre Gs (noté f(Gs)) ;
- le nouveau paramètre géologique associé au sous-ensemble 308 est fonction du paramètre G4 (noté f(G4)).
Il est également possible de réaliser d'autres décompositions de la colonne 351 . Ainsi, l'ensemble 353 correspond à une autre décomposition de la colonne 351 , en trois sous-ensembles (309 à 31 1 ). Les sous-ensembles 309 et 31 1 englobent respectivement les mailles 300 et 304, tandis que le sous-ensemble 310 englobe les trois mailles 301 , 302 et 303. Pour chacun de ces sous-ensembles (309 à 31 1 ), il est possible de déterminer un nouveau paramètre numérique géologique associé à ce sous-ensemble. Par exemple :
- le nouveau paramètre géologique associé au sous-ensemble 309 est fonction du paramètre Go (noté f(Go)) ;
- le nouveau paramètre géologique associé au sous-ensemble 310 est fonction des paramètres d, G2 et Gs (noté f(d, G2, G3)) ;
- le nouveau paramètre géologique associé au sous-ensemble 31 1 est fonction du paramètre G4 (noté f(G4)).
La fonction f permettant la détermination du nouveau paramètre peut être, par
Q
exemple une fonction moyenne exprimée de la manière suivantef(G1,...,Gn ) = ^—i- .
i f
La fonction f peut être également une fonction moyenne pondérée exprimée de la
Q
manière suivantef(G1,...,Gn ) = ^ki— i- , dans laquelle k, peut être un facteur i ^
géologique de la maille 30i comme la perméabilité. Il est également possible de noter G; ce nouveau paramètre associé au sous-ensemble d'indice j (en conférant, le cas échéant, à chaque sous-ensemble d'une décomposition, un indice). Par ailleurs, une valeur de dispersion peut être déterminée pour chacune des décompositions 352 et 353. La valeur de dispersion permet notamment de quantifier la perte d'hétérogénéité de la colonne lors d'une décomposition.
Ainsi, cette valeur de dispersion peut correspondre à la valeur de la variance interclasse au sein de la colonne. Ainsi, une telle variance peut être calculée de la
m [
manière suivante σ2 = (G| - Gtot )2— avec
i=i h tôt
- m le nombre de sous-ensembles déterminés (i.e. le nombre de mailles grossières cible),
- G; le nouveau paramètre géologique associé au sous-ensemble d'indice i (comme décrit précédemment par exemple),
- Gtot la valeur d'un nouveau paramètre géologique associé à un sous- ensemble comportant toutes les mailles de la colonne (par exemple, la moyenne des valeurs géologiques des mailles de la colonne),
- h| le nombre de mailles contenues dans le sous-ensemble d'indice i, et
- htot le nombre total de mailles dans la colonne.
Ainsi, plus la valeur de dispersion est élevée, plus l'hétérogénéité de la colonne est respectée. Nous pouvons écrire σ352 2 la variance interclasse de la décomposition
352 et σ353 2 la variance interclasse de la décomposition 353. Alors, dans le cas particulier de la figure 3, nous avons : σ352 2 o353 2
Figure imgf000015_0001
En supposant que Gi=G2=G3=1 et que Go=G4 =0 et que la fonction f correspond à la fonction moyenne, nous avons :
Figure imgf000015_0002
' 353 0 -- + - 1 -- + - 0 -- 0,36
La variance σ353 est la variance maximale possible à obtenir lors d'une décomposition avec ces hypothèses.
Ainsi, en ordonnant les différentes décompositions selon la valeur de dispersion associée, il est possible de connaître les décompositions respectant le plus l'hétérogénéité de la colonne (c'est-à-dire, ici, la décomposition 353).
La figure 4 illustre un organigramme d'un procédé possible dans un mode de réalisation de l'invention.
Sur réception d'un modèle géologique maillé (400), il est possible de déterminer un grand nombre de décompositions de colonnes en sous-ensembles de mailles adjacentes (étape 401 ).
Le modèle reçu comprend, par exemple, des mailles, chaque maille étant cubique et chaque côté de ce cube étant commun avec un autre cube. De plus, ce modèle peut également indiquer un ensemble de colonnes de mailles empilées les unes sur les autres, adjacentes et représentant des puits de forage géologique.
Ainsi, pour chacun de ces puits, il est possible de déterminer une pluralité de décompositions. L'objet de telles décompositions est de réduire le nombre de mailles du modèle géologique selon l'axe vertical (i.e. pour des couches successives du modèle). Si toutes les décompositions sont déterminées alors leur nombre est avec i le nombre initial de mailles dans la colonne, et f le nombre
Figure imgf000016_0001
de sous-ensembles dans la décomposition (i.e. le nombre de mailles dans le modèle cible).
S'il n'a pas été déterminé de valeur de dispersion pour chacune des décompositions déterminées précédemment (sortie OK du test 403), alors il est sélectionné une décomposition pour laquelle aucune valeur de dispersion n'a été déterminée. Alors, pour chacun des sous-ensembles de cette décomposition, un nouveau paramètre géologique associé au sous-ensemble est calculé (étape 403), par exemple, en réalisant une moyenne des paramètres géologiques des mailles contenues dans ce sous-ensemble.
Une valeur de dispersion associée à la décomposition sélectionnée est alors calculée (étape 404), par exemple en calculant la variance des moyennes précédemment calculées pour les sous-ensembles de cette décomposition.
Une fois toutes les valeurs de dispersion déterminées pour chacune des décompositions (sortie KO du test 403), il est possible alors d'ordonner ces décompositions en fonction de leur valeur respective de dispersion (étape 405).
Enfin cette liste ordonnée peut être fournie (message 406) à un utilisateur ou un autre système de traitement pour affichage ou de nouveaux calculs. En effet, l'utilisateur (par exemple, un géologue ou un ingénieur-réservoir) peut alors sélectionner la décomposition qui corresponde le plus à la réalité géologique constatée parmi les décompositions déterminées ayant les plus hautes valeurs de dispersion.
La figure 5 illustre des regroupements contraints par des valeurs de complétion de mailles dans un mode de réalisation de l'invention.
En effet, comme cela a été abordé précédemment, le nombre total de décompositions d'une colonne peut être relativement élevé et il peut être fastidieux de toutes les analyser. Ainsi, il est possible de fixer un certain nombre de contraintes lors de la décomposition afin de réduire le nombre de décompositions possibles.
A titre d'illustration, il est possible de fixer la taille maximale (en nombre de mailles) d'un sous-ensemble afin d'éviter notamment des phénomènes de « transfert de volume » lors de simulations. Il est également possible de « marquer » certaines mailles afin de rendre leur regroupement avec des mailles marquées de manière différente impossible ou difficile. Ce type de contrainte est le plus souvent souhaité dans le cas de contraintes dites de « complétion ». On entend par « contraintes de complétion » toute contrainte liée à l'exploitation du puits de forage, comme par exemple, l'ouverture dans le tubage (ou « casing » en anglais) d'ouvertures permettant l'exploitation de couches de pétrole, de gaz, etc.
Par exemple, lorsque le tubage est ouvert au niveau d'une maille donnée (i.e. C=1 en référence à la figure 5), il n'est pas souhaitable que cette maille soit regroupée avec une maille pour laquelle le tubage n'est pas ouvert (i.e. C=0). En référence au puits P1 de la figure 5, la décomposition proposée (R1 , R2, R3, R4) fait apparaître certains regroupements possibles tandis que d'autres regroupements ne sont pas possibles :
- le regroupement R1 est possible car les mailles 500 et 501 possèdent une même valeur de complétion C=1 ;
- le regroupement R2 est possible car les mailles 502 et 503 possèdent une même valeur de complétion C=0 ;
- le regroupement R3 n'est pas possible car les mailles 504 et 505 ne possèdent pas une valeur de complétion identique ;
- le regroupement R4 est possible car les mailles 506 et 507 possèdent une même valeur de complétion C=0 ;
En référence au puits P2 de la figure 5, il est possible de réaliser un raisonnement similaire :
- le regroupement R1 ' n'est pas possible car les mailles 510 et 51 1 ne possèdent pas une valeur de complétion identique ;
- le regroupement R2' est possible car les mailles 512 et 513 possèdent une même valeur de complétion C=0 ;
- le regroupement R3' est possible car les mailles 514 et 515 possèdent une même valeur de complétion C=1 ;
- le regroupement R4' est possible car les mailles 516 et 517 possèdent une même valeur de complétion C=0 ;
Ainsi, chaque maille des puits P1 et P2 est associée à une valeur numérique de complétion représentant un degré de complétion de la maille. Bien entendu, cette valeur de complétion n'est pas nécessairement 1 ou 0. Elle peut être également égale à 0,9 ou 0,2 ou tout autre valeur numérique.
Il est également possible d'interdire le regroupement de deux mailles si les valeurs de complétion de ces deux mailles diffèrent de plus d'un seuil Δ . En généralisant à un regroupement de plusieurs mailles, la valeur de complétion de chaque maille de ce regroupement (ou sous-ensemble) ne peut pas être différente des valeurs de complétion des autres mailles de ce regroupement d'au plus un seuil Δ .
Il peut être avantageux de réaliser les regroupements de mailles de façon identique entre les différentes colonnes afin de pouvoir construire un maillage régulier au sein d'un modèle grossier. Pour ce faire, les contraintes imposées au sein d'une colonne peuvent être répercutées sur les autres colonnes du modèle. Toujours en référence à la figure 5, les colonnes de puits P1 et P2 comportent des mailles (respectivement 500 à 507 et 510 à 517) associées chacune à une valeur de complétion. Il est possible de mettre en correspondance les mailles de ces deux puits en fonction de leur « profondeur » (i.e. leur position dans leur colonne respective). Ainsi, la maille 500 peut être mise en correspondance avec la maille 510, la maille 501 peut être mise en correspondance avec la maille 51 1 ... et la maille 507 peut être mise en correspondance avec la maille 517.
Afin d'assurer un maillage grossier régulier, il est possible d'associer aux mailles de la première colonne les valeurs de complétion de la deuxième colonne en fonction de la correspondance décrite ci-dessus. Ces valeurs sont alors appelées « valeur de complétion virtuelle ». A titre d'illustration, la valeur de complétion virtuelle de la maille 504 est C=1 (car la maille correspondante à la maille 504 est la maille 514 dans la deuxième colonne et que sa valeur de complétion vaut C=1 ). De même, la valeur de complétion virtuelle de la maille 501 est C=0 (car la maille correspondante à la maille 501 est la maille 51 1 dans la deuxième colonne et que sa valeur de complétion vaut C=0).
Le processus décrit ci-dessus peut alors être réalisé en transposant les valeurs de complétion avec les valeurs de complétion virtuelles.
Ainsi, pour vérifier si un regroupement est possible entre plusieurs mailles d'une colonne, il est possible de vérifier que la valeur de complétion virtuelle de chaque maille de ce regroupement est différente des valeurs de complétion virtuelle des autres mailles de ce même regroupement d'au plus un seuil Δ .
Le seuil Δ peut être fixé à l'avance dans un fichier de configuration ou dynamiquement en fonction des résultats d'un premier processus de décomposition. Une solution intermédiaire peut consister en une fixation du seuil Δ à une valeur numérique initiale assez basse (par exemple, 0) puis en une modification incrémentale à la hausse de sa valeur, si le nombre de décompositions déterminées n'est pas assez élevé.
La figure 6 représente un exemple de dispositif 600 de regroupements de mailles dans un modèle géologique maillé dans un mode de réalisation de l'invention.
Dans ce mode de réalisation, le dispositif 600 comporte un ordinateur, comprenant une mémoire 604 pour stocker des instructions permettant la mise en œuvre du procédé, les données de modèles reçues, et des données temporaires pour réaliser les différentes étapes du procédé tel que décrit précédemment.
L'ordinateur comporte en outre un circuit de regroupements de mailles 601 . Ce circuit peut être, par exemple :
- un processeur apte à interpréter des instructions sous la forme de programme informatique, ou
- une carte électronique dont les étapes du procédé de l'invention sont décrites dans le silicium, ou encore
- une puce électronique programmable comme une puce FPGA (pour « Field- Programmable Gâte Array » en anglais).
Cet ordinateur comporte une interface d'entrée 602 pour la réception de données de modèles géologiques maillés, et une interface de sortie 603 pour la fourniture de données de décompositions. Enfin, l'ordinateur comporte, pour permettre une interaction aisée avec un utilisateur, un écran 605 et un clavier 606.
Par ailleurs, le schéma fonctionnel présenté sur la figure 4 est un exemple typique d'un programme dont certaines instructions peuvent être réalisées auprès du dispositif décrit. A ce titre, la figure 4 peut correspondre à l'organigramme de l'algorithme général d'un programme informatique au sens de l'invention. Bien entendu, la présente invention ne se limite pas aux formes de réalisation décrites ci-avant à titre d'exemples ; elle s'étend à d'autres variantes.
D'autres réalisations sont possibles.
Par exemple, toutes les explications données et tous les enseignements fournis pour des modèles à une, deux ou trois dimensions sont transposables aux modèles de dimensions autres parmi les modèles à une, deux ou trois dimensions.

Claims

REVENDICATIONS
1 . Procédé d'assistance à la modélisation d'un sous-sol géologique dans un modèle géologique maillé de une à trois dimensions, le modèle géologique (100a) comportant une pluralité de mailles (101 , 102, 103, 104), chaque maille (101 , 102, 103, 104) dudit modèle étant associée à au moins un paramètre numérique (Goo, G01 , G02) représentant une propriété géologique, le procédé comportant les étapes : pour au moins une colonne (351 ) d'au moins n mailles du modèle :
a/ déterminer (401 ) une pluralité de décompositions (352, 353) de ladite colonne en m sous-ensembles (305, 306, 307, 308, 309, 310, 31 1 ) de mailles adjacentes, m étant un entier positif ou nul inférieur ou égal à n ; pour chaque décomposition déterminée :
a1/ pour chacun des m sous-ensembles, déterminer (403) un nouveau paramètre numérique associé audit sous-ensemble en fonction des paramètres numériques associés aux mailles dudit sous- ensemble ;
a2/ déterminer une valeur de dispersion (404) de ladite décomposition en fonction, au moins, des nouveaux paramètres numériques associés auxdits sous-ensembles déterminés à l'étape a1/ ;
bl fournir lesdites décompositions associées à leur valeur de dispersion, dans lequel n est un entier strictement positif.
2. Procédé selon la revendication 1 , dans lequel le procédé comprend en outre : c/ ordonner lesdites décompositions (405) fournies à l'étape b/ selon la valeur de dispersion déterminée pour chacune des décompositions ;
d/ fournir une liste ordonnée (406) desdites décompositions.
3. Procédé selon l'une des revendications précédentes, dans lequel, les étapes étant exécutées pour une pluralité de colonnes (P1 , P2), les colonnes étant distinctes et parallèles deux à deux, chaque décomposition déterminée pour une colonne de la pluralité de colonnes disposant d'une décomposition homologue dans chacune des autres colonnes ; le procédé comprend en outre :
el pour chaque décomposition des colonnes de la pluralité de colonnes, corriger les valeurs de dispersion déterminées en fonction de la valeur de dispersion de la décomposition et des valeurs de dispersion des décompositions homologues des autres colonnes.
4. Procédé selon l'une des revendications précédentes, dans lequel, chaque colonne comportant un paramètre d'importance, la valeur de dispersion de chaque décomposition est fonction dudit paramètre d'importance de la colonne.
5. Procédé selon l'une des revendications précédentes, dans lequel toutes les décompositions possibles de la colonne sont déterminées à l'étape a/.
6. Procédé selon l'une des revendications 1 à 4, dans lequel, chaque maille de la colonne étant associée à une valeur de complétion (C) représentant un degré de complétion de la maille, la pluralité de décompositions est déterminée à l'étape a/ de sorte que la valeur de complétion de chaque maille de chacun des m sous-ensembles est différente des valeurs de complétion des autres mailles du même sous-ensemble d'au plus un seuil déterminé.
7. Procédé selon l'une des revendications 1 à 4, dans lequel, le modèle géologique étant un modèle maillé à deux ou trois dimensions, chaque maille de chaque colonne étant associée à une valeur de complétion représentant un degré de complétion de la maille, chaque maille de la première colonne étant associée à une valeur de complétion associée à une maille de la deuxième colonne, dite valeur de complétion virtuelle, la pluralité de décompositions pour la première colonne est déterminée à l'étape a/ de sorte que la valeur de complétion virtuelle de chaque maille de chacun des m sous-ensembles est différente des valeurs de complétion virtuelle des autres mailles du même sous-ensemble d'au plus un seuil déterminé.
8. Procédé selon l'une des revendications 6 et 7, dans lequel le seuil déterminé est égal à 0.
9. Procédé selon la revendication 6, dans lequel le seuil déterminé est inférieur à la moitié de la différence entre la plus grande valeur de complétion et la plus petite valeur de complétion associées aux mailles de la colonne.
10. Procédé selon la revendication 7, dans lequel le seuil déterminé est inférieur à la moitié de la différence entre la plus grande valeur de complétion virtuelle et la plus petite valeur de complétion virtuelle associées aux mailles de la première colonne.
1 1 . Procédé selon l'une des revendications précédentes, dans lequel la valeur de dispersion d'une décomposition est déterminée en fonction de la variance interclasse des sous-ensembles déterminés.
12. Procédé selon l'une des revendications 1 à 4, dans lequel, chaque maille de la colonne étant associée à une valeur d'hétérogénéité représentant un paramètre géologique du sol, la pluralité de décompositions est déterminée à l'étape a/ de sorte que la valeur d'hétérogénéité de chaque maille de chacun des m sous-ensembles est différente des valeurs d'hétérogénéité des autres mailles du même sous-ensemble d'au plus un seuil déterminé.
13. Produit programme informatique comportant des instructions pour la mise en œuvre du procédé selon l'une des revendications 1 à 1 1 , lorsque ce programme est exécuté par un processeur.
14. Dispositif d'assistance à la modélisation d'un sous-sol géologique dans un modèle géologique maillé de une à trois dimensions, le modèle géologique comportant une pluralité de mailles, chaque maille dudit modèle étant associée à au moins un paramètre numérique représentant une propriété géologique, le dispositif comportant :
- un circuit électronique apte à, pour au moins une colonne d'au moins n mailles du modèle :
- déterminer une pluralité de décompositions de la colonne en m sous- ensembles de mailles adjacentes, m étant un entier positif ou nul inférieur ou égal à n ;
pour chaque décomposition déterminée :
a1/ pour chacun des m sous-ensembles, déterminer un nouveau paramètre numérique associé audit sous-ensemble en fonction des paramètres numériques associés aux mailles dudit sous-ensemble ;
a2/ déterminer une valeur de dispersion de ladite décomposition en fonction, au moins, des nouveaux paramètres numériques associés auxdits sous-ensembles ; - une interface de sortie pour fournir lesdites décompositions associées à leur valeur de dispersion,
dans lequel n est un entier strictement positif.
PCT/FR2013/051900 2012-08-06 2013-08-06 Procede d'assistance a la modelisation geologique par regroupements de mailles WO2014023916A1 (fr)

Priority Applications (3)

Application Number Priority Date Filing Date Title
GB1501860.9A GB2519037B (en) 2012-08-06 2013-08-06 Method of assistance in geological modeling by grouping meshes together
NO20150173A NO346799B1 (en) 2012-08-06 2013-08-06 Modelling a geological subsoil by grouping together mesh cells in a meshed geological model of one to three dimensions
US14/420,153 US10274641B2 (en) 2012-08-06 2013-08-06 Method of assistance in geological modeling by grouping meshes together

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1257647 2012-08-06
FR1257647A FR2994313B1 (fr) 2012-08-06 2012-08-06 Procede d'assistance a la modelisation geologique par regroupements de mailles

Publications (1)

Publication Number Publication Date
WO2014023916A1 true WO2014023916A1 (fr) 2014-02-13

Family

ID=47227960

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/FR2013/051900 WO2014023916A1 (fr) 2012-08-06 2013-08-06 Procede d'assistance a la modelisation geologique par regroupements de mailles

Country Status (5)

Country Link
US (1) US10274641B2 (fr)
FR (1) FR2994313B1 (fr)
GB (1) GB2519037B (fr)
NO (1) NO346799B1 (fr)
WO (1) WO2014023916A1 (fr)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111862326A (zh) * 2020-07-24 2020-10-30 咪咕文化科技有限公司 地质模型的存储方法和计算方法、电子设备及存储介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090299714A1 (en) * 2008-05-30 2009-12-03 Kelkar And Ass0Ciates, Inc. Dynamic updating of simulation models
US20120035896A1 (en) * 2010-08-09 2012-02-09 Conocophillips Company Reservoir upscaling method with preserved transmissibility

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2725794B1 (fr) * 1994-10-18 1997-01-24 Inst Francais Du Petrole Methode pour modeliser la distribution spatiale d'objets geometriques dans un milieu, tels que des failles dans une formation geologique
US6928399B1 (en) * 1999-12-03 2005-08-09 Exxonmobil Upstream Research Company Method and program for simulating a physical system using object-oriented programming
FR2962835B1 (fr) * 2010-07-16 2013-07-12 IFP Energies Nouvelles Methode pour generer un maillage hexa-dominant d'un bassin geometriquement complexe

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090299714A1 (en) * 2008-05-30 2009-12-03 Kelkar And Ass0Ciates, Inc. Dynamic updating of simulation models
US20120035896A1 (en) * 2010-08-09 2012-02-09 Conocophillips Company Reservoir upscaling method with preserved transmissibility

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
DAVID STERN: "Practical Aspects of Scaleup of Simulation Models", JOURNAL OF PETROLEUM TECHNOLOGY, vol. 57, no. 9, September 2005 (2005-09-01), XP055060786, ISSN: 0149-2136, DOI: 10.2118/89032-MS *

Also Published As

Publication number Publication date
FR2994313B1 (fr) 2014-08-29
GB2519037A (en) 2015-04-08
NO20150173A1 (en) 2015-03-06
US10274641B2 (en) 2019-04-30
US20150226877A1 (en) 2015-08-13
FR2994313A1 (fr) 2014-02-07
GB2519037B (en) 2018-08-29
GB201501860D0 (en) 2015-03-18
NO346799B1 (en) 2023-01-16

Similar Documents

Publication Publication Date Title
Hardebol et al. Multiscale fracture network characterization and impact on flow: A case study on the Latemar carbonate platform
EP2581767B1 (fr) Méthode pour construire un maillage d'un réservoir fracturé avec un nombre limité de noeuds dans le milieu matrice
EP2037080B1 (fr) Méthode pour estimer la preméabilité d'un réseau de fractures à partir d'une analyse de connectivité
US8301427B2 (en) Fracture network characterization method
CA2997608C (fr) Mise en correspondance d'historique de production d'hydrocarbures a partir de reservoirs heterogenes
EP2963235B1 (fr) Procede d'exploitation d'un gisement petrolier a partir d'une technique de positionnement des puits a forer
FR3034556A1 (fr) Cartographie d'ouverture correlee a la taille de la fracture pour une determination de porosite et de permeabilite localisees
WO2009098366A2 (fr) Methode de modelisation d'un milieu géologique poreux traversé par un réseau de fractures
EP2680045A1 (fr) Procédé d'exploitation d'un reservoir geologique à partir d'un modèle de réservoir calé au moyen d'un parametrage multi-echelles
FR3045868A1 (fr) Procede pour caracteriser et exploiter une formation souterraine comprenant un reseau de fractures
CA2821099C (fr) Procede d'exploitation d'un reservoir geologique a partir d'un modele de reservoir cale par le calcul d'une loi analytique de distribution conditionnelle de parametres incertains du modele
EP3164743B1 (fr) Procede de determination de caves geologiques
EP3146367B1 (fr) Procédé de détermination d'une carte de hauteur d'hydrocarbure liquide dans un réservoir
CN105093348A (zh) 一种多层油气资源量预测方法
CN107807396A (zh) 地层基质电阻率的确定方法和装置
FR3048300A1 (fr) Representation hybride geocellulaire en 3d de sous-ensembles de reseau de fracture naturelle choisis
Sylvester Turbidite bed thickness distributions: methods and pitfalls of analysis and modelling
FR3034894A1 (fr)
CA2892118A1 (fr) Procede pour determiner une representation d'un reservoir d'hydrocarbures
WO2014023916A1 (fr) Procede d'assistance a la modelisation geologique par regroupements de mailles
FR2994314A1 (fr) Procede de restauration de donnees de puits de forage
Martinius et al. Uncertainty analysis of fluvial outcrop data for stochastic reservoir modelling
FR3029664A1 (fr) Definition de petrofacies non-lineaires pour un modele de simulation de reservoir
WO2009090522A1 (fr) Procede, programme et systeme informatique de construction d'un modele geologique 3d
CA2892241A1 (fr) Procede pour evaluer une quantite d'hydrocarbures dans un gisement

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: 13758946

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 1501860

Country of ref document: GB

Kind code of ref document: A

Free format text: PCT FILING DATE = 20130806

WWE Wipo information: entry into national phase

Ref document number: 1501860.9

Country of ref document: GB

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 14420153

Country of ref document: US

122 Ep: pct application non-entry in european phase

Ref document number: 13758946

Country of ref document: EP

Kind code of ref document: A1