CN108873064A - Establishment method and system of lithofacies probability distribution model - Google Patents
Establishment method and system of lithofacies probability distribution model Download PDFInfo
- Publication number
- CN108873064A CN108873064A CN201810270086.8A CN201810270086A CN108873064A CN 108873064 A CN108873064 A CN 108873064A CN 201810270086 A CN201810270086 A CN 201810270086A CN 108873064 A CN108873064 A CN 108873064A
- Authority
- CN
- China
- Prior art keywords
- data
- petrofacies
- lithofacies
- wave impedance
- value
- 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
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000005070 sampling Methods 0.000 claims abstract description 6
- 238000007788 roughening Methods 0.000 claims description 49
- 241001269238 Data Species 0.000 claims description 12
- 238000000605 extraction Methods 0.000 claims description 9
- 238000010276 construction Methods 0.000 abstract description 4
- 238000012545 processing Methods 0.000 abstract description 3
- 238000005315 distribution function Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 238000006467 substitution reaction Methods 0.000 description 3
- 241001415846 Procellariidae Species 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000008676 import Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 230000001755 vocal effect Effects 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/282—Application of seismic models, synthetic seismograms
-
- 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
-
- 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6226—Impedance
-
- 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/66—Subsurface modeling
- G01V2210/665—Subsurface modeling using geostatistical modeling
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The application provides a method and a system for establishing a lithofacies probability distribution model, wherein the method comprises the following steps: loading the established three-dimensional construction model; importing preset lithofacies data into the three-dimensional structural model, and performing coarsening processing on the preset lithofacies data to obtain coarsened lithofacies data; sampling a preset wave impedance data volume into the three-dimensional construction model to obtain wave impedance data corresponding to the three-dimensional construction model; determining the probability value of the designated lithofacies data corresponding to each wave impedance value according to the wave impedance data and the coarsened lithofacies data so as to obtain probability distribution data of the designated lithofacies data; and respectively loading probability distribution data of each lithofacies data in the coarsened lithofacies data so as to display the probability distribution model of the coarsened lithofacies data in the same histogram window. According to the technical scheme, the lithofacies probability distribution model can be established by combining wave impedance and lithofacies data.
Description
Technical field
This application involves Review of geologic model building techniques field, in particular to the method for building up of a kind of lithofacies probability Distribution Model and it is
System.
Background technique
Lithofacies modeling is of great significance to Reservoir research, and the reasonability of phase model directly affects Reservoir Parameters prediction
Precision.High, the lateral abundant information of seismic data lateral resolution can reduce the uncertainty of phase model.Currently used rock
Phase modeling method is that the probability Distribution Model of lithofacies is established using wave impedance and individual well LITHOFACIES DATA, then in this phase probability distribution
Under the constraint of model, in conjunction with well point data and geological knowledge, lithographic model is established.However there is presently no open source informations to introduce benefit
The method detailed of lithofacies probability Distribution Model is established with wave impedance and LITHOFACIES DATA.
It should be noted that the above description of the technical background be intended merely to it is convenient to the technical solution of the application carry out it is clear,
Complete explanation, and facilitate the understanding of those skilled in the art and illustrate.Cannot merely because these schemes the application's
Background technology part is expounded and thinks that above-mentioned technical proposal is known to those skilled in the art.
Summary of the invention
The method for building up and system for being designed to provide a kind of lithofacies probability Distribution Model of the application embodiment, can
Lithofacies probability Distribution Model is established in conjunction with wave impedance and petrofacies data.
To achieve the above object, on the one hand the application provides the method for building up of a kind of lithofacies probability Distribution Model, the side
Method includes：It is loaded into established three-dimensional tectonic model；Default petrofacies data is imported in the three-dimensional tectonic model, and to described
Default petrofacies data carries out roughening treatment, obtains roughening petrofacies data；It include three kinds of petrofacies datas in the roughening petrofacies data；
By default Acoustic Impedance Data sampler body into the three-dimensional tectonic model, the corresponding wave impedance number of the three-dimensional tectonic model is obtained
According to；According to the Acoustic Impedance Data and the roughening petrofacies data, the corresponding specified petrofacies data of each wave impedance value is determined
Probability value, to obtain the probability distribution data of the specified petrofacies data；It is loaded into respectively each in the roughening petrofacies data
The probability distribution data of petrofacies data, to show the probability distribution mould of the roughening petrofacies data in the same histogram window
Type.
Further, after being loaded into established three-dimensional tectonic model, the method also includes：By the three-dimensional construction
All layers in model are set as identical web thickness.
Further, it is determined that the probability value that each wave impedance is worth corresponding specified petrofacies data includes：Construct the roughening
The histogram of Acoustic Impedance Data corresponding to petrofacies data is specified in petrofacies data；Extract the Acoustic Impedance Data in the histogram
Right, the Acoustic Impedance Data centering includes wave impedance value and corresponding number of grid；According to the Acoustic Impedance Data pair of extraction,
The probability value that each wave impedance is worth corresponding specified petrofacies data is calculated, the probability distribution number of the specified petrofacies data is obtained
According to.
Further, the histogram of Acoustic Impedance Data corresponding to petrofacies data is specified in constructing the roughening petrofacies data
When figure, the method also includes：Numerical filters are set for the histogram, and the numerical filters are applied to the finger
Determine in petrofacies data.
Further, it calculates each wave impedance and is worth the probability value of corresponding specified petrofacies data and include：It calculates described specified
Ratio of the number of grid of each wave impedance value corresponding to petrofacies data in the specified total number of grid of wave impedance value, and by institute
Ratio is stated as the corresponding probability value of each wave impedance value.
Further, the method also includes：When the corresponding wave impedance value of the specified petrofacies data is less than or equal to
When first threshold, first constant is set by the probability value that the wave impedance is worth the corresponding specified petrofacies data；When described
It is when the specified corresponding wave impedance value of petrofacies data is more than or equal to second threshold, wave impedance value is corresponding described specified
The probability value of petrofacies data is set as second constant.
Further, when wave impedance value is less than or equal to first threshold, 1 is set by the probability value of the first lithofacies,
0 is set by the second lithofacies and third lithofacies 0 probability value；When wave impedance value is more than or equal to second threshold, by the
The probability value of one lithofacies 0 is set as 0, sets 0 for the probability value of the second lithofacies, sets 1 for the probability value in third lithofacies.
To achieve the above object, what on the other hand the application also provided a kind of lithofacies probability Distribution Model establishes system, institute
The system of stating includes：Model is loaded into unit, for being loaded into established three-dimensional tectonic model；Roughening treatment unit, for described
Default petrofacies data is imported in three-dimensional tectonic model, and roughening treatment is carried out to the default petrofacies data, obtains roughening lithofacies
Data；It include three kinds of petrofacies datas in the roughening petrofacies data；Acoustic Impedance Data sampling unit is used for default wave impedance number
According to sampler body into the three-dimensional tectonic model, the corresponding Acoustic Impedance Data of the three-dimensional tectonic model is obtained；Probability distribution number
According to determination unit, for determining that each wave impedance value is corresponding according to the Acoustic Impedance Data and the roughening petrofacies data
The probability value of specified petrofacies data, to obtain the probability distribution data of the specified petrofacies data；Probability Distribution Model shows single
Member, for being loaded into the probability distribution data of each petrofacies data in the roughening petrofacies data respectively, in the same histogram
The probability Distribution Model of the roughening petrofacies data is shown in window.
Further, the probability distribution data determination unit includes：Histogram constructs module, for constructing the roughening
The histogram of Acoustic Impedance Data corresponding to petrofacies data is specified in petrofacies data；Data are described for extracting to extraction module
Acoustic Impedance Data pair in histogram, the Acoustic Impedance Data centering include wave impedance value and corresponding number of grid；Probability value
Computing module calculates each wave impedance and is worth corresponding specified petrofacies data for the Acoustic Impedance Data pair according to extraction
Probability value obtains the probability distribution data of the specified petrofacies data.
Further, the system also includes：First setup module, for being less than or equal to the first threshold when wave impedance value
When value, first constant is set by the probability value that the wave impedance is worth corresponding specified lithofacies；Second setup module, for working as wave
When impedance value is more than or equal to second threshold, second is set often for the probability value that the wave impedance is worth corresponding specified lithofacies
Number.
Therefore petrofacies data and Acoustic Impedance Data by being loaded onto the three-dimensional tectonic model pre-established by the application
In, so as to establish lithofacies probability Distribution Model in conjunction with petrofacies data and Acoustic Impedance Data.Specifically, it can determine specified
Wave impedance is worth corresponding various lithofacies probability values, to obtain the probability distribution data that each wave impedance is worth corresponding lithofacies.
Then it is loaded into the probability distribution data of each petrofacies data respectively.Technical solution provided by the present application can rapidly and accurately be built
Vertical lithofacies probability Distribution Model.
Referring to following description and accompanying drawings, specific implementations of the present application are disclosed in detail, specify the original of the application
Reason can be in a manner of adopted.It should be understood that presently filed embodiment is not so limited in range.In appended power
In the range of the spirit and terms that benefit requires, presently filed embodiment includes many changes, modifications and is equal.
The feature for describing and/or showing for a kind of embodiment can be in a manner of same or similar one or more
It uses in a other embodiment, is combined with the feature in other embodiment, or the feature in substitution other embodiment.
It should be emphasized that term "comprises/comprising" refers to the presence of feature, one integral piece, step or component when using herein, but simultaneously
It is not excluded for the presence or additional of one or more other features, one integral piece, step or component.
Detailed description of the invention
Included attached drawing is used to provide to be further understood from the application embodiment, and which constitute the one of specification
The principle of the application for illustrating presently filed embodiment, and with verbal description is come together to illustrate in part.It should be evident that
The accompanying drawings in the following description is only some embodiments of the application, for those of ordinary skill in the art, is not being paid
Out under the premise of creative labor, it is also possible to obtain other drawings based on these drawings.In the accompanying drawings：
Fig. 1 is the method for building up flow chart of lithofacies probability Distribution Model in the application embodiment；
Fig. 2 is the schematic diagram of probability distribution data in the application embodiment；
Fig. 3 is the schematic diagram of probability Distribution Model in the application embodiment；
Fig. 4 is that lithofacies probability Distribution Models establishes system structure diagram in the application embodiment.
Specific embodiment
In order to make those skilled in the art better understand the technical solutions in the application, below in conjunction with the application reality
The attached drawing in mode is applied, the technical solution in the application embodiment is clearly and completely described, it is clear that described
Embodiment is only a part of embodiment of the application, rather than whole embodiments.Based on the embodiment party in the application
Formula, all other embodiment obtained by those of ordinary skill in the art without making creative efforts, is all answered
When the range for belonging to the application protection.
Referring to Fig. 1, a kind of method for building up of lithofacies probability Distribution Model provided by the present application, includes the following steps.
S1：It is loaded into established three-dimensional tectonic model.
S2：Default petrofacies data is imported in the three-dimensional tectonic model, and the default petrofacies data is roughened
Processing obtains roughening petrofacies data；It include three kinds of petrofacies datas in the roughening petrofacies data.
S3：By default Acoustic Impedance Data sampler body into the three-dimensional tectonic model, the three-dimensional tectonic model pair is obtained
The Acoustic Impedance Data answered.
S4：According to the Acoustic Impedance Data and the roughening petrofacies data, determine that each wave impedance value is corresponding specified
The probability value of petrofacies data, to obtain the probability distribution data of the specified petrofacies data.
S5：It is loaded into the probability distribution data of the roughening petrofacies data and corresponding each petrofacies data respectively, same
The probability Distribution Model of the roughening petrofacies data is shown in one histogram window.
In the present embodiment, above-mentioned steps can execute in software " Petrel ".It specifically, can in " Petrel "
To load the three-dimensional tectonic model having had built up, which can embody the structure on purpose stratum.
It in the present embodiment, can will be in the three-dimensional tectonic model after being loaded with the three-dimensional tectonic model
All layers be set as identical web thickness.It specifically, can be by " Layering " control in software, by all layers
Web thickness is set as 5 meters.
In the present embodiment, after being loaded with three-dimensional tectonic model, can continue to import the lithofacies having had been prepared for
Data.The petrofacies data can carry out roughening treatment, to obtain roughening petrofacies data.It is had been prepared for further, it is also possible to import
Good Wave Impedance Data Volume.The Wave Impedance Data Volume can carry out sampling setting.Specifically, it can use earthquake weight in software
The mode for sampling (Seismic resampling) is obtained by the Acoustic Impedance Data sampler body into the three-dimensional tectonic model
The corresponding Acoustic Impedance Data of the three-dimensional tectonic model.
In the present embodiment, roughening petrofacies data may include these three petrofacies datas of high, inter, low.Locating
During reason, these three petrofacies datas can be handled respectively.By taking inter petrofacies data as an example, it can beat in software
Histogram window is opened, the corresponding relationship of Acoustic Impedance Data and number of grid can be shown in the histogram window.In practical application
In, numerical filters (value filter) can be set for the histogram, and by institute in the attribute column of histogram window
It states numerical filters to be applied in specified petrofacies data (inter), specifies rock in the roughening petrofacies data so as to construct
Histogram between phase data and the Acoustic Impedance Data.In the histogram, abscissa indicates Acoustic Impedance Data, ordinate table
Show the number of grid of specified petrofacies data.
In the present embodiment, after constructing the histogram, it can be based on the histogram, generate the specified rock
The corresponding wave impedance distribution function of phase data.Specifically, (can be divided by " the Distribution function " in software
Cloth function) function button, generate the corresponding wave impedance distribution function of the specified petrofacies data.
In the present embodiment, data pair, the data centering packet can be extracted from the wave impedance distribution function
Include wave impedance value and corresponding number of grid.Specifically, it can be popped up by " Spreadsheet " function button in software
The data pair of the wave impedance distribution function.It is then possible to by these data to copying in excel table.Wherein, a column can be with
Indicate the corresponding wave impedance value of inter petrofacies data, another column can indicate corresponding number of grid.
In the present embodiment, by repeating above-mentioned step, can successively extract high, inter and low this three
The data pair of kind of petrofacies data, these data are to may finally be summarized in excel table.According to the data pair of extraction,
So as to calculate the probability value that each wave impedance is worth corresponding various petrofacies datas, to obtain the probability of various petrofacies datas
Distributed data.Specifically, the number of grid of each wave impedance value corresponding to the specified petrofacies data can be calculated specified
Ratio in the total number of grid of wave impedance value, and it is worth corresponding probability value for the ratio as each wave impedance.For example,
Wave impedance be 1 when, lithofacies high has 1, and lithofacies inter has 1, and lithofacies low has 1, then respective probability be three/
One, when wave impedance is 2, lithofacies high has 1, and lithofacies inter has 2, and lithofacies low has 3, then respective probability is six points
One of, one third and half.Six column number evidences are shared in Fig. 2, every column number contains two groups of data in, wherein preceding
Three columns have separately included the corresponding wave impedance value (AI) of different petrofacies datas and number of grid (number), and then three columns separately include
The corresponding wave impedance value (AI) of different petrofacies datas and probability.
In the present embodiment, for different size of wave impedance value, some probability values artificially specified can be set.Tool
Body, when the wave impedance value is less than or equal to first threshold, the wave impedance is worth corresponding specified petrofacies data
Probability value is set as first constant；It is when the wave impedance value is more than or equal to second threshold, the wave impedance value is corresponding
The probability value of specified petrofacies data be set as second constant；Wherein, the first threshold is less than the second threshold.It please refers to
Fig. 2, for example, when wave impedance value is less than or equal to 12520,1 is set by the probability value of the first lithofacies (high), it will
The probability value of second lithofacies (inter) and third lithofacies (low) is set as 0；When wave impedance value is more than or equal to 18040,
0 is set by the probability value of the first lithofacies, sets 0 for the probability value of the second lithofacies, the probability value in third lithofacies is arranged
It is 1.
In this way, every kind of petrofacies data can obtain respective probability distribution data by above-mentioned processing.Actually answering
In, the corresponding probability distribution data of every kind of petrofacies data can respectively save as text file, then, it is loaded into respectively described thick
Change the probability distribution data of each petrofacies data in petrofacies data, it is described thick so as to be shown in the same histogram window
Change the corresponding probability Distribution Model of petrofacies data.Specifically, referring to Fig. 3, three text files can be loaded in software,
Then in histogram window, the probability Distribution Model being made of three kinds of curves can be shown.
Referring to Fig. 4, the application a kind of lithofacies probability Distribution Model is also provided establish system, the system comprises：
Model is loaded into unit 100, for being loaded into established three-dimensional tectonic model；
Roughening treatment unit 200, for importing default petrofacies data in the three-dimensional tectonic model, and to described default
Petrofacies data carries out roughening treatment, obtains roughening petrofacies data；It include three kinds of petrofacies datas in the roughening petrofacies data；
Acoustic Impedance Data sampling unit 300 is used for by default Acoustic Impedance Data sampler body into the three-dimensional tectonic model,
Obtain the corresponding Acoustic Impedance Data of the three-dimensional tectonic model；
Probability distribution data determination unit 400 is used for according to the Acoustic Impedance Data and the roughening petrofacies data, really
Fixed each wave impedance is worth the probability value of corresponding specified petrofacies data, to obtain the probability distribution number of the specified petrofacies data
According to；
Probability Distribution Model display unit 500, for being loaded into the probability distribution data of the roughening petrofacies data respectively, with
The probability Distribution Model of the roughening petrofacies data is shown in the same histogram window.
In the present embodiment, the probability distribution data determination unit includes：
Histogram constructs module, and wave impedance number corresponding to petrofacies data is specified in the roughening petrofacies data for constructing
According to histogram；
Data are to extraction module, for extracting the Acoustic Impedance Data pair in the histogram, the Acoustic Impedance Data centering
Including wave impedance value and corresponding number of grid；
It is corresponding to calculate each wave impedance value for the Acoustic Impedance Data pair according to extraction for probability value computing module
The probability value of specified petrofacies data, obtains the probability distribution data of the specified petrofacies data.
In the present embodiment, the system also includes：
First setup module, for when wave impedance value is less than or equal to first threshold, the wave impedance value to be corresponded to
The probability values of specified lithofacies be set as first constant.
Second setup module, for when wave impedance value is more than or equal to second threshold, the wave impedance value to be corresponded to
The probability values of specified lithofacies be set as second constant.
Therefore petrofacies data and Acoustic Impedance Data by being loaded onto the three-dimensional tectonic model pre-established by the application
In, so as to establish lithofacies probability Distribution Model in conjunction with petrofacies data and Acoustic Impedance Data.Specifically, it can determine specified
Wave impedance is worth corresponding various lithofacies probability values, to obtain the probability distribution data that each wave impedance is worth corresponding lithofacies.
Then it is loaded into the probability distribution data of each petrofacies data respectively.Technical solution provided by the present application can rapidly and accurately be built
Vertical lithofacies probability Distribution Model.
Those skilled in the art are supplied to the purpose described to the description of the various embodiments of the application above.It is not
It is intended to exhaustion or be not intended to and limit the invention to single disclosed embodiment.As described above, the application's is various
Substitution and variation will be apparent for above-mentioned technology one of ordinary skill in the art.Therefore, although specifically begging for
Some alternative embodiments are discussed, but other embodiment will be apparent or those skilled in the art are opposite
It is easy to obtain.The application is intended to include all substitutions of the invention discussed herein, modification and variation, and falls in
Other embodiment in the spirit and scope of above-mentioned application.
Each embodiment in this specification is described in a progressive manner, same and similar between each embodiment
Part may refer to each other, what each embodiment stressed is the difference with other embodiments.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that there are many deformations by the application
With variation without departing from spirit herein, it is desirable to which the attached claims include these deformations and change without departing from the application
Spirit.
Claims (9)
1. a kind of method for building up of lithofacies probability Distribution Model, which is characterized in that the method includes：
It is loaded into established three-dimensional tectonic model；
Default petrofacies data is imported in the three-dimensional tectonic model, and roughening treatment is carried out to the default petrofacies data, is obtained
To roughening petrofacies data；It include three kinds of petrofacies datas in the roughening petrofacies data；
By default Acoustic Impedance Data sampler body into the three-dimensional tectonic model, the corresponding wave resistance of the three-dimensional tectonic model is obtained
Anti- data；
According to the Acoustic Impedance Data and the roughening petrofacies data, the corresponding specified petrofacies data of each wave impedance value is determined
Probability value, to obtain the probability distribution data of the specified petrofacies data；
It is loaded into the probability distribution data of each petrofacies data in the roughening petrofacies data respectively, in the same histogram window
The middle probability Distribution Model for showing the roughening petrofacies data.
2. described the method according to claim 1, wherein after being loaded into established three-dimensional tectonic model
Method further includes：
Identical web thickness is set by all layers in the three-dimensional tectonic model.
3. the method according to claim 1, wherein determining that each wave impedance is worth corresponding specified petrofacies data
Probability value includes：
Construct the histogram that Acoustic Impedance Data corresponding to petrofacies data is specified in the roughening petrofacies data；
The Acoustic Impedance Data pair in the histogram is extracted, the Acoustic Impedance Data centering includes wave impedance value and corresponding grid
Quantity；
According to the Acoustic Impedance Data pair of extraction, the probability value that each wave impedance is worth corresponding specified petrofacies data is calculated, is obtained
To the probability distribution data of the specified petrofacies data.
4. according to the method described in claim 3, it is characterized in that, specifying petrofacies data in constructing the roughening petrofacies data
When the histogram of corresponding Acoustic Impedance Data, the method also includes：
Numerical filters are set for the histogram, and the numerical filters are applied in the specified petrofacies data.
5. according to the method described in claim 3, it is characterized in that, calculating each wave impedance is worth corresponding specified petrofacies data
Probability value includes：
The number of grid of each wave impedance value corresponding to the specified petrofacies data is calculated in the total grid number of specified wave impedance value
Ratio in amount, and it is worth corresponding probability value for the ratio as each wave impedance.
6. according to the method described in claim 5, it is characterized in that, the method also includes：
When the wave impedance value in the specified petrofacies data is less than or equal to first threshold, the wave impedance is worth corresponding
The probability value of first lithofacies is set as 1,0 is set by the probability value that the wave impedance is worth corresponding second lithofacies, by the wave
The probability value of the corresponding third lithofacies of impedance value is set as 0；
When the wave impedance value in the specified petrofacies data is more than or equal to second threshold, the wave impedance is worth corresponding
The probability value of third lithofacies is set as 1 and sets 0 for the probability value that the wave impedance is worth corresponding first lithofacies, by the wave resistance
The probability value of corresponding second lithofacies of anti-value is set as 0.
7. a kind of lithofacies probability Distribution Model establishes system, which is characterized in that the system comprises：
Model is loaded into unit, for being loaded into established three-dimensional tectonic model；
Roughening treatment unit, for importing default petrofacies data in the three-dimensional tectonic model, and to the default lithofacies number
According to roughening treatment is carried out, roughening petrofacies data is obtained；It include three kinds of petrofacies datas in the roughening petrofacies data；
Acoustic Impedance Data sampling unit, for default Acoustic Impedance Data sampler body into the three-dimensional tectonic model, to be obtained institute
State the corresponding Acoustic Impedance Data of three-dimensional tectonic model；
Probability distribution data determination unit, for determining each according to the Acoustic Impedance Data and the roughening petrofacies data
Wave impedance is worth the probability value of corresponding specified petrofacies data, to obtain the probability distribution data of the specified petrofacies data；
Probability Distribution Model display unit, for being loaded into the probability distribution of each petrofacies data in the roughening petrofacies data respectively
Data, to show the probability Distribution Model of the roughening petrofacies data in the same histogram window.
8. system according to claim 7, which is characterized in that the probability distribution data determination unit includes：
Histogram constructs module, specifies Acoustic Impedance Data corresponding to petrofacies data for constructing in the roughening petrofacies data
Histogram；
Data are to extraction module, and for extracting the Acoustic Impedance Data pair in the histogram, the Acoustic Impedance Data centering includes
Wave impedance value and corresponding number of grid；
It is corresponding specified to calculate each wave impedance value for the Acoustic Impedance Data pair according to extraction for probability value computing module
The probability value of petrofacies data obtains the probability distribution data of the specified petrofacies data.
9. system according to claim 8, which is characterized in that the system also includes：
First setup module, for when wave impedance value is less than or equal to first threshold, the wave impedance to be worth corresponding finger
The probability value for determining lithofacies is set as first constant；
Second setup module, for when wave impedance value is more than or equal to second threshold, the wave impedance to be worth corresponding finger
The probability value for determining lithofacies is set as second constant.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810270086.8A CN108873064B (en) | 2018-03-29 | 2018-03-29 | Establishment method and system of lithofacies probability distribution model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810270086.8A CN108873064B (en) | 2018-03-29 | 2018-03-29 | Establishment method and system of lithofacies probability distribution model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108873064A true CN108873064A (en) | 2018-11-23 |
| CN108873064B CN108873064B (en) | 2020-06-09 |
Family
ID=64326410
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810270086.8A Active CN108873064B (en) | 2018-03-29 | 2018-03-29 | Establishment method and system of lithofacies probability distribution model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108873064B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105301647A (en) * | 2014-06-10 | 2016-02-03 | 中国石油化工股份有限公司 | Method for distinguishing calcareous mudstone from sandstone |
| CN105738952A (en) * | 2016-02-02 | 2016-07-06 | 中国石油大学(北京) | Horizontal well region reservoir rock facies modeling method |
| CN106066493A (en) * | 2016-05-24 | 2016-11-02 | 中国石油大学(北京) | Bayes's petrofacies method of discrimination and device |
| EP3121622A1 (en) * | 2015-07-24 | 2017-01-25 | Bergen Teknologioverforing AS | Method of predicting parameters of a geological formation |
-
2018
- 2018-03-29 CN CN201810270086.8A patent/CN108873064B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105301647A (en) * | 2014-06-10 | 2016-02-03 | 中国石油化工股份有限公司 | Method for distinguishing calcareous mudstone from sandstone |
| EP3121622A1 (en) * | 2015-07-24 | 2017-01-25 | Bergen Teknologioverforing AS | Method of predicting parameters of a geological formation |
| CN105738952A (en) * | 2016-02-02 | 2016-07-06 | 中国石油大学(北京) | Horizontal well region reservoir rock facies modeling method |
| CN106066493A (en) * | 2016-05-24 | 2016-11-02 | 中国石油大学(北京) | Bayes's petrofacies method of discrimination and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108873064B (en) | 2020-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104750896B (en) | A kind of fractured-cavernous carbonate reservoir method for numerical simulation | |
| CN105913494B (en) | The modeling of multi-scale facture fine geology and method for numerical simulation and device | |
| Mallison et al. | Unstructured cut-cell grids for modeling complex reservoirs | |
| CN106227957A (en) | The method of equivalent fissure modeling | |
| WO2018093614A1 (en) | Parallel reservoir simulation with accelerated aquifer calculation | |
| CN105445785A (en) | Transverse wave offset velocity modeling method | |
| CN105205239A (en) | Method and device for modeling reservoir physical property parameter | |
| CN105549084A (en) | Method and system for building three-dimensional high-precision velocity model | |
| CN103984802A (en) | Finite element modeling simplification method of honeycomb sandwich structure | |
| CN106886046A (en) | Method and device for determining available reserves of unproductive blocks of fracture-cavity gas reservoir | |
| CN103969682A (en) | Drilling and logging data depth matching method and system | |
| CN104933246A (en) | Method for quantitatively representing contribution of gas field reserves parameters to reserves change | |
| CN114238488A (en) | Multi-professional digital fruit integration and data interaction method and system | |
| CN114547953B (en) | Fracturing construction parameter optimization method and system based on optimization design chart | |
| CN113552621A (en) | Shale gas ground stress determination method and device | |
| CN110390154A (en) | A method of improving Complex reservoir reservoir numerical simulation efficiency | |
| CN109558614B (en) | Simulation method and system for gas flow in shale gas reservoir multi-scale fracture | |
| CN110863818A (en) | Description method and device for residual oil/gas distribution | |
| CN111861129B (en) | Multi-fidelity injection-production optimization method based on multi-scale oil reservoir model | |
| CN108873064A (en) | Establishment method and system of lithofacies probability distribution model | |
| CN106662665A (en) | Re-ordered interpolation and convolution for faster staggered-grid processing | |
| CN106201655A (en) | Virtual machine distribution method and virtual machine distribution system | |
| CN112419493B (en) | Shale reservoir three-dimensional attribute model building method and device | |
| CN112507551B (en) | Unstructured dynamic mesh generation method and device | |
| CN106556870B (en) | A kind of well logging modeling method and system based on holding edge |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |