一种基于格网的快速等值线充填算法

本文提出一种基于格网的等值线跟踪,适用于任意边界分割的快速充填算法,实现充填的 量化效果。根据边界线与非封闭等值线间的关系,建立等值线间的拓扑关系,并以树结构方式存储,以准确快速地实现边界线、非封闭等值线之间的封闭和封闭等值线间的嵌套。该算法已成功应用于海量多波束地形数据的成图,克服了常用多波束后处理成图软件栅格充与等值线之间的失配。     

基于网格维数法的剑川县应急避难场所演化特征研究

采用分形理论中的网格维数法对1985、2000、2014年三个时期剑川县主城区应急避难场所的空间分布特征及相关维数进行测算,得出应急避难场所不同时期空间分布的容量维和信息维。研究发现:(1)剑川县城市应急避难场所的空间演化由集中分布逐步向均衡分布发展;(2)剑川县城市应急避难场所空间分布经历了由不符合分形结构特征到完全符合分形结构特征的过程。     

地震槽波的数学-物理模拟初探

针对地震槽波在低速层的传播特性,开展了煤层内地震槽波勘探的数值模拟和物理模拟研究的初探工作.在数值模拟研究方面,采用交错网格有限差分法对煤层中的地震槽波进行三分量全波场模拟.基于波场快照和人工合成地震记录研究了不同模型中的波场特征和各种波型的传播规律.在物理模拟方面,通过选用不同配比的环氧树脂和硅橡胶类材料构建地震槽波物理模型,利用透射法和反射法观测系统获得了清晰的地震槽波记录以研究槽波的地震学特征.研究表明,在煤层内槽波的地震波场中,Love型槽波的能量小于Rayleigh型槽波的SV分量,大于Rayleigh型槽波的SH分量.相对于Love型槽波和Rayleigh型槽波的SH分量,Rayleigh型槽波的SV分量在围岩中的泄露能量较强.在煤层界面附近的围岩中,地震波仍以槽波形式传播,随着距离的增加能量逐渐衰减.随着煤层变薄,煤层槽波主频向高频方向移动,频散现象增强,传播速度增大.     

下扬子及周边地区深部泊松比结构及深部动力过程约束研究

下扬子及周边地区存在着丰富的地质构造和多金属矿产资源,其深部结构和动力学过程已成为地学界研究的热点.为了更好地讨论下扬子及周边地区的深部动力学过程和岩浆活动机制,本文基于均匀网格层析成像方法提出了非均匀网格远震层析成像方法,利用大量的天然地震相对走时残差数据反演获得了下扬子及周边地区深至700 km范围内的三维S波和P波速度结构,并根据纵横波的比值关系计算出泊松比异常.由于S波速度比P波对流体的反应更加敏感,所以泊松比异常反映了物质是否包含流体或者物质的软硬、冷热程度.本研究结果显示：（1）长江中下游成矿带下方的上地幔内存在明显的高泊松比异常,而地幔过渡带内则存在明显的低泊松比异常;（2）大别造山带及其南侧的中扬子地块的上地幔中下部及地幔过渡带内都存在明显低泊松比异常,且呈现东深西浅的空间分布特征.结合已有的地质、地球物理和地球化学等资料,我们认为长江中下游成矿带下方的地幔过渡带内滞留着古太平洋俯冲板块,其上地幔内则赋存着软的上地幔热物质,为深部成矿提供了热量或幔源物质.因此,古太平洋板块的俯冲对长江中下游成矿带的形成发挥了至关重要的作用.     

基于非结构化网格的最小二乘逆时偏移

相对于传统的逆时偏移,最小二乘逆时偏移具有更高的成像质量,这种改善是通过迭代反演来获得的,另外其精度与效率依赖于求解波动方程算法的精度与效率.本文给出了基于非结构化网格的最小二乘逆时偏移的方法,该方法能够充分结合最小二乘算法与非结构化网格精细刻画地下界面以及随速度自适应剖分的优点;并采用带补偿的拉普拉斯滤波算法,来消除梯度计算中的低频噪声,从而加速目标函数的收敛速度.通过简单倾斜模型以及复杂Marmousi模型测试,显示了该方法的有效性和潜力.     

无网格局部Petrov-Galerkin法大地电磁场二维正演模拟

