核电站安全性分析与发展预测数学建模——陈宗则、熊桄清、张泽荣15142.docx

《核电站安全性分析与发展预测数学建模——陈宗则、熊桄清、张泽荣15142.docx》由会员分享，可在线阅读，更多相关《核电站安全性分析与发展预测数学建模——陈宗则、熊桄清、张泽荣15142.docx（18页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、核电站安安全性分分析与发发展预测测陈宗则 熊桄清清 张泽泽荣摘要本文意在在解决核核电站的的发展预预测问题题。发展展预测涉涉及到核核电站的的安全性性考核以以及核电电站经济济效益的的研究。为为此本文文一共分分为问题题提出、问问题分析析、模型型假设、模模型建立立与求解解、模型型评价、参参考文献献以及附附录七大大部分。其其中包含含三个模模型的研研究，核核泄漏放放射性物物质扩散散模型模型型I，核核电站边边际效益益评估模模型模型III以及及核电站站发展预预测模型型模型型IIII。模型I是是为了模模拟核泄泄漏的扩扩散范围围与路径径，运用用了烟雾雾扩散模模型，结结论是扩扩散面积积约为1173.57平平方公里里

2、。模型III是研究究核电站站的最大大经济效效益与技技术进步步的关系系，模型型为边际际效益模模型。先先建立经经济收益益Y的函函数，在在通过求求导得出出边际效效益函数数。边际际效益为为0时总总经济收收益最大大化，一一次可得得出经济济层面上上，为使使总经济济收益最最大化的的合理核核电站数数目与技技术进步步的配对对组合解解集。模型IIII在模模型I、模模型III的基础础上从总总经济收收益函数数出发，推推导出综综合考虑虑安全因因素与经经济收益益的发展展预测模模型。由由此得出出使中国国现有77座核电电站承担担起未来来20年年发生核核泄露风风险所要要求的技技术进步步的临界界值。如如若现实实的技术术进步大大于

3、临界界值，则则可以扩扩大建设设核电站站；反之之则不宜宜扩建，但但是核发发电技术术依然值值得发展展。三个模型型综合研研究了核核电站的的安全性性、效益益型以及及发展趋趋势。同同时在将将情景极极可能简简化的前前提下，得得出相对对可信的的预测结结果。具具有较强强的参考考价值。关键词：烟雾扩扩散模型型 边际际效益分分析模型型 量化化分析 控制变变量法一、问题题的提出出20111年3月月11日日，日本本遭受了了9级大大地震并并引发了了强烈的的海啸。这这次大地地震及其其引发的的海啸不不仅给日日本以重重创，而而且由此此造成的的福岛核核电站的的核泄漏漏更是引引起了全全世界对对核电站站及其安安全的重重新思考考。请

4、从互联联网或报报刊上搜搜集有关关数据，根根据这些些数据建建立评估估核电站站安全的的数学模模型。考考虑：（1） 随着人们们生活水水平的提提高，用用电量大大幅增加加，假设设不建设设核电站站，用电电量和发发电量之之间的差差距有多多大？建建设一个个某种规规格的核核电站能能提供多多少电力力？建设设核电站站的经济济成本和和效益如如何？目目前国内内有几个个核电站站、在建建或准备备建设的的有几个个？也就就是要求求建立建建设核电电站必要要性的数数学模型型并分析析；（2） 以秦山或或大亚湾湾核电站站为例（选选一个），如如果这些些地方出出现了严严重的自自然灾害害造成了了核泄漏漏（需要要你自己己作出合合理假设设），那

5、那么，在在一定气气象条件件（一定定风向、风风力、下下雨等）的的情况下下建立核核扩散的的数学模模型，并并讨论对对周围多多大范围围的居民民进行疏疏散以及及其他的的应对措措施和可可能的后后果；综合前两两问给出出核电站站的发展展前景预预测。二、问题题的分析析 自然灾灾害造成成核电站站的核泄泄漏问题题具有不不可预测测性与严严重性。不不可预测测性指的的是自然然灾害发发生的不不可预测测；严重重性指的的是核泄泄漏带来来的居民民生命财财产损失失的危害害严重。不不可预测测性与严严重性使使对于核核电站安安全的评评估问题题转化为为对从事事发之后后补救措措施实施施过程的的有效性性和成本本的评估估问题。而而对于事事后补救

6、救措施研研究的首首要问题题就是对对于核泄泄漏影响响区域的的模拟。为为此我们们建立了了核泄露露放射性性物质泄泄露模型型模型型I。模模型I采采用了烟烟雾扩散散模型，通通过量化化固定风风向风力力和降雨雨等因素素，定量量分析了了发射性性物质可可能的扩扩散范围围。基于于此，我我们定性性将讨论论在污染染范围内内所需采采取的补补救措施施的可行行性，已已经经济济成本。从从而对核核电站安安全性进进行评估估。对于核电电站的发发展预测测问题，我我们分为为三大步步：第一一步是利利用扩散散模型模模拟出核核泄漏可可能污染染的区域域，从而而计算损损失成本本；第二二步是建建立边际际效益模模型，并并且利用用分离变变量法，研研究

7、技术术进步对对于效益益的影响响，从而而得出效效益关于于技术进进步的函函数，由由此我们们建立了了核电站站边际效效益模型型模型型II；第第三步基基于扩散散模型得得出污染染范围以以及边际际效益模模型得出出效益技术进进步函数数，我们们建立了了发展预预测模型型模型型IIII，讨论论在不同同技术进进步的假假设下，核核电站的的发展趋趋势。三、模型型假设一、核泄泄露放射射性物质质泄露模模型模型II1、核泄泄漏看作作在地面面某一点点向四周周等强度度低随时时释放放放射性物物质，放放射性物物质在无无穷空间间扩赛，不不计地形形影响。2、放射射性物质质扩散服服从扩散散定律，即即单位时时间通过过单位法法向面积积的流量量与

8、它的的浓度梯梯度成正正比。3、放射射性物质质在穿过过降雨时时，其强强度由于于雨水的的吸收而而减少，减减少比率率为常数数。4、风的的影响只只考虑放放射原点点的某一一风向、某某一风速速的影响响。风的的影响范范围，随随着时间间而减小小。二、核电电站边际际效益评评估模型型模型型II1、排除除通货膨膨胀影响响因素，以以及假设设现金回回报率为为0。2、以大大亚湾核核电站为为标准核核电站3、由于于大亚湾湾核电站站截止220099年7月月31日日，大亚亚湾核电电站1号号机组实实现安 全运行行25339天，创创造了国国内核电电站单机机组安全全运行最最高纪录录。所以以假定安安全运行行是核电电站的常常规状态态，既不

9、不考虑发发生事故故。4、假设设目前只只有火电电站（化化石燃料料为能源源的发电电站）和和核电站站（新兴兴绿色能能源核电电站）5、假设设投入到到研究核核电技术术的资金金等于技技术进步步带来的的经济收收益三、核电电站发展展预测模模型模型IIII1、满足足模型II、III假设2、核电电站发生生事故的的概率与与全世界界核电站站的发生生泄露概概率频率率相等四、模型型的建立立与求解解核泄露放放射性物物质泄露露模型模型型I一、模型型设计放射性扩扩散的范范围和强强度主要要受到放放射源的的自然扩扩散、风风以及降降雨的影影响。自自然扩散散符合烟烟雾扩散散模型，在在给定放放射源浓浓度等数数据时，可可以求出出静态扩扩散

