The propagation of a single shear crack: A displacement discontinuity model

Numerical studies using the displacement discontinuity method show that a single shear crack under compression propagates in its own direction, because such propagation results in the maximum release of strain energy. The methods of linear elastic fracture mechanics may not be used for such a closed crack, and the stress intensity factors are meaningless in that case. Laboratory observations of propagation by means of kinks at an angle of approximately 70° to the crack may be due to heterogeneities, to the effect of a preexisting crack, to end effects, to microcracking, or to some combination of these factors. Such kinks may thus be local phenomena which cannot release most of the strain energy, and are not incompatible with our numerical results which are based on a global energy balance.     

About the specification of erosion flux for soft stratified cohesive sediments

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.     

Blind prediction of the site effects at Ashigara Valley,Japan, and its comparison with reality

Weak and strong ground motions were numerically predicted for three stations of the Ashigara Valley test site. The prediction was based on the records from a rock-outcrop station, one weak-motion record from a surface-sediments station, and the standard geotechnical model. The data were provided by the Japanese Working Group on the Effects of Surface Geology as a part of an international experiment. The finite-difference method for SH waves in a 2-D linear viscoelastic medium (a causalQ model) was employed.Comparison with the real records shows that at two stations the predictions fit better than at the third one. Strangely, the two better predictions were for stations situated at larger distances from the reference rock station (one station was on the surface, the other in a borehole). The strong ground motion (the peak acceleration of about 200 cm s^–2) was not predicted qualitatively worse than the weak motion (8 cm s^–2). A less sophisticated second prediction (not submitted during the experiment), in which we did not attempt to fit the available weak-motion record at the sedimentary station, agrees with the reality significantly better.     

Modelling of the vertical fluxes of nitric acid,ammonia, and ammonium nitrate

Results from numerical investigations regarding the exchange of HNO₃, NH₃, and NH₄NO₃ between the atmosphere and the biosphere are presented. The investigations were performed with a modified inferential method which is based on the generally accepted micrometeorological ideas of the transfer of momentum, sensible heat and matter near the Earth's surface and the chemical reactions among these nitrogen compounds. This modified inferential method calculates the micrometeorological quantities (such as the friction velocity and the fluxes of sensible and latent heat), the height-invariant fluxes of the composed chemically conservative trace species with group concentrationsc ₁=[HNO₃]+[NH₄NO₃] (total nitrate),c ₂=[NH₃]+[NH₄NO₃] (total ammonia), andc ₃=[HNO₃]-[NH₃] as well as the fluxes of the individual nitrogen compounds. The parameterization of the fluxes is based on the flux-gradient relationships in the turbulent region of the atmospheric surface layer. The modified inferential method requires only the data of wind velocity, temperature, humidity and concentrations (HNO₃, NH₃, and NH₄NO₃) measured at a reference height by stations of a monitoring network.     

重要找矿指示元素的野外快速分析法——金的野外快速分析

近年来,作者发展了七种以上重要找矿指示元素的野外快速分析法,此乃其一。文中提出一种新的金样冷浸技术,引进了新的光导光度检测方法,样品经溴化钠—硫酸—过氧化氢冷浸,用经丙酮—MIBK处理的泡沫塑料富集金,经混合解脱剂解脱并直接取解脱液用改进的微珠析出法或用光导光度法测定。取10g样品,可测定0.0005-50g/t范围内的金。方法简便、快速,可在普查分队驻地进行,每一工作日可测定80多个样品,由于废弃了王水分解样品和泡沫塑料的灰化,具有广阔的应用前景。     

山东省日照市蓝色经济区工程建设适宜性评价

