Hydraulic fracturing is an important stimulation method for the low-permeability oil and gas reservoirs，and an increasing number of studies in this aspect have been proposed. However，the heterogeneity of reservoirs has not been considered in the most studies，in spite of that the heterogeneity is one of the key geological factor affecting the fracturing performance. Therefore，in this study，a model was proposed based on the material point method to estimate the stimulation of heterogeneous reservoir，and the J-integral was used to estimate the stimulated reservoir area and to evaluate the heterogeneity. The simulated in-situ differential stress field was compared with the microseismic detection results so as to verify the accuracy of the model. The results show that the calculation results of model considering the heterogeneity highly matched with the actual microseismic signals in the field；a higher J-integral curve fluctuation indicated a stronger heterogeneity of reservoir；15-stage hydraulic fracture stimulation is the most reasonable with the ratio of stimulated reservoir area as 0.87 and the growth rate less than 0.01；the fractures at the two wings of the heterogeneous reservoir were quite different，and they are different from the designed half-length of the fractures.