10、范围围。风力力具有加加大扩散散面积的的作用，而而降雨具具有吸收收稀释放放射性元元素的作作用。二、变量量解释k:扩散散系数Q：放射射性物质质泄漏总总量C：某一一点的放放射性物物质浓度度t:扩散散的时间间r：降雨雨对于放放射性物物质的吸吸收比率率ws：风风速wa：风风向wd：风风力影响响范围每每天减小小量三、模型型建立与与求解1、建立立空间坐坐标轴将核泄漏漏时设为为t0=0，以以核电站站为坐标标原点，建建立以正正东为xx正方向向，正北北为y正正方向，建建立三维维直角坐坐标系。2、推导导放射源源自然辐辐射范围围半径：时刻t无无穷空间间中任意意一点（xx,y,z）的的放射性性物质浓浓度记为为C。根根据

11、假设设2，单单位时间间通过单单位法相相面的流流量为：q=-kkgraadC（11）k是扩散散系数，ggradd表示梯梯度，负负号表示示有浓度度高香浓浓度低的的地方扩扩散。考察空间间域，的体积积为V，包包围的曲面面为S，SS的外法法线向量量为n,则在t,tt+t内内通过的流量量为： Q1=tt+tSqnddt（22）而内放放射性物物质的增增量为：Q2=vC（tt）-CC(t+t)ddV（3）有质量守守恒定律律：Q1=QQ2（44）根据曲面面面积分分的奥氏氏公式：Sqqnd=vdivv q ddV（5）其中diiv是散散度记号号。有（1）（5）不不难得到到： Ct=kdiiv(ggraddC)=k

12、(2C/x2+2C/y2+2C/z2)，t0,- x,y,zz或t时，CC0通过将xx2+y2+z2=R2代入方程程（8）可可以解除除R：R2=(-4)ktln(C/Q)(4kt)3/22（9）R为在给给定扩散散系数kk，影响响时间tt，放射射性物质质浓度CC，以及及放射源源总泄漏漏量Q的的情况下下的受自自然扩散散影响范范围的半半径，即即在不考考虑风力力、降雨雨情况下下的影响响范围。通过maatlaab，可可以实现现该模型型的图像像生成，如如图一。在此假定定：k=0.000001(扩散系系数)Q=100000000（放放射源总总量）C=1000（放放射性物物质浓度度）t=100000000(秒

13、)（扩扩散时间间）（程序见见附录一一，此次次输入参参数为：exppandd(0.000001,100000000,1100,100000000,00.5,0,00,0)）计算机生生成图像像为：图一球面半径径R=77.122公里，即即距离放放射源77.122公里内内均有受受到污染染3、考虑虑风力与与降雨因因素的影影响:降雨因素素：假设降雨雨对于放放射性物物质的吸吸收率为为r,则则方程（99）变为为：R2=(-4)ktln(C/(11-r)Q)(4kt)3/22（10）风力因素素：只考虑产产生在辐辐射原点点某一角角度waa，某一一强度wws的风风力影响响。风力的影影响范围围从R开开始每天天减少ww

14、d，并并且对于于影响范范围内每每一点的的作用相相同。具体计算算机模拟拟实现见见附录一一。综上两个个因素，在在mattlabb上生成成考虑了了风速和和降雨因因素的影影响范围围图像（程序见见附录一一，输入入参数为为：expaand(0.0000001,1100000000,1000,1100000000,0.5,00.044,0.7855,0.3)）变量假设设：k=0.000001(扩散系系数)，参参考了其其他元素素粒子的的扩散系系数均为为此数量量级Q=100000000毫毫西弗（放放射源总总量），发发生七级级核泄漏漏时的放放射性物物质泄漏漏量C=1000毫西西弗（放放射性物物质浓度度），会会对人

15、体体产生危危害的浓浓度t=100000000(秒)（扩扩散时间间）r=0.5(假假设雨水水对于放放射性物物质的吸吸收率为为50%)ws=00.044(公里里每小时时)（假假设作用用在放射射源的风风速为440米/小时）wa=00.7885（假假设风向向与x00y面平平行，且且与x轴轴成正445度角角）wd=00.3（假假设风力力影响范范围从RR开始，每每天减弱弱3000米）计算机生生成图像像为：图二注：图形最大大半径球球的半径径R=7.12公公里，最最小半径径球的半半径R0=3.52公公里，图图形中分分线长度度为：77.122+111.111（球心心位移）+3.552=221.75（公公里）。即

