CN109725350A - A kind of method, apparatus and system of determining longitudinal wave azimuthal anisotropy parameter - Google Patents

Abstract

This specification embodiment discloses the method, apparatus and system of a kind of determining longitudinal wave azimuthal anisotropy parameter, the azimuthal anisotropy parameter includes an equivalent Circle Parameters and two elliptic parameters, wherein, the equivalent Circle Parameters include the radius of equivalent circle speed of velocity ellipse, and the elliptic parameter includes the short long axis speed ratio of velocity ellipse, the long axis angle of velocity ellipse.The described method includes: filtering out the channel set data of present analysis point from the prestack trace gather data in target work area；It is scanned according to the channel set data interaction and obtains normal-moveout spectrum, the speed parameter of present analysis point is determined according to the energy extreme point of the normal-moveout spectrum；Then, interaction scanning obtains angular spectrum, determines that the angle parameter of present analysis point, and interaction scanning obtain speed comparison according to the energy extreme point of the angular spectrum, the speed ratio parameter of present analysis point is determined according to the energy extreme point of the speed comparison.

Description

A kind of method, apparatus and system of determining longitudinal wave azimuthal anisotropy parameter

Technical field

The present invention relates to petroleum natural gas exploration technical fields, particularly, be related to a kind of determining longitudinal wave orientation respectively to The method, apparatus and system of Anisotropic parameter.

Background technique

In petroleum exploration field, need to study the generation, migration and preservation rule of petroleum resources.In oil-gas reservoir growth course In, oil gas can realize repeatedly migration by crack.Crack Detection is one of important topic of oil-gas exploration.Seismic prospecting utilizes people The geophysical property of work earthquake-wave-exciting Underground medium.In wide-azimuth seismic prospecting, due to fracture development, usually can Observe that the phenomenon that speed of seismic wave changes with observed azimuth, and the speed of this seismic wave changes with observed azimuth is known as Azimuthal anisotropy phenomenon.

Speed caused by vertical fracture is described with the regular available velocity ellipse of azimuthal variation, these describe velocity ellipse The parameter of geometrical characteristic is known as azimuthal anisotropy parameter.The velocity analysis technology of early stage directly ignores depositing for azimuthal anisotropy Wide-azimuth three dimensional seismic data is being handled under the premise of the assumption of isotropy in orientation.This simplified technology will lead to image blur. If significantly cutting off remote geophone offset seismic channel, the focusing effect of imaging can be improved to a certain extent；But significantly Excision necessarily reduces the degree of covering of imaging point again, to reduce the confidence level of imaging results.

The classical way for obtaining azimuthal anisotropy parameter is that point sector carries out velocity analysis.This method is needed earthquake Data are divided into multiple azimuth sectors by azimuthal distribution situation, then carry out velocity analysis to each sector, finally use institute There is the velocity fitting of sector to go out azimuthal anisotropy parameter.The shortcomings that this method is there are in terms of following three: each sector Data require to carry out Manual Speed's analysis respectively, a large amount of manpowers need to be expended；Due to the irregular azimuth distribution of seismic data, The degree of covering of certain sectors is lower, will affect the confidence level of final fitting result；Divide sector method by the azimuth information in sector It is averaged, significantly reduces bearing accuracy, to influence the precision of the azimuthal anisotropy parameter fitted.

In order to break through elliptical velocity it is assumed that calculate the speed of arbitrarily complicated changing rule, in above-mentioned point of sector speed On the basis of degree analysis and the fitting of more sector parameters, it is thus proposed that replaces the technical thought being fitted with interpolation.This technology tool Standby two features: first is that sector is divided to carry out velocity analysis；Second is that parabola interpolation is carried out with the velocity analysis result of adjacent sectors, Calculate the velocity amplitude of arbitrary orientation angle.The shortcomings that this technology, is: having used the velocity information interpolation meter of adjacent three sectors Calculate arbitrary orientation angular speed, the signal-to-noise ratio of these three sectors of the confidence level heavy dependence of calculated result.In addition to this, this skill Art still has the shortcomings that three of aforementioned point of sector art.

In order to solve the problems, such as that point sector velocity analysis and the fitting of more sector parameters encounter, develop again later comprehensive Anisotropic velocity inversion technique.This technology carries out orientation isotropism velocity analysis first, then with isotropism speed Degree carries out dynamic correction or pre-stack time migration operation to prestack trace gather；Residual move out time is counted with correlation method again；Then with minimum Quadratic method fits three elliptic parameters；Residual move out time is carried out to imaging trace gather with three elliptic parameters and isotropism speed again Correction.The problem of although this technology part solves the problems, such as that a point sector velocity analysis encounters, and there are still three aspects: When noise is relatively low, it is difficult to count residual move out time with correlation method.Because residual move out time statistics and elliptic parameter fitting operation be It is carried out before elliptic parameter application, application effect can not be predicted in statistics and fit procedure, it is difficult to determine that Optimal calculation is joined Number.When the application effect of elliptic parameter is undesirable, it is difficult to improve the time difference statistics of next round and parameter fitting operation.

Summary of the invention

The purpose of this specification embodiment is to provide a kind of method, apparatus of determining longitudinal wave azimuthal anisotropy parameter and is System, can more efficiently determine longitudinal wave azimuthal anisotropy parameter.

This specification provides the method, apparatus and system of a kind of determining longitudinal wave azimuthal anisotropy parameter, implementation May include:

A kind of method of determining longitudinal wave azimuthal anisotropy parameter, which comprises

The azimuthal anisotropy parameter includes an equivalent Circle Parameters and two elliptic parameters, wherein

The equivalent Circle Parameters include the radius of equivalent circle speed of velocity ellipse, and the radius of equivalent circle speed is for determining The amplitude of the low frequency component of dynamic correction value；

The elliptic parameter includes the short long axis speed ratio of velocity ellipse, the long axis angle of velocity ellipse；Wherein, short long axis Speed ratio is used to determine the amplitude of the high fdrequency component of dynamic correction value, and long axis angle is used to determine the phase of the high fdrequency component of dynamic correction value Position；

The longitudinal wave azimuthal anisotropy parameter is determined in a manner of interaction analysis, comprising:

The channel set data of present analysis point are filtered out from the prestack trace gather data in target work area；

Determine the isotropic radius of equivalent circle speed in orientation of the present analysis point, comprising: according to default scanning speed It spends section and isotropism velocity scanning in orientation is carried out to the channel set data, obtain the orientation of the present analysis point respectively to same Property normal-moveout spectrum, the orientation isotropic radius of equivalent circle speed of present analysis point is determined according to the energy extreme point of the normal-moveout spectrum Degree；

Based on the isotropic radius of equivalent circle speed in the orientation, the short long axis speed ratio of the present analysis point is determined With long axis angle, comprising: the initial short long axis speed ratio based on the preset present analysis point, according to default scanning angle area Between interaction scanning angle spectrum, the long axis angle of the present analysis point is determined according to the energy extreme point of the angular spectrum；And Based on the long axis angle of the present analysis point, according to default scanning speed than section interaction scanning speed comparison, according to described The energy extreme point of speed comparison determines the short long axis speed ratio of the present analysis point.

In another embodiment of the method that this specification provides, the short long axis speed ratio of the determination and major axis angle Degree, comprising:

Based on the isotropic radius of equivalent circle speed in the orientation, work as described in the method determination using interactive iteration analysis The short long axis speed ratio and long axis angle of preceding analysis site, wherein the method for interactive iteration analysis includes:

Based on the initial short long axis speed ratio of the preset present analysis point, swept according to the interaction of default scanning angle section Angular spectrum is retouched, the long axis angle of the present analysis point is determined according to the energy extreme point of the angular spectrum, and, it executes following Iterative step:

Based on the long axis angle for the present analysis point that last iteration determines, according to default scanning speed than section interaction Scanning speed comparison determines the short long axis speed ratio of the present analysis point according to the energy extreme point of the speed comparison；

Based on the short long axis speed ratio for the present analysis point that current iteration determines, handed over according to default scanning angle section Mutual scanning angle spectrum, the long axis angle of the present analysis point is determined according to the energy extreme point of the angular spectrum；

Until meeting default iteration requirement, the short long axle speed for the present analysis point that output last time iteration obtains When long axis angle.

In another embodiment of the method that this specification provides, the method also includes:

Determine the radius of equivalent circle speed of the azimuthal anisotropy of the present analysis point, and, it is each based on the orientation The radius of equivalent circle speed of anisotropy determines the short long axis speed ratio and long axis angle of the present analysis point, wherein

The radius of equivalent circle speed of the azimuthal anisotropy of the determination present analysis point, comprising:

According to default scanning speed section, the short long axis speed ratio and long axis angle of present analysis point, to described super Trace gather data carry out azimuthal anisotropy velocity scanning, obtain the azimuthal anisotropy normal-moveout spectrum of the present analysis point, according to The energy extreme point of the azimuthal anisotropy normal-moveout spectrum determines the radius of equivalent circle speed of the azimuthal anisotropy of present analysis point Degree.

In another embodiment of the method that this specification provides, the orientation of the determination present analysis point is each To the same sex radius of equivalent circle speed or determine the present analysis point azimuthal anisotropy radius of equivalent circle speed, Include:

It is determined according to the low frequency component amplitude size of the channel set residual normal moveout of the present analysis point described current The equivalent circular of the azimuthal anisotropy of the isotropic radius of equivalent circle speed in analysis site orientation or present analysis point half Diameter speed.

In another embodiment of the method that this specification provides, the short long axis of the determination present analysis point Speed ratio, comprising:

It is determined according to the high fdrequency component amplitude size of the channel set residual normal moveout of the present analysis point described current The short long axis speed ratio of analysis site.

This specification provide the method another embodiment in, the preset present analysis point it is initial Short long axis speed ratio is determined using following manner:

It is described current according to the setting of the high fdrequency component amplitude size of the channel set residual normal moveout of the present analysis point The initial short long axis speed ratio of analysis site.

In another embodiment of the method that this specification provides, the major axis angle of the determination present analysis point Degree, comprising:

Present analysis is determined according to the high fdrequency component phase size of the channel set residual normal moveout of the present analysis point The angle value of point.

In another embodiment of the method that this specification provides, the channel set number for filtering out present analysis point According to, comprising:

The channel set data that choice of parameters goes out present analysis point are sorted according to preset trace gather, the trace gather sorts parameter packet Include span parameter and parameters sortnig, wherein

The span parameter includes: adjacent lines number and adjacent C MP number；

The parameters sortnig includes: starting geophone offset, terminates geophone offset and geophone offset increment.

In another embodiment of the method that this specification provides, the channel set number for filtering out present analysis point According to, comprising:

The initial channel set data that parameter sub-elects present analysis point, the trace gather sorting ginseng are sorted according to preset trace gather Number includes span parameter and parameters sortnig, wherein the span parameter includes adjacent lines number and adjacent C MP number, the row Order parameter includes starting geophone offset, terminates geophone offset and geophone offset increment；

The initial channel set data sub-elected are grouped according to the parameters sortnig, by the trace gather data base in each group It is ranked up in the sequence that azimuth is incremented by, and each seismic channel data in initial channel set data is corrected to present analysis point On the datum level at place, the channel set data of present analysis point are obtained.

On the other hand, this specification embodiment also provides a kind of device of determining longitudinal wave azimuthal anisotropy parameter, comprising:

Channel set screening module, for filtering out the channel set of present analysis point from the prestack trace gather data in target work area Data；

Speed determination module, for carrying out orientation isotropism to the channel set data according to default scanning speed section Or azimuthal anisotropy velocity scanning, the orientation isotropism or azimuthal anisotropy normal-moveout spectrum of the present analysis point are obtained, The orientation isotropism of present analysis point or the equivalent circular of azimuthal anisotropy are determined according to the energy extreme point of the normal-moveout spectrum Radius speed；

Angle-determining module, for based on the orientation isotropism or azimuthal anisotropy radius of equivalent circle speed, And the initial short long axis speed ratio or the obtained present analysis of previous round iterative analysis of the preset present analysis point The short long axis speed ratio of point, according to default scanning angle section interaction scanning angle spectrum, according to the energy extreme value of the angular spectrum Point determines the long axis angle of the present analysis point；

Speed ratio determining module, for the radius of equivalent circle speed based on the orientation isotropism or azimuthal anisotropy The long axis angle of degree and present analysis point, according to default scanning speed than section interaction scanning speed comparison, according to institute The energy extreme point for stating speed comparison determines the short long axis speed ratio of the present analysis point.

On the other hand, this specification embodiment also provides a kind of equipment of determining longitudinal wave azimuthal anisotropy parameter, described Equipment includes Cooperation controlling element, data management element, channel set data sorting element, that channel set moves correcting element, orientation is each Anisotropy parameter initialization element, energy road computing element, normal-moveout spectrum computing element, speed comparison computing element, angular spectrum meter Element, speed spectrum analysis element, speed are calculated than spectrum analysis element, angular spectrum analysis element and channel set Dynamically Announce element； Wherein,

The Cooperation controlling element is used to coordinate and control the cooperation and operation of whole equipment；The Cooperation controlling element tool Standby user instruction interface, for receiving the instruction control of user and receiving the processing parameter of user's offer；

The data management element realizes data input, output and operational parameter under the control of Cooperation controlling element Management function；

The channel set data sorting element is used to receive the channel set data sort instruction that Cooperation controlling element is sent, from Cooperation controlling element receives input data and required parameter, realizes the sorting of channel set data；

The channel set, which moves correcting element, to be existed as independent component or exists as built-in element；The channel set is dynamic The function of correcting element includes: the dynamic correction instruction for receiving other elements and sending, and receives dynamic correction calculation data from other elements And parameter；Complete the dynamic correction calculation of channel set；Operation result is fed back into other elements, wherein

In the presence of the channel set moves correcting element as independent component, function includes: to receive Cooperation controlling element hair The dynamic correction instruction come receives dynamic correction calculation data and parameter from Cooperation controlling element, completes the dynamic correction calculation of channel set, Operation result is fed back into Cooperation controlling element；

When the dynamic correcting element of the channel set is built in inside other elements, function includes: to receive host's element to send Dynamic correction instruction, receive dynamic correction calculation data and parameter from host's element, complete channel set and move correction calculation, by operation knot Fruit feeds back to host's element；

The azimuthal anisotropy parameter initialization element, the azimuthal anisotropy sent for receiving Cooperation controlling element Parameter initialization instruction, completes azimuthal anisotropy parameter initialization；

Energy road computing element is built-in element, and function includes: the energy road computations for receiving host's element, The channel set data after dynamic correction are received from host's element, summation and the square operation of channel set is completed, obtains a superposition energy Road is measured, calculated result is fed back into host's element；

Include two built-in elements inside the normal-moveout spectrum computing element: channel set moves correcting element and energy road Computing Meta Part；The normal-moveout spectrum computing element receives the normal-moveout spectrum computations that Cooperation controlling element is sent, and completes normal-moveout spectrum and calculates, and will Calculated result feeds back to Cooperation controlling element；

Include two built-in elements inside the speed comparison computing element: channel set moves correcting element and energy road calculates Element；The speed comparison computing element receives the speed comparison computations that Cooperation controlling element issues, and completes speed comparison It calculates, and calculated result is fed back into Cooperation controlling element；

The angular spectrum computing element includes two built-in elements: channel set moves correcting element and energy road computing element； The angular spectrum computing element receives the instruction of Cooperation controlling element, completes angular spectrum and calculates, and calculated result is fed back to association Make control element；

