CN104570098A

CN104570098A - Method and device for acquiring common-azimuth-angle seismic data

Info

Publication number: CN104570098A
Application number: CN201510044378.6A
Authority: CN
Inventors: 熊登; 赵波; 程元方; 王仕俭; 尹天奎
Original assignee: NATIONAL ENGINEERING RESEARCH CENTER OF OIL AND GAS EXPLORATION SOFTWARE Co Ltd; BGP Inc
Current assignee: NATIONAL ENGINEERING RESEARCH CENTER OF OIL AND GAS EXPLORATION SOFTWARE Co Ltd; BGP Inc
Priority date: 2015-01-28
Filing date: 2015-01-28
Publication date: 2015-04-29
Anticipated expiration: 2035-01-28
Also published as: CN104570098B

Abstract

The invention discloses a method and a device for acquiring common-azimuth-angle seismic data. The method includes: building a coordinate system which is used for mapping the seismic data; building statistical points, which are arranged on the basis of a preset rule and used for storing seismic channel data mapped to the statistical points and mapped to positions near the statistical points, in the coordinate system; determining a to-be-extracted azimuth angle, and building the index of the statistical points located in the to-be-extracted azimuth angle; mapping the seismic data to the coordinate system, and by the statistical points, storing the seismic channel data mapped to the statistical points and mapped to positions near the statistical points; traversing the index to obtain the seismic channel data stored by the statistical points located in the to-be-extracted azimuth angle. By the method and the device for acquiring the common-azimuth-angle seismic data, fast common-azimuth-angle seismic data extraction can be performed during wide-azimuth exploration seismic data processing.

Description

A kind of acquisition methods of position angle geological data altogether and device thereof

Technical field

The application relates to field of seismic exploration, particularly relates to a kind of acquisition methods and device thereof of position angle geological data altogether.

Background technology

The meticulous depiction of underground structure and lithology is needed to the even multi-fold of better illumination and more high reps, therefore wide-azimuth exploration increasing prevalence.Wide-azimuth exploration increases reception channel by several times, especially rolls up in big gun and receives line arrangement, with obtain meticulousr bin cover and in bin evenly azimuthal distribution, to lifting seism processing effect.But wide-azimuth has highlighted the anisotropic tremendous influence of different azimuth, must take into full account that this type of affects when processing.So when processing wide-azimuth seismic data, existing a lot for the treatment of step needs point position angle to process, and adjacent seismic trace is divided in a certain size position angle, as common orientation angle gathers, then sends into treatment scheme.In same position angle, it has been generally acknowledged that the parameters such as such as seismic velocity are consistent, again result is combined after processing each position angle data respectively, just can reduce the impact that the seismic velocity etc. that causes due to anisotropy changes, and then significantly improve the data processed results, so extracting position angle geological data is altogether the routine operation in wide-azimuth seism processing.

The acquisition methods of tradition altogether position angle geological data is divided into three steps usually: 1, utilize geological data coordinate, with trigonometric function directly calculate each road geological data position angle and stored in data trace header; By the position angle size of specified design when 2, using, all seismic traces are screened and classifies as common orientation angle gathers; 3, sort according to nominal key in each angle gathers inside, common orientation, and key word can be offset distance, wire size or Taoist monastic name.

The acquisition methods that above-mentioned tradition is total to position angle geological data uses fairly simple, when data volume is little, uses traditional acquisition methods conveniently can calculate extraction.But when wide-azimuth is explored, the total number of channels of geological data sharply increases, after sometimes increasing to current hundred million road scales, classic method sharply increases with the time spending of trigonometric function computer azimuth angle in step sl.If remerge the time spending of step 2 and step 3, then azimuthal geological data sorting extraction time can reach a few hours or a couple of days altogether, take the seismic data processing time in large quantities, affect the treatment effeciency of seismic data, and then pole is unfavorable for efficiently processing at earthquake-capturing scene.

Summary of the invention

In view of the deficiencies in the prior art, the application provides a kind of acquisition methods and device thereof of position angle geological data altogether, can carry out the extraction work of common position angle geological data fast when processing the geological data of wide-azimuth exploration.

The acquisition methods of a kind of position angle geological data altogether that the application provides, comprising:

Set up coordinate system, described coordinate system is for mapping geological data;

At the statistics point of the built-in pre-defined rule arrangement that is based on of described coordinate system; Described statistics point is mapped on this statistics point and the seismic channel data be mapped near this statistics point for storing; Described pre-defined rule comprises the coordinate range that the described coordinate range adding up point is more than or equal to geological data, and the quantity of described statistics point is less than the quantity of the seismic channel data of geological data;

Determine position angle to be extracted and set up the index being positioned at the statistics point at position angle to be extracted;

By seismic data mapping to described coordinate system, described statistics point will be mapped on this statistics point and be mapped in the seismic channel data storage near this statistics point;

Travel through described index and obtain the seismic channel data that the statistics point that is positioned at position angle to be extracted stores.

Preferably, describedly set up the coordinate range that coordinate system is included according to geological data on described coordinate system and mark off mapping area;

Described statistics point of arranging at the built-in pre-defined rule that is based on of described coordinate system comprises: at the statistics point of the built-in pre-defined rule arrangement that is based on of described mapping area that described coordinate system divides.

Preferably, described statistics point comprises integer array and integer variable, and described integer array stores the Taoist monastic name of seismic channel data for storing described statistics point institute, described integer variable for store described add up put store the quantity of seismic channel data.

Preferably, the number of the described integer array storage Taoist monastic name of each described statistics point is less than or equal to 5.

Preferably, describedly on described coordinate system, mark off according to the coordinate range of geological data the coordinate range that mapping area is included according to geological data on described coordinate system mark off mapping area by mesh lines;

Described pre-defined rule comprises described statistics point and is positioned at grid node locations.

Preferably, described acquisition methods also comprises the track pitch and line-spacing that obtain geological data;

Described mesh lines comprises many parallels parallel to each other and many warps vertical with described parallel, the transverse axis of described coordinate system and the longitudinal axis are two mutually orthogonal described mesh lines, distance often between adjacent two described parallels is equal, and the distance often between adjacent two described warps is equal.

Preferably, describedly determine position angle to be extracted and set up the index being positioned at the statistics point at position angle to be extracted to comprise:

Region, position angle to be extracted is marked off in described coordinate system;

Set up the absorbing boundary equation in region, described position angle to be extracted;

Judge whether described statistics point is positioned at region, described position angle to be extracted according to described absorbing boundary equation;

Set up the index of the described statistics point being positioned at position angle to be extracted.

Calculate the position angle of described statistics point;

Position angle according to described statistics point judges whether this statistics point is positioned at region, described position angle to be extracted;

Preferably, described index is the coordinate of described statistics point, or described index is the order label of the distant relationships according to described statistics point range coordinate initial point.

Preferably, described seismic data mapping to be comprised to described coordinate system:

Taoist monastic name and the quantity of the seismic channel data stored by described statistics point are removed;

Each seismic channel data of described geological data is mapped in described mapping area according to its mapping point.

Preferably, described geological data comprises multiple roads collection, described geological data carry out described seismic data mapping to described coordinate system and described traversal described index obtained the seismic channel data that the statistics point that is positioned at position angle to be extracted stores time, described multiple roads collection one by one road collection carries out that described extremely described coordinate system and the described index of described traversal obtain the seismic channel data that the statistics point that is positioned at position angle to be extracted stores by seismic data mapping.

Preferably, described geological data to comprise in common big gun collection, CMP collection, common receiving point gather the combination of one or more; The receiving station coordinate that the mapping point of the seismic channel data of described altogether big gun collection is this seismic channel data deduct shot point coordinate draw the multiple of coordinate figure; The receiving station coordinate that the mapping point of the seismic channel data of described CMP collection is this seismic trace deduct center point coordinate draw the multiple of coordinate figure; The shot point coordinate that the mapping point of the seismic trace of described common receiving point gather is this seismic trace deduct acceptance point coordinate draw the multiple of coordinate figure.

Preferably, the described index of described traversal obtains seismic channel data that the statistics point that is positioned at position angle to be extracted stores and comprises the described index of traversal and collect seismic channel data that the statistics point that is positioned at position angle to be extracted stores and be stored as common orientation angle gathers.

The application is a kind of acquisition device being total to position angle geological data also, comprising:

First sets up module, and it is for setting up coordinate system; Described coordinate system is for mapping geological data;

Second sets up module, and it is for the statistics point in the built-in pre-defined rule arrangement that is based on of described coordinate system; Described statistics point is mapped on this statistics point and the seismic channel data be mapped near this statistics point for storing; The quantity of coordinate range and described statistics point that the coordinate range that described pre-defined rule comprises described statistics point is more than or equal to geological data is less than the quantity of the seismic channel data of geological data;

3rd sets up module, and it is for determining position angle to be extracted and setting up the index being positioned at the statistics point at position angle to be extracted;

Data mapping module, its for by seismic data mapping to described coordinate system, described statistics point will be mapped on this statistics point and the seismic channel data be mapped near this statistics point stores;

Data acquisition module, it obtains for traveling through described index the seismic channel data that the statistics point that is positioned at position angle to be extracted stores.

Preferably, described first sets up module comprises Region dividing module, and described Region dividing module marks off mapping area for the coordinate range according to geological data on described coordinate system;

Described second sets up module comprises the second son and sets up module, and the second son sets up the statistics point of described mapping area built-in be based on pre-defined rule arrangement of module for dividing at described coordinate system.

Preferably, described Region dividing module comprises Region dividing submodule, and the coordinate range that described Region dividing submodule is included according to geological data on described coordinate system marks off mapping area by mesh lines;

Described second sets up in module, and described pre-defined rule comprises described statistics point and is positioned at grid node locations.

Preferably, described acquisition device also comprises acquisition module, and described acquisition module is for obtaining track pitch and the line-spacing of geological data;

In described Region dividing submodule, described mesh lines comprises many parallels parallel to each other and many warps vertical with described parallel, the transverse axis of described coordinate system and the longitudinal axis are two mutually orthogonal described mesh lines, distance often between adjacent two described parallels is equal, and the distance often between adjacent two described warps is equal.

Preferably, the described 3rd set up module and comprise:

Position angle divides A module, and it for marking off region, position angle to be extracted in described coordinate system;

Establishing equation module, it is for setting up the absorbing boundary equation in region, described position angle to be extracted;

Judge A module, it is for judging according to described absorbing boundary equation whether described statistics point is positioned at region, described position angle to be extracted;

A module set up in index, and it is for setting up the index of the described statistics point being positioned at position angle to be extracted.

Preferably, the described 3rd set up module and comprise:

Position angle computing module, it is for calculating the position angle of described statistics point;

Position angle divides B module, and it for marking off region, position angle to be extracted in described coordinate system;

Judge B module, according to the position angle of described statistics point, it is for judging whether this statistics point is positioned at region, described position angle to be extracted;

B module set up in index, and it is for setting up the index of the described statistics point being positioned at position angle to be extracted.

Preferably, described data mapping module comprises:

Initialization module, it is for adding up a Taoist monastic name of the seismic channel data stored and quantity removing by described;

Mapping block, it is for mapping to described mapping area by each seismic channel data of described geological data according to its mapping point.

Preferably, described geological data comprises multiple roads collection, when described geological data is by described data mapping module and described data acquisition module process, described multiple roads collection one by one road collection by described data mapping module and described data acquisition module process.

Preferably, described data acquisition module comprises collection memory module, and described collection memory module is collected seismic channel data that the statistics point that is positioned at position angle to be extracted stores for traveling through described index and is stored as common orientation angle gathers.

In sum, can find out that the technical scheme that the application provides is passed through on coordinate system, to set up data-storable statistics point, and then make the seismic channel data in geological data according on its virtual borderlines to described coordinate system, then pass through set up index and the seismic channel data stored in statistics point is carried out collection storage, so, the acquisition methods of present embodiment is used only to need region, position angle to be extracted to divide out on described coordinate system, and then can evaluation work be greatly reduced, make the extraction work can carrying out common position angle geological data when processing the geological data of wide-azimuth exploration fast.

Accompanying drawing explanation

In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those skilled in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.

Fig. 1 is the schematic flow sheet of the acquisition methods of the common position angle geological data that a kind of embodiment of the application provides;

Fig. 2 be the application's embodiment the step S3 schematic flow sheet of the acquisition methods of common position angle geological data is provided;

Fig. 3 be another embodiment of the application the step S3 schematic flow sheet of the acquisition methods of common position angle geological data is provided.

Fig. 4 be the application's embodiment the step S4 schematic flow sheet of the acquisition methods of common position angle geological data is provided;

Fig. 5 be the application's embodiment provide the acquisition methods of common position angle geological data to set up coordinate system schematic diagram;

Fig. 6 is the module composition schematic diagram of the acquisition device of the common position angle geological data that a kind of embodiment of the application provides;

Fig. 7 be the application's embodiment provide the 3rd of the acquisition device of common position angle geological data the set up module module composition schematic diagram;

Fig. 8 be another embodiment of the application provide the 3rd of the acquisition device of common position angle geological data the set up module module composition schematic diagram;

Fig. 9 be the application's embodiment provide the data mapping module of the acquisition device of common position angle geological data module composition schematic diagram.

Embodiment

Technical scheme in the application is understood better in order to make those skilled in the art person, below in conjunction with the accompanying drawing in the embodiment of the present application, technical scheme in the embodiment of the present application is clearly and completely described, obviously, described embodiment is only some embodiments of the present application, instead of whole embodiments.Based on the embodiment in the application, the every other embodiment that those of ordinary skill in the art obtain under the prerequisite not making creative work, all should belong to the scope of protection of the invention.

Please refer to Fig. 1, a kind of embodiment of the application provides a kind of acquisition methods of position angle geological data altogether, comprises the following steps:

S1, set up coordinate system, described coordinate system is for mapping geological data.

Described geological data includes a large amount of seismic channel data, especially, when carrying out wide-azimuth exploration, its seismic channel data more, by setting up described coordinate system, described coordinate system can provide at most 360 degree of position angles, and then farthest can hold geological data.

Described coordinate system can be orthogonal coordinate system, and it comprises x-axis (also can be called transverse axis) and y-axis (also can be called the longitudinal axis).Described geological data comprises multiple roads collection usually, such as comprises multiple big gun collection altogether, multiple CMP collection, multiple common receiving point gather etc.The geophone station normally arranged during image data in seismic exploration carries out grid arrangement according to vertical survey line and cross-track, so the seismic trace that each geophone station gathers has its corresponding coordinate.There is its corresponding track pitch and line-spacing in each seismic trace of each road collection.So the transverse axis of described coordinate system (x-axis) represents the described vertical survey line concentrated, the longitudinal axis (y-axis) of described coordinate system represents the described cross-track concentrated, and then makes need not can directly mapping in converted coordinate orientation of described geological data.

S2, the statistics point of arranging at the built-in pre-defined rule that is based on of described coordinate system; Described statistics point is mapped on this statistics point and the seismic channel data be mapped near this statistics point for storing; Described pre-defined rule comprises the coordinate range that the described coordinate range adding up point is more than or equal to geological data, and the quantity of described statistics point is less than the quantity of the seismic channel data of geological data.

