无网格数值方法在结构物水动力研究中的应用

运用一种移动粒子半隐式法建立无网格水流数值模型,探讨该模型在结构物水动力特性研究中的可能性。在无网格条件下,粒子间的相互作用通过核子函数来表示,并通过建立物理量光滑模型、梯度模型和二阶导数模型来离散控制方程。流体的不可压缩性是通过引入粒子数密度概念实现的。该方法不受流体变形程度和固边界形状的限制,因而具有特定的适应性。运用该模型模拟了水柱倒塌后有无透空块体两种情况下右端直墙上的压力过程。计算结果表明透空块体可使直墙压力减小42%。     

Calibration of modelled mixing patterns in loess grain-size distributions: an example from the north-eastern margin of the Tibetan Plateau, China

Genetically meaningful decomposition (unmixing) of sediment grain-size distributions is accomplished with the end-member modelling algorithm. Unmixing of the loess grain-size distributions of a Late Quaternary loess–palaeosol succession from the north-eastern Tibetan Plateau indicates that the loess is a mixture of three end-members representing very fine sandy, coarse silty and medium silty loess. The unmixing approach potentially enables the unravelling of sediment fluxes from multiple dust sources, opening the way to significant advances in palaeoclimatic reconstructions from loess grain-size distribution data. However, as laser-diffraction size analysis is a volume-based technique, the proportional contributions of the modelled end-members might deviate (significantly) from weight proportions. Hence, calibration of the end-member volume proportions to weight proportions must be established before one can calculate the source-specific dust fluxes. This paper reports the findings of a sediment-mixing experiment which enables calibration of the modelled mixing patterns established for the Tibetan loess–palaeosol succession.     

极隙区纬度太阳活动高年背景电离层特征

本文利用 EISCAT-Svalbard 雷达1999年3月至2003年2月期间观测数据,对太阳活动高年前后极区电离层 F 区电子密度随高度和地方时的二维分布进行统计分析,研究极隙/极盖区背景电离层特征,并与 IRI-2001模式预测结果进行比较。分析结果表明,在 ESR 雷达所在纬度,不存在“冬季异常”现象;但首次发现,此处电离层冬季日变化存在磁午夜前主极大。另外,无论在什么季节,白天时段 ESR 雷达上空电子密度都在磁正午附近有一极大值。与 IRI 模式的比较表明,模式预测的 F 区电子密度与 ESR 实测结果存在较大偏差,在500km 高度以上的顶部电离层和冬季,差别尤其明显。简要讨论了软粒子沉降对造成电子密度日变化磁正午极大的作用和冬季磁午夜前主极大的形成机理。     

沁水盆地石炭-二叠纪岩相古地理及聚煤作用研究

晋东南沁水盆地是目前我国进行煤层气勘探开发的主要地区之一，煤层气储层主要是石炭系和二叠系的煤层，这些煤层气储层的厚度明显受沉积古地理控制。本文通过对露头及钻井剖面岩石学和沉积相研究,对该区主要含煤岩系-山西组和太原组进行了岩相古地理和煤储层聚集控制因素分析。太原组以石灰岩、铝土质泥岩、粉砂岩、粉砂质泥岩及砂岩为主，厚44.90～193.48 m，含多层可采煤层，总厚0～16.89 m，平均7.19 m，在太原组沉积期，研究区北部发育下三角洲平原相，煤层相对较厚，中部和南部为泻湖相，煤层相对较薄，东南角主要为滨外碳酸盐陆棚相，在晋城一带障壁砂坝相分布区，煤层亦较厚；山西组以砂岩、粉砂岩、泥岩为主，石灰岩仅在局部地区见到，该组厚18.60～213.25 m，含可采煤层1-2层，总厚0～10.0 m，平均4.20 m，在山西组沉积期，北部以下三角洲平原分流河道相为主，中部和南部以分流间湾相为主，东南部则以河口砂坝相为主，厚煤带都位于中部和南部三角洲分流间湾地区。总之，太原组富煤地带多与砂岩富集带相吻合，位于北部下三角洲平原和南部障壁砂坝地区，而山西组厚煤带大都位于南部三角洲分流间湾地区。     

伊朗中东部Ghaen-Birjand地区上地壳各向异性研究c

通过对1997年5月10日Zirkuh ( Ghaen-Birjand )破坏性地震若干余震横波分裂现象的观测，研究了该地区上地壳各向异性特征. 采用质点运动图和纵横比两种不同手段获取了横波分裂参数. 所选余震记录显示了明显的横波分裂，说明该地区介质具有强烈的各向异性. 通过质点运动的方法得出:快横波的偏振方向为22N19E；且快、慢横波的延时之差为6516 ms. 用纵横比的方法得到的结果为:快横波的偏振方向35N18E, 且快、慢横波的延时之差为(4910)ms. 对于厚约5 km各向异性均匀分布的水平层， 我们采用了一个沿应力分布的裂纹模型，并采用垂直于N22E排列, 且密度为0.01的裂纹情况来解释本文结果. 由于研究范围是上地壳，因此假设裂纹被水填充. 最后通过Hudson模型，对裂纹密度为0.005, 0.01和0.03的模型，给出横波速度随传播方向与裂纹对称轴夹角的变化曲线. 结果表明，当横波的偏振平行于裂纹表面时波速不随角度变化；而偏振垂直于裂纹表面时，速度曲线则有明显变化.      

