地球物理学报杂志社
分享到:

地球物理学报杂志

《地球物理学报》杂志在全国影响力巨大,创刊于1948年,公开发行的月刊杂志。创刊以来,办刊质量和水平不断提高,主要栏目设置有:学术论文、综述或评述、研究短论、科学论坛、学术动态等。
  • 主管单位:中国科学院
  • 主办单位:中国地球物理学会;中国科学院地质与地球物理研究所
  • 国际刊号:0001-5733
  • 国内刊号:11-2074/P
  • 出版地方:北京
  • 邮发代号:2-571
  • 创刊时间:1948
  • 发行周期:月刊
  • 期刊开本:B5
  • 复合影响因子:1.73
  • 综合影响因子:1.867
相关期刊
服务介绍

地球物理学报 2015年第06期杂志 文档列表

地球物理学报杂志空间物理学·大气物理学

顶部电离层离子密度经度结构的特征及其随季节、太阳活动和倾角的变化

摘要:本文利用DMSP卫星测量数据和傅里叶分解和重构方法,研究了地磁平静期顶部电离层总离子密度(Ni)经度结构的多重波数特征及波数4的年变化、逐年变化、地方时差异和随倾角的变化.傅里叶分解和重构的结果表明,顶部电离层平均Ni的经度结构中同时含有以波数1至波数4为主的多重分量,不同波数分量的幅度和相位各不相同.对波数4分量的分析表明,波数4的幅度在春秋季最强,北半球夏季高于冬季;随太阳活动水平增强,波数4分量的幅度增高,至太阳活动高年幅度达到最高,此后随太阳活动水平降低而减小,与F10.7呈正相关;春秋季和北半球夏季波数4分量在傍晚最强,晚上和上午次之,黎明最弱,从09LT到21LT,波数4的相位依次滞后,暗示向东移动.分析还发现,日落期间波数4幅度依赖倾角,春秋季随倾角的变化呈双峰结构,两个极大出现在倾角±18°附近,暗示赤道等离子体喷泉效应对顶部电离层经度结构的控制作用.
1843-1852

大功率电波加热电离层中热自聚焦不稳定性的理论研究和数值模拟

摘要:本文首先从电子密度及电子温度的输运方程和考虑自作用时的电磁波波动方程出发,利用简正模展开的方法推导出泵波在反射区域激发出热自聚焦不稳定性(thermal self-focusing instabilities,TSFI)所需电场阈值以及其增长率的完整数学表达式,并估算了TSFI激发阈值及所对应的有效辐射功率(ERP)的量级.随后利用三维垂直加热的理论模型,结合国际参考电离层(IRI-2012)和中性大气模型(MSIS-E-00)给出的背景参数,数值模拟了大功率高频泵波加热电离层时泵波反射区域电子密度及电子温度因TSFI而产生的变化及发展的过程,并对比分析了不同背景参数对较热效果的影响.结果表明:当高频泵波的加热阈值达到或超过百毫伏每米的量级时,即可激发TSFI,发展出大尺度电子密度及温度不均匀体,这些不均匀体内的密度耗空约为4%~10%,而电子温度剧烈增长,到达背景温度值的1.6~2.1倍;且在相当的加热条件下,背景电子温度越低、电子密度越小,加热效果越显著;电子密度及电子温度的扰动幅度随着加热时间的推移而逐渐减小,即扰动逐渐趋于饱和,且电子温度要快于电子密度达到饱和状态.本文还对泵波反射高度处的电子密度及电子温度变化率进行采样并求得其功率谱密度,分析结果表明:TSFI发展出的大尺度不均匀体满足幂律谱结构,谱指数随着加热的进行逐渐趋于稳定,白天与夜间的幂律谱指数区别不大,但电子密度与电子温度的幂律谱有所区别.
1853-1868
地球物理学报杂志地震学·海洋地球物理学

利用视震源时间函数反演尼泊尔M_S8.1地震破裂过程

