Continuous streamline trajectories on complex grids

Streamlines have been used for reservoir modeling and flow visualization in the petroleum industry and in computational fluid dynamics. When applied to the calculation of volumetric sweep and the identification of by-passed hydrocarbons for improved oil recovery, it is important that the velocity models that are used to trace trajectories across the cells of a grid are flux conservative. As such, the requirements on their tracing may be more stringent than in other disciplines. Flux conservation is also important at faults, at locally refined or coarsened embedded grid boundaries, and within unstructured grids, where the modeling of flow within a cell may not be consistent with the connection fluxes from adjacent cells. In such cases, additional degrees of freedom must be introduced to satisfy flux conservation. In this study, we introduce a flux conservative conforming cell face local boundary layer construction to resolve these inconsistencies. In contrast, solutions that rely upon spatial continuity of streamlines between elements are shown to not be flux conservative when these inconsistencies are present. The use of flux conservative conforming elements also allows the solution to be developed in local isoparametric coordinates, without explicit reference to cell or connection geometry. The solution has been implemented for both 3D corner point and for 2.5D PEBI grids. In all cases we utilize the lowest order Raviart–Thomas zeroth order velocity model, for which the trajectories and transit times may be obtained analytically. The results are demonstrated on a sequence of increasingly complex type, sector and full-field model applications.