CN106199707B - A kind of method and device of prediction sandbody distribution - Google Patents
A kind of method and device of prediction sandbody distribution Download PDFInfo
- Publication number
- CN106199707B CN106199707B CN201610431229.XA CN201610431229A CN106199707B CN 106199707 B CN106199707 B CN 106199707B CN 201610431229 A CN201610431229 A CN 201610431229A CN 106199707 B CN106199707 B CN 106199707B
- Authority
- CN
- China
- Prior art keywords
- seismic data
- data cube
- survey area
- plane layer
- phase
- 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
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. analysis, for interpretation, for correction
- 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. analysis, for interpretation, for correction
- 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6161—Seismic or acoustic, e.g. land or sea measurements
-
- 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
-
- 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase
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 present invention relates to technical field of geological exploration more particularly to it is a kind of prediction sandbody distribution method, the method includes:After exciting artificial earthquake, the first seismic data cube in survey area is obtained;Phase of first seismic data cube at 1/4 wavelength is inverted, the second seismic data cube is obtained;Extract the plane layer attribute of second seismic data cube;According to the plane layer attribute, the sandbody distribution of the survey area is predicted.The present invention is by inverting phase of first seismic data cube at 1/4 wavelength, obtain the second seismic data cube, to which the wave crest of the second seismic data cube no longer corresponds to the interface of mud stone and sandstone, the upper zero phase of second seismic data cube just corresponds to the reflecting interface of mud stone and sandstone, seismic reflection main lobe is transformed into thin layer center by the present invention by 90 ° of phase transitions, and the earthquake axis to be obtained according to the second seismic data cube is equal to actual geology axis.
Description
Technical field
The present invention relates to technical field of geological exploration more particularly to a kind of method and devices of prediction sandbody distribution.
Background technology
Seismic prospecting refer to earth's surface excite man-made explosion, from the vibrations caused by man-made explosion in the form of seismic wave to
Underground propagation, and earth's surface is passed in reflection back upwards under certain conditions, by the reflected seismic wave of the instrument record of earth's surface, from
And obtain seismic data.Due to seismic data be via obtained from underground medium, will necessarily be by the object of underground medium
The influence of rationality matter, this influence can be also saved in by the integrated information as underground medium in seismic data.Pass through specific aim
Seismic data is handled and is parsed, useful subsurface information can be obtained, and then geologic objective is instructed to explore.
With the progress of 3-d seismic exploration technology, and the deep application in oil field prospecting field, utilize seismic data
Seismic properties carry out predicting reservoir and become a kind of common application.The prior art, the ground that will usually be obtained according to seismic data
Shake axis is equal to geology axis, to predict the sandbody distribution in geology, still, and due to when exciting artificial earthquake, low speed mud stone
Stratum encounters high speed sandstone formation, and sound wave curve will will produce take-off, using standard zero phase seismic data be equal geology axis come
The result that prediction geology sandbody distribution obtains is different from practical geological structure, can not accurately predict geology sandbody distribution, no
Geological prospecting needs can be met.
Invention content
The present invention solves by providing a kind of method and device of prediction sandbody distribution and utilizes earthquake axis in the prior art
It is equal to geology axis to predict that the result that geology sandbody distribution obtains is different from practical geological structure, can not accurately predicts geology
The technical issues of sandbody distribution.
An embodiment of the present invention provides it is a kind of prediction sandbody distribution method, the method includes:
After exciting artificial earthquake, the first seismic data cube in survey area is obtained;
Phase of first seismic data cube at 1/4 wavelength is inverted, the second seismic data cube is obtained;
Extract the plane layer attribute of second seismic data cube;
According to the plane layer attribute, the sandbody distribution of the survey area is predicted.
Preferably, it before the sandbody distribution for predicting the survey area, is also wrapped according to the plane layer attribute described
It includes:
It is layered according to the drilling geology of the survey area, at least one stratum for including to second seismic data cube
The depth of each formation information is marked in information.
Preferably, described that the sandbody distribution of the survey area is predicted according to the plane layer attribute, including:
According to the depth and plane layer attribute of each formation information in second seismic data cube, the exploration is predicted
Relative position relation in region between the type and form and each geological stratification of geological stratification corresponding with each formation information.
Preferably, plane stratum attribute information includes amplitude and/or phase and/or frequency.
Based on same inventive concept, the embodiment of the present invention also provides a kind of device of prediction sandbody distribution, described device packet
It includes:
Acquisition module, for after exciting artificial earthquake, obtaining the first seismic data cube in survey area;
Reversal block obtains the second seismic data for inverting phase of first seismic data cube at 1/4 wavelength
Body;
Extraction module, the plane layer attribute for extracting second seismic data cube;
Prediction module, for according to the plane layer attribute, predicting the sandbody distribution of the survey area.
Preferably, described device further includes:
Mark module, for according to the layering of the drilling geology of the survey area, including to second seismic data cube
At least one formation information in the depth of each formation information be marked.
Preferably, the prediction module is specifically used for:
According to the depth and plane layer attribute of each formation information in second seismic data cube, the exploration is predicted
Relative position relation in region between the type and form and each geological stratification of geological stratification corresponding with each formation information.
Preferably, plane stratum attribute information includes amplitude and/or phase and/or frequency.
One or more of embodiment of the present invention technical solution, has at least the following technical effects or advantages:
Phase of the present invention by the first seismic data cube of reversion at 1/4 wavelength, obtains the second seismic data cube, to
The wave crest of second seismic data cube no longer corresponds to the interface of mud stone and sandstone, and the upper zero phase of the second seismic data cube is just right
The reflecting interface of mud stone and sandstone, lower zero phase is answered just to correspond to the reflecting interface of sandstone and mud stone, the present invention passes through 90 ° of phases
Seismic reflection main lobe is transformed into thin layer center by conversion, and the earthquake axis to be obtained according to the second seismic data cube is equal to reality
The geology axis on border.
Description of the drawings
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technology description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
The embodiment of invention for those of ordinary skill in the art without creative efforts, can also basis
The attached drawing of offer obtains other attached drawings.
Fig. 1 is a kind of method of prediction sandbody distribution in the embodiment of the present invention;
Fig. 2 is a kind of device of prediction sandbody distribution in the embodiment of the present invention.
Specific implementation mode
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art
The every other embodiment obtained without creative efforts, shall fall within the protection scope of the present invention.
The embodiment of the present invention provides a kind of method of prediction sandbody distribution, as shown in Figure 1, the method includes:
Step 101:After exciting artificial earthquake, the first seismic data cube in survey area is obtained.
Step 102:Phase of first seismic data cube at 1/4 wavelength is inverted, the second seismic data cube is obtained.
Step 103:Extract the plane layer attribute of second seismic data cube.
Step 104:According to the plane layer attribute, the sandbody distribution of the survey area is predicted.
Specifically, in embodiments of the present invention, man-made explosion is set in the earth's surface of survey area first, with manual method
Earthquake-wave-exciting, seismic wave encounters medium character different rock stratum interface when to underground propagation will occur to reflect and reflect.
Then, in a step 101, after the completion of artificial earthquake excitation, the first seismic data cube is obtained in survey area, specifically, logical
It crosses in the drilling well in earth's surface or survey area and receives seismic wave with wave detector, to obtain the first seismic data cube.
Further, take-off can occur since low speed mud stone encounters high speed sandstone, leads to Thin Sandbody seismic reflection and earthquake
It is not correspondence between lineups, the wave crest of the first seismic data cube corresponds to sandstone top, the trough pair of the first seismic data cube
Sandstone bottom is answered, therefore standard zero phase seismic data corresponded manner is not suitable for the lithologic interpretation of Thin Sandbody.Of the invention real
It applies in example, phase of the step 102 by the first seismic data cube of reversion at 1/4 wavelength obtains the second seismic data cube, to
The wave crest of second seismic data cube no longer corresponds to the interface of mud stone and sandstone, and the upper zero phase of the second seismic data cube is just right
The reflecting interface of mud stone and sandstone, lower zero phase is answered just to correspond to the reflecting interface of sandstone and mud stone, the present invention passes through 90 ° of phases
Seismic reflection main lobe is transformed into thin layer center by conversion, and the earthquake axis to be obtained according to the second seismic data cube is equal to reality
The geology axis on border.
Specifically, in embodiments of the present invention, after obtaining the second seismic data cube, to the second seismic data cube into
Row analysis, analyzes slice whens second seismic data cube etc., the plane layer attribute of the second seismic data cube is extracted, such as in data
Along layer or the plane layer attribute of interlayer, plane layer attribute includes amplitude and/or phase and/or frequency for extraction on body.
Since the information that seismic data cube includes belongs to body data, vertical resolution is relatively low, therefore, implements in the present invention
In example, before step 104, it is layered according to the drilling geology of the survey area, includes to second seismic data cube
The depth of each formation information is marked at least one formation information.Specifically, the drilling well in survey area, drilling well
Depth can be 3700m, acquire well data after drilling, can be obtained under different depth at drilling well by well data
The type of sand body structure obtains drilling geology layering, in turn, by drilling geology layering and the stratum in the second seismic data cube
Information carries out depth matching, and the depth of each formation information in the second seismic data cube is marked.Since well data has
There are the high-resolution on vertical and abundant high-frequency information, it is low in vertical upper resolution ratio so as to make up the second seismic data cube
Defect.
Further, step 104 specifically includes:According to the depth peace of each formation information in second seismic data cube
Face layer attribute predicts the type and form and each geology of geological stratification corresponding with each formation information survey area Nei
Relative position relation between layer.
Specifically, on the one hand, the formation information pair can be predicted according to the plane layer attribute of each formation information
The type and form for the practical geological stratification answered, on the other hand, according to the corresponding depth of each formation information, by believing all stratum
Breath is integrated, and can be obtained the position in survey area between the type and form and each geological stratification of all geological stratifications and be closed
System.
Citing illustrates the method for the prediction sandbody distribution of the present invention below:
After exciting artificial earthquake, the first seismic data cube in survey area is obtained, the first seismic data cube of reversion exists
The phase of quarter-wave strong point obtains the second seismic data cube.On the one hand, the drilling well in survey area, obtains according to drilling well
Drilling geology is layered, and drilling geology layering includes the first geological stratification, the second geological stratification and third geological stratification, and the first geological stratification corresponds to
Depth be earth's surface or less 1500m, the depth of the second geological stratification is earth's surface or less 1550m, and the corresponding depth of third geological stratification is
1600m below earth's surface, the first geological stratification,.On the other hand, corresponding in the second seismic data cube to include the first geological stratification information, the
Two geological stratification information and third geological stratification information, corresponding first ground of the first geological stratification information known to the second seismic data cube of analysis
Matter layer, corresponding second geological stratification of the second geological stratification information and the corresponding third geological stratification of third geological stratification information.Also, according to
The type of first geological stratification known to first geological stratification information is sandstone layer and the form of the first geological stratification.According to the second geology
The type of second geological stratification known to layer information is the form of shale layer and the second geological stratification.It can according to third geological stratification information
Know that the type of third geological stratification is similarly sandstone layer and the form of third geological stratification.But due to the second seismic data cube
Vertical resolution is relatively low, and the first geological stratification, the second geological stratification and third can not be accurately predicted according only to the second earthquake volume data
Depth location of the geological stratification in survey area, therefore, by drilling geology layering and each stratum in the second seismic data cube
Information is matched, so that it is determined that the depth of corresponding first geological stratification of the first formation information gone out in the second seismic data cube is
The depth of 1500m below earth's surface, corresponding second geological stratification of the second formation information in the second seismic data cube are earth's surface or less
The depth of 1550m, the corresponding third geological stratification of third formation information in the second seismic data cube are earth's surface or less 1600m, from
And the prediction to the sandbody distribution in survey area is realized, prediction result is accurate.
Based on same inventive concept, the embodiment of the present invention also provides a kind of device of prediction sandbody distribution, as shown in Fig. 2,
Described device includes:
Acquisition module 201, for after exciting artificial earthquake, obtaining the first seismic data cube in survey area;
Reversal block 202 obtains the second earthquake for inverting phase of first seismic data cube at 1/4 wavelength
Data volume;
Extraction module 203, the plane layer attribute for extracting second seismic data cube;
Prediction module 204, for according to the plane layer attribute, predicting the sandbody distribution of the survey area.
Further, described device further includes:
Mark module, for according to the layering of the drilling geology of the survey area, including to second seismic data cube
At least one formation information in the depth of each formation information be marked.
Further, prediction module 204 is specifically used for:
According to the depth and plane layer attribute of each formation information in second seismic data cube, the exploration is predicted
Relative position relation in region between the type and form and each geological stratification of geological stratification corresponding with each formation information.
Further, plane stratum attribute information includes amplitude and/or phase and/or frequency.
Technical solution in above-mentioned the embodiment of the present application, at least has the following technical effect that or advantage:
Phase of the present invention by the first seismic data cube of reversion at 1/4 wavelength, obtains the second seismic data cube, to
The wave crest of second seismic data cube no longer corresponds to the interface of mud stone and sandstone, and the upper zero phase of the second seismic data cube is just right
The reflecting interface of mud stone and sandstone, lower zero phase is answered just to correspond to the reflecting interface of sandstone and mud stone, the present invention passes through 90 ° of phases
Seismic reflection main lobe is transformed into thin layer center by conversion, and the earthquake axis to be obtained according to the second seismic data cube is equal to reality
The geology axis on border.
Since the second seismic data cube has the characteristics that range broad covered area, trace horizon information are continuous in the horizontal, and
Drilling geology is layered in the longitudinal direction not only high resolution, includes also abundant high-frequency information, therefore the present invention is by the second
The spread of the sand body in survey area is predicted in conjunction with drilling geology layering on the basis of shake data volume, prediction result is accurate
Really, high resolution.
Although preferred embodiments of the present invention have been described, it is created once a person skilled in the art knows basic
Property concept, then additional changes and modifications may be made to these embodiments.So it includes excellent that the following claims are intended to be interpreted as
It selects embodiment and falls into all change and modification of the scope of the invention.
Obviously, various changes and modifications can be made to the invention without departing from essence of the invention by those skilled in the art
God and range.In this way, if these modifications and changes of the present invention belongs to the range of the claims in the present invention and its equivalent technologies
Within, then the present invention is also intended to include these modifications and variations.
Claims (4)
1. a kind of method of prediction sandbody distribution, which is characterized in that the method includes:
After exciting artificial earthquake, the first seismic data cube in survey area is obtained;
Phase of first seismic data cube at 1/4 wavelength is inverted, the second seismic data cube is obtained;
Extract the plane layer attribute of second seismic data cube;
According to the plane layer attribute, the sandbody distribution of the survey area is predicted;
Wherein, before the sandbody distribution for predicting the survey area, further include according to the plane layer attribute described:
It is layered according to the drilling geology of the survey area, at least one formation information for including to second seismic data cube
In the depth of each formation information be marked;
Wherein, described that the sandbody distribution of the survey area is predicted according to the plane layer attribute, including:
According to the depth and plane layer attribute of each formation information in second seismic data cube, the survey area is predicted
Relative position relation between the type and form and each geological stratification of interior geological stratification corresponding with each formation information.
2. the method as described in claim 1, which is characterized in that plane stratum attribute information includes amplitude and/or phase
And/or frequency.
3. a kind of device of prediction sandbody distribution, which is characterized in that described device includes:
Acquisition module, for after exciting artificial earthquake, obtaining the first seismic data cube in survey area;
Reversal block obtains the second seismic data cube for inverting phase of first seismic data cube at 1/4 wavelength;
Extraction module, the plane layer attribute for extracting second seismic data cube;
Prediction module, for according to the plane layer attribute, predicting the sandbody distribution of the survey area;
Wherein, further include:
Mark module, for according to the drilling geology of the survey area be layered, to second seismic data cube include to
The depth of each formation information is marked in a few formation information;
Wherein, the prediction module is specifically used for:
According to the depth and plane layer attribute of each formation information in second seismic data cube, the survey area is predicted
Relative position relation between the type and form and each geological stratification of interior geological stratification corresponding with each formation information.
4. device as claimed in claim 3, which is characterized in that plane stratum attribute information includes amplitude and/or phase
And/or frequency.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610431229.XA CN106199707B (en) | 2016-06-16 | 2016-06-16 | A kind of method and device of prediction sandbody distribution |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610431229.XA CN106199707B (en) | 2016-06-16 | 2016-06-16 | A kind of method and device of prediction sandbody distribution |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106199707A CN106199707A (en) | 2016-12-07 |
CN106199707B true CN106199707B (en) | 2018-09-04 |
Family
ID=57460685
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610431229.XA Active CN106199707B (en) | 2016-06-16 | 2016-06-16 | A kind of method and device of prediction sandbody distribution |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106199707B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106940451B (en) * | 2017-05-08 | 2023-08-11 | 中油(天津)国际石油勘探开发技术有限公司 | Method for predicting plastic stratum deformation |
CN116047598B (en) * | 2022-11-07 | 2023-09-19 | 北京派特杰奥科技有限公司 | Oil exploration method and system based on artificial intelligence |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6654692B1 (en) * | 2002-11-21 | 2003-11-25 | Conocophillips Company | Method of predicting rock properties from seismic data |
CN102269824B (en) * | 2010-06-02 | 2014-02-05 | 中国石油天然气集团公司 | Phase conversion processing method for wavelet of seismic data |
CN103217714B (en) * | 2013-03-26 | 2015-04-15 | 中国石油大学(华东) | Seismic sedimentology interpretation method based on frequency-scale matching |
CN104991270A (en) * | 2015-06-08 | 2015-10-21 | 中国石油化工股份有限公司 | Superimposed sand body identification method based on seismic forward modeling |
-
2016
- 2016-06-16 CN CN201610431229.XA patent/CN106199707B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN106199707A (en) | 2016-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Toksöz et al. | Attenuation of seismic waves in dry and saturated rocks: I. Laboratory measurements | |
Castellaro | The complementarity of H/V and dispersion curves | |
US20130201795A1 (en) | Fracture identification from azimuthal migrated seismic data | |
Bouchaala et al. | Azimuthal investigation of compressional seismic-wave attenuation in a fractured reservoir | |
Naseer et al. | Porosity prediction of lower cretaceous unconventional resource play, south Indus Basin, Pakistan, using the seismic spectral decomposition technique | |
Brodic et al. | Three-component seismic land streamer study of an esker architecture through S-and surface-wave imaging | |
Martorana et al. | Seismo-stratigraphic model of “La Bandita” area in the Palermo Plain (Sicily, Italy) through HVSR inversion constrained by stratigraphic data | |
Ahmad et al. | High-resolution seismic and resistivity profiling of a buried Quaternary subglacial valley: Northern Alberta, Canada | |
CN106199707B (en) | A kind of method and device of prediction sandbody distribution | |
Karim et al. | Seismic reservoir characterization using model based post-stack seismic inversion: in case of Fenchuganj gas field, Bangladesh | |
CN104375166B (en) | Method for determining seismic exploration excitation depth | |
Noguera et al. | Multi-disciplinary approach for a landing point criteria in Vaca Muerta formation | |
Benfield et al. | Extracting geologic information directly from high-resolution full-waveform inversion velocity models—A case study from offshore Trinidad | |
Rauch-Davies et al. | Diffraction imaging applied to pre-existing 3D seismic data to map fracture corridors in an unconventional play | |
Shelley et al. | Outcrop-based reservoir characterization and modeling of an Upper Mississippian mixed carbonate–siliciclastic ramp, northeastern Oklahoma | |
Nworie et al. | Paleokarst features in the aptian carbonates of the barra velha formation, santos basin, brazil | |
Almalki et al. | Efficiency of seismic attributes in detecting near-surface cavities | |
Shelley et al. | Reservoir characterization and modeling of a subsurface Meramec analog from a quarry in northeastern Oklahoma | |
Tavella et al. | Integrated characterization of unconventional upper Jurassic reservoir in northern Mexico | |
Van Oorschot et al. | Development of the Brimmond Sand Fairway | |
Walker et al. | Straining at the leash: Understanding the full potential of the deep-water, subsalt Mad Dog field, from appraisal through early production | |
Williams-Rojas et al. | Geologic controls on reservoir performance in Muspac and Catedral gas fields, southeastern Mexico | |
Imposa et al. | Borehole seismic surveys for the mechanical characterization of a calcarenite cliff in the area of Ispica (southern Sicily) | |
Ogiesoba et al. | Diffraction imaging of lithology and fluid saturation in fault zones within the Austin Chalk and Eagle Ford Shale, Maverick Basin, South Texas | |
Luke et al. | Characterizing anomalous ground for engineering applications using surface-based seismic methods |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |