CN105158801B - The compression process method and apparatus of optical cable coupled noise - Google Patents
The compression process method and apparatus of optical cable coupled noise Download PDFInfo
- Publication number
- CN105158801B CN105158801B CN201510453845.0A CN201510453845A CN105158801B CN 105158801 B CN105158801 B CN 105158801B CN 201510453845 A CN201510453845 A CN 201510453845A CN 105158801 B CN105158801 B CN 105158801B
- Authority
- CN
- China
- Prior art keywords
- data
- noise
- curve
- optical cable
- cable coupled
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a kind of compression process method and apparatus of optical cable coupled noise, wherein this method includes：From the vertical seismic data gathered by distribution type fiber-optic, determine occur the interference range of optical cable coupled noise；It is determined that interference range in read noise periods number and corresponding time span, and calculate and obtain single noise periods, obtain nonmonotonic time curve；According to nonmonotonic time curve, amplitude attenuation curve is obtained；According to amplitude attenuation curve, calculating obtains average wavelet；Non-monotonic time curve and average wavelet are subjected to convolution operation, fitting noisy data is obtained；Fitting noisy data is subtracted by vertical seismic data, remaining data is obtained；Remaining data is filtered, the data after the compacting of optical cable coupled noise are obtained.The present invention solves the problems, such as in vertical seismic data to produce the compacting of interference noise because optic cable can not be close to the borehole wall, has reached the technique effect of the optical cable coupled noise in effective elimination vertical seismic data.
Description
Technical field
The present invention relates to oil exploration technology field, the compression process method and dress of more particularly to a kind of optical cable coupled noise
Put.
Background technology
Distribution type fiber-optic is progressively applied to multiple data acquisitions and transmission industry in recent years as a kind of emerging equipment,
Its operation principle is to be based on scattering effect, and union senses and is transmitted in one, and telemeasurement and monitoring can be achieved.For distribution
Formula optical fiber once determines the spatial parameter distribution map for being obtained with whole fiber distribution region, and can be with based on distribution type fiber-optic
Measure the information up to tens kilometers.Distribution type fiber-optic have sensitivity height, electromagnetism interference, good insulating, it is corrosion-resistant and
The advantages of being easy to networking and transmit over long distances, is especially suitable for applying the field being restricted in some traditional sensors.Because hanging down
Straight earthquake data acquisition usually requires going deep into underground number kilometer and temperature and pressure is higher, the borehole wall is in the environment of steel sleeve
Carry out, therefore, the advantage in this observed pattern with uniqueness of distribution type fiber-optic.
However, due to when current distribution type fiber-optic is gathered in underground, there is no backup equipment, therefore exist in record data
Optical cable coupled noise effects, and the noise is close with effective earthquake reflected wave apparent velocity, therefore, it is difficult to separate in record data
Optical cable coupled noise.
In view of the above-mentioned problems, not yet proposing effective solution at present.
The content of the invention
The embodiments of the invention provide a kind of compression process method of optical cable coupled noise, to reach that effectively compacting is distributed
The purpose of optical cable coupled noise in the geological data that optical fiber is measured, this method includes：
From the vertical seismic data gathered by distribution type fiber-optic, determine occur the interference range of optical cable coupled noise；
It is determined that interference range in read noise periods number and corresponding time span, and calculate and obtain single noise week
Phase；
By the interference range and single noise periods of determination, by the way of mathematical interpolation, away from song when obtaining nonmonotonic
Line；
According to the nonmonotonic time curve, amplitude attenuation curve is obtained；
The complete region of lineups waveform is selected from the interference range, according to the amplitude attenuation curve, calculating is obtained
Average wavelet；
The non-monotonic time curve and the average wavelet are subjected to convolution operation, fitting noisy data is obtained；
The fitting noisy data is subtracted by the vertical seismic data, remaining data is obtained；
The remaining data is filtered, the data after the compacting of optical cable coupled noise are obtained.
In one embodiment, according to the nonmonotonic time curve, amplitude attenuation curve is obtained, including：
According to the nonmonotonic time curve, remove in the seismic profile of the vertical seismic data with it is described non-monotonic
The corresponding sample value of time curve；
Function Fitting is carried out to removing the seismic profile after sample value, amplitude attenuation curve is obtained.
In one embodiment, single noise periods are calculated according to below equation：
Wherein, Δ T represents single noise periods, and N represents the noise periods number in the interference range, and T represents N number of noise week
Phase corresponding time span.
In one embodiment, the remaining data is filtered, including：
Bandpass filtering is carried out to the remaining data.
In one embodiment, the band logical used during bandpass filtering is carried out to the remaining data is scanned
Parameter is determined by the frequency range of the vertical seismic data.
The embodiment of the present invention additionally provides a kind of compression process device of optical cable coupled noise, to reach that effectively compacting is distributed
The purpose of optical cable coupled noise in the geological data that formula optical fiber is measured, the device includes：
Interference range determining module, for from the vertical seismic data gathered by distribution type fiber-optic, determining generation light
The interference range of cable coupled noise；
Single noise periods determining module, for it is determined that interference range in read noise periods number and the corresponding time it is long
Degree, and calculating obtains single noise periods；
Time curve curve determining module, for the interference range by determination and single noise periods, using mathematical interpolation
Mode, obtain nonmonotonic time curve；
Amplitude attenuation curve determining module, for according to the nonmonotonic time curve, obtaining amplitude attenuation curve；
Average wavelet computing module, for selecting the complete region of lineups waveform from the interference range, according to described
Amplitude attenuation curve, calculating obtains average wavelet；
Fitting noisy data acquisition module, for the non-monotonic time curve and the average wavelet to be carried out into convolution fortune
Calculate, obtain fitting noisy data；
Remaining data acquisition module, for subtracting the fitting noisy data by the vertical seismic data, is remained
Remainder evidence；
Filtration module, for being filtered to the remaining data, obtains the data after the compacting of optical cable coupled noise.
In one embodiment, amplitude attenuation curve determining module includes：
Removal unit, for according to the nonmonotonic time curve, removing the seismic profile of the vertical seismic data
In sample value corresponding with the non-monotonic time curve；
Amplitude attenuation curve determining unit, for carrying out Function Fitting to removing the seismic profile after sample value, is shaken
Width attenuation curve.
In one embodiment, the single noise periods determining module according to below equation specifically for calculating single
Noise periods：
Wherein, Δ T represents single noise periods, and N represents the noise periods number in the interference range, and T represents N number of noise week
Phase corresponding time span.
In one embodiment, the filtration module, specifically for by way of bandpass filtering to the remainder
According to being filtered.
In one embodiment, the filtration module is adopted during bandpass filtering is carried out to the remaining data
Band logical sweep parameter is determined by the frequency range of the vertical seismic data.
In embodiments of the present invention, the optical cable coupled noise in vertical seismic data is obtained by fitting, then by subtracting
The mode of method is suppressed noise, so as to solve in vertical seismic data to produce because optic cable can not be close to the borehole wall
The compacting problem of interference noise, has reached the technique effect of the optical cable coupled noise in effective elimination vertical seismic data.
Brief description of the drawings
Accompanying drawing described herein is used for providing a further understanding of the present invention, constitutes the part of the application, not
Constitute limitation of the invention.In the accompanying drawings：
Fig. 1 is the method flow diagram of the compression process method of optical cable coupled noise according to embodiments of the present invention；
Fig. 2 is the another method flow chart of the compression process method of optical cable coupled noise according to embodiments of the present invention；
Fig. 3 is that the vertical seismic data original record according to embodiments of the present invention disturbed by optical cable coupled noise is illustrated
Figure；
Fig. 4 is the time span T of N number of noise periods number according to embodiments of the present invention statistics schematic diagram；
Fig. 5 is that non-monotonic time curve according to embodiments of the present invention sets up schematic diagram；
Fig. 6 is optical cable coupled noise attenuation law statistics schematic diagram according to embodiments of the present invention；
Fig. 7 is that schematic diagram is chosen in the optical cable coupled noise interval of statistical wavelet according to embodiments of the present invention；
Fig. 8 is the average wavelet schematic diagram that superposition according to embodiments of the present invention is obtained；
Fig. 9 is the optical cable coupled noise schematic diagram obtained by fitting according to embodiments of the present invention；
Figure 10 is that initial data according to embodiments of the present invention subtracts the obtained remaining data of fitting optical cable coupled noise and shown
It is intended to；
Figure 11 is the remaining data schematic diagram obtained after bandpass filtering according to embodiments of the present invention；
Figure 12 is the structured flowchart of the compression process device of optical cable coupled noise according to embodiments of the present invention.
Embodiment
It is right with reference to embodiment and accompanying drawing for the object, technical solutions and advantages of the present invention are more clearly understood
The present invention is described in further details.Here, the exemplary embodiment of the present invention and its illustrating to be used to explain the present invention, but simultaneously
It is not as a limitation of the invention.
Inventor considers that optical cable coupled noise is regular stronger, it is easy to statistical simulation, therefore can be by being fitted
To the optical cable coupled noise in vertical seismic data, then noise is suppressed by way of subtraction, so as to solve
In vertical seismic data the compacting problem of interference noise is produced because optic cable can not be close to the borehole wall.
In embodiments of the present invention, it is proposed that a kind of compression process method of optical cable coupled noise, as shown in figure 1, including
Following steps：
Step 101：From the vertical seismic data gathered by distribution type fiber-optic, determine occur optical cable coupled noise
Interference range；
Step 102：It is determined that interference range in read noise periods number and corresponding time span, and calculate and obtain single
Noise periods；
In this step, single noise periods can be calculated according to below equation and obtained：
Wherein, Δ T represents single noise periods, and N represents the noise periods number in the interference range, and T represents N number of noise week
Phase corresponding time span.
Step 103：By the interference range and single noise periods of determination, by the way of mathematical interpolation, obtain non-monotonic
Time curve；
Step 104：According to the nonmonotonic time curve, amplitude attenuation curve is obtained；
Specifically, in this step, can be according to nonmonotonic time curve, the ground of the removal vertical seismic data
Shake sample value corresponding with the non-monotonic time curve in section；Then letter is carried out to removing the seismic profile after sample value again
Number fitting, so as to obtain amplitude attenuation curve；
Step 105：The complete region of lineups waveform is selected from the interference range, according to the amplitude attenuation curve,
Calculating obtains average wavelet；
Step 106：The non-monotonic time curve and the average wavelet are subjected to convolution operation, obtain being fitted noise number
According to；
Step 107：The fitting noisy data is subtracted by the vertical seismic data, remaining data is obtained；
Step 108：The remaining data is filtered, the data after the compacting of optical cable coupled noise are obtained；
The step mainly considers that remaining data is obtained in previous step can have exceptional value, accordingly, it would be desirable to it
It is filtered, while it is also contemplated that in subtraction, indivedual sampling point high-frequency anomalies may be caused, therefore can utilize band logical filter
Ripple carries out exceptional value elimination, to ensure the reliability of data.Specifically, band logical sweep parameter during bandpass filtering, can be by original
Frequency range in beginning data is determined, general main not injure the effective frequency composition in initial data as filtering principle.
In this example, it is proposed that then one kind is entered by being fitted optical cable coupled noise in vertical seismic data using subtraction
The method of row optical cable coupled noise compacting, is produced so as to solve in vertical seismic data because optic cable can not be close to the borehole wall
The compacting problem of interference noise, has reached the technique effect of optical cable coupled noise in effective elimination vertical seismic data, for pushing away
Wide application of the distribution type fiber-optic in geophysical exploration has very important meaning.
The compression process method of above-mentioned optical cable coupled noise is specifically described with reference to a specific embodiment, so
And it is worth noting that, the specific embodiment merely to the present invention is better described, does not constitute the improper limit to the present invention
It is fixed.
Generally speaking, the compression process method of the optical cable coupled noise in this example is gathered according to distribution type fiber-optic (DAS)
Vertical seismic data in optical cable coupled noise development characteristicses and design, mainly include：Set up non-monotonic time curve, meter
Calculate optical cable coupled noise attenuation law, extract optical cable coupled noise wavelet, fitting optical cable coupled noise data and subtraction compacting light
The several aspects of cable coupled noise.
Specifically, as shown in Fig. 2 including：
Step 201：Non-monotonic time curve is set up, including：
S1：Noise jamming area is determined
As shown in figure 3, abscissa is depth (unit：M), ordinate is time (unit：Ms), h1 and h2 is to be interfered
Two interval boundary depths, in the vertical seismic data S that distribution type fiber-optic is gathered, it is first determined occur optical cable coupled noise
Interference range, due to VSP data each trace record one depth of correspondence, therefore represent disturbed area with depth h1 and h2 respectively
Original position and end position.
S2：Cycle determines
As shown in figure 4, in the S1 h1 determined and h2 interval, reading the time span T of N number of noise periods number, calculate
Obtain single noise periods Δ T.It is worth noting that, in this step, periodicity N is under the premise of distinguishable, and value gets over the matter of fundamental importance
The precision of calculation is higher, in Fig. 4, and abscissa is depth (unit：M), ordinate is time (unit：Ms), h1 and h2 be by
To two boundary depths that interference is interval.
S3：Time curve is set up
According to the period Δ T determined in the border h1 and h2 and S2 defined in S1, by mathematical interpolation method, it can obtain
Nonmonotonic time curve (i.e. thinner broken line in Fig. 5) as shown in Figure 5.
Wherein, so-called non-monotonic time curve a, distance (referring to depth in this example) can correspond to multiple moment,
Different from the dull time curve at distance one moment of correspondence on ordinary meaning.
Further, since it is considered that the part before VSP first arrival times is disturbed in the absence of optical cable coupled noise, therefore should
Non-monotonic time curve should cut off the forward part of first arrival time.
Step 202：Calculate optical cable coupled noise attenuation law
According to the non-monotonic time curve obtained in step 201, the sample in the corresponding seismic profile of seismic is taken out
Point value, then does linear function fitting, can obtain amplitude attenuation curve as shown in Figure 6, wherein, so-called sample value refers to
Size of data (i.e. amplitude) in vertical section, any sample value corresponds to a time and a depth.
Step 203：Optical cable coupled noise wavelet is extracted, can be included：
S1：Preferred depth section
In the h1 shown in Fig. 3 and h2 interval, preferably go out the complete region of lineups waveform, avoid because lineups intersect
Superposition area, obtain geological data schematic diagram as shown in Figure 7.This mainly considers that time curve is broken line, in turning point
Annex waveform can produce superposition interference, not possess regularity, it is therefore desirable to avoid this region.
S2：According to the amplitude attenuation curve obtained in step 202, i.e. attenuation law, calculating obtains as shown in Figure 8 be averaged
Wavelet, specifically, is first compensated according to attenuation law so that the sizableness of all wavelet amplitudes for participating in statistics, then
The method for recycling statistical average, calculating obtains average wavelet.
Step 204：Optical cable coupled noise data are fitted, including：
Convolution fortune is carried out according to the average wavelet obtained in the non-monotonic time curve and step 203 obtained in step 201
Calculate, fitting obtains noise data Sn as shown in Figure 9.
Step 205：Subtraction suppresses optical cable coupled noise, can include：
S1：Subtracting techniques suppress noise
Fitting noise number is subtracted using the initial data S (i.e. the vertical seismic data S of cloth collecting fiber) in step 201
According to Sn, remaining data Sr as shown in Figure 10 is obtained；
S2：Eliminate exceptional value
Bandpass filtering is carried out to the remaining data Sr obtained in S1, the remainder after bandpass filtering as shown in figure 11 is obtained
According to Sr ', remaining data Sr ' is the result of final optical cable coupled noise subtraction compacting.
Why consider to eliminate exceptional value by the way of bandpass filtering herein, mainly in view of in subtraction
In, indivedual sampling point high-frequency anomalies may be caused, it is therefore desirable to carry out exceptional value elimination using bandpass filtering, can with ensure data
By property.Wherein, band logical sweep parameter during bandpass filtering, can be determined by the frequency range in initial data, mainly not hinder
And effective frequency composition is filtering principle in initial data.
Above-mentioned steps 201 to 205 are accomplished that the noise compacting to a region to be analyzed, repeat above-mentioned steps
201 to 205 can just realize the compression process to other region optical cable coupled noises.
Why do not use the methods such as pie slice to remove noise data in the method, be primarily due to optical cable coupling
Noise is approached with effective upgoing wave apparent velocity, removes noise effects by these modes in actual applications bad.
Based on same inventive concept, a kind of compression process dress of optical cable coupled noise is additionally provided in the embodiment of the present invention
Put, as described in the following examples.Because the principle that the compression process device of optical cable coupled noise solves problem is coupled with optical cable
The compression process method of noise is similar, thus the compression process device of optical cable coupled noise implementation may refer to optical cable coupling make an uproar
The implementation of the compression process method of sound, repeats part and repeats no more.Used below, term " unit " or " module " can be with
Realize the combination of the software and/or hardware of predetermined function.Although the device described by following examples preferably comes real with software
It is existing, but hardware, or the realization of the combination of software and hardware is also that may and be contemplated.Figure 12 is the embodiment of the present invention
A kind of structured flowchart of the compression process device of optical cable coupled noise, as shown in figure 12, including：Interference range determining module 1201,
It is single noise periods determining module 1202, time curve curve determining module 1203, amplitude attenuation curve determining module 1204, flat
Equal wavelet computing module 1205, fitting noisy data acquisition module 1206, remaining data acquisition module 1207 and filtration module
1208, the structure is illustrated below.
Interference range determining module 1201, for from the vertical seismic data gathered by distribution type fiber-optic, it is determined that setting out
The interference range of third contact of a total solar or lunar eclipse cable coupled noise；
Single noise periods determining module 1202, for it is determined that interference range in read noise periods number and it is corresponding when
Between length, and calculate obtain single noise periods；
Time curve curve determining module 1203, for the interference range by determination and single noise periods, using mathematics
The mode of interpolation, obtains nonmonotonic time curve；
Amplitude attenuation curve determining module 1204, it is bent for according to the nonmonotonic time curve, obtaining amplitude decay
Line；
Average wavelet computing module 1205, for selecting the complete region of lineups waveform from the interference range, according to
The amplitude attenuation curve, calculating obtains average wavelet；
Fitting noisy data acquisition module 1206, for the non-monotonic time curve and the average wavelet to be carried out into pleat
Product computing, obtains fitting noisy data；
Remaining data acquisition module 1207, for subtracting the fitting noisy data by the vertical seismic data, is obtained
To remaining data；
Filtration module 1208, for being filtered to the remaining data, obtains the data after the compacting of optical cable coupled noise.
In one embodiment, amplitude attenuation curve determining module 1204 can include：Removal unit, for according to described
Nonmonotonic time curve, removes sample corresponding with the non-monotonic time curve in the seismic profile of the vertical seismic data
Point value；Amplitude attenuation curve determining unit, for carrying out Function Fitting to removing the seismic profile after sample value, obtains amplitude and declines
Subtract curve.
In one embodiment, single noise periods determining module 1202 is specific can be used for calculating single according to below equation
Individual noise periods：
Wherein, Δ T represents single noise periods, and N represents the noise periods number in the interference range, and T represents N number of noise week
Phase corresponding time span.
In one embodiment, filtration module 1208, specifically can be used for by way of bandpass filtering to the residue
Data are filtered.
In one embodiment, filtration module 1208 is used during bandpass filtering is carried out to the remaining data
Band logical sweep parameter determined by the frequency range of the vertical seismic data.
In another embodiment, a kind of software is additionally provided, the software is used to perform above-described embodiment and preferred real
Apply the technical scheme described in mode.
In another embodiment, a kind of storage medium is additionally provided, be stored with above-mentioned software in the storage medium, should
Storage medium includes but is not limited to：CD, floppy disk, hard disk, scratch pad memory etc..
As can be seen from the above description, the embodiment of the present invention realizes following technique effect：Hung down by fitting
Optical cable coupled noise in straight geological data, is then suppressed noise by way of subtraction, vertical so as to solve
The compacting problem of interference noise is produced because optic cable can not be close to the borehole wall in geological data, effective elimination has been reached vertically
Shake the technique effect of the optical cable coupled noise in data.
Obviously, those skilled in the art should be understood that each module or each step of the above-mentioned embodiment of the present invention can be with
Realized with general computing device, they can be concentrated on single computing device, or be distributed in multiple computing devices
On the network constituted, alternatively, the program code that they can be can perform with computing device be realized, it is thus possible to by it
Store and performed in the storage device by computing device, and in some cases, can be to be held different from order herein
They, are either fabricated to each integrated circuit modules or will be multiple in them by the shown or described step of row respectively
Module or step are fabricated to single integrated circuit module to realize.So, the embodiment of the present invention is not restricted to any specific hard
Part and software are combined.
The preferred embodiments of the present invention are the foregoing is only, are not intended to limit the invention, for the skill of this area
For art personnel, the embodiment of the present invention can have various modifications and variations.Within the spirit and principles of the invention, made
Any modification, equivalent substitution and improvements etc., should be included in the scope of the protection.
Claims (10)
1. a kind of compression process method of optical cable coupled noise, it is characterised in that including：
From the vertical seismic data gathered by distribution type fiber-optic, determine occur the interference range of optical cable coupled noise；
It is determined that interference range in read noise periods number and corresponding time span, and calculate and obtain single noise periods；
By the interference range and single noise periods of determination, by the way of mathematical interpolation, nonmonotonic time curve is obtained；
According to the nonmonotonic time curve, amplitude attenuation curve is obtained；
The complete region of lineups waveform is selected from the interference range, according to the amplitude attenuation curve, calculating is averaged
Wavelet；
The non-monotonic time curve and the average wavelet are subjected to convolution operation, fitting noisy data is obtained；
The fitting noisy data is subtracted by the vertical seismic data, remaining data is obtained；
The remaining data is filtered, the data after the compacting of optical cable coupled noise are obtained.
2. according to the method described in claim 1, it is characterised in that according to the nonmonotonic time curve, obtain amplitude and decline
Subtract curve, including：
According to the nonmonotonic time curve, remove in the seismic profile of the vertical seismic data with it is described non-monotonic when away from
The corresponding sample value of curve；
Function Fitting is carried out to removing the seismic profile after sample value, amplitude attenuation curve is obtained.
3. according to the method described in claim 1, it is characterised in that calculate single noise periods according to below equation：
<mrow>
<mi>&Delta;</mi>
<mi>T</mi>
<mo>=</mo>
<mfrac>
<mi>T</mi>
<mi>N</mi>
</mfrac>
</mrow>
Wherein, Δ T represents single noise periods, and N represents the noise periods number in the interference range, and T represents N number of noise periods pair
The time span answered.
4. according to the method described in claim 1, it is characterised in that the remaining data is filtered, including：
Bandpass filtering is carried out to the remaining data.
5. method according to claim 4, it is characterised in that during bandpass filtering is carried out to the remaining data
The band logical sweep parameter used is determined by the frequency range of the vertical seismic data.
6. a kind of compression process device of optical cable coupled noise, it is characterised in that including：
Interference range determining module, occurs optical cable coupling for from the vertical seismic data gathered by distribution type fiber-optic, determining
Close the interference range of noise；
Single noise periods determining module, for it is determined that interference range in read noise periods number and corresponding time span,
And calculating obtains single noise periods；
Time curve curve determining module, for the interference range by determination and single noise periods, using the side of mathematical interpolation
Formula, obtains nonmonotonic time curve；
Amplitude attenuation curve determining module, for according to the nonmonotonic time curve, obtaining amplitude attenuation curve；
Average wavelet computing module, for selecting the complete region of lineups waveform from the interference range, according to the amplitude
Attenuation curve, calculating obtains average wavelet；
Fitting noisy data acquisition module, for the non-monotonic time curve and the average wavelet to be carried out into convolution operation,
Obtain fitting noisy data；
Remaining data acquisition module, for subtracting the fitting noisy data by the vertical seismic data, obtains remainder
According to；
Filtration module, for being filtered to the remaining data, obtains the data after the compacting of optical cable coupled noise.
7. device according to claim 6, it is characterised in that amplitude attenuation curve determining module includes：
Removal unit, for according to the nonmonotonic time curve, remove in the seismic profile of the vertical seismic data with
The corresponding sample value of the non-monotonic time curve；
Amplitude attenuation curve determining unit, for carrying out Function Fitting to removing the seismic profile after sample value, obtains amplitude and declines
Subtract curve.
8. device according to claim 6, it is characterised in that the single noise periods determining module specifically for according to
Below equation calculates single noise periods：
<mrow>
<mi>&Delta;</mi>
<mi>T</mi>
<mo>=</mo>
<mfrac>
<mi>T</mi>
<mi>N</mi>
</mfrac>
</mrow>
Wherein, Δ T represents single noise periods, and N represents the noise periods number in the interference range, and T represents N number of noise periods pair
The time span answered.
9. device according to claim 6, it is characterised in that the filtration module, specifically for passing through bandpass filtering
Mode is filtered to the remaining data.
10. device according to claim 9, it is characterised in that the filtration module is carrying out band to the remaining data
The band logical sweep parameter used during pass filter is determined by the frequency range of the vertical seismic data.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510453845.0A CN105158801B (en) | 2015-07-29 | 2015-07-29 | The compression process method and apparatus of optical cable coupled noise |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510453845.0A CN105158801B (en) | 2015-07-29 | 2015-07-29 | The compression process method and apparatus of optical cable coupled noise |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105158801A CN105158801A (en) | 2015-12-16 |
| CN105158801B true CN105158801B (en) | 2017-10-17 |
Family
ID=54799716
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510453845.0A Active CN105158801B (en) | 2015-07-29 | 2015-07-29 | The compression process method and apparatus of optical cable coupled noise |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105158801B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106842324B (en) * | 2016-12-08 | 2018-10-30 | 西安交通大学 | A kind of optical cable coupled noise drawing method based on the sparse optimization of time-frequency domain |
| CN107436451B (en) * | 2017-07-26 | 2019-10-11 | 西安交通大学 | A kind of amplitude spectral method of automatic calculating seismic data optical cable coupled noise degree of strength |
| CN111965700A (en) * | 2019-05-20 | 2020-11-20 | 中国石油天然气集团有限公司 | Method and system for eliminating zero wave number noise in optical fiber acoustic sensing seismic data |
| CN111965704B (en) * | 2019-05-20 | 2023-09-26 | 中国石油天然气集团有限公司 | Optical cable secondary oscillation noise suppression method and system for well seismic data |
| US11340371B2 (en) | 2019-07-11 | 2022-05-24 | Halliburton Energy Services, Inc. | Fiber optic noise mitigation for distributed acoustic sensing |
| CN116413803A (en) * | 2021-12-29 | 2023-07-11 | 中国石油天然气集团有限公司 | Method and device for improving DAS signal-to-noise ratio by utilizing local FK conversion |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6041018A (en) * | 1997-11-13 | 2000-03-21 | Colorado School Of Mines | Method for correcting amplitude and phase differences between time-lapse seismic surveys |
| CN102890291A (en) * | 2012-10-09 | 2013-01-23 | 中国海洋石油总公司 | Self-adapting suppressing method for external source interference in seismic data |
| CN103558634A (en) * | 2013-10-29 | 2014-02-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method and device for obtaining frequency attenuation gradient of seismic data |
| CN103954998A (en) * | 2014-05-09 | 2014-07-30 | 成都晶石石油科技有限公司 | Residual amplitude compensating method based on AVO |
| CN104656140A (en) * | 2014-12-01 | 2015-05-27 | 中国石油天然气集团公司 | Median filtering method and system based on vertical earthquake attenuation laws |
-
2015
- 2015-07-29 CN CN201510453845.0A patent/CN105158801B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6041018A (en) * | 1997-11-13 | 2000-03-21 | Colorado School Of Mines | Method for correcting amplitude and phase differences between time-lapse seismic surveys |
| CN102890291A (en) * | 2012-10-09 | 2013-01-23 | 中国海洋石油总公司 | Self-adapting suppressing method for external source interference in seismic data |
| CN103558634A (en) * | 2013-10-29 | 2014-02-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method and device for obtaining frequency attenuation gradient of seismic data |
| CN103954998A (en) * | 2014-05-09 | 2014-07-30 | 成都晶石石油科技有限公司 | Residual amplitude compensating method based on AVO |
| CN104656140A (en) * | 2014-12-01 | 2015-05-27 | 中国石油天然气集团公司 | Median filtering method and system based on vertical earthquake attenuation laws |
Non-Patent Citations (2)
| Title |
|---|
| "Walkaway_VSP技术及其实际应用";陈沅忠 等;《石油科技论坛》;20120229(第2期);第16-20页 * |
| "小波包分析法压制井间地震资料噪声";王立华 等;《勘探地球物理进展》;20050228;第28卷(第1期);第32-35页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105158801A (en) | 2015-12-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105158801B (en) | The compression process method and apparatus of optical cable coupled noise | |
| US20160252651A1 (en) | Trace downsampling of distributed acoustic sensor data | |
| CN101551465B (en) | Method for adaptively recognizing and eliminating seismic exploration single-frequency interference | |
| CN106707342B (en) | Common-shot-gather multistage surface wave pressing method and device | |
| CN103744114B (en) | Method and device for estimating quality factor on basis of zero offset VSP (vertical seismic profile) data | |
| CN109239780A (en) | Based on the synchronous method for squeezing wavelet transformation removal surface wave | |
| CN104570116A (en) | Geological marker bed-based time difference analyzing and correcting method | |
| CN106353817A (en) | Seismic exploration data multiple automatic search and suppression method | |
| CN103135133A (en) | Method and device of vector noise reduction of multi-component seismic data | |
| CN113189641B (en) | Two-channel multi-mode Rayleigh wave underground detection system and method | |
| CN112904415B (en) | Q-layer chromatography inversion method, system and processing equipment based on reconstructed source spectrum consistency | |
| CN102323618A (en) | Coherent noise suppression method based on fractional order Fourier transformation | |
| AU2012227284B2 (en) | Methods and systems for attenuating noise generated at fixed locations | |
| WO2021006948A1 (en) | Fiber optic noise mitigation for distributed acoustic sensing | |
| CN109212609A (en) | Near surface Noise Elimination method based on wave equation continuation | |
| CN113219539B (en) | Method and device for removing borehole wave interference in optical fiber acoustic sensing seismic data | |
| CN105301656A (en) | Method and device for suppressing ghost reflection signals | |
| CN112782763B (en) | Seismic quality factor estimation method, device, equipment and storage medium | |
| CN111505707B (en) | Method for extracting dispersion curve from vertical seismic profile data | |
| Feng et al. | Transform learning in the synchrosqueezing frequency domain—A novel denoising strategy for optical fiber seismic records | |
| CN106338769A (en) | Seismic data denoising method and system | |
| CN105866838B (en) | A kind of seismic data low-frequency information compensation method and device | |
| CN115236733A (en) | DAS-VSP data background noise suppression method based on deep learning | |
| CN105425297A (en) | Method and device of suppressing ghost signals | |
| CN112415601A (en) | Method and device for determining surface quality factor Q value |
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 |