地震背景噪声功率谱密度及概率密度函数的计算及其Matlab实现  

Seismic Background Noise Power Spectral Density and Probability Density Function Calculation Based on Matlab

作　　者：刘庆华[1,2] 鲁来玉 常利军[1] 李宗超[1] 侯博文[2] 谢佳兴 LIU Qinghua;LU Laiyu;CHANG Lijun;LI Zongchao;HOU Bowen;XIE Jiaxing(Institute of Geophysics,China Earthquake Administration,Beijing 100081,China;Henan Earthquake Agency,Zhengzhou 450016,China)

机构地区：[1]中国地震局地球物理研究所,北京100081 [2]河南省地震局,郑州450016

出　　处：《华北地震科学》2024年第3期102-107,共6页North China Earthquake Sciences

基　　金：国家自然科学基金地震联合基金(U1839209);中国地震局地球物理研究所基本科研业务项目(DQJB19B30);河南省科技攻关(242102320039);河南省地震局第四批科技创新团队(地磁观测台站网标准化建设探索创新团队)。

摘　　要：介绍了地震背景噪声功率谱密度及概率密度函数计算方法及涉及的一些概念,以台站实际数据为例,给出了应用Matlab计算噪声功率谱密度及其概率密度函数的程序算法,该方法直接输入文本格式数据,计算台基噪声功率谱概率密度函数;以洛阳地震台24小时的连续波形记录为例,详细叙述了原始数据去仪器响应、计算功率谱密度、倍频程计算和噪声功率谱密度的概率密度函数计算实现的全过程。该方法编程简单,适用于单台噪声水平评估,在野外台站的快速勘选中具有潜在应用价值。This paper introduces the calculation methods and concepts related to the power spectral density and probability density function of seismic background noise.Taking actual data from the station as an example,a program algorithm for calculating the power spectral density and probability density function of noise using Matlab is provided.This method directly inputs text format data to calculate the probability density function of the power spectral density of background noise in a station.Taking the 24-hour continuous waveform record at Luoyang Seismic Station as an example,the entire process of implementing the removing instrument response,power spectral density calculation,octave calculation,and probability density function calculation of noise power spectral density from raw data is described in detail.This method is simple to program and suitable for evaluating the noise level of a single station.It has potential application value in the rapid selection of field stations.

关 键 词：地震背景噪声 噪声功率谱 噪声概率密度函数 倍频程 仪器响应 

分 类 号：P315.63[天文地球—地震学]