Next Article in Journal
AR-AI Tools as a Response to High Employee Turnover and Shortages in Manufacturing during Regular, Pandemic, and War Times
Next Article in Special Issue
Minimum Safety Factor Evaluation of Slopes Using Hybrid Chaotic Sand Cat and Pattern Search Approach
Previous Article in Journal
Spatial Vegetation Patch Patterns and Their Relation to Environmental Factors in the Alpine Grasslands of the Qilian Mountains
Previous Article in Special Issue
Combining Numerical Simulations, Artificial Intelligence and Intelligent Sampling Algorithms to Build Surrogate Models and Calculate the Probability of Failure of Urban Tunnels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prediction of Permeability Coefficient k in Sandy Soils Using ANN

by
Grzegorz Wrzesiński
* and
Anna Markiewicz
Institute of Civil Engineering, Warsaw University of Life Sciences, Nowoursynowska 159 St., 02-776 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(11), 6736; https://0-doi-org.brum.beds.ac.uk/10.3390/su14116736
Submission received: 12 April 2022 / Revised: 16 May 2022 / Accepted: 29 May 2022 / Published: 31 May 2022

Abstract

:
The paper presents a method of application of an ANN (Artificial Neural Network) to predict the permeability coefficient k in sandy soils: FSa, MSa, CSa. To develop an ANN the results of permeability coefficients from pumping and consolidation tests were applied. The proposed ANN with an architecture 6-8-1 predicts the value of permeability coefficient k based on the following parameters: soil type, relative density ID, void ratio e and effective soil diameter d10. The mean relative error and single maximum value of the relative error for the proposed ANN are following: Mean RE = ±4%, Max RE = 7.59%. The use of the ANN to predict the soil permeability coefficient allows the reduction of the costs and time needed to conduct laboratory or field tests to determine this parameter.

1. Introduction

Due to the intensive development of the construction industry, areas with difficult subsoil conditions or high groundwater levels are more and more often built on. Therefore, before starting the preparation of an investment project, it is necessary to thoroughly identify the subsoil, i.e., to determine, apart from subsoil strength parameters, the permeability coefficient k. The parameter k describes the soil’s hydraulic permeability, i.e., the filtration ability of water through the subsoil. The permeability coefficient k in the subsoil depends on the following properties: type of soil, grain size and shape, porosity, soil structure and water viscosity [1]. Typical values of the parameter k in mineral soils are in the range 1 × 10−2 ÷ 1 × 10−10 m·s−1.
The permeability coefficient k is determined using many methods, ranging from estimations based on uncomplicated calculations to complex laboratory and field tests. Choosing the method of determining the parameter k depends mainly on the type of soil. Eurocode 7 distinguishes four methods of the evaluation of the parameter k: laboratory tests, field tests, estimation based on oedometer test and empirical correlations. Selection of the method for determining the parameter k should be preceded by a thorough analysis of the properties and homogeneity of the tested soil, as well as of the factors that may affect the tested parameter [2,3].
Permeability coefficient k in soils is often estimated based on empirical correlations, in which the content of grain size (most often the effective soil diameter d10), porosity and specific surface area are used [4,5]. However, existing empirical formulas ignore the impact of soil structure, permeability anisotropy and the soil grain shape. Research indicates that the parameter k of the same material may differs significantly if it is estimated based on different empirical formulas [6,7]. Therefore, it is crucial to determine the exact empirical formulas and the correctness of the values of the parameters used, which may affect the calculation results. It is also important that the empirical formulas depend on a greater number of parameters of the analysed soil and the exact determination of the parameters used in empirical formulas.
The most accurate way to evaluate the permeability coefficient k are field tests. In field tests the heterogeneity of the subsoil structure and hydraulic permeability anisotropy are best reflected. Obviously, the accuracy of the determined parameter depends on the correctness of the performed tests. The most common field test is pumping test. The main goal of the pumping test is to obtain, by pumping water from the well (piezometer), a hydrodynamic reaction of the subsoil, which enables the identification of the filtration parameters of the soil and the conditions of its supply [8,9,10]. The pumping test is most often used to estimate the parameter k in well-permeable soils, mainly sands. In soils with less permeability (k < 10−6 m·s−1), the BAT probe test is used. This test involves the connection of a piezometer to the probe measuring part with a glass water container in which pressure changes are registered. Parameter k is determined on the basis of pressure changes in a function of time. Depending on the degree of filling the pores with water in the soil, the test may be performed under the conditions of water supply (inflow) or outflow from the probe tip. In the laboratory, parameter k is determined using constant or variable gradient methods. Constant-gradient methods, including the Rowe chamber and Trautwein system and ZW-K2 apparatus, are most often used to evaluate the parameter k in well-permeable soils [11]. Variable-gradient methods, including the flow-pump method, Kamienski tube method and modified oedometer with a burette, are used to evaluate the parameter k in low-permeable soils. Of the above-mentioned variable-gradient methods, the flow-pump method is the most commonly used. In this method, the differences in pressure at the top and bottom of the soil sample are measured after establishing a constant water flow velocity in the test sample. The test is continued until the pressure stabilizes in the soil sample.
ANNs (Artificial Neural Networks) are a tool used in science for modelling complex phenomena [12,13,14,15]. In recent years, they have gained more and more popularity and are widely used for statistical analysis [16]. In recent years they have been successfully applied in civil engineering [17,18,19,20]. Waszczyszyn [21,22] used neural networks in the analysis of hybrid computational systems, regression problems, development of Recurent Cascade Neural Networks (RCNNs), Binarized Neural Networks (BNNs) and Principal Component Analysis. Rafiq et al. [23] applied three types of neural networks in engineering: multi-layer perceptron (MLP), normalized RBF (NRBF) and radial basis network (RBF). Research on the use of ANNs in construction, in particular in relation to the subsoil, was presented inter alia by: Shahin et al. [24,25], Dihoru et al. [26], Pichler et al. [27], Sulewska [28], Tian et al. [29], Zhou et al. [30], Ellis et al. [31], Sidarta and Ghaboussi [32], Penumadu and Zhao [33], Basheer [34], Najjar and Huang [35], Fu et al. [36], Lee et al. [37], Das and Basudhar [38], Byeon et al. [39], Wrzesinski [40], Wrzesinski et al. [41], Lin et al. [42,43], Liu et al. [44]. So far, many applications of ANNs in construction have been reported in the literature; however, there are no reports on the using of ANNs to evaluate the parameter k in non-cohesive soils.
The paper presents the using ANN (Artificial Neural Network) to predict the permeability coefficient k in sandy soils. In order to propose an ANN, a series of laboratory and field tests to determine parameter k were carried out and the basic index properties of tested soils were determined.

2. Materials and Methods

Research was performed on sands characterized by different grain size. The study started with the selection of 50 test sites in the field. The criterion for selecting the test sites was the subsoil with homogeneous permeable sandy soils and the possibility of performing field tests and collecting samples for laboratory tests. The conducted research to determine permeability coefficients was divided into field and laboratory tests. The field studies included pumping tests, while the laboratory studies included consolidometer tests and tests to determine basic index properties. A total of 50 field tests and 120 laboratory tests were performed, in which the permeability coefficients k were determined.
Pumping tests to determine permeability coefficients k were performed in accordance with the standard method [45,46,47]. In carried research 1 well and 2 piezometers were installed for each test. The pumping test started with the installation of a well in the subsoil. The diameters of all installed wells were 400 mm and they were installed at a depth of 1.10 ÷ 2.00 m. The installation depth depended on the test site. Then, 2 piezometers were installed in the subsoil. The first piezometer was installed at a distance of 2.0 m, the second at a distance of 4 or 5 m from the edge of the well. PVC wells and piezometers were used. The pumping tests were carried out by pumping water and simultaneous measurements of changes in the water level inside installed piezometers. The tests were repeated two or three times, at the test sites in particular, to verify the results. Differences in the permeability coefficients k did not exceed 5% (at the same test site). The water level in the piezometers was measured 28–30 days after the end of the studies. Parameters k in pumping tests were determined in accordance with the following formula:
k = Q π z 2 2 z 1 2 l n x 2 x 1    
where Q—water flow, z2—water level in piezometer 2, z1—water level in piezometer 1, x2—distance between well and piezometer 2 and x1—distance between well and piezometer 1.
A typical scheme of the performed pumping tests is presented in Figure 1.
Additionally, cone penetrometer tests (CPT) and electrical resistivity measurements were performed before pumping tests. Laboratory tests were performed on soil samples collected from the test sites. In laboratory the following tests were performed: tests on grain size distribution of soil, on the permeability coefficient in consolidometer and tests using a scanning electron microscope.
Grain size distribution of tested soils were carried out in accordance with EN ISO 14688-1: 2002 and EN ISO 14688-2: 2004 [48,49]. Based on these tests, the effective soil diameter d10 was evaluated. The results of the grain size distribution tests indicate that the analysed soils, according to EN ISO 14688-2: 2006 and EN ISO 14688-2: 2006-Ap2: 2012, are FSa-Fine Sand, MSa-Middle Sand and CSa-Coarse Sand. According to the EN ISO 14688-2: 2006 standard, the tested sands were characterized by a gravel fraction of 0 ÷ 19%, sand fraction of 80 ÷ 98%, silt fraction of 0 ÷ 12% and clay fraction of 0 ÷ 2%.
Relative densities ID of tested soils were determined using cone penetrometer tests (CPT). Void ratios e were determined using electrical resistivity measurement. After connecting the electrodes to the power source and the gauge, several dozen measurements were made with a frequency of 12 s for the same place. The arithmetic mean of the measurement results was taken as the final result of the electrical resistance.
For the analysed soils, tests of the structure and chemical composition were carried out using a scanning electron microscope (XL Series, FEG Quanta 250 model) equipped with a chemical composition analysis system based on EDS X-ray energy dispersion.
Permeability coefficient tests to determine parameters k were performed in a laboratory with the use of consolidometer (Figure 2). The tests started with the compaction of soil in the Proctor apparatus. The diameters of compacted samples were 150 mm and their heights were 60 mm. Soil samples were compacted until the values of the relative density ID determined in the field studies using the cone penetrometer tests (CPT) were reached [50,51]. The determination of the permeability coefficient was based on the analysis of the consolidation process in the uniaxial state of strain. For this purpose, consolidometer tests were performed using the flow pump method for soil saturation. Tests were performed with a continuous inflow of water at constant gradients of 0.50. Based on the obtained characteristics, the permeability coefficients k were determined by the Taylor method. The values of the parameter k for particular soils obtained with the same gradients did not differ from each other by more than 5%.
The grain size distribution, relative density ID, void ratio e and effective soil diameter d10 of tested soils can be seen in Table 1.

