WO2019231681A1 - Inverse stratigraphic modeling using a hybrid linear and nonlinear algorithm - Google Patents

Inverse stratigraphic modeling using a hybrid linear and nonlinear algorithm Download PDF

Info

Publication number
WO2019231681A1
WO2019231681A1 PCT/US2019/032489 US2019032489W WO2019231681A1 WO 2019231681 A1 WO2019231681 A1 WO 2019231681A1 US 2019032489 W US2019032489 W US 2019032489W WO 2019231681 A1 WO2019231681 A1 WO 2019231681A1
Authority
WO
WIPO (PCT)
Prior art keywords
optimized
computer
components
input parameters
different grain
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.)
Ceased
Application number
PCT/US2019/032489
Other languages
English (en)
French (fr)
Inventor
Ming Zhang
Mokhles Mustapha Mezghani
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.)
Saudi Arabian Oil Co
Aramco Services Co
Original Assignee
Saudi Arabian Oil Co
Aramco Services Co
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 Saudi Arabian Oil Co, Aramco Services Co filed Critical Saudi Arabian Oil Co
Priority to CA3101976A priority Critical patent/CA3101976A1/en
Priority to JP2020566913A priority patent/JP2021526634A/ja
Priority to EP19728231.2A priority patent/EP3803470A1/en
Priority to CN201980036566.0A priority patent/CN112689778A/zh
Publication of WO2019231681A1 publication Critical patent/WO2019231681A1/en
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • 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/282—Application of seismic models, synthetic seismograms
    • 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
    • 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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V20/00—Geomodelling in general
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12—Simultaneous equations, e.g. systems of linear equations
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • 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/61—Analysis by combining or comparing a seismic data set with other data
    • G01V2210/614—Synthetically generated data
    • 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/61—Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616—Data from specific type of measurement
    • G01V2210/6161—Seismic or acoustic, e.g. land or sea measurements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/66—Subsurface modeling
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/66—Subsurface modeling
    • G01V2210/661—Model from sedimentation process modeling, e.g. from first principles

