三维复杂山地多级次群推进迎风混合法多波型走时计算

三维复杂山地条件下的各种地震波型的走时计算技术,可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算各种地震波型的走时,提出一种多级次群推进迎风混合法.该算法利用不等距迎风差分法简洁稳定地处理三维复杂地表及附近的局部走时计算问题,利用计算精度不错的迎风双线性插值法处理绝大部分均匀正方体网格中的局部走时计算问题,利用群推进法模拟三维复杂地表条件下地震波前的扩展问题,利用多级次算法处理各种类型的地震波的走时计算问题.算法分析和计算实例表明:新方法具有很好的计算精度与效率,且能灵活稳定地处理三维复杂地表复杂介质条件下的多波型走时计算问题.     

甘肃地区一维速度模型计算研究

采用2009年1月1日—2014年6月甘肃省境内所有地震的震相数据,分类拟合计算该地区一维速度模型,计算结果显示,利用少量震级较大地震数据来拟合的速度结果近似于用所有地震拟合结果。拟合得出新模型初步结果后,再从所有地震资料中选定台站包围好、震相标识规范的142个地震,利用Hyposat定位方法来计算莫霍面深度,参考其他模型的莫霍面深度,综合分析后在一定的扰动范围内通过Hyposat批处理程序进行迭代定位,对8 640次计算模型做残差对比,最终选取残差最小模型作为新模型,此模型除部分参数有较小的差别外与甘青模型基本一致。     

杭州湾枯季水沙若干特征分析

根据2014年1月及2017年2月杭州湾大、小潮水沙资料,计算了流速、含沙量、单宽潮量、单宽输沙量,进而分析了潮周期断面净潮量、净输沙量和区域冲淤分布。研究发现杭州湾涨、落潮平均流速比值整体较20世纪八九十年代增大,绝大部分区域大于1.0,涨潮流相对增强,澉浦南岸和金山北岸尤为显著。含沙量平面上分布澉浦、杭州湾南岸两个高值区(大潮大于3 kg/m³,小潮大于1.3 kg/m³)以及北岸湾口至乍浦之间的低值区(大潮小于2 kg/m³,小潮小于0.9 kg/m³),随潮汛变化显著,最大含沙量浓度通常滞后于急流时刻。各测站涨、落潮量和输沙量呈现“大涨大落”和“大进大出”的特征,造成杭州湾短时间尺度内的“大冲大淤”。大潮两涨两落金山与乍浦、乍浦与澉浦之间区域净输沙量可达几千万吨,净冲淤则在几百万吨。     

Experimental study of aeolian sand ripples in a wind tunnel

The topographic parameters and propagation velocity of aeolian sand ripples reflect complex erosion, transport, and deposition processes of sand on the land surface. In this study, three Nikon cameras located in the windward (0–1 m), middle (4.5–5.5 m), and downwind (9–10 m) zones of a 10 m long sand bed are used to continuously record changes in sand ripples. Based on the data extracted from these images, this study reaches the following conclusions. (1) The initial formation and full development times of sand ripples over a flatbed decrease with wind velocity. (2) The wavelengths of full development sand ripples are approximately twice the wavelengths of initially formed sand ripples. Both wavelengths increase linearly with friction velocity. During the developing stage of sand ripples, the wavelength increases linearly with time. (3) The propagation velocity of full development sand ripples is approximately 0.6 times that of the initially formed sand ripples. The propagation velocity of both initial and full development of sand ripples increase as power functions with respect to friction velocity. During the developing stage of sand ripples, the propagation velocity decreases with time following a power law. These results provide new information for understanding the formation and evolution of aeolian sand ripples and help improve numerical simulations. Copyright © 2017 John Wiley & Sons, Ltd.     

DATABASE OF LANDSLIDES TRIGGERED BY 2015 GORKHA(NEPAL) MW7.8 EARTHQUAKE

