CN110542923A - 一种快速高精度叠后地震阻抗反演方法 - Google Patents

一种快速高精度叠后地震阻抗反演方法 Download PDF

Info

Publication number
CN110542923A
CN110542923A CN201910824513.7A CN201910824513A CN110542923A CN 110542923 A CN110542923 A CN 110542923A CN 201910824513 A CN201910824513 A CN 201910824513A CN 110542923 A CN110542923 A CN 110542923A
Authority
CN
China
Prior art keywords
longitudinal wave
seismic
reflection coefficient
inversion
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910824513.7A
Other languages
English (en)
Other versions
CN110542923B (zh
Inventor
聂文亮
文晓涛
何健
李波
黄伟
李垒
陈芊澍
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing Three Gorges University
Chengdu Univeristy of Technology
Original Assignee
Chongqing Three Gorges University
Chengdu Univeristy of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chongqing Three Gorges University, Chengdu Univeristy of Technology filed Critical Chongqing Three Gorges University
Priority to CN201910824513.7A priority Critical patent/CN110542923B/zh
Publication of CN110542923A publication Critical patent/CN110542923A/zh
Application granted granted Critical
Publication of CN110542923B publication Critical patent/CN110542923B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282—Application of seismic models, synthetic seismograms
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/622—Velocity, density or impedance
    • G01V2210/6226—Impedance

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种快速高精度叠后地震阻抗反演方法，包括：通过对地震褶积模型中的反射系数序列进行奇偶分解，得到反射系数序列分解后的地震褶积模型；利用所述反射系数序列分解后的地震褶积模型及拉格朗日乘子法，构建纵波低频模型约束的反演目标函数；利用L1‑2最小化算法对所述纵波阻抗低频模型同步约束的反演目标函数进行反演，得到奇偶分量系数；利用所述奇偶分量系数，得到纵波阻抗。

Description

一种快速高精度叠后地震阻抗反演方法
技术领域
本发明涉及地震阻抗反演技术领域，特别涉及一种快速高精度叠后地震阻抗反演方法。
背景技术
地震阻抗反演是叠后地震资料解释最为常用的方法之一，在储层描述方面得到了广泛的应用。目前已有反演方法的垂向分辨率不高、横向连续性较差，导致反演效果不好。针对这一问题，国内外学者开展了深入的研究。从本质上讲，现有的反演方法均是基于L1范数的凸优化算法，考虑到地震反演问题本身固有的“病态解”、“多解性”等问题，指出基于L1范数的凸优化算法会导致反演结果陷入次优稀疏解。为了得到最优稀疏解，非凸优化算法引起了众多学者的关注，比如，加权最小二乘迭代的Lp范数最小化算法，L1/L2和L1-2最小化算法。其中，有学者利用L1-2最小化算法对地震衰减进行补偿，并取得较好的实际应用效果。
发明内容
根据本发明实施例提供的方案解决的技术问题是横向连续性差，垂向分辨率低。
根据本发明实施例提供的一种快速高精度叠后地震阻抗反演方法，包括：
通过对地震褶积模型中的反射系数序列进行奇偶分解，得到反射系数序列分解后的地震褶积模型；
利用所述反射系数序列分解后的地震褶积模型及拉格朗日乘子法，构建纵波阻抗低频模型同步约束的反演目标函数；
利用L1-2最小化算法对所述纵波阻抗低频模型同步约束的反演目标函数进行反演，得到奇偶分量系数；
利用所述奇偶分量系数，得到纵波阻抗。
优选地，所述通过对地震褶积模型中的反射系数序列进行奇偶分解，得到反射系数序列分解后的地震褶积模型包括：
通过对所述发射系数序列进行奇偶分解，得到分解后的发射系数序列：r＝Dm；
通过将所述分解后的反射系数序列r＝Dm代入所述地震褶积模型s＝Wr+n中，得到反射系数序列分解后的地震褶积模型：s＝WDm+n；
其中，所述r是指反射系数序列；所述D是指反射系数序列分解算子；所述m是指奇偶分量对应的系数；所述s是指地震记录；所述W是指子波矩阵；所述n为地震噪声序列。
优选地，所述纵波阻抗低频模型同步约束的反演目标函数包括：
其中，所述ξlow是指纵波阻抗低频模型；所述L为低频滤波矩阵；所述C是指积分矩阵。
优选地，所述利用所述反射系数序列分解后的地震褶积模型及拉格朗日乘子法，构建纵波阻抗低频模型同步约束的反演目标函数包括：
根据测井资料，建立纵波阻抗低频模型；
根据所述拉格朗日乘子法以及所述纵波阻抗低频模型对基追踪标准公式进行同步约束，得到所述反演目标函数。
优选地，所述利用L1-2最小化算法对所述纵波阻抗低频模型同步约束的反演目标函数进行反演，得到奇偶分量系数包括：
通过对所述反演目标函数进行化简处理，得到L1-2最小化反演目标函数；
利用增广拉格朗日乘子法对所述L1-2最小化反演目标函数进行约束，得到最终目标函数；
利用交替方法乘子法ADMM对所述最终目标函数进行求解，得到奇偶分量系数。
优选地，所述利用所述奇偶分量系数，得到纵波阻抗包括：
根据所述奇偶分量系数，得到纵波反射系数；
通过对所述纵波反射系数进行积分处理，得到纵波阻抗。
根据本发明实施例提供的方案，使得反演方法的垂向分辨率高、横向连续性好，有效改善了反演效果。
附图说明
图1是本发明实施例提供的一种快速高精度叠后地震阻抗反演方法流程图；
图2是本发明实施例提供的反射系数奇偶分解示意图；
图3是本发明实施例提供的低通滤波矩阵效果示意图；
图4是本发明实施例提供的(a)L1和(b)L1-2范数单位球的示意图。
具体实施方式
以下结合附图对本发明的优选实施例进行详细说明，应当理解，以下所说明的优选实施例仅用于说明和解释本发明，并不用于限定本发明。
图1是本发明实施例提供的一种快速高精度叠后地震阻抗反演方法流程图，如图1所示，包括：
步骤S101：通过对地震褶积模型中的反射系数序列进行奇偶分解，得到反射系数序列分解后的地震褶积模型；
步骤S102：利用所述反射系数序列分解后的地震褶积模型及拉格朗日乘子法，构建纵波阻抗低频模型同步约束的反演目标函数；
步骤S103：利用L1-2最小化算法对所述纵波阻抗低频模型同步约束的反演目标函数进行反演，得到奇偶分量系数；
步骤S104：利用所述奇偶分量系数，得到纵波阻抗。
其中，所述通过对地震褶积模型中的反射系数序列进行奇偶分解，得到反射系数序列分解后的地震褶积模型包括：通过对所述发射系数序列进行奇偶分解，得到分解后的发射系数序列：r＝Dm；通过将所述分解后的反射系数序列r＝Dm代入所述地震褶积模型s＝Wr+n中，得到反射系数序列分解后的地震褶积模型：s＝WDm+n；其中，所述r是指反射系数序列；所述D是指反射系数序列分解算子；所述m是指奇偶分量对应的系数；所述s是指地震记录；所述W是指子波矩阵；所述n为地震噪声序列。
其中，所述纵波阻抗低频模型同步约束的反演目标函数包括：其中，所述ξlow是指纵波阻抗低频模型；所述L为低频滤波矩阵；所述C是指积分矩阵。
其中，所述利用所述反射系数序列分解后的地震褶积模型及拉格朗日乘子法，构建纵波阻抗低频模型同步约束的反演目标函数包括：根据测井资料，建立纵波阻抗低频模型；根据所述拉格朗日乘子法以及所述纵波低频模型对基追踪标准公式进行同步约束，得到所述反演目标函数。
其中，所述利用L1-2最小化算法对所述纵波阻抗低频模型同步约束的反演目标函数进行反演，得到奇偶分量系数包括：通过对所述反演目标函数进行化简处理，得到L1-2最小化反演目标函数；利用增广拉格朗日乘子法对所述L1-2最小化反演目标函数进行约束，得到最终目标函数；利用交替方法乘子法ADMM对所述最终目标函数进行求解，得到奇偶分量系数。
其中，所述利用所述奇偶分量系数，得到纵波阻抗包括：根据所述奇偶分量系数，得到纵波反射系数；通过对所述纵波反射系数进行积分处理，得到纵波阻抗。
本发明实施例首先根据奇偶分解理论将反射系数进行分解，同时采用拉格朗日乘子法，构造出低频模型约束的反演目标函数，使得反演结果符合实际地质背景，然后采用L1-2最小化算法进行反演；最后，通过模型试算和实际地震数据验证了准确性。
1、构建反演目标函数
根据Robinson地震褶积模型，地震记录中每一道地震信号由地震子波和相应反射系数序列褶积而成，写成矩阵形式，如下公式：
s＝Wr+n (1)
式中：s为地震记录；W为子波矩阵；r为反射系数序列；n为地震噪声序列。
1.1、反射系数奇偶双极子分解
根据奇偶分解理论，一个反射系数对可以唯一分解为一个反射系数偶分量和一个反射系数奇分量,如图2所示。图中c,d分别表示顶、底界面反射系数，δ(t)表示单位脉冲函数，Δt表示采样间隔；nΔt表示薄层时间厚度，n为采样点序号。
将反射系数r进行奇偶分解，得到：
r＝are+bro (2)
式中：re为反射系数偶分量；ro为反射系数奇分量；a为偶分量分解系数；b为奇分量分解系数；且a、b是唯一的。
对于任意一道反射系数序列均可表示为：
式中：i为顶层所在采样点位置；M为单道采样点个数；j为反射层厚度的采样点个数；N为反射层的最大厚度对应的采样点个数；ai,j为偶分量系数；bi,j为奇分量系数。
结合公式(2)，将公式(3)表达为矩阵形式：
r＝Dm (4)
式中：D表示反射系数分解算子，D＝[rero]；m表示奇偶分量对应的系数，即待反演量。
结合公式(1)和公式(4)，地震记录表达如下：
s＝WDm+n (5)
1.2、模型约束的反演目标函数
在奇偶反射系数分解的基础上，使用L1范数对反演目标函数进行稀疏约束，得到基追踪求解的标准形式：
式中，λ为正则化调节参数，用来控制待反演参数的稀疏度。在采用公式(6)直接进行反演时，横向连续性差，且反演结果缺乏实际工区地质情况。利用测井资料建立低频模型，对公式(6)进行约束，可有效补偿低频信息，使得反演结果包含来自实际工区的背景情况，可有效改善反演结果的横向连续性。为此，引入低频约束项，仅用测井低频信息来改善反演剖面的横向连续性。
对于叠后阻抗反演，反射系数与纵波阻抗的关系定义为：
式中，Ip(t)为t时刻的纵波阻抗，r(t)为t时刻的反射系数。公式(7)可近似表达为：
对公式(8)两边同时求积分，可得：
将公式(9)简写为矩阵形式有：
Cr＝ξ (10)
式中：C为积分矩阵，ξ为相对波阻抗序列，其表达形式分别为公式(11)和公式(12)：
结合公式(4)，可将(10)改写为：
CDm＝ξ (13)
为了使得反演结果符合实际工区的地质背景，将公式(13)变形为：
LCDm＝ξlow (14)
式中，ξlow为根据测井资料建立的低频模型，L为低频滤波矩阵。
根据拉格朗日乘子法，将公式(14)加入公式(6)中，最终得到低频模型约束的反演目标函数：
1.3、低通滤波矩阵的实现
为了使得反演结果更加贴近实际工区地质背景，采用在目标函数中添加低频地质背景进行约束，但适用于地震反演算法的低频滤波矩阵鲜有文献提及。已有的低通滤波矩阵在滤波的时会带来“截断效应”，为了消除这一影响，采用汉宁窗结合傅里叶正、反变换实现了低通滤波矩阵，但该方法将模型参数进行三倍延长，计算较复杂。为此根据DCT(Discrete Cosine Transform，离散余弦变换)设计了适用于地震反演的低通滤波矩阵，其效果如图2中点虚线所示。将该低通滤波矩阵定义如下：
式中：n、m代表低通滤波矩阵LLP中对应元素的行和列；N为待滤波信号的长度；M由截止频率决定，其定义为：
式中：ωc为截止频率。
2、L1-2最小化反演算法
指出块化的反演结果分辨率更高，对地层的解释能力更强。同时，指出基于L1范数的凸优化算法会导致反演结果陷入次优稀疏解。图4给出了L1和L1-2范数单位球示意图，由图4可知，L1-2最小化算法的求解结果更接近最优稀疏解，因此，本文采用L1-2最小化算法进行地震阻抗反演。
令将公式(15)可化简为L1-2最小化算法的求解形式：
式中，α为权重参数，采用DCA(Difference of Convex Algorithm，凸函数差异算法)将公式(18)分解为F(m)＝G(m)-H(m)，其中：
根据DCA迭代公式，将L1-2最小化反演目标函数简化为：
式中，yk为H(m)在mk处的梯度，其定义为：
引入增广拉格朗日乘数法对公式(20)进行约束，可得：
式中，z为辅助中间变量，w为拉格朗日乘子，ρ为惩罚参数。公式(20)由交替方法乘子法(Alternating Direction Method of Multipliers，ADMM)求解，其迭代递推公式如：
式中，z、m分别由软阈值函数和梯度下降法更新迭代，其求解公式分别为：
zk+1＝shrink(mk+wk/ρ,λ/ρ) (24)
mk+1＝(GTG+ρI)-1(GTs-yk+ρzk+1-wk) (25)
式中，shrink为软阈值函数。公式(23)迭代次数达到上限时，迭代终止。然后，对反演结果进行道积分，即可得出对应的纵波阻抗，其计算公式为：
根据本发明实施例提供的方案，使得反演方法的垂向分辨率高、横向连续性好，有效改善了反演效果。
尽管上文对本发明进行了详细说明，但是本发明不限于此，本技术领域技术人员可以根据本发明的原理进行各种修改。因此，凡按照本发明原理所作的修改，都应当理解为落入本发明的保护范围。

Claims (6)

1.一种快速高精度叠后地震阻抗反演方法，其特征在于，包括：
通过对地震褶积模型中的反射系数序列进行奇偶分解，得到反射系数序列分解后的地震褶积模型；
利用所述反射系数序列分解后的地震褶积模型及拉格朗日乘子法，构建纵波阻抗低频模型同步约束的反演目标函数；
利用L1-2最小化算法对所述纵波阻抗低频模型同步约束的反演目标函数进行反演，得到奇偶分量系数；
利用所述奇偶分量系数，得到纵波阻抗。
2.根据权利要求1所述的方法，其特征在于，所述通过对地震褶积模型中的反射系数序列进行奇偶分解，得到反射系数序列分解后的地震褶积模型包括：
通过对所述发射系数序列进行奇偶分解，得到分解后的发射系数序列：r＝Dm；
通过将所述分解后的反射系数序列r＝Dm代入所述地震褶积模型s＝Wr+n中，得到反射系数序列分解后的地震褶积模型：s＝WDm+n；
其中，所述r是指反射系数序列；所述D是指反射系数序列分解算子；所述m是指奇偶分量对应的系数；所述s是指地震记录；所述W是指子波矩阵；所述n为地震噪声序列。
3.根据权利要求2所述的方法，其特征在于，所述纵波阻抗低频模型同步约束的反演目标函数包括：
其中，所述ξlow是指纵波阻抗低频模型；所述L为低频滤波矩阵；所述C是指积分矩阵。
4.根据权利要求3所述的方法，其特征在于，所述利用所述反射系数序列分解后的地震褶积模型及拉格朗日乘子法，构建纵波阻抗低频模型同步约束的反演目标函数包括：
根据测井资料，建立纵波阻抗低频模型；
根据所述拉格朗日乘子法以及所述纵波阻抗低频模型对基追踪标准公式进行同步约束，得到所述反演目标函数。
5.根据权利要求4所述的方法，其特征在于，所述利用L1-2最小化算法对所述纵波阻抗低频模型同步约束的反演目标函数进行反演，得到奇偶分量系数包括：
通过对所述反演目标函数进行化简处理，得到L1-2最小化反演目标函数；
利用增广拉格朗日乘子法对所述L1-2最小化反演目标函数进行约束，得到最终目标函数；
利用交替方法乘子法ADMM对所述最终目标函数进行求解，得到奇偶分量系数。
6.根据权利要求5所述的方法，其特征在于，所述利用所述奇偶分量系数，得到纵波阻抗包括：
根据所述奇偶分量系数，得到纵波反射系数；
通过对所述纵波反射系数进行积分处理，得到纵波阻抗。
CN201910824513.7A 2019-09-02 2019-09-02 一种快速高精度叠后地震阻抗反演方法 Active CN110542923B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910824513.7A CN110542923B (zh) 2019-09-02 2019-09-02 一种快速高精度叠后地震阻抗反演方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910824513.7A CN110542923B (zh) 2019-09-02 2019-09-02 一种快速高精度叠后地震阻抗反演方法

Publications (2)

Publication Number Publication Date
CN110542923A true CN110542923A (zh) 2019-12-06
CN110542923B CN110542923B (zh) 2021-01-05

Family

ID=68712350

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910824513.7A Active CN110542923B (zh) 2019-09-02 2019-09-02 一种快速高精度叠后地震阻抗反演方法

Country Status (1)

Country Link
CN (1) CN110542923B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111366975A (zh) * 2020-03-21 2020-07-03 西华师范大学 基于交叉梯度正则化约束的叠前地震ava反演方法
CN112578443A (zh) * 2020-12-03 2021-03-30 成都理工大学 一种叠前地震反演的方法
CN112630834A (zh) * 2020-12-03 2021-04-09 重庆三峡学院 一种叠前地震直接反演的方法
CN112630835A (zh) * 2020-12-03 2021-04-09 重庆三峡学院 一种高分辨率的叠后地震波阻抗反演的方法
CN112859170A (zh) * 2021-03-24 2021-05-28 西南石油大学 一种基于时域高阶有限差分法的高精度波场数值模拟方法
CN113219530A (zh) * 2020-02-05 2021-08-06 中国石油天然气股份有限公司 非稳态盲反褶积方法及装置
CN113640871A (zh) * 2021-08-10 2021-11-12 成都理工大学 一种基于重加权l1范数稀疏约束的地震波阻抗反演方法
CN113885075A (zh) * 2021-09-15 2022-01-04 成都理工大学 一种叠后地震波相对波阻抗反演方法及系统
CN114076980A (zh) * 2020-08-17 2022-02-22 中国石油化工股份有限公司 一种针对薄层刻画的方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160109605A1 (en) * 2014-10-20 2016-04-21 Schlumberger Technology Corporation Fast Model Based Inversion of Acoustic Impedance of Annulus Behind Casing
CN106054247A (zh) * 2016-05-25 2016-10-26 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于转换波地震数据求取高精度反射系数的方法
CN107193046A (zh) * 2016-03-14 2017-09-22 中国石油化工股份有限公司 一种基于谱反演的砂体厚度预测方法及系统
CN107305600A (zh) * 2016-04-21 2017-10-31 新疆维吾尔自治区煤炭科学研究所 最小二乘法电阻率三维近似反演技术

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160109605A1 (en) * 2014-10-20 2016-04-21 Schlumberger Technology Corporation Fast Model Based Inversion of Acoustic Impedance of Annulus Behind Casing
CN107193046A (zh) * 2016-03-14 2017-09-22 中国石油化工股份有限公司 一种基于谱反演的砂体厚度预测方法及系统
CN107305600A (zh) * 2016-04-21 2017-10-31 新疆维吾尔自治区煤炭科学研究所 最小二乘法电阻率三维近似反演技术
CN106054247A (zh) * 2016-05-25 2016-10-26 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于转换波地震数据求取高精度反射系数的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘晓晶 等: "基于正交匹配追踪算法的叠前地震反演方法", 《石油地球物理勘探》 *
文晓涛 等: "低频稀疏双约束宽频带地震阻抗反演", 《地质勘探》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113219530A (zh) * 2020-02-05 2021-08-06 中国石油天然气股份有限公司 非稳态盲反褶积方法及装置
CN113219530B (zh) * 2020-02-05 2023-09-26 中国石油天然气股份有限公司 非稳态盲反褶积方法及装置
CN111366975A (zh) * 2020-03-21 2020-07-03 西华师范大学 基于交叉梯度正则化约束的叠前地震ava反演方法
CN114076980B (zh) * 2020-08-17 2024-04-02 中国石油化工股份有限公司 一种针对薄层刻画的方法及系统
CN114076980A (zh) * 2020-08-17 2022-02-22 中国石油化工股份有限公司 一种针对薄层刻画的方法及系统
CN112578443A (zh) * 2020-12-03 2021-03-30 成都理工大学 一种叠前地震反演的方法
CN112630834A (zh) * 2020-12-03 2021-04-09 重庆三峡学院 一种叠前地震直接反演的方法
CN112630835A (zh) * 2020-12-03 2021-04-09 重庆三峡学院 一种高分辨率的叠后地震波阻抗反演的方法
CN112630835B (zh) * 2020-12-03 2023-10-17 重庆三峡学院 一种高分辨率的叠后地震波阻抗反演的方法
CN112859170B (zh) * 2021-03-24 2022-08-02 西南石油大学 一种基于时域高阶有限差分法的高精度波场数值模拟方法
CN112859170A (zh) * 2021-03-24 2021-05-28 西南石油大学 一种基于时域高阶有限差分法的高精度波场数值模拟方法
CN113640871B (zh) * 2021-08-10 2023-09-01 成都理工大学 一种基于重加权l1范数稀疏约束的地震波阻抗反演方法
CN113640871A (zh) * 2021-08-10 2021-11-12 成都理工大学 一种基于重加权l1范数稀疏约束的地震波阻抗反演方法
CN113885075A (zh) * 2021-09-15 2022-01-04 成都理工大学 一种叠后地震波相对波阻抗反演方法及系统

Also Published As

Publication number Publication date
CN110542923B (zh) 2021-01-05

Similar Documents

Publication Publication Date Title
CN110542923B (zh) 一种快速高精度叠后地震阻抗反演方法
CN110542924B (zh) 一种高精度的纵横波阻抗反演方法
CN108919347A (zh) 基于vmd的地震信号随机噪声压制方法
CN110208862B (zh) 一种基于混合高阶分数阶ATpV稀疏正则化的地震反演方法
CN112630835B (zh) 一种高分辨率的叠后地震波阻抗反演的方法
CN111368247B (zh) 基于快速正交字典的稀疏表征正则化叠前avo反演方法
CN106772585B (zh) 一种基于弹性波解耦方程的优化拟解析方法及装置
CN104914466A (zh) 一种提高地震资料分辨率的方法
Yang et al. Seislet-based morphological component analysis using scale-dependent exponential shrinkage
CN110244353B (zh) 一种基于稀疏范数优化算法的地震数据规则化方法
CN101540039B (zh) 单帧图象超分辨方法
Gao et al. Deep learning vertical resolution enhancement considering features of seismic data
Sun et al. Full-waveform inversion using a learned regularization
Li et al. An attention‐guided convolution neural network for denoising of distributed acoustic sensing–vertical seismic profile data
Xie et al. A depth migration method based on the full-wave reverse-time calculation and local one-way propagation
CN111028160A (zh) 一种基于卷积神经网络的遥感图像噪声抑制方法
Gou et al. Complex seismic wavefield interpolation based on the Bregman iteration method in the sparse transform domain
CN109856673B (zh) 一种基于优势频率迭代加权的高分辨Radon变换数据分离技术
CN109085649B (zh) 一种基于小波变换优化的地震资料去噪方法
Gholami et al. Simultaneous constraining of model and data smoothness for regularization of geophysical inverse problems
US9014985B2 (en) System and method for compensating time and offset varying near-surface effects in seismic data background
CN109490954A (zh) 波场正演模拟方法及装置
CN113341460B (zh) 一种基于连续算子分裂的循环式极小化地震数据重建方法
CN116859461B (zh) 一种多次波成像方法及系统
Sang et al. Research Article Noise Attenuation of Seismic Data via Deep Multiscale Fusion Network

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant