文章快速检索     高级检索
  大地测量与地球动力学  2024, Vol. 44 Issue (9): 985-990  DOI: 10.14075/j.jgg.2023.12.599

引用本文  

吴雄伟, 杨江, 夏界宁, 等. 红沿河核电厂地震仪表系统震后数据分析[J]. 大地测量与地球动力学, 2024, 44(9): 985-990.
WU Xiongwei, YANG Jiang, XIA Jiening, et al. Post-Earthquake Data Analysis of Seismic Instrument System at Hongyanhe Nuclear Power Plant[J]. Journal of Geodesy and Geodynamics, 2024, 44(9): 985-990.

项目来源

中国地震局地震研究所和应急管理部国家自然灾害防治研究院基本科研业务费(IS202216317);湖北省重点研发计划(2022BAD059);中国地震局地震科技星火计划(XH231702YB)。

Foundation support

Scientific Research Fund of Institute of Seismology, CEA and National Institute of Natural Hazards, MEM, No. IS202216317; Key Research and Development Program of Hubei Province, No. 2022BAD059; The Spark Program of Earthquake Technology of CEA, No. XH231702YB.

第一作者简介

吴雄伟,工程师,主要从事核电地震仪表的设计研究,E-mail:wuxiongwei@whsii.com

About the first author

WU Xiongwei, engineer, majors in the development of seismic instrument of nuclear power station, E-mail: wuxiongwei@whsii.com.

文章历史

收稿日期:2023-12-27
红沿河核电厂地震仪表系统震后数据分析
吴雄伟1,2,3,4     杨江1,2,3,4     夏界宁1,2,3,4     范涛1,2,3,4     
1. 中国地震局地震研究所,武汉市洪山侧路40号,430071;
2. 地震预警湖北省重点实验室,武汉市洪山侧路40号,430071;
3. 武汉地震科学仪器研究院有限公司,湖北省咸宁市青龙路11号,437000;
4. 湖北省重大工程地震监测与预警处置工程技术研究中心,武汉市洪山侧路40号,430071
摘要:2023-08-23红沿河核电厂地震仪表系统成功记录到大连市普兰店区发生的4.6级地震,这是我国首次完整记录到核电厂附近中强震数据。对该系统记录到的地震数据进行分析,时程信号波形回放表明,本次地震记录数据清晰完整,自由场数据峰值与地震衰减经验公式计算结果相符。系统7个监测点加速度峰值对比分析表明,核电厂厂房对地震加速度信号具有放大效应,放大系数与建筑物标高正相关。加速度峰值数据频谱分析结果显示,核电厂厂房地震响应数据的卓越频率主要集中在10~20 Hz范围内,该范围内的地震加速度信号对厂房仍具有较大破坏性,这一结论与核电厂地震仪表准则NB/T 20076-2012中规定的地震触发频带(1~10 Hz)不符。鉴于我国核电厂地震仪表系统的地震触发滤波通频带设定范围为1~10 Hz,这一缺陷会降低地震仪表系统纵深安全防御的性能。
关键词核电厂地震仪表系统地震响应卓越频率

2023-08-23辽宁省大连市普兰店区发生4.6级地震,震中(39.42°N,122.2°E)距离红沿河核电厂厂区74 km,厂区有明显震感。地震发生后,核电厂3、4号机组及5、6号机组的地震仪表系统没有被触发启动记录,只有1、2号机组布设的地震仪表系统完整记录到地震的加速度时程数据,包含从地震开始到结束的整个过程。此次地震记录是我国核电厂首次记录到的中强震数据,具有重大的科研参考价值。本文分析了记录数据的完整性及合理性,并根据数据分析结果对核电厂配套地震仪表系统设计参数及安装方式等提出改进建议,以确保核电厂的地震防御安全。

1 地震仪表系统结构

根据核电厂地震规范要求,红沿河核电厂安装了3套地震仪表系统,用于对核电厂区的地震信息进行实时监测。当监测到破坏性地震事件时,向核电厂主控室发送实时地震报警信号,为核电厂提供纵深安全防御。3套系统的电气机柜分别安装在1、2号机组、3、4号机组及5、6机组的辅助厂房设备间(每2个机组共用一套地震仪表系统),3套系统在设备组成和监测点的布设方位上完全相同:系统主要由电气机柜和7台三轴向加速度传感器组成,电气机柜内安装了7台地震记录仪,7台加速度传感器分别布设在厂区自由场(001MZ)、反应堆厂房安全壳标高-3.5 m(002MZ)、4.65 m(007MZ)与20 m(003MZ)处及辅助厂房标高-3.45 m(004MZ)和0 m(008MZ)处,图 1为1、2号机组地震仪表系统各监测点的加速度传感器布设位置示意图[1-2]

图 1 核电厂加速度传感器布设 Fig. 1 Layout of acceleration sensors in nuclear power plants

这7台加速度传感器与电气机柜内部的7台地震记录仪一一对应,二者通过仪表电缆相连。其中,与002MZ和004MZ加速度传感器连接的2台地震记录仪被设置为地震触发器(触发阈值设置为10 Gal)[3],当这2台加速度传感器三轴向中任一轴向的加速度峰值超过地震报警触发阈值时,对应的2台地震记录仪自动通过网络接口发送触发信号,启动7台地震记录仪同步记录地震触发前10 s和触发条件失效后30 s地震数据,即核电厂地震仪表系统采用的是触发记录模式,与强震仪数据记录方式相同[4-5]。系统电气结构如图 2所示。

图 2 地震仪表系统电气结构 Fig. 2 Electrical structure of seismic instrument system
2 地震数据合理性分析

现阶段,我国民用核电厂均布设在东南沿海一带,选址前均进行了严格的地质调查。从秦山核电站开始运营至本次普兰店地震发生前,我国所有核电厂安装的地震监测设备均未记录到地震数据。地震发生后,红沿河核电厂3套地震仪表系统中仅有1套系统被触发,并记录到地震数据,另外2套系统未被触发,且核电厂厂区周边较近范围内也没有设置地震台站,无法提供参照数据进行对比。

对于红沿河核电厂3套地震仪表系统,其设计参数均满足强震监测规范要求,系统均按照生产工艺进行了严格的安装调试,并正常运行多年,期间也均依照维修保养流程进行了功能性测试与校验,对加速度传感器及地震记录仪的性能进行了周期性测试校验,设备各项指标均符合技术规范要求。可确认,核电厂3套地震仪表系统在硬件和软件功能上满足地震监测与记录要求。

2.1 地震数据完整性分析

地震结束后,通过专用软件对1、2号机组地震记录仪记录的地震数据进行数据回放,7个监测点的数据波形见图 3。结果显示,7台地震记录仪同步记录了地震触发前10 s及触发条件失效后30 s的连续时程数据,包含从地震开始到结束的整个过程,无数据丢失现象,也无异常数据突跳点,数据采集与存储完整。每台加速度传感器地震前与地震后的本底噪声相同,噪声峰值约为1~3 Gal,远小于地震触发报警阈值,与核电厂安装点的噪声环境要求相符[6]

图 3 地震加速度实时波形 Fig. 3 Real time waveform of seismic acceleration

地震仪表系统软件显示,触发地震事件记录的地震触发器对应监测点传感器为004MZ,该监测点Y轴向加速度峰值为13.1 Gal,超过了地震事件设定的触发阈值(10 Gal),满足地震事件记录触发条件;而另一监测点(002MZ)仅Y轴向加速度峰值超过触发阈值,2台传感器加速度峰值均超过设定阈值,地震触发记录合理有效。考虑到系统安装环境的差异性,3、4号机组及5、6号机组2套地震仪表系统未触发记录存在较大的合理性。

2.2 地震动加速度峰值衰减对比分析

红沿河核电厂每套地震仪表系统仅设置了单个自由场监测点,厂区周围无其他地震台站记录数据进行对比,因此有必要对自由场数据的合理性进行进一步分析。对于强震数据缺乏的地区,可依据震中区地震加速度记录数据的峰值和震级,选用对应的地震衰减关系经验公式,换算出该地区的地震加速度峰值[7-9]。本文选用文献[9]中经验公式对普兰店地震到达红沿河核电厂的加速度峰值进行计算分析:

$ \begin{gathered} \log _{10}\left(a_E\right)=c_1+c_2 \cdot M+c_3 \cdot M^2+ \\ c_4 \cdot \log _{10}\left(R+c_5 \cdot \mathrm{e}^{c_6} \cdot M\right) \end{gathered} $ (1)

式中,c1~c6为回归系数,由具体地质结构确定;R为震中距,单位km;M为震级。参照辽宁省地震局提供的普兰店地震信息可知,震级为4.6,震中距离红沿河核电厂74 km。鉴于红沿河核电厂核反应堆等建筑物均建立在基岩上,属于Ⅰ类场址,参考文献[9]依照表 1选取c1~c6参数值,代入式(1)可得:水平向加速度峰值为8.10 Gal,垂直向加速度峰值为3.45 Gal。与实地现场自由场数据对比可知,经衰减换算经验公式计算得到的加速度峰值与系统记录到的数值符合度较高,可确认地震仪表系统记录到的加速度数据合理有效。

表 1 加速度衰减计算公式参数 Tab. 1 Parameters of acceleration attenuation calculation formula
3 数据分析 3.1 放大效应分析说明

表 2为各监测点加速度峰值信息,表中X轴为平行于反应堆主轴方向,Y轴为垂直于反应堆主轴方向。由表可知,标高越高,加速度峰值越大。其中,安装在0 m标高附近的6台加速度传感器地震响应峰值数据接近,而核反应堆厂房20 m标高处监测点,无论是水平向还是垂直向的加速度峰值数据均较大,说明该监测点对地震加速度的放大效应明显,标高与放大效应的变化规律符合逻辑关系。放大系数为建筑物上布设的监测点加速度峰值与自由场监测点加速度峰值之比[10]。依照该定义可知,各监测点的放大系数如表 2所示,数据显示,核电厂建筑物对地震加速度响应具有放大效应。对于标高在0 m以下的监测点,其放大系数多数小于1,为负放大系数;对于0 m以上监测点,其放大系数均大于1,为正放大系数,这一结论符合建筑物对地震的放大效应特征。

表 2 监测点加速度峰值信息 Tab. 2 Peak acceleration information of monitoring points

由于核电厂建筑物的放大效应与标高为正相关,对于安装在安全壳基础上作为地震触发器的2台加速度传感器,其标高最低,地震放大效应最小,地震事件触发几率也最小。鉴于本次地震实际记录到的数据与地震触发阈值十分接近,同时考虑到3、4号机组与5、6号机组地震仪表系统安装环境的区别及监测点放大效应的差异化,极有可能导致系统记录到的加速度峰值小于触发阈值,这2套系统没有被触发启动记录存在一定的合理性。由于放大效应的影响,将低标高的监测点设置为地震触发器,会导致小震级地震被漏记的概率增加。后续系统改造升级时应考虑改进记录方式,增设高标高监测点作为地震触发器,或增加24 h连续记录功能,保证数据记录的可靠性。

3.2 傅里叶谱分析

