1、太阳能小屋的设计12012高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（
2、包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： 日期： 年 月 日赛区评阅编号（由赛区组委会评阅前进行编号）：2012高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：太阳
3、能小屋的设计摘要太阳能作为一种新兴的可再生能源，在生产和生活中具有广泛的应用，太阳能电池作为一种将太阳辐射能转换为电能的新型电源，具有永久性、清洁性和 灵活性三大优点，但其发电量受诸多因素影响，成本也较高，因此研究光伏电池在房表面的优化铺设是一个很重要的问题。一、 问题的重述在设计太阳能小屋时，需在建筑物外表面（屋顶及外墙）铺设光伏电池，光伏电池组件所产生的直流电需要经过逆变器转换成220V交流电才能供家庭使用，并将剩余电量输入电网。不同种类的光伏电池每峰瓦的价格差别很大，且每峰瓦的实际发电效率或发电量还受诸多因素的影响，如太阳辐射强度、光线入射角、环境、建筑物所处的地理纬度、地区的气候与气象
4、条件、安装部位及方式（贴附或架空）等。因此，在太阳能小屋的设计中，研究光伏电池在小屋外表面的优化铺设是很重要的问题。附件1-7提供了相关信息。请参考附件提供的数据，对下列三个问题，分别给出小屋外表面光伏电池的铺设方案，使小屋的全年太阳能光伏发电总量尽可能大，而单位发电量的费用尽可能小，并计算出小屋光伏电池35年寿命期内的发电总量、经济效益（当前民用电价按0.5元/kWh计算）及投资的回收年限。在求解每个问题时，都要求配有图示，给出小屋各外表面电池组件铺设分组阵列图形及组件连接方式（串、并联）示意图，也要给出电池组件分组阵列容量及选配逆变器规格列表。在同一表面采用两种或两种以上类型的光伏电池组件
5、时，同一型号的电池板可串联，而不同型号的电池板不可串联。在不同表面上，即使是相同型号的电池也不能进行串、并联连接。应注意分组连接方式及逆变器的选配。问题1：请根据山西省大同市的气象数据，仅考虑贴附安装方式，选定光伏电池组件，对小屋（见附件2）的部分外表面进行铺设，并根据电池组件分组数量和容量，选配相应的逆变器的容量和数量。问题2：电池板的朝向与倾角均会影响到光伏电池的工作效率，请选择架空方式安装光伏电池，重新考虑问题1。问题3：根据附件7给出的小屋建筑要求，请为大同市重新设计一个小屋，要求画出小屋的外形图，并对所设计小屋的外表面优化铺设光伏电池，给出铺设及分组连接方式，选配逆变器，计算相应结果
6、。二、 问题的分析该问题要求我们设计小屋外表面光伏电池的铺设方案，使得小屋的全年太阳能光伏发电总量尽可能大，而单位发电量的费用尽可能小。通过分析，我们发现以下三点：由于不同种类电池的发电量和价格不尽相同，因此我们需要对电池进行筛选，再进行铺设；总铺设费用包含电池的费用和逆变器的费用，为把问题简化，我们先考虑电池的选用与铺设，确定了电池的铺设方案后，再选择合适的逆变器，使得逆变器的费用尽可能小；由于小屋各表面形状、太阳辐射强度均不同，不同表面的电池也不能进行串并联，因此我们将各个表面分开考虑。针对问题一，电池仅考虑贴附安装方式，为寻找最优方案，在不考虑逆变器的情况下，我们首先引入电池单位面积收益
7、以及每块电池总收益这两个评价标准，并引入面积利用率这一概念，在进行铺设时，可以利用模拟退火算法的思想，先选择单位面积收益较高的电池铺设，观察面积的利用率，对于面积利用率不高的部分电池，用总收益更高的电池进行替换，最后用小面积电池填充，逐步优化，找到最佳铺设方案。确定了铺设方案后，再对光伏电池进行分组，确定串并联方式，同时考虑逆变器的分配，保证光伏分组阵列的端电压满足逆变器直流输入电压范围，光伏阵列的最大功率不超过逆变器的额定容量，综合考虑这几方面内容，进行最经济实用的分配,必要时电池可进行适当的替换。针对问题二，联系实际，我们不考虑墙面，仅考虑屋顶的电池架空安装。首先，我们确定电池架空时的倾斜
8、角，同时，考虑到电池架空后房屋表面会出现阴影，影响到后一排的电池受光，我们在安装电池时前后排应保持一定的间距，此间距可由电池架空时的倾斜角大致算得，在电池的选择和铺设方面仍然采用第一问的方法和思想，三、 模型假设与说明1、根据附件3，单晶硅电池的太阳光辐照阀值在太阳辐照强度低于200W/时， 电池转换效率转换效率的5%，假设此时电池转换效率为原转换效率的4%；2、根据相关资料，电池逆变器应避免安装在阳光直射处，否则可能会导致额外的逆变器内部温度，逆变器为保护内部元件而降额运行【1】，因此我们不必考虑逆变器所占用的面积；3、根据附件3，C薄膜电池C1、C2型号的太阳光辐照阀值在光照强度为200
9、W/是较1000 W/性能提高1%，我们将此理解为光照强度小于200W/且大于30W/时，电池转换效率提高1%；4、计算电池发电量时，当根据太阳辐射量计算出的电池发电量超过电池额定功率时，按照电池额定功率计算；5、假设任意放置的电池之间都可以进行串并联；6、假设屋顶前坡的梯形天窗为宽为1383.7mm的矩形；7、假设问题二中架空仅针对于屋顶，而墙面不考虑架空方式；8、假设问题二中太阳辐射不考虑散射量9、根据相关资料，倾斜面上的散射辐射量可由下式得到：其中，指倾斜面上散射辐射量，指水平面上散射辐射量，指倾斜面与水平面之间的夹角10、根据附件3的标注，所有光伏组件在0-10年效率按100%，10-
10、25年按照90%折算，25年后按80%折算，经换算，相当于光伏电池的寿命期为31.5年四、 符号说明及名词定义4.1 符号说明 ：房屋各面电池单位面积收益 ：房屋各面每块电池总收益 ：房屋各面电池每年单位面积发电量 ：房屋各面电池每年发电量 ：各面各时刻太阳辐射强度 ：各面全年太阳辐射总强度 ：单位功率电池价格 ：电池额定功率 ：电池转换效率 ：电池面积 ：小屋光伏电池35年寿命期内的发电总量 ：小屋光伏电池35年寿命期内的经济效益 ：每个面上每种电池的个数 ：电池投资的回收年限 ：电池架空时关于水平面的倾斜角 ：房屋屋顶前坡与水平面的夹角4.2 名词定义 1、电池收益：电池发电所产生的效益扣
11、除电池成本（不考虑逆变器）； 2、房屋各面电池单位面积收益：房屋各面电池单位面积发电量 民用电价-电池价格/电池面积； 3、房屋各面每块电池总收益：房屋各面电池发电量民用电价-电 池价格； 4、面积利用率：各表面上铺有电池的表面面积/表面总面积 五、 模型的建立与求解5.1 问题一 5.1.1 评价各类电池性能 首先，我们根据房屋各面的形状，太阳照射强度差异，将房屋表面分为东面，南面，西面，北面，屋顶前坡，后坡分别进行讨论。为简化模型，先不考虑逆变器的选用，只考虑电池的最优铺设。然后借助于C+，对附件所提供的数据进行分析和计算，求出各类电池每年在房屋各面单位面积收益及每块电池每年在房屋各面总收
12、益。具体算法如下：设在房屋各面电池每年单位面积发电量为，总发电量为，各面各时刻太阳辐射强度为，电池额定功率为，电池面积为则（单位） （单位）其中单晶硅电池当时，；单晶硅电池和多晶硅电池当时，；C薄膜电池C1、C2型号当时，；C薄膜电池当时，；所有型号电池当时，设房屋各面电池单位面积收益为，房屋各面每块电池总收益为，电池价格为（），民用电价为则 利用C+计算结果如下（具体程序见附录1）：产品型号每年单位面积发电量（)东南西北前坡后坡A168740.51474711176068307.7524588111936.3A267924.11457201162108209.0824296111794.5A
13、3763331637601305969225.3527303913254.6A467352.71444941152328140.0124091711695.3A561148.11311831046177390.1421872310617.9A661678.71323211055257454.2822062110710B18468116334812889121397.325232129323.7B285621.316516213032321634.925512329649.3B383479.416103112706321093.724874128907.6B477315.114914011768
14、019536.123037426773B583479.416103112706321093.724874128907.6B679404.715317112086120064.123660027496.6B778307.715105411919119786.923333127116.7C140571.773054.961133.617136.311041820681.9C235812.264484.853962151269746518255.7C3368576636655536.215567.310030818788.3C433896.861035.851075.81431792252.2172
15、79.3C537669.667829.256760.615910.510252019202.5C621069.437938.431747.58899.0957341.710740.4C721069.437938.431747.58899.0957341.710740.4C821243.638251.932009.98972.6457815.610829.1C921243.638251.932009.98972.6457815.610829.1C1023971.64316436120.410124.96524012219.8C1124784.244627.237344.810468.167451
16、.512634产品型号单位面积收益东南西北前坡后坡A1-1426.66 -186.65 -657.03 -2378.47 1363.30 -2321.32 A2-1428.39 -203.11 -667.89 -2368.91 1328.44 -2312.44 A3-1132.01 244.97 -277.37 -2188.95 1966.11 -2125.49 A4-1395.55 -180.58 -641.45 -2328.15 1338.09 -2272.15 A5-1269.43 -166.38 -584.80 -2116.12 1212.37 -2065.28 A6-1296.16
17、-183.54 -605.58 -2150.19 1207.18 -2098.91 B1-692.08 546.92 4.23 -1688.80 1948.25 -1563.96 B2-715.03 537.74 -10.97 -1722.81 1954.63 -1596.59 B3-470.74 750.70 215.70 -1453.31 2132.13 -1330.24 B4-626.31 504.93 9.44 -1536.33 1784.37 -1422.35 B5-489.00 732.44 197.45 -1471.57 2113.87 -1348.50 B6-649.80 51
18、2.01 3.13 -1584.42 1826.02 -1467.36 B7-640.15 505.60 3.76 -1561.86 1801.46 -1446.41 C1303.34 814.95 627.19 -65.77 1403.42 -9.92 C2267.63 719.22 553.49 -58.18 1238.66 -8.89 C3275.77 740.54 569.97 -59.54 1275.13 -8.81 C4253.36 680.79 523.92 -55.03 1172.45 -8.37 C5281.61 756.62 582.29 -61.10 1303.00 -9
19、.25 C6157.38 423.06 325.56 -34.31 728.67 -5.30 C7158.40 424.09 326.58 -33.28 729.69 -4.28 C8158.70 426.58 328.27 -34.57 734.71 -5.33 C9158.22 426.10 327.79 -35.04 734.23 -5.80 C10179.20 481.48 370.54 -38.89 829.18 -5.89 C11185.44 497.97 383.27 -40.04 857.45 -5.93 产品型号总收益东南西北前坡后坡A1-1821.33 -238.29 -8
20、38.79 -3036.46 1740.45 -2963.50 A2-2768.80 -393.71 -1294.64 -4591.88 2575.04 -4482.42 A3-1445.17 312.74 -354.10 -2794.51 2510.02 -2713.49 A4-2285.62 -295.74 -1050.56 -3813.03 2191.51 -3721.32 A5-2075.72 -272.06 -956.23 -3460.18 1982.41 -3377.05 A6-2512.47 -355.77 -1173.85 -4167.92 2340.00 -4068.53 B
21、1-1131.66 894.30 6.91 -2761.44 3185.68 -2557.31 B2-1386.00 1042.35 -21.27 -3339.49 3788.84 -3094.81 B3-692.05 1103.64 317.11 -2136.58 3134.54 -1955.65 B4-1018.93 821.47 15.36 -2499.42 2902.96 -2313.99 B5-948.82 1421.19 383.11 -2855.37 4101.66 -2616.57 B6-1260.85 993.49 6.08 -3074.33 3543.12 -2847.19
22、 B7-1067.78 843.34 6.27 -2605.18 3004.84 -2412.62 C1433.78 1165.38 896.88 -94.05 2006.89 -14.19 C2251.37 675.52 519.85 -54.64 1163.39 -8.35 C3434.40 1166.50 897.82 -93.79 2008.57 -13.87 C4390.17 1048.42 806.84 -84.74 1805.58 -12.89 C5433.68 1165.20 896.73 -94.09 2006.62 -14.24 C617.32 46.56 35.83 -3
23、.78 80.19 -0.58 C717.54 46.95 36.15 -3.68 80.78 -0.47 C834.65 93.13 71.67 -7.55 160.41 -1.16 C951.68 139.17 107.06 -11.45 239.80 -1.90 C1052.04 139.82 107.60 -11.29 240.78 -1.71 C11217.19 583.24 448.90 -46.89 1004.28 -6.94 5.1.2 确定电池铺设方案由此，我们可以筛选出收益较高的电池为各个表面铺设。我们利用模拟退火算法逐步找最优解的思想设计铺设方案，主要步骤如下：第一步，选
24、择几个单位面积收益相对较高的电池，分别对表面进行铺设，尽可能将表面铺满，计算每种电池铺设的面积利用率，选择其中面积利用率最高的；第二步，对面积利用率最高的铺设方案进行优化，将留有较大空隙的那部分电池换为每年总收益较之更高的电池进行铺设，将其尽可能铺满，若每年总收益更高的有几个，则按从大到小的次序分别尝试，选择面积利用率最高的；第三步，对铺设方案进行最终的优化，将表面上仍留有的空隙用面积较小的电池进行填充，使得表面利用率达到最大。各表面最终铺设方案如下：1屋顶前坡 由于屋顶前坡表面上有一个梯形天窗，但上下边长相差不大，因此我们将其看为宽为1383.7mm的矩形以简化问题。根据上述结果，前坡A3电
25、池的单位面积收益最高，因此，我们先用A3将屋顶前坡尽量铺满，其结果如下图1注：下列各图中打阴影部分均为门窗 图1此时屋顶表面左侧所留空隙较大，因此我们选择每年总收益比A3大的电池B1配合A3再重新进行铺设，结果如下图2（其中灰色矩形代表B1，白色矩形代表A3） 图2此时表面空隙较小，不能再铺设其他电池，为最优方案。2东面对于小屋东面，我们先选择单位面积收益最高的C1进行铺设，并与其他单位面积收益也较高的电池铺设方案进行对比，发现其面积利用率较高，定为初步方案，再用小面积电池C10填充，最优方案如下图3： 图33南面对于小屋南面，鉴于其表面形状的复杂性，我们将圆先看作是正方形，选择收益最高的电池
26、进行铺设，再经过层层优化，得到如下结果 图3然后，我们再将正方行替换为圆，继续填充空隙，得到最终结果如下图4 图44西面对于小屋西面，C1的单位面积收益最高，因此我们先选用C1对表面进行铺设，铺满后发现存在较多空隙，因此我们再用每年总收益更高的C2配合C1铺设，剩余部分再用面积较小的C7，C8，C10进行填补，其结果如下图5 图55北面及屋顶后坡根据附件所提供的数据，大同市位于北纬，小屋北面及屋顶后坡属于背阴面，受太阳辐射较少，由上述结果，我们也可以发现，各类电池受北向太阳辐射时，其单位面积每年收益及总收益均为负，因此北面和屋顶后坡并不适合安装光伏电池。 5.1.3 光伏电池分组及逆变器的选择
27、在确定了电池铺设方案后，我们对光伏电池进行分组和串并联，其中相同型号的电池进行串联，多个光伏组件串联后再进行并联（并联的光伏组件端电压相差不应超过10%），再为各组分配逆变器，并保证光伏分组阵列的端电压满足逆变器直流输入电压范围，光伏阵列的最大功率不超过逆变器的额定容量，综合考虑这几方面内容，进行最经济实用的分配,必要时电池可进行适当的替换。分配结果如下：1屋顶前坡经过分配，我们将A3分为4组，B1作为一组，每组电池先串联，再将5个组并联，此时光伏组件端电压相差不超过10%，满足条件，总容量为9250W，按此规格，我们可以选取一个型号为SN17的逆变器与之相连，具体结果见下表1组别电池型号个数连接方式端电压/V连接方式逆变器型号分组阵列容量/W屋顶前坡第一组A38串联368.8并联SN179250A38串联368.8A38串联368.8A39串联414.9B110串联379.1 表1分组阵列图形及组件连接方式如下图6 图62东面经分配，我们将C1分为13组，C10分为两组，并联，具体配备结果如下表2组别电池型号个数连接方式端电压/V连接方式逆变器型号分组阵列容量/W东侧第
