地球物理学报
阅读版


Earth and Planetary Physics (EPP)
地球物理学进展
地球物理学报翻译版
阅读版
  
首页 期刊介绍   编委会 期刊订阅 广告合作 留 言 板 联系我们                                                
地球物理学报
在线办公系统
作者在线投稿、查稿登录
专家在线审稿登录
编委办公登录
主编办公登录
编辑部登录
 
在线期刊
专题集New
热点文章
当期目录
下期目录
过刊浏览
ipad 移动阅读平台New
Email Alert
点击排行
下载排行
高被引论文
 
网站状态
 
  友情链接
  中国地球物理学会
  中国科学院地质与地球物理研究所
  中国地学期刊门户网
  中国知网(CNKI)
  万方数据
  中国科学D
  地球物理学进展
  Applied geophysics
  地震学报
  石油地球物理勘探
  石油物探
  勘探地球物理进展
  工程地球物理学报
  大地测量与地球动力学
  岩石学报
  地学前缘
  中国地质
  物探与化探
  物化探计算技术
  地球科学进展
  科学通报
  
中国科技期刊研究
 
文章快速检索     高级检索
 
2017 Vol.60 No.1
出版日期: 2017-01-05

目次
应用地球物理学
空间物理学·大气物理学·重力与大地测量学
地球动力学·地震学·地磁学
 
 
目次
0
2017年1月第60卷第1期目次
2017 Vol. 60 (1): 0-0 [摘要] ( 74 )  HTMLNew  PDF (684 KB)   ( 329 )
应用地球物理学
240 李振春, 黄金强, 黄建平, 朱峰
基于平面波加速的VTI介质最小二乘逆时偏移
地震各向异性集中表现为速度各向异性,势必影响地震波运动学特征.传统声波逆时偏移(RTM)和最小二乘逆时偏移(LSRTM)没有考虑介质各向异性特征,导致反射波不能正确归位、同相轴出现扭曲及寻优速度慢或不收敛等,VTI介质逆时偏移(VTI-RTM)矫正了声波成像的不足,但仍存在低频干扰严重、中深部成像不佳、振幅保持差等缺陷.为此,本文首先实现了VTI介质最小二乘逆时偏移(VTI-LSRTM)方法,为了节省I/O及内存需求并提高效率,进一步引入平面波编码技术,提出了一种基于平面波加速的VTI介质最小二乘逆时偏移(VTI-P-LSRTM)策略.在此基础上开展了简单模型及复杂Marmousi模型成像试验,并与标准逆时偏移剖面对比表明:本方法能够校正各向异性造成的相位畸变,且在迭代中自动压制串扰及低频噪声、补偿中深部能量,是一种兼具质量与效率的保幅成像策略;对速度误差的敏感性测试说明该方法需要相对正确的偏移速度及Thomsen参数模型.
2017 Vol. 60 (1): 240-257 [摘要] ( 238 )  HTMLNew  PDF (6890 KB)   ( 785 )
258 郭成锋, 杜启振, 张明强, 郑华灿, 韩冬
改进的TTI介质纯P波方程正演模拟与逆时偏移

逆时偏移作为一种高精度偏移方法已成为复杂构造成像的重要技术,描述纵波独立传播的延拓方程是各向异性介质逆时偏移的一个关键问题.在对VTI介质几个经典相速度近似公式回顾的基础上,针对常用于描述纯P波的Harlan近似公式在各向异性参数ε较大情况下近似精度较低的问题,本文对Harlan公式中的非椭圆项进行了修正,在非椭圆项前添加了一个与各向异性参数ε有关的修正系数,得到了三种改进型Harlan公式,并以近似精度最高的改进式为基础,推导了TTI介质纯P波方程.针对该伪微分方程,本文利用伪谱法和有限差分法联合实现波场延拓,对于常密度二阶方程,基于中心网格实现;对于一阶应力-速度方程则基于旋转交错网格实现.通过数值试验分析了TTI介质纯P波一阶应力-速度方程的近似精度,并以一阶纯P波方程为基础进行了TTI介质逆时偏移数值模拟试验.结果表明,本文给出的方法能够较准确地描述TTI介质纯P波波场特征,可以应用至各向异性介质逆时偏移.

2017 Vol. 60 (1): 258-270 [摘要] ( 241 )  HTMLNew  PDF (2724 KB)   ( 499 )
271 霍国栋, 杜启振, 王秀玲, 陈刚, 郑笑雪
纵向和横向同时约束AVO反演

为了提高AVO(amplitude versus offset)反演结果的精度和横向连续性,本文提出了一种新的AVO反演约束方法,该方法结合贝叶斯原理和卡尔曼滤波算法实现了对反演参数纵向和横向的同时约束.文章首先结合反演参数的纵向贝叶斯先验概率约束和反演参数的横向连续性假设建立了与卡尔曼滤波算法对应的AVO反演系统的数学模型,然后将该数学模型代入卡尔曼滤波算法框架,利用卡尔曼滤波算法实现了双向约束AVO反演.二维模型测试和实际数据测试结果表明,相对于单纯的纵向贝叶斯先验概率约束,双向约束能更准确地刻画参数的横向变化,得到更准确、横向连续性更好的反演结果.

2017 Vol. 60 (1): 271-282 [摘要] ( 146 )  HTMLNew  PDF (7649 KB)   ( 658 )
283 刁瑞, 吴国忱, 尚新民, 芮拥军, 宗兆云, 崔庆辉
三维地震与地面微地震联合校正方法
由于地面微地震监测台站布设在地表,会受到地表起伏、低降速带厚度和速度变化的影响,降低了微地震事件的识别准确度和定位精度,限制了地面微地震监测技术在复杂地表地区的应用.因此,将三维地震勘探技术的思路引入到地面微地震监测中,提出了三维地震与地面微地震联合校正方法,将油气勘探和开发技术更加紧密地结合在一起.根据三维地震数据和低降速带测量数据,通过约束层析反演方法建立精确的近地表速度模型,将地面微地震台站从起伏地表校正到高速层中的平滑基准面上,有效消除复杂近地表的影响.其次,根据射孔数据和声波测井速度信息,通过非线性反演方法建立最优速度模型,由于已经消除复杂近地表的影响,在进行速度模型优化时不需要考虑近地表的影响,因而建立的速度模型更加准确.最后,在精确速度模型的基础上,通过互相关方法求取剩余静校正量,进一步消除了复杂近地表和速度模型近似误差的影响.三维地震与地面微地震联合校正方法采用逐步校正的思路,能够有效消除复杂近地表的影响,提高微地震数据的品质和速度模型的精确度,保证了微地震事件的定位精度,具有良好的应用前景.
2017 Vol. 60 (1): 283-292 [摘要] ( 211 )  HTMLNew  PDF (3759 KB)   ( 546 )
293 谭玉阳, 李罗兰, 张鑫, 何川
一种改进的基于网格搜索的微地震震源定位方法
震源定位是微地震监测技术要解决的主要问题.目前,井下微地震监测多采用走时拟合法计算震源位置.常规方法受到环境噪声、初至拾取误差、速度模型误差等因素的影响,定位结果存在一定误差.为了提高定位精度,本文提出了一种改进的基于网格搜索的微地震震源定位方法.本文方法根据P波的偏振特征参数计算概率密度函数求取震源方位角,并采用改进的目标函数和搜索算法计算震源的径向距离和深度.模型数据和实际资料的处理结果表明,本文方法具有较强的抗噪性,计算得到的震源方位角更加接近真实值;与常规目标函数相比,本文方法采用的目标函数具有更好的收敛性,其定位结果受初至拾取误差和速度模型误差的影响更小;本文提出的搜索算法能够消除由于错误拾取造成的观测到时中的异常值对定位结果的影响.
2017 Vol. 60 (1): 293-304 [摘要] ( 235 )  HTMLNew  PDF (6411 KB)   ( 838 )
305 冯玮, 胡天跃, 姚逢昌, 张研, 崔永福, 彭更新
非稳态地震记录时变子波估计
地震子波估计是地震资料处理和解释中的一个关键问题,子波估计的可靠性会直接影响反褶积和反演的准确度.现有的子波估计方法分为确定型和统计型两种类型,本文通过结合这两类方法,利用确定型的谱分析法和统计型的偏度最大化方法,分别提取时变子波的振幅和相位信息,得到估计的时变子波.这种方法不需要对子波进行任何时不变或相位等的假设,具有对时变相位的估计能力.进而利用估计时变子波进行非稳态反褶积,提高地震记录的保真度,为精细储层预测和描述提供高质量的剖面.理论模型试算验证了方法的可行性,通过实际地震资料的处理应用,表明该方法能有效地提取出子波时变信息.
2017 Vol. 60 (1): 305-315 [摘要] ( 124 )  HTMLNew  PDF (3185 KB)   ( 671 )
316 杨春颖, 杨春, 陈双全, 李向阳, 王赟
各向异性衰减薄层地震响应特征研究
薄储层的叠前地震响应特征研究,特别是针对具有速度各向异性的含流体薄层,对储层描述具有十分重要的意义.文中基于波动方程数值模拟方法,正演得到各向同性弹性、各向同性衰减、速度各向异性、各向异性衰减模型的地震波场,并对比分析了四种模型的纵波(PP)和转换横波(PS)地震反射特征.研究结果表明:在衰减介质背景下,引入各向异性,PP和PS波的反射波振幅较弹性介质均减弱,且衰减因素对薄层振幅的影响强于各向异性.同时,VTI各向异性衰减在单频PS振幅曲线上表现出强差异性,而HTI各向异性衰减则会影响PP和PS波单频振幅曲线的极值点幅值和位置,通过分析单频振幅曲线的极值点振幅和极值点位置对各向异性衰减薄层预测有指导作用,尤其对平时较难分辨的VTI各向异性衰减薄层,单频分析方法的优势更明显.
2017 Vol. 60 (1): 316-326 [摘要] ( 182 )  HTMLNew  PDF (5702 KB)   ( 514 )
327 殷长春, 张博, 刘云鹤, 蔡晶
面向目标自适应三维大地电磁正演模拟

本文将面向目标的自适应算法应用于三维大地电磁数值模拟.使用基于非结构网格的矢量有限单元法对起伏地表大地电磁正演模拟问题进行求解.使用利用垂向电流密度在物性界面上的连续性对后验误差进行估算的算法指导网格优化.由于全局自适应算法针对观测点优化网格的能力较差,本文通过求解正演问题的对偶问题计算后验误差的加权系数,并对相关加权系数进行改进,从而实现了面向目标的自适应算法.与传统基于结构化网格的电磁正演算法相比,采用非结构网格能够更好地拟合起伏地表和地下不规则异常体.由于使用了面向目标的自适应算法,本文能够使用更少的网格达到较高的计算精度.通过对比本文模拟结果与半空间响应和全局自适应算法计算结果,并通过对比使用改进前和改进后加权系数得到的网格剖分结果验证了本文算法的有效性.

2017 Vol. 60 (1): 327-336 [摘要] ( 188 )  HTMLNew  PDF (5739 KB)   ( 819 )
337 薛帅, 白登海, 唐静, 马晓冰
耦合PML边界条件的三维大地电磁二次场有限差分法
本文将大地电磁场分解为一次场和二次场,应用交错网格有限差分法模拟计算大地电磁二次场,并引入各向异性最佳匹配层(PML)吸收边界条件作为二次场边界条件,实现了耦合PML吸收边界条件的三维大地电磁二次场有限差分正演模拟.为了确保正演的稳定性和效率,QMR求解器和磁感应矢量散度校正技术被用于PML吸收边界条件下系数矩阵的快速求解.三维模型正演响应表明,基于二次场的三维大地电磁有限差分算法具有较高的计算精度和可靠性.通过计算分析不同PML吸收因子条件的大地电磁正演结果,显示在适当的吸收因子下,PML吸收边界条件可较大幅度的减小外边界距离,从而有效的压缩模型求解空间,最终提高三维大地电磁正演模拟的效率.
2017 Vol. 60 (1): 337-348 [摘要] ( 159 )  HTMLNew  PDF (3820 KB)   ( 555 )
349 闫述, 薛国强, 邱卫忠, 周楠楠
CSAMT单分量数据解释方法
可控源音频大地电磁法(CSAMT)一直沿用大地电磁法(MT)的办法,通过计算电场分量与磁场分量的比值,求取卡尼亚视电阻率.而CSAMT场源已知,电场分量和磁场分量都与地下电阻率存在一定的关系,可以单独采用CSAMT电场分量或者磁场分量提取地下介质的视电阻率.本文通过分析电场分量与磁场分量的数据特性,提出利用CSAMT电场单分量数据进行视电阻率的计算,用改进的广义逆矩阵反演方法,使初始模型中的地电层数等于频道个数,克服了以往反演计算中层数较少的问题;实现全场区电场分量视电阻率曲线的拟合反演.同时对单分量视相位计算方法进行分析,结合山西大同地区积水采空区探测及数据解释结果,论证本文提出的单分量解释方法的有效性.
2017 Vol. 60 (1): 349-359 [摘要] ( 143 )  HTMLNew  PDF (2001 KB)   ( 461 )
360 于生宝, 郑建波, 高明亮, 栾卉
基于小波变换模极大值法和阈值法的CSAMT静态校正
采用小波分析方法进行CSAMT静态校正时,传统方法对视电阻率数据进行多尺度小波分解后,将所有尺度的细节系数设置为零,然后进行重构获得校正后的视电阻率.这使得在压制静态效应的同时,会损失一部分大构造异常的信息.针对传统小波分析方法在CSAMT静态校正中存在的问题,本文提出了基于小波变换模极大值法和阈值法相结合的静态效应校正方法.先对视电阻率数据进行多尺度小波变换,得到每一尺度上的模极大值,然后计算李氏指数进行静态效应的判断,通过设置合理的阈值,进行静态效应的压制.借助于正演模型和实测数据,验证了本方法的有效性.结果证明,本文提出的方法在压制静态效应的基础上,能够最大限度地保留大构造异常的信息.
2017 Vol. 60 (1): 360-368 [摘要] ( 111 )  HTMLNew  PDF (3936 KB)   ( 519 )
369 齐彦福, 殷长春, 刘云鹤, 蔡晶
基于瞬时电流脉冲的三维时间域航空电磁全波形正演模拟

