3D GIS下地下水流有限元数值模拟方法与应用

目前缺乏针对较强空间特性水文地质计算参数的高效采集方法与管理机制,而且现有孔隙地下水流有限元数值模拟在三维可视化与空间分析等方面仍存在很多不足。因此,本文以有限元数值计算方法和三维地理信息平台为基础,结合GIS空间分析算法和计算机图形学理论,针对孔隙承压地下水有限元分析过程的关键步骤（概念模型构建、空间离散、水文地质参数提取与赋值等）,提出了孔隙地下水有限元数值模拟过程在3D GIS下的实现方法和技术框架。基于水文地质钻孔数据和GTP（Generalized Tri-Prism）空间数据模型,以盐城市滨海平原水文地质区为研究区域,实现了地下水系统空间结构特征的三维表达,同时与地下水流有限元数值模型相耦合,提出了3D GIS模型与地下水非稳定流有限元数值模型耦合的核心处理算法,并在此基础上开发了3D GIS地下水流数值模拟系统。应用实例表明,该技术框架可有效简化有限元分析流程、优化模型计算效率、实现地下水有限元数值模拟过程及计算结果的三维可视化。     

基于三维粒子系统的雪灾模拟与实时绘制方法

利用遥感影像反演积雪的范围与雪深是当前的主流方法之一,但是由于卫星影像获取的时效性约束,存在难以表示积雪动态信息、模型通用性弱、对遥感资料质量依赖程度高等问题。三维粒子系统作为一种优秀的不规则模糊物体模拟方法,可以大大提高自然现象三维场景的真实感。本文提出了一种基于三维例子系统的雪灾模拟与实时绘制方法,将降雪的物理属性及几何属性转化为雪粒子属性,并以物理引擎为基础,建立基于三维粒子系统的积雪模型和融雪模型,模拟降雪过程并计算积雪范围和积雪深度,真实再现雪灾场景和发展趋势,方法简单、可操作且约束较少。实验表明,该方法的模拟结果与实际灾情具有较强的一致性,算法具有可靠性,并能再现三维雪灾场景,具有实时性和逼真性。     

泥石流数值模拟方法研究进展

泥石流作为一种多相混合介质,所包含的物理过程和动力学特征非常复杂。针对泥石流问题的数值模拟方法随着数值计算方法和物理计算模型的发展而发展,基于物理过程的数值模拟方法为探究泥石流复杂物理现象背后的机理提供了一种有效手段。回顾了求解泥石流动力问题的数值模拟方法,从连续介质计算方法、离散介质计算方法和混合介质计算方法3个方面分析了不同数值模拟方法的特点和适用情况,介绍了在泥石流分析中常用的数值模拟软件及其特点,展望了求解泥石流动力学问题的数值模拟方法的发展趋势。结果表明:传统的基于网格的计算方法已有长足发展和较长的应用历史,相对比较成熟,但是在处理大变形、快速运移的自由表面流问题时,存在网格容易畸变等问题;基于粒子的计算方法在处理上述问题时无需网格的划分和维护,易于确定自由表面位置和多相间的界面,但存在边界条件处理困难等问题;混合介质计算方法在较小尺度范围内对固体颗粒物质与液相相互作用机理进行探讨时具有重要作用。     

十二烷基羧基甜菜碱气液界面自组装行为的分子动力学模拟

应用分子动力学模拟方法研究十二烷基羧基甜菜碱(BS-12)在气/液界面的聚集行为,采用界面密度、径向分布函数、基团占有面积等参数分析BS-12气/液界面的聚集行为。结果表明:BS-12分子在气/液界面紧密排列形成紧密单分子层界面膜,界面膜厚度随着BS-12分子数增加而增大;BS-12分子亲水基团—N+—CH2—COO-中季氮阳离子与水分子相互作用高于亲水基团羧基与水分子间相互作用,加入Na+和Ca2+后羧基与水分子间的径向分布函数值没有变化,也未对扩散双电层产生较强的压缩作用,BS-12分子具有一定的抗盐能力,羧基阴离子与Ca2+间的相互作用高于与Na+间的相互作用,说明Ca2+影响BS-12分子的抗盐性能。     

有限元法与伪谱法混合求解弹性波动方程

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.     

神农架大九湖PAHs多介质归趋模拟

