CN104570072A - 一种粘弹性介质中的球面pp波反射系数建模方法 - Google Patents
一种粘弹性介质中的球面pp波反射系数建模方法 Download PDFInfo
- Publication number
- CN104570072A CN104570072A CN201310483343.3A CN201310483343A CN104570072A CN 104570072 A CN104570072 A CN 104570072A CN 201310483343 A CN201310483343 A CN 201310483343A CN 104570072 A CN104570072 A CN 104570072A
- Authority
- CN
- China
- Prior art keywords
- wave
- reflection coefficient
- spherical
- formula
- sphere
- 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
Landscapes
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种粘弹性介质中的球面PP波反射系数建模方法,属于石油地球物理勘探领域。本方法包括:(1)基于White模型计算纵波相速度υp和品质因子Q-1:从储层参数出发,基于White模型计算得到纵波相速度υp和品质因子Q-1;(2)求取频散孔隙介质中平面PP波反射系数;将第(1)步得到的纵波相速度υp和品质因子Q-1代入到频散介质的Zoeppritz方程中,计算得到频散孔隙介质中的平面PP波的反射系数(3)求取球面波反射系数:在第(2)步中求得平面PP波反射系数之后,根据球面波的平面波分解算法对频散孔隙介质中的球面PP波反射系数进行建模,求得频散孔隙介质中的球面PP波反射系数
Description
技术领域
本发明属于石油地球物理勘探领域,具体涉及一种粘弹性介质中的球面PP波反射系数建模方法,从储层参数出发,基于White模型和球面波的平面波分解算法计算粘弹性介质中球面PP波反射系数的建模。
背景技术
传统的AVO(振幅随偏移距)或AVA(振幅随入射角变化)是基于平面波(几何地震学)的假设。然而,当震源或检波器离反射界面很近或者需要研究临界角附近的反射波场的时候,几何地震学近似就变得不准确。这是由于实际地震勘探中使用的是点源,激发的是球面波。(1961)最早提出了球面波反射系数的概念,并通过解析解的结果分析了球面波反射波场(1959;1960)以及临界角附近的反射波和首波的动力学特征(1962;1965)。基于球面波的平面波分解算法,Haase(2004)给出了层状均匀介质中球面波反射系数的计算方法,并讨论了第1类和第3类AVO对应的球面PP和PS波的反射系数振幅和相位特点。Ursenbach(2007)通过引入一种特殊形式的子波,使得在计算球面波反射系数过程中对频率的积分可以解析求得,大大提高了计算效率。Ayzenberg等人(2007,2009)给出了利用有效反射系数(ERCs)表示的球面波反射系数形式。
传统AVO(AVA)技术的另一个假设就是地下介质完全弹性假设。而VSP记录(Hauge,1981;Hedl in等,2001)、测井数据(Schmitt,1999)以及实验室中的岩石物理测量(Behura等,2007;Winkler和Nur,1982)均显示地震波在传播过程中会出现衰减和频散现象,尤其是对于含有碳氢化合物的地层,地震波的衰减更为明显。根据地震波传播理论,国外部分学者将其归因于地下介质强吸收的反射系数。White(1975),De Hoop(1991),Ursin(2002)等在关于如何计算吸收反射系数方面给出了一系列理论上的讨论,Innanen和Weglein(2007)将一般的吸收逆散射问题简化为反演单一界面处的反射系数,Innanen(2011)详细推导了地震波由弹性非衰减介质入射到吸收介质时反射及透射系数的表达式。实际上,完整描述或反演地下介质中吸收的反射系数,既需要考虑振幅随角度的变化(AVA),同时也要考虑振幅随频率的变化关系(AVF)。
目前,国际上尚没有综合考虑了以上两种情况的关于粘弹性介质中球面波反射系数特征的研究。
发明内容
本发明的目的在于解决上述现有技术中存在的难题,提供一种粘弹性介质中的球面PP波反射系数建模方法,以实现粘弹性介质中球面波反射系数建模为目标,将White模型和球面波的平面波分解算法有效结合,实现了粘弹性介质中球面PP波反射系数建模。该建模方法物理意义明确且更加忠于实际地震勘探情况,对于更加精确地描述含油气储层有一定的指导意义。
本发明是通过以下技术方案实现的:
一种粘弹性介质中的球面PP波反射系数建模方法,包括:
(1)基于White模型计算纵波相速度υp和品质因子Q-1:
从储层参数出发,基于White模型计算得到纵波相速度υp和品质因子Q-1;所述储层参数包括岩石骨架特性参数和孔隙流体特性参数;
(2)求取频散孔隙介质中平面PP波反射系数;
将第(1)步得到的纵波相速度υp和品质因子Q-1代入到频散介质的Zoeppritz方程中,计算得到频散孔隙介质中的平面PP波的反射系数
(3)求取球面波反射系数;
在第(2)步中求得平面PP波反射系数之后,根据球面波的平面波分解算法对频散孔隙介质中的球面PP波反射系数进行建模,求得频散孔隙介质中的球面PP波反射系数
所述步骤(1)是利用公式(12)求得纵波相速度υp,利用公式(13)求得品质因子Q-1的:
其中,Re和Im分别代表实部和虚部;
υ为复速度,根据关系求得;
其中为平均密度,ρ1和ρ2分别为含气和含水孔隙介质的密度,p1和p2分别为权重因子;p1=d1/(d1+d2),p2=d2/(d1+d2),每个周期内包含两个孔隙介质层1和2,对应小标1和2,层厚为dl,l=1,2;
其中,
其中,
式中μdry为干岩石的剪切模量,KG为Gassmann体积模量:
KG=Kdry+α2M, (4)
其中
α为Biot系数,Ks是岩石矿物体积模量,Kf是流体体积模量,Kdry是干岩石体积模量,φ是孔隙度。
其中η为流体的粘滞性系数,κ为岩石骨架的渗透率,
所述步骤(2)是这样实现的:
所述频散介质的Zoeppritz方程如下:
其中c是速度,θ是角度,脚标p和s分别代表P波和S波,脚标1和2分别代表泥岩盖层和砂岩储层。依赖于频率的P波速度可以表示为:
V为相速度,通过公式(12)求得,α为吸收系数,其表达式为:
将步骤(1)得到的纵波相速度υp和品质因子Q-1代入公式(15)和公式(16),然后求解方程(14),得到频散孔隙介质中的平面PP波的反射系数是入射角θp1和频率ω的函数。
所述步骤(3)中所述根据球面波的平面波分解算法对频散孔隙介质中的球面PP波反射系数进行建模是通过下式实现的:
其中,
其中,S=Vp1/(Rω),Vp1为上层介质中的P波速度(Vp是已知的),R=(z+h)/cosθi,其中θi是球面波的入射角,θ是平面波的入射角,J1代表一阶Bessel函数,i是虚数单位;
h为震源到反射界面的深度,z为检波器到反射界面的深度;
通过求解公式(19)即得到频散孔隙介质中的球面PP波反射系数公式(19)中等号右边的所有参数都是已知的。
与现有技术相比,本发明的有益效果是:本发明相对于已有技术最大的创新点就是引入了地层的粘弹性信息。不仅考虑了地震波在地下传播过程中的球面波效应,而且考虑了含油气储层的频散和衰减效应,能够更加精确地反映储层的物理特性。
附图说明
图1 储层模型示意图。
图2 震源和接收点间界面反射射线的几何示意图。
图3-1 P波相速度随频率的变化关系。
图3-2 品质因子随频率的变化关系。
图4-1 平面PP波反射系数振幅随入射角和频率变化。
图4-2 平面PP波反射系数相位随入射角和频率变化。
图5-1 球面PP波反射系数振幅随入射角和频率变化。
图5-2 球面PP波反射系数相位随入射角和频率变化。
图6 本发明方法的步骤框图。
具体实施方式
下面结合附图对本发明作进一步详细描述:
本发明将White模型和球面波的平面波分解算法进行结合,综合考虑了粘弹性介质的衰减和频散效应以及地震波传播的球面波效应。首先根据White模型和粘弹性介质的Zeoppritz方程得到频散孔隙介质中的平面波反射系数,然后基于球面波的平面波分解算法得到频散孔隙介质中球面波反射系数的AVAF(振幅随频率和角度变化)特征。
如图6所示,本方法具体包括:
(1)基于White模型计算纵波相速度υp和品质因子Q-1
从储层参数出发,基于White模型的建模方法,可以得到纵波相速度υp和品质因子Q-1。储层参数包括岩石骨架特性参数和孔隙流体特性参数,所述岩石骨架特性参数包括岩石矿物体积模量Ks,干岩石体积模量Kdry,干岩石的剪切模量μdry,孔隙度φ,岩石骨架渗透率κ,所述孔隙流体特性参数包括流体体积模量Kf,流体粘滞性系数η,流体密度ρ。
White模型的具体建模方法如下:
考虑周期性薄互层,每个周期内包含两个孔隙介质层1和2(该方法的使用不受孔隙介质层数的限制),孔隙内分别含气和含水,层厚为dl,l=1,2(如图1所示)。根据White模型,Carcione和Picotti给出该模型的P波复模量表达式为:
其中,
式(2)中p1=d1/(d1+d2),l=1,2,忽略脚标,对于每层介质可以得到:
式中μdry为干岩石的剪切模量,KG为Gassmann体积模量:
KG=Kdry+α2M, (4)
其中
α也称Biot系数,Ks是岩石矿物体积模量,Kf是流体体积模量,Kdry是干岩石体积模量,φ是孔隙度。令
为快P波流体张力与总正应力的比值,
是与慢P波相关的阻抗,其中η为流体的粘滞性系数,κ为岩石骨架的渗透率,
是慢P波的复波数,
是有效模量,其中
若令ρ为体积密度,则根据关系可得到复速度υ,其中(P是权重因子,具体意义在公式(2)之后给予了说明,ρ1和ρ2分别为含气和含水孔隙介质的密度)是平均密度。这样,相速度和品质因子可表示为:
Re和Im分别代表实部和虚部。
(2)频散孔隙介质中平面PP波反射系数求取;
将第(1)步得到的纵波相速度υp和品质因子Q-1代入到频散介质的Zoeppritz方程中,可以得到频散孔隙介质中的平面PP波的反射系数
在White模型中,剪切波的吸收特性不受含气性影响,因此在大多数情况下可以忽略。由Hooke定律和界面处的连续性条件,可以得到频散介质对应的Zoeppritz方程:
其中c是速度,θ是角度,脚标p和s分别代表P波和S波,脚标1和2分别代表泥岩盖层和砂岩储层。依赖于频率的P波速度可以表示为:
V为相速度(就是公式(12)中的υp:),α为吸收系数,其表达式为
P波的相速度和品质因子可以通过式(12)和式(13)计算得到。
通过求解方程(14),可以得到频散孔隙介质中的平面PP波反射系数是入射角θp1和频率ω的函数。
(3)基于平面波分解算法的球面波反射系数求取;
在第(2)步中求得频散孔隙介质中平面PP波反射系数之后,根据球面波的平面波分解算法,最终就可以建模得到频散孔隙介质中的球面PP波反射系数具体方法如下:
根据球面波的平面波分解理论,反射球面谐波的位移势函数可以表示为(Haase,2004):
其中r,z,h的物理意义如图2所示,r为偏移距大小,h为震源到反射界面的深度,z为检波器到反射界面的深度,ω为球面谐波的角频率,p和ξ分别为上层介质中的水平和垂直慢度,Jo为零阶贝塞尔函数,Rpp为平面波的反射系数。在粘弹性介质中,平面波的反射系数Rpp不仅与入射角度有关,还与频率有关,所以将通过解方程(14)得到的代入式(17):
式(18)即为粘弹性介质中球面谐波势函数的表达式。为得到球面谐波的反射系数形式,对式(18)做如下处理:
①在检波点位置处求位移势函数在波传播方向的偏导数,得到位移函数;
②用均匀介质中的球面波场对反射波场校正,消除球面扩散对反射波振幅的影响;
③将积分变量由射线参数p变为cosθ。
最终得到频散孔隙介质中球面谐波反射系数的形式为:
其中,
上式中S=Vp1/(Rω),Vp1为上层介质中的P波速度,R=(z+h)/cosθi。
(19)式就是本发明方法最后得到的表达式,通过(19)式就完成了本发明的建模。
下面通过一个理论储层模型来说明本发明的效果。模型为具有低阻抗的浅层疏松含气砂岩储层,上覆有泥岩盖层。泥岩盖层的特性参数为:Vp=2190m/s,Vs=820m/s,ρ=2.16g/cm3。砂岩储层的特性参数为:Kdry=1.56GPa,μdry=1.10GPa,φ=0.33,K=2.00darcy,ρg=0.15g/cm3,ηg=0.01cP,Kw=2.42GPa,ρw=1.00g/cm3,ηw=1.00cP,d1=2m,d2=0.5m,h=500m。砂岩岩石颗粒的体积模量和密度分别为Ks=38GPa和ρs=2.65g/cm3。
图3-1和图3-2分别为由White模型得到的砂岩储层的P波相速度和品质因子随频率的变化关系。图4-1和图4-2分别为频散孔隙介质中平面PP波反射系数AVAF特征,即砂岩储层定界面对应的平面PP波反射系数(振幅和相位)随入射角和频率的变化关系。图5-1和图5-2为使用本发明的方法得到的砂岩储层定界面对应的球面PP波反射系数(振幅和相位)随入射角和频率的变化关系,即频散孔隙介质中球面PP波反射系数AVAF特征。
AVO技术作为一种较为成熟的烃类直接检测手段,在实际生产应用中已经取得了良好的效果。传统的AVO技术是基于平面波理论建立起来的,并总是假设地下介质是完全弹性的。然而,地震勘探中使用的是爆炸源,得到的是球面波记录,当震源或检波器离反射界面很近或者需要研究临界角附近的反射波场的时候,平面波近似就变得不准确。而且VSP资料、测井资料以及实验室中的岩石物理测量均显示地震波在实际传播过程中会发生衰减以及速度频散现象,尤其对于含有碳氢化合物的区域,衰减非常明显。所以,忽略地震波在地下传播的球面波效应及地下介质的频散和衰减效应会对AVO分析及储层预测带来巨大的风险。
本发明将White模型和球面波的平面波分解算法有效结合,综合考虑了地震波传播的球面波效应及介质的频散和衰减效应,得到了一种粘弹性介质中球面PP波反射系数的建模方法,并对其AVAF特征进行了研究。这项发明物理意义明确且更加忠于实际地震勘探情况,能够更加真实地反映含油气储层的物理特性,对于进行精确地储层预测工作具有一定的指导意义。
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。
Claims (4)
1.一种粘弹性介质中的球面PP波反射系数建模方法,其特征在于:所述方法包括:
(1)基于White模型计算纵波相速度υp和品质因子Q-1:
从储层参数出发,基于White模型计算得到纵波相速度υp和品质因子Q-1;所述储层参数包括岩石骨架特性参数和孔隙流体特性参数;
(2)求取频散孔隙介质中平面PP波反射系数;
将第(1)步得到的纵波相速度υp和品质因子Q-1代入到频散介质的Zoeppritz方程中,计算得到频散孔隙介质中的平面PP波的反射系数
(3)求取球面波反射系数;
在第(2)步中求得平面PP波反射系数之后,根据球面波的平面波分解算法对频散孔隙介质中的球面PP波反射系数进行建模,求得频散孔隙介质中的球面PP波反射系数
2.根据权利要求1所述的粘弹性介质中的球面PP波反射系数建模方法,其特征在于:所述步骤(1)是利用公式(12)求得纵波相速度υp,利用公式(13)求得品质因子Q-1的:
其中,Re和Im分别代表实部和虚部;
υ为复速度,根据关系求得;
其中为平均密度,ρ1和ρ2分别为含气和含水孔隙介质的密度,p1和p2分别为权重因子;p1=d1/(d1+d2),p2=d2/(d1+d2),每个周期内包含两个孔隙介质层1和2,对应小标1和2,层厚为dl,l=1,2;
其中,
其中,
式中μdry为干岩石的剪切模量,KG为Gassmann体积模量:
KG=Kdry+α2M, (4)
其中
α为Biot系数,Ks是岩石矿物体积模量,Kf是流体体积模量,Kdry是干岩石体积模量,φ是孔隙度;
其中η为流体的粘滞性系数,κ为岩石骨架的渗透率,
3.根据权利要求2所述的粘弹性介质中的球面PP波反射系数建模方法,其特征在于:所述步骤(2)是这样实现的:
所述频散介质的Zoeppritz方程如下:
其中c是速度,θ是角度,脚标p和s分别代表P波和S波,脚标1和2分别代表泥岩盖层和砂岩储层。依赖于频率的P波速度可以表示为:
V为相速度,通过公式(12)求得,α为吸收系数,其表达式为:
将步骤(1)得到的纵波相速度υp和品质因子Q-1代入公式(15)和公式(16),然后求解方程(14),得到频散孔隙介质中的平面PP波的反射系数是入射角θp1和频率ω的函数。
4.根据权利要求3所述的粘弹性介质中的球面PP波反射系数建模方法,其特征在于:所述步骤(3)中所述根据球面波的平面波分解算法对频散孔隙介质中的球面PP波反射系数进行建模是通过下式实现的:
其中,
其中,S=Vp1/(Rω),Vp1为上层介质中的P波速度,R=(z+h)/cosθi,其中θi是球面波的入射角,θ是平面波的入射角,J1代表一阶Bessel函数,i是虚数单位;
h为震源到反射界面的深度,z为检波器到反射界面的深度;
通过求解公式(19)即得到频散孔隙介质中的球面PP波反射系数
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310483343.3A CN104570072B (zh) | 2013-10-16 | 2013-10-16 | 一种粘弹性介质中的球面pp波反射系数建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310483343.3A CN104570072B (zh) | 2013-10-16 | 2013-10-16 | 一种粘弹性介质中的球面pp波反射系数建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104570072A true CN104570072A (zh) | 2015-04-29 |
CN104570072B CN104570072B (zh) | 2017-10-27 |
Family
ID=53086586
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310483343.3A Active CN104570072B (zh) | 2013-10-16 | 2013-10-16 | 一种粘弹性介质中的球面pp波反射系数建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104570072B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104931363A (zh) * | 2015-06-23 | 2015-09-23 | 江西理工大学 | 节理岩体变形模量测试方法 |
CN104950332A (zh) * | 2015-06-18 | 2015-09-30 | 河海大学 | 一种弹性多层介质中平面波反射系数的计算方法 |
CN106066490A (zh) * | 2016-05-25 | 2016-11-02 | 中国石油大学(北京) | 基于球面波的叠前密度反演方法及装置 |
CN106249281A (zh) * | 2015-06-10 | 2016-12-21 | 中国石油化工股份有限公司 | 叠前地震的粘弹性介质振幅随频率变化的反演方法和系统 |
CN106908839A (zh) * | 2017-05-05 | 2017-06-30 | 中国石油大学(华东) | 一种地震波精确反射系数解析方法及装置 |
CN107479090A (zh) * | 2016-06-08 | 2017-12-15 | 中国石油化工股份有限公司 | 计算含定向裂隙孔隙介质的频变avo响应的方法及系统 |
CN108181654A (zh) * | 2018-01-18 | 2018-06-19 | 中国石油大学(北京) | 基于多尺度岩石物理模型的avaf模拟方法及装置 |
CN110161562A (zh) * | 2019-04-29 | 2019-08-23 | 中国石油大学(华东) | 倾斜横向各向同性介质中的裂缝弱度反演方法及系统 |
CN110954959A (zh) * | 2018-09-27 | 2020-04-03 | 中国石油化工股份有限公司 | 球面透射波特性分析方法及计算机可读存储介质 |
CN111208312A (zh) * | 2019-12-30 | 2020-05-29 | 浙江大学 | 基于波的高速旋转管状结构的转速及应变监测方法 |
CN112130207A (zh) * | 2020-09-25 | 2020-12-25 | 中国科学院武汉岩土力学研究所 | 一种基于球形装药条件下由地面振动计算地下振动的方法 |
CN112505751A (zh) * | 2019-09-16 | 2021-03-16 | 中国石油化工股份有限公司 | 球面波ps反射系数计算方法及系统 |
CN112649845A (zh) * | 2019-10-10 | 2021-04-13 | 中国石油化工股份有限公司 | 球面波反射系数计算方法及系统 |
CN113805232A (zh) * | 2020-06-17 | 2021-12-17 | 中国石油化工股份有限公司 | 浅层地表的品质因子的估计方法、系统及存储介质 |
CN116879952A (zh) * | 2023-09-06 | 2023-10-13 | 自然资源部第一海洋研究所 | 点源弹性波海底反射系数的计算方法、存储介质和设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004088758A2 (fr) * | 2003-03-28 | 2004-10-14 | Centre National De La Recherche Scientifique- Cnrs | Procede et dispositif de mesure acoustique de caracteristique physique de materiaux poroelastique |
CN1987522A (zh) * | 2005-12-22 | 2007-06-27 | 陈信平 | 一种直接探测地下石油、天然气和煤层气的方法 |
CN102478667A (zh) * | 2010-11-30 | 2012-05-30 | 中国石油天然气集团公司 | 频散avo反演方法 |
CN103149586A (zh) * | 2013-02-04 | 2013-06-12 | 西安交通大学 | 一种倾斜层状粘弹性介质中波场正演模拟方法 |
-
2013
- 2013-10-16 CN CN201310483343.3A patent/CN104570072B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004088758A2 (fr) * | 2003-03-28 | 2004-10-14 | Centre National De La Recherche Scientifique- Cnrs | Procede et dispositif de mesure acoustique de caracteristique physique de materiaux poroelastique |
CN1987522A (zh) * | 2005-12-22 | 2007-06-27 | 陈信平 | 一种直接探测地下石油、天然气和煤层气的方法 |
CN102478667A (zh) * | 2010-11-30 | 2012-05-30 | 中国石油天然气集团公司 | 频散avo反演方法 |
CN103149586A (zh) * | 2013-02-04 | 2013-06-12 | 西安交通大学 | 一种倾斜层状粘弹性介质中波场正演模拟方法 |
Non-Patent Citations (1)
Title |
---|
王艳香等: "角道集波形校正方法", 《石油地球物理勘探》 * |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106249281A (zh) * | 2015-06-10 | 2016-12-21 | 中国石油化工股份有限公司 | 叠前地震的粘弹性介质振幅随频率变化的反演方法和系统 |
CN104950332B (zh) * | 2015-06-18 | 2018-02-02 | 河海大学 | 一种弹性多层介质中平面波反射系数的计算方法 |
CN104950332A (zh) * | 2015-06-18 | 2015-09-30 | 河海大学 | 一种弹性多层介质中平面波反射系数的计算方法 |
CN104931363A (zh) * | 2015-06-23 | 2015-09-23 | 江西理工大学 | 节理岩体变形模量测试方法 |
CN104931363B (zh) * | 2015-06-23 | 2017-06-30 | 江西理工大学 | 节理岩体变形模量测试方法 |
CN106066490A (zh) * | 2016-05-25 | 2016-11-02 | 中国石油大学(北京) | 基于球面波的叠前密度反演方法及装置 |
CN107479090B (zh) * | 2016-06-08 | 2019-11-12 | 中国石油化工股份有限公司 | 计算含定向裂隙孔隙介质的频变avo响应的方法及系统 |
CN107479090A (zh) * | 2016-06-08 | 2017-12-15 | 中国石油化工股份有限公司 | 计算含定向裂隙孔隙介质的频变avo响应的方法及系统 |
CN106908839B (zh) * | 2017-05-05 | 2018-09-28 | 中国石油大学(华东) | 一种地震波精确反射系数解析方法及装置 |
CN106908839A (zh) * | 2017-05-05 | 2017-06-30 | 中国石油大学(华东) | 一种地震波精确反射系数解析方法及装置 |
CN108181654A (zh) * | 2018-01-18 | 2018-06-19 | 中国石油大学(北京) | 基于多尺度岩石物理模型的avaf模拟方法及装置 |
CN108181654B (zh) * | 2018-01-18 | 2019-10-11 | 中国石油大学(北京) | 基于多尺度岩石物理模型的avaf模拟方法及装置 |
CN110954959A (zh) * | 2018-09-27 | 2020-04-03 | 中国石油化工股份有限公司 | 球面透射波特性分析方法及计算机可读存储介质 |
CN110954959B (zh) * | 2018-09-27 | 2021-11-05 | 中国石油化工股份有限公司 | 球面透射波特性分析方法及计算机可读存储介质 |
CN110161562A (zh) * | 2019-04-29 | 2019-08-23 | 中国石油大学(华东) | 倾斜横向各向同性介质中的裂缝弱度反演方法及系统 |
CN110161562B (zh) * | 2019-04-29 | 2020-09-18 | 中南大学 | 倾斜横向各向同性介质中的裂缝弱度反演方法及系统 |
CN112505751A (zh) * | 2019-09-16 | 2021-03-16 | 中国石油化工股份有限公司 | 球面波ps反射系数计算方法及系统 |
CN112505751B (zh) * | 2019-09-16 | 2023-03-21 | 中国石油化工股份有限公司 | 球面波ps反射系数计算方法及系统 |
CN112649845A (zh) * | 2019-10-10 | 2021-04-13 | 中国石油化工股份有限公司 | 球面波反射系数计算方法及系统 |
CN111208312A (zh) * | 2019-12-30 | 2020-05-29 | 浙江大学 | 基于波的高速旋转管状结构的转速及应变监测方法 |
CN113805232A (zh) * | 2020-06-17 | 2021-12-17 | 中国石油化工股份有限公司 | 浅层地表的品质因子的估计方法、系统及存储介质 |
CN113805232B (zh) * | 2020-06-17 | 2024-04-09 | 中国石油化工股份有限公司 | 浅层地表的品质因子的估计方法、系统及存储介质 |
CN112130207A (zh) * | 2020-09-25 | 2020-12-25 | 中国科学院武汉岩土力学研究所 | 一种基于球形装药条件下由地面振动计算地下振动的方法 |
CN112130207B (zh) * | 2020-09-25 | 2021-07-20 | 中国科学院武汉岩土力学研究所 | 一种基于球形装药条件下由地面振动计算地下振动的方法 |
CN116879952A (zh) * | 2023-09-06 | 2023-10-13 | 自然资源部第一海洋研究所 | 点源弹性波海底反射系数的计算方法、存储介质和设备 |
CN116879952B (zh) * | 2023-09-06 | 2023-11-21 | 自然资源部第一海洋研究所 | 点源弹性波海底反射系数的计算方法、存储介质和设备 |
US11966001B1 (en) | 2023-09-06 | 2024-04-23 | First Institute Of Oceanography, Mnr | Calculation method, storage medium and device for seabed reflection coefficient of point source elastic wave |
Also Published As
Publication number | Publication date |
---|---|
CN104570072B (zh) | 2017-10-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104570072A (zh) | 一种粘弹性介质中的球面pp波反射系数建模方法 | |
CN102478667B (zh) | 频散avo反演方法 | |
CN102466816B (zh) | 一种叠前地震数据地层弹性常数参数反演的方法 | |
CN105068117B (zh) | 用于裂缝性介质的avo反演方法、装置和设备 | |
Barbosa et al. | Estimation of fracture compliance from attenuation and velocity analysis of full‐waveform sonic log data | |
US8767508B2 (en) | Using seismic P and S arrivals to determine shallow velocity structure | |
Zheng et al. | Seismic characterization of fractured reservoirs by focusing Gaussian beams | |
CN102033242A (zh) | 一种深层倾斜裂缝储层地震振幅预测方法 | |
CN104375188A (zh) | 一种地震波透射衰减补偿方法及装置 | |
CN104316966B (zh) | 一种流体识别方法及系统 | |
CN104155693A (zh) | 储层流体流度的角道集地震响应数值计算方法 | |
Mahmoudian et al. | Estimation of elastic stiffness coefficients of an orthorhombic physical model using group velocity analysis on transmission data | |
CN103513277A (zh) | 一种地震地层裂隙裂缝密度反演方法及系统 | |
CN104459777A (zh) | 基于流体体积模量avo反演的流体识别的方法及系统 | |
CN103487831A (zh) | Avo地震正演计算方法 | |
Li et al. | Three‐dimensional passive‐source reverse‐time migration of converted waves: The method | |
CN104316965A (zh) | 一种裂缝方位和强度的预测方法及系统 | |
CN103576197B (zh) | 一种转换波角道集抽取方法 | |
Barley et al. | S to P scattering at the 650 km discontinuity | |
CN102230973A (zh) | 一种三维分步傅立叶粘滞声波深度偏移方法 | |
CN104007466B (zh) | 一种基于纵波振幅的无井约束叠前反演实现的储层与流体预测方法 | |
CN103558637B (zh) | 基于三分量传感器的远探测方法 | |
CN103364834A (zh) | 一种利用叠前地震频散分析预测储层渗透率的方法 | |
CN106353807A (zh) | 裂缝识别方法和装置 | |
CN103869366A (zh) | 一种确定裂隙裂缝走向的方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |