Hybrid seismic modeling based on discrete-wave number and finite-difference methods

Any calculation of seismic wave propagation comprising the seismic source, the travel path, and the receiver site in a single finite-difference (FD) model requires a considerable amount of computer time and memory. Moreover, the methods currently available for including point sources in the 2D FD calculations are far-field approximations only. Therefore we have developed a new hybrid method for treating the seismic wave fields at localized 2D near-surface structures embedded in a 1D background medium, and excited by a point source. The source radiation and propagation in the background model is solved by the discrete-wave number (DW) method, while the propagation in the local 2D structure is calculated by the FD method. The coupling between the two sets of calculations is performed on a rectangular excitation box surrounding the local structure. We show the usefulness of the method in ground-motion studies where both near-field source effects and local site effects are important. Technical problems connected with the inconsistency between the 3D source radiation and the 2D FD calculation are minor for the relatively distant in-plane point explosive sources, but are more serious for the in-plane dislocation sources.     

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%; 地下空间开发利用条件优越, 全区多数地区属于城市地下空间建设适宜区, 适宜性差及不适宜区主要分布在付疃河两岸附近和海岸附近、大型高速公路和铁路沿线一带区域及城区; 轨道交通建设适宜性良好, 小部分不适宜区主要集中在付疃河下游夹仓和小古镇一带及其沿海地区。评价结果可有效地服务于该地区规划布局、建设发展和工程选址, 具有重要的理论和实践指导意义。     

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

根据随钻测录井实时地质导向和大斜度、水平井评价成图技术需求，针对传统绘图方法存在的弊端，提出了将测录井信息、井眼轨迹和地质模型进行二维分解的实时绘图方法。针对二维分解绘图方法绘图时空复杂度较高的问题，给出了不同事件驱动下的局部实时计算和拷屏重绘算法，控制了对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 .