摘要 |
An apparatus and method are provided for modeling a reservoir simulation by solving a non-linear S-shaped function F = f(S) which is representative of a property S in a physical system, such saturation in a reservoir simulation. A Newton iteration (T) is performed on the function f(S) at Sto determine a next iterative value S. It is then determined whether Sis located on the opposite side of the inflection point Sfrom S. If Sis located on the opposite side of the inflection point from S, then Sis set to S, a modified new estimate. The modified new estimate, S, is preferably set to either the inflection point, S, or to an average value between Sand S, i.e., S= 0.5(S+S). The above steps are repeated until Sis within the predetermined convergence criteria. Also, solution algorithms are described for two-phase and three-phase flow with gravity and capillary pressure. |