CN105717538B - 起伏地表地震数据偏移基准面转换方法及装置 - Google Patents
起伏地表地震数据偏移基准面转换方法及装置 Download PDFInfo
- Publication number
- CN105717538B CN105717538B CN201410720019.3A CN201410720019A CN105717538B CN 105717538 B CN105717538 B CN 105717538B CN 201410720019 A CN201410720019 A CN 201410720019A CN 105717538 B CN105717538 B CN 105717538B
- Authority
- CN
- China
- Prior art keywords
- cmp
- data
- correcting value
- coordinate
- geophone station
- 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.)
- Active
Links
Abstract
本发明提供是一种起伏地表地震数据偏移基准面转换方法及装置，获取起伏地表的偏移基准面的数据信息，数据信息包括CMP点的X坐标和Y坐标及各CMP点的偏移基准面的第一校正量，获取炮点和检波点的坐标信息，根据炮点和检波点以及CMP点的物理位置，将与炮点和检波点相邻的CMP点对应的第一校正量插值到炮点和检波点上，得到炮点和检波点的偏移基准面的第二校正量，根据CMP点的第一校正量与炮点和检波点的第二校正量，计算出CMP局部平面与起伏地表的偏移基准面之间的转换校正量，使用转换校正量校正待转换数据得到目标数据。本发明基于转换校正量在不同系统或偏移方法中进行地震数据的转换，使用转换后的数据能获得正确的成像效果。
Description
技术领域
本发明涉及地震资料数字处理领域，尤其涉及一种起伏地表地震数据偏移基准面转换方法及装置。
背景技术
目前陆上反射地震勘探领域中获取最终地下成像的最佳方法是基于起伏地表条件的叠前深度偏移处理。由于地表的高低起伏需要根据地表及近地表情况选用起伏地表的偏移基准面。目前业界主要使用克希霍夫积分偏移、炮域单程波偏移和炮域逆时偏移等方法，在起伏地表偏移基准面上实现反射波地震资料叠前深度偏移。
实现这几种方法的地震数据处理软件系统有多套，不同方法与系统对输入数据的要求不同。现有的地震数据处理系统没有提供可适用不同软件系统的输入数据准备要求的数据转换方法。不同系统或者不同偏移方法之间进行起伏地表地震数据转换时容易出错，出现偏移后成像深度不一致、成像效果变差、存在偏移划弧等问题。
发明内容
本发明提供一种起伏地表地震数据偏移基准面转换方法及装置，用于解决现有不同系统或偏移方式进行地震数据转换时存在数据出错，导致成像深度不一致、成像效果较差以及偏移划弧等问题。
为了实现上述目的，本发明提供的一种起伏地表地震数据偏移基准面转换方法，包括：
获取起伏地表的偏移基准面的数据信息，所述数据信息包括各共中心点CMP点的X坐标和Y坐标以及各CMP点的偏移基准面的第一校正量；
获取炮点和检波点的坐标信息，所述坐标信息包括炮点的X坐标和Y坐标，以及检波点的X坐标和Y坐标；
根据所述炮点和所述检波点以及所述CMP点的物理位置，将与所述炮点和所述检波点相邻的所述CMP点对应的第一校正量插值到所述炮点和所述检波点上，得到所述炮点的偏移基准面的第二校正量和所述检波点的偏移基准面的第二校正量；其中所述物理位置是由X坐标和Y坐标确定的；
根据所述CMP点的所述第一校正量与所述炮点的所述第二校正量和所述检波点的所述第二校正量，计算得到CMP局部平面与起伏地表的偏移基准面之间的转换校正量；
使用所述转换校正量对待转换数据进行校正得到目标数据。
为了实现上述目的，本发明提供的一种起伏地表地震数据偏移基准面转换装置，包括：
第一获取模块，用于获取起伏地表的偏移基准面的数据信息，所述数据信息包括各共中心点CMP点的X坐标和Y坐标以及各CMP点的偏移基准面的第一校正量；
第二获取模块，用于获取炮点和检波点的坐标信息，所述坐标信息包括炮点的X坐标和Y坐标，以及检波点的X坐标和Y坐标；
插值模块，用于根据所述炮点和所述检波点以及所述CMP点的物理位置，将与所述炮点和所述检波点相邻的所述CMP点对应的第一校正量插值到所述炮点和所述检波点上，得到所述炮点的偏移基准面的第二校正量和所述检波点的偏移基准面的第二校正量；其中所述物理位置是由X坐标和Y坐标确定的；
计算模块，用于根据所述CMP点的所述第一校正量与所述炮点的所述第二校正量和所述检波点的所述第二校正量，计算得到CMP局部平面与起伏地表的偏移基准面之间的转换校正量；
校正转换模块，用于使用所述转换校正量对待转换数据进行校正得到目标数据。
本发明提供的起伏地表地震数据偏移基准面转换方法及装置，获取起伏地表的偏移基准面的数据信息，所述数据信息包括各共中心点CMP点的X坐标和Y坐标以及各CMP点的偏移基准面的第一校正量，获取炮点和检波点的坐标信息，所述坐标信息包括炮点的X坐标和Y坐标，以及检波点的X坐标和Y坐标，根据所述炮点和所述检波点以及所述CMP点的物理位置，将与所述炮点和所述检波点相邻的所述CMP点对应的第一校正量插值到所述炮点和所述检波点上，得到所述炮点的偏移基准面的第二校正量和所述检波点的偏移基准面的第二校正量；其中所述物理位置是由X坐标和Y坐标确定的，根据所述CMP点的所述第一校正量与所述炮点的所述第二校正量和所述检波点的所述第二校正量，计算得到CMP局部平面与起伏地表的偏移基准面之间的转换校正量，使用所述转换校正量对待转换数据进行校正得到目标数据。本发明基于CMP点的第一校正量获取炮点和检波点的第二校正量，从而得到获取CMP局部平面与偏移基准面之间的转换校正量，基于该转换校正量能在不同系统或偏移方法中进行地震数据的转换，使得转换数据过程简单准确，使用转换后的数据能获得正确的成像效果，而且成像深度能够保持一致。
附图说明
图1为本发明实施例一提供的一种起伏地表地震数据偏移基准面转换方法的流程示意图；
图2为起伏地表偏移基准面选取示意图；
图3为文本格式的起伏地表的偏移基准面数据信息；
图4为文本格式的炮点的坐标信息；
图5为插值后炮点的偏移基准面第二校正量；
图6为本发明实施例提供的地震数据校正到CMP局部平面与偏移基准面的区别；
图7为本发明实施例二提供的一种起伏地表地震数据偏移基准面转换装置的结构示意图；
图8为本发明实施例三提供的另一种起伏地表地震数据偏移基准面转换装置的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。
实施例一
图1为本发明实施例一提供的一种起伏地表地震数据偏移基准面转换方法的流程示意图。如图1所示，该起伏地表地震数据偏移基准面转换方法包括以下步骤：
101、获取起伏地表的偏移基准面的数据信息，所述数据信息包括各共中心点CMP点的X坐标和Y坐标以及各CMP点的偏移基准面的第一校正量。
图2为起伏地表偏移基准面选取示意图。如图2所示，真实地表曲线1地表起伏高低差异较大，为了对地表情况进行模拟，对真是地表进行平滑处理，得到起伏差异较小能基本反映出真实地表的偏移基准面曲线2。在地震勘查中，为了便于计算还会设置一个地震参考基准面曲线3。
基于选取的起伏地表的偏移基准面，对偏移基准面进行数据查询，得到该偏移基准面的数据信息，其中，该数据信息包括各共中心点(Common Mid Point，简称CMP)点的X坐和Y坐标，以及各CMP点的偏移基准面的第一校正量。其中，获取各CMP点的偏移基准面的第一校正量的过程包括：将各CMP点偏移基准面的高程与对应地震参考基准面高程求差值，使用该差值与替换速度做比值，该比值的2倍作为CMP点的第一校正量。
可选地，所述数据信息中还可以包括各CMP点的线号以及点号，其中线号和点号是为了使用的需要对CMP点做的二维编号。
例如，获取的起伏地表的偏移基准面的数据信息可以采用文本格式输出，如图3所示。
102、获取炮点和检波点的坐标信息，所述坐标信息包括炮点的X坐标和Y坐标，以及检波点的X坐标和Y坐标。
本实施例中，在进行地震勘查时，获取了一些炮点的大地坐标值，即炮点的X坐标和Y坐标。同样地，获取了一些检波点的大地坐标值，即检波点的X坐标和Y坐标。可选地，坐标信息中还可以包括各炮点以及检波点的标志号，例如顺序编写的序号。
例如，获取的炮点的坐标信息可以采用文本格式输出，如图4所示。
103、根据炮点和检波点以及CMP点的物理位置，将与炮点和检波点相邻的CMP点对应的第一校正量插值到炮点和检波点上，得到炮点的偏移基准面的第二校正量和检波点的偏移基准面的第二校正量。
其中，所述物理位置是由X坐标和Y坐标确定的。具体地，炮点的物理位置是由炮点的X坐标和Y坐标确定的，检波点的物理位置是由检波点的X坐标和Y坐标确定的，CMP点的物理位置是由CMP点的X坐标和Y坐标确定的。
在得知炮点的物理位置以及CMP点的物理位置后，能够根据炮点的物理位置获取与该物理位置相邻的CMP点。具体地，根据炮点的X坐标和Y坐标查询所述数据信息，获取与炮点相邻的CMP点。在获取到与炮点相邻的CMP点后，根据与炮点相邻的CMP点对应的第一校正量得到炮点的偏移基准面的第二校正量。即将与炮点相邻的CMP点对应的第一校正量插值到该炮点上，得到该炮点的偏移基准面的第二校正量。可选地，将与炮点相邻的CMP点对应的第一校正量进行加权平均，得到炮点的第二校正量。一般与炮点相邻的CMP点的个数为4个，对这4个CMP点的第一校正量进行加权平均，就可以得到炮点的第二校正量。
类似地，在得知检波点的物理位置以及CMP点的物理位置后，能够根据检波点的物理位置获取与该物理位置相邻的CMP点。具体地，根据检波点的X坐标和Y坐标查询所述数据信息，获取与检波点相邻的CMP点。在获取到与该检波点相邻的CMP点后，根据与检波点相邻的CMP点对应的第一校正量得到检波点的偏移基准面的第二校正量。即将与该检波点相邻的CMP点对应的第一校正量插值到该检波点上，得到该检波点的偏移基准面的第二校正量。可选地，将与检波点相邻的CMP点对应的第一校正量进行加权平均，得到检波点的第二校正量。一般与检波点相邻的CMP点的个数为4个，对这4个CMP点的第一校正量进行加权平均，就可以得到检波点的第二校正量。
例如，插值后炮的偏移基准面第二校正量采用文本格式输出，如图5所示。
104、根据CMP点的偏移基准面的第一校正量与炮点的偏移基准面的第二校正量和检波点的偏移基准面的第二校正量，计算得到CMP局部平面与起伏地表的偏移基准面之间的转换校正量。
具体地，在每一个地震数据道的道头信息内，将炮点的偏移基准面的第二校正量和检波点的偏移基准面的第二校正量相加，然后将相加的和一半与CMP点的偏移基准面的第一校正量做差，得到的差值为CMP局部平面与偏移基准面之间的转换校正量。
105、使用转换校正量对待转换数据进行校正得到目标数据。
一方面，如果待转换数据为CMP道集数据，所述使用所述转换校正量对待转换数据进行校正得到目标数据，包括：
将待转换CMP道集数据进行重排得到炮集数据，再使用获取到的转换校正量对重排得到的炮集数据进行校正，得到校正到偏移基准面上的目标数据。本实施例中，所述目标数据为校正到所述偏移基准面上的炮集数据。
另一方面，如果待转换数据为炮集数据，所述使用所述转换校正量对待转换数据进行校正得到目标数据，包括：
将待转换的炮集数据进行重排得到CMP道集数据，再使用转换校正量对CMP道集数据进行反向校正，得到校正到CMP局部平面上的目标数据。本实施例中，所述目标数据为校正到CMP局部平面上的CMP道集数据。
图6为本发明实施例提供的数据校正到CMP局部平面与偏移基准面的区别。图6中包括地震参考基准面、偏移基准面以及CMP局部平面。基于地震参考基准面的高程、偏移基准面的高程以及替换速度，能够获取炮点的偏移基准面的第二校正量，以及检波点的偏移基准面的第二校正量。基于地震参考基准面的高程、CMP局部平面的高程以及替换速度，能够获取CMP点的偏移基准面的第一校正量。
本实施例提供的起伏地表地震数据偏移基准面转换方法，获取起伏地表的偏移基准面的数据信息，其中数据信息包括各共中心点CMP点的X坐标和Y坐标以及各CMP点的偏移基准面的第一校正量，获取炮点和检波点的坐标信息，坐标信息包括炮点的X坐标和Y坐标以及检波点的X坐标和Y坐标，根据炮点和检波点以及CMP点的物理位置，将与炮点和检波点相邻的CMP点对应的第一校正量插值到炮点和检波点上，得到炮点的偏移基准面的第二校正量和检波点的偏移基准面的第二校正量，根据CMP点的第一校正量与炮点的第二校正量和检波点的第二校正量，计算得到CMP局部平面与起伏地表的偏移基准面之间的转换校正量，使用转换校正量对待转换数据进行校正得到目标数据。本实施例基于CMP点的第一校正量获取炮点和检波点的第二校正量，从而获取CMP局部平面与起伏地表的偏移基准面之间的转换校正量，基于该转换校正量能在不同系统或偏移方法中进行地震数据的转换，使得转换数据过程简单准确，使用转换后的数据能获得正确的成像效果，并且成像深度保持一致。
实施例二
图7为本发明实施例二提供的一种起伏地表地震数据偏移基准面转换装置的结构示意图。如图7所示，该起伏地表地震数据偏移基准面转换装置包括：第一获取模块21、第二获取模块22、插值模块23、计算模块24和校正转换模块25。
其中，第一获取模块21，用于获取起伏地表的偏移基准面的数据信息。
所述数据信息包括各共中心点CMP点的X坐标和Y坐标以及各CMP点的偏移基准面的第一校正量。
基于选取的起伏地表的偏移基准面，对偏移基准面进行数据查询，得到该偏移基准面的数据信息，其中，该数据信息包括CMP点的X坐和Y坐标，以及各CMP点的偏移基准面的第一校正量。其中，获取各CMP点的偏移基准面的第一校正量的过程包括：将各CMP点偏移基准面的高程与对应地震参考基准面高程求差值，使用该差值与替换速度做比值，该比值的2倍作为CMP点的第一校正量。可选地，所述数据信息中还可以包括各CMP点的线号以及点号，其中线号和点号是为了使用的需要对CMP点做的二维编号。
本实施例中偏移基准面的选取可参见上附图2以及实施例一中相关内容的记载，此处不再赘述。
第二获取模块22，用于获取炮点和检波点的坐标信息。
所述坐标信息包括炮点的X坐标和Y坐标，以及检波点的X坐标和Y坐标。
本实施例中，在进行地震勘查时，获取了一些炮点的大地坐标值，即炮点的X坐标和Y坐标。同样地，获取了一些检波点的大地坐标值，即检波点的X坐标和Y坐标。可选地，坐标信息中还可以包括各炮点以及检波点的标志号，例如顺序编写的序号。
与第一获取模块21和第二获取模块22连接的插值模块23，用于根据所述炮点和所述检波点以及所述CMP点的物理位置，将与所述炮点和所述检波点相邻的所述CMP点对应的第一校正量插值到所述炮点和所述检波点上，得到所述炮点的偏移基准面的第二校正量和所述检波点的偏移基准面的第二校正量。
其中所述物理位置是由X坐标和Y坐标确定的。
具体地，炮点的物理位置是由炮点的X坐标和Y坐标确定的，检波点的物理位置是由检波点的X坐标和Y坐标确定的，CMP点的物理位置是由CMP点的X坐标和Y坐标确定的。在得知炮点的物理位置以及CMP点的物理位置后，能够根据炮点的物理位置获取与该物理位置相邻的CMP点。在获取到与炮点相邻的CMP点后，根据与炮点相邻的CMP点对应的第一校正量得到炮点的偏移基准面的第二校正量。
类似地，在得知检波点的物理位置以及CMP点的物理位置后，能够根据检波点的物理位置获取与该物理位置相邻的CMP点。在获取到与该检波点相邻的CMP点后，根据与检波点相邻的CMP点对应的第一校正量得到检波点的偏移基准面的第二校正量。
与插值模块23连接的计算模块24，用于根据所述CMP点的所述第一校正量与所述炮点的所述第二校正量和所述检波点的所述第二校正量，计算得到CMP局部平面与偏移基准面之间的转换校正量。
具体地，计算模块24将炮点的偏移基准面的第二校正量和检波点的偏移基准面的第二校正量相加，然后将相加的和的一半与CMP点的偏移基准面的第一校正量做差，得到的差值为CMP局部平面与起伏地表的偏移基准面之间的转换校正量。
与计算模块24连接的校正转换模块25，用于使用所述转换校正量对待转换数据进行校正得到目标数据。
一方面，如果待转换数据为CMP道集数据，校正转换模块25用于将待转换CMP道集数据进行重排得到炮集数据，再使用获取到的转换校正量对重排得到的炮集数据进行校正，得到校正到偏移基准面上的目标数据。本实施例中，所述目标数据为校正到所述偏移基准面上的炮集数据。
另一方面，如果待转换数据为炮集数据，校正转换模块25用于将待转换的炮集数据进行重排得到CMP道集数据，再使用转换校正量对CMP道集数据进行反向校正，得到校正到CMP局部平面上的目标数据。本实施例中，所述目标数据为校正到CMP局部平面上的CMP道集数据。
本实施例提供的起伏地表地震数据偏移基准面转换装置，获取起伏地表的偏移基准面的数据信息，其中数据信息包括各共中心点CMP点的X坐标和Y坐标以及各CMP点的偏移基准面的第一校正量，获取炮点和检波点的坐标信息，坐标信息包括炮点的X坐标和Y坐标以及检波点的X坐标和Y坐标，根据炮点和检波点以及CMP点的物理位置，将与炮点和检波点相邻的CMP点对应的第一校正量插值到炮点和检波点上，得到炮点的偏移基准面的第二校正量和检波点的偏移基准面的第二校正量，根据CMP点的第一校正量与炮点的第二校正量和检波点的第二校正量，计算得到CMP局部平面与起伏地表的偏移基准面之间的转换校正量，使用转换校正量对待转换数据进行校正得到目标数据。本实施例基于CMP点的第一校正量获取炮点和检波点的第二校正量，从而得到获取CMP局部平面与偏移基准面之间的转换校正量，基于该转换校正量能在不同系统或偏移方法中进行地震数据的转换，使得转换数据过程简单准确，使用转换后的数据能获得正确的成像效果，成像深度保持一致。
实施例三
图8为本发明实施例三提供的另一种起伏地表地震数据转换装置的结构示意图。如图8所示，该起伏地表地震数据转换装置包括上述实施例二中的第一获取模块21、第二获取模块22、插值模块23、计算模块24和校正转换模块25。
本实施例中，插值模块23一种可选的结构方式包括：第一查询单元231、第一获取单元232、第二查询单元233和第二获取单元234。
其中，第一查询单元231，用于根据所述炮点的X坐标和Y坐标查询所述数据信息，获取与所述炮点相邻的所述CMP点。
与第一查询单元231连接的第一获取单元232，用于根据所述与所述炮点相邻的所述CMP点对应的所述第一校正量得到所述炮点的偏移基准面的第二校正量。
第一获取单元232将与炮点相邻的CMP点对应的第一校正量进行加权平均，得到炮点的第二校正量。一般与炮点相邻的CMP点的个数为4个，对这4个CMP点的第一校正量进行加权平均，就可以得到炮点的第二校正量。
第二查询单元233，用于根据所述检波点的X坐标和Y坐标查询所述数据信息，获取与所述检波点相邻的所述CMP点。
与第二查询单元233连接的第二获取单元234，用于根据所述与所述检波点相邻的所述CMP点对应的所述第一校正量得到所述检波点的偏移基准面的第二校正量。
第二获取单元234将与检波点相邻的CMP点对应的第一校正量进行加权平均，得到检波点的第二校正量。一般与检波点相邻的CMP点的个数为4个，对这4个CMP点的第一校正量进行加权平均，就可以得到检波点的第二校正量。
进一步地，校正转换模块25的一种可选的结构方式包括：重排单元251和校正单元252。
其中，重排单元251，用于在所述待转换数据为CMP道集数据时，将待转换CMP道集数据进行重排得到炮集数据。
与重排单元251连接的校正单元252，用于使用所述转换校正量对所述炮集数据进行校正，得到校正到偏移基准面上的所述目标数据，其中，所述目标数据为校正到所述偏移基准面上的炮集数据。
进一步地，重排单元251，还用于在所述待转换数据为炮集数据时，将待转换的所述炮集数据进行重排得到CMP道集数据。
校正单元252，还用于使用所述转换校正量对所述CMP道集数据进行反向校正，得到校正到CMP局部平面上的所述目标数据，其中，所述目标数据为校正到所述CMP局部平面上的CMP道集数据。
本实施例基于CMP点的第一校正量获取炮点和检波点的第二校正量，从而得到获取CMP局部平面与起伏地表的偏移基准面之间的转换校正量，基于该转换校正量能在不同系统或偏移方法中进行地震数据的转换，使得转换数据过程简单准确，使用转换后的数据能获得正确的成像效果，成像深度保持一致。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (10)
1.一种起伏地表地震数据偏移基准面转换方法，其特征在于，包括：
获取起伏地表的偏移基准面的数据信息，所述数据信息包括各CMP点的X坐标和Y坐标以及各CMP点的偏移基准面的第一校正量；
获取炮点和检波点的坐标信息，所述坐标信息包括炮点的X坐标和Y坐标，以及检波点的X坐标和Y坐标；
根据所述炮点和所述检波点以及所述CMP点的物理位置，将与所述炮点和所述检波点相邻的所述CMP点对应的第一校正量插值到所述炮点和所述检波点上，得到所述炮点的偏移基准面的第二校正量和所述检波点的偏移基准面的第二校正量；其中所述物理位置是由X坐标和Y坐标确定的；
根据所述CMP点的所述第一校正量与所述炮点的所述第二校正量和所述检波点的所述第二校正量，计算得到CMP局部平面与起伏地表的偏移基准面之间的转换校正量；
使用所述转换校正量对待转换数据进行校正得到目标数据。
2.根据权利要求1所述的起伏地表地震数据偏移基准面转换方法，其特征在于，所述根据所述炮点和所述检波点以及所述CMP点的物理位置，将与所述炮点和所述检波点相邻的所述CMP点对应的第一校正量插值到所述炮点和所述检波点上，得到所述炮点的偏移基准面的第二校正量和所述检波点的偏移基准面的第二校正量，包括：
根据所述炮点的X坐标和Y坐标查询所述数据信息，获取与所述炮点相邻的所述CMP点；
根据所述与所述炮点相邻的所述CMP点对应的所述第一校正量得到所述炮点的偏移基准面的第二校正量；
根据所述检波点的X坐标和Y坐标查询所述数据信息，获取与所述检波点相邻的所述CMP点；
根据所述与所述检波点相邻的所述CMP点对应的所述第一校正量得到所述检波点的偏移基准面的第二校正量。
3.根据权利要求2所述的起伏地表地震数据偏移基准面转换方法，其特征在于，所述根据所述与所述炮点相邻的所述CMP点对应的所述第一校正量得到所述炮点的偏移基准面的第二校正量，包括：
将所述与所述炮点相邻的所述CMP点对应的所述第一校正量进行加权平均，得到所述炮点的所述第二校正量；
所述根据所述与所述检波点相邻的所述CMP点对应的所述第一校正量得到所述检波点的偏移基准面的第二校正量，包括：
将所述与所述检波点相邻的所述CMP点对应的所述第一校正量进行加权平均，得到所述检波点的所述第二校正量。
4.根据权利要求1所述的起伏地表地震数据偏移基准面转换方法，其特征在于，如果所述待转换数据为CMP道集数据，所述使用所述转换校正量对待转换数据进行校正得到目标数据，包括：
将待转换CMP道集数据进行重排得到炮集数据；
使用所述转换校正量对所述炮集数据进行校正，得到校正到偏移基准面上的所述目标数据，其中，所述目标数据为校正到所述偏移基准面上的炮集数据。
5.根据权利要求1所述的起伏地表地震数据偏移基准面转换方法，其特征在于，如果所述待转换数据为炮集数据，所述使用所述转换校正量对待转换数据进行校正得到目标数据，包括：
将待转换的所述炮集数据进行重排得到CMP道集数据；
使用所述转换校正量对所述CMP道集数据进行反向校正，得到校正到CMP局部平面上的所述目标数据，其中，所述目标数据为校正到所述CMP局部平面上的CMP道集数据。
6.一种起伏地表地震数据偏移基准面转换装置，其特征在于，包括：
第一获取模块，用于获取起伏地表的偏移基准面的数据信息，所述数据信息包括各CMP点的X坐标和Y坐标以及各CMP点的偏移基准面的第一校正量；
第二获取模块，用于获取炮点和检波点的坐标信息，所述坐标信息包括炮点的X坐标和Y坐标，以及检波点的X坐标和Y坐标；
插值模块，用于根据所述炮点和所述检波点以及所述CMP点的物理位置，将与所述炮点和所述检波点相邻的所述CMP点对应的第一校正量插值到所述炮点和所述检波点上，得到所述炮点的偏移基准面的第二校正量和所述检波点的偏移基准面的第二校正量；其中所述物理位置是由X坐标和Y坐标确定的；
计算模块，用于根据所述CMP点的所述第一校正量与所述炮点的所述第二校正量和所述检波点的所述第二校正量，计算得到CMP局部平面与起伏地表偏移基准面之间的转换校正量；
校正转换模块，用于使用所述转换校正量对待转换数据进行校正得到目标数据。
7.根据权利要求6所述的起伏地表地震数据偏移基准面转换装置，其特征在于，所述插值模块，具体包括：
第一查询单元，用于根据所述炮点的X坐标和Y坐标查询所述数据信息，获取与所述炮点相邻的所述CMP点；
第一获取单元，用于根据所述与所述炮点相邻的所述CMP点对应的所述第一校正量得到所述炮点的偏移基准面的第二校正量；
第二查询单元，用于根据所述检波点的X坐标和Y坐标查询所述数据信息，获取与所述检波点相邻的所述CMP点；
第二获取单元，用于根据所述与所述检波点相邻的所述CMP点对应的所述第一校正量得到所述检波点的偏移基准面的第二校正量。
8.根据权利要求7所述的起伏地表地震数据偏移基准面转换装置，其特征在于，所述第一获取单元，具体用于将所述与所述炮点相邻的所述CMP点对应的所述第一校正量进行加权平均，得到所述炮点的所述第二校正量；
所述第二获取单元，具体用于将所述与所述检波点相邻的所述CMP点对应的所述第一校正量进行加权平均，得到所述检波点的所述第二校正量。
9.根据权利要求6所述的起伏地表地震数据偏移基准面转换装置，其特征在于，所述校正转换模块，具体包括：重排单元和校正单元；
所述重排单元，用于在所述待转换数据为CMP道集数据时，将待转换CMP道集数据进行重排得到炮集数据；
所述校正单元，用于使用所述转换校正量对所述炮集数据进行校正，得到校正到偏移基准面上的所述目标数据，其中，所述目标数据为校正到所述偏移基准面上的炮集数据。
10.根据权利要求9所述的起伏地表地震数据偏移基准面转换装置，其特征在于，所述重排单元，还用于在所述待转换数据为炮集数据时，将待转换的所述炮集数据进行重排得到CMP道集数据；
所述校正单元，还用于使用所述转换校正量对所述CMP道集数据进行反向校正，得到校正到CMP局部平面上的所述目标数据，其中，所述目标数据为校正到所述CMP局部平面上的CMP道集数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410720019.3A CN105717538B (zh) | 2014-12-02 | 2014-12-02 | 起伏地表地震数据偏移基准面转换方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410720019.3A CN105717538B (zh) | 2014-12-02 | 2014-12-02 | 起伏地表地震数据偏移基准面转换方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105717538A CN105717538A (zh) | 2016-06-29 |
| CN105717538B true CN105717538B (zh) | 2018-02-02 |
Family
ID=56146405
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410720019.3A Active CN105717538B (zh) | 2014-12-02 | 2014-12-02 | 起伏地表地震数据偏移基准面转换方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105717538B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109655888B (zh) * | 2017-10-11 | 2020-08-25 | 中国石油化工股份有限公司 | 地震数据处理中光滑浮动基准面的定量选择方法及系统 |
| CN109085644A (zh) * | 2018-07-30 | 2018-12-25 | 中国石油化工股份有限公司 | 基于双射线走时的真地表成像方法 |
| CN111999767B (zh) * | 2020-07-21 | 2023-09-26 | 中国石油天然气集团有限公司 | 起伏地表的偏移成像方法及装置 |
| CN113075732B (zh) * | 2021-03-25 | 2022-04-01 | 青岛海洋地质研究所 | 一种消除高分辨率小多道地震地层异常起伏的方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2691811A1 (fr) * | 1992-05-29 | 1993-12-03 | Total Sa | Procédé de détermination d'un modèle de distribution optimale des vitesses de propagation d'ondes sismiques dans le sous-sol. |
| US6810153B2 (en) * | 2002-03-20 | 2004-10-26 | Hitachi Software Global Technology, Ltd. | Method for orthocorrecting satellite-acquired image |
| GB2378507B (en) * | 2001-06-08 | 2005-08-17 | Pgs Americas Inc | Extraction of p-wave and s-wave velocities from multi-component seismic data by joint velocity inversion |
| CN1712992A (zh) * | 2004-06-25 | 2005-12-28 | 中国石油化工股份有限公司 | 一种在地震勘探中炮点或检波点的校正方法及其系统 |
| CN101900831A (zh) * | 2010-06-25 | 2010-12-01 | 恒泰艾普石油天然气技术服务股份有限公司 | 地震资料处理真地表条件下的椭圆展开成像方法和装置 |
| CN101984366A (zh) * | 2010-09-29 | 2011-03-09 | 北京吉星吉达科技有限公司 | 基于共孔径面的起伏地表叠前时间偏移方法及设备 |
| CN102879819A (zh) * | 2012-09-18 | 2013-01-16 | 中国石油天然气股份有限公司 | 保持地震波场运动学特征的地震数据处理方法及装置 |
| CN102914791A (zh) * | 2011-08-05 | 2013-02-06 | 中国石油天然气集团公司 | 一种起伏地表地震数据处理的克希霍夫叠前时间偏移方法 |
-
2014
- 2014-12-02 CN CN201410720019.3A patent/CN105717538B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2691811A1 (fr) * | 1992-05-29 | 1993-12-03 | Total Sa | Procédé de détermination d'un modèle de distribution optimale des vitesses de propagation d'ondes sismiques dans le sous-sol. |
| GB2378507B (en) * | 2001-06-08 | 2005-08-17 | Pgs Americas Inc | Extraction of p-wave and s-wave velocities from multi-component seismic data by joint velocity inversion |
| US6810153B2 (en) * | 2002-03-20 | 2004-10-26 | Hitachi Software Global Technology, Ltd. | Method for orthocorrecting satellite-acquired image |
| CN1712992A (zh) * | 2004-06-25 | 2005-12-28 | 中国石油化工股份有限公司 | 一种在地震勘探中炮点或检波点的校正方法及其系统 |
| CN101900831A (zh) * | 2010-06-25 | 2010-12-01 | 恒泰艾普石油天然气技术服务股份有限公司 | 地震资料处理真地表条件下的椭圆展开成像方法和装置 |
| CN101984366A (zh) * | 2010-09-29 | 2011-03-09 | 北京吉星吉达科技有限公司 | 基于共孔径面的起伏地表叠前时间偏移方法及设备 |
| CN102914791A (zh) * | 2011-08-05 | 2013-02-06 | 中国石油天然气集团公司 | 一种起伏地表地震数据处理的克希霍夫叠前时间偏移方法 |
| CN102879819A (zh) * | 2012-09-18 | 2013-01-16 | 中国石油天然气股份有限公司 | 保持地震波场运动学特征的地震数据处理方法及装置 |
Non-Patent Citations (3)
| Title |
|---|
| 变换坐标系下相移法起伏地表地震波场延拓;王祥春 等;《地球物理学进展》;20050930;第20卷(第3期);第677-680页 * |
| 基于起伏地表的混合法叠前深度偏移;王成祥;《石油地球物理勘探》;20020630;第37卷(第3期);第219-223页 * |
| 起伏地表二维声波方程地震波场延拓;王祥春 等;《中国矿业大学学报》;20090930;第38卷(第5期);第702-706页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105717538A (zh) | 2016-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104459782B (zh) | 使用薄层的沿层层析反演地震速度建模方法以及建模单元 | |
| CN105301639B (zh) | 基于vsp旅行时双加权层析反演速度场的方法及其装置 | |
| CN105717538B (zh) | 起伏地表地震数据偏移基准面转换方法及装置 | |
| CN113552625B (zh) | 一种用于常规陆域地震数据的多尺度全波形反演方法 | |
| CN104775811B (zh) | 一种地层各向异性信息的提取及校正方法和系统 | |
| CN105093319B (zh) | 基于三维地震数据的地面微地震静校正方法 | |
| CN102901985A (zh) | 一种适用于起伏地表的深度域层速度修正方法 | |
| CN104422963B (zh) | 一种变偏移距vsp资料时差校正方法 | |
| CN104360396B (zh) | 一种海上井间tti介质三种初至波走时层析成像方法 | |
| CN101984366A (zh) | 基于共孔径面的起伏地表叠前时间偏移方法及设备 | |
| CN102116869A (zh) | 高精度叠前域最小二乘偏移地震成像技术 | |
| CN104181593B (zh) | 一种三维无射线追踪回折波层析成像方法及装置 | |
| CN105301636A (zh) | 速度模型的建立方法和装置 | |
| CN105242315B (zh) | 一种地质构造成像的方法及装置 | |
| CN107894618A (zh) | 一种基于模型平滑算法的全波形反演梯度预处理方法 | |
| CN104865597A (zh) | 一种深度域层速度初始模型的建模方法 | |
| CN106405664A (zh) | 一种磁异常化极方法 | |
| CN107817516A (zh) | 基于初至波信息的近地表建模方法及系统 | |
| CN102338887B (zh) | 不规则尺寸空变网格层析成像静校正方法 | |
| CN104932021A (zh) | 基于反向射线追踪的约束层析速度建模方法 | |
| CN103852789A (zh) | 用于地震数据的非线性层析方法及其装置 | |
| CN107870361A (zh) | 一种地震回折波层析成像方法、装置及终端设备 | |
| CN104597496A (zh) | 一种二维地震资料中速度的三维空间归位方法 | |
| CN104950330B (zh) | 气顶油藏深度域成像的速度建模方法及系统 | |
| CN104570100B (zh) | 多子波克希霍夫地震数据偏移方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |