用平面波延拓方程进行地震数据的叠前速度反演

本文讨论地震勘探数据的叠前速度反演方法及其在海洋地震勘探数据上的反演试验.反演主要的计算步骤是:1.采用Fourier-Hankel变换把球面波分解为平面谐波;2.用平面谐波的延拓方程将上行波与下行波同时向下延拓,并计算每一层底部的反射系数和下一层的波阻抗;3.用最小二乘法从波阻抗中确定该层的声波速度.重复第2步与第3步,直到某一预定深度时结束.通过反演试验,对地震振幅比例的改变,子波变形,以及第1层速度和密度的误差对反演方法的稳定性及其精度的影响进行了分析.还通过实际海洋地震勘探数据的反演试验,对这一方法在地震勘探中的应用前景作了论述.     

基于人工地震资料的珠江口地区最小一维P波速度模型

珠江口地区位于南海北部大陆的边缘,具有洋陆过渡型地壳特征,且NE向滨海断裂带从其中穿过,强震风险不可忽视.文中基于2015年珠江口海陆联合三维人工地震探测数据,人工进行初至P波震相拾取,并使用VELEST程序分别反演了陆域和海域的最小一维P波速度模型(走时残差均方根最小)和台站校正结果.台站校正结果的空间分布与区域地形...     

利用爆破事件分析不同速度模型下常规定位方法定位结果

地震定位对速度模型的依赖性很强。四川地区地形复杂,常规工作中可选取多种速度模型进行定位。川西龙门山断裂带为东南部四川盆地和西北部青藏高原东部山区的明显分界线,近年在此断裂带上发生多次较大地震。对发生在该断裂带附近的6个爆破事件和15个天然地震重新定位,并对比结果。研究表明,相同台站包围情况下,川滇3D速度模型稳定性最好,但对浅表爆破不太准确。相比HypoSat（一维速度模型）组合,台站分布对Hypo2000（一维速度模型）和Hypo2000（赵珠速度模型）组合的定位结果影响较大。     

利用爆破事件分析不同速度模型下常规定位方法定位结果

地震定位对速度模型的依赖性很强.四川地区地形复杂,常规工作中可选取多种速度模型进行定位.川西龙门山断裂带为东南部四川盆地和西北部青藏高原东部山区的明显分界线,近年在此断裂带上发生多次较大地震.对发生在该断裂带附近的6个爆破事件和15个天然地震重新定位,并对比结果.研究表明,相同台站包围情况下,川滇3D速度模型稳定性最好,但对浅表爆破不太准确.相比HypoSat(一维速度模型)组合,台站分布对Hypo2000(一维速度模型)和Hypo2000(赵珠速度模型)组合的定位结果影响较大.     

利用DE算法反演地壳速度模型和地震定位

利用差异演化（Differential Evolution）非线性全局优化算法,设计了一种反演地壳速度模型和进行地震定位的方法,并给出了反演结果的具体分析.利用有限差分算法计算速度模型的走时场,可以节省大量的计算量,加快计算过程.反演得到的地壳速度模型和地震的震源参数可以直接用于地震层析成像研究,还可以利用地震层析成像得到的三维速度结构对地震重新定位,从而得到较为精确的震源参数.地壳速度模型的反演方法也可以用于三维速度结构的反演.     

Application of 3-D Crustal and Upper Mantle Velocity Model of North America for Location of Regional Seismic Events

—?Seismic event locations based on regional 1-D velocity-depth sections can have bias errors caused by travel-time variations within different tectonic provinces and due to ray-paths crossing boundaries between tectonic provinces with different crustal and upper mantle velocity structures. Seismic event locations based on 3-D velocity models have the potential to overcome these limitations. This paper summarizes preliminary results for calibration of IMS for North America using 3-D velocity model. A 3-D modeling software was used to compute Source-Station Specific Corrections (SSSCs(3-D)) for Pn travel times utilizing 3-D crustal and upper mantle velocity model for the region. This research was performed within the framework of the United States/Russian Federation Joint Program of Seismic Calibration of the International Monitoring System (IMS) in Northern Eurasia and North America.¶An initial 3-D velocity model for North America was derived by combining and interpolating 1-D velocity-depth sections for different tectonic units. In areas where no information on 1-D velocity-depth sections was available, tectonic regionalization was used to extrapolate or interpolate. A Moho depth map was integrated. This approach combines the information obtained from refraction profiles with information derived from local and regional network data. The initial 3-D velocity model was tested against maps of Pn travel-time residuals for eight calibration explosions; corrections to the 3-D model were made to fit the observed residuals. Our goal was to find a 3-D crustal and upper mantle velocity model capable predicting Pn travel times with an accuracy of 1.0–1.5 seconds (r.m.s.).¶The 3-D velocity model for North America that gave the best fit to the observed travel times, was used to produce maps of SSSCs(3-D) for seismic stations. The computed SSSCs(3-D) vary approximately from +5 seconds to ?5 seconds for the western USA and the Pre-Cambrian platform, respectively. These SSSCs(3-D) along with estimated modeling and measurement errors were used to relocate, using regional data, an independent set of large chemical explosions (with known locations and origin times) detonated within various tectonic provinces of North America. Utilization of the 3-D velocity model through application of the computed SSSCs(3-D) resulted in a substantial improvement in seismic event location accuracy and in a significant decrease of error ellipse area for all events analyzed in comparison both with locations based on the IASPEI91 travel times and locations based on 1-D regional velocity models.     

