Geophysical modelling performs to obtain subsurface structures in agreement with measured data. Freeware algorithms for geoelectrical data inversion have not been widely used in geophysical communities; however, different open-source modelling/inversion algorithms were developed in recent years. In this study, we review the structures and applications of openly Python-based inversion packages, such as pyGIMLi (Python Library for Inversion and Modelling in Geophysics), BERT (Boundless Electrical Resistivity Tomography), ResIPy (Resistivity and Induced Polarization with Python), pyres (Python wrapper for electrical resistivity modelling), and SimPEG (Simulation and Parameter Estimation in Geophysics). In addition, we examine the recovering ability of pyGIMLi, BERT, ResIPy, and SimPEG freeware through inversion of the same synthetic model forward responses. A versatile pyGIMLi freeware is highly suitable for various geophysical data inversion. The SimPEG framework is developed to allow the user to explore, experiment with, and iterate over multiple approaches to the inverse problem. In contrast, BERT, pyres, and ResIPy are exclusively designed for geoelectric data inversion. BERT and pyGIMLi codes can be easily modified for the intended applications. Both pyres and ResIPy use the same mesh designs and inversion algorithms, but pyres uses scripting language, while ResIPy uses a graphical user interface (GUI) that removes the need for text inputs. Our numerical modelling shows that all the tested inversion freeware could be effective for relatively larger targets. pyGIMLi and BERT could also obtain reasonable model resolutions and anomaly accuracies for small-sized subsurface structures. Based on the heterogeneous layered model and experimental target scenario results, the geoelectrical data inversion could be more effective in pyGIMLi, BERT, and SimPEG freeware packages. Moreover, this study can provide insight into implementing suitable inversion freeware for reproducible geophysical research, mainly for geoelectrical modelling.

We aim to develop a comprehensive tunnel lining detection method and clustering technique for semi-automatic rebar identification in order to investigate the ten tunnels along the South-link Line Railway of Taiwan (SLRT). We used the Ground Penetrating Radar (GPR) instrument with a 1000 MHz antenna frequency, which was placed on a versatile antenna holder that is flexible to the tunnel’s condition. We called it a Vehicle-mounted Ground Penetrating Radar (VMGPR) system. We detected the tunnel lining boundary according to the Fresnel Reflection Coefficient (FRC) in both A-scan and B-scan data, then estimated the thinning lining of the tunnels. By applying the Hilbert Transform (HT), we extracted the envelope to see the overview of the energy distribution in our data. Once we obtained the filtered radargram, we used it to estimate the Two-dimensional Forward Modeling (TDFM) simulation parameters. Specifically, we produced the TDFM model with different random noise (0–30%) for the rebar model. The rebar model and the field data were identified with the Hierarchical Agglomerative Clustering (HAC) in machine learning and evaluated using the Silhouette Index (SI). Taken together, these results suggest three boundaries of the tunnel lining i.e., the air–second lining boundary, the second–first lining boundary, and the first–wall rock boundary. Among the tunnels that we scanned, the Fangye 1 tunnel is the only one in category B, with the highest percentage of the thinning lining, i.e., 13.39%, whereas the other tunnels are in category A, with a percentage of the thinning lining of 0–1.71%. Based on the clustered radargram, the TDFM model for rebar identification is consistent with the field data, where k = 2 is the best choice to represent our data set. It is interesting to observe in the clustered radargram that the TDFM model can mimic the field data. The most striking result is that the TDFM model with 30% random noise seems to describe our data well, where the rebar response is rough due to the high noise level on the radargram.

This study aims to apply geophysical methods to determine the Specific Yield (Sy) and Groundwater Level (GWL) in an unconfined aquifer of the Pingtung Plain in South Taiwan. Sy is an important hydraulic parameter for assessing groundwater potential. Obtaining specific yield for a large area is impractical due to the limited coverage and the high cost of the pumping test, which limits the potential evaluation of regional groundwater. Therefore, we used time-lapse Electrical Resistivity Imaging (ERI) to determine the Sy and GWL. Seasonal variations were considered when measuring time-lapse resistivity for five different months in 2019. We calculated the Sy and GWL from inverted resistivity data using empirical formulas and the soil–water characteristic curve (SWCC). We first used Archie’s law to calculate the relative saturation change with depth for each ERI profile, and then we used the Van Genuchten (VG) and Brooks–Corey (BC) empirical equations to estimate Sy and GWL. Finally, we compared the obtained GWL to the existing observation well to verify the findings of our study. The results showed that the VG and BC are able to predict Sy and GWL; however, the BC result is less consistent with the observation well result. In the study area, the dry season GWL ranged from 24.5 m to 35.2 m for the VG results and from 25.7 m to 35.5 m for the BC results. The wet season GWL ranged from 26.5 m to 38.9 m for the VG and from 26.4 m to 38.2 m for the BC results. The spatial distribution of the GWL shows a high gradient of GWL in the northeastern region, induced by significant proximal fan recharge. The determined spatial distribution of Sy varies from 0.15 to 0.21 for the VG and 0.14 to 0.20 for the BC results, indicating the study area has significant potential for groundwater resources. Therefore, nondestructive resistivity imaging can be used to aid in the determination of hydraulic parameters.

In this study, we attempt to establish an alternative method for estimating the groundwater levels and the specific yields of an unconfined aquifer for the evaluation of potential groundwater reservoirs. We first converted the inverted resistivity into the normalized water content. Then, we inverted the parameters of the Brooks-Corey model from the vertical profiles of the water content by assuming that the suction head was in proportion to the elevation regarding a predefined base level. Lastly, we estimated the groundwater level, the theoretical specific yield, and the specific yield capacity from the Brooks-Corey parameters at every survey site in the study area. The contour maps of the time-lapse groundwater levels show that the groundwater flows downstream, with a higher hydraulic gradient near the river channel than in the area away from the main channel. We conclude that the estimated maximum specific yield capacities are consistent with that derived from the pumping tests in the nearby observation well. Additionally, the specific yield capacities are only three quarters to two thirds of the theoretical specific yields derived from the difference between the residual and saturated water contents in the Brooks-Corey model. We conclude that the distribution pattern of the specific yields had been subjected to the distribution of natural river sediments in the Minzu Basin, since the modern channel was artificially modified. Although we had to make some simple assumptions for the estimations, the results show that the surface resistivity surveys provide reasonable estimations of the hydraulic parameters for a preliminary assessment in an area with few available wells.

No abstract