2012年全国大学生数学建模大赛B题--论文.docx

2012高教社杯全国大学生数学建模竞赛B题 太阳能小屋的设计摘要：在太阳能小屋的设计中为实现太阳能光伏板最佳朝向、倾角及排布阵列设计及优化，通过建立倾斜放置的光伏板表面接收太阳辐射能模型，计算到达光伏板上的太阳辐射能量，推导出光伏板的最佳朝向及倾角。为使光伏板最大限度地接收太阳辐射的能量，在选择合适的朝向及倾角的基础上，对光伏电池排布阵列，建立目标规划，并通过与实际逆变器的相互匹配，不断对目标进行优化，最终得到一组最优解。通过上述研究，结合山西大同市本地情况，重新设计出一个更加适合当地地理及气象条件的太阳能光能房屋并为其选择最优的阵列排布方案。针对问题一：电池板只是铺设房屋的表面，没有涉及到电池板放的角度问题，先求算出房屋的角度为10.62度，再根据角度，建立模型算出光伏板上太阳能辐射量。并用目标规划阵列排列方案计算出电池的排布。再通过排布计算出经济效益，最后得出35年之内无法收回成本。针对问题二：通过对角度建立模型，计算得出最佳角度44.66度，通过排布计算出电池板排布最佳方案，建立模型计算出经济效益，在28.5年收回成本。如考虑货币时间价值，35年的经济效益是亏损的。针对问题三：要通过目标构建一个产电量尽量大，而成本尽量小的理想模型。假设小屋无挑檐、挑雨棚（即房顶的边投影与房体的长宽投影相等），建立模型计算出最佳的图形，并画出模型图。关键字：太阳能 太阳能辐射模型 最佳倾角 电池 模型 目标规划一.阐述问题太阳能作为迄今人类所认识的最清洁的可再生能源，其与建筑一体化将在建筑节能中起到十分重要的作用。屋顶在建筑外围结构中所接受的日照时间最长，接受的太阳辐射量也最大，具有利用太阳辐射的优越条件，同时，屋顶较开阔，便于大面积连续布置太阳能设备，因此，在城市中，建筑屋顶是太阳能利用的最佳场所。目前，许多国家已纷纷实施和推广“太阳能屋顶计划”，如有德国十万屋顶计划、美国百万屋顶计划以及日本的新阳光计划等。我国属于太阳能利用条件较好的地区，尤其是青藏高原地区太阳能。二.问题的分析图1 流程图首先，根据附件给出的数据进行分析和运算。已知山西大同市的地理条件（经纬度、海拔等等）、全年每个小时水平面总辐射强度、水平面散射辐射强度、法向直射辐射强度、各个方向总辐射强度等，建立太阳辐射几何模型，再由几何模型建立数学模型，得到各个表面每个小时太阳能光伏板表面的辐射量和中间参数的关系，将已知参数带入关系式得到各个表面每个小时太阳能光伏板表面的辐射量，又已知太阳能光伏板在不同年限内的转换效率（所有光伏组件在010年效率按100%，1025年按照90%折算，25年后按80%折算），由此得到各个表面每个小时太阳能光伏板的发电量。根据已知太阳能光伏板在各表面每个小时发电量的条件下，考虑小屋各外表面电池组件铺设分组阵列及组件连接方式（串、并联），做好初步预算并通过光伏电池和逆配器的参数等等数据对阵列排布进行目标规划，以求得最大的发电量及最小的费用，并给出电池组件分组阵列容量及选配逆变器规格列表。第一小题，仅仅需要考虑贴附安装方式，选定光伏电池组件，对小屋的部分外表面进行铺设。首先在已建立的数学关系模型中带入贴俯安装方式下各个变量的关系式和各个定量的值，得到此安装方式下各种型号的太阳能光伏板在各个表面的发电量，然后对各种型号的太阳能光伏板进行优化组合，得到一种最优的组合阵列，并得出电池组件分组阵列容量及选配逆变器规格列表。第二小题，着重考虑电池板的朝向与倾角，也就是在第一个问题的基础上，对太阳能光伏板的朝向与角度进行优化。首先考虑架空安装方式下太阳能伏板阵列排布的朝向，由于太阳总辐射中的散射部分与阵列朝向无关，所以只需要考虑阵列上太阳直射辐射强度随阵列面朝向的变化即可。光伏板的放置朝向可分为朝向赤道和任意方向两种，在相同倾角的情况下，赤道斜面接收的辐射能量要大于任意方向辐射接受的能量，所以在此仅讨论第一种情况，即电池板朝向赤道摆放时最佳倾角的选择及阵列的排布方式。题里给出地理纬度、地形高度等参数以后，倾角为的太阳能伏板表面1年内接受的总辐射Q是一个关于变量的函数Q（），对Q（）关于变量求导并取值为0，即dQ()d=0求解方程，即可得到年最佳倾角。至此，朝向及倾角问题已解决，对于太阳能伏板的选择及优化与第一题相同，对组合阵列进行目标规划，得到一种最优组合。第三题，根据附件7给出的小屋建筑要求，为大同市重新设计一个小屋。综合前前两道小题的结论，本题以原有小屋为基础，尽量提高房屋的有效利用面积，并依据已有模型给出小屋太阳能伏板的最佳朝向及最佳倾角，从而得到一个高效太阳能小屋。三.模型的假设1.模型的建立11几何关系建立及基本参数的确定S东=S总-S门=24.23m2S南=S总-S圆-S窗-S车库=21.78m2S西=26.98m2S北=S总-S窗-S门=28.119m2S南（顶）=60.785097m2S北（顶）=14.031368m2以上全部根据题中给出数据计算。四.模型的建立与求解1问题一的模型11几何关系建立及基本参数的确定首先对房屋顶面建立几何模型：图2图3由几何模型可以推导出 1=tan-1316cos=S1S2=90-1根据附件中山西大同市实地测量数据，结合已建立的几何模型，建立倾斜放置的光伏板表面太阳辐射量数学模型，系统框图如图1所示。图4 倾斜放置的光伏板表面太阳辐射量数学模型建立已知山西大同市的地理参数（纬度、地形高度等）以及中间参数（日序、光伏板倾角、方位角和时刻），可以得到逐时太阳能光伏板表面的辐射量和中间参数的关系。将逐时太阳能光伏板表面辐射量关于时间积分得到某一天的日辐射总量Wd（N,），再将Wd（N,）关于N累加得到太阳能光伏板表面的年累计辐射量Wy（）。计算地球表面任一点的太阳辐射量，首先确定一些基本的天文参数，主要包括地球表层大气外界上空的垂直太阳辐射强度、赤纬角、太阳高度角、太阳方位角和日出日落时刻等。由几何模型可以推导出 1=tan-1316cos=S1S2=90-1根据附件中山西大同市实地测量数据，结合已建立的几何模型，建立倾斜放置的光伏板表面太阳辐射量数学模型，系统框图如图1所示。图5 倾斜放置的光伏板表面太阳辐射量数学模型建立已知山西大同市的地理参数（纬度、地形高度等）以及中间参数（日序、光伏板倾角、方位角和时刻），可以得到逐时太阳能光伏板表面的辐射量和中间参数的关系。将逐时太阳能光伏板表面辐射量关于时间积分得到某一天的日辐射总量Wd（N,），再将Wd（N,）关于N累加得到太阳能光伏板表面的年累计辐射量Wy（）。计算地球表面任一点的太阳辐射量，首先确定一些基本的天文参数，主要包括地球表层大气外界上空的垂直太阳辐射强度、赤纬角、太阳高度角、太阳方位角和日出日落时刻等。1.1.1太阳赤纬角（）日地中心连线与赤道的夹角称为赤纬角，赤道以北为正、南为负，近似等于：=223.45360sin(2284+N365) (度)其中为日期序号，例如，1月1日为，3月22日为。（该定义摘自太阳能应用技术的第二章太阳辐射）1.1.2日出、日落时刻（Tsr，Tss）在北半球，纬度为、朝向正南、与水平面成倾斜角的太阳能电池板上，太阳光的入射角与纬度为（-）的水平面上太阳光入射角是相等的，即cos=sin（-）sin+cos（-）coscosee为倾斜面上日出、日落时角，e=cos-1（-tan（-）tan）由上可以得到倾斜面上日出日落时刻：Tsr=12（1-e15）Tss=12（1+e15）1.1.3太阳高度角太阳高度角是太阳相对于地平线的高度角，这是以太阳视盘面的几何中心和理想地平线所夹的角度。太阳高度角可以使用下面的算式，经由计算得到很好的近似值：其中为太阳高度角，为时角，为当时的太阳赤纬，为当地的纬度(大同的纬度为)。（该定义摘自维基百科）1.1.4太阳方位角（A）太阳方位角是太阳在方位上的角度，它通常被定义为从北方沿着地平线顺时针量度的角。它可以利用下面的公式，经由计算得到良好的近似值，但是因为反正弦值，也就是有两个以上的解。下面的两个公式也可以用来计算近似的太阳方位角，不过因为公式是使用余弦函数，所以方位角永远是正值，因此，角度永远被解释为小于180度，而必须依据时角来修正。当时角为负值时 (上午)，方位角的角度小于180度，时角为正值时 (下午)，方位角应该大于180度，即要取补角的值。 其中为太阳的方位角，为太阳高度角，为时角，为当时的太阳赤纬，为当地的地理纬度。（该定义摘自维基百科）1.2太阳能光伏板上太阳能总辐射量的计算光伏板的放置方式可分为朝向赤道和任意方向两种，在相同倾角的情况下，前者斜面接收的辐射能量要大于后者，所以在此仅讨论第一种情况。某一时刻，太阳能伏板上接收的总辐射能量主要由法向直接辐射能量（W1）和散射辐射（W2）组成，即：W总=W1+W2式中，电池板上太阳瞬时直接辐射能量:W1=Wocos电池板上瞬时太阳散射辐射:W2=Wpcos22sin一天内，太阳能伏板接收的太阳能总辐射能量为：Wd=TssTsrW总dT=Wd1+Wd2式中，Wd1、Wd2分别表示日直接辐射能量和日直接散射能量Wd1=TssTsrW1dTWd2=TssTsrW2dT年辐射总量为：Wy=1365Wd各电池光伏板发电量（w）：Q=W总10001000x1光伏电池实际价值： P=Q1000100010001p35年寿命期内的发电总量(Q35)根据折旧要求，所有光伏组件在010年效率按100%，1025年按照90%折算，25年后按80%折算，即：Q35=Qy10+90%Qy15+80%Qy1035年寿命期内的经济效益为应为发电总价值光伏阵列的总造价，即：S=Q3510000.5-P总x为电池光伏板的面积，p（元/Wp）光伏电池的价格，1为光伏电池的转换效率。1,3 贴附安装方式下太阳能伏板阵列排布方案及太阳总辐射能计算根据山西省大同市的气象数据，贴附安装方式下，=10.62。带入公式计算可得此条件下，房屋全年接收的总太阳辐射能量为W总=91607153.50W/m2。已知各型号光伏电池的面积和转换效率，则根据W总可以得到各电池光伏板的实际发电量以及实际成本。此时，将考虑如何排布太阳能伏板阵列。我们为了电池的排布达到最佳状态，建立阵列排布，来找电池最佳排布和逆电器的型号确定。maxz1=k=16j=111i=13WijkXijk minz2=k=16j=111i=13XijkPij j=111i=13XijkSij0.9Sk (k=1,2,3,4,5,6) Xijk为整数其中Xijk为第i种类型j 种规格电池在k 面上的数量 Wijk 为1块 i类型j 规格电池在k 面上的接收辐射量Pij 第i种类型j 种规格电池的价格 Sij第i种类型j 种规格电池的面积 表示Sk 表示k 面的面积 当k=1时 表示顶面当k=2时 表示北斜面 当k=3时 表示东面 当k=4时 表示南面 当k=5时 表示西面 当k=6时 表示北面 maxz1=k=16j=111i=13PijXijmm minz2=k=16fmymk 0.8PijXijmPmymk 0.8PijXijmPmymk 0.8XijmUijUm0.8XijmUijUmXijm为第i种类型j 种规格电池对应第m 种逆变器的电池数量 ；ymk 第m 种逆变器在k 面上的数量；fm第m 种逆变器的单价；Pij第i种类型j 种规格电池的功率；Pm 第m 种逆变器允许通过的最大功率；Pm 第m 种逆变器允许通过的最小功率；Uij 第i种类型j 种规格电池的电压；Um 第m 种逆变器允许通过的最大电压；Um 第m 种逆变器允许通过的最小电压；通过软件计算得出：南顶：使用40个A3单晶硅电池，四组10个A3单晶硅电池串联，四组串联单晶硅电池并联，并用一个SN17和SN1逆变器，如下图排列：图6南顶电池摆放位置北顶：使用7个A3单晶硅电池，并且每个电池都是并联的，用一个SN4逆变器如下图：图7北顶电池摆放位置东墙：使用12个C1薄膜电池，2个串联然后再把每个电池组并联，用一个SN5逆变器如下图：图8东墙电池摆放位置南墙：使用12个C7薄膜电池，4个串联然后再把每个电池并联，用一个SN3逆变器，如下图：图9 南墙电池摆放位置西墙：使用59 个C9薄膜电池，6个串联然后再把每个电池并联，用一个SN11逆变器，如下图：图10 西墙电池摆放位置北墙：使用17个C1薄膜电池，3个串联然后再把每个电池串联，用一个SN12逆变器，如下图：图11 北墙电池摆放位置1,3 经济效益35年寿命期内的发电总量根据折旧要求，所有光伏组件在010年效率按100%，1025年按照90%折算，25年后按80%折算，即:Q35=Qy10+90%Qy15+80%Qy10经济效益为：S=Q3510000.5-P总通过计算得到，35年寿命之内还是无法收回成本。2.问题二：架空安装方式下太阳能伏板阵列排布方案及太阳总辐射能计算2.1最佳方位角的确定对于全天无阴影遮盖的太阳能电池阵列，如果其倾角固定，则必然存在一个能够得到全天最多太阳总辐射能的最佳朝向，即最佳方位角。由于太阳总辐射中的散射部分与阵列朝向无关，所以只需要考虑阵列上太阳直射辐射强度随阵列面朝向的变化即可。2.2最佳倾角的确定给定地理纬度、地形高度等参数以后，倾角为的太阳能伏板表面1年内接受的总辐射Q是一个关于变量的函数Q（），对Q（）关于变量求导并取值为0，即：dQ()d=0求解方程，即可得到年最佳倾角。2.3最佳倾角的计算用Matlab编程将量化成0到90，增量为1，计算在各个值下的Q（）的值，从中挑选出使Q（）最大的为。下图所示为最佳倾角计算流程。图12根据计算从0到90的计算，最后算的最佳倾斜角度为44.66。2.4最佳朝向及倾角条件下的阵列排布方案在最佳倾角和最佳角度的时候，我们为了电池的排布达到最佳状态，建立阵列排布，来找电池最佳排布和逆电器的型号确定。maxz1=k=16j=111i=13WijkXijk minz2=k=16j=111i=13XijkPij j=111i=13XijkSij0.9Sk (k=1,2,3,4,5,6)Xijk为整数其中Xijk为第i种类型j 种规格电池在k 面上的数量； Wijk 为1块 i类型j 规格电池在k 面上的接收辐射量；Pij 第i种类型j 种规格电池的价格； Sij第i种类型j 种规格电池的面积； Sk 表示k面的面积，当k=1时， 表示顶面；当k=2时，表示北斜面； 当k=3时 表示东面；当k=4时，表示南面；当k=5时，表示西面；当k=6时，表示北面。maxz1=k=16j=111i=13PijXijmmminz2=k=16fmymk 0.8PijXijmPmymk 0.8PijXijmPmymk 0.8XijmUijUm0.8XijmUijUmXijm为第i种类型j 种规格电池对应第m 种逆变器的电池数量； ymk 第m 种逆变器在k 面上的数量；fm第m 种逆变器的单价；Pij第i种类型j 种规格电池的功率；Pm 第m 种逆变器允许通过的最大功率；Pm 第m 种逆变器允许通过的最大功率；Uij 第i种类型j 种规格电池的电压；Um 第m 种逆变器允许通过的最大电压；Um 第m 种逆变器允许通过的最大电压；通过计算的得到电池排布和第一问一样，但是只是屋顶电池板铺设时，角度与以前发生了改变。2,5经济效率35年寿命期内的发电总量根据折旧要求，所有光伏组件在010年效率按100%，1025年按照90%折算，25年后按80%折算，即:Q35=Qy10+90%Qy15+80%Qy10经济效益为：S=Q3510000.5-P总通过计算得到，28.5年收回成本。如考虑货币时间价值，35年之内是无法收回成本的，是亏损的，但是使用太阳能节约了能源减少了污染。3.问题三模型的建立与求解3.1.建立模型通过对问题三的分析，我们的最终目标是需要重新构建一个产电量尽量大，而成本尽量小的理想模型。假设所建小屋无挑檐、挑雨棚（即房顶的边投影与房体的长宽投影相等），之后我们可以把整个模型看成几何立体图形，不考虑其墙体的厚度进行模型的建立。a15 b3 ab（隐含条件：最短边一定比长边小）wl741）假设南北立面所映射的边为长边，东西立面所映射的边为短边。j=0xijyij0.35ah (i=1,2)j=0xijyij0.5bh (i=3)j=0xijyij0.3ah (i=4)2）假设南北立面所映射的边为短边，东西立面所映射的边为长边。j=0xijyij0.35bh (i=1,2)j=0xijyij0.5ah (i=3)j=0xijyij0.3bh (i=4) a ，b