选择美国环境保护署EPA优控污染物NaP（萘）、Phe（菲）、Fla（荧蒽）、BaA（苯并蒽）、Chr（屈）、Pyr（芘）和BaP（苯并芘）7种多环芳烃（polycyclic aromatic hydrocarbons,简称PAHs）作为研究对象,采用Level Ⅲ逸度模型对神农架大九湖地区PAHs多介质归趋进行了模拟,并利用神龙架大九湖大气、土壤、水体和沉积物中PAHs的实际监测数据对模型进行了验证;同时,将气候因素（温度、降水）对环境介质中PAHs转化的影响整合进模型中,探讨了气候变化对PAHs归趋的影响。结果表明:模型计算值与研究区PAHs监测值吻合良好,说明模型参数的选取合理可靠;不同PAHs的相间迁移通量具有一定差异,2环和3环PAHs则以沉积物向水体、水体向大气迁移为主,表现出由沉积物向水体和大气释放的特点,而4环和5环PAHs以大气向土壤和水体向沉积物迁移为主,体现出由大气向土壤和沉积物沉降的特点;温度与水-气、土-气和水-沉积物界面交换净通量存在明显的正相关关系,降水量与水-气、土-气界面交换净通量呈负相关关系,但与水-沉积物界面交换净通量呈正相关关系。本研究结果表明应重视因气候变化而引起的土壤、水体、沉积物的内源释放所造成的二次污染。      

CALPUFF模型模拟结果的时空多维可视化表达

大气污染扩散模型数值模拟结果的时空多维可视化表达,能有效展示大气污染扩散空间分布,为快速高效辅助后期污染分析与应急决策提供重要的参考。CALPUFF模型能高效模拟污染物在大气环境中时空变化的扩散过程,而现有大气污染扩散模型模拟结果可视化表达普遍存在不能有效集成地理信息、缺失部分重要维度特征等问题。论文首先分析了CALPUFF模型,及其与组件式GIS紧密集成的技术方案,为模型模拟结果与地理信息集成提供了技术支持。然后,基于CALPUFF模型的数值计算模拟结果,详细研究了风场数据预处理及其矢量绘制方法,浓度场的数据预处理及其栅格渲染方法。在三维地形场景下,探究了大气污染物的时空扩散模拟可视化表达方法,给出了CALPUFF模型模拟结果与三维地形场景的集成表达策略。最后,通过组件式GIS,构建了福州市三维地形场景,以大气污染CALPUFF扩散模拟结果为例,在WPF系统框架下,实现了风场和浓度场的时空动态模拟和空间层级多维可视化表达。实验结果表明,本文所采用的时空可视化技术方案能将CALPUFF模型模拟结果,有效集成三维地形场景,实现模型模拟结果在时间维度上动态扩散和在空间维度上的层级分布,可为大气污染模拟结果的时空多维可视化提供参考。     

浅谈VisualModflow和FEFLOW

地下水数值模拟是各种数值方法在水文地质计算中的应用,是目前定量研究地下水水资源量的重要手段。VisualModflow和FEFLOW都是当前世界上十分流行的可用于模拟三维地下水流和溶质运移模拟评价专业软件系统,属商业软件。它们都具有直观的、强有力的图形交互界面,模型剖分、输人参数和模拟结果,都可以用图形显示,并支持三维可视化,做到了真正的人机对话,在许多行业和部门内得到了广泛的应用。但同时,它们之间却存在着许多差别,各有千秋：     

生长模型驱动的单株杉木三维动态模拟

虚拟植物是一种潜在、有力的植物分析辅助工具,而单株木生长模拟是森林生态系统模拟的基础。为了动态模拟杉木生长发育过程,提出了参数化单株木三维形态结构建模和与距离无关的单木生长模型的集成方法,其主要过程为：首先,根据杉木的形态结构特征,采用参数化建模方法建立特定生长阶段的三维静态模型;其次,根据生长模型预测不同年龄杉木的树高、胸径、枝下高和分枝轮数,并与树木三维静态模型几何描述参数建立联动;最后,采用参数曲线调整干径和枝径变化、枝条长度、分枝角等形态结构参数,使模型形态随树龄增长而变化。在自主研发的ParaTree系统上,扩展杉木动态生长模拟模块,并以福建省漳平五一林场的二类调查数据为例,动态模拟了杉木的生长过程。模拟结果表明,本方法可直观表达林分中林木个体平均生长状况和大体形态结构特征。树木三维模型描述参数与传统树木统计生长模型结合,有利于重用林业领域淀积大量的生长模型。     

基于Visual ModFlow的苑庄水源地地下水位预测

根据苑庄水源地的水文地质条件,建立了水源地的水文地质概念模型。在此基础上,采用Visual ModFlow对水源地地下水流场进行了数值模拟。误差统计分析表明,数值计算模型在识别与验证过程中的误差较小,模拟效果较为理想。根据所建立的数值模型,并按照设计开采量,对水源地运行后区域地下水位的变化趋势进行了预测。     

