Search

Ground Penetrating Radar (GPR) is a non-destructive geophysical technique utilizing electromagnetic pulses to detect subsurface material properties. The analysis of regions of interest (ROIs) in GPR images often entails the identification of hyperbolic reflection regions of underground targets through accurate segmentation, a crucial preprocessing step. Currently, this represents a research gap. In the hyperbolic reflection region, manual segmentation not only demands professional expertise but is also time-consuming and error-prone. Automatic segmentation can aid in accurately determining the location and depth of the reflection region, thereby enhancing data interpretation and analysis. This study presents a deep residual Convolutional Neural Network (Res-CNN) that integrates skip connections within an encoder-decoder framework for ROI-binarized hyperbolic segmentation. The proposed framework includes designed downsampling and upsampling modules that facilitate feature computation sharing between these two modules through skip connections within network blocks. In the evaluation of both simple and complex models, our method attained PSNR, SSIM, and FSIM values of 57.1894, 0.9933, and 0.9336, and 58.4759, 0.9958, and 0.9677, respectively. Compared to traditional segmentation methods, the proposed approach demonstrated clearer segmentation results, enabling intelligent and effective identification of the ROI region containing abnormal hyperbolic reflection waves in GPR images. [ABSTRACT FROM AUTHOR]

Data assimilation for unconfined seepage analysis has faced significant challenges due to hybrid causes, such as sparse measurements, heterogeneity of porous media, and computationally expensive forward models. To address these bottlenecks, this paper introduces a physics-informed neural network (PINN) model to resolve the data assimilation problem for seepage analysis of unsaturated earth–rockfill dams. This strategy offers a solution that decreases the reliance on numerical models and enables an accurate and efficient prediction of seepage parameters for complex models in the case of sparse observational data. For the first attempt in this study, the observed values are obtained by random sampling of numerical solutions, which are then contributed to the synchronous constraints in the loss function by informing both the seepage control equations and boundary conditions. To minimize the effects of sharp gradient shifts in seepage parameters within the research domain, a residual adaptive refinement (RAR) constraint is introduced to strategically allocate training points around positions with significant residuals in partial differential equations (PDEs), which could facilitate enhancing the prediction accuracy. The model's effectiveness and precision are evaluated by analyzing the proposed strategy against the numerical solutions. The results indicate that even with limited sparse data, the PINN model has great potential to predict seepage data and identify complex structures and anomalies inside the dam. By incorporating coupling constraints, the validity of our PINN model could lead to theoretically viable applications of hydrogeophysical inversion or multi-parameter seepage inversion. The results show that the proposed framework can predict the seepage parameters for the entire research domain with only a small amount of observation data. Furthermore, with a small amount of observation data, PINNs are able to obtain more accurate results than purely data-driven DNNs. [ABSTRACT FROM AUTHOR]

For better stability, safety and water resource management in a dam, it is important to evaluate the amount of seepage from the dam body. This research is focused on machine learning approach to predict the amount of seepage from Pakistan’s Earth and rock fill Tarbela Dam during 2003 to 2015. The data of temperature, rainfall, water inflow, sediment inflow, reservoir level collected during 2003 to 2015 served as input while the seepage from dam during this period was the output. Artificial Neural Network (ANN), Random Forest (RF), Support Vector Machine (SVM), and CatBoost (CB), have been used to model the input-output relationship. The algorithms used to predict the dam seepage reported a high R2 scores between actual and predicted values of average seepage, suggesting their reliability in predicting the seepage in the Tarbela Dam. Moreover, the CatBoost algorithm outperformed, by achieving an R2 score of 0.978 in training, 0.805 in validation, and 0.773 in testing phase. Similarly, RMSE was 0.025 in training, 0.076 in validation, and 0.111 in testing phase. Furthermore, to understand the sensitivity of each parameter on the output (average seepage), Shapley Additive Explanations (SHAP), a model explanation algorithm, was used to understand the affect of each parameter on the output. A comparison of SHAP used for all the machine learning models is also presented. According to SHAP summary plots, reservoir level was reported as the most significant parameter, affecting the average seepage in Tarbela Dam. Moreover, a direct relationship was observed between reservoir level and average seepage. It was concluded that the machine learning models are reliable in predicting and understanding the dam seepage in the Tarbela Dam. These Machine Learning models address the limitations of humans in data collecting and analysis which is highly prone to errors, hence arriving at misleading information that can lead to dam failure.

