Brevet FR2765007A1 — publié le 19981224

PROCEDE ET APPAREIL DE QUADRILLAGE DE SIMULATION COMPORTANT UN DISPOSITIF DE QUADRILLAGE DE SURFACE STRUCTURE CONCU POUR ETRE UTILISE PAR UN SIMULATEUR DE RESERVOIR Arrière-plan de l'invention La présente invention concerne un appareil et un procédé logiciel fondés sur un poste de travail qui réagit à des données sismiques et à des données de diagraphie reçues d'une formation terrestre en réponse à une activation sismique, permettant de quadriller la formation terrestre en générant un quadrillage "structuré" et en plaçant ce "quadrillage structuré" sur la formation terrestre. Puisque le quadrillage structuré est composé d'une pluralité de blocs de quadrillage, lorsque le quadrillage structuré est placé sur la formation terrestre, le procédé et l'appareil logiciel fondés sur le poste de travail génèrent en outre une pluralité d'informations plus précises correspondant respectivement à la pluralité de blocs de quadrillage du quadrillage structuré. La pluralité d'informations plus précises concerne, par exemple, les propriétés de transmissibilité de la pluralité de blocs de quadrillage du quadrillage structuré. La pluralité d'informations plus précises est fournie en tant qu'entrée à un "simulateur de réservoir" logiciel fondé sur un poste de travail classique. En réponse à la pluralité d'informations plus précises, le "simulateur de réservoir" génère une pluralité correspondante de résultats de simulation (tels que des pressions et des saturations) appartenant respectivement à la pluralité de blocs de quadrillage du quadrillage structuré, la pluralité de résultats de simulation étant placée, SR 15546 US/RS respectivement, sur la pluralité de blocs de quadrillage du quadrillage structuré de sorte qu'un nouveau résultat de simulation est associé à chaque bloc de quadrillage du quadrillage structuré. Les blocs de quadrillage du quadrillage structuré et les nouveaux résultats de simulation associés à chaque bloc de quadrillage sont affichés sur le moniteur d'affichage du poste de travail pour être visualisés par un opérateur du poste de travail. Plus particulièrement, l'objet de la présente invention concerne un logiciel de quadrillage "structuré" amélioré destiné à être utilisé dans le cadre d'un programme de quadrillage de simulation, le programme de quadrillage de simulation étant ci-après appelé "Flogrid". Le logiciel de quadrillage "structuré", placé à l'intérieur du programme de quadrillage de simulation "Flogrid", est conçu pour être exécuté par le processeur du poste de travail. Lorsqu'il est exécuté par le processeur du poste de travail, le logiciel de quadrillage "structuré" génère finalement les "informations de propriété de bloc de quadrillage de formation terrestre plus précises" citées ci-dessus. Ces "informations de propriété de bloc de quadrillage plus précises" sont utilisées par un logiciel simulateur de réservoir pour générer les "résultats de simulation plus précis", tels que des pressions et des saturations, qui sont affichés sur le moniteur d'affichage du poste de travail. Les opérations sismiques sont effectuées près d'un ou de plusieurs puits de forage dans une formation terrestre et une pluralité de données sismiques est obtenue à partir de cette activation sismique. De plus, les opérations de diagraphie sont également effectuées SR 15546 US/RS dans un ou plusieurs puits de forage et des données de diagraphie sont également obtenues à partir des opérations de diagraphie. Les données sismiques et les données de diagraphie sont entrées dans un poste de travail d'un ordinateur dans lequel un programme d'interprétation est exécuté. Le programme d'interprétation de la technique antérieure était composé d'un premier programme parfois appelé "quadrillage" qui générait des données et d'un deuxième programme de simulation, réagissant au premier programme, qui recevait les données du premier programme "de quadrillage" et générait un ensemble de résultats de simulation et affichait des résultats de simulation sur l'affichage du poste de travail, les résultats de la simulation affichés permettant à un opérateur de déterminer les propriétés d'écoulement de la formation terrestre située près d'un ou plusieurs puits de forage forés dans la formation. En particulier, le premier programme "de quadrillage" établit un quadrillage pour chaque horizon de la formation terrestre proche des puits de forage, le quadrillage de chaque horizon comprenant une multitude de cellules individuelles. De plus, le premier programme "de quadrillage" génère des données et d'autres informations pour chacune des cellules individuelles de chaque horizon, les données et les autres informations de chaque cellule étant transmises au deuxième programme de simulation qui utilise des données et des informations pour chaque cellule de quadrillage, reçues du premier programme de quadrillage, pour générer un ensemble de résultats de simulation pour chaque cellule du quadrillage. Un SR 15546 US/RS résultat de simulation est affiché sur l'affichage du poste de travail pour chaque cellule du quadrillage ce qui permet à un opérateur du poste de travail de déterminer les propriétés de production de l'écoulement de chacune des cellules dans la formation terrestre quadrillée située près des puits de forage. Cependant, les efforts continus de mise au point sont centrés sur l'amélioration de la qualité et de la précision des données et d'autres informations générées par le premier programme "de quadrillage". Lorsqu'un ensemble de données plus précises améliorées est reçu par le deuxième programme de simulation, la fonction de simulation mise en pratique par le deuxième programme de simulation sera plus précise et plus complète; et, en conséquence, les résultats de simulation générés par le deuxième programme de simulation seront plus précis et plus complets. Par conséquent, au vu de l'ensemble des résultats de simulation plus précis et plus complets (par exemple les pressions et les saturations) générés par le deuxième programme de simulation, les propriétés d'écoulement associées à chaque cellule ou bloc de quadrillage du quadrillage placé sur la formation terrestre située près des puits de forage peuvent être déterminées avec plus de précision. En conséquence, il existe un besoin d'améliorer le premier programme "de quadrillage" de sorte que des données plus précises sont générées par le premier programme de quadrillage. La simulation informatique de procédés physiques utilisant des équations différentielles partielles nécessite un agencement en damier d'un volume spécifique d'espace en un ensemble de petits blocs. Ces SR 15546 US/RS agencements en damier peuvent être de différents types: 1. Quadrillages "structurés" qui sont des quadrillages "en forme de tartans" ou rectangulaires ou bien des versions de celles-ci ayant subi une distorsion; 2. Quadrillages "semi-structurés", dans lesquels le domaine est d'abord fractionné en un libre de sous- domaines et un quadrillage structuré est placé dans chacun de ceux-ci; 3. Quadrillages "non structurés" qui consistent en des agencements de triangles bidimensionnels ou de tétraèdres tridimensionnels; et 4. Quadrillages "structurés - non structurés" dans lesquels le domaine est d'abord fractionné en un ensemble libre de sous-domaines et dans lequel des quadrillages structurés sont placés dans certains sous- domaines et dans lesquels des quadrillages non structurés sont placés dans d'autres. La plupart des procédés de génération de quadrillages que l'on utilise aujourd'hui ont été mis au point dans le domaine de la conception assistée par ordinateur, en particulier dans l'industrie aérospatiale. Toutefois, en ce qui concerne les quadrillages "non structurés", les simulations d'écoulement sur des quadrillages fondés sur des triangles (c'est-à-dire pour des quadrillages "non structurés") ont été utilisées par plusieurs personnes dans le cadre de l'industrie pétrolière ou non. Des volumes de contrôle ont été formés autour de chaque noeud d'un quadrillage triangulaire en joignant les points bissecteurs des côtés au centre de gravité du triangle pour résoudre un SR 15546 US/RS problème magnétostatique à deux dimensions. Cette technique est connue sous le nom de procédé d'élément fini de volume de contrôle (CVFE). De plus, des volumes de contrôle ont été formés en joignant les bissectrices perpendiculaires des côtés du triangle d'une triangulation de Delaunay pour résoudre des équations de dispositifs semi-conducteurs. Cette technique a été appliquée à la simulation de réservoir, qui est connue sous le nom de méthode de Voronoi ou PEBI. Les problèmes hétérogènes ont été traités en faisant en sorte que la perméabilité soit constante sur un triangle. D'autres procédés différents ont été utilisés sur la base de la méthode CVFE dans laquelle les perméabilités sont définies comme étant constantes dans les volumes de contrôle. Cette approche traite mieux des limites de couches ayant des différences de perméabilité importantes que la méthode CVFE traditionnelle et comme avec la méthode traditionnelle, ceci donne un motif d'écoulement multipoint, appelé par conséquent procédé MPFA. Au contraire, la méthode PEBI se réduit à un motif d'écoulement de deux points. La méthode PEBI a été étendue pour traiter des systèmes hétérogènes anisotropes en définissant la perméabilité comme étant constante dans un triangle et en réglant l'angle entre les côtés du triangle et les limites de cellules. Cette approche présente deux inconvénients; tout d'abord, le fait de traiter des couches ayant des perméabilités différentes n'est pas souhaitable, ensuite, dans des systèmes hautement anisotropes, la condition d'angle entre les côtés du triangle et les limites de cellules peuvent devenir très défavorables de telle sorte que les cellules peuvent commencer à se SR 15546 US/RS chevaucher. Outre le fait d'utiliser des volumes de contrôle formés autour des noeuds des triangulations, il est possible d'utiliser les triangles eux-mêmes en tant que volumes de contrôle. Un inconvénient des volumes de contrôle triangulaires par rapport aux volumes de Voronoi est que le premier volume nécessite un nombre de cellules bien supérieur; pour des répartitions de points aléatoires, un facteur moyen de deux et cinq existe en deux et trois dimensions, respectivement. Un avantage des quadrillages triangulaires est la souplesse à s'adapter à des caractéristiques de production et géologiques complexes. Un logiciel de quadrillage "non structuré" amélioré est décrit dans la demande de brevet américain en suspens antérieure numéro de série 08/873.234 déposée le 11 juin 1997 et intitulée "Method and Apparatus for generating more accurate earth formation grid cell property information for use by a simulator to display more accurate simulation results of the formation near a wellbore" attribué à Dayal Gunasekera (appelée ci-après, la demande de Gunasekera). Le programme de quadrillage de simulation "Flogrid" de la présente invention décrit dans cette demande comporte: (1) le logiciel de quadrillage "non structuré" décrit dans la demande de Gunasekera, et (2) un logiciel de quadrillage "structuré" amélioré selon la présente invention qui comporte en outre un nouveau logiciel "de quadrillage de surface structuré" selon la présente invention. Traditionnellement, les modèles géologiques consistaient en des cartes et, en fonction d'un modèle géologique donné, un modèle de simulation a été élaboré SR 15546 US/RS à partir du modèle géologique. Toutefois, dans la technique antérieure, les ingénieurs s'intéressant au réservoir modifiaient directement le modèle de simulation plutôt que de mettre à jour le modèle géologique sous-jacent. Aujourd'hui, il y a une demande croissante pour une approche meilleure et plus intégrée de la modélisation du réservoir. Le logiciel de quadrillage "structuré" amélioré selon la présente invention aide le procédé itératif consistant à modifier les modèles géologiques sous-jacents puis à propager les modifications apportées au modèle de simulation plus rapidement que ce qui est possible actuellement. Plusieurs algorithmes différents ont été proposés pour accomplir automatiquement la tâche de quadrillage. Dans la pratique, il y a deux classes principales d'algorithmes de quadrillage de la technique antérieure: (1) des méthodes algébriques dans lesquelles une formule d'interpolation est utilisée pour interpoler les courbes limites; celles-ci peuvent fonctionner correctement si la région n'est pas trop distordue; et (2) des méthodes fondées sur la solution d'équations différentielles partielles. Les meilleures méthodes de cette classe résolvent une équation du type diffusion fondée sur la transformation de l'équation de Laplace dans l'espace physique en un problème non linéaire dans l'espace logique. Certaines variantes de cette approche donnent des quadrillages de bonne qualité mais peuvent rencontrer le problème posé par des cellules inversées près des courbures importantes au niveau des limites. Toutes ces méthodes requièrent que l'utilisateur fixe la répartition des points de quadrillage sur la limite de la région. SR 15546 US/RS Résumé de l'invention Par conséquent, un but principal de la présente invention est d'améliorer le premier programme "de quadrillage" de la technique antérieure et de surmonter les inconvénients du poste de travail de la technique antérieure et des programmes de quadrillage de simulation fondés sur un logiciel en proposant un programme de quadrillage de simulation amélioré appelé "Flogrid". Un aspect principal de la présente invention est de proposer le programme de quadrillage de simulation "Flogrid" amélioré, le programme de quadrillage de simulation "Flogrid amélioré comportant un dispositif de mise à une plus grande échelle et un nouveau "dispositif de quadrillage structuré" selon la présente invention, le nouveau "dispositif de quadrillage structuré" comportant en outre un dispositif de quadrillage de bloc et un nouveau "dispositif de quadrillage de surface structuré" appelé également "dispositif de quadrillage de point de coin" selon la présente invention, le "dispositif de quadrillage de surface structure" amélioré couplé au dispositif de quadrillage de bloc quadrillant une formation terrestre en plaçant une multitude de blocs de quadrillages d'un quadrillage structuré à l'intérieur d'une formation terrestre et en générant en outre une multitude d'informations de propriétés de blocs de quadrillages structurés plus précises qui correspond, respectivement, à la multitude de blocs de quadrillages, ces informations de propriété de blocs de quadrillages structurés étant utilisées par un poste de travail et un simulateur de réservoir fondé sur un SR 15546 US/RS logiciel pour générer une multitude de "résultats de simulation" plus précis qui correspondent, respectivement, à la multitude de blocs de quadrillages du quadrillage structuré, les résultats de simulation étant conçus pour être visualisés sur une visionneuse 3D. Un autre aspect de la présente invention est de proposer le programme de quadrillage de simulation "Flogrid" amélioré, dans lequel le "dispositif de quadrillage structuré" Flogrid comporte le "dispositif de quadrillage de surface structuré" et le "dispositif de quadrillage de surface structure" sert à quadriller un horizon faillé dans la formation terrestre qui comporte une ou plusieurs lignes d'intersection de failles en: (a) renfermant les lignes d'intersection de failles sur l'horizon avec une limite, (b) en élaborant une triangulation et recouvrant la triangulation sur la partie supérieure dudit horizon de façon à renfermer la limite et les lignes d'intersection de failles sur l'horizon, la triangulation ayant des noeuds, (c) en établissant une paire de champs de vecteur ayant une pluralité de vecteurs directionnels sur les noeuds de la triangulation, ce qui permet de produire une "triangulation vectorielle", en attribuant l'un des vecteurs directionnels des champs de vecteur à chaque noeud de la triangulation, chacun des vecteurs directionnels ayant un sens, (d) en réponse aux sens de la pluralité de vecteurs directionnels sur les noeuds de la triangulation, en élaborant un réseau "lignes de contrôle" qui rencontrent les côtés opposés de la limite en établissant une interconnexion entre les extrémités de chacune des lignes d'intersection de SR 15546 US/RS failles à l'intérieur de la limite et les côtés opposés de la limite, ce qui produit un "quadrillage de surface" qui recouvre la limite, (e) facultativement, en ajoutant des "lignes supplémentaires" au quadrillage de surface en établissant une interconnexion entre les lignes supplémentaires et les côtés opposés de la limite, ce qui produit un "quadrillage de surface structuré" ayant une pluralité de noeuds qui recouvre la limite de l'horizon faillé dans la formation terrestre, et (f) soit en plaçant à des intervalles réguliers les "lignes de contrôle" et les "lignes supplémentaires" du "quadrillage de surface structuré", soit en répartissant uniformément les "lignes de contrôle" et les "lignes supplémentaires" du quadrillage de surface structuré. Un autre aspect de la présente invention est de proposer le programme de quadrillage de simulation Flogrid amélioré dans lequel le "dispositif de quadrillage structuré" Flogrid comporte en outre le "dispositif de quadrillage de bloc", le "dispositif de quadrillage de bloc" servant à établir un "quadrillage structuré" dans la formation terrestre en: (g) connectant une pluralité de lignes de coordonnées à la pluralité de noeuds du "quadrillage de surface structuré" de l'étape mentionnée ci-dessus (e), et (h) réduisant la pluralité de lignes de coordonnées à partir des noeuds du "dispositif de quadrillage de surface structuré" et en permettant aux lignes de coordonnées réduites de former une intersection avec un ou plusieurs autres "horizons" orientés plus bas dans chaque formation terrestre, ce qui permet de produire un "quadrillage structuré", le "quadrillage structuré" comportant une multitude de "blocs de quadrillages", SR 15546 US/RS chacun des blocs de quadrillages étant placé entre des horizons adjacents de la formation terrestre. Un autre aspect de la présente invention est de proposer le programme de quadrillage de simulation Flogrid amélioré dans lequel le "dispositif de mise à une plus grande échelle" Flogrid reçoit le "quadrillage structuré" du "dispositif de quadrillage structuré" et attribue les "informations de propriété de bloc de quadrillage structuré plus précises" à chaque bloc de quadrillage de la multitude de blocs de quadrillages du "quadrillage structuré", un poste de travail et un simulateur de réservoir fondé sur un logiciel recevant les "informations de propriété de bloc de quadrillage structure" associées à la multitude de blocs de quadrillages du quadrillage structuré et générant une "multitude de résultats de simulation plus précis" également associée à la multitude de blocs de quadrillages, la multitude de "résultats de simulation" et la multitude correspondante de blocs de quadrillages étant conçue pour être visualisée sur une visionneuse 3D du poste de travail. Les buts mentionnés de la présente invention ci- dessus ainsi que d'autres sont atteints et accomplis grâce à un nouvel appareil et un nouveau procédé logiciel de quadrillage de simulation, appelé ci-après "Flogrid", qui comporte un "dispositif de quadrillage structure" amélioré selon la présente invention. Le "dispositif de quadrillage structuré amélioré" comporte en outre un "dispositif de quadrillage de bloc" et un "dispositif de quadrillage de surface structuré" également appelé "dispositif de quadrillage de point de coin" selon la présente invention. SR 15546 US/RS L'appareil et le procédé logiciel de quadrillage de simulation dénommé Flogrid, qui comporte le dispositif de quadrillage "Petragrid non structure" et le logiciel "dispositif de quadrillage structuré" de la présente invention (qui comporte en outre le "dispositif de quadrillage de surface structuré" de la présente invention), est initialement stocké dans un moyen de stockage comme un CD-Rom. Le moyen de stockage est inséré dans le poste de travail d'un ordinateur, et l'appareil et le procédé logiciel de quadrillage de simulation "Flogrid" de la présente invention est ensuite chargé et stocké dans la mémoire du poste de travail. Lorsque le logiciel de quadrillage de simulation "Flogrid" est exécuté par un processeur de ce poste de travail, les données entrées sont chargées dans un dispositif de stockage de données de réservoir dans Flogrid. Les résultats générés par le dispositif de stockage de données de réservoir sont fournis au canevas du réservoir de Flogrid. Le canevas du réservoir est composé d'un modèle structurel et d'un modèle de propriété. Le modèle structurel construit des modèles structuraux en 3D et le modèle de propriété crée des propriétés d'édition. Un dispositif de mise à une plus grande échelle Flogrid permet à l'utilisateur de regrouper des couches à l'intérieur d'un modèle à échelle précise en couches grossières. Les résultats générés par le modèle de propriété sont entrés dans le dispositif de quadrillage "Petragrid non structuré" Flogrid, qui est décrit en détail dans la "demande de Gunasekera". Les résultats générés par le modèle de propriété sont également entrés dans le dispositif de SR 15546 US/RS quadrillage "structuré" Flogrid de la présente invention. Le "dispositif de quadrillage structuré" de la présente invention quadrille une section d'une formation terrestre en élaborant un quadrillage "structuré" et en plaçant ce quadrillage "structuré" sur la formation terrestre. Un simulateur reçoit ce quadrillage "structuré" et, en réponse à cela, génère un ensemble de "résultats de simulation" plus précis qui sont visualisés par un opérateur sur une visionneuse 3D. Rappelons que la formation terrestre comporte une pluralité d'horizons faillés ayant des lignes d'intersection de failles qui les traversent et supposons qu'une limite sur l'un des horizons faillés renferme les lignes d'intersection de failles sur cet horizon, le "dispositif de quadrillage de surface structuré" (également appelé "dispositif de quadrillage de point de coin") placé à l'intérieur du "dispositif de quadrillage structuré" Flogrid élabore une triangulation sur cet horizon, la triangulation refermant la limite et les lignes d'intersection de failles à l'intérieur de la limite. La triangulation comporte une pluralité de noeuds. Le "dispositif de quadrillage de surface structuré" élabore un champ de vecteur composé d'une pluralité de vecteurs directionnels sur les noeuds de la triangulation, un vecteur directionnel étant attribué à chaque noeud de la triangulation. Le "dispositif de quadrillage de surface structuré" utilise alors le champ de vecteur de la triangulation pour élaborer une pluralité de "lignes de contrôle" qui traverseront les limites, de haut en bas et de gauche à droite. Une "ligne de contrôle" est élaborée et traverse la limite lorsqu'une paire de SR 15546 US/RS "lignes d'extension" contenues dans la limite (qui renferme les lignes d'intersection de failles) forme une interconnexion entre les deux extrémités d'une ligne d'intersection de failles sur cet horizon et les côtés opposés de la limite elle-même. Le "dispositif de quadrillage de surface structuré" ajoute également quelques "lignes supplémentaires" qui traversent également les limites de haut en bas et de gauche à droite. Le "dispositif de quadrillage de surface structuré" doit donc: (1) se trouver à une distance égale des lignes de contrôle et des lignes supplémentaires à l'intérieur de la limite, ou (2) répartir uniformément les lignes de contrôle et les lignes supplémentaires à l'intérieur de la limite. A ce stade, un "quadrillage de surface structuré" recouvre la limite sur cet horizon, en établissant une interconnexion entre les côtés opposés de la limite. Le "quadrillage de surface structuré" comporte une pluralité de "noeuds". Le "dispositif de quadrillage de bloc" du "dispositif de quadrillage structuré" sert à présent à relier les "lignes de coordonnées" aux "noeuds" du "quadrillage de surface structuré" et le "dispositif de quadrillage de bloc" réduit ces lignes de coordonnées à partir des noeuds du "quadrillage de surface structuré", dans la formation terrestre, ce qui permet aux lignes de coordonnées de former une intersection entre les horizons situés le plus bas à l'intérieur de la formation terrestre, ce qui produit une multitude de "blocs de quadrillage" d'un "quadrillage structuré". Les blocs de quadrillage du "quadrillage structuré" sont chacun situés entre les horizons adjacents de la formation terrestre, un "bloc SR 15546 US/RS de quadrillage unique" étant défini par quatre lignes de coordonnées adjacentes et une paire d'horizons adjacents de la formation terrestre limités au niveau des côtés opposés du bloc de quadrillage. Rappelons que le programme de quadrillage de simulation Flogrid comporte un dispositif de mise à une plus grande échelle, le dispositif de mise à une plus grande échelle attribuant à présent une propriété de blocs de quadrillage appropriée à chaque bloc de quadrillage du "quadrillage structuré". Lorsque la pluralité de propriété de blocs de quadrillages est attribuée à la pluralité respective des blocs de quadrillage du "quadrillage structuré" par le dispositif de mise à une plus grande échelle Flogrid, un "quadrillage structuré" se produit ainsi que des propriétés associées à celui- ci. Le "quadrillage structuré et les propriétés associées à celui-ci", qui représentent un "ensemble d'informations de propriétés de blocs de quadrillages plus précises", passent dans un simulateur de réservoir. En conséquence, rappelons que le procédé "Flogrid" de la présente invention comporte le dispositif de quadrillage "Petragrid non structuré" et le "dispositif de quadrillage structuré" de la présente invention, et le dispositif de quadrillage "Petragrid non structuré" ou le dispositif de quadrillage "structuré" de la présente invention produit, lorsqu'il est exécuté "l'ensemble d'informations de propriétés de blocs de quadrillages plus précises". Ces informations de propriétés de blocs de quadrillages plus précises sont reçues par le simulateur de réservoir qui est également stocké dans la mémoire du poste de travail. Le logiciel simulateur Eclipse répond aux informations SR 15546 US/RS de propriétés de blocs de quadrillages plus précises provenant de "Flogrid" en produisant un ensemble de résultats de simulation plus précis (par exemple des valeurs de pression et de saturation). Chacun des "résultats de simulation" plus précis est associé à l'un des blocs de quadrillage du quadrillage structuré qui est placé sur la formation terrestre et ces "résultats de simulation" sont affichés sur une visionneuse 3D du poste de travail pour être observés et analysés par un opérateur du poste de travail. Une plus grande étendue des capacités d'application de la présente invention apparaîtra plus clairement à la lecture de la description détaillée présentée ci-après. Toutefois, il faut comprendre que la description détaillée et les exemples spécifiques, bien qu'ils représentent un mode de réalisation préféré de la présente invention, ne sont donnés qu'à titre d'illustration, puisque plusieurs changements et modifications restant dans l'esprit et la portée de la présente invention apparaîtront évidents au spécialiste de la technique à la lecture de la description détaillée suivante. Brève description des dessins Une compréhension complète de la présente invention sera obtenue à la lecture de la description détaillée du mode de réalisation préféré présenté ci- après, et des dessins joints, qui sont donnés à titre d'illustration uniquement et qui ne limitent en aucun cas la présente invention, et parmi lesquels: les figures 1 à 23 sont utilisées conjointement avec la "description du mode de réalisation préféré", parmi lesquelles: SR 15546 US/RS les figures 1 à 13d2 sont des dessins de la "demande de Gunasekera" qui décrit le logiciel de quadrillage non structuré "Petragrid" (décrit dans la demande de brevet américain en suspens antérieure numéro de série 08/873.234 déposée le 11 juin 1997 et intitulée "Method and Apparatus for generating more accurate earth formation grid cell property information for use by a simulator to display more accurate simulation results of the formation near a wellbore" attribuée à Dayal Gunasekera, lesdits dessins étant fournis ici pour établir l'arrière-plan du logiciel de quadrillage "structuré" de la présente invention; les figures 14a à 14c illustrent le poste de travail de la figure 9 comportant la mémoire du poste de travail 40c, dans lequel la mémoire de poste de travail 40c est conçue pour recevoir et stocker un ensemble logiciel stocké sur un CD-Rom, le logiciel comportant l'appareil et le procédé de quadrillage de simulation "Flogrid" qui comportent le dispositif de quadrillage "structuré" de la présente invention, ainsi que le logiciel de quadrillage "non structuré" Petragrid, la figure 14c illustrant le nouveau contenu de la mémoire du poste de travail 40c une fois que le programme de quadrillage de simulation "Flogrid" stocké dans un CD-Rom est chargé dans le poste de travail, le nouveau contenu de la mémoire 40c comportant à la fois le logiciel de quadrillage "non structuré" Petragrid cl décrit dans la "demande de Gunasekera" et le logiciel de quadrillage "structuré" de la présente invention; la figure 15 illustre l'opération fonctionnelle effectuée par l'appareil et le procédé de quadrillage de simulation "Flogrid" et le logiciel simulateur SR 15546 US/RS "Eclipse" lorsque l'appareil et le procédé de quadrillage de simulation "Flogrid" de la présente invention et le logiciel simulateur "Eclipse" sont exécutés par le processeur du poste de travail en réponse à un enregistrement de sortie de diagraphie et un enregistrement de sortie de données sismiques réduit; La figure 16 illustre une construction plus détaillée de l'appareil et du procédé de quadrillage de simulation "Flogrid" de la figure 14c qui comporte à la fois le dispositif de quadrillage "non structuré" Petragrid qui est conçu pour générer les informations de propriétés de blocs de quadrillages de formation terrestre plus précises destinées à être utilisées par le simulateur de réservoir "Eclipse"; les figures 16a à 16c illustrent un quadrillage de surface structuré ayant une pluralité de noeuds, sur la figure 16 a, un ensemble de lignes de coordonnées réduites à partir des noeuds du quadrillage de surface structuré, sur la figure 16b, et la formation d'un "bloc de quadrillage" en coupant l'ensemble de lignes de coordonnées par une paire d'horizons, sur la figure 16c; les figures 17 à 23 sont utilisées conjointement avec une description fonctionnelle du fonctionnement du logiciel de quadrillage "structuré" de la présente invention conçu pour être utilisé à l'intérieur de l'appareil et du procédé de quadrillage de simulation Flogrid, parmi lesquels: la figure 17 illustre sept (7) étapes qui sont exécutées par le logiciel de quadrillage "structuré" Flogrid 68 de la présente invention au cours de la construction d'un "quadrillage structuré", en SR 15546 US/RS commençant par la construction d'un "quadrillage de surface structure"; la figure 18 illustre une limite, comportant deux failles i et une faille j unique; la figure 19 illustre des extensions uniformes de failles allant jusqu'à la limite pour élaborer des lignes de contrôle; la figure 20 illustre des lignes de contrôle ayant certaines lignes supplémentaires non uniformes; la figure 21 illustre des lignes de contrôle ayant certaines lignes supplémentaires - uniformes; la figure 22 illustre un exemple de quadrillage isotrope entre deux failles; la figure 23 illustre un exemple de quadrillage anisotrope entre deux failles; les figures 24 à 52 sont utilisées conjointement avec "la description détaillée du mode de réalisation préféré", parmi lesquelles: la figure 24 illustre une autre construction détaillée de l'appareil et du procédé de quadrillage de simulation "Flogrid" 60a de la présente invention comportant le "dispositif de quadrillage structuré" 68 de la présente invention; la figure 24a illustre une formation terrestre faillée comportant une pluralité d'horizons et une source de vibrations sonores qui sont réfléchies à l'extérieur des horizons afin de générer un enregistrement de sortie de données sismiques; la figure 25 illustre une construction du canevas structurel 149 de la figure 24; la figure 25a illustre un "modèle d'horizon faillé final" 151 (c'est-à-dire un modèle géologique SR 15546 US/RS 151) qui est généré par le modèle structurel 66a de la figure 25; la figure 26 illustre une construction du dispositif de quadrillage structuré 68 de la figure 24; la figure 27 illustre les étapes fonctionnelles exécutées par le "quadrillage de surface structuré" 68b de la figure 26; et les figures 28 à 52 sont des dessins qui sont utilisés conjointement avec une description détaillée de l'appareil et du procédé de quadrillage de simulation Flogrid 60a comportant une description détaillée du "dispositif de quadrillage de surface structuré" 68b de la figure 26 de la présente invention traitée dans la "description détaillée du mode de réalisation préféré"; parmi lesquelles: les figures 28a et 28b représentent les "données d'entrée" qui sont entrées dans Flogrid 60a; les figures 29 et 30 illustrent des horizons dans une formation terrestre faillée et la définition de "cartes de propriété"; la figure 31 illustre un horizon faillé 151 ayant des lignes d'intersection de failles 134 qui le traversent, et une "limite" 136 qui renferme les lignes d'intersection de failles 134. la figure 32 illustre la limite 136 de la figure 31 avant la distorsion; les figures 33 et 34 illustrent la limite 136 de la figure 31 après la distorsion et présentant un ensemble de lignes d'intersections de failles renfermées par la limite; la figure 35 illustre la limite 136 de la figure 34 après distorsion et ayant une pluralité de "lignes SR 15546 US/RS de contrôle" également appelées "lignes de flux" et "lignes supplémentaires" 150 placées sur la limite 136, ce qui permet de créer un "quadrillage de surface structuré"; les figures 35a et 36 illustrent la façon dont les "lignes de coordonnées" des figures 16a à 16c sont réduites à partir des noeuds du "quadrillage de surface structuré" pour créer une multitude de "blocs de quadrillages" du "quadrillage structuré"; les figures 37 et 38 illustrent la façon dont le dispositif de quadrillage structuré 68 gère les failles en zigzag 144; les figures 39 et 40 illustrent la fonction du "dispositif de mise à une plus grande échelle" Flogrid; les figures 41 et 42 illustrent des étiquettes de point sur le cube de référence et le bloc de quadrillage (figure 41) et l'interpolation d'un quadrillage précis 2xlx3 dans un bloc de quadrillage (figure 42); la figure 43 illustre la façon dont un utilisateur détermine une limite et choisit les failles i et j sur lesquelles appliquer un quadrillage; la figure 44 illustre la construction d'un quadrillage triangulaire régulier (une triangulation) sur la limite de la figure 43, la limite de la triangulation se trouvant bien à l'extérieur de la limite spécifiée par l'utilisateur; la figure 45a illustre les vecteurs "i et "j" appelés "vecteurs directionnels"; la figure 45bl illustre une triangulation 178a, 178b recouvrant une limite spécifiée par l'utilisateur SR 15546 US/RS 176 qui contient et renferme une pluralité de lignes d'intersection de failles 170, 172, 174; les figures 45b2 et 45b3 illustrent la façon dont le type "gelé" de "vecteurs directionnels" est attribué aux noeuds de la triangulation de la figure 45bl; les figures 46 et 47 illustrent les techniques d'interpolation anisotrope et isotrope qui sont des techniques permettant d'attribuer le type "non gelé" de "vecteurs directionnels" aux noeuds de la triangulation de la figure 45bl; la figure 48 illustre une limite spécifiée par l'utilisateur 176 renfermant les lignes d'intersection de failles 170, 172, 174 et ayant des "lignes de commande" également appelées "lignes de flux" allant du haut au bas et de la gauche à la droite de la limite 176, une telle "ligne de contrôle" étant identifiée par les numéros d'éléments suivants [200, 170, 200] et une autre "ligne de contrôle" étant identifiée par les numéros d'éléments suivants [200, 172, 200], et une autre "ligne de contrôle" étant identifiée par les numéros d'éléments suivants: [200, 174, 200]; la figure 49 illustre la façon dont une "ligne de contrôle" est élaborée en utilisant les vecteurs directionnels d'un champ de vecteur qui sont attribués à un noeud de la triangulation de la figure 45bl; la figure 50 utilise une limite spécifiée par l'utilisateur 176, les lignes d'intersection de failles , 172, 174, les lignes de contrôle de la figure 48 et un ensemble de "lignes supplémentaires" 240 qui s'ajoutent aux lignes de contrôle, ce qui permet de créer un "quadrillage de surface structuré" qui recouvre le limite 176 sur un horizon faillé dans une formation terrestre; SR 15546 US/RS la figure 51 illustre la façon dont l'utilisateur détermine le nombre de "lignes supplémentaires" qu'il souhaite avoir entre les "lignes de commande"; et la figure 52 illustre la façon dont les largeurs entre les "lignes de contrôle" et les "lignes supplémentaires" du "quadrillage de surface structuré" de la figure 50 peuvent varier progressivement, c'est- à-dire que les largeurs peuvent varier de façon "uniforme"; sur la figure 52, les largeurs ne sont pas espacées d'une distance égale, tel qu'illustré sur la figure 50. Description du mode de réalisation préféré La description du mode de réalisation préféré sera divisée en deux sections, une première section et une deuxième section. La première section de la description du mode de réalisation préféré de la demande, en référence aux figures 1 à 13d2, reprend certains aspects de la demande de Gunasekera identifiée cidessus, qui est la demande en suspens de la technique antérieure qui décrit le logiciel de quadrillage "non structuré" Petragrid qui est placé à l'intérieur du programme de quadrillage de simulation Flogrid. La deuxième section de la description du mode de réalisation préféré de la présente demande, en référence aux figures 14a à 23, traite du programme de quadrillage de simulation "Flogrid" comportant le logiciel de quadrillage "structuré", qui comporte en outre le "dispositif de quadrillage de surface structuré" de la présente invention qui est placé à SR 15546 US/RS l'intérieur du programme de quadrillage de simulation Flogrid. La description détaillée du mode de réalisation préféré, en référence aux figures 24 à 52, traitera en détail du programme de quadrillage de simulation "Flogrid" comportant un traitement détaillé du "dispositif de quadrillage de surface structuré" de la présente invention. Première section - Le logiciel de quadrillage "non structuré" Petragrid utilisé dans Flogrid (1) Informations d'arrière-plan En référence aux figures 1 à 13d2, en référence initiale aux figures 1 à 7, 11 et 12, les opérations de diagraphie pour produire des données de diagraphie et des opérations sismiques permettant de produire des données sismiques sont illustrées. Sur la figure 1, une pluralité d'horizons de formation terrestre 13 comprend une formation terrestre 15 et une paire de puits de forage 17 a été forée dans la formation 15 afin d'extraire finalement du fluide du puits de forage de la formation. La formation terrestre sera modélisée et simulée sur un ordinateur de poste de travail d'interprétation par le logiciel "Petragrid". Sur la figure 2, l'un des puits de forage 17 de la figure 1 est illustré et pénètre dans la formation terrestre 15 de la figure 1. Un camion chargé de la diagraphie 19 fait descendre un instrument de diagraphie 21 dans le puits de forage 17 et l'instrument de diagraphie 21 stimule et provoque la formation 15. En réponse à cela, des détecteurs de l'instrument de diagraphie 21 reçoivent des signaux de SR 15546 US/RS la formation 15 et, en réponse à cela, d'autres signaux représentant les données de diagraphie 27 propagent vers le haut du trou, entre l'instrument de diagraphie 21 et l'ordinateur du camion chargé de la diagraphie 23. Un enregistrement de sortie de diagraphie 25 est généré par l'ordinateur du camion chargé de la diagraphie 23 qui affiche les données de diagraphie 27. Sur la figure 3, une construction plus détaillée de l'ordinateur du camion chargé de la diagraphie 23 est illustrée. Un bus 23a reçoit les données de diagraphie 27 et, en réponse à cela, l'enregistrement de sortie de diagraphie 25 est généré par le processeur 23b, l'enregistrement de sortie de diagraphie 25 affichant et/ou stockant les données de diagraphie 27. L'enregistrement de sortie de diagraphie 25 est entré dans le poste de travail d'interprétation des figures 8 et 9. Sur les figures 4 et 5, un appareil et un procédé associé pour exécuter une opération sismique tridimensionnelle (3D) en un emplacement sur la surface de terre proche des puits de forage 17 de la figure 1 est illustré. Sur la figure 4, une source d'énergie acoustique ou explosive 10 située en dessous de la surface terrestre 12 détonne et génère une pluralité de vibrations sonores ou acoustiques 14 qui se propagent vers le bas et se réfléchissent hors d'une couche de l'horizon 16 à l'intérieur de la formation terrestre. La couche d'horizon 16 peut être une couche supérieure de roche ou de sable ou de schiste. Lorsque les vibrations sonores se réfléchissent hors de la couche de l'horizon 16, les vibrations sonores 14 se propagent vers le haut et sont reçues par une pluralité de SR 15546 US/RS récepteurs 18 appelés géophones 18 situés à la surface de la terre. La pluralité de géophones 18 génère chacun un signal électrique en réponse à la réception d'une vibration sonore et une pluralité de signaux électriques sont générés par les géophones 18, la pluralité de signaux étant reçue dans un camion d'enregistrement 20. La pluralité de signaux électriques provenant des géophones 18 représente un ensemble de caractéristiques de la formation terrestre comportant les horizons 16 situés à l'intérieur de la terre en dessous des géophones 18. Le camion d'enregistrement 20 contient un ordinateur 20a qui reçoit et stocke la pluralité de signaux reçus des géophones 18. Un enregistrement de sortie sismique 20a4 est généré à partir de l'ordinateur 20a dans le camion d'enregistrement 20 qui comporte et/ou affiche et/ou stocke la pluralité de signaux électriques qui sont représentatifs des caractéristiques de la formation terrestre comportant les horizons 16 situés à l'intérieur de la terre en dessous des géophones 18. Sur la figure 5, un autre procédé et appareil permettant d'exécuter une opération sismique en 3D est illustré. La figure 5 a été prise dans un livre intitulé "Seismic Velocity Analysis and the Convolutional Model" de Enders A. Robinson. L'opération sismique en 3D de la figure 4 est effectuée 10 fois. Par exemple, lorsque la source d'énergie explosive 10 est située dans la position 22 (première position ou position "0" le long de la surface de la terre) sur la figure 4, une première pluralité de signaux électriques provenant des géophones 18 est stockée dans l'ordinateur 20a dans le camion d'enregistrement 20. La source d'énergie explosive est déplacée jusqu'à la SR 15546 US/RS position 24. Lorsque la source d'énergie explosive 10 est située à la position 24, (deuxième position ou position "1" le long de la surface de la terre), une deuxième pluralité de signaux électriques est stockée dans l'ordinateur 20a dans le camion d'enregistrement 20. La source d'énergie explosive 10 est déplacée de façon répétitive et séquentielle des positions "2" à "9" sur la figure 5 jusqu'à ce qu'elle soit située en position 26 (c'est-à-dire la position "9" qui est la dixième position) sur la surface de la terre. Lorsque la source d'énergie explosive 10 est située en position 26 (dixième position le long de la surface de la terre), une dixième pluralité de signaux électriques est stockée dans l'ordinateur 20a dans le camion d'enregistrement 20. En conséquence, sur la figure 5, le camion d'enregistrement 20 enregistre 10 traces (dix ensembles de signaux électriques, chaque ensemble étant une pluralité de signaux électriques) entre la position 22 et la position 26 le long de la surface de la terre. Un enregistrement de sortie sismique 20a4 sera généré par l'ordinateur 20a dans le camion d'enregistrement 20 qui comporte les dix traces ou dix ensembles de signaux électriques reçues des géophones 18. Sur la figure 6, une construction plus détaillée de l'ordinateur du camion d'enregistrement 20a est illustrée. L'ordinateur du camion d'enregistrement 20a de la figure 4 comporte un processeur 20al et une mémoire 20a2 connectée à un bus de système. Les dix traces ou dix ensembles de signaux électriques (reçus des géophones 18 au cours de l'opération sismique 3D) sont reçus dans l'ordinateur du camion d'enregistrement a par l'intermédiaire du bloc "de données sismiques reçues" 20a3 de la figure 3 et sont stockés dans la SR 15546 US/RS mémoire 20a2 de l'ordinateur du camion d'enregistrement a. Si on le souhaite, un enregistrement de sortie sismique 20a4 est généré par l'ordinateur du camion d'enregistrement 20a, l'enregistrement de sortie sismique 20a4 étant conçu pour stocker et afficher "une pluralité de données sismiques" représentant les dix traces ou les dix ensembles de signaux électriques reçus par l'ordinateur du camion d'enregistrement 20a des géophones 18. Sur la figure 7, un schéma simplifié d'un ordinateur à unité centrale 30 est illustré et celui-ci utilise un "logiciel de réduction de données" stockées pour effectuer une opération "de réduction de données" sur la "pluralité de données sismiques" incluses dans l'enregistrement de sortie sismique 20a4 de la figure 6. L'ordinateur à unité centrale 30 produit un "enregistrement de sortie de données sismiques réduites" 30d de la figure 7 qui est conçu pour stocker et afficher des informations qui représentent des versions "réduites" de la "pluralité de données sismiques" incluses dans le moyen d'enregistrement de sortie 20a4 de la figure 3. L'ordinateur à unité centrale 30 de la figure 4 comporte un processeur d'unité centrale 30a connecté à un bus de système et à une mémoire 30b également connectée au bus du système qui stocke dans celui-ci un "logiciel de réduction de données". L'enregistrement de sortie sismique 20a4 de la figure 6, qui comporte la "pluralité de données sismiques", est connecté au bus du système de l'ordinateur à unité centrale 30 de la figure 7. En conséquence, la "pluralité de données sismiques" incluses dans le moyen d'enregistrement de sortie 20a4 de la figure 6, est maintenant entrée dans le SR 15546 US/RS processeur d'unité centrale 30a de la figure 7. Le processeur 30a de l'ordinateur à unité centrale 30 de la figure 7 exécute le "logiciel de réduction de données" stockées dans la mémoire 30b de l'ordinateur à unité centrale. Le "logiciel de réduction de données", qui est stocké dans la mémoire 30b de l'ordinateur à unité centrale 30 de la figure 4, peut être trouvé dans un livre intitulé "Seismic Velocity Analysis and the Convolutional Model", de Enders A. Robinson, dont la description est incorporée en tant que référence dans la présente demande. Lorsque le "logiciel de réduction de données" stockées dans la mémoire 30b est exécuté, le processeur de l'unité centrale 30a effectue une opération "de réduction de données" sur la "pluralité de données sismiques" qui est incluse dans l'enregistrement de sortie sismique 20a4 de la figure 7. Lorsque "l'opération de réduction de données" est effectuée, le processeur de l'unité centrale 30a génère un "enregistrement de sortie de données sismiques réduit" 30d qui stocke et est conçu pour afficher des informations: représentant une "version réduite" de la "pluralité de données sismiques" incluse dans l'enregistrement de sortie sismique 20a4 de la figure 7 et comportant un ensemble de caractéristiques appartenant à la formation terrestre située près du puits de forage 17 de la figure 1, les caractéristiques comportant l'emplacement et la structure des horizons 16 de la figure 4. Sur les figures 11 et 12, un organigramme du logiciel de réduction de données 30b stocké dans la mémoire 30b de l'ordinateur à unité centrale 30 de la figure 7 est illustré. L'organigramme du logiciel de réduction de données des figures 11 et 12 a été pris SR 15546 US/RS dans un livre intitulé "Seismic Velocity Analysis and the Convolutional Model", de Enders A. Robinson. L'organigramme du logiciel de réduction de données 30b comporte les blocs suivants: un bloc de démultiplexage 30bl relié à l'entrée, un bloc de tri 30b2, un bloc de retrait de gain 30b3, un bloc de filtrage de fréquence b4, un bloc de rééchantillonnage 30b5, un bloc de sélection de traces 30b6, une sortie 30b7 appelée regroupements choisis" (analyses de vitesse), correction d'amplitude 30b8, déconvolution 30b9, une deuxième sortie 30b10 appelée "traces triées CMP après déconvolution", un bloc de correction de temps 30bll, un bloc de commande automatique de gain 30b12, un bloc d'empilage 30b13, une troisième sortie 30b14 appelée "traces empilées (non filtrées)", un bloc de filtrage de fréquence 30b15, un autre bloc de commande automatique de gain 30b16, une quatrième sortie 30b17 appelée "traces empilées (filtrées)", une deuxième entrée appelée "informations en profondeur" 30b18, un bloc d'interpolation de traces 30b19, un bloc de migration 30b20, une cinquième sortie 30b21 appelée "traces migrées (non filtrées)", un bloc de filtrage de fréquence 30b22, un bloc de commande automatique de gain 30b23, une sixième sortie 30b24 appelée "traces migrées (filtrées)", un bloc de correction temps sur profondeur 30b25 et une septième sortie 30b26 appelée "traces migrées (migration en profondeur)". Dans l'organigramme des figures 11 et 12, l'une quelconque des sorties 30b7, 30b10, 30b14, 30b17, 30b21, 30b24 et 30b26 peut être utilisée comme entrée de l'ordinateur à unité centrale 30 ou du poste de travail SR 15546 US/RS d'interprétation mentionné ci-dessus en référence aux figures 8 et 9 des dessins. (2) Résumé de l'appareil et du procédé de quadrillage non structuré Petragrid En référence aux figures 8, 9, 10, 13a, 13b, 13cl et 13c2, un résumé de l'appareil et du procédé Petragrid, qui réagit à l'enregistrement de sortie de diagraphie 25 et l'enregistrement de sortie de données sismiques réduites 30d, produit par les opérations de diagraphie et des opérations sismiques des figures 1 à 7 est illustré. Sur la figure 8, l'enregistrement de sortie de diagraphie 25 de la figure 3 et l'enregistrement de sortie de données sismiques réduites 30d de la figure 7 sont tous deux entrés dans un poste de travail d'interprétation 40. Ce poste de travail 40 stocke le logiciel Petragrid de la présente invention qui, lorsqu'il est exécuté, génère des informations de propriété de cellule de quadrillage d'horizon plus précises qui sont conçues pour être utilisées par un programme de simulation, ledit programme de simulation, lorsqu'il est exécuté, modélisant et simulant les propriétés de production du fluide de forage de la formation terrestre 15 adjacente à un ou plusieurs puits de forage 17 de la figure 1. Sur la figure 9, le poste de travail d'interprétation 40 de la figure 8 est illustré plus en détail. Le poste de travail 40 comporte un bus 40a, un processeur de poste de travail 40b connecté au bus 40a, une mémoire de poste de travail 40c connectée au bus a et un affichage 40d connecté au bus 40a. L'enregistrement de sortie de diagraphie 35 et SR 15546 US/RS l'enregistrement de sortie de données sismiques réduites 30d sont tous deux connectés au bus 40a du poste de travail 40; en conséquence, les données de diagraphie et les données sismiques réduites stockées dans l'enregistrement de sortie de diagraphie 25 et l'enregistrement de sortie de données sismiques réduites, respectivement, sont entrés dans le poste de travail 40 et fournis au processeur du poste de travail b. Le contenu de la mémoire du poste de travail 40c est illustré sur la figure 10. Sur la figure 10, la mémoire de poste de travail c de la figure 9 stocke au moins deux blocs de logiciel: (1) un premier bloc de logiciel appelé ci-après le "logiciel Petragrid" 40cl qui génère des informations de propriété de cellule de quadrillage d'horizon plus précises nécessaires au logiciel simulateur Eclipse 40c2 traité ci- dessous; et (2) un deuxième bloc de logiciel appelé ci-après le "logiciel simulateur Eclipse" 40c2 qui reçoit les informations de propriété de cellule de quadrillage d'horizon plus précises qui sont générées par le logiciel Petragrid 40cl; le logiciel simulateur Eclipse 40c2 est traité dans un manuel intitulé "ECLIPSE 100 User Manual", que l'on peut se procurer auprès de "GeoQuest, division de Schlumberger Technology Corporation" située à Abingdon, Royaume-Uni. Le "ECLIPSE 100 User Manual". Le logiciel simulateur Eclipse 40c2 modélise et simule plus précisément les propriétés de production du fluide du puits de la formation terrestre 15 située près des puits de forage 17 de la figure 1 résultant des informations de SR 15546 US/RS propriété de cellule de quadrillage d'horizon plus précis qui sont reçues du logiciel Petragrid 40cl. En référence aux figures 13a et 13b, 13cl et 13c2, les paragraphes suivants résument le logiciel Petragrid 40cl de la figure 10 de la présente invention. Sur la figure 13a, la formation terrestre 15 de la figure 1 est de nouveau illustrée, la formation 15 comportant quatre (4) horizons 13 qui traversent l'étendue longitudinale de la formation 15 sur la figure 13a. Rappelons qu'un "horizon" 13 est défini comme étant la surface supérieure d'une couche de formation terrestre, la couche de formation terrestre comprenant par exemple du sable, du schiste ou du calcaire, etc. Toutefois, le logiciel Petragrid 40cl "quadrille" la formation 15 qui est située entre les couches d'horizons 13. C'est-à-direqu'entre les horizons 13 et sur l'horizon situé le plus haut 13 et en dessous de l'horizon situé le plus bas 13, le logiciel Petragrid cl "quadrille" la formation 15. Lors du quadrillage de la formation 15, la formation 15 est divisée en une multitude de cellules individuelles qui, lorsqu'elles sont reliées, comprennent le quadrillage. Sur la figure 13b, par exemple, la formation 15 comporte un horizon situé le plus haut 13a et un horizon situé le plus bas 13b qui est séparé de l'horizon 13a situé le plus haut par une couche de formation terrestre intermédiaire 15a. La couche de formation terrestre intermédiaire 15a comporte, par exemple, une couche de sable ou une couche de schiste ou une couche de calcaire, etc. Le programme de quadrillage de simulation Flogrid, comportant le SR 15546 US/RS logiciel Petragrid 40cl et le dispositif de quadrillage de surface structuré de la présente invention "quadrillent" la couche de formation terrestre l5a. C'est-à-dire que la couche de formation 15a est divisée, par le programme de quadrillage de simulation Flogrid (comportant le logiciel Petragrid 40cl), en une multitude de cellules individuelles 15al. Dans la technique antérieure, le produit logiciel connu sous le nom de "Grid" a été commercialisé par GeoQuest, division de Schlumberger Technology Corporation, Abingdon, Royaume-Uni. Le logiciel "Grid" divise les couches de formation 15a en une multitude de cellules. Toutefois, chaque cellule de la multitude de cellules était approximativement de coupe "rectangulaire". Sur la figure 13b, les cellules 15al sont représentées comme ayant une coupe approximativement "rectangulaire". Toutefois, sur la figure 13a, le logiciel Petragrid 40cl crée également la multitude de cellules al dans la formation terrestre 15 entre les horizons 13; toutefois, chaque cellule 15al peut présenter une forme en coupe transversale qui, en plus d'être approximativement "rectangulaire" (ci-après appelée quadrillage "structuré"), et soit approximativement polygonale soit en forme de coupe tétraédrique (appelé ci-après quadrillage "non structuré"). La figure 13a montre clairement une multitude de cellules 15al parmi lesquelles chaque cellule 15al présente une forme en coupe transversale qui est soit approximativement "polygonale" soit "tétraédrique" (c'est-à-dire "non structurée") en plus de "rectangulaire" (c'est-à-dire "structurée"). SR 15546 US/RS Sur la figure 13cl, après avoir décrit, en référence aux figures 13a et 13b, une fonction du logiciel Petragrid 40cl (c'est-à-dire "le quadrillage" de la formation avec une multitude de cellules) et décrit une première fonction potentiellement nouvelle du logiciel Petragrid 40cl (c'est-à-dire en créant des cellules individuelles du quadrillage qui présente une coupe transversale polygonale et tétraédrique en plus de rectangle), la figure 13cl illustre les "caractéristiques potentiellement nouvelles supplémentaires" du logiciel Petragrid 40cl de la figure 10. Sur la figure 13cl, le logiciel Petragrid 40cl reçoit l'enregistrement de sortie de diagraphie 25 de la figure 3 et l'enregistrement de sortie de données sismiques réduites 30d de la figure 7 et, en réponse à cela, le logiciel Petragrid 40cl génère un ensemble de données de sorties "plus précises" 40c3, les données de sortie "plus précises" 40c3 étant entrées dans le logiciel simulateur Eclipse 40c2. Le logiciel simulateur Eclipse 40c2 génère à son tour un ensemble "plus précis" de résultats de simulation qui sera affiché sur l'affichage du poste de travail 40d de la figure 9. En consequence, le logiciel simulateur Eclipse 40c2 modélise et simule "plus précisément" la formation terrestre 15 située près des puits de forage 17 de la figure 1 afin de déterminer les propriétés de production du fluide du puits de la formation terrestre 15. Sur la figure 13cl, le logiciel Petragrid 40cl comporte les blocs de code nouveaux suivants (blocs clA, 40clB, 40clC, 40clD et 40clE) qui représentent de nouvelles caractéristiques par rapport à celles du SR 15546 US/RS logiciel "Grid" de la technique antérieure et d'autres logiciels de la technique antérieure. Un organigramme plus complet du logiciel Petragrid 40cl sera traité ci- après en référence aux figures 13c2: 1. tétraédralisation/triangulation de base - bloc clA Sur la figure 13cl, avant de commencer, une définition du terme "perpendiculaire à K" est donnée. Un quadrillage est localement "perpendiculaire à K" si "Ka" est parallèle à "r", "k" étant le tenseur de perméabilité locale, "a" le vecteur de zone d'interface entre une paire de cellules, et "r" le vecteur allant du centre d'une cellule à l'autre. La génération d'une triangulation de base ou d'une tétraédralisation, comportant les coordonnées des sommets et du centre de chaque élément est traitée ici. Dans la technique antérieure, des quadrillages perpendiculaires à K sont générés en répartissant les points dans un espace à 2D ou 3D en fonction de la géométrie du modèle physique. Les points répartis subissent ensuite une triangulation (en 2D) ou une tétraédralisation (en 3D). Le centre de chaque triangle est calculé en faisant se croiser les lignes perpendiculaires à K et les points bissecteurs des côtés. Le centre de chaque tétraèdre est calculé en faisant se croiser les lignes perpendiculaires à K et les centres des cercles circonscrits des faces. Cette approche présente un inconvénient en ce qu'un centre se trouve loin du centre d'un élément lorsque la perméabilité présente un niveau élevé d'anisotropie, ce qui provoque les chevauchements des cellules de simulation. Selon un autre aspect du logiciel Petragrid concernant "la tétraédralisation/triangulation de base" *SR 15546 US/RS (bloc 40clA), un "nouvel algorithme" génère un quadrillage de base triangulaire "perpendiculaire à K" pour des quadrillages non structurés bidimensionnels (2D), et le nouvel algorithme génère également un quadrillage de base approximativement de forme tétraédrique "perpendiculaire à K" pour des quadrillages non structurés trimensionnels (3D). Toutefois, selon le nouvel algorithme qui génère le quadrillage de base ayant quasiment une forme de triangle "perpendiculaire à K" pour des quadrillages non structurés bidimensionnels "2D" et le quadrillage de base présentant une forme approximative de tétraèdre "perpendiculaire à K" pour des quadrillages non structurées trimensionnels "3D", effectuent cela en établissant l'échelle de la géométrie physique. Plus particulièrement, le "nouvel algorithme" mentionné ci- dessus associé à la " tétraédralisation/triangulation de base" (bloc 40cA de la figure 13cl) qui: (1) génère le quadrillage de base présentant une forme approximative de triangle "perpendiculaire à K" pour des quadrillages non structurés bidimensionnels (2D) en "établissant l'échelle" de la géométrie physique et (2) génère le quadrillage de base présentant une forme approximative de tétraèdre "perpendiculaire à K" pour des quadrillages non structurés tridimensionnels (3D) en "établissant l'échelle" de la géométrie physique, génère les deux types mentionnés ci-dessus de quadrillage en effectuant les étapes fonctionnelles suivantes: (1) on établit une carte de la géométrie physique dans un espace de calcul en utilisant une transformation fondée sur les composants du tenseur de perméabilité dans chaque couche du modèle; (2) des SR 15546 US/RS points sont répartis dans l'espace de calcul en fonction de la géométrie dessinée; (3) une triangulation de Delaunay ou une tétraédralisation des points est créée, dans lesquels: (3a) le centre de chaque élément (c'est-à-dire chaque "cellule") de la triangulation correspond au centre d'un cercle qui est tracé par les sommets de chaque cellule, et (3d) le centre de chaque cellule de la tétraédralisation correspond au centre (ou centre du cercle circonscrit) d'une sphère qui est tracée sur les sommets de chaque élément ou cellule; et (4) les sommets des éléments ou cellules et les centres des cellules sont ensuite transformés dans l'espace physique. Cette approche produit un quadrillage de meilleure qualité. 2. Agrégation triangle/tétraèdre - Bloc 40clB Sur la figure 13cl, avant de générer un quadrillage de simulation, tels que les quadrillages des figures 13a et 13b, les éléments ou cellules d'une triangulation de base (quadrillage triangulaire dans un espace bidimensionnel) ou une tétraédralisation (quadrillage de tétraèdre dans un espace tridimensionnel) peut être agrégée en des ensembles afin de réduire le nombre de cellules dans un quadrillage composite triangulaire/de tétraèdre et pour réduire le nombre de faces d'un quadrillage PEBI, avec une perte minime de précision des résultats de simulation suivants. Un nouvel algorithme associé à une "agrégation" triangle/tétraèdre (bloc 40clB de la figure 13cl) s'agrège en ensembles, une paire de cellules en forme de "triangle" ou une paire de cellules en forme de "tétraèdre" d'un quadrillage (comme les cellules 15al de la figure 13a) lorsque la SR 15546 US/RS distance entre le centre du cercle circonscrit de la paire de cellules est considéré comme petit par rapport à la taille de la paire de cellules. Cette opération est effectuée de façon itérative de sorte qu'une pluralité d'ensembles agrégés est créée, chaque ensemble de la pluralité d'ensembles agrégés comportant une ou plusieurs cellules. Si un ensemble de cellules de forme triangulaire d'un quadrillage ou bien un ensemble de cellules en forme de tétraèdre d'un quadrillage sont agglomérés, le nouvel algorithme associé à "l'agrégation triangle-tétraèdre" (bloc 40clB de la figure 13cl) regroupe cet ensemble de cellules. Ce groupement de cellules en forme de triangle de tétraèdre réduit le nombre de cellules à quadrillage sans réduire la précision de la simulation effectuée par le logiciel simulateur Eclipse 40c2 de la figure 10. 3. Génération de cellule - Bloc 40clC Sur la figure 13cl, la définition du terme "centre du cercle circonscrit" est donnée. Le "centre du cercle circonscrit" de chaque cellule d'une tétraédralisation (c'est-à-dire d'une cellule en forme de tétraèdre dans un quadrillage 3D) est le centre d'une sphère qui est tracée par les sommets de chaque cellule en forme de tétraèdre. Le "centre du cercle circonscrit" de chaque cellule d'une triangulation (c'est-à-dire d'une cellule triangulaire d'un quadrillage 2D) est le centre d'un cercle qui est tracé par les sommets de chaque cellule de forme triangulaire. Une cellule en forme de triangle ou une cellule en forme de tétraèdre ou bien encore une cellule de forme composite de triangle ou une cellule SR 15546 US/RS de forme composite de tétraèdre ou bien encore un agencement de l'une quelconque ou de plusieurs des cellules mentionnées ci-dessus correspondent chacune à une cellule de simulation perpendiculaire à K pourvu que le "centre du cercle circonscrit" soit utilisé comme centre de cellule de simulation. Une cellule en forme de triangle ou une cellule en forme de tétraèdre ou bien un agencement de celles-ci peuvent être utilisées comme cellules de simulation d'un quadrillage pourvu qu'on utilise les "centres des cercles circonscrits" comme centres de cellules de simulation. Lorsqu'un "assemblage" de cellules en forme de triangle ou de cellules en forme de tétraèdre d'un quadrillage est utilisé comme cellules de simulation, les centres des cercles circonscrits "combinés" de cet "assemblage" sont utilisés comme centres de cellules de simulation. 4. Calcul des propriétés des cellules - Bloc clD Sur la figure 13cl, un nouvel algorithme permettant de calculer la "transmissibilité" entre une paire de cellules approximativement perpendiculaires à K d'un quadrillage, par exemple entre les cellules adjacentes 15al du quadrillage représenté sur la figure 13a, est donné par la formule (17) reproduite ci- dessous. Transmissibilité -- La transmissibilité entre une paire de cellules i et j, pour des cellules TET 2D, PEBI 2D, TET 3D et PEBI 3D est donnée par les formules suivantes: SR 15546 US/RS T' = CD (17) l/t,+ t, o atK,r. (18) et o 0 NTG O K, = ONTG O Kit,,, (19) k 0 a = vecteur de surface de l'interface Ki = tenseur de perméabilité de la cellule i ajusté par net - brut ri = vecteur allant du centre de la cellule i à l'interface le long de la ligne des centres CD = constante de Darcy NTG = rapport d'épaisseur net - brut. La transmissibilité entre une paire de cellules radiales est calculée en utilisant la théorie de l'écoulement radial entre des centres moyens de pression. La transmissibilité entre une cellule radiale et une cellule linéaire est également calculée en utilisant l'écoulement radial entre les centres des cellules respectifs. 5. Numérotation de cellule - Bloc 40clE Sur la figure 13cl, le dispositif de résolution algébrique linéaire du logiciel simulateur Eclipse 40c2 fonctionne en utilisant des quadrillages "structurés" SR 15546 US/RS et non pas des quadrillages non structurés. Néanmoins, le quadrillage de la figure 13a utilisé par le logiciel Petragrid 40cl est un quadrillage non structuré. En conséquence, on a besoin d'un procédé permettant d'adapter le quadrillage non structuré à un quadrillage structuré destiné à être utilisé par le logiciel simulateur Eclipse 40c2. Par conséquent, selon un autre aspect de la présente invention, un nouvel algorithme associé au bloc de code "numérotation de cellule" clE "renumérote" chaque cellule d'un quadrillage non structuré et adapte chaque cellule du quadrillage non structuré numérotée dans un quadrillage structuré qui peut être utilisé par le logiciel simulateur Eclipse c2. En conséquence, le nouvel algorithme, appelé bloc de code "de renumérotation de cellule" 40c12 sur la figure 13cl "renumérote" chaque cellule du quadrillage non structuré en attribuant un indice (I, J, K) à "chaque cellule" du quadrillage non structuré. Un autre bloc de code appelé "exportateur de données" traité ci- dessous, adapte "chaque cellule" de quadrillage non structuré dans une cellule correspondante d'un quadrillage structuré, ladite cellule correspondante du quadrillage structuré se voyant attribuer une adresse composée dudit index (I, J, K) de chaque cellule. Dans la technique antérieure, le procédé correspondant comportait le placement d'un quadrillage rectangulaire présentant la même taille que la cellule sur la totalité de la géométrie et la numérotation des cellules dans le quadrillage non structuré en traversant le quadrillage rectangulaire dans un ordre de cellules logique. Toutefois, selon un autre aspect de la présente invention, un procédé de balayage de SR 15546 US/RS lignes basées sur le domaine est décrit pour adapter un quadrillage non structuré dans un quadrillage structuré pouvant être utilisé avec un simulateur de réservoir classique 40c2. La longueur de l'étape de balayage est déterminée de façon dynamique à partir des centres des cellules voisines. Rappelons qu'un quadrillage non structuré est composé d'un ensemble de domaines, dont certains sont structurés et certains non structurés, et le nouvel algorithme, associé à la "renumérotation des cellules" (bloc 40clE sur la figure 13cl) de la présente invention, numérote les domaines structurés dans un ordre logique. Les domaines non structurés sont numérotés comme une entité unique en balayant une ligne du quadrillage dans les sens X et Y. L'emplacement du centre de chaque cellule est utilisé au cours du balayage de ligne pour attribuer les indices I, J de cellules. K indices sont attribués en utilisant le nombre de cellules logiques, mais lorsqu'une numérotation logique n'est pas possible à effectuer, un balayage de ligne dans le sens Z peut être utilisé. Tous les domaines sont ensuite placés dans un quadrillage structuré global. Sur la figure 13c2, un organigramme plus complet du logiciel Petragrid 40cl des figures 10 et 13cl est illustré. Les nouvelles caractéristiques du logiciel Petragrid 40cl de la figure 13cl (blocs 40clA à 40clE) sont incluses dans l'organigramme de la figure 13c2. Sur la figure 13c2, le logiciel Petragrid 40cl comporte les blocs de codes suivants: un bloc de code appelé "dispositif d'élaboration de modèle interne" 40clF reçoit l'enregistrement de sortie de diagraphie 25 et l'enregistrement de sortie de données sismiques réduites 30d et, en réponse à cela, le dispositif SR 15546 US/RS d'élaboration du modèle interne 40clF génère une "structure de données de modèle interne 40clG". Le dispositif d'élaboration du modèle interne 40clF crée une "limite" dans laquelle les "points de la structure de données" de "points" 40clJ (générés par le "distributeur de points" 40clI) seront répartis. L'enregistrement de sortie de données sismiques réduites 30d définit l'endroit o les horizons 13 de la figure 1 et les failles sont situés à l'intérieur de la "limite" et l'enregistrement de sortie de diagraphie 25 définit l'endroit o les puits de forage 17 de la figure 1 sont situés à l'intérieur de la "limite". Toutefois, l'opérateur, au poste de travail 40 de la figure 9, définit la "limite" elle-même en définissant une ligne dans l'espace de surface. La structure de données du "modèle interne" 40clG comprend cette "limite". Un bloc de code appelé "distributeur de points selon les caractéristiques et contrôles" 40clI reçoit l'ensemble des données du modèle interne 40clG et répartit une pluralité de "points" à l'intérieur de la "limite" (par exemple un volume spatial en 3D pour les cellules en forme de tétraèdre); lorsque la pluralité de "points" est répartie de la "limite", une structure de données de "points" 40clJ est générée. Toutefois, un bloc de code appelé "éditeur de paramètre de contrôle" 40clH va: (1) contrôler la distance entre les "points" adjacents répartis à l'intérieur de la "limite" ce qui permet de définir un quadrillage précis ou un quadrillage grossier, (2) définir le nombre de ces "points" dans le sens X et le nombre de ces "points" dans le sens Y, et (3) définir le nombre de divisions radiales ou "anneaux" autour du puits de SR 15546 US/RS forage. Un bloc de code appelé "tétraédraliseur triangulateur automatique" 40clA (voir bloc 40clA de la figure 13cl) recevra la structure de données de "points" 40clJ qui contient une multitude de "points" à l'intérieur de la "limite" et connecte ensemble tous les "points" de la structure de données de "points" clJ de façon à créer une "pluralité de triangles" (pour un espace 2D) et "une pluralité de tétraèdres" (pour un espace 3D), ce qui permet de générer une structure de données "de tétraèdres et de triangles" clK. Un bloc de code appelé "agrégateur de triangles/tétraèdres" 40clB (voir bloc 40clB de la figure 13cl) va: (1) recevoir la "pluralité de triangles" et/ou la "pluralité de tétraèdres" dans la structure de données "de tétraèdres et de triangles" clK, et (2) produire une liste, appelée structure de données "de cartes d'agrégation" 40clL, qui établit quel ensemble de triangles parmi la "pluralité de triangles" et quel ensemble de tétraèdres de la "pluralité de tétraèdres" doivent être regroupés pour former une forme de polygone. Un bloc de code appelé "générateur de cellule" 40clC (voir bloc 401c de la figure 13cl) reçoit la structure de données "de cartes d'agrégation" 40clL et la structure de données "de modèle interne" 40clG et, en réponse à cela, le générateur de cellule 40clC fait le groupement mentionné ci-dessus; c'est-à-dire que le générateur de cellule 40clC regroupe "un premier ensemble de triangles" de la "pluralité de triangles" et "un premier ensemble de tétraèdres" de la "pluralité de tétraèdres" selon la liste de triangles et de tétraèdres mentionnés dans la structure de données de SR 15546 US/RS cartes d'agrégation 40clL, le "premier ensemble de triangles" et le "premier ensemble de tétraèdres" étant regroupés à l'intérieur de la "limite", c'est-à-dire définis par la structure de données de modèle interne 40clG, ce qui permet de définir une structure de données "de quadrillage non structuré" 40clM (voir le quadrillage de la figure 13a). Un bloc de code appelé "numéroteur de cellule" 40clE (voir bloc 40clE de la figure 13cl) attribue un indice (I, J, K) à chaque "cellule" dans la structure de données de quadrillage non structuré 40clM en réponse à la structure de données de modèle interne 40clG. Un bloc de code appelé "exportateur de données" 40clQ reçoit à la fois l'indice (I, J, K) mentionné ci-dessus pour chaque "cellule" du quadrillage non structuré et la structure de données de quadrillage non structuré 40clM et, en réponse à cela, l'exportateur de données 40clQ adapte chaque cellule de la structure de données du quadrillage non structuré 40clM dans un quadrillage "structuré, ce qui permet de produire un quadrillage structuré ayant une pluralité de cellules, chaque cellule du quadrillage structuré présentant une adresse d'indice (I, J, K). Le quadrillage structuré est utilisé par le logiciel simulateur Eclipse 40c2 (rappelons que le simulateur Eclipse 40c2 ne peut fonctionner que sur des quadrillages "structurés"). Un bloc de code appelé "calculateur de propriétés de cellules" 40clD (voir bloc 40clD de la figure 13cl) répond également à l'ensemble des données du modèles interne 40clG en calculant la "transmissibilité" et le "volume poreux" entre "chaque paire des cellules perpendiculaires à K du quadrillage non structuré SR 15546 US/RS clM. Rappelons que la "transmissibilité" est calculée en utilisant la formule suivante: Ti CD (17) t- l/t, + t Puisqu'il y a une multitude de paires de cellules perpendiculaires à K dans la structure de données du quadrillage non structuré 40clM, une multitude correspondante de "transmissibilité (T)" sera calculée par le calculateur de propriétés de cellules 40clD. Une "visionneuse 3D" 40clN est le moniteur d'affichage du poste de travail 40d de la figure 9 qui affiche la structure de données du quadrillage "non structuré" clM. Un bloc de code appelé "calculateur de connexion de puits" 40clP répond à la structure de données du quadrillage non structuré 40clM et à la structure de données du modèle interne 40clG en générant une structure de données "de connexion de puits" 40clR. Etant donné que chacun des puits de forage 17 de la figure 1 a une intersection avec une pluralité de cellules du quadrillage non structuré 40clM, le calculateur de connexion de puits 40clP calcule l'écoulement du fluide (appelé "coefficient de débit") entre chaque cellule du quadrillage non structuré 40clN qui a une intersection avec un puits de forage 17. La structure de données de connexion de puits 40clR comporte par conséquent, pour chaque puits de forage 17 de la figure 1, une liste particulière comprenant toutes les cellules du quadrillage non structuré 40clM qui sont coupées par et qui ont une intersection avec le puits de forage 17 de la figure 1, et une liste correspondante de coefficients de débit (appelés SR 15546 US/RS facteurs de connexion de puits) qui sont associés, respectivement, à cette liste particulière de cellules. L'exportateur de données 40clQ répond à la structure de données du quadrillage non structuré 40clM et à l'ensemble des données de connexion de puits 40clR en adaptant chaque cellule de la structure de données du quadrillage non structuré 40clM dans un quadrillage "structuré", ce qui permet de produire une structure "de données de sortie" 40c3 qui est composée du quadrillage structuré ayant une pluralité de cellules, chaque cellule du quadrillage structuré présentant une adresse d'indice (I, J, K) unique. Un élément très important de la structure "de données de sorties" 40c3 est la "transmissibilité" entre chaque paire de cellules perpendiculaires à K du quadrillage structuré et le volume poreux. Ces "données de sortie" 40c3 sont utilisées par le simulateur Eclipse 40c2. On peut faire référence à la figure 13dl ci-dessous pour comprendre mieux l'utilisateur des "données de sortie". (3) Résumé du logiciel simulateur Eclipse 40c2 Sur la figure 13dl, une construction plus détaillée du logiciel simulateur Eclipse 40c2 de la figure 13c2 est illustrée. Le logiciel simulateur Eclipse 40c2 de la figure 13dl est utilisé conjointement avec le mode de réalisation préféré de cette demande. Il faut comprendre, cependant, que d'autres simulateurs peuvent être utilisés conjointement avec le logiciel Petragrid 40cl de la présente invention. On peut se procurer le logiciel simulateur Eclipse 40c2 auprès de "GeoQuest", division de Schlumberger Technology Corporation", Houston, Texas ou Abingdon, Royaume-Uni. SR 15546 US/RS Sur la figure 13dl, les données de sortie 40c3 (comportant la transmissibilité traitée ci-dessus) sont entrées dans la partie "simulateur 40c2A du logiciel simulateur Eclipse 40c2. Cette partie du simulateur 40c2A reçoit également la sortie d'autres programmes, par exemple le programme PVT 40c2B, le programme VFP c2C, le programme Schedule 40c2D et d'autres programmes 40c2E, tous ces programmes pouvant être fournis par "GeoQuest", division de Schlumberger Technology Corporation", Houston, Texas. En réponse aux données de sortie 40c3, le simulateur 40c2A génère des résultats de simulation 40c2F pour chaque cellule du quadrillage structuré, par exemple pour chaque cellule al du quadrillage de la figure 13b. Parmi les exemples de "résultats de simulation" 40c2F, on peut citer une "pression" et une "saturation" dans chaque cellule de quadrillage structuré. Ces résultats de simulation 40c2F sont utilisés par des programmes séparés, par exemple le programme GRAF 50, le programme RTView 52 et d'autres programmes 54. Toutefois, ce qui est plus important, ces résultats de simulation 40c2F sont retransmis à la visionneuse 3D 40clN qui représente le moniteur d'affichage 40d du poste de travail 40 de la figure 9. Rappelons, d'après la figure 13c2, que le quadrillage non structuré clM (qui comprend une pluralité de cellules 15al) est transmis à la et affiché sur la visionneuse 3D 40clN. Par conséquent, le quadrillage non structuré comportant toutes ces cellules 15al sera affiché sur la visionneuse 3D 40clN (moniteur d'affichage 40d du poste de travail 40 de la figure 9). Sur la figure 13dl, les résultats de simulation 40c2F, pour chaque cellule 15al du quadrillage, sont également envoyés à la visionneuse SR 15546 US/RS 3D 40clN (affichage 40d). En conséquence, les résultats de simulation 40c2F pour "chaque cellule 15al sont placés sur et affichés dans "chaque cellule" 15al du quadrillage sur la visionneuse 3D 40clN (moniteur d'affichage 40d) du poste de travail 40. En fonction de la valeur de la pression ou de la saturation des résultats de simulation 40c2F, une couleur différente peut être attribuée aux résultats de simulation 40c2F; et, lorsque les résultats de simulation 40c2F sont placés sur "chaque cellule" particulière 15al sur la visionneuse 3D 40clN, la couleur attribuée associée aux résultats de simulation particuliers 40c2F (par exemple pression et saturation) sera utilisée pour afficher la valeur réelle des résultats de simulation 40c2F de "chaque cellule" 15al sur la visionneuse 3D 40clN (affichage 40d). Sur la figure 13d2, une "sortie typique" qui est affichée sur la visionneuse 3D 40clN des figures 13c2 et 13dl (c'est-à-dire sur "l'affichage" 40d de la figure 9) par le logiciel simulateur Eclipse 40c2est illustrée. Sur la figure 13d2, lorsque la couleur attribuée associée aux résultats de simulation particuliers c2F est utilisée pour afficher les résultats de simulation 40c2F de "chaque cellule" 15al sur la visionneuse 3D 40clN (affichage 40d), un "affichage de sortie" semblable à l'affichage de sortie 56 représenté sur la figure 13d2, est affiché sur le moniteur d'affichage 40d (visonneuse 3D 40clN) du poste de travail d'interprétation 40 de la figure 9. Une description fonctionnelle du fonctionnement du logiciel Petragrid 40cl de la présente invention, conjointement avec le logiciel simulateur Eclipse 40c2, SR 15546 US/RS sera décrite dans les paragraphes suivants en référence aux figures 1 à 13d2 des dessins. Les opérations de diagraphie des figures 1 et 2 sont exécutées et un enregistrement de sortie de diagraphie 25 est généré. De plus, les opérations sismiques des figures 4 et 5 sont effectuées dans la formation terrestre 15 près des puits de forage 17 de la figure 1 et l'enregistrement de sortie sismique 20a4 de la figure 6 est généré. L'enregistrement de sortie sismique 20a4 de la figure 6 subit une réduction de données dans l'ordinateur à unité centrale 30 de la figure 7 (en utilisant le logiciel représenté sur les figures 11 et 12) et un enregistrement de sortie de données sismiques réduites 30d de la figure 7 est généré. L'enregistrement de sortie de diagraphie 25 et l'enregistrement de sortie de données sismiques réduites 30d sont entrés dans un poste de travail d'interprétation 40 des figures 8 et 9. Un disque de stockage appelé "CD-Rom" peut, par exemple, stocker séparément le logiciel Petragrid 40cl et/ou le logiciel simulateur Eclipse 40c2. Sur la figure 9, ce CD-Rom est inséré par un opérateur, dans le poste de travail 40 de la figure 9 et, en conséquence, le logiciel Petragrid cl et/ou le logiciel simulateur Eclipse 40c2 est chargé dans la mémoire 40c du poste de travail 40 de la figure 9. La figure 10 illustre le logiciel Petragrid cl et le logiciel simulateur Eclipse 40c2 stockés dans la mémoire 40c du poste de travail 40. En réponse à une action de l'opérateur au poste de travail 40, le processeur du poste de travail 40b de la figure 9 commence à exécuter le logiciel Petragrid 40cl et le logiciel simulateur Eclipse 40c2 de la figure 10, le processeur du poste de travail 40b utilisant SR 15546 US/RS l'enregistrement de sortie de diagraphie 25 et l'enregistrement de sortie de données sismiques réduites 30d qui ont déjà été chargés dans le poste de travail 40. Lorsque cette exécution du logiciel Petragrid C1 et du logiciel simulateur Eclipse 40c2 est terminée, la formation terrestre 15 de la figure 1, y compris les horizons 13 encastrés dans la formation 15, est divisée en une multitude de "cellules" ou volumes en forme de tétraèdres interconnectés 15al qui sont représentés sur la figure 13a, et les "informations plus précises" associées à "chaque cellule" 15al, concernant les "propriétés de cellules" de "chaque cellule" 15al (y compris en particulier la "transmissibilité") seront générées. Les "informations plus précises" générées par le logiciel Petragrid 40cl sont utilisées par le logiciel simulateur Eclipse 40c2 pour générer les "résultats de simulation plus précis" c2F de la figure 13dl. La figure 13b illustre une vue éclatée d'une couche d'horizon de la figure 13a, illustrant plus en détail la multitude de cellules en forme de tétraèdre al. Bien que les cellules 15al de la figure 13b présentent une coupe approximativement "rectangulaire", lorsque l'exécution du logiciel Petragrid 40cl est achevée, les cellules 15al présentent une forme approximative triangulaire dans l'espace bidimensionnel ou une forme approximative de tétraèdre dans l'espace tridimensionnel, tel que représenté sur la figure 13a. Lorsque le processeur du poste de travail 40b de la figure 9 commence à exécuter le logiciel Petragrid cl et le logiciel simulateur Eclipse 40c2 de la figure 10, le dispositif d'élaboration du modèle SR 15546 US/RS interne clF de la figure 13c2 élabore une structure de données de modèle interne clG qui comprend une "limite" dans laquelle une multitude de "points" est répartie. Le répartiteur de points clI répartit la multitude de "points" dans la "limite" définie par la structure de données du modèle interne clG, la distance entre les "points" adjacents ainsi que le nombre de "points" dans le sens X et le nombre de "points" dans le sens Y étant contrôlé par l'éditeur de paramètre de contrôle 40clH. Le répartiteur de points 40clI génère une structure de données "de points" clJ qui est composée d'une multitude de points répartis dans la "limite". Le tétraédralisateur/triangulateur automatique 40clA répond à la multitude de points dans la structure de données de points 40clJ en établissant une interconnexion entre ces points de façon à former une multitude de cellules de forme triangulaire (pour un espace bidimensionnel) ou une multitude de cellules en forme de tétraèdre (pour un espace tridimensionnel), ce qui permet de générer la structure de données en tétraèdres/triangles 40clK. Plus particulièrement, le tétraédralisateur/triangulateur automatique 40clA de la figure 13c2 va: (1) générer le quadrillage de base approximativement triangulaire et "perpendiculaire à K" mentionné ci-dessus pour les quadrillages non structurés bidimensionnels (2D) en "mettant à l'échelle" la géométrie physique et (2) générer le quadrillage de base présentant approximativement la forme d'un tétraèdre et "perpendiculaire à K" pour des quadrillages non structurés tridimensionnels (3D) en "mettant à l'échelle" la géométrie physique et va générer les deux types de quadrillages mentionnés ci- SR 15546 US/RS dessus en exécutant les étapes fonctionnelles suivantes: (1) la géométrie physique est adaptée dans un espace de calcul en utilisant une transformation fondée sur les composants du tenseur de perméabilité dans chaque couche du modèle, (2) des points sont répartis dans l'espace de calcul en fonction de la géométrie adaptée; (3) une tétraédralisation ou une triangulation de Delaunay des points est créée, dans laquelle: (3a) le centre de chaque élément (c'est-à- dire chaque "cellule") de la triangulation correspond au centre d'un cercle qui est tracé par les sommets de chaque cellule, et (3b) le centre de chaque cellule de la tétraédralisation correspond au centre (ou centre du cercle circonscrit) d'une sphère qui est tracée par les sommets des éléments ou cellules, et (4) les sommets des éléments ou cellules ou centres des cellules sont ensuite retransformés dans l'espace physique. Lorsque ces étapes fonctionnelles sont achevées, la structure de données de tétraèdre/triangle C1K est créée. L'agrégateur de tétraèdre/triangle 40clB de la figure 13c2 reçoit la "pluralité de triangles" et/ou la "'pluralité de tétraèdres" dans la structure de données de tétraèdres et de triangles 40clK et produit une liste appelée "structure de données" de cartes d'agrégation 40clL, qui établit les ensembles de triangles de la "pluralité de triangles" et un ensemble de tétraèdres de la "pluralité de tétraèdres" doivent être regroupés pour former une forme de polygone. Le générateur de cellule 40clC de la figure 13c2 reçoit la structure de données de "carte d'agrégation" 40clL et 1 structure de données "de modèle interne" 40clG et en réponse à cela, le générateur de cellule 40clC regroupe SR 15546 US/RS "un premier ensemble de triangles" de la "pluralité de triangles" et "un premier ensemble de tétraèdres" de la "pluralité de tétraèdres" en fonction de la liste de triangles et de tétraèdres mentionnée dans la structure de données de carte d'agrégation 40clL, le "premier ensemble de triangles" et le "premier ensemble de tétraèdres" étant regroupés de la "limite" qui est définie par la structure de données du modèle interne clG, ce qui permet de définir une structure de données "de quadrillage non structuré" 40clM. Le numéroteur de cellule 40clE attribue un indice (I, J, K) à chaque "cellule" de la structure de données de quadrillage non structuré 40clM en réponse à la structure de données du modèle interne 40clG. L'exportateur de données 40clQ reçoit à la fois l'indice (I, J, K) mentionné ci-dessus pour chaque "cellule" de quadrillage non structuré et la structure de données de quadrillage non structuré 40clM, et, en réponse à cela, l'exportateur de données 40clQ adapte chaque cellule de la structure de données du quadrillage non structuré 40clM dans un quadrillage "structuré", ce qui permet de produit un quadrillage structuré et une pluralité de cellules, chaque cellule du quadrillage structuré présentant une adresse d'indice (I, J, K). Le calculateur de propriétés de cellules 40clD répond également à l'ensemble des données de modèle interne 40clG en calculant la "transmissibilité" et le "volume poreux" en "chaque paire" de cellules perpendiculaires à K du quadrillage non structuré 40clM. Rappelons que la "transmissibilité" est calculée en utilisant la formule suivante: SR 15546 US/RS T7= CD (17) A/t, + 1/ta Puisqu'il y a une multitude de paire de cellules perpendiculaires à K dans la structure de données de quadrillage non structuré clM, une multitude correspondante de "transmissibilité (T)" sera calculée par le calculateur de propriétés de cellules 40clD. La visionneuse 3D 40clN correspond au moniteur d'affichage du poste de travail 40d de la figure 9 qui affiche la structure de données de quadrillage "non structuré" clM. Le calculateur de connexion de puits 40clP répond à la structure de données de quadrillage non structuré 40clM et à la structure de données du modèles interne 40clG en générant la structure de données de connexion de puits 40clR. Etant donné que chacun des puits de forage 17 de la figure 1 présente une intersection avec une pluralité de cellules du quadrillage non structuré 40clM, le calculateur de connexion de puits 40clP calcule l'écoulement de fluide (appelé "coefficient de débit") entre chaque cellule du quadrillage non structuré clM qui présente une intersection avec un puits de forage 17. La structure de données de connexion de puits 40clR comporte par conséquent, pour chaque puits de forage 17 de la figure 1, une liste particulière de toutes les cellules du quadrillage non structuré 40clM qui sont coupées par et qui ont une intersection avec le puits de forage 17 de la figure 1 et une liste correspondante de coefficients de débit (appelé facteurs de connexion de puits) qui sont associés, respectivement, à chaque liste particulière de cellules. L'exportateur de données clQ répond à la structure de données du quadrillage SR 15546 US/RS non structuré 40clM et à l'ensemble des données des connexion de puits 40clR en adaptant chaque cellule de la structure de données du quadrillage non structuré clM dans un quadrillage "structure", ce qui permet de produire une structure "de données de sortie" 40c3 qui est composée du quadrillage structuré ayant une pluralité de cellules, chaque cellule du quadrillage structuré ayant une adresse d'indice (I, J, K) unique. La structure "de données de sortie" 40c3 comporte une "transmissibilité" entre chaque paire de cellules perpendiculaires à K du quadrillage structuré et le volume poreux. Les données de sortie 40c3 sont utilisées par le simulateur Eclipse 40c2. Sur la figure 13dl, le quadrillage non structuré 40clM est affiché sur la visionneuse 3D 40clN (affichage 40d de la figure 9). Le simulateur 40c2A reçoit les données de sortie 40c3 et en particulier les chiffres de transmissibilité et de volume poreux, et il génère ensuite un résultat de simulation 40c2F (par exemple la pression ou la saturation) pour "chaque cellule" du quadrillage affiché sur la visionneuse 3D clN. Les résultats de simulation 40c2F pour "chaque cellule" sont transmis à la visionneuse 3D 40clN et les résultats de simulation 40c2F pour "chaque cellule" sont placés sur "chaque cellule" sur la visionneuse 3D clN. Une couleur différente est attribuée à chaque valeur différente du résultat de simulation 40c2F et cette couleur est affichée à l'intérieur de "chaque cellule" sur la visionneuse 3D 40clN. Par exemple, pour une valeur de pression d'une cellule particulière, une première couleur est attribuée et pour une valeur différente de pression d'une autre cellule, une deuxième couleur est attribuée. La première couleur SR 15546 US/RS représente l'intérieur de la cellule particulière et la deuxième couleur représente l'intérieur de ladite autre cellule. En conséquence, la visionneuse 3D 40clN (ou moniteur d'affichage 40d de la figure 9) présente un affichage de sortie permettant d'être visualisé par un opérateur. Un exemple d'un tel affichage de sortie du simulateur d'Eclipse 40c2 est illustré sur la figure 13d2. Les données de sortie obtenues, qui sont visibles et qui peuvent être obtenues sur l'un quelconque des affichages de sortie du simulateur Eclipse 40c2, par exemple l'affichage de sortie de la figure 13d2, sont "plus précises" que celles de la technique antérieure, principalement du fait que les "données de sortie" 40c3 des figures 13c2 et 13cl, qui sont générées par le logiciel Petragrid 40cl de la présente invention, sont "plus précis". Deuxième section - Discussion générale sur le programme de quadrillage de simulation "Flogrid" comportant le logiciel "de quadrillage structuré" qui comporte en outre le "dispositif de quadrillage de surface structuré" de la présente invention En référence aux figures 14a à 14c de la figure 14a, le poste de travail d'interprétation 40 de la figure 9 est de nouveau illustré et comporte le moniteur, le processeur, le clavier et la souris. Sur la figure 14b, le poste de travail 40 est de nouveau représenté et comporte le processeur de station de base b qui est connecté au bus du système, l'affichage 40b connecté au bus et la mémoire du poste de travail 40c connectée au bus. Le bus reçoit l'enregistrement de SR 15546 US/RS sortie de diagraphie 25 et l'enregistrement de sortie de données sismiques réduites 30d. Un dispositif de stockage 62, conçu pour stocker les informations logicielles, appelé également CD-Rom 62, peut être inséré dans le poste de travail 40. Le bloc 60 de la figure 14b indique que le CD-Rom 62 est conçu pour stocker le programme de quadrillage de simulation "Flogrid" 60a, pour générer les informations de propriétés de cellules de quadrillage de formation terrestre plus précises et le logiciel simulateur "Eclipse" 40c2 qui reçoit les "informations de propriétés de cellules de quadrillage de formation terrestre plus précises" et, en réponse à cela, génère "l'ensemble plus précis" des résultats de simulation". Lorsque le CD-Rom 62 est inséré dans le poste de travail 40 de la figure 14a, le logiciel stocké dans le CD-Rom 62 comprenant le programme de quadrillage de simulation Flogrid qui comporte en outre le dispositif de quadrillage "non structuré" Petragrid et le logiciel de quadrillage de surface "structuré" de la présente invention, va être chargé dans la mémoire du poste de travail 40c. En conséquence, sur la figure 14c, la mémoire du poste de travail 40c des figures 9 et 14b stocke les paquets de logiciels suivants: le programme de quadrillage de simulation "Flogrid" 60a et le logiciel simulateur Eclipse 40c2 décrit dans la demande de Gunasekera. Toutefois, le programme de quadrillage de simulation "Flogrid" 60a comporte les blocs logiciels suivants: un dispositif de stockage de données de réservoir 64, un corps de réservoir 66, le logiciel de quadrillage "structuré" 68 selon la présente invention, le logiciel de quadrillage "non SR 15546 US/RS structuré" Petragrid 40cl décrit dans la demande de Gunasekera et un dispositif de mise à une plus grande échelle 70. En référence à la figure 15, une description fonctionnelle du fonctionnement du programme de quadrillage de simulation Flogrid 60a et du logiciel simulateur Eclipse 40c2 de la figure 14c est donnée dans le paragraphe suivant en référence à la figure 15 des dessins. Sur la figure 15, un ensemble de données de sortie 69, comportant l'enregistrement de sortie de diagraphie 25 et l'enregistrement de sortie de données sismiques réduites 30d de la figure 14b (en plus d'autres données de sorties comportant des affleurements et des dessins) est fourni en tant qu'entrée au programme de quadrillage de simulation "Flogrid" 60a. Plus particulièrement, les données d'entrée 69 sont entrées dans le canevas du réservoir Flogrid 66 de la figure 14c, le canevas du réservoir 66 comportant un paquet d'adaptations 65 et un géomodélisateur 67. Le canevas du réservoir 66 produit un modèle structurel 74 (qui comporte des cellules, des puits et des cartes d'horizons) qui est fourni comme entrée au dispositif de quadrillage "Flogrid" 68, 40cl. En réponse au modèle structurel 74, le dispositif de quadrillage "Flogrid" 68, 40cl (inclus dans le dispositif de quadrillage non structuré Petragrid 40cl et le dispositif de quadrillage structuré 68 de la présente invention) sert à produire un modèle de simulation quadrillé 76, le modèle de simulation 76 étant une section d'une formation terrestre quadrillée semblable à la section quadrillée représentée sur la figure 13b. Le modèle de simulation 76 est composé SR 15546 US/RS d'une pluralité de cellules de quadrillage (en coupe) présentant une forme approximativement rectangulaire, semblable aux cellules de quadrillage 15al représentées sur la figure 13b. Le modèle de simulation 76 est fourni comme entrée au simulateur Eclipse 40c2, avec d'autres données d'ingénierie 78. Le simulateur 40c2 génère les résultats de simulation 40c2F (voir figure 13dl) qui seront affichés sur la visionneuse 3D 40clN. Si les résultats de simulation 40c2F ne correspondent pas à certaines données historiques 80, les résultats de simulation 40c2F sont renvoyés, par l'intermédiaire de la ligne 82, au canevas du réservoir 66 afin de manipuler les résultats de simulation générés ci-après c2F par le simulateur 40c2. Le programme de quadrillage de simulation Flogrid a des figures 14c et 15 En référence à la figure 16, une construction plus détaillée du programme de quadrillage de simulation Flogrid 60a de la figure 14c et de la figure est illustrée. Sur la figure 16, les données d'entrée 69 de la figure 15 comprenant l'enregistrement de sortie de diagraphie 25 (piste de puits) et l'enregistrement de sortie de données sismiques réduites 30d (cartes et failles) ainsi que d'autres données d'entrée, par exemple des enregistrements SCAL, sont entrées et stockées dans le "dispositif de stockage de données de réservoir" Flogrid 64. On peut se procurer le programme de quadrillage de simulation Flogrid 60a de la figure auprès de "GeoQuest", division de Schlumberger Technology Corporation", situé à Houston, Texas. SR 15546 US/RS Le rôle du dispositif de stockage de données de réservoir 64 est de: importer, stocker, visualiser et classer des données pour qu'elles puissent être utilisées ultérieurement dans des modules en aval. Les classifications dépendent de la quantité exacte d'informations fournies par les données importées. Les options comprennent: (1) le fait de faire en sorte que les cartes correspondantes soient d'un type spécifique (par exemple porosité, perméabilité, horizon), (2) le fait d'établir ou de modifier les unités de cartes et de systèmes de coordonnées, (3) le fait de créer des surfaces - les surfaces représentent les couches supérieures, inférieures et intermédiaires et sont classées; l'utilisation de ces informations est important pour "l'automatisation" de la construction des unités de réservoir (strates) et des modèles structuraux en 3D obtenus ci-après; (4) le fait d'attribuer des cartes à une surface particulière (par exemple établir quelles sont les cartes qui représentent les couches supérieures, inférieures et intermédiaires), (5) le fait d'attribuer des cartes d'épaisseur à des surfaces qu'elles représentent, et (6) le fait d'attribuer des traces de failles à des surfaces particulières. Des facilités d'édition des données limitées sont également fournies dans le dispositif de stockage de données de réservoir 64. Les options se limitent en général à la création et à l'édition de traces de failles et de polygones dans la visionneuse 3D 40clN. Les données stockées dans le dispositif de stockage de données de réservoir 64 sont stockées dans le corps du réservoir 66. SR 15546 US/RS Le canevas du réservoir 66 de la figure 16 est composé de deux parties: (1) un modèle structurel 66a et (2) un modèle de propriété 66b. Le modèle structurel 66a va: construire "des modèles géologiques en 3D" à partir des données d'entrée 69, comportant des cartes (c'est-à-dire des horizons) et des failles. Le modèle de propriété 66b va: importer lesdits "modèles géologiques en 3D" provenant du modèle structurel 66a, créer et éditer les propriétés dans lesdits modèles géologiques en 3D et, en réponse à cela, le modèle de propriété 66b va créer et générer le "modèle structurel" mentionné ci-dessus 74. Le modélisateur structurel 66a de la figure 16 génère un "modèle géologique en 3D" 86. Le modélisateur structurel 66a est décrit dans les deux demandes de brevets américains en suspens: (1) la demande américaine antérieure en suspens numéro de série 08/823.107 attribuée à Abbott (ci-après appelée "demande d'Abbott"), déposée le 24 mars 1997 et intitulée "Method and Apparatus for Determining Geologic Relationships for Intersecting Faults". La demande d'Abbott décrit un appareil et un procédé fondés sur un logiciel et un poste de travail qui reçoit des informations concernant une multitude de failles dans une formation terrestre, détermine si chaque paire particulière de failles des données d'entrée "se rencontre" (c'est-à-dire si les paires de failles se coupent), détermine quelle faille des paires de failles qui se coupent constitue la faille "majeure" et quelle faille constitue la faille "mineure", détermine si la faille mineure coupe la faille majeure par le haut ou par le bas de la faille majeure, et, en réponse à cela, génère des "données de relation et de SR 15546 US/RS surface de faille"; et (2) la demande provisoire américaine antérieure en suspens numéro de série /029.524, déposée le 31 octobre 1996 intitulée "Method for 3D Modeling of Faulted Geologic Horizons", et la demande américaine antérieure en suspens numéro de série 08/916.841 attribuée à Graf et al (appelée ci- après "demande de Graf"), déposée le 22 août 1997 et intitulée "Automatic Non-Artificially Extended Fault Surface Based Horizon Modeling System". La demande de Graf décrit également un appareil et un procédé fondés sur un logiciel et un poste de travail qui reçoit les données "de relation et de surface de faille" générées par la "demande d'Abbott" identifiée ci-dessus et, en réponse à cela, elle génère le "modèle géologique en 3D" 86, également appelé "modèle d'horizon faillé final", qui représente un modèle structurel tridimensionnel élaboré à partir des données de cartes (c'est-à-dire des données d'horizons) et des données de failles. Plus particulièrement, le modélisateur structurel 66a reçoit les données choisies dans le dispositif de stockage de données de réservoir 64 et, en réponse à cela, construit des modèles géologiques en 3D stables à partir desquels les modèles de simulation sont élaborés. La construction d'un modèle géologique en 3D approprié est en général une tâche qui demande de nombreux calculs, et qui exige normalement que des surfaces importantes soient calculées et recalculées. Les surfaces peuvent devoir être calculées pour une variété de raisons. On peut citer: l'interpolation ou l'extrapolation des zones vides, le calcul de nouvelles surfaces à partir des données d'épaisseur ou le nouveau calcul de surface au voisinage de failles afin de SR 15546 US/RS corriger des surfaces incorrectes ou pour fournir des surfaces extrapolées afin d'échantillonner le mauvais côté de la faille. Les constructions d'un modèle structurel 66a impliquent trois étapes principales: définition de la limite, création de l'unité et fente du bloc de failles. La définition de la limite permet d'ignorer les parties inutiles des cartes, par exemple les zones ne contenant aucun hydrocarbure ou les grandes zones vides situées dans son périmètre, lors des étapes de construction et de nettoyage du modèle suivant. En ce qui concerne la création d'unité, elle permet d'élaborer des unités de réservoir à partir de représentations de surface et de surface choisies stockées dans le dispositif de stockage de données de réservoir 64. Les opérations utilisant des cartes maillées automatisées calculent des unités fondées sur l'ordre donné de surfaces, et permettent d'éroder les surfaces qui se chevauchent en fonction de leurs types de surface respectifs. En ce qui concerne le fractionnement de bloc de failles, cette étape achève le procédé de construction du modèle et donne une solution au problème complexe que pose l'échantillonnage des surfaces (et les propriétés qui en découlent) au voisinage de failles et de limites de couches. Le fractionnement de bloc de failles divise d'abord le modèle en un ensemble de blocs de failles, en extrapolant des failles si nécessaire. Ensuite, chaque bloc est en outre subdivisé en unités de blocs, en utilisant les unités créées précédemment. Enfin, les surfaces d'unités de blocs sont réparées et extrapolées pour produire un modèle qui supporte un échantillonnage SR 15546 US/RS rigoureux, comportant des valeurs d'échantillonnage provenant du "mauvais" côté de la faille. En conséquence, sur la figure 16, le "modèle géologique en 3D" 86 est fourni par le modélisateur structurel 66a et transmis au modélisateur de propriété 66b (le "modèle géologique 3D" 86 correspondant au "modèle d'horizon faillé final" généré par le logiciel de modélisation d'horizon décrit dans la "demande de Graf", identifiée ci-dessus). Le modélisateur de propriété 66b reçoit le modèle géologique 3D 86, construit un "canevas structurel" à partir du modèle géologique 3D 86 (c'est-à-dire à partir du "modèle d'horizon faillé final" décrit dans la demande de Graf) et, en réponse à cela, le modèle de propriété 66b génère un "modèle structurel" à partir du canevas structurel. Le "modèle structurel" fournit les connexions physiques explicites entre les surfaces de failles/horizons et leurs relations géologiques, l'existence de ces connexions physiques permettant à l'utilisateur d'éditer les étapes précédentes lors de la production du modèle structurel. Plus particulièrement, on propose un volume intéressant de canevas structurel, le canevas étant défini par des surfaces géologiques et sismiques. Les surfaces disponibles et les volumes existants du volume intéressant sont analysés et des surfaces candidates au fractionnement sont choisies pour une autre analyse. Les surfaces peuvent être analysées et les résultats validés. Les volumes sont sous-divisés en blocs qui dépendent de critères fournis par l'utilisateur. Le canevas structurel peut être édité par l'utilisateur, spécialement en sautant au moins une étape précédente. En particulier, le modélisateur de propriété 66b unifie SR 15546 US/RS des approches de modélisations géologiques en 3D et fondées sur des cartes pour donner une représentation de propriétés et structurelle en 3D utiles à partir de laquelle il est possible d'élaborer le dispositif de quadrillage de simulation. Ce module de modèle de propriété 66b a trois fonctions principales: (1) importer des modèles géologiques en 3D 86, (2) attribuer des cartes (horizons) à des unités créées dans un modèle structurel - ceci permet à l'utilisateur qui n'a pas indiqué que des cartes de propriété s'appliquent au- dessus ou en dessous une surface donnée dans le dispositif de stockage de données de réservoir 64 d'attribuer directement des cartes aux unités, (3) créer des propriétés supplémentaires sur un modèle de propriétés - les propriétés peuvent être définies par des expressions algébriques simples, par exemple des fonctions simples de propriétés existantes et avec des fonctions complexes définies par l'utilisateur en utilisant un calculateur. Le modélisateur de propriétés 66b répond à un dispositif de quadrillage progressif 84. C'est-à-dire que, tandis que l'on exécute le modèle de propriété 66b, le logiciel de quadrillage progressif 84, sur la figure 16, permet à l'utilisateur/opérateur de regrouper des couches (à l'intérieur d'un modèle à échelle précise) dans des couches grossières. Ceci est effectué manuellement ou automatiquement en utilisant des procédés qui limitent le flux total à travers la couche grossière ou limitent les variations de certaines propriétés à l'intérieur d'une couche. Les résultats générés par le modélisateur de propriétés 66b de la figure 16 sont entrés dans le dispositif de quadrillage "non structuré" 40cl. Le SR 15546 US/RS dispositif de quadrillage non structuré 40c1 de la figure 16 est le logiciel de quadrillage non structuré Petragrid qui est décrit dans la "demande de Gunasekera" identifiée ci-dessus et qui est décrit, ci- dessus dans la première section de la description du mode de réalisation préféré de la présente demande. Le dispositif de quadrillage non structuré 40cl de la figure 16 quadrille les caractéristiques du réservoir en utilisant des quadrillages PEBI et des quadrillages en forme de tétraèdre (et non pas de rectangle), des propriétés d'échantillons, des résultats d'importation, des propriétés d'édition et calcule des facteurs de connexion de puits. Les résultats générés par le modélisateur de propriétés de la figure 16 sont également fournis comme entrée au dispositif de quadrillage "structuré" 68 qui comporte le "dispositif de quadrillage de surface structuré" de la présente invention. Sur la figure 16, le dispositif de quadrillage "structuré" 68 reçoit des données du modélisateur de propriétés 66b. Les données de sortie 69 de la figure et de la figure 16 comportant des cartes et des traces de failles ainsi que des trajectoires de puits et des enregistrements SCAL, sont fournis comme entrée dans le dispositif de stockage de données de réservoir 64; et, en réponse à cela, le modélisateur de propriétés 66b génère un "modèle structurel", le "modèle structurel" étant entré dans le dispositif de quadrillage "structuré" 68 qui comporte le "dispositif de quadrillage de surface structuré" de la présente invention (également connu sous le nom de "dispositif de quadrillage de point de coin"). SR 15546 US/RS Dans la figure 16, le logiciel de quadrillage "structuré" 68 comporte un "dispositif de quadrillage de surface structuré" selon la présente invention, le "dispositif de quadrillage de surface structuré" représentant la nouvelle caractéristique du programme de quadrillage de simulation "Flogrid" 60a. Le logiciel de quadrillage structuré 68 est décrit plus en détail ci-dessous et est également décrit plus en détail dans la "description détaillée du mode de réalisation préféré" de la présente demande. Toutefois, en général, le dispositif de quadrillage structuré 68 de la présente invention est utilisé pour élaborer des quadrillages d'écoulement de fluide centrés sur des blocs et de point de coin structuré qui respectent les caractéristiques choisies d'un modèle de propriétés donné. Deux approches sont possibles: quadrillage et quadrillage progressif. Dans l'approche du quadrillage, l'utilisateur peut spécifier une nouvelle limite de quadrillage et contrôler la qualité et la taille du quadrillage vertical et "de surface". Une caractéristique importante du dispositif de quadrillage structuré 68 est la possibilité d'aligner i lignes (de contrôle) de quadrillage et j lignes (de contrôle) de quadrillage avec des failles choisies situées à l'intérieur d'une limite et de faire zigzaguer les failles restantes. Les paramètres de contrôle de quadrillage comportent une répartition globale Nx, Ny et Nz, la perpendicularité, l'uniformité et une distribution Nx et Ny des sous-quadrillages pour contrôler la répartition des lignes et des colonnes i et j dans différentes zones du quadrillage. On peut apporter une amélioration du quadrillage local SR 15546 US/RS cartésien. L'approche de quadrillage progressif (utilisant le dispositif de quadrillage progressif 84 mentionné ci-dessus) permet de construire des quadrillages de simulation en combinant des lignes, des colonnes et des couches d'un modèle géologique dans des cellules de simulation grossières. Le quadrillage progressif permet également à des caractéristiques géologiques significatives tels qu'une grande perméabilité des couches ou des contrastes importants dans les propriétés entre les couches géologiques, d'être saisies avec plus de précision dans le quadrillage de simulation obtenu. On peut utiliser des techniques de quadrillage manuel ou algorithmique (de variabilité et fondées sur l'écoulement). Sur la figure 16, le dispositif de quadrillage structuré 68 fournit ses résultats de "quadrillage de simulation structuré" à un dispositif de mise à une plus grande échelle 70. Le "quadrillage de simulation structuré" comporte une pluralité de "blocs de quadrillage". Le dispositif de mise à une plus grande échelle 70 échantillonne "des propriétés spécifiées" à l'intérieur de chaque "bloc de quadrillage" du "quadrillage de simulation structuré" et effectue un moyennage approprié pour attribuer une propriété de bloc de quadrillage représentatif à chaque "bloc de quadrillage" du "quadrillage de simulation structuré" qui est généré par le dispositif de quadrillage structuré 68. Le dispositif de mise à une plus grande échelle 70 sert à exécuter un "éventail de procédés de divisions progressives", chaque procédé servant à échantillonner les "propriétés spécifiées" à l'intérieur de chaque "bloc de quadrillage" du "quadrillage de simulation structuré" pour attribuer SR 15546 US/RS une "propriété de bloc de quadrillage" représentative à chaque "bloc de quadrillage" du "quadrillage de simulation structuré" généré par le dispositif de quadrillage structuré 68. Les "propriétés spécifiées" qui sont échantillonnées et, en conséquence, la "propriété de bloc de quadrillage" qui est attribuée à chaque "bloc de quadrillage" du "quadrillage de simulation structuré" comportent les propriétés suivantes: porosité, perméabilité absolue, perméabilité relative et pression capillaire. Le dispositif de mise à une plus grande échelle 70 comporte des instruments permettant de faire la moyenne de n'importe quelle valeur réelle. Le dispositif de mise à une plus grande échelle 70 effectue les divisions progressives en utilisant les procédés de phase simple à base d'écoulement et algébrique, des résultats d'importation et des propriétés d'édition. Plus particulièrement, le module du dispositif de mise à une plus grande échelle 70 permet aux utilisateurs de générer des propriétés de bloc de quadrillage de simulation pour chaque "bloc de quadrillage" des quadrillages structurés en échantillonnant et en faisant la moyenne des propriétés géologiques stockées dans un modèle de propriété. "L'éventail du procédé de division progressive" cité ci-dessus comporte les procédés suivants: (1) moyennage algébrique de la porosité, (2) moyennage géométrique et harmonique de Kx, Ky, Kz, (3) choix des types de roche, (4) division progressive fondée sur l'écoulement de Kx, Ky, Kz, et (5) division progressive algébrique simple des autres types de propriétés. En conséquence, sur la figure 16, un "quadrillage de simulation structuré et des propriétés associées à SR 15546 US/RS celui-ci" 40clD est généré par le dispositif de quadrillage non structuré (Petragrid) 40cl (rappelons que Petragrid génère un quadrillage non structuré mais il permet de convertir le dispositif de quadrillage non structuré en un dispositif de quadrillage structuré qui est utilisé par le simulation 40c2A). Il faut remarquer le calculateur de propriétés de cellules 40clD sur la figure 13c2 qui calcule certaines propriétés comme la transmissibilité et le volume poreux entre chaque paire de cellules d'un quadrillage non structuré. De plus, sur la figure 16, un "quadrillage de simulation structuré et des propriétés associées à celui-ci" 90 sont générés par le dispositif de mise à une plus grande échelle 70. En réalité, le dispositif de quadrillage structuré 68 génère le "quadrillage de simulation structuré" mais le dispositif de mise à une plus grande échelle 70 échantillonne des "propriétés spécifiées" à l'intérieur de chaque "bloc de quadrillage" de ce "quadrillage de simulation structuré" et effectue le moyennage approprié afin d'attribuer une "propriété de bloc de quadrillage" représentative à chaque "bloc de quadrillage" du "quadrillage de simulation structuré" qui est généré par le dispositif de quadrillage structuré 68. Sur la figure 16, en ce qui concerne les facteurs de connexion de puits 40clR, rappelons de ce qui vient d'être dit ci-dessus que la structure des données de connexion de puits 40clR comporte, pour chaque puits de forage 17 de la figure 1, une liste particulière de toutes les cellules du quadrillage non structuré 40clM qui sont coupées par et qui coupent le puits de forage 17 de la figure 1 et une liste correspondante de coefficients de débit (appelés facteurs de connexion de SR 15546 US/RS puits) qui sont associés, respectivement, à cette liste particulière de cellules. Sur la figure 16, le simulateur "Eclipse" 40c2A (de la figure 13dl) reçoit les "quadrillages de simulation structurés" et les propriétés associées pour chaque "bloc de quadrillage" 90 du dispositif de mise à une plus grande échelle 70 et le simulateur 40c2A reçoit également les "quadrillages de simulation structuré" et les propriétés associées pour chaque cellule 40clD et les facteurs de connexion de puits clR du dispositif de quadrillage non structuré (Petragrid) 40cl décrit dans la demande de Gunasekera, et, en réponse à cela, le simulateur 40c2A génère un ensemble de "résultats de simulation" 92, les résultats de simulation 92 étant fournis à la visionneuse 3D clN, lesdits résultats de simulation 92 étant également renvoyés au dispositif de quadrillage structuré 68 de la présente invention, au dispositif de mise à une plus grande échelle 70 et au dispositif de quadrillage non structuré (Petragrid) 40cl de la demande de Gunasekera afin de manipuler n'importe quel ensemble généré ci-après de résultats de simulation 92. La visionneuse 3D 40clN fournit une visualisation en 3D de toutes les données générées et entrées. Les caractéristiques comprennent: (1) des options spécialisées pour visualiser des cartes, des puits, des modèles géologiques, des quadrillages, des propriétés d'entrée de simulation et des résultats de simulation, (2) des contrôles interactifs pour changer de plan, effectuer une translation, une rotation et appliquer de la lumière, (3) un éditeur en 3D pour éditer des failles, des limites et des propriétés de quadrillage de simulation et (4) une légende de couleur qui SR 15546 US/RS supporte des objets multiples dans la visionneuse en 3D en même temps, etc. Dispositif de quadrillage structuré 68 comportant le quadrillage de surface structuré de la présente invention de la figure 16 En référence aux figures 16a à 23, une description détaillée du dispositif de quadrillage structuré 68 de la figure 16 est donnée dans les paragraphes suivants. * Des modèles de simulation en général sont composés de trois types: des modèles de quadrillage structuré (par exemple le dispositif de quadrillage structuré 68) qui effectue un "quadrillage de point de coin", également appelé "quadrillage de ligne de coordonnées" ou bien "quadrillage de surface", des modèles de quadrillage non structuré (par exemple le dispositif de quadrillage non structuré Petragrid 40cl) qui utilisent des quadrillages en forme au tétraèdre ou PEBI et des modèles structurés-non structurés qui utilisent à la fois des quadrillages en forme de tétraèdre et de point de coin. Le logiciel de quadrillage structuré "Flogrid" 68 de la figure 16 de la présente invention sert à effectuer le nouveau "quadrillage de point de coin" également appelé "quadrillage de ligne de coordonnées" ou bien "quadrillage de surface". Les paragraphes suivants, en référence aux figures 17 à 23 traiteront de la fonction de "quadrillage de surface" ou bien "de ligne de coordonnées" mentionnés ci-dessus de la présente invention qui est effectuée par le "dispositif de quadrillage structuré" 68. La fonction de "quadrillage SR 15546 US/RS de point de coin" ou bien "de quadrillage de surface" est la technologie centrale utilisée dans le "dispositif de quadrillage structuré" 68 de la figure 16 de la présente invention. Quadrillage de point de coin ou quadrillage de surface A ce stade, il est nécessaire de définir un "quadrillage de point de coin" et de souligner les étapes principales qui sont effectuées par le dispositif de quadrillage structuré 68 lors de l'élaboration d'un quadrillage de point de coin. Le "quadrillage de point de coin" comporte les composants principaux suivants: (1) une "ligne de coordonnées", (2), un "tube de coordonnées" qui présente des bords, chaque bord du tube de coordonnées étant composé d'au moins l'une des "lignes de coordonnées" et (3) un "bloc de quadrillage", un bloc de quadrillage étant obtenu en tranchant une paire d'horizons sur le tube de coordonnées. Les "lignes de coordonnées" sont définies comme étant un ensemble de lignes droites tracées vers le bas, à partir des noeuds d'un "quadrillage de surface structuré" de quadrilatères. Sur la figure 16a, un exemple de "quadrillage structuré" (ou quadrilatères) est illustré. Les "lignes de coordonnées" partent des noeuds du quadrillage structuré de la figure 16a et ces lignes de coordonnées entrent "dans la page" de la figure 16a. Les coins des quadrilatères sont situés en haut des lignes de coordonnées sur la figure 16a. Sur la figure 16e une vue en 3D du "tube de coordonnées" est illustrée. Le "tube de coordonnées" SR 15546 US/RS présente des bords et les "lignes de coordonnées" correspondent aux bords du tube de coordonnées, tel qu'illustré sur la figure 16b. Sur la figure 16b, ces "tubes de coordonnée" remplissent la région située à l'intérieur de la "limite définie par l'utilisateur" sans qu'il y ait de trous. Sur la figure 16c un "bloc de quadrillage" est obtenu en tranchant au moins une paire d'horizons à travers le tube de coordonnées. Le "tube de coordonnées" est tranché en "blocs de coordonnées" en plaçant des points sur chaque ligne de coordonnées, le même nombre de points étant placé sur chaque ligne de coordonnées. S'il y a N+l points sur chaque ligne de coordonnée, alors il y a "n" blocs de quadrillage dans le tube de coordonnées. Ceci est illustré sur la figure 16c pour n = 1. Les étapes principales d'élaboration d'un quadrillage de surface structuré ou de point de coin sont traitées dans le paragraphe suivant en référence aux figures 17 à 23 des dessins. Le dispositif de quadrillage "structuré" Flogrid 68 de la figure 14c et de la figure 16 comporte un "dispositif de quadrillage de surface structuré" connu également sous le nom "dispositif de quadrillage de point de coin" selon la présente invention. Le dispositif de quadrillage de surface structuré/dispositif de quadrillage de point de coin élabore un ensemble de lignes de contrôle qui forment un quadrillage de lignes de simulation Nx fois Ny fois 1 à l'intérieur d'une limite définie par l'utilisateur. Les lignes de commande de ce quadrillage de simulation respectent à la fois la "limite définie par l'utilisateur" et les failles i et j choisies à SR 15546 US/RS l'intérieur de la ligne. Les lignes de contrôle qui ne sont pas obligées de se tenir dans une limite ou une faille sont placées en fonction de paramètres de quadrillage contrôlés par l'utilisateur. Sur la figure 17, la construction, par le "dispositif de quadrillage structuré" 68 de la présente invention, d'un "quadrillage de surface structuré" ou bien un "quadrillage de point de coin" implique sept (7) étapes: (étape 1) spécifier une limite (sur un horizon d'une formation terrestre) et choisir des failles (qui se trouvent dans cet horizon et à l'intérieur de cette limite) sur lesquelles effectuer des quadrillages, bloc 96 sur la figure 17, (étape 2) élaborer une triangulation irrégulière interne qui respecte les failles choisies et la limite choisie, bloc 98, (étape 3) construire des champs de vecteur i et j, bloc 100, (étape 4) créer des lignes de commande en extrapolant des failles en utilisant des champs de vecteur, bloc 102, (étape 5) répartir des lignes de coordonnées le long des lignes de contrôle et de la limite, bloc 104, (étape 6) répartir les lignes de coordonnées restantes, bloc 106, et (étape 7) quadriller verticalement, bloc 107. L'utilisateur/opérateur, assis dans le poste de travail 40 des figures 14a et 14b, contrôle et exécute les étapes suivantes à partir de l'organigramme de la figure 17: étape 1 (spécifier une limite et choisir des failles sur lesquelles effectuer un quadrillage, bloc 96), étape 3 (construire des champs de vecteur i et j, bloc 100), étape 5 (répartir des lignes de coordonnées le long des lignes de commande et de la limite, bloc 104), et étape 7 (quadriller verticalement, bloc 107). SR 15546 US/RS Chacune de ces étapes (1) à (7), blocs 96 à 107 sur la figure 17, seront traitées plus en détail ci- dessous en référence aux figures 18 à 23 des dessins. Etape 1 - Spécifier une limite et choisir des failles sur lesquelles effectuer un quadrillage, bloc 96 de la figure 17 Supposons qu'un horizon d'une formation terrestre présente trois failles qui traversent l'horizon. Supposons en outre qu'une limite délimitée et approximativement rectangulaire et tracée sur cet horizon, la limite englobant les trois failles traversant l'horizon. Ensuite, supposons que les trois failles peuvent être classées en deux failles i et une faille j traversant l'horizon. La figure 18 illustre cette limite qui renferme les deux failles i et la faille j. Sur la figure 18, l'utilisateur doit définir cette limite mais la limite doit être soit un rectangle soit un rectangle déformé (un dispositif de quadrillage "structuré" doit avoir une limite en forme de rectangle ou en forme de rectangle déformé). Par conséquent, sur la figure 18, la "limite" est composée de deux "limites i" et deux "limites j". Une limite correcte doit suivre la tendance des failles particulières que l'utilisateur peut souhaiter respecter explicitement dans le quadrillage de simulation. Ceci est illustré sur la figure 18. Par conséquent, les limites i sont indiquées par le chiffre 108 sur la figure 18 et les limites j sont indiquées par le chiffre 110 sur la figure 18. A l'intérieur des limites 108, 110, il y a trois failles: une faille j 112 (une faille j étant dirigée dans le sens gauche- SR 15546 US/RS droite), et deux failles i 113 (une faille i étant dirigée dans le sens haut-bas). La limite 110 de la figure 18 suit la tendance de la faille j 112 et la limite 108 suit la tendance des failles i 113. Les failles sont classées par l'utilisateur et peuvent être des failles i, des failles j ou des failles en zigzag. Les failles i sont des failles qui présentent une valeur constante i dans le quadrillage de simulation généré (par exemple une faille peut se trouver dans la ligne i=3). Les failles j peuvent se trouver sur des lignes de quadrillage présentant une constante j). Les failles en zigzag ne sont alignées avec aucune limite ou bien sont trop distordues. Etape 2 - Elaborer une triangulation irrégulière interne qui respecte les failles choisies et la limite choisie, bloc 98 Le "dispositif de quadrillage de surface structuré" Flogrid placé à l'intérieur du dispositif de quadrillage structuré 68 élabore un quadrillage irrégulier de triangles qui respecte la limite et toutes les failles choisies (se référer à la figure 44 et au traitement détaillé en référence à la figure 44 mentionnée dans la description détaillée du mode de réalisation préféré ci-dessous). Ceci est effectué en modifiant séquentiellement une triangulation de départ et en ajoutant chaque segment de ligne provenant de la limite ou des failles tour à tour. Cette étape (bloc 98 de la figure 17) ne peut pas être vue par l'utilisateur/opérateur. Etape 3 - Construire des champs de vecteur i et j, bloc 100 SR 15546 US/RS Le "dispositif de quadrillage de surface structuré) 68b du dispositif de quadrillage structuré 68 de la figure 16 construit deux champs de vecteur (comme des champs magnétiques ou de vitesse), les vecteurs étant calculés sur chaque noeud de la triangulation (qui a été élaborée au bloc 98 de la figure 17 - voir figures 45bl, 45b2 et 49 de la "description détaillée du mode de réalisation préféré"). Les noeuds qui se trouvent sur la limite définie par l'utilisateur ou bien sur une faille présentent des vecteurs qui doivent respecter les tangentes aux failles et à la limite. Les vecteurs situés sur les noeuds restants sont interpolés. C'est-à-dire qu'une "pluralité de noeuds" de la triangulation (élaborée par l'intermédiaire du bloc 98) se trouve en réalité sur la partie supérieure des limites 108, 110 ou des failles 112, 113 sur la figure 18 et un "vecteur", appelé ci-après un "vecteur gelé" est connecté à chaque noeud de la "pluralité de noeuds". Chacun des "vecteurs gelés" présente un sens établi (ceci expliquant la raison de l'adjectif "gelé"). Le sens établi de chaque vecteur gelé est un sens qui est parallèle aux tangentes aux failles et/ou aux limites o les failles et les limites coupent réellement le noeud de la "pluralité de noeuds". Les vecteurs situés sur les noeuds restants (qui ne se trouvent pas sur les limites 108, 110 ou les failles 112, 113) sont appelés ci-après "vecteurs non gelés" puisque le sens de ces vecteurs n'est pas connu; par conséquent, le sens des vecteurs non gelés doit être déterminé par une technique d'interpolation, l'interpolation étant effectuée par rapport au sens des "vecteurs gelés". En conséquence, un ensemble de lignes SR 15546 US/RS de champ suit les lignes i et un ensemble de lignes de champ suit les lignes j. L'interpolation, qui permet de déterminer le sens des vecteurs "non gelés" à partir de la direction des vecteurs "gelés", peut être soit une interpolation "isotrope" soit une interpolation "anisotrope". Il y a également un contrôle d'orthogonalité qui influence les champs de vecteurs de sorte que les champs de vecteurs sont mutuellement orthogonaux. L'interpolation isotrope réduit l'influence de courbures importantes dans les failles ou sur les limites. Une interpolation anisotrope provoque la propagation de ces discontinuités. L'effet de ce choix sera illustré ci-après dans cette demande. Etape 4 - Créer des lignes de commande en extrapolant des failles en utilisant des champs de vecteur, bloc 102 de la figure 17 La figure 19 illustre des extensions uniformes de failles jusqu'à la limite pour élaborer des "lignes de contrôle". Par exemple, un ensemble d'extensions (également connu sous le nom de "courbes entières") , 117 crée une interconnexion entre la faille j 112 et la limite 108, ce qui permet de créer une "ligne de commande" [115, 112, 117]; un autre ensemble d'extensions (courbes entières) 119, 121 crée une interconnexion entre une faille i 113 et la limite 110, ce qui permet de créer une autre "ligne de contrôle" [119, 113, 121]; et un autre ensemble d'extensions (courbes entières) 123, 127 crée une interconnexion entre la faille i 113 et limite 110, ce qui permet de créer une autre "ligne de contrôle" [123, 113, 125]. Sur la figure 19, les "courbes entières" ou "extensions" [115, 117, 119, 121, 123, 125] sont SR 15546 US/RS construites loin des deux extrémités de chaque faille [112, 113]. Les courbes entières ressemblent à des lignes de flux dans un écoulement de fluide ou des lignes de champ de flux dans des systèmes électromagnétiques. Ces courbes entières et les failles interconnectées associées forment des "lignes de contrôle" qui (si l'utilisateur a effectué une classification correcte des failles et si une limite appropriée a été fournie) fractionnent la zone à l'intérieur de la limite définie par l'utilisateur en régions (ou sous-quadrillages) qui peuvent être déformées en continu dans un quadrillage "structuré" rectangulaire. Etape 5 - Distribuer des lignes de coordonnées le long des lignes de contrôle et de la limite, bloc 104 La figure 20 illustre les lignes de contrôle [115, 112, 117], [119, 113, 121], et [123, 113, 125], certaines "lignes supplémentaires" 127 (mais qui ne sont pas uniformes) étant disposées de façon adjacente aux lignes de contrôle. La figure 21 illustre également la même ligne de contrôle avec ces lignes supplémentaires 127 (mais ayant une uniformité finie). L'avant-dernière fonction, effectuée au cours de la construction d'un quadrillage structuré, implique la répartition de "points" le long des lignes de contrôle, entre les points o les lignes i coupent les lignes j. Le nombre de "points" à répartir est contrôlé par l'utilisateur de la façon suivante: (1) globalement, Nx et Ny sont spécifiés et Flogrid 60a décide de la répartition des points ou bien (2) localement, en spécifiant Nx entre une paire quelconque de lignes i ou Ny entre une paire quelconque de lignes j. Le paramètre SR 15546 US/RS d'uniformité contrôle la façon dont cela est effectué. Lorsqu'ils sont fixés à zéro, les points sont placés entre les intersections de façon équidistante. Si l'uniformité est augmentée jusqu'à un maximum de un, les points sont répartis d'une façon plus uniforme de sorte que la distance entre les paires successives ne change que progressivement. Ceci est illustré sur les figures 20 et 21 o l'une des lignes 127a a été déplacée de façon significative vers la droite puisque le sous-quadrillage adjacent présente une densité supérieure de lignes de quadrillage. Etape 6 - Répartir les lignes de coordonnées restantes, bloc 106 Le quadrillage final des lignes de coordonnées est construit en interpolant les points des limites de chaque sous-quadrillage en utilisant la méthoded'interpolation transfinie. Les lignes de coordonnées élaborées par les "dispositifs de quadrillage de surface structures" sont droites et verticales. Ainsi, les lignes de coordonnées sont complètement spécifiées par un quadrillage de surface 2D. Les lignes de coordonnées utilisées dans le quadrillage de point de coin sont obtenues en réduisant les lignes droites, verticalement à partir de noeuds du quadrillage de surface. Etape 7 - Quadriller verticalement, bloc 107 de la figure 17 Dans le traitement suivant de la présente invention, le terme "bloc de faille" est utilisé (pour être plus correct, il faudrait l'appeler "bloc" puisqu'il y a une différence entre un "bloc" du corps SR 15546 US/RS structurel et un "bloc de quadrillage" dans un quadrillage de simulation). Pour chaque "tube de coordonnées", le dispositif de quadrillage structuré Flogrid 60a 68 attribue ce "tube de coordonnées" à un "bloc de faille", pourvu que le point central du haut du tube de coordonnées occupe ce bloc de faille particulier. A ce stade, le terme "couche" dans une formation terrestre est défini par deux horizons (pas nécessairement adjacents) de formation terrestre. Au cours de cette étape 7, l'utilisateur peut spécifier le nombre de "couches de simulation" qui doivent être placées à l'intérieur de chaque "couche". Les horizons situés à l'intérieur du bloc de faille sont utilisés pour trancher le tube de coordonnées en "bloc de quadrillage". Ceci est effectué en plaçant les intersections des "lignes de coordonnées" avec chaque horizon de la formation, ce qui permet de produire quatre points d'intersection entre les quatre lignes de coordonnées et un horizon unique dans la formation terrestre, les quatre points étant utilisés pour définir quatre des points de coin du "bloc de quadrillage". Si une couche géologique doit devenir une "couche de simulation" unique, alors deux horizons adjacents de la formation terrestre, en séquence, définissent les huit points de coin d'un bloc de quadrillage de simulation unique. Si la "couche" doit être fractionnée en un "nombre donné de couches de simulation", et par conséquent en un "nombre de données de blocs de quadrillage", rappelons qu'une ligne de coordonnées particulière comprend une "pluralité de lignes" interposées, respectivement, entre une "pluralité de coins consécutifs" sur cette ligne de SR 15546 US/RS coordonnées particulière, lorsque la "pluralité de lignes" sur cette ligne de coordonnée particulière dans la "couche" est égal au "nombre donné de couches de simulation", alors la "couche" comporte le "nombre donné de couches de simulation" et le "nombre donné de couches de simulation" dans la "couche" comporte le "nombre donné de blocs de quadrillage". Les paragraphes suivants traitent des paramètres spécifiés par l'utilisateur et des choix faits par l'utilisateur suivants: orthogonalité, uniformité et quadrillage isotrope et anisotrope. L'orthogonalité - L'orthogonalité est un poids compris entre 0 et 1 qui influence les champs de vecteurs et par conséquent les lignes de contrôle pour qu'elles soient mutuellement orthogonales. Si le poids est égal à 0, les lignes suivent simplement les limites et leurs failles respectives. Si le poids est établi égal à 1, les lignes de champs tentent également de se couper à angle droit. Si la limite ou les failles ne sont pas mutuellement orthogonales, alors les vecteurs ne seront pas exactement orthogonaux. Les meilleures valeurs d'orthogonalité se trouvent quelque part entre ces valeurs. Une petite expérimentation peut être nécessaire de sorte que l'on puisse former un quadrillage répondant aux préférences de l'utilisateur. Si la limite ou les failles contiennent des coins très pointus, alors l'application d'une certaine orthogonalité est essentielle si des quadrillages inversés doivent être évités. Parfois, une faille peut ne pas être caractérisée comme une faille i ou j et est traitée de préférence comme un zigzag. De la même manière, si la limite est trop compliquée, elle peut mener à des quadrillages inversés. Le dispositif de SR 15546 US/RS quadrillage structuré 68 est fondé sur l'idée selon laquelle le quadrillage est un quadrillage rectangulaire distordu. Si les lignes cibles sont trop loin d'un quadrillage rectangulaire, quelque chose se cassera. Le logiciel avertit l'utilisateur de cela, et il doit ensuite se charger du réglage de la limite ou du changement de la classification des failles. Uniformité - L'uniformité du quadrillage de simulation est contrôlée par l'espacement des points requis par l'utilisateur sur les lignes de contrôle. Si l'uniformité est établie égale à 0, alors la distance qui sépare chaque paire de points le long de la ligne de contrôle entre deux lignes de contrôle sera la même. Ainsi, la distance qui sépare des paires de points peut présenter un saut distinct lorsqu'on se déplace entre les paires de lignes de contrôle. Si l'uniformité est établie à 1.0, alors la distance qui sépare les points change de façon uniforme d'une paire de lignes de commande à une autre. Ce contrôle est atteint grâce à un algorithme de spline de tension non linéaire. Reportons-nous maintenant aux figures 22 et 23. Quadrillage isotrope et anisotrope - Sur les figures 22 et 23, un exemple d'interpolation isotrope entre deux failles est illustré sur la figure 22 et un exemple d'interpolation anisotrope entre deux failles est illustré sur la figure 23. Les options principales du "quadrillage de surface" sont: (1) interpolation isotrope des champs de vecteurs (c'est-à-dire "quadrillage isotrope") et interpolation anisotrope des champs de vecteurs (c'est-à-dire "quadrillage anisotrope"). Dans une interpolation "isotrope", le dispositif de quadrillage structuré 68 établit le vecteur i "non SR 15546 US/RS gelé" sur le coin d'un triangle à: (1) la moyenne pondérée du vecteur orthogonal au vecteur j au niveau de ce coin et (2) la moyenne arithmétique des vecteurs i au niveau des coins de n'importe quel triangle partageant le premier coin. Ceci nécessite un procédé itératif et constitue la cause principale du laps de temps requis par le dispositif de quadrillage pour accomplir sa tâche. Le champ j est construit d'une façon similaire. Dans une interpolation "anisotrope", la moyenne arithmétique est remplacée par une moyenne pondérée de distance. Ceci prend plus de temps à calculer, ce qui explique le temps d'attente plus important pour ce type de quadrillage. Sur les figures 22 et 23, la différence entre les deux types d'interpolation (isotrope et anisotrope) est mieux illustrée par le résultat obtenu en réponse à la performance de l'interpolation isotrope par rapport à la performance de l'interpolation anisotrope. Dans la plupart des cas, l'interpolation isotrope et anisotrope produit des résultats similaires. Toutefois, dans certains cas, tel que cela est illustré sur les figures 22 et 23, les quadrillages sont quelque peu différents. Sur les figures 22 et 23, les lignes fines représentent les failles originales et les lignes en gras représentent les extensions dans des lignes de contrôle. En ce qui concerne la performance et la qualité, le dispositif de quadrillage structuré 68 de la présente invention fonctionne plus rapidement si on utilise l'interpolation isotrope sans orthogonalité. Le dispositif de quadrillage structuré 68 produit le quadrillage de meilleure qualité avec des points SR 15546 US/RS d'uniformité et d'orthogonalité élevés. Toutefois, la qualité est un concept subjectif et, par conséquent, une certaine expérimentation est suggérée. DESCRIPTION DETAILLEE DU MODE DE REALISATION PREFERE Une description détaillée du programme de quadrillage de simulation Flogrid 60a des figures 14c et 16, comportant une description détaillée du dispositif de quadrillage structuré 68 des figures 14c et 16 qui comporte le "dispositif de quadrillage de surface structuré" (connu également sous le nom de "dispositif de quadrillage de point de coin") de la présente invention et le "dispositif de quadrillage de bloc" va être énoncée dans les paragraphes suivants en référence aux figures 24 à 52 des dessins. Le programme de quadrillage de simulation "Flogrid" 60a En référence à la figure 24, une construction détaillée d'un système logiciel fondé sur un poste de travail est illustrée et ce système est conçu pour être stocké dans la mémoire du poste de travail 40c de la figure 14b, le système fondé sur le poste de travail de la figure 24 comportant le programme de quadrillage de simulation "Flogrid" 60a qui réagit aux données d'entrée 69 des figures 15 et 16 pour générer le modèle de simulation 76, un simulateur de logiciel (appelé "Eclipse") 40c2 réagissant au modèle de simulation 76 généré par le programme de quadrillage de simulation "Flogrid" 60a pour générer un ensemble de "résultats de simulation", et une visionneuse 3D 40clN conçue pour SR 15546 US/RS afficher les "résultats de simulation" générés par le simulateur "Eclipse" 40c2. Sur la figure 24, le programme de quadrillage de simulation Flogrid 60a reçoit les données d'entrée 69 des figures 15 et 16, les données d'entrée 69 comportant les données suivantes: données d'horizons dl, représentant une pluralité d'horizons dans une formation terrestre: données de failles 30d2 représentant une ou plusieurs failles traversant et coupant les horizons dans la formation terrestre; et une pluralité de modèles géocellulaires 147. Les données d'horizon 30dl et les données de failles 30d2 peuvent être mieux illustrées sur le dessin de la figure 24a. Sur la figure 24a, une pluralité d'horizons 129 dans une formation terrestre coupe une faille 131 qui traverse et coupe chacun des horizons 129. Les données d'horizons 30dl et les données de failles 30d2 sont recueillies de la manière suivante. Sur la figure 24a, une source 133 (telle qu'une source explosive) produit des vibrations sonores 135. Ces vibrations sonores 135 se propagent à l'intérieur d'une formation terrestre 153 et se réfléchissent à l'extérieur des horizons 129 dans la formation terrestre 153. Les vibrations sonores 135 se propagent vers le haut à partir des couches d'horizons 129 jusqu'à ce que les vibrations sonores soient reçues dans une pluralité de récepteurs géophones 137. Les récepteurs géophones 137 produisent un signal de sortie représentant les "données reçues" 139, les données reçues 139 étant entrées dans un camion d'enregistrement 141. Un ordinateur de camion d'enregistrement 143 reçoit les "données reçues" 139 et, en réponse à cela, il génère un "enregistrement de SR 15546 US/RS sortie de données sismiques" 145 qui contient des données sismiques représentant les horizons 129 et les failles 131 dans la formation terrestre. Les données d'horizons 30dl et les données de failles 30d2 de la figure 24 résultent des données sismiques stockées et contenues dans "l'enregistrement de sortie de données sismiques" 145. Sur la figure 24, le programme de quadrillage de simulation "Flogrid" 60a comporte les éléments structurels suivants. Un dispositif de stockage de données de réservoir 64 est un moyen de stockage conçu pour stocker provisoirement les données d'horizons 30dl, les données de failles 30d2 et les modèles géocellulaires 147. Un canevas structurel 149 reçoit les données stockées dans le dispositif de stockage de données de réservoir 64, élabore un "modèle géologique" et fait passer ce "modèle géologique" au modèle de propriété 66b. Le modèle de propriété 66b reçoit le "modèle géologique" qui provient du canevas structurel 149 et génère un "modèle structurel" qui sera envoyé au dispositif de quadrillage structuré 68. En utilisant un "dispositif de quadrillage progressif" 84, l'utilisateur peut fournir certaines entrées au modèle de propriété 66b pour simplifier le "modèle géologique" qui a été reçu dans le modèle de propriété 66b et pour regrouper certaines des couches du "modèle géologique" afin de générer le "modèle structurel". Le dispositif de quadrillage structuré 68 reçoit le "modèle structurel" provenant du modèle de propriété 66b et, en réponse à la réception du "modèle structurel", il élabore un "quadrillage de simulation structuré" dans la formation terrestre. C'est-à-dire que le dispositif de SR 15546 US/RS quadrillage structuré 68 de la présente invention comporte (1) un "dispositif de quadrillage de surface structuré" et (2) un "dispositif de quadrillage de bloc". Le "quadrillage de surface structuré" élabore un "quadrillage de surface structuré" sur un horizon faillé dans la formation terrestre et le "dispositif de quadrillage de bloc" réduit les lignes de coordonnées à partir des noeuds du "quadrillage de surface" afin de produire une pluralité de "blocs de quadrillage" d'un "quadrillage structuré" dans la formation terrestre. Le dispositif de quadrillage progressif 70 de la figure 24 attribue des propriétés (par exemple porosité, perméabilité, etc.) à chaque cellule du "dispositif de quadrillage de simulation structuré" généré par le dispositif de quadrillage structuré 68 et fait passer une "sortie de quadrillage de simulation" 90 comportant le quadrillage structuré et les propriétés associées à chaque cellule de celles-ci, au simulateur Eclipse c2. Dans le même temps, le dispositif de quadrillage non structuré Petragrid 40cl, (traité dans la première section de la description du mode de réalisation préféré de la présente demande) reçoit le "modèle structurel" provenant du modèle de propriétés 66b, et en réponse à cela, il génère sa propre sortie de quadrillage de simulation structuré 44clD que l'on fait passer au simulateur Eclipse 40c2. Le simulateur d'Eclipse 40c2 ne fait pas partie du programme de quadrillage de simulation Flogrid 60a mais sert à générer son propre ensemble de "résultats de simulation" qui sont affichés sur la visionneuse 3D clN. Se reporter à la figure 13d2 pour voir un exemple des "résultats de simulation" générés par le SR 15546 US/RS simulateur Eclipse 40c2 et affichés sur la visionneuse D 40clN. En référence à la figure 25, une construction détaillée du canevas structurel 149 de la figure 24 est illustrée. Sur la figure 25, le canevas structurel comprend le "dispositif d'élaboration du modèle interne" 40clF et le "modèle structurel" 66a. Rappelons, d'après la première section de la description du mode de réalisation préféré de la présente demande, lorsque le logiciel de quadrillage non structuré Petragrid 40cl a été traité, que le "dispositif d'élaboration du modèle interne" 40clF reçoit l'enregistrement de sortie de diagraphie 25 et l'ensemble de données sismiques réduites 30d et, en réponse à cela, le dispositif d'élaboration du modèle interne 40clF génère une "structure de données de modèle interne 40clG". Les dispositifs d'élaboration du modèle interne 40clF créent une "limite" dans laquelle les "points" de la structure de données "de points 40clJ" (générée par le "distributeur de points" 40cllI) vont être répartis. L'enregistrement de sortie de données sismiques réduites 30d définit l'endroit o les horizons 13 de la figure 1 et les failles sont situées à l'intérieur de la "limite", et l'enregistrement de sortie de diagraphie 25 définit l'endroit o les puits de forage 17 de la figure 17 sont situés à l'intérieur de la "limite". Toutefois, l'opérateur, au niveau du poste de travail 40 de la figure 9, définit la "limite" elle-même en définissant une ligne dans l'espace aérien. La structure de donnée du "modèle interne" clG comprend cette "limite". Rappelons de nouveau, d'après la description du mode de réalisation préféré de la présente demande, que SR 15546 US/RS le modèle structurel 66a est décrit dans les deux demandes de brevets américains en suspens suivantes: (1) demande américaine en suspens antérieure numéro de série 08/823.107 de Abbott (dénommée ci-après "la demande de Abbott"), déposée le 24 mars 1997 et intitulée "Method and Apparatus for Determining Geologic Relationships for Intersecting Faults"; et (2) demande provisoire américaine en suspens antérieure numéro de série 60/029.524, déposée le 31 octobre 1996, intitulée "Method for 3D Modeling of Faulted Geologic Horizons", et la demande américaine antérieure en suspens numéro de série 08/916.841 attribuée à Graf et al (dénommée ci-après la "demande de Graf"), déposée le 22 août 1997 et intitulée "Automatic Non-Artificially Extended Fault Surface Based Horizon Modeling System", correspond au numéro 94. 0008. La demande de Graf génère, en tant que produit final, un "modèle d'horizon faillé final". Par conséquent, sur les figures 25 et 25a, le modèle structurel 66a de la figure 25 génère le "modèle d'horizon faillé final" 151. Un exemple de "modèle d'horizon faillé final" 151 généré par le modèle structurel 66a est illustré sur la figure 25a. Sur la figure 25a, le "modèle d'horizon faillé final" 151 est une représentation tridimensionnelle de la formation terrestre 153 de la figure 24a qui a été soumise à une opération sismique traitée ci-dessus conjointement avec la figure 24a, et comporte une pluralité d'horizons 155 qui coupent une pluralité de failles 157. En référence à la figure 26, le dispositif de quadrillage 68 de la figure 24 comporte un "éditeur de limite" 68a; un "dispositif de quadrillage de surface structuré" (ou bien "dispositif de quadrillage de point SR 15546 US/RS de coin") 68b qui produit le "quadrillage de surface structuré" sur une limite spécifiée par l'utilisateur sur un horizon faillé; et un "dispositif de quadrillage de bloc" 68c qui réduit les lignes de coordonnées à partir des noeuds du "quadrillage de surface structuré" produit par le "quadrillage structuré" produisant ainsi le "quadrillage de simulation structuré" composé d'une pluralité de "blocs de quadrillage" (présentant une coupe transversale approximativement rectangulaire) qui sont entrés dans le dispositif de mise à une plus grande échelle 70. Les fonctions effectuées par chacun de ces éléments apparaîtront clairement à la lecture de la description suivante. En référence à la figure 27, le "dispositif de quadrillage de surface structuré" (ou bien "dispositif de quadrillage de point de coin") 68b de la figure 26 produit un "quadrillage de surface structuré" qui est disposé sur toute une limite spécifiée par l'utilisateur sur la surface supérieure d'un horizon faillé dans une formation terrestre. Le "dispositif de quadrillage de surface structuré" produit le "quadrillage de surface structuré" sur toute une limite sur ladite partie supérieure de l'horizon faillé dans la formation terrestre en exécutant les étapes fonctionnelles suivantes: (1) former des triangles d'une triangulation, bloc 68bl, (2) élaborer des vecteurs directionnels sur les noeuds de la triangulation, bloc 68b2, (3) en utilisant le sens des vecteurs directionnels, élaborer des lignes de flux (qui comportent des "lignes de contrôle" et "lignes supplémentaires") sur toute la limite spécifiée par l'utilisateur, bloc 68b3, (4) répartir à distance égale SR 15546 US/RS les lignes de flux, bloc 68b4, ce qui permet de produire un "quadrillage de surface structuré"; et (5) trouver des "intersections entre les lignes de flux", bloc 68b5. Le "quadrillage de surface structuré", généré par le dispositif de quadrillage de surface structuré 68b est reçu par le "dispositif de quadrillage de bloc" 68c. Le dispositif de quadrillage de bloc 68c réduit les lignes de coordonnées à partir des "intersections entre les lignes de flux" du quadrillage de surface structuré; et, lorsque les lignes de coordonnées coupent les lignes orientées situées en bas dans la formation terrestre, une multitude de "blocs de quadrillage", définis par leurs lignes de coordonnées et les deux horizons (le plus haut et le plus bas) sont formés, la multitude de "blocs de quadrillage" constituant ensemble un "quadrillage structure". Les éléments de Flogrid 68a traités ci-dessus en référence aux figures 24, 25, 26 et 27 seront traités plus en détail ci-dessous en référence aux figures 28a à 52. Sur les figures 28a à 28b, les données d'entrée 69 qui circulent dans le programme de quadrillage de simulation "Flogrid" 60a de la figure 24 (ci-après appelé "Flogrid 60a) comportent deux pistes de données: les données incluses dans la "piste 1" illustrée sur la figure 28a et les données incluses dans la "piste 2" illustrée sur la figure 28b. Dans la "piste 1" de la figure 28a, deux types de données d'entrée 69 sont entrés dans Flogrid 60a: un ensemble de "données de cartes" (c'est- à-dire les données d'horizons 30dl de la figure 24), et des "traces de failles" (c'est-à-dire les données de SR 15546 US/RS failles 30d2 de la figure 24). Pour la piste 1, l'utilisateur/opérateur doit choisir les données de cartes et les traces de failles 30bl, 30b2 pour élaborer un "modèle structurel". Toutefois, dans la "piste 2", de la figure 28b, les données entrées 69 comportent également un ensemble de données appelé "modèles géologiques de cellules fines" 147 (c'est-à-dire, les modèles géocellulaires 147 de la figure 24). Les "modèles géologiques de cellules fines" 147 de la figure 28b sont également entrés dans Flogrid 60a. Puisqu'un "modèle structurel" comprend un "modèle géologique" plus un "modèle de propriété", le "modèle géologique" est encore présent dans le "modèle géologique de cellules fines" de la piste 2. Par conséquent, pour la piste 2, il n'est pas nécessaire que l'utilisateur choisisse des données de cartes et des traces de failles (comme pour la piste 1) pour élaborer le modèle géologique. Par conséquent, les "données entrées" 69, entrées dans Flogrid 60a, comportent les "données (d'horizon) de cartes" et les "traces de failles" pour la piste 1, tel que représenté sur la figure 28a et les "modèles géologiques de cellules fines" pour la piste 2, tel que représenté sur la figure 28b. Le dispositif de stockage de données de réservoir 64 et l'élaboration du modèle structurel 66a de l'ossature structurelle le 149 de la figure 25 Sur les figures 24, 28a et 28b, les données entrées (c'est-à-dire les données d'entrée de la piste 1 de la figure 28a et les données d'entrée de la piste 2 de la figure 28b) sont chargées dans le dispositif de stockage de données de réservoir (RDS) 64 de Flogrid SR 15546 US/RS a. Les "données d'entrée" pour la piste 1 peuvent être composées de deux "cartes", chaque "carte" pouvant être une version différente du même horizon. L'utilisateur/opérateur, assis à son poste de travail 40 de la figure 14b, choisit les données, dans les "données entrées" des figures 28a et 28b, qui seront utilisées pour élaborer le modèle structurel 66a de la figure 25. Si l'utilisateur choisit les "données d'entrée" de la piste 2 de la figure 28b, les données d'entrée de la piste 2 du "modèle géologique de cellules fines" de la figure 28b comportent déjà le "modèle géologique". Rappelons qu'un "modèle structurel" 66a est réellement composé d'un modèle géologique plus un modèle de propriétés 66b. Par conséquent, dans la piste 2, nous disposons déjà d'un modèle géologique inclus dans le "modèle géologique de cellules fines" de la figure 28b. Le modèle structurel 66a définit l'endroit o les failles sont situées dans une formation terrestre, la façon dont les failles correspondent à des surfaces, illustre les horizons et définit les surfaces des horizons (c'est-à-dire la "géométrie" ou positions des surfaces dans la formation, la "topologie" de la formation ou bien le modèle et les réseaux contenus dans la formation). Un bon exemple d'un modèle structurel 66a est illustré sur la figure 25a. Par ailleurs, le modèle de propriétés 66b applique certaines propriétés à la roche de la formation terrestre, ces propriétés comportent la porosité et la perméabilité. Toutefois, dans la piste 1, o les données de cartes et les traces de failles sont entrées dans Flogrid 60a, Flogrid 60a élabore réellement un modèle structurel 66a à partir des SR 15546 US/RS données de cartes et des données de traces de failles. Les données entrées pour la piste 1 (données de cartes et données de traces de failles) sont supposées être des données correctes. Par conséquent, dans la piste 1, l'utilisateur choisit les "cartes" et les "traces de failles" des "données d'entrée" de la piste 1 (stockées dans le dispositif de stockage de données de réservoir 64), et les "cartes " et les "traces de failles" de la piste 1 (figure 28a) seront utilisées pour élaborer le modèle structurel 66a. Sur la figure 29, par exemple, l'opérateur assis au poste de travail 40 de la figure 14a élabore le modèle structurel 66a en choisissant une "carte 1" 120, une "carte 2" 122 et une "carte 3 124 et l'utilisateur/opérateur choisit alors les failles 126 qui coupent les cartes (ou horizons) 120, 122, 124. A ce stade, le modèle structurel 66a pour la piste 1 (figure 28a) est élaboré et défini par l'opérateur, en réponse aux données de cartes et aux données de traces de failles de la figure 28a. En fait, le modèle structurel 66a pour la piste 1 est défini dans une certaine mesure comme si le "modèle géologique de cellules fines" pour la piste 2 (figure 28b) avait été choisi par l'opérateur pour élaborer le modèle structurel 66a. Modélisateur de propriétés 66b Sur les figures 24 et 30, le modélisateur de propriétés 66b de la figure 24 applique certaines "propriétés" à une "couche" de formation terrestre, ces propriétés comportant, par exemple, la porosité, la perméabilité, etc. SR 15546 US/RS Pour la piste 1 de la figure 28a, le modélisateur de propriétés 66b utilise des "cartes de propriétés" pour attribuer certaines "propriétés" comme la porosité et la perméabilité, à chacune des couches de formation terrestre. Sur la figure 30, par exemple, il y a trois cartes (ou horizons), une première "carte 1" 120 représentant une première couche d'horizon, une deuxième "carte 2" 122 représentant une deuxième couche d'horizon et une troisième "carte 3" 124 représentant une troisième couche d'horizon. Une couche de "cartes de propriétés 1" 128 est définie par la "carte 1" 120 et la "carte 2" 122 et une couche de "carte de propriétés 2" 130 est définie par la "carte 2" 122 et la "carte 3" 124. Le modélisateur de propriétés 60b modélise (ou représente) certaines propriétés telles que la porosité et l'imperméabilité, par rapport à la couche de "carte de propriétés 1" 128 et le modélisateur de propriété 66a applique également d'autres propriétés comme la porosité et la perméabilité à la couche de "carte de propriétés 2" 130. Pour la piste 2, toutefois, par rapport au "modèle géologique de cellules fines", des propriétés sont déjà présentes dans chaque petite cellule du "modèle géologique de cellules fines" Parfois, pour la piste 1 ou la piste 2, l'utilisateur peut appliquer plus de propriétés au modèle en effectuant des calculs simples pour déduire une propriété d'une autre. A présent, puisqu'un "modèle géologique unique" 66a a été élaboré par l'ossature structurelle 149 en réponse aux données de la piste 1 ou de la piste 2, à ce stade les deux pistes (piste 1 et piste 2) SR 15546 US/RS fusionnent. Avant ce stade, il y avait deux façons principales de fourles données d'entrée 69 à Flogrid 60a, et ces deux façons principales nécessitaient deux pistes pour élaborer le modèle structurel. Le "modèle géologique simple" est à présent généré par l'ossature structurelle 149 et entré dans le modèle de propriétés 66b. Toutefois, avant de traiter du modèle de propriété 66b, remarquons la fonction du dispositif de quadrillage progressif 84 de la figure 24 qui est traitée ci-dessous. Dispositif de quadrillage progressif 84 de la figure 24 A ce stade, l'utilisateur peut vouloir simplifier le "modèle géologique unique" 66a, entré par l'ossature structurelle 142 dans le modèle de propriétés 66b de la figure 24, en utilisant le dispositif de quadrillage progressif 84 de la figure 24 pour regrouper les couches du "modèle géologique unique" (appelé "modèle géologique"). Par exemple, le "modèle géologique" peut être composé de 600 couches, et l'utilisateur peut utiliser le dispositif de quadrillage progressif 84 pour regrouper les couches du "modèle géologique", soit par la vue (visuellement) soit par un "algorithme". Lorsque les couches du "modèle géologique" sont regroupées par le dispositif de quadrillage progressif 84 en utilisant la méthode "d'algorithme", on peut utiliser deux méthodes: la méthode dite de regroupement de variations" (avec laquelle l'utilisateur peut regrouper des couches semblables) et la méthode dite de "groupement fondé sur le flux" (avec SR 15546 US/RS laquelle l'utilisateur peut pomper de l'eau à travers des couches et limiter la quantité de flux qu'une couche quelconque peut transporter). Considérons la description détaillée suivante du dispositif de quadrillage progressif 84 qui est décrit dans les paragraphes suivants. Le dispositif de quadrillage progressif 84 est utilisé pour regrouper des couches dans un modèle géologique en couches de simulation plus grossières. Il y a deux façons automatiques de le faire, en plus de la façon manuelle décrite dans le guide utilisateur: (1) une méthode fondée sur la variation qui regroupe les couches de sorte que la variation d'une couche est contrôlée, et (2) une méthode fondée sur l'écoulement qui regroupe les couches en contrôlant le flux total à l'intérieur d'une couche. Les interfaces des méthodes de variation et fondées sur l'écoulement ont été conçues pour être semblables, bien que les algorithmes sous-jacents soient différents. Il y a plusieurs options dans cette approche qui présentent une option correspondante dans l'autre approche. Les algorithmes traités ci-dessus sont décrits en référence aux couches de regroupement. Les lignes et les colonnes sont regroupées d'une façon directement analogue mais la méthode fondée sur l'écoulement présente des différences mineures qui seront détaillées ci-après. Dispositif de quadrillage progressif 84 - Quadrillage fondé sur la variation Cette méthode a été mise au point pour améliorer l'efficacité des méthodes de quadrillage fondées sur la SR 15546 US/RS variance décrites précédemment. Ces méthodes règlent les coins des blocs de quadrillage de simulation, de sorte que la variance de la propriété choisie par certains utilisateurs est répartie à distance égale sous la contrainte d'un nombre total spécifié de blocs de quadrillage. La technique peut fonctionner correctement. Toutefois, pour généraliser (lorsque l'on regroupe des cellules géologiques en cellules de simulation), elle est lente - principalement parce qu'un échantillonnage est nécessaire. Dans le cas spécial o les couches de quadrillage de simulation correspondent à un groupe exact de couches géologiques, l'algorithme est beaucoup plus rapide. Toutefois, on a toujours besoin de la solution d'un problème d'optimisation combinatoire compliqué et donc d'une entrée utilisateur obscure. La difficulté surgit de la nature non additive de la variance. Si on ajoute une couche de variance élevée à une couche de variance faible, la variance totale peut diminuer. Dans les paragraphes suivants, la notion de variabilité ou variation (qui est strictement croissante) est introduite. Considérons un modèle géologique ayant Ng couches et définissons la variation d'une propriété, p, entre deux couches en fonction de l'expression: Xo = 0 Xk = nInJ y Pjlk - P.jk+., Eq. 1 X Ng=O O Pi,j,k est la valeur de la propriété dans la cellule géologique (i,j, k), nInJ est le nombre total de cellules dans la couche et la somme se fait sur toutes SR 15546 US/RS les paires de cellules qui présentent un volume non nul. Définissons Vo et VNg pour qu'elles conviennent dans la définition de la variation bilatérale. La somme est effectuée pour toutes les valeurs de k entre 1 et Ng-1. La variation totale de tout le modèle est définie par la somme des variations entre les couches. Définissons la variation bilatérale d'une couche par l'expression: vtk =0,5(vk-l+Vk) Eq. 2 Ordonnons les variations bilatérales des couches de telle sorte que la couche ayant l'indice de couche tn soit la couche ayant la énième plus grande variation bilatérale. Si les couches géologiques sont regroupées, de façon consécutive, en couches de simulation, définissons la variation dans n'importe quelle couche grossière en fonction de l'équation suivante: vs = Z Xk Eq. 3 kas dans laquelle la somme se fait sur toutes les limites de couches géologiques à l'intérieur de la couche de simulation. L'utilisateur fixe ensuite le paramètre suivant: (1) une limite supérieure concernant la variation de la couche de simulation exprimée en tant que pourcentage de la variation totale du modèle, (2) une limite supérieure est fixée concernant le nombre de couches géologiques qui peuvent être groupées dans une couche de simulation et (3) le nombre de couches devant être conservées en tant que couches uniques. C'est-à-dire que des couches géologiques spécifiques doivent rester SR 15546 US/RS des couches de simulation. Si l'utilisateur décide que R couches doivent être conservées sous forme de couches uniques, ces R couches seront celles qui présentent les indices de couches t les plus élevés. Le regroupement commence à la base du modèle de sorte que le dernier groupe de couches, celui situé en haut, peut présenter une variation considérablement plus petite que le maximum spécifié par l'utilisateur. Dispositif de quadrillage progressif 84 - Quadrillage progressif fondé sur l'écoulement Dans cette approche, on résout deux problèmes d'écoulement incompressibles à phase unique sur le réservoir, de gauche à droite et du bas vers le haut. * On fera référence à ces expériences en les appelant expériences I et J. Lors des expériences, il est impératif d'imposer des limites de pression constante sur deux faces verticales opposées du réservoir et une limite d'écoulement nul sur les quatre autres. La pression sur une face est fixée à l'unité et sur la face opposée elle est fixée à zéro. Sans qu'il y ait de perte de généralité, on fixe la viscosité à l'unité. Les calculs d'écoulement sont effectués d'une manière complexe en fixant les transmissibilités perpendiculaires au sens de la chute de pression à zéro. Ceci accélère considérablement les calculs. Une fois que les calculs d'écoulement sont effectués, les flux des couches sont calculés en accumulant le flux volumétrique, indépendamment d'une couche à l'autre, dans chaque face entre toutes les paires de blocs de quadrillage et toutes les faces extérieures de la couche. Il faut remarquer que les SR 15546 US/RS flux de couche comportent les contributions des deux expériences I et J. Une fois que l'on a trouvé ces flux de couche, on ordonne les flux de couche pour obtenir l'indice de flux pour chaque couche. Le flux maximum est défini comme étant le flux de la couche présentant un indice de flux égal à un, c'est-à-dire le flux de la couche qui présente le flux maximum. L'utilisateur doit alors fixer trois paramètres: (1) une limite supérieure du flux total de la couche de simulation, exprimée en tant que pourcentage du flux maximum, (2) une limite supérieure du nombre de couches géologiques qui peuvent être regroupées en une couche de simulation, et (3) le nombre de couches à conserver comme couches uniques: c'est-à-dire que des couches géologiques spécifiques doivent rester des couches de simulation; si l'utilisateur décide que R couches doivent être maintenues sous la forme de couches uniques, ces R couches devront être celles qui présentent les R indices de flux les plus élevés. Une fois encore, le regroupement commence à la base du réservoir. Il faut remarquer que lors du regroupement de lignes ou de colonnes, à la place d'un flux de couche, on dispose d'un flux de ligne ou de colonne. Ces flux utilisent encore les expériences I et J décrites ci- dessus. Le flux de ligne (regroupant les lignes de cellules présentant la constante j) utilise uniquement les flux d'expérience I. Le flux de colonne (regroupant les colonnes de cellules présentant la constante I) utilise uniquement les flux d'expérience 1. SR 15546 US/RS Dispositif de quadrillage structuré 68 de la figure 24 Editeur de limite 68a de la figure 26 - Elaboration d'une limite En référence aux figures 31, 32 et 33, et en référence tout d'abord à la figure 31, supposons que la figure 31 représente une vue de dessus d'un horizon, tel que l'un des horizons 155 de la figure 25a, la vue de dessus de l'horizon 155 comportant des intersections 134 des failles 157 traversant l'horizon 155. Sur la figure 31, supposons qu'une limite 136 est tracée à l'intérieur de l'horizon 155, la limite 136 renfermant les intersections 134 des failles 157. Par conséquent, sur la figure 31, si l'on regarde l'horizon 155 par au-dessus (ce qui représente une "vue de surface"), on peut voir une limite spécifiée par l'utilisateur 136 contenant les lignes d'intersections 134 qui représentent les intersections des failles 157 qui traversent l'horizon 155 (voir figure 25a). Les lignes d'intersection des failles 134 de la figure 31 représentent les "traces de failles" d2 de la figure 28a et l'utilisateur effectue une numérisation dans une "limite" 136 autour des failles de traces 134, 30d2, et l'utilisateur envisage d'effectuer le "quadrillage de surface structuré" de la présente invention. Un quadrillage "structuré" nécessite une limite 136 ayant quatre (4) coins logiques et la forme de la limite 136 commence à ressembler à un "rectangle" ayant subi une distorsion. Par conséquent, un "quadrillage structuré" est une version ayant subi une distorsion d'un quadrillage rectangulaire. SR 15546 US/RS Dispositif de quadrillage de surface structuré 68b de la figure 26 Elaboration d'un "quadrillage de surface structuré" Sur la figure 32, pour élaborer un "quadrillage de surface structuré" sur l'horizon 155 de la figure 31, on commence par dessiner un ensemble complètement uniforme de rectangles à l'intérieur de la limite 136, tous les rectangles présentant la même forme et la même taille, tel qu'illustré sur la figure 32. On considère que les rectangles de la figure 32 sont élaborés dans un matériau présentant des caractéristiques comprises entre les caractéristiques d'une tige d'acier et d'une bande élastique; c'est-à-dire que le matériau des rectangles de la figure 32 est plus résistant qu'une bande élastique mais plus élastique qu'une tige en acier. Ensuite, on plie les rectangles de la figure 32 pour s'adapter à certaines caractéristiques qui nous intéressent. Après le pliage, la limite ayant subi une distorsion "136 ayant subi une distorsion" est illustrée sur la figure 33. N'importe quelle distorsion des rectangles de la figure 32 est valable tant qu'aucune des cellules de la figure 32 n'est "inversée". Toutefois, au cours du pliage des rectangles de la figure 32, un objectif principal est de faire subir aux rectangles de la figure 32 une distorsion afin de capturer le plus de lignes d'intersection de failles 134 possible situées dans l'horizon 155. Sur la figure 34, la limite 136 de la figure 31 est composée de deux parties: les lignes "i" et les lignes "j", les lignes "i" représentant le côté gauche et le côté droit de la limite 136 et les lignes "j" SR 15546 US/RS représentant le haut et le bas de la limite 136. Par exemple, sur la figure 34, la limite 136 présente une ligne "i" 136a qui forme le côté gauche de la limite 136, une autre ligne "i" 136b formant le côté droit de la limite 136, une ligne "j" 136c formant le haut de la limite 136 et une autre ligne "j" 136d formant le bas de la limite 136. Ensuite, sur la figure 34, on divise les lignes d'intersection de failles 134 de la figure 31 en trois classes: des failles "i", des failles "j" et des failles "en zigzag". On effectue un quadrillage sur les failles "i" et les failles "j" mais par sur les failles "de zigzag". Par exemple, sur la figure 34, la ligne d'intersection de faille 140 est une faille "i", la ligne de faille 142 est une faille "j", la ligne d'intersection de faille 143 est une autre faille "i" et la ligne d'intersection de faille 144 est une faille "en zigzag". Sur les figures 32, 34 et 35, en utilisant le quadrillage non déformé de la figure 32, on plie le quadrillage non déformé de la figure 32 pour s'adapter à la limite déformée 136 (composée de 136a, 136b, 136c, 136d) de la figure 34. Dans la figure 34, le terme "lignes de contrôle" est défini comme étant chacune des lignes limite 136 et chacune des "failles étendues". Une "faille étendue", cependant, est composée d'une "ligne d'intersection de faille" plus une ou plusieurs "extensions en pointillés" qui servent d'interconnexion entre la "ligne d'intersection de faille" et la limite 136. Sur la figure 35, les "lignes limite" 136 sont indiquées par les numéros de référence 136a, 136b, 136c et 136d, les "lignes d'intersection de failles" sont indiquées par les numéros de référence 140, 142 et 143 SR 15546 US/RS et les "extensions en pointillés" sont indiquées par les numéros de référence 146a, 146b, 146c, 146d, 146e et 146f. Par conséquent, une "ligne de contrôle" sur la figure 35 est composée de chacune des lignes limite 136a, 136b, 136c et 136b, et une "ligne de contrôle" sur la figure 35 est également composée de chacune des lignes d'intersection de failles "i" et "j" 140, 142 et 143 couplées avec ladite une ou lesdites plusieurs extensions en pointillés 146a à 146f connectées, respectivement, aux lignes d'intersection de failles , 142, 143. Par exemple, une première "ligne de contrôle" pourrait être la ligne d'intersection de faille "i" 140, l'extension en pointillés 146a pouvant connecter la ligne d'intersection de faille 140 à la limite supérieure 136c et l'extension en pointillés 146c connecter la ligne d'intersection de faille 140 à la limite inférieure 136d. Une deuxième "ligne de contrôle" pourrait être la ligne d'intersection de faille "j" 142, l'extension en pointillés 146e pouvant connecter la ligne d'intersection de faille 142 à la limite gauche 136a et l'extension en pointillés 146f connecter la ligne d'intersection de faille 142 à la limite droite 136b. Une troisième "ligne de contrôle" pourrait être la ligne d'intersection de faille "i" 143, l'extension en pointillés 146b pouvant connecter la ligne d'intersection de faille 143 à la limite supérieure 136c et l'extension en pointillés 146d connecter la ligne d'intersection de faille 143 à la limite inférieure 136d. Des quatrième, cinquième, sixième et septième "lignes de contrôle" peuvent être composées des limites 136a, 136b, 136c et 136d sur la figure 35. SR 15546 US/RS Dans la technique antérieure, l'utilisateur numérise et introduit manuellement les "lignes de contrôle". Toutefois, le dispositif de quadrillage de surface structuré 68b de Flogrid 60a de la présente invention élabore avec précaution et automatiquement les "failles étendues" (lignes d'intersection de failles 140, 142, 143 auxquelles sont connectées les extensions en pointillés 146a à 146f), de sorte que ces failles étendues apparaissent comme étant uniformes. Le dispositif de quadrillage de surface structuré 68b élabore alors un quadrillage à l'intérieur de la limite 136 de l'horizon 155 de la figure 35 en ajoutant des "lignes supplémentaires" 150 à l'intérieur de la limite 136, l'emplacement de chacune des "lignes supplémentaires" 150 à l'intérieur de la limite 136 étant déterminé en effectuant une interpolation entre les "lignes de contrôle" adjacentes à l'intérieur de la limite 136. Lorsqu'il aura élaboré automatiquement les "lignes de contrôle" sur la figure 35, l'utilisateur indique alors combien de "lignes supplémentaires" 150 il souhaite placer de façon adjacente aux ou entre les "lignes de contrôle". Par exemple, sur la figure 35, l'utilisateur indique que des "lignes supplémentaires" 150 seront placées de façon adjacent aux entre les "lignes de contrôle" de la figure 35. Les "lignes de contrôle" de la figure 35 sont placées dans une liste de la façon suivante: une première ligne de contrôle (140, 146a, 146c), une deuxième ligne de contrôle (143, 146b, 146d), une troisième ligne de contrôle (142, 146e, 146f), une quatrième ligne de contrôle 136a, une cinquième ligne de contrôle 136b, une sixième ligne de contrôle 136c et une septième ligne de contrôle 136d. SR 15546 US/RS Dispositif de quadrillage de bloc 68c de la figure 26 - Lignes de coordonnées réduites Sur les figure 35 et 35a, une fois que les lignes supplémentaires 150 sont ajoutées au quadrillage de la figure 35, on place chaque noeud (c'est-à-dire les sommets au point d'intersection) du quadrillage de la figure 35, on connecte une "ligne de coordonnées" à chaque noeud du quadrillage de la figure 35 et on réduit chaque "ligne de coordonnées", dans la formation terrestre (c'est-à-dire dans la page de la figure 35), en permettant aux "lignes de coordonnées" de couper d'autres cartes ou horizons de la formation terrestre et on détermine ensuite l'endroit o la "ligne de coordonnée" coupe les "cartes" ou horizons de la formation terrestre (rappelons que le terme "carte" fait également référence à un "horizon"). Un "quadrillage de simulation structuré" composé d'une pluralité de "blocs de quadrillage" est obtenu en réduisant les "lignes de coordonnées" à partir des noeuds ou des sommets du quadrillage de surface de la figure 35 et en permettant à ces lignes de coordonnées de couper les horizons ou cartes situés le plus en bas dans la formation terrestre, ce qui permet de former une pluralité de "blocs de quadrillage". Un "bloc de quadrillage" est défini en référence à la figure 35a. Sur les figures 35a et 36, par exemple, une coupe d'une formation terrestre, semblable à la coupe de la formation terrestre 153 de la figure 24a est illustrée. Sur les figures 35a et 36, en référence tout d'abord à la figure 36, une coupe d'une formation terrestre est représentée sur la figure 36 et une faille réelle 152 est représentée et coupe une SR 15546 US/RS 113 - pluralité d'horizons 154. Une pluralité de lignes de coordonnées 156 est réduite à partir des noeuds ou sommets du quadrillage de la figure 35, par le biais des horizons 154 et les intersections entre les lignes de coordonnées 156 et les horizons ou cartes 154 de la figure 36 sont indiqués. Le terme "sommets de quadrillage de la figure 35" se rapporte aux intersections, sur la figure 35, entre les "lignes de contrôle" et les "lignes supplémentaires" 150 ou bien des intersections entre les "lignes de contrôle" ou les intersections entre les "lignes supplémentaires" 150. Sur la figure 36, les intersections 158 sont déterminées comme existant entre chacune des lignes de coordonnées "réduites" 156 et chacun des horizons 154. Les intersections 158 de la figure 36 définissent quatre des huit coins et "blocs de quadrillage" qui vont être utilisés dans le simulateur 40c2A de la figure 16. Sur la figure 35a, une vue tridimensionnelle de la coupe de la formation terrestre de la figure 36 est illustrée. Sur la figure 35a, quatre lignes de coordonnées 156 sont connectées à quatre noeuds ou sommets 159 du quadrillage des figures 35 et 35a et les quatre lignes de coordonnées 156 sont réduites dans la formation terrestre, jusqu'à ce que les lignes de coordonnées 156 coupent deux autres horizons situés le plus bas 161 et 163 dans la formation terrestre. A ce stade, deux "blocs de quadrillage" 165 et 167 sont formés. Sur la figure 35a, le premier "bloc de quadrillage" 165 est limité par les lignes de coordonnées 156 et les noeuds 159 sur le "quadrillage de surface structuré" et les points d'intersection 156a, 156b, 156c et 156d de l'horizon 161; le deuxième SR 15546 US/RS "bloc de quadrillage" 167 est limité par les lignes de coordonnées 156 et les intersections 156a à 156d sur le premier horizon 161 et les intersections 163a à 163d sur le deuxième horizon 163. En utilisant la technique mentionnée ci-dessus en référence à la figure 35a, une multitude de ces "blocs de quadrillage" peut être mise au point en réduisant les lignes de coordonnées à partir de chacun des noeuds ou sommets du quadrillage de la figure 35 et en faisant se couper les lignes de coordonnées et les horizons situés le plus bas dans la formation terrestre. Sur la figure 37, en ce qui concerne les lignes d'intersection de failles en "zigzag" 144 de la figure , afin d'effectuer un quadrillage sur ces failles en "zigzag" 144, on les fait "s'engager" avec le noeud de quadrillage le plus proche. Sur la figure 37, la faille en "zigzag" 144 "s'engage avec le noeud de quadrillage le plus proche", ce qui permet de réduire les failles de la figure 37. En conséquence, on a modifié la forme de la faille en zigzag 144 de la figure 37 pour qu'elle corresponde aux sommets à 160a à 160d du quadrillage de la figure 37. Sur la figure 38, les lignes de coordonnées 162 sont représentées comme étant adjacentes à une faille en zigzag 164 qui coupe les horizons 166. Le dispositif de quadrillage de surface structuré 68b de la figure 26 s'étend hors des horizons 166 pour couper la ligne de coordonnées 162 afin de produire les intersections 168. Le dispositif de quadrillage de surface 68b utilise alors les intersections 168 entre l'horizon étendu 166 et les lignes de coordonnées 162 afin de former les coins des "blocs de quadrillage" du quadrillage structuré. SR 15546 US/RS Le dispositif de mise à une plus grande échelle Le quadrillage a été accompli et à présent il faut effectuer la mise à une plus grande échelle. En utilisant le dispositif de mise à une plus grande échelle 70 de la figure 24 pour effectuer la fonction des mises à une plus grande échelle, Flogrid 60a attribue certaines propriétés à chacun des "blocs de quadrillage" du modèle de simulation. Le dispositif de mise à une plus grande échelle 70 fonctionne uniquement avec le dispositif de quadrillage structuré 68 et non pas avec le dispositif de quadrillage non structuré cl de la figure 16, puisque le dispositif de quadrillage non structuré 40cl est capable d'attribuer seul les propriétés nécessaires au modèle de simulation. Sur la figure 39, traditionnellement, le dispositif de mise à une plus grande échelle 70 attribue des propriétés à chacun des "blocs de quadrillage" du modèle de simulation en échantillonnant la propriété au centre de chacun des "blocs de quadrillage". Par exemple, le dispositif de mise à une plus grande échelle 70 peut échantillonner la "porosité" ou la "perméabilité" au centre du "bloc de quadrillage". Toutefois, en pratique, les cartes (horizons) varient au centre du "bloc de quadrillage". Par conséquent, l'approche traditionnelle peut ne pas fonctionner correctement. Sur la figure 40, toutefois, dans la pratique réelle, le dispositif de mise à une plus grande échelle "quadrille avec précision" le "bloc de quadrillage". SR 15546 US/RS Si, par exemple, il y a une chute de pression dans le "bloc de quadrillage", le dispositif de mise à une plus grande échelle 70 estime l'écoulement réel dans le bloc de quadrillage puis, en tenant compte de toute la variabilité, le dispositif de mise à une plus grande échelle 70 définit la moyenne ou la valeur effective du débit d'écoulement dans le "bloc de quadrillage" comme étant égale à la valeur de débit d'écoulement, dans un bloc de quadrillage à pression uniforme, qui donnerait le même débit d'écoulement pour cette chute de pression. Considérons la description détaillée suivante du dispositif de mise à une plus grande échelle 70 qui est mentionnée dans les paragraphes suivants. Le dispositif de mise à une plus grande échelle exécute des procédés de mise à une plus grande échelle dans Flogrid 60a en fonction du concept selon lequel les propriétés mises à une plus grande échelle sont des valeurs semblables à celles mesurées en laboratoire. Ainsi, les "mesures" sont conçues pour calculer certaines propriétés mises à une plus grande échelle indépendamment des autres. On suppose que la porosité sera mise à une plus grande échelle indépendamment de toute autre propriété. De plus, les propriétés d'un "bloc de quadrillage" sont calculées au même moment, ce qui revient à découper un "bloc de quadrillage" provenant d'un réservoir de formation terrestre et à effectuer des expériences sur ce "bloc de quadrillage" dans un laboratoire. Dispositif de mise à une plus grande échelle 70 - Construction du quadrillage précis SR 15546 US/RS 117 - Chaque bloc de quadrillage 165, 167 est supposé être défini par ses huit points de coin. Sur la figure , afin d'effectuer un quadrillage précis du "bloc de quadrillage", un quadrillage "Nx" par "Ny" par "Nz" est construit à l'intérieur du "bloc de quadrillage" en utilisant une formule d'interpolation simple (se référer aux paragraphes suivants pour une description de cette formule d'interpolation). Dispositif de mise à une plus grande échelle 70 - Echantillonnage de propriétés Une fois que l'on a construit un "quadrillage précis" à l'intérieur du "bloc de quadrillage" (le quadrillage précis contenu à l'intérieur du bloc de quadrillage étant composé d'une pluralité de cellules), on peut échantillonner les différentes propriétés qu'il faut mettre une plus grande échelle au centre de chaque cellule de quadrillage précis. On suppose que les propriétés sont constantes pièce par pièce dans le quadrillage précis. Flogrid 60a choisit le nombre de cellules à l'intérieur de chaque "bloc de quadrillage quadrillé avec précision", de sorte que la taille de chaque "cellule précise" est environ celle de la cellule de modèle géologique la plus petite contenue dans le bloc du quadrillage grossier. On indiquera "Vb" du volume de bième quadrillage d'une "cellule de quadrillage fin". Pour des raisons d'efficacité, Flogrid 60a utilise une approximation de ces volumes de bloc. L'estimation du volume total de la cellule grossière obtenu en additionnant ces valeurs de SR 15546 US/RS quadrillage précis ont tendance à être la valeur exacte puisque le quadrillage précis est affiné. Dispositif de mise à une plus grande échelle 70 - Mise à une plus grande échelle des propriétés réelles non dimensionnelles Les propriétés inconnues, dont les dimensions (au sens de Darcy, mètres, etc.) ne sont pas celles de la perméabilité, sont supposées être desensembles de nombres réels sans dimension. La moyenne d'échelle grossière de chaque bloc de quadrillage grossier dispositif de mise à une plus grande échelle 70 est définie comme étant la moyenne pondérée du volume sur les cellules précises: Y VbPb P = b Eq. 1 b Pb étant la propriété dans le bième bloc précis et PB la valeur mise à une plus grande échelle dans le bloc grossier. Dispositif de mise à une plus grande échelle 70 - Mise à une plus grande échelle de la porosité La porosité est moyennée en utilisé le volume de cellule précis comme pondération. Ainsi: VbPb P b VEq. 2 Z Vb b Dispositif de mise à une plus grande échelle 70 - Mise à une plus grande échelle de l'imperméabilité absolue SR 15546 US/RS Dans les méthodes de "moyennage simple", on ne justifie pas la formule de moyennage en référence à une solution d'écoulement approximative. Dans le cas d'un moyennage arithmétique et harmonique, on pondère l'imperméabilité à l'aide de facteurs géométriques, de sorte que, dans les cas unidimensionnel dans lesquels le procédé est semblable à un procédé fondé sur l'écoulement, on arrive aux mêmes solutions. Ceci revient à prendre en considération la zone de coupe transversale et de longueur de chaque bloc de quadrillage à échelle précise. Ceci équivaut à utiliser le volume et la longueur au carré. Dans le cas d'un moyennage de puissance et géométrique, on n'inclut pas les facteurs géométriques car ils s'annulent. On présente uniquement la formule pour effectuer une mise à une plus grande échelle de la perméabilité dans le sens I. Les autres directions présentent la même forme. Dans les paragraphes suivants, les équations d'écoulement à phase unique dans un système de coordonnées trilinéaires général sont présentées. L'étude de ces équations peut aider les lecteurs qui souhaitent déduire les méthodes harmoniques- arithmétiques et arithmétiques-harmoniques. Les formules suivantes sont disponibles, dans lesquelles KI et la perméabilité absolue ayant subi une mise à une plus grande échelle dans le sens I, Vijk est une estimation du volume de la cellule, Diijk est une estimation de la longueur de bloc de quadrillage dans le sens I et Kiijk est la perméabilité absolue dans le ième bloc de quadrillage dans le sens I, pour la cellule (i,j,k): Moyennage arithmétique SR 15546 US/RS KIf VKk = E Vkl-k- Eq. 3 dIjk dIvk l'addition se faisant sur tous les blocs de quadrillage précis. Moyennage harmonique i di2 dl2 s EIE q. 4 KI,. V.jk Vjk kijk l'addition se faisant sur tous les blocs de quadrillage précis Moyennage géométrique N In K/= Ink kI Eq. 5 l'addition se faisant sur tous les blocs de quadrillage précis. Moyennage de puissance KPw =JN EkIk Eq. 6 l'addition se faisant sur tous les blocs de quadrillage précis. Dispositif de mise à une plus grande échelle 70 - Solutions unidimensionnelles composites Dans ces méthodes, on taille le quadrillage progressif dans un ensemble de tubes et on ignore l'écoulement transversal entre les tubes, ou bien on taille le quadrillage précis dans un ensemble linéaire de dalles en ignorant le flux transversal à l'intérieur de chaque dalle. Dans ce qui suit, on suppose que la viscosité est égale à l'unité; ceci simplifie les équations sans perte de généralité et sans affecter les valeurs de perméabilité ayant, été mise à une plus grande échelle. SR 15546 US/RS 121 - Dispositif de mise à une plus grande échelle 70 - Moyennage harmonique-arithmétique Dans le premier cas, on calcule le flux, uIjk, à travers chaque tube (j,k) séparément. Cela donne: UIjk Yev Uk = AP Eq. 7 ijk ijk Le flux total à travers le bloc de quadrillage est ensuite obtenu en additionnant uIjk, le flux à travers le tube fin (j,k) sur j et k. On calcule ensuite la perméabilité effective ayant été mise à une plus grande échelle, à partir de: 1> = dl,* E1 l2 Eq. 8 Z _2 2 Eq. 8 j dI k dlk ik i Vjk Hyk j Vijk Dispositif de mise à une plus grande échelle 70 - Moyennage arithmétique-harmonique Dans ce cas, on suppose une chute de pression dpi dans la ième dalle. Remarquons que, du fait de l'incompressibilité, le flux qui traverse chaque dalle peut être le même, le flux étant donné par la formule: u = pL-d2 kI,,k Eq. 9 jk ijk A présent, ajoutons les chutes de pression sur les dalles, et on obtient: PI =Po =-PL V=k Eq. 10 Z 2 j,k dl]k o po est la pression d'entrée et PL la pression de sortie. Mettons en équation la chute de pression en SR 15546 US/RS se basant sur l'expression similaire présentant la constante KI, la perméabilité ayant été mise à une plus grande échelle, et on obtient l'expression finale suivante: ! i\T ijk K= dIk Eq. 11 E _i; \';k Hijky dI2 j,k qk Dispositif de mise à une plus grande échelle 70 - Solutions tridimensionnelles complètes La solution numérique tridimensionnelle complète fondée sur l'écoulement est l'option de mise à une plus grande échelle à phase unique qui demande le plus de calculs. Elle est mise en oeuvre est utilisant aucune limite sans écoulement. Dans ce cas, on calcule le flux total qui traverse le système dans le sens de la chute de pression. La perméabilité effective est ensuite estimée en résolvant le même problème avec des perméabilités constantes choisies pour donner le même flux. Ceci est effectué en utilisant le procédé de moyennage harmonique-arithmétique. Une erreur est introduite par cela, à moins que les blocs de quadrillage soient rectangulaires. L'utilisation de la limite linéaire élimine le besoin de cette approximation. Dispositif de mise à une plus grande échelle 70 - Elaboration du quadrillage précis dans un bloc de quadrillage SR 15546 US/RS 123 - Sur les figures 41 et 42, après avoir choisi notre bloc de quadrillage, on divise le bloc en un ensemble de blocs précis Nx par NY par Nz. Ceci est effectué en utilisant une interpolation trilinéaire à partir de valeurs de coin. Une façon de considérer cela est d'imaginer que le bloc de quadrillage est une version déformée d'un cube unitaire. Sur un cube unitaire, on peut introduire un système de coordonnées (,n,) dans lequel chacune des variables se trouve entre 0 et 1. Ensuite on définit la déformation du cube unitaire jusqu'au bloc de quadrillage réel par la formule d'interpolation trilinéaire (voir figure 42): x(t,,4) = {[Xli - 4) + x24](1 - 'a) + [X3(1 - 4) + X44]rI}(1 -) + {[X5 (1 - + x64Xv1 - Ti) + [xX( + x Sur la figure 41, des étiquettes de coin sur le cube de référence et le bloc de quadrillage sont illustrées. Sur la figure 41, les points 1 à 4 se trouvent sur le dessus et les points 4 à 8 sur la base, tel que représenté sur la figure 41. On définit un quadrillage précis de blocs de quadrillage rectangulaires présentant chacun la même forme et la même taille, à l'intérieur du bloc de référence. La transformation déforme ensuite toutes ces cellules de référence précises en des blocs de quadrillage précis déformés. Sur la figure 42, l'étape d'interpolation d'un quadrillage précis 2xlx3 dans un bloc de quadrillage est illustrée. Dispositif de mise à une plus grande échelle 70 - Equation d'écoulement de phase unique avec coordonnées curvilignes SR 15546 US/RS La dérivation des équations d'écoulement en coordonnées générales est effectuée plus facilement en utilisant du matériel d'analyse de tenseur. On suppose, dans le cas normal, que le système de coordonnées trilinéaires est orthogonal, même lorsque cela n'est manifestement pas le cas. Dans ce qui suit, la convention d'addition concernant des indices répétés est appliquée. Pour un écoulement incompressible, la loi de l'équilibre des masses est donnée par: aui= 0 gu a= O o,' est la ième coordonnée curviligne (c,n, ), i = 1,2,3 respectivement. W est la racine carrée du déterminant du tenseur métrique et est le volume de l'espace physique par l'unité de volume d'espace de calcul en ( ,n, ). La loi de Darcy donne: ki ap u = 9- pgp o k' est la composante physique du tenseur de perméabilité dans le ième sens. L'approximation d'orthogonalité se produit lorsqu'on ignore les termes situés hors de la diagonale du tenseur de perméabilité transformé dans les coordonnées trilinéaires. gii est le ième composant du tenseur métrique et correspond à la longueur au carré de l'espace physique par unité de longueur au carré de l'espace de calcul en ( SR 15546 US/RS - niveau du centre du bloc donnée par (,n,) = (0,5, 0,5, 0,5). De cette manière, les équations dans le système curviligne présentent la même forme qu'un système équivalent en coordonnées rectangulaires mais présente des perméabilités variant différemment. Description détaillée du dispositif de quadrillage de surface structuré 68b de la figure 26 Le traitement ci-dessus représente une description détaillée du programme de quadrillage de simulation "Flogrid" 60a sur la figure 24, comportant une description détaillée du dispositif de quadrillage structuré 68 de la présente invention qui comporte le "dispositif de quadrillage de surface structuré" 68b de la présente invention. Les paragraphes suivants traitent, plus en détail, de la structure et de l'opération fonctionnelle du "dispositif de quadrillage de surface structuré" 68b de la figure 26 de la présente invention qui se trouve à l'intérieur de et fait partie du dispositif de quadrillage structuré 68 du programme de quadrillage de simulation "Flogrid" 60a de la figure 24. Sur la figure 43, en utilisant l'éditeur de limite 68a de la figure 26, l'utilisateur a déjà numérisé ou spécifié une limite 176 sur un horizon 155, ladite limite 176 entourant toutes les lignes d'intersection de failles de cet horizon. Une "limite" 176 est, en fait, l'enveloppe d'un horizon d'une formation terrestre à travers laquelle une pluralité de failles sont passées. Par exemple, sur la figure 43, une limite 76 renferme une pluralité de lignes d'intersection de failles 170, 172, 174, les SR 15546 US/RS lignes d'intersection de failles 170, 172, 174 existant éventuellement l'un ou plusieurs des horizons 155 de la formation terrestre représentés sur la figure 25a. De plus, l'utilisateur a déjà fait un choix des failles "i" et "j" (lignes d'intersection de failles) qui sont disposées à l'intérieur de cette limite 176. Tel qu'indiqué sur la figure 43, après avoir spécifié la limite 176 sur l'horizon 155, l'utilisateur choisit une paire de failles "i" 170 et 172 et l'utilisateur choisit une faille "j" 174. Dispositif de quadrillage de surface structuré 68b - Elaboration d'un quadrillage de surface structuré Le dispositif de quadrillage de surface structuré 68b du dispositif de quadrillage structuré 68 de la figure 26 élabore un quadrillage de surface structuré en exécutant les étapes fonctionnelles suivantes: Etape 1: construction d'un "quadrillage triangulaire régulier" dont la limite se trouve à l'extérieur de la limite de l'utilisateur 176 de la figure 43 Sur la figure 44, par exemple, un quadrillage de rectangles 178 est tracé sur la limite de l'utilisateur 176 sans aucune référence à la limite 176 autre que la taille générale de la limite 176. Ensuite, on divise chacun des rectangles 178 en "triangles" 178a, 178b, chaque rectangle 178 étant divisé en deux triangles, un premier triangle 178a et un deuxième triangle 178b. Le résultat de cette étape est qu'on a produit un "quadrillage triangulaire régulier", puisque chaque triangle 178a, 178b contenu à l'intérieur présente trois voisins. Sur la figure 44, ce "quadrillage SR 15546 US/RS triangulaire régulier" composé d'une pluralité de triangles 178a, 178b recouvrant la limite de l'utilisateur 176, sera appelé ci-après "triangulation". Etape 2: Elaboration d'une limite et d'une faille respectant un quadrillage non structuré (également appelée "technique de coupe d'isosurface") Sur la figure 44, considérons à leur tour chacune des lignes droites, ou segments de ligne droite de la figure 44 (c'est-à-dire considérons chacune des failles "i" et "j" 170, 172, 174 sur la figure 44 et la limite 176 sur la figure 44) qui se trouvent à l'intérieur du "quadrillage triangulaire régulier" de la figure 44. Sur la figure 44, à titre d'exemple, placer le coin 181 de la limite 176. Le traitement suivant portera sur la manière exacte dont le dispositif de quadrillage de surface structuré 68b effectuera un quadrillage sur ce coin 181. Afin d'effectuer un quadrillage sur le coin 181 de la figure 44, le dispositif de quadrillage de surface structuré 68b de la figure 26 va exécuter les étapes suivantes (1) à (4): (1) Déterminer les "triangles affectés" [c'est-à- dire déterminer ou identifier les triangles spécifiques du "quadrillage triangulaire régulier" de la figure 44 à l'intérieur duquel le segment de ligne (c'est-à-dire le coin 181 sur la figure 44) va être placé]; (2) marquer chaque noeud sur la triangulation existante avec une distance verticale entre le noeud et le segment de ligne, "+" sur un côté du segment de ligne et "-" sur l'autre côté du segment de ligne; en conséquence, le segment de ligne va comporter un SR 15546 US/RS "contour de valeur zéro", qui est défini comme étant la partie du segment de ligne o la distance entre le noeud de triangulation et de segment de ligne est "nulle", (3) si le noeud se trouve à environ une distance nulle du segment de ligne, on l'engage avec "zéro" (distance nulle entre le noeud et le segment de ligne), puis (4) on découpe les triangles le long du "contour de valeur nulle". L'extrémité des failles se trouve dans la triangulation et les points du quadrillage se trouve dans la triangulation. L'objectif majeur des étapes fonctionnelles citées ci-dessous (1) à (4) est le suivant: choisir la triangulation originale qui soit suffisamment précise pour que les triangles ne soient pas découpés plus de deux ou trois fois. Etape 3 - Lorsqu'une des lignes limite ou failles présentant un sens traverse un noeud d'un triangle d'une triangulation de la figure 4, attribuer un "vecteur directionnel" de longueur unitaire (également appelé vecteur directionnel "fixe") dont le sens se dirige vers le noeud, et relier ce vecteur directionnel à ce noeud du triangle de la triangulation Rappelons d'après la figure 43 que les segments de ligne comportent une ligne limite 176, une ligne de failles "i" 170, 172 et une ligne de faille "j" 174, les "vecteurs directionnels" comportant: un vecteur directionnel limite, un vecteur directionnel "i" et un vecteur directionnel "j". Les vecteurs directionnels "i" sont établis de façon à être alignés avec les lignes d'intersection de failles "i", les vecteurs directionnels "j" sont établis de façon à être alignés avec les lignes d'intersection de failles "j" et les SR 15546 US/RS vecteurs directionnels limites sont établis de façon à être alignés avec les lignes limites. Par exemple, sur la figure 45a, les vecteurs directionnels "i" de la figure 45a se dirigent vers le haut et les vecteurs directionnels "j" se dirigent vers la droite. Sur la figure 45bl, rappelons la règle générale concernant les vecteurs directionnels: un "vecteur directionnel" est connecté à la partie de chacun des segments de ligne (170, 172, 174, 176) de la figure bl o la partie du segment de ligne traverse un noeud de l'un des triangles 178a, 178b de la figure 45bl. Le sens de certains vecteurs directionnels est connu; et on appelle ces vecteurs des vecteurs "fixes". Toutefois, le sens de certains autres vecteurs directionnels n'est pas connu; et on peut appeler ces vecteurs des vecteurs "non fixes". Le vecteur directionnel qui est connecté au noeud du triangle de la triangulation à travers lequel une faille ou une limite passe est appelé un vecteur "fixe" car le sens de ce vecteur directionnel est connu et est le même que le sens de la faille ou de la limite d'intersection. Toutefois, le vecteur directionnel qui est connecté à un autre deuxième noeud d'un triangle à travers lequel une faille ou une limite ne passe pas est un vecteur directionnel "non fixe" car le sens de ce vecteur n'est pas connu puisqu'aucune faille ni aucune limite ne traverse ce second noeud. Sur la figure 45b2, un exemple est illustré et celui-ci démontre les définitions des vecteurs directionnels "fixes" et les vecteurs directionnels "non fixes". Une triangulation 183, sur la figure 45b2, comporte une pluralité de "noeuds" 185, un noeud ou SR 15546 US/RS sommet étant défini comme un point o des limites de la triangulation se coupent et forment une intersection l'une avec l'autre. Par exemple, sur la figure 45b2, une première limite de la triangulation 187 coupe et forme une intersection avec une deuxième limite de la triangulation 189 et un "noeud" ou "sommet" 191 est défini par l'intersection entre la première limite de la triangulation 187 et la deuxième limite de la triangulation 189. Un vecteur directionnel "fixe" apparaît au niveau d'un premier noeud particulier de la triangulation lorsqu'une ligne de faille ou une ligne limite forme une intersection avec le premier noeud particulier. Toutefois, un vecteur directionnel "non fixe" apparaît au niveau d'un deuxième noeud particulier de triangulation lorsqu'une ligne de faille ou une ligne limite ne forme pas d'intersection avec ce deuxième noeud particulier. Par exemple, sur la figure 45b2, les vecteurs directionnels 193, 195, 197 et 199 sont des vecteurs directionnels "fixes" puisque ces vecteurs directionnels sont connectés à et ont pour origine "certains noeuds particuliers" 201, 203, 205 et 207 sur la figure 45b2, lesdits "certains noeuds particuliers" formant une intersection avec les lignes d'intersection de failles 170, 172 (extrémité supérieure), 174 et 172 (extrémité inférieure), respectivement. Toutefois, les vecteurs directionnels 209 et 211 sont des vecteurs directionnels "non fixes" puisque ces vecteurs 209, 211 sont connectés à deux (2) noeuds de la triangulation 183 qui ne forment pas d'intersection avec une ligne d'intersection de faille quelconque 170, 172, 174 ou une ligne limite quelconque 176. SR 15546 US/RS Interpolation isotrope et anisotrope La figure 45b3 illustre la façon de déterminer les vecteurs directionnels "fixes" et la façon de déterminer les vecteurs directionnels "non fixes". Pour déterminer les vecteurs directionnels "fixes", lorsqu'un ligne limite ou faille traverse un noeud particulier dans la triangulation, on attribue un vecteur de longueur unitaire au noeud particulier dont le sens est le même que le sens de la ligne limite ou faille au niveau de ce noeud particulier. Pour déterminer les vecteurs directionnels "non fixes", on interpole les champs des vecteurs fixes (ayant un sens connu) sur les noeuds d'autres triangles pour déterminer le sens des vecteurs "non fixes". Puisqu'on connaît la direction des vecteurs directionnels "fixes", on peut utiliser une technique d'interpolation "isotrope" ou "anisotrope" (traitée ci- dessous) pour déterminer le sens des autres vecteurs directionnels "non fixes". Dans les paragraphes suivants, les deux méthodes permettant d'effectuer l'interpolation mentionnée ci- dessus, à savoir isotrope et anisotrope, seront traitées. Sur la figure 46, un ensemble de triangles est illustré et ceux-ci sont tracés autour d'un noeud afin d'illustrer "l'interpolation isotrope". Au cours de "l'interpolation isotrope" (absence d'orthogonalité), on inspecte chaque noeud du triangle, tour à tour. Les vecteurs directionnels de la figure 46 peuvent être "fixes" ou "non fixes". Le vecteur directionnel est "fixe" lorsque son sens est établi et fixé. Toutefois, le vecteur directionnel est "non fixe" SR 15546 US/RS (c'est-à-dire variable) lorsque son sens n'est pas connu ou pas établi ou calé. Si le vecteur est "non fixe", alors on utilise la méthode d'interpolation isotrope, on cale le vecteur "non fixe" égal à la valeur moyenne arithmétique de son voisin. C'est-à-dire qu'on calcule la moyenne arithmétique des six vecteurs voisins et on établit cette moyenne arithmétique égale à la valeur du vecteur non fixe. La méthode d'interpolation "isotrope" n'est pas une méthode d'interpolation moyenne pondérée de distance (comme la méthode d'interpolation anisotrope traitée ci-dessous). Lorsqu'on effectue une interpolation isotrope, on lisse les bosses sur une ou plusieurs failles. Sur la figure 47, en ce qui concerne l'interpolation "anisotrope", cette méthode d'interpolation utilise une technique "de moyenne pondérée de distance". En utilisant la méthode "anisotrope", en regardant les vecteurs "j", on détermine l'endroit o les vecteurs "j" forment une intersection avec les polygones de limite présents sur les triangles, on détermine les vecteurs "i" aux endroits o les vecteurs "j" forment une intersection avec les polygones limitant existant sur les triangles et on place une moyenne pondérée de distance au niveau du noeud au centre. L'effet de cette technique est de choisir les "bosses" qui sont lissées lorsqu'on utilise l'approche de quadrillage isotrope. Plus particulièrement, sur la figure 47, un ensemble de triangles entoure un point ou noeud d'une triangulation et un champ de vecteurs directionnels est connecté aux noeuds de ces triangles. D'abord, on regarde le vecteur "j" 192 de la figure 47 pour déterminer o il coupe ou SR 15546 US/RS forme une intersection avec les triangles voisins. Sur la figure 47, on remarque que le vecteur "j" 197 forme une intersection avec les triangles voisins au niveau de certains "points d'intersection". Au niveau de ces points d'intersection, on détermine le vecteur "i" 190 qui existe au niveau des points d'intersection mentionnés ci- dessus. Prenons une moyenne pondérée existante des deux vecteurs "i" 190, sur la figure 47. Etant donné que, sur la figure 47, le vecteur "i" 190a est "I-", le vecteur "i" 190b est "I" et le vecteur "i" c est "I+", alors une "moyenne pondérée de distance" est définie de la façon suivante: I = [(W+)(I-) + (W-)(I+)]/[W+) = (W-)], o "W+" est la distance existant entre le centre (vecteur "i" 190b) de la figure 47 et le vecteur "I+", et "W-" est la distance entre le centre (vecteur "i" b) de la figure 47 et le vecteur "I-". Le traitement ci-dessus décrit la façon d'interpoler les champs des vecteurs, en utilisant l'interpolation isotrope ou anisotrope, afin de déterminer le sens des vecteurs directionnels "non fixes" en fonction du sens connu des vecteurs directionnels "fixes". Le traitement mentionné ci-dessus décrit la façon d'élaborer le champ de vecteur. Le traitement suivant décrit la façon d'élaborer les "lignes de contrôle". Etape 4: Elaboration des lignes de contrôle Sur les figures 48 et 49, la limite 176 et les lignes d'intersection de failles "i" et "j" 170, 172, 174 ont déjà été établies d'après le traitement ci- dessus en référence à la figure 43. A présent, supposons de la même manière qu'un flux de fluide SR 15546 US/RS circule à une certaine "vitesse" et selon un "sens" le long des extrémités des failles 170, 172, 174. Un "vecteur directionnel" est connecté à l'extrémité des lignes d'intersection de failles 170, 172, 174, ce "vecteur directionnel" ayant un "sens" qui est de la même manière égal au "sens" mentionné précédemment du flux de fluide circulant au niveau de l'extrémité de la faille (le "vecteur directionnel" ayant une unité de longueur). Pour élaborer les "lignes de contrôle", en utilisant notre analogie, on se dirige vers l'extrémité des failles 170, 172, 174 de la figure 48, et on "saute par dessus" lesdites extrémités, dans le flux de fluide qui circule de la même manière avec une vitesse dans un sens le long des failles, et on dérive vers la limite 176. C'est-à-dire qu'on passe de l'extrémité de chacune des failles 170, 172, 174 dans le flux de fluide et on se déplace avec ce fluide pendant un certain moment en se dirigeant vers la limite 176. En utilisant la présente analogie, les "lignes de contrôle" présentent un pas particulier qui est traversé une fois que l'on est passé par dessus, les extrémités des failles 170, 172, 174 et qu'on est déplacé avec le flux de fluide vers la limite 176. En réalité, les "lignes de contrôle" présentent un pas particulier ou un sens particulier qui résulte directement des sens des "vecteurs directionnels" du champ de vecteur qui est connecté au noeud de la triangulation. Sur la figure 49, par exemple, afin d'élaborer une "ligne de contrôle" identifiée par le numéro d'élément "200-1", supposons qu'une faille 196 présente une extrémité 194 et que la faille est située de façon adjacente à une série de triangles 198 de la SR 15546 US/RS triangulation. Au niveau de l'extrémité de la faille 194, un premier vecteur directionnel 194a est connecté à un premier noeud de la triangulation, le premier vecteur directionnel 194a se dirigeant dans le sens de la faille, la longueur du premier vecteur directionnel 194 étant connotée à une "vitesse de fluide". Un deuxième vecteur directionnel 194b est connecté à un deuxième noeud de la triangulation de la figure 49. On passe par dessus l'extrémité 194 de la faille 196, dans le sens du premier vecteur directionnel 194a du champ de vecteur de la figure 49 qui est connecté au premier noeud de la triangulation puis on s'arrête, on détermine la vitesse du fluide en ce point (en utilisant la présente analogie). On se déplace dans le sens du premier vecteur directionnel 194a à cette vitesse pendant un instant. Ensuite, on détermine le sens du deuxième vecteur directionnel 194b du champ de vecteur de la figure 49 qui est connecté au deuxième noeud de la triangulation et on détermine la vitesse du fluide suivant du deuxième vecteur directionnel 194b en ce point. On se déplace dans le sens du deuxième vecteur directionnel 194b à la vitesse de fluide suivante pendant un instant, etc. On continue à dériver vers le haut jusqu'à atteindre la limite. En conséquence, la ligne de commande "200-1" de la figure 49 a été élaborée. * Dans la figure 50, les "lignes de contrôle" 200 sont élaborées en utilisant cette technique. En fait, cette technique est utilisée conjointement avec les lignes d'intersection de failles "i" 170, 172 et les lignes d'intersection de failles "j" 174. SR 15546 US/RS 136 - Par conséquent, sur les figures 48 et 49, lorsque la fonction décrite cidessus en référence à la figure 49 est achevée, les "lignes de contrôle" 200 (définies comme étant les "lignes en pointillés" 200 et les lignes de failles qui son fixées 170, 172, 174 de la figure 48) connectent les deux côtés de la ligne d'intersection de faille 170 à la limite 176, lesdites "lignes de contrôle" 200 connectent également les deux côtés de la ligne d'intersection de faille 172 à la limite 176, et lesdites "lignes de contrôle" 200 connectent également les deux côtés de la ligne d'intersection de faille 174 à la limite 176. Il faut garder à l'esprit le fait qu'une "ligne de contrôle" de la figure 48 est définie comme étant: la ligne en pointillés 200 et la ligne d'intersection de faille 170 (en utilisant la ligne d'intersection de faille 170 comme exemple). Toutefois, sur la figure 48, une ligne de contrôle 200 doit diviser une région en deux parties; c'est-à-dire qu'une ligne de contrôle de faille "i" 200 ne peut pas se terminer sur une limite de faille "j". Ceci signifie qu'un système de coordonnées peut être établi à l'intérieur de la région. Si on prend un point arbitraire (x, y) dans une région, on peut déterminer l'endroit o un point arbitraire est situé par rapport à une origine 202 sur la figure 48. C'est-à-dire qu'étant donné un point (x, y), on peut trouver d'o provient ce point amené depuis la limite par les vecteurs "i", ou d'o il provient après été amené de la limite par les vecteurs "j" (d'après un théorème, on sait que n'importe quel point doit être amené du haut d'une limite par les lignes "i" ou par la gauche par les lignes "j"). Ensuite, lorsqu'on connaît exactement SR 15546 US/RS la position du point arbitraire par rapport à l'origine 202 de la figure 48, on peut déterminer l'identité du triangle spécifique de la triangulation dans laquelle se trouve le point arbitraire. Ensuite, lorsqu'on connaît l'identité du triangle spécifique dans lequel est placé le point arbitraire, on peut déterminer les vecteurs directionnels voisins qui sont situés de façon adjacente au point arbitraire dans le triangle identifié. Ensuite, par interpolation, on peut déterminer le vecteur directionnel au niveau du point arbitraire. Etape 5 - Ajouter un quadrillage précis de lignes supplémentaires défini par les exigences des utilisateurs On demande à l'utilisateur combien il souhaite qu'il y ait de cellules ou "blocs de quadrillage" dans son quadrillage de simulation. S'il veut un quadrillage "50 par 50" (50 blocs de quadrillage sur un côté et 50 blocs de quadrillage sur l'autre côté), on élabore un quadrillage "50 par 50", avec un quadrillage de contrôle sous-jacent. S'il demande un quadrillage plus précis, par exemple "80 par 80", on ajoute un quadrillage précis de "lignes supplémentaires" défini par des points placés à intervalles réguliers sur les lignes de contrôle sous-jacentes. Ces chiffres, par exemple "50 par 50" sont répartis dans les régions fermées par la limite et les lignes de contrôle. Cette décision est prise sur la base des longueurs des côtés des sous-régions. Sur la figure 50, à l'intérieur de la limite 176, qui contient les failles "i" 170, 172 et la faille "j" SR 15546 US/RS 174, et les lignes de contrôle 200, on ajoute quelques "lignes supplémentaires" 204. Etape 6 - Placer à intervalles réguliers les lignes de contrôle et les lignes supplémentaires Sur la figure 51, l'utilisateur spécifie "Nx" et "Ny" entre des lignes de contrôle (c'est-à-dire que l'utilisateur indique combien de "lignes supplémentaires" il souhaite avoir entre les lignes de contrôle). Ce nombre de points est alors réparti, espacé régulièrement sur la longueur d'un arc le long des lignes de contrôle entre les intersections avec les autres lignes de contrôle ou avec la limite. Une fois que ceci a été effectué, les points du quadrillage restant à l'intérieur de chaque sous-région renfermée par les lignes de contrôle sont placés en utilisant une méthode de quadrillage adaptée à la limite standard (par exemple la méthode d'interpolation transfinie, la méthode de Winslow ou la méthode d'orthogonalité de zone de Knupp). En conséquence, on termine par des bandes ou des espaces espacés à intervalles réguliers entre des lignes de contrôle adjacentes 200 et entre des lignes supplémentaires adjacentes 204. Etape 7 - Avec uniformité Sur la figure 52, si l'utilisateur a besoin d'uniformité, on ne tente pas de placer à intervalles réguliers les lignes de contrôle 200 et les lignes supplémentaires 204 tel que cela a été indiqué ci- dessus à l'étape 6 en référence à la figure 51. C'est- à-dire qu'on ne divise pas la zone entre deux lignes de contrôle 200 et les lignes supplémentaires 204 avec des SR 15546 US/RS bandes de largeur égale; au lieu de cela, on règle la largeur entre les lignes de contrôle et les lignes supplémentaires, de telle sorte que, lorsque l'on se déplace à travers le quadrillage, ladite largeur change de façon uniforme. En résumé, le dispositif de quadrillage de surface structuré 68b de la figure 26 du dispositif de quadrillage structuré 68 des figures 14c, 16 et 24 de la présente invention sert à élaborer un "quadrillage de surface structuré" (également connu sous le nom "quadrillage de point de coin") en exécutant les étapes fonctionnelles suivantes: (1) Elaborer une triangulation qui absorbe une limite et un ensemble de lignes d'intersection de failles qui sont placées sur un horizon d'une formation terrestre; (2) Elaborer un champ de vecteur sur les noeuds de la triangulation, le champ de vecteur comportant un premier champ de vecteur comportant en outre des vecteurs directionnels et un deuxième champ de vecteur comportant en outre des vecteurs directionnels, le premier champ de vecteur ayant un premier sens de champ de vecteur qui correspond au sens des limites et des lignes d'intersection de failles qui forment une intersection avec les noeuds de la triangulation, le deuxième champ de vecteur ayant un deuxième sens de champ de vecteur qui est interpolé par rapport au premier sens du champ de vecteur et qui forme une intersection avec d'autres noeuds de la triangulation; (3) Elaborer un réseau de lignes de contrôle ayant un sens correspondant au sens des vecteurs directionnels dans les premier et deuxième champs de vecteurs, les lignes de contrôle comportant une ou SR 15546 US/RS plusieurs lignes d'intersection de failles et des extensions en pointillés qui forment une interconnexion avec les extrémités de chacune des lignes d'intersection de failles avec la limite; (4) Elaborer un réseau de lignes supplémentaires en plus du réseau de lignes de contrôle, ce qui permet de produit un "quadrillage de zone naïf"; et (5) Ensuite, lorsque le quadrillage de structure naïf a été élaboré, on post-traite le quadrillage de structure naïf de sorte que les lignes de contrôle et les lignes supplémentaires sont espacées à intervalles réguliers et réparties uniformément. L'invention a ainsi été décrite, et il apparaîtra évident que celle-ci peut être modifiée de différences façons. Ces variations ne doivent pas être considérées comme un éloignement de l'esprit et de la portée de l'invention et toutes ces modifications qui sembleront évidentes au spécialiste de la technique doivent être incluses dans la portée des revendications jointes. SR 15546 US/RS REVENDICATIONS 1. Appareil de quadrillage de simulation comportant un dispositif de quadrillage structuré, un horizon dans une formation terrestre comportant une ou plusieurs lignes d'intersection de failles, ledit dispositif de quadrillage structuré comprenant: un moyen renfermant une limite permettant de renfermer une limite autour desdites lignes d'intersection de failles sur ledit horizon, et un dispositif de quadrillage de surface structuré conçu pour élaborer un quadrillage de surface structuré sur ladite limite lorsque les dites lignes d'intersection de failles sont entourées par ladite limite par ledit moyen renfermant la limite, ledit quadrillage de surface structuré comportant deux lignes de flux qui forment l'intersection, lesdites lignes de flux qui forment l'intersection comportant une pluralité de noeuds, au moins certaines lignes de flux comportant des lignes de contrôle, lesdites lignes de contrôle comportant lesdites lignes d'intersection. 2. Appareil de quadrillage de simulation selon la revendication 1, dans lequel ledit dispositif de quadrillage structuré comprend en outre: un dispositif de quadrillage de bloc conçu pour connecter une pluralité de lignes de coordonnées à la pluralité de noeuds desdites lignes de flux qui forment l'intersection et de réduire ladite pluralité de lignes de coordonnées dans ladite formation terrestre jusqu'à ce que lesdites lignes de coordonnées forment l'intersection avec d'autres horizons dans ladite formation terrestre, ce qui permet de produire une SR 15546 US/RS pluralité de blocs de quadrillage d'un quadrillage de simulation structuré. 3. Appareil de quadrillage de simulation selon la revendication 2, dans lequel ledit dispositif de quadrillage de surface structuré comprend: un moyen permettant d'élaborer une triangulation et recouvrant ladite triangulation sur ladite limite et lesdites lignes d'intersection de failles renfermées par ladite limite, la triangulation présentant des noeuds. 4. Appareil de quadrillage de simulation selon la revendication 3, dans lequel ledit dispositif de quadrillage de surface structuré comprend en outre: un moyen pour élaborer un champ de vecteur comportant une pluralité de vecteurs directionnels sur les noeuds de ladite triangulation, ce qui permet de produire une triangulation de vecteurs, chacun des vecteurs directionnels du champ de vecteur est attribué à l'un des noeuds de la triangulation, chacun des vecteurs directionnels ayant un sens; et un moyen réagissant au sens de ladite pluralité de vecteurs directionnels sur lesdits noeuds de la triangulation pour élaborer un réseau de lignes de contrôle qui forment une intersection avec les côtés opposés de ladite limite, ce qui permet de produire un quadrillage préliminaire, chacune des lignes de contrôle étant formée lorsque les extrémités de l'une des lignes d'intersection de faille enfermées par ladite limite sont connectées auxdits côtés opposés de la dite limite. 5. Appareil de quadrillage de simulation selon la revendication 4, dans lequel ledit dispositif de quadrillage de surface structuré comprend en outre: SR 15546 US/RS un moyen pour ajouter des lignes supplémentaires audit quadrillage préliminaire en faisant interconnecter lesdites lignes supplémentaires avec les côtés opposés de ladite limite entre lesdites lignes de commande, ce qui permet de produire ledit quadrillage de surface structuré ayant une pluralité de noeuds qui recouvre ladite limite; et un moyen pour soit placer à intervalles réguliers des lignes de contrôle et des lignes supplémentaires, soit répartir uniformément lesdites lignes de contrôle et lesdites lignes supplémentaires dudit quadrillage de surface structure. 6. Appareil de quadrillage de simulation selon la revendication 2, dans lequel ledit appareil de quadrillage de simulation comporte en outre: un dispositif de mise à une plus grande échelle conçu pour recevoir ladite pluralité de blocs de quadrillage dudit dispositif de quadrillage de simulation structuré à partir du dudit dispositif de quadrillage de bloc et attribuer une pluralité d'informations de propriétés de blocs de quadrillage structuré plus précises à la pluralité respective de blocs de quadrillage dudit quadrillage de simulation structuré, grâce à quoi un simulateur reçoit ladite pluralité d'informations de propriétés de blocs de quadrillage plus précises associée à ladite pluralité de blocs de quadrillage du quadrillage de simulation structuré et génère une pluralité de résultats de simulation plus précis associée respectivement, à la pluralité de blocs de quadrillage. 7. Procédé de construction de quadrillage de simulation structuré conçu pour être placé sur une SR 15546 US/RS formation terrestre, un horizon contenu dans ladite formation terrestre comportant une pluralité de lignes d'intersection de failles, comprenant les étapes consistant à: (a) renfermer une limite autour des lignes d'intersection de failles de l'horizon; (b) élaborer une triangulation et recouvrir la triangulation sur ladite limite, ladite triangulation renfermant ladite limite et lesdites lignes d'intersection de failles à l'intérieur de ladite limite, la triangulation ayant des noeuds; et (c) élaborer un champ de vecteur ayant une pluralité de vecteurs directionnels sur le noeud de la triangulation, ce qui permet de produire une triangulation de vecteur, en attribuant l'un des vecteurs directionnels du champ de vecteur à chaque noeud de la triangulation, chacun des vecteurs directionnels ayant une direction. 8. Procédé selon la revendication 7, comprenant en outre les étapes consistant à: (d) en réponse au sens de la pluralité de vecteurs directionnels sur les noeuds de la triangulation, élaborer un réseau de lignes de contrôle qui forment une interconnexion avec les côtés opposés de la limite en interconnectant les extrémités de chacune des lignes d'intersection de la faille à l'intérieur de la limite avec les côtés opposés de la limite, ce qui permet de produire un quadrillage de surface qui recouvre la limite; (e) ajouter des lignes supplémentaires au quadrillage de surface en interconnectant les lignes supplémentaires avec les côtés opposés de la limite, ce qui permet de produire un quadrillage de surface SR 15546 US/RS structuré ayant une pluralité de noeuds à travers ladite limite, et (f) placer à intervalles réguliers les lignes de contrôle et les lignes supplémentaires de quadrillage de surface structuré ou bien répartir uniformément des lignes de contrôle et des lignes supplémentaires du quadrillage de surface structuré. 9. Procédé selon la revendication 8, comprenant en outre les étapes consistant à: (g) connecter une pluralité de lignes de coordonnées à la pluralité de noeuds du quadrillage de surface structuré de l'étape (e), et (h) réduire la pluralité de lignes de coordonnées à partir des noeuds du quadrillage de surface structuré et former une intersection avec lesdites lignes de coordonnées et un ou plusieurs autres horizons orientés vers le bas dans la formation terrestre, ce qui permet de produire ledit quadrillage de simulation structuré, ledit quadrillage de simulation structuré comportant une pluralité de blocs de quadrillage, chacun desdits blocs de quadrillage étant placé entre des horizons adjacents de la formation terrestre. 10. Appareil de quadrillage de simulation conçu pour générer un ensemble d'informations de propriétés de blocs de quadrillage plus précises conçu pour être utilisé par un simulateur, un horizon contenu dans une formation terrestre comportant une ou plusieurs lignes d'intersection de failles, comprenant: un dispositif de quadrillage structuré conçu pour générer un quadrillage de simulation structuré, ledit dispositif de quadrillage structuré comportant, SR 15546 US/RS un moyen renfermant une limite pour renfermer une limite autour desdites lignes d'intersection de failles sur ledit horizon, un dispositif de quadrillage de surface structuré conçu pour élaborer un quadrillage de surface structuré sur ladite limite une fois que lesdites lignes d'intersection de failles sont renfermées par ladite limite par ledit moyen renfermant une limite, ledit quadrillage de surface structuré comportant des lignes de flux formant une intersection, au moins certaines des lignes de flux comportant des lignes de contrôle, lesdites lignes de contrôle comportant lesdites lignes d'intersection de failles, lesdites lignes de flux qui forment l'intersection comportant une pluralité de noeuds; et un dispositif de quadrillage de bloc conçu pour connecter une pluralité de lignes de coordonnées à la pluralité de noeuds desdites lignes de flux formant une intersection et réduire ladite pluralité de lignes de coordonnées dans ladite formation terrestre jusqu'à ce que lesdites lignes de coordonnées forment une intersection avec d'autres horizons dans ladite formation terrestre, ce qui permet de produire une pluralité de blocs de quadrillage dudit quadrillage de simulation structuré. 11. Appareil de quadrillage de simulation selon la revendication 10, dans lequel ledit appareil de quadrillage de simulation comprend en outre: un dispositif de mise à une plus grande échelle conçu pour recevoir ladite pluralité de blocs de quadrillage dudit quadrillage de simulation structuré dudit dispositif de quadrillage de bloc et attribuer une pluralité d'informations de propriétés de blocs de SR 15546 US/RS quadrillage structuré plus précises à la pluralité de blocs de quadrillage dudit quadrillage de simulation structuré, grâce à quoi un simulateur reçoit ladite pluralité d'informations de propriétés de blocs de quadrillage structuré plus précises associée à ladite pluralité de blocs de quadrillage du quadrillage de simulation structuré et génère une pluralité de résultats de simulation plus précis associée, respectivement, à la pluralité de blocs de quadrillage. 12. Appareil de quadrillage de simulation selon la revendication 10, dans lequel ledit dispositif de quadrillage de surface structuré comprend en outre: un moyen pour élaborer une triangulation et recouvrir ladite triangulation sur ladite limite et les lignes d'intersection de failles renfermées par ladite limite, la triangulation ayant des noeuds. 13. Appareil de quadrillage de simulation selon la revendication 12, dans lequel ledit dispositif de quadrillage de surface structure comprend en outre: un moyen pour élaborer un champ de vecteur comportant une pluralité de vecteurs directionnels sur les noeuds de ladite triangulation, ce qui permet de produire une triangulation de vecteurs, chacun des vecteurs directionnels du champ de vecteurs étant attribué à l'un des noeuds de la triangulation, chacun des vecteurs directionnels ayant un sens; et un moyen réagissant au sens de ladite pluralité de vecteurs directionnels sur lesdits noeuds de la triangulation pour élaborer un réseau de lignes de contrôle qui forment une interconnexion avec les côtés opposés de ladite limite, ce qui permet de produire un quadrillage préliminaire, chacune desdites limites de SR 15546 US/RS contrôle étant formée lorsque les extrémités de l'une des lignes d'intersection de failles renfermée par ladite limite sont connectées auxdits côtés opposés de ladite limite. 14. Appareil de quadrillage de simulation selon la revendication 13, dans lequel ledit dispositif de quadrillage de surface structuré comprend en outre: un moyen pour ajouter des lignes supplémentaires audit quadrillage préliminaire en interconnectant lesdites lignes supplémentaires avec les côtés opposés de ladite limite entre lesdites lignes de contrôle, ce qui permet de produire ledit quadrillage de surface structuré ayant une pluralité de noeuds qui recouvrent ladite limite; et un moyen pour soit placer à intervalles réguliers lesdites lignes de contrôle et lesdites lignes supplémentaires, soit répartir uniformément lesdites lignes de contrôle et lesdites lignes supplémentaires dudit quadrillage de surface structuré. 15. Appareil de quadrillage de simulation selon la revendication 14, dans lequel ledit appareil de quadrillage de simulation comprend en outre: un dispositif de mise à une plus grande échelle conçu pour recevoir ladite pluralité de blocs de quadrillage dudit quadrillage de simulation structuré dudit dispositif de quadrillage de bloc et attribuer une pluralité d'informations de propriétés de blocs de quadrillage structuré plus précises à la pluralité de blocs de quadrillage dudit quadrillage de simulation structuré, grâce à quoi un simulateur reçoit ladite pluralité d'informations de propriétés de blocs de quadrillage structuré plus précise associée à ladite SR 15546 US/RS pluralité de blocs de quadrillage du quadrillage de simulation structuré et génère une pluralité de résultats de simulation plus précis associée, respectivement, à la pluralité de blocs de quadrillage. 16. Dispositif pour stocker des instructions destinées à élaborer un quadrillage de surface structuré conçu pour être placé sur une formation terrestre, dans une limite sur un horizon, lesdites instructions devant être exécutées par un processeur d'ordinateur et étant aptes à exécuter un procédé comprenant les opérations consistant à: (a) ren.ermer une limite autour des lignes d'intersection de failles de l'horizon; (b) élaborer une triangulation et recouvrir la triangulation sur ladite limite, ladite triangulation renfermant ladite limite et lesdites lignes d'intersection de failles à l'intérieur de ladite limite, la triangulation ayant des noeuds. 17>' Dispositif selon la revendication 16, dans lequel lesdites instructions conduisent à exécuter un procédé comprenant en outre les opérations suivantes: (c) élaborer un champ de vecteur ayant une pluralité de vecteurs directionnels sur le noeud de la triangulation, ce qui permet de produire une triangulation de vecteur, en attribuant l'un des vecteurs directionnels du champ de vecteur à chaque noeud de la triangulation, chacun des vecteurs directionnels ayant une direction; et (d) en réponse au sens de la pluralité de vecteurs directionnels sur les noeuds de la triangulation, élaborer un réseau de lignes de contrôle qui forment une interconnexion avec les côtés opposés de la limite en interconnectant les extrémités de SR 15546 US/RS chacune des lignes d'intersection de la faille à l'intérieur de la limite avec les côtés opposés de la limite, ce qui permet de produire un quadrillage de surface qui recouvre la limite. 18. Dispositif selon la revendication 17, dans lequel lesdites instructions conduisent à exécuter un procédé comprenant en outre les opérations suivantes: (e) ajouter des lignes supplémentaires au quadrillage de surface en interconnectant les lignes supplémentaires avec les côtés opposés de la limite, ce qui permet de -produire un quadrillage de surface structuré ayant- une pluralité de noeuds à travers ladite limite, et (f) placer à intervalles réguliers les lignes de contrôle et les lignes supplémentaires de quadrillage de surface structuré ou bien répartir uniformément des lignes de contrôle et des lignes supplémentaires du quadrillage de surface structuré. -S6 SR 15546 US/RS