城市垃圾填埋场振动台试验的数值模拟研究

对填埋场振动台模型试验进行非线性数值分析,并与试验结果相互补充、印证,进一步考察填埋场的地震响应规律。研究表明,特殊的夹层结构形式导致水平地震作用下城市垃圾填埋场存在两个主要的响应频率,在这两个主要频率附近模型顶部取得峰值加速度响应;不考虑其他因素的影响,填埋场防渗层的地震永久位移与输入地震动的水平位移振幅近似呈正比关系;坡比是影响填埋场覆盖层和衬垫层永久位移最主要的因素。     

模拟不同气候条件下碳酸盐岩风化作用的淋溶实验研究

通过模拟干热、湿热与干冷3种气候条件,以饱和CO2水作为淋溶液,对黔中岩溶区3条碳酸盐岩风化壳岩-土界面的岩粉层试样进行了淋溶实验(其中干热、湿热条件下淋溶到残余酸不溶物阶段),对淋出液的pH值以及主要造岩元素的浓度进行了动态分析。结果表明,碳酸盐岩风化壳岩-土界面由岩到土的转变过程中,伴随碳酸盐的溶蚀,酸不溶物已表现出明显的风化倾向。碳酸盐的溶蚀强度表现为干冷>干热>湿热的变化趋势。温度低,碳酸盐的溶解速率大;排水条件好,碳酸盐溶解释放的Ca、Mg易随风化流体排出体系。对于酸不溶物组分,淋溶实验中有:(1)K、Na、Mg、Si、P的载体矿物风化强度在干热条件下最大(至于Ca,由于方解石与白云石均是其主要的载体矿物,酸不溶物相中的Ca难以在淋出液中有效识别)。在干冷与湿热之间,K、Na、Mg等盐基离子的溶出能力大多表现为湿热>干冷,指示了温度对盐基离子释放强度的重要制约作用;而Si和P未表现出一致的变化趋势,可能源于淋溶体系微环境的差异。(2)Fe与Mn表现出弱迁移性。其中,Fe在干冷环境下淋出强度最弱,说明温度是制约含Fe矿物分解速率的重要因素。而Fe在干热与湿热之间,以及Mn在3种淋溶条件下,未呈现出一致的变化趋势。(3)Al和Ti在碳酸盐岩风化过程中表现出强烈的惰性。      

俯伪斜分段密集采煤法三维有限元数值模拟研究

俯伪斜分段密集采煤法,作为一种较新的采煤方法,在我国技术上还不很成熟,因此,首次针对京煤集团木城涧矿大台井的地质条件和赋煤状况,通过简化建立三维立体有限元模型,模拟计算了该矿煤层顶底板在开采前后运动变化情况。结果表明:煤层倾角、推采距离和开采巷道数三者是影响煤层顶底板运动规律的重要因素,得到了表征顶底板运动规律的顶底板法向最大应力和法向最大位移与三影响因素关系曲线,详细阐述了三影响因素下的煤层顶底板运动规律;同时,对现场开采中支撑压力大小,特别是在三影响因素下,超前支撑压力的作用范围与应力峰值位置进行了数值计算,并总结了其影响规律;最后,指出了该数值分析方法对现场开采将具有十分重要的应用价值和指导意义。     

Flexible membranes anchored to the ground for slope stabilisation: Numerical modelling of soil slopes using SPH

An alternative modelling for flexible membranes anchored to the ground for soil slope stabilisation is presented using Smoothed-Particle Hydrodynamics to model the unstable ground mass in a soil slope, employing a dynamic solve engine. A regression model of pressure normal to the ground, q_sim, and also membrane deflection, f_sim, have been developed using Design of Experiment. Finally, a comparison between the pressure obtained from numerical simulation and from a limit equilibrium analysis considering infinite slope has been carried out, showing differences in the results, mainly due to the membrane stiffness.     

Pull-down basin in the central part of Japan due to subduction-induced mantle flow

The driving force for the basin subsiding against isostatic balance in and around Lake Biwa in the Kinki district, Japan is discussed. The lake region is characterized by strong negative Bouguer anomalies, especially by a steep horizontal gradient zone of gravity anomaly running along the western margin of the lake. The large negative anomaly (>50 mgal) cannot be explained by low-density sediments beneath it. A down-warping structure extending to the Moho depth should be taken into account. This conjecture has been strongly supported by a short-period receiver function imaging, which shows a clear offset of about 8 km for the Moho discontinuity under the steep gravity gradient zone.A question arises as to what is the driving force to create such a large down-warping structure. We consider that the subduction of the shallow-dipping slab under the region (Philippine Sea Slab) may cause crustal deformation by dragging the viscous mantle downward. In order to verify this model, we simulated the induced mantle flow due to the subduction of the Philippine Sea Slab and the pressure distribution on the crust–mantle boundary. This numerical experiment showed that the induced flow makes a strong negative pressure zone under the lake region if the slab has a vertical offset along the direction of subduction. This offset of the slab is consistent with plate models deduced from hypocentral distributions and Sp phases of the deep-focus earthquakes.     

