WO2005033739A2 - Characterizing connectivity in reservoir models using paths of least resistance - Google Patents
Characterizing connectivity in reservoir models using paths of least resistance Download PDFInfo
- Publication number
- WO2005033739A2 WO2005033739A2 PCT/US2004/032942 US2004032942W WO2005033739A2 WO 2005033739 A2 WO2005033739 A2 WO 2005033739A2 US 2004032942 W US2004032942 W US 2004032942W WO 2005033739 A2 WO2005033739 A2 WO 2005033739A2
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- vertex
- vertices
- source
- paths
- path
- Prior art date
Links
- 238000000034 method Methods 0.000 claims abstract description 114
- 239000012530 fluid Substances 0.000 claims abstract description 57
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 50
- 230000005540 biological transmission Effects 0.000 claims abstract description 17
- 230000001186 cumulative effect Effects 0.000 claims description 8
- 238000005315 distribution function Methods 0.000 claims description 6
- 238000009472 formulation Methods 0.000 claims 4
- 239000000203 mixture Substances 0.000 claims 4
- 238000004088 simulation Methods 0.000 abstract description 64
- 229930195733 hydrocarbon Natural products 0.000 abstract description 13
- 239000004215 Carbon black (E152) Substances 0.000 abstract description 12
- 150000002430 hydrocarbons Chemical class 0.000 abstract description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 8
- 238000005094 computer simulation Methods 0.000 abstract description 2
- 239000012071 phase Substances 0.000 description 26
- 230000035699 permeability Effects 0.000 description 17
- 238000013459 approach Methods 0.000 description 15
- 230000008569 process Effects 0.000 description 13
- 238000004364 calculation method Methods 0.000 description 10
- 238000004519 manufacturing process Methods 0.000 description 10
- 230000037230 mobility Effects 0.000 description 9
- 238000005553 drilling Methods 0.000 description 8
- 239000011435 rock Substances 0.000 description 8
- 238000012512 characterization method Methods 0.000 description 7
- 238000004891 communication Methods 0.000 description 7
- 238000010168 coupling process Methods 0.000 description 7
- 238000005859 coupling reaction Methods 0.000 description 7
- 238000002347 injection Methods 0.000 description 7
- 239000007924 injection Substances 0.000 description 7
- 230000002829 reductive effect Effects 0.000 description 7
- 239000000243 solution Substances 0.000 description 7
- 230000008878 coupling Effects 0.000 description 6
- 230000000007 visual effect Effects 0.000 description 6
- 230000001419 dependent effect Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 239000011148 porous material Substances 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000009877 rendering Methods 0.000 description 4
- 230000036962 time dependent Effects 0.000 description 4
- 238000012546 transfer Methods 0.000 description 4
- 230000035508 accumulation Effects 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 3
- 230000001413 cellular effect Effects 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 3
- 230000000670 limiting effect Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 239000011347 resin Substances 0.000 description 3
- 229920005989 resin Polymers 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000006399 behavior Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 230000001788 irregular Effects 0.000 description 2
- 230000003340 mental effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000036961 partial effect Effects 0.000 description 2
- 230000002040 relaxant effect Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 240000000588 Hericium erinaceus Species 0.000 description 1
- 235000007328 Hericium erinaceus Nutrition 0.000 description 1
- 244000166071 Shorea robusta Species 0.000 description 1
- 235000015076 Shorea robusta Nutrition 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 239000008346 aqueous phase Substances 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000006735 deficit Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 125000001183 hydrocarbyl group Chemical group 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000008450 motivation Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- -1 physical Substances 0.000 description 1
- 230000000135 prohibitive effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000013341 scale-up Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000005204 segregation Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
- G01V2210/663—Modeling production-induced effects
Definitions
- the invention relates to numerical simulation of subsurface geological reservoirs. More specifically embodiments of the invention are related to computer modeling of the transmission of fluid (e.g., hydrocarbon natural resources and water) properties and/or reservoir properties within subsurface geological reservoirs.
- fluid e.g., hydrocarbon natural resources and water
- Reservoir simulation is used in the oil and gas production industry to assess the most profitable means to locate and operate wells (boreholes) to produce subsurface accumulations of oil and gas to the surface, where the fluids can be transported, refined, and sold.
- Reservoir simulation is of great interest because it infers the behavior of a real hydrocarbon-bearing reservoir from the performance of a mathematical or numerical model of that reservoir.
- the typical objective of reservoir simulation is to understand the complex chemical, physical, and fluid flow processes occurring in the reservoir sufficiently well to predict future behavior of the reservoir to maximize hydrocarbon recovery.
- Reservoir simulation calculations in such hydrocarbon systems are based on fluid flow through the reservoir being simulated. These calculations are performed with varying degrees of rigor, depending on the requirements of the particular simulation study and the capabilities of the simulation software being used.
- the physical system to be modeled is divided into cells (a set of which is called a grid or mesh), and the state variables that vary in space throughout the model are represented by sets of values for each cell.
- a reservoir simulation grid may be built by subdividing (discretizing or gridding) the volume of interest into a large number of polyhedral cells. The number of cells commonly ranges from a few tens of thousands to several hundreds of thousands.
- the volume of interest is defined areally and vertically by the extent of the oil and gas accumulation and of the water that is in pressure communication with the oil and gas. The area may be several square miles, and the thickness may be hundreds or even thousands of feet.
- the reservoir rock properties such as porosity and permeability are typically assumed to be constant inside a cell. Other variables such as fluid pressure and phase saturation are specified at specified points, sometimes called nodes, within the cell. A link between two nodes is called a "connection.” Fluid flow between two cells is typically modeled as flow along the connection between them.
- the state of a simulation cell may be defined by its pressure and its contents, i.e., the amounts of oil, gas, and water within the cell.
- the goal of simulation is to calculate the evolution through time of the states of all cells. This evolution is governed by the initial states and by the time-dependent removal of fluid from (production) or addition of fluid to (injection) the system by way of wells.
- the state of a cell changes in time because of fluid flow between pairs of neighboring cells or between a cell and a well. Fluid flows from high pressure to low pressure. Pressure gradients are induced by removing fluid from the reservoir (production) or adding fluid to the reservoir (injection) by way of wellbores that penetrate the porous and permeable rock. Within the reservoir, fluid converges on (flows toward) producing wellbores and diverges from (flows away from) injecting wellbores.
- the properties of cells and connections play a role in determining how the reservoir simulation model performs. These properties may be derived from a geologic reservoir characterization through processing techniques used in the art.
- Reservoirs are commonly segmented or compartmentalized by virtue of reduced permeability at fault surfaces and at stratigraphic sequence boundaries. These surfaces are often targeted for transmissibility adjustment and are just as often the source of inconsistent or erroneous adjustments that yield connectivity anomalies.
- the first method is finite difference simulation.
- the user runs a finite-difference reservoir simulation program and interrogates the results. This approach is slow (hours to tens of hours) and only informs the user about parts of the model that are well connected to wells.
- finite difference simulation there are three methods of interrogating results to assess connectivity.
- the first finite difference method for assessing connectivity is to visualize pressure results. This is the most common style of visual interrogation. Changes in pressure are associated with the flow of fluid. A common observation is that pressure depletion is observed in compartments previously thought to be isolated from (i.e., poorly connected to) producing wells. In such cases, the pressure field appears diffuse and does not offer any guidance as to the path that connects the pressure-depleted compartment to a producing well.
- the second finite difference method for assessing connectivity is to visualize tracer-to-component-ratio results.
- Reservoir simulation programs may offer the ability to associate named tracers with fluid components, e.g., water.
- fluid components e.g., water.
- tracer-ratio results are diffuse rather than specific. The fuzziness is exacerbated by the common situation that often only infrequent snapshots of tracerxomponent ratio results can be afforded. A low-frequency archive yields inadequately detailed information.
- the third finite difference method for assessing connectivity is to visualize flows.
- Customized software tools not available in standard software, make it possible to visualize statistics of internal flows. This method is still burdened by the need to actually conduct flow simulation and by the constraint of using actual wells. It works best for injector-producer pairs that are amenable to visualizing high flow rates.
- the second approach to assessing connectivity in numerical simulation uses mental tomography.
- the user builds a mental image of the permeability field based on viewing numerous slices through the model, where each slice displays the value of permeability using color.
- This technique may have been adequate for the coarsely gridded, nearly homogeneous, structured simulation grids that were built 20 years ago, but it is impractical for complex, unstructured, modern models.
- the third approach to assessing connectivity in numerical simulation uses streamline simulation. In this approach, the user runs a streamline flow simulation program and interrogates the results. This approach generally shows where fluid flows. However, it requires the a priori designation of source and/or sink wells, and it is slower than embodiments of the invention described herein to yield results.
- the fourth approach to assessing connectivity in numerical simulation uses geobody visualization.
- Customized software tools not available in standard software, make it possible to define and visualize geobodies with respect to user- defined criteria.
- a geobody is a connected component, which means for all possible pairs of cells that belong to the geobody, a path (sequence of connections) exists to connect the pair.
- User input for the "define geobodies" action consists of separate criteria for nodes and connections.
- the criteria are logical expressions (with a pass/fail outcome) that may range from simple to complex.
- the criteria are evaluated using node and connection properties. If a node fails the node criteria, the node is excluded from all geobodies, and the connections in which the failed node participates are treated as failed connections. If a connection fails the connection criteria, it is excluded from paths that may connect passing nodes. A node that participates in a failed connection is not marked as a failed node, because such a node may participate in other passing connections.
- the results of a "define geobodies" action are that each node is assigned values for two properties: a geobody membership index that identifies which geobody (connected component) the node belongs to and a geobody score.
- the geobody score is an aggregate measure of the size or importance of the geobody, such as the number of nodes or total pore volume of the geobody.
- the geobody membership index is commonly an integer rank, assigned in order of descending score, i.e., the highest scoring geobody is 1, the next lower scoring geobody is 2, etc.
- a special membership-index value such as zero or -1 marks failed nodes that belong to no geobodies. Connectivity among the failed nodes is not assessed.
- One embodiment of the invention includes a method of evaluating the transmission of a property within a subsurface geologic reservoir.
- the method may include providing a set of vertices representative of at least a portion of the reservoir and providing a plurality of edges representing property transmission paths between connected vertices within the set of vertices.
- the method may include associating with each edge a cost representative of the ability of a property to be transmitted across an edge from one vertex to another vertex within the set of vertices.
- the method may include selecting at least one source vertex from within the set of vertices and determining an extremum path between the source vertex and each of at least one and preferably at least two destination vertices that minimizes or maximizes the summed cost across the one or more edges included in the extremum path using a graph-theory single-source shortest-paths algorithm.
- One embodiment of the invention includes a computer-readable media tangibly embodying a program of instructions executable by a computer to evaluate the transmission of a property within a subsurface geologic reservoir.
- the program may include the following steps. a) reading a set of vertices representative of at least a portion of the reservoir; b) reading a plurality of edges, the edges representing property transmission paths between connected vertices within the set of vertices; c) associating with each edge a cost representative of the ability of a property to transmit across the edge from one vertex to another vertex within the set of vertices; d) reading a first vertex from within said set of vertices; and e) determining an extremum path between the first vertex and each of at least one and preferably at least two other vertices that minimizes or maximizes the summed connection cost across the one or more edges included in the extremum path using a graph-theory single-source shortest-paths algorithm.
- Fig. 1 is a graphical depiction of a simple set of two vertices and one edge.
- Fig. 2 is a block diagram of steps that may be completed according to some embodiments of the invention.
- Fig. 3 is a block diagram of steps that may be completed according to some embodiments of the invention.
- Fig. 4 is a depiction of a reservoir volume that has been discretized into a three by three by three group of volumetric cells.
- Fig. 5 is a depiction of the reservoir volume of Fig 4 that has been further abstracted into a set of vertices and edges thereby forming a graph-theory grid.
- Fig. 6 is a depiction of the grid of Fig 3 which includes a depiction of a lowest cost connected path from one vertex to a second vertex.
- Fig. 7 is a depiction of two paths from a source vertex to a destination vertex that have been generated on a 3D viewer.
- Fig. 8 is another depiction of two paths from a source vertex to a destination vertex that have been generated on a 3D viewer.
- Fig. 9 is a depiction that has been generated on a 3D viewer of two reservoir compartments that are separated by a fault.
- Figs. 10-12 are depictions of the two reservoir compartment/fault system of Fig. 9 that are depicted at varying threshold constraints.
- Figs. 13-16 are depictions of a source vertex and destination vertices at successively higher threshold values.
- Fig. 17 is a depiction of a shortest path from one vertex in one lens to a second vertex in a second lens.
- extremum means the maximum or minimum of some value or calculation.
- the extremum path of a group of paths may be either a path that has the highest cost or a path that has the lowest cost.
- the plural of extremum is extrema.
- Embodiments of the invention relate to reservoir simulation, including the numerical simulation of the transmission of fluid and/or reservoir properties within a geologic reservoir.
- Embodiments of the invention relate to numerical simulation of fluid flow in three dimensions (3D) of one, two, or three fluid phases in porous and permeable rock.
- the three fluid phases may be hydrocarbon liquid (oil), hydrocarbon vapor (gas), and an aqueous phase (water).
- Embodiments of the invention are also useful in simulating heat flow, chemical component transport and electrical current flow in a geologic reservoir.
- Embodiments of the invention apply graph-theory data structure and graph-theory algorithms to identify lowest-cost paths between any two locations in the model.
- these paths and any model properties on these paths can be visualized using a modeling software's 3D viewer. Numerical analysis of the paths can yield additional valuable information.
- a graph is a collection of points and lines connecting some subset of the points.
- the points of a graph are most commonly known as graph vertices, but may also be called “nodes” or simply “points.”
- the lines connecting the vertices of a graph are most commonly known as graph edges, but may also be called “arcs” or “lines.”
- Fig 1 depicting a first vertex 1 and a second vertex 2.
- the vertices in Fig 1 are connected by edge 3.
- the arrows on the ends of the edge 3 depict an edge that is capable of transmission from first vertex 1 to second vertex 2 or from second vertex 2 to first vertex 1.
- a graph is a mathematical abstraction that is useful for solving many kinds of problems.
- a graph consists of a set of vertices, and a set of edges, where an edge is something that connects two vertices in the graph. More precisely, a graph is a pair (V,E), where V is a finite set and E is a binary relation on V. V is called a vertex set whose elements are called vertices. E is a collection of edges, where an edge is a pair (u,v) with u,v in V.
- edges are ordered pairs, connecting a source vertex to a target vertex.
- edges are unordered pairs and connect the two vertices in both directions, hence in an undirected graph (u,v) and (v,u) are two ways of writing the same edge.
- Step 50 includes providing a set of vertices which may be representative of a geologic reservoir.
- Step 51 includes providing edges which may be representative of the ability of a property to transfer between connected vertices.
- Step 52 includes associating a cost with some or all of the edges. The cost may be representative of the ability of the property to transfer across an edge from one vertex that is connected to another vertex.
- Step 53 includes selecting a vertex from the set of vertices to act as a source vertex.
- Step 54 includes determining a low (or high) cost path to at least one other vertex using a graph theory single-source, shortest-paths algorithm.
- Figure 3 depicts another embodiment of the invention which includes more particular steps.
- This embodiment begins with the physical reservoir system 60 and from the reservoir system 60 determines a geologic reservoir characterization.
- the geologic reservoir characterization 60 may be in the form of a 3D geologic model or geologic maps.
- the gridding step 62 is accomplished in order to obtain a set of cells and connections that is representative of the physical reservoir system 60.
- a fine structured 3D geologic model grid may be translated into a coarser scale, unstructured grid that respects the geologic features (e.g. faults and shale lenses) of the geologic reservoir.
- An upscaling step 63 is completed to determine the cell properties and the connection properties 64.
- Cell properties may include, for example, porosity, permeability and facies, that is a characterization of the type of rock present.
- Connection properties may include, for example, transmissibility, the distance from a cell's center to a cell face, and a cell face's area. Transmissibility may be obtained from the cell's permeability, conductance and/or cell geometry.
- the system of cells and connections may be translated 65 into a system of vertices and edges to form a graph. Each edge or some number of the edges contained in the graph is assigned a cost or weight 66. The cost may be chosen to represent the ability of some property to move across the respective edge.
- One or more vertices are then selected 67 as a source vertex and the single- source, shortest path algorithm is then run to determine the lowest cost path from the source vertex to every other vertex which is connected to the source vertex 68.
- the data that is stored for each vertex on a shortest path may include the identity of the predecessor vertex and the summed total path cost up to that vertex 69.
- Various depictions of the paths may then be displayed on a 3-D viewer 70.
- a reservoir simulation grid 4 may be built by subdividing (discretizing or gridding) the reservoir volume of interest into a large number of polyhedral cells 5.
- the volume of interest may be defined areally and vertically by the extent of the oil and gas accumulation and of the water that is in pressure communication with the oil and gas.
- the area may be several square miles, and the thickness may be hundreds or even thousands of feet.
- the physical system to be simulated may be subdivided into a multiplicity of volumetric cells as in reservoir simulation. These cells can number from fewer than a hundred to millions.
- the discretization may be performed using finite difference, finite volume, finite element, or other methods known in the art that are based on dividing the physical system to be modeled into smaller units.
- the term "cell" is used, but it should be understood that if a simulation model uses the finite element method the term element would replace the term cell. Regardless of which of these methods is used, discretization reduces partial differential equations to a finite-dimensional system of algebraic equations.
- the cells can be of any geometric shape, such as parallelepipeds (or cubes) or hexahedra (typically having four approximately vertical corner edges which may vary in length), or tetrahedra, rhomboids, trapezoids, or triangles.
- the grid can comprise rectangular cells organized in a regular, structured pattern. Structured cells therefore have similar shape and the same number of sides or faces.
- the most commonly used structured grids are Cartesian or radial in which each cell has four sides in two dimensions or six faces in three dimensions.
- the grid can also comprise cells having a variety of shapes laid out in an irregular, unstructured pattern, or the grid can comprise a plurality of both structured and unstructured patterns.
- a grid is called flexible or unstructured when it is made up of polygons (polyhedra in three dimensions) having shapes, sizes, and number of sides or faces that can vary from place to place. Completely unstructured grids can be assembled that assume almost any shape. When viewed from above, all the cells are preferably, but not necessarily, boundary aligned with other cells, thereby avoiding having any side of a cell contacting the sides of two other cells.
- a mathematical construct called a connection may be defined to characterize the flow of a property between each pair of neighboring cells. If a property, for example fluid mass, can be transported directly between two cells, the cells are considered “connected,” and the fluid is represented as flowing through the connection between them. Alternatively a criterion may be selected to judge if two cells are connected. A cell may participate in from zero to many tens of connections, and most commonly participates in two to eight connections. Typically the total number of connections in a reservoir simulation model is three to four times the total number of cells. Properties are defined for these cell connections, just as properties are defined for the cells themselves.
- the properties of cells and connections may be derived, for example, from a geologic reservoir characterization through a two-step process called gridding and upscaling.
- a geologic reservoir characterization may be, for example, maps which include variations of geologic properties over the reservoir area or a 3D geologic model.
- the geologic reservoir characterization commonly is in the form of a 3D geologic model that discretizes the volume of interest into a large number of cells similar to a reservoir simulation model. Because geologic models are not intended for use in solving transient flow equations, they can accommodate much higher spatial resolution than a reservoir simulation model. 3D geologic models typically have many more cells than a reservoir simulation model (typically hundreds of times as many), and these geologic-model cells are usually much smaller than reservoir- simulation cells.
- the user may lay out the reservoir simulation grid on top of the 3D geologic model.
- the 3D geologic model is in the form of a structured grid while the reservoir simulation grid is an unstructured grid.
- the reservoir simulation grid may have larger cells with irregular shapes.
- the upscaling step undertakes to assign properties to the reservoir simulation cells based on the properties of the geologic-model cells. It may also be used to calculate connection properties for the reservoir simulation grid based on properties of the geologic-model cells.
- connection properties in the following way. First, potential connections could be identified between the geologic model's grid cells. Then a value for a property of interest, for example transmissibility, could be estimated based upon some characteristics of two neighboring cells. For example, transmissibility could be estimated based upon each cell's geometry (face area between neighboring cells and respective distance from each cell center to the face) and the permeability of each cell.
- a property of interest for example transmissibility
- the geologic model may itself contain connectivity anomalies that escape detection using current technology.
- Embodiments of the invention may be applied to 3D geologic models as well as reservoir simulation models. Explicitly controlling connectivity in the geologic model would result in higher-quality reservoir simulation models.
- connectivity anomalies are revealed by flow simulation using reservoir simulation models. After such anomalies are revealed the geologic model is rebuilt in order to cure the anomalies. Time would be saved if the anomalies were identified in the geologic model using embodiments of the invention.
- connection properties for neighboring cells as discussed above and then implement the following steps to identify anomalies.
- After identifying the source of the anomalies the geologic model could be corrected by exerting control over these shapes and placements, perhaps by adjusting the stratigraphic framework. This same procedure could also be used to identify and correct anomalies in a reservoir simulation model.
- a system of cells and connections may be translated into a graph-theory data structure 10 for use in some embodiments of the invention.
- Any 3D cellular model can be characterized as a collection of nodes (called vertices 11 in graph theory) linked where appropriate by connections (called edges 12 in graph theory). Vertices 11 correspond to cells, and edges 12, which correspond to connections, exist wherever two cells are connected, for example where faces of neighboring cells overlap.
- Edges may have properties that include, for example, transmissibility, which is a measure of the connection's conductance to fluid flow under the influence of pressure difference between the two connected vertices.
- This graph-theory data structure is native to modern reservoir simulation applications, but 3D geologic models may require significant computation to identify connections and to calculate transmissibilities as discussed previously herein.
- the graph-theory data structure may be extended from the reservoir grid to include well completions or other oil and gas production or processing equipment. This extension allows well vertices to be on paths. In actual applications of embodiments of the invention it is common for well vertices to be designated as a source vertex or as a destination vertex.
- the user may decide whether or not to permit a well vertex to be on the interior of a path, which is to say that fluid enters the well completion from one reservoir vertex and exits the well completion to a different reservoir vertex (crossflow).
- the well vertex may serve as the source or destination of a shortest path, but may not appear on the interior of any shortest path.
- the extension to include well vertices requires the host application to be aware of the time-dependent well-completion state, because the properties of connections between a well-completion and the reservoir depend on time.
- the graph-theory data structure may also be reduced by applying connection criteria to the connections to thereby reduce the number of edges.
- connections may be formed for any cells that have any overlap with neighboring cells or between cells that other wise have the potential for fluid flow.
- the resulting connections can be further evaluated in the process of translating connections into edges. For example, each connection may be tested against certain criteria. Connections that fail the criteria are not translated into edges while connections that pass the criteria are translated into edges.
- a non-limiting example of a connection-edge criterion includes the requirement that transmissibility exceeds a specified minimum value.
- connections could be tested by setting criteria for one or both of a connection's cells. If one or both of a connection's cells failed any of the criteria, then the associated connection would not translate into an edge in the graph-theory data set.
- the reservoir simulation model may be calibrated using historical performance data, which most often consist of measurements at regular intervals of produced fluid volumes and periodic measurements of pressures in wells. This calibration process is called history matching. History matching commonly requires revising volume, connectivity, and structure (elevations of key surfaces).
- Intra-compartment connectivity is relatively high, and inter-compartment connectivity is relatively low. Compartments may be delineated by faults (structural disruptions that tear stratigraphic surfaces) or by stratigraphic features such as laterally extensive shales. Modern reservoir-simulation grids commonly have tens to hundreds of faults and several stratigraphic units. The ability to understand and control intra- and inter-compartment connectivity is useful for successful building and history-matching reservoir simulation models. Embodiments of the invention may be used to complete history matching to adjust connection/edge properties, for example to reduce transmissibility of a connection/edge that connects cells/ vertices on different sides of a fault.
- the user can associate a cost or weight with some or all of the edges in the data structure.
- the cost of an edge may be defined by the user.
- resistance to flow is a cost that is well-suited to assessing connectivity and can be calculated using only input properties. Consequently, embodiments of the invention may be used as a tool for finding paths of least resistance in a reservoir.
- shortest path should be understood to mean “minimum-cost path” within this document.
- a path is a sequence of vertices and the path cost is the sum of the costs over the edges connecting successive vertices in the path. Using distance as the cost measure is useful only to illustrate that the path of least resistance is almost always longer (more tortuous) than the shortest path.
- the path cost can be interpreted as the pressure difference needed to induce a flow of a unit- mobility fluid at a volumetric rate of one standard unit.
- this means the path cost is the pressure difference in pounds per square inch needed to induce a flow of one barrel per day of a fluid with viscosity of one centipoise and a relative permeability of one. Therefore the path cost may have physical meaning and is not just a mathematical abstraction in embodiments of the invention.
- the best path could mean, for example, the path likely to conduct most of the fluid flow if a potential (pressure) difference were imposed between the two cells. If the two cells are not connected, there is no best path.
- Edge properties may include flow rates and resistance to flow for each fluid phase.
- the flow rate for a specific phase on a edge may be calculated as the product of three terms: transmissibility, phase mobility and phase potential difference.
- Transmissibility is the component of conductance that is independent of the fluid phases and depends only on permeability and grid geometry. The value is roughly the permeability times the face area divided by the distance between cell centers. More precisely, it attempts to account for the permeability field and the grid geometry so as to yield approximately accurate flow rates when multiplied by phase mobility and phase potential difference. Calculating reasonable connection/edge transmissibilities is not a trivial undertaking, but it is useful in reservoir simulation.
- One method of calculating transmissibility uses neighboring cell geometry and permeability as discussed earlier herein.
- connection/edge transmissibility can have a significant effect on the resulting model.
- Transmissibility can vary greatly depending on the methodology used during scale-up because it condenses so much information into so little. Transmissibilities are generally constant over time unless permeability or grid geometry depends on fluid pressure.
- Phase Mobility is the phase-dependent component of conductance. The value is the phase relative permeability divided by phase viscosity. It reflects how readily this phase flows in response to a pressure gradient (low- viscosity fluid flows more readily) and how successfully this phase competes with other phases (higher relative permeability is more successful). Phase mobility is normally evaluated using the properties of the "upstream" cell (the one with higher potential).
- Phase Potential Difference between the participating cells This is the driving force for flow. It depends on pressures, elevation of the cell (gravity), and relative amounts (saturations) of the phases.
- the heat flux, heat capacity or thermal conductivity of the fluid phase and/or the geologic strata could be used as edge properties.
- flux of a constituent, concentration differences, coefficient of diffusion or diffusivity could be used as edge properties.
- electrical conductivity of the fluid phase and or geologic strata could be used as edge properties.
- the invention uses graph-theory single-source, shortest paths algorithms to find the lowest cost path between a source vertex and one or more destination vertices.
- One problem in graph theory is to find the shortest path between two vertices in a graph.
- each edge is given a weight or cost.
- a shortest path is any path whose path cost is equal to the shortest path weight.
- FIG. 1 A graphical depiction of the shortest path solution is depicted in Fig 6.
- a user may wish to find the lowest cost path from first vertex 1 A to second vertex 2 A.
- each edge would have an associated cost value.
- the depicted solution, or path 13 would represent the lowest summed cost across any path from first vertex IA to second vertex IB.
- the solution path 13 travels through vertices IA, IB, 1C, ID, IE, IF to destination vertex 2A.
- the solution path 13 also travels through the edges contained between the above-listed vertices, that is through edges 12A through 12F.
- the total cost of the solution path 13 is found by summing the costs associated with each of the edges (12A-12F) located on the path 13.
- the invention uses a single-source, shortest paths (SSSP) algorithm, for example, Dijkstra's algorithm, which was published in 1959 and is widely used in such familiar applications as route-planning for a road trip. It finds the path from a source node (vertex) to a destination node (vertex) that minimizes the summation over the path (edges) of an edge property (the cost).
- SSSP single-source, shortest paths
- Dijkstra's algorithm A remarkable property of Dijkstra's algorithm is that it finds all shortest (minimum-cost) paths from a particular source vertex to all destination vertices just as quickly as it finds a single shortest path from that source to a particular destination vertex.
- Another example of a SSSP algorithm that may be used in the invention is the Bellman-Ford algorithm.
- the simplest illustrative cost property is the reciprocal of transmissibility, but could be any edge property (for example, the reciprocal of the product of transmissibility and phase mobility, where mobility is relative permeability divided by viscosity).
- the invention may implement Dijkstra's SSSP algorithm using a priority queue that in turn uses a binary heap.
- Dijkstra's algorithm maintains a set S of vertices whose final shortest- path weights from the source s have already been determined. The algorithm repeatedly selects the vertex w e V - S with the minimum shortest-path estimate d, inserts u into S, and relaxes all edges leaving u. This implementation maintains a priority queue that contains all the vertices in V - S, keyed by their d values. This implementation assumes that the graph G is represented by adjacency lists, which facilitate the identification of all neighboring vertices for any specific vertex. [0075] DIJKSTRA(G,W,S)
- Step 8 for each neighboring vertex v of the just-added vertex u, i.e., relax all edges (u, v) of vertex u.
- Relaxation decreases an upper bound d[v] on the actual shortest-path weight of a vertex ⁇ v] until the upper bound equals the shortest-path weight ⁇ v]. Relaxation is the only means by which shortest-path estimates and predecessors change.
- the process of relaxing an edge (u, v) consists of testing whether by going through u we can improve the shortest path to v found so far, and if so, updating v's shortest-path estimate d[v] and predecessor field ⁇ ]. In effect, vertex u is being evaluated as an improved predecessor for its neighboring vertex v.
- Steps 2-4 If the current shortest-path estimate d[v] for the weight of a shortest path from the source s to vertex v is greater than the shortest-path estimate d[u] for vertex u plus the (non-negative) weight w(u, v) for the edge (u, v) that leads from u to v, then implement Steps 2-4. That is, u is an improved predecessor for v because the path to v through u reduces the shortest-path estimate d[v] for v.
- a priority queue is a data structure for maintaining a set Q of elements, each with an associated value called a key.
- a priority queue supports the following operations (among others):
- INSERT(Q, x) inserts the element x into the set Q in a way that preserves sorting by the key value. This operation could be written as Q ⁇ — Q ⁇ x ⁇ .
- EXTRACT-MIN(Q) removes and returns the element of Q with the smallest key.
- DECREASE-KEY(Q, X) updates the priority queue (by preserving the sorting by key value) to account for a change in the key value for element x.
- the elements in the priority queue are vertices v in the set V - S, which is the set of vertices whose shortest paths from the source s have not yet been determined.
- the key is the shortest-path estimate d[v] for the weight of a shortest path from the source s to v.
- a priority queue may be implemented using a binary heap data structure.
- the procedures that update binary heaps preserve sorting by enforcing the "heap property."
- the heap property is that the key-field value of each of at most two child nodes is greater than or equal to the key-field value of a parent node.
- the priority-queue operations may be implemented using corresponding binary-heap procedure:
- INSERT (Q, X) may be implemented by HEAP-INSERT(A, k)
- EXTRACT-MIN(Q) may be implemented by HEAP-EXTRACT-MIN(A)
- DECREASE-KEY(Q, X) may be implemented by HEAP-DECREASE-KEY(A, i, k)
- A is an array that represents a heap and k is the key value for element x, and / is the index of a node.
- Both HEAP-INSERT(A, k) and HEAP-DECREASE-KEY(A, i, k) use iteration to enforce the heap property while updating the heap.
- HEAP-EXTRACT-MIN(A) invokes the heap operation HEAPIFY(A, i), which is usually implemented using recursion.
- An alternate implementation of INITIALIZE-SINGLE-SOURCE(G, S) may also use the recursive HEAPIFY(A, i) operation rather than the iterative Heap-Insert(A, k).
- Embodiments of the invention may substitute an iterative implementation of HEAPIFY(A, i) for the conventional recursive implementation to avoid the memory exhaustion and reduced speed associated with the recursive implementation.
- Dijkstra's algorithm solves the single-source shortest- paths problem on a weighted (i.e. containing costs) graph (directed or undirected), provided all edge-weights are nonnegative.
- Dijkstra's algorithm maintains a set of vertices whose final shortest-path weights from the source have already been calculated, along with a complementary set of vertices whose shortest-path weights have not yet been determined.
- the algorithm repeatedly selects the vertex with the minimum current shortest-path estimate among those whose shortest-path weights have not yet been determined. It updates the weight estimates to all vertices adjacent to the currently selected vertex (this updating is commonly referred to as "relaxation" of the edges between these vertices).
- the vertex is then added to the set of vertices whose final shortest-path weights have been calculated. It continues to do this until all vertices' final shortest-path weights have been calculated (until the set of vertices whose shortest-path weights have not been determined is empty).
- the SSSP algorithm may yield the following results for each vertex in the graph.
- the SSSP algorithm may yield the path cost from the source vertex to the current vertex. If the current vertex is not connected to the source vertex, a special value may be assigned to path cost, either negative or large and positive.
- the SSSP algorithm may yield the identity of the predecessor vertex on the shortest path.
- the source vertex has no predecessor, so its predecessor data item may be assigned a special value.
- Embodiments of the invention capture the identity of the edge between the predecessor vertex and the current vertex. This facilitates building an edge- view of shortest paths (as opposed to the conventional vertex-view).
- the vertex-view and edge-view are foundations for separate (vertex and edge) styles of 3D rendering and analysis.
- Embodiments of the invention may alternatively include a 3D vertex- viewer included in the cellular-modeling computing application.
- 3D viewers include any of the body of software methodologies for displaying a three-dimensional representation of a geologic reservoir.
- 3D viewers typically display the geologic reservoir in a multitude of colors and can be manipulated through a user interface to alter the content and quality of the graphical representation.
- the 3D viewer may be used to display two-dimensional renderings of the three-dimensional collections of cells that make up the cellular model.
- Each cell's color may be determined by the value of a user-selected cell property and by user-selected controls on the correspondence between color and value (the color legend).
- the 3D viewer may be used as the home for the dialog (user interface) through which the user controls and applies embodiments of the invention.
- a 3D-viewer-based dialog may support embodiments of the invention in various ways, for example, by setting controls for shortest-path computations, through input processing, by selecting computational tasks and by viewing results.
- a 3D-viewer-based dialog may be used in setting controls for shortest- path computations.
- the 3D-viewer-based dialog may be used in defining the cost to be minimized for the shortest paths computation. This cost may be evaluated for each edge and accumulated over the sequence of edges associated with the path. For example, the edge cost may be the reciprocal of transmissibility.
- Other illustrative examples include the connection length (distance) and the reciprocal of [transmissibility multiplied by phase mobility].
- a 3D-viewer-based dialog may be used in defining vertex and edge criteria. Vertices and edges that fail the criteria would be excluded from participating in any path. An example is, edge transmissibility must be at least 0.001 (RB-cp)/(d- psi), otherwise, the edge is ignored.
- a 3D-viewer-based dialog may be used in determining whether or not to extend the graph from reservoir vertices only by adding well vertices.
- a 3D-viewer-based dialog may be used in determining whether or not well vertices may be on the interior of a path or are restricted to participating only as source or destination vertices, i.e., whether or not to permit crossflow. This setting should be controllable for each individual well vertex.
- a 3D-viewer-based dialog may be used to accomplish input processing.
- Source vertices may be, for example, either reservoir vertices or well vertices.
- Reservoir vertices may be added to the list by exposing a cell (vertex) to view, selecting the cell with the mouse, then clicking.
- Well vertices may be added to the list by selecting them by name from a list.
- Some uses of the invention may use a list of one or more destination vertices. Composing this list of destination vertices may be accomplished using the 3D viewer in exactly the same manner as for composing the list of source vertices.
- a 3D-viewer-based dialog may be used in selecting the computational task to be preformed.
- Illustrative examples of such tasks include the following. Computing all shortest paths for each member in the current list of source vertices. Every vertex is a potential destination vertex. Using embodiments of the invention, the SSSP computations take about one second per source vertex. Another example includes identifying and displaying specific shortest paths for the current list of destination vertices. This process is effectively instantaneous. Another example includes conducting a multi-source competition. This means that each vertex is assigned an owner vertex from the current list of source vertices. The owner vertex is the source vertex that minimizes cost compared with other source vertices in the list. In effect, this process identifies for every potential destination vertex the source vertex to which the destination vertex is best connected (i.e., connected with minimum cost). This assignment takes about one second.
- Additional examples include computing user-defined cumulative measures as a function of cost for each of the source vertices and displaying the results graphically.
- An example measure is the volume of oil. The relative merits of drilling locations could be assessed by the accessibility of the location (the source vertex) to oil. Better locations have greater amounts of oil at lower cost (less resistance to flow).
- Another example includes managing the special task of determining shortest paths between all possible pairs of wells. Example 2 discusses this in detail.
- Another example includes exporting to a file results from any of the above computational tasks. The exported data may be supplemented by user-selected vertex properties and edge properties. Captured data are the source-vertex identifier, the current-vertex identifier, the predecessor-vertex identifier, the predecessor edge identifier, the path cost, and user-selected properties for the current vertex and for the predecessor edge.
- a 3D-viewer-based dialog may be used in viewing results.
- the computational tasks yield vertex properties that allow the user to view the results of the invention using the standard features of the 3D viewer.
- Paths are indicated by displaying the cells that correspond to vertices on the path rather than by rendering a 3D polyline whose segments connect cell centers and face centers along the path. The user chooses which property to display.
- the viewing of results may include viewing all shortest paths for a selected source vertex.
- the user will display the cost property, whose value is the cost of the shortest path from the selected source vertex.
- the source vertex has a cost of zero. Vertices that have no path connecting them to the source vertex have a special value for cost, either a negative value or a large, positive value.
- the threshold tool may be used to display only those vertices that satisfy user-specified criteria evaluated using the vertex's properties.
- the user specifies a list of values.
- the user defines a range of values by specifying minimum and maximum values.
- the user When viewing the all-shortest-paths results for a selected source vertex, the user thresholds on the cost property, specifying a range from zero to some small, positive value. The viewer displays all paths with costs in this range. By nudging the upper end of the cost range, the user can develop a feel for the way in which the source vertex is connected to its neighborhood. In many cases, the all-shortest-paths view is overwhelming, and it is necessary to view specific shortest paths.
- Another example of viewing results includes viewing specific shortest paths from a selected source vertex (or competition among a set of source vertices) to a list of selected destination vertices.
- a specific path consists of the sequence of vertices that starts at the source vertex and ends with the destination vertex, with ascending cost en route.
- the 3D viewer is capable of instantaneously displaying specific shortest paths to selected destination vertices. The viewer composes this display by identifying the vertices on the shortest path(s) and displaying only these vertices (in effect, imposing a hidden threshold condition).
- Fig 7 is a graphical depiction of a lowest cost connected path 23 from a source vertex 20 to a destination vertex 21.
- Fig 7 also depicts the minimum distance path 22 from source vertex 20 to destination vertex 21 as a contrast to show that the lowest cost path 23 is often less direct (more tortuous) than the shortest minimum distance path 22.
- Fig 8 depicts a more extreme example of the lowest cost path 23B being more tortuous than the minimum distance path 22B. Cost and distance values for the four minimum paths are given in Table 1.
- Another example (see Example 5) of viewing results includes viewing a collection of specific shortest paths between multiple sources and multiple destinations. Such displays suffer a loss of information (path cost and identity) whenever paths from different sources overlay, because the 3D viewer only accommodates a single value per vertex for each property. The user may accept this loss of information, or may select a feature that creates multiple views, each of which is guaranteed to have no overlapping paths. However, overlapping paths do not suffer a loss of information if the paths have the same source.
- Another example of viewing results includes viewing the results of a multi-source competition. A multi-source competition adds the notion of owner source- vertex.
- the user may use the 3D viewer's threshold tool to display a subset of owner source- vertices, as few as one owner at a time, or as many as all owners at once. For a particular selection of owner vertices, the user may choose to view all shortest paths (typically thresholding on maximum cost) or to view specific shortest paths to selected destinations.
- FIG. 1 Another example of viewing results includes viewing a line-graph (i.e., not a 3D rendering) of the cumulative distribution function (CDF) of all path costs for selected sources. This feature is described in Example 4.
- CDF cumulative distribution function
- Embodiments of the invention may be used to determine along which path would fluid flow if a pressure difference were imposed between vertex A and vertex B and if the fluid were required to choose a single, best path.
- Conventional solutions to this problem would define wells associated with vertices A and B, then solve flow equations to determine the velocity field. Defining wells is cumbersome, solving flow equations is slow, and interpreting the results is challenging.
- Embodiments of the invention may use existing cellular models that are readily cast as graphs (nodes (vertices) and connections (edges)), and graph-theory algorithms that efficiently solve the shortest-path problem.
- embodiments of the invention require only a point-and-click of the mouse (a handheld pointing device used to point to spots on the computer monitor) to mark vertices A and B.
- embodiments of the invention require minutes to tens of hours to solve flow equations, embodiments of the invention require about one second to find all shortest paths.
- embodiments of the invention show paths by displaying the sequence of cells that make up the path.
- Embodiments of the invention may increase the flexibility and speed of designating a source vertex and destination vertex. Any cell (or well) can play the role of source or destination and can be selected with a mouse click in the 3D viewer. No additional input is needed, in contrast with the complex input needed to describe a well for conventional flow simulation.
- Embodiments of the invention may display results very quickly.
- the single-source, shortest paths (SSSP) algorithm identifies all shortest paths from one source to all destinations in about one second. Displaying shortest paths from a source to any set of destinations may be instantaneous.
- Embodiments of the invention may be used to display visual results that are easily understood. Individual paths may be displayed as an unequivocal sequence of connected cells.
- One way to find a leak is to display the path of least resistance between a vertex in one compartment and a vertex in the other compartment.
- This task could be called, "single source, single path.” It is accomplished as follows: 1. In the 3D viewer, launch the shortest-path dialog, using the default definition of path cost (the reciprocal of edge transmissibility). 2. Pick a single source vertex either as a cell in Compartment 1 (pick by pointing at the cell with the mouse, then clicking) or as a well (-vertex) in Compartment 1 (point-and-click on the well, or select it from a list). A well vertex typically connects to multiple reservoir vertices, so the user must be confident that all such reservoir vertices are in Compartment 1. 3.
- the resulting 3D viewer display shows the path of least resistance between the two selected vertices.
- the shortest-path dialog may also provide the capability of managing (viewing and exporting to file) the shortest-path data, expressed either as a sequence of vertices (and user-selected vertex properties) or as a sequence of edges (and user-selected edge properties).
- the user could inspect the edge data for the purpose of choosing one edge as a target for transmissibility editing. There should be a geologic basis for picking this target edge, such as association with a fault or stratigraphic sequence boundary.
- the editing action could include setting transmissibility to a tiny value (even zero) on the target edge.
- One choice for the target edge is the one that connects a vertex in Compartment 1 to a vertex in Compartment 2.
- Figure 9 depicts three segmented compartments 40, 41 and 43, which are separated by faults 46 and 47 on a 3-D viewer.
- a source vertex 42 in first compartment 40 is selected to serve as a source vertex in one method of leak detection and the SSSP algorithm is run for this source vertex 42.
- Figures 10, 11 and 12 depict 3-D viewer representations of all the vertices that are connected to the source vertex 42 at three successively higher threshold values for cost. In these figures cost is represented as the summed cost of the reciprocal of transmissibility. At a threshold value of 0.005 only four vertices are displayed in Figure 10 and all such vertices are within the first compartment 40.
- Figure 12 depicts the two compartment system at a threshold value of 0.40 and shows the magnitude of the model anomaly.
- Another way to determine where best to impair or sever communication on a leak path is to display additional paths and then inspect the set of paths for commonality or overlap. This can be done by adding to the list of destination vertices in Compartment 2, yielding a display with multiple paths from the single source vertex in Compartment 1 to the selected destination vertices in Compartment 2. If there is a single (or “best") leak between the two compartments, the paths will overlay exactly in the source compartment (Compartment 1), and will diverge within the destination compartment (Compartment 2).
- 3D viewer's "Edit Property” tool to seal the leak by reducing the transmissibility of a selected edge on the leak path. This cycle of find-and-seal may have to be repeated (possibly many times, or perhaps not at all) to ensure the isolation of the two compartments, because it is unlikely that the original leakage resulted from a single edge.
- Embodiments of the invention make it possible to identify and display the path of least resistance between every pair of wells.
- the path cost is a natural measure of coupling, i.e., the lower the cost, the higher the coupling.
- the invention may be used to reveal in a quantitative and visual way the connectivity features that would otherwise be hidden.
- the distance between wells is often not a good measure of their coupling. That is, a well may be better connected to a remote well than to a closer well. This knowledge is useful for picking drilling locations and for adjusting grid properties to calibrate reservoir performance with historical data (history matching).
- the number of well pairs is N*(N-l)/2.
- the corresponding shortest paths could be determined in N-l SSSP calculations, each taking about one second. Since most reservoir models have a few tens of wells, the time required to determine shortest paths for all possible well pairs is not prohibitive.
- Well state affects any shortest-path assessment because properties of edges between well-vertices and reservoir vertices are time-dependent, reflecting real conditions in the oilfield (i.e., reflecting the life cycle of real wells: drilling, completion, stimulations, workovers, and abandonment). More accurate assessment of well coupling may be accomplished where the SSSP calculations be repeated every time any well-state changes.
- the user picks a particular time (date), the reservoir model knows all well-states at that time, and the invention evaluates shortest paths taking into account the current well-states.
- the number of paths displayed may be reduced so as to minimize visual overload.
- the number could be reduced by displaying only M nearest (measured by path-cost) neighbors, where M may be three or four.
- Another approach is to display only those well-pair paths with (total) cost below some threshold.
- any path-related property may lose information and be confusing. That is because where two or more paths overlay (i.e. share one or more vertices), the viewer can only display property- values for one of those paths. In practice, it is common that paths overlay, because connectivity tends to have an arterial structure. No such loss of information occurs when there is a single source and paths to multiple destinations overlap or when there are multiple sources and the user is displaying the results of a multi-source competition. Where the vertices and edges overlap, the path properties are identical, irrespective of the destination. Where there are multiple sources and the user is displaying the results of a multi-source competition, no two paths with different sources ever overlap (although two paths with the same source may overlap).
- One use of the invention for understanding well-coupling is to capture the all-well-pairs shortest paths results in a file for analysis, and to display well-pair shortest paths one source-well at a time. This display strategy avoids the issues of too many paths and multi-source path overlap.
- the user first launches the all-well-pairs shortest-paths task, then saves the file of results to serve as a guide as to which source-wells to visit for display purposes (e.g., those wells whose lowest-cost pairings have either very low cost or very high cost).
- the user picks a well of interest to serve as the source vertex, then launches the SSSP calculation for that source. Then the user requests to display specific shortest paths by multi-selecting some or all other well-vertices as destinations, guided by the all-well-pairs results This display strategy avoids the ambiguity associated with displaying overlapping paths from multiple sources.
- a multi-source competition in which the sources are all well-vertices assigns to each vertex an owner well-vertex.
- the multi-source competition in effect claims that all well- vertices act as producing wells with identical bottomhole pressures.
- the competition determines which well-vertices drain which reservoir vertices by assessing which well-vertex has the lowest-cost path to each reservoir vertex. Like any other task involving well vertices, this task may consider the time- dependent well state.
- a multi-source competition may use reservoir vertices as sources rather than well vertices.
- Different combinations of the settings expose different features of the grid's connectivity with respect to the sources (wells).
- Illustrative examples include the list of sources for which drainage volumes are to be displayed. This may be one, all, or a specific subset of the sources.
- the 3D viewer lets the user control time, because well states depend on time. At the start of the model, there may be only a few wells, and these have large drainage volumes. As time advances and additional wells are drilled and completed, these new wells claim drainage volume from older wells. When a well is out of service, its drainage volume will likely be claimed by its active neighbors.
- the 3D viewer lets the user control threshold (maximum value) on path cost (applies to all sources). If the user is displaying multiple sources, setting a low maximum-cost threshold yields some separation between adjacent drainage volumes. This makes it easier to understand the picture. If the user is displaying a single source, the maximum-cost can be nudged up from low values to high values to gain an appreciation for the shape of the drainage volume and to reveal obvious asymmetries.
- the maximum-cost can be nudged up from low values to high values to gain an appreciation for the shape of the drainage volume and to reveal obvious asymmetries.
- vertices In general there will be some vertices that are not connected to any of the selected, multiple sources. Such vertices may be identified as being owned by the source named "None,” and if present in sufficient number and if they contain a sufficient amount of hydrocarbon, may represent a target for drilling. Commonly these vertices are few and poorly connected among themselves, because most current reservoir simulation grids have relatively high connectivity.
- the drainage-volume view provided by the multi-source competition feature of embodiments of the invention yields a fast and information-rich view of the connectivity features of a reservoir simulation grid.
- the distinction between good wells and poor wells is readily apparent by comparing the size of their drainage volumes.
- the benefits of horizontal wells are evident from the size and shape of their drainage volumes. Oil that is difficult to drain from existing wells, hence represents a potential target for drilling and production, can be visualized by using the 3D viewer's threshold tool to display all vertices with cost higher than a user-specified value with all source vertices displayed.
- the multi-source competition results can be used analytically in at least the following ways.
- the multi-source competition results can be used in defining regions for averaging reservoir properties. Several situations arise in reservoir simulation in which it is necessary to identify the neighborhood (region) associated with each well. The region is needed to calculate average reservoir pressure, which may be used in history-matching to compare with field measurements of pressure in the well or used in prediction to control production and injection so as to conform to pressure guidelines. Current practice is to define well regions in an arbitrary (connectivity-blind) way, which is often at odds with the connectivity features revealed by the invention.
- the 3D-viewer interface could include a "Define Regions" dialog that would let the user mark sets of vertices as belonging to named regions based on those vertices in the view composed using a threshold on path-cost. These connectivity-based region definitions should be superior to regions that are defined without considering connectivity.
- the multi-source competition results can be used in calculating hydrocarbon volume versus cost. The results of a multi-source competition can be analyzed to yield for each source a curve of cumulative amount of hydrocarbon (oil and gas) within a given cost (resistance to flow).
- This curve is built for one source vertex by sorting the vertices in its drainage volume from lowest to highest cost, determining the hydrocarbon contents of each vertex, and accumulating the contents over this sorted list.
- This curve conveys the accessibility of this source vertex (well) to hydrocarbon. A more favorable well has more hydrocarbons at lower cost than a less favorable well.
- volume-versus-cost curves derived from multi-source competition offer a fast, quantitative way to assess drilling locations because the curves measure both contents and connectivity, and because they show to what extent a new well would harm (take production from) nearby existing wells.
- volume-versus-cost curves can also be used in history matching.
- transmissibilities should also be adjusted so as to preserve the original character evidenced in the volume-versus-cost curves. For example, if pore volume is significantly increased near a well, transmissibilities should be reduced. Otherwise, the volume-versus-cost curve may show an unnatural and exaggerated accessibility of the well to the added pore volume.
- launching the SSSP algorithm for a single source vertex takes about one second to identify the shortest paths from that source vertex to all other vertices in the model.
- the simplest way to exploit this information is to pick one or more destination vertices and display the shortest paths from the single source vertex to the set of selected destination vertices. In this case (meaning, for a single source vertex), the multiple paths commonly overlay near the source and diverge to seek out the destinations. There is no loss of information if the paths overlay near the source.
- Another approach to interpreting the connectivity with respect to a single source is to display all paths at once, using the 3D viewer's property- threshold tool to scan from destinations close to the source (low cost) toward destinations far from the source (high cost).
- the path-cost property is targeted for thresholding, the threshold minimum cost is set to zero (the source-vertex cost), and the maximum cost is set to a tiny value. Then the maximum value is nudged higher in small steps, at each step observing the shape of the growing collection of vertices.
- Figs 13-17 the source vertex 30 is depicted with a producing well 31 and the succeeding figures are completed at successively increasing values of summed reciprocal of transmissibility as the threshold value.
- the threshold value is 0.05.
- the threshold value is increased in the succeeding figures as follows: Fig 14, 0.20; Fig 15, 0.60; Fig 16, 2.60.
- Figure 14 the vertices that meet the threshold value form a nearly uniform circular lens 35.
- a second lens 36 begins to appear which indicates a path or chimney from the first lens 35 to the second lens 36.
- various lenses have been agglomerated to form an agglomerated mass of vertices 37.
- Figure 17 was generated by selecting the shortest path from source vertex 30 contained in first lens 35 to destination vertex 38 contained in the lowermost lens of the agglomerated mass 37 and connected to producing well 33 (32 and 34 are other producing wells).
- the shortest path uncovers the highly conductive chimney 39 that serves as the path by which the respective lenses are connected. This chimney 39 is most likely an anomaly that can be corrected in the model.
- the display could be interpreted as a "Medusa's head" of all shortest paths whose path-costs are confined within the current threshold.
- the advantage of seeing all paths is that shapes of the all-paths view often reveal important connectivity features.
- the all-paths view would present as a growing sphere (or cylinder of circular cross-section when confined between stratigraphic surfaces). But real reservoir grids generally present strikingly asymmetric (and informative) all-paths views. For example, if the source is on the edge of a sand body, the all-paths view would appear confined by the nearby bounding surfaces of the sand body but unconfined where the boundaries were distant.
- the view is so complex as to defy interpretation. Such cases may be better understood by displaying paths for a set of specific destinations or by using standard tools of the 3D viewer to slice the all-paths view.
- the implementation of the invention in the 3D viewer may provide a statistical-graph capability (the cumulative distribution function [CDF] of path-cost) that accommodates one or more single-source results.
- the CDF sorts the target property (path cost) from low to high values, and shows what proportion of all paths have a path-cost less than any particular value. Consequently, the CDF of path cost clearly reveals sharp increases in path-cost as steps in the otherwise smooth curve.
- a useful task to automate allows the user to generate and display in the
- N may be any integer, but would normally be from ten to 100.
- the user interface may be implemented in the 3D viewer, and it requests two input items from the user:
- the N random, shortest paths are identified by 2N random draws without replacement from the set of vertices, which yields 2N unique vertices, i.e., no vertex appears twice in this list of 2N vertices.
- the first vertex is treated as a source vertex, and an SSSP calculation is launched for that source vertex.
- the second vertex is treated as a destination vertex for that source vertex, and the shortest path is identified for this source-destination pair. This process is repeated (N-l) times, consuming the 2N randomly drawn vertices and yielding N random, independent, shortest paths.
- the resulting table of N source-destination vertex pairs may be viewed, and allows source (or destination) vertices of particular interest to be copied for re-use as source (or destination) vertices in other shortest-paths tools, such as displaying all shortest paths or multiple, specific paths for a source vertex.
- these shortest paths will also have an integer identifier, which allows subsets of them to be displayed (and the others hidden) using the 3D viewer's threshold tool.
- the user can choose to segregate the paths into subsets such that within each subset, no paths intersect.
- This segregation is implemented by an algorithm that identifies intersections between pairs of paths.
- Each candidate path visits each subset of non-intersecting paths and is tested for intersection with each resident path.
- the candidate path Upon encountering the first intersection, the candidate path is said to fail this subset.
- the candidate path then exits the current subset and moves to the next, where the process is repeated. If the candidate path survives (passes) the intersection test with all paths resident in a subset, the candidate path is added to that subset. If the candidate path fails all current subsets, a new subset is created containing only the candidate path.
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
Description
Claims
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/568,940 US7584086B2 (en) | 2003-09-30 | 2004-09-29 | Characterizing connectivity in reservoir models using paths of least resistance |
EP04789530A EP1668561A2 (en) | 2003-09-30 | 2004-09-29 | Characterizing connectivity in reservoir models using paths of least resistance |
CA2538577A CA2538577C (en) | 2003-09-30 | 2004-09-29 | Characterizing connectivity in reservoir models using paths of least resistance |
NO20061871A NO20061871L (en) | 2003-09-30 | 2006-04-27 | Characterization of connectivity in reservoir models using paths with least resistance |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US50733203P | 2003-09-30 | 2003-09-30 | |
US60/507,332 | 2003-09-30 |
Publications (2)
Publication Number | Publication Date |
---|---|
WO2005033739A2 true WO2005033739A2 (en) | 2005-04-14 |
WO2005033739A3 WO2005033739A3 (en) | 2006-04-06 |
Family
ID=34421609
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2004/032942 WO2005033739A2 (en) | 2003-09-30 | 2004-09-29 | Characterizing connectivity in reservoir models using paths of least resistance |
Country Status (6)
Country | Link |
---|---|
US (1) | US7584086B2 (en) |
EP (1) | EP1668561A2 (en) |
CN (1) | CN100590637C (en) |
CA (1) | CA2538577C (en) |
NO (1) | NO20061871L (en) |
WO (1) | WO2005033739A2 (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008021743A2 (en) * | 2006-08-09 | 2008-02-21 | Services Petroliers Schlumberger | Facilitating oilfield development with downhole fluid analysis |
EP1994488A2 (en) * | 2006-03-02 | 2008-11-26 | ExxonMobil Upstream Research Company | Method for quantifying reservoir connectivity using fluid travel times |
EP2098887A1 (en) | 2008-03-03 | 2009-09-09 | Exxonmobil Upstream Research Company | Determining connectivity architecture in 2-D and 3-D heterogeneous data |
US8731873B2 (en) | 2010-04-26 | 2014-05-20 | Exxonmobil Upstream Research Company | System and method for providing data corresponding to physical objects |
WO2015137916A1 (en) * | 2014-03-10 | 2015-09-17 | Landmark Graphics Corporation | Modeling geologic surfaces using unilateral non-node constraints from neighboring surfaces in the stratigraphic sequence |
NO337774B1 (en) * | 2009-05-07 | 2016-06-20 | Saudi Arabian Oil Co | Well pressure system and method |
EP2232406A4 (en) * | 2007-12-21 | 2017-10-18 | Exxonmobil Upstream Research Company | Method and apparatus for analyzing three-dimensional data |
Families Citing this family (71)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100489558C (en) * | 2004-06-07 | 2009-05-20 | 埃克森美孚上游研究公司 | Method for solving implicit reservoir simulation matrix equation |
US8139507B2 (en) * | 2004-09-15 | 2012-03-20 | Hewlett-Packard Development Company, L.P. | Network graph for alternate routes |
FR2886742B1 (en) * | 2005-06-02 | 2007-07-27 | Inst Francais Du Petrole | METHOD OF CHANGING SCALE OF ABSOLUTE PERMEABILITIES TO BUILD A FLOW SIMULATION MODEL |
US20070016389A1 (en) * | 2005-06-24 | 2007-01-18 | Cetin Ozgen | Method and system for accelerating and improving the history matching of a reservoir simulation model |
US7660711B2 (en) * | 2006-04-28 | 2010-02-09 | Saudi Arabian Oil Company | Automated event monitoring system for online reservoir simulation |
FR2918179B1 (en) * | 2007-06-29 | 2009-10-09 | Inst Francais Du Petrole | METHOD FOR ESTIMATING THE PERMEABILITY OF A FRACTURE NETWORK FROM A CONNECTIVITY ANALYSIS |
CN101896690B (en) * | 2007-12-13 | 2015-02-18 | 埃克森美孚上游研究公司 | Parallel adaptive data partitioning on a reservoir simulation using an unstructured grid |
WO2009076728A1 (en) * | 2007-12-17 | 2009-06-25 | Leximancer Pty Ltd | Methods for determining a path through concept nodes |
WO2009079123A1 (en) * | 2007-12-18 | 2009-06-25 | Exxonmobil Upstream Research Company | Determining connectivity architecture in 2-d and 3-d heterogeneous data |
US20100223039A1 (en) * | 2007-12-21 | 2010-09-02 | Serguei Maliassov | Modeling In Sedimentary Basins |
US20140313195A1 (en) * | 2008-02-29 | 2014-10-23 | Cherif Atia Algreatly | 3D Model Mapping |
CA2710809C (en) * | 2008-03-10 | 2017-06-13 | Exxonmobil Upstream Research Company | Method for determining distinct alternative paths between two object sets in 2-d and 3-d heterogeneous data |
CA2717514C (en) | 2008-05-05 | 2016-07-26 | Exxonmobil Upstream Research Company | Systems and methods for connectivity analysis using functional objects |
US8095349B2 (en) * | 2008-05-30 | 2012-01-10 | Kelkar And Associates, Inc. | Dynamic updating of simulation models |
US8355898B2 (en) * | 2008-06-17 | 2013-01-15 | Chevron U.S.A. Inc. | System and method for modeling flow events responsible for the formation of a geological reservoir |
US8560969B2 (en) * | 2008-06-26 | 2013-10-15 | Landmark Graphics Corporation | Systems and methods for imaging operations data in a three-dimensional image |
EP2342668B1 (en) * | 2008-09-19 | 2019-03-13 | Chevron U.S.A., Inc. | Computer-implemented systems and methods for use in modeling a geomechanical reservoir system |
BRPI0919456A2 (en) * | 2008-09-30 | 2015-12-22 | Exxonmobil Upstream Res Co | method for modeling fluid flow in a hydrocarbon reservoir |
US8311339B2 (en) * | 2008-11-05 | 2012-11-13 | Palo Alto Research Center Incorporated | Method to robustly match images with similar body and head pose from image sequence |
AU2009311619B2 (en) * | 2008-11-06 | 2015-10-01 | Exxonmobil Upstream Research Company | System and method for planning a drilling operation |
US8301426B2 (en) * | 2008-11-17 | 2012-10-30 | Landmark Graphics Corporation | Systems and methods for dynamically developing wellbore plans with a reservoir simulator |
US8350851B2 (en) * | 2009-03-05 | 2013-01-08 | Schlumberger Technology Corporation | Right sizing reservoir models |
US8339396B2 (en) * | 2009-03-05 | 2012-12-25 | Schlumberger Technology Corporation | Coarsening and splitting techniques |
CN102667804A (en) * | 2009-11-23 | 2012-09-12 | 埃克森美孚上游研究公司 | Method and system for modeling geologic properties using homogenized mixed finite elements |
FR2953039B1 (en) * | 2009-11-26 | 2012-01-13 | Inst Francais Du Petrole | METHOD OF OPERATING A PETROLEUM FACILITY BY RECONSTRUCTION OF TANK MODEL |
US8931580B2 (en) | 2010-02-03 | 2015-01-13 | Exxonmobil Upstream Research Company | Method for using dynamic target region for well path/drill center optimization |
EP2534605B1 (en) * | 2010-02-12 | 2020-06-17 | Exxonmobil Upstream Research Company | Method and system for partitioning parallel simulation models |
US9703006B2 (en) * | 2010-02-12 | 2017-07-11 | Exxonmobil Upstream Research Company | Method and system for creating history matched simulation models |
US8532968B2 (en) * | 2010-06-16 | 2013-09-10 | Foroil | Method of improving the production of a mature gas or oil field |
US20110310101A1 (en) * | 2010-06-22 | 2011-12-22 | Schlumberger Technology Corporation | Pillar grid conversion |
CN103210397B (en) | 2010-08-09 | 2015-12-09 | 科诺科菲利浦公司 | There is the oil reservoir method of roughening keeping transport |
US9229129B2 (en) * | 2010-12-10 | 2016-01-05 | Conocophillips Company | Reservoir geobody calculation |
EP2678802A4 (en) * | 2011-02-21 | 2017-12-13 | Exxonmobil Upstream Research Company | Reservoir connectivity analysis in a 3d earth model |
BR112013025220A2 (en) | 2011-04-12 | 2016-12-27 | Landmark Graphics Corp | porous media flow simulation method, computer program stored on a non-transient tangible computer readable storage medium |
EP2712440A4 (en) * | 2011-05-17 | 2016-05-25 | Exxonmobil Upstream Res Co | Method for partitioning parallel reservoir simulations in the presence of wells |
US10174593B2 (en) * | 2011-09-20 | 2019-01-08 | Landmark Graphics Corporation | System and method for coarsening in reservoir simulation system |
US9316761B2 (en) | 2012-01-25 | 2016-04-19 | Baker Hughes Incorporated | Determining reservoir connectivity using fluid contact gravity measurements |
US9535978B2 (en) | 2012-03-29 | 2017-01-03 | International Business Machines Corporation | Semantic mapping of topic map meta-models identifying assets and events to include weights |
US9158871B2 (en) | 2012-04-06 | 2015-10-13 | Ut-Battelle, Llc | Graph modeling systems and methods |
US9595129B2 (en) | 2012-05-08 | 2017-03-14 | Exxonmobil Upstream Research Company | Canvas control for 3D data volume processing |
AU2013350307A1 (en) * | 2012-11-20 | 2015-07-09 | Stochastic Simulation Limited | Method and system for characterising subsurface reservoirs |
AU2012396846B2 (en) | 2012-12-13 | 2016-12-22 | Landmark Graphics Corporation | System, method and computer program product for evaluating and ranking geobodies using a Euler Characteristic |
US10481297B2 (en) | 2013-01-28 | 2019-11-19 | Westerngeco L.L.C. | Fluid migration pathway determination |
WO2014117115A1 (en) * | 2013-01-28 | 2014-07-31 | Westerngeco Seismic Holdings Limited | Salt body extraction |
AU2013393346B2 (en) * | 2013-07-02 | 2017-06-29 | Landmark Graphics Corporation | Lofting algorithm for discrete network meshing |
SG11201510115UA (en) | 2013-07-02 | 2016-01-28 | Landmark Graphics Corp | 3d stadia algorithm for discrete network meshing |
WO2015035105A1 (en) * | 2013-09-05 | 2015-03-12 | Schlumberger Canada Limited | Integrated oilfield asset modeling using multiple resolutions of reservoir detail |
US10346745B2 (en) * | 2013-09-05 | 2019-07-09 | International Business Machines Corporation | Method of using graphical index maps to provide automated relationship discovery and impact analyses |
US9790770B2 (en) | 2013-10-30 | 2017-10-17 | The Texas A&M University System | Determining performance data for hydrocarbon reservoirs using diffusive time of flight as the spatial coordinate |
JP6502133B2 (en) * | 2014-03-28 | 2019-04-17 | 日本碍子株式会社 | Porous body, honeycomb filter, method of manufacturing porous body, and method of manufacturing honeycomb filter |
US10443358B2 (en) * | 2014-08-22 | 2019-10-15 | Schlumberger Technology Corporation | Oilfield-wide production optimization |
GB2550741B (en) * | 2015-03-04 | 2020-11-11 | Landmark Graphics Corp | Path optimization in production network systems |
US10042915B2 (en) | 2015-09-28 | 2018-08-07 | International Business Machines Corporation | Semantic mapping of topic map meta-models identifying assets and events to include directionality |
US10387476B2 (en) | 2015-11-24 | 2019-08-20 | International Business Machines Corporation | Semantic mapping of topic map meta-models identifying assets and events to include modeled reactive actions |
JP6605951B2 (en) * | 2015-12-25 | 2019-11-13 | 株式会社東芝 | Simulation apparatus and simulation method |
US10274331B2 (en) | 2016-09-16 | 2019-04-30 | Polaris Industries Inc. | Device and method for improving route planning computing devices |
WO2018160279A1 (en) * | 2017-02-28 | 2018-09-07 | National Technology & Engineering Solutions Of Sandia, Llc (Ntess) | Methods for preventing computationally explosive calculations for model parameters distributed on a hierarchy of geometric simplices |
US11263801B2 (en) | 2017-03-31 | 2022-03-01 | Schlumberger Technology Corporation | Smooth surface wrapping of features in an imaged volume |
CN111712823A (en) * | 2017-12-14 | 2020-09-25 | 斯伦贝谢技术有限公司 | System and method for simulating reservoir models |
CN109339772B (en) * | 2018-09-10 | 2022-01-04 | 中国石油天然气股份有限公司 | Method and device for judging gas reservoir communication relation of well region |
WO2021038261A1 (en) * | 2019-08-23 | 2021-03-04 | Total Se | Method for determining locations of wells in a field |
CN110796383A (en) * | 2019-11-01 | 2020-02-14 | 长沙理工大学 | Water system connectivity evaluation index system considering ecological base flow |
CN110837608B (en) * | 2019-11-07 | 2024-04-12 | 中科天玑数据科技股份有限公司 | Public opinion topic propagation path analysis system and method based on multi-source data |
CN111983678B (en) * | 2020-07-22 | 2024-02-23 | 中海石油(中国)有限公司深圳分公司 | Method for rapidly evaluating development potential of deep water sand body |
EP4244754A1 (en) * | 2020-11-13 | 2023-09-20 | Chevron U.S.A. Inc. | Well correlation through intermediary well |
EP4244755A1 (en) * | 2020-11-13 | 2023-09-20 | Chevron U.S.A. Inc. | Visualization of directional influence of wells on correlation |
US20240061146A1 (en) * | 2020-11-13 | 2024-02-22 | Chevron U.S.A. Inc. | Subsurface characterization based on multiple correlation scenarios |
CN113034686B (en) * | 2021-03-29 | 2023-10-13 | 中国科学院新疆生态与地理研究所 | Stratum sequence acquisition method based on topological sorting |
US20230119861A1 (en) * | 2021-10-20 | 2023-04-20 | Chevron U.S.A. Inc. | Acceleration of subsurface simulation using filtered grid connectivity graph |
CN113960288B (en) * | 2021-10-28 | 2023-09-19 | 中国石油大学(华东) | Quantitative evaluation method for source-sink system |
US20230376354A1 (en) * | 2022-05-19 | 2023-11-23 | The Boeing Company | Hybrid-Computing Resource Optimization |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6862531B2 (en) * | 2002-01-15 | 2005-03-01 | Westerngeco, L.L.C. | Layer stripping converted reflected waveforms for dipping fractures |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5501273A (en) * | 1994-10-04 | 1996-03-26 | Amoco Corporation | Method for determining the reservoir properties of a solid carbonaceous subterranean formation |
US5732192A (en) | 1994-11-30 | 1998-03-24 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Global qualitative flow-path modeling for local state determination in simulation and analysis |
US5757663A (en) | 1995-09-26 | 1998-05-26 | Atlantic Richfield Company | Hydrocarbon reservoir connectivity tool using cells and pay indicators |
US6106561A (en) * | 1997-06-23 | 2000-08-22 | Schlumberger Technology Corporation | Simulation gridding method and apparatus including a structured areal gridder adapted for use by a reservoir simulator |
US6662146B1 (en) | 1998-11-25 | 2003-12-09 | Landmark Graphics Corporation | Methods for performing reservoir simulation |
DE60044492D1 (en) | 1999-03-31 | 2010-07-15 | Exxonmobil Upstream Res Co | METHOD FOR SIMULATING A CHARACTERISTICS OF A PHYSICAL SYSTEM |
US6477515B1 (en) | 1999-08-11 | 2002-11-05 | The United States Of America As Represented By The Secretary Of The Navy | Efficient computation of least cost paths with hard constraints |
US6549879B1 (en) * | 1999-09-21 | 2003-04-15 | Mobil Oil Corporation | Determining optimal well locations from a 3D reservoir model |
WO2001027858A1 (en) | 1999-10-13 | 2001-04-19 | The Trustees Of Columbia University In The City Of New York | Petroleum reservoir simulation and characterization system and method |
US6928399B1 (en) | 1999-12-03 | 2005-08-09 | Exxonmobil Upstream Research Company | Method and program for simulating a physical system using object-oriented programming |
AU2001270265A1 (en) | 2000-06-29 | 2002-01-14 | Object Reservoir, Inc. | Method and system for coordinate transformation to model radial flow near a singularity |
US7277836B2 (en) | 2000-12-29 | 2007-10-02 | Exxonmobil Upstream Research Company | Computer system and method having a facility network architecture |
WO2002057901A1 (en) | 2001-01-17 | 2002-07-25 | Exxonmobil Upstream Research Company | Simulation method and system using component-phase transformations |
FR2823877B1 (en) | 2001-04-19 | 2004-12-24 | Inst Francais Du Petrole | METHOD FOR CONSTRAINING BY DYNAMIC PRODUCTION DATA A FINE MODEL REPRESENTATIVE OF THE DISTRIBUTION IN THE DEPOSIT OF A CHARACTERISTIC PHYSICAL SIZE OF THE BASEMENT STRUCTURE |
US7062719B2 (en) | 2001-08-15 | 2006-06-13 | National Instruments Corporation | Graphically defining a route through one or more switch devices |
-
2004
- 2004-09-29 WO PCT/US2004/032942 patent/WO2005033739A2/en active Application Filing
- 2004-09-29 CA CA2538577A patent/CA2538577C/en not_active Expired - Fee Related
- 2004-09-29 CN CN200480028500A patent/CN100590637C/en not_active Expired - Fee Related
- 2004-09-29 US US10/568,940 patent/US7584086B2/en not_active Expired - Fee Related
- 2004-09-29 EP EP04789530A patent/EP1668561A2/en not_active Withdrawn
-
2006
- 2006-04-27 NO NO20061871A patent/NO20061871L/en not_active Application Discontinuation
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6862531B2 (en) * | 2002-01-15 | 2005-03-01 | Westerngeco, L.L.C. | Layer stripping converted reflected waveforms for dipping fractures |
Non-Patent Citations (5)
Title |
---|
AKINSANMI O.B. ET AL: 'Application of Azimuthal Density While Drilling Images for Dips, Facies and Reservoir Characterization-Niger/Delta Experience' SPE/PETROLIUM SOCIETY OF CIM 65460 08 November 2000, pages 1 - 14, XP008115596 * |
C.V. DEUTSCH: 'Improved Reservoir Management Through Ranking Stochastic Reservoir Models' SPE/DOE 35411 24 April 1996, pages 105 - 113, XP008115597 * |
LI Y.S. ET AL: 'A Probabilistic Approach from Reliability Assesment of Offshore Structure Corrision and Fatique' SPE 59435 26 April 2000, XP008115595 * |
NORRENA K. ET AL: 'Automatic Determination of Well Placement Subject to Geostatistical Economic Constraints' SPE/PETROLIUM SOCIETY IF CIM/CHOA 78996 07 November 2002, pages 1 - 12, XP002499489 * |
TANGAS J. ET AL: 'Integrated Reservoir and Surface Network Simulation in Reservoir Management of Southern North Sea Gas Reservoirs' SPE 50635 22 October 1998, pages 51 - 62, XP000957854 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NO337471B1 (en) * | 2006-03-02 | 2016-04-18 | Exxonmobil Upstream Res Co | Method for quantifying reservoir connectivity using fluid travel time |
EP1994488A2 (en) * | 2006-03-02 | 2008-11-26 | ExxonMobil Upstream Research Company | Method for quantifying reservoir connectivity using fluid travel times |
EP1994488A4 (en) * | 2006-03-02 | 2011-07-13 | Exxonmobil Upstream Res Co | Method for quantifying reservoir connectivity using fluid travel times |
US8776895B2 (en) | 2006-03-02 | 2014-07-15 | Exxonmobil Upstream Research Company | Method for quantifying reservoir connectivity using fluid travel times |
WO2008021743A3 (en) * | 2006-08-09 | 2008-05-08 | Schlumberger Services Petrol | Facilitating oilfield development with downhole fluid analysis |
WO2008021743A2 (en) * | 2006-08-09 | 2008-02-21 | Services Petroliers Schlumberger | Facilitating oilfield development with downhole fluid analysis |
EP2232406A4 (en) * | 2007-12-21 | 2017-10-18 | Exxonmobil Upstream Research Company | Method and apparatus for analyzing three-dimensional data |
EP2098887A1 (en) | 2008-03-03 | 2009-09-09 | Exxonmobil Upstream Research Company | Determining connectivity architecture in 2-D and 3-D heterogeneous data |
NO337774B1 (en) * | 2009-05-07 | 2016-06-20 | Saudi Arabian Oil Co | Well pressure system and method |
US8731873B2 (en) | 2010-04-26 | 2014-05-20 | Exxonmobil Upstream Research Company | System and method for providing data corresponding to physical objects |
WO2015137916A1 (en) * | 2014-03-10 | 2015-09-17 | Landmark Graphics Corporation | Modeling geologic surfaces using unilateral non-node constraints from neighboring surfaces in the stratigraphic sequence |
AU2014386269B2 (en) * | 2014-03-10 | 2016-10-27 | Landmark Graphics Corporation | Modeling geologic surfaces using unilateral non-node constraints from neighboring surfaces in the stratigraphic sequence |
US9633475B2 (en) | 2014-03-10 | 2017-04-25 | Landmark Graphics Corporation | Modeling geologic surfaces using unilateral non-node constraints from neighboring surfaces in the stratigraphic sequence |
Also Published As
Publication number | Publication date |
---|---|
CN1860489A (en) | 2006-11-08 |
EP1668561A2 (en) | 2006-06-14 |
WO2005033739A3 (en) | 2006-04-06 |
US20070027666A1 (en) | 2007-02-01 |
NO20061871L (en) | 2006-04-27 |
CA2538577A1 (en) | 2005-04-14 |
CA2538577C (en) | 2010-08-03 |
US7584086B2 (en) | 2009-09-01 |
CN100590637C (en) | 2010-02-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2538577C (en) | Characterizing connectivity in reservoir models using paths of least resistance | |
US9874648B2 (en) | Reservoir connectivity analysis in a 3D earth model | |
US10838108B2 (en) | Method and system for interpolating discontinuous functions in a subsurface model | |
US9733388B2 (en) | Systems and methods for connectivity analysis using functional objects | |
US8386225B2 (en) | Method of optimizing the development of a fluid reservoir by taking into account a geologic and transient exchange term between matrix blocks and fractures | |
US20130138406A1 (en) | Method for constructing a fracture network grid from a voronoi diagram | |
EP1994488A2 (en) | Method for quantifying reservoir connectivity using fluid travel times | |
US10895131B2 (en) | Probabilistic area of interest identification for well placement planning under uncertainty | |
CA3012429C (en) | Hybrid 3d geocellular representation of selected natural fracture network subsets | |
Hazlett et al. | Optimal well placement in heterogeneous reservoirs through semianalytic modeling | |
US8718992B2 (en) | Method for history matching of a geological model comprising a sub-seismic fault network | |
US9329289B2 (en) | Method of constructing a grid representative of a property distribution by conditional multipoint statistical simulation | |
WO2016187048A1 (en) | Prospect assessment using structural frameworks | |
AlQassab et al. | Estimating the size and orientation of hydraulic fractures using microseismic events | |
Le Ravalec et al. | Integrating data of different types and different supports into reservoir models | |
Oliveira et al. | Well placement subclustering within partially oil-saturated flow units | |
CN112292714B (en) | Grid partition based on fault radiation | |
Litvak et al. | Uncertainty Estimation in Production Predictions Constrained by Production History and Time-Lapse Seismic in a GOM Oil Field | |
US10209380B2 (en) | Methods and systems for juxtaposition across geological discontinuities | |
EP3926368B1 (en) | Analyzing a hydrocarbon trap | |
Stright et al. | Coupled geological modeling and history matching of fine-scale curvilinear flow barriers | |
Nwachukwu | Model selection for CO₂ sequestration using surface deflection and injection data | |
AlQassab | Modeling hydraulic fractures using microseismic events | |
Madalimov et al. | A Method of Reservoir Connectivity Analysis for Single-and Dual-Porosity/Dual-Permeability Models | |
Lecour et al. | Consistent Structural Model Simulations for Volumetric Uncertainty Study |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 200480028500.0 Country of ref document: CN |
|
AK | Designated states |
Kind code of ref document: A2 Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW |
|
AL | Designated countries for regional patents |
Kind code of ref document: A2 Designated state(s): GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LU MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application | ||
WWE | Wipo information: entry into national phase |
Ref document number: 2007027666 Country of ref document: US Ref document number: 10568940 Country of ref document: US |
|
WWE | Wipo information: entry into national phase |
Ref document number: 1053/DELNP/2006 Country of ref document: IN |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2538577 Country of ref document: CA |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2004789530 Country of ref document: EP |
|
WWP | Wipo information: published in national office |
Ref document number: 2004789530 Country of ref document: EP |
|
WWP | Wipo information: published in national office |
Ref document number: 10568940 Country of ref document: US |