Here is a workaround:

load(grobner);
[x^2+y^2-d^2, (x-h)^2+(y-k)^2-s^2];
poly_reduced_grobner(%, [x,y]);
algsys(%, [x,y]);

Andrej