1:10000与1:50000比例尺基础测绘数据转换更新的探索

为进一步保持我国基础地理信息数据的现势性,开展并实现国、省两级联动更新及协同服务工作,通过探索地形数据转换更新的方法,构建数据转换模型、建立数据对应关系的方式,进行数据转换更新。本文以江苏省连云港市地级行政区域的现势性最新的省级1:10 000地形数据为主要数据源,依据《1:50 000地形数据库更新技术规定》和《江苏省1:10 000基础地理信息地形要素数据规范》开展研究,基于两种比例尺数据之间的关系,探索1:10 000地形数据转换更新1:50 000地形数据的方法,探讨其可能性与适用性,提升1:50 000地形数据更新效率,更好地为国民经济建设与社会发展提供可靠的测绘地理信息保障。     

联合钻孔资料及地形数据估算宝鸡地区vS30

将地形数据经过坡度计算,利用坡度与v_S30相关关系,计算宝鸡地区场地近似v_S30。为了得到接近实际的波速分布,借助宝鸡地区钻孔资料,研究选定区域内v_S30近似值与实测结果的差异,并对差异的空间变化趋势进行分析,在此基础上,对模型估算的v_S30进行修正,进而得到接近实际的v_S30分布,并给出基于NEHRP的宝鸡地区场地分类。     

汶川地震同震滑坡物质在震区的滞留和运移及其对龙门山地形演化的影响

2008年M_W7.9汶川地震导致龙门山断裂发生强烈地壳变形,同时引发的巨量同震滑坡加速了该地区的地表剥蚀和河流侵蚀.然而,目前尚缺少系统的数据定量研究滑坡物质的运移以及河流侵蚀速率随时间的演化规律,这些对理解龙门山前缘物质的再分配以及强震对活动造山带地形塑造的作用至关重要.为此,本研究在汶川地震后的6年间,对震区沱江上游3条支流湔江、石亭江、绵远河流域进行了多期次的定点现代河沙采样.通过系统测量河沙中的石英¹⁰Be浓度,并与震前已发表的数据进行对比,发现如下基本特点：（1）震后河沙¹⁰Be浓度均有明显降低,表明同震滑坡物质对河沙的稀释作用;（2）震后河流对河沙的运移量增加为震前的1.3~18.5倍,因此震后龙门山地区侵蚀速率短期显著增加;（3）初步估计得到汶川地震产生的滑坡物质被完全运移出造山带所需要的时间至少为100~4000年,接近龙门山地区强震复发周期;（4）震间和同震产生的构造变形和地表剥蚀在空间上具有互补性.考虑到地表剥蚀引起的地壳均衡反弹效应,认为类似汶川地震的强震有利于龙门山的隆升.认识震前、震时和震后的地壳变形及侵蚀过程有助于更好地理解单次强震事件对高原边界龙门山地形演化的作用.     

考虑地形起伏和障碍体破裂的汶川地震强地面运动数值模拟

2008年5月12日四川汶川地区发生M_W7.9地震,震中位置103.4°E,31.06°N.这次地震造成了以汶川、映秀为中心及其周边地域建筑物的严重破坏和人员的重大伤亡,且因为高山等地形复杂区域抢险救灾的艰巨性,为及时救援造成很大干扰.为更好理解地形因素对于强地面数值模拟结果的影响,建立了包含地形起伏影响及去除地形影响的两类模型.同时,依据震源破裂过程运动学反演结果,建立了包含障碍体破裂过程的震源滑动模型,实现断层分段、空间倾角以及滑移角的动态设定.基于动力学的地震动模拟方法,通过对地震波传播过程的数值计算和后处理分析,模拟由地震激发的区域强地面运动过程.结果显示：（1）强震动台站的断层距对地形效应具有放大或抑制作用,距离断层破裂带越近,地形效应越明显,反之,距离越远,则地形效应越微弱;（2）因为地形高差与障碍体的影响,地震造成的峰值可能出现在震中区域之外;（3）考虑地形影响模型的地表峰值速度（PGV）区域位于汶川与北川附近;而未考虑地形影响模型的PGV区域位于灌县—江油断层的后半段,处清平、安县附近;对汶川地震近实时强地面运动波场的模拟、峰值图谱的圈定及未来大地震强地面运动特征的预测都有重要指示意义.     

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

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

层状半空间中周期分布凸起地形对平面SH波的散射

提出了一种新的以层状半空间中周期分布斜线荷载动力格林函数为基本解的间接边界元方法,研究了周期分布凸起地形对平面SH波的散射问题.方法将散射波场分解为凸起内部散射波场和凸起外部散射波场.凸起内部散射波场通过在凸起闭合边界上施加虚拟斜线荷载产生的动力响应来模拟,而凸起外部散射波场则通过在凸起与半空间交界面上施加虚拟周期分布斜线荷载产生的动力响应来模拟.周期分布斜线荷载动力格林函数的引入,使得本文方法仅需针对一个凸起进行边界单元的离散和求解,便可完成问题的求解,避免了通过截断无限边界求解而引入的误差,方法具有较高精度的同时显著降低了求解自由度.文中通过与已有结果的比较,验证了方法的正确性,并以均匀半空间和基岩上单一土层中周期分布凸起为例进行了数值计算分析.研究表明,凸起间距对凸起地形间的动力相互作用有着显著的影响,同时层状半空间中周期分布凸起地形对SH波的散射与均匀半空间情况也有着显著的差别.     