16、即放射性性物质最最远可达达沿风向向距离放放射源221.75公里里处。图图形面积积为：(14.24+7.004)*7/22+pii*7.122/22+*3.522/22=1773.557（平平方公里里）四、模型型结果分分析模型的分分析结果果说明在在给定的的条件下下，扩散散的最远远距离为为21.75公公里。由由于现实实中风向向的不确确定性，出出于安全全考虑，在在核泄漏漏的中心心周边221.775公里里范围的的居民都都应该撤撤离。从从模型可可以看出出，主要要影响范范围的因因素主要要有核泄泄漏的总总量，以以及核泄泄漏在扩扩散过程程中的稀稀释。因因此控制制核泄漏漏主要有有两点可可以考虑虑。一个个是核泄泄

17、漏源，一一个是风风力对于于核泄漏漏影响过过程。从从核泄漏漏源考虑虑，可以以通过注注水，封封堆等控控制方法法减少核核泄漏的的量。比比如建立立一个检检测系统统，一旦旦发生核核泄漏，可可以自动动注水或或者封堆堆，以减减少核物物质的泄泄露总量量。从核核泄漏扩扩散过程程考虑，可可以通过过人工降降水来吸吸收泄露露出的核核物质，或或者在核核电站的的选址时时考虑地地形因素素，尽量量选择背背风点选选址。核电站边边际效益益评估模型型II：一、模型型设计本模型拟拟用核电电站的边边际经济济效用来来衡量核核电站的的存在必必要性。假设核电电站的总总经济效效用主要要受以下下因素影影响：1、建设设核电站站的投资资的产出出与成

