【技术实现步骤摘要】
快速编码免排序基准面校正方法及系统
本专利技术涉及地球物理勘探领域，更具体地，涉及一种快速编码免排序基准面校正方法及系统。
技术介绍
在反射波地震勘探中，叠前去噪、速度分析等处理方法是基于水平地表及水平层状近地表假设条件的，对于简单起伏地表、低降速带及速度均匀变化的近地表条件，可以通过常规静校正技术解决近地表问题，恢复反射波的几何形态。但是对于复杂近地表条件，如剧烈的起伏变化、强横向变化的低降速带，地表一致性的假设往往无法满足，如地表对同一位置造成的影响是随着地震波传播路径的不同而变化的，无法通过双程时差校正解决这一问题，此时，常规静校正技术无法完全解决近地表问题。基于波动方程的基准面校正，正是针对此类复杂近地表问题的有效处理手段，其思想是先将地震波场向下延拓到高速层的顶界面，再基于该工区给定的替换速度将地震波场向上延拓到最终的基准面，这个过程可以使得面波、反射波等地震波场在校正后的剖面中恢复线性传播特征，消除地表起伏、低降速带、复杂构造等复杂近地表速度分布带来的不利影响，为现有的基于水平观测面的叠前去噪技术提供优质...
【技术保护点】
1.一种快速编码免排序基准面校正方法，其特征在于，包括：/n步骤1：在共炮点记录中拾取全偏移距初至波，获得近地表速度模型，根据反演的射线密度确定波场延拓范围；/n步骤2：在所述波场延拓范围内，读入全部的共炮点地震数据，获得多个初始共炮点道集；/n步骤3：针对每一个初始共炮点道集，根据所述近地表速度模型，使用分步傅里叶方法将检波点延拓到高速层顶界面，进而根据替换速度将检波点向上延拓至基准面，获得多个共检波点道集；/n步骤4：针对每一个共检波点道集，使用分步傅里叶方法将炮点延拓到高速层顶界面，进而根据替换速度向上延拓至基准面，获得多个最终共炮点道集。/n
【技术特征摘要】
1.一种快速编码免排序基准面校正方法，其特征在于，包括：
步骤1：在共炮点记录中拾取全偏移距初至波，获得近地表速度模型，根据反演的射线密度确定波场延拓范围；
步骤2：在所述波场延拓范围内，读入全部的共炮点地震数据，获得多个初始共炮点道集；
步骤3：针对每一个初始共炮点道集，根据所述近地表速度模型，使用分步傅里叶方法将检波点延拓到高速层顶界面，进而根据替换速度将检波点向上延拓至基准面，获得多个共检波点道集；
步骤4：针对每一个共检波点道集，使用分步傅里叶方法将炮点延拓到高速层顶界面，进而根据替换速度向上延拓至基准面，获得多个最终共炮点道集。
2.根据权利要求1所述的快速编码免排序基准面校正方法，其中，所述步骤2还包括：
在所述波场延拓范围内，根据地震道中的炮、检点坐标对炮点和检波点进行统计，统计内容包括：
炮点编码，记录为Si，其中i＝1,2,3，…n，n为炮点总个数；
检波点编码，记录为Gj，其中j＝1,2,3,…m，m为检波点的总个数；
炮点覆盖次数，即每一炮涉及的地震道数，记录为SFi；
检波点覆盖次数，即每一个检波点涉及的炮点道数，记录为GFj；
根据统计的炮点编码和检波点编码，在磁盘上为每一个炮点和检波点创建一个空文件：为每一个炮点创建一个文件，命名为FSi；为每一个检波点创建一个文件，命名为FGj。
3.根据权利要求1所述的快速编码免排序基准面校正方法，其中，所述分步傅里叶方法为：
其中，为延拓前的波场，为延拓后的波场，Δz为延拓的步长，为……，u0为平均速度的倒数，kr为波数域中的水平波数向量，zn为第n步波场延拓时所在的深度，ω为波场在频率域中的一个角频率，i为虚数单位，为根据公式计算得到的垂直方向上的波数。
4.根据权利要求3所述的快速编码免排序基准面校正方法，其中，所述分步傅里叶方法还包括：
根据实际速度与平均速度的差值，通过公式(2)对所述延拓后的波场进行微调：
Pshift(r,zn+1,ω)＝e-iωΔu(r,z)ΔzP(r,zn,Δz,ω)(2)
其中，Pshift为微调后的波场，Δu(r,z)为实际速度的道数与平均速度的倒数的差值，r为水平位置向量，zn为第n步波场延拓时所在的深度，P表示在时间-空间域中zn深度上的波场。
5.根据权利要求1所述的快速编码免排序基准面校正方法，其中，将所述初始共炮点道集或所述共检波点道集分配给多个计算节点进行并行式计算。
6.根据权利要求5所述的快速编码免排序基准面校正方法，其中，进一步包括：
每个计算节点依次判断所述初始共炮点道集或所述共检波点道集是否在...
【专利技术属性】
技术研发人员：徐蔚亚，朱成宏，谢飞，赵艳平，魏哲枫，高鸿，张春涛，魏巍，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油勘探开发研究院，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。