CN108008465B - Signal source triggered time acquisition methods in engineering geophysics signal acquisition - Google Patents
Signal source triggered time acquisition methods in engineering geophysics signal acquisition Download PDFInfo
- Publication number
- CN108008465B CN108008465B CN201711015650.3A CN201711015650A CN108008465B CN 108008465 B CN108008465 B CN 108008465B CN 201711015650 A CN201711015650 A CN 201711015650A CN 108008465 B CN108008465 B CN 108008465B
- Authority
- CN
- China
- Prior art keywords
- time
- data
- signal
- engineering geophysics
- time point
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
- G01V11/002—Details, e.g. power supply systems for logging instruments, transmitting or recording data, specially adapted for well logging, also if the prospecting method is irrelevant
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention relates to signal source triggered time acquisition methods in engineering geophysics signal acquisition, including 1: setting physical field sensor-triggered threshold value；Physical prospecting instrument sampling time length L, pre-sampling time span M are set；2: engineering geophysics instrument acquires data, and data are stored in caching；3: physics field sensor starts synchronous acquisition signal, and corresponding physical quantity reading is compared by the moment with activation threshold value；4: recording temporal information t of the collected physics field intensity of physics field sensor greater than the activation threshold value when0；5: engineering geophysics instrument collects t0It when+(L-M) time point, stops working, retains t in caching0- M to t0The data at+(L-M) moment；6: Hilbert transform being carried out to collected data, finds out the mutation time point of data instantaneous frequency；Step 7: mutation time point in collected signal and mutation time point subsequent N number of time point being carried out curve fitting, the physical field practical triggered time is found out.
Description
Technical field
The present invention relates to Engineering geophysical exploration technical fields, and in particular to signal in a kind of engineering geophysics signal acquisition
Source triggered time acquisition methods.
Background technique
Engineering geophysical exploration, abbreviation engineering geophysics, its detection target are underground rock-soil layer or building structure, right
After detecting target load nature or artificial physics field, by its variation of Instrument observation, the spatial dimension to determine concealed target is (big
Small, shape, distribution), and the physical parameter of objective body is measured, solve a kind of physical prospecting method of geological problem.Make at this stage
In Engineering geophysical exploration instrument, mainly detected by the way of artificial field.
Then, in most of engineering geophysics operation, the excitation that prospecting instrument records artificial excitation's physical field is required
Time.There are mainly two types of methods at this stage: acquiring after first method triggering, i.e., one biography of setting near physics field excitation end
Sensor simultaneously sets a threshold value, after device ready, the neighbouring sensors measure physical field intensity in excitation end, when measured value is more than threshold value
When, prospecting instrument records the time, and in this, as the time of physics field excitation, synchronous signal receiving end acquires number from this moment on
According to；Second method is to acquire in advance, and sensor and given threshold, device ready are arranged equally near physics field excitation end
Afterwards, instrument receiving end first starts the cycle over acquisition data, while exciting end sensors measure physical field intensity nearby, when measured value is super
When crossing threshold value, prospecting instrument records the time, will prolong a set time before the record time, as acquisition starting point.
But both methods is problematic in that, for first method, since external interference and instrument bottom are made an uproar not
May be unlimited low, thus threshold value cannot be set to it is too low, that is to say, that the physical field firing time of instrument record actually compares object
The reason practical firing time in field wants late.For second method, complete physical field excitation process has been recorded really, still
Sensor sample rate is limited, can not accurately know the practical firing time of physical field, can only be with physics field intensity when handling signal
Maximum value is handled as time point.In this way handle the shortcomings that it is also obvious that signal pass through detection target after, the vibration of signal
The physical properties such as width, frequency, phase, velocity of wave can change, and be carried out with the time point of physical field maximum of intensity in signal
Processing, inevitably has error.
Summary of the invention
It is an object of the invention to, propose to believe in a kind of engineering geophysics signal acquisition aiming at the problems existing in the prior art
Number source triggered time acquisition methods, this method can accurately find the firing time of artificial physics field.
In order to solve the above technical problems, the signal source triggered time is obtained in a kind of engineering geophysics signal acquisition disclosed by the invention
Take method, which is characterized in that it includes the following steps:
Step 1: arranging corresponding physics field sensor at physics field excitation end, and set activation threshold value, when physical field passes
The correspondence physics field intensity of sensor is greater than the time corresponding when this activation threshold value for physical field firing time reference value；According to reality
It is L and pre-sampling time span M that border, which needs to be arranged engineering geophysics instrument sampling time length,；
Step 2: engineering geophysics instrument starts to acquire data, and data are stored in engineering geophysics instrument in a manner of recycling storage
Caching in；
Step 3: physics field sensor described in step 1 starts synchronous acquisition signal, and the moment will corresponding physical quantity reading
It is compared with the activation threshold value；
Step 4: when recording the collected physics field intensity of the physics field sensor greater than the activation threshold value when
Between information t0；
Step 5: engineering geophysics instrument collects t0It when+(L-M) time point, stops working, it is slow to retain engineering geophysics instrument
Deposit middle t0- M to t0The data at+(L-M) moment；
Step 6: Hilbert transform being carried out to collected data, the instantaneous frequency of data is calculated, finds out data wink
When frequency mutation time point t1；
Step 7: by mutation time point t in collected signal1And mutation time point t1Subsequent N number of time point carries out
Curve matching finds out the physical field practical triggered time by fit equation.
Beneficial effects of the present invention:
The present invention obtains physics field excitation time proximity point by the instantaneous frequency of signal, then approximate to physics field excitation
Sampled value near time point carries out curve fitting, calculated physical field firing time, very close with practical firing time,
The accuracy of physical field firing time can preferably be improved.
Detailed description of the invention
Fig. 1 is flow chart of the method for the present invention；
Fig. 2 is signal acquisition schematic diagram of the invention；
Fig. 3 is collected signal schematic representation in the present invention；
Fig. 4 is that instantaneous frequency abnormal point extracts schematic diagram in the present invention.
Specific embodiment
Below in conjunction with the drawings and specific embodiments, the present invention is described in further detail:
Signal source triggered time acquisition methods are as shown in Figure 1, it includes in a kind of engineering geophysics signal acquisition of the invention
Following steps:
Step 1: arranging that corresponding physics field sensor (can detect the physical field of transmitting at physics field excitation end
Sensor, such as wave detector, be just placed on transmitting end edge on), and set activation threshold value (i.e. critical value, default with sensor reading
Time greater than this value is physical field firing time reference value, but this time and the practical firing time of physical field have error, generally
Carry out given threshold, rule of thumb with the positive correlation of itself noise level of engineering geophysics instrument)；Engineering geophysics is set according to actual needs
Instrument sampling length is L and pre-sampling length M；The pre-sampling length M is greater than temporal information t0With the practical triggering of physical field
The difference of time；The pre-sampling time span M is 3%~5% that sampling length is L generally with empirical estimating；
The conventional method of geophysical exploration be by natural or artificial physical field (such as electric field, elastic wave field etc., here
Using artificial physics field) it is added on objective body to be explored, the variation of physical field is then measured, is derived by measurement result
Physical property to exploration targets body.Here excitation end refers to the device of transmitting physical field, such as the antenna of radar, elastic wave
Focus etc.；
Step 2: engineering geophysics instrument starts to acquire data, and data are stored in engineering geophysics instrument in a manner of recycling storage
Caching in (such as sensor cache size be 10 data units, after starting acquisition, data 1, data 2, data 3, data 4
Sequence deposit caching, caching is full when to the 10th data, and data 11 are stored in 1 number unit, overrides data 1, data 12 are covered
Lid data 2, and so on)；
Step 3: physics field sensor described in step 1 starts synchronous acquisition signal, and the moment will corresponding physical quantity reading
It is compared with the activation threshold value；
Step 4: when recording the collected physics field intensity of the physics field sensor greater than the activation threshold value when
Between information t0；
Step 5: engineering geophysics instrument collects t0It when+(L-M) time point, stops working, it is slow to retain engineering geophysics instrument
Deposit middle t0- M to t0The data at+(L-M) moment；As shown in Figure 2；
Step 6: Hilbert transform being carried out to collected data, the instantaneous frequency of data is calculated, finds out data wink
When frequency mutation time point t1；
Step 7: by mutation time point t in collected signal1With mutation time point t1Subsequent N number of time point march
Line fitting, finds out the physical field practical triggered time by fit equation.
In the step 2 of above-mentioned technical proposal, after data cached be filled with, new collected data are always replaced in caching
The time of storage earliest data.
N value is 5 in the step 7 of above-mentioned technical proposal.
Method of the invention is verified below,
Collected signal is not difficult to find out that physical field is practical as shown in figure 3, here by taking such as minor function as an example from function
Triggered time be 50.6.
In formula, tx is the preset physical field practical triggered time, and it is the time that value, which is 50.6, t, here, and e is natural constant,
A (t) is engineering geophysics instrument to collected data；
Then according to the present invention step 6 and 7 mode to A (t) carry out Hilbert transform, calculate the instantaneous frequency of signal
Rate finds out the catastrophe point t of signal transient frequency1, as shown in figure 4, calculated instantaneous frequency catastrophe point t here1It is 51.
By t in collected signal1, i.e., the point near 51 carries out curve fitting, and finds out physical field reality by fit equation
In the border triggered time, the calculated physical field triggered time is between 50.6001, with preset physical field practical triggered time here
It is very close.Hence it is demonstrated that the feasibility of the present invention program.
The content that this specification is not described in detail belongs to the prior art well known to professional and technical personnel in the field.
Claims (5)
1. signal source triggered time acquisition methods in a kind of engineering geophysics signal acquisition, which is characterized in that it includes the following steps:
Step 1: arranging corresponding physics field sensor at physics field excitation end, and set activation threshold value, when physics field sensor
Correspondence physics field intensity when being greater than this activation threshold value corresponding time be physical field firing time reference value；According to practical need
It is L and pre-sampling time span M that engineering geophysics instrument sampling time length, which is arranged,；
Step 2: engineering geophysics instrument starts to acquire data, and data are stored in the slow of engineering geophysics instrument in a manner of recycling storage
In depositing；
Step 3: physics field sensor described in step 1 starts synchronous acquisition signal, and the moment will corresponding physical quantity reading and institute
Activation threshold value is stated to be compared；
Step 4: recording time letter of the collected physics field intensity of the physics field sensor greater than the activation threshold value when
Cease t0；
Step 5: engineering geophysics instrument collects t0It when+(L-M) time point, stops working, retains in engineering geophysics instrument caching
t0- M to t0The data at+(L-M) moment；
Step 6: Hilbert transform being carried out to collected data, the instantaneous frequency of data is calculated, finds out the instantaneous frequency of data
The mutation time point t of rate1；
Step 7: by mutation time point t in collected signal1And mutation time point t1Subsequent N number of time point carries out curve
Fitting, finds out the physical field practical triggered time by fit equation.
2. signal source triggered time acquisition methods in engineering geophysics signal acquisition according to claim 1, it is characterised in that:
In the step 2, after data cached be filled with, new collected data replace the time stored in caching earliest data.
3. signal source triggered time acquisition methods in engineering geophysics signal acquisition according to claim 1, it is characterised in that:
N value is 5 in the step 7.
4. signal source triggered time acquisition methods in engineering geophysics signal acquisition according to claim 1, it is characterised in that:
The pre-sampling time span M is greater than temporal information t0With the difference in physical field practical triggered time.
5. signal source triggered time acquisition methods in engineering geophysics signal acquisition according to claim 1, it is characterised in that:
The pre-sampling time span M is 3%~5% that sampling time length is L.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711015650.3A CN108008465B (en) | 2017-10-25 | 2017-10-25 | Signal source triggered time acquisition methods in engineering geophysics signal acquisition |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711015650.3A CN108008465B (en) | 2017-10-25 | 2017-10-25 | Signal source triggered time acquisition methods in engineering geophysics signal acquisition |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108008465A CN108008465A (en) | 2018-05-08 |
| CN108008465B true CN108008465B (en) | 2019-06-14 |
Family
ID=62051946
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711015650.3A Active CN108008465B (en) | 2017-10-25 | 2017-10-25 | Signal source triggered time acquisition methods in engineering geophysics signal acquisition |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108008465B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN201681269U (en) * | 2010-07-27 | 2010-12-22 | 广州海洋地质调查局 | Multi-path submarine seismograph precision timer |
| CN103399508A (en) * | 2013-07-29 | 2013-11-20 | 电子科技大学 | Synchronous triggering management device and synchronous triggering management method |
| CN105300508A (en) * | 2015-10-30 | 2016-02-03 | 北京控制工程研究所 | Double-threshold detection circuit of ultrasonic wave signal propagation time, and detection method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10101495B2 (en) * | 2016-03-10 | 2018-10-16 | Exxonmobil Research And Engineering Company | Magneto-seismic exploration method and system |
-
2017
- 2017-10-25 CN CN201711015650.3A patent/CN108008465B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN201681269U (en) * | 2010-07-27 | 2010-12-22 | 广州海洋地质调查局 | Multi-path submarine seismograph precision timer |
| CN103399508A (en) * | 2013-07-29 | 2013-11-20 | 电子科技大学 | Synchronous triggering management device and synchronous triggering management method |
| CN105300508A (en) * | 2015-10-30 | 2016-02-03 | 北京控制工程研究所 | Double-threshold detection circuit of ultrasonic wave signal propagation time, and detection method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108008465A (en) | 2018-05-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Diehl et al. | Automatic S-wave picker for local earthquake tomography | |
| FI125828B (en) | Precipitation sensor | |
| AU2011100086A4 (en) | A method for measuring biomass of tree-roots by ground-penetrating radar | |
| CN107576674A (en) | A kind of method based on GPR measurement soil compression degree | |
| CN103196465B (en) | Method for separating noise of sensing signal and extracting signal of phase-sensitive optical time domain reflectometer | |
| Küperkoch et al. | Automated determination of S-phase arrival times using autoregressive prediction: application to local and regional distances | |
| WO2016112147A1 (en) | Gauge length optimization in distributed vibration sensing | |
| CN104483700B (en) | Formation fracture is monitored and early warning system and method | |
| NO20130782A1 (en) | IDENTIFICATION OF INVALID SEISMIC DATA | |
| CN104880728B (en) | A kind of method and relevant apparatus of the quality control for broadband scanning | |
| Miller et al. | The response of the Dines anemometer to gusts and comparisons with cup anemometers | |
| CN202420427U (en) | Signal detection system for electromagnetic thickness gauge | |
| CN106594529A (en) | Leakage signal propagation speed measurement method based on short-time energy and linear fitting | |
| CN109425894A (en) | A kind of seismic anomaly road detection method and device | |
| CA2818641C (en) | Methods of hydrocarbon detection using spectra dominant frequency and measures of energy decay on the low side and high side of spectra dominant frequency | |
| CN105866234A (en) | Eddy current detection-Barkhausen noise detection combined nondestructive testing apparatus and method for ferromagnetic material | |
| CN108008465B (en) | Signal source triggered time acquisition methods in engineering geophysics signal acquisition | |
| JP4065946B2 (en) | Initial motion detection method for multi-component AE waveform | |
| Abraham et al. | Use of a handheld broadband EM induction system for deriving resistivity depth images | |
| Blumentritt et al. | A comparison of magnetic susceptibility measurement techniques and ferrimagnetic component analysis from recent sediments in Lake Pepin (USA) | |
| Bohleber | Ground-penetrating radar assisted ice core research: the challenge of Alpine glaciers and dielectric ice properties | |
| CN106610509A (en) | Time domain processing method for transient electromagnetic data | |
| Vatandoost et al. | Automated multi-sensor petrophysical core logging | |
| RU2677837C1 (en) | Method of operational forecasting of earthquake hypocenter coordinates | |
| CN109061748A (en) | The method that Mine transient electromagnetic secondary electric potential relative error determines geologic body information |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |