CN106443775B - 高分辨率转换波裂缝预测方法 - Google Patents

高分辨率转换波裂缝预测方法 Download PDF

Info

Publication number
CN106443775B
CN106443775B CN201610352020.4A CN201610352020A CN106443775B CN 106443775 B CN106443775 B CN 106443775B CN 201610352020 A CN201610352020 A CN 201610352020A CN 106443775 B CN106443775 B CN 106443775B
Authority
CN
China
Prior art keywords
formula
indicate
wave
reflectance factor
resolution conversion
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.)
Active
Application number
CN201610352020.4A
Other languages
English (en)
Other versions
CN106443775A (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.)
China National Petroleum Corp
BGP Inc
Original Assignee
BGP Inc
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 BGP Inc filed Critical BGP Inc
Priority to CN201610352020.4A priority Critical patent/CN106443775B/zh
Publication of CN106443775A publication Critical patent/CN106443775A/zh
Application granted granted Critical
Publication of CN106443775B publication Critical patent/CN106443775B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/301Analysis for determining seismic cross-sections or geostructures

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范数联合优化算法求取高精度反射系数;构建高频子波同时保持该高频子波具有确定的频率宽度,用高频子波与高精度反射系数得到高分辨率转换波偏移剖面;利用改进的第三代特征值相干体计算法,对高分辨率转换波偏移剖面进行相干数据体切片提取,从而实现裂缝预测。本发明能够使得提取的裂缝预测剖面分辨率更高、裂缝走向更加清晰、细节刻画更加明显。

Description

高分辨率转换波裂缝预测方法
技术领域
本发明属于转换波地震资料处理技术领域,具体来讲,涉及一种高分辨率转换波裂缝预测方法。
背景技术
目前,相干体裂缝特征描述是探寻裂缝性油气藏的一项重要手段,而获取裂缝预测所需高分辨率数据体、裂缝参数的求取一直是裂缝油气藏勘探中重点研究的问题,尤其转换波低信噪比和分辨率的特点,需要提高其分辨率,便于后续解释的薄层分辨。但是,最小相位和白噪假设为前提的常规高分辨率处理方法,在主频带范围限制下,提高转换波分辨率的能力有限,无法达到后续裂缝预测分辨多薄层最小厚度的要求,所以需要采取新的高分辨率处理思路,即:转换波高精度反射系数求取与高频子波高分辨率重构的新思路,该思路不受任何假设和频带限制,能完成常规手段无法分辨多薄层最小厚度的要求,获取裂缝预测满意的高品质剖面。
通常,反射系数的求取需要构建多薄层目标函数,然而,由于常规方法构建的目标函数存在诸多不合理之处,故而,在求取过程中的反射系数存在极值和位置、大小上的误差,这种误差导致后续高分辨率剖面的错像。此外,目标函数的求解方法有很多,比如模拟退火法、共轭梯度法、匹配追踪法等这些算法,虽然都能求取目标函数的反射系数,但精度和稳定性较差,且对新构建的目标函数适应性较差。
三维相干体技术是20世纪90年代后期兴起的一项十分有效的地震解释技术,该技术主要从相邻地震道相互之间的相干性出发,给出一种定量描述,它能突出断层处地震波的变化情况,形成相干值异常区域,从而预测断裂展布。M.Bahorich等提出在振幅上进行互相关相干分析算法,但该算法计算出的断点非相干性点模糊,噪声重、能量弱。K.J.Marfurt年提出基于相似系数相干算法具有稳定性好、断层刻画精度高的特点,但横向分辨率低的劣势无法改变。Gersztenkorn等提出基于本征结构分析的算法,利用Hilbert变换计算相干体提高其稳定性和抗噪性,但该算法用于高陡复杂地区时,如不考虑倾角则不如第二代算法效果好,如考虑倾角则计算效率无法应用于生产中。
发明内容
本发明的目的在于解决现有技术存在的上述不足中的至少一项。例如,本发明的目的之一在于提供一种具有更高分辨率的转换波裂缝预测方法,以使得提取的裂缝预测剖面分辨率更高、裂缝走向更加清晰、细节刻画更加明显。
为了实现上述目的,本发明提供一种高分辨率转换波裂缝预测方法。所述方法包括:采集转换波地震数据,进行预处理和叠前偏移后得到转换波偏移剖面,在转换波偏移剖面上取n道地震记录x(x,t),x=1,2,…,n;在偏移剖面上建立稀疏约束和波阻抗约束的多薄层目标函数;对多薄层目标函数利用频域映射与L1范数联合优化算法求取高精度反射系数;构建高频子波同时保持该高频子波具有确定的频率宽度,用高频子波与高精度反射系数得到高分辨率转换波偏移剖面;利用改进的第三代特征值相干体计算法,对高分辨率转换波偏移剖面进行相干数据体切片提取,从而实现裂缝预测;其中,
所述目标函数为下式(1):
其中,r表示反射系数,ae代表反射系数的偶分量比例,ao代表反射系数的奇分量比例,re为反射系数的偶分量,ro为反射系数的奇分量,A为变换矩阵,b为地震数据和子波相关矩阵,λ为稀疏因子,为阻抗因子,C为积分算子矩阵,ξp为纵波波阻抗矩阵。
与现有技术相比,本发明的有益效果包括:能够基于转换波地震数据构建更加合理的多薄层目标函数,并能有效提高目标函数求解精度和抗噪性,能够获得高精度反射系数,能够使得提取的裂缝预测剖面分辨率更高、裂缝走向更加清晰、细节刻画更加明显。
附图说明
图1示出了通过现有的预处理和叠前偏移后得到的转换波偏移剖面。
图2示出了根据本发明的方法求取高精度反射系数对应的剖面。
图3示出了根据本发明的方法得到的高分辨率转换波偏移剖面。
图4A示出了图1的转换波偏移剖面的振幅谱数据;图4B示出了根据本发明的方法得到的高分辨率转换波偏移剖面的振幅谱数据。
图5A示出了利用改进的第三代特征值相干体计算法对图1的转换波偏移剖面进行相干切片提取的效果图;图5B示出了根据本发明的方法对高分辨率转换波偏移剖面进行相干切片提取的效果图。
具体实施方式
在下文中,将结合附图和示例性实施例来详细说明本发明的高分辨率转换波裂缝预测方法。
在本发明的一个示例性实施例中,高分辨率转换波裂缝预测方法可通过下面的步骤(1)至步骤(5)来实现:
(1)对采集到的转换波地震数据,经过常规预处理和叠前偏移后得到的转换波偏移剖面,取n道地震记录x(x,t),x=1,2,…,n;其中,x表示道数,t表示时间。
(2)在偏移剖面上建立由稀疏约束和波阻抗约束双重约束的多薄层目标函数;
步骤(2)可按下列方式完成:
(a)对于不含噪声的褶积模型,地震合成记录在频率域中可表示为:
S(f)=W(f)·R(f) (1)
其中,S(f)表示频率域下地震记录(即,上面的x(x,t)),W(f)表示频率域下地震子波,R(f)表示频率域下反射系数;f为频率,时间t经过傅里叶变换后变换到频率f中。
利用频率域下地震记录与地震子波的丰富信息,利用合适的常规反演算法(如,共轭梯度法、匹配追踪法、模拟退火等)消除子波(即,地震子波)影响(如,子波主频、旁瓣干扰等会影响反演求取反射系数精度的内容),从记录中得到所需的反射系数(这里,为了得到反射系数R(f))。利用反射系数奇偶分解原理,可得到频率域下目标函数表达式:
其中,ae代表反射系数R(f)的偶分量比例,ao代表反射系数R(f)的奇分量比例;Re表示R(f)的实部,Im表示R(f)的虚部。
(b)为建立更符合实际地层情况的反射系数模型,对步骤(a)中的多薄层反射系数模型下的目标函数进行推导或改进。
由奇偶反射系数分解原理,得多薄层反射系数奇偶分量可表示为:
re(t)代表反射系数r(t)的偶分量,ro(t)代表反射系数r(t)的奇分量,t表示时间,T表示时间厚度;r(t)为时间域反射系数,r(t)经过傅里叶变换后到上面的频率域反射系数R(f);g(t)当前分析位置反射系数,ge(t)代表g(t)的偶分量,go(t)代表g(t)的奇分量。
对式(3)进行傅里叶变换,得到其频率域实部与虚部表达式:
其中,π为圆周率。
利用式(2)和式(4)得到多薄层反射模型的目标函数表达式为:
其中,re代表反射系数r(t)的偶分量,ro代表反射系数r(t)的奇分量,Δt表示时移量;i无物理意义,是指数与正弦函数变化标识符。
为了便于求解,将式(5)可写成矩阵形式:
其中,
其中,t为时间,t1、t2、……、tN分别表示第1、2、……、N个时间,f为频率,f1、f2、……、fM分别表示第1、2、……、M个分析频率范围内的频率,T1、T2……、TN/2分别表示第1、2、……、N/2个时间厚度,Δt表示时移量
(c)由于该目标函数(6)(即,式(6))有多个全局最优解,为了降低求解问题时的多解性和不确定性,需要加上稀疏化约束,将目标函数改为范数求解形式为:
其中,λ为稀疏因子,λ一般在0-1范围内取值,|| ||1表示L1范数,表示欧几里得(Euclidean)范数的平方。
(d)由于波阻抗可看作是反射系数对时间的积分,因此发明人把它引入到目标函数中作为反射系数r的先验约束。
由于转换波相对于纵波来说其信噪比和分辨率较低,为了提高反射系数精度,本发明使用分辨率更高的纵波波阻抗作为转换波的约束条件。
若I(t)为纵波波阻抗值,则纵波相对波阻抗εt为:
其中,I(t0)表示纵波初始波阻抗值,r为反射系数,η表示在时间t0与时间t中取值符合,t0表示初始时间,tM-1表示结束时间。
把上式(8)写成矩阵形式:
ξp=Cr (9)
其中C为积分算子矩阵,其离散形式可表示为:
则由最小平方定义波阻抗约束得到新的目标函数为:
其中,r表示反射系数,ae和ao分别代表反射系数r(t)的偶奇分量比例,re是反射系数的偶分量,ro是反射系数的奇分量;A为变换矩阵,b是地震数据和子波相关矩阵,λ是稀疏因子,是阻抗因子(其大于零即可),C为积分算子矩阵,ξp是纵波波阻抗矩阵)
(3)对于步骤(2)中的式(11),利用频域映射与L1范数联合优化算法求取高精度反射系数。
本步骤(3)可按下列方式完成:
(a)对式(11)进行求解中,由于波阻抗约束项是由测井信号得到的已知先验约束条件,所以式(11)的求解类似于求矩阵方程Ax=b在含有噪声情况下的稀疏解。
可将目标函数(11)的求解问题转化为下式:
其中,υ为极小权重值,G为复数集合;本式(12)中,A表示Ae时b表示be,X表示re;A表示Ao时b表示bo,X表示ro
如果引入变量q∈Gm,那么式(12)等价于式(13):
式(13)对应的增广拉格朗日子问题可表示为:
其中,yT表示拉格朗日乘子的共轭转置运算,β为罚参数。
利用式(14)进行频域映射与L1范数联合优化算法实现步骤如下:
1)令k=0对qk,xk,yk赋初始值,并给定υ,β常数值(这里,υ,β为人为给定,可为大于零的常值),然后进行以下步骤2)至5)的算法运算,如果满足终止准则(例如,终止准则可以为:最优解值或者迭代次数),则算法完成运算;否则,进行步骤2);
2)令x=xk,y=yk,求解r的子问题得到qk+1
3)令q=qk+1,y=yk则关于x的极小化问题式(14)等价于式(16):
那么可以通过式(17)近似求解来完成对式(16)的精确求解:
其中,τ为大于零的邻近参数(这里,τ为人工给定,其大于零即可),gk为x=xk时二次项的梯度,那么式(17)通过下式(18)求解得到xk+1
其中,ο表示逐元素相乘,Shrink(,)表示一维收缩算子,sign表示符号函数。
4)令x=xk+1,q=qk+1,求解y的子问题得到yk+1
yk+1=yk-γβ(Axk+1`+qk+1-b) (19)
其中,γ为大于零的常数(这里,γ人工给定,其大于零即可)。
5)令k=k+1,重复步骤1)~步骤4)进行迭代运算。
(4)基于步骤(3)反射系数与高频子波重构,得到高分辨率转换波偏移剖面。
构建高频子波的同时需要保持该高频子波具有一定的频带宽度,以使得高分辨率转换波偏移剖面是宽频的。可利用常规的精度较高的俞式子波建立高频子波信号,亦可用诸如雷克子波等零相位子波建立高频子波信号。
时间域中俞氏子波表达式为:
其中q、p分别表示子波的低频值和高频值,t为时间。
用式(20)的子波与步骤(3)求取的高精度反射系数得到高分辨率转换波偏移剖面。
(5)对步骤(4)的高分辨率转换波偏移剖面进行改进的第三代特征值相干算法C3相干体切片裂缝预测。
相干数据体技术主要利用相邻地震道的相似性原理,描述地层及岩性的横向不均匀性。通过相干数据体算法,对地震数据体的不连续性进行系统分析,使解释人员能有效识别构造和断层的分布情况,避免解释的随意性。
(a)为了尽可能压制噪声干扰,提高相干体计算数据的分辨率,Gersztenkorn等提出了第三代特征值相干体计算法。
首先,利用Marfurt等在第二代多道相似性相干体计算中定义的椭圆和矩形分析窗口中,以坐标为(x,y)的分析点为中心,时窗中心t=nΔt,则时间窗口内含2k+1个元素的分析道与窗口内的J道数据排列成矩阵:
其中,表示第j道在时间t=mΔt-pΔxj-qyj沿着视倾角(p,q)插值后的振幅值。
那么矩阵U的第m列是分析体内每个地震道在时间t=mΔt-pΔxj-qyj(j=1,2,…,J)的振幅值,则沿着一对视倾角(p,q),时窗中心在t=mΔt的分析窗口内的2k+1个采样点的数据协方差矩阵为:
如λj为分析窗口内地震数据协方差矩阵式(22)的特征值,且有λ1≥λ2≥…≥λJ
则在以分析点为中心的分析窗口内,沿着视倾角(p,q)的J道c3特征值相干计算算法式可表示为:
则在以分析点为中心的分析窗口内,沿视倾角(p,q=0)的J道c3特征值相干计算算法式可表示为:
则沿着所有倾角的c3特征值相干计算算法式可表示为:
(b)因为第三代特征值相干算法C3是在子空间中计算,比前两代相干算法更能很好地压制噪声,但其计算量大、效率慢的缺点却影响了其在实际生产中的经济效益。
对该算法进行常规改进。例如,可利用基于相似性算法的高效率性估算各点的视倾角和方位角,然后在大于10倍的分析窗口内计算视倾角和方位角的平均值再在以分析点为中心的分析窗口内,沿着区域平滑后的视倾角计算地震数据协方差矩阵的特征值及相干值:
利用式(26)对步骤(4)的高分辨率剖面进行相干数据体切片提取,从而进行裂缝预测。
图1示出了通过现有的预处理和叠前偏移后得到的转换波偏移剖面,其横坐标表示CMP道号,图1中也示出了测井系数(即,图1中位于道号791和841之间的纵向曲线)。图2示出了根据本发明的方法求取高精度反射系数对应的剖面。由图2可以看出,本发明的方法求取高精度反射系数与测井系数(纵向曲线)匹配较好,各层位置和极性求取较佳。
图3示出了根据本发明的方法得到的高分辨率转换波偏移剖面。可以看出,图3明显比图1的分辨率高,且浅、中、深层的同相轴(例如,800ms附近的浅层、1000ms附近的中层、1150ms附近的深层的同相轴)连续性更好、断点更清晰、断层更明了,便于后续解释人员的裂缝预测处理。
图4A示出了图1的转换波偏移剖面的振幅谱数据;图4B示出了根据本发明的方法得到的高分辨率转换波偏移剖面的振幅谱数据。由图4A和图4B可知,经本发明的方法处理而得到的高分辨率转换波偏移剖面的振幅谱数据的频带得到明显拓宽,而且其振幅谱低频信息得到完整的保留。
图5A示出了利用改进的第三代特征值相干体计算法对图1的转换波偏移剖面进行相干切片提取的效果图;图5B示出了根据本发明的方法对高分辨率转换波偏移剖面进行相干切片提取的效果图。由图5A和图5B可知,经本发明的方法处理后的切片相干效果更好,裂缝刻画细节更清晰,避免了低分辨率数据为后续解释人员带来的假象。
综上所述,本发明通过利用稀疏约束和波阻抗双重约束构建多薄层目标函数,提高目标函数求解精度和抗噪性;利用频域映射与L1范数联合优化算法求取高精度反射系数,并通过与高频子波重构得到分辨率更高的新地震剖面,改善转换波的分辨率效果,利用改进的C3相干体切片进行裂缝预测,可有效提高地震勘探裂预测和描述地下构造的能力。也就是说,本发明的高分辨率转换波裂缝预测方法能够使得提取的裂缝预测剖面分辨率更高、裂缝走向更加清晰、细节刻画更加明显。
尽管上面已经结合附图和示例性实施例描述了本发明,但是本领域普通技术人员应该清楚,在不脱离权利要求的精神和范围的情况下,可以对上述实施例进行各种修改。

Claims (3)

1.一种高分辨率转换波裂缝预测方法,所述方法包括:采集转换波地震数据,进行预处理和叠前偏移后得到转换波偏移剖面,在转换波偏移剖面上取n道地震记录x(x,t),x=1,2,…,n,t表示时间;
在偏移剖面上建立稀疏约束和波阻抗约束的多薄层目标函数;
对多薄层目标函数利用频域映射与L1范数联合优化算法求取高精度反射系数;
构建高频子波同时保持该高频子波具有确定的频率宽度,用高频子波与高精度反射系数得到高分辨率转换波偏移剖面;
利用改进的第三代特征值相干体计算法,对高分辨率转换波偏移剖面进行相干数据体切片提取,从而实现裂缝预测;
其中,
所述目标函数为下式(1):
其中,r表示反射系数,ae代表反射系数的偶分量比例,ao代表反射系数的奇分量比例,re为反射系数的偶分量,ro为反射系数的奇分量,A为变换矩阵,b为地震数据和子波相关矩阵,λ为稀疏因子,为阻抗因子,C为积分算子矩阵,ξp为纵波波阻抗矩阵,其中,
其中,t为时间,t1、t2、……、tN分别表示第1、2、……、N个时间,f为频率,f1、f2、……、fM分别表示第1、2、……、M个分析频率范围内的频率,T1、T2……、TN/2分别表示1、2、……、N/2个时间厚度,Δt表示时移量,Re表示频率域下反射系数R(f)的实部,Im表示频率域下反射系数R(f)的虚部,S表示频率域下地震记录x(x,t),W表示频率域下地震子波;
其中,t0表示初始时间,tM-1表示结束时间,纵波相对波阻抗εt为:
I(t)为纵波波阻抗值,I(t0)表示纵波初始波阻抗值,r为反射系数,η表示在时间t0与时间t中取值符合;
C的离散形式表示为
2.如权利要求1所述的高分辨率转换波裂缝预测方法,其中,所述求取高精度反射系数的步骤包括:
A、将多薄层目标函数式(1)的求解问题转化为式(2):
其中,υ为极小权重值,G为复数集合;
B、引入变量q∈Gm,则式(2)等价于式(3):
C、式(3)对应的增广拉格朗日子问题表示为式(4):
其中,yT表示乘子的共轭转置运算,β为罚参数;
D、利用式(4)进行频域映射与L1范数联合优化算法求取高精度反射系数。
3.如权利要求2所述的高分辨率转换波裂缝预测方法,其中,所述步骤D包括步骤:
1)令k=0对rk,xk,yk赋初始值,并给定υ,β常数值,然后进行下面步骤2)至5)的算法运算,若满足终止准则,则完成运算,否则,进行步骤2);
2)令x=xk,y=yk,求解r的子问题得到是(5):
3)令r=rk+1,y=yk则关于x的极小化问题式(4)等价于式(6):
那么可以通过式(7)近似求解来完成对式(6)的精确求解,式(7)为:
其中,τ为大于零的邻近参数,gk为x=xk时二次项的梯度,那么式(7)通过式(8)求解得到xk+1,式(8)为:
其中,表示逐元素相乘,Shrink(,)表示一维收缩算子;
4)令x=xk+1,r=rk+1,求解y的子问题得到式(9):
yk+1=yk-γβ(Axk+1`+rk+1-b),
其中,γ为大于零的常数。
5)令k=k+1,重复步骤1)~步骤4)进行迭代运算。
CN201610352020.4A 2016-05-25 2016-05-25 高分辨率转换波裂缝预测方法 Active CN106443775B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610352020.4A CN106443775B (zh) 2016-05-25 2016-05-25 高分辨率转换波裂缝预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610352020.4A CN106443775B (zh) 2016-05-25 2016-05-25 高分辨率转换波裂缝预测方法

