基于弹性波场矢量分解与低秩分解的地震正演模拟方法
专利汇可以提供基于弹性波场矢量分解与低秩分解的地震正演模拟方法专利检索，专利查询，专利分析的服务。并且本 发明 属于 勘探地球物理学 领域，具体地，涉及一种弹性 地震 波 场正演模拟方法。首先，该方法基于波场矢量分解原理将弹性 地震波 场分解为纵波波场和横波波场，并分别计算相应的纵、横波波场传播矩阵，所得到的传播矩阵包含了对介质参数和数值模拟参数的补偿，因而具有较高的计算 精度 ；然后，基于低秩分解原理将传播矩阵进行分解得到矢量纵波和矢量横波的传播算子，从而降低数值模拟的计算量，提高计算效率；最后，将数值模拟得到的更新后的矢量纵波记录和矢量横波记录进行耦合，得到地震波场正演模拟结果。本发明构建的正演方法得到的地震波场记录几乎不存在数值频散，数值模拟 稳定性 较高，是一种高精度的弹性波场数值模拟方法。，下面是基于弹性波场矢量分解与低秩分解的地震正演模拟方法专利的具体信息内容。
1.一种弹性地震波场正演模拟方法。基于波场矢量分解原理将弹性地震波场分解为一个矢量纵波波场和一个矢量横波波场，并分别计算相应的纵、横波波场传播矩阵；然后，基于矩阵低秩分解原理将传播矩阵进行分解得到矢量纵波和矢量横波的传播算子；最后，将数值模拟得到的更新后的矢量纵波记录和矢量横波记录进行耦合，得到地震波场正演模拟结果，具体包括如下步骤：
步骤1，基于波场矢量分解原理将弹性地震波场分解为一个矢量纵波波场和一个矢量横波波场，根据递归时间积分原理分别计算矢量纵波和矢量横波的地震波场传播矩阵步骤2，基于矩阵低秩分解原理，分别将第1步计算得到的地震波传播矩阵进行分解得到矢量纵波和矢量横波的的地震波场传播算子
步骤3，基于第1步和第2步得到的地震波场传播算子进行弹性波场外推。在每一个时间步长中，分别应用矢量纵波和矢量横波的地震波传播算子计算得到更新后的纵、横波场步骤4，基于第1步、第2步和第3步得到的地震波场传播算子进行弹性波场外推。在每一个时间步长中，对于更新后的纵、横波场，通过纵、横波场耦合得到更新后的地震弹性波场。
2.根据权利要求1所述的基于波场矢量分解与耦合的弹性波正演模拟方法，其特征在于，步骤一具体为：
基于波场矢量分解原理将弹性地震波场分解为一个矢量纵波波场和一个矢量横波波场，根据递归时间积分原理分别计算矢量纵波和矢量横波的的地震波场传播矩阵；
波场矢量分解的公式如下：
P -1 P P
U ＝FFT [A(A·FFT[U])]
S -1 P P
U ＝-FFT [A×(A×FFT[U])]
P
其中， 表示波数域的纵波位移矢量、 表示波数域的横波位移矢量、U 表示空间域S -1 P
的纵波位移矢量、U 表示空间域的横波位移矢量、FFT和FFT 表示傅里叶正反变换、A 表示纵波的偏振矢量
递归时间积分的具体方法如下，以纵波为例：
第一步，假设该质点传播的地震波场可以近似成一个平面波：
其中， 是一个常数，表示P波振幅； 表示沿三
维空间不同方向的归一化波数；ω、t和j分别表示圆周频率，时间和复数单位，并且满足ω＝|k|vp；
第二步，平面P波沿空间的水平分量同样是一个平面波，通过坐标轴映射可以得到该分量的表达式为：
其中，角度θ和 的分别表示平面波传播方向与水平方向和垂直方向的夹角；
第三步，根据平面波公式，可以得到P波沿水平方向的分量的递归时间积分波场外推公式：
同理，横波的波场外推公式表示为：
P S
其中，u 和u 分别表示纵横波位移场；n+1、n、n-1分别表示下一个时刻、当前时刻和过去时刻；vP和vS分别表示纵横波的相速度；ω、Δt和|k|分别表示圆周频率、时间步长和波数。
3.根据权利要求1-2所述的基于波场矢量分解与耦合的弹性波正演模拟方法，其特征在于，步骤二具体为：
通过低秩分解分解原理，将递归时间积分波场外推矩阵进行分解，分别得到纵波和横波的外场外推算子。矩阵低秩分解后的子矩阵由原矩阵的部分行、部分列以及数个常数组成，可以表示为：
其中，矩阵W(x,k)是一个需要进行分解的目标矩阵；W(x,km)表示一个由矩阵W某几列组成的矩阵；W(xn,k)表示一个由矩阵W某几行组成的矩阵；amn表示一个常数矩阵；ML＝W(x,km)·amn，MR＝W(xn,k)。
4.根据权利要求1-3所述的基于弹性波场矢量分解和低秩分解的地震波场正演模拟，其特征在于，步骤三具体为：
第一步，在每一个时间步长中，对地震波场进行傅里叶正变换得到波数域的地震波场其中，U表示空间域的位移矢量、 表示波数域的位移矢量，FFT表示快速傅里叶变换；
第二步，对波数域的地震波场 分别乘以矢量纵波波场外推算子和矢 量横波波场外推算子：
其中， 表示波数域的位移矢量； 表示波数域的纵波位移矢量； 表示波数域的R R
横波位移矢量；n+1、n-1和n表示下一个时间时刻、过去时刻和当前时刻；MP 和MS 分别表示纵波和横波是地震波传播右矩阵；
第三步，对波数域的纵横波矢量 和 进行傅里叶反变换，并乘以地震波传播左矩阵得到更新后的纵横波波场：
P S
其中， 和 表示波数域的纵横波位移矢量，U 和U 表示空间域的纵横波位移矢量；
L L
MP 和MS 分别表示纵波和横波是地震波传播左矩阵；IFFT表示傅里叶反变换。
5.根据权利要求1-4所述的基于波场矢量分解与耦合的弹性波正演模拟方法，其特征在于，步骤四具体为：
将更新后的纵横波位移矢量耦合得到地震波正演模拟记录：
P S
U＝U+U
P S
其中，U 和U 表示空间域的纵横波位移矢量，U表示地震波正演模拟记录。
基于弹性波场矢量分解与低秩分解的地震正演模拟方法
技术领域
背景技术
发明内容
具体实施方式
＝W(x,km)·amn，M ＝W(xn,k)；
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 利用远偏移的增加的炮点间隔的地震勘探 | 2020-06-01 | 263 |
| 基于卷积神经网络的地震反演大数据生成方法 | 2020-05-15 | 823 |
| 陆地地震勘探资料的降维自适应层间多次波压制方法 | 2020-06-04 | 137 |
| 基于裂缝等效介质模型的地震响应模式建立方法 | 2020-05-19 | 234 |
| 直数 | 2020-06-07 | 164 |
| 一种基于伪随机相关辨识技术的电法勘探仪器及其方法 | 2020-06-02 | 879 |
| 一种自适应混合范数字典学习地震波阻抗反演方法 | 2020-05-20 | 229 |
| 基于最大相关熵的地震子波信号提取方法 | 2020-05-11 | 95 |
| 用于估测矿物勘探中的γ-γ测井数据的方法和系统 | 2020-05-31 | 782 |
| 一种减弱煤矿巷道电法体积效应干扰的电流导弹装置 | 2020-05-30 | 319 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
