专利汇可以提供基于弹性波场矢量分解与低秩分解的地震正演模拟方法专利检索,专利查询,专利分析的服务。并且本 发明 属于 勘探地球物理学 领域,具体地,涉及一种弹性 地震 波 场正演模拟方法。首先,该方法基于波场矢量分解原理将弹性 地震波 场分解为纵波波场和横波波场,并分别计算相应的纵、横波波场传播矩阵,所得到的传播矩阵包含了对介质参数和数值模拟参数的补偿,因而具有较高的计算 精度 ;然后,基于低秩分解原理将传播矩阵进行分解得到矢量纵波和矢量横波的传播算子,从而降低数值模拟的计算量,提高计算效率;最后,将数值模拟得到的更新后的矢量纵波记录和矢量横波记录进行耦合,得到地震波场正演模拟结果。本发明构建的正演方法得到的地震波场记录几乎不存在数值频散,数值模拟 稳定性 较高,是一种高精度的弹性波场数值模拟方法。,下面是基于弹性波场矢量分解与低秩分解的地震正演模拟方法专利的具体信息内容。
1.一种弹性地震波场正演模拟方法,其特征在于,包括如下步骤:
步骤1,基于波场矢量分解原理将弹性地震波场分解为一个矢量纵波波场和一个矢量横波波场,根据递归时间积分原理分别计算矢量纵波和矢量横波的地震波场传播矩阵;
步骤2,基于矩阵低秩分解原理,分别将步骤1计算得到的地震波场传播矩阵进行分解得到矢量纵波和矢量横波的地震波场传播算子;
步骤3,基于步骤1和步骤2得到的地震波场传播算子进行弹性波场外推,在每一个时间步长中,分别应用矢量纵波和矢量横波的地震波场传播算子计算得到更新后的纵、横波场;
步骤4,基于步骤1、步骤2和步骤3得到的地震波场传播算子进行弹性波场外推,在每一个时间步长中,对于更新后的纵、横波场,通过纵、横波场耦合得到更新后的地震弹性波场。
2.根据权利要求1所述的弹性地震波场正演模拟方法,其特征在于,
步骤1中,所述的基于波场矢量分解原理将弹性地震波场分解为一个矢量纵波波场和一个矢量横波波场,包括采用如下波场矢量分解公式:
UP=FFT-1[AP(AP·FFT[U])]
US=-FFT-1[AP×(AP×FFT[U])]
其中,UP表示空间域的纵波位移矢量、US表示空间域的横波位移矢量、FFT和FFT-1分别表示傅里叶正变换和反变换、AP表示纵波的偏振矢量、U表示空间域的位移矢量;
步骤1中,所述的根据递归时间积分原理分别计算矢量纵波和矢量横波的地震波场传播矩阵,包括:
采用如下公式计算矢量纵波沿水平方向的分量的递归时间积分波场外推公式:
以及,采用如下公式计算矢量横波的波场外推公式:
P S
其中,u和u分别表示纵横波位移场;n+1、n、n-1分别表示下一个时刻、当前时刻和过去时刻;vP和vS分别表示纵横波的相速度;ω、Δt和|k|分别表示圆周频率、时间步长和波数。
3.根据权利要求2所述的弹性地震波场正演模拟方法,其特征在于,步骤2具体为:
通过低秩分解原理,将递归时间积分波场外推矩阵进行分解,分别得到纵波和横波的波场外推算子,其中,矩阵低秩分解后的子矩阵由原矩阵的部分行、部分列以及数个常数组成,具体表示为:
其中,矩阵W(x,k)是需要进行分解的目标矩阵;W(x,km)表示由矩阵W的若干列组成的矩阵;W(xn,k)表示由矩阵W的若干行组成的矩阵;amn表示常数矩阵;ML=W(x,km)·amn,MR=W(xn,k)。
4.根据权利要求3所述的弹性地震波场正演模拟方法,其特征在于,步骤3包括:
第一步,在每一个时间步长中,对地震波场进行傅里叶正变换得到波数域的地震波场其中,U表示空间域的位移矢量、表示波数域的位移矢量,FFT表示快速傅里叶变换;
第二步,对波数域的地震波场 分别乘以矢量纵波波场外推算子和矢量横波波场外推算子:
其中, 表示波数域的位移矢量; 表示波数域的纵波位移矢量; 表示波数域的横波位移矢量;n+1、n-1和n分别表示下一个时间时刻、过去时刻和当前时刻;MPR和MSR分别表示纵波和横波的地震波传播右矩阵;
第三步,对波数域的纵横波矢量 和 进行傅里叶反变换,并乘以地震波传播左矩阵得到更新后的纵横波波场:
其中, 和 表示波数域的纵横波位移矢量,UP和US表示空间域的纵横波位移矢量;
MPL和MSL分别表示纵波和横波的地震波传播左矩阵;IFFT表示傅里叶反变换。
5.根据权利要求4所述的弹性地震波场正演模拟方法,其特征在于,步骤4具体为:
将更新后的纵横波位移矢量耦合得到地震波正演模拟记录:
U=UP+US
其中,UP和US表示空间域的纵横波位移矢量,U表示地震波正演模拟记录。
标题 | 发布/更新时间 | 阅读量 |
---|---|---|
基于卷积神经网络的地震反演大数据生成方法 | 2020-05-15 | 823 |
基于盆地建模对用于处理地震数据的速度模型的改进 | 2020-06-05 | 1004 |
弹性矢量逆时偏移成像方法 | 2020-05-16 | 747 |
基于裂缝等效介质模型的地震响应模式建立方法 | 2020-05-19 | 234 |
一种自适应混合范数字典学习地震波阻抗反演方法 | 2020-05-20 | 229 |
陆地地震勘探资料的降维自适应层间多次波压制方法 | 2020-06-06 | 907 |
一种基于伪随机相关辨识技术的电法勘探仪器及其方法 | 2020-06-03 | 924 |
一种裂缝性储层的粘弹性参数提取方法 | 2020-05-26 | 154 |
多域分频并行多尺度全波形反演方法 | 2020-05-21 | 579 |
一种减弱煤矿巷道电法体积效应干扰的电流导弹装置 | 2020-05-30 | 319 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。