Determination of multilayer soil parameters using Simulated Annealing algorithm
The paper presents a new method based on the Simulated Annealing Algorithm for the estimation of the parameters of multilayer soils. First the infinite integral of resistivity equation is approximated by use of complex image method. Then, Simulated Annealing Method is used to estimate the soil resistivity with the least error.
Parameter estimation of horizontal multilayer earth by complex image method
How to correctly obtain the structure parameters of earth is the precondition of grounding system design. Based on the local searching algorithm, an efficient method for estimating earth structure with any number of layers from the test data by Wenner configuration method is presented. In each iterative procedure, earth apparent resistivities are calculated by complex image method. A new method to directly obtain the partial derivatives of the earth apparent resistivities with respect to the earth parameters from the results of the complex image method is presented, which is more efficient than the conventional methods that arduously calculate the improper integrals with a Bessel's function and an infinite limit of integration. At the same time, a method to determine suitable initial earth parameters for the local searching algorithm is also proposed. The validity of the method is shown by several applications.
Archaeological Survey of Dangoyama Kofun by 3DGPR
We presented a three-dimensional ground penetrating radar (3DGPR) system for the archeological investigation in a harsh site, where we have to deal with many trees and/or great topographic variations. This system consists of a commercial GPR system with a laser positioning system. With this system, we can move the GPR antenna freely, on which a positioning detector can record accurate position for each GPR data trace. We applied our 3DGPR system at Dangoyama Kofun, a burial mound, in Fukushima prefecture, Japan. A lot of trees and bushes covered the survey area and cannot be cut. We acquired six datasets in two days in this area. We successfully imaged the boundary between the original surface of the burial mound and the covering soil clearly at about 0.3 m depth around the mound. Furthermore, we found some strong reflections appeared at 0.8 m depth in the top area of the mound, which can be from the artificial parts of the burial mound. The results would greatly aid the following archaeological excavation. In addition, we compared the non-interpolated depth slice with interpolated depth slice at same position and demonstrated that the interpolation process on threedimensional GPR data may introduce some artifacts which might mislead our interpretation.