山东烟台夹河中、下游地区海水入侵三维水质数值模拟研究

建立了三维变密度对流弥散水质数学模型来研究山东省烟台夹河中、下游地区咸淡水界面的运移规律。以四面体为基本离散单元 ,推导出三维海水入侵变密度水质模型求解的数值方法 ,其中水流方程求解时运用了迦辽金有限单元法。溶质运移方程求解时运用了欧拉拉格朗日混合方法 ,将对流项与弥散项分离 ,用传统迦辽金有限元方法求解弥散项 ;采用自适应MOC MMOC法求解对流项 ,以消除人工过量和数值弥散。根据地下水的潮汐效应观测信息 ,确定了含水系统的海底延伸边界 ;利用该地区地下水水头及水质长观资料识别了模型的水文地质参数 ,探讨了夹河地区海水入侵的原因 :认为夹河下游地区滨海地带地下水过量开采是造成烟台地区海水入侵的主要原因。此外 ,海水随潮定期地倒灌进入夹河 ,通过局部岩性天窗侵入淡水含水层加剧了沿夹河河床两侧地下水的咸化。同时还预测了几种情况下地下水的水质演化趋势 ,为防止和减轻夹河地区海水入侵提供合理、科学的依据。     

基于动态特征的强对流云团追踪

利用FY-2E卫星数据获取的强对流云团面积、重心、长短轴比、重心与形心距离、移动速度、移动角度和最低亮温等属性的变化可作为动态特征,利用慢特征分析方法提取云团中具有一定连续性和稳定性的动态特征对强对流云团不同阶段进行识别和追踪.结果表明,动态特征与强对流云团的不同发展阶段具有很好的对应关系:在初生阶段,云团的移动方向和速度不稳定,但是面积呈现出缓慢增长态势,云顶亮温缓慢下降,此时云团的慢特征为面积和云顶亮温;在成熟阶段,云团的移动路径趋于稳定,云顶亮温达到最低,云团重心和形心基本重合;在消散阶段,存在云团分裂和云团的重心与形心分离特征.云团长短轴比的变化与云团最低亮温的变化趋势一致,移速缓慢的对流云团更容易造成集中强降水,快速移动的对流云团大多造成地面大风.     

固体矿产资源储量估算应注意的问题

参与资源储量估算的探矿工程和原始资料质量应符合规范要求,并正确引用工业指标,根据地质规律和规范要求合理圈定和外推矿体;估算方法要符合矿床地质特征,能用简单方法的不用复杂方法;按照规范要求合理确定资源储量类型,块段划分要考虑资源储量类型、矿石类型和工程网度等;正确确定估算参数.从而获得可靠的资源储量结果.     

半隐式半拉格朗日动力框架的动能谱分析

大量的观测事实表明自由大气动能谱与波数之间满足如下关系:在大尺度区域满足E∝k-3关系,过渡到中尺度区域表现为E∝k-5/3关系.数值模式动能谱是测量模式动力框架的耗散的直接度量,而耗散对模式的性能有着重要影响,因此动能谱是研究和评估模式动力框架的非常规的有效方法.文中使用基于半隐式半拉格朗日动力框架的全球/区域一体化模式GRAPES进行数值模拟试验,然后计算GRAPES模式的动能谱并与实际观测得到的大气动能谱比较,得到GRAPES模式能够很好地复制出实际大气动能谱的分布特征,包括从大尺度区域的E∝k-3关系向中尺度的E∝k-5/3关系的过渡特征.并且发现GRAPES模式存在最大有效时间步长,当时间步长小于最大有效时间步长时,模式动能谱随时间步长增大而逐渐衰减;当时间步长大于最大有效时间步长时,模式动能谱随时间步长增大而虚假增长.同时通过与实际大气动能谱比较,发现模式动能谱在5Δx波长附近开始明显衰减,因此将5Δx波长定义为GRAPES模式的最高有效分辨尺度;当空间分辨率提高与时间步长等相协调时,中小尺度模式动能谱向中小尺度延伸而更接近实际大气动能谱;当空间分辨率提高与时间步长等不相协调时,中小尺度模式动能谱存在较大误差,相应的大尺度模式动能谱亦存在较大误差.此外,时间步长对模式spin up过程有着重要的影响,较小时间步长时,spin up过程能够很好发展出合理的动能谱结构,在物理空间上表现为模式能够在spin up时间内生成和发展出合理的中小尺度系统;而较大时间步长时,spin up过程很难发展出合理的动能谱结构,在物理空间上表现为模式未能在spin up时间内生成和发展出合理的中小尺度系统.最后,GRAPES模式动能谱与WRF模式动能谱具有一致性,GRAPES全球中期模式能够完美地模拟出大尺度的E∝k-3动能谱特征.综上所述,本文通过研究GRAPES模式动力框架的动能谱得到了一些有意义的结果,为进一步研究、完善、优化和应用模式提供了科学的指导,可见动能谱是评估模式动力框架的有效方法.     