本文以山东省日照市蓝色经济区为研究区, 通过野外地质调查和室内综合分析, 选取场地稳定性、工程地质条件、地下水条件, 地形排水条件、地质灾害及其他因素6个一级评价因子, 构造活动、场地土类型、建筑场地类别等13个二级评价因子, 利用专家-层次分析定权法进行区内工程建设适宜性评价, 得到了相应工程建设评价分区图。评价结果表明, 日照市场地工程建设适宜性总体较好, 大部分地区适宜进行高层建筑工程建设, 占84.50%; 地下空间开发利用条件优越, 全区多数地区属于城市地下空间建设适宜区, 适宜性差及不适宜区主要分布在付疃河两岸附近和海岸附近、大型高速公路和铁路沿线一带区域及城区; 轨道交通建设适宜性良好, 小部分不适宜区主要集中在付疃河下游夹仓和小古镇一带及其沿海地区。评价结果可有效地服务于该地区规划布局、建设发展和工程选址, 具有重要的理论和实践指导意义。     

南海海盆区莫霍面分布规律及其对深部钻探的意义

钻遇莫霍面是人类一直以来的梦想。深海海底是地球上离莫霍面最近的地方，目前有研究推测南海是世界上莫霍面深度最浅的海域之一，但缺乏足够的直接证据。深反射地震探测可以直接揭示岩石圈的构造形态，是莫霍面探测的重要手段。本文基于长达15000 km的深反射多道地震剖面的解释、处理、制图和分析，结合前人的研究，形成了南海海盆区莫霍面反射特征和空间分布的初步认识。① 南海东部次海盆南部早期经历了较快速扩张，岩浆供应充足，受扩张停止后岩浆活动影响较小，基底平坦，地质构造相对简单，同时洋壳地震速度结构不存在异常，且有较强的广角莫霍面反射波和可识别的地幔顶部折射波，具备莫霍面钻探的基本条件。② 南海海盆不同区域的莫霍面反射强度存在较大差异。其中东部次海盆莫霍面反射最为强烈且清晰，西北次海盆次之，西南次海盆仅有零星出现的清晰莫霍面反射且可信度不高。③ 识别南海海盆区莫霍面地震反射长度超过3500 km，首次形成了海盆区深度域莫霍面地震反射空间分布图。与重力反演的莫霍面深度相比，利用深反射多道地震计算的莫霍面深度细节更为丰富，并且可以在垂向上清晰刻画莫霍面的结构。整体上，南海海盆区莫霍面地震反射强烈和可信度高的区域中，深度较浅的区域之一是东部次海盆南部，最浅处仅约9. 5 km，其中水深4. 01 km，洋壳厚度仅5. 54 km。综合判断，东部次海盆南部是南海重要的莫霍面钻探备选区，这对南海莫霍面钻探选址具有重要意义。     

随钻测录井地质导向二维分解实时绘图方法

根据随钻测录井实时地质导向和大斜度、水平井评价成图技术需求，针对传统绘图方法存在的弊端，提出了将测录井信息、井眼轨迹和地质模型进行二维分解的实时绘图方法。针对二维分解绘图方法绘图时空复杂度较高的问题，给出了不同事件驱动下的局部实时计算和拷屏重绘算法，控制了对CPU和内存的消耗、提高了绘图效率，消除了实时绘图的闪烁和卡顿现象。应用实例表明，二维分解实时绘图方法能够实现大尺度随钻测录井地质导向图形的流畅、无卡顿实时绘图，可提高大斜度、水平井储层模型评价的刻画精度和时效。     

Numerical analysis of the symmetric methods

Aimed at the initial value problem of the particular second-order ordinary differential equations,y =f(x, y), the symmetric methods (Quinlan and Tremaine, 1990) and our methods (Xu and Zhang, 1994) have been compared in detail by integrating the artificial earth satellite orbits in this paper. In the end, we point out clearly that the integral accuracy of numerical integration of the satellite orbits by applying our methods is obviously higher than that by applying the same order formula of the symmetric methods when the integration time-interval is not greater than 12000 periods.     

Infinitesimal orbits around lagrange points in the elliptic,restricted three-body problem

This study presents a method of obtaining asymptotic approximations for motions near a Lagrange point in the planar, elliptic, restricted three-body problem by using a von Zeipel-type method. The calculations are carried out for a second-order escape solution in the proximity of the equilateral Lagrange point, L ₄, where the primaries' orbital eccentricity is taken as the small parameter .