江西石城国家地质公园龟裂地貌特征及成因分析

江西石城国家地质公园的龟裂地貌主要产于晚白垩世河口组红色砂岩中,造型奇特、多样,形成许多象形景观,如龟寿石、松果石、鳄鱼石、仙人犁田等,引人注目。但是,对这些龟裂地貌的微观研究不足,导致对其成因的解释比较模糊。本次在野外用罗盘测量裂隙走向,同时采集砂岩样品进行偏光显微镜和扫描电子显微镜实验观察,分析龟裂地貌形成的原因。野外观察发现,龟裂地貌主要发育在颗粒相对均匀的厚层细粒砂岩中,可产在砂岩层面或其它暴露面上,向下延伸一般不超过15 cm。在偏光显微镜下,砂岩样品的碎屑颗粒主要为次棱角状-次圆状,分选性较好,填隙物中含有方解石胶结物。在扫描电镜下,颗粒胶结疏松,颗粒表面的溶蚀小孔和裂隙比较常见。结合国内外类似研究,本文认为岩性是龟裂地貌发育的基本要素,裂隙走向的玫瑰花图显示,断裂构造控制了一部分龟裂的形成,但并不是唯一的原因。在风化作用过程中,由于硅、铁、锰氧化物沉淀而在砂岩表层形成薄层结壳,后由于应力拉张使结壳破裂,形成细小的裂隙,流水沿着裂隙侵蚀使其扩大,最终形成龟裂地貌。最后,提出龟裂地貌的四个阶段演化模式。     

基于网格细分技术的三维地质模型光滑方法研究

针对粗糙的三维地质模型,利用网格细分技术,生成加密光滑的三维模型。分析网格细分技术在三维地质建模中的应用需求以及Loop细分和改进的蝶形细分技术,为保持地质体之间公共面数据的一致性,对"改进的蝶形细分法"做了进一步改进,增加了对边界约束的处理。探讨利用细分技术生成多分辨率模型,通过试验验证了该方法的有效性。     

平面影像到QTM像元的转换算法及精度分析

平面影像数据到四元三角网(Quaternary Triangular Mesh,QTM)像元的转换是实现全球QTM格网影像无缝建模的首要任务。提出基于面积权的数据转换方法,给出不同类型的权重及格网值的计算,并利用ArcGIS中256级灰度的Wsiearth.tif图像数据分析了转换精度。结果表明:将4 km分辨率的影像数据转换到12层的QTM像元时,转换误差在0~2灰度内占96.27%,对一般应用的影响可以忽略不计。     

The effects of resolving the Canadian Arctic Archipelago in a finite element sea ice model

Though narrow straits may have a strong influence on the large-scale sea ice mass balance, they are often crudely represented in coarse resolution sea ice models. Unstructured meshes, with their natural ability to fit boundaries and locally increase the mesh resolution, propose an alternative framework to capture the complex oceanic areas formed by coasts and islands. In this paper, we develop a finite element sea ice model to investigate the sensitivity of the Arctic sea ice cover features to the resolution of the narrow straits constituting the Canadian Arctic Archipelago. The model is a two-level dynamic-thermodynamic sea ice model, including a viscous-plastic rheology. It is run over 1979–2005, forced by daily NCEP/NCAR reanalysis data. Confronting qualitatively numerical experiments with observations shows a good agreement with satellite and buoys measurements. Due to its simple representation of the oceanic interactions, the model overestimates the sea ice extent during winter in the southernmost parts of the Arctic, while the Baffin Bay and Kara Sea remain ice-covered during summer. In order to isolate the benefits from resolving the Canadian Arctic Archipelago, a numerical experiment is performed where we artificially close the archipelago. Focusing on the large-scale sea ice thickness pattern, no significant change is found in our model, except in the close surroundings of the archipelago. However, the local and short-term influences of the ice exchanges are nonnegligible. In particular, we show that the ice volume associated to the Canadian Arctic Archipelago represents 10% of the Northern Hemisphere sea ice volume and that the annual mean ice export towards Baffin Bay amounts to 125 km³ yr⁻¹, which may play an important role on the convective overturning in the Labrador Sea.     

Arbitrary Lagrangian–Eulerian method for large‐strain consolidation problems

In this paper, an arbitrary Lagrangian–Eulerian (ALE) method is generalized to solve consolidation problems involving large deformation. Special issues such as pore‐water pressure convection, permeability and void ratio updates due to rotation and convection, mesh refinement and equilibrium checks are discussed. A simple and effective mesh refinement scheme is presented for the ALE method. The ALE method as well as an updated‐Lagrangian method is then used to solve some classical consolidation problems involving large deformations with different constitutive laws. The results clearly show the advantage and efficiency of the ALE method for these examples. Copyright © 2007 John Wiley & Sons, Ltd.     

