The results of numerical simulations of hydraulic fracturing based on the concept of lattice, a grain based method, are presented. Five cases of numerical simulations were conducted to analyze the effects of changing rock properties and in-situ stresses on hydraulic fracture propagation and its containment. The study was done based on the data from one of the shale gas wells in the North Perth Basin (NPB) in Western Australia and the hydraulic fracturing lab experimental data on a 50 mm cubical sample. The hydraulic fracturing initiates and propagates within Carynginia Formation surrounded by the Irwin River Coal Measures (IRCM) Formation from the top and Kockatea Formation at the base. The lab experiments were upscaled to both intermediate and field scales to simulate fracture propagation in the viscosity dominated regime. Penetration parameters, including the depth and area of the fracture crossing the interface, were proposed to characterize fracture containment capacity. It was observed that the fracture is contained by both IRCM and Kockatea Formations due to the large contrast between the minimum horizontal stresses in these layers and that of the Carynginia formation. However, as quantified by the penetration parameters, the fracture propagates more easily into the IRCM Formation with higher brittleness than Carynginia Formation. Kockatea Formation with lower brittleness than that of Carynginia Formation inhibits vertical migration of the fracture, i.e. results in more visible fracture containment. The results suggest that decrease in rock brittleness inhibits fracture propagation and lower minimum horizontal stress contrast leads to fracture containment.