Publications (2)

Publication Number Publication Date
CN106443775A CN106443775A (zh) 2017-02-22
CN106443775B true CN106443775B (zh) 2018-09-28

Family

ID=58183493

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610352020.4A Active CN106443775B (zh) 2016-05-25 2016-05-25 高分辨率转换波裂缝预测方法

Country Status (1)

Country Link
CN (1) CN106443775B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107894616B (zh) * 2017-11-14 2020-01-10 中国石油集团川庆钻探工程有限公司地球物理勘探公司 多分量转换波裂缝预测方法
CN107942373B (zh) * 2017-11-14 2020-02-07 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于裂缝性油气储层断裂系统检测的相干算法
CN108535775B (zh) * 2018-03-30 2019-08-23 中国石油大学(北京) 非平稳地震资料声波阻抗反演方法及装置
CN111435170B (zh) * 2019-01-11 2022-11-04 中国石油天然气集团有限公司 薄层偏移剖面的计算方法及系统
CN111751877A (zh) * 2019-03-26 2020-10-09 中国石油天然气股份有限公司 一种地震数据多重积分相干裂缝预测方法与装置
CN109900805A (zh) * 2019-04-08 2019-06-18 大连理工大学 基于频域稀疏反演的tofd盲区内缺陷定量检测方法
CN112462442B (zh) * 2020-11-30 2022-04-08 山东大学 重磁位场场源位置估计方法、系统、介质及电子设备
CN112711861B (zh) * 2021-01-15 2022-05-17 青岛海洋地质研究所 一种确定高分辨率拖拽式浅剖多次波起跳时间的方法
CN114089428B (zh) * 2021-10-29 2022-08-12 中国矿业大学 一种时间域电磁张量测井方法
CN114002738B (zh) * 2021-11-02 2023-08-15 中海石油(中国)有限公司 高角度断裂约束的叠前裂缝预测方法、装置及存储介质
CN117055116B (zh) * 2023-10-12 2023-12-19 东北石油大学三亚海洋油气研究院 一种强稀疏约束Radon变换多次波压制方法及系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6820010B1 (en) * 2003-04-30 2004-11-16 Conocophillips Company Method for determining shear-wave velocity model for depth migration of mode-converted data
CN103076623B (zh) * 2011-10-25 2015-07-22 中国石油化工股份有限公司 一种基于叠前相干的裂缝检测方法
CN105445787B (zh) * 2014-08-11 2017-10-20 中国石油化工股份有限公司 一种最优方位子体相干的裂缝预测方法
CN104316958B (zh) * 2014-10-20 2017-12-05 中国石油天然气集团公司 一种识别不同尺度地层断裂的相干处理方法
CN104391325B (zh) * 2014-12-09 2017-02-22 中国石油天然气股份有限公司 不连续非均质地质体检测方法和装置

Also Published As

Publication number Publication date
CN106443775A (zh) 2017-02-22

Similar Documents

Publication Publication Date Title
CN106443775B (zh) 高分辨率转换波裂缝预测方法
Araya-Polo et al. Deep-learning tomography
CN108037526B (zh) 基于全波波场vsp/rvsp地震资料的逆时偏移方法
CN106405651B (zh) 一种基于测井匹配的全波形反演初始速度模型构建方法
CN103492910B (zh) 时域中的同步小波提取和反卷积
CN102298155B (zh) 一种基于高维小波变换的地震资料不连续性检测方法
CN108020863A (zh) 一种基于地震奇偶函数的碳酸盐岩薄储层孔隙度预测方法
Gong et al. Separation of prestack seismic diffractions using an improved sparse apex-shifted hyperbolic Radon transform
CN108845317B (zh) 一种基于分层介质格林函数的频域逆时偏移方法
CN107831542A (zh) Ddw高精度深度域井震匹配方法
CN103399346A (zh) 一种井震联合初始波阻抗建模方法
Li et al. Seismic coherence for discontinuity interpretation
CN107390213A (zh) 一种基于滑动时窗的探地雷达记录剖面的时延曲线提取方法
CN104730576A (zh) 基于Curvelet变换的地震信号去噪方法
CN109633593A (zh) 一种地质雷达信号定量分析方法及系统
CN106842299B (zh) 一种基于地震属性的裂缝定量化预测的方法
CN110261905A (zh) 基于倾角控制的复值相干微断层识别方法
CN107479091B (zh) 一种提取逆时偏移角道集的方法
Xu et al. Ground-roll separation of seismic data based on morphological component analysis in two-dimensional domain
CN107356963B (zh) 一种数据驱动的自适应的地震信号相干体属性分析方法
CN104422956A (zh) 一种基于稀疏脉冲反演的高精度地震谱分解方法
CN105005073A (zh) 基于局部相似度和评价反馈的时变子波提取方法
CN109738852A (zh) 基于低秩矩阵重建的分布式源二维空间谱估计方法
CN108549102A (zh) 联合梯度结构张量和多窗分析的地层结构曲率估计方法
Yablokov et al. Uncertainty quantification of multimodal surface wave inversion using artificial neural networks

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan

Applicant after: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd.

Address before: No. 216, No. 216, Huayang Avenue, Huayang Town, Shuangliu County, Shuangliu County, Sichuan

Applicant before: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd.

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20180330

Address after: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei

Applicant after: Dongfang Geophysical Exploration Co., Ltd., China Petrochemical Corp.

Address before: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan

Applicant before: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd.

GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20201117

Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No.

Patentee after: CHINA NATIONAL PETROLEUM Corp.

Patentee after: BGP Inc., China National Petroleum Corp.

Address before: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei

Patentee before: BGP Inc., China National Petroleum Corp.