CN105388513A - Method and device for establishing earthquake acquisition and observation system - Google Patents
Method and device for establishing earthquake acquisition and observation system Download PDFInfo
- Publication number
- CN105388513A CN105388513A CN201510671042.2A CN201510671042A CN105388513A CN 105388513 A CN105388513 A CN 105388513A CN 201510671042 A CN201510671042 A CN 201510671042A CN 105388513 A CN105388513 A CN 105388513A
- Authority
- CN
- China
- Prior art keywords
- orientation
- point
- recording geometry
- shot
- radius
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/02—Generating seismic energy
- G01V1/104—Generating seismic energy using explosive charges
- G01V1/13—Arrangements or disposition of charges to produce a desired pattern in space or time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/20—Arrangements of receiving elements, e.g. geophone pattern
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/14—Signal detection
- G01V2210/142—Receiver location
- G01V2210/1425—Land surface
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method and a device for establishing an earthquake acquisition and observation system, wherein the method comprises the steps of acquiring a pre-established first observation system which is orthogonally arranged on a to-be-observed work area; determining one or a plurality of azimuths with sparse shot-geophone distance distribution according to the surface element attribute of the first observation system; adding a shot-geophone point on each of the one or a plurality of azimuths, and furthermore performing shot arrangement on all shot points in the first observation system according to the rectangular arrangement dimension of the first observation system after the shot-geophone points are added, thereby obtaining a second observation system. The method and the device settle a problem of sparse shot-geophone distance distribution in an orthogonally configured observation system in prior art and realizes technical effects of relatively uniform azimuth and no abnormal sparse azimuth in the earthquake acquisition and observation system.
Description
Technical field
The present invention relates to technical field of geophysical exploration, particularly a kind of method for building up of earthquake-capturing recording geometry and device.
Background technology
Along with the development of computer equipment and seismic exploration equipment, earthquake data acquisition enters " high density wide-azimuth " epoch, and some Earthquakes list big guns have employed the record in 5 to 10 ten thousand roads, define high resolution space data collection.
Current, orthogonal laying recording geometry method for designing is a kind of main earthquake-capturing recording geometry method for designing, and this recording geometry is laid simple, is conducive to linear voice compacting, can realizes wide-azimuth earthquake-capturing simultaneously.But there is an obvious weak point in orthogonal laying recording geometry method for designing, is mainly manifested in: some position angle geophone offset distribution is comparatively sparse, and namely geophone offset cannot reach and be uniformly distributed in 360 ° of orientation.
But high density, wide-azimuth seismic data, in order to carry out reservoir prediction and reservoir description better, require that offset distance can be uniformly distributed in plane 360 ° of orientation.
Cannot ensure, in the equally distributed problem in each orientation, not yet to propose effective solution at present for geophone offset in above-mentioned seismic observation system.
Summary of the invention
Embodiments provide a kind of method for building up of earthquake-capturing recording geometry, to solve the problem that in prior art, in earthquake-capturing recording geometry, some orientation geophone offset is sparse, the method comprises:
Obtain the first recording geometry of the orthogonal laying in the work area to be measured set up in advance;
According to the bin attributes of described first recording geometry, determine that geophone offset distributes sparse one or more orientation;
Described one or more orientation increases big gun examine, and to increase big gun examine after the first recording geometry in all shot points blow out according to the rectangular arranged size of described first recording geometry, obtain the second recording geometry.
In one embodiment, described one or more orientation increases big gun and examines, comprising:
Determine radius of a circle;
Following operation is performed to each shot point in described first recording geometry:
Centered by current shot point, with described radius for radius draws circle;
By the point of each orientation on the circle drawn in described one or more orientation, as the shot point increased;
Following operation is performed to each geophone station in described first recording geometry:
Centered by current geophone station, with described radius for radius draws circle;
By the point of each orientation on the circle drawn in described one or more orientation, as the geophone station increased.
In one embodiment, determine radius of a circle, comprising:
Obtain acceptance point distance and the shot point distance of described first recording geometry;
The size of described acceptance point distance and described shot point distance is compared;
When described acceptance point is apart from being greater than described shot point distance, using the N of described shot point distance doubly as radius, wherein, N is positive integer;
When described acceptance point is apart from being less than described shot point distance, using the M of described acceptance point distance doubly as radius, wherein, M is positive integer.
In one embodiment, by the point of each orientation on the circle drawn in described one or more orientation, after the shot point increased, described method also comprises:
Determine whether the shot point quantity increased exceeds predetermined quantity;
If exceeded, then stress and strain model is carried out to the point of each orientation on the circle drawn in described one or more orientation;
Ask for the combination center point in each grid;
Using the combination center point in each grid of determining as the shot point that will increase;
By the point of each orientation on the circle drawn in described one or more orientation, after the geophone station increased, described method also comprises:
Determine whether the geophone station quantity increased exceeds predetermined quantity;
If exceeded, then stress and strain model is carried out to the point of each orientation on the circle drawn in described one or more orientation;
Ask for the combination center point in each grid;
Using the combination center point in each grid of determining as the geophone station that will increase.
In one embodiment, at the bin attributes according to described first recording geometry, determine that geophone offset distributes after sparse one or more orientation, described method also comprises:
Determine the orientation whether comprising the vertical and horizontal place of survey line in one or more orientation that the distribution that obtains is sparse;
If do not comprised, then using the orientation at the vertical and horizontal place of described survey line also as the sparse orientation of distribution.
The embodiment of the present invention additionally provides a kind of apparatus for establishing of earthquake-capturing recording geometry, and to solve the problem that in prior art, in earthquake-capturing recording geometry, some orientation geophone offset is sparse, this device comprises:
Acquisition module, for obtaining the first recording geometry of the orthogonal laying in the work area to be measured set up in advance;
Orientation determination module, for the bin attributes according to described first recording geometry, determines that geophone offset distributes sparse one or more orientation;
Add module, examine for increasing big gun in described one or more orientation, and to increase big gun examine after the first recording geometry in all shot points blow out according to the rectangular arranged size of described first recording geometry, obtain the second recording geometry.
In one embodiment, described interpolation module comprises:
Radius determining unit, for determining radius of a circle;
Add shot point unit, for performing following operation to each shot point in described first recording geometry: centered by current shot point, with described radius for radius draws circle; By the point of each orientation on the circle drawn in described one or more orientation, as the shot point increased;
Add geophone station unit, for performing following operation to each geophone station in described first recording geometry: centered by current geophone station, with described radius for radius draws circle; By the point of each orientation on the circle drawn in described one or more orientation, as the geophone station increased.
In one embodiment, described radius determining unit, comprising:
Point is apart from determining subelement, for obtaining acceptance point distance and the shot point distance of described first recording geometry;
Size compares subelement, for comparing the size of described acceptance point distance and described shot point distance;
First determines subelement, and for when described acceptance point distance is greater than described shot point distance, using the N of described shot point distance doubly as radius, wherein, N is positive integer;
Second determines subelement, and for when described acceptance point distance is less than described shot point distance, using the M of described acceptance point distance doubly as radius, wherein, M is positive integer.
In one embodiment, described interpolation shot point unit, also for by the point of each orientation on the circle drawn in described one or more orientation, after the shot point increased, determines whether the shot point quantity increased exceeds predetermined quantity; If exceeded, then stress and strain model is carried out to the point of each orientation on the circle drawn in described one or more orientation; Ask for the combination center point in each grid; Using the combination center point in each grid of determining as the shot point that will increase;
Described interpolation geophone station unit, also for by the point of each orientation on the circle drawn in described one or more orientation, after the geophone station increased, determines whether the geophone station quantity increased exceeds predetermined quantity; If exceeded, then stress and strain model is carried out to the point of each orientation on the circle drawn in described one or more orientation; Ask for the combination center point in each grid; Using the combination center point in each grid of determining as the geophone station that will increase.
In one embodiment, described orientation determination module is concrete also at the bin attributes according to described first recording geometry, determine that geophone offset distributes after sparse one or more orientation, determine the orientation whether comprising the vertical and horizontal place of survey line in one or more orientation that the distribution that obtains is sparse; If do not comprised, then using the orientation at the vertical and horizontal place of described survey line also as the sparse orientation of distribution.
In embodiments of the present invention, by the sparse orientation that distributes in the seismic observation system of determining built orthogonal laying, the orientation that the distribution determined is sparse increases big gun examine, and increase big gun is examined the recording geometry of rear generation as final recording geometry, thus efficiently solve the problem that in prior art, in earthquake-capturing recording geometry, some orientation geophone offset is sparse, reach and make each orientation in final earthquake-capturing recording geometry relatively uniform, there is not the technique effect in abnormal sparse orientation.
Accompanying drawing explanation
Accompanying drawing described herein is used to provide a further understanding of the present invention, forms a application's part, does not form limitation of the invention.In the accompanying drawings:
Fig. 1 is the method flow diagram of the method for building up of earthquake-capturing recording geometry according to the embodiment of the present invention;
Fig. 2 is the azimuthal distribution schematic diagram not limitting the orthogonal arranging system geophone offset of offset distance according to the embodiment of the present invention;
Fig. 3 is the azimuthal distribution schematic diagram of the orthogonal arranging system geophone offset of limit offset distance according to the embodiment of the present invention;
Fig. 4 is according to the point in the mapping on the circle of the formation of the embodiment of the present invention and circle;
Fig. 5 is the acceptance point and shot point distribution schematic diagram that are formed according to the combination center of the embodiment of the present invention;
Fig. 6 is the acceptance point and the shot point distribution schematic diagram that equal line-spacing according to the some distance of the embodiment of the present invention;
Fig. 7 is the azimuthal distribution schematic diagram generated apart from the acceptance point and shot point that equal line-spacing according to the point of the embodiment of the present invention;
Fig. 8 is the azimuthal distribution schematic diagram that the acceptance point that formed according to the combination center of the embodiment of the present invention and shot point generate;
Fig. 9 is according to the acceptance point of the embodiment of the present invention and shot point distribution and active arranged distribution schematic diagram;
Figure 10 is the azimuthal distribution schematic diagram laying geophone offset in bin according to the conventional orthorhombic of the embodiment of the present invention;
Figure 11 is acceptance point local distribution schematic diagram after the gridding according to the embodiment of the present invention;
Figure 12 is shot point local distribution schematic diagram after the gridding according to the embodiment of the present invention;
Figure 13 is whole shot point and acceptance point and active arranged distribution schematic diagram after the gridding according to the embodiment of the present invention;
Figure 14 is to compensating and the azimuthal distribution schematic diagram of geophone offset in bin that orthogonal laying recording geometry after gridding analyzes according to the bin of the employing 25 meters × 25 meters of the embodiment of the present invention;
Figure 15 is acceptance point local distribution schematic diagram when carrying out orthogonal laying according to the employing of the embodiment of the present invention identical shot point density and acceptance point density;
Figure 16 is shot point local distribution schematic diagram when carrying out orthogonal laying according to the employing of the embodiment of the present invention identical shot point density and acceptance point density;
Figure 17 carries out receiving according to the acceptance point of the employing same size of the embodiment of the present invention and adopts 25 meters × 25 rice and flour units to carry out analyzing the acceptance point of the orthogonal laying equal densities obtained and shot point and active arranged distribution schematic diagram;
Figure 18 is shot point according to the embodiment of the present invention and the orthogonal laying comparison diagram of acceptance point;
Figure 19 is the structured flowchart of the apparatus for establishing of earthquake-capturing recording geometry according to the embodiment of the present invention.
Embodiment
For making the object, technical solutions and advantages of the present invention clearly understand, below in conjunction with embodiment and accompanying drawing, the present invention is described in further details.At this, exemplary embodiment of the present invention and illustrating for explaining the present invention, but not as a limitation of the invention.
Inventor find in earthquake-capturing recording geometry the geophone offset in each orientation main cause pockety be some orientation distributes comparatively sparse, such as, what line direction generally can show is comparatively sparse, suppose that survey line is longitudinally 0 °, survey line is laterally 90 °, so on the position angle of 0 °, 90 °, 180 ° and 270 ° four direction, geophone offset distribution is poor, in rarefaction state.Certainly, line direction neither be only had to present rarefaction state, and other orientation also there will be sparse situation.
Based on this, since inventor expects that some orientation is sparse, in these orientation, so increase big gun examine the impact that just can reduce sparse distribution to a certain extent.For this reason, in inventive embodiments, provide a kind of method for building up of earthquake-capturing recording geometry, as shown in Figure 1, comprise the following steps:
Step 101: the first recording geometry obtaining the orthogonal laying in the work area to be measured set up in advance;
Step 102: according to the bin attributes of described first recording geometry, determines that geophone offset distributes sparse one or more orientation;
Step 103: increase big gun and examine in described one or more orientation, and to increase big gun examine after the first recording geometry in all shot points blow out according to the rectangular arranged size of described first recording geometry, obtain the second recording geometry.
In above-described embodiment mode, by the sparse orientation that distributes in the recording geometry of determining built orthogonal laying, the orientation that the distribution determined is sparse increases big gun examine, and using increase big gun examine after recording geometry as final recording geometry, thus efficiently solve the problem that in prior art, in recording geometry, some orientation geophone offset is sparse, reach and make each orientation in recording geometry relatively uniform, there is not the technique effect in abnormal sparse orientation.
Generally the recording geometry of orthogonal laying sets up according to the zone of interest position in current work area, the zone of interest degree of depth, resolution, the follow-up analysis requirement etc. to oily field and reservoir, that is, be need to set up respectively corresponding with this work area to establish recording geometry for different work areas.
In above-mentioned steps 103, consider that sparse degree has relation with some distance and line-spacing, therefore, big gun can be increased in such a way examine in the one or more orientation determined:
S1: determine radius of a circle;
Particularly, radius of a circle can be determined in such a way: the acceptance point distance and the shot point distance that obtain the first recording geometry (namely carrying out the orthogonal laying recording geometry before sparse filling); The size of acceptance point distance and shot point distance is compared, when acceptance point is apart from being greater than shot point distance, using the N of shot point distance doubly as radius; When acceptance point is apart from being less than shot point distance, using the M of acceptance point distance doubly as radius, namely using acceptance point distance and shot point apart from the multiple of smaller as radius, wherein, N and M can equal also can be unequal.
S2: following operation is performed to each shot point in the first recording geometry:
Centered by current shot point, with described radius for radius draws circle;
By the point of each orientation on the circle drawn in described one or more orientation, as the shot point increased;
S3: following operation is performed to each geophone station in the first recording geometry:
Centered by current geophone station, with described radius for radius draws circle;
By the point of each orientation on the circle drawn in described one or more orientation, as the geophone station increased.
That is, centered by each acceptance point (i.e. geophone station) of recording geometry and shot point (i.e. shot point), with the radius determined for radius draws circle, the point of each sparse orientation on circle is obtained.
Consider that the big gun increased in recording geometry is cautious more, so required cost is higher, for this reason, if need the cost of restriction, can by the point of each orientation in one or more orientation on the circle drawn, after the shot point increased: determine that the big gun increased is examined quantity and whether exceeded predetermined quantity; If exceeded, then carry out stress and strain model to the big gun of each orientation on the circle drawn in one or more orientation is cautious; Ask for the combination center point in each grid; Combination center point in each grid determined is examined as the big gun that will increase, the quantity that the big gun that effectively reducing in this way increases is examined, effectively can solve again the problem of sparse distribution.
Generally all there will be sparse on the vertical and horizontal of survey line, but sometimes may judge some orientation by accident when judgement, such as may think that the vertical and horizontal of survey line think not sparse orientation, to this, to be on the safe side, at every turn determining that geophone offset distributes behind sparse orientation, all judge the orientation whether comprising the vertical and horizontal place of survey line in these sparse orientation that distribute, if do not comprised, then using the orientation at the vertical and horizontal place of described survey line also as the sparse orientation of distribution, effectively can avoid the vertical and horizontal omitting survey line like this.
Such as: survey line is longitudinally 0 °, be laterally 90 °, the sparse orientation determined is: 45 °, 135 °, 225 ° and 315 °, so just using 0 °, 90 °, 180 ° and orientation, 270 ° of four direction also as sparse orientation.
Provide a specific embodiment in this example, but, it should be noted that this specific embodiment is only to better the present invention is described, not forming inappropriate limitation of the present invention.
In this example, propose a kind of three-dimensional seismic acquisition observation system method for building up, on the basis of orthogonal laying recording geometry, by increasing acceptance point and/or shot point in some orientation, the azimuthal distribution realizing CMP collection is more even, realizes being uniformly distributed of geophone offset in 360 degree of planar range, the bin generated by the method, geophone offset is evenly distributed, and in the bin that can overcome orthogonal recording geometry, sparse uneven problem appears in geophone offset distribution.
Particularly, can comprise the following steps:
Step 1: the recording geometry (G0) generating some distance, the satisfactory orthogonal laying of line-spacing according to work area geological tasks, and adopt G0 to lay shot point and the geophone station in work area.
Such as: the orthogonal parameter laying the recording geometry generated is:
Step 2: analyze the bin attributes of the recording geometry of orthogonal laying, determines that geophone offset distributes sparse orientation α
i, i=1,2 ... m, if sparse orientation α
ido not comprise survey line vertical and horizontal, then four horizontal and vertical orientation of survey line are also added as sparse orientation.
As shown in Figures 2 and 3, the azimuthal distribution of the recording geometry geophone offset of orthogonal laying is uneven, and wherein, Fig. 2 is the azimuthal distribution of the orthogonal arranging system geophone offset not limitting offset distance, and Fig. 3 is limited the azimuthal distribution of orthogonal arranging system geophone offset of offset distance.
As can be seen from Fig. 2 and 3: skewness is mainly manifested in 0 °, 90 °, 180 ° and 270 °, next shows 45 °, 135 °, 225 ° and 315 °.The position angle of geophone offset be employing 50 meters of acceptance points apart from, 50 meters of shot points apart from, 200 meters receive line-spacings and 200 meters of excitation line apart from orthogonal laying, adopt 40 line single line 160 roads to carry out single big gun and receive and generate bin and produce.
Step 3: select radius r, wherein, r can value be the k of smallest point distance doubly, wherein, smallest point apart from be generation orthogonal laying recording geometry acceptance point apart from or shot point apart from smaller, k is positive integer.
Step 4: with each acceptance point in recording geometry G0 and/or each shot point (i.e. shot point and geophone station) P
j(j=1,2 ... n-1, n) centered by, be that radius draws circle with r, obtain sparse orientation α
ipoint P on circle
ji.
As shown in Figure 4, be the circle formed, and the point in mapping on circle, wherein, the center of circle is that big gun is examined, and the small circle on circle represents the bearing point that formed in the sparse orientation of justifying.As seen from Figure 4, the point of sparse orientation and interpolation is all determined according to azimuthal distribution as a setting, and eight orientation are respectively 0 °, 45 °, 90 °, 135 °, 180 ° 225 °, 270 ° and 315 °.
Step 5: if the some P generated
jimeet cost requirement, then adopt the template length identical with recording geometry G0 and width to set up recording geometry, and simulation is blown out and analyzes, and determines shot point, geophone station position and recording geometry, otherwise enters step 6.
Step 6: by all acceptance points and/or sparse orientation α corresponding to each shot point
ipoint P on circle
jiaccording to the longitudinal size B of bin
inxand transverse dimension B
inycarry out stress and strain model, and obtain the some P in each grid
jicombination center point.
Step 7: combination center point acceptance point process formed is as new acceptance point, combination center point shot point process formed is as new shot point, the template length identical with recording geometry G0 and width is adopted to set up recording geometry, and simulation is blown out and analyzes, determine shot point, geophone station position and recording geometry, namely, with increase big gun examine after recording geometry in all shot points blow out according to the rectangular arranged size of the recording geometry of original orthogonal laying, obtain final earthquake-capturing recording geometry.
Be illustrated in figure 5 acceptance point and the shot point distribution of combination center formation, Fig. 6 is acceptance point and the shot point distribution that a distance equals line-spacing, carries out contrast as can be seen from Fig. 5 and Fig. 6: the acceptance point that combination center is formed, shot point density and point are apart from the acceptance point, the shot point equal density that equal line-spacing.
Be illustrated in figure 7 a distance and equal the acceptance point of line-spacing and the azimuthal distribution of shot point generation, Fig. 8 is the acceptance point of combination center formation and the azimuthal distribution of shot point generation, contrast is carried out: the two bin size adopted is identical, and degree of covering is identical as can be seen from Fig. 7 with Fig. 8.
As can be seen here, the orientation Sparse Problems that can overcome orthogonal laying recording geometry and exist examined by the big gun laid by the orthogonal laying of method establishment in this example recording geometry, to make in bin geophone offset distribution relatively more even, lay result and more level off to the laying result of random observation system.
Be described for a specific embodiment:
For the orthogonal laying recording geometry that parameter is following:
Its acceptance point and shot point distribution and active arranged distribution as shown in Figure 9, wherein, wherein comparatively weak color represents acceptance point and shot point distribution, and darker color represents active arranged distribution, and so-called activity distribution is after middle shot point is blown out has be in excited state for which region.
25 meters × 25 rice and flour units are adopted to analyze to the orthogonal laying recording geometry shown in Fig. 9, the conventional orthorhombic obtained lays the azimuthal distribution of geophone offset in bin as shown in Figure 10, as seen from Figure 10: compared with distributing with other orientation geophone offset, 0 °, orientation, 45 °, 90 °, 135 °, 180 °, 225 °, 270 °, 315 ° display geophone offsets distribute sparse.
The acceptance point lay above-mentioned conventional orthorhombic and shot point carry out the sparse compensation deals of geophone offset, with acceptance point and shot point for the center of circle, with 100 meters for radius, add 0 °, orientation, 45 °, 90 °, 135 °, 180 °, 225 °, 270 °, the compensation point of 315 °, acceptance point carpenters square cun and/or the gridding of shot point carpenters square little progress row can be adopted after whole acceptance point after being compensated and shot point, acceptance point local distribution figure as shown in figure 11 and shot point local distribution figure is as shown in figure 12 obtained after gridding, after gridding, whole shot point and acceptance point and active arrangement are as shown in figure 13.
Adopt the bin of 25 meters × 25 meters to analyze the orthogonal laying recording geometry after above-mentioned compensation and gridding, in the bin obtained, the azimuthal distribution of geophone offset as shown in figure 14.
As shown in figure 15 for adopting identical shot point density and acceptance point density to carry out orthogonal laying, acceptance point local distribution schematic diagram, Figure 16 carries out orthogonal laying for adopting identical shot point density and acceptance point density, shot point local distribution schematic diagram, as shown in figure 17 for adopting the acceptance point of same size to receive, and adopt 25 meters × 25 rice and flour units to analyze, the acceptance point of the orthogonal laying equal densities obtained and shot point and active arranged distribution schematic diagram.
As shown in figure 18, before the left side represents compensation, after intermediate representation compensation, the right represents equal densities shot point and the orthogonal laying of acceptance point, and by comparing Figure 18 and can finding, what provided by the application is set up orthogonal laying recording geometry by compensation way and can obtain good effect.
In this example, bin attributes analysis after laying according to orthogonal recording geometry, determines that geophone offset distributes sparse orientation, then with whole acceptance point and/or shot point be the center of circle, with the integral multiple of smallest point distance for radius, add sparse orientation compensation point, generate the acceptance point after compensating and/or shot point.And then, can require to determine whether to carry out lower step gridding according to cost of investment, wherein, gridding generates the acceptance point after compensating and shot point with acceptance point distance and/or shot point apart from carrying out gridding to the acceptance point after compensation and shot point, then carries out exciting and bin attributes analysis with the acceptance point of same size scope.
Based on same inventive concept, additionally provide a kind of apparatus for establishing of earthquake-capturing recording geometry in the embodiment of the present invention, as described in the following examples.The principle of dealing with problems due to the apparatus for establishing of earthquake-capturing recording geometry is similar to the method for building up of earthquake-capturing recording geometry, therefore the enforcement of the apparatus for establishing of earthquake-capturing recording geometry see the enforcement of the method for building up of earthquake-capturing recording geometry, can repeat part and repeats no more.Following used, term " unit " or " module " can realize the software of predetermined function and/or the combination of hardware.Although the device described by following examples preferably realizes with software, hardware, or the realization of the combination of software and hardware also may and conceived.Figure 19 is a kind of structured flowchart of the apparatus for establishing of the earthquake-capturing recording geometry of the embodiment of the present invention, can comprise as shown in figure 19: acquisition module 1901, orientation determination module 1902 and interpolation module 1903, be described this structure below.
Acquisition module 1901, for obtaining the first recording geometry of the orthogonal laying in the work area to be measured set up in advance;
Orientation determination module 1902, for the bin attributes according to described first recording geometry, determines that geophone offset distributes sparse one or more orientation;
Add module 1903, examine for increasing big gun in described one or more orientation, and to increase big gun examine after the first recording geometry in all shot points blow out according to the rectangular arranged size of described first recording geometry, obtain the second recording geometry.
In one embodiment, adding module 1903 can comprise: radius determining unit, for determining radius of a circle; Add shot point unit, for performing following operation to each shot point in described first recording geometry: centered by current shot point, with described radius for radius draws circle; By the point of each orientation on the circle drawn in described one or more orientation, as the shot point increased; Add geophone station unit, for performing following operation to each geophone station in described first recording geometry: centered by current geophone station, with described radius for radius draws circle; By the point of each orientation on the circle drawn in described one or more orientation, as the geophone station increased.
In one embodiment, radius determining unit can comprise: point is apart from determining subelement, for obtaining acceptance point distance and the shot point distance of described first recording geometry; Size compares subelement, for comparing the size of described acceptance point distance and described shot point distance; First determines subelement, and for when described acceptance point distance is greater than described shot point distance, using the N of described shot point distance doubly as radius, wherein, N is positive integer; Second determines subelement, and for when described acceptance point distance is less than described shot point distance, using the M of described acceptance point distance doubly as radius, wherein, M is positive integer.
In one embodiment, add shot point unit and for by the point of each orientation on the circle drawn in described one or more orientation, after the shot point increased, can also determine whether the shot point quantity increased exceeds predetermined quantity; If exceeded, then stress and strain model is carried out to the point of each orientation on the circle drawn in described one or more orientation; Ask for the combination center point in each grid; Using the combination center point in each grid of determining as the shot point that will increase; Add geophone station unit also for by the point of each orientation on the circle drawn in described one or more orientation, after the geophone station increased, determine whether the geophone station quantity increased exceeds predetermined quantity; If exceeded, then stress and strain model is carried out to the point of each orientation on the circle drawn in described one or more orientation; Ask for the combination center point in each grid; Using the combination center point in each grid of determining as the geophone station that will increase.
In one embodiment, orientation determination module 1902 specifically can also be used at the bin attributes according to described first recording geometry, determine that geophone offset distributes after sparse one or more orientation, determine the orientation whether comprising the vertical and horizontal place of survey line in one or more orientation that the distribution that obtains is sparse; If do not comprised, then using the orientation at the vertical and horizontal place of described survey line also as the sparse orientation of distribution.
In another embodiment, additionally provide a kind of software, this software is for performing the technical scheme described in above-described embodiment and preferred implementation.
In another embodiment, additionally provide a kind of storage medium, store above-mentioned software in this storage medium, this storage medium includes but not limited to: CD, floppy disk, hard disk, scratch pad memory etc.
From above description, can find out, the embodiment of the present invention achieves following technique effect: by the sparse orientation that distributes in the recording geometry of determining built orthogonal laying, the orientation that the distribution determined is sparse increases big gun examine, and using increase big gun examine after recording geometry as final earthquake-capturing recording geometry, thus efficiently solve the problem that in the recording geometry of orthogonal laying in prior art, some orientation geophone offset is sparse, reach and make each orientation in earthquake-capturing recording geometry relatively uniform, there is not the technique effect in abnormal sparse orientation.
Obviously, those skilled in the art should be understood that, each module of the above-mentioned embodiment of the present invention or each step can realize with general calculation element, they can concentrate on single calculation element, or be distributed on network that multiple calculation element forms, alternatively, they can realize with the executable program code of calculation element, thus, they can be stored and be performed by calculation element in the storage device, and in some cases, step shown or described by can performing with the order be different from herein, or they are made into each integrated circuit modules respectively, or the multiple module in them or step are made into single integrated circuit module to realize.Like this, the embodiment of the present invention is not restricted to any specific hardware and software combination.
The foregoing is only the preferred embodiments of the present invention, be not limited to the present invention, for a person skilled in the art, the embodiment of the present invention can have various modifications and variations.Within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (10)
1. a method for building up for earthquake-capturing recording geometry, is characterized in that, comprising:
Obtain the first recording geometry of the orthogonal laying in the work area to be measured set up in advance;
According to the bin attributes of described first recording geometry, determine that geophone offset distributes sparse one or more orientation;
Described one or more orientation increases big gun examine, and to increase big gun examine after the first recording geometry in all shot points blow out according to the rectangular arranged size of described first recording geometry, obtain the second recording geometry.
2. the method for claim 1, is characterized in that, described one or more orientation increases big gun and examines, comprising:
Determine radius of a circle;
Following operation is performed to each shot point in described first recording geometry:
Centered by current shot point, with described radius for radius draws circle;
By the point of each orientation on the circle drawn in described one or more orientation, as the shot point increased;
Following operation is performed to each geophone station in described first recording geometry:
Centered by current geophone station, with described radius for radius draws circle;
By the point of each orientation on the circle drawn in described one or more orientation, as the geophone station increased.
3. method as claimed in claim 2, is characterized in that, determine radius of a circle, comprising:
Obtain acceptance point distance and the shot point distance of described first recording geometry;
The size of described acceptance point distance and described shot point distance is compared;
When described acceptance point is apart from being greater than described shot point distance, using the N of described shot point distance doubly as radius, wherein, N is positive integer;
When described acceptance point is apart from being less than described shot point distance, using the M of described acceptance point distance doubly as radius, wherein, M is positive integer.
4. method as claimed in claim 2, is characterized in that:
By the point of each orientation on the circle drawn in described one or more orientation, after the shot point increased, described method also comprises:
Determine whether the shot point quantity increased exceeds predetermined quantity;
If exceeded, then stress and strain model is carried out to the point of each orientation on the circle drawn in described one or more orientation;
Ask for the combination center point in each grid;
Using the combination center point in each grid of determining as the shot point that will increase;
By the point of each orientation on the circle drawn in described one or more orientation, after the geophone station increased, described method also comprises:
Determine whether the geophone station quantity increased exceeds predetermined quantity;
If exceeded, then stress and strain model is carried out to the point of each orientation on the circle drawn in described one or more orientation;
Ask for the combination center point in each grid;
Using the combination center point in each grid of determining as the geophone station that will increase.
5. the method according to any one of Claims 1-4, is characterized in that, at the bin attributes according to described first recording geometry, determine that geophone offset distributes after sparse one or more orientation, described method also comprises:
Determine the orientation whether comprising the vertical and horizontal place of survey line in one or more orientation that the distribution that obtains is sparse;
If do not comprised, then using the orientation at the vertical and horizontal place of described survey line also as the sparse orientation of distribution.
6. an apparatus for establishing for earthquake-capturing recording geometry, is characterized in that, comprising:
Acquisition module, for obtaining the first recording geometry of the orthogonal laying in the work area to be measured set up in advance;
Orientation determination module, for the bin attributes according to described first recording geometry, determines that geophone offset distributes sparse one or more orientation;
Add module, examine for increasing big gun in described one or more orientation, and to increase big gun examine after the first recording geometry in all shot points blow out according to the rectangular arranged size of described first recording geometry, obtain the second recording geometry.
7. device as claimed in claim 6, it is characterized in that, described interpolation module comprises:
Radius determining unit, for determining radius of a circle;
Add shot point unit, for performing following operation to each shot point in described first recording geometry: centered by current shot point, with described radius for radius draws circle; By the point of each orientation on the circle drawn in described one or more orientation, as the shot point increased;
Add geophone station unit, for performing following operation to each geophone station in described first recording geometry: centered by current geophone station, with described radius for radius draws circle; By the point of each orientation on the circle drawn in described one or more orientation, as the geophone station increased.
8. device as claimed in claim 7, it is characterized in that, described radius determining unit, comprising:
Point is apart from determining subelement, for obtaining acceptance point distance and the shot point distance of described first recording geometry;
Size compares subelement, for comparing the size of described acceptance point distance and described shot point distance;
First determines subelement, and for when described acceptance point distance is greater than described shot point distance, using the N of described shot point distance doubly as radius, wherein, N is positive integer;
Second determines subelement, and for when described acceptance point distance is less than described shot point distance, using the M of described acceptance point distance doubly as radius, wherein, M is positive integer.
9. device as claimed in claim 7, is characterized in that:
Described interpolation shot point unit, also for by the point of each orientation on the circle drawn in described one or more orientation, after the shot point increased, determines whether the shot point quantity increased exceeds predetermined quantity; If exceeded, then stress and strain model is carried out to the point of each orientation on the circle drawn in described one or more orientation; Ask for the combination center point in each grid; Using the combination center point in each grid of determining as the shot point that will increase;
Described interpolation geophone station unit, also for by the point of each orientation on the circle drawn in described one or more orientation, after the geophone station increased, determines whether the geophone station quantity increased exceeds predetermined quantity; If exceeded, then stress and strain model is carried out to the point of each orientation on the circle drawn in described one or more orientation; Ask for the combination center point in each grid; Using the combination center point in each grid of determining as the geophone station that will increase.
10. the device according to any one of claim 6 to 9, it is characterized in that, described orientation determination module is concrete also at the bin attributes according to described first recording geometry, determine that geophone offset distributes after sparse one or more orientation, determine the orientation whether comprising the vertical and horizontal place of survey line in one or more orientation that the distribution that obtains is sparse; If do not comprised, then using the orientation at the vertical and horizontal place of described survey line also as the sparse orientation of distribution.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510671042.2A CN105388513B (en) | 2015-10-15 | 2015-10-15 | The method for building up and device of earthquake-capturing observation system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510671042.2A CN105388513B (en) | 2015-10-15 | 2015-10-15 | The method for building up and device of earthquake-capturing observation system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105388513A true CN105388513A (en) | 2016-03-09 |
| CN105388513B CN105388513B (en) | 2017-12-05 |
Family
ID=55420994
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510671042.2A Active CN105388513B (en) | 2015-10-15 | 2015-10-15 | The method for building up and device of earthquake-capturing observation system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105388513B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106019360A (en) * | 2016-05-19 | 2016-10-12 | 吉林大学 | Information driving optimization method for three-dimensional earthquake observation system |
| CN106970414A (en) * | 2017-05-15 | 2017-07-21 | 安徽省煤田地质局物探测量队 | A kind of 3-d seismic exploration Te Guan groups design system |
| CN110579799A (en) * | 2019-09-04 | 2019-12-17 | 中国科学院地质与地球物理研究所 | A method and system for seismic acquisition and observation with equal travel time intervals |
| CN111474574A (en) * | 2019-01-23 | 2020-07-31 | 中国石油天然气集团有限公司 | Compressed sensing-based generation method and device for seismic acquisition observation system |
| CN112394389A (en) * | 2019-08-15 | 2021-02-23 | 中国石油天然气集团有限公司 | Method and system for determining shot-geophone inspection point of two-dimensional seismic exploration curved line |
| CN112782758A (en) * | 2019-11-11 | 2021-05-11 | 中国石油天然气集团有限公司 | State determination method and device for sparse sampling observation system |
| CN113568033A (en) * | 2020-04-28 | 2021-10-29 | 中国石油天然气集团有限公司 | Design method and device of three-dimensional irregular sampling seismic acquisition observation system |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4476552A (en) * | 1972-02-24 | 1984-10-09 | Conoco Inc. | Geophysical prospecting methods |
| US5999487A (en) * | 1996-07-25 | 1999-12-07 | Elf Exploration Production | Seismic reflection data acquisition and processing method and device for prospecting in tectonically complex environments |
| CN102262240A (en) * | 2011-06-22 | 2011-11-30 | 中国石油集团川庆钻探工程有限公司 | Automatic layout method for full coverage area of three-dimensional observation system |
| CN102778689A (en) * | 2011-05-13 | 2012-11-14 | 中国石油化工股份有限公司 | Wide curved line seismic data underground reflection line building method |
| CN103454680A (en) * | 2013-08-27 | 2013-12-18 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for calculating vertical coverage times of Walk-away VSP observing system |
| CN102628962B (en) * | 2011-11-02 | 2014-12-17 | 中国科学院地质与地球物理研究所 | Seismic acquisition observation system for underground microseismic detection |
-
2015
- 2015-10-15 CN CN201510671042.2A patent/CN105388513B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4476552A (en) * | 1972-02-24 | 1984-10-09 | Conoco Inc. | Geophysical prospecting methods |
| US5999487A (en) * | 1996-07-25 | 1999-12-07 | Elf Exploration Production | Seismic reflection data acquisition and processing method and device for prospecting in tectonically complex environments |
| CN102778689A (en) * | 2011-05-13 | 2012-11-14 | 中国石油化工股份有限公司 | Wide curved line seismic data underground reflection line building method |
| CN102262240A (en) * | 2011-06-22 | 2011-11-30 | 中国石油集团川庆钻探工程有限公司 | Automatic layout method for full coverage area of three-dimensional observation system |
| CN102628962B (en) * | 2011-11-02 | 2014-12-17 | 中国科学院地质与地球物理研究所 | Seismic acquisition observation system for underground microseismic detection |
| CN103454680A (en) * | 2013-08-27 | 2013-12-18 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for calculating vertical coverage times of Walk-away VSP observing system |
Non-Patent Citations (2)
| Title |
|---|
| 卢湘鹏: "西部山前带观测系统设计方法", 《物探与化探》 * |
| 姜福豪 等: "《中国石油学会2015年物探技术研讨会论文集》", 30 April 2015 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106019360A (en) * | 2016-05-19 | 2016-10-12 | 吉林大学 | Information driving optimization method for three-dimensional earthquake observation system |
| CN106970414A (en) * | 2017-05-15 | 2017-07-21 | 安徽省煤田地质局物探测量队 | A kind of 3-d seismic exploration Te Guan groups design system |
| CN106970414B (en) * | 2017-05-15 | 2019-02-12 | 安徽省煤田地质局物探测量队 | A kind of 3-d seismic exploration Te Guan group designing system |
| CN111474574A (en) * | 2019-01-23 | 2020-07-31 | 中国石油天然气集团有限公司 | Compressed sensing-based generation method and device for seismic acquisition observation system |
| CN111474574B (en) * | 2019-01-23 | 2023-04-25 | 中国石油天然气集团有限公司 | Method and device for generating seismic acquisition observation system based on compressed sensing |
| CN112394389A (en) * | 2019-08-15 | 2021-02-23 | 中国石油天然气集团有限公司 | Method and system for determining shot-geophone inspection point of two-dimensional seismic exploration curved line |
| CN110579799A (en) * | 2019-09-04 | 2019-12-17 | 中国科学院地质与地球物理研究所 | A method and system for seismic acquisition and observation with equal travel time intervals |
| CN112782758A (en) * | 2019-11-11 | 2021-05-11 | 中国石油天然气集团有限公司 | State determination method and device for sparse sampling observation system |
| CN112782758B (en) * | 2019-11-11 | 2024-06-25 | 中国石油天然气集团有限公司 | State determination method and device for sparse sampling observation system |
| CN113568033A (en) * | 2020-04-28 | 2021-10-29 | 中国石油天然气集团有限公司 | Design method and device of three-dimensional irregular sampling seismic acquisition observation system |
| CN113568033B (en) * | 2020-04-28 | 2024-06-04 | 中国石油天然气集团有限公司 | Design method and device of three-dimensional irregular sampling seismic acquisition observation system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105388513B (en) | 2017-12-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105388513A (en) | Method and device for establishing earthquake acquisition and observation system | |
| US10467357B2 (en) | Geobody continuity in geological models based on multiple point statistics | |
| US10871590B2 (en) | Electromagnetic data inversion | |
| US9229129B2 (en) | Reservoir geobody calculation | |
| CN102057368B (en) | Distribution of properties in a 3D volumetric model using a maximum continuity field | |
| CN102243321B (en) | Method and system for processing seismic pre-stack time migration | |
| US20150100293A1 (en) | Method for efficient dynamic gridding | |
| Seers et al. | Probabilistic constraints on structural lineament best fit plane precision obtained through numerical analysis | |
| CN104155703A (en) | Method and device for evaluating three-dimensional observing system | |
| CN111474574B (en) | Method and device for generating seismic acquisition observation system based on compressed sensing | |
| Jahandari et al. | Forward modeling of direct-current resistivity data on unstructured grids using an adaptive mimetic finite-difference method | |
| Vachiratienchai et al. | A hybrid finite difference–finite element method to incorporate topography for 2D direct current (DC) resistivity modeling | |
| US20150317419A1 (en) | Local direct sampling method for conditioning an existing reservoir model | |
| CN118171586B (en) | Multi-beam survey line layout method, system, equipment and storage medium | |
| CN112379413B (en) | Irregular seismic source characterization method and device based on energy spectrum equivalence | |
| JP5435620B2 (en) | Sensitivity-compatible inverse analysis method in numerical analysis of resistivity method | |
| Loke et al. | The use of a semi‐structured finite‐element mesh in 3‐D resistivity inversion | |
| CN111123361B (en) | Method and device for regularly reconstructing seismic data of vertical seismic profile and storage medium | |
| CN113568033A (en) | Design method and device of three-dimensional irregular sampling seismic acquisition observation system | |
| Smaglichenko | Modification of Gaussian elimination for the complex system of seismic observations | |
| CN116089784A (en) | Data mining implementation method based on cloud computing | |
| CN111965694B (en) | Method and device for determining position of physical point of earthquake and earthquake observation system | |
| EP2904434B1 (en) | A method for determining a karstic region | |
| Pinnaduwa et al. | Three dimensional stochastic joint geometry modeling including a verification: A case study | |
| Kulatilake et al. | Stochastic three dimensional joint geometry modelling including a verification to an area in Stripa Mine, Sweden |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |