基于复模态的有限元模型修正算法

针对地下结构地震响应分析中无限地基辐射阻尼问题，引入复模态情况下的具有非简化的堆积阻尼矩阵的阻尼模型，并针对具有集中质量阵的阻尼模型提出了合并与质量有关的阻尼和堆积阻尼的思想，并据此提出了一种修正此类有限元模型的两步法，首先从复模态参数中提取实模态参数，采用基于模态残余力的识别算法修正刚度矩阵，然后根据复模态参数和已得的刚度矩阵来识别阻尼模型中的刚度参与系数和质量阻尼堆积阻尼联合矩阵。     

大地测量反演问题统一综述

大地测量反演在解决人类面临的资源短缺、自然灾害和生态环境退化问题显示出越来越重要的作用和地位。本文根据R.Parker的理论,首次尝试将近三十年来大地测量反演问题从解的存在性、模型构制、解的非唯一性、解的评价四个方面进行统一,促进了大地测量反演理论和应用的研究。     

基于ETOPO1数据的全球及重点海域水深统计与分析

由于不同的水深探测手段及海深反演方法适用于不同深度的海域,且全球海域水深分布不均,因此在探测及反演海深之前,需对目标海域的水深分布有一定认识。本文采用ETOPO1模型对全球及部分海域进行统计分析,结果表明,全球及大西洋、印度洋、太平洋区域的水深大多处于3 000—10 000 m,其中3 000—5 000 m水深海域占比较大,因此,在这些海域主要采用重力数据结合船测数据进行反演,同时应注意避开海岸带、岛屿较多、海脊存在的区域;对于水深小于1 000 m的海域,可以根据外部环境,采用多波束船测、水深遥感、机载激光测深等方式进行探测。通过分析及统计全球和重点海域的水深可以为海深探测及海底地形反演方法的选取提供参考。     

利用地震面波和重力资料联合反演地壳—上地幔三维密度结构的方法探讨

本文探讨了利用地震面波和重力资料联合反演地壳－地幔三维密度结构的反演问题。首先建立了地震面波和重力资料的观测方程,然后应用广义线性反线理论给出了反演问题的解。模型试验结果表明,与单纯面波反演的结果相比,联合反演的解在分析率和方差两个方面都有改善,当加入扰动重力数据时,在０－３００ｋｍ深度范围内联合反演的结果明显好于单纯面波反演的结果。     

基于城市地表温度与归一化植被指数对北京市建成区热环境分析

针对研究城市热环境的过程中,利用归一化植被指数(NDVI)进行地表温度(LST)反演,再将LST和NDVI结合说明地物变迁与城市热环境的影响的现状,利用Landsat-8多时相遥感影像、高分辨率影像、公开GIS等多源数据,通过人工交互判读和量化统计分析,实现了2013—2017年北京建成区NDVI变化及其对地表热环境影响分析,再对分析结果进行差值拟合评价。对NDVI阈值分割按照大小为LC1、LC2、LC3、LC4,对LST分为高温区(TH)、常温区(TR)和低温区(TL)。结果表明,2013年至2017年:1)建成区的平均NDVI增加0.03,其中LC1增加1.0%,LC2减少11.6%,LC3区域减少1.7%,LC4区域增加12.3%。2)建成区平均LST增加2.55 K,TH百分比增加0.6%;TR百分比减少1.1%,TL百分比增加0.5%。3)NDVI相对增加,地表温度相对下降以及NDVI相对减少,地表温度相对上升占60%,NDVI相对下降,地表温度相对下降以及NDVI相对增加,地表温度相对上升的占40%。     

海斗深渊区十字放炮与POBS精确位置校正

TS03航次创新性地投放了 6台被动源海底地震仪(POBS)到马里亚纳海沟南部"海斗深渊区",实施了万米级水深的天然地震观测,通过十字放炮获得POBS的精确坐底位置是后续研究工作中至关重要的一环.本文首先对POBS的SEG-Y数据进行裁截误差修复,然后采用直达水波走时反演的方法得到这6个POBS的精确坐底位置,并进行精...     

Imaging the ramp–décollement geometry of the Chelungpu fault using coseismic GPS displacements from the 1999 Chi-Chi, Taiwan earthquake

We use coseismic GPS data from the 1999 Chi-Chi, Taiwan earthquake to estimate the subsurface shape of the Chelungpu fault that ruptured during the earthquake. Studies prior to the earthquake suggest a ramp–décollement geometry for the Chelungpu fault, yet many finite source inversions using GPS and seismic data assume slip occurred on the down-dip extension of the Chelungpu ramp, rather than on a sub-horizontal décollement. We test whether slip occurred on the décollement or the down-dip extension of the ramp using well-established methods of inverting GPS data for geometry and slip on faults represented as elastic dislocations. We find that a significant portion of the coseismic slip did indeed occur on a sub-horizontal décollement located at 8 km depth. The slip on the décollement contributes 21% of the total modeled moment release. We estimate the fault geometry assuming several different models for the distribution of elastic properties in the earth: homogeneous, layered, and layered with lateral material contrast across the fault. It is shown, however, that heterogeneity has little influence on our estimated fault geometry. We also investigate several competing interpretations of deformation within the E/W trending rupture zone at the northern end of the 1999 ground ruptures. We demonstrate that the GPS data require a 22- to 35-km-long lateral ramp at the northern end, contradicting other investigations that propose deformation is concentrated within 10 km of the Chelungpu fault. Lastly, we propose a simple tectonic model for the development of the lateral ramp.     

北京大灰厂试验场重力变化特征及其解释

介绍了北京大灰厂试验场（简称京西试验场）的概况，对１９９３年３月以来重力重复观测资料进行了处理。计算了各期重力观测精度，结果显示１９９３年３月～１９９５年５月京西试验场存在一定重力变化。采用给定的模型反演了八宝山断层介质的密度变化，认为断层破碎带介质密度的变化可能是引起地表重力变化的主要因素     

Modelling and inversion of single-ended refraction data from the shot gathers of multifold deep seismic reflection profiling—an approach for deriving the shallow velocity structure

A multifold crustal-scale deep seismic near-vertical reflection profile generates a large number of single-ended shot gathers, which provide redundant data sets because of overlapping coverage of the shallow refractors. We present an approach for deriving the shallow velocity structure by modelling and inversion of single-ended seismic refraction first arrival traveltime data. We apply this method to a data set acquired with a 12-km long spread with 100 m spacing of shots and receivers, of the Neoproterozoic Marwar basin in the NW Indian shield. The approach is shown to be quite successful for delineating the shallow refractor depths, steep dips and velocities, even in the absence of regular reverse refraction profiles. The study reveals two-layered sedimentary formations, Malani volcanics and a complicated basement configuration of the Marwar basin, and provides a measure of resolution and uncertainty of the estimated model parameters. A seismic section of the near-trace gather is found to be qualitatively consistent with the derived structural features of the basin. The relative highs and lows, observed in the Bouguer gravity profile, further corroborate the derived velocity model. The present approach can be especially useful in offshore areas and elsewhere, where the single-ended multifold seismic profiles are the only available data sets.