传统时间域航空电磁全波形正演模拟主要采用间接法(褶积算法)和直接法(时域有限差分方法等),然而褶积算法需要获得精确的电流二阶导数,这给发射电流数据采集工作带来极大挑战;时域有限差分方法受到网格和时间步长的严格限制,缺乏灵活性.为解决这些问题,本文采用时域有限元方法,通过直接改变每个时间道上的瞬时电流强度模拟任意发射波形的电磁响应.由于无需计算电流二阶导数,大大提高了正演结果的精度.利用基于非结构四面体网格的矢量有限元方法和后推欧拉技术对时间域电场扩散方程进行空间和时间离散,实现三维航空电磁时间域全波形的直接正演模拟.由此不仅可以模拟复杂的地电结构,而且基于后推欧拉法的无条件稳定性,可以更加灵活地选取时间步长,提高计算效率.通过与1D数值模拟结果进行对比验证了该方法的准确性.本文对三维柱状体模型上HELITEM MULTIPULSE和VTEM系统实际发射波形电磁响应进行模拟,并与褶积算法的结果进行比较,验证了本文算法模拟实际发射波形电磁响应的优越性.对复杂三维地质体模型上不同发射波形电磁响应进行模拟,验证了时间域有限元算法可有效处理复杂地下地质结构.

2017 Vol. 60 (1): 369-382 [摘要] ( 177 )  HTMLNew  PDF (5858 KB)   ( 436 )
383 赵越, 李貅, 王祎鹏, 郭建磊, 曾友强
三维起伏地形条件下航空瞬变电磁响应特征研究

航空瞬变电磁法以其速度快、成本低、通行性好等的优势能够有效的应用于地质地形条件复杂的地区.目前对于航空瞬变电磁法的研究主要基于平坦地形的理想情况,对于地形效应的研究相对较少,然而实际应用中地形不可避免,若忽略地形影响将对资料解释造成较大的误差,从而制约航空电磁方法的进一步发展.本文基于交错网格的时域有限差分方法对三维起伏地形条件下航空瞬变电磁进行正演模拟,在保证算法准确性的前提下给出大量模型算例.以经典地形模型为例,利用所给方法计算三维正演响应,结果显示起伏对于航空瞬变电磁数据有着显著的影响且影响主要集中在早期.而后,以实际地质资料为基础,构建起伏地形条件下包含多个异常体的三维复杂模型,计算了复杂模型的航空瞬变电磁响应,并给出三维全域视电阻率曲线,从而对地形效应的影响有了更加直观的认知.最后,通过大量模型讨论了地形的尺寸参数、电性参数、飞行轨迹与飞行高度等因素变化对于航空瞬变电磁数据的影响情况,并得出有价值的结论.

