A Spatial-Based Integration Model for Regional Scale Solar Energy Technical Potential

: One of the main objectives of human society in the present century is to achieve clean and sustainable energy through utilization of renewable energy sources (RESs). In this paper, the main purpose is to identify the locations that are suitable for solar energy in the Kurdistan province of Iran. Initially, solar-related data are collected, and suitable criterion and assessment methods are chosen according to the available data. Then, the theoretical potential of solar energy is assessed and the solar radiation map is prepared. Moreover, the technical potential of various solar technologies is evaluated in that study area. These technologies include concentrating solar power (CSP) and photovoltaic (PV) in power plant applications, and rooftop PV panels and solar water heaters in general applications. The results show that the Kurdistan province has the potential capacity for 691 MW of solar photovoltaic power plants and 645 MW of CSP plants. In the case of using solar water heaters, 283 million cubic meters of natural gas and 1.2 million liters of gasoline could be saved in fuel consumption. The savings in the application of domestic PV will be 10.2 MW in power generation.


Introduction
Human society has experienced great development in recent years, but this growing development has been associated with high energy consumption, which has been mainly met by fossil fuels [1]. According to statistics from the International Energy Agency (IEA), world energy consumption in 2014 was equal to 9425 million tons of oil equivalent (Mtoe). However, the share of fossil fuels in the production of this energy was 84.5%. In addition, 32,381 million tons (Mt) of CO 2 were produced in this year, which has been nearly doubled compared to the level reported in 1973 [2]. These statistics show the importance of preventing the emission of greenhouse gasses and, consequently, the importance of finding appropriate solutions, such as optimization, increasing the efficiency of current energy systems, and moving towards the use of renewable energy sources (RESs) [3]. One of the main challenges of RESs is their high dependence on time and location. Accordingly, it is necessary to fully investigate renewable energies geographically and their distribution in different spatial and temporal scales [4].
Solar energy is one of the renewable energies that has attracted considerable attention in recent years. So far, several studies have been conducted in different countries on the solar energy resource Sarvabad, Saqqez, Sanandaj, Qorveh, Kamyaran, and Marivan. Based on the 2016 population and housing census, the province has a population of about 1.6 million-71% urban and 29% are rural dwellers. The relative density of the population is 51.2 people per km 2 . The location map of the study area shown in Error! Reference source not found..
In terms of climate, Kurdistan is a mountainous region, which has extensive plains and strewn valleys across the area. The lowest elevation is 900 meters above sea level (m.a.s.l) and highest is 3300 m.a.s.l, this height difference in elevation makes the creation of different climates. This province has been chosen as a case study and an example in the country and so the methods of research can be applied to other provinces of Iran too.

Methodology
In this paper, a GIS-based evaluation and integration method [29,30] is used for assessing the solar resource in Kurdistan. GIS is a powerful computer tool that enables the user to record, store, and organize databases, and analyze spatial and geographical data in different coordinate systems. Due to the nature of renewable energy, it is necessary to take several geographical analyses in their resources assessment [31]. In this paper, different data for the evaluation of solar resources in Kurdistan were obtained mostly from the Renewable Energy and Energy Efficiency Organization of Iran (Satba) [32] and also other national and International agencies and organizations. The solar irradiances were recorded in more than 150 stations in Iran during last 10 years by Satba and many other local power companies.
For solar resource assessment and in general RESs assessment, two types of potentials are usually considered-the geographical (theoretical) potential and technical potential. The technical potential is derived from the conversion potential, taking into account additional restrictions regarding the area that is realistically available for energy generation, while the theoretical potential identifies the physical upper limit of the energy available from a certain source [8].
The theoretical potential of solar energy is defined as the total amount of annual solar radiation in suitable areas for solar applications, which should take into account multiple constraints in the assessment phase to achieve the appropriate areas [8]. The technical potential is defined as the amount of the total theoretical potential that can be converted into electricity using existing technologies [8]. However, the technical potential can be limited more and more by economic assessments and evaluating competitiveness with non-renewable technologies [33]. In this paper, both the theoretical and technical potentials of the solar resources have been calculated.
In the solar power plant applications, solar radiation map in the province is prepared using data from meteorological stations, sunshine duration, and the Solar Radiation tool [34]. The technical, economic, environmental, and geographical constraints are used to eliminate unsuitable areas and identifying suitable sites for solar applications. According to the remaining locations and assessing In terms of climate, Kurdistan is a mountainous region, which has extensive plains and strewn valleys across the area. The lowest elevation is 900 meters above sea level (m.a.s.l) and highest is 3300 m.a.s.l, this height difference in elevation makes the creation of different climates. This province has been chosen as a case study and an example in the country and so the methods of research can be applied to other provinces of Iran too.

