Recent work on two-phase free boundary problems has led to the investigation of a new type of quadrature domain for harmonic functions. This paper develops a method of constructing such quadrature domains based on the technique of partial balayage, which has proved to be a useful tool in the study of one-phase quadrature domains and Hele-Shaw flows.