CN101545974A - Stratigraphic profile acoustic detection method for sub-bottom profiler - Google Patents

Stratigraphic profile acoustic detection method for sub-bottom profiler Download PDF

Info

Publication number
CN101545974A
CN101545974A CN200910097904A CN200910097904A CN101545974A CN 101545974 A CN101545974 A CN 101545974A CN 200910097904 A CN200910097904 A CN 200910097904A CN 200910097904 A CN200910097904 A CN 200910097904A CN 101545974 A CN101545974 A CN 101545974A
Authority
CN
China
Prior art keywords
sub
bottom profiler
signal
frequency response
frequency modulation
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.)
Pending
Application number
CN200910097904A
Other languages
Chinese (zh)
Inventor
施国全
王福林
韦俊霞
曹海林
郑军
范进良
赵爱君
刘畅
李小英
刘强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
HANGZHOU RESOUND MARINE INSTRUMENTS CO Ltd
Original Assignee
HANGZHOU RESOUND MARINE INSTRUMENTS CO Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by HANGZHOU RESOUND MARINE INSTRUMENTS CO Ltd filed Critical HANGZHOU RESOUND MARINE INSTRUMENTS CO Ltd
Priority to CN200910097904A priority Critical patent/CN101545974A/en
Publication of CN101545974A publication Critical patent/CN101545974A/en
Pending legal-status Critical Current

Links

Images

Abstract

The invention provides a stratigraphic profile acoustic detection method for a sub-bottom profiler, which comprises the following steps that: 1, the sub-bottom profiler firstly generates linear frequency modulation emission signals, and the linear frequency modulation emission signals are added with a sine-squared one eighth window function to form Chirp emission signals; and 2, the Chirp emission signals are transmitted to water body through a transmitting transducer, a receiving hydrophone receives echoed signals returned by stratum and then performs related processing on the echoed signals, and result output of the related processing is displayed on a display screen to give a string of pixels with different shading for reflecting interface characteristics of the stratum, wherein when the echoed signals are subjected to the related processing, frequency response compensation is performed on the sub-bottom profiler so as to detect a frequency response curve of the sub-bottom profiler, and the detected frequency response curve is utilized to perform frequency domain compensation on the echoed signals so as to remove frequency response interference of the sub-bottom profiler and improve signal correlation.

Description