Described statistics point is owing to being positioned at described coordinate system, so each described statistics point self has coordinate.Time on seismic data mapping to this coordinate system, the seismic trace of geological data is by himself virtual borderlines on the statistics corresponding with its coordinate point, and described seismic channel data stores by this statistics point.When the seismic channel data that has map upper less than statistics point and map between statistics point time, add up some meeting by nearby principle and the mapping position of this seismic trace are immediate this seismic channel data stored.Described nearby principle is for when a seismic channel data maps between several statistics point, stored by the statistics point nearest with described seismic channel data mapping point, if when occurring that multiple statistics point is identical apart from the distance of described seismic channel data, then in described multiple statistics point, random statistics point stores.Arranged by this kind, described coordinate system can adapt to the change of seismic data acquisition recording geometry left-right situs, survey line encryption, lacks road, lack the complicated situations of change such as arrangement, has extremely strong adaptability.

Described statistics point can comprise integer array and integer variable, and described integer array stores the Taoist monastic name of seismic channel data for storing described statistics point institute, described integer variable for store described add up put store the quantity of seismic channel data.Further, the number of the described integer array storage Taoist monastic name of each described statistics point is less than or equal to 5.First, described integer variable can limit described statistics point can store the maximum quantity of seismic trace, and can the quantity that records of described integer variable consistent with the number needs of the Taoist monastic name that described integer array stores, and then the accuracy of guarantee data.Meanwhile, geophone station coverage rate in this position angle can be counted by adding up described integer variable.

In a preferred embodiment, in the acquisition methods of the common position angle geological data that present embodiment provides, the coordinate range that described step S1 can also be included according to geological data on described coordinate system marks off mapping area.Accordingly, described step S2 comprises: at the statistics point of the built-in pre-defined rule arrangement that is based on of described mapping area that described coordinate system divides.

The shape of described mapping area can be rectangle or circle, and certain the application is total to the mapping area shape of the acquisition methods of position angle geological data not as limit.When described mapping area shape is rectangle, its length and width, by maximum coordinates scope in described geological data, be such as total to single big gun maximum coordinates scope of big gun collection, or CMP concentrate the maximum coordinates scope of single CMP.For common big gun collection, by adding up the maximum coordinates of each single big gun at vertical survey line and cross-track to determine the length and width of described rectangle mapping area.Certainly, when described mapping area is circular, the greater that its diameter is each single big gun in the maximum coordinates of vertical survey line and cross-track.

Further, describedly on described coordinate system, mark off according to the coordinate range of geological data the coordinate range that mapping area can be included according to geological data on described coordinate system mark off mapping area by mesh lines; Accordingly, in described step S2, described pre-defined rule comprises described statistics point and is positioned at grid node locations.

For adapting with the arrangement in actual seismic road, described mesh lines can comprise many parallels parallel to each other and many warps vertical with described horizontal line, and then a large amount of rectangular nodes is formed in described mapping area, described statistics point is positioned at four summits of described rectangular node.Further, the transverse axis of described coordinate system and the longitudinal axis can be two mutually orthogonal mesh lines.The origin position of described coordinate system can be the center of described mapping area.Arranged by this kind, can ensure that four quadrants that the mapping area of described coordinate system is separated to form by described coordinate system have identical area, and be distributed with the statistics point of common number, can ensure that again it can cover 360 degree of omnibearing position angles simultaneously.

Further, the distance often between adjacent two described parallels can be equal, and the distance often between adjacent two described warps can be equal.Certainly, the line-spacing that described distance often between adjacent two described parallels can be concentrated with actual seismic data track is multiple proportion, and the track pitch that the distance often between adjacent two described warps can be concentrated with actual seismic data track is multiple proportion.Concrete, the statistics point identical with its coordinate can be mapped directly to for making most of data of geological data, 0.5 times or 1 times of the line-spacing that described distance often between adjacent two described parallels can be concentrated for actual seismic data track, 0.5 times or 1 times of the track pitch that the distance often between adjacent two described warps can be concentrated for actual seismic data track.Distance between adjacent described parallel can be called mapping area line-spacing, and the distance between adjacent described warp can be called mapping area track pitch.

So the acquisition methods of the common position angle geological data that present embodiment provides can also comprise the track pitch and line-spacing that obtain geological data.When dividing described mapping area, described mesh lines can comprise many parallels parallel to each other and many warps vertical with described parallel, the transverse axis of described coordinate system and the longitudinal axis are two mutually orthogonal described mesh lines, distance often between adjacent two described parallels is equal, and the distance often between adjacent two described warps is equal.

S3, determine position angle to be extracted and set up to be positioned at the index of the statistics point at position angle to be extracted.

The described index set up in this step S3 be follow-up geological data is mapped on described coordinate system after, pass through set up index and can locate rapidly the statistics point being positioned at position angle to be extracted, described index can for described coordinate system being positioned at the coordinate of the statistics point at position angle to be extracted, or described index can also be the order label of the distant relationships according to described statistics point range coordinate initial point.Certainly, described index can artificially define, and for ease of quick indexing, the index number of each statistics point can be tactic natural number etc.

Meanwhile, consider that there is different to be extracted azimuthal demands usually, so many covers can be set up in described index.For different position angles to be extracted, the index set up is different.That is, multiple index has feature one to one for the position angle to be extracted corresponding to it.Certainly, exist in multiple azimuthal situation to be extracted, the index number of the lap in described multiple position angle to be extracted can be consistent.

In a feasible embodiment, please refer to Fig. 2, described step S3 can comprise following sub-step:

S31, in described coordinate system, mark off region, position angle to be extracted.

Described step S31 is for directly to reflect position angle to be extracted in a coordinate system, concrete, also can mark off region, position angle to be extracted in mapping area.Described position angle to be extracted can be determined by described step S31, such as described position angle to be extracted is 30 degree, north by east angular zone, can in described coordinate, the region rotating clockwise 30 degree by y-axis to x-axis in first quartile is region, position angle to be extracted, and the statistics point being positioned at region, described position angle to be extracted is required statistics point.

Certainly, described step S31 is used for determining position angle to be extracted, its execution sequence also performs after Without wishing and execution of step S1 and step S2, this step S31 can be able to perform after carrying out described step S1, certainly, the order between this step S31 and step S2 can not be fixed in the present embodiment.

S32, set up the absorbing boundary equation in region, described position angle to be extracted.

On described coordinate system, region, position angle to be extracted is determined by above-mentioned steps S31, and then need to obtain the statistics point being positioned at described position angle to be extracted, benchmark needs are now needed to go to judge whether described statistics point is positioned at region, described position angle to be extracted, consider based on this kind, can by setting up the absorbing boundary equation in region, described position angle to be extracted, then go to judge by described absorbing boundary equation, and then draw the statistics point being positioned at region, position angle to be extracted.

For the position angle described to be extracted in described step S31 for 30 degree, north by east angular zone, after marking off region, described position angle to be extracted, the absorbing boundary equation in region, position angle to be extracted can be: y >=2x; X >=0.

S33, judge whether described statistics point is positioned at region, described position angle to be extracted according to described absorbing boundary equation.

After obtaining described absorbing boundary equation by described step S32, can by the coordinate of each statistics point be substituting to described absorbing boundary equation, judge whether the orientation of each statistics point meets described absorbing boundary equation, and then know whether this statistics point is positioned at region, described position angle to be extracted.

Certainly, when described statistics point is positioned at the described mapping area divided by mesh lines, because the described parallel of described mesh lines and described warp are all arranged according to certain rule, and then after knowing described absorbing boundary equation, can also superpose according to the distance between described parallel and the distance between described warp the statistics point drawing and be positioned at region, described position angle to be extracted step by step.For the absorbing boundary equation of described step S32, the distance of definition parallel and warp is 1, and y≤100; And then with (0,0) be the first statistics point, then horizontal ordinate adds 1, show that the second statistics point is (1 according to y >=2x, 2), on the basis of the second statistics point, ordinate adds 1 until y≤100, and then show that horizontal ordinate is all statistics points being positioned at region, position angle to be extracted of 1; Then horizontal ordinate adds 1 again, show that the 3rd statistics point is for (2 according to y >=2x, 4), on the basis of the 3rd statistics point, ordinate adds 1 until y≤100, and then show that horizontal ordinate is all statistics points being positioned at region, position angle to be extracted of 2, carry out in turn until draw all statistics points meeting described absorbing boundary equation.

