利用GRACE数据检测日本Mw9.0地震同震和震后重力变化

基于德国地学中心（GFZ）发布的GRACE RL05月重力场模型数据,考虑全球陆地同化系统陆地水储量的影响,采用300 km的扇形滤波,利用叠积法提取了日本Mw9.0地震的同震和震后重力时变信号,并利用最小二乘拟合的方法计算了两个同震重力变化极值点在地震前后85个月的重力年变率。结果表明：土壤水分和雪水引起的重力变化为-0.34~0.09 μgal;利用GRACE数据检测到的同震重力变化为-5.3~4.2 μgal,与基于PSGRN/PSCMP模型计算的结果在空间分布和量级上具有较好的一致性;震后5 a震中附近区域的重力整体上呈现增加趋势,断层上盘所在的日本海域与下盘所在的太平洋区域重力增加的最大值分别为2.6、4.5 μgal,下盘重力增加较大,可能与断层下盘所在地层的黏滞性相对较低有关。     

利用剖面重力测量数据反演深圳市地下断层参数

探测深圳市地下断层参数对该地区的规划建设、抗震设防和可持续发展具有重要意义。利用人机联作选择法反演了深圳市两条测线的断层参数,反演结果表明:两条测线的正演模型与观测值之间相对误差分别为4.51%和4.26%,计算结果可靠;断层从地下4km附近开始出现破裂信号,至17km附近消失,倾角约为70°,属于倾角较大的正断层;断层位置与深圳市地形起伏可能存在一定相关性。研究结果与地质及地球物理勘探结果具有较好的一致性,验证了剖面重力测量数据反演深圳市地下断层参数的可行性和有效性。     

矩形断层向-位错引起的空间重力变化

针对位错理论难以解释断层面之间的旋转转动的问题,该文基于子断层叠加法,利用经典点源位错理论,得出了矩形断层向-位错引起的空间重力变化。利用C++语言编制了弹性半空间中的矩形断层向-位错正演程序,模拟计算并分析了倾角为90°的矩形断层滑动和转动引起的空间重力变化的分布特征。数值模拟结果显示:断层运动引起的空间重力变化与断层的滑动和转动量大小及断层产状有关,断层滑动和转动越大,引起的地表变形范围和量值就越大;断层转动引起的重力变化大小与断层埋深成反关系,即断层埋深越大,相应的重力变化量值就越小;宏观上,矩形断层位错运动引起的重力变化呈近似对称的特征,而断层向错打破了重力变化图像的对称性。     

联合地震位错模型和InSAR数据构建2017年九寨沟Mw 6.5地震同震三维形变场

利用哨兵（Sentinel）-1A卫星升、降轨影像,在地震位错模型约束下获取了2017年九寨沟Mw 6.5地震的高质量三维形变场。首先,利用合成孔径雷达干涉测量技术（interferometric synthetic aperture radar,InSAR）提取九寨沟地震升、降轨同震形变场;然后,通过“两步法”反演获取该地震发震断层的几何参数和分布式滑动模型,以此为约束,采用方差分量估计算法联合解算九寨沟地震三维形变场。结果表明,九寨沟地震同震三维形变场以水平位移为主,垂向形变较弱;南北向形变呈拉张趋势,断层上盘向南、下盘向北滑动,最大位移分别为－19.81 cm和14.38 cm;东西向形变不对称性明显,断层上盘西北部向东水平运动,最大位移为18.37 cm,下盘东南部向西运动,最大位移不足8 cm。将南北、东西向形变与6个全球导航卫星系统（global navigation satellite system,GNSS）台站观测数据进行比较,两者一致性较好且均方根误差较小,分别为1.44 cm和1.77 cm,表明联合升、降轨InSAR观测和地震位错模型约束构建同震三维形变场方法具有较高可行性,显著降低了大地测量数据不足、InSAR观测对南北向形变不敏感等问题的影响。     

2008年大柴旦M_w6.3级地震的InSAR同震形变观测及断层参数反演

2008年11月10日,青海省大柴旦地区发生了Mw6.3级地震。本文利用EnviSat卫星升降轨SAR数据和差分干涉测量技术提取了同震形变场,基于均匀位错模型反演确定了地震断层参数,然后利用格网迭代搜索法确定了较优断层倾角,同时基于非均匀位错模型获得了精细滑动分布。结果表明,地震使得上盘区域沿降轨、升轨视线向分别产生最大约8.5cm、10cm的抬升;较优断层倾角为47.9°;地震滑动未延伸至地表,主要发生在地下8.2~23.7km深度范围内,最大和平均滑动量为0.5m和0.19m,平均滑动角为104.9°。反演的地震矩为3.74×1018 N·m,矩震级为Mw6.35。     

基于参数不确定性分析的SAR土壤水分反演精度控制方法

