波场
- 基于单程波算子的全波场最小二乘偏移方法
地震成像应该是全波场（一次波、表面多次波和层间多次波）的贡献。为获得全面准确的地下信息，全波场成像是地震资料处理技术发展的必然趋势。野外采集的陆地资料相对复杂，一般包含直达波、面波、折射波、绕射波、一次波和多次波等，所以陆地资料的全波场成像非常困难。海洋地震资料不受近地表的影响，具有较高的信噪比，发育的波场主要是一次波和多次波，因此就全波场成像研究而言，海洋地震资料具有优越的数据基础。对于海洋深水地震资料，根据实现方法的不同，利用多次波成像方法主要有四类。
石油地球物理勘探 2023年5期2023-11-26
- 一次波和层间多次波联合成像方法
引入基于上、下行波场分离的成像条件提高成像精度；然后，针对联合成像过程中引入的串扰假象，结合逆时反传过程分析了串扰产生的根源，提出了基于缺失边界匹配补偿的解决策略，实现了波场的保真逆时延拓，有效避免了串扰假象的产生；最后，通过数值试算验证了方法的可行性、有效性和适用性。1 一次波和层间多次波联合成像方法原理与单程波偏移类似，逆时偏移也是先进行单炮偏移，然后通过多炮叠加实现最终成像。图1a展示了一次波逆时偏移过程，具体可描述为如下三步：①正演计算，地表激发地
石油地球物理勘探 2022年6期2022-12-09
- 三维TTI介质中的弹性波场分解
，重建P波和S波波场，然后运用成像条件分别得到PP波、PS波、SP波和SS波成像结果[5-6]。但是，地震波在实际地层中传播时，耦合的P波和S波波场会产生串扰噪声，严重降低成像剖面的质量，所以在应用成像条件之前，需要准确分离P波和S波波场。地面多分量地震记录可通过矢量合成法计算P波、S波的矢量方向[7]，把地震波振幅旋转到真正的纵、横波矢量轴，实现P波、S波波场分离[8-10]。另外，也可以在逆时偏移波场延拓过程中进行波场分离，获得纯P波、纯S波波场。基于
石油地球物理勘探 2022年6期2022-12-09
- 一种基于新差分模板和无穷范数的震源波场重建方法
随状态法通过震源波场和伴随波场间的相互运算(如相关)计算参数梯度/敏感核。但震源波场和伴随波场分别沿时间正向和逆向传播。为同时访问震源波场和伴随波场，震源波场需先保存在计算机内存中或在波场反向传播中进行重建。存储所有时刻、所有空间网格点的震源波场是不现实的，特别是三维情况。最优检测点方法(Optimal CheckPointing)在波场正向传播中仅需保存几个时刻(检测点)的震源波场，在波场反向传播中将存储的检测点波场作为初值重新计算相邻检测点之间的震源波
石油地球物理勘探 2022年6期2022-12-09
- 面向隧道仰拱的高频地震波逆时偏移成像
的传播规律，通过波场正传、波场重构与波场反传计算出准确的震源波场与检波点波场，然后对同一时刻的两个波场值，利用互相关成像条件提取该时刻的成像值，最后对所有时刻的成像值进行叠加，即可得到完整的成像剖面[5]。在二维各向同性声波介质中，一阶速度—应力方程可表示为如下形式：式中，Δx为x 方向网格间距，tΔ 为时间步长，为有限差分系数。式（2）即为震源波场正向时间延拓公式。检波点波场反向延拓就是正演的逆过程，其延拓公式与正向延拓公式类似：根据成像条件，震源波场与
中国水运 2022年10期2022-11-17
- 结合泊松算法和波场分解成像条件的各向异性优化纯声波方程逆时偏移
单、计算量较小、波场特征明确。目前在各向异性波场模拟和偏移中主要采用伪声波方程，该类方程的一种推导思路是将耦合的各向异性P-SV波频散关系中的横波速度设置为0简化而来[1]。在此基础上，发展了多种形式的伪声波方程[2-8]。另一种伪声波方程的推导思路是基于声学近似的胡克定律[9-10]。但是当给定的各向异性参数不满足伪声波假设时，现有伪声波方程在进行P波模拟时，容易出现伪横波干扰及数值不稳定。为了彻底解决伪声波存在的问题，Liu等[11]解耦了VTI介质P
石油地球物理勘探 2022年3期2022-06-11
- 隧道瞬变电磁克希霍夫偏移成像与地质灾害探测
精度的有效手段。波场变换为瞬变电磁拟波动场成像提供了前提基础。瞬变电磁扩散场方程与波动方程存在数学变换关系[10-12]，将满足扩散方程的瞬变电磁场转换为满足波动方程的虚拟波场后，可以借助地震勘探中发展起来的一些成像方法，解决复杂情况下瞬变电磁场的成像问题[13-14]；K.H.Lee 等[10]建立了满足时域扩散方程的电磁场与虚拟波场的关系式，通过这个变换式，可以将对应地电模型的波动方程模拟结果变换成时域电磁响应，但该变换的意义并不在于对电磁场的正演模拟
煤田地质与勘探 2022年5期2022-06-07
- 双检数据上下行波场分离技术研究进展
器接收的地震压力波场数据，能够很好地展现出地下地质构造和反射界面信息。可是，来自海面的虚反射噪声干扰总是如影随行，成为单一水检数据地震勘探一个难以克服的难题。海洋地震勘探中，海洋震源子波是一个远场子波，包括震源虚反射。Loewenthal等[1]最早认识到速度波场可用于地震勘探，提出使用压力检波器和速度检波器进行海上和陆上地震勘探。Barr等[2]首先在海底电缆(Ocean Bottom Cable，OBC)勘探中使用压力和速度检波器记录水、陆检数据，压制
石油地球物理勘探 2021年6期2021-12-04
- 煤层气富集矿井采空区瞬变电磁波场变换成像技术
瞬变电磁信号中的波场特性，基于速度分析反演成像方法，通过对比压裂前后的数据对煤层气水力压裂的效果进行了评价。段建华等[10]采用瞬变电磁法和微震监测技术联合对井下煤层气水力压裂效果及影响范围进行监测。综合以上研究成果可以看出，目前利用瞬变电磁法进行煤层气富集矿井低阻目标体的探测多为基于传统数据处理与解释方法，受体积效应影响，视电阻率断面探测成果可反映低阻区域，但对于最为重要的采空区及采动裂隙边界位置则难以区分，且基于井下全空间瞬变电磁法的含水采空区及采动裂
重庆大学学报 2021年10期2021-11-09
- 基于超大内存节点的波动方程逆时偏移高效实现
程波动方程延拓和波场成像两步构成,先对震源波场采用双程波方程正向外推,并保存外推波场数据;然后对接收波场采用双程波方程反向外推,并在外推过程中读取记录的震源外推波场快照,应用成像条件获取成像值,所有时间步求和得到单炮成像数据体。随着油气勘探向深层目标体和精细勘探方向发展,野外勘探普遍采用高精度与高密度地震采集技术,地震数据规模不断增大。对逆时偏移计算而言,地震数据规模增大,使得单炮数据成像空间规模增大,波场外推计算量与波场数据量增大。对大规模地震数据进行逆
石油物探 2021年5期2021-09-28
- 瞬变电磁扩散场到虚拟波场的精细积分变换算法
扩散场转换到虚拟波场.只有在此基础上，才能实现稳定高效的瞬变电磁法拟地震解释.瞬变电磁扩散场主要描述的是低频电磁场传播过程中的扩散和感应特征，具有较强的体积效应，因而对电性界面的分辨能力较差(Khan et al.，2018；底青云等, 2019; Di et al.，2020；Xue et al.，2019；Li et al.，2018)；而波动场能够比较好地描述波场的一系列特征，所以基于波动方程的拟地震成像可以很好地刻画地层的界面信息(Qi et al
地球物理学报 2021年9期2021-09-06
- 地下黏弹性介质波动方程及波场数值模拟*
全弹性假设的模拟波场和实际传播特征之间的差异明显, 波动理论在工程技术中的应用效果还有提升的空间. 在岩石物理中,品质因子Q是量化地震衰减强度的参数, 为了研究波在地下介质中的传播规律, 本文从常Q理论出发, 在黏弹性介质频散关系中, 利用多项式拟合和Taylor展开法将频率的分数阶转化为整数阶, 进而推导了时间域复数形式的地下黏弹性介质波动方程. 该近似处理避免了频散关系经域转换后出现分数阶时间微分项, 能有效地降低计算成本. 最后, 采用有限差分法联合
物理学报 2021年14期2021-08-05
- 纵横波波场分离效果对比
以得到完整的地震波场响应[1]。除了常规网格外，Madariaga提出的交错网格，它在不增加计算量的前提下，和常规网格相比局部精度提高了4倍，收敛速度也很快，很好地压制了数值频散[2]。弹性波纵横波场分离是进行弹性波逆时偏移的必要条件，如若不能正确分离纵横波场，则在偏移成像过程中由于纵横波的互相干扰，导致偏移成像失败。目前大多数学者都采用传统的Helmholtz原理分离纵横波场。Helmholtz分解基于弹性波混合波场是由一个无旋场和无散场组成的，这个无旋
化工设计通讯 2021年4期2021-04-30
- 基于CLG光学流和波场分解的逆时偏移角度道集提取方法研究
波分解方法需要将波场分解为局部平面波分量，应用角度域成像条件提取角度道集，计算代价较大．在逆时偏移角度道集中，算法的高效性和并行性是制约高效提取角度道集的重要因素．波矢量方向估计方法(Yoon et al.，2004；Yoon and Marfurt，2006；Dickens and Winbow，2011；Vyas et al.，2011a，b；Zhang and McMechan，2011a，b；王保利等，2013；Jin et al.，2014；Yo
地球物理学报 2021年4期2021-04-07
- 基于Poynting矢量行波分离的最小二乘逆时偏移
双程波动方程进行波场延拓的成像方法，逆时偏移具有诸多优点：可以利用反射波、多次波、回转波和棱柱波等各种波型，不受地层倾角的限制，并且对起伏地表和复杂速度分布有很好的适应性.但是相较于单程波偏移算子，逆时偏移面临着由背向反射所引起的低频噪声发育的问题.所谓背向反射，是指在震源波场或者检波器波场延拓过程中，当波遇到波阻抗界面时所发生的次生反射.背向反射波与正常传播的波之间的互相关在地震剖面上所表现出的结果即低频噪声(刘红伟等，2010；杜启振等，2013).一
地球物理学报 2021年2期2021-02-23
- 水陆检数据上下行波场分离方法
检标定数据（上行波场数据），消除了海面多次波干扰；将标定后陆检数据与对应水检数据相减，得到与海面相关的多次波数据（下行波场），用于后续上行波场与下行波场联合反褶积、去噪和偏移等处理，以进一步提高分辨率和信噪比［11-28］。地震数据大多包含噪声，因此常规处理方法在使用能量匹配因子分离上、下行波场时，分离出的上行波场中包含下行波场、下行波场中包含上行波场，即不能完全实现上、下行波场的完全分离。针对常规方法不能完全分离上、下行波场的窘况，本文研讨了水陆检数据上
石油地球物理勘探 2020年5期2020-10-17
- 应用GPU 的傅里叶有限差分逆时偏移
差分方法进行正反波场传播，利用互相关成像条件获取成像值，进而通过滤波方式去除低频噪声［8-9］。这类方式在应用中存在两个局限：一是实际偏移时计算网格较大，显式有限差分法存在较严重的数值频散［10-12］，且产生无规则噪声，不利于高频信息的保护与成像；二是成像后滤波在去除低频噪声的同时，会损伤有效低频信息，不利于低频信息的保护与成像。空间差分数值频散由地震波的传播方向、空间差分精度和空间采样密度等三个因素决定。其中：因地震波向四周传播，故传播方向无法控制；提
石油地球物理勘探 2020年5期2020-10-17
- 虚拟波场变换方法在电磁法中的进展
，由此开展了虚拟波场变换研究。虚拟波场变换方法是指通过数学上的积分变换，将麦克斯韦扩散方程转换为与电磁波类似的虚拟波动方程，在虚拟波场中进行计算和分析，该方法的初始研究是借鉴波动方程的性质，通过虚拟波场变换将地震中一些成熟的方法引入到电磁法的数据解释中，提高电磁法数据解释的分辨率，后来研究发现它在提高数值模拟的计算效率以及其他方面也可发挥重要作用。电磁法虚拟波场变换的方法研究开始于上世纪70年代。Kunetz[12]首先发现了扩散的电磁场和波动方程之间存在
物探化探计算技术 2020年1期2020-04-08
- 基于行波分解的绕射波成像方法研究
,从而实现对绕射波场单独成像的目的[9-10]。因此,绕射波成像方法主要分为:数据域绕射波成像方法与成像域绕射波成像方法。数据域中的分离方法有:基于绕射波旅行时方程构建绕射波时间剖面的D-Section分离方法[11-12];聚焦-反聚焦(Focusing-Defocusing)方法[8],其利用反射波具有固定镜像虚震源而绕射波在伪深度剖面上发散的特点,将聚焦的反射波切除再反聚焦从而在全波场中去除反射波;平面波域中的绕射波分离方法[13],其基本思想是在平
石油物探 2020年2期2020-03-30
- 煤矿采空区瞬变电磁超前探测波场成像研究
瞬变电磁超前探测波场成像技术研究，对于提高采空区水害探测预警的准确性具有重要的实际意义。笔者针对煤矿井下巷道掘进超前探测方法特点，采用波场变换数值计算与成像技术进行数据处理，显著提高了方法的纵向分辨率，确定了巷道掘进前方采空区的富水边界，可指导煤矿巷道超前探测工作。1 全空间波场变换瞬变电磁场分量与其相对应的波场之间的数学积分形式为[7，8]：将式(1)的积分方程离散化，写成矩阵形式：AU=P；U=(u1，…，uj，…，un)T，P=(p1，…，pj，…，
煤炭工程 2020年2期2020-02-24
- 基于能量范数成像条件的一阶方程弹性波逆时偏移
包括震源、检波点波场重建和成像条件应用。不少学者在弹性波RTM过程中先进行纵横波波场分离然后再成像[8-10]。采用的主要方法是基于Helmholtz分解的波场分离[11],可分为直接波场分离法[12-14]和间接波场分离法。直接波场分离法是在空间域对波场进行散度和旋度运算分离得到纵、横波波场,但存在振幅相位校正问题。AKI等[15]选择空间域矢量分解方法解决了成像时极性反转问题,但其波场的物理意义发生了改变;ZHU[16]尝试在波数域进行波场分离,但该方
石油物探 2020年1期2020-02-07
- 一种基于均匀稀疏采样的Lamb波场重构方法
式传感器阵列拾取波场信号以定位损伤，但此类方法具有以下缺点：① 传感器需通过耦合剂与试件黏接，采样信号易受耦合剂影响；② 阵列排布形式固定，操作不灵活；③ 空间分辨率受限，难以针对损伤的尺度、形状等特征作精确评估[3]。使用非接触式扫描激光多普勒测振仪(Scanning Laser Doppler Vibrometer,SLDV)能够获得板状结构被测区域内高空间分辨率Lamb波场数据，并利用波场分析技术对损伤进行精确评估[4-6]。“波场”主要指行波在被测
无损检测 2019年9期2019-09-24
- 巷道瞬变电磁超前探测波场成像技术
果。近年来发展的波场变换与合成孔径等高精度解译成像技术，为提高矿井瞬变电磁法对含水异常体探测的分辨率提供了新的思路，由此可见，开展巷道瞬变电磁超前探测波场成像技术研究，对于提高矿井水害探测预警的准确性具有重要的实际意义。笔者针对煤矿井下巷道掘进超前探测方法特点，采用波场变换数值计算与成像技术进行数据处理，显著提高了方法的纵向分辨率，确定了巷道掘进前方含水异常体的富水边界，并以此指导某矿3201回风巷超前探测工作，为巷道掘进前方水害防治工作提供依据。1 全空
中国煤炭地质 2019年7期2019-09-16
- 孙宇晨 拿波场币和巴菲特吃顿饭
前不久，波场（TRON）创始人孙宇晨在微博上宣布，自己以历史新高456.7888万美元成功拍下了巴菲特20周年慈善午宴。不过，显然大家对此并不买账。消息公布数小时后，波场官方代币波场币（TRX）开始暴跌，仅1小时就狂跌13.28%。“馬云门徒”“币圈贾跃亭”“世界经济论坛全球杰出青年”……孙宇晨被冠上了太多标签，毁誉参半。不过，这并不影响他成为“90后”创业的标杆性人物。2018年1月5日，距离创立仅半年，波场的流通市值达到130亿美元，总市值达到250亿
当代工人 2019年14期2019-09-09
- 基于复数域波场分解的保幅逆时偏移成像方法
差分方法进行震源波场正传与检波点波场反传,并将两种波场沿时间方向做零延迟互相关成像[3-4]。利用双程波波场模拟算法能够正确描述反射波、折射波、回转波、多次波等复杂波场现象,处理复杂构造地区的地震成像问题。然而,由于波场传播算子没有明确区分波的传播方向,引入了偏移噪声[5],因此当入射波场与反射波场传播方向相同时,将产生物理上不存在的成像噪声,表现为低频强振幅掩盖了真实的反射界面[6]。针对RTM成像低频噪声问题的处理方法主要可以分为两类:①减少背向反射波
石油物探 2019年2期2019-04-10
- 基于LSQR算法的二维声波方程频率域正演模拟与数值实现
少成就。大多数的波场模拟都是在时间域进行的，时间域计算方法是按时间片递推计算，由前一个时间片的波场值推出下一个时间片的波场值。因此前一时刻的计算误差会累加到下一时刻，如果计算的时间切片数目较多，计算误差将会累计，导致计算结果信噪比下降、精度太低而无法为接下来的波形反演提供研究基础。近三十年来，学者做了大量工作，频率域正演模拟最早由Lysmer等[1]提出；Shin等[2]进一步发展了这种方法，将频率域有限元法应用于地震波形反演; Jo等[3]提出了最优化9
物探化探计算技术 2019年1期2019-03-01
- 波场TRON主网络虚拟机TVM正式发布 网络设施已现清晰架构
多的是智能水平。波场——作为行业内的佼佼者，波场虚拟机TVM于2018年7月30日上线全新版本。什么是虚拟机？虚拟机，英文全名Virtual Machine，指通过软件模拟的具有完整硬件系统功能的、运行在一个完全隔离环境中的完整计算机系统。首先是虚拟系统，通过生成现有操作系统的全新虚拟镜像，具有真实windows系统完全一样的功能，进入虚拟系统后，所有操作都是在这个全新的独立的虚拟系统里面进行，可以独立安装运行软件，保存数据，拥有自己的独立桌面，不会对真正
投资者报 2018年30期2018-08-03
- 波场分解算法与逆时偏移角道集
成像是通过对震源波场和检测器波场应用传统的零延迟互相关成像条件[5]得到的,如果存在强反射,这个过程不仅在实际反射界面处产生图像,也会沿着波路径产生严重的低频噪声。这些噪声可能掩盖真实的构造,给解释工作造成困难。将这个传统的成像条件和真振幅偏移方法[6-11]做比较,我们发现主要区别在于依赖于震源和检波器波场夹角的权重函数。实际上,真振幅偏移中关于反射角的函数可以压制很多反向散射噪声,关键在于如何计算逆时偏移中的反射角。②快速生成偏移共成像点道集问题。共成
石油物探 2018年4期2018-08-01
- 弹性波波场分离方法对比及其在逆时偏移成像中的应用
holtz分解的波场分离法。波场分离法主要有两类：一类是直接分离法，另一类是间接分离法。直接分离法是直接对矢量波场进行散度、旋度计算，得到纵、横波波场，主要包括空间域分离及波数域分离两种。最早提出的空间域分离法是直接对波场进行散度、旋度运算得到纵波标量波场和横波矢量波场，这类方法得到的波场存在振幅和相位畸变，需做校正[6-8]。在此基础上，Aki等[9]提出空间域矢量分解法，得到的纵波与横波均为矢量，且在同一分量上与原波场分量极性一致，可直接解决成像时极性
石油地球物理勘探 2018年4期2018-07-16
- 非分裂完全匹配层边界存储时间域全波形反演
题，提出利用阻尼波场零频分量反演低频模型作为频率域波形反演的初始模型，即Laplace域全波形反演。频率域全波形反演由于其固有的多尺度特征使其理论研究和实际应用得到快速发展[6-10]。然而频域率反演最大的问题就是大型稀疏方程组的存储以及求解过程中巨大的内存需求，大多借用MUMPS线性方程组求解软件包进行求解。但对于大规模尤其是三维情况下，频率域正演对内存的超大要求，一般的计算机集群仍无法满足要求，因此近年来时间域全波形反演研究成为热点[11，12]，尤其
石油地球物理勘探 2018年4期2018-07-16
- 高效逆时偏移角度道集生成方法研究
[25-34]和波场局部平面波分解方法[17,35-39]。扩展成像条件方法首先生成地下局部偏移距道集,然后将其转换为角度域成像道集,该转换过程需要进行高维傅里叶变换和局部倾斜叠加计算,存在内存和I/O问题。在进行局部倾斜叠加时,时窗大小对生成道集的质量影响很大,选择的局部时窗小时,角度道集噪声少,但分辨率低;选择的时窗大时分辨率高,但是倾斜叠加产生的假象多,从而降低了角度域成像道集的质量,影响AVA特征提取。逆散射成像条件方法是利用方位角度滤波器对空间求
石油物探 2018年3期2018-05-31
- 基于解耦传播的波场分解方法在VTI介质弹性波逆时偏移中的应用
是利用模拟的震源波场和观测到的检波点波场对地下构造和界面进行成像。与声波逆时偏移相比，弹性波逆时偏移更符合实际情况，能够处理多分量数据，而且转换波的成像结果有更高的分辨率。因此，近些年有很多关于弹性波逆时偏移的研究[1-6]。还有一些研究涉及到横向各向同性(VTI)介质[7-8]。然而，对于弹性波逆时偏移，在源波场和检波点波场中都存在P波和S波。为了减少成像结果中的串扰假象，在逆时偏移之前进行P波和S波分解就变得非常必要。传统的波场分解方法，比如Helmh
吉林大学学报（地球科学版） 2018年3期2018-05-26
- 基于波场重建的时域全波形反演
成像方法时，震源波场的存储量决定着RTM和FWI方法对计算硬件存储性能的要求。同理，针对RTM方法的存储策略同样可借鉴应用至FWI方法中。Clapp提出在速度外围构建随机速度层作为边界条件[7]，将边界波场能量以随机噪声的形式保存于波场内部，进而通过波动方程的反向计算即可获得历史波场，该方法由于波场存储量较低，被广泛应用于RTM方法中；Shen等对上述方法进行了改进，加强了随机边界对低频数据的散射强度，并应用于波形反演算法中[8]。但随机边界条件产生的随机
中国锰业 2018年2期2018-05-23
- 随机边界的波场重建技术研究
偏移法采用的是全波场波动方程，对波动方程中的微分项采用差分离散(如有限元法、伪谱法、有限差分法等),其对波动方程的近似比较少(这主要由微分与差分方程的逼近程度来决定)，所以其不受偏移孔径和倾角的约束，在处理纵横向上有剧烈变化的地球介质物性特征(如密度、速度等)上很有成效。但是，计算量大和所需存储量高使其无法适用于大数据量的地震成像。利用随机边界波场重建技术，可以最大程度降低计算量和存储量，使得在实际操作过程中可以有效地利用地震波逆时偏移法。因此，对随机边界
中国锰业 2018年2期2018-05-23
- 基于检查点方法的各向异性介质逆时偏移
性,势必影响地震波场模拟精度和逆时偏移成像分辨率。虽然各向异性介质弹性波偏移理论已经取得了一定突破[5-7],但在应用过程中依然会遇到障碍。其一,弹性波是由一个三分量的位移矢量表征,包括两种模式即纵波和横波,如利用有限差分法求解波场,高昂的计算成本和巨大的存储量令人难以接受;其二,考虑横波速度较慢,为了满足波场稳定延拓条件,须采用较小时间步长和空间网格间距,必然会影响计算效率;其三,纵波与横波分离以及实际弹性参数的确定等都是其制约因素[8]。因此,对于地震
石油物探 2018年2期2018-04-04
- 基于波长平滑算子的逆时偏移回转波假象去除方法
自炮检点同向传播波场产生的低频噪声受到广泛关注，并涌现出多种消除此类低频噪声的方法。Loewenthal等[21]采用平滑速度方式消除背向散射波场进而压制低频噪声； Mulder等[22]通过带通滤波器对偏移后的成像结果进行滤波处理以消除低频假象； Valenciano等[23]利用反褶积成像条件压制低频噪声； Fletcher等[24]提出采用方向性衰减方法压制层间反射从而去除RTM中的低频假象； Yoon等[25]将Poynting矢量引入成像条件以消
石油地球物理勘探 2018年1期2018-03-10
- 基于稀疏表示的非零井源距VSP波场分离方法
非零井源距VSP波场分离方法赵彦青*萧蕴诗（同济大学电子与信息工程学院，上海201804）随着垂直地震剖面（VSP）技术的不断发展，高保真的波场分离成为VSP资料处理与应用的关键之一。引入欠定盲源分离思想，利用VSP波场的极化特征及多分量资料在时频域的稀疏性，提出一种基于稀疏表示的波场分离方法。该方法分为两步：第一步将地震信号变换到时频域，并估计极化矩阵；第二步在时频域分离各单一波场，再反变换回时间域。在第二步波场分离过程中，针对时频域波场稀疏性不足的情况
石油地球物理勘探 2017年1期2017-10-23
- 地震叠前逆时偏移衰减随机边界条件研究
边界条件,在真实波场外围依次加入衰减层和随机速度层,可有效降低随机反射对真实波场的能量干扰。对波动方程逆向运算重建激发点波场,结合互相关成像条件,完成逆时偏移计算。对比分析随机边界条件、完全匹配层边界条件和衰减随机边界条件的测试结果,与随机边界条件相比,衰减随机边界条件在牺牲较少存储效率的情况下,能够获得接近吸收边界条件的偏移效果。盐丘模型和大偏移距模型算例分析结果验证了以上结论的正确性。逆时偏移;随机边界;存储;波场重建逆时偏移[1]对复杂构造成像具有良
石油物探 2017年4期2017-08-30
- 基于随机采样的频率域多路径波场模拟与偏移成像
样的频率域多路径波场模拟与偏移成像李 博1，2（1．中国石油化工股份有限公司石油物探技术研究院，江苏南京211103；2．中国科学院地质与地球物理研究所油气资源研究重点实验室，北京100029）采用射线追踪方法的常规Kirchhoff深度偏移不能全面描述焦散区的波现象，有可能导致偏移成像产生畸变，利用此类方法进行复杂介质条件下的偏移成像，不可避免地存在焦散问题。详细推导了二次震源的波场模拟理论方法与区域分解的实施方案，在确保局部空间不存在射线交叉的前提下进
石油物探 2017年3期2017-06-29
- 双相介质波场数值模拟分析
制边界反射，通过波场数值模拟分析反射振幅的变化规律。关键词：双相介质；波场；模拟；振幅含流体双相介质是非常复杂的介质，双相介质中地震波传播理论是目前地震学研究的前沿和热门的问题，也是难点之一，其理论和应用正在进一步的完善和发展。理论描述了饱和流体孔隙介质中弹性波的传播，是双相介质地震波传播的理論基础。同时考虑了流动和喷射流动力学机制的模型是对理论的完善和发展。储层地球物理学表明，油气储层的地质模型主要可分为砂岩储层模型和碳酸盐岩储层模型，这两种模型均为含流
报刊荟萃（上） 2017年3期2017-06-26
- 瞬变电磁高分辨率三维成像技术研究
77）介绍了瞬变波场变换与高分辨率波场反变换，给出了三维成像的原理和实现过程。利用瞬变电磁三维成像技术进行模型数据以及数据的处理。利用三维成像技术在模型的处理中得到了很好的效果，并依此来探测矿区采空区分布，结果证明使用该方法进行数据处理达到了很好的效果。瞬变电磁法；波场变换；煤矿；偏移瞬变电磁法在水资源调查、矿产勘查、灾害地质调查、考古探测等诸多领域。目前，瞬变电磁数据处理基本停留在数据一维正演、一维反演、二维正演、二维反演阶段，对于数据的三维正演有一部分
西部探矿工程 2016年12期2016-12-21
- 多道瞬变电磁法(MTEM)虚拟波场提取技术
(MTEM)虚拟波场提取技术钟华森1,2， 薛国强1*， 李貅3， 智庆全4， 底青云11 中国科学院地质与地球物理研究所矿产资源研究重点实验室, 北京 100029 2 中国科学院大学， 北京 100049 3 长安大学地质工程与测绘学院, 西安 710054 4 中国地质科学院地球物理地球化学勘查研究所, 河北廊坊 065000多道瞬变电磁法(MTEM)是近年来发展起来的一种新的地球物理勘查技术，其数据采集方式与地震法类似，因此，采用瞬变电磁拟地震解释
地球物理学报 2016年12期2016-12-07
- 交错网格与旋转交错网格对VTI介质波场分离的影响分析
网格对VTI介质波场分离的影响分析1) 中国山东青岛266580中国石油大学(华东)地球科学与技术学院2) 中国天津300451中海油田服务股份有限公司利用交错网格有限差分和旋转交错网格有限差分进行各向异性介质弹性波场数值模拟时， 质点振动速度分量与应力张量的网格节点定义方式均不相同， 从而对各向异性波场分离效果产生不同的影响． 针对这一问题， 本文以具有垂直对称轴的横向各向同性(VTI)介质的波场分离为例， 首先分析了两种网格的参数定义方式以及VTI介质
地震学报 2016年1期2016-11-28
- 基于Hilbert变换的全波场分离逆时偏移成像
bert变换的全波场分离逆时偏移成像王一博1，郑忆康1，薛清峰1，常旭1，Tong W.Fei2，Yi Luo21 中国科学院地质与地球物理研究所，北京 1000292 EXPEC ARC,Saudi Aramco，Dhahran 31311，Saudi Arabia逆时偏移方法利用双程波算子模拟波场的正向和反向传播，通常采用互相关成像条件获得偏移剖面，是一种高精度的成像方法.但是传统的互相关成像条件会在偏移结果中产生低频噪声；此外，如果偏移速度中存在剧烈
地球物理学报 2016年11期2016-11-23
- 时间二阶积分波场的全波形反演
7时间二阶积分波场的全波形反演陈生昌， 陈国新浙江大学地球科学学院， 杭州310027通过对波场的时间二阶积分运算以增强地震数据中的低频成分，提出了一种可有效减小对初始速度模型依赖性的地震数据全波形反演方法—时间二阶积分波场的全波形反演方法.根据散射理论中的散射波场传播方程，推导出时间二阶积分散射波场的传播方程，再利用一阶Born近似对时间二阶积分散射波场传播方程进行线性化.在时间二阶积分散射波场传播方程的基础上，利用散射波场反演地下散射源分布，再利用波
地球物理学报 2016年10期2016-11-08
- 有限差分法在波场模拟中的分析研究
通过该方法能求得波场正演的数值解。关键词：有限差分；离散；近似；波场一、引言地震波是指从震源产生向四外辐射的弹性波。地震发生时，震源区的介质发生急速的破裂和运动，这种扰动构成一个波源。由于地球介质的连续性，这种波动就向地球内部及表层各处传播开去，形成了连续介质中的弹性波。波场、地震波场，指有地震波传播的空间。在这个空间的每一点上，一定时刻都有一定的波前通过，波的能量也按—定的规律传播。当我们按照顺时间轴方向研究波场的状态问题的时候，如果研究的就是波场的正演
大东方 2016年1期2016-05-30
- 基于多卡GPU集群的多次波逆时偏移成像技术
中的多次波预测、波场延拓等步骤的实现过程进行了研究。在此基础上，利用Sigsbee2B地质模型，应用基于反馈环理论的自由界面多次波衰减方法预测得到纯多次波炮集，实现了基于多卡GPU集群的多次波逆时偏移成像处理。处理结果显示，多次波逆时偏移的中浅层成像清晰，构造的成像精度明显高于常规逆时偏移。多卡GPU集群的应用可显著提高多次波逆时偏移的计算速度和效率，使得基于地震大数据体进行多次波逆时偏移成像处理成为可能。多次波波场延拓逆时偏移成像多卡GPU成像精度逆时偏
油气地质与采收率 2015年2期2015-10-18
- 基于解析时间波场外推与波场分解的逆时偏移方法研究
2基于解析时间波场外推与波场分解的逆时偏移方法研究胡江涛, 王华忠同济大学海洋与地球科学学院波现象与反演成像研究组， 上海 200092双程波方程逆时深度偏移是复杂介质高精度成像的有效技术，但其结果中通常包含成像方法引起的噪音和假象，一般的滤波方法会破坏成像剖面上的振幅，其中的假象也会给后续地质解释带来困扰.将波场进行方向分解然后实现入射波与反射波的相关成像能够有效地消除这类成像噪音，并提高逆时偏移成像质量.波传播方向的分解通常在频率波数域实现，它会占用
地球物理学报 2015年8期2015-03-01
- 地震波动方程方向行波波场分离正演数值模拟与逆时成像
的同时减少了虚假波场的引入，有效提高了计算结果的信噪比。地震波正演数值模拟技术的引入不仅大大降低了物理模拟技术和实际野外地震勘探的成本，同时还可有效指导实际地震资料的采集、处理、解释及方法验证等，因此，该技术在当前勘探地震学领域中发挥着至关重要的作用。随着地震波数值模拟技术的进步，如何在地震波正演模拟过程中实现不同类型波场的自动识别与分离问题逐渐得到广泛重视，有利于对复杂波场传播规律的认识和应用。其中，最具代表性的是弹性波波场分离数值模拟方法，目前该方法包
岩性油气藏 2014年4期2014-12-06
- 高分辨率拉东变换在阵列声波测井波场分离中的应用
记录的地层信息以波场的形式叠加在一起。波场分离效果直接影响到信号提取准确程度，关系到后续处理与解释的可靠性。波场分离技术主要包括中值滤波、f－k滤波、拉东变换等处理方法。阵列声波测井波场分离主要采用速度偏移法［1］，f－k滤波、拉东变换等技术应用较少。速度偏移法需要已知声波时差与窗长，条件相对苛刻；中值滤波［2］用于分离直达波与反射波，不易分离上、下行反射波；f－k滤波对滤波器依赖较大，滤波器带通过窄，易出现Rieber现象［3－4］，而带通过宽则压制噪声
测井技术 2014年4期2014-12-03
- 地震勘探基于波场分离的逆时偏移去噪方法
术一般分为有基于波场延拓的波动方程偏移方法和基于射线理论的Kirchhoff偏移方法两类。波动方程偏移方法是对波动方程的进行分解，能较好的适应速度模型的横向变化，但由于采用单程波的近似方程，在波场传播方向受到倾角限制，对陡倾角构造成像能力较差，尤其对于盐丘下回折波成像存在很大的困难。Kirchhoff偏移方法有很高的计算效率和实用性，在当前实际生产中广泛应用，但Kirchhoff偏移方法的分辨率会随着深度增加变差且缺乏相应的振幅信息。随着油气勘探难度的增大
云南地质 2014年2期2014-08-03
- 一种改进的中值滤波方法
包含了十分丰富的波场信息，它不仅能接收到上行纵波和上行转换波，也能接收到下行纵波、下行转换波及横波。在VSP资料处理中，有效地分离上行波和下行波，以及分离纵波和横波，是一项重要的任务。波场分离是VSP资料处理中至关重要的环节，常用的波场分离方法有F－K滤波和中值滤波。F－K滤波方法［2］是线性滤波，在频率－波数域中干扰波与有效波的视速度不同，若下行波的波数为正，则上行波的波数为负，下行波和上行波分别位于2个不同的区域，所以可以较好地衰减下行波从而突出上行波
成都理工大学学报（自然科学版） 2013年2期2013-08-01
- 逆时偏移计算中的边界处理分析及应用
低频成像噪音，和波场延拓方向不一致带来的存储等问题需要解决［6］.最近几年计算机技术的飞速发展，让逆时偏移重新获得生命力，它所存在的问题也得到了长足的改善.在处理计算量的问题上，通过计算机的并行计算［7］，应用GPU计算技术可以大幅提高计算效率［8－9］.在处理低频成像噪音的问题上，有在波传播过程中压制产生噪音的方法［10－12］，有对成像条件做出改进的方法［13－14］，有成像后滤波的方法［15－17］等.在解决由波场延拓方向不一致带来的波场传播历史的存
地球物理学报 2013年6期2013-04-11
- 基于GPU并行加速的叠前逆时偏移方法
双程声波方程进行波场延拓，突破倾角限制，借助于高阶有限差分方法实现叠前逆时偏移成像；利用GPU（Graphic Processing Unit）并行加速技术对波场延拓和成像进行计算，相比于传统算法，其计算效率有较大提高，可以解决叠前逆时偏移算法计算量过大问题；在获取波场信息过程中，也采用随机边界条件，实施以计算换存储策略，解决逆时偏移计算中的海量存储问题.模型测试结果表明，该方法能够高效和高精度地对地下复杂地质体成像.逆时偏移；GPU；加速；高阶有限差分；
东北石油大学学报 2012年4期2012-12-14
- 子波延迟对逆时偏移成像精度的影响
，这需要进行上行波场的反向外推。外推后求出的各点波场值，有的是来自本点的反射波，有的是该点下方许多点上的反射波叠加。因此，需要在外推波场中提取成像值。反射成像原则就是：反射面由反射点组成，在这些点上，其入射波的初至与反射波的产生时间相同，反射因数为反射波振幅u（x，z，td）（其中，x、z分别为水平坐标、深度坐标，m；td为反射时间，s）与入射波振幅d（x，z，td）的比值。因此，反射波成像可以用上下行波场进行互相关提取成像值。下行波初始时间的确定决定了成
石油天然气学报 2012年11期2012-08-20
- τ- p变换与偏振分析相结合波场分离方法的研究
与偏振分析相结合波场分离方法的研究张小慧，刘明洋(地球探测与信息技术教育部重点实验室(成都理工大学)，四川 成都 610059)王立娟(四川省安全科学技术研究院，四川 成都 610045)波场分离是多波多分量地震资料处理的重要环节之一，VSP三分量资料的各种类型的波相互叠加在一起，因此有效地分离出上下行波以及纵横波是VSP资料处理的关键。τ-p变换与偏振分析相结合的波场分离方法能够使地震数据按不同相轴呈点状收敛于不同的区域，用该方法进行VSP波场分离可使能
长江大学学报(自科版) 2011年10期2011-11-20
- 有限差分波动方程偏移波场特征定量分析
畸变。因为在实际波场成像过程中，波的动力学参数是随倾角变化的，改变各偏移处理参数也无法避免这种变化，从而导致偏移剖面各倾角成像波场的特征不同，这对构造解释影响不大，但在提取地震波动力学参数进行储层研究时，这种影响就要引起足够的重视，否则会在解释时无法正确判别参数的可靠性。针对上述情况，笔者在前人研究的基础上，通过点脉冲模型有限差分波动方程偏移处理，定量分析偏移后各倾角波场主振幅、频率和相位等参数的相对变化特征，讨论了这种变化对波场的影响程度。1 有限差分波
长江大学学报(自科版) 2011年28期2011-04-14
- 纵、横波分离的数值模拟研究
同性介质中的地震波场，并分别计算散度场和旋度场。结果表明，在满足差分范围内速度无剧烈变化的情况下，通过该方法能有效地分离出纵、横波场。另外，旋度运算的矢量性质造就了分离后的记录在震源两侧极性反向的特性，而这一特性与野外采集的横波数据在物理特性方面极其相似。弹性波方程；纵波；横波；散度；旋度通过弹性波方程数值模拟得到的是P波(纵波)和S波(横波)的混合波场。对于由P波震源激发记录中包含反射P波和转换S波的地震数据，可以通过散度和旋度的计算将P波和S波分离❶。
长江大学学报(自科版) 2010年1期2010-11-27