Definitions

  • a forward stratigraphic model is a digital representation of detailed internal geometry and rock properties of a subsurface earth volume, such as a petroleum reservoir or a sediment-filled basin, and is used as a tool to investigate formation of stratigraphic or sedimentary geologic layers on geological timescales.
  • Forward stratigraphic models provide geologic reservoir simulations used, for example, to select locations for new hydrocarbon wells, estimate hydrocarbon reserves, and plan hydrocarbon-reservoir-development strategies.
  • Three essential prerequisites are needed to consider stratigraphic forward simulation to be successful. First, a stratigraphic formation process recorded by outcrop, seismic, or other hard data must be well understood. Second, a forward model must be able to describe a stratigraphic process- response system. Third, an accurate estimation of many input parameters must be made. However, satisfying all three prerequisites has proven to be technically challenging.
  • the present disclosure describes inverse stratigraphic modeling using a hybrid linear and nonlinear algorithm.
  • a defined scope value is selected for each of a plurality of hydrodynamic input parameters.
  • a simulated topographical result is generated using the selected scope values and a forward model.
  • a detailed seismic interpretation is generated to represent specific seismic features or observed topography.
  • a calculated a misfit value representing a distance between the simulated topographical result and a detailed seismic interpretation is minimized.
  • An estimated optimized sand ratio and optimized hydrodynamic input parameters are generated.
  • a genetic algorithm is used to determine a proportion of each grain size in the estimated optimized sand ratio.
  • a misfit value is used that is calculated from thickness and porosity data extracted from well data and a simulation result generated by the forward model to generate optimized components of different grain sizes. Optimized hydrodynamic input parameters and optimized components of different grain sizes are generated.
  • Implementations of the described subject matter can be implemented using a computer- implemented method; a non-transitory, computer-readable medium storing computer- readable instructions to perform the computer-implemented method; and a computer- implemented system comprising one or more computer memory devices interoperably coupled with one or more computers and having tangible, non-transitory, machine- readable media storing instructions that, when executed by the one or more computers, perform the computer-implemented method/the computer-readable instructions stored on the non-transitory, computer-readable medium.
  • estimating and tuning values of the many required input parameters necessary for a successful forward stratigraphic model is a typically a manual process performed by geological engineers.
  • the manual estimation and tuning can be very time consuming and results are often inconsistent between different geological engineers.
  • the described methodology permits an automatic solution for calibrating a forward stratigraphic model with multiscale prior observation data (for example, well data (from well logs and well cores) and seismic data), which is more accurate and consistent than conventional methods of generating forward stratigraphic models.
  • the generated forward stratigraphic models can be used to improve overall accuracy of geological predictions.
  • the described methodology can be used by a real-time computing system to dynamically control, or direct control of, tangible equipment based on generated output data, such as a forward stratigraphic model.
  • the described methodology is less expensive from a computational standpoint and can improve the operation of a computer, at least in reducing computational requirements to generate a forward stratigraphic model.
  • FIG. 1 is a flowchart illustrating an example of a computer-implemented method for inverse stratigraphic modeling using a hybrid linear and nonlinear algorithm, according to an implementation of the present disclosure.
  • FIG. 2 A is a graph illustrating an example initial topography of a forward model, according to an implementation of the present disclosure.
  • FIG. 2B is a graph illustrating an example sinusoidal sea-level curve, according to an implementation of the present disclosure.
  • FIG. 3A is a data plot illustrating development of delta lobe along a centerline, according to an implementation of the present disclosure.
  • FIG. 3B is a graph illustrating travel distance of different grain sizes, according to an implementation of the present disclosure.
  • FIG. 3C is a data plot illustrating a 2D profile of a delta simulation without compaction, according to an implementation of the present disclosure.
  • FIG. 3D is a data plot illustrating a 2D profile of a delta simulation with compaction, according to an implementation of the present disclosure.
  • FIG. 4A is a graph illustrating that a cost function of grain size 1 exhibits non-convex features when multiple (more than two) grain sizes are considered, according to an implementation of the present disclosure.
  • FIG. 4B is a graph illustrating that a cost function of grain size 1 exhibits convex features when only a sand ratio is considered, according to an implementation of the present disclosure.
  • FIG. 5 are graphs illustrating a difference between observed data and an inverse result, according to an implementation of the present disclosure.
  • FIG. 6A is a graph illustrating a Genetic Algorithm (GA) implemented for grain size optimization on a high-frequency oscillation cost function, according to an implementation of the present disclosure.
  • GA Genetic Algorithm
  • FIG. 6B is a graph illustrating selection of global optima by using residual error and correlation, according to an implementation of the present disclosure.
  • FIG. 7 are graphs illustrating a comparison between simulated well data and real (observed) well data, according to an implementation of the present disclosure.
  • FIG. 8 A is a data plot of a location of an F3 block in the North Sea, according to an implementation of the present disclosure.
  • FIG. 8B is a Late Miocene to Pleistocene data plot of four system tracts, according to an implementation of the present disclosure.
  • FIG. 8C is a data plot of a location of an inline and a well, according to an implementation of the present disclosure.
  • FIG. 9A is a graph illustrating a difference between observed inline profile data and an inverse (simulation) result, according to an implementation of the present disclosure.
  • FIG. 9B is a data plot of an observed inline seismic profile, according to an implementation of the present disclosure.
  • FIG. 9C is a data plot illustrating an overlap of observed inline profiled data and an inverse (simulation) result, according to an implementation of the present disclosure.
  • FIG. 10A is a graph of a comparison of simulated well data and well porosity log, according to an implementation of the present disclosure.
  • FIG. 10B is a data plot of a well penetrating a deltaic system, according to an implementation of the present disclosure.
  • FIG. 11 is a block diagram illustrating an example of a computer- implemented system used to provide computational functionalities associated with described algorithms, methods, functions, processes, flows, and procedures, according to an implementation of the present disclosure.
  • a forward model for example, a two-dimensional (2D) hydrodynamic forward model designed for a fluvial-dominated delta system forward simulation
  • 2D two-dimensional hydrodynamic forward model designed for a fluvial-dominated delta system forward simulation
  • a forward model is a digital representation of detailed internal geometry and rock properties of a subsurface earth volume, such as a petroleum reservoir or a sediment-filled basin.
  • Forward models are used as forward engines for stratigraphic forward simulations. For example, calculations can be performed on a particular forward model to investigate/simulate formation of stratigraphic or sedimentary geologic layers on geological timescales and to provide geologic reservoir simulations for selection of locations for new hydrocarbon wells, estimation of hydrocarbon reserves, and planning hydrocarbon-reservoir- development strategies.
  • Stratigraphic forward models can be divided into four main sub- categories: 1) geometrical; 2) fuzzy logic; 3) diffusion; and 4) hydrodynamic process- response.
  • the geometrical model describes geometric results of a depositional process from geological aspect rather than the physical process itself.
  • the fuzzy logic model controls spatial and temporal distribution of sediments through a series of logical rules. Diffusion models assume that transportation and deposition of sediment could be described by diffusion processes with certain diffusion gradients on a pre-existing sediment mass/surface.
  • the hydrodynamic process-response model follows fundamental hydrodynamic laws and can numerically model different types of depositional system processes. [0035] Three prerequisites are typically needed to consider stratigraphic forward simulation to be successful.
  • an example of a standard input file that can be used with the described methodology can include the following input parameters.
  • the input file includes particular input parameters for the described inversion algorithm (that is, the Particle Swarm Optimization (PSO) and Genetic algorithm (GA) sections):
  • This input parameters list contains two parts, including forward model # parameters and inverse model parameters.
  • the leveraged inverse algorithm can only provide generally-optimized parameters based on the overall result of the forward simulation.
  • the generally-optimized parameters are normally a poor fit for each time step associated with the forward simulation.
  • the forward simulation is expensive with respect to time, and the information provided by the forward simulation is not be fully usable.
  • Variations of multiple input parameter values require consideration of many possible combinations. Since the forward model can be exercised multiple times in each iteration based on the possible input parameter value combinations, it is often necessary to artificially constrain the number of iterations, which limits inversion accuracy.
  • the current disclosure describes calibrating a forward model with multiscale prior observation data (for example, well data and seismic data), which is more accurate and consistent than conventional methods of generating input parameter values.
  • the generated forward model can be used to improve overall accuracy of geological predictions.
  • the described methodology improves upon building an objective function for calculation of a misfit and optimizing input parameters using the calculated misfit and an inverse algorithm.
  • a deltaic-type system is used, as formation processes for delta deposits have been well-studied and mathematically described by diffusion models, simplified models, and hydrodynamic process-response models with different fluid motion equations (for example, momentum, advection-diffusion, and stochastic parcel-based cellular routing).
  • an advection-diffusion equation is utilized as a forward engine.
  • PSO an Augmented Lagrange Multiplier (ALM) and a GA are also employed to build a two- step inverse algorithm.
  • An interpreted seismic surface, lithology thickness, and porosity extracted from well-logs are utilized for objective function establishment.
  • the described two-step inverse algorithm is used for input parameter optimization of a forward model designed for clastic deposits to quickly estimate some of the fundamental input parameters (for example, discharge, velocity, and components of different grain sizes) needed to conduct a detailed and full physical forward simulation.
  • Clastic deposits are composed of clasts (fragments) of pre-existing minerals and rock.
  • a clast can be considered to a fragment of geological detritus, such as chunks and smaller grains of rock broken off other rocks by physical weathering.
  • Clastic is used with reference to sedimentary rocks and particles in sediment transport, whether in a suspension or as bed load, and in sediment deposits.
  • the forward model is a 2D forward model designed for a fluvial-dominated delta system simulation.
  • the kernel engine of the forward model can be a steady-state 2D advection- diffusion equation in-plane and a linear model in the vertical direction, while compaction is considered.
  • the two-step inverse algorithm is designed based on forward model characteristics.
  • topography (geomorphological information) is used to build an objective function for a sand ratio and other parameter optimization (for example, hydrodynamic parameters) in each time-step of a forward simulation.
  • PSO and ALM are the engines of the piecewise inversion in the first step.
  • a global optimization algorithm for example, PSO
  • a linear optimization algorithm for example, ALM
  • ALM is then executed to solve components of each grain size and calculate misfits of topography produced from simulation and real data.
  • PSO or another global optimization algorithm is applied to update particles (that is, input parameters) and to search for an optimal solution according to the calculated misfits.
  • an optimal solution is found, different grain sizes are classified as sand and mud, and their components are added, respectively, to generate a sand ratio. Components of different grain sizes are calculated through linear optimization and other parameters are estimated through global optimization.
  • the second step of the inverse algorithm which generates a final optimization
  • thickness of different grain sizes and porosity (for example, extracted from well logs and well cores) is used to build an objective function for calibrating components of different grain sizes under a constraint of the sand ratio estimated from the first step.
  • a GA or another global algorithm is used as a kernel engine.
  • Correlation and residual error or other criterion are used to rank a similarity between simulation and observation.
  • model had been partly rewritten as a convex function in the first step, permitting convex function optimization in the first step.
  • a non-linear algorithm was utilized to optimize hydrodynamic parameters involved in a 2D advection-diffusion in-plane.
  • a quadratic optimization of a linear model was performed to determine components of different grain sizes, which control a vertical stacking pattern. The non-linear algorithm significantly improves the inversion efficiency and can be used on other forward models driven by hybrid linear and nonlinear algorithms.
  • different grain size components calculated were reduced dimensionally by calculating a sand ratio in the first step of the inversion.
  • the one-dimensional sand ratio is used to approximate an influence of high dimensional grain size components on topography, which is the low-frequency trend of forward model.
  • This dimension reduction process can be considered a smoothing of a cost function under a prior geological constraint.
  • the objective function is converted from a non-convex function to convex function for an accurate and stable solution.
  • Different grain size components are optimized by using well data in the second step of inversion under the constraint of the sand ratio calculated in the first step of the inversion.
  • the hierarchical/cascading inversion strategy reduces dimensions of parameters in each step, reduces the number of iterations needed, processing requirements, and increases solution stability. 3.
  • the application of a piecewise inversion results in the most fit parameters for each simulation step when compared to conventional approaches. Cumulative error is minimized and accuracy improved. Computation cost is also reduced since an optimization can be performed in each step rather than waiting until an entire forward model has been processed.
  • FIG. 1 is a flowchart illustrating an example of a computer-implemented method 100 for inverse stratigraphic modeling using a hybrid linear and nonlinear algorithm, according to an implementation of the present disclosure.
  • method 100 can be performed, for example, by any system, environment, software, and hardware, or a combination of systems, environments, software, and hardware, as appropriate.
  • various steps of method 100 can be run in parallel, in combination, in loops, or in any order.
  • method 100 contains a forward model and a two-step inverse engine.
  • the two-step inverse engine leverages an inverse algorithm to output a set of well-optimized parameters for subsequent forward modeling.
  • Seismic data is utilized in the first step (102) of the two-step inverse algorithm to catch a low-frequency trends (that is, topography of a delta deposit). Hydrodynamic input parameters (for example, discharge and velocity), and a sand ratio are optimized.
  • Well data obtained from penetrating the delta deposit are employed in the second step (104) of the two-step inverse algorithm to calibrate high-frequency input parameters (for example, porosity and reservoir thickness) of the delta deposit. Components of different grain sizes are optimized.
  • the resultant, well optimized input parameters can be used with a forward simulation generate an additional result.
  • a 2D hydrodynamic forward model is used simulate progradation and retrogradation of a fluvial-dominated delta system.
  • a relationship of liquid discharge, large-scale relief, and basin temperature has been used to calculate sediment concentration. The relationship reduces an amount of input parameters required for non-linear optimization, since sediment concentration is not considered an independent input parameter.
  • the kernel engine of the forward model is a steady- state 2D advection-diffusion equation in-plane, and a linear model in the vertical direction, while compaction is considered.
  • Suspended load is the primary mechanism of delta deposit development.
  • Three-dimensional (3D) suspended sediment motion can be expressed (as in Equation (1)) as sediment mass flux out of a control volume plus a mass variation rate in the control volume:
  • ps is the mass density of sediment
  • C is the instantaneous concentration in volume
  • u is instantaneous velocity in x direction
  • v is instantaneous velocity in y direction
  • w is instantaneous velocity in z direction.
  • Equation (2) a steady-state 2D advection- diffusion equation (as expressed in Equation (2)) is used to describe different grain size transportation in a lateral direction, without considering the suspended load variation in vertical direction:
  • q Sz is the suspended-load amount of different grain size in volume for each time step
  • pi is the components of different grain size
  • N is the number of grain size
  • C is instantaneous concentration in volume
  • u is instantaneous velocity in volume.
  • First Step (102) In the first step of the inverse algorithm, topography (geomorphological information) is used to build an objective function for a sand ratio and other parameter optimization (for example, hydrodynamic input parameters) in each time-step of a forward simulation using a forward model.
  • topography geomorphological information
  • other parameter optimization for example, hydrodynamic input parameters
  • a scope of input parameters is defined.
  • the scope can be defined by an input parameter (for example, as illustrated in the previously-defined input file).
  • the scope of the parameter“Discharge” can be defined by:
  • the parameter range will range from 100 to 500.
  • the discharge range is constrained by the size of a river channel and statistical result of known global modem river discharge amounts.
  • scope values can be entered by a user or dynamically determined by an automated software process. From 106, method 100 proceeds to 108.
  • the defined input parameter scopes are supplied to a PSO algorithm.
  • the PSO algorithm is configured to select a value randomly from the defined scope values for each input parameter with a possible range. For example, the PSO algorithm selects 300 in a first iteration for the Discharge input parameter. For each time step, a group of random particles are initialized to calculate a distribution of different grain sizes by solving a 2D advection-diffusion equation or other nonlinear equation. From 108, method 100 proceeds to 110.
  • the selected parameter values are combined with other fixed parameters (for example, from an input file) and supplied into a forward model to generate a simulated topographical result. From 110, method 100 proceeds to 112.
  • a detailed seismic data interpretation is generated to represent specific seismic features/observed topography (for example, Fig 8B is an illustration of a clinoform feature).
  • the detailed seismic data interpretation is different from a conventional seismic data interpretation, in that the detailed seismic data interpretation requires a depiction of not only atop surface, but also all time-equivalent seismic events within a target depositional unit constrained by interpreted seismic (strata) surfaces representing a geological surface or topography of each geological time step within the target depositional unit.
  • the interpreted seismic surfaces are preserved (for example, in an ASCII format). From 112, method 100 proceeds to 114.
  • a misfit value (Dl) is calculated.
  • Dl misfit value
  • L2-Norm of the simulated topography and interpreted seismic surfaces is utilized as an objective function for misfit calculation.
  • the misfit can be calculated by comparing the generated interpreted seismic surface of 112 with the simulated topographical result from 110:
  • Dl represents the Euler distance between simulated topography (Topography simulated, for example, at 1 10) and observed topography (Topography seismic, for example, at 1 12). From 1 14, method 100 proceeds to 1 16.
  • a linear optimization algorithm for example, an augmented
  • Lagrange multiplier (ALM) is applied to minimize the calculated misfit value of 114.
  • the resultant minimization value is used as an optimization in percentage (%) of grain size component values.
  • Optimized grain size components provided by the ALM are classified (for example, as sand and mud), and their components are added respectively to generate a sand ratio. From 116, method 100 proceeds to 118.
  • 108-116 is iterated one or more times to update input parameters until an optimal solution is found according to the ALM-calculated misfits. After all particles are calculated, the optimal solution has been found. With the optimal solution, different grain sizes are classified (for example, as sand and mud) and their components are added, respectively, to generate an estimated optimized sand ratio. From 118, method 100 proceeds to 120.
  • a data set of the generated estimated optimized sand ratio and hydrodynamic input parameters (for example, Discharge and velocity) is produced.
  • the estimated optimized sand ratio and hydrodynamic parameters will be used as fixed input parameters in the Second Step (104). From 120, method 100 proceeds to 122.
  • Second Step (104) In the second step of the inverse algorithm, components of different grain sizes are optimized.
  • the estimated optimized sand ratio and hydrodynamic parameters of 120 are fixed input parameters of the second step. Components of different grain sizes are required for optimization under the constraint of the sand ratio estimated from the first step (102).
  • GA is introduced as the engine of inversion. Correlation and residual error are used to rank similarity between simulation and observation on a 2D graph.
  • components of different grain sizes formed a high dimensional grid with specified intervals and are supplied to the GA as a first generation. In some implementations, only the top 40% seeds with high correlation and low residual error are preserved as parents for a subsequent generation. After four generations, the seed with a highest correlation and lowest residual error is recognized as the optimal seed.
  • the distribution scope of components of different grain size is constrained by the fixed input optimized sand ratio.
  • sand ratio there are four types of grain sizes, including coarse sand, fine sand, siltstone, and mudstone. Coarse sand and fine sand are classified as sand. Siltstone and mudstone are classified as mudstone.
  • the sand ratio is known to be 0.5, meaning 50% sandstone and 50% mudstone. From 122, method 100 proceeds to 124.
  • a GA is used to determine a proportion of each grain size.
  • the range of each grain size is from 0 to 50%.
  • the GA is performed to generate components of different grain sizes within the scope. From 124, method 100 proceeds to 126.
  • Thickness weiis represents lithology thickness extracted from wells
  • Thickness simulated represents lithology thickness extracted from simulation
  • Porosity weiis represents porosity extracted from wells
  • Porosity simulated represents porosity extracted from simulation. From 130, method 100 proceeds to 132. [0074] At 132, the misfit function (130) is used with the GA to iterate (124-130) one or more times until optimized components of different grain sizes are obtained. From 132, method 100 proceeds to 134.
  • the described methodology was tested with the described 2D hydrodynamic forward model on a synthetic, geologically- reasonable forward simulation result as observed data.
  • six input variables including discharge and components of five grain sizes, were calibrated.
  • the most optimal determined solution had a correlation of 97% and a residual error of 1.86, suggesting that the estimated input values were very close to true values.
  • FIG. 2A is a graph illustrating an example initial topography 200a of a forward model, according to an implementation of the present disclosure.
  • FIG. 2A has a horizontal axis 202a representing a Distance n (meters (m)) times 10 m and a vertical axis 204a representing topographical height (m). The topographical value is represent by 206a.
  • a forward simulation was conducted to provide test data at the beginning. The forward simulation was performed on a ramp 12.5 kilometers (km) long and 150 (m) high (as can be seen in FIG. 2A).
  • the knickpoint 208a was located about 2 km away from the left border.
  • the slope 2l0a is about 2 degrees above the knickpoint 208a.
  • the slope 2l2a is about 8 degrees below the knickpoint 208a.
  • FIG. 2B is a graph illustrating an example sinusoidal sea-level curve
  • FIG. 2B has a horizontal axis 202b representing Age in kilo annum (ka) and a vertical axis 204b representing Amplitude (m).
  • a sinusoidal fluctuation curve 206b ranges from 136 m to 149 m with a 100 year time-step and 2.0 ka duration used as sea level.
  • a stochastic discharge model consisting of a steady average discharge composition (500 m3/s) and a stochastic composition generated by a Markov process based climate-dependent discharge time series was used to produce variable discharge parameter for each step. Concentration was calculated by a relationship between catchment properties and river discharge series. Five classes of sediments, including coarse sand (1.2 millimeters (mm)), medium sand (0.3 mm), fine sand (0.15 mm), silt sand (0.06 mm), and clay (0.003 mm) with fixed components (35%, 10%, 15%, 20%, 20%, respectively) in each step were incorporated in the forward model.
  • FIG. 3 A is a data plot 300a illustrating development of delta lobe along a centerline, according to an implementation of the present disclosure.
  • FIG. 3A has a horizontal axis 302a representing Width (m) in a y direction times 30 m and a vertical axis 304a representing Distance (m) in an x direction times 10 m.
  • a total of 20 time steps were used to emulate the development of a fluvial-dominated delta system.
  • a delta lobe prograded 306a along a centerline 308a when sea level fell and retrograded 3l0a basinward when sea level rose.
  • FIG. 3B is a graph 300b illustrating travel distance of different grain sizes, according to an implementation of the present disclosure.
  • FIG. 3B has a horizontal axis 302b representing Distance (m) in an x direction times 10 m and a vertical axis 304b representing a Total Thickness (m). Travel distances of five grain sizes with different color (306b) are depicted. Grain sizes from 1 to 5, represent coarse sand to mud, respectively. The coarse grain size exhibits an s-shape 308b of depositional curve with short travel distance, while the mud exhibits a flat-shaped curve 3 lOb. This feature may be attributed to the rapid loss of momentum of grains the more course a grain size is and demonstrate that fine sized deposits travel much longer distances compared to coarse sized deposits.
  • FIG. 3C is a data plot 300c illustrating a 2D profile of a delta simulation without compaction, according to an implementation of the present disclosure.
  • FIG. 3C has a horizontal axis 302c representing Distance (m) in an x direction times 10 m and a vertical axis 304c representing Depth in decimeters (dm).
  • the illustrated 2D profile shows clinoform 306c (a sloping depositional surface of a major morphological feature giving seismic expression) evolution without compaction. Colors from black to grey (as in Grain size bar 308c) represent grain sizes from fine to coarse, respectively.
  • FIG. 3D is a data plot 300d illustrating a 2D profile of a delta simulation with compaction, according to an implementation of the present disclosure.
  • FIG. 3D has a horizontal axis 302d representing Distance (m) in an x direction times 10 m and a vertical axis 304d representing Depth (dm).
  • the illustrated 2D profile shows clinoform evolution 306d without compaction. Colors from black to grey (as in Grain size bar 308d) represent grain sizes from fine to coarse, respectively.
  • FIG. 4A is a graph 400a illustrating that a cost function of grain size 1 exhibits non-convex features when multiple (more than two) grain sizes are considered, according to an implementation of the present disclosure.
  • FIG. 4A has a horizontal axis 402a representing a Number of grain size sets and a vertical axis 404a representing Errors (x 10 4 ).
  • Line 406a represents Error of different grain size sets, while line 408a represents smoothed errors.
  • Components of different grain sizes were obtained through described linear optimization.
  • the cost function of each grain size is actually a non-convex function when multiple (more than two) grain sizes are considered.
  • the linear optimization method is inclined to generate a local minimum other than to find a global optimization.
  • FIG. 4B is a graph 400b illustrating that a cost function of grain size 1 exhibits convex features when only a sand ratio is considered, according to an implementation of the present disclosure.
  • FIG. 4B has a horizontal axis 402b representing a Number of sand ratio sets and a vertical axis 404b representing Errors (x 10 5 ).
  • Line 406b represents Error of sad ratio sets, while line 408b represents smoothed errors.
  • five grain sizes were classified into sand and mud according to geological classifications. Their components were added respectively to calculate a sand ratio and a mud ratio.
  • the described methodology squeezed the four dimensional cost function (that is, five components with four degrees- of-freedom) into one dimension under the strong constraint of geological concepts and generated a smooth cost function of a sand ratio.
  • the prediction of the sand ratio became stable and reliable.
  • FIG. 5 are graphs illustrating a difference 500 between observed data and an inverse result, according to an implementation of the present disclosure.
  • FIG. 5 has a horizontal axis 502 representing Distance (m) in an x direction and a vertical axis 504 representing Topographic height (m).
  • An original interpreted topography surface 506 is inset in the bottom left of FIG. 5.
  • a comparison between an inverse result 508 and observed data 510 against actual Topography 512 as in graph 514 shows that the simulation inverse result 508 almost overlaps with the interpreted topography surface 506, suggesting discharge and sand ratio were successfully optimized through the combination of PSO and ALM with more than 95% accuracy. The low frequency trend of the observation was well captured.
  • FIG. 6A is a graph 600a illustrating a GA implemented for grain size optimization on a high-frequency oscillation cost function, according to an implementation of the present disclosure.
  • FIG. 6A has a horizontal axis 602a representing a Number of datasets and a vertical axis 604a representing Error estimation as a scale representing the Euclidean distance of the simulation result and well data.
  • FIG. 6B is a graph 600b illustrating selection of global optima by using residual error and correlation, according to an implementation of the present disclosure.
  • FIG. 6B has a horizontal axis 602b representing Error between simulated well and real well data and a vertical axis 604b representing Correlation between simulated well and real well data.
  • the vertical axis represents Error estimation as a scale representing the Euclidean distance of the simulation result and well data.
  • the optimal solution 606b with a highest correlation (97%) and lowest residual error (1.86) was found.
  • FIG. 7 are graphs 700 illustrating a comparison between simulated well data and real (observed) well data, according to an implementation of the present disclosure.
  • FIG. 7 includes three wells (well 1 701 a, well 2 70 lb, and well 3 701 c. Each well is represented by two separate graphs of simulated well data and real well data (that is, Well 1 70la - 702a and 702b, respectively; Well 2 70lb - 704a and 704b, respectively; and Well 3 70lc - 706a and 706b, respectively).
  • Each graph (for example 702a) has a horizontal axis representing Lithology (for example, a category) and a vertical axis representing Thickness (m).
  • each well has a horizontal axis 708 representing Lithology and a vertical axis 710 representing Thickness (m).
  • FIG. 8A is a data plot 800a of a location of an F3 block in the North Sea, according to an implementation of the present disclosure.
  • the Netherlands Offshore F3 seismic survey was collected for inverse stratigraphic modeling verification. The data was derived from OPENDTECT free and public dataset (generated by dGB Earth Sciences, The Netherlands).
  • F3 Block North Sea 802a is a block covered by a 3D seismic survey in the Dutch sector of the North Sea.
  • FIG. 8B is a Late Miocene to Pleistocene data plot 800b of four system tracts, according to an implementation of the present disclosure.
  • the horizontal axis 80lb represents Distance (km), while the vertical axis 802b represents Time in milliseconds (ms).
  • TST transgressive
  • HST highstand
  • FSST falling stage
  • LST lowstand
  • HST 804b is dominated by a fluviodeltaic system (delta deposit 810b) with large scale sigmoidal clinoform.
  • the delta deposit 810b consists of sand and shale, with an overall high porosity. Bright spots observed on the profile may be attributed to biogenic gas pockets.
  • FSST 806b is dominated by mass transport deposits (MTD) 8l2b.“Target for Simulation” area 814 is dominated by a clinoform feature, which is an indicator of deltaic system.
  • Horizontal scale 816b and vertical scale 818b provide scale for the horizontal and vertical axes, respectively.
  • FIG. 8C is a data plot 800c of a location of an inline and a well, according to an implementation of the present disclosure.
  • the inline is identified as Inline 425 802c and the well as Well F03-4 804c.
  • a detailed interpretation of Inline 425 802c located in the center of the F3 survey was employed to provide seismic geomorphological information. The profile had been converted to a depth domain through seismic velocity field data. Since no grain size data was available for this verification, the delta deposit 810b was simply divided into sand and mud with assumed reasonable grain sizes (that is, 0.2 mm and 0.002 mm, respectively).
  • FIG. 9A is a graph 900a illustrating a difference between observed inline profile data and an inverse (simulation) result, according to an implementation of the present disclosure.
  • FIG. 9A has a horizontal axis 902a representing Distance (km) in an x direction and a vertical axis 904a representing a Topographic height (m).
  • a comparison of the inverse result 906a and observed seismic geomorphology data 908a as in overlapping data 910a exhibits a high similarity (approximately 90%).
  • the residual error (approximately 10%) can be attributed to the toe of the delta front. This misfit may be attributed to the limitation of seismic resolution. Note that some seismic events merge together at the base (for example, at 912a) since the deposits becomes too thin in this area.
  • FIG. 9B is a data plot 900b of an observed inline seismic profile, according to an implementation of the present disclosure.
  • the horizontal axis 902b represents Distance (km) and the vertical axis 904b represents Time (ms).
  • Horizontal scale 906b and vertical scale 908b provide scale for the horizontal and vertical axes, respectively.
  • the inline 910b is the Inline 425 (802c) of FIG. 8C.
  • the target area for the simulation is indicated approximately at the center of dashed circle 9l2b.
  • FIG. 9C is a data plot 900c illustrating an overlap of observed inline profiled data and an inverse (simulation) result, according to an implementation of the present disclosure.
  • the horizontal axis 902c represents Distance (km) and the vertical axis 904c represents Time (ms).
  • Horizontal scale 906c and vertical scale 908c provide scale for the horizontal and vertical axes, respectively.
  • Overlap (for example, at 9l0c) of the well-tuned inverse result on the observed seismic profile of FIG. 9B demonstrates that the inverse result is basically in accordance with the observed seismic profile of the deltaic deposit with about 90% accuracy.
  • FIG. 10A is a graph lOOOa of a comparison of simulated well data and well porosity log, according to an implementation of the present disclosure.
  • Graph lOOOa includes simulated well data l002a in a left pane and the well porosity log l004a of well F03-4 (FIG. 8C) in the right pane.
  • the comparison of extracted well data from an inverse simulation result and the porosity log of well F03-4 shows that the simulation could generally reflect the variation trend of real data.
  • well F03-4 was divided into four sections l006a, l008a, lOlOa, and l0l2a.
  • the top two sections (l006a and l008a) mainly consist of one reverse cycle with about a 0.5 to 0.6 sand ratio.
  • the third section (lOlOa) contains two sets of sand bodies and the lower sand bodies have larger thickness.
  • the bottom section (l0l2a) is characterized by several fine sand layers (laminations) in mud. These features are well reproduced by simulation result, suggesting that the well calibrated forward model do have a good prediction for the real case.
  • FIG. 10B is a data plot lOOOb of a well penetrating a deltaic system, according to an implementation of the present disclosure.
  • the horizontal axis l002b represents Distance (km) and the vertical axis l004b represents Time (ms).
  • the well is Well F03-4 (804c of FIG. 8C) and the deltaic system is delta deposit 8l0b (as in FIG. 8B).
  • the inverse result indicates a reasonable and quantitative geological interpretation of the sequence stratigraphic development and formation of the delta deposit 810b.
  • Horizontal scale l006b and vertical scale l008b provide scale for the horizontal and vertical axes, respectively.
  • the described methodology can be configured to send messages, instructions, or other communications to a computer-implemented controller, database, or other computer-implemented system to dynamically initiate control of, control, or cause another computer-implemented system to perform a computer-implemented or other function/operation.
  • operations based on data, operations, outputs, or interaction with a graphical user interface can be transmitted to cause operations associated with a computer, database, network, or other computer-based system to perform storage efficiency, data retrieval, or other operations consistent with this disclosure.
  • GUI graphical user interface
  • interacting with any illustrated GUI can automatically result in one or more instructions transmitted from the GUI to trigger requests for data, storage of data, analysis of data, or other operations consistent with this disclosure.
  • transmitted instructions can result in control, operation, modification, enhancement, or other operations with respect to a tangible, real-world piece of computing or other equipment.
  • the described GUIs can send a request to slow or speed up a computer database magnetic/optical disk drive, shut down/activate a computing system, cause a network interface device to disable, throttle, or increase data bandwidth allowed across a network connection, or sound an audible/visual alarm (for example, a mechanical alarm/light emitting device) as a notification of detected malicious behavior(s) with respect to a computing system(s) used with the described methodology or interacting with the computing system(s) used with the described methodology.
  • an audible/visual alarm for example, a mechanical alarm/light emitting device
  • the output of the described approach can be used to dynamically influence, direct, control, influence, or manage tangible equipment related to hydrocarbon production, analysis, and recovery.
  • real-time data received from an ongoing drilling operation can be analyzed using the described methodology.
  • a wellbore trajectory can be modified, a drill speed can be increased or reduced, a drill can be stopped, an alarm can be activated/deactivated (such as, visual, auditory, or voice alarms), refinery or pumping operations can be affected (for example, stopped, restarted, accelerated, or reduced).
  • Other examples can include alerting geosteering and directional drilling staff when incorrect directional survey data has been detected (such as, with a visual, auditory, or voice alarm).
  • the described approach can be integrated as part of a dynamic control system for any hydrocarbon- related equipment consistent with this disclosure.
  • FIG. 11 is a block diagram illustrating an example of a computer- implemented System 1100 used to provide computational functionalities associated with described algorithms, methods, functions, processes, flows, and procedures, according to an implementation of the present disclosure.
  • System 1100 includes a Computer 1102 and a Network 1130.
  • the illustrated Computer 1102 is intended to encompass any computing device such as a server, desktop computer, laptop/notebook computer, wireless data port, smart phone, personal data assistant (PDA), tablet computer, one or more processors within these devices, another computing device, or a combination of computing devices, including physical or virtual instances of the computing device, or a combination of physical or virtual instances of the computing device.
  • a server desktop computer
  • laptop/notebook computer wireless data port
  • smart phone smart phone
  • PDA personal data assistant
  • tablet computer one or more processors within these devices
  • another computing device or a combination of computing devices, including physical or virtual instances of the computing device, or a combination of physical or virtual instances of the computing device.
  • the Computer 1102 can include an input device, such as a keypad, keyboard, touch screen, another input device, or a combination of input devices that can accept user information, and an output device that conveys information associated with the operation of the Computer 1102, including digital data, visual, audio, another type of information, or a combination of types of information, on a graphical-type user interface (UI) (or GUI) or other UI.
  • UI graphical-type user interface
  • the Computer 1102 can serve in a role in a distributed computing system as a client, network component, a server, a database or another persistency, another role, or a combination of roles for performing the subject matter described in the present disclosure.
  • the illustrated Computer 1102 is communicably coupled with a Network 1130.
  • one or more components of the Computer 1102 can be configured to operate within an environment, including cloud-computing-based, local, global, another environment, or a combination of environments.
  • the Computer 1102 is an electronic computing device operable to receive, transmit, process, store, or manage data and information associated with the described subject matter. According to some implementations, the Computer 1102 can also include or be communicably coupled with a server, including an application server, e-mail server, web server, caching server, streaming data server, another server, or a combination of servers.
  • a server including an application server, e-mail server, web server, caching server, streaming data server, another server, or a combination of servers.
  • the Computer 1102 can receive requests over Network 1130 (for example, from a client software application executing on another Computer 1102) and respond to the received requests by processing the received requests using a software application or a combination of software applications.
  • requests can also be sent to the Computer 1102 from internal users (for example, from a command console or by another internal access method), external or third-parties, or other entities, individuals, systems, or computers.
  • Each of the components of the Computer 1102 can communicate using a System Bus 1103.
  • any or all of the components of the Computer 1102, including hardware, software, or a combination of hardware and software can interface over the System Bus 1103 using an application programming interface (API) 1112, a Service Layer 1113, or a combination of the API 1112 and Service Layer 1113.
  • the API 1112 can include specifications for routines, data structures, and object classes.
  • the API 1112 can be either computer-language independent or dependent and refer to a complete interface, a single function, or even a set of APIs.
  • the Service Layer 1113 provides software services to the Computer 1102 or other components (whether illustrated or not) that are communicably coupled to the Computer 1102.
  • the functionality of the Computer 1102 can be accessible for all service consumers using the Service Layer 1113.
  • Software services such as those provided by the Service Layer 1113, provide reusable, defined functionalities through a defined interface.
  • the interface can be software written in JAVA, C++, another computing language, or a combination of computing languages providing data in extensible markup language (XML) format, another format, or a combination of formats.
  • XML extensible markup language
  • alternative implementations can illustrate the API 1112 or the Service Layer 1113 as stand-alone components in relation to other components of the Computer 1102 or other components (whether illustrated or not) that are communicably coupled to the Computer 1102.
  • any or all parts of the API 1112 or the Service Layer 1113 can be implemented as a child or a sub-module of another software module, enterprise application, or hardware module without departing from the scope of the present disclosure.
  • the Computer 1102 includes an Interface 1104. Although illustrated as a single Interface 1104, two or more Interfaces 1104 can be used according to particular needs, desires, or particular implementations of the Computer 1102.
  • the Interface 1104 is used by the Computer 1102 for communicating with another computing system (whether illustrated or not) that is communicatively linked to the Network 1130 in a distributed environment.
  • the Interface 1104 is operable to communicate with the Network 1130 and includes logic encoded in software, hardware, or a combination of software and hardware.
  • the Interface 1104 can include software supporting one or more communication protocols associated with communications such that the Network 1130 or hardware of Interface 1104 is operable to communicate physical signals within and outside of the illustrated Computer 1102.
  • the Computer 1102 includes a Processor 1105. Although illustrated as a single Processor 1105, two or more Processors 1105 can be used according to particular needs, desires, or particular implementations of the Computer 1102. Generally, the Processor 1105 executes instructions and manipulates data to perform the operations of the Computer 1102 and any algorithms, methods, functions, processes, flows, and procedures as described in the present disclosure.
  • the Computer 1102 also includes a Database 1106 that can hold data for the Computer 1102, another component communicatively linked to the Network 1130 (whether illustrated or not), or a combination of the Computer 1102 and another component.
  • Database 1106 can be an in-memory, conventional, or another type of database storing data consistent with the present disclosure.
  • Database 1106 can be a combination of two or more different database types (for example, a hybrid in-memory and conventional database) according to particular needs, desires, or particular implementations of the Computer 1102 and the described functionality.
  • two or more databases of similar or differing types can be used according to particular needs, desires, or particular implementations of the Computer 1102 and the described functionality.
  • Database 1106 is illustrated as an integral component of the Computer 1102, in alternative implementations, Database 1106 can be external to the Computer 1102.
  • the Database 1106 holds a forward model 1116 and algorithms 1118 (for example, PSO, ALM, and GA).
  • the Computer 1102 also includes a Memory 1107 that can hold data for the Computer 1102, another component or components communicatively linked to the Network 1130 (whether illustrated or not), or a combination of the Computer 1102 and another component.
  • Memory 1107 can store any data consistent with the present disclosure.
  • Memory 1107 can be a combination of two or more different types of memory (for example, a combination of semiconductor and magnetic storage) according to particular needs, desires, or particular implementations of the Computer 1102 and the described functionality.
  • two or more Memories 1107 or similar or differing types can be used according to particular needs, desires, or particular implementations of the Computer 1102 and the described functionality.
  • Memory 1107 is illustrated as an integral component of the Computer 1102, in alternative implementations, Memory 1107 can be external to the Computer 1102.
  • the Application 1108 is an algorithmic software engine providing functionality according to particular needs, desires, or particular implementations of the Computer 1102, particularly with respect to functionality described in the present disclosure.
  • Application 1108 can serve as one or more components, modules, or applications.
  • the Application 1108 can be implemented as multiple Applications 1108 on the Computer 1102.
  • the Application 1108 can be external to the Computer 1102.
  • the Computer 1102 can also include a Power Supply 1114.
  • the Power Supply 1114 can include a rechargeable or non-rechargeable battery that can be configured to be either user- or non-user-replaceable.
  • the Power Supply 1114 can include power-conversion or management circuits (including recharging, standby, or another power management functionality).
  • the Power Supply 1114 can include a power plug to allow the Computer 1102 to be plugged into a wall socket or another power source to, for example, power the Computer 1102 or recharge a rechargeable battery.
  • Computers 1102 there can be any number of Computers 1102 associated with, or external to, a computer system containing Computer 1102, each Computer 1102 communicating over Network 1130.
  • client “user,” or other appropriate terminology can be used interchangeably, as appropriate, without departing from the scope of the present disclosure.
  • present disclosure contemplates that many users can use one Computer 1102, or that one user can use multiple computers 1102.
  • Described implementations of the subject matter can include one or more features, alone or in combination.
  • a computer-implemented method for inverse stratigraphic modeling comprising: a first step, further comprising: selecting a defined scope value for each of a plurality of hydrodynamic input parameters using a global optimization algorithm; generating a simulated topographical result using the selected scope values and a forward model; generating a detailed seismic interpretation to represent specific seismic features or observed topography; minimizing a calculated a misfit value, wherein the misfit value is calculated as a distance between the simulated topographical result and the detailed seismic interpretation; and; and generating an estimated optimized sand ratio and optimized hydrodynamic input parameters; and a second step, further comprising: determining, using a genetic algorithm, a proportion of each grain size in the estimated optimized sand ratio; using a misfit value calculated from thickness and porosity data extracted from well data and a simulation result generated by the forward model to generate optimized components of different grain sizes; and generating optimized hydrodynamic input parameters and optimized components of different grain sizes.
  • a first feature combinable with any of the following features, wherein the global optimization algorithm is configured to randomly select a value from each of the defined scopes.
  • a second feature combinable with any of the previous or following features, wherein the minimization of the calculated misfit value is optimized using a linear optimization algorithm.
  • a third feature combinable with any of the previous or following features, wherein the estimated optimized sand ratio is used as a fixed parameter to optimize components of different grain sizes within the estimated optimized sand ratio.
  • a fourth feature combinable with any of the previous or following features, wherein the distribution scope of components of different grain size is constrained by the estimated optimized sand ratio.
  • a fifth feature combinable with any of the previous or following features, wherein generating the estimated optimized sand ratio and optimized hydrodynamic input parameters and generating optimized hydrodynamic input parameters and optimized components of different grain sizes are performed by iterating, respectively, across operations of the first step and operations of the second step.
  • a sixth feature combinable with any of the previous or following features, further comprising performing forward modeling using the generated optimized hydrodynamic input parameters and optimized components of different grain sizes.
  • a non-transitory, computer-readable medium storing one or more instructions executable by a computer system to perform operations comprising: a first step, further comprising: selecting a defined scope value for each of a plurality of hydrodynamic input parameters using a global optimization algorithm; generating a simulated topographical result using the selected scope values and a forward model; generating a detailed seismic interpretation to represent specific seismic features or observed topography; minimizing a calculated a misfit value, wherein the misfit value is calculated as a distance between the simulated topographical result and the detailed seismic interpretation; and generating an estimated optimized sand ratio and optimized hydrodynamic input parameters; and a second step, further comprising: determining, using a genetic algorithm, a proportion of each grain size in the estimated optimized sand ratio; using a misfit value calculated from thickness and porosity data extracted from well data and a simulation result generated by the forward model to generate optimized components of different grain sizes; and generating optimized hydrodynamic input parameters and optimized components of different grain sizes.
  • a first feature combinable with any of the following features, wherein the global optimization algorithm is configured to randomly select a value from each of the defined scopes.
  • a second feature combinable with any of the previous or following features, wherein the minimization of the calculated misfit value is optimized using a linear optimization algorithm.
  • a third feature combinable with any of the previous or following features, wherein the estimated optimized sand ratio is used as a fixed parameter to optimize components of different grain sizes within the estimated optimized sand ratio.
  • a fourth feature combinable with any of the previous or following features, wherein the distribution scope of components of different grain size is constrained by the estimated optimized sand ratio.
  • a fifth feature combinable with any of the previous or following features, wherein generating the estimated optimized sand ratio and optimized hydrodynamic input parameters and generating optimized hydrodynamic input parameters and optimized components of different grain sizes are performed by iterating, respectively, across operations of the first step and operations of the second step.
  • a sixth feature combinable with any of the previous or following features, further comprising performing forward modeling using the generated optimized hydrodynamic input parameters and optimized components of different grain sizes.
  • a computer-implemented system comprising: one or more computers; and one or more computer memory devices interoperably coupled with the one or more computers and having tangible, non-transitory, machine- readable media storing one or more instructions that, when executed by the one or more computers, perform one or more operations comprising: a first step, further comprising: selecting a defined scope value for each of a plurality of hydrodynamic input parameters using a global optimization algorithm; generating a simulated topographical result using the selected scope values and a forward model; generating a detailed seismic interpretation to represent specific seismic features or observed topography; minimizing a calculated a misfit value, wherein the misfit value is calculated as a distance between the simulated topographical result and the detailed seismic interpretation; and generating an estimated optimized sand ratio and optimized hydrodynamic input parameters; and a second step, further comprising: determining, using a genetic algorithm, a proportion of each grain size in the estimated optimized sand ratio; using a misfit value calculated from thickness
  • a first feature combinable with any of the following features, wherein the global optimization algorithm is configured to randomly select a value from each of the defined scopes.
  • a second feature combinable with any of the previous or following features, wherein the minimization of the calculated misfit value is optimized using a linear optimization algorithm.
  • a third feature combinable with any of the previous or following features, wherein the estimated optimized sand ratio is used as a fixed parameter to optimize components of different grain sizes within the estimated optimized sand ratio.
  • a fourth feature combinable with any of the previous or following features, wherein the distribution scope of components of different grain size is constrained by the estimated optimized sand ratio.
  • a fifth feature combinable with any of the previous or following features, wherein generating the estimated optimized sand ratio and optimized hydrodynamic input parameters and generating optimized hydrodynamic input parameters and optimized components of different grain sizes are performed by iterating, respectively, across operations of the first step and operations of the second step.
  • a sixth feature combinable with any of the previous or following features, further comprising performing forward modeling using the generated optimized hydrodynamic input parameters and optimized components of different grain sizes.
  • Implementations of the subject matter and the functional operations described in this specification can be implemented in digital electronic circuitry, in tangibly embodied computer software or firmware, in computer hardware, including the structures disclosed in this specification and their structural equivalents, or in combinations of one or more of them.
  • Software implementations of the described subject matter can be implemented as one or more computer programs, that is, one or more modules of computer program instructions encoded on a tangible, non-transitory, computer-readable medium for execution by, or to control the operation of, a computer or computer-implemented system.
  • the program instructions can be encoded in/on an artificially generated propagated signal, for example, a machine-generated electrical, optical, or electromagnetic signal that is generated to encode information for transmission to a receiver apparatus for execution by a computer or computer-implemented system.
  • the computer-storage medium can be a machine-readable storage device, a machine-readable storage substrate, a random or serial access memory device, or a combination of computer-storage mediums.
  • Configuring one or more computers means that the one or more computers have installed hardware, firmware, or software (or combinations of hardware, firmware, and software) so that when the software is executed by the one or more computers, particular computing operations are performed.
  • the term“real-time,”“real time,”“realtime,”“real (fast) time (RFT),” “near(ly) real-time (NRT),”“quasi real-time,” or similar terms means that an action and a response are temporally proximate such that an individual perceives the action and the response occurring substantially simultaneously.
  • the time difference for a response to display (or for an initiation of a display) of data following the individual’s action to access the data can be less than 1 ms, less than 1 second (s), or less than 5 s.
  • data processing apparatus refers to data processing hardware.
  • Data processing hardware encompass all kinds of apparatuses, devices, and machines for processing data, including by way of example, a programmable processor, a computer, or multiple processors or computers.
  • the computer can also be, or further include special purpose logic circuitry, for example, a central processing unit (CPU), an FPGA (field programmable gate array), or an ASIC (application-specific integrated circuit).
  • CPU central processing unit
  • FPGA field programmable gate array
  • ASIC application-specific integrated circuit
  • the computer or computer-implemented system or special purpose logic circuitry can be hardware- or software-based (or a combination of both hardware- and software- based).
  • the computer can optionally include code that creates an execution environment for computer programs, for example, code that constitutes processor firmware, a protocol stack, a database management system, an operating system, or a combination of execution environments.
  • code that constitutes processor firmware, a protocol stack, a database management system, an operating system, or a combination of execution environments for example, code that constitutes processor firmware, a protocol stack, a database management system, an operating system, or a combination of execution environments.
  • the present disclosure contemplates the use of a computer or computer-implemented system with an operating system of some type, for example LINUX, UNIX, WINDOWS, MAC OS, ANDROID, IOS, another operating system, or a combination of operating systems.
  • a computer program which can also be referred to or described as a program, software, a software application, a unit, a module, a software module, a script, code, or other component can be written in any form of programming language, including compiled or interpreted languages, or declarative or procedural languages, and it can be deployed in any form, including, for example, as a stand-alone program, module, component, or subroutine, for use in a computing environment.
  • a computer program can, but need not, correspond to a file in a file system.
  • a program can be stored in a portion of a file that holds other programs or data, for example, one or more scripts stored in a markup language document, in a single file dedicated to the program in question, or in multiple coordinated files, for example, files that store one or more modules, sub-programs, or portions of code.
  • a computer program can be deployed to be executed on one computer or on multiple computers that are located at one site or distributed across multiple sites and interconnected by a communication network.
  • While portions of the programs illustrated in the various figures can be illustrated as individual components, such as units or modules, that implement described features and functionality using various objects, methods, or other processes, the programs can instead include a number of sub-units, sub-modules, third-party services, components, libraries, and other components, as appropriate. Conversely, the features and functionality of various components can be combined into single components, as appropriate. Thresholds used to make computational determinations can be statically, dynamically, or both statically and dynamically determined.
  • Described methods, processes, or logic flows represent one or more examples of functionality consistent with the present disclosure and are not intended to limit the disclosure to the described or illustrated implementations, but to be accorded the widest scope consistent with described principles and features.
  • the described methods, processes, or logic flows can be performed by one or more programmable computers executing one or more computer programs to perform functions by operating on input data and generating output data.
  • the methods, processes, or logic flows can also be performed by, and computers can also be implemented as, special purpose logic circuitry, for example, a CPU, an FPGA, or an ASIC.
  • Computers for the execution of a computer program can be based on general or special purpose microprocessors, both, or another type of CPU.
  • a CPU will receive instructions and data from and write to a memory.
  • the essential elements of a computer are a CPU, for performing or executing instructions, and one or more memory devices for storing instructions and data.
  • a computer will also include, or be operatively coupled to, receive data from or transfer data to, or both, one or more mass storage devices for storing data, for example, magnetic, magneto-optical disks, or optical disks.
  • mass storage devices for storing data, for example, magnetic, magneto-optical disks, or optical disks.
  • a computer need not have such devices.
  • a computer can be embedded in another device, for example, a mobile telephone, a personal digital assistant (PDA), a mobile audio or video player, a game console, a global positioning system (GPS) receiver, or a portable memory storage device.
  • PDA personal digital assistant
  • GPS global positioning system
  • Non-transitory computer-readable media for storing computer program instructions and data can include all forms of permanent/non-permanent or volatile/non-volatile memory, media and memory devices, including by way of example semiconductor memory devices, for example, random access memory (RAM), read-only memory (ROM), phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), and flash memory devices; magnetic devices, for example, tape, cartridges, cassettes, intemal/removable disks; magneto-optical disks; and optical memory devices, for example, digital versatile/video disc (DVD), compact disc (CD)-ROM, DVD+/-R, DVD-RAM, DVD-ROM, high-defmition/density (HD)-DVD, and BLU-RAY/BLU- RAY DISC (BD), and other optical memory technologies.
  • semiconductor memory devices for example, random access memory
  • the memory can store various objects or data, including caches, classes, frameworks, applications, modules, backup data, jobs, web pages, web page templates, data structures, database tables, repositories storing dynamic information, or other appropriate information including any parameters, variables, algorithms, instructions, rules, constraints, or references. Additionally, the memory can include other appropriate data, such as logs, policies, security or access data, or reporting files.
  • the processor and the memory can be supplemented by, or incorporated in, special purpose logic circuitry.
  • implementations of the subject matter described in this specification can be implemented on a computer having a display device, for example, a CRT (cathode ray tube), LCD (liquid crystal display), LED (Light Emitting Diode), or plasma monitor, for displaying information to the user and a keyboard and a pointing device, for example, a mouse, trackball, or trackpad by which the user can provide input to the computer.
  • a display device for example, a CRT (cathode ray tube), LCD (liquid crystal display), LED (Light Emitting Diode), or plasma monitor
  • a keyboard and a pointing device for example, a mouse, trackball, or trackpad by which the user can provide input to the computer.
  • Input can also be provided to the computer using a touchscreen, such as a tablet computer surface with pressure sensitivity, a multi-touch screen using capacitive or electric sensing, or another type of touchscreen.
  • Other types of devices can be used to interact with the user.
  • feedback provided to the user can be any form of sensory feedback (such as, visual, auditory, tactile, or a combination of feedback types).
  • Input from the user can be received in any form, including acoustic, speech, or tactile input.
  • a computer can interact with the user by sending documents to and receiving documents from a client computing device that is used by the user (for example, by sending web pages to a web browser on a user’s mobile computing device in response to requests received from the web browser).
  • GUI graphical user interface
  • a GUI can be used in the singular or the plural to describe one or more graphical user interfaces and each of the displays of a particular graphical user interface. Therefore, a GUI can represent any graphical user interface, including but not limited to, a web browser, a touch screen, or a command line interface (CLI) that processes information and efficiently presents the information results to the user.
  • a GUI can include a number of user interface (UI) elements, some or all associated with a web browser, such as interactive fields, pull down lists, and butons. These and other UI elements can be related to or represent the functions of the web browser.
  • UI user interface
  • Implementations of the subject mater described in this specification can be implemented in a computing system that includes a back-end component, for example, as a data server, or that includes a middleware component, for example, an application server, or that includes a front-end component, for example, a client computer having a graphical user interface or a Web browser through which a user can interact with an implementation of the subject mater described in this specification, or any combination of one or more such back-end, middleware, or front-end components.
  • the components of the system can be interconnected by any form or medium of wireline or wireless digital data communication (or a combination of data communication), for example, a communication network.
  • Examples of communication networks include a local area network (LAN), a radio access network (RAN), a metropolitan area network (MAN), a wide area network (WAN), Worldwide Interoperability for Microwave Access (WIMAX), a wireless local area network (WLAN) using, for example, 802.11 a/b/g/n or 802.20 (or a combination of 802.1 lx and 802.20 or other protocols consistent with the present disclosure), all or a portion of the Internet, another communication network, or a combination of communication networks.
  • the communication network can communicate with, for example, Internet Protocol (IP) packets, Frame Relay frames, Asynchronous Transfer Mode (ATM) cells, voice, video, data, or other information between network nodes.
  • IP Internet Protocol
  • ATM Asynchronous Transfer Mode
  • the computing system can include clients and servers.
  • a client and server are generally remote from each other and typically interact through a communication network.
  • the relationship of client and server arises by virtue of computer programs running on the respective computers and having a client-server relationship to each other.
  • any claimed implementation is considered to be applicable to at least a computer-implemented method; a non-transitory, computer-readable medium storing computer-readable instructions to perform the computer-implemented method; and a computer system comprising a computer memory interoperably coupled with a hardware processor configured to perform the computer-implemented method or the instructions stored on the non-transitory, computer-readable medium.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Geophysics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Geophysics And Detection Of Objects (AREA)
