While two-dimensional symmetry-enriched topological phases (π²π€π³s) have been studied intensively and systematically, three-dimensional ones are still open issues. We propose an algorithmic approach of imposing global symmetry Gs on gauge theories (denoted by π¦π³) with gauge group Gg. The resulting symmetric gauge theories are dubbed "symmetry-enriched gauge theories" (π²π€π¦), which may be served as low-energy effective theories of three-dimensional symmetric topological quantum spin liquids. We focus on π²π€π¦s with gauge group Gg=β€N1Γβ€N2Γβ― and on-site unitary symmetry group Gs=β€K1Γβ€K2Γβ― or Gs=U(1)Γβ€K1Γβ―. Each π²π€π¦(Gg,Gs) is described in the path integral formalism associated with certain symmetry assignment. From the path-integral expression, we propose how to physically diagnose the ground state properties (i.e., π²π€π³ orders) of π²π€π¦s in experiments of charge-loop braidings (patterns of symmetry fractionalization) and the \emph{mixed} multi-loop braidings among deconfined loop excitations and confined symmetry fluxes. From these symmetry-enriched properties, one can obtain the map from π²π€π¦s to π²π€π³s. By giving full dynamics to background gauge fields, π²π€π¦s may be eventually promoted to a set of new gauge theories (denoted by π¦π³β). Based on their gauge groups, π¦π³βs may be further regrouped into different classes each of which is labeled by a gauge group Gβg. Finally, a web of gauge theories involving π¦π³, π²π€π¦, π²π€π³ and π¦π³β is achieved. We demonstrate the above symmetry-enrichment physics and the web of gauge theories through many concrete examples.