The speed spectrum analysis element completes display and the radius of equivalent circle of normal-moveout spectrum under the control of Cooperation controlling element The pickup of speed parameter；

The speed than spectrum analysis element under the control of Cooperation controlling element, complete speed comparison display and short long axis Pickup of the speed than parameter；

The angular spectrum analysis element completes display and the angle parameter of angular spectrum under the control of Cooperation controlling element It picks up；

The channel set Dynamically Announce element realizes the Dynamically Announce of channel set under the control of Cooperation controlling element.

On the other hand, this specification embodiment also provides a kind of system of determining longitudinal wave azimuthal anisotropy parameter, including At least one processor and memory for storing computer executable instructions, when the processor executes described instruction, The structure of the function or the equipment of the step of realizing any one above-mentioned embodiment the method or described device.

A kind of method, apparatus for determining longitudinal wave azimuthal anisotropy parameter that this specification one or more embodiment provides And system, it can use radius of equivalent circle speed, three parameters of long axis angle and short long axis speed ratio describe the big of velocity ellipse Small, shape and orientation.And it can be by gradually scanning the Parameter Spectrums of above three anisotropic parameters, according to the energy of Parameter Spectrum Extreme's distribution determines corresponding parameter value, and then realizes the interactive iteration analysis of three azimuthal anisotropy parameters.It utilizes Each embodiment of this specification, can accurately determine azimuthal anisotropy parameter, and then more efficiently analyzes vertical fracture and lead The speed of cause with azimuthal variation rule.

Detailed description of the invention

In order to illustrate more clearly of this specification embodiment or technical solution in the prior art, below will to embodiment or Attached drawing needed to be used in the description of the prior art is briefly described, it should be apparent that, the accompanying drawings in the following description only needle To some embodiments recorded in this specification；For those of ordinary skill in the art, it is not making the creative labor Under the premise of, it can also be obtained according to these attached drawings other attached drawings.In addition, though all attached drawings of this specification provide only symbol The black and white or gray scale schematic diagram of patent application text specification are closed, still, possible embodiment is not limited to black and white or grayscale image Shape user interface.In the design of actual product, exploitation, test, implementation or application process, colour can be used.Obviously Ground, those of ordinary skill in the art understand that the actually described black and white or the corresponding color graphics of gray scale graphic user interface are used The implementation result of family skins is more excellent.In the accompanying drawings:

Fig. 1 is the azimuthal anisotropy parameter schematic diagram in one embodiment that this specification provides；

Fig. 2 is a kind of process signal of the embodiment of the method for determining longitudinal wave azimuthal anisotropy parameter that this specification provides Figure；

Fig. 3 is a kind of another embodiment of the method for determining longitudinal wave azimuthal anisotropy parameter that this specification provides Flow diagram；

Fig. 4 is a kind of another embodiment of the method for determining longitudinal wave azimuthal anisotropy parameter that this specification provides Flow diagram；

Fig. 5 is the azimuthal anisotropy parameter interactive iteration analysis process signal in the example that this specification provides Figure；

Fig. 6 is the orientation isotropism normal-moveout spectrum schematic diagram in the example that this specification provides；

Fig. 7 is the pickup result schematic diagram of the orientation isotropism normal-moveout spectrum in the example that this specification provides；

Fig. 8 is the channel set signal after the dynamic correction of orientation isotropism speed in the example that this specification provides Figure；

Fig. 9 is that the short long axis speed in the example that this specification provides adjusts schematic diagram than the initial value of parameter；

Figure 10 is the long axis angular spectrum schematic diagram in the example that this specification provides；

Figure 11 is the pickup result schematic diagram of the long axis angular spectrum in the example that this specification provides；

Channel set after moving correction in the long axis angle pick process in an example that Figure 12 provides for this specification shows It is intended to；

Figure 13 is the short long axis speed comparison schematic diagram in the example that this specification provides；

Figure 14 is pickup result schematic diagram of the short long axis speed in the example that this specification provides than parameter；

Figure 15 is that the short long axis speed in the example that this specification provides corrects the road Hou Chao than moving in pick process Collect schematic diagram；

Figure 16 is a kind of module knot of the Installation practice for determining longitudinal wave azimuthal anisotropy parameter that this specification provides Structure schematic diagram；

Figure 17 is the azimuthal anisotropy parameter interactive iteration analytical equipment in one embodiment that this specification provides Logical construction schematic diagram；

Figure 18 is the graphic user interface schematic diagram in one embodiment that this specification provides.

Specific embodiment

In order to make those skilled in the art more fully understand the technical solution in this specification, below in conjunction with this explanation Attached drawing in book one or more embodiment carries out the technical solution in this specification one or more embodiment clear, complete Site preparation description, it is clear that described embodiment is only specification a part of the embodiment, instead of all the embodiments.Based on saying Bright book one or more embodiment, it is obtained by those of ordinary skill in the art without making creative efforts all The range of this specification example scheme protection all should belong in other embodiments.

Although present description provides as the following examples or method operating procedure, apparatus function or equipment shown in the drawings Structure, but based on it is conventional or without creative labor in the method, device, equipment or system may include it is more or Less operating procedure or modular unit after person part merges.In the step of there is no necessary causalities in logicality or structure In, the modular structure of the execution of these steps sequence or device is not limited to this specification embodiment or shown in the drawings executes sequence Or modular structure.The method or modular structure device, server or end product in practice are in application, can be according to Embodiment or method shown in the drawings or modular structure carry out sequence execution or it is parallel execute (such as parallel processor or The environment of multiple threads, the even implementation environment including distributed treatment, server cluster).

In wide-azimuth seismic prospecting, due to fracture development, it can usually observe that the speed of seismic wave becomes with observed bearing Change.The phenomenon that this observation changes with observed azimuth referred to as azimuthal anisotropy phenomenon.

In general, speed caused by vertical fracture is with the available ellipse radii with azimuthal variation of azimuthal changing rule Length describes, it is this it is oval be known as velocity ellipse, longitudinal wave azimuthal anisotropy parameter can be defined to describe velocity ellipse Size, shape and orientation.Currently, usually using transverse speed, ellipse short shaft speed and ellipse short shaft azimuth as side Position anisotropic parameters；But through actual verification, the set parameter and formation fracture feature do not have specific corresponding relationship, and the parameter Also it is not appropriate for interacting iterative analysis.

Correspondingly, this specification embodiment provides a set of new longitudinal wave azimuthal anisotropy parameter, this specification is implemented The longitudinal wave azimuthal anisotropy parameter of example can be set as the combination of an equivalent Circle Parameters and two elliptic parameters；Correspondingly, The equivalent Circle Parameters may include the radius of equivalent circle speed of velocity ellipse, and the elliptic parameter may include velocity ellipse Short long axis speed ratio, the long axis angle of velocity ellipse.

By above-mentioned parameter that this specification embodiment is set physical significance analysis it is found that the velocity ellipse etc. Effect radius of circle speed can describe elliptical size, can be used to determine the amplitude of the low frequency component of dynamic correction value.The speed Elliptical short long axis speed ratio can describe elliptical shape, can be used to determine the amplitude of the high fdrequency component of dynamic correction value.Institute The long axis angle for stating velocity ellipse can describe elliptical orientation, can be used to determine the phase of the high fdrequency component of dynamic correction value. It can be easier to realize the high and low frequency of dynamic correction value using the above-mentioned longitudinal wave azimuthal anisotropy parameter that this specification embodiment defines Separation improves the working efficiency of parameter interaction analysis.

Meanwhile the it will be recognized by those skilled in the art that three azimuthal anisotropies ginseng provided in this specification embodiment On the basis of number (radius of equivalent circle speed, long axis angle and short long axis speed ratio), by simple mathematic(al) manipulation Other azimuthal anisotropy parameter combinations are obtained (as included but is not limited to: equivalent diameter speed, short axle angle, major and minor axis Speed ratio etc.) it should all be included within the protection scope of this specification.

Radius of equivalent circle speed, long axis angle and short are utilized for the convenience of description, only giving in this specification embodiment Long axis speed is than the description as azimuthal anisotropy parameter and corresponding embodiment.But theoretically, as long as by this specification After dynamic updating formula in embodiment does simple mathematic(al) manipulation, above-mentioned azimuthal anisotropy parameter is provided using this specification Other variants, still available and the same or similar result of this specification embodiment.It will be apparent that in this specification reality On the basis of applying example, the derivative parameter deformable body by similar simple transformation or derivation should all be included in this specification Within protection scope.

To simplify the description, in this specification embodiment, the radius of equivalent circle speed of velocity ellipse can be simplified and is described Short long axis speed for speed, velocity ellipse is described as speed ratio than simplifying, and the long axis angle simplification of velocity ellipse is described as angle Degree.

Fig. 1 is the azimuthal anisotropy parameter schematic diagram in this specification embodiment.As shown in Figure 1, vertical fracture causes Speed can be described with the length of the ellipse radii with azimuthal variation with azimuthal changing rule, it is this oval to be known as speed Degree is oval.InLine in Fig. 1 indicates that main profile, CrossLine indicate interconnection.In Fig. 1 it is shown in solid ellipse in, V_sIt is ellipse Round semi-minor axis length, the elliptical semi-minor axis speed of representation speed；V_fFor elliptical major semiaxis length, the elliptical length of representation speed Half axle speed；α is the angle of transverse and InLine (main profile) direction, the elliptical long axis angle of representation speed.

The short long axis speed ratio γ of velocity ellipse is defined as the ratio of elliptical semi-minor axis speed and major semiaxis speed:

This circle can be known as ellipse by the area of circle shown in dotted lines in Figure 1 and elliptical area equation shown in solid Equivalent circular, the V in Fig. 1_eFor the radius length of equivalent circular, it is properly termed as radius of equivalent circle speed.I.e. velocity ellipse is equivalent Radius of circle speed V_eIt can indicate are as follows:

Fig. 2 is a kind of embodiment of the method process signal for determining longitudinal wave azimuthal anisotropy parameter that this specification provides Figure.Specific one embodiment as shown in Fig. 2, the method for the determination longitudinal wave azimuthal anisotropy parameter that this specification provides one In a embodiment, the method may include:

S2: the channel set data of present analysis point are filtered out from the prestack trace gather data in target work area.

The prestack trace gather data can be original prestack common point (CMP) trace gather data, when being also possible to prestack Between the prestack CMP trace gather data that are obtained after reaction correction calculation of migration imaging trace gather.

The present analysis point may include current point to be analyzed.In some embodiments, the present analysis point information May include analysis site wire size and No. CMP.In order to simplify subsequent step and process loop iteration, one of this specification or In multiple embodiments, present analysis point can be set by the point of first iterative analysis, and can determine the line of present analysis point Number and No. CMP, the analysis of other points to be analyzed can carry out on the basis of the analysis site, to prevent from going out during iterative analysis Existing entanglement.

The channel set refers to be merged by several different seismic channel sets (different wire sizes or different No. CMP trace gather) Trace gather；After completing to merge, the wire size of all seismic channels in a channel set and No. CMP are respectively set to present analysis Point wire size and No. CMP.After merging, the seismic channel data in a channel set has common wire size and No. CMP.Certainly, originally Channel set data described in specification embodiment, can also be only from the same seismic channel set, at this time without merging.

When it is implemented, the channel set number for needing voluntarily to filter out present analysis point that can be analyzed according to actual parameter According to.In one embodiment of this specification, the channel set that choice of parameters goes out present analysis point can be sorted according to preset trace gather Data.The preset trace gather sorting parameter may include span parameter and parameters sortnig, and the span parameter may include phase Adjacent line number and adjacent C MP number；The parameters sortnig may include starting geophone offset, terminate geophone offset and geophone offset increment. Wherein, the adjacent lines number can refer to the item number for participating in the adjacent lines of present analysis point processing；The adjacent C MP number can To refer to adjacent C MP number for participating in present analysis point processing.Certainly, when the adjacent line number for participating in present analysis point processing When mesh and adjacent C MP number are all 1, the trace gather data sub-elected are still referred to as the channel set data of present analysis point.

Using preset trace gather sorting parameter to screen channel set, it can filter and drop off the farther away data of present analysis point, The data near analysis site are only allowed to participate in operation.It can not only guarantee computational accuracy in this way, but also operation efficiency can be improved.

It, can also be according to the position of present analysis point and above-mentioned trace gather sorting ginseng in another embodiment of this specification Number sub-elects trace gather data from prestack trace gather data, obtains the initial channel set data of present analysis point.Then, then to sorting Initial channel set data out are grouped according to the parameters sortnig, and the trace gather data in each group are incremented by based on azimuth Sequence is ranked up.And each seismic channel data in the initial channel set data is corrected to the benchmark where present analysis point On face, the channel set data of present analysis point are obtained.

When it is implemented, can according to the position of present analysis point and sorting parameter by required data from prestack trace gather data In sort out, merged with the trace gather data of present analysis point position, formed present analysis point initial channel set data. It is then possible to which each seismic channel obtained after merging is grouped according to above-mentioned parameters sortnig, i.e., seismic channel is pressed into geophone offset Size be divided into different groups.It is then possible to which each seismic channel sorts according to the sequence that azimuth is incremented by organizing again.Using above-mentioned Packet sequencing mode can be shown with further satisfaction channel set data interaction and (draw channel set data on the computer screen Image) needs.Meanwhile using above-mentioned packet sequencing mode, channel set residual normal moveout can also be more intuitively shown Low frequency and high fdrequency component, to interact for user, parameter is picked up and parameter adjustment provides guidance.

Finally, on datum level where each seismic channel data obtained after merging to be corrected to present analysis point.The base Quasi- face can refer to the elevation plane of reference corresponding to each seismic channel zero-time in channel set.Each channel set has respectively different Datum level, the datum level of each seismic channel is identical in channel set.Transformation of data can be prevented by carrying out base-level correction, and simplification is subsequent Dynamic correction and Parameter Spectrum calculate.And it can on this basis, construct the travel-time equation of dynamic correction, carry out dynamic correction calculation and meter Calculate stack power spectrum.

S4: the isotropic radius of equivalent circle speed in orientation of the present analysis point is determined.

Isotropism velocity scanning in orientation can be carried out to the channel set data according to default scanning speed section, obtained The orientation isotropism normal-moveout spectrum of the present analysis point；Present analysis point is determined according to the energy extreme point of the normal-moveout spectrum The isotropic radius of equivalent circle speed in orientation.

In some embodiments, the scanning range of speed parameter can be preset；Wherein, the scanning range of speed parameter Include: starting velocity, terminate speed and speed number.In other embodiments, can also preset speed, speed ratio, The iteration initial value of three parameters of angle.

When determining the isotropic radius of equivalent circle speed in orientation, speed can be set as 1.0, angle than initial value in advance Initial value is set as 0 degree.The short long axis speed of velocity ellipse is than describing elliptical shape, when the value of speed ratio γ is 1.0, quite It degenerates in the shape of velocity ellipse to be round, corresponding calculate is based under the isotropic assumed condition in orientation.

When it is implemented, can determine orientation using following manner on the basis of speed ratio γ=1.0, angle [alpha]=0 Isotropic radius of equivalent circle speed:

S402: interaction scan position isotropism normal-moveout spectrum.

Velocity scanning is carried out to the channel set data using the speed data in default sweep interval, is obtained described current The normal-moveout spectrum of analysis site determines the speed parameter of present analysis point according to the energy extreme point of the normal-moveout spectrum.

