Rapid Flash Calculations for Compositional Modeling

Author:

Li Y.1,Johns R.T.1

Affiliation:

1. The U. of Texas at Austin

Abstract

Abstract The computational time for conventional flash calculations increasessignificantly with the number of components making it impractical for use inmany finely-grided compositional simulations.Previous research toincrease flash calculation speed has been limited to those with zero binaryinteraction parameters (BIPs) or approximate methods based on an eigenvalueanalysis of the binary interaction matrix.Practical flash calculations, however, nearly always have some nonzero BIPs.Further, the accuracy andspeed of the eigenvalue methods varies depending on the choice and number ofthe dominant eigenvalues. This paper presents a new method for significantly increasing the speed offlash calculations for any number of nonzero BIPs.The approach requiresthe solution of up to six reduced parameters regardless of fluid complexity orthe number of components and is based on decomposing the BIPs into twoparameters using a simple quadratic expression.The new approach is exactin that the equilibrium phase compositions for the same BIPs are identical tothose with the conventional flash calculation; no eigenvalue analysis isrequired. Further, the new approach eliminates the Rachford-Riceprocedure and is more robust than the conventional flash calculationprocedure.We demonstrate the new approach for several example fluids, andshow that speed up by a factor of about 3 – 20 is obtained over conventionalflash calculations depending on the number of components. Introduction Gas injection into oil reservoirs results in complex interactions of flowwith phase behavior that are often not modeled accurately by black oilsimulation. This is especially true for miscible or nearly misciblefloods where significant mass transfer occurs between the hydrocarbonphases.Such floods are best modeled by compositionalsimulation. A significant disadvantage of compositional simulation, however, is that itis more computationally intensive than black-oil simulation.The primaryreason for the increased computational time is the result of solving repeatedflash calculations with cubic equations-of-state (EOS).Research has shownthat EOS flash calculations can occupy 50–70% of total computational time incompositional simulation.[1,2] The use of fewer pseudocomponents can reduce the flash computation time, butfewer components results in less accuracy.[3–5] This is especially truein MCM displacements, in which miscibility is developed by a combinedcondensing/vaporizing (CV) drive process.[6–9]Fluid characterization bypseudocomponent models can be improved by tuning to the analytical MME orMMP,[9] but those models still require significant computational time even forfewer pseudocomponents. Another way to reduce computation time is to reduce the number of gridblocks. With coarse grids, however, numerical dispersion is large, which maycloud the results.[10] More recently, methods have been examined to find reduced parameters forflash calculations.Michelsen[11–13] significantly increased flashcalculation speed by finding three reduced parameters, regardless of the numberof components.His method, however, assumes zero BIPs, which is toorestrictive for real fluid characterization. Michelsen also gave apractical method for stability calculations using the tangent plane distance(TPD).[12] Jensen and Fredenslund[14] extended Michelsen's research by the addition oftwo reduced parameters for every column of the BIP matrix with at least onenonzero element.Thus, there are a total of five reduced parameters forthe case when one column contains a single nonzero BIP. All other BIPs, however, must be zero in their approach, which again limits itsusefulness. In 1988, Hendricks[15] used an eigenvalue analysis to identify the dominantBIPs in flash calculations. The BIP matrix is recalculated by settingsmall eigenvalues to zero, using some predetermined criterion.Theprocedure was applied later by Hendricks and van Bergen[16] usingNewton-Raphson (NR) iteration. Although the method is faster thanconventional flash calculations, the reduction in eigenvalues can lead tonon-physical BIPs in that the diagonal elements of the BIP matrix can havesmall nonzero values. Their method is only an approximation to theoriginal phase behavior characterization.

Publisher

SPE

Cited by 2 articles. 订阅此论文施引文献 订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献

同舟云学术

1.学者识别学者识别

2.学术分析学术分析

3.人才评估人才评估

"同舟云学术"是以全球学者为主线,采集、加工和组织学术论文而形成的新型学术文献查询和分析系统,可以对全球学者进行文献检索和人才价值评估。用户可以通过关注某些学科领域的顶尖人物而持续追踪该领域的学科进展和研究前沿。经过近期的数据扩容,当前同舟云学术共收录了国内外主流学术期刊6万余种,收集的期刊论文及会议论文总量共计约1.5亿篇,并以每天添加12000余篇中外论文的速度递增。我们也可以为用户提供个性化、定制化的学者数据。欢迎来电咨询!咨询电话:010-8811{复制后删除}0370

www.globalauthorid.com

TOP

Copyright © 2019-2024 北京同舟云网络信息技术有限公司
京公网安备11010802033243号  京ICP备18003416号-3