CN105893674B - The method that geological property prediction is carried out using global covariance - Google Patents
The method that geological property prediction is carried out using global covariance Download PDFInfo
- Publication number
- CN105893674B CN105893674B CN201610196003.6A CN201610196003A CN105893674B CN 105893674 B CN105893674 B CN 105893674B CN 201610196003 A CN201610196003 A CN 201610196003A CN 105893674 B CN105893674 B CN 105893674B
- Authority
- CN
- China
- Prior art keywords
- geological
- data
- global covariance
- global
- geological property
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 83
- 239000011159 matrix material Substances 0.000 claims description 19
- 238000010276 construction Methods 0.000 claims description 4
- 230000035699 permeability Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 7
- 238000011161 development Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 239000003208 petroleum Substances 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
- 238000012916 structural analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Remote Sensing (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Computer Graphics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention proposes a kind of methods for carrying out geological property prediction using global covariance, include the following steps: to establish three-dimensional tectonic model according to the geologic data acquired, wherein, the input data of the three-dimensional tectonic model is the log of oil reservoir to be measured, output data is geological property prediction result, and the three-dimensional tectonic model includes multiple grids；The log of input is roughened into the grid, as the known data point of global covariance method, known to the geological property value of the known data point on the log；According to known data point, the geological property value of the oil reservoir is predicted using global covariance method.This method of the present invention in geological property predictive equation can be made to solve than original gram the speed of golden method for solving have and significantly promoted, solve geological property predictive equation speed and improve 3 times or more.
Description
Technical Field
The invention relates to the technical field of petroleum exploration and development, in particular to a method for predicting geological properties by using global covariance.
Background
Geological prediction is the early work in the petroleum exploration process, and the smooth petroleum exploration and development can be realized only on the basis of obtaining a relatively accurate geological prediction result.
In oil exploration and development, a kriging method is generally adopted in the existing geological property prediction method. Of these, the kriging method is a geostatistical gridding method that is useful in many fields. Kriging attempts to express trends implicit in the data, for example, high points would be connected along a ridge, rather than isolated by bullseye contours. The kriging process involves several factors: change map model, drift type and blockiness. The central idea of the kriging method is to describe the continuity and the anisotropy of data by a variogram, and to calculate the values of weights and covariances by the variogram, which measure the spatial correlation of an attribute at different positions. Is generally applied to any phenomenon that requires estimation of its spatial distribution with point data.
However, when the kriging method is adopted for geological prediction, the problem of poor continuity of a geological prediction result often exists, so that the geological prediction result is not in accordance with geological significance, and the subsequent petroleum exploration and development process is influenced.
Patent CN 104252549a discloses an analysis well placement method based on kriging interpolation, which adopts the kriging method to perform structural analysis on a space field, proposes a variation function model, and finally performs kriging calculation according to the variation graph function model to realize analysis on a well placement area. However, the crikin method is still used for analysis, and the continuity ratio of the analysis result is poor, so that the accuracy of the analysis result is not high.
Disclosure of Invention
The object of the present invention is to solve at least one of the technical drawbacks mentioned.
Therefore, the invention aims to provide a method for predicting the geological property by using the global covariance, which can greatly improve the speed of solving the geological property prediction equation compared with the original kriging solution method and improve the speed of solving the geological property prediction equation by more than 3 times.
In order to achieve the above object, an embodiment of the present invention provides a method for performing geologic attribute prediction by using global covariance, including the following steps:
step S1, establishing a three-dimensional structure model according to the collected geological data, wherein the input data of the three-dimensional structure model is a logging curve of an oil deposit to be detected, the output data is a geological attribute prediction result, and the three-dimensional structure model comprises a plurality of grids;
step S2, coarsening the input logging curve into the grid as the known data point of the global covariance method, wherein the geological attribute value of the known data point on the logging curve is known;
and step S3, according to the known data points in the step S2, predicting the geological attribute value of the oil deposit by using a global covariance method.
Further, in step S1, the building a three-dimensional structure model according to the collected geological data includes the following steps:
interpreting and analyzing the collected geological data to obtain corresponding interpreted geological data, wherein the interpreted geological data comprises: fault data, horizon data and layering data of the detected oil reservoir;
and establishing the three-dimensional construction model according to the interpreted geological data.
Further, the geological property prediction result comprises: and the porosity, saturation and permeability corresponding to the detected oil reservoir.
Further, in the step S3,
firstly, let the geological property prediction equation be as follows:
wherein u isiAnd ujThe geologic attribute values, C (u), for the ith and jth known data points, respectivelyi,uj) Is uiAnd ujCovariance of (d), wjFor weight value, M is the number of the known data points;
then, solving the geological property prediction equation by adopting a global covariance method, and rewriting the geological property prediction equation into the following form:
LUw＝r，
wherein, (LU)ij＝C(ui,uj)，ri＝C(u,ui) L and U are eachA lower triangular matrix and an upper triangular matrix;
then, the LU matrix is solved, and after the L matrix and the U matrix are obtained, the weight w is solvedj，
Finally, knowing the weight wjAnd predicting the geological property value of the unknown data point in the three-dimensional structure model, and recording as a geological property prediction result.
Further, the L and U matrices are:
wherein,
according to the method for predicting the geological property by using the global covariance, the geological property equation is solved by using the global covariance method, when the number of geological property grids is large, the speed of solving the geological property prediction equation is greatly improved compared with the speed of solving the original Krigin solution method, and the speed of solving the geological property prediction equation is improved by more than 3 times. The global covariance advancing method aims at the field of geological attribute prediction, the prediction effect is more geological significance than the original local Krigin, and the method can provide technical support for reservoir exploration and development research.
Additional aspects and advantages of the invention will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the invention.
Drawings
The above and/or additional aspects and advantages of the present invention will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which:
FIG. 1 is a flow chart of a method for geologic property prediction using global covariance, according to an embodiment of the invention;
FIG. 2 is a schematic diagram of a three-dimensional construction model according to an embodiment of the invention;
FIG. 3 is a schematic diagram illustrating the effect of a log upscaling to a three-dimensional formation model according to an embodiment of the present invention;
FIG. 4 is a schematic illustration of a geological property prediction according to an embodiment of the present invention;
FIGS. 5(a) and (b) are schematic diagrams of geological property predictions for the local kriging and global covariance methods, respectively;
fig. 6(a) and (b) are respectively an effect diagram of predicting geological properties by the local kriging method and the global covariance method when the variance is large.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the drawings are illustrative and intended to be illustrative of the invention and are not to be construed as limiting the invention.
The invention provides a method for predicting geological properties by using global covariance. The global covariance method can be understood as a solution method of a geological property prediction equation optimized by combining the geological property prediction field on the basis of a kriging algorithm.
As shown in fig. 1, the method for predicting geologic properties by using global covariance according to the embodiment of the present invention includes the following steps:
and step S1, establishing a three-dimensional structure model according to the collected geological data.
Specifically, the collected geological data is first explained and analyzed to obtain corresponding explained geological data. Wherein the interpreted geological data comprises: fault data, horizon data and stratification data of the measured reservoir. Then, a three-dimensional formation model is built based on the interpreted geological data.
As shown in FIG. 2, the three-dimensional model may represent the geometric features of a subsurface reservoir, including multiple networks. For example, three-dimensional construction models are typically composed of millions or tens of millions of meshes. Specifically, input data of the three-dimensional structure model is a logging curve of the oil reservoir to be tested, and output data is a geological property prediction result.
And step S2, coarsening the input logging curve into the grid of the three-dimensional structure model, wherein the geological property value of the known data point on the logging curve is known as the known data point of the global covariance method.
As shown in fig. 3, the characteristic points on each log are known data points, and the geologic property values of the known data points are known. In practice, often thousands of known data are available in reservoir models, and the number of grids can be as many as the order of millions.
And step S3, according to the known data points in the step S2, predicting the geological attribute value of the oil deposit by using a global covariance method.
Firstly, it is to be noted that the global covariance method optimizes the kriging algorithm as a solution method of the geological property prediction equation based on the kriging algorithm and in combination with the geological property prediction field, so that the algorithm efficiency is improved, and the effect is more significant than that of the original local kriging method in terms of the geological property prediction field. And predicting the geological attribute value of the unknown data point according to the known data point by using a global covariance method.
Step S31, the geological property prediction equation based on the kriging method is to adopt values of linear prediction properties at positions to be estimated of a series of properties, assuming that expected values of the properties at all positions are constants, and taking the variance of the minimized error as a target, so as to solve the weight of data, and through a series of deductions, the geological property prediction equation is set as follows:
wherein u isiAnd ujThe geologic attribute values, C (u), for the ith and jth known data points, respectivelyi,uj) Is uiAnd ujCovariance of (d), wjFor weight values, M is the number of known data points. To obtain M weights wj, an M-dimensional linear system of equations (1) needs to be solved.
The system of equations (1) is rewritten as the following matrix expression:
Aw＝r， (2)
wherein Aij ═ C (u)i,uj)，ri＝C(u,ui)。
For the solution of the geological property prediction equation (1) or (2), in order to better accord with the geological property rule, the global covariance method is adopted for the solution in the step. Since the global covariance method uses all known data to predict the data of the point to be estimated, i.e. when the known point is determined, the matrix a is invariant.
Step S32, solving the geological attribute prediction equation by adopting a global covariance method, and rewriting the geological attribute prediction equation into the following form:
LUw＝r， (3)
wherein, (LU)ij＝C(ui,uj)，ri＝C(u,ui) And L and U are a lower triangular matrix and an upper triangular matrix respectively.
Step S33, solving the LU matrix, and after obtaining the L and U matrices, solving the weight wj。
Specifically, the L and U matrices are:
wherein, U0j＝a0j，(j＝0，1，2，....n-1)，
It should be noted that, the calculation of the L and U matrices is a cyclic iterative process, the U matrix is iterated from row 1 to row n, the L matrix is iterated from column 1 to column n, and the U matrix is calculated before the L matrix in each iteration.
After calculating the L and U matrices according to the above process, the weight w is solved by equation (3)j
In step S34, the weight w is learnedjAnd predicting the geological property value of the unknown data point in the three-dimensional structure model, and recording as a geological property prediction result.
In one embodiment of the invention, the geological property prediction may include: and the porosity, saturation and permeability corresponding to the detected oil reservoir. For example, FIG. 4 shows a predictive value diagram of a reservoir porosity attribute.
In equation (3), compared to O (M) required for conventional solution Aw ═ r3) The operand, global covariance method solution can be in O (M)2) Is completed within the sub-operation. O (M) decomposing LU3) Taking time into consideration, the operation complexity of the geological attribute prediction equation is calculated by O (N.M) by using a global covariance method in a grid containing N unknown data points3) Reduced to O (M)3+N·M2)。
The following describes a comparison between the conventional local kriging method and the global covariance method of the present invention with reference to fig. 5(a) and 5(b), and fig. 6(a) and 6 (b).
Fig. 5(a) and (b) are schematic diagrams of geological property predictions for the local kriging and global covariance methods, respectively. Where fig. 5(a) is generated from local kriging using near 200 known data. As can be seen from fig. 5(a), there are many discontinuities (circles) in the three-dimensional structure model. The three-dimensional structure model shown in fig. 5(b) is obtained by the global covariance method, and the continuity is better than that of the local kriging method.
Fig. 6(a) and (b) are respectively an effect diagram of predicting geological properties by the local kriging method and the global covariance method when the variance is large. As can be seen by comparing the two graphs, the global covariance method adopted in FIG. 6(b) has stronger continuity of geological properties and better conforms to geological rules.
The following describes the comparison of the efficiency of the existing global and local kriging method and the global covariance method of the present invention with reference to tables 1 and 2. Wherein table 1 gives the complexity of the three methods. As can be seen from table 1, the global covariance approach of the present invention is the least complex.
| Global covariance method | Local kriging algorithm | Global kriging algorithm |
| O(NM) | O(NM log M+NK3) | O(NM3) |
TABLE 1
Table 2 gives the algorithm running times of the global and local kriging methods and the global covariance method on three-dimensional grids of different scales.
| Known point number | Size of the grid | Global covariance method | Local kriging algorithm | Global kriging algorithm |
| 281 | 40 × 64 × 8 ═ 2 ten thousand | <1sec | 3sec | 143sec |
| 584 | 80 x 128 x 16 ═ 16 ten thousand | 4sec | 59sec | 2.6hr |
| 1166 | 160 x 256 x 32 ═ 130 ten thousand | 64sec | 1076sec | 142 hr |
| 2315 | 320 × 512 × 64 ═ 1000 ten thousand | 944sec | 5.0hr | / |
TABLE 2
As can be seen from table 2, all the methods run with a linear increase in the number of grids, given that the number of data points remains the same. But under the same conditions, the global covariance method of the invention takes the shortest time. When the number of the three-dimensional grids is more than ten million, the kriging method takes more than 5 hours, the global covariance method only takes 944 seconds, and the algorithm efficiency is greatly improved. Therefore, the geological attribute prediction equation solved by using the global covariance method not only meets the requirement of geological attribute prediction practical problems, but also can greatly improve the solving speed.
According to the method for predicting the geological property by using the global covariance, the geological property equation is solved by using the global covariance method, when the number of geological property grids is large, the speed of solving the geological property prediction equation is greatly improved compared with the speed of solving the original Krigin solution method, and the speed of solving the geological property prediction equation is improved by more than 3 times. The global covariance method provided by the invention aims at the field of geological attribute prediction, the prediction effect is more geological significance than the original local kriging, and a technical support can be provided for reservoir exploration and development research.
In the description herein, references to the description of the term "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
Although embodiments of the present invention have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present invention, and that variations, modifications, substitutions and alterations can be made in the above embodiments by those of ordinary skill in the art without departing from the principle and spirit of the present invention. The scope of the invention is defined by the appended claims and their full range of equivalents.
Claims (4)
1. A method for predicting geological properties by using global covariance is characterized by comprising the following steps:
step S1, establishing a three-dimensional structure model according to the collected geological data, wherein the input data of the three-dimensional structure model is a logging curve of an oil deposit to be detected, the output data is a geological attribute prediction result, and the three-dimensional structure model comprises a plurality of grids;
step S2, coarsening the input logging curve into the grid as the known data point of the global covariance method, wherein the geological attribute value of the known data point on the logging curve is known;
step S3, according to the known data points in step S2, the geological attribute value of the oil reservoir is predicted by using a global covariance method, firstly, a geological attribute prediction equation is set as follows:
wherein u isiAnd ujThe geologic attribute values, C (u), for the ith and jth known data points, respectivelyi,uj) Is uiAnd ujCovariance of (d), wjFor weight value, M is the number of the known data points;
then, solving the geological property prediction equation by adopting a global covariance method, and rewriting the geological property prediction equation into the following form:
LUw＝r，
wherein, (LU)ij＝C(ui,uj)，ri＝C(u,ui) L and U are a lower triangular matrix and an upper triangular matrix respectively;
then, the LU matrix is solved, and after the L matrix and the U matrix are obtained, the weight w is solvedj，
Finally, knowing the weight wjAnd predicting the geological property value of the unknown data point in the three-dimensional structure model, and recording as a geological property prediction result.
2. The method for geologic property prediction using global covariance of claim 1 wherein, in step S1, the step of building a three-dimensional structure model from the collected geologic data comprises the steps of:
interpreting and analyzing the collected geological data to obtain corresponding interpreted geological data, wherein the interpreted geological data comprises: fault data, horizon data and layering data of the detected oil reservoir;
and establishing the three-dimensional construction model according to the interpreted geological data.
3. The method of claim 1, wherein the geologic property prediction result comprises: and the porosity, saturation and permeability corresponding to the detected oil reservoir.
4. The method for geologic property prediction employing global covariance of claim 1 wherein the L and U matrices are:
wherein, U0j＝a0j，(j＝0，1，2，....n-1)，
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610196003.6A CN105893674B (en) | 2016-03-31 | 2016-03-31 | The method that geological property prediction is carried out using global covariance |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610196003.6A CN105893674B (en) | 2016-03-31 | 2016-03-31 | The method that geological property prediction is carried out using global covariance |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105893674A CN105893674A (en) | 2016-08-24 |
| CN105893674B true CN105893674B (en) | 2019-10-25 |
Family
ID=57014225
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610196003.6A Active CN105893674B (en) | 2016-03-31 | 2016-03-31 | The method that geological property prediction is carried out using global covariance |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105893674B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109633781B (en) * | 2018-08-16 | 2021-03-02 | 清能艾科（深圳）能源技术有限公司 | Geological property acquisition method and device, electronic equipment and storage medium |
| CN110532519B (en) * | 2019-09-02 | 2021-06-11 | 中国矿业大学(北京) | Geological data generation method based on chicken flock algorithm and kriging method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101706768A (en) * | 2009-11-03 | 2010-05-12 | 上海第二工业大学 | Method for predicting information on same-position co-kriging by using Markov model |
| CN102693723A (en) * | 2012-04-01 | 2012-09-26 | 北京安慧音通科技有限责任公司 | Method and device for recognizing speaker-independent isolated word based on subspace |
| CN104252549A (en) * | 2013-06-28 | 2014-12-31 | 中国石油天然气股份有限公司 | Analytical well spacing method based on kriging interpolation |
| CN105118091A (en) * | 2015-08-26 | 2015-12-02 | 中国电建集团北京勘测设计研究院有限公司 | Method and system for constructing multi-precision non-uniform geological grid curved surface model |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2506649A (en) * | 2012-10-05 | 2014-04-09 | Bae Systems Plc | Hyperspectral image processing using estimated global covariance and mean |
-
2016
- 2016-03-31 CN CN201610196003.6A patent/CN105893674B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101706768A (en) * | 2009-11-03 | 2010-05-12 | 上海第二工业大学 | Method for predicting information on same-position co-kriging by using Markov model |
| CN102693723A (en) * | 2012-04-01 | 2012-09-26 | 北京安慧音通科技有限责任公司 | Method and device for recognizing speaker-independent isolated word based on subspace |
| CN104252549A (en) * | 2013-06-28 | 2014-12-31 | 中国石油天然气股份有限公司 | Analytical well spacing method based on kriging interpolation |
| CN105118091A (en) * | 2015-08-26 | 2015-12-02 | 中国电建集团北京勘测设计研究院有限公司 | Method and system for constructing multi-precision non-uniform geological grid curved surface model |
Non-Patent Citations (3)
| Title |
|---|
| Oil storage variability in Calanus finmarchicus;CB Miller, et al.;《Ices Journal of Marine Science》;20001231;第57卷(第57期);第1786-1799页 * |
| 协方差矩阵上-下三角分解法在区域土壤水盐条件模拟的应用;徐英，等;《水利学报》;20041130(第11期);第33-38页 * |
| 空间插值技术在地球化学图制作过程中的应用;姚卫华;《中国优秀博硕士学位论文全文数据库 (硕士) 基础科学辑》;20060415(第04期);第A008-79页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105893674A (en) | 2016-08-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Chen et al. | Assisted history matching of channelized models by use of pluri-principal-component analysis | |
| Chai et al. | An efficient method for fractured shale reservoir history matching: The embedded discrete fracture multi-continuum approach | |
| Arnold et al. | Hierarchical benchmark case study for history matching, uncertainty quantification and reservoir characterisation | |
| US20070055447A1 (en) | Method for updating a geological reservoir model by means of dynamic data | |
| US9805144B2 (en) | Method for exploiting a geological reservoir on the basis of a reservoir model matched by means of multiple-scale parameterization | |
| CN101266299B (en) | Method for forecasting oil gas utilizing earthquake data object constructional features | |
| J Al-Mudhafar | Integrating lithofacies and well logging data into smooth generalized additive model for improved permeability estimation: Zubair formation, South Rumaila oil field | |
| CN103364831A (en) | Physical property parameter quantification method based on neural network algorithm | |
| US11513255B2 (en) | Method of characterising a subsurface volume | |
| CN105893674B (en) | The method that geological property prediction is carried out using global covariance | |
| CN114114414A (en) | Artificial intelligence prediction method for 'dessert' information of shale reservoir | |
| CN109655394B (en) | Nuclear magnetic resonance T2 spectrum permeability calculation method under constraint of throat structure parameters | |
| NO20200978A1 (en) | Optimized methodology for automatic history matching of a petroleum reservoir model with ensemble kalman filter | |
| Abdideh | Estimation of permeability using artificial neural networks and regression analysis in an Iran oil field | |
| Masoudi et al. | Estimation of in place hydrocarbon volume in multilayered reservoirs using deterministic and probabilistic approaches | |
| US20200308934A1 (en) | Automatic calibration of forward depositional models | |
| GB2584449A (en) | Apparatus method and computer-program product for calculating a measurable geological metric | |
| Jahanbakhshi et al. | Joint estimation of facies boundaries and petrophysical properties in multi-facies channelized reservoirs through ensemble-based Kalman filter and level set parametrization | |
| CN111381279B (en) | Reservoir porosity quantitative prediction method and device | |
| Ginting et al. | Multiple Markov chains Monte Carlo approach for flow forecasting in porous media | |
| CN110794469B (en) | Gravity inversion method based on minimum geological feature unit constraint | |
| CN115880455A (en) | Three-dimensional intelligent interpolation method based on deep learning | |
| CN112612935A (en) | Complete logging data acquisition method based on self-reasoning model | |
| CN110570102A (en) | Reservoir evaluation method | |
| Cunha | Integrating static and dynamic data for oil and gas reservoir modelling |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 100094, room 4, building 3, No. 401, Feng Xiu Middle Road, Beijing, Haidian District Applicant after: LANDOCEAN ENERGY SERVICES CO.,LTD. Address before: 100094, Room 401, building 4, No. 3, Feng Xiu Middle Road, Beijing, Haidian District 100094, China Applicant before: Hengtai Aipu Oil and Gas Technology Service Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 100094 Room 401, building 4, courtyard 3, fengxiu Middle Road, Haidian District, Beijing Patentee after: New JCM Group Co.,Ltd. Address before: 100094 Room 401, building 4, courtyard 3, fengxiu Middle Road, Haidian District, Beijing Patentee before: LANDOCEAN ENERGY SERVICES CO.,LTD. |