摘要:以2015年4月26日MS7.1余震为经验格林函数事件,利用全国和全球的宽频带记录提取了2015年4月25日尼泊尔MS8.1地震的P波视震源时间函数和Rayleigh波视震源时间函数,并通过联合反演这些视震源时间函数获得了这次地震的时空破裂过程图像.无论是P波视震源时间函数还是Rayleigh波视震源时间函数都呈现出很强的方位依赖性,表明震源断层具有相当的尺度且破裂朝东南方向扩展.时空破裂过程图像清楚地证实了这一特征,并更清晰地显示,破裂几乎是纯粹的单侧破裂,从破裂起始点开始,沿断层面向东南方向扩展~100km,同时沿断层面向深部扩展~80km,形成~125°的破裂优势方向和~5.8m的最大位错.地震的破裂时间历史相对简单,呈非间断性扩展,持续时间约50s.
1881-1890

利用多台阵压缩传感方法反演尼泊尔M_W7.9地震破裂过程

摘要:2015年4月25日,尼泊尔地区发生MW7.9地震,震中位于28.1°N,84.7°E.为了详细地研究此次破坏性极强的地震的破裂过程,本文利用多台阵压缩传感方法,使用了阿拉斯加、欧洲和澳大利亚三个台网的共计179个台站的远场P波垂直分量的数据来反演,结果表明本次地震的破裂过程是一个清晰的南东东方向的单侧破裂,破裂尺度约为105km,整体持续时间约为58s.在破裂初始的前15s,能量辐射基本围绕在震源附近,16s后破裂开始向南东东方向以1.9km·s-1的速度破裂.释放能量最大的时间为第38s,位于距震中70km处.该位置从第29秒开始破裂,并持续释放能量长达30s之久.
1891-1899

2015年尼泊尔M_S8.1地震的地壳重力均衡背景与地表形变响应特征

摘要:对2015年尼泊尔MS8.1地震的地壳均衡背景及其引起的地表形变特征进行了研究,结果表明:(1)尼泊尔MS8.1地震震中以南的印度板块岩石圈有效弹性厚度大约为9km,加载主要来自地幔;地震以北的拉萨地块岩石圈有效弹性厚度大约为2km,加载主要来自地表.(2)尼泊尔MS8.1地震震中以南地区的地壳均衡异常大约为-100mGal(10-5 m·s-2),但其北部的地壳均衡异常则为300~400mGal,尼泊尔MS8.1地震发生在地壳均衡负异常向正异常过渡的高梯度带上.(3)尼泊尔MS8.1地震使震中周围地区的地壳整体向南运动,最大水平位移超过1.5m,分布在震中东南.震中以北的同震垂向位移总体为负值,最大下降幅度超过0.5m,同震重力变化总体为正值,最大超过60μGal(10-8 m·s-2);震中以南的垂向位移总体为正值,最大升幅超过0.7m,同震重力变化总体为负值,最大降幅超过-120μGal.(4)尼泊尔MS8.1地震使"世界屋脊"喜马拉雅山脉产生沉降,最大同震降幅超过120mm,震后松弛效应将使"世界屋脊"持续缓慢下降.该强震使世界最高峰珠穆朗玛峰降低了2~3mm,有可能被GPS、InSAR等现代大地测量工具检测到.
1900-1908

2013年7月22日甘肃岷县-漳县M_S6.6地震震源破裂过程

摘要:2013年7月22日,在甘肃岷县漳县交界处发生MS6.6地震,地震震中位置靠近临潭—宕昌断裂.本文通过构建有限断层模型,利用国家强震动台网中心提供的12条强地面运动三分量资料,通过波形反演方法来研究这次地震的震源破裂过程.结果显示这次地震是发生在甘东南地区岷县—宕昌断裂带东段附近的一次MW6.1级逆冲兼具左旋走滑破裂事件,最大滑动量约为80cm.发震断层走向及滑动性质与岷县—宕昌断裂吻合,推断本次地震与东昆仑断裂向北的扩展和推挤密切相关,是岷县—宕昌断裂进一步活动的结果.
1909-1918

基于经验模型和物理模型研究2013 M_S7.0芦山地震余震序列

摘要:基于中国地震台网中心2013 MS7.0芦山地震余震数据我们首先确定了余震空间分布范围并根据G-R关系计算了主震后半小时内的完备震级Mc=3.5,并且得到了ML≥3.5和ML≥3.0的地震在2001年至芦山地震前的背景场地震发生率.通过Omori-Ustu经验定律和两种Dieterich模型对芦山地震余震发生率的拟合,我们发现阶梯型Dieterich模型只能模拟p=1的情况,从而造成了模拟曲线与观测数据的差别;前人研究表明震后滑移同样是产生余震的原因,如果假设余震序列由主震静态剪应力Δτ和震后滑移共同作用所产生,我们数值模拟得到的对数型Dieterich模型能够较好地推断余震发生率R随时间t增加而衰减的趋势,能够从物理机制上解释MS7.0芦山地震余震序列衰减指数大于1这一现象.通过对数型Dieterich模型的拟合并结合Andrews的方法,我们还得到MS7.0芦山地震Aσ约为0.155 MPa,ta约为8.4年,这一值与前人研究结果十分接近.
1919-1930

2010年4月玉树M_S7.3地震序列的断层结构

摘要:利用双差定位方法对玉树地震序列2010年4月14日至10月31日间发生的ML≥1.0地震进行双差定位,得到1545个地震的重定位结果.综合分析地震双差定位结果和玉树地震序列中强地震震源机制解,发现玉树MS7.3地震发震构造由北西向和北东东向两条相交断层组成,主震发生在北西走向的甘孜—玉树断裂带上,5月29日的MS5.9余震序列发生在北东东走向的一条隐伏断裂上,两条断裂均接近直立.甘孜—玉树断裂是羌塘地块和巴彦喀拉地块的构造边界,由于羌塘地块和巴颜喀拉地块的差异运动使甘孜—玉树断裂强耦合段应力高度积累,在应变能超过岩石强度时破裂失稳发生了MS7.3地震.主震断层的左旋滑动导致北东东向断层的正应力减小,库伦应力增加,45天后触发了MS5.9余震序列的活动.
1931-1940

2013-2014年山东乳山地震序列发震构造初探

摘要:采用结合波形互相关技术的双差定位方法,对2013—2014年山东乳山地震序列重新定位,通过CAP及P波初动方法确定乳山序列较大地震的震源机制,在此基础上初步探讨乳山地震序列发震构造.结果显示,乳山序列呈现NW向展布,地震密集分布在8km×3km范围,震源深度分布在4~10km,4~7km区间相对集中.较大地震震源机制的节面Ⅰ方向与序列地震优势分布方向基本一致.综合考虑精确定位结果及较大地震震源机制,并结合震区附近地震资料,初步推测乳山地震序列发震断层为NW方向、近直立的走滑型隐伏断裂.
1954-1962

一种新型海底沉积物声学原位测量系统的研制及应用

摘要:为了能够精确地测量海底表层沉积物的声学参数,自主研制了一种新型海底沉积物声学原位测量系统,与国内外传统的声学原位测量系统相比,该系统能够实时显示声波波形,调整测量参数,其工作方式除了站位式测量之外,还实现了拖行式连续测量,极大地提高了工作效率.根据前期海试情况,对海底仪器结构进行了重新设计,使之可以同时测量海底沉积物及海底海水的声学参数,同时建立了双向数字信道,解决了测量过程中系统信号的干扰问题.该系统的结构分为两部分:甲板控制单元和水下测量单元,整套系统通过主机控制程序进行控制,采用GPS定位系统测定仪器的大地坐标.为了检验系统的稳定性及准确性,分别进行了实验室水槽实验和海上试验.利用水声测量设备对测量系统进行实验室水槽标定分析,实验结果表明系统测量值相对误差仅为0.04%,测量结果具有较高的精度.海上试验在青岛胶州湾和东海海域进行,获得了试验区域海底沉积物声速和声衰减系数的测量数据,将测量数据与他人的研究结果进行对比分析,结果表明测量数据与前人研究结果一致,较为准确.该原位测量系统在站位式测量和拖行式测量中都能够快速准确地测量出沉积物声速和声衰减系数,可以作为海底底质声学测量的调查设备.
1976-1984
地球物理学报杂志应用地球物理学

