一种基于长基距的采集参数确定方法和装置

摘要

本发明公开了一种基于长基距的采集参数确定方法，所述方法包括：获得待处理工区预置频率中每个频率所对应的长基距，其中，所述长基距是指基距长度大于25m的基距；针对该长基距设置所述每个频率所对应的组合单元的组内距，其中，所示组合单元为由n个检波器/震源组成的检波器/震源组合单元，n为大于1的正整数；根据所述每个预设频率所对应的组合单元的组内距，确定检波器/震源的空间分布；根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布；根据缩减的检波器/震源进行采集组合单元的分布进行地震数据采集。通过本发明的方案，在基于长基距地震数据采集时，可以缩减地震数据采集的工作量。

权利要求

1.一种基于长基距的采集参数确定方法，其特征在于，所述方法包括：

获得待处理工区预置频率中每个频率所对应的长基距，其中，所述长基距是指长度大于25m的基距；

针对该长基距设置所述每个频率所对应的组合单元的组内距，其中，所述组合单元为由n个检波器/震源组成的检波器/震源组合单元，n为大于1的正整数；

根据每个预设频率所对应的组合单元的组内距，确定检波器/震源的空间分布；

根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布，包括：获待处理工区预置频率的检波器/震源的空间分布后，将所有的频率所对应的检波器/震源组合单元的组内距进行排序；

选取检波器/震源组合单元的组内距中最小值和最大值所对应的组内距，按照缩减工作量的比例和缩减公式计算得到缩减的检波器/震源组合单元的分布。

2.根据权利要求1所述基于长基距的采集参数确定方法，其特征在于，所述获得待处理工区预置频率中每个频率所对应的长基距，包括：根据待处理工区的相关地质信息，确定待处理工区检波器的入射角范围；

根据所述检波器的入射角范围，采用预先设置的相控阵组合模型得到预置频率中每个频率所对应的基距。

3.根据权利要求2所述基于长基距的采集参数确定方法，其特征在于，所述相控阵组合模型，包括：其中， 为相控阵组合单元的响应；n为检波器/震源组合单元的数目；m为组合单元内检波器/震源的顺序号；组合单元为由n个检波器/震源组成的检波器/震源组合单元；Dm为检波器组合单元的能量强度；rm为每个检波器/震源组合单元到组合单元中心的距离；θ0为检波器的入射信号/出射信号的方位角； 为检波器的入射角/出射角定义；θm为检波器/震源的方位角； 为入射角/出射角；Vwater为水的速度，为组合响应 在‑π/2≤ψ0≤π/2范围内绝对值的和,Aη为在‑η≤ψ0≤η的入射角范围内的接收/激发的能量响应，Aη和 的比值作为组合响应能量AP，η为检波器的入射角范围。

4.根据权利要求3所述基于长基距的采集参数确定方法，其特征在于，所述采用预先设置的相控阵组合模型得到预置频率中每个频率所对应的基距，包括：根据所述待处理工区的相控阵组合模型，得到预置频率中每个频率相对应的能量曲线，其中，所述预置频率是根据所述待处理工区的地震资料的频率范围所设置的低频段频率；

根据每个频率的能量曲线确定该频率对应的长基距。

5.根据权利要求4所述基于长基距的采集参数确定方法，其特征在于，所述针对该长基距设置所述每个频率所对应的组合单元的组内距，包括：针对每个频率所对应的长基距，采用预先设置的检波器/震源的组合单元数量，分别设计每个频率所对应的检波器/震源组合单元的组内距；

结合每个频率所对应的检波器/震源组合单元的组内距，得到该待处理工区预置频率的检波器/震源的空间分布。

6.根据权利要求1所述基于长基距的采集参数确定方法，其特征在于，

其中，所述缩减公式，包括：

其中，rji为检波器/震源组合单元的原组内距，j为频率的编号，i为组合单元的编号，min(rji)为原组内距中的最小组内距，max(rji)为原组内距中的最大组内距，N为原检波器/震源组合单元的数目，1/M为缩减工作量的比例，M为2的整数倍，M＝2，4，8，16.......，k为缩减的检波器/震源组合单元的顺序号，k＝1,2,3,...........。

7.根据权利要求6所述基于长基距的采集参数确定方法，其特征在于，所述根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布，还包括：分别针对每个频率的检波器/震源的组合单元进行重排列，按照缩减工作量的比例和缩减公式计算得到对应的该频率的缩减的检波器/震源组合单元的分布；

获得所有频率缩减的检波器/震源组合单元的分布后，对所有频率的缩减后检波器/震源组合单元进行重新排列；

根据重新排列的缩减后的检波器/震源组合单元得到缩减的检波器/震源组合单元的分布。

8.一种基于长基距的采集参数确定装置，所述装置包括：存储器和处理器；其特征在于：所述存储器，用于保存用于长基距的采集参数确定的程序；

所述处理器，用于读取执行所述用于长基距的采集参数确定的程序，执行如下操作：获得待处理工区预置频率中每个频率所对应的长基距，其中，所述长基距是指长度大于25m的基距；

针对该长基距设置所述每个频率所对应的组合单元的组内距，其中，所述组合单元为由n个检波器/震源组成的检波器/震源组合单元，n为大于1的正整数；

根据每个预设频率所对应的组合单元的组内距，确定检波器/震源的空间分布；

根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布，包括：获待处理工区预置频率的检波器/震源的空间分布后，将所有的频率所对应的检波器/震源组合单元的组内距进行排序；

选取检波器/震源组合单元的组内距中最小值和最大值所对应的组内距，按照缩减工作量的比例和缩减公式计算得到缩减的检波器/震源组合单元的分布。

9.根据权利要求8所述基于长基距的采集参数确定装置，其特征在于，所述获得待处理工区预置频率中每个频率所对应的长基距，包括：根据待处理工区的相关地质信息，确定待处理工区检波器的入射角范围；

根据所述检波器的入射角范围，采用预先设置的相控阵组合模型得到预置频率中每个频率所对应的基距。

10.根据权利要求9所述基于长基距的采集参数确定装置，其特征在于，所述相控阵组合模型，包括：其中， 为相控阵组合单元的响应；n为检波器/震源组合单元的数目；m为组合单元内检波器/震源的顺序号；组合单元为由n个检波器/震源组成的检波器/震源组合单元；Dm为检波器组合单元的能量强度；rm为每个检波器/震源组合单元到组合单元中心的距离；θ0为检波器的入射信号/出射信号的方位角； 为检波器的入射角/出射角定义；θm为检波器/震源的方位角； 为入射角/出射角；Vwater为水的速度，为组合响应 在‑π/2≤ψ0≤π/2范围内绝对值的和,Aη为在‑η≤ψ0≤η的入射角范围内的接收/激发的能量响应，Aη和 的比值作为组合响应能量AP，η为检波器的入射角范围。

11.根据权利要求10所述基于长基距的采集参数确定装置，其特征在于，所述采用预先设置的相控阵组合模型得到预置频率中每个频率所对应的基距，包括：根据所述待处理工区的相控阵组合模型，得到预置频率中每个频率相对应的能量曲线，其中，所述预置频率是根据所述待处理工区的地震资料的频率范围所设置的低频段频率；

根据每个频率的能量曲线确定该频率对应的长基距。

12.根据权利要求11所述基于长基距的采集参数确定装置，其特征在于，所述针对该长基距设置所述每个频率所对应的组合单元的组内距，包括：针对每个频率所对应的长基距，采用预先设置的检波器/震源的组合单元数量，分别设计每个频率所对应的检波器/震源组合单元的组内距；

结合每个频率所对应的检波器/震源组合单元的组内距，得到该待处理工区预置频率的检波器/震源的空间分布。

13.根据权利要求8所述基于长基距的采集参数确定装置，其特征在于，其中，所述缩减公式，包括：

其中，rji为检波器/震源组合单元的原组内距，j为频率的编号，i为组合单元的编号，min(rji)为原组内距中的最小组内距，max(rji)为原组内距中的最大组内距，N为原检波器/震源组合单元的数目，1/M为缩减工作量的比例，M为2的整数倍，M＝2，4，8，16.......，k为缩减的检波器/震源组合单元的顺序号，k＝1,2,3,...........。

14.根据权利要求13所述基于长基距的采集参数确定装置，其特征在于，所述根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布，还包括：分别针对每个频率的检波器/震源的组合单元进行重排列，按照缩减工作量的比例和缩减公式计算得到对应的该频率的缩减的检波器/震源组合单元的分布；

获得所有频率缩减的检波器/震源组合单元的分布后，对所有频率的缩减检波器/震源组合单元进行重新排列；

根据重新排列的缩减的检波器/震源组合单元得到缩减的检波器/震源组合单元的分布。

