Techno-economic analysis of wind power integrated with both compressed air energy storage (CAES) and biomass gasification energy storage (BGES) for power generation

A techno-economic analysis of excess wind electricity powered adiabatic compressed air energy storage (A-CAES) and biomass gasification energy storage (BGES) for electricity generation is implemented to determine the performance of the system and the potential profitability of developing such a facility for distributed power generation in the UK by an investor, given the customer's demand for heat and electricity. The customers are considered to be about 1600 households in the Humber region, UK, who use heat generally for space heating and domestic hot water applications. The system is modelled using a developed Matlab computer code and its performance evaluated using total system efficiency (TSE), net present value (NPV) and cost of electricity (COE) as metrics. TSE of 36.8% is obtained for the system while the COE is found to be about £0.19 per kW h. In terms of profitability, the system returned a negative NPV of £2 144 062 signalling the non-profitability of the system in the proposed location. However, if 70% of total investment cost (TIC) of the system is provided for by means of a subsidy, the system becomes economically viable with positive NPV of £132 475 and COE of £0.10 per kW h respectively. The sensitivity study shows that the most significant factors swaying the NPV of the A-CAES-BMGES are TIC, O&M cost, excess wind electricity cost, electricity tariff and cost of diesel fuel.


Introduction
Global energy use is dominated by fossil fuelled sources. 1,2 In the recent past, the issues of climate change, depleting nature of fossil fuels and the yearly rise in global energy demand has led to a new global energy policy focusing on improvement of energy efficiency, promotion of energy conservation and an energy mix involving increased deployment of renewable energy (RE) sources [1][2][3] Wind power generation is one of the RE sources that is continuously becoming attractive for meeting current and future energy needs. 1,2 This is due to its established technology, wide availability, good scalability and relatively low cost with projected cost of electricity (COE) below 34 $ per MW h for onshore wind compared to 54 $ per MW h for large, groundmounted PV in the United States at 3% discount factor. 4 In 2016, wind power accounted for more than half of 14.1% global renewable energy use growth. 5 It has a signicant role to play in the quest to decarbonise the globe and switch from fossil based generation to renewable based generation. Wind is forecasted to have a good signicant future cost reduction potential as its technology improves. 4 However, wind suffer from a major limitation. It is unpredictable and can only provide energy intermittently. The intermittency leads to a mismatch between electricity demand and supply which consequently introduces an enormous challenge of safety and reliability to the grid operation. In addition, the energy density of wind is low in comparison to that found in conventional fossil fuel sources like coal, oil, gasoline and natural gas. 6 Energy storage (ES) has been identied as a major solution to the intermittency of wind. This is because of its favourable charging and discharging features, which makes it possible to overcome the uctuations. 7 ES has the capacity to even out the mismatch between energy supply and energy demand, thus playing a critical role in energy conservation. It also improves the efficiency and reliability of energy system. The increased efficiency would result to energy conservation and increased cost effectiveness.
Different energy storage (ES) technologies exist: Compressed Air Energy Storage (CAES), Pumped Hydro Storage (PHS), Battery Energy Storage (BES), Super Magnetic Energy Storage (SMES) and Flywheel Energy Storage (FES). 1,4,[8][9][10][11][12][13][14][15] Though each of the ES technologies offer exclusive advantages, they also suffer from special disadvantages. CAES is one of the mature and economically attractive ES technologies with capacity to enable a more efficient and exible energy system with a better use of the uctuating wind and other renewable energy sources. 3,16 In a CAES process, during the charge time at low-cost offpeak periods, ambient air is compressed and stored in large underground caverns or above ground storage vessels using excess electricity. During the discharge time at peak period, the stored compressed air is preheated and expanded through a turbine to produce electricity. 17 Two basic types of CAES system exist: diabatic CAES (D-CAES) and adiabatic CAES (A-CAES). 17 The D-CAES involves the combustion of fossil fuels with the compressed air (CA) to raise the temperature of the compressed air before expansion via gas turbines, leading to CO 2 emissions. The A-CAES stores the heat released during the compression process and reuse it to heat up the stored compressed air in the expansion process. 17 An alternative CAES system exists called the Isothermal CAES (I-CAES). Isothermal CAES (I-CAES) technology aims to address some of the diabatic and adiabatic CAES, by eradicating the need for fuel and high temperature heat energy storage, thus, offering an enhanced round trip efficiency in the range of 70-80%. [18][19][20] Being an isothermal process, the air is compressed without a change in temperature, thus reducing the work of compression while maximizing the work needed for expansion, through effective heat transfer with surroundings of the air vessel. An experimental study by Alami et al. 21 on low pressure, modular small scale compressed air energy storage (CAES) system for wind energy storage applications working on the I-CAEs principle found that an I-CAES CAES working with a series of off the shelf low pressure cylinders instead of pressurised containers or cavern could be preferable on smallmedium scale (1-10 MW), 21 and has the following advantages: enables the use cheaper, low pressure storage containers, and practically eliminates the need for heat removal considerations necessary in higher pressure systems to offset the temperature rise aer compression. Also, it allows close control of discharge rates according to power consumption needs while maintaining minimal losses. The authors also reported a maximum overall system efficiency of about 97.6% for the system, while its physical footprint is less than 0.6 m 3 .
However, while there is substantial research into nearisothermal compression for CAES (by companies like Lightsail and SustainX), 22 I-CAES is not yet commercially available and any currently available compression that approaches reversible isothermal compression is too slow for industrial use 20,23 due to the impractically small temperature differences required. 20 Consequently most reported A-CAES designs opt for a series of adiabatic or poly-tropic compressions, aer each of which the compressed air is cooled back to the ambient temperature in order to reduce the both the temperature and volume of the air. 22,24 There are different options for storing the thermal energy in a CAES process: sensible heat storage (SHS) in liquids and solids, latent heat storage (LHS) in phase change materials (PCMs) and thermochemical heat storage (THS) in chemical reactions. 17 Although ES can help solve the problem of variability of wind, it does not solve all the problems. Shortfalls resulting from losses in low roundtrip efficiencies of an ES system can be compensated by designing a hybrid ES systems involving the RE source with ES and other energy sources of generation that can be mobilised on short notice and whose generation features must complement each other. 25 For example, a combination of A-CAES and biomass gasication electrical storage (BMGES) system for meeting the power and heat requirements of a given community. In fact, wind powered A-CAES and BMGES system are complementary to some extent, since some of the heat from syngas produced during gasication can be used by the A-CAES process whereas the biomass gasier can use some of the cheap heat generated by the compression process of the A-CAES to dry its fuel before gasication. Moreover, combination of A-CAES and BMGES may act as a CHP system, by recovering heat from the dual fuel engine (DFE) during and aer the power generation phase. This can only be achieved by placing the A-CAES and BMGES close to the area of energy demand, given the technical difficulty to transfer thermal energy at large distances. 26 Such a plant allows to maximize and conjugate the extensive benets of distributed generation (DG) 11 with those of electricity storage.
Many hybrid energy storage systems involving A-CAES has been proposed and analysed in recent years: Garrison and Webber 27 proposed a novel hybrid conguration involving solar energy and excess wind electricity powered CAES system in which the air from the cavern is preheated by solar energy rather than natural gas prior to expansion in the turbine. They reported an overall efficiency of 46% for the coupled solar-CAES system. Zhang et al. 13 analysed a hybrid diesel DG system integrated with A-CAES and thermal energy storage (TES). Their results suggested that the hybrid system's exergy efficiency is 41.5%, and the primary fuel saving ratio is 23.13%. In another study by Singh and Baredar, 28 a techno-economic assessment of a solar PV, fuel cell, and biomass gasier hybrid energy system was reported. The cost of electricity (COE) and net present cost of the system was reported to be 0.23 $ per kW h (15.064 Rs per kW h) and $79 858.76 (Rs. 5 189 003) respectively.
A biomass-red combined cooling, heating and power system with TES system was studied by Caliano and Bianco. 29 Their result suggested that the combined use of a TES and cold TES during the hot season could represent a viable economical solution.
In another related study, Singh et al. 30 carried out optimal sizing and feasibility study of an island micro grid in rural area consisting of PV, wind, biomass and BES system using articial bee colony (ABC) algorithm. They reported the sizes of the system components as 250 kW solar PV, 19 kW wind turbines, 1400 batteries and 40 kW biomass gasier. 30 To verify the strength of the ABC technique, the results obtained were compared with that obtained from the standard soware tool, hybrid optimization model for electric renewable (HOMER) and particle swarm optimization (PSO) algorithm. It was reported that ABC provides an optimal conguration with least levelised COE of 0.173 $ per kW h.
A hybrid solar-biomass power plant without ES was studied by Srinivas and Reddy. 31 Their result indicated that the plant fuel energy efficiency increases from 16% to 29% with an increase in solar participation from 10% to 50% at the boiler pressure of 20 bar.
To provide a solution for managing excess heat production in tri-generation plant and thus increase the power plant annual efficiency, a hybrid optimization model of biomass trigeneration system combined with TES was studied by Dominković et al. 32 Three case studies with minimum yearly average power plant efficiencies of 50%, 65% and 75% respectively were conducted. It was reported that an increase in overall power plant efficiency from 50% to 65% in legislation, in order to be eligible for the maximal feed-in tariff, would not affect the economics of the system.
Hence, from the literatures reviewed so far and within the limits of the author's knowledge, no treatment exists in the literature currently that focused on the technical and economic analysis of wind power integrated with both Compressed Air Energy Storage (CAES) and Biomass Gasication Energy Storage (BGES) System for Power generation in the UK. Downdra gasier is suited for small and medium-sized applications. 33 This conversion technology has a fairly efficient biomass to gaseous fuel conversion and produce syngas with relatively low amounts of tar that is suitable for direct use in internal combustion engines. 33,34 Hence, a xed bed downdra gasier is selected for this work. Given the advantage of A-CAES system and the need to increase RE in the energy mix around the globe, it could be expected that the combination of wind powered A-CAES and BMGES hybrid system will be an important trend. Hence, in order to investigate this concept, we developed a mathematical model of A-CAES integrated with BMGES system in Matlab. The A-CAES is powered using off-peak electricity from wind turbine. Technical and economic analysis of the hybrid system is analysed using efficiency, Net Present Value (NPV) and Cost of Electricity (COE) as metrics. The impacts of some technical and economic factors on system efficiency and protability of the of hybrid system are subsequently investigated.

