Diffraction of SH-waves by topographic features in a layered transversely isotropic half-space

The scattering of plane SH-waves by topographic features in a layered transversely isotropic (TI) half-space is investigated by using an indirect boundary element method (IBEM). Firstly, the anti-plane dynamic stiffness matrix of the layered TI half-space is established and the free fields are solved by using the direct stiffness method. Then, Green’s functions are derived for uniformly distributed loads acting on an inclined line in a layered TI half-space and the scattered fields are constructed with the deduced Green’s functions. Finally, the free fields are added to the scattered ones to obtain the global dynamic responses. The method is verified by comparing results with the published isotropic ones. Both the steady-state and transient dynamic responses are evaluated and discussed. Numerical results in the frequency domain show that surface motions for the TI media can be significantly different from those for the isotropic case, which are strongly dependent on the anisotropy property, incident angle and incident frequency. Results in the time domain show that the material anisotropy has important effects on the maximum duration and maximum amplitudes of the time histories.     

汶川地震震前与同震断层滑移的比较分析

本文利用GPS观测的1999-2007年汶川震前3期地表变形数据和2008年汶川同震地表变形数据,结合地震位错理论,通过高斯变换和坐标旋转建立断层模型,运用遗传算法,反演了龙门山断裂带断层震前3期和同震滑动参数。结果表明龙门山断层震前3期平均走滑位移为-5.39mm,倾向位移为2.66mm,与同震断层滑移相比较,发现震前断层的滑移趋势与同震断层滑移一致,均为逆冲兼右旋的挤压运动。比较震前3期逆冲方向的滑移量,发现逆冲滑移有加速的现象。并根据震前和同震的断层滑动量估算了汶川地震复发周期。     

Retrieval of middle and upper atmospheric wind based on non-full circular fringe recorded by Fabry-Perot Inteferometer

Fabry-Perot Interferometer(FPI) has been used widely for wind measurements of the middle and upper atmosphere.To date, most of FPIs have been based on full-closed circular fringe, which needs 15–25 min to obtain a group of wind velocity(zonal and meridional). However, it is hard to improve the temporal resolution because full-closed circular fringe in several directions cannot be easily imaged onto the same Charge-Coupled Device(CCD) with enough airglow intensity. In this paper, a data processing method is proposed for non-full circular fringe of FPI, which can support CCD with enough area of observations in several directions simultaneously. The method is focused on the center determination of non-full fringe. It includes radial cross-section, peak coordinate determination, and center calculation. Based on the calculated center, the fringe is annular summed. Then its radius is determined subsequently using Gaussian fitting. Finally, the wind is retrieved from the fringe radius. For validation, fringes from two ground-based FPIs were used, which are deployed in Kelan(38.71°N, 111.58°E) and Xinglong(40.40°N, 117.59°E) in China. The results retrieved from non-full fringes of FPIs were compared with that from full-closed circular fringe. The averaged wind deviation between them demonstrates reasonable difference with 5.38 ms~-(1) for 892.0 nm airglow emission, 5.81 ms~-(1) for 630.0 nm emission, and 3.03 ms~-(1) for 557.7 nm emission. Besides, wind results of Xinglong FPI are compared roughly with measurements of meteor radar which is deployed in Ming Tombs of Beijing(40.3°N,116.2°E). Good agreement demonstrates that this method is robust enough for FPI wind retrieval of mesosphere and thermosphere.     

Evaluation of Hydraulic Properties of Aquitards Using Earthquake‐Triggered Groundwater Variation

