本发明专利技术提供了一种自动方位角道集部分叠加方法，属于地震勘探领域。所述自动方位角道集部分叠加方法包括：(0)设置输入数据参数和输出数据参数；(1)自动方位角地震数据抽取；(2)自动入射角地震数据抽取与部分叠加。利用本发明专利技术方法，能彻底解决地震方位角道集部分叠加过程中叠加次数的不均衡性，消除采集因素的影响，提高道集质量。
【技术实现步骤摘要】
本专利技术属于地震勘探领域，具体涉及，对叠前 道集进行自动分方位角、入射角部分叠加，获得规则化方位角道集。
技术介绍
由于地震方位AV0各向异性分析具有精度高、入射角范围要求小等特点，在裂缝 型储层地震预测中具有较重要的作用。 在地震方位AV0各向异性反演预测裂缝方位、密度中，需要利用叠前方位道集，道 集质量直接决定了各向异性参数反演精度。提高叠前方位道集质量的方法很多，其中通过 部分角道集叠加，增加叠加次数方法是最有效的方式之一。 部分角道集的现有部分叠加技术是： 首先，为了提高相邻道之间的信噪比，W减小炮检分布不匀带来的弊端，并保证在 不同方位上有足够密度的不同炮检距道集分布、较一致的叠加次数，通过扩大原CMP面元 的手段建立CMP宏面元。 其次，为使宏面元内各方位道数基本一致，W确保分方位叠加次数的一致性，同时 使各方位数据具有较统一的计算精度及误差，对十个宏面元数据进行了统计分析，W确定 各方位角范围。在方位角变化过程中，为使数据更为稳定，需要对方位角边界进行重叠。 然后，采用与方位角范围确定相同的方式确定入射角范围，并进行部分叠加。 最终，获得分方位、分入射角部分叠加道集。 现有技术中，确定方位角、入射角部分叠加范围是采用个别点统计，然后应用于全 区的方式。送种方式W角度划分为基础，能一定程度上解决部分叠加道集能量不均衡问题， 但对于复杂观测系统资料，还是不能彻底解决送种不均衡性、不可避免采集脚印的存在。
技术实现思路
本专利技术的目的在于解决上述现有技术中存在的难题，提供一种自动方位角道集部 分叠加方法，W部分叠加次数为基础，采用自动角度划分方式，彻底解决叠加次数的不均衡 性，消除采集因素的影响，提高道集质量。为获得高资料各向异性裂缝参数反演奠定基础。 本专利技术是通过W下技术方案实现的： -种自动方位角道集部分叠加方法，包括： (0)设置输入数据参数和输出数据参数； (1)自动方位角地震数据抽取； (2)自动入射角地震数据抽取与部分叠加。 所述步骤（0)中的输入数据参数包括： 宏面元中包含面元个数；unit ; 宏面元中每个面元所包含的地震道数；nstack ; 宏面兀中的总地震道数；Unstack = unit Xnstack ; 宏面元中每一个地震道的方位角值；azimuthj，其中j = 0,1，…，化stack-1表示 道序号； 宏面元中每个地震道的偏移距值；offsetj，其中j = 0,1，…，化stack-1表示道 序号； 宏面元中每一个地震道的采样点数；ns ; 宏面元中每一个地震道中每一个采样点的入射角值；incidencei, j，其中i = 0， 1，...，111131304-1表不道序号，j = 0,1，...，113-1表不义样点序号； 入射角 incidence^, j 的最大值为；incidence_max ;叠前地震数据为；data_ini,k，其中 i = 0,1，…，Unstack-1, k = 0,1，…，ns-1 ; 所述输出数据参数包括： 方位角道集方位数；azi_out_n ;[002引 方位角道集各方位的中必角；azi_outi = (；0. 5+;0 X 180/azi_out_n，其中i = 0， 1，…，azi_out_n-l为输出方位角道集方位角序号； 方位角道集入射角数：； 方位角道集的各入射中必角；inc_outi =化 5+i) X incidence_max/lnc_out_n, 其中i = 0,1，…，inc_out_n-l为输出方位角道集入射角序号； 初始化地震方位角道集数据；datai,k,.j = 0,其中i = 0,1，2，…，azi_out_n-l为 方位序号；k = 0，1，2，…，inc_out_n-l为入射角序号；j = 0，1，2，…，ns-1为采样点序 号。[003引所述步骤（1)包括 (a)计算输出方位角道集中每个方位的最大道数：(1)[003引 （b)初始化i = 0 ; (C)初始化 k = 0 ; (d)计算输出方位角扇区半径： (e)计算输出方位角扇区半径变化步长；azi_;r_st巧=azi_rk/l〇〇 ; (f)计算输出方位角扇区角度范围： (azi_outi-azi_rk，azi_outi+azi_rk] (g)统计，当j = 0,1，2，…，化stack-1时，宏面元中每一个地震道的方位角值 azimuthj值满足公式（3)的地震道总数，并记录该数目为Fold_Pa;rt_sumi; azi_outi_azi_;r_step < azimuth j《azi-OUti+azi-r-step (3) 记录所有满足公式（3)要求的地震道在宏面元中的位置序号azi_out_n〇i,t，其中 t = 0,1, ···, Fold_Part_sumi-l ； 化）如果 ρ·〇1(1_Ρ3；Γ?:_3?ιη?? > Fold_Pa;rt_sum_max-l，则跳转到步骤化）。 (;〇 如果 ρ·〇1(1_Ρ3；Γ?:_3?ιη?? < Fold_Pa;rt_sum_max-l，则： azi-fk+i = azi_rk+azi_r_step (4)[004引 ?·)计算k = k+1，返回步骤讯；[004引 似如果i > azi_out_n-l，自动方位抽取结束，固龄到步骤似甫则i = i+1，返 回步骤（d)。 所述步骤似包括： (A)计算输出方位角道集每个方位中每个入射角道集的最大道数：(目）[005引 度）对第i个输出方位角实施计算，初始化时i = 0 ;似当 t = 0,1，…，Fold_Pa;rt_sumi-l 时，每个 azi_out_n〇i, t 位置对应一个 巧嚴紐姑,_。",_嗎，值；依据鄉汾础_。气《扔从小到大的顺序，对azi_out_tn〇i,t进行排序；记录排 序后的宏面元中本方位的地震道位置；azi_out_no_soパi,t;[005引（D)对第j个采样点实施计算，初始化时j = 0 ;[005引 似对第k个入射角实施计算，初始化时k = 0 ; (巧搜索，当 t = 0，1，2，…，Fold_Pa;rt_sumi-l 时，是否有一个 。和。。娜值，满足公式化）：(6)[005引 如果当t = tt,且tt《Fold_Pa;rt_sumi-l时，满足公式化），跳转到步骤似；女口 果无法找到满足公式（6)条件的值，跳转到步骤（I); 似计算输出入射角扇区范围；虹，η) 其中 m = tt-Fold_Pa;rt_max/2 η = tt+Fold_Part_max/2[006引做计算第i输出方位角、第k入射角、第j采样点的输出地震数据：口）[00财 讯如果k《inc_out_n-l，则k = k+1，返回步骤胆）；否则跳转到步骤（J);[006引 α)如果j《ns-1，则j = j+1，返回步骤做；否则，固龄到步骤似； 似如果i《azi_out_n-l，则i = i+1，返回步骤做；否则，固龄到步骤（L); (L)输出的datai,k,.,即为所求的地震方位角道集。[006引所述步骤似中是采用冒泡排序算法依据0斯巧。从小到大的顺序，对本文档来自技高网...
【技术保护点】
一种自动方位角道集部分叠加方法，其特征在于：所述方法包括：(0)设置输入数据参数和输出数据参数；(1)自动方位角地震数据抽取；(2)自动入射角地震数据抽取与部分叠加。
【技术特征摘要】
【专利技术属性】
技术研发人员：肖鹏飞，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。