2017 Vol. 60 (1): 383-402 [摘要] ( 162 )  HTMLNew  PDF (9372 KB)   ( 769 )
403 陈曙东, 郭树旭, 张爽, 袁志文, 朱俊, 张海洋, 方东
几种典型近地表有限导体瞬变电磁特征响应研究
近地表有限导体的瞬变电磁响应中包含其形状、结构特征,如何建立这些特征与其响应之间的关系是近地表有限导体探测与识别的核心问题.本文建立了几种典型有限导体响应的三维感应偶极子模型,通过有限导体早期幂函数衰减的幅度与速率、晚期指数衰减的幅度与时间常数共四个模型参数表征有限导体的形状和结构特征,利用最小二乘拟合算法提取这些模型参数,揭示有限导体的形状、结构特征与模型参数之间的内在联系.理论与实验研究结果表明:具有绕球心旋转对称性的球形目标体,其三维感应偶极子响应完全相同,早期响应呈幂函数衰减,晚期呈e指数衰减,信号幅度和时间常数均能很好反映球体半径与材料特性;具有绕轴旋转对称性目标体,其三维感应偶极子响应部分相同,早期幂函数衰减的幅度与速率主要体现目标体的长度特征,晚期e指数衰减的时间常数体现了目标体的外径与壁厚特征;非轴对称性目标体,其三维感应偶极子响应均不相同,各个感应偶极子参数差异较大.本文的研究对提高瞬变电磁系统区分轴对称目标体与金属碎片能力,降低系统探测虚警率具有重要意义.
2017 Vol. 60 (1): 403-412 [摘要] ( 83 )  HTMLNew  PDF (1968 KB)   ( 890 )
413 冯德山, 王珣
基于卷积完全匹配层的非规则网格时域有限元探地雷达数值模拟
复频移完全匹配层(Complex Frequency-Shifted PML,CFS-PML)在长时间时域计算中对凋落波、倏失波具有好的吸收效果,并被广泛应用于时域有限差分模拟中.而本文采用卷积方法将CFS-PML应用于时域有限元求解GPR波动方程的数值模拟中.论文以TM波为例,推导了基于CPML(Convolutional PML)边界的时域有限元GPR波动方程求解公式,采用Newmark-β方法对时间导数进行离散,有效改善了时域有限元GPR数值计算程序的稳定性.并以狭长模型为例,开展了CPML边界中关键参数mRκ的选取实验,通过对比反射误差大小确定了综合最优参数组合.相同时刻UPML与CPML波场快照、3个检测点的反射误差比较,说明CPML较UPML具有更好的吸收效果.最后,采用非规则四边形网格对1个复杂GPR模型进行剖分,应用加载CPML边界条件的FETD程序对该模型进行了正演,得到了二维剖面法、宽角法正演GPR剖面图,说明非规则四边形对复杂模型的良好适应性,基于CPML边界条件的FETD可有效减少边界反射误差,能实现对任意复杂不规则模型的正演模拟.
2017 Vol. 60 (1): 413-423 [摘要] ( 164 )  HTMLNew  PDF (3773 KB)   ( 614 )
424 曹清华, 闫述, 薛国强, 朱娜
高分辨地电阻率探测与互联网络远程监测煤田火烧区
为发现、治理、监控煤田火烧区,将互联网的创新成果深度融合于传统地勘技术领域之中,开发了具有探测火烧区范围和中心温度、连续采集、高温传感、远程遥控等主要功能的永久性的无线传感器网络远程监测系统.研发出新型传感器节点,能够与单极-偶极装置互联储存并发送探测数据,同时将高温探头直接放置于地下火烧区巷道进行监测;通过带保护电路的太阳能供电装置,为现场数据采集网络提供了持续充足的能源;由移动通信网和互联网将信号传送到监测中心,实现远程控制监测;在新疆地区的12个火烧区同时部署了所研发的监测系统,至今已工作45个月,表明采取的技术措施可以保障监测系统的免维护长期稳定运行.乌鲁木齐监测中心收到110多万条现场信息,授权后可任意查询、永久保存,为新疆大面积煤田自燃火区的治理和保卫治理成果,提供基础数据.研究结果表明,充分发挥互联网在地球物理探测监测中的优化和集成作用,对提升煤田地质基勘查的创新和生产能力,具有重要的意义.
2017 Vol. 60 (1): 424-429 [摘要] ( 67 )  HTMLNew  PDF (1721 KB)   ( 346 )
430 啜晓宇, 滕吉文
强导(含)水隐伏陷落柱底板突水机理研究
华北型煤田采场底板隐伏强导(含)水岩溶陷落柱突水以其隐蔽性强,突水速度快,水量大,破坏性强,危害严重的特点已成为影响煤矿安全生产的重大隐患.但是,受奥陶系灰岩承压水和采动应力耦合作用下的岩溶陷落柱突水机理、突水过程和突水力学规律的研究尚不够系统和完善.本文以陷落柱上覆岩层层状介质假设和剪切破坏理论方法为出发点,建立了不同水平截面形状类型的强导(含)水岩溶陷落柱突水力学模型,并确立了奥灰承压水与采动应力共同作用诱发强导(含)水陷落柱突水的临界条件,并提出了突水危险性系数对陷落柱突水危险性进行定性和定量的评价.理论推导和数值分析表明,奥灰承压水和采动应力耦合作用下的强导(含)水陷落柱突水与其水平截面形状、隔水层厚度、岩石物理性质有关.当采场环境确定后,陷落柱水平截面形状则是决定其突水的关键因素,而且陷落柱水平截面形状是圆形时最容易发生突水.
2017 Vol. 60 (1): 430-440 [摘要] ( 102 )  HTMLNew  PDF (6364 KB)   ( 376 )
441 张超, 刘国强
基于周向磁场有限元方法的双侧向测井影响因素分析
电流线分布是评价双侧向测井仪器探测性能的一种辅助手段.通常采用基于电位的有限元数值模拟方法,通过对电位的微分获得电流线的分布,这种处理方法引入了额外误差.为解决此问题,本文提出了一种基于周向磁场强度的有限元方法,通过对视电阻率和电流线的数值模拟验证了新方法的正确性,且新方法在计算电流线方面优于电位法.在此基础上,针对泥浆低侵和高侵两种情况,模拟了深浅侧向测井的视电阻率响应和电流线的分布,分析了侵入带对双侧向测井响应的影响规律.
2017 Vol. 60 (1): 441-450 [摘要] ( 125 )  HTMLNew  PDF (3627 KB)   ( 450 )
451 国家自然科学基金委员会地球科学部
2016年国家自然科学基金资助项目——地球物理与空间物理部分
2017 Vol. 60 (1): 451-461 [摘要] ( 153 )  HTMLNew  PDF (540 KB)   ( 1053 )
空间物理学·大气物理学·重力与大地测量学
1 孙文杰, 宁百齐, 赵必强, 李国主, 胡连欢, 常首民
2015年3月磁暴期间中国中低纬地区电离层变化分析
2015年3月17日爆发了本太阳活动周最大的地磁暴,Dst指数达到-233 nT.本文利用电离层测高仪foF2hmF2、北斗同步卫星(BDS GEO)TEC以及GPS电离层闪烁S4指数对此次磁暴期间中国中低纬地区(北京、武汉、邵阳和三亚)的电离层变化进行分析,并对此次磁暴所引发电离层暴的可能机制进行了探讨.磁暴期间,中低纬电离层暴整体表现为正相暴之后长时间强的负相暴.3月17日白天中纬正相暴为风场抬升电离层所致,而驼峰区及低纬地区正相暴由东向穿透电场所引起;3月18日白天长时间的强负相暴为西向扰动发电机电场和成分扰动所引起;3月17和18日夜间的负相暴可能是日落东向电场受到抑制以及赤道向风场对扩散的抑制导致驼峰向赤道压缩所致,同时被抑制的日落东向电场强度不足以触发产生赤道扩展F,导致低纬三亚和邵阳夜间电离层闪烁在磁暴期间受到完全抑制.这是我们首次基于北斗同步卫星TEC组网观测开展的电离层暴研究.
2017 Vol. 60 (1): 1-10 [摘要] ( 219 )  HTMLNew  PDF (4560 KB)   ( 823 )
11 袁媛, 周宁芳, 李崇银
中国华北雾霾天气与超强El Niño事件的相关性研究
2015年11-12月,全国接连发生七次大范围、持续性雾霾天气过程,其中,11月27日-12月1日的雾霾天气过程持续时间长达五天,成为2015年最强的一次重污染天气过程;12月19-25日重度雾霾再次发展,影响面积一度达到35.2万km2.本文利用多种数据资料通过个例对比和历史统计详细分析了超强El Niño背景下雾霾天气频发的天气气候条件.其结果清楚表明:2015年11-12月欧亚中高纬度以纬向环流为主,东亚冬季风偏弱,使得影响我国的冷空气活动偏少,我国中东部大部地区对流层低层盛行异常偏南风,大气相对湿度明显偏大,并且大气层结稳定,对流层底层存在明显逆温.上述大气环流条件使得污染物的水平和垂直扩散条件差,因此在有一定污染排放的情况下,造成了重度雾霾天气过程的频发.由此,超强El Niño事件所导致的大尺度大气环流异常是我国中东部,尤其华北地区冬季雾霾天气频发的重要原因之一.
2017 Vol. 60 (1): 11-21 [摘要] ( 189 )  HTMLNew  PDF (6945 KB)   ( 1111 )
22 张亦洲, 苗世光, 李青春, 戴永久
北京城市下垫面对雾影响的数值模拟研究
为了探究北京城市下垫面对雾天气过程的影响,为北京地区雾数值预报水平的提高提供理论基础和科学依据,选取2011年10月29日北京地区雾天个例进行了数值模拟试验,通过对WRF/Noah/UCM模式系统中城市冠层参数的调整,显著改善了模式对此次雾天气过程的模拟效果.使用参数调整后的模式系统通过敏感性试验分析研究了北京城市下垫面对雾发生、发展和消散过程的影响.结果表明:参数调整后的WRF/Noah/UCM模式系统能够与实际观测较相符地模拟此次发生在北京地区的雾天气过程,北京城市下垫面主要通过对温度的改变对雾的形成、发展和消散产生显著影响,使雾不易在城市及其附近形成和发展,延后城市地区雾的形成,但城市的存在也使得城市地区及其附近雾不易消散,相较于没有城市时消散时间延后.
2017 Vol. 60 (1): 22-36 [摘要] ( 223 )  HTMLNew  PDF (9394 KB)   ( 546 )
37 杨帆, 王长青, 许厚泽, 钟敏, 周泽兵
利用运动学轨道提高GRACE时变重力场解算
基于变分方程法,本文利用GARCE高精度K波段星间测速数据KBRR,结合德国格拉茨大学发布的运动学轨道和GFZ发布的简动力学轨道作为两种伪观测值,分别解算了2005-2010年60阶全球时变重力场模型Hust-IGG01与Hust-IGG02.通过与GRACE官方机构发布的模型和其他国际主流权威模型进行对比,发现基于运动学轨道结合KBRR解算的模型Hust-IGG01优于基于简动力学轨道结合KBRR解算的模型Hust-IGG02:在重力场系数C20时间序列的统计数据上,Hust-IGG01比Hust-IGG02更接近SLR结果,在如C60C70C80以及C90等重力场低阶项上的数学统计均更接近CSR RL05;Hust-IGG01的重力场系数误差分布和GFZ RL05在同一水平,而Hust-IGG02的误差估计过于乐观;Hust-IGG02在主要质量变化区域上存在5%~10%信号低估,而Hust-IGG01能完全达到国际主流机构利用GPS观测数据的解算水平,Hust-IGG01与官方机构CSR、JPL和GFZ最新模型在格陵兰岛的冰川消融年际趋势分别是-125.4、-125.4、-127.3、-124.3 Gt·a-1,在亚马逊流域的平均等效水高周年振幅分别是17.56、17.40、17.46、17.22 cm,在撒哈拉沙漠的平均等效水高均方差分别是0.87、0.77、1.10、0.87 cm;另外在Hust-IGG01的实际应用上,本文分析了全球32个主要流域质量变化的年际趋势、周年振幅和半周年振幅三种信号模式,统计结果显示Hust-IGG01与CSR RL05结果基本吻合.
2017 Vol. 60 (1): 37-49 [摘要] ( 326 )  HTMLNew  PDF (4767 KB)   ( 968 )
50 苏勇, 于冰, 游为, 肖东升, 王继燕, 秦毅坤
基于重力卫星数据监测地表质量变化的三维点质量模型法
利用卫星重力测量手段监测全球质量变化取得了巨大成功,本文基于牛顿万有引力定律在三维空间直角坐标系中导出利用重力卫星观测数据监测全球质量变化的三维点质量模型法,该方法可直接利用重力卫星的轨道和星间观测数据或时变重力场模型计算全球质量变化,由于利用卫星观测数据计算地表质量变化的向下延拓过程以及观测数据噪声的影响,需要采用合适的空间约束方程或正则化技术对解算结果进行约束或平滑处理.利用合成全球质量变化模型模拟一个月的GRACE双星轨道和星间距离变率数据计算全球质量变化,对三维点质量模型法进行分析验证,采用零阶Tikhonov正则化技术处理病态问题.结果表明,三维点质量模型法可有效用于重力卫星观测数据监测全球质量变化,为利用重力卫星观测数据监测全球质量变化提供一种可选的途径.
2017 Vol. 60 (1): 50-60 [摘要] ( 190 )  HTMLNew  PDF (7890 KB)   ( 534 )
61 吴志露, 刘焱雄, 何秀凤, 彭琳, 郑剑
基于近岸海岛GPS数据反演海洋分潮负荷影响
本文探索了海潮负荷特征值的反演方法.基于中国近岸海岛GPS站(平潭与闸坡)观测数据,采用FFT方法提取了海潮负荷特征值,分析了反演特征值与全球海潮模型FES2004、NAO.99b和GOT4.7计算出的相关特征值之间的差异,评估了反演海潮分潮频率项的精度,并利用平潭站得到的海潮负荷特征值对附近的三沙站进行海潮改正以评价反演效果.结果表明:(1)频谱分析可精确提取4个半日分潮和4个全日分潮负荷的频率信息;与已知频率相比,P1分潮的反演频率误差为1.4%,其他7种分潮负荷反演频率误差均小于1%.(2)两个海岛分属不同的潮波入侵通道,反演分潮振幅和初相存在差异,但反演分潮频率几乎一致,间接证实它们属于同一潮波系统,也表明其潮差有别.(3)反演振幅与三种全球模型具有较好的一致性;其中S2、O1、P1、Q1四个分潮在水平方向互差为1~2 mm,高程方向上的互差均小于3 mm;K2、K1、M2、N2在水平方向振幅互差多数小于2 mm,个别差异高达4 mm,高程方向互差多为5~6 mm,个别超过10 mm.(4)反演得到海潮负荷改正模型相对于3种全球模型在三沙站的改正效果略佳,间接表明反演结果有效、可靠.(5)动态PPP结果中虽然存在多种误差,其时间序列仍可分离并提取海潮负荷的影响.
2017 Vol. 60 (1): 61-69 [摘要] ( 204 )  HTMLNew  PDF (4189 KB)   ( 341 )
地球动力学·地震学·地磁学
70 彭淼, 姜枚, Chen Youlin, 谭捍东, 李庆庆, 张立树, 许乐红
利用远震接收函数揭示的喜马拉雅东构造结台阵下方地壳结构及其动力学意义
在喜马拉雅造山带的东缘,雅鲁藏布江缝合带在这里发生急剧转折,南迦巴瓦变质体快速隆起,然而关于东构造结的形成机制一直未有定论.利用围绕南迦巴瓦峰的48个宽频带地震台站记录的远震数据提取P波接收函数,采用改进的H-κ叠加方法和共转换点叠加方法综合研究了东构造结的地壳厚度、波速比分布和地壳结构特征.结果表明:研究区平均地壳厚度为64.03 km,大部分台站介于60.48~66.55 km范围;平均波速比为1.728,主要集中范围为1.696~1.742.东构造结地壳厚度横向变化剧烈,构造结西端和北端厚而中间薄,东构造结核部Moho面呈现上隆的构造形态,东西向上隆幅度约为6~7 km,南北向的上隆超过9~10 km.东构造结核部地壳上隆减薄可能由高密度、高波速的岩石圈撕裂残片拆沉到上地幔软流圈后重力失衡所致.平均波速比超过1.8的高值异常展布于东构造结的两侧,推测为环东构造结的壳内部分熔融体.东构造结地壳上隆减薄和壳内部分熔融的存在很可能均与幔源热物质的上涌有关,而软流圈地幔的上涌则可能由印度板片的撕裂引起.
2017 Vol. 60 (1): 70-85 [摘要] ( 255 )  HTMLNew  PDF (7833 KB)   ( 1101 )
86 郭慧丽, 丁志峰, 徐小明
南北地震带北段的远震P波层析成像研究
本文利用“中国地震科学台阵”探测项目在南北地震带北段布设的678个流动地震台站在2013年10月至2015年4月期间记录到的远震波形数据,经过波形互相关拾取到473个远震事件共130309条P波走时残差数据,通过远震层析成像研究获得了该区(30°N-44°N,96°E-110°E)下方0.5°×0.5°的P波速度扰动图像.结果显示,研究区下方P波速度结构显示强烈的不均一性和显著的分区、分块特征.岩石圈速度结构具有显著的东西差异:祁连、西秦岭和松潘甘孜地块组成的青藏东北缘地区显示明显的低速异常,而属于克拉通性质的鄂尔多斯地块和四川盆地则显示高速异常,表明东部克拉通块体对青藏高原物质的东向挤出起到了强烈的阻挡作用.阿拉善地块显示出弱高速和局部弱低速的异常并存的特征.阿拉善地块西部显示低速异常,而东部与鄂尔多斯相邻的地区显示高速异常,可能表明该地区的岩石圈的变形主要受到青藏高原东北缘的挤压作用.在鄂尔多斯和四川盆地之间的秦岭下方100~250 km深度上表现为明显的低速异常,表明该处可能存在软流圈物质的运移通道.鄂尔多斯北部的河套裂陷盆地下方在100~500 km深度内低速异常表现明显,说明该区有深部热物质上涌且至少来源于地幔过渡带.青藏东北缘上地幔显示低速异常且地幔过渡带中出现明显的高速异常,这种结构模式暗示了在青藏高原东北缘可能发生了岩石圈拆沉作用,而高速异常体可能是拆沉的岩石圈地幔.
2017 Vol. 60 (1): 86-97 [摘要] ( 311 )  HTMLNew  PDF (8273 KB)   ( 1108 )
98 高雅健, 崔辉辉, 周元泽
中天山和塔里木盆地下方地幔转换带顶部P波速度结构探测
本文基于中国数字地震台网记录到的2016年1月12日发生于兴都库什-帕米尔地区的中源地震P波宽频带波形资料,利用P波三重震相方法研究了中天山和塔里木盆地北部410 km间断面起伏形态及间断面附近P波速度结构.通过波形拟合,我们发现中天山和塔里木盆地北部410 km间断面下沉约10 km;410 km间断面上方存在由塔里木盆地向中天山P波速度异常逐渐减小的约70 km厚的低速层(-3%~-1%).该低速层可能是地幔转换区物质上涌并部分脱水熔融及热效应共同作用的结果.另外我们发现中天山东北部约10~100 km深度范围内存在局部低速异常(约-6%),其很可能是上涌的软流圈物质.
2017 Vol. 60 (1): 98-111 [摘要] ( 190 )  HTMLNew  PDF (9413 KB)   ( 912 )
112 汪晟, 胥颐, 孟晓春
渤海海峡及周边区域地壳结构的层析成像特征
渤海海峡是连接中国东部山东半岛和辽东半岛的重要途径,其跨海通道的地壳稳定性研究受到高度关注.本文利用地震层析成像方法重建三维P波速度模型,揭示了渤海海峡及周边区域地壳和上地幔的构造特征.结果表明,渤海海峡的速度结构存在明显的非均匀性,海峡北部地壳速度较高,结构较为完整,断层活动不明显,与现今较弱的地震活动相吻合,但是地壳底部存在低速薄层,它有可能成为地壳和上地幔之间的滑脱带,需要开展进一步的研究加以确认.相比之下,海峡南部地壳速度偏低,附近区域地震活动频繁,与张家口-蓬莱断裂带通过于此有着密切的联系,该断裂持续不断的地震活动对海峡南部的地壳结构产生了较大的影响.在渤海南部,郯庐断裂带东、西两侧的地壳结构明显不同,西侧速度偏高,东侧至渤海海峡速度偏低,这一特征可能与此地区广泛发育的断层和地震活动有关.另外,受华北克拉通破坏及地幔上涌的影响,渤海地区地壳深部和上地幔速度偏低,郯庐断裂带及渤海海峡附近显示出深部热流的活动迹象,反映了岩石圈减薄和软流圈的局部抬升.
2017 Vol. 60 (1): 112-122 [摘要] ( 220 )  HTMLNew  PDF (6800 KB)   ( 918 )
123 滕吉文, 马学英, 张雪梅, 刘有山, 皮娇龙
2015年尼泊尔MS8.1大地震孕育的深层过程与发生的动力学响应
2015年4月25日发生在尼泊尔博克拉MS8.1大地震的深层动力过程与“地中海-喜马拉雅-南亚地震带”的中段喜马拉雅地震活动带密切相关.这次大地震是该地震强烈活动带上长期以来深部物质与能量强烈交换、运动,并导致构造活动和应力积累的产物.综合分析与研究提出:(1)博克拉MS8.1大地震的孕育、发生和发展具有长期活动和近年来相邻地带地震活动频繁的背景;这一地带自1505年-2015年,即500多年来相继发生多次MS≥8.0的大地震.(2)这一地带具有特异的地球物理边界场响应和深层动力过程,显示深部物质的重新分异、调整与能量交换.(3)大地震发生与周边地带应力场分布特异,壳、幔结构与介质属性变异及破裂响应与断层面解的属性相关.(4)喜马拉雅地带的三条北倾断裂带以不同角度向深部延伸、震源位置及浅表层的变形特征尚应深化理解.(5)MS8.1大地震的发生对相邻地带的波场影响强烈,故应强化高精度地球物理场的观测和探测,以“捕捉”未来可能大地震的孕育与发生.
2017 Vol. 60 (1): 123-141 [摘要] ( 133 )  HTMLNew  PDF (12435 KB)   ( 470 )
142 王晓欣, 丁志峰, 马延路
利用非线性台阵叠加方法快速追踪2015年4月25日尼泊尔MW7.9地震破裂能量中心运动轨迹

