首页 >
论文范文 >
工程工业论文 >
工程综合论文 >
矿业工程论文 >
频率—波数域水陆检数据上、下行波场分离方法
摘要：随着海洋油气勘探开发的发展，OBC/OBN水陆检海底接收的应用日趋广泛，上、下行波场分离作为水陆检数据处理的关键技术，决定了资料处理品质及应用效果。常规方法分离出来的上行波场中包含下行波场，下行波场中包含上行波场，即不能实现上行波场与下行波场的完全分离。为此，提出在频率—波数域利用包含直达波和切除直达波后的水陆检数据，直接计算水陆检数据标定滤波因子，实现水陆检数据标定和上、下行波场分离。分离后的上行波场数据中既消除了虚反射多次波干扰，又提高了地震数据信噪比和分辨率，为后续联合反褶积和偏移成像等处理提供了高保真的上行波场和下行波场数据。数据实例表明了所提方法的有效性和实用性。
加入收藏
拖缆是海洋地震勘探的主要采集方式，因海平面这个强反射界面产生的虚反射是拖缆数据中最严重干扰之一，造成地震数据频谱出现严重的陷波畸变，常规的数学处理方法难以恢复被陷波畸变了的频率成分。海底节点（OBN）采集技术将特性不同的水、陆两种检波器同时放置海底进行记录。水检记录的是压力波场，陆检垂直分量记录的是垂直速度波场，两种波场反射信号极性相同，虚反射和鸣震干扰极性相反。通过适当的数据合并与组合，物理上可有效消除虚反射和海水鸣震的影响。其关键是水检和陆检两种数据在振幅、相位等方面要尽量保持一致[1-6]。
上、下行波场分离技术是水陆检数据处理的基础和桥梁。高少武等[7-25]开展了海洋水陆检数据处理方法和技术的系统研究，首次提出了基于水检数据平均均方根振幅、陆检数据平均均方根振幅、标定因子、海水深度、海底反射系数、海平面反射系数、海水速度七个参数的上行波场分离公式，创立了七参数水陆检数据合并处理方法与技术[7-10]；首次提出了基于水检数据上、下行波场贡献因子和陆检数据上、下行波场贡献因子四个波场分离贡献因子参数的上、下行波场分离公式，创立了四参数水陆检数据上、下行波场分离方法与技术[20-23]；创造性提出了基于相位特征方程的相位匹配方法、匹配滤波器算子的滤波器匹配方法、最大方差模方程分子分母系数特征方程的标定因子匹配方法三大技术系列，组成了水陆检数据相位匹配、滤波器匹配、均方根振幅匹配和标定因子匹配的水陆检数据匹配方法与技术[10-13]。这些方法都是基于常数标定，不能满足水、陆检数据处理需要。
Amundsen[26]最早提出了频率—波数域水陆检数据上、下行波场分离计算公式。随后，Day等[27]、Amundsen等[28]和Caprioli等[29]相继发展了一系列频率—波数域水陆检数据标定和上、下行波场分离理论方法和技术，但都是使用直达波之下的水陆检数据进行标定因子计算。由于直达波之下水陆检数据信噪比低，特别是陆检数据中包含着较强的噪声，计算的标定因子存在误差，造成上、下行波场分离不彻底，上行波场中包含少量下行波场，下行波场中包含少量上行波场，无法满足水陆检数据处理需要[30-32]。
水陆检数据直达波部分信噪比高，为此，充分利用直达波部分数据，在频率—波数域直接计算水陆检数据标定因子，然后进行水陆检数据标定和上、下行波场处理，可实现频率—波数域水陆检数据上、下行波场分离，为双检数据的后续处理奠定良好的基础[33]。
1、方法原理
根据文献[26]，标定陆检数据之后，频率—波数域海底之上的上、下行波场数据可以表示为
式中：H (k,m)为水检数据；G (k,m)为陆检数据；Cz(k)为频率域标定滤波因子，k=1,2,⋯,NT,NT是频率采样点总个数；U (k,m)为标定处理后频率—波数域海底之上的上行波场；D (k,m)是标定处理后频率—波数域海底之上的下行波场，m=1,2,⋯,NX,NX是波数采样点总个数；ωk是第k个频率样点处的圆频率；kxm是第m个波数样点处空间方向圆波数；ρw为海水密度；cw是海水速度。
使用切除函数切除水陆检数据中直达波场数据。切除直达波场数据之后频率—波数域海底之上水陆检数据为
式中：Wd(k,m)为频率—波数域切除函数，用于从频率—波数域水检数据和陆检数据中去除直达波场数据；Hr(k,m)是切除直达波场之后频率—波数域水检数据；Gr(k,m)是切除直达波之后频率—波数域陆检数据。
切除直达波之后频率—波数域下行波场为
建立下列目标函数
式（7）两边对频率域标定滤波因子Cz[k]求导，并令导数为零，简化整理得
其中
式中上标“*”表示复共轭运算。
为了增加式（8）的稳定性，式（8）右边分子和分母乘以共轭因子Bz*[k]，则频率域标定滤波因子计算公式为
式中ε是一个非常小的正数（例如0.00000001），防止被零除产生的不稳定和数值溢出。详细推导过程见文献[30]。
2、模型数据处理
图1是水平层状深度—速度模型，炮点沉放深度为7.5 m，炮点间距为25 m；检波点沉放海底，检波点间距为25 m。图1模型合成水陆检共检波点道集数据如图2a、图2b。对比合成数据不同方法的处理结果可以看出：对于上行波场，常规方法（图2c）和本文方法（图2e）分离效果基本相当，但是常规方法结果中还包含少量下行波场（图2c中红色箭头所示）；常规方法分离的下行波场中（图2d）仍包含着较强的上行波场（图中蓝色和黄色箭头所示），而对应本文方法结果（图2f）中几乎看不到上行波场残留。说明本文方法处理后，上行波场和下行波场完全分离开来，特别是下行波场，几乎看不到上行波场残余。
图1 水平层状深度—速度模型
图2 图1模型合成水陆检数据及不同方法处理后的结果
图3是SEAM速度—深度模型，检波点位于海底1500 m处，检波点间距为10 m，分布范围如图中绿色方框所示；炮点位于水深15 m处，炮点间距为10 m，分布范围如图中红色方框所示。数据是SEG发布的标准海洋海底节点四分量模拟数据，水、陆检共炮点道集数据分别如图4a和图4b。本文方法与主流商业软件处理SEAM数据的结果对比显示：对于远炮检距数据，由于处理是在共检波点道集上进行，当共检波点道集数据道数较少时，主流商业软件不进行处理，直接将水检数据输出到上行波场，陆检数据输出到下行波场中，而本文方法使用相邻共检波点道集算子进行标定并进行上、下行波场分离处理（如图4红色箭头所示）；对于近炮检距数据，主流商业软件方法与本文方法分离效果总体相当，仅局部稍有差别；对于上行波场，主流商业软件方法（图4c）中还包含少量下行波场，而对应本文方法结果（图4e）中几乎看不到下行波场残留（图中蓝色箭头所示）；对于下行波场，主流商业软件方法（图4d）分离的下行波场中还包含着少量的上行波场，而对应本方法处理结果（图4f）中几乎看不到上行波场残留（图中黄色箭头所示）。说明本文方法处理后，上行波场和下行波场完全分离开来，不含剩余上、下行波场。
图3 SEAM速度—深度模型
图5是实际水陆检共炮点道集数据本文方法与主流商业软件方法处理结果对比，图6是实际水陆检叠后数据本文方法与主流商业软件方法处理结果对比。两种方法分离效果总体相当，但本文方法局部细节优于主流商业软件方法，如图6中圆圈区所示。
图4 SEAM模型水陆检共炮点道集数据及不同方法处理结果对比
图5 实际水陆检共炮点道集数据及不同方法处理结果对比
图6 实际水陆检叠后数据及不同方法处理结果对比
3、结束语
本文提出的频率—波数域水陆检数据上、下行波场分离方法，通过计算频率—波数域陆检数据加权自相关、频率—波数域水陆检数据加权互相关、切除直达波后频率—波数域陆检数据加权自相关、切除直达波后频率—波数域水陆检数据加权互相关，直接计算频率域水陆检数据标定滤波因子，实现水陆检数据上行波场与下行波场的分离处理，消除地震数据中海水多次波干扰影响，提高地震数据信噪比和分辨率，为后续联合反褶积和偏移成像处理提供高保真的上行波场与下行波场数据。模型数据和实际数据处理结果表明：
(1）本文方法直接计算频率域标定滤波因子，可实现水检数据上、下行波场的波场能量匹配和波场分离，是水陆检数据上、下行波场分离的有效方法；
(2）本文方法的水陆检数据上、下行波场的波场分离处理结果，上行波场有效地压制了虚反射多次波干扰，提高了反射数据的信噪比和分辨率；
(3）使用频率—波数域水陆检数据和切除直达波后频率—波数域水陆检数据，直接计算频率域标定滤波因子，使得频率域标定滤波因子参数的计算更加稳定和准确。
参考文献:
[7]高少武,钱忠平,马玉宁,等. OBC水陆检数据合并处理技术[J].石油地球物理勘探, 2018, 53(4):703-709.
[8]高少武,赵波,罗国安,等.一种水陆检波器地震数据合并处理的方法及装置:ZL201510556940. 3[P]. 2015-09-02.
[9]高少武,赵波,钱忠平,等.一种水陆检波器地震数据合并处理的方法及装置:ZL201710565184. X[P]. 2017-07-12.
[10]高少武,祝宽海.一种水陆检波器地震数据匹配方法:ZL201310531139. 4[P]. 2013-10-31.
[11]高少武,赵波,高轩,等. OBC水陆检数据匹配技术[J].石油地球物理勘探, 2015, 50(1):29-32.
[12]高少武,黄燕.一种水陆检波器地震数据标定因子反演方法:ZL201310736942. 1[P]. 2013-11-29.
[13]高轩,弓小平,高少武,等. OBC水陆检数据标定因子估算方法[J].石油地球物理勘探, 2016, 51(1):49-53.
[14]高少武,赵波,罗国安,等.一种水陆检波器数据海底反射系数反演方法:ZL201410379202. 1[P]. 2013-11-29.
[15]高少武,钱忠平,马玉宁,等.基于延迟波场特征法的海底反射系数估算技术[J].石油物探, 2018, 57(3):389-394.
[16]高少武,赵波,钱忠平,等.水陆检波器数据海水深度反演方法和装置:ZL201410377324. 7[P]. 2014-08-01.
[17]高少武,赵波,钱忠平,等.一种确定海水深度的方法和装置:ZL201510963991. 8[P]. 2015-12-21.
[18]高少武,钱忠平,孙鹏远,等.基于数据分段互双谱的海水深度计算方法[J].石油地球物理勘探, 2019,54(4):751-756.
[19]高少武,钱忠平,黄少卿,等. OBC水陆检数据海水深度估算方法[C].中国石油学会2017年物探技术研讨会论文集, 2017, 996-1000.
[20]高少武,钱忠平,孙鹏远,等.水陆检数据上下行波场分离方法[J].石油地球物理勘探, 2020, 55(5):991-996.
[21]高少武,钱忠平,李建峰,等.水陆检数据上下行波场分离[C].中国石油学会2019年物探技术研讨会论文集, 2019, 1086-1089.
[22]高少武,张少华,钱忠平,等.一种水陆检波器地震数据上下行波场分离方法及装置:ZL201810101188. 7[P]. 2018-02-01.
[23]高少武,张少华,王成祥,等.海底节点地震数据上下行波场分离的方法及装置:ZL201910144556. 0[P]. 2019-02-27.
[24]高少武.去除海洋拖缆数据接收点虚反射方法[C]. 2018勘探地球物理学进展学术研讨会论文集,2018, 253-263.
[25]高少武,孙鹏远,方云峰,等.双检数据上下行波场分离技术研究进展[J].石油地球物理勘探, 2021, 56(6):1419-1429.
[30]李映艳,张乐乐,陈刚,等.基于P/S波解耦和上下行波分离的二维VTI介质弹性波逆时偏移[J].石油地球物理勘探, 2023, 58(3):651-659.
[31]秦宁,王常波,梁鸿贤,等.一次波和层间多次波联合成像方法[J].石油地球物理勘探勘探, 2022, 57(6):1375-1383.
[32]青杰,顾汉明,王建花.频率—慢度域延拓关键参数优化的鬼波压制方法[J].石油地球物理勘探勘探,2023, 58(1):121-132.
[33]高少武,方云峰,马光凯,等.一种水陆检数据标定的方法:202111667515. 3[P]. 2021-12-30.
基金资助:中国石油集团公司项目“地震处理解释关键新技术与智能化软件研发”(2021ZG03)和“海洋节点(OBN)地震数据处理软件开发”(2021ZG02)联合资助;
文章来源:高少武,于万辉,张建磊,等.频率—波数域水陆检数据上､下行波场分离方法[J].石油地球物理勘探,2024,59(04):755-762.
分享：
围岩松动圈是指地下工程开挖后,由于应力重分布和爆破等因素影响,在洞室周围形成的具有一定厚度的破裂､松动区域[1]｡这一区域内的岩体力学性质发生显著变化,主要表现为强度降低､渗透性增加等特征｡围岩松动圈的形成和发展直接影响地下工程的稳定性和安全性｡
2025-09-07
钻孔探测作为矿井地质分析预测及灾害处理的重要手段,具有实揭性､灵活性及长距离等优点,常用于煤矿瓦斯抽放､探放水及地质填图等勘探工程｡然而,受限于软碎岩层､断层发育带等井下复杂地质条件,常规钻孔极易产生塌孔､缩径的问题,严重影响钻孔的成孔效率及施工效果｡
2025-08-03
断层是地壳中常见的地质构造之一,是由破断面两层岩体发生明显位移而形成的｡在采煤过程中,断层的发育增加了开采的难度,影响采煤方法的选择和回采速度｡断层的大小､走向､数量,均会加剧煤层漏风问题,进而增加煤层自燃的风险,对此,诸多学者开展了大量研究[1-2]｡
2025-08-03
贵州省具有丰富的煤炭资源,素有“西南煤海”之称,是我国华南型煤田中最大的产煤省区｡贵州省喀斯特矿区在高强度开采扰动后,极易造成矿井突水事故,给井下开采工作带来极大的安全隐患和经济损失[1]｡导水裂缝带,即垮落带和裂缝带,是贯通采空区与上部水体的通道,覆岩“两带”的高度及特征对于矿井井下水害防治具有重要的意义[2-3]｡
2025-08-03
采煤机自动拖缆是指采煤机在工作运行过程中,采用特定的拖曳装置移动电缆,替代人工来管理电缆,避免采煤机往复运行中,电缆因多次叠加引起电缆出槽､损伤等现象｡采煤机自动拖缆装置对于采煤机的正常运行至关重要,特别是对电缆槽高度较低的薄煤层工作面,效果更为显著,它确保了电缆不被损坏,从而保证了采煤机的电力供应｡
2025-08-03
光纤技术是基于光传播原理的高新技术创新,广泛应用于通信､传感､医疗等领域,光纤传感器利用光波在光纤内部的传播特性,这些特性使得光纤传感器在复杂环境中,特别是需要长期监测的工程领域中,应用更加广泛[1]｡在采矿作业中,随着深部采矿技术的不断发展,围岩的复杂性和变异性不断增加,对支护结构的稳定性提出了新的挑战[2]｡
2025-08-03
陶文斌[8]研究了大断面过破碎带失稳变形原因,提出了“预注浆+超前支护+架棚+二次注浆”的支护方案,有效保证了巷道围岩的稳定;伊丽娟[9]以某煤矿大巷过破碎带为工程背景,提出“锚杆+锚索+U型钢”的联合支护技术,成功控制了巷道围岩变形问题｡针对不同地质环境应确定相应的巷道支护技术｡
2025-08-03
煤炭企业的环境会计信息披露通常在社会责任报告中设专章,内容涉及环保管理体系､环保投资､污染物排放､重大环境问题及整改情况､节能减排等｡随着“双碳”目标的提出,许多煤炭企业深刻认识到低碳的重要性,开始在年度财务报告中披露环境会计信息｡
2025-07-12
针对荫营煤矿密闭区域检测环节操作繁琐，检测结果反馈时间长、效率低等问题，归纳了密闭“呼吸”现象的影响因素，基于现场测试分析了密闭“呼吸”现象的基本特征，对煤自燃气体产生的规律进行研究，结果表明：大气参数的变化导致密闭区出现“呼吸”效应，密闭内外压差的变化范围为-160～350 Pa,说明荫营煤矿密闭存在“呼吸”效应，且密闭内外气体交换量较大。研究结果对荫营煤矿快捷高效地实现井下密闭火情智能监测预警及风险管控以预防煤自燃具有重要指导和现实意义。
2025-07-07
煤炭服务企业与新一代信息技术的深度融合势在必行｡许多企业正在有序构建及优化管理､业务及服务三大核心系统,以实现信息化､数字化与智能化的融合升级｡管理层面涵盖企业资源规划(ERP)系统､办公自动化(OA)系统以及管理信息系统(MIS)等模块;业务层面涉及数字化协同设计系统､数字化设计交付系统以及工程采购与施工(EPC)总承包管理系统等。
2025-07-07我要评论
期刊名称：矿业工程研究
期刊人气：868
主管单位：湖南省教育厅
主办单位：湖南科技大学
出版地方：湖南
专业分类：煤矿
国际刊号：1674-5876
国内刊号：43-1493/TD
邮发代号：42-359
创刊时间：1980年
发行周期：季刊
期刊开本：16开
见刊时间：4-6个月
影响因子：0.463
影响因子：0.673
影响因子：0.642
影响因子：1.307
影响因子：0.480
科普杂志阅读、期刊信息查询、论文下载、文献查询、原创文章检测等多项服务。
服务热线
您的论文已提交，我们会尽快联系您，请耐心等待！
你的密码已发送到您的邮箱，请查看！
2024-08-05
46
上传者：管理员
