CN104111477B - Method and system for acquiring sliding time for harmonic suppression - Google Patents

Abstract

The invention discloses a method and a system for acquiring sliding time for harmonic suppression. The method comprises the following steps: the start frequency, the stop frequency and the scanning length of a vibroseis scanning signal are set, a plurality of pieces of vibration amplitude data are acquired by a seismic data acquisition device according to the signal; one piece of vibration amplitude data near a seismic source is selected for time-frequency conversion to generate a time-frequency correspondence, and the time-frequency correspondence is converted into a time-energy correspondence between a fundamental wave and each harmonic; the minimum sliding time is set, calculation is performed according to the start frequency, the stop frequency and the scanning length of the vibroseis scanning signal to obtain interference positions of harmonics of next shot to the fundamental wave of the current shot; the energy of the fundamental wave and the energy of each harmonic in each interference position are acquired according to the time-energy correspondence between the fundamental wave and the harmonic, and calculation is carried out to obtain the percentage of the energy of each harmonic from the energy of the fundamental wave; the harmonics of which the percentages exceeding a limit value are selected and suppressed, and the sliding time for harmonic suppression is calculated.

Description

A kind of sliding time acquisition methods for suppressing harmonic wave and system

Technical field

The present invention relates to petroleum exploration field, when being a kind of seismic prospecting to use controlled source slip scan method production The sliding time acquisition methods of compacting harmonic wave and system.

Background technology

At present, in oil seismic exploration, slip sweep is to become most in controlled source highly effective gathering technology in the world One of ripe, the most frequently used technology, the popularization and application the most at home of this technology also obtain good exploration effects and economic benefit.? In the application of controlled source slip sweep, harmonic wave interference is objective reality, needs special compact technique during indoor process, And can also effectively weaken harmonic wave by sliding time reasonable in design in the wild and disturb, and effect is than indoor pressing result more Good.Therefore, raw data quality when sliding time not only affects field process, or affect the most critical parameters of production efficiency, By quick and precisely choosing the sliding time parameter of beneficially harmonic wave compacting, to improving production efficiency, improving data quality has Significance.But, it is the heaviest in controlled source slip scan work progress for choosing sliding time parameter the most rapidly and accurately One of work wanted.

Summary of the invention

For choosing sliding time parameter the most rapidly and accurately, the present invention is by a kind of sliding time choosing method, knot Close this sliding time and can the most effectively carry out harmonic wave compacting, thus improve production efficiency, improve data quality.

For reaching above-mentioned purpose, the present invention proposes a kind of sliding time acquisition methods for suppressing harmonic wave, including with Lower step: set the initial frequency of controlled source scanning signal, terminate frequency, sweep length, and sweep according to described controlled source Retouching signal utilizes seismic data acquisition device to obtain multiple tracks Oscillation Amplitude data；Described multiple tracks Oscillation Amplitude data are chosen and leans on One Oscillation Amplitude data near earthquake source carry out time-frequency convert, generate corresponding T/F corresponding relation, and by described T/F corresponding relation is converted to the Time-energy corresponding relation of first-harmonic and each harmonic；Set the minimum slip time, Calculate according to described the minimum slip time, the initial frequency of described controlled source scanning signal, termination frequency, sweep length and obtain To each harmonic of next big gun each interference position to this big gun first-harmonic；The energy of the first-harmonic and each harmonic that set next big gun divides The most equal with the energy of this big gun first-harmonic corresponding and each harmonic, the Time-energy according to described first-harmonic and each harmonic is corresponding Relation, obtains first-harmonic energy and each harmonic energy of each interference position described, obtains each harmonic energy by calculating Account for the percentage ratio of the first-harmonic energy of correspondence；The harmonic wave selecting percentage ratio to exceed a certain limit value is suppressed, and calculates compacting institute State the sliding time of harmonic wave.

For reaching above-mentioned purpose, the present invention proposes a kind of sliding time for suppressing harmonic wave and obtains system, including: ground Shake data acquisition module, for setting the initial frequency of controlled source scanning signal, terminating frequency, sweep length, and according to institute State controlled source scanning signals collecting and obtain multiple tracks Oscillation Amplitude data；Corresponding relation modular converter, for shaking in described multiple tracks Dynamic amplitude data chooses one Oscillation Amplitude data near focus and carries out time-frequency convert, generate corresponding time m-frequency Rate corresponding relation, and described T/F corresponding relation is converted to the Time-energy corresponding relation of first-harmonic and each harmonic； Interference position acquisition module, for setting the minimum slip time, according to described the minimum slip time, the scanning of described controlled source The initial frequency of signal, terminate frequency, sweep length calculates and to obtain each harmonic of next big gun dry to each of this big gun first-harmonic Disturb position；Energy computation module, for set the energy of the first-harmonic of next big gun and each harmonic respectively with this big gun first-harmonic corresponding And the energy of each harmonic is equal, according to the Time-energy corresponding relation of described first-harmonic and each harmonic, obtain described each and do Disturb first-harmonic energy and each harmonic energy of position, obtain the hundred of the first-harmonic energy that each harmonic energy accounts for correspondence by calculating Proportion by subtraction；Sliding time computing module, suppresses for the harmonic wave selecting percentage ratio to exceed a certain limit value, and calculates compacting institute State the sliding time of harmonic wave.

Sliding time acquisition methods and the system for suppressing harmonic wave of the present invention, it is achieved that controlled source slip scan side Method construction in sliding time quick, accurately select, decrease field test workload, reduce field test expense, shorten Test period, effectively suppress harmonic wave interference, it is ensured that the reasonability that sliding time parameter is chosen.

Accompanying drawing explanation

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 flow chart of steps of the sliding time acquisition methods for suppressing harmonic wave of one embodiment of the invention.

Fig. 2 is the structural representation of the system that obtains for suppressing the sliding time of harmonic wave of one embodiment of the invention.

Fig. 3 is the F-T time-frequency figure of the present invention one specific embodiment.

Fig. 4 is first-harmonic and the Time-energy figure of each harmonic of the present invention one specific embodiment.

Fig. 5 is each harmonic of the present invention one specific embodiment schematic diagram to first-harmonic interference position.

Fig. 6 be the specific embodiment of the invention somewhere slip scan at the beginning of folded generalized section.

Detailed description of the invention

Hereinafter coordinate graphic and presently preferred embodiments of the present invention, the present invention is expanded on further for reaching predetermined goal of the invention institute The technological means taked.

Fig. 1 is the flow chart of steps of the sliding time acquisition methods for suppressing harmonic wave of one embodiment of the invention.Such as Fig. 1 Shown in, comprise the following steps:

Step S101, sets the initial frequency of controlled source scanning signal, terminates frequency, sweep length, and according to controlled Focus scanning signal utilizes seismic data acquisition device to obtain multiple tracks Oscillation Amplitude data；

In this step, the data acquisition unit used can be engineering seismograph conventional in existing seismic prospecting, inspection Ripple device, the equipment such as focus generating means.

Step S102, chooses one Oscillation Amplitude data near focus in multiple tracks Oscillation Amplitude data and carries out time-frequency and turn Change, generate corresponding T/F corresponding relation, and T/F corresponding relation is converted to first-harmonic and each harmonic Time-energy corresponding relation；

In this step, one Oscillation Amplitude data within focus 100 meters or 200 meters typically can be chosen, because Interference size is relevant with the distance of interference source, the most closely disturbs the strongest.

It addition, in this step, the T/F corresponding relation of generation can be represented by T/F figure；First-harmonic And the Time-energy corresponding relation of each harmonic can be represented by Time-energy figure；It is to say, time-frequency figure, time m- Energy diagram is a kind of form of expression of time-frequency corresponding relation, Time-energy corresponding relation respectively.

In this step, first, the T/F figure that time-frequency convert obtains is utilized, then, according to the wherein corresponding time With frequency abstraction first-harmonic and the energy peak of each harmonic, the Time-energy figure of first-harmonic and each harmonic can be obtained.

Step S103, sets the minimum slip time, according to the minimum slip time, the initial frequency of controlled source scanning signal Rate, termination frequency, sweep length, be calculated each harmonic each interference position to this big gun first-harmonic of next big gun in theory Put；

