EP0864882B1 - Procédé pour estimer ou simuler les paramètres d'une structure stratifiée - Google Patents

Procédé pour estimer ou simuler les paramètres d'une structure stratifiée Download PDF

Info

Publication number
EP0864882B1
EP0864882B1 EP98200712A EP98200712A EP0864882B1 EP 0864882 B1 EP0864882 B1 EP 0864882B1 EP 98200712 A EP98200712 A EP 98200712A EP 98200712 A EP98200712 A EP 98200712A EP 0864882 B1 EP0864882 B1 EP 0864882B1
Authority
EP
European Patent Office
Prior art keywords
parameters
determined
stratum structure
density function
model
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
Application number
EP98200712A
Other languages
German (de)
English (en)
Other versions
EP0864882A2 (fr
EP0864882A3 (fr
Inventor
Hendrik Willem Johan Debeye
Eyad Sabbah
Pieter Maarten Van Der Made
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fugro NV
Original Assignee
Jason Geosystems BV
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jason Geosystems BV filed Critical Jason Geosystems BV
Publication of EP0864882A2 publication Critical patent/EP0864882A2/fr
Publication of EP0864882A3 publication Critical patent/EP0864882A3/fr
Application granted granted Critical
Publication of EP0864882B1 publication Critical patent/EP0864882B1/fr
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/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/61—Analysis by combining or comparing a seismic data set with other data
    • G01V2210/614—Synthetically generated data

Definitions

  • the present invention provides a method for estimating or simulating one or more parameters of a stratum structure, comprising the steps of:
  • geostatistical techniques are combined with deterministic techniques which make use of seismic data, whereby the uncertainty in the solution of the model can be considerably reduced.
  • the stratum structure is usually built up of sandstones, shale and/or limestone-type rocks. Significant parameters are for instance the specific density, the porosity, the permeability and the like of such strata.
  • the first model of the stratum structure on which iterative (better) estimations are performed is determined for instance using seismic, geological and/or well data.
  • the variogram of the parameter can be determined on the basis of the experience of a geologist, measurements on a well and/or from samples obtained from a well.
  • the parameter for instance the porosity
  • U a schematically designated area U (fig. 1) of a stratum structure, in which for instance two boreholes V and W are arranged
  • the parameter for instance the porosity
  • the parameter must be determined of a large number of points - x, x i , x j and x k are shown - while the property of the structure at locations v and w at the respective boreholes V and W are known with a certain accuracy, for instance because measurements are carried out with probes along the boreholes and/or samples are taken.
  • a structural (geometric) model is first created on the basis of seismic measurement data at the earth's surface. Reflection surfaces which are considered important are formed to a structure of closed strata, for instance by the computer program EARTH MODEL commercially available from applicant, wherein a geologist can make determined choices.
  • the program EARTH MODEL defines the geological layers as well as the micro-structure in such a layer. It will be apparent that according to the present invention other methods may also be envisaged to obtain an initial simulated stratum structure.
  • Variograms are subsequently defined for a particular parameter, such as the porosity in determined layers.
  • a variogram such as in fig. 2 the correlation or variance is determined between two values of a parameter given at different spatial locations, i.e. as a function of the distance between two locations.
  • a variogram between two different parameters, for instance between porosity and permeability, is referred to as a cross-variogram.
  • the variogram becomes horizontal, above which distance the spatial correlation amounts to zero.
  • the variance therefore has a constant value at a distance beyond this distance.
  • this distance has a different value in each case for different directions, thus forming as it were an ellipsoid.
  • a number of points, for instance a, b, c and d, where the value of this parameter is known (or is already present in the model and has been determined previously in the same manner) are chosen in the vicinity of this point x 1 .
  • a value for the parameter is chosen randomly and stored for the point x 1 .
  • the thus simulated point is added to the known value and the procedure is repeated until all points x i in the space U have been simulated and a value assigned thereto.
  • the above stated preferred method differs from the so-called Kriging method (or in the case of cross-variograms co-Kriging) in that therein the most probable value is determined from the variograms (see C.V. Deutsch and A.G. Journel (1992), GSLIB: geostatistical software library and user's guide, Oxford University Press, New York).
  • the reservoir models fit available geostatistical information, such as variograms and histograms (i.e. wherein the data has a Gaussian distribution).
  • the models are subsequently checked against seismic data.
  • the parameters being estimated exclusively according to stochastic simulations, after obtaining a first stochastic model it is determined using seismic data whether or not a new iteration point fits better in the model than the previously found point.
  • a grid node (20) is selected at random, the conditional local probability density function is then determined subject to a particular lithology (indicator) (21), whereafter the realization lithology from the local probability density function is drawn (22) until the random path (23) is complete.
  • Fig. 6 shows the preferred embodiment, wherein the procedures of fig. 4 and 5 are combined and wherein the convergence step for the seismic data is shown in block 31.

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Claims (6)

  1. Procédé pour estimer ou simuler un ou plusieurs paramètres d'une structure stratifiée comprenant les étapes consistant à:
    a) obtenir des données sismiques de la structure stratifiée;
    b) déterminer à partir des données sismiques obtenues, un premier modèle de la structure stratifiée;
    c) déterminer un ou plusieurs variogrammes d'un ou de plusieurs paramètres dans une ou plusieurs strates de la structure stratifiée;
    d) mesurer une valeur pour le paramètre ou la pluralité de paramètres en une pluralité de points dans la structure stratifiée;
    e) estimer, sur la base du variogramme ou de la pluralité de variogrammes, une fonction de densité de probabilité (PDF) dudit paramètre ou de ladite pluralité de paramètres en d'autres points dans les strates adjacents ou distants des points où ledit paramètre ou ladite pluralité de paramètres ont été mesurés; et
    f) déterminer, en utilisant la fonction de densité de probabilité, une valeur pour ledit paramètre ou chacun de ladite pluralité de paramètres dans les strates au niveau des points où la fonction de densité de probabilité des paramètres a été déterminée;
    g) combiner les valeurs déterminées lors de l'étape f) audit paramètre correspondant ou à ladite pluralité de paramètres correspondants, déterminés à partir du modèle;
    h) obtenir à partir dudit paramètre ou de ladite pluralité de paramètres, déterminés lors de l'étape g), des valeurs acceptables pour ledit paramètre ou ladite pluralité de paramètres dans la structure stratifiée; et
    i) produire à partir des valeurs acceptables pour ledit paramètre ou ladite pluralité de paramètres un nouveau modèle;
    j) répéter les étapes (d)-(i) jusqu'à obtenir une estimation de la précision désirée.
  2. Procédé selon la revendication 1, selon lequel pour chaque emplacement, deux paramètres sont simulés en utilisant un variogramme transversal.
  3. Procédé selon la revendication 1 ou 2, selon lequel la fonction de densité de probabilité est sélectionnée sur la base d'un indicateur pour le type de lithologie à laquelle on doit s'attendre.
  4. Procédé selon les revendications 1 à 3, selon lequel le variogramme ou la pluralité de variogrammes sont déterminés sur la base des données disponibles d'une strate obtenues à partir d'échantillons de cette dernière.
  5. Procédé selon les revendications 1 à 4, selon lequel un ou plusieurs paramètres sont mesurés dans un trou de sondage.
  6. Procédé selon la revendication 1, selon lequel l'étape (h) comprend les étapes consistant à:
    former la convolution du paramètre ou de la pluralité de paramètres déterminés lors de l'étape (g) avec une ondelette pour obtenir un ou plusieurs paramètres traités par convolution;
    comparer le paramètre ou la pluralité de paramètres traités par convolution à un ou plusieurs paramètres correspondants déterminés à partir du modèle antérieur; et
    retenir le paramètre ou la pluralité de paramètres traités par convolution en tant que valeurs acceptables si le paramètre ou la pluralité de paramètres traités par convolution améliorent une probabilité cumulée de la fonction de densité de probabilité pour le paramètre ou la pluralité de paramètres à partir du modèle antérieur.
EP98200712A 1997-03-10 1998-03-06 Procédé pour estimer ou simuler les paramètres d'une structure stratifiée Expired - Lifetime EP0864882B1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US81427697A 1997-03-10 1997-03-10
US814276 1997-03-10

Publications (3)

Publication Number Publication Date
EP0864882A2 EP0864882A2 (fr) 1998-09-16
EP0864882A3 EP0864882A3 (fr) 2000-12-27
EP0864882B1 true EP0864882B1 (fr) 2003-07-09

Family

ID=25214600

Family Applications (1)

Application Number Title Priority Date Filing Date
EP98200712A Expired - Lifetime EP0864882B1 (fr) 1997-03-10 1998-03-06 Procédé pour estimer ou simuler les paramètres d'une structure stratifiée

Country Status (4)

Country Link
EP (1) EP0864882B1 (fr)
DE (1) DE69816152T2 (fr)
DK (1) DK0864882T3 (fr)
ES (1) ES2199406T3 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7373251B2 (en) 2004-12-22 2008-05-13 Marathon Oil Company Method for predicting quantitative values of a rock or fluid property in a reservoir using seismic data

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5915308A (en) * 1996-01-18 1999-06-29 Siemens Aktiengesellschaft Discharge apparatus
US6480790B1 (en) 1999-10-29 2002-11-12 Exxonmobil Upstream Research Company Process for constructing three-dimensional geologic models having adjustable geologic interfaces
FR2846767B1 (fr) * 2002-10-30 2004-12-24 Inst Francais Du Petrole Methode pour former plus rapidement un modele stochastique representatif d'un reservoir heterogene souterrain, contraint par des donnees dynamiques
GB0228884D0 (en) 2002-12-11 2003-01-15 Schlumberger Holdings Method and system for estimating the position of a movable device in a borehole
FR2869421B1 (fr) * 2004-04-27 2006-06-02 Inst Francais Du Petrole Methode de reconstruction d'un modele stochastique, representatif d'un milieu heterogene poreux, pour ameliorer son calage par les donnees de production
JP5379163B2 (ja) 2008-01-08 2013-12-25 エクソンモービル アップストリーム リサーチ カンパニー 地震探査データのスペクトルシェーピングインバージョン法及びマイグレーション法
WO2010151354A1 (fr) 2009-06-26 2010-12-29 Exxonmobil Upstream Research Company Construction de modèles de résistivité à partir d'une inversion stochastique
CA2855137A1 (fr) * 2011-11-15 2013-05-23 Chevron U.S.A. Inc. Systeme et procede d'utilisation de sous-ensembles spatialement independants de donnees pour determiner l'incertitude de la correction de biais de donnees calculees de distributions de proprietes pour des donnees de reservoirs spatialement correlees

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2710418B1 (fr) * 1993-09-21 1995-12-15 Inst Francais Du Petrole Méthode d'analyse des traces sismiques utilisant une technique de calibrage statistique pour en déduire des propriétés géologiques.
US5539704A (en) * 1995-06-23 1996-07-23 Western Atlas International, Inc. Bayesian sequential Gaussian simulation of lithology with non-linear data

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7373251B2 (en) 2004-12-22 2008-05-13 Marathon Oil Company Method for predicting quantitative values of a rock or fluid property in a reservoir using seismic data