说明书

一种基于长基距的采集参数确定方法和装置

技术领域

[0001] 本发明涉及地球物理勘探领域，尤指一种基于长基距的采集参数确定方法和装置。

背景技术

[0002] 海上空气枪震源常采用气枪组合方式来达到增强有效信号能量，压制气泡及噪音。传统地震勘探中也利用震源组合的方向性来增强地震波下传能量。但受限于震源组合
子阵基距，常规设计的基距为20m左右，对高频信号的定向增强明显，对10Hz以下的低频信号的定向性增强能力还是较差。

[0003] 高速层屏蔽区的地震勘探是地震勘探中的难点。高速层屏蔽区的地震勘探受作业环境的影响和采集方式的限制，使得高速层屏蔽区的勘探更加艰难。但经过不断地探索、研究发现，解决上述高速层屏蔽区的勘探所存在难题的有效手段是高覆盖、低频震源和宽方
位。传统的勘探中通过增大震源的总容量以及增加枪阵组合中的大枪比例来增强低频震源
能量，在实际的应用中往往不能达到期望的效果，在高速屏蔽区更是收效甚微。陆地的低速和高速屏蔽探区，受地表复杂的地质和环境背景噪音的影响，不利于震源的定向性设计；海域探区因海水层的存在，提供了一层相对简单的表层结构，海水速度的时间和空间变化相
对较小，有利于震源的定向性设计，因此实施针对低频段地震数据进行能量增强的地震数
据采集方法，但是基于低频段的长基距地震数据采集的工作量特别大，对于地震勘探增加
了实施的难度，因此，如何能够提供一种在基于长基距地震数据采集时，缩减地震数据采集的工作量的采集参数确定方法是亟待解决的问题。

发明内容

[0004] 为了解决上述技术问题，本发明提供了一种基于长基距的采集参数确定方法，在基于长基距地震数据采集时，可以缩减地震数据采集的工作量。

[0005] 为了达到本发明目的，本发明提供了一种基于长基距的采集参数确定方法，所述方法包括：

[0006] 获得待处理工区预置频率中每个频率所对应的长基距，其中，所述长基距是指长度大于25m的基距；

[0007] 针对该长基距设置所述每个频率所对应的组合单元的组内距，其中，所述组合单元为由n个检波器/震源组成的检波器/震源组合单元，n为大于1的正整数；

[0008] 根据所述每个预设频率所对应的组合单元的组内距，确定检波器/震源的空间分布；

[0009] 根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布。

[0010] 一种示例性的实施例中，所述获得待处理工区预置频率中每个频率所对应的长基距，包括：

[0011] 根据待处理工区的相关地质信息，确定待处理工区检波器的入射角范围；

[0012] 根据所述检波器的入射角范围，采用预先设置的相控阵组合模型得到预置频率中每个频率所对应的基距。

[0013] 一种示例性的实施例中，其特征在于，所述相控阵组合模型，包括：

[0014]

[0015] 其中， 为相控阵组合单元的响应；n为检波器/震源组合单元的数目；m为组合单元内检波器/震源的顺序号；组合单元为由n个检波器/震源组成的检波器/震源组合单元；
Dm为检波器组合单元的能量强度；rm为每个检波器/震源组合单元到组合单元中心的距离；
θ0为检波器的入射信号/出射信号的方位角；ψ0为检波器的入射角/出射角定义；θm为检波器/震源的方位角；ψm为入射角/出射角；Vwater为水的速度， 为组合响应 在‑π/2≤ψ0≤π/
2范围内绝对值的和,Aη为在‑η≤ψ0≤η的入射角范围内的接收/激发的能量响应，Aη和 的比值作为组合响应能量AP，η为检波器的入射角范围。

[0016] 一种示例性的实施例中，所述采用预先设置的相控阵组合模型得到预置频率中每个频率所对应的基距，包括：

[0017] 根据所述待处理工区的相控阵组合模型，得到预置频率中每个频率相对应的能量曲线，其中，所述预置频率是根据所述待处理工区的地震资料的频率范围所设置的低频段
频率；

[0018] 根据每个频率的能量曲线确定该频率对应的长基距。

[0019] 一种示例性的实施例中，所述针对该长基距设置所述每个频率所对应的组合单元的组内距，包括：

[0020] 针对每个频率所对应的长基距，采用预先设置的检波器/震源的组合单元数量，分别设计每个频率所对应的检波器/震源组合单元的组内距；

[0021] 结合每个频率所对应的检波器/震源组合单元的组内距，得到该待处理工区预置频率的检波器/震源的空间分布。

[0022] 一种示例性的实施例中，所述根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布，包括：

[0023] 获待处理工区预置频率的检波器/震源的空间分布后，将所有的频率所对应的检波器/震源组合单元的组内距进行排序；

[0024] 选取检波器/震源组合单元的组内距中最小值和最大值所对应的组内距，按照缩减工作量的比例和缩减公式计算得到缩减的检波器/震源组合单元的分布。

[0025] 一种示例性的实施例中，其中，所述缩减公式，包括：

[0026]

[0027] 其中，rji为检波器/震源组合单元的原组内距，j为频率的编号，i为组合单元的编号，min(rji)为原组内距中的最小组内距，max(rji)为原组内距中的最大组内距，N为原检波器/震源组合单元的数目，1/M为缩减工作量的比例，M为2的整数倍，M＝2，4，8，16.......，k为缩减的检波器/震源组合单元的顺序号，k＝1,2,3,...........。

[0028] 一种示例性的实施例中，所述根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布，还包括：

[0029] 分别针对每个频率的检波器/震源的组合单元进行重排列，按照缩减工作量的比例和缩减公式计算得到对应的该频率的缩减的检波器/震源组合单元的分布；

[0030] 获得所有频率缩减的检波器/震源组合单元的分布后，对所有频率的缩减检波器/震源组合单元进行重新排列；

[0031] 根据重新排列的缩减的检波器/震源组合单元得到缩减的检波器/震源组合单元的分布。

[0032] 为了解决上述问题，本发明还提供了一种长基距的采集参数确定装置，所述装置包括：存储器和处理器；

[0033] 所述存储器，用于保存用于长基距的采集参数确定的程序；

[0034] 所述处理器，用于读取执行所述用于长基距的采集参数确定的程序，执行如下操作：

[0035] 获得待处理工区预置频率中每个频率所对应的长基距，其中，所述长基距是指长度大于25m的基距；

[0036] 针对该长基距设置所述每个频率所对应的组合单元的组内距，其中，所述组合单元为由n个检波器/震源组成的检波器/震源组合单元，n为大于1的正整数；

[0037] 根据所述每个预设频率所对应的组合单元的组内距，确定检波器/震源的空间分布；

[0038] 根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布。

[0039] 一种示例性的实施例中，所述获得待处理工区预置频率中每个频率所对应的长基距，包括：

[0040] 根据待处理工区的相关地质信息，确定待处理工区检波器的入射角范围；

[0041] 根据所述检波器的入射角范围，采用预先设置的相控阵组合模型得到预置频率中每个频率所对应的基距。

[0042] 一种示例性的实施例中，所述相控阵组合模型，包括：

[0043]

[0044] 其中， 为相控阵组合单元的响应；n为检波器/震源组合单元的数目；m为组合单元内检波器/震源的顺序号；组合单元为由n个检波器/震源组成的检波器/震源组合单元；
Dm为检波器组合单元的能量强度；rm为每个检波器/震源组合单元到组合单元中心的距离；
θ0为检波器的入射信号/出射信号的方位角；ψ0为检波器的入射角/出射角定义；θm为检波器/震源的方位角；ψm为入射角/出射角；Vwater为水的速度， 为组合响应 在‑π/2≤ψ0≤π/
2范围内绝对值的和,Aη为在‑η≤ψ0≤η的入射角范围内的接收/激发的能量响应，Aη和 的比值作为组合响应能量AP，η为检波器的入射角范围。

[0045] 一种示例性的实施例中，所述采用预先设置的相控阵组合模型得到预置频率中每个频率所对应的基距，包括：

[0046] 根据所述待处理工区的相控阵组合模型，得到预置频率中每个频率相对应的能量曲线，其中，所述预置频率是根据所述待处理工区的地震资料的频率范围所设置的低频段
频率；

[0047] 根据每个频率的能量曲线确定该频率对应的长基距。

[0048] 一种示例性的实施例中，所述针对该长基距设置所述每个频率所对应的组合单元的组内距，包括：

[0049] 针对每个频率所对应的长基距，采用预先设置的检波器/震源的组合单元数量，分别设计每个频率所对应的检波器/震源组合单元的组内距；

[0050] 结合每个频率所对应的检波器/震源组合单元的组内距，得到该待处理工区预置频率的检波器/震源的空间分布。

[0051] 一种示例性的实施例中，所述根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布，包括：

[0052] 获待处理工区预置频率的检波器/震源的空间分布后，将所有的频率所对应的检波器/震源组合单元的组内距进行排序；

[0053] 选取检波器/震源组合单元的组内距中最小值和最大值所对应的组内距，按照缩减工作量的比例和缩减公式计算得到缩减的检波器/震源组合单元的分布。

[0054] 一种示例性的实施例中，其中，所述缩减公式，包括：

[0055]

[0056] 其中，rji为检波器/震源组合单元的原组内距，j为频率的编号，i为组合单元的编号，min(rji)为原组内距中的最小组内距，max(rji)为原组内距中的最大组内距，N为原检波器/震源组合单元的数目，1/M为缩减工作量的比例，M为2的整数倍，M＝2，4，8，16.......，k为缩减的检波器/震源组合单元的顺序号，k＝1,2,3,...........。

[0057] 一种示例性的实施例中，所述根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布，还包括：

[0058] 分别针对每个频率的检波器/震源的组合单元进行重排列，按照缩减工作量的比例和缩减公式计算得到对应的该频率的缩减的检波器/震源组合单元的分布；

[0059] 获得所有频率缩减的检波器/震源组合单元的分布后，对所有频率的缩减检波器/震源组合单元进行重新排列；

[0060] 根据重新排列的缩减的检波器/震源组合单元得到缩减的检波器/震源组合单元的分布。

[0061] 与现有技术相比，本发明公开了一种基于长基距的采集参数确定方法，所述方法包括：获得待处理工区预置频率中每个频率所对应的长基距，其中，所述长基距是指基距长度大于25m的基距；针对该长基距设置所述每个频率所对应的组合单元的组内距，其中，所示组合单元为由n个检波器/震源组成的检波器/震源组合单元，n为大于1的正整数；根据所述每个预设频率所对应的组合单元的组内距，确定检波器/震源的空间分布；根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单
元的分布。通过本发明的方案，在基于长基距地震数据采集时，可以缩减地震数据采集的工作量。

[0062] 本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。

附图说明

[0063] 附图用来提供对本发明技术方案的进一步理解，并且构成说明书的一部分，与本申请的实施例一起用于解释本发明的技术方案，并不构成对本发明技术方案的限制。

[0064] 图1为本发明实施例的基于长基距的采集参数确定方法流程图；

[0065] 图2是本发明实施例的基于长基距的采集参数确定装置示意图；

[0066] 图3是本发明实施例一示例的高速屏蔽层顶界与底界反射/透射示意图；

[0067] 图4是本发明实施例一示例的高速屏蔽层地质模型示意图；

[0068] 图5是本发明实施例一示例的震源及检波器组合基距设计示意图；

[0069] 图6是本发明实施例一示例的高速屏蔽区多层上覆地层震源同入射角的关系模型；

[0070] 图7是本发明实施例一示例的震源及检波器组合下传能量聚焦示意图；

[0071] 图8是本发明实施例一示例的能量聚焦宽频带示意图；

[0072] 图9是本发明实施例一示例的能量聚焦高维度示意图；

[0073] 图10是本发明实施例一示例的不同频率所对应的基距的能量聚焦曲线；

[0074] 图11是本发明实施例一示例的不同频率的组合单元示意图；

[0075] 图12是本发明实施例一示例的长基线相控阵组合单元采集方案示意图；

[0076] 图13是本发明实施例一示例的缩减工作量后组合单元的示意图。

具体实施方式

[0077] 为使本发明的目的、技术方案和优点更加清楚明白，下文中将结合附图对本发明的实施例进行详细说明。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互任意组合。

[0078] 在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行。并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。

[0079] 图1是本发明实施例的基于长基距的采集参数确定方法流程图。

[0080] 步骤101.获得待处理工区预置频率中每个频率所对应的长基距。

[0081] 在本实施例中，预置频率是根据待处理工区的地震资料的频率范围所设置的低频段频率。该低频段的频率根据待处理工区地震数据进行分析，现有技术中低频段频率一般
指低于10HZ以下的频率。根据待处理工区的需求，可以选择10HZ以下频率作为预置频率，例如：2HZ、4HZ、6HZ、8HZ等。也可以根据需求选择10HZ以上频率作为预置频率。对于该预置频率的选择根据需求进行设置，对此并不进行具体的限定。

[0082] 在本实施例中，长基距是指基距长度大于25m的基距。在常规地震采集中，震源组合常用的基距在18m左右，而检波器组合常用的基距在12.5m左右。而本实施例中，根据低频段所对应的基距，例如：748m，372m，187m，93m等等，该基距的长度要大于常规采集所采用的基距，因此，可以将基于该低频信号所确定的基距定义为长基距即指基距长度大于25m的基距。

[0083] 一种示例性的实施例中，所述获得待处理工区预置频率中每个频率所对应的长基距，包括：根据待处理工区的相关地质信息，确定待处理工区检波器的入射角范围；根据所述检波器的入射角范围，采用预先设置的相控阵组合模型得到预置频率中每个频率所对应
的基距。建立待处理工区的地质模型；根据所述地质模型确定待处理工区中检波器的入射
角范围。在本实施例中，创建地质模型是以Snell定律和Zoeppritz方程式作为理论基础。
Snell定律和Zoeppritz方程式是地震勘探中常用的两个重要工具。根据Snell定律，在存在阻抗差异界面上，入射角和透射角同界面上下两侧的介质速度呈正比例关系，地震波从低
速介质向高速介质传播时，在阻抗差异界面上产生反射和透射，当入射角大于临界角时发
生全反射，没有地震波透射到高速介质内部。临界角的大小同阻抗界面上下地层的速度比
存在严格的比例关系，临界角的变化范围在0～90°，随着速度差异的增大，临界角变小。以高速屏蔽区勘探为例，在高速层两侧存在两个速度差异很大的界面。在该高速层中，上界面阻碍地震波能量向下传播，下界面阻拦地震波能量的向上传播。上界面和下界面的传递的
信号被限制在各自的临界角范围之内。如图3所示，根据Snell关系在高速层上界面存在如
下关系，

[0084]

[0085] 假设高速层速度V1＝5000m/s，密度为ρ1＝2.5g/cm3，上覆地层速度V0＝2000m/s，3
密度为ρ0＝1.8g/cm。则当θ0＝21.80°时，θ1＝90.00，地震波发生全反射，因此，高速层上界面即高速层的顶界面透射的地震波入射角需限制在θ0＝21.80°之内。

[0086] 如图3所示，根据Snell关系在高速层底界面存在如下关系，

[0087]

[0088] 高速层速度V1＝5000m/s，密度ρ1＝2.5g/cm3，下部地层速度V2＝3000m/s，密度ρ2＝3
2.1g/cm。高速层下界面之下地层反射地震波向上传播时，当θ2＝30.96°时，θ1＝90.00°，地震波发生全反射，因此，高速层的下界面即底界面向上透射的地震波限制在θ2＝30.96°之内。

[0089] 通过上述分析，高速层顶界面的下传地震波有效入射角范围θ0＝21.80°之内，底界面上传地震波的有效入射角范围为θ2＝30.96°之内。也就是说，在本例的高速屏蔽区中，震源激发的地震波顶界面下传入射角有效范围在θ0＝21.80°之内，此时，对应的高速层底界面上传的地震波被控制在θ2＝30.96°之内，当高速层下伏地层倾角在0～30.96°之内时，下伏地层地震波能够反射回去，下伏地层倾角在大于30.96°时，地震波发生全反射，反射信息返回地下，此时通过反射波不能对下伏地层成像。高速层顶界面阻抗差异限定了入射角
的有效范围，底界面阻抗差异使得下伏地层的构造倾角限定在下界面的临界角内。

[0090] 根据上述所确定高速层顶底界面的传播规律建立高速屏蔽地质模型以用于震源和采集方案的设计，所建立的高速屏蔽地质模型如图4，A点为震源位置，B点为检波点位置，X为检波点与震源之间的偏移距；h0为高速层顶界深度，h1为高速层的厚度，h2为高速层下目的层的厚度；V0为高速层上覆地层的速度，V1为高速层的速度，V2为高速层下伏地层的速度；
ρ0为高速层上覆地层的速度，ρ1为高速层的速度，ρ2为高速层下伏地层的速度；x0为射线在高速层上覆地层传播时在地表的投影距离，x1为射线在高速层内传播时在地表的投影距
离，x2为射线在下伏地层内传播时在地表的投影距离；θ0为高速层顶界入射角，θ1为高速层底界的入射角，θ2为高速层下伏地层的入射角。

[0091] 根据所建立的高速屏蔽地质模型，根据Snell关系在高速层上界面存在如下关系，

[0092]

[0093] 如果待处理工区中高速层速度V1＝5000m/s，密度为ρ1＝2.5g/cm3，上覆地层速度3
V0＝2000m/s，密度为ρ0＝1.8g/cm 。则当θ0＝21.80°时，θ1＝90.00°，地震波发生全反射，因此得到上界面即高速层的顶界面透射的地震波入射角需限制在θ0＝21.80°之内。由此可以确定满足上覆地层即上界面入射的入射角有效范围是θ0要小于临界角。

[0094] 一种示例性的实施例中，所述相控阵组合模型，包括：

[0095]

[0096] 其中， 为相控阵组合单元的响应；n为检波器/震源组合单元的数目；m为组合单元内检波器/震源的顺序号；组合单元为由n个检波器/震源组成的检波器/震源组合单元；
Dm为检波器组合单元的能量强度；rm为每个检波器/震源组合单元到组合单元中心的距离；
θ0为检波器的入射信号/出射信号的方位角；ψ0为检波器的入射角/出射角定义；θm为检波器/震源的方位角；ψm为入射角/出射角；Vwater为水的速度，为组合响应 在‑π/2≤ψ0≤π/2范围内绝对值的和,Aη为在‑η≤ψ0≤η的入射角范围内的接收/激发的能量响应，Aη和 的比值作为组合响应能量AP，η为检波器的入射角范围。

[0097] 现有技术中对于震源组合或检波器组合方案中，一般采用以下的相控阵组合的基距设计方法；

[0098]

[0099] 公式3是现有的相控阵组合公式，m为组合内检波器/震源顺序号，Dm为检波器的灵敏度或者震源的能量大小，v*为视速度，v为入射地震波速度(相对检波器)或者出射地震波速度(相对震源)，α为入射角(相对检波器)或者出射角(相对震源)，xm为检波器/震源到组合对称中心的距离，结合相控阵技术对公式(1)进行了改造，将组合响应由二维拓展到三
维，公式中的xm替换为rm；将沿入射角α入射信号拓展到在三维空间由入射信号/出射信号的方位角θ0和仰角ψ0(入射角＝π/2‑ψ0或者出射角＝π/2‑ψ0)以及检波器/震源的方位角θm和仰角ψm决定；其次，将视速度v*用速度用v替代，速度v定义为水速Vwater，得到灵敏度的三维组合响应公式4。

[0100]

[0101] 结合在高速屏蔽区的震源激发的能量集中在有效入射角范围之内，将有效入射角范围以及下传能量作为限定约束条件，获得重新构建的相控阵组合模型如公式5所示：

[0102]

[0103] 一种示例性的实施例中，所述采用预先设置的相控阵组合模型得到预置频率中每个频率所对应的基距，包括：根据所述待处理工区的相控阵组合模型，得到预置频率中每个频率相对应的能量曲线，其中，所述预置频率是根据所述待处理工区的地震资料的频率范
围所设置的低频段频率；根据每个频率的能量曲线确定该频率对应的长基距。在本实施例
中，

[0104] 第一步：根据所述待处理工区的相控阵组合模型，得到基距和组合响应能量的关系曲线，其中，所述基距为各个组合单元到组合中心的距离。基于相控阵组合模型，结合待处理工区的地震资料情况，分析该地区地震资料的低频段范围，并选择相应的低频进行确
定该频率的能量曲线，以用于确定基距。例如：在火成岩，灰岩，岩盐等存在的高速层屏蔽区，在实际的资料分析中发现，高速层上下地震信号的频带发生突变，尤其是下伏地层的频带分布在3～24Hz，分析根本原因是因为高速屏蔽层对高频信号形成了遮挡。针对于这种对于高速屏蔽区，根据相控阵组合模型确定低频段内信号的能量曲线，在3～24Hz频段内可以选取四个低频段内的频率f1，f2，f3和f4进行基距确定。例如：低频频率选取的频率f1＝
3Hz，f2＝6Hz，f3＝12Hz和f4＝24Hz；对于该低频段频率的选择可以根据相应的待处理工区的地震资料的品质和待处理工区的需求选择相应的频率。在选取四个低频段内的频率f1，
f2，f3和f4后，根据待处理工区的相控阵组合模型求取四个低频频率f1，f2，f3和f4的聚焦能量随基距变化聚焦能量曲线如图10所示。

[0105] 第二步：根据每个频率的能量曲线，确定每个频率所对应的基距。具体的实现过程包括：根据计算得到预置频率中每个频率相对应的能量曲线，选择每个频率的能量曲线中多个波峰的位置所对应的基距作为待确定的基距；对多个待确定的基距进行对比，选择能
量聚焦性最优的基距作为该频率对应的基距。如图7所示，图7‑a为24Hz频率的地震信号对应着不同的基距的能量聚焦示意图，图7‑a中示意的是能量聚焦随着基距的变化，根据公式
5中，以第一个式子中的rm做为变量(rm为每个检波器/震源组合单元到组合单元中心的距
离)；其余参数做为常量求取K，然后按照公式5求取AP，得到组合响应随基距变化的曲线图
7‑a。图7‑a中示意的是能量聚焦随着基距的变化，限定的入射角范围η＝21.10°，基距从0开始增大，到达93m的时候，能量占比达到最大达到54％，随着基距的增大，能量聚焦性整体变弱，局部有振荡，存在局部能量聚焦高值点。选取第一和第二高值点的基距的组合方案进行比较，获得组合响应 随入射角的变化，如图7‑b，实线为第一高值点的组合响应随入射角的变化，虚线为第二高值点组合响应随入射角的变化，从图7‑b可以看出第一高值点对应的方案能量聚焦性最好。因此，选择第一高值点所对应的基距作为该频率的基距。

[0106] 如图10所示，针对于上述所建立的工区的相控阵组合模型确定低频信号的能量曲线，根据上述的确定基距的实现步骤，求得四个频率对应的基距分别为Lf1＝748m，Lf2＝
372m，Lf3＝187m，Lf4＝93m。在常规地震采集中，震源组合常用的基距在18m左右，而检波器组合常用的基距在12.5m左右。而本实施例中，根据低频能量聚焦设计的基距要大于常规采集所采用的基距，因此，可以将基于该低频信号所确定的基距定义为长基距。

[0107] 步骤102.针对该长基距设置所述每个频率所对应的组合单元的组内距。

[0108] 在本实施例中，组合单元为由n个检波器/震源组成的检波器/震源组合单元，n为大于1的正整数。

[0109] 一种示例性的实施例中，针对该长基距设置每个频率所对应的长基距的组合单元的组内距，包括：针对每个频率所对应的长基距，采用预先设置的检波器/震源的组合单元数量，分别设计每个频率所对应的检波器/震源组合单元的组内距；结合每个频率所对应的检波器/震源组合单元的组内距，得到该待处理工区预置频率的检波器/震源的空间分布。
在本实施例中，得到该待处理工区预置频率的检波器/震源的空间分布可以是优化的检波
器/震源的组合方案也可以是等间距的检波器/震源的组合方案。如果优化的检波器/震源
的组合方案可以通过包括：获取待处理工区中炮点的间距和检波点的间距；采用规则化方
法将所获取每个频率所对应的组合单元的组内距进行规则化；使用规则化后的组内距进行
地震数据采集。

[0110] 所述采用规则化方法将所获取每个频率所对应的组合单元的组内距进行规则化，包括：将每个频率所对应的组合单元的组内距进行微分计算，得到微分结果；对所述微分结果取整计算获得微分取整结果；对所述微分取整结果进行积分计算得到炮点的间距和检波
点的间距整数倍的组内距。其中，规则化的公式如下所示：

[0111]

[0112] 其中，rij为任一预置频率对应的组合单元的组内距，SI为炮点间距，RI为炮点间距。

[0113] 在本实施例中，当确定每个频率所对应的检波器/震源组合单元的组内距后，可以确定该待处理工区预置频率的检波器/震源的空间分布。

[0114] 步骤103.根据每个预设频率所对应的组合单元的组内距，确定检波器/震源的空间分布。

[0115] 在本实施例中，根据每个预设频率所对应的组合单元的组内距后，可以确定地震数据采集观测系统设计里面的炮线距、炮点距、检波线距以及检波点距。根据所确定的炮线距、炮点距、检波线距以及检波点距可确定检波器/震源的空间分布。

