甘青区域地震波走时表的地壳模型的选取

本文汇集了从六十年代末到八十年代初,地震勘探和测震方法给出的甘青地区地壳结构资料,根据地震观测和分析地震波的需求,简化出二十种地壳模型,由三千多个地震P波资料,进行观测走时与理论走时的对比,筛选出编算甘青区域地震波走时表的双层地壳模型:厚度H_1=22km,V_(P1)=6.10km/s,V_(S1)=3.55km/s Ha=29.5km,V_(P2)=6.47km/s,V_(S2)=3.81km/s H=51.5km,V_(Pn)=8.17km/s,V_(Sn)=4.62km/s     

Research on the Crustal Velocity Model of the Yunnan Region and Its Application

This paper selects the records of 7,412 earthquakes,each recorded by more than 10 stations in Yunnan between 2009 and 2014 to acquire the traveltime curves.Meanwhile,for improving precision,linear analysis,reduced traveltime curve and interval stability analysis are conducted focusing on the records of 83 earthquakes with M_L≥3.0 recorded each by≥80%of the stations,and by combining predecessors'research results,the initial crustal velocity model of the study area is obtained.By selecting 200 earthquakes with M≥3.0 occurring in Yunnan between 2010 and 2014,using the Hyposat batch location processing method to iterate the initial velocity model,and performing fitting to S waves layered velocity structure,we obtain the crustal velocity model for the Yunnan region,namely,the 2015 Yunnan model,with:v_(P1)=6.01km/s,v_(P2)=6.60km/s,v_(Pn)=7.89km/s,H_1=20km,H_2=21km,v_(S1)=3.52km/s,v_(S2)=3.86km/s,v_(Sn)=4.43km/s.Analysis on earthquake relocations based on the new model shows that most earthquakes occurring in Yunnan are at a depth of 10km-20km of the upper crust.The March 10,2011 M_S5.8Yingjiang and August 3,2014 M_S6.5 Ludian earthquakes are relocated,and the focal depths determined with the new model are respectively close to the precise positioning result and hypocentral distance to the strong motion stations at the epicenters,indicating that the new one-dimensional velocity model can better reflect the average velocity structure of the study area.     

新疆地区一维地壳速度模型研究

利用新疆区域地震台网观测到的2009年1月—2014年7月Pn、Sn、Pg和Sg震相数据,综合使用线性拟合、折合走时、PTD定深方法和HypoSAT定位方法确定该地区Pg、Pb和Pn的平均传播速度(v_(Pg),v_(Pb),v_(Pn))、康拉德界面和莫霍面的深度(H_(conr)和H_(moho))范围,以速度和深度步长分别为0.1km/s、1km精度迭代计算样本数据,通过对比分析计算结果与全国地震统一编目和3400模型下样本数据的定位结果后,确定RMS平均值最小的一维速度模型。在新模型中v_(Pg)、v_(Pb)和v_(Pn)分别为6.10km/s、6.70km/s和8.20km/s,H_(conr)和H_(moho)分别为26km、54km。通过检验对比,认为本文获得的新模型优于新疆地区现有的3400模型。     

内蒙古东部地区一维地壳速度结构分析

采用内蒙古测震台网2009—2016年记录的内蒙古东部地区131个地震资料,使用速度拟合、分区扫面、折合走时方法,反演得到该区域速度模型:v_1=6.10 km/s、v_(Pb)=6.72 km/s、v_n=8.05 km/s、H_1=23 km和H_2=16 km。东部模型检验结果显示,定位残差均值较华南模型和2015内蒙最优模型有明显的降低,且更加均匀稳定;东部模型与编目定位震中差较华南与编目、2015内蒙最优模型与编目和华南与编目震中差均值降低1 km左右。可见,东部模型更适合内蒙古东部地区。     

云南地区地壳速度结构和地震活动性研究

云南(21~29°N,97~106°E)地处印度板块与欧亚板块中国大陆碰撞带的东缘,是中国大陆最活跃的地震活动区域之一,其地震活动频度高、震级大、分布广,属于板缘、板内地震混合型地区.本研究采用中国地震科学台阵探测在南北地震带南段布置的367个流动地震台站记录到的波形数据,人工挑选出良好约束的,ML震级3.0的近震P波走时数据,采用VELEST算法,联合反演得到了研究区的一维P波速度结构和震源参数.结果表明,云南地区一维P波速度模型在0和30 km之间速度介于5.96~6.60 km/s之间;在30到40 km之间的下地壳P波速度值为6.60~6.80 km/s;上地幔顶部平均P波速度为7.80 km/s,低于全球大陆速度模型均值(8.04~8.10 km/s).基于新模型的地震活动性分析表明云南地区地震事件大多发生在中上地壳(20 km),约20 km处地震活动性最强且各个震级范围地震均有发生,可能指示该深度范围活跃的构造运动和复杂的应力状态.     

运用改进的地震定位交切法重定位2001—2010年云南地震

传统地震定位方法利用震源轨迹确定震源位置,但基于均匀或横向均匀介质模型必然导致定位误差。为此对传统方法进行改进,发展适用于三维复杂地壳速度模型的地震定位交切法。利用最小走时树射线追踪技术,以离散方式准确计算三维复杂地壳速度模型中的震源轨迹,将震源定位于震源轨迹交汇的密集点。将该方法应用于云南地区地震重定位,得到较高定位精度。     

青藏高原东北缘地壳及上地幔顶部速度结构研究

本文利用青藏高原东北缘71个固定台站与418个流动台站记录到的天然地震事件资料,采用双差层析成像方法对近震走时数据进行反演,获得了研究区高分辨率的三维P、S波速度结构和地震重定位结果.研究结果表明,本文给出的P、S波速度模型较已有的全球模型能更好的解释体波走时与面波相速度观测资料.松潘—甘孜和祁连构造带下方20～40 km深度范围表现为显著的P、S波低速异常,其中松潘—甘孜地块的壳内低速层可能与地壳部分熔融有关,而祁连构造带的壳内低速层则可能与地壳增厚有关.精定位后的岷漳6.7级地震和九寨沟7.0级地震震源深度都位于脆性的上地壳.两个地震的震源区地处不同块体的边界,均处在高、低速过渡带.震源区的壳内低速层可能处于部分熔融或易于蠕变的状态,脆性上地壳更容易积累应变能,从而导致地震的发生.     

小湾水库库区最小一维速度模型研究*

从小湾水库台网2005—2008年记录的众多地震中,挑选出最少被4个台站接收到的高质量地震数780个,一共5230条P波和4883条S波到时资料。利用Kissling方法得到了小湾水库库区最小一维P波和S波速度模型以及台站校正值。反演后的最小一维P波速度模型走时均方根残差从0.81s减少到0.12s,数据方差从1.64s²减少到0.04s²;地震震源深度比原来增加大约1公里,震源分布更加集中;不同台站的校正值差异表征了小湾库区速度结构存在横向不均匀性。最后利用得到的最小一维速度模型和台站校正值进行重定位,结果地震的走时均方根残差明显减少,表明得到的最小一维速度模型可信度较高。     

双差地震定位法误差定量分析

利用双层地壳模型,在下地壳中构造若干地震事件构成事件簇,分别在真实震源位置中添加不同程度的噪声后作为初始震源位置进行重定位,重定位中分别使用均匀分布和仅分布于事件簇一侧的台站,分析了初始震源位置和台站分布情况对双差重定位结果的影响.研究结果表明,在台站均匀分布情况下,当初始震中位置偏离真实震中小于22.3 km时,双差重定位后事件的震中位置平均偏离真实震中小于3.1km,事件间的相对震中距离为2.9km与真实值3.35 km接近,重定位后事件深度分布比真实事件集中;重定位后事件间的相对震中位置与添加的噪声水平关系不大,表明双差重定位法得到的事件间的相对震中位置更加稳定.当台站仅分布于震群一侧时,重定位结果与台站均匀分布时基本相同,表明台站分布对双差重定位结果影响不大.最后对双差重定位中最大震源对距离参数MAXSEP的分析表明,当事件簇中震源对平均距离比MAXSEP稍大时,重定位结果中部分事件的震源位置可能是不可靠的.     

小湾水库库区最小一维速度模型研究

从小湾水库台网2005—2008年记录的众多地震中,挑选出最少被4个台站接收到的高质量地震数780个,一共5 230条P波和4 883条S波到时资料。利用Kissling方法得到了小湾水库库区最小一维P波和S波速度模型以及台站校正值。反演后的最小一维P波速度模型走时均方根残差从0.81s减少到0.12s,数据方差从1.64s2减少到0.04s2;地震震源深度比原来增加大约1公里,震源分布更加集中;不同台站的校正值差异表征了小湾库区速度结构存在横向不均匀性。最后利用得到的最小一维速度模型和台站校正值进行重定位,结果地震的走时均方根残差明显减少,表明得到的最小一维速度模型可信度较高。     

The Crustal Velocity Structure of Western Inner Mongolia

The terrain of Inner Mongolia is long and narrow, and the geological structure is complicated. The South China crustal velocity model and Inner Mongolias optimal crustal velocity model 2015 cannot fully meet the earthquake location requirements of Inner Mongolia. Based on the seismological observations produced by Inner Mongolia Seismic Digital Network from 2009 to 2016,the initial model was obtained by using the linear fit of the seismic phases and the converted travel time curve. The Hyposat results of 225 earthquakes that occurred in western Inner Mongolia were scanned using this model,and the velocity model for western Inner Mongolia was determined as follows: V1= 6. 06 km/s;VPb= 6. 61 km/s; Vn= 8. 12 km/s; H1= 30 m; and the Moho depth H = 44 km. Comparison of the test results of the new model and the reference model shows that the residual error of the new model and the mean deviation of the epicenter location have obviously decreased.     

乳山震群震源区地壳速度结构初步研究

运用距离乳山震群最近的乳山台2011-2014年远震波形资料,计算接收函数,确定震群震源区及邻近区域的地壳厚度和波速比。结果显示：乳山台下方各个方位的接收函数差异比较大,地壳速度结构呈现横向非均匀性;震源区与邻近区域存在明显差别,邻近区域的中、下地壳存在明显的高速区,震源区中下地壳中存在明显的低速区;震群中M_L3.0以上地震基本发生在高低速交界处。据此推断,乳山震群可能是由于中下地壳小范围内的物质的不均匀性导致上地壳破裂。     

Relocation of the M 8.0 Wenchuan earthquake and its aftershock sequence

We relocated M8.0 Wenchuan earthquake and 2706 aftershocks with M⩾2.0 using double-difference algorithm and obtained relocations of 2553 events. To reduce the influence of lateral variation in crustal and upper mantle velocity structure, we used different velocity models for the east and west side of Longmenshan fault zone. In the relocation process, we added seismic data from portable seismic stations close to the shocks to constrain focal depths. The precisions in E-W, N-S, and U-D directions after relocation are 0.6, 0.7, and 2.5 km respectively. The relocation results show that the aftershock epi-centers of Wenchuan earthquake were distributed in NE-SW direction, with a total length of about 330 km. The aftershocks were concentrated on the west side of the central fault of Longmenshan fault zone, excluding those on the north of Qingchuan, which obviously deviated from the surface fault and passed through Pingwu-Qingchuan fault in the north. The dominant focal depths of the aftershocks are between 5 and 20 km, the average depth is 13.3 km, and the depth of the relocated main shock is 16.0 km. The depth profile reveals that focal depth distribution in some of the areas is characterized by high-angle westward dipping. The rupture mode of the main shock features reverse faulting in the south, with a large strike-slip component in the north. Supported by the Basic Research Project of Institute of Geophysics, China Earthquake Administration (Grant No. DQJB08Z03)     

Source Parameters of the Deadly M<Subscript>w</Subscript> 7.6 Kashmir Earthquake of 8 October, 2005

