CN105572723B - The design method of controlled source scanning signal based on autocorrelation wavelet - Google Patents
The design method of controlled source scanning signal based on autocorrelation wavelet Download PDFInfo
- Publication number
- CN105572723B CN105572723B CN201410542447.1A CN201410542447A CN105572723B CN 105572723 B CN105572723 B CN 105572723B CN 201410542447 A CN201410542447 A CN 201410542447A CN 105572723 B CN105572723 B CN 105572723B
- Authority
- CN
- China
- Prior art keywords
- mrow
- scanning signal
- frequency
- wavelet
- controlled source
- 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 present invention provides a kind of design method of the controlled source scanning signal based on autocorrelation wavelet, and this method includes：According to survey area frequency requirement, preferably go out the wavelet that wavelet form is good, frequency band is wider, and design the start-stop frequency, sweep length and start-stop slope length of scanning signal；Ask for wavelet spectrum；It should be allocated according to the sweep time of each frequency with the proportional relation of amplitude of frequency content requirement, obtain the function of time t (f) corresponding to each sample frequency；Function of time t (f) progress inverse transformation is asked for into time-frequency function f (t)；Instantaneous phase is asked for by being integrated to time-frequency function f (t), and then asks for sinusoidal controlled source scanning signal.The controlled source scanning signal that the present invention is produced overcomes that conventional sweep signal subspace wave morphology is poor, " rotatable phase, logarithm segmentation " scanning signal is relatively low in the energy of high band, the problem of single-shot frequency band is narrower is excited, controlled source seismic data quality can be greatly improved.
Description
Technical field
The present invention relates to oil field Seismic Exploration with Vibrator field, specifically refer to a kind of non-linear controlled source scanning letter
Number design method.
Background technology
Controlled source be dispersedly the incoming underground of energy within the long time, then will be scattered with correlation technique length
Scanning signal energy puts together and is collapsed into the short pulse of narrower in width, it be earth reflection and propagation characteristic finishing after
The set of the auto-correlation function (wavelet) of scanning signal, the quality of its autocorrelation wavelet just directly affects controlled source collection data
Quality.For a scanning signal, its auto-correlation function is that the center crest both sides of correlation wavelet must have many ripples
Valve.If these lobes decay it is very slow, then the correlation wavelet of the earthquake record middle-shallow layer reflected signal after correlation this
Class lobe will produce interference to deep reflex signal, that is correlation wavelet center crest both sides lobe plays noise contributions.It is this
The physical factor (vibrations and sound such as power network, wind, vehicle and personnel walking generation) that noise is not due to is caused, but by
The accessory produced in itself in related operation, so it is called correlated noise.The presence of the presence of secondary lobe, particularly the first secondary lobe
It is that we are undesirable, it can directly influence the resolution and detection of adjacent earth formations.
Scanning signal design method is a lot, and the most frequently used linear, nonlinear scanning, wherein linear signal application are the widest
General, it is easy to realize, the energy that each frequency distribution uniform scanning energy mean allocation, cell frequency are occupied is identical, to scanning
Signal component in signal frequency range does not have frequency filtering effect；And nonlinear scanning signal mainly due to stratum to earthquake
The attenuation by absorption of ripple is not linear and puts forward, and main application has two kinds, and exponential function scanning and logarithmic function scanning can
For supplement low frequency or radio-frequency component, but the linearity and non-linearity autocorrelation wavelet secondary lobe used at present is larger, produce compared with
Strong coherent noise, influences acquisition quality.
In order to reduce correlation side lobes, people have studied the scannings such as simulated frequency conversion scanning, array sweeping or pseudo random scanning
Method, but not extensive use.The characteristic of Zhang Hongle signal Analysis and signal correlation wavelet, draws " rotatable phase, logarithm point
Section " scanning signal meets this requirement.Show that " rotatable phase, logarithm segmentation " scanning signal is strictly one by analysis, demonstration
The scanning signal of correlation wavelet feature can be improved by planting, and in the theory analysis stage, (Zhang Hongle compiling one kind improves correlator baud
The scanning signal levied-" rotatable phase, logarithm segmentation " scanning signal geophysical equipments, 1999.08 (3)：17~20).In order to enter
The advantage Zhang Hongle, Wang Meisheng of one step card " rotatable phase, logarithm segmentation " scanning signal are to logarithm fractional scanning signal and general
Logical linear scan signal has carried out the experiment of seismic acquisition data comparison.Result of the test shows that S/N ratio of seismic records is obtained
It is obviously improved (a kind of scanning signal geophysical equipments of improvement correlation wavelet characteristic of Zhang Hongle, Wang Mei life, 2006.8 (S1)：33
~41)." rotatable phase, logarithm segmentation " scanning signal in order to design, east physical prospecting Cao business is auspicious et al. to be proposed a kind of shaping and calculates
Method, its principle designs scanning signal with Ricker wavelet wave character, and signal passes through the change of continuous phase in the design process
Change, the frequency spectrum and Ricker wavelet spectrums for obtaining signal carry out successive ignition fitting and realize that (Cao Wuxiang, Li Xianqing, Guo Hong open can
Control the plastics design method petroleum explorations of focus scanning signal, 2009.11 (6)：611~614).But this method is using scanning
Signal spectrum is consistent with Ricker wavelet spectrums, and its high band energy is relatively low, excites single-shot frequency band narrower, is unfavorable for high-resolution
Seismic prospecting.For this, we have invented a kind of design method of the new controlled source scanning signal based on autocorrelation wavelet, solution
The above of having determined technical problem.
The content of the invention
There is preferable correlation wavelet form it is an object of the invention to provide a kind of scanning signal, controllable shake is greatly improved
The design method of the controlled source scanning signal based on autocorrelation wavelet of focus earthquake data quality.
The purpose of the present invention can be achieved by the following technical measures：Controlled source scanning signal based on autocorrelation wavelet
Design method, should the design method of controlled source scanning signal based on autocorrelation wavelet include：Step 1, according to exploration area
Domain frequency requirement, preferably goes out the wavelet that wavelet form is good, frequency band is wider, and design start-stop frequency, the sweep length of scanning signal
And start-stop slope length；Step 2, wavelet spectrum is asked for；Step 3, should be with the frequency content according to the sweep time of each frequency
It is required that the proportional relation of amplitude be allocated, obtain the function of time t (f) corresponding to each sample frequency；Step 4, will
Function of time t (f) carries out inverse transformation and asks for time-frequency function f (t)；And step 5, asked by being integrated to time-frequency function f (t)
Instantaneous phase is taken, and then asks for sinusoidal controlled source scanning signal.
The purpose of the present invention can be also achieved by the following technical measures：
In step 2, the wavelet spectrum A (f) of scanning signal is obtained by fourier transform according to given wavelet σ：
In step 3, according to the wavelet spectrum curve in step 2, to sweep time of each frequency according to itself and this frequency
The relation that the amplitude of rate component requirements is proportional is allocated：
Dt (f)=kA (f) df (2)
Time of occurrence of each frequency content in scanning signal is sought, function of time t (f) is obtained：
In formula：
T---- controlled source scanning signal length；
K---- calculates use ratio constant；
The initial frequency of f1---- controlled source scanning signals；
F2---- controlled sources scanning signal, which rises, terminates frequency；
A (f) ----wavelet spectrum；
F---- controlled source scanning signal instantaneous frequencys；
Df---- controlled source scanning signal instantaneous frequency differential.
In step 4, curve t (f) during frequency is asked for the time-frequency curve f (t) of constant duration by inverse transformation.
In steps of 5, the phase of scanning signal, the sinusoidal controlled source scanning letter of output are asked for according to time-frequency function f (t)
Number:
In formula：
K=1,2 ..., Ni
In formula：
S (t) ----controlled source scanning signal；
B (t) ----Blacknam slope (Blackman) ramp function；
T1---- controlled source scanning signal the initial segments slope length；
T2---- controlled sources scanning signal terminates section slope length；
Δ t---- controlled source scanning signal sample rates.
The design method of the controlled source scanning signal based on autocorrelation wavelet in the present invention, is established based on correlator
The design cycle of the controlled source scanning signal of ripple, overcome conventional sweep signal subspace wave morphology it is poor, " rotatable phase, logarithm segmentation "
Scanning signal is relatively low in the energy of high band, excites the problem of single-shot frequency band is narrower, and the scanning signal of design has preferable phase
Climax wave morphology, to greatly improve controlled source seismic data quality.
Scanning signal method is designed based on autocorrelation wavelet in the present invention, pickup area can be analyzed to seismic spectrum
It is required that, the design of signal is scanned by preferred corresponding wavelet, the scanning signal autocorrelation wavelet designed has higher
Definition, and can be by the scanning signals of different wavelet design different frequency bands width.
With conventional linear compared with nonlinear scanning signal, scanning signal of the present invention not only has in earth-layer propagation process damages
Small feature, but also the autocorrelation wavelet with minimum dependence edge leaf are lost, therefore, the seismic data tool of signal excitation of the present invention
There is higher energy and signal to noise ratio, so as to greatly improve controlled source seismic data quality.
With the developing and development of controlled source Exploration Domain, survey area is prolonged to more complicated earth's surface, underground condition
Stretch, controlled source acquisition technique also further develops to high accuracy, high-resolution seismic survey, therefore controlled source scanning signal
Quality it is particularly important, it to improve Seismic Exploration Effect there is material impact.
The present invention is a kind of method that use wavelet carries out nonlinear properties design, from the successful application in In A Certain Place of Xinjiang area
From the point of view of, the design method of the controlled source scanning signal based on autocorrelation wavelet has good applicability and application prospect.
Brief description of the drawings
Fig. 1 is a specific embodiment of the design method of the controlled source scanning signal based on autocorrelation wavelet of the present invention
Flow chart；
Fig. 2 is the specific implementation based on Ricker wavelets with the scanning signal of Wide-band Ricker wavelet design of the invention
The signal of example；
Fig. 3 is linear scan, based on three kinds of specific embodiments of scanning signal one of Ricker wavelets and Wide-band Ricker wavelet
Time-frequency curve comparison；
Fig. 4 is linear scan, the specific implementation based on three kinds of scanning signals of Ricker wavelets and Wide-band Ricker wavelet
Example Spectrum Analysis Comparison；
Fig. 5 is linear scan, the specific implementation based on three kinds of scanning signals of Ricker wavelets and Wide-band Ricker wavelet
Example autocorrelation wavelet analysis contrast；
Fig. 6 compiles record for the single-shot solution of In A Certain Place of Xinjiang area linear signal；
Fig. 7 compiles record for the single-shot solution of In A Certain Place of Xinjiang area Ricker wavelets；
Fig. 8 is that In A Certain Place of Xinjiang area Wide-band Ricker wavelet signal single-shot solution compiles record；
Fig. 9 is In A Certain Place of Xinjiang area linear signal, Ricker wavelets and Wide-band Ricker wavelet signal single-shot target zone energy
Analysis；
Figure 10 is In A Certain Place of Xinjiang area linear signal, Ricker wavelets and Wide-band Ricker wavelet signal single-shot target zone noise
Than analysis.
Embodiment
For enable the present invention above and other objects, features and advantages become apparent, it is cited below particularly go out preferable implementation
Example, and coordinate institute's accompanying drawings, it is described in detail below.
It can be seen from controlled source scanning signal principle, the frequency spectrum of scanning signal and the frequency of autocorrelation wavelet are consistent
, therefore the design of signal can be scanned according to correlation wavelet.The design essence of scanning signal is exactly the design of time-frequency curve,
The sinusoidal signal for being available for controlled source to use is calculated further according to nonlinear scanning signal formula.
As shown in figure 1, knots of the Fig. 1 for the design method of the controlled source scanning signal based on autocorrelation wavelet of the present invention
Structure flow chart.
In step 101, according to survey area frequency requirement, preferably go out the wavelet σ that wavelet form is good, frequency band is wider, and set
Count start-stop frequency f1, f2, sweep length T and the start-stop slope length of scanning signal.In one embodiment, swept according to survey area
The requirement of frequency and energy is retouched, Ricker wavelets is have chosen and is grown with Wide-band Ricker wavelet, start-stop frequency 6-84Hz and scanning
26s is spent, start-stop slope is respectively 800ms, 500ms.Flow enters step 102.
In step 102, by the wavelet σ in step 101 by fourier transform, the frequency spectrum A (f) of scanning signal is obtained：
Flow enters step 103.
In step 103, by the spectrum curve in step 102, according to the sweep time of each frequency according to itself and the frequency
The relation that the amplitude of component requirements is proportional is allocated：
Dt (f)=kA (f) df (2)
Time of occurrence of each frequency content in scanning signal is sought, function of time t (f) is obtained：
In formula：
T---- controlled source scanning signal length；
K---- calculates use ratio constant；
The initial frequency of f1---- controlled source scanning signals；
F2---- controlled sources scanning signal, which rises, terminates frequency；
A (f) ----wavelet spectrum；
F---- controlled source scanning signal instantaneous frequencys；
Df---- controlled source scanning signal instantaneous frequency differential；
In one embodiment, according to the sweep time of each frequency should with the frequency content requirement amplitude it is proportional, ask
Time of occurrence of each frequency content in scanning signal is taken, relation during frequency to determine signal.Flow enters step 104.
In step 104, seek the instantaneous frequency in scanning signal, it is exactly t (f) inverse function f (t), be actually by
Above formula solves f (t) expression formula；In one embodiment, relation during frequency is asked for the time-frequency curve of constant duration by conversion
Data.Flow enters step 105.
In step 105, the phase of scanning signal, the sinusoidal controlled source scanning letter of output are asked for according to time-frequency function f (t)
Number:
In formula：
K=1,2 ..., Ni
In formula：
S (t) ----controlled source scanning signal；
B (t) ----Blacknam slope (Blackman) ramp function；
T1---- controlled source scanning signal the initial segments slope length；
T2---- controlled sources scanning signal terminates section slope length；
Δ t---- controlled source scanning signal sample rates；
In one embodiment, as shown in Fig. 2 having drawn based on Ricker wavelets and Wide-band Ricker wavelet, start-stop frequency
6-84Hz, sweep length 26s and start-stop slope are respectively 800ms, 500ms two controlled source scanning signals.
Performance of the present invention is as follows：
According to the requirement of Xinjiang exploration area scan frequency and energy, Ricker wavelets and broadband Ricker have chosen
Ripple, start-stop frequency 6-84Hz and sweep length 26s, start-stop slope are respectively 800ms, 500ms, and signal is carried out with the present invention
Design, implements step such as Fig. 1, the controlled source scanning signal of generation as shown in Fig. 2 and to linear signal, be based on
These three signals carry out time frequency analysis, such as Fig. 3 to Ricker wavelets with Wide-band Ricker wavelet；Carry out spectrum analysis, such as Fig. 4；Enter
Row autocorrelation wavelet is analyzed, such as Fig. 5, and Fig. 3-Fig. 5 shows that the sub- form of scanning signal auto-correlation that the present invention is produced has larger
Change.After the scanning signal test passes that the present invention is produced, it is input in controlled source with regard to operation can be carried out, such as Fig. 6 to Fig. 8 institutes
It is shown as linear signal, is remembered based on the earthquake single-shot obtained by Ricker wavelets and Wide-band Ricker wavelet these three signal scannings
Record, such as Fig. 9, Figure 10 are respectively to carry out the energy of target zone, Analysis signal-to-noise ratio (SNR) to three single-shots that focus is produced.Present invention scanning
The more conventional linear scan signal of signal is compared, the autocorrelation wavelet with minimum dependence edge leaf, the seismic data tool excited
There is higher energy and signal to noise ratio, more conventional signal scanning data is improved.
Claims (4)
1. the design method of the controlled source scanning signal based on autocorrelation wavelet, it is characterised in that autocorrelation wavelet should be based on
The design method of controlled source scanning signal include：
Step 1, according to survey area frequency requirement, select that wavelet form is good, bandwidth wavelet, and design the start-stop of scanning signal
Frequency, sweep length and start-stop slope length；
Step 2, wavelet spectrum is asked for；
Step 3, it should be allocated according to the sweep time of each frequency with the proportional relation of amplitude of frequency content requirement,
Obtain the function of time t (f) corresponding to each sample frequency；
Step 4, function of time t (f) progress inverse transformation is asked for into time-frequency function f (t)；
Step 5, instantaneous phase is asked for by being integrated to time-frequency function f (t), and then asks for sinusoidal controlled source scanning letter
Number；
In step 3, according to the wavelet spectrum curve in step 2, to sweep time of each frequency according to its with the frequency into
The relation for dividing desired amplitude proportional is allocated：
Dt (f)=kA (f) df (2)
Time of occurrence of each frequency content in scanning signal is sought, function of time t (f) is obtained：
<mrow>
<mi>T</mi>
<mo>=</mo>
<mi>k</mi>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mi>f</mi>
<mn>1</mn>
</mrow>
<mrow>
<mi>f</mi>
<mn>2</mn>
</mrow>
</msubsup>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>f</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>t</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>T</mi>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mi>f</mi>
<mn>1</mn>
</mrow>
<mi>f</mi>
</msubsup>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>f</mi>
</mrow>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mi>f</mi>
<mn>1</mn>
</mrow>
<mrow>
<mi>f</mi>
<mn>2</mn>
</mrow>
</msubsup>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>f</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula：
T---- controlled source scanning signal length；
K---- calculates use ratio constant；
The initial frequency of f1---- controlled source scanning signals；
F2---- controlled sources scanning signal, which rises, terminates frequency；
A (f) ----wavelet spectrum；
F---- controlled source scanning signal instantaneous frequencys；
Df---- controlled source scanning signal instantaneous frequency differential.
2. the design method of the controlled source scanning signal according to claim 1 based on autocorrelation wavelet, its feature exists
In in step 2, according to given wavelet σ by fourier transform, obtaining the wavelet spectrum A (f) of scanning signal：
<mrow>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mrow>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</msubsup>
<mi>&sigma;</mi>
<mo>&CenterDot;</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>j</mi>
<mi>&omega;</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
3. the design method of the controlled source scanning signal according to claim 1 based on autocorrelation wavelet, its feature exists
In in step 4, function of time t (f) being asked for into the time-frequency curve f (t) of constant duration by inverse transformation.
4. the design method of the controlled source scanning signal according to claim 1 based on autocorrelation wavelet, its feature exists
In in steps of 5, asking for the phase of scanning signal according to time-frequency function f (t), the sinusoidal controlled source scanning signal of output:
<mrow>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>B</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mo>&lsqb;</mo>
<mn>2</mn>
<mi>&pi;</mi>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>t</mi>
</msubsup>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula：
<mrow>
<mi>B</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<mfrac>
<mi>t</mi>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<mi>t</mi>
<mo>&le;</mo>
<mi>T</mi>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<mi>T</mi>
<mn>1</mn>
<mo><</mo>
<mi>t</mi>
<mo><</mo>
<mi>T</mi>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>T</mi>
<mo>-</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>T</mi>
<mn>2</mn>
<mo>&le;</mo>
<mi>t</mi>
<mo>&le;</mo>
<mi>T</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0.42</mn>
<mo>-</mo>
<mn>0.5</mn>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&pi;</mi>
<mfrac>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>0.08</mn>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mi>&pi;</mi>
<mfrac>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>T</mi>
<mi>i</mi>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula：
S (t) ----controlled source scanning signal；
B (t) ----Blacknam slope (Blackman) ramp function；
T1----- controlled source scanning signal the initial segments slope length；
T2----- controlled sources scanning signal terminates section slope length；
Δ t---- controlled source scanning signal sample rates；
N---- sampling numbers.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410542447.1A CN105572723B (en) | 2014-10-14 | 2014-10-14 | The design method of controlled source scanning signal based on autocorrelation wavelet |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410542447.1A CN105572723B (en) | 2014-10-14 | 2014-10-14 | The design method of controlled source scanning signal based on autocorrelation wavelet |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105572723A CN105572723A (en) | 2016-05-11 |
| CN105572723B true CN105572723B (en) | 2017-09-12 |
Family
ID=55883064
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410542447.1A Active CN105572723B (en) | 2014-10-14 | 2014-10-14 | The design method of controlled source scanning signal based on autocorrelation wavelet |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105572723B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106199684A (en) * | 2016-08-10 | 2016-12-07 | 中国石油天然气集团公司 | Scanning signal acquiring method and device |
| CN108181645A (en) * | 2017-11-24 | 2018-06-19 | 中国石油天然气集团公司 | A kind of controlled source frequency dividing while exciting method, apparatus and system |
| CN109975874B (en) * | 2017-12-28 | 2020-05-19 | 中石化石油工程技术服务有限公司 | Controllable seismic source scanning signal design method based on damping rake wavelets |
| CN112433250B (en) * | 2020-10-30 | 2024-08-27 | 中国石油天然气集团有限公司 | Method and device for separating orthogonal matching before frequency division aliasing scanning correlation |
| CN112764090B (en) * | 2020-12-24 | 2024-04-09 | 中国石油化工集团有限公司 | Seismic source nonlinear scanning signal design method based on forced vibration |
| CN113189642B (en) * | 2021-04-28 | 2024-04-05 | 中国石油化工集团有限公司 | Seismic source linear scanning signal design method based on forced vibration |
| CN114384580B (en) * | 2021-12-31 | 2023-05-02 | 同济大学 | Ideal wavelet customizing method based on controllable focus |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2219856A (en) * | 1986-09-29 | 1989-12-20 | Chevron Res | Method of producing and recording random vibratory seismic data |
| EP0953179B1 (en) * | 1997-11-14 | 2009-04-15 | Western Atlas International, Inc. | Seismic data acquisition and processing using non-linear distortion in a groundforce signal |
| CN101545981A (en) * | 2008-03-28 | 2009-09-30 | 中国石油天然气集团公司 | Method for realizing minimum phase of vibroseis seismic data zero-phase wavelet |
| CN102062872A (en) * | 2009-11-11 | 2011-05-18 | 中国石油天然气集团公司 | Multi-factor scanning method of controllable vibroseis |
| CN102692643A (en) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | Time varying controllable focal force signal deconvolution method |
| CN104035128A (en) * | 2013-03-06 | 2014-09-10 | 中国石油集团东方地球物理勘探有限责任公司 | Vibroseis pseudo random scanning signal generating method |
-
2014
- 2014-10-14 CN CN201410542447.1A patent/CN105572723B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2219856A (en) * | 1986-09-29 | 1989-12-20 | Chevron Res | Method of producing and recording random vibratory seismic data |
| EP0953179B1 (en) * | 1997-11-14 | 2009-04-15 | Western Atlas International, Inc. | Seismic data acquisition and processing using non-linear distortion in a groundforce signal |
| CN101545981A (en) * | 2008-03-28 | 2009-09-30 | 中国石油天然气集团公司 | Method for realizing minimum phase of vibroseis seismic data zero-phase wavelet |
| CN102062872A (en) * | 2009-11-11 | 2011-05-18 | 中国石油天然气集团公司 | Multi-factor scanning method of controllable vibroseis |
| CN102692643A (en) * | 2011-03-23 | 2012-09-26 | 中国石油天然气集团公司 | Time varying controllable focal force signal deconvolution method |
| CN104035128A (en) * | 2013-03-06 | 2014-09-10 | 中国石油集团东方地球物理勘探有限责任公司 | Vibroseis pseudo random scanning signal generating method |
Non-Patent Citations (1)
| Title |
|---|
| "一种改善相关子波特性的扫描信号";张宏生 等;《物探装备》;20061231(第S1（2006）期);第33-41页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105572723A (en) | 2016-05-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105572723B (en) | The design method of controlled source scanning signal based on autocorrelation wavelet | |
| RU2579164C1 (en) | Handling method for determining quality of geologic environment | |
| AU2009282411B2 (en) | Removal of surface-wave noise in seismic data | |
| CN104950326B (en) | The method for designing of controlled source nonlinear scanning signal based on target zone frequency spectrum | |
| US7917295B2 (en) | Modeling and filtering coherent noise in seismic surveying | |
| CN101598812B (en) | Method for removing abnormal noise in single-point reception of seismic record by digital detector | |
| CN109669212B (en) | Seismic data processing method, stratum quality factor estimation method and device | |
| CN102062872A (en) | Multi-factor scanning method of controllable vibroseis | |
| CN110018517A (en) | A kind of multiple dimensioned ground micro-seismic inverse time interference localization method | |
| CN103675899A (en) | Method for expanding post-stack seismic data frequency band based on wavelet compression | |
| CN106199532A (en) | Gpr Signal noise-reduction method based on mixing Fourier wavelet analysis | |
| CN108845357A (en) | A method of the equivalent quality factor in stratum is estimated based on the synchronous wavelet transformation that squeezes | |
| CN106680874A (en) | Harmonic noise suppression method based on waveform morphology sparse modeling | |
| CN102721979B (en) | Seismic data-based thin layer automatic interpretation and thickness prediction method and device | |
| CN113887398A (en) | GPR signal denoising method based on variational modal decomposition and singular spectrum analysis | |
| CN107436450A (en) | A kind of seismic signal bandwidth broadning method based on continuous wavelet transform | |
| CN106019376A (en) | Seismic wave compensation method constructed through frequency driving space-variant Q value model | |
| CN100501450C (en) | Method for performing waveform stretching compensation for seismic signal | |
| CN102323615A (en) | Method for reservoir prediction and fluid identification with earthquake data and device | |
| CN101874744A (en) | Ultrasonic guided wave parameter measurement method used for long bone analysis | |
| RU2412454C2 (en) | Method to process seismic data using discrete wavelet transform | |
| Shao et al. | Seismic data antialiasing interpolation using sparse Radon transform and dynamic mask function | |
| CN100412569C (en) | Method for conducting earthquake signal high frequency compensation utilizing earthquake micro metering | |
| CN109975874B (en) | Controllable seismic source scanning signal design method based on damping rake wavelets | |
| CN106324702A (en) | Quantitative evaluation method of seismic interferometry imaging observation system design |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190802 Address after: 100728 Beijing, Chaoyangmen, North Street, No. 22, No. Co-patentee after: SINOPEC OILFIELD SERVICE CORPORATION Patentee after: China Petroleum Chemical Industry Group Co Ltd Co-patentee after: SINOPEC PETROLEUM ENGINEERING GEOPHYSICS CO., LTD. Co-patentee after: Victory point company of petroleum works geophysics Co., Ltd of China Petrochemical Industry Address before: 257086 Shandong Province, Dongying city Dongying District Niuzhuang town before the Street No. 70 Patentee before: Victory point company of petroleum works geophysics Co., Ltd of China Petrochemical Industry |