System description
The system being proposed is as shown in Fig. 1. It is made of the following components: A-CAES powered by excess wind electricity, Hot air dryer (HAD) for drying biomass; Biomass gasier (BMG) for generating syngas from dried wood, dual fuel engine (DFE) running on syngas and diesel fuel and heat exchangers (HXs).
First the excess electricity of the wind turbine powers the compression train to compress air to a higher temperature and pressure. The heat of compression of the rst stage air compressor (AC1) is recovered by the air compressor heat exchanger (HX1) for biomass drying in the HAD, while the heat of compression of the remaining three stages (AC2, AC3, AC4) is recovered and stored as TES 1 , TES 2 and TES 3 respectively using the phase changes materials (PCM) in TES tanks, TEST 1 , TEST 2 and TEST 3 respectively. The cooled compressed air (CA) is stored in the air store (AS).
Subsequently, wood chips at initial moisture content (MC i ) of about 50% wt. are dried in the HAD to nal moisture content (MC f ) of about 10% wt. The dried wood is converted into syngas (sg) in the BMG. The heat carried by the generated syngas in recovered by means of the heat exchanger (HX2) and used to pre heat the temperature of the air exiting the last stage air compressor (AC4). The syngas is further cooled by recovering its heat using a heat exchanger (HX3). The recovered heat is used for biomass drying. Cooling the syngas helps to prevent preignition, to improve the engine volumetric efficiency, and to ease gas clean-up. 33 Aerwards, the cooled syngas which is limited at a temperature of 60 C, to minimise condensation, is cleaned and stored in the syngas store (SgS). Syngas from a gasier contain various contaminants and thus has to undergo several cleaning steps to reduce the concentrations of contaminants to level required by the downstream processes. For use in internal combustion engine, typical syngas cleaning requirements is particulate (<50 mg m À3 ), tars (<100 mg m À3 ). 35 Therefore, wet syngas cleaning, was used in this study. In such a scheme, the syngas is rst passed through a cyclone to reduce particulates. The syngas is then further cooled to a temperature below 100 C and passed through a spray scrubber to remove tar and ranges of other contaminants. 36 When electricity is required, the stored syngas is used to generate electricity in the DFE while the stored cool CA at high pressure exiting the AS is rst throttled down to the minimum operating pressure of the AS (40 bar) using a throttle valve (TV) and then heated using Q TES3 . The high temperature pressurised air leaving Q TES3 is subsequently heated up further by the exhaust heat from the DFE by means of the heat exchangers (HX4), and then expanded in AE1. Aer expansion in AE1, the temperature of air exiting the rst stage AE (TET AE1 ) is heated up rst by Q TES2 and further again by HX5 before it undergoes expansion in AE2. The process of heating and further heating by Q TES and HX is repeated until the last stage AE to produce electricity. Heat is recovered from the cooling water jacket of the DFE using cooling water heat exchanger (HX7).
A techno-economic analysis of the proposed hybrid system has been implemented to determine its value to the utility system and the potential protability of developing such a facility by an investor, given the customers demand for heat and electricity. The customer is considered to be about 1600 households in the Humber region UK who use heat generally for space heating and domestic hot water applications. For brevity, the proposed hybrid system is hereaer referred to as (A-CAES + BMGES)

