地震观测基础知识课件

CDSN 1960 1970 1980 1990 2000我国地震观测以模拟观测为主：400多有人值守台站 200个区域遥测子台我国数字地震观测迅速发展：1000个数据传输台站 200套应急流动观测设备 600套科学观测台阵数字地震观测起步数字观测快速发展模拟观测台网自由振荡地震面波地核地幔非均匀性各向异性全球地震活动性地方地震活动性近地震波勘探体波现代地震学的频率范围，以及有关研究对象的带宽和动态范围。J.Peterson 1993 海拉尔地震台(CDSN)台基噪声谱（J.Peterson 1993）白家疃地震台(CDSN)台基噪声谱（J.Peterson 1993）乌鲁木齐地震台(CDSN)台基噪声谱（J.Peterson 1993）地球噪声模型NLNM J.Peterson 1993NHNM J.Peterson 1993德国Gorleben附近在地表（上）和井下300米处（下）的噪声记录所得的速度功率谱密度环境背景噪声（包括海洋噪声）主要表现为面波，其幅度随深度呈现指数衰减的特点。在地表以及井下不同深度的地震计阵列所得到的短周期地震噪声记录（左）和信号记录（右）干扰源最小距离km最小距离比例系数级环境地噪声台站其他级别环境地噪声台站硬土和砂砾土基岩级（含级）以上铁路2.002.502.000.800.600.40县级以上（含县级）公路1.301.702.000.800.600.40飞机场3.005.002.000.800.600.40大型水库、湖泊10.0015.003.000.100.040.02海浪20.0020.008.000.200.100.05采石场、矿山2.503.002.000.800.600.40重型机械厂、岩石破碎机、火力发电厂、水泥厂2.503.002.000.800.600.40一般工厂、较大村落、旅游景点0.400.402.000.800.600.40大河流、江、瀑布2.503.004.000.600.400.20大型输油输气管道10.0010.002.000.600.400.2014层（含）以上高大建筑物0.200.202.000.500.300.106层（含）以下建筑物高大树木0.030.042.000.800.600.40高围栏、低树木、高灌木0.020.032.000.800.600.40注1：N级台站与干扰源之间最小距离=级台站与干扰源之间最小距离N级台站最小距离比例系数；注2：大型水库、湖泊：指库容量1010m3的水库、湖泊；注3：重型机械厂：指有大型机械、往复运动机械的工厂；注4：一般工厂：不产生明显振动感的工厂；注5：地震台站与713层建筑物的最小距离根据地震台站与6层和14层的最小距离按层数内插。宽频带或甚宽频带观测台站设计示意图 主体建设在地下，并采用了双层保温措施，具有风干扰小、温度变化小的特点型号频 带生产厂家备注FSS-31Hz 40Hz北京港震机电技术有限公司短周期地震计FBS-30.05Hz 40Hz宽频带地震计BBVS-600.017Hz 40Hz宽频带地震计DS-4A1Hz 50Hz北京市地震局短周期地震计CTS-10.0083Hz 50Hz中国地震局地震研究所甚宽频带地震计JCZ-10.0028Hz 50Hz超宽带地震计STS-20.0083Hz 50Hz甚宽频带地震计CMG-1T0.01Hz 50Hz英国GRALP公司可选规格：0.0027Hz10Hz或其它CMG-3T0.01Hz 50Hz可选规格：0.0027Hz50Hz或0.0083Hz50Hz或0.0333Hz50HzCMG-3ESP0.01Hz 50Hz可选规格：0.0083Hz50Hz或0.0333Hz50HzCMG-40T0.033Hz 50Hz可选规格：0.0083Hz50Hz或0.0188Hz50Hz或 0.05Hz50HzS-13J1Hz 100Hz美国Geotech公司机械摆，短周期地震计KS-20000.01Hz 50Hz宽频带地震计KS-540000.003Hz 16Hz加速度输出观测频带对记录信号的影响与速度成正比的宽带数字记录波形，为MOX台站记录的内华达核试验场的一次地下核试验。正的初动位移在记录BB中清晰可见，尽管信噪比很低。DAC ADC命令转发模拟信号采样示意图 频域时域频率混叠示意图ADC逐次比较式模数转换器 Delta-Sigma ADC1bit 码频谱与数字滤波器频率特性差分非线性积分非线性零位漂移增益误差ADC非量化误差人工值守地震台设备配置示意图 模拟遥测地震台设备配置示意图电信传输遥测台网中心设备配置示意图 电信传输遥测台网中心自动数据采集系统示意图 数字地震观测台网示意图 数字地震观测台站示意图 国家地震台网台站设备配置示意图 如果f(t)在区间（0，）上满足狄利克里条件（即f(t)存在有限个间断点和极限值），并且积分 收敛（或称f(t)绝对可积），则为f(t)的傅立叶变换。并有式中为角频率，傅立叶变换建立了时域和频域的对应关系。dtetfFtj)()(dttf)(),0()0(),()(21tftftfdeFtj在连续点在间断点0)()(dtetfsFstjjstdsesFjtf)(21)(正变换逆变换nnnnmmmmasasasabsbsbsbsXsYsH11101110)()()(对于一个稳定的线性动态系统，它的传递函数可以用一个有理分式形式的函数来表达：传递函数的分子多项式的根称为它的零点，分母多项式的根称为它的极点。传递函数零点、极点的分布表征了系统的动态特征。一个稳定系统的传递函数的非实极点和非实零点可能成对出现，互为共轭复数，极点位于s平面的左半平面。D=0.3D=0.707D=1.512)(22DssssH归一化的摆的二阶传递函数示例D=0.3D=0.707D=1.5121)(2DsssH归一化的摆的二阶传递函数示例nnznxzX)()(nnznx)(序列x(n)的z变换定义为 收敛的充分必要条件为 MkkNkkknxbknya00)()(MkkkNkkkzbzXzazY00)()(NkkkMkkkzazbzXzYzH00)()()(用常系数线性差分方程来表示系统的输入输出关系：当系统初始状态为零时，取z变换：系统的传递函数为：NkkMkkzdzcAzH1111)1()1()(传递函数的零极点表达：离散时间线性时不变系统的传递函数完全由它的零点、极点及常数A来决定。对于一个稳定系统，其极点应全部位于z平面单位圆内部。传递函数示例MkkzkhzH0)()(FIR数字滤波器：上图为线性相位系统；下图为最小相位系统。FIR数字滤波器示例数字滤波器示例1111zessT2216006)(ssksH309779.12542477796.9j200001sT2119990579662.09990185065.11)(zzkzHfr幅频特性dBHzHzdB示例：JCZ-1型超宽频带地震计BB通道的传递函数)()()()()(52512424123231222212121122KsKsKsKsKsKsKsKsKsKsKmssH式中，K=1000Vs/m，m=2.42061020，K11=0.024682，K12=0.00030462，K21=533.15，K22=142122，K31=667.60，K41=488.72，K51=178.88，K32=K42=K52=119422。传递函数的零点、极点表达：已知系统的传递函数H(s)可以直接令s=j代入H(s)来计算频率特性：jssHjH)()(f2)(jH)(argjH幅频特性：相频特性：IIR数字滤波器传递函数的一般表达：MkkkNkkkzazb)z(H101FIR数字滤波器传递函数的一般表达：10)()(NnnznhzH数字滤波器频率特性的计算在z平面的单位圆上进行：jezjzHeH)()(sff2)(jeH)(argjeH幅频特性：相频特性：