CN102435866A - Method for quickly identifying interference of ground object during archaeological detection of ground penetrating radar - Google Patents
Method for quickly identifying interference of ground object during archaeological detection of ground penetrating radar Download PDFInfo
- Publication number
- CN102435866A CN102435866A CN2011103701048A CN201110370104A CN102435866A CN 102435866 A CN102435866 A CN 102435866A CN 2011103701048 A CN2011103701048 A CN 2011103701048A CN 201110370104 A CN201110370104 A CN 201110370104A CN 102435866 A CN102435866 A CN 102435866A
- Authority
- CN
- China
- Prior art keywords
- source
- interference
- ground
- line
- interference 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.)
- Pending
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a method for quickly identifying the interference of a ground object during archaeological detection of a ground penetrating radar. The problem of electromagnetic interference caused by the ground object cannot be well solved by the conventional archaeological detection operation of the ground penetrating radar. The method comprises the following steps of: firstly, establishing a space coordinate system for a detection area, and establishing a rectangular coordinate system by using a detection line where a transmitting antenna is positioned as an x axis, a starting point of the detection line as an original point, a ground line which is vertical to the detection line as a y axis and a downward line which is vertical to the ground as a z axis; secondly, acquiring position parameters of an interference source, and calculating an interference time distance curve; and finally, charting data which is obtained by the archaeological detection of a field ground penetrating radar, comparing a chart of the data with the chart of the interference time distance curve which is obtained by calculation, and identifying ground interference from an actually-detected section. The method is applicable to calculation and identification of the interference of the ground object, and has the advantages of quickness, simplicity, convenience, intuitive result and the like.
Description
Technical field
The inventive method belongs to GPR archaeology exploration engineering field, relates in particular to the Electromagnetic Interference that the above object in a kind of ability quick identification ground produces.
Background technology
GPR becomes a kind of useful GEOPHYSICS FOR SHALLOW LAYER detection method in the science and technology archaeology with its high resolving power, advantage such as efficient, harmless.Because regular meeting adopted non-screened shielded antanna when the GPR archaeology was surveyed, the electromagnetic wave portion of energy of this kind antenna emission is aloft propagated.Electromagnetic wave attenuation is very little in the air.The Electromagnetic Interference of generations such as above object in ground such as iron block, electric wire, entanglement, buildings can be superimposed upon on the underground useful signal, and compares with useful signal, and interfering energy is strong, even might cover useful signal.Along with industrialization and quickening of urbanization process, the Electromagnetic Interference that the above object in ground produces is one of subject matter of facing of GPR archaeology investigation and prospecting.So need the technology of this interference of the easy quick identification of a kind of ability.
Summary of the invention
The purpose of this invention is to provide a kind of GPR and engage in archaeological studies when surveying, the method for the Electromagnetic Interference that the above interference source in ability quick identification ground produces.
Set up space coordinates to surveying the district; When coordinate system is set up with emitting antenna place survey line for
axle; The survey line starting point is an initial point; The vertical survey line in ground is
axle, vertical ground downwards for
axle set up rectangular coordinate system.
The practical implementation step of the inventive method:
Step (1). obtain the position of interference source parameter;
1-1. when interference source is the some Source Type; Point source interference source Q point coordinate is (
); Obtain point source interference source Q subpoint on survey line, and measure the position a of subpoint on survey line; Subpoint is to the distance b of survey line on the ground to measure point source interference source Q, and measurement point source interference source Q is to the vertical range c on ground;
1-2 When the interference source for the line source type, line source interference source
used at any point on a line source
(
) and with a non-zero vector
parallel straight expression; online source of interference sources
office to take some distance from two o'clock
and
, then Refer to Step 1-1 point source interference source Q parameter acquisition obtained coordinates of two points
(
) and
(
), a non-zero vector
, the
1-3. when interference source is the face Source Type; The equation of face source interference source S is
, and
is not 0 simultaneously; On the interference source S of face source, choose three the not point of conllinear
,
and
arbitrarily; Obtain coordinate
(
),
(
) and
(
) then with reference to step 1-1 point source interference source Q parameter acquiring mode at 3, and obtain 4 parameters of face source interference source S through following system of equations:
(2)
According to the relevant mathematical theorem of plane equation; Can establish
in the above-mentioned system of equations,
for arbitrarily not being 0 constant.
Step (2). calculate and disturb T-X curve;
The GPR archaeology is surveyed and is adopted profile method to measure often, and the relative position of emitting antenna and receiving antenna remains unchanged, along the same moved further of survey line; The antenna distance of emitting antenna and receiving antenna is
; And the line of emitting antenna and receiving antenna is
with
axle clamp angle;
; If the emitting antenna coordinate is (
), the receiving antenna coordinate is (
).
2-1. if interference source is a Source Type,
(3)
Wherein
is the skyborne velocity of propagation of electromagnetic wave,
the electromagnetic first break time of interference that produces for the point source interference source;
2-2. if interference source is the line source type,
Wherein: the electromagnetic first break time of interference that
produces for the line source interference source;
2-3. if interference source is the face Source Type,
When the both sides of emitting antenna and receiving antenna interference source or have one during in the face source in the face source;
When the homonymy of emitting antenna and receiving antenna interference source in the face source,
Wherein
is the electromagnetic first break time of interference that face source interference source produces;
Step (3). disturb on identification ground;
Open-air GPR archaeology detection gained data are carried out to figure disturb T-X curve to become figure to compare, and identify interference from measured section with the calculating gained.
The inventive method beneficial effect is following:
To calculating and the identification that ground object disturbs, have advantages such as quick, easy, visual result.Basically can when obtaining calculating parameter,, be easy to disturb identification at archaeology data acquisition scene with regard to the corresponding interference curve that calculates that can be fast and convenient.
Above-mentioned interferometer formula (3)~(5), in the time of can surveying for GPR archaeology, the compacting of the interference that the ground interference source produces provides the basis of usefulness.
Description of drawings
When Fig. 1 is image data, obtain the synoptic diagram of point source interference source parameter;
Fig. 2 is the situation of interference source in the instance 1;
Fig. 3 (a) is the section of archaeology GPR survey line gained among Fig. 2;
Fig. 3 (b) is according to formula (3) result calculated;
Fig. 4 is the situation of interference source in the instance 2;
Fig. 5 (a) is the section of archaeology GPR survey line gained among Fig. 4;
Fig. 5 (b) is according to formula (4) and formula (5) result calculated.
Embodiment
Below in conjunction with accompanying drawing the inventive method step is described further.
Set up space coordinates to surveying the district; When coordinate system is set up with emitting antenna place survey line for
axle; The survey line starting point is an initial point; The vertical survey line in ground is
axle, vertical ground downwards for
axle set up rectangular coordinate system.
Step (1). obtain the position of interference source parameter.
1-1. when interference source is the some Source Type; Can as shown in Figure 1ly obtain point source interference source Q location parameter; Point source interference source Q point coordinate is (
); Obtain point source interference source Q subpoint on survey line, and measure the position a of subpoint on survey line; Subpoint is to the distance b of survey line on the ground to measure point source interference source Q, and measurement point source interference source Q is to the vertical range c on ground;
1-2.? When the interference source is a line source type, you can get a line source as to the source of interference P positional parameters, line source interference source
used at any point on a line source
(
) and with a non-zero vector
parallel straight expression; Online source of interference sources
office to take some distance from two points
and
, and then refer to step 1-1 point source interference source Q parameter to obtain the coordinates of points obtained
(
) and
(
), a non-zero vector
, the
1-3. when interference source is the face Source Type; Can obtain face source interference source S location parameter as getting off; The equation of face source interference source S is
, and
is not 0 simultaneously; On the interference source S of face source, choose three the not point of conllinear
,
and
arbitrarily; Obtain coordinate
(
),
(
) and
(
) then with reference to step 1-1 point source interference source Q parameter acquiring mode at 3, and obtain 4 parameters of face source interference source S through following system of equations:
According to the relevant mathematical theorem of plane equation; Can establish
in the above-mentioned system of equations,
for arbitrarily not being 0 constant.
Step (2). calculate and disturb T-X curve.
The GPR archaeology is surveyed and is adopted profile method to measure often, and the relative position of emitting antenna and receiving antenna remains unchanged, along the same moved further of survey line; The antenna distance of emitting antenna and receiving antenna is
; And the line of emitting antenna and receiving antenna is
with
axle clamp angle;
; If the emitting antenna coordinate is (
), the receiving antenna coordinate is (
).
2-1. if interference source is a Source Type,
Wherein
is the skyborne velocity of propagation of electromagnetic wave,
the electromagnetic first break time of interference that produces for the point source interference source;
2-2. if interference source is the line source type,
Wherein: the electromagnetic first break time of interference that
produces for the line source interference source;
2-3. if interference source is the face Source Type,
When the both sides of emitting antenna and receiving antenna interference source or have one during in the face source in the face source;
When the homonymy of emitting antenna and receiving antenna interference source in the face source
Wherein
is the electromagnetic first break time of interference that face source interference source produces;
Step (3). disturb on identification ground.
Open-air GPR archaeology is surveyed the gained data be carried out to figure, like Fig. 3 (a) and Fig. 5 (a).To calculate gained and disturb T-X curve to become figure, like Fig. 3 (b) and Fig. 5 (b).
Open-air GPR archaeology detection gained data are carried out to figure disturb T-X curve to become figure to compare, and identify ground from measured section and disturb with the calculating gained.
Embodiment 1:
As shown in Figure 2 is the instance that ground object disturbs quick identification when carrying out the detection of GPR archaeology in the Liang Zhu ruins.There are two wire posts on ground, survey line region, and its distance to survey line is respectively 30.5m and 36m, because the distance of their size much smaller than them to survey line, so can regard the point source interference source as.Can obtain the location parameter of two point source interference sources with reference to step 1-1; The parameter of interference source 1 is
; The parameter of interference source 2 is
; Two curves have been produced through calculating two point source interference sources; Like Fig. 3 (b); And two curve lineups are also arranged in the GPR archaeology measured section relevant position of Fig. 3 (a), should be the Electromagnetic Interference that these two wire columns produce.
Embodiment 2:
As shown in Figure 4 be another carry out in the Liang Zhu ruins GPR archaeology when surveying ground object disturb the instance of quick identification.The survey line region, on have an electric wire (the substantially parallel face of land) in vain, can regard the line source interference source as; Also have a building on the ground, can regard face source interference source as.Can obtain the location parameter
of line source interference source with reference to step 1-2; Can obtain the location parameter
of face source interference source with reference to step 1-3; Produced curve through calculating the line source interference source; Like Fig. 5 (b) is the oblique line that face source interference source produces; And also can find curve lineups and oblique line lineups in the GPR archaeology measured section relevant position of Fig. 5 (a); Wherein the curve lineups should be the Electromagnetic Interference from aerial electric wire, and the oblique line lineups are the Electromagnetic Interference of surface structures.
Claims (1)
1. the GPR archaeology is surveyed ground object and is disturbed method for quickly identifying, it is characterized in that following steps:
Step (1). obtain the position of interference source parameter;
1-1. when interference source is the some Source Type; Point source interference source Q point coordinate is (
); Obtain point source interference source Q subpoint on survey line, and measure the position a of subpoint on survey line; Subpoint is to the distance b of survey line on the ground to measure point source interference source Q, and measurement point source interference source Q is to the vertical range c on ground;
1-2 When the interference source for the line type, line source interference
used at any point on a line source
(
) and with a non-zero vector
line parallel expression; online source interference sources
office to take some distance from two o'clock
and
, and then refer to step 1-1 point source interference source Q parameter acquisition obtained two point coordinate
(
) and
(
), a non-zero vector
, the
1-3. when interference source is the face Source Type; The equation of face source interference source S is
, and
is not 0 simultaneously; On the interference source S of face source, choose three the not point of conllinear
,
and
arbitrarily; Obtain coordinate
(
),
(
) and
(
) then with reference to step 1-1 point source interference source Q parameter acquiring mode at 3, and obtain 4 parameters of face source interference source S through following system of equations:
According to the relevant mathematical theorem of plane equation; Can establish
in the above-mentioned system of equations,
for arbitrarily not being 0 constant;
Step (2). calculate and disturb T-X curve;
The GPR archaeology is surveyed and is adopted profile method to measure often, and the relative position of emitting antenna and receiving antenna remains unchanged, along the same moved further of survey line; The antenna distance of emitting antenna and receiving antenna is
; And the line of emitting antenna and receiving antenna is
with
axle clamp angle;
; If the emitting antenna coordinate is (
), the receiving antenna coordinate is (
);
2-1. if interference source is a Source Type,
Wherein
is the skyborne velocity of propagation of electromagnetic wave,
the electromagnetic first break time of interference that produces for the point source interference source;
2-2. if interference source is the line source type,
Wherein: the electromagnetic first break time of interference that
produces for the line source interference source;
2-3. if interference source is the face Source Type,
When the both sides of emitting antenna and receiving antenna interference source or have one during in the face source in the face source;
When the homonymy of emitting antenna and receiving antenna interference source in the face source,
Wherein
is the electromagnetic first break time of interference that face source interference source produces;
Step (3). disturb on identification ground;
Open-air GPR archaeology detection gained data are carried out to figure disturb T-X curve to become figure to compare, and identify interference from measured section with the calculating gained.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011103701048A CN102435866A (en) | 2011-11-21 | 2011-11-21 | Method for quickly identifying interference of ground object during archaeological detection of ground penetrating radar |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011103701048A CN102435866A (en) | 2011-11-21 | 2011-11-21 | Method for quickly identifying interference of ground object during archaeological detection of ground penetrating radar |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102435866A true CN102435866A (en) | 2012-05-02 |
Family
ID=45984019
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011103701048A Pending CN102435866A (en) | 2011-11-21 | 2011-11-21 | Method for quickly identifying interference of ground object during archaeological detection of ground penetrating radar |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102435866A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103376443A (en) * | 2013-07-09 | 2013-10-30 | 浙江大学 | Ground penetrating radar terrestrial interference detecting and fast eliminating method |
CN107356985A (en) * | 2017-08-24 | 2017-11-17 | 浙江大学 | The method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion |
CN108646229A (en) * | 2018-06-14 | 2018-10-12 | 北京师范大学 | Underground column reflector inclination angle detection method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002086542A2 (en) * | 2001-04-20 | 2002-10-31 | Witten Alan J | A method and apparatus for identifying buried objects using ground penetrating radar |
WO2005093462A1 (en) * | 2004-03-24 | 2005-10-06 | Ids Ingegneria Dei Sistemi S.P.A. | Impulsive multi-channel ground penetrating radar |
CN102012509A (en) * | 2010-09-28 | 2011-04-13 | 浙江华东工程安全技术有限公司 | Radar detection method for occurrence of geological structure surface |
-
2011
- 2011-11-21 CN CN2011103701048A patent/CN102435866A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002086542A2 (en) * | 2001-04-20 | 2002-10-31 | Witten Alan J | A method and apparatus for identifying buried objects using ground penetrating radar |
WO2002086542A3 (en) * | 2001-04-20 | 2003-03-27 | Alan J Witten | A method and apparatus for identifying buried objects using ground penetrating radar |
WO2005093462A1 (en) * | 2004-03-24 | 2005-10-06 | Ids Ingegneria Dei Sistemi S.P.A. | Impulsive multi-channel ground penetrating radar |
CN102012509A (en) * | 2010-09-28 | 2011-04-13 | 浙江华东工程安全技术有限公司 | Radar detection method for occurrence of geological structure surface |
Non-Patent Citations (2)
Title |
---|
ANDREW T.CHAMBERLAIN, ECT: "Cave Detection in Limestone using Ground Penetrating Radar", 《JOURNAL OF ARCHAEOLOGICAL SCIENCE》 * |
林金鑫: "综合地球物理考古方法的应用可行性研究", 《中国博士学位论文全文数据库 哲学与人文科学辑(月刊)》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103376443A (en) * | 2013-07-09 | 2013-10-30 | 浙江大学 | Ground penetrating radar terrestrial interference detecting and fast eliminating method |
CN107356985A (en) * | 2017-08-24 | 2017-11-17 | 浙江大学 | The method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion |
CN108646229A (en) * | 2018-06-14 | 2018-10-12 | 北京师范大学 | Underground column reflector inclination angle detection method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11960047B1 (en) | Locating devices, systems, and methods using frequency suites for utility detection | |
US20140012505A1 (en) | Multi-component electromagnetic prospecting apparatus and method of use thereof | |
US11118912B2 (en) | Magnetic positioning system | |
CN102508308A (en) | Transient electromagnetic multipoint array-type detection method of tunnel | |
CN105093299A (en) | Observation system optimization method based on offset vector tile technology and apparatus thereof | |
US10222465B2 (en) | Magnetic field detector and ground-penetrating radar device with merged display | |
CN108076511A (en) | A kind of indoor wireless positioner and method based on low-frequency AC magnetic fields | |
Liu et al. | Estimation of subsurface cylindrical object properties from GPR full-waveform inversion | |
RU2656281C1 (en) | Method of applying swarm of unmanned aerial vehicles for remote determination of location of underground communications, their cross section and depth in ground | |
CN102435866A (en) | Method for quickly identifying interference of ground object during archaeological detection of ground penetrating radar | |
CN105388470B (en) | A kind of method of estimating target motion parameter | |
CN103376443B (en) | Ground penetrating radar terrestrial interference detecting and fast eliminating method | |
CN104391332A (en) | Shallow sea double-frequency controllable source electromagnetic prospecting method | |
CN113341467A (en) | Mine transient electromagnetic three-dimensional display method based on multi-interpolation method | |
CN102830424A (en) | Method for calculating combination parameters of detector | |
CN103969647A (en) | Sparse array architectural layout imaging method based on image domain filtering | |
Eröss et al. | Interpretation of very low frequency measurements carried out with an unmanned aerial system by 2D conductivity models | |
Pagani et al. | A study of HF transmitter geolocation through single-hop ionospheric propagation | |
WO2021098667A1 (en) | Full-coverage anti-interference controlled-source electromagnetic detection method | |
Yanju et al. | Interpretation research on electrical source of time domain ground-airborne electromagnetic data | |
RU2559565C2 (en) | Method of determining spatial position of extended objects located at depth, primarily under water, and electromagnetic line locator, primarily ship electromagnetic line locator for carrying out said method | |
CN108519622B (en) | Underground electric target detection method and device based on natural field source excitation | |
CN117741784B (en) | Ground-air cooperative electromagnetic exploration system and method with air magnetic reference channel | |
Novo | Ground-penetrating radar (GPR) | |
CN111856594B (en) | Tunnel three-direction video dispersion rate advanced detection method |
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: 20120502 |