单一松软煤层瓦斯抽采工艺数值模拟

借鉴油气藏数值模拟理论、技术和方法,对单一松软煤层的瓦斯抽采工艺进行了研究分析。根据单一松软煤层特点,提出了单一松软煤层的几何模型技术——单直径球型孔隙模型,在此基础上建立了煤层瓦斯抽采的地质和数理模型,并根据网格差分方法推导出数值解。以实际生产数据对新模型及其数值解进行验证,结果表明该模型可以用于指导单一松软煤层瓦斯抽采工艺的确定。     

The Zoulang Nanshan Caledonian Subduct ion Complex in the Northern Qilian Mountains and Its Dynamics

The front of the Zoulang Nanshan Caledonian volcanic island arc zone in the northern Qilian Mountains is a forearc accretionary terrane, composed of multiple accretionary volcanic island arcs, flysch accretionary wedges,high-pressure metamorphosed detachment zones and remnants of ophiolites. It resulted from the northeastward subduction of the Early Palaeozoic Qilan oceanic crust beneath the Alxa block. High-pressure metamorphism, which occurred during the subduction, progressed through three stages: the initial stage of medium T-high P,the main stage of temperature decrease and pressure increase, and the lag stage of pressure decrease and temperature increase. Finally the paper presents a retrotrench subduction dynamic model indicative of northward subduction of the central Qilian block and southward accretion of the Alxa block during the period of 450-500 Ma.     

Geostatistical modeling of the geological uncertainty in an iron ore deposit

This paper addresses the problem of quantifying the joint uncertainty in the grades of elements of interest (iron, silica, manganese, phosphorus and alumina), loss on ignition, granulometry and rock types in an iron ore deposit. Sampling information is available from a set of exploration drill holes. The methodology considers the construction of multiple rock type outcomes by plurigaussian simulation, then outcomes of the quantitative variables (grades, loss on ignition and granulometry) are constructed by multigaussian joint simulation, accounting for geological domains specific to each quantitative variable as well as for a stoichiometric closure formula linking these variables. The outcomes are validated by checking the reproduction of the data distributions and of the data values at the drill hole locations, and their ability to measure the uncertainty at unsampled locations is assessed by leave-one-out cross validation.Both the plurigaussian and multigaussian models offer much flexibility to the practitioner to face up to the complexity of the variables being modeled, in particular: (1) the contact relationships between rock types, (2) the geological controls exerted by the rock types over the quantitative variables, and (3) the cross-correlations and stoichiometric closure linking the quantitative variables. In addition to this flexibility, the use of efficient simulation algorithms turns out to be essential for a successful application, due to the high number of variables, data and locations targeted for simulation.     

非均匀月壤介质的被动微波辐射传输模拟

基于非均匀月壤物理模型和辐射传输方程,模拟月壤介质中的微波辐射传输特性,探讨频率、月壤厚度等与月表亮温的关系。结果表明:在低频段,月壤微波辐射亮温的动态变化范围较大,可探测的月壤厚度大,3 GHz时的最大可探测月壤厚度达12.4 m;在高频段对应的可探测月壤厚度较小,特别是从50GHz往后的频率段内,最大可探测月壤厚度均小于2 m。不同频率的亮温-厚度变化曲线没有交叉点,且频率越高,所能探测的月壤厚度越小。根据模拟结果,建立了月壤厚度与亮温的查找表。基于查找表,利用单个波段的亮温数据即可得到月壤厚度信息。     

液化土层地震动模拟计算方法及验证

针对目前土层反应分析方法难以模拟液化土层地震动时程的缺欠,在工程力学研究所原有有效应力分析程序基础上提出了一个改进的计算方法并进行了验证。改进方法中,通过每个应力循环模拟土的非线性进程以及液化导致的土刚度衰减过程,并引入了新建立的适于非均等固结随机地震荷载作用下的孔压增长模型。将计算方法与实际地震记录及大型振动台实验结果进行了对比,结果表明：在峰值加速度、液化后波形变化、时频曲线及加速度反应谱等主要特征上,计算结果均与现场实际记录及振动台实验基本一致;在液化引起地震动特征变化上,计算出的液化较非液化土层加速度反应谱的增量与实验结果一致。说明所提出的改进方法可以用于液化土层地震动的模拟。