1.一种基于Lq范数约束的用于地震勘探的速度分析方法,其特征在于,包括以下步骤:采集用于地震勘探的原始数据,进行第一坐标变换后,获取坐标变换数据;
对所述坐标变换数据进行双曲Radon变换,获取Radon域数据;
对所述坐标变换数据和所述Radon域数据,施加基于Lq范数的稀疏约束条件后,进行第二坐标变换,生成用于地震勘探的速度谱。
2.根据权利要求1所述一种基于Lq范数约束的用于地震勘探的速度分析方法,其特征在于:2
在进行第一坐标变换的过程中,对原始数据D(t,x)进行第一坐标变换t'=t ,获得其坐标变换后的数据D′(t′,x),其中,2
D′(t′,x)=D(t,x)。
3.根据权利要求2所述一种基于Lq范数约束的用于地震勘探的速度分析方法,其特征在于:在获取Radon域数据的过程中,所述Radon域数据表示为:
4.根据权利要求3所述一种基于Lq范数约束的用于地震勘探的速度分析方法,其特征在于:在施加基于Lq范数的稀疏约束条件的过程中,所述稀疏约束条件的目标函数表示为:其中,FFT和IFFT分别表示正、逆傅里叶变化,R表示Radon变换,λ为正则化参数,|| ||1为求取L1范数, 表示求取Lq范数,即 a,b分别代表D的总行数和列数。
5.根据权利要求4所述一种基于Lq范数约束的用于地震勘探的速度分析方法,其特征在于:在施加基于Lq范数的稀疏约束条件的过程中,所述稀疏约束条件的施加过程,包括以下步骤:(1)设置
(2)通过交替方向乘子算法进行求解目标函数获得约束后的Radon域数据
(3)判断是否满足条件 ε为截至阈值。如果满足,则停止循环,输出结果
6.根据权利要求5所述一种基于Lq范数约束的用于地震勘探的速度分析方法,其特征在于:在进行第二坐标变换的过程中,通过对约束后的输入结果 进行两次坐标变换,生成所述速度谱。
7.根据权利要求6所述一种基于Lq范数约束的用于地震勘探的速度分析方法,其特征在于:在生成速度谱的过程中,对约束后的输入结果 进行第一次坐标变换获得第一次坐标变换后的数据
对数据 进行第二次坐标变换 获得第二次坐标变换后的数据根据数据 获取所述速度谱。
8.根据权利要求7所述一种基于Lq范数约束的用于地震勘探的速度分析方法,其特征在于:在获取速度谱的过程中,所述速度谱表示为:‑1
其中, R、R 为正、逆Radon变换,w为时间窗口,N为总道数,S(t0,v)表示速度谱。
9.根据权利要求8所述一种基于Lq范数约束的用于地震勘探的速度分析方法,其特征在于:用于实现速度分析方法的速度分析系统包括:数据采集模块,用于采集所述原始数据;
第一数据处理模块,用于对所述原始数据,进行第一坐标变换,获取坐标变换数据;
第二数据处理模块,用于对所述坐标变换数据进行双曲Radon变换,获取Radon域数据;
速度分析模块,用于对所述坐标变换数据和所述Radon域数据,施加基于Lq范数的稀疏约束条件后,进行第二坐标变换,生成用于地震勘探的速度谱。
10.根据权利要求9所述一种基于Lq范数约束的用于地震勘探的速度分析方法,其特征在于:用于实现速度分析方法的速度分析系统,还包括:地震勘探模块,用于根据所述速度谱进行地震勘探。