根据GPS和InSAR数据反演2001年昆仑山口西地震同震破裂分布

利用干涉合成孔径雷达（interferometric synthetic aperture radar, InSAR）技术可以获取空间连续的地表形变监测结果,但仅能得到真实三维形变在雷达视线方向上的一维投影,降低了地表形变解译的可靠性。全球导航卫星系统（global navigation satellite system, GNSS）可提供高精度的三维形变结果,但其空间分辨率较低。因此可通过融合InSAR和GNSS来实现高精度高空间分辨率的三维地表形变监测。提出利用方差分量估计(variance component estimation, VCE)方法对经典的联合大地和卫星同步形变测量值的应变张量估计(simultaneous and integrated strain tensor estimation from geodetic and satellite deformation measurements, SISTEM)方法进行改进,发挥SISTEM方法可以考虑相邻点地表形变的空间相关性的优势,优化了函数模型中InSAR和GNSS观测的构成比例,引入VCE方法精确估计了各类观测值的后验方差,进而实现了高精度的三维地表形变测量。首先通过设计模拟实验对所提方法进行了验证,然后将该方法应用于2007年夏威夷基拉韦厄火山喷发案例中。相较于SISTEM方法,该方法在东西、南北、垂直向解算的形变精度均有大幅提升。     

轨迹大数据: 数据处理关键技术研究综述

针对轨迹大数据的高效点-轨迹k近邻（point to trajectory k nearest neighbor, P2T_kNN）查询处理需求,提出了一种融合时空剖分和轨迹分段的轨迹组织方法,其核心思想是在对轨迹作时间剖分的基础上,利用离散全球网格系统（discrete global grid system, DGGS）在空间上进行再次剖分,从而利用两次剖分得到的时空单元编码来索引落入其中的轨迹片段。在此基础上利用分布式列式存储技术设计了面向轨迹大数据的P2T_kNN查询处理框架,提出了一种顾及轨迹数据空间分布的自适应空间单元搜索算法,即通过分析轨迹数据在给定时间约束下的空间分异特征,动态调整空间单元的搜索步长,从而提升了轨迹稀疏区域的处理效率。针对亿级轨迹的实验结果表明,该方法适用于轨迹大数据的P2T_kNN查询处理,在轨迹稠密与稀疏区域的平均查询响应时间均小于1 s。     

省级测绘档案管理服务系统关键技术研究

结合山东省测绘档案管理服务系统的建设,理清了省级测绘档案接收、管理、服务的全业务流程,探讨了"图形、档案、库房"协同管理、基于模板的汇交管控等关键技术方法,实现了省级测绘档案"收、存、管、用"一站式服务.     

自然资源系统全流程业务平台研究与思考

分析了与测绘地理信息有关联的自然资源业务,了解自然资源业务对测绘地理信息的共性需求,分析测绘基础性公共服务平台与自然资源业务系统的关联性,梳理自然资源全流程业务关系,探索自然资源系统全流程业务平台的建设思路,为自然资源业务提供服务.     

Lippmann-Schwinger积分方程广义超松弛迭代解法及其收敛特性

为克服Born级数只适用于弱散射和短程传播的弱点,从光学散射理论中引入了针对强光学散射问题的Lippmann-Schwinger方程广义超松弛迭代解法,并对这种方法在反射地震条件下的算法表现、收敛特点、计算效率以及计算精度等问题进行了分析和讨论.理论分析和数值计算结果均表明：（1）与散射级数重整化相比,超松弛法具有更坚实的数学基础;（2）超松弛法可以有效地克服Born级数的弱点,得到与有限差分相当的数值模拟结果;（3）在反射地震条件下LS方程超松弛法表现优异,完全可以解决反射地震数值模拟中的强散射问题.

     

顾及共享交换的行业地理空间数据关联模式

在数字城市、"多规合一"、大数据建设推动下,政府管理部门和行业部门对地理空间数据共享有强烈需求.然而,其行业特殊性导致行业部门所拥有的数据无法进行有效融合与交换,数据共享成本大.以不改变行业原有的数据规范,实现行业之间地理空间数据的及时有效共享与交换为目标,通过分析行业地理空间数据特征,提炼数据的关联特征,认知并结构化...     

