概率安全评价(Probabilistic Safety Assessment, PSA)的不确定性分析和敏感性分析是研究非能动系统可靠性的有效方法[1-2]，非能动系统敏感性分析是评价系统输入参数的不确定性对系统输出不确定性的贡献率[3-4]，实际上是分析特定输入参数及其初始值发生变化对系统物理过程成功准则的影响[5-6]，其结果具有重要的参考价值，是工程决策的重要依据，也是可行性研究报告和工程评估报告不可缺少的组成部分[7-8]。倪超[9]采用Spearman秩相关系数作为敏感性度量手段，识别AP1000核电厂大破口失水事故下包壳峰值温度(Peak Cladding Temperature, PCT)的敏感参数；张中伟[10]和黄昌蕃[11]采用标准化回归系数表征各参数的敏感性影响，识别系统成功准则的敏感参数；文献[12-13]采用Pearson相关系数度量参数的敏感性，识别成功准则的敏感参数。

由于大型设备使用场景的专业性及特殊性，大型设备的管理都处于分散状态，设备基本信息及使用情况透明度都不高。同时，各单位作为资产占有使用者，管理站位高度不够，仅为满足了本单位事业发展需求占有使用，各自为政，导致整体管理水平低下，是加强大型设备使用绩效管理的“拦路虎”。

综上所述，目前刻画参数敏感性所采用的相关系数(Spearman秩相关系数、标准化回归系数、Pearson相关系数)，仅反映单个输入参数与系统物理过程成功准则的线性关系。而系统不是单输入参数，并且成功准则不止与一个输入参数具有线性关系，可能同时与多个参数存在相关性，即不只是单个参数影响成功准则而是多个参数同时制约其变化。成功准则与每个敏感参数既表现出相关联又有不确定的“特殊”关系，是多个输入参数共同作用的结果。至于这种“共同作用”如何协同显现，每个输入参数扮演什么角色、发挥多大作用，由输入参数组成的系统决定，是系统的“内在”规律。由于相关系数没有把输入参数置于多参数“共同作用”的重要前提下考虑，所以只能反映单个参数与成功准则的“外在”表征，无法表现输入参数在多参数“共同作用”前提下与成功准则的新相关程度。

针对相关系数在刻画参数敏感性方面的缺陷，本文基于偏相关系数提出迭代法，其综合三方面信息：输入参数与成功准则的相关系数、偏相关系数以及输入参数之间的相关系数，可全面准确地分析输入参数的敏感性，科学判定其与成功准则的相关性，合理解释敏感参数的判定结果。迭代法更为精确地判定参数敏感性，对降低参数识别的维数，提高识别效率具有重要意义。

1 基于相关系数识别PCT的敏感参数

本文以AP1000丧失正常给水事故为研究对象，设定PCT为成功准则[14-15]。分析非能动余热排出系统(Passive Residual Heat Removal System, PRHRS)物理过程的参数敏感性，主要研究输入参数对PCT准则的线性相关性，其分析过程如图1所示，重点研究内容集中在第六步敏感参数的识别。

  

图1 PRHRS物理过程成功准则的参数敏感性分析示意图Fig.1 Parameter sensitivity analysis diagram of PRHRS physical phenomena success criteria

式中，r是对应的偏相关系数，n为工况数，k为输入参数个数。

① 相关⟸P>P阈值>P⟹不相关

②

使用Pearson相关系数(简称相关系数)描述参数与PCT的线性关系，计算公式为：

(1)

式中，n为工况总数，xi为第i个工况的输入参数值，yi为第i个工况的PCT。Pearson相关系数又称线性相关系数、积差相关系数，适用于满足正态分布的数据，取值范围在(-1,1)之间，绝对值越大相关性越强，取值符号表示相关性方向。

根据文献[9，11,16-17]所采用的AP1000重要输入参数及不确定性分布，经与安全分析报告的参数设计值对比并征求多位专家意见，最终形成本文采用的输入参数及不确定性分布，见表1。

所谓制度追求价值的考量，完善制度最重要的是把社会主义的基本价值取向充分地反映在制度安排之中。马克思主义经典作家对社会主义社会的价值取向有过系统而丰富的论述，这些价值追求的实现必须通过制度安排才能做到。中国特色社会主义的理论和实践探索告诉我们，建设中国特色社会主义最重要的价值诉求是两条，一是人民当家作主，坚持人民主体地位，一是维护社会公平正义。只有坚持这两条价值原则，才能保证人民共享发展成果，走共同富裕道路，才能促进社会和谐，促进人的全面发展。党的领导、解放和发展社会生产力、改革开放等都是实现如上价值诉求的条件保证和路径选择。

[15] 林诚格. 非能动安全先进核电厂AP1000 [M]. 北京： 原子能出版社， 2008.

 

表1 输入参数及不确定性分布

 

Tab.1 Input parameters and their uncertainty distribution

  

序号参数单位分布1CMT2#隔离阀能阻系数正态(0.44,0.56)2HX水平管径m正态(0.01575,0.002)3二次侧(给水)初始温度K均匀(494.85,504.85)4HX有效总流通面积m2均匀(0.124168,0.144168)5反应堆总功率MWt均匀(3366,3434)6稳压器压力MPa均匀(15.19,15.81)7CMT水温K均匀(300,344.1)8IRWST温度-N①K均匀(300,344.1)9IRWST温度-F②K均匀(306,350.1)10HX水平管长度m均匀(3.0,3.4)11HX进口能阻系数均匀(0.24,0.3)12HX出口能阻系数均匀(0.24,0.3)13一次侧初始温度K均匀(569.87,578.21)14HX垂直管径m均匀(0.018,0.020)15HX热交换面积-外径水平段m2根据管数、长度、外径得到16HX热交换面积-外径垂直段m2根据管数、长度、外径得到17HX热交换面积-内径水平段m2根据管数、长度、内径得到18HX热交换面积-内径垂直段m2根据管数、长度、内径得到19轴向功率分布PBOT均匀(0.9,1.1)-变化因子20轴向功率分布PMID均匀(0.9,1.1)-变化因子21CMT管线能阻系数均匀(0.8,1.2)-变化因子22燃料导热率W/(m·K)正态(1,1.1)-变化因子23CMT1#隔离阀能阻系数正态(0.88,1.12)-变化因子24衰变热MWt正态(0.88,1.12)-变化因子25SG回路能阻系数均匀(0.9,1.1)-变化因子26HX水平管总长度m根据管数、单管长度得到27HX管数个均匀(620,750)

注：①IRWST温度-N：HX附近区域IRWST的温度；②IRWST温度-F：HX较远区域IRWST的温度。

 

表2 传统法判定与PCT存在相关性的参数

 

Tab.2 Parameters having correlation with PCT judging by traditional method

  

参 数相关系数IRWST温度-N0.116HX出口能阻系数0.132HX垂直管径-0.237HX热交换面积-外径垂直段-0.217HX热交换面积-内径垂直段-0.217轴向分布功率PBOT0.288轴向分布功率PMID0.215HX水平管总长度-0.105

PCT参数敏感性分析的目的是通过输入参数的敏感性系数值，删除与PCT弱相关的参数，重点考虑强相关的参数。如何根据敏感性系数确定哪些输入参数是强相关的，哪些又是弱相关的，选用什么标准或指标衡量敏感性的显著性，仍是有待深入研究的问题。目前此方面的研究较少[19]，基本以经验规则或人为设定相关系数阈值为标准，判定输入参数与PCT的相关性，传统方法过于简单粗略，有待细化、精准化。

2 基于偏相关系数识别PCT的敏感参数

2.1 迭代法

相关系数刻画参数敏感性时存在局限性，实际中不能仅根据相关系数盲目判定输入参数与PCT是否具有线性关系，必须考虑第三个参数或更多参数对所关注参数的影响。

偏相关系数描述在控制一个或几个输入参数的条件下，所关注的输入参数与物理过程成功准则之间的相关程度[20]，反映了输入参数与成功准则之间的净相关性。假设控制输入参数Z，则参数X与成功准则Y之间的偏相关系数为[21-22]：

 

(2)

[7] Hirschberg S. Human reliability analysis in probabilistic safety assessment for nuclear power plants[J].Nuclear Safety, 2004.

 

(3)

传统方法分析PCT输入参数敏感性，需根据相关系数判断参数是否与PCT存在线性相关性，唯一的判定依据是相关系数值，具体实施方式有以下两种。

虽然相关系数和偏相关系数均可以表示输入参数与成功准则的线性关系，但其物理意义有差异，如图2所示。

  

图2 相关系数与偏相关系数物理意义的差异Fig.2 Physical difference of correlation coefficient and partial correlation coefficient

图2中(a+b)/(a+b+c+d)表示由IRWST温度-N解释的PCT变异量的比例。若令ry1表示IRWST温度-N与PCT的相关系数，则另一方面，a是在IRWST温度-N和PCT的重叠面积中扣除HX水平管总长度重叠部分后的剩余面积，a/(a+b+c+d)是a+b扣除HX水平管总长度可以解释的部分后，单独由IRWST温度-N解释的PCT变异量的比例。若令sry1表示IRWST温度-N和PCT之间的偏相关系数，则

分析IRWST温度-N的偏相关性，结果见表3。

 

表3 IRWST温度-N的偏相关性分析

 

Tab.3 Partial correlation analysis of IRWST temperature-N

  

参数相关系数偏相关系数HX水平管总长度的相关系数IRWST温度-N0.1160.0310.100

[3] Saltelli A, Tarantola S, Campolongo F, et al. Sensitivity analysis in practice a guide to assessing scientific models [M]. USA: John Wiley& Sons 1nc, 2004.

本设计是一种新型电动车防盗装置，以STC12C5A60S2单片机作为控制器，利用震动传感器采集外界信号，并将信号反馈给单片机，单片机依据接收到的有效信号，控制蜂鸣器与发光二极管的开启，从而实现声光报警；利用红外遥控可对防盗装置进行布控或解控，当开启布控时，装置处于正常工作状态，解除布控则使装置处于休眠状态；此外，设计巧妙地将GPS与GSM功能融入到系统中，可以通过短信通讯实现对车辆的实时定位。系统结构设计图如图1所示。

综上所述，全面准确地分析输入参数的敏感性、科学判定与成功准则的相关性，必须综合三方面信息：输入参数与成功准则的相关系数；输入参数与成功准则的偏相关系数；输入参数之间的相关系数。

综合以上三方面信息提出迭代法解决基于相关系数识别敏感参数存在的问题，方法示意图如图3所示。

  

图3 迭代法示意图Fig.3 Schematic diagram of iterative method

2.2 迭代法识别PCT的敏感参数

设定阈值为0.1，判断参数之间的相关性，分别得到每一个参数与其余参数之间的相关性情况可知，HX水平管总长度与轴向分布功率PBOT存在线性关系，控制PBOT效应的同时分析HX水平管总长度与PCT之间的偏相关性，结果如表4所示。

 

表4 HX水平管总长度与PCT之间的偏相关性

 

Tab.4 Partial correlation between HX total length of horizontal pipe and PCT

  

控制参数参数PCT显著性无HX水平管总长度-0.1050.072PBOTHX水平管总长度-0.0160.761

观察表4可知，不控制轴向分布功率PBOT的效应时，HX水平管总长度与PCT之间的相关性呈负相关，相关系数为-0.105；但PBOT的效应后，HX水平管总长度的偏相关系数为-0.016，与PCT不存在相关性(p=0.761>0.05)。可推知尽管HX水平管总长度和PCT表面上负相关，但是由于它们的相关系数未能超过整个系统带来的阀值——PBOT与PCT的相关系数、两个参数间的相关系数，或者说，HX水平管总长度与PCT的相关程度不足以克服PBOT的相关影响，因而使表面的负相关在系统内部实际上表现为不相关的情形，这就是HX水平管总长度与PCT内在的本质联系。只有偏相关系数才能真正反映输入参数与成功准则的本质联系，所以不能简单地判断HX水平管总长度与PCT之间是否存在线性关系，结论应该为在控制PBOT不变的情况下，HX水平管总长度与PCT之间不存在显著的线性关系。

按示意图3应用迭代法识别与PCT均线性关系的参数，迭代结束后判定的与PCT存在线性关系的参数见表5。

 

表5 迭代结束后与PCT存在相关性的参数

 

Tab.5 Parameters existing correlation with PCT after iteration

  

编号参数相关系数偏相关系数1二次侧初始温度-0.095-0.1012反应堆总功率-0.046-0.0923HX垂直管径-0.237-0.1554HX热交换面积-外径垂直段-0.217-0.1585HX热交换面积-内径垂直段-0.217-0.1586轴向分布功率PBOT0.2880.3007轴向分布功率PMID0.2150.6158CMT1#隔离阀能阻系数0.0310.1479SG回路能阻系数0.0670.120

观察表5发现轴向分布功率PMID迭代结束前后的相关性指标相差较大，原因在于：HX垂直管径与PMID表现相关性，可能是PMID的压抑参数，压制了PMID与PCT的相关性；扣除HX垂直管径的影响后，分析PMID与PCT的偏相关性，相关性指标变大且相关性增强，可推知HX垂直管径的确是PMID的压抑参数。分析输入参数与成功准则的相关性时，不仅要考虑参数的偏相关系数，还应综合参数与成功准则及参数间的相关系数，只有综合三方面信息才能正确识别输入参数的敏感性。

就PCT而言，轴向分布功率PBOT、轴向分布功率PMID、CMT1#隔离阀能阻系数、SG回路能阻系数的偏相关系数为正，参数值越小，PCT值越低，PRHRS物理过程的可靠性越高；其余参数的偏相关系数为负，参数越大，PCT值越低，物理过程的可靠性越高。所有参数中只有轴向分布功率PBOT、轴向分布功率PMID等参数的偏相关系数绝对值大于0.2，其余参数均小于0.2，PMID的绝对值最大，为-0.615。

3 两种方法的对比

3.1 识别的参数结果对比

对比表2和表5可知，迭代法分析9个参数与PCT具有线性关系，比传统法多1个参数。迭代法相比传统法增加的参数为：二次侧初始温度、反应堆总功率、CMT1#隔离阀能阻系数和SG回路能阻系数；删除的参数有：IRWST温度-N、HX出口能阻系数、HX水平管总长度。迭代法增加或删除参数原因：其以参数之间的相关性为基础，分析参数的偏相关性，即控制附加参数效应的同时计算参数的偏相关系数，此系数能真正反映参数与PCT之间的线性关系；传统法独立计算单个参数对PCT的相关系数，而实际应用中往往因为第三个参数或更多参数的影响，使得相关系数不能真正反映指定参数与PCT之间的线性相关程度，所以传统法误认为与PCT存在相关性的参数，由迭代法将其删除，而传统法不能识别的与PCT存在相关性的参数，又由迭代法将其纳入到与PCT存在相关性的参数行列。

3.2 识别的参数效果对比

对比传统法和迭代法分析的与PCT存在相关性的参数效果，对比方案为：全参数样本中分别仅留下传统法和迭代法判定的与PCT存在相关性的参数，将两种方法对应的参数样本分别作为输入参数组合再次代入RELAP5程序，得到相应工况的PCT。对比指标为平均距离，其计算公式为：

 

(4)

式中，n为工况数，yi为传统方法或迭代法第i个工况下的PCT，y全i为全参数第i个工况下的PCT。

张兰花用一颗赤诚的心，连续五年参与了师团“5000致富带动工程”，与女职工郭艳英结成帮扶对子，使其成为团场“脱贫致富标兵”。在她的带动下，连队25名职工成为团场科技致富的“领头雁”。她通过自学获得了第二师颁发的“高级农艺工”证书。

两种方法的PCT代入式(4)计算平均距离，得到传统法距离为1.096，大于迭代法距离0.090，这说明传统法删除的参数中存在对PCT有影响的参数，以及将对PCT影响微弱或无影响的部分参数纳入与PCT存在相关性的行列，使得PCT与全参数偏差较大。相比传统法，这迭代法分析参数的偏相关性，即控制附加参数效应的同时计算参数的偏相关系数，能真正反映参数与PCT之间的线性关系，得到的平均距离小于传统法，这说明迭代法相比传统法可以更准确地识别与PCT存在相关性的参数，保留对PCT影响较大的参数，删除影响较小的参数，所以应采用迭代法分析参数偏相关性，判定与PCT的线性关系。

4 结论

本文分析核反应堆PRHRS物理过程的敏感性，研究输入参数对PCT准则的线性相关性，重点是敏感参数的识别。传统法采用相关系数刻画参数敏感性，只能反映单个参数与成功准则的“外在”表征，无法表现输入参数在多参数“共同作用”前提下与成功准则的新相关程度。针对相关系数的缺陷，本文基于偏相关系数提出迭代法分析输入参数对成功准则的敏感性。迭代法综合了输入参数与成功准则的相关系数、偏相关系数以及参数之间的相关系数三方面信息，在PCT准则下验证了迭代法的适用性，不仅可以全面分析输入参数的敏感性，还可以完整深入地解释敏感参数的判定结果。迭代法适合分析成功准则的参数敏感性，得到精确的敏感参数，对降低参数识别的维数，提高识别效率具有重要意义，是工程决策的重要依据，为将来提高核反应堆非能动系统物理过程可靠性提供可信依据。

参考文献(References)

[1] 扈本学， 王喆， 王国栋， 等. 基于抽样的不确定性及敏感性分析的方法在核电厂水膜蒸发试验误差分析中的应用[J]. 核安全， 2016， 15(1)： 84-94.

HU Benxue, WANG Zhe, WANG Guodong, et al. Application of the sampling-based uncertainty and sensitivity analysis method in error analysis of heated flat plate test at NPP [J]. Nuclear Safety, 2016, 15(1): 84-94. (in Chinese)

[2] 丁虎. 研发核电站PSA软件中不确定性和敏感性分析功能模块[D]. 合肥： 合肥工业大学， 2013： 1-2.

硬皮病以患者出现皮肤炎症、变性、纤维化继而硬化、萎缩等皮肤损害为标志性表现，能并发多器官功能障碍，属于风湿免疫系统疾病中较为顽固凶险的一种。目前国内学者大多在祖国医学“皮痹”的范畴对硬皮病的个别证型及治疗做过初步探讨，但未就硬皮病的中医证治规律独立进行系统性分析。笔者汇总并分析近年来关于硬皮病中医治疗方面的报道，试图整理出硬皮病的中医证治规律的初步思路。

DING Hu. Development of the uncertainty and the sensitivity analysis function module in the nuclear power station PSA software [D]. Hefei: Hefei University of Technology, 2013: 1-2. (in Chinese)

与IRWST温度-N存在相关性的参数不单有HX水平管总长度，但表3只分析与HX水平管总长度的相关性，原因：根据图2可知，只有PCT的敏感参数才会影响IRWST温度-N与PCT的重叠面积(相关性)，而与PCT不具有线性关系的参数不会影响IRWST温度-N与PCT的相关性，所以显著性检验法判定与PCT存在相关性符合要求的参数只有HX水平管总长度。由表3 可知，IRWST温度-N与PCT的相关系数为0.116，有较大负相关性，分析偏相关性可得偏相关系数为0.031，至此能够推断IRWST温度-N与PCT不存在相关性。相关性指标值出现巨变的原因：HX水平管总长度与IRWST温度-N具有较大重叠面积，相关系数为0.100；IRWST温度-N 与PCT的相关性受到HX水平管总长度的影响，当扣除IRWST温度-N与HX水平管总长度的重叠部分后，相关性指标值变小，与PCT的相关性减弱，由相关变为不相关。

[4] Helton J C, Davis F J. Illustration of sampling based methods for uncertainty [J]. Risk Analysis, 2002, 22(3): 591-622.

1995年，张培林在全国首创了分配“模糊弹性”理论。其核心内容为：从分配平等到机会平等；规则透明，数量同级间模糊；二八原则，向骨干倾斜；实行3～5倍级差。该机制既克服“大锅饭”计划经济问题，又有效驱动医院骨干的积极性、保持非骨干的稳定性。该项改革获得2009年度全国五一劳动奖状。

