US20160090825A1 - Method and System for Analyzing the Uncertainty of Subsurface Model - Google Patents
Method and System for Analyzing the Uncertainty of Subsurface Model Download PDFInfo
- Publication number
- US20160090825A1 US20160090825A1 US14/816,729 US201514816729A US2016090825A1 US 20160090825 A1 US20160090825 A1 US 20160090825A1 US 201514816729 A US201514816729 A US 201514816729A US 2016090825 A1 US2016090825 A1 US 2016090825A1
- Authority
- US
- United States
- Prior art keywords
- realizations
- conceptual
- instantiated
- subsurface
- instructions
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 150
- 238000004088 simulation Methods 0.000 claims description 41
- 238000004519 manufacturing process Methods 0.000 claims description 20
- 230000035699 permeability Effects 0.000 claims description 15
- 230000015572 biosynthetic process Effects 0.000 claims description 13
- 238000005755 formation reaction Methods 0.000 claims description 13
- 230000003993 interaction Effects 0.000 claims description 13
- 230000003750 conditioning effect Effects 0.000 claims description 11
- 230000003094 perturbing effect Effects 0.000 claims description 9
- 238000004891 communication Methods 0.000 claims description 7
- 229930195733 hydrocarbon Natural products 0.000 abstract description 22
- 150000002430 hydrocarbons Chemical class 0.000 abstract description 22
- 239000004215 Carbon black (E152) Substances 0.000 abstract description 11
- 238000011161 development Methods 0.000 abstract description 6
- 230000008569 process Effects 0.000 description 33
- 238000010586 diagram Methods 0.000 description 27
- 239000012530 fluid Substances 0.000 description 21
- 238000004458 analytical method Methods 0.000 description 17
- 239000011435 rock Substances 0.000 description 15
- 239000011159 matrix material Substances 0.000 description 14
- 230000008859 change Effects 0.000 description 13
- 230000000694 effects Effects 0.000 description 12
- 238000012937 correction Methods 0.000 description 11
- 238000012986 modification Methods 0.000 description 10
- 230000004048 modification Effects 0.000 description 10
- 238000013213 extrapolation Methods 0.000 description 9
- 230000002441 reversible effect Effects 0.000 description 9
- 230000009897 systematic effect Effects 0.000 description 9
- 239000003795 chemical substances by application Substances 0.000 description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 8
- 230000001143 conditioned effect Effects 0.000 description 6
- 238000009826 distribution Methods 0.000 description 6
- 208000035126 Facies Diseases 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 5
- 230000018109 developmental process Effects 0.000 description 5
- 230000007246 mechanism Effects 0.000 description 5
- 238000011084 recovery Methods 0.000 description 5
- 238000000151 deposition Methods 0.000 description 4
- 230000008021 deposition Effects 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 239000004576 sand Substances 0.000 description 4
- 238000003860 storage Methods 0.000 description 4
- 230000003190 augmentative effect Effects 0.000 description 3
- 230000006399 behavior Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000004590 computer program Methods 0.000 description 3
- 239000004020 conductor Substances 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 238000012804 iterative process Methods 0.000 description 3
- 238000005192 partition Methods 0.000 description 3
- 239000011148 porous material Substances 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000002829 reductive effect Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 238000005309 stochastic process Methods 0.000 description 3
- 230000007704 transition Effects 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 208000018747 cerebellar ataxia with neuropathy and bilateral vestibular areflexia syndrome Diseases 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000013400 design of experiment Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 230000035876 healing Effects 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000003012 network analysis Methods 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000008929 regeneration Effects 0.000 description 1
- 238000011069 regeneration method Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000002948 stochastic simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009885 systemic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- 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
-
- G01V99/005—
-
- 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
- E21B41/00—Equipment or details not covered by groups E21B15/00 - E21B40/00
Definitions
- This invention relates generally to the field of hydrocarbon exploration and extraction. Specifically, the invention is a method to examine uncertainty and risk associated with the development of a hydrocarbon resource by rapidly generating and analyzing variations of subsurface models realized from different scenarios.
- a geologic model is formed that includes various static properties.
- a reservoir model is created to model dynamic properties.
- the reservoir model is formed from geological horizons and faults.
- the reservoir model includes a framework that establishes the geometrical foundation for the three-dimensional grid and provides some of the boundaries for facies and petrophysical models that describe rock properties and contained fluids.
- the resulting reservoir model forms the basis for volumetric computations, reservoir simulations, facilities planning computations and well planning computations. While seismic and well data provide information regarding the reservoir model, considerable uncertainty may remain regarding the reservoir model.
- the effect of uncertainty is often examined through various means.
- the uncertainty may be examined by perturbing uncertain aspects or features of the reservoir model, recomputing the quantity of interest, and examining sensitivity of the quantity of interest with regard to the uncertain aspects.
- the problem with framework uncertainty related to the geometric foundation of the three-dimensional grid is that the steps between framework generation, definition of the grid, and computation of the quantity of interest are computationally and labor intensive, often requiring user input.
- Some conventional methods only perturb the depth of different model objects, such as faults and horizons.
- the depth perturbation may be spatially variable, for example allowing the flanks of an anticline to be pushed down or pulled up, while leaving the crest unperturbed.
- the resulting flexing of faults, horizons and other model objects occurs in the vertical direction, however, and the modeling grid is flexed simultaneously. That is, the geometry of the grid changes, but not its structure.
- Providing vertical and lateral movement of the model objects typically requires grid regeneration. Some conventional methods attempt to correct the existing grid. Other conventional methods move the model objects vertically and laterally and then adjust the intersections. The adjustment has two aspects, however, the removal of the original intersections and the implementation of the new intersections. For example, when a fault is shifted, a portion of a bisected horizon moves from the foot-wall side to the hanging-wall side or vice versa. The moving horizon piece has a previous displacement that should be undone. Typically, the undoing of previous displacements appears to be performed by deletion of horizons near faults and extrapolation towards the shifted fault.
- an enhancement to exploration and reservoir delineation techniques is needed to identify and recover hydrocarbons in light of imprecise data.
- the present techniques provide a streamlined method for generation of a perturbed framework and thus a perturbed reservoir model.
- the enhancements may provide a method for systematic removal of the effects of faults and folds to provide numerous realizations of the reservoir model in an efficient manner.
- the enhancements may provide a method to explore fault connectivity or check for geologic plausibility and technical validity which removes problematic model perturbation.
- a method includes analyzing uncertainty of subsurface formations and includes: creating a conceptual subsurface model, wherein the conceptual subsurface model is associated with a subsurface formation and comprises a plurality of objects; selecting parameter ranges for each of the plurality of objects and interactions between two or more of the plurality of objects; instantiating realizations for the plurality of objects based on the selected parameter ranges; and combining instantiated realizations of the plurality of objects into a reservoir model.
- a computer system for analyzing uncertainty of subsurface formations for production or exploration operations.
- the computer system may include a processor; memory in communication with the processor; and a set of instructions stored in memory and accessible by the processor.
- the set of instructions when executed by the processor, are configured to: create a conceptual subsurface model, wherein the conceptual subsurface model is associated with a subsurface formation and comprises a plurality of objects; select parameter ranges for each of the plurality of objects and interactions between two or more of the plurality of objects; instantiate realizations for the plurality of objects based on the selected parameter ranges; and combine instantiated realizations of these objects into a reservoir model.
- FIG. 1 is a flow chart for performing hydrocarbon exploration in accordance with an exemplary embodiment of the present techniques.
- FIG. 2 is a diagram of a realization of a geologic model.
- FIG. 3 is a diagram of a concept of the realization of FIG. 2 .
- FIG. 4 is a diagram of a schematic application of the undeforming and redeforming a framework in accordance with an exemplary embodiment of the present techniques.
- FIGS. 5A, 5B and 5C are diagrams of implausible realizations.
- FIGS. 6A, 6B and 6C are diagrams of unfaulting during concept creation.
- FIGS. 7A, 7B and 7C are diagrams of refaulting during instantiation of a realization.
- FIG. 8 is a diagram of the process from a base framework realization to an instantiated framework realization in accordance with an exemplary embodiment of the present technique.
- FIG. 9 is a block diagram of a computer system that may be used to perform any of the methods disclosed herein.
- Subsurface model is a reservoir model or a geologic model.
- “Geologic model”, as used herein, is three-dimensional model of the subsurface having static properties and includes faults, horizons, facies, lithology and properties such as porosity, permeability, or the proportion of sand and shale.
- Reservoir model is a three-dimensional model of the subsurface that in addition to static properties such as porosity and permeability also has dynamic properties that vary over the timescale of resource extraction such as fluid composition, pressure, and relative permeability.
- “Framework”, as used herein, is a geologic model formed from faults, horizons, model boundaries, and facies boundaries, e.g., a geologic model containing only surfaces and polylines.
- Concept is a model containing faults, horizons, and facies boundaries without geometry or location. Aspects of the geologic or reservoir models related to geometry and location are suppressed. For example, a concept is often conveyed as a sketch of only the objects deemed most pertinent, omitting anything else.
- “Scenario”, as used herein, is a concept or partial subsurface model in combination with select parameters and their ranges used to build realizations of subsurface models by deterministically or stochastically varying these parameters. Examples may be normal faults versus reverse faults, different time-depth conversions or environments of deposition, or high net-to-gross versus low-net-to-gross regions.
- Realization is a subsurface model (e.g., geologic model) created from a concept or scenario by assigning geometry and location to faults, horizons, and boundaries; and values to properties which may be utilized for computations and quantitative queries.
- geologic model e.g., geologic model
- “Instantiate”, as used herein, is the process of transforming a (qualitative) concept or object thereof to a (quantitative) realization or object thereof. Instantiation may be performed by interpolation or extrapolation of measurements or by application of a stochastic process.
- Simulation is the process of making a prediction related to the resource extraction based on the reservoir model.
- a simulation is typically performed by execution of a reservoir-simulator computer program on a processor, which computes composition, pressure, or movement fluid as function of time and space for a specified scenario of injection and production wells by solving a set of reservoir fluid flow equations.
- Unfaulting operation is the process of undoing the effects of a fault on the other objects of the model.
- a horizon typically exhibits a discontinuity where intersected by a fault.
- This fault is turned into a conceptual fault, its location and geometry are removed.
- Remaining model objects still exhibit the discontinuities caused by the original fault.
- the unfaulting operation removes these discontinuities, returning the remaining objects to a state unaffected by this fault, e.g., to a state where the fault never existed.
- faulting operation is the process of affecting the other model objects of a model when a fault is realized.
- the typical effect is the generation of discontinuities where the fault and existing objects intersect.
- Unfolding operation is process of at least partially removing location and geometry from a horizon transforming a horizon into a conceptual horizon.
- An unfolding operation removes the effects of a folding event (continuous deformation) from the model.
- Refolding operation is the process of applying a continuous deformation to the objects, either by assignment or perturbation of geometry and location to horizons and specified other objects.
- a subsurface model which may include a reservoir model or geologic model, is a computerized representation of a subsurface region based on geophysical and geological observations made on and below the surface of the Earth. It is the numerical equivalent of a three-dimensional geological map complemented by a description of physical quantities in the domain of interest.
- Subsurface models which are reservoir models, are often used as inputs to reservoir simulation programs that predict the behavior of rocks and fluids contained therein under various scenarios of hydrocarbon recovery. When producing an actual hydrocarbon reservoir, miscalculations or mistakes can be costly.
- Using subsurface models in simulations provides a mechanism to identify which recovery options offer the safest and most economic, efficient, and effective development plans for a particular reservoir.
- Construction of a subsurface model is typically a multistep process.
- First, a structural model or structural framework is created from surfaces that include faults, horizons, and if necessary, additional surfaces that bound the area of interest for the geologic model.
- the different surfaces define closed volumes often called zones.
- Second, each zone is meshed or partitioned into small cells defined by a three-dimensional grid.
- properties are assigned to surfaces (e.g., transmissibility) and individual cells (e.g., rock type, porosity, permeability, or oil saturation).
- each cell in the model is assigned a rock type.
- the cells may be beach sand, high water energy marine upper shoreface sand, intermediate water energy marine lower shoreface sand, and deeper low energy marine silt and shale.
- the distribution of these rock types within the model may be controlled by several methods, including map boundary polygons, rock type probability maps, or statistically emplaced based on concepts. Where available, rock type assignment may be conditioned to well data.
- Reservoir quality parameters typically include porosity and permeability, but may include measures of clay content, cementation factors, and other factors that affect the storage and deliverability of fluids contained in the pores of those rocks.
- Geostatistical techniques are typically used to populate the cells with porosity and permeability values that are appropriate for the rock type of each cell.
- Rock pores are saturated with groundwater, oil or gas. Fluid saturations may be assigned to the different cells to indicate which fraction of their pore space is filled with the specified fluids. Fluid saturations and other fluid properties may be assigned deterministically or geostatistically.
- Geostatistics is useful in modeling to interpolate observed data and to superimpose an expected degree of variability.
- kriging which uses the spatial correlation among data and intends to construct the interpolation via semi-variograms, may be used.
- geostatistical simulation is often used, for example based on variograms, training images, or parametric geological objects.
- Perturbing surface properties or cell properties, such as rock type, reservoir properties or fluid properties, is a conventional process, which may utilize deterministic or geostatistical methods to assign them. The assignment may include choosing a different variogram for kriging or a different seed for geostatistical simulation.
- a realization is instantiated from a concept.
- the difference between a concept and a realization is a complete geometry.
- objects such as points, polylines, polygons, horizons, faults, or compartments, have locations, relative positions with regard to each other, shapes, or sizes.
- the topology of the objects (e.g., their interactions) of the realization is defined.
- a property attached to an object has values at essentially every location of the object.
- a realization does not contain free parameters anymore.
- a concept contains free parameters relating to topology, geometry, and properties.
- a concept does not have geometry associated with its objects.
- At least some of the points, polylines, polygons, horizons, faults, or compartments do not have their locations, relative positions with regard to each other, shapes, or sizes defined.
- the topology of its objects can be completely specified, partially specified, or unspecified. For example, the order in which different horizons and faults terminate may be undefined. In some embodiments, the nature or interpretation of an object may be undefined.
- a horizon is (either implicitly or explicitly) typed conformable, basal, topal, erosional, or discontinuous while in a concept, the horizon type can but does not need to be specified.
- a fault is typed normal, reverse, or strike-slip, while in a concept, the fault type can but does not need to be typed.
- the cardinality of an object can be undefined.
- a single fault in the concept may be realized as multiple faults, for example, by realizing a fault as a set of relay faults or as a set of parallel faults.
- Framework uncertainty is caused by uncertainty in seismic migration, time-depth (or depth-true depth) conversion, structural interpretation, fault positions, well picks, horizon correlation and interpretation, and layer thicknesses.
- framework changes tend to create artifacts that have to be addressed. Accordingly, the present techniques provide a method and system for modifying an existing watertight (e.g., sealed) framework in such a manner that the resulting framework is yet again watertight without gaps or overlaps between model cells or reservoir compartments.
- the present techniques involve a workflow that may be used to analysis uncertainty and risk associated with the development of hydrocarbon resources by rapidly generating and analyzing variations of subsurface models realized from scenarios.
- structural artifacts may be removed once by systematic removal of the effects of faults and folds.
- numerous realizations of the subsurface model can then be generated by starting with an initial starting subsurface model and the corresponding concept and then, iteratively computing the effects of faults and folds on the revised starting subsurface model.
- the present techniques may also include the ability to omit a fault or fold from the perturbed model or to introduce additional faults and folds into the perturbed model. This capability allows the exploration of fault connectivity, for example the substitution of one long contiguous fault by a set of fault relays. Further still, the method may also include a performance of a check for geologic plausibility and technical validity which removes problematic model perturbations from the workflow or at least from the set of simulation results used for the final analysis.
- the workflow may include various steps. For example, as a first step, each horizon is unfaulted one fault at a time. Each unfaulting operation removes a tear in a horizon, bringing its edges back together. A fault causes a discontinuity in a given horizon. Unfaulting removes this discontinuity by adjusting the depths of the horizon near fault and propagating or extrapolating these adjustments away from the fault along the horizon. As a second step, the overarching remaining structure of now fault-free horizons is removed. This unfolding is performed by replacing the horizon with an approximately planar one while recording the necessary vertical depth adjustments. As a third step, a different overarching structure is imposed on the horizons.
- This refolding is performed by replacing the approximately planar horizons with differently shaped ones, preferably guided by the previously recorded adjustments.
- the horizons are refaulted.
- the refaulting is performed by moving the horizon on the different sides of the fault to the desired relative location and propagating these adjustments away from the fault along the horizon.
- Unfaulting (F ⁇ 1), refaulting (F), unfolding (S ⁇ 1), and refolding (S) may be viewed as mappings, transforms, and/or operators suggesting the notation F ⁇ 1, F, S ⁇ 1, S.
- the resulting framework should be identical to the existing framework. However, if the folding operator and/or the faulting operator are modified, then the resulting framework is a perturbation of the existing one.
- the modifications may include, but are not limited to: shift, rotate, scale or deform fault; change the throws or fault type; split one fault, combine two faults, or replicate a fault to distribute the throws; shift or deform a surface, change layer thickness or lateral change rates; change the thickness between two surfaces or lateral thickness changes; or change the topology.
- Unfaulting, refaulting, unfolding, and refolding can be performed with different methods depending on the desired degree of accuracy.
- the methods range from purely geometric methods; to kinematic methods that attempt to preserve distances, areas, and volumes; to geomechanical methods that model stresses, strains, elasticity, plasticity, failure, etc.
- the operators may vary for different embodiments. For example, purely geometrical operators may suffice because similar assumptions are made for unfaulting (e.g., unfolding) and refaulting (e.g., refolding). Artifacts introduced by application of simplistic inverse operators are largely removed when applying the similarly simplistic forward operators. Because the objective of the workflow is generation of a perturbed framework, any remaining artifacts may be considered part of the perturbation.
- the geometrical operators may include unfaulting geometrical operators, faulting geometrical operators, folding geometrical operators and refolding geometrical operators.
- the unfaulting geometrical operators may be constructed from fault-horizon-intersection polygons by estimating an intermediate polyline from the hanging- and footwall polylines. Given these three polylines and the constraint that perturbations in the far field are minimal, an operator in form of an elevation perturbation map or delta-z map can be constructed.
- Application of a map-based depth modification is conventional functionality in many commercially available geologic modeling packages. Removing a first of multiple faults changes the intersections between any given horizon and the remaining faults.
- fault-horizon-intersection polygons may need to be perturbed with the first unfaulting operator before removal of a second fault and repeated for other similar operations.
- Faulting geometrical operators may be constructed by specification of a throw distribution for a given fault, constructing the horizon fault-intersection polygon for the specified fault and any given horizon, constructing an elevation perturbation map conditioned on this polygon and constrained to minimize far field perturbation, and applying this perturbation map or delta-z map to the specified horizon.
- Folding and refolding geometrical operators may be constructed, for example, by geostatistics. Fault geometry may be changed deterministically or geostatistically.
- Any other object that is specified by coordinates may also be transferred from one framework to another by application of some or all of these operators.
- Examples may include well picks, well paths, generic polylines, or geobodies that may be used for conditioning objects of the subsurface model. For example, perturbed horizons may be conditioned to well picks that themselves may be perturbed to account for their uncertainty.
- the parameters (or geometry) of these operations is recorded to provide their application to some or all of the other objects in the model.
- the method may include instantiating a realization of a subsurface model (e.g., reservoir model).
- the method may involve the creation of a conceptual subsurface model, the selection of parameter ranges for the various objects of the conceptual model and their interactions, and combining instantiated realizations of these objects into a subsurface model.
- the conceptual model may be generated by an agent (a user or a computer program that acts on behalf of a user) using a concept editor; may be automatically created from an inputted base realization; and/or may be created from the base case by systematically undoing faults and/or faults.
- the parameter ranges may be estimated from the undoing of faults and/or folds. Also, the unfaulting and refaulting may be performed based on fault-horizon cutoff polygons.
- the instantiated realizations may be analyzed for geologic plausibility and, if warranted, rejected; may be analyzed for technical consistency and, if warranted, rejected; may be further augmented with properties such as porosity, permeability, or oil saturation; and/or may be simulated.
- the simulation may be performed using a simulation proxy method.
- a set of simulations related to different scenarios and/or realizations may be summarized with a statistics and/or may be analyzed to affect a decision.
- a connectivity measure may be used as a simulation proxy; and this connectivity measure may be based on a graph-based centrality measure. The centrality measure, which is described in U.S.
- 14/272,581 which is incorporated by reference, may include one or more of degree, betweenness, closeness, and eigenvector. Further, the method may include ranking the plurality of objects, instantiated realizations, or other items in order of the respective centralization measures.
- the method may include different concepts or a single concept.
- the method may involve creating a concept once and multiple realizations that are created from the concept.
- a concept may be a general layout of objects and relationships deemed to be of higher importance. The concept may not be drawn to scale, may not even show relative geometry or shape, and omit substantial amounts of detail.
- a realization is a model having more detail than a concept, which may be generated by a deterministic or stochastic process.
- the method may also include the creation of concepts that serve as scenarios.
- the method may also include multiple concepts that are created from multiple inputted base realizations to serve as different scenarios.
- FIGS. 1 to 10 Various aspects of the present techniques are described further in FIGS. 1 to 10 .
- FIG. 1 is a flow chart for performing hydrocarbon exploration or delineation of a potential hydrocarbon resource in accordance with an exemplary embodiment of the present techniques.
- the method may include creating a concept and forming various realizations based on the concept.
- some blocks may be omitted, repeated, performed in different order, or augmented with additional blocks not shown. That is, some blocks may be performed sequentially, while others are executed simultaneously in parallel.
- This flowchart 100 begins with inputting a base realization, as shown in block 101 .
- the base realization includes is a subsurface model having various properties.
- the concept is created.
- the creation of the concept may include reducing the base realization to a concept by systemic removal of all deformation.
- the removal may include unfaulting and unfolding.
- parameter ranges are selected.
- the selection of parameter ranges may be performed by an interpreter or computer algorithm.
- the selection may include determining bounds for deformations and other model parameters. Interactions between two or more of the plurality of objects may also be defined at this point in the process.
- the realizations are instantiated.
- the instantiating the realizations may include relying upon the parameter ranges (e.g., drawing parameters from these bounds) and selection of the interactions between two or more of the plurality of objects. Multiple realizations of the same concept or scenario may be obtained by systematic variation of parameters within their parameter ranges or by random sampling of the parameter ranges using a stochastic process.
- a plausibility or validity check may be performed, as shown in block 105 .
- the realization is checked for technical validity and/or geologic plausibility. For example, the realization may be verified by analyzing the validity of the resulting grid or by examination of fault polygons or thickness maps. Realizations that fail this check are either modified or discarded.
- the verification may include determining if the realization is geologically implausible and, if so, discarding the instantiated realizations that are geologically implausible.
- the verification may include determining if the realization is technically consistency and, if so, discarding the instantiated realizations that are technically inconsistency.
- geologically implausible realizations include those with stratigraphically older horizons conformably disposed over younger horizons, faults interpreted as belonging to a prior episode of deformation offset faults known to have moved during a later episode, and/or faults having displacement-to-map length ratios outside of bounds defined by interpretation of real faults systems.
- technical inconsistency include, gaps or holes in what should be continuous horizon representations, duplicated portions of horizons or faults, horizons or faults that loop back upon themselves, mesh triangles that face the wrong directions and/or grid cells that are inside-out, mesh triangles with four or more vertices, or mesh triangles that intersect each other without an explicitly represented intersection edge.
- the realization is then populated with properties, as shown in block 106 .
- the properties may include lithology, facies, porosity, permeability, fluid composition, and/or pressure.
- the realization is simulated, as shown in block 108 .
- the simulation of the realization may include pressure or saturation changes as functions of spatial position and time, bypassed or disconnected resources, or the fluid composition produced at a specified well.
- the determination may include creating a specific number of realizations, which may cover a specified part of the parameter space, cover a specified part of the response space, or design of experiments techniques to characterize the range of responses. If the determination is to perform another realization, the selection of the parameter ranges may be performed in block 103 . However, if another realization is not to be performed, a set of simulations of the realization(s) is analyzed, as shown in block 110 .
- the statistical or visual analysis may include ranking, whisper plots, box plots or other methods to review the different realizations.
- the analysis may include ranking the realizations based on a specified response, such as expected ultimate recovery or the maximum amount of oil, gas or water produced for a specified period of time.
- hydrocarbons are identified and produced, as shown in blocks 111 and 112 .
- hydrocarbons may be identified based at least partially on the analysis of the realizations.
- the realizations may be integrated with other measured data or subsurface models of the subsurface regions below the survey region. These different types of data may be integrated based on location information associated with the respective data to lessen uncertainty associated with the existence of hydrocarbons.
- the identified hydrocarbons may be produced, as shown in block 112 .
- various production operations may be performed to produce the hydrocarbons. For example, the operations may include drilling of a well to provide access to the hydrocarbon accumulation.
- the production may include installing a production facility configured to monitor and produce hydrocarbons from the production intervals that provide access to the hydrocarbons in the subsurface formation.
- the production facility may include one or more units to process and manage the flow of production fluids, such as hydrocarbons and/or water, from the formation.
- the production equipment and operations may be based on the realizations.
- the production facility may be coupled to a tree and various control valves via a control umbilical, production tubing for passing fluids from the tree to the production facility, control tubing for hydraulic or electrical devices, and a control cable for communicating with other devices within the wellbore.
- the present techniques provide a mechanism to provide a master subsurface model that may be used to generate other subsurface models.
- the present techniques involves generating perturbations from the concept to form various realizations. Then, the present techniques involve verification of the realizations (e.g., that the realizations are proper models). This provides a mechanism to lessen contamination by implausible models during subsequent statistical analysis.
- the geologic model may be divided into a framework and content.
- the framework is formed by collections of volumes (three-dimensional), their bounding surfaces as well as other surfaces (two-dimensional), polylines or curves (one-dimensional), and points (zero dimensional) that are embedded in a three-dimensional space.
- Surfaces relate to an area of interest, faults, and horizons.
- Curves relate to surface intersections, such as fault-horizon intersections, fault sticks, or polygons and polylines, separating gross geologic features, such as environments of deposition.
- Horizons partition the model into zones; while faults partition the model into segments. Faults, horizons, and polygons partition the model into compartments.
- Content refers to the properties associated with compartments (e.g., three dimensional distribution of properties), surfaces (e.g., two-dimensional distribution of properties), and polylines (e.g., one-dimensional distribution of properties).
- compartments e.g., three dimensional distribution of properties
- surfaces e.g., two-dimensional distribution of properties
- polylines e.g., one-dimensional distribution of properties.
- FIG. 2 and FIG. 3 These representations exemplify some of the differences between a realization, as shown in FIG. 2 , and a concept, as shown in FIG. 3 .
- FIG. 2 is a diagram of a realization 200 of a geologic model.
- a realization 200 includes a frame of reference or coordinate system, as indicated by the coordinate axes 201 and 202 .
- the realization 200 has an area of interest or region of interest 203 that specifies the spatial extent of the model. Inside this region of interest 203 , the realization 200 is completely quantified to enable numerical simulations. Outside this region of interest 203 , the realization 200 is not specified and/or quantified. There may or may not be any data or information available about the region outside of the region of interest 203 ; but the region outside of the region of interest 203 is irrelevant because it is not modeled.
- the realization 200 includes faults 204 and 205 (e.g., faults 205 a , 205 b , and 205 c ).
- faults 204 and 205 e.g., faults 205 a , 205 b , and 205 c .
- the minor fault 205 is realized by a set of parallel faults 205 a , 205 b , and 205 c.
- the realization 200 also includes horizons 206 and 207 . Both realized horizons 206 and 207 have attached geometries that describe depth, shape, and displacement caused by fault 204 .
- the horizons 206 and 207 create various zones, such as zones 208 , 209 , and 210 , which are bound by either the area of interest 203 or the horizons 206 or 207 .
- surfaces and zones, such as zones 208 , 209 , and 210 may have attached properties.
- zone 209 has an attached property 211 indicated by the gradual shading, such as porosity, net-to-gross ratio, or hydrocarbon saturation.
- the realized properties may be specified on a grid or mesh, such as mesh 212 .
- the realization contains the necessary information to perform a specified computation or simulation, such as the estimation of the gross rock volume (GRV), the stock tank original oil in place (STOOIP), the expected ultimate recovery (EUR) or the prediction of water cuts.
- GRV gross rock volume
- STOOIP stock
- FIG. 2 is a diagram of a concept 300 of the realization 200 of FIG. 2 .
- the concept 300 includes faults 304 and 305 , horizons 306 and 307 , and zones 308 , 309 and 310 . Because geometry in a concept 300 is typically unspecified, no frame of reference is utilized and horizons and faults are indicated by lines or planes, as shown by faults 304 and 305 and horizons 306 and 307 . In this diagram, fault 304 appears to be to the left of fault 305 , but without geometry, the spatial arrangement of the faults cannot be specified.
- constraints on the spatial arrangement may subsequently be imposed with the selection of parameter ranges, as noted in block 103 of FIG. 1 .
- horizon 306 appears to overlay horizon 307 , but again, without geometry, the spatial arrangement of the horizons cannot be specified.
- constraints on the spatial arrangement may subsequently be imposed with the selection of parameter ranges in block 103 of FIG. 1 .
- the topology is augmented with the concepts of younger (shallower) and older (deeper) to capture the relative order of horizons in the conceptual model.
- Horizons 306 and 307 are preferably typed as ‘base’, ‘top’ or ‘erosional’, ‘conformable’, or ‘unconformable’ or ‘discontinuous’.
- Zone 309 is bound, capped by horizon 306 , and based by horizon 307 .
- Zone 308 is unbound and based by horizon 306 .
- Zone 310 is unbound and capped by horizon 307 .
- Faults 304 and 305 are preferably typed as ‘normal’, ‘reverse’, or ‘strike-slip’.
- conceptual faults are further attributed with attributes such as ‘major’ or ‘minor’. If fault 304 is attributed with ‘major’ while fault 305 is attributed ‘minor’, then a realization of 304 truncates a realization of fault 305 in the event they intersect.
- the comparison of these diagrams exemplifies some of the differences between a realization and a concept.
- the conceptual faults 304 and 305 are realized as faults 204 and 205 (e.g., faults 205 a , 205 b , and 205 c ).
- the realization 200 there is a crosscutting relationship between the two faults where 204 is dominant or major, while, in the concept 300 , both faults 304 and 305 are typed ‘normal’.
- conceptual horizons 306 and 307 are realized as horizons 206 and 207 .
- Both realized horizons such as horizons 206 and 207 , have attached geometries that describe depth, shape, and displacement caused by fault 204 , while the conceptual horizons do not have such properties.
- zones 208 , 209 , and 210 are bound by either the area of interest 203 or the horizons 206 or 207 .
- surfaces and zones may have attached properties, while the conceptual zones do not include properties.
- the geologic model may be divided into a framework and content.
- the process of instantiating a realization is separated into two steps, which are i) instantiation of a framework realization and ii) instantiation of property realizations.
- the present techniques involve the concept creation (e.g., blocks 101 and 102 of FIG. 1 ) and instantiation of framework realizations (e.g., blocks 103 and 104 ).
- the concept is created directly with a suitable tool.
- the concept is created from a base realization (e.g., a realized geologic model that is stripped of geometry and possibly parts of its topology, meaning, and interpretation).
- this base realization is the most likely reservoir model, a model synthesized from the optimal interpretation, or the model is the statistically centralized (e.g., in the middle of the groupings) with regard to a specified prediction.
- FIG. 4 is a diagram 400 of a schematic application of the undeforming and redeforming a framework in accordance with an exemplary embodiment of the present techniques.
- the one (re)faulting operator and two (re)folding operators are utilized.
- the operators are perturbed, and their sequential order is commuted to create a realization caused by a different sequence of deformation events.
- Solid curves denote instantiated or realized objects, while dashed curves indicate conceptual objects.
- Block 402 depicts the base realization.
- Block 404 is the unfaulting operation F ⁇ 1 that removes the discontinuities imposed by the fault from the horizons and renders the fault conceptual. The result of this operation is shown in block 406 .
- Block 408 is a first unfolding operation S1 ⁇ 1 that removes the effect of one regional deformation and renders one horizon conceptual. The result of this operation is shown in block 410 .
- Block 412 is a second unfolding operation S2 ⁇ 1 that removes the effect of another regional deformation and renders the second horizon conceptual.
- Block 414 depicts the resulting conceptual model. Instantiating one realization, as shown in block 426 , begins with block 416 , the application of the faulting operator F that instantiates a different type of fault, a reverse fault. The result of this operation is shown in block 418 .
- a first refolding operator deforms both the already realized fault and instantiates a first horizon.
- the result of this operation is shown in block 422 .
- the second refolding operator deforms both the realized fault and horizon and instantiates the second horizon.
- the result of this operation is shown in block 426 , which is the instantiated realization.
- This sequence of operators as shown in block 428 , transformed the base realization 402 to another realization 426 .
- Unfaulting, refaulting, unfolding, and refolding can be performed with different methods depending on the desired degree of accuracy. The methods range from purely geometric methods; to kinematic methods that attempt to preserve distances, areas, and volumes; and to geomechanical methods that model stresses, strains, elasticity, plasticity, failure, etc.
- the modifications include but are not limited to: shift, rotate, scale or deform a fault; change the throws; change the fault type; split one fault into a set of parallel or echelon faults; combine multiple faults into one; or change the topology or interaction between faults.
- the modifications include but are not limited to: shift or deform a horizon; change zone thickness or lateral change rates; change the horizon type; or change the topology or interaction between horizons. Polylines may be shifted, scaled, or deformed.
- a base realization is inputted into the system that is systematically converted to a concept in block 102 by stripping properties and geometries from the base realization using a sequence of unfaulting and unfolding operations.
- at least some aspects of the stripped geometries and properties are retained to aid the selection of a parameter ranges in block 103 .
- an agent e.g., a user and/or a computer program that acts on behalf of a user selects a set of at least one modification from the exemplary set of modifications and specifies parameters for these modifications.
- a fault may be parameterized by its base location, its base shape, a perturbation of its base shape, a scale to increase or reduce its extent, and a throw profile.
- Other parameters may include type and its position within the deformation sequence.
- a horizon may be parameterized by its base location, its base shape, a perturbation of its base shape, a type and its positions within the stratigraphic sequence and deformation sequence.
- a polygon may be parameterized by its base location, its base shape, a perturbation of its base shape, a scale, a shift, or the zone(s) to which it is applied to impart an environment of deposition specified by an agent.
- deformation operators are used to instantiate a realization of the concept in block 104 .
- This may involve an agent specifying a sequence of continuous (e.g., refolding) and discontinuous (e.g., refaulting) deformations.
- the agent parameterizes the individual deformations and applies them to the concept objects of faults, horizons, and polylines.
- the realized faults and horizons may not intersect and truncate each other correctly. For example, horizons may need to be clipped or extended to the faults, while other horizons may need to clipped or extended to other horizons. Further, some faults may need to be clipped or extended to other faults.
- the instantiated framework realization may then be converted into a three-dimensional grid bound by the area of interest. Details of the gridding process are known to those skilled in the art. Examples may include U.S. Patent Application Publication Nos. 2013/218539 and 2012/265510 along with U.S. Pat. No. 7,248,259.
- Block 105 checks the geologic plausibility and technical validity of the instantiated realization. If the framework is found invalid or implausible, this framework realization is removed from the workflow and/or at least flagged. Judicious parameterization and narrow parameter ranges may limit the instantiation of unacceptable realizations. A tradeoff, however, exists between plausibility/validity and variety of realization. Narrow ranges may yield little variety between realizations, but a greater number of plausible or valid ones.
- parameters for the individual deformations are drawn independently for efficiency, but the resulting interaction of deformations may be far-fetched.
- An example of implausible realizations is shown below in FIGS. 5A, 5B and 5C .
- FIGS. 5A, 5B and 5C are diagrams of implausible realizations 500 , 520 and 540 . These realizations 500 , 520 and 540 involve two faults cut through each other multiple times or one fault cuts itself. Realization 500 changes polarity by slowly turning upside down. Realization 520 intersects itself and realization 540 contains two faults 542 and 544 that intersect each other multiple times.
- properties are instantiated, which is shown in block 106 .
- the realization is populated with properties such as porosity, net-to-gross ratio; oil, gas and water saturations, and horizontal and vertical permeabilities.
- the properties can be assigned deterministically, geostatistically, or by simulation; and conditioned or unconditioned with regard to other data, such as well logs or seismic data.
- Methods for instantiating properties in geologic models are known to those skilled in the art. An example is U.S. Pat. No. 7,415,401 to Calvert et al. entitled “Method for constructing 3-D geologic models by combining multiple frequency passbands”.
- properties are conditioned, as shown in block 107 .
- the property or data may be conditioned or at least guided by well markers, well logs, maps, or seismic horizon attributes and seismic volume attributes. All of these conditioning data have geometry, but the present techniques create realizations of geologic models by perturbing, distorting, or modifying geometry.
- Some preferred embodiments may involve modifying the geometry of the conditioning data in an approximately similar manner only to allow for geometric uncertainty in these data caused by data acquisition, data processing, or interpretation.
- some and/or all of the operators in the sequence may be modified or perturbed when applied to conditioning data.
- some additional operators are attached to the operator sequence for the conditioning data to model the geometric uncertainty of the conditioning data.
- the conditioning of the properties may include using undeformation and redeformations on the data to be used in the conditioning.
- the instantiated realization may then be simulated in block 108 and/or analyzed in block 110 to predict a specified quantity.
- Some predictions can be made directly from the subsurface model. Examples include but are not limited to gross rock volume (GRV), the stock tank original oil in place (STOOIP), or the expected ultimate recovery (EUR).
- Other predictions may involve additional financial assumption to calculate cash flow, discounted cash flow (DCF), discounted cash flow rate (DCFR), net present value (NPV), or return on capital employed (ROCE).
- Performing a reservoir simulation in block 108 provides a prediction of water cuts, flow streams, flow capacity, storage capacity, connectivity, or some other performance indicator.
- the analysis of the realization and/or the analysis of its simulation results may be performed, as noted in block 110 .
- multiple realizations are instantiated and simulated and included as part of the analysis.
- blocks 101 to 103 may be performed in an initial stage to convert the base realization to the concept with specified parameter ranges, while blocks 104 to 108 may be repeated through multiple iterations (e.g., multiple times or stages) to generate multiple realizations and/or simulations for analysis.
- the number of repetitions may be controlled by the user or an agent directly by specification or indirectly by selection of a stopping criterion to ensure an appropriate sampling of the parameter space.
- predictions may exhibit transitions between different behaviors where perturbing parameters up to a certain point yields similar results, but perturbing the parameters beyond this point yields a very different result (e.g., different regimes).
- This behavior may be likened to phase transitions in thermodynamic systems where the system can abruptly move to a different state with very different properties.
- the stopping criterion attempts to predict the number of ‘states’ and locate the transitions between the discovered states in the parameter space.
- the different regimes may involve predictions that impact the flow within the reservoir compartments and/or well.
- different regimes may include changes in flow that divide different compartments to adjust the amount of fluid communication between the compartments.
- each base realization corresponds to a different scenario.
- a scenario is an alternative working hypothesis; or in the context of this disclosure, a scenario is an alternative concept.
- the workflow uses at least one concept that in a preferred embodiment is generated from a base realization.
- multiple base realizations may be reduced to multiple concepts that differ from each other. Each of these different concepts may represent a different scenario.
- realizations are compared against each other by use of a metric that is used to group or cluster similar realizations (e.g., Suzuki et al, ‘Dynamic data integration for structural modeling: model screening approach using a distance-based model parameterization’, Computational Geosciences, 105-109, 2008).
- a metric that is used to group or cluster similar realizations (e.g., Suzuki et al, ‘Dynamic data integration for structural modeling: model screening approach using a distance-based model parameterization’, Computational Geosciences, 105-109, 2008).
- Techniques such as multi-dimensional scaling (MDS) may be used to group or cluster realizations and predictions.
- MDS multi-dimensional scaling
- Reservoir simulation can create large amounts of time-dependent results or time-series data.
- these time-series data are presented as contour boxplots (e.g., Sun and Genton, ‘Functional Boxplots’, Journal of Computational and Graphical Statistics, 20(2), 316-334, 2011).
- the inputting of base realization(s) in block 101 may be omitted.
- a user or agent may input or create at least one concept directly, for example, by using a concept editor in block 102 .
- a concept editor provides a mechanism for the creation of conceptual models by specifying at least the number of horizons, faults, and polylines.
- the concept editor may be used to specify certain attributes, such as fault type, horizon type, environment of depositions, and their interactions.
- the concept editor creates objects for the specified entities directly in a geologic modeling software package where they can be operated on with a sequence of deformation operators F and S.
- the concept editor creates objects for the specified entities either in memory, a file system, or in the cloud from where they can be imported by a geologic modeling software package to be operated on with a sequence of deformation operators F and S.
- FIGS. 6A, 6B and 6C are diagrams 600 , 620 and 640 of unfaulting during concept creation.
- the concept creation may be performed in block 102 of FIG. 1 .
- FIG. 6A is a diagram 600 having a horizon 601 that is bisected by normal fault 602 resulting in the foot-wall truncation 604 and hanging-wall truncation 603 .
- the truncations 604 and 603 are preferably represented as polylines forming a cutoff polygon for horizon 601 against fault 602 .
- a reference, such as reference line 605 is created from the foot-wall and hanging-wall polylines.
- One method for creating the reference is simply to average the depths or two-way travel times of the cutoff polygon.
- a local or floating reference is created for every polygon point.
- the laterally nearest point (e.g., neglecting vertical offset) of the hanging-wall polyline is determined and the local reference for the specified point is determined by averaging its depth with the depth of the nearest point on the hanging-wall polyline.
- the laterally nearest point (e.g., neglecting vertical offset) of the foot-wall polyline is determined and the local reference for the specified point is determined by averaging its depth with the depth of the nearest point on the foot-wall polyline.
- the dynamic-time-warping (DTW) algorithm may be an efficient method to determine corresponding points on foot-wall and hanging-wall polylines.
- the residual polygon consisting of residual polylines 603 ′ and 604 ′ is created by subtraction of the reference 605 from the cutoff polygon formed by polylines 603 and 604 .
- FIG. 6B is a diagram 620 having an unfaulting map or correction map 621 that is formed from the residual polylines 623 (corresponding to 603 ′ of FIG. 6A ) and 624 (corresponding to 604 ′).
- the map is formed by extrapolation from the residual polylines 623 and 624 .
- the extrapolation converges toward the level of zero, as shown by reference line 625 (corresponding to 605 ), at distance from the specified residual polylines.
- the extrapolator may involve regularization or other forms of extrapolation constraints. Minimal curvature may be a preferred regularization method.
- FIG. 6C is a diagram 640 having an unfaulted horizon 641 that is formed by subtracting the map 643 (corresponding 621 of FIG. 6B ) from the horizon 642 (corresponding to 601 of FIG. 6A ).
- the reference or level of zero is shown by reference line 645 .
- the effect of fault 644 is removed from horizon 641 .
- Fault 644 has been reduced to a concept as indicated by the dashes.
- the unfaulted horizon may contain gaps or artifacts that are preferably removed by filtering and interpolation.
- the base realization contains more than one fault
- removal of the first fault changes the horizon(s) and thus the cutoff polygon(s) for the remaining fault(s).
- the original cutoff polygons may be corrected by subtraction of the correction map. For example, if the base realization contains three faults, then removal of the first fault triggers the correction of the cutoff polygons around the second and third faults. Subsequent removal of the second fault triggers another correction of the cutoff polygon around the third fault. Subsequent removal of the third fault does not trigger any further corrections because the all faults are reduced to concepts.
- Realizing a fault has two aspects: (i) the specification of all its parameters and (ii) the redeformation of the other objects in the model. For computational efficiency, it may be advantageous to create an explicit fault object between the first and the second aspects.
- FIGS. 7A, 7B and 7C are diagrams 700 , 720 and 740 of refaulting during the instantiation of a realization.
- the instance realization may be performed in block 104 of FIG. 1 .
- FIG. 7A is a diagram 700 having a conceptual reverse fault 702 indicated by dashes that bisects horizon 701 .
- the process begins with parameterizing the fault 702 by specification of the fault geometry (e.g., location, orientation, shape, size, etc.). Some of these geometry parameters may be prescribed by the parameter ranges specified in block 103 of FIG. 1 , while others may be drawn at random from statistical distribution functions specified in block 103 of FIG. 1 .
- combinations of parameters may be selected by systematic sampling of the parameter ranges.
- a fault throw is also specified, for example at the intersection 708 of the realized fault with the horizon or preferably for every location on the fault by designating fault throw as a property attached to the fault.
- the fault-horizon intersection 708 is also used to define the local reference depth 705 .
- Fault type, reference, and throw are used to define the cutoff polygon consisting of foot-wall polyline 704 and hanging-wall polyline 703 .
- the foot-wall polyline 704 is determined by shifting the fault-horizon intersection 708 downwards along the realized fault 702 by half the local throw, while the hanging-wall polyline 703 is determined by shifting the fault-horizon intersection 708 upwards along the realized fault 702 by half the local throw.
- foot-wall and hanging-wall polylines are swapped: the foot-wall polyline 704 is determined by shifting the fault-horizon intersection 708 upwards along the realized fault 702 by half the local throw, while the hanging-wall polyline 703 is determined by shifting the fault-horizon intersection 708 downwards along the realized fault 702 by half the local throw.
- the polylines 703 and 704 are found by vertical shifting only of 708 , neglecting any lateral component introduced by shifting along the fault surface itself.
- the residual polylines 703 ′ and 704 ′ are determined from the polylines 703 and 704 by subtraction of the reference 705 .
- FIG. 7B is a diagram 720 having a refaulting map or correction map 721 (consisting of 721 ′ and 721 ′′) that is formed from the residual polylines 723 (corresponding to 704 ′ of FIG. 7A ) and 724 (corresponding to 703 ′ of FIG. 7A ).
- the map is formed by extrapolation from the residual polylines 723 and 724 .
- the extrapolation converges toward the level of zero 725 at distance from the specified residual polylines.
- the extrapolator may require regularization or another form of extrapolation constraint. Minimal curvature is a preferred regularization.
- a flat, single-valued horizon remains single valued; and the correction map 721 can be extrapolated from 723 and 724 directly without invoking 721 ′ and 721 ′′.
- a flat, single-valued horizon will become multi valued. In the region between the foot-wall cutoff and the hanging-wall cutoff, the horizon will be duplicated and overlapping itself.
- the refaulting map is multi valued between the foot-wall and the hanging-wall cutoffs shifting the meaning of map from a two-dimensional depiction of residual elevation toward a mathematical operator or transform. It may be advantageous to divide the multi-valued map 721 into the single-valued maps 721 ′ and 721 ′′.
- FIG. 7C is a diagram 740 having a refaulted horizon 741 consisting of 741 ′ and 741 ′′ that is formed by adding the map 742 consisting of 742 ′ (corresponding to 721 ′ of FIG. 7B ) and 742 ′′ (corresponding 721 ′′ of FIG. 7B ) to the horizon 744 (corresponding to 701 of FIG. 7A ), while providing multivaluedness or overlap in horizon 741 by use of an appropriate representation.
- the effect of fault 743 is thus imparted onto horizon 741 .
- Fault 743 has been realized from a concept as indicated by the solid line.
- the reference or level of zero may be shown by reference line 745 .
- the unfaulted (conceptual) horizon piece inside the cutoff polygon is used twice as it gets added both to 742 ′ and to 742 ′′.
- the unfaulted (conceptual) horizon piece inside the cutoff polygon would not be used at all.
- the refaulted horizon 741 may contain gaps or artifacts that are preferably removed by filtering and interpolation.
- a process such as disclosed in U.S. Pat. No.
- 7,756,694 is used to clean up the fault-horizon intersection by extrapolation of the horizon to the fault, cutback and truncation of the horizon by the fault, and creation of a watertight intersection.
- the process may also be used to clean up fault-fault or horizon-horizon intersections and to create cutoff polygons.
- realization of the first fault changes the horizon(s) and thus the intersections with the remaining fault(s).
- Either the original cutoff polygons (intersections between fault and horizon, e.g., fault-horizon intersection 708 ) are corrected by addition of the correction map, or preferably, cutoffs at the remaining faults are recreated.
- realization of the first fault triggers the correction of the cutoff polygons around the second and third faults.
- realization of the second fault triggers another correction of the cutoff polygon around the third fault.
- Subsequent realization of the third fault does not trigger any further corrections because the all faults have been realized.
- the map (e.g., map 742 ) is added only to the horizon (e.g., horizon 744 ), while in others the map is also added to some or all of the faults that are already realized to preserve their relative positions during the refaulting operation.
- FIG. 3 is a diagram 800 of the process from a base framework realization to an instantiated framework realization in accordance with an exemplary embodiment of the present techniques.
- Diagram 800 presents an exemplary application of the workflow from a base framework realization, such as base framework realization 810 , to an instantiated framework realization, such as such as instantiated framework realization 860 .
- the process begins with the base framework realization 810 .
- the base framework realization 810 includes three faults 811 , 812 , and 813 and one horizon 814 .
- the three faults are converted to conceptual faults by removing their geometry and healing their effects on the horizon 814 . This may be performed by applying the unfaulting operators, such as first unfaulting operator F1 ⁇ 1 for the first fault 811 , second unfaulting operator F2 ⁇ 1 for the second fault 812 , and third unfaulting operator F3 ⁇ 1 for the third fault 813 .
- These different unfaulting operators may be combined:
- the result is the intermediary model 820 that contains three conceptual faults 821 , 822 , and 823 and one horizon 824 .
- the horizon 824 is still realized, but healed.
- the horizon 824 does not exhibit any spatial discontinuities.
- Horizon 824 is continuous, while the base horizon 814 contained discontinuities at the fault locations.
- model 830 which is the conceptual model having three conceptual faults 821 , 822 , and 823 and the conceptual horizon 834 .
- the conceptual horizon 834 is reinstantiated creating the realized horizon 844 in the reinstantiated model 840 . This may be performed by applying the folding operator S1. Without any realized faults, the realized horizon 844 is continuous, but clearly different from horizon 824 .
- the conceptual fault 822 is reinstantiated creating the realized fault 852 in reinstantiated model 850 .
- This may be performed by applying the folding operator F2.
- both faults 821 and 823 were randomly determined to remain concepts and are not reinstantiated.
- the realization of fault 852 also introduced throws which that are applied to the continuous horizon 844 creating the faulted, discontinuous horizon 854 .
- the final framework realization 860 contains fault 852 and horizon 854 .
- Multiple framework realizations may be instantiated from the concept 830 by using different parameterizations for the conceptual faults 821 , 822 , and 823 and the conceptual horizon 834 . Creating multiple realizations may be useful to lessen uncertainty in the analysis of the realizations with regard to a specified problem, question, or decision.
- the flow chart of FIG. 1 may include various variations.
- the concept may be created in block 102 .
- an agent selects bounds for deformations and other model parameters. Then, parameters are selected from these bounds, and a realization of the concept is instantiated in block 104 .
- the concept is generated in block 102 by systematic removal of deformations from a base realization.
- the realization is checked for technical validity and/or geologic plausibility in block 105 because some combinations of parameters may result in realizations that are technically invalid or geologically implausible. Realizations that fail this test are either fixed or discarded outright. Then, the realizations are populated with properties in block 106 , simulated in block 108 , and analyzed in block 110 . The analysis results are summarized to facilitate business decisions and operations to produce hydrocarbons.
- FIG. 9 is a block diagram of a computer system 900 that may be used to perform any of the methods disclosed herein.
- a central processing unit (CPU) 902 is coupled to system bus 904 .
- the CPU 902 may be any general-purpose CPU, although other types of architectures of CPU 902 (or other components of exemplary system 900 ) may be used as long as CPU 902 (and other components of system 900 ) supports the inventive operations as described herein.
- the CPU 902 may execute the various logical instructions according to disclosed aspects and methodologies. For example, the CPU 902 may execute machine-level instructions for performing processing according to aspects and methodologies disclosed herein.
- the computer system 900 may also include computer components such as a random access memory (RAM) 906 , which may be SRAM, DRAM, SDRAM, or the like.
- the computer system 900 may also include read-only memory (ROM) 908 , which may be PROM, EPROM, EEPROM, or the like.
- RAM 906 and ROM 908 hold user and system data and programs, as is known in the art.
- the computer system 900 may also include an input/output (I/O) adapter 910 , a communications adapter 922 , a user interface adapter 924 , and a display adapter 918 .
- the I/O adapter 910 , the user interface adapter 924 , and/or communications adapter 922 may, in certain aspects and techniques, enable a user to interact with computer system 900 to input information.
- the I/O adapter 910 preferably connects a storage device(s) 912 , such as one or more of hard drive, compact disc (CD) drive, floppy disk drive, tape drive, etc. to computer system 900 .
- the storage device(s) may be used when RAM 906 is insufficient for the memory requirements associated with storing data for operations of embodiments of the present techniques.
- the data storage of the computer system 900 may be used for storing information and/or other data used or generated as disclosed herein.
- the communications adapter 922 may couple the computer system 900 to a network (not shown), which may enable information to be input to and/or output from system 900 via the network (for example, a wide-area network, a local-area network, a wireless network, any combination of the foregoing).
- User interface adapter 924 couples user input devices, such as a keyboard 928 , a pointing device 926 , and the like, to computer system 900 .
- the display adapter 918 is driven by the CPU 902 to control, through a display driver 916 , the display on a display device 920 .
- Information and/or representations of one or more 2D canvases and one or more 3D windows may be displayed, according to disclosed aspects and methodologies.
- system 900 may be varied as desired.
- any suitable processor-based device may be used, including without limitation personal computers, laptop computers, computer workstations, and multi-processor servers.
- embodiments may be implemented on application specific integrated circuits (ASICs) or very large scale integrated (VLSI) circuits.
- ASICs application specific integrated circuits
- VLSI very large scale integrated circuits
- the method may be implemented in machine-readable logic, set of instructions or code that, when executed, performs a method to analyzing uncertainty of subsurface formations.
- the code may be used or executed with a computing system such as computing system 900 .
- the computer system may be utilized to store the set of instructions that are utilized to manage the data and other aspects of the present techniques.
- a computer system 900 may be used to analyze uncertainty of subsurface formations for production or exploration operations.
- the computer system may include a processor; memory in communication with the processor; and a set of instructions stored in memory and accessible by the processor.
- the set of instructions when executed by the processor, are configured to: create a conceptual subsurface model, wherein the conceptual subsurface model is associated with a subsurface formation and comprises a plurality of objects; select parameter ranges for each of the plurality of objects and interactions between two or more of the plurality of objects; instantiate realizations for the plurality of objects based on the selected parameter ranges; and combine instantiated realizations of these objects into a reservoir model.
- the set of instructions are configured to create the conceptual subsurface model may be further configured to automatically create the conceptual subsurface model from an obtained base realization; may be further configured to undo one or more faults and folds in a sequential order; may be further configured estimate parameter ranges based on the undoing of one or more of faults and folds; and may be further configured to unfault an inputted base realization based on fault-horizon cutoff polygons to create the conceptual subsurface model. Further, the set of instructions may be configured to refault from the conceptual subsurface model based on fault-horizon cutoff polygons.
- the computer system may include other instructions to enhance efficiency of the operation of the present techniques.
- the set of instructions may be configured to analyze each of the instantiated realizations for geologic plausibility and, if one or more of the instantiated realizations are determined to be geologically implausible, discard the one or more of the instantiated realizations that are geologically implausible.
- the set of instructions may be configured to analyze each of the instantiated realizations for technical consistency and, if one or more of the instantiated realizations are determined to be technically inconsistency, discard the one or more of the instantiated realizations that are technically inconsistency.
- the computer system may include other enhancements.
- the set of instructions may be configured to instantiate properties into each of the instantiated realizations, wherein the properties comprise one or more of porosity, permeability and oil saturation; may be configured to condition the instantiating properties by perturbing, distorting, or modifying the geometry; and/or may be configured to condition the instantiating properties by applying a sequence of undeformations and redeformations that used to create the instantiated realizations.
- the set of instructions may be configured to simulate the instantiated realizations; may be configured to simulate using a simulation proxy method; may be configured to simulate proxy method using a connectivity measure as a simulation proxy for each of the instantiated realizations; may be configured to compute the connectivity measure based on graph based centrality measure; may be configured to rank the plurality of instantiated realizations in order of the respective centralization measures; and/or may be configured to simulate the instantiated realizations to create a set of simulations that are analyzed to affect a decision for production operations.
- the set of instructions configured to create the conceptual subsurface model may be further configured to create two or more instantiated realizations from the concept conceptual subsurface model; may be further configured to create two or more conceptual subsurface models to generate two or more scenarios; and/or may be further configured to create one or more conceptual subsurface models that are each based on different base realizations and are created to generate two or more scenarios.
- simulation may be approximated by a simulation proxy.
- a preferred simulation proxy is based on graph-based centrality.
- the centrality measure which is described in U.S. patent application Ser. No. 14/272,581, which is incorporated by reference, may include one or more of degree, betweenness, closeness, and eigenvector.
- a connectivity matrix expresses how well two neighboring grid cells are connected (transmissibility) or how similar a specified property is.
- a connection is weighted; while in the second case, each grid cell is associated with a label or index i and an attribute or property value vi.
- the two cases are not mutually exclusive: one definition of connection weight is the magnitude of their property or attribute difference. Another preferred definition of connection weight is their property or attribute average.
- an off-diagonal element of the connectivity matrix Cij for two neighboring grid cells i and j (where i ⁇ j) is set to ⁇ 1/2(vi+vj).
- a diagonal element Cii of the connectivity matrix is set to ⁇ 1/2 ⁇ ij (v i +v j ) where ⁇ ij is one when grid cells i and j are neighbors and zero when grid cells i and j not neighbors.
- the diagonal elements of the connectivity matrix are set to zero, effectively removing a self-interaction or self-connectivity.
- specified eigenvectors of the connectivity matrix are used to compute a connectivity measure for the grid cell.
- the first component of the specified eigenvectors defines the location of the first grid cell in a vector space.
- the second component of the specified eigenvectors defines the location of the second grid cell in said vector space, and so on for the remaining components and grid cells.
- the shortest distance to any other grid cell in said space defines a measure of connectivity indicating how connected the specified grid cell is to all others. Iterating this process over substantially all grid cells provides a connectivity measure for substantially every grid cell, resulting in a connectivity attribute.
- distance function Details of the distance function are irrelevant. Different distance functions result in different connectivity measures. Any metric or any generalized metric associated with said vector space results in a connectivity measure.
- distances may be computed directly from the connectivity matrix using either an iterative or algebraic process.
- the connectivity measure ci is computed iteratively as
- d is a small damping coefficient
- Mij 1/Cij if Cij ⁇ 0 and zero otherwise
- N refers to the number of grid cells
- 1 is a vector of dimension N containing only ones.
- An initial value for c may be 1/N.
- the connectivity matrix is normalized prior to the direct estimation of connectivity measures, for example by scaling each row sum, each columns sum, or each row sum and each column sum of the connectivity matrix C to one.
- centrality of a vertex measures its relative importance within a graph. Examples include how influential a person is within a social network, how well-used a road is within an urban network, or how well connected the grid cells are within their geobodies or connectivity structures. There are four main measures of centrality: degree, betweenness, closeness, and eigenvector.
- the connectivity measures disclosed with this invention are examples of eigenvector-based centrality measures.
- Degree centrality refers to the number of connections for a specified node, potentially weighted by the attribute value.
- degree centralities or degree-based connectivity measure may be computed by row sums, column sums, or row-column sums, preferably excluding elements on the matrix diagonals from the sum.
- a distance metric between any two grid cells belonging to this graph that is defined by the length of the shortest path between the two specified grid cells.
- the length of a path is defined by the number of connections linking the two specified grid cells, or in the attributed case, by the sum of the attributes along a path linking the specified grid cells.
- the farness of any grid cell is defined by the sum of its distances to all other grid cells of the graph. Closeness centrality is defined as the inverse of farness. The more central a grid cell is, the lower its total distances to all other grid cells. Closeness centrality can be viewed as a measure of how long it takes to spread information sequentially from a grid cell to all other grid cells belonging to the same graph.
- the grid cells with the largest closeness centralities or the largest closeness-based connectivity measures are the grid cells that provide fast drainage of a contiguous group of grid cells from their fluids.
- Extensions of closeness centrality account not only for the shortest path length but also for the number of paths.
- Betweenness centrality quantifies the number of times a grid cell acts as a bridge along the shortest path between any two grid cells of a subsurface model. It may be advantageous to scrutinize grid cells with high betweenness centrality because a small perturbation to the connectivity structure, permeability or transmissibility might dramatically alter the shortest paths and their spatial distributions.
- Eigenvector centrality is a measure of the influence of a grid cell in the connected graph of the subsurface model. Eigenvector centrality assigns a relative score to all grid cells based on the principle that connections from a specified grid cell to high-scoring grid cells contribute more to the score of the specified grid cell than connections to low-scoring grid cells.
- the dominant eigenvector associated with the largest eigenvalue is preferably obtained by an iterative process.
- a centralization measure is computed for a group of contiguous grid cells that have been attributed with a specified connectivity measure.
- Centralization for the specified group of grid cells measures how central its most central grid cell is in relation to all of its other grid cells, for example by computation of ⁇ c max ⁇ c i .
- this quantity is normalized by the number of grid cells or the theoretically largest sum of centrality differences for a graph of similar size. It may be advantageous to estimate the theoretically largest sum of centrality differences for a graph of similar size by constructing a compact group of grid cells with the same number of grid cells and maximal connectivity, for example in the shape of a ball. In the attributed case, every grid cell or connection of this ideal group is attributed with a maximal value in accordance to the specified attribute.
- groups of contiguous grid cells are ranked in order of their centralization measures.
- the group of contiguous grid cells is formed by thresholding, by definition of a spatial bounding box, or by any other method.
- a connectivity measure is assigned to groups of contiguous grid cells of the reservoir model.
- the connectivity measure serves as a proxy to a reservoir simulation or reservoir performance analysis. Proxy simulations for performance prediction are well known to practitioners of the art. Examples of such proxy simulations may be European Patent No. 1,994,488 to Li et al entitled “Method for Quantifying Reservoir Connectivity Using Fluid Travel Times”, U.S. Pat. No. 8,437,997 to Meurer et al entitled ‘Dynamic Connectivity Analysis’, U.S. Pat. No.
- proxies The purpose of these proxies is the analysis of different reservoir development or production scenarios to examine the connectivity between the oil-bearing reservoir and the producer wells or the connectivity between water-injection wells and hydrocarbon-production wells.
- the novel connectivity measures disclosed in this publication are independent of sources and targets. No well locations need to be specified. Grid cells do not need to be separated into mutually exclusive sources, sinks, and conductors. Instead, each grid cell is compared to all others. Each grid cell acts simultaneously as source, sink, and conductor.
- the disclosed connectivity measures allow examination of the model for highly connected regions, for disconnected compartments, for barriers, and regions where small even perturbations of connectivity and attributes or properties (porosity, permeability, or transmissibility) will change long-distance connectivity by disconnecting one region or compartment into multiple ones or connecting multiple regions or compartments into one, thus warranting additional scrutiny to analyze these sensitive regions.
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Description
- This application claims the benefit of U.S. Provisional Patent Application 62/057,797, filed Sep. 30, 2014, entitled METHOD AND SYSTEM FOR ANALYZING THE UNCERTAINTY OF SUBSURFACE MODEL, the entirety of which is incorporated by reference herein.
- This invention relates generally to the field of hydrocarbon exploration and extraction. Specifically, the invention is a method to examine uncertainty and risk associated with the development of a hydrocarbon resource by rapidly generating and analyzing variations of subsurface models realized from different scenarios.
- This section is intended to introduce various aspects of the art, which may be associated with exemplary embodiments of the present disclosure. This discussion is believed to assist in providing a framework to facilitate a better understanding of particular aspects of the disclosed methodologies and techniques. Accordingly, it should be understood that this section should be read in this light, and not necessarily as admissions of prior art.
- Typically, a geologic model is formed that includes various static properties. From the geologic model, a reservoir model is created to model dynamic properties. For example, the reservoir model is formed from geological horizons and faults. The reservoir model includes a framework that establishes the geometrical foundation for the three-dimensional grid and provides some of the boundaries for facies and petrophysical models that describe rock properties and contained fluids. The resulting reservoir model forms the basis for volumetric computations, reservoir simulations, facilities planning computations and well planning computations. While seismic and well data provide information regarding the reservoir model, considerable uncertainty may remain regarding the reservoir model.
- The effect of uncertainty is often examined through various means. For example, the uncertainty may be examined by perturbing uncertain aspects or features of the reservoir model, recomputing the quantity of interest, and examining sensitivity of the quantity of interest with regard to the uncertain aspects. The problem with framework uncertainty related to the geometric foundation of the three-dimensional grid is that the steps between framework generation, definition of the grid, and computation of the quantity of interest are computationally and labor intensive, often requiring user input.
- Some conventional methods only perturb the depth of different model objects, such as faults and horizons. The depth perturbation may be spatially variable, for example allowing the flanks of an anticline to be pushed down or pulled up, while leaving the crest unperturbed. The resulting flexing of faults, horizons and other model objects occurs in the vertical direction, however, and the modeling grid is flexed simultaneously. That is, the geometry of the grid changes, but not its structure.
- Providing vertical and lateral movement of the model objects typically requires grid regeneration. Some conventional methods attempt to correct the existing grid. Other conventional methods move the model objects vertically and laterally and then adjust the intersections. The adjustment has two aspects, however, the removal of the original intersections and the implementation of the new intersections. For example, when a fault is shifted, a portion of a bisected horizon moves from the foot-wall side to the hanging-wall side or vice versa. The moving horizon piece has a previous displacement that should be undone. Typically, the undoing of previous displacements appears to be performed by deletion of horizons near faults and extrapolation towards the shifted fault. See, e.g., Roe et al., ‘Flexible Simulation Of Faults’, SPE 134912, 2010; Cherpeau et al., ‘Stochastic Simulations of Fault Networks in 3D Structural Modeling’, Comptes Rendus Geoscience, 342, 687-694, 2010; Suzuki et al., ‘Dynamic Data Integration For Structural Modeling: Model Screening Approach Using Distance-Based Model Parameterization’, Computational Geosciences, 12, 105-119, 2008; Holden et al., ‘Stochastic Structural Modeling’, 35(8), 899-913, 2003; and Thore et al., ‘Structural Uncertainties: Determination, Management And Applications’, Geophysics, 67(3), 840-852, 2002.
- As a result, an enhancement to exploration and reservoir delineation techniques is needed to identify and recover hydrocarbons in light of imprecise data. The present techniques provide a streamlined method for generation of a perturbed framework and thus a perturbed reservoir model. Further, the enhancements may provide a method for systematic removal of the effects of faults and folds to provide numerous realizations of the reservoir model in an efficient manner. The enhancements may provide a method to explore fault connectivity or check for geologic plausibility and technical validity which removes problematic model perturbation.
- In one embodiment, a method is described. The method includes analyzing uncertainty of subsurface formations and includes: creating a conceptual subsurface model, wherein the conceptual subsurface model is associated with a subsurface formation and comprises a plurality of objects; selecting parameter ranges for each of the plurality of objects and interactions between two or more of the plurality of objects; instantiating realizations for the plurality of objects based on the selected parameter ranges; and combining instantiated realizations of the plurality of objects into a reservoir model.
- In yet another embodiment, a computer system for analyzing uncertainty of subsurface formations for production or exploration operations is described. The computer system may include a processor; memory in communication with the processor; and a set of instructions stored in memory and accessible by the processor. The set of instructions, when executed by the processor, are configured to: create a conceptual subsurface model, wherein the conceptual subsurface model is associated with a subsurface formation and comprises a plurality of objects; select parameter ranges for each of the plurality of objects and interactions between two or more of the plurality of objects; instantiate realizations for the plurality of objects based on the selected parameter ranges; and combine instantiated realizations of these objects into a reservoir model.
- The foregoing and other advantages of the present disclosure may become apparent upon reviewing the following detailed description and drawings of non-limiting examples of embodiments.
-
FIG. 1 is a flow chart for performing hydrocarbon exploration in accordance with an exemplary embodiment of the present techniques. -
FIG. 2 is a diagram of a realization of a geologic model. -
FIG. 3 is a diagram of a concept of the realization ofFIG. 2 . -
FIG. 4 is a diagram of a schematic application of the undeforming and redeforming a framework in accordance with an exemplary embodiment of the present techniques. -
FIGS. 5A, 5B and 5C are diagrams of implausible realizations. -
FIGS. 6A, 6B and 6C are diagrams of unfaulting during concept creation. -
FIGS. 7A, 7B and 7C are diagrams of refaulting during instantiation of a realization. -
FIG. 8 is a diagram of the process from a base framework realization to an instantiated framework realization in accordance with an exemplary embodiment of the present technique. -
FIG. 9 is a block diagram of a computer system that may be used to perform any of the methods disclosed herein. - In the following detailed description section, the specific embodiments of the present disclosure are described in connection with preferred embodiments. However, to the extent that the following description is specific to a particular embodiment or a particular use of the present disclosure, this is intended to be for exemplary purposes only and simply provides a description of the exemplary embodiments. Accordingly, the disclosure is not limited to the specific embodiments described below, but rather, it includes all alternatives, modifications, and equivalents falling within the true spirit and scope of the appended claims.
- Various terms as used herein are defined below. To the extent a term used in a claim is not defined below, it should be given the broadest definition persons in the pertinent art have given that term as reflected in at least one printed publication or issued patent.
- “Subsurface model”, as used herein, is a reservoir model or a geologic model.
- “Geologic model”, as used herein, is three-dimensional model of the subsurface having static properties and includes faults, horizons, facies, lithology and properties such as porosity, permeability, or the proportion of sand and shale.
- “Reservoir model”, as used herein, is a three-dimensional model of the subsurface that in addition to static properties such as porosity and permeability also has dynamic properties that vary over the timescale of resource extraction such as fluid composition, pressure, and relative permeability.
- “Framework”, as used herein, is a geologic model formed from faults, horizons, model boundaries, and facies boundaries, e.g., a geologic model containing only surfaces and polylines.
- “Concept”, as used herein, is a model containing faults, horizons, and facies boundaries without geometry or location. Aspects of the geologic or reservoir models related to geometry and location are suppressed. For example, a concept is often conveyed as a sketch of only the objects deemed most pertinent, omitting anything else.
- “Scenario”, as used herein, is a concept or partial subsurface model in combination with select parameters and their ranges used to build realizations of subsurface models by deterministically or stochastically varying these parameters. Examples may be normal faults versus reverse faults, different time-depth conversions or environments of deposition, or high net-to-gross versus low-net-to-gross regions.
- “Realization”, as used herein, is a subsurface model (e.g., geologic model) created from a concept or scenario by assigning geometry and location to faults, horizons, and boundaries; and values to properties which may be utilized for computations and quantitative queries.
- “Instantiate”, as used herein, is the process of transforming a (qualitative) concept or object thereof to a (quantitative) realization or object thereof. Instantiation may be performed by interpolation or extrapolation of measurements or by application of a stochastic process.
- “Simulate”, as used herein, is the process of making a prediction related to the resource extraction based on the reservoir model. A simulation is typically performed by execution of a reservoir-simulator computer program on a processor, which computes composition, pressure, or movement fluid as function of time and space for a specified scenario of injection and production wells by solving a set of reservoir fluid flow equations.
- “Unfaulting operation”, as used herein, is the process of undoing the effects of a fault on the other objects of the model. A horizon typically exhibits a discontinuity where intersected by a fault. When this fault is turned into a conceptual fault, its location and geometry are removed. Remaining model objects, however, still exhibit the discontinuities caused by the original fault. The unfaulting operation removes these discontinuities, returning the remaining objects to a state unaffected by this fault, e.g., to a state where the fault never existed.
- “Refaulting operation”, as used herein, is the process of affecting the other model objects of a model when a fault is realized. The typical effect is the generation of discontinuities where the fault and existing objects intersect.
- “Unfolding operation”, as used herein, is process of at least partially removing location and geometry from a horizon transforming a horizon into a conceptual horizon. An unfolding operation removes the effects of a folding event (continuous deformation) from the model.
- “Refolding operation”, as used herein, is the process of applying a continuous deformation to the objects, either by assignment or perturbation of geometry and location to horizons and specified other objects.
- To begin, a subsurface model, which may include a reservoir model or geologic model, is a computerized representation of a subsurface region based on geophysical and geological observations made on and below the surface of the Earth. It is the numerical equivalent of a three-dimensional geological map complemented by a description of physical quantities in the domain of interest. Subsurface models, which are reservoir models, are often used as inputs to reservoir simulation programs that predict the behavior of rocks and fluids contained therein under various scenarios of hydrocarbon recovery. When producing an actual hydrocarbon reservoir, miscalculations or mistakes can be costly. Using subsurface models in simulations provides a mechanism to identify which recovery options offer the safest and most economic, efficient, and effective development plans for a particular reservoir.
- Construction of a subsurface model is typically a multistep process. First, a structural model or structural framework is created from surfaces that include faults, horizons, and if necessary, additional surfaces that bound the area of interest for the geologic model. The different surfaces define closed volumes often called zones. Second, each zone is meshed or partitioned into small cells defined by a three-dimensional grid. Lastly, properties are assigned to surfaces (e.g., transmissibility) and individual cells (e.g., rock type, porosity, permeability, or oil saturation).
- The assignment of cell properties is often a two-step process where each cell is first assigned a rock type, and then each rock type is assigned spatially-correlated reservoir properties and/or fluid properties. Each cell in the model is assigned a rock type. For example, in a coastal clastic environment, the cells may be beach sand, high water energy marine upper shoreface sand, intermediate water energy marine lower shoreface sand, and deeper low energy marine silt and shale. The distribution of these rock types within the model may be controlled by several methods, including map boundary polygons, rock type probability maps, or statistically emplaced based on concepts. Where available, rock type assignment may be conditioned to well data.
- Reservoir quality parameters typically include porosity and permeability, but may include measures of clay content, cementation factors, and other factors that affect the storage and deliverability of fluids contained in the pores of those rocks. Geostatistical techniques are typically used to populate the cells with porosity and permeability values that are appropriate for the rock type of each cell. Rock pores are saturated with groundwater, oil or gas. Fluid saturations may be assigned to the different cells to indicate which fraction of their pore space is filled with the specified fluids. Fluid saturations and other fluid properties may be assigned deterministically or geostatistically.
- Geostatistics is useful in modeling to interpolate observed data and to superimpose an expected degree of variability. As an example, kriging, which uses the spatial correlation among data and intends to construct the interpolation via semi-variograms, may be used. To reproduce more realistic spatial variability and help assessing spatial uncertainty between data, geostatistical simulation is often used, for example based on variograms, training images, or parametric geological objects. Perturbing surface properties or cell properties, such as rock type, reservoir properties or fluid properties, is a conventional process, which may utilize deterministic or geostatistical methods to assign them. The assignment may include choosing a different variogram for kriging or a different seed for geostatistical simulation.
- For the purpose of this disclosure, a realization is instantiated from a concept. The difference between a concept and a realization is a complete geometry. In a realization, objects, such as points, polylines, polygons, horizons, faults, or compartments, have locations, relative positions with regard to each other, shapes, or sizes. The topology of the objects (e.g., their interactions) of the realization is defined. In a realization, a property attached to an object has values at essentially every location of the object. In summary, a realization does not contain free parameters anymore. On the other hand, a concept contains free parameters relating to topology, geometry, and properties. A concept does not have geometry associated with its objects. At least some of the points, polylines, polygons, horizons, faults, or compartments do not have their locations, relative positions with regard to each other, shapes, or sizes defined. The topology of its objects can be completely specified, partially specified, or unspecified. For example, the order in which different horizons and faults terminate may be undefined. In some embodiments, the nature or interpretation of an object may be undefined. In a realization, a horizon is (either implicitly or explicitly) typed conformable, basal, topal, erosional, or discontinuous while in a concept, the horizon type can but does not need to be specified. In a realization, a fault is typed normal, reverse, or strike-slip, while in a concept, the fault type can but does not need to be typed. In a concept, the cardinality of an object can be undefined. In some embodiments, a single fault in the concept may be realized as multiple faults, for example, by realizing a fault as a set of relay faults or as a set of parallel faults.
- One of the largest uncertainties in reservoir modeling relates to the framework formed by faults and certain horizons because this framework controls volumetrics and connectivity. Framework uncertainty is caused by uncertainty in seismic migration, time-depth (or depth-true depth) conversion, structural interpretation, fault positions, well picks, horizon correlation and interpretation, and layer thicknesses. Unfortunately, framework changes tend to create artifacts that have to be addressed. Accordingly, the present techniques provide a method and system for modifying an existing watertight (e.g., sealed) framework in such a manner that the resulting framework is yet again watertight without gaps or overlaps between model cells or reservoir compartments.
- The present techniques involve a workflow that may be used to analysis uncertainty and risk associated with the development of hydrocarbon resources by rapidly generating and analyzing variations of subsurface models realized from scenarios. Under the present techniques, structural artifacts may be removed once by systematic removal of the effects of faults and folds. In this manner, numerous realizations of the subsurface model can then be generated by starting with an initial starting subsurface model and the corresponding concept and then, iteratively computing the effects of faults and folds on the revised starting subsurface model.
- The present techniques may also include the ability to omit a fault or fold from the perturbed model or to introduce additional faults and folds into the perturbed model. This capability allows the exploration of fault connectivity, for example the substitution of one long contiguous fault by a set of fault relays. Further still, the method may also include a performance of a check for geologic plausibility and technical validity which removes problematic model perturbations from the workflow or at least from the set of simulation results used for the final analysis.
- In certain embodiments, the workflow may include various steps. For example, as a first step, each horizon is unfaulted one fault at a time. Each unfaulting operation removes a tear in a horizon, bringing its edges back together. A fault causes a discontinuity in a given horizon. Unfaulting removes this discontinuity by adjusting the depths of the horizon near fault and propagating or extrapolating these adjustments away from the fault along the horizon. As a second step, the overarching remaining structure of now fault-free horizons is removed. This unfolding is performed by replacing the horizon with an approximately planar one while recording the necessary vertical depth adjustments. As a third step, a different overarching structure is imposed on the horizons. This refolding is performed by replacing the approximately planar horizons with differently shaped ones, preferably guided by the previously recorded adjustments. In the fourth step, the horizons are refaulted. The refaulting is performed by moving the horizon on the different sides of the fault to the desired relative location and propagating these adjustments away from the fault along the horizon. Unfaulting (F−1), refaulting (F), unfolding (S−1), and refolding (S) may be viewed as mappings, transforms, and/or operators suggesting the notation F−1, F, S−1, S. If the workflow is performed with the refaulting operator being the inverse of the unfaulting operator (
F 1*F=I) and the refolding operator being the inverse of the refolding operator (S−1*S=I), then the resulting framework should be identical to the existing framework. However, if the folding operator and/or the faulting operator are modified, then the resulting framework is a perturbation of the existing one. - For faulting or folding, the modifications may include, but are not limited to: shift, rotate, scale or deform fault; change the throws or fault type; split one fault, combine two faults, or replicate a fault to distribute the throws; shift or deform a surface, change layer thickness or lateral change rates; change the thickness between two surfaces or lateral thickness changes; or change the topology.
- Unfaulting, refaulting, unfolding, and refolding can be performed with different methods depending on the desired degree of accuracy. The methods range from purely geometric methods; to kinematic methods that attempt to preserve distances, areas, and volumes; to geomechanical methods that model stresses, strains, elasticity, plasticity, failure, etc.
- The operators may vary for different embodiments. For example, purely geometrical operators may suffice because similar assumptions are made for unfaulting (e.g., unfolding) and refaulting (e.g., refolding). Artifacts introduced by application of simplistic inverse operators are largely removed when applying the similarly simplistic forward operators. Because the objective of the workflow is generation of a perturbed framework, any remaining artifacts may be considered part of the perturbation.
- The geometrical operators may include unfaulting geometrical operators, faulting geometrical operators, folding geometrical operators and refolding geometrical operators. The unfaulting geometrical operators may be constructed from fault-horizon-intersection polygons by estimating an intermediate polyline from the hanging- and footwall polylines. Given these three polylines and the constraint that perturbations in the far field are minimal, an operator in form of an elevation perturbation map or delta-z map can be constructed. Application of a map-based depth modification is conventional functionality in many commercially available geologic modeling packages. Removing a first of multiple faults changes the intersections between any given horizon and the remaining faults. Thus, the remaining fault-horizon-intersection polygons may need to be perturbed with the first unfaulting operator before removal of a second fault and repeated for other similar operations. Faulting geometrical operators may be constructed by specification of a throw distribution for a given fault, constructing the horizon fault-intersection polygon for the specified fault and any given horizon, constructing an elevation perturbation map conditioned on this polygon and constrained to minimize far field perturbation, and applying this perturbation map or delta-z map to the specified horizon. Folding and refolding geometrical operators may be constructed, for example, by geostatistics. Fault geometry may be changed deterministically or geostatistically.
- Any other object that is specified by coordinates may also be transferred from one framework to another by application of some or all of these operators. Examples may include well picks, well paths, generic polylines, or geobodies that may be used for conditioning objects of the subsurface model. For example, perturbed horizons may be conditioned to well picks that themselves may be perturbed to account for their uncertainty. Preferably, the parameters (or geometry) of these operations is recorded to provide their application to some or all of the other objects in the model.
- Each of these different aspects may be combined into a system that provides systematic examination of uncertainty by automatically perturbing selected objects of the geologic framework by selective perturbation of faults, horizons, boundaries, their topologies and their geometries.
- In the present disclosure, an enhancement to exploration and reservoir delineation is described. In one or more embodiments, the method may include instantiating a realization of a subsurface model (e.g., reservoir model). The method may involve the creation of a conceptual subsurface model, the selection of parameter ranges for the various objects of the conceptual model and their interactions, and combining instantiated realizations of these objects into a subsurface model. The conceptual model may be generated by an agent (a user or a computer program that acts on behalf of a user) using a concept editor; may be automatically created from an inputted base realization; and/or may be created from the base case by systematically undoing faults and/or faults. The parameter ranges may be estimated from the undoing of faults and/or folds. Also, the unfaulting and refaulting may be performed based on fault-horizon cutoff polygons.
- Further, other embodiments may include other features. For example, the instantiated realizations may be analyzed for geologic plausibility and, if warranted, rejected; may be analyzed for technical consistency and, if warranted, rejected; may be further augmented with properties such as porosity, permeability, or oil saturation; and/or may be simulated. Also, the simulation may be performed using a simulation proxy method. In addition, a set of simulations related to different scenarios and/or realizations may be summarized with a statistics and/or may be analyzed to affect a decision. A connectivity measure may be used as a simulation proxy; and this connectivity measure may be based on a graph-based centrality measure. The centrality measure, which is described in U.S. patent application Ser. No. 14/272,581, which is incorporated by reference, may include one or more of degree, betweenness, closeness, and eigenvector. Further, the method may include ranking the plurality of objects, instantiated realizations, or other items in order of the respective centralization measures.
- In one or more embodiments, the method may include different concepts or a single concept. For example, the method may involve creating a concept once and multiple realizations that are created from the concept. For example, a concept may be a general layout of objects and relationships deemed to be of higher importance. The concept may not be drawn to scale, may not even show relative geometry or shape, and omit substantial amounts of detail. As another example, a realization is a model having more detail than a concept, which may be generated by a deterministic or stochastic process. The method may also include the creation of concepts that serve as scenarios. The method may also include multiple concepts that are created from multiple inputted base realizations to serve as different scenarios. Various aspects of the present techniques are described further in
FIGS. 1 to 10 . -
FIG. 1 is a flow chart for performing hydrocarbon exploration or delineation of a potential hydrocarbon resource in accordance with an exemplary embodiment of the present techniques. The method may include creating a concept and forming various realizations based on the concept. As may be appreciated, some blocks may be omitted, repeated, performed in different order, or augmented with additional blocks not shown. That is, some blocks may be performed sequentially, while others are executed simultaneously in parallel. - This
flowchart 100 begins with inputting a base realization, as shown inblock 101. The base realization includes is a subsurface model having various properties. Atblock 102, the concept is created. The creation of the concept may include reducing the base realization to a concept by systemic removal of all deformation. The removal may include unfaulting and unfolding. Inblock 103, parameter ranges are selected. The selection of parameter ranges may be performed by an interpreter or computer algorithm. The selection may include determining bounds for deformations and other model parameters. Interactions between two or more of the plurality of objects may also be defined at this point in the process. Then, atblock 104, the realizations are instantiated. The instantiating the realizations may include relying upon the parameter ranges (e.g., drawing parameters from these bounds) and selection of the interactions between two or more of the plurality of objects. Multiple realizations of the same concept or scenario may be obtained by systematic variation of parameters within their parameter ranges or by random sampling of the parameter ranges using a stochastic process. - As some combinations of parameters may result in realizations that are technically invalid or geologically implausible, a plausibility or validity check may be performed, as shown in
block 105. Inblock 105, the realization is checked for technical validity and/or geologic plausibility. For example, the realization may be verified by analyzing the validity of the resulting grid or by examination of fault polygons or thickness maps. Realizations that fail this check are either modified or discarded. The verification may include determining if the realization is geologically implausible and, if so, discarding the instantiated realizations that are geologically implausible. Similarly, the verification may include determining if the realization is technically consistency and, if so, discarding the instantiated realizations that are technically inconsistency. Examples of geologically implausible realizations include those with stratigraphically older horizons conformably disposed over younger horizons, faults interpreted as belonging to a prior episode of deformation offset faults known to have moved during a later episode, and/or faults having displacement-to-map length ratios outside of bounds defined by interpretation of real faults systems. Examples of technical inconsistency include, gaps or holes in what should be continuous horizon representations, duplicated portions of horizons or faults, horizons or faults that loop back upon themselves, mesh triangles that face the wrong directions and/or grid cells that are inside-out, mesh triangles with four or more vertices, or mesh triangles that intersect each other without an explicitly represented intersection edge. After the plausibility or validity check, the realization is then populated with properties, as shown inblock 106. The properties may include lithology, facies, porosity, permeability, fluid composition, and/or pressure. Once populated, the realization is simulated, as shown inblock 108. The simulation of the realization may include pressure or saturation changes as functions of spatial position and time, bypassed or disconnected resources, or the fluid composition produced at a specified well. - Then, at
block 109, a determination is made whether to reiterate the process for another realization. That is, the process may be repeated to instantiate other realizations. The determination may include creating a specific number of realizations, which may cover a specified part of the parameter space, cover a specified part of the response space, or design of experiments techniques to characterize the range of responses. If the determination is to perform another realization, the selection of the parameter ranges may be performed inblock 103. However, if another realization is not to be performed, a set of simulations of the realization(s) is analyzed, as shown inblock 110. The statistical or visual analysis may include ranking, whisper plots, box plots or other methods to review the different realizations. The analysis may include ranking the realizations based on a specified response, such as expected ultimate recovery or the maximum amount of oil, gas or water produced for a specified period of time. - Once the realizations are analyzed, the hydrocarbons are identified and produced, as shown in
blocks block 111, hydrocarbons may be identified based at least partially on the analysis of the realizations. As an example, the realizations may be integrated with other measured data or subsurface models of the subsurface regions below the survey region. These different types of data may be integrated based on location information associated with the respective data to lessen uncertainty associated with the existence of hydrocarbons. Finally, the identified hydrocarbons may be produced, as shown inblock 112. With the identification of hydrocarbons, various production operations may be performed to produce the hydrocarbons. For example, the operations may include drilling of a well to provide access to the hydrocarbon accumulation. Further, the production may include installing a production facility configured to monitor and produce hydrocarbons from the production intervals that provide access to the hydrocarbons in the subsurface formation. The production facility may include one or more units to process and manage the flow of production fluids, such as hydrocarbons and/or water, from the formation. The production equipment and operations may be based on the realizations. To access the production intervals, the production facility may be coupled to a tree and various control valves via a control umbilical, production tubing for passing fluids from the tree to the production facility, control tubing for hydraulic or electrical devices, and a control cable for communicating with other devices within the wellbore. - Beneficially, by using a concept, the present techniques provide a mechanism to provide a master subsurface model that may be used to generate other subsurface models. As noted above (e.g., in block 104), the present techniques involves generating perturbations from the concept to form various realizations. Then, the present techniques involve verification of the realizations (e.g., that the realizations are proper models). This provides a mechanism to lessen contamination by implausible models during subsequent statistical analysis.
- For the purpose of the present disclosure, the geologic model may be divided into a framework and content. The framework is formed by collections of volumes (three-dimensional), their bounding surfaces as well as other surfaces (two-dimensional), polylines or curves (one-dimensional), and points (zero dimensional) that are embedded in a three-dimensional space. Surfaces relate to an area of interest, faults, and horizons. Curves relate to surface intersections, such as fault-horizon intersections, fault sticks, or polygons and polylines, separating gross geologic features, such as environments of deposition. Horizons partition the model into zones; while faults partition the model into segments. Faults, horizons, and polygons partition the model into compartments.
- Content refers to the properties associated with compartments (e.g., three dimensional distribution of properties), surfaces (e.g., two-dimensional distribution of properties), and polylines (e.g., one-dimensional distribution of properties). A distinction between content and framework is the existence of a mesh or grid used to discretize properties.
- As noted above, concepts and realizations relate to different aspects, which are further explained in
FIG. 2 andFIG. 3 . These representations exemplify some of the differences between a realization, as shown inFIG. 2 , and a concept, as shown inFIG. 3 . -
FIG. 2 is a diagram of arealization 200 of a geologic model. Arealization 200 includes a frame of reference or coordinate system, as indicated by the coordinateaxes realization 200 has an area of interest or region ofinterest 203 that specifies the spatial extent of the model. Inside this region ofinterest 203, therealization 200 is completely quantified to enable numerical simulations. Outside this region ofinterest 203, therealization 200 is not specified and/or quantified. There may or may not be any data or information available about the region outside of the region ofinterest 203; but the region outside of the region ofinterest 203 is irrelevant because it is not modeled. In this diagram, therealization 200 includesfaults 204 and 205 (e.g.,faults realization 200, there is a crosscutting relationship between the two faults, wherefault 204 is dominant or major as compared tofault 205. Theminor fault 205 is realized by a set ofparallel faults - The
realization 200 also includeshorizons horizons fault 204. Thehorizons zones interest 203 or thehorizons realization 200, surfaces and zones, such aszones zone 209 has an attachedproperty 211 indicated by the gradual shading, such as porosity, net-to-gross ratio, or hydrocarbon saturation. The realized properties may be specified on a grid or mesh, such asmesh 212. The realization contains the necessary information to perform a specified computation or simulation, such as the estimation of the gross rock volume (GRV), the stock tank original oil in place (STOOIP), the expected ultimate recovery (EUR) or the prediction of water cuts. - Unlike the
realization 200, aconcept 300 does not need a frame of reference and/or a region of interest. For example,FIG. 2 is a diagram of aconcept 300 of therealization 200 ofFIG. 2 . Theconcept 300 includesfaults horizons zones concept 300 is typically unspecified, no frame of reference is utilized and horizons and faults are indicated by lines or planes, as shown byfaults horizons fault 304 appears to be to the left offault 305, but without geometry, the spatial arrangement of the faults cannot be specified. If desired, constraints on the spatial arrangement may subsequently be imposed with the selection of parameter ranges, as noted inblock 103 ofFIG. 1 . Also,horizon 306 appears tooverlay horizon 307, but again, without geometry, the spatial arrangement of the horizons cannot be specified. If desired, constraints on the spatial arrangement may subsequently be imposed with the selection of parameter ranges inblock 103 ofFIG. 1 . Preferably, however, the topology is augmented with the concepts of younger (shallower) and older (deeper) to capture the relative order of horizons in the conceptual model.Horizons horizons zones Zone 309 is bound, capped byhorizon 306, and based byhorizon 307.Zone 308 is unbound and based byhorizon 306.Zone 310 is unbound and capped byhorizon 307.Faults fault 304 is attributed with ‘major’ whilefault 305 is attributed ‘minor’, then a realization of 304 truncates a realization offault 305 in the event they intersect. - The comparison of these diagrams exemplifies some of the differences between a realization and a concept. For example, the
conceptual faults faults 204 and 205 (e.g.,faults realization 200, there is a crosscutting relationship between the two faults where 204 is dominant or major, while, in theconcept 300, bothfaults conceptual horizons horizons horizons fault 204, while the conceptual horizons do not have such properties. Further,zones interest 203 or thehorizons realization 200, surfaces and zones may have attached properties, while the conceptual zones do not include properties. - Again, the geologic model may be divided into a framework and content. Thus, the process of instantiating a realization is separated into two steps, which are i) instantiation of a framework realization and ii) instantiation of property realizations.
- As noted above, the present techniques involve the concept creation (e.g., blocks 101 and 102 of
FIG. 1 ) and instantiation of framework realizations (e.g., blocks 103 and 104). In some embodiments, the concept is created directly with a suitable tool. Preferably, however, the concept is created from a base realization (e.g., a realized geologic model that is stripped of geometry and possibly parts of its topology, meaning, and interpretation). Preferably, this base realization is the most likely reservoir model, a model synthesized from the optimal interpretation, or the model is the statistically centralized (e.g., in the middle of the groupings) with regard to a specified prediction. Concept creation by removal of geometry can be seen as the process of systematic undeformation (e.g., unfaulting and unfolding). Unfaulting and (re)faulting are discontinuous deformations, while unfolding and (re)folding are continuous deformations. Instantiating a realization by attaching geometry may be the process of systematic deformation (e.g., refaulting and refolding). - As noted above, unfaulting F−1, refaulting F, unfolding S−1, and refolding S can be viewed as mappings, transforms, or operators. If the method is performed with the refaulting operator being the inverse of the unfaulting operator (e.g., F−1*F=F*F−1=1) and the refolding operator being the inverse of the unfolding operator (e.g., S−1*S=S*S−1=1), then the resulting framework realization may be substantially identical to the existing base framework. However, if the refolding operator and/or the refaulting operator are modified, then the resulting framework realization may be a perturbation of the base framework.
-
FIG. 4 is a diagram 400 of a schematic application of the undeforming and redeforming a framework in accordance with an exemplary embodiment of the present techniques. In this diagram 400, the one (re)faulting operator and two (re)folding operators are utilized. The operators are perturbed, and their sequential order is commuted to create a realization caused by a different sequence of deformation events. Solid curves denote instantiated or realized objects, while dashed curves indicate conceptual objects.Block 402 depicts the base realization.Block 404 is the unfaulting operation F−1 that removes the discontinuities imposed by the fault from the horizons and renders the fault conceptual. The result of this operation is shown inblock 406.Block 408 is a first unfolding operation S1−1 that removes the effect of one regional deformation and renders one horizon conceptual. The result of this operation is shown inblock 410.Block 412 is a second unfolding operation S2−1 that removes the effect of another regional deformation and renders the second horizon conceptual.Block 414 depicts the resulting conceptual model. Instantiating one realization, as shown inblock 426, begins withblock 416, the application of the faulting operator F that instantiates a different type of fault, a reverse fault. The result of this operation is shown inblock 418. A first refolding operator, as shown inblock 420, deforms both the already realized fault and instantiates a first horizon. The result of this operation is shown inblock 422. The second refolding operator, as shown inblock 424, deforms both the realized fault and horizon and instantiates the second horizon. The result of this operation is shown inblock 426, which is the instantiated realization. This sequence of operators, as shown inblock 428, transformed thebase realization 402 to anotherrealization 426. Unfaulting, refaulting, unfolding, and refolding can be performed with different methods depending on the desired degree of accuracy. The methods range from purely geometric methods; to kinematic methods that attempt to preserve distances, areas, and volumes; and to geomechanical methods that model stresses, strains, elasticity, plasticity, failure, etc. - For faulting, the modifications include but are not limited to: shift, rotate, scale or deform a fault; change the throws; change the fault type; split one fault into a set of parallel or echelon faults; combine multiple faults into one; or change the topology or interaction between faults. For folding, the modifications include but are not limited to: shift or deform a horizon; change zone thickness or lateral change rates; change the horizon type; or change the topology or interaction between horizons. Polylines may be shifted, scaled, or deformed. Those skilled in the art should recognize that the abovementioned lists of modifications are only meant to be exemplary and not meant to be exhaustive.
- As a specific example, returning to
FIG. 1 , atblock 101, a base realization is inputted into the system that is systematically converted to a concept inblock 102 by stripping properties and geometries from the base realization using a sequence of unfaulting and unfolding operations. Preferably, at least some aspects of the stripped geometries and properties are retained to aid the selection of a parameter ranges inblock 103. Inblock 103, an agent (e.g., a user and/or a computer program that acts on behalf of a user) selects a set of at least one modification from the exemplary set of modifications and specifies parameters for these modifications. A fault may be parameterized by its base location, its base shape, a perturbation of its base shape, a scale to increase or reduce its extent, and a throw profile. Other parameters may include type and its position within the deformation sequence. A horizon may be parameterized by its base location, its base shape, a perturbation of its base shape, a type and its positions within the stratigraphic sequence and deformation sequence. A polygon may be parameterized by its base location, its base shape, a perturbation of its base shape, a scale, a shift, or the zone(s) to which it is applied to impart an environment of deposition specified by an agent. - Analogous to using undeformation operators to convert a base realization to a concept, deformation operators are used to instantiate a realization of the concept in
block 104. This may involve an agent specifying a sequence of continuous (e.g., refolding) and discontinuous (e.g., refaulting) deformations. The agent parameterizes the individual deformations and applies them to the concept objects of faults, horizons, and polylines. The realized faults and horizons may not intersect and truncate each other correctly. For example, horizons may need to be clipped or extended to the faults, while other horizons may need to clipped or extended to other horizons. Further, some faults may need to be clipped or extended to other faults. It may also be advantageous to determine the intersection curves between faults and/or horizons and assign these curves to the geometries of the intersecting objects (e.g., by using known processes, such as the creation of a watertight framework). Methods for clipping and extending objects and the subsequent creation of a watertight framework are known to those skilled in the art. For example, one such disclosure is U.S. Pat. No. 7,756,694. - The instantiated framework realization may then be converted into a three-dimensional grid bound by the area of interest. Details of the gridding process are known to those skilled in the art. Examples may include U.S. Patent Application Publication Nos. 2013/218539 and 2012/265510 along with U.S. Pat. No. 7,248,259.
- Based on the automatic instantiation of a framework realization in
block 104 and the succeeding grid realization, the framework and/or the grid may be geologically implausible and/or technically invalid.Block 105 checks the geologic plausibility and technical validity of the instantiated realization. If the framework is found invalid or implausible, this framework realization is removed from the workflow and/or at least flagged. Judicious parameterization and narrow parameter ranges may limit the instantiation of unacceptable realizations. A tradeoff, however, exists between plausibility/validity and variety of realization. Narrow ranges may yield little variety between realizations, but a greater number of plausible or valid ones. In some embodiments, parameters for the individual deformations are drawn independently for efficiency, but the resulting interaction of deformations may be far-fetched. An example of implausible realizations is shown below inFIGS. 5A, 5B and 5C . -
FIGS. 5A, 5B and 5C are diagrams ofimplausible realizations realizations Realization 500 changes polarity by slowly turning upside down.Realization 520 intersects itself andrealization 540 contains twofaults - Technical validity refers to the ability to create a mesh or grid associated with the instantiated framework realization. Some realizations may simply violate assumptions made by the gridding algorithm leading to an abnormal termination of the gridding process. In the worst case, the gridding algorithm may even crash. Other realizations may stretch the gridding algorithm beyond its design specifications, causing the generation of poor grids with a large number of cells with high aspect ratios, highly obtuse angles, or negative areas and volumes. Analysis of the grid generation process and the resulting grid realization provides a mechanism for removal or at least flagging of invalid realizations from the remainder of the process.
- Typically after attaching a grid to the framework realization, properties are instantiated, which is shown in
block 106. The realization is populated with properties such as porosity, net-to-gross ratio; oil, gas and water saturations, and horizontal and vertical permeabilities. The properties can be assigned deterministically, geostatistically, or by simulation; and conditioned or unconditioned with regard to other data, such as well logs or seismic data. Methods for instantiating properties in geologic models are known to those skilled in the art. An example is U.S. Pat. No. 7,415,401 to Calvert et al. entitled “Method for constructing 3-D geologic models by combining multiple frequency passbands”. - Preferably, properties (e.g., within the model) are conditioned, as shown in
block 107. The property or data may be conditioned or at least guided by well markers, well logs, maps, or seismic horizon attributes and seismic volume attributes. All of these conditioning data have geometry, but the present techniques create realizations of geologic models by perturbing, distorting, or modifying geometry. In some embodiments, it may be advantageous to modify the geometry of the conditioning data in the same manner by application of the same sequence of undeformations and redeformations used to create the concept and instantiate the realization. Some preferred embodiments may involve modifying the geometry of the conditioning data in an approximately similar manner only to allow for geometric uncertainty in these data caused by data acquisition, data processing, or interpretation. In one embodiment, some and/or all of the operators in the sequence may be modified or perturbed when applied to conditioning data. In a preferred embodiment, some additional operators are attached to the operator sequence for the conditioning data to model the geometric uncertainty of the conditioning data. Further, the conditioning of the properties may include using undeformation and redeformations on the data to be used in the conditioning. - The instantiated realization may then be simulated in
block 108 and/or analyzed inblock 110 to predict a specified quantity. Some predictions can be made directly from the subsurface model. Examples include but are not limited to gross rock volume (GRV), the stock tank original oil in place (STOOIP), or the expected ultimate recovery (EUR). Other predictions may involve additional financial assumption to calculate cash flow, discounted cash flow (DCF), discounted cash flow rate (DCFR), net present value (NPV), or return on capital employed (ROCE). Performing a reservoir simulation inblock 108 provides a prediction of water cuts, flow streams, flow capacity, storage capacity, connectivity, or some other performance indicator. - Instead of computing a complete fluid-flow simulation based on full-physics models that include state equations for oil, gas and water, multiphase Navier-Stokes equations, and a complete development/production scenario with producer wells, injector wells, injection rates, and perforation zones, it may be advantageous to use reduced-order or reduced-physics model, also termed a proxy model, to achieve computational efficiency and to reduce complexity by suppressing needless detail. Examples of such proxy simulations may be European Patent No. 1,994,488; U.S. Pat. No. 8,437,997 to Meurer et al entitled ‘Dynamic Connectivity Analysis’, U.S. Pat. No. 7,164,990 to Bratvedt et al entitled ‘Method Of Determining Fluid Flow’, or Hirsch and Schuette, ‘Graph Theory Applications To Continuity And Ranking In Geologic Models’, Computers & Geosciences, 25(2), 127-139, 1999.
- Once realized and possibly simulated, the analysis of the realization and/or the analysis of its simulation results may be performed, as noted in
block 110. Preferably, multiple realizations are instantiated and simulated and included as part of the analysis. - In certain embodiments, blocks 101 to 103 may be performed in an initial stage to convert the base realization to the concept with specified parameter ranges, while
blocks 104 to 108 may be repeated through multiple iterations (e.g., multiple times or stages) to generate multiple realizations and/or simulations for analysis. The number of repetitions may be controlled by the user or an agent directly by specification or indirectly by selection of a stopping criterion to ensure an appropriate sampling of the parameter space. - Often, predictions may exhibit transitions between different behaviors where perturbing parameters up to a certain point yields similar results, but perturbing the parameters beyond this point yields a very different result (e.g., different regimes). This behavior may be likened to phase transitions in thermodynamic systems where the system can abruptly move to a different state with very different properties. In one preferred embodiment, the stopping criterion attempts to predict the number of ‘states’ and locate the transitions between the discovered states in the parameter space. For example, the different regimes may involve predictions that impact the flow within the reservoir compartments and/or well. In particular, different regimes may include changes in flow that divide different compartments to adjust the amount of fluid communication between the compartments.
- In yet other embodiments, it may also be advantageous to input not only one base realization into the process of
blocks 101 to 108, but to iterate over multiple base realizations. Each base realization corresponds to a different scenario. A scenario is an alternative working hypothesis; or in the context of this disclosure, a scenario is an alternative concept. The workflow uses at least one concept that in a preferred embodiment is generated from a base realization. Preferably, however, multiple base realizations may be reduced to multiple concepts that differ from each other. Each of these different concepts may represent a different scenario. - The analysis or simulation of multiple realizations of one or multiple scenarios creates large amounts of data that may be visualized or summarized. In one embodiment, realizations are compared against each other by use of a metric that is used to group or cluster similar realizations (e.g., Suzuki et al, ‘Dynamic data integration for structural modeling: model screening approach using a distance-based model parameterization’, Computational Geosciences, 105-109, 2008). Techniques such as multi-dimensional scaling (MDS) may be used to group or cluster realizations and predictions.
- Reservoir simulation can create large amounts of time-dependent results or time-series data. In one embodiment, these time-series data are presented as contour boxplots (e.g., Sun and Genton, ‘Functional Boxplots’, Journal of Computational and Graphical Statistics, 20(2), 316-334, 2011).
- In another embodiment, the inputting of base realization(s) in
block 101 may be omitted. Instead of concept creation by conversion of inputted base realizations by systematic removal of deformation (unfaulting and unfolding), a user or agent may input or create at least one concept directly, for example, by using a concept editor inblock 102. A concept editor provides a mechanism for the creation of conceptual models by specifying at least the number of horizons, faults, and polylines. Preferably, the concept editor may be used to specify certain attributes, such as fault type, horizon type, environment of depositions, and their interactions. In a preferred embodiment, the concept editor creates objects for the specified entities directly in a geologic modeling software package where they can be operated on with a sequence of deformation operators F and S. In another embodiment, the concept editor creates objects for the specified entities either in memory, a file system, or in the cloud from where they can be imported by a geologic modeling software package to be operated on with a sequence of deformation operators F and S. - A preferred method of unfaulting is presented in
FIG. 6 .FIGS. 6A, 6B and 6C are diagrams 600, 620 and 640 of unfaulting during concept creation. The concept creation may be performed inblock 102 ofFIG. 1 .FIG. 6A is a diagram 600 having ahorizon 601 that is bisected bynormal fault 602 resulting in the foot-wall truncation 604 and hanging-wall truncation 603. Thetruncations horizon 601 againstfault 602. A reference, such asreference line 605, is created from the foot-wall and hanging-wall polylines. One method for creating the reference is simply to average the depths or two-way travel times of the cutoff polygon. Preferably, a local or floating reference is created for every polygon point. For a specified point of the foot-wall polyline, the laterally nearest point (e.g., neglecting vertical offset) of the hanging-wall polyline is determined and the local reference for the specified point is determined by averaging its depth with the depth of the nearest point on the hanging-wall polyline. For a specified point of the hanging-wall polyline, the laterally nearest point (e.g., neglecting vertical offset) of the foot-wall polyline is determined and the local reference for the specified point is determined by averaging its depth with the depth of the nearest point on the foot-wall polyline. The dynamic-time-warping (DTW) algorithm may be an efficient method to determine corresponding points on foot-wall and hanging-wall polylines. The residual polygon consisting ofresidual polylines 603′ and 604′ is created by subtraction of thereference 605 from the cutoff polygon formed bypolylines -
FIG. 6B is a diagram 620 having an unfaulting map orcorrection map 621 that is formed from the residual polylines 623 (corresponding to 603′ ofFIG. 6A ) and 624 (corresponding to 604′). Preferably, the map is formed by extrapolation from theresidual polylines -
FIG. 6C is a diagram 640 having anunfaulted horizon 641 that is formed by subtracting the map 643 (corresponding 621 ofFIG. 6B ) from the horizon 642 (corresponding to 601 ofFIG. 6A ). The reference or level of zero is shown byreference line 645. The effect offault 644 is removed fromhorizon 641.Fault 644 has been reduced to a concept as indicated by the dashes. The unfaulted horizon may contain gaps or artifacts that are preferably removed by filtering and interpolation. - When the base realization contains more than one fault, then removal of the first fault changes the horizon(s) and thus the cutoff polygon(s) for the remaining fault(s). Either the fault cutoffs should be recreated, or preferably, the original cutoff polygons may be corrected by subtraction of the correction map. For example, if the base realization contains three faults, then removal of the first fault triggers the correction of the cutoff polygons around the second and third faults. Subsequent removal of the second fault triggers another correction of the cutoff polygon around the third fault. Subsequent removal of the third fault does not trigger any further corrections because the all faults are reduced to concepts.
- Realizing a fault has two aspects: (i) the specification of all its parameters and (ii) the redeformation of the other objects in the model. For computational efficiency, it may be advantageous to create an explicit fault object between the first and the second aspects.
- A preferred method of refaulting is presented in
FIGS. 7A, 7B and 7C .FIGS. 7A, 7B and 7C are diagrams 700, 720 and 740 of refaulting during the instantiation of a realization. The instance realization may be performed inblock 104 ofFIG. 1 .FIG. 7A is a diagram 700 having a conceptualreverse fault 702 indicated by dashes that bisectshorizon 701. The process begins with parameterizing thefault 702 by specification of the fault geometry (e.g., location, orientation, shape, size, etc.). Some of these geometry parameters may be prescribed by the parameter ranges specified inblock 103 ofFIG. 1 , while others may be drawn at random from statistical distribution functions specified inblock 103 ofFIG. 1 . In another embodiment, combinations of parameters may be selected by systematic sampling of the parameter ranges. A fault throw is also specified, for example at theintersection 708 of the realized fault with the horizon or preferably for every location on the fault by designating fault throw as a property attached to the fault. The fault-horizon intersection 708 is also used to define thelocal reference depth 705. Fault type, reference, and throw are used to define the cutoff polygon consisting of foot-wall polyline 704 and hanging-wall polyline 703. For a reverse fault, the foot-wall polyline 704 is determined by shifting the fault-horizon intersection 708 downwards along the realizedfault 702 by half the local throw, while the hanging-wall polyline 703 is determined by shifting the fault-horizon intersection 708 upwards along the realizedfault 702 by half the local throw. For a normal fault, foot-wall and hanging-wall polylines are swapped: the foot-wall polyline 704 is determined by shifting the fault-horizon intersection 708 upwards along the realizedfault 702 by half the local throw, while the hanging-wall polyline 703 is determined by shifting the fault-horizon intersection 708 downwards along the realizedfault 702 by half the local throw. In a preferred embodiment, thepolylines residual polylines 703′ and 704′ are determined from thepolylines reference 705. -
FIG. 7B is a diagram 720 having a refaulting map or correction map 721 (consisting of 721′ and 721″) that is formed from the residual polylines 723 (corresponding to 704′ ofFIG. 7A ) and 724 (corresponding to 703′ ofFIG. 7A ). Preferably, the map is formed by extrapolation from theresidual polylines - Under the process of normal faulting, a flat, single-valued horizon remains single valued; and the
correction map 721 can be extrapolated from 723 and 724 directly without invoking 721′ and 721″. Under the process of reverse faulting, however, a flat, single-valued horizon will become multi valued. In the region between the foot-wall cutoff and the hanging-wall cutoff, the horizon will be duplicated and overlapping itself. Thus for a reverse fault, the refaulting map is multi valued between the foot-wall and the hanging-wall cutoffs shifting the meaning of map from a two-dimensional depiction of residual elevation toward a mathematical operator or transform. It may be advantageous to divide themulti-valued map 721 into the single-valuedmaps 721′ and 721″. -
FIG. 7C is a diagram 740 having arefaulted horizon 741 consisting of 741′ and 741″ that is formed by adding themap 742 consisting of 742′ (corresponding to 721′ ofFIG. 7B ) and 742″ (corresponding 721″ ofFIG. 7B ) to the horizon 744 (corresponding to 701 ofFIG. 7A ), while providing multivaluedness or overlap inhorizon 741 by use of an appropriate representation. The effect offault 743 is thus imparted ontohorizon 741.Fault 743 has been realized from a concept as indicated by the solid line. The reference or level of zero may be shown byreference line 745. For refaulting with a reverse fault, the unfaulted (conceptual) horizon piece inside the cutoff polygon is used twice as it gets added both to 742′ and to 742″. For refaulting a horizon with a normal fault, the unfaulted (conceptual) horizon piece inside the cutoff polygon would not be used at all. In either case, therefaulted horizon 741 may contain gaps or artifacts that are preferably removed by filtering and interpolation. Preferably, a process, such as disclosed in U.S. Pat. No. 7,756,694, is used to clean up the fault-horizon intersection by extrapolation of the horizon to the fault, cutback and truncation of the horizon by the fault, and creation of a watertight intersection. Preferably, the process may also be used to clean up fault-fault or horizon-horizon intersections and to create cutoff polygons. - When the concept contains more than one fault, then realization of the first fault changes the horizon(s) and thus the intersections with the remaining fault(s). Either the original cutoff polygons (intersections between fault and horizon, e.g., fault-horizon intersection 708) are corrected by addition of the correction map, or preferably, cutoffs at the remaining faults are recreated. For example, if the concept contains three faults, then realization of the first fault triggers the correction of the cutoff polygons around the second and third faults. Subsequent realization of the second fault triggers another correction of the cutoff polygon around the third fault. Subsequent realization of the third fault does not trigger any further corrections because the all faults have been realized. In some embodiments, the map (e.g., map 742) is added only to the horizon (e.g., horizon 744), while in others the map is also added to some or all of the faults that are already realized to preserve their relative positions during the refaulting operation.
-
FIG. 3 is a diagram 800 of the process from a base framework realization to an instantiated framework realization in accordance with an exemplary embodiment of the present techniques. Diagram 800 presents an exemplary application of the workflow from a base framework realization, such asbase framework realization 810, to an instantiated framework realization, such as such as instantiatedframework realization 860. - The process begins with the
base framework realization 810. Thebase framework realization 810 includes threefaults horizon 814. First, the three faults are converted to conceptual faults by removing their geometry and healing their effects on thehorizon 814. This may be performed by applying the unfaulting operators, such as first unfaulting operator F1−1 for thefirst fault 811, second unfaulting operator F2−1 for thesecond fault 812, and third unfaulting operator F3−1 for thethird fault 813. These different unfaulting operators may be combined: -
F 3 −1 *F 2 −1 *F 1 −1 (e1) - meaning that first the unfolding operator F1 −1 is applied, then F2 −1, and lastly F3 −1.
- The result is the
intermediary model 820 that contains threeconceptual faults horizon 824. Thehorizon 824 is still realized, but healed. Thehorizon 824 does not exhibit any spatial discontinuities.Horizon 824 is continuous, while thebase horizon 814 contained discontinuities at the fault locations. - Then, a reduction of the
continuous horizon 824 to theconceptual horizon 834 by removal of its geometry is performed. This may be performed by applying the unfolding operator S1−1. The result ismodel 830, which is the conceptual model having threeconceptual faults conceptual horizon 834. - Based on the specified parameter ranges, the
conceptual horizon 834 is reinstantiated creating the realizedhorizon 844 in thereinstantiated model 840. This may be performed by applying the folding operator S1. Without any realized faults, the realizedhorizon 844 is continuous, but clearly different fromhorizon 824. - Then, based on the specified parameter ranges, the
conceptual fault 822 is reinstantiated creating the realizedfault 852 inreinstantiated model 850. This may be performed by applying the folding operator F2. In this specific example, bothfaults fault 852, however, also introduced throws which that are applied to thecontinuous horizon 844 creating the faulted,discontinuous horizon 854. - By suppressing the non-instantiated
conceptual faults final framework realization 860 containsfault 852 andhorizon 854. Multiple framework realizations may be instantiated from theconcept 830 by using different parameterizations for theconceptual faults conceptual horizon 834. Creating multiple realizations may be useful to lessen uncertainty in the analysis of the realizations with regard to a specified problem, question, or decision. - As may be appreciated, the flow chart of
FIG. 1 may include various variations. For example, the concept may be created inblock 102. Inblock 103, an agent selects bounds for deformations and other model parameters. Then, parameters are selected from these bounds, and a realization of the concept is instantiated inblock 104. - In certain embodiments, the concept is generated in
block 102 by systematic removal of deformations from a base realization. Preferably, the realization is checked for technical validity and/or geologic plausibility inblock 105 because some combinations of parameters may result in realizations that are technically invalid or geologically implausible. Realizations that fail this test are either fixed or discarded outright. Then, the realizations are populated with properties inblock 106, simulated inblock 108, and analyzed inblock 110. The analysis results are summarized to facilitate business decisions and operations to produce hydrocarbons. - As an example,
FIG. 9 is a block diagram of acomputer system 900 that may be used to perform any of the methods disclosed herein. A central processing unit (CPU) 902 is coupled tosystem bus 904. TheCPU 902 may be any general-purpose CPU, although other types of architectures of CPU 902 (or other components of exemplary system 900) may be used as long as CPU 902 (and other components of system 900) supports the inventive operations as described herein. TheCPU 902 may execute the various logical instructions according to disclosed aspects and methodologies. For example, theCPU 902 may execute machine-level instructions for performing processing according to aspects and methodologies disclosed herein. - The
computer system 900 may also include computer components such as a random access memory (RAM) 906, which may be SRAM, DRAM, SDRAM, or the like. Thecomputer system 900 may also include read-only memory (ROM) 908, which may be PROM, EPROM, EEPROM, or the like.RAM 906 andROM 908 hold user and system data and programs, as is known in the art. Thecomputer system 900 may also include an input/output (I/O)adapter 910, acommunications adapter 922, auser interface adapter 924, and adisplay adapter 918. The I/O adapter 910, theuser interface adapter 924, and/orcommunications adapter 922 may, in certain aspects and techniques, enable a user to interact withcomputer system 900 to input information. - The I/
O adapter 910 preferably connects a storage device(s) 912, such as one or more of hard drive, compact disc (CD) drive, floppy disk drive, tape drive, etc. tocomputer system 900. The storage device(s) may be used whenRAM 906 is insufficient for the memory requirements associated with storing data for operations of embodiments of the present techniques. The data storage of thecomputer system 900 may be used for storing information and/or other data used or generated as disclosed herein. Thecommunications adapter 922 may couple thecomputer system 900 to a network (not shown), which may enable information to be input to and/or output fromsystem 900 via the network (for example, a wide-area network, a local-area network, a wireless network, any combination of the foregoing).User interface adapter 924 couples user input devices, such as akeyboard 928, apointing device 926, and the like, tocomputer system 900. Thedisplay adapter 918 is driven by theCPU 902 to control, through adisplay driver 916, the display on adisplay device 920. Information and/or representations of one or more 2D canvases and one or more 3D windows may be displayed, according to disclosed aspects and methodologies. - The architecture of
system 900 may be varied as desired. For example, any suitable processor-based device may be used, including without limitation personal computers, laptop computers, computer workstations, and multi-processor servers. Moreover, embodiments may be implemented on application specific integrated circuits (ASICs) or very large scale integrated (VLSI) circuits. In fact, persons of ordinary skill in the art may use any number of suitable structures capable of executing logical operations according to the embodiments. - In one or more embodiments, the method may be implemented in machine-readable logic, set of instructions or code that, when executed, performs a method to analyzing uncertainty of subsurface formations. The code may be used or executed with a computing system such as
computing system 900. The computer system may be utilized to store the set of instructions that are utilized to manage the data and other aspects of the present techniques. - As an example, a
computer system 900 may be used to analyze uncertainty of subsurface formations for production or exploration operations. The computer system may include a processor; memory in communication with the processor; and a set of instructions stored in memory and accessible by the processor. The set of instructions, when executed by the processor, are configured to: create a conceptual subsurface model, wherein the conceptual subsurface model is associated with a subsurface formation and comprises a plurality of objects; select parameter ranges for each of the plurality of objects and interactions between two or more of the plurality of objects; instantiate realizations for the plurality of objects based on the selected parameter ranges; and combine instantiated realizations of these objects into a reservoir model. The set of instructions are configured to create the conceptual subsurface model may be further configured to automatically create the conceptual subsurface model from an obtained base realization; may be further configured to undo one or more faults and folds in a sequential order; may be further configured estimate parameter ranges based on the undoing of one or more of faults and folds; and may be further configured to unfault an inputted base realization based on fault-horizon cutoff polygons to create the conceptual subsurface model. Further, the set of instructions may be configured to refault from the conceptual subsurface model based on fault-horizon cutoff polygons. - The computer system may include other instructions to enhance efficiency of the operation of the present techniques. For example, the set of instructions may be configured to analyze each of the instantiated realizations for geologic plausibility and, if one or more of the instantiated realizations are determined to be geologically implausible, discard the one or more of the instantiated realizations that are geologically implausible. In addition to or alternatively, the set of instructions may be configured to analyze each of the instantiated realizations for technical consistency and, if one or more of the instantiated realizations are determined to be technically inconsistency, discard the one or more of the instantiated realizations that are technically inconsistency.
- In other embodiment, the computer system may include other enhancements. For example, the set of instructions may be configured to instantiate properties into each of the instantiated realizations, wherein the properties comprise one or more of porosity, permeability and oil saturation; may be configured to condition the instantiating properties by perturbing, distorting, or modifying the geometry; and/or may be configured to condition the instantiating properties by applying a sequence of undeformations and redeformations that used to create the instantiated realizations. The set of instructions may be configured to simulate the instantiated realizations; may be configured to simulate using a simulation proxy method; may be configured to simulate proxy method using a connectivity measure as a simulation proxy for each of the instantiated realizations; may be configured to compute the connectivity measure based on graph based centrality measure; may be configured to rank the plurality of instantiated realizations in order of the respective centralization measures; and/or may be configured to simulate the instantiated realizations to create a set of simulations that are analyzed to affect a decision for production operations. Further, the set of instructions configured to create the conceptual subsurface model may be further configured to create two or more instantiated realizations from the concept conceptual subsurface model; may be further configured to create two or more conceptual subsurface models to generate two or more scenarios; and/or may be further configured to create one or more conceptual subsurface models that are each based on different base realizations and are created to generate two or more scenarios.
- In some preferred embodiments, simulation may be approximated by a simulation proxy. A preferred simulation proxy is based on graph-based centrality. The centrality measure, which is described in U.S. patent application Ser. No. 14/272,581, which is incorporated by reference, may include one or more of degree, betweenness, closeness, and eigenvector.
- A connectivity matrix expresses how well two neighboring grid cells are connected (transmissibility) or how similar a specified property is. In the first case, a connection is weighted; while in the second case, each grid cell is associated with a label or index i and an attribute or property value vi. The two cases are not mutually exclusive: one definition of connection weight is the magnitude of their property or attribute difference. Another preferred definition of connection weight is their property or attribute average. With this definition of connection weight, an off-diagonal element of the connectivity matrix Cij for two neighboring grid cells i and j (where i≠j) is set to −1/2(vi+vj). A diagonal element Cii of the connectivity matrix is set to Σ1/2εij(vi+vj) where εij is one when grid cells i and j are neighbors and zero when grid cells i and j not neighbors.
- In some preferred embodiments of the inventive method, the diagonal elements of the connectivity matrix are set to zero, effectively removing a self-interaction or self-connectivity.
- In some embodiments of the inventive method, specified eigenvectors of the connectivity matrix are used to compute a connectivity measure for the grid cell. The first component of the specified eigenvectors defines the location of the first grid cell in a vector space. The second component of the specified eigenvectors defines the location of the second grid cell in said vector space, and so on for the remaining components and grid cells. For a specified grid cell in said vector space, the shortest distance to any other grid cell in said space defines a measure of connectivity indicating how connected the specified grid cell is to all others. Iterating this process over substantially all grid cells provides a connectivity measure for substantially every grid cell, resulting in a connectivity attribute. For computational efficiency, it may be advantageous to limit for a specified grid cell the search of its nearest grid cell in the vector space. Instead of computing the distance to every other grid cell in said vector space, it is preferable to compute only the distance in said vector space to its original neighbors as indicated by the connectivity matrix.
- Details of the distance function are irrelevant. Different distance functions result in different connectivity measures. Any metric or any generalized metric associated with said vector space results in a connectivity measure.
- Instead of explicitly computing all or a few specified eigenvectors from the connectivity matrix and using these eigenvectors to compute a distance between grid cells, distances may be computed directly from the connectivity matrix using either an iterative or algebraic process. In the iterative process, the connectivity measure ci is computed iteratively as
-
- until a specified (convergence) criteria is satisfied where d is a small damping coefficient, Mij=1/Cij if Cij≠0 and zero otherwise, N refers to the number of grid cells, and 1 is a vector of dimension N containing only ones. An initial value for c may be 1/N. In the algebraic process,
-
- where I is an identity matrix. For computational efficiency, the iterative process is preferably used. U.S. Pat. No. 6,285,999 to Page discloses a method for ranking linked web pages based on similar mathematical notions.
- Depending on the specifics of the connectivity matrix, in some embodiments of the inventive method the connectivity matrix is normalized prior to the direct estimation of connectivity measures, for example by scaling each row sum, each columns sum, or each row sum and each column sum of the connectivity matrix C to one.
- In graph theory and network analysis, centrality of a vertex measures its relative importance within a graph. Examples include how influential a person is within a social network, how well-used a road is within an urban network, or how well connected the grid cells are within their geobodies or connectivity structures. There are four main measures of centrality: degree, betweenness, closeness, and eigenvector. The connectivity measures disclosed with this invention are examples of eigenvector-based centrality measures.
- Degree centrality refers to the number of connections for a specified node, potentially weighted by the attribute value. For the disclosed connectivity matrices, degree centralities or degree-based connectivity measure may be computed by row sums, column sums, or row-column sums, preferably excluding elements on the matrix diagonals from the sum.
- In a connected graph, there is a distance metric between any two grid cells belonging to this graph that is defined by the length of the shortest path between the two specified grid cells. The length of a path is defined by the number of connections linking the two specified grid cells, or in the attributed case, by the sum of the attributes along a path linking the specified grid cells. The farness of any grid cell is defined by the sum of its distances to all other grid cells of the graph. Closeness centrality is defined as the inverse of farness. The more central a grid cell is, the lower its total distances to all other grid cells. Closeness centrality can be viewed as a measure of how long it takes to spread information sequentially from a grid cell to all other grid cells belonging to the same graph.
- When using permeability to compute the connectivity matrix, the grid cells with the largest closeness centralities or the largest closeness-based connectivity measures are the grid cells that provide fast drainage of a contiguous group of grid cells from their fluids.
- Extensions of closeness centrality account not only for the shortest path length but also for the number of paths.
- Betweenness centrality quantifies the number of times a grid cell acts as a bridge along the shortest path between any two grid cells of a subsurface model. It may be advantageous to scrutinize grid cells with high betweenness centrality because a small perturbation to the connectivity structure, permeability or transmissibility might dramatically alter the shortest paths and their spatial distributions.
- Eigenvector centrality is a measure of the influence of a grid cell in the connected graph of the subsurface model. Eigenvector centrality assigns a relative score to all grid cells based on the principle that connections from a specified grid cell to high-scoring grid cells contribute more to the score of the specified grid cell than connections to low-scoring grid cells. The centrality score or eigenvector-based connectivity measure c can be defined as solution to the eigenvector equation C c=λc. There will typically be multiple eigenvalues λ for which an eigenvector solution exists. The dominant eigenvector associated with the largest eigenvalue is preferably obtained by an iterative process.
- In some embodiments of the inventive method, a centralization measure is computed for a group of contiguous grid cells that have been attributed with a specified connectivity measure. Centralization for the specified group of grid cells measures how central its most central grid cell is in relation to all of its other grid cells, for example by computation of Σcmax−ci. Preferably, this quantity is normalized by the number of grid cells or the theoretically largest sum of centrality differences for a graph of similar size. It may be advantageous to estimate the theoretically largest sum of centrality differences for a graph of similar size by constructing a compact group of grid cells with the same number of grid cells and maximal connectivity, for example in the shape of a ball. In the attributed case, every grid cell or connection of this ideal group is attributed with a maximal value in accordance to the specified attribute.
- In some embodiments of the inventive method, groups of contiguous grid cells (e.g., compartments, segments, zones) are ranked in order of their centralization measures. In some embodiment of the inventive method, the group of contiguous grid cells is formed by thresholding, by definition of a spatial bounding box, or by any other method.
- In some preferred embodiments of the inventive method, a connectivity measure is assigned to groups of contiguous grid cells of the reservoir model. The connectivity measure serves as a proxy to a reservoir simulation or reservoir performance analysis. Proxy simulations for performance prediction are well known to practitioners of the art. Examples of such proxy simulations may be European Patent No. 1,994,488 to Li et al entitled “Method for Quantifying Reservoir Connectivity Using Fluid Travel Times”, U.S. Pat. No. 8,437,997 to Meurer et al entitled ‘Dynamic Connectivity Analysis’, U.S. Pat. No. 7,164,990 to Bratvedt et al entitled “Method Of Determining Fluid Flow”, or Hirsch and Schuette, “Graph Theory Applications To Continuity And Ranking In Geologic Models”, Computers & Geosciences, 25(2), 127-139, 1999. All these proxies, however, are source-target proxies where some grid cells or cells are designated to be sources or injectors and other grid cells are designated as targets, sinks, or producers. Sources, targets and conductors (i.e., grid cells that are neither sources nor sinks) are mutually exclusive. The purpose of these proxies is the analysis of different reservoir development or production scenarios to examine the connectivity between the oil-bearing reservoir and the producer wells or the connectivity between water-injection wells and hydrocarbon-production wells. The novel connectivity measures disclosed in this publication are independent of sources and targets. No well locations need to be specified. Grid cells do not need to be separated into mutually exclusive sources, sinks, and conductors. Instead, each grid cell is compared to all others. Each grid cell acts simultaneously as source, sink, and conductor. The disclosed connectivity measures allow examination of the model for highly connected regions, for disconnected compartments, for barriers, and regions where small even perturbations of connectivity and attributes or properties (porosity, permeability, or transmissibility) will change long-distance connectivity by disconnecting one region or compartment into multiple ones or connecting multiple regions or compartments into one, thus warranting additional scrutiny to analyze these sensitive regions.
- It should be understood that the preceding is merely a detailed description of specific embodiments of the invention and that numerous changes, modifications, and alternatives to the disclosed embodiments can be made in accordance with the disclosure here without departing from the scope of the invention. The preceding description, therefore, is not meant to limit the scope of the invention. Rather, the scope of the invention is to be determined only by the appended claims and their equivalents. It is also contemplated that structures and features embodied in the present examples can be altered, rearranged, substituted, deleted, duplicated, combined, or added to each other. The articles “the”, “a” and “an” are not necessarily limited to mean only one, but rather are inclusive and open ended so as to include, optionally, multiple such elements.
Claims (42)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14/816,729 US10995592B2 (en) | 2014-09-30 | 2015-08-03 | Method and system for analyzing the uncertainty of subsurface model |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201462057797P | 2014-09-30 | 2014-09-30 | |
US14/816,729 US10995592B2 (en) | 2014-09-30 | 2015-08-03 | Method and system for analyzing the uncertainty of subsurface model |
Publications (2)
Publication Number | Publication Date |
---|---|
US20160090825A1 true US20160090825A1 (en) | 2016-03-31 |
US10995592B2 US10995592B2 (en) | 2021-05-04 |
Family
ID=55583866
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/816,729 Active 2037-08-13 US10995592B2 (en) | 2014-09-30 | 2015-08-03 | Method and system for analyzing the uncertainty of subsurface model |
Country Status (1)
Country | Link |
---|---|
US (1) | US10995592B2 (en) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130338986A1 (en) * | 2010-12-21 | 2013-12-19 | Statoil Petroleum As | Quality control of sub-surface and wellbore position data |
US9733388B2 (en) | 2008-05-05 | 2017-08-15 | Exxonmobil Upstream Research Company | Systems and methods for connectivity analysis using functional objects |
CN107422378A (en) * | 2017-07-14 | 2017-12-01 | 中国石油大学(华东) | A kind of comprehensive discrimination method of polygonal fault system |
WO2018022649A1 (en) * | 2016-07-29 | 2018-02-01 | Exxonmobil Upstream Research Company | Method and system for generating a subsurface model |
WO2018080745A1 (en) * | 2016-10-28 | 2018-05-03 | Exxonmobil Upstream Research Company | Method and system for enhancing meshes for a subsurface model |
CN109239301A (en) * | 2018-09-12 | 2019-01-18 | 江苏科技大学 | A kind of anchor chain flash welding quality online evaluation method |
US20200211127A1 (en) * | 2018-12-31 | 2020-07-02 | Exxonmobil Upstream Research Company | Methods and Systems for Performing Decision Scenario Analysis |
US10890688B2 (en) * | 2015-11-20 | 2021-01-12 | Korea Institute Of Geoscience And Mineral Resources | Method for generating secondary data in geostatistics using observed data |
US11041969B2 (en) * | 2018-06-22 | 2021-06-22 | Exxonmobil Upstream Research Company | Methods and systems for modeling subsurfaces containing partial faults |
US11221427B2 (en) | 2017-09-28 | 2022-01-11 | Chevron U.S.A. Inc. | Systems and methods for estimating reservoir productivity as a function of depth in a subsurface volume of interest |
US11269110B2 (en) * | 2018-12-05 | 2022-03-08 | Schlumberger Technology Corporation | Computing system assessment of geological similarity of wells employing well-log data |
US11287541B2 (en) * | 2015-12-18 | 2022-03-29 | Exxonmobil Upstream Research Company | Method to design geophysical surveys using full wavefield inversion point- spread function analysis |
US11313993B2 (en) * | 2017-09-28 | 2022-04-26 | Gopro, Inc. | Systems and methods for estimating a likelihood of reservoir productivity as a function of position in a subsurface volume of interest |
WO2022093875A1 (en) * | 2020-10-26 | 2022-05-05 | Jerry Nims | Subsurface imaging and display of 3d digital image and 3d image sequence |
US11391854B2 (en) * | 2019-07-04 | 2022-07-19 | Chengdu University Of Technology | Optimization design method for volumetric fracturing construction parameters of infilled well of unconventional oil and gas reservoir |
US11917119B2 (en) | 2020-01-09 | 2024-02-27 | Jerry Nims | 2D image capture system and display of 3D digital image |
Family Cites Families (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0021284D0 (en) | 2000-08-30 | 2000-10-18 | Schlumberger Evaluation & Prod | Compositional simulation using a new streamline method |
US7248259B2 (en) | 2001-12-12 | 2007-07-24 | Technoguide As | Three dimensional geological model construction |
AU2007225358B2 (en) | 2006-03-02 | 2012-11-15 | Exxonmobil Upstream Research Company | Method for quantifying reservoir connectivity using fluid travel times |
US7756694B2 (en) | 2006-05-31 | 2010-07-13 | Schlumberger Technology Corporation | Method for interactive automation of fault modeling including a method for intelligently sensing fault-fault relationships |
US7585132B2 (en) * | 2006-06-27 | 2009-09-08 | James Imbrie | Method for remediating a contaminated site |
WO2009094064A1 (en) * | 2008-01-22 | 2009-07-30 | Exxonmobil Upstream Research Company | Dynamic connectivity analysis |
WO2009137150A1 (en) | 2008-05-09 | 2009-11-12 | Exxonmobil Upstream Research Company | Method for geophysical and stratigraphic interpretation using waveform anomalies |
WO2009142872A1 (en) | 2008-05-22 | 2009-11-26 | Exxonmobil Upstream Research Company | Seismic horizon skeletonization |
US8213261B2 (en) | 2008-05-22 | 2012-07-03 | Exxonmobil Upstream Research Company | Method for geophysical and geological interpretation of seismic volumes in the domains of depth, time, and age |
US9228415B2 (en) * | 2008-10-06 | 2016-01-05 | Schlumberger Technology Corporation | Multidimensional data repository for modeling oilfield operations |
AU2009314461B2 (en) * | 2008-11-14 | 2015-11-05 | Exxonmobil Upstream Research Company | Forming a model of a subsurface region |
CA2754690C (en) * | 2009-03-11 | 2017-01-24 | Exxonmobil Upstream Research Company | Gradient-based workflows for conditioning of process-based geologic models |
CA2764681C (en) | 2009-07-06 | 2017-11-21 | Exxonmobil Upstream Research Company | Method for seismic interpretation using seismic texture attributes |
CA2771865C (en) | 2009-10-23 | 2016-04-05 | Exxonmobil Upstream Research Company | Method for optimization with gradient information |
WO2011056347A1 (en) * | 2009-11-05 | 2011-05-12 | Exxonmobil Upstream Research Company | Method for creating a hierarchically layered earth model |
US9703006B2 (en) * | 2010-02-12 | 2017-07-11 | Exxonmobil Upstream Research Company | Method and system for creating history matched simulation models |
CA2799504C (en) * | 2010-05-28 | 2018-01-02 | Exxonmobil Upstream Research Company | Method for seismic hydrocarbon system analysis |
US20130151161A1 (en) * | 2010-08-27 | 2013-06-13 | Matthias G. Imhof | Seismic Horizon Skeletonization |
US20120099396A1 (en) * | 2010-10-22 | 2012-04-26 | Chevron U.S.A. Inc. | System and method for characterization with non-unique solutions of anisotropic velocities |
WO2012058029A1 (en) * | 2010-10-27 | 2012-05-03 | Exxonmobil Uspstream Research Comapny | Method and system for fracture stimulation |
CA2816931A1 (en) * | 2010-12-08 | 2012-06-14 | Exxonmobil Upstream Research Company | Constructing geologic models from geologic concepts |
US8965745B2 (en) | 2011-04-14 | 2015-02-24 | Schlumberger Technology Corporation | Grid from depositional space |
US9128204B2 (en) * | 2011-04-15 | 2015-09-08 | Exxonmobil Upstream Research Company | Shape-based metrics in reservoir characterization |
US8749549B2 (en) * | 2011-06-10 | 2014-06-10 | Schlumberger Technology Corporation | Prospect assessment and play chance mapping tools |
US8902221B2 (en) * | 2011-06-10 | 2014-12-02 | Schlumberger Technology Corporation | Prospect assessment and play chance mapping tools |
WO2013012469A1 (en) | 2011-07-20 | 2013-01-24 | Exxonmobil Upstream Research Company | Seismic pattern recognition using attributed relational graphs |
US9176247B2 (en) * | 2011-10-06 | 2015-11-03 | Exxonmobil Upstream Research Company | Tensor-based method for representation, analysis, and reconstruction of seismic data |
WO2013081708A1 (en) | 2011-11-29 | 2013-06-06 | Exxonmobil Upstream Research Company | Method for quantitative definition of direct hydrocarbon indicators |
AU2013219911B2 (en) | 2012-02-13 | 2014-11-06 | Exxonmobil Upstream Research Company | System and method for detection and classification of seismic terminations |
FR2987138A1 (en) | 2012-02-22 | 2013-08-23 | Schlumberger Technology Corp | TITLE NOT RENTED. |
US10422900B2 (en) | 2012-11-02 | 2019-09-24 | Exxonmobil Upstream Research Company | Analyzing seismic data |
WO2014074213A2 (en) | 2012-11-07 | 2014-05-15 | Exxonmobil Upstream Research Company | Method for knowledge capture and pattern recognition for the detection of hydrocarbon accumulations |
WO2014149344A1 (en) | 2013-03-15 | 2014-09-25 | Exxonmobil Upstream Research Company | Method and system for geophysical modeling of subsurface volumes |
US9952340B2 (en) | 2013-03-15 | 2018-04-24 | General Electric Company | Context based geo-seismic object identification |
US9995844B2 (en) | 2013-03-15 | 2018-06-12 | Exxonmobil Upstream Research Company | Method and system for geophysical modeling of subsurface volumes |
WO2014193619A2 (en) | 2013-05-30 | 2014-12-04 | Exxonmobil Upstream Research Company | Automated interetaton error correction |
US9824135B2 (en) | 2013-06-06 | 2017-11-21 | Exxonmobil Upstream Research Company | Method for decomposing complex objects into simpler components |
WO2016001752A2 (en) * | 2014-06-30 | 2016-01-07 | Cgg Services Sa | Ensemble-based multi-scale history-matching device and method for reservoir characterization |
-
2015
- 2015-08-03 US US14/816,729 patent/US10995592B2/en active Active
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9733388B2 (en) | 2008-05-05 | 2017-08-15 | Exxonmobil Upstream Research Company | Systems and methods for connectivity analysis using functional objects |
US20130338986A1 (en) * | 2010-12-21 | 2013-12-19 | Statoil Petroleum As | Quality control of sub-surface and wellbore position data |
US10890688B2 (en) * | 2015-11-20 | 2021-01-12 | Korea Institute Of Geoscience And Mineral Resources | Method for generating secondary data in geostatistics using observed data |
US11287541B2 (en) * | 2015-12-18 | 2022-03-29 | Exxonmobil Upstream Research Company | Method to design geophysical surveys using full wavefield inversion point- spread function analysis |
WO2018022649A1 (en) * | 2016-07-29 | 2018-02-01 | Exxonmobil Upstream Research Company | Method and system for generating a subsurface model |
WO2018080745A1 (en) * | 2016-10-28 | 2018-05-03 | Exxonmobil Upstream Research Company | Method and system for enhancing meshes for a subsurface model |
US10901118B2 (en) * | 2016-10-28 | 2021-01-26 | Exxonmobil Upstream Research Company | Method and system for enhancing meshes for a subsurface model |
CN107422378A (en) * | 2017-07-14 | 2017-12-01 | 中国石油大学(华东) | A kind of comprehensive discrimination method of polygonal fault system |
US11221427B2 (en) | 2017-09-28 | 2022-01-11 | Chevron U.S.A. Inc. | Systems and methods for estimating reservoir productivity as a function of depth in a subsurface volume of interest |
US11313993B2 (en) * | 2017-09-28 | 2022-04-26 | Gopro, Inc. | Systems and methods for estimating a likelihood of reservoir productivity as a function of position in a subsurface volume of interest |
US11041969B2 (en) * | 2018-06-22 | 2021-06-22 | Exxonmobil Upstream Research Company | Methods and systems for modeling subsurfaces containing partial faults |
CN109239301A (en) * | 2018-09-12 | 2019-01-18 | 江苏科技大学 | A kind of anchor chain flash welding quality online evaluation method |
US11269110B2 (en) * | 2018-12-05 | 2022-03-08 | Schlumberger Technology Corporation | Computing system assessment of geological similarity of wells employing well-log data |
US20200211127A1 (en) * | 2018-12-31 | 2020-07-02 | Exxonmobil Upstream Research Company | Methods and Systems for Performing Decision Scenario Analysis |
US11391854B2 (en) * | 2019-07-04 | 2022-07-19 | Chengdu University Of Technology | Optimization design method for volumetric fracturing construction parameters of infilled well of unconventional oil and gas reservoir |
US11917119B2 (en) | 2020-01-09 | 2024-02-27 | Jerry Nims | 2D image capture system and display of 3D digital image |
WO2022093875A1 (en) * | 2020-10-26 | 2022-05-05 | Jerry Nims | Subsurface imaging and display of 3d digital image and 3d image sequence |
Also Published As
Publication number | Publication date |
---|---|
US10995592B2 (en) | 2021-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10995592B2 (en) | Method and system for analyzing the uncertainty of subsurface model | |
US10838108B2 (en) | Method and system for interpolating discontinuous functions in a subsurface model | |
USRE49507E1 (en) | Faulted geological structures having unconformities | |
US10572611B2 (en) | Method and system for characterizing fractures in a subsurface region | |
EP3329307B1 (en) | Assignment of systems tracts | |
US10901118B2 (en) | Method and system for enhancing meshes for a subsurface model | |
EP3293552B1 (en) | System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles | |
Cherpeau et al. | Stochastic structural modelling in sparse data situations | |
US20180031719A1 (en) | Method and System for Generating a Subsurface Model | |
US11662501B2 (en) | Geologic modeling methods and systems having constrained restoration of depositional space | |
EP3526628B1 (en) | Geologic structural model generation | |
US11506807B2 (en) | Methods and systems for simulation gridding with partial faults | |
US11409024B2 (en) | Methods and systems for generating simulation grids via zone by zone mapping from design space | |
US11041976B2 (en) | Method and system for creating and using a subsurface model in hydrocarbon operations | |
US20220163692A1 (en) | Modeling and simulating faults in subterranean formations | |
CN112292714B (en) | Grid partition based on fault radiation | |
Kurianova et al. | Comprehensive Analysis of the Geological and Geophysical Data in the Study of the Upper Miocene Turbidite Systems of the Nam Con Son Basin, Vietnam |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: FINAL REJECTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: FINAL REJECTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT RECEIVED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT VERIFIED |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |