EP3683601B1 - Abbildung eines unterirdischen geologischen modells zu einer vergangenen wiederherstellungszwischenzeit - Google Patents

Abbildung eines unterirdischen geologischen modells zu einer vergangenen wiederherstellungszwischenzeit Download PDF

Info

Publication number
EP3683601B1
EP3683601B1 EP20150807.4A EP20150807A EP3683601B1 EP 3683601 B1 EP3683601 B1 EP 3683601B1 EP 20150807 A EP20150807 A EP 20150807A EP 3683601 B1 EP3683601 B1 EP 3683601B1
Authority
EP
European Patent Office
Prior art keywords
fault
time
model
restoration
past
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.)
Active
Application number
EP20150807.4A
Other languages
English (en)
French (fr)
Other versions
EP3683601A1 (de
Inventor
Jean-Laurent Mallet
Anne-Laure Tertois
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Aspen Paradigm Holding LLC
MALLET, JEAN-LAURENT
Original Assignee
Aspen Paradigm Holding LLC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Aspen Paradigm Holding LLC filed Critical Aspen Paradigm Holding LLC
Publication of EP3683601A1 publication Critical patent/EP3683601A1/de
Application granted granted Critical
Publication of EP3683601B1 publication Critical patent/EP3683601B1/de
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V20/00—Geomodelling in general
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures
    • G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05—Geographic models
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00—Manipulating 3D models or images for computer graphics
    • G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/51—Migration
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/64—Geostructures, e.g. in 3D data cubes
    • G01V2210/642—Faults
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/66—Subsurface modeling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2219/00—Indexing scheme for manipulating 3D models or images for computer graphics
    • G06T2219/20—Indexing scheme for editing of 3D models
    • G06T2219/2021—Shape modification

Definitions

  • Embodiments of the invention relate to the field of geological tomography for generating an image of the interior subsurface of the Earth based on geological data collected by transmitting a series of incident waves and receiving reflections of those waves across discontinuities in the subsurface.
  • the incident and reflected waves are reconstituted by a 3D model to generate an image of the reflecting surfaces interior to the Earth. Accordingly, geological tomography allows geophysicists to "see inside" the Earth.
  • Embodiments of the invention further relate to geological restoration in which the tomographic images of the present day geology are transformed into images of the past geology, as it was configured at an intermediate restoration time in the past ⁇ before the present day and after the start of deposition of the oldest subsurface layer being imaged.
  • New techniques are proposed herein to improve both the accuracy and computational speed of generating those images of the past restored geology.
  • Improved images may aid geoscientists exploring the subsurface geology for applications such as predicting tectonic motion or earthquakes, or by engineers in the mining or oil and gas industries.
  • the accuracy of a geological model of the present day configuration of the subsurface of the Earth may be improved by "restoring" the model to a past intermediate time ⁇ and checking model consistency at that time in the past.
  • restoring geological models is a complex task and current methods are typically inefficient, requiring extensive processing resources and time, as well as inaccurate, relying on over-simplifications that induce errors to moderate the complexity of the task.
  • US2017330373A1 discloses a geometric method for 3D structural restoration of a subsurface model including receiving data representative of a subsurface volume of interest including one or more chronohorizons and the geometry and topology of any faults of relevance; developing a fault framework model of the subsurface volume of interest; selecting a horizon, the deposition of which represents the geologic time to which the structural model should be restored; developing coordinate transformation constrained by a single datum horizon and, optionally, additional geologic constraints; applying the 3D transformation to all geologic features below and, optionally, above the datum surface; and scaling the vertical coordinates to accurately relate vertical and horizontal dimensions.
  • the method may be executed by a computer system.
  • US2013204598A1 discloses a method and system for modeling a subsurface structure at a time when the structure was originally formed.
  • a first model having non-planar horizons representing a current subsurface structure may be used to calculate a vector field based on the non-planar geometry of the horizons of the model.
  • the vector field may be non-uniform or uniform.
  • Geographic coordinates of the first model may be transformed to paleo-geographic coordinates of a model representing the subsurface structure in the past, where the non-planar horizons in the first model are transformed to planar horizons in the second model.
  • a set of points describing one or more fractures in the subsurface structure may be used to calculate a tuning parameter to correct a first set of paleo-geographic coordinates.
  • a second set of coordinates representing an improved prediction at a time period when the subsurface structure was originally formed may be generated.
  • the invention is directed to modeling restored geological models with ⁇ -active and ⁇ - inactive faults.
  • a system according to claim 12 and a computer-implemented method according to claim 1 are provided.
  • the computer-implemented method entails restoring a 3D model of the subsurface geology of the Earth from a present day geometry measured, using geological or seismic tomography of subsurface structures, at a present time to a predicted past geometry at a past restoration time.
  • the 3D model of the present day measured geometry comprising a network of faults is received, wherein a fault is a discontinuity that divides fault blocks that slide in opposite directions tangential to the surface of the fault as time approaches a modeled time.
  • a past restoration time ⁇ is selected that is prior to the present time and after a time when an oldest horizon surface in the 3D model was originally deposited.
  • the network of faults is divided into a subset of ⁇ -active faults and a subset of ⁇ -inactive faults, wherein a ⁇ -active fault is a fault that is active at the past restoration time ⁇ and a ⁇ -inactive fault is a fault that is inactive at the past restoration time ⁇ .
  • a fault is determined to be ⁇ -active when the fault intersects a horizon H ⁇ that was originally deposited at the past restoration time ⁇ and a fault is determined to be ⁇ -inactive when the fault does not intersect the horizon H ⁇ that was originally deposited at the past restoration time ⁇ .
  • the 3D model is restored from the present day measured geometry to the predicted past geometry at the past restoration time ⁇ by modeling each ⁇ -active and ⁇ -inactive fault differently.
  • Each ⁇ -active fault is modeled to join end points of a horizon H ⁇ separated on opposite sides of the fault in the present day model to merge into the same position in the restored model by sliding the end points towards each other in a direction tangential to the surface of the ⁇ -active fault.
  • Each ⁇ -inactive fault is modeled to keep collocated points on opposite sides of the fault together. Images of the interior subsurface of the Earth are generated for the past restored geology based on the restored 3D model.
  • Some embodiments that are not encompassed by the wording of the claims but are considered as useful for understanding the invention, comprise modeling restored geological models with new restoration coordinates u ⁇ , v ⁇ , t ⁇ .
  • a system and method is provided for restoring a 3D model of the subsurface geology of the Earth from a present day measured geometry to a predicted past geometry at a restoration time in the past ⁇ .
  • the 3D model of the present day geometry of the subsurface may be received, including one or more folded geological horizon surfaces.
  • a value may be selected of a restoration time in the past ⁇ before the present day and after a time an oldest horizon surface in the 3D model of the subsurface was deposited.
  • the 3D model may be restored from the present day measured geometry to the predicted past geometry at the restoration time in the past ⁇ using a 3D transformation.
  • the vertical component of the 3D transformation may restore the geometry to the vertical coordinate t ⁇ such that: points along a horizon surface H ⁇ modeling sediment that was deposited at the selected restoration time in the past ⁇ have a substantially constant value for the restored vertical coordinate t ⁇ ; and at any location in the 3D model, the restored vertical coordinate t ⁇ is equal to a sum of a first approximation t' ⁇ of the vertical coordinate and an error correction term ⁇ ⁇ , wherein the error correction term ⁇ ⁇ is computed by solving a linear relationship in which a variation in the sum of the first approximation t' ⁇ of the vertical coordinate and the error correction term ⁇ ⁇ between any two points separated by an infinitesimal difference in the direction of maximal variation of the sum is approximately equal to the distance between the points in the direction of maximal variation; and displaying an image of the restored 3D model of
  • Embodiments of the invention solve this problem by selectively activating and deactivating individual fault surfaces to be discontinuous or continuous, respectively, depending on the specific restoration geological-time ⁇ .
  • embodiments of the invention split faults into two complementary subsets of " ⁇ - active" faults and " ⁇ -inactive" faults, ⁇ -active faults are activated at restoration time ⁇ (e.g., a discontinuous fault surface along which fault blocks slide tangentially), whereas ⁇ -inactive faults are deactivated at restoration time ⁇ (e.g., a continuous surface that does not behave as a fault).
  • FIG. 8 shows the problem of a fault 300 being erroneously considered active at a restoration time before it formed (top image of FIG. 8 ) and the solution of modeling the fault as a ⁇ -inactive fault to correctly deactivate the fault at restoration time ⁇ according to embodiments of the invention (bottom image of FIG. 8 ).
  • top image of FIG. 8 when a horizon H ⁇ 210 is restored using u ⁇ v ⁇ t ⁇ -transform 201, fault block 800 (shaded region in the top-left image of FIG. 8 ) is bounded by an active fault 105 and an inactive fault 300. If however fault 300 is erroneously considered as an active fault then, after applying restoration 201:
  • fault 105 is modeled as a ⁇ -active fault (activating the fault)
  • fault 300 is modeled as a ⁇ -inactive fault (deactivating the fault). Accordingly, when u ⁇ v ⁇ t ⁇ -transform 201 is applied to fault block 800 (shaded region in the bottom-left image of FIG. 8 ), restored fault block 801 (shaded region in the bottom-right image of FIG. 8 ) is no longer bounded by an active fault ( ⁇ -inactive fault 300 is inactive at restored time ⁇ ).
  • the restored fault block 801 preserves volume and stays within (and does not cross) ⁇ -active fault 805 (because the deactivated boundary transformed from ⁇ -inactive fault 300 may shift to accommodate a shift in the restored ⁇ -active fault 805 ).
  • ⁇ -active and ⁇ -inactive faults produces more accurate results, e.g., even if there is no continuous path between (no way to continuously connect) a given fault block (e.g., 800 ) and the horizon H ⁇ (e.g., 210 ) deposited at geological time ⁇ , which typically requires additional processing that may induce errors.
  • a given fault block e.g., 800
  • the horizon H ⁇ e.g., 210
  • embodiments of the invention eliminate erroneous phantom faults and more accurately represent the faulted geology.
  • FIG. 16 is a flowchart of a method to restore a geological model using ⁇ -active and ⁇ -inactive faults, according to an embodiment of the invention.
  • a processor receives a 3D model of the present day measured geometry comprising a network of faults (e.g., present day model 202 ).
  • the present day model is measured tomographically by scanning the Earth's subsurface e.g., as described in reference to FIGS. 14 and 15 .
  • all faults in the model have already formed and so, represent active discontinuities that divide fault blocks which slide in opposite directions tangential to the surface of the fault as time approaches a modeled time.
  • a processor selects or receives a past restoration time ⁇ that is "intermediate” or prior to the present time and after the start of the subsurface's deposition (the time period when an oldest horizon surface in the 3D model was originally deposited).
  • a processor divides the network of faults into a subset of ⁇ -active faults and a subset of ⁇ -inactive faults
  • ⁇ -active faults may be faults that are active at the past restoration time ⁇
  • ⁇ -inactive faults are faults that are inactive at the past restoration time ⁇ .
  • a fault is determined to be ⁇ -active when the fault intersects a horizon H ⁇ that was originally deposited at the past restoration time ⁇ (e.g., see ⁇ -active faults 105 of FIG.
  • a fault is determined to be ⁇ -inactive when the fault does not intersect the horizon H ⁇ that was originally deposited at the past restoration time ⁇ (e.g., see ⁇ -inactive faults 300 of FIG. 3 ).
  • the processor may divide the fault network differently at different geological times.
  • a fault may be ⁇ - active at a first restoration time ⁇ ' (e.g., a time period during which the fault has formed) and ⁇ - inactive at a second restoration time ⁇ " (e.g., a time period different than that during which the fault has formed).
  • iso-value surfaces e.g., 208 of FIG.
  • each restoration coordinate e.g., u ⁇ , v ⁇ , and t ⁇
  • ⁇ -inactive faults e.g., 300 of FIG. 3
  • ⁇ -active faults e.g., 105 of FIG. 3
  • a processor restores the 3D model from the present day measured geometry to the predicted past geometry at the past restoration time ⁇ .
  • the processor may flatten a horizon H ⁇ (e.g., 210 of FIG. 4 ) that was originally deposited at time ⁇ to a substantially planar surface of approximately constant depth.
  • the processor may restore the horizons to non-planar surfaces, e.g., when the thickness of the layers is not constant.
  • restoring the 3D model to a past restoration time ⁇ may eliminate (e.g., removing or not displaying) all relatively shallower horizon surfaces that were originally deposited after the past restoration time ⁇ .
  • the processor treats ⁇ -active and ⁇ -inactive faults differently in operations 1650 and 1660, respectively.
  • a processor models the ⁇ -active fault as an active discontinuous fault surface and restore the horizon surface by removing or omitting the fault surface at the time of restoration.
  • the processor eliminates the ⁇ -active fault during restoration by sliding its adjacent fault blocks together. This joins end points of a horizon H ⁇ separated on opposite sides of the fault in the present day model to merge into the same position in the restored model by sliding the end points towards each other in a direction tangential to the surface of the ⁇ -active fault.
  • a processor models the ⁇ -inactive fault, not as a discontinuous fault surface, but as a continuous non-fault surface in the restoration transformation.
  • the ⁇ -inactive fault is modeled as a surface in which the discontinuity induced by the fault has been deactivated to prevent fault blocks from sliding in directions tangential to the surface of the fault as time approaches the restoration time ⁇ .
  • the processor models the ⁇ -inactive fault during restoration by keeping collocated points on opposite sides of the fault in the present day model together in the restored model.
  • the process may repeat to restore the model for a second different past restoration time ⁇ '.
  • the geological model may be sequentially restored to a sequence of multiple past restoration times ⁇ 1 , ⁇ 2 , ..., ⁇ n .
  • the fault network may be divided into a different subset of ⁇ - active and ⁇ -inactive faults, e.g., because different faults fracture the subsurface at different geological times.
  • a processor may play a moving image sequence in which the 3D model is iteratively restored in a forward or reverse order of the sequence of past restoration times ⁇ 1 , ⁇ 2 , ..., ⁇ n to visualize changes in the subsurface geology over the passage of time.
  • a processor displays a visualization of an image of the subsurface geology of the Earth overlaid with ⁇ -active faults and ⁇ -inactive faults in the restored model at past restoration time ⁇ .
  • the processor may display the ⁇ -active faults and the ⁇ -inactive faults with different visual identifiers, such as, different levels of translucency, different colors, different patterns, etc.
  • a restoration transformation may transform a geological image of the subsurface of the Earth from a present day space (e.g., x,y,z coordinates) to a restoration space (e.g., u ⁇ , v ⁇ , and t ⁇ coordinates) as it was formed at an intermediate restoration time in the past ⁇ (before the present-day but after the start of the subsurface deposition).
  • An ideal restoration should transform the vertical coordinate t ⁇ in a manner that strictly honors the thickness of layers, to preserve areas and volumes of the Earth, so that terrains are not stretched or squeezed over time in the vertical dimension.
  • conventional restoration transformations typically deform the vertical coordinates, forcing terrains to stretch and squeeze, resulting in errors in the restoration model.
  • Embodiments of the invention improve the accuracy of the restoration model by establishing a vertical restoration coordinate t ⁇ that preserves layer thickness.
  • This may be achieved by implementing a thickness-preserving constraint that sets a variation in the vertical restoration coordinate t ⁇ between any two points separated by an infinitesimal difference in the direction of maximal variation of the vertical coordinate t ⁇ to be approximately equal to the distance between the points in the direction of maximal variation.
  • Embodiments of the invention improve the accuracy of the restored model by establishing a new thickness-preserving constraint that introduces an error correction term ⁇ ⁇ .
  • the new thickness-preserving constraint sets the restored vertical coordinate t ⁇ to be equal to a sum of a first approximation t' ⁇ of the vertical coordinate and an error correction term ⁇ ⁇ , wherein the error correction term ⁇ ⁇ is computed by solving a relationship in which a variation in the sum of the first approximation t' ⁇ of the vertical coordinate and the error correction term ⁇ ⁇ between any two points separated by an infinitesimal difference in the direction of maximal variation of the sum is approximately equal to the distance between the points in the direction of maximal variation.
  • the new thickness-preserving constraint preserves layer thickness with greater accuracy as shown in histogram 503 of FIG. 5 as compared to conventional approximations shown in histograms 501 and 502 of FIG. 5 and minimizes volume variation with greater accuracy as shown in histogram 903 of FIG. 9 as compared to conventional approximations shown in histograms 901 and 902 of FIG. 9 , respectively.
  • Embodiments of the invention further improve the performance and computational speed of the computer generating the restored model by linearizing the new thickness-preserving constraint.
  • 1.
  • the error correction term ⁇ ⁇ may be generated such that the square of its spatial variation, ⁇ grad ⁇ ⁇ ⁇ 2 , is negligible.
  • the thickness-preserving constraint simplifies to a new linear thickness-preserving constraint of grad ⁇ ⁇ ⁇ grad t' ⁇ ⁇ 1 ⁇ 2 ⁇ 1- ⁇ grad t' ⁇ ⁇ ⁇ 2 (eqn. (37)).
  • This thickness-preserving constraint is linear because t' ⁇ is already known, so the constraint is a relationship between the gradient of the error ⁇ ⁇ and the gradient of the known first approximation of the vertical coordinate t ' ⁇ .
  • the computational complexity for performing the restoration transformation according to embodiments of the invention is significantly reduced compared to classical methods that are based on the mechanics of continuous media.
  • the modeling computer uses significantly less computational time and storage space to generate the inventive restoration model.
  • embodiments of the invention implement a new set of geometrical constraints and boundary conditions that preserve geological volumes and deformations while adhering to geological boundaries.
  • embodiments of the invention restore faults along fault striae (e.g., see FIG. 12 ) induced by the twin points associated with the paleo-geographic coordinates of a depositional (e.g., GeoChron) model, given as input of the restoration method.
  • a depositional e.g., GeoChron
  • An ideal restoration should also transform the horizontal coordinates u ⁇ and v ⁇ in a manner that strictly honors lateral spatial distribution, to preserve areas and volumes of the Earth, so that terrains are not stretched or squeezed over time in the horizontal dimensions.
  • conventional restoration transformations based on depositional coordinates e.g., paleo-geographic coordinates u and v
  • deform the horizontal coordinates forcing terrains to stretch and squeeze, resulting in errors in the restoration model.
  • Embodiments of the invention improve the accuracy of the restoration model at time ⁇ by establishing horizontal restoration coordinates u ⁇ and v ⁇ that restore the horizon surface H ⁇ deposited at time ⁇ consistently with horizontal depositional coordinates u and v whilst minimizing deformations.
  • the horizontal restoration coordinates u ⁇ and v ⁇ are equal to the depositional coordinates u and v (see e.g., equation (20)) and the spatial variations of the horizontal restoration coordinates u ⁇ and v ⁇ are preserved with respect to the horizontal depositional coordinates u and v (see e.g., equation (21)).
  • each restoration model at time ⁇ presents a horizon surface H ⁇ , as it was configured at that time ⁇ when it was originally deposited.
  • horizontal restoration coordinates u ⁇ and v ⁇ are modeled in a tectonic style (e.g., using constraints (22) or (23)) that is consistent with that of the horizontal coordinates u and v of the depositional model, which makes the restoration more accurate because the geological context is taken into account.
  • horizontal restoration coordinates u ⁇ and v ⁇ are modeled to minimize deformations induced by the restoration of horizon H ⁇ , rather than minimizing deformations in the whole volume G.
  • Horizontal restoration coordinates u ⁇ and v ⁇ may also be constrained only in G ⁇ , thereby only taking into account the part of the subsurface to be restored, not the entire model G.
  • horizontal restoration coordinates u ⁇ and v ⁇ may be constrained to be equal on opposite sides of ⁇ -active faults at twin point locations, where the twin points are computed from fault striae, which also ensures consistency with the depositional model (see e.g., equation (43)). Additionally or alternatively, horizontal restoration coordinates u ⁇ and v ⁇ are constrained to be equal on opposite sides of ⁇ -inactive faults at mate point locations to cancel the effect of inactive faults on the restoration model (see e.g., equation (43)).
  • FIG. 17 is a flowchart of a method to restore a geological model with improved accuracy using a new thickness-preserving constraint, according to an embodiment of the invention.
  • a processor may receive a 3D model of the present day measured geometry (e.g., present day model 202 ) comprising one or more folded (e.g., curvilinear or non-planar) geological horizon surfaces (e.g., 210 ).
  • the present day model may be measured tomographically by scanning the Earth's subsurface e.g., as described in reference to FIGS. 14 and 15 .
  • a processor may select or receive a past restoration time ⁇ that is "intermediate” or prior to the present time and after the start of the subsurface's deposition (the time period when an oldest horizon surface in the 3D model was originally deposited).
  • a processor may restore the 3D model from the present day measured geometry (e.g., present day model G ⁇ 202 in xyz-space G 220 ) to the predicted past geometry at the restoration time in the past ⁇ (e.g., restored model G ⁇ 203 in u ⁇ v ⁇ t ⁇ -space 219 ) using a 3D restoration transformation (e.g., u ⁇ v ⁇ t ⁇ -transform 201 ).
  • the geological layers above H ⁇ e.g., H ⁇ + 1 ...
  • the 3D restoration transformation includes a vertical component that restores the geometry to the vertical coordinate t ⁇ and two lateral or horizontal components that restore the geometry to the horizontal coordinates u ⁇ and v ⁇ .
  • the restored vertical coordinate t ⁇ and horizontal coordinates u ⁇ and v ⁇ represent the predicted vertical and horizontal positions, respectively, where particles in the subsurface were located in the Earth at the restoration time in the past ⁇ .
  • the processor may restore and compute coordinates for the part or subregion G ⁇ of the subsurface G that was deposited at a geological time of deposition t prior to or during the past restoration time ⁇ (e.g., deposited deeper than, or at the same layer in the subsurface as, H ⁇ ). Accordingly, the restored model eliminates or omits all relatively shallower or younger horizon surfaces or layers that were originally deposited after the past restoration time ⁇ .
  • the processor may restore the vertical coordinate t ⁇ such that points along a horizon surface H ⁇ (e.g., 210 ) modeling sediment that was deposited at the selected restoration time ⁇ have a substantially constant value for the restored vertical coordinate t ⁇ (see e.g., eqn. (19)).
  • the processor may restore the vertical coordinate t ⁇ such that at any location in the 3D model, the restored vertical coordinate t ⁇ is equal to a sum of a first approximation t' ⁇ of the vertical coordinate and an error correction term ⁇ ⁇ , wherein the error correction term ⁇ ⁇ is computed by solving a relationship in which a variation in the sum of the first approximation t' ⁇ of the vertical coordinate and the error correction term ⁇ ⁇ between any two points separated by an infinitesimal difference in the direction of maximal variation of the sum is approximately equal to the distance between the points in the direction of maximal variation.
  • the error correction term ⁇ ⁇ may correct errors in the first approximation t' ⁇ of the vertical coordinate. This constraint may be represented by a linear second order approximation (see e.g., eqn. (37)).
  • the processor computes the first approximation t' ⁇ of the vertical coordinate by solving a relationship in which the spatial variation of the vertical coordinate t' ⁇ is locally approximately proportional to the spatial variation of a geological time of deposition t.
  • the coefficient of proportionality is locally equal to the inverse of the magnitude of the maximal spatial variation of the geological time of deposition ( see e.g., eqn. (34)-(1)).
  • This relationship may give the vertical restoration coordinate t ⁇ the shape of the horizon H ⁇ because, on the horizon, the gradient of depositional time t is normal to the horizon surface.
  • follows the shape of the horizon.
  • the processor computes the first approximation t' ⁇ of the vertical coordinate by solving a relationship in which any infinitesimal displacement in the direction orthogonal to horizon surface H ⁇ results in a variation of the vertical coordinate t' ⁇ approximately equal to the length of the infinitesimal displacement for points on the horizon surface H ⁇ ( see e.g., eqn. (33)-(1)).
  • the processor computes the restored vertical coordinate t ⁇ in parts of the subsurface which are older than restoration time ⁇ such that iso-value surfaces of the restored vertical coordinate t ⁇ are parallel to the horizon surface H ⁇ and the difference in the restored vertical coordinate t ⁇ between two arbitrary iso-values is equal to the distance between the corresponding iso-surfaces (see e.g., eqn. (31)).
  • Parallel surfaces may be planar parallel in the restored model, and curved parallel (e.g., having parallel tangent surfaces) in present day model, such that the surfaces are non-intersecting at limits.
  • the error correction term ⁇ ⁇ is null at points along the horizon surface H ⁇ that was deposited at the selected restoration time in the past ⁇ so that the restored horizon surface H ⁇ is flat ( see e.g., eqn. (36)).
  • the restored horizontal coordinates u ⁇ and v ⁇ are constrained such that for each point along the horizon surface H ⁇ that was deposited at the selected restoration time in the past ⁇ : the restored horizontal coordinates u ⁇ and v ⁇ are equal to depositional horizontal coordinates u and v , respectively, and the spatial variations of the restored horizontal coordinates u ⁇ and v ⁇ are equal to the spatial variations of the depositional horizontal coordinates u and v, respectively (see e.g., eqns. (20)-(21)).
  • grad v ⁇ 1.
  • the processor may not constrain grad u ⁇ to be orthogonal to grad t ⁇ . This results from the boundary condition on H ⁇ and propagation through its constant gradient.
  • the restored horizontal coordinates u ⁇ and v ⁇ are constrained in parts of the subsurface which are older than restoration time ⁇ such that directions of maximal change of the restored horizontal coordinates u ⁇ and v ⁇ are linearly constrained by a local co-axis vector b ⁇ and a local axis vector a ⁇ , respectively ( see e.g., eqn. (41)).
  • the local axis vector a ⁇ is oriented approximately in the direction of maximal change of depositional horizontal coordinate u and orthogonal to the direction of maximal change of the vertical restoration coordinate t ⁇
  • the local co-axis vector b ⁇ is oriented orthogonal to the direction of the local axis vector a ⁇ and orthogonal to the direction of maximal change of the vertical restoration coordinate t ⁇ ( see e.g., eqn. (40)).
  • the restored horizontal coordinates u ⁇ and v ⁇ are computed over the part of the 3D model of the subsurface which is older than restoration time ⁇ such that the directions of maximal change of u ⁇ and v ⁇ are approximately orthogonal to the local co-axis vector b ⁇ and the local axis vector a ⁇ , respectively.
  • equation (40) constrains the local axis vector a ⁇ to be parallel to the gradient of u and the local co-axis vector b ⁇ to be orthogonal to the local axis vector a ⁇ , which means that the gradient of u is orthogonal to the local co-axis vector b ⁇ .
  • Equation (41) further constrains the gradient of u ⁇ to be approximately orthogonal to the local co-axis vector b ⁇ . Accordingly, the gradient of u ⁇ is approximately parallel to the gradient of u. The same logic implies the gradient of v ⁇ is approximately parallel to the gradient of v.
  • the restored horizontal coordinates u ⁇ and v ⁇ are computed over the part of the 3D model of the subsurface which is older than restoration time ⁇ such that projections of their directions of maximal change over the iso-value surfaces of the restored vertical coordinate t ⁇ are approximately orthogonal to local co-axis vector b ⁇ and the local axis vector a ⁇ , respectively (see e.g., eqn. (42)).
  • the values of the restored horizontal coordinates u ⁇ and v ⁇ are constrained in parts of the subsurface which are older than the restoration time ⁇ to be respectively equal on twin points on ⁇ -active faults, wherein twin points are points on opposite sides of a ⁇ -active fault that were collocated at the restoration time ⁇ and are located on the same fault stria in the present day model, to merge the twin points into the same position in the restored model by sliding the twin points towards each other in a direction tangential to the surface of the ⁇ -active fault (see e.g., eqn. (43)).
  • the values of the restored horizontal coordinates u ⁇ and v ⁇ are constrained in parts of the subsurface which are older than the restoration time ⁇ to be respectively equal on mate points on ⁇ -inactive faults, wherein mate points are points on opposite sides of a ⁇ -inactive fault that are collocated at present day time, to move mate points together on opposite sides of ⁇ -inactive faults (see e.g., eqn. (43)).
  • a processor may display an image of the restored 3D model of the subsurface geology of the Earth such that each point in the 3D model is positioned at the restored coordinates u ⁇ , v ⁇ , t ⁇ defining the location that a piece of sediment represented by the point was located at the restoration time in the past ⁇ .
  • the processor may receive an increasing chronological sequence of past restoration times ⁇ 1 , ⁇ 2 , ..., ⁇ n .
  • the processor may repeat operations 1720-1730 to compute a corresponding 3D restoration transformation R ⁇ i .
  • 3D restoration transformation R ⁇ i restores the part of the subsurface older than horizon H ⁇ i to its predicted past geometry at time ⁇ i , e.g., to 3D restored coordinates u ⁇ i , v ⁇ i , and t ⁇ i .
  • a processor may play a moving image sequence in which the 3D model is iteratively restored in a forward or reverse order of the sequence of past restoration times ⁇ 1 , ⁇ 2 , ..., ⁇ n to visualize changes in the subsurface geology over the passage of time.
  • the processor may edit the model in the restoration space and then reverse the restoration transformation to apply those edits in the present day space.
  • the processor may edit the depositional values u, v, and t associated with the restored 3D model, and then reverse transform the restored 3D model forward in time from the predicted past geometry at the restoration time in the past ⁇ to the present day measured geometry using an inverse of the 3D restoration transformation 200 to incorporate the edits from the restored model into the present day model.
  • FIGS. 16 and 17 may be performed for example using system 1505 of FIG. 15 , e.g., by one or more processor(s) 140 of FIG. 15 , or another suitable computing system.
  • the embodiments disclosed in reference to FIGS. 16 and 17 may be performed using other operations or orders of the operations, and the exact set of steps shown in the figures may be varied.
  • Embodiments of the invention provide a new, purely geometrical 3D restoration method based on the input of a depositional (e.g., GeoChron model).
  • a depositional e.g., GeoChron model
  • Embodiments of the invention are able to handle depositional models of any degree of geometrical and topological complexity, with both small and large deformations, do not assume elastic mechanical behavior, and do not require any prior knowledge of geo-mechanical properties.
  • Embodiments of the invention further reduce or eliminate gaps and overlaps along faults as part of the restoration transformation and do not resort to any post-processing to minimize such gaps and overlaps.
  • embodiments of the invention minimize deformations and volume variations induced by geological restoration with a higher degree of precision, unequaled so far ( see e.g., FIG. 5 and FIG. 9 ).
  • Embodiments of the invention further ensure that 2D deformations of horizon surfaces induced by the uvt -transform are kept coherent with 3D deformations of volumes
  • the set of faults is split into ⁇ -active and ⁇ -inactive subsets.
  • Embodiments of the invention input a 3D model of sedimentary terrains in the subsurface.
  • the input model may be the GeoChron TM model generated by SKUA ® software for use in mining and oil and gas industries.
  • Embodiments of the invention may build a 3D restoration transformation of this model in such a way that, after transformation, the new model represents terrains as they were at a given intermediate restoration-time ⁇ (where ⁇ 1 ⁇ ⁇ ⁇ ⁇ 2 , before the present day ⁇ 2 and after the time of the deposition of the oldest layer ⁇ 1 ).
  • G may represent the present day 3D geological domain of the region of the subsurface being modeled and G ⁇ 202 may represent the subset of G containing particles of sediment that were deposited at a time prior to or equal to ⁇ .
  • a depositional model may be generated by inputting a tomographic model of the present day subsurface geology of the Earth and transforming that geology to a past depositional time as each particle was configured when originally deposited in the Earth. Sedimentary particles are deposited over time in layers from deepest to shallowest from the earliest to the most recent geological time periods. Since various layers of terrain are deposited at different geological times, a depositional model does not image the geology at any one particular time period, but across many times periods, each layer modeled at the geological time when the layer was deposited. Accordingly, the vertical axis or depth in the depositional model may be a time dimension representing the time period of deposition, progressing from oldest to newest geological time as the model progresses vertically from deepest to shallowest layers.
  • the depositional model may be the GeoChron TM model, which is generated by SKUA TM software, that is routinely used by many oil & gas companies to build models of geologic reservoirs which help optimize hydrocarbon production and exploration.
  • GeoChron TM model is generated by SKUA TM software, that is routinely used by many oil & gas companies to build models of geologic reservoirs which help optimize hydrocarbon production and exploration.
  • An example implementation of the GeoChron model is discussed in U.S. Patent No. 8,600,708 : The depositional model is described in reference to the GeoChron model only for example, though any other depositional model may be used.
  • FIG. 7 schematically illustrates an example transformation from a present day model (upper-left image) to a depositional GeoChron model (bottom-right image), according to an embodiment of the invention.
  • the new depositional location r has a vertical coordinate that is the geological time t ( r ) when the particle at location r was deposited and has horizontal or paleo-geographic coordinates ⁇ u ( r ), v ( r ) ⁇ equal to the lateral spatial location where the particle at r was located at its depositional time t ( r ).
  • the paleo-geographic coordinates ⁇ u ( r ), v ( r ) ⁇ may be linked to the vertical time coordinate t ( r ) by different relationships (e.g., constrained by different systems of differential equations) depending on the structural style of their deposition (e.g., minimal deformation or flexural slip).
  • the uvt-transform 700 shown in FIG. 7 when the geological time coordinate t ( r ) is equal to the curvilinear distance to the top horizon H ⁇ 210 along curvilinear axis 222, the uvt -transform is a valid technique for imaging the depositional model.
  • the depositional (e.g., GeoChron) model includes the following data structures stored in a memory (e.g., memory 150 of FIG. 15 ) (see FIGS. 1 , 3 , 6 , and 7 ):
  • the depositional model may have the following properties:
  • GeoChron model and its features described herein are discussed only as an example of a depositional model and that these elements may differ in other models or implementations without changing the essence of the invention.
  • the restoration time ⁇ may be a given geological time in the past and subdomain G ⁇ 202 may be a part of a 3D present day geological domain G that has terrains older than (deposited at a time prior to) or equal to restoration time ⁇ and defined by a depositional model.
  • Embodiments of the invention provide a new unfolding technique that replaces the uvt -transform (converting the present day model to a depositional model rendering all layers at their many respective times of deposition) by a u ⁇ v ⁇ t ⁇ -transform 201 (converting the present day model to a restored model at a single restoration time ⁇ before the present day but after the earliest times of deposition of the deepest model layer): x y z ⁇ u ⁇ v ⁇ t ⁇ u ⁇ x y z , v ⁇ x y z , t ⁇ x y z Accordingly, present day geological space G ⁇ 202 is transformed into a restored geological space G ⁇ 203, such that:
  • Compaction may be handled in pre and post-restoration stages, as is known in the art. Thus, the model may be restored without taking compaction into account.
  • Some embodiments of the invention provide an inventive volume deformation with a new set of inventive geometric constraints on the depositional model to allow geologic layers to be restored at a given geological time ⁇ with a precision that has never before been reached.
  • paleo-geographic coordinates ⁇ u ( r ) and v ( r ) ⁇ and the geological time coordinate t ( r ) are replaced by new restoration coordinates respectively denoted ⁇ u ⁇ ( r ), v ⁇ ( r ) ⁇ and t ⁇ ( r ).
  • the coordinate frame unit vectors ⁇ r u ⁇ , r v ⁇ , r t ⁇ ⁇ 219 of the G ⁇ -space and its origin O u ⁇ v ⁇ t ⁇ may be equal to the coordinate frame unit vectors ⁇ r x , r y , r z ⁇ 220 of the G -space and its origin O xyz : r ⁇ u ⁇ ⁇ r x ; r ⁇ ⁇ ⁇ ⁇ r y ; r ⁇ t ⁇ ⁇ r z O ⁇ u ⁇ ⁇ ⁇ t ⁇ ⁇ O xyz
  • z ⁇ O may be assumed to be constant, e.g., equal to zero.
  • FIG. 2 shows the folded present-day volume G ⁇ 202 resulting from the deformation of restored volume G ⁇ 203 under tectonic forces following either a "minimal deformation” or a "flexural slip” tectonic style: G ⁇ ⁇ ⁇ tectonic forces ⁇ G ⁇ ⁇ G
  • the part of the subsurface observed today stratigraphically below H ⁇ 210 may be identified with the deformed volume G ⁇ 202, e.g., such that:
  • restoration coordinates ⁇ u ⁇ ( r ), v ⁇ ( r ), t ⁇ ( r ) ⁇ may be chosen in such a way that the u ⁇ v ⁇ t ⁇ -transform 201 of the present-day volume G ⁇ 202 into the restored volume G ⁇ 203 minimizes deformations and volume variations. This is achieved by constraining restoration coordinates ⁇ u ⁇ ( r ), v ⁇ ( r ), t ⁇ ( r ) ⁇ to honor the two following conditions in the present day G ⁇ domain:
  • the task of restoration includes:
  • S ⁇ (0) ( 209 ) is assumed to be a continuous, unfaulted surface whose altitude z ⁇ O is a given function z ⁇ O ( u , v ).
  • z ⁇ O may refer to a given function z ⁇ O ( u ( r ) , v ( r )) which may or may not be constant: ⁇ r ⁇ H ⁇ : z ⁇ o stands for z ⁇ o u r , v r
  • a given depositional (e.g., GeoChron) model may be received from storage in a digital device (e.g., from memory 150 of FIG. 15 ).
  • a geological time ⁇ may be selected that is associated with the given horizon H ⁇ 210 to be restored and the given altitude z ⁇ O of the surface S ⁇ (0) 209 onto which the horizon H ⁇ should be restored.
  • the region G ⁇ 202 may be retrieved as the part of the depositional model where geological time of deposition t ( r ) is less than or equal to ⁇ (subsurface regions deposited in a layer deeper than or equal to the layer deposited at time ⁇ ).
  • the set of faults may be split into a subset of ⁇ -active faults cutting H ⁇ 210 and a subset of ⁇ -inactive faults which do not cut H ⁇ .
  • a geologist or other user may decide to manually transfer some faults from the ⁇ - inactive fault set to the ⁇ -active set, or vice versa, which typically causes greater restoration deformations. For example, manually altering the set of automatically computed ⁇ -active and ⁇ -inactive faults typically makes the restoration less accurate.
  • r F ⁇ r F ⁇ ⁇ ( 304,306 ) represents a pair of "mate-points" collocated on both sides of F 300 and assigned to F + 103 and F - 104, respectively, and ⁇ ⁇ ( r ) represents an error correction constraint.
  • Constraints (25), (26), (27) and (28) may be referred to collectively as "fault transparency constraints.”
  • fault transparency constraints (25), (26), (27) and (28) may be locally installed at any point r F on a ⁇ -active fault F where fault throw is lower than TH min .
  • TH min may be equal to 1 meter.
  • Two new discrete vector fields r ⁇ and R ⁇ may be defined on 3D grid ⁇ 100.
  • the depositional uvt -transform 700 of G ⁇ 202 is typically correct when equation (3) is honored.
  • embodiments of the present invention adapt equation (3) for the inventive restoration technique, replacing the vertical depositional coordinate t ( r ) by a vertical restoration coordinate t ⁇ ( r ) and replacing equation (3) by the following inventive thickness-preserving constraint to ensure layer thickness is preserved and surfaces ⁇ S ⁇ ( d ) : d ⁇ 0 ⁇ are parallel: ⁇ grad t ⁇ r ⁇ ⁇ 1 ⁇ r ⁇ G ⁇
  • thickness-preserving equation (31) cannot be implemented as a DSI constraint, which must be linear.
  • various linear surrogates of equation (31) may be used to approximate t ⁇ ( r ) as follows:
  • Constraints (33) and (34) are only examples of possible surrogate-thickness-preserving constraints. Other examples of such surrogate thickness-preserving constraints may be used.
  • new inventive constraint (34)-(1) benefits from the geologic observation that, throughout the entire domain G ⁇ 202, surfaces ⁇ S ⁇ ( d ) : d ⁇ 0 ⁇ 208 generally have a shape roughly similar to the level sets of the geologic time of deposition t ( r );
  • a first approximation of vertical restoration coordinate t ' ⁇ ( r ) may be computed by running the DSI method on grid ⁇ 100.
  • FIG. 5 shows histograms 501 and 502 of the distributions of ⁇ grad t ⁇ ⁇ , where t ⁇ is approximated using constraints (33) or (34), respectively, in the example geological terrain G ⁇ 202 shown in FIG. 4 .
  • FIG. 5 shows that when t ⁇ is approximated by constraints (33) or (34), ⁇ grad t ⁇ ⁇ significantly differs from "1". Therefore, while constraints (34) provide a better approximation of the thickness-preserving equation (31) than constraints (33), both of these approximations are inaccurate.
  • FIG. 9 shows histograms 901 and 902 of relative variations of volume ⁇ V/V induced by the restoration of H ⁇ 210 over G ⁇ 202 shown in FIG. 4 , where t ⁇ is approximated using constraints (33) or (34), respectively.
  • a restoration transformation should minimize variations in volume ⁇ V /V from the present day to the restored model.
  • FIG. 9 shows that a restoration based on constraints (33) or (34) results in a volume variation ⁇ V /V that significantly differs from the ideal value of "0". While constraints (34) result in a smaller volume variation ⁇ V/V than constraints (33), both of these approximations induce a significant volume variation ⁇ V/V and are inaccurate.
  • horizontal restoration coordinates ⁇ u ⁇ ( r ), v ⁇ ( r ) ⁇ play a role similar to the one played by paleo-geographic coordinates ⁇ u ( r ), v ( r ) ⁇ with respect to horizons ⁇ H t : t ⁇ 0 ⁇ 216 of the depositional model provided as input.
  • horizontal restoration coordinates ⁇ u ⁇ ⁇ r ), v ⁇ ( r ) ⁇ may be generated as follows:
  • the restoration vector field R ⁇ (r) represents the field of deformation vectors from the present day (e.g., xyz) space to the restoration (e.g., u ⁇ v ⁇ t ⁇ ) space, e.g., computed from the u ⁇ v ⁇ t ⁇ - transform.
  • an additional restoration time ⁇ n +1 may be added to be associated with the horizontal plane H t n +1 located at a constant altitude z ⁇ n + 1 0 of the sea level.
  • Time ⁇ n +1 may be the present day geological time and, provided that ⁇ n +1 is greater than ⁇ n , any arbitrary value may be chosen for ⁇ n +1 .
  • Geological models are generated using geological or seismic tomography technology.
  • Geological tomography generates an image of the interior subsurface of the Earth based on geological data collected by transmitting a series of incident waves and receiving reflections of those waves across discontinuities in the subsurface.
  • a transmitter may transmit signals, for example, acoustic waves, compression waves or other energy rays or waves, that may travel through subsurface structures.
  • the transmitted signals may become incident signals that are incident to subsurface structures.
  • the incident signals may reflect at various transition zones or geological discontinuities throughout the subsurface structures, such as, faults or horizons.
  • the reflected signals may include seismic events.
  • a receiver may collect data, for example, reflected seismic events.
  • the data may be sent to a modeling mechanism that may include, for example, a data processing mechanism and an imaging mechanism.
  • FIG. 14 is a schematic illustration of a geological tomography technique in which a series of incident rays 111 and reflected rays 121 are propagated through a subsurface region of the Earth 30 to image the subsurface, according to an embodiment of the invention.
  • One or more transmitter(s) located at incident location(s) 60 may emit a series of incident rays 111.
  • Incident rays 111 may include for example a plurality of energy rays related to signal waves, e.g., sonic waves, seismic waves, compression waves, etc.
  • Incident rays 111 may be incident on, and reflect off of, a subsurface structure or surface 90 at a reflection point 50. Multiple reflection points 50 may be identified or imaged or displayed in conjunction to display, for example, a horizon.
  • One or more receiver(s) located at reflected location(s) 65 may receive the reflection rays 121.
  • Reflection rays 121 may be the reflected images of incident rays 111, for example, after reflecting off of image surface 90 at target point 50.
  • the angle of reflection 55 may be the angle between corresponding incident rays 111 and reflected rays 121 at reflection point 50.
  • An incident rays 111 and a corresponding reflected rays 121 may propagate through a cross-section of a subsurface structure 30.
  • Incident rays 111 may reflect off of a subsurface feature 90 at a reflection point 50, for example, a point on an underground horizon, the seafloor, an underground aquifer, etc.
  • One or more processor(s) may reconstitute incident and reflected rays 111 and 121 to generate an image the subsurface 30 using an imaging mechanism.
  • a common reflection angle migration (CRAM) imaging mechanism may image reflection points 50 by aggregating all reflected signals that may correspond to a reflection point, for example, reflected signals that may have the same reflection angle.
  • imaging mechanisms may aggregate reflected signals that may have the same reflection offset (distance between transmitter and receiver), travel time, or other suitable conditions.
  • the processor(s) may compose all of the reflection points 50 to generate an image or model of the present day underground subsurface of the Earth 30.
  • the processor(s) may execute a restoration transformation (e.g., u ⁇ v ⁇ t ⁇ - transform 201) to transform the present day model of subsurface 30 to a restored subsurface image 203 at a restoration time ⁇ .
  • One or more display(s) e.g., 180 of FIG. 15 ) may visualize the present day subsurface image 30 and/or the restored subsurface image 203.
  • FIG. 15 schematically illustrates a system including one or more transmitter(s), one or more receiver(s) and a computing system in accordance with an embodiment of the present invention. Methods disclosed herein may be performed using a system 1505 of Fig. 15 .
  • System 1505 may include one or more transmitter(s) 190, one or more receiver(s) 120, a computing system 130, and a display 180.
  • the aforementioned data e.g., seismic data used to form intermediate data and finally to model subsurface regions, may be ascertained by processing data generated by transmitter 190 and received by receiver 120.
  • Intermediate data may be stored in memory 150 or other storage units.
  • the aforementioned processes described herein may be performed by software 160 being executed by processor 140 manipulating the data.
  • Transmitter 190 may transmit signals, for example, acoustic waves, compression waves or other energy rays or waves, that may travel through subsurface (e.g., below land or sea level) structures.
  • the transmitted signals may become incident signals that are incident to subsurface structures.
  • the incident signals may reflect at various transition zones or geological discontinuities throughout the subsurface structures.
  • the reflected signals may include seismic data.
  • Receiver 120 may accept reflected signal(s) that correspond or relate to incident signals, sent by transmitter 190.
  • Transmitter 190 may transmit output signals.
  • the output of the seismic signals by transmitter 190 may be controlled by a computing system, e.g., computing system 130 or another computing system separate from or internal to transmitter 190.
  • An instruction or command in a computing system may cause transmitter 190 to transmit output signals.
  • the instruction may include directions for signal properties of the transmitted output signals (e.g., such as wavelength and intensity).
  • the instruction to control the output of the seismic signals may be programmed in an external device or program, for example, a computing system, or into transmitter 190 itself.
  • Computing system 130 may include, for example, any suitable processing system, computing system, computing device, processing device, computer, processor, or the like, and may be implemented using any suitable combination of hardware and/or software.
  • Computing system 130 may include for example one or more processor(s) 140, memory 150 and software 160.
  • Data 155 generated by reflected signals, received by receiver 120, may be transferred, for example, to computing system 130.
  • the data may be stored in the receiver 120 as for example digital information and transferred to computing system 130 by uploading, copying or transmitting the digital information.
  • Processor 140 may communicate with computing system 130 via wired or wireless command and execution signals.
  • Memory 150 may include cache memory, long term memory such as a hard drive, and/or external memory, for example, including random access memory (RAM), read only memory (ROM), dynamic RAM (DRAM), synchronous DRAM (SD-RAM), flash memory, volatile memory, non-volatile memory, cache memory, buffer, short term memory unit, long term memory unit, or other suitable memory units or storage units.
  • Memory 150 may store instructions (e.g., software 160) and data 155 to execute embodiments of the aforementioned methods, steps and functionality (e.g., in long term memory, such as a hard drive).
  • Data 155 may include, for example, raw seismic data collected by receiver 120, instructions for building a mesh (e.g., 100), instructions for partitioning a mesh, and instructions for processing the collected data to generate a model, or other instructions or data.
  • Memory 150 may also store instructions to divide and model ⁇ -active faults and ⁇ -inactive faults.
  • Memory 150 may generate and store the aforementioned constraints, restoration transformation (e.g., u ⁇ v ⁇ t ⁇ - transform 201), restoration coordinates (e.g., u ⁇ , v ⁇ , t ⁇ ), a geological-time and paleo-geographic coordinates (e.g., u, v, t ), a model representing a structure when it was originally deposited (e.g., in wvt-space), a model representing a structure at an intermediate restoration time (e.g., in u ⁇ , v ⁇ , t ⁇ -space), and/or a model representing the corresponding present day structure in a current time period (e.g., in xyz-space).
  • restoration transformation e.g., u ⁇ v ⁇ t ⁇ - transform 201
  • restoration coordinates e.g., u ⁇ , v ⁇ , t ⁇
  • Memory 150 may store cells, nodes, voxels, etc., associated with the model and the model mesh. Memory 150 may also store forward and/or reverse u ⁇ , v ⁇ , t ⁇ -transformations to restore present day models (e.g., in xyz-space) to restored models (e.g., in u ⁇ , v ⁇ , t ⁇ -space), and vice versa. Memory 150 may also store the three-dimensional restoration vector fields, which when applied to the nodes of the initial present day model, move the nodes of the initial model to generate one of the plurality of restored models.
  • Applying a restoration vector field to corresponding nodes of the present day model may cause the nodes to "move", “slide”, or “rotate”, thereby transforming modeled geological features represented by nodes and cells of the initial model.
  • Data 155 may also include intermediate data generated by these processes and data to be visualized, such as data representing graphical models to be displayed to a user.
  • Memory 150 may store intermediate data.
  • System 130 may include cache memory which may include data duplicating original values stored elsewhere or computed earlier, where the original data may be relatively more expensive to fetch (e.g., due to longer access time) or to compute, compared to the cost of reading the cache memory.
  • Cache memory may include pages, memory lines, or other suitable structures. Additional or other suitable memory may be used.
  • Computing system 130 may include a computing module having machine-executable instructions.
  • the instructions may include, for example, a data processing mechanism (including, for example, embodiments of methods described herein) and a modeling mechanism. These instructions may be used to cause processor 140 using associated software 160 modules programmed with the instructions to perform the operations described. Alternatively, the operations may be performed by specific hardware that may contain hardwired logic for performing the operations, or by any combination of programmed computer components and custom hardware components.
  • Embodiments of the invention may include an article such as a non-transitory computer or processor readable medium, or a computer or processor storage medium, such as for example a memory, a disk drive, or a USB flash memory, encoding, including or storing instructions, e.g., computer-executable instructions, which when executed by a processor or controller, carry out methods disclosed herein.
  • a non-transitory computer or processor readable medium or a computer or processor storage medium, such as for example a memory, a disk drive, or a USB flash memory
  • encoding including or storing instructions, e.g., computer-executable instructions, which when executed by a processor or controller, carry out methods disclosed herein.
  • Display 180 may display data from transmitter 190, receiver 120, or computing system 130 or any other suitable systems, devices, or programs, for example, an imaging program or a transmitter or receiver tracking device.
  • Display 180 may include one or more inputs or outputs for displaying data from multiple data sources or to multiple displays.
  • display 180 may display visualizations of subsurface models including subsurface features, such as faults, horizons and unconformities, as a present day subsurface image (e.g., 202 ), a restored subsurface image (e.g., 203 ) and/or a depositional model (e.g., 703 ).
  • Display 180 may display one or more present day model(s), depositional model(s), restoration model(s), as well as a series of chronologically sequential restoration models associated with a sequence of respective restoration times (e.g., ⁇ 1 ⁇ ⁇ 2 ⁇ ⁇ 3 ⁇ ⁇ 4 , as shown in FIG. 13 ).
  • the models may be displayed one at a time, two at a time, or many at a time (e.g., the number selected by a user or automatically based on the difference between models or the total number of models).
  • Display 180 may display the models in a sequence of adjacent models, through which a user may scan (e.g., by clicking a 'next' or 'previous' button with a pointing device such as a mouse or by scrolling through the models).
  • Input device(s) 165 may include a keyboard, pointing device (e.g., mouse, trackball, pen, touch screen), or cursor direction keys, for communicating information and command selections to processor 140.
  • Input device 165 may communicate user direction information and command selections to the processor 140.
  • a user may use input device 165 to select one or more preferred models from among the plurality of perturbed models, recategorize faults as ⁇ -active faults and ⁇ -inactive, or edit, add or delete subsurface structures.
  • Processor 140 may include, for example, one or more processors, controllers or central processing units ("CPUs").
  • Software 160 may be stored, for example, in memory 150.
  • Software 160 may include any suitable software, for example, DSI software.
  • Processor 140 may generate a present day subsurface image (e.g., 202 ), a restored subsurface image (e.g., 203 ) and/or a depositional model (e.g., 703 ), for example, using data 155 from memory 150.
  • a model may simulate structural, spatial or geological properties of a subsurface region, such as, porosity or permeability through geological terrains.
  • Processor 140 may initially generate a three dimensional mesh, lattice, grid or collection of nodes (e.g., 100 ) that spans or covers a domain of interest.
  • the domain may cover a portion or entirety of the three-dimensional subsurface region being modeled.
  • Processor 140 may automatically compute the domain to be modeled and the corresponding mesh based on the collected seismic data so that the mesh covers a portion or the entirety of the three-dimensional subsurface region from which geological data is collected (e.g., the studied subsurface region).
  • the domain or mesh may be selected or modified by a user, for example, entering coordinates or highlighting regions of a simulated optional domain or mesh.
  • the user may select a domain or mesh to model a region of the Earth that is greater than a user-selected subsurface distance (e.g., 100 meters) below the Earth's surface, a domain that occurs relative to geological features (e.g., to one side of a known fault or riverbed), or a domain that occurs relative to modeled structures (e.g., between modeled horizons H(t 1 ) and H(t 100 )).
  • Processor 140 may execute software 160 to partition the mesh or domain into a plurality of three-dimensional (3D) cells, columns, or other modeled data (e.g., represented by voxels, pixels, data points, bits and bytes, computer code or functions stored in memory 150 ).
  • 3D three-dimensional
  • the cells or voxels may have hexahedral, tetrahedral, or any other polygonal shapes, and preferably three-dimensional shapes.
  • data may include zero-dimensional nodes, one-dimensional segments, two-dimensional facet and three-dimensional elements of volume, staggered in a three-dimensional space to form three-dimensional data structures, such as cells, columns or voxels.
  • the cells preferably conform to and approximate the orientation of faults and unconformities.
  • Each cell may include faces, edges and/or vertices.
  • Each cell or node may correspond to one or more particles of sediment in the Earth (e.g., a node may include many cubic meters of earth, and thus many particles).
  • Data collected by receiver 120 after the time of deposition in a current or present time period include faults and unconformities that have developed since the original time of deposition, e.g., based on tectonic motion, erosion, or other environmental factors, may disrupt the regular structure of the geological domain. Accordingly, an irregular mesh may be used to model current geological structures, for example, so that at least some faces, edges, or surfaces of cells are oriented parallel to faults and unconformities, and are not intersected thereby.
  • a mesh may be generated based on data collected by receiver 120 , alternatively, a generic mesh may be generated to span the domain and the data collected by receiver 120 may be used to modify the structure thereof.
  • the data collected may be used to generate a set of point values at "sampling point".
  • the values at these points may reorient the nodes or cells of the mesh to generate a model that spatially or otherwise represents the geological data collected from the Earth.
  • Other or different structures, data points, or sequences of steps may be used to process collected geological data to generate a model.
  • the various processes described herein e.g., restoring a geological model using ⁇ -active and ⁇ - inactive faults, or restoring a geological model using a new thickness-preserving constraint
  • Restoration coordinates may be defined at a finite number of nodes or sampling points based on real data corresponding to a subsurface structure, e.g., one or more particles or a volume of particles of Earth. Restoration coordinates may be approximated between nodes to continuously represent the subsurface structure, or alternatively, depending on the resolution in which the data is modeled may represent discrete or periodic subsurface structures, e.g., particles or volumes of Earth that are spaced from each other.
  • the computing system of FIG. 15 may accept the data used in the operations of FIGS. 16 and 17 as for example a set of data generated by tomographic scanning of a subsurface geological region of the Earth as disclosed in reference to FIG. 14 , or such data augmented by another process.
  • the computing system may accept one or more of seismic and well data.
  • the computing device may generate one or more of seismic and well data.
  • Restoration or “intermediate” time ⁇ may refer to a time in the past before the present day and after a time when an oldest or deepest horizon surface in the 3D model was deposited.
  • “Restoration” or “intermediate” transformation or model may refer to a model or image of the surface as it was configured at the “intermediate” time in the past ⁇ .
  • An intermediate horizon may refer to a horizon that was deposited at the "intermediate" time ⁇ , which is located above the deepest horizon and below the shallowest horizon.
  • Time including the present-day, current or present time, the past restoration time ⁇ , and/or the depositional time t, may refer to geological time periods that span a duration of time, such as, periods of thousands or millions of years.
  • “Geological-time” t ( r ) may refer to the time of deposition when a particle of sediment represented by point r was originally deposited in the Earth.
  • the geological-time of the deposition may be replaced, e.g., by any arbitrary monotonic increasing function of the actual geological-time. It is a convention to use an monotonically increasing function, but similarly an arbitrary monotonic decreasing function may be used.
  • the monotonic function may be referred to as the "pseudo-geological-time”.
  • the geological-time of the deposition and restoration time of particles are predicted approximate positions since past configurations can not typically be verified.
  • “Current” or “present day” location for a particle (or data structure representing one or more particles) or subsurface feature may refer to the location of the item in the present time, as it is measured.
  • layers, horizons, faults and unconformities may be curvilinear surfaces which may be for example characterized as follows.
  • Terrain deformed in the neighborhood of a point r in the G-space may occur according to a "minimal deformation" tectonic style when, in this neighborhood:
  • Terrain deformed in the neighborhood of a point r in the G-space may occur according to a "flexural slip" tectonic style when, in this neighborhood:
  • DSI Discrete-Smooth-Interpolation
  • a function f ( x , y , z ) at nodes of a 3D grid or mesh ⁇ (e.g., 100 ), while honoring a given set of constraints.
  • the DSI method allows properties of structures to be modeled by embedding data associated therewith in a (e.g., 3D Euclidean) modeled space.
  • the function f ( x,y,z ) may be defined by values at the nodes of the mesh, ⁇ .
  • the DSI method allows the values of f ( x,y,z ) to be computed at the nodes of the mesh, ⁇ , so that a set of one or more (e.g., linear) constraints are satisfied.
  • DSI generally only applies linear constraints on the model.
  • bold symbols represent vectors or multi-dimensional (e.g., 3D) functions or data structures.
  • any other depositional model visualizing the predicted configuration of each particle, region or layer at its respective time of depositional may be used.
  • embodiments of the invention describe the present day coordinates as xyz, the restoration coordinates as u ⁇ v ⁇ t ⁇ , the depositional coordinates as uvt, the restoration transformation as a u ⁇ v ⁇ t ⁇ -transform, and the depositional transformation as a uvt-transform, any other coordinates or transformations may be used.
  • Embodiments of the invention may manipulate data representations of real-world objects and entities such as underground geological features, including faults and other features.
  • the data may be generated by tomographic scanning, as discussed in reference to FIG. 14 , e.g., received by for example a receiver receiving waves generated e.g., by an air gun or explosives, that may be manipulated and stored, e.g., in memory 150 of FIG. 15 , and data such as images representing underground features may be presented to a user, e.g., as a visualization on display 180 of FIG. 15 .
  • a subsurface image or model may refer to a computer-representation or visualization of actual geological features such as horizons and faults that exist in the real world.
  • Some features when represented in a computing device may be approximations or estimates of a real world feature, or a virtual or idealized feature, such as an idealized horizon as produced in a u ⁇ v ⁇ t ⁇ -transform.
  • a model, or a model representing subsurface features or the location of those features is typically an estimate or a "model", which may approximate or estimate the physical subsurface structure being modeled with more or less accuracy.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Graphics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Geophysics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Architecture (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Geophysics And Detection Of Objects (AREA)

Claims (14)

  1. Rechnerimplementiertes Verfahren zum Wiederherstellen eines 3D-Modells der unterirdischen Geologie der Erde aufgrund einer der Gegenwart zugehörigen Geometrie, die zu einer aktuellen Zeit gemessen wird, wobei eine geologische oder seismische Tomographie der unterirdischen Strukturen verwendet wird, um eine vorhergesagte der Vergangenheit zugehörige Geometrie zu einer vergangenen Wiederherstellungszeit zu erhalten, wobei das Verfahren wie folgt umfasst:
    Empfangen (1610) des 3D-Modells der der Gegenwart zugehörigen gemessenen Geometrie, die ein Netz von Verwerfungen umfasst, wobei eine Verwerfung eine Diskontinuität ist, die Verwerfungsblöcke teilt, die, insoweit wie die Zeit einer modellierten Zeit näher rückt, tangential zur Oberfläche der Verwerfung jeweils in gegenüber liegende Richtungen gleiten;
    Auswählen (1620) einer vergangenen Wiederherstellungszeit τ, die vor der aktuellen Zeit und nach einer Zeit liegt, als eine älteste Horizont-Oberfläche ursprünglich in dem 3D-Modell abgelagert wurde;
    wobei das Verfahren dadurch gekennzeichnet ist, dass es ferner wie folgt umfasst:
    Teilen (1630) des Netzwerks von Verwerfungen in eine Teilmenge τ-aktiver Verwerfungen und eine Teilmenge τ-inaktiver Verwerfungen, wobei eine τ-aktive Verwerfung eine Verwerfung ist, die zu der vergangenen Wiederherstellungszeit τ aktiv ist, und eine τ-inaktive Verwerfung eine Verwerfung ist, die zu der vergangenen Wiederherstellungszeit τ inaktiv ist, wobei eine Verwerfung als τ-aktiv bestimmt wird, wenn die Verwerfung sich mit einem Horizont H τ schneidet, der ursprünglich zu der vergangenen Wiederherstellungszeit τ abgelagert wurde, und eine Verwerfung wird als τ-inaktiv bestimmt wird, wenn die Verwerfung sich nicht mit dem Horizont H τ schneidet, der ursprünglich zu der vergangenen Wiederherstellungszeit τ abgelagert wurde;
    Wiederherstellen (1640) des 3D-Modells aufgrund der der Gegenwart zugehörigen gemessenen Geometrie, um die vorhergesagte der Vergangenheit zugehörige Geometrie zu der vergangenen Wiederherstellungszeit τ zu erhalten, indem wie folgt verfahren wird:
    Modellieren (1650) jeder τ-aktiven Verwerfung, um Endpunkte eines Horizonts Hτ , die in dem der Gegenwart zugehörigen Modell an gegenüber liegenden Seiten der Verwerfung getrennt sind, zu vereinen, sodass sie bei dem wiederhergestellten Modell in derselben Position zusammenlaufen, indem die Endpunkte in eine Richtung, die tangential zur der Oberfläche der τ-aktiven Verwerfung ist, aufeinander zugleiten; und
    Modellieren (1660) jeder τ-inaktiven Verwerfung, um jeweils zugewiesene Punkte an gegenüber liegenden Seiten der Verwerfung beizubehalten; und
    Erzeugen von Bildern des inneren Untergrunds der Erde für die der Vergangenheit zugehörige wiederhergestellte Geologie, die jeweils basierend auf dem 3D-Modell wiederhergestellt wurde.
  2. Verfahren nach Anspruch 1, umfassend Anzeigen einer Visualisierung eines Abbildes der unterirdischen Geologie der Erde, das mit τ-aktiven Verwerfungen und τ-inaktiven Verwerfungen in dem wiederhergestellten Modell zu der vergangenen Zeit τ in Deckung gebracht ist.
  3. Verfahren nach Anspruch 2, umfassend Anzeigen der τ-aktiven Verwerfungen mittels Transluzenz-Niveaus oder in Farben, die von denen der τ-inaktiven Verwerfungen verschieden sind.
  4. Verfahren nach Anspruch 1, wobei eine τ-inaktive Verwerfung eine Oberfläche ist, bei der die Diskontinuität, die von der Verwerfung induziert ist, deaktiviert wurde, um zu verhindern, dass Verwerfungsblöcke in Richtungen gleiten, die, insoweit wie die Zeit der Wiederherstellungszeit τ näher rückt, tangential zu der Oberfläche der Verwerfung sind.
  5. Verfahren nach Anspruch 1, wobei das der Gegenwart zugehörige Modell Iso-Wert-Oberflächen jeder Koordinate einer 3D-Restauration eines der Vergangenheit zugehörigen wiederhergestellten Modells der unterirdischen Geologie zu einer vergangenen Wiederherstellungszeit τ hat, wobei die Iso-Wert-Oberflächen überall auf τ-inaktiven Verwerfungen kontinuierlich und überall auf τ-aktiven Verwerfungen diskontinuierlich sind.
  6. Verfahren nach Anspruch 1, umfassend Ausflachen von Horizont H τ zur Wiederherstellungszeit τ, welcher ursprünglich zu der Zeit τ auf einer im Wesentlichen ebenen Oberfläche einer in etwa konstanten Tiefe abgelagert wurde.
  7. Verfahren nach Anspruch 6, wobei Horizonte, die älter sind als der Horizont H τ, als nicht ebene Oberflächen wiederhergestellt werden.
  8. Verfahren nach Anspruch 1, wobei Wiederherstellen des 3D-Modells zu dem zu einer vergangenen Wiederherstellungszeit τ ein Eliminieren von sämtlichen relativ flächigeren Horizontoberflächen, die ursprünglich nach der vergangenen Wiederherstellungszeit τ abgelagert wurden, umfasst.
  9. Verfahren nach Anspruch 1, wobei eine Verwerfung, die zu einer ersten Wiederherstellungszeit τ' τ-aktiv ist, zu einer zweiten Wiederherstellungszeit τ" τ-inaktiv ist.
  10. Verfahren nach Anspruch 1, umfassend sequenzielles Wiederherstellen des 3D-Models zu einer Sequenz von vergangenen Wiederherstellungszeiten τ1, τ2, ...τn, wobei jede einer Vielzahl der Sequenz von vergangenen Wiederherstellungszeiten eine verschiedene Teilmenge τ-aktiver und τ-inaktiver Verwerfungen hat.
  11. Verfahren nach Anspruch 10, umfassend Abspielen einer sich bewegenden Bildsequenz, bei der das 3D-Modell in einer Vorwärts- oder Rückwärts-Reihenfolge der Sequenz von vergangenen Wiederherstellungszeiten τ1, τ2, ...τn, iterativ wiederherstellt wird, um Änderungen in der unterirdischen Geologie im zeitlichen Verlauf zu visualisieren.
  12. System zum Wiederherstellen eines 3D-Modells der unterirdischen Geologie der Erde aufgrund einer der Gegenwart zugehörigen Geometrie, die zu einer aktuellen Zeit gemessen wird, um eine vorhergesagte der Vergangenheit zugehörige Geometrie zu einer vergangenen Wiederherstellungszeit zu erhalten, wobei das System einen oder mehrere Prozessoren umfasst, die dazu konfiguriert sind, das Verfahren nach einem der vorhergehenden Ansprüche auszuführen.
  13. System nach Anspruch 12, umfassend eine Anzeige, die dazu konfiguriert ist, ein Abbild der unterirdischen Geologie der Erde zu visualisieren, das mit τ-aktiven und τ-inaktiven Verwerfungen in einem wiederhergestellten Modell zu der Zeit τ in Deckung gebracht ist.
  14. System nach Anspruch 13, wobei die Anzeige dazu konfiguriert ist, die τ-aktiven Verwerfungen mittels Transluzenz-Niveaus oder in Farben, die von denen der τ-inaktiven Verwerfungen verschieden sind, zu visualisieren.
EP20150807.4A 2019-01-10 2020-01-08 Abbildung eines unterirdischen geologischen modells zu einer vergangenen wiederherstellungszwischenzeit Active EP3683601B1 (de)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US16/244,544 US10520644B1 (en) 2019-01-10 2019-01-10 Imaging a subsurface geological model at a past intermediate restoration time

Publications (2)

Publication Number Publication Date
EP3683601A1 EP3683601A1 (de) 2020-07-22
EP3683601B1 true EP3683601B1 (de) 2023-11-15

Family

ID=69057460

Family Applications (1)

Application Number Title Priority Date Filing Date
EP20150807.4A Active EP3683601B1 (de) 2019-01-10 2020-01-08 Abbildung eines unterirdischen geologischen modells zu einer vergangenen wiederherstellungszwischenzeit

Country Status (2)

Country Link
US (2) US10520644B1 (de)
EP (1) EP3683601B1 (de)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10918943B2 (en) * 2019-03-18 2021-02-16 Square Enix Ltd. Hexagonal fragmentation of terrain in computer game
US11446580B2 (en) 2019-03-18 2022-09-20 Square Enix Ltd. Rule-based level generation in computer game
EP4193197A1 (de) * 2020-09-08 2023-06-14 MATRIX JVCO LTD trading as AIQ Werkzeug mit künstlicher intelligenz zur geostress- und geostressanalyse
US11892580B2 (en) * 2021-01-13 2024-02-06 Chevron U.S.A. Inc. Decompaction of subsurface region using decompaction velocity
US11921250B2 (en) 2022-03-09 2024-03-05 Saudi Arabian Oil Company Geo-mechanical based determination of sweet spot intervals for hydraulic fracturing stimulation
CN117805900A (zh) * 2022-09-30 2024-04-02 中国石油化工股份有限公司 一种基于优势入射角-频率双域衰减的含气性预测方法

Family Cites Families (206)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3731269A (en) 1971-04-05 1973-05-01 Chevron Res Static corrections for seismic traces by cross-correlation method
US4063213A (en) 1976-06-28 1977-12-13 Texaco Inc. Methods for accurately positioning a seismic energy source while recording seismic data
US4504438A (en) 1981-12-07 1985-03-12 Levy Richard H Method and apparatus for determining the density characteristics of underground earth formations
US4715021A (en) 1984-08-08 1987-12-22 Exxon Production Research Co. Method for filtering and combining seismic data having different spectral characteristics
US4672546A (en) 1984-11-08 1987-06-09 Texas Instruments Incorporated Systems for extracting horizons contained in a 3-D network of turnings and converting the horizons to a 2-D format
US4821164A (en) 1986-07-25 1989-04-11 Stratamodel, Inc. Process for three-dimensional mathematical modeling of underground geologic volumes
FR2632734B1 (fr) 1988-06-10 1990-09-14 Schlumberger Prospection Procede d'etablissement d'un modele stratigraphique du sous-sol a partir d'un profil d'impedance acoustique et d'une section sismique
US4953140A (en) 1989-08-21 1990-08-28 Mobil Oil Corporation Method of subterranean mapping
FR2652180B1 (fr) 1989-09-20 1991-12-27 Mallet Jean Laurent Procede de modelisation d'une surface et dispositif pour sa mise en óoeuvre.
US5056066A (en) 1990-06-25 1991-10-08 Landmark Graphics Corporation Method for attribute tracking in seismic data
US5251184A (en) 1991-07-09 1993-10-05 Landmark Graphics Corporation Method and apparatus for finding horizons in 3D seismic data
GB9214482D0 (en) 1992-07-08 1992-08-19 Armitage Kenneth Sequence property interpretation & risk analysis link
US6005916A (en) 1992-10-14 1999-12-21 Techniscan, Inc. Apparatus and method for imaging with wavefields using inverse scattering techniques
US5537365A (en) 1993-03-30 1996-07-16 Landmark Graphics Corporation Apparatus and method for evaluation of picking horizons in 3-D seismic data
US5384752A (en) 1993-04-06 1995-01-24 Exxon Production Research Company Method for correcting a seismic source pulse waveform
US5570106A (en) 1993-07-28 1996-10-29 Landmark Graphics Corporation Method and apparatus for creating horizons from 3-D seismic data
USRE38229E1 (en) 1994-12-12 2003-08-19 Core Laboratories Global N.V. Method and apparatus for seismic signal processing and exploration
US5930730A (en) 1994-12-12 1999-07-27 Amoco Corporation Method and apparatus for seismic signal processing and exploration
US5594807A (en) 1994-12-22 1997-01-14 Siemens Medical Systems, Inc. System and method for adaptive filtering of images based on similarity between histograms
US5586082A (en) 1995-03-02 1996-12-17 The Trustees Of Columbia University In The City Of New York Method for identifying subsurface fluid migration and drainage pathways in and among oil and gas reservoirs using 3-D and 4-D seismic imaging
US5740342A (en) 1995-04-05 1998-04-14 Western Atlas International, Inc. Method for generating a three-dimensional, locally-unstructured hybrid grid for sloping faults
US5671136A (en) 1995-12-11 1997-09-23 Willhoit, Jr.; Louis E. Process for seismic imaging measurement and evaluation of three-dimensional subterranean common-impedance objects
FR2744224B1 (fr) 1996-01-26 1998-04-17 Inst Francais Du Petrole Methode pour simuler le remplissage d'un bassin sedimentaire
GB9619699D0 (en) 1996-09-20 1996-11-06 Geco Prakla Uk Ltd Seismic sensor units
US6014343A (en) 1996-10-31 2000-01-11 Geoquest Automatic non-artificially extended fault surface based horizon modeling system
US5862513A (en) 1996-11-01 1999-01-19 Western Atlas International, Inc. Systems and methods for forward modeling of well logging tool responses
US6018497A (en) 1997-02-27 2000-01-25 Geoquest Method and apparatus for generating more accurate earth formation grid cell property information for use by a simulator to display more accurate simulation results of the formation near a wellbore
US6106561A (en) 1997-06-23 2000-08-22 Schlumberger Technology Corporation Simulation gridding method and apparatus including a structured areal gridder adapted for use by a reservoir simulator
IT1294492B1 (it) 1997-09-16 1999-04-12 Enel Spa Metodo per la rappresentazione tramite griglia di calcolo semistrutturata di fenomeni fisici estendentisi in un dominio spaziale
US5987388A (en) 1997-12-26 1999-11-16 Atlantic Richfield Company Automated extraction of fault surfaces from 3-D seismic prospecting data
US5995907A (en) 1998-02-05 1999-11-30 Geoquest Seismic signal processing method and apparatus for generating time slice or horizon maps in response to seismic traces and quadrature traces to determine geologic features
JP2002503816A (ja) 1998-02-11 2002-02-05 アナロジック コーポレーション 対象を分類するコンピュータ断層撮影装置および方法
US6847737B1 (en) 1998-03-13 2005-01-25 University Of Houston System Methods for performing DAF data filtering and padding
RU2145100C1 (ru) 1999-02-16 2000-01-27 Миколаевский Эрнест Юлианович Способ поиска, разведки и исследования нефтегазовой залежи
US6950689B1 (en) 1998-08-03 2005-09-27 Boston Scientific Scimed, Inc. Dynamically alterable three-dimensional graphical model of a body region
US6138075A (en) 1998-08-05 2000-10-24 Landmark Graphics Corporation Methods and apparatus for analyzing seismic data
US6018498A (en) 1998-09-02 2000-01-25 Phillips Petroleum Company Automated seismic fault detection and picking
US6388947B1 (en) 1998-09-14 2002-05-14 Tomoseis, Inc. Multi-crosswell profile 3D imaging and method
US6236942B1 (en) 1998-09-15 2001-05-22 Scientific Prediction Incorporated System and method for delineating spatially dependent objects, such as hydrocarbon accumulations from seismic data
US6820043B2 (en) 1998-11-02 2004-11-16 T-Surf Corporation Modeling continuous properties on discontinuous surfaces
US6278949B1 (en) 1998-11-25 2001-08-21 M. Aftab Alam Method for multi-attribute identification of structure and stratigraphy in a volume of seismic data
US6246963B1 (en) 1999-01-29 2001-06-12 Timothy A. Cross Method for predicting stratigraphy
US6754588B2 (en) 1999-01-29 2004-06-22 Platte River Associates, Inc. Method of predicting three-dimensional stratigraphy using inverse optimization techniques
US6151555A (en) 1999-03-09 2000-11-21 Schlumberger Technology Corporation Seismic signal processing method and apparatus for generating a cube of variance values
US6169959B1 (en) 1999-05-27 2001-01-02 Baker Hughes Incorporated Method of predicting kinematics for surface multiples
DK1058131T3 (da) 1999-06-03 2004-09-20 Jason Geosystems B V Fremgangsmåde til bestemmelse af rumlige ændringer i underjordisk struktur, stratigrafi, litologi og fluidum-indhold og reduktion af seismisk stöj
DK1094338T3 (da) 1999-10-22 2006-12-27 Jason Geosystems B V Fremgangsmåde til vurdering af elasticitet og sammensætningsparametre fra seismiske og ekkoakustiske data
US6480790B1 (en) 1999-10-29 2002-11-12 Exxonmobil Upstream Research Company Process for constructing three-dimensional geologic models having adjustable geologic interfaces
FR2800881B1 (fr) 1999-11-05 2001-11-30 Inst Francais Du Petrole Methode pour realiser en 3d un modele cinematique de deformation d'un bassin sedimentaire
FR2802324B1 (fr) 1999-12-10 2004-07-23 Inst Francais Du Petrole Methode pour generer un maillage sur une formation heterogene traversee par une ou plusieurs discontinuites geometriques dans le but de realiser des simulations
MY125603A (en) 2000-02-25 2006-08-30 Shell Int Research Processing seismic data
US6370491B1 (en) 2000-04-04 2002-04-09 Conoco, Inc. Method of modeling of faulting and fracturing in the earth
FR2808336B1 (fr) 2000-04-26 2002-06-07 Elf Exploration Prod Methode d'interpretation chrono-stratigraphique d'une section ou d'un bloc sismique
US7369973B2 (en) 2000-06-29 2008-05-06 Object Reservoir, Inc. Method and system for representing reservoir systems
AU2001271629A1 (en) 2000-06-29 2002-01-14 Object Reservoir, Inc. Method and system for high-resolution modeling of a well bore in a hydrocarbon reservoir
GC0000235A (en) 2000-08-09 2006-03-29 Shell Int Research Processing an image
EP2362346B1 (de) 2001-04-18 2013-07-10 Landmark Graphics Corporation Verfahren und Computerprogrammprodukt zur Volumendarstellung
US20050114831A1 (en) 2001-04-18 2005-05-26 Andres Callegari Volume body renderer
EP1397703B1 (de) 2001-06-06 2008-02-27 Schlumberger Oilfield Assistance Limited Automatisches system zur modellierung von gestörten mehrwertigen horizonten
US6823266B2 (en) 2001-06-20 2004-11-23 Exxonmobil Upstream Research Company Method for performing object-based connectivity analysis in 3-D seismic data volumes
US6850845B2 (en) 2001-07-20 2005-02-01 Tracy Joseph Stark System for multi-dimensional data analysis
WO2003009004A1 (en) 2001-07-20 2003-01-30 Tracy Joseph Stark System for information extraction from geologic time volumes
US6853922B2 (en) 2001-07-20 2005-02-08 Tracy Joseph Stark System for information extraction from geologic time volumes
US7991717B1 (en) 2001-09-10 2011-08-02 Bush Ronald R Optimal cessation of training and assessment of accuracy in a given class of neural networks
US7187794B2 (en) 2001-10-18 2007-03-06 Research Foundation Of State University Of New York Noise treatment of low-dose computed tomography projections and images
GB0125713D0 (en) 2001-10-26 2001-12-19 Statoil Asa Method of combining spatial models
US6904169B2 (en) 2001-11-13 2005-06-07 Nokia Corporation Method and system for improving color images
FR2833384B1 (fr) 2001-12-10 2004-04-02 Tsurf Procede, dispositif et produit programme de modelisation tridimensionnelle d'un volume geologique
US7248259B2 (en) 2001-12-12 2007-07-24 Technoguide As Three dimensional geological model construction
US7069149B2 (en) 2001-12-14 2006-06-27 Chevron U.S.A. Inc. Process for interpreting faults from a fault-enhanced 3-dimensional seismic attribute volume
US7523024B2 (en) 2002-05-17 2009-04-21 Schlumberger Technology Corporation Modeling geologic objects in faulted formations
US7227983B1 (en) 2002-05-30 2007-06-05 The Regents Of The University Of California Automated macromolecular crystal detection system and method
US7744534B2 (en) 2002-06-07 2010-06-29 Verathon Inc. 3D ultrasound-based instrument for non-invasive measurement of amniotic fluid volume
US6791900B2 (en) 2002-06-13 2004-09-14 Exxonmobil Upstream Research Company Method of calculating a throw volume for quantitative fault analysis
US7184991B1 (en) 2002-07-12 2007-02-27 Chroma Energy, Inc. Pattern recognition applied to oil exploration and production
US7024021B2 (en) 2002-09-26 2006-04-04 Exxonmobil Upstream Research Company Method for performing stratigraphically-based seed detection in a 3-D seismic data volume
US6778909B1 (en) 2002-10-05 2004-08-17 3Dgeo Development, Inc. Seismic data processing systems and methods
AU2003275021A1 (en) 2002-10-18 2004-05-13 Exxonmobil Upstream Research Company A method for rapid fault interpretation of fault surfaces generated to fit three-dimensional seismic discontinuity data
WO2004044615A2 (en) 2002-11-09 2004-05-27 Geoenergy, Inc. Method and apparatus for seismic feature extraction
ATE445855T1 (de) 2002-12-13 2009-10-15 Schlumberger Technology Bv Verfahren und vorrichtung zur verbesserten tiefenanpassung von bohrlochbildern oder probenbildern
US7181380B2 (en) 2002-12-20 2007-02-20 Geomechanics International, Inc. System and process for optimal selection of hydrocarbon well completion type and design
FR2849211B1 (fr) 2002-12-20 2005-03-11 Inst Francais Du Petrole Methode de modelisation pour constituer un modele simulant le remplissage multilithologique d'un bassin sedimentaire
GB2396448B (en) 2002-12-21 2005-03-02 Schlumberger Holdings System and method for representing and processing and modeling subterranean surfaces
US6950751B2 (en) 2003-03-31 2005-09-27 Conocophillps Company Method and apparatus for the assimilation and visualization of information from 3D data volumes
US7248539B2 (en) 2003-04-10 2007-07-24 Schlumberger Technology Corporation Extrema classification
US6889142B2 (en) 2003-05-30 2005-05-03 Pgs Exploration (Uk) Limited Method of correcting for time shifts in seismic data resulting from azimuthal variation
DE10325738B3 (de) 2003-06-06 2004-12-02 Atlas Elektronik Gmbh Verfahren zur Generierung eines dreidimensionalen Geländemodells
US7298376B2 (en) 2003-07-28 2007-11-20 Landmark Graphics Corporation System and method for real-time co-rendering of multiple attributes
CA2455810C (en) 2004-01-23 2014-08-26 Tracy Joseph Stark System for information extraction from geologic time volumes
US7657414B2 (en) 2005-02-23 2010-02-02 M-I L.L.C. Three-dimensional wellbore visualization system for hydraulics analyses
FR2868848B1 (fr) 2004-04-07 2006-06-23 Earth Decision Sciences Sa Procede, dispositif et produit-programme de simulation stochastique de paleo chenaux ou de paleo structures geologiques
US7711532B2 (en) 2004-06-02 2010-05-04 Paradigm France Method for building a three dimensional cellular partition of a geological domain
RU2321064C2 (ru) 2004-06-03 2008-03-27 Мурманский государственный технический университет Способ построения обратимой трехмерной гидродинамической модели земли, калибруемой в реальном времени в процессе бурения
FR2871897B1 (fr) 2004-06-21 2006-08-11 Inst Francais Du Petrole Methode pour deformer une image sismique pour interpretation amelioree
EP1766439A2 (de) 2004-07-01 2007-03-28 ExxonMobil Upstream Research Company Auf hydrodynamik basierendes geologisches gridding-modellieren (hydro-gridding)
CA2616871A1 (en) 2004-07-30 2006-02-02 Algolith Inc. Apparatus and method for adaptive 3d noise reduction
US20060047429A1 (en) 2004-08-24 2006-03-02 Adams Steven L Method of estimating geological formation depths by converting interpreted seismic horizons from the time domain to the depth domain
US7925481B2 (en) 2004-09-10 2011-04-12 Van Wagoner John C Geologic models of subsurface sedimentary volumes
EA010085B1 (ru) 2004-09-10 2008-06-30 Эксонмобил Апстрим Рисерч Компани Способ оценивания свойств осадочного бассейна путем численного моделирования процессов осадконакопления
US7554883B2 (en) 2004-10-11 2009-06-30 Landmark Graphics Corporation Fault filter for seismic discontinuity data
EP1815272B1 (de) 2004-11-24 2021-07-28 Paradigm Geophysical Ltd. System und verfahren zur fehleridentifikation
AU2006214069B2 (en) 2005-02-18 2010-12-23 Bp Corporation North America Inc. System and method for using time-distance characteristics in acquisition, processing and imaging of t-CSEM data
WO2006097911A1 (en) 2005-03-17 2006-09-21 Algotec Systems Ltd. Bone segmentation
EP1869502A1 (de) 2005-04-15 2007-12-26 Shell Internationale Research Maatschappij B.V. Identifizieren einer belastungsanomalie in einer suboberflächenregion
EP1880327B1 (de) 2005-04-19 2021-09-08 D.E. Shaw Research, LLC Verfahren zur rechenlastverteilung in der berechnung von abstand-begrenzten wechselwirkungen zwischen teilchen
US7480205B2 (en) 2005-04-20 2009-01-20 Landmark Graphics Corporation 3D fast fault restoration
US7280918B2 (en) 2005-08-08 2007-10-09 Knowledge Systems, Inc. Method and system for combining seismic data and basin modeling
US7126340B1 (en) 2005-09-30 2006-10-24 Saudi Arabian Oil Company Method to characterize microfractured hydrocarbon reservoirs by artificially induced anisotropy of magnetic susceptibility
US7512529B2 (en) 2005-10-26 2009-03-31 Roxar Software Solutions A/S Analysis and characterization of fault networks
CA2625429C (en) 2005-11-03 2014-07-22 Saudi Arabian Oil Company Continuous reservoir monitoring for fluid pathways using 3d microseismic data
US7660481B2 (en) 2005-11-17 2010-02-09 Vital Images, Inc. Image enhancement using anisotropic noise filtering
US7486589B2 (en) 2006-02-09 2009-02-03 Schlumberger Technology Corporation Methods and apparatus for predicting the hydrocarbon production of a well location
US8605975B2 (en) 2006-02-13 2013-12-10 The University Of Chicago Image reconstruction from limited or incomplete data
US7539079B2 (en) 2006-03-29 2009-05-26 Pgs Geophysical As System and method for determining positions of towed marine source-array elements
US7565276B2 (en) 2006-04-05 2009-07-21 Seoul National University Industry Foundation Method of simulating detailed movements of fluids using derivative particles
EP2395375A3 (de) 2006-06-21 2012-04-11 Terraspark Geosciences, LLC Extraktion von Ablagerungssystemen
CN101583911B (zh) 2006-09-01 2013-01-02 哈利伯顿兰德马克绘图公司 用于对波形体进行成像的系统和方法
US20080088621A1 (en) 2006-10-11 2008-04-17 Jean-Jacques Grimaud Follower method for three dimensional images
US8120991B2 (en) 2006-11-03 2012-02-21 Paradigm Geophysical (Luxembourg) S.A.R.L. System and method for full azimuth angle domain imaging in reduced dimensional coordinate systems
FR2909185B1 (fr) 2006-11-27 2009-01-09 Inst Francais Du Petrole Methode d'interpretation stratigraphique d'images sismiques
GB2444506C (en) 2006-12-06 2010-01-06 Schlumberger Holdings Processing of stratigraphic data
US20090230969A1 (en) 2007-02-19 2009-09-17 Hall David R Downhole Acoustic Receiver with Canceling Element
CA2680021A1 (en) 2007-03-05 2008-09-12 Paradigm Geophysical (Luxembourg) S.A.R.L. Model-based time-preserving tomography
US20080232694A1 (en) 2007-03-21 2008-09-25 Peter Sulatycke Fast imaging data classification method and apparatus
FR2914434B1 (fr) 2007-03-30 2009-05-22 Inst Francais Du Petrole Methode de calage d'historique d'un modele geologique par modification graduelle des proportions des facies lithologiques
US8150663B2 (en) 2007-03-30 2012-04-03 Paradigm Geophysical (Luxembourg) S.A.R.L. Partitioning algorithm for building a stratigraphic grid
CA2690129C (en) 2007-06-07 2015-08-11 Paradigm Geophysical Corp. Device and method for displaying full azimuth angle domain image data
US8538702B2 (en) 2007-07-16 2013-09-17 Exxonmobil Upstream Research Company Geologic features from curvelet based seismic attributes
US8768672B2 (en) 2007-08-24 2014-07-01 ExxonMobil. Upstream Research Company Method for predicting time-lapse seismic timeshifts by computer simulation
JP5371221B2 (ja) 2007-09-11 2013-12-18 プロメテック・ソフトウェア株式会社 粒子法シミュレーションのためのスライスデータ構造、およびスライスデータ構造を利用した粒子法シミュレーションのｇｐｕへの実装方法
WO2009039422A1 (en) 2007-09-21 2009-03-26 Headwave, Inc. Seismic data processing and visualization
US8194068B1 (en) 2007-10-17 2012-06-05 Sandia Corporation Hexahedral finite element mesh coarsening using pillowing technique
US8140310B2 (en) 2007-11-01 2012-03-20 Schlumberger Technology Corporation Reservoir fracture simulation
GB0722469D0 (en) 2007-11-16 2007-12-27 Statoil Asa Forming a geological model
US7702463B2 (en) 2007-12-12 2010-04-20 Landmark Graphics Corporation, A Halliburton Company Systems and methods for enhancing a seismic data image
EP2232301A1 (de) 2007-12-21 2010-09-29 ExxonMobil Upstream Research Company Modellierung in sedimentären basins
US20090192717A1 (en) 2008-01-24 2009-07-30 Terje Iversen Methods and apparatuses for dynamic filtering of geometric primitives in 3d space
US7769546B2 (en) 2008-02-11 2010-08-03 Chevron U.S.A. Inc. Method for indexing a subsurface volume for the purpose of inferring geologic information
US8803878B2 (en) 2008-03-28 2014-08-12 Schlumberger Technology Corporation Visualizing region growing in three dimensional voxel volumes
CA2721008A1 (en) 2008-04-11 2009-10-15 Terraspark Geosciences, Llc Visulation of geologic features using data representations thereof
FR2930350B1 (fr) 2008-04-17 2011-07-15 Inst Francais Du Petrole Procede pour rechercher des hydrocarbures dans un bassin geologiquement complexe,au moyen d'une modelisation de bassin
CN102047294B (zh) 2008-06-06 2013-10-30 兰德马克绘图国际公司，哈里伯顿公司 对表示网格体的几何不规则网格数据的三维体成像的系统及方法
US9086507B2 (en) 2008-08-18 2015-07-21 Westerngeco L.L.C. Determining characteristics of a subterranean body using pressure data and seismic data
US9146329B2 (en) 2008-12-17 2015-09-29 Exxonmobil Upstream Research Company System and method for reconstruction of time-lapse data
KR101194605B1 (ko) 2008-12-22 2012-10-25 한국전자통신연구원 시간 연속적 텍스쳐 합성 장치 및 방법
WO2010082969A1 (en) 2009-01-13 2010-07-22 Exxonmobil Upstream Research Company Methods and systems to volumetrically conceptualize hydrocarbon plays
US8599644B2 (en) 2009-02-04 2013-12-03 Schlumberger Technology Corporation Velocity models for a single well and for a set of wells
US8165817B2 (en) 2009-03-09 2012-04-24 Schlumberger Technology Corporation Method for integrating reservoir charge modeling and downhole fluid analysis
US8085620B2 (en) 2009-03-27 2011-12-27 Westerngeco L.L.C. Determining a position of a survey receiver in a body of water
US20100256964A1 (en) 2009-04-07 2010-10-07 Schlumberger Technology Corporation System and technique to quantify a fracture system
US8219370B1 (en) 2009-05-20 2012-07-10 Adobe Systems Incorporated Simulation of shallow viscoelastic flows
US9418182B2 (en) 2009-06-01 2016-08-16 Paradigm Sciences Ltd. Systems and methods for building axes, co-axes and paleo-geographic coordinates related to a stratified geological volume
US9536022B1 (en) 2009-06-01 2017-01-03 Paradigm Sciences Ltd. Systems and methods for modeling faults in the subsurface
US8711140B1 (en) 2009-06-01 2014-04-29 Paradigm Sciences Ltd. Systems and methods for building axes, co-axes and paleo-geographic coordinates related to a stratified geological volume
US8600708B1 (en) 2009-06-01 2013-12-03 Paradigm Sciences Ltd. Systems and processes for building multiple equiprobable coherent geometrical models of the subsurface
FR2948215B1 (fr) 2009-07-16 2011-06-24 Inst Francais Du Petrole Methode pour generer un maillage hexa-dominant d'un milieu souterrain faille
EP2473947A4 (de) 2009-08-31 2015-11-04 Parx Ltd Vollautomatisches parksystem
US8655632B2 (en) 2009-09-03 2014-02-18 Schlumberger Technology Corporation Gridless geological modeling
US20110071799A1 (en) 2009-09-21 2011-03-24 Per Arne Slotte Grid models
US8743115B1 (en) 2009-10-23 2014-06-03 Paradigm Sciences Ltd. Systems and methods for coordinated editing of seismic data in dual model
EP2317348B1 (de) 2009-10-30 2014-05-21 Services Pétroliers Schlumberger Verfahren für den Aufbau eines Ablagerungsraumes entsprechend einem geologischen Bereich
US8437962B2 (en) 2009-11-25 2013-05-07 Halliburton Energy Services, Inc. Generating probabilistic information on subterranean fractures
US9176245B2 (en) 2009-11-25 2015-11-03 Halliburton Energy Services, Inc. Refining information on subterranean fractures
US20130144532A1 (en) 2009-12-21 2013-06-06 Schlumberger Technology Corporation Identification of reservoir geometry from microseismic event clouds
US8274859B2 (en) 2010-02-22 2012-09-25 Landmark Graphics Corporation Systems and methods for modeling 3D geological structures
US20120022837A1 (en) 2010-05-18 2012-01-26 Schlumberger Technology Corporation Smoothing Of Stair-Stepped Geometry In Grids
US8635052B1 (en) 2010-06-07 2014-01-21 Paradigm Ltd. Method and appartus for transforming a stratigraphic grid
US20110310101A1 (en) 2010-06-22 2011-12-22 Schlumberger Technology Corporation Pillar grid conversion
US8463586B2 (en) 2010-06-22 2013-06-11 Saudi Arabian Oil Company Machine, program product, and computer-implemented method to simulate reservoirs as 2.5D unstructured grids
FR2962835B1 (fr) 2010-07-16 2013-07-12 IFP Energies Nouvelles Methode pour generer un maillage hexa-dominant d'un bassin geometriquement complexe
US8731875B2 (en) 2010-08-13 2014-05-20 Exxonmobil Upstream Research Company System and method for providing data corresponding to physical objects
WO2012078217A1 (en) 2010-12-08 2012-06-14 Exxonmobil Upstream Research Company Constructing geologic models from geologic concepts
GB201103642D0 (en) 2011-03-03 2011-04-13 Univ Bradford Methods and systems for detection of liquid surface fluctuations
WO2012125558A2 (en) 2011-03-11 2012-09-20 Schlumberger Technology Corporation System and method for performing microseismic fracture operations
US9442205B2 (en) 2011-03-23 2016-09-13 Global Ambient Seismic, Inc. Method for assessing the effectiveness of modifying transmissive networks of natural reservoirs
US8965745B2 (en) 2011-04-14 2015-02-24 Schlumberger Technology Corporation Grid from depositional space
US9201993B2 (en) 2011-05-11 2015-12-01 Apple Inc. Goal-driven search of a stochastic process using reduced sets of simulation points
CA2743611C (en) 2011-06-15 2017-03-14 Engineering Seismology Group Canada Inc. Methods and systems for monitoring and modeling hydraulic fracturing of a reservoir field
WO2013015764A1 (en) 2011-07-22 2013-01-31 Landmark Graphics Corporation Mapping geologic features
FR2979152B1 (fr) 2011-08-17 2013-08-23 IFP Energies Nouvelles Methode pour construire un modele geologique comportant une mise en position de depot d'unites stratigraphiques.
CA2842398C (en) 2011-08-23 2017-11-28 Exxonmobil Upstream Research Company Estimating fracture dimensions from microseismic data
BR112014010082A8 (pt) 2011-10-28 2017-06-20 Landmark Graphics Corp método, sistema de computador e meio não transitório lido por computador armazenando um programa
US20130151215A1 (en) 2011-12-12 2013-06-13 Schlumberger Technology Corporation Relaxed constraint delaunay method for discretizing fractured media
US8744825B2 (en) 2012-02-13 2014-06-03 Livermore Software Technology Corp. Element refinement methods and systems in arbitrary lagrangian-eulerian (ALE) based finite element analysis
US20150009215A1 (en) 2012-02-17 2015-01-08 Schlumberger Technology Corporation Generating a 3d image for geological modeling
US10114134B2 (en) 2012-03-02 2018-10-30 Emerson Paradigm Holding Llc Systems and methods for generating a geological model honoring horizons and faults
FR2987903B1 (fr) 2012-03-09 2014-05-09 Schlumberger Services Petrol Structures de failles geologiques contenant des non-conformites.
US9759826B2 (en) 2012-04-03 2017-09-12 Paradigm Sciences Ltd. System and method for generating an implicit model of geological horizons
EP2901363A4 (de) 2012-09-28 2016-06-01 Exxonmobil Upstream Res Co Fehlerentfernung bei geologischen modellen
US9390555B2 (en) 2012-11-09 2016-07-12 International Business Machines Corporation Method to assess the impact of existing fractures and faults for reservoir management
US20140136171A1 (en) 2012-11-13 2014-05-15 Chevron U.S.A. Inc. Unstructured Grids For Modeling Reservoirs
US10012055B2 (en) 2013-01-24 2018-07-03 Schlumberger Technology Corporation Analysis of surface networks for fluids
EP2778725B1 (de) 2013-03-15 2018-07-18 Emerson Paradigm Holding LLC Systeme und Verfahren zum Aufbauen sedimentärer Attribute
EP2869096B1 (de) 2013-10-29 2019-12-04 Emerson Paradigm Holding LLC Systeme und verfahren zur multiskalaren vermaschung von geologischer zeitmodellierung
US10808501B2 (en) 2014-03-17 2020-10-20 Saudi Arabian Oil Company Modeling intersecting faults and complex wellbores in reservoir simulation
US10422923B2 (en) 2014-03-28 2019-09-24 Emerson Paradigm Holding Llc Systems and methods for modeling fracture networks in reservoir volumes from microseismic events
WO2016070073A1 (en) 2014-10-31 2016-05-06 Exxonmobil Upstream Research Company Managing discontinuities in geologic models
CA2963092C (en) 2014-10-31 2021-07-06 Exxonmobil Upstream Research Company Methods to handle discontinuity in constructing design space for faulted subsurface model using moving least squares
WO2016076847A1 (en) 2014-11-12 2016-05-19 Halliburton Energy Services, Inc. Reservoir mesh creation using extended anisotropic, geometry-adaptive refinement of polyhedra
US20160298427A1 (en) * 2015-04-09 2016-10-13 Schlumberger Technology Corporation Continuum sedimentary basin modeling using particle dynamics simulations
US9697621B2 (en) 2015-06-07 2017-07-04 Apple Inc. Unwrapping 3D meshes into low-distortion texture maps
US9690002B2 (en) 2015-06-18 2017-06-27 Paradigm Sciences Ltd. Device, system and method for geological-time refinement
US10400590B1 (en) 2015-10-16 2019-09-03 Emerson Paradigm Holding Llc Method and system for determining a distribution of rock types in geological cells around a wellbore
US20170184760A1 (en) 2015-12-23 2017-06-29 Paradigm Sciences Ltd. Device, system and method for a structure and stratigraphy preserving transformation of a geological model
US10685482B2 (en) * 2016-05-13 2020-06-16 Chevron U.S.A. Inc. System and method for 3D restoration of complex subsurface models
US10466388B2 (en) 2016-09-07 2019-11-05 Emerson Paradigm Holding Llc System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles

Also Published As

Publication number Publication date
EP3683601A1 (de) 2020-07-22
US10705254B1 (en) 2020-07-07
US10520644B1 (en) 2019-12-31
US20200225382A1 (en) 2020-07-16

Similar Documents

Publication Publication Date Title
EP3683601B1 (de) Abbildung eines unterirdischen geologischen modells zu einer vergangenen wiederherstellungszwischenzeit
US10598819B2 (en) Systems and methods to build sedimentary attributes
US10310137B1 (en) Systems and methods for building axes, co-axes and paleo-geographic coordinates related to a stratified geological volume
US11156744B2 (en) Imaging a subsurface geological model at a past intermediate restoration time
US9524590B1 (en) Systems and methods for building axes, co-axes and paleo-geographic coordinates related to a stratified geological volume
US10338252B1 (en) Systems and processes for building multiple equiprobable coherent geometrical models of the subsurface
US10331817B1 (en) Systems and methods for modeling faults in the subsurface
AU2011339017B2 (en) Constructing geologic models from geologic concepts
EP3293552B1 (de) System und verfahren zum editieren von geologischen modellen durch wechseln zwischen volumenbasierten modellen und oberflächenbasierten strukturmodellen mit stratigrafischen faserbündeln
US10995592B2 (en) Method and system for analyzing the uncertainty of subsurface model
Cacace et al. MeshIt—a software for three dimensional volumetric meshing of complex faulted reservoirs
Galera et al. Construction of coherent 3D geological blocks
Zhu et al. Coupled modeling between geological structure fields and property parameter fields in 3D engineering geological space
EP2778724B1 (de) Verfahren und Systeme zum Aufbau von Achsen, Co-Achsen und paleo-geographischen Koordinaten im Zusammenhang mit einem geschichteten geologischen Volumen
Ceccato et al. Constraints upon fault zone properties by combined structural analysis of virtual outcrop models and discrete fracture network modelling
WO2016001697A1 (en) Systems and methods for geologic surface reconstruction using implicit functions
Souche et al. Depositional space: Construction and applications to facies and petrophysical property simulations
EP3875995A1 (de) Abbildung eines unterirdischen geologischen modells zu einer vergangenen zwischenrestaurationszeit
Agi et al. Impact of geological interpretation on reservoir 3D static model: Workflow, methodology approach and delivery process
Westerteiger et al. Interactive Retro-Deformation of Terrain for Reconstructing 3D Fault Displacements
Mallet et al. 3d restoration of sedimentary terrains: The geochron approach

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION HAS BEEN PUBLISHED

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20200811

RBV Designated contracting states (corrected)

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

17Q First examination report despatched

Effective date: 20220125

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: GRANT OF PATENT IS INTENDED

INTG Intention to grant announced

Effective date: 20230615

P01 Opt-out of the competence of the unified patent court (upc) registered

Effective date: 20230601

RIN1 Information on inventor provided before grant (corrected)

Inventor name: TERTOIS, ANNE-LAURE

Inventor name: MALLET, JEAN-LAURENT

P02 Opt-out of the competence of the unified patent court (upc) changed

Effective date: 20230605

RAP3 Party data changed (applicant data changed or rights of an application transferred)

Owner name: MALLET, JEAN-LAURENT

Owner name: ASPEN PARADIGM HOLDING LLC

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE PATENT HAS BEEN GRANTED

RAP3 Party data changed (applicant data changed or rights of an application transferred)

Owner name: MALLET, JEAN-LAURENT

Owner name: ASPEN PARADIGM HOLDING LLC

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

Ref country code: GB

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602020020913

Country of ref document: DE

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20231228

Year of fee payment: 5

REG Reference to a national code

Ref country code: NL

Ref legal event code: FP

REG Reference to a national code

Ref country code: NO

Ref legal event code: T2

Effective date: 20231115

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: NL

Payment date: 20240112

Year of fee payment: 5

REG Reference to a national code

Ref country code: LT

Ref legal event code: MG9D

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20240216

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20240315

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20231115

REG Reference to a national code

Ref country code: AT

Ref legal event code: MK05

Ref document number: 1632276

Country of ref document: AT

Kind code of ref document: T

Effective date: 20231115

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20231115

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: ES

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20231115

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20231115

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20240315

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20240216

Ref country code: ES

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20231115

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20240215

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20231115

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20240315