CN102147479A - Modelling method of reservoir space physical property parameters - Google Patents
Modelling method of reservoir space physical property parameters Download PDFInfo
- Publication number
- CN102147479A CN102147479A CN 201110004610 CN201110004610A CN102147479A CN 102147479 A CN102147479 A CN 102147479A CN 201110004610 CN201110004610 CN 201110004610 CN 201110004610 A CN201110004610 A CN 201110004610A CN 102147479 A CN102147479 A CN 102147479A
- Authority
- CN
- China
- Prior art keywords
- data
- space
- phase
- mutually
- grid node
- 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.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to a modelling method of reservoir space physical property parameters, comprising the steps of: 1) digitalizing a deposition phase diagram of a reservoir space to be simulated, and performing integer coding identification, so as to obtain the coordinates of each grid node and the deposition phase identification parameter corresponding to the coordinates; 2) analyzing the correlation between the physical property parameters of the reservoir space to be simulated and the seismic data acting as the constraint data, and regarding the seismic data having the closest correlation with the physical property parameters as 'a second variable' for simulation; 3) performing variation function analysis of the known data by regarding the deposition phase identified in the step 1) and the seismic attributes on the reservoir space obtained in the step 2) as the constraints; 4) performing two-dimensional physical property modelling by regarding the deposition phase identified in the step 1) and the spatial seismic attributes as the constraints; and 5) building a three-dimensional model of the reservoir space physical property parameters according to various two-dimensional physical property models obtained in the step 4).
Description
Technical field
The invention belongs to the geostatistics modeling field of petroleum exploration and development, it particularly carries out the reservoir space physical parameter modeling method of (comprising factor of porosity, permeability, shale index etc.) about several data such as a kind of comprehensive earthquake that utilizes sedimentary facies constraint, well loggings.
Background technology
At present, no matter be abroad or domestic, utilizing geological data and sedimentary facies to carry out constraint modeling is the main flow direction of research at present.Utilize the sedimentary facies constraint to carry out modeling and be also referred to as " two-step approach modeling ".The first step is meant, with the discrete variable modeling method the caused non-average of different sedimentary facies is carried out modeling on a large scale, i.e. " phase modeling "; Second step was meant, and on the basis of first step analog result, promptly under the control of different sedimentary facies, utilized suitable geostatistics algorithm that various physical parameters are carried out modeling, i.e. " phased modeling ".
Particularly, conventional way is at first to adopt the stochastic simulation of based target in the geostatistics, perhaps adopts the multiple spot geostatistics method to generate phase model (or " training image "), carries out the rerum natura simulation then.The method of based target is used to represent the space distribution of geologic feature, and research object often is presumed to simple shape, and as square or spheroid, its position is at random in reservoir, distributes and simulates the geology variability by studying its shape and direction.And the multiple spot geostatistics is used the space structure that its " training image " replaces variation function expression Geological Variable, but the stationarity problem of " training image " is difficult to solve, when comprehensive earthquake information is studied, often earthquake information is made an explanation, be converted into a kind of training image, and then adopt the method for stochastic simulation to generate a plurality of realizations.Because the modeling method of the discrete variable in the stochastic modeling that the process of phase model or generation " training image " of setting up belongs to geostatistics, it is exactly to adopt the given data in the full work area to participate in analog computation that these two kinds of methods have individual common ground, thereby can not be to the single fine description that carries out mutually, so just can not give prominence to the DATA DISTRIBUTION feature of single phase, and can generate a plurality of realizations and select for use for researchists, therefore the rerum natura model that generates on this basis exists uncertainty and multi-solution, and artificial interpretive analysis is based on the realization that algorithm generates, be follow-up, therefore can not guarantee accuracy, will continue this uncertainty with this rerum natura simulation again, thereby make simulate effect not good as second step on basis.
Summary of the invention
At the problems referred to above, the purpose of this invention is to provide a kind of can be to the single fine description that carries out mutually, the modeling method of accuracy height, reservoir space physical parameter that simulate effect is good.
For achieving the above object, the present invention takes following technical scheme: the modeling method of a kind of reservoir space physical parameter, it may further comprise the steps: 1) with the deposition phasor digitizing in the reservoir space that will simulate, carry out the integer coding sign, obtain the coordinate and the sedimentary facies identification parameter corresponding of each grid node with it; 2) analyze the physical parameter in the reservoir space that will simulate and, " second variable " that will simulate with the geological data conduct the most closely of physical parameter correlativity as the correlativity between the geological data of bound data; 3) seismic properties on the reservoir space that obtains the sedimentary facies to identify in the step 1), and step 2) is carried out the variation Functional Analysis of given data as constraint; 4) with the sedimentary facies of step 1) sign and the seismic properties on the space as constraint, carry out two-dimentional rerum natura modeling; 5), set up the three-dimensional model of reservoir space physical parameter according to each the two-dimentional rerum natura model that obtains in the step 4).
The specific implementation process of described step 3) is: when 1. analyzing the variation function of a certain phase, with digitized sedimentary facies data, log data and geological data as the input data; 2. first grid node the deposition phasor after the digitizing in work area begins to scan one by one, and sign is had the coordinate of this phase, and the log data value and the geological data value record that project to behind the grid node get off; 3. log data, the variation function of geological data and their the intersection variation function of this phase of 2. being noted of analytical procedure, promptly obtain the data space distribution characteristics function of this phase, when space variation function is unstable, utilize the space similarity of seismic properties, directly calculate the well weighting coefficient λ of simulation
i:
Wherein, c
iDifference value inverse or the likeness coefficient of representing i mouth well and estimation point seismic properties, i are represented i mouth well, i=1, and 2,3,, n; c
jDifference value inverse or the likeness coefficient of representing j mouth well and estimation point seismic properties, j are represented j mouth well, j=1, and 2,3,, n; 4. went to for the 2. step, circulation finishes until the data of all phases are all analyzed by that analogy.
The specific implementation process of described step 4) is: 1. with digitized sedimentary facies data, geological data, log data as the input data, these three groups of data are begun scanning from first grid node, the geological data that obtains each node and log data be in respectively which mutually in; 2. to a certain when carrying out modeling mutually, begin scanning from first grid node in full work area, there is the coordinate of this phase all to note sign, write down the geological data and the log data of this phase simultaneously; 3. first grid node from this phase of noting begins to calculate, the data that participate in calculating are this interior mutually geological data and log data, until this grid node in mutually calculates and finishes with last, when space variation function is unstable, utilize the space similarity of seismic properties, promptly difference value is directly calculated the well weighting coefficient of simulation; 4. went to for the 2. step, carry out the calculating of next phase, and circulate, finish until all sedimentary facies are all calculated with this.
The concrete modeling method of described step 5) is: 1. deeply concern when setting up that according to earthquake-well logging demarcation utilize layer position, fault information in Depth Domain, figure sets up grid according to sedimentation model, forms the model screen work; If non-mode figure then control as screen work with layer position; 2. model meshes is transformed into time domain, sedimentation model figure or the phasor good according to digitizing determine that the grid of each sedimentary micro distributes, and obtain this interior mutually seismic properties and log data, ask for the weighting coefficient of this interior mutually each well to the simulation points analogue value, thus the physical parameter in the simulation mutually; 3. utilize each corresponding grid, carry out two-dimensional analog by the top by grid, each clathrum adopts the constraint mutually of its correspondence; 4. by determine in the well the time concern deeply, forward Depth Domain to by time domain, set up the three-dimensional oil reservoir model.
The present invention is owing to take above technical scheme, it has the following advantages: the deposition phasor that 1, the present invention studied is the phase distribution plan of determining that is provided by the researchist, rather than a plurality of realizations that draw by modeling algorithm, therefore, can reduce the uncertainty of the deposition phasor that modeling algorithm draws.2, the achievement that is provided by the explanation personnel is provided in the present invention, can determine that on the one hand the phase body of target area distributes, the manpower and the resource that are consumed in the time of can reducing the deposition phasor that draws at the screening modeling algorithm on the other hand.3, the present invention is after adopting constraint mutually, and the sharpness of border of each phase as constraint, is calculating a certain phase time with this border, filters out these interior mutually well point data and participates in calculating, and other well point data just do not participate in having calculated.So both followed explanation personnel's achievement in research, increased determinacy, guaranteed again that the sampling point data that participate in calculating can only comprise the information of this phase, had increased accuracy.4, utilization of the present invention distributes mutually directly to control and participates in rerum natura simulation well or condition data, has the consistance of height thereby make rerum natura distribute with distributing mutually.5, the present invention utilizes the weighting coefficient of the direct computer memory simulation of the space similarity of seismic properties, can be so that the space distribution of analog result and seismic properties be consistent, thus make final model to be consistent with data as much as possible.
Embodiment
Below in conjunction with embodiment the present invention is described in detail.
The present invention is based on following thought: adopt the deposition phasor of geological research and seismic properties to distribute, and the rerum natura simulation is carried out in phase-splitting on this basis, promptly adopting distributes mutually participates in space rerum natura simulation, make analog parameter distribution spatially with distribution is consistent mutually, so not only can guarantee the accuracy of single phase, and also can follow the understanding of researchist on the whole geologic feature.
The inventive method may further comprise the steps:
1) the deposition phasor in the reservoir space that will simulate that researchist in the work area is explained carries out digitizing, and the deposition phasor after the digitizing is carried out integer coding, so that sign; For example identifying the river course is integer 1, valley flat is an integer 2, by that analogy, so that obtain the modeling desired parameters, promptly obtain the coordinate of a certain grid node in the deposition phasor after the digitizing, and with the corresponding sedimentary facies identification parameter of this grid node coordinate, as " first variable ", so just can identify a certain grid node and belong to any sedimentary facies;
2) analyze the physical parameter in the reservoir space that will simulate and as the correlativity between the geological data of bound data, promptly physical parameter and the geological data (seismic properties) by well logging carries out correlation analysis, will with the physical parameter correlativity the most closely geological data as " second variable " of simulating; Wherein, geological data collects for adopting conventional method;
3) seismic properties on the reservoir space sedimentary facies and the step 2 to identify in the step 1)) is carried out the variation Functional Analysis of given data as constraint, and concrete implementation procedure is:
When 1. analyzing the variation function of a certain phase, with digitized sedimentary facies data, the log data that collects and geological data as the input data;
2. first grid node the deposition phasor after the digitizing in work area begins to scan one by one, and sign is had the coordinate of this phase, and the log data value and the geological data value record that project to behind the grid node get off;
3. log data, the variation function of geological data and their the intersection variation function of this phase of 2. noting of analytical procedure can obtain the data space distribution characteristics function of this phase, are used for the virtual space data; When space variation function is unstable, can utilize the space similarity of seismic properties, directly calculate the well weighting coefficient λ of simulation as difference value
i, a conditioned disjunction constraint function of computer memory parameter is used for the virtual space data.Wherein:
In the formula, c
iDifference value inverse or the likeness coefficient of representing i mouth well and estimation point seismic properties, i are represented i mouth well, i=1, and 2,3,, n; c
jDifference value inverse or the likeness coefficient of representing j mouth well and estimation point seismic properties, j are represented j mouth well, j=1, and 2,3,, n.
4. went to for the 2. step, circulation finishes until the data of all phases are all analyzed by that analogy;
4) seismic properties on the space that obtains with the sedimentary facies of step 1) sign, step 2), and the variation function in the step 3) or the likeness coefficient of seismic properties as constraint, are carried out two-dimentional rerum natura modeling, and concrete implementation procedure is:
1. with digitized sedimentary facies data, geological data, log data as the input data, these three groups of data are begun scanning from first grid node, the geological data that obtains each node and log data be in respectively which mutually in;
2. to a certain when carrying out modeling mutually, begin scanning from first grid node in full work area, there is the coordinate of this phase all to note sign, write down the geological data and the log data of this phase simultaneously, so just do not comprised the data of other phases;
3. first grid node from this phase of noting begins to calculate, the data that participate in calculating are this interior mutually geological data and log data, until being calculated, last grid node in mutually finishes, when space variation function is unstable, can utilize the space similarity of seismic properties, directly calculate the well weighting coefficient of simulation as difference value;
4. went to for the 2. step, carry out the calculating of next phase, and circulation according to this, until all being calculated, all sedimentary facies finish;
5) set up the three-dimensional model of reservoir space physical parameter, the method for employing is:
1. demarcate according to earthquake-well logging and concern deeply when setting up that utilize layer position, fault information in Depth Domain, figure sets up grid according to sedimentation model, forms the model screen work; If non-mode figure then control as screen work with layer position;
2. model meshes is transformed into time domain, sedimentation model figure or the phasor good according to digitizing determine that the grid of each sedimentary micro distributes, and obtain this interior mutually seismic properties and log data, ask for the weighting coefficient of this interior mutually each well to the simulation points analogue value, thus the physical parameter in the simulation mutually;
3. utilize each corresponding grid (corresponding a plurality of grids of possibility), carry out two-dimensional analog by the top by grid, each clathrum adopts the constraint mutually of its correspondence;
4. by determine in the well the time concern deeply, forward Depth Domain to by time domain, set up the three-dimensional oil reservoir model.
The various embodiments described above only are used to illustrate the present invention, and every equivalents of carrying out on the basis of technical solution of the present invention and improvement all should not got rid of outside protection scope of the present invention.
Claims (5)
1. the modeling method of a reservoir space physical parameter, it may further comprise the steps:
1) with the deposition phasor digitizing in the reservoir space that will simulate, carry out the integer coding sign, obtain the coordinate and the sedimentary facies identification parameter corresponding of each grid node with it;
2) analyze the physical parameter in the reservoir space that will simulate and, " second variable " that will simulate with the geological data conduct the most closely of physical parameter correlativity as the correlativity between the geological data of bound data;
3) seismic properties on the reservoir space that obtains the sedimentary facies to identify in the step 1), and step 2) is carried out the variation Functional Analysis of given data as constraint;
4) with the sedimentary facies of step 1) sign and the seismic properties on the space as constraint, carry out two-dimentional rerum natura modeling;
5), set up the three-dimensional model of reservoir space physical parameter according to each the two-dimentional rerum natura model that obtains in the step 4).
2. the modeling method of a kind of reservoir as claimed in claim 1 space physical parameter, it is characterized in that: the specific implementation process of described step 3) is:
When 1. analyzing the variation function of a certain phase, with digitized sedimentary facies data, log data and geological data as the input data;
2. first grid node the deposition phasor after the digitizing in work area begins to scan one by one, and sign is had the coordinate of this phase, and the log data value and the geological data value record that project to behind the grid node get off;
3. log data, the variation function of geological data and their the intersection variation function of this phase of 2. being noted of analytical procedure, promptly obtain the data space distribution characteristics function of this phase, when space variation function is unstable, utilize the space similarity of seismic properties, directly calculate the well weighting coefficient λ of simulation
i:
Wherein, c
iDifference value inverse or the likeness coefficient of representing i mouth well and estimation point seismic properties, i are represented i mouth well, i=1, and 2,3,, n; c
jDifference value inverse or the likeness coefficient of representing j mouth well and estimation point seismic properties, j are represented j mouth well, j=1, and 2,3,, n;
4. went to for the 2. step, circulation finishes until the data of all phases are all analyzed by that analogy.
3. the modeling method of a kind of reservoir as claimed in claim 1 space physical parameter, it is characterized in that: the specific implementation process of described step 4) is:
1. with digitized sedimentary facies data, geological data, log data as the input data, these three groups of data are begun scanning from first grid node, the geological data that obtains each node and log data be in respectively which mutually in;
2. to a certain when carrying out modeling mutually, begin scanning from first grid node in full work area, there is the coordinate of this phase all to note sign, write down the geological data and the log data of this phase simultaneously;
3. first grid node from this phase of noting begins to calculate, the data that participate in calculating are this interior mutually geological data and log data, until this grid node in mutually calculates and finishes with last, when space variation function is unstable, utilize the space similarity of seismic properties, promptly difference value is directly calculated the well weighting coefficient of simulation;
4. went to for the 2. step, carry out the calculating of next phase, and circulate, finish until all sedimentary facies are all calculated with this.
4. the modeling method of a kind of reservoir as claimed in claim 2 space physical parameter, it is characterized in that: the specific implementation process of described step 4) is:
1. with digitized sedimentary facies data, geological data, log data as the input data, these three groups of data are begun scanning from first grid node, the geological data that obtains each node and log data be in respectively which mutually in;
2. to a certain when carrying out modeling mutually, begin scanning from first grid node in full work area, there is the coordinate of this phase all to note sign, write down the geological data and the log data of this phase simultaneously;
3. first grid node from this phase of noting begins to calculate, the data that participate in calculating are this interior mutually geological data and log data, until this grid node in mutually calculates and finishes with last, when space variation function is unstable, utilize the space similarity of seismic properties, promptly difference value is directly calculated the well weighting coefficient of simulation;
4. went to for the 2. step, carry out the calculating of next phase, and circulate, finish until all sedimentary facies are all calculated with this.
5. as the modeling method of claim 1 or 2 or 3 or 4 described a kind of reservoir space physical parameters, it is characterized in that: the concrete modeling method of described step 5) is:
1. demarcate according to earthquake-well logging and concern deeply when setting up that utilize layer position, fault information in Depth Domain, figure sets up grid according to sedimentation model, forms the model screen work; If non-mode figure then control as screen work with layer position;
2. model meshes is transformed into time domain, sedimentation model figure or the phasor good according to digitizing determine that the grid of each sedimentary micro distributes, and obtain this interior mutually seismic properties and log data, ask for the weighting coefficient of this interior mutually each well to the simulation points analogue value, thus the physical parameter in the simulation mutually;
3. utilize each corresponding grid, carry out two-dimensional analog by the top by grid, each clathrum adopts the constraint mutually of its correspondence;
4. by determine in the well the time concern deeply, forward Depth Domain to by time domain, set up the three-dimensional oil reservoir model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201110004610 CN102147479B (en) | 2011-01-11 | 2011-01-11 | Modelling method of reservoir space physical property parameters |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201110004610 CN102147479B (en) | 2011-01-11 | 2011-01-11 | Modelling method of reservoir space physical property parameters |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102147479A true CN102147479A (en) | 2011-08-10 |
| CN102147479B CN102147479B (en) | 2013-05-29 |
Family
ID=44421845
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 201110004610 Active CN102147479B (en) | 2011-01-11 | 2011-01-11 | Modelling method of reservoir space physical property parameters |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102147479B (en) |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103617291A (en) * | 2013-12-14 | 2014-03-05 | 中国海洋石油总公司 | Equivalent characterization method for reservoir stratum cause unit interface |
| CN103765245A (en) * | 2011-08-25 | 2014-04-30 | 雪佛龙美国公司 | Hybrid deterministic-geostatistical earth model |
| CN104112285A (en) * | 2013-04-16 | 2014-10-22 | 北京金阳普泰石油技术股份有限公司 | Method of intelligently drawing sedimentary face map facing oilfield exploration and development, and system of intelligently drawing sedimentary face map facing oilfield exploration and development |
| CN104616353A (en) * | 2013-11-05 | 2015-05-13 | 中国石油天然气集团公司 | Modeling for random geologic model of reservoir and preferable method |
| CN104884974A (en) * | 2012-12-05 | 2015-09-02 | 兰德马克绘图国际公司 | Systems and methods for 3d seismic data depth conversion utilizing artificial neural networks |
| CN104991286A (en) * | 2015-07-16 | 2015-10-21 | 西南石油大学 | Sedimentary facies characterization method based on sedimentary modes |
| CN105204090A (en) * | 2015-10-09 | 2015-12-30 | 长江大学 | Boolean simulation method for sandbody with complicated morphology |
| CN105607120A (en) * | 2016-01-19 | 2016-05-25 | 中国海洋石油总公司 | Time-shifting-logging-based method for building initial model with seismic facies constraint |
| CN106522921A (en) * | 2016-11-10 | 2017-03-22 | 中国石油大学(北京) | Dynamic constraint random modeling method and device |
| CN106772587A (en) * | 2017-02-23 | 2017-05-31 | 河海大学 | Seismic elastic parameter Facies Control Modeling method based on same position multiphase collocating kriging |
| CN109072692A (en) * | 2016-05-04 | 2018-12-21 | 沙特阿拉伯石油公司 | Estimated using integrated static bottomhole pressure survey data and the two-dimentional reservoir pressure of simulation modelling |
| CN109425900A (en) * | 2017-09-05 | 2019-03-05 | 中国石油化工股份有限公司 | A kind of Seismic Reservoir Prediction method |
| CN109725348A (en) * | 2017-10-30 | 2019-05-07 | 中国石油化工股份有限公司 | A method of sedimentary facies is identified based on seismic data |
| CN112287532A (en) * | 2020-10-20 | 2021-01-29 | 中海石油(中国)有限公司 | Edge-controlled seismic driving modeling method |
| CN113740908A (en) * | 2020-05-29 | 2021-12-03 | 中国石油化工股份有限公司 | Two-dimensional variation analysis method for seismic slice, electronic device, and medium |
| CN115576007A (en) * | 2022-11-22 | 2023-01-06 | 西南石油大学 | Semi-deterministic fracture modeling method and system based on disorder matrix |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040122634A1 (en) * | 2002-12-19 | 2004-06-24 | Exxonmobil Upstream Research Company | Method of conditioning a random field to have directionally varying anisotropic continuity |
| CN100429528C (en) * | 2005-08-10 | 2008-10-29 | 大庆油田有限责任公司 | Method of deposition phase control for casting sandstone oil reservoir attribute |
-
2011
- 2011-01-11 CN CN 201110004610 patent/CN102147479B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040122634A1 (en) * | 2002-12-19 | 2004-06-24 | Exxonmobil Upstream Research Company | Method of conditioning a random field to have directionally varying anisotropic continuity |
| CN100429528C (en) * | 2005-08-10 | 2008-10-29 | 大庆油田有限责任公司 | Method of deposition phase control for casting sandstone oil reservoir attribute |
Non-Patent Citations (1)
| Title |
|---|
| 《中国地球物理》 20091231 周单等 地震和相约束的建模方法 201 1-5 , * |
Cited By (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103765245A (en) * | 2011-08-25 | 2014-04-30 | 雪佛龙美国公司 | Hybrid deterministic-geostatistical earth model |
| CN103765245B (en) * | 2011-08-25 | 2016-10-12 | 雪佛龙美国公司 | Hybrid Deterministic-Geostatistical Earth Model |
| CN104884974A (en) * | 2012-12-05 | 2015-09-02 | 兰德马克绘图国际公司 | Systems and methods for 3d seismic data depth conversion utilizing artificial neural networks |
| CN104112285A (en) * | 2013-04-16 | 2014-10-22 | 北京金阳普泰石油技术股份有限公司 | Method of intelligently drawing sedimentary face map facing oilfield exploration and development, and system of intelligently drawing sedimentary face map facing oilfield exploration and development |
| CN104616353A (en) * | 2013-11-05 | 2015-05-13 | 中国石油天然气集团公司 | Modeling for random geologic model of reservoir and preferable method |
| CN103617291B (en) * | 2013-12-14 | 2017-01-04 | 中国海洋石油总公司 | Equivalent characterization method for reservoir stratum cause unit interface |
| CN103617291A (en) * | 2013-12-14 | 2014-03-05 | 中国海洋石油总公司 | Equivalent characterization method for reservoir stratum cause unit interface |
| CN104991286A (en) * | 2015-07-16 | 2015-10-21 | 西南石油大学 | Sedimentary facies characterization method based on sedimentary modes |
| CN105204090A (en) * | 2015-10-09 | 2015-12-30 | 长江大学 | Boolean simulation method for sandbody with complicated morphology |
| CN105607120A (en) * | 2016-01-19 | 2016-05-25 | 中国海洋石油总公司 | Time-shifting-logging-based method for building initial model with seismic facies constraint |
| CN109072692A (en) * | 2016-05-04 | 2018-12-21 | 沙特阿拉伯石油公司 | Estimated using integrated static bottomhole pressure survey data and the two-dimentional reservoir pressure of simulation modelling |
| CN109072692B (en) * | 2016-05-04 | 2021-03-26 | 沙特阿拉伯石油公司 | Method for two-dimensional reservoir pressure estimation using integrated static bottom hole pressure survey data and simulation modeling |
| CN106522921A (en) * | 2016-11-10 | 2017-03-22 | 中国石油大学(北京) | Dynamic constraint random modeling method and device |
| CN106522921B (en) * | 2016-11-10 | 2017-10-27 | 中国石油大学(北京) | The stochastic modeling method and device of dynamic constrained |
| CN106772587A (en) * | 2017-02-23 | 2017-05-31 | 河海大学 | Seismic elastic parameter Facies Control Modeling method based on same position multiphase collocating kriging |
| CN109425900A (en) * | 2017-09-05 | 2019-03-05 | 中国石油化工股份有限公司 | A kind of Seismic Reservoir Prediction method |
| CN109725348A (en) * | 2017-10-30 | 2019-05-07 | 中国石油化工股份有限公司 | A method of sedimentary facies is identified based on seismic data |
| CN113740908A (en) * | 2020-05-29 | 2021-12-03 | 中国石油化工股份有限公司 | Two-dimensional variation analysis method for seismic slice, electronic device, and medium |
| CN113740908B (en) * | 2020-05-29 | 2024-05-07 | 中国石油化工股份有限公司 | Two-dimensional variogram analysis method, electronic equipment and medium for seismic slice |
| CN112287532A (en) * | 2020-10-20 | 2021-01-29 | 中海石油(中国)有限公司 | Edge-controlled seismic driving modeling method |
| CN112287532B (en) * | 2020-10-20 | 2024-03-05 | 中海石油(中国)有限公司 | Edge control earthquake driving modeling method |
| CN115576007A (en) * | 2022-11-22 | 2023-01-06 | 西南石油大学 | Semi-deterministic fracture modeling method and system based on disorder matrix |
| CN115576007B (en) * | 2022-11-22 | 2023-03-14 | 西南石油大学 | Semi-deterministic fracture modeling method and system based on disorder matrix |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102147479B (en) | 2013-05-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102147479B (en) | Modelling method of reservoir space physical property parameters | |
| Valentini et al. | Integrating faults and past earthquakes into a probabilistic seismic hazard model for peninsular Italy | |
| Taborda et al. | Large-scale earthquake simulation: computational seismology and complex engineering systems | |
| Lagacherie et al. | Geo-MHYDAS: A landscape discretization tool for distributed hydrological modeling of cultivated areas | |
| WO2017007924A1 (en) | Improved geobody continuity in geological models based on multiple point statistics | |
| CN104297785A (en) | Lithofacies constrained reservoir physical property parameter inversion method and device | |
| CN105701274A (en) | Generation method of three-dimensional local average random field samples of geotechnical parameters | |
| CN105205239A (en) | Method and device for modeling reservoir physical property parameter | |
| CN104360396B (en) | A kind of three kinds of preliminary wave Zoumaling tunnel methods of TTI medium between offshore well | |
| Turrini et al. | Three-dimensional seismo-tectonics in the Po Valley basin, Northern Italy | |
| CN110244363B (en) | Method for predicting fracture-cavity reservoir resource amount | |
| Al-Baldawi | Building a 3D geological model using petrel software for Asmari reservoir, south eastern Iraq | |
| Mojica et al. | Regularization parameter selection in the 3D gravity inversion of the basement relief using GCV: A parallel approach | |
| CN109598068B (en) | Ancient structure constraint modeling method, device and equipment | |
| CN107507179A (en) | Rock And Soil quantitative analysis method based on GOCAD | |
| GB2584449A (en) | Apparatus method and computer-program product for calculating a measurable geological metric | |
| CN103543478A (en) | Geologic morphological interpolation KM (Kriging and Multiple-point geostatistics) method | |
| CN117454710A (en) | Drilling sampling position determining method and modeling method for slope rock-soil body parameter random field simulation | |
| Lutome et al. | 3D geocellular modeling for reservoir characterization of lacustrine turbidite reservoirs: Submember 3 of the third member of the Eocene Shahejie Formation, Dongying depression, Eastern China | |
| GB2584447A (en) | Apparatus method and computer-program product for processing geological data | |
| Burs et al. | Developing a structural and conceptual model of a tectonically limited karst aquifer: a hydrogeological study of the Hastenrather Graben near Aachen, Germany | |
| Cafaro et al. | Landslide hazard assessment and judgment of reliability: a geomechanical approach | |
| Zuur et al. | Spatially continuous data analysis and modelling | |
| CN103901475A (en) | Attribute contour map drawing method and device | |
| Hou et al. | Entropy-based weighting in one-dimensional multiple errors analysis of geological contacts to model geological structure |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C56 | Change in the name or address of the patentee | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Patentee after: China National Offshore Oil Corporation Patentee after: CNOOC Research Institute Patentee after: Beijing Xuriaoyou Energy Technology Co., Ltd. Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Patentee before: China National Offshore Oil Corporation Patentee before: CNOOC Research Center Patentee before: Beijing Xuriaoyou Energy Technology Co., Ltd. |
|
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee after: CNOOC research institute limited liability company Patentee after: China Offshore Oil Group Co., Ltd. Co-patentee after: Beijing Xuriaoyou Energy Technology Co., Ltd. Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee before: CNOOC Research Institute Patentee before: China National Offshore Oil Corporation Co-patentee before: Beijing Xuriaoyou Energy Technology Co., Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20191209 Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee after: CNOOC research institute limited liability company Patentee after: China Offshore Oil Group Co., Ltd. Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee before: CNOOC research institute limited liability company Patentee before: China Offshore Oil Group Co., Ltd. Co-patentee before: Beijing Xuriaoyou Energy Technology Co., Ltd. |