南方海相地层近地表结构数据处理方法

摘要

提供了一种南方海相地层近地表结构数据处理方法，包括：对工区内的近地表结构调查点采集近地表结构数据；从近地表结构数据样本库中寻找工区内的未调查点的至少一个特征所对应的近地表结构数据的数据区间，并将未调查点的近地表结构数据约束在所述数据区间内；判断所述数据区间的变化范围是否小于预定值，如果所述数据区间的变化范围不小于预定值，则增加未调查点的特征以对未调查点的近地表结构数据进行进一步约束，直到约束后的数据区间小于预定值；使用调查点的近地表结构数据和未调查点的近地表结构数据的数据区间来获得工区内的近地表结构数据分布图。

说明书

技术领域

本申请涉及地震勘探领域，更具体地，涉及一种南方海相底层近地表结构数据处理方法。

背景技术

近地表结构调查是地震勘探的一个重要环节，在指导地震勘探野外生产的地震激发接收参数设计和室内地震资料静校正处理中都起到了重要的作用。目前采用的近地表结构调查的数据处理方法主要通过在调查区内每平方公里布设0.25～1个近地表结构调查点，采集近地表结构数据，并在处理时采用纯数学数据插值方法，对没有采集近地表结构数据的区域自动生成相似插值数据，最终制作生成整个调查区内的近地表结构图，其中，所述近地表结构数据主要包括近地表各层的速度数据和厚度数据，同时也可包括其它数据，诸如潜水面的深度数据、地表岩性数据、地形条件数据等。上述方法在地表地形变化不大的北方和岩石性质变化不大的陆相地层区，计算出的近地表结构数据的精度相对较高。

然而，南方海相地层具有如下特点：①海相沉积地层主要出露老地层，多为碳酸盐岩(即，灰岩)和砂岩地层出露地表，近地表速度和厚度变化大；②南方地形主要为山地，地表起伏大，近地表的结构变化快。因此，采用上述常规的近地表结构调查的数据处理方法计算出的近地表结构数据不能真实反映近地表的结构变化情况，误差较大，从而不能完全满足地震勘探的需要。

发明内容

本发明的目的在于提供一种南方海相地层近地表结构数据处理方法，包括：对工区内的近地表结构调查点采集近地表结构数据；从近地表结构数据样本库中寻找工区内的未调查点的至少一个特征所对应的近地表结构数据的数据区间，并将未调查点的近地表结构数据约束在所述数据区间内；判断所述数据区间的变化范围是否小于预定值，如果所述数据区间的变化范围不小于预定值，则增加未调查点的特征以对未调查点的近地表结构数据进行进一步约束，直到约束后的数据区间小于预定值；以及使用调查点的近地表结构数据和未调查点的近地表结构数据的数据区间来获得工区内的近地表结构数据分布图。

近地表结构数据可包括近地表各层的速度数据、厚度数据和潜水面深度数据中的至少一个。

可对速度数据设置100m/s的预定值，对厚度数据设置1m的预定值，并对潜水面深度设置0.5m的预定值。

可预先通过对各种近地表结构进行近地表结构调查采集来获得近地表结构数据样本库，或者将工区内的近地表结构调查点的近地表结构数据作为近地表结构数据样本库。

未调查点的特征可包括岩石性质、地形条件和潜水面深度。

地形条件可包括地表高程、坡度和离陡崖的距离。

工区内的近地表结构数据分布图可包括近地表各层的速度数据分布图、厚度数据分布图和潜水面深度数据分布图中的至少一个。

未调查点可基于数据处理的需要或近地表结构的变化情况而被分布。

数据区间可以是至少一个定值数据的集合。

可通过对小于预定值的约束后的数据区间的两个端点值取中值来获得定值数据，并使用调查点的近地表结构数据和所述定值数据来获得工区内的近地表结构数据分布图。

将在接下来的描述中部分阐述本发明另外的方面和/或优点，还有一部分通过描述将是清楚的，或者可以经过本发明的实施而得知。

附图说明

通过下面结合附图进行的详细描述，本发明的上述和其它目的和特点将会变得更加清楚，其中：

图1是示出根据本发明示例性实施例的南方海相地层近地表结构数据处理方法的流程图。

具体实施方式

现在，详细描述本发明的示例性实施例，其示例在附图中表示，其中，相同的标号始终表示相同的部件。

图1是示出根据本发明示例性实施例的南方海相地层近地表结构数据处理方法的流程图。