有限差分法和有限单元法在大地电磁场数值模拟中已经得到了广泛的应用,但其数值结果的精度在很大程度上依赖于网格的离散程度.当模拟起伏地形、弯曲界面等复杂地电模型大地电磁场响应时,常常需要花费大量的时间以便得到较合理的离散网格.无网格局部Petrov-Galerkin法(MLPG)不同于有限差分法和有限元法,其形函数和权函数脱离了网格的束缚.本文详细推导了二维大地电磁场边值问题的弱式形式,并将其离散为局部积分域内的表达形式.通过模拟二维海洋地电模型大地电磁场响应,并与结构网格有限元结果进行对比,验证了本文算法和程序的正确性及精度.设计了一个含有弯曲界面的二维地电模型,讨论了不同离散网格对MLPG无网格法模拟结果的影响,并与结构有限元法结果进行了比较,结果表明MLPG无网格法模拟结果受离散网格影响较小.最后利用MLPG无网格法计算了两个海洋起伏地形模型的大地电磁响应,讨论了海底起伏地形对大地电磁响应的影响.     

三维复杂山地多级次群推进迎风混合法多波型走时计算

三维复杂山地条件下的各种地震波型的走时计算技术,可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算各种地震波型的走时,提出一种多级次群推进迎风混合法.该算法利用不等距迎风差分法简洁稳定地处理三维复杂地表及附近的局部走时计算问题,利用计算精度不错的迎风双线性插值法处理绝大部分均匀正方体网格中的局部走时计算问题,利用群推进法模拟三维复杂地表条件下地震波前的扩展问题,利用多级次算法处理各种类型的地震波的走时计算问题.算法分析和计算实例表明:新方法具有很好的计算精度与效率,且能灵活稳定地处理三维复杂地表复杂介质条件下的多波型走时计算问题.     

波浪作用下锚链系泊浮式防波堤动力响应的数值模型研究

浮式防波堤充分利用波能在水深方向的分布特性,在满足工程消浪要求的同时对海域水沙交换影响较小,且能够快速布置,在某些实际工程有一定应用前景。为了深入了解波浪作用下浮式防波堤的动力响应,基于OpenFOAM标准求解器olaFlow,在刚体运动求解计算中植入锚链求解模块MOODY(mooring cable dynamics),实现了基于重叠网格方法的浮体运动与锚链受力耦合求解,建立了锚链系泊浮式防波堤动力响应的二维数值模型。利用该数值模型对锚链系泊单方箱浮式防波堤在波浪作用下的透射系数、运动响应、锚链张力进行了模拟,并和相关试验结果进行了比较。结果表明,模型能够准确模拟二维波浪和浮式防波堤的相互作用,并用于三维模型的改进。     

地震灾害风险普查中地震危险性计算效率

分别用120″、60″、30″、15″、9″和6″经纬度划分网格覆盖目标区并计算网格点地震危险性,进而利用空间插值得到较小网格的地震危险性结果,通过比较直接计算的小网格点地震危险性结果和通过较大空间网格插值得到的小网格的地震危险性结果,评估由较大空间网格插值所得结果的可靠性、适用性及使用要求。结果表明,30″经纬度划分网格计算值插值得到的各种小网格的地震危险性结果的误差都小于1 Gal,是兼顾精度和速度的网格间隔。     

基于EC细网格产品在乌鲁木齐机场低能见度预测中的释用

利用2015—2018年乌鲁木齐机场航空例行天气报告(METAR报)、ECMWF(European Centre for Medium-Range Weather Forecasting)细网格数值预报产品对影响能见度的主要因子进行分析,提取与低能见度相关性高的物理量作为预报因子,采用SVM方法,分别基于Poly、RBF核函数建立乌鲁木齐机场未来21 h能见度预报模型。结果表明:(1)基于预报因子区间分类的SVM模型物理意义明确,试验结果较好;以RBF为函数建立的SVM模型(SVM-RBF)预报能力更好,其训练样本预测的TS评分0.84,准确率89.20%。(2)SVM-RBF模型的检验样本中,预报准确样本的预报误差整体偏小;在漏报样本中则有能见度越低、预报误差越大的特点,模型的振荡性明显。(3)结合NCEP/NCAR再分析资料研究SVM-RBF模型对天气过程的预报表现,发现模型对于特定天气形势下引发的低能见度天气,预报误差较小且预报提前量较大。