MASW之APM告警

一、

   二所APM告警,和MASW告警同时计算,告警计算时不分优先级别,告警时显示“APM”。客户现场使用:MASW和APM重叠情况下,一般在低高告警对应APM划一块MASW抑制区,这样就只有APM生效。否则的话会同时有MASW和APM告警显示。


二、

   川大ATC系统APM告警使用适配参数dat_maswapplast,不使用ini配置文件apmparam.ini(已不用,可以删除),APM告警时在SDD显示MSAW(因为两项告警规范里最后进近告警APM归到MSAW告警一类中),MASW和APM区域重叠情况下,MASW预警、告警会失效,只按APM监测低高告警。

来自 “ ITPUB博客 ” ,链接:http://blog.itpub.net/7970627/viewspace-2157908/,如需转载,请注明出处,否则将追究法律责任。

转载于:http://blog.itpub.net/7970627/viewspace-2157908/

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
import numpy as np import matplotlib.pyplot as plt from obspy import read # 读取面波数据并画图。 st = read('MASW_DATA/Sample_Data/*.SAC') dt = st[0].stats.delta data = [] scale = 0.05 dx = 2 plt.figure(figsize=(8, 6)) for i, tr in enumerate(st): d = tr.data data.append(d) t = np.arange(len(d)) * dt plt.plot(t, d*scale+(i+1)*dx, lw=1, color='b') plt.xlabel('Time (s)') plt.ylabel('Offset (m)') plt.tight_layout() plt.savefig('Surface_wave.png') plt.show() # 二维FFT。 d = np.array(data) n = len(d[0]) # m为空间方向的采样点数,m增大可以让FK谱光滑一点,以达到插值效果。 m = len(d[:, 0]) * 5 D = np.zeros((m, n)) D[:len(d[:, 0])] = d # 时间采样率。 fs = 1 / dt # 空间采样率 xs = 1 / dx # 频率 (赫兹)。 f = np.arange(-n//2, n//2) * fs / (n-1) # 波数 (每米)。 k = 2 * np.pi * np.arange(-m//2, m//2) * xs / (m-1) # 二维FFT。 fk = np.fft.fft2(D) # 作图。 pmin = -10 P = abs(np.fft.fftshift(fk)); P /= P.max(); P = 10 * np.log10(P) P2 = abs(fk); P2 /= P2.max(); P2 = 10 * np.log10(P2) plt.figure(figsize=(11, 8)) plt.subplot(221) plt.pcolormesh(f, k, P2, cmap='magma', vmin=pmin, vmax=0) plt.xlabel('Frequency (s$^{-1}$)') plt.ylabel('Wave number (2$\pi$m$^{-1}$)') plt.subplot(222) plt.pcolormesh(f, k, P, cmap='magma', vmin=pmin, vmax=0) plt.plot([f[n//2], f[-1], f[-1], f[n//2], f[n//2]], [k[0], k[0], k[m//2], k[m//2], k[0]], lw=2, ls='--', color='r') plt.xlabel('Frequency (s$^{-1}$)') plt.ylabel('Wave number (m$^{-1}$)') plt.subplot(223) plt.pcolormesh(f[n//2:], k[:m//2], P[:m//2, n//2:], cmap='magma', vmin=pmin, vmax=0) plt.xlabel('Frequency (s$^{-1}$)') plt.ylabel('Wave number (m$^{-1}$)') plt.subplot(224) plt.pcolormesh(f[n//2:], abs(k[:m//2][::-1]), P[:m//2, n//2:][::-1], cmap='magma', vmin=pmin, vmax=0) cbar = plt.colorbar() cbar.set_label(r'FK spectra (dB)') plt.xlim(0, 100) plt.xlabel('Frequency (s$^{-1}$)') plt.ylabel('Wave number (m$^{-1}$)') plt.tight_layout() plt.show()
最新发布
06-13

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值