[0116] 通过上述的所确定的炮线距、炮点距、检波线距以及检波点距可确定检波器/震源的空间分布，采集观测系统如图12所示，上述方案设计的不等间距组合方案，如图12a中，采用表5中针对f3＝12Hz频率确定的检波器/震源的空间分布，布设了8条测线，测线间距呈不等间距，类似用上述参数定义的观测系统，也可以拓展到其它频率增强所采用的不等间距
的方案，如图12b，图中实线为检波线，点线为炮线，不仅限于此，也可以采用多个频率段相互累加的方案，图12c方案所示。

[0117] 步骤104.根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布。

[0118] 在本实施例中，考虑到采用长基距所确定的检波器/震源的空间分布，该采集的工作量较大，需要通过采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布。在得到缩减的检波器/震源组合单元的分布后，可以根据该检波器/震源组合单元
的分布进行地震数据采集。在基于表4理想方案，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布，如表13。

[0119] 表13基于表4理想方案工作量简化表

[0120]

[0121]

[0122] 一种示例性的实施例中，根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布，包括：

[0123] 获待处理工区预置频率的检波器/震源的空间分布后，将所有的频率所对应的检波器/震源组合单元的组内距进行排序；选取检波器/震源组合单元的组内距中最小值和最
大值所对应的组内距，按照缩减工作量的比例和缩减公式计算得到缩减的检波器/震源组
合单元的分布。

[0124] 一种示例性的实施例中，其中，所述缩减公式，包括：

[0125]

[0126] 其中，rji为检波器/震源组合单元的原组内距，j为频率的编号，i为组合单元的编号，min(rji)为原组内距中的最小组内距，max(rji)为原组内距中的最大组内距，N为原检波器/震源组合单元的数目，1/M为缩减工作量的比例，M为2的整数倍，M＝2，4，8，16.......，k为缩减的检波器/震源组合单元的顺序号，k＝1,2,3,...........。

[0127] 一种示例性的实施例中，所述根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布，还包括：分别针对每个频率的检波器/震源的组合单元进行重排列，按照缩减工作量的比例和缩减公式计算得到对应
的该频率的缩减的检波器/震源组合单元的分布；获得所有频率缩减的检波器/震源组合单
元的分布后，对所有频率的缩减检波器/震源组合单元进行重新排列；根据重新排列的缩减的检波器/震源组合单元得到缩减的检波器/震源组合单元的分布。

[0128] 针对于每个频率分别进行工作量缩减，而后在对组合方案进行重排，如以表3中的四个频率f1，f2，f3和f4所对应的工作量缩减为例，如表14为f1频率组合方案缩减方案，表
15为f2频率组合方案缩减方案，表16为f3频率组合方案缩减方案，表17为f4频率组合方案
缩减方案。

[0129] 表14 f1频率组合方案缩减方案

[0130]

[0131] 表15 f2频率组合方案缩减方案

[0132]

[0133]

[0134] 表16 f3频率组合方案缩减方案

[0135]

[0136] 表17 f4频率组合方案缩减方案

[0137]

[0138] 表14到表17组合方案中，则针对某一频率的检波器组合方案算上原方案有三种(原方案+工作量缩减1/2方案+工作量缩减1/4方案)，此时对多频率段进行增强组合设计，
如同时采用f1+f2+f3+f4组合方案，则根据局排列组合公式，可生成的低频拓展组合方案共
81种方案。实际中可根据频率增强的需求，在81种方案中选择满足要求的方案。如表18为根据表14到表17四个频率同时缩减组合工作量1/4的方案示意。

[0139] 表18 f4频率组合方案缩减方案

[0140]

[0141]

[0142] 图13为缩减工作量前后组合单元的空间位置及组合数目的变化，此方案中组合单元数由32降为8。

[0143] 为了解决上述问题，如图2所示，本发明还提供了一种长基距的采集参数确定装置，所述装置包括：存储器和处理器；

[0144] 所述存储器，用于保存用于长基距的采集参数确定的程序；

[0145] 所述处理器，用于读取执行所述用于长基距的采集参数确定的程序，执行如下操作：

[0146] 获得待处理工区预置频率中每个频率所对应的长基距，其中，所述长基距是指基距长度大于25m的基距；

[0147] 针对该长基距设置所述每个频率所对应的组合单元的组内距，其中，所示组合单元为由n个检波器/震源组成的检波器/震源组合单元，n为大于1的正整数；

[0148] 根据所述每个预设频率所对应的组合单元的组内距，确定检波器/震源的空间分布；

[0149] 根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布。

[0150] 一种示例性的实施例中，所述获得待处理工区预置频率中每个频率所对应的长基距，包括：

[0151] 根据待处理工区的相关地质信息，确定待处理工区检波器的入射角范围；

[0152] 根据所述检波器的入射角范围，采用预先设置的相控阵组合模型得到预置频率中每个频率所对应的基距。

[0153] 一种示例性的实施例中，所述相控阵组合模型，包括：

[0154]

[0155] 其中， 为相控阵组合单元的响应；n为检波器/震源组合单元的数目；m为组合单元内检波器/震源的顺序号；组合单元为由n个检波器/震源组成的检波器/震源组合单元；
Dm为检波器组合单元的能量强度；rm为每个检波器/震源组合单元到组合单元中心的距离；
θ0为检波器的入射信号/出射信号的方位角；ψ0为检波器的入射角/出射角定义；θm为检波器/震源的方位角；ψm为入射角/出射角；Vwater为水的速度，为组合响应 在‑π/2≤ψ0≤π/2范围内绝对值的和,Aη为在‑η≤ψ0≤η的入射角范围内的接收/激发的能量响应，Aη和 的比值作为组合响应能量AP，η为检波器的入射角范围。

[0156] 一种示例性的实施例中，所述采用预先设置的相控阵组合模型得到预置频率中每个频率所对应的基距，包括：根据所述待处理工区的相控阵组合模型，得到预置频率中每个频率相对应的能量曲线，其中，所述预置频率是根据所述待处理工区的地震资料的频率范
围所设置的低频段频率；根据每个频率的能量曲线确定该频率对应的长基距。

[0157] 一种示例性的实施例中，所述针对该长基距设置所述每个频率所对应的组合单元的组内距，包括：针对每个频率所对应的长基距，采用预先设置的检波器/震源的组合单元数量，分别设计每个频率所对应的检波器/震源组合单元的组内距；结合每个频率所对应的检波器/震源组合单元的组内距，得到该待处理工区预置频率的检波器/震源的空间分布。

[0158] 一种示例性的实施例中，所述根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布，包括：获待处理工区预置频率的检波器/震源的空间分布后，将所有的频率所对应的检波器/震源组合单元的组内距进
行排序；选取检波器/震源组合单元的组内距中最小值和最大值所对应的组内距，按照缩减工作量的比例和缩减公式计算得到缩减的检波器/震源组合单元的分布。

[0159] 一种示例性的实施例中，其中，所述缩减公式，包括：

[0160]

[0161] 其中，rji为检波器/震源组合单元的原组内距，j为频率的编号，i为组合单元的编号，min(rji)为原组内距中的最小组内距，max(rji)为原组内距中的最大组内距，N为原检波器/震源组合单元的数目，1/M为缩减工作量的比例，M为2的整数倍，M＝2，4，8，16.......，k为缩减的检波器/震源组合单元的顺序号，k＝1,2,3,...........。

[0162] 一种示例性的实施例中，所述根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布，还包括：分别针对每个频率的检波器/震源的组合单元进行重排列，按照缩减工作量的比例和缩减公式计算得到对应
的该频率的缩减的检波器/震源组合单元的分布；获得所有频率缩减的检波器/震源组合单
元的分布后，对所有频率的缩减检波器/震源组合单元进行重新排列；根据重新排列的缩减的检波器/震源组合单元得到缩减的检波器/震源组合单元的分布。

[0163] 一种示例性的实施例基于长基距的采集参数确定方法，具体实现过程如下所示：

[0164] 步骤1.建立待处理工区的地质模型。

[0165] 在本实施例中，创建地质模型依据Snell定律和Zoeppritz方程式作为理论基础。Snell定律和Zoeppritz方程式是地震勘探中常用的两个重要工具。根据Snell定律，在存在阻抗差异界面上，入射角和透射角同界面上下两侧的介质速度呈正比例关系，地震波从低
速介质向高速介质传播时，在阻抗差异界面上产生反射和透射，当入射角大于临界角时发
生全反射，没有地震波透射到高速介质内部。临界角的大小同阻抗界面上下地层的速度比
存在严格的比例关系，临界角的变化范围在0～90°，随着速度差异的增大，临界角变小。以高速屏蔽区勘探为例，在高速层两侧存在两个速度差异很大的界面。在该高速层中，上界面阻碍地震波能量向下传播，下界面阻拦地震波能量的向上传播。上界面和下界面的传递的
信号被限制在各自的临界角范围之内。如图3，根据Snell关系在高速层上界面存在如下关
系，

