模拟地震波传播的三维逐元并行谱元法

高效地震波场正演模拟对于复杂模型中地震波传播与成像研究至关重要.本文在谱元法原理框架内,对已有逐元谱元法改进,提出一种新的逐元并行谱元法求解三维地震波运动方程,并得到地震波场.逐元并行谱元法的核心思想在于在单元上进行质量矩阵与解向量的乘积运算,并将此运算平均分配至每一个CPU计算核心,此处理有利提升谱元法的并行计算效率...     

基于分数阶时间导数的黏弹性衰减VTI介质中平面波传播

目前在地震勘探频带范围内通常假设品质因子Q与频率无关,且呈衰减各向同性.事实上,相比较速度各向异性,介质的衰减各向异性同样不可忽视.本文将衰减各向异性和速度各向异性二者与常Q模型相结合,建立了黏弹性衰减VTI介质模型,并基于分数阶时间导数理论,给出了对应的本构关系和波动方程.利用均匀平面波分析和Poynting定理,推...     

FY-4A数据在2021年1月6—8日寒潮监测中应用

基于静止气象卫星图像解译理论和位涡理论,应用FY-4A水汽图像、干涉式大气垂直探测仪（GIIRS）产品和欧洲中心第五代大气再分析资料（ERA5）对2021年1月6—8日我国中东部地区寒潮天气进行分析。结果表明：高空冷涡和地面冷高压是影响此次寒潮天气的主要影响系统,冷空气5日08:00从贝加尔湖南部向东、向南加强发展,8日20:00 减弱;对比FY-4A/GIIRS与ERA5 850 hPa 24 h变温数据发现：两种数据反映的对流层低层冷空气东移南下特征和正负变温区域分布特征基本一致,FY-4A/GIIRS 850 hPa 24 h负变温中心和0℃等变温线的位置与ERA5数据接近,FY-4A/GIIRS 850 hPa 24 h变温中心值强度略大于ERA5。水汽图像暗区变暗附近的高层位涡的增加激发高层气旋性环流后侧下沉运动增强,干空气下沉造成了地面高压增强;500 hPa冷涡中心附近绝对涡度增加是引起冷空气增强的原因之一;冷空气急剧增强时段,中低层位涡的加大引起地面冷空气堆积,寒潮增强。FY-4A/GIIRS温度产品反映的冷空气移动特征与位涡理论对冷空气增强原因分析结果一致,说明应用FY-4A/GIIRS温度产品能够有效监测分析寒潮冷空气演变特征。     

2023年全球重大天气气候事件

2023年,全球平均气温比工业化前高出约1.45℃(±0.12℃),是有观测记录以来最热的一年。全球海平面继续上升,且全球平均海平面达到了有卫星记录(1993年至今)以来的最高水平,反映了持续的海洋变暖以及冰川和冰盖的融化。北极海冰面积仍远低于常年值,南极海冰面积创下历史新低。巴基斯坦、中国京津冀地区、意大利、巴西圣保罗州北部沿海地区、新西兰北岛等地遭受暴雨洪涝灾害,非洲西北部、中国云南、中美洲和南美洲北部发生严重干旱,南欧、北美、南美、东亚和南亚等地遭遇创纪录高温热浪,欧洲和北美等地遭遇寒流和暴风雪侵袭,强对流天气频繁袭击世界各处,全球热带气旋活动频繁。     

Shear wave propagation along infinite slopes: A theoretically based numerical study

In this paper, the seismic response of ‘infinitely’ long slopes is numerically analysed via the formulation of a 1D analytical/numerical model, in which the soil mechanical behaviour is assumed to be elasto‐perfectly viscoplastic and simple shear (SS) kinematical constraints are imposed. In order to simplify the problem, a theoretically based procedure to set up a fully 1D shear constitutive model is defined, within which the mechanical response of a multiaxial relationship is condensed. The use of a 1D shear constitutive model is aimed at reducing the number of unknowns and, therefore, the computational costs. In particular, the case of the Mohr–Coulomb yield criterion is considered, while an enhanced Taylor–Galerkin finite element algorithm is employed to simulate the seismic wave propagation within the soil stratum. The proposed ‘condensation’/calibration procedure captures both the ‘pseudo’‐hardening pre‐failure behaviour and the influence of dilation on the occurrence of strain‐localization, which characterize, under SS conditions, the static response of virgin perfectly plastic soils. The effectiveness of the conceived method is shown with reference to freshly deposited deposits, while, in the case of highly overconsolidated strata, some difficulties arise because of the brittle behaviour induced both by unloading and non‐associativeness. Copyright © 2011 John Wiley & Sons, Ltd.     

Falling sphere observations of anisotropic gravity wave motions in the upper stratosphere over Australia

A study of inertial scale gravity wave motions in the region of the atmosphere between 30 and 60 km has been undertaken, using wind and temperature data derived from rocket-borne falling sphere density experiments performed over Woomera, Australia between 1962 nad 1976. The gross features of the wave field compare favorably with those found in similar northern hemispheric studies. Wave propagation is found to be both vertically and horizontally anisotropic. A rotary spectral analysis indicates predominately upgoing wave energy, suggesting that the majority of sources of these waves lie below 30 km. A detailed statistical investigation of the waves, made using the Stokes parameters technique, reveals that phase progression is also highly directional in the horizontal, with a significant zonal component in summer, but with a strong meridional component in winter. Propagation towards the southeast is inferred in summer, with the waves possibly emanating from tropospheric sources in equatorial regions to the north of Australia. The technique also shows that, on average, the waves appear to have mean ellipse eccentricities (=f/) around 0.4–0.45. Indirect estimates of a number of important wave parameters are made. In particular,v andw flux estimates are made over several height intervals. The vertical gradient of density weighted flux implies wave-induced mean flow accelerations of the order 0.1–1 ms^–1day^–1. This suggests that dissipating gravity waves are a significant source of the momentum residuals that are encountered in studies of satellite data from this region.     

透射波波场重建和介质层厚反演研究

本文以P波和SV波波场比值,提出层状介质情况下透射波波场的反向重建算法原理,并据此给出在已知介质速度值时,求介质层厚度的反演方法。文中还以深源远震记录波形资料给出反演算例,证实本文的算法是可行的。     

应力波在波阻抗连续变化的弹性杆中的演化

利用特征线方法数值模拟计算了应力波在截面或弹性模量连续变化的弹性杆中的反射,同时也讨论了应力波脉冲宽度对反射波波形特征的影响,其结果对桩基测量有一定的参考价值。     

Oscillatory Rossby Solitary Waves in the Atmosphere

The linear Rossby wave frequency expression is expanded at higher accuracy based on the scale difference char?acteristics of atmospheric long waves in the and directions. That the nature of the waves represented by the expan?sion is identical to that of the original ones is demonstrated both in phase velocity and wave energy dispersion speed , followed by the derivation of the nonlinear expression describing atmospheric long wave behaviors with the associated approximate analytic solution obtained. Then, for the first time atmospheric’ oscillatory Rossby solitary wave’ with its dispersion relation is obtained by numerical calculation with the aid of physical parameters of the real atmosphere. The solitary wave is found to be very close to such longwave systems as blocking highs and cut-off de?pressions in the actual atmosphere.