3. Results

Performed research allowed the determination of the parameters k of tested soils using two methods: the pumping test and the consolidometer test. The parameters k, evaluated based on performed tests, are presented in Table 2.
On the basis of the conducted research, it can be stated that similar values of the permeability coefficient k were obtained in field and laboratory tests for soils from the same test site. The influence of the relative density ID on the permeability coefficients k is significant in the tested soils. Higher values of permeability coefficients k were in soils characterized by a lower relative density ID, higher void ratio e and lower value of effective soil diameter d10.
The tests carried out using the scanning electron microscope showed that the tested non-cohesive soils have different grain shapes, even if the sieve analysis shows a similar content of individual fractions and chemical composition. Sample SEM photos showing the structure and chemical composition spectra for selected sample of fine sand (FSa), characterized by similar relative density ID from two test sites (wells no. 1 and 6), are presented in Figure 3.
The conducted research shows that the shape of the soil grains and arrangement of grains has an impact on the parameter k. The performed tests confirmed that irregularly shaped soil grains hold more water in the micro-cavities in comparison with grains with more regular shapes. The differences in parameter k are significant for the same sand characterized by similar relative density ID, void ratio e and the effective soil diameter d10. For instance, the difference in parameter k is two-fold in FSa between wells no. 1 and 6. The impact of the analysed parameters on the parameter k is the highest in fine sands (FSa).

4. ANN (Artificial Neural Network) Analysis

4.1. Architecture of ANN

The class of ANNs used in the presented research is that of multilayer perceptrons (MLPs) with one hidden layer. The architecture of the ANNs N-H-M type is defined by: N—number of nodes in the input layer X1XN, H—number of nodes in the hidden layer, M = Y—number of nodes in the output layer.
ANNs were developed in accordance with the criterion of minimization of the error function. The criterion was the sum of squares of differences (SOS) evaluated using the equation:
E S O S = i = 1 P d i y i 2  
where P—number of cases of set P, di—known values of the tests, yi—predicted values using ANN.
Errors measures were separately calculated for subsets: training Tr, testing T and validation V. The assessment of the ANN’s predictive quality took into account the value of the error measures in the test subset T. The selected error measures for ANNs were determined using the following formulas:
-
Relative error for individual cases:
R E i = d i y i d i · 100 %    
-
Determination coefficient R2:
R 2 = 1 i = 1 P d i y i 2 i = 1 P d i d ¯ i 2  
-
Mean absolute error:
M A E = 1 i = 1 P d i y i P      
-
Root mean squared error:
R M S = 1 P · i = 1 P d i y i 2          
where di—measured value, yi—predicted value by ANN, d ¯ i —measured mean value in subset.

4.2. Data Sets, Training and Testing the ANN

In the performed ANN analysis, two data sets were used: set A and set B. Set A was obtained in field and laboratory tests and consisted of nA = 50 cases. Set A was described by five variables: X1 = soil type Є {FSa; MSa, CSa}; X2 = relative density ID Є {0.22 ÷ 0.92}; X3 = void ratio e Є {0.405 ÷ 0.728%}; and X4 = effective soil diameter d10 Є {0.04 ÷ 0.63} as well as Y = permeability coefficient k (Table 1). Variable X1 was treated as quality variable “one from N” type and required the use of many input neurons because three soil types were used. Set B consisted of data obtained in consolidometer tests and included nB = 120 cases. Set B was described by independent variables X1X4, and dependent variable Y = k. The development and training of the ANN was carried out with the use of set B. The selected ANN was used to predict permeability coefficient k on the basis of new data from field tests (set A).
Set B was randomly assigned to 60, 30 and 30 laboratory test results of the subsets: training Tr, testing T, and validation V, respectively. Subset testing T was used to periodically check the generalizability acquired by the ANN, while subset validation V was used for the final evaluation of the trained ANN.
Predictive quality of the neural regression model was estimated based on performed error analysis. Errors were determined in subset training Tr, testing T and validation V. ANNs were optimized in terms of the training method, the number of neurons in the hidden layer, the activation function of neurons in the hidden layer and the output layer. A Conjugate gradient (CG) and early stopping methods were used to train the network. The Conjugate gradient method is an algorithm for the numerical solution of particular systems of linear equations, namely those whose matrix is positive-definite, while the early stopping method is a form of regularization used to avoid overfitting when training a learner with an iterative method, such as gradient descent. Such methods update the learner so as to make it better fit the training data with each iteration. The training was finished after 28 training cycles. ANNs with the best predictive quality were evaluated on the basis of the lowest values of root mean squared errors RMS, the highest values of determination coefficient R2 and the lowest mean values of relative errors RE [52].
The performed neural network analysis allows the selection of ANN 6-8-1 (6 input neurons, 8 neurons in one hidden layer, 1 output neuron), which best predicts the parameter k. The selected ANN is presented in Figure 4. Activation functions were identified. For the hidden neurons, this was a tanh sigmoid curve, and for the output neuron, this was a logistic function. The errors measures (RMS, MAE, R2) for the developed ANN 6-8-1, in particular, and the subsets training Tr, testing T and validation V of set B are shown in Table 3. Prediction mean relative error Mean RE of the selected ANN 6-8-1 was a maximum of about ±4% in all subsets of data set B.

