浙江省梅汛期洪涝灾情分析和预测

通过所构设的雨涝强度、雨涝影响面积等统计量来估计灾情，有较高的正确率，能估计出极端情况，为灾害评估工作提供了新的思路。     

小兴安岭-张广才岭地区晚古生代至中生代花岗岩的成因及其地质意义

中国东北地区显生宙花岗岩的成因一直是中亚造山带东段研究的热点之一,尤其是小兴安岭-张广才岭地区的花岗岩其成因及形成的大地构造背景一直存在较大争议。本文新获得小兴安岭-张广才岭铁力和依兰地区的二长花岗岩LA-ICP-MS锆石U-Pb年龄分别为188±1 Ma和257±3 Ma。地球化学数据显示,两地区的二长花岗岩均为Ⅰ型花岗岩并且富集大离子亲石元素(Th和U等)和轻稀土元素,亏损高场强元素(Nb和Ta等)和重稀土元素具有弧岩浆岩的地球化学特征。锆石Hf同位素数据显示,铁力地区二长花岗岩的岩浆源区可能来自于中新元古代的下地壳部分熔融。综合前人在小兴安岭-张广才岭地区已发表的花岗岩类岩石的地质年代学和地球化学数据,初步推测小兴安岭-张广才岭地区在晚古生代至中生代期间处于活动大陆边缘环境。同时,小兴安岭-张广才岭地区晚古生代至中生代的岩浆岩具有自东向西形成时代逐渐变年轻的趋势,这可能是由于东侧的牡丹江洋在晚古生代至中生代期间俯冲角度逐渐变缓造成的。     

云南哀牢山-红河剪切带新生代富碱侵入岩中榍石年代学和地球化学特征

榍石作为副矿物在哀牢山-红河剪切带新生代富碱侵入岩中广泛存在。原位获取榍石矿物内部微量元素、U-Pb年龄和Sm-Nd同位素的空间变化对获取矿物和岩石的成因演化信息具有重要意义。本文使用四级杆/多接收电感耦合等离子体质谱(Quadrupole/Multi-Collector Inductively Coupled Plasma Mass Spectrometry,Q/MC-ICP-MS)与激光剥蚀系统(Laser Ablation,LA)联用,对哀牢山-红河剪切带5个富碱侵入岩体(桃花岩体、宁蒗-永胜岩体、哈播岩体、铜厂岩体和十里村岩体)中榍石开展了微区原位微量元素、U-Pb年代学和Sm-Nd同位素研究。微量元素分析结果表明,三江富碱侵入岩中榍石为岩浆成因,亏损Rb、Ba、Pb、Sr等大离子亲石元素,富集Th、U、Nb、Ta、Zr、Hf等高场强元素。榍石的稀土配分图均表现为明显右倾,不具有或具有弱的Eu负异常。与云南北衙、马厂箐矽卡岩矿床中的热液榍石相比,本文榍石在稀土元素组成上,具有较高的稀土元素总量、较高的Th/U、LREE/HREE和Ce/Ce*比值,具有较低的Eu/Eu*、Nb/Ta、Zr/Hf比值。微区原位LA-Q-ICP-MS U-Pb定年结果表明,研究区富碱岩体中榍石结晶年龄在32.5～37.9Ma之间,代表了岩体形成时代,与三江地区哀牢山-红河剪切带及其附近新生代富碱岩浆活动高峰期(～35Ma)一致,属于青藏高原晚碰撞期岩浆作用的产物。榍石微区原位LA-MCICP-MS Sm-Nd同位素分析结果显示,榍石颗粒的Nd同位素组成均一,表明榍石结晶过程中寄主岩浆的Nd同位素组成没有发生明显变化。各个富碱岩体之间的榍石Nd同位素组成变化范围在-6.8～-2.1之间,与全岩的同位素特征一致,表明榍石的原位微区Sm-Nd同位素可以作为富碱侵入岩研究中有效的示踪手段之一。     

山东荣成湾月湖地区的潟湖-潮汐汊道体系

泻湖－潮汐汊道体系为一复杂的反馈体系,其中任一参数的改变即可导致整个体系的复杂反映。以月湖为例探讨了泻湖地貌演化的过各,及人工作用泻湖各主要参数的变化规律;提出了月湖目前由于人工署堵特别是残余坝基对潮流的限制作用,致使Ｐ／Ａ值较高而饱满系数Ｋ值偏代,泻湖－潮汐汊道体系稳定性较差的问题。     

The sedimentary characteristics and formation mechanism of shell ridges along the southwest coast of Bohai Bay

The present paper studies the sedimentary characteristics and mechanism of the shell ridges on the southwestern coast of Bohai Bay, which are the largest and have the highest shell (or shell fragment) content in the world. These shell ridges are composed of two sedimentary subfacies: the ridge subfacies and the ridge infill subfacies with different inner textures and sedimentary structures. The ridge subfacies primarily consists of fresh shells and/or shell fragments with parallel beddings and high-angle oblique beddings. The ridge infill subfacies consists of finer shell fragments and silty sand with parallel beddings oblique to land. The evolution of the shell ridge is controlled by the accumulation of substantial shell material,the relative slow erosion of shoreline, storm waves, winds and the shift of river routes. Wind tunnel tests indicate that the critical movement velocity of shell is lower than that of quartz sand of the same grain size. Deltaic progradation alternates with shell ridge growth. While the Yellow River empties into the sea through this area, accompanied by deltaic progradating, the shell material is scarce, which is unfavorable to the growth of the shell ridge. Conversely, erosion occurs along an abandoned delta coast where a shell ridge may develop, e.g. the shell ridge on the southwest of Bohai Bay.     

用于颗粒土微观力学行为试验的微型三轴试验仪

颗粒土的微观力学行为（如土颗粒的运动与破碎等）决定着其宏观应力-应变行为,如应变局部化、应力硬化等。为研究颗粒土的微观力学行为,开发了一台微型三轴试验装置。它的轴向加载系统由伺服控制的步进转动马达与涡轮传动的减速器组成,围压由GDS压力控制器提供,压力室采用高透光率,高强度的轻质材料制成。借助于X射线显微CT及图像处理分析技术,该装置能实现对干砂土微尺寸试样（直径为8 mm,高度为16 mm）在三轴剪切条件下微观特性的无损检测。采用该三轴试验装置对粒径为0.60～1.18 mm 的LBS（Leighton Buzzard sand）试样在1.5 MPa的围压下进行了试验。结果显示,试验装置测得应力-应变曲线合理,显微CT 图像特征清晰,能够用于颗粒土体微观土力学行为的试验研究。     

Two Improvements on Numerical Simulation of 2-DOF Vortex-Induced Vibration with Low Mass Ratio

Till now, there have been lots of researches on numerical simulation of vortex-induced vibration. Acceptable results have been obtained for fixed cylinders with low Reynolds number. However, for responses of 2-DOF vortex-induced vibration with low mass ratio, the accuracy is not satisfactory, especially for the maximum amplitudes. In Jauvtis and Williamson''s work, the maximum amplitude of the cylinder with low mass ratio m*=2.6 can reach as large as 1.5D to be called as the “super-upper branch”, but from current literatures, few simulation results can achieve such value, even fail to capture the upper branch. Besides, it is found that the amplitude decays too fast in the lower branch with the RANS-based turbulence model. The reason is likely to be the defects of the turbulence model itself in the prediction of unsteady separated flows as well as the unreasonable setting of the numerical simulation parameters. Aiming at above issues, a modified turbulence model is proposed in this paper, and the effect of the acceleration of flow field on the response of vortex-induced vibration is studied based on OpenFOAM. By analyzing the responses of amplitude, phase and trajectory, frequency and vortex mode, it is proved that the vortex-induced vibration can be predicted accurately with the modified turbulence model under appropriate flow field acceleration.     

Simulation of interactions among multiple debris flows

Adjacent debris flows may interact in many ways: two or more concurrent debris flows may merge; one debris flow can run out over an existing debris flow fan. Such interactions may cause debris flow properties to change in the mixing process as well as more severe adverse effects than those caused by a single debris flow. This paper aims to investigate the interactions among channelized debris flows originated from adjacent catchments. Both concurrent and successive debris flows are considered. If several debris flows originate from different locations concurrently and merge, the volumetric sediment concentration (i.e., the ratio of the volume of solid material to the total volume of debris flow), C _v, is a good index to capture the mixing process of these debris flows. The change in C _v reflects where mixing occurs and the mixing degree. The debris flow properties (e.g., yield stress and dynamic viscosity) evolve in the mixing process and can be captured by the change in C _v. The debris flow with a larger volume dominates the mixing process, and the properties of the mixed debris flow are more similar to those of the larger debris flow. The inundated areas and runout distances of successive debris flows are smaller than those of concurrent debris flows of the same total volume due to the smaller scales of the individual events and blockage by the earlier debris flows. However, the deposit depth in the interacting part of the debris flow fans of successive debris flows can be much larger than that of concurrent debris flows, leading to more destructive cascading hazards (e.g., the formation of debris barrier lakes). The sequence of successive debris flows not only significantly influences the runout characteristics of the debris flows but also substantially affects the cascading hazards.     

Evaluation of CBM reservoir fracture characteristics by azimuthal AVA inversion

Exploitation of coalbed methane (CBM) is directly affected by the strike and density of fractures, particularly because fractures can act as conduits for water and gas (after desorption). To realize the azimuthal amplitude versus angle (AVA) inversion for a CBM reservoir, we treated reservoirs whose fractures are aligned vertically as being equivalent to a horizontal transverse isotropic medium and simplified the P-P wave reflection coefficient formula to an approximately linear equation. Then, we adopted a modified damping least squares algorithm to solve the azimuthal AVA inversion problem. With edited macro-bins, corrected residual moveout, and normalized amplitude, the azimuthal gathers exhibited unambiguous azimuthal AVA characteristics. The #3 CBM reservoir in the Qinshui Basin was selected for a case study. The density and strike of the reservoir fractures were obtained and analyzed, and the reservoir was divided into three classes according to fracture density. Finally, in comparison with the gas production data, we confirmed that the inversion results and reservoir classification were both reasonable and credible.     

Centroid sliding pyramid method for removability and stability analysis of fractured hard rock

This paper proposes a new method using centroid sliding pyramid (CSP) to identify the removability and stability of fractured hard rock in tunnel and slope engineering. The new method features two geometrical and topological improvements over the original key block method (KBM). Firstly, all the concave corners are considered as starting points of cutting process when a concave block is divided into a set of convex blocks in the original KBM. Only the concave corners formed by two joint planes are used for partitioning a concave block in the presented method and concave corners with free planes are excluded. Secondly, joint pyramid for removability computation in the original KBM is generated using all of the joint planes, while CSP is calculated only from the joint planes adjoining the free planes. The cone angle θ of CSP is the vectorial angle formed by the two candidate sliding surfaces of this CSP. Removability analysis of a block is transformed into calculating the cone angle of CSP. The geometrical relationship is simplified, and data size for removability computation is reduced compared with the original KBM. The provided method is implemented in a computer program and validated by examples of fractured rock slopes and tunnels.