一种信息驱动的三维地震观测系统的优化方法

文档序号:10652887阅读:279来源:国知局
一种信息驱动的三维地震观测系统的优化方法
【专利摘要】本发明涉及一种信息驱动的三维地震观测系统的优化方法,是以地下地质体非均匀性与尺度等先验信息作为驱动,通过稀疏表示理论,对地震的震源和检波器分布进行空间稀疏化采集和优化,以获得最优的地震采集观测系统。稀疏表示理论及地质体先验信息质量控制评价的优化地震观测系统优化方法,压缩了地震采集的冗余度,通过随机稀疏重建,减少了地震采集的震源和接收排列的数量,为三维高精度地震成像提供基础数据保障,提高了效率,降低了大量三维地震勘探采集成本。
【专利说明】
-种信息驱动的H维地震观测系统的优化方法
技术领域
[0001] 本发明设及一种地震勘探方法,尤其是信息驱动的=维地震观测系统的优化方 法。
【背景技术】:
[0002] =维地震勘探技术是获取地下地质体精细构造的有效方法技术,其对地下地质体 探测的分辨率与地震观测系统中的覆盖面元次数有关,即在规则地震观测系统中,与震源 间距和接收排列线距相关。一般来说,固定面积区域内,震源个数越多,接收排列线数越多, 则对地下地质体的分辨率越高,但运也会导致地震勘探采集有很高的成本。
[0003] 近年来,稀疏表示采样理论已逐渐代替有几十年历史的奈奎斯特采样定律,采样 间隔不再取决于信号的带宽,而是取决于信号的稀疏度及信息在信号中的结构与分布,并 且随机欠采样比规则稀疏采样更有利于信号的重建。实际规则网格地震勘探采集观测时, 震源和接收排列是存在较大的冗余度的,并且运种冗余度与勘探目标地质体的先验信息相 关,如地质体尺度大小、地下介质非均匀程度等相关。
[0004] 目前,复杂构造地质体需要=维地震观测和成像技术才能查明,其采集与处理的 成本制约着技术的应用。

【发明内容】

