WO2021130512A1 - Dispositif et procédé de prédiction de valeurs de porosité, de lithofaciès et de perméabilité dans un réservoir étudié de carbonates en fonction de données sismiques - Google Patents
Dispositif et procédé de prédiction de valeurs de porosité, de lithofaciès et de perméabilité dans un réservoir étudié de carbonates en fonction de données sismiques Download PDFInfo
- Publication number
- WO2021130512A1 WO2021130512A1 PCT/IB2019/001442 IB2019001442W WO2021130512A1 WO 2021130512 A1 WO2021130512 A1 WO 2021130512A1 IB 2019001442 W IB2019001442 W IB 2019001442W WO 2021130512 A1 WO2021130512 A1 WO 2021130512A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- values
- lithofacies
- porosity
- permeability
- prediction model
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 154
- 230000035699 permeability Effects 0.000 title claims abstract description 113
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 99
- 230000009466 transformation Effects 0.000 claims abstract description 21
- 238000000844 transformation Methods 0.000 claims abstract description 21
- 230000015654 memory Effects 0.000 claims description 34
- 238000012549 training Methods 0.000 claims description 31
- 238000004891 communication Methods 0.000 claims description 23
- 150000004649 carbonic acid derivatives Chemical class 0.000 claims description 4
- 238000012952 Resampling Methods 0.000 claims description 2
- 238000010801 machine learning Methods 0.000 description 38
- 238000007637 random forest analysis Methods 0.000 description 28
- 238000004422 calculation algorithm Methods 0.000 description 20
- 238000010586 diagram Methods 0.000 description 15
- 238000013528 artificial neural network Methods 0.000 description 14
- 239000011435 rock Substances 0.000 description 13
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 12
- 230000008569 process Effects 0.000 description 12
- 238000012706 support-vector machine Methods 0.000 description 12
- 235000019738 Limestone Nutrition 0.000 description 11
- 239000006028 limestone Substances 0.000 description 11
- 238000004458 analytical method Methods 0.000 description 9
- 238000012512 characterization method Methods 0.000 description 8
- 238000012545 processing Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 238000004590 computer program Methods 0.000 description 6
- 239000000377 silicon dioxide Substances 0.000 description 6
- 238000003860 storage Methods 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 5
- 238000009826 distribution Methods 0.000 description 5
- 238000005755 formation reaction Methods 0.000 description 5
- 208000035126 Facies Diseases 0.000 description 4
- 230000008901 benefit Effects 0.000 description 4
- 230000005251 gamma ray Effects 0.000 description 4
- 229930195733 hydrocarbon Natural products 0.000 description 4
- 150000002430 hydrocarbons Chemical class 0.000 description 4
- 238000012417 linear regression Methods 0.000 description 4
- 239000003208 petroleum Substances 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000002776 aggregation Effects 0.000 description 3
- 238000004220 aggregation Methods 0.000 description 3
- 229910052925 anhydrite Inorganic materials 0.000 description 3
- 238000013459 approach Methods 0.000 description 3
- OSGAYBCDTDRGGQ-UHFFFAOYSA-L calcium sulfate Chemical compound [Ca+2].[O-]S([O-])(=O)=O OSGAYBCDTDRGGQ-UHFFFAOYSA-L 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 230000004931 aggregating effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000002860 competitive effect Effects 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000003066 decision tree Methods 0.000 description 2
- 238000007636 ensemble learning method Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013145 classification model Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 238000001953 recrystallisation Methods 0.000 description 1
- -1 salt carbonates Chemical class 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000011282 treatment Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
-
- 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/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6244—Porosity
-
- 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/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6246—Permeability
Definitions
- the present invention relates to the field of exploration and appreciation of reservoirs.
- the invention relates to the field of exploration and appreciation of carbonate reservoirs using machine learning techniques on seismic data. More particularly, the invention provides a new method for predicting values of porosity, lithofacies and/or permeability, in a studied carbonate reservoir.
- Reservoir characterization has become increasingly important to hydrocarbon exploration
- Lithology and reservoir properties prediction from seismic data plays an essential role in reservoir quality evaluation, reservoir architecture delineation, and reservoir model building, which is of long-standing interest for petroleum reservoir exploration, development and production.
- Reservoir characterization attempts to describe petroleum deposits and the nature of the rocks that contain hydrocarbons, producing detailed geological reconstructions of both its geometry and of its lithological properties. This can provide important decision support, in particular in the highly competitive segment of carbonate reservoirs.
- Reservoir characterization relies on expertise from petroleum engineering, geology, and geophysics and can benefit from well logging.
- deterministic method and geostatistical method are both employed to estimate lithofacies and reservoir properties, such as porosity, from seismic data.
- Lithology and texture identification are usually performed using either a core sample or cutting analysis and the borehole log methods. Core samples are collected during the well-drilling process and are analyzed directly by experts, i.e., geologists, but this technique is costly. Indirect methods, i.e. well logs, provide information that can be used for rock characterization with relatively lower costs.
- the lithofacies and reservoir parameter determination problem entails some important issues to be considered: data are intrinsically noisy and imperfect; the physics between the reservoir properties and seismic responses are highly complex and nonlinear. This is especially true for the heterogeneous carbonate reservoirs.
- well log is at one-point location for exploration stage and is incapable of proposition a 3D analysis.
- Machine learning provides an intelligent and practical means of predicting geological features, or any spatially varying physical properties, from multi-dimensional geophysical data sets.
- the basic premise of supervised learning is that it requires training data containing labeled samples representing what is known about the inference target. T rained classification models are then applied to input variables with similar geological context to predict classes present within the training data.
- pre-salt oilfields For pre-salt oilfields, owing to the complex porosity distribution in carbonate reservoirs, predicting a reliable porosity is a fundamental step for reservoir modelling.
- pre salt carbonates are very heterogeneous reservoirs, in terms of facies, and consequently, in terms of porosity and permeability. Processes such as diagenesis and recrystallization can modify the primary porosity and make the lithofacies heterogeneous.
- applying a porosity volume derived from seismic data directly to the geological model is not an easy task since the difference between the seismic and geological grid could create upscaling issues both in vertical and lateral domains. Mori et al.
- the invention aims to overcome the disadvantages of the prior art.
- the invention proposes a method for predicting values of porosity, lithofacies and/or permeability in a studied carbonate reservoir based on seismic data, said method allowing to directly connect lithology and reservoir properties to seismic data.
- a method based on machine learning techniques allows to predict the spatial distribution of porosity, lithofacies and/or permeability from seismic data.
- the invention also proposes a computer device configured to predict values of porosity, lithofacies and/or permeability in a carbonate reservoir based on seismic data.
- a solution according to the invention can produce reliable values from a complex and heterogeneous carbonate reservoir.
- a method for predicting values of porosity, lithofacies and/or permeability in a studied carbonate reservoir based on seismic data comprising:
- an ensemble prediction model said ensemble prediction model having been trained with: o values of elastic attributes, said elastic attributes comprising V p , V s , density or their mathematic transformations, and o values of porosity, lithofacies and/or permeability;
- inverted elastic attributes calculated from seismic data generated on a studied carbonate reservoir, said inverted elastic attributes comprising V p , V s , density or their mathematic transformations;
- a method according to the invention is based on elastic attributes that are more sensitive to variations of lithofacies and reservoir properties.
- the ensemble prediction model used in this method has been trained with well logging data to generate a model that represent the complex nonlinear relationship between lithofacies/porosity/permeability and elastic attributes. Then the invention comprises the combination of the trained ensemble prediction model with inverted seismic cubes from seismic data (P-impedance and Vp/Vs ratio) to predict the spatial distribution of porosity, lithofacies and/or permeability from seismic data.
- the step of training the ensemble prediction model comprises a resampling of the values of porosity, lithofacies and/or permeability in order to reduce the frequency of data used in training.
- this comprises the use of substantially identical frequencies between data used in training and inverted elastic attributes generated from seismic data.
- a moving average filter to smooth the values of porosity, lithofacies and/or permeability.
- the elastic attributes comprise elastic attributes selected from the group consisting of: V p , V s , density, V p /V s ratio, shear modulus, bulk modulus, P-impedance, S- impedance, Poisson’s ratio and Lame’s coefficient;
- the ensemble prediction model includes any one of: a boosting such as gradient boosting or adaptative boosting, a bagging and a stacking.
- the ensemble prediction model includes any one of: random forest, adaboost and xgboost.
- the ensemble prediction model has been trained on one, preferably several, carbonate reservoirs that are different from the studied carbonate reservoir on which seismic data are generated.
- the method according to the invention do not need to realize log well on the studied reservoir to generate the predicted data.
- well log can be integrated on the studied carbonate reservoir to enhance the prediction model.
- the ensemble prediction model has been trained on one, preferably several, carbonate reservoirs that have a similar geological context.
- the step of retraining can comprise the use of key information on the formation of the reservoir such as water contact or gamma ray.
- the ensemble prediction model also includes reservoir type as input data, for example said reservoir type being selected from: lacustrine carbonates reservoir, cold-water carbonate reservoir, and warm-water carbonate reservoir.
- reservoir type being selected from: lacustrine carbonates reservoir, cold-water carbonate reservoir, and warm-water carbonate reservoir.
- the seismic data acquired for a carbonate reservoir correspond to 2D or 3D seismic data
- - generating predicted values of porosity, lithofacies and/or permeability comprises a generation of predicted values of porosity, lithofacies and/or permeability for all intersections of a resolution grid corresponding to 2D seismic data acquired,
- a computer device for predicting values of porosity, lithofacies and/or permeability in a studied carbonate reservoir based on seismic data comprising:
- a data memory configured to store an ensemble prediction model, said ensemble prediction model having been trained with: o values of elastic attributes, said elastic attributes comprising V p , V s , density or their mathematic transformations, and o values of porosity, lithofacies and/or permeability;
- a communication interface configured to acquire values of inverted elastic attributes calculated from seismic data generated on a studied carbonate reservoir, said inverted elastic attributes comprising V p , V s , density or their mathematic transformations;
- a processor configured to: o Load the ensemble prediction model, and o Generate predicted values of porosity, lithofacies and/or permeability in the studied carbonate reservoir, said predicted values being calculated from the acquired values of inverted elastic attributes and the ensemble prediction model.
- a non-transitory computer readable medium storing executable instructions which, when executed by a processor, implements a method according to the invention.
- a non-transitory computer readable medium storing executable instructions which, when executed by a processor of a computer device, implements a method for predicting values of porosity, lithofacies and/or permeability in a studied carbonate reservoir based on seismic data, the method comprising:
- an ensemble prediction model said ensemble prediction model having been trained with: o values of elastic attributes, said elastic attributes comprising V p , V s , density or their mathematic transformations, and o values of porosity, lithofacies and/or permeability;
- inverted elastic attributes calculated from seismic data generated on a studied carbonate reservoir, said inverted elastic attributes comprising V p , V s , density or their mathematic transformations; and Generating predicted values of porosity, lithofacies and/or permeability of the studied carbonate reservoir, said predicted values being calculated from the acquired values of inverted elastic attributes and the ensemble prediction model.
- FIG. 1 is a schematic view of a process flow diagram of the method according to an aspect of the invention.
- FIG. 2 is a schematic view of a specific step of the method according to an embodiment of the invention.
- FIG. 3 is a schematic view of a process flow diagram of a specific step of the method according to an embodiment of the invention.
- FIG. 4 is a schematic view showing a process flow diagram of a method according to another embodiment of the invention.
- FIG. 5 is a schematic block diagram of a computer device for predicting values according to an aspect of the invention.
- FIG. 6 is an illustration of a comparison of lithofacies prediction using different machine learning algorithms (6A: Fuzzy Logic, probabilistic neural network (PNN), support vector machine (SVM) - 6B : deep neural network (DNN), Naive Bayes, Random Forest, extreme gradient boosting (XGBOOST)). Each grey level corresponds to different lithofacies.
- machine learning algorithms 6A: Fuzzy Logic, probabilistic neural network (PNN), support vector machine (SVM) - 6B : deep neural network (DNN), Naive Bayes, Random Forest, extreme gradient boosting (XGBOOST)).
- PNN probabilistic neural network
- SVM support vector machine
- DNN deep neural network
- Naive Bayes Naive Bayes
- Random Forest Random Forest
- XGBOOST extreme gradient boosting
- FIG. 7 is an illustration of a confusing matrix of prediction results of five lithofacies (igneous, anhydrite, dolomitic limestone, silica-rich limestone, and limestone) using different machine learning algorithms: Random Forest (7A) ; extreme gradient boosting (7B); Fuzzy Logic (7C); Deep Neural Network (7D); Probabilistic Neural Network (7E); Naive Bayes (7F); Support Vector Machine (7G)).
- Random Forest (7A) ; extreme gradient boosting (7B); Fuzzy Logic (7C); Deep Neural Network (7D); Probabilistic Neural Network (7E); Naive Bayes (7F); Support Vector Machine (7G)).
- FIG. 8 is a schema of a comparison of different machine learning algorithm for porosity prediction: Support Vector Machine (8A); Deep Neural Network (8B) ; Linear regression (8C); Nonlinear regression (8D) ; Random Forest (8E) ; extreme gradient boosting (8F) of 3 wells (A1 , A2, A3).
- Support Vector Machine (8A) Deep Neural Network
- Linear regression (8C) Linear regression
- Nonlinear regression (8D) Nonlinear regression
- Random Forest (8E) extreme gradient boosting (8F) of 3 wells (A1 , A2, A3).
- each box in the flow diagrams or block diagrams may represent a system, a device, a module or code which comprises several executable instructions for implementing the specified logical function(s).
- the functions associated with the box may appear in a different order than indicated in the figures.
- two boxes successively shown may be executed substantially simultaneously, or boxes may sometimes be executed in the reverse order, depending on the functionality involved.
- Each box of flow diagrams or block diagrams and combinations of boxes in flow diagrams or block diagrams may be implemented by special systems that perform the specified functions or actions or perform combinations of special equipment and computer instructions.
- ensemble prediction method means a type of machine learning method that uses a series of learners to learn and uses some rules to integrate the learning results so as to achieve better learning effects than a single learner.
- the main idea of ensemble learning is to first generate a number of learners according to certain rules and then combine them by some integration strategies, and eventually output the final results by comprehensive judgment. Briefly, what ensemble learning does is to integrate multiple weak learners into one strong learner.
- Such ensemble prediction method can for example be selected from: Random Forest (RF, Breiman, 1996; 2001) and Extreme Gradient Boosting (XGBOOST, Chen and Guestrin, 2016).
- An ensemble prediction method is used to produce an ensemble prediction model.
- reservoir or “petroleum reservoir” can refers to a subsurface group of sedimentary rocks capable of storing a pool of hydrocarbons. This is commonly a porous sandstone or limestone.
- process By “process”, “compute”, “determine”, “display”, “extract”, “compare” or more broadly “executable operation” is meant, within the meaning of the invention, an action performed by a computer device or a processor unless the context indicates otherwise.
- the operations relate to actions and/or processes of a data processing system, for example a computer system or an electronic computing device, which manipulates and transforms the data represented as physical (electronic) quantities in the memories of the computer system or other devices for storing, transmitting or displaying information.
- calculation operations are carried out by the processor of the device, the produced data are entered in a corresponding field in a data memory and this field or these fields can be returned to a user for example through a Human Machine Interface formatting such data.
- These operations may be based on applications or software.
- application means any expression, code or notation, of a set of instructions intended to cause a data processing to perform a particular function directly or indirectly (for example after a conversion operation into another code).
- program codes may include, but are not limited to, a subprogram, a function, an executable application, a source code, an object code, a library and/or any other sequence of instructions designed for being performed on a computer system.
- processor is meant, within the meaning of the invention, at least one hardware circuit configured to perform operations according to instructions contained in a code.
- the hardware circuit may be an integrated circuit. Examples of a processor include, but are not limited to, a central processing unit, a graphics processor, an application-specific integrated circuit (“ASIC” according to Anglo-Saxon terminology), and a programmable logic circuit. A single processor or several other units may be used to implement the invention.
- a single processor or several other units may be used to implement the invention.
- coupled is meant, within the meaning of the invention, connected, directly or indirectly, with one or more intermediate elements. Two elements may be coupled mechanically, electrically or linked by a communication channel.
- human-machine interface corresponds to any element allowing a human being to communicate with a computer, in particular and without that list being exhaustive, a keyboard and means allowing in response to the commands entered on the keyboard to perform displays and optionally to select with the mouse or a touchpad items displayed on the screen.
- a touch screen for selecting directly on the screen the elements touched by the finger or an object and optionally with the possibility of displaying a virtual keyboard.
- computer device any device comprising a processing unit or a processor, for example in the form of a microcontroller cooperating with a data memory, possibly a program memory, said memories possibly being dissociated.
- the processing unit cooperates with said memories by means of internal communication bus.
- substantially refers to a majority of, or mostly, as in at least about 50%, 60%, 70%, 80%, 90%, 95%, 96%, 97%, 98%, 99%, 99.5%, 99.9%, 99.99%, or at least about 99.999% or more.
- carbonate reservoirs constitute a highly competitive segment for which physics between the reservoir properties and seismic responses are highly complex and nonlinear.
- Machine learning methods have been proposed for lithofacies prediction or porosity prediction. However, most methods are based on rock sample analysis or show a high level of confusion between lithofacies such as limestone classes. Moreover, no machine learning methods have been proposed for permeability prediction from inverted seismic data.
- the inventor developed solutions dedicated to carbonate reservoirs and based on machine learning to determine porosity, lithofacies and/or permeability from seismic data.
- well log data are used to build an ensemble prediction model to classify lithofacies using elastic attributes.
- Inverted elastic attributes coming from seismic data, preferably 2D or 3D, are processed by the trained prediction model to generate lithofacies, permeability and porosity prediction for example on the whole 3D area.
- the invention relates to a method 100 for predicting values of porosity, lithofacies and/or permeability in a studied carbonate reservoir based on seismic data.
- said method comprising the steps of: loading 130 an ensemble prediction model, acquiring 170 values of inverted elastic attributes calculated from seismic data generated on a studied carbonate reservoir, and generating 180 predicted values of porosity, lithofacies and/or permeability of the studied carbonate reservoir.
- a method 100 for predicting values of porosity, lithofacies and/or permeability according to the invention is mainly based on:
- This method is particularly adapted for carbonate reservoirs and, as illustrated by example, produces better results than other prediction methods already proposed on clastic or shale reservoirs for example.
- the ensemble prediction model has been previously trained.
- the invention set out herein can also comprise an embodiment where the ensemble prediction model is trained as part of the method according to the invention.
- a method 100 for predicting values of porosity, lithofacies and/or permeability according to the invention can comprise a facultative step of collecting 110 data from a well log.
- a method 100 for predicting values of porosity, lithofacies and/or permeability according to the invention can comprise a step of training 120 the ensemble prediction model.
- Machine learning is now widely adopted in various industrial fields. For example, Bayesian classification has been used for years in seismic data analysis. Prediction models can be divided into unsupervised learning methods and supervised learning methods.
- the unsupervised learning methods make it possible to determine groups of observations without a priori. Hence, those groups will be formed without a need for a label value on input data.
- the supervised learning methods link an input to an output based on example input-output pairs.
- a machine learning technique is used to build a supervised prediction model configured to estimate lithofacies, permeability and/or porosity from inverted elastic attributes.
- supervised learning methods neural networks, classification or regression trees, nearest neighbor search, and random forest are some of the most robust and efficient machine learning techniques according to the invention. For example, the Bayesian classification which has been used for years in seismic data analysis shows lower prediction result than ensemble methods evaluated.
- an ensemble prediction model that estimates lithofacies, permeability and/or porosity from elastic attributes, preferably inverted elastic attributes.
- the training of said ensemble prediction model is preferably done with values of elastic attributes, as input data, and with values of lithofacies, permeability and/or porosity, as target data.
- These values used to train the ensemble prediction model can be considered as reference data. They have advantageously been produced on a carbonate reservoir and have been inferred from a well log.
- Elastic attributes used to train the ensemble prediction model can be inverted elastic attributes, namely elastic attributes calculated with inversion of geophysical data, or elastic attributes generated from well log data.
- the ensemble prediction model has been trained on one, preferably several, carbonate reservoirs that are different from the studied carbonate reservoir on which seismic data is generated.
- the ensemble prediction model has been trained at least with the attribute to be predicted (i.e. porosity, lithofacies and permeability). For example, if the invention is used to predict porosity and lithofacies values, the ensemble model should have been trained at least with porosity and lithofacies values.
- the ensemble prediction model preferably can also include reservoir types as input data, for example said reservoir types being selected from: lacustrine carbonates formation, cold-water carbonate formation, and warm-water carbonate formation.
- reservoir types being selected from: lacustrine carbonates formation, cold-water carbonate formation, and warm-water carbonate formation.
- the ensemble prediction model is not calibrated/trained based on rock analysis of samples and preferably uses only well logs for training.
- the training 120 of the ensemble prediction model can be based on key information values 121 , elastic attributes values 122, porosity values 123, permeability values 125 and lithofacies values 124.
- the training 120 of the ensemble prediction model is at least based on elastic attributes values, porosity values, permeability values and lithofacies values.
- Key information values 121 can for example correspond to reservoirs typology type, a priori geological information, gamma ray values or water-contact values. Such key information without being essential can enhance the accuracy of the prediction.
- ensemble prediction models proved to be the best performers for permeability, lithofacies and porosity predictions on carbonate reservoirs. Indeed, while artificial neural networks have already been proposed for porosity prediction and lithofacies predictions, ensemble prediction models outperformed several other methods of machine learning in this context of carbonate reservoirs.
- Ensemble prediction model refers to a machine learning paradigm where multiple models (i.e. “weak learners”) are trained to solve the same problem and are combined to obtain better results. The main hypothesis is that when weak models are correctly combined, more accurate and/or robust models can be obtained.
- Ensemble prediction models can be combined in a “homogeneous ensemble prediction model”, homogeneous weak learners that are trained in different ways.
- an ensemble prediction model can use different types of base learning algorithms to combine heterogeneous weak learners into a “heterogeneous ensemble prediction model”.
- Ensemble prediction methods can also be divided into sequential ensemble prediction methods, where the base learners are generated sequentially (e.g. xgboost), and parallel ensemble prediction methods, where the base learners are generated in parallel (e.g. Random Forest).
- Stacking is an ensemble learning technique that combines multiple classification or regression models via a meta-classifier or a meta-regressor.
- the base level models are trained based on a complete training set, then the meta-model is trained on the outputs of the base level model as features.
- the base level often consists of different learning algorithms and therefore stacking ensembles are often heterogeneous.
- Bagging uses bootstrap sampling to obtain the data subsets for training the base learners.
- each model in the ensemble is trained on a random sample of the dataset, where preferably each random sample is the same size as the dataset and sampling with replacement is used.
- bagging uses voting for classification and averaging for regression.
- the random forest approach can be considered as being a bagging method where decision trees, fitted on bootstrap samples, are combined to produce an output with lower variance.
- Boosting refers to a family of algorithms that can fit a sequence of weak learners (such as small decision trees) to weighted versions of the data.
- each new model added to an ensemble is biased to pay more attention (more weight) to instances that previous model misclassified.
- the predictions are then combined through a weighted majority vote or a weighted sum to produce the final prediction.
- the principal difference between boosting and bagging, is that boosting works by iteratively creating models and adding them to the ensemble.
- Boosting comprise for example adaptative boosting and gradient boosting.
- the ensemble prediction model includes a method selected from: stacking, boosting, such as gradient boosting or adaptative boosting, and bagging such as random forest.
- the ensemble prediction model is selected from: random forest, xgboost, adaboost. Even more preferably, the ensemble prediction model is selected from xgboost and adaboost.
- the step of training 120 the ensemble prediction model can comprise the use of a moving average filter to smooth the values of permeability, lithofacies and/or porosity calculated from well logs. This will make it comparable to the scale of the inverted elastic attributes from pre-stack seismic data. Hence it significantly improves the resolution and prediction accuracy.
- values of lithofacies, permeability and porosity used for the training are preferably calculated using a petrophysics methodology, for example by log data directly (porosity) and a definition of lithofacies by the lithology content and petrophysics properties.
- a petrophysics methodology can correspond to a deterministic petrophysics methodology or a stochastic petrophysics methodology.
- a petrophysics methodology corresponds to a deterministic petrophysics methodology.
- Such methodologies are well known to the person skilled in the art, they also include knowhow and geoscientist expertise.
- the ensemble prediction model according to the invention is based at least partly on elastic attributes.
- Elastic attributes are preferably generated from well logs (in particular for the training) or seismic data with standard inversion methods (in particular for the prediction).
- elastic attributes comprise P-wave velocity V p , S-wave velocity V s , density or their mathematic transformations.
- the elastic attributes comprise elastic attributes selected from the group consisting of:
- Lame’s coefficient density * (V p 2 -4/3V s 2 ) - (2/3 * density * V s 2 ).
- the elastic attributes comprise elastic attributes selected from the group consisting of P-lmpedance and V p /V s ratio. Even more preferably, the elastic attributes consist of P-lmpedance and V p /V s ratio.
- a method 100 for predicting values of porosity, permeability and/or lithofacies comprises a step of loading 130 an ensemble prediction model.
- the loaded ensemble prediction model has been trained according to the step of training 120 describe above, more preferably with data collected according to the step of collecting 110 data from well logs.
- the ensemble prediction model includes any one of a stacking, a bagging such as random forest, a gradient boosting such as xgboost or an adaptative boosting such as adaboost.
- the ensemble prediction model includes any one of a random forest or a xgboost.
- the ensemble prediction model have been trained with values of elastic attributes and values of lithofacies, permeability and/or porosity. Such values can be considered as reference values.
- elastic attributes values are classically described as input data, whereas lithofacies, permeability and/or porosity values can be considered as target values.
- those reference values have been obtained from well logs. They are not based on rock sample analysis.
- a method 100 for predicting values of porosity, lithofacies and/or permeability comprises a step of acquiring 170 values of inverted elastic attributes calculated from seismic data generated on a studied carbonate reservoir.
- inverted elastic attributes are elastic attributes obtained from the inversion of seismic data.
- The preferably comprise V p , V s , density or their mathematic transformations.
- the inverted elastic attributes comprise inverted elastic attributes selected from the group consisting of: V p , V s , density, V p /V s ratio, shear modulus, bulk modulus, P- impedance, S-impedance, Poisson’s ratio and Lame’s coefficient.
- values of inverted elastic attributes include: P-lmpedance and V p /V s ratio. More preferably values of inverted elastic attributes consist in: P-lmpedance and V p /V s ratio.
- the inverted elastic attributes have been generated before the implementation of the method according to the invention, in some embodiments, the method according to the invention can include the generation of those inverted elastic attributes.
- a method 100 for predicting values of porosity, lithofacies and/or permeability according to the invention can comprise a step of collecting 140 seismic data.
- the seismic data generated on a carbonate reservoir preferably corresponds to 2D or 3D seismic data.
- the seismic data can be collected via several of transmitters and geophones.
- a method 100 for predicting values of porosity, lithofacies and/or permeability according to the invention can also comprise a step of inversion 150 of seismic data.
- seismic data do not contain direct facies information, it can be used to obtain P-wave velocity, S-wave velocity, density, and other reservoir properties using seismic inversion techniques.
- a method 100 for predicting values of porosity, lithofacies and/or permeability according to the invention can comprise a step of pretreatment 160 of inverted seismic data.
- inverted seismic data may be pretreated 160 so as to facilitate future operations.
- a pretreatment according to the invention may be: normalization of data, re-sampling, aggregation of data, and/or re-coding variables.
- inverted elastic attributes can comprise V p , V s , density or mathematic transformations based on these elastic attributes.
- a method 100 for predicting values of porosity, lithofacies and/or permeability comprises a step of generating 180 predicted values of porosity, lithofacies and/or permeability of the studied carbonate reservoir.
- This step can comprise treating or processing the acquired values of inverted elastic attributes with the ensemble prediction model.
- a method 100 for predicting values of porosity, lithofacies and/or permeability according to the invention can comprise a step of generating 180 predicted values of porosity, lithofacies or permeability of the studied carbonate reservoir.
- a method 100 for predicting values of porosity, lithofacies and/or permeability according to the invention preferably comprises a step of generating 180 predicted values of porosity and lithofacies of the studied carbonate reservoir. More preferably, the method according to the invention comprises a step of generating 180 predicted values of porosity and lithofacies, or permeability of the studied carbonate reservoir.
- a method 100 according to the invention can comprise a step of generating 180 predicted values of:
- the generation 180 of predicted values of porosity and lithofacies is based on a trained ensemble prediction model 181 , inverted elastic attributes 184 and it can also be based on key information 185 value.
- Key information 185 values can for example correspond to reservoirs typology type, a priori geological information, gamma ray values or water-contact values. Such key information without being essential can enhance the accuracy of the prediction.
- seismic data 182 and inverted seismic data 183 can be processed in a method according to the invention, such data are not injected in the trained ensemble prediction model 181 .
- generating 180 predicted values of porosity, lithofacies and/or permeability comprises a generation of predicted values of porosity, lithofacies and/or permeability for each pixel of a 2D model of the carbonate reservoir.
- generating 180 predicted values of porosity, lithofacies and/or permeability comprises a generation of predicted values of porosity, lithofacies and/or permeability for each 3D voxel of a 3D model of the carbonate reservoir.
- the invention generates lithological probability cubes from seismic
- the methods and devices of the present disclosure may facilitate determining how much hydrocarbon is in place in a subterranean carbonate reservoir.
- the methods and devices of the present disclosure may also help for positioning a well position.
- This method based on seismic data and a prediction model trained with well logs, can determine lithofacies collocated with the 3D seismic elastic attribute to build an accurate subsurface reservoir model and this without the need of any well log on this reservoir.
- the ensemble prediction model has been trained preferably from well logs generated from different carbonate reservoirs than the studied carbonate reservoir. Again, this is advantageous because without any well log data, the method according to the invention can produce relevant predicted values of porosity, lithofacies or permeability in similar geological context. However, the accuracy of prediction can be further enhanced.
- the method according to the invention can comprise a process 200 of retraining the trained ensemble prediction model.
- a process 200 of retraining the trained ensemble prediction model will be based on the trained 231 ensemble prediction model or on data used to train the trained ensemble prediction model (e.g. elastic attributes values, porosity values, lithofacies values, permeability values, reservoir typology values).
- data used to train the trained ensemble prediction model e.g. elastic attributes values, porosity values, lithofacies values, permeability values, reservoir typology values.
- those data will be advantageously completed with permeability, lithofacies and/or porosity 232 values and elastic attributes 233 generated from at least one well log 210 of the studied carbonate reservoir.
- well log will be analyzed 220, preferably according to a petrophysics methodology, more preferably a deterministic petrophysics methodology in order to determine lithofacies, permeability and porosity 232 values and elastic attributes 233.
- the value of the key information 234 such as reservoir typology can be advantageously used for the retraining 230 step.
- Key information values 234 can for example correspond to reservoirs typology types, a priori geological information, gamma ray values or water-contact values. Such key information without being essential can enhance the accuracy of the prediction.
- the well log received via a logging tool can be transmitted to a separate system used to perform the workflow steps dedicated to the elastic attributes, lithofacies and porosity calculation.
- the well log data may be input into a petrophysical model stored in a device according to the invention or another device.
- the petrophysical model may incorporate a multi-mineral analysis with classification techniques such as Bayesian classification.
- the method according to the invention can comprise the determination of lithofacies, permeability, porosity and elastic attributes values from a well log conducted on the studied carbonate reservoir and a step of retraining 230 the ensemble prediction model with a set of data comprising determined values of lithofacies, permeability, porosity and elastic attributes.
- one or more well logs of the studied carbonate reservoir can be available before the first implementation of the method according to the invention.
- the present invention uses an ensemble prediction model. As stated, whereas much of the focus of machine learning is on developing the single most accurate prediction model possible for a given task, those ensemble prediction model generate a set of models and then make prediction by aggregating the outputs of these models.
- Such a feature can be advantageously use when performing a retraining of the ensemble prediction model. Indeed, it is possible to give more importance to data generated from the studied carbonate reservoir than data generated from other carbonate reservoir when performing the step of 230 retraining the ensemble prediction model.
- the invention relates to a device 1 for predicting values of porosity, lithofacies and/or permeability in a studied carbonate reservoir based on seismic data.
- a device 1 may include any instrumentality or aggregate of instrumentalities operable to compute, classify, process, transmit, receive, retrieve, originate, switch, store, display, manifest, detect, record, reproduce, handle, or utilize any form of information, intelligence, or data.
- an information handling system may be a personal computer, a network storage device, or any other suitable device and may vary in size, shape, performance, functionality, and price.
- the information handling system may include random access memory (RAM), one or more processing resources such as a central processing unit (CPU) or hardware or software control logic, ROM, and/or other types of nonvolatile memory.
- Additional components of the information handling system may include one or more disk drives, one or more network ports for communication with external devices as well as various input and output (I/O) devices, such as a keyboard, a mouse, and a video display.
- the information handling system may also include one or more buses operable to transmit communications between the various hardware components.
- Figure 5 is a schematic block diagram illustrating various hardware components that may be utilized in the device 1 according to the invention to predict values of porosity, lithofacies and/or permeability in a studied carbonate reservoir based on seismic data.
- the device 1 comprises: one or more memory components 10 configured to store an ensemble prediction model, one or more communication interfaces 20 configured to acquire values of inverted elastic attributes calculated from seismic data; and one or more processors 30 configured to process the acquired values of inverted elastic attributes with the ensemble prediction model to generate predicted values of porosity, lithofacies and/or permeability in the studied carbonate reservoir. It can also comprise one or more user interfaces.
- a device 1 for predicting values of porosity, lithofacies and/or permeability in a studied carbonate reservoir based on seismic data may comprise a memory component 10.
- the memory component 10 may comprise any computer readable medium known in the art including, for example, a volatile memory, such as a static random access memory (SRAM) and a dynamic random access memory (DRAM) , and / or a non-volatile memory, such as read-only memory, flash memories, hard disks, optical disks and magnetic tapes.
- the memory component 10 may include a plurality of instructions or modules or applications for performing various functions.
- the memory component 10 can implement routines, programs, or matrix-type data structures.
- the memory component 10 may comprise a medium readable by a computer system in the form of a volatile memory, such as a random-access memory (RAM) and / or a cache memory.
- the memory component 10, like the other modules, can for example be connected with the other components of the device 1 via a communication bus and one or more data carrier interfaces.
- the memory component is preferably configured to store an ensemble prediction model.
- the ensemble prediction model has been in particular trained with: o values of elastic attributes, said elastic attributes comprising V p , V s , density or their mathematic transformations, and o values of porosity, lithofacies and/or permeability.
- the memory component 10 can be configured to store all data and values generated during the step of collecting 110 data from well logs.
- the memory component 10 can be configured to store all data and values such as reservoir typology values, values of elastic attributes and values of porosity, lithofacies and/or permeability, in particular those used to train the ensemble prediction model.
- the memory component 10 is preferably configured to store instructions capable of implementing the method according to the invention.
- the device 1 can also comprise a communication interface 20.
- the communication interface 20 is preferably configured to transmit data on at least one communication network and may implement a wired or wireless communication.
- the device 1 can communicate with other devices or computer systems and in particular with clients 2 thanks to the communication interface 20.
- communication interfaces 20 enable the device 1 to receive data from the various sensing components on location (e.g., 3D seismic data from geophones).
- the communication is operated via a wireless protocol such as Wi-Fi, 3G, 4G, and/or Bluetooth. These data exchanges may take the form of sending and receiving files.
- the communication interface 20 may be configured to transmit a printable file.
- the communication interface may in particular be configured to allow the communication with a remote terminal, including a client.
- the client is generally any hardware and/or software capable of communication with the device 1 .
- a communication interface 20 according to the invention is, in particular, configured to exchange data with third party devices or systems.
- a communication interface 20 configured to acquire values of inverted elastic attributes calculated from seismic data generated on a studied carbonate reservoir, said inverted elastic attributes comprising V p , V s , density or their mathematic transformations.
- the device 1 may be communicatively coupled to another computing system, such as the client 2 that is configured to execute an inversion modeling software.
- the device 1 may include a communication interface 20 through which another computing system, such as the client 2, sends the inverted elastic attributes as an input into the ensemble prediction model.
- another computing system such as the client 2
- the inversion modeling software may be included in the device 1 according to the invention, or the components of the illustrated control systems may be distributed throughout a greater number of control systems.
- Such devices may all be located at the reservoir site, or at a location remote from the reservoir site.
- a device 1 for predicting values of porosity, lithofacies and/or permeability in a studied carbonate reservoir based on seismic data may comprise a processor 30.
- the processor 30 may be operably coupled to the memory component 10 to execute instructions, encoded in programs, for carrying out the presently disclosed techniques, more particularly to perform the method according to the invention.
- the encoded instructions may be stored in any suitable article of manufacture (such as the memory component 10) that includes at least one tangible non-transitory, computer- readable medium that at least collectively stores these instructions or routines.
- the memory component 10 may contain a set of instructions that, when executed by the processor 30, performs the disclosed method.
- the memory component 10 may include any number of databases or similar storage media that can be queried from the processor 30 as needed to perform the disclosed method.
- the processor 30 is configured to: load the ensemble prediction model and process the acquired values of inverted elastic attributes with the ensemble prediction model to generate predicted values of porosity, lithofacies and/or permeability in the studied carbonate reservoir.
- modules or components are separated in Figure 5, but the invention may provide various types of arrangement, for example a single module cumulating all the functions described here. Similarly, these modules or components may be divided into several electronic boards or gathered on a single electronic board.
- a device 1 according to the invention can be incorporated into a computer system and able to communicate with one or several external devices such as a keyboard, a pointer device, a display, or any device allowing a user to interact with the device 1 .
- the device 1 may also be configured to communicate with or via a human-machine- interface.
- the device 1 can be coupled to a human interface machine (HMI).
- HMI human interface machine
- the HMI may be used to allow the transmission of parameters to the devices or conversely make available to the user the values of the data measured or calculated by the device.
- the HMI is communicatively coupled to a processor and includes a user output interface and a user input interface.
- the user output interface may include an audio and display output interface and various indicators such as visual indicators, audible indicators and haptic indicators.
- the user input interface may include a keyboard, a mouse, or another navigation module such as a touch screen, a touchpad, a stylus input interface, and a microphone for inputting audible signals such as a user speech, data and commands that can be recognized by the processor.
- the user interface may include various input/output devices that enable an operator to, for example, input values of inverted elastic attributes, or visualize an output image of the predicted values of porosity, lithofacies and/or permeability in the studied carbonate reservoir such as a 3D reservoir model on a digital display.
- the invention relates to a non-transitory computer readable medium storing executable instructions which, when executed by a processor, implement a method according to the invention.
- Computer-readable media may include any instrumentality or aggregation of instrumentalities that may retain data and/or instructions for a period of time.
- Computer-readable media may include, for example, without limitation, storage media such as a direct access storage device (e.g., a hard disk drive or floppy disk drive), a sequential access storage device (e.g., a tape disk drive), compact disk, CD-ROM, DVD, RAM, ROM, electrically erasable programmable read-only memory (EEPROM), and/or flash memory; as well as communications media such as wires, optical fibers, microwaves, radio waves, and other electromagnetic and/or optical carriers; and/or any combination of the foregoing.
- storage media such as a direct access storage device (e.g., a hard disk drive or floppy disk drive), a sequential access storage device (e.g., a tape disk drive), compact disk, CD-ROM, DVD, RAM, ROM, electrically erasable programmable read-only memory (EEPROM), and/or flash memory
- the invention relates to a non-transitory computer readable medium storing executable instructions which, when executed by a processor, implement a method for predicting values of porosity, lithofacies and/or permeability in a studied carbonate reservoir based on seismic data, the method comprising:
- an ensemble prediction model said ensemble prediction model having been trained with: o values of elastic attributes, said elastic attributes comprising V p , V s , density or their mathematic transformations, and o values of porosity, lithofacies and/or permeability;
- inverted elastic attributes calculated from seismic data generated on a studied carbonate reservoir, said inverted elastic attributes comprising Vp, Vs, density or their mathematic transformations;
- aspects of the present invention may be embodied as a device, system, method or computer program product. Accordingly, aspects of the present invention may take the form of a fully hardware embodiment, a fully software embodiment (including firmware, resident software, microcode, etc.) or a mode of operation. In addition, aspects of the present invention may take the form of a computer program product incorporated into one or more computer readable media having a computer readable program code embedded therein.
- a computer readable medium may be any tangible medium that may contain, or store, a program for use by or in connection with an instruction execution system, apparatus, or device.
- a computer readable medium may be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared or semiconductor system, apparatus or device, or any suitable combination of the foregoing. More specific examples (a non-exhaustive list) of the computer-readable storage medium would include: a hard disk, a random-access memory (RAM).
- Computer program code for performing operations for aspects of the present invention may be written in any combination of one or more programming languages, including an object- oriented programming language such as Java, C ++, or similar, the programming language "C” or similar programming languages, a scripting language such as Perl, or similar languages, and / or functional languages such as Meta Language.
- Program code can run entirely on a user's computer, partly on a user's computer, and partly on a remote computer or entirely on the computer or remote server. In the latter scenario, the remote computer can be connected to a user's computer by any type of network, including a local area network (LAN) or a wide area network (WAN).
- LAN local area network
- WAN wide area network
- These computer program instructions may be stored on a computer readable medium that can direct a computing device (i.e. computer, server ...), so that the instructions stored in the computer readable medium produce a computing device configured to implement the invention.
- the data used to test the machine learning workflow includes logging data from 7 wells labeled A1 to A7 and 2D inverted seismic data (P-impedance and V p /V s ratio) from presalt carbonate, offshore Brazil.
- Data has been prepared from well logging data.
- the input variable mainly P-impedance and Vp/V s ratio
- target variable i.e. permeability, lithofacies or porosity
- the input variable mainly P-impedance and Vp/V s ratio
- target variable i.e. permeability, lithofacies or porosity
- Machine learning models were used to determine patterns in the relationship between the input variable (P-impedance and V p /V s ratio) and the target variable (lithofacies and porosity) based on the labelled logging data. These patterns were trained to be applied on the seismic inversion results. Whereas seven machine learning models were trained (cf. Table 1), machine learning used in this study were mainly two ensemble learning approaches: Random Forest and XGBoost.
- Machine learning algorithm (MLA) parameter used for the training model of lithofacies prediction for well A6.
- Random Forest and XGBoost perform better with prediction accuracy above 80%.
- DNN Figure 7D
- Random Forest and XGBoost predicted accurately anhydrite with only 33.5% of accuracy while Random Forest and XGBoost predicted it respectively with 70.5% and 68.3% of accuracy
- Domolitic with 34% of accuracy while Random Forest and XGBoost were above 55%
- Silica-rich with 37% of accuracy while Random Forest and XGBoost predicted it respectively with 73.8% and 72.7% of accuracy.
- Naive Bayes Figure 7F predicted domolitic with only 34% of accuracy while Random Forest and XGBoost were above 55%
- Silica-rich with 32% of accuracy while Random Forest and XGBoost predicted it with 73.8% and 72.7% of accuracy respectively.
- Porosity Prediction The essence of porosity prediction is essentially a regression problem, while lithofacies prediction is essentially a classification problem. Hence there is a need to find a prediction model that will perform in both regression problem and classification problem. Therefore, not all machine learning algorithms can be employed to solve this regression problem.
- Ensemble learning methods were compared with four other methods for porosity prediction: SVM, DNN, linear and nonlinear regression. Here we took all the 7 wells into account. Similarly, 50% of the data from each well was randomly extracted to compose the training set which produces a generalized network. Then the generalized network was applied on each well to predict the porosity. We used correlation coefficients as the criteria to measure the performance of different machine learning algorithms.
- Table 3 A comparison of different machine learning algorithms and regression methods for the porosity prediction in terms of overall accuracy and time consumed.
- XGBOOST 120s 0.96 0.96 0.95 0.95 0.95 0.98 0.95 0.95 0.95
- the impact of data abundance on the lithofacies prediction accuracy of the each well have also been evaluated.
- the more data is accounted for classifier training the higher the validation accuracy for lithofacies prediction.
- the two ensemble learning classifiers can well predict the lithofacies (i.e. with 80 % of training data 94,5% for Random Forest and 89,2% for XGBoost).
- the predicted seismic lithofacies distribution at the well location roughly match with the lithofacies profile from wells A1 to A7.
- the silica-rich limestone is also fairly well delineated in the seismic profile.
- the section of silica-rich limestone partially occurring at wells A1 , A2, A4, and A6 can be captured by seismic prediction to a certain degree.
- the present invention describes how to combine, for carbonate reservoir characterization, machine learning with a successful ensemble prediction model, elastic attribute from well log, and inverted seismic attributes.
- ensemble train networks based on well logging data applied to seismic inversion results proved to be efficient in effectively mapping the spatial distribution of lithofacies and porosity.
- ensemble learning shows its potential in lithofacies, permeability and porosity prediction in terms of accuracy and efficiency, in comparison with other machine learning algorithm and industry tools.
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)
Abstract
L'invention concerne un procédé (100) de prédiction de valeurs de porosité, de lithofaciès et/ou de perméabilité dans un réservoir étudié de carbonates en fonction de données sismiques, le procédé consistant : - à charger (130) un modèle de prédiction d'ensemble, ledit modèle de prédiction d'ensemble ayant été formé : * par des valeurs d'attributs élastiques, lesdits attributs élastiques comprenant Vp, Vs, la densité ou leurs transformations mathématiques, et * par des valeurs de porosité, de lithofaciès et/ou de perméabilité ; - à acquérir (170) des valeurs d'attributs élastiques inversés, calculées à partir de données sismiques générées sur un réservoir étudié de carbonates, lesdits attributs élastiques inversés comprenant Vp, Vs, la densité ou leurs transformations mathématiques ; et à générer (180) des valeurs prédites de porosité, de lithofaciès et/ou de perméabilité du réservoir étudié de carbonates, lesdites valeurs prédites se calculant à partir des valeurs acquises d'attributs élastiques inversés et du modèle de prédiction d'ensemble.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/IB2019/001442 WO2021130512A1 (fr) | 2019-12-23 | 2019-12-23 | Dispositif et procédé de prédiction de valeurs de porosité, de lithofaciès et de perméabilité dans un réservoir étudié de carbonates en fonction de données sismiques |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/IB2019/001442 WO2021130512A1 (fr) | 2019-12-23 | 2019-12-23 | Dispositif et procédé de prédiction de valeurs de porosité, de lithofaciès et de perméabilité dans un réservoir étudié de carbonates en fonction de données sismiques |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2021130512A1 true WO2021130512A1 (fr) | 2021-07-01 |
Family
ID=70228375
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/IB2019/001442 WO2021130512A1 (fr) | 2019-12-23 | 2019-12-23 | Dispositif et procédé de prédiction de valeurs de porosité, de lithofaciès et de perméabilité dans un réservoir étudié de carbonates en fonction de données sismiques |
Country Status (1)
Country | Link |
---|---|
WO (1) | WO2021130512A1 (fr) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113933899A (zh) * | 2021-11-10 | 2022-01-14 | 中国石油大学(华东) | 基于地震属性约束指导下的砂砾岩储层建模方法和系统 |
CN113960659A (zh) * | 2021-10-14 | 2022-01-21 | 中国矿业大学 | 一种地震岩石物理驱动的煤层气储层含气量预测方法 |
CN114067214A (zh) * | 2022-01-17 | 2022-02-18 | 北京弘象科技有限公司 | 一种基于多模型融合卷积网络的暴雨识别方法和装置 |
CN114063162A (zh) * | 2021-11-05 | 2022-02-18 | 北京珠玛阳光科技有限公司 | 一种基于小样本机器学习SVM Vp/Vs预测方法 |
CN114167498A (zh) * | 2021-11-30 | 2022-03-11 | 中海石油(中国)有限公司 | 一种测井和地震属性联合的分步聚类分析方法 |
WO2022066447A1 (fr) * | 2020-09-22 | 2022-03-31 | Chevron U.S.A. Inc. | Systèmes et procédés destinés à générer des données de propriété de subsurface en fonction de la position et du temps dans un volume de subsurface à étudier |
CN114624773A (zh) * | 2022-04-29 | 2022-06-14 | 成都理工大学 | 一种确定致密储层含气性的方法及装置 |
CN115576028A (zh) * | 2022-12-01 | 2023-01-06 | 武汉盛华伟业科技股份有限公司 | 基于支持向量机的地质特征层预测方法及系统 |
CN115657132A (zh) * | 2022-12-29 | 2023-01-31 | 成都捷科思石油天然气技术发展有限公司 | 一种基于机器学习的储层预测方法及系统 |
CN116307871A (zh) * | 2023-03-06 | 2023-06-23 | 中国地质调查局油气资源调查中心 | 一种评价页岩微观储集空间有效性的方法 |
WO2023124912A1 (fr) * | 2021-12-31 | 2023-07-06 | 中国石油天然气股份有限公司 | Procédé et appareil de prédiction de catégorie de faciès sédimentaire de roche carbonatée |
US20230258838A1 (en) * | 2022-02-15 | 2023-08-17 | Chevron U.S.A. Inc. | Training data for machine learning seismic inversion |
US11733414B2 (en) | 2020-09-22 | 2023-08-22 | Chevron U.S.A. Inc. | Systems and methods for generating subsurface data as a function of position and time in a subsurface volume of interest |
US20230333277A1 (en) * | 2022-04-18 | 2023-10-19 | Saudi Arabian Oil Company | Lithofacies Guided Core Description Using Unsupervised Machine Learning |
CN117250658A (zh) * | 2023-11-17 | 2023-12-19 | 核工业北京地质研究院 | 建立研究区的地震数据集的方法 |
CN117332908A (zh) * | 2023-12-01 | 2024-01-02 | 长江水利委员会水文局 | 一种耦合集合预报的梯级水库多目标优化调度方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6374185B1 (en) * | 2000-02-18 | 2002-04-16 | Rdsp I, L.P. | Method for generating an estimate of lithological characteristics of a region of the earth's subsurface |
WO2012054487A2 (fr) * | 2010-10-18 | 2012-04-26 | Schlumberger Canada Limited | Estimation de production dans des formations souterraines |
CN109164491A (zh) * | 2018-10-15 | 2019-01-08 | 中国石油大学(北京) | 一种基于分类支持向量机的地震岩相识别方法及系统 |
US20190266501A1 (en) | 2018-02-27 | 2019-08-29 | Cgg Services Sas | System and method for predicting mineralogical, textural, petrophysical and elastic properties at locations without rock samples |
-
2019
- 2019-12-23 WO PCT/IB2019/001442 patent/WO2021130512A1/fr active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6374185B1 (en) * | 2000-02-18 | 2002-04-16 | Rdsp I, L.P. | Method for generating an estimate of lithological characteristics of a region of the earth's subsurface |
WO2012054487A2 (fr) * | 2010-10-18 | 2012-04-26 | Schlumberger Canada Limited | Estimation de production dans des formations souterraines |
US20190266501A1 (en) | 2018-02-27 | 2019-08-29 | Cgg Services Sas | System and method for predicting mineralogical, textural, petrophysical and elastic properties at locations without rock samples |
CN109164491A (zh) * | 2018-10-15 | 2019-01-08 | 中国石油大学(北京) | 一种基于分类支持向量机的地震岩相识别方法及系统 |
Non-Patent Citations (8)
Title |
---|
ABDULRAHMAN MOHAMMAD SALEH ALMOQBEL: "Carbonate Reservoir Characterization Based on Integration of 3-D Seismic Data and Well Logs Using Conventional and Artificial Intelligence Approaches", 14 September 2011 (2011-09-14), XP055730473, Retrieved from the Internet <URL:https://spiral.imperial.ac.uk/handle/10044/1/9516> [retrieved on 20200914] * |
BREIMAN, RANDOM FOREST, 1996 |
BRUNO DE RIBET ET AL: "Machine Learning Provides Higher-Quality Insights into Facies Heterogeneities over Complex Carbonate Reservoirs in a Recently Developed Abu Dhabi Oilfield, Middle East", ABU DHABI INTERNATIONAL PETROLEUM EXHIBITION & CONFERENCE, 1 January 2018 (2018-01-01), XP055730464, DOI: 10.2118/192944-MS * |
CHENGUESTRIN, EXTREME GRADIENT BOOSTING, 2016 |
MORI ET AL.: "Porosity Prediction of a Carbonate Reservoir", CAMPOS BASIN BASED ON THE INTEGRATION OF SEISMIC ATTRIBUTES AND WELL LOG DATA |
PEIPEI LIYURAN ZHANG: "Facies Characterization of a Reservoir in the North Sea Using Machine Learning Techniques", PROJECT FINAL REPORT FOR CS229 |
SILVA ET AL.: "Artificial neural networks to support petrographic classification of carbonate-siliciclastic rocks using well logs and textural information", JOURNAL OF APPLIED GEOPHYSICS, vol. 117, June 2015 (2015-06-01) |
T. MUKERIJ ET AL: "Statistical rock physics: Combining rock physics, information theory, and geostatistics to reduce uncertainty in seismic reservoir characterization", THE LEADING EDGE, 1 March 2001 (2001-03-01), pages 313 - 319, XP055240609, Retrieved from the Internet <URL:http://library.seg.org/doi/pdf/10.1190/1.1438938> [retrieved on 20160112], DOI: 10.1190/1.1438938 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022066447A1 (fr) * | 2020-09-22 | 2022-03-31 | Chevron U.S.A. Inc. | Systèmes et procédés destinés à générer des données de propriété de subsurface en fonction de la position et du temps dans un volume de subsurface à étudier |
US11733414B2 (en) | 2020-09-22 | 2023-08-22 | Chevron U.S.A. Inc. | Systems and methods for generating subsurface data as a function of position and time in a subsurface volume of interest |
CN113960659A (zh) * | 2021-10-14 | 2022-01-21 | 中国矿业大学 | 一种地震岩石物理驱动的煤层气储层含气量预测方法 |
CN114063162B (zh) * | 2021-11-05 | 2024-03-08 | 北京珠玛阳光科技有限公司 | 一种基于小样本机器学习SVM Vp/Vs预测方法 |
CN114063162A (zh) * | 2021-11-05 | 2022-02-18 | 北京珠玛阳光科技有限公司 | 一种基于小样本机器学习SVM Vp/Vs预测方法 |
CN113933899B (zh) * | 2021-11-10 | 2023-07-28 | 中国石油大学(华东) | 基于地震属性约束指导下的砂砾岩储层建模方法和系统 |
CN113933899A (zh) * | 2021-11-10 | 2022-01-14 | 中国石油大学(华东) | 基于地震属性约束指导下的砂砾岩储层建模方法和系统 |
CN114167498B (zh) * | 2021-11-30 | 2023-09-15 | 中海石油(中国)有限公司 | 一种测井和地震属性联合的分步聚类分析方法 |
CN114167498A (zh) * | 2021-11-30 | 2022-03-11 | 中海石油(中国)有限公司 | 一种测井和地震属性联合的分步聚类分析方法 |
WO2023124912A1 (fr) * | 2021-12-31 | 2023-07-06 | 中国石油天然气股份有限公司 | Procédé et appareil de prédiction de catégorie de faciès sédimentaire de roche carbonatée |
CN114067214B (zh) * | 2022-01-17 | 2022-05-13 | 北京弘象科技有限公司 | 一种基于多模型融合卷积网络的暴雨识别方法和装置 |
CN114067214A (zh) * | 2022-01-17 | 2022-02-18 | 北京弘象科技有限公司 | 一种基于多模型融合卷积网络的暴雨识别方法和装置 |
US20230258838A1 (en) * | 2022-02-15 | 2023-08-17 | Chevron U.S.A. Inc. | Training data for machine learning seismic inversion |
US20230333277A1 (en) * | 2022-04-18 | 2023-10-19 | Saudi Arabian Oil Company | Lithofacies Guided Core Description Using Unsupervised Machine Learning |
CN114624773A (zh) * | 2022-04-29 | 2022-06-14 | 成都理工大学 | 一种确定致密储层含气性的方法及装置 |
CN115576028A (zh) * | 2022-12-01 | 2023-01-06 | 武汉盛华伟业科技股份有限公司 | 基于支持向量机的地质特征层预测方法及系统 |
CN115657132A (zh) * | 2022-12-29 | 2023-01-31 | 成都捷科思石油天然气技术发展有限公司 | 一种基于机器学习的储层预测方法及系统 |
CN116307871A (zh) * | 2023-03-06 | 2023-06-23 | 中国地质调查局油气资源调查中心 | 一种评价页岩微观储集空间有效性的方法 |
CN116307871B (zh) * | 2023-03-06 | 2023-10-27 | 中国地质调查局油气资源调查中心 | 一种评价页岩微观储集空间有效性的方法 |
CN117250658A (zh) * | 2023-11-17 | 2023-12-19 | 核工业北京地质研究院 | 建立研究区的地震数据集的方法 |
CN117250658B (zh) * | 2023-11-17 | 2024-02-09 | 核工业北京地质研究院 | 建立研究区的地震数据集的方法 |
CN117332908A (zh) * | 2023-12-01 | 2024-01-02 | 长江水利委员会水文局 | 一种耦合集合预报的梯级水库多目标优化调度方法及系统 |
CN117332908B (zh) * | 2023-12-01 | 2024-02-13 | 长江水利委员会水文局 | 一种耦合集合预报的梯级水库多目标优化调度方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2021130512A1 (fr) | Dispositif et procédé de prédiction de valeurs de porosité, de lithofaciès et de perméabilité dans un réservoir étudié de carbonates en fonction de données sismiques | |
US11668853B2 (en) | Petrophysical inversion with machine learning-based geologic priors | |
US11520077B2 (en) | Automated reservoir modeling using deep generative networks | |
Hall | Facies classification using machine learning | |
CA3067013C (fr) | Modelisation de reservoir basee sur un apprentissage profond | |
US11435498B2 (en) | Subsurface models with uncertainty quantification | |
EP2761331B1 (fr) | Prédiction de propriétés de réservoir avec machine de vecteur de support de moindres carrés | |
Grana et al. | Probabilistic inversion of seismic data for reservoir petrophysical characterization: Review and examples | |
CN111596978A (zh) | 用人工智能进行岩相分类的网页显示方法、模块和系统 | |
Wang et al. | Data-driven S-wave velocity prediction method via a deep-learning-based deep convolutional gated recurrent unit fusion network | |
Wang et al. | An improved method for lithology identification based on a hidden Markov model and random forests | |
Tokpanov et al. | Deep-learning-based automated stratigraphic correlation | |
Ali et al. | Classification of reservoir quality using unsupervised machine learning and cluster analysis: Example from Kadanwari gas field, SE Pakistan | |
US20220178228A1 (en) | Systems and methods for determining grid cell count for reservoir simulation | |
Bhattacharya | A primer on machine learning in subsurface geosciences | |
Chaikine et al. | A new machine learning procedure to generate highly accurate synthetic shear sonic logs in unconventional reservoirs | |
Chen et al. | Prediction of Shear Wave Velocity Based on a Hybrid Network of Two‐Dimensional Convolutional Neural Network and Gated Recurrent Unit | |
Yang et al. | Deep-learning missing well-log prediction via long short-term memory network with attention-period mechanism | |
Bhattacharya | Unsupervised time series clustering, class-based ensemble machine learning, and petrophysical modeling for predicting shear sonic wave slowness in heterogeneous rocks | |
US20240127039A1 (en) | Geologic learning framework | |
Liu | Downscaling seismic data into a geologically sound numerical model | |
Zhao et al. | Rock-physics-guided machine learning for shear sonic log prediction | |
Lee et al. | Basin‐scale prediction of S‐wave Sonic Logs using Machine Learning techniques from conventional logs | |
Ore | A Machine Learning and Data-Driven Prediction and Inversion of Reservoir Brittleness from Geophysical Logs and Seismic Signals: A Case Study in Southwest Pennsylvania, Central Appalachian Basin | |
US20230323760A1 (en) | Prediction of wireline logs using artificial neural networks |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 19872257 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 19872257 Country of ref document: EP Kind code of ref document: A1 |