Methodology
In this paper, a GIS-based evaluation and integration method [29,30] is used for assessing the solar resource in Kurdistan. GIS is a powerful computer tool that enables the user to record, store, and organize databases, and analyze spatial and geographical data in different coordinate systems. Due to the nature of renewable energy, it is necessary to take several geographical analyses in their resources assessment [31]. In this paper, different data for the evaluation of solar resources in Kurdistan were obtained mostly from the Renewable Energy and Energy Efficiency Organization of Iran (Satba) [32] and also other national and International agencies and organizations. The solar irradiances were recorded in more than 150 stations in Iran during last 10 years by Satba and many other local power companies.
For solar resource assessment and in general RESs assessment, two types of potentials are usually considered-the geographical (theoretical) potential and technical potential. The technical potential is derived from the conversion potential, taking into account additional restrictions regarding the area that is realistically available for energy generation, while the theoretical potential identifies the physical upper limit of the energy available from a certain source [8].
The theoretical potential of solar energy is defined as the total amount of annual solar radiation in suitable areas for solar applications, which should take into account multiple constraints in the assessment phase to achieve the appropriate areas [8]. The technical potential is defined as the amount of the total theoretical potential that can be converted into electricity using existing technologies [8]. However, the technical potential can be limited more and more by economic assessments and evaluating competitiveness with non-renewable technologies [33]. In this paper, both the theoretical and technical potentials of the solar resources have been calculated.
In the solar power plant applications, solar radiation map in the province is prepared using data from meteorological stations, sunshine duration, and the Solar Radiation tool [34]. The technical, economic, environmental, and geographical constraints are used to eliminate unsuitable areas and identifying suitable sites for solar applications. According to the remaining locations and assessing the amount of producible power from CSP and PV power plants, the capacity for solar power plants in the study area can be calculated. For the general application of solar power, the available roof area for this application is calculated in residential buildings. According to this area, the production capacity of domestic PV can be assessed. The effect of using solar water heaters in fuel consumption reduction can be calculated by computing fuel consumption for hot water.

Solar Resource Assessment
In calculating the solar energy potential, according to the available technologies, it is necessary to determine which one of the technologies will be applicable. Therefore, the following technologies are selected among the existing technologies because they have been widely used worldwide and have been employed in various provinces of Iran. Rooftop PV panels Solar water heaters In the power plant applications, it is necessary to obtain the solar radiation map in the whole province. Therefore, the next section describes how to create the province's solar radiation map.

Solar Radiation Map
Radiation data can be obtained in several ways. Generally, these methods can be divided into two main categories:

•
The satellite data • Ground measured data At present, various meteorological satellites are active in the atmosphere, but due to lack of access, unfortunately, in this study, it is not possible to use the satellite data. The ground data is commonly recorded by weather stations or radiation survey meters such as pyrometers. In all meteorological stations, which are located in the Kurdistan province, the solar radiation is not recorded and archived. The only data that are recorded for a long time, and are associated with radiation, are the number of sunny hours per day, available at all stations in the province-considering that data of the sunny hours per day are available in the stations. Recorded data are used to extract the radiation data in the stations. For this purpose, an efficient Angstrom-Prescott model [35] is used, which has long been used in solar radiation research, and its efficiency has been proven. This model is a statistical model with two constant factors. Data for Sanandaj (the most reliable data) in recent years is used to determine the constants of the model. In summary, it can be said the Angstrom model uses the Equation (1): where H is actual solar radiation on a horizontal surface that has reached to ground. H 0 is solar radiation outside the atmosphere at the same point. Both are in W/m 2 . n is the number of sunny hours that is measured in the station, and N is the total number of hours per day. The Angstrom model is applied to Sanandaj station (airport), and the obtained coefficients will be applied to the other stations' data. H/H 0 in the vertical axis and n/N in the horizontal axis (average per month as one point) is drawn. This diagram is shown in Figure 2. radiation in all stations in the province can be calculated accordingly, as reported in Table 1. Similar results obtained in other studies [32,36] indicate the effectiveness and accuracy of the results obtained from the Angstrom method in this paper [37].  In the next step, it is necessary to consider the topographic conditions in providing a radiation map. By using the digital elevation model (DEM) [38] and Solar Radiation tool and by applying clearness and diffuse coefficients 0.5 and 0.3, respectively, using clearness index from Homer Energy [39] (in this map it is assumed that these coefficients are constant in the whole province) the radiation map can be obtained. Radiation at the stations is calculated from this map and radiation obtained from the Angstrom model is divided on it and the resulting number will be called the coefficient of each station. The value of this coefficient is interpolated in the whole area and the resulted map is called the factor map. By multiplying the factor map in the radiation map obtained from the Solar Radiation tool the radiation map of the province is achieved. DEM map and the final map of solar radiation for the whole province are shown Error! Reference source not found. and 4, respectively. As seen in Figure 1, coefficients "a" and "b" are 0.2798 and 0.4315, with standard errors of 2.52% and 3.61%, respectively, in Sanandaj station, and the R 2 is 0.6413. These factors are applied to the other stations in the province to calculate radiation value from the sunshine duration. Average radiation in all stations in the province can be calculated accordingly, as reported in Table 1. Similar results obtained in other studies [32,36] indicate the effectiveness and accuracy of the results obtained from the Angstrom method in this paper [37]. In the next step, it is necessary to consider the topographic conditions in providing a radiation map. By using the digital elevation model (DEM) [38] and Solar Radiation tool and by applying clearness and diffuse coefficients 0.5 and 0.3, respectively, using clearness index from Homer Energy [39] (in this map it is assumed that these coefficients are constant in the whole province) the radiation map can be obtained. Radiation at the stations is calculated from this map and radiation obtained from the Angstrom model is divided on it and the resulting number will be called the coefficient of each station. The value of this coefficient is interpolated in the whole area and the resulted map is called the factor map. By multiplying the factor map in the radiation map obtained from the Solar Radiation tool the radiation map of the province is achieved. DEM map and the final map of solar radiation for the whole province are shown Figures 3 and 4, respectively.

Site Selection
The main objective of the solar resource assessment is the evaluation of radiation values, as well as finding potential areas where solar power generation is possible with available technologies. In this case, according to the accuracy of the available data, narrowing the entire area to appropriate areas continues until the suitable areas are finally identified. The assessment process is based on a systematic survey process. At each step, less attractive areas are removed and the process continues on promising areas. Data and information management in the resource assessment process requires the integration and interpretation of results to identify suitable areas. A huge amount of technical and non-technical information is needed to identify suitable areas. The restrictive analytical method [40] is used to identify potentially suitable locations. In this regard, the study area can be divided into two parts based on the possibility and unfeasibility of applying solar technologies. The Boolean logic (restrictive) method uses pre-defined constraints (as shown in Figure 4) for identifying the suitable location for placing the solar technologies. Boolean logic uses a binary condition for input and checks a binary condition for the outputs. Logical math tools consider the value 0 for false conditions and the value of 1 for the true conditions. With respect to each criterion, the study area is divided into two discrete classes: 1 for areas with the possibility of constructing solar power plants, and the 0 for unfeasibility of the areas [41]. The main steps in the restrictive method are as follows: • creating a conceptual model; • determining and localization of the desired criterion for site selection; • collecting the required data; • assessment of the study area; • identification of promising and probable areas based on desired criteria for each data layer; • using a data integration method based on the conceptual model; • determining the suitable areas for the construction of the solar power plant (site selection).
A flow diagram (conceptual model) of data intergration method using restrictive data layers to select potential solar power plants sites is shown in Figure 5. Based on this integration method 12, data layers are applied to evaluate the study area for finding and defining suitable areas. The restrictive method's data layers are divided into four main groups: technical, economic, environmental, and geographical constraints. Previous studies have used many layers to apply their criteria and constraints. These indicators have been assessed and localized for Iran, according to national and local laws. Criteria and constraints related to solar site selection are given in Table 2.
By considering the economic indicators and given the need to create a temporary road from project sites to transport links, the distance to these links should not exceed a specified limit because it will increase operational costs. In this research, this distance is intended for up to 10 km.
Given the need to avoid comfort disturbance during construction, as well as the need to keep the industrial areas away from population centers, distance from residential areas should be observed. A 2000 m buffer zone for the cities and 1000 m for the villages are considered in this study.
Consideration of 250 m buffer zone for the transport links is due to property rights that are considered in the road property rights laws. Moreover, according to the power lines property rights laws, there should be a required minimum safety corridor around power lines with consideration to safety clearances. High-voltage transmission lines at different voltages have different limits, where the highest degree is taken into account here, which corresponds to 750 kV lines.
For wetlands, coastlines, forests, and faults a 500 m buffer zone is considered, while this is around 2000 m for environmental protected areas and 700 m for historical and cultural sites.  By applying these economic and environmental restrictions on radiation maps, unsuitable areas for the construction of solar power plants will be removed from the primary radiation map. A series of regions are also omitted due to technical reasons. According to the analysis carried out in [47], a minimum amount of radiation is required (here 5 kWh/m 2 /d) at a location to obtain power in an acceptable range from solar power plants. Accordingly, some parts are removed as well, and so radiation maps will be limited even more.
The photovoltaic panels with a nominal capacity of 1 kW, occupies an area of about 10 m 2 . In order to install a power plant with 2 MW capacity, an area of over 20,000 m 2 is required [47]. This also applies in relation to CSP. From the perspective of geographical constraints in solar power plants   [47] By applying these economic and environmental restrictions on radiation maps, unsuitable areas for the construction of solar power plants will be removed from the primary radiation map. A series of regions are also omitted due to technical reasons. According to the analysis carried out in [47], a minimum amount of radiation is required (here 5 kWh/m 2 /d) at a location to obtain power in an acceptable range from solar power plants. Accordingly, some parts are removed as well, and so radiation maps will be limited even more.
The photovoltaic panels with a nominal capacity of 1 kW, occupies an area of about 10 m 2 . In order to install a power plant with 2 MW capacity, an area of over 20,000 m 2 is required [47]. This also applies in relation to CSP. From the perspective of geographical constraints in solar power plants due to the construction costs and also for receiving maximum sunlight, the slope should be less than 3%.
As appears from the conceptual model in Figure 4, due to technical, economic, environmental, and geographical constraints with addressing the entire area of study, logical value for the areas that are not suitable for the construction of solar plants is zero, and these areas will be removed. Finally, all the appropriate technical, economic, environmental, and geographical layers are integrated to obtain final proper areas. It can be seen that the final areas are selected as suitable sites that satisfy all constraints, and in all criteria have a logical value of 1.
By applying economic and environmental constraints, suitable areas map for solar power plant construction will be as Figure 6. Province slope map was created using DEM maps. Geographically permissible area map (areas with slopes less than 3%) is shown in Figure 7. Finally, by applying the technical constraints, and the integration of four layers of suitable areas, in terms of technical, economic, environmental, and geographical constraints. Figure 8 is achieved as the final map of suitable areas for solar power plant applications. According to the map of Figure 8, it can be seen that an area of about 62 km 2 is available in the province for the construction of solar power plants.
Sustainability 2020, 12, x FOR PEER REVIEW 9 of 19 due to the construction costs and also for receiving maximum sunlight, the slope should be less than 3%.
As appears from the conceptual model in Figure 4, due to technical, economic, environmental, and geographical constraints with addressing the entire area of study, logical value for the areas that are not suitable for the construction of solar plants is zero, and these areas will be removed. Finally, all the appropriate technical, economic, environmental, and geographical layers are integrated to obtain final proper areas. It can be seen that the final areas are selected as suitable sites that satisfy all constraints, and in all criteria have a logical value of 1.
By applying economic and environmental constraints, suitable areas map for solar power plant construction will be as Error! Reference source not found.. Province slope map was created using DEM maps. Geographically permissible area map (areas with slopes less than 3%) is shown in Error! Reference source not found.. Finally, by applying the technical constraints, and the integration of four layers of suitable areas, in terms of technical, economic, environmental, and geographical constraints. Error! Reference source not found. is achieved as the final map of suitable areas for solar power plant applications. According to the map of Error! Reference source not found., it can be seen that an area of about 62 km 2 is available in the province for the construction of solar power plants.

