CN109884693A - Adaptively move towards normal-moveout spectrum acquiring method and system - Google Patents
Adaptively move towards normal-moveout spectrum acquiring method and system Download PDFInfo
- Publication number
- CN109884693A CN109884693A CN201910049867.9A CN201910049867A CN109884693A CN 109884693 A CN109884693 A CN 109884693A CN 201910049867 A CN201910049867 A CN 201910049867A CN 109884693 A CN109884693 A CN 109884693A
- Authority
- CN
- China
- Prior art keywords
- normal
- seismic
- spectrum
- moveout spectrum
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Normal-moveout spectrum acquiring method is adaptively moved towards, first step obtains seismic data, denoising；Second step carries out velocity analysis to the seismic data along earthquake main profile direction, obtains the first root mean sequare velocity；Third step forms normal-moveout spectrum and superposition seismic profile according to first root mean sequare velocity；Four steps establishes migration velocity field, and the superposition seismic profile or common midpoint gather are deviated, offset seismic profile is obtained；5th step analyzes the offset seismic profile and obtains obliquity structure coefficient, and the obliquity structure coefficient refers to inclination angle corresponding to each single-point of obliquity structure；6th step sets tilt threshold, the relationship of the slope coefficient and the tilt threshold is worked as according to the slope coefficient, it is determined whether turn azimuth.The present invention conveniently judges whether to need to carry out seismic data to turn azimuth processing；It is more accurate to make to turn azimuth processing；Obtain the channel set and supervelocity spectrum of high-resolution high s/n ratio.
Description
Technical field
This patent disclosure relates generally to field of seismic exploration, adaptively move towards normal-moveout spectrum acquiring method more particularly, to a kind of
And system.
Background technique
Complex geological structure seismic imaging is more demanding to velocity accuracy, seeks super road by correlation velocity spectrum calculation method
Collection and normal-moveout spectrum, can effectively improve normal-moveout spectrum quality.Existing seismic processes is generally carried out along acquisition azimuth direction, speed
Degree spectral method is by forming channel set and corresponding speed spectrum along the azimuth direction or the vertical azimuth directional combination mode, suddenly
The influence that channel set and normal-moveout spectrum are combined depending on the directional characteristic and complex geological structure of formation velocity.
It is disclosed application No. is 201210017112.9 patents and moves towards seismic velocity analysis methods along geological structure, including is true
Geological structure trend in constant speed degree analysis area keeps the distribution of velocity analysis face element consistent with geological structure trend, to single speed
Analysis site is spent centered on the velocity analysis point, is inclined to along geological structure and increases trace gather number, chooses best ellipsoid member combination
Minor axis length is moved towards to increase trace gather number along the geological structure, chooses best ellipsoid member combination long axis length.
But geological structure is moved towards to increase the difficulty to seism processing without differentiation in the prior art, it is small to incline
Corner structure influences very little to the signal-to-noise ratio of normal-moveout spectrum and resolution ratio, or is not necessarily to carry out turning azimuth processing to seismic data
When, it still will do it and turn azimuth processing, the introduction of extra data will affect signal-to-noise ratio and the resolution of normal-moveout spectrum to a certain extent
Rate.
Summary of the invention
The problem to be solved in the present invention is to provide one kind and adaptively moves towards normal-moveout spectrum acquiring method and system, distinguishes geology
Influence of the construction trend to speed spectral resolution and signal-to-noise ratio, seeks resolution ratio and the higher normal-moveout spectrum of signal-to-noise ratio and channel set.
The present invention provides adaptively move towards normal-moveout spectrum acquiring method, wherein including first step S1, obtains earthquake money
Material carries out denoising；Second step S2 carries out velocity analysis to the seismic data along earthquake main profile direction, obtains the
One root mean sequare velocity；Third step S3 forms normal-moveout spectrum and superposition seismic profile according to first root mean sequare velocity；4th step
Rapid S4, establishes migration velocity field, and the superposition seismic profile or common midpoint gather are deviated, and obtains offset earthquake and cuts open
Face；5th step S5 analyzes the offset seismic profile and obtains obliquity structure coefficient, and the obliquity structure coefficient refers to obliquity structure
Inclination angle corresponding to each single-point；6th step S6 sets tilt threshold, when the inclination angle overall coefficient is less than the tilt threshold
When, normal-moveout spectrum is formed according to the first root mean sequare velocity；When the obliquity structure coefficient is greater than the tilt threshold, incline along described
Angular direction redefines the corresponding root mean sequare velocity of obliquity structure, obtains the second root mean sequare velocity, according to second root mean square
Speed redefines normal-moveout spectrum.
According to embodiment of the present invention, the method includes determining seismic processing pitch region geological structure
Trend；The amendment obliquity structure coefficient is moved towards according to seismic processing pitch region geological structure.
According to embodiment of the present invention, the tilt threshold is obtained by single-spot testing: being obtained according to earthquake information
Obtain seismic profile and dip angle attribute section；Two or more test points are chosen from the obliquity structure region of seismic profile；To each
Test point forms single-point normal-moveout spectrum and corresponding trace gather along earthquake main profile direction；To each test point along earthquake main profile direction shape
Trace gather is composed and corresponded at supervelocity；The single-point normal-moveout spectrum and corresponding trace gather and supervelocity spectrum and corresponding trace gather are analyzed, is determined full
The inclination maximum that sufficient channel set and normal-moveout spectrum signal-to-noise ratio require, using the inclination maximum as tilt threshold.
It is according to embodiment of the present invention, described that along the inclination direction to redefine obliquity structure corresponding square
Root speed is obtained by following steps: when the slope coefficient is greater than the tilt threshold, determining the geology of the obliquity structure
Construction trend and earthquake main profile angular separation；According to formula Gyration seismic velocity is calculated,
The new azimuth seismic velocity direction made and geological structure trend are consistent.
According to embodiment of the present invention, the method includes, by each single-point of obliquity structure along new azimuth
Direction is grouped, and is adjusted weight relationship by different orientations grouping range, is increased new azimuth direction data weight, determine speed
Degree spectrum and channel set.
According to embodiment of the present invention, the method includes, to each single-point centered on the single-point, along new side
Parallactic angle direction increases CRP gather number, and selection forms channel set and normal-moveout spectrum signal-to-noise ratio and the optimal road of resolution ratio
Number；Increase common reflection point number along vertical new azimuth direction, selection forms channel set and normal-moveout spectrum signal-to-noise ratio and resolution ratio
Optimal road number；Choose the channel set and normal-moveout spectrum of optimal signal-to-noise ratio and resolution ratio.
According to an aspect of the invention, there is provided adaptively moving towards normal-moveout spectrum seeks system, including seismic data is preliminary
Processing module 1 carries out denoising for obtaining seismic data；Velocity analysis module 2, for along earthquake main profile direction pair
The seismic data carries out velocity analysis, obtains the first root mean sequare velocity；Normal-moveout spectrum and seismic profile construct module 3, are used for edge
Normal-moveout spectrum and superposition seismic profile are formed according to first root mean sequare velocity；Migrated section constructs module 4, for establishing offset
The superposition seismic profile is deviated, obtains offset seismic profile by velocity field；Obliquity structure analysis module 5, for analyzing
The offset seismic profile obtains obliquity structure coefficient；The obliquity structure coefficient refers to incline corresponding to each single-point of obliquity structure
Angle；Turn azimuth determination module 6, for setting tilt threshold, when the inclination angle overall coefficient is less than the tilt threshold, according to
First root mean sequare velocity forms normal-moveout spectrum；When the obliquity structure coefficient is greater than the tilt threshold, along the inclination direction
The corresponding root mean sequare velocity of obliquity structure is redefined, the second root mean sequare velocity is obtained, according to the second root mean sequare velocity weight
It is new to determine normal-moveout spectrum.
According to embodiment of the present invention, the system further includes Tectonic Analysis module, for definitely
Shake processing pitch region geological structure trend；The amendment inclination angle knot is moved towards according to seismic processing pitch region geological structure
Structure coefficient.
According to embodiment of the present invention, the system further includes tilt threshold determining module, the inclination angle threshold
Value determining module is used to form seismic profile and dip angle attribute section；Two or more are chosen from the obliquity structure region of seismic profile
A test point；Single-point normal-moveout spectrum and corresponding trace gather are formed along earthquake main profile direction to each test point；To each test point edge
Earthquake main profile direction forms supervelocity spectrum and corresponding trace gather；Analyze the single-point normal-moveout spectrum and corresponding trace gather and supervelocity spectrum and
Corresponding trace gather, determines the inclination maximum for meeting channel set and normal-moveout spectrum signal-to-noise ratio requirement, using the inclination maximum as inclination angle threshold
Value.
According to embodiment of the present invention, the system include turn Azimuth, Speed, Altitude spectrum determining module, for pair
Each single-point increases CRP gather number centered on the single-point, along new azimuth direction, and selection forms channel set and speed
Degree spectrum signal-to-noise ratio and the optimal road number of resolution ratio；Increase common reflection point number along vertical new azimuth direction, chooses and form super road
Collection and normal-moveout spectrum signal-to-noise ratio and the optimal road number of resolution ratio；Choose the channel set and speed of optimal signal-to-noise ratio and resolution ratio
Spectrum.
The present invention, which passes through, determines tilt threshold, can conveniently judge whether to need to carry out seismic data to turn orientation
Angle processing；Suitable test point is chosen by complex geological structure region, the geologicstructure area for having obliquity structure is determined
Suitable tilt threshold, make to whether carry out turn it is azimuthal judgement it is more accurate；By each single-point of obliquity structure along new side
Parallactic angle direction is grouped, make each grouping turn azimuth processing it is more accurate；By adaptively turning azimuth, high-resolution is obtained
The channel set and supervelocity of rate high s/n ratio are composed.
Detailed description of the invention
Fig. 1 is adaptively to move towards normal-moveout spectrum acquiring method step schematic diagram；
Fig. 2 is the step schematic diagram of determining tilt threshold；
Fig. 2-a to 2-e is the schematic diagram to the acquirement process of the tilt threshold；And
Fig. 3 is the schematic diagram for adaptively moving towards normal-moveout spectrum and seeking system.
Specific embodiment
The preferred embodiments of the present invention will be described in detail with reference to the accompanying drawing, and reference label refers to the group in the present invention
Part, technology, realizing under appropriate circumstances so as to advantages and features of the invention can be easier to be understood.Following description is pair
The materialization of the claims in the present invention, and other specific implementations not clearly stated relevant to claim also belong to power
The range that benefit requires.
Fig. 1, which is shown, adaptively moves towards normal-moveout spectrum acquiring method step schematic diagram.
As shown in Figure 1, adaptively moving towards normal-moveout spectrum acquiring method, including first step S1, seismic data is obtained, is gone
It makes an uproar；Second step S2 carries out velocity analysis to the seismic data along earthquake main profile direction, obtains the first root mean sequare velocity；
Third step S3 forms normal-moveout spectrum and superposition seismic profile according to first root mean sequare velocity；Four steps S4 establishes offset
The superposition seismic profile or common midpoint gather are deviated, obtain offset seismic profile by velocity field；5th step S5,
It analyzes the offset seismic profile and obtains obliquity structure coefficient, the obliquity structure coefficient refers to corresponding to each single-point of obliquity structure
Inclination angle；6th step S6 sets tilt threshold, square according to first when the inclination angle overall coefficient is less than the tilt threshold
Root speed forms normal-moveout spectrum；When the obliquity structure coefficient is greater than the tilt threshold, redefined along the inclination direction
The corresponding root mean sequare velocity of obliquity structure, obtains the second root mean sequare velocity, redefines speed according to second root mean sequare velocity
Degree spectrum.
Earthquake main profile direction is usually consistent with earthquake-capturing direction, that is, field seismometer puts direction.
Original earthquake single-shot has usually contained surface wave interference, linear disturbance wave etc., needs to carry out prestack denoising, prestack
Adaptive surface wave attenuation, subraction compacting linear disturbance, incorrect energy decaying etc. can be used in denoising module.
Mean-square-root velocity field is the common reflection point time-distance equation for HORIZONTAL LAYERED MEDIUM WITH HIGH ACCURACY, can use hyp time curve
Company approximately replaces.In seism processing, the various speed using speed from 2000 to 7000m/s are repeatedly to concentrically
Point trace gather carries out TEC time error correction and superposition, the obtained stack result juxtaposition of each speed in speed-round trip zero shot-geophone distance
In time plane, referred to as normal-moveout spectrum.We are picked up acquisition mean-square-root velocity field in normal-moveout spectrum, carry out TEC time error correction processing
So that the seismic channel of common midpoint gather is aligned in prestack, and then it is superimposed and obtains superposition seismic profile.
According to first mean-square-root velocity field, a smoother velocity field is obtained, as migration velocity field, in altogether
Heart point trace gather carries out offset and obtains offset domain CRP gather, then carries out normal-moveout spectrum calculating based on this, picks up again inclined
Velocity field is moved, successive ignition, obtains final migration velocity field and deviated repeatedly, exports final offset seismic profile and is total to
Reflection point trace gather.
Signal-to-noise ratio and resolution when inclination angle reaches certain degree in survey area geological structure, to the normal-moveout spectrum of synthesis
Influence caused by rate can be very big, so that influencing the later period to work such as the explanation of seismic data and seismic inversions.
When the inclination angle is smaller, it is several with the signal-to-noise ratio of normal-moveout spectrum and resolution ratio that are formed along line direction to turn azimuth
There is no difference, at this point, if will increase the calculating of data using azimuth processing seismic data is turned, reducing efficiency.
The tilt threshold refers in certain specific survey area, influences the inclination angle size of normal-moveout spectrum signal-to-noise ratio and resolution ratio
Critical value.When inclination angle is greater than tilt threshold, the signal-to-noise ratio and resolution ratio of the normal-moveout spectrum influence the processing of later period seismic data
Larger, when inclination angle is less than tilt threshold, influence can be ignored.
When the obliquity structure coefficient is greater than the tilt threshold, obliquity structure pair is redefined along the inclination direction
The root mean sequare velocity answered obtains the second root mean sequare velocity, redefines normal-moveout spectrum according to second root mean sequare velocity.Described in
Inclination direction redefines the corresponding root mean sequare velocity of obliquity structure and refers to root mean square to the geologicstructure area with inclination angle
Speed is redefined.
According to embodiment of the present invention, the method includes determining seismic processing pitch region geological structure
Trend；The amendment obliquity structure coefficient is moved towards according to seismic processing pitch region geological structure.
Pitch region geological structure trend can be primarily determined by the seismic profile that earthquake information is formed, so
It is revised afterwards using well logging information.
Based on seismic profile, dip angle attribute section is extracted, determines and is greater than tilt threshold region, when using seismic profile
Between be sliced or dip angle attribute section isochronous surface determine be greater than tilt threshold regional tectonics angle of strike.
Fig. 2 shows the step schematic diagrames of determining tilt threshold.
As shown in Fig. 2, the tilt threshold is obtained by single-spot testing: step 1 A1, according to earthquake information acquisition earthquake
Section and dip angle attribute section；Step 2 A2 chooses two or more test points from the obliquity structure region of seismic profile；Step
Three A3 form single-point normal-moveout spectrum and corresponding trace gather along earthquake main profile direction to each test point；To each test point along earthquake
Main profile direction forms supervelocity spectrum and corresponding trace gather；Step 4 A4 analyzes the single-point normal-moveout spectrum and corresponding trace gather and hypervelocity
Degree spectrum and corresponding trace gather, determine the inclination maximum for meeting channel set and normal-moveout spectrum signal-to-noise ratio requirement, using the inclination maximum as
Tilt threshold.
Dip angle attribute section is a kind of independently of seismic amplitude, and unaffected Seismic Attribute Parameters, it can be with
Reflect the change of pitch angle of stratum or reflecting interface.
Now it is illustrated with acquirement process of Fig. 2-a to 2-e to the tilt threshold.
Apparent dip in figure is the inclination angle in specification.
As shown in Fig. 2-a, left side is according to earthquake information acquisition seismic profile in figure, and right side is dip angle attribute section.Ground
It shakes and chooses No. 1 test point and No. 2 test points on section, such as the position of two vertical lines in Fig. 2-a.Simultaneously on dip angle attribute section
Also the position of No. 1 test point and No. 2 test points is marked.It determines that pitch region geology is moved towards, marks dip angle attribute and cut open
The range of face top rade.It for No. 1 test point and No. 2 test points, can be obtained from Fig. 2-a: the inclination angle of No. 1 test point on the whole
Range: 0 degree to 4 degree of apparent dip, 5 degree to 9 degree of apparent dip；Inclination angle range of No. 2 test points in obliquity structure geologicstructure area:
Apparent dip is greater than 10 degree (part that right figure is drawn a circle).
Fig. 2-b is the information in-a according to fig. 2, and No. 1 test point is along line direction (i.e. grid direction) according to seismic data
Form channel set and normal-moveout spectrum.It can be concluded that, compared with the single-point normal-moveout spectrum in left side and corresponding trace gather, surpassed by Fig. 2-b
Normal-moveout spectrum and corresponding channel set signal-to-noise ratio are increased.
Fig. 2-c is according to the information in such as 2-a, and No. 2 test points are along line direction (i.e. toward all directions) according to seismic data
The channel set and normal-moveout spectrum of formation.By in Fig. 2-c it can be concluded that, compared with the single-point normal-moveout spectrum in left side and corresponding trace gather,
Supervelocity spectrum and corresponding channel set signal-to-noise ratio and clarity have received large effect.
By the test point that selection is a series of, the tilt threshold in some pitch region is determined, it can be to seismic data
It carries out processing and retrieves signal-to-noise ratio and the high channel set of clarity and supervelocity spectrum.The quantity of the test point can be according to surveying
The architectonic complexity in exploratory area domain is selected.
Fig. 2-d is adaptively to turn azimuthal schematic diagram.
By taking No. 2 test points as an example, be determined by the geological structure trend to No. 2 test points, by No. 2 test points along
Construction trend forms new channel set by speed point.Furthermore it is also possible to by partial velocity on the direction of vertical configuration trend
The data of point are modified to the data being upwardly formed are walked along construction.
Fig. 2-e is No. 2 test points after adaptively turning azimuth, along the super road that construction trend is formed by seismic data
Collection and normal-moveout spectrum.Compared with single-point normal-moveout spectrum and corresponding trace gather, signal-to-noise ratio and clarity are greatly improved.
A kind of embodiment according to the present invention, described along the inclination direction to redefine obliquity structure corresponding square
Root speed is obtained by following steps:
When the slope coefficient is greater than the tilt threshold, the geological structure trend and earthquake of the obliquity structure are determined
Main profile angular separation；
According to formula
Gyration seismic velocity is calculated, the new azimuth seismic velocity direction made and geological structure trend are consistent.
Each single-point of obliquity structure is grouped along new azimuth direction, presses by a kind of embodiment according to the present invention
Different orientations are grouped range and adjust weight relationship, increase new azimuth direction data weight, determine normal-moveout spectrum and channel set.
When determine carry out turning azimuth re-forming channel set when, can be true according to the difference of geological structure trend
The azimuth of fixed different single-points, can take mean value also can be set as needed different weights, determines and adaptively turns azimuth
Size, according to seismic data obtain supervelocity spectrum and channel set.
According to embodiment of the present invention, to each single-point centered on the single-point, increase altogether along new azimuth direction
Reflection point trace gather number, selection form channel set and normal-moveout spectrum signal-to-noise ratio and the optimal road number of resolution ratio；Along vertical new orientation
Angular direction increases common reflection point number, and selection forms channel set and normal-moveout spectrum signal-to-noise ratio and the optimal road number of resolution ratio；It chooses
The channel set and normal-moveout spectrum of optimal signal-to-noise ratio and resolution ratio.
It can successively be chosen properly in azimuth direction and vertical orientations angular direction durings forming new channel set
Trace gather, form new channel set.
Fig. 3, which is shown, adaptively moves towards the schematic diagram that normal-moveout spectrum seeks system.
As shown in figure 3, adaptively turning Azimuth, Speed, Altitude spectrum seeks system, including, seismic data preliminary treatment module 1 is used
In obtaining seismic data, denoised；Velocity analysis module 2, for being carried out along earthquake main profile direction to the seismic data
Velocity analysis obtains the first root mean sequare velocity；Normal-moveout spectrum and seismic profile construct module 3, for along square according to described first
Root speed forms normal-moveout spectrum and superposition seismic profile；Migrated section constructs module 4, will be described folded for establishing migration velocity field
Add seismic profile to be deviated, obtains offset seismic profile；Obliquity structure analysis module 5 is cutd open for analyzing the offset earthquake
Face obtains obliquity structure coefficient, and the obliquity structure coefficient refers to inclination angle corresponding to each single-point of obliquity structure；Turn azimuth judgement
Module 6, for setting tilt threshold, when the inclination angle overall coefficient is less than the tilt threshold, according to the first root mean sequare velocity
Form normal-moveout spectrum；When the obliquity structure coefficient is greater than the tilt threshold, inclination angle is redefined along the inclination direction and is tied
The corresponding root mean sequare velocity of structure, obtains the second root mean sequare velocity, redefines normal-moveout spectrum according to second root mean sequare velocity.
When the obliquity structure coefficient is less than the tilt threshold, no longer carry out turning azimuth operation, when the inclination angle
When structural coefficient is greater than the tilt threshold, need to by the processed seismic data weight of seismic data preliminary treatment module 1
It is new to carry out velocity analysis, different road formation channel sets is extracted to needing to turn azimuthal each speed single-point.
A kind of embodiment according to the present invention, the system further include Tectonic Analysis module, for definitely
Shake processing pitch region geological structure trend；The amendment inclination angle knot is moved towards according to seismic processing pitch region geological structure
Structure coefficient.
A kind of embodiment according to the present invention, the system further include tilt threshold determining module, the inclination angle threshold
Value determining module is used to form seismic profile and dip angle attribute section；Two or more are chosen from the obliquity structure region of seismic profile
A test point；Single-point normal-moveout spectrum and corresponding trace gather are formed along earthquake main profile direction to each test point；To each test point edge
Earthquake main profile direction forms supervelocity spectrum and corresponding trace gather；Analyze the single-point normal-moveout spectrum and corresponding trace gather and supervelocity spectrum and
Corresponding trace gather, determines the inclination maximum for meeting channel set and normal-moveout spectrum signal-to-noise ratio requirement, using the inclination maximum as inclination angle threshold
Value.
The tilt threshold determining module can choose two or more test points, and the position of selection and quantity can be with
It is determined according to architectonic complexity.
A kind of embodiment according to the present invention, the system, including turn Azimuth, Speed, Altitude spectrum determining module, for pair
Each single-point increases CRP gather number centered on the single-point, along new azimuth direction, and selection forms channel set and speed
Degree spectrum signal-to-noise ratio and the optimal road number of resolution ratio；Increase common reflection point number along vertical new azimuth direction, chooses and form super road
Collection and normal-moveout spectrum signal-to-noise ratio and the optimal road number of resolution ratio；Choose the channel set and speed of optimal signal-to-noise ratio and resolution ratio
Spectrum.
The present invention, which passes through, determines tilt threshold, can conveniently judge whether to need to carry out seismic data to turn orientation
Angle processing；Suitable test point is chosen by complex geological structure region, the geologicstructure area for having obliquity structure is determined
Suitable tilt threshold, make to whether carry out turn it is azimuthal judgement it is more accurate；By each single-point of obliquity structure along new side
Parallactic angle direction is grouped, make each grouping turn azimuth processing it is more accurate；By adaptively turning azimuth, high-resolution is obtained
The channel set and supervelocity of rate high s/n ratio are composed.
It should be noted that the above-mentioned embodiments illustrate rather than limit the invention, and this
Field technical staff can be designed alternative embodiment without departing from the scope of the appended claims.In claim
In, any reference symbol between parentheses should not be configured to limitations on claims.
Claims (10)
1. adaptively moving towards normal-moveout spectrum acquiring method, wherein including
First step (S1) obtains seismic data, carries out denoising；
Second step (S2) carries out velocity analysis to the seismic data along earthquake main profile direction, obtains the first root mean square speed
Degree；
Third step (S3) forms normal-moveout spectrum and superposition seismic profile according to first root mean sequare velocity；
Four steps (S4), establishes migration velocity field, and the superposition seismic profile or common midpoint gather are deviated, obtained
Offset seismic profile；
5th step (S5) analyzes the offset seismic profile and obtains obliquity structure coefficient, and the obliquity structure coefficient refers to inclination angle
Inclination angle corresponding to each single-point of structure；
6th step (S6) sets tilt threshold, square according to first when the inclination angle overall coefficient is less than the tilt threshold
Root speed forms normal-moveout spectrum；
When the obliquity structure coefficient is greater than the tilt threshold, it is corresponding that obliquity structure is redefined along the inclination direction
Root mean sequare velocity obtains the second root mean sequare velocity, redefines normal-moveout spectrum according to second root mean sequare velocity.
2. according to the method described in claim 1, including determining that seismic processing pitch region geological structure moves towards；
The amendment obliquity structure coefficient is moved towards according to seismic processing pitch region geological structure.
3. according to the method described in claim 1, the tilt threshold is obtained by single-spot testing:
According to earthquake information acquisition seismic profile and dip angle attribute section；
Two or more test points are chosen from the obliquity structure region of seismic profile；
Single-point normal-moveout spectrum and corresponding trace gather are formed along earthquake main profile direction to each test point；
Supervelocity spectrum and corresponding trace gather are formed along earthquake main profile direction to each test point；
The single-point normal-moveout spectrum and corresponding trace gather and supervelocity spectrum and corresponding trace gather are analyzed, determines and meets channel set and normal-moveout spectrum letter
It makes an uproar than desired inclination maximum, using the inclination maximum as tilt threshold.
4. according to the method described in claim 1, described redefine the corresponding root mean square of obliquity structure along the inclination direction
Speed is obtained by following steps:
When the slope coefficient is greater than the tilt threshold, the geological structure trend and the main survey of earthquake of the obliquity structure are determined
Line angular separation；
According to formulaCalculate gyration seismic velocity, the new azimuth seismic velocity direction made with
Geological structure trend is consistent.
5. according to the method described in claim 4, include being grouped each single-point of obliquity structure along new azimuth direction,
Weight relationship is adjusted by different orientations grouping range, increases new azimuth direction data weight, determines normal-moveout spectrum and super road
Collection.
6. according to the method described in claim 5, including to each single-point centered on the single-point, increasing along new azimuth direction
CRP gather number, selection form channel set and normal-moveout spectrum signal-to-noise ratio and the optimal road number of resolution ratio；
Increase common reflection point number along vertical new azimuth direction, selection forms channel set and normal-moveout spectrum signal-to-noise ratio and resolution ratio
Optimal road number；
Choose the channel set and normal-moveout spectrum of optimal signal-to-noise ratio and resolution ratio.
7. adaptively moving towards normal-moveout spectrum seeks system, including
Seismic data preliminary treatment module (1) is denoised for obtaining seismic data；
Velocity analysis module (2) obtains first for carrying out velocity analysis to the seismic data along earthquake main profile direction
Root speed；
Normal-moveout spectrum and seismic profile building module (3), for forming normal-moveout spectrum and superposition along according to first root mean sequare velocity
Seismic profile；
Migrated section constructs module (4), for establishing migration velocity field, by the superposition seismic profile or common midpoint gather into
Line displacement obtains offset seismic profile；
Obliquity structure analysis module (5) obtains obliquity structure coefficient for analyzing the offset seismic profile；The obliquity structure
Coefficient refers to inclination angle corresponding to each single-point of obliquity structure；
Turn azimuth determination module (6), for setting tilt threshold, when the inclination angle overall coefficient is less than the tilt threshold,
Normal-moveout spectrum is formed according to the first root mean sequare velocity；When the obliquity structure coefficient is greater than the tilt threshold, along the inclination angle
Direction redefines the corresponding root mean sequare velocity of obliquity structure, obtains the second root mean sequare velocity, according to the second root mean square speed
Degree redefines normal-moveout spectrum.
8. system according to claim 7 further includes Tectonic Analysis module, for determining seismic processing pitch region
Geological structure trend；
The amendment obliquity structure coefficient is moved towards according to seismic processing pitch region geological structure.
9. system according to claim 7 further includes tilt threshold determining module, the tilt threshold determining module is used for
Form seismic profile and dip angle attribute section；
Two or more test points are chosen from the obliquity structure region of seismic profile；
Single-point normal-moveout spectrum and corresponding trace gather are formed along earthquake main profile direction to each test point；
Supervelocity spectrum and corresponding trace gather are formed along earthquake main profile direction to each test point；
The single-point normal-moveout spectrum and corresponding trace gather and supervelocity spectrum and corresponding trace gather are analyzed, determines and meets channel set and normal-moveout spectrum letter
It makes an uproar than desired inclination maximum, using the inclination maximum as tilt threshold.
10. system according to claim 7, including turn Azimuth, Speed, Altitude spectrum determining module, it is used for each single-point with the list
Centered on point, increase CRP gather number along new azimuth direction, selection formed channel set and normal-moveout spectrum signal-to-noise ratio and
The optimal road number of resolution ratio；
Increase common reflection point number along vertical new azimuth direction, selection forms channel set and normal-moveout spectrum signal-to-noise ratio and resolution ratio
Optimal road number；
Choose the channel set and normal-moveout spectrum of optimal signal-to-noise ratio and resolution ratio.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910049867.9A CN109884693B (en) | 2019-01-18 | 2019-01-18 | Self-adaptive trend velocity spectrum solving method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910049867.9A CN109884693B (en) | 2019-01-18 | 2019-01-18 | Self-adaptive trend velocity spectrum solving method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109884693A true CN109884693A (en) | 2019-06-14 |
| CN109884693B CN109884693B (en) | 2021-10-15 |
Family
ID=66926293
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910049867.9A Active CN109884693B (en) | 2019-01-18 | 2019-01-18 | Self-adaptive trend velocity spectrum solving method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109884693B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111239825A (en) * | 2020-03-09 | 2020-06-05 | 辽宁工程技术大学 | Method for determining seismic wave velocity by ellipse method |
| CN113514880A (en) * | 2020-04-10 | 2021-10-19 | 中国石油化工股份有限公司 | Construction method for superposition speed at intersection point of two-dimensional measuring lines |
| CN115840874A (en) * | 2023-02-21 | 2023-03-24 | 中海石油（中国）有限公司北京研究中心 | Method and device for determining source control forward product inclination angle, computer storage medium and equipment |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5097452A (en) * | 1991-05-21 | 1992-03-17 | Western Atlas International, Inc. | Analysis of migration velocity by migration of velocity spectra |
| CN102221709A (en) * | 2011-06-01 | 2011-10-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Velocity analysis and dynamic correction method based on formation parameter information |
| CN102879814A (en) * | 2011-07-15 | 2013-01-16 | 中国石油天然气集团公司 | Accurate depth domain layer speed updating method |
| CN103064113A (en) * | 2011-10-19 | 2013-04-24 | 恒泰艾普石油天然气技术服务股份有限公司 | Two-parameter expansion method and device confirming underground velocity structure |
| CN103217712A (en) * | 2012-01-19 | 2013-07-24 | 中国石油化工股份有限公司 | Earthquake velocity analysis method along geologic structure trend |
| WO2014191650A1 (en) * | 2013-05-30 | 2014-12-04 | Total Sa | Methods and devices for acquiring seismic data relative to an underground area beneath the sea |
| CN104614762A (en) * | 2014-12-17 | 2015-05-13 | 中国石油天然气股份有限公司 | Loose sandstone gas reservoir boundary determination method and device |
| CN105676278A (en) * | 2016-01-15 | 2016-06-15 | 中国海洋石油总公司 | Method and device for processing seismic wave data |
| CN107179546A (en) * | 2017-05-11 | 2017-09-19 | 中国石油天然气股份有限公司 | A kind of method and device for determining velocity field |
| CN107884827A (en) * | 2016-09-29 | 2018-04-06 | 中国石油化工股份有限公司 | Earthquake overlap speed turns the method and device of interval velocity |
| CN107957594A (en) * | 2017-11-15 | 2018-04-24 | 中国石油集团东方地球物理勘探有限责任公司 | The oval bearing calibration of seismic data, dynamic bearing calibration and normal-moveout spectrum computational methods |
-
2019
- 2019-01-18 CN CN201910049867.9A patent/CN109884693B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5097452A (en) * | 1991-05-21 | 1992-03-17 | Western Atlas International, Inc. | Analysis of migration velocity by migration of velocity spectra |
| CN102221709A (en) * | 2011-06-01 | 2011-10-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Velocity analysis and dynamic correction method based on formation parameter information |
| CN102879814A (en) * | 2011-07-15 | 2013-01-16 | 中国石油天然气集团公司 | Accurate depth domain layer speed updating method |
| CN103064113A (en) * | 2011-10-19 | 2013-04-24 | 恒泰艾普石油天然气技术服务股份有限公司 | Two-parameter expansion method and device confirming underground velocity structure |
| CN103217712A (en) * | 2012-01-19 | 2013-07-24 | 中国石油化工股份有限公司 | Earthquake velocity analysis method along geologic structure trend |
| WO2014191650A1 (en) * | 2013-05-30 | 2014-12-04 | Total Sa | Methods and devices for acquiring seismic data relative to an underground area beneath the sea |
| CN104614762A (en) * | 2014-12-17 | 2015-05-13 | 中国石油天然气股份有限公司 | Loose sandstone gas reservoir boundary determination method and device |
| CN105676278A (en) * | 2016-01-15 | 2016-06-15 | 中国海洋石油总公司 | Method and device for processing seismic wave data |
| CN107884827A (en) * | 2016-09-29 | 2018-04-06 | 中国石油化工股份有限公司 | Earthquake overlap speed turns the method and device of interval velocity |
| CN107179546A (en) * | 2017-05-11 | 2017-09-19 | 中国石油天然气股份有限公司 | A kind of method and device for determining velocity field |
| CN107957594A (en) * | 2017-11-15 | 2018-04-24 | 中国石油集团东方地球物理勘探有限责任公司 | The oval bearing calibration of seismic data, dynamic bearing calibration and normal-moveout spectrum computational methods |
Non-Patent Citations (4)
| Title |
|---|
| EVGENY LANDA 等: "separation，imaging，and velocity analysis of seismic diffractions using migrated dip-angle gathers", 《SEG TECHNICAL PROGRAM EXPANDED ABSTRACTS》 * |
| 徐常练: "三参量速度分析中的构造形态约束", 《石油地球物理勘探》 * |
| 缪志伟 等: "三江盆地浓江凹陷烃源岩热成熟度定量预测", 《物探化探计算技术》 * |
| 项龙云: "基于体模式的速度分析方法研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111239825A (en) * | 2020-03-09 | 2020-06-05 | 辽宁工程技术大学 | Method for determining seismic wave velocity by ellipse method |
| CN113514880A (en) * | 2020-04-10 | 2021-10-19 | 中国石油化工股份有限公司 | Construction method for superposition speed at intersection point of two-dimensional measuring lines |
| CN113514880B (en) * | 2020-04-10 | 2024-02-20 | 中国石油化工股份有限公司 | Method for constructing superposition speed at intersection point of two-dimensional measuring lines |
| CN115840874A (en) * | 2023-02-21 | 2023-03-24 | 中海石油（中国）有限公司北京研究中心 | Method and device for determining source control forward product inclination angle, computer storage medium and equipment |
| CN115840874B (en) * | 2023-02-21 | 2023-04-28 | 中海石油（中国）有限公司北京研究中心 | Source control forward product dip angle determining method, device, computer storage medium and equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109884693B (en) | 2021-10-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109738945B (en) | Method for directly generating construction diagram by using prestack depth migration result | |
| CN104678434B (en) | Method for predicting storage layer crack development parameters | |
| CN102841379B (en) | Method for analyzing pre-stack time migration and speed based on common scatter point channel set | |
| CN103733089B (en) | For including the system and method that the underground of uncertainty estimation characterizes | |
| CN109884693A (en) | Adaptively move towards normal-moveout spectrum acquiring method and system | |
| CN104280777B (en) | Method for suppressing interference of seismic data multiples on land | |
| CN104297784A (en) | Primary wave azimuthal anisotropy based fracture predicting method | |
| CN105093301B (en) | The generation method and device of common imaging point angle of reflection angle gathers | |
| CN104570076A (en) | Automatic seismic wave first-arrival picking method based on dichotomy | |
| CN109839660A (en) | A method of velocity depth model is established using prestack trace gather data | |
| CN107728204A (en) | Based on the anisotropic crack prediction method of prestack compressional wave and system | |
| CN109164487A (en) | A kind of method based on model foundation average velocity field and finely at the method for figure | |
| CN111983683B (en) | Prediction method and system for lake-facies limestone reservoir under low-well condition | |
| CN106443777B (en) | The acquisition methods and converted-wave velocity analysis method of converted wave anisotropic parameters spectrum | |
| CN106324682B (en) | A kind of surface structure survey method applied to frozen soil layer area | |
| CN106199704B (en) | A kind of Three-dimendimal fusion submarine cable seismic data velocity modeling method | |
| CN107656308B (en) | A kind of common scattering point pre-stack time migration imaging method based on time depth scanning | |
| CN103076628B (en) | The disposal route of the pre-stack time migration that a kind of aperture is optimized | |
| CN105093318B (en) | A kind of adaptive wave equation wave field extrapolation static correcting method | |
| CN106443791B (en) | The method for asking for tilted stratum or anisotropic formation shear wave Value of residual static correction | |
| CN106990434B (en) | Ellipse expansion converted wave imaging method and system | |
| CN106226815B (en) | The acquisition methods and converted-wave velocity analysis method of converted wave velocity equivalent comparison | |
| CN106125136B (en) | A kind of common point main road set creation method of relief surface | |
| CN109143398B (en) | Modeling method for automatic grid chromatography depth domain speed | |
| CN111399031A (en) | Method for acquiring and processing mountain land seismic data |
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 |