[5] Lofaro R J. Understanding aging in containment cooling systems[J]. Nuclear Engineering and Design, 1996, 163(3): 301-314.

[6] Borgonovo E, Apostolakis G E, Tarantola S, et al. Comparison of global sensitivity analysis techniques and importance measures in PSA [J]. Reliability Engineering and System Safety, 2003, 79(2): 175-185.

式中，rXY,Z表示固定输入参数Z时，参数X与成功准则Y之间的偏相关系数，rXZ、rYZ分别表示参数X与参数Z以及参数Z与成功准则之间的相关系数。对偏相关系数进行显著性检验，检验统计量为：

[8] 陈炼， 胡啸， 邓程程， 等. 最佳估算方法中不确定度评估关键问题分析[J]. 原子能科学技术， 2016， 50(5)： 851-858.

CHEN Lian, HU Xiao, DENG Chengcheng, et al. Analysis on key issue of uncertainty evaluation in best estimate method [J]. Atomic Energy Science and Technology, 2016, 50(5): 851-858. (in Chinese)

HUANG Changfan. Preliminary investigation on reliability assessment of passive safety system [D]. Shanghai: Shanghai Jiao Tong University, 2012: 12-13. (in Chinese)

NI Chao. Modeling of AP1000 nuclear power plant LB-LOCA best estimate analysis and uncertainty study [D]. Shanghai: Shanghai Jiao Tong University, 2011: 47-80. (in Chinese)

推荐理由:今年正值知青上山下乡50周年，留守知青年纪大的已年过古稀，年纪小的也将近花甲，作者从采访过的数百位对象中选出19位最有代表性的知青，这19条命运的曲线构成了16个真实而感人的故事，述说着从下乡到坚守的过往岁月。北大荒留守知青用50年的生命给我们提供了一个人生的坐标，让我们找到自己的位置，让我们懂得什么叫责任，什么叫担当，什么叫承诺，什么叫人生……

[10] 张中伟. 保守分析与最佳估算相结合之LOCA认证分析方法论[D]. 上海： 上海交通大学， 2011： 46-65.

供应不足的大环境，为整个尿素行情提供了长久而稳固利好支撑。随时间的推移，环保力度不断升级，尿素企业的生产阻力不断增大，开工率也难有起色。一些前期检修的企业，仍旧没有复产的消息，而山东、河北、内蒙古等地区部分企业已经出现停产、减产的情况。

ZHANG Zhongwei. Deterministic-realistic hybrid methodology (DRHM) for LOCA licensing analysis [D]. Shanghai: Shanghai Jiao Tong University, 2011: 46-65. (in Chinese)

[11] 黄昌蕃. 非能动安全系统可靠性评估方法初步研究[D]. 上海： 上海交通大学， 2012： 12-13.

[9] 倪超. AP1000核电厂大破口失水事故最佳估算分析建模与不确定性研究[D]. 上海： 上海交通大学， 2011： 47-80.

[12] 王成龙， HU Linwen， 秋穗正， 等. 可移动式熔盐冷却高温堆(TFHR)热工水力特性不确定性分析[C]//第十四届热工流体会议论文集， 2015.

WANG Chenglong, HU Linwen, QIU Huizheng, et al. Uncertainty analysis of transportable fluoride-salt-cooled high temperature reactor (TFHR) [C]//Proceeding of the 14th Thermal Engineering Fluid Conference, 2015. (in Chinese)

[13] Gertman A E, Mesina G L. Uncertainty analysis of RELAP 5-3D: INL/EXT-12-26182[R]. Idaho National Laboratory, 2012.

[14] 孙汉虹. 第三代核电技术AP1000[M].2版. 北京： 国家电力出版社， 2016： 34-78.

试验动物ICR小鼠(雄性，体重18～22 g)和豚鼠、绵羊血红细胞(SRBC)和S180肿瘤细胞均购于医学部实验动物中心；鸡血红细胞购于中国农业大学动物医学院；医药级香菇多糖购于山西泰盛制药有限公司；胎牛血清购于杭州四季青公司；其他常用化学试剂购于北京化学试剂公司(分析纯)。