A quasi single-phase model for debris flows and its comparison with a two-phase model

A depth-averaged quasi single-phase mixture model is proposed for debris flows over inclined bed slopes based on the shallow water hydrosediment-morphodynamic theory with multi grain sizes. The stresses due to fluctuations are incorporated based on analogy to turbulent flows, as estimated using the depth-averaged k-? turbulence model and a modification component. A fully conservative numerical algorithm, using wellbalanced slope limited centred scheme, is deployed to solve the governing equations. The present quasi single-phase model using four closure relationships for the bed shear stresses is evaluated against USGS experimental debris flow and compared with traditional quasi single-phase models and a recent physically enhanced two-phase model. It is found that the present quasi single-phase model performs much better than the traditional models, and is attractive in terms of computational cost while the two-phase model performs even better appreciably.     

Simulation of two-phase debris flow scouring bridge pier

Debris flows are typical two-phase flows, which commonly accompany erosion in mountainous areas, and may destroy bridge engineering by scouring. In this study, a physically-based two-phase model is applied for the simulation of debris flow scouring of bridge pier. In this model, the shear stress of debris flow on an erodible bed is considered to be a function of the solid shear stress, fluid shear stress, and volume fraction; accordingly, the erosion is incorporated into the two-phase model. Using a highaccuracy computational scheme based on the finite volume method, the model is employed for simulating a dynamic debris flow over an erodible bed. The numerical results are consistent with the experimental data, and verify the feasibility of the two-phase model. Moreover, a simple numerical test is performed to exhibit the fundamental behaviour of debris flow scouring of bridge pier, which shows that the degree of erosion on each side of the pier is higher compared to other areas. The scouring depth is influenced by the variations of solid volume fraction and velocity of debris flow and pier width.     

视频监视系统中运动目标检测的新方法

运动目标检测是计算机视觉监视系统的核心。对于采用固定摄像机监控视频运动目标的检测,利用Affine Lucas Kanade特征跟踪算法的图像金字塔模型,根据模型的不同层数,将对应不同层数分辨率下的差分图像与Susan算子结合,通过形态学处理实现对复杂背景下运动目标的检测。仿真实验结果表明,采用的算法有效地抑制了噪声对运动目标检测效果的影响,且计算量随金字塔层数的增加而成倍减少。     

Dynamic process simulation with a Savage-Hutter type model for the intrusion of landslide into river

Natural damming of rivers by mass movements is a very common and potentially dangerous phenomena which has been documented all over the world. In this paper, a two-layer model of Savage-Hutter type is presented to simulate the dynamic procedure for the intrusion of landslide into rivers. The two-layer shallow water system is derived by depth averaging the incompressible Navier-Stokes equations with the hydrostatic assumption. A high order accuracy scheme based on the finite volume method is proposed to solve the presented model equations. Several numerical tests are performed to verify the realiability and feasibility of the proposed model. The numerical results indicate that the proposed method can be competent for simulating the dynamic process of landslide intrusion into the river. The interaction effect between both layers has a significant impact on the landslide movement, water fluctuation and wave propagation.     

Interaction of Waves, Surface Currents, and Turbulence： the Application of Surface-Following Coordinate Systems

Surface waves comprise an important aspect of the interaction between the atmosphere and the ocean, so a dynamically consistent framework for modelling atmosphere-ocean interaction must take account of surface waves, either implicitly or explicitly. In order to calculate the effect of wind forcing on waves and currents, and vice versa, it is necessary to employ a consistent formula- tion of the energy and momentum balance within the airflow, wave field, and water column. It is very advantageous to apply sur- face-following coordinate systems, whereby the steep gradients in mean flow properties near the air-water interface in the cross-interface direction may be resolved over distances which are much smaller than the height of the waves themselves. We may account for the waves explicitly by employing a numerical spectral wave model, and applying a suitable theory of wave–mean flow interaction. If the mean flow is small compared with the wave phase speed, perturbation expansions of the hydrodynamic equations in a Lagrangian or generalized Lagrangian mean framework are useful: for stronger flows, such as for wind blowing over waves, the presence of critical levels where the mean flow velocity is equal to the wave phase speed necessitates the application of more general types of surface-following coordinate system. The interaction of the flow of air and water and associated differences in temperature and the concentration of various substances (such as gas species) gives rise to a complex boundary-layer structure at a wide range of vertical scales, from the sub-millimetre scales of gaseous diffusion, to several tens of metres for the turbulent Ekman layer. The bal- ance of momentum, heat, and mass is also affected significantly by breaking waves, which act to increase the effective area of the surface for mass transfer, and increase turbulent diffusive fluxes via the conversion of wave energy to turbulent kinetic energy.     