Simulating the hydraulic characteristics of the lower Yellow River by the finite‐volume technique

The finite‐volume technique is used to solve the two‐dimensional shallow‐water equations on unstructured mesh consisting of quadrilateral elements. In this paper the algorithm of the finite‐volume method is discussed in detail and particular attention is paid to accurately representing the complex irregular computational domain. The lower Yellow River reach from Huayuankou to Jiahetan is a typical meandering river. The generation of the computational mesh, which is used to simulate the flood, is affected by the distribution of water works in the river channel. The spatial information about the two Yellow River levee, the protecting dykes, and those roads that are obviously higher than the ground, need to be used to generate the computational mesh. As a result these dykes and roads locate the element interfaces of the computational mesh. In the model the finite‐volume method is used to solve the shallow‐wave equations, and the Osher scheme of the empirical function is used to calculate the flux through the interface between the neighbouring elements. The finite‐volume method has the advantage of using computational domain with complex geometry, and the Osher scheme is a method based on characteristic theory and is a monotone upwind numerical scheme with high resolution. The flood event with peak discharge of 15 300 m³/s, occurring in the period from 30 July to 10 August 1982, is simulated. The estimated result indicates that the simulation method is good for routing the flood in a region with complex geometry. Copyright © 2002 John Wiley & Sons, Ltd.     

内陆单斜构造内咸水入侵淡水含水层三维数值模拟：——以山西柳林？…

本文建立了内陆单斜构造内咸水入侵淡水含水层的三维可混溶溶质运移模型。同时,为了尽可能提高模拟精度,还建立了区域二维水流模型和水质模型,为三维模型提供边界条件。模型应用于描述柳林电厂水源地咸－淡水界面运移特征,模拟效果良好,还据此预报了咸水入侵的可能性并探讨了淡水流的强弱对咸－淡水界面运移的影响。     

多边形有限单元形函数的几何构造方法

根据非均质材料的细观结构,划分多边形的计算网格。采用多边形单元进行有限元分析,实现了基于材料真实结构的数值模拟。给出了多边形有限单元形函数的几何构造方法,构造了一个辅助多边形,采用散度定理推导出多边形单元形函数的表达式。分析总结了多边形Wachspress插值、Laplace插值和平均值插值的构造方法和性质。     

Numerical Prediction of Hydrodynamic Forces on A Ship Passing Through A Lock

While passing through a lock, a ship usually undergoes a steady forward motion at low speed. Owing to the size restriction of lock chamber, the shallow water and bank effects on the hydrodynamic forces acting on the ship may be remarkable, which may have an adverse effect on navigation safety. However, the complicated hydrodynamics is not yet fully understood. This paper focuses on the hydrodynamic forces acting on a ship passing through a lock. The unsteady viscous flow and hydrodynamic forces are calculated by applying an unsteady RANS code with a RNG k?ε turbulence model. User-defined function (UDF) is compiled to define the ship motion. Meanwhile, the grid regeneration is dealt with by using the dynamic mesh method and sliding interface technique. Numerical study is carried out for a bulk carrier ship passing through the Pierre Vandamme Lock in Zeebrugge at the model scale. The proposed method is validated by comparing the numerical results with the data of captive model tests. By analyzing the numerical results obtained at different speeds, water depths and eccentricities, the influences of speed, water depth and eccentricity on the hydrodynamic forces are illustrated. The numerical method proposed in this paper can qualitatively predict the ship-lock hydrodynamic interaction. It can provide certain guidance on the manoeuvring and control of ships passing through a lock.     

Numerical prediction of hydrodynamic forces on a ship passing through a lock

While passing through a lock, a ship usually undergoes a steady forward motion at low speed. Owing to the size restriction of lock chamber, the shallow water and bank effects on the hydrodynamic forces acting on the ship may be remarkable, which may have an adverse effect on navigation safety. However, the complicated hydrodynamics is not yet fully understood. This paper focuses on the hydrodynamic forces acting on a ship passing through a lock. The unsteady viscous flow and hydrodynamic forces are calculated by applying an unsteady RANS code with a RNG k-ε turbulence model. User-defined function(UDF) is compiled to define the ship motion. Meanwhile, the grid regeneration is dealt with by using the dynamic mesh method and sliding interface technique. Numerical study is carried out for a bulk carrier ship passing through the Pierre Vandamme Lock in Zeebrugge at the model scale. The proposed method is validated by comparing the numerical results with the data of captive model tests. By analyzing the numerical results obtained at different speeds, water depths and eccentricities, the influences of speed, water depth and eccentricity on the hydrodynamic forces are illustrated. The numerical method proposed in this paper can qualitatively predict the ship-lock hydrodynamic interaction. It can provide certain guidance on the manoeuvring and control of ships passing through a lock.