The speed parameter in default sweep interval can be determined according to default starting velocity, termination speed and speed number, To give a series of scanning speed values, dynamic correction is carried out to channel set respectively with each scanning speed.It is dynamic to correct the travelling used Shi Fangcheng can be indicated are as follows:

Wherein, t is two-way travel time, t₀For the self excitation and self receiving time, x is geophone offset, V_eFor radius of equivalent circle speed.

It theoretically, is orientation isotropism normal-moveout spectrum using the normal-moveout spectrum that formula (3) is calculated.It can be according to given Parameter can use formula (3) and calculate each t₀It is worth corresponding t value, then the sampling point of t moment is moved to t₀Moment.Wherein, In the given parameter of dynamic timing, V_eFor each scanning speed value, geophone offset, azimuth are according to the channel set data of present analysis point It determines.

By each trace-stacking after dynamic correction at together, after square operation, the corresponding superposition of each scanning speed is obtained Energy road.Correspondingly, stack power can be calculated according to following formula:

Wherein, i is the serial number for inputting seismic channel, and N is seismic channel number in channel set, and j is the serial number of sampling point in seismic channel, M For number of samples in seismic channel, A_ijFor j-th of sample value of i-th of seismic channel after dynamic correction；S_jFor the stack power road of output J-th of sample value.

For any one scanning speed value, the S of its corresponding all sampled point is calculated_j, that is, obtain the scanning speed value Stack power road.

The corresponding stack power road of each scanning speed is arranged according to the sequence of speed from small to large again, obtains present analysis The normal-moveout spectrum of point.

S404: the isotropic radius of equivalent circle speed in orientation is picked up in interaction.

It is referred to the energy extreme point distribution situation of normal-moveout spectrum, adds new velocity pick point or the existing speed of adjustment Pickup point makes the longitudinally varying trend of velocity pick point meet the regularity of distribution of normal-moveout spectrum energy extreme point.The pickup point can To refer to a point of the interaction pickup on normal-moveout spectrum, each velocity pick o'clock is described by a time value and a velocity amplitude. It is thus possible to which the distribution according to energy extreme point on normal-moveout spectrum, obtains the speed ginseng of present analysis point in such a way that interaction is picked up Number.During the pickup of speed parameter and adjustment, formula (3) can also be used to carry out dynamic correction in real time to channel set, to sentence Whether the disconnected speed parameter picked up is reasonable.

Scheme in above-described embodiment, normal-moveout spectrum calculates and the dynamic correction of channel set is all based on the orientation assumption of isotropy；When The speed parameter of preceding acquisition is only capable of substantially evening up the reflection line-ups of channel set, and caused by azimuthal anisotropy when travelling Disturbance quantity exists as residual normal moveout.

It, can be according to the low of the channel set residual normal moveout of the present analysis point in one embodiment of this specification Frequency component amplitude size determines the isotropic radius of equivalent circle speed in present analysis point orientation.

, as specified parameter, other specified parameters can be kept not for according to the velocity amplitude of normal-moveout spectrum energy extreme value dot picking Become, dynamic correction is carried out using formula (3) to the channel set of present analysis point, the channel set data after obtaining dynamic correction.Then, may be used To analyze the low frequency component amplitude size of the channel set data after dynamic correction, determine whether the velocity amplitude picked up is reasonable.If super The low frequency component amplitude of the residual normal moveout of trace gather is larger, then illustrating the velocity amplitude of current pickup, there are large error, Ying Chong The velocity amplitude that new adjustment is picked up.

Wherein, if determining that the standard of the low frequency component amplitude size of the residual normal moveout of channel set may include: dynamic The reflection interval of reflection line-ups after correction shows as one close to the horizontal curve of cyclical fluctuations, then low frequency component amplitude has been Very little.Conversely, then there is biggish low frequency component.

Using the above scheme of the present embodiment, by during speed parameter picks up and adjusts, to channel set data into Row dynamic correction in real time, and the amplitude size of the residual normal moveout low frequency component by analysis channel set, to determine the speed of pickup Whether reasonable spend parameter value.The above scheme of the present embodiment can be reduced because energy extreme value point focusing is not particularly evident or pole The not enough accurate bring velocity pick error that value point position determines, to further increase the equivalent circular half of present analysis point The accuracy that diameter speed determines improves data-handling efficiency.

S6: it is based on the isotropic radius of equivalent circle speed in the orientation, determines the short long axle speed of the present analysis point Degree ratio and long axis angle.

Under the premise of the isotropic radius of equivalent circle speed in orientation that can be determined in keeping above-mentioned steps is constant, really Fixed short long axis speed ratio and long axis angle parameter.

In some embodiments, speed ratio, the scanning range of angle parameter can be preset.Wherein, speed ratio parameter Scanning range include: starting velocity ratio, terminate speed ratio and speed than three parameters of number；The scanning range packet of angle parameter Include three start angle, termination point and angle number parameters.

When it is implemented, can be swept based on the initial short long axis speed ratio of the preset present analysis point according to default Angular interval interaction scanning angle spectrum is retouched, the major axis angle of the present analysis point is determined according to the energy extreme point of the angular spectrum Degree.Then, then the long axis angle based on the present analysis point, according to default scanning speed than section interaction scanning speed ratio Spectrum, the short long axis speed ratio of the present analysis point is determined according to the energy extreme point of the speed comparison.

It is then possible to determine long axis angle parameter using following manner:

S602: interaction scanning angle spectrum.

When it is implemented, the initial short long axis speed ratio of the present analysis point can be preset, in orientation respectively to same Property radius of equivalent circle speed and pre-set initial short long axis speed ratio on the basis of, according to default scanning angle area Between interaction scanning angle spectrum, the long axis angle of the present analysis point is determined according to the energy extreme point of the angular spectrum.

The initial short long axis speed is than that can be set greater than 0 and the value less than 1.0 previously according to actual needs.So It afterwards, can be according to the initial velocity of setting than the angle parameter in parameter, the speed parameter of present analysis point, default sweep interval Angle scanning is carried out to the channel set data, obtains the angular spectrum of the present analysis point.

Scanning step is as follows: according to angle scanning range, a series of scanning angles is given, with each scanning angle respectively to super Trace gather carries out dynamic correction, and dynamic correction equation is as follows accordingly:

Wherein, t is two-way travel time, t₀For the self excitation and self receiving time, x is geophone offset, V_eFor speed, γ is speed ratio, α For angle, β is the azimuth of seismic channel.

Each t can be calculated according to given parameter₀It is worth corresponding t value, then the sampling point of t moment is moved to t₀When It carves.Wherein, implement in the given parameter of dynamic timing, V_eSpeed parameter for the present analysis point determined in above-described embodiment, γ is initial velocity than parameter, and α is each scan angle angle value, and geophone offset, azimuth are true according to the channel set data of present analysis point It is fixed.

It is then possible to which each earthquake trace-stacking after dynamic correction after square operation, is obtained each scan angle at together Spend corresponding stack power road.When it is implemented, can use formula (4) calculates stack power.It later, can be again by each scanning The corresponding stack power road of angle is arranged according to the sequence of scanning angle from small to large, obtains the angular spectrum of present analysis point.

S604: interaction pick-up angles parameter.

It is referred to the distribution situation of angular spectrum energy extreme point, adds new angular pickup point or the existing angle of adjustment Pickup point makes the longitudinally varying trend of angular pickup point meet the regularity of distribution of angular spectrum energy extreme point.With reference to above-mentioned implementation Example, the angular pickup point can be by a time values and an angle value description.To according to the energy extreme point of angular spectrum Distributed interactive picks up the angle parameter for obtaining present analysis point.

Further, during the pickup of angle parameter and adjustment, formula (5) can be used, channel set is carried out in real time Dynamic correction, to judge whether the angle value for picking up or editing is reasonable.

It, can also be according to the channel set residual normal moveout of the present analysis point in one embodiment of this specification High fdrequency component phase size determines the angle value of present analysis point.

Other specified parameters can be kept not according to the angle value of angular spectrum energy extreme value dot picking as specified parameter Become, dynamic correction is carried out using formula (5) to the channel set of present analysis point, the channel set data after obtaining dynamic correction.Then, may be used To analyze the high fdrequency component phase size of the residual normal moveout of the channel set data after dynamic correction, determine that the angle value picked up is It is no reasonable.If the high fdrequency component phase of the residual normal moveout of channel set is larger, illustrate that the angle value of current pickup exists Large error should readjust the angle value of pickup.

Wherein, if determining that the standard of the high fdrequency component phase size of the residual normal moveout of channel set may include: dynamic The reflection time curve of reflection line-ups after correction close to Yu Xuanbo, then the high fdrequency component phase of residual normal moveout compared with Greatly；, whereas if the time curve of the reflection line-ups after dynamic correction is close to dextrorotation wave, then the high frequency division of residual normal moveout It is smaller to measure phase.

Using the above scheme of the present embodiment, by during angle parameter picks up and adjusts, to channel set data into Row dynamic correction in real time, and the phase size of the residual normal moveout high fdrequency component by analysis channel set, to determine the angle of pickup Whether reasonable spend parameter value.The above scheme of the present embodiment can be reduced because energy extreme value point focusing is not particularly evident or pole Value point position determine not enough precisely and bring angular pickup error, thus further increase present analysis point angle determine Accuracy, improve data-handling efficiency.

In some embodiments, short long axis speed ratio parameter can be determined using following manner:

S606: interaction scanning speed comparison.

The speed parameter and angle parameter of the present analysis point determined in available above-described embodiment, it is then possible to According to pre-set velocity than scanning range, the speed ratio parameter of default sweep interval is determined, to obtain series of scans speed ratio ginseng Number.

It is then possible to according to the speed ratio in the speed parameter of present analysis point, angle parameter and default sweep interval Argument sequence carries out speed ratio scanning to the channel set data, the speed comparison of the present analysis point is obtained, according to described The energy extreme point of speed comparison determines the speed ratio parameter of present analysis point.

Dynamic correction is carried out to channel set when it is implemented, formula (5) can be based respectively on each scanning speed ratio.It can be with According to given parameter, each t is calculated₀It is worth corresponding t value, then the sampling point of t moment is moved to t₀Moment.Wherein, specific real It applies in the given parameter of timing, V_eSpeed parameter for the present analysis point determined in above-described embodiment, α are above-mentioned implementation The angle parameter of the present analysis point determined in example, γ are each scanning speed ratio, and geophone offset, azimuth are according to present analysis point Channel set data determine.

It is then possible to which each earthquake trace-stacking after dynamic correction after square operation, is obtained each scanning speed at together Degree is than corresponding stack power road.When it is implemented, can use formula (4) calculates stack power.Later, it will can respectively sweep again It retouches speed and is arranged according to scanning speed than sequence from small to large than corresponding stack power road, obtain the speed of present analysis point Comparison.

S608: interaction pickup velocity compares parameter.

Distribution situation referring to speed than spectrum energy extreme point adds new speed speed more existing than pickup point or adjustment Than pickup point, make speed than the regularity of distribution that the longitudinally varying trend of pickup point meets speed than spectrum energy extreme point.With reference to upper Embodiment is stated, the speed can be by a time value and a velocity ratio description than pickup point.To according to speed comparison Energy extreme point distributed interactive pick up obtain present analysis point speed ratio parameter.

Further, during pickup and adjustment of the speed than parameter, formula (5) can be used, channel set is carried out in fact When dynamic correction, judge to pick up or whether the velocity ratio of editor reasonable.

In one embodiment of this specification, the short long axis speed ratio of the determination present analysis point can also be wrapped It includes:

Described in being determined according to the high fdrequency component amplitude size of the channel set residual normal moveout of the present analysis point The short long axis speed ratio of present analysis point.

Speed in formula (5) can be appointed as to the speed ratio according to speed than spectrum energy extreme value dot picking than parameter Value keeps other specified parameter constants, carries out dynamic correction using formula (5) to the channel set of present analysis point, obtains dynamic correction Channel set data afterwards.It is then possible to analyze the high fdrequency component amplitude of the residual normal moveout of the channel set data after dynamic correction Size determines whether the velocity ratio picked up is reasonable.If the high fdrequency component amplitude of the residual normal moveout of channel set is larger, Illustrating the velocity ratio of current pickup, there are large errors, should readjust the velocity ratio of pickup.

Wherein, if determining that the standard of the high fdrequency component amplitude size of the residual normal moveout of channel set may include: dynamic The high fdrequency component amplitude of the reflection time curve of reflection line-ups after correction is larger, then the high fdrequency component of residual normal moveout Amplitude is larger；Conversely, then the high fdrequency component amplitude of residual normal moveout is smaller.

Using the above scheme of the present embodiment, by during speed is picked up and adjusted than parameter, to channel set data Dynamic correction in real time, and the amplitude size of the residual normal moveout high fdrequency component by analysis channel set are carried out, to determine pickup Whether speed is more reasonable than parameter value.The above scheme of the present embodiment can reduce because energy extreme value point focusing be not it is particularly evident or Person's extreme point position it is determining not enough precisely and bring speed ratio picking error, to further increase the speed of present analysis point Degree improves data-handling efficiency than determining accuracy.

In another embodiment of this specification, the initial of the present analysis point can be preset using following manner Short long axis speed ratio:

It is described current according to the setting of the high fdrequency component amplitude size of the channel set residual normal moveout of the present analysis point The initial short long axis speed ratio of analysis site.

It can be first according to the current orientation isotropism radius of equivalent circle speed of present analysis point, initial velocity ratio (1.0), initial angle parameter (0 degree) carries out dynamic correction to the channel set data of present analysis point.When it is implemented, can first by As specified parameter, reference is above-mentioned for current orientation isotropism radius of equivalent circle speed, initial velocity ratio, initial angle parameter Formula (5) carries out dynamic correction to channel set data.It is then possible to analyze the residue of lineups in the channel set data after dynamic correction The initial short long axis of the present analysis point is arranged according to high fdrequency component amplitude size for the high fdrequency component amplitude size of dynamic correction value Speed ratio.When it is implemented, can be by the value for being set as less than 1.0 and being greater than 0 by current value 1.0 of parametric speed ratio γ.Its In, determine that the mode of the high fdrequency component amplitude size of channel set residual normal moveout can be implemented with reference to above-described embodiment.It can be with According to the residual normal moveout high fdrequency component amplitude size of channel set, the size of γ value is set.The residual normal moveout of channel set High fdrequency component amplitude it is bigger, γ value is smaller.

Initial velocity ratio can also be replaced as specified ginseng using each speed ratio adjusted during adjustment Number keeps other specified parameter constants, carries out dynamic correction in real time using formula (5), analyzes corresponding residual normal moveout high frequency The size of component, so that the initial value of parametric speed ratio γ is gradually determined as a suitable value.

Before the present embodiment implementation, because the current value of parameter γ is 1.0 in formula (5), normal-moveout spectrum calculates and surpasses The dynamic correction of trace gather is all based on the orientation assumption of isotropy, and the speed parameter currently obtained is only capable of the reflection line-ups of channel set Disturbance quantity is as residual NMO correction when substantially evening up and (reduce dynamic correction value low frequency component), and travelling caused by azimuthal anisotropy Amount exists.It, can be preliminary through this embodiment by the value that current short long axis speed ratio settings are less than 1.0 and greater than 0 The high fdrequency component interference for reducing residual normal moveout caused by azimuthal anisotropy is conducive to improve subsequent angle and speed ratio ginseng The efficiency and accuracy of number adjustment.

