CN105807319A - Near-field wavelet bubble period picking up method and system - Google Patents
Near-field wavelet bubble period picking up method and system Download PDFInfo
- Publication number
- CN105807319A CN105807319A CN201610141270.3A CN201610141270A CN105807319A CN 105807319 A CN105807319 A CN 105807319A CN 201610141270 A CN201610141270 A CN 201610141270A CN 105807319 A CN105807319 A CN 105807319A
- Authority
- CN
- China
- Prior art keywords
- seismic wave
- wave intensity
- main peak
- bubble
- wavelet
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000005070 sampling Methods 0.000 claims abstract description 62
- 241000239290 Araneae Species 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 6
- 238000003908 quality control method Methods 0.000 description 6
- 229910002056 binary alloy Inorganic materials 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 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
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
The invention provides a near-field wavelet bubble period picking up method and system. The method comprises the steps of: obtaining a plurality of pieces of seismic wave sampling point information collected by a detector in a preset time period; traversing a two-dimensional array created by the plurality of pieces of seismic wave sampling point information, obtaining a maximum seismic wave intensity recorded in the two-dimensional array, recording the maximum seismic wave intensity as a first maximum seismic wave intensity, and determining the wave peak where the first maximum seismic wave intensity is arranged to be a wavelet pulse main peak; obtaining a reversing point of the wavelet pulse main peak; starting from the reversing point, searching the two-dimensional array for a bubble main peak; and determining the bubble period according to the wavelet pulse main peak and the bubble pulse main peak. By adopting the scheme provided by the invention, the near-field wavelet bubble period picking up speed and accuracy are substantially improved.
Description
Technical field
The present invention relates to marine seismic prospectiong field, particularly relate to a kind of near-field wavelet bubble period pick-up method and system.
Background technology
Marine seismic prospectiong focus is made up of several high pressure airguns, and in oil-gas seismic exploration process, data quality and efficiency of construction are had decisive influence by focus.Before project kickoff, it is necessary to all air cannons carry out single rifle Quality Control, being contrasted with theoretical bubble value by the actual wavelet data of the SEGY form of near field cymoscope record, whether checking device parameter meets detail design.Tradition quality detecting method relies on business software, carries out man-machine interaction pickup under graphic user interface, as it is shown in figure 1, Quality Control timeliness is low, and accuracy rate is low.
Summary of the invention
In order to solve the problems referred to above, the present invention proposes a kind of near-field wavelet bubble period pick-up method and system, it is possible to be greatly improved 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 the method includes:
Obtain multiple seismic wave sampling point information that in preset time period, cymoscope collects.
Travel through the two-dimensional array according to multiple seismic wave sampling point information creatings, obtain the maximum seismic wave intensity of record in two-dimensional array, and this maximum seismic wave intensity is recorded as the first maximum seismic wave intensity, the crest at the first maximum seismic wave intensity place is defined as wavelet pulse main peak.
Obtain the reversal point of wavelet pulse main peak.
By this reversal point, two-dimensional array is searched for bubble main peak.
Bubble period is determined according to wavelet pulse main peak and bubble pulse main peak.
Preferably,
Obtain multiple seismic wave sampling point information that in preset time period, cymoscope collects to include:
Within the default time period, by each passage of cymoscope respectively according to default detection cycle locality seismic wave information;Wherein, each passage is respectively provided with default channel number.
Seismic wave intensity and generation time point corresponding to seismic wave intensity is obtained from the seismic wave information gathered.
Using each seismic wave intensity generation time point corresponding with seismic wave intensity of obtaining as a seismic wave sampling point information;Wherein, each seismic wave sampling point information is corresponding with the channel number gathering sampling point information.
Preferably, the method also includes:
Using the channel number of passages different on cymoscope and the multiple seismic wave sampling point information corresponding with each channel number as dimension, create two-dimensional array.
Wherein, different channel numbers is the first row of two-dimensional array;The multiple seismic wave sampling point information corresponding with each channel number are the row of two-dimensional array.
Seismic wave sampling point information in every string of two-dimensional array arranges downwards successively according to the time point order from front to back that produces of each seismic wave intensity.
Preferably, the reversal point of wavelet pulse main peak refers to: on wavelet pulse main peak, occurs in the point that first seismic wave intensity is zero that the first maximum seismic wave intensity is later.
The reversal point obtaining wavelet pulse main peak includes:
The row at maximum seismic wave intensity place are determined according to two-dimensional array.
From the first maximum seismic wave intensity, along the leu time search downwards at the first maximum seismic wave intensity place, the point that first seismic wave intensity is zero on wavelet pulse main peak occurred after obtaining the first maximum seismic wave intensity.
By reversal point, two-dimensional array is searched for bubble main peak and includes:
By reversal point, along the leu time search downwards at the first maximum seismic wave intensity place until searching second crest, using second crest as bubble pulse main peak.
Preferably, determine that bubble period includes according to wavelet pulse main peak and bubble pulse main peak:
Generation time point corresponding for first maximum seismic wave intensity is recorded as the first generation time point.
Determine the maximum seismic wave intensity on bubble pulse main peak and be recorded as the second maximum seismic wave intensity;Generation time point corresponding for second maximum seismic wave intensity is recorded as the second generation time point.
Calculate the first generation time point and second and produce the time difference of time point, this time difference is defined as bubble period.
In order to achieve the above object, the invention allows for a kind of near-field wavelet bubble period picking up system, this system includes: the first acquisition module, spider module, the second acquisition module, search module and determine module.
First acquisition module, for obtaining multiple seismic wave sampling point information that cymoscope in preset time period collects.
Spider module, for traveling through the two-dimensional array according to multiple seismic wave sampling point information creatings, obtain the maximum seismic wave intensity of record in this two-dimensional array, and this maximum seismic wave intensity is recorded as the first maximum seismic wave intensity, the crest at the first maximum seismic wave intensity place is defined as wavelet pulse main peak.
Second acquisition module, for obtaining the reversal point of wavelet pulse main peak.
Search module, for, by this reversal point, searching for bubble main peak in two-dimensional array.
Determine module, for determining bubble period according to wavelet pulse main peak and bubble pulse main peak.
Preferably, multiple seismic wave sampling point information that in the first acquisition module acquisition preset time period, cymoscope collects include:
Within the default time period, by each passage of cymoscope respectively according to default detection cycle locality seismic wave information;Wherein, each passage is respectively provided with default channel number.
Seismic wave intensity and generation time point corresponding to this seismic wave intensity is obtained from the seismic wave information gathered.
Using each seismic wave intensity generation time point corresponding with seismic wave intensity of obtaining as a seismic wave sampling point information;Wherein, each seismic wave sampling point information is corresponding with the channel number gathering this sampling point information.
Preferably, this system also includes creation module.
Creation module for: using the channel number of passages different on cymoscope and the multiple seismic wave sampling point information corresponding with each channel number as dimension, create two-dimensional array.
Wherein, different channel numbers is the first row of two-dimensional array;The multiple seismic wave sampling point information corresponding with each channel number are the row of two-dimensional array.
Seismic wave sampling point information in every string of two-dimensional array arranges downwards successively according to the time point order from front to back that produces of each seismic wave intensity.
Preferably,
The reversal point of wavelet pulse main peak refers to: on wavelet pulse main peak, occurs in the point that first seismic wave intensity is zero that the first maximum seismic wave intensity is later.
Second acquisition module obtains the reversal point of wavelet pulse main peak and includes:
The row at maximum seismic wave intensity place are determined according to two-dimensional array.
From the first maximum seismic wave intensity, along the leu time search downwards at the first maximum seismic wave intensity place, the point that first seismic wave intensity is zero on wavelet pulse main peak occurred after obtaining the first maximum seismic wave intensity.
Search module, by reversal point, is searched for bubble main peak in two-dimensional array and is included:
By reversal point, along the leu time search downwards at the first maximum seismic wave intensity place until searching second crest, using this second crest as bubble pulse main peak.
Preferably, it is determined that according to wavelet pulse main peak and bubble pulse main peak, module determines that bubble period includes:
Generation time point corresponding for first maximum seismic wave intensity is recorded as the first generation time point.
Determine the maximum seismic wave intensity on bubble pulse main peak and be recorded as the second maximum seismic wave intensity;Generation time point corresponding for second maximum seismic wave intensity is recorded as the second generation time point.
Calculate the first generation time point and second and produce the time difference of time point, this time difference is defined as bubble period.
Compared with prior art, the present invention includes: obtain multiple seismic wave sampling point information that in preset time period, cymoscope collects.Travel through the two-dimensional array according to multiple seismic wave sampling point information creatings, obtain the maximum seismic wave intensity of record in two-dimensional array, and this maximum seismic wave intensity is recorded as the first maximum seismic wave intensity, the crest at the first maximum seismic wave intensity place is defined as wavelet pulse main peak.Obtain the reversal point of wavelet pulse main peak.By this reversal point, two-dimensional array is searched for bubble main peak.Bubble period is determined according to wavelet pulse main peak and bubble pulse main peak.Pass through the solution of the present invention, it is possible to be greatly improved near-field wavelet bubble period pickup velocity and accuracy rate.
Accompanying drawing explanation
Below the accompanying drawing in the embodiment of the present invention being illustrated, the accompanying drawing in embodiment is for a further understanding of the present invention, is used for explaining the present invention, is not intended that limiting the scope of the invention together with description.
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.
Detailed description of the invention
For the ease of the understanding of those skilled in the art, below in conjunction with accompanying drawing, the invention will be further described, can not 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 it is shown in figure 1, the method includes:
Multiple seismic wave sampling point information that in S101, acquisition preset time period, cymoscope collects.
Preferably,
Obtain multiple seismic wave sampling point information that in preset time period, cymoscope collects to include:
S1011, within the default time period, by each passage of cymoscope respectively according to default detection cycle locality seismic wave information;Wherein, each passage is respectively provided with default channel number.Here this time period preset can need to carry out according to user self-defined.
S1012, from the seismic wave information gathered, obtain seismic wave intensity and generation time point corresponding to seismic wave intensity.
S1013, using each seismic wave intensity generation time point corresponding with seismic wave intensity of obtaining as a seismic wave sampling point information;Wherein, each seismic wave sampling point information is corresponding with the channel number gathering sampling point information.Namely, the seismic wave information of each channel acquisition of cymoscope can record under corresponding channel number, multiple seismic wave information that each channel number can collect according to time sequencing corresponding record, equally, each seismic wave sampling point information also can record under corresponding channel number.
S102, traversal are according to the two-dimensional array of multiple seismic wave sampling point information creatings, obtain the maximum seismic wave intensity of record in two-dimensional array, and this maximum seismic wave intensity is recorded as the first maximum seismic wave intensity, the crest at the first maximum seismic wave intensity place is defined as wavelet pulse main peak.
In embodiments of the present invention, before step of the present invention starts, it is necessary to beforehand through multiple one two-dimensional array of seismic wave sampling point information creating gathered.
Concrete creation method includes: using the channel number of passages different on cymoscope and the multiple seismic wave sampling point information corresponding with each channel number as dimension, create two-dimensional array.
Wherein, different channel numbers is the first row of two-dimensional array;The multiple seismic wave sampling point information corresponding with each channel number are the row of two-dimensional array.
Seismic wave sampling point information in every string of two-dimensional array arranges downwards successively according to the time point order from front to back that produces of each seismic wave intensity.
On the basis creating this two-dimensional array, according to energy spherical diffusion principle, and receive, apart from the cymoscope that single rifle is nearest, the principle that energy is the strongest, as long as we find the maximum seismic wave intensity recorded in two-dimensional array, the channel number gathering this maximum seismic wave intensity will be found, specifically, it is determined that the row at this maximum seismic wave intensity place, find, according to the row at its place, the channel number that these row are corresponding.
S103, obtain wavelet pulse main peak reversal point.
Preferably, the reversal point of wavelet pulse main peak refers to: on wavelet pulse main peak, occurs in the point that first seismic wave intensity is zero that the first maximum seismic wave intensity is later.
Polarity inversion point represents the end of this channel reception wavelet main pulse, is the important separation of wavelet main pulse and bubble pulse.
The reversal point obtaining wavelet pulse main peak includes:
S1031, determine the row at maximum seismic wave intensity place according to two-dimensional array.
S1032, from the first maximum seismic wave intensity, along the leu time search downwards at the first maximum seismic wave intensity place, the point that first seismic wave intensity is zero on wavelet pulse main peak occurred after obtaining the first maximum seismic wave intensity.
S104, by this reversal point, in two-dimensional array search for bubble main peak.
Preferably, by reversal point, two-dimensional array is searched for bubble main peak and includes:
By reversal point, along the leu time search downwards at the first maximum seismic wave intensity place until searching second crest, using this second crest as bubble pulse main peak.
In embodiments of the present invention, owing to bubble pulse main peak is the next crest following wavelet pulse main peak closely, therefore, as long as along bubble pulse main peak temporally by first sequential search backward, having searched the next crest following wavelet pulse main peak closely and just searched bubble pulse main peak.
S105, determine bubble period 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, generation time point corresponding for the first maximum seismic wave intensity is recorded as the first generation time point.
S1052, the maximum seismic wave intensity determined on bubble pulse main peak are also recorded as the second maximum seismic wave intensity;Generation time point corresponding for second maximum seismic wave intensity is recorded as the second generation time point.
S1053, calculating the first generation time point and second produce the time difference of time point, and this time difference is defined as bubble period.
So far, we just can obtain a bubble period, and circulation performs above-mentioned steps S101 to S105, until all bubble period have been picked up, obtains performance data.This performance data is preserved and exports in the form of text.After obtaining this performance data, just by performance data compared with the Ct value preset, can determine that whether bubble rifle is qualified according to comparative result.Here, this Ct value preset is the parameter value that dispatches from the factory robbed by inspection bubble that person skilled pre-saves.
In order to achieve the above object, the invention allows for a kind of near-field wavelet bubble period picking up system 01, as in figure 2 it is shown, this system includes: the first acquisition module 02, spider module the 03, second acquisition module 04, search module 05 and determine module 06.
First acquisition module 02, for obtaining multiple seismic wave sampling point information that cymoscope in preset time period collects.
Spider module 03, for traveling through the two-dimensional array according to multiple seismic wave sampling point information creatings, obtain the maximum seismic wave intensity of record in this two-dimensional array, and this maximum seismic wave intensity is recorded as the first maximum seismic wave intensity, the crest at the first maximum seismic wave intensity place is defined as wavelet pulse main peak.
Second acquisition module 04, for obtaining the reversal point of wavelet pulse main peak.
Search module 05, for, by this reversal point, searching for bubble main peak in two-dimensional array.
Determine module 06, for determining bubble period according to wavelet pulse main peak and bubble pulse main peak.
Preferably, the first acquisition module 02 obtains multiple seismic wave sampling point information that cymoscope in preset time period collects and includes:
Within the default time period, by each passage of cymoscope respectively according to default detection cycle locality seismic wave information;Wherein, each passage is respectively provided with default channel number.
Seismic wave intensity and generation time point corresponding to this seismic wave intensity is obtained from the seismic wave information gathered.
Using each seismic wave intensity generation time point corresponding with seismic wave intensity of obtaining as a seismic wave sampling point information;Wherein, each seismic wave sampling point information is corresponding with the channel number gathering this sampling point information.
Preferably, this system also includes creation module 07.
Creation module 07 for: using the channel number of passages different on cymoscope and the multiple seismic wave sampling point information corresponding with each channel number as dimension, create two-dimensional array.
Wherein, different channel numbers is the first row of two-dimensional array;The multiple seismic wave sampling point information corresponding with each channel number are the row of two-dimensional array.
Seismic wave sampling point information in every string of two-dimensional array arranges downwards successively according to the time point order from front to back that produces of each seismic wave intensity.
Preferably, the reversal point of wavelet pulse main peak refers to: on wavelet pulse main peak, occurs in the point that first seismic wave intensity is zero that the first maximum seismic wave intensity is later.
Second acquisition module 04 obtains the reversal point of wavelet pulse main peak and includes:
The row at maximum seismic wave intensity place are determined according to two-dimensional array.
From the first maximum seismic wave intensity, along the leu time search downwards at the first maximum seismic wave intensity place, the point that first seismic wave intensity is zero on wavelet pulse main peak occurred after obtaining the first maximum seismic wave intensity.
Search module 05, by reversal point, is searched for bubble main peak in two-dimensional array and is included:
By reversal point, along the leu time search downwards at the first maximum seismic wave intensity place until searching second crest, using this second crest as bubble pulse main peak.
Preferably, it is determined that according to wavelet pulse main peak and bubble pulse main peak, module 06 determines that bubble period includes:
Generation time point corresponding for first maximum seismic wave intensity is recorded as the first generation time point.
Determine the maximum seismic wave intensity on bubble pulse main peak and be recorded as the second maximum seismic wave intensity;Generation time point corresponding for second maximum seismic wave intensity is recorded as the second generation time point.
Calculate the first generation time point and second and produce the time difference of time point, this time difference is defined as bubble period.
Compared with prior art, the present invention includes: obtain multiple seismic wave sampling point information that in preset time period, cymoscope collects.Travel through the two-dimensional array according to multiple seismic wave sampling point information creatings, obtain the maximum seismic wave intensity of record in two-dimensional array, and this maximum seismic wave intensity is recorded as the first maximum seismic wave intensity, the crest at the first maximum seismic wave intensity place is defined as wavelet pulse main peak.Obtain the reversal point of wavelet pulse main peak.By this reversal point, two-dimensional array is searched for bubble main peak.Bubble period is determined according to wavelet pulse main peak and bubble pulse main peak.Pass through the solution of the present invention, it is possible to be greatly improved near-field wavelet bubble period pickup velocity and accuracy rate.
In a word, the method can high timeliness output Quality Control achievement: on the south exploratory area, Hypon island be example, 120 minutes conventional human's interactive mode pick-up times, an accuracy rate 94%;Adopting the quality control method of the present invention, the pick-up time is only 0.05 minute, and an accuracy rate 100%.
Wavelet bubble period automatic pick method based on matrix search relies on data-driven, directly directly picks up all bubble data from binary system wavelet information, it is to avoid the error of the mutual Quality Control of conventional human, also significantly have compressed single rifle Quality Control time.
It should be noted that; embodiment described above is for only for ease of those skilled in the art and understands; it is not limited to protection scope of the present invention; under the premise without departing from the inventive concept of the present invention, those skilled in the art to the made any apparent replacement and improvement etc. of the present invention all within protection scope of the present invention.
Claims (10)
1. a near-field wavelet bubble period pick-up method, it is characterised in that described method includes:
Obtain multiple seismic wave sampling point information that in preset time period, cymoscope collects;
Travel through the two-dimensional array according to the plurality of seismic wave sampling point information creating, obtain the maximum seismic wave intensity of record in described two-dimensional array, and described maximum seismic wave intensity is recorded as the first maximum seismic wave intensity, the crest at described first maximum seismic wave intensity place is defined as wavelet pulse main peak;
Obtain the reversal point of described wavelet pulse main peak;
By described reversal point, described two-dimensional array is searched for bubble main peak;
Bubble period is determined according to described wavelet pulse main peak and described bubble pulse main peak.
2. near-field wavelet bubble period pick-up method as claimed in claim 1, it is characterised in that
Multiple seismic wave sampling point information that in described acquisition preset time period, cymoscope collects include:
Within the described default time period, by each passage of cymoscope respectively according to default detection cycle locality seismic wave information;Wherein, described each passage is respectively provided with default channel number;
Seismic wave intensity and generation time point corresponding to described seismic wave intensity is obtained from the described seismic wave information gathered;
Using each seismic wave intensity generation time point corresponding with described seismic wave intensity of obtaining as a seismic wave sampling point information;Wherein, each seismic wave sampling point information is corresponding with the channel number gathering described sampling point information.
3. near-field wavelet bubble period pick-up method as claimed in claim 2, it is characterised in that described method also includes:
Using the channel number of passages different on described cymoscope and the multiple seismic wave sampling point information corresponding with each channel number as dimension, create two-dimensional array;
Wherein, described different channel number is the first row of described two-dimensional array;The row that multiple seismic wave sampling point information be described two-dimensional array corresponding with each channel number;
Seismic wave sampling point information in every string of described two-dimensional array arranges downwards successively according to the time point order from front to back that produces of each seismic wave intensity.
4. near-field wavelet bubble period pick-up method as claimed in claim 3, it is characterised in that
The reversal point of described wavelet pulse main peak refers to: on described wavelet pulse main peak, occurs in the point that first seismic wave intensity is zero that described first maximum seismic wave intensity is later;
The reversal point of the described wavelet pulse main peak of described acquisition includes:
The row at described maximum seismic wave intensity place are determined according to described two-dimensional array;
From described first maximum seismic wave intensity, along the leu time search downwards at described first maximum seismic wave intensity place, the point that first seismic wave intensity is zero on described wavelet pulse main peak occurred after obtaining described first maximum seismic wave intensity;
By described reversal point, described two-dimensional array is searched for bubble main peak and includes:
By described reversal point, along the leu time search downwards at described first maximum seismic wave intensity place until searching second crest, using described second crest as bubble pulse main peak.
5. near-field wavelet bubble period pick-up method as claimed in claim 2, it is characterised in that described determine that bubble period includes according to described wavelet pulse main peak and described bubble pulse main peak:
Generation time point corresponding for described first maximum seismic wave intensity is recorded as the first generation time point;
Determine the maximum seismic wave intensity on described bubble pulse main peak and be recorded as the second maximum seismic wave intensity;Generation time point corresponding for described second maximum seismic wave intensity is recorded as the second generation time point;
Calculate the described first time difference producing time point and described second generation time point, described time difference is defined as described bubble period.
6. a near-field wavelet bubble period picking up system, it is characterised in that described system includes: the first acquisition module, spider module, the second acquisition module, search module and determine module;
Described first acquisition module, for obtaining multiple seismic wave sampling point information that cymoscope in preset time period collects;
Described spider module, for traveling through the two-dimensional array according to the plurality of seismic wave sampling point information creating, obtain the maximum seismic wave intensity of record in described two-dimensional array, and described maximum seismic wave intensity is recorded as the first maximum seismic wave intensity, the crest at described first maximum seismic wave intensity place is defined as wavelet pulse main peak;
Described second acquisition module, for obtaining the reversal point of described wavelet pulse main peak;
Described search module, for, by described reversal point, searching for bubble main peak in described two-dimensional array;
Described determine module, for determining bubble period according to described wavelet pulse main peak and described bubble pulse main peak.
7. near-field wavelet bubble period picking up system as claimed in claim 6, it is characterised in that described first acquisition module obtains multiple seismic wave sampling point information that in preset time period, cymoscope collects and includes:
Within the described default time period, by each passage of cymoscope respectively according to default detection cycle locality seismic wave information;Wherein, described each passage is respectively provided with default channel number;
Seismic wave intensity and generation time point corresponding to described seismic wave intensity is obtained from the described seismic wave information gathered;
Using each seismic wave intensity generation time point corresponding with described seismic wave intensity of obtaining as a seismic wave sampling point information;Wherein, each seismic wave sampling point information is corresponding with the channel number gathering described sampling point information.
8. near-field wavelet bubble period pick-up method as claimed in claim 7, it is characterised in that described system also includes creation module;
Creation module for: using the channel number of passages different on described cymoscope and the multiple seismic wave sampling point information corresponding with each channel number as dimension, create two-dimensional array;
Wherein, described different channel number is the first row of described two-dimensional array;The row that multiple seismic wave sampling point information be described two-dimensional array corresponding with each channel number;
Seismic wave sampling point information in every string of described two-dimensional array arranges downwards successively according to the time point order from front to back that produces of each seismic wave intensity.
9. near-field wavelet bubble period picking up system as claimed in claim 8, it is characterised in that
The reversal point of described wavelet pulse main peak refers to: on described wavelet pulse main peak, occurs in the point that first seismic wave intensity is zero that described first maximum seismic wave intensity is later;
The reversal point that described second acquisition module obtains described wavelet pulse main peak includes:
The row at described maximum seismic wave intensity place are determined according to described two-dimensional array;
From described first maximum seismic wave intensity, along the leu time search downwards at described first maximum seismic wave intensity place, the point that first seismic wave intensity is zero on described wavelet pulse main peak occurred after obtaining described first maximum seismic wave intensity;
Described search module, by described reversal point, is searched for bubble main peak in described two-dimensional array and is included:
By described reversal point, along the leu time search downwards at described first maximum seismic wave intensity place until searching second crest, using described second crest as bubble pulse main peak.
10. near-field wavelet bubble period picking up system as claimed in claim 7, it is characterised in that described determine according to described wavelet pulse main peak and described bubble pulse main peak, module determines that bubble period includes:
Generation time point corresponding for described first maximum seismic wave intensity is recorded as the first generation time point;
Determine the maximum seismic wave intensity on described bubble pulse main peak and be recorded as the second maximum seismic wave intensity;Generation time point corresponding for described second maximum seismic wave intensity is recorded as the second generation time point;
Calculate the described first time difference producing time point and described second generation time point, described time difference is defined as described 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 true CN105807319A (en) | 2016-07-27 |
CN105807319B 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) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106443763A (en) * | 2016-08-31 | 2017-02-22 | 中国石油天然气集团公司 | Near-field signal obtaining method and device for single gun |
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 |
---|
张正和 等: "蒙特卡洛法速度谱的自动拾取", 《大庆石油地质与开发》 * |
李长勤 等: "气枪组合震源模拟分析及其应用", 《海洋石油》 * |
杨博: "气枪震源机理与子波处理", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106443763A (en) * | 2016-08-31 | 2017-02-22 | 中国石油天然气集团公司 | Near-field signal obtaining method and device for single gun |
CN106443763B (en) * | 2016-08-31 | 2018-06-01 | 中国石油天然气集团公司 | A kind of the near-field signals acquisition methods and device of list rifle |
Also Published As
Publication number | Publication date |
---|---|
CN105807319B (en) | 2018-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102928873B (en) | Method for positioning ground micro-seismic based on four-dimensional energy focusing | |
CN110058305A (en) | A kind of DAS seismic data noise-reduction method based on convolutional neural networks | |
CN102565857B (en) | Automatic remaining dynamic correction method | |
CN103399300A (en) | Wave packet superposition microseism ground location method | |
CN104570076A (en) | Automatic seismic wave first-arrival picking method based on dichotomy | |
CN105301641A (en) | Azimuthal anisotropy speed inversion method and apparatus | |
CN106777060A (en) | Node type seismic detector data file management system and method | |
CN106959466A (en) | Marine seismic acquisition system and method | |
CN103163567A (en) | Air gun wavelet bubble ratio analysis method and device based on energy | |
CN104182651B (en) | For the automatic quality control method in micro-seismic event azimuth that three-component geophone is received | |
CN104570087A (en) | Method for performing seismic data first break and event pickup extraction through instantaneous attributes | |
CN104133244B (en) | A kind of near surface static corrections processing method of focus shear wave exploration | |
CN105807319A (en) | Near-field wavelet bubble period picking up method and system | |
CN111308548A (en) | High-precision microseism data first arrival pickup device, system and method | |
CN108375789B (en) | Synchronous matching method for jointly acquiring seismic data | |
CN104199105A (en) | Pickup method and device of first arrivals in high-density three-dimensional exploration | |
CN113534259A (en) | Vibroseis efficient acquisition real-time prestack time migration imaging method | |
CN106837313B (en) | LWF memory-type acoustic logging slowness extracting method | |
CN105738949B (en) | A kind of nine bin uniformity method for parallel processing for time-lapse seismic | |
CN106646592B (en) | The real-time stacking method of single shot record | |
CN104280766B (en) | A kind of direct offset method of utilization local data lineups slope | |
CN104865601B (en) | The method for determining the type of the earthquake record gathered in seismic prospecting | |
CN105572733A (en) | Seismic velocity spectrum automatic pick-up method | |
CN110161560B (en) | Method and device for detecting microseismic event | |
CN110703319B (en) | Microseismic event detection method and system based on Khichin-Einstein theorem |
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 |