The surrogate model predicts the hydraulic conductivity tensor Keq from voxelized 3D domains representing tensor‑valued random fields of DFM. Replacing expensive 3D numerical homogenization, it cuts simulation time for fractured crystalline media by an order of magnitude. The method integrates into a multilevel Monte Carlo framework, enabling repeated evaluations. Engineers can now perform groundwater flow studies with less computational overhead.