Traditional seismic statistics are usually limited to regular statistical scope, or require a lot of manual intervention somehow. To solve this problem, in this paper we introduce the spatial statistics method based on geographic information system, and apply it to seismic spatial statistics, which can efficiently and rapidly excavate the interrelation among various elements containing geographic information. Take a case-study to demonstrate the basic and detailed technological process of spatial statistics, we analyzed the seismic data(M_L ≥ 1.0) of the Sichuan and Yunnan region. The results show that 0.5°×0.5° is the most suitable scale size for this case-study, and the region with most earthquake events in 2019 is the Weiyuan-Changning region on the southwestern margin of the Sichuan Basin. In addition, the Longmenshan region, junction area of the Longmenshan fault zone and Xianshuihe-Xiaojiang fault zones, Lijiang-Xiaojinhe region, and the southern margin of the Yunnan Province have relative high rate of earthquake occurrence in 2019. Besides, we applied the method to statistically study the earthquake energy released of every province during 2015-2019. We found that the Sichuan Province has released the largest amount of seismic energy, Xinjiang and Xizang Autonomous Region has gone through the second largest amount energy. Besides the Sichuan Province, the land area of Taiwan, Sichuan Province and Xizang Autonomous Region have experienced very high seismic energy per 10⁴ square kilometers by area or per 10⁴ people by population. Through above two case-studies we show the preliminary application of GIS in the earthquake statistical analysis.