18、本本2、用于于技术研研发的资资金3、发电电量供需需缺口4、环保保优势二、变量量解释Y：所有有核电站站总经济济效益产产出D：全国国对于用用供电经经济效用用总需求求缺口 N：现有有核电站站总数目目I：建设设一座核核电站所所需投资资的金额额（元）Qn：一一座核电电站年发发电量（千千瓦时）Qc：全全国每年年的用电电量（千千瓦时）Pe：电电价(元元/每千千瓦时）R：真实实GDPP增长率率r:环境境治理投投资占GGDP的的比重Ep:环环保系数数Cn:核核电站碳碳排放量量Cf:火火电站太太排放量量Td：投投资于单单个核电电站促进进技术进进步的资资金所带带来的经经济效益益t:研究究的时间间范围k：中国国用于研

19、研究开发发领域经经费占国国内生产产总值的的比重三、模型型建立本模型通通过研究究边际分分析，研研究核电电站最优优数量，如如果最有有数量等等于或大大于目前前我国已已有核电电站数目目+在建建核电站站数目+计划建建设核电电站数目目，那么么我们认认为核电电在有继继续建设设的需求求；反之之，则没没有必要要性。模型为：一、本模模型认为为所有核核电站的的总经济济效益分分三部分分影响因因素：1、核电电站在研研究时间间范围内内通过发发电产生生的经济济收益，即即：(QnNPet (1-r) -Epp（11.1）其中QnnN为所所有核电电站总发发电量，PPe为电价价，t为为有效时时间；因为考虑虑到核电电站为绿绿色能源

20、源，解约约了治理理环境的的成本，所所以出去去了治理理环境的的成本。在在这里假假设核电电站用于于治理环环境的资资金占总总产出的的比例与与中国投投入环境境治理的的费用占占GDPP的比重重相同，即即r。与与此同时时，引入入判断核核电站是是否为绿绿色电站站的指标标环保系系数Epp，如果果Cn(Cn+Cf)00.5,则Epp=0，反反之为11。2、建设设时投入入的建设设成本，即即：IN（1.2）其中I为为投资成成本，NN为核电电站的数数量3、由于于技术进进步带来来的经济济效益，即即：TdNN（11.3）其中Tdd为一座座核电站站由于技技术进步步产生的的经济收收益综合（11.1）、（11.2）、（11.3

21、），得得出所有有核电站站的总经经济效益益函数YY：Y= (QnPe(1+r) EptNIN+TTdN（11）二、本模模型技术术进步主主要受国国内对于于发展该该技术的的投资的的影响：1、假设设国内对对于发展展核电站站的需求求是国内内总发电电量与所所有核电电站发电电量得差差值的货货币表示示，即：C=( Qc- QQnN) Pe (11+R)t（2.1）其中，QQc- QQnN为国内内总发电电量与所所有核电电站发电电量得差差值，PPe为电价价，目的的在于将将该差值值用货币币表示；与此同同时，考考虑到对对于技术术进步的的需求会会随着经经济增长长而增加加，所以以在此假假设对于于技术进进步的需需求每年年以

22、真实实GDPP增长率率增长，即即：(11+R)t。2、对于于技术进进步的投投资为对对于技术术进步的的需求的的函数，即即：kC（2.2）在此假定定对于核核电站技技术进步步的投资资额占需需求额的的比重与与我国投投资于研研究开发发领域资资金占GGDP的的比重相相等，即即为k。3、技术术进步带带来的单单位核电电站经济济产出等等于技术术进步的的投资，即即：Td=kkC（2）在此假设设技术是是可以无无成本复复刻的，所所以一旦旦技术进进步了，所所有的核核电站的的经济效效益都同同时提升升相同量量，且等等于投资资额。三、将所所有核电电站的总总经济效效益Y变变形:将（2.1）、（22）代入入（1）中中，得：Y=

23、(QnPe(1-r)-1EptI+ k QcPe (11+R)t)N- kQnPe (11+R)tN2（33.1）令:a=(QQnPe(1-r)-1EptI+ k QcPe (11+R)t);b= kkQnPe (11+R)t则（3.1）变变为：Y=aN-bbN2（33.2）对Y关于于N求导导，得出出所有新新建核电电站边际际经济效效益产出出：dYdN= a-2bN（33.3）有边际效效益理论论知，当当边际效效益为00时，总总效益达达到最大大化。令核电站站的边际际经济效效益产出出为0，即即dYdN=0,得得最大核核电站数数量N*为：N*=a2b（33）其中a=QnnPe(1-r)-1EptI+

24、k QcPe (11+R)t;b= kkQnPe (11+R)t那么出于于（3）上上的每一一点都是是可以使使总效益益最大化化的，或或者说是是边际效效益为00的点。我我们认为为这就是是最合理理的核电电站数目目。通过计算算机实现现模型III在这里，我我们打算算用控制制变量法法，隔离离出技术术进步因因素K对于合合理核电电站数目目的影响响。以k为自自变量，研研究合理理核电厂厂数N：Qn=1150000000000000（千千瓦时）单单座核电电厂年发发电量（数数据来自自中国核核能协会会网关于于大亚湾湾核电站站年上网网电量统统计）Pe=全全国电价价并不唯唯一，但但大多分分布在00.20.88之间，故故在此

25、假假设为00.6元元/千瓦瓦时r=1.49% 为国国家环境境治理投投资占GGDP的的比重（数数据来自自中国环环境保护护部官网网，20008年年环境统统计年报报）Ep=11，由之之前的定定义知，由由于Cnn(Cnn+Cff)=000.5（数数据来自自中广核核集团官官网）I=26646880000000000（元元人民币币），单单座核电电厂总投投资额。折折合美元元40.72亿亿美元（数数据来源源中广核核集团官官网）Qc=4470000000000000000（千千瓦时）全全国年用用电量（数数据来自自中国电电力联合合会官网网中国电电力联合合会）t=为考考察年限限，在此此定200年。R=为真真实GDD

26、P增长长率。真真实GDDP增长长率每年年略有不不同，在在此设为为8%。其函数图图象为：（具体程程序见附附录二，测试试数据为为：nucllearr(1550000000000000,00.6,0.001499,1,26446800000000000,4470000000000000000,220，00.088）)图四注：1（1）平平行于KK轴的是是N=00,直线线一；（2）平平行于NN轴的是是K=0.662（%），直线线二；（3）曲曲线为：N*=a2b ，其其中R为为自变量量。设为为曲线11.其中：aa=QnnPe(1-r)-1EptI+ k QcPe (11+R)t;b= kQnPe (11+

27、R)t结果分析析（2）如图四，在在给定时时间和其其他因素素的基础础上，单单一研究究技术进进步K对于核核电站数数量的影影响，可可见核电电站数量量呈现增增长趋势势，但是是边际增增长递减减。当技技术进步步率在00.622%水平平左右时时，核电电站的合合理数目目为0，也也就是说说在此点点之前的的点建设设核电站站出于亏亏损状态态；此点点之后建建设核电电站可以以盈利。在在技术进进步水平平出于00.622%至100%之间间时，随随着技术术进步水水平的增增加，合合理核电电站数目目不断上上升。当当技术进进步突破破10%点点时，核核发电技技术足够够成熟，核核电站数数目达到到饱和为为1566座左右右。核电站发发展预

28、测测模型模型型IIII一、模型型设计从模型II我们得得出了核核泄漏的的扩散路路径和范范围，即即1733.577平方公公里；从从边际效效益模型型我们得得出N=N（kk）是单单调递增增函数，进进而得到到核电站站最大经经济收益益为技术术进步的的单调增增函数（即即在边际际收益为为0的前前提下）。在此，我们的模型III从边际效益模型开始推导出综合考虑经济效益和安全性的发展预测模型。二、变量量解释Y：所有有核电站站总经济济效益产产出N=7座座（我国国已建成成），在在建9座座，计划划建设22座I=26646880000000000（元元）Qn=115000000000000000（千瓦瓦时）Qc=44700

29、00000000000000（千千瓦时）Pe=00.6元元每千瓦瓦时电价价，电价价R=8%，真实实GDPP增长率率r=1.49%，环境境治理投投资占GGDP的的比重Ep=11，环保保系数t=200年，研研究的时时间范围围k：中国国用于研研究开发发领域经经费占国国内生产产总值的的比重P0：一一座核电电站一年年内核泄泄漏发生生概率P：在研研究时间间范围内内发生核核泄漏的的概率C：核泄泄漏带来来单位面面积的经经济损失失A：搬迁迁费用B：污染染地区在在未被污污染前每每年的GGDPS=1773.557平方方公里，污污染面积积detaal:核核电站产产生的总总经济收收益与核核泄漏带带来的期期望损失失的差值

30、值三、模型型建立与与求解模型为：模型III的（33.1）式式，将tt=200等参数数代入，得得到最大大总经济济收益YY关于技技术进步步k的函函数。Y=-55.7337010122+8.99952210144k（11）全世世界4338座核核电站，550年内内共发生生10起起事故，故故每年发发生的先先验概率率为：P0=110/(4388*500)= 0.000044.56662（2）那么200年故事事发概率率为：P=200*100/(4438*50)=0.00991（33）查得得单位面面积的GGDP损损失加上上迁移费费（假设设20年年无法恢恢复生产产）为：C=A+B201、假设设单位GGDP产产值

31、与深深圳单位位面积GGDP产产值的数数量相等等，即为为4.888077亿元/平方公公里，故故设B=1亿元元/平方方公里=10992、假设设搬迁费费用与三三峡移民民单位面面积搬迁迁费用数数量级相相等，即即5100亿/333899.255万平方方公里或或1555050000元元/平方方公里。故A=11000000000S=11735570000000=1.73557109元C=1773577000000+100000000000020=2.11736610100(元)（44）由此此，我们们可以退退出核电电站的总总经济收收益与核核泄漏期期望损失失的差值值：deltta=YY- PPSC =-5.77

32、370010122+8.99952210144k-0.009911733.5772.11736610100= -55.77713ee+0112+8.99952210114k（5）然后我们们假定差差值为00，即核核电站的的总收益益与核泄泄漏期望望损失相相等，可可以解出出：K*=00.00064（6）K*的含含义是，为为使现有有的7座座核电站站承担起起未来220年核核泄漏安安全事故故的风险险所需要要的核发发电技术术进步的的投入经经费占全全国对于于核电站站发电需需求的比比例。所以：当kKK*时表明明核电站站可以承承担更大大的风险险，所以以可以继继续增加加建设。当kK*时表明核电站承担风险的能力已不足

33、，不宜继续增加建设。五、模型型评价模型I：优点：此模型在在把问题题最大程程度简化化的基础础上，解解决了放放射性元元素扩散散的影响响范围模模拟。有有助于对对情况有有直观的的认识。对对于最快快反应出出应对措措施有着着借鉴意意义。缺点：此模型将将风速、风风向、风风力以及及降雨都都假定为为定值或或者简单单的线性性函数，与与现实情情况有一一定差距距。如果果能够这这些因素素改为更更贴近现现实的函函数，将将有助于于模型的的推广。模型III：优点：此模型充充分考虑虑了时间间，发电电量，电电力供应应缺口，科科技进步步以及核核电站的的环保性性等因素素，起到到了确立立技术进进步与基基于经济济层面合合理核电电站的函函

34、数关系系。缺点：对于技术术进步、环环保性的的量化不不够精准准。如果果可以有有更加贴贴切的函函数模拟拟，可以以使模型型结果更更为准确确。模型IIII：优点：模型综合合考虑了了核电站站的安全全性能，以以及不同同技术进进步带来来的不同同总收益益，继而而从风险险综合收收益的角角度对在在不同技技术进步步前提下下核发电电产业未未来可能能的走向向。缺点：对于迁移移费用、单单位面积积损失不不够具体体落实到到不同核核电站的的具体经经济状况况。六、参考考文献1 王玉磊磊 邱罡罡 编著著 从从零开始始学MAATLAAB 中国铁铁道出版版社2 大亚湾湾核电运运营管理理有限责责任公司司官网 n/n2290336/iin

35、deex.hhtmll3中中广核官官网 wwww.cggnpcc.coom.ccn/nn288819559/iindeex.hhtmll4国国家电力力信息网网 htttp:/n/5中中国电力力企业联联合会官官网 hhttpp:/tj.cecc.orrg.ccn/6中中国核能能协会官官网 hhttpp:/wwww.chhinaa-neea.ccn/7中中华人民民共和国国统计局局官网 htttp:/8中中华人民民共和国国环境保保护部官官网 hhttpp:/七、附录录附录一：核泄露放放射性物物质泄露露模型模型型I计算算机实现现:funcctioon =eexpaand(k,QQ,C,t,rr,wss

36、,waa,wdd)R=sqqrt(-44)*kk*t*logg(CC/(1-rr)*QQ)*(4*pi*k*tt)(1.55)for i=00:244:t/36000if (wd*(i/24)=R R=RR-wdd; R0=R;t=0:0.11:0.5*ppi;p=0:0.11:2*pi;wsx=ws*coss(waa)*ii;wsy=ws*sinn(waa)*ii;theeta,phii=mmeshhgriid(tt,p);x=R*sinn(thhetaa).*coss(phhi)+wsxx;y=R*sinn(thhetaa).*sinn(phhi)+wsyy;z=R*coss(thhetaa);holdd onnsurffc(xx,y,z)endendR0=RR0vieww(200,700)titlle(核泄漏漏有风向向影响);xlabbel(x);ylabbel(y);zlabbel(z);axiss(-10,15,-100,255,0,10)x0=wws*ccos(wa)*100000000/36000y0=wws*ssin(wa)*100000000/36000R1=ssqrtt(xx0)2+(y0)2)end附录二核电站关关于技术术进步的的边际效效益模型型模型型II：