对地震数据进行傅里叶谱分析,结果如图 4所示。自由场监测点傅里叶谱(图 4(a))分析结果表明,本次记录到的地震波频率主要集中在1~10 Hz频段,10~20 Hz以上频段迅速衰减,即本次地震的加速度信号频率主要集中在1~10 Hz频段范围内;安装在厂房内的6个监测点傅里叶谱(图 4(b)~4(g))分析结果表明,本次记录到的地震波频率主要集中在10~20 Hz,10~20 Hz频段高频信号区域衰减较慢。对比自由场数据同频段数据可知,10~20 Hz频段内核电厂建筑物有共振现象,即本次地震信号的能量主要集中在1~10 Hz频段范围内,但与核电厂建筑物存在共振效应,需考虑地震信号10~20 Hz频段高频信号的破坏性影响;另外,7个监测点均存在50 Hz频率突跳现象,应为工频干扰所致[11-12]

图 4 各监测点傅里叶谱 Fig. 4 Fourier spectra of monitoring points
3.3 加速度反应谱分析

对地震数据进行反应谱计算分析,阻尼比选取5%,反应谱如图 5所示,各监测点的卓越频率点见表 3。各监测点反应谱对比结果表明,反应谱峰值主要集中在1~20 Hz范围内,20 Hz以上迅速衰减,与傅里叶谱分析结论一致;50 Hz频率点处无峰值突跳出现,表明系统已对工频噪声的影响进行了抑制;2个标高-3.5 m及1个标高-3.45 m的监测点水平方向的卓越频率相同,表明数据采样的一致性较好[13-14]。对各监测点的卓越频率峰值进行对比可以看出,监测点的标高越高,其反应谱峰值越大,即核电厂建筑物对地震波具有明显的放大效应,且与标高正相关。

图 5 各监测点加速度反应谱 Fig. 5 Acceleration response spectra of monitoring points

表 3 反应谱卓越频率 Tab. 3 Response spectrum predominant frequency list

对比安装在厂房内的监测点反应谱可知,002MZ、004MZ、006MZ及008MZ的卓越频率点在10~20 Hz范围内,该频段反应谱值并不比1~10 Hz频段反应谱值小,表明本次地震的地震波在10~20 Hz频率范围内依旧具有较大的破坏性。这一结论与NB/T 20076-2012地震仪表准则[15]中规定的地震触发器的地震滤波范围(1~10 Hz)并不相符。依照设计方案,地震仪表系统对地震安全停堆保护报警(SSE)的判别依旧沿用1~10 Hz的滤波通带范围,这会导致10~20 Hz范围内的地震加速度采样信号被衰减,使破坏性地震的漏报率增加。后续在地震仪表系统改造过程中,应考虑适度展宽地震触发器的滤波通带范围。

4 结语

对红沿河核电厂地震仪表系统记录到的地震数据进行分析,结果显示,记录的地震数据完整、合理,首次直接证明红沿河核电厂地震仪表系统具备中强震监测能力。记录数据显示,核电厂厂房对地震加速度信号具有放大效应,放大系数与标高正相关。地震加速度对核电厂的破坏性信号卓越频率主要集中在10~20 Hz范围内,表明10~20 Hz以内的地震加速度信号对核电厂建筑物同样具备破坏性,这一结论与相关规范给出的地震数据采集滤波通带范围不符,有可能导致地震仪表系统对地震停堆识别的漏报。建议后期在对地震仪表系统进行改造的过程中,将高标高的监测点设定为地震触发器,或添加系统24 h数据记录功能,增强系统地震数据记录能力;同时,对地震触发阈值判别的滤波通带范围进行适度展宽,降低地震报警判别的漏报率。