Fig. 3 is a kind of another embodiment stream of the method for determining longitudinal wave azimuthal anisotropy parameter that this specification provides Journey schematic diagram.As shown in figure 3, it is isotropic equivalent to be also based on the orientation in another embodiment of this specification Radius of circle speed determines the short long axis speed ratio and long axis angle of the present analysis point in the way of interactive iteration, comprising:

Firstly, based on the initial short long axis speed ratio for presetting the present analysis point described in previous embodiment, according to Default scanning angle section interaction scanning angle spectrum, determines the present analysis point according to the energy extreme point of the angular spectrum Long axis angle；

Then, following iterative steps are executed:

Based on the long axis angle for the present analysis point that last iteration determines, according to default scanning speed than section interaction Scanning speed comparison determines the short long axis speed ratio of the present analysis point according to the energy extreme point of the speed comparison；

Based on the short long axis speed ratio for the present analysis point that current iteration determines, handed over according to default scanning angle section Mutual scanning angle spectrum, the long axis angle of the present analysis point is determined according to the energy extreme point of the angular spectrum；

Until meeting default iteration requirement, the short long axle speed for the present analysis point that output last time iteration obtains When long axis angle.

It can use the scheme of above-described embodiment, first primarily determine the short long axis speed ratio and long axis angle of present analysis point Parameter, then, then progressive alternate optimization, until result meets preset requirement.Wherein, the preset requirement may include: Speed comparison and angular spectrum energy extreme value point focusing are good and are reasonably distributed, and pickup point is located on energy extreme point, and channel set is same The high fdrequency component of the residual normal moveout of phase axis is minimized.

Wherein, if judging that the standard of the high fdrequency component size of the residual normal moveout of channel set may include: dynamic correction Reflection line-ups afterwards reflection time curve high fdrequency component amplitude it is larger, then the high fdrequency component of residual normal moveout compared with Greatly；Conversely, then the high fdrequency component of residual normal moveout is smaller.

The short long axis speed ratio of present analysis point and long axis angle determination method for parameter are with reference to upper in iterative optimization procedure Embodiment implementation is stated, is not described herein.

Fig. 4 is a kind of another embodiment stream of the method for determining longitudinal wave azimuthal anisotropy parameter that this specification provides Journey schematic diagram.As shown in figure 4, in another embodiment of this specification, the method can also include:

S8: the radius of equivalent circle speed of the azimuthal anisotropy of the present analysis point is determined, comprising:

According to default scanning speed section, the short long axis speed ratio and long axis angle of present analysis point, to described super Trace gather data carry out azimuthal anisotropy velocity scanning, obtain the azimuthal anisotropy normal-moveout spectrum of the present analysis point, according to The energy extreme point of the azimuthal anisotropy normal-moveout spectrum determines the radius of equivalent circle speed of the azimuthal anisotropy of present analysis point Degree；

S10: the short long axle speed of the present analysis point is determined based on the radius of equivalent circle speed of the azimuthal anisotropy Degree ratio and long axis angle.

It is possible, firstly, to the short long axis speed when major axis angle degree parameter determined in obtaining step S6, and according to default scanning Speed interval carries out dynamic correction to channel set using formula (5).I.e. according to the short long axis speed determined in S6 when major axis angle degree A series of scanning speed values in parameter and default scanning speed section are as given parameter, using formula (5) to super road Collection carries out dynamic correction.It is then possible to corresponding stack power be calculated using formula (4), to obtain azimuthal anisotropy speed Spectrum.

Then, pickup is interacted according to the azimuthal anisotropy normal-moveout spectrum, obtain the azimuthal anisotropy etc. Imitate radius of circle speed.Specific embodiment can be implemented with reference to the radius of equivalent circle velocity pick method in step S4, here It does not repeat them here.

It is then possible to keep the radius of equivalent circle speed of azimuthal anisotropy constant, determined in the way of interactive iteration The short long axis speed ratio and long axis angle of the present analysis point.Specific embodiment can be real with reference to the method in step S6 Apply, when angular spectrum or speed comparison scan, the corresponding isotropic radius of equivalent circle speed in orientation be substituted for orientation respectively to Anisotropic radius of equivalent circle speed.

Using above scheme provided in this embodiment, short long axis speed when major axis angle degree anisotropy ginseng can be being determined On the basis of number, radius of equivalent circle speed is advanced optimized, and using the radius of equivalent circle speed after optimization, further determine that pair The short long axis speed answered when major axis angle degree parameter.To improve the accuracy of finally determining parameter.

In some embodiments, it can determine the need for carrying out by analyzing the residual normal moveout of channel set data The optimization of radius of equivalent circle speed, and optimization number.Judge whether the residual normal moveout of channel set data meets default want The standard asked may include: that the low frequency component of the residual normal moveout of channel set lineups is minimized.

Later, can according to actual needs, be arranged next analysis site wire size and No. CMP, and the analysis site is arranged It for present analysis point, repeats the above steps, determines its corresponding anisotropic parameters.

Certainly, above-mentioned speed, the speed sequence that when angle parameter determines are only for example, on when it is implemented, It states the sequence that parameter determines not limit, the given parameters that corresponding dynamic correction uses also accordingly carry out proper transformation.For example, It can also first determine speed parameter, initial angle parameter is then given according to the trend of work area vertical fracture, is then moved Correction and speed comparison, which scan, determines speed ratio parameter, finally carries out dynamic correction than parameter according to new speed again and angular spectrum is swept It retouches, determines final angle parameter etc..

Fig. 5 indicates the detailed process schematic diagram of the realization azimuthal anisotropy parameter in one example of this specification.Such as Fig. 5 Shown, in order to better illustrate the practicability and feasibility of this specification scheme provided by the above embodiment, this specification is also mentioned A kind of specific example of scheme using above-described embodiment is supplied.

1, the position of first analysis site is set.

The wire size of first analysis site is set as 100, No. CMP is set as 200.And first analysis site (is used into wire size 100 and CMP Numbers 200 indicate) it is set as present analysis point.

2, setting trace gather sorts parameter.

It includes two classifications of span parameter and parameters sortnig that trace gather, which sorts parameter,.Span parameter includes that adjacent lines number (is set 1) and two parameters of adjacent C MP number (being set as 1) for；Parameters sortnig includes starting geophone offset (being set as 0 meter), terminates geophone offset (being set as 4500 meters) and three parameters of geophone offset increment (being set as 500 meters).

3, channel set data are sorted.

Parameter is sorted according to the position of present analysis point and trace gather, prestack trace gather data are sorted, are sorted and benchmark Face correction, forms the channel set data of present analysis point.

Steps are as follows: firstly, according to the position (wire size 100 and No. CMP 200) of analysis site and sorting parameter (adjacent line number Mesh 1,0 meter of geophone offset of starting, terminates 500 meters of 4500 meters of geophone offset and geophone offset increment at adjacent C MP number 1) by required data from Merging is sorted out in prestack trace gather data.Then, each seismic channel after merging is grouped according to above-mentioned packet parameters, Each seismic channel sorts according to the sequence that azimuth is incremented by organizing again.Finally, by each earthquake channel correction sub-elected to current point On datum level where analysis point, the channel set data of present analysis point are formed.

4, parameter scanning range is set.

Need to be arranged speed (the radius of equivalent circle speed of velocity ellipse), speed ratio (the short long axis speed ratio of velocity ellipse) and The scanning range of angle (the long axis angle of velocity ellipse) three parameters.Wherein, the scanning range of speed parameter includes: starting speed It spends (being set as 1800 meter per seconds), terminate speed (being set as 5400 meter per seconds) and speed number (being set as 121)；Scanning of the speed than parameter Range includes: starting velocity ratio (being set as 0.7), terminates speed ratio (being set as 1.0) and speed ratio number (being set as 31)；Angle parameter Scanning range include start angle (be set as -90 degree), termination point (being set as 90 degree), angle number (being set as 91).

5, the initial value of azimuthal anisotropy parameter is set.

Setting including following parameters initial value: velocity original value be set as 3600 meter per seconds, speed than initial value be set as 1.0, Angle initial value is set as 0 degree.

6, first scan normal-moveout spectrum.

(1800 meter per second of starting velocity terminates 5400 meter per second of speed, speed to the velocity scanning range set according to step 4 Number 121), a series of scanning speed (1800 meter per seconds, 1830 meter per seconds, 1860 meter per second ..., 5400 meter per second) is given, with first Beginning speed ratio, initial angle and each scanning speed carry out dynamic correction to channel set respectively.It is calculated when the travelling that dynamic correction uses public Formula reference formula (5).Because the current value of parameter γ is 1.0 in formula (5), the calculating of present speed spectrum is based on orientation respectively to same Property assume；At this point, theoretically, formula (5) can degenerate for formula (3).Therefore, optionally, dynamic to correct the trip used Calculation formula is also possible to formula (3) when row.

The way of dynamic correction is to calculate each t according to given parameter₀It is worth corresponding t value, then by t₀The sampling point at moment It moves to t moment.Each earthquake trace-stacking after dynamic correction after square operation, is obtained into each scanning speed at together again Corresponding stack power road.Energy balane is overlapped according to formula (4).

When having calculated all S_jAfterwards, the corresponding stack power road of a certain given scanning speed has just been obtained.

Again by the corresponding stack power road of each scanning speed according to scanning speed from small to large sequence (1800 meter per seconds, 1830 meter per seconds, 1860 meter per second ..., 5400 meter per second) arrangement, the normal-moveout spectrum of present analysis point is just obtained, as shown in Figure 6.

Abscissa is radius of equivalent circle speed in Fig. 6, and unit is meter per second；Ordinate is t₀Time, unit are millisecond；Fig. 6 The background graphics of middle irregular distribution indicate different time and the corresponding stack power of speed, and the cross wires in figure indicates part energy Measure extreme point position.Dot on white vertical line indicates the pickup point of radius of equivalent circle speed, is now 3600 meter per second of initial value. Because set when first scan normal-moveout spectrum the value of short long axis speed ratio γ as normal-moveout spectrum shown in 1, Fig. 6 be orientation isotropism Normal-moveout spectrum.

7, pickup velocity parameter.

The pickup result schematic diagram of Fig. 7 expression speed parameter.As shown in fig. 7, referring to the energy extreme's distribution of normal-moveout spectrum Situation adds new velocity pick point or the existing velocity pick point of adjustment, meets the longitudinally varying trend of velocity pick point The regularity of distribution of normal-moveout spectrum energy extreme point.Dot in Fig. 7 on white broken line indicates the pickup point of radius of equivalent circle speed, this The speed of secondary pickup is orientation isotropism speed.

During pickup velocity parameter, public affairs can be used under the premise of keeping the value of short long axis speed ratio γ to be 1 Formula (5) carries out dynamic correction in real time to channel set, to judge whether pickup velocity parameter is reasonable；Optionally, can also make Dynamic correction in real time is carried out to channel set with formula (3), to judge whether pickup velocity parameter is reasonable.Fig. 8 shows based on initial speed Degree ratio, initial angle and the orientation isotropism speed parameter of pickup carry out the channel set schematic diagram after dynamic correction to channel set. Abscissa is geophone offset in Fig. 8, and unit is rice；Ordinate is t₀Time, unit are millisecond.

Analysis chart 8 is it is found that t₀Though the back wave at 1100 milliseconds to 1600 milliseconds of the time is substantially evened up, but still is had similar The residual disturbance amount of dextrorotation wave.This explanation, the dynamic correction of orientation isotropism can even up the low frequency part of dynamic correction value, and with split The high frequency section for stitching related dynamic correction value exists as residual normal moveout.

8, default speed ratio parameter is adjusted.

Can regulate the speed ratio according to the size of the residual normal moveout high fdrequency component of channel set.By parameter γ from initial Value 1.0 is adjusted to less than 1.0 and to be greater than 0 value.The standard of adjustment may include: the high frequency division of the residual normal moveout of channel set Amount is bigger, and γ value is smaller.

If judging that the standard of the high fdrequency component amplitude size of the residual normal moveout of channel set may include: after moving correction Reflection line-ups reflection time curve high fdrequency component amplitude it is larger, then the high fdrequency component amplitude of residual normal moveout It is larger；, whereas if reflection interval shows as a more smooth curve or straight line, then the high fdrequency component vibration of residual normal moveout Width is smaller.At this point, moving calculation formula when the travelling that correction uses is formula (5).

γ value is adjusted to 0.98 herein, as shown in Figure 9.Fig. 9 indicates that initial value adjustment of the short long axis speed than parameter is shown It is intended to.Abscissa is short long axis speed ratio in Fig. 9, and unit is one thousandth；Ordinate is t₀Time, unit are millisecond；It does not advise Then background graphics indicate different time and speed than corresponding stack power, and cross wires indicates part energy level value point position.Speed Dot at than 0.98 on white vertical line indicates short long axle speed speed than value adjusted.The speed comparison of the Fig. 9 is basis (0 degree) of orientation isotropism speed parameter and initial angle angle value of pickup calculates determining.

Because being the orientation isotropism speed parameter picked up used in given parameters in the calculating process of this speed comparison With initial long axis angle (0 degree), the scanning result error of speed comparison is larger, and stack power is not restrained, and therefore, directly adjusts herein Whole (rather than being picked up referring to speed comparison) speed ratio parameter.

Disturbance quantity corresponding residual normal moveout when travelling caused by crack can be tentatively reduced using the scheme of this step, It improves successive iterations and optimizes the efficiency of determining angle and speed than parameter.

9, scanning angle is composed.

According to angle scanning range (degree of start angle -90,90 degree of termination point, angle number 91), a series of sweep is given It retouches angle (- 90 degree, -88 degree, -86 90 degree of degree ...), after being adjusted in the radius of equivalent circle velocity amplitude picked up based on step 7, step 8 Velocity ratio and above-mentioned each scanning angle dynamic correction is carried out to channel set respectively.It is calculated when the travelling that dynamic correction uses public Formula is formula (5).

The way of dynamic correction is to calculate each t according to given parameter₀It is worth corresponding t value, then by t₀The sampling point at moment It moves to t moment.After square operation, it is corresponding to be obtained into each scanning angle at together for each trace-stacking after dynamic correction again Stack power road.The calculation formula of stack power is formula (4).

Again by each scanning angle (- 90 degree, -88 degree, -86 90 degree of degree ...) corresponding stack power road according to angle from it is small to Big sequence arrangement, just obtains the angular spectrum of present analysis point.As shown in Figure 10.Figure 10 indicates long axis angular spectrum schematic diagram.Figure Abscissa is long axis angle in 10, and unit is degree；Ordinate is t₀Time, unit are millisecond；Irregular background graphics indicate not With time and the corresponding stack power of angle, cross wires indicates part energy level value point position；Dot on white vertical line indicates length The initial value (0 degree) of shaft angle degree.

10, pick-up angles parameter.

The energy extreme point distribution situation of angular referenced modal data adds new angular pickup point or the existing angle of adjustment Pickup point makes the longitudinally varying trend of angular pickup point meet the regularity of distribution of angular spectrum energy extreme point.Figure 11 indicates angle The pickup result schematic diagram of spectrum.Abscissa is long axis angle in Figure 11, and unit is degree；Ordinate is t₀Time, unit are millisecond； Irregular background graphics indicate that different time and the corresponding stack power of angle, cross wires indicate part energy level value point position；It is white Dot on color broken line indicates the pickup value of long axis angle.

