CN110161562A - 倾斜横向各向同性介质中的裂缝弱度反演方法及*** - Google Patents
倾斜横向各向同性介质中的裂缝弱度反演方法及*** Download PDFInfo
- Publication number
- CN110161562A CN110161562A CN201910354239.1A CN201910354239A CN110161562A CN 110161562 A CN110161562 A CN 110161562A CN 201910354239 A CN201910354239 A CN 201910354239A CN 110161562 A CN110161562 A CN 110161562A
- Authority
- CN
- China
- Prior art keywords
- weakness
- crack
- inversion
- isotropy
- titled
- 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 37
- 239000011159 matrix material Substances 0.000 claims abstract description 45
- 238000009826 distribution Methods 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 238000000926 separation method Methods 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 4
- 238000011160 research Methods 0.000 abstract description 4
- 238000004891 communication Methods 0.000 abstract description 3
- 230000002265 prevention Effects 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 11
- 239000013598 vector Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000004590 computer program Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000010287 polarization Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002250 progressing effect Effects 0.000 description 1
- 238000012360 testing method 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
-
- 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
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
本公开提供了一种倾斜横向各向同性介质中的裂缝弱度反演方法及***,使用垂直横向各向同性介质和Bond矩阵的刚度矩阵,根据裂缝弱度和倾角推导出倾斜横向各向同性介质的刚度矩阵,结合刚度扰动和散射理论,分离两个弱各向异性倾斜横向各向同性介质的弱对比界面,得到该界面下PP波反射系数的方位角,以倾角作为先验信息,利用部分入射角叠加方位地震数据实现方位角叠前反演,根据反演结果估计贝叶斯框架中的裂缝弱度,合理的估计P波和S波模量和裂缝弱度,为地震波的传播研究和地震灾害预防提供了技术支持。
Description
技术领域
本公开属于物理反演领域,具体涉及一种倾斜横向各向同性介质中的裂缝 弱度反演方法及***。
背景技术
本部分的陈述仅仅是提供了与本公开相关的背景技术信息,不必然构成在 先技术。
在具有倾角的横向各向同性介质中传播的地震波的速度随着不同的方向而 变化,水平横向各向同性(HTI)或垂直横向各向同性(VTI)可以被认为是倾 斜横向各向同性(TTI)的特殊情况。地震波在这种TTI介质中传播过程中表现 出方位角速度变化,由裂缝弱度定义的各向异性可以表征方位角速度变化。
发明内容
本公开为了解决上述问题,提出了一种倾斜横向各向同性介质中的裂缝弱 度反演方法及***,本公开在分离两个倾斜横向各向同性介质的界面情况下PP 波反射系数的方位角具有偏移参数化,以倾角作为先验信息,利用部分入射角 叠加方位地震数据实现方位角叠前反演,能够有效估计/预测贝叶斯框架中的裂 缝弱度,即能够合理的估计P波和S波模量和裂缝弱度,为地震波的传播研究 和地震灾害预防提供了技术支持。
根据一些实施例,本公开采用如下技术方案:
一种倾斜横向各向同性介质中的裂缝弱度反演方法,包括以下步骤:
使用垂直横向各向同性介质和Bond矩阵的刚度矩阵,根据裂缝弱度和倾角 推导出倾斜横向各向同性介质的刚度矩阵,结合刚度扰动和散射理论,分离两 个弱各向异性倾斜横向各向同性介质的弱对比界面,得到该界面下PP波反射系 数的方位角,以倾角作为先验信息,利用部分入射角叠加方位地震数据实现方 位角叠前反演,根据反演结果估计贝叶斯框架中的裂缝弱度。
作为进一步的限定,根据法向和切向过剩弱度以及Bond变换对弱各向异性 垂直横向各向同性介质的刚度矩阵进行积分,根据法向和切向过剩弱度以及倾 角推导倾斜横向各向同性介质的刚度矩阵。
作为更进一步的限定,考虑由平面水平弱相对界面分离的两个均匀弱各向 异性倾斜横向各向同性介质,并且通过使用P波和S波速度作为坐标轴来参数 化半空间中的每个倾斜横向各向同性介质,融合密度、法向和切向裂缝弱度参 数、以及对称轴和z轴之间的角度差和方位角,确定倾斜横向各向同性介质的 刚度矩阵。
作为进一步的限定,将散射理论与倾斜横向各向同性介质刚度矩阵的扰动 相结合,根据背景P-和S-波模量、密度、裂缝弱度和倾斜角得出线性化PP波反 射系数,分离两个弱各向异性倾斜横向各向同性介质的弱对比界面。
作为进一步的限定,在贝叶斯框架下,进行方位角的反演,得到贝叶斯框 架下的裂缝弱度。
作为进一步的限定,利用Cauchy稀疏正则化和低频信息正则化以及迭代加 权最小二乘法优化算法的约束来实现方位角的叠前反演,根据反演结果估计裂 缝弱度。
作为更进一步的限定,使用Cauchy函数作为先验分布,使用高斯函数作为 似然函数,通过Cauchy先验分布和高斯似然函数联合解决后验分布,得到分布 最大化方程,将低频信息正则化与最大化方程相结合,使用迭代重加权最小二 乘算法来求解,得到模型参数的最终结果,即反演结果。
一种倾斜横向各向同性介质中的裂缝弱度反演***,被配置为运行在服务 器或处理器上,执行以上倾斜横向各向同性介质中的裂缝弱度反演方法。
一种计算机可读存储介质,其中存储有多条指令,所述指令适于由终端设 备的处理器加载并执行所述的倾斜横向各向同性介质中的裂缝弱度反演方法。
一种终端设备,包括处理器和计算机可读存储介质,处理器用于实现各指 令;计算机可读存储介质用于存储多条指令,所述指令适于由处理器加载并执 行所述的倾斜横向各向同性介质中的裂缝弱度反演方法。
与现有技术相比,本公开的有益效果为:
为了表征倾斜的横向各向同性,提出了使用可观察的PP波方位角反射数据 的方位AVO参数化和裂缝弱度的反演。根据法向和切向过剩弱度以及Bond变换 对弱各向异性VTI介质的刚度矩阵进行积分,首先根据法向和切向过剩弱度以 及倾角推导出TTI介质的刚度矩阵。将散射理论与TTI介质刚度矩阵的扰动相 结合,然后根据背景P-和S-波模量,密度,裂缝弱度和倾斜角得出线性化PP 波反射系数。分离两个弱各向异性TTI介质的弱对比界面。在贝叶斯框架内, 利用Cauchy稀疏正则化和低频信息正则化以及IRLS优化算法的约束来实现方 位AVO反演。能够合理的估计P波和S波模量和裂缝弱度,为地震波的传播研 究和地震灾害预防提供了技术支持。
附图说明
构成本公开的一部分的说明书附图用来提供对本公开的进一步理解,本公 开的示意性实施例及其说明用于解释本公开,并不构成对本公开的不当限定。
图1是本公开的流程图;
具体实施方式:
下面结合附图与实施例对本公开作进一步说明。
应该指出,以下详细说明都是例示性的,旨在对本公开提供进一步的说明。 除非另有指明,本文使用的所有技术和科学术语具有与本公开所属技术领域的 普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图 限制根据本公开的示例性实施方式。如在这里所使用的,除非上下文另外明确 指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说 明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、 组件和/或它们的组合。
正如背景技术中所述的水平横向各向同性(HTI)或垂直横向各向同性(VTI) 可以被认为是倾斜横向各向同性(TTI)的特殊情况。地震波在这种TTI介质中 传播过程中表现出方位角速度变化。裂缝弱度使我们能够表征各向异性和裂缝 发展。
为了预测这种TTI介质中的裂缝弱度,在本实施例中,使用可观察的反射 数据提出了具有偏移(AVO)参数化和反演方法的新的方位角幅度变化。首先使 用VTI介质和Bond矩阵的刚度矩阵,根据裂缝弱度和倾角推导出TTI介质的刚 度矩阵。结合刚度扰动和散射理论,在分离两个TTI介质的界面情况下PP波反 射系数的方位AVO参数化。以倾角作为先验信息,最终利用部分入射角叠加方 位地震数据实现方位AVO反演,以估计贝叶斯框架中的裂缝弱度。对合成数据 和实际数据的测试表明,对于TTI介质的表征,可以合理地和地质地估计P波 和S波模量和裂缝弱度。
具体的,倾斜横向各向同性介质的弹性刚度矩阵的建立:
为了描述TTI介质,我们考虑由平面水平弱相对界面分离的两个均匀弱各 向异性TTI介质,并且通过使用P波和S波速度αb和βb作为坐标轴来参数化半 空间中的每个TTI介质,并且密度ρb,法向和切向裂缝弱度参数δN和δT,以及 两个角度例如θ°,即倾斜(对称轴和z轴之间的角度差)和方位角φ(接收器线 的方位角φobs和对称轴φsym之间的方位角差)。入射P波位于平面[x,z]中。
通过在[x,z]平面中旋转VTI介质的对称轴,可以获得TTI介质的弹性刚度矩 阵。我们假设θ°在顺时针旋转期间倾斜为正,在逆时针旋转期间倾斜为负。
弱向各向异性VTI介质的弹性刚度矩阵在正常和切向过剩弱度方面的线性 化表达式由下式给出:
当
和
密度ρb,法向和切向裂缝弱度参数δN和δT,角度θ°,即倾斜(对称轴和z 轴之间的角度差)方位角φ(接收器线的方位角φobs和对称轴φsym之间的方位角差)
表示VTI介质的有效弹性刚度分量;δN和δT是正的并且具有无量纲法向和 切向裂缝弱度的物理意义
EN和ET表示裂缝的非负的无量纲的法向和切向柔度,即裂缝中的柔度与非 裂缝介质中的柔度的比;ZN和ZT也是非负的并且表示由裂缝添加到非裂缝介质 的单位宽度的总法向和切向柔度的物理意义
Mb和μb是非裂缝介质的压缩(P波)和剪切(S波)模量,并且λb=Mb-2μb是第一拉梅常数。
在等式(1),(2a)和(2b)中,表示VTI介质的有效弹性刚度分量;δN和 δT是正的并且具有无量纲法向和切向裂缝弱度的物理意义,从0到1不等;EN和 ET表示裂缝的非负的无量纲的法向和切向柔度,即裂缝中的柔度与非裂缝介质 中的柔度的比率;ZN和ZT也是非负的并且表示由裂缝添加到非裂缝介质的单位 宽度的总法向和切向柔度的物理意义,其具有尺寸1/应力;和是非裂缝介质的压缩(P波)和剪切(S波)模量,并且λb=Mb-2μb是第一拉梅 常数。
因此,在执行旋转之后导出TTI介质的弹性刚度矩阵(Bond,1943),其表 示为
其中Bond矩阵围绕z轴旋转的可写为
等式(3)中TTI介质的非零刚度分量如下:
等式(5a)-(5m)是TTI介质的弹性刚度矩阵,其可用于驱动弱对比界面 分离两个弱各向异性TTI介质的刚度分量的扰动。
倾斜横向各向同性介质中的线性PP波反射系数
对于由弱对比界面分离的两个均匀弱各向异性TTI介质的情况,我们假设 两个半空间中弱各向异性介质的刚度分量cij满足一阶扰动条件,由下式给出
在等式(6)中的上标0表示背景量的情况下,符号Δ表示在背景介质中的 小扰动或与它们的值的偏差,并且模|·|可以被定义为
考虑到界面上弱对比度的条件和小的弱度,我们可以导出刚度分量中扰动 的弱对比度和弱各向异性近似,并将它们写为
Δc12≈Δλb-λb(χbcos2θ°+sin2θ°)ΔδN, (7b)
Δc22≈ΔMb-χb 2MbΔδN, (7e)
Δc23≈Δλb-λb(cos2θ°+χbsin2θ°)ΔδN, (7f)
Δc44≈Δμb-μbcos2θ°ΔδT, (7j)
Δc66≈Δμb-μbsin2θ°ΔδT, (7m)
其中符号Δ还分别表示刚度分量,弹性和弱度参数中的小扰动。
将散射理论与TTI介质刚度分量中的导出扰动相结合,Shaw和Sen (2004,2006)提出了一种推导弱各向异性介质的线性化PP波反射系数的方法, 该方法由下式给出:
其中θ代表入射角;ρ表示各向同性背景的密度,Δρ表示密度项的小扰动; ΔC表示刚度矩阵中的扰动(如公式7a-7m所示),参数ξ与η慢度和极化矢量相 关,如附录A所示。
我们根据背景P-和S-波模量,密度,法向和切向裂缝弱度以及从两个弱各 向异性TTI介质分离的弱对比界面情况下的倾角的线性PP波反射系数,如下
其中这五个系数代表模型参数的权重系数,其表达式见附录A;ΔMb/MbΔμb/μbΔρb/ρb的分母代表上层和下层背景弹性参数的平均数量;并分别表示正常 和切向裂缝弱度的扰动。
如附录A所示,当θ°=90°TTI介质的反射系数与HTI介质的反射系数完全一 致,这进一步验证了TTI介质中理论推导的正确性。
法向和切向裂缝弱度项的倾角在模型参数对PP波反射系数的影响中起关键 作用。因此,包含倾斜角θ°的反射系数方程(9)极大地有助于在TTI介质中应 用方位AVO参数化和反演。
方位AVO反演用于倾斜横向各向同性介质中的裂缝弱度
相对ΔMb/Mb,Δμb/μb,和Δρb/ρb可以被ΔMb/Mb≈Δ(lnMb),Δμb/μb≈Δ(lnμb)和Δρb/ρb≈Δ(lnρb)代替,可以用符号ln(·)代表自然对数来代替,并且符号Δ表示界面 上弹性参数的小对比,即|ΔMb/Mb|<<1,|Δμb/μb|<<1和|Δρb/ρb|<<1。此外,通过使用连续变化假设,对数弹性参数和裂缝弱度的对比可以近似为Δ(lnMb)≈d(lnMb), Δ(lnμb)≈d(lnμb)和Δ(lnρb)≈d(lnρb)。
方程(9)的PP波反射系数与地震子波进行卷积,得到正方程,由下式给 出
其中SPP代表合成PP波数据,W是小波矩阵,D是差分算子的矩阵。
等式(10)可以以矩阵形式写成
[S]MNL×1=[G]MNL×5L[m]5L×1, (11)
其中M是入射角数N是方位角数,L是样本数;S给出了可观测方位地震数据 的向量,G表示与子波矩阵相关的前向算子矩阵、差分算子矩阵和系数加权矩阵, m=[lnMb,lnμb,lnρb,δN,δT]T表示模型参数的向量,表示为
其中
[SPP(θi,φj;θ°)]L×1=[SPP(t1,θi,φj;θ°),...,SPP(tL,θi,φj;θ°)]T;
此外,考虑模型参数之间的去相关过程,并且去相关矩阵是和并且其中矩阵E表示去相关矩阵,并且矩阵E-1表示矩阵的逆。
因此,去相关之后的等式(11)可以写成
我们使用Cauchy函数作为先验分布,使用高斯函数作为似然函数,然后可 以通过Cauchy先验分布和高斯似然函数(Pan,2018)联合解决后验分布,由 下式给出:
是地震噪声的方差,是模型参数的方差;p(·)代表概率密度函数。将低 频信息正则化与最大化方程(15)相结合,我们得到了
其中表示模型参数的正则化系数,在和中表示模型参数的初始值。
然后求解函数方程(16),我们得到了
其中
和
其中QCauchy表示Cauchy稀疏矩阵。
我们最终使用迭代重加权最小二乘(IRLS)算法来求解方程(17)并得到 模型参数的最终结果
根据最终结果可以合理地和地质地估计P波和S波模量和裂缝弱度。
其中,附录A:
对于PP波入射和反射的情况,极化和慢度向量由(Shaw和Sen,2006)给 出。
t=[sinθcosφ,sinθsinφ,cosθ], (A-1a)
t'=[-sinθcosφ,-sinθsinφ,cosθ], (A-1b)
其中φ代表方位角;p和t这两个参数是慢度矢量和极化矢量;ξ和η这两 个参数与慢度矢量p和极化矢量t相关联,表示为
和
相关参数ηmn如下
将方程(A-3a)和(A-3b)代入方程(8),我们推导出PP波反射系数的表 达式为
推导后,计算式(A-4)得到
其中
和g=μb/Mb.
可以推导出HTI介质中的PP波反射系数。当θ°=90°,我们可以得到
和
根据Pan和Zhang(2018)的研究,我们知道反射系数与HTI介质完全一致。 说明本实施例提供的方法的结果的有效性。
综上,在具有倾角的横向各向同性介质中传播的地震波的速度随着不同的 方向而变化,并且由裂缝弱度定义的各向异性可以表征方位角速度变化。为了 表征倾斜的横向各向同性,我们提出了使用可观察的PP波方位角反射数据的方 位AVO参数化和裂缝弱度的反演。根据法向和切向过剩弱度以及Bond变换对弱 各向异性VTI介质的刚度矩阵进行积分,我们首先根据法向和切向过剩弱度以 及倾角推导出TTI介质的刚度矩阵。将散射理论与TTI介质刚度矩阵的扰动相 结合,然后根据背景P-和S-波模量,密度,裂缝弱度和倾斜角得出线性化PP 波反射系数。分离两个弱各向异性TTI介质的弱对比界面。在贝叶斯框架内,我们接下来利用Cauchy稀疏正则化和低频信息正则化以及IRLS优化算法的约 束来实现方位AVO反演。我们最终使用合成和实际数据集来确认所提出的方位 AVO参数化和TTI介质中的裂缝弱度的反演的有效性。
本领域内的技术人员应明白,本公开的实施例可提供为方法、***、或计 算机程序产品。因此,本公开可采用完全硬件实施例、完全软件实施例、或结 合软件和硬件方面的实施例的形式。而且,本公开可采用在一个或多个其中包 含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、 CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本公开是参照根据本公开实施例的方法、设备(***)、和计算机程序产品 的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/ 或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或 方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式 处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机 或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流 程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备 以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的 指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流 程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使 得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理, 从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程 或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本公开的优选实施例而已,并不用于限制本公开,对于本领 域的技术人员来说,本公开可以有各种更改和变化。凡在本公开的精神和原则 之内,所作的任何修改、等同替换、改进等,均应包含在本公开的保护范围之 内。
上述虽然结合附图对本公开的具体实施方式进行了描述,但并非对本公开 保护范围的限制,所属领域技术人员应该明白,在本公开的技术方案的基础上, 本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本公开 的保护范围以内。
Claims (10)
1.一种倾斜横向各向同性介质中的裂缝弱度反演方法,其特征是:包括以下步骤:
使用垂直横向各向同性介质和Bond矩阵的刚度矩阵,根据裂缝弱度和倾角推导出倾斜横向各向同性介质的刚度矩阵,结合刚度扰动和散射理论,分离两个弱各向异性倾斜横向各向同性介质的弱对比界面,得到该界面下PP波反射系数的方位角,以倾角作为先验信息,利用部分入射角叠加方位地震数据实现方位角叠前反演,根据反演结果估计贝叶斯框架中的裂缝弱度。
2.如权利要求1所述的一种倾斜横向各向同性介质中的裂缝弱度反演方法,其特征是:根据法向和切向过剩弱度以及Bond变换对弱各向异性垂直横向各向同性介质的刚度矩阵进行积分,根据法向和切向过剩弱度以及倾角推导倾斜横向各向同性介质的刚度矩阵。
3.如权利要求2所述的一种倾斜横向各向同性介质中的裂缝弱度反演方法,其特征是:考虑由平面水平弱相对界面分离的两个均匀弱各向异性倾斜横向各向同性介质,并且通过使用P波和S波速度作为坐标轴来参数化半空间中的每个倾斜横向各向同性介质,融合密度、法向和切向裂缝弱度参数、以及对称轴和z轴之间的角度差和方位角,确定倾斜横向各向同性介质的刚度矩阵。
4.如权利要求1所述的一种倾斜横向各向同性介质中的裂缝弱度反演方法,其特征是:将散射理论与倾斜横向各向同性介质刚度矩阵的扰动相结合,根据背景P-和S-波模量、密度、裂缝弱度和倾斜角得出线性化PP波反射系数,分离两个弱各向异性倾斜横向各向同性介质的弱对比界面。
5.如权利要求1所述的一种倾斜横向各向同性介质中的裂缝弱度反演方法,其特征是:在贝叶斯框架下,进行方位角的反演,得到贝叶斯框架下的裂缝弱度。
6.如权利要求1所述的一种倾斜横向各向同性介质中的裂缝弱度反演方法,其特征是:利用Cauchy稀疏正则化和低频信息正则化以及迭代加权最小二乘法优化算法的约束来实现方位角的叠前反演,根据反演结果估计裂缝弱度。
7.如权利要求6所述的一种倾斜横向各向同性介质中的裂缝弱度反演方法,其特征是:使用Cauchy函数作为先验分布,使用高斯函数作为似然函数,通过Cauchy先验分布和高斯似然函数联合解决后验分布,得到分布最大化方程,将低频信息正则化与最大化方程相结合,使用迭代重加权最小二乘算法来求解,得到模型参数的最终结果,即反演结果。
8.一种倾斜横向各向同性介质中的裂缝弱度反演***,其特征是:被配置为运行在服务器或处理器上,执行权利要求1-7中任一项所述的倾斜横向各向同性介质中的裂缝弱度反演方法。
9.一种计算机可读存储介质,其特征是:其中存储有多条指令,所述指令适于由终端设备的处理器加载并执行权利要求1-7中任一项所述的倾斜横向各向同性介质中的裂缝弱度反演方法。
10.一种终端设备,其特征是:包括处理器和计算机可读存储介质,处理器用于实现各指令;计算机可读存储介质用于存储多条指令,所述指令适于由处理器加载并执行权利要求1-7中任一项所述的倾斜横向各向同性介质中的裂缝弱度反演方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910354239.1A CN110161562B (zh) | 2019-04-29 | 2019-04-29 | 倾斜横向各向同性介质中的裂缝弱度反演方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910354239.1A CN110161562B (zh) | 2019-04-29 | 2019-04-29 | 倾斜横向各向同性介质中的裂缝弱度反演方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110161562A true CN110161562A (zh) | 2019-08-23 |
CN110161562B CN110161562B (zh) | 2020-09-18 |
Family
ID=67633109
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910354239.1A Expired - Fee Related CN110161562B (zh) | 2019-04-29 | 2019-04-29 | 倾斜横向各向同性介质中的裂缝弱度反演方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110161562B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110687601A (zh) * | 2019-11-01 | 2020-01-14 | 中南大学 | 一种正交各向异性介质流体因子与裂缝参数反演方法 |
CN113009571A (zh) * | 2021-02-18 | 2021-06-22 | 中国矿业大学(北京) | 双相介质中的水平裂缝的反射系数和透射系数确定方法 |
CN113176612A (zh) * | 2021-03-26 | 2021-07-27 | 中国石油大学(华东) | 一种高角度裂缝密度反演方法 |
CN113189648A (zh) * | 2021-04-30 | 2021-07-30 | 西南石油大学 | 一种正交各向异性页岩脆性指数预测方法 |
CN113189647A (zh) * | 2021-04-30 | 2021-07-30 | 西南石油大学 | 一种横观各向同性页岩地层脆性指数预测方法 |
CN113568045A (zh) * | 2021-07-26 | 2021-10-29 | 中南大学 | 一种tti介质裂缝倾角的反演方法 |
CN114063163A (zh) * | 2021-11-11 | 2022-02-18 | 中南大学 | 一种裂缝型储层单斜等效介质地震表征与反演方法及*** |
CN114384067A (zh) * | 2021-12-27 | 2022-04-22 | 中国科学院福建物质结构研究所 | 一种各向同性激光介质中弱各向异性的测量装置及测量方法和应用 |
CN114488302A (zh) * | 2022-01-25 | 2022-05-13 | 中南大学 | 一种原位各向异性地应力场预测方法及*** |
CN115061194A (zh) * | 2022-06-02 | 2022-09-16 | 中国地质大学(北京) | 一种裂缝发育状态的识别方法、装置、电子设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050135189A1 (en) * | 2003-10-28 | 2005-06-23 | Baker Hughes Incorporated | Vector 3-component 3-dimensional kirchhoff prestack migration |
US20120185171A1 (en) * | 2008-12-03 | 2012-07-19 | Yonghe Sun | Multiple anisotropic parameter inversion for a tti earth model using well data |
CN104570072A (zh) * | 2013-10-16 | 2015-04-29 | 中国石油化工股份有限公司 | 一种粘弹性介质中的球面pp波反射系数建模方法 |
CN108627871A (zh) * | 2017-03-24 | 2018-10-09 | 中国石油化工股份有限公司 | 一种tti介质裂缝性质参数的反演方法 |
CN108897043A (zh) * | 2018-09-03 | 2018-11-27 | 中国石油天然气股份有限公司 | 基于弱度参数检测裂缝型储层流体的方法及装置 |
-
2019
- 2019-04-29 CN CN201910354239.1A patent/CN110161562B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050135189A1 (en) * | 2003-10-28 | 2005-06-23 | Baker Hughes Incorporated | Vector 3-component 3-dimensional kirchhoff prestack migration |
US20120185171A1 (en) * | 2008-12-03 | 2012-07-19 | Yonghe Sun | Multiple anisotropic parameter inversion for a tti earth model using well data |
CN104570072A (zh) * | 2013-10-16 | 2015-04-29 | 中国石油化工股份有限公司 | 一种粘弹性介质中的球面pp波反射系数建模方法 |
CN108627871A (zh) * | 2017-03-24 | 2018-10-09 | 中国石油化工股份有限公司 | 一种tti介质裂缝性质参数的反演方法 |
CN108897043A (zh) * | 2018-09-03 | 2018-11-27 | 中国石油天然气股份有限公司 | 基于弱度参数检测裂缝型储层流体的方法及装置 |
Non-Patent Citations (3)
Title |
---|
潘新朋 等: "裂缝-孔隙型含气储层流体与裂缝参数贝叶斯地震反演方法", 《中国科学:地球科学》 * |
潘新朋 等: "非均质HTI介质裂缝弱度参数地震散射反演", 《石油地球物理勘探》 * |
王涛: "各向同性/各向异性介质中弹性波场数值模拟方法技术研究", 《中国博士学位论文全文数据库》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110687601A (zh) * | 2019-11-01 | 2020-01-14 | 中南大学 | 一种正交各向异性介质流体因子与裂缝参数反演方法 |
CN113009571A (zh) * | 2021-02-18 | 2021-06-22 | 中国矿业大学(北京) | 双相介质中的水平裂缝的反射系数和透射系数确定方法 |
CN113176612A (zh) * | 2021-03-26 | 2021-07-27 | 中国石油大学(华东) | 一种高角度裂缝密度反演方法 |
CN113189648B (zh) * | 2021-04-30 | 2022-03-11 | 西南石油大学 | 一种正交各向异性页岩脆性指数预测方法 |
CN113189648A (zh) * | 2021-04-30 | 2021-07-30 | 西南石油大学 | 一种正交各向异性页岩脆性指数预测方法 |
CN113189647A (zh) * | 2021-04-30 | 2021-07-30 | 西南石油大学 | 一种横观各向同性页岩地层脆性指数预测方法 |
CN113189647B (zh) * | 2021-04-30 | 2022-03-11 | 西南石油大学 | 一种横观各向同性页岩地层脆性指数预测方法 |
CN113568045A (zh) * | 2021-07-26 | 2021-10-29 | 中南大学 | 一种tti介质裂缝倾角的反演方法 |
CN114063163A (zh) * | 2021-11-11 | 2022-02-18 | 中南大学 | 一种裂缝型储层单斜等效介质地震表征与反演方法及*** |
CN114063163B (zh) * | 2021-11-11 | 2022-11-11 | 中南大学 | 一种裂缝型储层单斜等效介质地震表征与反演方法及*** |
CN114384067A (zh) * | 2021-12-27 | 2022-04-22 | 中国科学院福建物质结构研究所 | 一种各向同性激光介质中弱各向异性的测量装置及测量方法和应用 |
CN114384067B (zh) * | 2021-12-27 | 2023-02-10 | 中国科学院福建物质结构研究所 | 一种各向同性激光介质中弱各向异性的测量装置及测量方法和应用 |
CN114488302A (zh) * | 2022-01-25 | 2022-05-13 | 中南大学 | 一种原位各向异性地应力场预测方法及*** |
CN115061194A (zh) * | 2022-06-02 | 2022-09-16 | 中国地质大学(北京) | 一种裂缝发育状态的识别方法、装置、电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN110161562B (zh) | 2020-09-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110161562A (zh) | 倾斜横向各向同性介质中的裂缝弱度反演方法及*** | |
US20180259664A1 (en) | Automated system and methods for adaptive robust denoising of large-scale seismic data sets | |
Kim et al. | Discrimination of earthquakes and explosions in southern Russia using regional high-frequency three-component data from the IRIS/JSP Caucasus network | |
US9121968B2 (en) | Extracting geologic information from multiple offset stacks and/or angle stacks | |
CN105137486B (zh) | 各向异性介质中弹性波逆时偏移成像方法及其装置 | |
US8972195B2 (en) | Extracting geologic information from multiple offset stacks and/or angle stacks | |
US8838391B2 (en) | Extracting geologic information from multiple offset stacks and/or angle stacks | |
Sambridge et al. | Seismic event location: nonlinear inversion using a neighbourhood algorithm | |
CN107153216B (zh) | 确定地震波场的坡印廷矢量的方法、装置以及计算机存储介质 | |
MX2007012979A (es) | Restoracion de falla rapida 3d. | |
Trifu et al. | A fast evaluation of the seismic moment tensor for induced seismicity | |
CN109164483A (zh) | 多分量地震数据矢量去噪方法及多分量地震数据矢量去噪装置 | |
Diallo et al. | Characterization of polarization attributes of seismic waves using continuous wavelet transforms | |
Jin et al. | Reflection and transmission approximations for weak contrast orthorhombic media | |
Zhang et al. | Prompt elasto-gravity signals (PEGS) and their potential use in modern seismology | |
Grechka | Shear-wave group-velocity surfaces in low-symmetry anisotropic media | |
EP2671100A2 (en) | Extracting geologic information from multiple offset stacks and/or angle stacks | |
WO2014102616A2 (en) | Efficient wavefield extrapolation in anisotropic media | |
Jin et al. | Reflection and transmission approximations for monoclinic media with a horizontal symmetry plane | |
Tavakoli F et al. | Matrix-free anisotropic slope tomography: Theory and application | |
CN109541682A (zh) | 各向同性弹性参数保幅反演方法及装置 | |
Fu et al. | Time-frequency analysis of seismic wave propagation across a rock mass using the discontinuous deformation analysis method | |
Stovas | On parameterization in monoclinic media with a horizontal symmetry plane | |
You et al. | Multiscale ambient noise tomography of short-period Rayleigh waves across northern Taiwan | |
Silver et al. | A source extent analysis of the Imperial Valley earthquake of October 15, 1979, and the Victoria earthquake of June 9, 1980 |
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 | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20191106 Address after: 410083 Hunan province Changsha Lushan Road No. 932 Applicant after: Central South University Applicant after: China University of Petroleum (East China) Address before: 266580 Qingdao economic and Technological Development Zone, Changjiang Road, No. 66, Shandong Applicant before: China University of Petroleum (East China) |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200918 |
|
CF01 | Termination of patent right due to non-payment of annual fee |