  Simulating 2D flow in fractured crystalline rock requires 2D stochastic
discrete-fracture matrix (DFM) models. To obtain the simulation statistics of
interest at an affordable computational cost, we aim to use the multilevel
Monte Carlo method. To use this multiscale approach, one needs to upscale the
hydraulic conductivity of the fractures by numerical homogenization. In this
work, we substitute numerical homogenization with a surrogate model to speed up
the computations. In particular, we resort to a deep convolutional neural
network (CNN) connected to a deep feed-forward neural network. The equivalent
hydraulic conductivity tensor $K_{eq}$ is predicted based on an input spatial
random field (SRF) of hydraulic conductivity tensors, cross-section, and
hydraulic conductivity of fractures. Three independent surrogates with the same
architecture are trained using data from DFM models with three different ratios
of hydraulic conductivities of fracture and bulk $K_f/K_b$. As the $K_f/K_b$
ratio increases, the multivariate $K_{eq}$ distribution becomes more complex,
and thus, the prediction accuracy of the trained surrogates deteriorates.
Regardless of $K_f/K_b$, however, an improvement in the prediction accuracy of
the trained surrogates is noted as the considered fracture density of the
modeling setup decreases. We also investigate prediction accuracy on input SRFs
of different correlation lengths. Upscaling by numerical homogenization and by
surrogate modeling is compared on two practical problems: upscaling of the
hydraulic conductivity tensor and groundwater flow through a given surface. We
obtained equally accurate results for the equivalent hydraulic tensor
calculation of upscaled DFM models regardless of the upscaling method. For the
groundwater flow problem, the accuracy of quantity of interest imitates the
accuracy of $K_{eq}$ predictions.