PCT/US2019/032489 2018-05-31 2019-05-15 Inverse stratigraphic modeling using a hybrid linear and nonlinear algorithm Ceased WO2019231681A1 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CA3101976A CA3101976A1 (en) 2018-05-31 2019-05-15 Inverse stratigraphic modeling using a hybrid linear and nonlinear algorithm
JP2020566913A JP2021526634A (ja) 2018-05-31 2019-05-15 線形及び非線形のハイブリッドアルゴリズムを用いた逆層序モデル化
EP19728231.2A EP3803470A1 (en) 2018-05-31 2019-05-15 Inverse stratigraphic modeling using a hybrid linear and nonlinear algorithm
CN201980036566.0A CN112689778A (zh) 2018-05-31 2019-05-15 使用混合线性和非线性算法的反演地层建模

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US15/994,586 US11506804B2 (en) 2018-05-31 2018-05-31 Inverse stratigraphic modeling using a hybrid linear and nonlinear algorithm
US15/994,586 2018-05-31

Publications (1)

Publication Number Publication Date
WO2019231681A1 true WO2019231681A1 (en) 2019-12-05

Family

ID=66691015

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2019/032489 Ceased WO2019231681A1 (en) 2018-05-31 2019-05-15 Inverse stratigraphic modeling using a hybrid linear and nonlinear algorithm

