Xiaoyi Liu - Academia.edu (original) (raw)
Papers by Xiaoyi Liu
Water Resources Research, 2009
1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m i... more 1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m in deep boreholes at the Mizunami Underground Research Laboratory (MIU) construction site in central Japan. During these two tests, induced groundwater responses were monitored at many observation intervals at various depths in different boreholes at the site. We analyze the two cross-hole pumping tests using transient hydraulic tomography (THT) based on an efficient sequential successive linear estimator to compute the hydraulic conductivity (K) and specific storage (S s ) tomograms, as well as their uncertainties in three dimensions. The equivalent K and S s estimates obtained using asymptotic analysis treating the medium to be homogeneous served as the mean parameter estimates for the 3-D stochastic inverse modeling effort. Results show several, distinct, high K and low S s zones that are continuous over hundreds of meters, which appear to delineate fault zones and their connectivity. The THT analysis of the tests also identified a low K zone which corresponds with a known fault zone trending NNW and has been found to compartmentalize groundwater flow at the site. These results corroborate well with observed water level records, available fault information, and coseismic groundwater level responses during several large earthquakes. The successful application of THT to cross-hole pumping tests conducted in fractured granite at this site suggests that THT is a promising approach to delineate large-scale K and S s heterogeneities, fracture connectivity, and to quantify uncertainty of the estimated fields.
Journal of Hydrology, 2007
Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hy... more Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hydraulic heterogeneity. To date, a comprehensive validation of hydraulic tomography has not been done either at the laboratory or field scales. The main objective of this paper is to examine the accuracy of hydraulic conductivity (K) tomograms obtained from the steady-state hydraulic tomography algorithm of [Yeh, T.-C. J., Liu, S., 2000. Hydraulic tomography: development of a new aquifer test method. Water Resources Research 36, 2095-2105]. We first obtain a reference K tomogram through the inversion of synthetic cross-hole test data generated through numerical simulations. The purpose of reference K tomogram generation is to examine the ability of the algorithm to image the heterogeneity pattern under optimal conditions without experimental errors and with full control of forcing functions (initial and boundary conditions as well as source/sink terms). Parallel to the generation of synthetic data, we conduct hydraulic tests at multiple scales in a laboratory aquifer with deterministic heterogeneity to generate data that are used to validate K tomograms from hydraulic tomography. Measurements include multiple K estimates from core, slug, single-hole and cross-hole tests as well as several unidirectional, flow-through experiments conducted on the sandbox under steady-state conditions. Validation of K tomograms involved a multi-method and multiscale approach proposed herein which include: (1) visual comparisons of K tomograms to the true sand distributions and the reference K tomogram; (2) testing the ability of 0022-1694/$ -see front matter Published by Elsevier B.V. a v a i l a b l e a t w w w . s c i e n c e d i r e c t . c o m j o u r n a l h o m e p a g e : w w w . e l s e v i e r . c o m / l o c a t e / j h y d r o l
Water Resources Research, 2009
1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m i... more 1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m in deep boreholes at the Mizunami Underground Research Laboratory (MIU) construction site in central Japan. During these two tests, induced groundwater responses were monitored at many observation intervals at various depths in different boreholes at the site. We analyze the two cross-hole pumping tests using transient hydraulic tomography (THT) based on an efficient sequential successive linear estimator to compute the hydraulic conductivity (K) and specific storage (S s ) tomograms, as well as their uncertainties in three dimensions. The equivalent K and S s estimates obtained using asymptotic analysis treating the medium to be homogeneous served as the mean parameter estimates for the 3-D stochastic inverse modeling effort. Results show several, distinct, high K and low S s zones that are continuous over hundreds of meters, which appear to delineate fault zones and their connectivity. The THT analysis of the tests also identified a low K zone which corresponds with a known fault zone trending NNW and has been found to compartmentalize groundwater flow at the site. These results corroborate well with observed water level records, available fault information, and coseismic groundwater level responses during several large earthquakes. The successful application of THT to cross-hole pumping tests conducted in fractured granite at this site suggests that THT is a promising approach to delineate large-scale K and S s heterogeneities, fracture connectivity, and to quantify uncertainty of the estimated fields.
Water Resources Research, 2009
1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m i... more 1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m in deep boreholes at the Mizunami Underground Research Laboratory (MIU) construction site in central Japan. During these two tests, induced groundwater responses were monitored at many observation intervals at various depths in different boreholes at the site. We analyze the two cross-hole pumping tests using transient hydraulic tomography (THT) based on an efficient sequential successive linear estimator to compute the hydraulic conductivity (K) and specific storage (S s ) tomograms, as well as their uncertainties in three dimensions. The equivalent K and S s estimates obtained using asymptotic analysis treating the medium to be homogeneous served as the mean parameter estimates for the 3-D stochastic inverse modeling effort. Results show several, distinct, high K and low S s zones that are continuous over hundreds of meters, which appear to delineate fault zones and their connectivity. The THT analysis of the tests also identified a low K zone which corresponds with a known fault zone trending NNW and has been found to compartmentalize groundwater flow at the site. These results corroborate well with observed water level records, available fault information, and coseismic groundwater level responses during several large earthquakes. The successful application of THT to cross-hole pumping tests conducted in fractured granite at this site suggests that THT is a promising approach to delineate large-scale K and S s heterogeneities, fracture connectivity, and to quantify uncertainty of the estimated fields.
Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hy... more Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hydraulic heterogeneity. To date, a comprehensive validation of hydraulic tomography has not been done either at the laboratory or field scales. The main objective of this paper is to examine the accuracy of hydraulic conductivity (K) tomograms obtained from the steady-state hydraulic tomography algorithm of [Yeh, T.-C. J., Liu, S., 2000. Hydraulic tomography: development of a new aquifer test method. Water Resources Research 36, 2095-2105]. We first obtain a reference K tomogram through the inversion of synthetic cross-hole test data generated through numerical simulations. The purpose of reference K tomogram generation is to examine the ability of the algorithm to image the heterogeneity pattern under optimal conditions without experimental errors and with full control of forcing functions (initial and boundary conditions as well as source/sink terms). Parallel to the generation of synthetic data, we conduct hydraulic tests at multiple scales in a laboratory aquifer with deterministic heterogeneity to generate data that are used to validate K tomograms from hydraulic tomography. Measurements include multiple K estimates from core, slug, single-hole and cross-hole tests as well as several unidirectional, flow-through experiments conducted on the sandbox under steady-state conditions. Validation of K tomograms involved a multi-method and multiscale approach proposed herein which include: (1) visual comparisons of K tomograms to the true sand distributions and the reference K tomogram; (2) testing the ability of 0022-1694/$ -see front matter Published by Elsevier B.V. a v a i l a b l e a t w w w . s c i e n c e d i r e c t . c o m j o u r n a l h o m e p a g e : w w w . e l s e v i e r . c o m / l o c a t e / j h y d r o l
Journal of Hydrology, 2007
Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hy... more Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hydraulic heterogeneity. To date, a comprehensive validation of hydraulic tomography has not been done either at the laboratory or field scales. The main objective of this paper is to examine the accuracy of hydraulic conductivity (K) tomograms obtained from the steady-state hydraulic tomography algorithm of [Yeh, T.-C. J., Liu, S., 2000. Hydraulic tomography: development of a new aquifer test method. Water Resources Research 36, 2095-2105]. We first obtain a reference K tomogram through the inversion of synthetic cross-hole test data generated through numerical simulations. The purpose of reference K tomogram generation is to examine the ability of the algorithm to image the heterogeneity pattern under optimal conditions without experimental errors and with full control of forcing functions (initial and boundary conditions as well as source/sink terms). Parallel to the generation of synthetic data, we conduct hydraulic tests at multiple scales in a laboratory aquifer with deterministic heterogeneity to generate data that are used to validate K tomograms from hydraulic tomography. Measurements include multiple K estimates from core, slug, single-hole and cross-hole tests as well as several unidirectional, flow-through experiments conducted on the sandbox under steady-state conditions. Validation of K tomograms involved a multi-method and multiscale approach proposed herein which include: (1) visual comparisons of K tomograms to the true sand distributions and the reference K tomogram; (2) testing the ability of 0022-1694/$ -see front matter Published by Elsevier B.V. a v a i l a b l e a t w w w . s c i e n c e d i r e c t . c o m j o u r n a l h o m e p a g e : w w w . e l s e v i e r . c o m / l o c a t e / j h y d r o l
Water Resources Research, 2009
1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m i... more 1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m in deep boreholes at the Mizunami Underground Research Laboratory (MIU) construction site in central Japan. During these two tests, induced groundwater responses were monitored at many observation intervals at various depths in different boreholes at the site. We analyze the two cross-hole pumping tests using transient hydraulic tomography (THT) based on an efficient sequential successive linear estimator to compute the hydraulic conductivity (K) and specific storage (S s ) tomograms, as well as their uncertainties in three dimensions. The equivalent K and S s estimates obtained using asymptotic analysis treating the medium to be homogeneous served as the mean parameter estimates for the 3-D stochastic inverse modeling effort. Results show several, distinct, high K and low S s zones that are continuous over hundreds of meters, which appear to delineate fault zones and their connectivity. The THT analysis of the tests also identified a low K zone which corresponds with a known fault zone trending NNW and has been found to compartmentalize groundwater flow at the site. These results corroborate well with observed water level records, available fault information, and coseismic groundwater level responses during several large earthquakes. The successful application of THT to cross-hole pumping tests conducted in fractured granite at this site suggests that THT is a promising approach to delineate large-scale K and S s heterogeneities, fracture connectivity, and to quantify uncertainty of the estimated fields.
Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hy... more Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hydraulic heterogeneity. To date, a comprehensive validation of hydraulic tomography has not been done either at the laboratory or field scales. The main objective of this paper is to examine the accuracy of hydraulic conductivity (K) tomograms obtained from the steady-state hydraulic tomography algorithm of [Yeh, T.-C. J., Liu, S., 2000. Hydraulic tomography: development of a new aquifer test method. Water Resources Research 36, 2095-2105]. We first obtain a reference K tomogram through the inversion of synthetic cross-hole test data generated through numerical simulations. The purpose of reference K tomogram generation is to examine the ability of the algorithm to image the heterogeneity pattern under optimal conditions without experimental errors and with full control of forcing functions (initial and boundary conditions as well as source/sink terms). Parallel to the generation of synthetic data, we conduct hydraulic tests at multiple scales in a laboratory aquifer with deterministic heterogeneity to generate data that are used to validate K tomograms from hydraulic tomography. Measurements include multiple K estimates from core, slug, single-hole and cross-hole tests as well as several unidirectional, flow-through experiments conducted on the sandbox under steady-state conditions. Validation of K tomograms involved a multi-method and multiscale approach proposed herein which include: (1) visual comparisons of K tomograms to the true sand distributions and the reference K tomogram; (2) testing the ability of 0022-1694/$ -see front matter Published by Elsevier B.V. a v a i l a b l e a t w w w . s c i e n c e d i r e c t . c o m j o u r n a l h o m e p a g e : w w w . e l s e v i e r . c o m / l o c a t e / j h y d r o l
Ground Water, 2007
Hydraulic tomography has been developed as an alternative to traditional geostatistical methods t... more Hydraulic tomography has been developed as an alternative to traditional geostatistical methods to delineate heterogeneity patterns in parameters such as hydraulic conductivity (K) and specific storage (Ss). During hydraulic tomography surveys, a large number of hydraulic head data are collected from a series of cross-hole tests in the subsurface. These head data are then used to interpret the spatial distribution of K and Ss using inverse modeling. Here, we use the Sequential Successive Linear Estimator (SSLE) of Yeh and Liu (2000) to interpret synthetic pumping test data created through numerical simulations and real data generated in a laboratory sandbox aquifer to obtain the K tomograms. Here, we define “K tomogram” as an image of K distribution of the subsurface (or the inverse results) obtained via hydraulic tomography. We examine the influence of signal-to-noise ratio and biases on results using inverse modeling of synthetic and real cross-hole pumping test data. To accomplish this, we first show that the pumping rate, which affects the signal-to-noise ratio, and the order of data included into the SSLE algorithm both have large impacts on the quality of the K tomograms. We then examine the role of conditioning on the K tomogram and find that conditioning can improve the quality of the K tomogram, but can also impair it, if the data are of poor quality and conditioning data have a larger support volume than the numerical grid used to conduct the inversion. Overall, these results show that the quality of the K tomogram depends on the design of pumping tests, their conduct, the order in which they are included in the inverse code, and the quality as well as the support volume of additional data that are used in its computation.
Ground Water, 2007
Hydraulic tomography has been developed as an alternative to traditional geostatistical methods t... more Hydraulic tomography has been developed as an alternative to traditional geostatistical methods to delineate heterogeneity patterns in parameters such as hydraulic conductivity (K) and specific storage (Ss). During hydraulic tomography surveys, a large number of hydraulic head data are collected from a series of cross-hole tests in the subsurface. These head data are then used to interpret the spatial distribution of K and Ss using inverse modeling. Here, we use the Sequential Successive Linear Estimator (SSLE) of Yeh and Liu (2000) to interpret synthetic pumping test data created through numerical simulations and real data generated in a laboratory sandbox aquifer to obtain the K tomograms. Here, we define “K tomogram” as an image of K distribution of the subsurface (or the inverse results) obtained via hydraulic tomography. We examine the influence of signal-to-noise ratio and biases on results using inverse modeling of synthetic and real cross-hole pumping test data. To accomplish this, we first show that the pumping rate, which affects the signal-to-noise ratio, and the order of data included into the SSLE algorithm both have large impacts on the quality of the K tomograms. We then examine the role of conditioning on the K tomogram and find that conditioning can improve the quality of the K tomogram, but can also impair it, if the data are of poor quality and conditioning data have a larger support volume than the numerical grid used to conduct the inversion. Overall, these results show that the quality of the K tomogram depends on the design of pumping tests, their conduct, the order in which they are included in the inverse code, and the quality as well as the support volume of additional data that are used in its computation.
Water Resources Research, 2009
1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m i... more 1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m in deep boreholes at the Mizunami Underground Research Laboratory (MIU) construction site in central Japan. During these two tests, induced groundwater responses were monitored at many observation intervals at various depths in different boreholes at the site. We analyze the two cross-hole pumping tests using transient hydraulic tomography (THT) based on an efficient sequential successive linear estimator to compute the hydraulic conductivity (K) and specific storage (S s ) tomograms, as well as their uncertainties in three dimensions. The equivalent K and S s estimates obtained using asymptotic analysis treating the medium to be homogeneous served as the mean parameter estimates for the 3-D stochastic inverse modeling effort. Results show several, distinct, high K and low S s zones that are continuous over hundreds of meters, which appear to delineate fault zones and their connectivity. The THT analysis of the tests also identified a low K zone which corresponds with a known fault zone trending NNW and has been found to compartmentalize groundwater flow at the site. These results corroborate well with observed water level records, available fault information, and coseismic groundwater level responses during several large earthquakes. The successful application of THT to cross-hole pumping tests conducted in fractured granite at this site suggests that THT is a promising approach to delineate large-scale K and S s heterogeneities, fracture connectivity, and to quantify uncertainty of the estimated fields.
Journal of Hydrology, 2007
Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hy... more Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hydraulic heterogeneity. To date, a comprehensive validation of hydraulic tomography has not been done either at the laboratory or field scales. The main objective of this paper is to examine the accuracy of hydraulic conductivity (K) tomograms obtained from the steady-state hydraulic tomography algorithm of [Yeh, T.-C. J., Liu, S., 2000. Hydraulic tomography: development of a new aquifer test method. Water Resources Research 36, 2095-2105]. We first obtain a reference K tomogram through the inversion of synthetic cross-hole test data generated through numerical simulations. The purpose of reference K tomogram generation is to examine the ability of the algorithm to image the heterogeneity pattern under optimal conditions without experimental errors and with full control of forcing functions (initial and boundary conditions as well as source/sink terms). Parallel to the generation of synthetic data, we conduct hydraulic tests at multiple scales in a laboratory aquifer with deterministic heterogeneity to generate data that are used to validate K tomograms from hydraulic tomography. Measurements include multiple K estimates from core, slug, single-hole and cross-hole tests as well as several unidirectional, flow-through experiments conducted on the sandbox under steady-state conditions. Validation of K tomograms involved a multi-method and multiscale approach proposed herein which include: (1) visual comparisons of K tomograms to the true sand distributions and the reference K tomogram; (2) testing the ability of 0022-1694/$ -see front matter Published by Elsevier B.V. a v a i l a b l e a t w w w . s c i e n c e d i r e c t . c o m j o u r n a l h o m e p a g e : w w w . e l s e v i e r . c o m / l o c a t e / j h y d r o l
Water Resources Research, 2009
1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m i... more 1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m in deep boreholes at the Mizunami Underground Research Laboratory (MIU) construction site in central Japan. During these two tests, induced groundwater responses were monitored at many observation intervals at various depths in different boreholes at the site. We analyze the two cross-hole pumping tests using transient hydraulic tomography (THT) based on an efficient sequential successive linear estimator to compute the hydraulic conductivity (K) and specific storage (S s ) tomograms, as well as their uncertainties in three dimensions. The equivalent K and S s estimates obtained using asymptotic analysis treating the medium to be homogeneous served as the mean parameter estimates for the 3-D stochastic inverse modeling effort. Results show several, distinct, high K and low S s zones that are continuous over hundreds of meters, which appear to delineate fault zones and their connectivity. The THT analysis of the tests also identified a low K zone which corresponds with a known fault zone trending NNW and has been found to compartmentalize groundwater flow at the site. These results corroborate well with observed water level records, available fault information, and coseismic groundwater level responses during several large earthquakes. The successful application of THT to cross-hole pumping tests conducted in fractured granite at this site suggests that THT is a promising approach to delineate large-scale K and S s heterogeneities, fracture connectivity, and to quantify uncertainty of the estimated fields.
Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hy... more Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hydraulic heterogeneity. To date, a comprehensive validation of hydraulic tomography has not been done either at the laboratory or field scales. The main objective of this paper is to examine the accuracy of hydraulic conductivity (K) tomograms obtained from the steady-state hydraulic tomography algorithm of [Yeh, T.-C. J., Liu, S., 2000. Hydraulic tomography: development of a new aquifer test method. Water Resources Research 36, 2095-2105]. We first obtain a reference K tomogram through the inversion of synthetic cross-hole test data generated through numerical simulations. The purpose of reference K tomogram generation is to examine the ability of the algorithm to image the heterogeneity pattern under optimal conditions without experimental errors and with full control of forcing functions (initial and boundary conditions as well as source/sink terms). Parallel to the generation of synthetic data, we conduct hydraulic tests at multiple scales in a laboratory aquifer with deterministic heterogeneity to generate data that are used to validate K tomograms from hydraulic tomography. Measurements include multiple K estimates from core, slug, single-hole and cross-hole tests as well as several unidirectional, flow-through experiments conducted on the sandbox under steady-state conditions. Validation of K tomograms involved a multi-method and multiscale approach proposed herein which include: (1) visual comparisons of K tomograms to the true sand distributions and the reference K tomogram; (2) testing the ability of 0022-1694/$ -see front matter Published by Elsevier B.V. a v a i l a b l e a t w w w . s c i e n c e d i r e c t . c o m j o u r n a l h o m e p a g e : w w w . e l s e v i e r . c o m / l o c a t e / j h y d r o l
Ground Water, 2007
Hydraulic tomography has been developed as an alternative to traditional geostatistical methods t... more Hydraulic tomography has been developed as an alternative to traditional geostatistical methods to delineate heterogeneity patterns in parameters such as hydraulic conductivity (K) and specific storage (Ss). During hydraulic tomography surveys, a large number of hydraulic head data are collected from a series of cross-hole tests in the subsurface. These head data are then used to interpret the spatial distribution of K and Ss using inverse modeling. Here, we use the Sequential Successive Linear Estimator (SSLE) of Yeh and Liu (2000) to interpret synthetic pumping test data created through numerical simulations and real data generated in a laboratory sandbox aquifer to obtain the K tomograms. Here, we define “K tomogram” as an image of K distribution of the subsurface (or the inverse results) obtained via hydraulic tomography. We examine the influence of signal-to-noise ratio and biases on results using inverse modeling of synthetic and real cross-hole pumping test data. To accomplish this, we first show that the pumping rate, which affects the signal-to-noise ratio, and the order of data included into the SSLE algorithm both have large impacts on the quality of the K tomograms. We then examine the role of conditioning on the K tomogram and find that conditioning can improve the quality of the K tomogram, but can also impair it, if the data are of poor quality and conditioning data have a larger support volume than the numerical grid used to conduct the inversion. Overall, these results show that the quality of the K tomogram depends on the design of pumping tests, their conduct, the order in which they are included in the inverse code, and the quality as well as the support volume of additional data that are used in its computation.
Water Resources Research, 2009
1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m i... more 1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m in deep boreholes at the Mizunami Underground Research Laboratory (MIU) construction site in central Japan. During these two tests, induced groundwater responses were monitored at many observation intervals at various depths in different boreholes at the site. We analyze the two cross-hole pumping tests using transient hydraulic tomography (THT) based on an efficient sequential successive linear estimator to compute the hydraulic conductivity (K) and specific storage (S s ) tomograms, as well as their uncertainties in three dimensions. The equivalent K and S s estimates obtained using asymptotic analysis treating the medium to be homogeneous served as the mean parameter estimates for the 3-D stochastic inverse modeling effort. Results show several, distinct, high K and low S s zones that are continuous over hundreds of meters, which appear to delineate fault zones and their connectivity. The THT analysis of the tests also identified a low K zone which corresponds with a known fault zone trending NNW and has been found to compartmentalize groundwater flow at the site. These results corroborate well with observed water level records, available fault information, and coseismic groundwater level responses during several large earthquakes. The successful application of THT to cross-hole pumping tests conducted in fractured granite at this site suggests that THT is a promising approach to delineate large-scale K and S s heterogeneities, fracture connectivity, and to quantify uncertainty of the estimated fields.
Journal of Hydrology, 2007
Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hy... more Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hydraulic heterogeneity. To date, a comprehensive validation of hydraulic tomography has not been done either at the laboratory or field scales. The main objective of this paper is to examine the accuracy of hydraulic conductivity (K) tomograms obtained from the steady-state hydraulic tomography algorithm of [Yeh, T.-C. J., Liu, S., 2000. Hydraulic tomography: development of a new aquifer test method. Water Resources Research 36, 2095-2105]. We first obtain a reference K tomogram through the inversion of synthetic cross-hole test data generated through numerical simulations. The purpose of reference K tomogram generation is to examine the ability of the algorithm to image the heterogeneity pattern under optimal conditions without experimental errors and with full control of forcing functions (initial and boundary conditions as well as source/sink terms). Parallel to the generation of synthetic data, we conduct hydraulic tests at multiple scales in a laboratory aquifer with deterministic heterogeneity to generate data that are used to validate K tomograms from hydraulic tomography. Measurements include multiple K estimates from core, slug, single-hole and cross-hole tests as well as several unidirectional, flow-through experiments conducted on the sandbox under steady-state conditions. Validation of K tomograms involved a multi-method and multiscale approach proposed herein which include: (1) visual comparisons of K tomograms to the true sand distributions and the reference K tomogram; (2) testing the ability of 0022-1694/$ -see front matter Published by Elsevier B.V. a v a i l a b l e a t w w w . s c i e n c e d i r e c t . c o m j o u r n a l h o m e p a g e : w w w . e l s e v i e r . c o m / l o c a t e / j h y d r o l
Water Resources Research, 2009
1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m i... more 1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m in deep boreholes at the Mizunami Underground Research Laboratory (MIU) construction site in central Japan. During these two tests, induced groundwater responses were monitored at many observation intervals at various depths in different boreholes at the site. We analyze the two cross-hole pumping tests using transient hydraulic tomography (THT) based on an efficient sequential successive linear estimator to compute the hydraulic conductivity (K) and specific storage (S s ) tomograms, as well as their uncertainties in three dimensions. The equivalent K and S s estimates obtained using asymptotic analysis treating the medium to be homogeneous served as the mean parameter estimates for the 3-D stochastic inverse modeling effort. Results show several, distinct, high K and low S s zones that are continuous over hundreds of meters, which appear to delineate fault zones and their connectivity. The THT analysis of the tests also identified a low K zone which corresponds with a known fault zone trending NNW and has been found to compartmentalize groundwater flow at the site. These results corroborate well with observed water level records, available fault information, and coseismic groundwater level responses during several large earthquakes. The successful application of THT to cross-hole pumping tests conducted in fractured granite at this site suggests that THT is a promising approach to delineate large-scale K and S s heterogeneities, fracture connectivity, and to quantify uncertainty of the estimated fields.
Water Resources Research, 2009
1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m i... more 1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m in deep boreholes at the Mizunami Underground Research Laboratory (MIU) construction site in central Japan. During these two tests, induced groundwater responses were monitored at many observation intervals at various depths in different boreholes at the site. We analyze the two cross-hole pumping tests using transient hydraulic tomography (THT) based on an efficient sequential successive linear estimator to compute the hydraulic conductivity (K) and specific storage (S s ) tomograms, as well as their uncertainties in three dimensions. The equivalent K and S s estimates obtained using asymptotic analysis treating the medium to be homogeneous served as the mean parameter estimates for the 3-D stochastic inverse modeling effort. Results show several, distinct, high K and low S s zones that are continuous over hundreds of meters, which appear to delineate fault zones and their connectivity. The THT analysis of the tests also identified a low K zone which corresponds with a known fault zone trending NNW and has been found to compartmentalize groundwater flow at the site. These results corroborate well with observed water level records, available fault information, and coseismic groundwater level responses during several large earthquakes. The successful application of THT to cross-hole pumping tests conducted in fractured granite at this site suggests that THT is a promising approach to delineate large-scale K and S s heterogeneities, fracture connectivity, and to quantify uncertainty of the estimated fields.
Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hy... more Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hydraulic heterogeneity. To date, a comprehensive validation of hydraulic tomography has not been done either at the laboratory or field scales. The main objective of this paper is to examine the accuracy of hydraulic conductivity (K) tomograms obtained from the steady-state hydraulic tomography algorithm of [Yeh, T.-C. J., Liu, S., 2000. Hydraulic tomography: development of a new aquifer test method. Water Resources Research 36, 2095-2105]. We first obtain a reference K tomogram through the inversion of synthetic cross-hole test data generated through numerical simulations. The purpose of reference K tomogram generation is to examine the ability of the algorithm to image the heterogeneity pattern under optimal conditions without experimental errors and with full control of forcing functions (initial and boundary conditions as well as source/sink terms). Parallel to the generation of synthetic data, we conduct hydraulic tests at multiple scales in a laboratory aquifer with deterministic heterogeneity to generate data that are used to validate K tomograms from hydraulic tomography. Measurements include multiple K estimates from core, slug, single-hole and cross-hole tests as well as several unidirectional, flow-through experiments conducted on the sandbox under steady-state conditions. Validation of K tomograms involved a multi-method and multiscale approach proposed herein which include: (1) visual comparisons of K tomograms to the true sand distributions and the reference K tomogram; (2) testing the ability of 0022-1694/$ -see front matter Published by Elsevier B.V. a v a i l a b l e a t w w w . s c i e n c e d i r e c t . c o m j o u r n a l h o m e p a g e : w w w . e l s e v i e r . c o m / l o c a t e / j h y d r o l
Journal of Hydrology, 2007
Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hy... more Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hydraulic heterogeneity. To date, a comprehensive validation of hydraulic tomography has not been done either at the laboratory or field scales. The main objective of this paper is to examine the accuracy of hydraulic conductivity (K) tomograms obtained from the steady-state hydraulic tomography algorithm of [Yeh, T.-C. J., Liu, S., 2000. Hydraulic tomography: development of a new aquifer test method. Water Resources Research 36, 2095-2105]. We first obtain a reference K tomogram through the inversion of synthetic cross-hole test data generated through numerical simulations. The purpose of reference K tomogram generation is to examine the ability of the algorithm to image the heterogeneity pattern under optimal conditions without experimental errors and with full control of forcing functions (initial and boundary conditions as well as source/sink terms). Parallel to the generation of synthetic data, we conduct hydraulic tests at multiple scales in a laboratory aquifer with deterministic heterogeneity to generate data that are used to validate K tomograms from hydraulic tomography. Measurements include multiple K estimates from core, slug, single-hole and cross-hole tests as well as several unidirectional, flow-through experiments conducted on the sandbox under steady-state conditions. Validation of K tomograms involved a multi-method and multiscale approach proposed herein which include: (1) visual comparisons of K tomograms to the true sand distributions and the reference K tomogram; (2) testing the ability of 0022-1694/$ -see front matter Published by Elsevier B.V. a v a i l a b l e a t w w w . s c i e n c e d i r e c t . c o m j o u r n a l h o m e p a g e : w w w . e l s e v i e r . c o m / l o c a t e / j h y d r o l
Water Resources Research, 2009
1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m i... more 1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m in deep boreholes at the Mizunami Underground Research Laboratory (MIU) construction site in central Japan. During these two tests, induced groundwater responses were monitored at many observation intervals at various depths in different boreholes at the site. We analyze the two cross-hole pumping tests using transient hydraulic tomography (THT) based on an efficient sequential successive linear estimator to compute the hydraulic conductivity (K) and specific storage (S s ) tomograms, as well as their uncertainties in three dimensions. The equivalent K and S s estimates obtained using asymptotic analysis treating the medium to be homogeneous served as the mean parameter estimates for the 3-D stochastic inverse modeling effort. Results show several, distinct, high K and low S s zones that are continuous over hundreds of meters, which appear to delineate fault zones and their connectivity. The THT analysis of the tests also identified a low K zone which corresponds with a known fault zone trending NNW and has been found to compartmentalize groundwater flow at the site. These results corroborate well with observed water level records, available fault information, and coseismic groundwater level responses during several large earthquakes. The successful application of THT to cross-hole pumping tests conducted in fractured granite at this site suggests that THT is a promising approach to delineate large-scale K and S s heterogeneities, fracture connectivity, and to quantify uncertainty of the estimated fields.
Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hy... more Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hydraulic heterogeneity. To date, a comprehensive validation of hydraulic tomography has not been done either at the laboratory or field scales. The main objective of this paper is to examine the accuracy of hydraulic conductivity (K) tomograms obtained from the steady-state hydraulic tomography algorithm of [Yeh, T.-C. J., Liu, S., 2000. Hydraulic tomography: development of a new aquifer test method. Water Resources Research 36, 2095-2105]. We first obtain a reference K tomogram through the inversion of synthetic cross-hole test data generated through numerical simulations. The purpose of reference K tomogram generation is to examine the ability of the algorithm to image the heterogeneity pattern under optimal conditions without experimental errors and with full control of forcing functions (initial and boundary conditions as well as source/sink terms). Parallel to the generation of synthetic data, we conduct hydraulic tests at multiple scales in a laboratory aquifer with deterministic heterogeneity to generate data that are used to validate K tomograms from hydraulic tomography. Measurements include multiple K estimates from core, slug, single-hole and cross-hole tests as well as several unidirectional, flow-through experiments conducted on the sandbox under steady-state conditions. Validation of K tomograms involved a multi-method and multiscale approach proposed herein which include: (1) visual comparisons of K tomograms to the true sand distributions and the reference K tomogram; (2) testing the ability of 0022-1694/$ -see front matter Published by Elsevier B.V. a v a i l a b l e a t w w w . s c i e n c e d i r e c t . c o m j o u r n a l h o m e p a g e : w w w . e l s e v i e r . c o m / l o c a t e / j h y d r o l
Ground Water, 2007
Hydraulic tomography has been developed as an alternative to traditional geostatistical methods t... more Hydraulic tomography has been developed as an alternative to traditional geostatistical methods to delineate heterogeneity patterns in parameters such as hydraulic conductivity (K) and specific storage (Ss). During hydraulic tomography surveys, a large number of hydraulic head data are collected from a series of cross-hole tests in the subsurface. These head data are then used to interpret the spatial distribution of K and Ss using inverse modeling. Here, we use the Sequential Successive Linear Estimator (SSLE) of Yeh and Liu (2000) to interpret synthetic pumping test data created through numerical simulations and real data generated in a laboratory sandbox aquifer to obtain the K tomograms. Here, we define “K tomogram” as an image of K distribution of the subsurface (or the inverse results) obtained via hydraulic tomography. We examine the influence of signal-to-noise ratio and biases on results using inverse modeling of synthetic and real cross-hole pumping test data. To accomplish this, we first show that the pumping rate, which affects the signal-to-noise ratio, and the order of data included into the SSLE algorithm both have large impacts on the quality of the K tomograms. We then examine the role of conditioning on the K tomogram and find that conditioning can improve the quality of the K tomogram, but can also impair it, if the data are of poor quality and conditioning data have a larger support volume than the numerical grid used to conduct the inversion. Overall, these results show that the quality of the K tomogram depends on the design of pumping tests, their conduct, the order in which they are included in the inverse code, and the quality as well as the support volume of additional data that are used in its computation.
Ground Water, 2007
Hydraulic tomography has been developed as an alternative to traditional geostatistical methods t... more Hydraulic tomography has been developed as an alternative to traditional geostatistical methods to delineate heterogeneity patterns in parameters such as hydraulic conductivity (K) and specific storage (Ss). During hydraulic tomography surveys, a large number of hydraulic head data are collected from a series of cross-hole tests in the subsurface. These head data are then used to interpret the spatial distribution of K and Ss using inverse modeling. Here, we use the Sequential Successive Linear Estimator (SSLE) of Yeh and Liu (2000) to interpret synthetic pumping test data created through numerical simulations and real data generated in a laboratory sandbox aquifer to obtain the K tomograms. Here, we define “K tomogram” as an image of K distribution of the subsurface (or the inverse results) obtained via hydraulic tomography. We examine the influence of signal-to-noise ratio and biases on results using inverse modeling of synthetic and real cross-hole pumping test data. To accomplish this, we first show that the pumping rate, which affects the signal-to-noise ratio, and the order of data included into the SSLE algorithm both have large impacts on the quality of the K tomograms. We then examine the role of conditioning on the K tomogram and find that conditioning can improve the quality of the K tomogram, but can also impair it, if the data are of poor quality and conditioning data have a larger support volume than the numerical grid used to conduct the inversion. Overall, these results show that the quality of the K tomogram depends on the design of pumping tests, their conduct, the order in which they are included in the inverse code, and the quality as well as the support volume of additional data that are used in its computation.
Water Resources Research, 2009
1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m i... more 1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m in deep boreholes at the Mizunami Underground Research Laboratory (MIU) construction site in central Japan. During these two tests, induced groundwater responses were monitored at many observation intervals at various depths in different boreholes at the site. We analyze the two cross-hole pumping tests using transient hydraulic tomography (THT) based on an efficient sequential successive linear estimator to compute the hydraulic conductivity (K) and specific storage (S s ) tomograms, as well as their uncertainties in three dimensions. The equivalent K and S s estimates obtained using asymptotic analysis treating the medium to be homogeneous served as the mean parameter estimates for the 3-D stochastic inverse modeling effort. Results show several, distinct, high K and low S s zones that are continuous over hundreds of meters, which appear to delineate fault zones and their connectivity. The THT analysis of the tests also identified a low K zone which corresponds with a known fault zone trending NNW and has been found to compartmentalize groundwater flow at the site. These results corroborate well with observed water level records, available fault information, and coseismic groundwater level responses during several large earthquakes. The successful application of THT to cross-hole pumping tests conducted in fractured granite at this site suggests that THT is a promising approach to delineate large-scale K and S s heterogeneities, fracture connectivity, and to quantify uncertainty of the estimated fields.
Journal of Hydrology, 2007
Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hy... more Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hydraulic heterogeneity. To date, a comprehensive validation of hydraulic tomography has not been done either at the laboratory or field scales. The main objective of this paper is to examine the accuracy of hydraulic conductivity (K) tomograms obtained from the steady-state hydraulic tomography algorithm of [Yeh, T.-C. J., Liu, S., 2000. Hydraulic tomography: development of a new aquifer test method. Water Resources Research 36, 2095-2105]. We first obtain a reference K tomogram through the inversion of synthetic cross-hole test data generated through numerical simulations. The purpose of reference K tomogram generation is to examine the ability of the algorithm to image the heterogeneity pattern under optimal conditions without experimental errors and with full control of forcing functions (initial and boundary conditions as well as source/sink terms). Parallel to the generation of synthetic data, we conduct hydraulic tests at multiple scales in a laboratory aquifer with deterministic heterogeneity to generate data that are used to validate K tomograms from hydraulic tomography. Measurements include multiple K estimates from core, slug, single-hole and cross-hole tests as well as several unidirectional, flow-through experiments conducted on the sandbox under steady-state conditions. Validation of K tomograms involved a multi-method and multiscale approach proposed herein which include: (1) visual comparisons of K tomograms to the true sand distributions and the reference K tomogram; (2) testing the ability of 0022-1694/$ -see front matter Published by Elsevier B.V. a v a i l a b l e a t w w w . s c i e n c e d i r e c t . c o m j o u r n a l h o m e p a g e : w w w . e l s e v i e r . c o m / l o c a t e / j h y d r o l
Water Resources Research, 2009
1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m i... more 1] Two large-scale cross-hole pumping tests were conducted at depths of 191-226 m and 662-706 m in deep boreholes at the Mizunami Underground Research Laboratory (MIU) construction site in central Japan. During these two tests, induced groundwater responses were monitored at many observation intervals at various depths in different boreholes at the site. We analyze the two cross-hole pumping tests using transient hydraulic tomography (THT) based on an efficient sequential successive linear estimator to compute the hydraulic conductivity (K) and specific storage (S s ) tomograms, as well as their uncertainties in three dimensions. The equivalent K and S s estimates obtained using asymptotic analysis treating the medium to be homogeneous served as the mean parameter estimates for the 3-D stochastic inverse modeling effort. Results show several, distinct, high K and low S s zones that are continuous over hundreds of meters, which appear to delineate fault zones and their connectivity. The THT analysis of the tests also identified a low K zone which corresponds with a known fault zone trending NNW and has been found to compartmentalize groundwater flow at the site. These results corroborate well with observed water level records, available fault information, and coseismic groundwater level responses during several large earthquakes. The successful application of THT to cross-hole pumping tests conducted in fractured granite at this site suggests that THT is a promising approach to delineate large-scale K and S s heterogeneities, fracture connectivity, and to quantify uncertainty of the estimated fields.
Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hy... more Hydraulic tomography potentially is a viable technology that facilitates subsurface imaging of hydraulic heterogeneity. To date, a comprehensive validation of hydraulic tomography has not been done either at the laboratory or field scales. The main objective of this paper is to examine the accuracy of hydraulic conductivity (K) tomograms obtained from the steady-state hydraulic tomography algorithm of [Yeh, T.-C. J., Liu, S., 2000. Hydraulic tomography: development of a new aquifer test method. Water Resources Research 36, 2095-2105]. We first obtain a reference K tomogram through the inversion of synthetic cross-hole test data generated through numerical simulations. The purpose of reference K tomogram generation is to examine the ability of the algorithm to image the heterogeneity pattern under optimal conditions without experimental errors and with full control of forcing functions (initial and boundary conditions as well as source/sink terms). Parallel to the generation of synthetic data, we conduct hydraulic tests at multiple scales in a laboratory aquifer with deterministic heterogeneity to generate data that are used to validate K tomograms from hydraulic tomography. Measurements include multiple K estimates from core, slug, single-hole and cross-hole tests as well as several unidirectional, flow-through experiments conducted on the sandbox under steady-state conditions. Validation of K tomograms involved a multi-method and multiscale approach proposed herein which include: (1) visual comparisons of K tomograms to the true sand distributions and the reference K tomogram; (2) testing the ability of 0022-1694/$ -see front matter Published by Elsevier B.V. a v a i l a b l e a t w w w . s c i e n c e d i r e c t . c o m j o u r n a l h o m e p a g e : w w w . e l s e v i e r . c o m / l o c a t e / j h y d r o l
Ground Water, 2007
Hydraulic tomography has been developed as an alternative to traditional geostatistical methods t... more Hydraulic tomography has been developed as an alternative to traditional geostatistical methods to delineate heterogeneity patterns in parameters such as hydraulic conductivity (K) and specific storage (Ss). During hydraulic tomography surveys, a large number of hydraulic head data are collected from a series of cross-hole tests in the subsurface. These head data are then used to interpret the spatial distribution of K and Ss using inverse modeling. Here, we use the Sequential Successive Linear Estimator (SSLE) of Yeh and Liu (2000) to interpret synthetic pumping test data created through numerical simulations and real data generated in a laboratory sandbox aquifer to obtain the K tomograms. Here, we define “K tomogram” as an image of K distribution of the subsurface (or the inverse results) obtained via hydraulic tomography. We examine the influence of signal-to-noise ratio and biases on results using inverse modeling of synthetic and real cross-hole pumping test data. To accomplish this, we first show that the pumping rate, which affects the signal-to-noise ratio, and the order of data included into the SSLE algorithm both have large impacts on the quality of the K tomograms. We then examine the role of conditioning on the K tomogram and find that conditioning can improve the quality of the K tomogram, but can also impair it, if the data are of poor quality and conditioning data have a larger support volume than the numerical grid used to conduct the inversion. Overall, these results show that the quality of the K tomogram depends on the design of pumping tests, their conduct, the order in which they are included in the inverse code, and the quality as well as the support volume of additional data that are used in its computation.