In this talk, first, we show that the color N-dependent area law falloffs of the double-winding Wilson loop averages for the SU(N) lattice gauge model are reproduced from the Z(N) lattice Abelian gauge model due to the center group dominance in quark confinement. Next, we discuss lattice
gauge-scalar models which allow analytic continuation for gauge invariant operators between confinement region and Higgs region. Applying the cluster expansion, we try to understand nontrivial contribution from scalar field in quark confinement mechanism. In order to understand quark confinement further, moreover, we study double-winding Wilson loop averages in the analytical region on the phase diagram.