Also Published As

Publication number Publication date
DE69816152T2 (de) 2004-05-27
EP0864882A2 (fr) 1998-09-16
ES2199406T3 (es) 2004-02-16
EP0864882A3 (fr) 2000-12-27
DK0864882T3 (da) 2003-08-18
DE69816152D1 (de) 2003-08-14

Similar Documents

Publication Publication Date Title
US11519262B2 (en) Systematic evaluation of shale plays
EP0674189B1 (fr) Méthodes de génération des modèles numériques à partir de l'information sismique
EP0750203B1 (fr) Modélisation de subsurface à partir de données sismiques et de mesures secondaires
RU2271023C2 (ru) Способ предсказания количественных значений свойства скальной породы или текучей среды в продуктивном пласте с помощью сейсмических данных
CA2543801C (fr) Procedes de fabrication de modele de reservoir
US8498848B2 (en) Method for upscaling a reservoir model using deep reading measurements
CN111273348B (zh) 基于更新概率比率恒定理论的多点地质统计叠前反演方法
JP3390445B2 (ja) 地質学的特性を推定するために統計学的較正技術を使用する地震トレース解析方法
CN104755960B (zh) 基于盆地建模对用于处理地震数据的速度模型的改进
US20150247940A1 (en) History matching of time-lapse crosswell data using ensemble kalman filtering
US8874419B2 (en) Method of developing a petroleum reservoir from a facies map construction
AU6257198A (en) Method for determining barriers to reservoir flow
US8965744B2 (en) Method of developing a petroleum reservoir by reservoir model reconstruction
EP0864882B1 (fr) Procédé pour estimer ou simuler les paramètres d'une structure stratifiée
US20210405250A1 (en) Conditioning of Surface-Based Geologic Models
WO2016187238A1 (fr) Auto-validation d'un système d'interprétation et de modélisation terrestre
Noroozi et al. 3D Geometrical-Stochastical modeling of rock mass joint networks: case study of the right bank of Rudbar Lorestan Dam plant
Vidal et al. Geostatistical simulations of geothermal reservoirs: two-and multiple-point statistic models
GhojehBeyglou et al. A geostatistical analysis approach for facies and porosity modelling of a heterogeneous sandstone reservoir to compare four practical stochastic methods.
US20230359793A1 (en) Machine-learning calibration for petroleum system modeling
Madani et al. Improvement of fracture network modeling in fractured reservoirs using conditioning and geostatistical method
Almeida et al. An integrated approach to reservoir studies using stochastic simulation techniques
Niven Geostatistics for naturally fractured reservoirs
Bisht et al. Generation of regional trend from Boolean log of lithofacies for building a realistic fine scale geological model using stochastic facies modelling–A case study
Fitzgerald et al. Delineate Three-Dimensional Iron Ore Geology and Resource Models Using the Potential Field Method

