数值模拟在沂水铁矿区地下水环境影响评价中的应用

地下水数值模拟在矿区地下水环境评价中应用较少,通过对山东沂水铁矿矿区进行水文地质详细勘察,建立了尾矿库及矿区地下水水流及水质数学模型,预测和评价了该矿建设实施过程中对地下水环境可能造成的直接影响和间接危害。结论为矿山建设阶段和运行阶段正常情况下,尾矿区水位水质变化对地下水环境没有明显的影响,非正常情况或者事故状态下,预测污染因子在泄漏点附近一定范围出现超标现象;采矿区矿井排水对当地浅层地下水位影响较大。针对这种影响和危害提出防治措施,为沂水铁矿建设项目选址决策、工程设计和环境管理提供科学依据。     

西沙湾工程波浪推算及波浪场数值模拟

本研究根据崇武水文站1976—2007年共31 a的波浪观测资料,使用皮尔逊Ⅲ型频率曲线进行拟合分析,推算出崇武在不同重现期条件下的特征波浪要素,再将其作为外海波浪要素输入,通过基于椭圆型缓坡方程的CGWAVE近岸波浪数值模型,模拟西沙湾的重现期波浪场,得到其设计波浪要素,得到较为合理的计算结果。西沙湾的海底地形强烈地影响近岸波浪,导致了波高和波向的不均匀分布特征,岛屿和礁石都是天然屏障,东侧水域受龟屿及闽台码头的掩护,波浪反射、绕射明显,对避风坞起到了有效的保护作用。此外,对于不同地质的边界,反射系数不同使得近岸波高分布有较大差异。     

三亚湾海岸动力演变的数值模拟研究

近年来,海岸侵蚀呈现不断加剧的趋势,给沿海地区经济社会的发展造成了极大危害。关于三亚湾海岸侵蚀过程及其动力学研究多停留在定性或半定量阶段。从海岸动力地貌学角度出发,基于沉积物迁移理论,利用数值模拟方法,定量研究了三亚湾海岸在波浪和潮流作用下的动力演变过程。计算结果重现了海水波浪的传播形态及动态特征:涨潮时,波高极值相对较大,水位随之升高,波动现象明显,波动形态复杂;落潮时,波高极值相对减小,水位也随之降低,海面相对较平静,波形态在离岸区域较单一,仅在近岸区域相对多样化。同时,数值模拟结果也很好地显示了三亚湾岸线附近的沉积或侵蚀过程特征:三亚湾西段海洋动力条件较强,沉积或侵蚀现象明显,沉积/侵蚀能力由三亚湾西段逐渐向东段减小。此外,以海洋动力学机制为基础,定量模拟了三亚湾海岸侵蚀与演变物理过程,为该地区海岸开发利用、海岸环境保护提供科学依据。     

A semi-analytical solution for frost heave prediction of clay soil

Frost heave is one of the main freezing problems for construction in permafrost regions.The Konrad-Morgenstern segregation potential（SP） model is being used in practice for frost heave using numerical techniques.However,the heat release from in-situ and migrated water in the freezing zone could result in some numerical instability,so the simulation of frost fringe is not ideal.In this study,a semi-analytical solution is developed for frost heave prediction of clay soil.The prediction results to the two tests with different freezing mode with clay soil agree well with the tested behavior,which indicates the feasibility of the solution.     

下地壳流变与造山带同挤压期地壳伸展的动力学关系

应用不可压缩非牛顿粘性流体的本构关系,对造山带同挤压期下地壳流变及其与上地壳构造伸展的动力学关系进行了二维有限元数值模拟. 结果表明,在板块侧向挤压下,当造山带山根下陷和地表隆起达一定程度后,地壳不同层圈岩石将发生复杂的粘性流变. 流变的运动学方式和分布范围不仅与时间有关、同时还受地壳厚度转变带形态的制约. 在构造挤压和山体荷载达到弹性平衡状态后,地壳流变首先发生在造山带下地壳山根,但经一定的Maxwell时间后,流变将不断局限于造山带前缘的厚度转变带. 这一流变方式的变化是导致造山带浅部地壳动力学转变的主要原因. 它造成造山带内上地壳最小主应力从近水平挤压不断转化为近水平拉张,由此使造山带前陆发生挤压冲断的同时,山体的核部发生上地壳的拉张伸展. 最后,应用这一结果讨论了青藏高原南缘南北向地壳伸展的动力学性质.     

考虑风能输入的抛物型缓坡方程

在Radder和Kirby发展的波浪折射绕射缓坡方程抛物型模型基础上,对这种模型进行了改进,改进后的模型除可以考虑波浪传播过程中的底摩阻损耗、非线性作用外,加入了风能输入对波浪传播的影响。基于风能输入项的波浪模型数值计算结果表明,在纯风浪情况下的计算结果与传统的风浪计算方法结果一致,在波浪传播过程中由于风的作用,将导致波高比无风作用下计算的波高大。     

胶州湾大沽河口洪水期三维水沙数值模拟研究

大沽河是胶州湾最大的入海河流,每年都会给胶州湾带来大量泥沙。本文利用三维的ECOMSED水动力-泥沙数值模型对洪水期的大沽河进行了数值模拟,得到了洪水存在期间以及洪水结束后胶州湾的余流分布、大沽河河道上中下游各层盐度、泥沙浓度随时间变化的趋势以及入海冲淡水和泥沙对胶州湾的影响。     

海坛海峡二维潮流场数值模拟

海坛海峡为南北狭长型海峡，海峡内潮波属于前进波．本文建立了平面二维浅水波数学方程，利用欧拉-拉格朗日差分方法得到数值解，模型采用随时间变化的动边界技术，成功地模拟了海坛海峡的前进波特征，并根据实测数据进行了验证．同时计算了同潮时线和等振幅线，不同时刻的潮流场和潮流平均流速分布．计算结果表明，北部湾口M2分潮高潮时间比南部湾口早约5～6min，等振幅线范围约为2．12～2．15m．海峡内流速分布呈南北强、中间弱的特点，最大流速1m／s左右．     

陆态网VLBI站坐标解算

介绍“中国大陆构造环境监测网络”VLBI站坐标解算的方法和过程.采用1994-2009年共50次观测数据解算台站及基线的年变化率,得到了台站的坐标及基线的年变化,并给出了VLBI站坐标解算的统计分析结果.     

广义凸向量极值问题

提出了Ⅰ型、Ⅱ型广义凸最优化问题,给出了广义ｋｕｈｎ－Ｔｕｃｋｅｒ点和广义ｋｕｈｎ－Ｔｕｃｋｅｒ鞍点的概念,并讨论了广义凸最优化问题是Ⅱ型凸的充要条件。