As a fundamental part of ground penetrating radar (GPR) data processing, reverse time migration (RTM) can correctly position reflection waves and focusing diffraction waves on the proper spatial position. Least-squares reverse-time migration (LSRTM) is widely used in the seismic field for its ability to suppress artifacts and generate high-resolution images in comparison to conventional RTM. However, in the particular case of GPR detection methods, LSRTM is extremely susceptible to aliasing artifacts caused by under-sampling. In pursuit of enhanced precision in underground structure characterization, this paper presents the development of a new LSRTM based on modified total variation (MTV) regularization to improve imaging resolution. Initially, the objective function of LSRTM is derived by combining the Born approximation in 2-D transversal magnetic mode. Next, the adjoint equations and their gradients are solved using the Lagrange multiplier method. The objective function is then constrained by MTV regularization to ensure the precision and convergence of the LSRTM, which delivers a refined edge with reconstruction details. In the numerical experiments, in comparison to the conventional LSRTM method, the LSRTM-MTV algorithm demonstrated a 30.4% increase in computational speed and a 21.1% reduction in mean squared error (MSE). The outperformance of the proposed method is verified in detail through the image resolution and amplitude preservation in the test of synthetic data and laboratory data. Future research efforts will center on applying the proposed method to models featuring dispersive or anisotropic media that closely mimic real-world conditions and extending the application to various imaging techniques involving objective function minimization. [ABSTRACT FROM AUTHOR]

Seepage analysis is the central issue in solving the unconfined seepage problems from complex earth‐rock dams. The precision and efficiency of seepage simulation are the major technical barriers to exploring the heterogeneities of hydrological state variables in dams. This paper addresses the heterogeneity of the saturated‐unsaturated zones of dams and develops the smoothed finite element method (S‐FEM) capable of solving the seepage problem in domains with arbitrary geometry and continuously varied permeability. The gradient smoothing technique, in which the area integrals are transformed into the line integrals around edges of smoothing cells, is used to obtain the element matrices, which is used to simplify the solution of seepage problem. To be specific, we employ Van Genuchten method to modify the hydraulic conductivity of unsaturated zones medium and Signorini condition to convert the outflow boundary into the water head boundaries. The precision of the proposed method together with the importance of considering the unsaturated zones are first deliberated by the verification tests. Then, to quantify the impacts on the hydrological state variables, we focus on the heterogeneity of the entire zones through exploratory tests. The results show that the presence of leakage anomalies and auxiliary dam structures has a significant effect on the redistribution of hydrological state variables, which convinces that the findings may lay the foundation for the development of multi‐parameter seepage inversion and the derivation of coupled hydro‐geophysical inversion using the most appropriate starting model. [ABSTRACT FROM AUTHOR]

The essentiality of diffraction wavefront in structural mapping has been perceived, mainly because of its imaging potential for small structural variations and microscopic discontinuities. However, full extraction of ground penetrating radar (GPR) diffractions cannot be done perfectly even in theory as the faint components are typically masked by the more dominant reflections, resulting in a resolution loss of minutiae features. To address these issues, this letter develops a semblance-based coherence analysis framework for GPR diffraction wavefront reconstruction and full waveform inversion (FWI). Specifically, adaptive amplitude restoration is initially performed using the mean instant amplitude for analytic attenuation function derivation, which is followed by modified multitrace coherence stacking. Diffraction reconstruction is then executed with corrective factors determined by EM semblance attributes, while being assessed by quantified indices for reconstruction performance. The proposed framework is tested for validity and practicality with a synthetic model and laboratory experiment by comparing the conventional FWI and the diffraction driven FWI. The results present a high-resolution image of diffraction occurrences even with a homogeneous initial model, and more intuitively show a higher-resolution image even with a suboptimal initial model by diffraction-driven permittivity analysis. These encouraging results could potentially enable a hybrid inversion which benefits from reconstructing both GPR reflections and diffractions. [ABSTRACT FROM AUTHOR]

