Brevet FR3012899A1 — publié le 20150508

Contexte [0001] La modélisation de bassins est utilisée pour évaluer la génération, la migration et l'accumulation d'hydrocarbures dans les bassins sédimentaires. En modélisation de bassins, la perméabilité peut couvrir plusieurs ordres de grandeur. Les grès fortement perméables présents dans les roches réservoirs ont des perméabilités pouvant atteindre des valeurs de l'ordre de plusieurs Darcy (D), alors que pour du schiste tassé, la perméabilité tombe à quelques nanodarcy (nD), c.-à-d. approximativement neuf ordres de grandeur de différence. En utilisant ces perméabilités, les échelles de temps sur lesquelles a lieu l'écoulement d'hydrocarbures peuvent être estimées. Ces différences de perméabilité engendrent des difficultés dans la modélisation de la migration. [0002] Un traitement exact de l'écoulement d'hydrocarbures à l'intérieur d'un bassin sédimentaire à l'aide de la loi de Darcy sans approximation supplémentaire impose des pas de temps de l'ordre d'une année ou plus petits. Toutefois, il se peut que cela ne soit pas réalisable, car les temps d'exécution de simulation augmentent lorsque la durée du pas de temps est diminuée. Par conséquent, des approximations supplémentaires peuvent être employées. [0003] Pour traiter numériquement la loi de Darcy, la perméabilité peut être limitée de sorte que l'échelle de temps des processus d'écoulement soit au moins de l'ordre de la durée des pas de temps, limitant ainsi l'écoulement à l'intérieur de roches fortement perméables. En utilisant l'algorithme d'écoulement de Darcy, la migration dans des roches faiblement perméables est prise en considération de manière plus exacte, surtout lorsque les effets de chronologie de l'expulsion des roches sources sont traités avec exactitude. Cependant, par comparaison à d'autres procédés de migration, des pas de temps plus petits sont utilisés, ce qui introduit là encore des coûts de calcul élevés. Résumé [0004] Un procédé de modélisation d'une migration de fluides d'hydrocarbures est décrit. Le procédé peut comprendre les étapes consistant à décomposer, par l'utilisation d'un processeur, un domaine d'un modèle de bassin pour donner une première pluralité de mailles présentant une perméabilité inférieure à une valeur seuil et une deuxième pluralité de mailles présentant une perméabilité supérieure à la valeur seuil. Une première analyse de migration peut être effectuée pour analyser une migration relativement lente en commençant dans la première pluralité de mailles. Une deuxième analyse de migration peut être effectuée pour analyser une migration relativement plus rapide en commençant dans la deuxième pluralité de mailles. [0005] Un support lisible par ordinateur est également décrit. Sur le support lisible par ordinateur peuvent être conservées des instructions qui, lorsqu'elles sont exécutées par un processeur, sont configurées pour faire en sorte que le processeur effectue des opérations. Les opérations peuvent comprendre l'étape consistant à décomposer un domaine d'un modèle de bassin pour donner une première pluralité de mailles présentant une perméabilité inférieure à une valeur seuil et une deuxième pluralité de mailles présentant une perméabilité supérieure à la valeur seuil. Les opérations peuvent également comprendre l'étape consistant à réaliser un algorithme d'écoulement de Darcy pour analyser une migration relativement lente en commençant dans la première pluralité de mailles. Les opérations peuvent comprendre en outre l'étape consistant à réaliser un algorithme de percolation d'invasion pour analyser une migration relativement plus rapide en commençant dans la deuxième pluralité de mailles. [0006] Un système informatique est également décrit. Le système informatique peut comprendre un processeur et un système de mémoire comprenant un ou plusieurs supports lisibles par ordinateur sur lesquels sont conservées des instructions qui, lorsqu'elles sont exécutées par le processeur, sont configurées pour faire en sorte que le système informatique effectue des opérations. Les opérations peuvent comprendre l'étape consistant à décomposer un domaine d'un modèle de bassin pour donner une première pluralité de mailles présentant une perméabilité inférieure à une valeur seuil et une deuxième pluralité de mailles présentant une perméabilité supérieure à la valeur seuil. Les opérations peuvent également comprendre l'étape consistant à réaliser un algorithme d'écoulement de Darcy pour analyser une migration relativement lente en commençant dans la première pluralité de mailles. Les opérations peuvent comprendre en outre l'étape consistant à réaliser un algorithme de percolation d'invasion pour analyser une migration relativement plus rapide en commençant dans la deuxième pluralité de mailles.30 [0007] On aura conscience que le résumé qui précède vise simplement à introduire un sous-ensemble de l'objet exposé ci-dessous et n'est donc pas limitatif. Description succincte des dessins [0008] Les dessins joints, qui sont incorporés dans et font partie de cette spécification, illustrent des modes de réalisation des présents enseignements et, conjointement à la description, servent à expliquer les principes des présents enseignements. Dans les figures: [0009] La Figure 1 illustre une vue d'un domaine d'un modèle de bassin, indiquant la perméabilité, selon un mode de réalisation. [0010] La Figure 2 illustre une vue d'une décomposition du domaine du modèle de bassin de la Figure 1, selon un mode de réalisation. [0011] La Figure 3 illustre une vue d'une partie agrandie du modèle de bassin de la Figure 2, selon un mode de réalisation. [0012] La Figure 4 illustre une vue d'une décomposition d'un domaine d'un modèle de bassin présentant une perforation de bouchon, selon un mode de réalisation. [0013] La Figure 5 illustre un organigramme d'un procédé de modélisation d'une migration de fluides d'hydrocarbures, selon un mode de réalisation. [0014] La Figure 6 illustre une vue schématique d'un système informatique, selon un mode de 25 réalisation. Description détaillée [0015] La description détaillée qui suit fait référence aux dessins joints. Chaque fois que possible, les mêmes numéros de repères sont utilisés dans les dessins et la description qui suit 30 pour faire référence à des pièces identiques ou similaires. Bien que plusieurs modes de réalisation et traits caractéristiques de la présente invention soient décrits ici, des modifications, des adaptations et d'autres mises en oeuvre sont possibles, sans s'écarter de l'esprit et du champ de la présente invention. [0016] Des modes de réalisation de la présente invention concernent un procédé de migration, qui peut être basé sur une décomposition en domaines d'un modèle de bassin, et des systèmes destinés à réaliser de tels procédés. La Figure 1 illustre un modèle 100 de bassin doté d'une grande étendue de perméabilité parmi des couches, selon un mode de réalisation. Comme représenté, le modèle 100 de bassin comprend sept couches 101 à 107; toutefois, comme on peut s'en rendre compte, le nombre de couches dans le modèle 100 peut varier en fonction du domaine faisant l'objet de la cartographie, les dimensions du domaine, etc. [0017] La perméabilité dans le modèle 100 est mesurée en [log(mD)], car plus de 12 ordres de grandeur sont couverts. Dans le modèle 100 de bassin illustré, la première couche 101, dite "inférieure", peut présenter une perméabilité d'environ -6 à -4 [log(mD)]. La deuxième couche 102 peut présenter une perméabilité d'environ -1 à 3 [log(mD)]. La troisième couche 103 peut présenter une perméabilité d'environ -7 à -3 [log(mD)]. La quatrième couche 104 peut présenter une perméabilité d'environ 1 à 4 [log(mD)]. La cinquième couche 105 peut présenter une perméabilité d'environ -7 à -2 [log(mD)]. La sixième couche 106 peut présenter une perméabilité d'environ 2 à 5 [log(mD)]. La septième couche 107, dite "supérieure", peut présenter une perméabilité d'environ -2 à -1 [log(mD)]. Comme on peut s'en rendre compte, les couches 101 à 107 et la perméabilité de chaque couche 101 à 107 sont illustratives et peuvent varier d'un modèle à un autre. [0018] La Figure 2 illustre le modèle 100 de bassin de la Figure 1, mais décomposé (par ex. partitionné ou divisé) en deux parties 201, 202 sur la base de la perméabilité, selon un mode de réalisation. Une valeur seuil pour la perméabilité peut être sélectionnée. La valeur seuil peut être d'environ -5 log(mD) à environ -3 log(mD), d'environ -3 log(mD) à environ -1 log(mD), d'environ -1 log(mD) à environ 1 log(mD), d'environ 1 log(mD) à environ 3 log(mD), ou d'environ 3 log(mD) à environ 5 log(mD). Sur la Figure 2, la valeur seuil sélectionnée est d'environ 0 [log(mD)]. Une fois la valeur seuil sélectionnée, le modèle 100 de bassin peut être décomposé en une première partie 201 présentant une perméabilité inférieure à la valeur seuil (par ex. les couches 101, 103, 105, 107) et une deuxième partie 202 présentant une perméabilité supérieure à la valeur seuil (par ex. les couches 102, 104, 106). Dans d'autres modes de réalisation, un seuil de n'importe quelle autre valeur peut être sélectionné, sans limitation, en fonction de divers facteurs relatifs au modèle 100 de bassin en cours d'analyse. [0019] La partie 201 présentant la perméabilité la plus basse peut être constituée de roches sources, de schistes, ou de siltites schisteuses. On entend par roches sources des roches à partir desquelles un fluide d'hydrocarbure a été généré ou peut être généré. On entend par schistes des roches à grain fin comprenant un mélange de minéraux argileux et d'autres minéraux comme du quartz, de la calcite, etc. La partie 202 présentant la perméabilité la plus élevée peut être constituée de roches réservoirs, de grès ou de calcaires perméables (par ex. des carbonates poreux). On entend par roches réservoirs des roches qui présentent une porosité suffisante pour contenir des fluides d'hydrocarbures. [0020] La Figure 3 illustre une vue d'une partie agrandie du modèle 100 de bassin après la décomposition, selon un mode de réalisation. Chaque couche (par ex. 103, 104) dans le modèle 100 de bassin peut être fragmentée en une pluralité de mailles de taille similaire. Comme représenté sur la Figure 3, chaque maille s'étend sur environ 0,2 km (suivant l'axe X) et environ 120 m (suivant l'axe Y); toutefois, la taille des mailles peut varier dans une ou les deux directions. [0021] Comme représenté dans le coin inférieur gauche de la Figure 3, un fluide d'hydrocarbure est généré dans la partie 201 de plus faible perméabilité (par ex. la couche 103) du modèle 100 de bassin. Le fluide d'hydrocarbure s'écoule ou migre vers le haut à travers la couche 103 et des points d'injection (c.-à-d. l'interface entre les couches 103, 104) jusqu'à la partie 202 de plus forte perméabilité (par ex. la couche 104), comme indiqué par les flèches 210. Le fluide d'hydrocarbure s'écoule vers le haut et sur le côté droit de la partie 202 de plus forte perméabilité (par ex. la couche 104), comme indiqué par les carrés 212. Ceci peut être appelé "chapelet de migration". Le fluide d'hydrocarbure s'accumule ensuite dans le coin supérieur droit, comme indiqué par les cercles 214. [0022] L'écoulement d'hydrocarbures à l'intérieur du modèle 100 de bassin représenté sur la Figure 3 peut être évalué par pas de temps. Les pas de temps peuvent être de même durée ou ils peuvent varier. Les pas de temps peuvent s'échelonner d'un minimum d'environ 10 ans, environ 100 ans, environ 1000 ans, ou environ 10 000 ans à un maximum d'environ 100 000 ans, environ 500 000 ans, environ 1 000 000 ans, ou davantage. Pour chaque pas de temps, une première analyse de migration, comme un algorithme d'écoulement de Darcy, peut être effectuée sur chacune des mailles dans le modèle 100 de bassin (c.-à-d. à la fois les parties 201 de plus faible perméabilité et les parties 202 de plus forte perméabilité). Dans un autre mode de réalisation, l'algorithme d'écoulement de Darcy peut être réalisé sur un sous-ensemble ciblé des mailles. Par exemple, l'algorithme d'écoulement de Darcy peut être réalisé sur les mailles dans les parties 201 de plus faible perméabilité. [0023] L'algorithme d'écoulement de Darcy servant à calculer l'écoulement d'un fluide visqueux (par ex. un liquide ou un gaz d'hydrocarbure) à travers des milieux poreux peut être décrit par l'équation suivante : = où f) dénote la vitesse de l'écoulement, y dénote le tenseur de mobilité qui est donné par le quotient de la perméabilité de la matrice rocheuse, k, et la viscosité du liquide / du gaz, y, comme = -k. La force motrice Vu dénote un gradient, qui est orienté dans la direction de la diminution la plus abrupte du champ de potentiel u. Le gradient Vu est une formulation mathématique de différences de potentiel de pression sur des distances de taille infinitésimalement petites dans des directions spatiales. [0024] Lorsque plus d'un composant est pris en considération, des effets supplémentaires peuvent être pris en compte, comme des pressions capillaires, des perméabilités relatives, etc. Des pressions capillaires peuvent jouer un rôle dans les systèmes eau-pétrole-gaz et introduire des effets non linéaires. En général, il se peut que des équations caractérisant de tels modèles qui incluent la prise en considération d'effets multiples ne puissent pas être résolues analytiquement. Par conséquent, un traitement numérique peut être utilisé, qui est basé sur une discrétisation à la fois de l'espace (c.-à-d. que le modèle est décomposé en un ensemble de mailles spatiales de propriétés égales) et du temps (introduction de pas de temps). Diverses méthodes numériques peuvent ainsi être utilisées afin de traiter l'écoulement d'hydrocarbures à travers des milieux poreux, comme des méthodes d'éléments finis ou des méthodes de volumes finis. [0025] La taille spatiale des mailles décrivant le processus d'écoulement est donnée par l'échelle à laquelle les propriétés physiques (par ex. la perméabilité, la viscosité et / ou la porosité) varient. L'échelle à laquelle les saturations varient de façon significative doit être plus grande que la taille des mailles de la grille. Des tailles typiques de mailles sont de l'ordre de 1 km horizontalement et 100 m verticalement. Une durée maximale raisonnable du pas de temps peut être estimée par le débit obtenu par la loi de Darcy. Lorsqu'un algorithme d'écoulement de Darcy pur est envisagé, on peut choisir un pas de temps tel que la composition en hydrocarbures d'une maille n'est pas significativement modifiée dans la limite d'un pas de temps. Si les pas de temps sont trop grands, on peut limiter le débit dans les parties de plus forte perméabilité pour stabiliser l'algorithme, par exemple en réduisant la perméabilité dans ces zones. [0026] Après que l'algorithme d'écoulement de Darcy a été réalisé, une deuxième analyse de migration, comme un algorithme de percolation d'invasion, peut être effectuée sur les mailles des parties 202 de plus forte perméabilité (par ex. la couche 104) pour le premier pas. Comme l'algorithme d'écoulement de Darcy, l'algorithme de percolation d'invasion peut être utilisé pour calculer l'écoulement d'un fluide visqueux (par ex. un liquide ou un gazeux d'hydrocarbure) à travers des milieux poreux. [0027] La modélisation de l'écoulement des fluides par percolation d'invasion peut négliger la chronologie pendant le mouvement du fluide. Au lieu de cela, un bilan des forces de flottabilité d'hydrocarbure avec des seuils de pression capillaire, (par ex. au dessous d'un bouchon) peut être recherché. On peut supposer que dans des roches perméables, un équilibre sera atteint au bout d'un temps suffisant. Les échelles de temps géologiques dans le cadre de la modélisation de bassins, qui couvrent l'évolution globale de bassins sédimentaires, dépassent habituellement de tels laps de temps. L'algorithme de percolation d'invasion est décrit plus en détail dans: T. Hantschel & A. Kauerauf (2009): "Fundamentals of Basin and Petroleum Systems Modeling," Chapitre 6.8, Springer, qui est incorporé ici par référence dans son intégralité. [0028] L'algorithme de percolation d'invasion est ainsi basé sur deux approximations. Premièrement, tout effet visqueux des hydrocarbures peut être négligé et ainsi, la migration peut être vue comme se produisant instantanément, c.-à-d. que tous effets de chronologie de la migration peuvent ne pas être pris en compte. Deuxièmement, il suppose que la migration a lieu dans des parcours de migration, également appelés "chapelets." Par exemple, les parcours de migration peuvent être des parcours ou chapelets unidimensionnels. [0029] L'algorithme de percolation d'invasion peut prendre en considération avec exactitude la migration rapide dans des roches fortement perméables. Comme la perméabilité est élevée, tous les processus d'écoulement à l'intérieur de ces roches peuvent être considérés comme instantanés aux échelles de temps géologiques. Ainsi, le fait de négliger tous les effets visqueux est une bonne approximation pour les roches fortement perméables. En outre, la résolution spatiale peut être accrue, donnant une représentation exacte de géométries complexes, comprenant des structures de failles complexes. On peut supposer que des perforations, c.-à-d. des pertes de fluide d'hydrocarbure à travers des couches de bouchon, ont lieu dans des parcours de migration en chapelets, ce qui peut être décrit avec plus d'exactitude par l'algorithme de percolation d'invasion que par l'algorithme d'écoulement de Darcy. [0030] L'algorithme de percolation d'invasion peut être réalisé sur un maillage plus fin (c.-à-d. des mailles plus petites) que l'algorithme d'écoulement de Darcy. Plus particulièrement, l'algorithme de percolation d'invasion peut être traité sur un maillage d'une résolution supérieure d'environ 10 à environ 100 fois à celle du maillage de l'algorithme d'écoulement de Darcy. Par exemple, les mailles du modèle 100 de bassin pour l'algorithme de percolation d'invasion peuvent être d'environ 25 m, d'environ 50 m ou d'environ 100 m (suivant l'axe X horizontal) et d'environ 5 m, d'environ 10 m ou d'environ 20 m (suivant l'axe Y vertical). [0031] Comme observé plus haut, l'écoulement analysé par l'algorithme de percolation d'invasion dans la partie 202 de plus forte perméabilité peut se produire instantanément (c.-à-d. être approché comme étant instantané) du point de vue de l'algorithme d'écoulement de Darcy, car la perméabilité dans la partie 202 de plus forte perméabilité est supérieure à la valeur seuil. Ainsi, l'échelle de temps sur laquelle la migration a lieu peut être plus courte que la durée des pas de temps employés par l'analyse numérique du modèle 100 de bassin. [0032] Comme la migration d'hydrocarbures dans les parties de plus forte perméabilité 202 est traitée par l'algorithme de percolation d'invasion, il est possible d'augmenter la durée des pas de temps de l'algorithme d'écoulement de Darcy, qui sont dans ce cas limités uniquement par l'écoulement dans les parties 201 de plus faible perméabilité. La durée des pas de temps peut être accrue d'un facteur d'environ 1 jusqu'à un facteur d'environ 1000, voire au-delà, permettant ainsi d'utiliser des pas de temps de 1000 ans jusqu'à 10 million d'années et plus. Cela pourrait se traduire par des améliorations significatives de performances d'exécution. [0033] La Figure 4 illustre une vue d'une décomposition d'un domaine d'un modèle 400 de bassin présentant une ou plusieurs perforations 406 de bouchon, selon un mode de réalisation. Comme sur la Figure 3, les flèches 410 indiquent une migration d'une partie 401 de plus faible perméabilité vers une partie 402 de plus forte perméabilité. Cette migration peut être effectuée par l'algorithme d'écoulement de Darcy. La résolution de l'écoulement calculé par l'algorithme de percolation d'invasion est élevée et des accumulations peuvent apparaître comme un nuage de points reliés. [0034] Les perforations 406 de bouchon peuvent être mieux décrites comme des parcours de migration en chapelets. Par conséquent, de tels parcours peuvent être pris en considération par l'algorithme de percolation d'invasion. A cette fin, le procédé peut comprendre l'étape consistant à imposer une condition supplémentaire pour l'algorithme d'écoulement de Darcy. Plus précisément, l'écoulement de parties 402 de plus forte perméabilité jusque dans des parties 401 de plus faible perméabilité (c.-à-d. du réservoir jusque dans la couche de bouchon), comme représenté sur la Figure 4, peut être empêché. Des accumulations de fluides d'hydrocarbures dans les parties 401 de plus faible perméabilité, créées par des perforations 406 de bouchon, peuvent alors être analysées ou réinjectées en utilisant l'algorithme de percolation d'invasion lors de pas de temps ultérieurs.30 [0035] La Figure 5 illustre un organigramme d'un procédé 500 de modélisation d'une migration de fluides d'hydrocarbures dans un bassin sédimentaire, selon un mode de réalisation. Le procédé 500 peut comprendre l'étape consistant à décomposer un modèle de bassin pour donner des première et deuxième pluralités de mailles, comme représenté en 502. Le modèle de bassin peut être décomposé en sélectionnant d'abord une valeur seuil de perméabilité. La première pluralité de mailles peut présenter une perméabilité inférieure à la valeur seuil et la deuxième pluralité de mailles peut présenter une perméabilité supérieure à la valeur seuil. [0036] Un algorithme d'écoulement de Darcy peut ensuite être réalisé pour analyser une migration relativement lente en commençant dans la première pluralité de mailles, comme représenté en 504. L'algorithme d'écoulement de Darcy peut également être réalisé en commençant dans la deuxième pluralité de mailles. L'algorithme d'écoulement de Darcy peut être réalisé selon le choix de pas de temps de la première pluralité de mailles uniquement (par ex. des pas de temps moins nombreux et plus longs). Une fois que l'algorithme d'écoulement de Darcy s'est achevé, un algorithme de percolation d'invasion peut être réalisé pour analyser une migration relativement plus rapide en commençant dans la deuxième pluralité de mailles (c.-à-d. présentant la perméabilité supérieure à la valeur seuil), comme représenté en 506. En outre, des structures de réservoirs complexes comprenant des répartitions complexes de faciès et des failles peuvent être modélisées en utilisant l'algorithme de percolation d'invasion. Le procédé 500 peut favoriser un rendement et des performances améliorés. Par comparaison à un algorithme d'écoulement de Darcy pur, les pas de temps peuvent être plus grands (les réservoirs fortement perméables n'étant pas pris en compte pour l'extension de l'algorithme d'écoulement de Darcy), diminuant ainsi les exigences de calcul, tout en conservant une précision suffisante. [0037] Les modes de réalisation de l'invention peuvent également comprendre un ou plusieurs systèmes destinés à mettre en oeuvre un ou plusieurs modes de réalisation du procédé de la présente invention. La Figure 6 illustre une vue schématique d'un tel système 600 informatique ou à processeur, selon un mode de réalisation. Le système 600 à processeur peut comprendre un ou plusieurs processeurs 602 dont les configurations de coeurs (y compris à coeurs multiples) et les fréquences d'horloge varient. Le ou les processeurs 602 peuvent être utilisables pour exécuter des instructions, appliquer une logique, etc. On aura conscience que ces fonctions peuvent être assurées par des processeurs multiples ou des coeurs multiples sur une même puce, fonctionnant en parallèle et / ou liés de façon à communiquer entre eux. [0038] Le système 600 à processeur peut également comprendre un système de mémoire, qui peut être ou comprendre un ou plusieurs dispositifs de mémoire et/des ou supports 604 lisibles par ordinateur (par ex. des supports non transitoires lisibles par ordinateur) dont les dimensions physiques, l'accessibilité, les capacités de stockage, etc. varient, comme des unités Flash, des disques durs, des disques, une mémoire vive, etc., servant à stocker des données, comme des images, des fichiers et des instructions de programme en vue de leur exécution par le processeur 602. Dans un mode de réalisation, les supports 604 lisibles par ordinateur peuvent conserver des instructions qui, lorsqu'elles sont exécutées par le processeur 602, sont configurées pour faire en sorte que le système 600 à processeur effectue des opérations. Par exemple, l'exécution de ces instructions peut faire en sorte que le système 600 à processeur mette en oeuvre une ou plusieurs parties et / ou modes de réalisation du ou des procédés 500 décrits plus haut. [0039] Le système 600 à processeur peut également comprendre une ou plusieurs interfaces 606 de réseau. Les interfaces 606 de réseau peuvent comprendre n'importe quels matériels, applications et / ou autres logiciels. En conséquence, les interfaces 606 de réseau peuvent comprendre des adaptateurs Ethernet, des émetteurs-récepteurs sans fil, des composants de réseau à interface PCI et / ou série, destinés à communiquer sur des supports filaires ou sans fil en utilisant des protocoles comme Ethernet, Ethernet sans fil, etc. [0040] Le système 600 à processeur peut comprendre en outre une ou plusieurs interfaces 608 de périphériques, servant à une communication avec un écran d'affichage, un projecteur, des claviers, des souris, des pavés tactiles, des capteurs, d'autres types de périphériques d'entrée et / ou de sortie, et / ou d'autres similaires. Dans certaines réalisations, les composants du système 600 à processeur ne sont pas nécessairement englobés à l'intérieur d'un boîtier unique ou même situés à proximité immédiate les uns des autres, mais dans d'autres réalisations, les composants et / ou d'autres peuvent être placés dans un boîtier unique.30 [0041] Les supports 604 lisibles par ordinateur peuvent être physiquement ou logiquement disposés ou configurés pour conserver des données sur un ou plusieurs dispositifs 610 de stockage. Le dispositif 610 de stockage peut comprendre un ou plusieurs systèmes de fichiers ou bases de données dans un quelconque format approprié. Le dispositif 610 de stockage peut également comprendre un ou plusieurs programmes logiciels 612, qui peuvent contenir des instructions interprétables ou exécutables destinées à réaliser un ou plusieurs des processus décrits. Sur demande du processeur 602, un ou plusieurs des programmes logiciels 612, ou une partie de ceux-ci, peuvent être chargés des dispositifs 610 de stockage vers les supports 604 lisibles par ordinateur en vue de leur exécution par le processeur 602. [0042] Les personnes qualifiées dans le métier auront conscience que les composants décrits ci-dessus ne constituent qu'un exemple d'une configuration matérielle, car le système 600 à processeur peut comprendre n'importe quel type de composants matériels, y compris un quelconque microprogramme ou logiciel éventuel d'accompagnement nécessaire, pour réaliser les mises en oeuvre décrites. Le système 600 à processeur peut également être mis en oeuvre en partie ou en totalité par des composants de circuits électroniques ou des processeurs, tels que des circuits intégrés à application spécifique (ASIC) ou des circuits intégrés prédiffusés programmables (FPGA). [0043] La description qui précède de la présente invention, ainsi que les modes de réalisation et exemples qui lui sont associés, a été présentée uniquement à des fins d'illustration. Elle n'est pas exhaustive et ne limite pas la présente invention à la forme précise décrite. Les personnes qualifiées dans le métier auront conscience, d'après la description qui précède, que des modifications et des variations sont possibles à la lumière des enseignements ci-dessus ou peuvent être acquises à partir de la mise en pratique des modes de réalisation décrits. [0044] Par exemple, les mêmes techniques décrites ici en faisant référence au système 600 à processeur peuvent être utilisées pour exécuter des programmes selon des instructions reçues en provenance d'un autre programme ou d'un tout autre système à processeur. De façon analogue, des commandes peuvent être reçues, exécutées, et leurs sorties renvoyées entièrement à l'intérieur du traitement et / ou de la mémoire du système 600 à processeur. En conséquence, ni un terminal de commande à interface visuelle ni même un quelconque terminal ne sont strictement nécessaires pour concrétiser les modes de réalisation décrits. [0045] De même, il n'est pas nécessaire que les étapes décrites soient effectuées dans le même ordre que celui évoqué ou avec le même degré de séparation. Diverses étapes peuvent être éliminées, répétées, combinées ou divisées, le cas échéant, pour atteindre des objectifs ou des améliorations identiques ou similaires. En conséquence, la présente invention n'est pas limitée aux modes de réalisation décrits ci-dessus, mais est plutôt définie par les revendications jointes au vu de toute leur gamme d'équivalents. En outre, dans la description ci-dessus et dans les revendications ci-dessous, sauf indication contraire, le terme "exécuter" et ses variantes doivent être interprétés comme se rapportant à toute utilisation de code de programme ou d'instructions sur un dispositif, qu'ils soient compilés, interprétés ou exécutés en utilisant d'autres techniques. REVENDICATIONS1 Système informatique comportant: un processeur; et un système de mémoire comportant un ou plusieurs supports lisibles par ordinateur sur lesquels sont conservées des instructions qui, lorsqu'elles sont exécutées par le processeur, sont configurées pour faire en sorte que le système informatique effectue des opérations, les opérations comportant les étapes consistant à: décomposer un domaine d'un modèle de bassin pour donner une première pluralité de mailles présentant une perméabilité inférieure à une valeur seuil et une deuxième pluralité de mailles présentant une perméabilité supérieure à la valeur seuil; réaliser un algorithme d'écoulement de Darcy pour analyser une migration relativement lente en commençant dans la première pluralité de mailles; et réaliser un algorithme de percolation d'invasion pour analyser une migration relativement plus rapide en commençant dans la deuxième pluralité de mailles. 2. Système informatique selon la revendication 1, dans lequel les opérations comportent en outre l'étape consistant à réaliser l'algorithme d'écoulement de Darcy en commençant dans la deuxième pluralité de mailles. 3. Système informatique selon la revendication 1 ou la revendication 2, dans lequel les opérations comportent en outre l'étape consistant à réaliser l'algorithme de percolation d'invasion après que l'algorithme d'écoulement de Darcy est réalisé. 4. Système informatique selon l'une quelconque des revendications 1 à 3, dans lequel la deuxième pluralité de mailles comporte une perforation de bouchon à travers une couche comprenant la première pluralité de mailles. 5. Support lisible par ordinateur sur lequel sont conservées des instructions qui, lorsqu'elles sont exécutées par un processeur, sont configurées pour faire en sorte que le processeur effectue des opérations, les opérations comportant les étapes consistant à:décomposer un domaine d'un modèle de bassin pour donner une première pluralité de mailles présentant une perméabilité inférieure à une valeur seuil et une deuxième pluralité de mailles présentant une perméabilité supérieure à la valeur seuil; réaliser un algorithme d'écoulement de Darcy pour analyser une migration relativement lente en commençant dans la première pluralité de mailles; et réaliser un algorithme de percolation d'invasion pour analyser une migration relativement plus rapide en commençant dans la deuxième pluralité de mailles.