[发明专利]一种保护激发频率的延时震源组合参数计算方法

摘要

本发明涉及一种保护激发频率的延时震源组合参数计算方法，该方法首先需要建立地质模型，在后期优化可能用到的坐标位置上布设震源，在以最浅目的层深度为半径组合中心为球心的球面上等角度布设检波器，完成所有震源的单炮正演；然后在给定组合个数和组合形态的基础上通过数值寻优算法或贪婪算法计算得到最佳组合参数，计算过程需要满足给定的主频约束条件。本发明克服了当前震源组合参数设计上的某些不足，寻找到了一种快速的组合参数定量化计算方法，同时具有良好地可操作性，可以为野外施工提供方案指导。

主权项

1.一种保护激发频率的延时震源组合参数计算方法，其特征在于：包括以下步骤：步骤(1)、根据工区实际情况，建立三维地质模型；步骤(2)、基于三维地质模型，确定待优化的组合中心，围绕该中心在三维地质模型的地质网格节点上布设一定数量的震源；步骤(3)、获取最浅目的层深度信息以步骤(2)确定的组合中心为球心，以最浅目的层深度为半径做球面，该球面与地质网格相交，在球面上以等倾角和水平角度间隔布设检波器，以获取最浅目的层深度信息；步骤(4)、对步骤(2)中所有震源依次进行单程波动方程正演，记录下各检波器接收到的首波波形，存储待用；步骤(5)、根据检波器接收排列长度和最浅勘探目标层的深度确定有效入射角α，根据组合激发图形中心点的深度和最浅勘探目标层的深度确定非有效入射角β；步骤(6)、统计球面上投射至目的层锥角范围内的能量和其它方向能量的比值，以能量比值最大化为判别标准进行优化计算根据步骤(5)得到的有效入射角α和非有效入射角β，将给定震源组合个数和平面组合形态作为已知量，统计球面上投射至目的层锥角范围内的能量和其它方向能量的比值，以能量比值最大化为判别标准，采用粒子群优化算法或贪婪算法计算延时震源组合参数，延时震源组合参数包括各个震源的三维空间坐标(xi,yi,zi)、组合加权系数Ci和组合激发延迟时间△t延迟i。