The study presented in this manuscript aimed to relate the sedimentary strata imaged by the ground penetrating radar(GPR)method through numerical modeling with the mapping of sedimentary strata acquired through geotec...The study presented in this manuscript aimed to relate the sedimentary strata imaged by the ground penetrating radar(GPR)method through numerical modeling with the mapping of sedimentary strata acquired through geotechnical surveys.The study aimed to expose how obtaining subsoil information through noninvasive/destructive electromagnetic waves is beneficial,as they are reliable and less costly than drilling holes beyond what is necessary to have a subsurface mapping.In this sense,physical-geological modeling was carried out.The information on the type of sediments,acquired through simple recognition surveys carried out in the city of Belém-PA,helped to create a model of a sedimentary package with its respective intrinsic physical properties.The result shows that the GPR recovered with good vertical and horizontal resolution at the beginning and end of the layers of the sedimentary package studied,proving to be very effective for locating geotechnical sounding points and safely reducing costs.展开更多
In this paper,a layer-constrained triangulated irregular network( LC-TIN) algorithm is proposed for three-dimensional( 3 D) modelling,and applied to construct a 3 D model for geological disease information based on gr...In this paper,a layer-constrained triangulated irregular network( LC-TIN) algorithm is proposed for three-dimensional( 3 D) modelling,and applied to construct a 3 D model for geological disease information based on ground penetrating radar( GPR) data. Compared with the traditional TIN algorithm,the LCTIN algorithm introduced a layer constraint to the discrete data points during the 3 D modelling process,and it can dynamically construct networks from layer to layer and implement 3 D modelling for arbitrary shapes with high precision. The experimental results validated this method,the proposed algorithm not only can maintain the rationality of triangulation network,but also can obtain a good generation speed. In addition,the algorithm is also introduced to our self-developed 3 D visualization platform,which utilized GPR data to model geological diseases. Therefore the feasibility of the algorithm is verified in the practical application.展开更多
Ground penetrating radar (GPR) and the portable seismic property analyzer (PSPA) have been extensively used in the past two decades for monitoring, quantifying, and mapping the deterioration of bridge decks. Using PSP...Ground penetrating radar (GPR) and the portable seismic property analyzer (PSPA) have been extensively used in the past two decades for monitoring, quantifying, and mapping the deterioration of bridge decks. Using PSPA and GPR ensures regular monitoring of bridge conditions, leads to the early detection of deterioration. This research is to address the condition of August A. Busch bridge deck owned by the Missouri Department of Conservation. Visual inspection, GPR, and PSPA data were acquired on the bridge deck. Over 90% of the bridge deck was in fair to good condition with an average compressive strength of over 2500 psi. GPR data showed no indication of significant deterioration. The overall bridge deck was determined to be in fair to good condition.展开更多
Kom Ombo temple is one of temples which were belted over high plateau close to the River Nile, near to Aswan in Egypt in the Greek-Roman period. The expected archaeological remains in the selected area are the hidden ...Kom Ombo temple is one of temples which were belted over high plateau close to the River Nile, near to Aswan in Egypt in the Greek-Roman period. The expected archaeological remains in the selected area are the hidden tunnels of the mummified crocodiles. The aim of the present work is to detect any of these tunnels by the application of the (GPR) and (SP) methods. The interpretation of the 10 GPR profiles revealed some locations of possible hidden tunnels. These locations show different contrasts and high amplitudes of the reflected signals, compared to the enclosing soil;also the scattering of the signals is higher than the bed layer in these locations, which may reveal the possible buried mummified crocodile tunnels in the study area. The depths of the possible targets range from 2.0 m to 2.5 m. The SP electric map shows that the study area possesses a range of about 135 mV of the potential differences between the measured stations. The positive response of the SP data is mainly concentrated at the central part of the study area. The relatively weak, negative SP anomalies may be related to moisture in the soil. The positive SP anomalies on the SP electric map display possible significant correlation between them and the inferred tunnel locations from the GPR data. The calculated depths from the SP profiles show significant agreement with that estimated from GPR data depths, which indicate that the SP electric method can be used as a successful tool in detecting buried archaeological remains in support of GPR.展开更多
China’s Chang’e 3(CE-3)and Chang’e 4(CE-3)missions made historic progress by sending rovers equipped with scientific instruments to the new sites on the nearside and farside of the moon,respectively.The same lunar ...China’s Chang’e 3(CE-3)and Chang’e 4(CE-3)missions made historic progress by sending rovers equipped with scientific instruments to the new sites on the nearside and farside of the moon,respectively.The same lunar penetrating radar(LPR)that uses pulses of electromagnetic energy to reveal the underground structure and properties of the lunar soil"regolith"covering most of the lunar surface was carried by the rovers.It provided for the first time the opportunity for in situ LPR measurements of the subsurface substrate in two geologically different places on the moon.At present,the Yutu rover of CE-3 mission traveled along a path of total length of about 114 m while the Yutu 2 rover of CE-4 has traversed over 1000 m and keeps going.This paper summarizes the fruitful results so far obtained by LPRs,including the physical properties and layered structure of the lunar regolith and shallow crust beneath the two landing sites.The regolith layer thickness at the CE-3 site is thinner than that at the CE-4 site due to its relatively young age.The penetration depth below CE-4 site is about 2.85 times(in terms of the forward and return path delay)deeper than CE-3 as indicated by their different loss tangent values(0.0039±0.0002 vs.0.013),which is probably due to the differences in abundance of ilmenite and rocks in the regolith.Other physical parameters including dielectric permittivity vs.depth profile,bulk density and electrical conductivity have been estimated using various methods.Thanks to the low signal loss,CE-4 LPR is able to present clear cross section views of two buried craters and the paleo-surface(ancient surface)of the landing site.The multiple stratums observed by the low frequency channel of LPRs indicate several episodes of lava eruptions occurred in the late stage of the formation of the nearside Imbrium Basin and the farside Von Kármán crater despite the asymmetric volcanisms distributed on two sides.The discoveries revealed by LPRs advance our knowledges on the formation process and properties of lunar regolith,thickness of ejecta deposits caused by lunar impact events,the evolution of the nearside and farside volcanisms,etc.They also demonstrate the feasibility of applying ground penetrating radar for non-terrestrial explorations such as China’s first Martian mission,Tianwen 1 mission.展开更多
Ground penetrating radar (GPR) surveys have being applied to investigate very near-surface stratification of sedimentary units in coastal plains and to define their depositional conditions. This paper presents, howeve...Ground penetrating radar (GPR) surveys have being applied to investigate very near-surface stratification of sedimentary units in coastal plains and to define their depositional conditions. This paper presents, however, low-frequency GPR survey to investigate fault-related depositional systems at greater depths. The Quinta-Cassino area in the Rio Grande do Sul Coastal Plain (RGSCP, Brazil) shows a wide strandplain that is made off by very long, continuous, and linear geomorphic features (beach ridges). This strandplain extends for ~70 km southward. The beach ridges show low-angle truncations against the Quinta escarpment, and also truncations in the strandplain. The traditional approach points that RGSCP was developed by juxtaposition of four lagoons/barrier systems as consequence of sea level changes;previous model assumes that no deformational episode occurred in RGSCP. The geophysical and geological surveys carried out in this area showed the existence of listric fault controlling the beach ridges in the escarpments and hanging-wall blocks. The radargrams could distinguish Pleistocene basement unit anticlockwise rotation, thickening of beach ridges radarfacies close to listric normal faults, and horst structures. These deformational features indicate that the extensional zone of a large-scale gravity-driven structure controlled the mechanical subsidence, the Holocene sedimentation and its stratigraphic and geomorphic features in the Quinta-Cassino area to build up an asymmetric delta. The results point to a new approach in dealing with RGSCP Holocene evolution.展开更多
Soil water content(SWC)is a crucial parameter in ecology,agriculture,hydrology,and engineering studies.Research on non-invasive monitoring of SWC has been a long-lasting topic in these fields.Ground penetrating radar(...Soil water content(SWC)is a crucial parameter in ecology,agriculture,hydrology,and engineering studies.Research on non-invasive monitoring of SWC has been a long-lasting topic in these fields.Ground penetrating radar(GPR),a non-destructive geophysical technique,has the advantages of high resolution,deep detection depth,and high efficiency in SWC measurements at medium scale.It has been successfully applied in field investigations.This paper summarizes the recent progress in developing GPR-based SWC measurement methods,including reflected wave,ground wave,surface reflection,borehole GPR,full waveform inversion,average envelope amplitude,and frequency shift methods.The principles,advantages,limitations,and applications of these methods are described in detail.A comprehensive technical framework,which comprises the seven methods,is proposed to understand their principles and applicability.Two key procedures,namely,data acquisition and data processing,are emphasized as crucial to method applications.The suitable methods that will satisfy diverse application demands and field conditions are recommended.Future development,potential applications,and advances in hardware and data processing techniques are also highlighted.展开更多
Ground Penetrating Radar(GPR) is an effective Non-Destructive Testing(NDT) technique for highway pavement surveys, which is able to acquire continuous pavement data compared with traditional core drilling method. In t...Ground Penetrating Radar(GPR) is an effective Non-Destructive Testing(NDT) technique for highway pavement surveys, which is able to acquire continuous pavement data compared with traditional core drilling method. In this study, we proposed an accurate and efficient method to estimate the thickness of each pavement layer using an air-coupled GPR system. For this work, the main difficulties are estimating each pavement layer's time delay and dielectric constant. We first give the basic signal model for pavement evaluation, and then present an Intrinsic Mode Functions(IMFs) product detector to determine each pavement layer's time delay. This method is based on Empirical Mode Decomposition(EMD), which is an adaptive signal decomposition procedure and proved to be suitable for suppressing noises in GPR signal. The dielectric constant was determined by metal reflection measurement. The laboratory and highway experiments illustrate that the proposed thickness estimation method yields reasonable result, thus meets the requirements of practical highway pavement survey with massive GPR data.展开更多
Garbage disposal has gradually become a key issue faced by environmental governance. What must be controlled in garbage disposal is the calculation of garbage volume. In this paper, according to the different landfill...Garbage disposal has gradually become a key issue faced by environmental governance. What must be controlled in garbage disposal is the calculation of garbage volume. In this paper, according to the different landfill conditions of a landfill in Anhui, including: covering soil, overlying HDPE membrane, overlying concrete, etc., combined with the physical differences between the rock mass and the garbage at the bottom of the landfill. The landfill covered with HDPE membrane adopts the transient electromagnetic method, the landfill covered with soil layer adopts the high-density electrical method, and the landfill covered with concrete adopts the combination of transient electromagnetic method and ground penetrating radar. At the same time, Combines the borehole data to determine that the resistivity value of the interface between the garbage and the soil was 29 Ω⋅m. Finally, the sections of multiple survey lines can accomplish three-dimensional stereoscopic that calculate the distribution and actual volume of the garbage. This method is used to determine the final investment cost.展开更多
This paper presents an overview on the wide-ranging track structure studies at the Tampere University of Technology(TUT),Finland dealing with the key aspects of track geotechnics related to high-speed passenger traffi...This paper presents an overview on the wide-ranging track structure studies at the Tampere University of Technology(TUT),Finland dealing with the key aspects of track geotechnics related to high-speed passenger traffic on ballasted tracks.Special attention is paid to ballast and sub-ballast,while also considering frost action,embankment stability,track stiffness,track geometry and transition zones.As a result,this paper states that understanding the ballast degradation mechanism and its consequences and assessment of its condition occupy an important role in the construction and maintenance of a smooth high-speed rail line.The choices related to building the sub-ballast also have a dramatic impact on later track deformations and maintenance needs.In cold climate,especially where seasonal frost occurs,understanding and taking into account the frost action mechanism is crucial.Especially in the maintenance and rehabilitation planning of existing tracks,high-class analyses of ground penetrating radar data and its integrated analysis with other data can yield considerable benefits.展开更多
Finite-Difference Time-Domain(FDTD)is the most popular time-domain approach in computational electromagnetics.Due to the Courant-Friedrich-Levy(CFL)condition and the perfect match layer(PML)boundary precision,FDTD can...Finite-Difference Time-Domain(FDTD)is the most popular time-domain approach in computational electromagnetics.Due to the Courant-Friedrich-Levy(CFL)condition and the perfect match layer(PML)boundary precision,FDTD cannot simulate soil medium whose surface is connected by multiple straight lines or curves(multi-scale)accurately and efficiently,which greatly limits the application of FDTD method to simulate buried objects in soils.Firstly,this study proposed the absorption boundary and adopted two typical perfect matching layers(UPML and CPML)to compare their absorption effects,and then using the three forms of improved Yee-FDTD algorithm,alternating-direction implicit(ADI-FDTD),unconditionally stable(US-FDTD)and hybrid implicit explicit finite time-domain(HIE-FDTD)to divide and contrast the boundary model effects.It showed that the HIE-FDTD was suitable for inversion of multi-scale structure object modeling,while ADI-FDTD and US-FDTD were ideal for single-boundary objects in both uniaxial perfectly matched layer(UMPL)and convolution perfectly matched layer(CPML)finite element space.After that,all the models were tested by computer performance for their simulated efficiency.When simulating single boundary objects,UPML-US-FDTD and ADI-FDTD could achieve the ideal results,and in the boundary inversion of multi-scale objects,HIE-FDTD modeling results and efficiency were the best.Test modeling speeds of CPML-HIE-FDTD were compared with three kinds of waveform sources,Ricker,Blackman-Harris and Gaussian.Finally,under the computer condition in which the CPU was i5-8250,the HIE-FDTD model still had better performance than the traditional Yee-FDTD forward modeling algorithm.For modeling multi-scale objects in farmland soils,the methods used CPML combined with the HIE-FDTD were the most efficient and accurate ways.This study can solve the problem that the traditional FDTD algorithm cannot construct non-mesh objects by utilizing the diversity characteristics of Yee cell elements.展开更多
文摘The study presented in this manuscript aimed to relate the sedimentary strata imaged by the ground penetrating radar(GPR)method through numerical modeling with the mapping of sedimentary strata acquired through geotechnical surveys.The study aimed to expose how obtaining subsoil information through noninvasive/destructive electromagnetic waves is beneficial,as they are reliable and less costly than drilling holes beyond what is necessary to have a subsurface mapping.In this sense,physical-geological modeling was carried out.The information on the type of sediments,acquired through simple recognition surveys carried out in the city of Belém-PA,helped to create a model of a sedimentary package with its respective intrinsic physical properties.The result shows that the GPR recovered with good vertical and horizontal resolution at the beginning and end of the layers of the sedimentary package studied,proving to be very effective for locating geotechnical sounding points and safely reducing costs.
基金Supported by the National Science Foundation of China(61302157)the National High Technology Research and Development Program of China(863 Program)(2012AA12A308)the Yue Qi Young Scholars Project of China University of Mining&Technology(Beijing)(800015Z1117)
文摘In this paper,a layer-constrained triangulated irregular network( LC-TIN) algorithm is proposed for three-dimensional( 3 D) modelling,and applied to construct a 3 D model for geological disease information based on ground penetrating radar( GPR) data. Compared with the traditional TIN algorithm,the LCTIN algorithm introduced a layer constraint to the discrete data points during the 3 D modelling process,and it can dynamically construct networks from layer to layer and implement 3 D modelling for arbitrary shapes with high precision. The experimental results validated this method,the proposed algorithm not only can maintain the rationality of triangulation network,but also can obtain a good generation speed. In addition,the algorithm is also introduced to our self-developed 3 D visualization platform,which utilized GPR data to model geological diseases. Therefore the feasibility of the algorithm is verified in the practical application.
文摘Ground penetrating radar (GPR) and the portable seismic property analyzer (PSPA) have been extensively used in the past two decades for monitoring, quantifying, and mapping the deterioration of bridge decks. Using PSPA and GPR ensures regular monitoring of bridge conditions, leads to the early detection of deterioration. This research is to address the condition of August A. Busch bridge deck owned by the Missouri Department of Conservation. Visual inspection, GPR, and PSPA data were acquired on the bridge deck. Over 90% of the bridge deck was in fair to good condition with an average compressive strength of over 2500 psi. GPR data showed no indication of significant deterioration. The overall bridge deck was determined to be in fair to good condition.
文摘Kom Ombo temple is one of temples which were belted over high plateau close to the River Nile, near to Aswan in Egypt in the Greek-Roman period. The expected archaeological remains in the selected area are the hidden tunnels of the mummified crocodiles. The aim of the present work is to detect any of these tunnels by the application of the (GPR) and (SP) methods. The interpretation of the 10 GPR profiles revealed some locations of possible hidden tunnels. These locations show different contrasts and high amplitudes of the reflected signals, compared to the enclosing soil;also the scattering of the signals is higher than the bed layer in these locations, which may reveal the possible buried mummified crocodile tunnels in the study area. The depths of the possible targets range from 2.0 m to 2.5 m. The SP electric map shows that the study area possesses a range of about 135 mV of the potential differences between the measured stations. The positive response of the SP data is mainly concentrated at the central part of the study area. The relatively weak, negative SP anomalies may be related to moisture in the soil. The positive SP anomalies on the SP electric map display possible significant correlation between them and the inferred tunnel locations from the GPR data. The calculated depths from the SP profiles show significant agreement with that estimated from GPR data depths, which indicate that the SP electric method can be used as a successful tool in detecting buried archaeological remains in support of GPR.
文摘China’s Chang’e 3(CE-3)and Chang’e 4(CE-3)missions made historic progress by sending rovers equipped with scientific instruments to the new sites on the nearside and farside of the moon,respectively.The same lunar penetrating radar(LPR)that uses pulses of electromagnetic energy to reveal the underground structure and properties of the lunar soil"regolith"covering most of the lunar surface was carried by the rovers.It provided for the first time the opportunity for in situ LPR measurements of the subsurface substrate in two geologically different places on the moon.At present,the Yutu rover of CE-3 mission traveled along a path of total length of about 114 m while the Yutu 2 rover of CE-4 has traversed over 1000 m and keeps going.This paper summarizes the fruitful results so far obtained by LPRs,including the physical properties and layered structure of the lunar regolith and shallow crust beneath the two landing sites.The regolith layer thickness at the CE-3 site is thinner than that at the CE-4 site due to its relatively young age.The penetration depth below CE-4 site is about 2.85 times(in terms of the forward and return path delay)deeper than CE-3 as indicated by their different loss tangent values(0.0039±0.0002 vs.0.013),which is probably due to the differences in abundance of ilmenite and rocks in the regolith.Other physical parameters including dielectric permittivity vs.depth profile,bulk density and electrical conductivity have been estimated using various methods.Thanks to the low signal loss,CE-4 LPR is able to present clear cross section views of two buried craters and the paleo-surface(ancient surface)of the landing site.The multiple stratums observed by the low frequency channel of LPRs indicate several episodes of lava eruptions occurred in the late stage of the formation of the nearside Imbrium Basin and the farside Von Kármán crater despite the asymmetric volcanisms distributed on two sides.The discoveries revealed by LPRs advance our knowledges on the formation process and properties of lunar regolith,thickness of ejecta deposits caused by lunar impact events,the evolution of the nearside and farside volcanisms,etc.They also demonstrate the feasibility of applying ground penetrating radar for non-terrestrial explorations such as China’s first Martian mission,Tianwen 1 mission.
文摘Ground penetrating radar (GPR) surveys have being applied to investigate very near-surface stratification of sedimentary units in coastal plains and to define their depositional conditions. This paper presents, however, low-frequency GPR survey to investigate fault-related depositional systems at greater depths. The Quinta-Cassino area in the Rio Grande do Sul Coastal Plain (RGSCP, Brazil) shows a wide strandplain that is made off by very long, continuous, and linear geomorphic features (beach ridges). This strandplain extends for ~70 km southward. The beach ridges show low-angle truncations against the Quinta escarpment, and also truncations in the strandplain. The traditional approach points that RGSCP was developed by juxtaposition of four lagoons/barrier systems as consequence of sea level changes;previous model assumes that no deformational episode occurred in RGSCP. The geophysical and geological surveys carried out in this area showed the existence of listric fault controlling the beach ridges in the escarpments and hanging-wall blocks. The radargrams could distinguish Pleistocene basement unit anticlockwise rotation, thickening of beach ridges radarfacies close to listric normal faults, and horst structures. These deformational features indicate that the extensional zone of a large-scale gravity-driven structure controlled the mechanical subsidence, the Holocene sedimentation and its stratigraphic and geomorphic features in the Quinta-Cassino area to build up an asymmetric delta. The results point to a new approach in dealing with RGSCP Holocene evolution.
基金supported by the National Natural Science Foundation of China(Grant No.41571404)on project of State Key Laboratory of Earth Surface Processes and Resource Ecology.
文摘Soil water content(SWC)is a crucial parameter in ecology,agriculture,hydrology,and engineering studies.Research on non-invasive monitoring of SWC has been a long-lasting topic in these fields.Ground penetrating radar(GPR),a non-destructive geophysical technique,has the advantages of high resolution,deep detection depth,and high efficiency in SWC measurements at medium scale.It has been successfully applied in field investigations.This paper summarizes the recent progress in developing GPR-based SWC measurement methods,including reflected wave,ground wave,surface reflection,borehole GPR,full waveform inversion,average envelope amplitude,and frequency shift methods.The principles,advantages,limitations,and applications of these methods are described in detail.A comprehensive technical framework,which comprises the seven methods,is proposed to understand their principles and applicability.Two key procedures,namely,data acquisition and data processing,are emphasized as crucial to method applications.The suitable methods that will satisfy diverse application demands and field conditions are recommended.Future development,potential applications,and advances in hardware and data processing techniques are also highlighted.
基金Supported by the 863 National High Technology Research and Development Program(No.2012AA121901)
文摘Ground Penetrating Radar(GPR) is an effective Non-Destructive Testing(NDT) technique for highway pavement surveys, which is able to acquire continuous pavement data compared with traditional core drilling method. In this study, we proposed an accurate and efficient method to estimate the thickness of each pavement layer using an air-coupled GPR system. For this work, the main difficulties are estimating each pavement layer's time delay and dielectric constant. We first give the basic signal model for pavement evaluation, and then present an Intrinsic Mode Functions(IMFs) product detector to determine each pavement layer's time delay. This method is based on Empirical Mode Decomposition(EMD), which is an adaptive signal decomposition procedure and proved to be suitable for suppressing noises in GPR signal. The dielectric constant was determined by metal reflection measurement. The laboratory and highway experiments illustrate that the proposed thickness estimation method yields reasonable result, thus meets the requirements of practical highway pavement survey with massive GPR data.
文摘Garbage disposal has gradually become a key issue faced by environmental governance. What must be controlled in garbage disposal is the calculation of garbage volume. In this paper, according to the different landfill conditions of a landfill in Anhui, including: covering soil, overlying HDPE membrane, overlying concrete, etc., combined with the physical differences between the rock mass and the garbage at the bottom of the landfill. The landfill covered with HDPE membrane adopts the transient electromagnetic method, the landfill covered with soil layer adopts the high-density electrical method, and the landfill covered with concrete adopts the combination of transient electromagnetic method and ground penetrating radar. At the same time, Combines the borehole data to determine that the resistivity value of the interface between the garbage and the soil was 29 Ω⋅m. Finally, the sections of multiple survey lines can accomplish three-dimensional stereoscopic that calculate the distribution and actual volume of the garbage. This method is used to determine the final investment cost.
基金Project supported by the Finnish Transport Agency
文摘This paper presents an overview on the wide-ranging track structure studies at the Tampere University of Technology(TUT),Finland dealing with the key aspects of track geotechnics related to high-speed passenger traffic on ballasted tracks.Special attention is paid to ballast and sub-ballast,while also considering frost action,embankment stability,track stiffness,track geometry and transition zones.As a result,this paper states that understanding the ballast degradation mechanism and its consequences and assessment of its condition occupy an important role in the construction and maintenance of a smooth high-speed rail line.The choices related to building the sub-ballast also have a dramatic impact on later track deformations and maintenance needs.In cold climate,especially where seasonal frost occurs,understanding and taking into account the frost action mechanism is crucial.Especially in the maintenance and rehabilitation planning of existing tracks,high-class analyses of ground penetrating radar data and its integrated analysis with other data can yield considerable benefits.
基金This work was financially supported by the State Key Research Program of China(Grant No.2016YFD0700101)the State Key Research Program of China(Grant No.2017YFD0700404)+1 种基金the Guangdong Provincial Department of Agriculture’s Specialized Program for Rural Area Rejuvenation(Grant No.2019KJ129)and the Guangdong Provincial Department of Agriculture’s Modern Agricultural Innovation Team Program for Animal Husbandry Robotics(Grant No.200-2018-XMZC-0001-107-0130).
文摘Finite-Difference Time-Domain(FDTD)is the most popular time-domain approach in computational electromagnetics.Due to the Courant-Friedrich-Levy(CFL)condition and the perfect match layer(PML)boundary precision,FDTD cannot simulate soil medium whose surface is connected by multiple straight lines or curves(multi-scale)accurately and efficiently,which greatly limits the application of FDTD method to simulate buried objects in soils.Firstly,this study proposed the absorption boundary and adopted two typical perfect matching layers(UPML and CPML)to compare their absorption effects,and then using the three forms of improved Yee-FDTD algorithm,alternating-direction implicit(ADI-FDTD),unconditionally stable(US-FDTD)and hybrid implicit explicit finite time-domain(HIE-FDTD)to divide and contrast the boundary model effects.It showed that the HIE-FDTD was suitable for inversion of multi-scale structure object modeling,while ADI-FDTD and US-FDTD were ideal for single-boundary objects in both uniaxial perfectly matched layer(UMPL)and convolution perfectly matched layer(CPML)finite element space.After that,all the models were tested by computer performance for their simulated efficiency.When simulating single boundary objects,UPML-US-FDTD and ADI-FDTD could achieve the ideal results,and in the boundary inversion of multi-scale objects,HIE-FDTD modeling results and efficiency were the best.Test modeling speeds of CPML-HIE-FDTD were compared with three kinds of waveform sources,Ricker,Blackman-Harris and Gaussian.Finally,under the computer condition in which the CPU was i5-8250,the HIE-FDTD model still had better performance than the traditional Yee-FDTD forward modeling algorithm.For modeling multi-scale objects in farmland soils,the methods used CPML combined with the HIE-FDTD were the most efficient and accurate ways.This study can solve the problem that the traditional FDTD algorithm cannot construct non-mesh objects by utilizing the diversity characteristics of Yee cell elements.
