A computerized geologic mapping system based on logical models of geologic structures

In this article we describe the basic framework of the computerized geologic mapping system cigma. The system, whic is based on a mathematical formulation of geologic concepts, consists of the following six subsystems: (1) input of geologic data set; (2) inference of stratigraphic sequence; (3) construction of logical models of geologic structures; (4) determination of three-dimensional geologic boundary surfaces; (5) construction of three-dimensional solid model of geologic structures; and (6) graphical presentation. Geologic structures are summarized in several tables called logical models of geologic structures. Each model is constructed automatically from input data on structural relations between geologic bodies. The model interprets the data automatically to create data files necessary to determine the shapes of geologic boundaries; it also provides a threedimensional solid model of geologic structures referring to the shapes of boundaries. As a prototype, we introduce two types of contacts corresponding to conformity and unconformity into the logical model and show that it is possible to draw a geologic map automatically. More complex geologic structures can be introduced into the geologic mapping system through further formulation of geologic structures.     

Reseach on three-dimensional tide current mathematical model of coast and its application

Ｒｅｓｅａｃｈ　ｏｎ　ｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌ　ｔｉｄｅ　ｃｕｒｒｅｎｔ　ｍａｔｈｅｍｔｉｃａｌ　ｍｏｄｅｌ　ｏｆ　ｃｏａｓｔ　ａｎｄ　ｉｔｓ　ａｐｐｌｉｃａｔｉｏｎTXＲｅｓｅａｃｈｏｎｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌｔｉｄｅｃｕｒｒｅｎ...     

Cohesive sediment transport in the 3D-hydrodynamic-baroclinic circulation model,: study case for idealized tidal inlet

This work provides a general hydrodynamic circulation model that can be used to understand density driven flows, which may arise in the case of suspension of fine-grained materials. The research is expected to provide a better understanding of the characteristics of spatial and temporal variability of current, which is associated with the period of ebb and flood tidal cycles.The model development includes extending the existing three-dimensional (3D) ADCIRC model with (1) baroclinic forcing term and (2) transport module of suspended and soluble materials. The transport module covers the erosion, material suspension and deposition processes for cohesive type sediment. In the case of an idealized tidal inlet in stratified water, the inclusion of baroclinic term can demonstrate the prevailing longshore sediment transport. It is shown that the model has application to the transport of the cohesive sediments from the mouth of the Mississippi River along the north shore of the Gulf of Mexico towards and along the Texas coast.     

Three-dimensional hydrodynamic model, study cases for quarter annular and idealized ship channel problems

An extended version of the three-dimensional hydrodynamic model, ADCIRC 3D-DSS, was utilized to simulate both horizontal and vertical flows in a (quarter) annular harbor (QATP and ATP) and rectangular basin with an idealized ship channel (RBSC). Comparison of horizontal and vertical solutions to the analytical solution and results of other researchers are in good agreement. The vertical velocity solution is highly sensitive to the horizontal velocity solutions. The presence of the sidewall boundary may also affect the vertical solutions. Around the sloping bank of RBSC channel with one-third gradient, the vertical velocity becomes important. The maximum vertical velocity approaches ±50% of the sediment fall velocity of fine sand.     

一种二维和三维嵌套海洋流体动力学数值模式及其在北部湾潮汐和潮流数值模拟中的应用

研究设计了一种二维和三维嵌套、外模态和内模态分离的水动力学数值模式,既可用于潮汐,也可用于风暴潮的数值计算。该模型对全部海区进行二维计算,对其中重点关心的海区同时进行三维计算。在三维计算区域,采用了内、外模态既分离又耦合的计算技术。数值格式采用全部交错的网格结构,三维模型中垂直方向采用σ-坐标代替通常的z-坐标,垂直涡动粘性系数由混合长度理论确定,垂直粘性项采用隐式差分格式。作为算例,本文对南海北部湾潮波进行了细网格的数值计算。     

A three-dimensional modeling of multidirectional random wave diffraction by rectangular submarine pits

A three-dimensional modeling of multidirectional random-wave diffraction by a group of rectangular submarine pits is presented in this paper. The fluid domain is divided into N interior regions representing the pit area and an overall exterior region separated by the imaginary pit boundaries. In the interior region, the analytical expressions of the Fourier series expansion for velocity potentials in the pit regions have been derived with the unknown coefficients determined from a series of Green's function based boundary integral equations. The boundary integral approach has also been applied to obtain the velocity potential and free-surface elevation in the exterior region. The Pierson–Moskowitz (P–M) frequency spectrum was selected for the random wave simulation using the superposition of solutions of a finite number of decomposed wave components. Numerical results for the cases of regular waves and random waves are presented to examine the influences of the pit geometry and incident wave condition on the overall wave field. The general diffraction pattern of alternate bands of increase and decrease of relative wave height in front of the pit system can be observed. It is found that, in the shadow region, the relative wave height is reduced. As the number of pit increases, the effectiveness of reducing the relative wave height behind the multiple-pit system increases. However, the relative wave height within the pit area and in front of the leading pit shows increasing trend. It is noticed that under the random-wave condition, the level of increase and decrease of the relative wave height due to the existence of submarine pits is less pronounced than that observed from results in regular-wave condition.     

