Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Chapitre II Analyse Nodale

Télécharger au format doc, pdf ou txt
Télécharger au format doc, pdf ou txt
Vous êtes sur la page 1sur 13

Chapitre II Analyse nodale

II.1. Introduction :

L'analyse nodale est ut il is ée pour évaluer un système de production global


(commençant par la pression statique du réservoir, et finissant avec celle de séparateur).
C'est une technique d'optimisation employée pour analyser les problèmes de production et
pour améliorer les performances du système de plusieurs de puits. Cette méthode est
utilisée intensivement dans les gisements d'huile et de gaz, depuis qu'elle a été introduite
par Gilbert dans les années 50.

II.2. Importance de l’analyse Nodale :


N'importe quel puits de production est foré et complété pour extraire le pétrole, le gaz
ou l'eau du réservoir à la surface. Pour surmonter des pertes de charge dans le système de
transport pendant l'extraction de ces fluides, de l'énergie est exigée. Le système de
production peut être simple ou peut inclure beaucoup de composants dans lesquels les
pertes de charge se produisent.
La figure (2.1) représente un schéma simplifié de l’écoulement de l’effluent lors de la
production et les différentes pertes de charges qui peuvent avoir lieu dans tout le système du
réservoir jusqu'au séparateur. Il peut être subdivisé comme suit :
1. L’écoulement dans le milieu poreux,
2. Complétion (stimulation, perforation, et gravel pack),
3. L’écoulement dans le tubing vertical ou dirigé (restriction, vanne de Sécurité),
4. L’écoulement en surface dans les réseaux de collecte (duse, pipes, vannes, etc..).

16
Chapitre II Analyse nodale

Figure 2.1: Pertes de charge possible dans un système de production

II.3. Procédure d’application de l’analyse nodale :

L’analyse nodale est appliquée pour analyser les performances des systèmes qui se
composent de plusieurs éléments agissants les uns sur les autres. Le procédé consiste à
choisir un nœud dans le puits et diviser le système au niveau de ce nœud. Les nœuds
utilisés sont montrés sur la figure 2.2
Tous les composants en amont du nœud composent la section Inflow, alors que la section
Outflow est compose par tous les éléments en aval du nœud.
Une relation entre le débit et la chute de pression doit être établie pour chaque élément
du système. Le débit à travers le système est déterminé une fois que les conditions
suivantes sont satisfaites au niveau du nœud :
 Le débit entrant égal à celui sortant.
 Une seule pression peut exister.
Une fois qu’on sélectionne le nœud, la pression à ce dernier est déterminée par :

Inflow : 𝑃node = 𝑃r − ∆p…………………………………………………. (2.1)

17
Chapitre II Analyse nodale

Outflow : 𝑃node = 𝑃sep + ∆p…………………………………………………. (2.2)


La chute de pression dans n’importe quelle composante varie avec le débit q une
représentation de la pression en fonction du débit produit deux courbes dont
l’intersection donnera un point qui vérifie les deux conditions citées au-dessus ; c’est le
point de fonctionnement du system.

L’effet du changement dans n’importe quelle composante peut être analysé en


recalculant la pression du nœud en fonction du débit en utilisant les nouvelles
caractéristiques de la composante. Si le changement se fait au niveau des composantes
ascendantes, la courbe d’outflow ne change pas. Avec le changement de l’une des deux
courbes, on aura un autre point de fonctionnement c’est à dire une nouvelle capacité
d’écoulement même si les pressions fixées se change à cause de l’épuisement ou un
changement des conditions de séparation.

La procédure est la suivante :


 Choisir les composants à optimiser,

 Sélectionner l'endroit du nœud qui ressentira l'effet du changement dans le composant

 Développer les expressions pour l'inflow et l’outflow,

 Obtenir les données nécessaires pour la construction des IPR,

 Déterminer l'effet du changement des caractéristiques des composants choisis


en traçant l'inflow ou l'outflow.

Figure. 2.2: Courbe d’outflow et inflow.

18
Chapitre II Analyse nodale

II.4. Les différentes positions du nœud :

Figure. 2.3: Les différentes positions du nœud

 Nœud 1:( séparateur)


Le choix du nœud au niveau du séparateur permet d'étudier l'effet de la pression de
séparateur sur le fonctionnement du puits.

 Nœud 2: (duse)

L'emplacement de nœud 2 (duse) nous permet d'étudier L'effet de la duse et de contrôler


