近地表模型构建方法与静校正量求取方法技术领域
本发明涉及地球物理勘探技术领域，具体来讲，涉及一种能够提高近地
表结构模型构建以及地震数据炮点与接收点静校正量计算精度和效率的方
法。
背景技术
目前，采用层析反演方法获得近地表速度、厚度信息并建立近地表速度
模型，进而求取炮点和接收点静校正量的方法已经成为当前地震勘探的主流
技术。
现有的层析反演技术方法通常是利用地震单炮记录上的初至波旅行时，
估计地下地层速度、厚度结构。在现有技术的常规地震勘探中，地震波激发
点间距80～100米，接收点间距一般在30～60米之间，所得的数据即为常规
的地震数据。由于常规地震数据空间采样间隔较大，地震波传播的射线密度
较为稀疏，对于近地表(例如，10米以内)的浅层地层速度和厚度的反演结
果误差较大，降低了静校正处理的精度。因此，现有技术中也有采用表层层
析勘探(SurfaceTomographicProspecting，简称STP)方法的初至波信息建立
表层模型，即采用小药量、小接收点距(例如，5～10米)和大激发点距(例
如，100～300米，为了减少钻井成本)，获得表层层析勘探数据，这样能够
通过加密接收点采样间隔，即增加地震波传播的射线密度，降低近地表浅层
结构模型的反演误差。
目前，在现有技术中，为了应用常规地震数据和表层层析勘探数据获得
高精度的炮点和接收点静校正量，需要采用以下技术步骤：
(1)分别拾取表层层析勘探数据初至波旅行时和常规地震数据初至波旅
行时；
(2)根据表层层析勘探数据初至波旅行时应用层析反演方法，从表层层
析勘探数据中获得高精度的近地表结构模型；
(3)根据常规地震数据资料的初至波旅行时，并将上述获得的高精度近
地表结构模型作为初始模型用于常规地震数据的层析反演，以获得最终的用
于炮点和接收点静校正量计算的近地表模型。
然而，现有技术中的上述方法主要存在以下三点不足：一是，需要进行
两次层析反演，因为层析反演属于非线性反演方法，故而导致计算量大，耗
时多；二是，虽然表层层析勘探获得的高精度模型可以作为初始模型，但是
由于常规地震数据空间采样间隔较大，地震波传播的射线密度并没有增加，
最终用常规地震数据通过层析反演获得的近地表模型的精度会有所降低；三
是，由于常规地震数据采集与表层层析勘探数据采集所用的观测系统不同，
两种观测系统的炮点和接收点不会完全重合，因此，仅仅用步骤(2)得到的
高精度模型不能获得常规地震数据处理所需的炮点和接收点的静校正量，需
要通过上述步骤(3)获得，或者将步骤(2)的结果加载到常规地震数据中，
处理过程较为繁琐。
为了解决现有技术的上述不足，本发明提出了一种简化而高效的近地表
结构模型层析反演与静校正方法。
发明内容
本发明的目的在于解决现有技术存在的上述不足中的至少一项。
为了实现上述目的，本发明的一方面提供了一种近地表模型构建方法，
所述方法包括以下步骤：分别拾取表层层析勘探数据的初至波旅行时和常规
地震数据的初至波旅行时；在同一层析反演方法中同时应用表层层析勘探数
据的初至波旅行时和常规地震数据的初至波旅行时，获得高精度的近地表结
构模型。
本发明的另一方面提供了一种近地表模型构建方法，所述方法包括以下
步骤：分别拾取第一初至波旅行时和第二初至波旅行时，所述第一初至波旅
行时为表层层析勘探数据的初至波旅行时，所述第二初至波旅行时为常规地
震数据的初至波旅行时；读入第一初至波旅行时、第二初至波旅行时、以及
给定的初始近地表速度-深度模型；采用层析反演技术，通过正演模拟获得第
三初至波旅行时，将第三初至波旅行时同时分别与所述第一初至波旅行时和
第二初至波旅行时进行对比，通过修正所述初始近地表速度-深度模型，以使
第三初至波旅行时与第一、第二初至波旅行时之间的总误差达到精度要求(例
如，误差平方和最小或尽量小)，然后，输出最终的近地表速度-深度模型。
本发明的另一方面提供了一种静校正量求取方法，所述方法包括以下步
骤：分别拾取表层层析勘探数据的初至波旅行时和常规地震数据的初至波旅
行时；在同一层析反演方法中同时应用表层层析勘探数据的初至波旅行时和
常规地震数据的初至波旅行时，获得高精度的近地表结构模型；根据所述高
精度的近地表结构模型进一步获得表层层析勘探数据的炮点的静校正量和接
收点的静校正量、以及处理所述常规地震数据时所需的炮点的静校正量和接
收点的静校正量。
本发明的另一方面还提供了一种静校正量求取方法，所述方法包括以下
步骤：分别拾取第一初至波旅行时和第二初至波旅行时，所述第一初至波旅
行时为表层层析勘探数据的初至波旅行时，所述第二初至波旅行时为常规地
震数据的初至波旅行时；读入第一初至波旅行时、第二初至波旅行时、以及
给定的初始近地表速度-深度模型；采用层析反演技术，通过正演模拟获得第
三初至波旅行时，将第三初至波旅行时同时分别与所述第一初至波旅行时和
第二初至波旅行时进行对比，通过修正所述初始近地表速度-深度模型，以使
第三初至波旅行时与第一、第二初至波旅行时之间的误差达到精度要求(例
如，误差平方和最小或尽量小)，然后，输出最终的近地表速度-深度模型；
利用所述最终的近地表速度-深度模型和给定的地震数据处理基准面，分别计
算出表层层析勘探数据的炮点的静校正量和接收点的静校正量、以及常规地
震数据的炮点的静校正量和接收点的静校正量。
与现有技术相比，本发明的有益效果包括以下方面：
1)简化了技术流程，将以往两次层析反演改进为一次层析反演，极大地
提高了处理效率；
2)在层析反演中同时应用表层层析勘探数据和常规地震数据初至波旅行
时，增加了地震波在近地表浅-中层地层中传播的射线密度，显著高于表层层
析勘探数据或常规地震数据的射线密度，进而提高了近地表结构模型的反演
精度，尤其是对地震数据静校正处理至关重要的10米以内的浅层；
3)同时应用表层层析勘探数据和常规地震数据初至波旅行时，能够同时
获得表层层析勘探和常规地震采集的炮点和接收点处的近地表模型数据和静
校正量值，也简化了地震数据静校正处理流程。
附图说明
图1示出了根据本发明的近地表模型构建方法与静校正量求取方法的一
个示意流程图。
图2示出了合并数据集的时-距关系图。
具体实施方式
在下文中，将结合示例性实施例来详细说明本发明的近地表模型构建方
法与静校正量求取方法。
在本发明的一个示例性实施例中，本发明的近地表模型构建方法是通过
在一次层析反演算法中同时应用表层层析勘探数据的初至波旅行时和常规地
震数据的初至波旅行时来构建近地表模型，其步骤可以为：(1)分别拾取表
层层析勘探数据初至波旅行时和常规地震数据初至波旅行时；(2)在层析反
演方法中同时应用表层层析勘探数据初至波旅行时和常规地震数据初至波旅
行时，从而获得高精度的近地表结构模型。
本发明的近地表模型构建方法也可通过以下步骤来实现：
1、分别拾取表层层析勘探数据和常规地震数据初至波旅行时(直接在两
种数据上读取)。
2、读入上述两种初至波旅行时数据和给定的初始近地表速度-深度模型。
具体而言，可通过解释由上述两种初至波数据合并后的数据集来获得初始近
地表速度-深度模型，即使所获得的初始近地表速度-深度模型的精度很低亦
可。例如，根据上述两种初至波数据合并后的数据集内个数据的时-距关系图
(横坐标为激发点和接收点之间的距离，纵坐标为处置时间)，判断出图中的
拐点位置及拐点个数n，并沿横坐标以拐点为端点划分出n+1层(相当于，
n+1条线段连成的折线)，每一层的速度为其线段斜率的倒数，从而以每一层
速度与时间的乘积构成每一层的初始近地表速度-深度模型，各层厚度之和即
为给定的初始近地表速度-深度模型。
3、采用层析反演技术，将通过正演模拟获得的理论初至波旅行时同时与
读取的表层层析勘探数据初至波旅行时和常规地震数据初至波旅行时进行对
比，通过修正初始近地表速度-深度模型，使正演获得的理论初至波旅行时与
表层层析勘探数据初至波旅行时和常规地震数据初至波旅行时之间的误差达
到精度要求(例如，误差平方和(即所有误差的平方和)最小或尽可能小，
或者小于某一给定值)时，输出最终的近地表速度-深度模型。这里，对初始
近地表速度-深度模型的修正属于现有层析反演技术的常规做法，通常可进行
多次修正。理论初至波旅行时与表层层析勘探数据初至波旅行时和常规地震
数据初至波旅行时之间的误差平方和越小，则表示反演得到的近地表速度-深
度模型与实际模型越接近。
在本发明的另一个示例性实施例中，本发明的静校正量求取方法是通过
在一次层析反演算法中同时应用表层层析勘探数据的初至波旅行时和常规地
震数据的初至波旅行时来构建近地表模型，随后获得地震数据的炮点和接收
点静校正量，其步骤可以为：(1)分别拾取表层层析勘探数据初至波旅行时
和常规地震数据初至波旅行时；(2)在层析反演方法中同时应用表层层析勘
探数据初至波旅行时和常规地震数据初至波旅行时，从而获得高精度的近地
表结构模型，进一步获得地震数据处理时所需的炮点的静校正量和接收点的
静校正量。
本发明的静校正量求取方法也可通过以下步骤来实现：
1、分别拾取表层层析勘探数据和常规地震数据初至波旅行时。
2、读入上述两种初至波旅行时数据和给定的初始近地表速度-深度模型。
3、采用层析反演技术，将通过正演模拟获得的理论初至波旅行时同时与
读取的表层层析勘探数据初至波旅行时和常规地震数据初至波旅行时进行对
比，通过修正初始近地表速度-深度模型，使正演获得的理论初至波旅行时与
表层层析勘探数据初至波旅行时和常规地震数据初至波旅行时之间的误差达
到精度要求(例如，误差平方和(即所有误差的平方和)最小或尽可能小)
时，输出最终的近地表速度-深度模型。
4、利用输出的近地表速度-深度模型和给定的地震数据处理基准面及替
换速度，分别计算出表层层析勘探数据的炮点和接收点静校正量、以及地震
勘探数据的炮点和接收点静校正量。例如，给定处理基准面和替换速度属于
本领域普通技术人员的常规技能。通过输出的近地表速度-深度模型，计算静
校正量的过程是属于本领域普通技术人员的基本技能。
下文将结合附图对本发明的一个实施例进行详细描述。
图1示出了本发明的一个示例性实施例的示意流程图。图2示出了合并
数据集的时-距关系图，其中，横坐标为激发点和接收点之间的距离，纵坐标
为初至时间。
第一步，以本领域公知的方法进行表层层析勘探采集，获得表层层析勘
探数据，并以相同的方法获取该数据的初至数据集，该数据集特征是由表层
层析勘探激发点到接收点之间的距离和地震波传播时间(简称为时-距关系)
来描述。
第二步，以本领域公知的方法进行野外地震勘探采集，获得地震勘探数
据，并以本领域公知的方法获取该数据的初至数据集，该数据集特征由地震
勘探激发点到接收点之间的距离和地震波传播时间(简称为时-距关系)来描
述。表层层析勘探与常规地震勘探采用不同的观测方式，通常只有部分接收
点和激发点重合。与常规地震勘探技术相比，表层层析勘探通常采用较稀疏
的炮点激发和较密集的检波点接收，表层层析勘探的初至波含有更多近地表
浅层结构的信息，常规地震勘探初至波含有更多深层结构信息。由于地震勘
探和表层层析勘探均采用多次覆盖技术(即，采用多点激发多道接收技术)，
因此，每一个初至数据集的数据量非常大，一般可达十万-百万道级。
第三步，分别读入的两种初至数据集并合并为一个数据集，数据集内各
数据的时-距关系保持不变，合并数据集更全面地反映了近地表浅层-深层地层
结构，使处理人员更容易识别出拐点位置(图2中的A点)。处理人员根据拐
点解释出折线AB和AC，AB的斜率的倒数表示近地表第一层的速度V1，AC
的斜率的倒数表示近地表第二层的速度V2，第一层的厚度H1可由B点时间T1与
第一层的速度V1相乘获得。至此，建立了初始的近地表速度-深度模型，该模
型为水平层状均匀速度模型，不能真实地反应变化剧烈的近地表速度结构。
图2中有一个拐点，初始模型为两层结构，如有2个拐点，则可以解释成3
层结构模型，3个拐点，可获得4层结构模型，……，以此类推。
第四步，利用合并后的初至数据集和初始速度-深度模型，即可进行层析
反演。层析反演能够精细估计地质模型，是地震勘探领域常用的技术。反演
结果具有多解性是层析反演算法的固有问题。在本发明中通过输入合并的初
至数据集和初始模型来增加层析反演的输入数据(或信息)，从而有利于降低
反演结果的多解性。采用层析反演技术，将通过正演模拟获得的理论初至波
旅行时同时与读取的表层层析勘探数据初至波旅行时和常规地震数据初至波
旅行时进行对比，通过修正初始近地表速度-深度模型，使正演获得的理论初
至波旅行时与表层层析勘探数据初至波旅行时和常规地震数据初至波旅行时
之间的误差达到精度要求(例如，误差平方和(即所有误差的平方和)最小
或尽可能小)时，输出最终的近地表速度-深度模型。
第五步，层析反演结果是获得了近地表速度-深度模型，该模型能够精细
描述近地表速度的纵横向变化规律，划分近地表较浅处的低降速层和较深处
的高速层，更加符合实际的地质情况。处理人员输入模型中低降速层的速度
值(比如，800米/秒)，即可获得低降速层底界面的深度。表层层析勘探和常
规地震勘探的所有激发点和接收点在该模型的相应位置处均有速度和深度
值。
第六步，求取地震勘探炮点、检波点静校正量是地震数据处理的必需环
节，也是本领域的常规技术。处理人员根据近地表模型，输入低降速层的速
度、低降速层底界面的深度以及给定的处理基准面、替换速度，即可获得地
震勘探和表层层析勘探所有炮点与接收点的静校正量。
综上所述，本发明不仅简化了应用表层层析勘探数据和常规地震数据建
立近地表结构模型以及求取炮点和接收点静校正量的技术流程，而且提高了
近地表结构模型的精度，进而提高了地震数据静校正处理的效率。本发明的
适用范围包括但不限于表层层析勘探数据和常规地震数据的层析反演与静校
正应用，还可适用于多种不同观测系统的地震数据的联合应用，建立近地表
结构模型，求取所有参与反演的地震数据的炮点和接收点静校正量。本发明
可广泛用于需要对近地表和深层地层结构进行精细了解的地震勘探、矿山开
采、建筑、水利等工程物探领域。
尽管上面已经结合附图和示例性实施例描述了本发明，但是本领域普通
技术人员应该清楚，在不脱离权利要求的精神和范围的情况下，可以对上述
实施例进行各种修改。