Seismogram Analysis and Fitting of Three Earthquakes,West Papua,Recorded at Stations ENH and KMI,China

The S wave velocity structure of the earth below Eastern Southeast Asia has been investigated by analyzing the seismogram from surface wave to multiple depth waves in the time domain and three Cartesian components simultaneously. The wave passes across the front area of subduction zone between the Philippine plate and the Asian plate. The main data are waveform comparisons, instead of the arrival times. The synthetic seismogram is calculated using the GEMINI method. The synthetic seismogram constructed by PREMAN global earth model deviates greatly from the measured one. To solve this problem, corrections are needed for the β speed structure. Corrections cover the gradient change of βh, which turns from negative to positive in upper mantle layers as in the PREMAN, change of earth crust depth and change of zero order coefficients of β velocity function in all earth mantle layers. So, the fitting is obtained, as well as the arrival time or the waveform of Love and Rayleigh surface waves, the S wave and the repetitive depth waves ScS2 and ScS3. This result reveals that the Southeast Asia, being stretched due to tectonic release, has a mantle in some parts with negative anomaly of S wave velocity and vertical anisotropy in all earth mantle layers.     

利用数字地震记录研究唐山震区台下的P、S波速度结构

直接采集数字地震记录的原始 P波、S波记录波形 ,与利用 Haskell矩阵传递方法得到的理论综合地震图进行对比 ,通过不断调整介质模型 ,最后得到台站下方的体波速度结构模型 ,此方法中得到的 S波速度结构不同于通过其他波形转换而来的结果 ,而是同 P波一样从原始波形拟合而来 ,同时获得地震台下方 P、S速度介质模型 ,这为我们讨论地下介质的性质提供了更多的信息。数字地震波形资料取自于唐山市遥测地震台网陡河和滦县的短周期地震记录 ,共选用了震中距范围为 30～ 60°、震源深度为 1 0 0～ 60 0 km的 5个不同地点的深源远震记录资料。结果表明 ,台站下方的纵波和横波速度分层结构相当吻合 ,说明拟合结果相当可靠。两个台站下的地壳厚度均为 38km,速度结构在 1 3km的低速层处出现差异 ,滦县的低速层比陡河要更厚一些。两个台站的地壳速度结构都是高低相同 ,而且陡河的高低差异较滦县更为显著。唐山震区滦县部分的震源比陡河附近的浅 ,而 1 976年唐山 7.8级地震发生在陡河附近 ,可见地震的发生与壳内低速层的差异幅度有关。从泊松比高值来看 ,地幔顶部存在低速层 ,该区域是地幔上隆区 ,可能有部分熔融     

3D P-WAVE VELOCITY STRUCTURE AT THE NORTHEASTERN MARGIN OF ORDOS BLOCK

Using the 7 100 absolute first arrivals of P waves and 91 513 relative P arrival times of 726 events at the northeastern margin of the Ordos block since 2009, the 3D fine structure of P wave velocity within the depth of 15km in the crust was inverted by the double difference seismic tomography method. The results show that there exist obvious high-speed continuous bodies in the northwest of the study area, and their lateral areas increase gradually with depth, while the velocity of east and south is relatively low. The velocity inhomogeneity exists and differs at different depths. The lateral differences of velocity are related to seismicity and faults. The 5~15km depth profile shows that earthquakes tend to occur in the area with relatively high velocity or high speed transition zones, which to some extent reflects the fragility of regional crustal media and the strong differential movement of faults in vertical and horizontal directions where the crust body is easy to absorb and store strain energy and generate major earthquakes. A "Y"-shape low-velocity channel is present in the lower crust around Liangcheng, corresponding to the NW-trending Heilaoyao-Shahukou fault set, which may reveal the migration path of the Late Tertiary-Quaternary basalt eruption. The Helingeer M6.2 earthquake in 1976 was related to the formation of the locking section of the thermal welding in this area. The three-dimensional fine structure of P wave velocity presented in this paper provides intuitive seismological evidence for physical and chemical properties of crustal media and the deep tectonic environment of earthquake preparation.     