SUN Hanhong. The third generation nuclear technology AP1000[M].2nd ed. Beijing: China Electric Power Press, 2016: 34-78. (in Chinese)

利用拉丁超立方抽样方法[18]，在输入参数取值范围内抽取188组参数样本，得到稳态或瞬态模型不确定性传递计算的输入参数组合。其中每一组由27个输入参数组成，对应一组RELAP5程序输入。参数样本带入RELAP5程序运行，得到对应工况下的PCT，参数组合与对应的PCT代入式(1)，得到输入参数的相关系数。应用传统法判定参数与PCT是否存在相关性，设定相关系数阈值为0.1，最终判定9个参数与PCT具有线性关系，其结果如表2所示。

LIN Chengge. Passive safety advanced nuclear plant AP1000 [M]. Beijing: Atomic Energy Press, 2008. (in Chinese)

[16] 张顺香. AP1000电厂状态参数不确定性对LBLOCA影响的量化分析[D]. 上海： 上海交通大学， 2013： 41-51.

ZHANG Shunxiang. Uncertainty analysis of the effect of plant state parameter on LB-LOCA in AP1000 [D]. Shanghai: Shanghai Jiao Tong University, 2013： 41-51. (in Chinese)

[17] 刘强. 基于神经网络方法的AP1000 非能动系统可靠性分析[D]. 北京： 清华大学， 2014： 48-61.

LIU Qiang. Reliability analysis of AP1000 passive system based on artificial neural networks [D]. Beijing: Tsinghua University, 2014: 48-61. (in Chinese)

“毫无疑问，对焦在这里是个难点，况且我用的还是一款为胶片机设计的古董镜头。但好在最后拍摄效果不错，这让我不得不对这台相机刮目相看。”

[18] Imam R L, Conover W J. Small sample sensitivity analysis techniques for computer models with an application to risk assessment [J]. Communications in Statistics-Theory and Methods, 1980, 9(17): 1749-1842.

[19] 蔡毅, 邢岩, 胡丹. 敏感性分析综述 [J]. 北京师范大学学报(自然科学版)， 2008， 44(1)： 9-16.

CAI Yi, XING Yan, HU Dan. On sensitivity analysis [J]. Journal of Beijing Normal University(Natural Science), 2008, 44(1): 9-16. (in Chinese)

[20] 王黎明， 陈颖， 杨楠. 应用回归分析[M]. 上海： 复旦大学出版社， 2008： 65-66.

WANG Liming, CHEN Ying, YANG Nan. Applied regression analysis [M]. Shanghai: Fudan University Press, 2008: 65-66. (in Chinese)

工艺过程：主反应锅中加入A相搅拌，加热至70～75 ℃，加入混合好的B相，搅拌溶解均匀后加入C相，搅拌溶解。冷却至40～45 ℃时加入D相，混合均匀，继续冷却至30～35 ℃，加入E相。搅拌均匀，至25～30 ℃时加入F相、G相，混合均匀。

[21] 罗伯特S.平狄克， 丹尼尔L.鲁宾费尔德.计量经济模型与经济预测[M].4版.钱小军, 译. 北京：机械工业出版社， 2005： 57-65.

Pindyck R S, Rubinfeld D L. Econometric Models and Economic Forecasts [M].4th ed. Translated by QIAN Xiaojun. Beijing: China Machine Press, 2005: 57-65. (in Chinese)

种箱的结构示意图如图4所示。排种器壁上开有滑槽，侧壁可以用手工进行转动固定，取种盘是依附在侧壁的内侧，距离侧壁内面尽可能保持在5mm左右，以减小与侧壁的摩擦。

[22] 何晓群. 多元统计分析 [M]. 2版. 北京： 中国人民大学出版社， 2010： 12-15.

HE Xiaoqun. Multivariate statistical analysis [M].2nd ed. Beijing: China Renmin University Press, 2010: 12-15. (in Chinese)

2）摘袋方法一定要正确。首先避免高温、阳光强烈时段摘袋。其次摘袋不能一次性完成。单层袋先撕开底部呈伞状遮挡住果实，过5天左右再全部去除。双层袋先去外袋，过3～5个晴天再去内袋。