低渗透压力条件下砂岩渗透性质的CT试验

运用岩石高压三轴加载装置和渗透压加载装置,对砂岩进行固定应力变渗透水压力试验,同时借助SOMATOMPLUS螺旋CT扫描机进行实时观测。通过试验结果分析,推出了基于CT数的岩石空隙率公式,从CT尺度上分析了岩石空隙率、渗透速度、渗流速度、微孔隙直径、渗透率等随渗透水压力的变化规律。结果表明:岩石的渗透参数随渗透压力的增加而增大,渗透参数随渗透压力的变化呈对数关系,这是由于渗透水压力使岩石内部空隙发生变形,这一变形为弹性变形;而渗流速度与渗透水压力呈线性关系,符合宏观尺度渗流的达西定律。      

福建省地倾斜固体潮特征及其动力学分析

利用福建地区厦门、漳州、莆田、南平、龙岩5个台站地倾斜固体潮观测资料,分别计算了地震平静期的倾斜固体潮潮汐变化特征量,并进行了对比分析;除此之外还对固体潮特征进行了地球动力学分析。结果表明,沿海地区海潮是决定倾斜固体潮特征的主要因素。     

湖泊磷的核磁共振(NMR)分析——1环境标准样品的31P-NMR表征

通过对环境标准样进行31P-NMR测试,研究了磷的化学结构与化学位移的关系.结果表明:无机正磷酸盐(PO43-)化学位移在6.412×10-6,焦磷酸盐的化学位移在-4.353×10-6;多聚磷的化学位移分布特点是:主链末端磷化学位移在-3.53×10-6至-4.40×10-6之间,中间磷在-18.159×10-6至-20.246×10-6之间,有机磷中单酯和二醑的化学位移分别在4.3540×10-6至6.3652×10-6和-0.107×10-6至2.195×10-6之间,部分二脂不稳定,会降解为单脂;膦酸盐的化学位移在18.670×10-6至21.460×10-6之间,磷酸肌酸位于-0.380×10-6,有机多聚磷中α-P的化学位移在-8.833×10-6至-9.449×10-6之间,由于偶合作用产生P峰的裂分,P峰的偶合符合核磁中n 1的规则.结构环境影响31P的化学位移值的偏移.     

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

对昆明地震台山洞内放置的宽频带地震仪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.