CN110879414A - 一种成像振幅一致性校正方法及校正系统 - Google Patents
一种成像振幅一致性校正方法及校正系统 Download PDFInfo
- Publication number
- CN110879414A CN110879414A CN201811033328.8A CN201811033328A CN110879414A CN 110879414 A CN110879414 A CN 110879414A CN 201811033328 A CN201811033328 A CN 201811033328A CN 110879414 A CN110879414 A CN 110879414A
- Authority
- CN
- China
- Prior art keywords
- transformed
- transformation
- coordinate
- point
- imaging
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 126
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000012937 correction Methods 0.000 title claims description 31
- 230000005012 migration Effects 0.000 claims abstract description 52
- 238000013508 migration Methods 0.000 claims abstract description 52
- 230000001788 irregular Effects 0.000 claims abstract description 50
- 230000001131 transforming effect Effects 0.000 claims abstract description 9
- 230000009466 transformation Effects 0.000 claims description 166
- 230000006870 function Effects 0.000 claims description 61
- 238000013507 mapping Methods 0.000 claims description 52
- 239000000126 substance Substances 0.000 claims description 18
- 238000005070 sampling Methods 0.000 claims description 17
- 238000011426 transformation method Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 abstract description 10
- 230000008569 process Effects 0.000 abstract description 9
- 238000004458 analytical method Methods 0.000 abstract description 5
- 230000009286 beneficial effect Effects 0.000 abstract description 4
- 239000000523 sample Substances 0.000 description 45
- 238000001514 detection method Methods 0.000 description 15
- 230000008859 change Effects 0.000 description 12
- 238000009792 diffusion process Methods 0.000 description 6
- 230000010354 integration Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种成像振幅一致性校正方法及校正系统,包括:将不规则的观测系统变换为规则的观测系统,获得雅克比行列式;计算雅克比行列式的值;针对当前地震道,将雅克比行列式的值作为加权因子调整Kirchhoff积分偏移公式的积分系数,获得变换后的Kirchhoff积分偏移公式;根据变换后的Kirchhoff积分偏移公式对当前地震道成像振幅进行调整。本发明通过雅克比行列式作为加权函数调整偏移过程中每一道的成像权重值,从而调整每一地震道的成像振幅对最终叠加结果的贡献权重,最终消除了观测系统不均匀导致的成像振幅不一致现象,实现了保幅成像目的,有利于后续的AVO分析及其它地震解释工作,且该方法操作方便,效率更高,不会引入过多的计算量。
Description
技术领域
本发明属于地震勘探领域,具体涉及一种成像振幅一致性校正方法及校正系统。
背景技术
Kirchhoff积分法偏移(包括时间偏移和深度偏移)是工业界中广泛应用的地震成像处理技术。该技术采用逐道成像、积分累加的方式进行。而在Kirchhoff积分法成像过程中假设观测系统是规则的,即炮点间隔以及检波点间隔(或CMP间隔)为一常数,最终叠加结果一般是对各个地震道的成像结果的简单叠加,不做加权区分。但实际采集过程中很难保证等间隔采样,实际观测系统一般是不规则的。因此,直接采用原始的Kirchhoff积分方程进行偏移得到的成像振幅一致性遭到破坏,而这种振幅不一致性是由于观测系统不均匀导致的,不能反应地下的真实反射情况,这些不准确的振幅能量关系会误导后继的AVO分析及各种地震解释工作。
现有做法一般是对叠前数据进行规则化处理再进行常规的成像,但是对叠前数据进行规则化处理的计算量较大且会生成更大空间的数据量,不利于目标快速处理。因此,特别需要一种效率更高的校正方法来实现成像振幅一致性。
发明内容
本发明的目的是提出一种效率更高的成像振幅一致性校正方法及校正系统。
根据本发明的一方面,提出了一种成像振幅一致性校正方法,包括:将不规则的观测系统变换为规则的观测系统,获得雅克比行列式;计算所述雅克比行列式的值;针对当前地震道,将所述雅克比行列式的值作为加权因子调整Kirchhoff积分偏移公式的积分系数,获得变换后的Kirchhoff积分偏移公式;根据所述变换后的Kirchhoff积分偏移公式对所述当前地震道成像振幅进行调整。
优选的,采用线性坐标变换方法将不规则的观测系统变换为规则的观测系统,获得雅克比行列式,所述雅克比行列式为:
其中,|J|为雅克比行列式,为变换后的CMP点,为变换后的半偏移距,为变换后的CMP点的x坐标,为变换后的CMP点的y坐标,为变换后的半偏移距的x分量,为变换后的半偏移距的y分量,f为CMP点x坐标变换前与变换后的映射函数,g为CMP点y坐标变换前与变换后的映射函数,k为半偏移距x分量变换前与变换后的映射函数,l为半偏移距y分量变换前与变换后的映射函数。
优选的,所述计算所述雅克比行列式的值,包括:构建变换前当前地震道所覆盖的至少一个多边形;计算所述至少一个多边形的总面积,所述总面积即当前地震道对应的所述雅克比行列式的值。
优选的,根据以下原则构建变换前当前地震道所覆盖的至少一个多边形:对于每个所述多边形,所述多边形内部仅有多边形所在平面内的一个CMP点,且所述多边形内任意一点p到所述CMP点的距离小于点p到其它多边形内的CMP点的距离。
优选的,所述变换后的Kirchhoff积分偏移公式为:
其中,为成像空间坐标;z为成像深度;为成像振幅,为变换后的地震道数据,为变换后的对成像点有贡献的地震道的集合,t为采样时间,为变换后的CMP点的x、y坐标,为变换后的CMP点半偏移距的x、y分量,|J|为雅克比行列式。
根据本发明的另一方面,提出了一种成像振幅一致性校正系统,该系统包括:存储器,存储有计算机可执行指令;处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:将不规则的观测系统变换为规则的观测系统,获得雅克比行列式;计算所述雅克比行列式的值;针对当前地震道,将所述雅克比行列式的值作为加权因子调整Kirchhoff积分偏移公式的积分系数,获得变换后的Kirchhoff积分偏移公式;根据所述变换后的Kirchhoff积分偏移公式对所述当前地震道成像振幅进行调整。
优选的,采用线性坐标变换方法将不规则的观测系统变换为规则的观测系统,获得雅克比行列式,所述雅克比行列式为:
其中,|J|为雅克比行列式,为变换后的CMP点,为变换后的半偏移距,为变换后的CMP点的x坐标,为变换后的CMP点的y坐标,为变换后的半偏移距的x分量,为变换后的半偏移距的y分量,f为CMP点x坐标变换前与变换后的映射函数,g为CMP点y坐标变换前与变换后的映射函数,k为半偏移距x分量变换前与变换后的映射函数,l为半偏移距y分量变换前与变换后的映射函数。
优选的,所述计算所述雅克比行列式的值,包括:构建变换前当前地震道所覆盖的至少一个多边形;计算所述至少一个多边形的总面积,所述总面积即当前地震道对应的所述雅克比行列式的值。
优选的,根据以下原则构建变换前当前地震道所覆盖的至少一个多边形:对于每个所述多边形,所述多边形内部仅有多边形所在平面内的一个CMP点,且所述多边形内任意一点p到所述CMP点的距离小于点p到其它多边形内的CMP点的距离。
优选的,所述变换后的Kirchhoff积分偏移公式为:
其中,为成像空间坐标,z为成像深度,为成像振幅,为变换后的地震道数据,为变换后的对成像点有贡献的地震道的集合,t为采样时间,为变换后的CMP点的x、y坐标,为变换后的CMP点半偏移距的x、y分量,|J|为雅克比行列式。
本发明的有益效果在于:本发明通过将不规则的观测系统变换为规则的观测系统,获取了一个雅克比行列式,并计算雅克比行列式的值,针对当前地震道,用雅克比行列式的值作为加权函数调整Kirchhoff积分偏移公式的积分系数,获得变换后的Kirchhoff积分偏移公式,根据变换后的Kirchhoff积分偏移公式对当前地震道成像振幅进行调整,正是通过雅克比行列式作为加权函数调整偏移过程中每一道的成像权重值,从而调整每一地震道的成像振幅对最终叠加结果的贡献权重,最终消除了观测系统不均匀导致的成像振幅不一致现象,实现了保幅成像目的,有利于后续的AVO分析及其它地震解释工作,且该方法操作方便,效率更高,不会引入过多的计算量。
本发明具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的,或者将在并入本文中的附图和随后的具体实施例中进行详细陈述,这些附图和具体实施例共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施方式进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施方式中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的一个实施例的不规则观测点对应的多边形覆盖范围。
图2示出了根据本发明的一个实施例的一种成像振幅一致性校正方法的流程图。
图3示出了根据本发明的一个实施例的不规则数据、无校正的偏移剖面。
图4示出了根据本发明的一个实施例的一种成像振幅一致性校正方法的偏移剖面。
图5示出了根据本发明的一个实施例的原始规则数据的偏移剖面。
图6示出了根据本发明的一个实施例的不规则数据、无校正的成像道集。
图7示出了根据本发明的一个实施例的一种成像振幅一致性校正方法的成像道集。
图8示出了根据本发明的一个实施例的原始规则数据的成像道集。
具体实施方式
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
根据本发明的一种成像振幅一致性校正方法,包括:将不规则的观测系统变换为规则的观测系统,获得雅克比行列式;计算雅克比行列式的值;针对当前地震道,将雅克比行列式的值作为加权因子调整Kirchhoff积分偏移公式的积分系数,获得变换后的Kirchhoff积分偏移公式;根据变换后的Kirchhoff积分偏移公式对当前地震道成像振幅进行调整。
具体的,通过将不规则的观测系统变换为规则的观测系统,获取了一个雅克比行列式,并计算雅克比行列式的值,针对当前地震道,用雅克比行列式的值作为加权函数调整Kirchhoff积分偏移公式的积分系数,获得变换后的Kirchhoff积分偏移公式,根据变换后的Kirchhoff积分偏移公式对当前地震道成像振幅进行调整。
根据示例性的成像振幅一致性校正方法通过雅克比行列式作为加权函数调整偏移过程中每一道的成像权重值,从而调整每一地震道的成像振幅对最终叠加结果的贡献权重,最终消除了观测系统不均匀导致的成像振幅不一致现象,实现了保幅成像目的,有利于后续的AVO分析及其它地震解释工作,且该方法操作方便,效率更高,不会引入过多的计算量。
作为优选方案,采用线性坐标变换方法将不规则的观测系统变换为规则的观测系统,获得雅克比行列式,雅克比行列式为:
其中,|J|为雅克比行列式,为变换后的CMP点,为变换后的半偏移距,为变换后的CMP点的x坐标,为变换后的CMP点的y坐标,为变换后的半偏移距的x分量,为变换后的半偏移距的y分量,f为CMP点x坐标变换前与变换后的映射函数,g为CMP点y坐标变换前与变换后的映射函数,k为半偏移距x分量变换前与变换后的映射函数,l为半偏移距y分量变换前与变换后的映射函数。
具体的,建立不规则观测系统变换为规则观测系统的坐标变换关系:
假设可以找到如下的线性坐标变换将不规则的观测系统变为规则的观测系统:
其中,变量的上标i为变换前的不规则观测系统,r为变换后的规则观测系统,为变换前的炮点的x坐标,为变换后的炮点的x坐标,为变换后的炮点的y坐标,为变换前的炮点的y坐标,为变换前的检波点的x坐标,为变换后的检波点的x坐标,为变换后的检波点的y坐标,为变换前的检波点的y坐标,f为CMP点x坐标变换前与变换后的映射函数,g为CMP点y坐标变换前与变换后的映射函数,k为半偏移距x分量变换前与变换后的映射函数,l为半偏移距y分量变换前与变换后的映射函数。
其中,为变换前的炮点,为变换前的检波点,为变换前的炮点的x坐标,为变换前的炮点的y坐标,为变换前的检波点的x坐标,为变换前的检波点的y坐标,为变换后的炮点的x坐标,为变换后的炮点的y坐标,为变换后的检波点的x坐标,为变换后的检波点的y坐标,|J|为雅克比行列式,为变换后的炮点,为变换后的检波点,f为CMP点x坐标变换前与变换后的映射函数,g为CMP点y坐标变换前与变换后的映射函数,k为半偏移距x分量变换前与变换后的映射函数,l为半偏移距y分量变换前与变换后的映射函数。
作为优选方案,计算所述雅克比行列式的值,包括:构建变换前当前地震道所覆盖的至少一个多边形;计算至少一个多边形的总面积,总面积即当前地震道对应的所述雅克比行列式的值。
具体的,积分系数表示的几何意义是当前道所覆盖的面积,其中,为变换前的炮点,为变换前的检波点,为变换前的炮点的x坐标,为变换前的炮点的y坐标,为变换前的检波点的x坐标,为变换前的检波点的y坐标,雅克比行列式|J|则表示坐标变换前后的面积比。既然坐标变换后的采样是等间隔的,即是个常数,其中,为变换后的炮点,为变换后的检波点,为变换后的炮点的x坐标,为变换后的炮点的y坐标,为变换后的检波点的x坐标,为变换后的检波点的y坐标,|J|为雅克比行列式,不放设为1,则雅克比行列式|J|的值就是坐标变换前当前道所覆盖的面积。
作为优选方案,根据以下原则构建变换前当前地震道所覆盖的至少一个多边形:对于每个多边形,多边形内部仅有多边形所在平面内的一个CMP点,且多边形内任意一点p到CMP点的距离小于点p到其它多边形内的CMP点的距离。
具体的,假设平面上每一道的CMP点坐标为ai,1≤i≤n,n为总共的CMP点数,该道所覆盖的范围是一个不规则的多边形vi,1≤i≤n,采用如下原则构建:
(1)多边形vi内部有且仅有平面内的一个点ai;
(2)多边形内任意一点p到该CMP点坐标ai的距离小于到其它多边形内的CMP点坐标距离,即满足:||(p,ai)||2<||(p,aj)||2,1≤j≤n,j≠i。
如图1所示,根据上述两个条件构建不规则观测点对应的多边形覆盖范围:点代表不规则观测地震道对应的CMP点;多边形对应该道的有效成像贡献范围,满足上式两个条件所构建的多边形所围成的面积即当前道对应的雅克比行列式值。
作为优选方案,变换后的Kirchhoff积分偏移公式为:
其中,为成像空间坐标,z为成像深度,为成像振幅,为变换后的地震道数据,为变换后的对成像点有贡献的地震道的集合,t为采样时间,为变换后的CMP点的x、y坐标,为变换后的CMP点半偏移距的x、y分量,|J|为雅克比行列式。
具体的,常规Kirchhoff积分法偏移(包括时间偏移和深度偏移)是通过如下方程实现的:
其中,为成像空间坐标,z为成像深度,为成像振幅,函数为整个地震数据,为地震道空间坐标(xs,ys,xg,yg),t为地震波旅行时,即采样时间,A为与几何扩散及传播方向有关的校正系数,为对成像点有贡献的地震道的集合,可以理解成偏移孔径内地震道的集合,分别为炮点和检波点的坐标。
其中,为成像空间坐标,z为成像深度,为成像振幅,函数为整个地震数据,为地震道空间坐标(xs,ys,xg,yg),t为地震波旅行时,即采样时间,A为与几何扩散及传播方向有关的校正系数,为对成像点有贡献的地震道的集合,可以理解成偏移孔径内地震道的集合,为CMP点的x、y坐标,为半偏移距的x、y分量,即
常规Kirchhoff积分法偏移假设观测系统采样均匀,认为每一地震道所覆盖的地表面积是相等的,即是个常数,其中,为变换前的炮点,为变换前的检波点,为变换前的炮点的x坐标,为变换前的炮点的y坐标,为变换前的检波点的x坐标,为变换前的检波点的y坐标,因此对公式(5)的计算一般不考虑积分系数
进一步地,将Kirchhoff积分方程(5)写成如下形式:
其中,为成像空间坐标,z为成像深度,为成像振幅,为变换前的地震道数据,为变换前的地震道空间坐标(xs,ys,xg,yg),为变换前的CMP点的x、y坐标,为变换前的半偏移距的x、y分量,为变换后的地震道数据,为变换后的对成像点有贡献的地震道的集合,t为采样时间,为变换后的CMP点的x、y坐标,为变换后的CMP点半偏移距的x、y分量,|J|为雅克比行列式。
根据本发明的一种成像振幅一致性校正系统,该系统包括:存储器,存储有计算机可执行指令;处理器,处理器运行存储器中的计算机可执行指令,执行以下步骤:将不规则的观测系统变换为规则的观测系统,获得雅克比行列式;计算雅克比行列式的值;针对当前地震道,将雅克比行列式的值作为加权因子调整Kirchhoff积分偏移公式的积分系数,获得变换后的Kirchhoff积分偏移公式;根据变换后的Kirchhoff积分偏移公式对当前地震道成像振幅进行调整。
具体的,通过将不规则的观测系统变换为规则的观测系统,获取了一个雅克比行列式,并计算雅克比行列式的值,针对当前地震道,用雅克比行列式的值作为加权函数调整Kirchhoff积分偏移公式的积分系数,获得变换后的Kirchhoff积分偏移公式,根据变换后的Kirchhoff积分偏移公式对当前地震道成像振幅进行调整。
根据示例性的成像振幅一致性校正系统通过雅克比行列式作为加权函数调整偏移过程中每一道的成像权重值,从而调整每一地震道的成像振幅对最终叠加结果的贡献权重,最终消除了观测系统不均匀导致的成像振幅不一致现象,实现了保幅成像目的,有利于后续的AVO分析及其它地震解释工作,且该方法操作方便,效率更高,不会引入过多的计算量。
作为优选方案,采用线性坐标变换方法将不规则的观测系统变换为规则的观测系统,获得雅克比行列式,雅克比行列式为:
其中,|J|为雅克比行列式,为变换后的CMP点,为变换后的半偏移距,为变换后的CMP点的x坐标,为变换后的CMP点的y坐标,为变换后的半偏移距的x分量,为变换后的半偏移距的y分量,f为CMP点x坐标变换前与变换后的映射函数,g为CMP点y坐标变换前与变换后的映射函数,k为半偏移距x分量变换前与变换后的映射函数,l为半偏移距y分量变换前与变换后的映射函数。
具体的,建立不规则观测系统变换为规则观测系统的坐标变换关系:
假设可以找到如下的线性坐标变换将不规则的观测系统变为规则的观测系统:
其中,变量的上标i为变换前的不规则观测系统,r为变换后的规则观测系统,为变换前的炮点的x坐标,为变换后的炮点的x坐标,为变换后的炮点的y坐标,为变换前的炮点的y坐标,为变换前的检波点的x坐标,为变换后的检波点的x坐标,为变换后的检波点的y坐标,为变换前的检波点的y坐标,f为CMP点x坐标变换前与变换后的映射函数,g为CMP点y坐标变换前与变换后的映射函数,k为半偏移距x分量变换前与变换后的映射函数,l为半偏移距y分量变换前与变换后的映射函数。
其中,为变换前的炮点,为变换前的检波点,为变换前的炮点的x坐标,为变换前的炮点的y坐标,为变换前的检波点的x坐标,为变换前的检波点的y坐标,为变换后的炮点的x坐标,为变换后的炮点的y坐标,为变换后的检波点的x坐标,为变换后的检波点的y坐标,|J|为雅克比行列式,为变换后的炮点,为变换后的检波点,f为CMP点x坐标变换前与变换后的映射函数,g为CMP点y坐标变换前与变换后的映射函数,k为半偏移距x分量变换前与变换后的映射函数,l为半偏移距y分量变换前与变换后的映射函数。
作为优选方案,计算所述雅克比行列式的值,包括:构建变换前当前地震道所覆盖的至少一个多边形;计算至少一个多边形的总面积,总面积即当前地震道对应的所述雅克比行列式的值。
具体的,积分系数表示的几何意义是当前道所覆盖的面积,其中,为变换前的炮点,为变换前的检波点,为变换前的炮点的x坐标,为变换前的炮点的y坐标,为变换前的检波点的x坐标,为变换前的检波点的y坐标,雅克比行列式|J|则表示坐标变换前后的面积比。既然坐标变换后的采样是等间隔的,即是个常数,其中,为变换后的炮点,为变换后的检波点,为变换后的炮点的x坐标,为变换后的炮点的y坐标,为变换后的检波点的x坐标,为变换后的检波点的y坐标,|J|为雅克比行列式,不放设为1,则雅克比行列式|J|的值就是坐标变换前当前道所覆盖的面积。
作为优选方案,根据以下原则构建变换前当前地震道所覆盖的至少一个多边形:对于每个多边形,多边形内部仅有多边形所在平面内的一个CMP点,且多边形内任意一点p到CMP点的距离小于点p到其它多边形内的CMP点的距离。
具体的,假设平面上每一道的CMP点坐标为ai,1≤i≤n,n为总共的CMP点数,该道所覆盖的范围是一个不规则的多边形vi,1≤i≤n,采用如下原则构建:
(1)多边形vi内部有且仅有平面内的一个点ai;
(2)多边形内任意一点p到该CMP点坐标ai的距离小于到其它多边形内的CMP点坐标距离,即满足:||(p,ai)||2<||(p,aj)||2,1≤j≤n,j≠i。
如图1所示,根据上述两个条件构建不规则观测点对应的多边形覆盖范围:点代表不规则观测地震道对应的CMP点;多边形对应该道的有效成像贡献范围,满足上式两个条件所构建的多边形所围成的面积即当前道对应的雅克比行列式值。
作为优选方案,变换后的Kirchhoff积分偏移公式为:
其中,为成像空间坐标,z为成像深度,为成像振幅,为变换后的地震道数据,为变换后的对成像点有贡献的地震道的集合,t为采样时间,为变换后的CMP点的x、y坐标,为变换后的CMP点半偏移距的x、y分量,|J|为雅克比行列式。
具体的,常规Kirchhoff积分法偏移(包括时间偏移和深度偏移)是通过如下方程实现的:
其中,为成像空间坐标,z为成像深度,为成像振幅,函数为整个地震数据,为地震道空间坐标(xs,ys,xg,yg),t为地震波旅行时,即采样时间,A为与几何扩散及传播方向有关的校正系数,为对成像点有贡献的地震道的集合,可以理解成偏移孔径内地震道的集合,分别为炮点和检波点的坐标。
其中,为成像空间坐标,z为成像深度,为成像振幅,函数为整个地震数据,为地震道空间坐标(xs,ys,xg,yg),t为地震波旅行时,即采样时间,A为与几何扩散及传播方向有关的校正系数,为对成像点有贡献的地震道的集合,可以理解成偏移孔径内地震道的集合,为CMP点的x、y坐标,为半偏移距的x、y分量,即
常规Kirchhoff积分法偏移假设观测系统采样均匀,认为每一地震道所覆盖的地表面积是相等的,即是个常数,其中,为变换前的炮点,为变换前的检波点,为变换前的炮点的x坐标,为变换前的炮点的y坐标,为变换前的检波点的x坐标,为变换前的检波点的y坐标,因此对公式(5)的计算一般不考虑积分系数
进一步地,将Kirchhoff积分方程(5)写成如下形式:
其中,为成像空间坐标,z为成像深度,为成像振幅,为变换前的地震道数据,为变换前的地震道空间坐标(xs,ys,xg,yg),为变换前的CMP点的x、y坐标,为变换前的半偏移距的x、y分量,为变换后的地震道数据,为变换后的对成像点有贡献的地震道的集合,t为采样时间,为变换后的CMP点的x、y坐标,为变换后的CMP点半偏移距的x、y分量,|J|为雅克比行列式。
实施例
图2示出了根据本发明的一个实施例的一种成像振幅一致性校正方法的流程图。
如图2所示,一种成像振幅一致性校正方法,包括:
S102:将不规则的观测系统变换为规则的观测系统,获得雅克比行列式;
其中,采用线性坐标变换方法将不规则的观测系统变换为规则的观测系统,获得雅克比行列式,雅克比行列式为:
其中,|J|为雅克比行列式,为变换后的CMP点,为变换后的半偏移距,为变换后的CMP点的x坐标,为变换后的CMP点的y坐标,为变换后的半偏移距的x分量,为变换后的半偏移距的y分量,f为CMP点x坐标变换前与变换后的映射函数,g为CMP点y坐标变换前与变换后的映射函数,k为半偏移距x分量变换前与变换后的映射函数,l为半偏移距y分量变换前与变换后的映射函数。
具体的,建立不规则观测系统变换为规则观测系统的坐标变换关系:
假设可以找到如下的线性坐标变换将不规则的观测系统变为规则的观测系统:
其中,变量的上标i为变换前的不规则观测系统,r为变换后的规则观测系统,为变换前的炮点的x坐标,为变换后的炮点的x坐标,为变换后的炮点的y坐标,为变换前的炮点的y坐标,为变换前的检波点的x坐标,为变换后的检波点的x坐标,为变换后的检波点的y坐标,为变换前的检波点的y坐标,f为CMP点x坐标变换前与变换后的映射函数,g为CMP点y坐标变换前与变换后的映射函数,k为半偏移距x分量变换前与变换后的映射函数,l为半偏移距y分量变换前与变换后的映射函数。
其中,为变换前的炮点,为变换前的检波点,为变换前的炮点的x坐标,为变换前的炮点的y坐标,为变换前的检波点的x坐标,为变换前的检波点的y坐标,为变换后的炮点的x坐标,为变换后的炮点的y坐标,为变换后的检波点的x坐标,为变换后的检波点的y坐标,|J|为雅克比行列式,为变换后的炮点,为变换后的检波点,f为CMP点x坐标变换前与变换后的映射函数,g为CMP点y坐标变换前与变换后的映射函数,k为半偏移距x分量变换前与变换后的映射函数,l为半偏移距y分量变换前与变换后的映射函数。
S104:计算雅克比行列式的值;
其中,计算所述雅克比行列式的值,包括:构建变换前当前地震道所覆盖的至少一个多边形;计算至少一个多边形的总面积,总面积即当前地震道对应的所述雅克比行列式的值。
具体的,积分系数表示的几何意义是当前道所覆盖的面积,其中,为变换前的炮点,为变换前的检波点,为变换前的炮点的x坐标,为变换前的炮点的y坐标,为变换前的检波点的x坐标,为变换前的检波点的y坐标,雅克比行列式|J|则表示坐标变换前后的面积比。既然坐标变换后的采样是等间隔的,即是个常数,其中,为变换后的炮点,为变换后的检波点,为变换后的炮点的x坐标,为变换后的炮点的y坐标,为变换后的检波点的x坐标,为变换后的检波点的y坐标,|J|为雅克比行列式,不放设为1,则雅克比行列式|J|的值就是坐标变换前当前道所覆盖的面积。
其中,根据以下原则构建变换前当前地震道所覆盖的至少一个多边形:对于每个多边形,多边形内部仅有多边形所在平面内的一个CMP点,且多边形内任意一点p到CMP点的距离小于点p到其它多边形内的CMP点的距离。
具体的,假设平面上每一道的CMP点坐标为ai,1≤i≤n,n为总共的CMP点数,该道所覆盖的范围是一个不规则的多边形vi,1≤i≤n,采用如下原则构建:
(1)多边形vi内部有且仅有平面内的一个点ai;
(2)多边形内任意一点p到该CMP点坐标ai的距离小于到其它多边形内的CMP点坐标距离,即满足:||(p,ai)||2<||(p,aj)||2,1≤j≤n,j≠i。
如图1所示,根据上述两个条件构建不规则观测点对应的多边形覆盖范围:点代表不规则观测地震道对应的CMP点;多边形对应该道的有效成像贡献范围,满足上式两个条件所构建的多边形所围成的面积即当前道对应的雅克比行列式值。
S106:针对当前地震道,将雅克比行列式的值作为加权因子调整Kirchhoff积分偏移公式的积分系数,获得变换后的Kirchhoff积分偏移公式;
其中,变换后的Kirchhoff积分偏移公式为:
其中,为成像空间坐标,z为成像深度,为成像振幅,为变换后的地震道数据,为变换后的对成像点有贡献的地震道的集合,t为采样时间,为变换后的CMP点的x、y坐标,为变换后的CMP点半偏移距的x、y分量,|J|为雅克比行列式。
具体的,常规Kirchhoff积分法偏移(包括时间偏移和深度偏移)是通过如下方程实现的:
其中,为成像空间坐标,z为成像深度,为成像振幅,函数为整个地震数据,为地震道空间坐标(xs,ys,xg,yg),t为地震波旅行时,即采样时间,A为与几何扩散及传播方向有关的校正系数,为对成像点有贡献的地震道的集合,可以理解成偏移孔径内地震道的集合,分别为炮点和检波点的坐标。
其中,为成像空间坐标,z为成像深度,为成像振幅,函数为整个地震数据,为地震道空间坐标(xs,ys,xg,yg),t为地震波旅行时,即采样时间,A为与几何扩散及传播方向有关的校正系数,为对成像点有贡献的地震道的集合,可以理解成偏移孔径内地震道的集合,为CMP点的x、y坐标,为半偏移距的x、y分量,即
常规Kirchhoff积分法偏移假设观测系统采样均匀,认为每一地震道所覆盖的地表面积是相等的,即是个常数,其中,为变换前的炮点,为变换前的检波点,为变换前的炮点的x坐标,为变换前的炮点的y坐标,为变换前的检波点的x坐标,为变换前的检波点的y坐标,因此对公式(5)的计算一般不考虑积分系数
进一步地,将Kirchhoff积分方程(5)写成如下形式:
其中,为成像空间坐标,z为成像深度,为成像振幅,为变换前的地震道数据,为变换前的地震道空间坐标(xs,ys,xg,yg),为变换前的CMP点的x、y坐标,为变换前的半偏移距的x、y分量,为变换后的地震道数据,为变换后的对成像点有贡献的地震道的集合,t为采样时间,为变换后的CMP点的x、y坐标,为变换后的CMP点半偏移距的x、y分量,|J|为雅克比行列式。
S108:根据变换后的Kirchhoff积分偏移公式对当前地震道成像振幅进行调整。
图3示出了根据本发明的一个实施例的不规则数据、无校正的偏移剖面。图4示出了根据本发明的一个实施例的一种成像振幅一致性校正方法的偏移剖面。图5示出了根据本发明的一个实施例的原始规则数据的偏移剖面。
如图3、图4和图5所示,采用一个二维正演数据:共有2200个CMP道集,每个道集内有60个地震道,且偏移距分布均匀。图3的不规则数据体是通过在原始数据体上进行抽道形成的:首先在2200个CMP道集中随机挑选800个CMP道集,然后在每个被挑中的CMP道集中随机去掉2/3道。将图3和图4的数据与图5的原始规则数据做对比,可以看出图3的数据振幅不一致,图4的数据振幅一致。在图3-图5,横坐标为CMP,纵坐标为时间(s)
图6示出了根据本发明的一个实施例的不规则数据、无校正的成像道集。图7示出了根据本发明的一个实施例的一种成像振幅一致性校正方法的成像道集。图8示出了根据本发明的一个实施例的原始规则数据的成像道集。
如图6、图7和图8所示,采用一个二维正演数据:共有2200个CMP道集,每个道集内有60个地震道,且偏移距分布均匀。图6的不规则数据体是通过在原始数据体上进行抽道形成的:首先在2200个CMP道集中随机挑选800个CMP道集,然后在每个被挑中的CMP道集中随机去掉2/3道。将图6和图7的数据与图8的原始规则数据做对比,可以看出图6的数据振幅不一致,图7的数据振幅一致。在图6-图8,横坐标为道数,纵坐标为时间(s)。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (10)
1.一种成像振幅一致性校正方法,其特征在于,包括:
将不规则的观测系统变换为规则的观测系统,获得雅克比行列式;
计算所述雅克比行列式的值;
针对当前地震道,将所述雅克比行列式的值作为加权因子调整Kirchhoff积分偏移公式的积分系数,获得变换后的Kirchhoff积分偏移公式;
根据所述变换后的Kirchhoff积分偏移公式对所述当前地震道成像振幅进行调整。
3.根据权利要求1所述的成像振幅一致性校正方法,其特征在于,所述计算所述雅克比行列式的值,包括:
构建变换前当前地震道所覆盖的至少一个多边形;
计算所述至少一个多边形的总面积,所述总面积即当前地震道对应的所述雅克比行列式的值。
4.根据权利要求3所述的成像振幅一致性校正方法,其特征在于,根据以下原则构建变换前当前地震道所覆盖的至少一个多边形:
对于每个所述多边形,所述多边形内部仅有多边形所在平面内的一个CMP点,且所述多边形内任意一点p到所述CMP点的距离小于点p到其它多边形内的CMP点的距离。
6.一种成像振幅一致性校正系统,其特征在于,该系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
将不规则的观测系统变换为规则的观测系统,获得雅克比行列式;
计算所述雅克比行列式的值;
针对当前地震道,将所述雅克比行列式的值作为加权因子调整Kirchhoff积分偏移公式的积分系数,获得变换后的Kirchhoff积分偏移公式;
根据所述变换后的Kirchhoff积分偏移公式对所述当前地震道成像振幅进行调整。
8.根据权利要求6所述的成像振幅一致性校正系统,其特征在于,所述计算所述雅克比行列式的值,包括:
构建变换前当前地震道所覆盖的至少一个多边形;
计算所述至少一个多边形的总面积,所述总面积即当前地震道对应的所述雅克比行列式的值。
9.根据权利要求8所述的成像振幅一致性校正系统,其特征在于,根据以下原则构建变换前当前地震道所覆盖的至少一个多边形:
对于每个所述多边形,所述多边形内部仅有多边形所在平面内的一个CMP点,且所述多边形内任意一点p到所述CMP点的距离小于点p到其它多边形内的CMP点的距离。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811033328.8A CN110879414A (zh) | 2018-09-05 | 2018-09-05 | 一种成像振幅一致性校正方法及校正系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811033328.8A CN110879414A (zh) | 2018-09-05 | 2018-09-05 | 一种成像振幅一致性校正方法及校正系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110879414A true CN110879414A (zh) | 2020-03-13 |
Family
ID=69727801
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811033328.8A Pending CN110879414A (zh) | 2018-09-05 | 2018-09-05 | 一种成像振幅一致性校正方法及校正系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110879414A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050143923A1 (en) * | 2003-12-19 | 2005-06-30 | Henk Keers | Processing seismic data |
CN102914791A (zh) * | 2011-08-05 | 2013-02-06 | 中国石油天然气集团公司 | 一种起伏地表地震数据处理的克希霍夫叠前时间偏移方法 |
CN104142518A (zh) * | 2013-10-29 | 2014-11-12 | 中国石油化工股份有限公司 | 一种地震观测系统叠前时间偏移响应的分析方法 |
CN104407380A (zh) * | 2014-11-27 | 2015-03-11 | 中国石油化工股份有限公司 | 一种处理叠前偏移距分组地震数据的方法 |
-
2018
- 2018-09-05 CN CN201811033328.8A patent/CN110879414A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050143923A1 (en) * | 2003-12-19 | 2005-06-30 | Henk Keers | Processing seismic data |
CN102914791A (zh) * | 2011-08-05 | 2013-02-06 | 中国石油天然气集团公司 | 一种起伏地表地震数据处理的克希霍夫叠前时间偏移方法 |
CN104142518A (zh) * | 2013-10-29 | 2014-11-12 | 中国石油化工股份有限公司 | 一种地震观测系统叠前时间偏移响应的分析方法 |
CN104407380A (zh) * | 2014-11-27 | 2015-03-11 | 中国石油化工股份有限公司 | 一种处理叠前偏移距分组地震数据的方法 |
Non-Patent Citations (6)
Title |
---|
ANAT CANNING ET AL.: "Reducing 3-D acquisition footprint for 3-D DMO and 3-D prestack migration", 《GEOPHYSICS》 * |
HERMAN JARAMILLO ET AL.: "Practical aspects of Voronoi-based area weights for Kirchhoff migration.", 《SEG/HOUSTON 2005 ANNUAL MEETING》 * |
NORMAN BLEISTEIN ET AL.: "On the imaging of reflectors in the earth", 《GEOPHYSICS》 * |
XIONGWEN WANG ET AL.: "Reducing the acquisition footprint using coordinate transformation", 《SEG HOUSTON 2009 INTERNATIONAL EXPOSITION AND ANNUAL MEETING》 * |
中国石化石油物探技术研究院: "《油气地球物理技术新进展—第79届SEG年会论文概要》", 31 March 2011, 石油工业出版社 * |
林琳等: "《地理信息系统基础及应用》", 31 May 2018, 中国矿业大学出版社 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105005042B (zh) | 一种探地雷达地下目标定位方法 | |
KR101475382B1 (ko) | 광학적 3차원 측량의 자기 적응 윈도우 푸리에 위상추출방법 | |
CN110866545A (zh) | 一种探地雷达资料中管线目标的自动识别方法及系统 | |
Terrasse et al. | Application of the curvelet transform for clutter and noise removal in GPR data | |
CN104715491B (zh) | 一种基于一维灰度矩的亚像素边缘检测方法 | |
CN102750705A (zh) | 基于图像融合的光学遥感图像变化检测 | |
CN112258542A (zh) | 红外图像目标检测方法、装置、电子设备及存储介质 | |
CN110764087A (zh) | 一种基于干涉成像高度计的海面风向反加权反演方法 | |
CN110907995A (zh) | 井中vsp地震数据的逆时偏移方法及装置 | |
CN107610174B (zh) | 一种鲁棒的基于深度信息的平面检测方法及系统 | |
Zhang et al. | An improved BM3D algorithm based on anisotropic diffusion equation | |
CN110879414A (zh) | 一种成像振幅一致性校正方法及校正系统 | |
CN109461171A (zh) | 基于多通道改进DoG滤波的红外弱小目标检测算法 | |
CN111965645A (zh) | 一种顾及几何约束的多基线植被高度反演方法及装置 | |
CN115507752B (zh) | 基于平行环境元素的单目视觉测距方法和系统 | |
CN111461999A (zh) | 一种基于超像素相似性测量的sar图像相干斑抑制方法 | |
CN108428216B (zh) | 基于散布矩阵特征的二阶偏微分方程遥感图像去噪方法 | |
CN111126508A (zh) | 一种基于hopc改进的异源图像匹配方法 | |
CN115598704A (zh) | 一种基于最小二乘逆时偏移生成保幅角道集的方法、设备及可读存储介质 | |
CN110223250B (zh) | 基于单应变换的sar几何校正方法 | |
Xu et al. | Evaluation of a canonical image representation for sidescan sonar | |
CN109975397B (zh) | 基于多频涡流复信号的传热管损伤信息高保真提取方法 | |
CN109557541B (zh) | 一种全息穿透成像雷达极坐标数据处理方法 | |
CN109655888B (zh) | 地震数据处理中光滑浮动基准面的定量选择方法及系统 | |
CN111751874B (zh) | 一种变偏移距vsp叠后变覆盖次数校正方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200313 |