Characteristics of horizontal ground motion measures along principal directions

<正>Ground motion records are often used to develop ground motion prediction equations(GMPEs) for a randomly oriented horizontal component,and to assess the principal directions of ground motions based on the Arias intensity tensor or the orientation of the major response axis.The former is needed for seismic hazard assessment,whereas the latter can be important for assessing structural responses under multi-directional excitations.However,a comprehensive investigation of the pseudo-spectral acceleration(PSA) and of GMPEs conditioned on different axes is currently lacking.This study investigates the principal directions of strong ground motions and their relation to the orientation of the major response axis, statistics of the PSA along the principal directions on the horizontal plane,and correlation of the PSA along the principal directions on the horizontal plane.For these,three sets of strong ground motion records,including intraplate California earthquakes,inslab Mexican earthquakes,and interface Mexican earthquakes,are used.The results indicate that one of the principal directions could be considered as quasi-vertical.By focusing on seismic excitations on the horizontal plane,the statistics of the angles between the major response axis and the major principal axis are obtained;GMPEs along the principal axes are provided and compared with those obtained for a randomly oriented horizontal component;and statistical analysis of residuals associated with GMPEs along the principal directions is carried out.     

Blasius方程中临界值的分析估计

证明了流体力学边界层理论中著名的Blasius方程的一个新性质,它在二维平板流的研究中具有重要的作用,由此建立它的积分形式,并给出方程中的临界值βc的分析估计。     

考虑软土结构性的折线型强度公式及应用

结合上海某植物园堆山工程,通过直剪试验及三轴剪切试验,得到考虑软土结构性的折线型强度公式和不考虑软土结构性的直线型强度公式的参数,采用这两种方法对山体荷载作用下不同位置处的软土地基的强度增长进行了预测,并与现场十字板剪切试验结果进行了对比。结果表明：由于该工程的软土地基具有明显的结构性,采用直线型强度公式得到的地基强度增长量要大于采用采用折线型的增长量;另外,在应用考虑土的结构性的折线型强度公式时,应该尽量采用三轴试验参数而不是直剪试验的参数,这样得到的结果更接近十字板试验结果。     

基于逐减随机震源采样法的频率域二维黏滞声波方程全波形反演

全波形反演方法利用叠前地震波场的运动学和动力学信息重建地下速度结构,具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而,巨大的计算量是阻碍其发展的一个瓶颈问题。为此,研究者们提出了震源编码技术来减少计算量,但是此方法在模型更新过程中会引进随机串扰噪声,降低反演结果准确性。所以,在保证计算精度的情况下,本文提出了采用逐减随机震源采样的方法来高效计算全波形反演问题。笔者将此方法应用于频率域二维黏滞声波波动方程全波形反演,开始了在频率域进行随机震源采样类方法的研究,计算过程中共使用了依次增大的8个频率段;并应用Overthrust模型来验证此类随机震源采样法的正确性。实验结果表明：基于逐减随机震源采样法的反演结果与实际Overthrust模型的拟合误差为0.065 65,而应用基于全部震源的全波形反演方法得到的反演结果与实际Overthrust模型的拟合误差为0.064 64,两者差别不大;但计算用时由740 min减少到291.2 min,即计算效率提高了2.54倍。为了更好地确定方法的有效性,将其应用于Marmousi模型进行试算。模型试算结果表明：基于逐减随机震源和基于全部震源得到的反演结果与实际Marmousi模型的拟合误差分别为0.080 12和0.078 97,相差不大;但计算用时由1 218.9 min减少到274.4 min,计算效率提高了4.44倍。综上,在保证反演精度的情况下,基于逐减随机震源采样法的频率域全波形反演方法大大减少了计算量,具有不可替代的计算优势,并且没有引进随机串扰噪声。     

暴雨内涝危险性情景模拟方法研究——以上海中心城区为例

本文采用地理信息系统技术（GIS）与水文/水动力学相结合的方法,使用修正的SCS模型进行产流模拟,利用局部等体积法和水动力模型进行汇流模拟,建立了基于情景的城市暴雨内涝危险性模拟工具,并对相同雨强情景下不同汇流模型的模拟结果与精度进行了对比分析。结果表明：（1）局部等体积法计算简单,对降雨的时程分配较为敏感,峰前历时和雨强对结果影响显著,较适合于雨型确定的城市暴雨内涝危险性快速模拟;水动力模型计算复杂,与整个降雨历时存在明显关系,模拟精度较高。（2）利用台风“麦莎”带来强降雨导致的积水实测数据,对模拟结果进行对比和验证,表明本文建立的水动力模型法模拟结果与实测降雨积水更为接近,更适合上海中心城区暴雨内涝危险性情景模拟。     

Identification of diffusion parameters in a nonlinear convection–diffusion equation using the augmented Lagrangian method

Numerical identification of diffusion parameters in a nonlinear convection–diffusion equation is studied. This partial differential equation arises as the saturation equation in the fractional flow formulation of the two-phase porous media flow equations. The forward problem is discretized with the finite difference method, and the identification problem is formulated as a constrained minimization problem. We utilize the augmented Lagrangian method and transform the minimization problem into a coupled system of nonlinear algebraic equations, which is solved efficiently with the nonlinear conjugate gradient method. Numerical experiments are presented and discussed. This work was partially supported by the Research Council of Norway (NFR), under grant 128224/431.     

On the scavenging of SO2 by cloud and raindrops: II. An experimental study of SO2 absorption and desorption for water drops in air

For the purpose of testing our previously described theory of SO₂ scavenging a laboratory investigation was carried out in the UCLA 33 m long rainshaft. Drops with radii between 250 and 2500 m were allowed to come to terminal velocity, after which they passed through a chamber of variable length filled with various SO₂ concentrations in air. After falling through a gas separating chamber consisting of a fluorocarbon gas the drops were collected and analyzed for their total S content in order to determine the rate of SO ₂ absorption.The SO₂ concentration in air studied ranged between 1 and 60% (v). Such relatively large concentrations were necessary due to the short times the drops were exposed to SO₂ in the present setup. The present experimental results were therefore not used to simulate atmospheric conditions but rather to test our previously derived theory which is applicable to any laboratory or atmospheric condition. Comparison of our studies with the results from our theory applied to our laboratory conditions led to predicted values for the S concentration in the drops which agreed well with those observed if the drops had radii smaller than 500 m. In order to obtain agreement between predicted and observed S concentrations in larger drops, an empirically derived eddy diffusivity for SO₂ in water had to be included in the theory to take into account the effect of turbulent mixing inside such large drops.In a subsequent set of experiments, drops initially saturated with S (IV) were allowed to fall through S-free air to determine the rate of SO ₂ desorption. The results of these studies also agreed well with the results of our theoretical model, thus justifying the reversibility assumption made in our theoretical models.In a final set of experiments, the effects of oxidation on SO₂ absorption was studied by means of drops containing various amounts of H₂O₂. For comparable exposure times to SO₂, the S concentration in drops with H₂O₂ was found to be up to 10 times higher than the concentration in drops in which no oxidation occurred.     

缓变地形下β效应的Rossby代数孤立波

正压流体中,从准地转位涡方程出发采用摄动方法和时空伸长变换推导了在缓变地形下β效应的Rossby代数孤立波方程,得到Rossby波振幅满足带有缓变地形非齐次Benjamin-Davis-Ono(BDO)方程的结论.通过分析孤立Rossby波振幅的演变,指出了β效应、地形效应是诱导Rossby孤立波产生的重要因素,说明了在缓变地形强迫效应和非线性作用相平衡的假定下,Rossby孤立波振幅的演变满足非齐次BDO方程,给出在切变基本气流下缓变地形和正压流体中Rossby波的相互作用.     

Perfectly matched layer-absorbing boundary condition for finite-element time-domain modeling of elastic wave equations

The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second-order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite-element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.