Horizontal and vertical components of earthquake ground motions at liquefiable sites

Field observations on ground motions from recent earthquakes imply that current knowledge is limited with regard to relating vertical and horizontal motions at liquefiable sites. This paper describes a study with the purpose of clarifying this emerging issue to some extent. A series of numerical analyses is carried out on a liquefiable soil deposit with a verified, fully coupled, nonlinear procedure. It is shown that the transformation of vertical motions in the deposit differs considerably from the transformation of horizontal motions. Both the amplitude and frequency content of the horizontal motions are strongly dependent on the shaking level or the associated nonlinear soil behavior. The transfer function for vertical motions is however likely to be independent of the intensity of input motions; no reduction in the amplitude occurs even in the case of strong shaking. The results are shown to be in consistence with the laboratory observations on shaking table tests and recent field observations that less nonlinearity exists for vertical motions. It is also shown that the possibility exists for using information on spectral ratios between the horizontal and vertical surface motions to quickly identify in situ soil behavior and liquefaction that are not readily covered by conventional field or laboratory experimentation procedures.     

Computational modeling of cyclic mobility and post-liquefaction site response

Under seismic excitation, liquefied clean medium to dense cohesionless soils may regain a high level of shear resistance at large shear strain excursions. This pattern of response, known as a form of cyclic mobility, has been documented by a large body of laboratory sample tests and centrifuge experiments. A plasticity-based constitutive model is developed with emphasis on simulating the cyclic mobility response mechanism and associated pattern of shear strain accumulation. This constitutive model is incorporated into a two-phase (solid–fluid), fully coupled finite element code. Calibration of the constitutive model is described, based on a unique set of laboratory triaxial tests (monotonic and cyclic) and dynamic centrifuge experiments. In this experimental series, Nevada sand at a relative density of about 40% is employed. The calibration effort focused on reproducing the salient characteristics of dynamic site response as dictated by the cyclic mobility mechanism. Finally, using the calibrated model, a numerical simulation is conducted to highlight the effect of excitation frequency content on post-liquefaction ground deformations.     

PRELIMINARY NUMERICAL STUDY OF TOPOGRAPHIC EFFECTS OF THE TIBETAN PLATEAU ON SURFACE DIRECT RADIATION*

To improve the simulating ability of a model,this paper presents a scheme of calculating direct radiation at land surface with topography in the model.A numerical study is conducted for the topographic effects of the Tibetan Plateau on the direct radiation using NCEP terrain data.Results show that,after taking account into the topographic radiation effect,the regional average of the radiation over the Plateau obviously increases in the local early morning and late afternoon,but changes less around noon.The effect is stronger in winter than that in summer.And heterogeneous topography has also affected the distribution of the radiation in this area.A simple numerical experiment shows that considering the effect will lead ground temperature to increase on the slope having more sunshine,and vice versa.     

DSQ型水管倾斜仪电压灵敏度的调整技巧

DSQ型水管倾斜仪是“九五”期间推出的数字化地震前兆仪器，本文介绍了一些如何调整换能器的电压灵敏度及调整中的技巧。     

由GPS观测结果推导中国大陆现今水平应变场

以中国大陆及周边近400个GPS测站的水平运动速率为基础，给出了现今地壳水平应变场结果表明:①中国大陆水平应变为西强东弱，剪应变数值大于正应变数值(绝对值)，应变量级一般为10-8/a，局部区域达到10-7/a，但应变分布不均匀；②南北向应变最突出的部位为中国西南部西段的喜马拉雅条带、西北部的36N~42N段及柴达木断块的北缘；③东西向应变西边缘变化最大．此外，由西向东还具有正负交替的变化特征；④REN(东-北向剪切应变)与Rmax(最大剪切应变)数值较大的区域分别是喜马拉雅条带、西北部的36N～42N段、柴达木断块的西部、川滇菱形块体，以及阿拉善、祁连及塔里木断块的交界区；⑤青藏块体周边以面收缩为主，内部则以面膨胀为主．其以北的地区以面收缩为主．西界数值最大，东部数值最小(除燕山构造带外)；⑥西部西区主压应变为南北向，主张应变为东西向．西部东缘区主压应变为近东西向，主张应变为近南北向．川滇菱形块体主应变的方向发生了很大的变化，北部地区为东西压南北张，南部地区则恰好相反；⑦中国大陆的应变模式可能是断块模式与连续模式的组合．此外，小尺度优势应变可能是剪切应变．造成上述结果除与印度板块的碰撞及边界耦合有关外，还与深部物质的活动及地壳介质的物性有密切的关系．必须指出，由于GPS测站在空间上分布的不均匀性，那么，由此而来的应变场，其应变尺度也不一样．      

新疆伽师强震群区的横波分裂与应力场特征