During the last six years, National Geophysical Research Institute, Hyderabad has established a semi-permanent seismological network of 5–8 broadband seismographs and 10–20 accelerographs in the Kachchh seismic zone, Gujarat with a prime objective to monitor the continued aftershock activity of the 2001 M_w 7.7 Bhuj mainshock. The reliable and accurate broadband data for the 8 October M_w 7.6 2005 Kashmir earthquake and its aftershocks from this network as well as Hyderabad Geoscope station enabled us to estimate the group velocity dispersion characteristics and one-dimensional regional shear velocity structure of the Peninsular India. Firstly, we measure Rayleigh-and Love-wave group velocity dispersion curves in the period range of 8 to 35 sec and invert these curves to estimate the crustal and upper mantle structure below the western part of Peninsular India. Our best model suggests a two-layered crust: The upper crust is 13.8 km thick with a shear velocity (Vs) of 3.2 km/s; the corresponding values for the lower crust are 24.9 km and 3.7 km/sec. The shear velocity for the upper mantle is found to be 4.65 km/sec. Based on this structure, we perform a moment tensor (MT) inversion of the bandpass (0.05–0.02 Hz) filtered seismograms of the Kashmir earthquake. The best fit is obtained for a source located at a depth of 30 km, with a seismic moment, Mo, of 1.6 × 10²⁷ dyne-cm, and a focal mechanism with strike 19.5°, dip 42°, and rake 167°. The long-period magnitude (M_A ~ M_w) of this earthquake is estimated to be 7.31. An analysis of well-developed sPn and sSn regional crustal phases from the bandpassed (0.02–0.25 Hz) seismograms of this earthquake at four stations in Kachchh suggests a focal depth of 30.8 km.     

Rupture process of the M L=4.1 earthquake in Huailai Basin on July 20, 1995

On July 20, 1995, an earthquake of M _L=4.1 occurred in Huailai basin, northwest of Beijing, with epicenter coordinates 40.326°N, 115.448°E and focal depth 5.5 km. Following the main shock, seismicity sharply increased in the basin. This earthquake sequence was recorded by Sino-European Cooperative Huailai Digital Seismograph Network (HDSN) and the hypocentres were precisely located. About 2 hours after the occurrence of the main shock, a smaller event of M _L=2.0 took place at 40.323°N, 115.447°E with a focal depth of 5.0 km, which is very close to the main shock. Using the M _L=2.0 earthquake as an empirical Green’s function, a regularization method was applied to retrieve the far-field source-time function (STF) of the main shock. Considering the records of HDSN are the type of velocity, to depress high frequency noise, we removed instrument response from the records of the two events, then integrated them to get displacement seismogram before applying the regularization method. From the 5 field stations, P phases in vertical direction which mostly are about 0.5 s in length were used. The STFs obtained from each seismic phases are in good agreement, showing that the M _L=4.1 earthquake consisted of two events. STFs from each station demonstrate an obvious “seismic Doppler effect”. Assuming the nodal plane striking 37° and dipping 40°, determined by using P wave first motion data and aftershock distribution, is the fault plane, through a trial and error method, the following results were drawn: Both of the events lasted about 0.1 s, the rupture length of the first one is 0.5 km, longer than the second one which is 0.3 km, and the rupture velocity of the first event is 5.0 km/s, larger than that of the second one which is about 3.0 km/s; the second event took place 0.06 s later than the first one; on the fault plane, the first event ruptured in the direction γ=140° measured clockwise from the strike of the fault, while the second event ruptured at γ=80°, the initial point of the second one locates at γ=−100° and 0.52 km from the beginning point of the first one. Using far-field ground displacement spectrum measurement method, the following source parameters about the M _L=4.1 earthquake were also reached: the scalar earthquake moment is 3.3×10¹³ N·m, stress drop 4.6 MPa, rupture radius 0.16 km. Contribution No. 99FE2022, Institute of Geophysics, China Seismological Bureau. This study is supported by the Chinese Joint Seismological Science Foundation (95-07-411).     

Source Parameters of the 8 October, 2005 M<Subscript>w</Subscript>7.6 Kashmir Earthquake

