Material point method modeling in oil and gas reservoirs

Patent metadata

Field	Value
Publication number	US-9378310-B2
Application number	US-201213649655-A
Country	US
Kind code	B2
Filing date	Oct 11, 2012
Priority date	Oct 13, 2011
Publication date	Jun 28, 2016
Grant date	Jun 28, 2016

Abstract

A computer system and method of simulating the behavior of an oil and gas reservoir including changes in the margins of frangible solids. A system of equations including state equations such as momentum, and conservation laws such as mass conservation and volume fraction continuity, are defined and discretized for at least two phases in a modeled volume, one of which corresponds to frangible material. A material point model technique for numerically solving the system of discretized equations, to derive fluid flow at each of a plurality of mesh nodes in the modeled volume, and the velocity of at each of a plurality of particles representing the frangible material in the modeled volume. A time-splitting technique improves the computational efficiency of the simulation while maintaining accuracy on the deformation scale. The method can be applied to derive accurate upscaled model equations for larger volume scale simulations.

First claim

What is claimed is:

1. A method of operating a computer system to simulate the fluid and structural behavior of a volume of the earth near a wellbore, comprising the steps of:
retrieving, from a memory resource in the computer system, parameters defining properties of at least one fluid phase at each of a plurality of mesh nodes in a grid representative of a sub-surface formation in the volume to be simulated, the parameters comprising a velocity parameter for fluid at each mesh node;
retrieving, from a memory resource in the computer system, parameters defining properties for each of a plurality of particles of a frangible material within the volume, the parameters for each particle comprising a velocity parameter, and a damage parameter having possible states comprising undamaged and damaged;
selecting initial pressure conditions at each of the plurality of mesh nodes;
then, at each of a plurality of time steps over a simulation time period, solving a system of equations comprising a state equation for fluid flow at each of the plurality of mesh nodes and a state equation for particle behavior for each of the plurality of particles; and
processing the results of the solving steps to estimate fluid flow and changes in the sub-surface formation over the simulation time period;
wherein the solving step comprises:
for each of a first plurality of time steps separated by a first time period:
setting the velocity parameter to zero for each particle having a damage parameter in the undamaged state;
setting an elastic stress term to a constant for each particle having a damage parameter in the undamaged state;
then solving the state equation for fluid flow at each of the plurality of mesh nodes; and
solving the state equation for particle behavior for each of the plurality of particles; and

after the first plurality of time steps, for each of a second plurality of time steps separated by a second time period, the second time period being shorter than the first time period:
setting the velocity parameter for fluid in each mesh node to a constant value corresponding to a recent value of the velocity parameter in that mesh node;
then solving the state equation for particle behavior for each of a subset of the plurality of particles within a selected portion of the volume; and
then, for each particle in the subset, evaluating a damage model responsive to one or more stress values determined in the solving step, to determine the state of the damage parameter for that particle.

2. The method of claim 1 , wherein the retrieved parameters defining at least one fluid phase represent hydraulic fracturing fluid;
and wherein the retrieved parameters defining particles of a frangible material represent a gas shale.

3. The method of claim 1 , wherein the retrieved parameters defining at least one fluid phase represent heavy oil;
wherein the retrieved parameters defining particles of a frangible material represent unconsolidated sand;
wherein the state equation for fluid flow expresses momentum of the heavy oil at a mesh node in terms of a pressure gradient at the mesh node, viscous stress of the heavy oil at the mesh node, and material interaction with the unconsolidated sand at the mesh node;
wherein the state equation for particle behavior expresses momentum of unconsolidated sand in terms of a pressure gradient at the particle, elastic stress at the particle, viscous stress at the particle, and material interaction of the particle with heavy oil;
and wherein the system of equations further comprises equations corresponding to mass conservation laws for the heavy oil and the unconsolidated sand, and an equation corresponding to a volume fraction continuity constraint.

4. The method of claim 3 , wherein the processing step comprises:
volume averaging the system of equations to a first volume scale;
executing a simulation of the behavior of a volume corresponding to the first volume scale using the volume averaged system of equations;
setting values of constants in the volume averaged equations responsive to a comparison of the results of the executing step and a first instance of the step of solving the system of equations; and
upscaling the volume averaged system of equations to a larger volume scale; and
executing a simulation of the behavior of a volume corresponding to the larger volume scale using the upscaled system of equations.

5. The method of claim 1 , wherein the retrieved parameters defining at least one fluid phase represent heavy oil;
and wherein the retrieved parameters defining particles of a frangible material represent unconsolidated sand.

6. The method of claim 5 , wherein the step of evaluating the damage model comprises:
determining a shear stress value and a normal stress value resulting from the solving step;
responsive to the shear stress value exceeding a limit corresponding to the normal stress value, setting the damage parameter of the particle to the damaged state.

7. The method of claim 5 , further comprising:
after the second plurality of time steps, repeating the setting and solving steps for a third plurality of time steps separated by the first time period.

8. The method of claim 1 , wherein the step of evaluating the damage model comprises:
determining a shear stress value and a normal stress value resulting from the solving step;
responsive to the shear stress value exceeding a limit corresponding to the normal stress value, setting the damage parameter of the particle to the damaged state.

9. The method of claim 1 , further comprising:
after the second plurality of time steps, repeating the setting and solving steps for a third plurality of time steps separated by the first time period.

10. The method of claim 1 , wherein the retrieved parameters defining at least one fluid phase represent heavy oil;
and wherein the retrieved parameters defining particles of a frangible material represent unconsolidated sand;
wherein the state equation for fluid flow expresses momentum of the heavy oil at a mesh node in terms of a pressure gradient at the mesh node, viscous stress of the heavy oil at the mesh node, and material interaction with the unconsolidated sand at the mesh node;
and wherein the state equation for particle behavior expresses momentum of an unconsolidated sand particle in terms of a pressure gradient at the particle, elastic stress at the particle, viscous stress at the particle, and material interaction of the particle with heavy oil.

11. A computer system for simulating the fluid and structural behavior of a volume of the earth near a wellbore, comprising:
a processing unit for executing program instructions;
a memory resource, coupled to the processing unit, for storing data representative of properties of at least one fluid phase at each of a plurality of mesh nodes in a grid representative of a sub-surface formation in the volume to be simulated, and data representative of properties for each of a plurality of particles of a frangible material within the volume, each particle associated with a damage parameter having possible states comprising undamaged and damaged; and
program memory, coupled to the processing unit, for storing a computer program including program instructions that, when executed by the one or more processing units, causes the computer system to perform a sequence of operations comprising:
retrieving, from the memory resource, at least a portion of the data defining properties of at least one fluid phase at each of a plurality of mesh nodes in a grid representative of a sub-surface formation in the volume to be simulated;
retrieving, from the memory resource, at least a portion of the data definin

Assignees

• Vanderheyden William Brian
• Zhang Duan
• Los Alamos Nat Security Llc
• Bp Corp North America Inc

Inventors

• Vanderheyden William Brian
• Zhang Duan

Classifications

• G01V1/00

  Seismology; Seismic or acoustic prospecting or detecting · CPC title

• E21B43/00 Primary

  Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells (applicable only to water E03B) · CPC title

• E21B43/20

  Displacing by water · CPC title

• G06F2113/08

  Fluids · CPC title

• G06F30/20 Primary

  Design optimisation, verification or simulation (optimisation, verification or simulation of circuit designs G06F30/30) · CPC title