基于波场延拓原理的起伏地表组合震源波场定向方法
【专利摘要】本发明提供了一种基于波场延拓原理的起伏地表组合震源波场定向方法。首先，针对复杂起伏地表设计一个水平基准面，将不同高程位置的震源，投影至此水平基准面；其次，根据水平地表组合震源的方向因子公式，设计组合震源的延时激发参数；再次，对于每个延时激发的震源波场记录以自激自收方式在频率波数域正向外推至实际起伏地表高程；最后，将延拓后的各个震源波场进行干涉叠加，得到定向传播的地震波场。经过数值测试，本方法可以得到稳定良好的定向波场，起伏地表对定向波场影响很小，可作为起伏地表条件下组合震源数值计算的技术方法。
【专利说明】
基于波场延拓原理的起伏地表组合震源波场定向方法
技术领域
[0001] 本发明专利属于勘探地震领域，涉及一种地震勘探组合震源的定向方法，尤其是 在复杂起伏地表条件下的地震波场定向。
【背景技术】：
[0002] 地震勘探方法是探测地下地质体构造的一种地球物理勘探方法，其重要的环节之 一是地震数据的采集。而数据采集质量的好坏与震源激发有密切关系，野外勘探时常用组 合震源的方式以提高震源激发的能量和地下目标体反射信号的信噪比。
[0003] 目前组合震源可考虑以下几种情况:水平地表条件下简单线性组合震源可用组合 原理的方向因子公式来确定震源的激发延时;倾角固定的倾斜地表条件下有类似的方向因 子公式;复杂起伏地表条件下可用旋转坐标方法的几何路径计算激发延时。如巩向博等，地 球物理学报2014, "起伏地表组合震源地震波场定向方法"。然而，对于非均匀复杂起伏地表 情况，通过几何关系计算的组合震源激发延时有一定的误差。
【发明内容】
：
[0004] 本发明针对现有技术的不足，提供了一种基于波场延拓原理的起伏地表条件下组 合震源激发定向波场的方法。
[0005] 为了实现本发明的上述目的，基于波场延拓原理的起伏地表组合震源定向波场方 法包括以下步骤：
[0006] a、建立起伏地表地球物理参数模型，将模型参数离散化至每个网格节点；
[0007] b、记录此模型的地表水平位置与相应高程坐标，沿起伏地表设置N个水平间隔一 致的点震源Si, S2，…，Sn，分别记录其尚程坐标为ZSl，ZS2,…ZSn;
[0008] C、将震源的高程进行排序，取其中的最大高程，作为水平基准面高程ZH = maX(Zsi， zs2,…zsn)，max是取最大值函数；
[0009] d、从模型最左网格至最右网格，将起伏地表模型向上填充等效介质至水平基准 面，等效介质物理参数为起伏地表表层位置处的地球物理参数；
[0010] e、将起伏地表的N个地震点震源投影至在水平基准面上，记录起伏地表至水平基 准面尚程 Δ ZS1 = Z『ZS1，Δ ZS2 = Z『ZS2，…，Δ ZSn = Z『ZSn;
[0011] f、根据水平地表组合震源的方向因子公式计算相邻震源的延时时差参数τ;
[0012] g、按照延时参数依次激发N个点震源，得到N个时间域波场，分别为Ul(Xl，zH，t)，u 2 (χ2,ΖΗ,τ+?) ,··· ,U2(X2,ZH, (η-1)*τ+?);
[0013] h、将此N个时间域波场分别对时间和水平坐标两个维度进行傅里叶变换，得到N个 频率波数域波场，分别为U1 (kx，Zh，ω )，U2 (kx，Zh，ω )，…，Un(kx，Zh，ω )，其中kx是水平波数；
[0014] i、根据频率波数域的下行波动方程，将N个频率波数域波场以自激自收方式正向 外推（向下延拓)至实际起伏地表高程，结果分别为U 1 (kx，ZS1，ω )，U2(kx，ZS2，ω )，…，Un(kx， Z Sn , ^ )
[0015] j、将延拓后N个频率波数域波场叠加，得到频率域起伏地表组合震源定向波场记 录；
[0016] k、将上步结果进行时间和水平坐标的二维反傅里叶变换，得到时间域定向传播的 波场。
[0017] 有益效果:经模拟试验，本发明公开的起伏地表组合震源波场定向的方法，计算准 确，有效得到定向激发的地震波场，并改善了由于地表非均匀性引起的波场计算误差，在实 际应用中更具有普遍性。从图中明显看出，本发明通过a~k步骤能够得到任意方向角的地 震波场，且起伏地表对定向传播的地震波场的影响很小。
【附图说明】
[0018] 图1复杂起伏地表地质模型图
[0019] 图2组合震源定向波场在0.3秒处的波场快照图。
[0020] 2a是方向角Θ为30度，2b是方向角Θ为〇度，2c是方向角Θ为-30度
【具体实施方式】
[0021]下面结合附图和实例对本发明进一步的详细说明。
[0022] -种起伏地表条件下组合震源激发定向波场的方法，在图1所示的非均匀起伏地 表模型情况下，包括如下步骤：
[0023] a、建立起伏地表地球物理参数模型，将模型参数离散化至每个网格节点；
[0024] b、记录此模型的地表水平位置与相应高程坐标，沿起伏地表设置N个水平间隔一 致的点震源Si, S2，…，Sn，分别记录其尚程坐标为ZSl，ZS2,…ZSn，设震源间距为d;
[0025] c、将震源的高程进行排序，取其中的最大高程，作为水平基准面高程ZH = max(Zsl， zs2,…zsn)，max是取最大值函数；
[0026] d、从模型最左网格至最右网格，将起伏地表模型向上填充等效介质至水平基准 面，等效介质物理参数为起伏地表表层位置处的地球物理参数；
[0027] e、将起伏地表的N个地震点震源投影至在水平基准面上，记录起伏地表至水平基 准面尚程 Δ ZSl = ZH-ZSl，Δ ZS2 = ZH_ZS2，…，Δ ZSn = ZH-ZSn;
[0028] f、设计地震波传播的方向角为Θ，根据水平地表组合震源的方向因子公式计算相 邻震源的相位延迟β，根据相位延迟与延时关系得到震源的延时时差参数τ = β/ω，其中ω 是角频率；
[0029] g、按照延时延时依次激发N个点震源，得到N个时间域波场，分别为m(X1，ZH，t)，u 2 (X2，ZH，T+t)，···，U2(X2，ZH, (n_l)*T+t);
[0030] h、将此N个时间域波场分别对时间和水平坐标两个维度进行傅里叶变换，得到N个 频率波数域波场，分别为U1 (kx，ZH，ω )，U2 (kx，ZH，ω )，…，Un(kx，ZH，ω )，其中kx是水平波数； [0031 ] i、根据频率波数域的下行波动方程，将N个频率波数域波场以自激自收方式正向 外推（向下延拓)至实际起伏地表高程，结果分别为U1 (kx，ZS1，ω )，U2(kx，ZS2，ω )，…，Un(kx， ZSn, ω )，其中波场延拓公式为
[0032] U(kx，zsi，ω ) = -U(kx，ZH, ω )exp(ikz Δ Zsi)，
[0033] U(kx，zS2, ω ) = -U(kx，ZH, ω )exp(ikz Δ zS2)，
[0034] …，
[0035] U(kx，zsn，ω ) =-U(kx，ZH，ω )exp(ikz Δ zsn) ·
[0036] 其中kz是垂直波数，符合频散关系 ，V代表地下介质速度，是空间变 量的函数v(x，z);
[0037] j、将延拓后N个频率波数域波场叠加，得到频率域起伏地表组合震源定向波场记 录；
[0038] k、将上步结果进行时间和水平坐标的二维反傅里叶变换，得到时间域定向传播的 波场。
【主权项】
1. 一种基于波场延拓原理的起伏地表组合震源波场定向方法，其特征在于，包括以下 步骤： a、 建立起伏地表地球物理参数模型，将模型参数离散化至每个网格节点； b、 记录此模型的地表水平位置与相应高程坐标，沿起伏地表设置N个水平间隔一致的 点震源Si，S2，…，Sn，分别记录其尚程坐标为ZS1，ZS2，'"ZSi!; C、将震源的高程进行排序，取其中的最大高程，作为水平基准面高程ZH = max(Zsl， zs2,…zsn)，max是取最大值函数； d、 从模型最左网格至最右网格，将起伏地表模型向上填充等效介质至水平基准面，等 效介质物理参数为起伏地表表层位置处的地球物理参数； e、 将起伏地表的N个地震点震源投影至在水平基准面上，记录起伏地表至水平基准面 尚程 Δ ZS1 = ZH-ZS1，Δ ZS2 = ZH-ZS2，…，Δ ZSn = ZH-ZSn; f、 根据水平地表组合震源的方向因子公式计算相邻震源的延时时差参数τ; g、 按照延时参数依次激发Ν个点震源，得至ijN个时间域波场，分别为m(X1，ZH，t)， U2(X2， ΖΗ,τ+t)，···，U2(X2，zh，（n_l)*T+t); h、 将此N个时间域波场分别对时间和水平坐标两个维度进行傅里叶变换，得到N个频率 波数域波场，分别为Ui(kx，ZH, ω )，U2(kx，ZH, ω )，'"，Un(kx，ZH, ω )，其中kx是水平波数； i、 根据频率波数域的下行波动方程，将N个频率波数域波场以自激自收方式正向外推 至实际起伏地表高程，结果分别为iMk^zsi，ω )，U2(kx，ZS2, ω )，···，Un(kx，ZSn，ω ) j、 将延拓后Ν个频率波数域波场叠加，得到频率域起伏地表组合震源定向波场记录； k、 将j步骤结果进行时间和水平坐标的二维反傅里叶变换，得到时间域定向传播的波 场。
【文档编号】G01V1/28GK105842731SQ201610308447
【公开日】2016年8月10日
【申请日】2016年5月11日
【发明人】巩向博, 王升超, 于晨霞, 韩立国
【申请人】吉林大学