中国东南沿海地区近震走时表

本文根据浙江、福建两省地震台网所记录的浙、闽、台三省及邻近海域地震的观测资料,以单层均匀地壳模型,利用统计和速度正反演方法,得到P_n、P_a、S_n、S_a波的速度值和地壳厚度H,编制了中国东南沿海地区近震走时表。统计表明,本走时表比J-B表及中国地震走时表更适合于我国东南沿海地区地震观测的实际,可供该区地震工作者使用,也可作为有关人员的参考资料。     

利用兰州小孔径地震台阵资料叠加观测走时曲线

使用兰州小孔径地震台阵记录的近10年地震观测垂直分量波形数据,采用长、短时间平均数比值方法(STA/LTA)叠加出适用于青藏高原东北缘地区的观测走时曲线.结果表明,兰州小孔径地震台阵独特的地理位置,基本上能够记录到不同震中距(0—180°)和不同方位的地震事件;使用不同频率滤波处理之后的地震数据记录叠加出相应的观测走时曲线,从观测走时曲线中可以识别出不同体波震相(P,PKIKP,PKP,PP,PPP,PcP,ScP,S,SS等)的到时及其观测走时曲线特征.这对识别地震各种震相,认识和研究地球内部精细结构等具有非常重要的科学意义.     

新型二进制走时表文件的编制与查询

为将新型二进制结构走时表文件实际应用到测震分析软件,利用Matlab软件编写二进制结构综合走时表文件的通用程序,给出查询震相走时的算法,并成功配置到地震分析软件中用于地震定位.使用者只要仿照该通用程序的结构框架,根据需要适当修改程序,就能快速生成适合地震台站定位的走时表文件.     

地震分析软件震相走时表文件结构研究

通过对现有地震分析软件所配置的震相走时表文件的结构类型及优劣进行分析,将震相走时表分为4种组合结构并分析其实用性。在IASPEI91震相走时表文件结构的基础上进行修改、补充和完善,从而提出了一种安全性高、适用面广、灵活性强、能存储各种结构走时表的新型走时表文件结构。应用这种新型结构震相走时表,地震分析软件可以灵活配置不同的走时表,实现精细定位,并为数字化测震台网快速利用测震资料编制和应用高精度的区域震相走时表创造了条件,对于地球深部构造研究具有一定的促进意义。     

山西地区分区近震走时表的编制与检验

主要以山西太原、大同、临汾３个遥测台网的记录，结合山西省内及邻省部分台站的地震观测资料，用勘探、爆破与测震方法得出山西地区地壳模型结构。用观测走时与理论走时对比，划分出６个区域，计算出＂山西地区近震走时表＂及＂山西南部地区近震走时表＂。检验结果表明，发震时刻、震中位置及深度等各项参数较Ｊ－Ｂ走时表均有极大的提高。经试用确认其适用于山西地区。     

五大连池地区区域走时表研究

对黑龙江省区域数字地震台网建立以来五大连池火山区域所发生的天然地震,挑选出其中定位精度较高的85个地震,利用观测的地震波走时反演地下介质速度结构,在此基础上建立基于单层平均速度模型的五大连池火山区域走时表.试用结果表明,使用新走时表处理资料精度较J-B走时表有明显提高.     

数字仪与模拟仪地震处理的对比分析及近震走时表的修改

利用都兰地震台数字地震仪记录资料，在地震处理中使用修改后的走时表，可以实现从模拟到数字的顺利过渡，更好的使用数字化食品进行地震观测。通过传递函数仿真计算震级的对比以及走时表修改后与修改前震中距的对比分析，找出了两的差异和结合点。     

编制区域地震波走时表方法探讨

根据区域地震台网的地震观测资料，采用网格的方法划分等区域，将观测走时和理论走时对比，计算出适合本地区的区域走时表。     

华南地区近震走时表（速度模型）的精度检验

利用１６６次地震资料，通过３种地壳速度模型的对比，按照规定的评判原则，对《华南地区近震走时表》（速度模型）对福建及其邻区的地震定位精度进行了检验．结果表明，华南地区地壳模型（走时表）比Ｊ－Ｂ模型（走时表）更适合于本区地震的定位．     

利用J-B地震走时表和IASPEI走时表计算P波残差分析黑龙江省国家测震台的区域地震误差成因

利用黑龙江省国家测震台接收到的日本本岛及周边海域、台湾岛东北部及附近海域、中国大陆东北深震区以及中国东北与俄罗斯交界处部分地震,使用J-B地震走时表和IASPEI走时表分别计算P波走时残差,探讨黑龙江省测震台站接收研究区域地震事件误差成因。     

核幔界面衍射波Pdiff震相最大 记录震中距探讨

智利M_S8.8特大地震在甘肃数字测震台网记录的最大震中距达179.7°, 通过读取该地震初至震相的观测走时, 计算该震相的平均观测慢度为4.104 s/°, 与P_diff震相的理论慢度4.439 s/°基本一致. 研究结果表明: 由于甘肃数字测震台网和甘东南野外观测流动台阵独特的地理位置, 记录到了至今全球最大震中距(179.7°)的核幔界面衍射波P_diff震相; P_diff震相起始观测走时明显滞后于IASP91理论走时; 核幔界面复杂的形态对不同方位传播的P_diff震相观测走时产生了一定的影响.      

Numerical modeling of global seismic phases and its application in seismic phase identification

Earthquake data include informative seismic phases that require identification for imaging the Earth's structural interior. In order to identify the phases, we created a numerical method to calculate the traveltimes and raypaths by a shooting technique based upon the IASP91 Earth model, and it can calculate the traveltimes and raypaths for not only the seismic phases in the traditional traveltime tables such as IASP91, AK135, but also some phases such as pPcP, pPKIKP, and PPPPP. It is not necessary for this method to mesh the Earth model, and the results from the numerical modeling and its application show that the absolute differences between the calculated and theoretical traveltimes from the ISAP91 tables are less than 0.1 s. Thus, it is simple in manipulation and fast in computation, and can provide a reliable theoretical prediction for the identification of a seismic phase within the acquired earthquake data.     

2021年11月17日江苏大丰海域MS5.0地震震源深度对比分析

选取江苏测震台网震中距2°以内18个地震台站记录清晰的波形资料,使用单纯形、LocSAT、HypoSAT地震定位方法,分别采用华南速度模型、郯庐断裂带中南段速度模型以及IASP91全球速度模型,重新测定江苏大丰海域M_S5.0地震震源深度,对获取的深度值进行对比分析。结果表明:3种地震定位方法测定的震源深度结果接近,约为11 ± 2 km;HypoSAT定位方法计算结果较为稳定,IASP91、郯庐断裂带中南段模型均适用于该区域震源深度计算。      

Fresnel volumes corresponding to PKP waves in the IASP91 model

Fresnel volumes provide a suitable tool for estimating the resolution limits or accuracy of tomographic methods based only on the arrival times of these seismic waves, which penetrate the deep interior of the Earth. The standard definition of the Fresnel volume was modified to include effects of multipathing, i.e., multivalued travel times. Computations were performed in the IASP91 model for three representative epicentral distances and two prevailing signal periods. Preliminary results have been obtained. They place a rather strong restriction on the resolution limits of tomographic methods, depending on the prevailing period, due to the relatively large fatness of PKP Fresnel volumes.     

Travel time source-specific station corrections related to lithospheric structures in the Mediterranean region

We compare the locations obtained from arrival times collected by the International Seismological Centre from a network of regional and teleseismic stations for a cluster of Italian earthquakes with the locations of the same events obtained by the dense national seismic network operated by the Istituto Nazionale di Geofisica e Vulcanologia. We find mislocations on the order of 15 km for epicentral coordinates and on the order of 25 km for depths calculated from the regional and teleseismic network and using the standard IASP91 travel times. These mislocations are generally larger than the sizes of the respective error ellipse semi-axes. We then show that systematic shifts of hypocentral coordinates can be substantially reduced by applying source-specific station corrections. Moreover, we find that the size of error ellipses characterizing the teleseismic locations is significantly reduced by the application of such corrections. Our travel time corrections are compared and found fairly consistent with information available in the literature on tomographic studies on the crust and upper mantle in the European-Mediterranean region.     

Commentary on “Earthquake Depths in the Crimea–Black Sea Region” by V.E. Kulchitsky,B.G. Pustovitenko and V.A. Svidlova

It is well known that the results of determining earthquake parameters depend to a large extent on data processing algorithms and velocity models of the seismic wave propagation medium used in solving hypocenter problems. In 1992, V.Yu. Burmin developed a hypocentric algorithm that minimizes the functional of distances between the points corresponding to the theoretical and observed travel times of seismic waves from an earthquake source to recording stations. The determination of the coordinates of earthquake hypocenters in this case is much more stable than for the commonly used minimization of the functional of discrepancies in the seismic wave arrival times at a station. Using this algorithm and the refined velocity model of the medium, V.Yu. Burmin and L.A. Shumlyanskaya reinterpreted the earthquake parameters for the Crimea–Black Sea region. The most important result of this reinterpretation was the conclusion about the occurrence of deep earthquakes with a source depth of more than 60 km in the region. This result contradicts the conventional beliefs about the seismicity of the region and therefore aroused strong criticism from experts directly involved in compiling the existing catalogs of regional earthquakes. These comments and criticisms are presented by V.E. Kulchitsky with coauthors in a work published in this issue of the journal. In the present paper, we analyze the comments in detail and respond. In particular, we show that the previously used methods of seismic data processing made it highly unlikely by default that deep earthquakes would appear in the results. As an example, we refer to the use of travel-time curves for depths down to 35 km. It is clear that deep earthquakes could not have been found with this approach.     

华东地区台站偏差的研究

用华东地区十年地震资料和迭代方法，在华东走时表的基础上，计算出百余地震台站的作为距离的函数的走时偏差值，将这些结果当做台站校正值用于地震定位，定位精度得到明显提高。     

兰州台阵响应函数及不同方位小地震事件FK分析结果

计算了兰州台阵不同频率的台阵响应函数,对台阵附近不同方位的4个小震进行了FK分析。结果显示兰州台阵2 Hz以下的台阵响应函数在原点附近呈现出对称的单主峰图样,而高于2Hz时台阵响应函数出现较多旁瓣,表明该台阵对于2 Hz以下的信号具有较好的分辨能力。小震的FK分析结果表明,根据甘肃省地震局地震目录和IASP91速度模型计算得到P波方位角和慢度与FK分析的结果总体上比较一致,但也存在一定差别,这种差别可能由所使用的理论模型和当地模型有差别以及未做台阵校正所引起。     

2010年4月14日青海玉树7.1级地震序列中小地震辐射能量的估计

本文利用震后架设的4个流动台站和震源区周围3个固定台站记录到的青海玉树7.1级地震序列波形资料,资料时间截止到2010年5月19日.经过传播路径效应和场地响应校正后,估算了玉树地震序列中80个中小地震的辐射能量.比较了辐射能量估计中不同方法之间结果存在的差异性,首次考虑了由于仪器频带宽度的限制而造成的中小地震辐射能量的低估与补偿问题,得到以下主要结论:(1)分别采用Andrews方法、震源实际谱方法(区分能量补偿前和后二种情况)和根据理论模型并将积分上限外推到所计算地震拐角频率的10倍的方法(理论谱方法)计算出的中小地震辐射能量4种结果基本上位于同一个量级水平之上,表明结果具有一定的可靠性和稳定性.这其中,Andrews方法得到的辐射能量最低,理论谱方法得到的结果最高,而由实际谱方法得到的地震辐射能量介于两者之间.(2)由古登堡经验公式得到的地震辐射能量远高于上述4种方法计算出的结果,能够达到理论谱方法结果的4倍,而且随震级的增大,由它带来的辐射能量的高估趋势更加明显.(3)由于没有考虑地震能量的补偿,Andrews方法得到的地震辐射量明显偏低.用两种不同方法进行地震辐射能量补偿后得到的结果比较接近,但考虑到仪器频带上限值设定的复杂性,认为由理论谱方法计算中小地震的辐射能量更加合适.(4)经过能量补偿后地震辐射能量要高于未作能量补偿的结果,能量补偿带来的辐射能量增加的作用对于小地震更加明显.在本文研究范围之内,不同方法得到的折合能量均明显存在随地震矩的增大而增大的现象,反映出玉树余震序列中大震地震辐射效率高于小震.