In the present embodiment, generally also can set the record time when step S101, can by this record time as The minimum slip time, Theoretical Calculation is gone out next big gun each harmonic interference position to this big gun first-harmonic, then is obtained by subsequent step The size of interference decides whether to avoid this subharmonic.

Step S104, the energy setting the first-harmonic of next big gun and each harmonic is humorous with this big gun first-harmonic corresponding and each time respectively The energy of ripple is equal, according to the Time-energy corresponding relation of first-harmonic and each harmonic, obtains the first-harmonic energy of each interference position And each harmonic energy, by calculating the percentage ratio obtaining the first-harmonic energy that each harmonic energy accounts for correspondence；

In this step, owing to before and after slip scan, two emplacements are put adjacent, shooting parameter is consistent, it is possible to think adjacent The harmonic wave of two big guns and first-harmonic energy are substantially suitable, namely the energy of the first-harmonic of next big gun and each harmonic respectively with corresponding The energy of big gun first-harmonic and each harmonic is equal.

In this step, owing to first-harmonic energy and harmonic energy can not the most simply be divided by calculating percentage ratio, institute Can first it be normalized, obtain each interference position normalized first-harmonic energy value and normalized respectively Subharmonic energy value, calculates the percentage ratio of harmonic wave and first-harmonic energy the most again；Wherein,

The normalization formula utilized is:

$\mathrm{SPL}=20\lg \frac{P_e}{P_{\mathrm{ref}}};$

Wherein, SPL is first-harmonic energy or each harmonic energy of each interference position；

P_refValue is 1；

P_eFor normalized first-harmonic energy value or normalized each harmonic energy value.

Step S105, the harmonic wave selecting percentage ratio to exceed a certain limit value is suppressed, and is calculated the slip of compacting harmonic wave Time.

Wherein, a certain limit value is 5%, 10% or 20%, can be adjusted according to actual needs.

In this step, the formula of the sliding time utilization calculating compacting harmonic wave is:

$T_i=\frac{(k-1)f_E}{k(f_E-f_B)}T;$

Wherein, T_iFor suppressing the sliding time of harmonic wave；

K is the ripple of the harmonic wave exceeding a certain restriction percentage ratio；

f_BInitial frequency for controlled source scanning signal；

f_ETermination frequency for controlled source scanning signal；

T is the sweep length of controlled source scanning signal.

In actual applications, can be to the sliding time T of compacting harmonic wave_iWhen obtaining the slip that rounds after rounding up Between.

Based on same inventive concept, the embodiment of the present invention additionally provides a kind of sliding time for suppressing harmonic wave and obtains System, as described in the following examples.Owing to the principle of this acquisition system solution problem is similar to this acquisition methods, thus be accordingly used in The enforcement of the sliding time acquisition system of compacting harmonic wave may refer to the enforcement of the sliding time acquisition methods for suppressing harmonic wave, Repeat no more in place of repetition.Used below, term " unit " or " module " can realize predetermined function software and/or The combination of hardware.Although the device described by following example preferably realizes with software, but hardware, or software is with hard The realization of the combination of part also may and be contemplated.

Fig. 2 is the structural representation of the system that obtains for suppressing the sliding time of harmonic wave of the embodiment of the present invention, such as Fig. 2 institute Show, including: earthquake data acquisition module 11, corresponding relation modular converter 12, interference position acquisition module 13, energy computation module 14, sliding time computing module 15, is specifically described this structure below.

Earthquake data acquisition module 11, for setting the initial frequency of controlled source scanning signal, terminating frequency, scanning length Degree, and obtain multiple tracks Oscillation Amplitude data according to controlled source scanning signals collecting；

Wherein, earthquake data acquisition module 11 can according to controlled source scanning signals collecting acquisition multiple tracks Oscillation Amplitude data To utilize the equipment such as engineering seismograph conventional in existing seismic prospecting, cymoscope, focus generating means.

Corresponding relation modular converter 12, for choosing one Oscillation Amplitude near focus in multiple tracks Oscillation Amplitude data Data carry out time-frequency convert, generate corresponding T/F corresponding relation, and are converted to by T/F corresponding relation The Time-energy corresponding relation of first-harmonic and each harmonic；

Interference position acquisition module 13, for setting the minimum slip time, sweeps according to the minimum slip time, controlled source Retouch the initial frequency of signal, terminate frequency, sweep length and calculate and obtain each harmonic of next big gun to each of this big gun first-harmonic Interference position；

Energy computation module 14, for set the energy of the first-harmonic of next big gun and each harmonic respectively with this big gun base corresponding The energy involving each harmonic is equal, according to the Time-energy corresponding relation of first-harmonic and each harmonic, obtains each interference position First-harmonic energy and each harmonic energy, obtain each harmonic energy account for the percentage ratio of first-harmonic energy of correspondence by calculating；

Sliding time computing module 15, suppresses for the harmonic wave selecting percentage ratio to exceed a certain limit value, and calculates The sliding time of compacting harmonic wave.

In the present embodiment, corresponding relation modular converter 12 can be additionally used according to T/F corresponding relation, according to wherein Corresponding time and frequency abstraction first-harmonic and the energy peak of each harmonic, generate first-harmonic and the Time-energy pair of each harmonic Should be related to.

In the present embodiment, it is possible to computing module 14 sets the energy of the first-harmonic of next big gun and each harmonic respectively with phase This big gun first-harmonic answered and the energy of each harmonic equal be owing to this big gun is adjacent with the position of next big gun, setup parameter is consistent.

In the present embodiment, it is possible to computing module 14 is additionally operable to the first-harmonic energy according to each interference position and each time humorous Wave energy, is normalized it, obtain each interference position normalized first-harmonic energy value and normalized respectively Subharmonic energy value；Wherein,

The normalization formula utilized is:

$\mathrm{SPL}=20\lg \frac{P_e}{P_{\mathrm{ref}}};$

Wherein, SPL is first-harmonic energy or each harmonic energy of each interference position；

P_refValue is 1；

P_eFor normalized first-harmonic energy value or normalized each harmonic energy value.

In the present embodiment, a certain limit value in sliding time computing module 15 is 5%, 10% or 20%.

In the present embodiment, sliding time computing module 15, enters for the harmonic wave selecting percentage ratio to exceed a certain limit value Row is suppressed, and the formula that the sliding time calculating compacting harmonic wave utilizes is:

$T_i=\frac{(k-1)f_E}{k(f_E-f_B)}T;$

Wherein, T_iFor suppressing the sliding time of harmonic wave；

K is the ripple of the harmonic wave exceeding a certain restriction percentage ratio；

f_BInitial frequency for controlled source scanning signal；

f_ETermination frequency for controlled source scanning signal；

T is the sweep length of controlled source scanning signal.

When reality is applied, can be to the sliding time T of compacting harmonic wave_iWhen obtaining the slip that rounds after rounding up Between.

In order to the above-mentioned sliding time acquisition methods for suppressing harmonic wave and system are carried out apparent explanation, below Illustrate in conjunction with a specific embodiment, however, it should be noted that this embodiment is merely to be better described this Bright, it is not intended that and the present invention is limited improperly.

1, first, start-stop frequency f of controlled source scanning signal is set_B-f_E(Hz), sweep length T (s), it is also possible to set Just listen the time (t).Concrete value is: start-stop frequency 6-84Hz, sweep length 14s.

2, scene gathers controlled source initial data, these data according to this signal by seismic data acquisition device in the wild Include multiple tracks Oscillation Amplitude data.

3, choose a track data near earthquake source and carry out time-frequency convert and generate T/F figure, as it is shown on figure 3, for this tool The F-T time-frequency figure of body embodiment.Figure can recognize that 9 subharmonic, according to the characteristic of focus, 13 subharmonic can be produced in theory, But higher hamonic wave energy is the most weak.

4, then, this time-frequency figure is changed, according to wherein corresponding time and frequency abstraction first-harmonic and each harmonic Energy peak, generate first-harmonic and the Time-energy figure of each harmonic.Wherein,

First-harmonic start-stop frequency 6-84Hz, second harmonic start-stop frequency is 12-168Hz, and triple-frequency harmonics start-stop frequency is 18- 252Hz, other subharmonic start-stop frequencies are multiplied by number of times by fundamental frequency successively；

