We present a comparison of the Solar Dynamics Observatory (SDO) analysis of NOAA Active Region (AR) 11158 and numerical simulations of flux-tube emergence, aiming to investigate the formation process of the flareproductive AR. First, we use SDO/Helioseismic and Magnetic Imager (HMI) magnetograms to investigate the photospheric evolution and Atmospheric Imaging Assembly (AIA) data to analyze the relevant coronal structures. Key features of this quadrupolar region are a long sheared polarity inversion line (PIL) in the central δ-sunspots and a coronal arcade above the PIL. We find that these features are responsible for the production of intense flares including an X2.2-class event. Based on the observations, we then propose two possible models for the creation of AR 11158 and conduct flux emergence simulations of the two cases to reproduce this AR. Case 1 is the emergence of a single flux tube, which is split into two in the convection zone and emerges at two locations, while Case 2 is the emergence of two isolated but neighboring tubes. We find that, in Case 1, a sheared PIL and a coronal arcade are created in the middle of the region, which agrees with the AR 11158 observation. However, Case 2 never build a clear PIL, which deviates from the observation. Therefore, we conclude that the flare-productive AR 11158 is, between the two cases, more likely to be created from a single split emerging flux than two independent flux bundles.