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 PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims description 30
- 230000001186 cumulative effect Effects 0.000 claims description 2
- GNFTZDOKVXKIBK-UHFFFAOYSA-N 3-(2-methoxyethoxy)benzohydrazide Chemical compound COCCOC1=CC=CC(C(=O)NN)=C1 GNFTZDOKVXKIBK-UHFFFAOYSA-N 0.000 claims 1
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 4
- 238000002948 stochastic simulation Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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)
- 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; etf) 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; eti) 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.
- Procédé selon la revendication 1, selon lequel pour chaque emplacement, deux paramètres sont simulés en utilisant un variogramme transversal.
- 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.
- 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.
- Procédé selon les revendications 1 à 4, selon lequel un ou plusieurs paramètres sont mesurés dans un trou de sondage.
- 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; etretenir 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.
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)
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)
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)
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 |
-
1998
- 1998-03-06 DE DE69816152T patent/DE69816152T2/de not_active Expired - Lifetime
- 1998-03-06 ES ES98200712T patent/ES2199406T3/es not_active Expired - Lifetime
- 1998-03-06 DK DK98200712T patent/DK0864882T3/da active
- 1998-03-06 EP EP98200712A patent/EP0864882B1/fr not_active Expired - Lifetime
Cited By (1)
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 |