[0005] 本发明的目的就在于针对上述现有技术的不足,提供一种提高地震勘探采集效 率、优化地震观测系统的一种信息驱动的地震观测系统优化方法。
[0006] 本发明的主要思想是稀疏表示采样理论,对信号随机稀疏采样通过数据重建,亦 能保真重建恢复原信号。通过对地震震源和接收排列的随机欠采样(减少一半的采样点), 及地下模型信息驱动的地震观测系统评价方法,W减少地震采集的震源个数与接收排列 数,节约=维地震勘探的采集成本。
[0007] 本发明的目的是通过W下技术方案实现的:
[000引信息驱动的地震观测系统优化方法,包括下述步骤:
[0009] a、根据地调基础数据,包括地质体尺度信息、非均匀性信息,建立随机函数的=维 地下介质地球物理参数模型;
[0010] b、在勘探区域内,设计地面规则密集网格震源和接收排列,设规则震源个数为化, 规则接收排列条数为Nr,通过波动方程和步骤a建立的地球物理模型,数值模拟规则密集网 格下的地震数据;
[OOW C、稀疏化采集接收排列:将规则排列的检波器随机缺失50%,即接收排列条数降 低为round (Nr/2)条,round取整数函数;
[0012] d、通过波动方程和步骤a建立的地球物理模型,数值模拟计算步骤C稀疏化接收排 列的地震数据;
[0013] e、将步骤d数据进行稀疏约束重建,得到重建后的密集规则网格地震数据,W先验 信息驱动的地质体地震照明分析强度为质量评价标准,质量控制重建后的地震数据,若符 合地质体的勘探目标要求,则重复步骤C,直至不符合地震数据质量控制标准,记录最终的 接收排列条数N rsparse ;
[0014] f、同时稀疏化采集震源:将规则排列的震源随机缺失50%,即激发震源个数降低 为round (化/2)个,round取整数函数;
[0015] g、通过波动方程和步骤a建立的地球物理模型,数值模拟计算步骤f稀疏化震源的 地震数据;
[0016] h、将步骤g数据进行稀疏约束重建,得到重建后的密集规则网格地震数据,W先验 信息驱动的地质体地震照明分析强度为质量评价标准,质量控制重建后的地震数据,若符 合地质体的勘探目标要求,则重复步骤f,直至不符合地震数据质量控制标准,记录最终的 震源个数N SSparse ;
[0017] i、按照新的震源个数化sparse和接收排列条数Nrsparse,设计新的地震观测系统,其中 震源和接收排列随机分布。
[0018] 有益效果:经试验,稀疏表示理论及地质体先验信息质量控制评价的优化地震观 测系统优化方法,压缩了地震采集的冗余度,通过随机稀疏重建,减少了地震采集的震源和 接收排列的数量,为=维高精度地震成像提供基础数据保障,提高了效率,降低了大量=维 地震勘探采集成本。
【附图说明】
[0019] 图1 一种信息驱动的地震观测系统优化方法技术路线图流程图
【具体实施方式】
[0020] 下面结合附图对本发明进一步的详细说明。
[0021] a、根据基础地调数据,包括地质体尺度信息、非均匀性信息等,建立随机函数的= 维地下介质地球物理参数模型;
[0022] b、勘探区域内,设计地面规则密集网格的震源和接收排列,设规则震源个数为化, 规则接收排列条数为Nr,波动方程和步骤a建立的地球物理模型,数值模拟规则密集网格下 的地震数据;
[0023] C、稀疏化采集接收排列:将规则排列的接收排列随机缺失50%,即接收排列条数 降低为round(Nr/2)条,其中round为取整数函数;
[0024] d、通过波动方程和步骤a建立的地球物理模型,数值模拟计算步骤C稀疏化接收排 列的地震数据;
[0025] e、将步骤d数据进行稀疏约束重建,得到重建后的密集规则网格地震数据,W先验 信息驱动的地质体地震照明分析强度为质量评价标准,质量控制重建后的地震数据,若符 合地质体的勘探目标要求,则重复步骤C,直至不符合地震数据质量控制标准,记录最终的 接收排列条数N rsparse ;
[0026] f、同时稀疏化采集震源:将规则排列的震源随机缺失50%,即激发震源个数降低 为round(化/2)个,其中round为取整数函数;
[0027] g、通过波动方程和步骤a建立的地球物理模型,数值模拟计算步骤f稀疏化震源的 地震数据;
[0028] h、将步骤g数据进行稀疏约束重建,得到重建后的密集规则网格地震数据,W先验 信息驱动的地质体地震照明分析强度为质量评价标准,质量控制重建后的地震数据,若符 合地质体的勘探目标要求,则重复步骤f,直至不符合地震数据质量控制标准,记录最终的 震源个数N SSparse ;
[0029] i、按照新的震源个数化sparse和接收排列条数Nrsparse,设计新的地震观测系统,其中 震源和接收排列随机分布。
[0030] 图1是本发明实施方案的技术路线图,其中重点描述随机函数的地球物理参数模 型建立,稀疏约束重建地震波场,W及地震照明分析评价先验信息的地质体探测效果。信息 驱动的地震观测系统优化方法,包括如下步骤:
[0031] a、根据基础地调数据,包括地质体尺度信息、非均匀性信息等,建立随机函数的= 维地下介质地球物理参数模型。S维地球物理模型可由大、小两种尺度参数表示,其中大尺 度参数为地质体的背景介质参数,小尺度参数为地质体的随机扰动参数,其速度v(x,y,z) 和密度P(x,y,z)函数可表示为;
[0032] v(x,y ,Z) =vo+5v(x,y ,Z),
[0033] p(x,y ,z) =po+5p(x,y ,z).
[0034] 其中VQiPQ为背景介质参数,S为随机扰动参数。
[0035] b、勘探区域内,设计地面规则密集网格的震源和接收排列,设规则震源个数为化, 规则接收排列条数为Nr,波动方程和步骤a建立的地球物理模型,数值模拟规则密集网格下 的地震数据;
[0036] C、稀疏化采集接收排列:将规则排列的接收排列随机缺失50%,即接收排列条数 降低为round(Nr/2)条,其中round为取整数函数;
[0037] d、通过波动方程和步骤a建立的地球物理模型,数值模拟计算步骤C稀疏化接收排 列的地震数据;
[0038] e、将步骤d数据进行稀疏约束重建,得到重建后的密集规则网格地震数据,W先验 信息驱动的地质体地震照明分析强度为质量评价标准,质量控制重建后的地震数据,若符 合地质体的勘探目标要求,则重复步骤C,直至不符合地震数据质量控制标准,记录最终的 接收排列条数N rsparse O
[0039] 规格网格数据的地震数据稀疏约束数据重建,可W表达为下式;
[0040] da1:aQbs = T ? da1:a
[0041] 其中data是密集规则测网数据,dataobs是观测数据,T为稀疏采样矩阵。W较稀疏 观测数据来预测密集规则的测网数据,是个欠定问题,一方面控制重建后的数据误差足够 小,另一方面约束变换域数据m是稀疏的,即求解混合范数最优化问题,建立如下稀疏约束 目标函数:
[0042]
[0043] 其中A是控制变换域稀疏度与重建反演误差的折中参数,变量P和q是测量反演误 差和模型稀疏化的范数,常用的提升变换域稀疏性的方式是令,p = 2和q=l。运里面L是任 一种稀疏变换算子。
[0044] 地震照明技术评价信息驱动的地质体的探测效果:地震观测系统的照明强度DU, y,z)可用下式计算,即时间域地震波场沿震源和接收排列维度求和:
[0045]
[0046] 其中G(x,y,z,t)是数值计算的时间域地震波场,At是时间采集间隔,化是采样点 数,Xsi表示第i个震源位置,Xu表示第j条接收排列。地下目标地质体区域内的地震照明分 布均匀,且强度值大于模型平均值,则判断地震观测系统对目标体符合质量要求,可W进一 步的随机缺失震源数量;
[0047] f、同时稀疏化采集震源:将规则排列的震源随机缺失50%,即激发震源个数降低 为round(化/2)个,其中round为取整数函数;
[0048] g、通过波动方程和步骤a建立的地球物理模型,数值模拟计算步骤f稀疏化震源的 地震数据;
[0049] h、将步骤g数据进行稀疏约束重建,得到重建后的密集规则网格地震数据,W先验 信息驱动的地质体地震照明分析强度为质量评价标准,质量控制重建后的地震数据,若符 合地质体的勘探目标要求,则重复步骤f,直至不符合地震数据质量控制标准,记录最终的 震源个数N SSparse O 具体实施过程同步骤e;
[0050] i、按照新的震源个数化sparse和接收排列条数Nrsparse,设计新的地震观测系统,其中 震源和接收排列随机分布。
【主权项】
1. 一种信息驱动的地震观测系统优化方法,其特征在于,包括以下步骤: a、 根据地调基础数据,包括地质体尺度信息、非均匀性信息,建立随机函数的三维地下 介质地球物理参数模型; b、 在勘探区域内,设计地面规则密集网格震源和接收排列,设规则震源个数为Ns,规则 接收排列条数为Nr,通过波动方程和步骤a建立的地球物理模型,数值模拟规则密集网格下 的地震数据; c、 稀疏化采集接收排列:将规则排列的检波器随机缺失50%,即接收排列条数降低为 round (Nr/2)条,round取整数函数; d、 通过波动方程和步骤a建立的地球物理模型,数值模拟计算步骤c稀疏化接收排列的 地震数据; e、 将步骤d数据进行稀疏约束重建,得到重建后的密集规则网格地震数据,以先验信息 驱动的地质体地震照明分析强度为质量评价标准,质量控制重建后的地震数据,若符合地 质体的勘探目标要求,则重复步骤c,直至不符合地震数据质量控制标准,记录最终的接收 排列条数N rsparse ; f、 同时稀疏化采集震源:将规则排列的震源随机缺失50%,即激发震源个数降低为 round (Ns/2)个,round取整数函数; g、 通过波动方程和步骤a建立的地球物理模型,数值模拟计算步骤f稀疏化震源的地震 数据; h、 将步骤g数据进行稀疏约束重建,得到重建后的密集规则网格地震数据,以先验信息 驱动的地质体地震照明分析强度为质量评价标准,质量控制重建后的地震数据,若符合地 质体的勘探目标要求,则重复步骤f,直至不符合地震数据质量控制标准,记录最终的震源 个数N ssparse ; i、 按照新的震源个数Nssparse和接收排列条数Nrsparse,设计新的地震观测系统,其中震源 和接收排列随机分布。
【文档编号】G01V1/00GK106019360SQ201610333130
【公开日】2016年10月12日
【申请日】2016年5月19日
【发明人】巩向博, 王升超, 于晨霞, 韩立国
【申请人】吉林大学
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1