In this study, a detailed database of landslides triggered by the 25 April 2015 Gorkha (Nepal)M_W7.8 earthquake is constructed based on visual interpretation of pre- and post-earthquake high-resolution satellite images and field reconnaissance. Results show the earthquake triggered at least 47 200 landslides, which have a NWW direction spatial distribution, similar with the location and strike of the seismogenic fault. The landslides are of a total area about 110km² and an oval distribution area about 35 700km². On the basis of a scale relationship between landslide area (A)and volume (V), V=1.314 7×A^{1.208 5}, the total volume of the coseismic landslides is estimated to be about 9.64×10⁸m³. In the oval landslide distribution area, the landslide number density, area density, and volume density were calculated and the results are 1.32km^-2, 0.31%, and 0.027m, respectively. This study provides a detailed and objective inventory of landslides triggered by the Gorkha earthquake, which provides very important and essential basic data for study of mechanics of coseismic landslides, spatial pattern, distribution law, and hazard assessment. In addition, the landslide database related to an individual earthquake also provides an important earthquake case in a subduction zone for studying landslides related to multiple earthquakes from a global perspective.     

Current Vertical Motion Analysis of Northwestern Margin of Ordos Based on Precise Leveling

Ordos block was squeezed by the Qinghai-Tibetan block and North China block, and the tectonic activity was intense. In the periphery of Ordos block, there was a series of folds zones and compressed faults with complicate structures. This paper used three-phase data of 1980,1990 and 2014 to calculate vertical velocity of Northwestern margin of Ordos and the analytical results indicated that ① the Hetao basin between the rise of Yinshan fault block and Ordos fault block showed relatively subsidence, in which Linhe basin was the most evident and the subsidence rate was about 2-4mm/a. The subsidence rate of Jartai -Yinchuan rift zone on the western margin of Ordos block was about 2mm/a; ② the whole testing zone exhibited the evident inherited movement characterized by mountain rise and basin subsidence; ③the two leveling section through the northern margin fault and Dengkou-Benjing fault showed that the difference between vertical velocities on two sides of the fault was less than 0.5mm/a.     

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.     

利用GPS数据研究南北地震带北段近期地壳水平形变特征

采用中国地震局"陆态网络"GPS基准站和流动站观测数据,利用MIT的开源软件GAMIT/GLOBK解算南北地震带北段(33°~42°N,97°~110°E),特别是六盘山地区、祁连地块、银川地堑1999年以来在欧亚框架下、以鄂尔多斯块体作为参考的GPS速度场。计算南北地震带北段5条断裂两侧50~350km范围内各站垂直/平行于断裂走向的速度投影,以及14个基准站18条基线2010年以来的长度时间序列,试图进一步估计关心区域的应力积累情况,探索断层闭锁段。结合地震活动性,认为南北地震带北段未来50年有发生M8地震的可能性;该段1900年以来的第4个M6、M7地震活跃期可能已经拉开序幕,未来2年有发生M7以上地震的可能性;祁连断裂带、海原断裂、西秦岭北缘地震带具备发生M6和M7地震的优先条件,尤其是海原断裂。     

福建及其邻区背景噪声面波群速度层析成像

基于福建及邻区108个宽频带地震台站2016年6月到7月两个月垂直和水平分量波形连续记录, 利用相位加权叠加算法提高信噪比, 计算得到108个台站对的高质量经验格林函数。 对所获取的经验格林函数, 采用时频分析的方法在1~20 s频段内量取了大量高质量的Rayleigh波和Love波群速度频散数据。 在此基础上, 采用基于射线追踪的二维层析成像方法反演得到了福建及其邻区1~20 s的Rayleigh波和Love波群速度分布。 分辨率测试结果表明群速度分布的分辨率能达50 km。 成像结果显示1~10 s的群速度分布与地壳中上部地质特征有很好的一致性, 区内的福州盆地和漳州盆地在浅层结构中表现出明显的低速异常。 长周期的群速度则揭示了漳州西北的高地热区内中下地壳低速体, 政和-大埔断裂两侧的速度差异, 表明其可能是一个深大断裂, 并呈现明显的东西差异。     

温州深软场地剪切波速实测误差及其对PGA影响研究

该研究在温州两个深软土层场地上完成了剪切波速测试误差专项实验,统计拟合了误差分布规律并研究了其对PGA的影响,得到结论:每米深度上的实测误差基本符合标准正态分布,温州地区深软场地的测试误差与我国其他几个地区覆盖土层较浅场地上的每米实测误差程度基本相同;两个场地上各深度的误差标准差大致相同,约为16%;实测误差对PGA影响明显,误差取1倍标准差时PGA的变化程度可能达到30%,取2倍标准差时可能高达45%;PGA的变化程度受输入地震动的频率特性制约,长周期频率成分对深软场地计算结果的影响明显。