The travel time and amplitude of ground-penetrating radar （GPR） waves are closely related to medium parameters such as water content, porosity, and dielectric permittivity. However, conventional estimation methods, ...The travel time and amplitude of ground-penetrating radar （GPR） waves are closely related to medium parameters such as water content, porosity, and dielectric permittivity. However, conventional estimation methods, which are mostly based on wave velocity, are not suitable for real complex media because of limited resolution. Impedance inversion uses the reflection coefficient of radar waves to directly calculate GPR impedance and other parameters of subsurface media. We construct a 3D multiscale stochastic medium model and use the mixed Gaussian and exponential autocorrelation function to describe the distribution of parameters in real subsurface media. We introduce an elliptical Gaussian function to describe local random anomalies. The tapering function is also introduced to reduce calculation errors caused by the numerical simulation of discrete grids. We derive the impedance inversion workflow and test the calculation precision in complex media. Finally, we use impedance inversion to process GPR field data in a polluted site in Mongolia. The inversion results were constrained using borehole data and validated by resistivity data.展开更多
The englacial structures and ice thickness of the Laohugou No. 12 （L12） Glacier in the Qilian Mountains, China, were retrieved from ground-penetrating radar （GPR） profile dzta acquired in August of 2007. Here the ...The englacial structures and ice thickness of the Laohugou No. 12 （L12） Glacier in the Qilian Mountains, China, were retrieved from ground-penetrating radar （GPR） profile dzta acquired in August of 2007. Here the interpretation of a typical GPR image is validated using two-dimensional, Finite-Difference Time-Domain （FDTD） numerical modeling. Data analyses revealed many en- glacial characteristics, such as temperate ice, crevasses, and cavities at the position of convergence between the eastern and west- ern glacial branches of L12, and at an altitude between 4,600 and 4,750 m a.s.1, on the east branch. Combining ice thickness, en- glacial structures, subglacial topography, and surface flow velocities of this glacier, we analyzed the reasons for the distribution of temperate ice. The results show that greater englacial water content is associated with englacial crevassing and surface moulins, which allow water to be channeled to the temperate ice aquifer beneath the surface cold ice layer. Analysis of air temperature data shows that as more meltwater imports into the ice body, this has a great effect on water conservation and dynamics conditions. With climate warming, and under the influence of crevasses, subglacial structures, and ice thickness, ice thickness reduction on the L12 east branch is more rapid than that on the west branch.展开更多
Ground-penetrating radar (GPR) has been used predominantly for environments with low electrical conductivity like freshwater aquifers, glaciers, or dry sandy soils. The objective of the present study was to explore it...Ground-penetrating radar (GPR) has been used predominantly for environments with low electrical conductivity like freshwater aquifers, glaciers, or dry sandy soils. The objective of the present study was to explore its application for mapping in subsurface agricultural soils to a depth of several meters. For a loamy sand and a clayey site on the North China Plain, clay inclusions in the sand were detected; the thickness, inclination, and continuity of the confining clay and silt layers was assessed; and a local water table was mapped. Direct sampling (soil coring and profiling) in the top meter and independent measurement of the water table were utilized to confirm the findings. Also, effective estimates of the dielectric number for the site with the dielectric number of moist clayey soils depending strongly on frequency were obtained. Thus, important properties of soils, like the arrangement and type of layers and in particular their continuity and inclination, could be explored with moderate efforts for rather large areas to help find optimal locations for the time-consuming and expensive measurements which would be necessary to detail a model of the subsurface.展开更多
In ground-penetrating radar （GPR） imaging, it is common for the depth of investigation to be on the same order as the variability in surface topography, In such cases, migration fails when it is carried out from a d...In ground-penetrating radar （GPR） imaging, it is common for the depth of investigation to be on the same order as the variability in surface topography, In such cases, migration fails when it is carried out from a datum after the application of elevation statics, We introduce a reverse-time migration （RTM） algorithm based on the second-order decoupled form of Maxwell＇s equations, which requires computation of only the electric field, The wavefield extrapolation is computed directly from the acquisition surface without the need for datuming, In a synthetic case study, the algorithm significantly improves image accuracy over a processing sequence in which migration is performed after elevation statics, In addition, we acquired a field dataset at the Coral Pink Sand Dunes （CPSD） in Utah, USA, The data were acquired over rugged topography and have the complex internal stratigraphy of multiply eroded, modern, and ancient eolian deposits, The RTM algorithm significantly improves radar depth images in this challenging environment,展开更多
In order to achieve a rapid and accurate identification of soil stratification information and accelerate the development of smart agriculture,this paper conducted soil stratification experiments on agricultural soils...In order to achieve a rapid and accurate identification of soil stratification information and accelerate the development of smart agriculture,this paper conducted soil stratification experiments on agricultural soils in the Mollisols area of Northeast China using Ground Penetrating Radar(GPR)and obtained different types of soil with frequencies of 500 MHz,250 MHz,and 100 MHz antennas.The soil profile data were obtained for 500 MHz,250 MHz,and 100 MHz antennas,and the dielectric properties of each type of soil were analyzed.In the image processing procedure,wavelet analysis was first used to decompose the pre-processed radar signal and reconstruct the high-frequency information to obtain the reconstructed signal containing the stratification information.Secondly,the reconstructed signal is taken as an envelope to enhance the stratification information.The Hilbert transform is applied to the envelope signal to find the time-domain variation of the instantaneous frequency and determine the time-domain location of the stratification.Finally,the dielectric constant of each soil horizon is used to obtain the propagation velocity of the electromagnetic wave at the corresponding position to obtain the stratification position of each soil horizon.The research results show that the 500 MHz radar antenna can accurately delineate Ap/Ah,horizon and the absolute accuracy of the stratification is within 5 cm.The effect on the soil stratification below the tillage horizon is not apparent,and the absolute accuracy of the 250 MHz and 100 MHz radar antennas on the stratification is within 9 cm.The overwhelming majority of the overall calculation errors are kept to within 15%.Based on the three central frequency antennas,the soil horizon detection rate reaches 93.3%,which can achieve accurate stratification of soil profiles within 1 m.The experimental and image processing methods used are practical and feasible;however,the GPR will show a missed detection for soil horizons with only slight differences in dielectric properties.Overall,this study can quickly and accurately determine the information of each soil stratification,ultimately providing technical support for acquiring soil configuration information and developing smart agriculture.展开更多
A comprehensive Ground Penetration Radar (GPR) investigations and hazard assessment for the rehabilitation and strengthening of Habib Sakakini’s Palace in Cairo is presented herein, which is considered one of the mos...A comprehensive Ground Penetration Radar (GPR) investigations and hazard assessment for the rehabilitation and strengthening of Habib Sakakini’s Palace in Cairo is presented herein, which is considered one of the most significant architectural heritage sites in Egypt. The palace located on an ancient water pond at the eastern side of Egyptian gulf besiding Sultan Bebris Al-Bondoqdary mosque is a place also called “Prince Qraja al-Turkumany pond”. That pond had been filled down by Habib Sakakini at 1892 to construct his famous palace in 1897. The integrated geophysical survey of the palace allowed the identification of several targets of potential archaeological and geotechnical engineering interest buried in fill and silty clay in the depth range between 100 - 700 cm. the methodological development focused on Multi-Fold (MF) Ground Penetrating Radar (GPR) imaging and subsurface characterization based on integrated velocity and attenuation analysis. Eight hundred sqm of Ground penetration Radar (GPR) profiling have been conducted to monitor the subsurface conditions. 600 meters are made in the surrounding area of the Palace and 200 sqm at the basement. The aim is to monitor the soil conditions beneath and around the Palace and to identify potential geological discontinuities, or the presence of faults and cavities. A suitable single and dual antenna are used (500 - 100 MHZ) is used to penetrate the desired depth of 7 meters (ASTM D6432). The GPR is used also detect the water table. At the building basement the GPR is used to identify the foundation thickness and soil-basement interface. As well as the inspection of cracks in some supporting columns, piers and masonry walls. The GPR also was used to investigate the floors and ceilings conditions and structural mapping. The results were validated by the geotechnical and structural surveys. All these results together with the seismic hazard analysis will be used for the complete analysis of the palace in the framework of the rehabilitation and strengthening works foreseen in a second stage.展开更多
探讨了多小波函数及其预处理方法对探地雷达(ground penetrating radar—GPR)图象去噪性能的影响,在Donoho D L和Johnstone I M提出的小波阈值去噪方法的基础上提出了一个改进的阈值函数,并对实际的GPR图象进行阈值化处理和对比分析,结...探讨了多小波函数及其预处理方法对探地雷达(ground penetrating radar—GPR)图象去噪性能的影响,在Donoho D L和Johnstone I M提出的小波阈值去噪方法的基础上提出了一个改进的阈值函数,并对实际的GPR图象进行阈值化处理和对比分析,结果表明选取合适的预处理方法,采用DGHM和STT多小波对GPR图象去噪可获得比其他方法更好的效果。展开更多
基金supported by the Doctoral Fund Project of the Ministry of Education(No.20130061110060 class tutors)the Post-Doctoral Fund Project(No.2015M571366)+1 种基金the National Natural Science Foundation of China(No.41174097)US DoD ARO Project"Advanced Mathematical Algorithm"(No.W911NF-11-2-0046)
文摘The travel time and amplitude of ground-penetrating radar （GPR） waves are closely related to medium parameters such as water content, porosity, and dielectric permittivity. However, conventional estimation methods, which are mostly based on wave velocity, are not suitable for real complex media because of limited resolution. Impedance inversion uses the reflection coefficient of radar waves to directly calculate GPR impedance and other parameters of subsurface media. We construct a 3D multiscale stochastic medium model and use the mixed Gaussian and exponential autocorrelation function to describe the distribution of parameters in real subsurface media. We introduce an elliptical Gaussian function to describe local random anomalies. The tapering function is also introduced to reduce calculation errors caused by the numerical simulation of discrete grids. We derive the impedance inversion workflow and test the calculation precision in complex media. Finally, we use impedance inversion to process GPR field data in a polluted site in Mongolia. The inversion results were constrained using borehole data and validated by resistivity data.
基金funded by a project of the Chinese Glacier Inventory(41071047 and KZCX2-YW-GJ04)
文摘The englacial structures and ice thickness of the Laohugou No. 12 （L12） Glacier in the Qilian Mountains, China, were retrieved from ground-penetrating radar （GPR） profile dzta acquired in August of 2007. Here the interpretation of a typical GPR image is validated using two-dimensional, Finite-Difference Time-Domain （FDTD） numerical modeling. Data analyses revealed many en- glacial characteristics, such as temperate ice, crevasses, and cavities at the position of convergence between the eastern and west- ern glacial branches of L12, and at an altitude between 4,600 and 4,750 m a.s.1, on the east branch. Combining ice thickness, en- glacial structures, subglacial topography, and surface flow velocities of this glacier, we analyzed the reasons for the distribution of temperate ice. The results show that greater englacial water content is associated with englacial crevassing and surface moulins, which allow water to be channeled to the temperate ice aquifer beneath the surface cold ice layer. Analysis of air temperature data shows that as more meltwater imports into the ice body, this has a great effect on water conservation and dynamics conditions. With climate warming, and under the influence of crevasses, subglacial structures, and ice thickness, ice thickness reduction on the L12 east branch is more rapid than that on the west branch.
基金Project supported in part by the Deutsche Forschungsgemeinschaft (DFG), Germany (No. RO 1080/8-1) jointly by Max-Planck Gesellschaft and the Chinese Academy of Sciences through a travel grant to the first author.
文摘Ground-penetrating radar (GPR) has been used predominantly for environments with low electrical conductivity like freshwater aquifers, glaciers, or dry sandy soils. The objective of the present study was to explore its application for mapping in subsurface agricultural soils to a depth of several meters. For a loamy sand and a clayey site on the North China Plain, clay inclusions in the sand were detected; the thickness, inclination, and continuity of the confining clay and silt layers was assessed; and a local water table was mapped. Direct sampling (soil coring and profiling) in the top meter and independent measurement of the water table were utilized to confirm the findings. Also, effective estimates of the dielectric number for the site with the dielectric number of moist clayey soils depending strongly on frequency were obtained. Thus, important properties of soils, like the arrangement and type of layers and in particular their continuity and inclination, could be explored with moderate efforts for rather large areas to help find optimal locations for the time-consuming and expensive measurements which would be necessary to detail a model of the subsurface.
基金The Herbette Foundation at the University of Lausanne provided support for the development of the RTM algorithm
文摘In ground-penetrating radar （GPR） imaging, it is common for the depth of investigation to be on the same order as the variability in surface topography, In such cases, migration fails when it is carried out from a datum after the application of elevation statics, We introduce a reverse-time migration （RTM） algorithm based on the second-order decoupled form of Maxwell＇s equations, which requires computation of only the electric field, The wavefield extrapolation is computed directly from the acquisition surface without the need for datuming, In a synthetic case study, the algorithm significantly improves image accuracy over a processing sequence in which migration is performed after elevation statics, In addition, we acquired a field dataset at the Coral Pink Sand Dunes （CPSD） in Utah, USA, The data were acquired over rugged topography and have the complex internal stratigraphy of multiply eroded, modern, and ancient eolian deposits, The RTM algorithm significantly improves radar depth images in this challenging environment,
基金Under the auspices of the National Key R&D Program of China(No.2021YFD1500100)the Strategic Priority Research Program of the Chinese Academy of Sciences(No.XDA28100000)。
文摘In order to achieve a rapid and accurate identification of soil stratification information and accelerate the development of smart agriculture,this paper conducted soil stratification experiments on agricultural soils in the Mollisols area of Northeast China using Ground Penetrating Radar(GPR)and obtained different types of soil with frequencies of 500 MHz,250 MHz,and 100 MHz antennas.The soil profile data were obtained for 500 MHz,250 MHz,and 100 MHz antennas,and the dielectric properties of each type of soil were analyzed.In the image processing procedure,wavelet analysis was first used to decompose the pre-processed radar signal and reconstruct the high-frequency information to obtain the reconstructed signal containing the stratification information.Secondly,the reconstructed signal is taken as an envelope to enhance the stratification information.The Hilbert transform is applied to the envelope signal to find the time-domain variation of the instantaneous frequency and determine the time-domain location of the stratification.Finally,the dielectric constant of each soil horizon is used to obtain the propagation velocity of the electromagnetic wave at the corresponding position to obtain the stratification position of each soil horizon.The research results show that the 500 MHz radar antenna can accurately delineate Ap/Ah,horizon and the absolute accuracy of the stratification is within 5 cm.The effect on the soil stratification below the tillage horizon is not apparent,and the absolute accuracy of the 250 MHz and 100 MHz radar antennas on the stratification is within 9 cm.The overwhelming majority of the overall calculation errors are kept to within 15%.Based on the three central frequency antennas,the soil horizon detection rate reaches 93.3%,which can achieve accurate stratification of soil profiles within 1 m.The experimental and image processing methods used are practical and feasible;however,the GPR will show a missed detection for soil horizons with only slight differences in dielectric properties.Overall,this study can quickly and accurately determine the information of each soil stratification,ultimately providing technical support for acquiring soil configuration information and developing smart agriculture.
文摘A comprehensive Ground Penetration Radar (GPR) investigations and hazard assessment for the rehabilitation and strengthening of Habib Sakakini’s Palace in Cairo is presented herein, which is considered one of the most significant architectural heritage sites in Egypt. The palace located on an ancient water pond at the eastern side of Egyptian gulf besiding Sultan Bebris Al-Bondoqdary mosque is a place also called “Prince Qraja al-Turkumany pond”. That pond had been filled down by Habib Sakakini at 1892 to construct his famous palace in 1897. The integrated geophysical survey of the palace allowed the identification of several targets of potential archaeological and geotechnical engineering interest buried in fill and silty clay in the depth range between 100 - 700 cm. the methodological development focused on Multi-Fold (MF) Ground Penetrating Radar (GPR) imaging and subsurface characterization based on integrated velocity and attenuation analysis. Eight hundred sqm of Ground penetration Radar (GPR) profiling have been conducted to monitor the subsurface conditions. 600 meters are made in the surrounding area of the Palace and 200 sqm at the basement. The aim is to monitor the soil conditions beneath and around the Palace and to identify potential geological discontinuities, or the presence of faults and cavities. A suitable single and dual antenna are used (500 - 100 MHZ) is used to penetrate the desired depth of 7 meters (ASTM D6432). The GPR is used also detect the water table. At the building basement the GPR is used to identify the foundation thickness and soil-basement interface. As well as the inspection of cracks in some supporting columns, piers and masonry walls. The GPR also was used to investigate the floors and ceilings conditions and structural mapping. The results were validated by the geotechnical and structural surveys. All these results together with the seismic hazard analysis will be used for the complete analysis of the palace in the framework of the rehabilitation and strengthening works foreseen in a second stage.
文摘探讨了多小波函数及其预处理方法对探地雷达(ground penetrating radar—GPR)图象去噪性能的影响,在Donoho D L和Johnstone I M提出的小波阈值去噪方法的基础上提出了一个改进的阈值函数,并对实际的GPR图象进行阈值化处理和对比分析,结果表明选取合适的预处理方法,采用DGHM和STT多小波对GPR图象去噪可获得比其他方法更好的效果。
