  We propose a non-perturbative method to determine the mixing coefficients of
$\Delta s=2$ four-quark operators for the Wilson quark action using chiral Ward
identities. The method is applied to calculate B_K in quenched QCD.