[0166]

[0167] 假设高速层速度V1＝5000m/s，密度为ρ1＝2.5g/cm3，上覆地层速度V0＝2000m/s，3
密度为ρ0＝1.8g/cm。则当θ0＝21.80°时，θ1＝90.00°，地震波发生全反射，因此得到上界面即高速层的顶界面透射的地震波入射角需限制在θ0＝21.80°之内。

[0168] 如图3，根据Snell关系在高速层底界面存在如下关系，

[0169]

[0170] 高速层速度V1＝5000m/s，密度ρ1＝2.5g/cm3，下部地层速度V2＝3000m/s，密度ρ2＝3
2.1g/cm。高速层下界面之下地层反射地震波向上传播时，则当θ2＝30.96°时，θ1＝90.00°，地震波发生全反射，因此，高速层的下界面即底界面向上透射的地震波限制在θ2＝30.96°之内。

[0171] 通过上述分析，高速层顶界面的下传地震波有效入射角范围θ0＝21.80°之内，底界面上传地震波的有效入射角范围为θ2＝30.96°之内。也就是说，在本例的高速屏蔽区中，震源激发的地震波顶界面下传入射角有效范围在θ0＝21.80°之内，此时对应的高速层底界面上传的地震波被控制在θ2＝30.96°之内，当高速层下伏地层倾角在0～30.96°之内时，下伏地层地震波能够反射回去，下伏地层倾角在大于30.96°时，地震波发生全反射，反射信息返回地下，此时通过反射波不能对下伏地层成像。高速层顶界面阻抗差异限定了入射角的
有效范围，底界面阻抗差异使得下伏地层的构造倾角限定在下界面的。

[0172] 根据上述所确定高速层顶底界面的传播规律建立高速屏蔽地质模型以用于震源和采集方案的设计，如图4，A点为震源位置，B点为检波点位置，X为检波点与震源之间的偏移距；h0为高速层顶界深度，h1为高速层的厚度，h2为高速层下目的层的厚度；V0为高速层上覆地层的速度，V1为高速层的速度，V2为高速层下伏地层的速度；ρ0为高速层上覆地层的速度，ρ1为高速层的速度，ρ2为高速层下伏地层的速度；x0为射线在高速层上覆地层传播时在地表的投影距离，x1为射线在高速层内传播时在地表的投影距离，x2为射线在下伏地层内传播时在地表的投影距离；θ0为高速层顶界入射角，θ1为高速层底界的入射角，θ2为高速层下伏地层的入射角。

[0173] 步骤2.根据所述地质模型确定待处理工区检波器的入射角范围。

[0174] 在本实施例中，根据步骤1中所建立的高速屏蔽地质模型，根据Snell关系在高速层上界面存在如下关系，

[0175]

[0176] 如果待处理工区中高速层速度V1＝5000m/s，密度为ρ1＝2.5g/cm3，上覆地层速度3
V0＝2000m/s，密度为ρ0＝1.8g/cm 。则当θ0＝21.80°时，θ1＝90.00°，地震波发生全反射，因此得到上界面即高速层的顶界面透射的地震波入射角需限制在θ0＝21.80°之内。由此可以确定确定满足上覆地层即上界面入射的入射角有效范围是θ0要小于临界角。

[0177] 步骤3.根据上述建立的高速屏蔽地质模型以及地震波在地下传播的几何关系，建立高速屏蔽区的有效偏移距评价模型。

[0178] 在本步骤中，根据上述建立的高速屏蔽地质模型以及地震波在地下传播的几何关系，可以得到检波点与震源之间的偏移距X的公式，

[0179] X＝2(h0 tanθ0+h1 tanθ1+h2 tanθ2)                     (6)

[0180] 并且结合三角关系，得到公式7，cscθ0，cscθ1，cscθ2为各入射角对应的余割函数，[0181]

[0182] 另外，则根据Snell定律存在如下的关系，

[0183]

[0184] 联合公式6和公式7进一步进行简化，将公式中θ1和θ2消掉，则得到公式9，

[0185]

[0186] 根据上述的所确定满足上覆地层即上界面入射的入射角有效范围是θ0要小于临界角。定义顶界面的临界角为θ0c，底界面的临界角为θ2c。则公式9进一步简化为公式10，[0187]

[0188] 由于在实际勘探中速度大小随着深度的埋深加大而增大，一般情况下，高速屏蔽层上覆地层速度要低于下伏地层的速度，即存在关系V0

[0189]

[0190] 定义11式右侧为检波点与震源之间的偏移距X的最小偏移距Xmin设计准则，建立高速屏蔽区的有效偏移距评价模型如公式12

[0191]

[0192] 步骤4.确定待处理工区中的震源出射角及检波器入射角范围。

[0193] 公式12中，第一个式子为高速屏蔽区最小最大偏移距的设计公式，θ0c为高速层顶界的临界角，限定了有效透射信息的范围；底界的临界角为θ2c限定了高速屏蔽层下伏地层3
的有效倾角范围。以图4中给定高速屏蔽模型为例V1＝5000m/s，ρ1＝2.5g/cm ，上覆地层速
3 3
度V0＝2000m/s，ρ0＝1.8g/cm ，下伏地层速度V2＝3000m/s，ρ2＝2.1g/cm。h0＝1200m，h1＝
500m h2＝200m。在此参数条件下，则可根据公式12制作对应的高速屏蔽的偏移距设计参考量版，制作过程包括：第一步，根据12式中第一式建立高速屏蔽区的有效偏移距模型；第二步，选取实际生产中常用的偏移距范围，例如：1000m，2000m，3000m,4000m,5000m,6000m,
7000m,800m,9000m,10000m,11000m,12000m，则在图5建立的量版中可以找到各个偏移距所
对应的入射角范围，如表1。

[0194] 表1实际勘探常用偏移距对应入射角

[0195]

[0196] 将图5中的模型进一步拓展到多层介质，建立起高速屏蔽区多层上覆地层震源同入射角的关系模型如图6所示，根据Snell定律存在如下关系，假设水速V水，

[0197]

[0198] 则由公式13中提取水层与高速层上覆地层对应的公式如公式14，

[0199]

[0200] 假设水速V水＝5000m/s，ρ水＝1.03g/cm3，则可根据公式12换算出表1中各个偏移距对应的震源对应的出射角度范围，如表2。

[0201] 表2震源出射角

[0202]

[0203] 步骤5.建立待处理工区的相控阵组合模型。

[0204] 在本实施例中，现有技术中对于震源组合或检波器组合方案中，一般采用以下的相控阵组合的基距设计方法；

[0205]

[0206] 公式3是常用的组合公式，m为组合内检波器/震源顺序号，式中Dm为检波器的灵敏度或者震源的能量大小，v*为视速度，v为入射地震波速度(相对检波器)或者出射地震波速度(相对震源)，α为入射角(相对检波器)或者出射角(相对震源)，xm为检波器/震源到组合对称中心的距离，结合相控阵技术对公式(1)进行了改造，将组合响应由二维拓展到三维，公式中的xm替换为rm；将沿入射角α入射信号拓展到在三维空间由入射信号/出射信号的方位角θ0和仰角ψ0(入射角＝π/2‑ψ0或者出射角＝π/2‑ψ0)以及检波器/震源的方位角θm和仰角ψm决定；其次，将视速度v*用速度用v替代，速度v定义为水速Vwater，得到灵敏度的三维响应公式4。

[0207]

[0208] 为了满足在高速屏蔽区的震源激发的能量集中在有效入射角范围之内，加入有效入射角范围以及下传能量，获得重新构建的相控阵组合模型如公式5，

[0209]

[0210] 定义η为有效入射角范围，公式中n为组合单元的数目；Dm为组合单元的灵敏度或者能量强度；rm为各个组合单元到组合中心得距离；θ0入射信号/出射信号的方位角；ψ0为入射角或者出射角定义；以及检波器/震源的的方位角θm和入射角/出射角ψm；Vwater为水速，V0为高速层上覆地层的速度，为组合相应 的在‑π/2≤ψ0≤π/2的均匀半空间接收到/激发的总能量响应,Aη为在‑η≤ψ0≤η的入射角范围内的接收到/激发的能量响应，Aη和 的比值构成了组合响应聚焦的判别公式。如图7为24Hz信号对应着不同的基距的能量聚焦示意图，图7a中示意的是能量聚焦随着基距的变化，根据公式5中，以第一个式子中的rm做为变量，其余参数做为常量求取K，然后按照下式求取AP，得到组合响应随基距变化的曲线图7a。图
7a中示意的是能量聚焦随着基距的变化，限定的入射角范围η＝21.10°，基距从0开始增大，到达93m的时候，能量占比达到最大达到54％，随着基距的增大，能量聚焦性整体变弱，局部有振荡，存在局部能量聚焦高值点。选取第一和第二高值点的基距的组合方案进行比较，看组合响应 随入射角的变化，如图7b，实线为第一高值点的组合响应随入射角的变化，虚线为第二高值点组合响应随入射角的变化，明显看到第一高值点对应的方案能量聚焦性最
好。上面的示例是针对单一频率的信号的，实际地震勘探中，常用的频带范围分布在0～
200Hz，则根据公式5也可以拓展到宽频带的设计中。如图8所示能量聚焦宽频带设计示意
图，则可以针对实际的不同的勘探频率需求设计基距组合方案。

