【自动p相位到达时间选择器】具有SNR输出的自动P相到达时间选择器研究（Matlab代码实现）

 👨‍🎓个人主页

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 文献来源

🌈4 Matlab代码及文章下载

---

💥1 概述

摘要提出了一种无需设置检测间隔或阈值即可提取单分量加速度或宽带速度记录中p相到达时间的新方法。phasepicker算法将信号转换为具有粘性阻尼的单自由度振荡器的响应域，然后跟踪耗散阻尼能量的变化率来拾取p波相位。

SDOF振荡器具有短的自然周期和相应的高谐振频率，比地震波中的大多数频率都要高。它还具有高阻尼比(60%的临界)。在此阻尼水平下，频率响应接近巴特沃斯最大平坦幅度滤波器，相位角保持不变。输入信号传递给振荡器的相对输入能量被转换为弹性应变能，然后作为阻尼能被阻尼元件耗散。随着时间的推移，阻尼能量产生一个平滑的包络;它在信号开始时为零，在P相到达之前为零或接近于零，并随着P波迅速增加。由于阻尼能量函数在P波开始时变化很大，因此它被用作跟踪和选择P相到达时间的度量。

PPHASEPICKER采用直方图法检测p期开始。将其性能与使用短期平均与长期平均之比的拾取技术和使用赤池信息准则寻找第一个p相位到达时间的拾取方法进行了比较。使用大量不同强度和信噪比的记录对PPHASEPICKER进行了测试，结果表明，PPHASEPICKER能够更准确地从背景噪声中挑选出真实信号的起始点，并正确区分首次到达的信号是P波(突发或脉冲)还是信号来自故障传感器。

原文摘要：

Abstract Presented is a new approach for picking P-phase arrival time in single-component acceleration or broadband velocity records without requiring detection in-
terval or threshold settings. The algorithm PPHASEPICKER transforms the signal into a
response domain of a single-degree-of-freedom (SDOF) oscillator with viscous damping
and then tracks the rate of change of dissipated damping energy to pick P-wave phases.
The SDOF oscillator has a short natural period and a correspondingly high resonant
frequency, which is higher than most frequencies in a seismic wave. It also has a high
damping ratio (60% of critical). At this damping level, the frequency response
approaches the Butterworth maximally flat magnitude filter, and phase angles are pre-
served. The relative input energy imparted to the oscillator by the input signal is con-
verted to elastic strain energy and then dissipated by the damping element as damping
energy. The damping energy yields a smooth envelope over time; it is zero in the begin-
ning of the signal, zero or near zero before the P-phase arrival, and builds up rapidly
with the P wave. Because the damping energy function changes considerably at the
onset of the P wave, it is used as a metric to track and pick the P-phase arrival time.
The PPHASEPICKER detects P-phase onset using the histogram method. Its performance
is compared with picking techniques using short-term-average to long-term-average ra-
tio, and a picking method that finds the first P-phase arrival time using the Akaike in-
formation criterion. A large set of records with various intensities and signal-to-noise
ratios is used for testing the PPHASEPICKER, and it is demonstrated that PPHASEPICKER
is able to more accurately pick the onset of genuine signals against the background noise
and to correctly distinguish between whether the first arrival is a P wave (emergent or
impulsive) or whether the signal is from a faulty sensor.

📚2 运行结果

 

 

 

 部分代码：

% Load strong-motion acceleration waveform (component HNN of station
% CE.68150 record from the 2014 M6.0 South Napa earthquake in California).

load ('strong-motion.mat');    

%% Run PphasePicker with default picking parameters
type = 'SM'; % Input type is strong motion  
pflag = 'Y'; % To plot waveform and P-phase onset

[loc] = PphasePicker(x, dt, type, pflag);
pause(2);

% Run PphasePicker with optional picking parameters
Tn = 0.01;      % undamped natural period of oscillator in second 
xi = 0.6;       % damping ratio 
nbins = 200;    % histogram bin size
o = 'to_peak'   % 'to_peak' to take segment of waveform from beginning to
                % absolute peak value
[loc, snr_db] = PphasePicker(x, dt, type, pflag, Tn, xi, nbins, o);

pause(2); 

%% Example 2 (VELOCITY WAVEFORM WITH HIGH SAMPLING RATE)
% Load micro seismic velocity waveform with 6000 samples-per-second

load ('weak-motion.mat');    

%% Run PphasePicker with optional picking parameters
type = 'na'; % no bandpass filtering  
pflag = 'Y'; % To plot waveform and P-phase onset

Tn = 0.01;      % undamped natural period of oscillator in second 
xi = 0.6;       % damping ratio 
nbins = 200;    % histogram bin size
o = 'to_peak'   % 'to_peak' to take segment of waveform from beginning to
                % absolute peak value
[loc, snr_db] = PphasePicker(x, dt, type, pflag, Tn, xi, nbins, o);

%% Example 3 (ACCELERATION WAVEFORM WITH LOW SAMPLING RATE)
% Load acceleration waveform with 80 samples-per-second

load ('acc-motion_low_sampling.mat');    

%% Run PphasePicker with optional picking parameters
type = 'na'; % no bandpass filtering  
pflag = 'Y'; % To plot waveform and P-phase onset

% USE Tn = 0.1 for records with low sampling rate lower than 100 sps 
Tn = 0.1;       % undamped natural period of oscillator in second 
xi = 0.6;       % damping ratio 
nbins = 200;    % histogram bin size
o = 'to_peak'   % 'to_peak' to take segment of waveform from beginning to
                % absolute peak value
[loc, snr_db] = PphasePicker(x, dt, type, pflag, Tn, xi, nbins, o);

🎉3 文献来源

部分理论来源于网络，如有侵权请联系删除。

🌈4 Matlab代码及文章下载