发明内容
本发明实施例提供了一种地震资料处理真地表条件下的振幅补偿方法和装置。所述技术方案如下:
一种地震资料处理真地表条件下的振幅补偿方法,包括:
获取真地表条件下的地震资料;
根据获取的所述地震资料计算得到检波点的出射角;
根据所述出射角和所述检波点处的接收能量对所述接收能量进行补偿;
根据所述补偿后的接收能量和成像时间信息计算得到零偏移距时间剖面。
所述获取的地震资料包含经过处理的地表函数,相应地,根据获取的所述地震资料计算得到检波点的出射角,具体包括:
根据所述经过处理的地表函数计算检波点的法线斜率;
根据所述地震资料计算所述检波点与其相应的反射点之间连线的斜率;
根据所述法线斜率和所述连线的斜率计算得到检波点的出射角。
根据所述出射角和所述检波点处的接收能量对所述接收能量进行补偿,具体包括:
将所述检波点处的接收能量除以所述出射角的余弦值,得到振幅补偿后的接收能量。
一种地震资料处理真地表条件下的振幅补偿装置,包括:
地震资料获取模块,用于获取真地表条件下的地震资料;
第一计算模块,用于根据获取的所述地震资料计算得到检波点的出射角;
能量补偿模块,用于根据所述出射角和所述检波点处的接收能量对所述接收能量进行补偿;
第二计算模块,用于根据所述补偿后的接收能量和成像时间信息计算得到零偏移距时间剖面。
所述获取的地震资料包含经过处理的地表函数,相应地,所述第一计算模块包括:
第一斜率计算单元,用于根据所述经过处理的地表函数计算检波点的法线斜率;
第二斜率计算单元,用于根据所述地震资料计算所述检波点与其相应的反射点之间连线的斜率;
出射角计算单元,用于根据所述法线斜率和所述连线的斜率计算得到检波点的出射角。
所述能量补偿模块具体用于将所述检波点处的接收能量除以所述出射角的余弦值,得到振幅补偿后的接收能量。
本发明实施例提供的技术方案的有益效果是:
通过地震波反射射线出射角度和出射能量计算出补偿后的接收能量,能够有效补偿真地表、特别是起伏地表、高速岩石直接出露等不满足地表一致性和垂直出射假设条件地区采集的地震资料因出射角导致的能量损失,提高振幅处理保真度和成像质量。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
实施例1
为了有效补偿地震资料振幅保持中不满足地表一致性和垂直出射假设条件地区因出射角导致的能量损失,提高地震资料处理振幅保真度和成像质量,本发明实施例提供了一种地震资料处理真地表条件下的振幅补偿方法,参见图1,该方法包括:
101:获取真地表条件下的地震资料;
102:根据获取的所述地震资料计算得到检波点的出射角;
103:根据所述出射角和所述检波点处的接收能量对所述接收能量进行补偿;
104:根据所述补偿后的接收能量和成像时间信息计算得到零偏移距时间剖面。
通过地震波反射射线出射角度和出射能量计算出补偿后的接收能量,能够有效补偿真地表、特别是起伏地表、高速岩石直接出露等不满足地表一致性和垂直出射假设条件地区采集的地震资料因出射角导致的能量损失,提高振幅处理保真度和成像质量。
实施例2
为了有效补偿地震资料振幅保持中不满足地表一致性和垂直出射假设条件的地区因出射角导致的能量损失,提高地震资料处理振幅保真度和成像质量,本发明实施例提供了一种地震资料处理真地表条件下的振幅补偿方法,该方法以上述真地表条件下的椭圆展开成像理论和真地表条件下的参数展开成像理论为基础,由于在两个成像理论中的补偿实现方法是相同的,此处以前者为例,参见图2,该方法包括:
步骤201:采集地震原始资料;
在地震勘探领域中,该地震原始数据的采集主要是利用地震勘探设备,在野外地面设计的炮线上人为激发地震波,在预勘探区域设置多个接收点,该多个接收点散布在包括炮点在内的地面上一个勘探目标区域范围的矩形网格点上,使用地震检波设备获取爆炸波即地震波的波形和记录的炮点、检波点的大地坐标及炮点、检波点的相对关系等地震原始资料。
步骤202:对采集到的地震资料进行室内解编处理;
本领域技术人员可以获知,室内解编处理是指将在野外接收到的地震数据拿到室内后,将这些以segd或segy格式存储的数据转换成便于地震资料处理软件处理的内部格式数据。
步骤203:加载观测系统,预处理解编后的资料;
其中,加载观测系统的目的是利用观测系统信息,建立接收到的地震数据与大地坐标的关系,得到目前所处理的资料所属的具体位置等信息。
其中,预处理是指对地震资料进行前期处理,包括去噪,去多次波,能量补偿,反褶积等,目的是去伪存真,使希望得到的信息更清晰可靠,为后续处理准备好数据。
其中,加载观测系统是对原始资料进行观测系统定义的过程,在这个过程中,具体地,将地震采集时记录的炮点、检波点的大地坐标及炮点、检波点的相对关系置于采集到的单炮的每一个记录道的道头中。
步骤204:根据预处理的结果求取地表函数;
具体地,利用观测系统信息,特别是记录到的炮点和检波点的坐标及高程信息通过插值和平滑得到地表函数z=f(x)。
步骤205:计算地表检波点位置的法线的斜率;
具体地,求出步骤204中得到的地表函数的在每个检波点处的一阶导数,该一阶导数即为每个地表检波点处的法线的斜率kn,其具体计算方法可以为kn=-dx/dz。
图3a所示为真地表椭圆展开成像算法基本原理图。均匀介质,波的传播速度为V,法线在S-R线上的出露点到炮点距离为10,对t时刻的信号进行椭圆展开变换可以根据已知公式求出法线的长度V*t0/2。如果V*t、V*t0/2、10和1为已知,那么就能知道可能的反射点O的位置。过S、R和O三个点做一个圆,圆的半径为r。做S-R的中垂线,它与圆的交点分别为B、E。可以证明,O点的法线(即SOR的角平分线)与圆的交点也在B点。因为SB弧和BR弧相等,反射角与入射角相等。
本领域技术人员可以获知,斜率的计算可以通过导数计算或者对坐标进行运算得到,本发明实施例对此不做具体限定。
步骤206:根据反射点O和相应的检波点R位置,求取OR线斜率k1;
图3b是真地表地震资料采集时地震波射线传播示意图。其中,S为炮点,R为检波点,O为反射点,地震波上行射线与其在地表出射位置的法线的夹角组成出射角i。出射角越大,地震波出射能量损失越大,需要补偿的能量也就越多。若地表起伏变化剧烈,则出射角也变化莫测,必然导致接收能量极不均匀,常规振幅补偿方法很难做到较为准确、合理的保真振幅补偿处理,严重影响成像质量。要做到接收能量保真处理,需要针对每一反射点以及对应的炮检对计算出出射角。这在常规以CMP(Common midpoint,共中心点)理论和水平地表假设条件为基础的处理方法很难实现。而依托地震资料处理真地表条件下的椭圆展开成像方法,可以很容易的计算出每一炮检对和相应地下反射点所对应的出射角,进行相应的接收能量补偿。
步骤207:根据求取的斜率,计算相应地出射角余弦;
优选地,出射角i的余弦可以由下式得到,
步骤208:根据该检波点处的接收能量和出射角余弦,计算振幅补偿后的接收能量。
具体地,将该检波点处的接收能量Eu除以出射角余弦cosi,得到振幅补偿后的接收能量Eu=Er/cosi,再通过干涉叠加,最终得到零偏移距时间剖面的振幅特征。
对于非三分量地震检波器,它们在地表接收到的能量Er是上行射线的能量Eu在法向的投影分量Er=Eu×cosi,则补偿后的接收能量应为Eu=Er/cosi。
本领域技术人员可以获知,本发明实施例中所述的振幅补偿是对接收能量的补偿。
步骤209:根据补偿后的接收能量和成像时间信息计算得到零偏移距时间剖面;
其中,成像时间信息由地震资料根据地震资料处理真地表条件下的椭圆展开成像方法计算得到,地震资料处理真地表条件下的椭圆展开成像方法包括:根据炮点的位置坐标、检波点的位置坐标和反射点的位置坐标,计算得到该反射点的法线在炮检线上的出露点的位置坐标;根据预设的基准面和该反射点的法线方程,计算得到成像位置点的位置坐标;根据该反射点的位置坐标,该成像位置点的位置坐标,以及该反射点与该出露点之间的双程旅行时间,计算得到该出露点与该成像位置点之间的时间校正量;根据该出露点与该成像位置点之间的时间校正量,考虑接收能量补偿,进行椭圆展开相切干涉叠加,得到零偏移距时间剖面。
该地震资料处理真地表条件下的椭圆展开成像方法,通过获取到出露点与成像位置点之间的时间校正量,根据出露点与成像位置点之间的时间校正量,考虑接收能量补偿,进行椭圆展开相切干涉叠加,得到零偏移距时间剖面,可以处理真地表条件下采集的地震资料,处理结果客观真实,当地表和地下都比较复杂时,也能得到较好的保幅零偏移距时间剖面,可以满足复杂地表和复杂构造成像的问题。并且可有效的避免常规处理方法存在的导致构造畸变、不保幅等问题,对起伏地表地区的油气和矿产资源勘探等有着重要的实际应用价值。另外该方法事先无须对地震资料进行任何的静校正处理,直接从起伏地表开始展开处理,将“静”校正量隐含地包括在其中,这种时间域校正不但包含了旅行时的纵向分量,同时还包含了旅行时的横向分量。
通过地震波反射射线出射角度和出射能量计算出补偿后的接收能量,能够有效补偿真地表、特别是起伏地表、高速岩石直接出露等不满足地表一致性和垂直出射假设条件地区采集的地震资料因出射角导致的能量损失,提高处理保真度和成像质量。
实施例3
为了有效补偿地震资料振幅保持中不满足地表一致性和垂直出射假设条件的地区因出射角导致的能量损失,提高地震资料处理振幅保真度和成像质量,本发明实施例提供了一种地震资料处理真地表条件下的振幅补偿装置,参见图4,该装置包括:
地震资料获取模块401,用于获取真地表条件下的地震资料;
第一计算模块402,用于根据获取的该地震资料计算得到检波点的出射角;
能量补偿模块403,用于根据该出射角和该检波点处的接收能量对该接收能量进行补偿;
第二计算模块404,用于根据该补偿后的接收能量和成像时间信息计算得到零偏移距时间剖面。
该获取的地震资料包含经过处理的地表函数,相应地,参见图5,该第一计算模块402包括:
第一斜率计算单元402a,用于根据该经过处理的地表函数计算检波点的法线斜率;
第二斜率计算单元402b,用于根据该地震资料计算该检波点与其相应的反射点之间连线的斜率;
出射角计算单元402c,用于根据该法线斜率和该连线的斜率计算得到检波点的出射角。
该能量补偿模块403具体用于将该检波点处的接收能量除以该出射角的余弦值,得到振幅补偿后的接收能量。
本实施例提供的装置,与方法实施例属于同一构思,其具体实现过程详见方法实施例,这里不再赘述。
为了测试本发明实施例提供的振幅补偿方法的实用性,本发明实施例进行了实验,实验结果如下:
参见图6,图6为一个地表起伏很大的简单地质构造模型示意图,其中,地表高程最低640m,最高1800m,高差达到近1200m;模型宽25km;地下有两个反射层,第一反射层纵波速度4000m/s,第二反射层纵波速度6000m/s。图7a所示为某炮地震波射线追踪图。对该模型进行高斯射线束正演得到的单炮地震记录如图7b。反射同相轴不再是类双曲线,而是发生了严重的畸变,这是由于炮点和检波点位于剧烈起伏的地表上。从图中不难看出,在同相轴上斜率较小的位置能量强,而斜率较大的位置能量弱。充分体现了出射角度对接收能量的影响。图中的红色同相轴表示补偿的能量,可以看到斜率大的位置能量补偿多,斜率小的位置能量补偿少。黑色和红色叠加后,接收能量基本均匀,说明我们所采用的接收能量补偿方案是真实有效的。
图8a和图8b为真地表椭圆展开成像(图8a)和常规CMP全部静校正处理后的最终叠加效果对比图(图8b),基准面高程定在800m。常规静校正方法所用的近地表替代速度为4000m/s。但是在成像剖面上几乎看不到有效信号。其可能原因是:由于该测线段地表地形极为复杂,低降速带较薄,多有高速刚性岩层直接出露地表,近地表速度横向变化较大,若假定射线在地表垂直出射,单纯采用高程静校正的办法做处理会带来很大误差。尽管后续做了折射静校正、剩余静校正等处理,但由于该地区完全不符合CMP理论的处理和成像条件,而且检波器的特征是(方向)选择性接收,接收的能量与出射角有关等,导致不能成像。
真地表椭圆展开成像方法直接对未做静校正的数据进行处理,可保证在射线追踪的过程中采用正确的炮点和检波点位置,而且成像过程中加入了与真地表相关的接收能量补偿,消除了地表起伏对接收能量的影响,大大提高了成像剖面的保幅性和质量。成像剖面上可以看到同相轴,说明采集到的地震资料中确实有有效信号。
本发明实施例提供的上述技术方案的全部或部分可以通过程序指令相关的硬件来完成,所述程序可以存储在可读取的存储介质中,该存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。