CN107329169A - The extracting method and device of a kind of angular-trace gather - Google Patents
The extracting method and device of a kind of angular-trace gather Download PDFInfo
- Publication number
- CN107329169A CN107329169A CN201710628826.6A CN201710628826A CN107329169A CN 107329169 A CN107329169 A CN 107329169A CN 201710628826 A CN201710628826 A CN 201710628826A CN 107329169 A CN107329169 A CN 107329169A
- Authority
- CN
- China
- Prior art keywords
- mrow
- incident angle
- ray
- msup
- relation
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
Abstract
The application embodiment discloses the extracting method and device of a kind of angular-trace gather, wherein, methods described includes：Geological data and speed data are obtained, and determines based on the geological data and the speed data incident angle of ray；Based on the speed data and default ray parameter, the Fresnel-zone width of the ray is determined；Enter row interpolation weighting to the geological data in the Fresnel-zone width range, and the result that interpolation is weighted is used as the corresponding angular-trace gather of the incident angle.The extracting method and device for a kind of angular-trace gather that the application is provided, it is possible to increase the extraction quality of angular-trace gather, so as to improve the imaging effect of angular stack section and improve AVA analysis precision.
Description
Technical field
The application is related to seismic data processing technology field, the extracting method and device of more particularly to a kind of angular-trace gather.
Background technology
AVO (Amplitude variation with offset, amplitude with offset distance change) technology is to utilize amplitude
With offset distance Variations and the seismic exploration technique for distinguishing lithology and oil-gas reservoir.AVO analyses are carried out, key is prestack
The fine processing of geological data and recovery and holding relative amplitude information.Conventional AVO invertings are using anti-on angular-trace gather
Penetrate coefficient approximate formula and obtain AVO attributes, the amplitude of the technical requirements geological data is capable of inverting underground medium reflectance factor
Strong or weak relation.It has studied Different Strata inclination angle lower angle in the prior art with offset distance to change, it is indicated that AVA (Amplitude
Variation with Angle, amplitude with incident angle change) analysis must take into consideration stratigraphic dip influences, when there is tomography
When the back wave of the section ripple, diffracted wave and the target zone that are produced Deng stratum cross directional variations is interfered, on this basis AVA points
Analysis accuracy is poor, therefore can not be useful in the complicated area of strata condition, and angle calculation is strictly limited to rate pattern, accurate
True degree is relatively low.In order to overcome the defect of above-mentioned angular-trace gather, prior art is proposed extracts angular-trace gather based on migration before stack
Algorithm.Prestack migration method is both the important means for solving engineering construction system and velocity analysis, is again complicated structure region
AVA analyses, inverting provide strong support.
Migration before stack can be imaged while enable reflection line-ups to correct and playback to complex geologic body, utilize prestack
The angular-trace gather that the CMP or CRP trace gathers of skew are extracted, which carries out AVA invertings, can obtain higher precision.At present, a kind of method is carried
Having gone out the Amplitude preserving algorithm based on round trip wave equation enables migration imaging amplitude to reflect reflectance factor power change.In addition,
There is a method in which introducing the concept of the common picture point trace gather of underground offset distance, and angular-trace gather is converted to using slant stack, but
It is more suitable for carrying out migration velocity analysis.There is a method in which proposing angle road is extracted using double square root equation Harmonic Decomposition method
Collection, improves the extraction accuracy of angular-trace gather.
The characteristics of propagation of the seismic wave in complex dielectrics has multipath, can not with the common imaging gather of migration in offset domain
Correct to determine reflector position, the approach of solution is the common imaging gather replacement with angle domain.For seismic attribute abstraction, angle
The change (AVA) of amplitude can more reflect the change of lithology than amplitude with the change (AVO) of geophone offset in degree trace gather.As oil gas is surveyed
That visits continues to develop, and many scholars are being directed to the research to fluid Forecasting recognition.Poststack data mainly has time frequency analysis side
Method and the wave impedance inversion fluid identification based on seismic data.Hydrocarbon potential detection is carried out using AVO technologies, forefathers have done
Substantial amounts of research work, and achieve good effect in the prediction of different types of oil and gas reservoir, it is seen that angular-trace gather oil
Important function in gas exploration.
Angular-trace gather acquiring method mainly has two classes, and a class is wave equation theoretical method, and another kind of is ray theory side
Method.Angular-trace gather extracting method efficiency high based on ray theory, is the main flow of current industrial application.But, it is currently based on and penetrates
The angular-trace gather extracting method of lineation opinion, the quality of the angular-trace gather extracted is generally relatively low, so as to cause angular stack section
Imaging effect it is not good.
The content of the invention
The purpose of the application embodiment is to provide the extracting method and device of a kind of angular-trace gather, it is possible to increase angle road
The extraction quality of collection, so as to improve the imaging effect of angular stack section.
To achieve the above object, the application embodiment provides a kind of extracting method of angular-trace gather, and methods described includes：
Geological data and speed data are obtained, and determines based on the geological data and the speed data incident angle of ray；
Based on the speed data and default ray parameter, the Fresnel-zone width of the ray is determined；To the Fresnel bandwidth
Geological data in the range of degree enters row interpolation weighting, and the result that interpolation is weighted is used as the corresponding angle road of the incident angle
Collection.
Further, determine that the incident angle of ray includes based on the geological data and the speed data：
Set up default the first relation between ray parameter and incident angle；Wherein, layer is also included in first relation
Speed；
The second relation when setting up the default ray parameter and offset distance, travelling, between layer root mean sequare velocity；
Based on first relation and second relation, the incident angle of ray is calculated.
Further, first relation is determined according to the following equation：
Wherein, P represents the default ray parameter, and θ represents the incident angle, VintRepresent the interval velocity.
Further, second relation is determined according to the following equation：
Wherein, P represents the default ray parameter, and X represents the offset distance, when t represents the travelling, VrmsRepresent institute
State a layer root mean sequare velocity.
Further, the incident angle of ray is calculated according to the following equation：
Wherein, θ represents the incident angle, and X represents the offset distance, VintRepresent the interval velocity, VrmsRepresent described
Layer root mean sequare velocity, t0When representing zero shot-geophone distance bilayer travelling.
Further, the Fresnel-zone width of the ray is determined according to the following equation：
Wherein, W represents the Fresnel-zone width, VrmsThe layer root mean sequare velocity is represented, P represents the default ray
Parameter, t0When representing zero shot-geophone distance bilayer travelling, f0Represent the basic frequency of the geological data.
Further, row interpolation weighting is entered to the geological data in the Fresnel-zone width range to be included：
For all recorded traces in the range of outgoing in t0The amplitude at moment, by between the recorded trace and eye point away from
From the corresponding weight factor of amplitude as the recorded trace reciprocal, and based on the amplitude and corresponding weight factor
It is weighted superposition.
To achieve the above object, the application also provides a kind of extraction element of angular-trace gather, and described device includes：
Incident angle determining unit, for obtaining geological data and speed data, and based on the geological data and institute
State the incident angle that speed data determines ray；
Fresnel-zone width determining unit, for based on the speed data and default ray parameter, it is determined that described penetrate
The Fresnel-zone width of line；
Angular-trace gather determining unit, for entering row interpolation weighting to the geological data in the Fresnel-zone width range,
And the result for weighting interpolation is used as the corresponding angular-trace gather of the incident angle.
Further, the incident angle determining unit includes：
First relation sets up module, for setting up the first relation between default ray parameter and incident angle；Wherein, institute
Stating also includes interval velocity in the first relation；
Second relation sets up module, during for setting up the default ray parameter and offset distance, travelling, layer root mean sequare velocity
Between the second relation；
Incident angle computing module, for based on first relation and second relation, calculating the incidence angle of ray
Degree.
Further, the incident angle computing module calculates the incident angle of ray according to the following equation：
Wherein, θ represents the incident angle, and X represents the offset distance, VintRepresent the interval velocity, VrmsRepresent described
Layer root mean sequare velocity, t0When representing zero shot-geophone distance bilayer travelling.
Therefore, for any incident angle, inserted using the geological data in Fresnel-zone width range
Value weighting, is used as the angular-trace gather of the incident angle.Fresnel-zone width is the limit for differentiating geologic body interval, limited angle model
Partial stack in enclosing is corresponding with the width of Fresnel-zone, it is ensured that the resolution ratio of angular-trace gather；Meanwhile, utilize the luxuriant and rich with fragrance alunite of ray
Your band is theoretical, can accurately calculate the scope of ray incident angle, improves the quality of angular-trace gather.Model test and actual number
Show according to the study, the angular-trace gather obtained based on the theoretical angular-trace gather conversion method of Fresnel-zone, its AVA curve conversion trend
Accurately, it is bigger compared to the angle codomain obtained in the prior art in the angular range of deep layer, angular stack section imaging effect
More preferably.
Brief description of the drawings
, below will be to embodiment in order to illustrate more clearly of the application embodiment or technical scheme of the prior art
Or the accompanying drawing used required in description of the prior art is briefly described, it should be apparent that, drawings in the following description are only
It is some embodiments described in the application, for those of ordinary skill in the art, is not paying creative labor
On the premise of, other accompanying drawings can also be obtained according to these accompanying drawings.
A kind of extracting method flow chart for angular-trace gather that Fig. 1 provides for the application embodiment；；
Fig. 2 (a) is the geological data schematic diagram after denoising in the application；
Fig. 2 (b) is the angular-trace gather schematic diagram that extracts in the application；
Fig. 3 (a) is the stacked section schematic diagram of CRP trace gathers in the application；
Fig. 3 (b) is the angular-trace gather stacked section schematic diagram that extracts in the application；
Fig. 4 is the functional block diagram of the extraction element of angular-trace gather in the application.
Embodiment
In order that those skilled in the art more fully understand the technical scheme in the application, it is real below in conjunction with the application
The accompanying drawing in mode is applied, the technical scheme in the application embodiment is clearly and completely described, it is clear that described
Embodiment is only a part of embodiment of the application, rather than whole embodiments.Based on the embodiment party in the application
Formula, the every other embodiment that those of ordinary skill in the art are obtained under the premise of creative work is not made all should
When the scope for belonging to the application protection.
Referring to Fig. 1, the application provides a kind of extracting method of angular-trace gather, methods described includes：
S1：Geological data and speed data are obtained, and ray is determined based on the geological data and the speed data
Incident angle.
In the present embodiment, the geological data can be CMP (Common Middle Point, CMP) road
Collect data or CRP (Common Reflection Point, common reflection point) trace gather data.Specifically, the CMP trace gathers number
According to or the CRP trace gathers data can utilize Prestack Migration Technology so that amplitude playback.In addition, the CMP trace gathers data or institute
State the data that CRP trace gather data can also be amplitude distortion caused by corrected inhomogeneous formation.In the geological data,
The offset distance of each seismic channel can be included.
In the present embodiment, the speed data can be the speed mould obtained after being modeled in advance to geological data
Data in type.Specifically, the rate pattern on stratum can be obtained by geological data by technologies such as tomographies, so as to
To obtain the speed data.In the speed data, each layer of interval velocity can be included in stratum and to a certain layer
Layer root mean sequare velocity.
In the present embodiment, the incident angle of ray is determined based on the geological data and the speed data, can be with
Comprise the following steps：
S11：Set up default the first relation between ray parameter and incident angle；Wherein, also wrapped in first relation
Include interval velocity；
S12：The second relation when setting up the default ray parameter and offset distance, travelling, between layer root mean sequare velocity；
S13：Based on first relation and second relation, the incident angle of ray is calculated.
Specifically, in the present embodiment, it is possible to use Walden equations, the incident angle of ray is calculated.Specifically,
Assuming that in the case of multilayer dielectricity, first relation can be determined according to the following equation：
Wherein, P represents the default ray parameter, and θ represents the incident angle, VintRepresent the interval velocity.
In the present embodiment, examined when can set up travelling according to the following equation with offset distance, layer root mean sequare velocity, zero big gun
Contact between during away from bilayer travelling：
In addition, default ray parameter can be determined according to following differential formulas：
It therefore, it can determine second relation according to the following equation：
Wherein, P represents the default ray parameter, and X represents the offset distance, when t represents the travelling, VrmsRepresent institute
State a layer root mean sequare velocity.
So, first relation and the second relation are subjected to simultaneous, so as to calculate ray according to the following equation
Incident angle：
Wherein, θ represents the incident angle, and X represents the offset distance, VintRepresent the interval velocity, VrmsRepresent described
Layer root mean sequare velocity, t0When representing zero shot-geophone distance bilayer travelling.
S2：Based on the speed data and default ray parameter, the Fresnel-zone width of the ray is determined.
In the present embodiment, because seismic reflection is a kind of wave phenomenon, the energy of back wave is from a table
Face rather than a point, the width on this surface are generally determined by Fresnel-zone.In order to realize the portion based on Fresnel-zone
Divide superposition, it usually needs provide the scope of angular stack.Specifically, each layer can be included in stratum in the speed data
Layer root mean sequare velocity, then the Fresnel-zone width of the ray can be determined according to the following equation：
Wherein, W represents the Fresnel-zone width, VrmsThe layer root mean sequare velocity is represented, P represents the default ray
Parameter, t0When representing zero shot-geophone distance bilayer travelling, f0Represent the basic frequency of the geological data.
S3：Enter row interpolation weighting, and the result that interpolation is weighted to the geological data in the Fresnel-zone width range
It is used as the corresponding angular-trace gather of the incident angle.
Specifically, when entering row interpolation weighting, for all recorded traces in the range of outgoing in t0The amplitude at moment, by institute
Recorded trace weight factor corresponding with the amplitude reciprocal as the recorded trace of the distance between eye point is stated, and based on institute
State amplitude and corresponding weight factor is weighted superposition.
In actual applications, the geological data after denoising can be as shown in Fig. 2 (a), and the angular-trace gather of extraction can be such as Fig. 2
(b) shown in.In addition, the stacked section of CRP trace gathers can be as shown in Fig. 3 (a), the angular-trace gather stacked section of extraction can be such as figure
Shown in 3 (b).
Referring to Fig. 4, the application also provides a kind of extraction element of angular-trace gather, described device includes：
Incident angle determining unit 100, for obtaining geological data and speed data, and based on the geological data and
The speed data determines the incident angle of ray；
Fresnel-zone width determining unit 200, for based on the speed data and default ray parameter, it is determined that described
The Fresnel-zone width of ray；
Angular-trace gather determining unit 300, adds for entering row interpolation to the geological data in the Fresnel-zone width range
Power, and the result that interpolation is weighted is used as the corresponding angular-trace gather of the incident angle.
In the present embodiment, the incident angle determining unit 100 includes：
First relation sets up module, for setting up the first relation between default ray parameter and incident angle；Wherein, institute
Stating also includes interval velocity in the first relation；
Second relation sets up module, during for setting up the default ray parameter and offset distance, travelling, layer root mean sequare velocity
Between the second relation；
Incident angle computing module, for based on first relation and second relation, calculating the incidence angle of ray
Degree.
In the present embodiment, first relation, which sets up module, can determine first relation according to the following equation：
Wherein, P represents the default ray parameter, and θ represents the incident angle, VintRepresent the interval velocity.
In the present embodiment, second relation, which sets up module, can determine second relation according to the following equation：
Wherein, P represents the default ray parameter, and X represents the offset distance, when t represents the travelling, VrmsRepresent institute
State a layer root mean sequare velocity.
In the present embodiment, the incident angle computing module calculates the incident angle of ray according to the following equation：
Wherein, θ represents the incident angle, and X represents the offset distance, VintRepresent the interval velocity, VrmsRepresent described
Layer root mean sequare velocity, t0When representing zero shot-geophone distance bilayer travelling.
In the present embodiment, the Fresnel-zone width determining unit 200 can determine described penetrate according to the following equation
The Fresnel-zone width of line：
Wherein, W represents the Fresnel-zone width, VrmsThe layer root mean sequare velocity is represented, P represents the default ray
Parameter, t0When representing zero shot-geophone distance bilayer travelling, f0Represent the basic frequency of the geological data.
Therefore, for any incident angle, inserted using the geological data in Fresnel-zone width range
Value weighting, is used as the angular-trace gather of the incident angle.Fresnel-zone width is the limit for differentiating geologic body interval, limited angle model
Partial stack in enclosing is corresponding with the width of Fresnel-zone, it is ensured that the resolution ratio of angular-trace gather；Meanwhile, utilize the luxuriant and rich with fragrance alunite of ray
Your band is theoretical, can accurately calculate the scope of ray incident angle, improves the quality of angular-trace gather.Model test and actual number
Show according to the study, the angular-trace gather obtained based on the theoretical angular-trace gather conversion method of Fresnel-zone, its AVA curve conversion trend
Accurately, it is bigger compared to the angle codomain obtained in the prior art in the angular range of deep layer, angular stack section imaging effect
More preferably.
In the 1990s, for a technology improvement can clearly distinguish be on hardware improvement (for example,
Improvement to circuit structures such as diode, transistor, switches) or software on improvement (for the improvement of method flow).So
And, with the development of technology, the improvement of current many method flows can be considered as directly improving for hardware circuit.
Designer nearly all obtains corresponding hardware circuit by the way that improved method flow is programmed into hardware circuit.Cause
This, it cannot be said that the improvement of a method flow cannot be realized with hardware entities module.For example, PLD
(Programmable Logic Device, PLD) (such as field programmable gate array (Field Programmable Gate
Array, FPGA)) it is exactly such a integrated circuit, its logic function is determined by user to device programming.By designer
Voluntarily programming comes a digital display circuit " integrated " on a piece of PLD, without asking chip maker to design and make
Special IC chip.Moreover, nowadays, substitution manually makes IC chip, and this programming is also used instead mostly " patrols
Volume compiler (logic compiler) " software realizes that software compiler used is similar when it writes with program development,
And the source code before compiling also write by handy specific programming language, this is referred to as hardware description language
(Hardware Description Language, HDL), and HDL is also not only a kind of, but have many kinds, such as ABEL
(Advanced Boolean Expression Language)、AHDL(Altera Hardware Description
Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL
(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby
Hardware Description Language) etc., VHDL (Very-High-Speed are most generally used at present
Integrated Circuit Hardware Description Language) and Verilog2.Those skilled in the art
It will be apparent to the skilled artisan that only needing to method flow slightly programming in logic and being programmed into integrated circuit with above-mentioned several hardware description languages
In, it is possible to it is readily available the hardware circuit for realizing the logical method flow.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
Realized by the mode of software plus required general hardware platform.Understood based on such, the technical scheme essence of the application
On the part that is contributed in other words to prior art can be embodied in the form of software product, the computer software product
It can be stored in storage medium, such as ROM/RAM, magnetic disc, CD, including some instructions are to cause a computer equipment
(can be personal computer, server, or network equipment etc.) performs each embodiment of the application or embodiment
Method described in some parts.
Each embodiment in this specification is described by the way of progressive, identical similar between each embodiment
Part mutually referring to what each embodiment was stressed is the difference with other embodiment.Especially, pin
For the embodiment of device, the introduction control for being referred to the embodiment of preceding method is explained.
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes performing particular task or realizes routine, program, object, the group of particular abstract data type
Part, data structure etc..The application can also be put into practice in a distributed computing environment, in these DCEs, by
Remote processing devices connected by communication network perform task.In a distributed computing environment, program module can be with
Positioned at including in the local and remote computer-readable storage medium including storage device.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application has many deformations
With change without departing from spirit herein, it is desirable to which appended claim includes these deformations and changed without departing from the application
Spirit.
Claims (10)
1. a kind of extracting method of angular-trace gather, it is characterised in that methods described includes：
Geological data and speed data are obtained, and determines based on the geological data and the speed data incidence angle of ray
Degree；
Based on the speed data and default ray parameter, the Fresnel-zone width of the ray is determined；
Enter row interpolation weighting to the geological data in the Fresnel-zone width range, and the result that interpolation is weighted is as described
The corresponding angular-trace gather of incident angle.
2. according to the method described in claim 1, it is characterised in that determine to penetrate based on the geological data and the speed data
The incident angle of line includes：
Set up default the first relation between ray parameter and incident angle；Wherein, interval velocity is also included in first relation；
The second relation when setting up the default ray parameter and offset distance, travelling, between layer root mean sequare velocity；
Based on first relation and second relation, the incident angle of ray is calculated.
3. method according to claim 2, it is characterised in that determine first relation according to the following equation：
<mrow>
<mi>P</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
</mrow>
<msub>
<mi>V</mi>
<mi>int</mi>
</msub>
</mfrac>
</mrow>
Wherein, P represents the default ray parameter, and θ represents the incident angle, VintRepresent the interval velocity.
4. method according to claim 2, it is characterised in that determine second relation according to the following equation：
<mrow>
<mi>P</mi>
<mo>=</mo>
<mfrac>
<mi>X</mi>
<mrow>
<msubsup>
<mi>tV</mi>
<mrow>
<mi>r</mi>
<mi>m</mi>
<mi>s</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
Wherein, P represents the default ray parameter, and X represents the offset distance, when t represents the travelling, VrmsRepresent the layer
Root mean sequare velocity.
5. method according to claim 2, it is characterised in that calculate the incident angle of ray according to the following equation：
<mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>V</mi>
<mi>int</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
<msubsup>
<mi>V</mi>
<mrow>
<mi>r</mi>
<mi>m</mi>
<mi>s</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<mrow>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mrow>
<mi>r</mi>
<mi>m</mi>
<mi>s</mi>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>t</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
Wherein, θ represents the incident angle, and X represents the offset distance, VintRepresent the interval velocity, VrmsRepresent that the layer is equal
Root speed, t0When representing zero shot-geophone distance bilayer travelling.
6. according to the method described in claim 1, it is characterised in that determine the Fresnel bandwidth of the ray according to the following equation
Degree：
<mrow>
<mi>W</mi>
<mo>=</mo>
<msub>
<mi>V</mi>
<mrow>
<mi>r</mi>
<mi>m</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<mi>P</mi>
<mo>&CenterDot;</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
Wherein, W represents the Fresnel-zone width, VrmsThe layer root mean sequare velocity is represented, P represents the default ray parameter,
t0When representing zero shot-geophone distance bilayer travelling, f0Represent the basic frequency of the geological data.
7. method according to claim 6, it is characterised in that enter to the geological data in the Fresnel-zone width range
Row interpolation weighting includes：
For all recorded traces in the range of outgoing in t0The amplitude at moment, by the distance between the recorded trace and eye point
Amplitude corresponding weight factor reciprocal as the recorded trace, and carried out based on the amplitude and corresponding weight factor
Weighted superposition.
8. a kind of extraction element of angular-trace gather, it is characterised in that described device includes：
Incident angle determining unit, for obtaining geological data and speed data, and based on the geological data and the speed
Degrees of data determines the incident angle of ray；
Fresnel-zone width determining unit, for based on the speed data and default ray parameter, determining the ray
Fresnel-zone width；
Angular-trace gather determining unit, for entering row interpolation weighting to the geological data in the Fresnel-zone width range, and will
The result of interpolation weighting is used as the corresponding angular-trace gather of the incident angle.
9. device according to claim 8, it is characterised in that the incident angle determining unit includes：
First relation sets up module, for setting up the first relation between default ray parameter and incident angle；Wherein, described
Also include interval velocity in one relation；
Second relation sets up module, during for setting up the default ray parameter and offset distance, travelling, between layer root mean sequare velocity
The second relation；
Incident angle computing module, for based on first relation and second relation, calculating the incident angle of ray.
10. device according to claim 9, it is characterised in that the incident angle computing module is counted according to the following equation
Calculate the incident angle of ray：
<mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>&theta;</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>V</mi>
<mi>int</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
<msubsup>
<mi>V</mi>
<mrow>
<mi>r</mi>
<mi>m</mi>
<mi>s</mi>
</mrow>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mfrac>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<mrow>
<msup>
<mi>X</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>V</mi>
<mrow>
<mi>r</mi>
<mi>m</mi>
<mi>s</mi>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>t</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
Wherein, θ represents the incident angle, and X represents the offset distance, VintRepresent the interval velocity, VrmsRepresent that the layer is equal
Root speed, t0When representing zero shot-geophone distance bilayer travelling.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710628826.6A CN107329169A (en) | 2017-07-28 | 2017-07-28 | The extracting method and device of a kind of angular-trace gather |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710628826.6A CN107329169A (en) | 2017-07-28 | 2017-07-28 | The extracting method and device of a kind of angular-trace gather |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107329169A true CN107329169A (en) | 2017-11-07 |
Family
ID=60200320
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710628826.6A Pending CN107329169A (en) | 2017-07-28 | 2017-07-28 | The extracting method and device of a kind of angular-trace gather |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107329169A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110531414A (en) * | 2019-08-21 | 2019-12-03 | 中国矿业大学 | A kind of reflection seismic tomography fine detection method at high dip angle multilayer interface |
| CN112748466A (en) * | 2019-10-30 | 2021-05-04 | 中国石油天然气集团有限公司 | Travel time field data processing method and device based on Fresnel body |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020042678A1 (en) * | 2000-08-07 | 2002-04-11 | Dimitri Bevc | Typing picks to horizons in migration velocity analysis |
| CN104268412A (en) * | 2014-09-29 | 2015-01-07 | 中国石油天然气股份有限公司 | Angle gather ray chromatography migration velocity analyzing method and device |
| CN104360387A (en) * | 2014-10-20 | 2015-02-18 | 李晓峰 | Diffracted wave separating and imaging method based on reflected energy prediction |
| US20160061976A1 (en) * | 2014-04-17 | 2016-03-03 | Saudi Arabian Oil Company | Generating Subterranean Imaging Data Based on Vertical Seismic Profile Data |
-
2017
- 2017-07-28 CN CN201710628826.6A patent/CN107329169A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020042678A1 (en) * | 2000-08-07 | 2002-04-11 | Dimitri Bevc | Typing picks to horizons in migration velocity analysis |
| US20160061976A1 (en) * | 2014-04-17 | 2016-03-03 | Saudi Arabian Oil Company | Generating Subterranean Imaging Data Based on Vertical Seismic Profile Data |
| CN104268412A (en) * | 2014-09-29 | 2015-01-07 | 中国石油天然气股份有限公司 | Angle gather ray chromatography migration velocity analyzing method and device |
| CN104360387A (en) * | 2014-10-20 | 2015-02-18 | 李晓峰 | Diffracted wave separating and imaging method based on reflected energy prediction |
Non-Patent Citations (3)
| Title |
|---|
| 大港油田科技丛书编委会 编: "《地震勘探资料处理和解释技术》", 30 September 1999, 石油工业出版社 * |
| 张广智 等: ""基于射线理论的角度道集提取方法综述"", 《地球物理学进展》 * |
| 王文龙 等: ""基于射线理论的角度道集提取方法研究"", 《中国地球物理2011》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110531414A (en) * | 2019-08-21 | 2019-12-03 | 中国矿业大学 | A kind of reflection seismic tomography fine detection method at high dip angle multilayer interface |
| CN110531414B (en) * | 2019-08-21 | 2020-10-30 | 中国矿业大学 | Reflection seismic fault fine detection method for high-dip-angle multilayer interface |
| CN112748466A (en) * | 2019-10-30 | 2021-05-04 | 中国石油天然气集团有限公司 | Travel time field data processing method and device based on Fresnel body |
| CN112748466B (en) * | 2019-10-30 | 2024-03-26 | 中国石油天然气集团有限公司 | Fresnel-based travel time field data processing method and device |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105675635B (en) | Tight rock component relative amount and brittleness index determine method and apparatus | |
| CN104237940B (en) | A kind of diffraction wave imaging method based on dynamic characteristic and device | |
| CN104375188A (en) | Seismic wave transmission attenuation compensation method and device | |
| CN108415073B (en) | Angle domain back scattering offset imaging method and device | |
| CN107255831A (en) | A kind of extracting method of prestack frequency dispersion attribute | |
| CN107765298A (en) | A kind of method and device for determining near-surface velocity model | |
| CN104459768A (en) | Three-dimensional space target geologic body tracking method based on visualization | |
| CN107817526A (en) | Prestack seismic gather segmented amplitude energy compensation method and system | |
| CN107728204A (en) | Based on the anisotropic crack prediction method of prestack compressional wave and system | |
| CN109188520A (en) | Thin reservoir thickness prediction method and device | |
| CN104375178B (en) | A kind of carbonate rock fractured cave method for predicting reservoir and device | |
| CN103869362A (en) | Method and equipment for obtaining body curvature | |
| EP3749987B1 (en) | Systems and methods to enhance 3-d prestack seismic data based on non-linear beamforming in the cross-spread domain | |
| US11561312B2 (en) | Mapping near-surface heterogeneities in a subterranean formation | |
| CN102866426A (en) | Method for analyzing oil gas information of rock body by applying amplitude versus offset (AVO) large-angle trace gathers | |
| US10969507B2 (en) | Seismic channel geobody extraction for basin fairway simulation | |
| CN107329169A (en) | The extracting method and device of a kind of angular-trace gather | |
| CN107576985A (en) | A kind of method and apparatus of seismic inversion | |
| CN105093300B (en) | A kind of boundary recognition of geological body method and device | |
| CN106707345B (en) | A kind of angle elastic parameter Lithology Identification Methods and device | |
| Liu et al. | Combined petrophysics and 3D seismic attributes to predict shale reservoirs favourable areas | |
| CN106772599B (en) | A kind of method and device calculating formation shear speed | |
| CN106501851A (en) | A kind of optimum methods of seismic attributes and device | |
| CN107991705B (en) | Log curve correction method and apparatus based on Two-dimensional Statistical feature | |
| CN105467452B (en) | A kind of method and device that static correction is carried out to transitional region |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20171107 |
|
| WD01 | Invention patent application deemed withdrawn after publication |