[0211] 公式5为一高维度设计公式，参变量多，可以进一步扩展公式5的应用维度，例如：可以将检波器组合单元数目N的设计考虑进去，得到图9的高维度设计示意，假设检波器组
合单元数目N的变化范围[4，34]，从图9能量聚焦高维度示意可以看出检波器数N＝6的情况下，检波器组合方案的整体性能最优。另外，针对公式5中某些参数可以选取某些参数进行联合设计，得到高维度的设计方案。

[0212] 步骤6.获得待处理工区预置频率中每个频率所对应的长基距。

[0213] 在本实施例中，长基距是指长度大于25m的基距。获得待处理工区预置频率中每个频率所对应的长基距的实现过程包括：

[0214] 第一步：根据所述待处理工区的相控阵组合模型，计算得到预置频率相对应的每个频率的能量曲线，其中，所述预置频率是根据所述待处理工区的地震资料分析设置的低
频段频率。

[0215] 在本步骤中，基于步骤5所创建的相控阵组合模型，结合待处理工区的地震资料情况，分析该地区地震资料的低频段范围，并选择相应的预置频率进行确定该预置频率相对
应的能量曲线，并进一步确定预置频率中所选择的每个频率所对应的基距。例如：在火成
岩，灰岩，岩盐等存在的高速层屏蔽区，在实际的待处理工区的地震资料分析发现，高速层上下地震信号的频带发生突变，尤其是下伏地层的频带分布在3～24Hz，导致频带发生突变的根本原因是由于高速屏蔽层对高频信号形成了遮挡。在本步骤中，采用所建立的相控阵
组合模型确定预置频率所对应的能量曲线。在3～24Hz频段内，选取四个低频频率f1，f2，f3和f4进行基距设计，例如：低频频率选取的频率f1＝3Hz，f2＝6Hz，f3＝12Hz和f4＝24Hz；对于该低频频率的选择可以根据相应的待处理工区的地震资料的品质进行分析选择相应的
频率。根据待处理工区的相控阵组合模型求取四个低频频率f1，f2，f3和f4的能量随基距变化曲线如图10所示。

[0216] 第二步：根据能量曲线确定每个频率对应的长基距。

[0217] 在本步骤中，根据第一步所获得的每个预置频率中每个频率所对应的能量曲线，选择每个频率的能量曲线中多个波峰的位置所对应的基距作为待确定的基距；对多个待确
定的基距进行对比，选择能量聚焦的基距作为该频率对应的基距。具体实现过程如下：根据第一步中获得的每个频率的能量曲线，选取第一和第二高值点的基距的组合相应进行比
较，分析组合响应 随入射角的变化，如图7b所示，在图7b中实线为第一高值点(第一高值
点对应第一个波峰)的组合响应随入射角的变化，虚线为第二高值点(第二高值点对应第二
个波峰)的组合响应随入射角的变化，从对比分析来看，第一高值点对应的组合响应能量聚焦性最好。根据上述的确定过程，分别确定四个频率对应的基距为Lf1＝748m，Lf2＝372m，Lf3＝187m，Lf4＝93m。在常规采集中震源和检波器采用的组合方式所使用的基距，震源组合常用的基距在18m左右，而检波器组合常用的基距在12.5m左右。根据低频能量聚焦设计的基
距要大于常规采集所采用的基距，因此，针对低频段所确定的基距定义为长基距。长基距是指基距长度大于25m的基距。

[0218] 步骤7.针对该长基距设置所述每个频率所对应的组合单元的组内距。

[0219] 在本实施例中，组合单元为由n个检波器/震源组成的检波器/震源组合单元，n为大于1的正整数。

[0220] 一种示例性的实施例中，针对每个频率所对应的长基距，采用预先设置的检波器/震源的组合单元数量，分别设计每个频率所对应的检波器/震源组合单元的组内距。在本实施例中，得到该待处理工区预置频率的检波器/震源的空间分布可以是优化的检波器/震源
的组合方案即不等间距的检波器/震源的组合单元，也可以是等间距的检波器/震源的组合
方案即等间距的检波器/震源的组合单元。

[0221] 在本实施例中，根据所确定的每个频率的基距可以采用最优组合设计和等间距组合设计两种方案，具体的实现过程如下：

[0222] 步骤71.采用最优组合方案。

[0223] 对每个频率对应的基距，采用预先设置的震源/检波器的组合单元数量，分别设计每个频率所对应的组合单元的组内距；结合每个频率所对应的组合单元的组内距，得到该
待处理工区预置频率的震源/检波器的空间分布。

[0224] 在本实施例中，是根据所确定的每个频率分别设置组合方案中每个组合单元的组内距即选择采用最优组合方案，具体的实现过程如下所示：以检波器组合为例，假设采用组合单元数N＝8，则先针对四个基距Lf1，Lf2，Lf3和Lf4分别设计各个组合单元的组内距r。f1频率对应的各检波器组合单元中心距为r11，r12，r13，r14，r15，r16，r17，r18；f2频率对应的各检波器组合单元中心距为r21，r22，r23，r24，r25，r26，r27，r28；f3频率对应的各检波器组合单元中心距为r31，r32，r33，r34，r35，r36，r37，r38；f4频率对应的各单元中心距为r41，r42，r43，r44，r45，r46，r47，r48；根据所确定的每个频率的基距计算得到每个频率对应的各检波器组合组合单元中心距，如表3。

[0225] 表3基距设计表

[0226]

[0227] 由于针对不同频率所对应的基距及组内距离是随机分布的，按照由大到小的顺序进行重排组内距，通过重排组内距会合并重复出现的组内距的方案，减少重复的工作量。因此，在得到每个频率对应的各组合单元的组内距后，对表3基距进行重排得到表4。

[0228] 表4基距重排表

[0229]

[0230]

[0231] 步骤72.采用等间距设计。

[0232] 在本实施例中，在实际地震数据采集实施过程中，图11所设计的方案为各组合单元的不等间距方案，即检波器组合中各个单元的组内距是变化的。在实际的生产中，一般是采用等间距的设计，即检波器组合中各个单元的组内距是固定。因此，在得到每个低频的长基距后，需要结合实际生产中的参数来实现炮点长基距组合和检波点长基距组合，可以根
据炮点间距SI和检波器间距RI对地震数据采集参数进行优化，使得基距为SI/RI的整数倍。
具体的实现过程包括：

[0233] 步骤720.获取待处理工区中炮点间距SI和检波器间距RI；

[0234] 步骤721.采用基距规则化方法将所获取每个频率所对应的组合单元的组内距规则化，得到炮点间距SI和检波器间距RI整数倍的组内距。

[0235] 在该步骤中，将每个频率所对应的组合单元的组内距进行微分计算，得到微分结果；对所述微分结果取整计算获得微分取整结果；对所述微分取整结果进行积分计算得到
规则化的基距。具体实现过程例如：使得基距的变化为SI/RI的整数倍。具体公式如下：

[0236]

[0237] 先对基距求微分，接着对微分结果取整计算，对微分取整结果积分重构基距r'ij。在本公式中，对于炮点间距SI和检波器间距RI进行分别计算。针对表4进行基距规则化，则分别得到SI＝RI＝25m和SI＝RI＝12.5m的基距规则化后的设计表，表5和表6。

[0238] 表5 SI＝RI＝25m基距规则化设计表

[0239]

[0240] 表6 SI＝RI＝12.5m基距规则化设计表

[0241]

[0242] 对表5和表6的参数进行重排得到新的组合设计，如表7和表8所示。

[0243] 表7 SI＝RI＝25m基距规则化重排表

[0244]

[0245]

[0246] 表8 SI＝RI＝12.5m基距规则化重排表

[0247]

[0248]

