太阳能小屋的设计(最终).doc

2016重庆邮电大学大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 重庆邮电大学 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： 日期：2015年8月24日赛区评阅编号（由赛区组委会评阅前进行编号）：2016重庆邮电大学大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：摘要本文用EXCEL软件对给出的山西大同典型气象年逐时参数进行全面性分析，进而计算出各个类型的光伏太阳能电池板的各项参数，采用模糊综合评价的模型在光电池的功率，转换效率，工作时长以及价格进行比较，选择出最佳的光电池问题一：以各光伏太阳能电池板的额定功率为阀值，筛选出以额定功率工作的时长和低于额定功率状态时所做的功，通过模糊综合评价的模型对各电池板的性能进行综合性评价，再计算出各光电池一年内所获得利润大小，最后选出合适的电池板为B2和A3。根据小屋各个面的面积确定出电池板的数量，进而选出合适的逆变器。在35年使用寿命内，经济效益约为36.5271元，投资回收年限为28年。问题二：在第一问的基础上，考虑到地理纬度，电池板倾斜角度等因素的影响，我们对太阳方位角、太阳高度角、太阳赤纬角、太阳时角进行了量化处理，通过月总辐射量在全年范围内求和，利用matlab工具采取计算机循环寻优算法，计算出电池板的最佳倾角为44，沿用解决问题一的思路对逆变器进行了选择。在35年使用寿命内，经济效益约为45.1639万元，投资回收年限为22.17年问题三：基于之前的计算结果和结论，并对小屋的建筑要求进行了线性规划，用LINGO软件进行处理，找到了小屋面积，朝向及其屋顶倾角的最合理的设计方法，选出了相应的逆变器。在35年使用寿命内，经济效益约为49.0331万元，投资回收年限为21.12年。 关键词：太阳能光伏电池板；模糊综合评价；投资回收年限；最佳倾角；计算机循环寻优；一、 问题重述在设计太阳能小屋时，需在建筑物外表面（屋顶及外墙）铺设光伏电池，光伏电池组件所产生的直流电需要经过逆变器转换成220V交流电才能供家庭使用，并将剩余电量输入电网。不同种类的光伏电池每峰瓦的价格差别很大，且每峰瓦的实际发电效率或发电量还受诸多因素的影响，如太阳辐射强度、光线入射角、环境、建筑物所处的地理纬度、地区的气候与气象条件、安装部位及方式（贴附或架空）等。因此，在太阳能小屋的设计中，研究光伏电池在小屋外表面的优化铺设是很重要的问题。参考附件提供的数据，对下列三个问题，我们分别给出小屋外表面光伏电池的铺设方案，使小屋的全年太阳能光伏发电总量尽可能大，而单位发电量的费用尽可能小，并计算出小屋光伏电池35年寿命期内的发电总量、经济效益（当前民用电价按0.5元/kWh计算）及投资的回收年限。问题1：根据山西省大同市的气象数据，仅考虑贴附安装方式，选定光伏电池组件，对小屋（见附件2）的部分外表面进行铺设，并根据电池组件分组数量和容量，选配相应的逆变器的容量和数量。问题2：电池板的朝向与倾角均会影响到光伏电池的工作效率，选择架空方式安装光伏电池，重新考虑问题1。问题3：根据附件7给出的小屋建筑要求，为大同市重新设计一个小屋，要求画出小屋的外形图，并对所设计小屋的外表面优化铺设光伏电池，给出铺设及分组连接方式，选配逆变器，计算相应结果。二、 问题分析对于问题一,我们认为在仅考虑贴附安装方式的情况下，只要在光伏太阳能电池板的收入大于其需要的成本的条件下，在给定的面积中让其利用率达到最大即可保证最大收益，从而建立模型求解。对于问题二，我们认为地理纬度，光伏电池板倾斜角度以及安装方位角等因素均会影响到电池板的工作效率，且以架空方式进行铺设，所以我们对一些主要的影响因素进行比较后建立模型，并结合问题一中得出的结论求出更加合理的铺设方案。对于问题三，我们将采用前两问中的计算结果，本着成本最低，可实际操作，环境友好的原则，用太阳能电池板利用率最高的方式来设计太阳能小屋。综合上述三个问题的理解和分析，本文着重考虑了光伏太阳能电池板现有技术条件下的可用性与实用性，使这种环境友好型的新型技术可以得到推广和优化应用，给我们的生活带来更多便利。三、符号说明符号说明W功率（瓦）l收益（元）t工作时长（小时）电池的工作效率F太阳辐照值（Wm2）s面积（m2）太阳方位角（）h太阳高度角（）太阳赤纬角（）当地纬度（）太阳时角（）光伏电池板倾斜角度（）光伏电池安装方位角（）电池工作因子四、问题假设1.电池的工作因子=70%为常数。2.由于屋顶与水平的夹角cos =0.983，近似于1。近似将屋顶作为水平面处理。3.逆变器在屋里，不占用光伏电池空间。4.在一月、二月、十一月、十二月份屋顶有积雪取0.7，其他月份取0.22。5.光电池不受天气状况影响。五、模型的建立与求解5.1问题一5.1.1模型的准备首先对三种电池进行分析，通过附件1我们得知不同型号的电池对太阳辐照阈值要求不同，同时结合附件4，我们对一年中各个型号电池的工作时长作统计。得到下表：不同型号工作条件和时长对比表型号工作条件东墙一年工作时长（小时）西墙一年工作时长（小时）北墙一年工作时长（小时）南墙一年工作时长（小时）屋顶一年工作时长（小时）A200Wm21168147017121241860B80Wm22133314289931422480C30Wm232413713308837992817（表1）5.1.2模型的建立运用模糊综合评价法，首先确定隶属函数。24种电池的因素评论域为功率、转换效率、工作时间和价格。进行单因素评价，功率的区分为4个模糊集合。W1、W2、W3、W4.将其用柯西型模糊分布表示为：W1=0 x28011+2280(x-280)-2 x280W2=11+(x-210)22210W3=11+(x-140)22140W4=1 x7011+(x-70)2270 x70转换效率也分为4个模糊集合。1、2、3、4，也用柯西模糊分布表示为：1=0 x1611+(x-16)-2 x162=11+(x-12)23=11+(x-8)24=1 x411+(x-4)2 x4以此类推价格Q的4个模糊集合为：Q4=0 x350011+23500(x-3500)-2 x3500Q3=11+(x-2500)222500Q2=11+(x-1500)221500Q1=1 x50011+(x-500)22500 x500工作时间t（为每年每种电池各个面的平均工作时长）的4个模糊集合为：t1=0 x300011+23000(x-3000)-2 x3000t2=11+(x-2400)222400t3=11+(x-1800)221800t4=1 x120011+(x-1200)221200 x1200计算以上公式，建立相关4个评价指标功率、转换效率、工作时间和价格的模糊关系矩阵R。R=r1,1r1,24r4,1r4,24,其中ri,j为第j个电池种类对第i个评价的隶属程度。对于每个R矩阵乘一个14的评语量化集D=(0.4,0.3,0.2,0.1)，即每项评价的得分状况。从而得到综合评判矩阵4个B=DOR=(b1,b2,b3,b4,b5,b24)。将这4个B矩阵先合成再分为24个41的C矩阵，代表了每种电池对于4个指标的综合评价。确定评价因素权向量，对于工作功率，转换效率，工作时间和价格，选择评价因素权向量A=a1,a2,a3,a4=(0.4,0.35,0.15,0.1)。每种电池的综合评价矩阵C与A作乘积，得到每种电池的最终评价结果为矩阵S。WA、WB、WC分别为A、B、C种电池的功率。那么它们在额定功率下工作一年所产生的利润是：lA=WA10-3tA0.5lB=WB10-3tB0.5lC=WC10-3tC0.5以屋顶为例，结合附件3的数据带入计算在辐射度F=1000W/,光谱为AM1.5，温度为25标准情况下。且实际发电量计算时可以不考虑AM值的影响，根据电池表面接收到的太阳总辐射量参数进行计算得出下表：不同型号电池利润表（表5-2）电池型号利润平均单位利润A1199.95156.6220704A2302.25155.9278909A3186145.6949492A4251.1153.316172A5227.85139.3450142A6274.85141.5345471B1328.6200.9601566B2396.8204.7053337B3260.4177.1255061B4297.6182.9268293B5347.2178.9366053B6365.8188.5224949B7310185.8513189C1110.8585.26923077C264.29368.45280837C3110.8570.37219495C499.76564.78246753C5110.8571.98051948C64.43440.29077692C74.43440.05420054C88.86840.61834421C913.30240.72872015C1013.30245.80736251C1155.42547.32164202如果考虑不在额定功率工作的情况，结合表5-3，对于每种电池，FW额s时在额定功率工作。如果FF）假设选择n块i号电池和j号逆变器，且这n块电池有串联a个在并联b串。那么j号逆变器关于容量W、允许输入电压V输入的约束条件为：nWiWjVjminaViVjmaxab=n5.1.3模型的求解求解5.1.2我们可以得到最优为A3其次为B2最优且相差不大。从表5-2得出B2利润要大于A3。所以在条件允许的情况下，优先选择利润最大的。得出关于逆变器和B2，A3的串并联情况为下表所示：逆变器关于两种光伏电池的排布表逆变器种类B2A3B2单位价格A3单位价格个数排列方式个数排列方式SN10-0-SN20-0-SN32并联4并联22501125SN45并联8并联1380862.5SN57并联12并联1457.143850SN612并联20并联1250750SN77最多串联3个12最多串联3个1457.143850SN812最多串联3个20最多串联3个1275765SN925最多串联3个40最多串联3个1400875SN1050最多串联3个80最多串联3个1276797.5SN112可串联4-6个4可串联4-6个22501125SN125可串联4-6个8可串联4-6个1380862.5SN137可串联4-6个12可串联4-6个1471.429858.333SN1412可串联4-6个20可串联4-6个1275765SN1518可串联4-6个30可串联4-6个1222.222733.333SN1625可串联4-6个40可串联4-6个1400875SN1731可串联6-17个50可串联6-17个1411.29875SN1837可串联6-17个60可串联6-17个1478.378911.6667（表3）从表5-3中可以得出不同逆变器和不同电池的最佳组合。1. 屋顶B2在F165Wm2时在额定功率下工作。由附件2筛选得到符合165Wm2F80Wm2的值有444个结果，且F=54482.08Wm2.l1=16.39%195699110-654482.0810-30.50.7=6.0582筛选得到F165Wm2的值有2036个结果。即一年中按照额定功率工作时长t=2036小时。l2=32010-320360.50.7=228.032l=l1+l2=234.09两个顶面一共用29块B2， 成本价为400029+6900+35000=157900l顶=234.0929=6788.61根据题意要求，在电压和功率的制约下，本着价格低廉的方式，两个顶面一共有29块B2，用1个16号逆变器分管较大顶面上的24块4个串联6组并联和1个4号逆变器并联5块电池在小顶面。屋顶排列光电池（图1）2. 北墙附件2中北墙辐射强度符合70Wm2F30Wm2的共有1956个结果，且F=93800Wm2l1=3.281588C1在F70Wm2，有1132个小时。l2=68.46l=l1+l2=71.74由于背面北向阳光，且收益为五面中最少的，所以我们将其忽略不计，即不在北墙做任何修改和变动。3南墙由于B2面积过大，在南墙能铺设的较少，所以这里我们南墙选用光伏电池A3，此种类在F200Wm2正常工作。共有2124个小时。l=159.831经过排列，南墙用八块比较合适，具体如图2所示。l南=1278.648南墙根据表5-3，得出南墙使用1个4号逆变器。成本价为3203.58+6900=32528元。南墙排列光电池（图2）4西墙对于A1来说，有1470个小时工作。lA=110.103相对于B2，165Wm2F80Wm2的值有593个结果。F=69530.53Wm2。lB1=7.73在额定功率下工作1642个小时。lB2=183.9lB=191.63具体排列如图3所示l西=2107.93西墙可以安装11个光电池，用6号逆变器，一起并联。成本价为114000+15000=59000西墙排列光电池（图3）5东墙对于A1来说，有1168个小时。lA=87.89对于B2，165Wm2F80Wm2的值有815个结果，F=91719.61Wm2lB1=10.245在额定功率下工作1318个小时。lB2=147.62lB=157.86l东=1420.74东墙一共有9个光电池，用6号逆变器，都并联。 成本价为94000+15000=51000元东墙排列光电池（图4）所以总成本价为157900+51000+59000+32528=300428元l总=6788.61+1278.648+2107.93+1420.74=11595.928十年收益115959.28，35年收益365271.732元。10-25年每年收益下降为10436.3352元，25以后收益下降为92767.424元共28年可以回本。5.2问题二5.2.1模型的准备太阳方位角和高度角h示意图其中太阳高度角h=90-（-2326sinr365.2422）其中r是具体日期距离春分（3月20日）的天数。各月代表日距离的高度角月份123456789101112月日序151515151515151515151515距离6433533930927824821718615612595高度角29.1660837.3990347.9027960.1092669.089773.1069370.8619862.8316951.2135939.5879930.424326.68562赤纬角8.911842-15.538-24.2626-9.9134514.0096724.2055710.13968-14.4878-24.1344-9.6006714.9574224.16607（表4）太阳方位角cos=sinhsin+sincoshcos太阳赤纬角=23.45sin360284+n365（n为日序）各月代表日日序如下表所示：各月代表日日序月份123456789101112月日序171616151511171615151410年日序174775105135162198228258288318344（表5）太阳直射平面地轴赤道平面太阳时角，0为水平面日出时太阳时角，cs为倾斜面上日落时太阳时角，cr为倾斜面上日出时太阳时角。光伏电池安装方位角为光伏电池与正南方向的夹角。5.2.2模型的建立对于非屋顶的东西两墙，在调整角度和使增大的基本思想是使太阳辐射角度基本接近于法向，但是这样会最大面积的遮挡其它光伏电池的采光，反而使利益最小，所以在东西墙面我们仍然采用贴附安装方式，对于背向阳光的北墙，基本忽略在其以架空方式安装所增加的收益，与第一问中的安排保持一致。仅在南墙屋顶考虑最佳倾斜角度和最佳方位角。查阅资料1得到水平面日出时太阳时角0=arccos(-tgtg)对于冬半年，cs=-cr=0每月的总辐射量HT=A1sin+B1cos+C1，简化得到HT=A12+B12sin+x+C1其中A1=HbDsinsin0cos-0sincos，B1=HbDsin0coscos+0sinsin+13Hd-12G，C1=23Hd+12G, D=sin0coscos+0sinsin，x=tan-1B1A1那么HTmax=A12+B12+C1，此时+x=90x确定便能得到佳角度。类比于冬半年，下半年的变化为将0替换为c=cos-1-tan-tanHT=HbDsinccoscoscos+sinsin+csincos-cossinsin+Hdcos+23+121-cosG由上式可以得到每个月的最佳倾角。并用计算机算出年最佳倾角=44。同时年太阳辐照强度最大，F=1.29441011查阅资料2得到倾角对太阳辐照值的影响为：倾角上的太阳辐照值F倾=F法=coscosh+sincoshcos, cos=sinhsin+sincoshcos。5.2.3模型的求解由于目标地理位置处于北半球，且由附件2易得出南墙的太阳辐射强度最强，所以在这里我们令光伏电池的安装方位角=0。由5.2.2公式可计算得出表6：相关量的值（表6）根据上式得到相关因子约为0.87，并据其修改南墙和屋顶的太阳辐照强度。由5.1.2的公式进行效率的计算。得到12个月份南墙和屋顶的辐照量如下表所示：在最佳角度下各个月份的太阳辐照强度南墙屋顶F200F load sj.txt for i = 1:1:90for j = 1:1:12sum(j,i) = sj(j,1)*sin(i*pi/180)+sj(j,2)*cos(i*pi/180)+sj(j,3);endendfunction d= d(a,b,c)x = a*sin(b)*sin(c);y = cos(b)*cos(c)*sin(a);d = x+y;function text= juli(b,c,d,e,f)u = b;t = c;w = sin(d)*sin(e)*cos(f);x = e*sin(f)*cos(d);r = u/t;s = w - x;text = r*s;endm = 10(5*1.4317)+10(5*1.3832)+10(5*2.2221)+10(5*1.6605)+10(5*1.4717)+10(5*1.447)+10(5*1.4751)+10(5*1.5792)+10(5*1.4366)+10(5*1.5683)+10(5*1.2629)+10(5*1.4403)function t = synt(b,a)m = size(b,1);n = size(a,2);for i = 1:m for j = 1:n t(i,j) = max(min(b(i,:);a(:,j); endendfunction c = ex4(b,a)t = a ;for i =1:1:24 c(1,i) = b*t(:,i);endlingo：max = b*a/cos(x);3a;a15;3b;b15;a*b74;x0.716;b*tan(x)2.6;附录R1，4*24=0 0.7834 0 0 0 0.2866 0 0.7407 0 0 0 0.2866 0 0 0 0 0 0 0 0 0 0 0 00.9438 0.0308 0.8077 0.1045 0.2553 0.0549 0.1219 0.0335 1.0000 0.3182 0.0789 0.0549 0.2079 0.0335 0.0179 0.0335 0.0283 0.0335 0.0098 0.0098 0.0102 0.0106 0.0106 0.01610.0474 0.0081 0.0722 0.0163 0.0248 0.0115 0.0176 0.0086 0.0541 0.0272 0.0141 0.0115 0.0226 0.1489 0.0400 0.1489 0.1007 0.1489 0.0149 0.0149 0.0158 0.0168 0.0168 0.03340.0066 0.0021 0.0082 0.0035 0.0046 0.0028 0.0037 0.0022 0.0071 0.0048 0.0032 0.0028 0.0043 0.1346 1.0000 0.1346 0.2593 0.1346 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000R2，4*24=0 0 0 0 0 0 0 0 0 0 0 0 0 0.9481 0.9481 0.9481 0.9481 0.9481 0.9481 0.9481 0.9481 0.9481 0.9481 0.94810.0044 0.0044 0.0044 0.0044 0.0044 0.0044 0.7406 0.7406 0.7406 0.7406 0.7406 0.7406 0.7406 0.0055 0.0055 0.0055 0.0055 0.0055 0.0055 0.0055 0.0055 0.0055 0.0055 0.00550.0181 0.0181 0.0181 0.0181 0.0181 0.0181 0.0114 0.0114 0.0114 0.0114 0.0114 0.0114 0.0114 0.0015 0.0015 0.0015 0.0015 0.0015 0.0015 0.0015 0.0015 0.0015 0.0015 0.00150.0877 0.0877 0.0877 0.0877 0.0877 0.0877 0.0018 0.0018 0.0018 0.0018 0.0018 0.0018 0.0018 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005R3,4*24=0.4137 0.2906 0.8794 0.2000 0 0 0.0422 0.1320 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00.0409 0.0444 0.0218 0.0471 0.1012 0.0937 0.0534 0.0493 0.0594 0.1131 0.0594 0.0890 0.1006 0.0383 0.0286 0.0304 0.0257 0.0319 0.0141 0.0141 0.0142 0.0142 0.0159 0.01650.0126 0.0132 0.0087 0.0137 0.0201 0.0194 0.0146 0.0140 0.0155 0.0212 0.0155 0.0189 0.0201 0.4950 0.2299 0.2686 0.1765 0.3049 0.0498 0.0498 0.0504 0.0504 0.0626 0.06710.0060 0.0062 0.0046 0.0064 0.0082 0.0080 0.0067 0.0065 0.0069 0.0085 0.0069 0.0079 0.0082 0.1006 0.1752 0.1533 0.2280 0.1389 1.0000 1.0000 1.0000 1.0000 0.9834 0.9321R4,4*24=0.4137 0.2906 0.8794 0.2000 0 0 0.0422 0.1320 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 00.0409 0.0444 0.0218 0.0471 0.1012 0.0937 0.0534 0.0493 0.0594 0.1131 0.0594 0.0890 0.1006 0.0383 0.0286 0.0304 0.0257 0.0319 0.0141 0.0141 0.0142 0.0142 0.0159 0.01650.0126 0.0132 0.0087 0.0137 0.0201 0.0194 0.0146 0.0140 0.0155 0.0212 0.0155 0.0189 0.0201 0.4950 0.2299 0.2686 0.1765 0.3049 0.0498 0.0498 0.0504 0.0504 0.0626 0.06710.0060 0.0062 0.0046 0.0064 0.0082 0.0080 0.0067 0.0065 0.0069 0.0085 0.0069 0.0079 0.0082 0.1006 0.1752 0.1533 0.2280 0.1389 1.0000 1.0000 1.0000 1.0