FR2811430A1 - Methode de modelisation permettant de predire en fonction du temps la composition detaillee de fluides porudits par un gisement souterrain en cours de production - Google Patents
Methode de modelisation permettant de predire en fonction du temps la composition detaillee de fluides porudits par un gisement souterrain en cours de production Download PDFInfo
- Publication number
- FR2811430A1 FR2811430A1 FR0009008A FR0009008A FR2811430A1 FR 2811430 A1 FR2811430 A1 FR 2811430A1 FR 0009008 A FR0009008 A FR 0009008A FR 0009008 A FR0009008 A FR 0009008A FR 2811430 A1 FR2811430 A1 FR 2811430A1
- Authority
- FR
- France
- Prior art keywords
- thermodynamic
- fluids
- phases
- input data
- phase
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 45
- 239000000203 mixture Substances 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 11
- 239000012071 phase Substances 0.000 claims abstract description 66
- 238000004088 simulation Methods 0.000 claims abstract description 31
- 229930195733 hydrocarbon Natural products 0.000 claims abstract description 26
- 150000002430 hydrocarbons Chemical class 0.000 claims abstract description 26
- 239000004215 Carbon black (E152) Substances 0.000 claims abstract description 24
- 239000007788 liquid Substances 0.000 claims abstract description 12
- 239000008346 aqueous phase Substances 0.000 claims abstract description 5
- 238000002347 injection Methods 0.000 claims abstract description 5
- 239000007924 injection Substances 0.000 claims abstract description 5
- 238000006243 chemical reaction Methods 0.000 claims abstract description 3
- 238000006073 displacement reaction Methods 0.000 claims abstract description 3
- 230000009172 bursting Effects 0.000 claims description 9
- 239000007791 liquid phase Substances 0.000 abstract description 4
- 238000004880 explosion Methods 0.000 abstract 5
- 239000007792 gaseous phase Substances 0.000 abstract 1
- 239000007789 gas Substances 0.000 description 23
- 239000000470 constituent Substances 0.000 description 12
- 230000006870 function Effects 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 4
- 238000004090 dissolution Methods 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000004514 thermodynamic simulation Methods 0.000 description 3
- 230000008016 vaporization Effects 0.000 description 3
- 238000009834 vaporization Methods 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 101100365516 Mus musculus Psat1 gene Proteins 0.000 description 2
- RAHZWNYVWXNFOC-UHFFFAOYSA-N Sulphur dioxide Chemical compound O=S=O RAHZWNYVWXNFOC-UHFFFAOYSA-N 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 238000004513 sizing Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 239000012808 vapor phase Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/26—Oils; Viscous liquids; Paints; Inks
- G01N33/28—Oils, i.e. hydrocarbon liquids
- G01N33/2823—Raw oil, drilling fluid or polyphasic mixtures
Landscapes
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medicinal Chemistry (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- General Chemical & Material Sciences (AREA)
- Food Science & Technology (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
- Méthode utilisant une modélisation de type " Black Oil " pour prédire en fonction du temps la composition détaillée de fluides produits par un gisement souterrain en cours de production, combinée avec une étape d'éclatement ou " delumping " permettant une représentation thermodynamique détaillée des fluides de gisement. Comme données d'entrée pour le modèle, on introduit non seulement les paramètres thermodynamiques des fluides tels que viscosité, facteur de volume, densité, rapport de dissolution gaz-huile, etc. (sous forme de tables, et/ ou par corrélation, comme fonction de la pression, de la température si celle-ci varie) et, le cas échéant un paramètre additionnel conservant une mémoire de la composition du gaz tel par exemple que la densité du gaz), mais également des données concernant les variations en fonction des mêmes " abscisses ", des paramètres de phase nécessaires au " delumping ", sans que ceux-ci soient utilisés durant la simulation " Black Oil " des écoulements.- Application par exemple à l'établissement de profils prévisionnels de la composition détaillée des hydrocarbures produits par un gisement.
Description
La présente invention concerne une méthode de modélisation permettant de
prédire en fonction du temps la composition détaillée de fluides produits par un gisement souterrain en cours de production, et notamment d'un gisement d'hydrocarbures en cours de production. Un tel modèle permet à l'ingénieur de gisement d'établir des profils compositionnels détaillés qui, sont nécessaires en particulier au dimensionnement prévisionnel et à la gestion d'installations de surface, telles que séparateurs, unités de traitement, lignes de transport, etc., donc utiles en ingénierie de surface et de
procédés.
Etat de la technique La modélisation des écoulements dans un réservoir pétrolier est fondée essentiellement sur l'application au réservoir préalablement maillé (ou à une portion de celui-ci) de la loi de Darcy bien connue décrivant l'écoulement de fluides en milieux poreux, de la loi de conservation des masses dans chaque unité de volume, des relations thermodynamiques régissant l'évolution des paramètres de phase des fluides tels que la viscosité, la masse volumique, sur des conditions initiales, sur des
conditions aux limites de fermeture de la structure, et sur des conditions aux puits.
Le modèle dit " Black Oil " désigné ci-après en abrégé par B.O., est l'un des plus utilisés en simulation pétrolière. Il permet de décrire un écoulement tridimensionnel et triphasique (eau-huile-gaz) compressible. Les effluents pétroliers intervenant dans ce modèle sont décrits généralement par un constituant eau, et deux constituants pour le fluide de gisement, le mot de constituant recouvrant ici la notion de composant (comme H20 pour l'eau) et la notion de pseudo-composant (regroupement de composants). Les constituants intervenant dans ce modèle sont au nombre de trois: un constituant eau (E), un constituant hydrocarbure lourd (L), et un constituant hydrocarbure léger (V). Dans un modèle de type B.O. dit " strict ", le constituant (E) n'est présent que dans la phase eau, le constituant (L) n'est présent que dans la phase liquide hydrocarbure (appelée huile ou condensat), le constituant
(V) se partage entre les phases hydrocarbures liquide et vapeur (phase appelée gaz).
Même si l'usage des modèles B.O. est à déconseiller dans certains cas de gaz à condensats, il n'en reste pas moins applicable à un grand nombre de cas industriels. Un autre modèle de simulation bien connu dit " compositionnel " est aussi utilisé, o les fluides hydrocarbures sont représentés par un nombre plus élevé de composants, au moins trois, souvent davantage. La modélisation des flux de ces fluides mieux détaillés se traduit par des temps de calcul très importants (beaucoup plus longs que ceux requis pour une modélisation de type B.O.) du fait du plus grand nombre de constituants mais aussi parce qu'il est souvent nécessaire de réduire la taille des mailles pour limiter l'erreur numérique et par conséquent d'augmenter le
nombre de mailles.
Pour des raisons pratiques, les fluides en place sont décrits comme constitués d'un nombre de composants ou pseudo-composants beaucoup plus réduit que le nombre réel de composants, de façon que les calculs de modélisation puissent être
menés dans un délai raisonnable. Une composition réduite à quelque 5 ou 10 pseudo-
composants s'avère en général suffisante pour représenter le comportement des
fluides dans le gisement.
Par la demande de brevet WO 99/42937 et aussi le papier de C. Leibovici et J. Barker "A Method for Delumping the Results of a Compositional Reservoir Simulation" SPE 49068, présenté à la SPE Annual Technical Conference and Exhibition New Orleans, 27-30 Septembre 1998, on connaît une méthode de prévision de l'évolution de la composition détaillée au cours du temps à partir des calculs effectués dans une simulation de type compositionnel de fluides décrits par un certain nombre réduit de pseudo- composants (principe de représentation 'groupée' ou 'regroupée' obtenue par une opération dite de 'lumping'), le nombre de composants étant au moins de trois. La méthode permet donc de prédire quels auraient été les résultats d'une simulation de réservoir si l'on avait utilisé une modélisation finement détaillée o les fluides sont représentés par un nombre plus élevé de composants. C'est une opération connue des gens de l'art sous l'appellation
de "delumping", que l'on désignera par la suite aussi par "éclatement".
Le principe de l'étape d'éclatement antérieure est de calculer le coefficient ADo et les n coefficients ADp (soit n+l coefficients, n étant le nombre de paramètres de l'équation d'état) d'une équation générale connue, publiée antérieurement dans un papier de C.F. Leibovici, E.H. Stenby, K. Knudsen, "A Consistent Procedure for Pseudo-Component Delumping", Fluid Phase Equilibria, 1996, 117, 225-232.: n (1) Ln(ki) = ADo + _ADp rlpi p=l o les IIpi sont des paramètres fixes de caractérisation du constituant i dans l'équation d'état pour une représentation thermodynamique donnée, à partir des constantes d'équilibre ki de chaque constituant de la représentation thermodynamique regroupée calculées pendant la simulation compositionnelle dans chaque maille et à chaque pas de temps. Si Nrg est le nombre de composants de la représentation thermodynamique regroupée, on dispose donc de Nrg équations pour déterminer n+1
coefficients. Une condition nécessaire est donc que Nrg soit au moins égal à n+l.
Pour l'équation d'état de Peng-Robinson à deux paramètres, il faut donc une
représentation thermodynamique regroupée à au moins trois composants.
Une fois le coefficient ADo et les n coefficients ADp calculés, on les utilise pour calculer les constantes d'équilibre des composants de la représentation thermodynamique détaillée (Nrd composants) en appliquant l'équation (1) aux Nrd composants avec leurs propres paramètres fixes de caractérisation dans la
représentation thermodynamique détaillée.
En utilisant a) les constantes d'équilibre pour la représentation thermodynamique détaillée ainsi déterminées, b) les flux entre chaque maille et aux puits, c) la fraction vapeur dans chaque maille issus de la simulation compositionnelle regroupée, et, d) la composition détaillée globale dans chaque maille et aux puits injecteurs en début de pas de temps, on estime alors dans chaque maille la composition détaillée de chaque phase hydrocarbure au pas de temps t et la composition détaillée globale de chaque maille au pas de temps (t+1) suivant. L'un des intérêts de cette méthode est qu'il n'est pas nécessaire, dans l'étape d'éclatement (" delumping ") à chaque pas de temps, de résoudre l'équation d'état, que ce soit pour la représentation regroupée ou la représentation détaillée, ce qui constitue une économie de temps de calcul. Un inconvénient de la méthode est qu'elle n'est pas applicable aux simulations de type B.O, puisque les équations d'état
pratiques ont au moins deux paramètres.
Par la demande de brevet WO 98/5710726, on connaît également une méthode pour prédire l'évolution de la composition détaillée au cours du temps à partir des calculs de flux effectués dans une simulation de type B.O. o les phases
hydrocarbures sont décrites par seulement deux composants (L) et (V).
L'inconvénient de la méthode est qu'elle requiert l'utilisation de l'équation d'état de la représentation détaillée à chaque pas de temps, et donc d'être coûteuse en temps de
calcul pendant l'étape de "delumping".
La méthode selon l'invention La méthode selon l'invention combine une simulation prévisionnelle du comportement du gisement par une modélisation de type B.O. avantageuse puisqu'elle ne nécessite aucun recours à une équation d'état, et une étape d'éclatement ou de "delumping" qui ne nécessite à aucun pas de temps une
quelconque résolution d'une équation d'état.
Dans une modélisation de type B.O. on introduit, entre autres paramètres, les paramètres thermodynamiques des fluides (tels que viscosité, facteur de volume, densité, rapport de dissolution gaz-huile, etc.) dans un fichier des données d'entrée du simulateur (sous forme de tables, et ou par corrélation, comme fonction de la pression, de la température si celle-ci varie) et, le cas échéant (pour le modèle de type B.O. étendu), d'un paramètre additionnel conservant une mémoire de la
composition du gaz tel par exemple que la densité du gaz).
L'idée essentielle est ici d'introduire simplement dans les données d'entrée les (n+l) colonnes ou corrélations supplémentaires correspondant aux différences ADp, par exemple en les incluant dans un fichier annexe au fichier principal utilisé par la simulation. I est également possible d'introduire 2 x (n+1) colonnes ou corrélations supplémentaires correspondant aux paramètres Dp de chaque phase. Une condition pratique pour que des données puissent être entrées sous forme de table est que les
dépendances fonctionnelles introduites soient suffisamment continues.
Par là, on s'affranchit de la nécessité du calcul des constantes d'équilibre par
le modèle.
La méthode de simulation selon l'invention permet de prédire en fonction du temps et dans au moins une zone thermodynamique, une composition détaillée d'un fluide produit par un réservoir et plus particulièrement une composition détaillée d'un fluide contenu dans et produit par un gisement pétrolier dans lequel sont
implantés un ou plusieurs puits de production.
Elle est caractérisée en ce que: - on discrétise le gisement par un maillage, chaque maille contenant une ou plusieurs phases dont au moins une phase non aqueuse; - on détermine la variation de paramètres thermodynamiques des phases non aqueuses nécessaires à une modélisation de type " Black Oil " (B.O.) au cours d'étapes d'un chemin thermodynamique suivi par les fluides dans le gisement (tels par exemple que la viscosité des phases, le rapport de dissolution du gaz dans l'huile, les facteurs de volume de l'huile ou de compressibilité du gaz, etc.); - on définit les fluides par une représentation détaillée à Nrd composants et/ou pseudo-composants; - on détermine au moins une équation d'état à n paramètres permettant de simuler au moins durant la phase de préparation des données d'entrée (simulation qui est nécessaire strictement que durant cette phase de préparation des données d'entrée, et donc non nécessaire de façon répétitive au cours ou après la simulation de type B.O.), le comportement thermodynamique des fluides au cours des étapes d'un chemin thermodynamique suivi, permettant de générer, par zone ou domaine thermodynamique, n+1 relations fonctionnelles supplémentaires (sous forme de tables ou de corrélations); - on traduit le comportement thermodynamique de chaque phase non aqueuse en des données d'entrée convenant pour une représentation thermodynamique de type B.O., données d'entrée que l'on complète par les données d'entrée supplémentaires convenant à l'opération d'éclatement - on réalise une modélisation de type B.O., permettant de déterminer dans chaque maille et à des pas de temps consécutifs, des caractéristiques thermodynamiques de chaque phase non aqueuse et des données représentatives de déplacement des phases dans le gisement; et - on réalise une opération d'éclatement pour obtenir la composition détaillée des
fluides dans chacune des mailles.
Pour caractériser le chemin thermodynamrnique, on détermine par exemple les
variations de la composition des fluides et des propriétés des différentes phases.
Au cours de la simulation de type B.O., on détermine par exemple dans chaque maille des paramètres thermodynamiques tels que la pression (pj), la température (Tj) si celle-ci varie, les saturations des phases hydrocarbure liquide (Soj) et vapeur (Sgj), les débits d'injection ou de production, et, pour chaque paire
de mailles(j,h), les débits volumiques des phases liquide (uojh) et vapeur (Ugjh).
L'opération d'éclatement comporte par exemple une détermination des constantes d'équilibre à partir de données d'entrées spécifiques à l'opération d'éclatement, et une conversion de résultats exprimés en volume par des résultats
utilisables dans des équations de conservation des quantités molaires ou massiques.
Présentation sommaire des figures
D'autres caractéristiques et avantages de la méthode selon l'invention,
apparaîtront à la lecture de la description ci-après d'un exemple non limitatif de
réalisation, en se référant aux dessins annexés o: - les Fig. 1 à 3 montrent respectivement les variations en fonction de la pression de trois coefficients ADo, AD1, AD2, au cours de la vaporisation différentielle d'une huile de gisement simulée à l'aide de l'équation de Peng-Robinson à deux paramètres.
Description détaillée
Les différentes étapes de la méthode de modélisation permettant notamment de prédire en fonction du temps la composition détaillée de fluides produits par un gisement d'hydrocarbures en cours de production sont les suivantes: a) On représente tout d'abord le réservoir sous la forme d'un réseau de mailles (j) dont chacune constitue un volume élémentaire rempli de fluide(s) sous forme de une ou plusieurs phases, avec au moins une phase non aqueuse. Les phases non aqueuses sont dites phases hydrocarbure même si elles peuvent contenir certains composants autres que des hydrocarbures tel que l'azote, le dioxyde de carbone,
l'anhydride sulfureux.
b) On établit, d'une façon connue en soi, le comportement thermodynamique de la (des) phase(s) non aqueuse(s) nécessaire à une représentation thermodynamique de type B.O.. Il s'agit de décrire la dépendance de paramètres de phase avec la pression, la température si celle varie en cours de simulation et, éventuellement, pour les modélisations de type B. O. non strict, un indicateur d'évolution de la composition (par exemple le contenu en condensat dissous dans la phase gaz) au cours d'un chemin (évolution) thermodynamique représentatif de celui que va suivre le fluide dans le gisement. Le chemin thermodynamique le plus courant est celui de l'opération dite différentielle au cours de laquelle la pression du mélange est progressivement abaissée à température constante: sous le point de saturation, l'une des phases est évacuée à chaque palier de pression (pour la vaporisation différentielle d'une phase liquide, la phase gaz est évacuée); durant cette opération, les compositions et propriétés des phases (viscosités, masses volumiques, facteur de volume, rapport de dissolution gaz-huile
de la phase liquide, etc.) sont déterminées.
Il est à noter que l'on peut également traiter les cas qui requièrent l'usage de plusieurs représentations thermodynamiques " Black Oil " par exemple, si la composition en début de simulation varie avec la profondeur ou latéralement à l'intérieur du réservoir, ou si l'on peut distinguer plusieurs chemins thermodynamiques locaux durant la modélisation. Plusieurs zones ou domaines de variations des grandeurs thermodynamiques ou de composition, désignés souvent par les spécialistes comme zones thermodynamiques, pourront donc être définis et utilisés. L'origine des données sur le comportement du fluide de gisement peut être de nature expérimentale (expériences de laboratoire), elles peuvent être issues de mesures sur champ, ou prédites à partir de données frustres et d'hypothèses complémentaires en simulant numériquement le comportement du fluide de gisement. Parmi les propriétés de phase ainsi habituellement décrites en entrées des modélisations de type B.O. on trouve la viscosité des phases, le rapport de dissolution du gaz dans l'huile, le facteur de volume de l'huile, le facteur de
compressibilité du gaz.
c) On définit, avant la simulation B.O. et pour chaque zone ou domaine thermodynamique, le fluide par une représentation détaillée à Nrd composants et/ou pseudo-composants, d) Par zone thermodynamique, on détermine une équation d'état permettant de reproduire préalablement, avant la simulation B.O., le comportement thermodynamique du fluide au cours du chemin thermodynamique représentatif de
celui que va suivre le fluide dans le gisement.
e) On reproduit avec l'équation d'état les étapes (souvent des paliers de pression) du chemin thermodynamique de chaque zone, en particulier sous la pression de saturation o deux phases sont présentes. Lors de cette simulation thermodynamique préalable, on utilise, à la pression de saturation du fluide et en dessous, à chaque étape du chemin thermodynamique, les constantes d'équilibre obtenues à partir des résultats de la simulation, pour calculer le coefficient ADo et les n coefficients ADp de l'équation (1) correspondant à l'équation d'état choisie pour décrire le comportement du fluide hydrocarbure par exemple en minimisant la fonction: Nrd [ n (2) O(ADo, ADI,..,ADp,..,ADn) - = DO + ZADp 1-pi- Ln(ki) i=l p=l L'équation (1) ayant été obtenue à partir d'une équation sur la fugacité du constituant i dans une phase o apparaissent les coefficients Do et les n coefficients Dp propres à la phase: È(3) L{ =nDo+Z Dpfpn (3) Ln Xi= Do + Dp FI"pi, Ix+p) p=l on peut alternativement déterminer ces coefficients qui apparaissent comme des paramètres de phases, l'inconvénient étant d'augmenter le volume des données d'entrées. On dispose alors à chaque étape du chemin thermodynamique, prévu comme représentatif de celui que va suivre le fluide dans le gisement, de la dépendance en fonction de la pression, de la température et, le cas échéant, d'autres indicateurs nécessaires, des propriétés des phases nécessaires à la simulation B.O., et, sous la pression de saturation, du coefficient ADo et des n coefficients ADp (ou des n+l paramètres Do et Dp, p variant de 1 à n, par phase) qui seront utilisés pour calculer
les constantes d'équilibre lors de l'étape de "delumping" décrite ciaprès.
Les variations des différences ADp (p= 0 à 2) pour une opération
différentielle sur une huile de gisement simulée à l'aide de l'équation de Peng-
Robinson à deux paramètres, sont illustrées aux Fig. 1 à 3. On peut constater que le comportement de ces paramètres est suffisamment continu pour être aisément introduit sous forme de table, ou alternativement par une corrélation simple en fonction de la pression. Par exemple, dans le cas présent pour ADo: ADo = 1.406231 ( 1 - LN(P/Psat) [1 + 0.038 x LN(P/Psat)]) Au stockage pour utilisation ultérieure de l'évolution des valeurs de ces coefficients en fonction des abscisses de pression, température, et indicateur compositionnel, il peut être pratique d'ajouter l'évolution des valeurs de masse volumique de l'huile dans les conditions thermodynamiques de chaque étape du
chemin thermodynamique suivi, tant au-dessus de la pression de saturation que au-
dessous. On peut choisir également ou alternativement de stocker d'autres paramètres
des phases huile et/ou gaz, mais ce n'est pas indispensable.
f) On traduit, d'une façon connue en soi, le comportement thermodynamique de la (des) phase(s) hydrocarbure(s) en une ou plusieurs représentation(s) thermodynamique(s) de type B.O. dans les données d'entrée de la modélisation. Ces dépendances sont décrites au-dessus de la pression de saturation du mélange d'hydrocarbures et au-dessous de la pression de saturation, soit point par point dans des tables avec une méthode d'interpolation et, le cas échéant, une méthode d'extrapolation, soit par des corrélations analytiques, soit par un mélange de tables et de corrélations. g) On réalise, de façon connue en soi, une simulation de type B.O., ladite simulation permettant de calculer au moins dans chaque maille (j) et à des pas de m temps consécutifs (m, m+1, etc.), la pression (pj), la température (T7) (si celle-ci varie), les saturations des phases hydrocarbure liquide (Soj) et vapeur (Sgj), les débits d'injection ou de production, et, pour chaque paire de mailles (j,h), les débits m Il
volumiques des phases liquide (uojh) et vapeur (Ugih).
h) On détermine, à chaque pas de temps (m) et dans chaque maille (j), m les constantes d'équilibre (kij) du composant i à partir des dépendances des coefficients (ADo, AD...ADp...ADn) de l'équation (1), ou des coefficients (Do, Dl,.,Dp,..Dn) de l'équation (3) écrite par phase, établies en e) avec la pression (pj), la température (Tj) et avec, le cas échéant en cas de modélisation B.O. non stricte un indicateur de composition comme le contenu en condensat dissous dans la phase gaz. i) On détermine, à chaque pas de temps (m), la fraction vaporisée (0j) dans chaque maille (j), soit à partir des résultats de la simulation de type B. O. soit, pour une meilleure précision en cas de coexistence des phases huile et gaz, en résolvant l'équation de Rachford-Rice connue des spécialistes qui est appliquée par exemple dans le document suivant: Rachford, H.H.Jr et Rice, J.D.; "Procedure for Use of Electronic Digital
Computers in Calculating Flash Vaporization Hydrocarbon Equilibrium", J. Pet.
Technol., 1952, 14, 19, à partir des fractions molaires de chaque composant i dans la composition détaillée m globale (zij) du fluide hydrocarbure dans la maille (j) au pas de temps (m) m m Nrd zii (kii-1) (4) m m- 0 i= + (kij -1) oj j) On estime la composition détaillée de chaque phase hydrocarbure, à chaque pas de temps (m) et dans chaque maille (j), en utilisant les relations suivantes m m (5) Xm Z1ii (kii - 1) Pur laphasehuile = 1 + (kj -1) ou m m m m kii Zii (kii - 1) (6) Yij - + (kj -1) O pour la phase gaz k) On estime la densité molaire de la phase huile,oj dans la maille (j) au pas de temps (m) par exemple à partir de la masse volumique de l'huile poj, dans la maille (j) au pas de temps (m), par la relation: m (7) oj = Mo, MMo étant la masse molaire de la phase liquide
hydrocarbure (dite 'huile').
m Dans l'équation (7), la masse volumique poj peut donc être issue des résultats au pas de temps (m) de la simulation B.O. ou obtenues à partir des données d'entrée spécifiques à l'étape d'éclatement stockées lors de la simulation thermodynamique préalable. MMo peut être calculée à partir des résultats de la simulation B.O. si la masse molaire du constituant hydrocarbure lourd est définie par défaut ou dans les données d'entrées utilisateur, ou à partir de: Nid (8) MMo7 = x MM x,étant la masse molaire du composant ou i=l m
pseudo-composant i, et les quantités xij ayant été calculées dans l'étape précédente.
1) On estime la densité molaire de la phase gaz ogj dans chaque maille (j) à chaque pas de temps (m), soit à partir de données d'entrée spécifiques à l'étape d'éclatement stockées lors de la simulation thermodynamique préalable, soit à partir m m des résultats de la simulation B.O. en pression Pj, température Tj, par exemple avec, lorsque le facteur de compressibilité du gaz fait partie des données d'entrée de la représentation thermodynamique de type B.O. la relation suivante: m (9), gj = Zgj R T, Zgj étant le facteur de compressibilité du gaz, R la constante des gaz parfaits ou, par exemple, si c'est le facteur de volume du gaz qui fait partie des données d'entrée de la représentation thermodynamique de type B.O. avec la relation suivante: m 1 (10) g j = m, Vmolst étant le volume molaire du gaz en conditions Bgj Vmolst
standard et BgJ étant le facteur de volume du gaz.
m) On évalue dans chaque maille (j) à chaque pas de temps (m) les débits molaires des phases liquide (Uoh) et vapeur (u.h) à partir des débits volumiques des phases liquide (uojh) et vapeur (ugh) obtenus lors de la simulation B.O. et les densités
molaires obtenues dans les étapes k) et 1).
n) On estime au pas de temps m+l la fraction molaire de chaque composant ou pseudo composant i dans la composition détaillée globale (z I) du composant ou pseudo composant i dans la composition détaillée globale (z ij) du fluide hydrocarbure dans la maille (j) connaissant la fraction molaire (zij) au pas de temps m à partir des équations ci-après: Zij Nj At (yij' Sg + xij, Soj)- At A (Yij' ugjh + xj uojh) h+1 J(j) (11)zi = m+l mI m m (12)NJ = NJ - At (S + Soj)- At A (Ugh + uOjh) h J(j) m l'équation (12) traduisant le bilan molaire total hydrocarbure sur la maille (j), NJ, étant la quantité molaire hydrocarbure contenue dans la maille (j) au pas de temps (m) en tenant compte des échanges de matière avec toutes les mailles (h) voisines de
(j) qui forment l'ensemble J(j).
m Dans l'équation (11), l'écriture des termes (yij) et (xij), - dans lesquels j'=j pour un débit de la maille (j) vers la maille (h) ou dans le puits, et j'=h pour un débit de la maille (h) vers la maille (j), et j' correspondant au fluide injecté dans le cas des puits d'injection, S étant alors négatif -, suppose implicitement l'utilisation d'un schéma simple amont pour les flux compositionnels. Une écriture plus générale de ces termes est (yijh) et (Xijh) o xijh et Yih décrivent les compositions des phases liquide et gaz, obtenues d'une façon connue en soi, s'écoulant entre les mailles (j) et (h). Les termes source (débits molaires) (Sm) et (Somj) sont obtenus à partir des termes source (débits volumiques) (Si) et (Soj) dans la simulation B.O. en
reproduisant pour les termes source les étapes j) à m).
Les étapes h) à n) décrivent l'opération d'éclatement ou "delumping". Les détails fournis sont basés sur des équations de conservation en moles; on aurait pu
tout aussi bien se baser sur des équations de conservation en masse.
Dès lors que l'on sait décrire la composition détaillée dans chaque maille du modèle de simulation de type B.O. à l'instant t, qui peut être en particulier l'instant initial, on sait grâce à la présente méthode, décrire l'évolution de la composition détaillée dans chaque maille au cours du procédé d'exploitation modélisé dans un modèle de type " Black Oil >".
Claims (4)
1) Méthode de modélisation permettant de prédire en fonction du temps la composition détaillée des fluides produits par un gisement souterrain en cours de production, dans au moins une zone thermodynamique définie dans la modélisation, caractérisée en ce que: - on discrétise le gisement par un maillage, chaque maille contenant une ou plusieurs phases dont au moins une phase non aqueuse; - on détermine la variation de paramètres thermodynamiques des phases non aqueuses nécessaires à une modélisation de type " Black Oil " (B.O.) au cours d'étapes d'un chemin thermodynamique suivi par les fluides dans le gisement; - on définit les fluides par une représentation détaillée à Nrd composants et/ou pseudo- composants; - on détermine au moins une fonction d'état à n paramètres permettant de simuler au moins dans la phase de préparation des données d'entrée, le comportement thermodynamique des fluides au cours des étapes du chemin thermodynamique suivi de façon à générer, dans chaque dite zone thermodynamique o l'on souhaite effectuer une opération d'éclatement, des données d'entrée supplémentaires; - on traduit le comportement thermodynamique de chaque phase non aqueuse en des données d'entrée convenant pour une représentation thermodynamique de type B.O., données d'entrée que l'on complète par les dites données d'entrée supplémentaires convenant à une opération d'éclatement; - on réalise une modélisation de type B.O., permettant de déterminer dans chaque maille et à des pas de temps consécutifs, des caractéristiques thermodynamiques de chaque phase non aqueuse et des données représentatives de déplacement des phases dans le gisement; et - on réalise une opération d'éclatement pour obtenir la composition détaillée des
fluides dans chacune des mailles.
2) Méthode selon la revendication 1, caractérisée en ce que l'on caractérise le chemin thermodynamique en déterminant les variations de la composition des fluides et des propriétés des différentes phases.
3) Méthode selon la revendication 1 ou 2, caractérisée en ce que l'on détermine dans chaque maille au cours de la simulation de type B.O., des paramètres thermodynamiques tels que la pression (pj), la température (Tj) , si celle-ci varie, les saturations des phases hydrocarbure liquide (Soj) et vapeur (Sgj), les débits d'injection ou de production, et, pour chaque paire de mailles(j,h), les débits m m
volumiques des phases liquide (uojh) et vapeur (Ugih).
4) Méthode selon l'une des revendications précédentes, caractérisée en ce que
l'opération d'éclatement comporte une détermination des constantes d'équilibre à partir de données d'entrées spécifiques à l'opération d'éclatement, et une conversion de résultats exprimés en volume par des résultats utilisables dans des équations de
conservation des quantités molaires ou massiques.
) Méthode selon l'une des revendications précédentes, caractérisée en ce que
l'on utilise chaque fonction d'état pour générer n+1 relations fonctionnelles supplémentaires sous forme de tables de données ou de corrélations, que l'on inclut
dans les données d'entrée.
Priority Applications (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0009008A FR2811430B1 (fr) | 2000-07-10 | 2000-07-10 | Methode de modelisation permettant de predire en fonction du temps la composition detaillee de fluides porudits par un gisement souterrain en cours de production |
NL1018475A NL1018475C2 (nl) | 2000-07-10 | 2001-07-06 | Modelleerwerkwijze die het als functie van de tijd voorspellen van de gedetailleerde samenstelling van flu´da mogelijk maakt, die voortkomen uit een ondergronds reservoir onder productie. |
GB0116518A GB2369217B (en) | 2000-07-10 | 2001-07-06 | Modelling method allowing to predict as a function of time the detailed composition of fluids produced by an underground reservoir under production |
US09/899,105 US6985841B2 (en) | 2000-07-10 | 2001-07-06 | Modelling method allowing to predict as a function of time the detailed composition of fluids produced by an underground reservoir under production |
CA2352621A CA2352621C (fr) | 2000-07-10 | 2001-07-09 | Methode de modelisation permettant de predire en fonction du temps de la composition detaillee de fluides produits par un gisement souterrain en cours de production |
NO20013391A NO319817B1 (no) | 2000-07-10 | 2001-07-09 | Modelleringsmetode som tillater a forutsi som en funksjon av en tid den detaljerte sammensetning av fluider produsert av et undergrunnsreservoar under produksjon |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0009008A FR2811430B1 (fr) | 2000-07-10 | 2000-07-10 | Methode de modelisation permettant de predire en fonction du temps la composition detaillee de fluides porudits par un gisement souterrain en cours de production |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2811430A1 true FR2811430A1 (fr) | 2002-01-11 |
FR2811430B1 FR2811430B1 (fr) | 2002-09-06 |
Family
ID=8852325
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR0009008A Expired - Fee Related FR2811430B1 (fr) | 2000-07-10 | 2000-07-10 | Methode de modelisation permettant de predire en fonction du temps la composition detaillee de fluides porudits par un gisement souterrain en cours de production |
Country Status (6)
Country | Link |
---|---|
US (1) | US6985841B2 (fr) |
CA (1) | CA2352621C (fr) |
FR (1) | FR2811430B1 (fr) |
GB (1) | GB2369217B (fr) |
NL (1) | NL1018475C2 (fr) |
NO (1) | NO319817B1 (fr) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109184677A (zh) * | 2018-09-26 | 2019-01-11 | 西南石油大学 | 用于非均质互层砂体的储层评价方法 |
Families Citing this family (46)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9269043B2 (en) | 2002-03-12 | 2016-02-23 | Knowm Tech, Llc | Memristive neural processor utilizing anti-hebbian and hebbian technology |
US8909580B2 (en) | 2011-01-26 | 2014-12-09 | Knowmtech, Llc | Methods and systems for thermodynamic evolution |
WO2004049216A1 (fr) * | 2002-11-23 | 2004-06-10 | Schlumberger Technology Corporation | Procede et systeme pour simulations de reseaux d'installations de surface et de reservoirs integres |
GB2434235A (en) * | 2002-11-23 | 2007-07-18 | Schlumberger Technology Corp | Method and System of Black Oil Delumping |
AU2006344398B2 (en) * | 2005-10-06 | 2011-05-19 | Logined B.V. | Method, system and apparatus for numerical black oil delumping |
US20100132450A1 (en) * | 2007-09-13 | 2010-06-03 | Pomerantz Andrew E | Methods for optimizing petroleum reservoir analysis |
WO2009035918A1 (fr) * | 2007-09-13 | 2009-03-19 | Schlumberger Canada Limited | Procédé pour optimiser l'analyse de réservoir de pétrole |
CA2702965C (fr) | 2007-12-13 | 2014-04-01 | Exxonmobil Upstream Research Company | Partitionnement parallele adaptatif de donnees sur une simulation de reservoir utilisant une grille non structuree |
US8365831B2 (en) * | 2007-12-18 | 2013-02-05 | Exxonmobil Upstream Research Company | Determining connectivity architecture in 2-D and 3-D heterogeneous data |
CA2705340C (fr) | 2007-12-21 | 2016-09-27 | Exxonmobil Upstream Research Company | Procede et appareil pour analyser des donnees tridimensionnelles |
CN101903805B (zh) * | 2007-12-21 | 2013-09-25 | 埃克森美孚上游研究公司 | 沉积盆地中的建模 |
WO2009094064A1 (fr) * | 2008-01-22 | 2009-07-30 | Exxonmobil Upstream Research Company | Analyse de connectivité en dynamique |
US7920970B2 (en) * | 2008-01-24 | 2011-04-05 | Schlumberger Technology Corporation | Methods and apparatus for characterization of petroleum fluid and applications thereof |
US7822554B2 (en) * | 2008-01-24 | 2010-10-26 | Schlumberger Technology Corporation | Methods and apparatus for analysis of downhole compositional gradients and applications thereof |
US9026418B2 (en) | 2008-03-10 | 2015-05-05 | Exxonmobil Upstream Research Company | Method for determining distinct alternative paths between two object sets in 2-D and 3-D heterogeneous data |
US8428919B2 (en) * | 2008-03-28 | 2013-04-23 | Exxonmobil Upstream Research Company | Computing a consistent velocity vector field from a set of fluxes |
US9733388B2 (en) | 2008-05-05 | 2017-08-15 | Exxonmobil Upstream Research Company | Systems and methods for connectivity analysis using functional objects |
US8352228B2 (en) * | 2008-12-23 | 2013-01-08 | Exxonmobil Upstream Research Company | Method for predicting petroleum expulsion |
US9552462B2 (en) * | 2008-12-23 | 2017-01-24 | Exxonmobil Upstream Research Company | Method for predicting composition of petroleum |
CA2750926A1 (fr) * | 2009-01-30 | 2010-08-05 | Chevron U.S.A. Inc. | Systeme et methode permettant de predire l'ecoulement d'un fluide dans des reservoirs souterrains |
EP2406663A1 (fr) | 2009-03-13 | 2012-01-18 | Exxonmobil Upstream Research Company | Procédé de prédiction d'un écoulement de fluide |
US8359163B2 (en) * | 2009-06-01 | 2013-01-22 | Schlumberger Technology Corporation | Estimating fluid flow in a reservoir |
EP2491431A1 (fr) | 2009-10-20 | 2012-08-29 | Exxonmobil Upstream Research Company | Procédé d'évaluation quantitative de la connectivité pour des paires de puits à des fréquences variables |
WO2011159372A1 (fr) | 2010-06-15 | 2011-12-22 | Exxonmobil Upstream Research Company | Procédé et système permettant de stabiliser des procédés de formulation |
GB2498117B (en) * | 2010-09-28 | 2015-07-08 | Schlumberger Holdings | Methods for reservoir evaluation employing non-equilibrium compositional gradients |
US8972232B2 (en) | 2011-02-17 | 2015-03-03 | Chevron U.S.A. Inc. | System and method for modeling a subterranean reservoir |
WO2012141740A1 (fr) | 2011-04-14 | 2012-10-18 | Exxonmobil Upstream Research Company | Méthode de préparation d'un échantillon à base de pétrole pour analyse élémentaire et isotopique |
CA2853297C (fr) | 2011-11-11 | 2019-12-24 | Exxonmobil Upstream Research Company | Procede de determination de la position, de la dimension et de la composition de fluide d'une accumulation d'hydrocarbures de sous-sol |
WO2013071189A1 (fr) | 2011-11-11 | 2013-05-16 | Exxonmobil Upstream Research Company | Procédé et système pour surveillance de gisement utilisant des données d'isotopes et/ou de gaz nobles agglomérés |
US10677022B2 (en) | 2012-06-15 | 2020-06-09 | Landmark Graphics Corporation | Systems and methods for solving a multi-reservoir system with heterogeneous fluids coupled to a common gathering network |
US9031822B2 (en) | 2012-06-15 | 2015-05-12 | Chevron U.S.A. Inc. | System and method for use in simulating a subterranean reservoir |
CN103902752B (zh) * | 2012-12-27 | 2018-04-03 | 中国地质大学(北京) | 一种新型黑油、组分一体化混合数值模拟方法及装置 |
CN103760306B (zh) * | 2014-01-15 | 2016-01-20 | 中国人民解放军后勤工程学院 | 一种受限空间油气热爆炸发生的预警和控制方法 |
RU2634793C1 (ru) | 2014-03-07 | 2017-11-03 | Эксонмобил Апстрим Рисерч Компани | Способ разведки и система для обнаружения углеводородов по водяному столбу |
US10472960B2 (en) * | 2014-12-30 | 2019-11-12 | Schlumberger Technology Corporation | Estimating contamination during focused sampling |
WO2017058738A1 (fr) * | 2015-09-28 | 2017-04-06 | Schlumberger Technology Corporation | Déroulement des opérations de simulation basée sur un réseau |
US10571605B2 (en) | 2016-03-28 | 2020-02-25 | Exxonmobil Upstream Research Company | Method and system for modeling and simulating a fractured reservoir |
US10329905B2 (en) | 2016-04-07 | 2019-06-25 | Baker Hughes, A Ge Company, Llc | Method to estimate the influence of pore-size distribution on phase equilibrium of multi-component hydrocarbon systems in unconventional shale gas and oil reservoirs |
US10846447B2 (en) | 2016-04-29 | 2020-11-24 | Exxonmobil Upstream Research Company | Method and system for stacking fracture prediction |
US10572611B2 (en) | 2016-04-29 | 2020-02-25 | Exxonmobil Upstream Research Company | Method and system for characterizing fractures in a subsurface region |
US10712472B2 (en) | 2016-04-29 | 2020-07-14 | Exxonmobil Upstresm Research Company | Method and system for forming and using a subsurface model in hydrocarbon operations |
CA3051877C (fr) | 2017-02-28 | 2021-07-06 | Exxonmobil Upstream Research Company | Applications d'isotopes metalliques dans l'exploration, le developpement et la production d'hydrocarbures |
CN107145671B (zh) * | 2017-05-09 | 2019-04-12 | 中国石油大学(华东) | 一种油藏数值模拟方法及系统 |
CN107605474B (zh) * | 2017-09-15 | 2020-11-13 | 中石化石油工程技术服务有限公司 | 一种随钻预测气层产量的方法及装置 |
CN107908842B (zh) * | 2017-11-06 | 2021-01-29 | 齐鲁工业大学 | 润滑油动力黏度-温度模型建立方法 |
CN108229095A (zh) * | 2018-02-02 | 2018-06-29 | 华北理工大学 | 油中溶解气体体积分数的预测方法及终端设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5710726A (en) * | 1995-10-10 | 1998-01-20 | Atlantic Richfield Company | Semi-compositional simulation of hydrocarbon reservoirs |
FR2756044A1 (fr) * | 1996-11-18 | 1998-05-22 | Inst Francais Du Petrole | Methode pour constituer un modele representatif d'ecoulements polyphasiques dans des conduites de production petroliere |
FR2775094A1 (fr) * | 1998-02-18 | 1999-08-20 | Elf Exploration Prod | Methode de simulation pour predire en fonction du temps une composition detaillee d'un fluide produit par un reservoir |
WO1999042937A1 (fr) * | 1998-02-18 | 1999-08-26 | Elf Exploration Production | Methode de simulation pour predire en fonction du temps une composition detaillee d'un fluide produit par un reservoir |
WO2000037898A1 (fr) * | 1998-12-18 | 2000-06-29 | Exxonmobil Upstream Research Company | Procede d'evaluation des proprietes d'un fluide a plusieurs composants a l'aide de pseudocomposants |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6101447A (en) * | 1998-02-12 | 2000-08-08 | Schlumberger Technology Corporation | Oil and gas reservoir production analysis apparatus and method |
-
2000
- 2000-07-10 FR FR0009008A patent/FR2811430B1/fr not_active Expired - Fee Related
-
2001
- 2001-07-06 NL NL1018475A patent/NL1018475C2/nl not_active IP Right Cessation
- 2001-07-06 GB GB0116518A patent/GB2369217B/en not_active Expired - Fee Related
- 2001-07-06 US US09/899,105 patent/US6985841B2/en not_active Expired - Lifetime
- 2001-07-09 NO NO20013391A patent/NO319817B1/no not_active IP Right Cessation
- 2001-07-09 CA CA2352621A patent/CA2352621C/fr not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5710726A (en) * | 1995-10-10 | 1998-01-20 | Atlantic Richfield Company | Semi-compositional simulation of hydrocarbon reservoirs |
FR2756044A1 (fr) * | 1996-11-18 | 1998-05-22 | Inst Francais Du Petrole | Methode pour constituer un modele representatif d'ecoulements polyphasiques dans des conduites de production petroliere |
FR2775094A1 (fr) * | 1998-02-18 | 1999-08-20 | Elf Exploration Prod | Methode de simulation pour predire en fonction du temps une composition detaillee d'un fluide produit par un reservoir |
WO1999042937A1 (fr) * | 1998-02-18 | 1999-08-26 | Elf Exploration Production | Methode de simulation pour predire en fonction du temps une composition detaillee d'un fluide produit par un reservoir |
WO2000037898A1 (fr) * | 1998-12-18 | 2000-06-29 | Exxonmobil Upstream Research Company | Procede d'evaluation des proprietes d'un fluide a plusieurs composants a l'aide de pseudocomposants |
Non-Patent Citations (4)
Title |
---|
C LEIBOVICI ET AL: "A Consistent Procedure for Pseudo-Component Delumping", FLUID PHASE EQUILIBRIA,NL,ELSEVIER SCIENTIFIC PUBLISHING COMPANY, AMSTERDAM, vol. 117, 1996, pages 225 - 232, XP002083980 * |
J A TRANGENSTEIN ET AL: "Mathematical Structure of Compositional Reservoir Simulation", SIAM JOURNAL ON SCIENTIFIC AND STATISTICAL COMPUTING,US,PHILADELPHIA, PA, vol. 10, no. 5, 1 September 1989 (1989-09-01), pages 817 - 845, XP002083981, ISSN: 0196-5204 * |
ODEH A S: "COMPARISON OF SOLUTIONS TO A THREE-DIMENSIONAL BLACK-OIL RESERVOIR SIMULATION PROBLEM", JOURNAL OF PETROLEUM TECHNOLOGY,US,AIME, DALLAS, TX, vol. 33, no. 1, 1981, pages 13 - 25, XP000910660, ISSN: 0149-2136 * |
WANG P ET AL: "COMPARISON OF FLASH CALCULATIONS FOR IMPES COMPOSITIONAL RESERVOIR SIMULATORS", IN SITU,DEKKER, NEW YORK, NY,,US, vol. 21, no. 3, August 1997 (1997-08-01), pages 223 - 237, XP000910535, ISSN: 0146-2520 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109184677A (zh) * | 2018-09-26 | 2019-01-11 | 西南石油大学 | 用于非均质互层砂体的储层评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CA2352621A1 (fr) | 2002-01-10 |
GB2369217A (en) | 2002-05-22 |
US6985841B2 (en) | 2006-01-10 |
CA2352621C (fr) | 2011-10-18 |
NO20013391L (no) | 2002-01-11 |
FR2811430B1 (fr) | 2002-09-06 |
GB0116518D0 (en) | 2001-08-29 |
NO20013391D0 (no) | 2001-07-09 |
NO319817B1 (no) | 2005-09-19 |
US20020016703A1 (en) | 2002-02-07 |
GB2369217B (en) | 2004-07-21 |
NL1018475C2 (nl) | 2002-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR2811430A1 (fr) | Methode de modelisation permettant de predire en fonction du temps la composition detaillee de fluides porudits par un gisement souterrain en cours de production | |
Killough et al. | Fifth comparative solution project: evaluation of miscible flood simulators | |
Al-Meshari et al. | New Strategic Method to Tune Equation-of-state for compositional simulation | |
EP2775092B1 (fr) | Procédé d'exploitation d'un gisement d'hydrocarbures contenant des composés organo-soufrés au moyen d'un modèle thermo-cinétique et d'une simulation de réservoir compositionnelle | |
EP1462605A1 (fr) | Méthode de pseudoisation et d'éclatement pour décrire des fluides hydrocarbones | |
Barroux et al. | Forecasting of H2S Production due to Aquathermolysis reactions | |
El-Banbi et al. | New modified black-oil correlations for gas condensate and volatile oil fluids | |
EP3118758B1 (fr) | Procede d'exploitation d'un gisement d'hydrocarbures contenant des composes organo-soufres au moyen d'un modele thermo-cinetique et d'une simulation de reservoir compositionnelle | |
Yang et al. | Predict Reservoir Fluid Properties from Advanced Mud Gas Data | |
Leibovici et al. | Method for delumping the results of compositional reservoir simulation | |
Gaganis et al. | An accurate method to generate composite PVT data for black oil simulation | |
Al-Meshari et al. | Validation of splitting the hydrocarbon plus fraction: first step in tuning equation of state | |
CA2421096C (fr) | Methode pour modeliser des conditions de floculation d'asphaltenes dans des fluides hydrocarbones apparentes a un fluide de reference | |
Nederlof et al. | Understanding the H2S variations in the Upper Arab Formation in Abu Dhabi | |
Faissat et al. | Fluid modelling consistency in reservoir and process simulations | |
Gessner et al. | An internally consistent procedure to characterize single carbon number fractions for phase equilibrium of petroleum mixtures: Application to Brazilian pre-salt reservoir fluids | |
Levaché et al. | Deep offshore fluid evaluation/connectivity study | |
Sedarat et al. | A quality control protocol for gas condensate fluid samples | |
Rubin | Simulating gravity drainage and reinfiltration with a subdomain-dual-permeability hybrid fracture model | |
Promzelev et al. | Peculiarities of Identification of Reservoir Fluids Properties of Two-Phase With Oil Rim and Gas Cap Deposits: From Sampling to Justification of the Parameters for Calculating Reserves and PVT Data for Hydrodynamic Simulation of Field Development | |
Imo-Jack et al. | Splitting the Heptanes-Plus fraction of Condensate and Volatile systems: Application of the Gamma distribution Model | |
AL-Muftah | Unified Field Compositional Fluid Model for the Bahrain Field | |
Beilby et al. | Distribution of hydrogen sulphide in the Buzzard Field | |
Canel et al. | Compositional material balance: its application to the development of an oil and gas field with retrograde condensation | |
Johnson-Achilike et al. | Thermodynamic Modelling for Mapping Wax Deposition Envelope for an Identified Niger Delta Oilfield |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
CD | Change of name or company name | ||
PLFP | Fee payment |
Year of fee payment: 17 |
|
PLFP | Fee payment |
Year of fee payment: 18 |
|
PLFP | Fee payment |
Year of fee payment: 19 |
|
ST | Notification of lapse |
Effective date: 20200305 |