le débit de production. Les données nécessaires sont:
IPR.
Chute de pression en fonction de débit.
L'équation 𝑃𝑤f = f (𝐺𝐿��,q , ��).
Avec :
𝜕 : Diamètre de la duse 1/64 in
q : Débit de production.
 Nœud 3: (tête du puits)

19
Chapitre II Analyse nodale

Le choix du nœud au niveau de la tête de puits permet d'étudier l'effet du diamètre des
flowlines sur la performance du puits.
Les données nécessaires sont:
 IPR (Inflow Performance Relationship) du puits.
 Chute de pression dans le tubing en fonction de débit.
 Chute de pression dans la collecte en fonction de débit.
 Pression de séparation.
La figure (2.4) montre l'effet de trois diamètres de collecte sur la performance du puits ;
L’augmentation du diamètre de la collecte (𝐷1 → 𝐷2 → 𝐷3) est suivi par une augmentation
de débit de production (q1 → q2 → q3)

Figure. 2.4: L'effet du diamètre du flowline.


 Nœud 6: (fond du puits)
Le choix du nœud dans le fond du puits nous permet d'étudier l'effet de l'IPR (nœud au
fond du puits) et de diamètre de tubing sur la performance du puits qui est l’un des
objectives de ce travail. Les données nécessaires sont:
 Courbe d’IPR mesurée dans le trou de puits.
 Chute de pression dans le tubing en fonction de débit.
 Chute de pression dans le réseau de collecte en fonction de débit.
 Pression de séparation.
 Nœud 8 : (réservoir)
Le choix du nœud dans le réservoir permet de savoir l'effet de l'épuisement du
réservoir sur la performance du puits. Les données nécessaires sont:

20
Chapitre II Analyse nodale

 IPR.
 Chute de pression dans le tubing en fonction du débit.
 Chute de pression dans le réseau de collecte en fonction du débit.
 Pression de séparation.
La figure (2.5) montre l'effet du déclin de la pression dans le réservoir. Le déclin de
pression dans le réservoir s'accompagne avec un déclin de la production.

Figure. 2.5: L'effet du déclin de la pression du réservoir.

II.5. Objectifs de l'analyse nodale:


Les objectifs de l'analyse nodale sont :
- Déterminer le débit avec lequel un puits d'huile ou du gaz produira avec la
considération des limitations de la géométrie et de la complétion du puits (d'abord
par déplétion naturelle).
- Déterminer sous quelles conditions d'écoulement (qui peuvent être liées au temps)
un puits cessera de produire.
- Définir le moment opportun pour l'installation d’un mécanisme d’activation
artificielle et aider dans le choix de la méthode d’activation.
- Optimiser le système pour produire avec un débit planifié.
- Vérifier chaque composante dans le système de production (déterminer si elle
affecte négativement ou positivement débit de production).

21
Chapitre II Analyse nodale

- Permettre l'identification rapide par (le personnel de la gestion et l'intérieur) des


manières avec lesquelles on peut améliorer la production.

II.6. Inflow performance relationship (IPR) :


L’écoulement du réservoir vers le fond est une partie essentielle du système
d’écoulement dans le puits.
Cet écoulement dépend de la chute de pression dans le réservoir. La relation entre le
débit et la chute de pression dans le milieu poreux peut être très complexe et dépend de
paramètres, tels que les propriétés pétro physiques de la roche et les propriétés des fluides,
régime d’écoulement, la saturation des roches en fluide, la compressibilité des fluides, la
formation endommagée ou stimulée….
L’IPR (Inflow Performance Relationship) montre donc la relation entre la pression de
fond dynamique et le débit de production provenant du réservoir.

Plusieurs méthodes ont été élaborées pour construire la courbe IPR des puits à huile et à
gaz, parmi ces méthodes on a :

II.6.1. LOI DE DARCY :


La loi fondamentale du mouvement liquide dans les milieux poreux est la loi de Darcy.
L'expression mathématique développée par Darcy en 1956 déclare que la vitesse d'un
fluide homogène dans un milieu poreux est proportionnelle au gradient de pression, et
inversement proportionnel à la viscosité liquide.
Cette loi est donnée par l'expression suivante:

V= q /A = (k/µ) ∂P/∂ x………..………………………….……………(2.3)

La loi de Darcy s'applique seulement quand les conditions suivantes existent :


 Écoulement à faible vitesse (laminaire).

 Écoulement permanent.

 Formation homogène.

 Pas de réaction entre le fluide et la formation.


22
Chapitre II Analyse nodale

