利用FCT方法压制弹性波数值模拟中的数值频散

在基于波动方程的有限差分数值模拟中,会不可避免地出现数值频散（也称网格发散）问题。数值频散问题通常会给数值模拟的结果造成严重影响,因此在数值模拟中,应尽量设法消除这种现象。这里在前人的基础上,在基于各向同性介质的弹性波方程的数值模拟过程中,通过引入通量校正方法来解决数值频散问题。由数值模拟结果表明,该方法可以有效地消除数值频散现象,大大改善数值模拟的结果。     

中国县级以上行政中心迁移的特征、效应及其影响因素

行政中心设置是否科学合理,直接影响着资源要素的高效配置、基本公共服务均等化、区域协调高质量发展,以及国家治理体系与治理能力现代化。论文采用空间分析、双重差分模型等方法,对改革开放以来中国县级以上行政中心迁移的基本特征、经济效应及其影响因素进行了分析。研究表明：在迁移特征方面,行政级别越高,迁移距离相对越远,集中出现在2000年以后,以华北地区最为频繁,西部地区迁移距离普遍比其他地区更远;在迁移效应方面,新行政中心周边人口、经济和建设用地增速都明显快于原行政中心,行政中心迁移对经济发展起到了一定的促进作用,但也存在着明显的差异性;在影响因素方面,主要受资源环境条件制约、旧城规划改造推动、新区开发建设引导、行政体制改革和城市发展战略等因素的影响,行政中心迁移通常是多重因素共同作用的结果。     

青藏高原西部地表反射率的合成分析

利用中(国)日(本)青藏高原陆面过程合作试验(GAME/Tibet)设立于青藏高原西部狮泉河和改则两地自动气象站(AW S)1997~1998年观测的大气、辐射和土壤等资料,计算了两地的逐时地表反射率,并进行了日平均、月平均和日合成分析,据此讨论了青藏高原西部地表反射率的若干变化特征,并在计算了对反射率影响最大的因子——太阳高度角的基础上,着重讨论了青藏高原西部地表反射率的日变化特征。结果表明:夏半年,由于降水使土壤湿度增大,植被生长、地表反射率值较低;冬半年,受土壤湿度减小、雪盖覆盖影响,地表反射率值较高,12月平均值可达0.5以上。而地表反射率的日变化特征表现为,清晨、黄昏地表反射率高,中午地表反射率低,大致呈U形曲线,与太阳高度角的变化相反。该工作有助于深入了解高原地区地表反射率的平均状况及其变化特征。     

最优插值校准法定量估测降水及效果分析

天气雷达可以监测大范围降水以及降水分布情况,但是由于雷达估测的是离地面某一高度处的降水,会造成降水估算的误差。地面雨量站测量的是地面上某个点累计的雨量,单点精度较高,但往往漏掉强降水,暴雨中心。结合两者的优点,研究了雷达和雨量计联合估测降水方法中的最优插值校准法,并和Z—I关系估测的降水结果进行了对比分析。在此基础上还比较了不同雷达初值场对降水精度的影响。结果表明,最优插值校准法定量估测降水可以提高雷达估算降水的精度,同时雷达初值场的代表性好坏直接关系到降水估测的效果。     

Landsat遥感影像正规化处理的模型比较研究

卫星遥感影像在接收的过程中由于受到许多因素的干扰,使影像的质量受到一定程度的影响。而不同时相的影像也会因为成像时光照和大气条件的不同而出现差异。因此,在使用遥感影像之前往往需要进行必要的辐射校正。目前,主要用于辐射校正的两种基于影像的正规化模型分别是Markham和Irish等人提出的校正日照差异的正规化模型(ICM)和美国尤他州立大学遥感和地理信息系统实验室新近提出的日照大气综合校正模型(IACM/IACMt)。本文以三对不同时相的Landsat TM/ETM+影像为例,利用不变特征地物法求出各影像的相对噪音、RMSE和斜率,并用这些指标来对模型进行定量对比和评价。结果表明,当不同时相影像之间的季相、太阳天顶角和大气条件有较大差别时,经过正规化处理的遥感影像的噪音和RMSE比未经正规化处理的影像得以明显降低,且斜率更接近于1,镶嵌接缝淡化,效果更好。但如果遥感影像之间的季相一致、成像时的太阳天顶角接近,且影像晴空无云的话,则经过正规化处理的影像和未经该处理的影像之间的区别并不大。就正规化模型之间的对比而言,当影像无明显的大气影响时,ICM要好于IACM;但当受到明显的大气影响时,带有τ值的IACMt模型更好。     

随机模型对北斗DCB估计和电离层建模的影响

