Flow to wells in nonuniform geologic formations is of central interest to hydrogeologists and petroleum engineers. There are, however, very few mathematical analyses of such flow. We present analytical expressions for leading statistical moments of vertically averaged hydraulic head and flux under steady-state flow to a well that pumps water from a bounded, randomly heterogeneous aquifer. Like in the widely used Thiem equation, we prescribe a constant pumping rate deterministically at the well and a constant head at a circular outer boundary of radius L. We model the natural logarithm Y = 1nT of aquifer transmissivity T as a statistically homogeneous random field with a Gaussian spatial correlation function. Our solution is based on exact nonlocal moment equations for multidimensional steady state flow in bounded, randomly heterogeneous porous media. Perturbation of these nonlocal equations leads to a system of local recursive moment equations that we solve analytically to second order in the standard deviation of Y. In contrast to most stochastic analyses of flow, which require that log transmissivity be multivariate Gaussian, our solution is free of any distributional assumptions. It yields expected values of head and flux, and the variance-covariance of these quantities, as functions of distance from the well. It also yields an apparent transmissivity, T-a, defined as the negative ratio between expected flux and head gradient at any radial distance. The solution is supported by numerical Monte Carlo simulations, which demonstrate that it is applicable to strongly heterogeneous aquifers, characterized by large values of log transmissivity variance. The two-dimensional nature of our solution renders it useful for relatively thin aquifers in which vertical heterogeneity tends to be of minor concern relative to that in the horizontal plane. It also applies to thicker aquifers when information about their vertical heterogeneity is lacking, as is commonly the case when measurements of head and flow rate are done in wells that penetrate much of the aquifer thickness. Potential uses include the analysis of pumping tests and tracer test conducted in such wells, the statistical delineation of their respective capture zones, and the analysis of contaminant transport toward fully penetrating wells.