The index that the described statistics that S34, foundation are positioned at position angle to be extracted is put.

On the basis being positioned at the statistics point at described position angle to be extracted obtained in described step S33, index is set up to each described statistics point being positioned at position angle to be extracted.In the present embodiment, index is an order label for the distant relationships according to described statistics point range coordinate initial point, as by the natural numbers such as 1,2,3 to each described statistics point label in turn, until marked the described statistics point being positioned at described position angle to be extracted.

In the embodiment that another is feasible, please refer to Fig. 3, described step S3 can comprise following sub-step:

S31 ＇, calculate described statistics point position angle.

Number due to described statistics point is far smaller than the quantity of the seismic channel data in described geological data, so the computing velocity of this step S31 ＇ extremely can't affect arithmetic speed.Because each statistics point all has the coordinate of self, calculated the position angle on described coordinate system of each described statistics point by this coordinate, this position angle is measured with 0 to 360 degree certainly.

S32 ＇, in described coordinate system, mark off region, position angle to be extracted.

Described step S32 ＇ is for directly to reflect position angle to be extracted in a coordinate system, concrete, also can mark off region, position angle to be extracted in mapping area.Described position angle to be extracted can be determined by described step S31, such as described position angle to be extracted is 30 degree, north by east angular zone, can in described coordinate, the region rotating clockwise 30 degree by y-axis to x-axis in first quartile is region, position angle to be extracted, and the statistics point being positioned at region, described position angle to be extracted is required statistics point.

Certainly, described step S32 ＇ is used for determining position angle to be extracted, its execution sequence being not limited to performs after execution of step S1 and step S2, this step S32 ＇ can be able to perform after carrying out described step S1, certainly, the order between this step S31 and step S2 can not be fixed in the present embodiment.

S33 ＇, according to described statistics point position angle judge whether this statistics point is positioned at region, described position angle to be extracted.

After carrying out described step S31 ＇ and described S32 ＇, know position angle and position angle to be extracted that to add up described in each a little, the position angle added up a little and region, described position angle to be extracted are compared can judge to show whether this statistics point is positioned at region, described position angle to be extracted described in each.Such as, the position angle of a certain statistics point is 30 degree, and region, position angle to be extracted is 0 to 45 degree, and namely its region is 0 degree≤position angle≤45 degree to be extracted, judges 0 degree≤30 degree≤45 degree, and then show that this statistics point is positioned at region, described position angle to be extracted.

The index that the described statistics that S34 ＇, foundation are positioned at position angle to be extracted is put.

On the basis being positioned at the statistics point at described position angle to be extracted obtained in described step S33 ＇, this step S34 ＇ sets up index to each described statistics point being positioned at position angle to be extracted.In the present embodiment, index is an order label for the distant relationships according to described statistics point range coordinate initial point, as by the natural numbers such as 1,2,3 to each described statistics point label in turn, until marked the described statistics point being positioned at described position angle to be extracted.

S4, by seismic data mapping to described coordinate system, described statistics point will be mapped on this statistics point and the seismic channel data be mapped near this statistics point stores.

Described geological data can comprise multiple roads collection, described geological data carry out described seismic data mapping to described coordinate system and the described index of described traversal obtained the seismic channel data that the statistics point that is positioned at position angle to be extracted stores time, described multiple roads collection one by one road collection carries out described seismic channel data seismic data mapping to described coordinate system and the described index of described traversal being obtained the statistics point that is positioned at position angle to be extracted and store, so described coordinate system and described statistics point only need to set up once in the extraction work of whole data, for whole geological data multiple roads collection described in coordinate system and described statistics point can Reusability.

In a specific embodiment, please refer to Fig. 4, described step S4 can comprise following sub-step:

S41, the Taoist monastic name described statistics being put the seismic channel data stored and quantity are removed.

Described step S41 is initialization operation, this step is considering that geological data to be extracted has a large amount of road collection usually, all needed to carry out initialization operation, to prevent the extraction work of the data of this road collection of data influence of collection together before each road collection maps.Described seismic trace quantity is set to zero, and namely described integer variable is set to zero, shows now do not have mapping (enum) data in described coordinate system.Then when each road collection carries out extracted data work, Counts is carried out in circulation, and then ensures that the data counts of this road collection is not by the impact of upper collection together.

S42, each seismic channel data of described geological data to be mapped in described mapping area according to its mapping point.

This step S42 is that each seismic channel data all can have mapping point when mapping, and described mapping point can be identical with the described coordinate added up a little by seismic data mapping to described coordinate system.Consider that described geological data to generally include in common big gun collection, CMP collection, common receiving point gather the combination of one or more.The receiving station coordinate that the mapping point of the seismic channel data of described altogether big gun collection is this seismic channel data deduct shot point coordinate draw the multiple of coordinate figure; The receiving station coordinate that the mapping point of the seismic channel data of described CMP collection is this seismic trace deduct center point coordinate draw the multiple of coordinate figure; The shot point coordinate that the mapping point of the seismic trace of described common receiving point gather is this seismic trace deduct acceptance point coordinate draw the multiple of coordinate figure.

Can find out, when the data pick-up work carrying out common big gun collection, the shot point of described big gun collection is altogether the initial point of described coordinate system; When the data pick-up work carrying out CMP collection, the central point of described CMP collection is the initial point of described coordinate system; When the data pick-up work carrying out common receiving point gather, the acceptance point of described common receiving point gather is the initial point of described coordinate system.

Certainly, because the distance between the distance between parallel adjacent on described coordinate system and adjacent warp can be self-defined, the actual relative coordinate of described mapping point and seismic trace to be mapped is (for common big gun collection, actual relative coordinate is that the receiving station coordinate of this seismic trace deducts shot point coordinate) there are multiple mapping relations, described mapping point is reflected on described coordinate system and is the coordinate that seismic trace to be mapped maps to statistics point.

For single big gun collection altogether, the actual relative coordinate of seismic trace to be mapped is (x1, y1), and the track pitch of described geological data is L, and line-spacing is H; Described mapping point horizontal ordinate can be (x1/L, y1/H), finds corresponding statistics point on coordinate system according to this mapping point (x1/L, y1/H).Leading portion situation and this section of situation are only two specific embodiments of present embodiment, because the mapping relations of described mapping point are various, so the application is not as limit.

The coordinate of road collection belonging to it at the coordinate of described coordinate system and the seismic channel data of mapping of the statistics point on described coordinate system is consistent, and, vertical survey line, the cross-track of described coordinate system and road collection are consistent, so after mapping, the concrete orientation of described seismic trace does not change, and then accurately can reflect the affiliated azimuth information of described seismic trace.

S5, travel through described index and obtain the seismic channel data that the statistics point that is positioned at position angle to be extracted stores.

After geological data maps, the statistics being positioned at position angle to be extracted according to setting up at described step S3 puts corresponding index, travel through region, described position angle to be extracted, and then the data that the statistics point described coordinate system being positioned at region, described position angle to be extracted stores are collected, the common orientation angle gathers needed for seismic trace composition formed after collecting.So, the described index of traversal can be comprised in described step S5 and collect seismic channel data that the statistics point that is positioned at position angle to be extracted stores and be stored as common orientation angle gathers.

Described orientation angle gathers altogether can be two-dimentional orientation angle gathers altogether also can be three-dimensional orientation angle gathers altogether.Described position angle to be extracted can be a certain angle, also can be a certain angular range.When described position angle to be extracted is a certain angle, described orientation angle gathers is altogether two dimension; When described position angle to be extracted is a certain angular range, described orientation angle gathers is altogether three-dimensional.Illustrate, when described position angle to be extracted is 30 degree of position angles, north by east, the seismic channel data that the data of described extraction store for the statistics point dropped on the angle line at 30 degree of angles, north by east; When described position angle to be extracted is 0 degree, north by east to 30 degree of position angles, the seismic channel data that the data of described extraction store to the statistics point in region, 30 degree of position angles for dropping on 0 degree, north by east.Meanwhile, find the statistics point in region, described position angle to be extracted by described index, by extracting the Taoist monastic name of the seismic channel data that described statistics point stores, collecting described seismic channel data according to described Taoist monastic name and storing.

In sum, can find out that the technical scheme that present embodiment provides is passed through on coordinate system, to set up data-storable statistics point, and then make the seismic channel data in geological data according on its virtual borderlines to described coordinate system, then pass through set up index and the seismic channel data stored in statistics point is carried out collection storage, so, the acquisition methods of present embodiment is used only to need region, position angle to be extracted to divide out on described coordinate system, and then can evaluation work be greatly reduced, make the extraction work can carrying out common position angle geological data when processing the geological data of wide-azimuth exploration fast.

Meanwhile, the acquisition methods of the common position angle geological data that present embodiment provides can extract orientation angle gathers altogether to the geological data of arbitrary size efficiently, can replace the method that common orientation angle gathers is extracted in traditional computer azimuth angle again.Three advantages of the technical scheme that present embodiment provides are: the acquisition methods of the common position angle geological data that (1) present embodiment provides directly does not calculate geological data position angle, but set up can full azimuth coordinate system repeatedly, and map seismic trace on described coordinate system, extract a position angle geological data demand altogether and get seismic trace in sector that on described coordinate system, each position angle is corresponding and order is preserved; (2) acquisition methods of common position angle geological data that present embodiment provides once can set up the common orientation angle gathers index of many cover different orientations sizes, position angle splitting schemes corresponding different respectively; (3) situations are seen in the various changes in the acquisition methods energy natural adaptation seismic data acquisition of common position angle geological data that provides of present embodiment, such as arrange number of channels increase and decrease, unidirectional spread and bilateral arrangement mixing, local and receive line increase and decrease encryption etc., unrestricted condition, succinctly efficiently.

Please continue to refer to Fig. 5, be described in detail the principle of the abstracting method of the common position angle geological data that the application provides below with a specific embodiment, the application not with content shown in this embodiment for restriction.

The present embodiment is for the common big gun collection in geological data.Obtain prestack list big gun data, add up the maximum coordinates of each single big gun at vertical survey line and cross-track, be respectively 30000 meters and 2000 meters, and, obtain geological data track pitch and reception line-spacing simultaneously, be respectively 30 meters and 100 meters.

First, set up coordinate system 1, and mark off mapping area 2 according to each single big gun described in the maximum coordinates of vertical survey line and cross-track on described coordinate system 1, described mapping area 2 entirety is rectangle.

Using the maximum coordinates scope of the maximum coordinates in single big gun data as described coordinate system 1, distance between the adjacent warp getting described coordinate system 1 is the distance between 15 meters and adjacent parallel is 50 meters, by described coordinate system 1 coordinate range respectively divided by described 15 meters and described 50 meters, and then by a large amount of parallel and warp, mapping area 2 is carried out gridding, set up the described mapping area 2 divided by mesh lines.The transverse axis (x-axis) of described coordinate system 1 parallels with the vertical survey line of described big gun collection altogether, and the longitudinal axis (y-axis) of described coordinate system 1 parallels with the cross-track of described big gun collection altogether.Certainly, described warp, parallel are equivalent to the described cross-track, the described vertical survey line that are total to big gun collection.

Statistics point 3 in described mapping area 2 on single parallel adds up to 2000, and statistics point 3 number on single warp is 41.Getting coordinate system 1 center is sp location, and described mapping area 2 comprises the identical quadrant of scope, and covers the position angle of 360 degree.Each statistics point 3 in described mapping area 2 is provided with the integer array that can store 5 seismic channel data Taoist monastic names, is used for recording at most the Taoist monastic name of 5 track datas, arranges the seismic channel data number of channels that this statistics of integer variable record puts 3 actual storage simultaneously.

It is be less than track pitch to process two earthquake track pitches that each statistics point 3 is set up integer array, with the situation receiving line encryption, now having multiple seismic trace to be mapped in described mapping area 2 on same statistics point 3, now, described statistics point 3 can by fall thereon and fall in its vicinity and the Taoist monastic name of seismic trace do not dropped on statistics point 3 all preserved.Because coordinate system 1 contains 360 position angles centered by road collection central point, and coordinate system 1 range set is to hold whole seismic traces of any single road collection, therefore the left and right permutations of this coordinate system 1 energy natural adaptation seismic data acquisition recording geometry, survey line are encrypted, lack road, are lacked the complicated situations of change such as arrangement, have extremely strong adaptability.

Coordinate system 1 schematic diagram that Fig. 5 sets up for the present embodiment.In the mapping area 2 of described coordinate system 1, circle and black splotch are the statistics point 3 after coordinate system 1 grid discretization.The asterisk at described mapping area 2 center is sp location, and its coordinate is (0,0).Described mapping area 2 contains 360 degree of orientation.The horizontal direction of described coordinate system 1 is transverse axis, i.e. vertical line direction, coordinate range be (-4,0) to (4,0), vertical direction is the longitudinal axis, i.e. cross-track direction, and coordinate range is that (0 ,-2) arrive (0,2).

Then, the index of the statistics point 3 described coordinate system 1 being positioned at position angle to be extracted is set up.Specify position angle to be extracted to be 0 to 360 degree, angle intervals is 45, equidistant 8 position angles to be extracted, namely needs 8 orientation angle gathers altogether.Dotted line in described mapping area 2 is angular dividing line, is positioned at 45 degree and 45 degree of integral multiple positions, coordinate system 1 is divided into totally 8 position angles to be extracted, sector (1) to sector (8).Extract between separatrix, that is to say that the statistics in the mapping area 2 of the sector region that position angle is corresponding puts 3 coordinates as correspondence azimuthal index to be extracted, store for subsequent use.This index will be used for setting up three-dimensional orientation angle gathers altogether.

Three-dimensional as shown in Figure 5 altogether azimuthal coordinate be 1 index comprise coordinate system 1 coordinate of sector (1) to sector (8) in totally 8 position angles to be extracted.Index includes the coordinate of following point in sector (1): (1,0), (2,0), (3,0), (4,0), (2,1), (3,1), (4,1), (4,2); Following point coordinate is included: (0,1), (1,1), (0,2), (1,2), (2,2), (3,2) in sector (2).Index in other sector internal coordinate values by that analogy.Certainly, the statistics point 3 be positioned on separatrix belongs to the position angle to be extracted comprising this angle value.

Again by described seismic data mapping to described extraction coordinate system 1, the seismic channel data mapped on this statistics point 3 and near this statistics point 3 stores by the statistics point 3 on described extraction coordinate system 1.

First the integer variable on statistics point 3 each in the full azimuth coordinate system 1 of above-mentioned foundation is set to 0, represents now there is no mapping (enum) data in coordinate system 1.

Then, the acceptance point coordinate of each seismic trace concentrated by single big gun altogether deducts the mapping point that shot point coordinate obtains this seismic trace, be mapped in described mapping area 2 again according to described mapping point, and the integer variable of the statistics of mapping position point 3 is added 1, in the integer array of this statistics point 3, record the Taoist monastic name of current mapping simultaneously.After mapping terminates, the quantity of the seismic trace that integer variable value i.e. this statistics point 3 of the statistics point 3 in described mapping area 2 stores, the Taoist monastic name being the seismic trace that this statistics point 3 stores that the integer array of statistics point 3 records.

