基于快速分布式目标探测的时序雷达干涉测量方法:以Lost Hills油藏区为例

针对当前分布式目标雷达干涉测量运算效率低、选点困难等问题,本文提出了一种建立在快速同质点选取下的干涉数据处理框架.相比之前的时序数据处理方法,新方法具有选点快速、自适应性强的特点,能在保留影像分辨率基础之上增加空间点密度.另外,在统计推断的基础上,提出基于无偏空间相干性估计的分布式目标选择方法,进而弥补了传统经验阈值设定的缺陷.本文以美国加州Lost Hills油田区为例,在论证数据处理框架的可行性基础之上,分析了因孔隙流体萃取和孔隙压力降低引起的地表变形.     

2011年9月10日瑞昌-阳新地震发震构造初探

基于密集流动台阵,联合使用主事件法和 sPL 深度震相,对2011年9月10日瑞昌-阳新MS4．6地震及其余震序列进行重新定位,获得更为准确的震源位置。结果显示：重新定位后主震的水平位置为29．70°N,115．47°E,误差约1 km;震源深度为15 km,误差约2 km。结合地震序列定位结果、主震震源机制解和震区NE向断层发育的构造背景,初步探讨主震可能的发震构造。     

一种InSAR大气相位建模与估计方法

为了削弱大气延迟对干涉结果的影响以提高InSAR的测量能力,本文在InSAR大气相位特征分析的基础上,研究了一种新的InSAR大气相位建模与估计方法.首先采用稳健估计确定大气垂直分层部分的模型参数,然后利用基于Matern模型的Kriging插值估计大气紊流部分,最后应用估计的大气垂直分层和紊流资料改正InSAR测量结果.利用覆盖河南义马地区的ASAR数据对本文提出的方法进行了验证,结果表明去除大气影响后,InSAR重建的DEM与参考DEM的高程差异的均方误差由19.5m降至5.3m,精度提高了约72%.同时,改正后的干涉图更合理地揭示了义马矿区的沉降漏斗情况,进一步验证了本文方法的有效性.     

黄海及其邻近地区的Pn波速度与各向异性

利用中国东部地震台网和ISC 报告1980~2004年的地震走时数据,反演了黄海及其邻近地区的Pn波速度和各向异性,根据岩石层地幔的横向非均匀性分析了区域地质构造的深部特点.Pn波速度的变化与区域地质构造有一定的对应关系,黄海地区上地幔顶部的P波平均速度较高,没有发现明显的低速异常,表明上地幔顶部不存在大范围的地幔扰动.速度异常的分布表明,南黄海东部和西部有可能分属于不同的构造块体,其间的分界大致对应于南北走向的黄海东部断裂带,具有相对较低的Pn波速度.边界东、西两侧的Pn波各向异性存在明显的差异：南黄海西部Pn波的快波方向以北东—北北东方向为主,反映了海区内部扬子块体向北运动产生的构造变形;南黄海东部Pn波的快波方向为南北方向,与黄海东部断裂带的走向基本一致,说明黄海东部和西部之间存在一个深达岩石层地幔的南北向转换边界.结合相关资料估计黄海东部断裂带在中生代时期发生了右旋走滑运动,以响应中国东部郯庐断裂带的大规模左旋剪切以及南黄海扬子块体的向北嵌入.     

南海东北部及邻区深部结构的综合地球物理研究

南海东北部及邻区,特别是洋陆转换带地区的复杂地壳结构特征一直是南海岩石层结构研究的热点.本文在Pn波地震层析成像结果和深部地震探测剖面的约束下,利用重力数据建立该区两条剖面的密度模型.两条重力剖面二度半密度正反演的拟合结果支持琼粤隆起至吕宋岛弧区一带的地壳结构中存在下地壳高速层的观点,同时认为台西南盆地的拟合结果表明南部凹陷区仍属于过渡型地壳.本文认为剖面AA′和剖面BB′的构造属性虽然均总体倾向于火山岩型,但二者的地质结构并不完全一致,表明了北部陆缘深部结构的横向差异与构造属性的复杂.     

利用DE算法反演地壳速度模型和地震定位

利用差异演化（Differential Evolution）非线性全局优化算法,设计了一种反演地壳速度模型和进行地震定位的方法,并给出了反演结果的具体分析.利用有限差分算法计算速度模型的走时场,可以节省大量的计算量,加快计算过程.反演得到的地壳速度模型和地震的震源参数可以直接用于地震层析成像研究,还可以利用地震层析成像得到的三维速度结构对地震重新定位,从而得到较为精确的震源参数.地壳速度模型的反演方法也可以用于三维速度结构的反演.     

