Synthesising Residential Electricity Load Profiles at the City Level Using a Weighted Proportion (Wepro) Model

: It is important to understand residential energy use as it is a large energy consumption sector and the potential for change is of great importance for global energy sustainability. A large energy-saving potential and emission reduction potential can be achieved, among others, by understanding energy consumption patterns in more detail. However, existing studies show that it requires many input parameters or disaggregated individual end-uses input data to generate the load proﬁles. Therefore, we have developed a simpliﬁed approach, called weighted proportion (Wepro) model, to synthesise the residential electricity load proﬁle by proportionally matching the city’s main characteristics: Age group, labour force and gender structure with the representative households proﬁles provided in the load proﬁle generator. The ﬁndings indicate that the synthetic load proﬁles can represent the local electricity consumption characteristics in the case city of amsterdam based on time variation analyses. The approach is in particular advantageous to tackle the drawbacks of the existing studies and the standard load model used by the utilities. Furthermore, the model is found to be more e ﬃ cient in the computational process of the residential sector’s load proﬁles, given the number of households in the city that is represented in the local proﬁle.


Introduction
The residential energy sector plays a crucial role in achieving greater energy efficiency and emissions reduction goals. Studies have suggested that residential energy use is of great importance in ensuring global energy sustainability, given its energy-saving potential [1,2]. The International Energy Agency (IEA) has calculated that the residential sector contributes about 25% of energy consumption and 17% of carbon dioxide (CO 2 ) emissions globally. It is therefore, essential to understand the residential energy consumption patterns locally to allow for an assessment of the energy-saving potential in the sector. However, lack of accessibility to measured high-resolution electricity consumption data at the city level such as smart-meter data and time use survey (TUS) data makes it difficult to understand the characteristics of electricity consumption locally. Research into this aspect will improve our understanding of residential electricity load profiles, which can be used to achieve improvements in energy efficiency as the residential sector has a major potential for energy savings [3]; to reduce CO 2 emissions as extensive studies have identified that household behaviour has a significant impact on consumption and emission [3,4]; and to optimise energy management [5] as these types of studies have supported transmission grid planning for better energy management [5,6]. This suggests that energy  [12,34,35,37] Furthermore, some case studies have employed data-mining techniques to identify residential electricity load profiles [4,8,34,35]. Recent studies have proposed data-mining-based methods such as K-means [4,29,34], hierarchical [29,35] and fuzzy algorithms for purposes of electricity load profile modelling [29]. A clustering-based framework to analyse household electricity consumption patterns using a k-means algorithm has been proposed for a study conducted in China. The clustering method was selected since the electricity consumption patterns in the data were relatively smooth. A k-means algorithm was applied because it works considerably faster than other cluster algorithms, and it was easier to interpret the clustering results. The analysis was conducted in three consecutive stages: Bottom-up [3,[14][15][16][17][18][19][20][21][22]28] Energies 2020, 13, x FOR PEER REVIEW 4 of 28 Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Others [12,34,35,37] Furthermore, some case studies have employed data-mining techniques to identify residential electricity load profiles [4,8,34,35]. Recent studies have proposed data-mining-based methods such as K-means [4,29,34], hierarchical [29,35] and fuzzy algorithms for purposes of electricity load profile modelling [29]. A clustering-based framework to analyse household electricity consumption patterns using a k-means algorithm has been proposed for a study conducted in China. The clustering method was selected since the electricity consumption patterns in the data were relatively smooth. A k-means algorithm was applied because it works considerably faster than other cluster algorithms, and it was easier to interpret the clustering results. The analysis was conducted in three consecutive stages: Top-down [12,13] Methods Energies 2020, 13, x FOR PEER REVIEW 4 of 28 Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works.  [12,34,35,37] Furthermore, some case studies have employed data-mining techniques to identify residential electricity load profiles [4,8,34,35]. Recent studies have proposed data-mining-based methods such as K-means [4,29,34], hierarchical [29,35] and fuzzy algorithms for purposes of electricity load profile modelling [29]. A clustering-based framework to analyse household electricity consumption patterns using a k-means algorithm has been proposed for a study conducted in China. The clustering method was selected since the electricity consumption patterns in the data were relatively smooth. A k-means algorithm was applied because it works considerably faster than other cluster algorithms, and it was easier to interpret the clustering results. The analysis was conducted in three consecutive stages: Stochastics model [18][19][20] Energies 2020, 13, x FOR PEER REVIEW 4 of 28 Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works.  [12,34,35,37] Furthermore, some case studies have employed data-mining techniques to identify residential electricity load profiles [4,8,34,35]. Recent studies have proposed data-mining-based methods such as K-means [4,29,34], hierarchical [29,35] and fuzzy algorithms for purposes of electricity load profile modelling [29]. A clustering-based framework to analyse household electricity consumption patterns using a k-means algorithm has been proposed for a study conducted in China. The clustering method was selected since the electricity consumption patterns in the data were relatively smooth. A k-means algorithm was applied because it works considerably faster than other cluster algorithms, and it was Machine learning techniques [4,8,34,35] Energies 2020, 13, x FOR PEER REVIEW 4 of 28 Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works.  [12,34,35,37] Furthermore, some case studies have employed data-mining techniques to identify residential electricity load profiles [4,8,34,35]. Recent studies have proposed data-mining-based methods such as K-means [4,29,34], hierarchical [29,35] and fuzzy algorithms for purposes of electricity load profile modelling [29]. A clustering-based framework to analyse household electricity consumption patterns using a k-means algorithm has been proposed for a study conducted in China. The clustering method was selected since the electricity consumption patterns in the data were relatively smooth. A k-means Others Load profile's aspect Energies 2020, 13, x FOR PEER REVIEW 4 of 28 Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works.  [12,34,35,37] Furthermore, some case studies have employed data-mining techniques to identify residential electricity load profiles [4,8,34,35]. Recent studies have proposed data-mining-based methods such as K-means [4,29,34], hierarchical [29,35] and fuzzy algorithms for purposes of electricity load profile modelling [29]. A clustering-based framework to analyse household electricity consumption patterns Behavioural [3,[14][15][16][17][19][20][21] Energies 2020, 13, x FOR PEER REVIEW 4 of 28 Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Furthermore, some case studies have employed data-mining techniques to identify residential electricity load profiles [4,8,34,35]. Recent studies have proposed data-mining-based methods such as K-means [4,29,34], hierarchical [29,35] and fuzzy algorithms for purposes of electricity load profile Behavioural-psychology [14,15] Energies 2020, 13, x FOR PEER REVIEW 4 of 28 Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Temporal profile and spatial profile [12,21] Measured data: Smart-meter data, utilities data [14][15][16][17][18] Energies 2020, 13, x FOR PEER REVIEW 4 of 28 Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Transmission System Operators (ENTSO-E) represents most of the electricity Transmission System Operators (TSOs) across Europe. In fact, the data consists of all sectors: residential, industrial and others which is also required to be synthesised if we want to take the residential part of this national load profile. In fact, Amsterdam might have a different residential profile load profile than the national's residential profile. The third dataset is the residential electricity load profiles dataset provided by NEDU [33], which will be presented in the Section 4. Therefore, a future study would be followed to improve current work when there is more data available. In addition, Table 1 provides an initial overview of the important categories in the load profiles studies based on the discussion in the related works. Others [12,34,35,37] Furthermore, some case studies have employed data-mining techniques to identify residential electricity load profiles [4,8,34,35]. Recent studies have proposed data-mining-based methods such as K-means [4,29,34], hierarchical [29,35] and fuzzy algorithms for purposes of electricity load profile modelling [29]. A clustering-based framework to analyse household electricity consumption patterns using a k-means algorithm has been proposed for a study conducted in China. The clustering method was selected since the electricity consumption patterns in the data were relatively smooth. A k-means Energies 2020, 13, 3543 5 of 29 algorithm was applied because it works considerably faster than other cluster algorithms, and it was easier to interpret the clustering results. The analysis was conducted in three consecutive stages: holidays, seasonal and shifting phenomena [34]. Similarly, our study also clusters the load-profile analysis into three stages: seasonal, monthly and typical seasonal days. Another case study in China employed hierarchical clustering, which is widely recognised in the context of pattern recognition, because it is easy to operate, efficient and practical [35]. A quantitative analysis approach based on association rule mining (ARM) was proposed in [4] in order to identify the impacts of household characteristics (HCs) on residential electricity consumption patterns. In any case it is assumed that the load profile data on weekdays are somehow more typical and significant than those on weekend days, while our work has covered both the weekday load profiles and weekend day load profiles through selected typical days [4].
A range of statistical analysis methods have also been applied in order to model residential electricity load profiles [6,28,[37][38][39], including determination of the key drivers of residential peak electricity demand. Some studies provided panel datasets including data from smart-meters [24,26,40]. A model was developed using Australian data for the greater Sydney region to analyse and model residential peak demand by providing both daily and seasonal patterns [37]. The analysis was in line with the results of multiple studies showing that peak residential electricity consumption was significantly influenced by the climate and the demand for cooling. In another study, hourly residential electricity consumption was used to estimate the Monte Carlo stochastic building-stock energy model of the dwellings in the sample and the climate data sources [28]. An error analysis was performed using normalised root mean square error (NRMSE), normalised mean absolute error (NMAE), maximum absolute difference (MAXAD) and maximum relative difference (MAXRD). The results from the modelling were validated using the hourly energy equations and electricity consumption data and the uncertainty of the Monte Carlo model was calculated using multiple runs as a sample. When combined with knowledge of user behaviour, this bottom-up building-stock approach, which uses energy performance certificate (EPC) databases, can be used to estimate aggregate mean hourly electricity consumption. In this case, calibration was required to develop urban energy models. This also indicated that the outdoor air temperature had a significant influence on the model [28].

Electricity Consumption Studies in The Netherlands
As an overview, some studies in relation with the electricity consumption in the case study's country are provided. The household electricity consumption constitutes approximately twenty percent of the total energy consumption in The Netherlands [41]. Behavioural profiles of electricity consumption can be determined according to Dutch household and dwelling characteristics [16,17,42]. A study based on collected questionnaires relating to the dwellings above in winter 2008 showed that household size, dwelling type, use of dryers, washing cycles and number of showers influence electricity consumption significantly [43]. Furthermore, a model-based analysis [41] has been performed to explore the effects of smart-meter adoption, occupant behaviour and appliance efficiency on reducing electricity consumption in relation to CO 2 emissions in The Netherlands. The paper looked at electricity consumption by end-users, projecting the best-and worst-case scenarios for carbon intensity annually. All cases assumed that carbon intensity would not increase in the future under current Dutch and European policies [41]. A real-life assessment of the effect of smart electrical appliances was conducted among Dutch households with a dynamic electricity tariff, an energy management system and a smart washing machine [29]. The results showed changes in laundry behaviour and thus electricity usage. The households regularly used the automation that came with smart washing machines [44]. The results of the study are interesting and could be a focus in our future work.
In relation to the residential Dutch load profiles, a recent study includes the local impact of an increasing penetration of photovoltaic (PV) panels and heat pumps (HPs) using the load measurements from three Dutch areas. It shows that the average daily load profile, without photovoltaics (PVs) and heatpumps (HPs) in all areas resembles the standard residential load profile. However, because of Energies 2020, 13, 3543 6 of 29 a shift from gas to electric stoves the time of peak load occurs earlier in the day [11]. We have also mentioned another profile generator ALPG [16,17] in our review, which is applied in our model. It is an open source load profile generator developed based on Dutch dwelling setting. The generated load profile is compared with measured data in Lochem (The Netherlands) over one year. It indicates a similar statistical trend, although some minor differences were identified, for instance the static standby power usage from the ALPG is too flat [16].
In brief, our contributions in this paper are the following: (1) We have developed a simplified method for modelling residential electricity load profiles in cities using Weighted proportion (Wepro) model that reflects local characteristics. (2) We introduce a practical and efficient approach to synthesize electricity load profiles, which does not require many input parameters or disaggregated individual end-uses input data to generate the load profiles. (3) We assess residential electricity load profiles based on time-division concepts: seasonal variation, monthly variation, typical seasonal days and hourly variation. The approach adopted here is illustrated the application in the case study with simple examples of the proportion adjustments of the city's profiles and household's profiles.
The rest of the paper is structured as follows: Section 2 describes the research design; Section 3 presents the results, which is the application of the method for the amsterdam case; Section 4 evaluates and discusses the results; and Section 5 concludes the paper and present the research implications for future work.

Materials and Methods
The proposed method consists of four phases: data collection, data pre-processing, data-modelling and load-profile analyses. Data collection can be challenging, frustrating and time-consuming, especially when we want to acquire high-resolution time series data. In order to generate the hourly profile of residential electricity consumption in cities, it is required to provide city's main input data on population information such as on gender, age groups and labour force. Furthermore, it is essential to identify the required dataset or information such as national holidays per year, solar irradiation dataset and outdoor temperature dataset. All these data should cover the same periods of time. In this work, the proposed model is validated by the case-study city of the H2020 ClairCity project presented here, namely amsterdam (The Netherlands). ClairCity is a research project modelling air pollution and carbon emissions. The project identifies current air emissions or pollutant concentrations by technology and citizens' activities, behaviour and practices in six pilot cities or regions: amsterdam, Bristol, Aveiro, Liguaria, Ljubljana and Sosnowiec. The aim is to develop locally specific policy packages in which clean-air, low-carbon, healthy futures are quantified, modelled and analysed [45][46][47][48][49][50][51].
In data collection and pre-processing phases, it is important to study the latter comprehensively, as it can improve data quality and the accuracy of the result [52]. Data corruption, missing values and outliers are the commonest problems in data-processing [52,53]. In general, there are four tasks in data pre-processing: cleaning, transformation, integration and reduction [52,54,55]. Table 2 summarises the common problems of data pre-processing tasks and their solutions: In this work, the data collected from amsterdam (The Netherlands), are in the form of a panel dataset, which is a cross-sectional data sample at specific point in time [52]. The panel dataset consists of information on age groups, the gender structure of each age group, the labour force, national holidays, solar irradiation and temperature datasets. The information on age groups, gender and the labour force are obtained indirectly [56][57][58] from Central Bureau Statistics (CBS, The Netherlands). In this case, we have elected to model the load profile for 2015. The population age is grouped into three groups: 0-17 years old, 18-64 years old and above 64 years old. The unemployment rate is recorded as 6.7% [56]. The labour force and age groups data are not in the form of datasets. Both of them provide information on the share of employment and unemployment, and the share of population's age groups and gender structure in the city, during the selected period. Therefore, there is no pre-processing technique is required in this case as well as for the solar irradiation dataset provided in ALPG. Data on public holidays are integrated into LPG's model as one of the independent inputs, Energies 2020, 13, 3543 7 of 29 like the temperature dataset. The temperature dataset and solar irradiation dataset are retrieved from the Royal Netherlands Meteorological Institute (KNMI), the official Dutch national weather service. More specifically for temperature, we selected the data from the 240 Schiphol weather station, which is the nearest station to amsterdam and is in the same region of Noord-Holland. In this dataset, there is no missing values, noisy or inconsistent data. A reduction technique is applied, since the station code variable is not required in the modelling tool. Furthermore, due to the different standards between the data source and LPG's format. We transformed the dataset from .txt to .csv by reducing the first variable, station code, and normalising the temperature value. As mentioned, we have done data pre-processing tasks and documenting our specific work in relation to the data pre-processing steps in more details is in preparation. Table 2. Data pre-processing: The tasks perform in data pre-processing include their common problems and solutions of these problems [52,55].

Data Modelling
In the data modelling we will apply the Wepro model to synthesise the residential electricity load profile at the city level through the household profile generators namely LPG and ALPG.

Weighted Proportion (Wepro) Model
The Wepro model is a simplified approach to model residential electricity load profiles in cities by adjusting and matching the proportion of city's weighted profiles with the households' profiles through the existing household profile generators. First, it is necessary to collect information on the citizens' age groups (AG), gender (GD) and labour force (LF). In this case, a figure for annual electricity consumption is not required, since we only focus on providing the share of hourly electricity load profiles. Second, we coupled the share of age groups and labour force and applied this share to proportionally fit the total population. The population is categorised into three groups by age: 0-17 years old, 18-64 years old and over 64 years old. Thus, the sum of the composition of these age groups represents the city's population by age group is expressed in Equation (1): where Tag is the total share of the age groups' share in the city. AG 1 is the age group for people aged 0 to 17, and AG 2 for people aged 18 to 64 and AG 3 for people over the age of 64. In more detail, each age group has gender information, although we can also identify gender information at the higher level of the age groups, giving totals for each gender in the city. In this model, more details on the gender composition of each age group is required as expressed in Equation (2): where Tmf is the total share of male's share and female's share in the city. Ml is Male and Fm is Female. We also need to identify the city's labour force composition. The shares of employment and unemployment represent the city's labour force is formed in Equation (3): where Tlf is the total share of employment's share and unemployment's share in the city. Em is Employment and Un = Unemployment. The labour force data are measured on the basis of the labour force population, which is only derived from one of the age groups. In this case, the labour force is included in AG 2 = 18-64 years old. Here the labour force is the proper set of age groups, labour force being an aspect of the age groups but not equal to age groups as shown in Equation (4): As mentioned, we employ the household profile generators in this case LPG and ALPG to generate the household load profiles. The first step is to select the household profiles to be modelled by the profile generators. The fundamental consideration is that the selected household profiles in the profile generators should represent the city's characteristics in term of age groups, gender structure and labour force, this being the focus of our study. This means that the selected household profiles should represent the city's profiles proportionally as depicted in Figure 1.
where Tlf is the total share of employment's share and unemployment's share in the city. Em is Employment and Un = Unemployment. The labour force data are measured on the basis of the labour force population, which is only derived from one of the age groups. In this case, the labour force is included in AG2 = 18-64 years old. Here the labour force is the proper set of age groups, labour force being an aspect of the age groups but not equal to age groups as shown in Equation (4): As mentioned, we employ the household profile generators in this case LPG and ALPG to generate the household load profiles. The first step is to select the household profiles to be modelled by the profile generators. The fundamental consideration is that the selected household profiles in the profile generators should represent the city's characteristics in term of age groups, gender structure and labour force, this being the focus of our study. This means that the selected household profiles should represent the city's profiles proportionally as depicted in Figure 1. • Capacity, fairness of allocation and rounding number We apply the capacity model based on the Amsterdam's age groups share in Figure 1 for selecting which household profiles to be modelled. The main goal is to determine the number of the occupants's profiles to be modelled as shown in the following expression of Equation (5): where Tamt is the total number of the occupants' profiles. AG1wt is the number in age group 1 based on it's weight. AG2wt is the number in age group 2 based on it's weight and AG3wt is the number in age group 3 based on it's weight. The share of the occupants for each group are converted to decimal form to provide the results of the total number of occupant-profiles from each age group. Furthermore, the capacity model can also be extended to determine the gender of the selected profiles as expressed in Equation (6) if it is supported by the profile generators. In this case, it is applicable to LPG, since LPG provides detail characters of the occupants' gender information: Here Tg is the total number of combinations of the occupants' gender. AG1m is the share of males in age group 1. AG1f is the share of females in age group 1. AG2m is the share of males in age • Capacity, fairness of allocation and rounding number We apply the capacity model based on the amsterdam's age groups share in Figure 1 for selecting which household profiles to be modelled. The main goal is to determine the number of the occupants's profiles to be modelled as shown in the following expression of Equation (5): where Tamt is the total number of the occupants' profiles. AG1wt is the number in age group 1 based on it's weight. AG2wt is the number in age group 2 based on it's weight and AG3wt is the number in age group 3 based on it's weight. The share of the occupants for each group are converted to decimal form to provide the results of the total number of occupant-profiles from each age group. Furthermore, the capacity model can also be extended to determine the gender of the selected profiles as expressed in Equation (6) if it is supported by the profile generators. In this case, it is applicable to LPG, since LPG provides detail characters of the occupants' gender information: Here Tg is the total number of combinations of the occupants' gender. AG1m is the share of males in age group 1. AG1f is the share of females in age group 1. AG2m is the share of males in age group 2. AG2f is the share of females in age group 2. AG3m is the share of males in age group 3 and AG3f is the share of females in age group 3. In this case a widely used fairness sharing technique called max-min fairness can be applied in sharing the allocations if it is required.
Therefore, the application of the Wepro model to the case-study city namely amsterdam is as follows: First, the city's population is represented by the sum of the composition of age groups in amsterdam. We grouped the city's age groups into three categories: 0-17 years old = 17.5%; 18-64 years old = 70.3%; and above 64 years old = 12.2% [57,58] using the formula in Equation (1): In more detail, the gender structure is classified into three age groups. For the age group of 0 to 17-year-olds, 51.58% are male and 48.42% female. In the age group of 18-to 64-year-olds, 50.24% are male and 49.75% female. Finally, for the age group above 65, we identified 46.24% male and 53.75% female [57,58]. Therefore, Equation (2) is presented to identify the gender at the city level: Tm f = 49.5% + 50.5% Tm f = 100% Furthermore, the labour force data are measured on the basis of the labour force population, which is only derived from age group among 18-to 64-year-olds. The unemployment rate is recorded as 6.7% [56]. In this case, Equation (3) is used to identify the employment and unemployment shares.

Tl f = 100%
Here, Equation (4) is used where the labour force is the proper set of age groups, labour force being an aspect of the age groups but not equal to age groups: We coupled the share of age groups and labour force and applied the Proportional matched profile to the total population as the city's main characteristics. Therefore, as displayed in Figure 2, the Amsterdam's main profile should reflect: The age groups, labour force and gender classes. This means that from the age group percentage: The aged 0-17 group is nearly 20%, aged 18-64 is 70% and the rest 10% is for aged over 64. From this 70% where the aged 18-64, there is about 93% of this age group are people with work and the rest are not working. Furthermore, each age group illustrates a slight difference in the share of gender information, except for the aged over 64, where the female populations are slightly more dominant than the male populations. We coupled the share of age groups and labour force and applied the Proportional matched profile to the total population as the city's main characteristics. Therefore, as displayed in Figure 2, the Amsterdam's main profile should reflect: The age groups, labour force and gender classes. This means that from the age group percentage: The aged 0-17 group is nearly 20%, aged 18-64 is 70% and the rest 10% is for aged over 64. From this 70% where the aged 18-64, there is about 93% of this age group are people with work and the rest are not working. Furthermore, each age group illustrates a slight difference in the share of gender information, except for the aged over 64, where the female populations are slightly more dominant than the male populations.

•
Capacity, fairness of allocation and rounding number Furthermore, Equation (5) is presented, where the weighted city's age group values are applied into a simple capacity model, in order to determine the capacity of the allocation. Therefore, based on the weighted values, we have ten capacity of the households profiles. It means, we can only select maximum ten occupants from the household profiles generators: Tamt = 17.5% + 70.3% + 12.2% Tamt = 1.75 + 7.03 + 1.2 Furthermore, if it is supported by the profile generators, the capacity model can also be extended to determine the gender of the selected profiles as expressed in Equation (6). In this case, it is applicable to LPG, since LPG provides detail characters of the occupants gender information: Tg = (1.03 + 0.96) + (3.51 + 3.48) + (0.46 + 0.54) As shown in the Equation (5), age group 1 has two allocations, age group 2 has seven allocations and age group 3 has one allocation. Thus, there are currently two resources for two allocations, which after the division between them, resulting in 1. Furthermore, AG1m has an excess of 0.03, where the excess can be taken and divided among the remaining demands, which is only AG1f. Therefore, AG1f = 1. As a result of the capacity and fairness of allocation model depicted in Figure 3, there will be two occupants: one male and one female in age group 0-17. Furthermore for the case Age group 2 and Age group 3, we cannot fully apply the max-min fairness. We simply apply rounding number because we have only two resources per age group. For instance, for age group 3, there are two resources for only one allocation. Therefore, rounding number is applied to the highest weight between the resources. working people and one person not working. Considering the gender share is quite balance in this age group, then it is either four females and three males, or four males and three females in the occupants' list. Lastly, for the aged over 64, which has only one allocation, based on the results of the capacity and fairness of allocation model, we apply rounding value to the one which has the highest share to represent the senior age group. Therefore, we selected a female senior to represent this age group.  Figure 3. The application of the Wepro model's structure to the case-study city, namely Amsterdam, The Netherlands. It consists of the Amsterdam age group share, labour force composition share and gender share of each age group and their capacity of the occupants to be modelled. As a result, the age group 18-64 should consist of seven adult occupants with six of them working people and one person not working. Considering the gender share is quite balance in this age group, then it is either four females and three males, or four males and three females in the occupants' list. Lastly, for the aged over 64, which has only one allocation, based on the results of the capacity and fairness of allocation model, we apply rounding value to the one which has the highest share to represent the senior age group. Therefore, we selected a female senior to represent this age group.

Profile Generators: LPG and ALPG
To produce the load profiles of the selected households profiles as the result of the Wepro model between the city's main characteristics and the households occupants, we use LPG and ALPG as the load profiles generators. Thus, in this case we optimise a bottom-up approach provided by the generators, scale-up from the household level to the city level based on the down-scale task perform previously in the weighting model, and employ the profile generator's model at the former level.
The main reason of choosing LPG and ALPG because both of them are developed based on behavioral model, which is in line with ClairCity project's goal to model the citizen's behaviour. LPG's model has been selected for use in our model, as it offers a mature model with which to synthesise household energy load profiles based on various occupants' profiles. Pflugrandt has developed the model with a strong focus on modelling the behavioural aspect. The basic elements for modelling a single household in Figure 4 are the desire to do so and expressions of the need to do something. The model specifies weight, threshold and decay time as desired properties [15].
Weight is the relative weight of a need compared to all a person's other needs. In selecting for the next action, the minimisation of the deviation requirement is used as a criterion, the weighting acting as a multiplier in this calculation. Threshold determines when the person really feels a need, that is, when it is included in the next action selection of the calculation. For example, in reality there is usually no eating after lunch because only 10% of the hunger sensation is evident. Instead, one generally waits until a noticeable feeling of hunger has built up before having dinner. Finally, Decay Time describes the half-life, until 50% of the requirement is reached. It has been found that activities at 50% threshold mostly after the two to three times the decay time, depending on the weighting and the other available activities. The decay constant is calculated from the decay time by which the current value of the need is multiplied in each time step [15]. When creating households, it has been found that activities at the 50% threshold are usually executed after two to three times the decay time, depending on the weighting and the other available activities. Furthermore, besides desire, it is also essential to identify the individual's properties (age, gender, sick leave in the year, average duration of illness, needs when healthy, needs when ill) and load type, which in this case is electricity [14,15]. LPG provides various pre-defined German household profiles. To produce the load profiles of the selected households profiles as the result of the Wepro model between the city's main characteristics and the households occupants, we use LPG and ALPG as the load profiles generators. Thus, in this case we optimise a bottom-up approach provided by the generators, scale-up from the household level to the city level based on the down-scale task perform previously in the weighting model, and employ the profile generator's model at the former level.  The main reason of choosing LPG and ALPG because both of them are developed based on behavioral model, which is in line with ClairCity project's goal to model the citizen's behaviour. LPG's model has been selected for use in our model, as it offers a mature model with which to synthesise household energy load profiles based on various occupants' profiles. Pflugrandt has developed the model with a strong focus on modelling the behavioural aspect. The basic elements for modelling a single household in Figure 4 are the desire to do so and expressions of the need to do something. The model specifies weight, threshold and decay time as desired properties [15].
Weight is the relative weight of a need compared to all a person's other needs. In selecting for the next action, the minimisation of the deviation requirement is used as a criterion, the weighting acting as a multiplier in this calculation. Threshold determines when the person really feels a need, that is, when it is included in the next action selection of the calculation. For example, in reality there is usually no eating after lunch because only 10% of the hunger sensation is evident. Instead, one generally waits until a noticeable feeling of hunger has built up before having dinner. Finally, Decay Time describes the half-life, until 50% of the requirement is reached. It has been found that activities at 50% threshold mostly after the two to three times the decay time, depending on the weighting and the other available activities. The decay constant is calculated from the decay time by which the current value of the need is multiplied in each time step [15]. When creating households, it has been The second profile generator used in our work is ALPG. ALPG employs household occupancy profiles generated by a simple behavioural model, which creates consistent profiles for the devices. The devices' flexibility is specified through four classes: timeshiftables, buffer-timeshiftables, buffers and curtailable. The inflexible electricity profiles are grouped into the following categories: stand-by load, consumer electronics, lighting, inductive devices, fridges, and other. Furthermore, to show annual electricity consumption, the individual profiles are scaled in magnitude, making it easier to alter the profile if there is a change in electricity usage by the external factors. An example of such a change could be the adoption of a new technology, for instance, light-emitting diode (LED) lights. Moreover, the following classes in Figure 5 are implemented in the simulation model: neighbourhood, household, person, device, house, writer and ALPG. Electricity usage in a typical Dutch setting is the focus of ALPG, which is also in line with our work in modelling residential electricity load profiles, with Amsterdam as the case-study city [16,17].
Furthermore, we after applying the capacity allocation into LPG and ALPG the following are closest profiles that reflect the city's proportion of the age groups, gender and labour force mentioned above.
• LPG The following are the simplified Wepro-based selected pre-defined households profiles in LPG although there could be also several other options that may fulfill the Wepro model composition: found that activities at the 50% threshold are usually executed after two to three times the decay time, depending on the weighting and the other available activities. Furthermore, besides desire, it is also essential to identify the individual's properties (age, gender, sick leave in the year, average duration of illness, needs when healthy, needs when ill) and load type, which in this case is electricity [14,15]. LPG provides various pre-defined German household profiles.
The second profile generator used in our work is ALPG. ALPG employs household occupancy profiles generated by a simple behavioural model, which creates consistent profiles for the devices. The devices' flexibility is specified through four classes: timeshiftables, buffer-timeshiftables, buffers and curtailable. The inflexible electricity profiles are grouped into the following categories: stand-by load, consumer electronics, lighting, inductive devices, fridges, and other. Furthermore, to show annual electricity consumption, the individual profiles are scaled in magnitude, making it easier to alter the profile if there is a change in electricity usage by the external factors. An example of such a change could be the adoption of a new technology, for instance, light-emitting diode (LED) lights. Moreover, the following classes in Figure 5 are implemented in the simulation model: neighbourhood, household, person, device, house, writer and ALPG. Electricity usage in a typical Dutch setting is the focus of ALPG, which is also in line with our work in modelling residential electricity load profiles, with Amsterdam as the case-study city [16,17].  [17] that shows the cardinality of a class in relation to another. The example of one-to-one (1..1) relationship is depicted between Household and House, where a household lives in a house and a house belongs to a household. The one-to-many(1..*) relationship is shown between Household and Device, where a household has one or more devices, and each device belongs to a household. Each class from these multiple classes represents a part of the model, which makes the software flexible to be extended in the future work.
Furthermore, we after applying the capacity allocation into LPG and ALPG the following are closest profiles that reflect the city's proportion of the age groups, gender and labour force mentioned above.
• LPG The following are the simplified Wepro-based selected pre-defined households profiles in LPG although there could be also several other options that may fulfill the Wepro model composition: Couple, both of whom work, with one child Couple, one at work, one at home, with one child Couple both of whom work Single with work Senior at home  [17] that shows the cardinality of a class in relation to another. The example of one-to-one (1..1) relationship is depicted between Household and House, where a household lives in a house and a house belongs to a household. The one-to-many(1..*) relationship is shown between Household and Device, where a household has one or more devices, and each device belongs to a household. Each class from these multiple classes represents a part of the model, which makes the software flexible to be extended in the future work. Furthermore, we should select some options such as which pre-defined profile to be modelled, geographic location and temperature profile based on temperature dataset that we input before, if the temperature dataset is not provided yet by LPG. That is why we need to pre-process our input data such as temperature dataset in order to be matched with LPG's format. Then we can calculate the household profile one by one which may require a computational processing time and the result is generated in comma-separated values (.CSV) file. Furthermore, we can insert these occupant's list to the Wepro composition in order to validate the model. As illustrated in Figure 6, the selected household profiles can fulfill the Wepro's model composition. Then, we generate these LPG's pre-defined households' load profiles one by one. The LPG can be downloaded free from https://www.loadprofilegenerator.de/. In generating one predefined household's load profile, after we download and open the windows program, we can go to "calculation" menu.  Figure 6. The application of the Wepro model's structure for Amsterdam's household occupancy profiles in LPG. It consists of the Amsterdam's age group share, labour force composition share and gender share of each age group, their capacity of the occupants to be modelled and the selected gender character provided in LPG.
Furthermore, we should select some options such as which pre-defined profile to be modelled, geographic location and temperature profile based on temperature dataset that we input before, if the temperature dataset is not provided yet by LPG. That is why we need to pre-process our input data such as temperature dataset in order to be matched with LPG's format. Then we can calculate the household profile one by one which may require a computational processing time and the result is generated in comma-separated values (.CSV) file.
• ALPG As shown in Table 4, the pre-defined households profiles in ALPG are not as detailed as in LPG, but they simply can fulfill the Wepro model. The pre-defined households class contains seven types of households: Single worker, dual worker, family dual worker, family single parent, dual retired and single retired. Dual profile means a couple. In this case, each type of household corresponds to a category of electricity annual consumption in Kilowatt hour and amount of occupants or persons. Figure 6. The application of the Wepro model's structure for amsterdam's household occupancy profiles in LPG. It consists of the amsterdam's age group share, labour force composition share and gender share of each age group, their capacity of the occupants to be modelled and the selected gender character provided in LPG.
• ALPG As shown in Table 4, the pre-defined households profiles in ALPG are not as detailed as in LPG, but they simply can fulfill the Wepro model. The pre-defined households class contains seven types of households: Single worker, dual worker, family dual worker, family single parent, dual retired and single retired. Dual profile means a couple. In this case, each type of household corresponds to a category of electricity annual consumption in Kilowatt hour and amount of occupants or persons. Table 4. Pre-defined households configurations in ALPG based [17]. To fulfill the Wepro model and simplify the process, we selected: one single worker, one single retired, two dual worker and one family dual worker. ALPG is open-source code and the code is available at it's github page. Figure 7 shows the snipped code of the households profiles selection, where the ALPG program runs by executing profilegenerator.py.: Furthermore, the same procedure with LPG, in Figure 8 we inserted these occupant's list to the Wepro composition in order to validate the model. Accordingly, it indicates a different result in comparison with LPG because in ALPG there is no need to identify the gender characteristics as it has simplified and consistent pre-defined profiles list as provided in Table 4. Moreover, these selected occupancy's list in ALPG may fulfill the Wepro model regardless the gender detail. In consequences, there are five generated households load profiles both in LPG and ALPG. We used the average load profile's value of these generated load profiles in the analysis. To fulfill the Wepro model and simplify the process, we selected: one single worker, one single retired, two dual worker and one family dual worker. ALPG is open-source code and the code is available at it's github page. Figure 7 shows the snipped code of the households profiles selection, where the ALPG program runs by executing profilegenerator.py.:    available at it's github page. Figure 7 shows the snipped code of the households profiles selection, where the ALPG program runs by executing profilegenerator.py.:   Figure 8. The application of the Wepro model's structure for Amsterdam's household occupancy profiles in ALPG. It consists of the Amsterdam's age group share, labour force composition share and gender share of each age group, their capacity of the occupants to be modelled and the selected predefined profiles provided in LPG. Figure 8. The application of the Wepro model's structure for amsterdam's household occupancy profiles in ALPG. It consists of the amsterdam's age group share, labour force composition share and gender share of each age group, their capacity of the occupants to be modelled and the selected pre-defined profiles provided in LPG.

Load Profile Analyses Based on Time-Division
As the model has produced residential electricity load profiles in high resolution, we focus on analysis of the load profile results in this section. The visualisation, charting and plotting of the hourly resolution are executed in Python. The residential electricity load profile is analysed at four levels: seasonal analysis, monthly analysis, days analysis and hourly analysis.

Seasonal Analysis
Electricity consumption patterns based on the seasons is interesting to distinguish, as the temperature influences the interval of the seasons. A time-division according to the four seasons in the year has been defined based on the meteorological concept (Table 5).

Monthly Analysis
Monthly characteristics are examined through the monthly average share to identify the monthly pattern of residential electricity consumption in the city. The analyses cover the 12 months electricity data from the load profiles results. Besides to identify the potential energy savings, the monthly analysis is beneficial to plan the generation and distribution of power utilities.  The meteorological concept (Table 5) is quite simple and is the most widely used, being broken down into four three-month periods. Winter has the three coldest months in the Northern Hemisphere, namely December, January and February. Spring runs from March to May, summer from June to August, and the other months belong to autumn [59]. Hence, a seasonal electricity load profiles model is proposed as follows in Equation (7): where Ts% is the total of seasons share, Sw% is the share of Winter season, Ssp% is the share of Spring season, Ssm% is the share of Summer season and Sa is the share of Autumn season.

Days Analysis
The days analysis is provided based on the hourly average share load of the days in each season and the typical days share in each season. The typical days are the selected days, of one weekday and one weekend day in each season as listed in Table 6. None of the selected days listed below in Table 6 is a national holiday in the selected city, meaning that the selected days represent people's normal daily activities on a weekday and at weekends.

Hourly Analysis
In this part, we will process the output from the load profile generators as one minute resolution and create a time series with an hourly resolution. The objective is to show the hourly characteristics of electricity consumption on the hourly average load profiles in the selected year and the hourly average of the seasonal load profiles. Hence, the average load profiles model is proposed as follows in Equation (8): where h n is the hourly average. n is the number of hours which is in the range of 0 to 23. d is the number of days in a year, which is in the range of 1 to 365 for the selected year 2015. x i is the data i-th.
To create the hourly average of the seasonal load profiles, we replace d with the number of the days in each season.

Results: Load Profile Analyses in amsterdam as the Case Study
The analyses of the generated load profiles by the model will be presented based on the time variation of the case-study city, namely amsterdam (The Netherlands). The hourly temporal results will be moderately validated by the standard average of Dutch household load profile.

Load-Profile Analyses Based on Time-Division
The generated load profiles produced by the model will be analysed at four levels: seasonal analysis, monthly analysis, days analysis and hourly analysis, where hourly resolution is the core output of our temporal profile.

Seasonal Analysis
The seasonal variations based on the meteorological concept (Table 5) are provided in two generated models: Wepro-LPG and Wepro-ALPG, where we used Equation (7). In this concept, Winter has the fewest days, Autumn has one day more compared to winter, while spring and summer have two days more compared to winter. Based on the generated load profiles of Wepro-LPG and Wepro-ALPG in Figure 9, it indicates that Winter is the highest consumption share, slightly followed by Autumn, and Spring. Both models show Summer as the lowest consumption share.

Monthly Analysis
Amsterdam's monthly electricity load share is illustrated in Figure 10, which shows a distinct profile in the Wepro-LPG model. This demonstrates that December has the highest consumption share compared to the other months, which may exhibit seasonal variations. Surprisingly, this load profile identifies May as having the second highest electricity share in 2015, followed by October, August and January. The lowest monthly consumption share is in July, which concurs with the seasonal analysis result that summer has the lowest consumption share in all load profiles. The second lowest consumption share is in March, followed by September.
The Wepro-ALPG model indicates December as having the highest consumption share, the same as in Wepro-LPG model. The second highest consumption share is in January, followed by October and November. The lowest consumption share is in June, followed by February and April.  Wepro-LPG Ts% = 26.02% + 25.22% + 23.46% + 25.28% while for Wepro-ALPG, the highest consumption share occurred in winter, slightly followed by autumn, then spring. Summer is recorded as the lowest consumption period. Wepro-ALPG Ts% = 25.91% + 24.67% + 23.91% + 25.49%

Monthly Analysis
Amsterdam's monthly electricity load share is illustrated in Figure 10, which shows a distinct profile in the Wepro-LPG model. This demonstrates that December has the highest consumption share compared to the other months, which may exhibit seasonal variations. Surprisingly, this load profile identifies May as having the second highest electricity share in 2015, followed by October, August and January. The lowest monthly consumption share is in July, which concurs with the seasonal analysis result that summer has the lowest consumption share in all load profiles. The second lowest consumption share is in March, followed by September.
Energies 2020, 13, x FOR PEER REVIEW 18 of 28 Figure 10. Monthly average electricity load share based on the results of the generated Wepro-LPG and Wepro-ALPG models.

Days Analysis
The days analysis is provided based on the hourly average share load of the days in each season which are depicted in Figure 11, and the daily load share of the typical selected days, which is illustrated in Figure 12. The Wepro-ALPG model indicates December as having the highest consumption share, the same as in Wepro-LPG model. The second highest consumption share is in January, followed by October and November. The lowest consumption share is in June, followed by February and April.

Days Analysis
The days analysis is provided based on the hourly average share load of the days in each season which are depicted in Figure 11, and the daily load share of the typical selected days, which is illustrated in Figure 12.  The days analysis is provided based on the hourly average share load of the days in each season which are depicted in Figure 11, and the daily load share of the typical selected days, which is illustrated in Figure 12.  For the winter period in Wepro-LPG, the top morning peak appears around 7 am for all weekdays. For some weekdays, the morning peak then continues with a light peak between 9 am and 10 am. The following daylight peaks are identified around lunchtime, from 12 am to 1 pm, while for weekend days, the curve show several daylight peaks from 7 am to late lunchtime, around 1 pm. The morning and daylight peaks appear around 7 am, 9 am, 11 am and 1 pm, among which 11 am is identified as the top morning peak, which may be associated with the brunch time. Furthermore, the evening peaks for all days started from about 5 pm to 8 pm and mostly reach the top value around 7 pm, while for the weekend days, the curves show a longer evening peak from 6 pm to around 7:30 pm. The weekend days show a quite higher daylight share than the weekdays' daylight share. Furthermore, for Wepro-ALPG, the top morning peak is shown around 7 am for all weekdays, followed by a light peak around 8 am on some weekdays. For some weekdays, the morning peak then continues with a light peak between 10 am to 12 am. The load share continues to increase until 2 pm, while for the weekend days, the morning peaks are characterized by a light peak around 7 am that continues to increase to 10 am. The top morning peak is around 12 am. The weekend days show a higher daylight share than the weekdays' daylight share. Most of the days indicate 6 pm as the top daylight share. Most of the days indicate 6 pm as the top evening peak, two weekdays show 5 pm as the top evening peak and the weekend days shows 7 pm as the top evening peak. Furthermore, the hourly average load share in Autumn period for Wepro-LPG illustrates the top morning peak on 6 am on most weekdays and one weekday has 7 am as the top morning peak. The following morning peak occurs on 10 am on most weekdays, followed by another daylight peak at 2 pm, while for the weekend days, the curve identifies several peaks at 6 am, 8 am, 9 am and 11 am. The morning peaks are seen around 6 am, 9 am, 11 am and an afternoon one at 2 pm, where 11 am is identified as the top morning peak, which may be associated with brunch. Furthermore, the evening peaks for all days start to increase from about 5 pm or 6 pm and all days reach the maximum evening peak on 7 pm. It is obvious that the weekend days show a quite significantly higher daylight share than the weekdays' daylight share. Lastly, the hourly average load share in Autumn period for Wepro-ALPG shows the top morning peak around 7 am for all weekdays, with a slight peak occurring before at 6 am. The next peak happens at 10 am, with the load share increasing gradually from 11 am to 3 pm on all weekdays, while for the weekend days, the morning peaks start with a slight peak at 7 am, and the load share keeps increasing until it reaches another peak at 10 am. It then increases further to reach the top morning peak at 12 am. The weekend days show a higher daylight share than the weekdays' daylight share. Most of the days indicate 7 pm as the top evening peak, one day identifies 5 pm as the top evening peak, another day shows 6 pm as the top evening peak and another day has 8 pm as the top evening peak.
In addition, as an overview of the daily total share load, we present the selected typical days analysis in Figure 12, where the selected days represent weekdays and weekend days of each season. For the winter period in Wepro-LPG, the top morning peak appears around 7 am for all weekdays. For some weekdays, the morning peak then continues with a light peak between 9 am and 10 am. The following daylight peaks are identified around lunchtime, from 12 am to 1 pm, while for weekend days, the curve show several daylight peaks from 7 am to late lunchtime, around 1 pm. The morning and daylight peaks appear around 7 am, 9 am, 11 am and 1 pm, among which 11 am is identified as the top morning peak, which may be associated with the brunch time. Furthermore, the evening peaks for all days started from about 5 pm to 8 pm and mostly reach the top value around 7 pm, while for the weekend days, the curves show a longer evening peak from 6 pm to around 7:30 pm. The weekend days show a quite higher daylight share than the weekdays' daylight share. Furthermore, for Wepro-ALPG, the top morning peak is shown around 7 am for all weekdays, followed by a light peak around 8 am on some weekdays. For some weekdays, the morning peak then continues with a light peak between 10 am to 12 am. The load share continues to increase until 2 pm, while for the weekend days, the morning peaks are characterized by a light peak around 7 am that continues to increase to 10 am. The top morning peak is around 12 am. The weekend days show a higher daylight share than the weekdays' daylight share. Most of the days indicate 6 pm as the top evening peak, some weekdays identify 5 pm as the top evening peak and one weekend day shows 7 pm as the top evening peak.
The hourly average load share in Spring period for Wepro-LPG illustrates the top morning peak around 6 am for all weekdays. The following morning and daylight peaks occur around lunchtime from 12 am to 2 pm, while for weekend days, the curves show several morning peaks from 6 am to late lunchtime around 2 pm. The morning and daylight peaks are shown around 6 am, 9 am, 11 am and 2 pm, where 11 am is identified as the top morning peak, which may be associated with brunch. Furthermore, the evening peaks for all days start from about 6 pm to 8 pm and mostly reach the top peak around 7 pm. It is obvious that the weekend days show a quite significant higher daylight share than the weekdays' daylight share. Furthermore, the hourly average load share in Spring period for Wepro-ALPG shows the top morning peak around 7 am for all weekdays, followed by a light peak around 8 am on some weekdays. For some weekdays, the morning and daylight peaks then continue with a light peak at 10 am, 11 am, 1 pm and 2 pm. After 7 am, the curves are continually declining until 9 am. The load share continues to increase again with a slight share from 9 am to 2 pm, while for the weekend days, the morning peaks start with a slight peak from 7 am, and gradually increase to reach the top on 10 am. It then increases slightly at 12 am, which is identified as the top morning peak in the weekend days. The weekend days show a higher daylight share than the weekdays' daylight share. Most of the days indicate 6 pm as the top evening peak, one weekday identifies 5 pm as the top evening peak and some days shows 7 pm as the top evening peak.
Furthermore, the Wepro-LPG in Summer period indicates the top morning peak on 6 am and 7 am for all weekdays. The following morning peak is around 10 am to 2 pm, while for weekend days, the curve shows several morning and daylight peaks from 6 am to late lunchtime, around 2 pm. These peaks are evident at 6 am, 7 am, 8 am, 11 am and 2 pm, where 11 am is identified as the top morning peak, which may be associated with brunch. Furthermore, the evening peaks for all days have started from about 6 pm to 8 pm, which all days reach the top peak on 7 pm. The weekend days show a quite significant higher daylight share than the weekdays' daylight share, while for Wepro-ALPG in Summer period, the top morning peak is shown around 7 am for all weekdays, followed by a light peak around 11 am for most weekdays. The load share continues to increase again with a slight share from 11 am to 2 pm, while for the weekend days, the morning peaks start with a slight peak at 7 am, then a gradual slight one increasing each next hour and reaching the top on 10 am. It then slightly increases further at 12 am. The weekend days show a higher daylight share than the weekdays' daylight share. Most of the days indicate 6 pm as the top evening peak, two weekdays show 5 pm as the top evening peak and the weekend days shows 7 pm as the top evening peak.
Furthermore, the hourly average load share in Autumn period for Wepro-LPG illustrates the top morning peak on 6 am on most weekdays and one weekday has 7 am as the top morning peak.
The following morning peak occurs on 10 am on most weekdays, followed by another daylight peak at 2 pm, while for the weekend days, the curve identifies several peaks at 6 am, 8 am, 9 am and 11 am. The morning peaks are seen around 6 am, 9 am, 11 am and an afternoon one at 2 pm, where 11 am is identified as the top morning peak, which may be associated with brunch. Furthermore, the evening peaks for all days start to increase from about 5 pm or 6 pm and all days reach the maximum evening peak on 7 pm. It is obvious that the weekend days show a quite significantly higher daylight share than the weekdays' daylight share. Lastly, the hourly average load share in Autumn period for Wepro-ALPG shows the top morning peak around 7 am for all weekdays, with a slight peak occurring before at 6 am. The next peak happens at 10 am, with the load share increasing gradually from 11 am to 3 pm on all weekdays, while for the weekend days, the morning peaks start with a slight peak at 7 am, and the load share keeps increasing until it reaches another peak at 10 am. It then increases further to reach the top morning peak at 12 am. The weekend days show a higher daylight share than the weekdays' daylight share. Most of the days indicate 7 pm as the top evening peak, one day identifies 5 pm as the top evening peak, another day shows 6 pm as the top evening peak and another day has 8 pm as the top evening peak.
In addition, as an overview of the daily total share load, we present the selected typical days analysis in Figure 12, where the selected days represent weekdays and weekend days of each season.
For the weekdays, we selected 1:

Hourly Analysis
The hourly average load profiles in a year are provided in Figure 13 based on the expression in Equation (8). The Wepro-ALPG curve illustrates an increasing load from 5 am and reaches the morning peak around 7 am. Then, the load decreases gradually until 9 am and increases again to 10 am. After that, the curve remains flat from 10 am to 1 pm although there is a subtle peak in between around 12 am or during lunchtime. It is then increases slightly from 1 pm to 3 pm, and after 3 pm the curve is increasing significantly and is reaching the evening peak around 6 pm. Furthermore, after 6 pm the curve is decreasing gradually until midnight. After that, the curve remains flat until 5 am. While for the Wepro-LPG, the morning peak starts to increase significantly from 5 am and reaches the peak at 6 am. After 6 am, the curve is gradually decreasing until 4 pm although there is a subtle peak around 11 am. It starts to increase again significantly and reaches the evening peak at 7 pm. The curve decreases significantly after 7 pm to midnight. Furthermore, it remains quite flat until 3 am. There is a slightly increase load between 3 am to 5 am before the morning peak.
Energies 2020, 13, x FOR PEER REVIEW 21 of 28 increasing significantly and is reaching the evening peak around 6 pm. Furthermore, after 6 pm the curve is decreasing gradually until midnight. After that, the curve remains flat until 5 am. While for the Wepro-LPG, the morning peak starts to increase significantly from 5 am and reaches the peak at 6 am. After 6 am, the curve is gradually decreasing until 4 pm although there is a subtle peak around 11 am. It starts to increase again significantly and reaches the evening peak at 7 pm. The curve decreases significantly after 7 pm to midnight. Furthermore, it remains quite flat until 3 am. There is a slightly increase load between 3 am to 5 am before the morning peak. Furthermore, the hourly average load per season for both models indicate the consistent curve shape within the season and model, either Wepro LPG or Wepro ALPG. The winter curve indicates the highest load profile among the hourly season curves, and it is more obvious for Wepro-ALPG whereas the winter's curve is shown as the highest load and slightly followed by autumn's curve for Wepro-LPG. The demand peaks show similarity with the hourly average in a year as the top peak is in the evening, followed by the morning peak and a subtle peak during the lunch time, with both models identifying the same peak hours as the peak hours for the hourly average in a year. In general, Wepro-LPG's load profiles for all seasons have higher load share than the Wepro-ALPG's load share for all seasons during the day. Conversely, the Wepro-ALPG's load indicate longer peak share than the Wepro-LPG's load during the night time.

Validation with Case Study's Measured Data
Validation can be done by comparing the data generated by Wepro with the city's measured data such as smart-meter data, TUS data or data from the utilities. In this case, we cannot make an in-depth validation as the measured city's data is unavailable. Therefore, a future study would follow to improve our current work when the city's measured data is available.
In practice, we can still compare our model with the standard load profiles for Dutch households published by the Energy Data Services Netherlands (EDSN) to validate whether our hourly average generated load profiles have the same trends the standard Dutch residential load profile. The average normalised standard household load in The Netherlands based on EDSN is provided in Figure 1 of [11]. It is shown that the morning peak starts to increase from 5 am, similar to both our generated models. It then reaches the peak around 10 am, while both our models identify the morning peak around 6 am to 7 am. The EDSN's load remains flat from 10 am to 13 pm, although there is a subtle peak at 12 am during lunchtime. This curve from 10 am to 13 pm is quite similar to the Wepro-ALPG Furthermore, the hourly average load per season for both models indicate the consistent curve shape within the season and model, either Wepro LPG or Wepro ALPG. The winter curve indicates the highest load profile among the hourly season curves, and it is more obvious for Wepro-ALPG whereas the winter's curve is shown as the highest load and slightly followed by autumn's curve for Wepro-LPG. The demand peaks show similarity with the hourly average in a year as the top peak is in the evening, followed by the morning peak and a subtle peak during the lunch time, with both models identifying the same peak hours as the peak hours for the hourly average in a year. In general, Wepro-LPG's load profiles for all seasons have higher load share than the Wepro-ALPG's load share for all seasons during the day. Conversely, the Wepro-ALPG's load indicate longer peak share than the Wepro-LPG's load during the night time.

Validation with Case Study's Measured Data
Validation can be done by comparing the data generated by Wepro with the city's measured data such as smart-meter data, TUS data or data from the utilities. In this case, we cannot make an in-depth validation as the measured city's data is unavailable. Therefore, a future study would follow to improve our current work when the city's measured data is available.
In practice, we can still compare our model with the standard load profiles for Dutch households published by the Energy Data Services Netherlands (EDSN) to validate whether our hourly average generated load profiles have the same trends the standard Dutch residential load profile. The average normalised standard household load in The Netherlands based on EDSN is provided in Figure 1 of [11]. It is shown that the morning peak starts to increase from 5 am, similar to both our generated models. It then reaches the peak around 10 am, while both our models identify the morning peak around 6 am to 7 am. The EDSN's load remains flat from 10 am to 13 pm, although there is a subtle peak at 12 am during lunchtime. This curve from 10 am to 13 pm is quite similar to the Wepro-ALPG model one. Furthermore, like the Wepro-LPG model, the ESDN's load is decreasing gradually to 4 pm. After that, the curve starts to increase significantly like the curves of both our generated models. The EDSN's model reaches the peak between 19:00 and 19:30 similar to the Wepro-LPG. Furthermore, the load is decreasing quite significantly until 2 am. It remains flat from 2 am to 5 am which is similar to our Energies 2020, 13, 3543 23 of 29 generated models. In general, it can be concluded that the generated hourly average share of the Wepro models have similar curve trends as the EDSN's load trend, although the morning peaks in the generated Wepro models have different time characteristics from EDSN's morning peak. Furthermore, both our models and the EDSN's model show a subtle peak during lunchtime. The evening peak occurs after dinner in the Wepro-LPG and EDSN model, while the the evening peak is occurring exactly at the dinner time in Wepro-ALPG.
As an update, the load profile data in EDSN have been moved to de Vereniging Nederlandse Energie Data Uitwisseling's (NEDU) page [33]. The data provided in NEDU's page start from year 2016, therefore data 2016 are used in this initial validation. Smart-meter data is used as a basis for the consumption/production profiles as described in 'Profielenmethodiek elektriciteit', where the documentation is available in Dutch. The raw data are provided in 15 min resolution, which show how much electricity is allocated in that 15 min. The data are obtained from 3,002,450 households type E1A in 2016. The comparison of the standard average Dutch household load profile in 2016, which has similar trends with Figure 1 of [11] for E1A residential type and the hourly average load profile in a year of Wepro-LPG and Wepro-ALPG is plotted in Figure 14.
Energies 2020, 13, x FOR PEER REVIEW 22 of 28 model one. Furthermore, like the Wepro-LPG model, the ESDN's load is decreasing gradually to 4 pm. After that, the curve starts to increase significantly like the curves of both our generated models.
The EDSN's model reaches the peak between 19:00 and 19:30 similar to the Wepro-LPG. Furthermore, the load is decreasing quite significantly until 2 am. It remains flat from 2 am to 5 am which is similar to our generated models. In general, it can be concluded that the generated hourly average share of the Wepro models have similar curve trends as the EDSN's load trend, although the morning peaks in the generated Wepro models have different time characteristics from EDSN's morning peak. Furthermore, both our models and the EDSN's model show a subtle peak during lunchtime. The evening peak occurs after dinner in the Wepro-LPG and EDSN model, while the the evening peak is occurring exactly at the dinner time in Wepro-ALPG. As an update, the load profile data in EDSN have been moved to de Vereniging Nederlandse Energie Data Uitwisseling's (NEDU) page [33]. The data provided in NEDU's page start from year 2016, therefore data 2016 are used in this initial validation. Smart-meter data is used as a basis for the consumption/production profiles as described in 'Profielenmethodiek elektriciteit', where the documentation is available in Dutch. The raw data are provided in 15 minutes resolution, which show how much electricity is allocated in that 15 min. The data are obtained from 3,002,450 households type E1A in 2016. The comparison of the standard average Dutch household load profile in 2016, which has similar trends with Figure 1 of [11] for E1A residential type and the hourly average load profile in a year of Wepro-LPG and Wepro-ALPG is plotted in Figure 14.
This simple validation is an initial check to see whether our generated load profiles resemble the standard Dutch's household load profile characteristics before going into an in-depth validation with the city's measure data. Moreover, further study in the future is required.

Discussion
Based on the results, our weighted proportion (Wepro) model can be applied to generate the residential electricity load profiles at the city level by utilising the exisiting household profile generators, either LPG or ALPG, which we have employed here, given that they both have specific behavioural profile models. The seasonal share analysis based on Wepro-LPG and Wepro-ALPG, shows each season's consumption share is in the range of 23% to 26%. The 1% share consists of approximately 80 hours of load or about 3 days of load when calculated on the basis of the hourly This simple validation is an initial check to see whether our generated load profiles resemble the standard Dutch's household load profile characteristics before going into an in-depth validation with the city's measure data. Moreover, further study in the future is required.

Discussion
Based on the results, our weighted proportion (Wepro) model can be applied to generate the residential electricity load profiles at the city level by utilising the exisiting household profile generators, either LPG or ALPG, which we have employed here, given that they both have specific behavioural profile models. The seasonal share analysis based on Wepro-LPG and Wepro-ALPG, shows each season's consumption share is in the range of 23% to 26%. The 1% share consists of approximately 80 h of load or about 3 days of load when calculated on the basis of the hourly dataset. For instance, if we compare the winter and summer seasons to the whole year in Wepro-LPG as shown in Figure 9, where winter is 26.02% and summer is 23.46%, it indicates that the electricity load in winter is almost 3% higher than in summer to the whole year, which is equal to approximately 240 h or about 9 days. In addition, both seasonal profiles indicate Winter as having the highest consumption share, which concurs with the known seasonal pattern in energy demands studies [60][61][62]. In addition, the seasonal analyses based on meteorological is important to be mentioned as some studies did not mention which time-division concept they used for analysing the seasonal electricity profile. Furthermore, the monthly analysis results illustrate that December is having the highest consumption share, which accords with the result of some monthly electricity studies [61,63].
The hourly average share based on the days in each season show that the weekend days indicate a higher daylight share than the weekdays' daylight share in both models. The result of the daily share of the selected typical days for all models indicates that most weekend days have a higher consumption share than weekdays in the same season. It concurs with an analysis of weekday and weekend variation, where weekend days show slightly more electricity use than weekdays [64]. Exception found in the Wepro-ALPG model's selected days in winter, where weekday consumption is higher than at weekend.
The hourly average load profiles identify the morning and evening peaks in Wepro-ALPG and Wepro-LPG, where the Wepro-LPG model has a higher load than the Wepro-ALPG model for both peaks. It is also identified that the evening peak has a significant higher load value than the morning peak load value in both models. All the hourly average loads in a year and per season demonstrate a consistent curve shape within season and model, either Wepro LPG or Wepro ALPG. The consistency is also shown within the curve shape of the hourly seasonal average load share with the hourly seasonal load share based on the days within the model.
As a consequence, the application of our model requires a profile generator as an external tool to match the weighted city's profile with the representative occupants' profiles at the household level, since we are not building our own profile generator. It also influences the results of the generated load profiles where they will be based on the characteristics of the developed model in profile generator, include relying on the few selected input profiles as a result of the approach taken in this study. The issue of relying on the few selected input profiles may result in the less fluctuations load profiles as shown in the Wepro-ALPG load profiles for the morning curves. The main difference of the hourly average in a year between the models is shown in the morning curve, where for Wepro-LPG after reaches the peak on 6am, the load share is declined gradually until 4 pm, with some light peaks in between, while for Wepro-ALPG, the curve declines slightly until 9 am after reaching a peak at 7 am. It increases again at 10 am and remains stable until 1 pm. This issue is also has been initially identified in [16] where the generated profiles show less fluctuations on the single household level, while the fluctuations at the neighbourhood level matched with the measured values. We assume that the less fluctuations during the morning period generated in Wepro-ALPG might be caused by the consistent pre-defined profiles in ALPG, where they are developed based on the simple behavioural model of an occupancy profile. The occupancy model for general events in ALPG is configured using mean times to change the state of a person. In this case, it is limited to the three person's states: active (being home), inactive (e.g., sleeping) and away (e.g., to work) [16], while in the generated load profiles of Wepro-LPG, the fluctuations are obviously shown during the morning period which might be caused by the detailed behavioural model that emphasised on the person's desire developed in LPG model. Although, it requires a future analysis. In general, the Wepro-ALPG has more aligned curve shape with the average standard Dutch residential load profile as illustrated in Figure 14 than the Wepro-LPG, where it could be because ALPG model is built based on Dutch dwelling setting. Moreover, a measured dataset that adequately represents the case study is required for validation purposes although in general our generated hourly average load profiles have similar curve's trend with the standard Dutch residential load profile provided by NEDU.
In addition, our model is found to be more efficient in respect of its computational time. In this case processing the load profiles of the city's residential sector, which consists of a large number of households is more efficient rather than generating each household in LPG or a certain number of city's households in ALPG. It takes about ten minutes computation to generate a single-person household load profile in LPG and about fifteen minutes computation to generate a multi-occupants household load profile in LPG, for instance profile: Family, 3 children, both adults at work. Thus, in takes 60-min to generate the Wepro's selected five profiles of Table 3 in LPG, where we used LPG version 8.9.0. Furthermore, the simulation of the current configuration that consists of five households from four types of pre-defined profiles in ALPG takes about eight minutes. We use Python 3.7 (64-bit) to run this configuration. All of these simulations either LPG or ALPG were conducted on a computer using an Intel core i5-5300U CPU processor @2.3 GHz and 8 GB of installed memory (RAM). Thus, the computation will take much longer than our approach to generate a single or the few load profiles at the city level. In this case, our approach to model the residential sector at the city level has also tackle the limitation addressed by the ALPG that the tool is aimed at small group of houses which is maximum about 100 households per simulation. Consequently, our approach also creates efficiencies in the size and storage of the generated files. For instance, the output folder of one "single with work" profile generated in LPG has 2.6 GB size and the output folder of our selected pre-defined profiles in ALPG has 1.5 GB size.

Conclusions
This work has developed a simplified and practical approach to model residential electricity load profiles where the model can match the main city's characteristics with the representative pre-defined households profiles proportionally. The Wepro model is advantageous as an efficient approach to develop the residential electricity load profiles at the city level, especially where survey data, smart-meter data or any other local temporal profiles dataset are unavailable. The findings concur with some load profile studies from the similar climate profile which indicate Winter as the highest consumption share and illustrate either December or January is having the highest consumption share. The results of the selected typical days for all load profiles indicate that most weekend days have a higher consumption share than weekdays in the same season. Moreover, all the hourly average load profiles in a year and per season demonstrate the consistent curve shapes, demand peaks and the peak hours within the season and model, either using Wepro-LPG or Wepro-ALPG. In terms of the curve shape and daylight characteristics between the models, the hourly average in a year of Wepro-ALPG is preferred to be used because it also shows a high similarity with the shape of the standard Dutch household provided by NEDU or previously EDSN, although the Wepro-ALPG load profiles illustrate less morning fluctuations as a result of the few input profiles taken by the approach. In addition, in terms of the evening peak, the hourly average in a year of Wepro-LPG is preferable to be used, because it resembles the evening peak time of the Dutch household characteristics, where the evening peak takes place after dinner time, which concurs with a Dutch load profile study that the evening peak takes place after dinnertime when e.g., TV, dishwasher, etc., are on because within the average Dutch household, cooking is done using gas instead of electricity.
Moreover, our work contributes by evaluating the characteristics of residential electricity load profiles based on time variation analyses: seasonal analysis, monthly analysis, days analysis and hourly analysis. In addition, this method is applicable to model previous year, current year and future year, where for current year and future year are used city's projected numbers.
Furthermore, the few selected household profiles which are the representative of the city's profile may dominate the shape of the output profiles where all of input have represented the city's age group, labour force composition and gender share. Although the few selected profiles may dominantly influence the output profile, based on the results, they still resemble the Dutch average household profile and concur with the common peak demands characteristics. In addition, although the Wepro model depends on external household profile generators such as LPG and ALPG, the Wepro model is found to be more efficient in storage capacity and computational process of the residential sector's load profiles, given the number of households in the city that can represent the local profile.
In future work, it would be interesting to identify the potential of energy savings based on the generated load profiles using a relevant machine-learning technique. We also look forward to add more main input parameter to the model and compare with the case study's measured data. Further work might also be conducted to extend residential electricity temporal profiles into spatial profiles.
Author Contributions: The idea, method and analysis of the study is designed by A.K. A.K., P.D.K.M. and P.S.N. wrote the paper. P.D.K.M. performed the complex pre-processing and modelling tasks in Python. P.S.N. reviewed and proofread the paper. All authors have read and agreed to the published version of the manuscript.