添加正则化项的黏声逆时偏移成像方法研究

在实际地球介质中传播的地震波会产生衰减和频散现象,因此其更接近于黏弹性介质,在地震处理中补偿这些黏性影响是十分必要的。基于波动方程的叠前深度偏移中进行吸收衰减补偿更准确,也更有物理意义,直接求解双程波动方程的逆时偏移(RTM)能够成像大倾角复杂构造,具有诸多优势。然而当考虑吸收衰减补偿时通常会产生不稳定现象,大部分研究都是在逆时偏移的波场延拓中进行波数域的低通滤波来解决这个问题。本文采用广义标准线性固体的黏声波动方程进行吸收衰减补偿的Q--RTM方法,通过添加正则化项的方式来稳定延拓过程。添加正则化项本质上是低通滤波,滤波窗口是指数形式,在时空域有明确的形式,可以阻止发生高频不稳定。与直接滤波相比,正则化参数可以是空变的,因此比较适合剧烈变化的区域,我们还发现震源归一化互相关成像条件更适合Q--RTM方法。      

Seismic wavefield modeling based on time-domain symplectic and Fourier finite-difference method

Seismic wavefield modeling is important for improving seismic data processing and interpretation. Calculations of wavefield propagation are sometimes not stable when forward modeling of seismic wave uses large time steps for long times. Based on the Hamiltonian expression of the acoustic wave equation, we propose a structure-preserving method for seismic wavefield modeling by applying the symplectic finite-difference method on time grids and the Fourier finite-difference method on space grids to solve the acoustic wave equation. The proposed method is called the symplectic Fourier finite-difference (symplectic FFD) method, and offers high computational accuracy and improves the computational stability. Using acoustic approximation, we extend the method to anisotropic media. We discuss the calculations in the symplectic FFD method for seismic wavefield modeling of isotropic and anisotropic media, and use the BP salt model and BP TTI model to test the proposed method. The numerical examples suggest that the proposed method can be used in seismic modeling of strongly variable velocities, offering high computational accuracy and low numerical dispersion. The symplectic FFD method overcomes the residual qSV wave of seismic modeling in anisotropic media and maintains the stability of the wavefield propagation for large time steps.     

伸缩缝刚度对大跨度悬索桥动力特性的影响

伸缩缝作为大跨度桥梁与引桥之间的重要连接构件,其抗推刚度及可能存在的变异性对主桥及引桥动力特性的影响不可忽略。本文建立了大跨度悬索桥及引桥的有限元模型,采用弹簧单元模拟加劲梁与引桥箱梁之间的伸缩缝,分析伸缩缝刚度对悬索桥及引桥自振特性及其地震响应的影响规律。分析结果表明:伸缩缝刚度对加劲梁的横弯振型、竖弯与纵飘耦合振型的频率有明显的影响;伸缩缝刚度的变化会导致加劲梁与引桥的振型相互耦合,同时这些振型的频率发生相应的突变,当伸缩缝刚度较大时,加劲梁两个竖弯与纵飘的耦合振型解耦成为独立的竖弯和纵飘振型;当引桥与悬索桥加劲梁的纵飘振型发生耦合时,在纵向和竖向地震作用下的悬索桥及引桥的地震响应达到最小。伸缩缝刚度对悬索桥动力特性影响的分析可为悬索桥的模态参数确认、损伤识别、抗震性能分析提供有价值的借鉴。     

Stagnant low-energy ions in the near cusp region observed by Cluster

Measuring the low-energy ions in the Earth's magnetotail lobes is difficult, because a spacecraft becomes positively charged in a sunlit and tenuous plasma environment. Recent studies have introduced a new method, making use of the positive electric potential on the Cluster spacecraft, to measure the low-energy ions(less than a few tens of electronvolts) in the polar caps/magnetotail lobes in the years 2001–2010. With the measured velocities, we are able to study the trajectories of these low-energy ions. Particle tracing has been used in previous studies, confirming that ions of ionospheric origin are the dominant contributor to the ion population in the Earth's magnetotail lobes. In this work, we continue to study the source of low-energy ions measured in the lobes. We found that not all of the low-energy ions in the lobes come directly from the ionosphere. Particle tracing infers that some of the low-energy ions start to move tailward from the cusp/near-cusp region with a zero parallel velocity. In the following, we refer to these low-energy ions as stagnant low-energy ions. On the other hand, the in situ measurements by Cluster show a population of low-energy ions in the cusp/near-cusp region with pitch angles near 90°(i.e., no significant parallel velocity).The locations of stagnant low-energy ions are determined by particle tracing and in situ measurements. Similar ion energies and spatial distributions determined by these two methods confirm the presence of the stagnant low-energy ion population.     

2017年四川九寨沟MS7.0强震的余震重定位及主震震源机制反演

本文采用双差定位法对2017年8月8日至10月31日期间四川九寨沟M_S7.0主震及5200个余震序列进行相对定位,得到4036个重定位地震事件.采用中国区域地震台网观测到的宽频带垂直分向波形数据和W震相反演方法,得到了主震震源机制解.重定位结果显示,余震序列分别沿NNW和SSE两个方向扩展,展布长度约58 km,且这些余震主要集中在22 km深度之上.余震分布的另一个重要特点是具有分区特性,即在主震NNW方向约5 km处存在明显的西北和东南两区余震活动分界线;西北区的余震由深至浅具有较好连续性,而东南区却在约10 km深度处存在不连续性.余震分布的这种分区特征,说明九寨沟地震震源区的地壳结构存在强烈的不均匀性.余震分布与主震破裂特征的一致性,证实了我们定位结果的可靠性.主震的震源机制解展示出节面Ⅰ的走向/倾角/滑动角分别为246°/83.7°/-177°,而节面Ⅱ的走向/倾角/滑动角为155.7°/87.1°/-6.3°,最佳质心深度为15.5 km,矩震级M_W为6.5.根据余震分布较为垂直和主震震源机制解两节面的倾角均在80°以上,并结合野外地质调查结果,推测此次九寨沟地震为与节面Ⅱ参数相近的一次高角度的左旋走滑型事件.     

基于偏微分方程的模型光滑处理在偏移成像中的应用

本文针对射线类偏移成像当中的速度模型光滑处理问题,借鉴数字图像处理当中的偏微分方程法,基于能量泛函,应用变分方法导出基于速度模型的偏微分方程实现射线类偏移成像当中的速度模型的光滑处理.由于偏微分方程法具有线性叠加特性、模型解的唯一性和局部特征保持性,因此,应用该算法可以实现基于原始速度模型空间结构的模型光滑处理.通过在原始速度模型以及光滑处理后的速度模型上计算速度的空间分布以及地震波走时、射线路径可以得出,偏微分方程法对速度模型的光滑处理能够很好地保持原始模型的空间结构,偏移成像结果也证明了该方法的实用性.     

NEW EVIDENCES FOR LATE QUATERNARY ACTIVITY IN THE SOUTHERN SEGMENT OF THE YISHU-TANGTOU FAULT,THE TAN-LU FAULT ZONE,AND ITS TECTONIC IMPLICATION

The Tan-Lu Fault Zone(TLFZ), a well-known lithosphere fault zone in eastern China, is a boundary tectonic belt of the secondary block within the North China plate, and its seismic risk has always been a focus problem. Previous studies were primarily conducted on the eastern graben faults of the Yishu segment where there are historical earthquake records, but the faults in western graben have seldom been involved. So, there has been no agreement about the activity of the western graben fault from the previous studies. This paper focuses on the activity of the two buried faults in the western graben along the southern segment of Yishu through combination of shallow seismic reflection profile and composite drilling section exploration. Shallow seismic reflection profile reveals that the Tangwu-Gegou Fault(F₄)only affects the top surface of Suqian Formation, therefore, the fault may be an early Quaternary fault. The Yishui-Tangtou Fault(F₃)has displaced the upper Pleistocene series in the shallow seismic reflection profile, suggesting that the fault may be a late Pleistocene active fault. Drilling was implemented in Caiji Town and Lingcheng Town along the Yishui-Tangtou Fault(F₃)respectively, and the result shows that the latest activity time of Yishui-Tangtou Fault(F₃)is between(91.2±4.4)ka and(97.0±4.8)ka, therefore, the fault belongs to late Pleistocene active fault. Combined with the latest research on the activity of other faults along TLFZ, both faults in eastern and western graben were active during the late Pleistocene in the southern segment of the Yishu fault zone, however, only the fault in eastern graben was active in the Holocene. This phenomenon is the tectonic response to the subduction of the Pacific and Philippine Sea Plate and collision between India and Asian Plate. The two late Quaternary active faults in the Yishu segment of TLFZ are deep faults and present different forms on the surface and in near surface according to studies of deep seismic reflection profile, seismic wave function and seismic relocation. Considering the tectonic structure of the southern segment of Yishu fault zone, the relationship between deep and shallow structures, and the impact of 1668 Tancheng earthquake(M=8(1/2)), the seismogenic ability of moderate-strong earthquake along the Yishui-Tangtou Fault(F₃)can't be ignored.     

松嫩地块东缘和佳木斯地块西缘电性结构

横过松嫩地块东缘和佳木斯地块西缘的大地电磁测深剖面揭示了两块体结合带附近的深部电性结构.本文对剖面测点做了标准化数据处理,并对二维偏离度、构造走向进行了计算和分析,采用非线性共轭梯度(NLCG)二维反演方法对TM模式的数据进行了反演,获得了该剖面的地壳、上地幔电性结构模型,划分出三个典型构造单元:松嫩地块东缘、碰撞拼合带和佳木斯地块西缘.研究结果表明,研究区上地壳基本呈高阻特征,可能为岩浆岩,代表其经历了多期次岩浆作用,而松嫩地块东缘和佳木斯地块西缘的中下地壳的高导体可能与地幔物质的上涌有关;拼合带下方存在西倾的高导体和高阻体,可能是佳木斯地块向西俯冲到松嫩地块下方的构造遗迹;研究区可能发生了拆沉作用,与之伴随的地幔物质上涌可能是后期伸展作用的一个动力.     

2017年四川九寨沟MS7.0地震前区域重力场变化

利用南北地震带2014-2017年期间的流动重力观测资料,系统分析了区域重力场变化及其与2017年8月8日四川九寨沟7.0级地震发生的关系.结果表明：①区域重力场异常变化与北西西向塔藏断裂和南北向岷江断裂带在空间上关系密切,反映了沿控震断裂在2013-2017年期间发生了引起地表重力变化效应的地壳变形和构造活动.②九寨沟地震前,测区内出现了大空间范围的区域性重力异常,而震源区附近产生了局部重力异常,沿塔藏断裂带形成了重力变化高梯度带,其中,甘肃玛曲、迭部、青海河南蒙古族自治县、四川若尔盖、九寨沟一带重力差异变化达100×10^-8m·s^-2以上;这些可能反映九寨沟地震前,区域及震源区附近均产生与该地震孕育、发生有关的构造运动或应力增强作用.③九寨沟地震震中位于重力差异运动剧烈的鞍部等值线附近,与断裂走向基本一致的重力变化高梯度带零值线上,这一观测事实进一步佐证了重力场动态变化图像对强震地点预测具有重要的指示意义.     

层状半空间中周期分布凸起地形对平面SH波的散射

提出了一种新的以层状半空间中周期分布斜线荷载动力格林函数为基本解的间接边界元方法,研究了周期分布凸起地形对平面SH波的散射问题.方法将散射波场分解为凸起内部散射波场和凸起外部散射波场.凸起内部散射波场通过在凸起闭合边界上施加虚拟斜线荷载产生的动力响应来模拟,而凸起外部散射波场则通过在凸起与半空间交界面上施加虚拟周期分布斜线荷载产生的动力响应来模拟.周期分布斜线荷载动力格林函数的引入,使得本文方法仅需针对一个凸起进行边界单元的离散和求解,便可完成问题的求解,避免了通过截断无限边界求解而引入的误差,方法具有较高精度的同时显著降低了求解自由度.文中通过与已有结果的比较,验证了方法的正确性,并以均匀半空间和基岩上单一土层中周期分布凸起为例进行了数值计算分析.研究表明,凸起间距对凸起地形间的动力相互作用有着显著的影响,同时层状半空间中周期分布凸起地形对SH波的散射与均匀半空间情况也有着显著的差别.