1) There indeed exist theories in which the gauge vector-like fields coupled to the axial fermion current enter the game. The most familiar example is, of course, the Standard model, in which there is the local $SU_{L}(2)$ symmetry (left fermion doublets interact with 3 gauge fields $W_{\mu}^{a}$). It is possible to rewrite the theory in terms of vector-axial basis,
$$
\psi_{L} \equiv \frac{1}{2}\psi - \frac{\gamma_{5}}{2}\psi,
$$
and corresponding vector-axial gauge fields are
$$
V_{\mu}^{a} = \frac{W_{\mu}^{a}}{2}, \quad A_{\mu}^{a} = -\frac{W_{\mu}^{a}}{2}
$$
2) There also exist many realistic field theories in which the axial-vector fields exist, but where they are not the gauge fields; typically there are fields representing some particles or the symmetries of the theory.
The familiar example is axial-vector mesons in the QCD near and below the global chiral symmetry $SU_{L}(3)\times SU_{R}(3)$ breaking scale. An approach to introduce these mesons is following: one may gauge this symmetry by introducing massless axial-vector fields, and after that to break it explicitly by adding the mass terms. Although this approach looks unnatural, in fact it has some theoretical origin (the action of chiral perturbation theory has hidden local gauge $SU_{L}(3)\times SU_{R}(3)$ symmetry), phenomenological origin (axial-vector mesons is of course the part of the QCD which respects the approximate chiral symmetry) deep historical roots (the so-called vector meson dominance model) and more or less succesfully describes the data.
Another example of such effective theory is the Weyl semimetal near the bands crossing point. It is given by the theory of massless chiral fermions with non-zero distance in momentum and energy space between their spectrum (being the Dirac cones), parametrized by $b_{\mu}$. It is local because of tensions and dislocations into the semimetal. The lagrangian of such model effectively coincides with
$$
L = \bar{\psi}(i\gamma_{\mu}\partial^{\mu} -\gamma^{\mu}\gamma_{5}b_{\mu})\psi
$$
3) Also, background axial gauge fields are sometime introduced when we need to define the axial current through the action. Precisely, if we introduce such field $A_{\mu}$, then the corresponding current is given by
$$
J_{\mu}^{A}(x) = \frac{\delta \Gamma[A]}{\delta A^{\mu}(x)}
$$
After calculating the current and its properties (often related to the chiral anomaly) the coupling to $A$ is set to zero. This trick was used by Bardeen when he calculated the anomaly in $SU_{L}(3)\times SU_{R}(3)$ theory similar to the QCD.
Note about the mixed anomalies
Of course, if You have both vector and axial-vector gauged symmetries, then the gauge anomaly appears. If You have only one fermion specie, then in general the theory will be inconsistent. Other fermions required in order to cancel this anomaly. However, You see from the text written above that often axial (or vector) gauge fields are in fact fictive (or may even correspond to physical particle), and in these cases You have not to worry about the anomalies. Indeed, suppose the theory in which both vector $V_{\mu}$ and axial-vector $A_{\mu}$ fields are present; however, the axial field is physical. Then the (consistent) anomaly for vector and axial currents is
$$
\tag 1 \partial_{\mu}J^{\mu}_{V} = \frac{1}{48\pi^{2}}\epsilon^{\mu\nu\alpha\beta}F_{\mu\nu}^{A}F_{\alpha\beta}^{V},
$$
$$
\tag 2 \partial_{\mu}J^{\mu}_{A} = \frac{1}{96\pi^{2}}\epsilon^{\mu\nu\alpha\beta}(F_{\mu\nu}^{A}F_{\alpha\beta}^{A} + F_{\mu\nu}^{V}F_{\alpha\beta}^{V})
$$
In order to remove the anomaly from the vector current conservation (and corresponding Ward identities, of course) You have to add the local counter-term (in literature it is called the Bardeen counter-term), which cancels the right hand-side of $(1)$. The expression $(2)$ instead takes the form
$$
\partial_{\mu}J^{\mu}_{A} = \frac{1}{96\pi^{2}}\epsilon^{\mu\nu\alpha\beta}(F_{\mu\nu}^{A}F_{\alpha\beta}^{A} + 3F_{\mu\nu}^{V}F_{\alpha\beta}^{V})
$$
This counter-term even generates the additional physical part of the vector current $J_{\mu}$, namely
$$
\Delta J_{\mu}^{V} \sim \epsilon_{\mu\nu\alpha\beta}A^{\nu}F^{\alpha\beta}_{V}
$$