[1] Periodic pumping tests, in which a fluid is extracted during half a period, then reinjected, have been used historically to estimate effective aquifer properties. In this work, we suggest a modified approach to periodic pumping test analysis in which one uses several periodic pumping signals of different frequencies as stimulation, and responses are analyzed through inverse modeling using a ''steady-periodic'' model formulation. We refer to this strategy as multifrequency oscillatory hydraulic imaging. Oscillating pumping tests have several advantages that have been noted, including no net water extraction during testing and robust signal measurement through signal processing. Through numerical experiments, we demonstrate additional distinct advantages that multifrequency stimulations have, including: (1) drastically reduced computational cost through use of a steady-periodic numerical model and (2) full utilization of the aquifer heterogeneity information provided by responses at different frequencies. We first perform fully transient numerical modeling for heterogeneous aquifers and show that equivalent results are obtained using a faster steady-periodic heterogeneous numerical model of the wave phasor. The sensitivities of observed signal response to aquifer heterogeneities are derived using an adjoint state-based approach, which shows that different frequency stimulations provide complementary information. Finally, we present an example 2-D application in which sinusoidal signals at multiple frequencies are used as a data source and are inverted to obtain estimates of aquifer heterogeneity. These analyses show the different heterogeneity information that can be obtained from different stimulation frequencies, and that data from several sinusoidal pumping tests can be rapidly inverted using the steady-periodic framework.

We discuss efficient solutions to systems of shifted linear systems arising in computations for oscillatory hydraulic tomography (OHT). The reconstruction of hydrogeological parameters such as hydraulic conductivity and specific storage using limited discrete measurements of pressure (head) obtained from sequential oscillatory pumping tests, leads to a nonlinear inverse problem. We tackle this using the quasi-linear geostatistical approach [15]. This method requires repeated solution of the forward (and adjoint) problem for multiple frequencies, for which we use flexible preconditioned Krylov subspace solvers specifically designed for shifted systems based on ideas in [13]. The solvers allow the preconditioner to change at each iteration. We analyze the convergence of the solver and perform an error analysis when an iterative solver is used for inverting the preconditioner matrices. Finally, we apply our algorithm to a challenging application taken from oscillatory hydraulic tomography to demonstrate the computational gains by using the resulting method.

No abstract

Oscillatory pumping tests were conducted at the Boise Hydrological Research Site. A periodic pressure signal is generated by pumping and injecting water into the aquifer consecutively and the pressure response is recorded at many points around the source. We present and analyze the data from the

Abstract. An implementation of GMRES with multiple preconditioners (MPGMRES) is proposed for solving shifted linear systems with shift-and-invert preconditioners. With this type of preconditioner, the Krylov subspace can be built without requiring the matrix-vector product with the shifted matrix. Furthermore, the multipreconditioned search space is shown to grow only linearly with the number of preconditioners. This allows for a more efficient implementation of the algorithm. The proposed implementation is tested on shifted systems that arise in computational hydrology and the evaluation of different matrix functions. The numerical results indicate the effectiveness of the proposed approach.