基于分离变量的地震PSA相关性分析方法研究.docx

? ? 基于分离变量的地震PSA相关性分析方法研究 ? ? 冯丙辰，杨建峰，张晓明，王 照 (苏州热工研究院有限公司，广东 深圳 518000) 福岛事故后，核电厂外部灾害风险分析越来越受到业界的重视[1-2]。地震概率安全评价(PSA)作为一种定量评价核电厂地震风险的分析方法，已逐步在核电厂的设计和运行中得到应用[3-4]。随着不断地研究和实践，地震PSA方法已逐渐发展成熟，国内外相关技术导则[5]也陆续发布，用于指导地震PSA的开发和应用。 与内部事件PSA相比，地震PSA有两个显著特点[6]：1) 构筑物、系统和部件(SSC)的失效概率与地震的级别有关，需考虑核电厂所有可能发生的地震动范围及其影响后果；2) 地震可能同时对多个冗余或相似的SSC造成破坏，地震PSA定量化过程中需合理考虑这种相关性的影响。对于第1个特点，通过开发体现SSC不同地震动下条件失效概率的易损度曲线和体现厂址不同等级地震超越频率的危险性曲线，对两者进行卷积分[7]计算实现。但对于第2个特点，目前仍缺乏一套完善的处理方法。现有地震PSA中SSC相关性处理一般采用简化的假设[8]：位于相同位置(如相同厂房、高度、系统)相互冗余的同类设备相关性为1，其余相关性为0。该方法将系统冗余列中的设备按完全相关处理，通常情况下被认为是保守的。然而，对于可能存在部分相关性的情况，如不同系统中的相同设备，该方法认为其完全不相关，由此得到的结果将不保守甚至错误。 为更好地体现地震情况下SSC相关性的影响，本文对基于分离变量的易损度相关性分析方法开展研究，并对其在地震PSA中的应用进行讨论。 1 地震易损度评价 1.1 易损度模型 SSC的地震易损度定义为在确定的失效模式和给定地震动参数值(通常为峰值地面加速度PGA)下的条件失效概率。地震易损度评价的目的是确定SSC的抗震能力。抗震能力是一个随机变量，由概率分布来描述[9]： A=AmeReU (1) 式中：A为SSC的抗震能力；Am为抗震能力中值；eR和eU是中值为1的随机变量，分别代表随机性和不确定性。易损度采用双对数正态分布模型[10]，即假设eR和eU分别服从对数标准差为βR和βU的对数正态分布。 在给定失效模式和Am的情况下(仅考虑随机性βR)，对于给定的加速度a，SSC的条件失效概率f0可表示为： f0=Φ[ln(a/Am)/βR] (2) 式中，Φ[·]为标准正态累积分布函数。 当包含不确定性βU时，易损度成为一不确定的随机变量。对于给定的加速度a，易损度f可用主观概率密度函数表示。不超越易损度f′的主观概率Q(也称为置信度)和f′的关系如下： (3) 式中：Q=P[ff′|a]，表示在给定的a下f小于f′的主观概率；Φ-1[·]为标准正态累积分布函数的反函数。a[11]可表示为： p A=AmeΦ-1(f′)βRe-Φ-1(Q)βU (4) 1.2 易损度参数估计 地震易损度评价要确定SSC的易损度参数Am、βR和βU，从而开发易损度曲线。计算易损度参数一般用中间随机变量安全因子F[12]，SSC的抗震能力可表示为： A=FASSE (5) 式中，ASSE为设计要求的安全停堆地震地面加速度能力。 设备的安全因子F可表示为： F=FCFSRFRE (6) 式中：FC为能力因子；FSR为结构响应因子；FRE为设备响应因子。 影响结构响应变化的变量包括谱形(SA)、地震动空间变化(GMI)、阻尼(δ)、模型(M)、模态组合(MC)和土壤与结构相互作用(SSI)。利用这些变量可计算结构响应因子及相应的随机性和不确定性[13]： FSR=FSAFGMIFδFMFMCFSSI (7) (8) (9) FC和FRE也可类似地分解为不同的变量来计算。 对于抗震能力中值Am，有： Fm=Am/ASSE (10) 式中，Fm为安全因子中值。由此，可通过计算F的中值Fm及随机性和不确定性获得Am、βR和βU。 2 基于分离变量的相关性分析方法 2.1 方法原理 目前地震PSA中采用简化假设的主要原因是难以获取相关性参数，缺乏相关性分析的输入。基于分离变量的相关性分析方法(以下称分离变量法)能避免直接求解相关性参数的问题，而是将其转化为易损度分析中相关变量的识别，易损度分析人员完全有能力对此进行判断[14]。 基于分离变量的地震相关性分析原理由Reed和McCann于1985年提出[14]，通过在地震易损度评价中分离独立和相关变量来实现，与易损度评价过程紧密结合。这些变量为影响FC、FSR和FRE的各参数。通过识别这些参数是否具有相同的来源，来确定其为独立变量或相关变量。例如，假设两个不同设备位于相同厂房的相同位置，则它们具有相同的FSR及其随机性和不确定性，具有不同的FRE及其随机性和不确定性。 2.2 计算模型 利用