Estimating the quantity of seepage through the foundation and body of a dam using proper health and safety monitoring is critical to the effective management of disaster risk in a reservoir downstream of the dam. In this study, a deep learning model was constructed to predict the extent of seepage through Pakistan's Tarbela dam, the world's second largest clay and rock dam. The dataset included hydro-climatological, geophysical, and engineering characteristics for peak-to-peak water inflows into the dam from 2014 to 2020. In addition, the data are time series, recurring neural networks (RNN), and long short-term memory (LSTM) as time series algorithms. The RNN–LSTM model has an average mean square error of 0.12, and a model performance of 0.9451, with minimal losses and high accuracy, resulting in the best-predicted dam seepage result. Damage was projected using a deep learning system that addressed the limitations of the model, the difficulties of calculating human activity schedules, and the need for a different set of input data to make good predictions. [ABSTRACT FROM AUTHOR]

This study aims to develop a novel wavelet neural-network (WNN) model for solving electrical resistivity imaging (ERI) inversion with massive amounts of measured data in control and measurement fields. In the proposed method, we design a mixed multilayer WNN (MMWNN) which uses Morlet and Mexican wavelons as different activation functions in a cascaded hidden layer structure. Meanwhile, a hybrid STGWO-GD learning approach is used to improve the learning ability of the MMWNN, which is a combination of the self-tuning grey wolf optimizer (STGWO) and the gradient descent (GD) algorithm adopting the advantages of each other. Moreover, updating formulas of the GD algorithm are derived, and a Gaussian updating operator with weighted hierarchical hunting, a chaotic oscillation equation, and a nonlinear modulation coefficient are introduced to improve the hierarchical hunting and the control parameter adjustment of the modified STGWO. Five examples are used with the aim of assessing the availability and feasibility of the proposed inversion method. The inversion results are promising and show that the introduced method is superior to other competitors in terms of inversion accuracy and computational efficiency. Furthermore, the effectiveness of the proposed method is demonstrated over a classical benchmark successfully. [ABSTRACT FROM AUTHOR]

Subsurface mapping can be implemented by the migration imaging of ground penetrating radar (GPR) data. Finite difference (FD)-based migration can effectively achieve fast imaging in media with lateral heterogeneity. However, this migration inherently suffers from insufficient accuracy in determining the true morphology of any steep-obliquity interface. In this letter, we propose a new algorithm based on crossed differential operator (CDO) and reverse-time migration (RTM) strategy to refocus the reflections back into its precise location. By employing the implicit Crank–Nicolson scheme without any operator approximations, wave field can be extrapolated nearly along the characteristic directions of wave propagation based on the nonapproximate equations. We determined the potential of the algorithm using synthetic and field data. For synthetic data, we consider a multigeometry scenario, including different angle interfaces, different-sized spheres, and lateral heterogeneity. In particular, we focus on extrapolation accuracy, amplitude–frequency focusing ability, entropy analysis, and imaging errors. In the field data test, we found that the finer morphology of the steeper-obliquity interfaces can be determined, the interpretability and amount of discernible details can also be improved. This improvement makes CDO-based migration a very promising tool for the finer imaging of shallow subsurface. [ABSTRACT FROM AUTHOR]

