Approximation of the exponential integral (well function) using sampling methods

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

Exponential integral (also known as well function) is often used in hydrogeology to solve Theis and Hantush equations. Many methods have been developed to approximate the exponential integral. Most of these methods are based on numerical approximations and are valid for a certain range of the argument value.This paper presents a new approach to approximate the exponential integral. The new approach is based on sampling methods. Three different sampling methods; Latin Hypercube Sampling (LHS), Orthogonal Array (OA), and Orthogonal Array-based Latin Hypercube (OA-LH) have been used to approximate the function. Different argument values, covering a wide range, have been used. The results of sampling methods were compared with results obtained by Mathematica software, which was used as a benchmark.All three sampling methods converge to the result obtained by Mathematica, at different rates. It was found that the orthogonal array (OA) method has the fastest convergence rate compared with LHS and OA-LH. The root mean square error RMSE of OA was in the order of 1E-08. This method can be used with any argument value, and can be used to solve other integrals in hydrogeology such as the leaky aquifer integral.

Original languageEnglish
Pages (from-to)278-282
Number of pages5
JournalJournal of Hydrology
Volume523
DOIs
Publication statusPublished - 1 Apr 2015

Fingerprint

sampling
hydrogeology
method
aquifer
software
rate

Keywords

  • Exponential integral
  • Latin hypercube sampling
  • Orthogonal array
  • Orthogonal array-based latin hypercube sampling
  • Well function

ASJC Scopus subject areas

  • Water Science and Technology

Cite this

Approximation of the exponential integral (well function) using sampling methods. / Baalousha, Husam.

In: Journal of Hydrology, Vol. 523, 01.04.2015, p. 278-282.

Research output: Contribution to journalArticle

@article{f518d7fb771842e38a6d843cab8cee3b,
title = "Approximation of the exponential integral (well function) using sampling methods",
abstract = "Exponential integral (also known as well function) is often used in hydrogeology to solve Theis and Hantush equations. Many methods have been developed to approximate the exponential integral. Most of these methods are based on numerical approximations and are valid for a certain range of the argument value.This paper presents a new approach to approximate the exponential integral. The new approach is based on sampling methods. Three different sampling methods; Latin Hypercube Sampling (LHS), Orthogonal Array (OA), and Orthogonal Array-based Latin Hypercube (OA-LH) have been used to approximate the function. Different argument values, covering a wide range, have been used. The results of sampling methods were compared with results obtained by Mathematica software, which was used as a benchmark.All three sampling methods converge to the result obtained by Mathematica, at different rates. It was found that the orthogonal array (OA) method has the fastest convergence rate compared with LHS and OA-LH. The root mean square error RMSE of OA was in the order of 1E-08. This method can be used with any argument value, and can be used to solve other integrals in hydrogeology such as the leaky aquifer integral.",
keywords = "Exponential integral, Latin hypercube sampling, Orthogonal array, Orthogonal array-based latin hypercube sampling, Well function",
author = "Husam Baalousha",
year = "2015",
month = "4",
day = "1",
doi = "10.1016/j.jhydrol.2015.01.062",
language = "English",
volume = "523",
pages = "278--282",
journal = "Journal of Hydrology",
issn = "0022-1694",
publisher = "Elsevier",

}

TY - JOUR

T1 - Approximation of the exponential integral (well function) using sampling methods

AU - Baalousha, Husam

PY - 2015/4/1

Y1 - 2015/4/1

N2 - Exponential integral (also known as well function) is often used in hydrogeology to solve Theis and Hantush equations. Many methods have been developed to approximate the exponential integral. Most of these methods are based on numerical approximations and are valid for a certain range of the argument value.This paper presents a new approach to approximate the exponential integral. The new approach is based on sampling methods. Three different sampling methods; Latin Hypercube Sampling (LHS), Orthogonal Array (OA), and Orthogonal Array-based Latin Hypercube (OA-LH) have been used to approximate the function. Different argument values, covering a wide range, have been used. The results of sampling methods were compared with results obtained by Mathematica software, which was used as a benchmark.All three sampling methods converge to the result obtained by Mathematica, at different rates. It was found that the orthogonal array (OA) method has the fastest convergence rate compared with LHS and OA-LH. The root mean square error RMSE of OA was in the order of 1E-08. This method can be used with any argument value, and can be used to solve other integrals in hydrogeology such as the leaky aquifer integral.

AB - Exponential integral (also known as well function) is often used in hydrogeology to solve Theis and Hantush equations. Many methods have been developed to approximate the exponential integral. Most of these methods are based on numerical approximations and are valid for a certain range of the argument value.This paper presents a new approach to approximate the exponential integral. The new approach is based on sampling methods. Three different sampling methods; Latin Hypercube Sampling (LHS), Orthogonal Array (OA), and Orthogonal Array-based Latin Hypercube (OA-LH) have been used to approximate the function. Different argument values, covering a wide range, have been used. The results of sampling methods were compared with results obtained by Mathematica software, which was used as a benchmark.All three sampling methods converge to the result obtained by Mathematica, at different rates. It was found that the orthogonal array (OA) method has the fastest convergence rate compared with LHS and OA-LH. The root mean square error RMSE of OA was in the order of 1E-08. This method can be used with any argument value, and can be used to solve other integrals in hydrogeology such as the leaky aquifer integral.

KW - Exponential integral

KW - Latin hypercube sampling

KW - Orthogonal array

KW - Orthogonal array-based latin hypercube sampling

KW - Well function

UR - http://www.scopus.com/inward/record.url?scp=84922636311&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84922636311&partnerID=8YFLogxK

U2 - 10.1016/j.jhydrol.2015.01.062

DO - 10.1016/j.jhydrol.2015.01.062

M3 - Article

AN - SCOPUS:84922636311

VL - 523

SP - 278

EP - 282

JO - Journal of Hydrology

JF - Journal of Hydrology

SN - 0022-1694

ER -