对黄海-东海研究区深部结构的一些新认识

综合地震层析成像与重磁数据的处理结果,选择26°N～36°N,120°E～130°E的范围作为研究区,讨论了黄、东海研究区的深部结构特点及其与周边各地质单元的相互关系,完成了研究区两条剖面的密度结构反演,认为东海陆架地区地壳厚度变化与大陆地区相比并不明显,显著减薄开始于冲绳海槽地区,中地壳消失;琉球岛弧处地壳厚度明显再度增加,特别是上地壳的厚度增加最大,推断其原因应与俯冲作用及俯冲带附近板块与地幔的运动速率之差有关.地球物理场“东西分带”是黄海—东海地区壳内结构从西向东变化的反映,但随着深度的增加,研究区的岩石层结构出现以近EW向为优势的构造格局.因此推断深部近EW向的异常是三叠纪时期南北板块碰撞、挤压所致,浅部的NE向条带异常是后期构造运动在岩石层较浅部位构造效应的反映.黄海—东海地区岩石层结构存在浅部与深部优势构造方向不协调的现象.层析成像结果证实了南黄海东缘断层的存在,还勾绘出绍兴—十万大山碰撞带为以40°左右的倾角向NW方向倾斜的高速带,另一条倾向基本相同的高速带则是南、北扬子块体结合带在深部的反映.     

环渤海地区的地震层析成像与地壳上地幔结构

利用环渤海地区的天然地震P波到时资料,采用纬度和经度方向分别为05°×06°的网格划分,反演了该地区地壳上地幔的三维P波速度结构.初步结果表明,环渤海地区地壳上地幔的速度结构具有明显的横向不均匀性：京津唐地区地壳中上部的速度异常反映了浅表层的地质构造特征,造山带和隆起区对应于高速异常,坳陷区和沉积盆地对应于低速异常;地壳下部出现大规模的低速异常与华北地区广泛存在的高导层相对应,估计与壳内的滑脱层和局部熔融、岩浆活动有关;莫霍面附近的速度异常反映了地壳厚度的变化及壳幔边界附近热状态的差异;上地幔顶部大范围的低速异常可能是上地幔软流层热物质大规模上涌所致.     

基于Pn/Pg相对定位方法研究2017年8月8日九寨沟M7.0地震起始破裂深度

本文利用Pn/Pg相对定位方法,测定了2017年8月8日四川九寨沟M7.0主震及部分余震的起始破裂深度.地震的起始破裂深度是理解地震孕震机理的重要参数,而九寨沟地区台网稀疏,地壳速度结构复杂,基于传统的到时定位方法测定地震起始破裂深度误差较大.Pn/Pg相对定位方法首先基于流动观测近台记录对2017年8月10日M4.1和11月7日M4.5余震震源位置进行测定,选择其为参考事件,再利用Pg校正主震水平位置,Pn约束震源起始深度,基于参考事件可以有效降低速度模型对震源位置测定的影响.结果显示：九寨沟主震起始破裂深度约9 km,早期余震的震源深度分布在7~13 km,主要集中在主震起始破裂深度附近.     

南海地区地震背景噪声成像与壳幔深部结构

利用南海地区28个陆地地震台站和2个布设于太平岛和东沙岛的新增海岛地震台站2011—2016年间的连续地震背景噪声波形数据,使用互相关方法计算得到了台站间的互相关函数,并提取出Rayleigh面波群速度和相速度频散曲线.采用快速行进和子空间方法反演获得了南海及周边地区12～40s周期的Rayleigh面波群速度和相速度图像,并联合反演得到了研究区深至60km的三维S波速度结构.考虑到南海数千米厚海水层对于面波频散反演的严重影响,本文在反演模型中加入了水层,显著提高了反演结果的可靠性.成像结果表明:南海及周边地区地壳上地幔顶部S波速度结构存在显著的横向不均匀性,并与这一区域的主要构造单元具有较好的空间对应关系.在5～10km深度,莺歌海—宋红盆地区的低速异常特征可能与盆地较厚的沉积层有关.在5～15km深度,海域高速异常区与海盆空间位置具有高度一致性,推测与海盆区地壳厚度相对陆缘区明显偏薄有关.当深度从20km增加至30km,海盆区的高速特征扩展至了陆缘地区,反映了地壳厚度从海盆至陆缘逐渐增厚的趋势,与OBS(海底地震仪)深地震剖面给出的地壳精细结构结果一致.至35～60km深度,海盆的高速异常特征依然明显,且速度值随深度增加整体呈现上升的趋势,推测南海海盆区的岩石圈厚度应该大于60km.