CN102323623A - GPR Attribute Analysis Method in Archaeological Survey - Google Patents
GPR Attribute Analysis Method in Archaeological Survey Download PDFInfo
- Publication number
- CN102323623A CN102323623A CN201110232249A CN201110232249A CN102323623A CN 102323623 A CN102323623 A CN 102323623A CN 201110232249 A CN201110232249 A CN 201110232249A CN 201110232249 A CN201110232249 A CN 201110232249A CN 102323623 A CN102323623 A CN 102323623A
- Authority
- CN
- China
- Prior art keywords
- attribute
- data
- gpr
- archaeological
- extracts
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 10
- 238000011835 investigation Methods 0.000 claims abstract description 5
- 238000005457 optimization Methods 0.000 claims abstract 2
- 239000000284 extract Substances 0.000 claims description 12
- 238000000034 method Methods 0.000 claims description 11
- 238000000605 extraction Methods 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 4
- 238000011160 research Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 2
- 238000001303 quality assessment method Methods 0.000 claims description 2
- 238000011157 data evaluation Methods 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 3
- 238000009412 basement excavation Methods 0.000 abstract description 2
- 230000000149 penetrating effect Effects 0.000 abstract 4
- 230000000704 physical effect Effects 0.000 abstract 1
- 238000005516 engineering process Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 2
- 238000011002 quantification Methods 0.000 description 2
- 238000013480 data collection Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a ground penetrating radar attribute analysis method in archaeological investigation. Attribute characteristics of describable quantitative electromagnetic wave geometry, kinematics, dynamics, statistics and the like are extracted from ground penetrating radar data, detailed information contained in the data is subjected to deep excavation, information of structure, physical property and the like of an underground archaeological object body are depicted and described. According to the invention, through extracting ground penetrating radar attributes of optimization combination, physical and mathematical characteristics of an electromagnetic wave recorded by a ground penetrating radar can be combined with different archaeological objects well, the archaeological objects are subjected to three-dimensional depiction, and a directive function is provided for detection of a cultural sites unknown object body.
Description
Technical field
The invention belongs to the near surface geophysics field, relate to the GPR property analysis method in a kind of archaeological investigation.
Background technology
GPR is a high-resolution near surface geophysics technology; But further requirement along with archaeology detection of a target fineness; Only cook time slice at original 3-D data volume; Can not satisfy the demand, this just requires us from original 3-D data volume, further to extract the information of hiding.Many similaritys according to GPR and seismic method data acquisition and data interpretation technology; Can use for reference the theory of seismic properties; But extract the attributive character of describable quantification in the Radar Data, and then the structure of portrayal and description buried target body is with information such as transitivities.
Summary of the invention
The objective of the invention is deficiency, the property analysis method of the GPR in a kind of archaeological investigation is provided to prior art.This method is to utilize GPR three-dimensional properties figure that the structure of underground archaeology objective body is carried out meticulous depiction with transitivity.
After the present invention carries out quality assessment and handles field data; But from the GPR data, extract the attributive character such as electromagnetic wave geometry, kinematics, dynamics and statistics of describable quantification, the detailed information of containing in the further excavation data.Make full use of Given information; Especially open-air results of drilling; Multiple attribute is carried out analysis-by-synthesis; And then the GPR attribute is optimized combination, carry out qualitative to the multiple attribute after the optimal combination at last and demarcate quantitatively, thereby realize utilizing GPR three-dimensional properties figure the meticulous depiction of underground archaeology objective body structure with transitivity.
The present invention can extract from the original 3-D data volume of GPR and hide Info, and the archaeology target is carried out three-dimensional meticulous depiction, and then help the detection to site of ancient culture remains unknown object body.
Embodiment
The present invention needs high-quality field data; And the field data collection of standard is the prerequisite that guarantees quality data; In addition; Before the Application Radar attribute technology, the reply data are carried out multianalysis and evaluation (like signal to noise ratio (S/N ratio), borehole data, geological information), and the reliability of data should obtain enough attention.Before data processing, the arrangement that will standardize usually comprises a series of previous works such as unified road header, polar switching, zero correction, merged file.Conventional processing comprises dewow filtering, digital filtering, velocity analysis, skew playback etc.
The object of attributes extraction is divided into two-dimensional radar data and 3-D data volume.The method for distilling that carries out attribute mainly is divided into the extraction of three-dimensional properties body, the reflection horizon bit attribute extracts two aspects.It is to be the basis with the 3-D data volume that described three-dimensional properties body extracts; Extract on the space not people having a common goal's GPR electromagnetic wave attribute; As extracting traditional three wink attributes instantaneous amplitude, instantaneous phase and instantaneous frequency one by one through multiple recording channel analysis, the mathematical principle through infinitesimal analysis extracts road integration, line integral and road differential information.Described reflection horizon bit attribute extracts and is divided into two kinds of methods: a kind of method is that zone of interest is carried out tracing of horizons, opens window for the moment up and down along zone of interest, and then the reflection horizon attribute in the window when extracting; Another kind method is to the reflection horizon, the attribute during set time position, and then extraction up and down in the window.For solving different problems, should be to different archaeology objective bodies, different attributes extraction method, the different attributive character of extraction of application.The attributive character of different target body is represented with attributed graph.
Different objective body reflection wave characteristic is different, so when before concrete attributes extraction, the information of study area being carried out comprehensive as far as possible understanding, also will be based on full appreciation to the reflection wave characteristic.A large amount of attributes are extracted out; But some attribute effect in practical application is unsatisfactory, and also comprises the composition that is relative to each other in a large amount of attributes, causes the repetition and the waste of bulk information; Need be optimized attribute, and optimize the attribute the most responsive finding the solution problem.According to known achievement in research, study area geology characteristic, combine the rerum natura and the structural information of the physics of GPR institute recording electromagnetic wave and mathematical characteristic and goal in research body at last, carry out qualitative and demarcation quantitatively the multiple attribute of GPR.
Claims (1)
1. the GPR property analysis method in the archaeological investigation is characterized in that this method may further comprise the steps:
Step 1. data evaluation and routine data are handled; Specifically: the data of gathering are carried out quality assessment according to signal to noise ratio (S/N ratio), borehole data and geological information; Carry out corresponding routine data according to demand then and handle, comprise unified road header, polar switching, zero correction, merged file, dewow filtering, digital filtering, velocity analysis and skew playback successively;
The attributes extraction of step 2. GPR, specifically: the object of attributes extraction is divided into two-dimensional radar data and 3-D data volume, and the method for distilling that carries out attribute information is divided into the extraction of three-dimensional properties body, the reflection horizon bit attribute extracts two aspects; It is to be the basis with the 3-D data volume that described three-dimensional properties body extracts, and extracts on the space not people having a common goal's GPR electromagnetic wave attribute; Described reflection horizon bit attribute extracts and is divided into two kinds of methods: a kind of method is that zone of interest is carried out tracing of horizons, opens window for the moment up and down along zone of interest, and then the reflection horizon attribute in the window when extracting; Another kind method is to the reflection horizon, the attribute during set time position, and then extraction up and down in the window;
The attribute optimization of step 3. GPR and attribute are demarcated, specifically: based on study area information, according to Principle of Statistics, from the attribute that is extracted, optimize the attribute the most responsive to finding the solution problem; According to known achievement in research, study area geology characteristic, goal in research physical difference, carry out qualitative and demarcation quantitatively to the multiple attribute of GPR.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110232249A CN102323623A (en) | 2011-08-15 | 2011-08-15 | GPR Attribute Analysis Method in Archaeological Survey |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201110232249A CN102323623A (en) | 2011-08-15 | 2011-08-15 | GPR Attribute Analysis Method in Archaeological Survey |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN102323623A true CN102323623A (en) | 2012-01-18 |
Family
ID=45451402
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201110232249A Pending CN102323623A (en) | 2011-08-15 | 2011-08-15 | GPR Attribute Analysis Method in Archaeological Survey |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102323623A (en) |
Cited By (4)
| 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 |
| CN107167845A (en) * | 2017-05-23 | 2017-09-15 | 浙江大学 | Improve the system and method for magnetic archaeology Data Gathering efficiency |
| CN107356985A (en) * | 2017-08-24 | 2017-11-17 | 浙江大学 | The method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion |
| CN108287369A (en) * | 2017-12-08 | 2018-07-17 | 中国煤炭地质总局勘查研究总院 | A kind of isochronous surface acquisition methods and device based on three-dimensional Ground Penetrating Radar attribute volume |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1191594A (en) * | 1995-06-19 | 1998-08-26 | 弗米尔制造公司 | Excavator data acquisition and control system and process |
| US20040133347A1 (en) * | 2003-12-08 | 2004-07-08 | Britt John T. | Automated resource management system (ARMSTM) |
| WO2005093462A1 (en) * | 2004-03-24 | 2005-10-06 | Ids Ingegneria Dei Sistemi S.P.A. | Impulsive multi-channel ground penetrating radar |
| CN1896768A (en) * | 2005-07-14 | 2007-01-17 | 中国石油化工股份有限公司 | Treatment for displaying three-dimensional earthquake data reflecting structure |
| WO2008053489A2 (en) * | 2006-11-03 | 2008-05-08 | Paradigm Geophysical (Luxembourg) S.A.R.L. | System and method for full azimuth angle domain imaging in reduced dimensional coordinate systems |
| CN101669043A (en) * | 2007-03-08 | 2010-03-10 | 格库技术有限公司 | Systems and methods for seismic data acquisition employing asynchronous, decoupled data sampling and transmission |
-
2011
- 2011-08-15 CN CN201110232249A patent/CN102323623A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1191594A (en) * | 1995-06-19 | 1998-08-26 | 弗米尔制造公司 | Excavator data acquisition and control system and process |
| US20040133347A1 (en) * | 2003-12-08 | 2004-07-08 | Britt John T. | Automated resource management system (ARMSTM) |
| WO2005093462A1 (en) * | 2004-03-24 | 2005-10-06 | Ids Ingegneria Dei Sistemi S.P.A. | Impulsive multi-channel ground penetrating radar |
| CN1896768A (en) * | 2005-07-14 | 2007-01-17 | 中国石油化工股份有限公司 | Treatment for displaying three-dimensional earthquake data reflecting structure |
| WO2008053489A2 (en) * | 2006-11-03 | 2008-05-08 | Paradigm Geophysical (Luxembourg) S.A.R.L. | System and method for full azimuth angle domain imaging in reduced dimensional coordinate systems |
| CN101669043A (en) * | 2007-03-08 | 2010-03-10 | 格库技术有限公司 | Systems and methods for seismic data acquisition employing asynchronous, decoupled data sampling and transmission |
Non-Patent Citations (4)
| Title |
|---|
| 《中国地球物理2010--中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集》 20101231 赵文轲，田钢 探地雷达三维属性技术在考古调查中的应用研究 , * |
| 林金鑫,田钢等: "良渚遗址古水系调查中的综合地球物理方法", 《浙江大学学报(工学版)》 * |
| 谢昭晖等: "我国探地雷达的应用现状及展望", 《工程勘察》 * |
| 赵文轲，田钢: "探地雷达三维属性技术在考古调查中的应用研究", 《中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集》 * |
Cited By (4)
| 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 |
| CN107167845A (en) * | 2017-05-23 | 2017-09-15 | 浙江大学 | Improve the system and method for magnetic archaeology Data Gathering efficiency |
| CN107356985A (en) * | 2017-08-24 | 2017-11-17 | 浙江大学 | The method that earthen ruins archaeological investigation is carried out using GPR wave impedance inversion |
| CN108287369A (en) * | 2017-12-08 | 2018-07-17 | 中国煤炭地质总局勘查研究总院 | A kind of isochronous surface acquisition methods and device based on three-dimensional Ground Penetrating Radar attribute volume |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104502969A (en) | Channel sandstone reservoir identification method | |
| Wood et al. | An interpretation of surface and borehole seismic surveys for mine planning at the Millennium uranium deposit, northern Saskatchewan, Canada | |
| Fernández-Álvarez et al. | Combined GPR and ERT exploratory geophysical survey of the Medieval Village of Pancorbo Castle (Burgos, Spain) | |
| CN101000378A (en) | Method for determining air layer using acoustic wave split-double pole transverse wave well-logging data | |
| CN101625413A (en) | Oil-gas exploration method | |
| Farrugia et al. | Noninvasive techniques for site characterization of Alberta seismic stations based on shear‐wave velocity | |
| Khalil et al. | Engineering applications of seismic refraction method: A field example at Wadi Wardan, Northeast Gulf of Suez, Sinai, Egypt | |
| Yaede et al. | A geophysical strategy for measuring the thickness of the critical zone developed over basalt lavas | |
| CN102323623A (en) | GPR Attribute Analysis Method in Archaeological Survey | |
| Chandler et al. | Utility of the horizontal-to-vertical spectral ratio passive seismic method for estimating thickness of Quaternary sediments in Minnesota and adjacent parts of Wisconsin | |
| CN106443791B (en) | The method for asking for tilted stratum or anisotropic formation shear wave Value of residual static correction | |
| CN102096098A (en) | Seismic attribute extracting method | |
| CN106033127A (en) | Ground stress orientation seismic prediction method based on shear wave speed variation rate | |
| Goldie | A comparison between conventional and distributed acquisition induced polarization surveys for gold exploration in Nevada | |
| Lowe et al. | Exploring ground-penetrating radar and sediment magnetic susceptibility analyses in a sandstone rockshelter in northern Australia | |
| CN113487734B (en) | Three-dimensional geological model construction method based on field outcrop and ground penetrating radar | |
| Urban et al. | Geophysical investigation of a Middle Holocene archaeological site along the Yukon River, Alaska | |
| Shi et al. | Identification and distribution of fractures in the Zhangjiatan shale of the Mesozoic Yanchang Formation in Ordos Basin | |
| Wang et al. | 3D amplitude-preserved full-azimuth anisotropic imaging-offering more reliable azimuthal AVO for fracture detection | |
| Łyskowski et al. | Georadar investigations and geochemical analysis in contemporary archeological studies | |
| Menezes et al. | High-resolution facies modeling of presalt lacustrine carbonates reservoir analog: Morro do Chaves Formation example, Sergipe-Alagoas Basin, Brazil | |
| Cheng et al. | Geophysical case study of the Iso and new Insco deposits, Québec, Canada, part I: data comparison and analysis | |
| Reilly et al. | Data report: paleomagnetic directions from IODP Expedition 354, Hole U1451A, Cores 23H and 24H | |
| Mustaqeem et al. | Application of the Gray Level Co-occurrence Matrix (GLCM) Technique for Analyzing Fracture Networks in the Granitic Basement of Cuu Long Basin, Vietnam | |
| Strobach et al. | Cross well radar and vertical radar profiling methods for time lapse monitoring of rainfall infiltration |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20120118 |
|
| RJ01 | Rejection of invention patent application after publication |