昆明地震台地面与井下地震观测背景噪声对比

对昆明地震台山洞内放置的宽频带地震仪BBVS-60及观测井内放置的GL-S60B深井宽频带地震仪所记录的资料,进行背景噪声对比分析,得到基岩地区地面及井下背景噪声干扰,为今后在云南地区推广井下观测提供参考依据。结果显示,在1—20 Hz频段内,井下台背景噪声RM_S值比地表台小,降噪效果UD向达30%以上,EW向达50%,NS向达60%。在100—30 s频带内噪声不降反升,原因有待进一步分析。     

甘肃省地震局高性能计算系统

高性能计算是现代科学研究、工程技术开发和大规模数据处理的关键支撑技术,在地震系统发挥着日益重要的作用。甘肃省地震局高性能计算集群系统由60个节点、采用40 GB InfiniBand网络互连组成,峰值计算速度达9.2万亿次/秒。作为地震系统重要的科技服务平台之一,为地震反演、数据分析、数值模拟等领域以及科研创新平台提供高性能计算环境。     

Three-dimensional forward modeling of DC resistivity using the aggregation-based algebraic multigrid method

To speed up three-dimensional (3D) DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method (AGMG). We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner (AGMG-CG) to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling (3DDCXH). In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.     

A method of reconstructing complex stratigraphic surfaces with multitype fault constraints

The construction of complex stratigraphic surfaces is widely employed in many fields, such as petroleum exploration, geological modeling, and geological structure analysis. It also serves as an important foundation for data visualization and visual analysis in these fields. The existing surface construction methods have several deficiencies and face various difficulties, such as the presence of multitype faults and roughness of resulting surfaces. In this paper, a surface modeling method that uses geometric partial differential equations (PDEs) is introduced for the construction of stratigraphic surfaces. It effectively solves the problem of surface roughness caused by the irregularity of stratigraphic data distribution. To cope with the presence of multitype complex faults, a two-way projection algorithm between threedimensional space and a two-dimensional plane is proposed. Using this algorithm, a unified method based on geometric PDEs is developed for dealing with multitype faults. Moreover, the corresponding geometric PDE is derived, and an algorithm based on an evolutionary solution is developed. The algorithm proposed for constructing spatial surfaces with real data verifies its computational efficiency and its ability to handle irregular data distribution. In particular, it can reconstruct faulty surfaces, especially those with overthrust faults.     

华北克拉通地壳结构及动力学机制分析

本文对布设在华北克拉通三个陆块的199个宽频带台站记录的远震数据进行了接收函数计算.利用H-κ迭代方法获得了该区域基岩地区的地壳结构,平滑处理后作为背景结构模型中的基岩地区地壳结构;利用相邻算法对沉积层地区的接收函数进行了波形拟合计算,获得了沉积层结构,平滑后作为背景结构模型中的沉积层结构;结合前人的研究成果,完善了研究区域的背景结构模型.以此模型为基础,对接收函数进行了CCP（Common Conversion Point,共转换点）叠加成像,获得了Moho面成像结果,对比沉积层的成像结果发现：西部陆块中鄂尔多斯块体东部地区地壳厚度较大,约为42 km,泊松比较低,小于0.24,为长英质含量较多的地壳层;位于中部陆块的山西地堑地壳厚度小于鄂尔多斯块体,且变化较大,西侧地壳厚度约为40 km,东侧重力梯度带附近地壳厚度迅速减薄至36 km左右,张家口-怀来-大同一带出现了地壳的局部抬升,地壳厚度等值线基本以北北东方向为主,与构造带方向基本一致,地堑内泊松比约为0.26~0.28,前人对此区域的层析成像研究结果表明太行山隆起和阴山隆起存在壳内低速层,推测为地壳部分熔融以及上地幔物质上涌造成的;东部陆块中渤海湾盆地的地壳厚度较薄,约为32 km,部分地区小于30 km,其中冀中坳陷带地壳厚度最薄,约为28 km,沉积层基底分布与Moho面分布呈镜像对称趋势,沉积层较厚地区的地壳较薄,推测东部陆块在太平洋板块俯冲作用下,存在北西-南东向的拉张作用,使其内发育了大量断陷盆地.     

