CN110857997A - 一种基于横向约束的分步叠前弹性参数反演方法及系统 - Google Patents
一种基于横向约束的分步叠前弹性参数反演方法及系统 Download PDFInfo
- Publication number
- CN110857997A CN110857997A CN201810966875.5A CN201810966875A CN110857997A CN 110857997 A CN110857997 A CN 110857997A CN 201810966875 A CN201810966875 A CN 201810966875A CN 110857997 A CN110857997 A CN 110857997A
- Authority
- CN
- China
- Prior art keywords
- wave
- inversion
- change rate
- transverse
- constraint
- 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.)
- Pending
Links
Images
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. analysis, for interpretation, for correction
-
- 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. analysis, for interpretation, for correction
- 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
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/324—Filtering
-
- 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/624—Reservoir parameters
- G01V2210/6242—Elastic parameters, e.g. Young, Lamé or Poisson
Abstract
公开了一种基于横向约束的分步叠前弹性参数反演方法及系统,该方法及系统包括:基于经过同相轴匹配一致的PP波和PS波角道集,获得PP波和PS波角度域反射系数;利用Stewart反射系数公式进行多波联合反演,并且利用卡尔曼滤波算法对反演方程进行横向约束,获得纵波速度变化率和横波速度变化率;基于所述纵波速度变化率和横波速度变化率,代入Aki‑Richards三参数反射系数公式,并同时利用卡尔曼滤波算法对反演方程进行横向约束,得到密度变化率;对获得的反演结果进行低频补偿获得最终的叠前弹性参数的绝对值。本发明通过基于横向约束的分步AVO反演方法,对地下储层识别及流体预测有重要意义。
Description
技术领域
本发明涉及地震资料AVO叠前弹性参数反演与流体识别领域,更具体地,涉及一种基于横向约束的分步叠前弹性参数反演方法及系统。
背景技术
叠前地震反演是储层预测和流体识别的主要技术之一,纵、横波速度和密度等弹性参数可以研究储层岩性和孔隙流体性质,其中由反演得到的密度参数是检测储层含油气性的重要参数。常规的纵波叠前地震反演方法仅使用了反射纵波地震数据,由纵波叠前地震数据提取纵、横波速度和密度参数存在多解性、反演过程不适定等问题,因此人们试图寻找一种稳定的高精度叠前地震反演方法。
AVO反演作为一种叠前反演方法已经成为复杂油气藏储层预测的关键技术,弹性参数反演结果的准确性是AVO叠前反演技术的关键难题。部分实际地震资料信噪比、分辨率不高,PP波反演结果不稳定、反演过程不适定等问题亟需解决。因此,有必要开发一种基于横向约束的分步叠前弹性参数反演方法及系统。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种基于横向约束的分步叠前弹性参数反演方法及系统,其能够通过对密度反演结果的准确性问题,提出了基于卡尔曼滤波约束下的分步多波AVO反演方法研究,解决部分实际地震资料信噪比、分辨率不高,PP波反演结果不稳定、反演过程不适定等问题。
根据本发明的一方面,提出了一种基于横向约束的分步叠前弹性参数反演方法。所述方法可以包括:
1)基于经过同相轴匹配一致的PP波和PS波角道集,获得PP波和PS波角度域反射系数;
2)利用Stewart反射系数公式进行多波联合反演,并且利用卡尔曼滤波算法对反演方程进行横向约束,获得纵波速度变化率和横波速度变化率;
3)基于所述纵波速度变化率和横波速度变化率,代入Aki-Richards三参数反射系数公式,并同时利用卡尔曼滤波算法对反演方程进行横向约束,得到密度变化率;
4)对步骤2)和步骤3)中获得的反演结果进行低频补偿获得最终的叠前弹性参数的绝对值。
其中,RPP、RPS分别为纵波、转换波反射系数;VP=(VP1+VP2)/2,ΔVP=VP2-VP1,VP1、VP2分别为界面上、下层的纵波速度;VS=(VS1+VS2)/2,ΔVS=VS2-VS1,VS1、VS2分别为界面上、下层的SV波速度;ρ=(ρ1+ρ2)/2,Δρ=ρ2-ρ1,ρ1、ρ2分别为界面上、下层的密度;θ=(θ1+θ2)/2,θ1、θ2分别为P波入射角和透射角; 分别为转换波反射角和透射角。
其中,R'PP为新推导得到的反射系数,VP=(VP1+VP2)/2,ΔVP=VP2-VP1,VP1、VP2分别为界面上、下层的纵波速度;VS=(VS1+VS2)/2,ΔVS=VS2-VS1,VS1、VS2分别为界面上、下层的SV波速度,θ=(θ1+θ2)/2,θ1、θ2分别为P波入射角和透射角。
优选地,在步骤2)中,通过如下步骤进行横向约束:
2-1)根据第k道反演的纵、横波速度变化率预测第k+1道的反演结果;
2-2)基于卡尔曼增益利用测量方程对步骤3-1)中预测的第k+1道反演结果进行修正进而获得最优化反演结果。
优选地,在步骤3)中,通过如下步骤进行横向约束:
3-1)根据第k道反演的密度变化率预测第k+1道的反演结果;
3-2)基于卡尔曼增益利用测量方程对步骤3-1)中预测的第k+1道反演结果进行修正进而获得最优化反演结果。
优选地,通过如下公式预测第k+1道的反演结果:
mk+1,k=Ak,kmk,k+rk (4)
其中,mk+1,k为由第k道预测的第k+1道的反演参数,Ak为状态变换矩阵,mk,k为由第k道反演的纵、横波速度变化率和密度变化率的最优结果,rk为mk,k和mk+1,k之间的状态转换误差。
优选地,所述测量方程为:
dk=Gmk+lk (5)
其中,dk、mk分别为AVO反演方程的地震记录与弹性参数,lk为地震数据反演结果的误差值,G为系数矩阵。
优选地,所述卡尔曼增益表达式为:
Zk+1=Yk+1,kGk+1 T(Gk+1Yk+1,kGk+1 T+B)-1 (6)
其中,Zk+1是第k+1道的卡尔曼增益,Yk+1,k为mk+1,k的预测误差协方差矩阵,B是lk的协方差矩阵,Gk+1为系数矩阵。
优选地,通过如下公式对获得最优化反演结果:
mk+1,k+1=mk+1,k+Zk+1(dk+1-Gk+1mk+1,k) (7)
其中,mk+1,k+1为第k+1道的最终优化反演结果,Z为第k+1道卡尔曼增益,Gk+1为系数矩阵,mk+1,k为由第k道预测的第k+1道的反演参数,dk+1是第k+1道地震数据。
根据本发明的另一方面,提出了一种基于横向约束的分步叠前弹性参数反演系统,其上存储有计算机程序,其中,所述程序被处理器执行时实现以下步骤:
步骤1:基于经过同相轴匹配一致的PP波和PS波角道集,获得PP波和PS波角度域反射系数;
步骤2:利用Stewart反射系数公式进行多波联合反演,并且利用卡尔曼滤波算法对反演方程进行横向约束,获得纵波速度变化率和横波速度变化率;
步骤3:基于所述纵波速度变化率和横波速度变化率,代入Aki-Richards三参数反射系数公式,并同时利用卡尔曼滤波算法对反演方程进行横向约束,得到密度变化率;
步骤4:对步骤2和步骤3中获得的反演结果进行低频补偿获得最终的叠前弹性参数的绝对值。
本发明具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的,或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述,这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的一种基于横向约束的分步叠前弹性参数反演方法的步骤的流程图;
图2(a)示出了深度道集CDP 1000点的纵波角道集;
图2(b)示出了深度道集CDP 1000点的转换横波角道集;
图3(a)示出了AVO反演后经过低频补偿的M测线目标层纵波速度剖面;
图3(b)示出了AVO反演后经过低频补偿的M测线目标层横波速度剖面;
图3(c)示出了AVO反演后经过低频补偿的M测线目标层密度剖面。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的一种基于横向约束的分步叠前弹性参数反演方法的步骤的流程图。
在该实施例中,根据本发明的一种基于横向约束的分步叠前弹性参数反演方法的可以包括:
步骤101,基于经过同相轴匹配一致的PP波和PS波角道集,获得PP波和PS波角度域反射系数;
步骤102,利用Stewart反射系数公式进行多波联合反演,并且利用卡尔曼滤波算法对反演方程进行横向约束,获得纵波速度变化率和横波速度变化率;
其中,RPP、RPS分别为纵波、转换波反射系数;VP=(VP1+VP2)/2,ΔVP=VP2-VP1,VP1、VP2分别为界面上、下层的纵波速度;VS=(VS1+VS2)/2,ΔVS=VS2-VS1,VS1、VS2分别为界面上、下层的SV波速度;ρ=(ρ1+ρ2)/2,Δρ=ρ2-ρ1,ρ1、ρ2分别为界面上、下层的密度;θ=(θ1+θ2)/2,θ1、θ2分别为P波入射角和透射角; 分别为转换波反射角和透射角。
在一个示例中,通过如下步骤进行横向约束:
2-1)根据第k道反演的纵、横波速度变化率预测第k+1道的反演结果;
2-2)基于卡尔曼增益利用测量方程对步骤2-1)中预测的第k+1道反演结果进行修正进而获得最优化反演结果。
具体地,在步骤2-1)中通过如下公式预测第k+1道的反演结果:
mk+1,k=Ak,kmk,k+rk (4)
其中,mk+1,k为由第k道预测的第k+1道的反演参数,Ak为状态变换矩阵,mk,k为由第k道反演的纵、横波速度变化率最优结果,rk为mk,k和mk+1,k之间的状态转换误差。
步骤2-2)中所述测量方程为:
dk=Gmk+lk (5)
其中,dk、mk分别为AVO反演方程的地震记录与弹性参数,lk为地震数据反演结果的误差值,G为公式(1)、(2)、(3)式中方程的系数矩阵。
在步骤2-2)所述卡尔曼增益表达式为:
Zk+1=Yk+1,kGk+1 T(Gk+1Yk+1,kGk+1 T+B)-1 (6)
其中,Zk+1是第k+1道的卡尔曼增益,Yk+1,k为mk+1,k的预测误差协方差矩阵,B是lk的协方差矩阵,Gk+1为(1),(2),(3)式中方程的系数矩阵。
在步骤2-2)中通过如下公式对获得最优化反演结果:
mk+1,k+1=mk+1,k+Zk+1(dk+1-Gk+1mk+1,k) (7)
其中,mk+1,k+1为第k+1道的最终优化反演结果,Zk+1为第k+1道卡尔曼增益,Gk+1为系数矩阵,mk+1,k为由第k道预测的第k+1道的反演参数,dk+1是第k+1道地震数据。
步骤103,基于所述纵波速度变化率和横波速度变化率,代入
Aki-Richards三参数反射系数公式,并同时利用卡尔曼滤波算法对反演方程进行横向约束,得到密度变化率;
其中,R'PP为推导得到的新反射系数方程,VP=(VP1+VP2)/2,ΔVP=VP2-VP1,VP1、VP2分别为界面上、下层的纵波速度;VS=(VS1+VS2)/2,ΔVS=VS2-VS1,VS1、VS2分别为界面上、下层的SV波速度,θ=(θ1+θ2)/2,θ1、θ2分别为P波入射角和透射角。
在一个示例中,通过如下步骤进行横向约束:
3-1)根据第k道反演的密度变化率预测第k+1道的反演结果;
3-2)基于卡尔曼增益利用测量方程对步骤3-1)中预测的第k+1道反演结果进行修正进而获得最优化反演结果。
具体地,在步骤3-1)中通过如下公式预测第k+1道的反演结果:
mk+1,k=Ak,kmk,k+rk (4)
其中,mk+1,k为由第k道预测的第k+1道的反演参数,Ak为状态变换矩阵,mk,k为由第k道反演的密度变化率最优结果,rk为mk,k和mk+1,k之间的状态转换误差。
在步骤3-2)中所述测量方程为:
dk=Gmk+lk (5)
其中,dk、mk分别为AVO反演方程的地震记录与弹性参数,lk为地震数据反演结果的误差值,G为公式(1)、(2)、(3)式中方程的系数矩阵。
在步骤3-2)所述卡尔曼增益表达式为:
Zk+1=Yk+1,kGk+1 T(Gk+1Yk+1,kGk+1 T+B)-1 (6)
其中,Zk+1是第k+1道的卡尔曼增益,Yk+1,k为mk+1,k的预测误差协方差矩阵,B是lk的协方差矩阵,Gk+1为公式(1)、(2)、(3)式中方程的系数矩阵。
在步骤3-2)中通过如下公式对获得最优化反演结果:
mk+1,k+1=mk+1,k+Zk+1(dk+1-Gk+1mk+1,k) (7)
其中,mk+1,k+1为第k+1道的最终优化反演结果,Z为第k+1道卡尔曼增益,Gk+1为系数矩阵,mk+1,k为由第k道预测的第k+1道的反演参数,dk+1是第k+1道地震数据。
步骤104,对步骤102和步骤103中获得的反演结果进行低频补偿获得最终的叠前弹性参数的绝对值。
在一个示例中,低频补偿是通过对反演获得的弹性参数结果的反射率计算即可获取参数的绝对值,这个计算的过程将低频成分包含进来。
应用示例
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
本发明选取了中国LJ地区A油田M测线的纵波和转换横波地震数据进行测试。该区在地质上位于凹陷南缓斜坡,受燕山运动和喜马拉雅运动的影响,形成了向北倾没的大型鼻状构造,油气藏类型主要为受构造控制的断块油气藏,兼有岩性、泥岩裂缝等特殊油气藏。如图2(a)所示,为原始地震记录的深度道集CDP 1000点的纵波角道集,图2(b)为深度道集CDP 1000点的转换波角道集,采用的地震资料PP波主频接近35Hz,PS波主频接近15Hz。在时间域内将纵波地震数据与转换横波地震数据匹配一致得到纵波角道集和转换横波角道集。
经过本发明中的基于横向约束的基于横向约束的分步叠前弹性参数反演方法反演后获得如图3(a)中所示的经过低频补偿的M测线目标层段纵波剖面,如图3(b)中所示的经过低频补偿的M测线目标层段横波剖面和如图3(c)中所示的经过低频补偿的M测线目标层段密度剖面。根据图中可以清晰的看出反演的目标层,经与生产井信息进行对比发现,纵波速度,即图3(a)中黑色线段标注的生产井段内的油层和密度反演结果,即图3(c)中黑色线段标注的生产井段内的油层接近真实值范围,实现了较为准确的预测。
本发明还提出了一种基于横向约束的分步叠前弹性参数反演系统,其上存储有计算机程序,其中,所述程序被处理器执行时实现以下步骤:
步骤1:基于经过同相轴匹配一致的PP波和PS波角道集,获得PP波和PS波角度域反射系数;
步骤2:利用Stewart反射系数公式进行多波联合反演,并且利用卡尔曼滤波算法对反演方程进行横向约束,获得纵波速度变化率和横波速度变化率;
步骤3:基于所述纵波速度变化率和横波速度变化率,代入Aki-Richards三参数反射系数公式,并同时利用卡尔曼滤波算法对反演方程进行横向约束,得到密度变化率;
步骤4:对步骤2和步骤3中获得的反演结果进行低频补偿获得最终的叠前弹性参数的绝对值。
综上所述,本发明通过针对部分实际地震资料信噪比、分辨率不高,PP波反演结果不稳定、反演过程不适定等问题,本发明提出了一种基于横向约束的分步叠前弹性参数反演方法及系统,能够较好的反演地下储层中密度等弹性参数,其基于卡尔曼滤波算法则增强横向连续性,理论严谨可靠且操作流程简单实用,为AVO反演识别提供了帮助。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (10)
1.一种基于横向约束的分步叠前弹性参数反演方法,其特征在于,包括:
1)基于经过同相轴匹配一致的PP波和PS波角道集,获得PP波和PS波角度域反射系数;
2)利用Stewart反射系数公式进行多波联合反演,并且利用卡尔曼滤波算法对反演方程进行横向约束,获得纵波速度变化率和横波速度变化率;
3)基于所述纵波速度变化率和横波速度变化率,代入Aki-Richards三参数反射系数公式,并同时利用卡尔曼滤波算法对反演方程进行横向约束,得到密度变化率;
4)对步骤2)和步骤3)中获得的反演结果进行低频补偿获得最终的叠前弹性参数的绝对值。
2.根据权利要求1所述的基于横向约束的分步叠前弹性参数反演方法,其中,在步骤2)中利用以下Stewart纵波反射系数公式和Stewart转换波反射系数公式获得纵波速度变化率和横波速度变化率
4.根据权利要求1中所述的基于横向约束的分步叠前弹性参数反演方法,其中,在步骤2)中,通过如下步骤进行横向约束:
2-1)根据第k道反演的纵、横波速度变化率预测第k+1道的反演结果;
2-2)基于卡尔曼增益利用测量方程对步骤2-1)中预测的第k+1道反演结果进行修正进而获得最优化反演结果。
5.根据权利要求1中所述的基于横向约束的分步叠前弹性参数反演方法,其中,在步骤3)中,通过如下步骤进行横向约束:
3-1)根据第k道反演的密度变化率预测第k+1道的反演结果;
3-2)基于卡尔曼增益利用测量方程对步骤3-1)中预测的第k+1道反演结果进行修正进而获得最优化反演结果。
6.根据权利要求4或5中所述的基于横向约束的分步叠前弹性参数反演方法,其中,通过如下公式预测第k+1道的反演结果:
mk+1,k=Ak,kmk,k+rk (4)
其中,mk+1,k为由第k道预测的第k+1道的反演参数,Ak为状态变换矩阵,mk,k为由第k道反演的纵、横波速度变化率和密度变化率的最优结果,rk为mk,k和mk+1,k之间的状态转换误差。
7.根据权利要求4或5中所述的基于横向约束的分步叠前弹性参数反演方法,其中,所述测量方程为:
dk=Gmk+lk (5)
其中,dk、mk分别为AVO反演方程的地震记录与弹性参数,lk为地震数据反演结果的误差值,G为系数矩阵。
8.根据权利要求4或5中所述的基于横向约束的分步叠前弹性参数反演方法,其中,所述卡尔曼增益表达式为:
Zk+1=Yk+1,kGk+1 T(Gk+1Yk+1,kGk+1 T+B)-1 (6)
其中,Zk+1是第k+1道的卡尔曼增益,Yk+1,k为mk+1,k的预测误差协方差矩阵,B是lk的协方差矩阵,Gk+1为系数矩阵。
9.根据权利要求4或5中所述的基于横向约束的分步叠前弹性参数反演方法,其中,通过如下公式对获得最优化反演结果:
mk+1,k+1=mk+1,k+Zk+1(dk+1-Gk+1mk+1,k) (7)
其中,mk+1,k+1为第k+1道的最终优化反演结果,Z为第k+1道卡尔曼增益,Gk+1为系数矩阵,mk+1,k为由第k道预测的第k+1道的反演参数,dk+1是第k+1道地震数据。
10.一种基于横向约束的分步叠前弹性参数反演系统,其上存储有计算机程序,其中,所述程序被处理器执行时实现以下步骤:
步骤1:基于经过同相轴匹配一致的PP波和PS波角道集,获得PP波和PS波角度域反射系数;
步骤2:利用Stewart反射系数公式进行多波联合反演,并且利用卡尔曼滤波算法对反演方程进行横向约束,获得纵波速度变化率和横波速度变化率;
步骤3:基于所述纵波速度变化率和横波速度变化率,代入Aki-Richards三参数反射系数公式,并同时利用卡尔曼滤波算法对反演方程进行横向约束,得到密度变化率;
步骤4:对步骤2和步骤3中获得的反演结果进行低频补偿获得最终的叠前弹性参数的绝对值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810966875.5A CN110857997A (zh) | 2018-08-23 | 2018-08-23 | 一种基于横向约束的分步叠前弹性参数反演方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810966875.5A CN110857997A (zh) | 2018-08-23 | 2018-08-23 | 一种基于横向约束的分步叠前弹性参数反演方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110857997A true CN110857997A (zh) | 2020-03-03 |
Family
ID=69635141
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810966875.5A Pending CN110857997A (zh) | 2018-08-23 | 2018-08-23 | 一种基于横向约束的分步叠前弹性参数反演方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110857997A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113900140A (zh) * | 2021-09-30 | 2022-01-07 | 中国石油大学(北京) | 一种基于时空组合的地震数据优化的方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103293552A (zh) * | 2013-05-24 | 2013-09-11 | 中国石油天然气集团公司 | 一种叠前地震资料的反演方法及系统 |
CN103675907A (zh) * | 2012-09-20 | 2014-03-26 | 中国石油化工股份有限公司 | 一种基于岩相约束的avo反演油气检测方法 |
US20140324354A1 (en) * | 2013-04-29 | 2014-10-30 | King Fahd University Of Petroleum And Minerals | Transmission coefficient method for avo seismic analysis |
CN104297800A (zh) * | 2014-10-09 | 2015-01-21 | 中国海洋石油总公司 | 一种自相控叠前反演方法 |
-
2018
- 2018-08-23 CN CN201810966875.5A patent/CN110857997A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675907A (zh) * | 2012-09-20 | 2014-03-26 | 中国石油化工股份有限公司 | 一种基于岩相约束的avo反演油气检测方法 |
US20140324354A1 (en) * | 2013-04-29 | 2014-10-30 | King Fahd University Of Petroleum And Minerals | Transmission coefficient method for avo seismic analysis |
CN103293552A (zh) * | 2013-05-24 | 2013-09-11 | 中国石油天然气集团公司 | 一种叠前地震资料的反演方法及系统 |
CN104297800A (zh) * | 2014-10-09 | 2015-01-21 | 中国海洋石油总公司 | 一种自相控叠前反演方法 |
Non-Patent Citations (3)
Title |
---|
郑笑雪,等: "横向约束分步叠前弹性参数反演", 《石油地球物理勘探》 * |
闫惠中: "多波联合AVO反演及流体预测", 《中国优秀硕士学位论文全文数据库》 * |
霍国栋,等: "纵向和横向同时约束AVO反演", 《地球物理学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113900140A (zh) * | 2021-09-30 | 2022-01-07 | 中国石油大学(北京) | 一种基于时空组合的地震数据优化的方法及装置 |
CN113900140B (zh) * | 2021-09-30 | 2022-08-23 | 中国石油大学(北京) | 一种基于时空组合的地震数据优化的方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106842313B (zh) | 基于方位叠前地震数据的各向异性参数反演方法 | |
CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
US6014342A (en) | Method of evaluating a subsurface region using gather sensitive data discrimination | |
US7768870B2 (en) | Method for adjusting a seismic wave velocity model according to information recorded in wells | |
US5696735A (en) | Seismic migration using offset checkshot data | |
CN103733089B (zh) | 用于包括不确定性估计的地下表征的系统和方法 | |
AU2004220765B2 (en) | Method for stable estimation of anisotropic parameters for P-wave prestack imaging | |
CN103097914A (zh) | 用于表征储集层的演变的方法 | |
CN102692645A (zh) | 利用纵波、转换波数据联合反演储层纵横波速度比的方法 | |
Guo et al. | Becoming effective velocity-model builders and depth imagers, Part 2—The basics of velocity-model building, examples and discussions | |
CN110187389A (zh) | 一种基于薄层反射理论的ava反演方法 | |
CN110857997A (zh) | 一种基于横向约束的分步叠前弹性参数反演方法及系统 | |
CN109324344A (zh) | 基于纯纵波和拟声波反演的页岩厚度预测方法及系统 | |
CN106597547A (zh) | 薄储层地震准确描述方法 | |
CN107238859A (zh) | 计算虚拟井合成地震记录的方法及系统 | |
CN112684498A (zh) | 一种基于宽方位地震数据的储层裂缝预测方法及系统 | |
CN105445795A (zh) | 倾斜界面的速度分析方法及其系统 | |
CN109581521B (zh) | Tti各向异性的局部层析方法及系统 | |
CN104965226A (zh) | 一种岩体中油气信息的分析方法 | |
CN113589375A (zh) | 基于倾斜层约束旅行时计算的vsp层速度反演方法 | |
CN110618448B (zh) | 多波动态匹配方法及系统 | |
CN113075734B (zh) | 一种基于信噪比约束的剩余曲率谱计算方法及装置 | |
CN105467448B (zh) | 高速层覆盖下基于初至拾取的近地表参数反演方法 | |
CN113009579B (zh) | 地震数据反演方法及装置 | |
CN110850504B (zh) | 一种基于铀曲线拟阻抗约束的页岩密度参数叠前反演方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200303 |
|
RJ01 | Rejection of invention patent application after publication |