多震相走时和波形三维地震成象方法及在祁连山地壳结构研究中的应用

本文的理论方法是以几何射线理论为基础发展起来的、天然地震走时反演技术及天然地震层析成像技术。它存在着震源函数与介质参数的解耦问题。本研究采取了五种方法来改善反演结果。包括,利用,Pg,Pn等震相增大约束条件;用已有精度较高的人工地震测深结果作速度约柬：用波形反演来修改模型,把诸多物理量开发出来互为约束,以修改后的模型再作反演,使解的稳定性大大提高：采用最优化过程,选择遗传算法。可以进行震源定位,走时反演,波形反演：得到任意深度的速度分布及从地表到Moho面的速度剖面。用于在祁连山地区的结果表明,这些层析剖面对认识大地构造、重大深部事件动力学是很有益的。     

Three—dimensional P velocity structure in Beijing area

A detail three-dimensional P wave velocity structure of Beijing,Tianjin and Tangshan area(BTT area)was deter-mined by inverting local earthquake data.In total 16 048 Pwave first arrival times from 16048 shallow and mid-depth crustal earthquakes,which occurred in and around the BTT area from 1992to 1999were used.The first arrival times are recorded by Northern China Unived Telemetry Seismic Network and Yanqing-Huailai Digital Seismic Network.Hypocentral parameters of 1 132 earthquakes with magnitude ML=1.7-6.2 and the three-dimensional P wave velocity structure were obtained simultaneously.The inversion result reveals the com-plicated lateral heterogeneity of P wave velocity structure around BTT area.The tomographic images obtained are also found to explain other seismological observations well.     

基于小波分析的场地波速测量方法

在场地波速测量中,由于噪声等因素的影响很难准确识别P、S波的初至时刻,致使波速结果存在很大的误差。本文给出了一种基于小波变换的波速测量的新方法。该方法利用波动信号的小波变换与弹性波群速度的关系准确识别弹性波初至时刻。弹性波小波变换的峰值时刻代表着以群速度传播的弹性波的初至时刻,使P波、S波的初至时刻的确定具有明确的物理意义,波速的结果准确、可靠、稳定。此外,波动信号的小波多尺度分析还可以确定地层中传播的弹性波的频散特性。最后,该方法在场地波速测量的实测信号的应用表明该方法可准确确定P、S波速。     

Tomographic inversion of P410s delay times for simultaneous determination of P and S velocities of the upper mantle beneath the Baltic Shield

Teleseismic data recorded by stations in the Swedish National Seismic Network (SNSN) are used for a study of upper mantle structure beneath the Baltic Shield using the receiver function technique. The data show very clear conversions from the 410 and 660 km discontinuities. The signals associated with P to S conversions at these discontinuities arrive 1-2 s earlier than predicted by global models such as IASP91 or PREM. We interpret this as a manifestation of higher than average velocities in the mantle beneath the shield, consistent with lower than average global temperatures. For a 1400 km profile along the network, we observe variations of around 1 second in delay times of P410s and slightly less for P660s. Under the assumption that the mantle discontinuities are at a given constant depth, the delay times of the mantle converted phases are tomographically inverted to reveal P and S velocity structure below the stations. Synthetic tests show that this tomographic inversion has the potential to resolve P and S velocity variations at structural scales adequate for upper mantle studies. Results from application to real data appear to be consistent with independently produced mantle velocity structures deduced from normal tomographic arrival time data. For the P velocity model, a north-dipping body of (relatively) low velocity is found for the central part of the profile at 58-64°N. A sharp contrast from low to high velocities that may be associated with the Proterozoic-Archean boundary is found at 66°N.     

联合利用走时与波形反演技术研究地壳三维速度结构(Ⅱ)——应用