差分码偏差(DCB)作为电离层建模和导航定位中一项重要的误差源,对其进行估计求解至关重要.为提高北斗卫星导航系统(BDS) DCB估计和电离层建模精度,提出了一种综合高度角、卫地距和测站纬度多因素的随机模型,并对比分析了不同随机模型对BDS DCB估计和电离层垂直总电子含量(VTEC)建模精度的影响.结果表明：不同随机模型对卫星端DCB解算产生约0.2 ns差异.相较于高度角随机模型,采用高度角、卫地距组合模型测站DCB估计精度平均提高0.13 ns,电离层建模精度提高了约0.2 TECU.新提出的随机模型,在低纬度测站DCB解算精度上差于高度角模型和高度角、卫地距组合模型,但在高纬度测站DCB解算结果上更优,且对电离层VTEC建模精度提升效果明显,与前两种随机模型相比分别提升了0.88 TECU和0.68TECU.     

基于不同GNSS增强系统的激光巡线数据分析

无人直升机搭载激光雷达(LiDAR)进行输电线路巡检具有高精度、低成本的优点,为了获取高精度输电线路三维点云高精度地理空间坐标,采用基于基站差分动态后处理(PPK)、江苏省连续运行参考站(JS-CORS)、千寻位置(QX-CORS)等三种不同定位服务方式进行输电线路巡检,分析了三种定位下获取的点位定位精度以及点云质量....     

液化场地水库土石坝地震动力响应研究

基于有限差分法基本原理,采用FLAC3D软件对水库大坝在地震作用下的动力响应进行分析,研究地震强度、水库蓄水高度等对坝体稳定性的影响.结果表明:随着地震强度增加,坝基液化程度加剧,诱发的坝顶竖向沉陷增大;随着蓄水高度的增加,坝基液化程度和范围增大,坝顶竖向沉陷变形增大.当蓄水高度为10 m,地震峰值加速度为0.4 g时...     

时间-空间高阶精度矩形交错网格隐式有限差分声波正演模拟

有限差分方法因其操作简单、计算消耗低而成为地震勘探领域中最为常用的数值模拟方法之一, 然而用离散的显式差分算子数值逼近地震波动方程中的连续导数容易导致数值频散, 并且基于正方形网格离散形式的有限差分方法对不同地质模型的适应性较低.针对一阶变密度声波方程的数值模拟, 本文发展了一种适用于矩形网格离散形式的时间高阶空间隐式有限差分格式, 可以有效压制时间和空间频散, 同时灵活的网格剖分增强了其应用的广泛性.基于本文矩形交错网格时间高阶空间隐式有限差分格式的时空域频散关系和变量替换的思想, 首先采用泰勒级数展开方法求解不同方向的非轴上时间差分系数及轴上空间差分系数, 使本文差分格式可以获得任意偶数阶时间和空间精度.为了进一步提高本文差分格式在更大波数区域的空间模拟精度, 我们采用线性优化方法来求取新的轴上空间差分系数用于一阶变密度声波方程的波场迭代求解中.频散、稳定性分析及数值模拟算例表明: 相比于传统十字形空间域隐式有限差分格式, 本文矩形交错网格时间高阶空间隐式有限差分格式在精度、稳定性和效率方面均具有优势.

     

基于纵横波分离FCT弹性波正演频散压制

有限差分方法被广泛应用于地震波数值模拟和传播.传统有限差分法采用Taylor级数展开实现空间偏导数的差分, 但该方法会因为网格离散化而产生数值频散, 降低地震波模拟的精度.优化差分系数正演方法能在一定程度上压制部分频散, 然而纵、横波速度取值差异较大, 在弹性波有限差分正演模拟中, 在满足纵波最大速度确定的稳定性条件下, 浅层低速横波波场往往会产生明显的频散现象.为了削弱弹性波场正演数值频散, 提高数值模拟精度, 本文首先采用优化差分网格系数降低数值频散, 然后再采用通量校正传输(Flux-Correction Transport, FCT)法来进一步压制弹性波场有限差分数值频散.常规的FCT法是对弹性波场直接进行频散压制, 但由于弹性波场中纵、横波速度差异明显, 横波波场频散明显强于纵波, 为了压制横波波场的数值频散, 往往需要选取较大的频散压制参数, 但这会使频散较弱的纵波产生假象.因此本文提出基于纵横波分离FCT弹性波正演频散压制方法, 对分离之后的纵横波场分别选择合适的频散压制参数进行通量校正, 可以有效压制数值频散, 削弱纵波FCT产生的假象.通过理论分析和数值算例发现, 本文方法能有效削弱弹性波场有限差分数值频散, 相对于常规FCT方法没有假象产生.