Country Status (6)

Country Link
US (1) US11506804B2 (https=)
EP (1) EP3803470A1 (https=)
JP (1) JP2021526634A (https=)
CN (1) CN112689778A (https=)
CA (1) CA3101976A1 (https=)
WO (1) WO2019231681A1 (https=)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP4016137A1 (en) 2020-12-21 2022-06-22 Abu Dhabi National Oil Company Method and system for determining a location of hydrocarbon reservoirs

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11055908B2 (en) * 2018-08-27 2021-07-06 Exxonmobil Research And Engineering Company Inversion of large, nearly-homogeneous geobodies via ultra low-dimensional shape representation using orthogonal basis functions
US11507054B2 (en) * 2018-12-31 2022-11-22 Palo Alto Research Center Incorporated Method and system for hierarchical multi-scale part design with the aid of a digital computer
US20220236448A1 (en) * 2019-05-10 2022-07-28 Totalenergies Se Method for modelling a water current induced by a river in a geological gridded model
US12560740B2 (en) * 2019-05-10 2026-02-24 Totalenergies One Tech Method for modelling the formation of a sedimentary basin using a stratigraphic forward modeling program
CN114945939A (zh) * 2020-01-03 2022-08-26 国际壳牌研究有限公司 润滑剂图像处理和分析
CN115704914B (zh) * 2021-08-17 2026-03-27 中国石油天然气股份有限公司 页岩气ⅰ类储层品质参数确定方法、装置以及存储介质
US12298457B2 (en) * 2022-08-31 2025-05-13 Saudi Arabian Oil Company Method of stratigraphic calibration
GB2639790A (en) * 2023-03-14 2025-10-01 Landmark Graphics Corp Diffusion modeling based subsurface formation evaluation
WO2025175409A1 (en) * 2024-02-19 2025-08-28 Saudi Arabian Oil Company Method and system for rock physics-based grid refinement of grid in forward stratigraphic modeling
CN118228565A (zh) * 2024-04-09 2024-06-21 西安交通大学 一种减少厨房油烟机叶轮内油烟颗粒沉积的方法
CN119006740B (zh) * 2024-07-25 2025-04-15 成都理工大学 一种基于粒子群优化算法的地质模型四边形网格优化方法
CN119622169B (zh) * 2024-08-09 2025-09-30 广州海洋地质调查局 一种基于遗传算法的地形处理方法、装置、设备及介质
CN120296443B (zh) * 2025-04-03 2025-09-02 北京建筑大学 一种适用于跨界水污染的边界探测与定位方法
CN120298617B (zh) * 2025-06-12 2025-08-22 中南大学 一种非线性约束耦合驱动的隐式建模优化方法、系统及存储器
CN120911269B (zh) * 2025-07-23 2026-03-17 延安大学 一种基于大数据的大气污染物快速溯源系统