基于快速推进迎风双线性插值法的三维地震波走时计算

摘要:三维地震波走时计算技术是三维地震反演、层析成像、偏移成像等诸多地震数据处理技术中非常重要的正演计算工具.为了获得精度高且兼顾效率的三维走时计算方法:首先,在常规双线性插值公式推导过程中,充分利用平面波双线性假设的结论,获得了二元极小值超越方程的解析解,进而推导出了准确的局部走时计算公式,同时构造性地证明了该计算公式满足地震波的传播规律和Eikonal方程;其次,引入迎风差分的基本思想,提出迎风双线性插值的局部走时计算策略,该计算策略能简化算法、提高效率且保证无条件稳定性;然后,将上述计算公式和迎风双线性插值策略与常规快速推进法中的窄带技术结合,获得了一种新的基于快速推进迎风双线性插值法的三维地震波走时计算方法;最后,通过精度和效率分析检验了新算法的精度、效率和正确性,并通过计算实例验证了算法在面对复杂介质时的稳定性和有效性.
2011-2023

叠前地震数据特征波场分解、偏移成像与层析反演

摘要:本文提出了一套叠前地震数据稀疏表达(特征波场合成)、深度偏移成像和层析成像的处理流程.不同于传统的变换域中的数据稀疏表达理论,本文利用局部平面波的传播方向(慢度矢量),在中心炮检点处同时进行波束合成,从而将地震数据投影到局部平面波域(高维空间)中.由于波束合成后的地震数据描述了局部平面波的方向特征,因此称之为特征波场.然而波束合成算法需要估计局部平面波的慢度矢量.当地震数据受噪声干扰时,难以在常规τ-p谱中自动估计局部平面波的射线参数(慢度矢量).本文提出了基于反演理论的特征波场合成方法,可以同时反演局部平面波及其传播方向,从而提高特征波合成的自动化程度并保持方法的稳健性.通过特征波场合成,可以将地震数据分解为单独的震相(波形).这样的数据可以直接用来成像及反演.在局部平面波域中,由于局部平面波的入射与出射射线参数已知,传统的Kirchhoff叠前深度偏移(PSDM)和高斯束/控制束PSDM可以实现从"沿等时面的画弧"到"向反射点(段)的直接投影"的转变,叠前偏移的效率以及成像质量可以同时提高.此外,特征波场与地下反射点(段)的一对一映射关系使得叠前深度偏移与层析成像融为一体,可以极大地提高速度反演的效率.数值试验证明了特征波场合成、叠前深度成像以及层析反演的有效性.
2024-2034

各向异性TI介质qP反射波走时层析成像

摘要:地震走时层析成像是反演地层各向异性参数分布的有效方法,但是关于地震各向异性介质走时层析成像的研究并不多,其技术远远没有达到成熟的阶段.在野外数据采集时,地表反射波观测方式相对井间和垂直地震剖面观测方式的成本更低,利用qP反射波走时反演各向异性参数具有更加广泛的实用价值.本文实现的TI介质地震走时层析成像方法结合了TI介质反射波射线追踪算法、走时扰动方程和非线性共轭梯度算法,它可以对任意强度的TI介质模型进行反演,文中尝试利用qP反射波走时重建TI介质模型的参数图像.利用qP反射波对层状介质模型和块状异常体模型进行走时反演,由于qP波相速度对弹性模量参数和Thomsen参数的偏微分不同,所以可以分别反演弹性模量参数和Thomsen参数.数值模拟结果表明:利用qP反射波可以反演出TI介质模型的弹性模量参数与Thomsen参数,不同模型的走时迭代反演达到了较好的收敛效果,与各向同性介质走时反演结果相比较,各向异性介质走时反演结果具有较好的识别能力.
2035-2045

一种基于平面波静态编码的最小二乘逆时偏移方法

摘要:平面波偏移是一种面炮偏移方法,相对于常规逐炮偏移,其具有较高的计算效率.然而常规平面波偏移方法成像精度低,且成像时会产生串扰噪音.为此,本文在实现常规平面波偏移算法基础上,引入反演思想实现了基于静态平面波编码的最小二乘偏移理论方法及处理流程,在优化算法基础上对平层模型和复杂砂砾断块模型进行了成像测试并与其他成像策略进行对比.研究结果表明:基于时移编码的平面波最小二乘偏移能有效抑制低频成像噪音和串扰噪音,补偿中深部成像能量,是一种较为有效的保幅成像策略.
2046-2056

基于Duffing振子混沌系统的地震速度分析方法

摘要:强随机噪声干扰是导致地震勘探资料低信噪比的主要原因,如何在强随机噪声干扰下获取有效的信息是值得关注的问题.Duffing振子混沌系统是一个非线性的动力学系统,其对强随机噪声具有免疫能力,而对特定的周期性信号具有敏感性.本文提出一种基于Duffing振子混沌系统的速度分析方法.对CMP道集按照时距曲线关系进行移动窗口截取,将所截取的信号构建为待测信号加入Duffing振子混沌系统,通过相图网格分割方法(GPM)判断系统状态的改变,从而在强随机噪声背景下获得高分辨率的速度谱.理论模型和实际资料的处理结果表明,与传统的水平叠加速度分析方法相比,本方法能够在强随机噪声背景下获得更准确的速度分析结果.
2057-2068

基于Curvelet变换与POCS方法的三维数字岩心重建

摘要:随着页岩气勘探与开发的深入,研究页岩裂隙的三维空间展布成为页岩岩石物理研究的必要步骤之一.但由于仪器的限制,页岩切片在深度上具有不连续性,以及数字岩心纵向上成像最小间隔与横向分辨率的不一致成为影响裂隙表征和数字岩石物理模拟精度提高的重要因素.为了更好的研究裂隙在三维的空间展布,本文将curvelet稀疏变换与凸集投影(POCS)迭代算法有效结合,实现三维数字岩心重建.首先对X射线扫描砂岩得到的三维数据体进行隔片抽稀,利用本文方法实现三维数据体重建,重建结果与完整数据体具有很好的一致性,且优于现有方法(spgl1),验证了新方法的有效性与先进性.其次对聚焦离子束扫描电镜(FIB-SEM)得到的纳米级页岩二维切片在深度上进行了加密重建,获得纵向上成像最小间隔与横向分辨率基本一致的三维数字岩心,由于仪器限制引起的页岩切片深度上的不连续性得到减弱,裂隙展布更加清晰.砂岩CT图像以及页岩FIB-SEM成像数据的重建结果验证了本文方法的有效性与先进性.
2069-2078

储层重力密度反演后验约束正则化方法

摘要:本文针对蒸汽辅助重力泄油(SAGD)生产中开发监测问题,发展了综合应用地震及重力数据反演储层密度的联合反演算法.通过测井数据建立纵波阻抗与密度的直接关系,并推导出这种关系下重力与纵波阻抗数据联合反演的计算方法,从而计算出蒸汽腔体密度分布规律.文中应用密度反演后验约束正则化方法,采用Tikhonov正则化模型,通过波阻抗数据作为约束进行联合反演,在算法上提高了稳定性,同时得到较高的反演精度.文中对SAGD生产中的理论模型进行了方法试算,并分析了算法的误差,最终应用于SAGD生产的实际数据中,通过最终反演结果分析,该方法取得了很好的应用效果.
2079-2086

3D DEMAP数值模拟与观测实验

摘要:为推动电磁剖面(EMAP)技术向找矿勘探领域发展,本文针对密集阵列剖面(DEMAP)测量方式,利用三维积分方程法对层状介质中赋存三维地质异常体的电性结构进行了数值模拟,并在青海野马泉铁锌矿集区西部进行了大地电磁(MT)与DEMAP的野外对比观测实验.数值模拟结果显示,DEMAP观测方式获得的结果受偏移距影响,造成的视电阻率和相位误差与偏移距大小和地下电性结构的复杂性密切相关,但在整体上,统计误差≤10%.通过对野马泉矿集区的实验结果分析可得,DEMAP与MT的视电阻率断面和相位断面具有较好的一致性.本文从理论和实验上都表明DEMAP观测技术是一种有潜力的、经济、快速、有效的找矿手段.
2103-2111