Pour l'écoulement turbulent, qui se produit à des vitesses plus élevées, une modification
spéciale de l'équation de Darcy est nécessaire. Quand l'écoulement turbulent existe,
l'application de l'équation de Darcy peut avoir comme conséquence des erreurs très graves.

II.6.2. Méthode de Vogel:


J.V Vogel a élaboré, en 1968, une formule empirique pour prévoir l’indice de
productivité pour des puits ou la pression de fond statique est inférieure à la pression de
bulle. Il a calculé l’IPR à partir de plusieurs puits qui ont des pressions de fonds
différentes et des formations de différents paramètres pétro-physiques contenant des
fluides différents en PVT, et il a tracé la courbe IPR, à partir de cette courbe il a fait
une approche numérique est-il a déduit son équation :
……………..(2.4)

q0 : Débit entrant (inflow rate) correspondant à Pwf .

q0max  : Débit entrant (inflow rate) maximum, correspondant à une pression dynamique

Pwf : Pression de fond dynamique.


PR : Pression du réservoir approximative.

II.6.3. Les principales corrélations utilisées pour l'IPR:


II.6.3.1. Corrélation de Jones:
L’équation de Jones pour le gaz est une forme modifiée de l’équation de Darcy qui prend
en considération les pertes de charge dans l’écoulement laminaire et turbulent, elle est
exprimée sous la forme :
𝑃𝑟 − 𝑃𝑤𝑓 = 𝐴 𝑄2 + 𝐵𝑄
"A" et " B" : sont calculés à partir les propriétés de réservoir, comme ils peuvent être
déterminés à partir les essais multi-rate ; en traçant la courbe ΔP /Q =f (Q) en coordonnées
cartésiennes, on obtient une ligne droite ayant une pente " A " et une intersection "B" quand le
débit (Q) égale à zéro.

23
Chapitre II Analyse nodale

II.6.3.2. Corrélation de Forcheimer :


Cette corrélation exprime l’IPR en fonction des coefficients de perte de charge dans
l’écoulement turbulent et l’écoulement laminaire.
𝑃𝑟 − 𝑃𝑤𝑓 = 𝐴 𝑄2 + 𝐵𝑄
Tel que :
A : coefficient de perte de charge due à l’écoulement turbulent.
B : coefficient de perte de charge due à l’écoulement laminaire.

II.6.3.3. Back pressure :


Sa forme est :
𝑄 = 𝑐(𝑃𝑟2 − 𝑃𝑤𝑓2)𝑛
Où "c" et "n" sont déterminés à partir de pression et des propriétés de réservoir.

II.6.3.4. C and N:
C’est la même forme que l’équation back pressure mais elle diffère dans la manière de
l'obtention .
des constante c et n :
𝑄 = 𝑐(𝑃𝑟2 − 𝑃𝑤𝑓2)𝑛

Où "c" et "n" sont déterminés en traçant la courbe de Q en fonction de (𝑃𝑟2 − 𝑃𝑤𝑓2)


sur un papier log-log .Et "n" varie de 0.5 pour un écoulement turbulent à 1 pour un
écoulement laminaire.
II.6.3.5. Multi-rate Jones :
Plus de10 points de test peuvent être introduits pour adaptés la corrélation de Jones pour
les gaz qui est exprimé sous la forme suivante :
𝑃𝑟2 − 𝑃𝑤𝑓2 = 𝐴 𝑄2 + 𝐵𝑄

24
Chapitre II Analyse nodale

II.6.3.6. Modified isochronal :


Les tests isochrones modifiés sont généralement utilisés pour les puits à gaz parce qu’ils ont
besoin moins de temps et d’argents pour produire les même résultats que ceux du test
isochrone. Ce modèle est basé sur celui de back pressure :
𝑄 = 𝑐(𝑃𝑟2 − 𝑃𝑤𝑓2)𝑛

𝑠𝑖 𝑛 = 1 𝑎𝑙𝑜𝑟𝑠 𝑐=0,703 𝑘 /𝑢𝑔 𝑇 𝑍 ln⁡(0,606 𝑟𝑒 / 𝑟𝑤)

II.6.3.7. Forcheimer avec pseudo pressure:


C’est un modèle récent qu’est basé sur le modèle standard de Forcheimer, en remplaçant le
terme pression par la fonction pseudo pression :
Ѱr − Ѱwf = 𝐴 𝑄2 + 𝐵𝑄

II.7. Pertes de charge dans le tubing (Outflow) :


L’équation de perte de charge est composée par trois termes dans n’importe quel
écoulement :

 Le terme d’élévation où le terme statique.


 Le terme de frottement (friction).
 Le terme d’accélération