Demand proles
To simulate the proposed system performance for the assumed 1600 households, the electrical and heat demand proles were generated as shown in Fig. 2. The electrical prole was obtained from the data of half-hourly average daily electrical energy consumption in kW h for the period 01/05/2011 to 31/05/2012 for 5554 households from the UK Northern power grid region. 37 The data was manipulated to give hourly average data based on the average daily electricity consumption of 4115 KW h per household obtained from a DECC report. 38 The heat demand proles are based on real hourly average data from a 100 home community housing development in the UK, obtained from a study by Wood. 39 The data was manipulated to give electricity to heat ratio equal to the national average for UK social housing of 0.28. 39

Modelling of the A-CAES + BMGES
The modelling of the proposed system is presented hereunder. To simplify the analysis, the following assumptions is made: The system and its sub components in any operating mode attain steady state. 24,40 Pressure losses in heat exchangers and pipes is ignored. 24,40 There is no change in enthalpy during the throttling process of CA. 24,40 Temperature of the CA in the AS is assumed equal to ambient temperature. 24,40 Air and the DFE exhaust gas is treated as ideal gas. Space heating water enters HX7 at 35 C and leaves at 55 C. 12 The changes in kinetic and potential energy are neglected. Isentropic efficiencies of the compressor and turbine are 0.85 and 0.88 respectively. 13 The motor (generator) efficiency of AC (AE) are 0.99 and 0.97 respectively. 13 The heat exchangers (HXs) has effectiveness (3) of 0.7.

