Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.     

MODELING OF THE DIVERSION CHANNEL CLOSURE FOR THE THIRD STAGE OF THE THREE GORGES PROJECT

The closure of third stage diversion channel for the Three Gorges Dam is characterized by large closure discharge, large drop in water surface, and high gap velocity. 1D and 2D flow mathematical models were used in this paper to simulate the hydraulic conditions during the closure. The variation of discharge in the diversion channel and the drops in water surface shared by the upstream and downstream cofferdams were computed using the ID model, and the detailed hydraulic patterns in the diversion channel were simulated using the 2D model. The computed results indicate that the designed closure scheme for discharge of 9,010 m^3/s was feasible for construction, while the designed closure scheme for discharge of 12,200 m^3/s was inapplicable.     

Integrated spatial sampling modeling of geospatial data

Sampling is to, by efficient selection of samples, acquire the accurate information about the population (the research object) at less cost. Spatial sampling is a kind of sampling toward geospatial objects or features with spatial correlation. The differences between effi-cient sampling and completely universal survey lie in quality, time and cost. Sampling provides a kind of economical, prompt and accurate survey[13]. Efficient spatial sampling can be regarded as the optimization of the sampl…     

A patch hierarchy approach to modeling surface and subsurface hydrology in complex flood‐plain environments

Geomorphology interacts with surface‐ and ground‐water hydrology across multiple spatial scales. Nonetheless, hydrologic and hydrogeologic models are most commonly implemented at a single spatial scale. Using an existing hydrogeologic computer model, we implemented a simple hierarchical approach to modeling surface‐ and ground‐water hydrology in a complex geomorphic setting. We parameterized the model to simulate ground‐ and surface‐water ?ow patterns through a hierarchical, three‐dimensional, quantitative representation of an anabranched montane alluvial ?ood plain (the Nyack Flood Plain, Middle Fork Flathead River, Montana, USA). Comparison of model results to ?eld data showed that the model provided reasonable representations of spatial patterns of aquifer recharge and discharge, temporal patterns of ?ood‐water storage on the ?ood plain, and rates of ground‐water movement from the main river channel into a large lateral spring channel on the ?ood plain, and water table elevation in the alluvial aquifer. These results suggest that a hierarchical approach to modeling ground‐ and surface‐water hydrology can reproduce realistic patterns of surface‐ and ground‐water ?ux on alluvial ?ood plains, and therefore should provide an excellent ‘quantitative laboratory’ for studying complex interactions between geomorphology and hydrology at and across multiple spatial scales. Copyright © 2004 John Wiley & Sons, Ltd.     

大冶铁矿露天转地下开采的离散元数值模拟研究

离散单元方法特别适合于节理比较发育的岩体。对于模拟和分析矿山中因采矿引起的围岩松动和冒落等力学行为是比较有效的。用离散单元法研究了大冶铁矿狮子山采区露天转地下开采，采用充填法、崩落法采掘对露天边坡和围岩的影响及围岩变形规律，并得出了相应的结论，为设计和方案的选用提供了一些有益的根据。     

某深埋长隧地应力演化及围岩应力位移模拟研究

结合西部地区某深埋长大公路隧道信息化施工 ,对深埋长隧地应力演化及围岩应力位移进行了弹塑性有限元数值模拟研究。研究结果表明 ,隧道轴线现代地应力状况与隧道埋深、地层岩性及构造发育程度有关 ,最大地应力为40 0MPa左右 ;隧道周边围岩应力在曲边墙底部最大 ,约 3 8 0MPa ,隧道开挖引起的围岩应力影响范围约 2 5 0m ;隧道水平收敛和拱顶下沉位移与隧道埋深近于成直线关系。这些研究结果为深埋长隧信息化设计和施工以及围岩稳定性分析提供了科学依据。     

储层建模中几种原型模型的建立

储层地质知识库的建立是储层建模过程中极为重要的一步，它为建模提供基本的条件数据及各种统计参数。本文总结了建立地质知识库的主要步骤及基本内容。在资料较少的情况下，需要借助原型模型来提供一些无法直接获得但对建模十分重要的参数，如砂体宽厚比、水平方向变差函数等。文中用3个实例较详细的介绍了建立原型模型的3种主要方法，分别是利用露头、成熟油田和沉积模拟实验建立原型模型，并比较了各自的优缺点。     

广东花岗岩类温泉水-汽转化过程中氡浓度的变化及其危害

广东省是我国温泉出露最多的省份之一。在温泉附近建立温泉宾馆,是利用温泉资源大力发展本地旅游业的有效途径之一。本文利用新型NR 667A连续测氡仪测量了广东省4处温泉宾馆室内外空气中氡浓度以及室内温泉水氡浓度。温泉水氡浓度变化范围为53 4～292 5Bq L,室外氡浓度范围是17～48Bq m3,而客房内平均氡浓度水平比没有用温泉水时高出10 9%～813%,浴室内氡浓度比没有用温泉水洗浴时高出13 8%～516%。在宾馆室内使用温泉水有较高的氡转换系数。某些宾馆室内氡浓度超标,可能对宾馆工作人员造成潜在健康危害,应采取一些管理措施或技术手段降低其暴露时间或程度。     

基于共反射点叠加的观测参数选择

在地震勘探数据处理中，通常采用的共中心点(CMP)叠加技术对于复杂的地质构造难以取得好的效果。若通过从道间距、覆盖次数等观测参数与目标层地层倾角、共反射点离散度的关系出发，并结合模型正演，对复杂构造地层的观测参数选取进行优化，使之能对目标层尤其是逆断层下降盘的阴影层段进行较充分的“照明”，即可得到共反射点叠加剖面，这对共反射点叠加是一项有益的研究。     

Seismic horizontal pullout capacity of shallow vertical anchors in sand

The concept of the equivalent free surface has been extended to determine the seismic horizontal pullout capacity of shallow vertical strip plate anchors buried in sand. The analysis has been done rigorously by using the method of stress characteristics. The results have been expressed in the form of non-dimensional charts. The pullout resistance has been found to reduce quite extensively with increase in the magnitude of horizontal earthquake acceleration. The results were compared with the previously published data, and it was seen that the computed pullout resistance with the proposed method was found to be lowest. This revised version was published online in July 2006 with corrections to the Cover Date.