地球物理学报杂志 北大期刊 CSCD期刊 统计源期刊

地球物理学报 2016年第12期杂志 文档列表

《地球物理学报》稿约简27-

电磁测深数据地形影响的快速校正简4408-4413

摘要：地形起伏会对电磁法的数据产生一定影响,尤其会影响浅部地层电性结果的准确性.本文通过对地形影响基本规律、经典比值校正原理的分析,认为可以用实测标准电阻率替代数值模拟中均匀半空间电阻率,提出一种新的地形影响快速校正方法,即采用小极矩直流电阻率法获得无地形影响的表层电阻率值,作为地形校正的标准电阻率,以此构造一个新的校正公式.分别对CSAMT和TEM仿真和实际测量资料进行地形校正处理,效果较好.说明新的比值校正公式,是一种快速、有效和实用的校正算法.

伪随机编码源激发下的时域电磁信号合成简4414-4423

摘要：将伪随机编码技术引入到人工源电磁法后,可以通过加大发射功率以及应用后续的相关处理技术来达到压制噪声、加大探测深度及提高分辨率的目的,因此引起了越来越多学者的关注及研究,但大多数研究集中在资料处理的相关技术上,对模拟电磁信号关注较少.然而,资料处理工作大多是从电磁信号出发的,模拟伪随机编码源激发下的电磁信号不但可以为资料处理环节提供理论数据,而且可以为检测资料处理的效果提供中间结果,因此,模拟电磁信号工作必不可少.本文根据获得接收信号的物理过程来实现伪随机编码源激发下的电磁信号合成.首先用解析公式获得特定地电结构的大地频率域响应,然后通过余弦变换得到时间域阶跃响应,接下来用阶跃响应的时间导数得到大地脉冲响应,通过将大地脉冲响应与伪随机编码源的褶积得到理想接收信号,最后,用低通滤波器来模拟发射设备和接收设备的频带限制,将之和噪声一起加到理想接收信号上,最终模拟出仿真的合成信号.通过和野外实际接收信号对比发现本文合成信号仿真度较高,可以服务于后续的数据处理环节.

多道瞬变电磁法(MTEM)虚拟波场提取技术简4424-4431

摘要：多道瞬变电磁法(MTEM)是近年来发展起来的一种新的地球物理勘查技术,其数据采集方式与地震法类似,因此,采用瞬变电磁拟地震解释方法对MTEM数据进行处理解释具有一定的优越性.研究MTEM虚拟波场有效提取方法具有重要意义.在以往奇异值分解法、预条件正则化共轭梯度法两种不同的常规提取方法研究基础上,本文提出采用相关叠加法提取MTEM虚拟波场信息.首先采用预条件正则化共轭梯度法对全时段MTEM数据进行虚拟波场提取,然后采用同种方法对划分的各时间段数据进行虚拟波场提取,最后对全时段提取结果与各个时间段提取结果进行相关性叠加,叠加结果作为最终的提取结果.实测MTEM数据虚拟波场提取结果表明,采用相关叠加法可以得到稳定、光滑的虚拟波场波形曲线,且抗干扰能力强.

高密度电阻率法在海底金矿含水构造探测中的应用简4432-4438

摘要：海底金矿上覆岩层中含水直接对矿床的开采构成威胁,选择山东三山岛金矿新立矿区海底-135m水平沿脉巷道,采用地球物理高密度电阻率法对该巷道635m测线以下200m深度范围内岩体中的含水构造进行了坑道探测,利用三种装置(温纳装置(Wenner)、偶极装置(Dipole-Dipole)、Schlumberger装置)进行实施,相互验证,取得了很好的探测结果.温纳、偶极和Schlumberger三种装置的视电阻率反演结果反映在测线(-135m沿脉巷道)以下至30m深度为一高阻层,表明自-135~-165m的岩体中已不含水或含少量的水;测线以下30~60m段为低阻层,反映-165~-200m的岩体中含基岩裂隙水;测线以下大于60m的地段为特高阻层,反映-200m以下岩体含水性逐渐变差.该探测结果与矿区水文地质结构调查分析结论具有很好的一致性,表明了坑道高密度电阻率法探测含水构造是可行的.

多道瞬变电磁法共中心点道集数据联合反演简4439-4447

摘要：多道瞬变电磁法是目前地面资源勘查领域的研究热点,开展该方法观测数据的一维反演对完成装备研发及技术推广具有积极的推动作用.针对多道瞬变电磁法特殊的源波形和数据采集方式,本文对多道瞬变电磁数据的响应提取和一维反演进行了研究.首先,采用反卷积方法从伪随机二进制序列全波形响应数据中提取大地脉冲响应.为了避免脉冲响应中空气波的干扰,提出采用反向积分算法,把脉冲响应转换成下降沿阶跃响应,并将其整理成共中点道集.最后,为充分利用各个偏移距下的响应数据对不同深度范围目标体的分辨能力,采用基于OCCAM算法的联合反演方法对共中心点道集数据进行联合反演.结果表明在不同目标体深度以及有噪声干扰的情况下,共中心点道集数据联合反演计算均可获得较好的结果.

瞬变电磁探测系统(CASTEM)试验对比——安徽颖上大王庄铁矿简4448-4456

摘要：为推进国产电磁勘探设备的高端化进程,自主研发了CASTEM瞬变电磁系统.在前期系统开发、测试及集成基础上,开展野外试验研究,以检验CASTEM系统的稳定性及探测能力.选择已取得显著勘探成果的安徽颖上县大王庄铁矿进行试验,试验内容包括大功率发射机性能、接收系统噪声性能、实测响应曲线和反演效果,并与国际高端瞬变电磁仪在上述几方面进行了对比.试验结果表明CASTEM系统的各项指标性能稳定、可靠,已达到国际先进水平;探测成果与高端瞬变电磁仪器探测成果以及钻孔资料相吻合,表明CASTEM系统可以用于深部勘探.

磁性层状介质的CSAMT电场分量响应特征简4457-4463

摘要：在传统的勘探电磁学理论中,往往把岩层磁导率近似看成空气中的磁导率,从而得到简化后的电磁场响应式,主要用于指导非磁性介质情况下的电磁探测.本文开展了磁性介质的电磁响应理论、数值模拟和特征分析等研究,以便实现对磁性矿体的电磁精细探测.首先推导出磁性层情况下的电磁响应表达式;然后,计算了几种典型磁性介质地质断面的响应曲线,并对其响应特征进行了分析;研究结果表明:当磁性层位于地表时,当表层相对磁导率μr>1.2,磁性层对电磁测深曲线具体明显的影响;当岩层为中间层时,只要μr>1.2,对H型、A型曲线中段会有影响.对H型曲线,主要表现在使极小值处的曲线抬高、变宽,以致使曲线形态发生变化;当磁性层处于最后一层时,不论岩层的磁导率是多少,对各类二、三层曲线来讲,不会影响曲线的形态,只是使渐近线前第一个极值点的位置右移,极值点至渐近线之间的线段变得较为陡峭.

基于交叉梯度约束的地震初至纵波与瑞雷面波联合反演简4464-4472

摘要：本文提出了一种初至纵波(P波)与瑞雷面波的交叉梯度联合反演策略.通过对初至P波进行全波形反演可以获得近地表P波速度结构;通过对仅含瑞雷面波信息的地震数据转换到频率-波数域进行加窗振幅波形反演(Windowed-Amplitude Waveform Inversion,w-AWI)可获得近地表横波(S波)速度结构.在二者反演的目标函数中均加入P波速度和S波速度的交叉梯度作为正则化约束项,使得在反演过程中P波速度和S波速度相互制约,相互约束,从而实现对地震初至P波与瑞雷面波的联合反演.数值模拟结果表明交叉梯度联合反演可以提高S波速度反演分辨率,而P波速度反演结果并没有得到提高.实际资料的反演结果表明,交叉梯度联合反演能够获得更加可信的近地表速度结构.

不依赖源子波的跨孔雷达时间域波形反演简4473-4482

摘要：波形反演是近年来较热门的反演方法,其分辨率可以达到亚波长级别.在波形反演的实际应用中,源子波的估计十分重要.传统方法使用反褶积来估计源子波并随着反演过程更新,该方法在合成数据波形反演中效果较好,但在实际数据反演过程中存在一系列的问题.由于实际数据信噪比较低,在源子波估计过程中需要大量的人为干涉,且结果并不一定可靠.本文使用一种基于褶积波场的新型目标函数,令反演过程不再依赖源子波.详细推导了针对跨孔雷达波形反演的梯度及步长公式,实现介电常数和电导率的同步反演.针对一个合成数据模型同时反演介电常数和电导率,结果表明该方法能够反演出亚波长尺寸异常体的形状和位置.接着,将该方法应用到两组实际数据中,并与基于估计源子波的时间域波形反演结果进行比较.结果表明不依赖源子波的时间域波形反演结果分辨率更高,也更准确.

起伏地形下各向异性的2D大地电磁无网格法数值模拟简4483-4493

摘要：作为一种基于节点的计算方法,无网格法具有构造高阶导数方便,自适应分析便利的优点,特别适合复杂地质构造的数值模拟.本文针对实际地球物理勘探中存在的起伏地形和各向异性的地电结构,提出用无网格法来模拟大地电磁响应,采用复合二次径向基函数构造了形函数,推导了大地电磁无网格法等价线性方程组,研究了系数矩阵的压缩存储方法以及大型稀疏复线性方程组快速求解算法,实现了起伏地形下各向异性的2D大地电磁高精度数值模拟.基于层状模型验证了算法的正确性,计算结果表明:无论是TM模式还是TE模式,计算相对误差均小于1%;通过对地垒和地堑模型的模拟,得出了起伏地形对视电阻率和相位的影响规律;对起伏地形下含有不同各向异性系数异常体的模型进行了数值模拟,为开展复杂地质结构的电磁场特征研究以及地形校正奠定了理论基础.随着计算科学的快速发展,无网格法必将成为新的高精度电磁场数值模拟方法.

富含煤层气储层的叠前地震反演预测方法研究——以山西晋东南测区为例简4494-4504

摘要：叠前地震反演技术能够将地震中振幅、偏移距和入射角等多种有效信息与油气的敏感参数结合起来,在有效识别油气藏“甜点”方面虽发挥了重要的作用,但对于煤层气的有效预测仍需研究.针对煤层气与石油天然气的赋存地质条件不同,本文以山西沁水盆地煤层气地震资料为例,借鉴叠前地震反演技术的多参数预测思想,通过模型正演和实测资料处理,得到纵波、横波速度和密度等不同的数据体,利用数据体和弹性参数关系式,采用Aki&Richard近似公式与叠前宽角度反演联合的方法,寻找出富含煤层气储层的敏感参数.从而实现富含煤层气储层的有效预测.经与测井曲线结果对比,吻合率较高,说明了该方法具有可行性与实用性.

深海热液硫化物矿体3D瞬变电磁正演简4505-4512

摘要：深海热液硫化物矿体瞬变电磁的正演是考虑深海环境的全空间条件下三维体的涡流电磁响应.采用全空间矢量有限元法模拟计算深海热液硫化物矿的三维瞬变电磁响应,对硫化物矿体采用矩形单元模型剖分,应用Galerkin法推导有限元方程,先计算频率域响应,再通过Fourier反变换将其转换至时间域,得出深海热液硫化物矿矿体的瞬变电磁响应.并用双半空间模型的解析解检验了全空间矢量有限元法模拟计算算法和程序的正确性,最后按照等比例缩小电磁物理实验原则,比对数值计算和物理实验结果论证了全空间3D模型数值的正确性.结果表明:对于海水、矿体以及围岩复杂电磁边界,应用全空间矢量有限元法模拟计算深海热液硫化物矿瞬变电磁响应异常与物理模拟结果一致,而且计算方法简单精确,异常幅值明显,边界清晰.

煤矿井下微震震源高精度定位研究简4513-4520

摘要：煤矿井下微震震源准确定位,对于动力灾害监测预警具有重要意义.由于微震震源需要通过井下传感器接收信息反演确定,传感器的安装位置限制于煤矿井下巷道周围,传感器沿巷道面的不合理布置将大大降低震源定位精度.针对由传感器信息反演震源位置引起的病态问题,本文提出了基于微震监测测点优化布置的震源高精度定位算法.首先通过计算系数矩阵条件数,判定病态问题;然后利用中心化法和行平衡法联合进行病态矩阵预处理.对预处理后的矩阵A、b利用L曲线法计算正则参数,结合Tikhonov正则化算法计算得到震源坐标正则解.研究结果表明,中心化法有效降低了矩阵数量级,行平衡预处理降低了病态条件数,预处理后Tikhonov正则解的震源坐标误差最小可以达到3.09m,与预处理前的高斯消去解相比误差大大降低.通过上述优化处理,实现了井下受限空间微震监测震源高精度定位.

一种刻槽式随钻声波测井隔声体的理论与实验研究简4521-4528

摘要：隔声体的设计与效果评价对随钻声波测井仪器的研发至关重要.针对刻槽式隔声体结构,利用数值模拟的方法分析了钻铤波在其传播过程中隔声阻带形成、拓宽并达到稳定的物理过程.在此基础上优化设计并加工了一种可用于随钻声波测井实验样机的刻槽式隔声体.在实验室内利用声波频率扫描方法测量了该隔声体的钻铤波声衰减,并进一步在实验井中测试验证了隔声体在井孔中工作时的隔声效果.最后,开展了随钻测井环境下的现场实验,测到了高质量的地层波信号和对应的纵波时差.理论分析和实验对比研究结果表明,刻槽隔声体有效地降低了钻铤直达波的干扰.本文的工作为随钻声波测井刻槽式隔声体的设计及评价测试提供了一套系统的理论及实验方法.

2015年尼泊尔M_S8.1地震震源区S波三维速度结构与强震发生机理研究简4529-4543

摘要：利用2001—2003年期间在2015年4月12日尼泊尔MS8.1级强震震源区流动地震观测记录到的连续波形数据,提取了5~25s周期的瑞利波相速度频散曲线,并构建了尼泊尔地震震源区二维瑞利波相速度分布图像.以0.5°×0.5°为网格大小将研究区网格化,采用NA算法反演得到尼泊尔地震震源地区三维S波速度结构.结果显示,在上地壳,以主前锋逆冲断裂带(MFT)为界,其以北地区为高波速异常,而其以南为明显低波速异常;在中地壳,以藏南拆离系(STDS)为界,南北两侧速度结构也存在明显差别,以南地区为明显高波速异常,而以北地区为明显低波速异常.这些结构特征说明,印度板块与欧亚板块碰撞挤压作用形成地幔热物质上涌并造成地壳物质部分熔融,并由此形成了东西向拉张的南北向裂谷.2015年尼泊尔MS8.1级主震和最大余震均发生于高低波速异常过渡区且偏向高波速异常区,暗示了这样的波速异常区易于积累能量孕育强震.主震和最大余震的南侧均存在明显的低波速异常,与主喜马拉雅滑脱断裂带(MHT)相对应,可能代表部分熔融或深部流体作用于主边界断裂带(MBT)附近的MHT断裂带,降低断层面上的有效正应力,从而触发尼泊尔强震及最大余震的发生.主震与最大余震之间的余震分布于高低波速异常变化较为明显的地区,说明研究区内地震的发生受震源区附近的速度结构控制.

地壳分层各向异性介质接收函数及其粒子群反演简4544-4559

摘要：地壳不同深度介质的地震各向异性是研究地壳不同深度范围变形方式的重要依据.鉴于地壳介质的复杂性,如何从远震体波接收函数中提取不同深度的各向异性参数仍是一个有待深入研究的课题.在已有研究的基础上,本文利用广义反射-透射系数矩阵方法计算的合成地震图,研究了复杂地壳分层各向异性介质的接收函数随反方位角(back azimuth)变化及不同层位各向异性参数对接收函数波场的影响,为各向异性介质接收函数的解释提供了新的理论依据.通过引入粒子群优化理论,发展了分层各向异性介质接收函数全局反演算法.数值及观测数据的验证结果表明,在各向同性速度模型确定的前提下,我们的方法能够可靠地提取地壳分层各向异性参数;在反演中引入曲波变换去噪技术,对于正确解析不同层位的各向异性参数具有重要价值.

柴达木盆地东部新生代盆地结构与演化简4560-4569

摘要：本文对柴达木盆地东部新生代盆地结构和构造演化进行了研究.地震剖面揭示柴达木盆地东部新生界分别向南、北盆地边缘变薄和尖灭,盆地北部被欧龙布鲁克山和埃姆尼克山隆起强烈改造.通过对新生代地层厚度横向变化以及地层剖面分析,确定欧北断裂自中新世晚期开始向北逆冲,导致欧龙布鲁克山发生隆升和德令哈凹陷的形成.埃北断裂从上新世开始活动,与欧北断裂同时向北逆冲,导致德令哈凹陷进一步沉降,形成厚度达2600m的狮子沟组.埃南断裂在第四纪开始大规模向南逆冲,不仅造成其北侧的欧龙布鲁克山和埃姆尼克山隆起强烈抬升和向南推覆,而且导致南侧霍布逊凹陷的形成,成为柴达木盆地第四纪沉积中心.早期提出的前陆盆地和背驮式盆地模型显然不能解释柴达木盆地东部新生代构造格架和演化历史.本次研究认为柴达木盆地东部的形成是强烈的水平挤压作用导致地壳发生大规模褶皱的结果,即柴达木盆地东部新生代是一个大规模向斜.该向斜盆地模型很好地解释了新生代地层向盆地边缘减薄以及沉积中心主要位于盆地中部等现象.了解柴达木盆地东部构造发展对了解青藏高原侧向扩展具有重要意义.