本文应用一种非线性台阵叠加方法,通过对地震记录瞬时相位信号的延迟叠加来估计信号间的一致性程度,用于在震后快速追踪大地震破裂能量中心相对于震中运动轨迹.基于此方法,文中利用欧洲地震台网不同频率段的数据对2015年4月25日尼泊尔MW7.9地震的破裂轨迹进行了追踪,并对成像结果的不确定程度进行了估计,且通过理论地震图计算,验证了该结果的可靠性.研究结果表明,该事件为东南向发展的单侧破裂,破裂尺度至少达到120 km.整个破裂过程主要分为三个阶段,第一阶段破裂在震中附近释放小部分能量;第二阶段破裂向东南方向发展,表现出明显的沿断层倾向辐射频率变化的现象,相对更高频的能量辐射在更深处发生;第三阶段破裂继续向东南方向发展,不同频段得到的能量中心运动轨迹趋于一致,这可能提示了滑移范围的逐渐变窄.同时,推断在主喜马拉雅逆冲断裂的更深部(断层深度大于等于20 km处)也存在破裂的可能.

2017 Vol. 60 (1): 142-150 [摘要] ( 241 )  HTMLNew  PDF (4707 KB)   ( 1285 )
151 屈春燕, 左荣虎, 单新建, 张国宏, 宋小刚, 刘云华, 余露
尼泊尔MW7.8地震InSAR同震形变场及断层滑动分布
采用DInSAR技术和欧空局2014年新发射的Sentinel-1A/IW数据,获取了2015年4月25日尼泊尔MW7.8地震的InSAR同震形变场.所用InSAR数据扫描范围东西长约500 km,南北宽约250 km,覆盖了整个变形区域,揭示了形变场的全貌及其空间连续变化形态.此次地震造成的地表形变场总体呈现为中部宽两端窄的纺锤形,从震中向东偏南约20°方向延伸,主要形变区东西长约160 km,南北宽约110 km,由规模较大的南部隆升区和规模较小的北部沉降区组成,南部最大LOS向隆升量达1.1 m,北部最大LOS向沉降量约在0.55 m.在隆升和沉降区之间干涉纹图连续变化,没有出现由于形变梯度过大或地表破裂而导致的失相干现象,表明地震断层未破裂到地表.基于InSAR形变场和部分GPS观测数据,利用弹性半空间低倾角单一断层面模型进行了滑动分布单独反演和联合反演,三种反演结果均显示出一个明显的位于主震震中以东的滑动分布集中区,向外围衰减很快,主要滑动发生于地下7~23 km的深度范围内.InSAR单独反演的破裂范围,特别是东西向破裂长度大于GPS单独反演的破裂长度,而InSAR单独反演的最大滑动量则低于GPS单独反演的滑动量.因此认为联合反演结果更为可靠.联合反演的破裂面长约150 km,沿断层倾向宽约70 km,最大滑移量达到4.39 m,矩震级为MW7.84,与之前用地震波数据和GPS数据反演的结果一致.
2017 Vol. 60 (1): 151-162 [摘要] ( 307 )  HTMLNew  PDF (9744 KB)   ( 824 )
163 黄禄渊, 张贝, 程惠红, 罗纲, 瞿武林, 石耀霖
2015年智利Illapel MW8.3地震同震效应及其对南美大陆地震危险性影响
基于高性能并行有限元方法,建立含地表地形和Moho面起伏的大规模非均匀椭球地球模型,计算了2015年9月17日智利Illapel MW8.3地震同震效应,同时讨论2010年智利Maule MW8.8地震和2014年智利Iquique MW8.1地震对2015年智利Illapel地震的加载作用,进一步分析这3个地震共同作用对南美大陆及周围断层的地震危险性影响.结果表明:对此次智利Illapel地震CEA(中国地震局地球物理研究所)滑动模型计算的地表最大水平位移约为3 m,USGS(美国地调局)模型的最大水平位移约6 m.特大地震影响范围广,南美大陆接近一半区域同震水平位移量级达到0.5 mm;受3次地震共同作用,智利2010年Maule地震南部俯冲带、2015年Illapel地震和2014年Iquique地震之间俯冲带,以及2014年Iquique地震以北俯冲带,受到10 kPa东西向压应力加载;2010年Maule地震和2014年Iquique地震共同加速了2015年Illapel地震滑动;1960年Valdivia MW9.5地震破裂区、2014年Iquique地震和1995年Antofagasta地震之间区域以及2010年Maule地震和2015年Illapel地震的未破裂区重合处库仑应力变化均超过10 kPa,有可能加速强震的发生.
2017 Vol. 60 (1): 163-173 [摘要] ( 281 )  HTMLNew  PDF (6885 KB)   ( 488 )
174 邓园浩, 程惠红, 张怀, 瞿武林, 张贝, 石耀霖
2016年3月2日苏门答腊MS7.8地震同震位移和应力场数值模拟研究
苏门答腊地区发育多条左旋走滑性质断层,地震活动活跃.2006年3月2日该区西南海域发生了MS7.8大地震.大地震的发生常常会引发区域位移场和应力场及周围断层应力状态发生变化.本文建立全球PREM有限元地球模型,据已有的断层滑动模型计算了此次苏门答腊地震引发的同震位移和应力及库仑应力变化,并进一步讨论了此次地震对周围断层的影响以及区域构造应力场对库仑应力变化计算的影响.初步结果表明此次苏门答腊MS7.8地震造成较大的南北向水平位移且集中在探测者破裂区(Investigator Fracture Zone),最大水平位移量约6.74 m,断层倾角接近垂直,下盘向北运动而上盘向南,进一步表明MS7.8地震为典型的左旋走滑为主的地震,发生海啸的可能性较低;库仑应力变化达MPa量级的区域集中在震中,但近场大部分余震分布在库仑应力减小区域,有效摩擦系数变化和区域构造应力场的耦合作用可能是其原因;利用改进的库仑应力变化计算方法和最优破裂方向计算得出的结果显示库仑应力触发理论可较好地解释余震分布.
2017 Vol. 60 (1): 174-186 [摘要] ( 290 )  HTMLNew  PDF (5455 KB)   ( 537 )
187 郑露露, 林建民, 倪四道, 祝捍皓, 郑红
台风激发的第二类地脉动特征及激发模式分析
0.003~1 Hz频段的地脉动主要来源于海浪运动与固体地球的耦合作用,台风引起的强烈海浪运动往往可使地脉动能量显著增强.由于涉及大气-海洋-固体地球三个圈层之间的复杂动量传递与耦合过程,迄今为止,关于台风激发地脉动的具体源区位置及激发机制尚存在争议.本文选取日本、中国东南沿海及台湾地区的地震台站波形连续记录,研究了2008年台风“森拉克”和“黑格比”激发地脉动的时频特征,开展相应数值模拟,并与观测数据进行了对比分析研究.结果表明台风激发第二类地脉动存在两种主要模式:(1)近岸源区激发,即台风引起波浪入射至海岸反射并与后续来波相互作用形成驻波作用于海底而激发;(2)台风中心附近源区激发,即台风中心移动过程中不同时期激发的同频率波浪相向传播、相互作用产生驻波作用于海底而激发,源区位置主要集中于台风中心左后方.此外,结合波浪再分析数据、台风风场特征,我们进一步对第二类地脉动激发过程中的影响因素进行了分析,发现:第一种模式激发的地脉动与近岸源区波浪场强度、观测点至源区距离及台风中心至海岸线距离等因素相关;而第二种模式激发的地脉动则主要受台风中心附近波浪场的频率成分与传播方向影响.
2017 Vol. 60 (1): 187-197 [摘要] ( 194 )  HTMLNew  PDF (4770 KB)   ( 1461 )
198 李红玉, 于湘伟, 章文波
场地效应研究的传递台站线性反演法——以芦山地震为例
采用线性反演法研究了地震动的场地效应,并针对线性反演法的局限性,提出了“传递台站”的方法来扩展其应用.通过传递台站,地震事件不再局限于必须被参考台记录到,台站也不再局限于须与参考台站同时获取到地震记录,改进后可大大增加数据样本,使得结果更趋稳定.同时,对于原来难以参与分析的台站,也可得到其场地反应.研究中,还采用了基于遗传算法的非参考台方法,进一步分析了场地反应和验证所提方法的合理性.将所提方法应用于2013年芦山7.0级地震周边区域,分析了这一地区的强震台站的场地反应和地震波衰减因子Qs值.研究结果显示,当采用一个和二个传递台站进行分析时,所能利用的地震事件和所能求解场地反应的台站数目逐步增加,随着数据样本的增加,分析结果趋于稳定.研究区域的Qs值拟合结果为Qsf)=334.16f0.73.遗传算法结果表明,线性反演法中所选取的参考台站有其自身的场地反应.当利用遗传算法获得的参考场的场地反应去校正线性反演法的结果时,对于多数场地,两者的结果符合得很好,说明了遗传算法的结果可以合理地认为是场地较为真实的反应.
2017 Vol. 60 (1): 198-211 [摘要] ( 170 )  HTMLNew  PDF (6482 KB)   ( 438 )
212 袁杰, 朱守彪
断层阶区对产生超剪切地震破裂的促进作用
地震时若断层发生超剪切破裂,地震灾害会显著加剧.因此研究超剪切破裂的形成机理有着非常重要的科学意义.本文利用动力有限单元方法,模拟断层破裂从初始断层跳跃传播到另一条平行的次级断层(断层阶区)时破裂速度的变化情况,并分析断层阶区几何特征等物理参数对产生超剪切地震破裂的促进作用.计算结果表明,断层阶区的诸多物理因素(如:重叠长度、相隔距离以及摩擦系数等)都会对破裂的传播速度产生影响.在一定条件下,破裂传播速度会由在初始断层上的亚剪切波速度,转换为在次级断层上的超剪切波速度.在破裂速度转换过程中,断层间隔起着重要作用,当断层阶区中两断层垂直间隔距离小到一定程度时,破裂跳跃阶区后,破裂速度不会发生变化.所以对于分段断层(可视为一种特殊的断层阶区),由于其断层垂直间隔为0,也就不会出现破裂速度变化的现象,模拟结果对此也进行了证实.然而,若断层间隔太大,当其距离超过一定的限度后,破裂通常无法跨越断层阶区继续传播,而是终止在初始断层上.模拟结果还表明,初始断层与次级断层之间的重叠距离也十分重要,只有当断层阶区中两平行断层之间的重叠部分达到一定长度后,断层的破裂速度才有可能发生转换.此外,计算结果显示,破裂过程中断层面上的应力变化可能是破裂传播速度发生转换的直接原因.最后,模拟还发现,当破裂跨越断层阶区发生速度转换时,破裂需要停顿一定的时间,以便积聚足够的能量来实现破裂速度的增快.
2017 Vol. 60 (1): 212-224 [摘要] ( 176 )  HTMLNew  PDF (4414 KB)   ( 462 )
225 张蕾, 孙知明, 李海兵, 赵来时, 曹勇, 叶小舟, 王雷振, 王焕, 何祥丽, 韩帅, 白明坤, 葛成隆, 赵越
龙门山构造带WFSD-2钻孔岩心磁化率特征及其对大地震活动的响应

岩石磁化率特征可以帮助判断岩石的形成环境,对地震过程中滑动摩擦伴随高温的物理-化学变化具有显著反应.本研究以钻穿龙门山中段构造带的汶川地震断裂科学钻探2号孔(WFSD-2)岩心为研究对象,使用Bartington MS2K磁化率仪对500~2283.56 m深度的岩心进行高分辨率无损磁化率测试,并结合岩性特征和显微结构探讨了龙门山构造带主要岩石单元的磁化率特征及其地震断裂活动的磁学响应.磁化率测试结果表明,由花岗岩和火山碎屑岩组成的彭灌杂岩体的磁化率值(数十到数千个10-6SI)普遍高于上三叠统须家河组沉积岩的磁化率值(数个到数十个10-6SI).从WFSD-2岩性分布来看,彭灌杂岩上下出露四段,其磁化率值特征反映它们属于不同的岩石单元,它们与下伏须家河组地层呈断层接触,构成叠瓦状构造,指示了龙门山构造带具有强烈的地壳缩短作用.断裂带中处于滑动带的断层泥和假玄武玻璃具有高磁化率特征,而断层角砾岩和碎裂岩不具有高磁化率值特征,表明断层岩磁化率增高的原因可能主要与地震断裂滑动摩擦过程中高温作用下发生的磁性矿物转换有关,断层岩中高磁化率异常可作为大地震活动的证据.WFSD-2岩心中的映秀-北川断裂带(600~960 m)可识别出约80条高磁化率异常的断层岩带,揭示映秀-北川断裂带是一条长期活动的断裂带,龙门山构造带形成演化过程中伴随着大地震活动.

2017 Vol. 60 (1): 225-239 [摘要] ( 229 )  HTMLNew  PDF (6728 KB)   ( 1138 )
地球物理学报
   作者服务中心
关于在学术论文署名中常见问题或错误的诚信提醒
我要投稿
投稿须知
论文模板
版权转让协议
标准地图底图下载

   最新地球物理学图书
刘光鼎选集

 
 网络资源
相关地球物理学期刊
相关数据库
相关机构
相关工具软件
版权所有 © 2010 《地球物理学报》编辑部
通信地址:北京9825信箱《地球物理学报》编辑部 (100029)
电 话:010-82998105,82998113
Email:actageop@mail.igcas.ac.cn; geophy@163bj.com
本系统由北京玛格泰克科技发展有限公司设计开发  技术支持:support@magtech.com.cn
京ICP备13017565号-7