System.ArgumentOutOfRangeException: 索引和长度必须引用该字符串内的位置。 参数名: length 在 System.String.Substring(Int32 startIndex, Int32 length) 在 zhuanliShow.Bind()
【技术实现步骤摘要】
本申请属于地震数据处理领域,具体涉及一种地震数据层析静校正方法、装置、电子设备及存储介质。
技术介绍
1、复杂地表引起的静校正问题对地震采集数据和数据处理都有严重的影响,静校正问题的存在使得数据处理中的部分地球物理假设不满足,给数据处理带来很大的困难。在复杂探区,近地表低降速带巨厚、测点高程起伏大、速度横向差异大,静校正问题比较难解决。针对二维测线的静校正问题处理中,主要的方法有模型法静校正、高程静校正和层析静校正,在生产应用中主要采用的是二维测线层析静校正技术,该方法能充分利用各类型波的初至信息,可以得到相对比较准确、稳定、符合射线理论的单线近地表速度模型,计算得到静校正量,从而能够的解决单条二维测线的静校正问题。
2、在沙漠工区存在巨大沙丘、信噪比低、高程起伏大、地表条件复杂、观测系统变化大等问题,因此采用层析静校正的方法处理静校正问题。常规的二维测线层析静校正处理流程中,分别对每条线单独进行反演、建立近地表速度模型、计算静校正量、计算生成浮动基准面,但由于空间采样不足,仅考虑沿线方向的信息,增加了反演结果的不确定性,特别是缺少近偏移距的数据时,由于层析方法固有的缺陷,导致反演结果不可靠,而且每条测线单独静校正得到的浮动基准面有差异,不能做到全区统一浮动基准面处理,影响处理效果;因反演缺少纵向信息,对交点处的成像也有影响。
技术实现思路
1、基于以上技术问题,本申请提出一种地震数据层析静校正方法、装置、电子设备及存储介质。
2、第一方面,本申请提出一种地
3、步骤s1:建立目标工区的近地表速度模型体,并得到地下介质慢度场;
4、步骤s2:利用所述近地表速度模型体计算初至时间矩阵;
5、步骤s3:利用初至时间残差矩阵进行层析反演,得到地下介质慢度场的更新量,所述初至时间残差矩阵为所述初至时间矩阵与实测的初至时间矩阵之间的差值;
6、步骤s4:如果更新量更新后的地下介质慢度场不满足精度要求,则利用所述更新量更新近地表速度模型体,转到步骤s2;
7、步骤s5:如果更新量更新后的地下介质慢度场满足精度要求,则根据当前的近地表速度模型体,计算静校正量。
8、所述地下介质慢度场是将所述近地表速度模型体中的每一个位置点的速度取倒数得到。
9、所述利用所述近地表速度模型体计算初至时间矩阵,包括:
10、根据射线路径矩阵以及所述地下介质慢度场,计算初至时间矩阵。
11、所述初至时间残差矩阵的计算方式如下:
12、将初至时间的二维测线合并为实测的初至时间矩阵;
13、将所述初至时间矩阵与所述实测的初至时间矩阵之间的差值,作为初至时间残差矩阵。
14、所述利用所述更新量更新近地表速度模型体,包括:
15、将所述更新量与所述近地表速度模型体相加,得到更新后的近地表速度模型体。
16、所述根据射线路径矩阵以及所述地下介质慢度场,计算初至时间矩阵,计算式如下:
17、t=ab
18、其中,a为射线路径矩阵,b为地下介质慢度场,t为初至时间矩阵。
19、所述利用初至时间残差矩阵进行层析反演,得到地下介质慢度场的更新量,计算式如下:
20、aδb=δt
21、其中,a为射线路径矩阵,δb为地下介质慢度场的更新量,δt初至时间残差矩阵。
22、所述根据当前的近地表速度模型体,计算静校正量,计算式如下:
23、
24、其中,st(j)为第j个检波点或者炮点的静校正量,hj为第j点位置处的低降速带厚度,从当前的近地表速度模型体中提取,vj为第j点位置处的低降速带速度,从当前的近地表速度模型体中提取,ed为处理基准面高程,ej为第j个检波点或者炮点的高程,v0为处理替换速度。
25、所述地震数据层析静校正方法,在计算静校正量之后,还包括:
26、在目标工区内定义网格面;
27、将所有静校正量组成的曲线进行平滑处理,得到浮动基准面;
28、在所述网格面中每一个网格中,基于所述浮动基准面对所述静校正量进行高低频量分离,得到所述静校正量的高频量和低频量;
29、基于浮动基准面,利用所述静校正量的低频量,将所述静校正量校正到预定义的最终基准面。
30、所述静校正量的低频量为静校正量对应浮动基准面垂点到预定义的最终基准面垂直距离值;静校正量的高频量为静校正量对应浮动基准面的垂点与所述静校正量之间的距离值。
31、第二方面,本申请提出一种地震数据层析静校正装置,包括:模型体建立模块、初至时间计算模块、层析反演模块、模型体更新模块、静校正量计算模块;
32、所述模型体建立模块用于建立目标工区的近地表速度模型体,并得到地下介质慢度场;
33、所述初至时间计算模块用于:利用所述近地表速度模型体计算初至时间矩阵;
34、所述层析反演模块用于利用初至时间残差矩阵进行层析反演,得到地下介质慢度场的更新量,所述初至时间残差矩阵所述初至时间矩阵与实测的初至时间矩阵之间的差值;
35、所述模型体更新模块用于如果更新量更新后的地下介质慢度场不满足精度要求,则利用所述更新量更新近地表速度模型体,转到初至时间计算模块;
36、所述静校正量计算模块用于如果更新量更新后的地下介质慢度场满足精度要求,则根据当前的近地表速度模型体,计算静校正量。
37、第三方面,本申请提出一种电子设备,包括:一个或多个处理器,以及存储器,所述存储器存储指令,当所述指令被所述一个或多个处理器执行时,使得所述一个或多个处理器执行所述的地震数据层析静校正方法。
38、第四方面,本申请提出一种计算机可读存储介质,其存储有可执行指令,所述指令当被执行时使得机器执行所述的地震数据层析静校正方法。
39、有益效果:
40、本申请通过是在二维层析静校正的基础上提出了拟三维层析静校正,充分利用交点处不同方向的初至信息以及考虑到炮点和检波线不在同一直线的情况,做到全区统一数据库、统一分解得到浮动基准面,解决了全区统一浮动基准面的问题,同时改善了观测系统变化处和交点处的成像。
本文档来自技高网...【技术保护点】
1.一种地震数据层析静校正方法,其特征在于,包括:
2.如权利要求1所述的地震数据层析静校正方法,其特征在于,所述地下介质慢度场是将所述近地表速度模型体中的每一个位置点的速度取倒数得到;所述利用所述近地表速度模型体计算初至时间矩阵,包括:
3.如权利要求1所述的地震数据层析静校正方法,其特征在于,所述初至时间残差矩阵的计算方式如下:
4.如权利要求1所述的地震数据层析静校正方法,其特征在于,所述利用所述更新量更新近地表速度模型体,包括:
5.如权利要求2所述的地震数据层析静校正方法,其特征在于,所述根据射线路径矩阵以及所述地下介质慢度场,计算初至时间矩阵,计算式如下:
6.如权利要求1所述的地震数据层析静校正方法,其特征在于,所述利用初至时间残差矩阵进行层析反演,得到地下介质慢度场的更新量,计算式如下:
7.如权利要求1所述的地震数据层析静校正方法,其特征在于,所述根据当前的近地表速度模型体,计算静校正量,计算式如下:
8.如权利要求1所述的地震数据层析静校正方法,其特征在于,所述地震数据层析静校
9.如权利要求8所述的地震数据层析静校正方法,其特征在于,所述静校正量的低频量为静校正量对应浮动基准面垂点到预定义的最终基准面垂直距离值;静校正量的高频量为静校正量对应浮动基准面的垂点与所述静校正量之间的距离值。
10.一种地震数据层析静校正装置,其特征在于,包括:模型体建立模块、初至时间计算模块、层析反演模块、模型体更新模块、静校正量计算模块;
11.一种电子设备,其特征在于,包括:一个或多个处理器,以及存储器,所述存储器存储指令,当所述指令被所述一个或多个处理器执行时,使得所述一个或多个处理器执行权利要求1至9中任一项所述的地震数据层析静校正方法。
12.一种计算机可读存储介质,其特征在于,其存储有可执行指令,所述指令当被执行时使得机器执行权利要求1至9中任一项所述的地震数据层析静校正方法。
...【技术特征摘要】
1.一种地震数据层析静校正方法,其特征在于,包括:
2.如权利要求1所述的地震数据层析静校正方法,其特征在于,所述地下介质慢度场是将所述近地表速度模型体中的每一个位置点的速度取倒数得到;所述利用所述近地表速度模型体计算初至时间矩阵,包括:
3.如权利要求1所述的地震数据层析静校正方法,其特征在于,所述初至时间残差矩阵的计算方式如下:
4.如权利要求1所述的地震数据层析静校正方法,其特征在于,所述利用所述更新量更新近地表速度模型体,包括:
5.如权利要求2所述的地震数据层析静校正方法,其特征在于,所述根据射线路径矩阵以及所述地下介质慢度场,计算初至时间矩阵,计算式如下:
6.如权利要求1所述的地震数据层析静校正方法,其特征在于,所述利用初至时间残差矩阵进行层析反演,得到地下介质慢度场的更新量,计算式如下:
7.如权利要求1所述的地震数据层析静校正方法,其特征在于,所述根据当前的近地表速度模型体,计算静校正量,计...
【专利技术属性】
技术研发人员:王豆豆,刘燕峰,葛中慧,杨丽莹,邹少峰,
申请(专利权)人:中国石油化工股份有限公司,
类型:发明
国别省市:
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。