参数不确定性是SAR反演土壤水分的重要不确定性来源,为控制土壤水分反演精度,提出一种基于参数不确定性的有效控制土壤水分反演精度的方法,使用该方法可以控制参数的误差范围。首先使用全局敏感性分析方法,确定后向影响散射系数输出的主要参数;在不同量级高斯噪声随机扰动下,将大量各参数采值输入AIEM模型中,得到带噪声的后向散射系数集合;再使用LUT法反演土壤水分,计算反演结果满足误差量级控制范围。以此为基础,利用ENVISAT ASAR双极化数据(VV、VH)和实测土壤水分数据进行验证,利用LUT法反演得到带噪声的土壤水分,计算ASAR影像中采样点土壤水分反演值RMSE0.04cm3/cm3。结果表明各影响参数误差量级控制范围可有效控制土壤水分反演精度,在较大的入射角范围内都适用。     

利用STLN和InSAR数据反演2008年青海大柴旦Mw6.3级地震断层参数

利用Envisat/ASAR雷达影像处理后获取的2008年11月10日青海大柴旦Mw6.3级地震InSAR同震形变场数据和LN算法反演观测值含有粗差情况下的地震Okada矩形位错模型断层参数;断层深度为19.691 5 km,倾角为56.892 3°,长度为16.499 3km,宽度为6.854 8km,倾滑量为0.832 6 m,中心经度为95.884 7°,中心纬度为37.529 5°,走向方位角为116.411 5°,地震矩为3.108 4×1018 Nm(Mw6.295 0)。反演结果表明,当观测值含有粗差时,LN算法在1范数下可以有效抵抗粗差的干扰,充分利用已有先验信息的同时顾及观测、线性化等原因引起的系数矩阵误差的影响。     

利用带倾角断层形变公式研究川滇块体东边界断裂带形变特征

在经典的走滑位移公式的基础上进行改进, 引入断层倾角信息, 使其反映了走滑断层倾角对地表形变的影响。由于断层倾角的存在使得震间位移曲线的中心是以断层面闭锁-滑动分界在地表的投影为中心, 并且此处形变变化最大; 公式中的倾角项表示拟合的断层倾角, 体现了断层出露地表处与断层闭锁-滑动分界的整体偏移, 其值往往小于近地表倾角测量值。以改进的公式作为断层剖面分析的拟合公式, 应用于汶川地震前后川滇块体东边界主要走滑断裂形变特征的剖面分析研究中, 得到的拟合结果与其他方法得到的结果具有较好的一致性。改进公式比原公式更客观地反映了断层的形变状态, 但由于针对单条断层的剖面拟合受多种因素的影响, 拟合结果需结合实际情况作进一步分析。     

大气水分含量的空间异质性及其对地表温度反演的影响

大气水分含量是影响单波段热红外地表温度反演的关键因素之一。研究大气水分含量空间异质性对地表温度反演的影响误差有助于单波段热红外数据的合理应用。以MODIS MOD 05_L2产品为数据源,探讨区域大气水分含量的空间异质性及其对Landsat地表温度反演的影响。在推导大气水分含量对地表温度反演误差模型的基础上,针对Landsat5热红外波段,结合不同区域的大气水分含量分布状况,对大气水分含量的空间异质性可能引起的地表温度反演误差进行了模拟计算与对比分析。结果显示,大气水分含量空间异质性对地表温度反演的影响与大气水分含量及其异质性程度有关:当大气水分含量的空间异质性相对较大时,误差值和误差范围都比较大,达到0.5~2℃,相对误差为2%~20%;当大气水分含量空间分布较均匀时,误差值和误差范围为0.1~0.2℃,相对误差为0.5%~2%。可见,对于空间异质性程度较大的区域,区域面上大气参数的获得对地表温度的准确反演具有重要意义。     

空域最小二乘法用于重力卫星误差分析

重力测量卫星性能不仅与轨道参数、载荷误差、数据分辨率等因素密切相关，也与反演算法有关。传统的分析方法如动力学法、短弧法等用于误差分析，不可避免将算法误差引入分析结果，使得分析结论确定性不足。为解决这一问题，提出了空域最小二乘分析法，用空域格网重力扰动数据替代重力卫星载荷数据反演地球重力场，有效避免了算法误差对于分析结果的影响。分析结果表明，重力卫星在500 km轨道高度、一次数据覆盖条件下，测量重力场最高阶数约为240阶，载荷误差为1×10-10 m·s-2·Hz-1/2水平时，测量重力场最高阶数为136阶，其累积重力异常误差为2.7 mGal，累积大地水准面误差为14 cm。要达到最优测量能力，轨道倾角通常不小于89°。为减小地球引力高频信号对于地球重力场低阶位系数估计值的影响，估计位系数最高阶数需大于240阶。     

汶川地震前后川滇地区主要断裂活动趋势分析

