In this paper we use gravity data to study fine crustal structure and seismogenic environment beneath Beijing and its surrounding regions. Multi-scale wavelet analysis method is applied to separating gravity fields. Logarithmic power spectrum method is also used to calculate depth of gravity field source. The results show that the crustal structure is very complicated beneath Beijing and its surrounding areas. The crustal density exhibits laterally inhomogeneous. There are three large scale tectonic zones in North China, i.e., WNW-striking Zhangjiakou-Bohai tectonic zone (ZBTZ), NE-striking Taihang piedmont tectonic zone (TPTZ) and Cangxian tectonic zone (CTZ). ZBTZ and TPTZ intersect with each other beneath Beijing area and both of them cut through the lithosphere. The upper and middle crusts consist of many small-scale faults, uplifts and depressions. In the lower crust, these small-scale tectonic units disappear gradually, and they are replaced by large-scale tectonic units. In surrounding regions of Beijing, ZBTZ intersects with several other NE-striking tectonic units, such as Cangxian uplift, Jizhong depression and Shanxi Graben System (SGS). In west of Taihangshan uplift, gravity anomalies in upper and middle crusts are correlated with geological and topographic features on the surface. Compared with the crust, the structure is comparatively simple in uppermost mantle. Earthquakes mainly occurred in upper and middle crusts, especially in transitional regions between high gravity anomaly and low gravity anomaly. Occurrence of large earthquakes may be related to the upwelling of upper mantle and asthenosphere heat flow materials, such as Sanhe earthquake (MS8.0) and Tangshan earthquake (MS7.8).