CN105807319B - A kind of near-field wavelet bubble period pick-up method and system - Google Patents
A kind of near-field wavelet bubble period pick-up method and system Download PDFInfo
- Publication number
- CN105807319B CN105807319B CN201610141270.3A CN201610141270A CN105807319B CN 105807319 B CN105807319 B CN 105807319B CN 201610141270 A CN201610141270 A CN 201610141270A CN 105807319 B CN105807319 B CN 105807319B
- Authority
- CN
- China
- Prior art keywords
- seismic wave
- main peak
- bubble
- wave intensity
- maximum
- 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
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000005070 sampling Methods 0.000 claims abstract description 58
- 238000001514 detection method Methods 0.000 claims description 7
- 241000239290 Araneae Species 0.000 claims description 5
- 210000001367 artery Anatomy 0.000 claims 1
- 210000003462 vein Anatomy 0.000 claims 1
- 238000003908 quality control method Methods 0.000 description 6
- 238000010276 construction Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 229910002056 binary alloy Inorganic materials 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (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
A kind of near-field wavelet bubble period pick-up method of present invention proposition and system, including：Obtain the collected multiple seismic wave sampling point information of wave detector in preset time period.Traverse the two-dimensional array according to multiple seismic wave sampling point information creatings, obtain the maximum seismic wave intensity recorded in two-dimensional array, and the maximum seismic wave intensity is recorded as the first maximum seismic wave intensity, the wave crest where the first maximum seismic wave intensity is determined as wavelet pulse main peak.Obtain the rollback point of wavelet pulse main peak.By the rollback point, bubble main peak is searched in two-dimensional array.Bubble period is determined according to wavelet pulse main peak and bubble pulse main peak.Scheme through the invention can greatly improve near-field wavelet bubble period pickup velocity and accuracy rate.
Description
Technical field
The present invention relates to marine seismic prospectiong field more particularly to a kind of near-field wavelet bubble period pick-up method and it is
System.
Background technology
Marine seismic prospectiong focus is made of several high pressure airguns, and during oil-gas seismic exploration, focus is to data
Quality has decisive influence with construction efficiency.Before project kickoff, needs to carry out single rifle Quality Control to all air cannons, be examined by near field
The practical wavelet data of the SEGY formats of wave device record are compared with theoretical bubble value, and whether verification device parameter meets construction
Design.Traditional quality detecting method relies on business software, and human-computer interaction pickup is carried out under graphic user interface, and Quality Control timeliness is low, and
And accuracy rate is low.
Invention content
To solve the above-mentioned problems, the present invention proposes a kind of near-field wavelet bubble period pick-up method and system, can
Greatly improve near-field wavelet bubble period pickup velocity and accuracy rate.
In order to achieve the above object, the present invention proposes a kind of near-field wavelet bubble period pick-up method, and this method includes：
Obtain the collected multiple seismic wave sampling point information of wave detector in preset time period.
The two-dimensional array according to multiple seismic wave sampling point information creatings is traversed, the maximum earthquake recorded in two-dimensional array is obtained
Intensity of wave, and the maximum seismic wave intensity is recorded as the first maximum seismic wave intensity, it will be where the first maximum seismic wave intensity
Wave crest be determined as wavelet pulse main peak.
Obtain the rollback point of wavelet pulse main peak.
By the rollback point, bubble main peak is searched in two-dimensional array.
Bubble period is determined according to wavelet pulse main peak and bubble pulse main peak.
Preferably,
Obtaining the collected multiple seismic wave sampling point information of wave detector in preset time period includes：
Within the preset period, seismic wave is acquired according to the preset detection period by each channel of wave detector respectively
Information；Wherein, each channel is respectively provided with preset channel number.
Seismic wave intensity and the corresponding generation time point of seismic wave intensity are obtained from the seismic wave information of acquisition.
By each seismic wave intensity of acquisition generation time o'clock corresponding with seismic wave intensity as a seismic wave sample
Point information；Wherein, each seismic wave sampling point information is corresponding with the acquisition channel number of sampling point information.
Preferably, this method further includes：
The channel number in different channels on wave detector and multiple seismic wave sampling point information corresponding with each channel number are made
For dimension, two-dimensional array is created.
Wherein, different channel numbers is the first row of two-dimensional array；Multiple seismic wave sampling points corresponding with each channel number
Information is the row of two-dimensional array.
Seismic wave sampling point information in each row of two-dimensional array according to each seismic wave intensity generation time point in the past
Sequence backward arranges downwards successively.
Preferably, the rollback point of wavelet pulse main peak refers to：On wavelet pulse main peak, the first maximum seismic wave is appeared in
The point that first later seismic wave intensity of intensity is zero.
Obtain wavelet pulse main peak rollback point include：
The row where maximum seismic wave intensity are determined according to two-dimensional array.
Since the first maximum seismic wave intensity, along the Leie time search downwards where the first maximum seismic wave intensity, obtain
Take the point that first seismic wave intensity on the wavelet pulse main peak occurred after the first maximum seismic wave intensity is zero.
By rollback point, search bubble main peak includes in two-dimensional array：
By rollback point, along the Leie time search downwards where the first maximum seismic wave intensity until searching second
Wave crest, using second wave crest as bubble pulse main peak.
Preferably, determine that bubble period includes according to wavelet pulse main peak and bubble pulse main peak：
The corresponding generation time point of first maximum seismic wave intensity is recorded as the first generation time point.
It determines the maximum seismic wave intensity on bubble pulse main peak and is recorded as the second maximum seismic wave intensity；Most by second
The corresponding generation time point of violent earthquake intensity of wave is recorded as the second generation time point.
The time difference for calculating the first generation time point and the second generation time point, which is determined as bubble period.
In order to achieve the above object, the invention also provides a kind of near-field wavelet bubble period picking up system, the system packets
It includes：First acquisition module, spider module, the second acquisition module, search module and determining module.
First acquisition module, for obtaining the collected multiple seismic wave sampling point information of wave detector in preset time period.
Spider module obtains the two-dimensional array for traversing the two-dimensional array according to multiple seismic wave sampling point information creatings
The maximum seismic wave intensity of middle record, and the maximum seismic wave intensity is recorded as the first maximum seismic wave intensity, most by first
Wave crest where violent earthquake intensity of wave is determined as wavelet pulse main peak.
Second acquisition module, the rollback point for obtaining wavelet pulse main peak.
Search module searches for bubble main peak for by the rollback point in two-dimensional array.
Determining module, for determining bubble period according to wavelet pulse main peak and bubble pulse main peak.
Preferably, the first acquisition module obtains the collected multiple seismic wave sampling point packets of wave detector in preset time period
It includes：
Within the preset period, seismic wave is acquired according to the preset detection period by each channel of wave detector respectively
Information；Wherein, each channel is respectively provided with preset channel number.
Seismic wave intensity and the corresponding generation time point of the seismic wave intensity are obtained from the seismic wave information of acquisition.
By each seismic wave intensity of acquisition generation time o'clock corresponding with seismic wave intensity as a seismic wave sample
Point information；Wherein, each seismic wave sampling point information is corresponding with the acquisition channel number of the sampling point information.
Preferably, which further includes creation module.
Creation module is used for：By the channel number in different channels and multiple earthquakes corresponding with each channel number on wave detector
Wave sampling point information creates two-dimensional array as dimension.
Wherein, different channel numbers is the first row of two-dimensional array；Multiple seismic wave sampling points corresponding with each channel number
Information is the row of two-dimensional array.
Seismic wave sampling point information in each row of two-dimensional array according to each seismic wave intensity generation time point in the past
Sequence backward arranges downwards successively.
Preferably,
The rollback point of wavelet pulse main peak refers to：On wavelet pulse main peak, appear in the first maximum seismic wave intensity with
The point that first seismic wave intensity afterwards is zero.
Second acquisition module obtain wavelet pulse main peak rollback point include：
The row where maximum seismic wave intensity are determined according to two-dimensional array.
Since the first maximum seismic wave intensity, along the Leie time search downwards where the first maximum seismic wave intensity, obtain
Take the point that first seismic wave intensity on the wavelet pulse main peak occurred after the first maximum seismic wave intensity is zero.
Search module is by rollback point, search bubble main peak includes in two-dimensional array：
By rollback point, along the Leie time search downwards where the first maximum seismic wave intensity until searching second
Wave crest, using second wave crest as bubble pulse main peak.
Preferably, determining module determines that bubble period includes according to wavelet pulse main peak and bubble pulse main peak：
The corresponding generation time point of first maximum seismic wave intensity is recorded as the first generation time point.
It determines the maximum seismic wave intensity on bubble pulse main peak and is recorded as the second maximum seismic wave intensity；Most by second
The corresponding generation time point of violent earthquake intensity of wave is recorded as the second generation time point.
The time difference for calculating the first generation time point and the second generation time point, which is determined as bubble period.
Compared with prior art, the present invention includes：Obtain the collected multiple seismic wave samples of wave detector in preset time period
Point information.The two-dimensional array according to multiple seismic wave sampling point information creatings is traversed, the maximum earthquake recorded in two-dimensional array is obtained
Intensity of wave, and the maximum seismic wave intensity is recorded as the first maximum seismic wave intensity, it will be where the first maximum seismic wave intensity
Wave crest be determined as wavelet pulse main peak.Obtain the rollback point of wavelet pulse main peak.By the rollback point, in two-dimensional array
Search for bubble main peak.Bubble period is determined according to wavelet pulse main peak and bubble pulse main peak.Scheme through the invention, can
Greatly improve near-field wavelet bubble period pickup velocity and accuracy rate.
Description of the drawings
The attached drawing in the embodiment of the present invention is illustrated below, the attached drawing in embodiment be for the present invention into one
Step understands, for explaining the present invention together with specification, does not constitute limiting the scope of the invention.
Fig. 1 is the near-field wavelet bubble period pick-up method flow chart of the present invention；
Fig. 2 is the near-field wavelet bubble period picking up system composition frame chart of the present invention.
Specific implementation mode
For the ease of the understanding of those skilled in the art, the invention will be further described below in conjunction with the accompanying drawings, not
It can be used for limiting the scope of the invention.
In order to achieve the above object, the present invention proposes a kind of near-field wavelet bubble period pick-up method, as shown in Figure 1,
This method includes：
S101, the collected multiple seismic wave sampling point information of wave detector in preset time period are obtained.
Preferably,
Obtaining the collected multiple seismic wave sampling point information of wave detector in preset time period includes：
S1011, within the preset period, acquired respectively according to the preset detection period by each channel of wave detector
Seismic wave information；Wherein, each channel is respectively provided with preset channel number.Here the preset period can be according to user's need
It carries out self-defined.
S1012, seismic wave intensity and the corresponding generation time of seismic wave intensity are obtained from the seismic wave information of acquisition
Point.
S1013, by each seismic wave intensity of acquisition generation time o'clock corresponding with seismic wave intensity as a ground
Seismic wave sampling point information；Wherein, each seismic wave sampling point information is corresponding with the acquisition channel number of sampling point information.That is, wave detector
The seismic wave information of each channel acquisition can be recorded under corresponding channel number, and each channel number can be corresponded to according to time sequencing
Collected multiple seismic wave information are recorded, equally, each seismic wave sampling point information can be also recorded under corresponding channel number.
The two-dimensional array of S102, traversal according to multiple seismic wave sampling point information creatings, obtains and is recorded most in two-dimensional array
Violent earthquake intensity of wave, and the maximum seismic wave intensity is recorded as the first maximum seismic wave intensity, the first maximum seismic wave is strong
Wave crest where degree is determined as wavelet pulse main peak.
In embodiments of the present invention, before step of the present invention starts, the pre- multiple seismic wave samples for first passing through acquisition are needed
Point one two-dimensional array of information creating.
Specifically creation method includes：By the channel number in different channel on wave detector and corresponding with each channel number multiple
Seismic wave sampling point information creates two-dimensional array as dimension.
Wherein, different channel numbers is the first row of two-dimensional array；Multiple seismic wave sampling points corresponding with each channel number
Information is the row of two-dimensional array.
Seismic wave sampling point information in each row of two-dimensional array according to each seismic wave intensity generation time point in the past
Sequence backward arranges downwards successively.
On the basis of creating the two-dimensional array, according to energy spherical diffusion principle, and the inspection nearest apart from single rifle
Wave device receives the strongest principle of energy and will be looked for as long as we find the maximum seismic wave intensity recorded in two-dimensional array
To the channel number for acquiring the maximum seismic wave intensity, specifically, it is determined that the row where the maximum seismic wave intensity, according to its place
Row find the corresponding channel number of the row.
S103, the rollback point for obtaining wavelet pulse main peak.
Preferably, the rollback point of wavelet pulse main peak refers to：On wavelet pulse main peak, the first maximum seismic wave is appeared in
The point that first later seismic wave intensity of intensity is zero.
Polarity rollback point represents the end of the channel reception wavelet main pulse, is the important of wavelet main pulse and bubble pulse
Separation.
Obtain wavelet pulse main peak rollback point include：
S1031, row where maximum seismic wave intensity are determined according to two-dimensional array.
S1032, since the first maximum seismic wave intensity, searched downwards along the Leie where the first maximum seismic wave intensity
Rope obtains the point that first seismic wave intensity on the wavelet pulse main peak occurred after the first maximum seismic wave intensity is zero.
S104, by the rollback point, in two-dimensional array search for bubble main peak.
Preferably, by rollback point, search bubble main peak includes in two-dimensional array：
By rollback point, along the Leie time search downwards where the first maximum seismic wave intensity until searching second
Wave crest, using second wave crest as bubble pulse main peak.
In embodiments of the present invention, since bubble pulse main peak is the next wave crest for following wavelet pulse main peak closely,
As long as having searched temporally by sequential search backward first along bubble pulse main peak and having followed the next of wavelet pulse main peak closely
Wave crest has just searched bubble pulse main peak.
S105, bubble period is determined according to wavelet pulse main peak and bubble pulse main peak.
Preferably, determine that bubble period includes according to wavelet pulse main peak and bubble pulse main peak：
S1051, the corresponding generation time point of the first maximum seismic wave intensity is recorded as the first generation time point.
S1052, it determines the maximum seismic wave intensity on bubble pulse main peak and is recorded as the second maximum seismic wave intensity；It will
The corresponding generation time point of second maximum seismic wave intensity is recorded as the second generation time point.
S1053, the time difference for calculating the first generation time point and the second generation time point, which is determined as bubble
Period.
So far, we can obtain a bubble period, and cycle executes above-mentioned steps S101 to S105, until all gas
Dwell cycle pickup is completed, and performance data is obtained.The performance data is preserved and is exported in the form of text.Obtaining the achievement number
After according to this, performance data can be compared with preset Ct value, determines whether bubble rifle is qualified according to comparison result.
Here, which is the manufacture parameter value robbed by inspection bubble that related technical personnel pre-save.
In order to achieve the above object, the invention also provides a kind of near-field wavelet bubble period picking up system 01, such as Fig. 2 institutes
Show, which includes：First acquisition module 02, spider module 03, the second acquisition module 04, search module 05 and determining module
06。
First acquisition module 02, for obtaining the collected multiple seismic wave sampling point information of wave detector in preset time period.
Spider module 03 obtains the two-dimemsional number for traversing the two-dimensional array according to multiple seismic wave sampling point information creatings
The maximum seismic wave intensity recorded in group, and the maximum seismic wave intensity is recorded as the first maximum seismic wave intensity, by first
Wave crest where maximum seismic wave intensity is determined as wavelet pulse main peak.
Second acquisition module 04, the rollback point for obtaining wavelet pulse main peak.
Search module 05 searches for bubble main peak for by the rollback point in two-dimensional array.
Determining module 06, for determining bubble period according to wavelet pulse main peak and bubble pulse main peak.
Preferably, the first acquisition module 02 obtains the collected multiple seismic wave sampling point information of wave detector in preset time period
Including：
Within the preset period, seismic wave is acquired according to the preset detection period by each channel of wave detector respectively
Information；Wherein, each channel is respectively provided with preset channel number.
Seismic wave intensity and the corresponding generation time point of the seismic wave intensity are obtained from the seismic wave information of acquisition.
By each seismic wave intensity of acquisition generation time o'clock corresponding with seismic wave intensity as a seismic wave sample
Point information；Wherein, each seismic wave sampling point information is corresponding with the acquisition channel number of the sampling point information.
Preferably, which further includes creation module 07.
Creation module 07 is used for：By the channel number in different channel on wave detector and corresponding with each channel number multiple
Seismic wave sampling point information creates two-dimensional array as dimension.
Wherein, different channel numbers is the first row of two-dimensional array；Multiple seismic wave sampling points corresponding with each channel number
Information is the row of two-dimensional array.
Seismic wave sampling point information in each row of two-dimensional array according to each seismic wave intensity generation time point in the past
Sequence backward arranges downwards successively.
Preferably, the rollback point of wavelet pulse main peak refers to：On wavelet pulse main peak, the first maximum seismic wave is appeared in
The point that first later seismic wave intensity of intensity is zero.
Second acquisition module 04 obtain wavelet pulse main peak rollback point include：
The row where maximum seismic wave intensity are determined according to two-dimensional array.
Since the first maximum seismic wave intensity, along the Leie time search downwards where the first maximum seismic wave intensity, obtain
Take the point that first seismic wave intensity on the wavelet pulse main peak occurred after the first maximum seismic wave intensity is zero.
Search module 05 is by rollback point, search bubble main peak includes in two-dimensional array：
By rollback point, along the Leie time search downwards where the first maximum seismic wave intensity until searching second
Wave crest, using second wave crest as bubble pulse main peak.
Preferably, determining module 06 determines that bubble period includes according to wavelet pulse main peak and bubble pulse main peak：
The corresponding generation time point of first maximum seismic wave intensity is recorded as the first generation time point.
It determines the maximum seismic wave intensity on bubble pulse main peak and is recorded as the second maximum seismic wave intensity；Most by second
The corresponding generation time point of violent earthquake intensity of wave is recorded as the second generation time point.
The time difference for calculating the first generation time point and the second generation time point, which is determined as bubble period.
Compared with prior art, the present invention includes：Obtain the collected multiple seismic wave samples of wave detector in preset time period
Point information.The two-dimensional array according to multiple seismic wave sampling point information creatings is traversed, the maximum earthquake recorded in two-dimensional array is obtained
Intensity of wave, and the maximum seismic wave intensity is recorded as the first maximum seismic wave intensity, it will be where the first maximum seismic wave intensity
Wave crest be determined as wavelet pulse main peak.Obtain the rollback point of wavelet pulse main peak.By the rollback point, in two-dimensional array
Search for bubble main peak.Bubble period is determined according to wavelet pulse main peak and bubble pulse main peak.Scheme through the invention, can
Greatly improve near-field wavelet bubble period pickup velocity and accuracy rate.
In short, this method can high timeliness output Quality Control achievement：On the south for the exploratory area of Hypon island, conventional human's interactive mode
120 minutes pick-up times, an accuracy rate 94%；Quality control method using the present invention, pick-up time are only 0.05 minute, and
Accuracy rate 100%.
Wavelet bubble period automatic pick method based on matrix search relies on data-driven, directly believes from binary system wavelet
All bubble data are directly picked up in breath, are avoided conventional human from interacting the error of Quality Control, are also substantially had compressed single rifle Quality Control time.
It should be noted that embodiment described above be merely for convenience of it will be understood by those skilled in the art that, and
It is not used in and limits the scope of the invention, under the premise of not departing from the inventive concept of the present invention, those skilled in the art couple
Any obvious replacement and improvement that the present invention is made etc. are within protection scope of the present invention.
Claims (10)
1. a kind of near-field wavelet bubble period pick-up method, which is characterized in that the method includes：
The collected multiple seismic wave sampling point information of wave detector in preset time period are obtained, by each seismic wave intensity of acquisition
Generation time o'clock corresponding with the seismic wave intensity is as a seismic wave sampling point information；
The two-dimensional array according to the multiple seismic wave sampling point information creating is traversed, the maximum recorded in the two-dimensional array is obtained
Seismic wave intensity, and the maximum seismic wave intensity is recorded as the first maximum seismic wave intensity, by the described first maximum earthquake
Wave crest where intensity of wave is determined as wavelet pulse main peak, wherein in the two-dimensional array, different channel numbers is on wave detector
The first row of the two-dimensional array；The multiple seismic wave sampling point information corresponding with each channel number is the two-dimensional array
Row；
Obtain the rollback point of the wavelet pulse main peak；
By the rollback point, bubble main peak is searched in the two-dimensional array；
Bubble period is determined according to the wavelet pulse main peak and the bubble pulse main peak.
2. near-field wavelet bubble period pick-up method as described in claim 1, which is characterized in that
The collected multiple seismic wave sampling point information of wave detector include in the acquisition preset time period：
Within the preset period, seismic wave is acquired according to the preset detection period by each channel of wave detector respectively
Information；Wherein, each channel is respectively provided with preset channel number；
Seismic wave intensity and the corresponding generation time point of the seismic wave intensity are obtained from the seismic wave information of acquisition；
By each seismic wave intensity of acquisition generation time o'clock corresponding with the seismic wave intensity as a seismic wave sample
Point information；Wherein, each seismic wave sampling point information is corresponding with the channel number of sampling point information is acquired.
3. near-field wavelet bubble period pick-up method as claimed in claim 2, which is characterized in that the method further includes：
The channel number in different channels on the wave detector and multiple seismic wave sampling point information corresponding with each channel number are made
For dimension, two-dimensional array is created；
Wherein, the different channel number is the first row of the two-dimensional array；Multiple seismic waves corresponding with each channel number
Sampling point information is the row of the two-dimensional array；
Seismic wave sampling point information in each row of the two-dimensional array according to each seismic wave intensity generation time point in the past
Sequence backward arranges downwards successively.
4. near-field wavelet bubble period pick-up method as claimed in claim 3, which is characterized in that
The rollback point of the wavelet pulse main peak refers to：On the wavelet pulse main peak, the described first maximum earthquake is appeared in
The point that first later seismic wave intensity of intensity of wave is zero；
The rollback point for obtaining the wavelet pulse main peak includes：
The row where the maximum seismic wave intensity are determined according to the two-dimensional array；
Since the described first maximum seismic wave intensity, searched downwards along the Leie where the described first maximum seismic wave intensity
Rope, first seismic wave intensity obtained on the wavelet pulse main peak occurred after the described first maximum seismic wave intensity are
Zero point；
By the rollback point, search bubble main peak includes in the two-dimensional array：
By the rollback point, along the search downwards of the Leie time where the described first maximum seismic wave intensity until searching the
Two wave crests, using second wave crest as bubble pulse main peak.
5. near-field wavelet bubble period pick-up method as claimed in claim 2, which is characterized in that described according to the wavelet arteries and veins
It rushes main peak and the bubble pulse main peak determines that bubble period includes：
The corresponding generation time point of the first maximum seismic wave intensity is recorded as the first generation time point；
It determines the maximum seismic wave intensity on the bubble pulse main peak and is recorded as the second maximum seismic wave intensity；By described
The corresponding generation time point of two maximum seismic wave intensities is recorded as the second generation time point；
The time difference is determined as described by the time difference for calculating the first generation time point and the second generation time point
Bubble period.
6. a kind of near-field wavelet bubble period picking up system, which is characterized in that the system comprises：First acquisition module, traversal
Module, the second acquisition module, search module and determining module；
First acquisition module will for obtaining the collected multiple seismic wave sampling point information of wave detector in preset time period
Each seismic wave intensity generation time o'clock corresponding with the seismic wave intensity obtained is as a seismic wave sampling point information；
The spider module obtains described two for traversing the two-dimensional array according to the multiple seismic wave sampling point information creating
The maximum seismic wave intensity recorded in dimension group, and the maximum seismic wave intensity is recorded as the first maximum seismic wave intensity,
Wave crest where described first maximum seismic wave intensity is determined as wavelet pulse main peak, wherein in the two-dimensional array, detection
Different channel numbers is the first row of the two-dimensional array on device；The multiple seismic wave sampling point letter corresponding with each channel number
Breath is the row of the two-dimensional array；
Second acquisition module, the rollback point for obtaining the wavelet pulse main peak；
Described search module searches for bubble main peak for by the rollback point in the two-dimensional array；
The determining module, for determining bubble period according to the wavelet pulse main peak and the bubble pulse main peak.
7. near-field wavelet bubble period picking up system as claimed in claim 6, which is characterized in that first acquisition module obtains
The collected multiple seismic wave sampling point information of wave detector in preset time period are taken to include：
Within the preset period, seismic wave is acquired according to the preset detection period by each channel of wave detector respectively
Information；Wherein, each channel is respectively provided with preset channel number；
Seismic wave intensity and the corresponding generation time point of the seismic wave intensity are obtained from the seismic wave information of acquisition；
By each seismic wave intensity of acquisition generation time o'clock corresponding with the seismic wave intensity as a seismic wave sample
Point information；Wherein, each seismic wave sampling point information is corresponding with the channel number of sampling point information is acquired.
8. near-field wavelet bubble period picking up system as claimed in claim 7, which is characterized in that the system also includes establishments
Module；
Creation module is used for：By the channel number in different channels and multiple earthquakes corresponding with each channel number on the wave detector
Wave sampling point information creates two-dimensional array as dimension；
Wherein, the different channel number is the first row of the two-dimensional array；Multiple seismic waves corresponding with each channel number
Sampling point information is the row of the two-dimensional array；
Seismic wave sampling point information in each row of the two-dimensional array according to each seismic wave intensity generation time point in the past
Sequence backward arranges downwards successively.
9. near-field wavelet bubble period picking up system as claimed in claim 8, which is characterized in that
The rollback point of the wavelet pulse main peak refers to：On the wavelet pulse main peak, the described first maximum earthquake is appeared in
The point that first later seismic wave intensity of intensity of wave is zero；
The rollback point that second acquisition module obtains the wavelet pulse main peak includes：
The row where the maximum seismic wave intensity are determined according to the two-dimensional array；
Since the described first maximum seismic wave intensity, searched downwards along the Leie where the described first maximum seismic wave intensity
Rope, first seismic wave intensity obtained on the wavelet pulse main peak occurred after the described first maximum seismic wave intensity are
Zero point；
Described search module is by the rollback point, search bubble main peak includes in the two-dimensional array：
By the rollback point, along the search downwards of the Leie time where the described first maximum seismic wave intensity until searching the
Two wave crests, using second wave crest as bubble pulse main peak.
10. near-field wavelet bubble period picking up system as claimed in claim 7, which is characterized in that the determining module according to
The wavelet pulse main peak and the bubble pulse main peak determine that bubble period includes：
The corresponding generation time point of the first maximum seismic wave intensity is recorded as the first generation time point；
It determines the maximum seismic wave intensity on the bubble pulse main peak and is recorded as the second maximum seismic wave intensity；By described
The corresponding generation time point of two maximum seismic wave intensities is recorded as the second generation time point；
The time difference is determined as described by the time difference for calculating the first generation time point and the second generation time point
Bubble period.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610141270.3A CN105807319B (en) | 2016-03-11 | 2016-03-11 | A kind of near-field wavelet bubble period pick-up method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610141270.3A CN105807319B (en) | 2016-03-11 | 2016-03-11 | A kind of near-field wavelet bubble period pick-up method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105807319A CN105807319A (en) | 2016-07-27 |
| CN105807319B true CN105807319B (en) | 2018-11-09 |
Family
ID=56468246
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610141270.3A Active CN105807319B (en) | 2016-03-11 | 2016-03-11 | A kind of near-field wavelet bubble period pick-up method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105807319B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106443763B (en) * | 2016-08-31 | 2018-06-01 | 中国石油天然气集团公司 | A kind of the near-field signals acquisition methods and device of list rifle |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1794268A (en) * | 2005-12-21 | 2006-06-28 | 兰州大学 | Method of abstracting data from curve function recorded on cooraination peper |
| WO2010109280A2 (en) * | 2009-03-27 | 2010-09-30 | Geco Technology B.V. | Processing seismic data |
| CN103808492A (en) * | 2014-02-13 | 2014-05-21 | 上海温光自动化技术有限公司 | Curve display method of OTDR collected data |
-
2016
- 2016-03-11 CN CN201610141270.3A patent/CN105807319B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1794268A (en) * | 2005-12-21 | 2006-06-28 | 兰州大学 | Method of abstracting data from curve function recorded on cooraination peper |
| WO2010109280A2 (en) * | 2009-03-27 | 2010-09-30 | Geco Technology B.V. | Processing seismic data |
| CN103808492A (en) * | 2014-02-13 | 2014-05-21 | 上海温光自动化技术有限公司 | Curve display method of OTDR collected data |
Non-Patent Citations (3)
| Title |
|---|
| 气枪组合震源模拟分析及其应用;李长勤 等;《海洋石油》;20140930;第34卷(第3期);第101-105页 * |
| 气枪震源机理与子波处理;杨博;《中国优秀硕士学位论文全文数据库 基础科学辑》;20150615(第06期);第A011-108页 * |
| 蒙特卡洛法速度谱的自动拾取;张正和 等;《大庆石油地质与开发》;20020630;第21卷(第3期);第79-80页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105807319A (en) | 2016-07-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110058305A (en) | A kind of DAS seismic data noise-reduction method based on convolutional neural networks | |
| CN103758511B (en) | Method and device for identifying hidden reservoir through underground reverse time migration imaging | |
| CN103742131B (en) | With the time difference real time extracting method boring sound wave underground signal Acquire and process system | |
| CN102928873B (en) | Method for positioning ground micro-seismic based on four-dimensional energy focusing | |
| CN106154324A (en) | Down-hole micro-seismic event automatic identifying method based on multiple tracks scanning superposition | |
| CN111175815A (en) | Method and system for solving micro-seismic monitoring crack seismic source mechanism in oil reservoir transformation | |
| CN111538082B (en) | Automatic first arrival picking method for seismic wave time-frequency domain | |
| CN106291678A (en) | A kind of seismic data acquisition method and system | |
| CN103645248A (en) | High-temperature alloy grain size evaluation method based on ultrasonic phase velocity | |
| CN105807319B (en) | A kind of near-field wavelet bubble period pick-up method and system | |
| CN104182651B (en) | For the automatic quality control method in micro-seismic event azimuth that three-component geophone is received | |
| Li et al. | Automatic arrival-time picking of P-and S-waves of microseismic events based on object detection and CNN | |
| CN113534259B (en) | Real-time prestack time migration imaging method for efficient collection of controllable seismic source | |
| CN104570087A (en) | Method for performing seismic data first break and event pickup extraction through instantaneous attributes | |
| Ma et al. | Machine learning-assisted processing workflow for multi-fiber DAS microseismic data | |
| CN106593420A (en) | Multisource combined sound reflection logger and well side interface detection method | |
| CN114842280A (en) | Automatic micro-seismic signal identification algorithm based on convolutional neural network | |
| CN107229069B (en) | The method that a kind of pair of common receiver data trace gather carries out velocity analysis | |
| CN1465989A (en) | Multiwave multi-component earthquake data collecting system | |
| CN104865601B (en) | The method for determining the type of the earthquake record gathered in seismic prospecting | |
| CN113253335B (en) | Ground micro-seismic observation device and method | |
| CN111929721A (en) | Seismic information acquisition system and method | |
| CN111679321A (en) | Microseismic signal reconstruction method and system by using generalized gradient | |
| CN105911138A (en) | Impact imaging method and system based on average response energy | |
| CN112179566A (en) | Gun array seismic source air leakage early warning method and system for marine streamer acquisition |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing Applicant after: China Offshore Oil Group Co., Ltd. Applicant after: China Oilfield Services Limited Address before: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing Applicant before: China National Offshore Oil Corporation Applicant before: China Oilfield Services Limited |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |