Sharp price fluctuations and increasing environmental and distributional concerns, among other issues, have led to a renewed academic interest in energy demand. In this paper we estimate, for the first time in Spain, an energy demand system with household microdata. In doing so, we tackle several econometric and data problems that are generally recognized to bias parameter estimates. This is obviously relevant, as obtaining correct price and income responses is essential if they may be used for assessing the economic consequences of hypothetical or real changes. With this objective, we combine data sources for a long time period and choose a demand system with flexible income and price responses. We also estimate the model in different sub-samples to capture varying responses to energy price changes by households living in rural, intermediate and urban areas. This constitutes a first attempt in the literature and it proved to be a very successful choice.