基于GIS和SD方法的社会经济发展与生态环境响应动态模拟预测研究

在着重于空间现象分析处理的GIS与动态模拟的SD方法结合的模式下，动态模拟预测了不同社会经济发展策略对生态环境的影响，揭示了生态经济系统社会经济发展与生态环境相互影响关系，为区域可持续发展提供决策支持服务。     

水下无源导航系统仿真匹配算法研究

针对水下无源导航(地形辅助导航)系统的特点，充分利用各种海洋地理特征传感器获取信息，进行数据预处理、算法分析、性能评估，然后进行智能数据融合，进而完成水下无源导航系统模块化的仿真设计，为水下无源导航系统的进一步研究提供可行的方法和技术。仿真结果表明，采用水下地形匹配辅助导航，能够有效地修正惯导误差达到设计要求，为真实海洋环境下的精确导航提供可靠的技术依据。     

Determination of silica sand stiffness

Site investigation and evaluation of properties of soil or rock are important aspects of geotechnical design. Determination of the ground stiffness is one of the important parameters in geotechnical engineering. Since the measurement of shear modulus is very sensitive to soil disturbance, especially for sand, determination of the stiffness of soil in the field is more reliable than in laboratory tests on sampled specimens. Measurement of shear modulus is one of the most common applications of self-boring pressuremeter testing. As an in situ device, the pressuremeter provides a unique method for assessing directly the in situ shear modulus of a soil. This paper describes a laboratory study of silica sand stiffness, which varies with stress level and strain amplitude. The results show that the elastic shear modulus value is mainly dependent on the value of the mean effective stress and relative density.     

An experimental investigation of LNAPL migration in an unsaturated/saturated sand

Accidental spills of hydrocarbons, such as Light Non-Aqueous Phase Liquids (LNAPLs), are one of the most common sources of subsurface contamination. Migration of LNAPL in a porous medium is influenced by various factors such as the number of fluids present in the unsaturated/saturated zones and the proportion of pores occupied by each fluid. The results for relationship between matric suction and degree of saturation are presented in this paper for water–air, water–LNAPL and LNAPL–air systems in a sand. A simple and reliable setup using Buchner funnel was designed to obtain these relations. It was found that the relationship between matric suction head and degree of saturation is hysteretic for all the fluid systems (water–air, water–LNAPL and LNAPL–air). Furthermore, the amount of hysteresis depended upon the fluid system, with the maximum hysteresis occurring for water–air system. The results suggest that the amount of trapped air depends upon the reversal degree of saturation from drying to wetting.     

Atomistic simulation of the structure and elastic properties of pyrite (FeS2) as a function of pressure

Interatomic potential parameters have been derived at simulated temperatures of 0 K and 300 K to model pyrite FeS₂. The predicted pyrite structures are within 1% of those determined experimentally, while the calculated bulk modulus is within 7%. The model is also able to simulate the properties of marcasite, even though no data for this phase were included in the fitting procedure. There is almost no difference in results obtained for pyrite using the two potential sets; however, when used to model FeS₂ marcasite, the potential fitted at 0 K performs better. The potentials have also been used to study the high-pressure behaviour of pyrite up to 44 GPa. The calculated equation of state gives good agreement with experiment and shows that the Fe–S bonds shorten more rapidly that the S–S dimer bonds. The behaviour of marcasite at high pressure is found to be similar to that of pyrite.     

Wave-influenced deltas: geomorphological implications for facies reconstruction

ABSTRACT A process‐based facies model for asymmetric wave‐influenced deltas predicts significant river‐borne muds with potentially lower quality reservoir facies in prodelta and downdrift areas, and better quality sand in updrift areas. Many ancient barrier‐lagoon systems and ‘offshore bars’ may be better reinterpreted as components of large‐scale asymmetric wave‐influenced deltaic systems. The proposed model is based on a re‐evaluation of several modern examples. An asymmetry index A is defined as the ratio between the net longshore transport rate at the mouth (in m³ year^?1) and river discharge (in 10⁶ m³ month^?1). Symmetry is favoured in deltas with an index below ≈ 200 (e.g. Tiber, lobes of the Godavari delta, Rosetta lobe of the Nile, Ebro), whereas deltas with a higher index are asymmetric (e.g. Danube – Sf. Gheorghe lobe, Brazos, Damietta lobe of the Nile). Periodic deflection of the river mouth for significant distances in the downdrift direction occurs in extreme cases of littoral drift dominance (e.g. Mahanadi), resulting in a series of randomly distributed, quasi‐parallel series of sand spits and channel fills. Asymmetric deltas show variable proportions of river‐, wave‐ and tide‐dominated facies both among and within their lobes. Bayhead deltas, lagoons and barrier islands form naturally in prograding asymmetric deltas and are not necessarily associated with transgressive systems. This complexity underlines the necessity of interpreting ancient depositional systems in a larger palaeogeographic context.     

Are catch-up reefs an artefact of coring?

Drill cores through modern coral reefs commonly show a time lag in reef initiation followed by a phase of rapid accretion to sea level from submerged foundations – the so-called ‘catch-up response’. But because of the difficulty of drilling in these environments, core distribution is usually restricted to accessible areas that may not fully represent reef history, especially if the reef initiated in patches or developed with a prograde or retrograde geometry. As a consequence, core data have the potential to give a misleading impression of reef development, particularly with respect to the timing of initiation and response to sea-level rise. Here, we use computer models to simulate keep-up reef development and, from them, quantify variations in the timing of reef initiation and accretion rate using mock cores taken through the completed simulations. The results demonstrate that cores consistently underestimate the timing of reef initiation and overestimate the reef accretion rate so that, statistically, a core through a keep-up reef will most likely produce a catch-up pattern – an initiation lag followed by a phase of rapid accretion to sea level. This implies that catch-up signatures may be an artefact of coring and that keep-up reefs are significantly more common than previous core studies claim.     

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.     

新疆某水利枢纽工程软弱结构面抗剪强度参数的试验研究

某水利枢纽工程位于伊犁地区,影响工程的主要问题为坝基岩体软弱结构面和断层的抗滑稳定问题。根据工程实际条件和软弱夹层特性,采取现场断层软弱结构面试验和室内模拟试验相结合,获得软弱夹层的抗剪参数。实践证明,把两种方法结合起来的试验成果可靠,尤其对一些不具备进行现场大型岩体软弱夹层抗剪试验的工程,采用室内中型直剪仪进 行剪切试验具有实用意义。     

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

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