Method and system for interpolating discontinuous functions in a subsurface model
A method and system are described for creating a subsurface model. In this method, discontinuous functions are interpolated in a subsurface model. The method extrapolates specified values over discontinuous manifolds by embedding the manifold in a higher-dimensional space that amplifies distances across discontinuities and thus eliminates special consideration to prevent extrapolation across discontinuities. The resulting subsurface model may be used in reservoir simulations and hydrocarbon operations.
1. A method for generating a subsurface model for a subsurface region comprising:
creating a mesh for a subsurface region in an original domain;
computing connectivity and compartmentalization of the mesh;
transforming the mesh into a transform domain wherein the transform domain is configured to exaggerate distance across discontinuities, wherein dimensionality of the created mesh plus the number of discontinuities is an upper bound on the dimensionality of the transform domain;
transforming measured data into the transform domain;
interpolating the measured data to the mesh in the transform domain;
transforming the interpolated measurements to the original domain; and
outputting the interpolated mesh and/or measurements in the original domain.
2. The method of claim 1 , wherein the transforming mesh into the transform domain further comprises preserving local distances between surface nodes in the mesh and exaggerating distances when nodes are separated by a discontinuity.
3. The method of claim 2 , wherein the distance is computed based on a Manhattan distance that is computed by propagating the distance from a first surface node to a second surface node.
4. The method of claim 2 , wherein the distance is computed based on a Euclidian distance that is computed between a first surface node to a second surface node.
5. The method of claim 1 , further comprising building a reservoir model from the interpolated mesh and assigning properties to the mesh elements of the interpolated mesh.
6. The method of claim 5 , wherein the properties comprise one or more of porosity, permeability and any combination thereof.
7. The method of claim 5 , further comprising simulating fluid flow within the reservoir model to create simulation results.
8. The method of claim 1 , wherein the measured data are either at least one component of the mesh coordinates or a categorical, discrete, or continuous property attributed to mesh nodes.
9. The method of claim 7 , further comprising managing hydrocarbon operations based on the simulation results.
10. A system for generating a subsurface model having one or more objects associated with a subsurface region, comprising:
a processor;
an input device in communication with the processor and configured to receive input data associated with a subsurface region;
memory in communication with the processor, the memory having a set of instructions, wherein the set of instructions, when executed, are configured to:
create a mesh for a subsurface region in an original domain;
compute connectivity and compartmentalize of the mesh;
transform mesh into a transform domain wherein the transform domain is configured to exaggerate distance across discontinuities, wherein dimensionality of the created mesh plus the number of discontinuities is an upper bound on the dimensionality of the transform domain;
transform measured data into the transform domain;
interpolate the measured data to a mesh in the transform domain;
transform the interpolated measurements to the original domain; and
output the interpolated mesh and/or measurements in the original domain.
11. The system of claim 10 , wherein the set of instructions are further configured to:
preserve local distances between surface nodes in the mesh and to exaggerate distances when nodes are separated by a discontinuity.
12. The system of claim 11 , wherein the set of instructions, when executed by the processor, are configured to: compute the distance based on a Manhattan distance that is computed by
propagating the distance from a first surface node to a second surface node.
13. The system of claim 11 , wherein the set of instructions, when executed by the processor, are configured to: compute the distance based on a Euclidian distance that is computed between a first surface node to a second surface node.
14. The system of claim 10 , wherein the set of instructions, when executed by the processor, are configured to: build a reservoir model from the interpolated mesh and assign properties to the mesh elements of the interpolated mesh.
15. The system of claim 14 , wherein the properties comprise one or more of porosity, permeability and any combination thereof.
16. The system of claim 14 , further comprising simulating fluid flow within the reservoir model to create simulation results.
17. The system of claim 16 , wherein the set of instructions, when executed by the processor, are further configured to: manage hydrocarbon operations based on the simulation results.