基于开源数据和条件生成对抗网络的地形重建方法

如何使用少量的地形特征复原地形地貌一直为地学领域的难题。本文使用开源数据集提取地形特征要素,使用地形特征要素作为约束条件构建了用于生成DEM的条件生成对抗网络(Conditional Generative Adversarial Networks, CGAN),设计了基于开源DEM、开源DEM与遥感影像组合、以及5m高精度DEM提取地形特征要素生成DEM的对比实验,并对结果进行视觉效果、相关性分析以及地形因子的对比与评价。结果表明:(1)在视觉效果上,3种不同方式生成的DEM在视觉效果上均十分逼近原始5 m DEM,都远好于传统插值方法生成DEM,基于开源12.5m DEM提取要素和1m遥感影像的重建效果最接近于原始5 m DEM;(2)在相关性上,三种不同方式生成的DEM与原始5m DEM相关性均能达到0.75以上,组合开源数据提取要素重建DEM与原始5 m DEM相关性可达到0.85以上;(3)在地形因子方面,基于开源12.5 m DEM和1 m遥感影像提取要素重建DEM的坡度和坡向的分布趋势与原始5 m DEM最为一致。本文为高精度DEM建模提供了新的思路,在高精度DEM难以获取...     

局部山体地形对断层动力学破裂过程的影响研究

基于曲线网格有限差分方法研究了垂直走滑断层在不同山体地形情况下的动力学破裂传播,分析并讨论了局部山体地形对断层破裂过程及相应地面地震动的影响,得到了各模型断层面的动力学破裂过程及相应的地表峰值速度特征。研究结果表明,山体地形尺度（山体高度及底部延展距离等）对断层动力学破裂过程影响较大,进而影响到相应的地面地震动分布。当山体地形处于自由地表上亚剪切向超剪切转换的位置附近时,山体地形会阻碍断层面上自由地表超剪切的产生。一般而言,对于具有一定埋深的断层,当山体地形底部延展距离一定时,山体高度越高,其对自由地表超剪切的阻碍程度越大;当山体高度一定时,地形底部延展距离越大,越会阻碍自由地表超剪切的产生,这种破裂过程的变化会导致相应地面地震动呈现不同特征的分布。此外,还探讨了断层破裂过程及相应地震动对成核区外初始剪切应力变化的响应,结果显示,当初始剪切应力较高时,高应力降引起的超剪切破裂会对断层破裂及相应的地震动分布起主导作用。      

基于二次场的海洋可控源电磁法三维有限元正演

根据库伦规范势的定义,推导出关于磁矢量势和电标量势的偏微分方程,为了克服由电流源引起的奇异性和数值模拟计算困难,将电磁总场分解为一次场和二次场,一次场由基于Schelkunoff势函数的一维正演算法得到,二次场由有限元法计算得到,实现了海洋可控源电磁法三维有限元正演算法。通过一维数值模拟实例,验证该算法的计算精度。然后,利用该算法对带海底地形的三层储层模型进行正演,分析了海底地形对海洋控源电磁场各分量产生的影响。      

Gakkel洋中脊基底隆起的非对称地壳结构

超慢速扩张的北冰洋Gakkel洋中脊具有六个沿扩张方向的线性基底隆起(本文编号为A—F).这些线性基底隆起在中轴两侧的地球物理场和地壳结构呈现不同程度的非对称性.本文利用Gakkel洋中脊的地形、空间重力异常(FAA)和航空磁力数据,计算了它的扩张速率、剩余地幔布格重力异常(RMBA)、地壳厚度和非均衡地形.根据中轴两侧地形和地壳厚度的对称关系,我们将六个基底隆起分为对称型和非对称型两种类型.整体上,B、D和F区基底隆起在中轴两侧的地形和地壳厚度的非对称幅值(两侧差值的绝对值)较小,其中地形的非对称幅值分别为～157m、～125m、～208m,地壳厚度的非对称幅值分别为～1km、～0.06km、～0.3km;而A、C和E区的非对称幅值较大,其中地形的非对称幅值分别为～510m、～410m、～673m,地壳厚度的非对称幅值分别为～2km、～2.5km、～1.1km.我们因此推断B、D和F区具有相对对称的地壳结构,而A、C和E区具有非对称的地壳结构.根据A、C和E区中轴两侧非均衡地形的对称关系和非对称地形的补偿状态,推测A区的非对称性可能是由岩浆分配不均所导致;而C区和E区的非对称性可能是由构造断层作用使断层下盘向上抬升变薄所导致.我们进一步推测洋中脊走向的改变可能使得构造作用更易集中于基底隆起的一侧.