CA2362285C - Uncertainty constrained subsurface modeling - Google Patents
Uncertainty constrained subsurface modeling Download PDFInfo
- Publication number
- CA2362285C CA2362285C CA002362285A CA2362285A CA2362285C CA 2362285 C CA2362285 C CA 2362285C CA 002362285 A CA002362285 A CA 002362285A CA 2362285 A CA2362285 A CA 2362285A CA 2362285 C CA2362285 C CA 2362285C
- Authority
- CA
- Canada
- Prior art keywords
- model
- data
- measurement data
- updating
- subsurface area
- 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.)
- Expired - Lifetime
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. analysis, for interpretation, for correction
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- 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/667—Determining confidence or uncertainty in parameters
Abstract
A method, apparatus, and article of manufacture are provided that use measurement data to create a model of a subsurface area. The method includes creating an initial parameterized model having an initial estimate of model parameter uncertainties; considering measurement data from the subsurface area; updating the model and its associated uncertainty estimate; and repeating the considering and updating steps with additional measurement data. A computer-based apparatus and article of manufacture for implementing the method are also disclosed. The method, apparatus, and article of manufacture are particularly useful in assisting oil companies in making hydrocarbon reservoir data acquisition, drilling and field development decisions.
Description
UNCERTAINTY CONSTRAINED SUBSURFACE MODELING
Related Applications This application claims the benefit of U.S. Provisional Application No. 60/119,821, filed February 12, 1999, and U.S.
Provisional Application No. 60/165,333, filed November 12, 1999.
Background of the Invention This invention is related to subsurface modeling, and is more particularly concerned with a parametric subsurface modeling method, apparatus, and article of manufacture that use uncertainty estimates of subsurface model parameters.
Subsurface models are typically created by geoscientists and engineers to allow development strategies for the subsurface area to be evaluated. Models of this type are commonly created in connection with the development of hydrocarbon reservoirs and mining sites, but they can also used during drilling and related activities where the physical properties of the subsurface area are important. This patent application will focus on the process of creating and updating a model of a subsurface hydrocarbon reservoir, but it should be understood that this merely represents one specific example of how a model of any subsurface area may be created and updated.
Currently, hydrocarbon reservoir modeling is performed most commonly in high-risk, high-profile situations. Typical applications include discoveries in new areas, deepwater exploration, fields in which production surprises or drilling hazards have been encountered, fields in which secondary and tertiary recovery activities are planned, and fields which are being considered for sale or abandonment. The failure to adequately model hydrocarbon reservoirs can have numerous adverse financial consequences, including inaccurate reserve calculations, drilling or completion problems, improper production facility sizing, and suboptimal well placement.
The general problem addressed by this invention is how to construct a model of a subsurface area that is in agreement with multiple sets of measurement data. A model that is in agreement with all of the measurement data obtained from the reservoir can help address. many of the problems noted above.
By 'reservoir model' we mean a quantitative parameterized representation of the subsurface in terms of geometries and material properties. The geometrical model parameters will typically identify geological boundaries, such as contacts between different geologic layers, faults, or fluid/fluid interfaces. The material model parameters will typically identify properties of distributed subsurface materials, such as seismic wave velocities, porosities, permeabilities, fluid saturations, densities, fluid pressures, or temperatures.
By 'agreement' we mean that the data predicted from the reservoir model fit measurements made on the actual reservoir (seismic data, drilling data, well logging data, well test data, production history data, permanent monitoring data, ground penetrating radar data, gravity measurements, etc.).
Virtually all types of measurement data have quantifiable uncertainties and the reservoir model agrees with the measurement data when the difference between data predicted by the reservoir model and measurement data obtained from the reservoir is less than this inherent measurement uncertainty.
While creating a reservoir model that fits one data set is a relatively straightforward task, it is much more difficult to ensure that the model is in agreement with multiple data sets, particularly if the data sets consist of different types of data.
A reservoir model, however, is nonunique even if it is made to fit a variety of data, because different values of material properties and geometries within the model can result in similar predicted measurement values. In other words, the reservoir model has inherent uncertainties: each of the numerical parameters in the reservoir model (e. g., values of material properties within a layer) can take a range of values while the model remains in agreement with the data. This range in parameter values is the uncertainty associated with the reservoir model. The invention described herein is a method to integrate information from multiple measurements and to obtain a reservoir model with quantitative uncertainties in the model parameters. A model of the reservoir that fits the data and has quantified uncertainties can be used to assess the risk inherent in reservoir development decisions (e.g., deciding on the location of additional wells) and to demonstrate the value of additional measurements by showing how these measurements decrease uncertainties in model parameters of interest (e. g., the location of a drilling target or hazard).
A Shared Earth Model (SEM) is a geometrical and material property model of a subsurface area. The model is 'shared' in the sense that it integrates the work of several experts (geologists, geophysicists, well log analysts, reservoir engineers, etc.) who use information from a variety of measurements and interact with the model through different application programs. Ideally, the SEM contains all available information about a reservoir and thus is the basis to make forecasts and plan future actions.
Yet, in any practical case, the information in the measurements is not sufficient to uniquely constrain the parameters (geometries and material properties) of a SEM. As noted above, any SEM has an associated uncertainty, defined here as the range that model parameters can take while fitting available measurements.
Related Applications This application claims the benefit of U.S. Provisional Application No. 60/119,821, filed February 12, 1999, and U.S.
Provisional Application No. 60/165,333, filed November 12, 1999.
Background of the Invention This invention is related to subsurface modeling, and is more particularly concerned with a parametric subsurface modeling method, apparatus, and article of manufacture that use uncertainty estimates of subsurface model parameters.
Subsurface models are typically created by geoscientists and engineers to allow development strategies for the subsurface area to be evaluated. Models of this type are commonly created in connection with the development of hydrocarbon reservoirs and mining sites, but they can also used during drilling and related activities where the physical properties of the subsurface area are important. This patent application will focus on the process of creating and updating a model of a subsurface hydrocarbon reservoir, but it should be understood that this merely represents one specific example of how a model of any subsurface area may be created and updated.
Currently, hydrocarbon reservoir modeling is performed most commonly in high-risk, high-profile situations. Typical applications include discoveries in new areas, deepwater exploration, fields in which production surprises or drilling hazards have been encountered, fields in which secondary and tertiary recovery activities are planned, and fields which are being considered for sale or abandonment. The failure to adequately model hydrocarbon reservoirs can have numerous adverse financial consequences, including inaccurate reserve calculations, drilling or completion problems, improper production facility sizing, and suboptimal well placement.
The general problem addressed by this invention is how to construct a model of a subsurface area that is in agreement with multiple sets of measurement data. A model that is in agreement with all of the measurement data obtained from the reservoir can help address. many of the problems noted above.
By 'reservoir model' we mean a quantitative parameterized representation of the subsurface in terms of geometries and material properties. The geometrical model parameters will typically identify geological boundaries, such as contacts between different geologic layers, faults, or fluid/fluid interfaces. The material model parameters will typically identify properties of distributed subsurface materials, such as seismic wave velocities, porosities, permeabilities, fluid saturations, densities, fluid pressures, or temperatures.
By 'agreement' we mean that the data predicted from the reservoir model fit measurements made on the actual reservoir (seismic data, drilling data, well logging data, well test data, production history data, permanent monitoring data, ground penetrating radar data, gravity measurements, etc.).
Virtually all types of measurement data have quantifiable uncertainties and the reservoir model agrees with the measurement data when the difference between data predicted by the reservoir model and measurement data obtained from the reservoir is less than this inherent measurement uncertainty.
While creating a reservoir model that fits one data set is a relatively straightforward task, it is much more difficult to ensure that the model is in agreement with multiple data sets, particularly if the data sets consist of different types of data.
A reservoir model, however, is nonunique even if it is made to fit a variety of data, because different values of material properties and geometries within the model can result in similar predicted measurement values. In other words, the reservoir model has inherent uncertainties: each of the numerical parameters in the reservoir model (e. g., values of material properties within a layer) can take a range of values while the model remains in agreement with the data. This range in parameter values is the uncertainty associated with the reservoir model. The invention described herein is a method to integrate information from multiple measurements and to obtain a reservoir model with quantitative uncertainties in the model parameters. A model of the reservoir that fits the data and has quantified uncertainties can be used to assess the risk inherent in reservoir development decisions (e.g., deciding on the location of additional wells) and to demonstrate the value of additional measurements by showing how these measurements decrease uncertainties in model parameters of interest (e. g., the location of a drilling target or hazard).
A Shared Earth Model (SEM) is a geometrical and material property model of a subsurface area. The model is 'shared' in the sense that it integrates the work of several experts (geologists, geophysicists, well log analysts, reservoir engineers, etc.) who use information from a variety of measurements and interact with the model through different application programs. Ideally, the SEM contains all available information about a reservoir and thus is the basis to make forecasts and plan future actions.
Yet, in any practical case, the information in the measurements is not sufficient to uniquely constrain the parameters (geometries and material properties) of a SEM. As noted above, any SEM has an associated uncertainty, defined here as the range that model parameters can take while fitting available measurements.
The invention has two primary aspects. The first aspect is a method to quantify and update model parameter uncertainties based on available measurements. One embodiment of this method is based on Bayes' rule, with SEM uncertainty quantified by a posterior probability density function (PDF) ~of the model parameters, conditioned on the measurements used to constrain the model. This posterior PDF may be approximated by a multivariate normal distribution, which is fully described by the posterior mean and covariance matrix of the SEM parameters. Alternatively, one can use a Monte Carlo method to obtain a sample of models drawn from the posterior PDF. This sample of models spans the uncertainty implied by the measurements.
The second aspect is how such a measure of uncertainty acts as a 'memory' of the SEM and can be used for consistent model updating. Quantified uncertainties provide a mechanism to ensure that updates of the SEM based on new data (e. g., well data) are consistent with information provided by data examined previously (2.g., surface seismic data). In particular, we show through a simple example how the effects of a local update of the model can be propagated using the posterior covariance matrix of the SEM parameters. We also show how to update a sample of models obtained by the Monte Carlo method to include new information.
The ideal of a SEM is that all specialists should be able to interact with a common geometry and material property model of the reservoir, incorporating changes into the model using measurements from their own domain of expertise, while maintaining model consistency with previous measurements.
This SEM representation would always be consistent with all available information and should be easy to update as soon as new measurements become available (e. g., from additional wells). Model building would not be a task done episodically, but instead the reservoir model would evolve incrementally as more and more information becomes available during development and production.
V~Thile acquiring more measurements can reduce uncertainty, it is important to weigh the cost of data acquisition against the benefits of reducing uncertainty. This can be done using the tools of decision theory, where different decisions are compared given their associated gains/costs and current uncertainties. A consistent quantification of uncertainties can assist oil companies in making data acquisition, drilling, or development decisions.
Currently, reservoir models are simply modified to fit new data and confirming that the modification is not inconsistent with the previously obtained measurement data is left up to the discretion of the user. The reservoir model may be the result of years of effort and may incorporate measurement data from a wide variety of sources. A user will often only confirm that the change made is not inconsistent with the measurement data within his or her area of expertise (a well log analyst may confirm, for instance, that the change made is consistent with the other well logging data, but may not determine whether the change has introduced an inconsistency with the seismic or geologic data from the area). Many reservoir simulations rely heavily on production data from wells and only four types of geological or geophysical reservoir information: structure of the top of the reservoir, reservoir thickness, porosity, and the ratio of net pay to gross pay. These maps are often constructed from seismic and well log data alone. Incorporating all available data, such as core analyses, seismic-guided reservoir property distributions and fluid analyses, and making certain that the reservoir model is consistent with these different types of data is a cost-effective way to stregthen and validate reservoir models across disciplines.
An iterative method to obtain a model that fits some of the data has been described by George W. Celniker in commonly-assigned U.S. Patent No. 5,905,657, issued March 18, 1999 and incorporated herein by reference. In the Celniker method, the user examines the difference between predicted and measured data, modifies the model attempting to decrease this difference, and repeats the procedure until the fit is satisfactory. This procedure may be adequate if all data sets are considered simultaneously, which may be impractical for diverse and large data sets. If instead the model is modified to fit each of the N data sets in turn (say, from ~(~i to t~(~' ), there is no guarantee that the modifications made to fit the ith data set d~~#do not make the model inconsistent with any of the data sets examined previously ( ~~~' , ~~Z) , . - - , d~~'~?). The model can be assured to be consistent with all data sets only by repeating the comparisons with each data set. Also, if a new data set is acquired and the model is modified to fit it, all other data sets must be examined again to ensure consistency. These repeated checks can make the method time-consuming and inefficient in practice.
Moreover, an iterative comparison of predicted and measured data does not by itself quantify the uncertainty in the model (defined, e.g., as the range that the model parameters can span while still fitting the measured data).
Brief Summary of the Invention The invention comprises a parametric subsurface modeling method, apparatus, and article of manufacture that use measurement data to create a model of a subsurface area. The method includes the steps of: creating a parameterized model having an initial estimate of model parameter uncertainties;
considering measurement data from the subsurface area;
updating the model to fit the measurement data, the updated model having an updated estimate of model parameter uncertainties; and repeating the considering and updating steps with additional measurement data. A computer-based apparatus and article of manufacture for implementing the method are also disclosed. The method, apparatus, and article of manufacture are particularly useful in assisting oil companies in making hydrocarbon reservoir data acquisition and field development decisions. Features of the invention, preferred embodiments and variants thereof, possible applications and their advantages will become appreciated and understood by those skilled in the art from the following detailed description and drawings.
Brief Description of the Drawings FIG. 1 is a flowchart showing steps associated with the present method, apparatus, and article of manufacture;
FIG. 2 is a schematic illustration of computer hardware associated with the apparatus and article of manufacture;
FIG. 3 is a related group of diagrams used to describe a first embodiment of the inventive method;
FIG. 4 is a related group of diagrams used to describe the first embodiment of the inventive method;
FIG. 5 is a related group of diagrams used to describe a second embodiment of the inventive method.
Detailed Description of the Invention Figure 1 shows several steps associated with the present method, apparatus and article of manufacture and provides a general overview of the invention. In the Create Subsurface Model Step 10, a parameterized model of a subsurface area is created using initial reservoir information 12 and/or initial measurement data 14. The initial parameterized model will have an associated initial estimate of model parameter uncertainties.
The model of the subsurface area will typically have geometrical model parameters representing geological boundaries and material parameters representing properties of distributed subsurface materials. The model of the subsurface area may be, for instance, a layered medium representing a layered earth with material properties that are constant or variable within each layer; a geocellular model having material property values defined on a regular or irregular three-dimensional grid; or may be a geometry-based model having material property values defined on a plurality of discrete geometrical sub-regions within the subsurface area.
The initial information may consist of prior knowledge of the spatial distribution of material properties in the subsurface, e.g., the increase of seismic velocity with depth. The initial information may come from physical laws or measurements made in subsurface areas other than the one being modeled. The initial measurement data may consist of seismic data, drilling data, well logging data, well test data, production history data, permanent monitoring data, ground penetrating radar data, gravity measurements, etc. or various combinations of these types of data.
The initial estimate of model parameter uncertainties will typically consist of probability density functions, and preferably consist of multivariate normal/lognormal probability density functions definable by mean vectors and covariance matrices.
In the Consider Measurement Data Step 16, Measurement Data 18 that provides information regarding the subsurface model parameters is examined. In certain cases, the Measurement Data 18 will provide information that may be compared directly to one or more of the subsurface model parameters.
Well logging data may, for instance, give direct measurements of the thickness or the compressional seismic wave velocity of a given geologic layer. In other cases, a Prediction Algorithm 20 will be used to compute data predicted by the model that can be compared to Measurement Data 18, because the Measurement Data 18 only indirectly measures one or more of the subsurface model parameter values. Identify Inconsistencies 22 is used to flag those occasions when the Measurement Data 18 being considered is inconsistent with the model and its associated uncertainty estimate.
In the Update Model and Uncertainty Estimate Step 24, the model is updated to fit the Measurement Data 18 (typically within the model uncertainty constraints). The estimate of model parameter uncertainties is updated as well using the Measurement Data 18. The Consider Measurement Data Step 16, may for instance, produce a likelihood function that is combined with the initial estimate of model parameter uncertainties to produce both the updated parameterized model as well as the updated estimate of model parameter uncertainties. For example, the mean and covariance matrix of the updated estimate of model parameter uncertainties may be computed using deterministic optimization. Alternatively, one may use Monte Carlo sampling to obtain a number of models that are consistent with the Measurement Data 18. These samples may be used to compute the mean and covariance matrix of the updated model parameter uncertainties.
The Consider Measurement Data Step 16 and the Update Model and Uncertainty Estimate Step 24 are repeated for Additional Measurement Data 26 to produce a further updated parameterized model having a further updated estimate of model parameter uncertainties.
The Measurement Data 18 and the Additional Measurement Data 26 may, for instance, consist of different types of data, such as seismic data and well logging data. The Measurement Data 18 and the Additional Measurement Data 26 may, alternatively, consist of the same type of data that has been acquired from the subsurface area at different times to measure changes in reservoir, such as time-lapse/4D surface seismic data.
At least two additional alternative methods for managing the model update process are possible. In some cases, it may be preferable to allow the model parameters to be changed only when the change is consistent with the earlier estimate of model parameter uncertainties. In other cases, Identify Inconsistencies 22 may be used to identify changes in the model parameters that appear to be inconsistent with the initial estimate of model parameter uncertainties derived from the previously considered measurement data. This is generally an indication that model assumptions and data quality assumptions need to be re-examined.
The Consider Measurement Data Step 16 and the Update Model and Uncertainty Estimate Step 24 may be repeated as desired using Additional Measurement Data 26 to further update the model of the subsurface area and its associated estimate of model parameter uncertainties. This may be repeated, for instance, whenever a new set of Additional Measurement Data 26 from the subsurface area becomes available.
Figure 2 schematically illustrates computer hardware that may be used to implement the inventive method. Computer 30 has a media reading device, such as floppy disk device 32, a CD-ROM
Reader or a ZIP drive. The media reading device may also be capable of recording the output of the program the computer 30 is running. A user of the computer 30 may enter commands using a user input device, such as a keyboard 34 or a mouse, may view output of the program code on a visual display device, such as monitor 36, and may make hardcopies of output using an output device, such as printer 38. When properly configured, computer 30 (and its associated peripheral devices) is an apparatus for creating a model of a subsurface area in accordance with the present invention. Computer media, such as floppy disk 40, a CD-ROM or a ZIP disk, may have computer readable program code that allows the computer 42 to create a model of a subsurface area in accordance with the inventive method.
A preferred embodiment of the inventive method will now be described in substantially more detail. The inventive method addresses two primary issues: how to quantify uncertainties in a SEM given measurements and how to use these uncertainties to ensure consistent model updating. The latter is an important issue because in a SEM environment one should be able to continuously update the model; however, model updates based on a set of new data must be consistent with the information provided by data examined previously.
We will now show how to generally address these issues using two simple examples where model uncertainties are calculated and updated using seismic and well data.
The first example will be used to illustrate the quantification of uncertainty in a multivariate normal distribution and consistent model updating environment and will use a simple two-dimensional SEM containing three layers (see Figure 3, diagram 42). This model has seven parameters:
the thicknesses hi of the two top layers at two locations define the SEM geometry, and three compressional wave velocities VPi are material properties. We use layer thicknesses here rather than the absolute depths of horizons because thicknesses are better suited for optimization.
Uncertainties in layer thickness can be easily transformed to uncertainties in absolute depths.
We suppose that at the outset there are measured surface seismic data in the form of four traces recorded at different locations (Figure 3, diagram 44). These four traces were computed by convolving a 50-Hz seismic wavelet with the normal-incidence reflection coefficients of a 'true' model.
The layer thicknesses and velocities in the true model are in parentheses in diagram 42. The densities needed to compute the reflection coefficients were obtained using Gardner's law, which gives rock density as a function of compressional wave velocity. To make the data more realistic, we added band-limited noise for a signal-to-noise ratio of 10.
The second example, see Figure 5, consists of a layered model and will be used to illustrate the quantification of uncertainty using a Monte Carlo method. The data consists of a single seismic reflection trace, and the uncertainty quantification problem is to infer how much variability in the layered medium parameters is consistent with the data given an assumed signal-to-noise ratio.
Quantifying SEM uncertainty from measurements To quantify uncertainty, we use the Bayesian approach widely adopted in the statistical and geophysical inversion literature. Denote the SEM parameters with a vector 111 = [m1, m2, ... , mM]T and the seismic data with a vector d = [dl, d2, ... , dN]T (we use bold lowercase letters for vectors and bold uppercase letters for matrices; the superscript T denotes the transpose). The final uncertainty is quantified by a posterior probability density function (PDF) of the model parameters in 111, which can be obtained from Bayes' rule as the product of a prior PDF and a likelihood function (Figure 4, diagrams 60-74):
p(111 I d, I) G~ p(111 I I) L(111 I d, I) , where p(111 ~ d, I) is the posterior PDF (the PDF of the SEM
parameters 111 given the data d, p(111 ~ I) is the prior PDF, L(1'1'1 I d, I) is the likelihood function, and c~ denotes proportionality. It might seem an impossible task to write an expression for the model PDF given the poor state of knowledge of the reservoir. This is possible by making the PDFs conditional on some prior information I, which includes the parametric form of the model, the noise model, and the accuracy of the physics and simulators used to predict the measurements for a given value of the model parameters.
The prior PDF quantifies what is known about the model parameters from the prior information only, i.e., independently of the information provided by the data. For the initial SEM of diagram 42, we used as a prior PDF a uniform distribution for the layer thicknesses (between a minimum of 1 m and a maximum of 400 m) and a normal distribution for the compressional wave velocities (with a mean of 2500 m/s, a standard deviation of 500 m/s, a minimum of 1500 m/s and a maximum of 5000 m/s). Prior uncertainties in the SEM parameters were also assumed to be independent, so that the prior PDF of all parameters is just the product of the prior PDFs of each:
p (m ~ I) - p (ml ~ I) P (.mz ~ I) ... p (m~ ~ I) This prior PDF represents an initial state of information where the layer thicknesses are unknown, while the prior PDF
of the velocities reflects what is expected for sedimentary rocks.
Information provided by measurements is quantified by the likelihood function, which is formally the probability of observing the data d when the model parameters equal m. For the data of diagram 44, we assumed from prior information that the measured data contained measurement errors that were normally distributed with a standard deviation o'. The likelihood function is then L(m ~ d, I) a exp{-f1/ (2~» fd - g(m) lT fd - g(m» }, where g(m) is a forward modeling (i.e. simulation, prediction) operator that returns the value of the data calculated for a given value of the SEM parameters. In our case, this operator gives the data computed by convolving a seismic wavelet (assumed known) with the reflection coefficient sequence corresponding to the parameters in the SEM. The combination of the prior and likelihood tells us what we know about the model parameters a posteriori.
Diagrams 60 through 74 show how to update a posterior PDF
when new data become available. In the first row, a posterior PDF is obtained when there are only data d~'~. When more data become available, indicated by d~2~, the posterior PDF as constrained by all the data can be obtained by applying Bayer' rule again while using as the prior PDF the posterior obtained previously. This simple update can be done if errors in the data d~'~ and d~2~ are independent, which is a reasonable assumption for data of different types.
Diagrams 60 through 74 also illustrate how uncertainty can be quantified. At any stage in the process, the posterior PDF
can be approximated by a multivariate normal distribution as in p(m ~ d, z) - const. expf-(1/2) (m - ~.)T C' (m - ~) 7 .
In this approximation, the posterior PDF is fully described by the mean value of the model parameters ~.,~, and by the covariance matrix C. The mean (indicated by the white triangles 62, 68, and 76 in diagrams 60, 66, and 74) gives a most probable, 'best' value of the SEM parameters; the covariance matrix defines an ellipsoid in parameter space describing the shape of the posterior PDF (the ellipses 70 and 78 in diagrams 66 and 74), an indication of model parameter uncertainty.
A general practical recipe to evaluate ~u and C uses nonlinear optimization of an objective function x(111) - log p(t11 I d, I) .
A generic optimizer (e. g., quasi-Newton) can be used to find the maximum of the objective function, and the value of m at the maximum can be taken to correspond to the posterior mean The posterior covariance matrix C can be computed as the inverse of the Hessian matrix (the matrix of second derivatives) of the objective function evaluated by finite differences around ~.
The result of applying the nonlinear optimization procedure to the SEM and data of diagrams 42 and 44 is illustrated in diagrams 46 and 48. While the seismic data predicted by the best values of the SEM match closely the measured data, there are large uncertainties in the layer thicknesses and wave velocity in the top layer. This illustrates the fundamental non-uniqueness of time-depth conversion: increasing simultaneously thickness and wave velocity in the top layer does not change the travel time of the reflection. The data do not constrain well the combination of layer thickness and wave velocity, and this uncertainty has been captured in the posterior covariance.
Generic optimization algorithms are typically 'local', in the sense that they find a maximum by moving toward higher values of the objective function from a starting point. Therefore, if the objective function has multiple maxima (as often is the case for band-limited seismic data), the optimizer may converge to a meaningless local maximum. If there are multiple local maxima, these optimizers will converge to the global maximum only if they start from a value of m that is close to (in the sense of being downhill from) that maximum.
For these optimizers to be useful in practice, the user should have the capability to search for a reasonably good starting point by trial-and-error interactions with the SEM.
The uncertainties computed from the Hessian matrix are also 'local' because they are obtained from the local curvature of the objective function near its maximum. The uncertainties computed in this fashion will be accurate only if the objective function is well approximated by a quadratic, i.e., if the posterior PDF is well approximated by a multivariate normal distribution. An alternative is to use a Monte Carlo sampling strategy where values of the model parameter vector 111 are sampled from the posterior PDF. While uncertainties computed from the Hessian matrix are likely to be useful in many instances, there may be cases where a Monte Carlo approach is necessary to obtain a sufficiently accurate uncertainty quantification.
The Monte Carlo approach is illustrated in the example of Figure 5. The data d 86 are seismic reflection data, and the model parameters in the vector m are m = (n, t, v) , where n is the number of layers, t is a vector of travel times to the layer interfaces, and v is a vector of compressional wave velocities or acoustic impedances in each of the layers.
The three profiles 88 are obtained by sampling the posterior PDF of the parameters using a Monte Carlo algorithm. In practice, this can be done in two steps: first obtain a sample of layered media in travel time from the posterior PDF, and then convert each of the sampled layered media from travel time to depth starting from a travel time-depth tie (a point where the travel time and depth are known or have been measured). The image 90 is obtained by superimposing a large number of layered media sampled from the posterior PDF and gives an image of the uncertainty in compressional wave velocity with depth for a given travel time-depth tie.
A method that may be used to obtain a sample from the posterior PDF is the Metropolis-Hastings algorithm. Each step of the algorithm consists of choosing a "candidate" layered medium by perturbing the current one (e.g., by adding a layer, deleting a layer, or changing the travel time to a layer interface). This amounts to choosing a candidate parameter vector m' from a "candidate PDF" q(m'~m), that is the distribution of possible candidates given the current value m of the parameter vector. In the Metropolis-Hastings algorithm, the candidate is accepted with probability alpha = min {1, [p (m' ~d, I) /p(m~d, I) ] [q(m~m') / q(m' ~m) ] } .
It can be shown that a sequence of models obtained with this simple algorithm will eventually sample the posterior PDF
p(mld, I). Once a sample of layered media is obtained, it is easy to convert travel time to depth by starting from a known travel time-depth tie and computing the thickness of the i-th layer hi as h~ = delta ti v~ /2, where deltati is the thickness of the i-th layer in travel time and vi is the velocity of the i-th layer.
The Monte Carlo approach provides a more detailed and accurate quantification of uncertainty compared to the multivariate normal distribution method described above. A
Monte Carlo approach such as that shown in Figure 5 accounts for the possibility of having different number of layers in the reservoir model and captures the uncertainty of posterior PDFs that are not well approximated by a normal distribution.
For example, Figure 5 shows that the PDF of compressional wave velocity at a given depth may be multimodal. If a description of uncertainty in terms of a multivariate normal distribution is needed, however, it is easy to compute a posterior mean ~..~ and a posterior covariance matrix C from the result of Monte Carlo sampling.
Using SEM Uncertainty for Consistent Model Updating The uncertainties computed and displayed in diagram 46 are obviously useful in that they quantify how well the seismic data constrain each parameter of the SEM. Uncertainties also provide a 'memory' to the model: if information on the covariance matrix were stored with the model, it would effectively remember the constraints placed by all the data types that the model was tested against.
The most general way to update the model and its uncertainty when new data become available is illustrated in diagrams 60 through 74, and consists of a repeated application of Bayes' rule. In this section we describe two methods that can be used to quickly accomplish model updating and that can be used when new data becomes available (e.g., from a well) to update a SEM whose uncertainty has been quantified using a multivariate normal distribution or the Monte Carlo approach.
An important use of quantified uncertainty is in propagating model updates to ensure consistency of the SEM with previously examined data. We illustrate the first method for consistent model updating, which is applicable when the uncertainty is quantified in a multivariate normal PDF form, by considering what happens when new information from a well becomes available. The simplest update is illustrated in diagrams 50 and 52: when observations from a well drilled at the left edge of the SEM show that the thicknesses of the top two layers are 200 m and 20 m, the model is simply updated to reflect the new information.
If layer thicknesses in the model of diagrams 46 and 48 are changed to reflect values observed at a well without changing any other SEM parameters, the predicted seismic data do not match anymore the measured data. To ensure consistency, one would have to match the model to the seismic data again after the update.
In the case of time-depth conversion this is not unrealistic, and the consequences of changes in layer thicknesses can be calculated in a straightforward way. But in the general case where a SEM has been constrained by a variety of data, there is no obvious way to compute how the effects of a local update should be propagated to maintain consistency.
We now describe a simple mechanism where quantified uncertainties are used to propagate the effects of a local model update to ensure consistency of a SEM with data examined previously. The mechanism uses the transmission of information provided by Bayes' rule (diagrams 60 through 74) and is illustrated in diagram 80. After considering only the seismic data, diagram 80 shows as a shaded ellipse the posterior PDF of the compressional wave velocity VP and thickness h of a layer; the best value of VP and h is shown as a white triangle 82. Suppose now that, as shown in diagrams 50 and 52, new information from a well becomes available and it is now known that the layer thickness must be within a narrow interval, indicated by the narrow cross-hatched rectangle in chart 80. The best value for VP and h given the well data is now located at dot 84, which is the point within the narrow rectangle where the posterior PDF
given the seismic data only is greatest. Note that the information from the seismic data specifies that if the best value of the layer thickness changes, the best value of the layer wave velocity must change as well.
At the limit where the layer thickness is fixed with absolute certainty by the well observations (i.e., the rectangle in chart 80 shrinks to a vertical line), the operation illustrated in chart 80 is the calculation of the conditional mean of VP for a given value of h. In the general case where one wants to compute the mean of some model parameters conditional on other model parameters that have been updated, this calculation is straightforward. Start by rearranging the vector of SEM parameters t11 into two vectors ml m= , m2 where the vector m~ contains the parameters whose conditional mean is to be computed, and the vector 1112 contains the parameters that are fixed. In the example of Figure 3, diagrams 50 and 52, 1112 would contain the two layer thicknesses constrained by the well observations and 111 the two other layer thicknesses and the three wave velocities.
After considering some data (seismic data in the example of Figure 3, diagrams 50 and 52), suppose we have a multivariate normal posterior PDF for 111 with mean ~ and covariance C.
Rearrange the mean vector as done for 11'1 in two vectors ~~ and ~,2 and the covariance matrix as C - Cm Ciz T
Ciz C2z It can be shown that the mean of the parameters in conditional on the other parameters having a value m2 is and the conditional covariance of the parameters in 111 is -, T
" - C,2 C22 ,2 This approach can be easily applied to the example of Figure 3, diagrams 50 and 52. Figure 3, diagrams 54 and 56 show the result of adding the well data using a 'consistent' model update, i.e., computing the mean and covariance of the seismically-defined SEM parameters conditional on the two values of thickness observed at the well. In simple terms, this is a consistent update because it uses the information provided by the seismic data to propagate the effects of a local model update done on the basis of well data to other model parameters in order to preserve consistency with the seismic data. For example, the seismic data imply that if the thickness of the top layer is modified, the wave velocity must change as well. The conditional value computed from the uncertainty for the top layer wave velocity is 1997 m/s, which is very close to the true value (2000 m/s; see diagrams 42 and 54).
Diagrams 54 and 56 also show that the data predicted by the updated SEM match the measurements well, in contrast with the SEM of diagrams 50 and 52. It is important to stress that this close match has not been achieved by comparing the SEM
predictions with the seismic data again after the well data were incorporated in the model. The good fit to the measurements is simply a result of using the posterior PDF to compute a consistent model update. This is the key advantage of quantifying uncertainties; at first approximation, the effects of any update to the SEM can be propagated using the posterior uncertainties without needing to re-examine all data previously incorporated into the model. The posterior uncertainties are a 'memory' mechanism that allows the SEM to remember how closely data constrain its geometries and material property distributions.
The limitations of this mechanism follow from the limitations on the quantification of uncertainty noted in the previous section. If the multivariate normal representation we use here does not approximate the posterior PDF closely, the computed consistent update will suffer. In our example, it is clear that while the seismic data predicted by the consistently updated SEM are close to the measurements, the fit may be improved (compare diagrams 54 and 56 with diagrams 46 and 48). Also, some uncertainties are smaller than they should be; for example, the uncertainty in the thickness of the second layer should be greater than 0.8 m.
On the other hand, the updated model should be close to the best value and thus an automated optimization applied at this stage should have a good chance of succeeding. Starting from the model of diagrams 50 and 52, an optimization converges to a meaningless local maximum; if one starts instead from the consistently updated model of diagrams 54 and 56, the optimization quickly converges to a best model with an improved fit to the seismic data. Once the best model is found, the posterior covariance matrix may be evaluated again for more accurate uncertainty quantification.
In the previous section we also described a Monte Carlo approach to quantify uncertainty. If the model uncertainty is quantified in a multivariate normal distribution whose mean ~..~, and covariance matrix C are evaluated from the result of Monte Carlo sampling, the consistent update method we have illustrated above can be used to compute a quick update, i.e.
without requiring the Monte Carlo sampling process to be repeated. We now illustrate a second method that uses additional information to directly update the sample from the posterior PDF obtained by a Monte Carlo method.
The method we describe can be applied when the additional information that becomes available directly controls the properties of the posterior PDF sample. In the example of Figure 5, diagrams 86-90, a Monte Carlo method was used to obtain a sample from the posterior PDF of the thicknesses and compressional wave velocities in a layered medium. Because the seismic reflection data are in travel time, obtaining a layered medium description in depth requires having a travel time-depth tie, that is, having a point where the absolute depth and the corresponding travel time are known or measured. This travel time-depth tie may be the Earth's surface, where depth and travel time are zero. For offshore reservoirs, the travel time-depth tie may be at the sea bottom. If the seismic data were acquired using receivers or sources in a borehole, the travel time-depth tie may be obtained at the depth of a receiver/source in the borehole.
As will be clear to those skilled in the art, the uncertainty in the depth of a layer as determined from seismic reflection data is greater the farther. away the layer is from the deepest travel time-depth tie. This is because depth is the sum of layer thicknesses starting from the deepest travel time-depth tie; the uncertainty in depth is the cumulative uncertainty of the layer thicknesses.
Diagrams 92 and 94 illustrate our method of updating a layered model and its uncertainties as new information on travel time-depth ties becomes available. Diagram 92 shows the results of Monte Carlo sampling obtained for a relatively shallow travel time-depth tie (e.g. the Earth's surface, the sea floor, the depth of a seismic receiver/source in a borehole). The results are presented in the form of a posterior image (as in diagram 90) computed by superimposing all the layered media sampled by the Monte Carlo method. As noted in the previous section, in practice the Monte Carlo sampling obtains a sample of layered media in travel time, rather than depth; each of these layered media is then converted to depth using a known/measured travel time-depth tie.
Diagram 94 illustrates how the uncertainty in depth changes if a deeper travel time-depth tie is measured. This additional travel time-depth tie will typically be provided by a seismic receiver/source placed in a borehole, and it could be acquired during or after drilling the well. Diagram 94 shows how the uncertainty in deep layers is reduced by using the deeper travel time-depth tie. It should be stressed that the application of this method does not require sampling the posterior PDF again, but simply using the same sample of layered media in travel time obtained earlier while recomputing the travel time to depth conversion for the additional, deeper travel time-depth tie. This reduction in uncertainty can be extremely important in making decisions during the drilling of a well (e. g., on determining the depth of zones that have anomalous pore pressures and are drilling hazards).
Uncertainty quantification and consistent model updating can improve significantly the efficiency of constructing and modifying a SEM. Because industry-standard interpretation workflows don't account for uncertainty, model consistency is maintained by making elements of the model 'rigid' as the interpretation progresses down the workflow. The term 'rigid' means here that once a domain expert has set the optimal values for parameters under his control (e. g., a geophysicist interpreting the model framework), these values are not changed by later experts (e.g., a reservoir engineer) for fear that the model will no longer be consistent with the previous data. For example, once the model framework is fixed, the reservoir engineer is only left flow parameters to adjust during history matching. Explicitly accounting for uncertainty would allow that same reservoir engineer to adjust all parameters of the model within their acceptable ranges to obtain the best history match. The remaining uncertainty in the model then provides input to decision-making tools which, for example, can be used to estimate a PDF of the net present value of a reservoir. Thus, uncertainty quantification allows us to approach the ideal of a SEM that is constrained by as many data as possible, can be easily updated, and ties directly into decision-making tools.
The foregoing description of the preferred and alternate embodiments of the present invention have been presented for purposes of illustration and description. It is not intended to be exhaustive or limit the invention to the precise examples described. Many modifications and variations will be apparent to those skilled in the art. The embodiments were chosen and described in order to best explain the principles of the invention and its practical application, thereby enabling others skilled in the art to understand the invention for various embodiments and with various modifications as are suited to the particular use contemplated. It is intended that the scope of the invention be defined by the accompanying claims and their equivalents.
The second aspect is how such a measure of uncertainty acts as a 'memory' of the SEM and can be used for consistent model updating. Quantified uncertainties provide a mechanism to ensure that updates of the SEM based on new data (e. g., well data) are consistent with information provided by data examined previously (2.g., surface seismic data). In particular, we show through a simple example how the effects of a local update of the model can be propagated using the posterior covariance matrix of the SEM parameters. We also show how to update a sample of models obtained by the Monte Carlo method to include new information.
The ideal of a SEM is that all specialists should be able to interact with a common geometry and material property model of the reservoir, incorporating changes into the model using measurements from their own domain of expertise, while maintaining model consistency with previous measurements.
This SEM representation would always be consistent with all available information and should be easy to update as soon as new measurements become available (e. g., from additional wells). Model building would not be a task done episodically, but instead the reservoir model would evolve incrementally as more and more information becomes available during development and production.
V~Thile acquiring more measurements can reduce uncertainty, it is important to weigh the cost of data acquisition against the benefits of reducing uncertainty. This can be done using the tools of decision theory, where different decisions are compared given their associated gains/costs and current uncertainties. A consistent quantification of uncertainties can assist oil companies in making data acquisition, drilling, or development decisions.
Currently, reservoir models are simply modified to fit new data and confirming that the modification is not inconsistent with the previously obtained measurement data is left up to the discretion of the user. The reservoir model may be the result of years of effort and may incorporate measurement data from a wide variety of sources. A user will often only confirm that the change made is not inconsistent with the measurement data within his or her area of expertise (a well log analyst may confirm, for instance, that the change made is consistent with the other well logging data, but may not determine whether the change has introduced an inconsistency with the seismic or geologic data from the area). Many reservoir simulations rely heavily on production data from wells and only four types of geological or geophysical reservoir information: structure of the top of the reservoir, reservoir thickness, porosity, and the ratio of net pay to gross pay. These maps are often constructed from seismic and well log data alone. Incorporating all available data, such as core analyses, seismic-guided reservoir property distributions and fluid analyses, and making certain that the reservoir model is consistent with these different types of data is a cost-effective way to stregthen and validate reservoir models across disciplines.
An iterative method to obtain a model that fits some of the data has been described by George W. Celniker in commonly-assigned U.S. Patent No. 5,905,657, issued March 18, 1999 and incorporated herein by reference. In the Celniker method, the user examines the difference between predicted and measured data, modifies the model attempting to decrease this difference, and repeats the procedure until the fit is satisfactory. This procedure may be adequate if all data sets are considered simultaneously, which may be impractical for diverse and large data sets. If instead the model is modified to fit each of the N data sets in turn (say, from ~(~i to t~(~' ), there is no guarantee that the modifications made to fit the ith data set d~~#do not make the model inconsistent with any of the data sets examined previously ( ~~~' , ~~Z) , . - - , d~~'~?). The model can be assured to be consistent with all data sets only by repeating the comparisons with each data set. Also, if a new data set is acquired and the model is modified to fit it, all other data sets must be examined again to ensure consistency. These repeated checks can make the method time-consuming and inefficient in practice.
Moreover, an iterative comparison of predicted and measured data does not by itself quantify the uncertainty in the model (defined, e.g., as the range that the model parameters can span while still fitting the measured data).
Brief Summary of the Invention The invention comprises a parametric subsurface modeling method, apparatus, and article of manufacture that use measurement data to create a model of a subsurface area. The method includes the steps of: creating a parameterized model having an initial estimate of model parameter uncertainties;
considering measurement data from the subsurface area;
updating the model to fit the measurement data, the updated model having an updated estimate of model parameter uncertainties; and repeating the considering and updating steps with additional measurement data. A computer-based apparatus and article of manufacture for implementing the method are also disclosed. The method, apparatus, and article of manufacture are particularly useful in assisting oil companies in making hydrocarbon reservoir data acquisition and field development decisions. Features of the invention, preferred embodiments and variants thereof, possible applications and their advantages will become appreciated and understood by those skilled in the art from the following detailed description and drawings.
Brief Description of the Drawings FIG. 1 is a flowchart showing steps associated with the present method, apparatus, and article of manufacture;
FIG. 2 is a schematic illustration of computer hardware associated with the apparatus and article of manufacture;
FIG. 3 is a related group of diagrams used to describe a first embodiment of the inventive method;
FIG. 4 is a related group of diagrams used to describe the first embodiment of the inventive method;
FIG. 5 is a related group of diagrams used to describe a second embodiment of the inventive method.
Detailed Description of the Invention Figure 1 shows several steps associated with the present method, apparatus and article of manufacture and provides a general overview of the invention. In the Create Subsurface Model Step 10, a parameterized model of a subsurface area is created using initial reservoir information 12 and/or initial measurement data 14. The initial parameterized model will have an associated initial estimate of model parameter uncertainties.
The model of the subsurface area will typically have geometrical model parameters representing geological boundaries and material parameters representing properties of distributed subsurface materials. The model of the subsurface area may be, for instance, a layered medium representing a layered earth with material properties that are constant or variable within each layer; a geocellular model having material property values defined on a regular or irregular three-dimensional grid; or may be a geometry-based model having material property values defined on a plurality of discrete geometrical sub-regions within the subsurface area.
The initial information may consist of prior knowledge of the spatial distribution of material properties in the subsurface, e.g., the increase of seismic velocity with depth. The initial information may come from physical laws or measurements made in subsurface areas other than the one being modeled. The initial measurement data may consist of seismic data, drilling data, well logging data, well test data, production history data, permanent monitoring data, ground penetrating radar data, gravity measurements, etc. or various combinations of these types of data.
The initial estimate of model parameter uncertainties will typically consist of probability density functions, and preferably consist of multivariate normal/lognormal probability density functions definable by mean vectors and covariance matrices.
In the Consider Measurement Data Step 16, Measurement Data 18 that provides information regarding the subsurface model parameters is examined. In certain cases, the Measurement Data 18 will provide information that may be compared directly to one or more of the subsurface model parameters.
Well logging data may, for instance, give direct measurements of the thickness or the compressional seismic wave velocity of a given geologic layer. In other cases, a Prediction Algorithm 20 will be used to compute data predicted by the model that can be compared to Measurement Data 18, because the Measurement Data 18 only indirectly measures one or more of the subsurface model parameter values. Identify Inconsistencies 22 is used to flag those occasions when the Measurement Data 18 being considered is inconsistent with the model and its associated uncertainty estimate.
In the Update Model and Uncertainty Estimate Step 24, the model is updated to fit the Measurement Data 18 (typically within the model uncertainty constraints). The estimate of model parameter uncertainties is updated as well using the Measurement Data 18. The Consider Measurement Data Step 16, may for instance, produce a likelihood function that is combined with the initial estimate of model parameter uncertainties to produce both the updated parameterized model as well as the updated estimate of model parameter uncertainties. For example, the mean and covariance matrix of the updated estimate of model parameter uncertainties may be computed using deterministic optimization. Alternatively, one may use Monte Carlo sampling to obtain a number of models that are consistent with the Measurement Data 18. These samples may be used to compute the mean and covariance matrix of the updated model parameter uncertainties.
The Consider Measurement Data Step 16 and the Update Model and Uncertainty Estimate Step 24 are repeated for Additional Measurement Data 26 to produce a further updated parameterized model having a further updated estimate of model parameter uncertainties.
The Measurement Data 18 and the Additional Measurement Data 26 may, for instance, consist of different types of data, such as seismic data and well logging data. The Measurement Data 18 and the Additional Measurement Data 26 may, alternatively, consist of the same type of data that has been acquired from the subsurface area at different times to measure changes in reservoir, such as time-lapse/4D surface seismic data.
At least two additional alternative methods for managing the model update process are possible. In some cases, it may be preferable to allow the model parameters to be changed only when the change is consistent with the earlier estimate of model parameter uncertainties. In other cases, Identify Inconsistencies 22 may be used to identify changes in the model parameters that appear to be inconsistent with the initial estimate of model parameter uncertainties derived from the previously considered measurement data. This is generally an indication that model assumptions and data quality assumptions need to be re-examined.
The Consider Measurement Data Step 16 and the Update Model and Uncertainty Estimate Step 24 may be repeated as desired using Additional Measurement Data 26 to further update the model of the subsurface area and its associated estimate of model parameter uncertainties. This may be repeated, for instance, whenever a new set of Additional Measurement Data 26 from the subsurface area becomes available.
Figure 2 schematically illustrates computer hardware that may be used to implement the inventive method. Computer 30 has a media reading device, such as floppy disk device 32, a CD-ROM
Reader or a ZIP drive. The media reading device may also be capable of recording the output of the program the computer 30 is running. A user of the computer 30 may enter commands using a user input device, such as a keyboard 34 or a mouse, may view output of the program code on a visual display device, such as monitor 36, and may make hardcopies of output using an output device, such as printer 38. When properly configured, computer 30 (and its associated peripheral devices) is an apparatus for creating a model of a subsurface area in accordance with the present invention. Computer media, such as floppy disk 40, a CD-ROM or a ZIP disk, may have computer readable program code that allows the computer 42 to create a model of a subsurface area in accordance with the inventive method.
A preferred embodiment of the inventive method will now be described in substantially more detail. The inventive method addresses two primary issues: how to quantify uncertainties in a SEM given measurements and how to use these uncertainties to ensure consistent model updating. The latter is an important issue because in a SEM environment one should be able to continuously update the model; however, model updates based on a set of new data must be consistent with the information provided by data examined previously.
We will now show how to generally address these issues using two simple examples where model uncertainties are calculated and updated using seismic and well data.
The first example will be used to illustrate the quantification of uncertainty in a multivariate normal distribution and consistent model updating environment and will use a simple two-dimensional SEM containing three layers (see Figure 3, diagram 42). This model has seven parameters:
the thicknesses hi of the two top layers at two locations define the SEM geometry, and three compressional wave velocities VPi are material properties. We use layer thicknesses here rather than the absolute depths of horizons because thicknesses are better suited for optimization.
Uncertainties in layer thickness can be easily transformed to uncertainties in absolute depths.
We suppose that at the outset there are measured surface seismic data in the form of four traces recorded at different locations (Figure 3, diagram 44). These four traces were computed by convolving a 50-Hz seismic wavelet with the normal-incidence reflection coefficients of a 'true' model.
The layer thicknesses and velocities in the true model are in parentheses in diagram 42. The densities needed to compute the reflection coefficients were obtained using Gardner's law, which gives rock density as a function of compressional wave velocity. To make the data more realistic, we added band-limited noise for a signal-to-noise ratio of 10.
The second example, see Figure 5, consists of a layered model and will be used to illustrate the quantification of uncertainty using a Monte Carlo method. The data consists of a single seismic reflection trace, and the uncertainty quantification problem is to infer how much variability in the layered medium parameters is consistent with the data given an assumed signal-to-noise ratio.
Quantifying SEM uncertainty from measurements To quantify uncertainty, we use the Bayesian approach widely adopted in the statistical and geophysical inversion literature. Denote the SEM parameters with a vector 111 = [m1, m2, ... , mM]T and the seismic data with a vector d = [dl, d2, ... , dN]T (we use bold lowercase letters for vectors and bold uppercase letters for matrices; the superscript T denotes the transpose). The final uncertainty is quantified by a posterior probability density function (PDF) of the model parameters in 111, which can be obtained from Bayes' rule as the product of a prior PDF and a likelihood function (Figure 4, diagrams 60-74):
p(111 I d, I) G~ p(111 I I) L(111 I d, I) , where p(111 ~ d, I) is the posterior PDF (the PDF of the SEM
parameters 111 given the data d, p(111 ~ I) is the prior PDF, L(1'1'1 I d, I) is the likelihood function, and c~ denotes proportionality. It might seem an impossible task to write an expression for the model PDF given the poor state of knowledge of the reservoir. This is possible by making the PDFs conditional on some prior information I, which includes the parametric form of the model, the noise model, and the accuracy of the physics and simulators used to predict the measurements for a given value of the model parameters.
The prior PDF quantifies what is known about the model parameters from the prior information only, i.e., independently of the information provided by the data. For the initial SEM of diagram 42, we used as a prior PDF a uniform distribution for the layer thicknesses (between a minimum of 1 m and a maximum of 400 m) and a normal distribution for the compressional wave velocities (with a mean of 2500 m/s, a standard deviation of 500 m/s, a minimum of 1500 m/s and a maximum of 5000 m/s). Prior uncertainties in the SEM parameters were also assumed to be independent, so that the prior PDF of all parameters is just the product of the prior PDFs of each:
p (m ~ I) - p (ml ~ I) P (.mz ~ I) ... p (m~ ~ I) This prior PDF represents an initial state of information where the layer thicknesses are unknown, while the prior PDF
of the velocities reflects what is expected for sedimentary rocks.
Information provided by measurements is quantified by the likelihood function, which is formally the probability of observing the data d when the model parameters equal m. For the data of diagram 44, we assumed from prior information that the measured data contained measurement errors that were normally distributed with a standard deviation o'. The likelihood function is then L(m ~ d, I) a exp{-f1/ (2~» fd - g(m) lT fd - g(m» }, where g(m) is a forward modeling (i.e. simulation, prediction) operator that returns the value of the data calculated for a given value of the SEM parameters. In our case, this operator gives the data computed by convolving a seismic wavelet (assumed known) with the reflection coefficient sequence corresponding to the parameters in the SEM. The combination of the prior and likelihood tells us what we know about the model parameters a posteriori.
Diagrams 60 through 74 show how to update a posterior PDF
when new data become available. In the first row, a posterior PDF is obtained when there are only data d~'~. When more data become available, indicated by d~2~, the posterior PDF as constrained by all the data can be obtained by applying Bayer' rule again while using as the prior PDF the posterior obtained previously. This simple update can be done if errors in the data d~'~ and d~2~ are independent, which is a reasonable assumption for data of different types.
Diagrams 60 through 74 also illustrate how uncertainty can be quantified. At any stage in the process, the posterior PDF
can be approximated by a multivariate normal distribution as in p(m ~ d, z) - const. expf-(1/2) (m - ~.)T C' (m - ~) 7 .
In this approximation, the posterior PDF is fully described by the mean value of the model parameters ~.,~, and by the covariance matrix C. The mean (indicated by the white triangles 62, 68, and 76 in diagrams 60, 66, and 74) gives a most probable, 'best' value of the SEM parameters; the covariance matrix defines an ellipsoid in parameter space describing the shape of the posterior PDF (the ellipses 70 and 78 in diagrams 66 and 74), an indication of model parameter uncertainty.
A general practical recipe to evaluate ~u and C uses nonlinear optimization of an objective function x(111) - log p(t11 I d, I) .
A generic optimizer (e. g., quasi-Newton) can be used to find the maximum of the objective function, and the value of m at the maximum can be taken to correspond to the posterior mean The posterior covariance matrix C can be computed as the inverse of the Hessian matrix (the matrix of second derivatives) of the objective function evaluated by finite differences around ~.
The result of applying the nonlinear optimization procedure to the SEM and data of diagrams 42 and 44 is illustrated in diagrams 46 and 48. While the seismic data predicted by the best values of the SEM match closely the measured data, there are large uncertainties in the layer thicknesses and wave velocity in the top layer. This illustrates the fundamental non-uniqueness of time-depth conversion: increasing simultaneously thickness and wave velocity in the top layer does not change the travel time of the reflection. The data do not constrain well the combination of layer thickness and wave velocity, and this uncertainty has been captured in the posterior covariance.
Generic optimization algorithms are typically 'local', in the sense that they find a maximum by moving toward higher values of the objective function from a starting point. Therefore, if the objective function has multiple maxima (as often is the case for band-limited seismic data), the optimizer may converge to a meaningless local maximum. If there are multiple local maxima, these optimizers will converge to the global maximum only if they start from a value of m that is close to (in the sense of being downhill from) that maximum.
For these optimizers to be useful in practice, the user should have the capability to search for a reasonably good starting point by trial-and-error interactions with the SEM.
The uncertainties computed from the Hessian matrix are also 'local' because they are obtained from the local curvature of the objective function near its maximum. The uncertainties computed in this fashion will be accurate only if the objective function is well approximated by a quadratic, i.e., if the posterior PDF is well approximated by a multivariate normal distribution. An alternative is to use a Monte Carlo sampling strategy where values of the model parameter vector 111 are sampled from the posterior PDF. While uncertainties computed from the Hessian matrix are likely to be useful in many instances, there may be cases where a Monte Carlo approach is necessary to obtain a sufficiently accurate uncertainty quantification.
The Monte Carlo approach is illustrated in the example of Figure 5. The data d 86 are seismic reflection data, and the model parameters in the vector m are m = (n, t, v) , where n is the number of layers, t is a vector of travel times to the layer interfaces, and v is a vector of compressional wave velocities or acoustic impedances in each of the layers.
The three profiles 88 are obtained by sampling the posterior PDF of the parameters using a Monte Carlo algorithm. In practice, this can be done in two steps: first obtain a sample of layered media in travel time from the posterior PDF, and then convert each of the sampled layered media from travel time to depth starting from a travel time-depth tie (a point where the travel time and depth are known or have been measured). The image 90 is obtained by superimposing a large number of layered media sampled from the posterior PDF and gives an image of the uncertainty in compressional wave velocity with depth for a given travel time-depth tie.
A method that may be used to obtain a sample from the posterior PDF is the Metropolis-Hastings algorithm. Each step of the algorithm consists of choosing a "candidate" layered medium by perturbing the current one (e.g., by adding a layer, deleting a layer, or changing the travel time to a layer interface). This amounts to choosing a candidate parameter vector m' from a "candidate PDF" q(m'~m), that is the distribution of possible candidates given the current value m of the parameter vector. In the Metropolis-Hastings algorithm, the candidate is accepted with probability alpha = min {1, [p (m' ~d, I) /p(m~d, I) ] [q(m~m') / q(m' ~m) ] } .
It can be shown that a sequence of models obtained with this simple algorithm will eventually sample the posterior PDF
p(mld, I). Once a sample of layered media is obtained, it is easy to convert travel time to depth by starting from a known travel time-depth tie and computing the thickness of the i-th layer hi as h~ = delta ti v~ /2, where deltati is the thickness of the i-th layer in travel time and vi is the velocity of the i-th layer.
The Monte Carlo approach provides a more detailed and accurate quantification of uncertainty compared to the multivariate normal distribution method described above. A
Monte Carlo approach such as that shown in Figure 5 accounts for the possibility of having different number of layers in the reservoir model and captures the uncertainty of posterior PDFs that are not well approximated by a normal distribution.
For example, Figure 5 shows that the PDF of compressional wave velocity at a given depth may be multimodal. If a description of uncertainty in terms of a multivariate normal distribution is needed, however, it is easy to compute a posterior mean ~..~ and a posterior covariance matrix C from the result of Monte Carlo sampling.
Using SEM Uncertainty for Consistent Model Updating The uncertainties computed and displayed in diagram 46 are obviously useful in that they quantify how well the seismic data constrain each parameter of the SEM. Uncertainties also provide a 'memory' to the model: if information on the covariance matrix were stored with the model, it would effectively remember the constraints placed by all the data types that the model was tested against.
The most general way to update the model and its uncertainty when new data become available is illustrated in diagrams 60 through 74, and consists of a repeated application of Bayes' rule. In this section we describe two methods that can be used to quickly accomplish model updating and that can be used when new data becomes available (e.g., from a well) to update a SEM whose uncertainty has been quantified using a multivariate normal distribution or the Monte Carlo approach.
An important use of quantified uncertainty is in propagating model updates to ensure consistency of the SEM with previously examined data. We illustrate the first method for consistent model updating, which is applicable when the uncertainty is quantified in a multivariate normal PDF form, by considering what happens when new information from a well becomes available. The simplest update is illustrated in diagrams 50 and 52: when observations from a well drilled at the left edge of the SEM show that the thicknesses of the top two layers are 200 m and 20 m, the model is simply updated to reflect the new information.
If layer thicknesses in the model of diagrams 46 and 48 are changed to reflect values observed at a well without changing any other SEM parameters, the predicted seismic data do not match anymore the measured data. To ensure consistency, one would have to match the model to the seismic data again after the update.
In the case of time-depth conversion this is not unrealistic, and the consequences of changes in layer thicknesses can be calculated in a straightforward way. But in the general case where a SEM has been constrained by a variety of data, there is no obvious way to compute how the effects of a local update should be propagated to maintain consistency.
We now describe a simple mechanism where quantified uncertainties are used to propagate the effects of a local model update to ensure consistency of a SEM with data examined previously. The mechanism uses the transmission of information provided by Bayes' rule (diagrams 60 through 74) and is illustrated in diagram 80. After considering only the seismic data, diagram 80 shows as a shaded ellipse the posterior PDF of the compressional wave velocity VP and thickness h of a layer; the best value of VP and h is shown as a white triangle 82. Suppose now that, as shown in diagrams 50 and 52, new information from a well becomes available and it is now known that the layer thickness must be within a narrow interval, indicated by the narrow cross-hatched rectangle in chart 80. The best value for VP and h given the well data is now located at dot 84, which is the point within the narrow rectangle where the posterior PDF
given the seismic data only is greatest. Note that the information from the seismic data specifies that if the best value of the layer thickness changes, the best value of the layer wave velocity must change as well.
At the limit where the layer thickness is fixed with absolute certainty by the well observations (i.e., the rectangle in chart 80 shrinks to a vertical line), the operation illustrated in chart 80 is the calculation of the conditional mean of VP for a given value of h. In the general case where one wants to compute the mean of some model parameters conditional on other model parameters that have been updated, this calculation is straightforward. Start by rearranging the vector of SEM parameters t11 into two vectors ml m= , m2 where the vector m~ contains the parameters whose conditional mean is to be computed, and the vector 1112 contains the parameters that are fixed. In the example of Figure 3, diagrams 50 and 52, 1112 would contain the two layer thicknesses constrained by the well observations and 111 the two other layer thicknesses and the three wave velocities.
After considering some data (seismic data in the example of Figure 3, diagrams 50 and 52), suppose we have a multivariate normal posterior PDF for 111 with mean ~ and covariance C.
Rearrange the mean vector as done for 11'1 in two vectors ~~ and ~,2 and the covariance matrix as C - Cm Ciz T
Ciz C2z It can be shown that the mean of the parameters in conditional on the other parameters having a value m2 is and the conditional covariance of the parameters in 111 is -, T
" - C,2 C22 ,2 This approach can be easily applied to the example of Figure 3, diagrams 50 and 52. Figure 3, diagrams 54 and 56 show the result of adding the well data using a 'consistent' model update, i.e., computing the mean and covariance of the seismically-defined SEM parameters conditional on the two values of thickness observed at the well. In simple terms, this is a consistent update because it uses the information provided by the seismic data to propagate the effects of a local model update done on the basis of well data to other model parameters in order to preserve consistency with the seismic data. For example, the seismic data imply that if the thickness of the top layer is modified, the wave velocity must change as well. The conditional value computed from the uncertainty for the top layer wave velocity is 1997 m/s, which is very close to the true value (2000 m/s; see diagrams 42 and 54).
Diagrams 54 and 56 also show that the data predicted by the updated SEM match the measurements well, in contrast with the SEM of diagrams 50 and 52. It is important to stress that this close match has not been achieved by comparing the SEM
predictions with the seismic data again after the well data were incorporated in the model. The good fit to the measurements is simply a result of using the posterior PDF to compute a consistent model update. This is the key advantage of quantifying uncertainties; at first approximation, the effects of any update to the SEM can be propagated using the posterior uncertainties without needing to re-examine all data previously incorporated into the model. The posterior uncertainties are a 'memory' mechanism that allows the SEM to remember how closely data constrain its geometries and material property distributions.
The limitations of this mechanism follow from the limitations on the quantification of uncertainty noted in the previous section. If the multivariate normal representation we use here does not approximate the posterior PDF closely, the computed consistent update will suffer. In our example, it is clear that while the seismic data predicted by the consistently updated SEM are close to the measurements, the fit may be improved (compare diagrams 54 and 56 with diagrams 46 and 48). Also, some uncertainties are smaller than they should be; for example, the uncertainty in the thickness of the second layer should be greater than 0.8 m.
On the other hand, the updated model should be close to the best value and thus an automated optimization applied at this stage should have a good chance of succeeding. Starting from the model of diagrams 50 and 52, an optimization converges to a meaningless local maximum; if one starts instead from the consistently updated model of diagrams 54 and 56, the optimization quickly converges to a best model with an improved fit to the seismic data. Once the best model is found, the posterior covariance matrix may be evaluated again for more accurate uncertainty quantification.
In the previous section we also described a Monte Carlo approach to quantify uncertainty. If the model uncertainty is quantified in a multivariate normal distribution whose mean ~..~, and covariance matrix C are evaluated from the result of Monte Carlo sampling, the consistent update method we have illustrated above can be used to compute a quick update, i.e.
without requiring the Monte Carlo sampling process to be repeated. We now illustrate a second method that uses additional information to directly update the sample from the posterior PDF obtained by a Monte Carlo method.
The method we describe can be applied when the additional information that becomes available directly controls the properties of the posterior PDF sample. In the example of Figure 5, diagrams 86-90, a Monte Carlo method was used to obtain a sample from the posterior PDF of the thicknesses and compressional wave velocities in a layered medium. Because the seismic reflection data are in travel time, obtaining a layered medium description in depth requires having a travel time-depth tie, that is, having a point where the absolute depth and the corresponding travel time are known or measured. This travel time-depth tie may be the Earth's surface, where depth and travel time are zero. For offshore reservoirs, the travel time-depth tie may be at the sea bottom. If the seismic data were acquired using receivers or sources in a borehole, the travel time-depth tie may be obtained at the depth of a receiver/source in the borehole.
As will be clear to those skilled in the art, the uncertainty in the depth of a layer as determined from seismic reflection data is greater the farther. away the layer is from the deepest travel time-depth tie. This is because depth is the sum of layer thicknesses starting from the deepest travel time-depth tie; the uncertainty in depth is the cumulative uncertainty of the layer thicknesses.
Diagrams 92 and 94 illustrate our method of updating a layered model and its uncertainties as new information on travel time-depth ties becomes available. Diagram 92 shows the results of Monte Carlo sampling obtained for a relatively shallow travel time-depth tie (e.g. the Earth's surface, the sea floor, the depth of a seismic receiver/source in a borehole). The results are presented in the form of a posterior image (as in diagram 90) computed by superimposing all the layered media sampled by the Monte Carlo method. As noted in the previous section, in practice the Monte Carlo sampling obtains a sample of layered media in travel time, rather than depth; each of these layered media is then converted to depth using a known/measured travel time-depth tie.
Diagram 94 illustrates how the uncertainty in depth changes if a deeper travel time-depth tie is measured. This additional travel time-depth tie will typically be provided by a seismic receiver/source placed in a borehole, and it could be acquired during or after drilling the well. Diagram 94 shows how the uncertainty in deep layers is reduced by using the deeper travel time-depth tie. It should be stressed that the application of this method does not require sampling the posterior PDF again, but simply using the same sample of layered media in travel time obtained earlier while recomputing the travel time to depth conversion for the additional, deeper travel time-depth tie. This reduction in uncertainty can be extremely important in making decisions during the drilling of a well (e. g., on determining the depth of zones that have anomalous pore pressures and are drilling hazards).
Uncertainty quantification and consistent model updating can improve significantly the efficiency of constructing and modifying a SEM. Because industry-standard interpretation workflows don't account for uncertainty, model consistency is maintained by making elements of the model 'rigid' as the interpretation progresses down the workflow. The term 'rigid' means here that once a domain expert has set the optimal values for parameters under his control (e. g., a geophysicist interpreting the model framework), these values are not changed by later experts (e.g., a reservoir engineer) for fear that the model will no longer be consistent with the previous data. For example, once the model framework is fixed, the reservoir engineer is only left flow parameters to adjust during history matching. Explicitly accounting for uncertainty would allow that same reservoir engineer to adjust all parameters of the model within their acceptable ranges to obtain the best history match. The remaining uncertainty in the model then provides input to decision-making tools which, for example, can be used to estimate a PDF of the net present value of a reservoir. Thus, uncertainty quantification allows us to approach the ideal of a SEM that is constrained by as many data as possible, can be easily updated, and ties directly into decision-making tools.
The foregoing description of the preferred and alternate embodiments of the present invention have been presented for purposes of illustration and description. It is not intended to be exhaustive or limit the invention to the precise examples described. Many modifications and variations will be apparent to those skilled in the art. The embodiments were chosen and described in order to best explain the principles of the invention and its practical application, thereby enabling others skilled in the art to understand the invention for various embodiments and with various modifications as are suited to the particular use contemplated. It is intended that the scope of the invention be defined by the accompanying claims and their equivalents.
Claims (21)
1. A method of using measurement data to create a model of a subsurface area comprising the steps of:
creating a parameterized model of the subsurface area having an initial estimate of model parameter uncertainties;
considering measurement data from the subsurface area;
updating said parameterized model to fit said measurement data, said updated parameterized model having an updated estimate of model parameter uncertainties; and repeating said considering and updating steps with additional measurement data to produce a further updated parameterized model having a further updated estimate of model parameter uncertainties.
creating a parameterized model of the subsurface area having an initial estimate of model parameter uncertainties;
considering measurement data from the subsurface area;
updating said parameterized model to fit said measurement data, said updated parameterized model having an updated estimate of model parameter uncertainties; and repeating said considering and updating steps with additional measurement data to produce a further updated parameterized model having a further updated estimate of model parameter uncertainties.
2. A method according to claim 1, wherein said model of the subsurface area has geometrical model parameters representing geological boundaries.
3. A method according to claim 1, wherein said model of the subsurface area has material parameters representing properties of distributed subsurface materials.
4. A method according to claim 1, wherein said measurement data comprise seismic data, drilling data, well logging data, well test data, production history data, or permanent monitoring data.
5. A method according to claim 1, wherein said measurement data and said additional measurement data comprise different types of data.
6. A method according to claim 5, wherein said measurement data comprises seismic data and said additional measurement data comprises well logging data.
7. A method according to claim 1, wherein said parameterized model of the subsurface area comprises a layered model having material properties associated with each layer.
8. A method according to claim 1, wherein said parameterized model of the subsurface area comprises a geocellular model having material property values defined on a regular or irregular three-dimensional grid.
9. A method according to claim 1, wherein said parameterized model of the subsurface area comprises a geometry-based model having material property values defined on a plurality of discrete geometrical sub-regions within the subsurface area.
10. A method according to claim 1, wherein said initial and updated estimates of model parameter uncertainties comprise probability density functions computed using deterministic optimization or Monte Carlo sampling.
11. A method according to claim 10, wherein said probability density functions are multivariate normal probability density functions definable by mean vectors and covariance matrices.
12. A method according to claim 11, wherein a mean and covariance matrix of said updated estimate of model parameter uncertainties is computed using deterministic optimization or Monte Carlo sampling.
13. A method according to claim 1, wherein said considering step comprises comparing said measurement data to a model parameter, said updating step comprises updating said parameterized model to fit said measurement data, and said updating step further comprising updating said initial estimate of model parameter uncertainties to be consistent with uncertainties associated with said measurement data.
14. A method according to claim 1, wherein said considering step comprising comparing said additional measurement data to predicted data produced by a prediction algorithm operating on said initial parameterized model.
15. A method according to claim 1, wherein said considering step produces a likelihood function that is combined with said initial estimate of model parameter uncertainties to produce said updated estimate of model parameter uncertainties.
16. A method according to claim 1, wherein said updating step comprises directly updating a sample from a posterior PDF obtained by a Monte Carlo method using said measurement data.
17. A method according to claim 17, wherein said updating step comprises directly updating said sample using travel time-depth measurement data obtained while drilling.
18. A method according to claim 1, wherein said parameterized model parameters may be changed during said updating step only when said change is consistent with said initial estimate of model parameter uncertainties.
19. A method according to claim 1, wherein any changes to said parameterized model parameters that are inconsistent with said initial estimate of model parameter uncertainties are identified during said updating step.
20. An apparatus for creating a model of a subsurface area comprising:
means for creating a parameterized model of the subsurface area having an initial estimate of model parameter uncertainties;
means for considering measurement data from the subsurface area;
means for updating said parameterized model to fit said measurement data, said updated parameterized model having an updated estimate of model parameter uncertainties;
and means for repeating said considering and updating steps with additional measurement data to produce a further updated parameterized model having a further updated estimate of model parameter uncertainties.
means for creating a parameterized model of the subsurface area having an initial estimate of model parameter uncertainties;
means for considering measurement data from the subsurface area;
means for updating said parameterized model to fit said measurement data, said updated parameterized model having an updated estimate of model parameter uncertainties;
and means for repeating said considering and updating steps with additional measurement data to produce a further updated parameterized model having a further updated estimate of model parameter uncertainties.
21. An article of manufacture, comprising:
a computer usable medium having a computer readable program code means embodied therein for creating a model of a subsurface area, the computer readable program code means in said article of manufacture comprising:
computer-readable program means for creating a parameterized model of the subsurface area having an initial estimate of model parameter uncertainties;
computer-readable program means for considering measurement data from the subsurface area;
computer-readable program means for updating said parameterized model to fit said measurement data, said updated parameterized model having an updated estimate of model parameter uncertainties; and computer-readable program means for repeating said considering and updating steps with additional measurement data to produce a further updated parameterized model having a further updated estimate of model parameter uncertainties.
a computer usable medium having a computer readable program code means embodied therein for creating a model of a subsurface area, the computer readable program code means in said article of manufacture comprising:
computer-readable program means for creating a parameterized model of the subsurface area having an initial estimate of model parameter uncertainties;
computer-readable program means for considering measurement data from the subsurface area;
computer-readable program means for updating said parameterized model to fit said measurement data, said updated parameterized model having an updated estimate of model parameter uncertainties; and computer-readable program means for repeating said considering and updating steps with additional measurement data to produce a further updated parameterized model having a further updated estimate of model parameter uncertainties.
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11982199P | 1999-02-12 | 1999-02-12 | |
US60/119,821 | 1999-02-12 | ||
US16533399P | 1999-11-12 | 1999-11-12 | |
US60/165,333 | 1999-11-12 | ||
PCT/US2000/003615 WO2000048022A1 (en) | 1999-02-12 | 2000-02-11 | Uncertainty constrained subsurface modeling |
Publications (2)
Publication Number | Publication Date |
---|---|
CA2362285A1 CA2362285A1 (en) | 2000-08-17 |
CA2362285C true CA2362285C (en) | 2005-06-14 |
Family
ID=26817738
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CA002362285A Expired - Lifetime CA2362285C (en) | 1999-02-12 | 2000-02-11 | Uncertainty constrained subsurface modeling |
Country Status (6)
Country | Link |
---|---|
US (1) | US6549854B1 (en) |
EP (1) | EP1151326B1 (en) |
AU (1) | AU3229900A (en) |
CA (1) | CA2362285C (en) |
NO (1) | NO20013894L (en) |
WO (1) | WO2000048022A1 (en) |
Families Citing this family (200)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2354852B (en) * | 1999-10-01 | 2001-11-28 | Schlumberger Holdings | Method for updating an earth model using measurements gathered during borehole construction |
US6826486B1 (en) | 2000-02-11 | 2004-11-30 | Schlumberger Technology Corporation | Methods and apparatus for predicting pore and fracture pressures of a subsurface formation |
AU2001247576A1 (en) * | 2000-03-17 | 2001-10-03 | Optim, L.L.C. | Optimization apparatus, system, and method of use doing business |
US7035790B2 (en) * | 2000-06-02 | 2006-04-25 | Canon Kabushiki Kaisha | Speech processing system |
US20020026253A1 (en) * | 2000-06-02 | 2002-02-28 | Rajan Jebu Jacob | Speech processing apparatus |
US7010483B2 (en) * | 2000-06-02 | 2006-03-07 | Canon Kabushiki Kaisha | Speech processing system |
US6954745B2 (en) * | 2000-06-02 | 2005-10-11 | Canon Kabushiki Kaisha | Signal processing system |
US7072833B2 (en) * | 2000-06-02 | 2006-07-04 | Canon Kabushiki Kaisha | Speech processing system |
US7003439B2 (en) * | 2001-01-30 | 2006-02-21 | Schlumberger Technology Corporation | Interactive method for real-time displaying, querying and forecasting drilling event and hazard information |
US7526413B2 (en) * | 2001-01-31 | 2009-04-28 | Exxonmobil Upstream Research Company | Volumetric laminated sand analysis |
FR2823877B1 (en) | 2001-04-19 | 2004-12-24 | Inst Francais Du Petrole | METHOD FOR CONSTRAINING BY DYNAMIC PRODUCTION DATA A FINE MODEL REPRESENTATIVE OF THE DISTRIBUTION IN THE DEPOSIT OF A CHARACTERISTIC PHYSICAL SIZE OF THE BASEMENT STRUCTURE |
US6581198B1 (en) * | 2001-06-13 | 2003-06-17 | Cadence Design Systems, Inc. | Method and arrangement for extracting capacitance in integrated circuits having non Manhattan wiring |
US6892366B1 (en) * | 2001-08-28 | 2005-05-10 | Cadence Design Systems, Inc. | Method and apparatus for performing extraction using a model trained with Bayesian inference via a Monte Carlo method |
US7051293B1 (en) | 2001-08-28 | 2006-05-23 | Cadence Design Systems, Inc. | Method and apparatus for creating an extraction model |
US7103524B1 (en) | 2001-08-28 | 2006-09-05 | Cadence Design Systems, Inc. | Method and apparatus for creating an extraction model using Bayesian inference implemented with the Hybrid Monte Carlo method |
GB0125713D0 (en) | 2001-10-26 | 2001-12-19 | Statoil Asa | Method of combining spatial models |
US6880138B1 (en) | 2002-01-31 | 2005-04-12 | Cadence Design Systems, Inc. | Method and apparatus for creating a critical input space spanning set of input points to train a machine learning model for extraction |
FR2837947B1 (en) * | 2002-04-02 | 2004-05-28 | Inst Francais Du Petrole | METHOD FOR QUANTIFYING THE UNCERTAINTIES RELATED TO CONTINUOUS AND DESCRIPTIVE PARAMETERS OF A MEDIUM BY CONSTRUCTION OF EXPERIMENT PLANS AND STATISTICAL ANALYSIS |
WO2004022905A1 (en) * | 2002-04-10 | 2004-03-18 | Schlumberger Technology Corporation | Method, apparatus and system for pore pressure prediction in presence of dipping formations |
US7899657B2 (en) * | 2003-01-24 | 2011-03-01 | Rockwell Automoation Technologies, Inc. | Modeling in-situ reservoirs with derivative constraints |
US6810332B2 (en) * | 2003-01-31 | 2004-10-26 | Chevron U.S.A. Inc. | Method for computing complexity, confidence and technical maturity indices for reservoir evaluations |
MXPA04003024A (en) * | 2003-04-03 | 2004-10-06 | Schlumberger Technology Bv | Methods and apparatus for predicting pore and fracture pressures of a subsurface formation. |
US7835893B2 (en) * | 2003-04-30 | 2010-11-16 | Landmark Graphics Corporation | Method and system for scenario and case decision management |
US20050004969A1 (en) * | 2003-07-03 | 2005-01-06 | Reaction Design, Llc | System for information exchange for integration of multiple data sources |
JP2007531068A (en) * | 2003-07-03 | 2007-11-01 | リアクション デザイン | Method and system for integrated uncertainty analysis |
US20050004833A1 (en) * | 2003-07-03 | 2005-01-06 | Reaction Design, Llc | Method and system for integrated uncertainty analysis |
US7286214B2 (en) * | 2003-08-25 | 2007-10-23 | University Of South Florida | Method and program product for determining a radiance field in an optical environment |
NL1024444C2 (en) * | 2003-10-03 | 2005-04-08 | J O A Beheer B V | Method, device, computer program and data carrier for modeling a multi-dimensional heterogeneous structure with a digital processing unit. |
US7379819B2 (en) * | 2003-12-04 | 2008-05-27 | Schlumberger Technology Corporation | Reservoir sample chain-of-custody |
GB2410550B8 (en) * | 2003-12-04 | 2008-10-01 | Schlumberger Holdings | Fluids chain-of-custody |
US7359844B2 (en) | 2004-01-20 | 2008-04-15 | Saudi Arabian Oil Company | Real time earth model for collaborative geosteering |
EP1714230B1 (en) | 2004-01-30 | 2010-04-14 | ExxonMobil Upstream Research Company | Reservoir model building methods |
CN1898640A (en) * | 2004-01-30 | 2007-01-17 | 埃克森美孚上游研究公司 | Reservoir evaluation methods |
US20050171700A1 (en) * | 2004-01-30 | 2005-08-04 | Chroma Energy, Inc. | Device and system for calculating 3D seismic classification features and process for geoprospecting material seams |
US7054750B2 (en) * | 2004-03-04 | 2006-05-30 | Halliburton Energy Services, Inc. | Method and system to model, measure, recalibrate, and optimize control of the drilling of a borehole |
CN1299127C (en) * | 2004-03-26 | 2007-02-07 | 中国石油天然气集团公司 | Laminated dielectric double-focusing method for earthquake observation system optimized design |
AU2005272112B2 (en) * | 2004-07-07 | 2010-04-22 | Exxonmobil Upstream Research Company | Predicting sand-grain composition and sand texture |
CA2572981A1 (en) * | 2004-07-07 | 2006-10-26 | Exxonmobil Upstream Research Company | Bayesian network applications to geology and geophysics |
US20060025969A1 (en) * | 2004-07-30 | 2006-02-02 | Conocophillips Company | Method for converting a geologic panel into a simulator format |
FR2874706B1 (en) * | 2004-08-30 | 2006-12-01 | Inst Francais Du Petrole | METHOD OF MODELING THE PRODUCTION OF A PETROLEUM DEPOSITION |
CA2586603C (en) * | 2004-11-18 | 2014-04-01 | Exxonmobil Upstream Research Company | Method for quantification and mitigation for dip-induced azimuthal avo |
US8209202B2 (en) * | 2005-04-29 | 2012-06-26 | Landmark Graphics Corporation | Analysis of multiple assets in view of uncertainties |
FR2886740B1 (en) * | 2005-06-03 | 2007-09-28 | Inst Francais Du Petrole | METHOD FOR UPDATING A GEOLOGICAL MODEL WITH SEISMIC AND PRODUCTION DATA |
EP1922663A4 (en) | 2005-07-27 | 2015-11-04 | Exxonmobil Upstream Res Co | Well modeling associated with extraction of hydrocarbons from subsurface formations |
MX2007016586A (en) * | 2005-07-27 | 2008-03-04 | Exxonmobil Upstream Res Co | Well modeling associated with extraction of hydrocarbons from subsurface formations. |
EA200800436A1 (en) * | 2005-07-27 | 2008-08-29 | Эксонмобил Апстрим Рисерч Компани | MODELING OF A WELL, ASSOCIATED WITH THE PRODUCTION OF HYDROCARBONS FROM UNDERGROUND FORMATIONS |
GB0524134D0 (en) * | 2005-11-26 | 2006-01-04 | Univ Edinburgh | Improvements in and relating to hydrocarbon recovery from a hydrocarbon reservoir |
US7366616B2 (en) * | 2006-01-13 | 2008-04-29 | Schlumberger Technology Corporation | Computer-based method for while-drilling modeling and visualization of layered subterranean earth formations |
US7486589B2 (en) | 2006-02-09 | 2009-02-03 | Schlumberger Technology Corporation | Methods and apparatus for predicting the hydrocarbon production of a well location |
CA2664352C (en) * | 2006-09-28 | 2011-09-27 | Exxonmobil Upstream Research Company | Iterative inversion of data from simultaneous geophysical sources |
AU2007339997A1 (en) * | 2006-12-28 | 2008-07-10 | Chevron U.S.A. Inc. | History matching and forecasting of hydrocarbon-bearing reservoirs utilizing proxies for likelihood functions |
US7467044B2 (en) * | 2007-01-15 | 2008-12-16 | Chevron U.S.A. Inc | Method and system for assessing exploration prospect risk and uncertainty |
US7636275B2 (en) * | 2007-02-06 | 2009-12-22 | Conocophillips Company | Direct time lapse inversion of seismic data |
US7529624B2 (en) * | 2007-02-21 | 2009-05-05 | Geomechanics International, Inc. | Method and apparatus for remote characterization of faults in the vicinity of boreholes |
US8548782B2 (en) | 2007-08-24 | 2013-10-01 | Exxonmobil Upstream Research Company | Method for modeling deformation in subsurface strata |
CA2690991C (en) * | 2007-08-24 | 2013-12-24 | Exxonmobil Upstream Research Company | Method for multi-scale geomechanical model analysis by computer simulation |
US8768672B2 (en) * | 2007-08-24 | 2014-07-01 | ExxonMobil. Upstream Research Company | Method for predicting time-lapse seismic timeshifts by computer simulation |
CA2690992C (en) * | 2007-08-24 | 2014-07-29 | Exxonmobil Upstream Research Company | Method for predicting well reliability by computer simulation |
MX2010005116A (en) * | 2007-11-10 | 2010-09-09 | Landmark Graphics Corp | Systems and methods for workflow automation, adaptation and integration. |
US7886844B2 (en) * | 2007-11-12 | 2011-02-15 | Schlumberger Technology Corporation | Borehole survey method and apparatus |
GB0722469D0 (en) * | 2007-11-16 | 2007-12-27 | Statoil Asa | Forming a geological model |
FR2923855B1 (en) * | 2007-11-21 | 2009-12-11 | Inst Francais Du Petrole | METHOD FOR OPTIMIZING THE TRACK OF A DRILL BY UPDATING STRUCTURAL INTERPRETATION |
WO2009075961A1 (en) * | 2007-12-12 | 2009-06-18 | Exxonmobil Upstream Research Company | Method and apparatus for evaluating submarine formations |
US9638830B2 (en) | 2007-12-14 | 2017-05-02 | Westerngeco L.L.C. | Optimizing drilling operations using petrotechnical data |
US8417497B2 (en) * | 2008-01-18 | 2013-04-09 | Westerngeco L.L.C. | Updating a model of a subterranean structure using decomposition |
US7869955B2 (en) | 2008-01-30 | 2011-01-11 | Chevron U.S.A. Inc. | Subsurface prediction method and system |
BRPI0908156A2 (en) * | 2008-03-21 | 2015-08-11 | Exxonmobil Upstream Company | Methods for invoking measured geophysical data and for producing hydrocarbons from a subsurface region |
US8793111B2 (en) * | 2009-01-20 | 2014-07-29 | Schlumberger Technology Corporation | Automated field development planning |
US8527248B2 (en) * | 2008-04-18 | 2013-09-03 | Westerngeco L.L.C. | System and method for performing an adaptive drilling operation |
WO2009137181A1 (en) * | 2008-05-05 | 2009-11-12 | Exxonmobil Upstream Research Company | Modeling dynamic systems by visualizing and narrowing a parameter space |
US8527203B2 (en) | 2008-05-27 | 2013-09-03 | Schlumberger Technology Corporation | Method for selecting well measurements |
GB2474157B (en) * | 2008-07-01 | 2012-10-17 | Schlumberger Holdings | Effective hydrocarbon reservoir exploration decision making |
US8892410B2 (en) * | 2008-08-11 | 2014-11-18 | Exxonmobil Upstream Research Company | Estimation of soil properties using waveforms of seismic surface waves |
BRPI0914510B1 (en) | 2008-10-14 | 2019-07-16 | Prad Research And Development Limited | METHODS TO AUTOMATE OR PARTIALLY AUTOMATE OPTIMIZATION OF AN AUTOMATED OR PARTIALLY AUTOMATED DRILLING OPERATION, SYSTEM TO AUTOMATE AUTOMATICALLY TO PRODUCT PERFORMANCE, AUTOMATED PROCESSING PERFORMANCE DRILLING CONTROL |
CA2737415C (en) * | 2008-11-06 | 2017-03-28 | Exxonmobil Upstream Research Company | System and method for planning a drilling operation |
WO2010071701A1 (en) * | 2008-12-16 | 2010-06-24 | Exxonmobil Upstream Research Company | Systems and methods for hydrocarbon reservoir development and management optimization |
CN102257407B (en) * | 2008-12-18 | 2014-05-21 | 埃克森美孚上游研究公司 | Sediment transport by fully developed turbulent flows |
BRPI1006862B1 (en) | 2009-01-13 | 2020-03-17 | Exxonmobil Upstream Research Company | METHOD AND SYSTEM FOR OPTIMIZING DECISION-MAKING FOR A HYDROCARBONET WELL, AND, SYSTEM ASSOCIATED WITH HYDROCARBON PRODUCTION |
GB0905261D0 (en) * | 2009-03-27 | 2009-05-13 | Geco Technology Bv | Processing seismic data |
US20100257004A1 (en) * | 2009-04-01 | 2010-10-07 | Chervon U.S.A. Inc. | Method and system for conducting geologic basin analysis |
US8341984B2 (en) * | 2009-05-27 | 2013-01-01 | Westerngeco L.L.C. | Estimating velocities with uncertainty |
US8494827B2 (en) * | 2009-09-25 | 2013-07-23 | Exxonmobil Upstream Research Company | Method of predicting natural fractures and damage in a subsurface region |
NO20101371A1 (en) * | 2009-10-05 | 2011-04-06 | Logined Bv | Method, system and apparatus for modeling uncertainty in production system networks |
US8234073B2 (en) * | 2009-10-23 | 2012-07-31 | Chevron U.S.A. Inc. | System and method for estimating geological architecture of a geologic volume |
US20110098996A1 (en) * | 2009-10-26 | 2011-04-28 | David Nichols | Sifting Models of a Subsurface Structure |
WO2011062794A1 (en) * | 2009-11-18 | 2011-05-26 | Conocophillips Company | Attribute importance measure for parametric multivariate modeling |
BR112012012817A2 (en) | 2009-11-30 | 2020-08-11 | Exxonmobil Upstream Research Company | adaptive newton method for reservoir simulation |
AU2011213261B2 (en) * | 2010-02-02 | 2015-04-02 | Conocophillips Company | Multilevel percolation aggregation solver for petroleum reservoir simulations |
EP2531694B1 (en) | 2010-02-03 | 2018-06-06 | Exxonmobil Upstream Research Company | Method for using dynamic target region for well path/drill center optimization |
US8537638B2 (en) * | 2010-02-10 | 2013-09-17 | Exxonmobil Upstream Research Company | Methods for subsurface parameter estimation in full wavefield inversion and reverse-time migration |
WO2011100002A1 (en) | 2010-02-12 | 2011-08-18 | Exxonmobil Upstream Research Company | Method and system for partitioning parallel simulation models |
US9367564B2 (en) | 2010-03-12 | 2016-06-14 | Exxonmobil Upstream Research Company | Dynamic grouping of domain objects via smart groups |
US8223587B2 (en) * | 2010-03-29 | 2012-07-17 | Exxonmobil Upstream Research Company | Full wavefield inversion using time varying filters |
CA2795835C (en) | 2010-04-30 | 2016-10-04 | Exxonmobil Upstream Research Company | Method and system for finite volume simulation of flow |
US8694299B2 (en) | 2010-05-07 | 2014-04-08 | Exxonmobil Upstream Research Company | Artifact reduction in iterative inversion of geophysical data |
US8756042B2 (en) | 2010-05-19 | 2014-06-17 | Exxonmobile Upstream Research Company | Method and system for checkpointing during simulations |
US10087721B2 (en) | 2010-07-29 | 2018-10-02 | Exxonmobil Upstream Research Company | Methods and systems for machine—learning based simulation of flow |
WO2012015521A1 (en) | 2010-07-29 | 2012-02-02 | Exxonmobil Upstream Research Company | Method and system for reservoir modeling |
CA2805446C (en) | 2010-07-29 | 2016-08-16 | Exxonmobil Upstream Research Company | Methods and systems for machine-learning based simulation of flow |
US8767508B2 (en) | 2010-08-18 | 2014-07-01 | Exxonmobil Upstream Research Company | Using seismic P and S arrivals to determine shallow velocity structure |
AU2011295892A1 (en) | 2010-09-03 | 2013-05-02 | Chevron U.S.A. Inc. | Iterative method and system to construct robust proxy models for reservoir simulation |
CA2807300C (en) | 2010-09-20 | 2017-01-03 | Exxonmobil Upstream Research Company | Flexible and adaptive formulations for complex reservoir simulations |
US8437998B2 (en) | 2010-09-27 | 2013-05-07 | Exxonmobil Upstream Research Company | Hybrid method for full waveform inversion using simultaneous and sequential source method |
EP2622457A4 (en) * | 2010-09-27 | 2018-02-21 | Exxonmobil Upstream Research Company | Simultaneous source encoding and source separation as a practical solution for full wavefield inversion |
CA2814669A1 (en) | 2010-11-23 | 2012-05-31 | Exxonmobil Upstream Research Company | Variable discretization method for flow simulation on complex geological models |
CN103238158B (en) | 2010-12-01 | 2016-08-17 | 埃克森美孚上游研究公司 | Utilize the marine streamer data source inverting simultaneously that mutually related objects function is carried out |
AU2010365379B2 (en) * | 2010-12-16 | 2015-07-23 | Landmark Graphics Corporation | Method and system of plotting correlated data |
GB201022128D0 (en) | 2010-12-31 | 2011-02-02 | Foster Findlay Ass Ltd | Active contour segmentation |
CN103348265B (en) * | 2011-01-13 | 2016-10-19 | 界标制图有限公司 | Update the method and system of geologic element model |
WO2012102784A1 (en) | 2011-01-26 | 2012-08-02 | Exxonmobil Upstream Research Company | Method of reservoir compartment analysis using topological structure in 3d earth model |
US9146330B2 (en) * | 2011-03-29 | 2015-09-29 | Westerngeco L.L.C. | Selecting a survey setting for characterizing a target structure |
CN103703391B (en) | 2011-03-30 | 2017-05-17 | 埃克森美孚上游研究公司 | System of full wavefield inversion using spectral shaping and computer implementing method |
CN103460074B (en) | 2011-03-31 | 2016-09-28 | 埃克森美孚上游研究公司 | Wavelet estimators and the method for many subwaves prediction in full wave field inversion |
US8749549B2 (en) * | 2011-06-10 | 2014-06-10 | Schlumberger Technology Corporation | Prospect assessment and play chance mapping tools |
US8902221B2 (en) * | 2011-06-10 | 2014-12-02 | Schlumberger Technology Corporation | Prospect assessment and play chance mapping tools |
US9846256B2 (en) | 2011-08-09 | 2017-12-19 | Schlumberger Technology Corporation | Interactive display of results obtained from the inversion of logging data |
US8694262B2 (en) * | 2011-08-15 | 2014-04-08 | Chevron U.S.A. Inc. | System and method for subsurface characterization including uncertainty estimation |
ES2640824T3 (en) | 2011-09-02 | 2017-11-06 | Exxonmobil Upstream Research Company | Use of projection on convex assemblies to limit the inversion of the entire wave field |
EP2756382A4 (en) | 2011-09-15 | 2015-07-29 | Exxonmobil Upstream Res Co | Optimized matrix and vector operations in instruction limited algorithms that perform eos calculations |
BR112014009734A2 (en) * | 2011-10-20 | 2017-04-18 | Prad Res & Dev Ltd | method for controlling hydrocarbon recovery equipment from a reservoir |
US20130223187A1 (en) * | 2011-11-11 | 2013-08-29 | International Geophysical Company, Inc. | Geological Structure Contour Modeling and Imaging |
US9176930B2 (en) * | 2011-11-29 | 2015-11-03 | Exxonmobil Upstream Research Company | Methods for approximating hessian times vector operation in full wavefield inversion |
US10030499B2 (en) | 2011-12-06 | 2018-07-24 | Bp Corporation North America Inc. | Geological monitoring console |
CA2861863A1 (en) | 2012-03-08 | 2013-09-12 | Exxonmobil Upstream Research Company | Orthogonal source and receiver encoding |
EP2831647A4 (en) * | 2012-03-28 | 2016-02-24 | Geco Technology Bv | Providing an objective function based on variation in predicted data |
US20130261981A1 (en) * | 2012-04-03 | 2013-10-03 | Westerngeco L.L.C. | Covariance estimation using sparse wavelet representation |
WO2013169429A1 (en) | 2012-05-08 | 2013-11-14 | Exxonmobile Upstream Research Company | Canvas control for 3d data volume processing |
NO342738B1 (en) * | 2012-12-06 | 2018-08-06 | Roxar Software Solutions As | Procedure and a system for presenting seismic information |
US20130341093A1 (en) * | 2012-06-21 | 2013-12-26 | Stuart Inglis Jardine | Drilling risk avoidance |
EP2877690B1 (en) | 2012-08-10 | 2023-03-01 | Halliburton Energy Services Inc. | Methods and systems for borehole image forward modeling of formation properties |
AU2013324162B2 (en) | 2012-09-28 | 2018-08-09 | Exxonmobil Upstream Research Company | Fault removal in geological models |
CN102937724B (en) * | 2012-11-20 | 2016-06-08 | 中国神华能源股份有限公司 | A kind of detection method of stope bottom rock stratum of open coal mine |
WO2014084945A1 (en) | 2012-11-28 | 2014-06-05 | Exxonmobil Upstream Resarch Company | Reflection seismic data q tomography |
NO20121473A1 (en) * | 2012-12-06 | 2014-06-09 | Roxar Software Solutions As | System for modeling geological structures |
US9575205B2 (en) * | 2013-01-17 | 2017-02-21 | Pgs Geophysical As | Uncertainty-based frequency-selected inversion of electromagnetic geophysical data |
BR112015021666B1 (en) * | 2013-03-05 | 2021-01-26 | Technological Resources Pty. Limited | computer-implemented method for updating an estimate for a material property of a volume and computer system for updating an estimate for a material property of a volume |
US20140257700A1 (en) * | 2013-03-08 | 2014-09-11 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | System and method for estimating uncertainty for geophysical gridding routines lacking inherent uncertainty estimation |
CN105308481B (en) * | 2013-04-22 | 2019-08-02 | 斯塔特伊石油公司 | The method for determining temperature under earth's surface |
CA2909105C (en) | 2013-05-24 | 2018-08-28 | Ke Wang | Multi-parameter inversion through offset dependent elastic fwi |
US10459117B2 (en) | 2013-06-03 | 2019-10-29 | Exxonmobil Upstream Research Company | Extended subspace method for cross-talk mitigation in multi-parameter inversion |
MX350181B (en) | 2013-06-03 | 2017-08-29 | Exxonmobil Upstream Res Co | Uncertainty estimation of subsurface resistivity solutions. |
CA2907728C (en) | 2013-06-10 | 2021-04-27 | Exxonmobil Upstream Research Company | Interactively planning a well site |
US9702998B2 (en) | 2013-07-08 | 2017-07-11 | Exxonmobil Upstream Research Company | Full-wavefield inversion of primaries and multiples in marine environment |
DK3036566T3 (en) | 2013-08-23 | 2018-07-23 | Exxonmobil Upstream Res Co | SIMILAR SOURCE APPLICATION DURING BOTH SEISMIC COLLECTION AND SEISMIC INVERSION |
SG11201600091WA (en) * | 2013-08-23 | 2016-02-26 | Landmark Graphics Corp | Local updating of 3d geocellular model |
US10036818B2 (en) | 2013-09-06 | 2018-07-31 | Exxonmobil Upstream Research Company | Accelerating full wavefield inversion with nonstationary point-spread functions |
US9864098B2 (en) | 2013-09-30 | 2018-01-09 | Exxonmobil Upstream Research Company | Method and system of interactive drill center and well planning evaluation and optimization |
US10400570B2 (en) * | 2013-11-13 | 2019-09-03 | Schlumberger Technology Corporation | Automatic wellbore condition indicator and manager |
US9910189B2 (en) | 2014-04-09 | 2018-03-06 | Exxonmobil Upstream Research Company | Method for fast line search in frequency domain FWI |
AU2015256626B2 (en) | 2014-05-09 | 2017-10-12 | Exxonmobil Upstream Research Company | Efficient line search methods for multi-parameter full wavefield inversion |
US10185046B2 (en) | 2014-06-09 | 2019-01-22 | Exxonmobil Upstream Research Company | Method for temporal dispersion correction for seismic simulation, RTM and FWI |
RU2016150545A (en) | 2014-06-17 | 2018-07-17 | Эксонмобил Апстрим Рисерч Компани | FAST VISCOACOUSTIC AND VISCOELASTIC INVERSION OF A FULL WAVE FIELD |
GB2528130A (en) | 2014-07-11 | 2016-01-13 | Total E&P Uk Ltd | Method of constraining an inversion in the characterisation of the evolution of a subsurface volume |
GB2528129A (en) | 2014-07-11 | 2016-01-13 | Total E&P Uk Ltd | Method for obtaining estimates of a model parameter so as to characterise the evolution of a subsurface volume |
US10838092B2 (en) | 2014-07-24 | 2020-11-17 | Exxonmobil Upstream Research Company | Estimating multiple subsurface parameters by cascaded inversion of wavefield components |
US10422899B2 (en) | 2014-07-30 | 2019-09-24 | Exxonmobil Upstream Research Company | Harmonic encoding for FWI |
WO2016018723A1 (en) | 2014-07-30 | 2016-02-04 | Exxonmobil Upstream Research Company | Method for volumetric grid generation in a domain with heterogeneous material properties |
US10359523B2 (en) | 2014-08-05 | 2019-07-23 | Exxonmobil Upstream Research Company | Exploration and extraction method and system for hydrocarbons |
US10386511B2 (en) | 2014-10-03 | 2019-08-20 | Exxonmobil Upstream Research Company | Seismic survey design using full wavefield inversion |
AU2015337108B2 (en) | 2014-10-20 | 2018-03-01 | Exxonmobil Upstream Research Company | Velocity tomography using property scans |
US10605939B2 (en) | 2014-10-27 | 2020-03-31 | Cgg Services Sas | Multi-vintage energy mapping |
US11598195B2 (en) | 2014-10-27 | 2023-03-07 | Baker Hughes, A Ge Company, Llc | Statistical approach to incorporate uncertainties of parameters in simulation results and stability analysis for earth drilling |
AU2015339883B2 (en) | 2014-10-31 | 2018-03-29 | Exxonmobil Upstream Research Company | Methods to handle discontinuity in constructing design space for faulted subsurface model using moving least squares |
WO2016069171A1 (en) | 2014-10-31 | 2016-05-06 | Exxonmobil Upstream Research Company | Handling domain discontinuity in a subsurface grid model with the help of grid optimization techniques |
CA2961923A1 (en) | 2014-10-31 | 2016-05-06 | Exxonmobil Upstream Research Company | Managing discontinuities in geologic models |
EP3234659A1 (en) | 2014-12-18 | 2017-10-25 | Exxonmobil Upstream Research Company | Scalable scheduling of parallel iterative seismic jobs |
US10190406B2 (en) * | 2014-12-23 | 2019-01-29 | Baker Hughes, A Ge Company, Llc | Formation fracturing potential using surrounding pore pressures |
WO2016110660A1 (en) | 2015-01-06 | 2016-07-14 | Total E&P Uk Limited | Method for obtaining estimates of a model parameter so as to characterise the evolution of a subsurface volume over a time period |
US10520618B2 (en) | 2015-02-04 | 2019-12-31 | ExxohnMobil Upstream Research Company | Poynting vector minimal reflection boundary conditions |
CA2972028C (en) | 2015-02-13 | 2019-08-13 | Exxonmobil Upstream Research Company | Efficient and stable absorbing boundary condition in finite-difference calculations |
RU2017132164A (en) | 2015-02-17 | 2019-03-18 | Эксонмобил Апстрим Рисерч Компани | MULTI-STAGE PROCESS OF INVERSION OF A FULL WAVE FIELD WHEN EXECUTING WHICH ARE ARRAYING AN ARRAY FREE OF MULTIPLE DATA WAVES |
SG11201708665VA (en) | 2015-06-04 | 2017-12-28 | Exxonmobil Upstream Res Co | Method for generating multiple free seismic images |
US10838093B2 (en) | 2015-07-02 | 2020-11-17 | Exxonmobil Upstream Research Company | Krylov-space-based quasi-newton preconditioner for full-wavefield inversion |
CN108139499B (en) | 2015-10-02 | 2020-02-14 | 埃克森美孚上游研究公司 | Q-compensated full wavefield inversion |
US20170103144A1 (en) * | 2015-10-08 | 2017-04-13 | Schlumbeger Technology Corporation | Well trajectory adjustment |
WO2017065889A1 (en) | 2015-10-15 | 2017-04-20 | Exxonmobil Upstream Research Company | Fwi model domain angle stacks with amplitude preservation |
FR3043227A1 (en) * | 2015-11-04 | 2017-05-05 | Services Petroliers Schlumberger | |
CN105653808A (en) * | 2016-01-05 | 2016-06-08 | 成都阳光铝制品有限公司 | Specified plasticity extension strength uncertainty evaluation method based on Monte Carlo |
US10768324B2 (en) | 2016-05-19 | 2020-09-08 | Exxonmobil Upstream Research Company | Method to predict pore pressure and seal integrity using full wavefield inversion |
US10605937B2 (en) * | 2016-05-26 | 2020-03-31 | Cgg Services Sas | Device and method for smart picking surface waves dispersion curves |
GB2556621B (en) * | 2016-09-30 | 2020-03-25 | Equinor Energy As | Improved structural modelling |
WO2018069744A1 (en) | 2016-10-14 | 2018-04-19 | Schlumberger Technology Corporation | Geologic structural model generation |
EP3559401B1 (en) | 2016-12-23 | 2023-10-18 | ExxonMobil Technology and Engineering Company | Method and system for stable and efficient reservoir simulation using stability proxies |
US11644589B2 (en) | 2017-05-23 | 2023-05-09 | Schlumberger Technology Corporation | Analogue facilitated seismic data interpretation system |
US11041976B2 (en) | 2017-05-30 | 2021-06-22 | Exxonmobil Upstream Research Company | Method and system for creating and using a subsurface model in hydrocarbon operations |
EP3645834B1 (en) | 2017-06-27 | 2024-04-10 | Services Pétroliers Schlumberger | Real-time well construction process inference through probabilistic data fusion |
CN110412530B (en) * | 2018-04-27 | 2021-09-17 | 深圳市优必选科技有限公司 | Method and device for identifying charging pile and robot |
US10557345B2 (en) | 2018-05-21 | 2020-02-11 | Saudi Arabian Oil Company | Systems and methods to predict and inhibit broken-out drilling-induced fractures in hydrocarbon wells |
US10753203B2 (en) | 2018-07-10 | 2020-08-25 | Saudi Arabian Oil Company | Systems and methods to identify and inhibit spider web borehole failure in hydrocarbon wells |
RU2718137C1 (en) * | 2018-12-05 | 2020-03-30 | Общество с ограниченной ответственностью "Газпромнефть Научно-Технический Центр" (ООО "Газпромнефть НТЦ") | Method of estimating wave field attenuation parameter for determining hydrocarbon saturation of formation in borehole space in constructing geological model |
WO2020163372A1 (en) | 2019-02-05 | 2020-08-13 | Motive Drilling Technologies, Inc. | Downhole display |
RU2720859C1 (en) * | 2019-05-29 | 2020-05-13 | Общество с ограниченной ответственностью "Газпром добыча Кузнецк" | Method of horizontal methane coal well design selection |
EP3990949A4 (en) * | 2019-06-28 | 2023-07-26 | Services Pétroliers Schlumberger | Field data acquisition and virtual training system |
CN110322563B (en) * | 2019-07-09 | 2020-11-10 | 吉林大学 | Supplementary drilling site selection method and system for building geological model |
US20220243575A1 (en) * | 2019-09-04 | 2022-08-04 | Schlumberger Technology Corporation | Autonomous operations in oil and gas fields |
US11348218B2 (en) | 2020-03-30 | 2022-05-31 | Halliburton Energy Services, Inc. | Hybrid interpretation approach for borehole imaging |
US11572785B2 (en) | 2021-01-26 | 2023-02-07 | Saudi Arabian Oil Company | Drilling uncertainty real time updates for accurate well placement |
EP4330732A1 (en) * | 2021-04-30 | 2024-03-06 | Zed Depth Exploration Data GmbH | Geophysical data-space mining using a trans-dimensional algorithm |
CN114137612A (en) * | 2021-11-30 | 2022-03-04 | 沈阳建筑大学 | Conditional mean spectrum generation method for vertical seismic oscillation considering spectral correlation |
WO2023214196A1 (en) * | 2022-05-02 | 2023-11-09 | Abu Dhabi National Oil Company | Method and system for exploiting a subterranean reservoir |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4340934A (en) * | 1971-09-07 | 1982-07-20 | Schlumberger Technology Corporation | Method of generating subsurface characteristic models |
US4679174A (en) * | 1984-04-26 | 1987-07-07 | Western Geophysical Company Of America | Method for seismic lithologic modeling |
US5103920A (en) * | 1989-03-01 | 1992-04-14 | Patton Consulting Inc. | Surveying system and method for locating target subterranean bodies |
US5321613A (en) * | 1992-11-12 | 1994-06-14 | Coleman Research Corporation | Data fusion workstation |
US5416750A (en) * | 1994-03-25 | 1995-05-16 | Western Atlas International, Inc. | Bayesian sequential indicator simulation of lithology from seismic data |
US5648937A (en) * | 1995-01-18 | 1997-07-15 | Atlantic Richfield Company | Method and apparatus for correlating geological structure horizons from velocity data to well observations |
US5539704A (en) | 1995-06-23 | 1996-07-23 | Western Atlas International, Inc. | Bayesian sequential Gaussian simulation of lithology with non-linear data |
US5838634A (en) * | 1996-04-04 | 1998-11-17 | Exxon Production Research Company | Method of generating 3-D geologic models incorporating geologic and geophysical constraints |
US6014343A (en) * | 1996-10-31 | 2000-01-11 | Geoquest | Automatic non-artificially extended fault surface based horizon modeling system |
US5905657A (en) | 1996-12-19 | 1999-05-18 | Schlumberger Technology Corporation | Performing geoscience interpretation with simulated data |
US6128577A (en) | 1996-12-19 | 2000-10-03 | Schlumberger Technology Corporation | Modeling geological structures and properties |
-
2000
- 2000-02-11 US US09/502,920 patent/US6549854B1/en not_active Expired - Lifetime
- 2000-02-11 CA CA002362285A patent/CA2362285C/en not_active Expired - Lifetime
- 2000-02-11 EP EP00910159A patent/EP1151326B1/en not_active Expired - Lifetime
- 2000-02-11 WO PCT/US2000/003615 patent/WO2000048022A1/en active IP Right Grant
- 2000-02-11 AU AU32299/00A patent/AU3229900A/en not_active Abandoned
-
2001
- 2001-08-09 NO NO20013894A patent/NO20013894L/en not_active Application Discontinuation
Also Published As
Publication number | Publication date |
---|---|
AU3229900A (en) | 2000-08-29 |
CA2362285A1 (en) | 2000-08-17 |
US6549854B1 (en) | 2003-04-15 |
NO20013894L (en) | 2001-10-12 |
WO2000048022A1 (en) | 2000-08-17 |
EP1151326B1 (en) | 2005-11-02 |
NO20013894D0 (en) | 2001-08-09 |
EP1151326A1 (en) | 2001-11-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2362285C (en) | Uncertainty constrained subsurface modeling | |
AU2013230933B2 (en) | Model predicting fracturing of shale | |
US6662109B2 (en) | Method of constraining by dynamic production data a fine model representative of the distribution in the reservoir of a physical quantity characteristic of the subsoil structure | |
US20200095858A1 (en) | Modeling reservoir permeability through estimating natural fracture distribution and properties | |
EP2810101B1 (en) | Improving efficiency of pixel-based inversion algorithms | |
US8359184B2 (en) | Method, program and computer system for scaling hydrocarbon reservoir model data | |
Dobróka et al. | Interval inversion approach for an improved interpretation of well logs | |
EP3788411B1 (en) | System and method for mapping hydrocarbon source rock using seismic attributes | |
EP4042211B1 (en) | Modeling reservoir permeability through estimating natural fracture distribution and properties | |
US20110315376A1 (en) | Gradient-Based Workflows For Conditioning Of Process-Based Geologic Models | |
US20180372896A1 (en) | Method of predicting parameters of a geological formation | |
US20220341835A1 (en) | Methods and Systems for Determining Reservoir and Fracture Properties | |
Gervais et al. | Probability maps of reservoir presence and sensitivity analysis in stratigraphic forward modeling | |
Orellana et al. | Influence of variograms in 3D reservoir-modeling outcomes: An example | |
EP3931600B1 (en) | Iterative stochastic seismic inversion | |
Poplavskii et al. | Two‐dimensional inverse modeling of sedimentary basin subsidence | |
US11927717B2 (en) | Optimized methodology for automatic history matching of a petroleum reservoir model with Ensemble Kalman Filter (EnKF) | |
US20100017131A1 (en) | Method for evaluating measured electromagnetic data relating to a subsurface region | |
Gutierrez et al. | Rock physics workflows for exploration in frontier basins | |
Jahani et al. | Real-Time Ensemble-Based Well-Log Interpretation for Geosteering | |
Cunha | Integrating static and dynamic data for oil and gas reservoir modelling | |
Gaillot et al. | Uncertainty analysis in formation evaluation: Rationale, methods and examples | |
US20230417949A1 (en) | Method for determining uncertainties associated with a model of a sedimentary basin | |
Ma et al. | Uncertainty Analysis | |
Stephen et al. | A multiple model approach to history matching and uncertainty analysis using time-lapse seismic |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
EEER | Examination request | ||
MKEX | Expiry |
Effective date: 20200211 |