利用GPS和GRACE分析四川地表垂向位移变化

陆地水储量的季节性变化是导致地表周期性负荷形变位移的主要因素,有效地剔除地表位移中的陆地水储量影响,是获取地壳构造垂向运动的必要过程.四川地处青藏高原东边缘,地形分区明显,境内以长江水系为主,水资源丰富,研究四川地区地表负荷形变位移,有助于分析陆地水储量的时空分布特性及地壳构造形变信息.本文利用研究区域内59个CORS站的GPS观测数据,计算了CORS站点的垂向位移,并将其与GRACE所得相应结果进行对比分析.结果显示,GPS和GRACE所得垂向位移时间序列的振幅大小整体相符,但存在明显的相位差.GPS站点振幅最大值为12.7 mm,对应HANY站,最小值为1.5 mm,对应SCMX站.GRACE所得的地表垂向位移振幅大小均为3~4 mm,且最大位移集中出现在7-9月份;而GPS站点出现最大位移的月份和地形相关,东部盆地、西北部高原和南部山地分别出现在7-8月份、10-11月份和10月份.GPS站点时间序列中的周年项与陆地水的季节性变化强相关,为了讨论陆地水储量对GPS站点位移的影响,本文利用改进的总体经验模态分解方法（MEEMD：Modified Ensemble Empirical Mode Decomposition）,从GPS垂向位移时间序列中提取出周年项及约2年的年际变化项.发现利用MEEMD获取的周年项改正原始GPS时间序列时可使其WRMS（Weight Root Mean Square）减少量减小约26%,结果优于最小二乘拟合方法提取的GPS周年项改正效果,验证了MEEMD方法在GPS坐标时间序列处理中的可行性及有效性.     

Variation in reach‐scale thalweg‐migration intensity in a braided reach of the lower Yellow River in 1986–2015

Thalweg migration of an alluvial river plays a key role in channel evolution, which may influence the effect of existing river training works and biodiversity on floodplains, and cause losses in riparian land and property. The braided reach of the Lower Yellow River underwent continuous channel aggradation during the period from 1986 to 1999, and then remarkable channel degradation in 1999–2015 owing to the state of operation of the Xiaolangdi Reservoir in 1999. Here we quantify associated thalweg migration changes and identify the key influencing factor in the braided reach. Thalweg‐migration distances and intensities at section‐ and reach‐scales were calculated during the past 30 years from 1986 to 2015, in order to investigate the characteristics of thalweg migration in the reach. There was a 47% reduction in the reach‐scale thalweg‐migration distance and a 35% reduction in the corresponding migration intensity after the reservoir operation. It is also revealed that fluvial erosion intensity is a dominant factor in controlling the thalweg migration, based on the investigation into various influencing factors in the study reach. The thalweg‐migration intensity of the braided reach can be expressed as a power function of the previous four‐year average fluvial erosion intensity. The calculated thalweg‐migration intensities in 1986–2015 using the proposed relation generally agree with the observed data. Copyright © 2017 John Wiley & Sons, Ltd.     

宾川地区的地壳厚度和泊松比分布特征研究

This paper presents the changes of crust thickness and Poissons ratios distribution in the Binchuan region,where the first air-gun transmitting station and its a small dense array were deployed.From September 2011 to January 2014,more than 239 teleseismic events of M≥6.0 were recorded in 16 stations in the Binchuan region.Their P-wave receiver functions were analyzed respectively.The high spatial resolution result shows that the average crust thickness of Binchuan region is 45.3km,it follows the rule of"deeper in the north and east part,shallower in the south and west part."The deepest region is in Xiaoyindian Station;the crust thickness is 47.9km;the shallowest region is in Paiying Station,it has the thickness of 42.1km.It shows that the deeper Moho surface nearby the Chenghai fault and shallower nearby the Honghe fault;the isoline distribution of thickness changes greatly nearby the Chenghai fault and slowly nearby the Honghe fault.From the distribution of Poissons ratios,it is unevenly in the study area with a great difference from the north part to the south part,which shows a characteristic of"lower in the north,higher in the south".The Poissons ratio nearby the Honghe fault is medium too high(0.26≤σ≤0.29);lower nearby the Chenghai fault(σ≤0.26).This paper concludes the possible reason of different characteristic between Poissons ratio and crust thickness is thicker in the upper crust in the Binchuan region.     

不同尺度裂缝对弹性波速度和各向异性影响的实验研究

裂缝广泛分布于地球介质中并且具有多尺度的特点,裂缝尺度对于油气勘探和开发有着重要的意义.本文制作了一组含不同长度裂缝的人工岩样,其中三块含裂缝岩样中的裂缝直径分别为2 mm、3 mm和4 mm,裂缝的厚度都约为0.06 mm,裂缝密度大致相同(分别为4.8%、4.86%和4.86%).在岩样含水的条件下测试不同方向上的纵横波速度,实验结果表明,虽然三块裂缝岩样中的裂缝密度大致相同,但是含不同直径裂缝岩样的纵横波速度存在差异.在各个方向上,含数量众多的小尺度裂缝的岩样中纵横波速度都明显低于含少量的大尺度裂缝的岩样中纵横波速度.尤其是对纵波速度和SV波速度,在不同尺度裂缝岩样中的差异更明显.在含数量多的小尺度裂缝的岩样中纵波各向异性和横波各向异性最高,而含少量的大尺度的裂缝的岩样中的纵波各向异性和横波各向异性较低.实验测量结果与Hudson理论模型预测结果进行了对比分析,结果发现Hudson理论考虑到了裂缝尺度对纵波速度和纵波各向异性的影响,但是忽略了其对横波速度和横波各向异性的影响.     

汶川地震科学钻探3号井孔附近微震震源参数确定及其意义

基于微震监测仪器频带及波形记录的时频特征,我们筛选出汶川地震科学钻探3号井孔周围微震台阵2012年记录的218个t_s-t_p1 s的微震.通过盖戈法与和达法相结合确定微震震源的几何参数,发现这些微震分布大体呈NE-SW展布,与龙门山断裂带的走向基本一致.无论是采用《地震台站观测规范》中的量规函数(量规函数GF)还是李学政等(2003)的量规函数(量规函数LXZ),近震震级均与矩震级呈现出较好的线性关系,向震级小的一端延伸时都表现为M_LM_W,但采用李学政等的量规函数时该趋势更加明显.同时采用Brune和Boatwright震源谱衰减模型对观测震源谱的拟合表明拐角频率具有模型依赖性:基于Brune模型拟合得到的拐角频率大于基于Boatwright模型得到的拐角频率.无论基于哪一种模型,矩震级与拐角频率、破裂半径的对数线性关系均较弱,与应力降、视应力的对数则存在较好的线性关系.这些关系不支持应力降、折合能量、视应力为常数的观点,表明微震与大地震的震源物理过程存在差异.视应力与应力降成比例,比例系数小于0.5,表明破裂动力学模式符合SavageWood模式.近震震级与矩震级拟合关系Mw=a+bM_L中b的大小与应力降和地震矩的关系有关,△σ∝M_0~γ,则b=1/(1+γ),因此从b的大小可以粗略地判断应力降与地震矩的关系.本文对应于γ≈1的情况,与基于量规函数LXZ得到的近震震级与矩震级关系中系数b=0.53吻合,这说明仅从辐射能量的角度考虑量规函数LXZ较量规函数GF准确.