As shown in Figure 4, for first-harmonic and the Time-energy figure of each harmonic of this specific embodiment, figure has been got 5 times humorous Ripple.Owing to harmonic energy gradually reduces, therefore it is to separate on this figure.

5, the minimum slip time is set as 6s, in conjunction with start-stop frequency 6-84Hz of controlled source scanning signal, sweep length 14s, when the minimum slip time 6s, each harmonic is to first-harmonic interference position, as it is shown in figure 5, be this specific embodiment for Theoretical Calculation The each harmonic schematic diagram to first-harmonic interference position.

Wherein, this big gun sends when 0s, and its first-harmonic start-stop frequency is 6-84hz, and sweep length is 14s, is also recorded in figure The second harmonic of this big gun, triple-frequency harmonics, four-time harmonic；Next big gun sends at 6s, have recorded the first-harmonic of next big gun, secondary humorous in figure Ripple, triple-frequency harmonics, four-time harmonic.

So, the two of next big gun, three, four-time harmonic be second harmonic interference position to the interference position of this big gun first-harmonic 10.8s, third harmonic interference position 8s, four-time harmonic interference position 7s.

6, due to before and after slip scan two emplacements put adjacent, shooting parameter is consistent, it is taken as that adjacent big gun harmonic wave and first-harmonic Energy is the most suitable.So, in conjunction with the Time-energy figure shown in Fig. 4, it is possible to obtain each harmonic that interference position is corresponding and base Wave energy.

Wherein, the decibel value of second harmonic interference position ENERGY E 2 (value position is 10.8-6=4.8s) is-20dB, right The decibel value answering first-harmonic ENERGY E 12 (value position 10.8s) is-10dB.

The decibel value of third harmonic interference position ENERGY E 3 (value position 8-6=2s) is-35dB, corresponding first-harmonic energy The decibel value of E13 (value position 8s) is-9.5dB.

The decibel value of four-time harmonic interference position ENERGY E 4 (value position 7-6=1s) is-45dB, corresponding first-harmonic energy The decibel value of E14 (value position 7s) is-6dB.

7, respectively above decibel value is normalized；

The normalized value of calculating second harmonic interference position ENERGY E 2:Wherein, SPL=-20, P_ref=1, can be in the hope of P_e=0.1；

So, successively can be in the hope of: P (E2)=0.1；P (E12)=0.316；P (E3)=0.017783；P (E13)= 0.334965；P (E4)=0.0056；P (E14)=0.501.

Then the percentage ratio that harmonic energy each time accounts for the first-harmonic energy of its correspondence is asked for:

$R2=\frac{P\left(E2\right)}{P\left(E12\right)}\×100\%=31.6\%;$

$R3=\frac{P\left(E3\right)}{P\left(E13\right)}\×100\%=5.3\%;$

$R4=\frac{P\left(E4\right)}{P\left(E14\right)}\×100\%=1.1\%.$

8, according to practical operation requirement, the harmonic wave selecting more than 10% is suppressed, so this time selecting compacting secondary humorous Ripple.In conjunction with computing formulaWherein, k=2, T=14, f_B=6, f_E=84；It is calculated T₂=7.5s； For 8s after rounding up, so this time desirable 8s of sliding time, the impact of second harmonic interference can be avoided in the wild.

If in practical operation, will select to suppress two, third harmonic interference, then k desirable 3, are calculated T₃> 10.05s, takes sliding time 11s.

In one embodiment, sliding time acquisition methods and the system slip scan examination in somewhere of the present invention are utilized The sliding time tested has carried out Preliminary Applications on selecting.Expert advice sliding time takes 14s, and reason is that fearness is domestic removes harmonic management Technology does not passes a test, and affects live data quality；The suggestion of chief engineer is that sliding time takes 10s, brings into and serves harmonic wave interference, test GeoEast processes noise removal capability, and facilitated technique research is with progressive.By advising after the technique study of the present invention that sliding time takes 8s, suppresses second harmonic；Or 11s, compacting two, triple-frequency harmonics.Finally using 10s, actual count is the most sliding in the wild The dynamic time is 10.88s, and major part big gun avoids two, triple-frequency harmonics.Fig. 6 is the somewhere slip scan of the specific embodiment of the invention The most folded generalized section, from the point of view of Fig. 6, after adjusting sliding time, substantially avoided harmonic wave interference so that harmonic wave obtains Good compacting.

