用于计算双向层替换浮动基准面的海拔高程的方法
【技术领域】
[0001]本发明涉及复杂探区石油和天然气地震勘探技术领域，特别涉及一种用于计算双向层替换浮动基准面的海拔高程的方法。
【背景技术】
[0002]在地震资料处理中，进行速度分析时引进一个中间参考面，该参考面将CDP道集内所有炮点和检波点的基准面静校正量进行平均，得到静校正量的低频分量，该低频分量就是浮动基准面。浮动基准面是一个时间面，是速度谱的起始点。
[0003]在进行地震资料变速成图时，解释人员一般都是在固定基准面成图选择浮动基准面作为构造解释的基准面，原因是解释人员所用的速度谱的起点就是浮动基准面。但是，解释人员实现上述目的时需要获得浮动基准面对应的海拔高程。浮动基准面的海拔高程计算是复杂探区变速成图的关键环节。
[0004]在实际应用中，有的解释人员将地形的平滑面作为浮动基准面的海拔高程，有的解释人员则根据固定基准面的海拔高程和替换速度反向推算浮动基准面的海拔高程。当近地表地质情况简单即地形平坦时，地震波沿横向的传播速度稳定，上述两种方法产生的误差较小，可以忽略不计。但是，对于我国西部地区的一些探区，诸如山地、山前带、沙漠、黄土塬和戈壁滩等复杂探区，其近地表地形起伏剧烈，地震波沿近地表横向的速度变化很大，使得高速顶界面不稳定，如果按照上述两种计算方法进行计算，将会产生很大的误差，所产生的误差无法再忽略不计，特别是对于低幅度构造勘探来说，上述两种计算方法所产生的误差足以导致勘探的失败。因此，上述两种计算方法无法较好地适用于地形起伏剧烈的复杂探区O
[0005]另外，在现有技术中，用于计算浮动基准面的海拔高程的方法还存在着空间不闭合的问题。因为在计算过程中是按照每个CDP道集所涉及到的炮点和检波点进行计算的，不同测线在交点处是不闭合的，这也会导致误差的产生。
[0006]综上所述，现有技术中用于计算浮动基准面的海拔高程的方法存在如下不足:一是仅仅能适用于地形平坦的探区，无法适用于近地表地形起伏剧烈的复杂探区；二是空间不闭合。
【发明内容】
[0007]本发明的目的是针对现有技术的上述缺点，提供一种用于计算双向层替换浮动基准面的海拔高程的方法。
[0008]本发明提供的用于计算双向层替换浮动基准面的海拔高程的方法，该方法包括:
[0009]确定多个观测点，获得每一个观测点的大地坐标和实际高程；
[0010]根据每一个观测点的大地坐标，通过坐标变换得到该观测点的相对坐标；
[0011]根据每一个观测点的相对坐标，获得低速层的速度和高速层顶界面的高程；
[0012]由低速层的速度和高速层顶界面的高程计算得到从实际地表下剥到高速层顶界面的下剥时间；
[0013]用替换速度替代实际的低速层，保持下剥时间不变，由下剥时间、高速层顶界面的高程和替换速度得到每一个观测点的拟地形高程；
[0014]按照排列长度对每一个观测点的拟地形高程进行平滑处理，得到该观测点的相对于拟地形的浮动基准面的海拔高程；
[0015]按照排列长度对每一个观测点的实际高程进行平滑处理，得到该观测点的相对于实际地形的浮动基准面的海拔高程；
[0016]对工区内所有的二维测线的相对于实际地形的浮动基准面进行三维空间的平滑处理，以保证相交的测线闭合；
[0017]根据传播时间相等原则，由替换速度、高速层顶界面的高程和相对于拟地形的浮动基准面的海拔高程，计算得到实际地形的浮动基准面的海拔高程对应的低速层的速度。
[0018]优选地，观测点包括炮点和检波点。
[0019]优选地，排列长度为:对于每一个炮点，该炮点与最远的检波点之间的沿测线的距离。
[0020]优选地，传播时间相等原则为:对于任意一个观测点，从实际地形的高程以实际低速层的速度传播到高速层顶界面的时间与从拟地形高程以替换速度传播到高速层顶界面的时间相等;从相对于实际地形的浮动基准面以低速层的速度传播到高速层顶界面的时间与从相对于拟地形的浮动基准面以替换速度传播到高速层顶界面的时间相等。
[0021 ] 优选地，大地坐标为以参考椭球面为基准面的坐标。
[0022]上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂，以下特举本发明的【具体实施方式】。
[0023]本发明具有如下有益效果:
[0024](I)与现有技术的计算浮动基准面的海拔高程的方法相比，本发明的方法通过双向层替换能够明显降低计算误差，从而得到精确的浮动基准面的海拔高程；
[0025](2)本发明提供的方法能够实现相交测线的闭合，从而解决浮动基准面的海拔高程计算过程中的空间不闭合的问题；
[0026](3)本发明提供的方法不仅能够适用于地形平坦的探区，而且还能够适用于近地表地形起伏剧烈的复杂探区。
【附图说明】
[0027]通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中:
[0028]图1为本发明实施例提供的用于计算双向层替换浮动基准面的海拔高程的方法的流程示意图。
【具体实施方式】
[0029]下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例，然而应当理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本公开，并且能够将本公开的范围完整的传达给本领域的技术人员。
[0030]图1示出了本发明实施例提供的用于计算双向层替换浮动基准面的海拔高程的方法的流程示意图，如图1所示，该方法包括如下步骤:
[0031 ] 步骤S1:确定多个观测点，获得每