Plant capacity
The electrical capacity of the proposed system is sized to match the maximum electricity demand of 1300 MW. Any excess/ decit electricity/heat is sold/bought to/from the grid. The DFE is sized to supply 25% of the total power while the A-CAES system supplies the rest. This journal is © The Royal Society of Chemistry 2018

Gasier and dual fuel engine model
The design requirement is to generate a rated output (P e,DFE ) in the dual fuel mode using syngas mixed with diesel. The performance of the dual-fuelled engine is characterised by a factor called syngas fraction (sf). It indicates the amount of diesel fuel (df) input that is replaced by the syngas (sg) in the dual fuel mode. Diesel replacement up to 80% is possible. 33 The thermal power input of the syngas (P th,sg ) and diesel (P th,df ) to the dual-fuel engine is determined as follows: The expected feed rate of biomass on dry basis (d.b) to achieve P th,sg is calculated as follows:

Dryer
According to the HAD, as in Fig. 3, the mass balances in the dryer (D) can be written for the dry biomass (B), air (a) and water (w) as follows: where _ m is mass ow rate, 4 is the humidity ratio and subscripts a, D, B, i, e, w and ev denotes air, dryer, biomass, in, exit/nal, water and evaporation for specie involved respectively. The mass ow rate of air entering the HAD (ṁa ,D ) is obtained as follows: where ṁD ,AC and ṁD ,sg denote mass ow rate of air recovered from AC1 and HX3 respectively. The initial feed rate of wt wood chips to the dryer (ṁB ,i ) at initial moisture content (MC i ) on wet basis is derived from the feed rate of wood in the biomass gasier (ṁB G ) as follows: The heat required for drying (Q D ) was obtained from energy balance as follows: 41,42 Q ev is the heat required to evaporate water, Q l,D is heat loss from dryer structures and Q s is sensible heat required to heat the biomass and drying structure to drying temperature. where: Heat loss (Q l,D ) is expressed as 5% of the heat input (Q D ). The heat input to the dryer, Q D is supplied by the hot air as it is cooled from its initial state ( _ m a,D , T a,i , RH i ) to nal state (ṁa ,D , T a,e , RH e ). Thus the energy supplied by the hot air (HA) is as follows: where In most hot air dryers, the exit temperature of the dried solid is a few degrees lower than the nal temperature of the hot air. 43 In this model, the difference between T a,e and T B,e is set to be 5 C.
With known ṁa ,D determined from eqn (7), eqn (6) and (9) are solved to obtain the nal humidity of air at the exit of the HAD and the temperature of the exiting air and biomass respectively.
The fan electrical power required for circulation of air in the dryer is calculated as follows: j is the specic fan power. In this model j is set to 1.5 kW m À3 s À1 as found by Nilsson. 44