Experimental investigation of the velocity of a sand cloud blowing over a sandy surface

The velocity of a wind‐blown sand cloud is important for studying its kinetic energy, related erosion, and control measures. PDA (particle dynamics analyser) measurement technology is used in a wind tunnel to study the probability distribution of particle velocity, variations with height of the mean velocity and particle turbulence in a sand cloud blowing over a sandy surface. The results suggest that the probability distribution of the particle velocity in a blowing sand cloud is stochastic. The probability distribution of the downwind velocity complies with a Gaussian function, while that of the vertical velocity is greatly complicated by grain impact with the bed and particle–particle collisions in the air. The probability distribution of the vertical velocity of ?ne particles (0·1–0·3 mm sands) can be expressed as a Lorentzian function while that of coarse particles (0·3–0·6 mm sands) cannot be expressed by a simple distribution function. The mean downwind velocity is generally one or two orders greater than the mean vertical velocity, but the particle turbulence in the vertical direction is at least two orders greater than that in the downwind direction. In general, the mean downwind velocity increases with height and free‐stream wind velocity, but decreases with grain size. The variation with height of the mean downwind velocity can be expressed by a power function. The particle turbulence of a blowing sand cloud in the downwind direction decreases with height. The variations with height of the mean velocity and particle turbulence in the vertical direction are very complex. It can be concluded that the velocity of a sand cloud blowing over a sandy surface is mainly in?uenced by wind velocity, grain impact with the bed and particle–particle collisions in the air. Wind velocity is the primary factor in?uencing the downwind velocity of a blowing sand cloud, while the grain impact with the bed and particle–particle collisions in the air are the primary factors responsible for the vertical velocity. Copyright © 2004 John Wiley & Sons, Ltd.     

用多阶Mueller矩阵解对非均匀地表植被生物量和土壤湿度的迭代反演

根据一层随机非球形粒子矢量辐射传输方程的一阶与高阶Mueller矩阵解的数值模拟及其差异 ,提出了在一些参数已知的条件下 ,用水平与垂直同极化和交叉极化的后向散射测量反演植被覆盖下地表湿度和植被散射元单位面积占空比的迭代方法。利用模拟测量数据以及机载合成孔径雷达 (AirSAR)实测数据 ,对不同土壤湿度和单位面积占空比的情况进行了反演试验     

Numerical study of characteristics of underground blast induced surface ground motion and their effect on above-ground structures. Part I. Ground motion characteristics

Current codes of practice in assessing the blast ground motion effect on structures are mainly based on the ground peak particle velocity (PPV) or PPV and the principal frequency (PF) of the ground motion. PPV and PF of ground motion from underground explosions are usually estimated by empirical formulae derived from field blast tests. Not many empirical formulae for PF, but many empirical formulae for PPV are available in the literature. They were obtained from recorded data either on ground surface or in the free field (inside the geological medium). Owing to the effect of surface reflection, blast motions on ground surface and in the free field are very different. But not many publications in the open literature discuss the differences of blast motions on ground surface and in the free field. Moreover, very few publications discuss the blast ground motion spatial variation characteristics. As ground motion directly affects structural responses, it is very important to study its characteristics in order to more reliably assess its effects on structures. In this paper, a validated numerical model is used to simulate stress wave at a granite site owing to explosion in an underground chamber. Using the simulated stress wave, the relations such as PPV and PF attenuation as well as spatial variation of motions on ground surface and in the free field are derived. Discussions on the differences of the characteristics of surface and free field motions are made. Results presented in this paper can be used in a more detailed assessment of ground motion effect on structures.     

新疆库车坳陷侏罗系碎屑岩中粘土矿物粒径研究

选取库车坳陷侏罗系阿合组—齐古组的碎屑岩样品,采用跟进与离心相结合的方式将岩样进行5级分离,对分离后岩样SEM(EDS)分析观察到不同粒径粘土矿物的形貌特征。在20～2μm范围内的样品中,可见片状碎屑粘土矿物;在<2μm的样品中可见自生粘土矿物,以集合体形式存在。对不同粒级的粘土XRD分析表明,该区碎屑岩的主要粘土矿物类型为伊利石、高岭石、绿泥石和伊蒙混层,伊蒙混层主要出现在<0.2μm范围内。高岭石的含量随着岩样粒度的变粗呈现上升的趋势,而伊利石则相反。在气候温暖潮湿的早侏罗世和中侏罗世早期,高岭石的含量相对高;在气候干旱的中侏罗世晚期及晚侏罗世,高岭石的含量减少,伊利石则呈现相反趋势。区分不同成因的粘土矿物可以为进一步研究这一地区的沉积环境、成岩演化历史提供有力依据。     

A NUMERICAL MODEL OF MIXED CONVECTIVESTRATIFORM CLOUD

A 2-D slab-symmetric model of mixed convective-stratiform cloud is developed bysuperimposing convective cloud-size field on the convergence field,in order to simulate and studythe mixed clouds consisting of stratiform cloud and convective cloud.A deep convective,anelasticand conservative system of equations with basic variables(V,θ,π')is solved by a new method tocalculate dynamic field.The water substance in the cloud is divided into 6 categories and themicrophysical processes are described in spectrum with two variable parameters and morereasonable particle number/size distributions.To compare with measured radar echo intensity andstructure,the model may calculate echo intensity of the model cloud observed by radar.