黏弹性介质中瑞雷波有限差分数值模拟与波场分析

瑞雷面波经常被用来反演地表浅层横波速度,受到越来越广泛的关注。对瑞雷波的研究一般都基于完全弹性介质,而实际地层更接近黏弹性介质,对黏弹性介质中的瑞雷面波进行模拟更具实际意义。本文采用广义标准线性体模型来描述黏弹性介质,并采用交错网格有限差分法对考虑水平自由表面的黏弹介质进行正演模拟,再与弹性介质中的结果进行对比分析。首先采用非线性最优化算法根据期望常数品质因子直接求取松弛时间来拟合常Q模型,并给出广义标准线性固体的具体算例,实施自由表面条件时采用声学-弹性边界近似法,通过剪切模量不变来考虑自由表面上、下横向应力保持连续的条件。对于非自由表面,采用非分裂的多轴卷积完全匹配层来吸收波场。然后对几种典型的数值模型进行正演模拟计算,数值解与解析解的对比验证了本文方法的准确性与有效性,正演结果的对比表明波场尤其是面波频散会受黏弹性影响,因此有必要在面波勘探中考虑黏弹性因素。     

某地区黑云角闪岩应力松弛特性研究

岩体工程失稳往往由工程岩体蠕变、应力松弛等时效力学行为下的损伤演化所致。本文以某坝基黑云角闪岩为研究对象,进行了2组单轴应力松弛试验,获取了不同应力条件下的应力松弛及其相应的应力松弛速率曲线。试验结果表明：根据应力松弛曲线形态可分为快速松弛,减速松弛及应力稳定3个阶段,应力松弛的演化与岩石裂隙的萌生和流变变形密切相关。在此基础上,提出了采用非线性Maxwell模型描述岩石应力松弛特性,模拟结果与试验结果吻合较好。本文研究可为类似工程的长期安全稳定性研究提供理论依据。     

图像重建块迭代算法中松弛参数选取的研究

在图像重建中,Landweber迭代算法是图像重建算法中的重要方法。本文将针对Landweber分块迭代算法中松弛参数的选取进行研究。在重建过程中采取对投影矩阵按投影角度分块的方法,选取特定的松弛参数。通过数值实验得出结论:对于按角度分块的块迭代算法,松弛参数选取为λ乘以块矩阵与其共轭转置矩阵乘积的最大特征值分之一,当采集完全投影数据,且λ接近(1/6)~(1/7)时效果最好。另外,本文按角度分块的做法和松弛参数的选取方法对于有限角度图像重建问题也是可行的。     

Design of controlled elastic and inelastic structures

One of the founders of structural control theory and its application in civil engineering,Professor Emeritus Tsu T.Soong,envisioned the development of the integral design of structures protected by active control devices.Most of his disciples and colleagues continuously attempted to develop procedures to achieve such integral control.In his recent papers published jointly with some of the authors of this paper,Professor Soong developed design procedures for the entire structure using a design-redesign proce...     

地面核磁共振找水方法在花岗岩地区的应用

文章简要地阐述了地面NMR方法的直接找水原理、特点以及仪器组成、测量质量的控制.通过在花岗岩地区的找水实践,说明了该方法的可行性和优点,最后总结了影响NMR信号采集和NMR信号解释结果的因素及解决办法,并提出了几点建议及今后的研究方向.     

基于松弛法的影像边缘提取

提出一种基于松弛法的影像边缘提取算法.本文详述了边缘标号的设置及初始概率的计算,并根据邻域内的信息分布,综合考虑了边缘的先验知识,设计出自适应的支持函数.     

Using NWP to Assess the Influence of the Arctic Atmosphere on Midlatitude Weather and Climate

The influence of the Arctic atmosphere on Northern Hemisphere midlatitude tropospheric weather and climate is explored by comparing the skill of two sets of 14-day weather forecast experiments using the ECMWF model with and without relaxation of the Arctic atmosphere towards ERA-Interim reanalysis data during the integration. Two pathways are identified along which the Arctic influences midlatitude weather: a pronounced one over Asia and Eastern Europe, and a secondary one over North America. In general, linkages are found to be strongest(weakest) during boreal winter(summer) when the amplitude of stationary planetary waves over the Northern Hemisphere is strongest(weakest). No discernible Arctic impact is found over the North Atlantic and North Pacific region, which is consistent with predominantly southwesterly flow. An analysis of the flow-dependence of the linkages shows that anomalous northerly flow conditions increase the Arctic influence on midlatitude weather over the continents. Specifically, an anomalous northerly flow from the Kara Sea towards West Asia leads to cold surface temperature anomalies not only over West Asia but also over Eastern and Central Europe. Finally, the results of this study are discussed in the light of potential midlatitude benefits of improved Arctic prediction capabilities.     

Multi-scale Incremental Analysis Update Scheme and Its Application to Typhoon Mangkhut (2018) Prediction

In the traditional incremental analysis update(IAU) process, all analysis increments are treated as constant forcing in a model’s prognostic equations over a certain time window. This approach effectively reduces high-frequency oscillations introduced by data assimilation. However, as different scales of increments have unique evolutionary speeds and life histories in a numerical model, the traditional IAU scheme cannot fully meet the requirements of short-term forecasting for the damping of hig...     

Propagation of seismic waves in a viscoelastic medium and the effect of topography

In this paper we have simulated the propagation of seismic waves in viscoelastic medium and calculated the effect of irregular surface topography by using finite element method. Several types of elastic and viscoelastic medium models, such as block structure, vidges and/or valley with and without a soft layer underneath have been studied. The distribution of maximum amplitude of displacement and acceleration on surface has been calculated in the case of vertically incident SH and P waves from the basement. The possible application in earthquake engineering is also discussed. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 442–449, 1991.     

An efficient finite–discrete element method for quasi‐static nonlinear soil–structure interaction problems

An efficient finite–discrete element method applicable for the analysis of quasi‐static nonlinear soil–structure interaction problems involving large deformations in three‐dimensional space was presented in this paper. The present method differs from previous approaches in that the use of very fine mesh and small time steps was not needed to stabilize the calculation. The domain involving the large displacement was modeled using discrete elements, whereas the rest of the domain was modeled using finite elements. Forces acting on the discrete and finite elements were related by introducing interface elements at the boundary of the two domains. To improve the stability of the developed method, we used explicit time integration with different damping schemes applied to each domain to relax the system and to reach stability condition. With appropriate damping schemes, a relatively coarse finite element mesh can be used, resulting in significant savings in the computation time. The proposed algorithm was validated using three different benchmark problems, and the numerical results were compared with existing analytical and numerical solutions. The algorithm performance in solving practical soil–structure interaction problems was also investigated by simulating a large‐scale soft ground tunneling problem involving soil loss near an existing lining. Copyright © 2011 John Wiley & Sons, Ltd.