Pour un écoulement vertical et incliné, le terme d’élévation est le plus important.
L’équation d’écoulement des fluides, qui est valable pour n’importe quel fluide
(monophasique où poly phasique) et pour n’importe quelle inclinaison du pipe

Perte de charge totale = Perte de charge par friction + Perte de charge par élévation +
Perte de charge par accélération
II.7.1. Développement des Méthodes de Calcul des Pertes de Charge dans le
tubing:
Plusieurs méthodes empiriques existent pour le calcul des pertes de charges liées au tubing
de production. Parmi ces méthodes on peut citer les modèles suivants :
 Modèle homogène.
 Modèle de glissement.

25
Chapitre II Analyse nodale

 Modèle généralisé
a. Modèles homogènes
Ne tient pas compte de HL pour le calcul de la densité, qui est évaluée par simple
reconstitution des phases. Et considère une même vitesse d'écoulement gaz et liquide. Ces
modèles ont des domaines de validité assez réduit, et celui de Poettmann et Carpenter est
limité à une concentration de WOR= 20% lors de la présence de trois (03) phases.

b. Modèles de glissement

On les considéré plus précis que les précédent, et admettent que le liquide et le gaz se
déplacent à des vitesses différentes.

Ce modèle comprend ceux qui calculent la densité du mélange en fonction de HL.

Celui-ci peut être lié ou indépendant du coefficient de pertes de charge qui dépend lui-
même des propriétés des fluides.

La configuration d'écoulement est prise en compte, la masse volumique est calculée en


fonction de HL et le facteur de frottement dépend surtout de la phase continue. Dans ce sens,
plusieurs théories ont vue le jour : DUNS et ROS (1962),

ROS (1961), ORKISZEWSKI (1967), BEGGS et BRILL (1973), ALL (1974).

c. Modèles généralisés:

Sont apparus pour combler les lacunes laissées par le modèle de glissement. Dans ces
modèles dits à deux fluides, chacune des phases est considérée séparément, ce qui conduit à
deux systèmes d'équations, concernant la conservation de masse, des quantités de mouvement,
et de l'énergie de chaque phase.

II.7.2. Les principales corrélations utilisées pour l'écoulement vertical :


II.7.2.1. Duns & Ros :
C'est le résultat du travail de laboratoire où le holdup de liquide et le régime d'écoulement
sont observés. Ils ont utilisés une carte du modèle d'écoulement pour déterminer la vitesse de
glissement (et par conséquent le holdup de liquide) et le facteur de frottement. Cette
corrélation est recommandée pour les puits où les grands rapports gaz-liquide et les vitesses
d'écoulement ont induit le comportement de régime d'écoulement.

II.7.2.2. Petroleum Expert :


Cette corrélation combine les meilleurs dispositifs des corrélations existantes.

26
Chapitre II Analyse nodale

Elle emploie la carte d'écoulement de Gould et Al et la corrélation de Hagedorn Brown dans


le slug flow, et des Duns & ROS pour le mist flow. Dans le régime transitoire, une
combinaison des résultats de ces modes d’écoulement sont employées.
II.7.2.3. Duns & Ros original :
C’est la méthode éditée originale, sans perfectionnements appliqués dans la corrélation
primaire de Duns & ROS. Cette dernière a été développée et optimisée pour l’utilisation avec
le condensat.
II.7.2.4. Duns & Ros modified :
C'est le résultat du travail de laboratoire où le holdup de liquide et le régime d'écoulement
sont observés. Ils ont utilisés une carte du modèle d'écoulement pour déterminer la vitesse de
glissement (et par conséquent le holdup de liquide) et le facteur de frottement. Cette
corrélation est recommandée pour les puits où les grands rapports gaz-liquide et les vitesses
d'écoulement ont induit le comportement de régime d'écoulement.

Tableau 2.1: Classification des corrélations de perte de charge

Méthode categorie
a
Poettmann and Carpenter
Baxendell and Thomas a
Fancher and Brown a
Hagedorn and Brown b
Gray b
Ashiem b
Duns and Ros c
Orkiszewski c
Aziz et al c
Chierici et al c
Beggs and Brill c

Conclusion :

27
Chapitre II Analyse nodale

Pour les réservoirs riches en gaz en condensat dont la pression du réservoir est inférieure à
la pression de rosée on utilise les modèles: Multi rate Jones, C and N pour tracer l'IPR et les
corrélations de Duns and Ros, Petroleum Expert établir le VLP.

28

Vous aimerez peut-être aussi