基于Low rank近似的黏弹性介质衰减补偿微地震逆时定位与裂缝成像

真实地下介质具有黏弹性,地震波在传播过程中会发生耗散与频散.忽视黏弹性介质的吸收衰减效应,逆时延拓过程中地震波将会出现振幅减弱、相位失真等现象,无法准确定位震源真实位置,因此需要对黏弹性介质中传播的波场进行衰减补偿,并通过采用合适的成像算子对微地震震源进行定位与裂缝成像.本文基于耗散与频散解耦的分数阶黏弹性波动方程模拟波场,采用low rank分解近似混合域算子,分离衰减相关项并反转耗散项符号,并在补偿的衰减项波场的波数域中进行低通滤波,压制噪声的影响;使用优化后的成像算子进行微地震震源定位,并通过分离散射波场,对散射波进行逆时反传寻找裂缝.数值实验证明,本文方法通过low rank近似有效提高了计算效率,衰减补偿算子在滤波器约束下能够稳定地补偿反向延拓的波场,优化后的成像算子能够在压制随机噪声的同时进一步提高计算效率和定位分辨率.

     

一种利用Nyström离散与FFT快速褶积的散射地震波并行计算方法

利用数值方法解Lippermann-Schwinger (L-S)方程的主要困难在于系数矩阵存储和线性方程组求解.这主要是因为L-S方程的积分部分是一个空间褶积,在离散后将导致一个满秩矩阵,进而形成一个大型或超大型代数方程组.因此,在利用L-S解决地震波散射问题时,一般是利用散射级数法而非数值方法.然而,散射级数法的计算精度和收敛性强烈地依赖于速度扰动的强度,而克服这种依赖性的一个可能的途径就是对现有的数值方法进行改进或是建立新的数值求解方案.在这种思想指导下,首先对L-S方程进行改写,得到一个与原L-S方程等价的积分方程(等价L-S方程).然后,对等价L-S方程进行逐点归一化处理,并利用Nyström法对经归一化处理的等价L-S方程(归一化等价L-S方程)进行离散,并用FFT计算空间褶积.之所以这样选择是由于归一化等价L-S方程经Nyström法离散生成的系数阵为一个Toeplitz阵,可利用其Toeplitz性质降低存储空间;而FFT可以将矩矢空间褶积转化为乘积,且积分核部分只要计算一次即可.进一步,为节约正演计算时间,设计了进程级和线程级相结合的MPI+OpenMP并行模式.数值试验表明,与传统的积分方程数值算法相比,利用等价L-S方程、Nyström离散和FFT快速褶积的计算方案可极大地降低存储需求,进而在保证精度的同时提高计算效率.

     

多尺度阵列嵌套组合反演宾川气枪源区横波速度结构

本文基于多尺度阵列嵌套组合的方式,利用频率-贝塞尔变换法（Frequency-Bessel,F-J方法）提取背景噪声面波频散信息,通过多个阵列融合的频散曲线反演得到宾川气枪发射台周边不同深度的横波速度结构.结果显示：浅层一阶面波频散信息的加入,使得基阶反演结果更加收敛,反演深度加深到8 km;深度在8 km以下的结构的研究利用多尺度阵列（密集台阵-宾川气枪台网-云南区域地震台网）嵌套组合的方式,面波基阶低频信息从0.55 Hz拓宽到0.008 Hz,使横波速度结构的反演深度显著增加,同时对反演过程提供约束,使得70 km深度以上的横波速度更收敛.由此本文所得的横波速度结构为该区地下结构的探测提供基础,多尺度阵列嵌套组合频散谱的研究方式也为以后区域结构的研究提供一种新的方法和思路.

     

地震地球物理数据监控软件设计与实现

台站设备故障响应和风险排除是地震监测工作的基本内容,针对其便捷性和高时效性要求,采用Java语言,开发基于Android平台的地震地球物理数据监控软件.该软件整合重庆市地震局监测工作中的各项业务需求,采用云服务、Tomcat服务器及百度地图等技术,实现了数据实时查看及历史数据波形浏览等功能,及时获取仪器设备工作状态,及...