[0249] 上述的炮点间距SI和检波器间距RI采用了25m间距和12.5m间距，炮点间距SI和检波器间距RI采用了25m间距和12.5m间距是目前海上常用的两种间距，假设间距为：SI＝RI
＝6.25m，SI＝RI＝3.125m，SI＝RI＝1.5625m和SI＝RI＝1m三种间距的长基距组合方式。所采用的频率采用和上述步骤中的例子保持一致即：f1＝3Hz，f2＝6Hz，f3＝12Hz和f4＝
24Hz。将上述4个方案基距规则化重排表按照上述的规则化方法获得，SI＝RI＝6.25m基距
规则化重排表、SI＝RI＝3.125m重排表、SI＝RI＝3.125m重排表、SI＝RI＝1.5625m重排表和SI＝RI＝1m重排表。重排的信息如下表所示：表9为SI＝RI＝6.25m重排表，表10为SI＝RI＝3.125m重排表，表11为SI＝RI＝1.5625m重排表，表12为SI＝RI＝1m重排表。

[0250] 表9 SI＝RI＝6.25m基距规则化重排表

[0251]

[0252]

[0253] 表10 SI＝RI＝3.125m基距规则化重排表

[0254]

[0255]

[0256] 表11 SI＝RI＝1.5625m基距规则化重排表

[0257]

[0258]

[0259] 表12 SI＝RI＝1m基距规则化重排表

[0260]

[0261]

[0262] 对于上述所提到的方案表4中理想的组合方案，及SI＝RI＝25m，SI＝RI＝12.5m，SI＝RI＝6.25m，SI＝RI＝3.125m，SI＝RI＝1.5625m和SI＝RI＝1m情况下的组合方案。这里的连等不只限定同时采用的参数，炮点和检波点组合方案分别的计算和设计也包含在内，
针对于炮点和检波器的间距可以分别采用不同的参数设计。

[0263] 步骤8.根据所述每个预设频率所对应的组合单元的组内距，确定检波器/震源的空间分布。

[0264] 在本步骤中，根据所确定的每个预设频率所对应的组合单元的组内距，确定检波器/震源的空间分布。例如：表4即为兼顾f1，f2，f3和f4四个频率理想的组合方案，在获得表
4的组合方案后，获得如图11所示兼顾f1，f2，f3和f4四个频率的理想的组合方案的空间分布状态。

[0265] 步骤9.根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布。

[0266] 采集方法实现主要包括实际采集中，按照上述提出的组合方案的设计实现的采集方法(主要包括按照所提出的方案设计的炮线距，炮点距，检波线距以及检波点距)，以及通过上述的基础方案进行空间内插得到的简化工作量的方案。采集方法限制的主要包括按照
上述方案设计的不等间距组合方案，如图12a中，采用表5中针对f3＝12Hz频率增强的方案，布设了8条测线，测线间距呈不等间距，类似用上述参数定义的观测系统，也可以拓展到其它频率增强所采用的不等间距的方案，如图12b，图中实线为检波线，点线为炮线，不仅限于此，也可以采用重排后多个频率段相互累加的方案，图12c方案所示。采集数据后，对应8条线的炮点数据可采用炮点垂直叠加和炮点静校正后垂直叠加重构单炮数据。

[0267] 针对于上述步骤所得到的基于长基距的采集系统的实施工作量较大，需要通过简化工作量，基于所提出的长基距的采集系统方案进行空间内插重新得到的组合方案。

[0268] 步骤91.获待处理工区预置频率的检波器/震源的空间分布后，将所有的频率所对应的检波器/震源组合单元的组内距进行排序；选取检波器/震源组合单元的组内距中最小
值和最大值所对应的组内距，按照缩减工作量的比例和缩减公式计算得到缩减的检波器/
震源组合单元的分布。

[0269] 步骤911.获待处理工区预置频率的检波器/震源的空间分布后，将所有的频率所对应的检波器/震源组合单元的组内距进行排序；

[0270] 步骤912.选取检波器/震源组合单元的组内距中最小值和最大值所对应的组内距，按照缩减工作量的比例和缩减公式计算得到缩减的检波器/震源组合单元的分布；

[0271] 其中，所述缩减公式，包括：

[0272]

[0273] 其中，rji为原基距，其中j为频率的编号，i为组合单元的编号，min(rji)为原基距中的最小基距，max(rji)为原基距中的最大基距，N为原组合单元的数目，缩减工作量的比例1/M，M为2的整数倍，M＝2，4，8，16.......，k为缩减组合单元的顺序号，k＝1,2,
3,...........。采用上述缩减的步骤，图11示意的是在基于表4理想方案，进行空间内插得到的简化工作量的方案，如表13所示：

[0274] 表13基于表4理想方案工作量简化表

[0275]

[0276]

[0277] 根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布，还包括：步骤92.分别针对每个频率的检波器/震源的组合单元进行重排列，按照缩减工作量的比例和缩减公式计算得到对应的该频率的缩减的检波
器/震源组合单元的分布；获得所有频率缩减的检波器/震源组合单元的分布后，对所有频
率的缩减检波器/震源组合单元进行重新排列；根据重新排列的缩减的检波器/震源组合单
元得到缩减的检波器/震源组合单元的分布。

[0278] 步骤921.分别针对每个频率的检波器的组合单元进行重排列，按照缩减工作量的比例和缩减公式计算得到对应的该频率的缩减组合单元的分布；

[0279] 步骤922.获得所有预置频率的缩减组合单元的分布后，对所有预置频率的缩减组合单元进行重新排列。

[0280] 采用上述步骤92对于组合方案进行缩减工作量，首先对于每个频率分别进行工作量缩减，而后在对组合方案进行重排，如以表3中的四个频率f1，f2，f3和f4所对应的工作量缩减为例，如表14为f1频率组合方案缩减方案，表15为f2频率组合方案缩减方案，表16为f3频率组合方案缩减方案，表17为f4频率组合方案缩减方案。

[0281] 表14 f1频率组合方案缩减方案

[0282]

[0283] 表15 f2频率组合方案缩减方案

[0284]

[0285] 表16 f3频率组合方案缩减方案

[0286]

[0287] 表17 f4频率组合方案缩减方案

[0288]

[0289] 表14到表17组合方案中，则针对某一频率增强的组合方案算上原方案有三种(原方案+工作量缩减1/2方案+工作量缩减1/4方案)，此时对多频率段进行增强组合设计，如同时采用f1+f2+f3+f4组合方案，则根据排列组合公式，可生成的检波器组合方案81种方案。
实际中可根据低频段频率增强的需求，在81中方案中选择满足要求的方案。如表18为根据
表14到表17四个频率同时缩减组合工作量1/4的方案示意。

[0290] 表18 f4频率组合方案缩减方案

[0291]

[0292] 图13为缩减工作量前后组合单元的空间位置及组合单元数目的变化，此方案中组合单元数由32降为8。在得到缩减的检波器/震源组合单元的分布后，可以根据该缩减的检
波器/震源组合单元的分布进行地震数据采集。

[0293] 通过本发明方法，根据所述检波器/震源的空间分布，采用预定规则缩减组合单元的个数，得到缩减的检波器/震源组合单元的分布；根据缩减的检波器/震源进行采集组合
单元的分布进行地震数据采集。通过本发明的方案，在基于长基距地震数据采集时，可以缩减地震数据采集的工作量。

[0294] 本领域普通技术人员可以理解，上文中所公开方法中的全部或某些步骤、系统、装置中的功能模块/单元可以被实施为软件、固件、硬件及其适当的组合。在硬件实施方式中，在以上描述中提及的功能模块/单元之间的划分不一定对应于物理组件的划分；例如，一个物理组件可以具有多个功能，或者一个功能或步骤可以由若干物理组件合作执行。某些组件或所有组件可以被实施为由处理器，如数字信号处理器或微处理器执行的软件，或者被
实施为硬件，或者被实施为集成电路，如专用集成电路。这样的软件可以分布在计算机可读介质上，计算机可读介质可以包括计算机存储介质(或非暂时性介质)和通信介质(或暂时
性介质)。如本领域普通技术人员公知的，术语计算机存储介质包括在用于存储信息(诸如
计算机可读指令、数据结构、程序模块或其他数据)的任何方法或技术中实施的易失性和非易失性、可移除和不可移除介质。计算机存储介质包括但不限于RAM、ROM、EEPROM、闪存或其他存储器技术、CD‑ROM、数字多功能盘(DVD)或其他光盘存储、磁盒、磁带、磁盘存储或其他磁存储装置、或者可以用于存储期望的信息并且可以被计算机访问的任何其他的介质。此
外，本领域普通技术人员公知的是，通信介质通常包含计算机可读指令、数据结构、程序模块或者诸如载波或其他传输机制之类的调制数据信号中的其他数据，并且可包括任何信息
递送介质。