武汉和阿德莱德大气太阴N2潮及M2潮的对比分析
摘要：利用阿德莱德（35°S,138°E）和武汉（30.6°N,114.5°E）的流星雷达观测数据首次给出了80-100 km高度上周期为12.66太阳时的大气太阴半日潮汐（N2潮）的季节、高度、年度变化及其与周期为12.42太阳时的大气太阴半日潮汐（M2潮）的对比分析.分析结果表明：武汉和阿德莱德的N2潮和M2潮均有明显的季节、高度和年度变化.N2潮与M2潮的幅度比值大于其引力势之比0.191,在某些年份的不同季节和高度上,N2潮的幅度甚至大于M2潮的幅度.大多幅度之比接近或超过N2潮和M2潮引力势之比的2倍.中低热层的大气太阴N2潮汐值得关注.基于重力地质法的南中国海海底地形反演
摘要：根据重力地质法（GGM）,利用南中国海海域内63179个船测控制点水深将测高自由空间重力异常划分为长波参考场和短波残差场,并反演出了该海域112°E—119°E,12°N—20°N范围的1′×1′海底地形模型,该过程中使用的海水和海底洋壳密度差异常数1.32 g·cm-1通过实测水深估计得到.利用反演得到的GGM模型对剩余的10529个检核点船测水深插值计算后与实测水深进行比较,其较差结果的均值为-1.64 m,标准差为76.95 m,相对精度为4.06%.此外,根据船测点数量、分布和海底地形的不同,选择了三个海域进行统计,结果表明：在船测控制点分布均匀的海域,GGM模型精度优于ETOPO1模型,在控制点过于分散的海域其精度会有所下降,但好于船测水深的直接格网化结果.为进一步探究检核点的较差结果中出现较大数值的成因,本文对精度较差的点位进行了单独分析,选择了两条船测航迹剖面进行了研究,并分析了检核点的水深较差、相对精度与水深和重力异常的关系,结果表明：GGM模型精度受水深和重力异常的相关性影响较小,受海底地形复杂程度影响较大,地形坡度变化平缓海域的预测精度明显高于海山地区.最后,综合GGM模型和ETOPO1模型优势,利用所有船测水深作为控制,生成了综合的海底地形模型.青藏高原Moho界面起伏样式对通道流模型的动力学响应——基于数值模拟
摘要：青藏高原的高精度人工源深部地震探测发现：壳-幔边界（Moho）并非是一物理学上的＂刚性＂界面,它不仅起伏变化强烈,极为凹凸不平,而且被一系列规模不一,产状各异的深大断裂所切割,故必然会导致复杂的地表和深部物质运移动力学响应的复杂化.在常见的通道流模型中,一般均设定下地壳与上地幔之间为一平缓界面,在数值模拟工作中亦常简化为平滑的约束界面.为此,基于青藏高原实际资料提取的壳、幔介质平均速度模型,采用黏弹性数值模型讨论Moho界面起伏变换样式对通道流模型产生的响应.研究结果表明：（1）通道流效应的影响仅限于小区域内,当Moho面存在起伏样式变化时,确会对通道流产生影响,Moho界面的起伏增强了下地壳和岩石圈地幔的同步运动效应,但是其影响范围是有限的;（2）Moho界面起伏形态变化对地表和Moho界面水平位移产生的影响各异,在Moho界面发生错断的地方,呈现为地表水平位移开始发生明显加速减小的地方,即地表与深部介质水平位移解耦,模型深部动力学效应与地表的响应并非为局部性效应,而至少体现出区域性的响应.利用瑞利面波相速度和方位各向异性研究鄂尔多斯块体的岩石圈变形特征
摘要：利用在鄂尔多斯块体内部布设的45个宽频带流动台站和固定台站的资料,用双平面波方法反演了20-143 s共12个周期的基阶瑞利面波的平均相速度和方位各向异性,并反演了一维S波速度结构.反演结果显示50-100 s中长周期的瑞利面波相速度高于AK135速度模型的相速度。为高速异常,S波速度显示高速异常主要位于180 km深度范围内,表明鄂尔多斯块体保留有厚的高速岩石圈.20-111 s周期的方位各向异性强度小于1%,较小的各向异性表明鄂尔多斯块体岩石圈变形较弱.20-50 s周期的平均快波方向为近EW向,67-143 s周期的平均快波方向为NW-SE向,相对发生了整体改变,快波方向的转变约开始于80-100 km深度范围,这表明岩石圈上下部存在着由不同变形机制导致的各向异性.上部岩石圈中各向异性可能主要为残留的＂化石＂各向异性,而下部岩石圈各向异性可能是现今板块构造运动导致的变形而形成.鄂尔多斯块体岩石圈垂向上的变形差异可能主要与岩石圈温度随深度的变化以及青藏高原NE-NNE向挤压引起的上部岩石圈逆时针旋转有关.蒙古中南部地区的上地幔P波速度结构
摘要：利用中蒙国际科技合作专项项目在蒙古中南部地区布设的69个宽频带地震仪在2011年8月—2013年7月所记录的远震波形数据,采用波形相关方法在0.02-0.1 Hz的频段内拾取了18551个有效的P震相相对走时残差数据,进一步采用有限频走时层析成像的方法,反演获取了蒙古中南部下方深达800 km的P波速度结构.结果显示：肯特山下方的低速异常极有可能是推测中的肯特山地幔柱的反映;呼斯坦瑙鲁和曼达尔戈壁西边的低速异常可能与杭爱山高原下的地幔柱或地幔对流有关联;戈壁滩的低速异常带可能和达里甘嘎火山有相同的热来源。可能是该区火成岩存在的深部原因,也可能是戈壁滩及其周边地区地幔柱或下地幔热物质上涌的表现.中国大陆及周缘地震目录完整性统计分析
摘要：地震目录资料是进行地震预测、地震危险性分析、地震活动性研究等地震学研究的基础资料.对地震目录资料的完整性可靠性分析是地震学研究的基础工作之一.近几十年来我国积累的大量的仪器地震记录,历史地震也在不断更新,为地震活动性研究提供了更为丰富的样本,因此有必要对新的地震目录进行完整性分析.本文采用了新的统计方法对我国仪器记录地震目录和历史破坏性地震目录进行了完整性分析.结果表明,对于仪器记录地震目录（M≥3.0）我国东部地区1975年后基本完整,西部地区1980年后基本完整.对于历史地震目录（M≥5.0）,东部地区1500年后基本完整,西部地区1950年后基本完整.考虑到华北地震区历史破坏性地震（M≥4^3/4）活动的时间非平稳性,我们采用突变点分析法研究了华北地震区历史破坏性地震目录的时间特征,得到了华北地震区地震活动周期性变化的突变点位置,并揭示了华北地震区历史破坏性地震目录完整起始时间约为公元1500年前后.高频GNSS形变波的震相识别：模拟实验与实例分析
摘要：高频GNSS震时形变波震相及识别是GNSS地震学的重要内容.在实时数据处理基础上,本文利用振动台的高频GNSS观测实验,并结合近期部分大震的高频GNSS形变波震相特征进行研究.数据处理结果表明,实时处理与事后处理的精度在同一量级,且与采样率无关.通过与同址观测强震仪和地震计记录的对比及特定震相的频谱分析,发现高频GNSS可完整记录P波、S波、Love波及Rayleigh波震相,影响震相记录的主要因素是GNSS定位精度与震级,而仅当震中距很小时,采样率将产生一定影响.研究结果得出：基于地震波传播规律,利用高频GNSS台阵记录的形变波空间分段特征,结合震相运动轨迹及其他地震波记录,可实现实时高频GNSS形变波的震相识别.雅鲁藏布大峡谷格嘎古堰塞湖沉积岩石磁学性质指示的溃决事件
摘要：末次冰期时,在雅鲁藏布大峡谷入口处发育了一个古冰川堰塞湖,称为格嘎古堰塞湖.该堰塞湖是认识冰川阻江地貌效应的重要地质证据.古冰川堰塞湖可能会发生多次溃决洪水,通过堰塞湖沉积的分析可以获取堰塞期间的水文事件的信息.本文通过对古堰塞湖沉积39个样品多种磁学参数的详细分析发现,湖相沉积中的砂层低频磁化率（χlf）为291×10-8m3/kg,泥质沉积的χlf仅为42×10-8m3/kg;而非磁滞剩磁磁化率（χARM）和饱和等温剩磁（SIRM）的差别则较小.进一步的χ-T、Loop和IRM研究以及χARM/χlf表明,砂层中包含较多的粗颗粒磁铁矿是χlf较高的主要原因,说明砂层中的磁性矿物没有经过风化破坏,应为近源沉积,并且快速埋藏,指示了溃决事件.本文的结果表明,利用岩石磁学性质以检测古堰塞湖溃决事件是一种可行的新方法.双相介质中地震波的频率-空间域数值模拟
摘要：本文利用优化的25点频率一空间域有限差分算法对基于BISQ模型双相各向同性介质中的地震波进行了数值模拟.通过与经典的Biot模型理论模拟结果进行对比,分析了Biot流动（宏观流体流动）和Squirt流动（微观流体流动）耦合作用对地震波在孔隙介质中传播特性的影响.数值模拟在地震频段进行,结果显示：在理想相界和黏滞相界情况下,Squirt流动机制都比Biot流动机制产生了更大的速度频散和能量衰减.其中,在Biot流动和Squirt流动耦合作用下的快P波的速度和振幅小于仅考虑Biot流动影响下快P波速度和振幅,而且慢P波的衰减也更加强烈.本文还研究了地震波在双层双相各向同性介质分界面处的反射和透射特征,双相介质中波的反射与透射现象类似于单相介质的情况.模拟结果表明,利用优化25点频率-空间域有限差分法模拟双相孔隙介质中的地震波场是可行的,这为开展双相孔隙介质全波形反演问题的研究提供了可能.基于Born波路径的高斯束初至波波形反演
摘要：为了提高表层速度反演精度,本文提出了一种新的波形反演方法.该方法只利用初至波波形信息以减少波形反演对初始模型的依赖性,降低反演多解性与稳定性.由于只利用初至波波形信息,所以该方法利用高斯束计算格林函数和正演波场,以减少正演计算量.为了避免庞大核函数的存储,该方法基于Born波路径,利用矩阵分解算法实现方向与步长的累加计算.将此基于Born波路径的初至波波形反演方法应用于理论模型实验,并与声波方程全波形反演和初至波射线走时层析方法相对比,发现该方法的反演效果略低于全波形反演方法,但明显优于传统初至波射线走时层析方法,而计算效率却与射线走时层析相当.同时,相对于全波形反演,本文方法对初始模型的依赖性也有所降低.起伏地表条件下各向异性地震波最短路径射线追踪
摘要：在地震波正反演研究中,考虑起伏地表和地震各向异性具有非常重要的理论意义和实际应用价值.本文在前人研究的基础上,将最短路径追踪算法引入到起伏地表各向异性介质模型的地震波走时计算中.模型剖分时,整体模型划分成正方形单元,起伏边界附近以不规则网格逼近,进而采用非规则节点布置实现非规则网格处的最短路径计算.追踪计算中采用Sena群速度近似公式,得到各向异性地震波的走时,实现了复杂地表情况下各向异性介质模型中地震波的射线追踪.理论模型计算结果显示,本文方法能够可靠地应用于复杂各向异性介质模型,具有较高的计算精度.3D高阶抛物Radon变换地震数据保幅重建
摘要：本文结合传统3D抛物Radon变换（PRT）和AVO数据正交多项式拟合,给出了3D高阶抛物Radon变换方法（HOPRT）.该变换增加了描述AVO数据变化的梯度信息和曲率信息,拓展了传统3D抛物Radon变换方法,使其在具有AVO特征的数据重建中具有更高的准确度,从而提高AVO分析的可靠性.文中给出了3D高阶抛物Radon变换进行地震数据保幅重建的流程.理论模型和实际地震资料的重建结果显示了本文方法的优点.各向异性Radon变换及其在多次波压制中的应用
摘要：即使采用分辨率很高的双曲Radon变换,对速度各向异性发育介质及长偏移距情况下的地震数据,其Radon域内能量仍不收敛.为了克服此难题,我们在Radon变换的积分路径中考虑了非双曲走时的影响,通过引入非双曲时差公式中的各向异性非椭圆率η参数,可以准确描述出长偏移距条件下来自同一层位的时距曲线,并推导了由偏移距、慢度、非椭圆率三参数控制的积分曲线正反变换公式,我们称之为各向异性Radon变换.离散化求解时,各向异性Radon变换是时变的,频率域快速算法已不适用,本文采用了最优相似系数加权Gauss-Seidel迭代算法,保持其计算精度的同时也有较高的计算效率.将此方法应用在模型数据以及实际长偏移距海上地震数据的多次波压制处理中。收到了较好的处理效果.基于压缩感知的高分辨率平面波分解方法研究
摘要：平面波分解方法是地震资料处理中的一项关键技术,它广泛地运用在平面波偏移、各类Beam偏移等成像方法中.平面波分解不仅可以提高偏移成像的效率,而且可以压制地震资料中的随机噪音,提高地震数据的信噪比.线性Radon变换（LRT）是一种常见的实现平面波分解的方法,但常规的LRT存在以下两个缺点：（1）分辨率受测不准原理限制;（2）变换结果存在很多噪音和空间假频.为了克服LRT的上述缺点,本文提出一种基于压缩感知的高分辨率平面波分解方法,并利用加权匹配追踪（WMP）技术实现了该方法.该方法将LRT视为一个参数估计问题,并将LRT结果的稀疏性作为约束条件,在压缩感知理论的指导下利用WMP方法得到高分辨率、高信噪比的平面波分解结果.另外,该方法还可以用于提取地震数据的线性信号、压制随机噪音、实现高维地震数据规则化等地震资料处理技术.数值实验结果证明：WMP方法可以有效地提取地震数据中的线性信号,提高LRT的分辨率和信噪比,从而改善平面波分解的质量.东海陆架盆地雁荡低凸起综合地球物理解释及其成因探讨
摘要：通过对东海陆架盆地西部地震和重磁资料的综合地球物理解释,对雁荡低凸起展布形态进行了细致刻画,凸起呈NE方向不连续展布于瓯江凹陷和闽江凹陷之间,长约170 km、宽约15-50 km.地震资料揭示雁荡低凸起上广泛发育了侏罗纪与白垩纪地层。厚度约为500-1500 m,展布面积约5000 km2,局部缺失中生界地层.凸起两侧中生代盆地结构差异明显,西侧瓯江凹陷为典型的断陷盆地,东断西超、断裂发育,半地堑、掀斜断块等中生界构造样式发育;东侧闽江凹陷为坳陷型盆地,断裂、火成岩不发育,挤压背斜、断背斜、反转构造等中生界构造样式发育.自由空间重力异常图与剩余重力异常图上,凸起表现为一系列NE向团块状重力高值区,而磁力异常△T、图上则表现为深部磁场强度低的特征,火成岩影响部位可见磁力高值异常.综合凸起及邻域重磁震、莫霍面深度等地质地球物理资料,认为雁荡低凸起为一元古界组成的古隆起,区别于东部的台北低凸起.同时,结合区域构造演化及沉积特征,推测侏罗纪时期雁荡低凸起与浙闽隆起区可能连为一体,晚白垩世近东西向伸展作用下浙闽隆起发生裂陷进而形成了雁荡低凸起.西南印度洋脊中段Indomed-Gallieni洋中脊岩浆-构造动力模式
摘要：利用西南印度洋脊中段Indomed-Gallieni洋段49—51°E区段全覆盖高分辨率多波束水深地形资料,应用构造地貌学分析方法,结合区域地形及其他地球物理等资料,在分段分析49—51°E区段岩浆一构造动力学模式的基础上,进一步探讨了约10 Ma以来Indomed-Gallieni洋段的演化史.28、29洋段目前岩浆供应不足,在轴部不对称深断层的控制之下不对称扩张。属于超慢速扩张洋脊较常见的演化方式.轴部火山建造主要向北翼增生,发育与火山脊相关的火山地貌;南翼构造拉张作用强烈,地貌上可观察到大量断块,拆离断层可能大量存在.而27洋段水深浅、火山密集、轴部缺失裂谷,超慢速扩张下却具有较高的岩浆通量.Indomed-Gallieni洋段地形高地建造于一次岩浆增强事件,但应该不是因为Crozet热点的影响.27洋段为目前仍受该岩浆增强事件影响的唯一区段,但其强度和规模也在逐渐减小;包括28、29洋段在内的Indomed-Gallieni段其他部分,已重新恢复到岩浆供应不足的正常超慢速扩张洋脊演化模式.28、29洋段和27洋段岩浆供应均存在岩浆通量由多至少的周期,周期内岩浆供应较多时期轴部建脊,减少时期轴部火山建造裂离.但27洋段由于仍受岩浆增强事件的影响,与28、29洋段表现形式不同,主要表现为火山建造裂离方式、岩浆供应周期长短以及构造活动强烈程度的不同.基于勒夫波的鲁甸Ms6.5地震震源复杂性分析
摘要：2014年8月3日,云南省鲁甸县发生了Ms6.5地震,造成600余人死亡,震害之大出人预料.我们利用经验格林函数技术从国家地震台网记录的宽频带勒夫波提取了视震源时间函数,通过分析其方位依赖性获得了关于震源结构与能量辐射复杂性的认识,进而讨论了它们与震害之间的关系.对视震源时间函数的分析表明,能量释放的时间过程比较简短,83%的能量集中释放于前10 s,但这次地震发生于两个走向分别为60°和150°左右的相互交叉的断层,破裂开始于60°取向的断层,起初以双侧破裂为主,后来向60°方位的传播占优势,但很快触动了150°取向的断层,同样,起初以双侧破裂为主,后来向150°方位的传播占优势.由于两条断层同时参与了这次地震,近场能量辐射具有较强的方位复杂性,其中60°-210°方位能量辐射较强.所以,能量释放时间过程之简短但交叉的震源结构和依赖于方位的能量辐射应该是造成如此震害的主要原因.2014年8月3日云南鲁甸Ms6.5级地震序列重定位与震源机制研究
摘要：本研究采用双差定位法对2014年8月3日至7日期间鲁甸Ms6.5级主震及647个余震序列进行重新定位,得到471个重定位结果.结果显示,主震的震源深度为13.3 km,与破裂过程显示的初始破裂深度较为接近,余震序列呈现出近东西向-北西向的不对称共轭状分布。近东西向长约17 km,而北西向长约22 km,小震优势分布深度为10 km以上,且由主震处沿共轭断层分别向东南向和近东西向逐渐往10 km深度以上的浅部迁移.小震分布还展示出发震断层高倾角分布,且与昭通-鲁甸断裂分支断裂包谷垴一小河断裂活动相关.由于主震破裂的质心深度可为深入认识本次地震灾害严重提供重要证据,为此我们采用gCAP（generalized Cut And Paste）方法反演了包括主震在内共5个4.0级以上地震的震源机制解,结果显示主震质心深度仅约5.0 km,与已有破裂过程显示的较大滑移量处于2-8 km之间的深度一致.本次主震错断了互为共轭的两条断裂,这种共轭破裂模式与矩心深度较浅,可能为本次地震致灾严重的重要原因.
