Abstract: Yangjiaopu protective district located in Fengfeng mining site is a large water source district. About90 thousand m3 of karst water from Ordovician is explored in Yangjiaopu every day.“Sub-divisional Technological Regulation for Source Protective Zones of Drinking Water(H J/T338-200)”published by State Environmental Protection Administration in2007 states that division of the large scale water protective zone should use numerical model to simulate and calculate. According this regulation, this paper simulate and forecast the solute transport of Yangjiaopu karst groundwater system using numerical simulation method with visual modflow software, to divide the extent of the water source protective zones in different grade. The division result is as follow: The area centered on the groundwater exploitation well and delimited by the radius of solutes transporting for 100 day, can be regarded as the first-grade protection zone. Excepting the first-grade protection zone,the area delimited by the radius of solutestransportfor1000days,is consider as the second-grade protection zone. The bare Ordovices Limestone area recharging the karst water subsystem could be divided into prospective protection zone.