During pick-up angles parameter, formula (5) can be used, dynamic correction in real time is done to channel set, to judge to pick up Or whether the angle value of editor is reasonable.Figure 12 indicates the radius of equivalent circle velocity amplitude picked up based on step 7, in step 8 after adjustment Velocity ratio and the angle parameter of above-mentioned pickup carry out the channel set schematic diagram after dynamic correction.Abscissa is big gun inspection in Figure 12 Away from unit is rice；Ordinate is t₀Time, unit are millisecond；Waveform indicates the back wave after dynamic correction.Compared with Fig. 8, figure T in 12₀Back wave at 1100 milliseconds to 1600 milliseconds of time is evened up substantially, and residual normal moveout is smaller.

11, scanning speed comparison.

According to speed ratio scanning range (starting ratio 0.7 terminates ratio 1.0 and ratio number 31), a series of sweep is given Retouching speed ratio, (0.7,0.71,0.72 ... 1.0), the long axis of the radius of equivalent circle velocity amplitude, step 10 pickup that are picked up based on step 7 Angle value and above-mentioned scanning speed compare channel set and carry out dynamic correction.Calculation formula is formula when the travelling that dynamic correction uses (5)。

The way of dynamic correction is: according to given parameter, calculating each t₀It is worth corresponding t value, then by t₀The sampling point at moment T moment is arrived in value removal of home.

Each trace-stacking after dynamic correction after square operation, is obtained into each scanning speed than corresponding at together again Stack power road.The calculation formula of stack power is formula (4).

Again by each scanning speed ratio (0.7,0.71,0.72 ... 1.0) corresponding stack power road according to speed ratio from it is small to Big sequence arrangement, just obtains the speed comparison of present analysis point, as shown in figure 13.

Figure 13 indicates the pickup result schematic diagram of short long axis speed comparison.Abscissa is short long axis speed ratio in Figure 13, single Position is one thousandth；Ordinate is t₀Time, unit are millisecond；Irregular background graphics indicate that different time and speed ratio are corresponding Stack power, cross wires indicate part energy level value point position.Speed than 0.98 at dot on white vertical line indicate short long axis Speed is than first value adjusted (value before adjustment is 1).Pickup has been used in the calculating process of the speed comparison of this step The long axis angle of radius of equivalent circle speed and pickup, the scanning result of this speed comparison is (compared with speed comparison shown in Fig. 9 Compared with) error significantly reduces, stack power convergence and reliable.

12, pickup velocity compares parameter.

Referring to the energy extreme point distribution situation of speed comparison, new speed speed more existing than pickup point or adjustment is added Than pickup point, make present speed than the regularity of distribution that the longitudinally varying trend of pickup point meets speed than spectrum energy extreme point.Figure 14 indicate the pickup result schematic diagram of short long axis speed ratio.Abscissa is short long axis speed ratio in Figure 14, and unit is one thousandth； Ordinate is t₀Time, unit are millisecond；Irregular background graphics indicate different time and speed than corresponding stack power, fork Line indicates part energy level value point position；Dot on white broken line indicates the current value of short long axis speed ratio.Relative to Figure 13, Only part pickup point is adjusted in Figure 14, so as to illustrate, after radius of equivalent circle velocity pick, and first just successive step Then speed ratio, then carries out the pickup of subsequent angle and speed than parameter, angle and speed can be improved than determining precision, And then subsequent angle, speed can be improved than efficiency that iteration updates.

During pickup and adjustment, real-time dynamic correction is carried out to channel set using formula (5), is picked up or is adjusted with judgement Velocity ratio it is whether reasonable.Figure 15 indicates the major axis angle that the radius of equivalent circle velocity amplitude picked up based on step 7, step 10 are picked up Angle value and the speed of above-mentioned pickup are than the channel set schematic diagram after carrying out dynamic correction.Abscissa is geophone offset in Figure 15, and unit is Rice；Ordinate is t₀Time, unit are millisecond；Waveform indicates to carry out the reflection after dynamic correction using azimuthal anisotropy parameter Wave.Comparison diagram 12 and Figure 15 can have found, channel set improvement effect of the channel set in Figure 15 compared with Figure 12 very little.

13, iteration renewal speed ratio and angle parameter.

9 to 12 steps, progressive alternate renewal speed ratio and angle parameter are repeated, until result is satisfactory.More by iteration After new, the energy extreme value point focusing of speed comparison and angular spectrum is good and is reasonably distributed, and pickup point is located on energy extreme point, remains The high fdrequency component of remaining dynamic correction value very little.

14, scanning speed is composed again.

It is first determined whether needing scanning speed spectrum again.Criterion may include: the residual normal moveout of channel set Low frequency component it is whether sufficiently small.If low frequency component is larger, normal-moveout spectrum, pickup velocity value are rescaned；And iteration again Renewal speed ratio and angle parameter value.As shown in figure 14, in this example, the low frequency component of residual normal moveout very little, because This, without the spectrum of scanning speed again and parameter of regulating the speed, without iteration renewal speed ratio and angle parameter again.

15, next analysis site is analyzed.

Repeatable 2 to 14 steps, complete the azimuthal anisotropy Parameter analysis of next analysis site.

The scheme that this specification said one or multiple embodiments provide gives a set of interactive iteration that is suitable for and analyzes Azimuthal anisotropy parameter, comprising: radius of equivalent circle speed, short long axis speed ratio and long axis angle.These three parameters have Determining physical significance, (comprehensive speed, intensity and crack including fracture development stratum are walked with the physical features in crack To) there is specific corresponding relationship.

Meanwhile in determining the azimuthal anisotropy parametric procedure, longitudinal wave can also be realized using above three parameter The high and low frequency of dynamic correction value separates, and then realizes from orientation isotropism velocity analysis to azimuthal anisotropy Parameter analysis Smooth transition has established solid theoretical basis for alternating direction iteration method azimuthal anisotropy Parameter analysis.

On this theoretical basis, the side for implementing the analysis of longitudinal wave azimuthal anisotropy parameter interactive iteration is further illustrated Method.Radius of equivalent circle velocity parameter values are determined under the premise of the orientation assumption of isotropy first；It is then possible to use interactive iteration The short long axis speed of method progressive updating is than parameter and major axis orientation angular dimensions.Situation can also be evened up according to channel set lineups, Decide whether to correct radius of equivalent circle speed parameter, and whether carries out short long axis speed than parameter and major axis orientation angular dimensions Next round iteration updates.Using the scheme of this specification above-described embodiment, the successive optimization and whole process matter of parameter may be implemented Amount control ensure that even if under Low SNR, can also obtain preferably analytical effect by the method for progressive alternate.

All the embodiments in this specification are described in a progressive manner, the same or similar between each embodiment Part may refer to each other, and the highlights of each of the examples are differences from other embodiments.It is specific to implement carefully Section is referred to the description of preceding aspect embodiments, does not do repeat one by one herein.

The embodiment above is described only for this specification specific embodiment, and other embodiments are wanted in appended right In the range of seeking book.In some cases, the movement recorded in detail in the claims or step can be according to different from embodiments In sequence execute and still may be implemented desired result.In addition, process depicted in the drawing not necessarily require by Desired result is just able to achieve according to particular order or the consecutive order implementation shown.In some embodiments, at multitask Reason and parallel processing are also possible or may be advantageous.

A kind of method for determining longitudinal wave azimuthal anisotropy parameter that this specification one or more embodiment provides, can be with The size, shape and side of velocity ellipse are described using three radius of equivalent circle speed, long axis angle and short long axis speed ratio parameters Position.And it can be by gradually scanning the Parameter Spectrums of above three anisotropic parameters, according to the energy extreme's distribution of Parameter Spectrum It determines corresponding parameter value, and then realizes the interactive iteration analysis of three azimuthal anisotropy parameters.It is each using this specification Embodiment can accurately determine azimuthal anisotropy parameter, so more efficiently analyze speed caused by vertical fracture with The rule of azimuthal variation

Based on the method for determination longitudinal wave azimuthal anisotropy parameter described above, this specification one or more embodiment A kind of device of determining longitudinal wave azimuthal anisotropy parameter is also provided.The device may include having used this specification implementation System, software (application), module, component, server of example the method etc. simultaneously combine the necessary device for implementing hardware.It is based on Same innovation thinking, the device in one or more embodiments that this specification embodiment provides is as described in the following examples. Since the implementation that device solves the problems, such as is similar to method, the implementation of the specific device of this specification embodiment can be joined See the implementation of preceding method, overlaps will not be repeated.Used below, term " unit " or " module " may be implemented pre- Determine the combination of the software and/or hardware of function.Although device described in following embodiment is preferably realized with software, The realization of the combination of hardware or software and hardware is also that may be contemplated.Specifically, Figure 16 indicates one that specification provides Kind determines the modular structure schematic diagram of the Installation practice of longitudinal wave azimuthal anisotropy parameter, and as shown in figure 16, described device can To include:

Channel set screening module 102 can be used for filtering out present analysis point from the prestack trace gather data in target work area Channel set data；

Speed determination module 104 can be used for carrying out orientation to the channel set data according to default scanning speed section Isotropism or azimuthal anisotropy velocity scanning obtain the orientation isotropism or azimuthal anisotropy of the present analysis point Normal-moveout spectrum determines the orientation isotropism or azimuthal anisotropy of present analysis point according to the energy extreme point of the normal-moveout spectrum Radius of equivalent circle speed；

Angle-determining module 106 can be used for the equivalent circular based on the orientation isotropism or azimuthal anisotropy half The present analysis point of the initial short long axis speed ratio or pickup of diameter speed and the preset present analysis point Current short long axis speed ratio, according to default scanning angle section interaction scanning angle spectrum, according to the energy extreme value of the angular spectrum Point determines the long axis angle of the present analysis point；

Speed can be used for the equivalent circular based on the orientation isotropism or azimuthal anisotropy than determining module 108 The long axis angle of radius speed and present analysis point interacts scanning speed comparison than section according to default scanning speed, The short long axis speed ratio of the present analysis point is determined according to the energy extreme point of the speed comparison.

It should be noted that device described above can also include other embodiment party according to the description of embodiment of the method Formula.Concrete implementation mode is referred to the description of related method embodiment, does not repeat one by one herein.

A kind of device for determining longitudinal wave azimuthal anisotropy parameter that this specification one or more embodiment provides, can be with The size, shape and side of velocity ellipse are described using three radius of equivalent circle speed, long axis angle and short long axis speed ratio parameters Position.And it can be by gradually scanning the Parameter Spectrums of above three anisotropic parameters, according to the energy extreme's distribution of Parameter Spectrum It determines corresponding parameter value, and then realizes the interactive iteration analysis of three azimuthal anisotropy parameters.It is each using this specification Embodiment can accurately determine azimuthal anisotropy parameter, so more efficiently analyze speed caused by vertical fracture with The rule of azimuthal variation.

Method or apparatus described in above-described embodiment that this specification provides can realize that business is patrolled by computer program It collects and records on a storage medium.The storage medium may include: the physical unit for storing information, usually will letter After breath digitlization, the media of the modes such as electricity, magnetic or optics is recycled to be stored.The storage medium may include: to utilize Electric energy mode stores the device of information such as, various memory, such as RAM, ROM；The device of information is stored in the way of magnetic energy such as, Hard disk, floppy disk, tape, core memory, magnetic bubble memory, USB flash disk；Using optical mode storage information device such as, CD or DVD.Certainly, there are also readable storage medium storing program for executing of other way, such as quantum memory, graphene memory etc..

The computer program of the storage medium storage can be read and executed by a computer, and realize this specification The effect of scheme described by embodiment.Therefore, this specification also provides a kind of equipment of determining longitudinal wave azimuthal anisotropy parameter, Memory including processor and storage processor executable instruction, described instruction may be implemented down when being executed by the processor The service logic step stated:

Firstly, filtering out the channel set data of present analysis point from the prestack trace gather data in target work area；

Then, it is determined that the isotropic radius of equivalent circle speed in orientation of the present analysis point, comprising: swept according to default It retouches speed interval and isotropism velocity scanning in orientation is carried out to the channel set data, the orientation for obtaining the present analysis point is each To same sex normal-moveout spectrum, the isotropic equivalent circular in orientation half of present analysis point is determined according to the energy extreme point of the normal-moveout spectrum Diameter speed；

Then, it is based on the isotropic radius of equivalent circle speed in the orientation, determines the short long axis of the present analysis point Speed ratio and long axis angle, comprising:

Based on the initial short long axis speed ratio of the preset present analysis point, swept according to the interaction of default scanning angle section Angular spectrum is retouched, the long axis angle of the present analysis point is determined according to the energy extreme point of the angular spectrum；And based on described The long axis angle of present analysis point, according to default scanning speed than section interaction scanning speed comparison, according to the speed comparison Energy extreme point determine the short long axis speed ratio of the present analysis point.

Figure 17 indicates the determination longitudinal wave azimuthal anisotropy parameter in one or more embodiment that this specification provides Equipment logical construction schematic diagram.The logical construction can both be realized with the mode of electronic component, can also be used tricks The mode of calculation machine software is realized；It should be noted that when the logic for realizing the equipment by the way of the computer software When structure, the element for forming the equipment can be implemented as forming the component of the computer software.

As shown in figure 17, a kind of equipment of determining longitudinal wave azimuthal anisotropy parameter can be relatively independent by 12 Element composition.Wherein, centrally located Cooperation controlling unit plays maincenter control action, has user instruction interface function.Top (speed spectrum analysis element, speed are than spectrum analysis element, angular spectrum analysis element, channel set Dynamically Announce member for 4 elements in portion Part) have graphical user interface function.3 composite components (normal-moveout spectrum computing element, speed comparison computing element and the angle of bottom Degree spectrum computing element) all built-in respective dedicated channel set moves correcting element and energy road computing element, speed is undertaken respectively Spectrum, the calculating task of speed comparison and angular spectrum.

The data management element in left side is responsible for the management work of the input of data, output and calculating parameter.The super road in left side Collection data sorting element is responsible for the sorting operation of channel set data.The channel set on right side moves correcting element and is responsible for channel set data Dynamic correction calculation.Three azimuthal anisotropy parameters (equivalent circular half is arranged in the azimuthal anisotropy parameter initialization element on right side Diameter speed, short long axis speed ratio and long axis angle) initial value.

The data management element can realize data input, output and operation ginseng under the control of Cooperation controlling element Several management functions.Its function may include:

Receive the data input instruction that Cooperation controlling element is sent；Longitudinal wave azimuthal anisotropy parameter is inputted from device external Prestack trace gather data needed for analysis；Give prestack trace gather data transmission to Cooperation controlling element.

Receive the data output instruction that Cooperation controlling element is sent；Receive the result data that Cooperation controlling element transmits；It will Result data is output to device external.Be output to device external data be analysis obtain azimuthal anisotropy parameter (including Radius of equivalent circle speed, short long axis speed ratio and long axis angle).

It receives the parameter that Cooperation controlling element is sent and saves instruction；Operational parameter (operational parameter packet is saved inside element Include present analysis point position, trace gather sorting parameter and scanning range parameter；Present analysis point position includes: the line of present analysis point Number and No. CMP of present analysis point；It includes adjacent lines number and adjacent C MP number that trace gather, which sorts parameter,；Scanning range parameter packet It includes: the scanning range of the scanning range of speed parameter, the scanning range of ratio parameter and angle parameter.Wherein, speed parameter Scanning range includes: starting velocity, terminates speed and speed number；The scanning range of ratio parameter includes: starting ratio, terminates Three parameters of ratio and ratio number；The scanning range of angle parameter includes start angle, termination point, angle number three ginsengs Number).