Reducing horizontal diffusion errors in σ-coordinate coastal ocean models with a second-order Lagrangian-interpolation finite-difference scheme

When a steep bottom slope exists, it is well known that conventional methods for calculating horizontal diffusion in sigma-coordinate coastal ocean models causes spurious transport (e.g. salinity, temperature, and sediments) and currents. In this study, a second-order accurate finite-difference algorithm and program have been developed to reduce the spurious numerical diffusion errors. In the proposed algorithm, the finite differencing is performed in the x–z coordinate system to approximate the horizontal gradient. Each variable in the finite differential formation is calculated in the sigma-coordinate grid cells using a second-order Lagrangian interpolation polynomial. In conjunction with a stepwise bottom boundary condition, numerical experiments show that the proposed finite-difference scheme considerably reduces numerical errors compared to conventional approaches when dealing with horizontal diffusion over steep topography, which often occurs in coastal oceans and navigation channels.     

三维有云大气辐射变温率数值模拟研究——邻云辐射效应

大气数值模式普遍采用一维辐射传输模型,无法表征有云大气中的三维辐射传输过程,从而影响数值模拟乃至数值天气预报的准确性。为评估这种不确定性,研究了水平云体间的三维辐射相互作用及其改变云体热力结构的规律,力求为改进数值模式辐射计算方案提供理论依据。选取I3RC PhaseⅡ的典型积云场和层积云场作为试验对象,将云场中心区域云体视为目标云体,周围云体为邻近云体,采用宽带三维辐射传输模式SHDOM模拟长波和短波辐射变温率的空间分布,定量阐明邻近云体对目标云体热力结构的影响。结果表明,邻云在长波区域对目标云体主要起辐射保温的作用,目标云体增温区域集中于邻云一侧的云体表面层,增温强度与云覆盖率和云间距离倒数成正比,最高可达3.08K/h,云体增温的厚度与目标云体液态水含量成反比;在短波区域,邻云同时起散射增温和遮蔽降温的作用,太阳垂直入射时,散射增温效应较弱,变温率空间差异小;当太阳天顶角增大后,遮蔽降温效应逐步起主导作用,造成目标云体被邻云遮挡一侧的表面层明显降温,峰值可达-1.72K/h,数值上甚至超过邻云长波增温效应。总之,邻近云体可以明显改变目标云体的变温率空间分布,引入三维邻云辐射效应对改进大气数值模式辐射计算方案具有重要意义。     

A symmetrisation method for non-associated unified hardening model

This paper presents a simple method for symmetrising the asymmetric elastoplastic matrix arising from non-associated flow rules. The symmetrisation is based on mathematical transformation and does not alter the incremental stress–strain relationship. The resulting stress increment is identical to that obtained using the original asymmetrized elastoplastic matrix. The symmetrisation method is applied to integrate the Unified Hardening (UH) model where the elastoplastic matrix is asymmetric due to stress transformation. The performance of the method is verified through finite element analysis (FEA) of boundary value problems such as triaxial extension tests and bearing capacity of foundations. It is found that the symmetrisation method can improve the convergence of the FEA and reduce computational time significantly for non-associated elastoplastic models.     

Soil nailed slope by strength reduction and limit equilibrium methods

A detailed study of nailed slopes under different conditions is reported in this paper. No major difference is found in terms of safety and slip surfaces between the strength reduction method (SRM) and the limit equilibrium method (LEM) in general cases. Appreciable differences between the SRM and LEM appear, however, if the nail load is controlled by the overburden stress. Some special slip surfaces from the SRM obtained by using a very fine mesh are discussed. Field tests demonstrated that the nail head is important in determining the failure mode and the factor of safety of a nailed slope, while the effect of the nail elastic modulus is more noticeable only when the slope is very steep. The optimum layout of the soil nail was found to be longer at bottom and shorter at the top, which is contrary to some engineers’ guidelines for soil nail design during top-down construction. The distribution of tensional force along the soil nail is influenced by the state of the slope (service state, limit state) and the failure modes (external failure, internal failure). In general, the line of maximum tension may not correspond to the critical slip surface as commonly believed, except for the case where the failure mode is an internal tensile failure.