针对川滇地区地震活动频繁的现状,该文首先对1998—2008年和2008—2015年该地区两期GNSS连续观测数据进行解算(时间节点为汶川地震);然后结合前人对川滇区域棱形块体的划分方案,基于欧拉矢量计算方法求出各块体的欧拉矢量以及旋转极经度、旋转极纬度、旋转角速度,并通过反欧拉矢量方法反演出鲜水河-小江断裂带和龙门山断裂带的错动趋势变化。结果表明:汶川地震造成两条断裂带在错动方向和错动速率方面均有变化,并且主要集中在川滇区域北部和西部的断裂带上。     

接收机自主完好性监测的算法分析

RAIM技术在接收机端实现卫星完好性监测,因此,监测算法的原理、正确性和灵敏度值得分析与研究。文中给出了RAIM故障检测和识别的数学模型,同时给出了两者的完好性保证算法。对告警门限值、误警率、漏检率之间的关系进行了较为全面的分析,分析了故障检测和故障识别功能与误警率、漏检率之间的关系。分析结果表明:可视卫星数相同时,误警概率越低,FD和FI的精度越低;误警概率相同时,可视卫星数越多,FD的灵敏度越高,FI的灵敏度越低;RAIM可用性一定时,同时满足低漏检率和低误警率是不可能的。     

一种等高线自动绘制中断层处理的新方法

提出了一种断层处理新方法--分块法与断层恢复法合成算法,这种方法兼有分块法与断层恢复法的优点,又克服了其不足,较好地解决了含逆断层时煤层底板等高线自动绘制及任意边界的构网问题,使绘出的等值线图误差较小,更加符合实际.     

DEM构建中的断裂线处理

根据断裂线在地表曲面中的数学性质,将断裂线分为硬断裂线和软断裂线,处理软断裂线时采用修改采样点权值的方法,处理硬断裂线时采用直接划分两侧采样点的方法。实验结果表明,该方法思路简单、容易实现,计算结果精度较高。     

鲜水河与龙门山断层CR布设与识别

介绍了CR的设计,并在鲜水河和龙门山断层进行了布设。针对单幅影像中CR点识别的困难,在采用GPS测量辅助CR点识别的基础上,提出利用安装CR点前后多幅影像的幅度平均对比法来识别CR点。     

小江断裂带地壳运动特征及地震危险性研究

利用1999~2013年间川滇地区的GPS速度场观测数据和块体模型研究小江断裂带地壳运动特征及地震危险性。结果表明:小江断裂带闭锁程度较高,应变量积累较大,在地表以下6km的闭锁程度约为0.94,滑动亏损速率约为(9.2±0.4)mm/a,整条断裂带左旋走滑速率约为(10.8~12.1)mm/a。结合小江断裂带上的历史地震数据估算,小江断裂带未来地震最大震级为7.3级,与近年来断裂带活动较为活跃相吻合。     

武汉地区断裂构造研究

本文运用数字图象处理技术,统计分析方法、构造地质、地震地质等多种定量、定性分析方法对武汉地区断裂构造特征及其与地震活动性的关系进行了深入研究。文中编制了断裂构造系列等密度图,指出了断裂构造的网络特征,进而对隐伏断裂进行了讨论,首次提出了汉阳——梁子湖北西向隐伏断裂的概念,最后对断裂构造与地震活动的关系进行了探讨。     

雅玛里克断裂南支断层形变监测与分析

利用跨断层短水准测量,建立了跨雅玛里克断裂南支断层形变监测网,通过资料分析得出该断裂带断层活动特征, 确定了雅玛里克断裂为全新世早期活动断裂.     

软件模块故障倾向预测方法研究

研究了在区分故障严重程度下的软件模块故障倾向预测方法,将故障分为高严重程度和低严重程度两种类型,用统计分析和机器学习方法分析静态代码度量与故障倾向之间的关系。以公开和私有两种类型的失效数据集作为实验数据,分析发现,故障的严重程度影响预测性能,预测不同严重程度的故障需要选择不同的度量和分类模型,预测低严重程度故障的性能好于预测高严重程度故障的性能。     

红河断裂带闭锁程度和滑动亏损分布特征研究

利用1999~2013年青藏高原东南缘的GPS速度场观测数据,采用DEFNODE程序反演红河断裂带走滑速率、三维闭锁程度和滑动亏损分布。反演结果表明:红河断裂带北、中、南段右旋走滑速率分别约为(5.9±1.2)mm/a,(4.8±0.6)mm/a和(4.3±0.4)mm/a,在地表以下6km的闭锁程度分别为0.43,0.22和0.25,其滑动亏损速率分别为3.3mm/a,2.6mm/a和2.3mm/a,红河断裂带闭锁程度和应变积累程度都比较低,与近年来红河断裂带整体活动微弱的现状相吻合。