核电安全分析与预测论文(完整版).doc

核电安全分析与预测论文摘要本文基于核电站的安全问题，建立了发展预测模型。该模型主要考虑到核电站的需求供给关系、技术转化率和环境效益以及核电安全风险。其中包含两个模型的研究，核泄漏放射性物质连续点源高斯优化模型模型A，核电站供求均衡评估模型模型B。模型A关键在于建立空间中任意一点的浓度的模型，确定Z=0时的模型，并根据人体承受浓度的标准，考虑地面反射、泄漏源的有效高度、放射性物质自身的沉降作用和雨水吸附作用等对放射性物质浓度影响，确定关于x ,y的函数模型，最后通过编程进行测算和检验。模型B是从核电站资本成本、风险成本的投入得出核电站资本存量，利用柯布道格拉斯生产模型，结合技术转化率建立出核电站供给模型；又利用核电发电缺口，结合技术转化率、GDP发展水平和环保效益建立核电站需求模型。利用供需均衡，求解相应的稳态数量。模型A解决问题二。模型B具有一定普适性，结合中国实际情况，可以有效解决问题一和问题三。可以直观看出建立核电站必要性，并对核电安全发展进行相关预测。关键词：连续点源高斯优化模型 核电站供求均衡评估模型 柯布-道格拉斯函数目录1、问题的提出42、问题的分析53、模型假设63.1、核泄露放射性物质泄露模型模型A63.2、核电站供需均衡模型模型B64、模型的建立与求解74.1.1、模型A设计74.1.2、模型A的建立84.1.2.1、模型A的建立84.1.2.2、模型A的修正94.2.1、模型B设计154.2.2、模型B变量解释164.2.3、模型B建立174.2.3.1、供给模型174.2.3.2、需求模型184.2.3.3、均衡模型184.2.4、模型B分析184.2.4.1三维模型分析184.2.4.2控制时间因素的技术转化率因素分析204.2.4.3控制技术转化率因素的时间因素分析204.3、模型求解214.3.1问题一214.3.2问题二224.3.3问题三285、模型评价30模型A：30模型B：316、参考文献327、附录33附录一，模型B33附录二，模型A341、问题的提出我国是一个电力消耗大国，经济的快速发展和人们生活水平的日益提高，对我国的电力生产提出了严峻的挑战。中电联副秘书长欧阳昌裕在“2011中国清洁电力论坛上”表示，“十二五”期间，中国电力需求将保持8.5%的年均增长速度，每年要增加5000亿度电，而未来二十年乃至三十年，中国都将保持较高的电力需求增长。由于煤炭资源的不可再生性以及煤电的污染问题，核电似乎是我们为满足电力需求的所应作出的比较理智的选择。然而，今年日本福岛核电站的核泄漏问题表明核电并不是一只温顺的小绵羊，该事故引起了全世界对核电站及其安全的重新思考。如德国已经宣布2022年前关闭所有核电站。请从互联网或报刊上搜集有关数据，根据这些数据建立评估核电站安全的数学模型。考虑：（1）随着人们生活水平的提高，用电量大幅增加，假设不建设核电站，用电量和发电量之间的差距有多大？建设一个某种规格的核电站能提供多少电力？建设核电站的经济成本和效益如何？目前国内有几个核电站、在建或准备建设的有几个？也就是要求建立建设核电站必要性的数学模型并分析；（2）以秦山或大亚湾核电站为例，如果这些地方出现了严重的自然灾害造成了核泄漏（需要你自己作出合理假设），那么，在一定气象条件（一定风向、风力、下雨等）的情况下建立核扩散的数学模型，并讨论对周围多大范围的居民进行疏散以及其他的应对措施和可能的后果；（3）综合前两问给出核电站的发展前景预测。2、问题的分析核电虽然是高效的清洁能源，但一旦发生核电站泄漏事故，其带来的经济损失和环境损失是不可估量的。对于其损失的大小，我们利用其各种风险概率计算出相应的事故率和期望损失值。而对其事故造成影响的研究，首要问题就是对于核泄漏影响区域的模拟。为此我们建立了核泄漏放射性物质连续点源高斯优化模型模型A。模型I采用了烟雾扩散模型，通过量化固定风向风力和降雨等因素，定量分析了发射性物质可能的扩散范围。在建模的过程中，根据题设和合理的假设，运用概率学知识，通过图形解析和数学理论推导得出“连续点源放射性物质高斯扩散模型”。并进一步考虑地面反射、泄漏源的有效高度、放射性物质自身的沉降作用和雨水吸附作用等对放射性物质浓度影响，进行模型建立。基于此，我们定性将讨论在污染范围内所需采取的补救措施的可行性，已经经济成本。从而对核电站安全性进行评估。对于核电站的安全和发展预测问题，我们分为三个层次：首先是利用“核泄漏放射性物质连续点源高斯优化模型”模拟出核泄漏可能污染的区域，从而得出损失成本；其次，建立核电供需均衡模型，并且利用分离变量法，和别研究技术转化率和时间对于核电站稳态数量的影响模型B;最后基于模型B，结合中国实际国情，探讨并预测核电未来的发展趋势。我们认为，随着技术进步，科技的发展，核电技术会经历成熟到衰落的阶段。在较高的科技进步水平下，最终有可能会被其他更加清洁和安全的能源所取代。3、模型假设3.1、核泄露放射性物质泄露模型模型A1、放射性气体初始泄漏时可看作在空中某一个连续点源向四周等强度的瞬时释放，在于风向垂直的三维坐标的两个方向上轴和轴的分布呈高斯分布。2、放射性物质在无穷空间扩散过程中不发生化学反应，且不计地形影响。3、扩散系数呈各向同性，且放射性物质扩散服从扩散定律。4、放射性物质在穿过降雨区域时，其强度由于雨水的吸收而减少，减少比率为常数。5、假设地面对放射性物质起不完全反射作用，存在放射性物质的沉积消耗。6、假设风向为水平风向，且风向、风速不随时间和地点变化。风速大于1m/s。7、扩散过程中不考虑泄漏源点内部温度的变化对烟气气体扩散的影响。3.2、核电站供需均衡模型模型B1、排除通货膨胀影响因素，以及假设现金回报率为0。2、以大亚湾核电站的运营数据为标准核电站数据。3、技术进步表现为技术转化率，体现核电能源和新兴绿色能源的替代转化。4、核电的发展符合规模经济。4、模型的建立与求解核泄露放射性物质泄露模型模型A4.1.1、模型设计空间任意一点的放射性物质的扩散系数空间任意一点的放射性物质浓度泄漏源泄漏的放射性物质总量源强，用浓度标准差表示的轴上的扩散参数任意扩散时刻泄漏源距地面的实际高度烟云抬升高度泄漏源的有效高度泄漏源出口的有效内径泄漏源高度处的平均风速，m/ s放射性物质的扩散速度, m/ s沉降速度，空气的动力粘性系数地面反射系数大气压力，雨水吸附系数放射性核素的半衰期太阳高度角浮力通量， m4 / s泄漏源的热排放率，泄漏源出口处温度，环境大气平均温度，泄漏源出口温度与环境温度差，放射性物质粒子密度，重力加速度，9. 806 54.1.2模型的建立4.1.2.1模型的建立当风速为时，利用连续点源高斯扩散模型分析核电站周边放射性物质浓度的变化情况。此泄漏点源是有边界点源，泄漏点源的实际高度为。以泄漏点源在地面的投影点为坐标原点，以风向方向为轴，铅直方向为轴，与轴水平面垂直方向为轴建立三维坐标系，由于扩散过程中浓度在、轴上的变化分布符合高斯分布，所以下风向的任意一点的浓度函数为：（1）根据概率统计我们可以得出方差的表达式为：（2）进而，源强的积分公式可以根据假设得出：（3）把（13）式代入（14）积分可以得出：（4）将(13)式和（16）式代入（15）式可以得出：（5） 最后再将把（16）、（17）式代入（13）式可以得出：（6）4.1.2.2模型的修正（1）、地面对放射性物质的发射作用考虑到地面会对放射性烟流的扩散有一定的反射作用，同时扩散的核素粒子受沉降等作用，核扩散物质又不可能被全部反射回去，因而核辐射物质只能是部分反射回大气，为了便于描述和模拟，设地面反射系数为。这样进入大气的核辐射物质可以看成是两个部分：一是从泄漏源直接扩散到空间A点；二是从地面反射进入空间A点（见图1）。.图1高架连续点源扩散的地面反射图从几何物理学分析可知，通过地面反射进入大气的核扩散物质相当于虚拟泄漏源泄露的核辐射物质在原来空间的一个浓度叠加。如果所设核泄漏点源在距有效地面的高度为的地方，实际情况中下风向任一点的浓度值：1、实际泄漏源的贡献为： （7）2、虚像泄漏源的贡献为： （8）则该处的实际浓度为： （9）综合上面的公式得到连续点源高斯烟羽扩散模型： （10）（2）、泄漏源的有效高度对模型的修正1、 考虑热力抬升作用对模型的修正如图 2 所示，为核泄漏点源的有效源高。它是由两部分构成：一是核泄漏口的有效高度；二是在实际核扩散中核泄漏气团从泄漏口排出时，由于受到热力抬升和本身动力抬升，进而产生的一个附加高度，有。图2 热力抬升示意图对于，主要由浮升力和泄漏的初始动量决定，同时还要受到泄漏口温度、大气温度、风速、地形地貌等多种因素的影响。我们直接引用气体污染扩散学中应用较广泛的，有关烟气抬升高度的综合分析公式：由 Briggs 抬升公式知，浮力通量: 并且，根据该结论的分析，带入风速及放射性气体的扩散速度有三种情况：1） ， 适用于k=1.0m / s且=35K2）， 适用于k=1.0m / s且=35K3） , 适用于k不限（k 1.0时，取k =1.0），=1.0m / s且=35K2） ， 适用于k=1.0m / s且=35K3） , 适用于k不限（k 1.0时，取k =1.0），在计算是，也将注意x的符号），则上风向的放射性气体的扩散浓度为其他具体的参数表达形式可参考上页。核电站供求均衡评估模型B：4.2.1、模型设计本模型拟用核电站的稳态数量与现实情况来衡量核电站的存在必要性。假设核电站的供给主要受以下因素影响：1、建设核电站的投资成本、运营成本、废弃成本、风险成本2、技术转化率3、核电站数量4、观察时间核电站需求主要受以下因素影响：1、电力缺口（标准发电量、GDP增速、核电站数量）2、环保效益3、技术转化率4、观察时间4.2.2、变量解释C：核电投资资本存量I：单位核电站投资成本c：单位核电站运营成本w：观察期内单位核电站废弃成本p：观察期内单位核电站风险成本N：观察期内核电站数量t：观察期：资本存量的弹性系数k：观察期内的年均技术转化率AS：核电总供给Qn：单位核电站发电量（千瓦时）Pe：单位电价（元/千瓦时）i：观察期内GDP年均增长率r:环境治理投资占GDP的比重Ep:环保系数D：不建核电站的电力缺口AD：核电总需求4.2.3、模型建立本模型通过研究核电需求函数与核电供给函数，分析得出其稳态数量。如果稳态数量等于或大于目前我国已建、在建和筹建的核电站数目，那么认为核电站有继续建设的必要；反之，则没有必要。模型为：4.2.3.1、供给模型本模型认为所有核电站供给符合柯布道格拉斯的生产函数模型。（1）、资本存量C：资本存量等于投资成本和运营成本、废弃成本和风险成本之和。运营成本、废弃成本和风险成本与时间有关。其中，废弃成本是将废弃处理费用摊销到观察期所得，风险成本根据模型A与事故概率计算得出期望损失值。（2）、总供给AS：柯布道格拉斯生产函数，其中，将技术水平因素替换为，资本存量替换为C，可根据美国等核电技术成熟国家通过拟合算出，约为0.866。由此，总供给函数整理得：4.2.3.2、需求模型本模型考虑到发电缺口与技术转化的关系（1）不建核电站的发电缺口（2）总需求AD：技术进步会产生核能的不断被替代，所以在总需求模型中考虑到该部分被转化替代的核能。4.2.3.3 均衡模型本模型认为核电稳态均衡符合AS=AD，即：整理出关于N的二元函数，得.（*）（*）式是求出在不同观察期和技术转化率水平下的均衡N值，下面我们将利用计算机实现模型B，考察其函数特点。然后再用控制变量法，分别隔离时间因素和技术进步因素，计算不同控制条件下的核电站的稳态数目。4.2.4、模型分析4.2.4.1三维模型分析该模型的求解通过计算机实现，得出稳态的核电站数目。通过查阅资料得到以下数据：（1）Qn=30000000000（千瓦时）单座核电厂年发电量；（2）工业交通建筑用电占84.3%，均价0.65元/千瓦时；居民用电占7.7%,均价0.5元/千瓦时；农业用电占6.8%，均价0.5元/千瓦时；商业用电占1.2%，均价0.8元/千瓦时则用电均价Pe为：0.650.843+0.50.077+0.50.068+0.80.0120.6（3）r=2.00% 为国家环境治理投资占GDP的比重，Ep=1（4）I= 27900000000（元人民币）建设一座核电站前期所需投资的成本（元）（5）w=3500000000（元人民币）一座核电站处理核废料所需投资成本（元）（6）c=20000000（元人民币）一座核电站正常运营成本（元）（7）p=15000000（元人民币）一座核电站事故风险成本（元）（8）GDP平均增长率为9%（9）=0.866用matlab编程，代入以上数据进行模型求解，具体程序见附录模型求解的图像如下：由图形可以直观看出，当技术进步率水平不高时，可以看出需要核电站数量较大。就目前实际来看，现阶段技术进步率水平处于较低水平，可以预见到核电的缺口较大，所以建设核电站是必要的。技术进步率水平适中的时候，核电站的数量水平比较稳定。当技术水平处于较高阶段，核电有可能被其他能源替代。4.2.4.2控制时间因素的技术转化率因素分析假设期为30年，技术转化率在0-9%，该模型的求解通过计算机实现，得出不同技术转化率下稳态的核电站数目的图像：由图直观看出，稳态核电站数目是单调递减函数，在技术转化率低的时候，用电量缺口较大，只能用现有的核电技术解决较大的用电缺口问题；在技术转化率较高的情况下，用电缺口将由更加高效清洁的能源所替代。现在全球的技术转化率大约为1%左右，从图中可以看到，稳态核电站数量为603座，而2006年实际核电站数量为428座，所以有必要建设。而近几年来，各国都在大力发展核电事业，从这点，模型与事实具有较高的拟合度。4.2.4.3控制技术转化率因素的时间因素分析假定的考察国家为美国，其技术转化率为2.5%，该模型的求解通过计算机实现，得出不同时间下稳态的核电站数目的图像：由图直观看出，稳态核电站数量呈先逐渐减少后逐渐稳定的图像。在开始的10年中，随着时间的增加，在一定的技术转化率水平下，一部分核电需求逐渐减少，被新能源所替代。新能源最终与核电达到平衡，达到了稳态数量的稳定。美国现阶段技术转化率大约为2.5%，其实际核电站数量为114座，并且暂缓核电站的建设。根据模型的预测约为112座，已经达到了该技术水平下的稳态数量，与美国实际数量较接近，不宜再增加核电站数量，该情况较符合事实，模型具有较高的拟合度。4.3、模型求解4.3.1问题一（1）对于用电量和发电量差距的问题，我们通过查找相应的数据和资料。可以看出：我国供电量和用电量基本保持均衡。随着我国经济的发展，GDP增速的提高，我国的用电需求在逐年增大。除去核电站的发电量，我国的用电缺口量将呈上升趋势。（2）通过一定的搜索和整理，我们找到了目前我国目前已有7个投入运营的核电，11个在建的核电站，25个筹建的核电站。通过计算得到平均的发电能力为594.24万千瓦。（3）建设核电站的经济成本和效益，根据模型B的供给模型中可以求得其成本函数，为具体含义和变量分析见4.2模型分析。其效益函数等同于总需求函数，即效益分析见4.2模型分析。（4）我国核电站建设的必要性。现阶段我国的技术转化率约为3%，通过（*）式，可以计算出得我国核电站需求稳态数目为65座。而现阶段我国已建、在建和筹建的核电站数量为43座，可以看出，还有较大的缺口电量，所以我国有必要建设数量规模适中的核电站。4.3.2问题二（一）、模型参数的确定模型所需参数的选取对模型的求解至关重要，通常情况下气象参数的选取是利用该地区多年气象资料，采取工业安全与环保统计的方法进行有关参数的确定，而其他扩散参数则以实际测定为准。1、大气稳定度的计算：根据我国标准(GB/ T 1320191) )制订地方大气污染物排放标准的技术方法的规定，大气稳定度分为6级，分别为A极不稳定、B不稳定、C弱不稳定、D中性、E弱稳定、F稳定。该方法的技术路线是：根据核泄漏源所在地的经度和纬度以及泄漏的日期和时间计算当时的太阳高度角，利用天气条件确定辐射等级，然后利用辐射等级和风速确定大气稳定度，最后查扩散参数幂函数表，确定扩散参数，基本过程，如图3。 图3-中国国家标准规定大气稳定度基本理念首先，然后，由太阳高度角和云量查出太阳辐射等级；最后，再根据地面风速确定当时的大气稳定度。天空云层的情况日照角60度日照角35度日照角15度天空云量为4/8或高空有薄云强中等弱天空云量为5/8-7/8，云层高度为2134-4877m中等弱弱天空云量为5/8-7/8，或云层高度2134m弱弱弱表一：日照强度的确定地面风速白天日照夜间条件强中等弱云量为4/8云量为3/86CDDDD表二： Pasquill大气稳定度确定2、扩散系数的计算利用上述途径确定大气稳定度，进而确定扩散参数的计算方程式。扩散参数的确定，采用给出一套扩散参数幂函数表，如表三和表四：大气稳定度ABCDEF表三： 扩散参数（平原地区和城市远郊区）大气稳定度A-BCDE-F表四： 扩散参数（工业区和城市中心区）（二）相关数据的确定通过对官方网站上的统计数据和相关的新闻的整理，可以直接和间接得出如下数据：假设1、放射性浓度c=100毫西弗（对人体会产生危害临界值）2、假设天空云量为5/87/8 3、云层高度（21344877）m4、35日照角605、地面风速为4m/s，6、扩散系数计算公式：=；=7、泄漏源出口温度与环境温度差大于35摄氏度（=30=80）确定的计算公式：8、吸附系数，。预期放射性物质含碘，则。（三）疏散范围的计算方法1、根据客观上的一些基础数据，代入建立的模型当中，计算污染物扩散的地面浓度，即设立Z=0时，平面上的污染物浓度的分布，得到浓度C与（x , y）之间的模型关系。2、以放射性浓度c=100毫西弗（对人体会产生危害临界值）做为临界值，使核泄漏污染物的到达某点的地面浓度不超过该临界值，确定一个二元方程式，并通过matlab编程进行直观图的制作和临界区域的图形描述，确定需要及时疏散人群的区域范围。3、根据2中所得二元方程式，计算疏散区域的面积以及受核污染不同损害程度的区域面积，以对核电发展必要性分析中的数据进行检验。生成图形为：图三：核泄漏模拟图因为数据的准确性有待验证，所以经过初步估计，下风向范围，最大半径球的半径R=7.12公里。上风向上，最小半径球的半径R0=3.52公里。图形中分线长度为：7.12+11.11（球心位移）+3.52=21.75（公里）。即放射性物质最远可达沿风向距离放射源21.75公里处。图形面积约为173.4（平方公里）。并以此确定疏散范围。（四）应对措施对周围居民来说，采取辐射防护措施的依据是辐射监测情报。应根据不同的辐射水平采取不同的防护措施。如果环境辐射水平低于通常的限制水平时，就不必采取防护措施。如果环境辐射水平达到或者超过所规定的警报水平时，就应当迅速通知周围居民关闭门窗，这可以使居民受到的辐射剂量减少到1/5；或者通知居民执行撤离计划。 1、自我防护方法 空气有放射性污染时，屋外人员可利用口罩、手帕、毛巾、纸等捂住口鼻，对呼吸道进行防护，可使吸入的放射性核素所致剂量减少到1/10。 2、服用稳定性碘保护甲状腺 成人一次服用100稳定性碘碘（相当于130KI或170KIO3），一般在530分钟内就可阻止甲状腺对放射性碘的吸收，大约在一周后对碘的吸收恢复正常。3、撤离 从核泄漏现场迅速撤离是避免或减少各种途径照射的最有效的防护对策。 4、永久性重新定居 如果预计在1年或2年之内，月累积剂量不会降低到10 mSv以下，则应考虑实施不再返回原来家园的再定居。5、对食品和水的去污 控制饮食：碘131在释放后48小时就牛奶中出现。受污染的食物，如牛奶、水果、蔬菜、谷物，可采用加工、洗消、去皮等方法除污染，也可在低温下保存，使短寿命的核素自行蜕变。 6、加强辐射防护意识或知识 在核和辐射突发事件的处置过程中，政府和救援人员加强与公众的沟通十分重要。7、清除污染 辐射防护除了保护工作人员和公众及其他们的后代以外，还有一个重要的责任就是保护我们人类和其他生物赖以生存唯一家园地球环境。（五）核泄漏的主要的后果 核泄漏导致的后果与泄漏方式、量和范围以及环境情况、气象条件等诸多因素有关。严重的核泄漏事故，除了直接造成人员伤亡外，还会引起潜伏期较长的癌症为主的恶性疾病和遗传性疾病、以及人们心理恐慌和社会经济秩序混乱。 1、污染我们的生存环境 较大的核泄漏事故造成的环境污染所产生的后果，无论是在持续时间上还是在波及范围上都比它直接照射产生的后果更严重。2、引起有害的辐射生物效应 核泄漏事故发生后控制放射性物质扩散很难，可引起受照者造血功能障碍、神经系统损伤等近期效应；也可诱发可怕的癌症或恶性肿瘤等持久的威胁健康和生命的远期效应；还导致破坏遗传基因、先天性 3、可造成巨大的经济损失 较大的核事故所造成的经济损失巨大。为处理切尔诺贝利核电站事故前苏联政府动用了全苏必要的科学、技术及经济方面的应变力量。4、对公众的社会心理影响更大 几次重大核事故的经验证明，核事故也可造成很大社会心理影响，引起人群心理紊乱、焦虑、恐慌和长期慢性心理应激。这种不良的社会心理效应，其危害可能比辐射本身导致的后果更严重。由于核武器首次直接用于战争，已经留下了可怕的印象和后5、公众对核和辐射的重新审视 核事故和核泄露等事件使人们重新审视核电与核设施的安全问题。一旦发生事故，造成的危害是长期的甚至是全球性的。在美国，核武器试验引起长期的公众恐惧和社会动荡，并对政府进行谴责，包括诉讼。切尔诺贝利核电站事故进一步促使公众对国家4.3.3问题三对于我国未来核电站的前景预测，我们将从两方面考虑。首先，我们假定我国在30年内技术转化率现在水平不变，即我国未来经济形势状态良好。现阶段，我国的技术转化率约为3%，通过（*）式，结合计算机我们得到图像如下：通过稳态图像，我们可以看到，我国的核电站稳态数量呈递减趋势，最终在30年后达到65座的稳态水平。但是现阶段，我国核电需求量远远大于供给量，所以现阶段，需要大力发展核电事业，填补较大的用电缺口。在10年之后的时间段，应该减缓核电的发展速度，提高核电的规模和技术水平，保持较少的核电数量。同时，随着核燃料原料成本的加大，核电站成本的上升，也会导致核电站数量的减少。其次，我们给定30年的考察期，探讨我国对于科技的投资水平对于核电站发展前景的影响。通过（*）式，结合计算机我们已经得到图像如下：可以看出，如果在未来我国一段时间内，我国逐步加大对科技研究投资的比例，加大新能源的开发力度，则核电需求数量会逐渐减少，而逐渐被其他绿色新能源所替代。5、模型评价模型A：1、 模型优点第一：根据实际情况，参考高斯模型的已受到广泛运用并得到很好的认可的模型。并加以考虑了地面反射、烟云抬升高度，、放射性物质自身的沉降作用、雨水的吸附作用以及放射性衰变等因素对模型进行修正，得到了更优化的连续点源高斯扩散模型。此模型能较好的求出核电站周边放射性物质浓度，并与模型一求出的浓度进行比较，从而给出核电站周边放射性物质浓度的变化情况。第二：在高架连续点源核污染扩散模型高斯模型的基础上，尝试加以更全面的预测下风向的污染物扩散情况。继承了模型二的优点，使放射性物质浓度的预测更加精确。第三：我们建立的以上模型在核事故的应急救援过程中，对救援人员划定警戒区和确定周围居民的疏散范围具有重要意义，并可为制定救援方案和应急决策提供科学依据。对放射性云团在空中迁移和扩散提供浓度的定量描述，也对可能发生的核事件的放射性核素浓度监测及监测时间范围提供相关信息。2、模型的缺点由于扩散是一个非常复杂，影响因素众多的过程，本文中只考虑了热力抬升、核衰变等五个方面对模型的修正，难免有