The data management element also can receive the parameter acquisition instruction that Cooperation controlling element is sent；It will be saved inside equipment Operational parameter feed back to Cooperation controlling element.

The channel set data sorting element can receive the channel set data sort instruction that Cooperation controlling element is sent, from Cooperation controlling element receives input data and required parameter, completes the sorting of channel set data.Its function may include:

Receive the channel set data sort instruction that Cooperation controlling element is sent；Present analysis point is received from Cooperation controlling element Position and channel set sort parameter；Prestack trace gather data are sorted, are sorted and base-level correction, present analysis point is formed Channel set data；Separation results are fed back into Cooperation controlling element.

The azimuthal anisotropy parameter initialization element can receive the azimuthal anisotropy that Cooperation controlling element is sent Parameter initialization instruction, completes azimuthal anisotropy parameter initialization.Its function may include:

Receive the azimuthal anisotropy parameter initialization instruction that Cooperation controlling element is sent；Setting radius of equivalent circle speed, The initial value of three short long axis speed ratio, long axis angle sweep parameters (physical significance of each parameter is as shown in Fig. 1).Its is medium Radius of circle speed parameter is set as 3600 meter per seconds to effect, short long axis speed is set as 1.0, major axis orientation angular dimensions than parameter and is set as 0 degree； Each parameter value after initialization is fed back into Cooperation controlling element.

The channel set, which moves correcting element, can not only be used for independent individual presence, also built-in in three kinds of spectrum computing elements (packet Include normal-moveout spectrum computing element, speed comparison computing element and angular spectrum computing element).

When the channel set, which moves correcting element, to be individually present, function may include: to receive Cooperation controlling element to send Dynamic correction instruction；Dynamic correction calculation data and parameter are received from Cooperation controlling element；Complete the dynamic correction calculation of channel set；It will Operation result feeds back to Cooperation controlling element.

When the channel set, which moves correcting element, to be built in inside other elements, function may include: to receive host's member The dynamic correction instruction that part is sent；Dynamic correction calculation data and parameter are received from host's element；It completes channel set and moves correction calculation.It will Operation result feeds back to host's element.

Orientation isotropism move correction calculation use travelling when calculation formula can with reference formula (3), and orientation respectively to Calculation formula can be with reference formula (5) when the travelling that the dynamic correction calculation of the opposite sex uses.When short long axis speed is equal to 1.0 than parameter When, when travelling used in the dynamic correction of formula (5) computer azimuth isotropism also can be used.

Energy road computing element can be built in three kinds of spectrum computing elements (including normal-moveout spectrum computing element, speed Comparison computing element and angular spectrum computing element).Its function may include: the energy road computations for receiving host's element；From Host's element receives the channel set data after dynamic correction；Summation and the square operation for completing channel set, obtain a stack power Road；Calculated result is fed back into host's element.The formula that energy road calculates can be with reference formula (4).

After having calculated all sample values in the energy road, just obtained a certain sweep parameter (velocity ellipse Radius of equivalent circle speed, short long axis speed ratio, long axis angle thrin) corresponding stack power road.

May include two built-in elements inside the normal-moveout spectrum computing element: channel set moves correcting element and energy road meter Calculate element.After normal-moveout spectrum computing element is connected to the normal-moveout spectrum computations that Cooperation controlling element is sent, connect from Cooperation controlling element Channel set data and normal-moveout spectrum after contracture choosing calculate required operational parameter, complete normal-moveout spectrum and calculate, and calculated result is anti- It feeds Cooperation controlling element.

The function of normal-moveout spectrum computing element may include: the normal-moveout spectrum computations for receiving Cooperation controlling element and issuing；From Cooperation controlling element receives channel set data；Scanning speed range parameter is received from Cooperation controlling element；According to velocity scanning model It encloses and (includes: starting velocity, terminates speed and speed number), determine scanning speed sequence；Call the dynamic correction of built-in channel set The function of element carries out dynamic correction to channel set respectively with each scanning speed；Built-in energy road computing element is recycled to calculate The corresponding stack power road of each speed；The corresponding stack power road of each scanning speed is arranged according to the sequence of speed from small to large again Column, obtain normal-moveout spectrum；Calculated result is fed back into Cooperation controlling element.

May include two built-in elements inside the speed comparison computing element: channel set moves correcting element and energy road Computing element.Speed comparison computing element can receive the speed comparison computations of Cooperation controlling element sending, complete speed Comparison calculates, and calculated result is fed back to Cooperation controlling element.

The function of speed comparison computing element may include: to receive the speed comparison calculating that Cooperation controlling element issues to refer to It enables；Channel set data are received from Cooperation controlling element；Speed ratio range parameter is received from Cooperation controlling element；It is swept according to speed ratio Range (include: starting ratio, terminate ratio and ratio number) is retouched, determines that scanning speed compares sequence；Call built-in channel set The function of dynamic correcting element, with each scanning speed than carrying out dynamic correction to channel set respectively；Utilize built-in energy road Computing Meta Part calculates each speed than corresponding stack power road；By each speed than corresponding stack power road according to speed ratio from small to large Sequence arranges, and obtains short long axis speed comparison；Calculated result is fed back into Cooperation controlling element.

May include two built-in elements inside the angular spectrum computing element: channel set moves correcting element and energy road meter Calculate element.Angular spectrum computing element can receive the instruction of Cooperation controlling element, completes angular spectrum and calculates, and calculated result is anti- It feeds Cooperation controlling element.

The function of angular spectrum computing element may include: the angular spectrum computations for receiving Cooperation controlling element and issuing；From Cooperation controlling element receives channel set data；(include: start angle, terminate from Cooperation controlling element receiving angle range parameter Angle and angle number)；According to angle scanning range, scanning angle sequence is determined；Built-in channel set is called to move correcting element Function, dynamic correction is carried out to channel set respectively with each scanning angle；Each scanning is calculated using built-in energy road computing element The corresponding stack power road of angle；The corresponding stack power road of each scanning angle is arranged according to the sequence of angle from small to large again Column, obtain angular spectrum；Calculated result is fed back into Cooperation controlling element.

The speed spectrum analysis element can complete the display and equivalent circular of normal-moveout spectrum under the control of Cooperation controlling element The pickup of radius speed parameter.Its function may include: to receive the instruction for the renewal speed spectrum display that Cooperation controlling element issues； New speed modal data is received from Cooperation controlling element；Show speed modal data；The interactive action of user is responded, equivalent circular is completed The pickup of radius speed parameter；By newest pickup as a result, real-time informing Cooperation controlling element.

The speed can complete the display of speed comparison and short under the control of Cooperation controlling element than spectrum analysis element Pickup of the long axis speed than parameter.Its function may include: to receive the renewal speed comparison that Cooperation controlling element issues to show Instruction；New speed ratio modal data is received from Cooperation controlling element；Show speed ratio modal data；The interactive action of user is responded, Complete pickup of the short long axis speed than parameter；And by newest pickup as a result, real-time informing Cooperation controlling element.

The angular spectrum analysis element can complete the display and major axis angle of angular spectrum under the assistance of Cooperation controlling element Spend the pickup of parameter.Its function may include: the instruction for receiving the update angular spectrum that Cooperation controlling element issues and showing；From cooperation Control element receives new angle modal data；Angles of display modal data；The interactive action of user is responded, long axis angle parameter is completed Pickup；And by newest pickup as a result, real-time informing Cooperation controlling element.

The dynamic that the channel set Dynamically Announce element under the control of Cooperation controlling element, can complete channel set is aobvious Show.Its function may include: the instruction for receiving the update channel set data that Cooperation controlling element issues and showing；Receive new dynamic school Channel set data after just；Show new channel set data.

The Cooperation controlling element is the maincenter control element of this equipment, and function is to coordinate and control the association of other elements Make and run, multiple relatively independent elements is made to form a co-operating entirety, it is common to realize azimuthal anisotropy parameter The function of interactive iteration analysis.The element has user instruction interface, can receive the instruction control of user, receive user and provide Processing parameter.Cooperation controlling element and other elements are connected with each other, and are played intermediation, are realized inter-component data, information With the access of instruction transmitting.

The function of Cooperation controlling element may include:

User instruction interface is provided, receive user instruction control (starting and termination including process, whether carry out it is next Whether wheel iteration continues to analyze next analysis site instruction)；Receive user provide processing parameter (including analysis site position, ginseng Number scanning range, azimuthal anisotropy initial parameter value, short long axis speed ratio parameter adjusted)；The finger provided according to user Order and processing parameter, establishing the process that azimuthal anisotropy parameter interactive iteration is analyzed, (process step is referring to method of the invention Part).

Data input instruction is sent to data management element；Receive the prestack trace gather data of data management cell feeds back.To Data management element sends data output instruction；Data to be output are sent to data management element；It is sent out to data management element Send parameter acquisition instruction；Receive the operational parameter of parameter management cell feeds back；Parameter, which is sent, to data management element saves instruction； Parameter to be saved is sent to data management element.

Channel set data sort instruction is sent to channel set sorting element；Prestack trace gather number is sent to channel set sorting element According to；Channel set data, which are sent, to channel set sorting element sorts parameter；Receive the separation results number of channel set sorting cell feeds back According to.

The initialization directive of orientation anisotropic parameters is sent to azimuthal anisotropy initialization element；Orientation is received respectively to different Property initialization cell feeds back azimuthal anisotropy initial parameter value.

Correcting element, which is moved, to channel set sends dynamic correction instruction；Correcting element, which is moved, to channel set sends dynamic correction calculation data With parameter (radius of equivalent circle speed, short long axis speed ratio and long axis angle that parameter includes: velocity ellipse)；It is dynamic to receive channel set The operation result of correcting element feedback.

Computations are composed to normal-moveout spectrum computing element transmission speed；Channel set data are sent to normal-moveout spectrum computing element；To Normal-moveout spectrum computing element sends scanning speed range parameter (include: starting velocity, terminate speed and speed number)；Receive speed Compose the normal-moveout spectrum calculated result of computing element feedback.

To speed comparison computing element transmission speed comparison computations；Channel set number is sent to speed comparison computing element According to；Scanning speed is sent to speed comparison computing element (to include: starting ratio, terminate ratio and ratio than range parameter Number)；Receive the speed comparison calculated result of speed comparison computing element feedback.

Angular spectrum computations are sent to angular spectrum computing element；Channel set data are sent to angular spectrum computing element；To Angle span parameter (including: start angle, termination point and angle number) is swept in the transmission of angular spectrum computing element；Receiving angle spectrum The angular spectrum calculated result of computing element feedback.

The instruction of renewal speed spectrum display is sent to speed spectrum analysis element；New speed is sent to speed spectrum analysis element Modal data；The speed parameter for receiving speed spectrum analysis cell feeds back picks up result.

The instruction that renewal speed comparison is shown is sent than spectrum analysis element to speed；It is sent newly to speed than spectrum analysis element Speed ratio modal data；The speed that speed is received than spectrum analysis cell feeds back picks up result than parameter.

The instruction for updating angular spectrum and showing is sent to angular spectrum analysis element；New angle is sent to angular spectrum analysis element Modal data；The angle parameter of receiving angle spectrum analysis cell feeds back picks up result.

The instruction for updating channel set data and showing is sent to channel set Dynamically Announce element；To channel set Dynamically Announce element Send channel set data to be shown.

Figure 18 indicates the graphic user interface schematic diagram of four visualization elements.Four operation windows in Figure 18 are all can Example depending on changing the graphic user interface of element, it is corresponding with four visualization elements at the top of equipment shown in Figure 17 respectively.From Three graphical windows of left-to-right be respectively speed spectrum analysis element, speed than spectrum analysis element and angular spectrum analysis element figure Shape user interface；The graphical window of the rightmost side is the graphic user interface of channel set Dynamically Announce element.Four windows are lateral simultaneously Column arrange, longitudinal direction is aligned by same time.This layout shows convenient for spectrum scanning and the comparison of dynamic correction result, also reduces friendship The difficulty of mutual pickup work.

In order to better illustrate the operability and operating mechanism of above equipment, this specification additionally provides above-mentioned one kind Determine the dynamic duty mechanism of the equipment of longitudinal wave azimuthal anisotropy parameter, as follows:

1, the position of first analysis site is set

The wire size of first analysis site can be set as 200 by Cooperation controlling element, be set as 100 for its No. CMP；And it will be first A analysis site (wire size 100 and No. CMP 100) is set as present analysis point.Cooperation controlling element is sent to data management element joins Number saves instruction and parameter value；After data management element is connected to parameter preservation instruction, parameter is received and saved, and store and arrive data Manage element.

2, setting trace gather sorts parameter

Trace gather can be set by Cooperation controlling element and sort parameter (including two classifications of span parameter and parameters sortnig). Span parameter includes two parameters of adjacent lines number (being set as 1) and adjacent C MP number (being set as 1)；Parameters sortnig includes starting big gun Inspection is away from (being set as 0 meter), termination geophone offset (being set as 4500 meters) and three parameters of geophone offset increment (being set as 500 meters).Cooperation controlling Element sends parameter to data management element and saves instruction and parameter value to be saved；Data management element is connected to parameter preservation and refers to After order, parameter to be saved is received, and store and arrive data management element.

3, channel set data are sorted

Cooperation controlling element sends the input instruction of prestack trace gather to data management element, and data management element is connected to the finger After order, according to the position (wire size 100 and No. CMP 200) of present analysis point and sorting parameter (adjacent lines number 1, adjacent C MP number Mesh 1, terminates 500 meters of 4500 meters of geophone offset and geophone offset increment at 0 meter of geophone offset of starting) input prestack trace gather data, and by data Cooperation controlling element is given with sorting parameter feedback.Cooperation controlling element sends channel set data point to channel set data sorting element Choosing instruction, data to be sorted and sorting parameter.

After channel set data sorting element receives the channel set data sort instruction of Cooperation controlling element transmission, from cooperation Control element receives prestack trace gather data and sorting parameter, is sorted, is sorted and base-level correction, and present analysis point is formed Channel set data, and feed back to Cooperation controlling element.The channel set data of present analysis point are sent to number by Cooperation controlling element According to management element, after data management element receives channel set data, it is spare to carry out inner buffer.

4, parameter scanning range is set

Parameter scanning range (the scanning range including speed parameter: starting velocity can be set by Cooperation controlling element 1800 meter per seconds, termination 5400 meter per second of speed, speed number 121；Scanning range of the speed than parameter: starting velocity than 0.7, Speed ratio 1.0 is terminated, speed is than number 31；The scanning range of angle parameter: start angle -90 is spent, and 90 degree of termination point, angle Number 91).Cooperation controlling element is sent to data management element saves scanning range instruction and parameter scanning value range.Data pipe After reason element is connected to the preservation scanning range instruction that Cooperation controlling element is sent, parameter to be saved is received from Cooperation controlling element It is spare to carry out inner buffer for value.

5, the initial value of azimuthal anisotropy parameter is set

