Enhanced geothermal systems (EGSs) are a promising concept to make geothermal power generation available in many regions worldwide. However, in current EGSs generally only a fraction of the geothermal reservoir is effectively accessed by the production fluid, resulting in suboptimal performance. Recent studies in the literature on a two-dimensional Darcy flow in a circular reservoir driven by reoriented injector–producer wells demonstrated that well configurations and pumping schemes designed on the basis of chaos theory enable efficient distribution of the production fluids throughout the entire reservoir. Key to this is accomplishment of chaotic advection, i.e. the rapid dispersion and stretching of material fluid regions, by a “proper” flow forcing. However, these studies concern isotropic reservoirs, while geothermal reservoirs typically are highly anisotropic. Our theoretical/computational study expands on said studies by investigating the impact of anisotropy on fundamental aspects of the Lagrangian transport of production fluids. This reveals that anisotropy generically eliminates key organizing mechanisms in the Lagrangian transport, viz. symmetries, and thus tends to promote disorder and, inherently, chaotic advection. However, symmetries are partially preserved—and thus order and coherence partially restored—in non-generic cases such as pumping schemes employing an even number of injector–producer wells and well configurations aligned with the anisotropy. Symmetry associated with well alignment in fact appears crucial to an intriguing “order within chaos” observed only in such cases: prolonged confinement of fluid to subregions of chaotic areas.