Method for multiphase flow upscaling

A method of multiphase flow upscaling of a fine scale geologic model, comprising:

• from the fine scale geologic model, defining a local domain for generating a dynamic pseudo-function, the local domain including an upstream coarse grid domain adjacent a downstream coarse grid domain, and a concatenated domain adjacent the upstream coarse grid domain, the local domain being configured such that fluid flow to be modeled therethrough flows first through the concatenated domain, then through the upstream coarse grid domain, and lastly through the downstream coarse grid domain, wherein a coarse grid interface is defined between the upstream coarse grid domain and the downstream coarse grid domain;
• imposing a set of boundary conditions on an upstream face of the concatenated domain and a downstream face of the downstream coarse grid domain;
• imposing a no-flow boundary condition on faces of the concatenated domain, the upstream coarse grid domain, and the downstream coarse grid domain parallel to a direction of fluid flow;
• simulating fluid flow through the local domain for a plurality of time steps until the local domain is completely flooded;
• at each time step, recording a fine scale solution of saturation and flow rate of individual fluid phases across the coarse grid interface;
• constructing a pseudo-function by post-processing results of the simulation of fluid flow;
• and
• upscaling the fine scale geologic model to a coarse scale geologic model using the pseudo-function.

A system for multiphase flow upscaling of a fine scale geologic model, the system comprising:

• a processor;
• a storage medium;
• and
• a computer application accessible by the processor and stored on at least one of the storage medium and the processor, the computer application adapted to:
• from the fine scale geologic model, define a local domain to generate a dynamic pseudo-function, the local domain including an upstream coarse grid domain adjacent a downstream coarse grid domain, and a concatenated domain adjacent the upstream coarse grid domain, the local domain being configured such that fluid flow to be modeled therethrough flows first through the concatenated domain, then through the upstream coarse grid domain, and lastly through the downstream coarse grid domain, wherein a coarse grid interface is defined between the upstream coarse grid domain and the downstream coarse grid domain;
• impose a constant pressure boundary condition of P=1 and a saturation condition of Sd=1 on an upstream face of the concatenated domain;
• impose a constant pressure boundary condition of P=0 and a saturation condition of Sd=0 on a downstream face of the downstream coarse grid domain;
• impose a no-flow boundary condition on faces of the concatenated domain, the upstream coarse grid domain, and the downstream coarse grid domain parallel to a direction of fluid flow;
• simulate fluid flow through the local domain for a plurality of time steps until the local domain is completely flooded;
• at each time step, record a fine scale solution of saturation and flow rate of individual fluid phases across the coarse grid interface;
• construct a pseudo-function by post-processing results of the simulation of fluid flow;
• and
• upscale the fine scale geologic model to a coarse scale geologic model using the pseudo-function.

A computer program product having computer executable logic recorded on a tangible, machine readable medium, the computer program product comprising:

• code for defining, from the fine scale geologic model, a local domain for generating a dynamic pseudo-function, the local domain including an upstream coarse grid domain adjacent a downstream coarse grid domain, and a concatenated domain adjacent the upstream coarse grid domain, the local domain being configured such that fluid flow to be modeled therethrough flows first through the concatenated domain, then through the upstream coarse grid domain, and lastly through the downstream coarse grid domain, wherein a coarse grid interface is defined between the upstream coarse grid domain and the downstream coarse grid domain;
• code for imposing a constant pressure boundary condition of P=1 and a saturation condition of Sct=1 on an upstream face of the concatenated domain;
• code for imposing a constant pressure boundary condition of P=0 and a saturation condition of Sct=0 on a downstream face of the downstream coarse grid domain;
• code for imposing a no-flow boundary condition on faces of the concatenated domain, the upstream coarse grid domain, and the downstream coarse grid domain parallel to a direction of fluid flow;
• code for simulating fluid flow through the local domain for a plurality of time steps until the local domain is completely flooded;
• code for recording, at each time step, a fine scale solution of saturation and flow rate of individual fluid phases across the coarse grid interface;
• code for constructing a pseudo-function by post-processing results of the simulation of fluid flow;
• and
• code for upscaling the fine scale geologic model to a coarse scale geologic model using the pseudo-function.

A method of managing hydrocarbon resources, comprising:

• from a fine scale geologic model, defining a local domain for generating a dynamic pseudo-function, the local domain including an upstream coarse grid domain adjacent a downstream coarse grid domain, and a concatenated domain adjacent the upstream coarse grid domain, the local domain being configured such that fluid flow to be modeled therethrough flows first through the concatenated domain, then through the upstream coarse grid domain, and lastly through the downstream coarse grid domain, wherein a coarse grid interface is defined between the upstream coarse grid domain and the downstream coarse grid domain;
• imposing a constant pressure boundary condition of P=1 and a saturation condition of Sd=1 on an upstream face of the concatenated domain;
• imposing a constant pressure boundary condition of P=0 and a saturation condition of Sd=0 on a downstream face of the downstream coarse grid domain;
• imposing a no-flow boundary condition on faces of the concatenated domain, the upstream coarse grid domain, and the downstream coarse grid domain parallel to a direction of fluid flow;
• simulating fluid flow through the local domain for a plurality of time steps until the local domain is completely flooded;
• at each time step, recording a fine scale solution of saturation and flow rate of individual fluid phases across the coarse grid interface;
• constructing a pseudo-function by post-processing results of the simulation of fluid flow;
• upscaling the fine scale geologic model to a coarse scale geologic model using the pseudo-function;
• and
• managing hydrocarbon resources using the coarse scale geologic model.