1 / 42

PETE 310 Lectures # 32 to 34

PETE 310 Lectures # 32 to 34. Cubic Equations of State. …Last Lectures. Instructional Objectives. Know the data needed in the EOS to evaluate fluid properties Know how to use the EOS for single and multicomponent systems

arich
Download Presentation

PETE 310 Lectures # 32 to 34

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. PETE 310Lectures # 32 to 34 Cubic Equations of State …Last Lectures

  2. Instructional Objectives • Know the data needed in the EOS to evaluate fluid properties • Know how to use the EOS for single and multicomponent systems • Evaluate the volume (density, or z-factor) roots from a cubic equation of state for • Gas phase (when two phases exist) • Liquid Phase (when two phases exist) • Single phase when only one phase exists

  3. Equations of State (EOS) • Single Component Systems Equations of State (EOS) are mathematical relations between pressure (P) temperature (T), and molar volume (V). • Multicomponent Systems For multicomponent mixtures in addition to (P, T & V) , the overall molar composition and a set of mixing rules are needed.

  4. Uses of Equations of State (EOS) • Evaluation of gas injection processes (miscible and immiscible) • Evaluation of properties of a reservoir oil (liquid) coexisting with a gas cap (gas) • Simulation of volatile and gas condensate production through constant volume depletion evaluations • Recombination tests using separator oil and gas streams

  5. Equations of State (EOS) • One of the most used EOS’ is the Peng-Robinson EOS (1975). This is a three-parameter corresponding states model.

  6. T T c c CP e e sur sur T Pres Pres 2 v P 1 T 1 L 2 - P has es V L V Mo Mo la la r r V V o o lu lu m m e e PV Phase Behavior Pressure-volume behavior indicating isotherms for a pure component system

  7. Equations of State (EOS) • The critical point conditions are used to determine the EOS parameters

  8. Equations of State (EOS) • Solving these two equations simultaneously for the Peng-Robinson EOS provides • and

  9. Equations of State (EOS) Where and with

  10. EOS for a Pure Component CP CP e e e sur sur sur T T Pres Pres Pres 2 2 4 4 T T 1 1 v v 3 3 5 5 A A A 2 2 2 1 1 P P 1 1 - - L L 1 1 1 A A A 1 1 1 ¶ ¶ æ æ ö ö P P 1 1 1 0 0 2 2 2 > > ç ç ÷ ÷ 0 0 V V ~ ~ 0 ¶ ¶ è è ø ø V V T T 2 2 - - P P has has es 2 2 L L V V 7 7 6 6 Mo Mo la la r r V V o o lu lu m m e e

  11. Equations of State (EOS) • PR equation can be expressed as a cubic polynomial in V, density, or Z. with

  12. Equations of State (EOS) • When working with mixtures (aa) and (b) are evaluated using a set of mixing rules • The most common mixing rules are: • Quadratic for a • Linear for b

  13. Quadratic MR for a • where the kij’s are called interaction parameters and by definition

  14. Linear MR for b

  15. Example • For a three-component mixture (Nc = 3) the attraction (a) and the repulsion constant (b) are given by

  16. Equations of State (EOS) • The constants a and b can be evaluated using • Overall compositions zi with i = 1, 2…Nc • Liquid compositions xi with i = 1, 2…Nc • Vapor compositions yi with i = 1, 2…Nc

  17. Equations of State (EOS) • The cubic expression for a mixture is then evaluated using

  18. Analytical Solution of Cubic Equations • The cubic EOS can be arranged into a polynomial and be solved analytically as follows.

  19. Analytical Solution of Cubic Equations • Let’s write the polynomial in the following way Note: “x” could be either the molar volume, or the density, or the z-factor

  20. Analytical Solution of Cubic Equations • When the equation is expressed in terms of the z factor, the coefficients a1 to a3 are:

  21. Procedure to Evaluate the Roots of a Cubic Equation Analytically • Let

  22. Procedure to Evaluate the Roots of a Cubic Equation Analytically • The solutions are,

  23. Procedure to Evaluate the Roots of a Cubic Equation Analytically • If a1, a2 and a3 are real and if D = Q3 + R2 is the discriminant, then • One root is real and two complex conjugate if D > 0; • All roots are real and at least two are equal if D = 0; • All roots are real and unequal if D < 0.

  24. Procedure to Evaluate the Roots of a Cubic Equation Analytically where

  25. Procedure to Evaluate the Roots of a Cubic Equation Analytically where x1, x2 and x3 are the three roots.

  26. Procedure to Evaluate the Roots of a Cubic Equation Analytically • The range of solutions that are used for the engineer are those for positive volumes and pressures, we are not concerned about imaginary numbers.

  27. Solutions of a Cubic Polynomial From the shape of the polynomial we are only interested in the first quadrant.

  28. Solutions of a Cubic Polynomial • http://www.uni-koeln.de/math-nat-fak/phchem/deiters/quartic/quartic.html contains Fortran codes to solve the roots of polynomials up to fifth degree.

  29. Web site to download Fortran source codes to solve polynomials up to fifth degree

  30. EOS for a Pure Component CP CP e e e sur sur sur T T Pres Pres Pres 2 2 4 4 T T 1 1 v v 3 3 5 5 A A A 2 2 2 1 1 P P 1 1 - - L L 1 1 1 A A A 1 1 1 ¶ ¶ æ æ ö ö P P 1 1 1 0 0 2 2 2 > > ç ç ÷ ÷ 0 0 V V ~ ~ 0 ¶ ¶ è è ø ø V V T T 2 2 - - P P has has es 2 2 L L V V 7 7 6 6 Mo Mo la la r r V V o o lu lu m m e e

  31. Parameters needed to solve EOS • Tc, Pc, (acentric factor for some equations I.e Peng Robinson) • Compositions (when dealing with mixtures) • Specify P and T  determine Vm • Specify P and Vm  determine T • Specify T and Vm  determine P

  32. Tartaglia: the solver of cubic equations http://es.rice.edu/ES/humsoc/Galileo/Catalog/Files/tartalia.html

  33. Cubic Equation Solver http://www.1728.com/cubic.htm

  34. Equations of State (EOS) • Phase equilibrium for a single component at a given temperature can be graphically determined by selecting the saturation pressure such that the areas above and below the loop are equal, these are known as the van der Waals loops.

  35. Two-phase VLE • The phase equilibria equations are expressed in terms of the equilibrium ratios, the “K-values”.

  36. Dew Point Calculations • Equilibrium is always stated as: (i = 1, 2, 3 ,…Nc) • with the following material balance constrains

  37. Dew Point Calculations • At the dew-point (i = 1, 2, 3 ,…Nc)

  38. Dew Point Calculations • Rearranging, we obtain the Dew-Point objective function

  39. Bubble Point Equilibrium Calculations • For a Bubble-point

  40. Flash Equilibrium Calculations • Flash calculations are the work-horse of any compositional reservoir simulation package. • The objective is to find the fv in a VL mixture at a specified T and P such that

  41. Evaluation of Fugacity Coefficients and K-values from an EOS • The general expression to evaluate the fugacity coefficient for component “i” is

  42. Evaluation of Fugacity Coefficients and K-values from an EOS • The final expression to evaluate the fugacity coefficient using an EOS is.

More Related