4.3. ANN Estimation

The proposed ANN 6-8-1 was tested by applying it to predict the permeability coefficient k on the basis of field data from set A. Variables X1, X2, X3 and X4 (Table 1) were introduced and Y = k values were estimated using the proposed ANN. Measured values of parameter k determined in pumping tests and values of parameter k predicted using developed ANN 6-8-1 are presented in Table 4. The value of the maximum single relative error Max RE was equal to 7.59%.

5. Conclusions

The paper concerns the determination of the permeability coefficient k using an ANN (Artificial Neural Network). In order to develop an ANN, field and laboratory tests of the parameter k were performed and index properties of the tested soils were determined. The research was performed in following soils: FSa, MSa and CSa. On the basis of the conducted research, it can be stated that similar values of the permeability coefficient k were obtained in field and laboratory tests from soils from the same test site. The influence of the relative density ID on the permeability coefficients k is significant in the tested soils. Higher values of permeability coefficients k were found in soils characterized by a lower relative density ID, higher void ratio e and lower value of effective soil diameter d10.
The proposed ANN with the 6-8-1 architecture predicts the real value of the permeability coefficient k based on the following data: soil type, relative density ID, void ratio e and effective soil diameter d10. The presented ANN estimates the permeability coefficients k with values of determination coefficient R2 = 0.97, mean relative error RE = ±4% and single maximum relative error Max RE = 7.59%.
The presented research frames the work as an experimental campaign and characterization of sandy soils. Further field and laboratory tests of the permeability coefficient k and the consideration of additional factors influencing the value of the tested parameter should be performed because the conclusions of the analysis refer only to the tested soils. The presented ANN model may be upgraded, new functions can be entered, including other types of soil, and can thus be extended to more types of non-cohesive and cohesive soils. The development of an ANN for the prediction of the soil permeability coefficient allows the reduction of the costs and time needed to conduct laboratory or field tests to determine this parameter.

Author Contributions

G.W. and A.M. prepared the research program, performed laboratory and field tests, developed an ANN and prepared the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets generated and analyzed during the current study are available from the authors upon reasonable request.

Conflicts of Interest

Authors declare no conflict of interest.

