CN103605163B - A kind of static correction value method for quantitatively evaluating and device - Google Patents
A kind of static correction value method for quantitatively evaluating and device Download PDFInfo
- Publication number
- CN103605163B CN103605163B CN201310516700.1A CN201310516700A CN103605163B CN 103605163 B CN103605163 B CN 103605163B CN 201310516700 A CN201310516700 A CN 201310516700A CN 103605163 B CN103605163 B CN 103605163B
- Authority
- CN
- China
- Prior art keywords
- arrival
- grid
- static correction
- correction value
- matching
- 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.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a kind of static correction value method for quantitatively evaluating and device, wherein, the method comprises: select the first arrival in described grid, the first arrival of selecting is carried out to pretreatment; Wherein, described pretreated method is: the static correction value of the static correction value of shot point corresponding described first arrival, acceptance point and first arrival time are added; Pretreated first arrival in each grid is carried out to matching and obtain matching first arrival accordingly; Obtain the time difference root-mean-square value of each grid according to described matching first arrival and pretreated first arrival; The time difference root-mean-square value of cumulative all grids obtains the quantitative error of the static correction value of survey area; According to the quantitative error of the static correction value of described survey area, static correction value is evaluated.
Description
Technical field
The present invention relates to geophysical prospecting for oil static correction field, particularly a kind of first arrival time of utilizing seismic waveTo static correction value method for quantitatively evaluating and device.
Background technology
The computational methods of static correction value are a lot, have interlayer relation Y-factor method Y, refraction process, Relative refraction statics,Chromatography etc., various static correcting methods have its scope of application and application conditions separately. Practically in seismic data processing,By multiple static correcting method deriving static correction values, sometimes even need comprehensive multiple static correction side in complicated earth surface areaMethod obtains the static correction value of a set of optimum, and effect how to evaluate a set of static correction value is a very important job.
As shown in Figure 1, work area g does not add static correction profile; As shown in Figure 2, work area g elevation static correction stack is cutd openFace figure; As shown in Figure 3, work area g chromatography static correction stacked profile map; As shown in Figure 4, work area gGMG refractionStatic correction stacked profile map; As shown in Figure 5, work area gKL refraction static correction stacked profile map. Current evaluation sideMethod is all to be done qualitatively and to be evaluated by personnel's observation, by static correction value is applied to various roads collection (Gong Bao road collection,Gong Jian road collection, common midpoint gather, common offset road collection) the upper effect of evening up that contrasts first arrival, or at stacked sectionThe quality of static correction value is evaluated in upper these measures of contrast image quality. Increase like this workload for the treatment of people, andAnd sometimes this evaluate qualitatively can not objective the quality that reflects static correction value, especially for static correction valueLittle or the lower situation of data signal to noise ratio of difference, for example Fig. 4 and Fig. 5 are just difficult to find out quiet school from stacked sectionThe quality of plus effect, so industry needs a kind of method to do quantitative assessment to static correction value, this is also static correction technologyThe demand of development.
Summary of the invention
For addressing the above problem, the little or lower situation of data signal to noise ratio for static correction value difference, the present invention proposesA kind of static correction value method for quantitatively evaluating and device, to overcome the static correction effect to above-mentioned two situations in prior artCannot evaluate accurately.
For achieving the above object, the invention provides a kind of static correction value method for quantitatively evaluating, the method comprises:
Select the first arrival in grid, the first arrival of selecting is carried out to pretreatment; Wherein, described pretreated method is:The static correction value of the static correction value of shot point corresponding described first arrival, acceptance point and first arrival time are added;
Pretreated first arrival in each grid is carried out to matching and obtain matching first arrival accordingly;
Obtain the time difference root-mean-square value of each grid according to described matching first arrival and pretreated first arrival;
The time difference root-mean-square value of cumulative all grids, and accumulated value is averaging to the determining of static correction value that obtains survey areaAmount error;
According to the quantitative error of the static correction value of described survey area, static correction value is evaluated.
Optionally, in an embodiment of the present invention, the step of the first arrival in the described grid of described selection comprises:
Obtain shot point coordinate, acceptance point coordinate and the first arrival time of every one first arrival in survey area;
Described survey area is arranged to grid distance, and grid division, according to shot point coordinate, the reception of described every one first arrivalPoint coordinates obtains the middle point coordinates of described every one first arrival, according to the first arrival time of described every one first arrival, middle point coordinatesObtain the first arrival in corresponding grid;
From first arrival corresponding to each grid, select the first arrival within the scope of suitable geophone offset.
Optionally, in an embodiment of the present invention, describedly from first arrival corresponding to each grid, select suitable geophone offset modelThe method of the first arrival in enclosing is refraction top and bottom process.
Optionally, in an embodiment of the present invention, described matching adopts fitting a straight line or multinomialMatching.
For achieving the above object, the present invention also provides a kind of static correction value quantitative evaluation device, and this device comprises:
Pretreatment unit, for selecting the first arrival of grid, carries out pretreatment by the first arrival of selecting; Wherein, described inPretreated method is: by the static correction value of the static correction value of shot point corresponding described first arrival, acceptance point and first arrival timeBe added;
Matching unit, for to each grid pretreated first arrival carry out matching and obtain matching first arrival accordingly;
Time difference root-mean-square value acquiring unit, for obtaining each grid according to described matching first arrival and pretreated first arrivalTime difference root-mean-square value;
The quantitative error acquiring unit of static correction value, for the time difference root-mean-square value of cumulative all grids, and to accumulated valueBe averaging the quantitative error of the static correction value that obtains survey area;
Evaluation unit, for evaluating static correction value according to the quantitative error of the static correction value of described survey area.
Optionally, in an embodiment of the present invention, described pretreatment unit comprises:
First arrival parameter acquisition module, for obtaining shot point coordinate, the acceptance point coordinate of every one first arrival in survey areaAnd first arrival time;
Grid first arrival acquisition module, for described survey area is arranged to grid distance, grid division, according to described eachThe shot point coordinate of road first arrival, the middle point coordinates that acceptance point coordinate obtains described every one first arrival, at the beginning of described every oneFirst arrival time, middle point coordinates extremely obtain the first arrival in corresponding grid;
Select module, for selecting the first arrival within the scope of suitable geophone offset from first arrival corresponding to each grid.
Optionally, in an embodiment of the present invention, the method that described selection module adopts is refraction top and bottom process.
Optionally, in an embodiment of the present invention, described matching unit adopts fitting a straight line or fitting of a polynomial.
Technique scheme has following beneficial effect: the technical program can be evaluated the quality of static correction value rapidly,Have more intuitive advantage with respect to traditional qualitative static correction evaluation method, this algorithm has reduced static correction greatlyThe workload of amount contrast.
Brief description of the drawings
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existingHave the accompanying drawing of required use in technical description to be briefly described, apparently, the accompanying drawing in the following describes onlySome embodiments of the present invention, for those of ordinary skill in the art, in the prerequisite of not paying creative workUnder, can also obtain according to these accompanying drawings other accompanying drawing.
Fig. 1 is that work area g does not add static correction profile;
Fig. 2 is work area g elevation static correction stacked profile map;
Fig. 3 is work area g chromatography static correction stacked profile map;
Fig. 4 is work area gGMG refraction static correction stacked profile map;
Fig. 5 is work area gKL refraction static correction stacked profile map;
Fig. 6 is a kind of static correction value method for quantitatively evaluating flow chart that the present invention proposes;
Fig. 7 is a kind of static correction value quantitative evaluation device block diagram that the present invention proposes;
Fig. 8 is the cautious elevation curve map of work area g big gun in embodiment;
Fig. 9 is bin number 101 effective geophone offset scope and fit line schematic diagrames in embodiment;
Figure 10 is the various static correcting method root-mean-square error of work area g block diagram in embodiment.
Detailed description of the invention
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, completeDescribe wholely. Obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiment.Based on the embodiment in the present invention, those of ordinary skill in the art obtain under creative work prerequisite not makingEvery other embodiment, belongs to the scope of protection of the invention.
As shown in Figure 6, a kind of static correction value method for quantitatively evaluating flow chart proposing for the present invention. The method comprises:
Step 601): select the first arrival in grid, the first arrival of selecting is carried out to pretreatment; Wherein, described pre-placeThe method of reason is: the static correction value of the static correction value of shot point corresponding described first arrival, acceptance point and first arrival time are added;
Step 602): pretreated first arrival in each grid is carried out to matching and obtain matching first arrival accordingly;
Step 603): the time difference root-mean-square value that obtains each grid according to described matching first arrival and pretreated first arrival;
Step 604): the time difference root-mean-square value of cumulative all grids, and accumulated value is averaging and obtains survey areaThe quantitative error of static correction value;
Step 605): according to the quantitative error of the static correction value of described survey area, static correction value is evaluated.
Optionally, in an embodiment of the present invention, the step of the first arrival in the described grid of described selection comprises:
Obtain shot point coordinate, acceptance point coordinate and the first arrival time of every one first arrival in survey area;
Described survey area is arranged to grid distance, and grid division, according to shot point coordinate, the reception of described every one first arrivalPoint coordinates obtains the middle point coordinates of described every one first arrival, according to the first arrival time of described every one first arrival, middle point coordinatesObtain the first arrival in corresponding grid;
From first arrival corresponding to each grid, select the first arrival within the scope of suitable geophone offset.
Optionally, in an embodiment of the present invention, describedly from first arrival corresponding to each grid, select suitable geophone offset modelThe method of the first arrival in enclosing is refraction top and bottom process.
Optionally, in an embodiment of the present invention, described matching adopts fitting a straight line or multinomialMatching.
As shown in Figure 7, a kind of static correction value quantitative evaluation device block diagram proposing for the present invention. This device comprises:
Pretreatment unit 701, for selecting the first arrival of grid, carries out pretreatment by the first arrival of selecting; Wherein,Described pretreated method is: by the static correction value of the static correction value of shot point corresponding described first arrival, acceptance point and first arrivalTime is added;
Matching unit 702, at the beginning of to each grid, pretreated first arrival is carried out matching and is obtained matching accordinglyExtremely;
Time difference root-mean-square value acquiring unit 703, for obtaining each according to described matching first arrival and pretreated first arrivalThe time difference root-mean-square value of grid;
The quantitative error acquiring unit 704 of static correction value, for the time difference root-mean-square value of cumulative all grids, and to tiredThe value added quantitative error that is averaging the static correction value that obtains survey area;
Evaluation unit 705, for commenting static correction value according to the quantitative error of the static correction value of described survey areaValency.
Optionally, in an embodiment of the present invention, described pretreatment unit 701 comprises:
First arrival parameter acquisition module 7011, for obtaining shot point coordinate, the reception of every one first arrival in survey areaPoint coordinates and first arrival time;
Grid first arrival acquisition module 7012, for described survey area is arranged to grid distance, grid division, according to instituteState the shot point coordinate of every one first arrival, the middle point coordinates that acceptance point coordinate obtains described every one first arrival, according to described everyFirst arrival time, the middle point coordinates of one first arrival obtain the first arrival in corresponding grid;
Select module 7013, for selecting the first arrival within the scope of suitable geophone offset from first arrival corresponding to each grid.
Optionally, in an embodiment of the present invention, the method that described selection module 7013 adopts is refraction top and bottom process.
Optionally, in an embodiment of the present invention, described matching unit 702 adopts fitting a straight line or fitting of a polynomial.Embodiment:
As shown in Figure 8, be the cautious elevation curve map of work area, mountain region g two dimension survey line big gun. Elevation changes at 210m～1500mBetween, observation system parameter: track pitch 20m, shooting distance 80m, totally 609 big guns, it is fixed that the present embodiment calculates static correctionThe step of amount error is as follows:
1) adopt usual way to obtain shot point and acceptance point coordinate, adopt common geological data pick-up method to obtainFirst arrival time;
2) grid is set apart from 80m, grid division, according to the middle point coordinates of per pass first arrival shot point and acceptance point, calculatesGo out the first arrival in each grid.
3) adopt common refraction layered approach, select the first arrival within the scope of suitable geophone offset.
4) geophone offset sequence is pressed in the first arrival in each grid, static correction value is applied to just supreme.
5) select a grid, grid number is 101(Gi, i=101), according to the actual conditions of data to grid GiIn within the scope of suitable geophone offset first arrival carry out fitting a straight line, as shown in Figure 9, calculate in geophone offset per pass first arrival andThe time difference of fit line is designated as Δ tj, calculates its root mean square and is designated as ti.
Wherein, i is grid sequence number, and j is the first arrival sequence number that participates in calculating in this grid, and n participates in grid calculatingFirst arrival number.
6) root mean square of cumulative all grids, asks its average to be designated as T, as the quantitative error of this cover static correction value.
Wherein, i is grid sequence number, and m is the meshes number that participates in calculating.
Fig. 1～Fig. 5 compares and can find out, Fig. 4 and Fig. 5 section imaging effect are best, be better than above three cover quietCorrecting value. But these two only there is nuance, be difficult to distinguish intuitively, illustrate GMG refraction statics andThe overall imaging of KL refraction statics is on close level.
As shown in figure 10, the quantitative error block diagram of the various static correcting methods of work area g calculating for this method. In figureThe quantitative error minimum of KL refraction static correction, illustrates the best results of KL refraction static correction.
Finally it should be noted that: above-mentioned only in order to the present invention to be described and unrestricted technical scheme described in the invention;Although this description has been described in detail to the present invention,, those skilled in the art still can be to thisInvention is modified or is equal to replacement, and all do not depart from technical scheme and the improvement thereof of the spirit and scope of the present invention, itsAll should be encompassed in claim scope of the present invention.
Claims (8)
1. a method that obtains static correction value quantitative error, is characterized in that, the method comprises:
Select the first arrival in grid, the first arrival of selecting is carried out to pretreatment; Wherein, described pretreated method is:
The static correction value of the static correction value of shot point corresponding described first arrival, acceptance point and first arrival time are added;
Pretreated first arrival in each grid is carried out to matching and obtain matching first arrival accordingly;
Obtain the time difference root-mean-square value of each grid according to described matching first arrival and pretreated first arrival;
The time difference root-mean-square value of cumulative all grids, and accumulated value is averaging to the determining of static correction value that obtains survey areaAmount error.
2. the method for claim 1, is characterized in that, the step of the first arrival in the described grid of described selectionComprise:
Obtain shot point coordinate, acceptance point coordinate and the first arrival time of every one first arrival in survey area;
Described survey area is arranged to grid distance, and grid division, according to shot point coordinate, the reception of described every one first arrivalPoint coordinates obtains the middle point coordinates of described every one first arrival, according to the first arrival time of described every one first arrival, middle point coordinatesObtain the first arrival in corresponding grid;
From first arrival corresponding to each grid, select the first arrival within the scope of suitable geophone offset.
3. method as claimed in claim 2, is characterized in that, describedly from first arrival corresponding to each grid, selectsThe method of the first arrival within the scope of suitable geophone offset is refraction top and bottom process.
4. the method for claim 1, is characterized in that, described matching adopts fitting a straight line or multinomialMatching.
5. a device that obtains static correction value quantitative error, is characterized in that, this device comprises:
Pretreatment unit, for selecting the first arrival of grid, carries out pretreatment by the first arrival of selecting; Wherein, described inPretreated method is: by the static correction value of the static correction value of shot point corresponding described first arrival, acceptance point and first arrival timeBe added;
Matching unit, for to each grid pretreated first arrival carry out matching and obtain matching first arrival accordingly;
Time difference root-mean-square value acquiring unit, for obtaining each grid according to described matching first arrival and pretreated first arrivalTime difference root-mean-square value;
The quantitative error acquiring unit of static correction value, for the time difference root-mean-square value of cumulative all grids, and to accumulated valueBe averaging the quantitative error of the static correction value that obtains survey area.
6. device as claimed in claim 5, is characterized in that, described pretreatment unit comprises:
First arrival parameter acquisition module, for obtaining shot point coordinate, the acceptance point coordinate of every one first arrival in survey areaAnd first arrival time;
Grid first arrival acquisition module, for described survey area is arranged to grid distance, grid division, according to described eachThe shot point coordinate of road first arrival, the middle point coordinates that acceptance point coordinate obtains described every one first arrival, at the beginning of described every oneFirst arrival time, middle point coordinates extremely obtain the first arrival in corresponding grid;
Select module, for selecting the first arrival within the scope of suitable geophone offset from first arrival corresponding to each grid.
7. device as claimed in claim 6, is characterized in that, the method that described selection module adopts is refraction pointLayer method.
8. device as claimed in claim 5, is characterized in that, described matching unit adopts fitting a straight line or multinomialFormula matching.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310516700.1A CN103605163B (en) | 2013-10-28 | 2013-10-28 | A kind of static correction value method for quantitatively evaluating and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310516700.1A CN103605163B (en) | 2013-10-28 | 2013-10-28 | A kind of static correction value method for quantitatively evaluating and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103605163A CN103605163A (en) | 2014-02-26 |
| CN103605163B true CN103605163B (en) | 2016-05-11 |
Family
ID=50123403
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310516700.1A Active CN103605163B (en) | 2013-10-28 | 2013-10-28 | A kind of static correction value method for quantitatively evaluating and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103605163B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112099091B (en) * | 2019-06-17 | 2024-03-01 | 中国石油天然气集团有限公司 | Static correction method and device for seismic exploration point data |
| CN112053030A (en) * | 2020-07-28 | 2020-12-08 | 中国石油天然气集团有限公司 | Method and device for evaluating positioning accuracy of detector |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4069471A (en) * | 1976-03-12 | 1978-01-17 | Geophysical Systems Corporation | Method of determining weathering corrections in seismic record processing |
| CN101363916A (en) * | 2007-08-10 | 2009-02-11 | 中国石油集团东方地球物理勘探有限责任公司 | High precision refraction static correction data inversion method |
| CN102692648A (en) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | Refraction wave residual static correction method based on cannon first arrival |
-
2013
- 2013-10-28 CN CN201310516700.1A patent/CN103605163B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4069471A (en) * | 1976-03-12 | 1978-01-17 | Geophysical Systems Corporation | Method of determining weathering corrections in seismic record processing |
| CN101363916A (en) * | 2007-08-10 | 2009-02-11 | 中国石油集团东方地球物理勘探有限责任公司 | High precision refraction static correction data inversion method |
| CN102692648A (en) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | Refraction wave residual static correction method based on cannon first arrival |
Non-Patent Citations (1)
| Title |
|---|
| 静校正量的评价思路;Adnan Al-maghouth 等;《石油物探》;20100930;第49卷(第5期);483-488 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103605163A (en) | 2014-02-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104142514A (en) | Three-dimensional earthquake observing system quantitative designing method | |
| Perez-Astudillo et al. | DNI, GHI and DHI ground measurements in Doha, Qatar | |
| CN104483705B (en) | A kind of three-dimensional residual static corrections | |
| GB2586403A (en) | Method, system, and apparatus for correcting translucency artifacts in data representing a support structure | |
| CN105224908A (en) | A kind of roadmarking acquisition method based on orthogonal projection and device | |
| CN102279973A (en) | Sea-sky-line detection method based on high gradient key points | |
| CN103605163B (en) | A kind of static correction value method for quantitatively evaluating and device | |
| GB2527247A (en) | Method and system for detecting heartbeat irregularities | |
| CN103576213A (en) | Method for determining uniformity of properties of shot-geophone distances | |
| CN103713318A (en) | Three-dimensional stacked seismic data splicing method | |
| CN103217710A (en) | Method and system for acquiring common offset distance and common azimuth angle domain gathers | |
| CN104091329A (en) | Calibration method and device for CT image as well as CT system | |
| CN104769459A (en) | High-efficient velocity analysis method for pre-stack time migration | |
| CN103576197A (en) | Method for extracting converted wave angle channel set | |
| CN105319576A (en) | Concave polygon region observation system automatic full coverage quick layout method | |
| CN106054252B (en) | A kind of method and device of pre-stack time migration | |
| CN103744116A (en) | Prestack gather all-time-domain phase-equalization correction method | |
| CN102692651A (en) | Preliminary wave residual static correction method with space-variant velocity | |
| CN104502963B (en) | Detection method of amplitude energy abnormality phenomenon of post-stack data | |
| CN103645501A (en) | A method for determining a largest offset | |
| CN106802962B (en) | A method of calculating and correct ocean essential mean value | |
| CN104730217B (en) | The defect distribution display methods of a kind of glass substrate and display device | |
| CN104502962B (en) | A kind of method and device of design shot point | |
| CN104121882B (en) | Be in a pout detection method and the device of steel pipe seam | |
| CN111047083A (en) | Processing method of stepping type vorticity covariance observation data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |