We consider the following repulsive-productive chemotaxis model: Let p ∈ (1, 2), find u ≥ 0, the cell density, and v ≥ 0, the chemical concentration, satisfyingin a bounded domain Ω ⊆ R d , d = 2, 3. By using a regularization technique, we prove the existence of solutions for problem (1). Moreover, we propose three fully discrete Finite Element (FE) nonlinear approximations of problem (1), where the first one is defined in the variables (u, v), and the second and third ones introduce σ = ∇v as auxiliary variable.We prove some unconditional properties such as mass-conservation, energy-stability and solvability of the schemes. Finally, we compare the behavior of these schemes throughout several numerical simulations and give some conclusions.