Sliding time acquisition methods and the system for suppressing harmonic wave of the present invention, it is achieved that controlled source slip scan side Method construction in sliding time quick, accurately select, decrease field test workload, reduce field test expense, shorten Test period, effectively suppress harmonic wave interference, it is ensured that the reasonability that sliding time parameter is chosen.

Particular embodiments described above, has been carried out the purpose of the present invention, technical scheme and beneficial effect the most in detail Describe in detail bright, be it should be understood that the specific embodiment that the foregoing is only the present invention, the guarantor being not intended to limit the present invention Protect scope, all within the spirit and principles in the present invention, any modification, equivalent substitution and improvement etc. done, should be included in this Within the protection domain of invention.

Claims (10)

1. the sliding time acquisition methods being used for suppressing harmonic wave, it is characterised in that comprise the following steps:

Set the initial frequency of controlled source scanning signal, terminate frequency, sweep length, and according to described controlled source scanning letter Seismic data acquisition device number is utilized to obtain multiple tracks Oscillation Amplitude data；

In described multiple tracks Oscillation Amplitude data, choose one Oscillation Amplitude data near focus carry out time-frequency convert, generate with The T/F corresponding relation of its correspondence, according to described T/F corresponding relation, according to wherein corresponding time and frequency Extract first-harmonic and the energy peak of each harmonic, generate first-harmonic and the Time-energy corresponding relation of each harmonic；

Set the minimum slip time, scan the initial frequency of signal, end according to described the minimum slip time, described controlled source Only frequency, sweep length calculate and obtain each harmonic each interference position to this big gun first-harmonic of next big gun；

The energy setting the first-harmonic of next big gun and each harmonic is equal with the energy of this big gun first-harmonic corresponding and each harmonic respectively, Time-energy corresponding relation according to described first-harmonic and each harmonic, obtain each interference position described first-harmonic energy and Each harmonic energy, according to first-harmonic energy and each harmonic energy of each interference position described, is normalized place to it Reason, obtains the normalized first-harmonic energy value of each interference position described and normalized each harmonic energy value, then passes through Calculate the percentage ratio obtaining the first-harmonic energy that each harmonic energy accounts for correspondence；Wherein,

The normalization formula utilized is:

$SPL=20\lg \frac{P_e}{P_{ref}};$

Wherein, SPL is first-harmonic energy or each harmonic energy of each interference position described；

P_refValue is 1；

P_eFor normalized first-harmonic energy value or normalized each harmonic energy value；

The harmonic wave selecting percentage ratio to exceed a certain limit value is suppressed, and calculates the sliding time suppressing described harmonic wave.

Sliding time acquisition methods for suppressing harmonic wave the most according to claim 1, it is characterised in that under described setting The first-harmonic of one big gun and the energy of each harmonic are equal with the energy of this big gun first-harmonic corresponding and each harmonic respectively to be included:

Owing to this big gun is adjacent with the position of next big gun, setup parameter is consistent, so setting the first-harmonic of next big gun and each harmonic Energy is equal with the energy of this big gun first-harmonic corresponding and each harmonic respectively.

Sliding time acquisition methods for suppressing harmonic wave the most according to claim 1, it is characterised in that described selection hundred Proportion by subtraction exceedes the harmonic wave of a certain limit value and suppresses, and calculates and suppress the sliding time of described harmonic wave and include:

Wherein, described a certain limit value is 5%, 10% or 20%.

Sliding time acquisition methods for suppressing harmonic wave the most according to claim 1, it is characterised in that described selection hundred Proportion by subtraction exceedes the harmonic wave of a certain limit value and suppresses, and the formula calculating the sliding time utilization suppressing described harmonic wave is:

$T_i=\frac{(k-1)f_E}{k(f_E-f_B)}T;$

Wherein, T_iSliding time for described compacting harmonic wave；

K exceedes the ripple time of harmonic wave of a certain restriction percentage ratio described in being；

f_BInitial frequency for described controlled source scanning signal；

f_ETermination frequency for described controlled source scanning signal；

T is the sweep length of described controlled source scanning signal.

Sliding time acquisition methods for suppressing harmonic wave the most according to claim 4, it is characterised in that described selection hundred Proportion by subtraction exceedes the harmonic wave of a certain limit value and suppresses, and calculates and suppress the sliding time of described harmonic wave and also include:

Sliding time T to described compacting harmonic wave_iA sliding time rounded is obtained after rounding up.

6. one kind obtains system for suppressing the sliding time of harmonic wave, it is characterised in that including:

Earthquake data acquisition module, for setting the initial frequency of controlled source scanning signal, terminating frequency, sweep length, and Multiple tracks Oscillation Amplitude data are obtained according to described controlled source scanning signals collecting；

Corresponding relation modular converter, for choosing one Oscillation Amplitude number near focus in described multiple tracks Oscillation Amplitude data According to carrying out time-frequency convert, generate corresponding T/F corresponding relation, according to described T/F corresponding relation, according to Wherein corresponding time and frequency abstraction first-harmonic and the energy peak of each harmonic, generate first-harmonic and the time m-energy of each harmonic Amount corresponding relation；

Interference position acquisition module, for setting the minimum slip time, according to described the minimum slip time, described controlled source The initial frequency of scanning signal, terminate frequency, sweep length calculate and obtain each to this big gun first-harmonic of each harmonic of next big gun Individual interference position；

Energy computation module, for set the energy of the first-harmonic of next big gun and each harmonic respectively with this big gun first-harmonic corresponding and respectively The energy of subharmonic is equal, according to the Time-energy corresponding relation of described first-harmonic and each harmonic, obtains each interference position described The first-harmonic energy put and each harmonic energy, according to first-harmonic energy and each harmonic energy of each interference position described, It is normalized, obtain each interference position described normalized first-harmonic energy value and normalized each time humorous Wave energy value, then by calculating the percentage ratio obtaining the first-harmonic energy that each harmonic energy accounts for correspondence；Wherein,

The normalization formula utilized is:

$SPL=20\lg \frac{P_e}{P_{ref}};$

Wherein, SPL is first-harmonic energy or each harmonic energy of each interference position described；

P_refValue is 1；

P_eFor normalized first-harmonic energy value or normalized each harmonic energy value；

Sliding time computing module, suppresses for the harmonic wave selecting percentage ratio to exceed a certain limit value, and calculates compacting institute State the sliding time of harmonic wave.

Sliding time for suppressing harmonic wave the most according to claim 6 obtains system, it is characterised in that described energy meter Calculate module, set the energy of the first-harmonic of next big gun and each harmonic respectively with this big gun first-harmonic corresponding and the energy phase of each harmonic Etc. including:

Owing to this big gun is adjacent with the position of next big gun, setup parameter is consistent, so setting the first-harmonic of next big gun and each harmonic Energy is equal with the energy of this big gun first-harmonic corresponding and each harmonic respectively.

Sliding time for suppressing harmonic wave the most according to claim 6 obtains system, it is characterised in that during described slip Between computing module, suppress for the harmonic wave that selects percentage ratio to exceed a certain limit value, wherein, described a certain limit value is 5%, 10% or 20%.

Sliding time for suppressing harmonic wave the most according to claim 6 obtains system, it is characterised in that during described slip Between computing module, the harmonic wave exceeding a certain limit value for described selection percentage ratio suppresses, and calculates the described harmonic wave of compacting Sliding time utilize formula be:

$T_i=\frac{(k-1)f_E}{k(f_E-f_B)}T;$

Wherein, T_iSliding time for described compacting harmonic wave；

K exceedes the ripple time of harmonic wave of a certain restriction percentage ratio described in being；

f_BInitial frequency for described controlled source scanning signal；

f_ETermination frequency for described controlled source scanning signal；

T is the sweep length of described controlled source scanning signal.

Sliding time for suppressing harmonic wave the most according to claim 9 obtains system, it is characterised in that described slip Time Calculation module, suppresses for the harmonic wave selecting percentage ratio to exceed a certain limit value, and calculates the described harmonic wave of compacting Sliding time also includes:

Sliding time T to described compacting harmonic wave_iA sliding time rounded is obtained after rounding up.