Crustal structure in Dabieshan UHP metamorphic belt and its tectonic implication

IntroductionGeologistsfirstlyfoundcoesite-bearingecologitesattheendof1980'sandthenthemicrodiamond(Xu,elal,1992)inDabieshanarea.Theultra-highpressure(UHP)metamorphismandthegeodynamicprocessesofDabieorogenhaveattractedmanygeoscientists(Wang,etal,1989,O...     

Regional minimum 1-D P-wave velocity model for a new seismicity catalogue with precise and consistent earthquake locations in southern Iran

For faster and more robust ray tracing in 1-D velocity models and also due to the lack of reliable 3-D models, most seismological centers use 1-D models for routine earthquake locations. In this study, as solution to the coupled hypocenter-velocity problem, we compute a regional P-wave velocity model for southern Iran that can be used for routine earthquake location and also a reference initial model for 3-D seismic tomography. The inversion process was based on travel time data from local earthquakes paired reports obtained by merging the catalogues of Iranian Seismic Center (IRSC, 6422 events) and by the Broadband Iranian National Seismic Network (BIN, 4333 events) for southern Iran in the period 2006 through July 2017. After cleaning the data set from large individual reading errors and by identifying event reports from both networks belonging to same earthquake (a process called event pairing), we obtained a data set of 1115 well-locatable events with a total number of 24,606 P-wave observations. This data set was used to calculate a regional minimum 1-D model for southern Iran as result of an extensive model search by trial-and-error process including several dozens of inversions. Significantly different from previous models, we find a smoothly increasing P-velocity by depth with velocities of 5.8 km/s at shallow and velocities of 6.4 km/s at deepest crustal levels. For well-locatable events, location uncertainties are estimated in the order of ±?3 km for epicenter and double this uncertainty for hypocentral depth. The use of the minimum 1-D model with appropriate station delays in routine hypocenter location processing will yield a high-quality seismic catalogue with consistent uncertainty estimates across the region and it will also allow detection of outlier observations. Based on the two catalogues by IRSC and BIN and using the minimum 1-D model and station delays for all stations in the region, we established a new combined earthquake catalogue for southern Iran. While the general distribution of the seismicity corresponds well with that of the two individual catalogues by IRSC and BIN, the new catalogue significantly enhances the correlation of seismicity with the regional fault systems within and between the major crustal blocks that as an assembly build this continental region. Furthermore, the unified seismic catalogue and the minimum 1-D model resulting from this study provide important ingredients for seismic hazard studies.     

2-D crustal structure of Changbaishan-Tianchi volcanic region determined by seismic traveltime inversion

IntroductionTrial-and-error forward modeling of wide-angle seismic reflection/refraction traveltimes for 2-D velocity structure is extremely time-consuming, even for experienced data interpreters. For wide-angle seismic reflection/refraction experiments that consist of numerous shots along a single line, it is quite difficult through repeated trial-and-error forward modeling to construct a 2-D model that fits the data within acceptable limits (Cerveny, et al, 1977; ZHANG, et al, 200 . In ad…     

A multi-scale 3-D crust velocity model in the Hefei-Chao Lake area around the southern segment of Tanlu Fault Zone

Regional high-precision velocity models of the crust are an important foundation for examining seismic activity, seismogenic environments, and disaster distribution characteristics. The Hefei-Chao Lake area contains the main geological units of Hefei Basin, with thick sediments and the Chao Lake depression. Several major concealed faults of the southern NNE-trending Tanlu Fault Zone cross this area. To further explore the underground distribution characteristics of the faults and their tectonic evolutionary relationship with adjacent tectonic units, this study used ambient noise data recorded by a seismic array deployed in Hefei City and Chao Lake, constructing a 3-D velocity model at the depth of 1–8 km. Then a multi-scale high-resolution 3-D velocity model of this area was constructed by this new upper crustal velocity model with the previous middle and lower crustal model. The new model reveals that a high-velocity belt is highly consistent with the strike of the Tanlu Fault Zone, and a low-velocity sedimentary characteristic is consistent with the Hefei Basin and Chao Lake depression. The distribution morphology of high and low velocity bodies shows that the sedimentary pattern of Hefei-Chao Lake area is closely related to the tectonic evolution of the Tanlu Fault Zone since the Mesozoic. This study also identifies multiple low-velocity anomalies in the southeastern Hefei City. We speculate that strong ground motion during the 2009 Feidong earthquake (magnitude of 3.5) was related to amplification by the thick sediments in the Hefei Basin. We also discuss further applications of multi-scale high-resolution models of the shallow layer to strong ground motion simulations in cities and for earthquake disaster assessments.     

地壳三维构造反演和速度层析成像

本文给出利用空间深地震测深反射资料重建地壳三维构造和速度分布的方法，适用于任意分层并且每层可由若干断块组成的地壳模型。通过采用适当的反演方法获得地壳界面的三维构造，并在此结果基础之上，进一步将剩余走时残差归因于地壳速度的不均匀性，采用东分块的反演方法，重建地壳三维速度分布，数值模拟的结果表明算法是有效的。     

青藏高原东部及其边缘地区的地壳上地幔三维速度结构

利用在青藏高原东部及其邻近地区记录到的１万余条近震到时资料,反演该地区的地壳上地幔三维速度结构。采用网格点模型描述三维速度结构,模型维数为２２２２６,网格点间距水平向为１００ｋｍ,垂直向为２０ｋｍ,网格点之间的速度值通过线性插值给出。采用改进了的快速三维射线追踪方法,确定三维非均匀介质中的地震射线路径和理论走时。反演结果显示,青藏高原南部的上地壳中（３０ｋｍ左右的深度）存在一低速区,这和面波反演的结果一致,羌塘块体下地壳有明显的低速异常带,青藏公路沿线的垂直速度剖面显示出岩石层受挤压增厚的构造特征。     

上海及其邻近地区三维地壳结构层析成像

利用上海及其邻近地区地震台网的地震记录及人工爆破资料,采用地震层析成像方法反演研究区地壳三维P波、S波速度结构模型。成像结果表明:不同深度的P波、S波速度扰动呈北西向展布特征;地震大多数分布在低速区的块体内,说明地震不仅与断裂活动有关,有可能还与物质结构和性质有关;上海及其邻近地区上地壳速度结构,断裂活动和物质性质三者之间存在密切关系。     

用地震走时反演长白山天池火山地区的二维地壳结构

应用一种同时获得地壳上地幔二维速度结构和界面形态的地震走时反演算法，对在长白山天池火山区L1剖面获得的人工地震折射和宽角反射走时数据进行反演，得到了该剖面的地壳上地幔二维速度结构和构造．结果表明:长白山天池火山口下方中下地壳存在一倒三角状的低速体，壳内反射界面和莫霍面一致呈下凹形态，幅度2～6km不等，形成一个地壳根，推测是长白山天池火山系统．最后将所得结果与SEIS83正演试错模拟结果进行了比较．      

南海西沙地块岛屿地震观测和海陆联测初步结果

为了研究南海西沙地块下方的地壳结构,我们在岛屿区架设了流动地震台站进行天然地震观测和海底地震仪的人工地震探测.本文利用远震接收函数方法和射线追踪方法对琛航台的远震数据和海陆联测数据进行了初步处理和分析,建立了琛航岛下方的一维横波速度模型以及横穿琛航岛的二维地壳结构模型.琛航岛顶部存在2 km厚的新生代低速沉积层,下地壳...     

地球横向不均匀结构对地表以及空间固定点同震重力变化的影响

本文提出一个新算法,用来高精度计算三维不均匀地球模型中地震位错引起的地表以及空间固定点同震重力变化.具体地说,我们首先把实际三维不均匀地球分解成球对称地球模型和对应的横向不均匀增量,分别进行计算,二者对应的计算结果分别称为球对称解和三维响应.由于球对称解可直接利用球对称地球模型位错理论计算得到,本文的目标是计算三维响应即地球的横向不均匀结构对同震重力变化的影响.然后,我们把三维响应再分为震源的响应和地球横向不均匀结构的响应,它们可分别借助对震源函数的扰动以及对平衡方程式的变分求解.本文推导出六个特殊点源位错引起的地表以及空间固定点同震重力变化计算公式(一个垂直走滑位错,两个相互垂直的倾滑位错,三个开裂位错),对这些公式进行适当组合就可以计算任意位置任意类型位错产生的同震重力变化,对应的计算公式同步给出.接着,依据36阶P波速度模型,我们利用岩石试验经验关系式推导出三维S波速度模型,密度模型,位场模型以及重力模型.最后,利用上述三维模型,本文计算出三种典型类型的点源位错产生的同震重力变化,结果显示三维响应与位错类型,震源深度都有关系,其最大响应占球对称解的0.5%左右,且在所有影响因素中S波速度模型影响最大.数值结果同时表明,三维响应中震源的响应与地球横向不均匀构造的响应处于同一量级.本文给出的地表和空间固定点同震重力变化计算公式可分别高精度解析地表重力和卫星重力观测数据(GRACE、GOCE等),提高大地测量数据理论解析水平.     

用可视化技术分析虚拟现实环境下三维地震数据场的研究

目的 用可视化方法分析三维地震数据，直接观察数据中隐藏的地质现象。方法：将第每个反射子波用一个粒子模型表示，来自同一反射面的子波粒子相互叠合成一个同相面。为加快显示速度，将一个三维地震数据体转换为三维二值图像，用有序扩展算法追踪同相面所在的三维连通域，细化后得到反射面的坐标数据。再在虚拟现实系统下的对反射面进行立体显示，直观地分析地下地层的结构。结果 实际数据试验得到了反射面的三维形态及断层的分布规律，立体显示更深刻地揭示了复杂地质现象及其规律。结论 只有在虚拟现实系统中进行三维地震数据的可视化，才能直观地看到地下地层的真实场景，真正实现三维数据的全三维解释，所提出的方法是有效。     

BARENTS16: a 1-D velocity model for the western Barents Sea

A minimum 1-D seismic velocity model for routine seismic event location purposes was determined for the area of the western Barents Sea, using a modified version of the VELEST code. The resulting model, BARENTS16, and corresponding station corrections were produced using data from stations at regional distances, the vast majority located in the periphery of the recorded seismic activity, due to the unfavorable land–sea distribution. Recorded seismicity is approached through the listings of a joint bulletin, resulting from the merging of several international and regional bulletins for the region, as well as additional parametric data from temporary deployments. We discuss the challenges posed by this extreme network-seismicity geometry in terms of velocity estimation resolution and result stability. Although the conditions do not facilitate the estimation of meaningful station corrections at the farthermost stations, and even well-resolved corrections do not have a convincing contribution, we show that the process can still converge to a stable velocity average for the crust and upper mantle, in good agreement with a priori information about the regional structure and geology, which reduces adequately errors in event location estimates.     

单源单缆方式采集的天然气水合物三维地震数据处理技术

作为潜在的能源储备,天然气水合物的重要性持续上升.针对天然气水合物的地震勘探也逐步由二维方式向三维方式转变.本文以单源单缆方式采集的天然气水合物三维地震数据为基础,探讨了地震数据成像的关键技术,在获得准确的三维地震成像数据体的同时,还得到三维地层速度数据体.三维地震数据更有利于描述含天然气水合物地层的地震反射特征,并刻画可能含天然气水合物地层的分布,从而拓宽对天然气水合物储层的认识.     

Modeling of Wide-Angle Seismic Attributes Using the Gaussian Beam Method

In this study, observed seismic attributes from shot gather 11 of the SAREX experiment are used to derive a preliminary velocity and attenuation model for the northern end of the profile in southern Alberta. Shot gather 11 was selected because of its prominent Pn arrivals and good signal to noise ratio. The 2-D Gaussian beam method was used to perform the modeling of the seismic attributes including travel times, peak envelope amplitudes and pulse instantaneous frequencies for selected phases. The preliminary model was obtained from the seismic attributes from shot gather 11 starting from prior tomographic results. The amplitudes and instantaneous frequencies were used to constrain the velocity and attenuation structure, with the amplitudes being more sensitive to the velocity gradients and the instantaneous frequencies more sensitive to the attenuation structure. The resulting velocity model has a velocity discontinuity between the upper and lower crust, and lower velocity gradients in the upper and lower crust compared to earlier studies. The attenuation model has Q _p ^-1 values between 0.011 and 0.004 in the upper crust, 0.0019 in the lower crust and a laterally variable Q _p ^-1 in the upper mantle. The Q _p ^-1 values are similar to those found in Archean terranes from other studies. Although the results from a single gather are non-unique, the initial model derived here provides a self-consistent starting point for a more complete seismic attribute inversion for the velocity and attenuation structure.