Просто разбейте конкатенацию, как предложил @brezniczky. Или используйте hstack:
import numpy as np
a = np.ones((508,12))
b = np.zeros((508,13))
col = 3
final = np.hstack((b[:,0:col],a,b[:,col:]))
print(final[0])
[0. 0. 0. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.]
В этом сценарии hstack и конкатенация с осью = 1 одинаковы, я просто предпочитаю hstack для лучшей читаемости