关于地动噪声功率谱的讨论.ppt

关于测定台站噪声功率谱密度的建议 1 给出约定，明确用于计算噪声功率谱的信号长度、分段数、相邻段重叠百分比、窗函数、以及谱平滑点数，以使各台站噪声功率谱估计有统一标准。 2 统一计算功率谱和RMS振幅时对PSD定义的约定，例如，采用单边PSD定义，计算RMS振幅时根号内不再乘以2，对《GB/T 19531.1-2004 地震台站观测环境技术要求 第1部分：测震》中的描述和公式加以小修，予以明确。 3 对计算的PSD，明确给出95%（或90%）置信度误差估计值（如Peterson，1993）或在PSD曲线图上给出给定置信度的误差限线段长度（如Holcomb，1989）。 关于地动噪声功率谱的讨论 赵仲和 2010-12-03 主要内容 噪声功率谱密度（PSD）定义的一致性 PSD估计值的误差估计 减小PSD估计值方差的措施 关于测定台站噪声功率谱密度（PSD）的建议 NORSAR 功率谱密度(http://www.norsardata.no/NDC/spectraplot/) 使用Welch法对每个通道计算每小时功率谱密度 (PSD) : 把1小时时间窗分成200秒长的时段，相邻段重叠100秒 (即总共35段). 计算每段的周期图（periodogram） (即 FFT并取振幅平方). 对35个周期图取平均，得到每小时的 PSD. 没有对地震图进行仪器校正，所以以康特（ count）为单位.因此， PSD 的单位是counts**2/Hz We display the logarithm (log10) of the PSD in monthly panels (see example for NAO00_SZ above). The x-axis is time in units of days the y-axis is frequency in [Hz] and the PSD is color-coded (blue: low amplitude, red. high amplitude). In order to increase the dynamic resolution of the plots the average PSD (computed from all hourly PSDs in the panel) has been substracted. The average PSD is overlaid on the plot as a red curve. Its scale is indicated with red markers on top and on the bottom of the plot. 讨论1：噪声功率谱密度定义的一致性 数学定义 自功率谱密度函数由数据时间序列的付里叶变换与它的共轭相乘得到。…基本运算如下： x(nΔt)→X(mΔf) （1） Z(mΔf)=X(mΔf)? X*(mΔf) （2） （GB 10084-88《振动、冲击数据分析和表示方法》） 其中x(nΔt)为信号，n 为数据点序号， Δt为采样间隔， X(mΔf)为其付里叶变换，m为频率序号， Δf为频率间隔， Z(mΔf)为功率谱，*表示复数的共轭 GB/T 19531.1-2004 《在地震台站观测环境技术要求 第1部分：测震》中引用这一定义 地动噪声的均方根振幅 根据上述定义，得到地动噪声的均方根（RMS）振幅 aRMS={2P x (fu-fl)}1/2 (NMSOP 4.11) 其中P为（1）式定义的功率谱密度，fu和fl分别为计算均方根振幅的频段的两端频率。 在GB/T 19531.1-2004 《在地震台站观测环境技术要求 第1部分：测震》附录A中采用这一定义。 （NMSOP：P. Bormann主编《新地震观测实践手册》） 地球正常噪声新模型（Peterson,1993） 由美国USGS（美国地质调查局）的J.Peterson及其所领导的研究小组，观测和研究了全世界各地正常地球噪声，得到了新的全球公认的地球正常噪声新模型，包括地球高噪声新模型（NHNM）和地球低噪声新模型（NLNM）。正常地噪声水平在高噪声模型（NHNM）和低噪声模型（NLNM）之间，（见图A.1）。图A.1是Peterson用全球范围75个不同地点的地面运动加速度的功率谱密度曲线集，由它的包络得到地球正常噪声模型。 （GB/T 1953