CN107356985A - The method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion - Google Patents

The method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion Download PDF

Info

Publication number
CN107356985A
CN107356985A CN201710736517.0A CN201710736517A CN107356985A CN 107356985 A CN107356985 A CN 107356985A CN 201710736517 A CN201710736517 A CN 201710736517A CN 107356985 A CN107356985 A CN 107356985A
Authority
CN
China
Prior art keywords
gpr
wave impedance
dielectric constant
record
core
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
CN201710736517.0A
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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201710736517.0A priority Critical patent/CN107356985A/en
Publication of CN107356985A publication Critical patent/CN107356985A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/885Radar or analogous systems specially adapted for specific applications for ground probing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/887Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a kind of method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion.It carries out the data acquisition of GPR common offset and drilling extracting core first, and dielectric constant values are measured and record at regular intervals to obtain well data to the core of taking-up;Then, Coherent Noise in GPR Record is pre-processed, and well data is processed to obtain a low-frequency impedance body;Wave impedance inversion is done using the GPR after processing and well data, the wave impedance result of Whole frequency band is obtained, wave impedance is converted into dielectric constant, obtains a high-resolution underground dielectric constant profile;Contrast the dielectric constant values of the core of some positions of record, it can be deduced that the distribution characteristics of the different archaeology targets in underground.The invention provides the detailed distributed intelligence of the high-resolution physical parameter of underground medium, the distribution situation of dirt archaeology target can be accurately inferred to using physical parameter information, so as to improve detectivity of the GPR to target of engaging in archaeological studies.

Description

The method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion
Technical field
The present invention relates to Agrocybe praecox to investigate field, more particularly to one kind utilizes GPR wave impedance inversion to carry out native something lost The method of location archaeological investigation.
Background technology
Ground penetrating radar exploration has the advantages of quick nondestructive and higher resolution ratio, afield wide in archaeological investigation at present General application.However, simply carrying out conventional treatment to Coherent Noise in GPR Record at present, its result can only provide the position form such as interface So limited information, and formation physical parameters distribution can not be provided, seriously constrain the differentiation to dirt archaeology destination properties.
Drilling extracting core is the effective earthen ruins archaeological investigation means of a tradition, and drilling extracting core construction is carried out using Luoyang Spade Simply and to ruins destructiveness very little, the core taken out by analyzing drilling well, mesh of engaging in archaeological studies can be accurately determined in vertical direction Target property, but be often a peephole view, laterally complete earthen ruins archaeological investigation result can not be obtained.
Therefore, either GPR or drilling extracting core method, application effect is all by serious restriction.
The content of the invention
The purpose of the present invention is overcome the deficiencies in the prior art, there is provided a kind of to carry out soil using GPR wave impedance inversion The method of ruins archaeological investigation.
The technical solution adopted in the present invention is as follows:It is a kind of to carry out earthen ruins archaeology tune using GPR wave impedance inversion The method looked into, this method comprise the following steps:
(1) GPR and drilling data acquisition:
According to archaeology target design survey line, Coherent Noise in GPR Record is gathered using common offset mode, if being chosen on survey line Dry position carries out drilling extracting core, measures and record the dielectric constant values of core at regular intervals;
(2) GPR and well data processing:
The Coherent Noise in GPR Record collected in step (1) is pre-processed, the dielectric of the core measured in step (1) is normal Number is converted to wave impedance, and a low-frequency impedance body is obtained by interpolation processing;
(3) wave impedance inversion and data interpretation:
Pretreated Coherent Noise in GPR Record and low-frequency impedance body are merged using wave impedance inversion, obtain one entirely The wave impedance result of frequency band, dielectric constant is converted to by wave impedance, you can obtains a high-resolution underground dielectric constant point Cloth section, the dielectric constant values of the core of some positions of contrast step (1) record, you can analyze the different archaeology targets in underground Distribution situation.
Further, the pretreatment in the step (2) includes dewow, amplitude compensation, horizontal filtering, deconvolution, band logical Filtering and migration processing.
Compared with prior art, the present invention is merged GPR and well data using wave impedance inversion, as a result The detailed distributed intelligence of high-resolution physical parameter can be provided, so as to add the differentiation to earthen ruins archaeology destination properties Ability;The function of boring coring is extended, it is organically combined with Coherent Noise in GPR Record, solves its transverse direction The problem of loss of learning;Whole flow process is more time saving and energy saving, and place destructiveness is smaller, improves GPR and afield engages in archaeological studies The detectivity in investigation field.
Brief description of the drawings
Fig. 1 is the flow chart that archaeological investigation is carried out using GPR wave impedance inversion;
Fig. 2 is the schematic diagram for designing survey line and bore position;
Fig. 3 is the result schematic diagram of wavelet transform;
Fig. 4 is the result schematic diagram of well data processing;
Fig. 5 is the result schematic diagram of GPR wave impedance inversion.
Embodiment
The present invention is described further with reference to the accompanying drawings and examples.
As shown in figure 1, the present embodiment offer is a kind of to carry out earthen ruins archaeological investigation using GPR wave impedance inversion Method, this method comprise the following steps:
1) GPR and drilling data acquisition
This example archaeology target is an ancient stream channel, designs what a 25m vertical with ancient stream channel direction grew according to Given information Survey line, as shown in Figure 2.Spy land mine is carried out in a manner of common offset using PulseEKKO Pro systems and supporting 100MHz antennas Up to data acquisition, antenna spacing 1m, step-length 0.1m, sampling interval 0.4ns, 64 superpositions of per pass;CMP is done on survey line Velocity analysis, it is 0.075m/ns to obtain reference velocity;According to the Coherent Noise in GPR Record feature collected, 5m, 15m, 25m are chosen Three positions carry out drilling extracting core using Luoyang Spade, drilling depth 2m, have all detected raw soil layer, and archaeology point is done to core Analysis, the dielectric constant values and record of core are V.7 measured using Percometer according to 0.1m interval;
2) GPR and well data processing
Dewow and amplitude compensation, use level are done to Coherent Noise in GPR Record using the Reflexw softwares that version number is 5.0 Filtering removes direct wave, then the resolution ratio and signal to noise ratio of data are lifted with Least square deconvolution and bandpass filtering, is offset with f-k Reflective information is accurately playbacked;Time and depth transfer is done according to reference velocity, and cuts off 2m data below, remaining data is changed It is as shown in Figure 3 to return the result that time-domain obtains;The dielectric constant data that well logging obtains is converted into wave impedance data, to every mouthful of well Data are done one-line interpolation and obtained with GPR per track data identical sampled point, then do two-dimensional linear to three mouthfuls of well data Interpolation extrapolation, and do using in soil deeply and obtain low-frequency impedance body such as Fig. 4 in time-domain;
3) wave impedance inversion and data interpretation
The low-frequency impedance body obtained using the Coherent Noise in GPR Record after processing and well data, using wave impedance inversion to visiting Ground radar and well data are merged, and obtain the wave impedance result of a Whole frequency band;Wave impedance is converted into dielectric constant, can be obtained To a high-resolution underground dielectric constant profile result such as Fig. 5;Analyzed and measured according to core before, learn river course Alluvial has higher dielectric constant than the raw-soil, so as to clearly recognize the lower interface of ancient stream channel such as the dotted line institute in Fig. 5 Show, and the vertically and horizontally detailed distribution situation of ancient stream channel alluvial can be recognized accurately.
To sum up apparently, compared to traditional GPR archaeological investigation method for only providing the limited informations such as interface location, sheet Method is merged GPR and well data using wave impedance inversion, as a result can provide high-resolution physical parameter Detailed distributed intelligence, so as to add to earthen ruins archaeology destination properties discriminating power;The function of boring coring is carried out Extension, makes it organically combine with Coherent Noise in GPR Record, solves the problems, such as its horizontal loss of learning;Whole flow process compared with To be time saving and energy saving, and place destructiveness is smaller, improves the detectivity in GPR afield archaeological investigation field.

