CN104330827B - 一种地表模型静校正的处理方法 - Google Patents

Abstract

本发明涉及一种模型静校正的处理方法，属于石油天然气地震勘探技术领域。本发明通过工区的地表高程、微测井和小折射数据，建立地表模型。根据该模型计算工区内已测定的每个炮点、检波点的低、降速层的厚度与速度，以及高速层顶界面的高程。由任意炮点、检波点的地表高程和高速层顶界面高程之差，结合经验公式可计算每个炮点和检波点的表层静校正量值。然后将得到的每个炮点、检波点的静校正量加载到单炮记录上。该方法能够控制好整体地表表层顶底界面的形态和厚度，合理消除地表表层结构变化引起的静校正问题，弥补地表表层静校正处理难题的技术空白，并且能极大的降低成本，提高工作效率，具有广阔的市场应用前景。

Description

一种地表模型静校正的处理方法

技术领域

本发明涉及一种地表模型静校正的处理方法，属于石油天然气地震勘探技术领域。

背景技术

静校正是地震资料处理与解释的基础步骤，关系着最终成果的可靠性甚至成败。如果建立的近地表模型无法有效控制静校正量值，地震资料处理和解释结果的精度都将受到影响。在沙漠、戈壁、黄土塬、山地等复杂地表区，利用地表高程、微测井、小折射等数据，建立地表模型，利用发明的经验公式计算静校正量，可解决表层静校正问题。该经验公式的特点，在微测井、小折射采集数据密集的地区，效果较好，且效率较高，在沙漠厚度较深且高速层较平缓的地区，其效果更加明显。

发明内容

本发明的目的是提供一种地表模型静校正的处理方法，以解决石油地震资料处理中低降速带变化所引起的静校正问题。

本发明的技术方案是：

1)根据已有炮点、检波点高程，计算出工区内所有炮点检波点的地表高程；

2)利用工区的地表高程、微测井和小折射数据，建立地表模型，根据该模型计算工区内已测定的每个炮点、检波点的低降速层速度与厚度以及高速层顶界面高程；

3)根据得到的每个炮点、检波点的地表高程和高速层顶界面高程之差利用经验公式计算每个炮点、检波点的静校正量；

4)将得到的每个炮点、检波点的静校正量加载到单炮记录上。

所述步骤3)中的经验公式为：

t＝h/(V0+C*ln h)

其中t为静校正量，h为地表表层厚度，V0＝350，C＝39。

所述步骤2)是采用克里格内插法实现的，所采用的克里格公式为：

其中z(x_i)为观测值，i＝1,2…n，它们分别位于区域内x_i的位置，x₀是一个未采样点，λ_i为权重,且

本发明的有益效果是：本发明通过工区的地表高程、微测井和小折射数据，建立地表模型，根据该模型计算工区内已测定的每个炮点、检波点的低、降速层的厚度与速度，以及高速层顶界面的高程。由任意炮点、检波点的地表高程和高速层顶界面高程之差，结合经验公式可计算每个炮点和检波点的表层静校正量值。然后将得到的每个炮点、检波点的静校正量加载到单炮记录上。该方法能够控制好整体地表表层顶底界面的形态和厚度，合理消除地表表层结构变化引起的静校正问题，弥补地表表层静校正处理难题的技术空白，并且能极大的降低成本，提高工作效率，具有广阔的市场应用前景。

附图说明

图1某测线地表表层厚度示意图；

图2本发明得到的炮点和检波点静校正量曲线示意图；

图3本发明实施例中施加静校正量前后的单炮记录示意图；

图4本发明实施例中施加静校正量前后的叠加剖面示意图。

具体实施方式

下面结合附图对本发明的具体实施方式做进一步的说明。

本发明所提供的表层静校正量的计算方法是利用地表高程、微测井和小折射等数据，建立地表模型，然后根据经验公式计算静校正量，从而解决表层静校正问题，该方法具体实施步骤如下：

1.根据待处理地震信号工区的已有炮点和检波点高程，得到工区内地表高程，通过每个小折射、微测井点的低速层厚度和该点的地表表层高程，计算出该点的高速层顶界面高程。地表高程＝野外解释的高速层顶界面的高程+低降速带厚度，野外解释的高速层顶界面的高程和低降速带厚度由微测井和小折射解释后得到。

2.利用克里格内插法获得工区内每个炮点、检波点的低降速层速度与速层厚度以及高速层顶界面高程，所采用的克里格公式如下：

其中z(x_i)为观测值，i＝1,2…n，它们分别位于区域内x_i的位置，x₀是一个未采样点，λ_i为权重,且

克里格方法通过引进以距离为自变量的变异差函数来计算权值，由于变异函数既可以反映变量的空间结构特性，又可以反映变量的随机分布特性，因此利用克里格方法进行空间数据插值能够取得理想的效果。

3.利用经验公式计算每个炮点、检波点的静校正量，该公式如下：

t＝h/(V0+C*ln h)

其中t为静校正量，h为地表表层厚度，V0＝350，C＝39。

4.将计算得到的每个炮点、检波点的静校正量加载到单炮记录上，从而达到解决由于表层变化引起的静校正问题。

下面以某工区为例说明本发明的效果。图1为某测线的地表表层厚度示意图。通过本发明的方法计算出每个炮点和检波点的静校正量如图2所示，将得到的校正量加载到单炮记录上，如图3所示，其中左边为施加校正量前的单炮记录，右边为施加校正量后的单炮记录，对比图3可以看出，施加静校正量后，单炮初至更加光滑，同时部分时间段上隐约可见双曲特征。施加静校正量前后的叠加剖面如图4所示，其中左边为施加校正量前，右边为施加静校正量后，对比图4左右叠加剖面可以看出：施加静校正量后，叠加剖面同相轴连续性明显增强，资料信噪比得到大幅度提高，地震剖面处理效果明显得到改善。

Claims (2)

1.一种地表模型静校正的处理方法，其特征在于，该方法包括以下步骤：

1)根据已有炮点、检波点高程，计算出工区内所有炮点检波点的地表高程；

2)利用工区的地表高程、微测井和小折射数据，建立地表模型，根据该模型计算工区内已测定的每个炮点、检波点的低降速层速度与厚度以及高速层顶界面高程；

3)根据得到的每个炮点、检波点的地表高程和高速层顶界面高程之差利用经验公式计算每个炮点、检波点的静校正量；

4)将得到的每个炮点、检波点的静校正量加载到单炮记录上；

步骤3)中的经验公式为：

t＝h/(V0+C*ln h)

其中t为静校正量，h为地表表层厚度，V0＝350，C＝39。

2.根据权利要求1所述的地表模型静校正的处理方法，其特征在于，所述步骤2)是采用克里格内插法实现的，所采用的克里格公式为：

$x_0={\mathrm{\Σ}}_{i=1}^n{\mathrm{\λ}}_{i}z\left(x_i\right)$

其中z(x_i)为观测值，i＝1,2…n，它们分别位于区域内x_i的位置，x₀是一个未采样点，λ_i为权重,且