使用阻尼最小二乘法进行震源参数和地壳三维速度结构的走时联合反演．所用资料为Ｓ波和Ｐ波到时差,并用人工地震资料的二维解释结果作为三维速度模型的特定约束条件．为建立初始模型,又利用天然地震构成了准二维剖面．在走时反演基础上,利用遗传算法进行了几个地震事件的波形反演尝试,并对走时反演获得的地壳速度结构模型的局部进行了修正．以３４°～４２°Ｎ,９４°～１１２°Ｅ作为研究区域,在该区域中收集了１９８６年以来大量地震的Ｓ波和Ｐ波到时差资料,７条人工地震二维速度剖面资料和２个数字化地震台的几个地震的三分向记录资料．对这些资料进行了处理,最后得出了０～２５ｋｍ深度不同截面的速度分布,并对所得结果进行了分析．     

芦山震区地壳三维P波速度精细结构及地震重定位研究

联合芦山地震序列5285个地震的50711条P波初至绝对到时数据及7294691条高质量的相对到时数据,利用双差地震层析成像方法联合反演了芦山震源区高分辨率的三维P波速度精细结构及5115个地震震源参数.反演结果表明,芦山主震震中为30.28°N,103.98°E,震源深度为16.38km,主震南西段余震扩展长度约23km,余震前缘倾角较和缓,主震北东段余震扩展长度约12km,余震前缘呈铲形,倾角较陡.芦山震源区P波三维速度结构表现出明显的横向不均匀性,近地表处的P波速度异常与地形起伏及地质构造密切相关:宝兴杂岩对应明显的高速异常,此异常由地表延伸到地下15km深度附近,而中新生代岩石表现为低速异常;大兴附近区域亦显示出小范围的大幅度高速异常,宝兴高速异常与大兴高速异常在10km深度附近相连,进而增加了芦山震源区的高低速异常对比幅度.在芦山主震的南西、北东两段速度结构存在着较大差异,芦山主震在水平向位于宝兴及大兴高速异常所包围的低速异常的前缘.主震南西段余震主要发生在倾向北西的高低速异常转换带上并靠近低速一侧,其下盘为低速异常,上盘为高速异常.而芦山主震北东段的余震主要分布在宝兴高速体与大兴高速体之间,主发震层向北西倾斜,主发震层上方的宝兴高速异常下边界出现一条南东倾向的反冲地震带,两地震带呈"y"型分布.     

北京地区上地壳三维细结构层析成像

利用三维不分块地震层析成像技术对北京东北部地区的三维人工地震莫霍反射波走时数据进行反演,得到了该区上地壳的三维P波、S波慢度及波速比结构图像,三种结构不同深度的切片图像揭示出慢度和波速比分布有呈北东向和北西向展布的特征.断裂活动引起沉积不均匀造成现今的速度分布,地震大多分布在波速比较低的块体内,说明地震不仅与断裂活动有关,有可能还与物质结构和性质相关;北京地区上地壳速度结构、断裂活动和物质性质三者之间存在密切关系.     

山东地区地壳三维纵横波速度和泊松比结构及其对地震活动性的影响

基于山东及邻区丰富的P波和S波到时数据反演获得了研究区内高精度的三维纵横波速度结构和泊松比异常分布形态。结果表明:2020年济南长清M_S4.1地震震中位于P波、S波高低速异常和高低泊松比异常过渡带,可能是区域构造应力下长清断裂发生左旋走滑运动的结果;2003年青岛崂山M_L4.1地震、崂山震群、乳山震群和长岛震群等的发生可能都受到流体的强烈影响,流体沿已有的较大断裂或相对完整岩体内的裂隙侵入,诱发断裂活动或裂隙破裂,从而导致中强地震或震群活动的发生。      

汶川地震余震区波速比与强余震的关系

利用2008年5月12日汶川地震余震的P波和S波到时数据,对龙门山地区波速比进行了计算。结果表明,该地区的强余震全部发生在波速比相对较高的地区,这应该是地下介质的破碎原因造成的。     