Claims (2)

  1. A kind of 1. method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion, it is characterised in that this method bag Include following steps:
    (1) GPR and drilling data acquisition:
    According to archaeology target design survey line, Coherent Noise in GPR Record is gathered using common offset mode, some positions are chosen on survey line Carry out drilling extracting core is put, measures and record the dielectric constant values of core at regular intervals.
    (2) GPR and well data processing:
    The Coherent Noise in GPR Record collected in step (1) is pre-processed, the dielectric constant of the core measured in step (1) is turned Wave impedance is changed to, a low-frequency impedance body is obtained by interpolation processing.
    (3) wave impedance inversion and data interpretation:
    Pretreated Coherent Noise in GPR Record and low-frequency impedance body are merged using wave impedance inversion, obtain a Whole frequency band Wave impedance result, wave impedance is converted into dielectric constant, you can obtain a high-resolution underground dielectric constant distribution and cut open Face, the dielectric constant values of the core of some positions of contrast step (1) record, you can analyze point of the different archaeology targets in underground Cloth situation.
  2. 2. a kind of method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion according to claim 1, Characterized in that, pretreatment in the step (2) include dewow, amplitude compensation, horizontal filtering, deconvolution, bandpass filtering and Migration processing.
CN201710736517.0A 2017-08-24 2017-08-24 The method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion Pending CN107356985A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710736517.0A CN107356985A (en) 2017-08-24 2017-08-24 The method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710736517.0A CN107356985A (en) 2017-08-24 2017-08-24 The method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion

Publications (1)

Publication Number Publication Date
CN107356985A true CN107356985A (en) 2017-11-17

Family

ID=60289257

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710736517.0A Pending CN107356985A (en) 2017-08-24 2017-08-24 The method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion

Country Status (1)

Country Link
CN (1) CN107356985A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111239726A (en) * 2020-02-26 2020-06-05 山东超越数控电子股份有限公司 Method, device and readable storage medium for identifying underground pollutant area
CN113376629A (en) * 2021-05-13 2021-09-10 电子科技大学 In-well radar least square inversion method based on non-uniform input parameter grid

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005093462A1 (en) * 2004-03-24 2005-10-06 Ids Ingegneria Dei Sistemi S.P.A. Impulsive multi-channel ground penetrating radar
CN102323623A (en) * 2011-08-15 2012-01-18 浙江大学 Ground penetrating radar attribute analysis method in archaeological investigation
CN102435866A (en) * 2011-11-21 2012-05-02 浙江大学 Method for quickly identifying interference of ground object during archaeological detection of ground penetrating radar

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005093462A1 (en) * 2004-03-24 2005-10-06 Ids Ingegneria Dei Sistemi S.P.A. Impulsive multi-channel ground penetrating radar
CN102323623A (en) * 2011-08-15 2012-01-18 浙江大学 Ground penetrating radar attribute analysis method in archaeological investigation
CN102435866A (en) * 2011-11-21 2012-05-02 浙江大学 Method for quickly identifying interference of ground object during archaeological detection of ground penetrating radar

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
ZENG ZHAO-FA 等: "Recursive impedance inversion of ground-penetrating radar data in stochastic media", 《APPLIED GEOPHYSICS》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111239726A (en) * 2020-02-26 2020-06-05 山东超越数控电子股份有限公司 Method, device and readable storage medium for identifying underground pollutant area
CN113376629A (en) * 2021-05-13 2021-09-10 电子科技大学 In-well radar least square inversion method based on non-uniform input parameter grid
CN113376629B (en) * 2021-05-13 2022-08-05 电子科技大学 In-well radar least square inversion method based on non-uniform input parameter grid

Similar Documents

Publication Publication Date Title
Zenone et al. Preliminary use of ground-penetrating radar and electrical resistivity tomography to study tree roots in pine forests and poplar plantations
Aizebeokhai et al. The use of the multiple-gradient array for geoelectrical resistivity and induced polarization imaging
Xiao et al. Permafrost subgrade condition assessment using extrapolation by deterministic deconvolution on multifrequency GPR data acquired along the Qinghai-Tibet railway
Wood et al. An interpretation of surface and borehole seismic surveys for mine planning at the Millennium uranium deposit, northern Saskatchewan, Canada
Wang et al. The automatic interpretation of structural plane parameters in borehole camera images from drilling engineering
CN108241180A (en) Ion type rareearth ore bottom plate exploration method
CN107356985A (en) The method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion
Anbazhagan et al. Identification of karstic features in lateritic soil by an integrated geophysical approach
Zhao et al. Application of 3D GPR attribute technology in archaeological investigations
Adepelumi et al. Geophysical mapping of subsurface stratigraphy beneath a river bed using ground penetrating radar: Lagos Nigeria case study
Afshar et al. Using ground‐penetrating radar to explore the cemented soil horizon in an arid region in Iran
Shaikh et al. Frequency band adjustment match filtering based on variable frequency GPR antennas pairing scheme for shallow subsurface investigations
Sarhan Possibility of intrusive igneous body beneath the Cretaceous sequence in Abu Gharadig Basin, Egypt: integration of geophysical data interpretations
Dhakate et al. Application of ground-penetrating radar for identification of groundwater resources in a coastal terrain
CN115857047B (en) Comprehensive prediction method for earthquake reservoir
CN102323623A (en) Ground penetrating radar attribute analysis method in archaeological investigation
Coutinho et al. Applications of geophysical techniques to improve a groundwater conceptual model in an outcrop area of the Guarani Aquifer System, in Brazil
CN104863578A (en) Electrical conductive rate well measuring method of radar in well
Knoll et al. Characterization of a sand and gravel aquifer using ground-penetrating radar, Cape Cod, Massachusetts
Seshunarayana et al. Multichannel analysis of surface waves (MASW) for mapping shallow subsurface layers—a case study, Jabalpur, India
Zeng et al. Three-dimensional seismic analysis of a polygonal fault system (PFS) in the Northern Carnarvon Basin, Australia: Implications for fluid flow migration and gas hydrate system
Pei et al. Contribution of WBH acquisition for reservoir geophysical prospecting on Songliao Basin
Burlakovs et al. Ground-penetrating Radar (GPR) Geoenvironmental Screening in Lakes of Latvia-Challenges and Outcomes
Cheng et al. Geophysical case study of the Iso and new Insco deposits, Québec, Canada, part I: data comparison and analysis
Febriansanu et al. Recent Sedimentary Deposit and Dynamics Sedimentation Analysis at Teluk Betung Using GPR Interpretation and Core Data

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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20171117

WD01 Invention patent application deemed withdrawn after publication