Efficient and accurate simulation of ground penetrating radar (GPR) in the open region helps immensely in both grasping the features of echoes and facilitating the interpretation of real GPR data. Due to the limitation of the computer model, however, the strong artificial boundary reflections, especially the low-frequency propagating waves encountered at the late stage of simulation greatly affect the simulation accuracy of GPR. This paper presents an innovative optimized unsplit-field convolutional perfectly matched layer (UCPML) based on rotated staggered grid (RSG) scheme to truncate the finite-difference time-domain (FDTD) lattices. Rather than obey the sharp variation based on an ${m}$ th-order polynomial, the optimized approach employs a novel optimized term and an adjustment factor to seek a gentle variation on optimal constitutive coefficients. This guarantees that the determination of optimal constitutive coefficients can be less influenced by the order of polynomial and especially, to improve the absorptive performance on low-frequency propagating waves. The calculating efficiency and accuracy of the RSG-FDTD scheme, as well as the absorbing performance of the optimized UCPML, are verified by two numerical examples. In particular, the analysis of the amplitude-frequency features of low-frequency clutters at steady state of the electromagnetic (EM) field and the corresponding global reflection error in the time–frequency domain is also presented. [ABSTRACT FROM AUTHOR]

An integrated geophysical investigation was performed at S dam located at Dadu basin in China to assess the condition of the dam curtain. The key methodology of the integrated technique used was flow-field fitting method, which allowed identification of the hydraulic connections between the dam foundation and surface water sources (upstream and downstream), and location of the anomalous leakage outlets in the dam foundation. Limitations of the flow-field fitting method were complemented with resistivity logging to identify the internal erosion which had not yet developed into seepage pathways. The results of the flow-field fitting method and resistivity logging were consistent when compared with data provided by seismic tomography, borehole television, water injection test, and rock quality designation. [ABSTRACT FROM AUTHOR]

We introduce a novel implementation of a rotated staggered grid (RSG) finite difference (FD) method based on the unsplit convolutional perfectly matched layer (UCPML). The method is carried out using rotated staggered FD operators and a general complex frequency shifted stretching factor. It can offer a number of advantages over the RSG with split convolutional perfectly matched layer (SCPML), such as improving computational efficiency, overcoming the numerical dispersions resulting from the high-frequency propagation and strong contrasts in material, and specifically absorbing the evanescent modes at grazing incidence. A half-space numerical example with a lossy medium attached demonstrates that large-scale spurious reflections and evanescent modes at grazing incidence are absorbed efficiently by the UCPML in RSG scheme. Meanwhile, the comparison of the reflection error between the proposed algorithm and traditional approaches [using time-dependent reflection error (TDRE) and frequency-dependent reflection error (FDRE)] suggests that the proposed algorithm of RSG provides a preferable improvement in absorbing low-frequency coda wave of evanescent modes, particularly in the lossy medium that possesses inhomogeneous discontinuities. [ABSTRACT FROM PUBLISHER]

The chaos oscillation particle swarm optimization (COPSO) algorithm is prone to binge trapped in the local optima when dealing with certain complex models in ground-penetrating radar (GPR) data inversion, because it inherently suffers from premature convergence, high computational costs, and extremely slow convergence times, especially in the middle and later periods of iterative inversion. Considering that the bilateral connections between different particle positions can improve both the algorithmic searching efficiency and the convergence performance, we first develop a fast single-trace-based approach to construct an initial model for 2-D PSO inversion and then propose a TV-regularization-based improved PSO (TVIPSO) algorithm that employs total variation (TV) regularization as a constraint technique to adaptively update the positions of particles. B by adding the new velocity variations and optimal step size matrices, the search range of the random particles in the solution space can be significantly reduced, meaning blindness in the search process can be avoided. By introducing constraint-oriented regularization to allow the optimization search to move out of the inaccurate region, the premature convergence and blurring problems can be mitigated to further guarantee the inversion accuracy and efficiency. We report on three inversion experiments involving multilayered, fluctuated terrain models and a typical complicated inner-interface model to demonstrate the performance of the proposed algorithm. The results of the fluctuated terrain model show that compared with the COPSO algorithm, the fitness error (MAE) of the TVIPSO algorithm is reduced from 2.3715 to 1.0921, while for the complicated inner-interface model the fitness error (MARE) of the TVIPSO algorithm is reduced from 1.9539 to 1.5674. [ABSTRACT FROM AUTHOR]