<p>The Zhugongtang lead-zinc deposit is a super-large lead-zinc deposit controlled by the fault-fold structure in the Western Guizhou, China. In this paper, the complex migration process of ore-bearing fluid in fault-fold structures is reproduced by numerical simulation based on a multi-field coupled framework combining thermodynamics and fluid mechanics. The results show that ore-bearing fluid first migrates vertically rapidly along the fault and then transfers horizontally along the stratum to the SE flank of the fold. Negative pressure anomalies appear at the narrowing position in the deep part of the fault and at the intersection of faults and folds. The spatial distribution of Zn²⁺ is related to the known mineralization distribution. Under the influence of heat flow field and pressure, ore-bearing fluid tends to the abnormal area of negative stress. Comparative experiments with different interlimb angle show that for smaller interlimb angle, mineralization tends to occur in faults, and conversely, mineralization tends to occur in strata. This is consistent with the structural control of mineralization observed in other deposits in northwestern Guizhou. This research provides a new insight into the understanding of fault ore control model.</p>