Balance-equation method for simulating terahertz quantum-cascade lasers using a wave-function basis with reduced dipole moments of tunnel-coupled states