The stratigraphic section acoustic detection method that is used for sub-bottom profiler
Technical field
The present invention relates to a kind of stratigraphic section echoed signal that is used for oceanographic instrumentations such as sub-bottom profiler and handle particularly a kind of stratigraphic section acoustic detection method that is used for sub-bottom profiler
Background technology
In the matched filter processing of sub-bottom profiler, for suppressed sidelobes, carry out windowing process to transmitting, conventional sub-bottom profiler adopts Blackman-Harris window or Dolph-Chebyshev window to carry out windowing process to transmitting usually.
Said method has following shortcoming:
1, the side lobe levels of auto-correlation output is higher, and stratigraphic resolution is reduced;
2, cause the bandwidth that transmits to narrow down after the windowing;
3, forward and backward in windowing, the transmit signal energy loss is bigger, exports at amplifier under the situation of same power, and the stratum penetration depth decreases.
Summary of the invention
Fundamental purpose of the present invention is to provide a kind of stratigraphic section acoustic detection method that is used for sub-bottom profiler, and the transmit signal energy loss that exists in the solution prior art is bigger, and fuzzy relatively shortcoming and defect is differentiated on the stratum.
The present invention is achieved through the following technical solutions.A kind of stratigraphic section acoustic detection method that is used for sub-bottom profiler may further comprise the steps:
One, sub-bottom profiler at first produces linear frequency modulation and transmits
x(t)=Re{exp[j2πt(f1+Kt/2)]} (1)
In the formula: x (t) transmits for linear frequency modulation; Re Re{exp[j2 π t (f1+Kt/2)] } for it is got real part;
f 1Initial frequency for FM signal; Then chirp rate K is: K=Δ f/ Δ t
△ f is the frequency span of FM signal;
△ t is the time span of FM signal;
Linear Frequency Modulation transmits and adds sinusoidal square 1/8 windowed function w (t) among the x (t), forms Chirp and transmits:
s(t)=x(t)*w(t); (2)
Sinusoidal square 1/8 a windowed function w (t) is:
w ( t ) = 1 N / 8 - 1 < t &le; 7 N / 8 - 1 sin 2 ( 4 t&pi; / N ) 0 &le; t &le; N / 8 - 1 sin 2 ( &pi; / 2 - 4 ( t - 7 N / 8 + 1 ) &pi; / N 7 N / 8 - 1 < t &le; N - 1 - - - ( 3 )
In the formula, N is the sampling number that linear frequency modulation transmits;
Two, Chirp s emission signal s (t) is launched to water body through transmitting transducer, receiving hydrophone receives the echoed signal that the stratum is returned, then echoed signal is done relevant treatment, result's output of relevant treatment is shown on the display screen, provides the pixel of a string deep or light reflection bed boundary feature that differs.
Wherein, when echoed signal is carried out relevant treatment, sub-bottom profiler is carried out the frequency response compensation, record the frequency response curve of sub-bottom profiler, the frequency response curve that utilization records carries out frequency domain compensation to echoed signal, remove the frequency response of sub-bottom profiler and disturb, improve the correlativity of signal.
The stratigraphic section acoustic detection method that is used for sub-bottom profiler of the present invention, adopt sinusoidal square 1/8 windowing process technology, when reducing the transmit signal energy loss, the relevant output of stratum echoed signal has less side lobe levels, it is darker to make sub-bottom profiler and section plotter in the past compare earth penetrating, and the stratigraphic section layering of measurement is more clear.
Description of drawings
Fig. 1 is the signal graph of the Blackman-Harris window of traditional sub-bottom profiler employing;
Fig. 2 is the auto-correlation output map of the Blackman-Harris window of traditional sub-bottom profiler employing;
Fig. 3 adopts the signal graph of sinusoidal square 1/8 windowed function for the present invention;
Fig. 4 adopts the auto-correlation output map of sinusoidal square 1/8 windowed function for the present invention;
The echoed signal figure that Fig. 5 returns for stratum in the specific embodiment of the invention;
Fig. 6 is that the echoed signal relevant treatment is exported figure as a result in the specific embodiment of the invention.
Embodiment
The present invention is described in detail below in conjunction with accompanying drawing.
If the frequency range that linear frequency modulation transmits is 3kHz~8kHz, the initial frequency f1=3kHz of FM signal, the frequency span △ f=5kHz of FM signal, the time span △ t=10ms of FM signal, it is as follows that signals sampling frequency f s=100kHz is rewritten into the discrete digital form with (1) formula:
x(n)=Re{exp[j2πt(f1+Knt/2)]} (4)
In the formula: n=0~N-1.
Can produce the digitizing linear frequency modulation time series signal of 3kHz~8kHz by formula (4).
It is as follows that sinusoidal square 1/8 windowed function is rewritten into the discrete digital form:
w ( n ) = 1 N / 8 - 1 < n &le; 7 N / 8 - 1 sin 2 ( 4 n&pi; / N ) 0 &le; n &le; N / 8 - 1 sin 2 ( &pi; / 2 - 4 ( n - 7 N / 8 + 1 ) &pi; / N 7 N / 8 - 1 < n &le; N - 1 - - - ( 5 )
In the formula, w (n) is the time series of windowed function, n=0~N-1.
The discrete digital time series of sinusoidal square 1/8 windowing linear FM signal
s(n)=x(n)*w(n) (6)
In the formula, s (n) is the discrete digital time series of windowed function linear FM signal, n=0~N-1.
For example, when n=0~9, the linear frequency modulation s emission signal s (n) after the windowing is respectively 0.0002,0.0005,0.0008,0.0008 ,-0.0001 ,-0.0019 ,-0.0047 ,-0.0084 ,-0.0123 ,-0.0157.
Linear frequency modulation after the windowing transmits and launches to water body through transmitting transducer, after the emission, the echoed signal that the stratum is returned as shown in Figure 5, the result of echoed signal relevant treatment is as shown in Figure 6.
At last, it is also to be noted that what more than enumerate only is specific embodiments of the invention.Obviously, the invention is not restricted to above embodiment, many distortion can also be arranged.All distortion that those of ordinary skill in the art can directly derive or associate from content disclosed by the invention all should be thought protection scope of the present invention.

Claims (1)

1, a kind of stratigraphic section acoustic detection method that is used for sub-bottom profiler is characterized in that, may further comprise the steps:
One, sub-bottom profiler at first produces the linear frequency modulation x (t) that transmits:
x(t)=Re{exp[j2πt(f1+Kt/2)]} (1)
In the formula:
X (t) transmits for linear frequency modulation; Re{exp[j2 π t (f1+Kt/2)] } for it is got real part;
F1 is the initial frequency of FM signal; Then chirp rate K is: K=Δ f/ Δ t
Δ f is the frequency span of FM signal;
Δ t is the time span of FM signal;
Linear Frequency Modulation transmits and adds sinusoidal square 1/8 windowed function w (t) among the x (t), forms linear frequency modulation s emission signal s (t):
s(t)=x(t)*w(t); (2)
In the formula, w (t) is the time series of windowed function, t=0~N-1;
Sinusoidal square 1/8 a windowed function w (t) is:
w ( t ) = 1 N / 8 - 1 < t &le; 7 N / 8 - 1 sin 2 ( 4 t&pi; / N ) 0 &le; t &le; N / 8 - 1 sin 2 ( &pi; / 2 - 4 ( t - 7 N / 8 + 1 ) &pi; / N 7 N / 8 - 1 < t &le; N - 1 - - - ( 3 )
In the formula, N is the sampling number that linear frequency modulation transmits;
Two, linear frequency modulation s emission signal s (t) is launched to water body through transmitting transducer, receiving hydrophone receives the echoed signal that the stratum is returned, then echoed signal is done relevant treatment, result's output of relevant treatment is shown on the display screen, provides the pixel of a string deep or light reflection bed boundary feature that differs;
Wherein, when echoed signal is carried out relevant treatment, sub-bottom profiler is carried out the frequency response compensation, record the frequency response curve of sub-bottom profiler, the frequency response curve that utilization records carries out frequency domain compensation to echoed signal, remove the frequency response of sub-bottom profiler and disturb, improve the correlativity of signal.
CN200910097904A 2009-04-23 2009-04-23 Stratigraphic profile acoustic detection method for sub-bottom profiler Pending CN101545974A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200910097904A CN101545974A (en) 2009-04-23 2009-04-23 Stratigraphic profile acoustic detection method for sub-bottom profiler

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200910097904A CN101545974A (en) 2009-04-23 2009-04-23 Stratigraphic profile acoustic detection method for sub-bottom profiler

Publications (1)

Publication Number Publication Date
CN101545974A true CN101545974A (en) 2009-09-30

Family

ID=41193212

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200910097904A Pending CN101545974A (en) 2009-04-23 2009-04-23 Stratigraphic profile acoustic detection method for sub-bottom profiler

Country Status (1)

Country Link
CN (1) CN101545974A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103884625A (en) * 2014-04-15 2014-06-25 黄河水利委员会黄河水利科学研究院 Data-fusion-based reservoir bottom mud density detection method and system
CN106597546A (en) * 2016-12-08 2017-04-26 黄河水利委员会黄河水利科学研究院 Method for detecting river bottom surface sediment physical properties
CN106814360A (en) * 2015-11-30 2017-06-09 江苏中海达海洋信息技术有限公司 A kind of multibeam sounding system based on linear FM signal
CN113589294A (en) * 2021-07-20 2021-11-02 武汉船舶通信研究所(中国船舶重工集团公司第七二二研究所) Method, system, equipment and medium for processing underwater shallow stratum profile data

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103884625A (en) * 2014-04-15 2014-06-25 黄河水利委员会黄河水利科学研究院 Data-fusion-based reservoir bottom mud density detection method and system
CN106814360A (en) * 2015-11-30 2017-06-09 江苏中海达海洋信息技术有限公司 A kind of multibeam sounding system based on linear FM signal
CN106814360B (en) * 2015-11-30 2019-07-09 江苏中海达海洋信息技术有限公司 A kind of multibeam sounding system based on linear FM signal
CN106597546A (en) * 2016-12-08 2017-04-26 黄河水利委员会黄河水利科学研究院 Method for detecting river bottom surface sediment physical properties
CN106597546B (en) * 2016-12-08 2019-04-02 黄河水利委员会黄河水利科学研究院 A method of detection river bed surface layer silt physical property
CN113589294A (en) * 2021-07-20 2021-11-02 武汉船舶通信研究所(中国船舶重工集团公司第七二二研究所) Method, system, equipment and medium for processing underwater shallow stratum profile data

Similar Documents

Publication Publication Date Title
Park et al. Multimodal analysis of high frequency surface waves
US9268022B2 (en) Underwater detection device and underwater detecting method
EP2624015B1 (en) Sweep design for seismic sources
CN101627322A (en) Vibroseis acquisition method
CA2778760A1 (en) Methods and apparatus to process time series data for propagating signals in a subterranean formation
CN103698807B (en) Scalariform two-dimensional wide-band observation system design method
CN102279394B (en) Low-frequency narrow-band high-resolution ultrasonic detection imaging method
CN107436451B (en) A kind of amplitude spectral method of automatic calculating seismic data optical cable coupled noise degree of strength
CN101545974A (en) Stratigraphic profile acoustic detection method for sub-bottom profiler
CN106526678B (en) A kind of wave field separation method and device of reflected acoustic wave well logging
Brooks et al. Green's function approximation from cross-correlations of 20–100hz noise during a tropical storm
JP5710289B2 (en) Underground radar equipment
CN1981212B (en) Method for attenuating noise in seismic data using complex trace diversity filter
CN100485413C (en) Frequency-scanning signal time-frequency correlation detection and time delay estimating method
Zhang et al. Lithospheric imaging through reverberant layers: sediments, oceans, and glaciers
Yamada et al. Limitation of the predominant-period estimator for earthquake early warning and the initial rupture of earthquakes
Zou et al. Enhancing low-frequency water-column acoustic reflections in marine multichannel seismic data for seismic oceanography
CN104502968B (en) Controlled source geological data detection method based on threshold value multistage median filtering
RU2517983C1 (en) Method of profiling bottom deposits
Zheng et al. The application of a homemade boomer source in offshore seismic survey: From field data acquisition to post-processing
Woolfe et al. Optimized extraction of coherent arrivals from ambient noise correlations in a rapidly fluctuating medium
KR101839236B1 (en) Subbottom profiler and subbottom profiling method
CN104345347A (en) Logging curve recovery method used for compact gas-containing sandstone reservoir prediction
CN117872478A (en) Weak reflection signal extraction method and device for array acoustic logging data
CN113672840B (en) Stoneley wave attenuation extraction method and stoneley wave attenuation extraction device

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20090930