A novel approach to the stability of clathrate hydrates: grandcanonical MC simulation