Dual fuel engine (DFE) exhaust ow
From the known fuel input energy, the ow rate of the syngas and diesel fuel is calculated as follows: The exhaust gas ow rate of the DFE is estimated using: Table 1 shows the HAD and DFE input design parameters. The relevant characteristic parameters of the selected diesel engine has been take from the manufacturer's specication. 45 3.5 A-CAES system 3.5.1 Compressor. Assuming maximum work per compressor stage, the mass ow rate of air in each compressor stage (mȦC) is calculated as follows: The pressure ratio of each compressor (b AC ) stage is a function of the number of stages (n s,AC ) and overall pressure ratio (b O,AC ) of the compressor. It is calculated using the following equation: The outlet temperature from the compressor is obtained using the relation: The heat recovered from AC1 for biomass drying (RH AC1 ) was obtained as follows: The mass ow rate of air from AC1 for biomass drying is obtained as follows: The syngas (sg) produced from a downdra BMG is drawn from the high temperature reduction zone operating at about 815-1200 C. 48 At the point of leaving the gasier, the temperature of the hot syngas (T sg,i ) lies in the range of about 500-816 C. 48,49 For a given heat exchanger effectiveness (h recup ), the air temperature leaving the Cc HX2 and entering the TEST is obtained as follows: where C min is min [ _ m a,AC C p,a , _ m sg , C p,sg ] The temperature of the syngas leaving the Cc HX (T sg,O ) is obtained by energy balance. As the temperature of the air exiting the compressor is somewhat high, all the heat contained in the syngas from the gasier is not fully recovered. If the temperature of the syngas exiting the HX2 is more than 200 C, the heat it contains is recovered further using syngas cooler HX3 for use in generating hot air at a temperature of about 80 C for biomass drying while the syngas is cooled to a nal temperature not less than 60 C The heat recovered from syngas cooler HX3 for biomass drying (RH D,sg ) was expressed as follows: The mass ow rate of air from HX3 for biomass drying is obtained as follows: 3.5.2 Thermal energy storage tank (TEST). The thermal storage capacity of the TEST is expressed as follows: PCM thermal storage media has been considered in the present analysis because among the various heat storage techniques, latent heat (LH) TES has the capacity to provide large heat storage capacity and isothermal behaviour during the charging and discharging processes. The other benets of LH TES include compactness in sizes and low weight per unit storage capacity. 50 To achieve a suitable TES, the highest melting point of the PCM in each TEST must match the required turbine inlet temperature for each stage (TIT sn ). The chosen melting point of PCM (T m ) in each TEST is related to the entering temperature of the air (T TEST,in )  80 33 according to the relation T m ¼ T TEST,i À 10 C. during the discharge mode, the temperature of the PCM when cold (T PCM,c ) must always be above the temperature of the cold pressurised air entering the TEST. In the model, T PCM,c is set as 50 C, 70 C and 50 C for TEST 1 , TEST 2 and TEST 3 respectively. Table 2 shows the thermo physical properties and cost of the PCMs. The mass (m) of the PCM in each TEST is estimated as follows: 51 where subscript m stands melting, sol for solid and liq for liquid. The magnitude of the recovered energy depends on the roundtrip efficiency (RTE) of the charge/discharge cycles. The outlet temperature of the air from the TEST (T TEST,O ) is calculated thus: The temperature of the air leaving HX4, HX5 and HX6 and entering the AE1, AE2 and AE3 respectively (TIT) in each stage is obtained by energy balance as follows: The outlet temperature of the DFE exhaust gas (T ex,O ) is set to the stack temperature of the selected dual fuel engine which is 120 C to prevent the corrosive effects of condensation in the exhaust piping. 52 3.5.3 Air store volume. The output power of the A-CAES system depends on the volume of air compressed during charging. The volume of the storage chamber (V AS ) is estimated as follows: 3.5.4 Air expander (AE). The outlet temperature of air exiting the AE in each stage T AE,o is calculated as follows: The charge and discharge time ratio the AC to the AE is directly dependent on the discharge to charge mass ratio and its density and can be represented as: 53 The generated power output by the A-CAES (P AE ) and the total power output by the hybrid system (P T ) at a given instant is calculated by the following equations: 3.5.5 Waste heat recovery (WH) for district heating (DH). In the present analysis, the supply of medium low pressure hot water at a network supply temperature of 55 C for domestic hot water or district heating (DH)/low grade heat process requirements 55,56 is considered. Due to the low temperatures, Counter current (Cc) shell-and-tube type heat exchangers with high efficiency and low operating temperature difference is considered. 57 Heat can be recovered from three sources in the system (if available). These are the AE exhaust, the cooling water jacket and the hot exhaust gas. 58,59 In typical engines, 30% of the fuel energy input is contained in the jacket water and is capable of producing 90 to 99 C hot water. 52 Based on the thermodynamic analysis using energy balances, the following equations can be used to estimate the recovered thermal energy from the cooling water jacket. The heat content of the cooling water is estimated based on the assumption that 30% of the fuel input is lost in the cooling water. 52 RH cw ¼ 0.3h HX (P th,sg + P th,df ) Counter current ow conguration is considered for all HXs and the area of all heat exchangers is modelled using LMTD method. 60 The energy balance equations applied in the HX is as follows: where the subscript c and h denote cold and hot respectively. The heat transfer area of the HX is obtained as follows: In the analysis, U value of 120W m À2 K À1 , is assumed. 61 Similar procedure is followed for the calculation of the area of other HXs. 3.6 System performance criteria 3.6.1 Energy efficiency. The system energy efficiency is analysed with two parameters: the system electrical efficiency (SEE) and the total system efficiency (TSE). SEE is the ratio of the useful net energy output from the system to the total energy input. TSE is the ratio of the sum of the net useful electric output and net useful recovered heat from system, to the total fuel energy input. 40,53,62 TSE ¼ (P T + RH cw )t dch /(P WT + E fuel )t ch (40) where E fuel is the total fuel input consisting of biomass energy input and diesel fuel energy inputs to the DFE system The term useful output (for both electricity and heat) covers all the electricity or heat produced, no matter what the demand in the study area is. Table 3 shows the values for the input parameters used in the A-CAES model.

Equipment cost estimation
where is cf $À£ conversion factor for dollar to pounds; BY is base year and CY is currency year. The capital cost of the TEST including the AS is made up of the capital cost of the steel tanks, the PCM lling it and insulation as follows: 67 In the model, an insulation thickness of 0.038 m and density of steel of 7900 kg m À3 is used. 67 The cost of stainless steel (c steel ) is obtained as 3.18 £ per kg from. 68 The capital cost of the A-CAES assembly is estimated with the following equation: The specic capital cost of biomass gasier coupled to an internal combustion engine electric generator set (BMG + EGS) for power generation and CHP was quoted by IEA 69 to lie in the range of 3000 to 4000 ($ per kW e ) in year 2007. Because of technology learning, the lower value of 3000 $ per kW e is used in the analysis and increased by 5% for contingencies and BOP.
The purchase equipment cost of the hot air dryer (HAD), cyclone and scrubber was obtained by scaling using the "sixtenths rule". 65 In this approach, the costs (c G ) of a given plant size ( _ m G ) is obtained from a known or reference plant size (ṁr ef ) and cost (c ref ) using the following equation: 65 where 4 is a scaling factor. Table 4 shows the reference cost factors used in the computation. An average 4 value of 0.6 has been used in the analysis. 65 Finally, cost of the hot water tank (c HWT ) vary widely depending on its volume, quality and its capabilities. In the model, a modest low cost value of 925 (£ per m 3 ) has been used. 70 Cost escalation factor (e) also called ination rate, is used to account for the increase in cost of components and services over time. Naturally, these cost parameters are all subjected to strong uncertainties; therefore, in the present study all item costs were escalated from their base year to 2017 using only one ination rate e.