Family Cites Families (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4783751A (en) 1983-08-17 1988-11-08 University Of South Carolina Analysis of pore complexes
US4868883A (en) 1985-12-30 1989-09-19 Exxon Production Research Company Analysis of thin section images
FR2776393B1 (fr) 1998-03-20 2001-09-07 Inst Francais Du Petrole Methode pour former automatiquement un modele simulant la structure stratigraphique d'une zone souterraine
US6754588B2 (en) 1999-01-29 2004-06-22 Platte River Associates, Inc. Method of predicting three-dimensional stratigraphy using inverse optimization techniques
US6246963B1 (en) 1999-01-29 2001-06-12 Timothy A. Cross Method for predicting stratigraphy
BRPI0511511B1 (pt) * 2004-05-27 2017-01-31 Exxonmobil Upstream Res Co método para prever litologia e porosidade a partir de dados sísmicos de reflexão
AU2005285360B2 (en) * 2004-09-10 2011-02-03 Exxonmobil Upstream Research Company A method for constructing geologic models of subsurface sedimentary volumes
KR100772506B1 (ko) 2005-12-08 2007-11-01 한국전자통신연구원 영상처리를 이용한 지질재료 분류 방법 및 그 장치
TW200834459A (en) 2007-02-05 2008-08-16 Huper Lab Co Ltd Video object segmentation method applied for rainy situations
US20090164186A1 (en) * 2007-12-20 2009-06-25 Bhp Billiton Innovation Pty Ltd. Method for determining improved estimates of properties of a model
WO2009137181A1 (en) * 2008-05-05 2009-11-12 Exxonmobil Upstream Research Company Modeling dynamic systems by visualizing and narrowing a parameter space
MX2012010271A (es) * 2010-03-05 2012-12-17 Vialogy Llc Cálculos de inyección de ruido activo para predictibilidad mejorada en descubrimiento y caracterización de depósitos de petróleo y gas.
US8583410B2 (en) 2010-05-28 2013-11-12 Ingrain, Inc. Method for obtaining consistent and integrated physical properties of porous media
CN102243678A (zh) * 2011-07-19 2011-11-16 北京师范大学 一种基于沉积动力学反演的储集砂体分析方法
GB2511228B (en) 2011-12-20 2017-06-14 Shell Int Research Automated calibration of a stratigraphic forward modelling (sfm) tool using a neighborhood algorithm with explicit escape clauses
CA2911247C (en) * 2013-08-13 2023-01-10 Schlumberger Canada Limited Digital core sensitivity analysis
CN103454677B (zh) * 2013-09-18 2016-03-02 中海石油(中国)有限公司湛江分公司 基于粒子群与线性加法器结合的地震数据反演方法
US9430499B2 (en) 2014-02-18 2016-08-30 Environmental Systems Research Institute, Inc. Automated feature extraction from imagery
EP2930543B1 (de) 2014-04-10 2021-09-15 Müller-Elektronik GmbH Fallrohrsensor und Verfahren zur Einzelkornerkennung
US9626771B2 (en) 2015-05-20 2017-04-18 Saudi Arabian Oil Company Image-based analysis of a geological thin section
US10614063B2 (en) 2015-10-01 2020-04-07 Microsoft Technology Licensing, Llc. Streaming records from parallel batched database access
WO2017188858A1 (en) * 2016-04-28 2017-11-02 Schlumberger Canada Limited Reservoir performance system
US10713398B2 (en) * 2016-05-23 2020-07-14 Saudi Arabian Oil Company Iterative and repeatable workflow for comprehensive data and processes integration for petroleum exploration and production assessments
CN107085638B (zh) * 2017-04-17 2019-11-12 西南石油大学 一种水力压裂支撑剂参数优化方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DUAN TAIZHONG ED - SONG YAN: "Similarity measure of sedimentary successions and its application in inverse stratigraphic modeling", PETROLEUM SCIENCE, CHINA UNIVERSITY OF PETROLEUM (BEIJING), HEIDELBERG, vol. 14, no. 3, 24 July 2017 (2017-07-24), pages 484 - 492, XP036294635, ISSN: 1672-5107, [retrieved on 20170724], DOI: 10.1007/S12182-017-0174-1 *
HUSSAIN ALQATTAN ET AL: "Moving Models into Reality: An Automated Workflow to Calibrate Forward Stratigraphic Modeling; Application to Hanifa and Arab-D in Central Saudi Arabia", ABU DHABI INTERNATIONAL PETROLEUM EXHIBITION & CONFERENCE, 1 January 2017 (2017-01-01), XP055610443, DOI: 10.2118/188920-MS *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP4016137A1 (en) 2020-12-21 2022-06-22 Abu Dhabi National Oil Company Method and system for determining a location of hydrocarbon reservoirs
US12111436B2 (en) 2020-12-21 2024-10-08 Abu Dhabi National Oil Company Method and system for determining a location of hydrocarbon reservoirs

Also Published As

Publication number Publication date
US20190369278A1 (en) 2019-12-05
US11506804B2 (en) 2022-11-22
EP3803470A1 (en) 2021-04-14
JP2021526634A (ja) 2021-10-07
CA3101976A1 (en) 2019-12-05
CN112689778A (zh) 2021-04-20

Similar Documents

Publication Publication Date Title
US11506804B2 (en) Inverse stratigraphic modeling using a hybrid linear and nonlinear algorithm
RU2669948C2 (ru) Оптимизация многоступенчатого проекта нефтяного месторождения в условиях неопределенности
CN109478208B (zh) 用于石油勘探和生产评估的综合数据和过程集成的迭代且可重复的工作流程
RU2602752C1 (ru) Анализ микросейсмических данных от операции по разрыву пласта
US12182715B2 (en) Proxy modeling workflow and program for reactive transport modeling
Hui et al. A robust embedded discrete fracture modeling workflow for simulating complex processes in field-scale fractured reservoirs
US12189072B2 (en) System and method to identify high-impact discrete fracture model realizations for accelerated calibration of reservoir simulation models
US11493654B2 (en) Construction of a high-resolution advanced 3D transient model with multiple wells by integrating pressure transient data into static geological model
CA3134777A1 (en) Automatic calibration of forward depositional models
US11899150B2 (en) Velocity model for sediment-basement interface using seismic and potential fields data
Zhang et al. 3D parametric modeling of complex geological structures for geotechnical engineering of dam foundation based on t‐splines
US11391856B2 (en) Stochastic dynamic time warping for automated stratigraphic correlation
Koray et al. Improving subsurface characterization utilizing machine learning techniques
US20230393303A1 (en) Integrated diagenetic-depositional facies (iddf) characterization and 3d geomodeling
US11754734B2 (en) Placing wells in a hydrocarbon field based on seismic attributes and quality indicators
Davudov et al. Machine learning assisted well placement optimization
US12385392B2 (en) Determining a three-dimensional fracability index for identifying fracable areas in a subsurface region
US11733415B2 (en) Parallelization of seismic data-related modelling
US20250189687A1 (en) Method for smoothly combining overlapping mistied structural depth grids
US20240240554A1 (en) Training dataset generation process for moment tensor machine learning inversion models
US20230252202A1 (en) Quantification of expressive experimental semi-variogram ranges uncertainties
WO2025068736A1 (en) Method and system for predicting well production rate
US20230168409A1 (en) Hydrocarbon phase behavior modeling for compositional reservoir simulation
Stephen Seismic history matching with saturation indicators combined with multiple objective function optimization
Zhen et al. Method and application of sand body thickness prediction based on virtual sample machine learning

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 19728231

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 3101976

Country of ref document: CA

ENP Entry into the national phase

Ref document number: 2020566913

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2019728231

Country of ref document: EP

Effective date: 20210111