[发明专利]一种基于大炮初至的折射波剩余静校正方法

说明书

技术领域

本发明涉及地球物理勘探技术，是一种基于大炮初至的折射波剩余静校正方法。

背景技术

静校正是对地震资料所作的校正，用于补偿由高程、风化层厚度以及风化层速度产生的影响，把资料校到一个指定的基准面上。其目的通常是获得在一个平面上进行采集，且没有风化层或低速介质存在时的反射波到达时间。它是一种地震勘探中消除近地表影响从而提高地震反射波成像质量的主要数据处理方法。

2004年井西利等人公开了一种自适应混合反演剩余静校正的方法，针对剩余静校正问题的非线性特点，分析了已有的线性反演及非线性反演方法的优缺点，在此基础上提出了一种计算效率较高、搜寻全局最优解的混合反演方法。即在求解过程中每次扰动采用线性寻优方法搜寻局部最优解，在迭代过程中采用模拟退火法随机搜寻，就可以用较少的迭代次数搜寻全局最优解。还提出了反演参数的选取方法，使得该法具有自适应的特点。理论分析和数值计算表明，该方法不仅解决了剩余静校正较强的非线性全局优化问题，而且提高了计算效率。

2005年李全胜等人公开了一种基于模型道的地表一致性剩余静校正方法，论述了模型道的构造方法，将剩余静校正问题作为广义线性反问题来求解炮点和检波点剩余静校正量。通过实际资料应用效果表明，这种方法能够很好地解决较复杂地区地震资料的剩余静校正问题。

2006年段云卿公开了一种折射波剩余静校正方法，利用折射波信噪比高的特点，将反射波剩余静校正原理应用于折射波处理，通过交互手段，逐段估算折射波的速度，用合适的速度对地震记录进行线性动校正，在共炮点或共中心点道集上，用相关方法计算各道与模型道时差，再用统计方法计算出各炮点及检波点的剩余静校正量。

上述技术都依赖地震处理系统，即需要在对地震资料进行一定处理后才能应用，所要计算的数据量大，要求地震资料具有较高的信噪比，以便于模型道的制作和相关计算的时差准确，三种方法对软件和硬件都提出了较高要求，如需要配置工作站和资料处理系统，一般只能在资料处理中心进行，野外现场无法在微机上实际应用。

发明内容

本发明的目的是提供一种针对在复杂地区折射层复杂和地震资料信噪比较低的情况下，在现场微机上计算处理，提高基准面静校正精度的基于大炮初至的折射波剩余静校正方法。

本发明通过以下具体步骤实现：

1)利用表层资料或者大炮初至资料进行基准面静校正；

步骤1)所述的基准面静校正对于表层资料是采用小折射或微测井对勘探工区进行表层结构调查，得到各调查点的表层结构数据，采用层间关系系数建立表层模型；对于大炮初至资料是先拾取野外生产的大炮初至时间，采用折射反演或层析反演方法建立表层模型，在统一定义水平基准面和填充速度后，由建立的表层模型数据计算基准面静校正。

2)对拾取的所有地震道的初至时间数据施加该地震道对应的炮点和检波点的基准面静校正，生成新的地震道的初至时间数据T_s，j；

步骤2)所述的基准面静校正是将每一地震道的初至时间加上该地震道对应的检波点的静校正和炮点的静校正。

3)生成共检波点道集对应的地震道的初至时间数据T_r，k；

步骤3)所述的生成地震道的初至时间数据是按炮点的地面位置进行顺序排序，抽取该炮中含有某指定检波点的地震道的初至时间，按正负炮检距的大小进行排序，得到某共检波点道集的初至时间。

4)在所有共炮点道集的正负偏移距上根据折射速度的大小和分布范围分别划分出一个折射层的范围，对每个正负偏移距上指定折射层范围内各检波点的偏移距与初至时间分别进行线性拟合或二次曲线拟合，得到每个炮点道集内各检波点的拟合初至时间拟合前初至时间为T_s，j，由用拟合初至时间减去拟合前的初至时间得到每个道集指定折射层范围内各检波点处的初至拟合时差；

步骤4)所述的各检波点处的初至拟合时差为Δt_s，j＝T′_s，j-T_s，j，

Δt_s，j为s炮点时j检波点处地震道的初至拟合时差，用拟合初至时间减去拟合前的初至时间得到；

T′_s，j为s炮点时j检波点处采用数据拟合计算的地震道的拟合初至时间；

T_s，j为s炮点时j检波点处地震道拟合前的初至时间。