Deposition Morphology of Non-homogeneous Debris Flow and Its Energy Characteristics

Non-homogeneous two-phase debris flows are widely found in the western mountainous regions of China. To investigate the characteristics of the debris flow deposition process related to the morphology and extent of the debris fan, a series of physical experiments were carried out using an experimental flume. Some useful relationships were obtained to link the flow velocity with the geometric characteristics of deposition morphology and the corresponding area or volume. Based on these, some expressions about energy dissipation process in both the transport-deposition zone and deposition zone are presented, and improved equations describing solidliquid two-phase energy transformations in the specific deposition zone are also established. These results provide a basis for further investigating the underlying mechanisms of non-homogeneous debris flows, based upon which effective disaster control measures can be undertaken.     

Modeling of low-frequency seismic waves in a shallow sea using the staggered grid difference method

Elastic waves in the seabed generated by low-frequency noise radiating from ships are known as ship seismic waves and can be used to detect and identify ships.To obtain the propagation characteristics of ship seismic waves,an algorithm for calculating seismic waves at the seafl oor is presented based on the staggered-grid fi nite diff erence method.The accuracy of the algorithm was tested by comparison with analytical solutions.Numerical simulation of seismic waves generated by a low-frequency point sound source in a typical shallow sea environment was carried out.Using various source frequencies and locations in the numerical simulation,we show that the seismic waves in the near fi eld are composed mostly of transmitted S-waves and interface waves while transmitted P-waves are weak near the seafl oor.However,in the far fi eld,the wave components of the seismic wave are mainly normal modes and interface waves,with the latter being relatively strong in the waveforms.As the source frequency decreases,the normal modes become smaller and the interface waves dominate the time series of the seismic waves.     

New Damage-locating Method for Bridges Subjected to a Moving Load

A new damage-locating method for bridges subjected to a moving load is presented, and a new ‘moving load damage-locating indicator＇ （MLDI） is introduced. A vehicle is modeled as a moving load, the bridge is simplified as an Euler-Bernoulli beam, and the damage is simulated by a reduction of stiffness properties of the elements. The curvature and MLDI values at each node of the baseline model （undamaged） and the damage model are computed respectively. Then the damage or damages can be located from a sudden change of the MLDI value. The feasibility and effectiveness of the proposed method are validated by numerical simulation. The results indicate that the method is effective, being able to not only locate a single damage accurately, but also locate multiple damages in simply-supported bridges, including multiple damages in continuous bridges. The results also indicate that the MLDI can accurately locate damages under 5% measurement noise.     

Debris flows introduced in landslide deposits under rainfall conditions: The case of Wenjiagou gully

Debris flows often occur in landslide deposits during heavy rainstorms. Debris flows are initiated by surface water runoff and unsaturated seepage under rainfall conditions. A physical model based on an infinitely long, uniform and void-rich sediment layer was applied to analyze the triggering of debris-flow introduced in landslide deposits. To determine the initiation condition for rainfall-induced debris flows, we conducted a surface water runoff and saturated-unsaturated seepage numerical program to model rainfall infiltration and runoff on a slope. This program was combined with physical modeling and stability analysis to make certain the initiation condition for rainfall-introduced debris flows. Taking the landslide deposits at Wenjiagou gully as an example, the initiation conditions for debris flow were computed. The results show that increase height of surface-water runoff and the decrease of saturated sediment shear strength of are the main reasons for triggering debris-flows under heavy rainfall conditions. The debris-flow triggering is affected by the depth of surface-water runoff, the slope saturation and shear strength of the sediment.     

Comparative Analysis of Seismic Response Characteristics of Pile-Soil-Stnicture Interaction System

The study on the earthquake-resistant performance of a pile-soil-structure interaction system is a relatively complicated and primarily important issue in civil engineering practice. In this paper, a computational model and computation procedures for pile-supported structures, which can duly consider the pile-soil interaction effect, are established by the finite element method. Numerical implementation is made in the time domain. A simplified approximation for the seismic response analysis of pile-soil-structure systems is briefly presented. Then a comparative study is performed for an engineering example with numerical results computed respectively by the finite element method and the simplified method. Through comparative analysis, it is shown that the results obtained by the simplified method well agree with those achieved by the finite element method. The numerical results and findings will offer instructive guidelines for earthquake-resistant analysis and design of pile-supported structures.