Вы можете просто зациклить,
import numpy as np
x1 = 0.; y1=0.; z1 = 1.
x2 = 1.; y2=0.; z2 = 1.
centers=[[x1, y1], [x2, y2]]
stds=[z1, z2]
np.random.seed(1)
smpl = []
for c, std in zip(centers, stds):
smpl.append(np.random.normal(loc=c, scale=std))
print(smpl)
, но передача в виде списков также, кажется, работает и, вероятно, будет более эффективной,
np.random.seed(1)
smpl = np.random.normal(loc=centers, scale=std)
print(smpl)