Economic analysis of the integrated system
To evaluate the proposed integrated system economic performance, the cost of electricity (COE), total life cycle cost (TLCC) and net present value (NPV) were computed. The COE is the minimum price at which energy must be sold for an energy project to break even.  First an estimate of the TLCC for the system is carried out as follows: 71 .
The total investment cost of the integrate system is the sum of the total capital cost (c T x ) of the various components that make up the system as follows: The remaining direct costs (c o ) and indirect costs (c i x ) are related with equipment installation, piping, instrumentation, controls, electrical equipment and materials, land, civil structural and architectural work, and service facilities. These costs cannot be estimated directly. In the model, they are calculated as 5% of the purchased equipment costs i.e. (c 0 + c i x ¼ 0.5 x ) The annual operation and maintenance (O&M) cost can be expressed as a percentage of the capital cost. 4,11 In the model, the O&M cost is estimated as 5% of the TIC. 4 The annual cost of the input energy (c En,i ) i.e. the cost of biomass (c w,BM ), diesel (c df ) and excess wind energy (P WT ) cost (c r,e ) used by the integrated system is estimated as follows: The cost of wet biomass feed (c w,BM ) at any moisture content (MC) for industry is estimated as follows: Data for tting c w,BM above has been taken from the report AEBIOM-wood fuels handbook. 72 Finally the COE is calculated as follows: 71 .
The ultimate capital recovery factor (UCRF) is calculated using.
The real discount (d r ) rate is calculated from the nominal discount rate (r) using the following formula: 71 .
In the constant dollar analysis, it is necessary to convert cash ows for depreciation from current to constant dollars using the ination rate (e). The following expression shows the conversion from current depreciation (D) to constant depreciation (D r ): The present value (PV) of the O&M, energy cost (c En ) and depreciation (c D ) is obtained by multiplying the appropriate energy cost, with the present value factor as follows.
The net present value (NPV) which represents the present value of all the expenses and revenues streams over the life of the system is evaluated with the following: 11,71 where R is the revenue generated by selling electricity and heat produced by the system. R is estimated as follows: In the model, the excess wind electricity price (c r,e ) while charging is assumed to be 42 £ MW per h, the 2017 average offpeak spot price. 73 For discharging (selling), the electricity consumer tariff/price (c ect ) is assumed to be 14.05 p per kW h, the UK national average price per kW h for electricity. 74 The electricity buy back tariff (c ebt ) and the electricity sell back tariff (c est ) is assumed as 43 £ MW per h, the average 2017 wholesale electricity price in the UK. 75 Table 5 shows the representative values for the main cost parameters of the integrated system. All results have been reported in the 2017 constant dollar year. The methodology described in the preceding section is applied to  34 71 analyse the performance of the system. Fig. 4 shows a physical set up of the algorithm for the implementation of the model in Matlab.  25.76%, and the A-CAES system about 74.24%. It is interesting to observe that about 86% of the excess wind electricity used during the charging mode is produced in the A-CAES during the discharging mode. This characteristic demonstrates the hybrid nature of the A-CAES + BMGES system and it is so because of the added syngas and exhaust heat energy which is used to further heat up the air to a higher temperature aer compression in the compressor and during expansion mode respectively. The total heat recovered from the system for district heating is about 586.37 MW h th per year representing only about 0.16% of the total input fuel energy. This recovered heat, is 100% from the DFE cooling system jacket water.