CSP Plants
As mentioned, in this study two application of solar power plant is intended. The first case relates to a CSP system with the following characteristics [47]:

CSP Plants
As mentioned, in this study two application of solar power plant is intended. The first case relates to a CSP system with the following characteristics [47]: • dry-cooled; • six hours of storage;

CSP Plants
As mentioned, in this study two application of solar power plant is intended. The first case relates to a CSP system with the following characteristics [47]: • dry-cooled; • six hours of storage; • the amount of solar multiple equal to 2; • the power density of 32.895 MW/km 2 .
The performance of the system will be different for different environmental conditions. According to an analysis done with the System Advisor Model (SAM), various coefficients values (called capacity factors) are obtained for different radiation values [47]. By categorizing solar radiation to different classes, Table 3 is achieved. Finally, the potential of solar energy (in MWh) in the CSP system will be obtained from the following Equation.
where A is the available area for the construction of solar power in each class. P max is the maximum power produced by this system, which is equal to 32.895 MW/km 2 . h yr is the number of hours in a year, which is equal to 8760, z is the capacity factor in each class and is obtained from Table 3. Therefore, from Equation (2), the potential of solar energy in each class will be achieved, and the total potential of different classes will be equal to the potential of this application. The amount of radiation at selected suitable locations for solar power plant applications (Figure 7) are placed in class 1 of Table 3. Therefore, the capacity factor (z) will be equal to 0.315. According to the map of Figure 7, the amount of available area (A) is equivalent to 62.289 km 2 . By inserting these numbers in Equation (2), the evaluated power in the application of CSP plants will be equal to 645 MW.

PV Farms
PV power plants usually use a lot of same panels in a wide range of ground, so each panel reveals the function of the whole system. The considered system, which its analysis is available in [47] has the following characteristics: • single-axis tracking, with a north-south axis of rotation; • the slope angle of zero degrees; • maximum power of 48 MW/km 2 .
These systems will create a full-function only in designed conditions. According to the analysis conducted by NREL, the performance of these systems under various conditions is a coefficient of maximum power. These coefficients in different climatic conditions, and in various US states, have been presented in [47]. The capacity factor of 0.231 is considered for Kurdistan province 3030. Therefore, the power generated by the small PV power plant systems, in Kurdistan province is calculated from Equation (2), and with a capacity factor of (Z) 0.231, maximum power of 48 MW/km 2 , and an available area of 62.289 km 2 . By inserting these numbers in Equation (2), the evaluated power in the solar PV power plants application will be equal to 691 MW. This is the potential power of the solar PV power plants in Kurdistan province.
In the final evaluation of solar power plants in Kurdistan province from PV and CSP plants, it should be mentioned that solar power plant potential is not the sum of these two numbers, because if you use one of them, another cannot be used, and this amount depends on the amount of the use of these two applications based on desired scenarios.

Residential Applications
Another application in the use of solar energy is related to small-scale applications. In this regard, two applications are considered in this study: Rooftop PV panels and solar water heaters. In the residential areas, the installation location of this equipment is the flat and non-shadow surface of rooftops, which is obtained as fallows.

Available Area for Residential Applications
In accordance with the Master Plan maps of urban areas, and after analyzing in the ArcMap software, and separation of residential buildings, the surface areas of these sites will be achieved ( Figure 9). Unfortunately, Master Plan maps in Marivan and Sarvabad cities, and urban areas that are not located in the center of the counties were not available. Thus, from the population in cities with Master Plan maps, we tried to estimate the other cities' residential areas. Thus, by calculating the per capita residential areas in cities with maps, and their average, and then by multiplying this average by the population of cities without Master Plan map, the residential area can be calculated in these cities.