参考文献
[1]
陈志高, 吴雄伟, 姚立民. 岭澳核电站地震监测系统[J]. 大地测量与地球动力学, 2012, 32(增1): 170-172 (Chen Zhigao, Wu Xiongwei, Yao Limin. Earthquake Monitoring System of Ling'ao Nuclear Power Station[J]. Journal of Geodesy and Geodynamics, 2012, 32(S1): 170-172) (0)
[2]
肖志慧. CPR1000核电地震仪表系统大修年检工艺优化[J]. 中国设备工程, 2018(14): 56-57 (Xiao Zhihui. Optimization of Annual Inspection Process for Overhaul of CPR1000 Nuclear Power Seismic Instrument System[J]. China Plant Engineering, 2018(14): 56-57 DOI:10.3969/j.issn.1671-0711.2018.14.033) (0)
[3]
陈志高, 黄俊. 核电站地震监测系统及其开发进展[J]. 大地测量与地球动力学, 2015, 35(6): 1065-1068 (Chen Zhigao, Huang Jun. Development and Processes of Seismic Monitoring Systems for Nuclear Power Plants[J]. Journal of Geodesy and Geodynamics, 2015, 35(6): 1065-1068) (0)
[4]
王雷, 周琪, 高峰, 等. 新型数字强震仪关键技术分析[J]. 地震工程与工程振动, 2011, 31(6): 153-158 (Wang Lei, Zhou Qi, Gao Feng, et al. New Technology Used in Strong Motion Seismograph[J]. Journal of Earthquake Engineering and Engineering Vibration, 2011, 31(6): 153-158) (0)
[5]
杨江, 范涛, 黄蔚, 等. 基于MEMS传感器的工程强震动加速度计研制[J]. 自然灾害学报, 2022, 31(4): 181-189 (Yang Jiang, Fan Tao, Huang Wei, et al. Development of Engineering Strong Vibration Accelerometer Based on MEMS Sensor[J]. Journal of Natural Disasters, 2022, 31(4): 181-189) (0)
[6]
The American Society of Mechanical Engineers OM Committee. Operation and Maintenance of Nuclear Power Plants: ASME OM-2022[S]. New York: The American Society of Mechanical Engineers, 2023 (0)
[7]
汪素云, 俞言祥, 高阿甲, 等. 中国分区地震动衰减关系的确定[J]. 中国地震, 2000, 16(2): 99-106 (Wang Suyun, Yu Yanxiang, Gao Ajia, et al. Development of Attenuation Relations for Ground Motion in China[J]. Earthquake Research in China, 2000, 16(2): 99-106 DOI:10.3969/j.issn.1001-4683.2000.02.001) (0)
[8]
俞言祥, 李山有, 肖亮. 为新区划图编制所建立的地震动衰减关系[J]. 震灾防御技术, 2013, 8(1): 24-33 (Yu Yanxiang, Li Shanyou, Xiao Liang. Development of Ground Motion Attenuation Relations for the New Seismic Hazard Map of China[J]. Technology for Earthquake Disaster Prevention, 2013, 8(1): 24-33 DOI:10.3969/j.issn.1673-5722.2013.01.003) (0)
[9]
席远, 耿淑伟, 董满生, 等. 基于现行规范场地划分标准的地震动峰值加速度衰减关系研究[J]. 合肥工业大学学报: 自然科学版, 2012, 35(9): 1264-1268 (Xi Yuan, Geng Shuwei, Dong Mansheng, et al. Attenuation Relation of Strong Ground Motion PGA Based on Classified Site Criterion in Current Seismic Design Code[J]. Journal of Hefei University of Technology(Natural Science), 2012, 35(9): 1264-1268 DOI:10.3969/j.issn.1003-5060.2012.09.026) (0)
[10]
吴亦恺, 丁海平. 利用强震记录估计场地放大效应的一种方法[J]. 苏州科技大学学报: 工程技术版, 2023, 36(1): 15-21 (Wu Yikai, Ding Haiping. A Method for Estimating Site Amplification Effect Using Strong Earthquake Records[J]. Journal of Suzhou University of Science and Technology(Engineering and Technology Edition), 2023, 36(1): 15-21) (0)
[11]
孙常青, 朱振家, 潘章容, 等. 甘肃地震预警测站典型异常波形的频谱分析[J]. 地震工程学报, 2023, 45(3): 614-624 (Sun Changqing, Zhu Zhenjia, Pan Zhangrong, et al. Spectrum Analysis of the Typical Abnormal Waveforms of Earthquake Early Warning Stations in Gansu Province[J]. China Earthquake Engineering Journal, 2023, 45(3): 614-624) (0)
[12]
席继楼, 赵家骝, 关歆莹, 等. 地电场观测中的工频干扰抑制方法研究[J]. 地震, 2015, 35(4): 53-63 (Xi Jilou, Zhao Jialiu, Guan Xinying, et al. Methods for Restraining Power Frequency Interferences in Geo-Electric Field Observation[J]. Earthquake, 2015, 35(4): 53-63) (0)
[13]
陈永新, 迟明杰, 李小军. 基于强震动记录确定的场地卓越周期[J]. 地震学报, 2016, 38(1): 138-145 (Chen Yongxin, Chi Mingjie, Li Xiaojun. Determination of Site Dominant Period Based on Strong Motion Records[J]. Acta Seismologica Sinica, 2016, 38(1): 138-145) (0)
[14]
侯瑞彬. 基于观测记录的场地卓越周期与场地基本周期对地震动放大模型的影响研究[J]. 地震工程与工程振动, 2021, 41(6): 202-210 (Hou Ruibin. The Effect of Site Predominant Period Derived from Seismic Records and Site Fundamental Period on Seismic Site Amplification Model[J]. Earthquake Engineering and Engineering Dynamics, 2021, 41(6): 202-210) (0)
[15]
国家能源局. 核电厂地震仪表准则: NB/T 20076-2012[S]. 北京: 原子能出版社, 2012 (National Energy Bureau of the People's Republic of China. Earthquake Instrumentation Criteria for Nuclear Power Plants: NB/T 20076-2012[S]. Beijing: Atomic Energy Press, 2012) (0)
Post-Earthquake Data Analysis of Seismic Instrument System at Hongyanhe Nuclear Power Plant
WU Xiongwei1,2,3,4     YANG Jiang1,2,3,4     XIA Jiening1,2,3,4     FAN Tao1,2,3,4     
1. Institute of Seismology, CEA, 40 Hongshance Road, Wuhan 430071, China;
2. Hubei Key Laboratory of Earthquake Early Waring, 40 Hongshance Road, Wuhan 430071, China;
3. Wuhan Institute of Seismic Scientific Instruments Co Ltd, 11 Qinglong Road, Xianning 437000, China;
4. Engineering Technology Research Center for Earthquake Monitoring and Early Warning Disposal of Major Projects in Hubei Province, 40 Hongshance Road, Wuhan 430071, China
Abstract: On August 23, 2023, the seismic instrument system of the Hongyanhe Nuclear Power Plant successfully recorded a M4.6 earthquake in Pulandian, Dalian. This is the first complete strong earthquake record around the nuclear power plant area in China. We analyze the seismic data recorded by the system, and the data playback shows that the seismic recorded data was clear and complete, and the peak value of the free field monitoring data is consistent with the calculation result of the seismic attenuation empirical formula. Comparing the peak acceleration values of seven monitoring points data, the result shows that the nuclear power plant building has an amplification effect on seismic acceleration signals, and the amplification coefficient is positively correlated with building elevation. The conclusion of spectral analysis of peak acceleration data shows that predominant frequency points of seismic response data in nuclear power plant buildings are mainly concentrated in the range of 10 to 20 Hz, and seismic acceleration signals in the range of 10 to 20 Hz also have certain destructive effects on the building. This conclusion is inconsistent with the seismic triggering filtering range (1 to 10 Hz) specified in NB/T 20076-2012 seismic instrument criteria, due to the fact that the seismic trigger recognition algorithm of the seismic instrument system is set in a filtering range of 1 to 10 Hz. This defect will reduce the performance of the seismic instrument system's depth safety defense.
Key words: nuclear power plant; seismic instrument system; seismic response; predominant frequency