2. 青岛海洋科学与技术国家实验室海洋矿产资源评价与探测技术功能实验室, 青岛 266061;
3. 中交天航(大连)浚航工程有限公司, 大连 116039
2. Evaluation and Detection Technology Laboratory of Marine Mineral Resources, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266061, China;
3. CCCC-TDC(DaLian) Channel Dredging Engineering Co., Ltd., Dalian 116039, China
在地震层析成像中,远偏移距折射初至对提高速度模型深部结构的成像精度具有重要作用(Mooney and Weaver, 1989;Operto and Charvis, 1996;Funck et al., 2008;王帅军等, 2014).然而,随着偏移距的增加,折射波的信噪比较低.这个问题在噪声强烈的海洋广角海底地震仪(OBS)观测中尤为突出(卫小冬等, 2011;吴振利等, 2011;阮爱国等, 2011;Alshuhail et al., 2012;Al-Shuhail, 2015;Liu et al., 2015;刘俊等, 2015;Zou et al., 2016).因此,有必要改善折射波的信噪比,以提高地震初至拾取的准确度.
折射波干涉法是一种较为有效的折射波增强方法,它通过对不同炮的地震道做互相关和叠加,达到压制噪声的效果(Dong et al., 2006).但该方法只能得到某个虚炮点激发、各检波点接收的数据,而无法得到真实炮点激发、各检波点接收的数据.为解决此问题,Bharadwaj等(2011)和Hanafy等(2011)提出了超虚折射干涉法(Supervirtual Refraction Interferometry,SRI),该方法在互相关构建虚拟道的基础上,将虚拟道与对应近道做卷积和叠加实现了对原始数据中折射波的增强.初步的应用显示,使用SRI方法增强后的折射波进行地震初至拾取时,所拾取初至的稳定性和精度显著提高(Mallinson et al., 2011;Al-Hagan et al., 2014),层析速度建模的准确度具有明显改善(Mallinson et al., 2011).为进一步提高折射波增强的稳定性和精度,前人在SRI方法基础上进行了一系列的改进,所解决的问题包括远近道叠加次数不均匀(乔宝平等, 2014;An et al., 2017a)、叠加道相似度低(吕雪梅等, 2018)、三维折射波的增强(Lu et al., 2014;An et al., 2017b)等.SRI方法还被拓展应用于深海OBS数据(Bharadwaj et al., 2013)、地核折射波数据(Bharadwaj et al., 2014)等,取得了良好效果.
然而,SRI方法效果依赖于叠加炮点数量,参与互相关产生虚拟道的炮点数量越少,虚拟道叠加次数越少、抗噪能力越差、准确度越低,最终的增强效果越差.在广角OBS探测中,由于OBS数量有限,导致长测线中的OBS间距通常在数公里甚至十几公里(Yan et al., 2001;吴振利等, 2011;阮爱国等, 2011;卫小冬等, 2011;Zou et al., 2016).这种情况下,同一炮所激发的折射波仅能被少数几个OBS观测到.在应用SRI方法时,通常基于互易性原理把台站与炮点互换,因此较大的台站间距导致可参与虚拟道叠加的台站数较少.同时,强烈的海洋噪声也使远偏移距折射波的信噪比较低.这些原因共同导致了SRI方法对广角OBS数据远偏移距折射波的增强能力不足.
针对此问题,本文提出基于相邻虚拟道叠加的超虚折射干涉法(Supervirtual Refraction Interferometry based on the Stacking of Neighbor Virtual-traces, SRI-SNV),通过引入相邻虚拟道叠加,以有效提高虚拟道叠加次数.本文将在介绍SRI-SNV方法基本原理的基础上进行理论实验和实际资料应用,并将SRI-SNV方法结果与SRI方法结果进行对比分析,以显示SRI-SNV方法对OBS数据折射波的增强效果.
1 基于相邻虚拟道叠加的超虚折射干涉法原理SRI方法的目标是重构高信噪比的折射波.该方法先通过互相关建立参考地震道(高信噪比地震道)与目标地震道(低信噪比地震道)之间的关系,然后在此关系基础上将参考地震道数据叠加以重构目标地震道(Hanafy et al., 2011;Bharadwaj et al., 2012).对于广角OBS观测系统,SRI方法实现折射波增强的基本步骤是:首先,在共台站集选择一个参考地震道(高信噪比近道)和一个目标地震道(低信噪比远道),将不同的台站集里这两道的折射波形做时窗截取,并通过互相关产生虚拟道;其次,叠加不同台站所对应的虚拟道;之后,将叠加虚拟道与对应的参考地震道做卷积,生成新的目标地震道;然后,选择一个新的参考地震道重新进行上述操作,生成新的目标地震道;最后,叠加所有的新目标地震道,从而实现一个目标地震道的增强.循环上述步骤即可实现对所有地震道上折射波的增强.具体流程如图 1所示.
上述步骤显示,SRI方法在增强折射波信噪比过程中有两次叠加操作,一次是虚拟道叠加,另一次是新目标地震道叠加.虚拟道叠加次数越多,降噪效果越好.由于SRI方法的虚拟道叠加次数依赖于台站数量,因此参与同一折射波增强的台站数量越多,对地震道信噪比的改善效果越好.但在实际OBS数据采集中台站分布稀疏,同一炮所激发折射波的有效接收台站数量较少,导致可参与虚拟道叠加的台站数较少,使虚拟道产生过程中的抑噪效果降低.因此,提高虚拟道叠加次数是广角OBS数据折射波增强的关键.
对于广角OBS数据,SRI方法选择共台站集上的两道进行互相关,从而产生虚拟道.由于不同台站所对应虚拟道的峰值应该在相同时刻,因此可以通过叠加来提高虚拟道的信噪比.若近道和远道的折射波前近似平行,选择近道的邻近道和远道的邻近道,并保持这两个邻近道具有与原来远近道相同的道间距差,那么这两个邻近道的折射波互相关所产生的虚拟道与SRI所产生虚拟道的峰值应该在相同位置,因此同样可以参与叠加以提高虚拟道信噪比.需要注意的是,由于实际地下界面存在起伏,上述假设只在小的空间范围内才能近似成立,所以只有与SRI所选择地震道相邻的地震道才能产生与SRI虚拟道足够相近的虚拟道,我们称之为相邻虚拟道.相邻虚拟道参与叠加后,虚拟道的叠加次数成倍增加,有助于提高增强后折射波的信噪比.该操作可以添加到SRI方法互相关产生虚拟道的步骤中(图 1),我们将新方法称为基于相邻虚拟道叠加的超虚折射干涉法(SRI-SNV).下面介绍SRI-SNV方法中相邻虚拟道叠加的原理.
图 2展示了SRI-SNV方法的相邻虚拟道叠加过程,它包含两个基本步骤:相邻地震道互相关生成相邻虚拟道和不同台站所对应的相邻虚拟道叠加.图 2a显示,在mi、mj点激发和n1点接收的两道做互相关后,生成了虚拟道ψ(mi, mj)n1.选取与mi和mj同侧相邻且具有相同偏移距的地震道重复上述互相关操作,就生成了相邻虚拟道.相邻虚拟道参与叠加后,在不增加台站数量的情况下,虚拟道叠加次数会成倍增加,如图 2b所示.所得到的叠加虚拟道记为Φ(mi, mj),此即对应于mi和mj的最终虚拟道.按照图 1中的流程,该虚拟道被用于和原始记录卷积来生成增强后的数据.下面介绍SRI-SNV方法的具体数学过程.
SRI方法的互相关和叠加过程可以用如下公式表示(安圣培等, 2015):
(1) |
(2) |
其中,*为复共轭;m为整数,代表台站的序号,M为台站的总个数;φ(m, i)和φ(m, j)对应于炮点i和j激发,m点接收的折射波记录;虚拟道ψ(i, j)m对应于第m个台站在i、j两点的折射波互相关,Φ(i, j)对应于叠加虚拟道;SRI方法生成的虚拟道叠加次数λ=M.
SRI-SNV方法将目标虚拟道和相邻虚拟道进行叠加,生成叠加虚拟道.其中,目标虚拟道由所选参考地震道和目标地震道在不同台站集中的折射波记录经互相关后产生,相邻虚拟道则由与所选道同侧等距相邻道的折射波经互相关后产生.然后将目标虚拟道与相邻虚拟道进行叠加.针对共台站集记录,假设目标虚拟道对应于不同台站集中第i、j两点的折射波互相关,则相邻虚拟道ψ(i-n, j-n)m和叠加虚拟道Φ(i, j)的表达式可以写为:
(3) |
(4) |
其中,m为台站的序号;M为台站的总个数;n为整数,代表单个台站所对应目标虚拟道和相邻虚拟道的参考地震道序号间隔;N为相邻虚拟道总数,本文选择的相邻虚拟道对称分布,故N为偶数,N/2为整数;φ(m, i-n),φ(m, j-n)分别代表生成虚拟道的不同偏移距折射波记录.当N=0时,无相邻虚拟道叠加,对应于SRI方法叠加生成虚拟道的过程,此时目标虚拟道ψ(i, j)m和叠加虚拟道Φ(i, j)的表达式与SRI方法一致,且虚拟道叠加次数λ=M.当N≠0时,表示有相邻虚拟道叠加,此时虚拟道叠加次数λ的取值与M和N的取值有关,且满足关系式:
(5) |
由于叠加次数增加,SRI-SNV方法对噪声的压制能力优于SRI方法,下文将通过数值实验和实际数据测试来展示这一点.
2 数值测试 2.1 模型与数据为了测试SRI-SNV方法对折射波的增强效果,本文设计了折射界面起伏的速度模型,如图 3所示.模型的上层和下层的纵波速度分别为2500 m·s-1和4000 m·s-1,纵横波速度比设定为常数1.8.模型中的速度界面分别在20.56~30.72 km和32.80~49.44 km之间存在两个高度不同的凸起.图 3还展示了观测系统,三个台站间距为0.8 km,631个炮点间距为0.08 km.
正演采用弹性波伪谱法(马锐和邹志辉, 2017),震源子波采用主频为5 Hz的雷克子波.正演模拟所生成的理论地震数据如图 4a所示.为模拟实际情况,理论数据添加了75%的随机噪声,以测试本文方法在低信噪比情况下的效果,如图 4b所示.其中,75%指的是噪声最大振幅与位于34.8 km处地震道折射波绝对值的最大振幅之比.
干涉法互相关步骤的输入数据要使用时窗截取,每道数据的时窗设置按如下步骤执行:首先,粗略估计初至时间,设定时窗位置;然后,利用估计的初至时间直接截取一定长度的折射波,截取时窗宽度的选择要足够宽,以保证折射波在截取的范围之内,本文选择的截取时窗宽度为1.4 s;最后,在截取的边缘处以高斯函数作为衰减函数平滑过渡到0值,高斯函数f(t)的表达式如下:
(6) |
其中,a为常数1;e为自然常数;b为高斯系数,本文取100;t为时间;t0为截取边缘处的时间值.由于边缘两端都需要施加高斯函数,故每道数据都需选取两个t0.实际应用中,我们将粗略初至作为前端边缘的t0,而将粗略初至与时窗宽度之和作为后端边缘的t0.
需要注意的是,本文干涉法的卷积步骤的输入均采用了原始地震数据,并未进行时窗截取,以保证不会遗漏重要的震相信息.
2.2 增强效果分析将图 4黑色虚矩形框中的折射波作为输入数据,利用SRI和SRI-SNV分别增强后的结果如图 5所示.图 5a显示,使用SRI方法对折射波进行增强后,所生成的折射波虽然得到一定的增强,但其附近仍然存在较强噪声.相比之下,图 5b中SRI-SNV方法增强后的折射波形的连续程度明显改善,并与图 4a中的无噪声波形接近.这种增强效果的差异主要来源于虚拟道叠加次数的不同.
图 6展示了SRI-SNV方法和SRI方法的虚拟道叠加效果.在SRI方法中虚拟道仅有3次叠加(图 6a),而在SRI-SNV方法中虚拟道叠加次数达到39次(图 6c),叠加次数是SRI方法的13倍.虚拟道叠加次数λ的取值需要选择恰当的邻近道范围,过小或者过大的λ值均不利于增强效果的改善.这是因为,当λ的取值偏小时,虚拟道叠加次数仍然较低,导致虚拟道叠加过程的抑噪效果有限;而当λ的取值偏大时,虽然虚拟道叠加次数较多,但由于用于叠加的虚拟道距离原始虚拟道较远、差异较大,难以满足相邻虚拟道叠加的假设条件,这些较大差异的虚拟道参与叠加后将导致波形的叠加相消不充分而引入更多噪声.对于本文情况,选择λ=39较为恰当.对比SRI方法(图 6b)和SRI-SNV方法(图 6d)的叠加结果可见,充足的虚拟道叠加次数使最终生成的虚拟道信噪比明显改善,其在峰值处与理论虚拟道形态相似,更加接近于脉冲函数.由于虚拟道代表了高信噪比参考地震道与低信噪比目标地震道之间的关系,因此准确地构建虚拟道将提高重构目标地震道的精度.
为显示折射波增强对初至拾取精度的影响,采用STA/LTA与多道互相关结合的方法拾取增强后的初至,并与无噪声理论数据的初至拾取结果作对比,如图 7a和7b所示.可以看出,SRI方法增强后的初至拾取结果与无噪声数据拾取结果的吻合度较差,有相当一部分的初至发生了偏离.而SRI-SNV方法增强后的初至拾取结果与无噪声数据匹配较好,初至基本上能与无噪声数据初至重合.进一步将两种方法增强后的初至拾取结果与对应无噪声数据初至拾取结果作差,拾取误差如图 8所示.一般而言,T/4(T为主周期)是判断初至拾取误差精度的标准,如果拾取误差小于T/4,则认为该初至能应用于层状模型中的走时层析成像(Mallinson et al., 2011).由图 8可知,SRI-SNV方法增强后的初至走时拾取误差有97.6%小于50 ms,基本都能满足小于T/4(主周期T=0.2 s)的要求.拾取误差的大小也可以通过其与理论初至之间的标准差衡量,SRI方法和SRI-SNV方法所对应拾取误差的标准差分别为0.151 s和0.028 s,说明SRI-SNV方法增强后的折射波有利于提高初至拾取的精度.另外,图 8还显示,相比于SRI方法,SRI-SNV方法对应的远偏移距处的拾取误差更小,表明SRI-SNV方法更适合于低信噪比远道的初至波增强.
为进一步分析方法的抗噪性,我们采用信噪比不同的理论数据对SRI方法和SRI-SNV方法进行了测试.将不含噪声数据施加振幅由低至高的随机噪声,分别生成五组数据.图 9显示了不同方法在不同程度信噪比情况下的折射波增强效果.从图中可以看出,当施加噪声水平较低(如25%)时,SRI方法和SRI-SNV方法对折射波的增强效果较好,说明高信噪比情况下两种方法都可行;当施加噪声水平达到一定程度(如75%、125%)时,SRI方法对折射波的增强效果较差,而SRI-SNV方法对折射波的增强效果较好,易于进行初至拾取;当施加噪声水平过高(如175%、225%)时,SRI方法的结果中完全无法识别折射波,而SRI-SNV方法的结果中仍然可以看到折射波的趋势,但两种方法增强后的折射波信噪比都过低而难以用于拾取初至.上述对比分析表明,相比于SRI方法,本文的SRI-SNV方法更适用于低信噪比情况下的折射波增强.
为检验SRI-SNV方法在实际资料中的增强效果,本文将其应用到南黄海北部某区域采集的广角OBS数据.野外采集的炮间距为125 m,台站间距为6 km,具体观测系统见Zou等(2016).我们选择了三个相邻的OBS台站,每个台站截取其接收的1401个炮的数据.将原始广角OBS数据做带通滤波以滤除一部分噪声,滤波频带范围为1~20 Hz,所得到数据如图 10所示.可见,滤波后折射波的信噪比仍较低,尤其在大于28 km的偏移距,波形已较难识别,如图 10a所示.我们使用宽度为1.2 s的时窗对滤波后三个台站的数据分别进行截取,并将截取后的折射波作为输入,以进行初至波增强.图 11显示了C1号台站的折射波增强结果,增强区域对应图 10a黑色虚矩形框中的数据.图 11a显示,SRI方法增强后的折射波连续性不足,而且部分偏移距还出现了明显的噪声.这些噪声的振幅与增强后的折射波的振幅相当,有可能像理论测试那样产生拾取误差.图 11b是使用SRI-SNV方法增强后的结果.与数值测试相同,该处同样选择39次虚拟道叠加.结果显示,SRI-SNV方法结果中的折射波连续性优于SRI方法,且信噪比也较高.
地震波形增强效果与互相关虚拟道的准确程度直接相关.虚拟道准确程度的一个衡量标准是其与脉冲函数的相似程度.图 12对比了SRI方法和本文SRI-SNV方法的虚拟道叠加效果.结果显示,SRI-SNV方法所得到的虚拟道在42.5 s附近的峰值更尖锐,更接近脉冲样式.这种脉冲形式的虚拟道在用于卷积生成目标地震道的过程中会减少旁瓣的产生,就像图 11b中显示的那样.
为对比SRI和SRI-SNV两种方法对初至拾取的影响,我们采用STA/LTA与多道互相关结合的自动拾取方法(林凡生和邹志辉,2016)对它们所增强波形进行初至拾取,所获得初至如图 13中红点所示.图 13a显示,SRI方法增强后的初至容易产生更多的离群值,这些离群值很明显偏离了初至的趋势,若用于层析将会在速度反演中带来噪声.与此相比,SRI-SNV方法增强后的初至拾取结果与数据的匹配程度更好,稳定性更高,更符合真实情况下的拾取结果,如图 13b所示.同时,我们使用层析反演得到的速度模型(Zou et al., 2016)和最短路径射线追踪(Moser, 1991)计算了初至到时,并将其分别投影到增强后的数据上,如图 13中绿点所示.可见,虽然根据层析速度模型计算的初至在大部分位置上与SRI-SNV方法增强后初至吻合,但仍有一些细节与该方法增强后初至符合得不够好,比如30~40 km之间和55 km附近等.若使用增强后波形拾取的初至进行层析反演,将有助于反演出速度结构的更多细节特征.
针对广角OBS数据远偏信噪比较低和台站分布稀疏的问题,本文提出了基于相邻虚拟道叠加的超虚折射干涉法(SRI-SNV方法).该方法通过叠加相邻虚拟道大幅增加了生成虚拟道过程中的叠加次数,提高了虚拟道计算的准确度,可以用于增强较大台站间距和较低信噪比的广角OBS折射波.理论和实际计算结果均表明,相比于SRI方法,SRI-SNV方法的虚拟道更加接近脉冲函数的形式,减少了增强后折射波的旁瓣,更有利于低信噪比情况下的初至增强.SRI-SNV方法比SRI方法更有利于提高折射初至拾取的稳定性和精度,使用SRI-SNV方法增强后的波形所拾取初至误差的标准差远低于SRI方法,且这种差异在远偏移距处更为明显,表明SRI-SNV方法更有利于低信噪比的远偏移距地震折射波的增强.
致谢 感谢两位匿名审稿人对本文提出的建设性意见和建议.
Al-Hagan O, Hanafy S M, Schuster G. 2014. Iterative supervirtual refraction interferometry. Geophysics, 79(3): Q21-Q30. DOI:10.1190/geo2013-0210.1 |
Alshuhail A, Aldawood A, Hanafy S. 2012. Application of super-virtual seismic refraction interferometry to enhance first arrivals:A case study from Saudi Arabia. The Leading Edge, 31(1): 34-39. DOI:10.1190/1.3679326 |
Al-Shuhail A A. 2015. Improving automatic first-arrival picking by supervirtual interferometry:examples from Saudi Arabia. Arab. J. Geosci., 8(10): 8731-8740. DOI:10.1007/s12517-015-1804-9 |
An S P, Liang X H, Peng G X, et al. 2015. Recognization and identification of first-arrival refractions in low SNR prestack seismic data. OGP (in Chinese), 50(3): 451-459, 468. |
An S P, Hu T Y, Liu Y M, et al. 2017a. Automatic first-arrival picking based on extended super-virtual interferometry with quality control procedure. Explor. Geophys., 48(2): 124-130. DOI:10.1071/EG14120 |
An S P, Hu T Y, Peng G X. 2017b. Three-dimensional cumulant-based coherent integration method to enhance first-break seismic signals. IEEE Transactions on Geoscience and Remote Sensing, 55(4): 2089-2096. DOI:10.1109/TGRS.2016.2636336 |
Bharadwaj P, Schuster G T, Mallinson I. 2011. Super-virtual refraction interferometry: Theory.//81st SEG Annual Meeting, Expanded Abstracts. SEG, 3809-3813.
|
Bharadwaj P, Schuster G, Mallinson I, et al. 2012. Theory of supervirtual refraction interferometry. Geophys. J. Int., 188(1): 263-273. DOI:10.1111/j.1365-246x.2011.05253.x |
Bharadwaj P, Wang X, Schuster G, et al. 2013. Increasing the number and signal-to-noise ratio of OBS traces with supervirtual refraction interferometry and free-surface multiples. Geophys. J. Int., 192(3): 1070-1084. DOI:10.1093/gji/ggs087 |
Bharadwaj P, Nissen-Meyer T, Schuster G, et al. 2014. Enhancing core-diffracted arrivals by supervirtual interferometry. Geophys. J. Int., 196(2): 1177-1188. DOI:10.1093/gji/ggt452 |
Dong S Q, Sheng J M, Schuster J T. 2006. Theory and practice of refraction interferometry.//76th SEG Annual Meeting, Expanded Abstracts. New Orleans, Louisiana: SEG, 3021-3025.
|
Funck T, Andersen M S, Neish J K, et al. 2008. A refraction seismic transect from the Faroe Islands to the Hatton-Rockall Basin. J. Geophys. Res., 113(B12): 178-196. DOI:10.1029/2008jb005675 |
Hanafy S M, Al-Hagan O, Al-Tawash F. 2011. Super-virtual refraction interferometry: Field data example over a colluvial wedge.//81st SEG Annual Meeting, Expanded Abstracts. San Antonio, Texas: SEG, 3814-3818.
|
Lin F S, Zou Z H. 2016. Automatic seismic first-arrival picking method based on the combination of STA/LTA and multi-channel cross correlation.//12st National Safety Geophysical Series-Geophysics and Information Perception (in Chinese). Beijing: Chinese Geophysical Society, 74-80.
|
Liu J, Wu S Y, Shi J. 2015. Application of pre-stack noise suppression technique in the South Yellow Sea. Marine Geology and Quaternary Geology (in Chinese), 35(6): 165-173. DOI:10.16562/j.cnki.0256-1492.2015.06.017 |
Liu L H, Hao T Y, Lü C C, et al. 2015. Crustal structure of Bohai Sea and adjacent area (North China) from two onshore-offshore wide-angle seismic survey lines. J. Asian Earth Sci., 98: 457-469. DOI:10.1016/j.jseaes.2014.11.034 |
Lu K, Altheyab A, Schuster G T. 2014. 3D super-virtual refraction interferometry.//84th SEG Annual Meeting, Expanded Abstracts. SEG, 4203-4207.
|
Lv X M, An S P, Hu T Y, et al. 2018. Similarity-weighted super-virtual interferometry to enhance first breaks. Acta Scientiarum Naturalium Universitatis Pekinensis (in Chinese), 54(1): 87-93. DOI:10.13209/j.0479-8023.2017.071 |
Mallinson I, Bharadwaj P, Schuster G, et al. 2011. Enhanced refractor imaging by supervirtual interferometry. The Leading Edge, 30(5): 546-550. DOI:10.1190/1.3589113 |
Ma R, Zou Z H. 2017. Composite absorption boundary conditions of pseudo spectral elastic wave simulation. 2017 Geophysical Technology Symposium of Chinese Petroleum Society (in Chinese). Tianjin, 1052-1055.
|
Mooney W D, Weaver C S. 1989. Regional crustal structure and tectonics of the Pacific Coastal States; California, Oregon, and Washington. Memoir of the Geological Society of America, 172(1): 129-161. DOI:10.1130/MEM172-p129 |
Moser T J. 1991. Shortest path calculation of seismic rays. Geophysics, 56(1): 59-67. DOI:10.1190/1.1442958 |
Operto S, Charvis P. 1996. Deep structure of the southern Kerguelen Plateau (southern Indian Ocean) from ocean bottom seismometer wide-angle seismic data. J. Geophys. Res., 101(B11): 25077-25103. DOI:10.1029/96JB01758 |
Qiao B P, Guo P, Wang P, et al. 2014. Effectively picking weak seismic signal near the surface based on reverse virtual refraction interferometry. Chinese J. Geophys. (in Chinese), 57(6): 1900-1909. DOI:10.6038/cjg20140621 |
Ruan A G, Niu X W, Qiu X L, et al. 2011. A wide angle Ocean Bottom Seismometer profile across Liyue Bank, the southern margin of South China Sea. Chinese J. Geophys. (in Chinese), 54(12): 3139-3149. DOI:10.3969/j.issn.0001-5733.2011.12.014 |
Wang S J, Wang F Y, Zhang J S, et al. 2014. The P-wave velocity structure of the lithosphere of the North China Craton-Results from the Wendeng-Alxa Left Banner deep seismic sounding profile. Sci. China:Earth Sci., 57(9): 2053-2063. DOI:10.1007/s11430-014-4903-7 |
Wei X D, Ruan A G, Zhao M H, et al. 2011. A wide-angle OBS profile across Dongsha Uplift and Chaoshan Depression in the mid-northern South China Sea. Chinese J. Geophys. (in Chinese), 54(12): 3325-3335. DOI:10.3969/j.issn.0001-5733.2011.12.030 |
Wu Z L, Li J B, Ruan A G, et al. 2011. Crustal structure of the northwestern sub-basin, South China Sea:results from a wide-angle seismic experiment. Sci. China:Earth Sci., 55(1): 159-172. |
Yan P, Zhou D, Liu Z S. 2001. A crustal structure profile across the northern continental margin of the South China Sea. Tectonophysics, 338(1): 1-21. DOI:10.1016/s0040-1951(01)00062-2 |
Zou Z H, Liu K, Zhao W N, et al. 2016. Upper crustal structure beneath the northern South Yellow Sea revealed by wide-angle seismic tomography and joint interpretation of geophysical data. Geol. J., 51(S1): 108-122. DOI:10.1002/gj.2847 |
安圣培, 梁向豪, 彭更新, 等. 2015. 低信噪比地震数据折射初至的判定与识别. 石油地球物理勘探, 50(3): 451-459, 468. |
林凡生, 邹志辉. 2016.基于STA/LTA与多道互相关结合的地震波初至自动拾取方法.//国家安全地球物理丛书(十二)——地球物理与信息感知.北京: 中国地球物理学会, 74-80.
|
刘俊, 吴淑玉, 施剑. 2015. 南黄海地震资料叠前去噪技术应用. 海洋地质与第四纪地质, 35(6): 165-173. DOI:10.16562/j.cnki.0256-1492.2015.06.017 |
吕雪梅, 安圣培, 胡天跃, 等. 2018. 相似度加权的超虚干涉法加强初至波信号. 北京大学学报(自然科学版), 54(1): 87-93. DOI:10.13209/j.0479-8023.2017.071 |
马锐, 邹志辉. 2017.伪谱法弹性波模拟的复合吸收边界条件.//中国石油学会2017年物探技术研讨会论文集.天津, 1052-1055.
|
乔宝平, 郭平, 王璞, 等. 2014. 基于逆虚折射干涉法有效提取近地表弱地震信号. 地球物理学报, 57(6): 1900-1909. DOI:10.6038/cjg20140621 |
阮爱国, 牛雄伟, 丘学林, 等. 2011. 穿越南沙礼乐滩的海底地震仪广角地震试验. 地球物理学报, 54(12): 3139-3149. DOI:10.3969/j.issn.0001-5733.2011.12.014 |
王帅军, 王夫运, 张建狮, 等. 2014. 华北克拉通岩石圈二维P波速度结构特征——文登-阿拉善左旗深地震测深剖面结果. 中国科学:地球科学, 44(12): 2697-2708. |
卫小冬, 阮爱国, 赵明辉, 等. 2011. 穿越东沙隆起和潮汕坳陷的OBS广角地震剖面. 地球物理学报, 54(12): 3325-3335. DOI:10.3969/j.issn.0001-5733.2011.12.030 |
吴振利, 李家彪, 阮爱国, 等. 2011. 南海西北次海盆地壳结构:海底广角地震实验结果. 中国科学:地球科学, 41(10): 1463-1476. |