CN106842310A - 叠前地震四参数同步反演方法 - Google Patents
叠前地震四参数同步反演方法 Download PDFInfo
- Publication number
- CN106842310A CN106842310A CN201510883629.XA CN201510883629A CN106842310A CN 106842310 A CN106842310 A CN 106842310A CN 201510883629 A CN201510883629 A CN 201510883629A CN 106842310 A CN106842310 A CN 106842310A
- Authority
- CN
- China
- Prior art keywords
- parameter
- matrix
- wave
- velocity
- density
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 239000012530 fluid Substances 0.000 claims abstract description 13
- 238000010521 absorption reaction Methods 0.000 claims abstract description 10
- 230000006641 stabilisation Effects 0.000 claims abstract description 5
- 238000011105 stabilization Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 33
- 230000005540 biological transmission Effects 0.000 claims description 12
- 238000011156 evaluation Methods 0.000 claims description 4
- 230000015556 catabolic process Effects 0.000 claims description 3
- 238000000205 computational method Methods 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- AYEKOFBPNLCAJY-UHFFFAOYSA-O thiamine pyrophosphate Chemical compound CC1=C(CCOP(O)(=O)OP(O)(O)=O)SC=[N+]1CC1=CN=C(C)N=C1N AYEKOFBPNLCAJY-UHFFFAOYSA-O 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 abstract description 11
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000009671 shengli Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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/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/624—Reservoir parameters
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
本发明公开了一种叠前地震四参数同步反演方法。该方法在粘弹介质理论基础下,建立了地震波反射系数与地层纵波速度、横波、密度和地层吸收参数Q间的直接关系,通过稳定的叠前地震反演算法实现地层纵波速度、横波、密度及地层吸收参数Q四参数同步反演。这对于利用地层吸收衰减特征识别储层的含油气性,刻画储层含流体性和空间变化,降低勘探风险具有十分重要的意义。
Description
技术领域
该技术属于油气勘探地震资料处理技术领域,特别是一种叠前地震四参数同步反演方法。
现有技术
叠前地震反演是获取地层参数、进行有利储层评价的主要方法。目前常用的叠前地震反演方法有AVO/AVA反演和EI反演,这两种反演方法的理论基础是Zoeppritz方程,所关注的是地震波在弹性介质分界面上的反射系数。常规叠前地震反演建立了地震纵波反射系数与地层纵波速度、横波速度、密度等三个弹性参数之间的关系,利用该关系式,通过叠前反演方法即可同时反演这三个参数,因此,常规叠前反演技术反演得到的是纵波速度、横波速度、密度三个弹性参数。然而,地下介质并非完全弹性,地震波在实际介质中传播时,会因岩石内部流体及颗粒的相对运动及地层本身的非弹性性质产生衰减,使用弹性介质理论不能充分描述地震波在实际地层中的传播特征。
发明目的
本发明的目的是针对现有技术存在的问题,提出一项新的叠前地震四参数同步反演方法。该方法在粘弹介质理论基础下,建立了地震波反射系数与地层纵波速度、横波、密度和地层吸收参数Q间的直接关系,通过稳定的叠前地震反演算法实现地层纵波速度、横波、密度及地层吸收参数Q四参数同步反演。这对于利用地层吸收衰减特征识别储层的含油气性,刻画储层含流体性和空间变化,降低勘探风险具有十分重要的意义。
发明内容
一种叠前地震四参数同步反演方法,首先建立地震纵波反射系数与地层纵波速度、横波速度、密度和品质因子Q四个参数的理论关系,然后基于贝叶斯方法,利用叠前不同角度道集资料的联合,同时反演这四个参数。
上述方案进一步包括:
建立地震纵波反射系数与地层纵波速度、横波速度、密度和品质因子Q四个参数的理论关系中粘弹性地震纵波反射系数近似公式为:
所述同步反演四参数是结合贝叶斯理论,建立了稳定的叠前地震四参数同步反演目标函数:
从测井数据获得低频信息,根据低频信息约束的计算方法,利用测井数据得到低频约束表示为:
其中,I0和IN分别是第1个点和第N个点的阻抗信息或者表示纵横波速度、密度及Q,Ri是第i个界面的反射系数,用矩阵表示为:
将低频约束考虑到目标函数中,将目标函数求极值,得到反射系数所满足的方程,通过求解方程,得到反演的纵波速度、横波速度、密度及品质因子Q参数。
粘弹性地震纵波反射系数近似公式推导过程为:
根据介质分解理论,对粘弹性Zoeppritz方程取近似,粘弹性介质精确的Zoeppritz方程整理为下式的形式:
MR=N (1)
其中R=(RPP,RPS,TPP,TPS)T是带吸收的反射系数和透射系数组成的向量,M和N是与弹性介质参数、粘弹性介质参数和入射角有关的系数矩阵和向量;
M中的分量用介质的参数和入射角表示出来,M11为例:
根据散射理论,在弱粘弹性近似和相似介质近似下,将粘弹性介质当作背景,将弹性参数的变化和粘弹性参数均看成扰动,因此可将矩阵M,R,N分解为背景矩阵和扰动矩阵;其中矩阵M分解为:
M=Mu+ΔM (4)
其中背景矩阵Mu为(1)式所示的矩阵。
向量N分解为:
N=Nu+ΔN (5)
其中背景矩阵Nu为:
Nu=[N11 N21 N31 N41]T (6)
扰动向量ΔN为:
ΔN=[ΔN11 ΔN21 ΔN31 ΔN41]T (7)
粘弹性介质反射透射系数向量R可分解为:
R=Ru+ΔR (8)
其中Ru表示粘弹性介质背景反射透射系数,ΔR表示粘弹性参数和扰动反射透射系数;
此时粘弹性介质弹性波拟Zoeppritz方程变为:
(Mu+ΔM)(Ru+ΔR)=(Nu+ΔN) (9)
背景矩阵满足MuRu=Nu
利用Mu可求得其逆(Mu)-1,
经过整理可得粘弹性介质弹性参数和扰动反射系数矩阵:
所以反射系数RPP为:
横波品质因子Qs与饱和度无关,它是完全独立于孔隙流体的,因此在化简的过程中不考虑横波品质因子,同时在该近似中,把看作是小量,在实际资料中Qp的取值是几到几百,的取值是小于10-2,弹性参数的变化与的乘积也是小于10-2,因此在方程近似中假定相邻的上下两层介质中Qp的值相等,即Qp1=Qp2=Qp,最后取反射吸收的实部,得到下式:
进一步合并,简化得到
进一步忽略高阶弱项,即可将上式简化为如下形式:
前述叠前地震四参数同步反演方法,还包括应用以下三种方法提取Q值曲 线:(1)基于井旁地震道的Q值提取方法;(2)经验公式Q值估算方法;(3)基于测井资料的Q值提取方法。
利用该方法进行数据测试,在测试之前,需要进行相应的数据准备工作,包括地震资料、测井资料及与测井对应的Q曲线。其中地震数据既可以输入叠前CMP道集,也可以输入叠前部分角度叠加道集。
发明效果
叠前地震四参数同步反演方法能够同时从地震数据中直接反演得到地层纵波速度、横波速度、密度和地层吸收有关的参数Q,其中Q对含气储层具有很高的敏感性,纵横波速度和密度可以进一步获取储层流体敏感参数,进行流体识别。与基于波动方程波形反演等方法相比,本方法具有计算效率高且算法稳定、便于实现的特点。
附图说明
图1所示一维地层模型。
图2是利用单程波法正演模拟得到的叠前CMP道集。
图3是利用叠前地震四参数同步方法反演的纵波速度、横波速度、密度和地层Q值曲线。
图4是合成的Marmousi模型叠前CMP道集。
图5(a)(b)(c)分别是已知的纵波速度、横波速度和密度剖面,图5(d)是利用李氏经验公式计算的纵波品质因子剖面。
图6(a)(b)(c)(d)分别为利用本方法反演的纵波速度、横波速度、密度和地层Q值剖面。
图7是提取的第52道反演结果与原始数据的对比图。
具体实施方式
常规叠前地震反演建立了地震纵波反射系数与地层纵波速度、横波速度、密 度等三个弹性参数之间的关系,利用该关系式,通过叠前反演方法即可同时反演这三个参数。
本实施例则是首先建立地震纵波反射系数与地层纵波速度、横波速度、密度和品质因子Q等四个参数的理论关系,然后基于贝叶斯方法,利用叠前不同角度道集资料的联合,同时反演这四个参数。
(1)粘弹性地震纵波反射系数近似公式推导
根据介质分解理论,对粘弹性Zoeppritz方程取近似。粘弹性介质精确的Zoeppritz方程整理为下式的形式:
MR=N (1)
其中R=(RPP,RPS,TPP,TPS)T是带吸收的反射系数和透射系数组成的向量,M和N是与弹性介质参数、粘弹性介质参数和入射角有关的系数矩阵和向量。
M中的分量都可以用介质的参数和入射角表示出来,M11为例:
根据散射理论,在弱粘弹性近似和相似介质近似下,将粘弹性介质当作背景,将弹性参数的变化和粘弹性参数均看成扰动,因此可将矩阵M,R,N分解为背景矩阵和扰动矩阵。其中矩阵M可分解为:
M=Mu+ΔM (4)
其中背景矩阵Mu为(1)式所示的矩阵。
向量N可以分解为:
N=Nu+ΔN (5)
其中背景矩阵Nu为
Nu=[N11 N21 N31 N41]T (6)
扰动向量ΔN为:
ΔN=[ΔN11 ΔN21 ΔN31 ΔN41]T (7)
粘弹性介质反射透射系数向量R可分解为:
R=Ru+ΔR (8)
其中Ru表示粘弹性介质背景反射透射系数,ΔR表示粘弹性参数和扰动反射透射系数。
此时粘弹性介质弹性波拟Zoeppritz方程变为:
(Mu+ΔM)(Ru+ΔR)=(Nu+ΔN) (9)
背景矩阵满足MuRu=Nu
利用Mu可求得其逆(Mu)-1,
经过整理可得粘弹性介质弹性参数和扰动反射系数矩阵:
所以反射系数RPP为:
横波品质因子Qs与饱和度无关,它是完全独立于孔隙流体的,因此在化简的过程中不考虑横波品质因子。同时在该近似中,把看作是小量,在实际资料中Qp的取值一般是几到几百,的取值一般是小于10-2,弹性参数的变化与的乘积一般也是小于10-2,因此在方程近似中假定相邻的上下两层介质中Qp的值相等,即Qp1=Qp2=Qp。反射系数中含有虚部项,不便于求取。最后取反射吸收的 实部,得到下式:
进一步合并,简化得到
进一步忽略高阶弱项,即可将上式简化为如下形式:
该表达式建立了地震纵波反射系数与纵波速度、横波速度、密度和品质因子等四个参数之间的理论关系。基于该表达式,即可利用叠前不同角度道集资料同时反演这四个参数。
(2)四参数同步反演框架建立
叠前地震四参数同步反演在贝叶斯框架下进行。结合先验信息和与正演模型有关的似然函数,取最大后验概率作为反演问题的解。基于贝叶斯理论,可以将先验信息与观测数据结合起来,充分考虑地震数据信噪比对反演的影响,通过均衡采用的地震数据与测井数据的信息量,将少参数反演的误差,提高反演结果的可信程度。结合贝叶斯理论,建立了稳定的叠前地震四参数同步反演目标函数。
由于地震数据是带限的,缺少低频信息,因此需要从测井数据获得低频信息,从而使反演的结果更符合地质特征。根据Yin等(2008)的低频信息约束的计算方法,我们利用测井数据得到低频约束,可以表示为:
其中,I0和IN分别是第1个点和第N个点的阻抗信息(也可以表示纵横波速度、密度及Q),Ri是第i个界面的反射系数,用矩阵可以表示为:
将低频约束考虑到目标函数中,将目标函数求极值,就可以得到反射系数所满足的方程,通过求解方程,得到反演的纵横波速度、密度及Q参数。
(3)Q曲线求取
在叠前地震四参数同步反演中,除了需要纵横波速度、密度,还需要相应的Q曲线。目前主要应用以下三种方法提取Q值曲线:(1)基于井旁地震道的Q值提取方法。(2)经验公式Q值估算方法。(3)基于测井资料的Q值提取方法。
下面结合附图对上述实施例做进一步说明。
图1所示一维地层模型,共分为三层,第二层含气,地层Q值比较低,为50;第一层和第三层Q值比较高,为500。图中从左到右分别是纵波速度、横波速度、密度和地层Q值曲线。
图2是利用单程波法正演模拟得到的叠前CMP道集,1毫秒采样,炮间距10-250m,覆盖次数25次。
图3是利用叠前地震四参数同步方法反演的纵横波速度、密度和地层Q值 曲线,其中点线是反演的结果,实线是原始数据,从图中可以看出,反演的总趋势与模型一致。以Marmousi2模型为例来验证方法的正确性和适用性。
图4是合成的叠前CMP道集,图5(a)(b)(c)分别是已知的纵波速度、横波速度和密度剖面,图5(d)是利用李氏经验公式计算的纵波品质因子剖面。模型中有一个气层和一个油层,气层位置大约在时间0.2s、CDP30~CDP70之间,油层位置大约在时间0.8s、CDP20~CDP150之间,它们表现出低纵波速度、密度和品质因子,而横波速度没有降低的特征。
图6(a)(b)(c)(d)为利用本章的方法反演的纵波速度、横波速度、密度和地层Q值剖面,从反演的结果可以看出地层特征得到了比较完整的恢复,从而说明了方法的适用性。
图7是提取的第52道反演结果与原始数据的对比图,从左到右分别是纵波速度、横波速度、密度和地层Q值的倒数,其中淡色是模型数据,深色是反演结果,从图中可以看出,反演结果的总体趋势与模型数据曲线一致。
应用本方法这些程序对胜利地区典型模型进行了成功的反演,检验了方法的正确性、有效性和稳定性。对具有较复杂构造的胜利油田地区实际资料进行了联合反演,获得了较好的储层预测结果。
Claims (4)
1.一种叠前地震四参数同步反演方法,其特征在于首先建立地震纵波反射系数与地层纵波速度、横波速度、密度和品质因子Q四个参数的理论关系,然后基于贝叶斯方法,利用叠前不同角度道集资料的联合,同时反演这四个参数。
2.根据权利要求1所述的叠前地震四参数同步反演方法,其特征在于,
建立地震纵波反射系数与地层纵波速度、横波速度、密度和品质因子Q四个参数的理论关系中粘弹性地震纵波反射系数近似公式为:
所述同步反演四参数是结合贝叶斯理论,建立了稳定的叠前地震四参数同步反演目标函数:
从测井数据获得低频信息,根据低频信息约束的计算方法,利用测井数据得到低频约束表示为:
其中,I0和IN分别是第1个点和第N个点的阻抗信息或者表示纵横波速度、密度及Q,Ri是第i个界面的反射系数,用矩阵表示为:
将低频约束考虑到目标函数中,将目标函数求极值,得到反射系数所满足的方程,通过求解方程,得到反演的纵波速度、横波速度、密度及品质因子Q参数。
3.根据权利要求2所述的叠前地震四参数同步反演方法,其特征在于粘弹性地震纵波反射系数近似公式推导过程为:
根据介质分解理论,对粘弹性Zoeppritz方程取近似,粘弹性介质精确的Zoeppritz方程整理为下式的形式:
MR=N (1)
其中R=(RPP,RPS,TPP,TPS)T是带吸收的反射系数和透射系数组成的向量,M和N是与弹性介质参数、粘弹性介质参数和入射角有关的系数矩阵和向量;
M中的分量用介质的参数和入射角表示出来,M11为例:
根据散射理论,在弱粘弹性近似和相似介质近似下,将粘弹性介质当作背景,将弹性参数的变化和粘弹性参数均看成扰动,因此可将矩阵M,R,N分解为背景矩阵和扰动矩阵;其中矩阵M分解为:
M=Mu+ΔM (4)
其中背景矩阵Mu为(1)式所示的矩阵。
向量N分解为:
N=Nu+ΔN (5)
其中背景矩阵Nu为:
Nu=[N11 N21 N31 N41]T (6)
扰动向量ΔN为:
ΔN=[ΔN11 ΔN21 ΔN31 ΔN41]T (7)
粘弹性介质反射透射系数向量R可分解为:
R=Ru+ΔR (8)
其中Ru表示粘弹性介质背景反射透射系数,ΔR表示粘弹性参数和扰动反射透射 系数;
此时粘弹性介质弹性波拟Zoeppritz方程变为:
(Mu+ΔM)(Ru+ΔR)=(Nu+ΔN) (9)
背景矩阵满足MuRu=Nu
利用Mu可求得其逆(Mu)-1,
经过整理可得粘弹性介质弹性参数和扰动反射系数矩阵:
所以反射系数RPP为:
横波品质因子Qs与饱和度无关,它是完全独立于孔隙流体的,因此在化简的过程中不考虑横波品质因子,同时在该近似中,把看作是小量,在实际资料中Qp的取值是几到几百,的取值是小于10-2,弹性参数的变化与的乘积也是小于10-2,因此在方程近似中假定相邻的上下两层介质中Qp的值相等,即Qp1=Qp2=Qp,最后取反射吸收的实部,得到下式:
进一步合并,简化得到
进一步忽略高阶弱项,即可将上式简化为如下形式:
4.根据权利要求1或2、3所述的叠前地震四参数同步反演方法,其特征在于还包括应用以下三种方法提取Q值曲线:(1)基于井旁地震道的Q值提取方法;(2)经验公式Q值估算方法;(3)基于测井资料的Q值提取方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510883629.XA CN106842310B (zh) | 2015-12-04 | 2015-12-04 | 叠前地震四参数同步反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510883629.XA CN106842310B (zh) | 2015-12-04 | 2015-12-04 | 叠前地震四参数同步反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106842310A true CN106842310A (zh) | 2017-06-13 |
CN106842310B CN106842310B (zh) | 2020-10-13 |
Family
ID=59149707
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510883629.XA Active CN106842310B (zh) | 2015-12-04 | 2015-12-04 | 叠前地震四参数同步反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106842310B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108627872A (zh) * | 2018-07-04 | 2018-10-09 | 中国石油大学(华东) | 一种油气储层密度高可靠性预测方法 |
CN110161563A (zh) * | 2019-06-12 | 2019-08-23 | 中国石油大学(华东) | 一种深度域地震流体分析方法、装置、系统及存储介质 |
CN110261897A (zh) * | 2019-04-26 | 2019-09-20 | 中国石油化工股份有限公司 | 基于组稀疏的叠前四参数反演方法 |
CN111060964A (zh) * | 2018-10-16 | 2020-04-24 | 中国石油天然气股份有限公司 | 地层弹性参数的确定方法及装置 |
CN111239811A (zh) * | 2018-11-28 | 2020-06-05 | 中国石油天然气股份有限公司 | 一种判断储层流体趋势与压实趋势的方法及装置 |
CN112147683A (zh) * | 2019-06-28 | 2020-12-29 | 中国石油化工股份有限公司 | 基于岩石物理关系约束的叠前稀疏层反演方法及系统 |
CN112505771A (zh) * | 2020-10-16 | 2021-03-16 | 中国石油大学(华东) | 一种叠前统计特征参数估计方法、装置、处理器及存储介质 |
CN112578458A (zh) * | 2020-10-19 | 2021-03-30 | 中国石油大学(华东) | 叠前弹性阻抗随机反演方法、装置、存储介质及处理器 |
CN112630840A (zh) * | 2020-10-16 | 2021-04-09 | 中国石油大学(华东) | 基于统计特征参数的随机反演方法及处理器 |
CN113311483A (zh) * | 2020-02-27 | 2021-08-27 | 中国石油天然气集团有限公司 | 基于页岩油的叠前弹性参数组合反演方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104316966A (zh) * | 2014-11-12 | 2015-01-28 | 中国石油大学(华东) | 一种流体识别方法及系统 |
-
2015
- 2015-12-04 CN CN201510883629.XA patent/CN106842310B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104316966A (zh) * | 2014-11-12 | 2015-01-28 | 中国石油大学(华东) | 一种流体识别方法及系统 |
Non-Patent Citations (1)
Title |
---|
刘浩杰等: "《地层吸收衰减反演方法》", 30 September 2015, 北京:石油工业出版社 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108627872A (zh) * | 2018-07-04 | 2018-10-09 | 中国石油大学(华东) | 一种油气储层密度高可靠性预测方法 |
CN111060964B (zh) * | 2018-10-16 | 2022-05-10 | 中国石油天然气股份有限公司 | 地层弹性参数的确定方法及装置 |
CN111060964A (zh) * | 2018-10-16 | 2020-04-24 | 中国石油天然气股份有限公司 | 地层弹性参数的确定方法及装置 |
CN111239811B (zh) * | 2018-11-28 | 2022-07-05 | 中国石油天然气股份有限公司 | 一种判断储层流体趋势与压实趋势的方法及装置 |
CN111239811A (zh) * | 2018-11-28 | 2020-06-05 | 中国石油天然气股份有限公司 | 一种判断储层流体趋势与压实趋势的方法及装置 |
CN110261897B (zh) * | 2019-04-26 | 2021-07-20 | 中国石油化工股份有限公司 | 基于组稀疏的叠前四参数反演方法 |
CN110261897A (zh) * | 2019-04-26 | 2019-09-20 | 中国石油化工股份有限公司 | 基于组稀疏的叠前四参数反演方法 |
CN110161563A (zh) * | 2019-06-12 | 2019-08-23 | 中国石油大学(华东) | 一种深度域地震流体分析方法、装置、系统及存储介质 |
CN112147683A (zh) * | 2019-06-28 | 2020-12-29 | 中国石油化工股份有限公司 | 基于岩石物理关系约束的叠前稀疏层反演方法及系统 |
CN113311483A (zh) * | 2020-02-27 | 2021-08-27 | 中国石油天然气集团有限公司 | 基于页岩油的叠前弹性参数组合反演方法及系统 |
CN113311483B (zh) * | 2020-02-27 | 2024-04-30 | 中国石油天然气集团有限公司 | 基于页岩油的叠前弹性参数组合反演方法及系统 |
CN112505771A (zh) * | 2020-10-16 | 2021-03-16 | 中国石油大学(华东) | 一种叠前统计特征参数估计方法、装置、处理器及存储介质 |
CN112630840A (zh) * | 2020-10-16 | 2021-04-09 | 中国石油大学(华东) | 基于统计特征参数的随机反演方法及处理器 |
CN112578458A (zh) * | 2020-10-19 | 2021-03-30 | 中国石油大学(华东) | 叠前弹性阻抗随机反演方法、装置、存储介质及处理器 |
Also Published As
Publication number | Publication date |
---|---|
CN106842310B (zh) | 2020-10-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106842310A (zh) | 叠前地震四参数同步反演方法 | |
CN100557464C (zh) | 基于叠前波场模拟的地震勘探层位标定方法 | |
CN102466816B (zh) | 一种叠前地震数据地层弹性常数参数反演的方法 | |
CN106443768B (zh) | 一种叠前深度域合成地震记录制作方法 | |
CN104237945B (zh) | 一种地震资料自适应高分辨处理方法 | |
Yanhu et al. | A method of seismic meme inversion and its application | |
CN103527184A (zh) | 一种白云岩储层的预测方法和系统 | |
CN105089652A (zh) | 一种拟声波曲线重构与稀疏脉冲联合反演方法 | |
CN104297800B (zh) | 一种自相控叠前反演方法 | |
CN107132573A (zh) | 一种应用子波分解重构技术识别强阻抗屏蔽下岩性油藏的方法 | |
CN105938203A (zh) | 一种储层特性的检测方法及装置 | |
CN106842326B (zh) | 无横波速度测井时砂泥互层co2地质封存时移地震正演模拟方法 | |
Leaney | Microseismic source inversion in anisotropic media | |
CN106501872A (zh) | 一种裂缝储层地应力特征的计算方法及装置 | |
CN106353807A (zh) | 裂缝识别方法和装置 | |
CN104424393A (zh) | 一种基于主成分分析的地震数据储层反射特征加强方法 | |
Tuppen | Nine-component seismic amplitude inversion: A case study in the Eagle Ford shale | |
Fernandez-Concheso | Characterizing an unconventional reservoir with conventional seismic data: A case study using seismic inversion for the Vaca Muerta Formation, Neuquen Basin, Argentina | |
Ningkai et al. | Stepped and detailed seismic prediction of shallow-thin reservoirs in Chunfeng oilfield of Junggar Basin, NW China | |
CN109239777B (zh) | 一种利用联合反演方法检测构造煤发育的方法 | |
Tavella et al. | Integrated characterization of unconventional upper Jurassic reservoir in northern Mexico | |
Huang et al. | Amplitude variation with angle and azimuth inversion to estimate fracture properties in shale-gas reservoirs | |
CN106353796A (zh) | 一种提高地面地震资料分辨率方法 | |
CN104297781A (zh) | 一种射线参数域无约束弹性参数反演方法 | |
CN104516021A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |