An incremental variational gradient damage model for saturated poroelastic media incorporating cohesive zone effects