基于环形扫面测量的三维直流电阻率法任意各向异性模型响应特征

高精度地下任意各向异性介质电性分布特征识别是当前的热门研究课题。本文以国内外前人工作为基础,基于非结构化网格,利用基于梯度恢复的后验误差估计指导网格自适应细化过程,实现了直流电阻率法三维有限元数值模拟。通过与一维模型半解析解的对比,验证了本文算法的有效性。考虑到电各向异性介质在观测中存在视电阻率反常现象,本文采用了环形扫面测量方法。通过对几种典型各向异性模型的模拟分析,得到了相应的各向异性影响规律和识别特征,各向异性主轴电阻率之间的比值大小决定了椭圆型视电阻率极性曲线长轴与短轴的比值大小,主轴电阻率的旋转方向会改变视电阻率极性曲线的形状。本文的算法研究及数值模拟技术可为直流电法数据精细处理和解释提供技术支撑。     

海洋地理信息系统中动态海面仿真方法研究

海洋地理信息系统(MGIS)中加入动态水面效果,可以有效地提升海洋环境仿真的真实感。本文从水波模型的构建出发,介绍了利用正弦波动方程模拟海面的方法;针对动态构建水面网格时高计算性能的需求,设计了基于GPU并行计算的架构;最后利用程序开发实现本文的设计思路,经过实验分析表明,该方法能够有效地模拟三维地理环境中动态海面的效果。     

利用三维高斯射线束成像进行地震定位

常规的地震定位方法通常需要拾取地震记录的初至,当初至不明显或被较高水平的噪声淹没时精度较低.本文采用基于三维高斯射线束的偏移成像方法对震源进行定位,较好地解决了该问题.通过三维高斯射线束对台站记录进行偏移归位,并将各台站成像结果的交点作为地震能量释放的中心位置;当各台站成像结果不能交于一点时,采用三维空间高斯滤波方法可实现震源位置的自动获取.提出的变网格计算方案极大地减少了计算量,显著地提高了成像精度和计算效率.利用首都圈地震台网数据,对涿鹿、滦县以及房山三个地震事件进行试算,结果表明:基于变网格三维高斯束偏移成像的地震定位方法自动化程度很高,而且具有较好的抗噪能力,特别适合处理低信噪比资料的地震定位问题.     

电导率各向异性的海洋电磁三维有限单元法正演

本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究.     

一种基于全局优化的交错网格有限差分法

在地震波场数值模拟中, 交错网格有限差分技术得到了广泛的应用, 但是在弹性模量变化较大时, 通常会因插值而导致模拟误差增大. 旋转交错网格可以很好地克服这个缺点, 因而适合于各向异性介质正演模拟. 但是对于同样大小的网格单元, 旋转交错网格需要的步长比常规交错网格要大, 这会使梯度和散度算子的误差增大因而更易产生空间数值频散. 针对这些问题, 本文提出了旋转交错网格与紧致有限差分相结合的方法, 并基于模拟退火算法进行全局优化, 压制数值频散, 拓宽波数范围. 数值模拟结果表明, 此方法可以有效地压制数值频散, 且具有较高的模拟精度.      

基于SGOG瓦块的数字地球真三维可视化技术与应用

以Google Earth为代表的第一代数字地球平台,利用球面剖分瓦片,构建了2.5维静态的数字地球表面,实现了地表空间数据的集成、建模和应用,但没有涉及地表上下的空间。本文阐述了球体大圆弧QTM八叉树剖分的规则、瓦块体系的几何特征及其编码原理,设计并实现了瓦块编码到空间直角坐标之间的互换算法,并利用C++和OSG等主流语言,构建了真三维数字地球可视化平台的雏形框架;实现了球面和球体的任意剖分、地表上下空间实体的可视化建模及简单三维空间分析。初步展示了球体瓦块的新一代数字平台在天地一体化空间大数据集成、管理、建模与应用中的优势及前景。研究表明,球体大圆弧QTM八叉树网格,具有简单、规整、几何特征明晰、适宜作为全球离散空间的坐标基准,利于空间实体建模与可视化等特点,可作为构建新一代数字地球平台的数据模型基础。     

