CN102254321B - Method for automatically identifying polarity-reversed channel based on first arrival wave - Google Patents
Method for automatically identifying polarity-reversed channel based on first arrival wave Download PDFInfo
- Publication number
- CN102254321B CN102254321B CN201110130303.1A CN201110130303A CN102254321B CN 102254321 B CN102254321 B CN 102254321B CN 201110130303 A CN201110130303 A CN 201110130303A CN 102254321 B CN102254321 B CN 102254321B
- Authority
- CN
- China
- Prior art keywords
- frontier point
- seismic
- seismic trace
- arrival
- polarity
- 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.)
- Expired - Fee Related
Links
Images
Abstract
The invention provides a method for automatically identifying a polarity-reversed channel based on a first arrival wave. The method comprises the following steps of: converting a seismic record oscillogram of a plurality of seismic channels into a grayscale, and performing edge detection on the grayscale to obtain boundary points of each seismic channel; determining uncertain boundary points in the obtained boundary points and deleting the uncertain boundary points; performing interpolation on the rest of boundary points of the seismic channels to obtain the boundary points of all the seismic channels; mapping the smoothed boundary points to the seismic record oscillogram so as to pick up a first takeoff of each seismic channel; and judging the polarity of the seismic channels according to whether a first arrival peak is behind the first arrival takeoff of the seismic channels.
Description
Technical field
The present invention relates to geological data process field, more particularly, relate to the method based on primary wave identification reversed polarity road.
Background technology
In petroleum prospecting, extremely important for the evaluation of field monitor record quality, reversed polarity often occurs, therefore, the evaluation of reversed polarity and analysis is seemed to particularly important in single big gun data.
Yet, the identification that is confined to naked eyes of current recognition method more, so workload is large, recognition efficiency is low, needs high efficiency recognition methods.
Summary of the invention
In order to improve the recognition efficiency of the seismic trace in seismologic record data, the invention provides a kind of method based on primary wave identification reversed polarity road.
A kind of method based on primary wave identification reversed polarity road according to the present invention comprises step: the seismic wave figure of a plurality of seismic traces is converted into gray-scale map and gray-scale map is carried out to the frontier point that rim detection obtains each seismic trace; Determine the incorrect frontier point in the frontier point obtaining and delete; Frontier point by remaining seismic trace carries out the frontier point that interpolation obtains all seismic traces; The frontier point of smoothing processing is mapped to the first arrival take-off that seismic wave figure picks up each seismic trace; And, according to the first arrival take-off of seismic trace, whether be the polarity that first arrival crest judges seismic trace afterwards.
In addition, can carry out conic fitting to the frontier point obtaining, the distance that departs from fit line be surpassed to predetermined value frontier point and be defined as incorrect frontier point and delete.
In addition, determine that the incorrect frontier point in frontier point the step of deleting can comprise: (a) to the frontier point obtaining by rim detection, judge the seismic trace before each seismic trace and this seismic trace frontier point time difference absolute value and with this seismic trace after the absolute value of difference of time of frontier point of seismic trace at least one whether to be greater than Preset Time poor, if it is poor to be greater than Preset Time, the frontier point of this seismic trace is confirmed as incorrect frontier point deleted; (b) the remaining frontier point to passing through step (a), calculate the mean value of absolute value of difference of time of the frontier point of all adjacent seismic traces, judge the seismic trace before each seismic trace and this seismic trace frontier point time difference absolute value and with this seismic trace after the absolute value of difference of time of frontier point of seismic trace at least one whether be greater than 4 times of described mean value, if be greater than 4 times, the frontier point of this seismic trace is confirmed as incorrect frontier point and deletes; (c) repeating step (b), until the difference absolute value of the time of the frontier point of the seismic trace before remaining each seismic trace and this seismic trace and be all less than 1.5 times of described mean value with the absolute value of difference of time of frontier point of seismic trace after this seismic trace.
In addition, by judging that the preset distance place, first arrival jump space of each seismic trace is the polarity that first arrival trough or first arrival crest are determined this seismic trace.
According to the present invention, utilize Image Edge-Detection to carry out automatic decision to the first arrival of single big gun data, on the basis of automatic decision, recycling polarity amplitude characteristic is judged the polarity of seismic trace, thereby has improved the efficiency of the polarity of judgement seismic trace.
Accompanying drawing explanation
By the description of embodiment being carried out below in conjunction with accompanying drawing, these and/or other aspect of the present invention and advantage will become clear and be easier to and understand, wherein:
Fig. 1 illustrates the process flow diagram of the method based on primary wave identification reversed polarity road according to an exemplary embodiment of the present invention.
Fig. 2 illustrates the definite incorrect frontier point of another exemplary embodiment according to the present invention and the process flow diagram of the step of deleting.
Fig. 3 illustrates to use according to of the present invention and automatically identifies based on primary wave the reversed polarity road that the method automatic decision in reversed polarity road goes out.
Embodiment
Now the embodiment of the present invention is described in detail, in the accompanying drawings, wherein, identical label represents same parts to its example shown all the time.Below with reference to the accompanying drawings embodiment is described to explain the present invention.
In geological data process field, primary wave typically refers to and is excited by shot point, arrives at first the seismic signal of acceptance point.The source of primary wave is very extensive, has almost contained the seismic event of all kinds, comprises direct wave, reflection wave, refraction wave, diffracted wave, diving Wave etc.Due in arriving the seismic signal of same point, the travel-time of primary wave is the shortest, and energy is conventionally stronger.
Below, with reference to Fig. 1, describe the method based on primary wave identification reversed polarity road according to an exemplary embodiment of the present invention in detail.Fig. 1 illustrates the process flow diagram of the method based on primary wave identification reversed polarity road according to an exemplary embodiment of the present invention.
At step S101, the seismic wave figure of a plurality of seismic traces is converted into gray-scale map and gray-scale map is carried out to the frontier point that rim detection obtains each seismic trace.
Primary wave shows as the boundary of effective reception signal in seismologic record, and energy is strong, take-off is obvious, and in seismic trace set, the first arrival take-off time between each seismic trace is conventionally more continuous, wherein, and the beginning of first arrival take-off indication primary wave signal.Therefore, can regard the line of the time of the first arrival take-off of each seismic trace in seismic trace set as border between disturbance and seismologic record data.Therefore, can extract frontier point by rim detection, that is, and the first arrival take-off of primary wave.
Because the frontier point that uses edge detection method to obtain image is common practise, at this, omit the description to it.
In step 102, determine the incorrect frontier point in the frontier point obtaining and delete.
Owing to there being the reasons such as error in edge detection algorithm, need to determine whether the frontier point of each seismic trace obtaining is correct frontier point.Due to the Time Continuous of the first arrival take-off of the primary wave between each seismic trace, therefore, the frontier point of acquisition should be also continuous in time.So frontier point that can be by judging seismic trace whether with its near the frontier point of seismic trace whether continuous in time, determine whether the frontier point of this seismic trace is correct frontier point.If not correct frontier point, delete the frontier point of this seismic trace and by interpolation method, obtain the frontier point of this seismic trace.
According to an exemplary embodiment of the present, the frontier point obtaining in step 101 is carried out to conic fitting, judge whether the frontier point of seismic trace and the distance of matched curve are greater than predetermined value, if the predetermined value of being greater than, is defined as the frontier point of this seismic trace incorrect frontier point and deletes.
Another exemplary embodiment according to the present invention, by the continuity of the judgement time of frontier point of seismic trace and near the time of the frontier point of the seismic trace this seismic trace, determines whether this seismic trace is correct frontier point.The time of frontier point refers to the time that seismic trace detects the first arrival take-off of primary wave.
Below, with reference to Fig. 2, describe the definite incorrect frontier point of another exemplary embodiment according to the present invention the step of deleting in detail.Fig. 2 illustrates the definite incorrect frontier point of another exemplary embodiment according to the present invention and the process flow diagram of the step of deleting.
In step 201, frontier point to the seismic trace obtaining in step 101, judge the seismic trace before the frontier point of each seismic trace and this seismic trace frontier point time difference absolute value and with this seismic trace after the absolute value of difference of time of frontier point of seismic trace at least one whether to be greater than Preset Time poor, if it is poor that at least one in two absolute values is greater than Preset Time, the frontier point of this seismic trace is confirmed as incorrect frontier point deleted.
Suppose, T (n) represents the time (that is, the time of the first arrival take-off of primary wave) of the frontier point of n seismic trace.Judgement | T (n)-T (n-1) | and | T (n+1)-T (n) | at least one whether be greater than the poor T of Preset Time
preset.If | T (n)-T (n-1) | and | T (n+1)-T (n) | at least one be greater than T
preset, this seismic trace, the frontier point of n seismic trace is incorrect frontier point and deletes.
When the frontier point of n seismic trace is confirmed as incorrect frontier point, delete the frontier point of n seismic trace.Therefore,, when whether the frontier point of determining (n+1) seismic trace is correct frontier point, the frontier point of (n+1) seismic trace was compared with the time of the frontier point of (n-1) seismic trace and (n+2) seismic trace.
Whether the frontier point that judges successively as stated above all seismic traces is incorrect frontier point and deletes incorrect frontier point.
In step 202, the frontier point of remaining seismic trace to passing through step 201, calculate the mean value of absolute value of mistiming of the frontier point of all adjacent seismic traces, then judge the seismic trace before the frontier point of each seismic trace and this seismic trace frontier point time difference absolute value and with this seismic trace after the absolute value of difference of time of frontier point of seismic trace at least one whether be greater than 4 times of described mean value, if be greater than 4 times, the frontier point of this seismic trace is confirmed as incorrect frontier point and deletes.
Suppose, through after step 201, the seismic trace quantity corresponding with remaining frontier point is N, and respectively be the 1st seismic trace, the 2nd seismic trace, the 3rd seismic trace ..., N-1 seismic trace and N seismic trace.
Calculate the mean value of absolute value of mistiming of the frontier point of adjacent seismic trace, that is,
Wherein, T (i) represents the time of the frontier point of i seismic trace.
Judgement | T (n)-T (n-1) | and | T (n+1)-T (n) | at least one whether be greater than 4T
ave.
If | T (n)-T (n-1) | and | T (n+1)-T (n) | at least one be greater than 4T
ave, the frontier point of this seismic trace, that is, the frontier point of n seismic trace is incorrect frontier point and deletes.Whether the frontier point of all seismic traces of judgement residue is as stated above incorrect frontier point and deletes incorrect frontier point.
In step 203, judgement is through step 202 and in remaining seismic trace, the absolute value of the time value difference of the frontier point of the seismic trace before the frontier point of each seismic trace and this seismic trace and with the absolute value of the difference of time of the frontier point of this seismic trace seismic trace afterwards at least one whether be greater than 1.5T
avedoubly, that is, and judgement | T (n)-T (n-1) | and | T (n+1)-T (n) | at least one whether be greater than 1.5T
ave.If at least one in remaining seismic trace meets above-mentioned condition, remaining seismic trace is re-executed to the method as described in step 202, until remaining frontier point meets above-mentioned condition.
In step 103, by the frontier point of remaining seismic trace is carried out to the frontier point that interpolation obtains all seismic traces;
In step 104, the frontier point of smoothing processing is mapped to the first arrival take-off that seismic wave figure picks up each seismic trace.
In step 105, according to the first arrival take-off of seismic trace, whether be the polarity that first arrival crest judges seismic trace afterwards.For each seismic trace, from frontier point, start scanning, when first scanning first arrival crest after frontier point, judge that this seismic trace is reversed polarity road, while first scanning first arrival trough on the contrary after frontier point, judge that this seismic trace is polarity road.
Alternatively, can be the polarity that first arrival trough or first arrival crest are determined this seismic trace with preset distance place, first arrival jump space by judgement.For example, if play preset distance place, jump space with first arrival, be crest (that is, first arrival crest), this seismic trace is reversed polarity road; If playing preset distance place, jump space with first arrival is trough (that is, first arrival trough), this seismic trace is confirmed as polarity road.
Fig. 3 illustrates to use according to of the present invention and automatically identifies based on primary wave the reversed polarity road that the method automatic decision in reversed polarity road goes out.In Fig. 3, X-axis represents seismic trace, and Y-axis represents the time.Square frame in Fig. 3 partly represents reversed polarity road.
According to the present invention, utilize Image Edge-Detection to carry out automatic decision to the first arrival of single big gun data, on the basis of automatic decision, recycling polarity amplitude characteristic is judged the polarity of seismic trace, thereby has improved the efficiency of the polarity of judgement seismic trace.
Although represented and described some embodiments of the present invention, it should be appreciated by those skilled in the art that in the situation that do not depart from principle of the present invention and the spirit that is limited its scope by claim and equivalent thereof, can modify to these embodiment.
Claims (3)
1. the method based on primary wave identification reversed polarity road, is characterized in that, comprises step:
The seismic wave figure of a plurality of seismic traces is converted into gray-scale map and gray-scale map is carried out to the frontier point that rim detection obtains each seismic trace;
Determine the incorrect frontier point in the frontier point obtaining and delete;
Frontier point by remaining seismic trace carries out the frontier point that interpolation obtains all seismic traces;
The frontier point of all seismic traces that obtain by interpolation is mapped to the first arrival take-off that seismic wave figure picks up each seismic trace; And
According to the first arrival take-off of seismic trace, whether be the polarity that first arrival crest judges seismic trace afterwards.
2. method according to claim 1, is characterized in that, the frontier point obtaining is carried out to conic fitting, and the frontier point that the distance that departs from fit line is surpassed to predetermined value is defined as incorrect frontier point and deletes.
3. method according to claim 1, is characterized in that, by judgement, is the polarity that first arrival trough or first arrival crest are determined this seismic trace with the preset distance place, first arrival jump space of each seismic trace.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110130303.1A CN102254321B (en) | 2011-05-19 | 2011-05-19 | Method for automatically identifying polarity-reversed channel based on first arrival wave |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110130303.1A CN102254321B (en) | 2011-05-19 | 2011-05-19 | Method for automatically identifying polarity-reversed channel based on first arrival wave |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102254321A CN102254321A (en) | 2011-11-23 |
| CN102254321B true CN102254321B (en) | 2014-03-26 |
Family
ID=44981562
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201110130303.1A Expired - Fee Related CN102254321B (en) | 2011-05-19 | 2011-05-19 | Method for automatically identifying polarity-reversed channel based on first arrival wave |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102254321B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106168679B (en) * | 2015-05-18 | 2018-03-09 | 中国石油化工股份有限公司 | Seismic acquisition records the processing method of polarity |
| CN108445533B (en) * | 2018-03-08 | 2019-06-07 | 青岛海洋地质研究所 | The method of long offset building phytal zone rate pattern based on OBS |
| CN110927782A (en) * | 2018-09-20 | 2020-03-27 | 中国石油化工股份有限公司 | Effective signal determination method and system for UNITE seismic acquisition instrument data |
| CN113589377B (en) * | 2020-04-30 | 2023-04-28 | 中国石油化工股份有限公司 | Seismic data polarity quantitative discrimination method based on seismic amplitude |
| CN112162321B (en) * | 2020-09-17 | 2023-03-14 | 中海油田服务股份有限公司 | Detection method and device for reverse-polarity seismic channel |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1395114A (en) * | 2001-07-05 | 2003-02-05 | 中国石油天然气集团公司 | Three-D earthquake prediction and observation system restoring technology |
| CN101663596A (en) * | 2006-11-03 | 2010-03-03 | 帕拉迪姆科学有限公司 | System and method for full azimuth angle domain imaging in reduced dimensional coordinate systems |
| CN101688926A (en) * | 2007-06-26 | 2010-03-31 | 申昌秀 | Method for velocity analysis using waveform inversion in laplace domain for geophysical imaging |
| CN102004264A (en) * | 2010-10-18 | 2011-04-06 | 中国石油化工股份有限公司 | Quantitative analysis and evaluation method for quality of acquired seismic data |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2549213C2 (en) * | 2007-11-14 | 2015-04-20 | ТЕРРАСПАРК ДЖИОСАЙЕНСИЗ, ЭлЭлСи | Processing seismic data |
| US8095345B2 (en) * | 2009-01-20 | 2012-01-10 | Chevron U.S.A. Inc | Stochastic inversion of geophysical data for estimating earth model parameters |
-
2011
- 2011-05-19 CN CN201110130303.1A patent/CN102254321B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1395114A (en) * | 2001-07-05 | 2003-02-05 | 中国石油天然气集团公司 | Three-D earthquake prediction and observation system restoring technology |
| CN101663596A (en) * | 2006-11-03 | 2010-03-03 | 帕拉迪姆科学有限公司 | System and method for full azimuth angle domain imaging in reduced dimensional coordinate systems |
| CN101688926A (en) * | 2007-06-26 | 2010-03-31 | 申昌秀 | Method for velocity analysis using waveform inversion in laplace domain for geophysical imaging |
| CN102004264A (en) * | 2010-10-18 | 2011-04-06 | 中国石油化工股份有限公司 | Quantitative analysis and evaluation method for quality of acquired seismic data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102254321A (en) | 2011-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102254321B (en) | Method for automatically identifying polarity-reversed channel based on first arrival wave | |
| KR101951034B1 (en) | Method for improving processing speed of OS-CFAR detection | |
| NO342539B1 (en) | Consistent fall angle estimation for seismic imaging | |
| CN102279973B (en) | Sea-sky-line detection method based on high gradient key points | |
| EP3059615B1 (en) | Amplitude-versus-angle analysis for quantative interpretation | |
| Forte et al. | Automated phase attribute-based picking applied to reflection seismics | |
| MX2010007253A (en) | Method for dynamic aperture determination for three-dimensional surface-related multiple elimination. | |
| CN103487788A (en) | Method for rapidly and automatically extracting train pulse signals | |
| CN104656139A (en) | Static correction method and equipment for converted wave | |
| US9523782B2 (en) | System and method for detection and classification of seismic terminations | |
| CN104463057B (en) | A kind of magnetic card coding/decoding method and device | |
| CN106353805A (en) | Intuitive microearthquake event classification method | |
| CN114841927A (en) | Shale reservoir fracture identification method, device, equipment and storage medium | |
| CN101923651A (en) | Calibration method of automatic recognition accuracy of corrupted data based on pavement images | |
| CN114236612A (en) | Signal processing method for earthquake background noise of water area in land | |
| CN103487804A (en) | Semi-automatic antarctic subglacial lake recognition method based on ice radar technology | |
| CN102338885A (en) | First-arrival-time automatic picking method of tricomponent VSP (Vertical Seismic Profile) data | |
| CN101813601B (en) | Sediment settling velocity measuring device with automatic measuring function and method thereof | |
| CN104570113A (en) | Method for self-adaptively removing strong reflection of earthquake | |
| US20120059592A1 (en) | Tracing seismic sections to convert to digital format | |
| CN114779211A (en) | Laser pulse radar equipment and point cloud density improving method and equipment | |
| CN103837890B (en) | Obtain the method and apparatus of geological data | |
| CN109655885A (en) | A kind of automatic first break pickup method and system | |
| CN103258196B (en) | In complex environment, be close to the character separation method of character string | |
| US10534102B2 (en) | Deblended and deghosted seismic data |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140326 Termination date: 20180519 |