确定海底地震仪的水平方位Determine the horizontal orientation of OBS

引言

陆地地震仪在布设时，需要将其N分量指向地理北向。然而海底地震仪（ocean bottom seismograph）在布设时是自由沉入海底，其N分量（以下用BHN表示）的指向是随机的，因此需要通过数据处理的方法来确定海底地震仪的水平方位。

方法

这里我们使用(Stachnik et al., 2012)的方法来确定OBS的方位角，该方法是基于瑞利波（Rayleigh wave）的极化实现的。
瑞利波在垂向(Z)和径向(R)上存在90°的相位差，因此我们可以通过旋转水平分量来与垂直分量的希尔伯特变换(Z Hilbert)求相关性，从而确定径向的方向，再由径向(baz)和旋转的角度(θ)确定BHN的方向(X)，如下图所示。

对于一个三分量地震，通过地震和传感器的位置可以求出反方位角(baz)。通过旋转水平分量与垂直分量的希尔伯特变换求相关性，相关性最大的角度θ即为径向。通过反方位角(baz)和旋转角(θ)即可求出BHN的方位(X)。

例子

我们使用上述方法来确定台站B01的方位角。如下图，首先挑选台站记录到的一个三分量地震波形，通过地震和台站的位置算出反方位角(baz)为262°，接着，从波形中截取瑞利波，对水平分量的瑞利波进行旋转。

数据来自Mantle serpentinization and water cycling through the Mariana Trench and Forearc: XF_2012

对每一个旋转的水平分量与垂直分量的希尔伯特变换(Z Hilbert)求相关性，如下图所示，相关性最大的方向即为径向。对于这个地震，得到径向的旋转角(θ)为241°，即从BHN的方向顺时针旋转241°得到径向的方向：$BHR=BHN\times cos(\theta )+BHE\times sin(\theta )$

这里我们使用两种方法求相关性：
1. Normalized cross correlation
$$C_{zr}=\frac{S_{zr}}{\sqrt{S_{zz}{S}_{zr}}}.$$
2. Second normalized cross correlation
$$C_{zr}^{\ast }=\frac{S_{zr}}{S_{zz}}.$$

• $C_{zr}$是一个真正的归一化互相关，它的值在-1到1之间，但是其在最大值附近具有相似值，难以确定最大值。
• $C_{zr}^{\ast }$具有明确定义的最大值，但它的值是无限的。

我们将三分量地震波形，径向(R)以及垂向的希尔伯特变换(Z Hilbert)对比可以发现，垂向(Z)与垂向的希尔伯特变换(Z Hilbert)之间存在90°相位差，而Z的希尔伯特变换(Z Hilbert)与径向(R)相关性非常高，说明我们找到了正确的径向。

因此，通过这个地震确定的水平方位角为：

$$X=baz-\theta =262°-241°=21°$$

结果

我们可以通过台站布设期间记录的多个地震来求得更准确的结果：

项目地址

该程序命名为orient，托管在Github上：https://github.com/SeisPiano/orient
希望可以帮到你，欢迎点击star🌟，谢谢(◍•ᴗ•◍)

引用

If you use orient in your work, please cite:
Liu Dan, Ting Yang, Ba Manh Le, Yuechu Wu, Yizhi Wang, Xinfeng Huang, Haoran Du, Jian Wang, Yongshun Chen. (2022). Seismometer-detached broadband ocean bottom seismograph (OBS): development, test, and data quality analysis. Chinese Journal of Geophysics, 65(7): 2560-2572, https://doi.org/10.6038/cjg2022P0441

参考文献

• Stachnik, J. C., Sheehan, A. F., Zietlow, D. W., Yang, Z., Collins, J., Ferris, A. (2012) Determination of New Zealand Ocean Bottom Seismometer Orientation via Rayleigh-Wave Polarization. Seismological Research Letters, 83(4): 704–713. https://doi.org/10.1785/0220110128