CN102540255B - Method for predicting extended elastic impedance inversion fluidity-containing probability - Google Patents
Method for predicting extended elastic impedance inversion fluidity-containing probability Download PDFInfo
- Publication number
- CN102540255B CN102540255B CN201110425742.5A CN201110425742A CN102540255B CN 102540255 B CN102540255 B CN 102540255B CN 201110425742 A CN201110425742 A CN 201110425742A CN 102540255 B CN102540255 B CN 102540255B
- Authority
- CN
- China
- Prior art keywords
- eei
- elastic impedance
- probability
- gas
- bearing formation
- 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 27
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 33
- 238000004445 quantitative analysis Methods 0.000 claims abstract description 17
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 32
- 238000004141 dimensional analysis Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 abstract description 6
- 238000005192 partition Methods 0.000 abstract 1
- 239000004744 fabric Substances 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000012113 quantitative test Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000000554 iris Anatomy 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
The invention aims at solving the problem that partition of air layer or water layer distributing areas is largely interfered by human factors and providing a method for predicting the extended elastic impedance inversion fluidity-containing probability. The method comprises the steps: calculating extended elastic impedance energy efficiency index (EEI) values of different incident angles of each sample according to logging information to obtain EEI (theta) j; calculating EEI differential values delta EEI j of different angles of each sample; applying a column diagram to analyze a distribution range of the delta EEI j to obtain distribution ranges of the probabilities of air layers or water layers so as to establish quantitative analysis templates of probability distribution of the air layers or the air layers and the water layers; inverting seismic data to obtain extended elastic impedance inversion cross sections of the different angles, removing non-reservoir sections and only reserving inversion data of reservoir sections; utilizing the inversion data of the reservoir sections to calculate extended elastic impedance differential value data; and applying established quantitative analysis templates to the extended elastic impedance differential value data to obtain air-containing probability data, and performing time window analysis to object layers to obtain air-containing probability plane analysis results.
Description
Technical field
The present invention relates to the reservoir prediction field of geophysical survey, more particularly, relate to a kind of Forecasting Methodology of extended elastic impedance inversion fluidity-containing probability.
Background technology
Prestack elastic impedance (EI) inverting is to utilize the well-log informations such as different geophone offsets (or incident angle) road collection data and shear wave, compressional wave, density, joint inversion goes out and lithology, multiple elastic parameter that oil-gas possibility is relevant, comprehensively judges reservoir properties and oil-gas possibility.There is following problem in conventional elastic impedance method: on yardstick, great changes will take place with angle variation for the value of asking for, also less stable of the reflection coefficient of asking for.Because seismic amplitude is with the impact of variation (AVO, the Amplitude versus Offset) effect of offset distance, the EI value of inverting changes acute variation occurs with incident angle.Closely, in, the magnitude of the elastic wave impedance value of different angles far away not in a rank, although elastic impedance value changes, be extremely small with respect to the variation of magnitude, this contrast of intersection to the elastic impedance under different incidence angles brings difficulty.
The proposition of expansion elastic impedance (EEI), has eliminated the impact of incident angle on EI, has improved the scope of application of incident angle, has improved the sensitivity to the reflection of air water difference.But at present to the analytical approach of inversion result mainly by the mode of X plot, there is the problems such as artificial property is large, analysis result poor stability.
Fig. 1 is the process flow diagram that the Forecasting Methodology of the extended elastic impedance inversion fluidity-containing probability of prior art is shown.Method shown in Fig. 1 is the analytical approach to inversion result of the prior art, main by the mode of X plot, describes the Forecasting Methodology of the extended elastic impedance inversion fluidity-containing probability of prior art below with reference to Fig. 1.
As shown in Figure 1, at step S101, calculate gas-bearing formation, water layer different angles EEI value by well-log information.
At step S102, the EEI value of the different angles to the gas-bearing formation calculating, water layer is made X plot.Because gas-bearing formation, water layer are with the variation of incident angle, the Changing Pattern of its EEI value is variant, is reflected in the distributing position difference on X plot, analyzes the main distributed areas of gas-bearing formation (or water layer), irises out Favorable Areas on X plot, sets up and explains template.
At step S103, seismic data is carried out to inverting.Obtain the expansion elastic impedance inverting data volume of different angles by inverting, and remove non-Reservoir Section, only retain the inverting value of Reservoir Section.
At step S104, utilize well-log information to make X plot and explain template, the inversion result of seismic data is made an explanation, obtain gas-bearing property and predict the outcome.
But, set up the mode of explaining template by X plot, the division of gas-bearing formation (or water layer) distributed areas is subject to interference from human factor larger, is difficult to set up unified standard, and the poor stability of explanation results is larger on the accuracy impact of explanation results.
Therefore, need a kind of Forecasting Methodology of the extended elastic impedance inversion fluidity-containing probability that is not subject to interference from human factor.
Summary of the invention
Division for the gas-bearing formation existing in prior art or water layer distributed areas is subject to the problem that interference from human factor is larger, the present invention proposes a kind of Forecasting Methodology of extended elastic impedance inversion fluidity-containing probability.The Forecasting Methodology of described extended elastic impedance inversion fluidity-containing probability can comprise the following steps: calculate the expansion elastic impedance EEI value of the different angles of each sample by well-log information, obtain EEI (θ)
j, wherein, θ represents the angle of sample, j represents catalogue number(Cat.No.); Calculate the EEI difference DELTA EEI of the different angles of each sample
j, wherein, j represents catalogue number(Cat.No.); Application histogram analysis Δ EEI
jdistribution range, obtain the probability distribution scope of gas-bearing formation and water layer, to set up the quantitative analysis template of gas-bearing formation probability distribution, or set up the quantitative analysis template of gas-bearing formation and water layer probability distribution; Seismic data is carried out to inverting, obtain the expansion elastic impedance inverting section of different angles, and remove non-Reservoir Section, only retain the inverting data of Reservoir Section; Utilize the inverting data of Reservoir Section to calculate expansion elastic impedance difference data body; The quantitative analysis template of foundation is applied to expansion elastic impedance difference data body, obtains gassiness probability data body, window analyzes to obtain gassiness probability two dimensional analysis result when zone of interest is carried out.
Obtain EEI (θ)
jstep can comprise: in target area, choose typical well, gas-bearing formation, water layer for objective interval in the typical well of choosing are determined analyzing samples, the expansion elastic impedance EEI value of calculating different angles for definite analyzing samples, obtains EEI (θ)
j.
The Forecasting Methodology of extended elastic impedance inversion fluidity-containing probability is mainly set up quantitative analysis template by the calculating of EEI difference and histogram according to an exemplary embodiment of the present invention, quantitative analysis template is applied to real data body (expansion elastic impedance difference data body), finally obtains the probability distribution graph of fluid (gassiness or moisture).Apply after the method for this quantitative test expansion elastic impedance inversion result, can reduce the artificial property that adopts traditional X plot delineation range of profitability analytical approach, improve precision and the efficiency of analysis result.
Accompanying drawing explanation
By the detailed description of carrying out below in conjunction with accompanying drawing, above and other objects of the present invention, feature and advantage will become apparent, wherein:
Fig. 1 is the process flow diagram that the Forecasting Methodology of the extended elastic impedance inversion fluidity-containing probability of prior art is shown;
Fig. 2 illustrates the process flow diagram of the Forecasting Methodology of extended elastic impedance inversion fluidity-containing probability according to an exemplary embodiment of the present invention;
Fig. 3 is that the incident angle of the gas-bearing formation (representing with square) that calculated by well-log information, water layer (representing with pentagram) is the EEI value X plot of 5 ° and 25 °;
Fig. 4 is the histogram of far away, the near incident angle expansion elastic impedance difference cloth of statistics gas-bearing formation and water layer;
Fig. 5 is the sectional view that certain well has been removed two sections of expansion elastic impedance inverting differences of palpus of non-reservoir;
Fig. 6 is the sectional view containing fluidity probability that utilization fractional analysis template obtains;
Fig. 7 is the planimetric map of the gassiness probability of the objective interval that obtains by probability statistics.
Embodiment
Now, describe more fully with reference to the accompanying drawings according to example embodiment of the present invention.
Fig. 2 illustrates the process flow diagram of the Forecasting Methodology of extended elastic impedance inversion fluidity-containing probability according to an exemplary embodiment of the present invention.
With reference to Fig. 2, at step S201, calculate expansion elastic impedance (EEI) value of the different angles of each sample by well-log information, obtain EEI (θ)
j, wherein, θ represents the angle of sample, j represents catalogue number(Cat.No.).Specifically, can in target area, choose typical well, gas-bearing formation, water layer for objective interval in the typical well of choosing are determined analyzing samples, calculate the expansion elastic impedance EEI value of different angles for definite analyzing samples, obtain EEI (θ)
j.
Because gas-bearing formation, water layer are with the variation of incident angle, the Changing Pattern of its EEI value is variant.Therefore,, at step S202, can calculate the EEI difference DELTA EEI of the different angles of each sample
j, j represents catalogue number(Cat.No.);
At step S203, set up the quantitative analysis template of gas-bearing formation (or water layer) probability distribution.Specifically, application histogram analysis Δ EEI
jdistribution range, obtain the probability distribution scope of gas-bearing formation and water layer, to set up the quantitative analysis template of gas-bearing formation probability distribution, or set up the quantitative analysis template of gas-bearing formation and water layer probability distribution.
At step S204, seismic data is carried out to inverting, obtain the expansion elastic impedance inverting section of different angles, and remove non-Reservoir Section, only retain the inverting data of Reservoir Section.
At step S205, utilize the inverting data of Reservoir Section to calculate expansion elastic impedance difference data body.
At step S206, the quantitative analysis template that will set up at step S203 be applied to expansion elastic impedance difference data body, obtain gassiness probability data body, and window analyzes to obtain gassiness probability two dimensional analysis result when zone of interest is carried out.
To describe by Fig. 3 to Fig. 7 by the various diagrams of the Forecasting Methodology of the extended elastic impedance inversion fluidity-containing probability according to an exemplary embodiment of the present invention shown in Fig. 2 below.
Fig. 3 is that the incident angle of the gas-bearing formation (representing with square) that calculated by well-log information, water layer (representing with pentagram) is the EEI value X plot of 5 ° and 25 °.With reference to Fig. 3, in the situation that little incident angle elastic impedance is identical, water layer is obviously greater than gas-bearing formation the EEI value of 25 °, this shows, incident angle expansion elastic impedance difference far away by analyzing, near, can reach the object of distinguishing water layer and gas-bearing formation.
Fig. 4 is the histogram of far away, the near incident angle expansion elastic impedance difference cloth of statistics gas-bearing formation and water layer, wherein, and Δ EEI
j=EEI (25 °)
j-EEI (5 °)
j, transverse axis is Δ EEI value, the longitudinal axis on the top of Fig. 4 is that the accumulation frequency appears in sampling point, i.e. and distribution probability, the longitudinal axis of the bottom of Fig. 4 is the frequency that the sampling point of statistics occurs.The Δ EEI of the gas-bearing formation in Fig. 4
jmainly be distributed in below 650, more than 80% be distributed in below 400, and more than 80% water layer is distributed in the 450～1100,450～650th, the interval that gas-bearing formation and water layer all occur, but the probability occurring is different.Can obtain the distribution probability relation of gas-bearing formation, water layer according to the statistics shown in Fig. 4.
Fig. 5 be certain well (for example, Peng Lai 2 wells) remove the sectional view of two sections of the palpuses expansion elastic impedance inverting differences of non-reservoir, near well in Fig. 5, palpus two epimere elastic impedance differences are relatively little, mostly below 700, and palpus two hypomere elastic impedance differences are large, take more than 800 high values as main, this well practical logging explanation reservoir is mainly grown at palpus two epimeres and hypomere, wherein epimere is mainly gas-bearing formation, and hypomere is water layer.
Fig. 6 is the sectional view containing fluidity probability that utilization fractional analysis template obtains, and the upper graph of Fig. 6 represents the sectional view of gassiness probability, and the lower graph of Fig. 6 represents the sectional view of moisture probability.As can be seen from Figure 6, must two epimere gassiness probability higher, most of more than 60%, must two hypomere gassiness probability lower, major part be less than the probability shown in the upper graph of 20%, Fig. 6 and the probability shown in lower graph and be complementary trend.
Fig. 7 is the planimetric map of the gassiness probability of the objective interval that obtains by probability statistics, the most shallow region representation gassiness probability of color is high, the darkest region representation gassiness probability of color is lower, color region representation gassiness probability placed in the middle is low, from color region placed in the middle to the darkest region of color again to the most shallow region of color, gassiness probability is increased to 1 gradually, can identify quantitatively thus gassiness Favorable Areas.
As can be seen from the above description, the Forecasting Methodology of extended elastic impedance inversion fluidity-containing probability is mainly set up quantitative analysis template by the calculating of EEI difference and histogram according to an exemplary embodiment of the present invention, quantitative analysis template is applied to real data body (expansion elastic impedance difference data body), the final probability distribution graph obtaining containing fluid (gassiness or moisture).Apply after the method for this quantitative test expansion elastic impedance inversion result, can reduce the artificial property that adopts traditional X plot delineation range of profitability analytical approach, improve precision and the efficiency of analysis result.In addition, the Forecasting Methodology of extended elastic impedance inversion fluidity-containing probability according to an exemplary embodiment of the present invention, predict the gassiness probability distribution of continental clastic gas reservoir, predict the outcome and show that high gas rate well and air water are arranged in height, gassiness Probability Region mostly with the industrial gas well producing, its gas-bearing property prediction coincidence rate reaches 81%, there is higher precision of prediction, continental clastic gas reservoir is had a good application prospect.
Although shown with reference to certain exemplary embodiments of the present invention and described the present invention, but it should be appreciated by those skilled in the art, in the case of not departing from the spirit and scope of the present invention of claim and equivalent restriction thereof, can carry out in form and details various changes.
Claims (1)
1. a Forecasting Methodology for extended elastic impedance inversion fluidity-containing probability, is characterized in that, described Forecasting Methodology comprises the following steps:
The expansion elastic impedance EEI value of calculating the different incidence angles degree of each sample by well-log information, obtains EEI (θ)
j, wherein, θ represents the angle of sample, j represents catalogue number(Cat.No.);
Calculate the EEI difference DELTA EEI of the different angles of each sample
j, wherein, j represents catalogue number(Cat.No.);
Application histogram analysis Δ EEI
jdistribution range, obtain the probability distribution scope of gas-bearing formation and water layer, to set up the quantitative analysis template of gas-bearing formation probability distribution, or set up the quantitative analysis template of gas-bearing formation and water layer probability distribution;
Seismic data is carried out to inverting, obtain the expansion elastic impedance inverting section of different angles, and remove non-Reservoir Section, only retain the inverting data of Reservoir Section;
Utilize the inverting data of Reservoir Section to calculate expansion elastic impedance difference data body;
The quantitative analysis template of the gas-bearing formation of the quantitative analysis template of the gas-bearing formation probability distribution of foundation or foundation and water layer probability distribution is applied to expansion elastic impedance difference data body, obtain gassiness probability data body, window analyzes to obtain gassiness probability two dimensional analysis result when zone of interest is carried out
Wherein, the step that obtains EEI (θ) j comprises: in target area, choose typical well, gas-bearing formation, water layer for objective interval in the typical well of choosing are determined analyzing samples, the expansion elastic impedance EEI value of calculating different angles for definite analyzing samples, obtains EEI (θ)
j.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110425742.5A CN102540255B (en) | 2011-12-19 | 2011-12-19 | Method for predicting extended elastic impedance inversion fluidity-containing probability |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110425742.5A CN102540255B (en) | 2011-12-19 | 2011-12-19 | Method for predicting extended elastic impedance inversion fluidity-containing probability |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102540255A CN102540255A (en) | 2012-07-04 |
| CN102540255B true CN102540255B (en) | 2014-05-21 |
Family
ID=46347565
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201110425742.5A Active CN102540255B (en) | 2011-12-19 | 2011-12-19 | Method for predicting extended elastic impedance inversion fluidity-containing probability |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102540255B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104656137B (en) * | 2015-02-27 | 2017-07-28 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Elastic wave impedance linear combination gas-oil detecting method |
| CN105301643A (en) * | 2015-11-13 | 2016-02-03 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Geological reservoir detection method based on extended elastic impedance |
| CN112180442A (en) * | 2019-07-03 | 2021-01-05 | 中国石油天然气集团有限公司 | Lithology inversion method and system based on CRP (common reflection point) gather |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101907729A (en) * | 2010-06-11 | 2010-12-08 | 北京诺克斯达石油科技有限公司 | Ray elastic parameter inversion method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007001759A2 (en) * | 2005-06-24 | 2007-01-04 | Exxonmobil Upstream Research Company | Method for obtaining porosity and shale volume from seismic data |
-
2011
- 2011-12-19 CN CN201110425742.5A patent/CN102540255B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101907729A (en) * | 2010-06-11 | 2010-12-08 | 北京诺克斯达石油科技有限公司 | Ray elastic parameter inversion method |
Non-Patent Citations (2)
| Title |
|---|
| 扩展弹性阻抗反演在储层预测中的应用;牛聪等;《石油地球物理勘探》;20111130;第46卷(第增刊1期);67-71 * |
| 牛聪等.扩展弹性阻抗反演在储层预测中的应用.《石油地球物理勘探》.2011,第46卷(第增刊1期),67-71. |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102540255A (en) | 2012-07-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102033242B (en) | Deep inclined fractured reservoir earthquake amplitude prediction method | |
| CN108931814B (en) | Multi-attribute fusion based bedrock fracture prediction method | |
| CN102466815B (en) | Triassic clastic rock petroleum reservoir identification method | |
| CN104155701B (en) | A kind of multi-scale facture Forecasting Methodology utilizing Prestack seismic data and well information | |
| CN111506861B (en) | Method for calculating crack strength of favorable region of target layer | |
| CN104678434A (en) | Method for predicting storage layer crack development parameters | |
| CN102109611A (en) | Fast and convenient method for predicting high-quality petroleum reservoir in virtue of seism attributes | |
| CN112505754B (en) | Method for collaborative partitioning sedimentary microfacies by well-seismic based on high-precision sequence grid model | |
| CN103777245A (en) | Method of quantitative evaluation on oil and gas accumulation conditions based on seismic data | |
| CN103869362A (en) | Method and equipment for obtaining body curvature | |
| CN104570086B (en) | Method for predicting pre-stack cracks in common offset and common azimuth angle domain | |
| CN102540255B (en) | Method for predicting extended elastic impedance inversion fluidity-containing probability | |
| CN114994758A (en) | Wave impedance extraction and structure characterization method and system for carbonate fracture control reservoir | |
| CN109283577B (en) | Seismic horizon calibration method | |
| CN113109875B (en) | Inversion method of carbonate rock reservoir under full waveform velocity field constraint | |
| CN117930380A (en) | Carbonate reservoir prediction method and system based on seismology | |
| CN111352154B (en) | Reservoir prediction method based on wide-azimuth earthquake | |
| CN108931813B (en) | Method for identifying reef cover, reef core and reef base | |
| CN114488305B (en) | A fine calibration method for geological horizons of seismic data in new exploration areas without wells | |
| CN115267937A (en) | Braided river delta sand body configuration prediction method and system under offshore thin well pattern condition | |
| CN112198554B (en) | Method for establishing high-precision inversion initial model driven by seismic waveform | |
| CN113589365B (en) | Reservoir pinch-out line description method based on time-frequency domain information | |
| CN112444857A (en) | Paste-salt earth physical prediction method based on post-stack seismic data | |
| CN115113266B (en) | Volcanic depth domain phased velocity modeling method, storage medium and computer equipment | |
| CN114861515B (en) | Method, device, equipment and medium for calculating layer speed data body |
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 | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20120704 Assignee: SICHUAN JISAITE TECHNOLOGY CO.,LTD. Assignor: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. Contract record no.: 2016510000019 Denomination of invention: Method for predicting extended elastic impedance inversion fluidity-containing probability Granted publication date: 20140521 License type: Exclusive License Record date: 20160920 |
|
| LICC | Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180131 Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189 Patentee after: BGP INC., CHINA NATIONAL PETROLEUM Corp. Address before: 610213 No. 1, No. 1, No. 1, Huayang Avenue, Huayang Town, Shuangliu County, Chengdu, Sichuan Patentee before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200921 Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No. Patentee after: CHINA NATIONAL PETROLEUM Corp. Patentee after: BGP Inc., China National Petroleum Corp. Address before: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189 Patentee before: BGP Inc., China National Petroleum Corp. |
|
| TR01 | Transfer of patent right |