The hydraulic properties of aquitards are not easily obtained because monitoring wells are usually installed in aquifers for groundwater resources management. Earthquake‐induced crust stress (strain) triggers groundwater level variations over a short period of time in a large area. These groundwater anomalies can be used to investigate aquifer systems. This study uses a poroelastic model to fit the postseismic variations of groundwater level triggered by the Chi‐Chi earthquake to evaluate the hydraulic properties of aquitards in the Jhoushuei River alluvial fan (JRAF), Taiwan. Six of the adopted eight wells with depths of 70 to 130 m showed good agreement with the recovery theory. The mean hydraulic conductivities (K) of the aquifers for the eight wells are 1.62 × 10^?4 to 9.06 × 10^?4 m/s, and the thicknesses are 18.8 to 46.1 m. The thicknesses of the aquitards are 11.3 to 42.0 m. Under the isotropic assumption for K, the estimated values of K for the aquitards are 3.0 × 10^?8 to 2.1 × 10^?6 m/s, corresponding to a silty medium. The results match the values obtained for the geological material of the drilling core and those reported in previous studies. The estimated values were combined with those given in previous studies to determine the distribution of K in the first two aquitards in the JRAF. The distribution patterns of the aquitards reflect the sedimentary environments and fit the geological material. The proposed technique can be used to evaluate the K value of aquitards using inverse methods. The inversion results can be used in hydrogeological analyses, contaminant modeling, and subsidence evaluation.     

Origin and Evolution of Aquitard Porewater in the Western Coastal Plain of Bohai Bay,China

High‐salinity paleowater from low‐permeability aquitards in coastal areas can be a major threat to groundwater resources; however, such water has rarely been studied. The chemical and isotopic compositions of porewater extracted from a 200‐m‐thick Quaternary sedimentary sequence in the western coastal plain of Bohai Bay, China, were analyzed to investigate the salinity origin and chemical evolution of porewater in aquitards. Porewater samples derived at depths shallower than 32 m are characterized by Cl‐Na type saline water (total dissolved solids [TDS], 10.9–84.3 g/L), whereas those at depths greater than 32 m comprise Cl·SO₄‐Na type brackish water (TDS, 2.2–6.3 g/L). Saline porewater is interpreted as evaporated seawater prior to halite saturation, as evidenced by Cl‐Br relationships. Although substantial dilution of saline porewater with meteoric water is supported by a wider Cl^? range and δ²H‐δ¹⁸O covariance, the original marine waters were not completely flushed out. The deeper brackish porewater is determined to be a mixture of fresher porewater and brine groundwater and had a component of old brine of less than 10%, as indicated by a mixing model defined using δ²H and Cl^? tracers. Porewater δ²H‐δ¹⁸O relationships and negative deuterium excess ranging from ?25.9‰ to ?2.9‰ indicate the existence of an arid climate since Late Pleistocene in Tianjin Plain. The aquitard porewaters were chemically modified through water‐rock interactions due to the long residence time.     

天津地区土动力学参数变异性对地表地震动参数的影响

收集天津地区十多年来的83份地震安全性评价报告,统计202个钻孔的1 650组动三轴数据,分别给出不同种类土在不同区间深度下的统计代表值及其标准差。以某典型Ⅲ类场地为例,用等效线性化方法进行多种地震动强度及相位输入下的水平成层场地地震反应分析计算,详细研究该地区覆盖土层动剪模量比和阻尼比变异性对地表峰值加速度及其反应谱的影响。结果表明:场地地表峰值加速度和反应谱随土动剪切模量比增大或减小而增大或减小,随土动阻尼比增大或减小而减小或增大;在大震输入条件下,地表峰值加速度和地表反应谱的中、高频段随土动剪切模量比减小而减小的变化尤为明显,动阻尼比变化也有一定影响,但不如土动剪切模量比减小变化时影响明显;在中震、小震地震动输入条件下,场地土的动剪模量比和阻尼比变异性对地表峰值加速度和反应谱影响并不显著。     

沉积物粒度组分空间预测的神经网络残余kriging方法

针对近海表层沉积物粒度组分空间变异的尺度差异性,提出了基于广义回归神经网络残余kriging的沉积物粒度组分空间预测方法,并以海州湾沉积物粒度数据为基础,分析了其在沉积物粒度组分空间预测和底质类型制图中的应用效果。结果表明,广义回归神经网络残余kriging方法能够获得比普通kriging方法更高的沉积物粒度组分空间预测精度,并且其底质类型的总体空间预测精度达到85%以上,相应的Kappa系数也在0.8以上,显示底质制图的预测类型与样本的实测类型具有较高的一致性。新方法对于开展定量化的沉积物粒度组分空间预测和底质类型制图具有参考价值。     

填充墙刚度影响的实测高层建筑增量动力分析

高层结构中的填充墙在地震作用下与周围结构构件之间的相互作用十分复杂,对建筑结构的整体抗震性能具有较大影响。然而我国规范在设计阶段通常不考虑填充墙对结构抗震性能的影响,统一采用周期折减系数来考虑其刚度变化引起的内力变化,因此准确评估填充墙对结构抗震性能的影响具有重要意义。本文在对一栋框-剪结构和一栋剪力墙高层建筑进行随机振动测试的基础上,利用Perform-3D对每栋高层分别建立了3种分析模型。其中对未考虑填充墙作用的结构模型,分别采取规范建议值和实测结果值两种方式进行周期折减。对通过添加斜撑单元来考虑填充墙作用的结构模型,利用环境激励测试识别获得的结构模态信息进行模型修正。对该3种模型进行了增量动力分析,探讨潜在危险性水平地震作用下填充墙对高层建筑抗震性能的影响。结果表明,填充墙增加了结构在弹性阶段的整体初始刚度,但随着地震动强度的增加逐渐丧失对结构刚度的贡献作用。相比考虑了填充墙作用的模型计算结果,规范建议的周期折减系数较为保守。同时研究发现,填充墙对高层框-剪结构的影响程度要比剪力墙结构大。     

黄土液化微细观特性试验研究

黄土液化演化过程的微观机理分析是液化防御的科学问题之一。通过微细观及动力学试验探索黄土液化的本质和影响因素。首先用CT细观扫描实验探索黄土渗透液化的细观变化,研究表明土体液面上升的根本原因是弱碱性盐类胶结物的吸水作用导致土样含水面整体上升;试样达到高饱和度,大孔隙周围颗粒间胶结物质破坏后有效应力为零,土层液化。粉土的孔隙尺寸和特殊的胶结物质导致高饱和度。土样微观结构的差异也会影响土的液面上升和破坏强度。针对低黏性粉土、粉质砂土及粉质黏土的三类黄土液化实验分析表明,低黏性粉土动荷加载时间更短,更易于液化,即低粘性粉土液化最为严重,粉质砂土为中等液化,粉质黏土相比其他黄土类别不易液化。电镜扫描土样微观结构参数分析表明,土颗粒周围胶结物质的化学元素比值(Ca/Fe),以及土颗粒粒径分布和孔隙尺寸(孔隙与颗粒比)均影响液化等级,可初步判断液化的强弱。     

Non-symmetrical levee breaching processes in a channel bend due to overtopping

A series of experiments were done to reveal the overtopping breaching process of non-cohesive and cohesive levees in a U-bend flume. The flood hydrograph and breaching geometry were measured and analyzed in detail. The results show that the levee breaching processes can be briefly divided into four stages: slope erosion, longitudinal headward gully-cutting, lateral erosion, and relative stabilization. For non-cohesive levees, non-symmetrical lateral development of the breach occurs throughout the four stages, and the final non-symmetrical coefficient is approximately 2.2–2.6. Larger flow discharge or higher water level can accelerate the breaching process, while coarser sands tend to accelerate the process initially but depress the process at the end. The fluvial erosion rate of a non-cohesive breach shows a power-function relation with the excess wall shear stress. For cohesive levees, a plateau forms in the breach partially blocking the flow in the first two stages. The breach flow is approximately perpendicular to the levee body, and, thus, the erosion rates of the two breach sides are almost the same. Non-symmetrical lateral development mainly occurs in the third stage when the deep gully forms. The final non-symmetrical coefficient is approximately 2.7–3.3. It is expected that these findings can provide a valuable experimental dataset and a theoretical basis for breach closure and flood alleviation.