Use Z/(2)[a,b,c,d];

K:=Ideal(a^2+a, b^2+b, c^2+c, d^2+d);

I:=Ideal(ab+a, bc+b, cd+c, bcd, abd+ad+bd+d) + K;

ReducedGBasis(I);