Numerical prediction of blast‐induced stress wave from large‐scale underground explosion

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.     

非线性Rossby波所满足的三阶Zakharov方程与低频振荡

本文在基本气流具有水平切变的情况下,利用摄动法导出了非线性Rossby波所满足的三阶Zakharov方程,然后,考虑了基流具有弱切变的情况,通过使用三阶Zakharov方程研究了Rossby波列的第一类不稳定性问题。结果表明:通过非线性作用,大气中的Rossby波列可产生调制不稳定。同时,本文对这种不稳定的区域,增长率和周期进行了详细的计算,并讨论了波振幅、波数、纬度和基流切变对它们的影响,指出Rossby波列的调制不稳定可以激发30～60天的低频振荡。     

利用高光谱数据进行植被生化成分反演方法研究

高光谱数据包含着丰富光谱信息,能够定量地分析物质成分^[1]。由高光谱数据,可以运用多元逐步线性回归方法反演植被生化组分含量,从而达到监测植被生长状况的目的^[2]。本文具体介绍了多元逐步线性回归方法,及其在氮、木质素和纤维素含量估算中的应用。     

Observations of High Resolution Spectra of Comet Hale-Bopp (C/1995 O1) at the SAO of the RAS

We present the results of the preliminary study of the comet Hale-Bopp spectrum obtained April 17, 1997 by K. Churyumov and F. Mussayev with the help of the 1-meter Zeiss reflector and the echelle spectrometer (spectral resolutionλ/Δ λ ≈ 50000), CCD and the long slit, oriented along the radius-vector(“Sun-comet direction”). Energy distributions for three selected regions including the C₃, C₂ (0-0) and CN(Δ ν = 0) molecules emissions of the comet Hale-Bopp spectrum were built. The rotational lines of the CN(Δ ν = 0) band were identified. The nature of the high emission peak near λ 4020 Å in the C₃ band is discussed. The presence of the cometary continuum of the nonsolar origin is assumed.     

Seismic rays in an ellipsoidal earth model: theoretical analysis and estimation of deviations due to the ellipticity

Summary. A first-order form of the Euler's equations for rays in an ellipsoidal model of the Earth is obtained. The conditions affecting the velocity law for a monotonic increase, with respect to the arc length, in the angular distance to the epicentre, and in the angle of incidence, are the same in the ellipsoidal and spherical models. It is therefore possible to trace rays and to compute travel times directly in an ellipsoidal earth as in the spherical model. Thus comparison with the rays of the same coordinates in a spherical earth provides an estimate of the various deviations of these rays due to the Earth's flattening, and the corresponding travel-time differences, for mantle P -waves and for shallow earthquakes. All these deviations are functions both of the latitude and of the epicentral distance. The difference in the distance to the Earth's centre at points with the same geocentric latitude on rays in the ellipsoidal and in the spherical model may reach several kilometres. Directly related to the deformation of the isovelocity surfaces, this difference is the only cause of significant perturbation in travel times. Other differences, such as that corresponding to the ray torsion, are of the first order in ellipticity, and may exceed 1 km. They induce only small differences in travel time, less than 0.01s. Thus, we show that the ellipticity correction obtained by Jeffreys (1935) and Bullen (1937) by a perturbational method can be recovered by a direct evaluation of the travel times in an ellipsoidal model of the Earth. Moreover, as stated by Dziewonski & Gilbert (1976), we verify the non-dependence of this correction on the choice of the velocity law.     

密度和压缩系数的散射层析成像法

本文在速度成像的基础上研究了同时对密度和压缩系数成像的散射波层析成像法.对不同散射角度的计算可以得到一系列反演图像,拟合这些图像,从而可以有效地达到对密度和压缩系数（或速度）成像的目的.与单纯的速度成像相比,增加了反演的难度.首先是对资料的方位性要求增加;其次是对资料的利用率下降.即便如此,从对较少量的炮点和检波点资料的数值计算来看,仍取得了满意的成像结果.我们对组成字母“A”的散射体结构进行了成像计算,结果能够同时再现密度和压缩系数,成像清晰,表明了方法的可行性,并能应用于复杂结构的成像问题.     

青藏高原及邻近区域的S波三维速度结构

本文收集了WWSSN台网和我国台网中13个地震台站的长周期地震记录,用140条10-90s瑞利波频散曲线和作者提出的Tarantola-Backus面波频散层析成象方法,作了青藏高原及邻区的速度反演,得出该地区岩石层速度结构的三维图象.结果表明,1.在10-110km深度范围内,速度结构出现与大地构造特征相关的分区性,显示出四个构造单元:青藏块体、柴达木-巴颜喀拉-三江块体、塔里木块体和印度块体.2.高原内部,深度为10-70km内速度较低,莫霍界面呈不对称盆形分布,藏北那曲附近地壳厚度超过70km,高原边缘壳厚为45-50km,90-110km为高速异常,表明高原内部存在上地幔盖层.3.高原北部的班公湖断裂和东部的三江断裂系是该区重要的分界线,是岩石层结构存在明显差异的重要接触部位,可能是冈瓦纳古陆与欧亚古陆的缝合带.4.柴达木-巴颜喀拉-三江块体内部速度分布不均匀,地壳厚度由北向南从45km加深到60km;在深度90-110km存在一低速层.5.塔里木地块内速度随深度均匀增加,从地壳到上地幔110km内没有发现低速层.地壳厚度约50km.     

论裸眼井中的“共振纵波”和“共振横波”

一、引言 文献[1]和[2]指出,充流体裸眼井中存在频率确定、波速恒定、且沿传播方向不衰减的“共振纵波”和“共振横波”。本文的目的在于论证并指出,这两种波是不存在的。 二、论证 设一点源位于z=0的井轴上,则井轴上z处的声波频率响应为     

Predicting ground-motion variations at the Turkey-flat test site,California

Results of a single group participating in an international experiment are analyzed. The experiment served to verify computational predictions of the ground-motion variations due to near-surface geological effects at a site established for that purpose by the California Department of Conservation. Based on an acceleration record at a rock location, and geotechnical model of medium, records at the other locations of a nearby sedimentary deposit were predicted. A 2-D finite-difference sensitivity analysis suggested that the lateral wave-propagation effects are negligibly small, and locally 1-D computations are sufficient for the present site. Those computations are compared with observations not available to the authors during the blind prediction. Peak accelerations, peak velocities and RMS accelerations were predicted with errors less than 159%, 114% and 62%, respectively. Maxima of the response spectra were fitted within a factor of 2. The predicted and observed Husid's plots (i.e., the normalized cumulative plots of the acceleration squared) have the correlation coefficients 0.98. The detected misfits do not show any simple relation to the instrument location, component, frequency, or time.