Results and discussion
Although about 60 to 70% of the energy input of the fuel in the DFE is contained in the thermal energy contained in the exhaust gas and cooling systems, 52 the total percentage of the recovered heat is low because part of the heat content of the exhaust gas is used to heat up the air aer the TEST before expansion in the AE stages. The system electrical efficiency is about 28.97%; the total system efficiency for the production of heat and power is about 36.81%. Table 7 shows the technical performance of the overall integrated A-CAES + BMGES system. The thermal storage capacity of the thermal energy storage tanks, TES 1 , TES 2 and TES 3 is estimated as 0.19, 0.16 and 0.36 MW respectively. As can be observed, the thermal capacity of TES 3 is about double that of TES 1 and TES 2 . This can be explained by higher temperature difference in TEST 3 made possible by highest and lowest temperature of HTF entering the TEST during the charging and discharging modes respectively when compared with TEST 1 and TEST 2 .
The analysed A-CAES + BMGES system is not a single ES technology but a hybrid ES system. Thus comparing its efficiency with an A-CAES system can be misleading. However, the efficiency values computed (both SEE and TSE) assume fairly moderate values when compared with data found in the technical literature for both A-CAES and similar hybrid systems. As a matter of fact, Garrison and Webber 27 reported an overall efficiency of 46% for a coupled solar-CAES system. Zhao et al. 78 reported RTE ranging from 70-74 for a hybrid A-CAES and ywheel energy storage system (FESS). In another study, electrical efficiencies ranging from 24.79-62.09 was reported by Karellas and Tzouganatos 79 for an A-CAES system of different conguration (i.e. different stages with and without air preheating). Liu and Wang 24 presented RTEs that vary between 53.9% and 66.98%.
Facci et al. 26 estimated system electrical efficiencies of 30.2% for tri-generative compressed air storage (T-CAES) compatible for a small to medium size civil application. It was reported that the efficiency could vary between 20 to 40% depending on the operating storage pressure (10-100 bar). 80 In another study Kim and Favrat 81 presented efficiencies that vary from 23.6% to 71.6%. The wide variation in the reported efficiencies of CAES emanates from the fact that different congurations, operating parameters and conditions are applied by the various authors in their analysis. Such parameters include (a) the average storage pressure, (b) the compressors and turbines efficiencies, (c) the number of stages of AEs and ACs. Table 8 shows the important  It can be observed that the specic cost of TES 1 and TES 2 (£ per kW h) are more than that of TES 3 by over 40%. This is because, more heat is stored in TEST 3 owing to higher HTF temperature entering it when compared with TEST1 and TEST 2 . In addition, the cold PCM temperature in TEST 2 is higher than in TEST 3 since the temperature of the compressed air from the air store entering TEST 3 is at ambient. In TEST 1 and TEST 2 , the cold temperature of the PCMs is set by the minimum TET from the AEs which is higher than ambient temperature. It is fascinating to assess the shares of cost of the different components of the A-CAES + BMGES system in the TIC as in Table 8. The CAES has the highest share of all the components of the system. It accounted for about 59.8% of the TIC consisting of 15.4% for AC, 13.5% for AE, 5.1% for TES and 25.8% for AS. The biomass gasier with cleaning system and the DFE has the second highest share of 35% with the syngas cleaning system contributing 6.3%. The recuperator accounted for 4.45%. The NPV for the system is found to be negative with a value of £2 144 062 indicating the non-protability of the system in the selected location.

Sensitivity analyses
4.1.1 Impact of cost factors. The performance of the integrated system depends on a number of technical and cost factors. Cost factors include the total investment cost (TIC), the discount rate (d), the fuel price, the O&M cost (c O&M ), ination rate (e), and cost of excess electricity (sold and bought). The technical factors are the syngas fraction (sf), capacity factor (CF), round trip efficiency (RTE) of the A-CAES TES system, the syngas temperature (SGT) and the exhaust gas temperature (EGT) of the DFE. In the sensitivity analyses, the effect of changes in the technical and economic parameters were studied and analysed so as to evaluate their effects on the overall specications of the system. Specically, in the cost sensitivity, the base line value of each cost parameter input into the COE calculation is adjusted by multiplying each of the base line cost parameters by a base line cost component multiplier (BCCM) that ranged from a smallest value of 0.3 to highest value of 1.5 in accordance with IEA methodology for energy plants. 4 Then a sensitivity plot (SP) of COE for the variation of the factors is illustrated. The greater the slope of the COE is with respect to the relative change in the cost factor, the bigger the inuence on the COE and vice versa. Fig. 5 shows the results of this sensitivity analysis for the COE. As it can be observed, the TIC has the steepest slope of the factors and thus has the strongest inuence on the COE. This is because TIC remain the single biggest element of the system's cost. If the TIC increases from its base line value by a BCCM of 1.3 (30% increase) while all other variables are held constant, the COE increases to 0.23 £ per kW h. This represents a percentage increase in the COE of about 19.5%. However, at a BCCM of 0.7 (30% reduction), the COE reduces by about 19.5% to 15.25p per kW h which is outside the range of current UK national average price for electricity which is 14.05p per kW h. 74 This fairly high percentage change in COE with variation in the TIC highlights the signicant sensitivity of the COE to change in the TIC.
The discount rate (d) has the second highest impact on COE. The COE is found to increase with increasing discount rate. This is because the discount rate inuences the costs of capital. The higher the discount rate, the higher the cost of capital and vice versa. The COE for the discount rate varied from a low value of about 0.14 £ per kW h at a low BCCM of 0.3 (d ¼ 3%) to a high value of about 0.23 £ per kW h at a BCCM of 1.5 (d ¼ 15%), representing about 21.7% increase in COE from the base value.  Fig. 9 shows the EGT sensitivity to the COE of the A-CAES + BMGES system. The most immediate and direct impact of EGT on the system is the variation of the TIT and hence TSE. In dual fuel mode, the EGT is usually higher than the diesel only mode at the same load. 83 This may be due to slow combustion ame speed resulting from the incomplete combustion features of the syngas which leads to a fairly higher exergy loss through exhaust gas. The typical EGT of dual fuel compression ignition (CI) engine operating on diesel fuel and rened rice bran oil blends, having a compression ratio of 18.4 : 1 can vary from a low value of 341 C at 63% load and 50% sf to a high value of 524 C at 98% load and 75% sf, depending on the load, DF, compression ratio and ambient conditions. 84 Based on 120 C as the nal stack temperature, varying the starting temperature of the DFE exhaust gas temperature from 361 C (634 K) up to 541 C (814 K) results in about 75% increase in temperature difference available and thus, available energy for recovery. The degree to which the energy obtainable in the higher EGT can be utilised depends on the characteristics of the system. Within the EGT range of studied, it is observed that increase in EGT improved the performance of the system. The TSE increased from 36.2 to 37.3%. The increase in the TSE is due to increase in available heat for recovery at high EGT which leads to increase in TIT as can be seen in Fig. 8. As expected the increase in the overall performance of the system at higher EGT led to improved COE of the system. Fig. 10 shows the sensitivity of COE of the system on the temperature of syngas (SGT) exiting the gasier during the charging mode. As expected, the COE of the system varies inversely with the SGT. As the SGT reduces from the base case value of 973 K to 773 K, the COE increased from 0.1895 to 0.1903 £ per kW h. However, as the SGT is increased from the base value to 1073 K, the COE reduced to 0.1891 £ per kW h. The trend can be explained by the increase in the TSE with SGT. With increase in SGT from the base value to 1073 K, the higher TSE reveals the benet of reheating on efficiency. TES efficiency rises with air re-heating since the increased turbine inlet temperature (TIT) in the AE stages leads to increased power outputs and hence TSE as can be seen in Fig. 9. Fig. 11 shows the percentage variation of the NPV of the system to the variation of the sensitivity factors examined above from their minimum and maximum values. It can be seen clearly that the most signicant factors swaying the NPV of the A-CAES + BMGES are TIC, cost of O&M, excess wind electricity cost, electricity tariff, and cost of diesel fuel

