首页 >
论文范文 >
工程工业论文 >
工程综合论文 >
矿业工程论文 >
频率—波数域水陆检数据上、下行波场分离方法
0
给我们留言吧!
摘要：随着海洋油气勘探开发的发展，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.
分享：
黄金是全球战略性关键矿产资源,为我国经济社会发展提供了重要支撑｡卡林型金矿是全球分布较广泛的微细浸染型金矿床,主要分布于美国内华达州､犹他州以及我国的滇黔桂､川陕甘两个金三角内｡卡林型金矿是一种难采选矿床,金主要以细微粒浸染状产于碳酸岩中,必须通过加压预氧化等方式使金暴露出来｡
2025-05-07
煤炭在国内能源结构中占据着主导地位,对于保障国家能源安全和经济发展具有不可替代的作用｡煤炭地质勘探是煤炭资源开发的前期基础工作,其目的是查明煤炭资源储量､质量､分布规律以及开采技术条件等,为煤炭资源的合理开发和利用提供科学依据｡
2025-05-07
根据对矿山中岩土､矿石种类､分布以及矿山内部具体形状等多方面情况来进行勘察钻探,对地质构造､岩层赋存情况､矿体规模等进行综合分析,本文通过实例阐述利用岩土工程勘察及钻探技术对某矿区的水文､工程､环境地址请进行分析,不仅能够促进该矿区开采计划的顺利制定,而且还能够在矿山开采的后续工作中将矿石的开发利用率最大化,从根本上有利于推动矿山开采相关行业的快速发展｡
2025-05-07
地采矿山的井下通风管理一直是一个复杂而关键的问题｡由于井下空间相对狭小,自然光照不足,以及可能出现的各种地质灾害,如坍塌､涌水等,都给通风管理带来了极大的困难[1]｡此外,随着矿山开采深度的不断增加,供氧量的需求也在逐步上升,这使得矿井通风管理面临着前所未有的挑战｡
2025-05-07
在当前采矿充填工艺中,充填材料和充填工艺的选择,给矿山企业实际操作造成了诸多困扰,采空区充填充实率和充填工作效率难题凸显,鉴于这种情况,结合矿山回采工艺和流程以及一般顶板垮落特性,对充填材料和支护工艺进行改良｡
2025-05-07
当前多数选矿厂的一段磨矿分级系统由球磨机､渣浆泵和水力旋流器组成｡球磨机磨矿过程中产生的碎钢球除了可利用圆筒筛部分排除外,通常随矿浆进入泵池经泵进入水力旋流器分级｡由于碎钢球粒度粗､硬度大,容易造成渣浆泵磨损和故障,加剧水力旋流器磨损｡
2025-05-05
对于含向斜工作面开采,大部分煤矿采用跳采方式,即上下2个切眼的刀把式布置,等工作面正常后再布置全长工作面｡该方式避开了恶劣的向斜构造,有利于设备正常运转及煤炭安全开采,但也丢失了大量煤炭资源,同时刀把式工作面衔接时准备工序繁多､对接工艺复杂,对技术要求较高[1-3]｡
2025-04-09
倾斜煤层一般是指煤层倾角大于25°的煤层｡倾斜煤层在我国分布广泛,特别是在我国西部部分矿区,其储量占比达到煤炭资源总量的60%以上,具有巨大的开采价值｡近年来,随着中东部地区赋存条件较好的煤炭资源枯竭殆尽,为了保持矿区的可持续发展,就必须开采倾斜煤层｡
2025-04-09
在煤矿现场,瓦斯抽采管网复杂,积水的地点分散,大部分煤矿靠人工定时巡检放水,存在工作量大､放水不及时､放水器放不出水等问题｡许多专家学者对自动排出瓦斯抽采管路中的积水进行了大量的研究,祝钊等[3]通过对瓦斯抽采管路智能放水系统研究,采用PLC､电磁阀､浮球液位计等设计了一种自动放水系统,能够有效排出瓦斯管路中的积水。
2025-04-09
蚁群算法凭借其卓越的鲁棒性､显著的正反馈效应及良好的兼容性,在路径规划领域展现了广泛的应用潜力｡蚁群算法存在着收敛速度慢､易陷入局部最优的缺点｡为克服这些不足,王志中[3]将启发式信息矩阵引入,兼顾了前期大范围搜索和后期快速收敛;WenxiangGao等[4]设计一种新的信息素扩散梯度公式,增强信息素作用｡
2025-04-07我要评论
期刊名称：矿业工程研究
期刊人气：739
主管单位：湖南省教育厅
主办单位：湖南科技大学
出版地方：湖南
专业分类：煤矿
国际刊号：1674-5876
国内刊号：43-1493/TD
邮发代号：42-359
创刊时间：1980年
发行周期：季刊
期刊开本：16开
见刊时间：4-6个月
影响因子：0.463
影响因子：0.673
影响因子：0.642
影响因子：1.307
影响因子：0.480
学术期刊咨询、科普杂志阅读、专利申请方案、课题研究设计、专著挂名咨询、论文下载、文献查询、原创文章检测、SCI实验设计方案等多项服务。
服务热线
您的论文已提交，我们会尽快联系您，请耐心等待！
你的密码已发送到您的邮箱，请查看！
2024-08-05
25
上传者：管理员
