基于Sentinel-1A的2021年双湖MW5.7地震同震破裂模型

利用Sentinel-1A SAR数据提取2021年西藏双湖县M_W5.7地震同震形变场及2.5D形变场,反演断层滑动分布模型。计算不同节面解为接收断层产生的库仑应力变化差异确定发震构造,并结合余震分布信息评估未来地震风险性。结果表明,地震震中为34.37°N、87.71°E,震源深度6.51 km,发震断层倾向东、走向33°、倾角50°、平均滑动角-74°,以倾滑为主兼有少量左旋走滑分量,最大滑动量0.26 m。短时间内,震区南部地震风险较小,北部则需要结合更多资料进一步分析。本次地震是在羌塘块体持续向东扩张的背景下,受EW向拉伸作用使得黄水湖正断层发生的一次弥散型变形活动,SN向地堑得到进一步扩张。     

衍射地震CT技术的数值研究

本文首先解决了声波方程的非Born近似的正演计算问题,从而获得理论上不带近似的正演数据;然后,推导了井间（CBP）、垂直地震剖面（VSP）和地面反射（SRP）三种不同的数据采集方式下的衍射CT的重建公式;利用这些重建算法和正演数据,系统地研究了影响到地球物理CT成象质量的三种因素,即:（1）数据采集方式,（2）异常程度和（3）成象区域的尺寸,对重建图象的影响;并比较了衍射地震CT和射线地震CT的成象质量。     

双侧向电阻率测井反演算法研究

双侧向测井仪因其激励源的复杂性给其测量数据的反演带来了极大的困难．本文利用变形玻恩（Born）迭代法对轴对称二维非均匀介质中双侧向测井仪的位场数据进行了严格反演．在每一次迭代过程中采用了高效数值模式匹配的正演计算方法,它可半解析地求出反演中所需计算的格林（Green）函数及其偏导数,同时利用所推导出的位场非线性积分方程中积分运算的半解析形式,大大提高了计算效率和反演质量．通过反演实例表明,反演结果与电导率原状剖面吻合较好,并且受噪声影响较小．     

Generalized Born inversion of seismic reflection data

Summary . Born inverse methods give accurate and stable results when the source wavelet is impulsive. However, in many practical applications (reflection seismology) an impulsive source cannot be realized and the inversion needs to be generalized to include an arbitrary source function. In this paper, we present a Born solution to the seismic inverse problem which can accommodate an arbitrary source function and give accurate and stable results. It is shown that the form of the generalized inversion algorithm reduces to a Wiener shaping ***filter, which is solved efficiently using a Levinson recursion algorithm. Numerical examples of synthetic and real field data illustrate the validity of our method.     

On the spatial distribution of seismicity and the 3D tectonic stress field in western Greece

We analyzed a large number of focal mechanisms and relocated earthquake hypocenters to investigate the geodynamics of western Greece, the most seismically active part of the Aegean plate-boundary zone. This region was seismically activated multiple times during the last decade, providing a large amount of enhanced quality new information that was obtained by the Hellenic Unified Seismological Network (HUSN). Relocated seismicity using a double-difference method appears to be concentrated above ∼35 km depth, exhibiting spatial continuity along the convergence boundary and being clustered elsewhere. Earthquakes are confined within the accreted sediments escarpment of the down-going African plate against the un-deformed Eurasian hinterland. The data arrangement shows that Pindos constitutes a seismic boundary along which large stress heterogeneities occur. In Cephalonia no seismicity is found to be related with the offshore Cephalonia Transform Fault (CTF). Onshore, NS crustal extension dominates, while in central and south Peloponnesus the stress field appears rotated by 90°. Shearing-stress obliquity by 30° is indicated along the major strike-slip faults, consistent with clockwise crustal rotation. Within the lower crust, the stress field appears affected by plate kinematics and distributed deformation of the lower crust and upper mantle, which guide the regional geodynamics.     

青藏高原及邻近地区地壳与上地幔剪切波三维速度结构

本文利用30个基准台所记录的238条长周期面波资料,经过适配滤波和分格频散反演,得到中国大陆及邻区147个分格10—105s的纯路径频散,进而反演出青藏高原及邻近地区深至170km的剪切波三维速度结构.研究表明,青藏高原中西部地区和东部地区的地壳平均厚度分别为70±7km和65±7km,地壳平均剪切波速度分别为3.55和3.62km/s,上地幔顶盖平均速度分别为4.63和4.61km/s; 岩石层厚度均为120±10km;东部地区下地壳内30—40km深度处普遍存在低速层;青藏高原及其东侧的上地幔低速层内有横贯东西且明显向上隆起的低速腔.滇西缅北地区的地壳厚45±5km,上地壳及下地壳内都有低速层;上地幔顶盖的速度为4.42km/s,比青藏高原本体及恒河平原都低.恒河平原地壳厚34±2km,速度平均为3.45km/s;上地幔顶盖厚86±10km,速度平均为4.63km/s,顶盖内55—83km深处有一个低速夹层.     

青藏高原大尺度地表流体的重力效应特征分析

基于2002~2017年三维大气再分析模型、水文模式及GRACE卫星重力数据,采用负荷格林函数和卫星重力反演方法分析青藏高原大尺度地表流体的重力效应特征。结果表明,青藏高原地区大气扣除季节项前后的重力效应峰对峰变化介于1.9~10.7 μGal和1.1~4.6 μGal,其负荷效应呈现出显著的季节性特征;土壤水和积雪扣除季节项前后的重力效应峰对峰变化介于0.9~9.9 μGal和0.7~8.3 μGal,其负荷效应呈现出显著的季节性和年际变化特征;GRACE时变重力场扣除季节项前后的重力变化峰对峰变化介于4.4~24.1 μGal和3.4~18.3 μGal,考虑观测误差、构造运动和信号泄漏误差的影响,其显著的季节性和年际变化反映了多种地表流体的综合影响。     

三河-平谷8.0级地震区浅层三维密度结构反演研究

为研究三河-平谷地震区浅层构造背景、地震孕育机理及地震与构造的关系,利用高精度重力异常数据,采用基于块体生长模式的重力三维反演算法对地震区浅层三维密度结构进行反演,并通过模拟试算验证基于块体生长模式反演方法的有效性和稳定性。高精度布格重力异常显示,三河-平谷8.0级地震位于大兴重力局部高、三河-马坊重力局部高与大厂重力低之间的交汇过渡低值区域。研究区浅层三维密度结构反演结果表明,1679年三河-平谷M8.0地震明显受NE向夏垫断裂控制,断裂两侧密度差异明显且向下延伸约10 km,推测发震部位深约10 km。     

基于非结构化有限元的三维井地电阻率法约束反演

电磁探测反演是典型的不适定问题,易造成反演结果的多解性,不适定性是反演自身固有的特征,没有求解的附加信息这一本质困难是很难克服的,解决该问题的有效方法是研究约束反演。本文采用目前较为主流的高斯牛顿—共轭梯度法(GN-CG),在反演目标函数中直接施加约束条件,将介质电阻率的取值范围作为先验信息和约束条件以外点罚函数法的方式引入到反演目标函数中,与常规三维电阻率反演目标函数相比,增加了不等式约束项的目标函数,理论上可以压制反演的多解性。通过多种理论模型的测试结果表明,本文基于不等式约束的三维井地电阻率反演算法有效地改善了反演结果的精度,以惩罚函数法施加不等式约束条件的方式是现实可行及有效的。