In mapping area 2 as shown in Figure 5, black color dots indicates that seismic trace maps, and in hollow dots denotation coordination system 1, this position does not have seismic trace to map.In the present embodiment, black splotch (having seismic trace statistics point 3) upper seismic trace maps number of channels and is 1, therefore do not mark especially in figure.

Finally, travel through described index the seismic channel data that the statistics point 3 being positioned at described position angle to be extracted stores is carried out collecting and is stored as common orientation angle gathers.

In single road collection, all seismic traces are after coordinate system 1 is shown and penetrated, according to 8 position angles to be extracted of 45 degree of spacing corresponding to set up index, described mapping area 2 is traveled through by angle, find the statistics point 3 of the mapping (enum) data in position angle to be extracted that in described mapping area 2, described index is corresponding, geological data is extracted with the Taoist monastic name that the statistics point 3 of mapping (enum) data stores, seismic trace corresponding for each position angle to be extracted is stored together, becomes the common orientation angle gathers that this position angle is corresponding.After being disposed to position angles all in index, the common bearing data that this road is concentrated has extracted.In mapping area 2 as shown in Figure 1, sector (4) map without seismic trace, show that the position angle at (4) place is without geological data.There is geological data at other position angle.

When described geological data comprises multiple road collection, to each the road collection in whole geological data, traversal step S4, to step S5, completes the screening that common orientation angle gathers is concentrated in whole road.

Please refer to Fig. 6, a kind of embodiment of the application also provides a kind of acquisition device of position angle geological data altogether, comprises with lower module:

First sets up module 100, and it is for setting up coordinate system; Described coordinate system is for mapping geological data.

Described first sets up in module 100, and described coordinate system can be orthogonal coordinate system, and it comprises x-axis (also can be called transverse axis) and y-axis.Described geological data comprises multiple roads collection usually, such as comprises multiple big gun collection altogether, multiple CMP collection, multiple common receiving point gather etc.The geophone station normally arranged during image data in seismic exploration carries out grid arrangement according to vertical survey line and cross-track, so the seismic trace that each geophone station gathers has its corresponding coordinate.There is its corresponding track pitch and line-spacing in each seismic trace of each road collection.So the transverse axis of described coordinate system (x-axis) represents the described vertical survey line concentrated, the longitudinal axis (y-axis) of described coordinate system represents the described cross-track concentrated, and then makes need not can directly mapping in converted coordinate orientation of described geological data.

Second sets up module 200, and it is for the statistics point in the built-in pre-defined rule arrangement that is based on of described coordinate system; Described statistics point is mapped on this statistics point and the seismic channel data be mapped near this statistics point for storing; The quantity of coordinate range and described statistics point that the coordinate range that described pre-defined rule comprises described statistics point is more than or equal to geological data is less than the quantity of the seismic channel data of geological data.

Set up in module 200 described second, described statistics point is owing to being positioned at described coordinate system, so each described statistics point self has coordinate.Time on seismic data mapping to this coordinate system, the seismic trace of geological data is by himself virtual borderlines on the statistics corresponding with its coordinate point, and described seismic channel data stores by this statistics point.When the seismic channel data that has map upper less than statistics point and map between statistics point time, add up some meeting by nearby principle and the mapping position of this seismic trace are immediate this seismic channel data stored.Described nearby principle is for when a seismic channel data maps between several statistics point, stored by the statistics point nearest with described seismic channel data mapping point, if when occurring that multiple statistics point is identical apart from the distance of described seismic channel data, then in described multiple statistics point, random statistics point stores.Arranged by this kind, described coordinate system can adapt to the change of seismic data acquisition recording geometry left-right situs, survey line encryption, lacks road, lack the complicated situations of change such as arrangement, has extremely strong adaptability.

In a preferred embodiment, in the acquisition device of the common position angle geological data that present embodiment provides, described first sets up module 100 can also comprise Region dividing module, and described Region dividing module marks off mapping area for the coordinate range according to geological data on described coordinate system.Accordingly, described second sets up module 200 comprises the second son and sets up module, and described second son sets up the statistics point of described mapping area built-in be based on pre-defined rule arrangement of module for dividing at described coordinate system.

In described Region dividing module, the shape of described mapping area can be rectangle or circle, and certain the application is total to the mapping area shape of the acquisition methods of position angle geological data not as limit.When described mapping area shape is rectangle, its length and width, by maximum coordinates scope in described geological data, be such as total to single big gun maximum coordinates scope of big gun collection, or CMP concentrate the maximum coordinates scope of single CMP.For common big gun collection, by adding up the maximum coordinates of each single big gun at vertical survey line and cross-track to determine the length and width of described rectangle mapping area.Certainly, when described mapping area is circular, the greater that its diameter is each single big gun in the maximum coordinates of vertical survey line and cross-track.

Further, described Region dividing module can comprise Region dividing submodule, and described Region dividing submodule marks off mapping area for the coordinate range according to geological data on described coordinate system by mesh lines; Accordingly, set up in module 200 described second, described pre-defined rule comprises described statistics point and is positioned at grid node locations.

In described Region dividing submodule, consider as adapting with the arrangement in actual seismic road, described mesh lines can comprise many parallels parallel to each other and many warps vertical with described horizontal line, and then a large amount of rectangular nodes is formed in described mapping area, described statistics point is positioned at four summits of described rectangular node.Further, the transverse axis of described coordinate system and the longitudinal axis can be two mutually orthogonal mesh lines.The origin position of described coordinate system can be the center of described mapping area.Arranged by this kind, can ensure that four quadrants that the mapping area of described coordinate system is separated to form by described coordinate system have identical area, and be distributed with the statistics point of common number, can ensure that again it can cover 360 degree of omnibearing position angles simultaneously.

Further, the distance often between adjacent two described parallels can be equal, and the distance often between adjacent two described warps can be equal.Certainly, the line-spacing that described distance often between adjacent two described parallels can be concentrated with actual seismic data track is multiple proportion, and the track pitch that the distance often between adjacent two described warps can be concentrated with actual seismic data track is multiple proportion.Concrete, the statistics point identical with its coordinate can be mapped directly to for making most of data of geological data, 0.5 times or 1 times of the line-spacing that described distance often between adjacent two described parallels can be concentrated for actual seismic data track, 0.5 times or 1 times of the track pitch that the distance often between adjacent two described warps can be concentrated for actual seismic data track.

So the acquisition device of the common position angle geological data that present embodiment provides can also comprise acquisition module, described acquisition module is for obtaining track pitch and the line-spacing of geological data.In described Region dividing submodule, described mesh lines can comprise many parallels parallel to each other and many warps vertical with described parallel, the transverse axis of described coordinate system and the longitudinal axis are two mutually orthogonal described mesh lines, distance often between adjacent two described parallels is equal, and the distance often between adjacent two described warps is equal.

3rd sets up module 300, and it is for determining position angle to be extracted and setting up the index being positioned at the statistics point at position angle to be extracted.

Described 3rd sets up in module 300, the described index set up be follow-up geological data is mapped on described coordinate system after, pass through set up index and can locate rapidly the statistics point being positioned at position angle to be extracted, described index can for described coordinate system being positioned at the coordinate of the statistics point at position angle to be extracted, or described index can also be the order label of the distant relationships according to described statistics point range coordinate initial point.Certainly, described index can artificially define, and for ease of quick indexing, the index number of each statistics point can be tactic natural number etc.

In a feasible embodiment, please refer to Fig. 7, the described 3rd sets up module 300 can comprise following submodule:

Position angle divides A module 301, and it for marking off region, position angle to be extracted in described coordinate system.