Cooperation controlling element sends the initialization of orientation anisotropic parameters to azimuthal anisotropy parameter initialization element and refers to It enables, after data management element is connected to the azimuthal anisotropy parameter initialization instruction that Cooperation controlling element is sent, is arranged each to different The initial value (3600 meter per second of radius of equivalent circle speed, short long axis speed ratio 1.0,0 degree of long axis angle) of property parameter, by each to different The initial value of property parameter feeds back to Cooperation controlling element.Cooperation controlling element receives the initial value of azimuthal anisotropy parameter Afterwards, data management element is forwarded it to.After data management element receives the initial value of azimuthal anisotropy parameter, save as The current value of azimuthal anisotropy parameter.

6, first scan normal-moveout spectrum

Cooperation controlling element from data management element obtain channel set data, velocity scanning range (1800 meters of starting velocity/ Second, terminate 5400 meter per second of speed, speed number 121) and azimuthal anisotropy parameter current value (radius of equivalent circle speed 3600 meter per seconds, short long axis speed ratio 1.0, long axis angle 0), and it is sent to normal-moveout spectrum computing element.Normal-moveout spectrum computing element root According to the velocity scanning range received, a series of scanning speeds (1800 meter per seconds, 1830 meter per seconds, 1860 meter per seconds ... 5400 are given Meter per second), it calls built-in channel set to move correcting element respectively and dynamic correction calculation is carried out to the channel set data received, then Built-in energy road computing element is called, the corresponding energy road of each scanning speed is calculated.Again by the corresponding superposition of each scanning speed Energy road is arranged according to the sequence (1800 meter per seconds, 1830 meter per seconds, 1860 meter per second ..., 5400 meter per second) of speed from small to large, just Obtain the normal-moveout spectrum of present analysis point.

Then, the normal-moveout spectrum data feedback of present analysis point is given to Cooperation controlling element.Cooperation controlling element receives speed After spending modal data, speed spectrum analysis element is forwarded it to.After speed spectrum analysis element receives speed modal data, renewal speed The normal-moveout spectrum of spectrum analysis element shows figure.

7, pickup velocity parameter

Be referred to the energy extreme point distribution situation of speed spectrum analysis element speeds spectrogram shape, add new pickup point or Existing pickup point is adjusted, the longitudinally varying trend of pickup point is made to meet the regularity of distribution of normal-moveout spectrum energy extreme point.It picks up During speed parameter, the velocity amplitude of pickup is sent to Cooperation controlling element by speed spectrum analysis element in real time, then by cooperating The velocity amplitude notification data of pickup is managed element by control element, by data management element by current radius of equivalent circle velocity amplitude It is updated to the velocity amplitude of current pickup, and is cached.

Cooperation controlling element is by current azimuthal anisotropy parameter (radius of equivalent circle speed, short long axis speed ratio, long axis Angle) it is transmitted to the dynamic correcting element of channel set.Then, Cooperation controlling element obtains the super of present analysis point from data management element Trace gather data, and be transmitted to channel set and move correcting element.Channel set moves correcting element and receives the side that Cooperation controlling element is sent After position anisotropic parameters and channel set data, dynamic correction calculation carried out to channel set, and by the channel set data after dynamic correction Feed back to Cooperation controlling element.

Finally, Cooperation controlling element gives the channel set data forwarding after dynamic correction to channel set Dynamically Announce element, by surpassing Trace gather Dynamically Announce element real-time update shows the channel set data after dynamic correction, the speed ginseng picked up for user's real-time judge Whether number is reasonable.

Because the value of current anisotropic parameters γ is 1.0 (setting in step 5), normal-moveout spectrum calculates (step 6) and channel set Dynamic correction is all based on the orientation assumption of isotropy, and current speed parameter is only capable of substantially drawing the reflection line-ups of channel set It is flat, and high frequency components amount exists as residual normal moveout when travelling caused by crack.

8, default speed ratio parameter is adjusted

It can notify Cooperation controlling unit according to the situation of evening up of current channel set lineups, present speed is compared into parameter Value is set as 0.98.After Cooperation controlling element is connected to instruction, reinforms data management element and be adjusted to speed than parameter value 0.98, and by data management element cache.

9, scanning angle is composed

Cooperation controlling element from data management element obtain channel set data, velocity scanning range (start angle -90 degree, 90 degree of termination point, angle number 91) and azimuthal anisotropy parameter current value (step 7 determine radius of equivalent circle speed, Short long axis speed ratio 0.98, long axis angle 0), and it is sent to angular spectrum computing element.

Angular spectrum computing element gives a series of scanning angles (- 90 degree, -88 according to the angle scanning range received Degree, -86 90 degree of degree ...), it calls built-in channel set to move correcting element respectively and dynamic correction is carried out to the channel set data received Operation.Then, built-in energy road computing element is called, the corresponding energy road of each scanning angle is calculated.Again by each scanning angle Corresponding stack power road according to angle from small to large sequence (- 90, -88, -86 ... 90) arrange, just obtain present analysis point Angular spectrum.

Then, the angular spectrum data feedback of present analysis point is given to Cooperation controlling element.Cooperation controlling element receives angle After spending modal data, angular spectrum analysis element is forwarded it to.After angular spectrum analysis element receives angle modal data, angle is updated The display figure of spectrum analysis element angle spectrum.

10, pick-up angles parameter

Be referred to the energy extreme point distribution situation of angular spectrum analysis element angular spectrum figure, add new pickup point or Existing pickup point is adjusted, the longitudinally varying trend of pickup point is made to meet the regularity of distribution of angular spectrum energy extreme point.It picks up During angle parameter, the angle value of pickup is sent to Cooperation controlling element by angular spectrum analysis element in real time, then by cooperating The angle value of pickup is sent to data management element by control element, is updated current long axis angle value by data management element For the angle value of current pickup, and cached.

Cooperation controlling element is by current azimuthal anisotropy parameter (radius of equivalent circle speed, short long axis speed ratio, long axis Angle) it is transmitted to the dynamic correcting element of channel set.Then, Cooperation controlling element obtains the super of present analysis point from data management element Trace gather data, and be transmitted to channel set and move correcting element.Channel set moves correcting element and receives the side that Cooperation controlling element is sent After position anisotropic parameters and channel set data, dynamic correction calculation carried out to channel set, and by the channel set data after dynamic correction Feed back to Cooperation controlling element.

Finally, Cooperation controlling element gives the channel set data forwarding after dynamic correction to channel set Dynamically Announce element, by surpassing Trace gather Dynamically Announce element real-time update shows the channel set data after dynamic correction, the angle ginseng picked up for user's real-time judge Whether number is reasonable.

11, scanning speed comparison

Cooperation controlling element obtains channel set data, speed ratio scanning range (starting velocity ratio from data management element 0.7, terminate speed is than 1.0 and speed is than number 31) and azimuthal anisotropy parameter current value (step 7 determination equivalent circular The long axis angle that radius speed, short long axis speed ratio 0.98, step 10 determine), and speed is sent to than spectrum computing element.

Speed comparison computing element according to the speed ratio scanning range received, give a series of scanning speed ratios (0.7, 0.71,0.72 ... 1.0), call built-in channel set to move correcting element respectively and carry out dynamic correction to the channel set data received Then built-in energy road computing element is called in operation, calculate each scanning speed than corresponding energy road.Again by each scanning speed Than corresponding stack power road, according to speed ratio sequence from small to large, (0.7,0.71,0.72 ... 1.0) is arranged, and just obtains current The speed comparison of analysis site.

Then, the speed of present analysis point is fed back into Cooperation controlling element than modal data.Cooperation controlling element receives After speed is than modal data, speed is forwarded it to than spectrum analysis element.Speed receives speed ratio modal data than spectrum analysis element Afterwards, the display figure of renewal speed comparison analysis element speed comparison.

12, pickup velocity compares parameter

It is referred to energy extreme point distribution situation of the speed than spectrum analysis element speeds comparison figure, adds new pickup Point adjusts existing pickup point, the regularity of distribution for making the longitudinally varying trend of pickup point meet speed than spectrum energy extreme point. During pickup velocity is than parameter, the velocity ratio of pickup is sent to Cooperation controlling member than spectrum analysis element by speed in real time Part, then the velocity ratio notification data of pickup is managed into element by Cooperation controlling element, by data management element by current speed Degree ratio is updated to the velocity ratio of current pickup, and is cached.

Cooperation controlling element is by current azimuthal anisotropy parameter (radius of equivalent circle speed, short long axis speed ratio, long axis Angle) it is transmitted to the dynamic correcting element of channel set.Then, Cooperation controlling element obtains the super of present analysis point from data management element Trace gather data, and be transmitted to channel set and move correcting element.Channel set moves correcting element and receives the side that Cooperation controlling element is sent After position anisotropic parameters and channel set data, dynamic correction calculation carried out to channel set, and by the channel set data after dynamic correction Feed back to Cooperation controlling element.

Finally, Cooperation controlling element gives the channel set data forwarding after dynamic correction to channel set Dynamically Announce element, by surpassing Trace gather Dynamically Announce element real-time update shows the channel set data after dynamic correction, the speed ratio picked up for user's real-time judge Whether parameter is reasonable.

13, iteration renewal speed ratio and angle parameter

9 to 12 steps, progressive alternate renewal speed ratio and angle parameter are repeated, until result is satisfactory.Such as pass through iteration After update, speed comparison (as shown in figure 14) and angular spectrum (as shown in figure 11) energy extreme value point focusing are good and are reasonably distributed； Pickup point is located on energy extreme point；As shown in figure 15, the high fdrequency component of the residual normal moveout of channel set lineups drops to very It is small.

14, scanning speed composes (optional) again

It is first determined whether needing scanning speed spectrum again.Criterion is: the low frequency of the residual normal moveout of channel set Whether component is sufficiently small；As shown in figure 15, the low frequency component of the residual normal moveout of channel set lineups very little, therefore, Without the spectrum of scanning speed again and parameter of regulating the speed, without iteration renewal speed ratio and angle parameter again.

15, next velocity analysis point (optional) is analyzed

Repeatable 2 to 14 steps, complete the azimuthal anisotropy Parameter analysis of next analysis site.The present embodiment in order to avoid Repeated explanation decides not to the azimuthal anisotropy Parameter analysis for carrying out next analysis site.

The equipment of this specification longitudinal wave azimuthal anisotropy parameter interactive iteration analysis provided by the above embodiment, by 12 co-operating element compositions, centered on Cooperation controlling element, using maincenter formula logical construction.And establish cooperation Special information transfer channel between control element and other in non-built element, for instructing, parameter or data unidirectional or Transmitted in both directions.In addition to Cooperation controlling element, without direct communication between other in non-built elements.In order to improve work effect Rate, respectively a built-in channel set inside normal-moveout spectrum computing element, speed comparison computing element and angular spectrum computing element Dynamic correcting element and an energy road computing element, to mitigate the burden of Cooperation controlling element.The maincenter formula logical construction, is convenient for Each performs its own functions and interelement cooperates for element, is particularly suitable for realizing complicated multi-parameter interactive iteration analysis process, thus The efficiency of above-mentioned parameter analysis can be greatly improved.

It should be noted that the description of specific embodiment involved in above equipment can also include other embodiment party Formula.Concrete implementation mode is referred to the description of related method embodiment, does not repeat one by one herein.

A kind of equipment of determining longitudinal wave azimuthal anisotropy parameter, can use radius of equivalent circle described in above-described embodiment Three speed, long axis angle and short long axis speed ratio parameters describe the size, shape and orientation of velocity ellipse.And can by by The Parameter Spectrum of step scanning above three anisotropic parameters, determines corresponding parameter according to the energy extreme's distribution of Parameter Spectrum Value, and then realize the interactive iteration analysis of three azimuthal anisotropy parameters.It, can be accurate using each embodiment of this specification Ground determines azimuthal anisotropy parameter, and then more efficiently analyzes speed caused by vertical fracture with the rule of azimuthal variation Rule.

This specification also provides a kind of system of determining longitudinal wave azimuthal anisotropy parameter, and the system can be individual Longitudinal wave azimuthal anisotropy parametric analysis system, can also apply in other earthquake data processing systems.The system can Think individual server, also may include the one or more the methods or one or more implementations for having used this specification Server cluster, system (including distributed system), the software (application), practical operation device, logic gate of example device or equipment Circuit device, quantum computer etc. simultaneously combine the necessary terminal installation for implementing hardware.The determining longitudinal wave azimuthal anisotropy The system of parameter may include at least one processor and the memory for storing computer executable instructions, and the processor is held The mould of the step of method described in above-mentioned any one or multiple embodiments is realized when row described instruction or described device The element of block or the equipment.

It should be noted that system described above can also include according to the description of method, apparatus or apparatus embodiments Other embodiments, concrete implementation mode are referred to the description of related method embodiment, do not repeat one by one herein.

A kind of system of determining longitudinal wave azimuthal anisotropy parameter, can use radius of equivalent circle described in above-described embodiment Three speed, long axis angle and short long axis speed ratio parameters describe the size, shape and orientation of velocity ellipse.And can by by The Parameter Spectrum of step scanning above three anisotropic parameters, determines corresponding parameter according to the energy extreme's distribution of Parameter Spectrum Value, and then realize the interactive iteration analysis of three azimuthal anisotropy parameters.It, can be accurate using each embodiment of this specification Ground determines azimuthal anisotropy parameter, and then more efficiently analyzes speed caused by vertical fracture with the rule of azimuthal variation Rule.

It should be noted that this specification device, equipment or system described above is according to related method embodiment Description can also include other embodiments, and concrete implementation mode is referred to the description of embodiment of the method, does not make herein It repeats one by one.All the embodiments in this specification are described in a progressive manner, same and similar between each embodiment Part may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for hardware For+program class, storage medium+program embodiment, since it is substantially similar to the method embodiment, so describing simpler Single, the relevent part can refer to the partial explaination of embodiments of method.

Although referring to anisotropic parameters, Parameter Spectrum etc. in this specification embodiment content to obtain, definition, interaction, count Calculate, judgement etc. operations and data description, still, this specification embodiment be not limited to must be comply with standard data model/ Situation described in template or this specification embodiment.Certain professional standards use customized mode or embodiment to describe On practice processes embodiment modified slightly also may be implemented above-described embodiment it is identical, it is equivalent or it is close or deformation after can The implementation result of expectation.Using the implementation of the acquisitions such as these modifications or deformed data acquisition, storage, judgement, processing mode Example, still may belong within the scope of the optional embodiment of this specification.On the basis of this specification embodiment, pass through The variant of above-mentioned simple transformation, deformation or derivation and derivative parameter, method, apparatus, equipment or system, should all be included in this Within the protection scope of patent.

This specification specific embodiment is described in the embodiment above.Other embodiments are in the appended claims In the range of.In some cases, the movement recorded in detail in the claims or step can be according to different from embodiments Sequence executes and still may be implemented desired result.In addition, process depicted in the drawing is not necessarily required and is shown Particular order or consecutive order are just able to achieve desired result.In some embodiments, multitasking and parallel processing It is also possible or may be advantageous.

System, device, equipment, module or the unit that above-described embodiment illustrates, can be soft by computer chip or computer Part is realized, or is realized by the product with certain function.It is a kind of typically to realize that equipment is computer.Specifically, calculating Machine for example can be in personal computer, laptop computer, vehicle-mounted human-computer interaction device, tablet computer or these equipment Any equipment combination.

