Abstract:To overcome the shortcomings of the energy consumption prediction models in the application during the design stage, a quick prediction model for energy consumption is proposed based on the decoupling method. Taking typical residential and office buildings in hot summer and cold winter zones as research objects, the influence factors on building energy consumption are classified into intrinsic factors and operational factors on the basis of the heat transfer principle. Then, using the intrinsic factors as the fundamental variables and operational factors as the modified variables, the quick prediction model for the buildings in typical cold and hot zones is proposed based on the decoupling method and the accuracy of the proposed model is verified. The results show that compared to the simulation results of EnergyPlus, the relative error of the prediction model is less than 1.5%; compared with the real operating data of the building, the relative error is 13.14% in 2011 and 8.56% in 2012 due to the fact that the coincidence factor becomes larger than the design value about 16% in 2011 and 13% in 2012. The finding reveals that the proposed model has the advantages of rapid calculation compared with EnergyPlus and Design Builder when predicting building energy consumption in building designs. The energy consumption prediction model is of great practical value in optimal operation and building designs.
Keywords:building energy consumption; energy conservation; load prediction; EnergyPlus
Biography:Xu Jianqun(1959—), male, doctor, professor, qlj1062@163.com.
Nowadays, with the increasing population, rapid urbanization and high thermal comfort requirement, building energy consumption increases sharply and it will catch up with and even exceed the industrial and transportation energy consumption. In the future, the building will be the most energy-consuming[1-3]. To realize an efficient energy management system, it is crucially important to establish a reliable and accurate prediction model for building energy consumption, which causes the construction of prediction models for building energy consumption to become an essential issue[4-5].
It is difficult to precisely predict the building energy consumption, since there are many factors influencing energy consumption, such as the change in weather conditions, structures of the building, geographic location, layout and the operation of internal facilities[6-8]. The prediction accuracy directly affects design, configuration and operation. During the past several decades, there have been various forecasting techniques that are applied to the building energy consumption calculation, such as engineering methods, statistical methods, artificial neural networks, support vector machines and grey model techniques, etc[9], to name only a few. Stephenson et al.[10-12]put forward the reaction coefficient to discuss the heat transfer and heat load. Spitler et al.[13]utilized periodic response factors to reduce the unnecessary computational quantity. The frequency-domain regression method was introduced by Wang et al.[14]to calculate transient heat flow. A report sponsored by the agency of the United States government proposed the key criteria that was used for precise heating and cooling load calculations[15]. Softwares Design Builder and EnergyPlus are usually adopted to calculate energy consumption, the accuracy of which has been validated by the American HVAC Engineering Association[16-17].
These methods mainly concentrate on reducing energy consumption and optimizing the building operation. However, there is little research on how to predict building energy consumption at the stage of building designs. The usage of new materials and new technology can reduce the building’s energy consumption, and thus it is important to comprehensively calculate the building energy consumption in construction and operation[8,18]. On the other hand, sometimes these methods present a hysteretic nature from the perspective of engineering application and require elaborate construction drawings and parameters, whose data is difficult to obtain in the stage of building designs. Nevertheless, designers cannot have the comprehensive understanding, for instance, when the government wants to convert a coal-fired power plant into a CCHP (combined cooling, heating with power) power plant. They should find some consumers to utilize the cooling and heating supply produced by the plant. Under this circumstance, designers hope to know the dynamic characteristics of building load and afterwards configure different types of buildings in a region, so that they can make the peak-valley difference as small as possible. Moreover, with the development of the economy, occupancy density, coincidence factor and schedule may change, which may be different from the design conditions[19-20]. Therefore, the influential factors of building should be grasped and the prediction model that satisfies accuracy requirements should be established.
This work develops an energy consumption prediction model for typical buildings based on the heat transfer principle. This work also analyzes and distinguishes the influential factors of building energy consumption to simplify the prediction model. Finally, the accuracy of the proposed model is verified.
1.1 Heat transfer model
The building energy system is fairly complex; therefore, there are four assumptions when predicting the energy consumption. 1) The building consists of four parts: vertical walls, windows, ceilings, and floor. The heat storage of the building can be ignored when calculating building energy consumption[16]. 2) Only the internal heat generation of the building which includes electrical equipment (including computers, lighting, etc.) and body heat dissipation is taken into account. 3) The windows and doors are shut when air conditioning is running. 4) Indoor heat transfers promptly so that temperature distribution is uniform.
There are many factors influencing the building energy consumption as shown in Fig.1, including solar radiation, air circulation, internal heat generation and so on. As can be seen in Fig.1, conductive heat transfer via the walls and ceiling, convective heat transfer with air circulation, radiation through the windows, and solar radiation absorption are the main heat transfer mechanisms in buildings. In this section, we concisely discuss the building model. The cooling and heating load of the building can be demonstrated in the following equation:
(1)
Fig.1 Heat transfer model of a building
whereqcandqhare the cooling load and the heating load;q1,q2andq3are the heat flux between indoors and outdoors, which is triggered by temperature difference, the solar radiation and the air circulation, respectively; and the internal heat generation of the building is denoted byq4. The total building energy consumption can be obtained from the time integration ofqcandqh.
According to the Chinese National Standard (GB 50189—2005)[19]and the basic principle of heat transfer[20],q1is described by the following general equation:
q1=
(2)
whereTo,TiandTgare the ambient temperature, the indoor temperature and the ground temperature;Awa,AwiandAaare the areas of wall, window and floor, respectively;Awi/(Awa+Awi) is defined as the window-wall ratior;RwaandRwiare the thermal resistances of the wall and the window between indoors and outdoors, respectively; andkgis the heat transfer coefficient. The ambient temperature and the ground temperature are obtained by consulting meteorological data.
q2is calculated by the following equation:
q2=
(3)
whereq2consists of the direct radiative heat flux density of architectural enclosure and the indoor radiative heat flux density;εstandεabare the correction coefficients of the enclosure structure and solar radiation absorption;kstandkanare the heat transfer coefficients of the enclosure structure and solar radiation angle. The parametersεst,εab,kstandkanare related to building orientation denoted byθ[19].
q3is calculated as
q3=
(4)
whereCpis the average specific heat capacity at constant pressure;NandVdenote the rate and volume of air exchange per unit time, respectively.Nis closely affected by indoor heat distribution, proportional toq4below.
The ambient temperature in Eqs.(2) to (4) is presented in Fig.2, which is the hourly air temperature of Nanjing, a typical city in the hot summer and cold winter zone in China. The date of Fig.2 can be obtained from Autodesk Ecotect Analysis.
q4is calculated asgreat influence onqcandqh. In a typical building, the internal heat generation is caused by electrical equipment and body heat dissipation.eelandeboare the electrical equipment density and body heat dissipation,ebo=50.peis the occupancy density.
q4=(eel+ebope)
(5)
wherehis the height of a single-story and its value is 3.6 m. The internal heat generation of the building has a
Fig.2 Ambient temperature of Nanjing
1.2 Calculation process
In the proposed model, the influential factors of the cooling and heating load can be divided into two categories: intrinsic and operational parameters. Their concrete explanations are as follows:
1) Intrinsic parameters contain orientation, thermal resistance of walls and windows, window-wall ratio, and so on. When the building has been built, these parameters will not transform basically.
2) Operational parameters include schedules, the staff and their behavior, the comfortable indoor temperature, electrical equipment, and so on. In an engineering application, these parameters will vary with users’ requirements frequently.
The proposed model is a multi-input and multi-output system. In order to simplify the calculation, the multi-variable system is divided into several subsystems, and each subsystem has a function matrix, as shown in Fig.3. The cooling and heating loads are the coupling consequence of multiple factors. The research strategy is demonstrated in Fig.4.
Fig.3 Division of multivariable process
Fig.4 Research strategy of the predictive model
In the model, the detailed procedure can be conducted as follows:
Step1 Operational parameters are considered to be the most important factors in this prediction model, which directly determine the application of the building.
Step2 When the application of the building has been determined, intrinsic parameters can be used to revise the load of building. Intrinsic parameters will produce an opposite effect on the load of building when the internal heat generation is greater than a certain value. Therefore, the coupling relationship between intrinsic and operational parameters is taken into account.
Step3 The temperature difference between indoors and outdoors and the radiative heat flux absorbed by the building change with time. Therefore, a temperature difference coefficient matrix (TDCM) and radiation coefficient matrix (RCM) are employed when building the mathematical models.
2.1 Influence analysis of operational attribute parameters
When studying the impact of operational parameters, a group of intrinsic parameters are assumed unchanged: site=Nanjing,θ=0,Rwa=2 (m2·K)/W,Rwi=0.5 (m2·K)/W,r=0.5.
It is assumed that the indoor comfortable temperature is 26 ℃ in summer and 20 ℃ in winter. Figs.5 to 7 show the change of cooling load in the typical day (July 15th) and cooling months, and the relationship between cooling load andq4in cooling months.
Fig.5 displays that the cooling load in the typical day in summer has a certain proportional relationship withq4. The coefficient of cooling loadacand the coefficient of heating loadahof each duration whenq4=25 W/m2are shown in Tab.1. Figs.6 and 7 indicate that the cooling load of different months or typical days has a similar proportional relationship withq4. Based on the above-mentioned analysis, we introduce baseline values, proportionalcoefficients, and correction coefficients to express forming energy consumption. Baseline values can be determined byq4using the following equation:
Fig.5 Influence of time andq4on cooling density
Fig.6 Influence of month andq4on cooling density
Fig.7 Influence ofq4on cooling density in typical day
Bc=f1(q4),Bh=f2(Bc)
(6)
whereBcandBhare the average value of the daily cooling and heating load from June to August, and the calculation period is from December to February.BcandBhcan be solved by the polynomial function. The coefficients of the polynomial are listed in Tab.2, whererepresents the coefficient of one degree term ofBcand the rest can be deduced by analogy.
Tab.1 Cooling (heating) load ratio, radiation coefficient and temperature difference coefficient
Duration/hθcθhtcthacah0⁃1000.01940.04560.03530.04521⁃2000.01610.04610.03420.04602⁃3000.01380.04650.03320.04663⁃4000.01170.04690.03240.04714⁃5000.01000.04720.03170.04755⁃60.046200.01040.04710.03120.04776⁃70.04750.09300.01390.04650.03210.04747⁃80.05140.09070.02310.04490.03480.04628⁃90.05650.08650.03580.04270.03820.04419⁃100.06190.08200.04810.04050.04190.041810⁃110.06680.07780.05900.03860.04520.039611⁃120.07060.07430.06720.03720.04780.037912⁃130.07330.07180.07330.03610.04960.036613⁃140.07590.06980.07750.03540.05130.035614⁃150.07810.06870.07880.03530.05290.035015⁃160.07900.06880.07620.03560.05340.035116⁃170.07840.07010.07120.03650.05300.035717⁃180.07620.07200.06480.03760.05160.036718⁃190.07180.07470.05570.03920.04860.038019⁃200.066300.04700.04070.04480.039720⁃21000.04030.04190.04210.041021⁃22000.03400.04300.04000.042222⁃23000.02870.04390.03820.043323⁃24000.02380.04480.03670.0443
Tab.2 Polynomial function coefficient of baseline value
CoefficientB1cB0cB6hB5hB4hB3hB2hB1hB0hValue1.1855.814×101-1.511×10-111.822×10-8-8.758×10-62.134×10-3-2.761×10-11.837×101-4.390×102
2.2 Influence analysis of intrinsic parameters
Intrinsic parameters directly affectq1,q2andq3. The results obtained from partial derivatives ofqcandqh, respectively, with respect toθ,r,Rwa, andRwi, can be expressed in the following terms:
(7)
The polynomial coefficients of the above-mentioned equations are presented in Tab.3, whererepresents the quadratic coefficient ofGc(θ);
represents the quadratic coefficient ofGc(Rwa); other terms can be obtained similarly.
2.3 Influence analysis of coupling factors
The coupling relationship between intrinsic and operational parameters includes three aspects: coupling terms ofq4, TDCM, and RCM.
1) Coupling terms ofq4. Whenq4is very small, the cooling and heating load of the building will become smaller when the wall or window resistance becomes larger.
Tab.3 Polynomial function coefficients of self-attribute parameters
CoefficientValueCoefficientValueCoefficientValueCoefficientValueθ4c 4.563×10-9θ4h 6.864×10-10R3c,wa 1.480×10-3R3h,wa-4.983×10-3θ3c-1.062×10-6θ3h-2.034×10-7R2c,wa-6.902×10-3R2h,wa 4.692×10-2θ2c 7.369×10-5θ2h 2.240×10-5R1c,wa 8.912×10-4R1h,wa-1.653×10-1θ1c-6.544×10-4θ1h-5.529×10-4R0c,wa 1.270×10-2R0h,wa 1.771×10-1θ0c 3.636×10-5θ0h 7.158×10-5R2c,wi-9.988×10-2R2h,wi 4.749×10-1r1c 4.676×10-1r1h-7.373×10-2R1c,wi 1.696×10-1R1h,wi-7.405×10-1r0c-2.059×10-1r0h 3.266×10-2R0c,wi-5.303×10-2R0h,wi 2.249×10-1
However, whenq4is greater than a certain value, an opposite effect may occur.
2) The calculation of TDCM is carried out on the basis of the temperature difference between indoors and outdoors of a building in Nanjing.
3) In order to obtain RCM, software Autodesk Ecotect Analysis is introduced to simulate sunshine, shadow, and lighting of a building in Nanjing, and the results are demonstrated in Figs.8(a) to (d). Figs.8(a) and (b) display the light and shadow of the identical building (θis 0°) at 8:00 and 17:00 on July 15th. Fig.8(c) indicates the light and shadow of a building (θis 40° south by east, other parameters are the same as those shown in Figs.8(a) and (b)) at 8:00. Fig.8(d) represents natural lighting and lamp illumination inside the building.θcandθhare the radiation coefficients.tcandthare the temperature difference coefficients. The values of these coefficients at different moments are given in Tab.1.
(a)(b)
(c)(d)
Fig.8 Analysis on RCM. (a) 8:00 on July 15th (θ=0); (b) 5:00 on July 15th (θ=0°); (c) 8:00 on July 15th (θ=40°); (d) Indoor natural lighting
2.4 Mathematical expression of the prediction model
Based on the above analysis, the following equation is used to calculate the average cooling and heating load of the prediction model:
(8)
whereand
are the sum of the factor of cooling and heating capacity during the working period, respectively;Sc(θ) andSh(θ) are the sum of the radiation coefficient in cooling and heating time, respectively;Sc(t) andSh(t) are the sum of the temperature coefficient in cooling and heating time, respectively.
2.5 Accuracy verification of the model by software
Software EnergyPlus can be used to simulate and analyze building energy consumption comprehensively[21], which is chosen to verify the accuracy of the prediction model proposed in this paper. Two types of typical buildings (Buildings No.1 and No.2) are established and simulated. Figs.9(a), (b), (c) and (d) show the overall appearance and floor plan of the two buildings in turn.
(a)(b)
(c)(d)
Fig.9 Overall appearances and floor plans of the two buildings. (a) Appearance of Building No.1; (b) Appearance of Building No.2; (c) Floor plan of Building No.1; (d) Floor plan of Building No.2
The main parameters are given in Tab.4. In order to investigate the performance of the prediction model, EnergyPlus is also used for performance comparisons and the results are shown in Tab.5. It can be seen from the above investigation that the prediction model generally achieves similar accuracy compared with EnergyPlus and the relative error is permissible.
Tab.4 Key parameters of the two buildings
ParameterBuildingNo.1BuildingNo.2θ/°32.20Aa/m2403067180.11Rwa/(m2·K·W-1)1.852.52Rwi/(m2·K·W-1)0.4350.427r0.540.57N/s-14.97WorkingperiodAllday6:00⁃9:00,11:00⁃13:00,17:00⁃19:00pe/(person·m-2)0.080.26eel/(W·m-2)6789q4/(W·m-2)5892
Tab.5 Comparison of calculation results
ParameterBuildingNo.1BuildingNo.2qcqc,a/(W·m-2)140.0176.1qc,b/(W·m-2)141.3177.0Relativeerror/%0.920.46qhqh,a/(W·m-2)65.691.5qh,b/(W·m-2)67.387.0Relativeerror/%2.52-5.17Notes:qx,aaretheresultscalculatedbysoftwareEnergy⁃Plus;qx,baretheresultscalculatedbythepredictivemodelinthispaper.
2.6 Accuracy verification of the model by literature
In Ref.[22], the energy consumption of a teaching building in Donghua University was simulated by the software Design Builder and the actual building energy consumptions in the year of 2011 and 2012 were measured. The building is located in Shanghai, which also belongs to a hot summer and cold winter zone and the working time is from 7:00 to 18:00, May 20 to October 1 in summer and November 20 to March 20 in winter. The parameters of the building are listed in Tab.6.
Tab.6 Parameters of the teaching building
Parameter/constantValueRwa/(m2·K·W-1)0.6097Rwi/(m2·K·W-1)0.1719r0.5134Aa/m215989.06θ/(°)0N/s-13Coincidencefactor/%45pe/(person·m-2)0.08eel/(W·m-2)35q4/(W·m-2)35.4
According to the values listed in Tab.6, the building load is calculated by the prediction model proposed in this paper. Comparisons of the results are detailed in Tab.7.
It is surprising to find that the relative error between the data in Tab.7 and the results calculated by the prediction model is merely 4.38%, which indicates that the resultsof the prediction model are close to the results of Design Builder. Nevertheless, there is a significant relative error when compared to the actual operating data. Ref.[22] showed that the coincidence factor becomes larger than the design value with the development of Donghua University, about 16% in 2011 and 13% in 2012, which results in a higher energy consumption.
Tab.7 Comparisons of calculation results of the teaching building load
ParameterValueCoolingdemandcalculatedbypredictionmodel/(MW·h)979.68Heatingdemandcalculatedbypredictionmodel/(MW·h)585.73Totaldemandcalculatedbypredictionmodel/(MW·h)1565.41TotaldemandcalculatedbyRef.[22]/(MW·h)1496.67Actualconsumptionof2011/(MW·h)1802.22Actualconsumptionof2012/(MW·h)1711.91RelativeerrorbetweenpredictionmodelandRef.[22]/%4.38Relativeerrorbetweenpredictionmodelandactualoperatingdatain2011/%13.14Relativeerrorbetweenpredictionmodelandactualoperatingdatain2012/%8.56
In this paper, the prediction model of typical buildings is proposed based on the heat transfer principle. First, the influential factors of the cooling and heating load can be divided into two categories: intrinsic and operational parameters. In order to simplify the calculation, the multi-variable system is divided into several subsystems. Then, the influences of operational and intrinsic parameters and coupling factors on building energy consumption are analyzed, and the mathematical expression of the prediction model is presented. Finally, the high accuracy of the prediction model is verified by comparing the simulation results of EnergyPlus and Design Builder and the measured results of other literature. Furthermore, the results satisfy the requirements of engineering application.
Our findings demonstrate that the proposed model has the advantages of rapid calculation and high accuracy compared with other models when predicting building energy consumption for building designs. The prediction model can be used to guide building energy conservation design and provide the basis for the optimization of building operation, which has a high practical value in engineering applications.
[1]Lin B Q, Liu H X. China’s building energy efficiency and urbanization [J].EnergyandBuildings, 2015,86: 356-365. DOI:10.1016/j.enbuild.2014.09.069.
[2]Zhang Y, He C Q, Tang B J, et al. China’s energy consumption in the building sector: A life cycle approach [J].EnergyandBuildings, 2015,94: 240-251. DOI:10.1016/j.enbuild.2015.03.011.
[3]Wu D W, Wang R Z. Combined cooling, heating and power: A review [J].ProgressinEnergyandCombustionScience, 2006,32(5): 459-495. DOI:10.1016/j.pecs.2006.02.001.
[4]Fumo N. A review on the basics of building energy estimation [J].RenewableandSustainableEnergyReviews, 2014,31: 53-60. DOI:10.1016/j.rser.2013.11.040.
[5]Li D H W, Yang L, Lam J C. Zero energy buildings and sustainable development implications—a review [J].Energy, 2013,54: 1-10. DOI:10.1016/j.energy.2013.01.070.
[6]Feng X H, Yan D, Hong T Z. Simulation of occupancy in buildings [J].EnergyandBuildings, 2015,87: 348-359. DOI:10.1016/j.enbuild.2014.11.067.
[7]Ahmad A S, Hassan M Y, Abdullah M P, et al. A review on applications of ANN and SVM for building electrical energy consumption forecasting [J].RenewableandSustainableEnergyReviews, 2014,33: 102-109. DOI:10.1016/j.rser.2014.01.069.
[8]Zhao H X, Magoules F. A review on the prediction of building energy consumption [J].RenewableandSustainableEnergyReviews, 2012,16(6): 3586-3592. DOI:10.1016/j.rser.2012.02.049.
[9]Mitalas G P, Stephenson D G. Room thermal response factors [J].ASHRAETransactions, 1967,73(1): 261-266.
[10]Stephenson D G, Mitalas G P. Cooling load calculation by thermal response factor method [J].ASHRAETransactions, 1967,73(1): 152-165.
[11]Stephenson D G, Mitalas G P. Calculation of heat conduction transfer functions for mult-layer slabs [J].ASHRAETransactions, 1971,77(2): 117-126.
[12]Mitalas G P, Kimura K. A calorimeter to determine cooling load caused by lights [J].ASHRAETransactions, 1971,77(2): 65-72.
[13]Spitler J D, Fisher D E. Development of periodic response factors for use with the radiant time series method [J].ASHRAETransactions, 1999,105(1): 491-509.
[14]Wang S W, Chen Y M. Transient heat flow calculation for multilayer constructions using frequency-domain regression method [J].BuildingandEnvironment, 2003,38(1): 45-61. DOI:10.1016/s0360-1323(02)00024-0.
[15]Arlan Burdick IBACOS, Inc. Strategy guideline: Accurate heating and cooling load calculations[R]. Pittsburgh, USA: Arlan Burdick IBACOS, Inc. 2011.
[16]EnergyPlusTM. Energy plus engineering reference—The reference to energy plus calculation, 867-1008[R/OL]. (2007-08-01)[2016-09-03]. http://appsl.eere.energy.gov/buildings/energyplus/pdfs/engineering reference.
[17]ANSI/ASHRAE 2007. Standard method of test for the evaluation of building energy analysis computer programs[R]. Atlanta, GA, USA: America Society of Heating, Refrigerating and Air-Conditioning Engineers, 2007.
[18]Mihalakakou G, Santamouris M, Tsangrassoulis A. On the energy consumption in residential buildings [J].EnergyandBuildings, 2002,34(7): 727-736. DOI:10.1016/s0378-7788(01)00137-2.
[19]Ministry of Construction. GB 50189—2005 Design standard for energy efficiency of public buildings [S]. Beijing: China Building Industry Press, 2005. (in Chinese)
[20]Rohsenow W M, Hartnett J P, Cho Y I.Handbookofheattransfer[M]. 3rd ed. New York: McGraw-Hill Professional, 1998.
[21]Crawley D B, Lawrie L, Winkelman F C, et al. Energy Plus: Creating a new-generation building energy simulation program [J].EnergyandBuildings, 2001,33(4): 319-331. DOI:10.1016/s0378-7788(00)00114-6.
[22]Zhou Jian. Research on building energy consumption simulation and energy-saving reform plan [D]. Shanghai: College of Environmental Science and Engineering of Donghua University, 2014: 19-27. (in Chinese)
References
摘要:为了克服建筑能耗预测模型在建筑设计阶段的应用不足,通过解耦方法建立了建筑能耗快速预测模型.以中国冬冷夏热地区典型建筑为研究对象,依据传热学的基本原理,将影响建筑能耗的因素分为运行属性因素和自有属性因素两类,以运行属性参数为基准,自有属性参数为修正,依据主元解耦方法,建立了典型建筑冷热负荷计算快速预测模型,并验证了该数学模型的准确性.研究结果表明:与软件EnergyPlus相比,该模型计算结果的相对误差在1.5%以下;与实际运行数据相比,该模型计算结果的相对误差在2011年为13.14%,2012年为8.56%,其主要原因是建筑的同时率分别高于设计值16%和13%;与软件EnergyPlus和Design Builder相比,该模型计算快速、简捷且能够用于设计阶段.该建筑能耗快速预测模型对建筑的设计和运行优化具有一定工程价值.
关键词:建筑能耗;能量守恒;负荷预测;EnergyPlus
中图分类号:TK11
JournalofSoutheastUniversity(EnglishEdition) Vol.33,No.3,pp.355⁃361Sept.2017 ISSN1003—7985
DOI:10.3969/j.issn.1003-7985.2017.03.015
Citation:Xu Jianqun, Zhang Fang, Chen Feixiang, et al. Energy consumption prediction model of typical buildings in hot summer and cold winter zone of China[J].Journal of Southeast University (English Edition),2017,33(3):348-354.
DOI:10.3969/j.issn.1003-7985.2017.03.015.