By using the Wilf-Zeilberger method, we prove a novel finite combinatorial identity related to a bivariate generating function for $ζ(2+r+2s)$ (an extension of a Bailey-Borwein-Bradley Apery-like formula for even zeta values). Such identity is then applied to show several supercongruences.