《反射波广义拉东谱法绕射波场分离技术.pdf》由会员分享，可在线阅读，更多相关《反射波广义拉东谱法绕射波场分离技术.pdf（8页珍藏版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 102854529 A (43)申请公布日 2013.01.02 CN 102854529 A *CN102854529A* (21)申请号 201210241784.8 (22)申请日 2012.07.13 G01V 1/28(2006.01) (71)申请人 孙赞东 地址 北京市昌平区府学路 18 号中国石油大 学 ( 北京 ) 申请人 白英哲 (72)发明人 孙赞东 白英哲 (74)专利代理机构 北京国林贸知识产权代理有 限公司 11001 代理人 李桂玲 李富华 (54) 发明名称 反射波广义拉东谱法绕射波场分离技术 (57) 摘要 地震勘探通过地震波对地下介。
2、质的响应来推 知地下信息， 当地下存在介质突变时， 往往响应为 绕射波。在塔里木碳酸盐岩地区， 主要的油气储 集空间是次生的溶蚀空洞和裂缝， 这些孔洞和裂 缝的地震响应以绕射为主。因此有必要研究这些 绕射波， 但绕射波与反射波特征近似， 往往难以识 别。 因此有人提出了通过倾角域偏移成像的方法， 在倾角域实现绕射波场的分离。我们沿着这一思 路， 通过绕射波和反射波在倾角域的不同响应， 提 出一种新的倾角域绕射波场分离方法反射波 广义拉东谱的绕射波场分离方法， 成功地在倾角 道集中将绕射波和反射波进行分离， 并对分离后 的绕射波场单独成像。 并通过实际资料上的应用， 验证了方法的有效性和实用性。。
3、 (51)Int.Cl. 权利要求书 1 页 说明书 4 页 附图 2 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书 1 页 说明书 4 页 附图 2 页 1/1 页 2 1. 本发明意在通过提供一种新的倾角域绕射波场分离技术， 通过反演的方法预测反射 波， 在预测反射波的过程中通过拾取反射波广义拉东谱对反射波方程进行线性化， 从而求 取倾角道集中的反射波， 以便更好的分区倾角道集中的绕射波和反射波。 本专利使用倾角域偏移得到的倾角域共反射点道集作为输入， 实现上述目的采取的技 术方案如下 ： 步骤 1 ： 对输入的倾角域共反射角道集按不同倾角不同深度的反射曲线进。
4、行叠加， 将 叠加的数据按叠加使用的倾角和深度排布， 形成一张谱， 称作反射波广义拉东谱。 步骤 2 ： 拾取步骤 1 得到的反射波广义拉东谱上的能量团， 获得反射界面的深度和倾 角。 步骤 3 ： 通过褶积理论和由步骤 2 求得的反射界面的深度和倾角预测反射波。 步骤 4 ： 将步骤 3 预测的反射波从总场中减去， 就得到绕射波。 2. 根据权利要求书 1 所述的反射波广义拉东谱法绕射波场分离， 其特征在于， 所述的 步骤 1 中的生成反射波广义拉东谱的方法为 ： 将数据按反射波规律进行叠加， 式中 S 为反射波广义 Radon 谱 A 为界面的视倾角， Z 为界面的是深度， 倾角道集的倾 。
5、角维度， z 为倾角道集的深度维度。 反射波广义 Radon 变换是沿着不同深度， 不同倾角的反射曲线叠加， 沿这些路径叠加 的数据在反射波广义 Radon 谱中表示为一个点。当叠加路径与反射同相轴重合时， 叠加能 量最强， 因此反射同相轴在反射波广义 Radon 谱上会收敛为一个能量团， 能量团的位置指 示着反射界面的深度和倾角。 通过自动拾取反射波能量团就可以获得反射界面的视深度和 视倾角。 3. 根据权利要求书 1 所述的反射波广义拉东谱法绕射波场分离， 其特征在于， 所述的 步骤 3 中预测反射波使用下述方程 ： 式中 R 为反射波场， w 为地震子波， 是狄拉克函数， ri、 Zi和。
6、 Ai分别为第 i 个界面的 反射系数、 深度和倾角， 为观测倾角， z 为观测深度。 认为在全波场中反射波是主要的， 这假设在通常的地震记录中都能满足， 则可以写出 误差函数 E |Q-R|2 式中 Q 为观测数据， 即倾角道集， E 为观测误差。式中的 Zi和 Ai分别通过反射波广义 拉东谱的方法求取， ri通过广义最小二乘法求取。求取后可以重构反射波场 R。 权 利 要 求 书 CN 102854529 A 2 1/4 页 3 反射波广义拉东谱法绕射波场分离技术 技术领域 0001 本发明涉及油气田勘探技术领域， 属于地震资料处理范畴， 具体的说是在倾角域 实现绕射波与反射波分离的技术。。
7、 背景技术 0002 对于复杂介质， 如塔里木盆地， 油气储集在古老的寒武系和奥陶系碳酸盐岩地层 中， 这些地层的埋深在 5000 米以下， 原生孔隙早已被压实、 充填殆尽， 如今的储集空间是靠 溶蚀、 断裂等作用后天形成的次生孔隙。 这些孔隙的分布不受沉积相带的控制， 且分布的规 模及其的不均匀。因此， 想要寻找这些洞缝型储层只能通过地震的手段。 0003 溶蚀孔洞的尺寸往往较之地震波的波长要小得多， 可以认为是岩性的突变点。在 岩性的突变点上地震波就会发生绕射， 只要探测到这些绕射波， 通过地震偏移技术， 绕射波 收敛， 就可以确定溶蚀孔洞的位置。然而事实并不这么简单， 在偏移后的剖面上，。
8、 绕射波往 往被收敛成一个点， 从而淹没在众多的反射同相轴之间， 不易分辨。 在通常的叠前偏移后的 共反射点道集 (CRP) 中， 绕射点和反射点又同样表现为同相轴的拉平， 因此， 通过叠前偏移 的共反射点道集 (CRP) 也不能区分反射波与绕射波。 0004 目前国内外在绕射波场分离方面的研究还处于起步阶段， 还没有成型的研究成 果。在勘探现场， 绕射波仍是处理的难题。在很多情况下， 由于没有有效的分离手段， 只能 将绕射波当成一种干扰进行压制。 但这样一来， 地下的介质突变信息就丢失了， 地震勘探的 分辨率也随之下降。因此， 石油勘探行业尤其是在面对具有溶蚀孔洞的复杂储层时迫切需 要一种行。
9、之有效的绕射波场分离方法。 0005 因此国外有学者提出应该将绕射波从反射波中分离出来， 单独成像研究。传统的 绕射波场分离只是在叠后剖面上分离 “串珠状” 反射， 精度和准度受到很大制约。 Landa等人 于 1998 年提出了 D-section 的概念， 提出用绕射波来检测局部非均质性。Vermeulen 等人 于 2006 提出通过相干异常的方法来分离断点的绕射波。于此同时， Taner 等人提出在平面 波震源下绕射波场与平面波场的差别， 并主张用拉东变换分离绕射波。Landa 等人于 2008 年正式提出了倾角域绕射波场分离的概念。在此基础上白英哲等人又做了大量的研究工 作， 推导了。
10、不同维数下绕射波与反射波的倾角域响应， 将倾角域成像推广到三维。 在此同时 Klokov 等人提出了一套新的类拉东变换的倾角域绕射波场分离方法， 并开始了对断棱绕射 的讨论。 发明内容 0006 本发明意在通过提供一种新的倾角域绕射波场分离技术， 通过反演的方法预测反 射波， 在预测反射波的过程中通过拾取反射波广义拉东谱对反射波方程进行线性化， 从而 求取倾角道集中的反射波， 以便更好的分区倾角道集中的绕射波和反射波。 0007 本专利使用倾角域偏移得到的倾角域共反射点道集作为输入， 实现上述目的采取 的技术方案如下 ： 说 明 书 CN 102854529 A 3 2/4 页 4 0008 。
11、步骤 1 ： 对输入的倾角域共反射角道集按不同倾角不同深度的反射曲线进行叠 加， 将叠加的数据按叠加使用的倾角和深度排布， 形成一张谱， 称作反射波广义拉东谱。 0009 将数据按反射波规律进行叠加， 0010 0011 式中 S 为反射波广义拉东谱 A 为界面的视倾角， Z 为界面的是深度， 倾角道集的 倾角维度， z 为倾角道集的深度维度。 0012 反射波广义拉东变换是沿着不同深度， 不同倾角的反射曲线叠加， 沿这些路径叠 加的数据在反射波广义拉东谱中表示为一个点。当叠加路径与反射同相轴重合时， 叠加能 量最强， 因此反射同相轴在反射波广义拉东谱上会收敛为一个能量团， 能量团的位置指示 。
12、着反射界面的深度和倾角。 通过自动拾取反射波能量团就可以获得反射界面的视深度和视 倾角。 0013 步骤 2 ： 拾取步骤 1 得到的反射波广义拉东谱上的能量团， 获得反射界面的深度和 倾角。 0014 步骤 3 ： 通过褶积理论和由步骤 2 求得的反射界面的深度和倾角预测反射波。 0015 预测反射波使用下述方程 ： 0016 0017 式中 R 为反射波场， w 为地震子波， 是狄拉克函数， ri、 Zi和 Ai分别为第 i 个界 面的反射系数、 深度和倾角， 为观测倾角， z 为观测深度。 0018 认为在全波场中反射波是主要的， 这假设在通常的地震记录中都能满足， 则可以 写出误差函数。
13、 0019 E |Q-R|2 0020 式中 Q 为观测数据， 即倾角道集， E 为观测误差。式中的 Zi和 Ai分别通过反射波 广义拉东谱的方法求取， ri通过广义最小二乘法求取。求取后可以重构反射波场 R。 0021 步骤 4 ： 将步骤 3 预测的反射波从总场中减去， 就得到绕射波。 0022 令反射波为 D 则可以通过下式表示绕射波的分离过程， 0023 D Q-R 0024 以上具体实施方式仅用于说明本发明， 而非用于限定本发明。 附图说明 0025 图 1(a) 是一个合成倾角道集， 由一个向上开口的反射同相轴和一个水平的绕射 同相轴构成 ； 0026 图 1(b) 是图 1(a)。
14、 所示倾角道集的反射波广义拉东谱 ； 0027 图 1(c) 是图 1(a) 所示倾角道集由反射波广义拉东谱法重构的反射波道集 ； 0028 图 1(d) 是图 1(a) 所示倾角道集中总场减去反射波得到的绕射波道集。 0029 图 2 是中值滤波结果和反射波广义拉东谱法分离的绕射波在 30 度倾角处地震道 的比较。 0030 图 3(a) 是实际数据中的一个倾角道集 ； 说 明 书 CN 102854529 A 4 3/4 页 5 0031 图 3(b) 是图 3(a) 所示倾角道集的反射波广义拉东谱 ； 0032 图 3(c) 是图 3(a) 所示倾角道集中以反射波广义拉东谱法分离的绕射波。
15、场。 具体实施方式 0033 通过一个合成倾角域共反射点道集来说明 ： 0034 图 1(a) 是一个合成的倾角道集， 图 1(b) 为该倾角道集的反射波广义拉东谱。按 照步骤 1 所述， 反射波广义拉东变换是沿着不同深度， 不同倾角的反射曲线叠加， 道集上的 几条虚实不同的曲线是几个代表的叠加路径， 沿这些路径叠加的数据在反射波广义拉东谱 中表示为一个点， 在谱上十字表示的点对应着相应颜色的路径。当叠加路径与反射同相轴 重合时， 叠加能量最强， 因此反射同相轴在反射波广义拉东谱上会收敛为一个能量团， 能量 团的位置指示着反射界面的深度和倾角。按步骤 2 所述， 通过自动拾取反射波能量团就可 。
16、以获得反射界面的深度 Z 和倾角 A。 0035 由步骤 3 可知， 当通过反射波广义拉东谱求取了各个反射界面的深度和倾角后， 反演反射波的问题就简化成了一个线性反演的问题， 未知参数就只剩下的各个界面的反射 系数 r， 通过线性的最小二乘法可以轻松的求出各层的反射系数， 进而重构反射波 ( 如图 1(c) 所示 )。最后按照步骤 4， 在全波场中减去反射波就可以得到绕射波场， 如图 1(d) 所 示。 0036 按照上述方法得到的结果可以认为没有振幅上的损失， 这比传统的中值滤波法有 明显的进步。图 2 是中值滤波结果和反演分离的绕射波在 30 度倾角处地震道的比较， 图中 的实线是原始记录。
17、， 虚线是中值滤波的记录， 点是反演得到的绕射波记录。 通过对比可以看 处， 在水平绕射同相轴的情况下， 中值滤波和反演法都能很好的保持绕射振幅， 但当绕射点 非正上方观测， 绕射同相轴的斜率有变化时， 反演法可以很好的恢复绕射同相轴， 其振幅与 原始记录非常接近， 但中值滤波会有较大偏差。 0037 参考文献 0038 Landa E， Keydar S.Seismic monitoring of diffraction images for detection of local heterogeneities.Geophysics， 1998， 63 ： 1093-1100. 0039 V。
18、ermeulen J， Gurevich B， Urosevic M， et al.Enhancing Coherency analysis for fault detection and mapping using 3D diffraction imaging.76th SEG Annual Meeting， 2006， 1108-1112. 0040 Taner M T， Fomel S， Landa E.Prestack separation of seismic diffractions using plane-wave decomposition.76th Annual Meetin。
19、g.SEG Annual Meeting， 2006， 2401-2404. 0041 Landa E.， Fomel S， Reshef M.Separation imaging and velocity analysis of seismic diffractions using migrated dip-angle gathers.78th SEG Annual Meeting， 2008， 2176-2180. 0042 Bai Y Z， Sun Z D， Yang H J， et al.Seismic diffraction separation in 2D and 3D space。
20、s.73rd EAGE annual meeting， 2011. 0043 Klokov A， Baina R， Landa E.Separation and imaging of seismic diffractions in dip angle domain.72nd EAGE annual meeting， 2010. 说 明 书 CN 102854529 A 5 4/4 页 6 0044 Klokov A， Baina R， Landa E.Point and edge diffractions in three dimensions.73rd EAGE annual meeting.2011。 说 明 书 CN 102854529 A 6 1/2 页 7 图 1 说 明 书 附 图 CN 102854529 A 7 2/2 页 8 图 2 图 3 说 明 书 附 图 CN 102854529 A 8 。