During the last six years, the National Geophysical Research Institute, Hyderabad has established a semi-permanent seismological network of 5 broadband seismographs and 10 accelerographs in the Kachchh seismic zone, Gujarat, with the prime objective to monitor the continued aftershock activity of the 2001 M_w7.7 Bhuj mainshock. The reliable and accurate broadband data for the M_w 7.6 (8 Oct., 2005) Kashmir earthquake and its aftershocks from this network, as well as from the Hyderabad Geoscope station, enabled us to estimate the group velocity dispersion characteristics and the one-dimensional regional shear-velocity structure of peninsular India. Firstly, we measure Rayleigh- and Love-wave group velocity dispersion curves in the range of 8 to 35 sec and invert these curves to estimate the crustal and upper mantle structure below the western part of peninsular India. Our best model suggests a two-layered crust: The upper crust is 13.8-km thick with a shear velocity (Vs) of 3.2 km/s; the corresponding values for the lower crust are 24.9 km and 3.7 km/sec. The shear velocity for the upper mantle is found to be 4.65 km/sec. Based on this structure, we perform a moment tensor (MT) inversion of the bandpass (0.05–0.02 Hz) filtered seismograms of the Kashmir earthquake. The best fit is obtained for a source located at a depth of 30 km, with a seismic moment, Mo, of 1.6 × 10²⁷ dyne-cm, and a focal mechanism with strike 19.5°, dip 42°, and rake 167°. The long-period magnitude (M_A ~ M_w) of this earthquake is estimated to be 7.31. An analysis of well-developed sPn and sSn regional crustal phases from the bandpassed (0.02–0.25 Hz) seismograms of this earthquake at four stations in Kachchh suggests a focal depth of 30.8 km.     

Earthquake relocation and 3-dimensional crustal structure of P-wave velocity in cen-tral-western China

IntroductionUnderstandingthemechanismofcontinentalearthquakesisveryimportantforseismichaz-ardpreventionandearthquakeprediction.Themodernseismotectonictheoryandtheideaofearthquakepredictionaredevelopedmainlyfromthestudiesoninterplateearthquakes,whicharedifficulttoexplainthephenomenaofintraplateearthquakes,suchasthecontinentalearthquakesoccurredinChinesemainland.Whiletheinterplateearthquakesoccurredalongtheplatebounda-ries,theintraplateearthquakesdistributediffuselyintheinterioroftheplates.Thus…     

Rupture process of November 6, 1988, Lancang-Gengma, Yunnan, China, earthquake of M s=7.6 using empirical Green’s function deconvolution method

The empirical Green’s function deconvolution technique is applied to the November 6, 1998,M _s=7.6, Yunnan, China earthquake to study the rupture process of this event. An aftershock ofM _s=6.3 is taken as the empirical Green’s function. Recordings of these two events at six stations of China Digital Seismographic Network (CDSN) are employed. Deconvolution results show that this event is a relatively simple event. Directivity analysis indicates that the rupture was initiated at hypocenter and propagated bilaterally and early symmetrically towards the northwest and southeast directions with a total length of 70 km and a time duration of 19 s. The rupture velocity is estimated to be about 2.0 km/s. Contribution No. 98A01004, Institute of Geophysics, State Seismological Bureau, China. This study is supported by the Chinese Joint Seismological Science Foundation(95-07-411).     

Near surface velocity and QS structure of the Quaternary sediment in Bohai basin,China

Heavily populated by Beijing and Tianjin cities, Bohai basin is a seismically active Cenozoic basin suffering from huge lost by devastating earthquakes, such as Tangshan earthquake. The attenuation (Q_P and Q_S) of the surficial Quaternary sediment has not been studied at natural seismic frequency (1?10 Hz), which is crucial to earthquake hazards study. Borehole seismic records of micro earthquake provide us a good way to study the velocity and attenuation of the surficial structure (0?500 m). We found that there are two pulses well separated with simple waveforms on borehole seismic records from the 2006 M_W4.9 Wen'an earthquake sequence. Then we performed waveform modeling with generalized ray theory (GRT) to confirm that the two pulses are direct wave and surface reflected wave, and found that the average v_P and v_S of the top 300 m in this region are about 1.8 km/s and 0.42 km/s, leading to high v_P/v_S ratio of 4.3. We also modeled surface reflected wave with propagating matrix method to constrain QS and the near surface velocity structure. Our modeling indicates that Q_S is at least 30, or probably up to 100, much larger than the typically assumed extremely low Q (~10), but consistent with Q_S modeling in Mississippi embayment. Also, the velocity gradient just beneath the free surface (0?50 m) is very large and velocity increases gradually at larger depth. Our modeling demonstrates the importance of borehole seismic records in resolving shallow velocity and attenuation structure, and hence may help in earthquake hazard simulation.