CN111158048B - Analysis method for improving reservoir prediction precision through seismic waveform envelope interpretation - Google Patents
Analysis method for improving reservoir prediction precision through seismic waveform envelope interpretation Download PDFInfo
- Publication number
- CN111158048B CN111158048B CN202010007364.8A CN202010007364A CN111158048B CN 111158048 B CN111158048 B CN 111158048B CN 202010007364 A CN202010007364 A CN 202010007364A CN 111158048 B CN111158048 B CN 111158048B
- Authority
- CN
- China
- Prior art keywords
- seismic
- interpretation
- waveform
- time window
- layer
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 24
- 238000006243 chemical reaction Methods 0.000 claims abstract description 7
- 238000000605 extraction Methods 0.000 claims abstract description 7
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 6
- 238000010586 diagram Methods 0.000 claims abstract description 5
- 238000004519 manufacturing process Methods 0.000 claims description 6
- 238000005553 drilling Methods 0.000 claims description 4
- 239000004576 sand Substances 0.000 claims description 4
- 239000003208 petroleum Substances 0.000 abstract description 4
- 238000011160 research Methods 0.000 abstract description 4
- 238000000034 method Methods 0.000 abstract description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- 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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides an analysis method for improving reservoir prediction precision by seismic waveform envelope interpretation, which comprises the following steps: 1) Calibrating a seismic waveform main lobe corresponding to a research target interval on a seismic reflection waveform section by finely synthesizing seismic record calibration and seismic phase 90-degree conversion by using seismic data interpretation software; 2) Defining an interpretation horizon top layer and an interpretation bottom layer, respectively carrying out horizon tracking interpretation along the top and bottom of a seismic waveform main lobe corresponding to a target layer, and interpreting density 2x2; 3) Explaining a time window between the top layer and the bottom layer of the horizon to extract a time window for seismic waveform attributes and determining the time window as a non-isochronous optimal window for waveform seismic attribute extraction analysis; 4) Extracting and analyzing the arc length, the deflection and the torsion of the earthquake waveform in the time window in the non-isochronous optimal window, drawing an attribute plane diagram and carrying out reservoir prediction; the prediction precision is higher, the method is closer to the actual geological reality reflected by the actual target interval, and better foundation geological basis can be provided for petroleum exploration decision-making.
Description
Technical Field
The invention relates to the technical field of oil exploration, in particular to an analysis method for improving reservoir prediction precision through seismic waveform envelope interpretation.
Background
The reservoir prediction is carried out by utilizing the seismic waveform attribute extraction, and the method is widely applied to oil exploration and development. The basic principles of seismic stratigraphy assume that the event of a reflected wave on a seismic section has the meaning of a geologic time interface. The reservoir prediction by using the seismic waveform attribute must give an attribute analysis time window, and the conventional attribute analysis time window is determined by firstly tracing and explaining a seismic horizon along the seismic waveform peak homophase axis, then respectively drifting upwards and downwards for a certain time along the horizon, and using the time window to perform the seismic waveform attribute analysis. However, the seismic waveform information in the time window cannot truly reflect the real information of the target layer of the research, and often contains seismic waveform information of a non-target layer, or the seismic waveform information of the target layer is not contained in the analysis time window because the definition of the time window is not proper, so that the reservoir prediction precision by utilizing the seismic waveform attribute is influenced; and an 'isochronous window' generated by the isochronous drift along the layer is used for seismic waveform attribute prediction, the analysis time window is not set accurately, and non-target layer seismic waveform attribute information is contained in the analysis window.
Disclosure of Invention
In view of the technical deficiencies, the invention aims to provide an analysis method for improving the reservoir prediction precision by seismic waveform envelope interpretation, which carries out 'seismic waveform envelope interpretation' on seismic waveform information reflecting a target interval, determines an optimal time window reflecting the seismic waveform information of the target interval to the maximum extent, is more accurate than a time window which is drifted up and down along the interval by conventional interpretation, has higher reservoir prediction precision by using the method, is closer to the geological reality reflected by the actual target interval, and can provide better foundation geological basis for petroleum exploration decision.
In order to solve the technical problems, the invention adopts the following technical scheme:
an analysis method for improving reservoir prediction accuracy by seismic waveform envelope interpretation, which utilizes a computer provided with seismic data interpretation software to perform seismic waveform analysis, and is characterized by comprising the following steps:
1) Calibrating a seismic waveform main lobe corresponding to a research target interval on a seismic reflection waveform section by finely synthesizing seismic record calibration and seismic phase 90-degree conversion by using seismic data interpretation software;
2) Defining an interpretation horizon top layer and an interpretation bottom layer, respectively performing horizon tracking interpretation along the top and bottom of a main lobe of a seismic waveform corresponding to a target layer, interpreting density 2x2, and performing 1x1 interpolation after the interpretation is completed;
3) Explaining a time window between the top layer and the bottom layer of the horizon, extracting the time window for seismic waveform attributes, and determining the time window as a non-isochronous optimal window for waveform seismic attribute extraction analysis;
4) And extracting and analyzing the arc length, the deflection and the torsion of the earthquake waveform in the time window in the non-isochronous optimal window, drawing an attribute plane diagram and carrying out reservoir prediction.
The invention has the beneficial effects that: the invention carries out 'seismic waveform enveloping interpretation' on the seismic waveform information reflecting the target interval, determines the optimal time window reflecting the seismic waveform information of the target interval to the maximum extent, is more accurate than the time window which is opened by the conventional interpretation along the up-and-down drift of the layer, has higher reservoir prediction precision by utilizing the invention, is more close to the geological reality reflected by the actual target interval, and can provide better foundation geological basis for petroleum exploration decision.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a flow chart of an analysis method for improving reservoir prediction accuracy through seismic waveform envelope interpretation according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 1, an analysis method for improving reservoir prediction accuracy by seismic waveform envelope interpretation utilizes a computer installed with seismic data interpretation software to perform seismic waveform analysis, and comprises the following steps:
1) Calibrating a seismic waveform main lobe corresponding to a research target interval on a seismic reflection waveform section by finely synthesizing seismic record calibration and seismic phase 90-degree conversion by using seismic data interpretation software;
2) Defining an interpretation horizon top layer and an interpretation bottom layer, respectively performing horizon tracking interpretation along the top and bottom of a main lobe of a seismic waveform corresponding to a target layer, interpreting density 2x2, and performing 1x1 interpolation after the interpretation is completed;
3) Explaining a time window between the top layer and the bottom layer of the horizon, extracting the time window for seismic waveform attributes, and determining the time window as a non-isochronous optimal window for waveform seismic attribute extraction analysis;
4) And extracting and analyzing the arc length, the deflection and the torsion of the earthquake waveform in the time window in the non-isochronous optimal window, drawing an attribute plane diagram and carrying out reservoir prediction.
When the seismic interpretation data processing method is used, firstly, the obtained original seismic data (sgy standard data), well drilling and logging curve data, well drilling and well deviation data and well position coordinate data are sorted and standardized, and a seismic interpretation work area is established; determining a target layer through stratum division and comparison, synthesizing a seismic record manufacturing module by using software, manufacturing a fine synthetic seismic record to carry out the corresponding relation of the target layer section and a seismic reflection waveform in-phase axis, corresponding the target layer (a single sand body or a sand layer group combination) to a seismic reflection main lobe through 90-degree phase conversion, and synthesizing seismic record manufacturing and 90-degree phase conversion.
Secondly, defining an interpretation horizon top layer (top) and an interpretation bottom layer (bot), and respectively carrying out horizon tracking interpretation along the top and bottom of the seismic waveform main lobe corresponding to the target layer, wherein the interpretation density is 2x2, and 1x1 is preferred. If the interpretation is 2x2, 1x1 interpolation is performed after the interpretation is finished.
And after the top layer (top) and the bottom layer (bot) are interpreted, the time window between the two layers is the 'non-isochronous optimal window' of the waveform seismic attribute extraction analysis through the horizon interpolation calculation. Compared with an 'isochronous window' obtained by shifting up and down for a certain millisecond number after the homophase axis tracking of a conventional peak, the 'non-isochronous optimal window' can completely contain the seismic information corresponding to the target interval in an analysis window.
And finally, extracting and analyzing seismic waveform arc length, deflection and torsion seismic attributes in the time window by utilizing a software corresponding attribute extraction module in the non-isochronous optimal window, drawing a corresponding attribute plane diagram, and developing reservoir prediction.
The invention carries out 'seismic waveform enveloping interpretation' on the seismic waveform information reflecting the target interval, determines the optimal time window reflecting the seismic waveform information of the target interval to the maximum extent, is more accurate than the time window which is opened by the conventional interpretation along the up-and-down drift of the layer, has higher reservoir prediction precision by utilizing the invention, is more close to the geological reality reflected by the actual target interval, and can provide better foundation geological basis for petroleum exploration decision.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include such modifications and variations.
Claims (1)
1. An analysis method for improving reservoir prediction accuracy by seismic waveform envelope interpretation, which utilizes a computer provided with seismic data interpretation software to perform seismic waveform analysis, and is characterized by comprising the following steps:
1) Arranging and standardizing the obtained original seismic data volume, well drilling and logging curve data, well drilling well deviation data and well position coordinate data, establishing an earthquake interpretation work area, determining a target layer through stratum division and comparison, synthesizing a seismic record manufacturing module by using software, manufacturing a fine synthetic seismic record, carrying out the same-phase axis corresponding relation between a target layer and a seismic reflection waveform, enabling the target layer, namely a single sand body or sand layer group combination to correspond to a seismic reflection main lobe through 90-degree phase conversion, manufacturing the synthetic seismic record and carrying out 90-degree phase conversion, and calibrating the seismic waveform main lobe corresponding to the target layer to be researched on a seismic reflection waveform section;
2) Defining an interpretation horizon top layer and an interpretation bottom layer, respectively carrying out horizon tracking interpretation along the top and bottom of a seismic waveform main lobe corresponding to a target layer, interpreting density 2x2, and carrying out 1x1 interpolation after the interpretation is finished;
3) Explaining a time window between the top layer and the bottom layer of the horizon to extract a time window for seismic waveform attributes and determining the time window as a non-isochronous optimal window for waveform seismic attribute extraction analysis;
4) And extracting and analyzing the arc length, the deflection and the torsion of the earthquake waveform in the time window in the non-isochronous optimal window, drawing an attribute plane diagram and carrying out reservoir prediction.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010007364.8A CN111158048B (en) | 2020-01-04 | 2020-01-04 | Analysis method for improving reservoir prediction precision through seismic waveform envelope interpretation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010007364.8A CN111158048B (en) | 2020-01-04 | 2020-01-04 | Analysis method for improving reservoir prediction precision through seismic waveform envelope interpretation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111158048A CN111158048A (en) | 2020-05-15 |
CN111158048B true CN111158048B (en) | 2022-10-04 |
Family
ID=70561257
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010007364.8A Active CN111158048B (en) | 2020-01-04 | 2020-01-04 | Analysis method for improving reservoir prediction precision through seismic waveform envelope interpretation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111158048B (en) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6463387B1 (en) * | 2001-01-31 | 2002-10-08 | Phillips Petroleum Company | 3-D seismic event tracking |
CN101408624A (en) * | 2007-10-08 | 2009-04-15 | 陶庆学 | Forecasting and evaluating technologies of three-dimensional earthquake optimum time window river course sand body storage layer |
CN102650701A (en) * | 2012-04-27 | 2012-08-29 | 中国石油天然气股份有限公司 | Reservoir wave impedance prediction method based on dipole wave |
CN104898164A (en) * | 2015-03-23 | 2015-09-09 | 中国石油天然气股份有限公司 | thin tight reservoir earthquake prediction method based on earthquake phase tiny change analysis |
CN105093306A (en) * | 2014-05-15 | 2015-11-25 | 中国石油化工股份有限公司 | Method for automatic interpretation and thickness calculation of reservoir in geophysical exploration |
CN107132573A (en) * | 2017-02-28 | 2017-09-05 | 中国石油化工股份有限公司 | A kind of method that application wavelet decomposition reconfiguration technique recognizes the lower lithological pool of strong impedance shielding |
CN109188520A (en) * | 2018-09-17 | 2019-01-11 | 中国石油天然气股份有限公司 | Thin reservoir thickness prediction method and device |
CN109725348A (en) * | 2017-10-30 | 2019-05-07 | 中国石油化工股份有限公司 | A method of sedimentary facies is identified based on seismic data |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7800977B2 (en) * | 2004-06-01 | 2010-09-21 | Westerngeco L.L.C. | Pre-stack combining of over/under seismic data |
FR2909185B1 (en) * | 2006-11-27 | 2009-01-09 | Inst Francais Du Petrole | METHOD OF STRATIGRAPHIC INTERPRETATION OF SEISMIC IMAGES |
US7577529B2 (en) * | 2007-05-22 | 2009-08-18 | Chevron U.S.A. Inc. | Method for determining attributes associated with net-sand thickness |
CN103412332B (en) * | 2013-01-22 | 2016-05-25 | 中国地质大学(北京) | A kind of method of determining thin reservoir thickness |
CN110068859A (en) * | 2019-04-12 | 2019-07-30 | 中国石油大学(华东) | A kind of horizontal well region braid deltas seismic reservoir sedimentology quantitatively characterizing method |
-
2020
- 2020-01-04 CN CN202010007364.8A patent/CN111158048B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6463387B1 (en) * | 2001-01-31 | 2002-10-08 | Phillips Petroleum Company | 3-D seismic event tracking |
CN101408624A (en) * | 2007-10-08 | 2009-04-15 | 陶庆学 | Forecasting and evaluating technologies of three-dimensional earthquake optimum time window river course sand body storage layer |
CN102650701A (en) * | 2012-04-27 | 2012-08-29 | 中国石油天然气股份有限公司 | Reservoir wave impedance prediction method based on dipole wave |
CN105093306A (en) * | 2014-05-15 | 2015-11-25 | 中国石油化工股份有限公司 | Method for automatic interpretation and thickness calculation of reservoir in geophysical exploration |
CN104898164A (en) * | 2015-03-23 | 2015-09-09 | 中国石油天然气股份有限公司 | thin tight reservoir earthquake prediction method based on earthquake phase tiny change analysis |
CN107132573A (en) * | 2017-02-28 | 2017-09-05 | 中国石油化工股份有限公司 | A kind of method that application wavelet decomposition reconfiguration technique recognizes the lower lithological pool of strong impedance shielding |
CN109725348A (en) * | 2017-10-30 | 2019-05-07 | 中国石油化工股份有限公司 | A method of sedimentary facies is identified based on seismic data |
CN109188520A (en) * | 2018-09-17 | 2019-01-11 | 中国石油天然气股份有限公司 | Thin reservoir thickness prediction method and device |
Also Published As
Publication number | Publication date |
---|---|
CN111158048A (en) | 2020-05-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2018340369B2 (en) | Method and device for determining thin interlayer | |
CN109188520A (en) | Thin reservoir thickness prediction method and device | |
CN102866426B (en) | A kind of method utilizing AVO wide-angle road set analysis rock mass hydrocarbon information | |
CN108121009B (en) | Development late stage complex fault block velocity-variant depth mapping method | |
CN105549084A (en) | Method and system for building three-dimensional high-precision velocity model | |
CN107132574B (en) | Prediction method for lithologic combination of marine stratum clastic rocks | |
CN111158048B (en) | Analysis method for improving reservoir prediction precision through seismic waveform envelope interpretation | |
CN105093266A (en) | Seismogeological comprehensive description method for stratigraphic overlap line | |
CN110727030B (en) | Reservoir prediction method based on wavelet multi-extreme characteristic joint analysis | |
Liu et al. | Wave-equation diffraction imaging using pseudo dip-angle gather | |
CN113031070B (en) | Method for making depth domain synthetic seismic record | |
CN110632660B (en) | Thin sand body characterization method and device based on seismic data body | |
Galloway et al. | Seismic stratigraphic model of depositional platform margin, eastern Anadarko basin, Oklahoma | |
Zhang et al. | Identifying minor faults on top of coalfield Ordovician limestone stratum using seismic attributes derived from azimuthally stacked data | |
CN112147676A (en) | Method for predicting thickness of coal bed and gangue | |
CN112394398A (en) | Stratum pinch-out line identification method based on lithofacies and seismic facies | |
CN111665558A (en) | Oil gas distribution detection method and system | |
CN109143354B (en) | Method and device for decomposing seismic waveform characteristics | |
Wang et al. | Lithological Reservoir Exploration Technology | |
CN114722852A (en) | CNN-based ablation point identification method | |
Rigatti et al. | 3-D PSDM case history in a thrust belt: Quiriquire Block, Eastern Venezuela Basin | |
CN116338786A (en) | Earthquake inversion method and system for eliminating strong reflection shielding effect | |
Fei et al. | Seismic characterization of individual geologic factors with disentangled features | |
CN118011487A (en) | Beach dolomite reservoir earthquake prediction method and earthquake prediction equipment | |
CN117950017A (en) | Reservoir prediction parameter fusion method and device |
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 |