Legal Events

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

Free format text: ORIGINAL CODE: 0009012

AK Designated contracting states

Kind code of ref document: A2

Designated state(s): DE DK ES FR GB IT NL

AX Request for extension of the european patent

Free format text: AL;LT;LV;MK;RO;SI

PUAL Search report despatched

Free format text: ORIGINAL CODE: 0009013

AK Designated contracting states

Kind code of ref document: A3

Designated state(s): AT BE CH DE DK ES FI FR GB GR IE IT LI LU MC NL PT SE

AX Request for extension of the european patent

Free format text: AL;LT;LV;MK;RO;SI

17P Request for examination filed

Effective date: 20010530

AKX Designation fees paid

Free format text: DE ES FR GB IT NL

RBV Designated contracting states (corrected)

Designated state(s): DE DK ES FR GB IT NL

GRAH Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOS IGRA

GRAH Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOS IGRA

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

AK Designated contracting states

Designated state(s): DE DK ES FR GB IT NL

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

REF Corresponds to:

Ref document number: 69816152

Country of ref document: DE

Date of ref document: 20030814

Kind code of ref document: P

REG Reference to a national code

Ref country code: DK

Ref legal event code: T3

REG Reference to a national code

Ref country code: ES

Ref legal event code: FG2A

Ref document number: 2199406

Country of ref document: ES

Kind code of ref document: T3

ET Fr: translation filed
PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

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

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

26N No opposition filed

Effective date: 20040414

REG Reference to a national code

Ref country code: GB

Ref legal event code: 732E

NLS Nl: assignments of ep-patents

Owner name: FUGRO N.V.

Effective date: 20070403

REG Reference to a national code

Ref country code: FR

Ref legal event code: TP

REG Reference to a national code

Ref country code: ES

Ref legal event code: PC2A

REG Reference to a national code

Ref country code: DE

Ref legal event code: R082

Ref document number: 69816152

Country of ref document: DE

Representative=s name: DREISS PATENTANWAELTE PARTG MBB, DE

Ref country code: DE

Ref legal event code: R082

Ref document number: 69816152

Country of ref document: DE

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 18

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

Ref country code: NL

Payment date: 20150213

Year of fee payment: 18

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

Ref country code: DK

Payment date: 20150219

Year of fee payment: 18

Ref country code: DE

Payment date: 20150311

Year of fee payment: 18

Ref country code: IT

Payment date: 20150325

Year of fee payment: 18

Ref country code: ES

Payment date: 20150324

Year of fee payment: 18

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

Ref country code: FR

Payment date: 20150313

Year of fee payment: 18

Ref country code: GB

Payment date: 20150316

Year of fee payment: 18

REG Reference to a national code

Ref country code: DE

Ref legal event code: R081

Ref document number: 69816152

Country of ref document: DE

Owner name: CGG SERVICES (NL) B.V., NL

Free format text: FORMER OWNER: FUGRO N.V., LEIDSCHENDAM, NL

REG Reference to a national code

Ref country code: GB

Ref legal event code: 732E

Free format text: REGISTERED BETWEEN 20160114 AND 20160120

REG Reference to a national code

Ref country code: DE

Ref legal event code: R119

Ref document number: 69816152

Country of ref document: DE

REG Reference to a national code

Ref country code: DK

Ref legal event code: EBP

Effective date: 20160331

REG Reference to a national code

Ref country code: NL

Ref legal event code: MM

Effective date: 20160401

GBPC Gb: european patent ceased through non-payment of renewal fee

Effective date: 20160306

REG Reference to a national code

Ref country code: FR

Ref legal event code: ST

Effective date: 20161130

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

Ref country code: DE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20161001

Ref country code: GB

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20160306

Ref country code: NL

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20160401

Ref country code: FR

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20160331

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

Ref country code: IT

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20160306

REG Reference to a national code

Ref country code: ES

Ref legal event code: FD2A

Effective date: 20170428

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

Ref country code: DK

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20160331

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

Ref country code: ES

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20160307