山西临汾流动地震观测台址P／S波速度结构测试

用单孔检层法测定山西临汾流动地震观测台址两个钻孔的P／S波速度结构,为研究地震观测台址速度结构对地震波接收效果的影响提供基础资料。与常规P／S波测井不同的是,在P波速度测试中采用高频24道水听器串共炮点接收,有效降低单点多炮接收因重复性问题带来的零时误差,提高P波速度测定的精度。该方法可为类似工程场址的高精度P／S波速度测试提供新途径。     

云南地区地壳中上部横波速度结构研究

根据云南地区的基阶瑞利波相速度频散资料，用面波层析成像方法反演得到该区域中上地壳S波速度结构. 给出了研究区域内在4个深度上的S波速度水平分布图像和沿100.5°E、24°N、25°N、26°N及27°N的S波速度-深度剖面图. 结果表明：在小江断裂与红河断裂围成的川滇菱形块体内，26～30km深度处的速度明显低于周边地区，其南段从地表到15km深度均为明显的低速区域. 云南地区的强震（M＞6.0）震中位置与S波速度分布图像具有明显的相关性，主要分布于高速与低速的过渡区域.     

多震相走时和波形三维地震成像方法及在朝鲜半岛和祁连山地壳结构研究的应用

本项成果包括：提出天然地震走时反演层析成像技术,采用下列方法使得处理结果得以改善：1)利用Pg,Sg,Pm,Sm,Pn,Sn等震相增大约束条件;2)用已有精度较高的人工地震测深结果作速度约束;3)用波形反演来修改模型,把波源,介质吸收,散射等全部物理特征集中反映在记录中,把诸多物理量开发出来互为约束,以修改后的模型再作反演,使解的稳定性大大提高;4)采用最优化过程,选择遗传算法。可以进行震源定位,走时反演,波形反演;5)得到任意深度的速度分布及从地表到Moho面的速度剖面。主要应用结果：对于朝鲜半岛南部,划分为5(沿纬度)*６(沿经度)*8(沿深度)块,对于中部分288块.得到从地表到M面的8个水平切面;中部地区沿纬度13个二维剖面及其Moho面深度分布。上述方法也用于祁连山中东段地壳三维结构成像加上地震台网数字记录,反演。该区属塔里木-阿拉善地块走廊过渡带与北祁连褶皱带;从剖面可看出该地块上地壳低速层厚,下地壳有-低速层。北祁连褶皱带盆地与隆起构造之间的起伏差异,显示古浪断裂与金强河断裂之间的深部差异与界线。两地区的结果表明,这些剖面对认识大地构造、地质结构的稳定性,深部事件的性质是很有益的,对地球动力学研究也有重要意义。     

Vp和Vp/Vs联合反演:在日本东北地区深部结构成像(英文)

本文提出了一种新的反演方法:通过采用纵、横波走时数据对(从相同的震源产生的P和S波被同一台站记录)来联合反演纵波速度(Vp)和纵、横波速度比(Vp/Vs),然后单独反演横波速度Vs,在反演过程中同时对地震参数进行定位。该方法不需要假设P和S波的射线路径一致,它是沿着P和S波射线路径计算相对慢度扰动值。该方法直接把Vp/Vs作为一个模型参数,由此能获得比采用从独立反演获得的Vp和Vs计算出Vp/Vs的方法更精确的速度比值。该新方法被应用到反演日本东北地区的壳幔速度及波速比结构的研究中,获得了较好的效果。反演结果表明,在日本东北地区,太平洋俯冲板块为一高Vp,高Vs和低Vp/Vs异常区,而在活火山下方的浅部地幔楔以及背弧深部地区为低Vp,低VS和高Vp/VS异常。虽然这些特征在前人的研究中已经报道过,但与前人的研究结果相比,本次研究所获得的Vp/Vs的空间分布具有较小的分散性,同时,它的分布特征能较好的与地震波速度结构相吻合。     

Locating microseismic events using borehole data