Described position angle divides A module 301 for directly to be reflected in a coordinate system at position angle to be extracted, concrete, also can mark off region, position angle to be extracted in mapping area.Described position angle to be extracted can be determined by described step S31, such as described position angle to be extracted is 30 degree, north by east angular zone, can in described coordinate, the region rotating clockwise 30 degree by y-axis to x-axis in first quartile is region, position angle to be extracted, and the statistics point being positioned at region, described position angle to be extracted is required statistics point.

Certainly, described position angle divides A module 301 for determining position angle to be extracted, its execution sequence Without wishing with execute described first and set up after module 100 and described second sets up module 200 and perform, this position angle divides after A module 301 can set up module 100 carrying out described first and can perform, certainly, this position angle division A module 301 and the second order set up between module 200 can not be fixed in the present embodiment.

Establishing equation module 302, it is for setting up the absorbing boundary equation in region, described position angle to be extracted.

Divide A module 301 by above-mentioned position angle and on described coordinate system, determine region, position angle to be extracted, and then need to obtain the statistics point being positioned at described position angle to be extracted, benchmark needs are now needed to go to judge whether described statistics point is positioned at region, described position angle to be extracted, consider based on this kind, described establishing equation module 302 can by setting up the absorbing boundary equation in region, described position angle to be extracted, then go to judge by described absorbing boundary equation, and then draw the statistics point being positioned at region, position angle to be extracted.

For the position angle described to be extracted in described position angle division A module 301 for 30 degree, north by east angular zone, after marking off region, described position angle to be extracted, the absorbing boundary equation in region, position angle to be extracted can be: y >=2x; X >=0.

Judge A module 303, it is for judging according to described absorbing boundary equation whether described statistics point is positioned at region, described position angle to be extracted.

After obtaining described absorbing boundary equation by described establishing equation module 302, can by the coordinate of each statistics point be substituting to described absorbing boundary equation, judge whether the orientation of each statistics point meets described absorbing boundary equation, and then know whether this statistics point is positioned at region, described position angle to be extracted.

Certainly, when described statistics point is positioned at the described mapping area divided by mesh lines, because the described parallel of described mesh lines and described warp are all arranged according to certain rule, and then after knowing described absorbing boundary equation, can also superpose according to the distance between described parallel and the distance between described warp the statistics point drawing and be positioned at region, described position angle to be extracted step by step.For the absorbing boundary equation of described establishing equation module 302, the distance of definition parallel and warp is 1, and y≤100; And then with (0,0) be the first statistics point, then horizontal ordinate adds 1, show that the second statistics point is (1 according to y >=2x, 2), on the basis of the second statistics point, ordinate adds 1 until y≤100, and then show that horizontal ordinate is all statistics points being positioned at region, position angle to be extracted of 1; Then horizontal ordinate adds 1 again, show that the 3rd statistics point is for (2 according to y >=2x, 4), on the basis of the 3rd statistics point, ordinate adds 1 until y≤100, and then show that horizontal ordinate is all statistics points being positioned at region, position angle to be extracted of 2, carry out in turn until draw all statistics points meeting described absorbing boundary equation.

A module 304 set up in index, and it is for setting up the index of the described statistics point being positioned at position angle to be extracted.

What obtain in described judgement A module 303 is positioned on the basis of statistics point at described position angle to be extracted, sets up index to each described statistics point being positioned at position angle to be extracted.In the present embodiment, index is an order label for the distant relationships according to described statistics point range coordinate initial point, as by the natural numbers such as 1,2,3 to each described statistics point label in turn, until marked the described statistics point being positioned at described position angle to be extracted.

In the embodiment that another is feasible, please refer to Fig. 8, the described 3rd sets up module 300 can comprise following submodule:

Position angle computing module S301 ＇, it is for calculating the position angle of described statistics point.

Number due to described statistics point is far smaller than the quantity of the seismic channel data in described geological data, so the computing velocity of this position angle computing module S301 ＇ extremely can't affect arithmetic speed.Because each statistics point all has the coordinate of self, calculated the position angle on described coordinate system of each described statistics point by this coordinate, this position angle is measured with 0 to 360 degree certainly.

Position angle divides B module S302 ＇, and it for marking off region, position angle to be extracted in described coordinate system.

Described position angle divides B module S302 ＇ for directly to be reflected in a coordinate system at position angle to be extracted, concrete, also can mark off region, position angle to be extracted in mapping area.Described position angle to be extracted can be determined by described step S31, such as described position angle to be extracted is 30 degree, north by east angular zone, can in described coordinate, the region rotating clockwise 30 degree by y-axis to x-axis in first quartile is region, position angle to be extracted, and the statistics point being positioned at region, described position angle to be extracted is required statistics point.

Certainly, described position angle divides B module S302 ＇ and is used for determining position angle to be extracted, its execution sequence Without wishing with execute described first and set up after module 100 and described second sets up module 200 and perform, this position angle divides after B module S302 ＇ can set up module 100 carrying out described first and can perform, certainly, this position angle division B module S302 ＇ and the described second order set up between module 200 can not be fixed in the present embodiment.

Judge B module S303 ＇, according to the position angle of described statistics point, it is for judging whether this statistics point is positioned at region, described position angle to be extracted.

Carrying out having known after described position angle computing module S301 ＇ and position angle divide B module S302 ＇ position angle and position angle to be extracted that to add up described in each a little, the position angle added up a little and region, described position angle to be extracted being compared can judge to show whether this statistics point is positioned at region, described position angle to be extracted described in each.Such as, the position angle of a certain statistics point is 30 degree, and region, position angle to be extracted is 0 to 45 degree, and namely its region is 0 degree≤position angle≤45 degree to be extracted, judges 0 degree≤30 degree≤45 degree, and then show that this statistics point is positioned at region, described position angle to be extracted.

B module S304 ＇ set up in index, and it is for setting up the index of the described statistics point being positioned at position angle to be extracted.

On the basis being positioned at the statistics point at described position angle to be extracted obtained in described judgement B module S303 ＇, this index is set up B module S304 ＇ and is set up index to each described statistics point being positioned at position angle to be extracted.In the present embodiment, index is an order label for the distant relationships according to described statistics point range coordinate initial point, as by the natural numbers such as 1,2,3 to each described statistics point label in turn, until marked the described statistics point being positioned at described position angle to be extracted.

Data mapping module S402400, its for by seismic data mapping to described coordinate system, described statistics point will be mapped on this statistics point and the seismic channel data be mapped near this statistics point stores.

In described data mapping module S402400, described geological data can comprise multiple roads collection, described geological data carry out described seismic data mapping to described coordinate system and the described index of described traversal obtained the seismic channel data that the statistics point that is positioned at position angle to be extracted stores time, described multiple roads collection one by one road collection carries out described seismic channel data seismic data mapping to described coordinate system and the described index of described traversal being obtained the statistics point that is positioned at position angle to be extracted and store, so described coordinate system and described statistics point only need to set up once in the extraction work of whole data, for whole geological data multiple roads collection described in coordinate system and described statistics point can Reusability.

In a specific embodiment, please refer to Fig. 9, described data mapping module S402400 can comprise following submodule:

Initialization module S401, it is for adding up a Taoist monastic name of the seismic channel data stored and quantity removing by described.

Described initialization module S401 performs initialization operation to statistics point, this step is considering that geological data to be extracted has a large amount of road collection usually, all needed to carry out initialization operation, to prevent the extraction work of the data of this road collection of data influence of collection together before each road collection maps.Described seismic trace quantity is set to zero, and namely described integer variable is set to zero, shows now do not have mapping (enum) data in described coordinate system.Then when each road collection carries out extracted data work, Counts is carried out in circulation, and then ensures that the data counts of this road collection is not by the impact of upper collection together.

Mapping block S402, it is for mapping to described mapping area by each seismic channel data of described geological data according to its mapping point.

This mapping block S402 is that each seismic channel data all can have mapping point when mapping, and described mapping point can be identical with the described coordinate added up a little by seismic data mapping to described coordinate system.Consider that described geological data to generally include in common big gun collection, CMP collection, common receiving point gather the combination of one or more.The receiving station coordinate that the mapping point of the seismic channel data of described altogether big gun collection is this seismic channel data deduct shot point coordinate draw the multiple of coordinate figure; The receiving station coordinate that the mapping point of the seismic channel data of described CMP collection is this seismic trace deduct center point coordinate draw the multiple of coordinate figure; The shot point coordinate that the mapping point of the seismic trace of described common receiving point gather is this seismic trace deduct acceptance point coordinate draw the multiple of coordinate figure.

Data acquisition module 500, it obtains for traveling through described index the seismic channel data that the statistics point that is positioned at position angle to be extracted stores.

After geological data maps, the statistics being positioned at position angle to be extracted according to setting up at described step S3 puts corresponding index, travel through region, described position angle to be extracted, and then the data that the statistics point described coordinate system being positioned at region, described position angle to be extracted stores are collected, the common orientation angle gathers needed for seismic trace composition formed after collecting.So described data acquisition module 500 can comprise collection memory module, described collection memory module is collected seismic channel data that the statistics point that is positioned at position angle to be extracted stores for traveling through described index and is stored as common orientation angle gathers.

More than show and describe ultimate principle of the present invention, principal character and advantage of the present invention.The technician of the industry should understand; the present invention is not restricted to the described embodiments; what describe in above-described embodiment and instructions just illustrates principle of the present invention; without departing from the spirit and scope of the present invention; the present invention also has various changes and modifications, and these changes and improvements all fall in the scope of protection of present invention.Application claims protection domain is defined by appending claims and equivalent thereof.

Claims

1. be total to an acquisition methods for position angle geological data, it is characterized in that, comprising:

2. the acquisition methods of position angle geological data altogether as claimed in claim 1, is characterized in that: describedly set up the coordinate range that coordinate system is included according to geological data on described coordinate system and mark off mapping area;

3. the acquisition methods of position angle geological data altogether as claimed in claim 2, it is characterized in that: described statistics point comprises integer array and integer variable, described integer array stores the Taoist monastic name of seismic channel data for storing described statistics point institute, described integer variable for store described add up put store the quantity of seismic channel data.

4. the acquisition methods of position angle geological data altogether as claimed in claim 3, is characterized in that: the number of the described integer array storage Taoist monastic name of each described statistics point is less than or equal to 5.

5. the acquisition methods of position angle geological data altogether as claimed in claim 2, is characterized in that: describedly on described coordinate system, mark off the coordinate range that mapping area is included according to geological data on described coordinate system according to the coordinate range of geological data mark off mapping area by mesh lines;

6. the acquisition methods of position angle geological data altogether as claimed in claim 5, is characterized in that: described acquisition methods also comprises the track pitch and line-spacing that obtain geological data;

7. the acquisition methods of position angle geological data altogether as claimed in claim 1, is characterized in that, describedly determines position angle to be extracted and set up the index being positioned at the statistics point at position angle to be extracted to comprise:

8. the acquisition methods of position angle geological data altogether as claimed in claim 1, is characterized in that, describedly determines position angle to be extracted and set up the index being positioned at the statistics point at position angle to be extracted to comprise:

Calculate the position angle of described statistics point;

9. the acquisition methods of position angle geological data altogether as claimed in claim 1, is characterized in that: described index is the coordinate of described statistics point, or described index is the order label of the distant relationships according to described statistics point range coordinate initial point.

10. the acquisition methods of position angle geological data altogether as claimed in claim 3, is characterized in that, is describedly comprised by seismic data mapping to described coordinate system:

The acquisition methods of 11. common position angle geological datas as described in claim 1 to 10 any one, it is characterized in that: described geological data comprises multiple roads collection, described geological data carry out described seismic data mapping to described coordinate system and described traversal described index obtained the seismic channel data that the statistics point that is positioned at position angle to be extracted stores time, described multiple roads collection one by one road collection carries out that described extremely described coordinate system and the described index of described traversal obtain the seismic channel data that the statistics point that is positioned at position angle to be extracted stores by seismic data mapping.

The acquisition methods of 12. as claimed in claim 10 altogether position angle geological datas, is characterized in that: described geological data to comprise in common big gun collection, CMP collection, common receiving point gather one or more combination; The receiving station coordinate that the mapping point of the seismic channel data of described altogether big gun collection is this seismic channel data deduct shot point coordinate draw the multiple of coordinate figure; The receiving station coordinate that the mapping point of the seismic channel data of described CMP collection is this seismic trace deduct center point coordinate draw the multiple of coordinate figure; The shot point coordinate that the mapping point of the seismic trace of described common receiving point gather is this seismic trace deduct acceptance point coordinate draw the multiple of coordinate figure.

The acquisition methods of 13. as claimed in claim 1 altogether position angle geological datas, is characterized in that: described traversal described index obtains seismic channel data that the statistics point that is positioned at position angle to be extracted stores and comprises and travel through described index and collect seismic channel data that the statistics point that is positioned at position angle to be extracted stores and be stored as common orientation angle gathers.

The acquisition device of 14. 1 kinds of common position angle geological datas, is characterized in that, comprising:

The acquisition device of 15. position angle geological datas altogether as claimed in claim 14, it is characterized in that: described first sets up module comprises Region dividing module, and described Region dividing module marks off mapping area for the coordinate range according to geological data on described coordinate system;

The acquisition device of 16. position angle geological datas altogether as claimed in claim 15, it is characterized in that: described statistics point comprises integer array and integer variable, described integer array stores the Taoist monastic name of seismic channel data for storing described statistics point institute, described integer variable for store described add up put store the quantity of seismic channel data.

The acquisition device of 17. position angle geological datas altogether as claimed in claim 15, it is characterized in that: described Region dividing module comprises Region dividing submodule, the coordinate range that described Region dividing submodule is included according to geological data on described coordinate system marks off mapping area by mesh lines;

The acquisition device of 18. position angle geological datas altogether as claimed in claim 17, it is characterized in that: described acquisition device also comprises acquisition module, described acquisition module is for obtaining track pitch and the line-spacing of geological data;

The acquisition device of 19. position angle geological datas altogether as claimed in claim 14, it is characterized in that, the described 3rd sets up module comprises:

The acquisition device of 20. position angle geological datas altogether as claimed in claim 14, it is characterized in that, the described 3rd sets up module comprises:

The acquisition device of 21. position angle geological datas altogether as claimed in claim 16, it is characterized in that, described data mapping module comprises:

The acquisition device of 22. common position angle geological datas as described in claim 14 to 21 any one, it is characterized in that: described geological data comprises multiple roads collection, when described geological data is by described data mapping module and described data acquisition module process, described multiple roads collection one by one road collection by described data mapping module and described data acquisition module process.

The acquisition device of 23. as claimed in claim 21 altogether position angle geological datas, is characterized in that: described geological data to comprise in common big gun collection, CMP collection, common receiving point gather one or more combination; The receiving station coordinate that the mapping point of the seismic channel data of described altogether big gun collection is this seismic channel data deduct shot point coordinate draw the multiple of coordinate figure; The receiving station coordinate that the mapping point of the seismic channel data of described CMP collection is this seismic trace deduct center point coordinate draw the multiple of coordinate figure; The shot point coordinate that the mapping point of the seismic trace of described common receiving point gather is this seismic trace deduct acceptance point coordinate draw the multiple of coordinate figure.

The acquisition device of 24. position angle geological datas altogether as claimed in claim 14, it is characterized in that: described data acquisition module comprises collection memory module, described collection memory module is collected seismic channel data that the statistics point that is positioned at position angle to be extracted stores for traveling through described index and is stored as common orientation angle gathers.