Limitations of the study
The result of the study should not be considered as representative of the UK or the world in any way. Electricity and heat Fig. 9 Impact of EGT on COE.   Paper demand including energy prices depends on a range of different supply and demand situations, including the geopolitical location, local climate, the national energy mix, network costs and taxation. Hence, the case study using Hull should only be considered as an example of the applicability of the methodology. The methodology can be scaled up to any capacity and used for any location or country while factoring in labour costs, fuels and electricity prices. Moreover, it should be noted the details of the hybrid systems cost and economics are valid elsewhere. The primary differences between the systems would be ensuring they meet the standards of the local regulators in the region of interest.

Conclusions and recommendations
A wind powered A-CAES system and a downdra biomass gasier with dryer coupled to a dual fuel engine are integrated in this paper to produce heating and power simultaneously for supplying some of the energy needs of 1600 households in Hull humber region.
Techno-economic analysis is carried out to investigate the performance of the system. The technical assessment focused mainly on electrical power generation, fuel consumption and efficiencies. The economic analysis is carried out based on total life cycle cost (TLCC), cost of electricity (COE), and net present value (NPV). Moreover a sensitivity analysis is carried out to appraise the effects of important technical and cost parameters on system performance. The main conclusions are summarized as follows: The results obtained indicates that the proposed system is technically possible for commercial distributed power and heat generation. It uses renewable energy and can simultaneously produce electrical and heating energy.
Within the baseline design conditions, the system has the capability to produce a total electric power output of 1.16 MW (869.51 kW from A-CAES and 0.3 from DFE), 0.315 kW heat for sale and 318.03 kg of dry wood chips per hour with electric efficiency of about 29% and total system efficiency for the simultaneous production of heat and electricity of 36.8%. To achieve this, the system requires wind turbine excess electricity of 1 MW, a biomass gasier with capacity to process 318.03 kg of wood per hour and a 0.3 MWe dual fuel engine with syngas and diesel consumption of 23.74 and 881.77 kg h À1 respectively.
Results of economic analysis shows that the COE of the system (0.19 £ kW h) is outside the range of the average cost of electricity for a medium user home in the UK (2500 to 5000 kW h a year), including taxes which is 15.2p per kW h.
Results of the sensitivity analysis of the technical parameters showed that increasing the EGT, RTE of A-CAES TEST and SGT increases the TSE of the system. In addition, it is concluded that since the A-CAES is interconnected with BMG and DFE, any variation in the parameter of BMG or DFE, will cause a variation in the performance of the A-CAES + BMGES system.
The economic results of the system and the sensitivity (Fig. 5-10) provide investors with clear guides to the COE of the hybrid systems for distributed generation. The negative NPV of the system even with modest TIC reduction coupled with high COE cost as compared to global average COE of about 0.1 $ per kW h shows that the integrated system is not economically viable for distributed electricity generation in the selected location unless about 70% the TIC of the system is met by a government subsidy.
The economic position of the system may be further enhanced in the future by: Reducing CAES and BMG plant cost, which can occur through efficiency improvements as well as commercialisation.
Increasing sell back to the grid tarrifs for power exported by private commercial generators.
Granting high subsidy to investment in hybrid biomass and wind generation systems.
Possible reductions in O&M, fuel, and other BOP costs have a moderately small effect, but of course should not be ignored as part of overall cost reductions.
It should be taken note of, that these results are independent of renewable energy certicates, feed-in tariffs or other government subsidies. The expected reduction in A-CAES and BMG costs due to commercialisation, standardisation, lower cost of capital and technology learning in the future is likely to improve the economy of the overall system.

Conflicts of interest
The authors declare no conict of interest.