CN109143348A - 3D seismic data tomography enhanced processing method - Google Patents
3D seismic data tomography enhanced processing method Download PDFInfo
- Publication number
- CN109143348A CN109143348A CN201710506012.5A CN201710506012A CN109143348A CN 109143348 A CN109143348 A CN 109143348A CN 201710506012 A CN201710506012 A CN 201710506012A CN 109143348 A CN109143348 A CN 109143348A
- Authority
- CN
- China
- Prior art keywords
- seismic data
- seismic
- cube
- tomography
- data
- 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.)
- Pending
Links
- 238000003325 tomography Methods 0.000 title claims abstract description 58
- 238000003672 processing method Methods 0.000 title claims abstract description 17
- 230000001427 coherent effect Effects 0.000 claims abstract description 63
- 238000000034 method Methods 0.000 claims abstract description 48
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 38
- 238000001914 filtration Methods 0.000 claims abstract description 29
- 230000001133 acceleration Effects 0.000 claims abstract description 28
- 238000010276 construction Methods 0.000 claims abstract description 25
- 230000008569 process Effects 0.000 claims abstract description 16
- 238000012545 processing Methods 0.000 claims abstract description 12
- 230000002708 enhancing effect Effects 0.000 claims abstract description 5
- 230000004927 fusion Effects 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 12
- 238000004458 analytical method Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 7
- 238000003384 imaging method Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 6
- 238000011161 development Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 241000500042 Eogenes Species 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 238000004088 simulation Methods 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/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/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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
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
Disclose a kind of 3D seismic data tomography enhanced processing method.This method may include: to determine target zone window scope based on seismic data cube, pre-process to seismic data cube;Seismic data cube attribute is counted, construction guide way is calculated;Based on construction guide way and seismic data cube, the Acceleration Algorithm in Seismic Coherence Cube in the window scope of target zone is calculated;Attribute based on Acceleration Algorithm in Seismic Coherence Cube and seismic data cube obtains the coherent value of each seismic data of seismic data cube, sets threshold value;It is less than the seismic data of threshold value for coherent value, it is handled using diffusing filter, obtain diffusing filter data volume, it is greater than the seismic data of threshold value for coherent value, using median filter process, median filtering data volume is obtained, fusion is carried out to diffusing filter data volume and median filtering data volume and obtains tomography enhancing data volume.The present invention combines diffusing filter processing with median filter process, improves fault recognizing precision, solves the problems, such as that section information and layer position continuity are strengthened.
Description
Technical field
The present invention relates to oil-gas explorations and tectonic geology to investigate field, disconnected more particularly, to a kind of 3D seismic data
Layer enhanced processing method.
Background technique
Eastern China Bohai gulf basin morning third generation and Song-liao basin Yanshan orogen (Fault-subsidence Stage) period, tectonic activity
Strongly, a large amount of tomographies are formed, it is oil gas that tomography, which plays important control action to sedimentation, oil-gas migration and structural trap etc.,
Emphasis describes element in exploration.Bohai gulf basin Eogene and Song-liao basin sieve dwarf-Cretaceous System rift layer buried depth are big, earthquake money
Expect that frequency, noise are relatively low, breakpoint, cross-sectional imaging effect are often bad.Especially Song-liao basin Fault-subsidence Stage, tectonic movement are living
Jump, colcanism are frequent, and clastic deposited sediments are influenced by volcanic eruption, intrusion, and seismic imaging precision reduces, and breakpoint, tomography are anti-
Feature Fuzzy is penetrated, identification is explained difficult.
To improve fault interpretation precision, generally requires to carry out Acceleration Algorithm in Seismic Coherence Cube processing, improve tomographic imaging.From nineteen ninety-five
Since Bahorich delivers one text of Acceleration Algorithm in Seismic Coherence Cube analytical technology, coherent body technique has been widely used in fault recognizing explanation,
Increasingly important role is played in structure interpretation.By development in more than 20 years, four generation coherent body skills are successively had developed
Art.
First generation coherent body technique is the coherent body technique based on cross-correlation, and the algorithm based on cross-correlation has only been used 3 genuine
Data are shaken, but the practical application effect in 3-d seismic data set is significant, than seismic properties such as seismic amplitude slices to tomography
Portray and become apparent from, be used widely in the explanation of the area three-dimensionals earthquake fault such as European North Sea.
Second generation coherent body technique mainly have, variance similar based on multiple tracks, manhatton distance waveform similarity estimation or
The coherent body technique of complex seismic trace, it is more common with first two technology, routine is formed in based on the similar coherent body technique of multiple tracks
Need to define the aperture of a room and time in speed and Tomography Velocity analysis foundation, in calculating, i.e., when three-dimensional data is analyzed
The size of window, while its orientation and inclination angle need to be also defined in each data point of data volume, increase the ground for calculating coherent body
Road quantity is shaken, increases relevant volume data calculation amount, but improve the signal-to-noise ratio of coherent body;Variance is per pass and mean-trace
Difference valuation, if all roads are all identical, the variance valuation of coherent body is zero, and similarity factor estimated value is 1, however,
Even if all seismic channel waveforms are identical, if their amplitude is different, but the phase responsibility at this time based on cross-correlation estimation
Number is also 1, that is the calculation of coherence cube based on variance only considers waveform, does not consider the other information such as amplitude, earthquake square margin
Have become the important foundation that fault interpretation or tomography automatic tracing are explained at present according to body.
Third generation coherent body technique is the coherent body technique based on characteristic value structure, and what this method needed initially to estimate inclines
Angle or azimuth, the estimation of both occurrence can be scanned realization at regular intervals, and the coherent algorithm of characteristic value structure
Inherently may be implemented to carry out inclination angle or azimuth sweep at regular intervals, first from data volume when extraction and analysis in window one
Group seismic channel data generate sample vector, sample vector analysis time t ± k it is vertical when window in, by time t depending on inclining
It in the plane that angle p and q are determined, is generated by data sample interpolation, each sample vector forms the row of new data matrix.Again by
The column of newly-generated data matrix carry out auto-correlation and cross-correlation, generate covariance matrix according to covariance matrix and calculate it
Value indicative and its corresponding eigenvector.If the waveform in all roads is all consistent in analysis window, characteristic value coherence factor is equal to
1, if waveform is consistent, but amplitude is different, and similitude coherence factor and variance type coherence factor are all not equal to 0 at this time.It is being based on
Data coherency imaging effect can be improved using complex seismic trace technology in the coherent body of characteristic value structure, in current producing scientific research
This coherent body technique obtains most commonly used application, and some maturation business interpretation softwares have had based on characteristic value structure
Calculation of coherence cube module.
The coherent body technique of a new generation is the coherent body technique based on geometric structural tensor, and this geometric structural tensor includes
The inclination angle of reflecting interface and azimuth information, when can steadily estimate in window the reflecting interface of analysis site t inclination angle and side
Parallactic angle is very suitable to the automatic interpretation and the analysis of 3-D seismics structure attribute of earthquake data volume.Firstly, defining a three-dimensional data
Seismic signal body, it has fixed an inclination angle and orientation, the directional derivative vector of seismic amplitude sampling point by three directions ground
The covariance matrix that amplitude sampling point directional derivative generates the local direction derivative of a 3x3 is shaken, by available 3 sheets of the matrix
Value indicative and its corresponding eigenvector.Using 3 characteristic values, the metric attribute for being " chaos " is defined, for describing ground
The internal structure of layer, the method can show the spatial of seismic data interrupting layer well.The coherent body technique exists at present
Both at home and abroad in the stage is studied and improved, business software is not yet formed.
The appearance and rapid development of coherent body technique substantially increase 3-D seismics interrupting layer image quality, it has also become ground
Important fault recognizing technology in shake explanation.Numerous research and development, the coherent body technique article that substantial amounts have been write using personnel.
Yuan Shujin has combed the generation and development course of Acceleration Algorithm in Seismic Coherence Cube Technique, describes every generation coherent body technique in detail
Principle and algorithm, applicable elements, briefly list application example, effect, looked forward to the development trend of Acceleration Algorithm in Seismic Coherence Cube Technique,
Analyze the application conditions of coherent body technique in the factor and seismic interpretation that influence coherent body resolution ratio.
Zhang Junhua et al. proposes improved coherence algorithm on the basis of three generations's coherent body technique theory analysis, and
Corresponding software is developed, in the area Dagang Oilfield TJH, the coherent body technique of application enhancements algorithm has studied target zone stratum knot
Structure and reservoir change, and feature tomography.
Wang Xiwen et al. has researched and proposed the computation of seismic coherence cube based on wavelet transformation, proposes with simulation seismic wavelet
Wavelet function wavelet transformation, obtain frequency dividing instantaneous phase, then calculate the algorithm 1 of coherent body and obtained with wavelet transformation
Real, imaginary part calculates the algorithm 2 of coherent body, in order to protrude craven fault feature, is reconstructed with the coherent body that Color seperation grating goes out.It is practical
Material computation shows that 2 anti-noise ability of algorithm is strong, and than algorithm 2, the small advertisement in craven fault is portrayed becomes apparent from algorithm 1.
But being met certain condition using needs for coherent body technique, broken and distribution of strata complex area is being constructed, such as
Fruit seismic data quality is not again high, and the application effect of coherent body technique will give a discount, it is difficult to clearly symbolize fault feature.?
In the Cretaceous System rift layer seismic interpretation of Song-liao basin south, due to target zone bury depth, occurrence variation greatly, volcanic rock divides extensively
The factors such as cloth, multi-stage tectonic movement superposition and mature fault, seismic data noise is relatively low, target zone reflection wave groups
And section wave characteristic is unclear.Using third generation coherent body technique, tomographic imaging quality is improved, but low by earthquake signal-to-noise ratio,
The unfavorable conditions such as attitude of stratum is changeable, lithology cross directional variations are big restrict, coherent body data outage Surface wave characteristic is fully aware of,
It is prominent, moreover, the not instead of fault information of the irrelevant reflection in part, lithology cross directional variations produce by comparing with seismic data
Raw seismic reflection dies down or discontinuously, in this way, relevant volume data just includes stronger multi-solution and uncertain, it is difficult to improve
Fault recognizing Explanation Accuracy.Therefore, it is necessary to develop a kind of tomography intensive treatment suitable for low signal-to-noise ratio 3D seismic data
Method, to improve cross-sectional imaging and fault recognizing Explanation Accuracy.
The information for being disclosed in background of invention part is merely intended to deepen the reason to general background technique of the invention
Solution, and it is known to those skilled in the art existing to be not construed as recognizing or imply that the information is constituted in any form
Technology.
Summary of the invention
The invention proposes a kind of 3D seismic data tomography enhanced processing methods, can be by diffusing filter processing in
Value filtering processing combines, and improves fault recognizing precision, solves the problems, such as that section information and layer position continuity are strengthened.
The invention proposes a kind of 3D seismic data tomography enhanced processing methods.The method may include: based on ground
Data volume is shaken, determines the window scope of target zone, and pre-process to the seismic data cube;Count the seismic data cube
Attribute, and then calculate construction guide way;Based on the construction guide way and the seismic data cube, the target zone is calculated
Acceleration Algorithm in Seismic Coherence Cube in window scope;Attribute based on the Acceleration Algorithm in Seismic Coherence Cube Yu the seismic data cube obtains seismic data
The corresponding coherent value of each seismic data in body, and then set the relevant threshold value of tomography;It is less than for the coherent value described
The seismic data of the relevant threshold value of tomography is handled the seismic data using diffusing filter, obtains diffusing filter data volume, needle
It is greater than the seismic data of the relevant threshold value of the tomography to the coherent value, median filter process is used to the seismic data,
Median filtering data volume is obtained, fusion is carried out to the diffusing filter data volume and median filtering data volume and obtains tomography enhancing number
According to body.
Preferably, carrying out pretreatment to the seismic data cube includes: to remove relevant noise and random noise.
Preferably, the attribute of the seismic data cube includes: the stratigraphic dip of the waveform of the seismic data cube, axis of reflection
With azimuth.
Preferably, calculating construction guide way include: based on the stratigraphic dip and azimuthal angle calculation inclination angle body,
As the construction guide way.
Preferably, calculating the Acceleration Algorithm in Seismic Coherence Cube in the window scope of the target zone includes: based on the construction guide way
With the seismic data cube, the similitude of waveform between the neighboring track in the window scope of target zone is calculated, obtains the earthquake
Coherent body.
It preferably, include: to be estimated describedly by gradient-structure tensor using diffusing filter processing to the seismic data
Shake the discontinuity surface of data volume;Judge with the presence or absence of discontinuity surface in window scope, if there is no discontinuity surface, then along institute
It states stratigraphic dip and azimuth and carries out smooth, acquisition diffusing filter data volume to the seismic data cube.
It preferably, include: according to the stratigraphic dip, along pendage using median filter process to the seismic data
Direction carries out median filtering to the seismic data cube, obtains median filtering data volume.
Preferably, determine the seismic profile breakpoint in the seismic data cube, read in Acceleration Algorithm in Seismic Coherence Cube with seismic profile
The corresponding coherent value of breakpoint sets the coherent value as the relevant threshold value of tomography.
The beneficial effects of the present invention are: for guarantor's width with higher, fidelity, target zone has certain signal-to-noise ratio
Original earthquake data improves fault recognizing precision, realizes by combining diffusing filter processing with median filter process
Fine tomography, layer position are explained and construction is implemented, and establish solid foundation for evaluation of trap and Efficient Exploration exploitation.
Methods and apparatus of the present invention has other characteristics and advantages, these characteristics and advantages are attached from what is be incorporated herein
It will be apparent in figure and subsequent specific embodiment, or will be in the attached drawing being incorporated herein and subsequent specific reality
It applies in mode and is stated in detail, the drawings and the detailed description together serve to explain specific principles of the invention.
Detailed description of the invention
Exemplary embodiment of the present is described in more detail in conjunction with the accompanying drawings, of the invention is above-mentioned and other
Purpose, feature and advantage will be apparent, wherein in exemplary embodiments of the present invention, identical reference label is usual
Represent same parts.
Fig. 1 shows the flow chart of the step of 3D seismic data tomography enhanced processing method according to the present invention.
Fig. 2 a and Fig. 2 b respectively illustrate the section and seismic facies of seismic data cube according to an embodiment of the invention
The schematic diagram of the section of stem body.
Fig. 3 a and Fig. 3 b respectively illustrate the section of seismic data cube according to an embodiment of the invention and diffusion is filtered
The schematic diagram of the section of wave data volume.
Fig. 4 a and Fig. 4 b show the section and median filtering number of seismic data cube according to an embodiment of the invention
According to the schematic diagram of the section of body.
Fig. 5 a and Fig. 5 b respectively illustrate the section of seismic data cube according to an embodiment of the invention and tomography increases
The schematic diagram of the section of strong data volume.
Fig. 6 a and Fig. 6 b, which respectively illustrate Fault profile according to an embodiment of the invention and tomography, enhances data volume
Section schematic diagram.
Fig. 7 a and Fig. 7 b respectively illustrate the slice of seismic data cube according to an embodiment of the invention and tomography increases
The schematic diagram of the slice of strong data volume.
Specific embodiment
The present invention will be described in more detail below with reference to accompanying drawings.Although showing the preferred embodiment of the present invention in attached drawing,
However, it is to be appreciated that may be realized in various forms the present invention and should not be limited by the embodiments set forth herein.On the contrary, providing
These embodiments are of the invention more thorough and complete in order to make, and can will fully convey the scope of the invention to ability
The technical staff in domain.
Fig. 1 shows the flow chart of the step of 3D seismic data tomography enhanced processing method according to the present invention.
3D seismic data tomography enhanced processing method according to the present invention may include:
Step 101, it is based on seismic data cube, determines the window scope of target zone, and pre-process to seismic data cube;
The window scope of target zone is determined according to synthetic seismogram calibration and target zone top, bottom reflecting layer primary explanation achievement;
Step 102, the attribute of seismic data cube is counted, and then calculates construction guide way;
Step 103, based on construction guide way and seismic data cube, the seismic coherence in the window scope of target zone is calculated
Body;Calculate the similitude of waveform between the neighboring track in the window scope of target zone, obtained data volume i.e. Acceleration Algorithm in Seismic Coherence Cube;
Step 104, the attribute based on Acceleration Algorithm in Seismic Coherence Cube and seismic data cube obtains each earthquake number in seismic data cube
According to corresponding coherent value, and then set the relevant threshold value of tomography;Assuming that completely relevant, i.e. seismic data waveform, inclination angle and orientation
The attributes such as angle are completely the same, and defining coherent value is 1, then, the smaller expression incoherence of coherent value is stronger;
Step 105, the seismic data of the relevant threshold value of tomography is less than for coherent value, diffusing filter is used to seismic data
Processing obtains diffusing filter data volume, and the seismic data of the relevant threshold value of tomography is greater than for coherent value, is used to seismic data
Median filter process obtains median filtering data volume, carries out fusion acquisition to diffusing filter data volume and median filtering data volume
Tomography enhances data volume;Median filter process increases the continuity of lineups, and diffusing filter processing highlights fault information, makes
Near Faults lineups sharpen, and for the ease of structure interpretation, need seismic data interrupting layer development region lineups to sharpen, tomography
Agensis region lineups are continuous;It is concerned with threshold value according to tomography, coherent value is less than tomography and is concerned with the data of threshold value, to earthquake
Data are handled using diffusing filter, obtain diffusing filter data volume, coherent value is greater than the data of the relevant threshold value of tomography, to earthquake
Data use median filter process, obtain median filtering data volume, and then merge and obtain tomography enhancing data volume, can solve simultaneously
The problem of section information and layer position continuity are strengthened.
In one example, carrying out pretreatment to seismic data cube includes: to remove relevant noise and random noise.For example,
The relevant noise and random noise of seismic data cube can be removed by median filtering and bandpass filtering, median filtering generally can be absolutely
Noise peak is prevented, local smoothing method is carried out to data sequence, but sometimes cannot be completely smooth to the small amplitude occurred at random, because
This needs to do after median filtering bandpass filtering, and the small amplitude in smooth pulse improves target zone signal-to-noise ratio.
In one example, the attribute of seismic data cube include: the waveform of seismic data cube, axis of reflection stratigraphic dip with
Azimuth.
In one example, calculating construction guide way includes: based on stratigraphic dip and azimuthal angle calculation inclination angle body, as structure
Make guide way.It is successively counted in the main profile direction of seismic data cube with cross-track direction by stratigraphic dip with azimuth
It calculates, obtains inclination angle body, as construction guide way.
In one example, calculate the Acceleration Algorithm in Seismic Coherence Cube in the window scope of target zone include: based on construction guide way with
Seismic data cube calculates the similitude of waveform between the neighboring track in the window scope of target zone, obtains Acceleration Algorithm in Seismic Coherence Cube.It is based on
Guide way and seismic data cube are constructed, the similitude of waveform between the neighboring track in the window scope of target zone is calculated, according to phase
It is Acceleration Algorithm in Seismic Coherence Cube like the data volume that property obtains, reduces irrelevant caused by stratigraphic dip, azimuthal variation, highlights
Section generates irrelevant.
In one example, it determines the seismic profile breakpoint in seismic data cube, reads in Acceleration Algorithm in Seismic Coherence Cube and cutd open with earthquake
Breakpoint corresponding coherent value in face sets coherent value as the relevant threshold value of tomography.Position is explained according to seismic profile breakpoint, definitely
The seismic profile breakpoint in data volume is shaken, since the data of seismic data cube and Acceleration Algorithm in Seismic Coherence Cube are correspondingly, can to read
Coherent value corresponding with seismic profile breakpoint in Acceleration Algorithm in Seismic Coherence Cube is taken, sets the coherent value as the relevant threshold value of tomography.
It in one example, include: that gradient knot is passed through based on construction guide way using diffusing filter processing to seismic data
The discontinuity surface of structure tensor estimation seismic data cube;Judge with the presence or absence of discontinuity surface in window scope, if there is no not
Continuous surface then carries out median filter smoothness of image to seismic data cube along stratigraphic dip and azimuth, obtains diffusing filter data volume.One
As random noise and relevant noise in seismic data can be removed by 1-5 iteration, enhance target zone it is irrelevant at
Image quality amount.
It in one example, include: to be based on construction guide way and stratigraphic dip using median filter process to seismic data,
Median filtering is carried out to seismic data cube along pendage direction, obtains median filtering data volume.By along pendage direction
Median filtering is carried out to seismic data cube, the amplitude inconsistency of seismic reflection lineups caused by non-tomography factor can be reduced
Weaken with continuity, improves the continuity and wave group feature of normal stratum reflection.
The present invention has the original earthquake data of certain signal-to-noise ratio for guarantor's width with higher, fidelity, target zone, leads to
It crosses and combines diffusing filter processing with median filter process, improve fault recognizing precision, realize fine tomography, layer position solution
It releases and constructs and implement, establish solid foundation for evaluation of trap and Efficient Exploration exploitation.
Using example
A concrete application example is given below in the scheme and its effect of the embodiment of the present invention for ease of understanding.This field
It should be understood to the one skilled in the art that the example is only for the purposes of understanding the present invention, any detail is not intended to be limited in any way
The system present invention.
Fig. 2 a and Fig. 2 b show the section of seismic data cube according to an embodiment of the invention and Acceleration Algorithm in Seismic Coherence Cube
Section schematic diagram.
Based on seismic data cube, the window scope of target zone is determined, earthquake number is removed by median filtering and bandpass filtering
According to the relevant noise and random noise of body.The attribute for counting seismic data cube, including waveform, the stratigraphic dip of axis of reflection and orientation
Angle is successively calculated with cross-track direction in the main profile direction of seismic data cube with azimuth by stratigraphic dip, is obtained
Inclination angle body is obtained, as construction guide way.Based on construction guide way and seismic data cube, the phase in the window scope of target zone is calculated
The similitude of waveform between neighboring trace, is Acceleration Algorithm in Seismic Coherence Cube according to the data volume that similitude obtains, as shown in Figure 2 b, with Fig. 2 a's
Seismic data cube Profile Correlation, Acceleration Algorithm in Seismic Coherence Cube reduce irrelevant caused by stratigraphic dip, Orientation differences, highlight section
What is generated is irrelevant.
Fig. 3 a and Fig. 3 b show the section and diffusing filter number of seismic data cube according to an embodiment of the invention
According to the schematic diagram of the section of body.
Fig. 4 a and Fig. 4 b show the section and median filtering number of seismic data cube according to an embodiment of the invention
According to the schematic diagram of the section of body.
Fig. 5 a and Figure 5b shows that the section of seismic data cube according to an embodiment of the invention and tomographies to enhance number
According to the schematic diagram of the section of body.
Fig. 6 a and Fig. 6 b, which show Fault profile according to an embodiment of the invention and tomography, enhances cuing open for data volume
The contrast schematic diagram in face.
Fig. 7 a and Fig. 7 b, which show the slice of seismic data cube according to an embodiment of the invention and tomography, enhances number
According to the schematic diagram of the slice of body.
Attribute based on Acceleration Algorithm in Seismic Coherence Cube and seismic data cube, each seismic data obtained in seismic data cube are corresponding
Coherent value, it is assumed that completely relevant, i.e. the attributes such as seismic data waveform, inclination angle and azimuth are completely the same, and defining coherent value is 1,
So, the smaller expression incoherence of coherent value is stronger.Position is explained according to seismic profile breakpoint, determines the ground in seismic data cube
Section breakpoint is shaken, since the data of seismic data cube and Acceleration Algorithm in Seismic Coherence Cube are correspondingly, can to read in Acceleration Algorithm in Seismic Coherence Cube
Coherent value corresponding with seismic profile breakpoint is 0.5, then sets 0.5 for the relevant threshold value of tomography, for coherent value less than 0.5
Data are handled using diffusing filter, are calculated construction guide way according to seismic data cube, are counted the stratigraphic dip and orientation of axis of reflection
Angle, and then based on construction guide way, the discontinuity surface of seismic data cube is estimated by gradient-structure tensor;Judge in window scope
It is interior then seismic data cube to be carried out along stratigraphic dip and azimuth flat with the presence or absence of discontinuity surface if there is no discontinuity surface
It is sliding, diffusing filter data volume is obtained, as shown in Figure 3b, the seismic data cube Profile Correlation with Fig. 3 a, diffusing filter data volume is sharp
Change hiding fault information, improves tomographic imaging precision;It is greater than 0.5 data for coherent value, at median filtering
Reason is carried out median filtering to seismic data cube along pendage direction, is obtained intermediate value filter based on construction guide way and stratigraphic dip
Wave data volume, as shown in Figure 4 b, the seismic data cube Profile Correlation with Fig. 4 a, median filtering data volume reduce non-tomography not phase
Stemness enhances the continuity and consistency of wave group, and then merges and obtain tomography enhancing data volume, such as Fig. 5 b, Fig. 6 b and Fig. 7 b
It is shown, it is compared with Fig. 5 a, Fig. 6 a and Fig. 7 a, the present invention improves fault tectonic explanation and mapping precision.
In conclusion the present invention, for guarantor's width with higher, fidelity, target zone has certain signal-to-noise ratio primitively
Data are shaken, by combining diffusing filter processing with median filter process, fault recognizing precision is improved, realizes fine disconnected
Layer, layer position are explained and construction is implemented, and establish solid foundation for evaluation of trap and Efficient Exploration exploitation.
It will be understood by those skilled in the art that above to the purpose of the description of the embodiment of the present invention only for illustratively saying
The beneficial effect of bright the embodiment of the present invention is not intended to limit embodiments of the invention to given any example.
Various embodiments of the present invention are described above, above description is exemplary, and non-exclusive, and
It is not limited to disclosed each embodiment.Without departing from the scope and spirit of illustrated each embodiment, for this skill
Many modifications and changes are obvious for the those of ordinary skill in art field.The selection of term used herein, purport
In the principle, practical application or improvement to the technology in market for best explaining each embodiment, or make the art
Other those of ordinary skill can understand each embodiment disclosed herein.
Claims (8)
1. a kind of 3D seismic data tomography enhanced processing method, comprising:
Based on seismic data cube, the window scope of target zone is determined, and the seismic data cube is pre-processed;
The attribute of the seismic data cube is counted, and then calculates construction guide way;
Based on the construction guide way and the seismic data cube, the seismic coherence in the window scope of the target zone is calculated
Body;
Attribute based on the Acceleration Algorithm in Seismic Coherence Cube Yu the seismic data cube obtains each seismic data pair in seismic data cube
The coherent value answered, and then set the relevant threshold value of tomography;
It is less than the seismic data of the relevant threshold value of the tomography for the coherent value, diffusing filter is used to the seismic data
Processing obtains diffusing filter data volume, the seismic data of the relevant threshold value of the tomography is greater than for the coherent value, to described
Seismic data uses median filter process, median filtering data volume is obtained, to the diffusing filter data volume and median filtering number
Fusion, which is carried out, according to body obtains tomography enhancing data volume.
2. 3D seismic data tomography enhanced processing method according to claim 1, wherein to the seismic data cube into
Row pretreatment includes: to remove relevant noise and random noise.
3. 3D seismic data tomography enhanced processing method according to claim 1, wherein the category of the seismic data cube
Property includes: the waveform of the seismic data cube, the stratigraphic dip of axis of reflection and azimuth.
4. 3D seismic data tomography enhanced processing method according to claim 3, wherein the calculating constructs guide way
It include: based on the stratigraphic dip and azimuthal angle calculation inclination angle body, as the construction guide way.
5. 3D seismic data tomography enhanced processing method according to claim 1, wherein calculate the target zone when
Acceleration Algorithm in Seismic Coherence Cube within the scope of window includes: to calculate the when window of target zone based on the construction guide way and the seismic data cube
The similitude of waveform between neighboring track in range, obtains the Acceleration Algorithm in Seismic Coherence Cube.
6. 3D seismic data tomography enhanced processing method according to claim 3, wherein used to the seismic data
Diffusing filter is handled
The discontinuity surface of the seismic data cube is estimated by gradient-structure tensor;
Judge with the presence or absence of discontinuity surface in window scope, if there is no discontinuity surface, then along the stratigraphic dip and side
Parallactic angle carries out smooth, acquisition diffusing filter data volume to the seismic data cube.
7. 3D seismic data tomography enhanced processing method according to claim 3, wherein used to the seismic data
Median filter process includes:
According to the stratigraphic dip, median filtering is carried out to the seismic data cube along pendage direction, obtains median filtering
Data volume.
8. 3D seismic data tomography enhanced processing method according to claim 1, wherein determine the seismic data cube
In seismic profile breakpoint, read Acceleration Algorithm in Seismic Coherence Cube in coherent value corresponding with seismic profile breakpoint, set the coherent value as
Tomography is concerned with threshold value.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710506012.5A CN109143348A (en) | 2017-06-28 | 2017-06-28 | 3D seismic data tomography enhanced processing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710506012.5A CN109143348A (en) | 2017-06-28 | 2017-06-28 | 3D seismic data tomography enhanced processing method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109143348A true CN109143348A (en) | 2019-01-04 |
Family
ID=64805407
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710506012.5A Pending CN109143348A (en) | 2017-06-28 | 2017-06-28 | 3D seismic data tomography enhanced processing method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109143348A (en) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109884708A (en) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | Lower Paleozoic buried hill reservoir seismic prediction method based on construction Steerable filter |
CN111142158A (en) * | 2020-01-15 | 2020-05-12 | 西安石油大学 | Calibration method for improving seismic interpretation precision of fault-order structure |
CN111257936A (en) * | 2020-02-29 | 2020-06-09 | 中国石油大学(华东) | Fracture characteristic re-enhancement method based on FEF effective signal compensation |
CN111582114A (en) * | 2020-04-29 | 2020-08-25 | 中国海洋石油集团有限公司 | Seismic fault identification method, device, equipment and storage medium |
CN111751871A (en) * | 2019-03-29 | 2020-10-09 | 中国石油天然气集团有限公司 | Three-dimensional seismic edge layer coherent slice processing method and device and electronic equipment |
CN111796323A (en) * | 2019-04-09 | 2020-10-20 | 中国石油化工股份有限公司 | Method and system for judging sliding fracture boundary and main section |
CN111929726A (en) * | 2019-09-04 | 2020-11-13 | 中国石油天然气股份有限公司 | Seismic coherent data volume processing method and device |
CN112114358A (en) * | 2019-06-20 | 2020-12-22 | 中国石油化工股份有限公司 | Underground volcanic channel identification method based on three-dimensional seismic data representation |
CN112213777A (en) * | 2019-07-12 | 2021-01-12 | 中国石油化工股份有限公司 | Geosteering phase interpretation method and device |
CN112394394A (en) * | 2019-08-13 | 2021-02-23 | 中国石油化工股份有限公司 | Identification method of ultra-deep volcanic rock |
CN112444880A (en) * | 2019-08-30 | 2021-03-05 | 中国石油化工股份有限公司 | Fast filtering method for suppressing clutter and storage medium |
CN113093274A (en) * | 2020-01-08 | 2021-07-09 | 中国石油天然气股份有限公司 | Method, device, terminal and storage medium for low-level sequence fault recognition |
CN113325470A (en) * | 2020-05-28 | 2021-08-31 | 中国海洋石油集团有限公司 | Self-adaptive multi-azimuth coherent fusion method |
CN114137613A (en) * | 2020-09-03 | 2022-03-04 | 中国石油化工股份有限公司 | Stratum fracture identification method and system, storage medium and electronic equipment |
CN114252913A (en) * | 2020-09-25 | 2022-03-29 | 中国石油天然气股份有限公司 | Method and device for identifying plane fault information |
CN114355449A (en) * | 2022-01-05 | 2022-04-15 | 电子科技大学 | Vector median constrained structure-oriented three-dimensional seismic image enhancement method |
CN113093274B (en) * | 2020-01-08 | 2024-05-28 | 中国石油天然气股份有限公司 | Method, device, terminal and storage medium for identifying low-order faults |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675902A (en) * | 2012-09-07 | 2014-03-26 | 中国石油化工股份有限公司 | Optimal direction edge monitoring method |
CN104459801A (en) * | 2014-12-10 | 2015-03-25 | 中国石油天然气集团公司 | Coherence enhancement processing method used for recognizing fault |
CN104570096A (en) * | 2015-01-22 | 2015-04-29 | 中国海洋石油总公司 | Granite buried hill top surface fault interpretation method |
-
2017
- 2017-06-28 CN CN201710506012.5A patent/CN109143348A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675902A (en) * | 2012-09-07 | 2014-03-26 | 中国石油化工股份有限公司 | Optimal direction edge monitoring method |
CN104459801A (en) * | 2014-12-10 | 2015-03-25 | 中国石油天然气集团公司 | Coherence enhancement processing method used for recognizing fault |
CN104570096A (en) * | 2015-01-22 | 2015-04-29 | 中国海洋石油总公司 | Granite buried hill top surface fault interpretation method |
Non-Patent Citations (2)
Title |
---|
张少青 等: "滤波技术在煤田断层精细检测中的应用", 《物探化探计算技术》 * |
田涛 等: "基于地层倾角的相干体技术——以渤海湾盆地A油田为例", 《中国石油学会2017年物探技术研讨会论文集》 * |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109884708A (en) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | Lower Paleozoic buried hill reservoir seismic prediction method based on construction Steerable filter |
CN109884708B (en) * | 2019-03-20 | 2021-09-07 | 中国石油化工股份有限公司 | Construction oriented filtering-based ancient boundary buried hill reservoir earthquake prediction method |
CN111751871B (en) * | 2019-03-29 | 2023-02-07 | 中国石油天然气集团有限公司 | Three-dimensional seismic edge layer coherent slice processing method and device and electronic equipment |
CN111751871A (en) * | 2019-03-29 | 2020-10-09 | 中国石油天然气集团有限公司 | Three-dimensional seismic edge layer coherent slice processing method and device and electronic equipment |
CN111796323A (en) * | 2019-04-09 | 2020-10-20 | 中国石油化工股份有限公司 | Method and system for judging sliding fracture boundary and main section |
CN112114358A (en) * | 2019-06-20 | 2020-12-22 | 中国石油化工股份有限公司 | Underground volcanic channel identification method based on three-dimensional seismic data representation |
CN112213777A (en) * | 2019-07-12 | 2021-01-12 | 中国石油化工股份有限公司 | Geosteering phase interpretation method and device |
CN112394394A (en) * | 2019-08-13 | 2021-02-23 | 中国石油化工股份有限公司 | Identification method of ultra-deep volcanic rock |
CN112444880A (en) * | 2019-08-30 | 2021-03-05 | 中国石油化工股份有限公司 | Fast filtering method for suppressing clutter and storage medium |
CN111929726B (en) * | 2019-09-04 | 2023-02-10 | 中国石油天然气股份有限公司 | Seismic coherent data volume processing method and device |
CN111929726A (en) * | 2019-09-04 | 2020-11-13 | 中国石油天然气股份有限公司 | Seismic coherent data volume processing method and device |
CN113093274A (en) * | 2020-01-08 | 2021-07-09 | 中国石油天然气股份有限公司 | Method, device, terminal and storage medium for low-level sequence fault recognition |
CN113093274B (en) * | 2020-01-08 | 2024-05-28 | 中国石油天然气股份有限公司 | Method, device, terminal and storage medium for identifying low-order faults |
CN111142158A (en) * | 2020-01-15 | 2020-05-12 | 西安石油大学 | Calibration method for improving seismic interpretation precision of fault-order structure |
CN111257936B (en) * | 2020-02-29 | 2021-06-01 | 中国石油大学(华东) | Fracture characteristic re-enhancement method based on FEF effective signal compensation |
CN111257936A (en) * | 2020-02-29 | 2020-06-09 | 中国石油大学(华东) | Fracture characteristic re-enhancement method based on FEF effective signal compensation |
CN111582114A (en) * | 2020-04-29 | 2020-08-25 | 中国海洋石油集团有限公司 | Seismic fault identification method, device, equipment and storage medium |
CN111582114B (en) * | 2020-04-29 | 2023-09-01 | 中国海洋石油集团有限公司 | Earthquake fault identification method, device, equipment and storage medium |
CN113325470A (en) * | 2020-05-28 | 2021-08-31 | 中国海洋石油集团有限公司 | Self-adaptive multi-azimuth coherent fusion method |
CN113325470B (en) * | 2020-05-28 | 2022-08-12 | 中国海洋石油集团有限公司 | Self-adaptive multi-azimuth coherent fusion method |
CN114137613A (en) * | 2020-09-03 | 2022-03-04 | 中国石油化工股份有限公司 | Stratum fracture identification method and system, storage medium and electronic equipment |
CN114137613B (en) * | 2020-09-03 | 2024-05-17 | 中国石油化工股份有限公司 | Formation fracture identification method, system, storage medium and electronic equipment |
CN114252913A (en) * | 2020-09-25 | 2022-03-29 | 中国石油天然气股份有限公司 | Method and device for identifying plane fault information |
CN114355449B (en) * | 2022-01-05 | 2023-04-25 | 电子科技大学 | Structure-oriented three-dimensional seismic image enhancement method based on vector median constraint |
CN114355449A (en) * | 2022-01-05 | 2022-04-15 | 电子科技大学 | Vector median constrained structure-oriented three-dimensional seismic image enhancement method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109143348A (en) | 3D seismic data tomography enhanced processing method | |
Faleide et al. | Impacts of seismic resolution on fault interpretation: Insights from seismic modelling | |
Sava | Micro-earthquake monitoring with sparsely sampled data | |
WO2008111037A2 (en) | A method for identifying and analyzing faults/fractures using reflected and diffracted waves | |
CN104360388B (en) | A kind of 3 D seismic observation system evaluation method | |
CN104268412B (en) | A kind of angle gathers ray tomography migration velocity analysis method and device | |
Justice et al. | Multidimensional attribute analysis and pattern recognition for seismic interpretation | |
CN109001813A (en) | A kind of method, apparatus and system of multiple suppression | |
US6925387B2 (en) | Method and apparatus for kinematically linking multiple seismic domains | |
CN102565852B (en) | Angle domain pre-stack offset data processing method aiming to detect oil-gas-bearing property of reservoir | |
AU2014200562B2 (en) | Systems and methods for detecting swell noise in a seismic gather | |
Chen et al. | A compact program for 3D passive seismic source‐location imaging | |
Wei et al. | Deep learning-based P-and S-wave separation for multicomponent vertical seismic profiling | |
EA030770B1 (en) | System and method for seismic adaptive optics | |
CN105093300B (en) | A kind of boundary recognition of geological body method and device | |
CN109116419A (en) | Based on the geologic anomaly conformation identification method and device that frequency-shared phase guiding is relevant | |
Carpentier et al. | Semiautomated suppression of above-surface diffractions in GPR data | |
Rybakov et al. | 3-D gravity and magnetic interpretation for the Haifa Bay area (Israel) | |
CN105954800A (en) | Method for utilizing earthquake texture characteristics to identify minor fault | |
CN109782344A (en) | Depositional sequence Boundary Recognition method and device | |
Faleide et al. | Exploring seismic detection and resolution thresholds of fault zones and gas seeps in the shallow subsurface using seismic modelling | |
CN109143398A (en) | A kind of modeling method of automatic mesh chromatography Depth Domain speed | |
Murty et al. | Mapping the sub-trappean Mesozoic sediments in the western part of Narmada–Tapti region of Deccan Volcanic Province, India | |
Qin | Full-waveform inversion of ground-penetrating radar data and its indirect joint petrophysical inversion with shallow-seismic data | |
CN109490952A (en) | Acceleration Algorithm in Seismic Coherence Cube analysis method and system |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190104 |