FR3036820A1 - Modelisation de la saturation et permeabilite de reservoir de champ petrolifere - Google Patents

Abstract

Système et procédé pour corriger des courbes de pression capillaire, comprenant la création (610) d'une courbe de pression capillaire en utilisant de multiples tangentes hyperboliques liées, la détermination (620) d'une pression de coupure de correction de fermeture de la courbe de pression capillaire, et la correction (630) de la courbe de pression capillaire. La correction peut comprendre la normalisation (640) de la courbe de pression capillaire et l'extrapolation (660) de la courbe de pression capillaire.

Description

MODÉLISATION DE SATURATION ET DE PERMÉABILITÉ DE RÉSERVOIR DE CHAMP PÉTROLIFÈRE Arrière-plan [0001] Pour créer des modèles précis de réservoir de champ pétrolifère, une saturation en eau et d'hydrocarbures peuvent être prédites à un point donné dans le réservoir de champ pétrolifère. [0002] Des données de saturation peuvent être disponibles à l'échelle correcte, où elles peuvent être dérivées avec précision à partir de données de diagraphie de puits pétrophysiques en utilisant différent(e)s flux de traitement et normes. Toutefois, il peut s'avérer souhaitable de calculer une saturation à l'échelle du réservoir, où peu de propriétés du réservoir sont connues. Dans pareils cas, un modèle de saturation peut être obtenu en utilisant une fonction de hauteur de saturation. Toutefois, des modèles de saturation peuvent être basés sur des fonctions de hauteur de saturation pour des systèmes de gorge à pore unique, ou si une modélisation de gorge à pores multiples est possible, sur des modèles instables qui dépendent du nombre de points de données utilisés et de la sélection des meilleurs intervalles d'ajustement. Résumé [0003] Des modes de réalisation de la divulgation peuvent fournir un système de calcul, un support lisible par ordinateur non transitoire et un procédé pour corriger une courbe de pression capillaire.

Par exemple, le procédé comprend la création d'une courbe de pression capillaire en utilisant de multiples tangentes hyperboliques liées, et la détermination d'une pression de coupure de correction de fermeture de la courbe de pression capillaire. Le procédé peut comprendre en outre la correction de la courbe de pression capillaire, où la correction de la courbe de pression capillaire comprend la normalisation de la courbe de pression capillaire ou l'extrapolation de la courbe de pression capillaire. [0004] Dans un autre mode de réalisation, le procédé peut comprendre en outre la détermination d'une pression d'entrée pour la courbe de pression capillaire après la correction de la courbe de pression capillaire. [0005] Dans un autre mode de réalisation, la détermination de la pression de coupure de correction de fermeture de la courbe de pression capillaire comprend la détermination d'un minimum local de la dérivée seconde de la courbe de pression capillaire, et la pression de coupure de correction de fermeture correspond au minimum local où la pression capillaire est la plus basse. 3036820 2 [0006] Dans un autre mode de réalisation, la normalisation de la courbe de pression capillaire comprend le réglage de la pression de coupure de correction de fermeture afm de correspondre à une saturation de phase de mouillage de 100 %, et la pression d'entrée correspond à la pression de coupure de correction de fermeture après le réglage de la pression de coupure de correction de 5 fermeture. [0007] Dans un autre mode de réalisation, l'extrapolation de la courbe de pression capillaire comprend la détermination d'une première pression de coupure dérivée, l'exécution d'une analyse de régression entre la première pression de coupure dérivée et la pression de coupure de correction de fermeture, et l'extrapolation de la courbe de pression capillaire au point de saturation de phase 10 de mouillage de 100 %, dans lequel la première pression de coupure dérivée correspond au premier minimum local dérivé où la pression capillaire est la plus basse, et dans lequel la pression d'entrée correspond à une pression au point de saturation de phase de mouillage de 100 %. [0008] Dans un autre mode de réalisation, la correction de la courbe de pression capillaire comprend le rognage de la courbe de pression capillaire afin d'éliminer les données de pression 15 capillaire en dessous de la pression de coupure de correction de fermeture. [0009] Dans un autre mode de réalisation, la courbe de pression capillaire est défmie par l'équation suivante: f(P, an, Wn, tn) = al + aN Eri=i (an+i - ).tanh(w.(P -ta)) 20 avec les contraintes suivantes: wn > 0, tin E [1, /V] n, N E N an+1 < an, Vn E [1, N -1] n,N E N 25 où P représente une transformée logarithmique d'une pression capillaire normalisée et N représente le nombre de tangentes hyperboliques. [0010] Dans un autre mode de réalisation, le support lisible par ordinateur non transitoire stocke des instructions qui, lorsqu'elles sont exécutées par un ou plusieurs processeur(s) d'un système de 30 calcul, amènent le système de calcul à exécuter des opérations. Par exemple, les opérations peuvent comprendre la création d'une courbe de pression capillaire en utilisant de multiples tangentes hyperboliques liées, et la détermination d'une pression de coupure de correction de fermeture de la courbe de pression capillaire. Les opérations peuvent comprendre en outre la correction de la courbe de pression capillaire, et la détermination d'une pression d'entrée pour la courbe de 3036820 3 pression capillaire après la correction de la courbe de pression capillaire, où la correction de la courbe de pression capillaire comprend la normalisation de la courbe de pression capillaire ou l'extrapolation de la courbe de pression capillaire. [0011] Dans un autre mode de réalisation, le système de calcul peut comprendre un ou plusieurs 5 processeur(s), et un système de mémoire comprenant un ou plusieurs support(s) lisible(s) par ordinateur non transitoire(s) pour stocker des instructions qui, lorsqu'elles sont exécutées par ledit/lesdits un ou plusieurs processeur(s) du système de calcul, amènent le système de calcul à exécuter des opérations. Par exemple, les opérations peuvent comprendre la création d'une courbe de pression capillaire en utilisant de multiples tangentes hyperboliques liées, et la détermination 10 d'une pression de coupure de correction de fermeture de la courbe de pression capillaire. Les opérations peuvent comprendre en outre la correction de la courbe de pression capillaire, et la détermination d'une pression d'entrée pour la courbe de pression capillaire après la correction de la courbe de pression capillaire, où la correction de la courbe de pression capillaire comprend la normalisation de la courbe de pression capillaire ou l'extrapolation de la courbe de pression 15 capillaire. [0012] Ce résumé est fourni pour introduire une sélection de concepts qui sont décrits en détail ci-dessous dans la description détaillée. Ce résumé n'a pas pour objectif d'identifier des caractéristiques clés ou essentielles de l'objet revendiqué, et ne doit pas non plus être utilisé comme aide pour limiter la portée de l'objet revendiqué.

20 Brève description des dessins [0013] Les dessins annexés, qui sont incorporés dans et constituent une partie du présent document, illustrent des modes de réalisation des présents enseignements. Ces aspects et ces avantages, et/ou d'autres, dans les modes de réalisation de la divulgation, apparaîtront clairement et 25 seront plus facilement appréciés à la lecture de la description qui suit des différents modes de réalisation, effectuée en conjonction avec les dessins annexés, dans lesquels: [0014] la Figure 1 illustre un exemple d'un système qui comprend différents composants de gestion afin de gérer différents aspects d'un environnement géologique selon un mode de réalisation; [0015] la Figure 2 est un organigramme d'un procédé pour modéliser une saturation dans un 30 réservoir selon un mode de réalisation; [0016] la Figure 3 montre un modèle de tangentes hyperboliques dans un système de pression capillaire et de saturation en eau selon un mode de réalisation; 3036820 4 [0017] la Figure 4 montre un modèle de tangentes hyperboliques dans un système de pression capillaire et de saturation en eau selon un mode de réalisation; [0018] la Figure 5 montre un modèle de tangentes hyperboliques dans un système de pression capillaire et de saturation en eau selon un mode de réalisation; 5 [0019] la Figure 6 montre des données de pression capillaire à partir d'un système de gorge à pores multiples selon un mode de réalisation; [0020] la Figure 7 montre un ajustement de courbe à des données de pression capillaire selon un mode de réalisation; [0021] la Figure 8 montre des tangentes hyperboliques qui correspondent à des gorges à pores 10 selon un mode de réalisation; [0022] la Figure 9 est un organigramme d'un procédé pour corriger une courbe de pression capillaire selon un mode de réalisation; [0023] la Figure 10 montre une courbe de pression capillaire selon un mode de réalisation; [0024] la Figure 11 montre une dérivée seconde d'une courbe de pression capillaire selon un mode 15 de réalisation; [0025] la Figure 12 montre une courbe de pression capillaire rognée selon un mode de réalisation; [0026] la Figure 13 montre une courbe de pression capillaire normalisée selon un mode de réalisation; [0027] la Figure 14 montre une dérivée première d'une courbe de pression capillaire selon un mode 20 de réalisation; [0028] la Figure 15 montre une courbe de pression capillaire extrapolée selon un mode de réalisation; et [0029] la Figure 16 est une vue schématique d'un système de calcul selon un mode de réalisation. [0030] On notera que certains détails des dessins ont été simplifiés et sont dessinés de manière à 25 faciliter la compréhension des présents enseignements plutôt qu'à maintenir une précision structurelle, un détail et une échelle stricts. Ces dessins/figures sont fournis dans un but explicatif et non limitatif. Description détaillée 30 [0031] Nous allons maintenant nous référer en détail aux différents modes de réalisation de la présente divulgation, dont des exemples sont illustrés dans les dessins et les figures annexés. Les modes de réalisation sont décrits ci-dessous dans le but d'assurer une compréhension plus complète 3036820 5 des composants, procédés et appareils divulgués ici. Tous les exemples fournis ont un but illustratif, et non limitatif. Il apparaîtra cependant clairement à l'homme du métier que l'invention peut être mise en pratique sans ces détails spécifiques. Dans d'autres exemples, des procédés, procédures, composants, circuits, et réseaux bien connus n'ont pas été décrits en détail de manière à 5 ne pas obscurcir inutilement des aspects des modes de réalisation. [0032] D'un bout à l'autre du fascicule et des revendications, les termes suivants ont les significations explicitement associées ici, sauf si le contexte indique clairement autre chose. Les expressions "dans certains modes de réalisation" et "dans un mode de réalisation", telles qu'elles sont employées ici, ne font pas nécessairement référence au(x) même(s) mode(s) de réalisation, 10 bien que cela puisse être le cas. En outre, les expressions "dans un autre mode de réalisation" et "dans certains autres modes de réalisation", telles qu'elles sont employées ici, ne font pas nécessairement référence à un mode de réalisation différent, bien que cela puisse être le cas. Comme cela est décrit ci-dessous, différents modes de réalisation peuvent être facilement combinés, sans sortir de la portée ou de l'esprit de la présente divulgation. 15 [0033] Tel qu'il est employé ici, le terme "ou" est un opérateur inclusif, et est équivalent au terme "et/ou," sauf si le contexte indique clairement autre chose. Le terme "basé sur" n'est pas exclusif et permet un fondement sur des facteurs supplémentaires, comme cela est décrit, sauf si le contexte indique clairement autre chose. Dans le fascicule, l'énumération de "au moins un parmi A, B, et C," comprend des modes de réalisation contenant A, B, ou C, de multiples exemples de A, B, ou C, ou 20 des combinaisons de A/B, A/C, BIC, A/B/B/ B/B/C, A/B/C, etc. En outre, d'un bout à l'autre du fascicule, la signification de "un", "une" et "le/la" inclut plusieurs références. La signification de "dans" englobe "dans" et "sur". [0034] On comprendra également que, bien que les termes premier, deuxième, etc. puissent être employés ici pour décrire des éléments différents, ces éléments ne seront pas limités par ces 25 termes. Ces termes sont employés dans le but de distinguer un élément d'un autre. Par exemple, un premier objet ou une première étape pourrait être appelé(e) un deuxième objet ou une deuxième étape et, de façon similaire, un deuxième objet ou une deuxième étape pourrait être appelé(e) un premier objet ou une première étape, sans sortir de la portée de l'invention. Le premier objet ou la première étape, et le deuxième objet ou la deuxième étape, sont tous les deux des objets ou des 30 étapes, respectivement, mais ils ne doivent pas être considérés comme étant le même objet ou la même étape. On comprendra en outre que les termes "comprend," "comprenant," "inclut" et/ou "incluant," lorsqu'ils sont employés dans le présent fascicule, spécifient la présence de 3036820 6 caractéristiques, entiers, étapes, opérations, éléments, et/ou composants évoqués, mais n'excluent pas la présence ou l'ajout d'un(e) ou de plusieurs autre(s) caractéristique(s), entiers(s), étape(s), opération(s), élément(s), composant(s) et/ou groupe(s) de ceux-ci. En outre, tel qu'il est employé ici, le terme "si" peut être interprété comme signifiant "lorsque" ou "lors" ou "en réponse à la 5 détermination" ou "en réponse à la détection," en fonction du contexte. [0035] Lorsqu'il est fait référence ici à une plage numérique de valeurs, ces plages doivent être interprétées comme englobant chacun et tous les nombres et/ou fractions entre les plages minimum et maximum indiquées. Par exemple, une plage de 0,5 % à 6 % inclurait expressément les valeurs intermédiaires de 0,6 %, 0,7 % et 0,9 %, jusqu'à et y compris 5,95 %, 5,97 % et 5,99 %. Le même 10 principe s'applique à chaque autre propriété numérique et/ou plage élémentaire précisée ici, sauf si le contexte indique clairement autre chose. [0036] L'attention est maintenant portée sur des procédures, procédés, techniques et flux de traitement de traitement qui sont présentés selon certains modes de réalisation. Certaines opérations dans les procédures, procédés, techniques et flux de traitement de traitement divulgués ici peuvent 15 être combinés et/ou l'ordre de certaines opérations peut être changé. [0037] La Figure 1 montre un exemple d'un système 100 qui comprend différents composants de gestion 110 pour gérer différents aspects d'un environnement géologique 150 (par exemple, un environnement qui comprend un bassin sédimentaire, un réservoir 151, une ou plusieurs faille(s) 153-1, un ou plusieurs corps géologique(s) 153-2, etc.). Par exemple, les composants de gestion 20 110 peuvent permettre une gestion directe ou indirecte de la détection, du forage, de l'injection, de l'extraction, etc., par rapport à l'environnement géologique 150. Ensuite, des informations supplémentaires relatives à l'environnement géologique 150 peuvent devenir disponibles en retour 160 (par exemple, optionnellement comme une entrée dans une ou plusieurs des composants de gestion 110). 25 [0038] Dans l'exemple qui est montré dans la Figure 1, les composants de gestion 110 comprennent un composant de données sismiques 112, un composant d'informations supplémentaires 114 (par exemple, des données de puits/diagraphie), un composant de traitement 116, un composant de simulation 120, un composant d'attributs 130, un composant d'analyse/visualisation 142 et un composant de flux de traitement 144. Lors du fonctionnement, des 30 données sismiques et d'autres informations fournies par les composants 112 et 114 peuvent être entrées dans le composant de simulation 120. 3036820 7 [0039] Dans un exemple de mode de réalisation, le composant de simulation 120 peut être fondé sur des entités 122. Les entités 122 peuvent comprendre des entités de terre ou des objets géologiques tels que des puits, des surfaces, des corps, des réservoirs, etc. Dans le système 100, les entités 122 peuvent comprendre des représentations virtuelles d'entités physiques réelles qui sont 5 reconstituées à des fins de simulation. Les entités 122 peuvent comprendre des entités basées sur des données acquises par l'intermédiaire d'une détection, d'une observation, etc. (par exemple, les données sismiques 112 et les autres informations 114). Une entité peut être caractérisée par une ou plusieurs propriété(s) (par exemple, une entité de grille de pilier géométrique d'un modèle terrestre peut être caractérisée par une propriété de porosité). Ces propriétés peuvent représenter un(e) ou 10 plusieurs mesure(s) (par exemple, des données acquises), calcul(s), etc. [0040] Dans un exemple de mode de réalisation, le composant de simulation 120 peut opérer en conjonction avec une infrastructure logicielle telle qu'une infrastructure orientée objet. Dans une telle infrastructure, les entités peuvent comprendre des entités basées sur des classes prédéfinies afin de faciliter la modélisation et la simulation. Un exemple disponible dans le commerce d'une 15 infrastructure basé sur un objet est l'infrastructure MICROSOFT® .NET® (Redmond, Washington), qui fournit un ensemble de classes d'objets extensibles. Dans l'infrastructure .NET®, une classe d'objets englobe un module de code réutilisable et de structures de données associées. Les classes d'objet peuvent être utilisées pour instancier des instances d'objet à utiliser par un programme, un script, etc. Par exemple, des classes de trous de forage peuvent défmir des objets pour représenter 20 des trous de forage basés sur des données de puits. [0041] Dans l'exemple qui est montré dans la Figure 1, le composant de simulation 120 peut traiter des informations pour correspondre à un ou plusieurs attribut(s) spécifié(s) par le composant d'attributs 130, qui peut comprendre une bibliothèque d'attributs. Ce traitement peut être effectué avant l'entrée dans le composant de simulation 120 (par exemple, considérons le composant de 25 traitement 116). A titre d'exemple, le composant de simulation 120 peut exécuter des opérations sur des informations d'entrée basées sur un ou plusieurs attribut(s) spécifié(s) par le composant d'attributs 130. Dans un exemple de mode de réalisation, le composant de simulation 120 peut construire un ou plusieurs modèle(s) de l'environnement géologique 150, sur lequel (lesquels) il est possible de se baser pour simuler un comportement de l'environnement géologique 150 (par 30 exemple, sensible à une ou plusieurs action(s), qu'elles soient naturelles ou artificielles). Dans l'exemple qui est montré dans la Figure 1, le composant d'analyse/visualisation 142 peut permettre une interaction avec un modèle ou des résultats basés sur un modèle (par exemple, des résultats 3036820 8 d'une simulation, etc.). A titre d'exemple, une sortie du composant de simulation 120 peut être entrée dans un ou plusieurs autres flux de traitement, comme cela est indiqué par un composant de flux de traitement 144. [0042] A titre d'exemple, le composant de simulation 120 peut comprendre un(e) ou plusieurs 5 caractéristique(s) d'un simulateur tel que le simulateur de réservoir ECLIPSETm (Schlumberger Limited, Houston, Texas), le simulateur de réservoir INTERSECTTm (Schlumberger Limited, Houston, Texas), etc. A titre d'exemple, un composant de simulation, un simulateur, etc. peuvent comprendre des caractéristiques pour mettre en oeuvre une ou plusieurs techniques sans maillages (par exemple, pour résoudre une ou plusieurs équation(s), etc.). A titre d'exemple, un réservoir ou 10 des réservoirs peu(ven)t être simulé(s) par rapport à une ou plusieurs technique(s) de récupération améliorée(s) (par exemple, considérons un traitement thermique tel que SAGD, etc.). [0043] Dans un exemple de mode de réalisation, les composants de gestion 110 peuvent comprendre des caractéristiques d'une infrastructure disponible dans le commerce tel que l'infrastructure logicielle de simulation sismique PETREL® (Schlumberger Limited, Houston, 15 Texas). L'infrastructure PETREL® fournit des composants qui permettent une optimisation des opérations d'exploration et de développement. L'infrastructure PETREL® comprend des composants de logiciel de simulation sismique qui peuvent sortir des informations à utiliser pour augmenter la performance du réservoir, par exemple en améliorant la productivité de l'équipe active. En utilisant une infrastructure de ce type, différents professionnels (par exemple, des 20 géophysiciens, des géologues et des ingénieurs de réservoir) peuvent développer des flux de traitement collaboratifs et intégrer des opérations à des procédés rationnels. Une telle infrastructure peut être considérée comme une application et peut être considérée comme une application dirigée par les données (par exemple, dans laquelle des données sont entrées à des fins de modélisation, de simulation, etc.). 25 [0044] Dans un exemple de mode de réalisation, différents aspects du composant de gestion 110 peuvent comprendre des modules additionnels ou des modules d'extension qui opèrent selon des spécifications d'un environnement d'infrastructure. Par exemple, un environnement d'infrastructure disponible dans le commerce commercialisé sous le nom d'environnement d'infrastructure OCEAN® (Schlumberger Limited, Houston, Texas) permet d'intégrer des modules additionnels (ou 30 des modules d'extension) dans un flux de traitement d'infrastructure PETREL®. L'environnement d'infrastructure OCEANe influence des outils .NET® (Microsoft Corporation, Redmond, Washington) et procure des interfaces stables et conviviales pour un développement efficace. Dans 3036820 9 un exemple de mode de réalisation, des composants différentes peuvent être mises en oeuvre comme des modules additionnels (ou des modules d'extension) qui se conforment à et opèrent suivant des spécifications d'un environnement d'infrastructure (par exemple, selon des spécifications d'interface de programmation d'application (API), etc.). 5 [0045] La Figure 1 montre également un exemple d'une infrastructure 170 qui comprend une couche de simulation de modèle 180 de concert avec une couche de services d'infrastructure 190, une couche centrale d'infrastructure 195 et une couche de modules 175. L'infrastructure 170 peut comprendre l'infrastructure OCEAN® disponible dans le commerce où la couche de simulation de modèle 180 est le pack logiciel centré sur le modèle PETREL® disponible dans le commerce qui 10 contient les applications d'infrastructure OCEAN®. Dans un exemple de mode de réalisation, le logiciel PETREL® peut être considéré comme une application dirigée par les données. Le logiciel PETREL® peut comprendre une infrastructure pour la construction et la visualisation de modèles. [0046] A titre d'exemple, une infrastructure peut comprendre des caractéristiques pour mettre en oeuvre une ou plusieurs technique(s) de génération de maillages. Par exemple, une infrastructure 15 peut comprendre un composant d'entrée pour recevoir des informations provenant de l'interprétation de données sismiques, un ou plusieurs attribut(s) basé(s) au moins en partie sur des données sismiques, des données de diagraphie, des données d'image, etc. Une infrastructure de ce type peut comprendre un composant de génération de maillages qui traite des informations d'entrée, optiormellement en conjonction avec d'autres informations, pour générer un maillage. 20 [0047] Dans l'exemple qui est montré dans la Figure 1, la couche de simulation de modèle 180 peut fournir des objets de domaine 182, agir en tant que source de données 184, alimenter un rendu 186 et alimenter différentes interfaces d'utilisateur 188. Le rendu 186 peut fournir un environnement graphique dans lequel des applications peuvent afficher leurs données, alors que les interfaces d'utilisateur 188 peuvent apporter l'apparence et la convivialité pour des composants d'interface 25 d'utilisateur. [0048] A titre d'exemple, les objets de domaine 182 peuvent comprendre des objets d'entité, des objets de propriété et optionnellement d'autres objets. Les objets d'entité peuvent être utilisés pour représenter de façon géométrique des puits, des surfaces, des corps, des réservoirs, etc., alors que les objets de propriété peuvent être utilisés pour fournir des valeurs de propriété aussi bien que des 30 paramètres de version et d'affichage de données. Par exemple, un objet d'entité peut représenter un puits où un objet de propriété peut fournir des informations de diagraphie aussi bien que des 3036820 10 informations de version et des informations d'affichage (par exemple, pour afficher le puits comme faisant partie d'un modèle). [0049] Dans l'exemple qui est montré dans la Figure 1, des données peuvent être stockées dans une ou plusieurs source(s) de données (ou magasin(s) de données, d'une manière générale des 5 dispositifs de stockage de données physiques), qui peuvent se trouver sur le même site physique ou sur des sites physiques différents et qui sont accessibles par l'intermédiaire d'un ou de plusieurs réseau(x). La couche de simulation de modèle 180 peut être configurée de manière à modéliser des projets. De ce fait, un projet particulier peut être stocké, où des informations de projet stockées peuvent comprendre des entrées, des modèles, des résultats et des cas. Donc, lors de 10 l'accomplissement d'une session de modélisation, un utilisateur peut stocker un projet. Il est possible ultérieurement d'accéder au projet et de le récupérer en utilisant la couche de simulation de modèle 180, qui est capable de recréer des exemples des objets de domaine pertinents. [0050] Dans l'exemple qui est montré dans la Figure 1, l'environnement géologique 150 peut comprendre des couches (par exemple, une stratification) qui comprennent un réservoir 151 et une 15 ou plusieurs autre(s) caractéristique(s), telles que la faille 153-1, le corps géologique 153-2, etc. A titre d'exemple, l'environnement géologique 150 peut être équipé de l'un quelconque parmi une variété de capteurs, de détecteurs, d'actionneurs, etc. Par exemple, l'équipement 152 peut comprendre un circuit de communication pour recevoir et transmettre des informations par rapport à un ou plusieurs réseau(x) 155. Ces informations peuvent comprendre des informations associées 20 à un équipement de forage 154, qui peut être un équipement pour acquérir des informations, aider à récupérer des ressources, etc. Un autre équipement 156 peut être situé à distance d'un site de puits et comprendre un circuit de captation, de détection, d'émission ou autre. Un tel équipement peut comprendre un circuit de stockage et de communication pour stocker et communiquer des données, des instructions, etc. A titre d'exemple, un ou plusieurs satellite(s) peu(ven)t être prévu(s) des fms 25 de communication, d'acquisition de données, etc. Par exemple, la Figure 1 montre un satellite en communication avec le réseau 155 qui peut être configuré pour des communications, notons que le satellite peut comprendre un circuit destiné à une imagerie (par exemple, spatiale, spectrale, temporelle, radiométrique, etc.). [0051] La Figure 1 montre également l'environnement géologique 150 comme comprenant 30 optionnellement un équipement 157 et 158 associé à un puits qui présente une partie sensiblement horizontale qui peut intersecter une ou plusieurs fracture(s) 159. Par exemple, considérons un puits dans une formation de schiste pouvant comporter des fractures naturelles, des fractures artificielles 3036820 11 (par exemple, des fractures hydrauliques) ou une combinaison de fractures naturelles et artificielles. A titre d'exemple, un puits peut être foré pour un réservoir qui est latéralement étendu. Dans un exemple de ce type, des variations latérales des propriétés, des contraintes, etc. peuvent exister, où une évaluation de ces variations peut faciliter la planification, les opérations, etc. pour 5 développer un réservoir latéralement étendu (par exemple, par l'intermédiaire d'une fracturation, d'une injection, d'une extraction, etc.). A titre d'exemple, l'équipement 157 et/ou 158 peut comprendre des composants, un système, des systèmes, etc. pour la fracturation, la détection sismique, l'analyse de données sismiques, l'évaluation d'une ou de plusieurs fracture(s), etc. [0052] Comme cela est mentionné, le système 100 peut être utilisé pour exécuter un ou plusieurs 10 flux de traitement. Un flux de traitement peut être un procédé qui comprend un certain nombre d'étapes de traitement. Une étape de traitement peut opérer sur des données, par exemple, pour créer de nouvelles données, pour mettre à jour des données existantes, etc. A titre d'exemple, elle peut opérer sur une ou plusieurs entrée(s) et générer un ou plusieurs résultat(s), par exemple, basé(s) sur un ou plusieurs algorithme(s). A titre d'exemple, un système peut comprendre un 15 éditeur de flux de traitement pour la création, l'édition, l'exécution, etc. d'un flux de traitement. Dans un exemple de ce type, l'éditeur de flux de traitement peut assurer la sélection d'une ou plusieurs étape(s) de travail prédéfinie(s), d'une ou de plusieurs étapes de traitement personnalisée(s), etc. A titre d'exemple, un flux de traitement peut être un flux de traitement qui peut être mis en oeuvre dans le logiciel PETREL®, par exemple, qui opère sur des données 20 sismiques, un ou plusieurs attribut(s) sismique(s), etc. A titre d'exemple, un flux de traitement peut être un procédé qui peut être mis en oeuvre dans l'infrastructure OCEAN®. A titre d'exemple, un flux de traitement peut comprendre une ou plusieurs étape(s) de travail qui accède(nt) à un module tel qu'un module d'extension (par exemple, un code externe exécutable, etc.). [0053] Comme cela est décrit ci-dessus, le système 100 peut être utilisé pour simuler ou modéliser 25 un environnement géologique 150 et/ou un réservoir 151. Les modèles de réservoir se basent souvent sur des données de saturation comme composant. Dans certains modes de réalisation, le système 100 peut être basé sur un modèle de saturation comme composant du modèle de réservoir 151. [0054] La Figure 2 est un organigramme d'un procédé 200 pour modéliser une saturation dans un 30 réservoir, selon un mode de réalisation. Comme cela est illustré dans la Figure 2, le procédé 200 peut démarrer par l'obtention de données pétrophysiques lors de l'opération 210. Par exemple, lors de l'opération 210, les données pétrophysiques en provenance du réservoir peuvent être collectées 3036820 12 ou reçues. Les données pétrophysiques peuvent comprendre des données de pression capillaire et des données de perméabilité. Dans certains modes de réalisation, les données pétrophysiques peuvent également comprendre des données de porosité, de hauteur au-dessus du niveau d'eau libre et des données de type de roche. 5 [0055] A l'opération 220, un nombre de gorges à pores peut être déterminé à partir des données pétrophysiques obtenues. Par exemple, un nombre de gorges à pores peut être déterminé à partir des données de pression capillaire obtenues. Dans d'autres modes de réalisation, le nombre de gorges à pores dans le système peut être prédéterminé. [0056] Une fois que le nombre de gorges à pores a été déterminé, un ensemble de tangentes 10 hyperboliques dont le nombre est égal au nombre de gorges à pores peut être établi lors de l'opération 230. [0057] A l'opération 240, l'ensemble de tangentes hyperboliques peut être utilisé pour créer une courbe pour approcher les données pétrophysiques obtenues et défmir un ensemble de paramètres de tangente hyperbolique. Par exemple, l'ensemble de tangentes hyperboliques peut être utilisé 15 pour créer une courbe pour approcher les données de pression capillaire obtenues et défmir un ensemble de paramètres de tangente hyperbolique associées à ladite courbe. [0058] Dans certains modes de réalisation, la courbe créée pour approcher les données de pression capillaire peut être corrigée afin de tenir compte de la manière avec laquelle les données de pression capillaire ont été obtenues. Par exemple, la courbe peut être corrigée pour des effets de 20 conformité ou de fermeture. Le procédé 200 peut exécuter automatiquement une correction de fermeture sur l'ensemble de tangentes hyperboliques utilisées pour créer une courbe pour approcher les données de pression capillaire obtenues, et peut également déterminer la pression d'entrée qui correspond à la courbe de pression capillaire. [0059] Une fois que les paramètres de tangente hyperbolique ont été définis, au moins un 25 paramètre de tangente hyperbolique peut être combiné avec les données pétrophysiques obtenues afin de définir des dépendances pour une fonction de hauteur de saturation lors de l'opération 250. Par exemple, au moins un paramètre de tangente hyperbolique peut être combiné avec les données de perméabilité obtenues afin de définir une dépendance de perméabilité pour certains des paramètres qui définissent une fonction de hauteur de saturation. 30 [0060] A l'opération 260, la fonction de hauteur de saturation peut être combinée avec des données pétrophysiques afin de modéliser une saturation dans le réservoir. Par exemple, la saturation en eau et d'hydrocarbures dans un réservoir peut être calculée à partir de la fonction de hauteur de 3036820 13 saturation en utilisant des données de perméabilité, des données de porosité, et une hauteur au-dessus du niveau d'eau libre. Dans certains modes de réalisation, la fonction de hauteur de saturation peut également être combinée avec des données de type de roche. Par exemple, la fonction de hauteur de saturation peut être limitée à un seul type de roche ou un seul type de roche 5 peut être supposé pour le modèle de réservoir. [0061] A l'opération 270, le modèle de saturation peut être affiché. Par exemple, lors de l'opération 270, le modèle de saturation ou des changements du modèle de saturation peu(ven)t être affiché(s). Dans d'autres modes de réalisation, le modèle de saturation peut être affiché comme une partie du modèle de réservoir plus grand. 10 [0062] Comme cela est décrit ci-dessus, un modèle de données de saturation peut être utilisé pour prédire une saturation en eau et d'hydrocarbures à un point donné dans un réservoir de champ pétrolifère. Par exemple, un modèle de données de saturation peut être créé en utilisant des propriétés de réservoir telles que la perméabilité, la porosité, la hauteur au-dessus du niveau d'eau libre, et une fonction de hauteur de saturation. Dans certains modes de réalisation, des données de 15 porosité, de perméabilité et de type de roche peuvent être obtenues à partir de données sismiques et/ou de données de puits. De façon similaire, la fonction de hauteur de saturation peut être une fonction des données de pression capillaire, de saturation en eau et de perméabilité. Dans certains modes de réalisation, les données pétrophysiques pour ces propriétés de champ pétrolifère sont obtenues à partir d'une analyse d'échantillons de bouchon de carotte représentatifs du réservoir de 20 champ pétrolifère. [0063] Tel qu'il est employé ici, le terme "pression capillaire" fait référence à la différence dans les forces capillaires générées par deux ou plus de deux fluides immiscibles à l'intérieur de vides d'une roche. Les données de pression capillaire peuvent être mesurées par l'intermédiaire d'une expérimentation ou peuvent être reçues dans le modèle. Par exemple, la pression capillaire peut 25 être mesurée par l'intermédiaire d'expériences par plaque poreuse, centrifugation ou injection de mercure. Les données de pression capillaire peuvent comprendre la mesure d'une saturation à un niveau de pression et/ou de hauteur différent. Dans certains modes de réalisation, un registre de données de pression capillaire de laboratoire par rapport à une saturation de phase de mouillage ou à une saturation de phase de non mouillage est obtenu et est utilisé pour établir la fonction de 30 hauteur de saturation. Dans un autre mode de réalisation, les données de pression capillaire obtenues par expérimentation sont normalisées avant que les données de pression capillaire soient utilisées pour établir la fonction de hauteur de saturation. La normalisation peut permettre 3036820 14 l'utilisation de la fonction de hauteur de saturation avec un réservoir avec différents systèmes de fluide, tels que gaz/eau, pétrole/eau, et pétrole/eau/gaz. Dans un mode de réalisation, les données de pression capillaire mesurées sont représentatives de la pression capillaire du réservoir de champ pétrolifère ou d'une partie de celui-ci. Par exemple, des données de pression capillaire en termes de 5 hauteur peuvent représenter une épaisseur maximum du réservoir à modéliser. [0064] Tel qu'il est employé ici, le terme "saturation en eau" fait référence à une partie d'une porosité de substrat remplie d'eau. Dans un mode de réalisation, des données de saturation en eau peuvent être obtenues par expérimentation. Par exemple, la saturation en eau peut être obtenue à partir d'expériences sur la pression capillaire: la saturation de phase de non mouillage (en cas 10 d'injection de mercure) peut être calculée comme étant le volume occupé par la phase de non mouillage (mesure du volume injecté pendant l'expérience) sur le volume total de pores. Dans certains modes de réalisation, les données de saturation en eau sont normalisées. Dans un mode de réalisation, les données de saturation en eau sont représentatives de la saturation en eau du réservoir de champ pétrolifère ou d'une partie de celui-ci. 15 [0065] Tel qu'il est employé ci, le terme "perméabilité" fait référence à la capacité d'un substrat à transmettre un fluide. Dans un mode de réalisation, les données de perméabilité peuvent être obtenues par expérimentation. Par exemple, des données de perméabilité peuvent être dérivées à partir de pressions mesurées avant d'entrer un échantillon de substrat et après la sortie du substrat en utilisant un fluide de viscosité connue. Dans le cas d'un gaz, des corrections, telles qu'une 20 correction pour l'effet de Klinkenberg, peuvent être incluses. Dans un mode de réalisation, les données de perméabilité mesurées sont représentatives de la perméabilité du réservoir de champ pétrolifère ou d'une partie de celui-ci. [0066] Dans un mode de réalisation, la fonction de hauteur de saturation est basée sur deux équations pour approcher des données de pression capillaire mesurées à partir du réservoir: une 25 première équation qui résout un ensemble de paramètres inconnus en utilisant des données de pression capillaire mesurées, et une deuxième équation utilisant les paramètres inconnus résolus pour appliquer un ensemble de tangentes hyperboliques afin d'approcher des données de pression capillaire obtenues à partir d'un système de gorge à pore unique ou à pores multiples. Dans un mode de réalisation, ces équations ajustent des données de pression capillaire mesurées à partir du 30 réservoir en utilisant un procédé des moindres carrés non linéaire contraints. Dans un autre mode de réalisation, ces équations ajustent la pression capillaire et les données de saturation mesurées à partir du réservoir en utilisant un procédé des moindres carrés non linéaire contraints. 3036820 15 [0067] Par exemple, une première équation (Équation 1) peut utiliser un ensemble M de données mesurées de saturation en eau et de pression capillaire. Dans un mode de réalisation, les données de saturation en eau et de pression capillaire sont obtenues par analyse et expérimentation basées sur des échantillons de bouchon de carotte provenant du réservoir. Dans un autre mode de réalisation, 5 les données de saturation en eau et de pression capillaire sont normalisées, et la pression capillaire normalisée est transformée en logarithme de la pression capillaire avant son incorporation dans l'équation 1. [0068] Dans un mode de réalisation, l'équation 1 utilise l'ensemble M de données mesurées de saturation en eau et de pression capillaire dans un procédé des moindres carrés non linéaire pour 10 trouver des paramètres inconnus (an, wn, tu) d'un modèle qui minimise une erreur E entre les données et un modèle de pression capillaire f. Dans un mode de réalisation, la première équation correspond à l'équation suivante: Équation 1: E = Ei,_1(Srneas, - f (Pmeasi, an tn))2 15 où Smeas et Pmeas représentent les données de saturation en eau et de pression capillaire et an, wn, tu sont les paramètres inconnus à résoudre. [0069] Dans un autre mode de réalisation, une deuxième équation incorpore les paramètres résolus préalablement inconnus (an, wn, hi) dans un modèle qui définit un ensemble N de tangentes hyperboliques. Par exemple, dans un mode de réalisation, la deuxième équation correspond à 20 l'équation suivante: Équation 2: f(P) a Wn, tn) = aN EnN,i an+1 .tanh(w P n)) avec les contraintes suivantes 25 > 0, Vn c [1,N] n, N E N an+i < an,\In E [1, N - 1] n,N E N 30 où P est la transformée logarithmique de la pression capillaire normalisée et N est le nombre de tangentes hyperboliques établies pour le modèle. 3036820 16 [0070] Dans un mode de réalisation, le nombre de tangentes hyperboliques du modèle dans l'équation 2 est prédéterminé. Par exemple, la Figure 6 montre des données de pression capillaire à partir d'un système de gorge à 3 pores, par conséquent, les équations 1 et 2 seraient établies à N=3. [0071] Dans un mode de réalisation, les facteurs d'échelle (an+i-an) de chaque tangente 5 hyperbolique dans l'ensemble N sont liés ensemble de telle sorte que la somme des tangentes hyperboliques soit limitée entre 2a1 et 2aN. La liaison peut forcer la séparation des tangentes hyperboliques entre différentes gorges à pores. Par exemple, forcer une tangente hyperbolique par gorge à un pore au lieu d'une tangente hyperbolique sur une gorge à 3 pores et deux autres tangentes hyperboliques sans aucune contribution. En d'autres termes, comme cela est illustré dans 10 la Figure 8, chaque tangente hyperbolique peut être limitée à une gorge à un pore. [0072] Dans un mode de réalisation, les contraintes présentes dans l'équation 2 sont configurées de manière à limiter les tangentes hyperboliques à des courbes de pression capillaire réalistes et à améliorer la stabilité du modèle. Par exemple, les tangentes hyperboliques peuvent être triées par le nombre de gorges à pores dans le système, avec la "première" tangente hyperbolique qui 15 commence sur la gauche. Chaque gorge à pores et la tangente hyperbolique combinée correspondante peuvent être établies comme des fonctions décroissantes monotones. Par exemple, les Figures 3, 4, et 5 illustrent un modèle de tangentes hyperboliques dans un système de pression capillaire et de saturation en eau selon un mode de réalisation. La Figure 3 illustre une seule tangente hyperbolique 310 dans un système de pression capillaire et de saturation en eau créée en 20 utilisant l'équation 2 ci-dessus avec les contraintes dans celle-ci. L'axe x représente la pression capillaire et l'axe y représente la saturation en eau. La Figure 4 illustre deux tangentes hyperboliques 320 et 330 créées en utilisant l'équation 2 ci-dessus avec les contraintes dans celle-ci. Comme cela est illustré dans la Figure 4, une troisième tangente hyperbolique 340 est la somme des tangentes hyperboliques 320 et 330 et représente un système de gorge à deux pores. 25 [0073] La Figure 5 montre deux tangentes hyperboliques 350 et 360 créées sans les contraintes dans l'équation 2 ci-dessus, et une troisième tangente hyperbolique 370 qui est la somme des tangentes hyperboliques 350 et 360. Comme cela est illustré dans la Figure 5, la troisième tangente hyperbolique 370 ne peut pas représenter une courbe de pression capillaire réaliste parce que les tangentes hyperboliques non contraintes sous-jacentes 350 et 360 se dirigent dans des directions 30 différentes. Une tangente hyperbolique peut également ne pas représenter une courbe de pression capillaire réaliste si elle résulte en une fonction décroissante non monotone. 3036820 17 [0074] Dans un mode de réalisation, une routine d'optimisation non linéaire est utilisée pour trouver les paramètres de meilleur ajustement. Par exemple, une routine d'optimisation non linéaire configurée pour gérer des contraintes d'inégalités linéaires, telle qu'une programmation quadratique séquentielle, peut être utilisée pour trouver les paramètres de meilleur ajustement. 5 [0075] Les Figures 6, 7 et 8 montrent un modèle de pression capillaire selon des modes de réalisation de la présente divulgation. La Figure 6 illustre des données de pression capillaire à partir d'un système de gorge à pores multiples. La Figure 7 illustre une courbe de meilleur ajustement 410 sur les données de pression capillaire. Comme cela est illustré dans la Figure 7, la courbe de meilleur ajustement 410 est la somme de trois tangentes hyperboliques 420, 430 et 440. La Figure 10 8 illustre les trois tangentes hyperboliques 420, 430, et 440 décalées qui montrent quelle tangente hyperbolique correspond à quelle gorge à pores. [0076] Comme cela est illustré dans les Figures 6 à 8, un modèle de pression capillaire qui incorpore les équations 1 et 2 montre un bon ajustement aux données de pression capillaire mesurées de puits, et un nombre de tangentes hyperboliques N peut être établi pour approcher le 15 nombre de gorges à pores dans le système. Dans certains modes de réalisation, un bon ajustement est déterminé par la quantité d'erreur dans l'équation 1: l'erreur moindre sur l'équation 1 signifiant le meilleur ajustement, tandis qu'une valeur d'erreur plus élevée indique une qualité moindre de l'ajustement. [0077] Dans un mode de réalisation, une fonction de hauteur de saturation est créée en combinant 20 le modèle de pression capillaire des équations 1 et 2 de concert avec des équations qui incorporent d'autres propriétés physiques du réservoir. Par exemple, un modèle de pression capillaire peut être créé en utilisant les équations 1 et 2 pour approcher des données de pression capillaire mesurées tout en incorporant simultanément des données de perméabilité afm de créer une fonction de hauteur de saturation. 25 [0078] Comme cela est décrit ci-dessus, le modèle de pression capillaire des équations 1 et 2 peut être soumis à une correction supplémentaire. Des données de pression capillaire obtenues à partir d'échantillons de bouchon de carotte provenant du réservoir peuvent être soumises à des effets de conformité ou de fermeture selon la manière avec laquelle le bouchon de carotte est obtenu. Lorsqu'un bouchon de carotte est immergé dans du mercure pour une procédure HPMI (injection 30 de mercure à haute pression), avec des pressions aussi élevées, par exemple, qu'environ 60000 psi, ou une procédure MICP (pression capillaire d'injection de mercure), avec des pressions inférieures, par exemple, à environ 2000 psi, l'injection de mercure commence d'abord à une pression très 3036820 18 basse. Toutefois, l'interface du bouchon de carotte peut présenter une certaine rugosité de bord etiou d'autres altérations qui ne sont pas représentatives du système de pores interne du réservoir. Par exemple, l'interface du bouchon de carotte peut être plus ou moins friable, ou peut présenter une rugosité altérée par le procédé de coupe du bouchon de carotte ainsi que par la manutention, le 5 vieillissement et la manipulation avant les expériences d'injection de mercure, avec comme résultat que les pores exposés sur l'interface de l'échantillon de bouchon de carotte sont plus grands qu'ils ne l'étaient à l'origine. Lorsque la pression de mercure augmente, le mercure commencera à se répandre autour du bouchon de carotte, et pénètrera dans l'interface rugueuse endommagée à une certaine pression dépendant de l'étendue de l'endommagement du bouchon de carotte. Par 10 exemple, plus il est endommagé, plus les ouvertures créées sont grandes, et plus la pression à laquelle le mercure commencera à pénétrer dans le bouchon sera basse. Ces pores endommagés ne peuvent pas être représentatifs du système de pores du réservoir, et par conséquent cela doit être corrigé. [0079] Tel qu'il est employé ici, le terme "pression d'entrée" fait référence au point sur la courbe 15 de pression capillaire auquel le mercure entre initialement dans l'échantillon lors d'une expérience d'injection de mercure. La pression d'entrée est également quelquefois appelée pression de déplacement, et peut être une caractéristique de la roche qui est liée à la force requise dans le réservoir pour que le pétrole se trouve au-dessus du niveau d'eau libre. [0080] Dans un mode de réalisation, un procédé peut être utilisé pour exécuter automatiquement 20 une correction de fermeture sur des courbes de pression capillaire en utilisant de multiples tangentes hyperboliques liées. La correction peut être réalisée par extrapolation pour des échantillons de bouchon de carotte homogènes, tels que du grès, ou par normalisation pour des échantillons de bouchon de carotte hétérogènes, tels que du carbonate de roche à structure cellulaire sans fractures. En outre, le procédé 600 peut également être utilisé pour déterminer la 25 pression d'entrée qui correspond à la courbe de pression capillaire. [0081] La Figure 9 est un organigramme d'un procédé pour corriger une courbe de pression capillaire et obtenir la pression d'entrée pour la courbe capillaire. Comme cela est illustré dans la Figure 9, un procédé 600 peut commencer avec la création d'une courbe de pression capillaire lors de l'opération 610. Le procédé 600 peut ensuite continuer avec la détermination d'une pression de 30 coupure de correction de fermeture lors de l'opération 620 et la correction de la courbe de pression capillaire lors de l'opération 630. Dans un mode de réalisation, la courbe de pression capillaire est corrigée par l'intermédiaire d'une normalisation lors de l'opération 640 et la pression d'entrée est 3036820 19 déterminée lors de l'opération 650. Dans un autre mode de réalisation, la courbe de pression capillaire est corrigée par l'intermédiaire d'une extrapolation lors de l'opération 660 et la pression d'entrée est déterminée lors de l'opération 670. [0082] La Figure 10 montre une courbe de pression capillaire avec 2 gorges à pores. Comme cela 5 est illustré dans la Figure 10, lors de l'opération 610, une courbe de pression capillaire 700 est créée en liant de multiples tangentes hyperboliques 710 et 720 en utilisant l'équation 2 pour approcher les données de pression capillaire obtenues 750, de façon similaire à ce qui a été décrit ci-dessus en relation avec les Figures 2 à 8. [0083] A l'opération 620, la pression de coupure de correction de fermeture est déterminée. Dans 10 un mode de réalisation, la pression de coupure de correction de fermeture correspond au minimum local de la dérivée seconde de la courbe de pression capillaire où la pression capillaire est la plus basse. Par exemple, la courbe de pression capillaire 700 peut être créée en utilisant l'équation 2: Équation 2: f (P, an, w, t) = ai + a N zEri=i (an+i - an).tanh(wn.(P ta)) 15 dans laquelle N = 2 pour représenter un système de gorge à deux pores. [0084] La dérivée seconde de la courbe de pression capillaire 700 peut alors être calculée en utilisant l'équation suivante: Équation 5: (17 2n-1- 1 e2z (n\ tanh.z ( 1)k. 2 dzn (1+e2.)1+1.kz (n 20 dans laquelle k) est le nombre eulérien et n = 2. [0085] La Figure 11 montre une dérivée seconde 800 de la courbe de pression capillaire 700. Comme cela est illustré dans la Figure 11, la dérivée seconde présente deux minima locaux 801 et 802, avec la pression la plus basse au minimum local 802. Par conséquent, la pression de coupure de correction de fermeture 803 correspond à la pression au minimum local 802 le long de la ligne 25 "A". [0086] La Figure 12 montre une courbe de pression capillaire rognée 770. Comme cela est illustré dans la Figure 12, toutes les données de pression capillaire inférieures à la pression de coupure de correction de fermeture 803 sont supprimées. Comme cela est décrit ci-dessus, les données de 3036820 20 pression capillaire supprimées peuvent correspondre à des données capillaires soumises à des effets de conformité ou de fermeture. [0087] Une fois que la pression de coupure de correction de fermeture 803 a été déterminée, la courbe de pression capillaire 770 peut être corrigée et la pression d'entrée peut être déterminée lors 5 de l'opération 630. [0088] La correction peut être réalisée par normalisation lors de l'opération 640 ou par extrapolation lors de l'opération 660. [0089] Dans les opérations 640 et 650, la correction est réalisée par normalisation et la pression de coupure de correction de fermeture correspond à la pression d'entrée. 10 [0090] La Figure 13 montre une courbe de pression capillaire normalisée. Comme cela est illustré dans la Figure 13, une courbe de pression capillaire normalisée 900 est créée en normalisant la courbe de pression capillaire rognée 770 de telle sorte que la pression de coupure de correction de fermeture 803 corresponde à une saturation de phase de mouillage de 100 % lors de l'opération 640. La pression de coupure de correction de fermeture 803 correspondra à la pression d'entrée 850 15 le long de la ligne "A" lors de l'opération 650. [0091] La normalisation de la courbe de pression capillaire rognée 770 peut être effectuée par l'intermédiaire de procédés de normalisation standard dans l'industrie pétrolière. Par exemple, la normalisation peut être effectuée en prenant la valeur de saturation qui correspond à la pression de coupure de correction de fermeture 803 et en divisant la totalité de la saturation de la pression 20 capillaire par cette valeur, de telle sorte que le dernier point de la courbe de pression capillaire normalisée correspondra à la saturation à la pression de coupure de correction de fermeture 803 et/ou à la saturation à la pression de coupure de correction de fermeture 803 = 1. [0092] A l'opération 660, la correction est réalisée par extrapolation en utilisant la pression de coupure de correction de fermeture et une première pression de coupure dérivée, et la pression 25 d'entrée est la pression à la saturation de phase de mouillage de 100 %. [0093] Dans un mode de réalisation, la première pression de coupure dérivée correspond au minimum local de la dérivée première de la courbe de pression capillaire 700 où la pression capillaire est la plus basse. Par exemple, la dérivée première de la courbe de pression capillaire 700 peut être calculée en utilisant l'équation suivante: 30 Équation 6: N _ dP Ln=1(an+1 an).wn..( tanh(wn.(P a)).tarth(w P a)) 3036820 21 [0094] La Figure 14 montre une dérivée première 880 et la dérivée seconde 800 d'une courbe de pression capillaire rognée 770. Comme cela est illustré dans la Figure 14, la dérivée première comprend deux minima locaux 881 et 882, avec la pression la plus basse au minimum local 882. Par conséquent, la première pression de coupure dérivée 883 correspond à la pression au minimum 5 local 882 le long de la ligne "B." Comme cela est illustré dans la Figure 14, la pression de coupure de correction de fermeture 803 correspond au minimum local 802 de la dérivée seconde 800, et la courbe de pression capillaire 770 est rognée afin de supprimer toutes les données de pression capillaire inférieures à la pression de coupure de correction de fermeture 803. [0095] La Figure 15 montre une courbe de pression capillaire extrapolée. Comme cela est illustré 10 dans la Figure 15, une courbe de pression capillaire extrapolée 980 est créée par une analyse de régression entre la première pression de coupure dérivée 883 et la pression de coupure de correction de fermeture 803. Par exemple, comme cela est illustré dans la Figure 15, une régression des moindres carrés linéaire entre la première pression de coupure dérivée 883 et la pression de coupure de correction de fermeture 803 peut être calculée lors de l'opération 660. 15 [0096] La régression des moindres carrés linéaire peut être utilisée pour prédire la saturation en dessous de la pression de coupure de correction de fermeture 803 jusqu'à l'obtention d'une saturation de phase de mouillage de 100 %. Par exemple, la régression des moindres carrés linéaire peut être utilisée pour extrapoler la courbe de pression capillaire extrapolée 980 au point de saturation de phase de mouillage de 100 % en utilisant une équation linéaire, voir "D". 20 [0097] A l'opération 670, comme cela est illustré dans la Figure 15, la pression d'entrée 850 correspond à une pression à une saturation de phase de mouillage de 100 % le long de la ligne "C." [0098] Par conséquent, dans certains modes de réalisation, le procédé 600 peut exécuter automatiquement la correction de fermeture et déterminer la pression d'entrée de données de pression capillaire pour des échantillons de bouchon de carotte homogènes et hétérogènes en 25 utilisant une analyse de normalisation ou d'extrapolation. [0099] Dans certains modes de réalisation, les procédés de la présente divulgation peuvent être exécutés par un système de calcul. La Figure 16 illustre un exemple d'un tel système de calcul 500, selon certains modes de réalisation. Le système de calcul 500 peut comprendre un ordinateur ou un système d'ordinateur 501A, qui peut être un système d'ordinateur individuel 501A ou un 30 agencement de systèmes d'ordinateur distribués. Le système d'ordinateur 501A comprend un ou plusieurs module(s) d'analyse 502 qui est (sont) configuré(s) pour exécuter différents tâches selon certains modes de réalisation, telles qu'un ou plusieurs procédé(s) divulgué(s) ici. Pour exécuter 3036820 22 ces différentes tâches, le module d'analyse 502 opère de façon indépendante, ou en coordination avec un ou plusieurs processeur(s) 504, qui est (sont) connecté(s) à un ou plusieurs support(s) de stockage 506. Le(s) processeur(s) 504 est (sont) également connecté(s) à une interface réseau 507 afin de permettre au système d'ordinateur 501A de communiquer sur un réseau de données 509 5 avec un ou plusieurs système(s) d'ordinateur et/ou système(s) de calcul supplémentaire(s), tels que 501B, 501C et/ou 501D (notons que les systèmes d'ordinateur 501B, 501C et/ou 501D peuvent partager ou non la même architecture que le système d'ordinateur 501A, et peuvent être situés à des emplacements physiques différents, par exemple, les systèmes d'ordinateur 501A et 501B peuvent être situés dans une installation de traitement, tout en étant en communication avec un ou 10 plusieurs système(s) d'ordinateur tels que 501C et/ou 501D qui est (sont) situé(s) dans un ou plusieurs centre(s) de données, et/ou situé(s) dans différents pays sur des continents différents). [00100] Un processeur peut comprendre un microprocesseur, un micro-dispositif de commande, un module ou un sous-système de processeur, un circuit intégré programmable, un réseau de portes programmable ou un autre dispositif de commande ou de calcul. 15 [00101] Les supports de stockage 506 peuvent être mis en oeuvre sous la forme d'un ou de plusieurs support(s) de stockage lisible(s) par un ordinateur ou lisible(s) par une machine. Notons que bien que dans l'exemple de mode de réalisation qui est montré dans la Figure 16 les supports de stockage 506 soient représentés comme étant situés à l'intérieur du système d'ordinateur 501A, dans certains modes de réalisation, les supports de stockage 506 peuvent être distribués à l'intérieur 20 et/ou en travers de multiples coffrets internes ou externes du système de calcul 501A et/ou de systèmes de calcul supplémentaires. Les supports de stockage 506 peuvent comprendre une ou plusieurs forme(s) différente(s) de mémoire comprenant des dispositifs de mémoire à semiconducteur tels que des mémoires vives dynamiques ou statiques (DRAM ou SRAM), des mémoires mortes effaçables et programmables (EPROM), des mémoires mortes électriquement 25 effaçables et programmables (EEPRO) et des mémoires flash, des disques magnétiques tels que des disques fixes, des disquettes et des disques amovibles, d'autres supports magnétiques comprenant des bandes, des supports optiques tels que des disques compacts (CD) ou des disques vidéo numériques (DVD), des disques BLUERAY® ou d'autres types de stockage optique, ou d'autres types de dispositifs de stockage. Notons que les instructions discutées ci-dessus peuvent être 30 fournies sur un support de stockage lisible par un ordinateur ou lisible par une machine, ou peuvent être fournies sur de multiples supports de stockage lisibles par un ordinateur ou lisibles par une machine distribués dans un grand système comprenant éventuellement plusieurs noeuds. Ce(s) 3036820 23 support(s) lisible(s) par un ordinateur ou lisible(s) par une machine est (sont) considéré(s) comme faisant partie d'un article (ou d'un article de manufacture). Un article ou un article de manufacture peut faire référence à n'importe quel(s) composant unique ou composants multiples. Le(s) support(s) de stockage peu(ven)t être situé(s) soit dans la machine qui exécute les instructions 5 lisibles par une machine, soit situé(s) à un site distant à partir duquel les instructions lisibles par une machine peuvent être téléchargées par l'intermédiaire d'un réseau pour être exécutées. [00102] Dans certains modes de réalisation, le système de calcul 500 contient un ou plusieurs module(s) de modélisation 508. Dans l'exemple de système de calcul 500, le système d'ordinateur 501A comprend le module de modélisation 508. Dans certains modes de réalisation, 10 un seul module de modélisation peut être utilisé pour exécuter au moins certains aspects d'un ou de plusieurs modes de réalisation des procédés divulgués ici. Dans des modes de réalisation alternatifs, une pluralité de modules de modélisation peut être utilisée pour exécuter au moins certains aspects des procédés décrits ici. [00103] On appréciera que le système de calcul 500 est un exemple d'un système de calcul, 15 et que le système de calcul 500 peut comprendre plus ou moins de composants que ce qui est montré, peut combiner des composants supplémentaires qui ne sont pas décrits dans l'exemple de mode de réalisation de la Figure 16, et/ou le système de calcul 500 peut présenter une configuration ou un agencement différent(e) des composants représentés dans la Figure 16. Les différentes composants montrés dans la Figure 16 peuvent être mises en oeuvre dans un matériel, un logiciel ou 20 une combinaison à la fois de matériel et de logiciel, comprenant un ou plusieurs circuit(s) intégré(s) spécifique(s) de traitement et/ou d'application de signal. [00104] En outre, des aspects des procédés de traitement décrits ici peuvent être mis en oeuvre en exécutant un ou plusieurs module(s) fonctionnel(s) dans un appareil de traitement d'informations tel que des processeurs universels ou des puces spécifiques d'application, telles que 25 des ASIC, des FPGA, des PLD, ou d'autres dispositifs appropriés. Ces modules, combinaisons de ces modules et/ou leur combinaison avec un matériel universel sont inclus à l'intérieur de la portée de la protection de l'invention. [00105] Des interprétations géologiques, des modèles géologiques et/ou d'autres aides d'interprétation peuvent être affinés d'une façon itérative; ce concept est applicable aux procédés 30 discutés ici. Ceci peut comprendre l'utilisation de boucles de rétroaction exécutées sur une base algorithmique, tel qu'un dispositif de calcul (par exemple, un système de calcul 500, Figure 16), et/ou par l'entremise d'une commande manuelle par un utilisateur qui peut produire des 3036820 24 déterminations permettant de savoir si une étape, une action, un gabarit, un modèle ou un ensemble de courbes donné(e) est devenu suffisamment précis(e) pour évaluer une formation géologique tridimensionnelle souterraine considérée.

1001061 La présente divulgation a été décrite en se référant aux modes de réalisation. Bien 5 que quelques modes de réalisation aient été présentés et décrits, l'homme du métier appréciera que des changements peuvent être apportés à ces modes de réalisation sans trahir les principes et l'esprit de la description détaillée qui précède. Il est prévu que la présente divulgation soit interprétée comme englobant ces modifications et variations dans la mesure où elles rentrent à l'intérieur de la portée des revendications annexées ou des équivalents de celles-ci. 10

Claims (14)

1. REVENDICATIONS1. Procédé pour corriger une courbe de pression capillaire, comprenant les étapes suivantes: obtention de données de pression capillaire représentant une pression capillaire dans un réservoir, création (610) d'une courbe de pression capillaire en utilisant de multiples tangentes hyperboliques liées pour approcher les données de pression capillaire obtenues; détermination (620) d'une pression de coupure de correction de fermeture de la courbe de pression capillaire; et correction (630) de la courbe de pression capillaire, dans lequel la correction de la courbe de pression capillaire comprend au moins une étape parmi une normalisation (640) de la courbe de pression capillaire et une extrapolation (660) de la courbe de pression capillaire.
2. 2. Procédé selon la revendication 1, comprenant en outre les étapes suivantes: détermination d'une pression d'entrée pour la courbe de pression capillaire après la correction de la courbe de pression capillaire, dans lequel la détermination de la pression de coupure de correction de fermeture de la courbe de pression capillaire comprend une détermination d'un minimum local de la dérivée seconde de la courbe de pression capillaire, et dans lequel la pression de coupure de correction de fermeture correspond au minimum local où la pression capillaire est la plus basse.
3. 3. Procédé selon la revendication 2, dans lequel la normalisation de la courbe de pression capillaire comprend un réglage de la pression de coupure de correction de fermeture afin de correspondre à une saturation de phase de mouillage de 100 %, et dans lequel la pression d'entrée correspond à la pression de coupure de correction de fermeture après le réglage de la pression de coupure de correction de fermeture.
4. 4. Procédé selon la revendication 2, dans lequel l'extrapolation de la courbe de pression capillaire comprend les étapes suivantes: détermination d'une première pression de coupure dérivée; 3036820 26 exécution d'une analyse de régression entre la première pression de coupure dérivée et la pression de coupure de correction de fermeture; et extrapolation de la courbe de pression capillaire au point de saturation de phase de mouillage de 100%, 5 dans lequel la première pression de coupure dérivée correspond au premier minimum local dérivé où la pression capillaire est la plus basse, et dans lequel la pression d'entrée correspond à une pression au point de saturation de phase de mouillage de 100 %. 10
5. 5. Procédé selon la revendication 1, dans lequel la courbe de pression capillaire est définie par l'équation suivante: f(P,an,vin,tn) -= ai + aN + E_i (an+i an).tanh(wn,.(P te)) 15 avec les contraintes suivantes: ton. > 0, E [1, N] n,N E N an+1< an,Vn E [1,A 1] nNeN 20 où P représente une transformée logarithmique d'une pression capillaire normalisée et N représente le nombre de tangentes hyperboliques.
6. 6. Support lisible par ordinateur non transitoire (506) pour stocker des instructions qui, 25 lorsqu'elles sont exécutées par un ou plusieurs processeur(s) (504) d'un système de calcul (501A), amènent le système de calcul à exécuter des opérations, les opérations comprenant: obtention de données de pression capillaire représentant une pression capillaire dans un réservoir, création (610) d'une courbe de pression capillaire en utilisant de multiples tangentes 30 hyperboliques liées pour approcher les données de pression capillaire obtenues; détermination (620) d'une pression de coupure de correction de fermeture de la courbe de pression capillaire; correction (630) de la courbe de pression capillaire; et 3036820 29469/302FR I 27 détermination (650, 670) d'une pression d'entrée pour la courbe de pression capillaire après la correction de la courbe de pression capillaire, dans lequel la correction de la courbe de pression capillaire comprend au moins une opération parmi une normalisation (640) de la courbe de pression capillaire et une extrapolation 5 (660) de la courbe de pression capillaire.
7. 7. Support lisible par ordinateur non transitoire selon la revendication 6, dans lequel la détermination de la pression de coupure de correction de fermeture de la courbe de pression capillaire comprend une détermination d'un minimum local de la dérivée seconde de la courbe de 10 pression capillaire, dans lequel la pression de coupure de correction de fermeture correspond au minimum local où la pression capillaire est la plus basse, et dans lequel la correction de la courbe de pression capillaire comprend un rognage de la courbe de pression capillaire afin d'éliminer les données de pression capillaire en dessous de la pression de coupure de correction de fermeture. 15
8. 8. Support lisible par ordinateur non transitoire selon la revendication 7, dans lequel la normalisation de la courbe de pression capillaire comprend un réglage de la pression de coupure de correction de fermeture afin de correspondre à une saturation de phase de mouillage de 100 %, et dans lequel la pression d'entrée correspond à la pression de coupure de correction de fermeture après le réglage de la pression de coupure de correction de fermeture, et 20 dans lequel l'extrapolation de la courbe de pression capillaire comprend une détermination d'une première pression de coupure dérivée, une exécution d'une analyse de régression entre la première pression de coupure dérivée et la pression de coupure de correction de fermeture, et une extrapolation de la courbe de pression capillaire au point de saturation de phase de mouillage de 100 %, dans lequel la première pression de coupure dérivée correspond au premier minimum local 25 dérivé où la pression capillaire est la plus basse, et dans lequel la pression d'entrée correspond à une pression au point de saturation de phase de mouillage de 100 %.
9. 9. Support lisible par ordinateur non transitoire selon la revendication 6, dans lequel la courbe de pression capillaire est définie par l'équation suivante: 30 f(I), wn, tn) = a + (an+i - a ).ta h(wie - tn)) 3036820 28 avec les contraintes suivantes: wn, > 0, Vn E [1, N] rt, N E N 5 an,+1 < an, Vn E [11N -1] ral N E N où P représente une transformée logarithmique d'une pression capillaire normalisée et N représente le nombre de tangentes hyperboliques. 10
10. 10. Système de calcul (501A), comprenant: un ou plusieurs processeur(s) (504); et un système de mémoire comprenant un ou plusieurs support(s) lisible(s) par ordinateur non transitoire(s) (506) pour stocker des instructions qui, lorsqu'elles sont exécutées par ledit/lesdits un ou plusieurs processeur(s) du système de calcul, amènent le système de calcul à exécuter des 15 opérations, les opérations comprenant: obtention de données de pression capillaire représentant une pression capillaire dans un réservoir, création (610) d'une courbe de pression capillaire en utilisant de multiples tangentes hyperboliques liées pour approcher les données de pression capillaire obtenues; 20 détermination (620) d'une pression de coupure de correction de fermeture de la courbe de pression capillaire; correction (630) de la courbe de pression capillaire; et détermination (650, 670) d'une pression d'entrée pour la courbe de pression capillaire après la correction de la courbe de pression capillaire, 25 dans lequel la correction de la courbe de pression capillaire comprend au moins une opération parmi une normalisation (640) de la courbe de pression capillaire et une extrapolation (660) de la courbe de pression capillaire.
11. 11. Système de calcul selon la revendication 10, dans lequel la détermination de la pression de 30 coupure de correction de fermeture de la courbe de pression capillaire comprend une détermination d'un minimum local de la dérivée seconde de la courbe de pression capillaire, dans lequel la pression de coupure de correction de fermeture correspond au minimum local où la pression capillaire est la plus basse, et dans lequel la correction de la courbe de pression capillaire comprend 3036820 29 un rognage de la courbe de pression capillaire afin d'éliminer les données de pression capillaire en dessous de la pression de coupure de correction de fermeture.
12. 12. Système de calcul selon la revendication 11: 5 dans lequel la normalisation de la courbe de pression capillaire comprend un réglage de la pression de coupure de correction de fermeture afin de correspondre à une saturation de phase de mouillage de 100 %, et dans lequel la pression d'entrée correspond à la pression de coupure de correction de fermeture après le réglage de la pression de coupure de correction de fermeture. 10
13. 13. Système de calcul selon la revendication 11, dans lequel l'extrapolation de la courbe de pression capillaire comprend: détermination d'une première pression de coupure dérivée; exécution d'une analyse de régression entre la première pression de coupure dérivée et la pression de coupure de correction de fermeture; et 15 extrapolation de la courbe de pression capillaire au point de saturation de phase de mouillage de 100 %, dans lequel la première pression de coupure dérivée correspond au premier minimum local dérivé où la pression capillaire est la plus basse, et dans lequel la pression d'entrée correspond à une pression au point de saturation de phase 20 de mouillage de 100 %.
14. 14. Système de calcul selon la revendication 10, dans lequel la courbe de pression capillaire est définie par l'équation suivante: 25 f(P)art,Ivn)tn vq't1 =- + aN an).ta, h(tiin,.(P e)) 30 avec les contraintes suivantes: > 0, \in E [1, IV] n,IV E N an+i < a,,Vra E [1, N - 1] n, N E N où P représente une transformée logarithmique d'une pression capillaire normalisée et N représente le nombre de tangentes hyperboliques.