半隐式半拉格朗日平方守恒计算格式的构造

在显式半拉格朗日完全平方守恒格式基础上, 构造出半隐式半拉格朗日完全平方守恒计算格式, 它继承了显式半拉格朗日完全平方守恒格式的优点, 并突破计算不稳定柯朗条件对时间步长的约束, 使时间步长大为增长.此外, 还给出这种新的计算格式在一维原始方程上的应用.     

交叉相关算法在强对流天气临近预报中的应用

利用"雨燕"临近预报系统的雷达回波外推算法,对三次强对流个例的3 km CAPPI反射率因子、回波顶和垂直累积液态含水量进行0～1小时外推预测,将预测结果和组合阈值比较进行强对流天气预警区域预报试验。资料来源于CINRAD-SB雷达基数据经RPG算法处理得到的雷达产品数据。所用的方法是在扩展后的交叉相关追踪算法基础上实现雷达回波移动预测。试验结果表明:3 km CAPPI反射率因子、回波顶和垂直累积液态含水量雷达产品与实况有较好的吻合,能较好地预报出回波的形状、变化趋势和移动方向,回波的范围、位置和强度中心与实况相似,外推时间越短,预报效果越好。交叉相关追踪算法外推得到的雷达回波结果在0～1 h内是可用的;外推预测结果经过组合阈值的过滤后,得出的强对流天气预警区域预警结果也是可行的,对强对流天气临近预报有一定指示意义。     

利用矩阵分解理论的双程波方程叠前深度偏移方法

叠前深度偏移理论及方法一直是地震数据成像中研究的热点问题.业界对单程波叠前深度偏移方法和逆时深度偏移开展了深入的研究,但对双程波方程波场深度延拓理论及成像方法的研究还鲜有报道.本文以地表记录的波场值为基础,利用单程波传播算子估计波场对深度的偏导数,为在深度域求解双程波方程提供充分的边界条件,并提出利用矩阵分解理论实现双程波方程的波场深度外推.通过对强速度变化介质中传播波场的计算,与传统的单程波偏移方法相比,本文提出的偏移方法计算的波场与常规有限差分技术计算的波场相一致,证明了本方法计算的准确性.通过对SEAM模型的成像,在相同的成像参数下,与传统的单程波偏移算法和逆时深度偏移算法方法相比,本文提出的偏移方法能够提供更少的虚假成像和更清晰的成像结果.本文所提偏移算法具有深度偏移和双程波偏移的双重特色,推动和发展了双程波叠前深度偏移的理论和实践.     

Application of the Characteristic CIP Method to a Shallow Water Model on the Sphere

Semi-implicit algorithms are popularly used to deal with the gravitational term in numerical models. In this paper, we adopt the method of characteristics to compute the solutions for gravity waves on a sphere directly using a semi-Lagrangian advection scheme instead of the semi-implicit method in a shallow water model, to avoid expensive matrix inversions. Adoption of the semi-Lagrangian scheme renders the numerical model always stable for any Courant number, and which saves CPU time. To illustrate the efficiency of the characteristic constrained interpolation profile (CIP) method, some numerical results are shown for idealized test cases on a sphere in the Yin-Yang grid system.     

一种基于多项式外推的局部透射边界位移解(外行波为非平面波情形)

一般情况下,通过人工边界向外透射的波动常为非平面波。对于近场非平面波在人工边界上的透射问题,不能采用简单的平面波透射方法。本文针对近场非平面波在人工边界上的透射特点,分析了非平面波沿人工边界法向视波速的变化规律及其近场失稳机制,推导了适合近场非平面波的、人工边界点位移解的不等步距线性外推公式。     

两种北斗卫星精密星历外推方法的精度分析

国际全球卫星导航系统服务(IGS)中心只提供当天00:00:00—23:45:00的精密星历数据,若要获取当天完整的星历数据,往往采用坐标外推的方法得到23:45:00—24:00:00时段的星历数据。本文采用三角函数插值法和广义延拓逼近法对一段时间不同类型和不同时间间隔的北斗卫星精密星历轨道坐标进行外推。结果表明,GEO卫星精密星历轨道坐标外推精度优于IGSO卫星和MEO卫星;两种外推方法对于5 min时间间隔精密星历外推精度明显高于15 min时间间隔精密星历,且三角函数插值法外推精度优于广义延拓逼近法1～2个数量级。因此,在对北斗卫星精密星历进行轨道坐标外推时,可以优先采用精度更高的三角函数插值法和使用5 min时间间隔的星历数据进行外推计算。