参照图1，在步骤S110，对工区内的近地表结构调查点采集近地表结构数据。这里，所述近地表结构数据可包括近地表各层的速度数据、厚度数据和潜水面深度数据中的至少一个，但应理解本发明不限于此，所述近地表结构数据还可包括其它数据，诸如地表岩性数据、地形条件数据等；可在调查区内每平方公里布设0.25～1个近地表结构调查点，但本发明不限于此，本领域技术人员完全可按照实际情况和需要对调查点进行相应的调整。此外，可通过诸如小折射、微测井等现有方法来对近地表结构调查进行采集，这里将不再赘述。

在步骤S120，从近地表结构数据样本库中寻找工区内的未采集数据区中的未调查点的至少一个特征所对应的近地表结构数据的数据区间，并将未调查点的近地表结构数据约束在所述数据区间内。这里，所述未调查点的特征可包括但不限于岩石性质、地形条件和潜水面深度，所述地形条件可包括但不限于地表高程、坡度和离陡崖的距离，并且未调查点的这些特征可由本领域技术人员通过地质调查等方法而被确定；所述数据区间可以是多个定值数据的集合，并且本领域技术人员可基于数据处理的需要或近地表结构的变化情况来分布所述未调查点。此外，可按照岩石性质、地形条件和潜水面深度的顺序来对近地表结构数据进行约束，但这仅是示例，本领域技术人员完全可根据需要对约束的顺序进行调整和改变；仅作为示例，如果将某一未调查点的岩石性质作为特征，并且所述岩石性质是灰岩，则从近地表结构数据样本库中得到灰岩所对应的近地表结构数据的数据区间(诸如，速度数据区间、厚度数据区间和潜水面的深度数据区间等)，并将所述未调查点的近地表结构数据约束在获得的数据区间内。

在步骤S130，判断步骤S120中的所述数据区间的变化范围是否小于预定值。这里，本领域技术人员可根据实际需要来设置所述预定值，诸如，可对速度数据设置100m/s(米/秒)的预定值，对厚度数据设置1m(米)的预定值，并对潜水面深度设置0.5m(米)的预定值。

如果所述数据区间的变化范围小于预定值，则在步骤S140，使用调查点的近地表结构数据和未调查点的近地表结构数据的数据区间来获得工区内的近地表结构数据分布图。这里，工区内的近地表结构数据分布图可包括但不限于近地表各层的速度数据分布图、厚度数据分布图和潜水面的深度数据分布图中的至少一个。此外，可通过诸如数据成图等现有方法来获得近地表结构数据分布图，这里将不再赘述。

如果所述数据区间的变化范围不小于预定值，则在步骤S150增加未调查点的特征，并且操作返回步骤S120，即，从近地表结构数据样本库中寻找未调查点的所述特征和所述增加的特征所共同对应的近地表结构数据的数据区间，以对未调查点的近地表结构数据进行进一步约束，直到约束后的数据区间的变化范围小于预定值。

仅作为示例，如果上述某一未调查点的通过岩石性质约束后，获得了灰岩所对应的近地表结构数据的数据区间并且该数据区间的变化范围不小于预定值，则可对所述未调查点继续进行约束，诸如，如果所述未调查点的高程为800m，则可从近地表结构数据样本库中获得岩石性质是灰岩并且高程是800m所对应的近地表结构数据的数据区间(诸如，速度数据区间、厚度数据区间和潜水面的深度数据区间等)，如果所述数据区间的变化范围仍然不小于预定值，则继续进行约束，诸如，如果所述未调查点的潜水面深度为10m，则可从近地表结构数据样本库中获得岩石性质是灰岩、高程是800m并且潜水面深度为10m所对应的近地表结构数据的数据区间，等等，直到约束后的数据区间的变化范围小于预定值小于预定值。

应该了解，上面所述的未调查点的近地表结构数据的数据区间可以是至少一个定值数据的集合，并且仅作为示例的是：如果未调查点的近地表结构数据的数据区间不是定值数据的集合，则可通过诸如对小于预定值的约束后的数据区间的两个端点值取中值来获得定值数据，并使用调查点的近地表结构数据和所述定值数据来获得工区内的近地表结构数据分布图。

根据本发明的示例性实施例，通过采用带地质条件约束的插值方法代替传统的纯数学数值插值方法，同时引入了统计学的原理，从而使地质相关性大大加强，提高了南方海相地层近地表结构调查数据处理精度，解决了南方海相地层近地表结构调查数据处理数据不准的问题，也间接提高了地震勘探后期野外施工激发接收参数设计精度和地震资料室内处理静校正的精度，具有可操作性强，精度高的特点；此外，本发明不仅适用于南方海相地区，还可以在其它类似地区进行推广。

虽然已经参照特定示例性实施例示出和描述了本发明，但是本领域的技术人员将理解，在不脱离范围由权利要求及其等同物限定的本发明的精神和范围的情况下可作出形式和细节上的各种改变。