It for convenience of description, is that various modules or element are retouched respectively by function division when describing apparatus above or equipment It states.Certainly, in the one or more embodiments for implementing this specification, can each module or the function of element same or It is realized in multiple software and/or hardwares, the module of same function can also will be realized by the combination of multiple submodule or subcomponent Realize etc..Device or apparatus embodiments described above is only schematical, for example, the division of the module or element, Only a kind of logical function partition, there may be another division manner in actual implementation, such as multiple modules or element can be with In conjunction with or be desirably integrated into another system, or some features can be ignored or not executed.Another point, it is shown or discussed Data, parameter or transmission channel between equipment component is instructed to can be through some interfaces, between device, equipment or element It connects coupling or communication connection is realized, can be electrical property, mechanical or other forms.

It is also known in the art that other than realizing controller in a manner of pure computer readable program code, it is complete Entirely can by by method and step carry out programming in logic come so that controller with logic gate, switch, specific integrated circuit, programmable Logic controller realizes identical function with the form for being embedded in microcontroller etc..Therefore this controller is considered one kind Hardware component, and to the inside that can also be considered as hardware component for realizing the device or equipment of various functions that its inside includes Structure.Or even, can will be considered as realizing the device of various functions or equipment either the software of implementation method is real Body can be the internal structure of hardware component again.

The present invention referring to according to the method for the embodiment of the present invention, device, equipment, system and computer program product Flowchart and/or the block diagram describes.It should be understood that can be realized by computer program instructions every in flowchart and/or the block diagram The combination of process and/or box in one process and/or box and flowchart and/or the block diagram.It can be by the computer journey Sequence instructs the processing for being supplied to general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices Device, so that the instruction that the computer or the processor of other programmable data processing devices execute may be implemented in flow chart The service logic of function is specified in one box or multiple boxes in one process or multiple processes and/or block diagrams.

These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates, Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or The function of being specified in multiple boxes.

These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one The step of function of being specified in a box or multiple boxes.

In a typical configuration, the system of a set of determining longitudinal wave azimuthal anisotropy parameter may include one or more A processor (CPU), input/output interface, network interface and memory, and for realizing part described in this specification or entirely One or more sets computer softwares of portion's method.The method, apparatus and equipment that this specification one or more embodiment provides can To be embodied as the function, structure and component of the computer software.

It should also be noted that, the terms "include", "comprise" or its any other variant are intended to nonexcludability It include so that the process, method, commodity or the equipment that include a series of elements not only include those elements, but also to wrap Include other elements that are not explicitly listed, or further include for this process, method, commodity or equipment intrinsic want Element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including described want There is also other identical elements in the process, method or equipment of element.

It will be understood by those skilled in the art that this specification one or more embodiment can provide as method, system or calculating Machine program product.Therefore, this specification one or more embodiment can be used complete hardware embodiment, complete software embodiment or The form of embodiment combining software and hardware aspects.Moreover, this specification one or more embodiment can be used at one or It is multiple wherein include computer usable program code computer-usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) on the form of computer program product implemented.

This specification one or more embodiment can describe in the general context of computer executable instructions, such as journey Sequence module.Generally, program module include routines performing specific tasks or implementing specific abstract data types, programs, objects, Component, data structure etc..This specification one or more embodiment can also be practiced in a distributed computing environment, at these In distributed computing environment, by executing task by the connected remote processing devices of communication network.In distributed computing In environment, program module can be located in the local and remote computer storage media including storage equipment.

All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method Part explanation.In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", The description of " specific example " or " some examples " etc. means specific features described in conjunction with this embodiment or example, structure, material Or feature is contained at least one embodiment or example of this specification.In the present specification, to the signal of above-mentioned term Property statement must not necessarily be directed to identical embodiment or example.Moreover, specific features, structure, material or the spy of description Point may be combined in any suitable manner in any one or more of the embodiments or examples.In addition, without conflicting with each other, Those skilled in the art can be by different embodiments or examples described in this specification and different embodiments or examples Feature is combined.

The foregoing is merely the embodiments of this specification, are not limited to this specification.For art technology For personnel, this specification can have various modifications and variations.It is all made any within the spirit and principle of this specification Modification, equivalent replacement, improvement etc., should be included within the scope of the claims of this specification.

Claims (12)

1. a kind of method of determining longitudinal wave azimuthal anisotropy parameter, which is characterized in that the azimuthal anisotropy parameter includes One equivalent Circle Parameters and two elliptic parameters, wherein

The equivalent Circle Parameters include the radius of equivalent circle speed of velocity ellipse, and the radius of equivalent circle speed is for determining dynamic school The amplitude of the low frequency component of positive quantity；

The elliptic parameter includes the short long axis speed ratio of velocity ellipse, the long axis angle of velocity ellipse；Wherein, short long axle speed Than the amplitude for determining the high fdrequency component of dynamic correction value, long axis angle is used to determine the phase of the high fdrequency component of dynamic correction value；

The longitudinal wave azimuthal anisotropy parameter is determined in the method for interaction analysis, comprising:

The channel set data of present analysis point are filtered out from the prestack trace gather data in target work area；

Determine the isotropic radius of equivalent circle speed in orientation of the present analysis point, comprising: according to default scanning speed area Between isotropism velocity scanning in orientation is carried out to the channel set data, obtain the orientation isotropism speed of the present analysis point Degree spectrum, the isotropic radius of equivalent circle speed in orientation of present analysis point is determined according to the energy extreme point of the normal-moveout spectrum；

Based on the isotropic radius of equivalent circle speed in the orientation, the short long axis speed ratio and length of the present analysis point are determined Shaft angle degree, comprising: the initial short long axis speed ratio based on the preset present analysis point is handed over according to default scanning angle section Mutual scanning angle spectrum, the long axis angle of the present analysis point is determined according to the energy extreme point of the angular spectrum；And it is based on The long axis angle of the present analysis point, according to default scanning speed than section interaction scanning speed comparison, according to the speed The energy extreme point of comparison determines the short long axis speed ratio of the present analysis point.

2. the method for determining longitudinal wave azimuthal anisotropy parameter according to claim 1, which is characterized in that the determination is short Long axis speed ratio and long axis angle, comprising:

Based on the isotropic radius of equivalent circle speed in the orientation, described current point is determined using the method that interactive iteration is analyzed Analyse the short long axis speed ratio and long axis angle of point, wherein the method for interactive iteration analysis includes:

Based on the initial short long axis speed ratio of the preset present analysis point, according to default scanning angle section interaction scan angle Degree spectrum, the long axis angle of the present analysis point is determined according to the energy extreme point of the angular spectrum, and, execute following iteration Step:

Based on the long axis angle for the present analysis point that last iteration determines, scanned according to default scanning speed than section interaction Speed comparison determines the short long axis speed ratio of the present analysis point according to the energy extreme point of the speed comparison；

Based on the short long axis speed ratio for the present analysis point that current iteration determines, swept according to the interaction of default scanning angle section Angular spectrum is retouched, the long axis angle of the present analysis point is determined according to the energy extreme point of the angular spectrum；

Until meeting default iteration requirement, the short long axis speed for the present analysis point that output last time iteration obtains is when Long axis angle.

3. the method for determining longitudinal wave azimuthal anisotropy parameter according to claim 1 or 2, which is characterized in that the side Method further include:

Determine the radius of equivalent circle speed of the azimuthal anisotropy of the present analysis point, and, based on the orientation respectively to different The radius of equivalent circle speed of property determines the short long axis speed ratio and long axis angle of the present analysis point, wherein

The radius of equivalent circle speed of the azimuthal anisotropy of the determination present analysis point, comprising:

According to default scanning speed section, the short long axis speed ratio and long axis angle of present analysis point, to the channel set Data carry out azimuthal anisotropy velocity scanning, the azimuthal anisotropy normal-moveout spectrum of the present analysis point are obtained, according to described The energy extreme point of azimuthal anisotropy normal-moveout spectrum determines the radius of equivalent circle speed of the azimuthal anisotropy of present analysis point.

4. the method for determining longitudinal wave azimuthal anisotropy parameter according to claim 3, which is characterized in that the determining institute It states the isotropic radius of equivalent circle speed in orientation of present analysis point or determines the orientation of the present analysis point respectively to different The radius of equivalent circle speed of property, comprising:

According to the low frequency component amplitude size of the channel set residual normal moveout of the present analysis point, the present analysis is determined The radius of equivalent circle of the azimuthal anisotropy of the isotropic radius of equivalent circle speed in orientation or present analysis point of point Speed.

5. the method for determining longitudinal wave azimuthal anisotropy parameter according to claim 1 or 2, which is characterized in that described true The short long axis speed ratio of the fixed present analysis point, comprising:

The present analysis is determined according to the high fdrequency component amplitude size of the channel set residual normal moveout of the present analysis point The short long axis speed ratio of point.

6. the method for determining longitudinal wave azimuthal anisotropy parameter according to claim 1 or 2, which is characterized in that described pre- If the present analysis point initial short long axis speed ratio using following manner determine:

According to the high fdrequency component amplitude size of the channel set residual normal moveout of the present analysis point, the present analysis is set The initial short long axis speed ratio of point.

7. the method for determining longitudinal wave azimuthal anisotropy parameter according to claim 1 or 2, which is characterized in that described true The long axis angle of the fixed present analysis point, comprising:

Present analysis point is determined according to the high fdrequency component phase size of the channel set residual normal moveout of the present analysis point Angle value.

8. the method for determining longitudinal wave azimuthal anisotropy parameter according to claim 1, which is characterized in that described to filter out The channel set data of present analysis point, comprising:

Sort choice of parameters according to preset trace gather and go out the channel set data of present analysis point, the trace gather sorting parameter include across Spend parameter and parameters sortnig, wherein

The span parameter includes: adjacent lines number and adjacent C MP number；

The parameters sortnig includes: starting geophone offset, terminates geophone offset and geophone offset increment.

9. the method for determining longitudinal wave azimuthal anisotropy parameter according to claim 1, which is characterized in that described to filter out The channel set data of present analysis point, comprising:

The initial channel set data that parameter sub-elects present analysis point are sorted according to preset trace gather, the trace gather sorts parameter packet Include span parameter and parameters sortnig, wherein the span parameter includes adjacent lines number and adjacent C MP number, the sequence ginseng Number includes starting geophone offset, terminates geophone offset and geophone offset increment；

The initial channel set data sub-elected are grouped according to the parameters sortnig, the trace gather data in each group are based on side The incremental sequence of parallactic angle is ranked up；

Each seismic channel data in the initial channel set data is corrected on the datum level where present analysis point, is worked as The channel set data of preceding analysis site.

10. a kind of device of determining longitudinal wave azimuthal anisotropy parameter, which is characterized in that described device includes:

Channel set screening module, for filtering out the channel set number of present analysis point from the prestack trace gather data in target work area According to；

Speed determination module, for carrying out orientation isotropism or side to the channel set data according to default scanning speed section Position anisotropic velocity scanning, obtains the orientation isotropism or azimuthal anisotropy normal-moveout spectrum of the present analysis point, according to The energy extreme point of the normal-moveout spectrum determines the orientation isotropism of present analysis point or the radius of equivalent circle of azimuthal anisotropy Speed；

Angle-determining module, the equivalent circular half for orientation isotropism or azimuthal anisotropy based on the present analysis point The initial short long axis speed ratio or previous round iterative analysis of diameter speed and the preset present analysis point are obtained The short long axis speed ratio of present analysis point, according to default scanning angle section interaction scanning angle spectrum, according to the angular spectrum Energy extreme point determines the long axis angle of the present analysis point；

Speed ratio determining module, for based on the orientation isotropism or azimuthal anisotropy radius of equivalent circle speed, with And the long axis angle of the present analysis point, according to default scanning speed than section interaction scanning speed comparison, according to the speed The energy extreme point of degree comparison determines the short long axis speed ratio of the present analysis point.

11. a kind of equipment of determining longitudinal wave azimuthal anisotropy parameter, which is characterized in that the equipment includes: Cooperation controlling member Part, data management element, channel set data sorting element, channel set move correcting element, azimuthal anisotropy parameter initialization member Part, energy road computing element, normal-moveout spectrum computing element, speed comparison computing element, angular spectrum computing element, speed spectrum analysis member Part, speed are than spectrum analysis element, angular spectrum analysis element and channel set Dynamically Announce element；Wherein,

The Cooperation controlling element is used to coordinate and control the cooperation and operation of whole equipment；The Cooperation controlling element has use Family instruction interface, for receiving the instruction control of user and receiving the processing parameter of user's offer；

The data management element realizes the management of data input, output and operational parameter under the control of Cooperation controlling element Function；

The channel set data sorting element is for receiving the channel set data sort instruction that Cooperation controlling element is sent, from cooperation Control element receives input data and required parameter, realizes the sorting of channel set data；

The channel set, which moves correcting element, to be existed as independent component or exists as built-in element；The dynamic correction of the channel set The function of element includes: the dynamic correction instruction for receiving other elements and sending, and receives dynamic correction calculation data and ginseng from other elements Number；Complete the dynamic correction calculation of channel set；Operation result is fed back into other elements, wherein

In the presence of the channel set moves correcting element as independent component, function includes: to receive Cooperation controlling element to send Dynamic correction instruction receives dynamic correction calculation data and parameter from Cooperation controlling element, completes the dynamic correction calculation of channel set, will transport It calculates result and feeds back to Cooperation controlling element；

The channel set moves correcting element when being built in inside other elements, function include: receive host's element send it is dynamic Correction instruction receives dynamic correction calculation data and parameter from host's element, completes channel set and moves correction calculation, and operation result is anti- It feeds host's element；

The azimuthal anisotropy parameter initialization element, the azimuthal anisotropy parameter sent for receiving Cooperation controlling element Azimuthal anisotropy parameter initialization is completed in initialization directive；

Energy road computing element is built-in element；Its function includes: the energy road computations for receiving host's element, from place Major component receives the channel set data after dynamic correction, completes summation and the square operation of channel set, obtains a stack power road, Calculated result is fed back into host's element；

Include two built-in elements inside the normal-moveout spectrum computing element: channel set moves correcting element and energy road computing element； The normal-moveout spectrum computing element receives the normal-moveout spectrum computations that Cooperation controlling element is sent, and completes normal-moveout spectrum and calculates, and will meter It calculates result and feeds back to Cooperation controlling element；

Include two built-in elements inside the speed comparison computing element: channel set moves correcting element and energy road Computing Meta Part；The speed comparison computing element receives the speed comparison computations that Cooperation controlling element issues, and completes speed comparison meter It calculates, and calculated result is fed back into Cooperation controlling element；

The angular spectrum computing element includes two built-in elements: channel set moves correcting element and energy road computing element；It is described Angular spectrum computing element receives the instruction of Cooperation controlling element, completes angular spectrum and calculates, and calculated result is fed back to cooperation control Element processed；

The speed spectrum analysis element completes the display and radius of equivalent circle speed of normal-moveout spectrum under the control of Cooperation controlling element The pickup of parameter；

The speed under the control of Cooperation controlling element, completes the display and short long axle speed of speed comparison than spectrum analysis element Than the pickup of parameter；

The angular spectrum analysis element completes the display of angular spectrum and picking up for angle parameter under the control of Cooperation controlling element It takes；

The channel set Dynamically Announce element realizes the Dynamically Announce of channel set under the control of Cooperation controlling element.

12. a kind of system of determining longitudinal wave azimuthal anisotropy parameter, which is characterized in that including at least one processor and deposit The memory of computer executable instructions is stored up, the processor realizes any one of claim 1-9 when executing described instruction The step of the method.