鱼眼镜头全景影像可测量关键算法研究

鱼眼镜头具有短焦距、大视场的优点,在虚拟实景、视频监控、智能交通、机器人导航等领域得到广泛的应用.随着测绘科学技术的发展,希望在广视角图片上有更多的信息,同时通过这些信息的获取能够更好地服务社会.鱼眼镜头摄像机拍摄的影像具有非常严重的变形,需要通过一定的方法将其纠正.本文提出了一种基于球面模型鱼眼影像畸变差修正、共线条件方程及相似三角形原理的全景影像可量测性方法,使全景图像在测绘领域的应用更加广泛.     

最大准则优化技术在贴体网格中的应用

贴体网格在地质数值模拟中具有广阔的应用前景,为解决贴体网格生成时边界离散问题,提出了最大长度准则和最大面积准则,把曲线逼近和曲面网格优化问题转化为数学优化问题,为求解该问题,提出了改进的单粒子优化算法.试验表明,最大长度准则和最大面积准则的优化效果好于常规方法;以改进的单粒子优化算法求解该问题时,计算效率是智能单粒子优化算法的30倍左右（节点量为200）,从而实现最大长度准则和最大面积准则在贴体网格生成中的应用.针对最大面积准则优化曲面网格不能控制网格步长的情况,提出了限定步长的网格优化算法,使网格步长合理化,并通过实例验证了该算法的有效性.研究成果提供了生成贴体网格时边界优化准则和求解方法,对今后复杂边界的贴体网格生成具有重要意义.     

城市地理国情监测统计单元网格划分方法研究

地理国情监测是我国测绘和地理信息发展的重大战略之一。其中地理国情监测单元网格划分是进行统计分析得到地理国情监测成果的基础。针对城市若以区为统计单位,统计单元网格过大,或以社区为单位,社区界线在图上位置不是很明确。本文提出在行政区划界线（区界）的基础上,采用道路网、河流、湖泊等构建网格,采用渐进式的综合方法对网格进行聚合,得到城市统计单元网格。试验证明,该方法可作为城市统计单元网格的划分方法。     

三角网格谱元法地震波场数值模拟

谱元法结合了有限元法的灵活性和谱方法的指数收敛性,高效且高精度,是近年来发展的一种重要的地震波场数值模拟方法.经典的谱元法采用四边形(六面体)网格,利用一维Gauss-Legendre-Lobatto(GLL)积分的张量积得到对角的质量矩阵,以大大提高计算效率,但是四边形(六面体)网格不能够灵活地刻画复杂的几何模型的弯曲界面.为此,在谱元法中引入三角形(四面体)网格到二维(三维)是十分必要的.不同于经典的谱元法,在非结构化网格中不能使用GLL积分的张量积,使得非结构化网格的谱元法的实现存在着诸多的困难.目前,比较流行的三角网格谱元法,通过使用KoornwinderDubiner(KD)正交多项式,并正交化这些KD多项式构建基函数,同时利用重合的插值节点和积分节点以获取对角的质量矩阵;它所使用的积分点为优化的点集——Fekete点,且这些积分点能与四边形网格完全耦合.相比于四边形,三角网格谱元法能显著提高复杂模型的描述能力,对起伏地表模型有很大优势.本文引入高效的最佳匹配层(PML)吸收边界条件,并通过数值试验将三角网格谱元法与经典的谱元法进行对比研究.相比于经典的谱元法,三角网格谱元法显著缺点为较低的计算精度.对于7阶谱元,为了能够精确地模拟面波,三角网格谱元法需要在每个最短的面波波长内至少有11个采样点,然而经典的谱元法仅需4个采样点,并且前者所需的内存量约为后者的5.5倍.