Low-frequency 3-D seismic survey of upper crustal magmatic intrusions in the northeastern Pannonian basin of Hungary

A specially designed 700-km² grid survey, deploying 1000 regularly distributed low-frequency seismic recording systems, successfully investigated one of the most complex geologic environments of the Pannonian basin. The wide-angle signals penetrated through over 1000 m of multi-phase igneous lithology and recognized, for the first time, the underlying enigmatic Permian to Early Triassic basement rocks. Tomographic inversion of the first arrival grid data resulted in determination of an accurate three-dimensional (3-D) velocity field, to a depth of 4 km. The anomalous changes of the spatial velocity data outline the regional extent of the Late Miocene magmatic intrusions, which are covered by over 2000 m of Mid-Miocene to Pleistocene clastics. Complex relationship was found between the surface potential data and the intrusive bodies. This multi-faceted geophysical data analysis established a functional technique for mapping a subsurface with intricate acoustic and structural complexity.     

TRMM星载测雨雷达和地基雷达反射率因子数据的三维融合

TRMM卫星上的测雨雷达(TRMMPR)探测资料分布均匀且具有很高的垂直分辨率,但灵敏度较低;地基雷达(GR)水平分辨率较高且具有较高的灵敏度,但其垂直分辨率低。通过将TRMM PR与GR反射率因子数据的三维数据融合,得到了更优的反射率因子图像。测雨雷达与地基雷达三维数据融合主要分为以下几步:测雨雷达与地基雷达数据预处理——如去杂波、衰减校正;测雨雷达与地基雷达时空匹配;选取和应用合适的三维图像融合算法;对融合后的图像进行效果评估。试验结果表明:融合后的图像不仅增大了信息量,更好地检测弱降水,还提高了空间三维(3D)分辨率,能更好地反映降水区域细节,且使得数据总体上具有更高的完整性和可靠性。此外,还将基于雷达资料估测的降水数据与地面雨量计数据进行对比,估计反射率因子数据融合在降水测量上的有效性。      

月球构造新认识：基于60m高精度全月月形图的月球构造形态学研究

一直以来月球上都没有发现像地球上一样的全球板块构造现象,被认为是属于单板块构造的行星体并且是不活动的。随着月球地形探测精度的提高,人们逐渐发现了一些地质现象表明月球并不像之前想象的一样是不活动的。高精度的月形数据显示出了月球上除了存在众所周知的(撞击作用产生的)复杂的撞击构造外,仍然存在许多可能由内生作用产生的构造特征:如裂谷、山脊、火山群等。本文利用最新高精度月球激光高度数据(LOLA数据)建立了全月60m(相当于7秒网格)高精度月形图。由于数据量超大(TB级别),作者编写了基于并行架构的集群计算处理代码,其计算效率提高了近300倍。基于这些图件,作者在全球范围内找出许多重要的构造现象,并根据其构造及地质特征进行分类研究。研究结果表明并不像之前认识一样,月球上的确存在许多地质构造体是由月球的内生地质作用所产生,并建立对于月球结构和构造的新认识。此外,通过对比发现,像地球上全球尺度的线性山脉只是板块构造运动的结果;另一方面,在无板块构造的行星上,就像月球一样,会产生大面积的高地(月陆)和低地(月海地区),这种地质现象很有可能是行星地壳岩浆演化的直接结果。在文章最后,作者还将制作出的60m分辨率的全月月形系列图件附后并公开,以方便其他科研人员进一步开展相关研究;并以此敬献中国地质科学院60周年庆典。     

球面距离计算方法及精度比较

球面距离(角间距)计算是天文或地理学中极常用的计算之一,也是目标查找、锥形检索、交叉证认等方法的基础。数学上,通过球面几何可以直接计算出两点的距离,前人已经推导出了多个复杂程度不一的计算方法。但是由于计算机的精度有限,在进行数值计算时有舍入误差,导致公式计算结果出现偏差。对几个常用的球面距离计算公式进行了考察,测试并对比它们在不同计算环境下的精度与优缺点。此外还展示并比较了几种常用天文软件包、数据库的球面距离计算方法,以期有助于天文工作者选择适合自己当前需要的计算方法。     

九寨沟Ms7.0级地震三维同震形变场提取

针对D-InSAR地震形变监测中存在LOS(line of sight)视线向模糊的问题,该文构建了一种融合升降轨不同视线向干涉测量数据获取三维同震形变场的方法。以2017年九寨沟M_S7.0级地震为例,基于研究区的Sentine1_1A数据和SRTM(1″)DEM数据,采用二轨差分的D-InSAR技术,融合升降轨LOS向以及自定义视线向的干涉测量数据进行联合解算,获取了研究区三维同震形变场,解决了单一轨道雷达LOS提取结果不能准确反映地表三维形变场的问题。实验表明,此次地震造成了震中附近一定范围地表的沉降和东南向的滑移,并通过震区形变剖面图和形变等值线图,分析了形变场的空间分布以及此次地震所造成地表断裂的位置,认为此次地震与塔藏断裂、虎牙断裂存在关联性,其运动形式为主动盘逆冲的走滑型地震。     

科学计算程序库WHIGG F90 LIB及其地学应用

WHIGG F90LIB(WFL)是一款用FORTRAN90语言编写的用于地球学科的计算程序库。应用WFL可以减少重复性编程,方便地解决科学计算中数据的预处理问题,结合多种分析方法,对科学问题进行定量分析。利用WFL在地球重力场、气候涛动和数字信号处理3个方面的应用实例,以验证WFL在科学计算中应用的有效性。     

地震层析成像LSQR算法的并行化

讨论了地震层析成像的LSQR算法(最小二乘QR分解). 在建立偏导数矩阵方程组时,对区内地震在方程中保留震源项,引入正交投影算子进行参数分离,对区外远震采用传统的平滑处理方式,用LSQR法求解联立的方程组. 由于区内地震的正交分解处理和区外远震的平滑处理,使得偏导数矩阵中的非零元素成倍增加,对于大型反演问题,这些非零元素常常达到几十GB到几百GB的数量级,巨量的内存占用成为LSQR算法的瓶颈. 针对这一问题,本文研究了偏导数矩阵中非零元素的分布规律,设计出合理的存储结构,采用分布式存储进行矩阵计算,提出了LSQR算法的并行化方案,并在联想深腾6800超级计算机上实现. 导出了LSQR算法的并行效率估算公式. 对两个地区的实际地震层析成像数据进行了效率测试.     

多视密集匹配并行传播GPU-PatchMatch算法

针对多视密集匹配的效率较低的问题,提出了GPU-PatchMatch多视密集匹配算法。该算法使用GPU提高PatchMatch的计算效率;同时充分利用稀疏场景信息,对深度信息进行规则初始化;为提高传播效率,使用了金字塔红黑板并行传播深度信息。最后在DTU、Strecha和Vaihigen数据集上进行了试验,并与常用的多视密集匹配算法进行对比。试验结果表明,本文算法在重建效率上有较大提高,与CPU算法(PMVS、MVE、OpenMVS)相比有7倍以上提升,与GPU算法相比也有2.5倍以上提升,表明本文算法的有效性。     

胶州湾三维变动边界的潮流数值模拟

基于Blumberg等（1986）的河口、陆架和海洋模式，引入变边界处理技术，建立了胶州湾三维变动边界的潮流模型，模型以σ-正交曲线坐标下三维非线性潮波方程为基本方程，引入2.5阶瑞封闭方程组，采用分裂算子法数值求解方程组，利用湍封闭方程求解垂直紊动粘滞系数，采用干湿网格方法模拟潮流漫滩过程，三维变动边界潮流模型计算结果与实测值吻合良好。     

基于MPI的LAGFD-WAM海浪数值模式并行算法研究

在LAGFD-WAM海浪串行数值模式基础上,利用MPI信息传递机制实现其并行化。通过对模拟区域合理划分,对数据采取分块加载,实现了各个节点的负载平衡;通过对算法的改进实现了粗粒度计算,大幅度减少了通信量,从而提高了程序的执行效率。对串行计算和并行计算的效率比较表明,本文建立的方法能够得到较高的加速比。对全球海浪模式,加速比和CPU数目大致呈线性关系。对高分辨率的区域海浪模式,在128 CPU条件下加速比可以达到91.9。