Residential Applications
Another application in the use of solar energy is related to small-scale applications. In this regard, two applications are considered in this study: Rooftop PV panels and solar water heaters. In the residential areas, the installation location of this equipment is the flat and non-shadow surface of rooftops, which is obtained as fallows.

Available Area for Residential Applications
In accordance with the Master Plan maps of urban areas, and after analyzing in the ArcMap software, and separation of residential buildings, the surface areas of these sites will be achieved (Error! Reference source not found.). Unfortunately, Master Plan maps in Marivan and Sarvabad cities, and urban areas that are not located in the center of the counties were not available. Thus, from the population in cities with Master Plan maps, we tried to estimate the other cities' residential areas. Thus, by calculating the per capita residential areas in cities with maps, and their average, and then by multiplying this average by the population of cities without Master Plan map, the residential area can be calculated in these cities. In rural areas, a roof area of 100 m 2 can be considered for each family [48]. According to the statistics of rural households in each city, we can also calculate the rural residential area. Population, occupied area by residents and residential area per capita in cities with Master Plan in Kurdistan province, is presented in Table 4.
The average residential area is 27.79 m 2 /person. Using this average, the residential area is calculated in other urban areas (urban areas without a Master Plan map). By gathering the urban areas of each county together, the urban residential areas in each county will be achieved (the second column of Table 5). Furthermore, with regard to the consideration of 100 m 2 for each rural household, the approximated rural residential area can be obtained for each county (the third column of Table 5) and by summing these two values, the entire residential area in the counties can be achieved (fourth column of Table 5). In rural areas, a roof area of 100 m 2 can be considered for each family [48]. According to the statistics of rural households in each city, we can also calculate the rural residential area. Population, occupied area by residents and residential area per capita in cities with Master Plan in Kurdistan province, is presented in Table 4.
The average residential area is 27.79 m 2 /person. Using this average, the residential area is calculated in other urban areas (urban areas without a Master Plan map). By gathering the urban areas of each county together, the urban residential areas in each county will be achieved (the second column of Table 5). Furthermore, with regard to the consideration of 100 m 2 for each rural household, the approximated rural residential area can be obtained for each county (the third column of Table 5) and by summing these two values, the entire residential area in the counties can be achieved (fourth column of Table 5). According to the analysis conducted by NREL [49], it can be assumed that in lands with residential applications, 8% of the area is flat and on the roof; thus, can be considered as a potential for installing solar equipment. In addition, other factors, such as shading, roof structure, and so on will restrict this amount. According to this analysis, in the cold conditions (similar to Kurdistan province) 22% of the previous value is appropriate for using solar energy [49]. Therefore, available areas for solar energy in general applications in each county, obtained by multiplying the coefficients of 0.08 and 0.22 in the residential area of the counties, which is shown in the 5 th column of Table 5. After identifying this area, we must determine the potential solar energy production from the two above-mentioned applications.

Rooftop PV Panels
In the use of solar PV panels, the density of power generation in the area is considered around 110 W/m2. On average, the efficiency of small-scale PV systems can also be considered as 13.5% [47]. So by multiplying these two amounts in the available area, the potential of the producible power of photovoltaic panels can be obtained. This amount is presented in the last column of Table 5, for each county. This amount is more than 10.2 MW in total for Kurdistan province.

Solar Water Heater
In Iran, natural gas and diesel are mainly used for water heating. Therefore, by having gas and diesel consumption in the residential sector in each county, the amount of fuel needed to generate hot water can be achieved. For this purpose, the following assumptions are made in this study.

•
About 23% of fuel consumption in the residential sector is for hot water generation [25]; • solar water heater systems will reduce about 50% of fuel consumption for hot water [25]; • the efficiency of water heating systems based on gas and diesel is expected to be 50% [50]; • the energy content of natural gas and diesel used in residential sector is 37.2 MJ/m 3 and 36.6 MJ/l respectively.
According to the above cases, we can consider the following equations: where, V 1 is the volume of natural gas consumed in the residential sector (m 3 ), V 2 is the volume of diesel consumed in the residential sector (L), and V 11 is energy from the consumed natural gas (MJ), V 21 is energy from consumed diesel (MJ) and P is the potential of solar energy for domestic hot water supply (MJ). Hence, the potential for energy saving will be achieved in each county. Statistics related to fuel consumption has been extracted from 3737, and the percentage of consumption in the residential section for 2011 has been extracted from 2727. Iran's energy balance report for 2011 indicates that 70.18% of natural gas and 1.23% of diesel consumption is in the residential sector (as well as gas and diesel consumption in each county of the province is shown in Table 6, see the second and third columns). Given that solar water heaters will reduce about 23% of fuel consumption in the residential sector, natural gas and gasoline savings potentials by application of solar water heater in each county have been computed and are shown in Table 6, (fourth and fifth columns). Finally, the amount of energy saving in fuel consumption, as a result of using solar water heaters in each county can also be calculated according to Equation (5). These results are shown in Table 6, in the sixth and seventh columns, separately for each county.

Evaluation of Total Solar Energy Potential
As mentioned, the available area will be for the use of both systems (not simultaneously); and so it cannot be represented total solar energy potential in general application. This will be possible by considering more scenarios, close to reality forecasts, and more detailed assessment. Therefore, solar potential map in these two general applications is presented separately for each application. The results of this section are shown in the map of Figure 10 as a potential use of solar water heater and Figure 11 shows domestic PV potential. The results of this section are shown in the map of Error! Reference source not found. as a potential use of solar water heater and Error! Reference source not found. shows domestic PV potential.   The results of this section are shown in the map of Error! Reference source not found. as a potential use of solar water heater and Error! Reference source not found. shows domestic PV potential.   The results show that in terms of the potential use of solar water heater, Sanandaj is the highest and Sarvabad is the lowest priority. The same is true for domestic PV potential. In fact, in the general applications of solar potential, there is a direct proportion to the population. These results also show that fuel consumption saving by using solar water heaters is 1.2 million liters for diesel and 282.8 million cubic meters of natural gas. As well as by implementation of the domestic PV systems, power consumption is reduced more than 10 MW. The technical potential of solar energy in Kurdistan province, using various solar technologies are summarized in Table 7. As can be seen in the table, the results indicated that the Kurdistan province has capability to produce 691 MW of power from PV panels, and 645 MW from CSP plant, which are more than total power demand of the province. It means that by developing full capacity of solar potential in the Kudrestan total power demand can be supplied in the future. In the case of using solar energy for water heating, totally, 283 million cubic meters of natural gas and 1.2 million liters of gasoline could be saved. By developing solar energy, system energy in the area the energy supply security could be increased in wintertime when the high snowfall blocks the roads for fuel transportation.

Conclusions
In this paper, solar energy resources in Iran's Kurdistan province were identified and evaluated. At first, solar resource-related data were collected in Kurdistan province and an appropriate assessment method was selected based on available data. Then, the theoretical potential of the solar energy in the province was evaluated and solar energy atlas was prepared. In the following, the technical potential of various solar technologies was studied. These technologies include CSP plant and PV farms in power plant applications, and rooftop PV panels and solar water heaters in general domestic applications.
In solar power plant applications, primarily using data from meteorological stations, sunlight duration, and Solar Radiation tool, radiation map was prepared in Kurdistan province. Then, using the technical, economic, environmental, and geographical constraints, some areas were discarded and finally, the suitable areas were identified. The available technical potential in the use of solar power plants was evaluated by considering the suitable areas and the recoverable amount of power in the application of PV and CSP plants. For general application, the available area in residential rooftops was calculated. Finally, with respect to this area, the production capacity of domestic PV was evaluated. The effect of solar water heater implementation in gas and diesel consumption value was calculated using fuel consumption for hot water and resulted in saving in the use of this application.
The results show that the Kurdistan province is capable of achieving, 691 MW of PV power generation, and 645 MW of CSP plant. In the case of using solar water heaters, 283 million cubic meters of natural gas and 1.2 million liters of gasoline could be saved. The savings in the application of domestic PV panels will be 10.2 MW in power generation.
Future works will be dedicated to the assessment of the potential of other renewable energy sources (RESs), comprehensive evaluation of the use of RESs, and finally, evaluation of different scenarios to increase the share of these sources in the consumption pattern of the Kurdistan province.