References

  1. Todd, D. Groundwater Hydrology, 2nd ed.; John Wiley & Sons: Chichester, UK, 1980. [Google Scholar]
  2. Wrzesiński, G. Anisotropy of soil shear strength parameters caused by the principal stress rotation. Arch. Civ. Eng. 2021, 67, 163–187. [Google Scholar] [CrossRef]
  3. Wrzesiński, G.; Lechowicz, Z. Testing of undrained shear strength in a Hollow Cylinder Apparatus. Studia Geotech. Mech. 2015, 37, 69–73. [Google Scholar] [CrossRef] [Green Version]
  4. Szymkiewicz, A.; Kryczałło, A. Calculation of permeability coefficient of sands and gravel based on grain size distribution curve: Review of empirical relations. Inżynieria Morska Geotech. 2011, 2, 110–121. [Google Scholar]
  5. Twardowski, K.; Drożdżak, R. Pośrednie metody oceny właściwości filtracyjnych gruntów. Wiertnictwo Nafta Gaz 2006, 23, 477–486. [Google Scholar]
  6. Parylak, K.; Zięba, Z.; Bułdys, A.; Witek, K. The verification of determining a permeability coefficient of non-cohesive soil based on empirical formulas including its microstructure. Acta Sci. Pol. Archit. 2013, 12, 43–51. [Google Scholar]
  7. Wrzesiński, G. Permeability coefficient tests in non-cohesive soils. Sci. Rev. Eng. Environ. Sci. 2020, 29, 72–80. [Google Scholar] [CrossRef] [Green Version]
  8. MacDonald, A.; Barker, J.; Davies, J. The bailer test: A simple effective pumping test for assessing borehole success. Hydrogeol. J. 2008, 16, 1065–1075. [Google Scholar] [CrossRef] [Green Version]
  9. Pawluk, K.; Połoński, M.; Lendo-Siwicka, M.; Wrzesiński, G. Two-objective optimization for optimal design of the multi-layered permeable reactive barriers. IOP Conference Series: Materials Science and Engineering. IOP Conf. Ser. Mater. Sci. Eng. 2019, 471, 1–11. [Google Scholar] [CrossRef]
  10. Polak, K.; Kaznowska-Opala, K.; Pawlecka, K.; Klich, J. Interpretation of pumping tests results on the basis of examination of AGH-1 well. Przegląd górniczy 2014, 10, 106–111. [Google Scholar]
  11. Head, K.; Epps, R. Manual of Soil Laboratory Testing. Vol. 2. Permeability, Shear Strength and Compressibility Test; Whittles Publishing: Scotland, UK, 2011. [Google Scholar]
  12. Jang, D.; Park, H.; Choi, G. Estimation of Leakage Ratio Using Principal Component Analysis and Artificial Neural Network in Water Distribution Systems. Sustainability 2018, 10, 750. [Google Scholar] [CrossRef] [Green Version]
  13. Kang, P.-S.; Lim, J.-S.; Huh, C. Artificial Neural Network Model to Estimate the Viscosity of Polymer Solutions for Enhanced Oil Recovery. Appl. Sci. 2016, 6, 188. [Google Scholar] [CrossRef] [Green Version]
  14. Li, W.; Cui, L.; Zhang, Y.; Cai, Z.; Zhang, M.; Xu, W.; Zhao, X.; Lei, Y.; Pan, X.; Li, J.; et al. Using a Backpropagation Artificial Neural Network to Predict Nutrient Removal in Tidal Flow Constructed Wetlands. Water 2018, 10, 83. [Google Scholar] [CrossRef] [Green Version]
  15. Zhu, X.; Fan, Y.; Zhang, F.; Ye, X.; Chen, C.; Yue, H. Multiple-Factor Based Sparse Urban Travel Time Prediction. Appl. Sci. 2018, 8, 279. [Google Scholar] [CrossRef] [Green Version]
  16. Paez, T.L. Neural networks in mechanical systems simulation, identification and assessment. Shock. Vib. 1993, 1, 177–199. [Google Scholar] [CrossRef]
  17. Adeli, H. Neural Networks in Civil Engineering: 1989–2000. Comput.-Aided Civ. Inf. 2001, 16, 126–142. [Google Scholar] [CrossRef]
  18. Ghaboussi, J. Biologically inspired soft computing methods in structural mechanics and engineering. Struct. Eng. Mech. 2001, 11, 485–502. [Google Scholar] [CrossRef]
  19. Rojas, R. Neural Networks: A Systematic Introduction; Springer: Berlin/Heidelberg, Germany, 1996. [Google Scholar]
  20. Wrzesiński, G.; Sulewska, M.J.; Lechowicz, Z. Evaluation of the Change in Undrained Shear Strength in Cohesive Soils due to Principal Stress Rotation Using an Artificial Neural Network. Appl. Sci. 2018, 8, 781. [Google Scholar] [CrossRef] [Green Version]
  21. Waszczyszyn, Z. Advances of Soft Computing in Engineering CISM Courses and Lectures; Springer: Vienna, Austria, 2010. [Google Scholar]
  22. Waszczyszyn, Z. Artificial neural networks in civil engineering: Another five years of research in Poland. Comp. Assist. Mech. Eng. Sci. 2011, 18, 131–146. [Google Scholar]
  23. Rafiq, M.Y.; Bugmann, G.; Easterbrook, D.J. Neural network design for engineering applications. Comput. Struct. 2001, 79, 1541–1552. [Google Scholar] [CrossRef]
  24. Shahin, M.A.; Jaksa, M.B.; Maier, H.R. Artificial Neural Network Applications in Geotechnical Engineering. Aust. Geomech. 2001, 36, 49–62. [Google Scholar]
  25. Shahin, M.A.; Jaksa, M.B.; Maier, H.R. State of the Art of Artificial Neural Networks in Geotechnical Engineering. Elect. J. Geotech. Eng. 2008, 8, 1–26. [Google Scholar]
  26. Dihoru, L.; Muir Wood, D.; Sadek, T.; Lings, M. A neural network for error prediction in a true triaxial apparatus with flexible boundaries. Comput. Geotech. 2005, 32, 59–71. [Google Scholar] [CrossRef]
  27. Pichler, B.; Lackner, R.; Mang, H.A. Back analysis of model parameters in geotechnical engineering by means of soft computing. Int. J. Numer. Meth. Eng. 2003, 57, 1943–1978. [Google Scholar] [CrossRef]
  28. Sulewska, M.J. Applying artificial neural networks for analysis of geotechnical problems. Comp. Assist. Mech. Eng. Sci. 2011, 18, 231–241. [Google Scholar]
  29. Tian, J.; Li, C.; Liu, J.; Yu, F.; Cheng, S.; Zhao, N.; Wan Jaafar, W.Z. Groundwater Depth Prediction Using Data-Driven Models with the Assistance of Gamma Test. Sustainability 2016, 8, 1076. [Google Scholar] [CrossRef] [Green Version]
  30. Zhou, T.; Wang, F.; Yang, Z. Comparative Analysis of ANN and SVM Models Combined with Wavelet Preprocess for Groundwater Depth Prediction. Water 2017, 9, 781. [Google Scholar] [CrossRef] [Green Version]
  31. Ellis, G.W.; Yao, C.; Zhao, R.; Penumadu, D. Stress-strain modeling of sands using artificial neural network. J. Geotech. Eng. 1995, 121, 429–435. [Google Scholar] [CrossRef]
  32. Sidarta, D.E.; Ghaboussi, J. Constitutive modeling of geomaterials from non-uniform material test. Comput. Geotech. 1998, 22, 53–71. [Google Scholar] [CrossRef]
  33. Penumadu, D.; Zhao, R. Triaxial compression behavior of sand and gravel using artificial neural networks (ANN). Comput. Geot. 1999, 24, 207–230. [Google Scholar] [CrossRef]
  34. Basheer, I.A. Selection of methodology for neural network modeling of constitutive hysteresis behavior of soils. Comput.-Aided Civ. Inf. 2000, 15, 440–458. [Google Scholar] [CrossRef]
  35. Najjar, Y.M.; Huang, C. Simulating the stress-strain behavior of Georgia kaolin via recurrent neuronet approach. Comput. Geot. 2007, 34, 346–361. [Google Scholar] [CrossRef]
  36. Fu, Q.; Hashash, Y.M.A.; Hung, S.; Ghaboussi, J. Integration of laboratory testing and constitutive modeling of soils. Comput. Geotech. 2007, 34, 330–345. [Google Scholar] [CrossRef]
  37. Lee, S.J.; Lee, S.R.; Kim, Y.S. An approach to estimate unsaturated shear strength using artificial neural network and hyperbolic formulation. Comput. Geotech. 2003, 30, 489–503. [Google Scholar] [CrossRef]
  38. Das, S.K.; Basudhar, P.K. Undrained lateral load capacity of piles in clay using artificial neural network. Comput. Geotech. 2006, 33, 454–459. [Google Scholar] [CrossRef]
  39. Byeon, W.Y.; Lee, S.R.; Kim, Y.S. Application of flat DMT and ANN to Korean soft clay deposits for reliable estimation of undrained shear strength. Int. J. Offshore Polar Eng. 2006, 16, 73–80. [Google Scholar]
  40. Wrzesiński, G. Stability Analysis of an Embankment with Influence of the Principal Stress Rotation on the Shear Strength of Subsoil. Ph.D. Thesis, Warsaw University of Life Sciences–SGGW, Warsaw, Poland, 2016. [Google Scholar]
  41. Wrzesiński, G.; Kowalski, J.; Miszkowska, A. Numerical analysis of dewatering process of deep excavation. In Proceedings of the 18th International Multidisciplinary Scientific Geoconference SGEM 2018: Hydrogeology, Engineering Geology and Geotechnics, Albena, Bulgaria, 30 June 2018; Issue 1.2, Science and Technologies in Geology, Exploration and Mining. Volume 18, pp. 497–504. [Google Scholar]
  42. Lin, P.Y.; Ni, P.P.; Guo, C.C.; Mei, G.X. Mapping Soil Nail Loads Using FHWA Simplified Models and Artificial Neural Network Technique. Can. Geotech. J. 2020, 57, 1453–1471. [Google Scholar] [CrossRef]
  43. Lin, P.Y.; Chen, X.Y.; Jiang, M.J.; Song, X.G.; Xu, M.J.; Huang, S. Mapping shear strength and compressibility of soft soils with artificial neural networks. Eng. Geol. 2022, 300, 106585. [Google Scholar] [CrossRef]
  44. Liu, D.; Lin, P.Y.; Zhao, C.Y.; Qiu, J.J. Mapping horizontal displacement of soil nail walls using machine learning approaches. Acta Geotech. 2021, 16, 4027–4044. [Google Scholar] [CrossRef]
  45. Driscoll, F. Groundwater and Wells, 2nd ed.; Johnson Filtration Systems Inc.: St Paul, MN, USA, 1986. [Google Scholar]
  46. ICRC. Technical Review: Practical Guidelines for Test Pumping in Water Wells; International Committee of the Red Cross: Geneva, The Switzerland, 2011. [Google Scholar]
  47. Kruseman, G.P.; de Ridder, N.A. Analysis and Evaluation of Pumping Test Data, 2nd ed.; Publication 47; International Institute for Land Reclamation and Improvement: Wageningen, The Netherlands, 1994. [Google Scholar]
  48. EN ISO 14688-1; Geotechnical Investigation and Testing-Identification and Classification of Soil-Part 1: Identification and Description. International Organization for Standardization: Geneva, Switzerland, 2002.
  49. EN ISO 14688-2; Geotechnical Investigation and Testing-Identification and Classification of Soil-Part 2: Principles for a Classification. International Organization for Standardization: Geneva, Switzerland, 2004.
  50. Head, K. Manual of Soil Laboratory Testing. Vol. 1. Soil Classification and Compaction Test; Pentech Press: London, UK, 1996. [Google Scholar]
  51. Tymosiak, D.; Sulewska, M.J. The study of compactibility parameters in non-cohesive soils by Proctor compaction test. Acta Sci. Pol. Archit. 2016, 15, 43–54. [Google Scholar]
  52. Bishop, C.M. Neural Networks for Pattern Recognition; Oxford University Press: Oxford, UK, 1995. [Google Scholar]
Figure 1. Pumping test system scheme.
Figure 1. Pumping test system scheme.
Sustainability 14 06736 g001
Figure 2. Consolidometer used in tests.
Figure 2. Consolidometer used in tests.
Sustainability 14 06736 g002
Figure 3. Chemical composition spectra and sample SEM photos showing the structure for FSa, characterized by similar relative density ID from two test sites (wells no. 1 and 6): (a) FSa from well no. 1, (b) FSa from well no. 6.
Figure 3. Chemical composition spectra and sample SEM photos showing the structure for FSa, characterized by similar relative density ID from two test sites (wells no. 1 and 6): (a) FSa from well no. 1, (b) FSa from well no. 6.
Sustainability 14 06736 g003
Figure 4. Architecture of the developed ANN 6-8-1.
Figure 4. Architecture of the developed ANN 6-8-1.
Sustainability 14 06736 g004
Table 1. Grain size distribution, relative densities, void ratios and effective soil diameters d10 of soils.
Table 1. Grain size distribution, relative densities, void ratios and effective soil diameters d10 of soils.
No. of WellSoilFraction-EN ISO 14688-1: 2002; EN ISO 14688-2: 2004
(%)
Relative Density
ID (-)
Void Ratio
e (-)
Effective Soil Diameter
d10 (mm)
GrSaSiCl
1FSa091900.490.5420.07
2FSa190900.670.5910.07
3FSa092800.610.4980.08
4FSa192700.640.5230.09
5FSa290800.410.6560.08
6FSa094600.540.5990.10
7FSa193600.510.5870.10
8FSa097300.560.4860.17
9FSa095320.390.7050.17
10FSa095500.500.5890.12
11FSa0861220.350.7120.04
12FSa2861020.680.5060.063
13FSa189820.390.7010.09
14FSa1881100.690.5200.06
15FSa094600.220.7280.11
16FSa1851130.540.5990.06
17FSa290800.390.6800.08
18FSa089920.420.7010.07
19FSa290800.800.4910.08
20FSa1871020.450.5930.063
21MSa099100.480.6030.25
22MSa097210.410.5890.23
23MSa196300.580.5210.20
24MSa097300.520.5800.20
25MSa098200.610.5360.21
26MSa294310.650.5140.21
27MSa292420.700.5450.18
28MSa195400.670.5130.17
29MSa098200.560.5630.21
30MSa198100.780.4560.25
31MSa097210.370.6310.23
32MSa095320.850.4050.20
33MSa197200.780.4150.24
34MSa296200.630.5200.24
35MSa198100.560.5470.26
36MSa098200.820.4360.21
37MSa197200.3906530.24
38MSa097210.330.6060.23
39MSa096400.860.4100.25
40MSa098200.610.5350.21
41MSa296200.540.5180.23
42MSa098200.490.5520.21
43CSa892000.710.4700.50
44CSa1287100.680.4560.46
45CSa1981000.590.5070.61
46CSa1882000.400.5320.58
47CSa1780210.540.5280.45
48CSa1882000.280.5900.57
49CSa1882000.380.5760.60
50CSa1682200.920.4230.63
Table 2. Permeability coefficients k, evaluated based on pumping and consolidometer tests.
Table 2. Permeability coefficients k, evaluated based on pumping and consolidometer tests.
No. of WellSoilPermeability Coefficient k (m/s)
Pumping
Test
Consolidometer Test
1FSa2.32 × 10−52.20 × 10−5
2FSa3.69 × 10−53.42 × 10−5
3FSa2.10 × 10−52.00 × 10−5
4FSa1.24 × 10−51.33 × 10−5
5FSa5.77 × 10−55.64 × 10−5
6FSa4.68 × 10−54.33 × 10−5
7FSa3.79 × 10−53.66 × 10−5
8FSa4.40 × 10−53.97 × 10−5
9FSa4.78 × 10−54.64 × 10−5
10FSa5.59 × 10−55.27 × 10−5
11FSa9.32 × 10−59.05 × 10−5
12FSa3.85 × 10−53.64 × 10−5
13FSa8.48 × 10−58.50 × 10−5
14FSa4.54 × 10−54.63 × 10−5
15FSa9.86 × 10−59.65 × 10−5
16FSa5.08 × 10−54.96 × 10−5
17FSa7.20 × 10−57.32 × 10−5
18FSa8.64 × 10−58.73 × 10−5
19FSa5.30 × 10−54.98 × 10−5
20FSa6.75 × 10−56.14 × 10−5
21MSa1.69 × 10−41.57 × 10−4
22MSa2.97 × 10−42.93 × 10−4
23MSa2.28 × 10−42.12 × 10−4
24MSa1.49 × 10−41.45 × 10−4
25MSa1.32 × 10−41.32 × 10−4
26MSa1.35 × 10−41.25 × 10−4
27MSa1.48 × 10−41.50 × 10−4
28MSa1.36 × 10−41.22 × 10−4
29MSa2.20 × 10−42.03 × 10−4
30MSa1.17 × 10−41.16 × 10−4
31MSa2.78 × 10−42.65 × 10−4
32MSa1.45 × 10−41.26 × 10−4
33MSa1.63 × 10−41.55 × 10−4
34MSa2.08 × 10−41.92 × 10−4
35MSa2.23 × 10−42.05 × 10−4
36MSa1.85 × 10−41.64 × 10−4
37MSa2.89 × 10−42.57 × 10−4
38MSa2.54 × 10−42.36 × 10−4
39MSa1.29 × 10−41.16 × 10−4
40MSa1.98 × 10−41.97 × 10−4
41MSa1.75 × 10−41.63 × 10−4
42MSa1.70 × 10−41.61 × 10−4
43CSa3.73 × 10−43.68 × 10−4
44CSa4.14 × 10−43.84 × 10−4
45CSa4.85 × 10−44.78 × 10−4
46CSa6.28 × 10−45.89 × 10−4
47CSa5.84 × 10−45.80 × 10−4
48CSa7.05 × 10−46.98 × 10−4
49CSa6.97 × 10−47.02 × 10−4
50CSa3.24 × 10−43.29 × 10−4
Table 3. Errors measures for ANN 6-8-1 in the subsets training Tr, testing T, validation V of set B.
Table 3. Errors measures for ANN 6-8-1 in the subsets training Tr, testing T, validation V of set B.
ErrorsSubset TrSubset TSubset V
RMS0.00980.00960.0084
MAE0.02150.02040.0119
R20.9760.9760.976
Table 4. Measured values of permeability coefficient k from field tests and predicted values of k using developed ANN 6-8-1.
Table 4. Measured values of permeability coefficient k from field tests and predicted values of k using developed ANN 6-8-1.
No. of WellsSoilMeasured
Values of k
in Pumping Tests
di (m/s)
Predicted
Values of k Based on ANN 6-8-1
yi (m/s)
Relative
Errors
of Individual
Case
REi (%)
1FSa2.32 × 10−52.25 × 10−53.02
2FSa3.69 × 10−53.58 × 10−52.98
3FSa2.10 × 10−52.23 × 10−56.19
4FSa1.24 × 10−51.31 × 10−55.65
5FSa5.77 × 10−55.70 × 10−51.21
6FSa4.68 × 10−54.65 × 10−50.64
7FSa3.79 × 10−53.80 × 10−50.26
8FSa4.40 × 10−54.28 × 10−52.73
9FSa4.78 × 10−54.69 × 10−51.88
10FSa5.59 × 10−55.45 × 10−52.50
11FSa9.32 × 10−59.28 × 10−50.43
12FSa3.85 × 10−53.84 × 10−50.26
13FSa8.48 × 10−58.48 × 10−50
14FSa4.54 × 10−54.53 × 10−50.22
15FSa9.86 × 10−59.82 × 10−50.41
16FSa5.08 × 10−55.03 × 10−50.98
17FSa7.20 × 10−56.98 × 10−53.06
18FSa8.64 × 10−58.63 × 10−50.12
19FSa5.30 × 10−55.46 × 10−53.02
20FSa6.75 × 10−56.84 × 10−51.33
21MSa1.69 × 10−41.65 × 10−42.37
22MSa2.97 × 10−43.05 × 10−42.69
23MSa2.28 × 10−42.22 × 10−42.63
24MSa1.49 × 10−41.49 × 10−40
25MSa1.32 × 10−41.31 × 10−40.76
26MSa1.35 × 10−41.33 × 10−41.48
27MSa1.48 × 10−41.47 × 10−40.68
28MSa1.36 × 10−41.30 × 10−44.41
29MSa2.20 × 10−42.18 × 10−40.91
30MSa1.17 × 10−41.16 × 10−40.85
31MSa2.78 × 10−42.73 × 10−41.80
32MSa1.45 × 10−41.34 × 10−47.59
33MSa1.63 × 10−41.60 × 10−41.84
34MSa2.08 × 10−42.05 × 10−41.44
35MSa2.23 × 10−42.28 × 10−42.24
36MSa1.85 × 10−41.84 × 10−40.54
37MSa2.89 × 10−42.89 × 10−40
38MSa2.54 × 10−42.53 × 10−40.39
39MSa1.29 × 10−41.21 × 10−46.20
40MSa1.98 × 10−41.99 × 10−40.51
41MSa1.75 × 10−41.75 × 10−40
42MSa1.70 × 10−41.68 × 10−41.18
43CSa3.73 × 10−43.71 × 10−40.54
44CSa4.14 × 10−44.25 × 10−42.66
45CSa4.85 × 10−44. 73 × 10−42.47
46CSa6.28 × 10−46.62 × 10−45.41
47CSa5.84 × 10−45.84 × 10−40
48CSa7.05 × 10−47.03 × 10−40.28
49CSa6.97 × 10−46.75 × 10−43.16
50CSa3.24 × 10−43.20 × 10−41.23
Max
RE32 = 7.59%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wrzesiński, G.; Markiewicz, A. Prediction of Permeability Coefficient k in Sandy Soils Using ANN. Sustainability 2022, 14, 6736. https://0-doi-org.brum.beds.ac.uk/10.3390/su14116736

AMA Style

Wrzesiński G, Markiewicz A. Prediction of Permeability Coefficient k in Sandy Soils Using ANN. Sustainability. 2022; 14(11):6736. https://0-doi-org.brum.beds.ac.uk/10.3390/su14116736

Chicago/Turabian Style

Wrzesiński, Grzegorz, and Anna Markiewicz. 2022. "Prediction of Permeability Coefficient k in Sandy Soils Using ANN" Sustainability 14, no. 11: 6736. https://0-doi-org.brum.beds.ac.uk/10.3390/su14116736

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop