<p>This research presented the direction and distribution of the in situ stress field along three metallogenic fault zones in the northwestern Jiaodong Peninsula of China by analyzing an extensive dataset of 164 in situ stress measurements from six major gold mines. A geological modeling and numerical simulation coupling technology was developed, along with an interface program and grid optimization techniques to facilitate the transition from geological models to numerical computation models. Employing the BP neural network model coupled with inversion methods, the study has successfully derived the optimal boundary conditions for simulating the stress field around the fault zones. A composite metallogenic fault zone comprising two faults F1 and F3 in the Sanshandao gold mine was selected to conduct inversion simulation. These findings reveal that the tectonic stress controls the evolution process of metallogenic faults, and the fault evolution promotes the formation of the present-day stress field. The maximum horizontal stress (σ_H) in the northwestern Jiaodong Peninsula is predominantly oriented in an NWW-SEE direction, which is nearly perpendicular to the strike of the faults. The interaction between the in situ stress field and the composite fault zone demonstrates that the presence of a fault results in a discontinuous transfer of stress across its boundaries. The maximum principal stress (σ_max) at the intersection of two faults and the surrounding area is lower than that at equivalent depths elsewhere. When a fault is horizontal with broken rock, stress release occurs within the fault, leading to stress concentration along its sides. Additionally, the orientation and shape of a fault significantly influence the magnitude and direction of stress deflection around it. Variations in these characteristics can substantially alter the distribution of the crustal stress field. Conversely, the stress field may influence the shape of the fault.</p>