利用在新疆伽师地区布设的流动台阵记录到的地方震波形数据,研究了伽师强震群附近各台站横波分裂现象,给出了相应的快波偏振方向的平面场分布. 发现在台阵的塔里木盆地一侧,波偏振方向为近SN向,与塔里木盆地的区域主压应力方向一致,但在塔里木盆地北部边缘的褶皱变形带内,快波偏振方向变为近EW向,特别是在柯坪断裂附近,快波偏振方向与阿图什地震的震源断层方向基本一致. 由于快波偏振方向平行于主压应力方向,给出的快波偏振方向反映了相应的主压应力场特征. 结果表明,伽师强震群的成因很可能是塔里木盆地北缘横向非均匀变形造成的局部张性剪切应变能的释放.     

南向行星际磁场事件与磁暴关系的研究

利用172-182年IMP-8飞船的太阳风观测资料和相应地磁活动性指数Dst和AE,研究了43个南向行星际磁场事件期间太阳风和磁层的耦合问题. 与这43个事件对应的地磁暴是中等的和强的磁暴(Dst＜-50nT). 结果表明:(1) 在43个事件中有11个(约占25.6髎)紧随激波之后,18个处于激波下游流场中(占42髎),其余14个(占33髎)和激波没有关连. 绝大多数事件都伴有太阳风动压和总磁场强度的增加;(2) 当行星际晨昏向电场强度EI＞-4mV/m时,只引起磁亚暴,对Dst指数没有明显影响. 仅当EI＜-5mV/m时,磁亚暴和磁暴才会同时出现;(3) 太阳风动压的增加会增强能量向环电流的输入,但不是密度和速度单独起作用,而是以PK=ρV2的组合形式影响能量的输入;(4) 虽然行星际磁场(IMF)南向分量BZ对太阳风和磁层的耦合起着关键作用,但IMF的BX和BY分量相对于BZ的大小对太阳风向磁层的能量传输也有一定影响. 当BX、BY相对BZ较大时能量耦合加强.     

中朝地台东北缘地区的地震层析成像

根据中朝地台东北缘地区 (东经 1 1 7°0 0′— 1 2 6°0 0′ ,北纬 36°0 0′— 44°0 0′) 1 980— 1 997年的 380 0 0余条P波走时数据 ,利用正交投影法重建了该区地壳和上地幔的三维速度结构 .通过分析及同人工地震测深剖面的详细对比 ,证明了成像结果的可靠性 .结果表明 :中朝地台东北缘地区地壳上地幔介质存在显著的横向不均匀性 ,直至 1 2 0km深度处依然明显 ;地壳上部的速度图像清楚地反映了不同岩石单元的分布与该区不同性质的基岩分布基本吻合 ;从上、中地壳的速度图像中发现了研究区存在海城、朝阳、义县、丹东南、唐山等几个低速异常区(即速度逆反层区 ) ,其中海城、唐山、朝阳等地区的壳内低速层已由深地震测深资料所证实 ;研究区陆地发生的几次强震均发生于壳内低速层上方的高速脆性介质内 ,而渤海发生的强震 ,此现象不明显 ,但都发生于横向介质速度显著突变的位置 ;在地壳不同深度上发现了普兰店至山海关横跨渤海的北西向低速异常带 ;地震层析二维速度图像与深地震测深资料的对比表明 ,研究区利用地震层析成像技术 ,在一定条件下可以获得与人工地震测深相似的效果     

剥蚀及地幔作用下青藏高原隆升过程的数值模拟

修改了England和Mckenzie的黏性薄层流变模型中控制大陆形变的连续性方程,将剥蚀作用对高原隆升演化的影响直接引入该方程,并考虑下伏地幔小尺度对流对增厚岩石层的搬离作用对高原隆升演化后期的影响,用有限差分法直接模拟青藏高原隆升过程. 数值模拟结果所显示的高原隆升演化过程与实际观测资料吻合较好,揭示了高原隆升演化过程的非平稳和多阶段的特性;同时还表明上地幔小尺度对流对岩石层底部的搬离作用可能是最近8Ma以来高原快速隆升的主导机制.     

三维地震数据离散光滑插值的共轭梯度法

针对三维地震数据插值,提出采用Laplacian算子进行光滑约束的插值方法,并借鉴Mallet研究的离散光滑插值思路,采用预条件共轭梯度法,直接生成网格节点上的值,从而回避寻求满足插值方程的函数. 为了实现其中Laplacian算子的快速求逆,文中引入Claerbout螺旋坐标系谱因式分解理论. 在螺旋坐标系下,Laplacian算子的表示矩阵具有Toeplitz结构,其快速求逆可由谱法LU分解实现. 基于二维离散光滑插值,文中还给出共轭梯度法与NMO相结合的沿时间切片逐层处理的离散光滑插值流程. 最后,应用该方法对模型数据和实际三维地震数据进行了处理.