Constraining microseismic hypocentres in and around hydrocarbon reservoirs and their overburdens is essential for the monitoring of deformation related to hydraulic fracturing, production and injection and the assessment of reservoir security for CO₂ and wastewater storage. Microseismic monitoring in hydrocarbon reservoirs can be achieved via a variety of surface and subsurface acquisition geometries. In this study we use data from a single, subsurface, vertical array of sensors. We test an existing technique that uses a 1D velocity model to constrain locations by minimizing differential S‐to‐P arrival times for individual sensors. We show that small errors in either arrival time picks or the velocity model can lead to large errors in depth, especially near velocity model discontinuities where events tend to cluster. To address this issue we develop two methods that use all available arrival times simultaneously in the inversion, thus maximizing the number of potential constraints from to N, where N is the number of phase picks. The first approach minimizes all available arrival time pairs whilst the second approach, the equal distance time (EDT) method defines the hypocentre as the point where the maximum number of arrival time surfaces intersect. We test and compare the new location procedures with locations using differential S‐to‐P times at each individual sensor on a microseismic data set recorded by a vertical array of sensors at the Ekofisk reservoir in the North Sea. Specifically, we test each procedure's sensitivity to perturbations in measured arrival times and the velocity model using Monte Carlo analysis. In general, location uncertainties increase with increasing raypath length. We show that errors in velocity model estimates are the most significant source of uncertainty in source location with these experiments. Our tests show that hypocentres determined by the new procedures are less sensitive to erroneous measurements and velocity model uncertainties thus reducing the potential for misinterpretation of the results.     

中国大陆岩石标本在高压下的弹性波速的初步研究

对10种岩石的20个试件和1个氯化银试件测得了高压下纵波和横波的速度,其中3个试件的纵波测至20,000公斤/厘米~2。发现压力高于4000公斤/厘米~2以后,速度基本上随压力呈直线变化,但当压力高于10,000公斤/厘米~2时直线的斜率逐级降低。在横波速度测定中,为提高脉冲初至到时的测定精度,在发射和接收端分别加入耦合电路。利用S波波列的多点测量外推方法有效地提高了测量精度。最后对结果作了细致的误差分析。     

首都圈地区地壳P波和S波三维速度结构及其与大地震的关系

选用华北地震遥测台网和首都圈数字化地震遥测台网1993～2005年记录的2866个地震事件中的33487条P波和31822条S波的到时资料,计算得到了水平分辨率25km到50km之间的首都圈地区（385°N～41°N, 114°E～120°E）地壳三维P波和S波速度结构,并进一步获得泊松比分布.研究表明,首都圈地区P波和S波速度分布表现出强烈的横向不均匀性,浅层速度分布同地表地质结构分布相一致.分析得出震区强震多发生在低速体与高速体之间、低泊松比地区,且震源下部存在低速、高泊松比异常体.经过与其他地区研究结果相比较,对地震触发与流体的关系进行了探讨.流体在地震孕育及触发的过程中可能起着重要作用.     

河北省强震区内外深部S波速度结构特征研究

利用短周期记录的S波记录反演了沙城老震区及其邻区共3个台站下方的剪切波速度结构，并进行了对比分析。结合以前得到的邢台、唐山震区及其附近5个台站的速度结构，发现强震区具有几个大致相同的构造特征，即：（1）地壳厚度低于周围台站；（2）S波低速层速度低于周围地区，高低速层差异大；（3）具有较低的上地幔S波速度；（4）强震往往发生在低速层上面的高速层内或高低速层交界面上；（5）易震层往往是两个低速层夹持的高速层等，并对这种结构易于发震的机理给出了分析。总之强震的发生是与低速层紧密相磁的，而剪切波对于反映低速体构造非常敏感，因此S波速度结构研究对于揭示强震的深部构造背景，深入研究震源机制都是很有意义的。而且，强震我所共同具有的这种特征为我们今后强震危险区的划分和强震地点的预测提供了较为科学的依据。