### Abstract

Phase equilibrium calculations become computationally intensive in compositional simulation as the number of components and phases increase. Reduced methods were developed to address this problem where the binary interaction parameter (BIP) matrix is either approximated by spectral decomposition (Hendriks and van Bergen, 1992), or by using a two-parameter BIP formula (Li and Johns, 2006). Haugen and Beckner (2011) and Yan et al. (2011), however, recently stated that the spectral decomposition method, and by reference all reduced methods, are not as fast as previously reported in the literature. In this paper we present the first study that compares all eight reduced and conventional methods published to date using optimized code and compilers. The results show that the spectral decomposition method and its variants are not as fast as other reduced methods, and can be slower than the conventional approach when fewer than 20 components are used. The reason for the slow speed is the additional number of nested loops required in the code and the requirement that the code must allow for a variable number of eigenvalues. We show that the reduced method of Li and Johns (2006) and its variants, however, are significantly faster. Significant speed up in flash calculations is achieved for all fluids studied when more than six components are used. For example, for ten-component fluids, a speed up of 2 to 3 in the computational time for Newton-Raphson iterations are obtained compared to the conventional method based on minimization of Gibbs energy. We also demonstrate that the reduced method by Gorucu and Johns (2011, 2013), an offshoot of the Li and Johns (2006) reduced approach, is more robust than all other methods.

Original language | English (US) |
---|---|

Title of host publication | Society of Petroleum Engineers - SPE Reservoir Simulation Symposium 2013 |

Pages | 17-32 |

Number of pages | 16 |

State | Published - Aug 8 2013 |

Event | SPE Reservoir Simulation Symposium 2013 - The Woodlands, TX, United States Duration: Feb 18 2013 → Feb 20 2013 |

### Publication series

Name | Society of Petroleum Engineers - SPE Reservoir Simulation Symposium 2013 |
---|---|

Volume | 1 |

### Other

Other | SPE Reservoir Simulation Symposium 2013 |
---|---|

Country | United States |

City | The Woodlands, TX |

Period | 2/18/13 → 2/20/13 |

### Fingerprint

### All Science Journal Classification (ASJC) codes

- Geochemistry and Petrology
- Modeling and Simulation

### Cite this

*Society of Petroleum Engineers - SPE Reservoir Simulation Symposium 2013*(pp. 17-32). (Society of Petroleum Engineers - SPE Reservoir Simulation Symposium 2013; Vol. 1).

}

*Society of Petroleum Engineers - SPE Reservoir Simulation Symposium 2013.*Society of Petroleum Engineers - SPE Reservoir Simulation Symposium 2013, vol. 1, pp. 17-32, SPE Reservoir Simulation Symposium 2013, The Woodlands, TX, United States, 2/18/13.

**Comparison of reduced and conventional phase equilibrium calculations.** / Gorucu, S. E.; Johns, Russell Taylor.

Research output: Chapter in Book/Report/Conference proceeding › Conference contribution

TY - GEN

T1 - Comparison of reduced and conventional phase equilibrium calculations

AU - Gorucu, S. E.

AU - Johns, Russell Taylor

PY - 2013/8/8

Y1 - 2013/8/8

N2 - Phase equilibrium calculations become computationally intensive in compositional simulation as the number of components and phases increase. Reduced methods were developed to address this problem where the binary interaction parameter (BIP) matrix is either approximated by spectral decomposition (Hendriks and van Bergen, 1992), or by using a two-parameter BIP formula (Li and Johns, 2006). Haugen and Beckner (2011) and Yan et al. (2011), however, recently stated that the spectral decomposition method, and by reference all reduced methods, are not as fast as previously reported in the literature. In this paper we present the first study that compares all eight reduced and conventional methods published to date using optimized code and compilers. The results show that the spectral decomposition method and its variants are not as fast as other reduced methods, and can be slower than the conventional approach when fewer than 20 components are used. The reason for the slow speed is the additional number of nested loops required in the code and the requirement that the code must allow for a variable number of eigenvalues. We show that the reduced method of Li and Johns (2006) and its variants, however, are significantly faster. Significant speed up in flash calculations is achieved for all fluids studied when more than six components are used. For example, for ten-component fluids, a speed up of 2 to 3 in the computational time for Newton-Raphson iterations are obtained compared to the conventional method based on minimization of Gibbs energy. We also demonstrate that the reduced method by Gorucu and Johns (2011, 2013), an offshoot of the Li and Johns (2006) reduced approach, is more robust than all other methods.

AB - Phase equilibrium calculations become computationally intensive in compositional simulation as the number of components and phases increase. Reduced methods were developed to address this problem where the binary interaction parameter (BIP) matrix is either approximated by spectral decomposition (Hendriks and van Bergen, 1992), or by using a two-parameter BIP formula (Li and Johns, 2006). Haugen and Beckner (2011) and Yan et al. (2011), however, recently stated that the spectral decomposition method, and by reference all reduced methods, are not as fast as previously reported in the literature. In this paper we present the first study that compares all eight reduced and conventional methods published to date using optimized code and compilers. The results show that the spectral decomposition method and its variants are not as fast as other reduced methods, and can be slower than the conventional approach when fewer than 20 components are used. The reason for the slow speed is the additional number of nested loops required in the code and the requirement that the code must allow for a variable number of eigenvalues. We show that the reduced method of Li and Johns (2006) and its variants, however, are significantly faster. Significant speed up in flash calculations is achieved for all fluids studied when more than six components are used. For example, for ten-component fluids, a speed up of 2 to 3 in the computational time for Newton-Raphson iterations are obtained compared to the conventional method based on minimization of Gibbs energy. We also demonstrate that the reduced method by Gorucu and Johns (2011, 2013), an offshoot of the Li and Johns (2006) reduced approach, is more robust than all other methods.

UR - http://www.scopus.com/inward/record.url?scp=84881047759&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84881047759&partnerID=8YFLogxK

M3 - Conference contribution

AN - SCOPUS:84881047759

SN - 9781627480246

T3 - Society of Petroleum Engineers - SPE Reservoir Simulation Symposium 2013

SP - 17

EP - 32

BT - Society of Petroleum Engineers - SPE Reservoir Simulation Symposium 2013

ER -