High-resolution crustal structure in the Songliao basin
-
摘要: 基于东北地区已有的宽频带流动台阵远震数据，利用波场延拓和分解的H-β网格搜索法，对松辽盆地的沉积层及地壳结构进行了深入分析。结果显示：松辽盆地的沉积层厚度为0.2—2.5 km，整体呈现中央坳陷区厚、边缘薄且西南地区最薄的分布特征；研究区地壳较薄，厚度介于24—34 km之间，其横向变化特征与沉积层厚度分布具有一定的对应性。依据沉积层和地壳的厚度计算了地壳伸展系数，其平均值接近于以往接收函数研究估测的岩石圈伸展因子。因此，本文推测松辽盆地在伸展构造过程中，其地壳和岩石圈的减薄以纯剪切模式为主。此外，松辽盆地具有较高的地壳平均波速比vP/vS，暗示盆地下方岩石圈地幔的减薄过程中可能存在岩浆的底侵作用。Abstract: High-resolution shallow crustal structure beneath the Songliao basin of Northeast （NE） China has obvious economic and scientific significance. To constrain the sediment and crustal structure of the Songliao basin, H-β grid search method based on wavefield downward continuation and decomposition is used with teleseismic data recorded from portable broadband seismic arrays in the NE China. The results show that the estimated sediment thickness is 0.2−2.5 km, and becomes thinner from the central depression toward the margin of the basin, with the thinnest sediment in the southwestern region. The crustal thickness varies from 24 km to 34 km, and its lateral variation correlates with the distribution of sedimentary thicknesses to a certain extent. The crustal stretching factor is calculated from the sedimentary and crustal thicknesses, with an average close to the lithospheric stretching estimation from previous receiver function studies. Thus, we infer that the thinning of the crust and lithosphere is dominated by pure shear mode during the extensional process of the Songliao basin. Moreover, the Songliao basin has a high crustal vP/vS ratio, indicating possible magmatic underplating during the lithospheric extension beneath the Songliao basin.
-
Keywords:
- Songliao basin /
- H-β technique /
- sediment thickness /
- crustal structure /
- vP/vS ratio
-
引言
相对于陆地，海洋的地震观测要困难得多，而海洋面积占地球表面积71%，这极大限制了全球的天然地震观测。近30年来，尽管海底地震仪（ocean bottom seismograph，缩写为OBS）的应用在很大程度上弥补了海底地震观测这一空缺（Moore et al，1981；Collins et al，2001；Stephen et al，2003），但OBS数据相对陆地台站通常表现出更高的噪声水平，尤其是在对反演海洋地下深部结构非常重要的低频段（Webb，1988，1998；Webb，Crawford，1999）。海洋的复杂环境特别是各种水波，包括底流、潮汐、内波、表面水波、海洋生物等在内的各种因素共同导致了OBS的高噪声数据。降低OBS数据的噪声水平、提高其记录地震波形的信噪比对于地震学及海洋地球物理学的发展具有重要意义。
一方面，可以依靠设计和研发更高性能的OBS来提高海底地震数据的质量。例如：将传感器模块和仪器本体分离的分体式OBS技术（Duennebier et al，1981）；将传感器模块浅埋入海底沉积层的技术（Wooding et al，2001；Shiobara et al，2013）；屏蔽OBS地震传感器（即给地震传感器增加外部屏障，使水流不能直接作用于传感器）的技术（刘丹等，2022）等。这些技术的进步无疑推动了OBS的发展，使得一些大型被动源OBS观测计划得以实施。分体式设计使得地震传感器模块相对独立，可降低仪器本体噪声对地震数据的干扰。这一设计现已经被广泛采用，成功的被动源OBS观测计划，如Cascadia Initiative （Toomey et al，2014），RHUM-RUM （Stähler et al，2016）和ORCA （Eilon et al，2022）等都是应用这种类型的仪器。利用重物或遥控无人潜水器（remote operated vehicle，缩写为ROV）将地震传感器埋入沉积层虽然可以显著降低OBS的噪声水平，但同时大大增加了布放成本，因此这种技术目前只适用于单个OBS实验（如：Wooding et al，2001；Shiobara et al，2013），难以应用于大型被动源OBS观测台阵。利用自埋式机器人技术，使OBS的传感器模块可以自动浅埋入沉积层的自埋式OBS具有广阔的前景，可以实现在成本可控的前提下有效提高OBS的数据质量，但目前，该技术还尚未成熟。实验研究显示，给地震传感器增加外部保护壳几乎可以完全屏蔽水流对地震传感器的影响，是一种低成本且比较有效的降噪技术。
另一方面，可以利用数据处理方法去除OBS数据的噪声，提高地震信号的信噪比，但目前这种方法仅限于去除OBS垂直分量上的噪声。OBS的低频噪声有两个主要来源：一个是底流作用导致的倾斜噪声（tilt noise），底流会导致OBS水平分量产生低频（＜0.1 Hz）噪声，当传感器尚未完美调平时，水流产生的噪声会从水平分量“泄漏”到垂直分量；另一个是次重力波（infragravity waves）对海底作用产生的屈服噪声（compliance noise），当海洋表面水波的波长超过水深时，长周期的水波可以穿透水层造成海底的垂直运动进而导致OBS数据上的屈服噪声。因此屈服噪声的高频极限与水深相关，可以通过一个简单的公式$f{\text{≈}} \sqrt{{g}/ （ {2\pi d} ） }$来近似计算（Bell et al，2015）。目前普遍使用的去除OBS噪声的数据处理方法也是针对这两种噪声，即分别利用水平、垂直噪声数据和压力、垂直噪声数据来计算水平到垂直分量以及压力到垂直分量的传递函数，并分别去除垂直分量上的倾斜噪声与屈服噪声，从而提高垂直分量上地震信号的信噪比，得到更高质量的地震波形。这种方法已被广泛应用于最近的科学实验中（Bell et al，2016；Bowden et al，2016；Lin et al，2016；Agius et al，2017；Tian，Ritzwoller，2017；Cai et al，2018；Doran，Laske，2019；Janiszewski et al，2019）。An等（2021）的研究表明，OBS水平分量的低频（＜0.1 Hz）噪声存在一个几乎不变的主导方向，建议通过旋转水平分量数据来有效降低垂直主导方向上的噪声水平。
基于上述问题，本文拟利用2019—2020年磐鲲OBS南海测试的数据，采用Bell等（2015）的方法去除海底地震仪垂直分量上的倾斜噪声，以提高垂直分量上地震波形的信噪比，为海洋地下深部结构反演等研究提供可靠数据基础。
1. 数据和方法
1.1 数据
本文使用2019—2020年磐鲲OBS在南海测试（图1）的数据。磐鲲海底地震仪（OBS）是南方科技大学自主研发的分体式宽频带OBS，关于2019—2020年南海测试和磐鲲OBS更详细的介绍可参考刘丹等（2022）、Wang等（2022）和Liu 等（2023）的文章。在用于测试的6台仪器（5台返回了有效数据，另1台也成功回收）中，只有台站K02产生了明显的倾斜噪声，因此本文的分析是基于台站K02在南海深水区（深度为3 896 m）自2019年10月至2020年5月期间的连续7个月的观测数据。
1.2 方法
在海底底流的作用下，OBS会产生水平方向上的噪声，而如果仪器倾斜，噪声会从水平分量“泄漏”到垂直分量。为了去除倾斜噪声，需要计算水平与垂直分量之间的传递函数。传递函数作为频率的函数显示了从一个分量传递到另一个分量的相移和比例系数。如果传递函数的相移为0，比例系数为常数且与频率不相关，那么去除倾斜噪声就可以简化为旋转地震计的坐标系，即将垂直分量旋转到重力的方向。然而，传递函数存在频率相关性，因此，使用传递函数去除倾斜噪声的方法更加可靠。另外，OBS在海底的倾斜特征可能会随着时间发生变化，为了使去除倾斜噪声的效果最优，需要至少每天计算一个传递函数（Bell et al，2015）。
本文研究使用Bell等（2015）的方法来计算传递函数。这一方法来源于Bendat和Piersol （1986）的理论。假设$ x （ t ） $和$ y （ t ） $均为具有正态分布误差的平稳随机过程（即地震仪记录中的理想噪声数据），将$ x （ t ） $和$ y （ t ） $分解为$ {n} $个子过程，则互谱密度函数${G}_{ xy} （ f ） \mathrm{可定}\mathrm{义}$为
$$ {G}_{ xy} （ f ） ＝\frac{1}{{n}}\sum\limits _{i＝1}^{{n}}{X}_{i}^{*} （ f ） {Y}_{i} （ f ） ，$$ (1) 式中，$ f $为谱频率，$ X （ f ） $和$ Y （ f ） $分别为$ x （ t ） $和$ y （ t ） $的傅里叶变换，*表示共轭。将X和Y分解为实部和虚部，分别定义实部和虚部的互谱密度函数$ {C}_{xy} （ f ） $和$ {Q}_{ xy} （ f ） $为：
$$ {C}_{xy} （ f ） ＝\frac{1}{{n}}\sum\limits _{i＝1}^{{n}} （ {X}_{{\rm{R}}}{Y}_{{\rm{R}}}＋{X}_{{\rm{I}}}{Y}_{{\rm{I}}} ） ，$$ (2) $$ {Q}_{xy} （ f ） ＝\frac{1}{{n}}\sum\limits _{i＝1}^{{n}} （ {X}_{{\rm{R}}}{Y}_{{\rm{I}}}－{X}_{{\rm{I}}}{Y}_{{\rm{R}}} ） {\text{.}} $$ (3) 定义$ x （ t ） $相对于自身的谱密度函数为$ {G}_{xx} $，则X到Y的最优线性传递函数$ {H}_{xy} （ f ） $为：
$$ {H}_{xy} （ f ） ＝\frac{{G}_{xy} （ f ） }{{G}_{xx} （ f ） } {\text{.}} $$ (4) 结合上述内容可以得到描述复传递函数的三个参数，即：增益系数A，相关性γ2和相位ϕ。其中，增益系数表示传递函数的振幅；相关性给出了$ x （ t ） $在指定频率与$ y （ t ） $匹配的程度，它是相关系数的平方，表示可以从$ x （ t ） $线性预测出$ y （ t ） $的幂的部分；相位则表示$ x （ t ） $与$ y （ t ） $之间的相移。描述复传递函数的增益系数、相关性和相位三个参数可表示如下：
$$ {A}_{xy} （ f ） ＝\frac{\left|{G}_{xy} （ f ） \right|}{{G}_{xx} （ f ） } ，$$ (5) $$ {\gamma }_{xy}^{2} （ f ） ＝\frac{{\left|{G}_{xy} （ f ） \right|}^{2}}{{G}_{xx} （ f ） {G}_{yy} （ f ） } ，$$ (6) $$ {\phi }_{xy} （ f ） ＝\mathrm{arctan}\frac{{Q}_{xy} （ f ） }{{C}_{xy} （ f ） } {\text{.}} $$ (7) 利用传递函数，可以计算出从水平分量“泄漏”到垂直分量的噪声，从而自垂直分量中去除这部分噪声。在早期的噪声去除方法（Crawford，Webb，2000）中，先后对两个水平分量各自“泄漏”到垂直分量的噪声进行去除。Bell等（2015）改进了这一方法，首先旋转水平分量以找到与垂直分量具有最大相关性的水平方向，接着利用旋转的水平分量一次性去除垂直分量的倾斜噪声。这种方法在最大相关性的方向上给出了地震计倾斜的物理方向，而倾斜角度可以通过计算倾斜频段上增益系数的反正切得到，并且最大化相关性可以将传递函数中的不确定性降至最低。
2. 处理过程
传递函数量化了两个分量在频域上的增益系数、相关性和相位关系。传递函数的计算依赖于信号平稳的统计特性，因此在计算传递函数前需要剔除地震等瞬变信号。通过在频域上对多个噪声时间窗中计算的增益系数、相关性和相位关系进行平均，可以得到合适的传递函数。选择噪声数据的方法包括：使用较长的时间窗口，使地震等瞬变信号占据较短的时间窗长度（Yang et al，2012）；利用地震目录选择相对没有地震发生的时间窗口（Bell et al，2015）；手动挑选时间窗口（Bell et al，2015）；利用程序自动剔除含有瞬变信号的窗口（Janiszewski et al，2019）等。本文采用Janiszewski等（2019）的方法选择噪声数据，因为地震目录并不能完全描述OBS记录中的瞬变信号，例如一些较小的近震不会出现在地震目录上，但它们会对平稳的噪声信号造成干扰。
选择噪声数据：首先需去除原始数据的仪器响应，设置时间窗口长度为7 200 s，相邻时间窗重叠30%，基于每天的噪声数据设置阈值以去除能量相对较高的噪声窗口（图2）；然后，利用Bell等（2015）的方法对水平分量（1和2分量）进行旋转并计算与垂直分量（Z分量）的相关性，具有最大相关性的水平分量（H分量）的方向即为地震计倾斜的方向。对每天的数据重复此处理过程，结果如图3所示，可以看出，台站K02地震计的倾斜方向在布设期间一直很稳定，大约在水平分量1顺时针旋转27°的方向。
![]()
图 2 2020年3月18日垂直（a）和水平（b，c）分量的频谱及噪声时间窗选择（d）图（d）中的方块表示窗口被接受，而方块的缺失（箭头所示）表示窗口未通过基于日平均值设置的阈值Figure 2. Spectrogram for the vertical （a） and horizontal （b，c） components，and the noise time windows selection （d） on 18th March，2020A square in fig. （d） indicates that the window is accepted， while the absence of a square （indicated by the arrow） indicates that the window does not pass the threshold set based on the daily average![]()
图 3 台站K02布设期间水平与垂直分量的最大相关性方向（a）和对应的相关大小（b）图（a）中的倾斜方向为相对水平分量1顺时针转过的角度Figure 3. Maximum coherence direction （a） and corresponding coherence value （b） of the horizontal to vertical component during station K02 deploymentThe tilt direction in fig.（a） represents the degrees turned clockwise relative to the horizontal component 1OBS在布设期间的倾斜特征可能会持续变化，对每天的噪声数据计算传递函数。利用图3中得到的H分量计算水平到垂直分量传递函数的增益系数、相关性和相位，过程示例如图4所示。为了使处理过程更加严谨，考虑到地震计的频带宽度（100 Hz—120 s），暂不考虑频率低于0.008 Hz的数据。在相关性最大的倾斜频段（0.008—0.01 Hz）上，增益系数基本恒定，相移近似为0，与倾斜效应一致。计算倾斜频段上增益系数的反正切得到倾斜角度，对每天的数据重复此过程即可得到布设期间台站的倾斜特征（图5）。可以看出，布设初期台站的倾斜角度出现一些小的变化，可能反映了OBS在海底软质沉积层中的不均匀沉降。在随后的绝大部分时间里台站K02的倾斜角度一直稳定在1°左右，对应水平与垂直分量的相关性在大部分时间里均较高，说明计算结果相对可靠。
![]()
图 4 水平到垂直分量传递函数的增益系数（a）、相关性（b）和相位（c）Figure 4. Admittance （a），coherence （b） and phase （c） of the horizontal to vertical component transfer function![]()
图 5 台站K02布设期间的倾斜角度（a）和其对应的相关性（b）Figure 5. Tilt angle （a） and corresponding coherence value （b） during station K02 deployment3. 结果
Peterson （1993）根据全球地震台网75个宽频带陆地地震台站的地震数据，挑选出了近
2000 条相对平静的噪声记录并分析了它们的噪声功率谱密度（power spectral density，缩写为PSD），给出了全球陆地台站的新高噪声模型（new high noise model，缩写为NHNM）和新低噪声模型（new low noise model，缩写为NLNM），如今这一模型已成为衡量地震台站噪声水平的重要依据。将传递函数应用于对应的水平分量地震数据，并将传递的那一部分噪声从垂直数据中去除，得到倾斜校正的地震波形。某两个地震波形去除倾斜噪声前后的对比，如图6所示。
![]()
图 6 去除倾斜噪声前后波形（滤波至0.01—0.1 Hz）及其功率谱密度 （PSD） 对比（a） 2020年3月18日巴厘岛南部附近发生的MW6.2地震；（b） 2019年12月25日东太平洋海隆南部附近发生的MW6.1地震Figure 6. Comparison of waveforms （filtered at 0.01−0.1 Hz） and its PSDs before and after tilt noise removed（a） MW6.2 earthquake occurred near the south of Bali on March 18，2020；（b） MW6.1 earthquake occurred near the southern East Pacific Rise on December 25，2019图6a所示为台站K02记录的一个发生在巴厘岛附近的MW6.2近震，在原始数据中就可以看到较为清晰的地震波形，但仍受到明显的长周期噪声的干扰。经过倾斜噪声校正后，长周期噪声的干扰几乎被完全去除了，在地震发生前的相对安静时期波形振幅在0附近无明显波动，去除倾斜噪声后的地震波形几乎达到了与附近陆地台站相当的波形质量（图7）。图6b是台站K02记录的一个发生在东太平洋海隆附近的MW6.1较远地震，它的振幅比几乎同样震级的近震（图6a）小了一个数量级，在原始波形中几乎看不到地震信号，但进行倾斜噪声校正后，各震相变得较为清晰，这表明去除倾斜噪声可以提高OBS对远震的探测能力。同时，从去除倾斜噪声前后地震数据的功率谱密度曲线可以看出，倾斜噪声主要存在于低频段（＜0.03 Hz），去除倾斜噪声后，地震波形的PSD曲线在低频段均有不同程度的降低，表明去除倾斜噪声可以提高OBS记录的地震波形在低频段的信噪比。
![]()
图 7 巴厘岛附近MW6.2地震和台站位置（a）及去除倾斜噪声后垂直分量波形与附近陆地台站对比（b）（滤波至0.01—0.1 Hz）Figure 7. The MW6.2 earthquake near Bali and station locations （a） and comparison of vertical component waveforms after tilt noise removed with nearby land stations （b）（filtered at 0.01−0.1 Hz）由于海底地震台站比陆地台站稀疏得多，远震体波在海洋岩石圈深部结构成像中的作用受到较大限制，而地震面波横跨海盆传播，因此地震面波尤其是低频面波对于洋壳和海洋岩石圈的探测至关重要。我们对台站K02接收到的部分地震做了群速度的频率−时间分析得到其频散曲线，图8给出了一个地震波形（图6b中的地震）去除倾斜噪声前后的瑞雷波频散特征。可以看出，在去除倾斜噪声前，频散曲线（图8a）受到噪声影响无法反映海洋岩石圈的结构变化特征；而去除倾斜噪声后，典型的、平滑的频散曲线（图8b）呈现出面波群速度随岩石圈及上地幔结构在深度方向上的变化规律；所以，去除倾斜噪声可以提高OBS记录的低频面波数据质量，使其更有利于进行海洋岩石圈的速度结构成像。
![]()
图 8 原始地震波形（a）与去除倾斜噪声后的地震波形（b）的频散特征Figure 8. Dispersion characteristics of the raw seismic waveform （a） and the seismic waveform after tilt noise removed （b）4. 讨论与结论
4.1 地震计倾斜角度对倾斜噪声的影响
地震计的倾斜角度对OBS垂直分量数据的倾斜噪声有着很大影响。台站K03与台站K02具有相似的倾斜特征（图9），但台站K03的倾斜角度（约0.5°）要小于台站K02 （约1.0°）。如图10所示，台站K03的数据并未表现出明显的倾斜噪声。另一台站K08，其水平分量和垂直分量噪声的相关性一直小于0.4 （刘丹等，2022），表明这一台站的地震计没有倾斜，其数据也未受到倾斜噪声的影响。由于这几台OBS布放在同一海区的相近位置（图1），受到底流等因素的影响相似，因此它们的倾斜噪声水平取决于地震计的倾斜角度，当倾斜角度约小于0.5°时，其倾斜噪声水平就很微弱，可以忽略不计。
![]()
图 9 台站K03布设期间的倾斜角度（a）及对应的相关性大小（b）Figure 9. Tilt angle （a） and corresponding coherence value （b） during station K03 deployment![]()
图 10 巴厘岛附近MW6.2地震去除倾斜噪声前后波形及其功率谱密度（PSD）对比Figure 10. Comparison of the waveforms and its PSDs before and after tilt noise removed of the MW6.2 earthquake occurred near Bali上述分析表明OBS的调平精度严重影响了仪器的数据质量。在OBS布放初期地震计调平并锁紧后，由于海底的复杂环境OBS很可能会再次倾斜，例如，地震计在海底沉积物中逐渐发生不均匀沉降。因此，对于实际观测的OBS来说，除了精准的调平系统以外，还需要定期监测地震计的倾斜角度，当倾斜角度超过设定阈值时重新调平是降低OBS数据倾斜噪声的有效方案（刘丹等，2022；Liu et al，2023）。
4.2 海底底流对倾斜噪声的影响
屏蔽OBS的地震计模块，使得底流不直接作用于地震计，可以大大降低OBS数据的噪声水平（Bell et al，2015），最近作者的OBS水槽实验也证实了这一点。而南科大的磐鲲正是这种将地震计屏蔽在其内部的OBS，磐鲲OBS的地震计不应受到底流的影响。然而，通过对台站K02的分析，磐鲲OBS还是产生了较为明显的倾斜噪声。
这一现象产生的原因可能是由于底流与实际海底之间相互作用并不像刚性基底水槽中那样简单。实验室水槽中的水流厚度较低，且水流速度在各深度基本相同，而实际海底会形成厚度为几米到几十米的底流边界湍流层。因此，海底底流中的地震计会记录两个分量，一个是水流直接作用在地震计上绕流产生的波动，另一个是地震计附近更大范围内湍流产生的波动。越接近海底，底流流速越低，因此OBS记录的因绕流产生的波动较小，但与海底边界层中的湍流涡旋相关的波动不会朝着海底逐渐减小，因此，海底地震仪记录到的压力波动可能占主要地位（Webb，1988）。对于磐鲲OBS来说，绕流产生的噪声基本被完全屏蔽了，但湍流涡旋在地震计周边产生的压力波动并没有降低，这种压力波动可能通过海底或者屏蔽空间内的水体传至地震计，产生了较明显的倾斜噪声。
此外，底流一直被认为是倾斜噪声的重要来源，但尚无观测数据证明底流是倾斜噪声的唯一来源。当仪器发生倾斜时，其它因素导致的水平噪声也有可能“泄漏”至垂直分量从而产生倾斜噪声，例如，次重力波也会产生水平分量的噪声（Webb，1998）。因此，磐鲲OBS的倾斜噪声来源也可能是底流之外的其它因素。
若要进一步研究倾斜噪声与海底底流的关系，还需要在OBS附近同时观测海底底流以及设置地震计有无屏蔽的对照实验。
4.3 结论
受多个因素影响，OBS的地震计在海底并不总是能够实现理想的调平（如Liu et al，2014；Hung et al，2019），在底流作用下，OBS数据中常常含有不同程度的倾斜噪声。利用旋转的水平分量（H分量）与垂直分量（Z分量）之间的传递函数可以得到OBS的倾斜特征，并且可以在很大程度上去除垂直分量上的倾斜噪声，提高OBS的数据质量，尤其可以提高地震波形在低频段的信噪比。
对于原始记录中震相较为清晰的地震事件，去除倾斜噪声可以进一步提高地震波形的信噪比，使其达到与附近陆地台站相当的波形质量。而对于在原始记录中震相不够清晰的地震事件，去除倾斜噪声可以大幅提高地震波形的质量，使得地震波形震相较为清晰，其瑞雷波的频散特征更有利于进行深部岩石圈和上地幔的速度结构成像。
尽管去除倾斜噪声可以提高OBS记录的垂直分量地震波形的信噪比，但是对于水平分量上的噪声却无能为力。如何提高OBS水平分量数据的信噪比，是一个更大的挑战。这需要通过OBS与海底更好地耦合，降低底流对OBS的影响等技术问题的解决来实现。未来OBS数据质量的进一步提高离不开技术的进步，而高质量的OBS数据也会持续推动地震学和海洋地球物理学的进一步发展。
-
![]()
图 3 松辽盆地沉积层厚度（a）、地壳厚度（b）、地壳平均vP/vS （c）以及地壳伸展系数（d）分布图
图（b）中，黑色方块为采用H-β搜索方法获得有效结果的台站，三角形为H-β搜索结果无效的台站
Figure 3. Map of sediment thickness （a），crustal thickness （b），crustal vP/vS ratio （c） and stretching factor （d） beneath the Songliao basin
The black squares in Fig. （b） represent the seismic stations with the valid H-β values，and triangles represent those without valid H-β values
![]()
图 1 东北地区地质构造和地震台站分布示意图
黄色实心圆为满洲里—绥芬河线性流动台阵（SM台阵），蓝色实心圆为NECESSArray流动台阵
Figure 1. Map showing major geological tectonic settings and portable seismic arrays deployed in Northeast China
Yellow and blue solid circles represent two portable seismic arrays，the Suifenhe to Manzhouli （SM） and NECESSArray，respectively. The purple solid line outlines the Songliao basin
![]()
图 2 远震地震事件分布
红色三角形为研究区域中心，蓝色圆点为地震事件的震中分布
Figure 2. Epicentral distribution of teleseismic events used in this study
Red triangle represents the center of the studied area，blue dots denote seismic events
![]()
图 4 地壳厚度与沉积层厚度（a）和壳内平均波速比 vP/vS （b）的相关性
Figure 4. Relationship between crustal thicknesses and sediment thicknesses （a） and vP/vS ratios （b）
表 1 本文H-β网格搜索法中所用模型
Table 1 Model used in H-β grid search method
模型分层 vP/（km·s−1） ρ/（10 3 kg·m−3） 厚度H/km vS/（km·s−1） 沉积层 2.1 1.97 0—7.0 0—2.5 地壳 6.4 2.7 25—40 3.0—5.0 地幔 8.0 3.3 ∞ 4.5 
下载: 导出CSV
-
包汉勇,郭战峰,张罗磊,黄亚平. 2013. 盆地伸展系数求取方法与评价:以苏北盆地为例[J]. 石油实验地质,35(3):331–338. doi: 10.11781/sysydz201303331 Bao H Y,Guo Z F,Zhang L L,Huang Y P. 2013. Calculating methods and assessment of stretching factor:A case study of northern Jiangsu basin[J]. Petroleum Geology &Experiment,35(3):331–338 (in Chinese).
窦立荣. 1992. 东北含油气区中新生代断陷盆地石油地质特征及资源远景[J]. 江汉石油学院学报,14(1):1–8. Dou L R. 1992. Characteristics of petroleum geology and petroleum potential in Meso-Cenozoic fault basins in northeastern petroleum-bearing provinces,China[J]. Journal of Jianghan Petroleum Institute,14(1):1–8 (in Chinese).
符伟. 2019. 深反射剖面揭示的松辽盆地北部深部结构、动力学背景与油气远景[D]. 长春: 吉林大学: 67–82. Fu W. 2019. Deep Structure, Dynamic Background and Hydrocarbon Prospect of Northern Songliao Basin: Revealed by Deep Seismic Reflection Profile[D]. Changchun: Jilin University: 67–82 （in Chinese）.
傅维洲,杨宝俊,刘财,Кгылов С B. 1998. 中国满洲里—绥芬河地学断面地震学研究[J]. 长春科技大学学报,28(2):87–93. Fu W Z,Yang B J,Liu C,Кгылов C B. 1998. Study on the seismology in Manzhouli-Suifenhe geoscience transect of China[J]. Journal of Changchun University of Science and Technology,28(2):87–93 (in Chinese).
高立新,戴勇. 2020. 现今中国东北地区地球动力学环境[J]. 大地测量与地球动力学,40(11):1101–1107. Gao L X,Dai Y. 2020. The present geodynamic environment of Northeast China[J]. Journal of Geodesy and Geodynamics,40(11):1101–1107 (in Chinese).
高延光,李永华. 2014. 中国东北-华北地区地壳厚度与泊松比及其地质意义[J]. 地球物理学报,57(3):847–857. doi: 10.6038/cjg20140314 Gao Y G,Li Y H. 2014. Crustal thickness and vP/vS in the Northeast China-North China region and its geological implication[J]. Chinese Journal of Geophysics,57(3):847–857 (in Chinese).
高占永. 2015. 中国东北地区地壳上地幔结构的接收函数研究[D]. 北京: 中国地震局地球物理研究所: 31–46. Gao Z Y. 2015. The Study of the Crustal and Upper Mantle Structure in Northeast China From Teleseismic Receiver Function[D]. Beijing: Institute of Geophysics, China Earthquake Administration: 31–46 （in Chinese）.
葛荣峰,张庆龙,王良书,解国爱,徐士银,陈娟,王锡勇. 2010. 松辽盆地构造演化与中国东部构造体制转换[J]. 地质论评,56(2):180–195. Ge R F,Zhang Q L,Wang L S,Xie G A,Xu S Y,Chen J,Wang X Y. 2010. Tectonic evolution of Songliao basin and the prominent tectonic regime transition in eastern China[J]. Geological Review,56(2):180–195 (in Chinese).
侯贺晟,王成善,张交东,马峰,符伟,王璞珺,黄永建,邹长春,高有峰,高远,张来明,杨瑨,国瑞. 2018. 松辽盆地大陆深部科学钻探地球科学研究进展[J]. 中国地质,45(4):641–657. doi: 10.12029/gc20180401 Hou H S,Wang C S,Zhang J D,Ma F,Fu W,Wang P J,Huang Y J,Zou C C,Gao Y F,Gao Y,Zhang L M,Yang J,Guo R. 2018. Deep continental scientific drilling engineering in Songliao basin:Progress in earth science research[J]. Geology in China,45(4):641–657 (in Chinese).
胡望水,吕炳全,张文军,毛治国,冷军,官大勇. 2005. 松辽盆地构造演化及成盆动力学探讨[J]. 地质科学,40(1):16–31. Hu W S,Lü B Q,Zhang W J,Mao Z G,Leng J,Guan D Y. 2005. An approach to tectonic evolution and dynamics of the Songliao basin[J]. Chinese Journal of Geology,40(1):16–31 (in Chinese).
嵇少丞,王茜,杨文采. 2009. 华北克拉通泊松比与地壳厚度的关系及其大地构造意义[J]. 地质学报,83(3):324–330. doi: 10.3321/j.issn:0001-5717.2009.03.002 Ji S C,Wang Q,Yang W C. 2009. Correlation between crustal thickness and Poisson’s ratio in the North China craton and its implication for lithospheric thinning[J]. Acta Geologica Sinica,83(3):324–330 (in Chinese).
李恩泽,刘财,张良怀,曾昭发. 2012. 松辽盆地地震构造与地震活动相关性研究[J]. 地球物理学进展,27(4):1337–1349. doi: 10.6038/j.issn.1004-2903.2012.04.007 Li E Z,Liu C,Zhang L H,Zeng Z F. 2012. The correlation of structure and earthquake in Songliao basin[J]. Progress in Geophysics,27(4):1337–1349 (in Chinese).
李国良. 2016. 瑞雷波椭圆率的测定与在反演S波速度结构中的应用[D]. 北京: 中国石油大学（北京）: 14–28. Li G L. 2016. Measurement of Rayleigh Wave Ellipticity and Its Application to the Joint Inversion of High-Resolution S-Wave Velocity Structure[D]. Beijing: China University of Petroleum （Beijing）: 14–28 （in Chinese）.
李国良. 2019. 利用被动源数据联合反演盆地3D速度结构[D]. 北京: 中国石油大学（北京）: 47–69. Li G L. 2019. Joint Inversion of Basin-Wide 3D Sedimentary Structure With Passive Seismic Data[D]. Beijing: China University of Petroleum （Beijing）: 47–69 （in Chinese）.
李英宾,李毅,魏滨,刘波,张占彬,杨明. 2019. CSAMT和浅层地震在松辽盆地西南部铀矿勘查中的应用[J]. 地质与勘探,55(6):1442–1451. Li Y B,Li Y,Wei B,Liu B,Zhang Z B,Yang M. 2019. Application of CSAMT and shallow seismic reflection to uranium exploration in southwestern Songliao basin[J]. Geology and Exploration,55(6):1442–1451 (in Chinese).
林畅松,张燕梅. 1995. 拉伸盆地模拟理论基础与新进展[J]. 地学前缘,2(3/4):79–88. Lin C S,Zhang Y M. 1995. Quantitative stretching models and computer simulation of rift basin[J]. Earth Science Frontiers,2(3/4):79–88 (in Chinese).
刘德来,陈发景,关德范,唐建人,刘翠荣. 1996. 松辽盆地形成、发展与岩石圈动力学[J]. 地质科学,31(4):397–408. Liu D L,Chen F J,Guan D F,Tang J R,Liu C R. 1996. A study on lithospheric dynamics of the origin and evolution in the Songliao basin[J]. Chinese Journal of Geology,31(4):397–408 (in Chinese).
马海超,储日升,盛敏汉,危自根. 2020. 利用深源近震高频Ps转换波震相研究松辽盆地沉积层结构[J]. 大地测量与地球动力学,40(2):214–220. Ma H C,Chu R S,Sheng M H,Wei Z G. 2020. Sedimentary structures of the Songliao basin using high-frequency Ps converted wave from local deep earthquakes[J]. Journal of Geodesy and Geodynamics,40(2):214–220 (in Chinese).
王璞珺,刘海波,任延广,万晓樵,王树学,瞿雪姣,蒙启安,黄永建,黄清华,高有峰,王成善. 2017. 松辽盆地白垩系大陆科学钻探“松科2井”选址[J]. 地学前缘,24(1):216–228. Wang P J,Liu H B,Ren Y G,Wang X Q,Wang S X,Qu X J,Meng Q A,Huang Y J,Huang Q H,Gao Y F,Wang C S. 2017. How to choose a right drilling site for the ICDP Cretaceous Continental Scientific Drilling in the Songliao basin (SK2),Northeast China[J]. Earth Science Frontiers,24(1):216–228 (in Chinese).
王仁涛,李志伟,包丰,谢军,赵建忠. 2019. 松辽盆地沉积层结构的短周期地震背景噪声成像研究[J]. 地球物理学报,62(9):3385–3399. doi: 10.6038/cjg2019M0144 Wang R T,Li Z W,Bao F,Xie J,Zhao J Z. 2019. S-wave velocity structure of sediment in Songliao basin from short-period ambient noise tomography[J]. Chinese Journal of Geophysics,62(9):3385–3399 (in Chinese).
危自根,陈凌. 2012. 东北地区至华北北缘地壳结构的区域差异:地壳厚度与波速比的联合约束[J]. 地球物理学报,55(11):3601–3614. doi: 10.6038/j.issn.0001-5733.2012.11.009 Wei Z G,Chen L. 2012. Regional differences in crustal structure beneath northeastern China and northern North China Craton:Constraints from crustal thickness and vP/vS ratio[J]. Chinese Journal of Geophysics,55(11):3601–3614 (in Chinese).
危自根,储日升,陈凌. 2015. 华北克拉通地壳结构区域差异的接收函数研究[J]. 中国科学:地球科学,45(10):1504–1514. Wei Z G,Chu R S,Chen L. 2015. Regional differences in crustal structure of the North China Craton from receiver functions[J]. Science China Earth Sciences,58(12):2200–2210. doi: 10.1007/s11430-015-5162-y
武岩. 2011. 利用接收函数方法研究华北克拉通地壳上地幔结构[D]. 北京: 中国地震局地球物理研究所: 25–39. Wu Y. 2011. The Structure of the Crust and Upper Mantle in North China Craton From Teleseismic Receiver Function[D]. Beijing: Institute of Geophysics, China Earthquake Administration: 25–39 （in Chinese）.
谢振新,吴庆举,周仕勇,朱敏. 2018. 兴蒙造山带诺敏河火山群地壳厚度与波速比研究[J]. 地球物理学报,61(12):4805–4816. doi: 10.6038/cjg2018M0197 Xie Z X,Wu Q J,Zhou S Y,Zhu M. 2018. Study of crustal thickness and vP/vS ratio beneath the Nuomin River volcanoes[J]. Chinese Journal of Geophysics,61(12):4805–4816 (in Chinese).
杨宝俊,穆石敏,金旭,刘财. 1996. 中国满洲里—绥芬河地学断面地球物理综合研究[J]. 地球物理学报,39(6):772–782. doi: 10.3321/j.issn:0001-5733.1996.06.007 Yang B J,Mu S M,Jin X,Liu C. 1996. Synthesized study on the geophysics of Manzhouli-Suifenhe geoscience transect,China[J]. Chinese Journal of Geophysics,39(6):772–782 (in Chinese).
姚志祥,王椿镛,曾融生,楼海,周民都. 2014. 利用接收函数方法研究西秦岭构造带及其邻区地壳结构[J]. 地震学报,36(1):1–19. doi: 10.3969/j.issn.0253-3782.2014.01.001 Yao Z X,Wang C Y,Zeng R S,Lou H,Zhou M D. 2014. Crustal structure in western Qinling tectonic belt and its adjacent regions deduced from receiver functions[J]. Acta Seismologica Sinica,36(1):1–19 (in Chinese).
余嘉顺,曹俊兴,鲍新毅,黄跃. 2003. 表面低速层对勘探地震横波波形影响的模拟研究[J]. 成都理工大学学报(自然科学版),30(6):583–587. doi: 10.3969/j.issn.1671-9727.2003.06.006 Yu J S,Cao J X,Bao X Y,Huang Y. 2003. A modeling of the effects of ground surface weathering layers on reflected shear waves in seismic exploration[J]. Journal of Chengdu University of Technology (Science &Technology Edition)
,30(6):583–587 (in Chinese). 张广成,吴庆举,潘佳铁,张风雪,余大新. 2013. 利用H-K叠加方法和CCP叠加方法研究中国东北地区地壳结构与泊松比[J]. 地球物理学报,56(12):4084–4094. doi: 10.6038/cjg20131213 Zhang G C,Wu Q J,Pan J T,Zhang F X,Yu D X. 2013. Study of crustal structure and Poisson ratio of NE China by H-K stack and CCP stack methods[J]. Chinese Journal of Geophysics,56(12):4084–4094 (in Chinese).
张毅. 2019. 应用接收函数方法研究中国东部壳幔间断面结构[D]. 北京: 中国地质大学（北京）: 29–41. Zhang Y. 2019. Study of the Crust-Mantle Discontinuity Structure in Eastern China With Receiver Function Method[D]. Beijing: China University of Geosciences （Beijing）: 29–41 （in Chinese）.
周庆华,冯子辉,门广田. 2007. 松辽盆地北部徐家围子断陷现今地温特征及其与天然气生成关系研究[J]. 中国科学:D辑,37(增刊2):177–188. Zhou Q H,Feng Z H,Men G T. 2008. Present geotemperature and its suggestion to natural gas generation in Xujiaweizi fault-depression of the northern Songliao basin[J]. Science in China:Series D,51(1):207–220.
朱洪翔,田有,刘财,冯晅,杨宝俊,刘才华,刘廷,马锦程. 2017. 中国东北地区高分辨率地壳结构:远震接收函数[J]. 地球物理学报,60(5):1676–1689. doi: 10.6038/cjg20170506 Zhu H X,Tian Y,Liu C,Feng X,Yang B J,Liu C H,Liu T,Ma J C. 2017. High-resolution crustal structure of Northeast China revealed by teleseismic receiver functions[J]. Chinese Journal of Geophysics,60(5):1676–1689 (in Chinese).
朱洪翔,田有,刘财,冯晅. 2018. 沉积盆地地区地壳结构估计:预测反褶积方法消除接收函数多次波混响[J]. 地球物理学报,61(9):3664–3675. doi: 10.6038/cjg2018L0152 Zhu H X,Tian Y,Liu C,Feng X. 2018. Estimation of the crustal structure beneath the sedimentary basin:Predictive deconvolution method to remove multiples reverberations of the receiver function[J]. Chinese Journal of Geophysics,61(9):3664–3675 (in Chinese).
Bao Y F,Niu F L. 2017. Constraining sedimentary structure using frequency-dependent P wave particle motion:A case study of the Songliao basin in NE China[J]. J Geophys Res,122(11):9083–9094. doi: 10.1002/2017JB014721
Christensen N I,Mooney W D. 1995. Seismic velocity structure and composition of the continental crust:A global view[J]. J Geophys Res,100(B6):9761–9788. doi: 10.1029/95JB00259
Gilbert F,Backus G E. 1966. Propagator matrices in elastic wave and vibration problems[J]. Geophysics,31(2):326–332. doi: 10.1190/1.1439771
Graves R W,Pitarka A,Somerville P G. 1998. Ground-motion amplification in the Santa Monica area:Effects of shallow basin-edge structure[J]. Bull Seismol Soc Am,88(5):1224–1242.
Guo Z,Chen Y J,Ning J Y,Feng Y G,Grand S P,Niu F L,Kawakatsu H,Tanaka S,Obayashi M,Ni J. 2015. High resolution 3-D crustal structure beneath NE China from joint inversion of ambient noise and receiver functions using NECESSArray data[J]. Earth Planet Sci Lett,416:1–11. doi: 10.1016/j.jpgl.2015.01.044
Haskell N A. 1953. The dispersion of surface waves on multilayered media[J]. Bull Seismol Soc Am,43(1):17–34. doi: 10.1785/BSSA0430010017
He J,Wu Q J,Sandvol E,Ni J,Gallegos A,Gao M T,Ulziibat M,Demberel S. 2016. The crustal structure of south central Mongolia using receiver functions[J]. Tectonics,35(6):1392–1403. doi: 10.1002/2015TC004027
Kennett B L N,Kerry N J,Woodhouse J H. 1978. Symmetries in the reflection and transmission of elastic waves[J]. Geophys J Int,52(2):215–229. doi: 10.1111/j.1365-246X.1978.tb04230.x
Leahy G M,Saltzer R L,Schmedes J. 2012. Imaging the shallow crust with teleseismic receiver functions[J]. Geophys J Int,191(2):627–636. doi: 10.1111/j.1365-246X.2012.05615.x
Li G L,Chen H C,Niu F L,Guo Z,Yang Y J,Xie J. 2016. Measurement of Rayleigh wave ellipticity and its application to the joint inversion of high-resolution S wave velocity structure beneath Northeast China[J]. J Geophys Res,121(2):864–880. doi: 10.1002/2015JB012459
Lin C S,Li S T,Zhang Q M. 1997. Lithospheric stretching,subsidence and thermal history modeling:Application to Yinggehai,Qiongdongnan and Songliao basins in East China[J]. J China Univ Geosci,8(1):83–89.
McKenzie D. 1978. Some remarks on the development of sedimentary basins[J]. Earth Planet Sci Lett,40(1):25–32. doi: 10.1016/0012-821X(78)90071-7
Ren J Y,Tamaki K,Li S T,Zhang J X. 2002. Late Mesozoic and Cenozoic rifting and its dynamic setting in eastern China and adjacent areas[J]. Tectonophysics,344(3/4):175–205.
Tao K,Liu T Z,Ning J Y,Niu F L. 2014a. Estimating sedimentary and crustal structure using wavefield continuation:Theory,techniques and applications[J]. Geophys J Int,197(1):443–457. doi: 10.1093/gji/ggt515
Tao K,Niu F L,Ning J Y,Chen Y J,Grand S,Kawakatsu H,Tanaka S,Obayashi M,Ni J. 2014b. Crustal structure beneath NE China imaged by NECESSArray receiver function data[J]. Earth Planet Sci Lett,398:48–57. doi: 10.1016/j.jpgl.2014.04.043
Wei H H,Liu J L,Meng Q R. 2010. Structural and sedimentary evolution of the southern Songliao basin,Northeast China,and implications for hydrocarbon prospectivity[J]. AAPG Bull,94(4):531–564.
Wernicke B. 1985. Uniform-sense normal simple shear of the continental lithosphere[J]. Can J Earth Sci,22(1):108–125. doi: 10.1139/e85-009
Xiong X S,Gao R,Li Y K,Hou H S,Liang H D,Li W H,Guo L H,Lu Z W. 2015. The lithosphere structure of the Great Xing’an Range in the eastern Central Asian Orogenic Belt:Constrains from the joint geophysical profiling[J]. J Asian Earth Sci,113:481–490. doi: 10.1016/j.jseaes.2015.06.006
Yeck W L,Sheehan A F,Schulte-Pelkum V. 2013. Sequential H‐κ stacking to obtain accurate crustal thicknesses beneath sedimentary basins[J]. Bull Seismol Soc Am,103(3):2142–2150. doi: 10.1785/0120120290
Yu Y Q,Song J G,Liu K H,Gao S S. 2015. Determining crustal structure beneath seismic stations overlying a low-velocity sedimentary layer using receiver functions[J]. J Geophys Res:Solid Earth,120(5):3208–3218. doi: 10.1002/2014JB011610
Zhang R Q,Wu Q J,Sun L,He J,Gao Z Y. 2014. Crustal and lithospheric structure of Northeast China from S-wave receiver functions[J]. Earth Planet Sci Lett,401:196–205. doi: 10.1016/j.jpgl.2014.06.017
Zheng C,Zhang R Q,Wu Q J,Li Y H,Zhang F X,Shi K X,Ding Z F. 2019. Variations in crustal and uppermost mantle structures across eastern Tibet and adjacent regions:Implications of crustal flow and asthenospheric upwelling combined for expansions of the Tibetan Plateau[J]. Tectonics,38(7):3167–3181.
Zhu L,Kanamori H. 2000. Moho depth variation in southern California from teleseismic receiver functions[J]. J Geophys Res:Solid Earth,105(B2):2969–2980.
-
期刊类型引用(1)
1. 刘丹，杨挺，吴越楚，王宜志，林俊国. 应用于宽频带海底地震仪的低频差分压力传感器的研制. 地球物理学报. 2025(04): 1367-1380 . 
百度学术
其他类型引用(0)
计量
- 文章访问数: 1073
- HTML全文浏览量: 403
- PDF下载量: 226
- 被引次数: 1
