CN104730579A - 一种基于表层横波速度反演的纵横波联合静校正方法 - Google Patents
一种基于表层横波速度反演的纵横波联合静校正方法 Download PDFInfo
- Publication number
- CN104730579A CN104730579A CN201310520027.9A CN201310520027A CN104730579A CN 104730579 A CN104730579 A CN 104730579A CN 201310520027 A CN201310520027 A CN 201310520027A CN 104730579 A CN104730579 A CN 104730579A
- Authority
- CN
- China
- Prior art keywords
- wave
- partiald
- velocity
- dispersion curve
- inversion
- 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
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于表层横波速度反演的纵横波联合静校正方法,该方法通过高精度的波场插值,联合应用纵波和转换波资料面波信息,并通过加权叠加,有效提高了直接反演横波速度的精度,有效的解决了常规方法通过间接求取横波静校正量,存在静校正不准的问题,从而有效提高多波资料的成像水平。经过对胜利油田罗家地区三分量地震数据处理中的应用试验,取得了较好的处理效果,证明基于表层横波速度反演的纵横波联合静校正方法是一可行的方法技术。
Description
技术领域
本发明涉及地震勘探或勘测领域,具体涉及一种基于表层横波速度反演的纵横波联合静校正方法。
背景技术
由于地形的起伏变化或浅层存在横向的速度不均匀体,造成到达地表的反射地震信号在不同的检波器上除存在正常的偏移距和深度变化引起的正常时差外,还会出现由于地表与浅层低、降速带异常所引起的走时变化。这一走时变化必须在反射波信号的处理中予以消除,才能获得地下介质的正常的成像。这一时差的校正在人工地震信号处理领域被称为静校正。
静校正一般分为二大类方法:一类方法是通过获得地表的速度结构,将横向速度不均匀体产生的异常时差予以校正,这种方法一般又统称为一次静校正或长波长静校正;另一类方法是利用相邻地震道上信号的相关性,通过相关分析求取各道间的时差以达到消除这种时差的方法。这种方法在地震信号处理领域一般又称为剩余静校正,即这一静校正方法的应用一般在一次静校正之后使用,目的是消除各道间变化不大的短波长校正量。
纵波的静校正方法目前相对来说比较成熟,有层析静校正、折射波静校正等,但是由于转换波特有的反射特征使得这些纵波静校正方法无法应用于转换波,而由于转换波静校正需要炮点和检波点分开处理,使得其静校正方法更加复杂和困难,尤其是检波点的静校正问题更为突出。目前也出现了一些像共检波点叠加相关法静校正,基于构造时间控制的转换波静校正等常规转换波静校正方法,这些方法能解决部分静校正问题,但是由于间接应用纵波速度或者速度比,以及由于存在的纵横波匹配的等难题,使得这些方法无法准确求取横波静校正量,影响静校正效果。
公开号为CN101907727A的专利公开了“一种面波多分量转换波静校正方法”,该方法基于炮记录中的瑞雷面波,通过面波的提取、频散曲线反演求取浅层的地层分层及横波速度,进而精确计算转换波检波点静校正量。该方法和本发明相比,本发明建立了基于小波变换方法的高精度的道间插值,进一步提高了资料的信噪比,综合应用纵横波资料的面波信息,并进行能谱加权叠加,充分利用了纵波和横波的高频和低频优势,提高能谱拾取精度,应用了阻尼最小二乘法进行反演,进一步增加了计算的稳定性和效率。
芦俊等的“面波法转换波静校正研究”(《中国地球物理2006》.2006.25),简单的介绍了利用面波进行浅层横波速度反演的一些思路,而且提出只是利用垂直分量的面波信息,整个技术方法还没有实现。本发明提出了综合利用纵横波资料的优势面波信息,进一步提高了反演的精度,降低了多解性。而且在反演方法等方面采用了更加先进的方法,并进行了可视化实现。
另外,像姚姚等提出的利用转换折射波和反射波进行转换波静校正的方法也取得了一定的效果,但因为转换折射波是续至波不是初至波不易拾取以及拾取反射波时存在不确定性,给计算静校正量造成了一定的困难。
发明内容
静校正的效果会严重影响成像的质量,但是一般情况下,由于很难得到近地表的横波速度信息,使得转换波静校正处理中检波点静校正量难以准确计算,影响了多波资料的成像质量。本法的目的在于提供一种基于表层横波速度反演的纵横波联合静校正方法,通过多波资料中信息丰富的面部信息来直接反演横波速度信息,进而直接得到横波静校正量,减少中间过程的累积误差,提高静校正的精度。
一种基于表层横波速度反演的纵横波联合静校正方法,包括如下步骤:
第一步:波场高精度插值
所述波场高精度插值的具体方法是:
(a)对面波数据进行一定程度的线性动校(LNMO);
(b)进行小波变换插值;
(c)反线性动校(RLNMO)恢复信号。
第二步:面波频散曲线提取
提取频散曲线可以通过两个线性变换实现,具体方法是:
(a)对原始数据做倾斜叠加,将数据变换到了慢度-截距域(p-τ);所述对原始数据做倾斜叠加采用线性Radon变换来完成;
(b)沿截距τ方向做一维傅里叶变换得到慢度-频率域(p-f);
(c)将慢度转换成速度即可获得速度-频率域。
第三步:横波速度反演
横波速度的反演需要经过面波正演和面波反演两个方面;
1)面波正演
面波正演主要是根据弹性波动理论建立描述面波相速度C与波动频率f以及其它弹性参数之间关系的面波频散方程F(C,f),并解方程F(C,f)=0的根,即计算不同频率时的面波传播速度,由此了解特定介质模型中的面波频散特征。
2)面波反演
建立的面波频散方程:
其中:其中:[J]是由传播矩阵表示的连乘形式,△n为纵波体积膨胀系数,Wn为横波角位移,C为相速度,μ0为介质的拉梅常数,w0为法向位移。
因为[J]中的元素是相速度C和常数R的函数,所以式(1)表示了C与R的关系,从而导出了相速度C的频散函数
F(C,f)=0;
通过二分法求解得到不同频率下的速度值,最终得到面波速度随频率的变化——频散曲线;
面波反演的主要包括频散曲线C-f的求取和对频散曲线的解释,
反演过程中的目标函数为
Φ(X)=∑[di-Ci(xj)]2 (2)
其中:Ci为相速度的理论值,di(i=1,2,...,M)为面波相速度的实测值,xj(j=1,2,...,N)为近地表模型参数;
目标函数待到最小的充要条件是
根据式(3)可以求出
ATAΔX=AT·ΔC (4)
式(4)中的A为Jacobi矩阵,C为相速度,其具体形式为:
面波反演是对介质各层厚度和横波速度同时反演的过程,具体方法是:
(a)给定初值:根据实测频散曲线的变化规律,做出定性解释,即确定出可能的层数及各层的横波速度和厚度的变化范围,选择反演开始时各层横波速度和厚度的初值,也就是建立一个初始地层模型;
(b)将这个初始地层模型的理论频散曲线与实测瑞波频散曲线相比较,求得理论频散曲线与实测频散曲线之间的偏差ΔC(f);
(c)采用阻尼最小二乘法进行误差修正,调整初始地层模型的参数,重新计算调整后的地层模型的理论频散曲线,再与实测频散曲线比较,得到新的偏差ΔC(f);
(d)重复上述过程,反复调整地层模型参数,是理论频散曲线不断逼近实测频散曲线,直至偏差ΔC(f)小于给定的量值标准;
第四步:静校正量的计算
近地表静校正,使用替换速度静校正方法,利用直接反演的横波速度,对炮点和检波点分开求取静校正量。
进一步,在反演过程中,还包括将实验室测出的横波速度VS、纵波速度VP、密度ρ数据,根据公式
VP=F(VS)
ρ=G(VP)
反求出的函数关系引入到反演过程中的步骤,可减少未知数的个数,提高反演的精度。
本发明通过高精度的波场插值,联合应用纵波和转换波资料面波信息,并通过加权叠加,有效提高了直接反演横波速度的精度,有效的解决了常规方法通过间接求取横波静校正量而存在静校正不准的问题,从而有效提高多波资料的成像水平。
附图说明
图1是反演与理论频散曲线示意图。
图2a是罗家地区某区块三分量资料中的X分量资料图。
图2b罗家地区某区块三分量资料中的Y分量资料图。
图2c是罗家地区某区块三分量资料中的Z分量资料图。
图3a是罗家地区某区块X分量数据瑞雷面波频散曲线的提取图。
图3b是罗家地区某区块Y分量数据瑞雷面波频散曲线的提取图。
图3c是罗家地区某区块Z分量数据瑞雷面波频散曲线的提取图。
图4a是罗家地区某区块XZ分量叠加瑞雷面波频散曲线的提取图。
图4b是罗家地区某区块XZ分量加权叠加瑞雷面波频散曲线的提取图。
图5是罗家地区某区块反演的横波速度剖面图。
图6a是罗家地区某区块静校正前剖面图。
图6b是罗家地区某区块静校正后叠加剖面。
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
作为体纵波与横波的复合波,瑞雷面波的速度与纵波、横波速度有关,因此可以通过瑞雷面波反演近地表的横波速度结构。而瑞雷面波作为陆地地震勘探中的常见噪音,一直是被作为强线性干扰来压制的。如果能利用提取的面波实现转换波的静校正,就相当于“变废为宝”了。
一种基于表层横波速度反演的纵横波联合静校正方法,包括如下步骤:
第一步:波场高精度插值
在地震勘探中,因受野外施工条件、成本、效益等种因素的限制,野外采集道间距过大,容易产生空间假频。而利用地震道插值能加密空间采样率,防止偏移频散出现,提高信噪比。且由于面波具有速度低的特点,其数据往往会出现超周期现象,直接对这样的数据进行插值效果很差。而面波又有线性特征,首先对面波数据进行一定程度的线性动校(LNMO),然后进行小波变换插值,最后再反线性动校(RLNMO)恢复信号。
“小波”就是小区域、长度有限、均值为0的波形。所谓“小”是指它具有衰减性;而称之为“波”则是指它的波动性,其振幅正负相间的震荡形式。与Fourier变换相比,小波变换是时间(空间)频率的局部化分析,它通过伸缩平移运算对信号(函数)逐步进行多尺度细化,最终达到高频处时间细分,低频处频率细分,能自动适应时频信号分析的要求,从而可聚焦到信号的任意细节,解决了Fourier变换的困难问题,成为继Fourier变换以来在科学方法上的重大突破。小波变换是一种新的变换分析方法,它继承和发展了短时傅立叶变换局部化的思想,同时又克服了窗口大小不随频率变化等缺点,能够提供一个随频率改变的时间一频率窗口,是进行信号时频分析和处理的理想工具。
第二步:面波频散曲线提取
提取频散曲线可以通过两个线性变换实现:首先对原始数据做倾斜叠加(即线性Radon变换),将数据变换到了慢度-截距域(p-τ);其次,沿截距τ方向做一维傅里叶变换得到慢度-频率域(p-f);最后将慢度转换成速度即可获得速度-频率域(v-f)。
第三步:横波速度反演
横波速度的反演需要经过面波正演和反演两个方面,面波正演主要是根据弹性波动理论建立描述面波相速度C与波动频率f以及其它弹性参数之间关系的面波频散方程F(C,f),并解方程F(C,f)=0的根,即计算不同频率时的面波传播速度,由此了解特定介质模型中的面波频散特征。
建立的面波频散方程为:
其中:[J]是由传播矩阵表示的连乘形式,△n为纵波体积膨胀系数,Wn为横波角位移,C为相速度,μ0为介质的拉梅常数,w0为法向位移。
因为[J]中的元素是相速度C和常数R的函数,所以式(1)表示了C与R的关系,从而导出了相速度C的频散函数F(C,f)=0。
该方程是一个非线性方程,通过二分法求解得到不同频率下的速度值。最终得到面波速度随频率的变化——频散曲线。
面波反演的主要包括频散曲线C-f的求取和对频散曲线的解释,
反演过程中的目标函数为
Φ(X)=∑[di-Ci(xj)]2 (2)
其中Ci为相速度的理论值,di(i=1,2,...,M)为面波相速度的实测值,xj(j=1,2,...,N)为近地表模型参数。
目标函数待到最小的充要条件是
根据式(3)可以求出
ATAΔX=AT·ΔC (4)
式(4)中的A为Jacobi矩阵,C为相速度,其具体形式为:
一般而言,横波速度VS与纵波速度VP的关系、纵波速度VP与密度ρ的关系可以从实验室中得出,即
VP=F(VS)、ρ=G(VP)
在反演过程中如果知道F、G则可减少未知数的个数。另外,面波反演是对介质各层厚度和横波速度同时反演的过程,具体方法是:
a)给定初值:根据实测频散曲线的变化规律,做出定性解释,即确定出可能的层数及各层的横波速度和厚度的变化范围,选择反演开始时各层横波速度和厚度的初值,也就是建立一个初始地层模型;
(b)将这个初始地层模型的理论频散曲线与实测瑞波频散曲线相比较,求得理论频散曲线与实测频散曲线之间的偏差ΔC(f);
(c)采用阻尼最小二乘法进行误差修正,调整初始地层模型的参数,重新计算调整后的地层模型的理论频散曲线,再与实测频散曲线比较,得到新的偏差ΔC(f);
(d)重复上述过程,反复调整地层模型参数,是理论频散曲线不断逼近实测频散曲线,直至偏差ΔC(f)小于给定的量值标准。
第四步:静校正量的计算
近地表静校正,使用替换速度静校正方法,利用直接反演的横波速度,对炮点和检波点分开求取静校正量。实验验证
为对本方法这些程序应用了理论模型进行试验,检验方法的正确性、有效性和稳定性,对胜利油田罗家地区某区块的三维多分量资料进行了静校正处理,获得了较好的静校正效果。
图1是对理论模型反演计算得到的理论频散曲线与正演频散曲线进行对比,反演结果与模型参数吻合性较好,误差较小。
具体情况见下表
层数 | 反演厚度(m) | 反演横波速度(m/s) | 模型厚度(m) | 模型横波速度(m/s) |
1 | 10.14 | 99.67 | 10 | 90 |
2 | 100.87 | 148.75 | 100 | 160 |
3 | ∞ | 239.88 | ∞ | 240 |
对罗家三分量资料进行静校正处理,图2a、图2b、图2c为罗家多分量资料。如图3a、图3b、图3c所示,对其单分量数据提取频散曲线,可以看到谱能量不是很集中。如图4b所示,对其多分量数据纵波和转换波联合应用经过加权叠加提取的频散曲线,可以看到加权叠加后谱能量集中性明显变好,而如图4a所示,普通叠加的效果则没有加权叠加的效果好。如图5所示,为反演的罗家一条线的横波速度剖面。如图6a、图6b所示,通过静校正前后的叠加剖面对比,可以看到静校正后剖面同相轴的能量和连续性变好,取得了不错的静校正效果。
Claims (2)
1.一种基于表层横波速度反演的纵横波联合静校正方法,包括如下步骤:
第一步:波场高精度插值
所述波场高精度插值的具体方法是:
(a)对面波数据进行线性动校(LNMO);
(b)进行小波变换插值;
(c)反线性动校(RLNMO)恢复信号;
第二步:面波频散曲线提取
提取频散曲线可以通过两个线性变换实现,具体方法是:
(a)对原始数据做倾斜叠加,将数据变换到了慢度-截距域(p-τ);所述对原始数据做倾斜叠加采用线性Radon变换来完成;
(b)沿截距τ方向做一维傅里叶变换得到慢度-频率域(p-f);
(c)将慢度转换成速度即可获得速度-频率域;
第三步:横波速度反演
横波速度的反演需要经过面波正演和面波反演两个方面;
1)面波正演
面波正演主要是根据弹性波动理论建立描述面波相速度C与波动频率f以及其它弹性参数之间关系的面波频散方程F(C,f),并解方程F(C,f)=0的根,即计算不同频率时的面波传播速度,由此了解特定介质模型中的面波频散特征;
2)面波反演
建立的面波频散方程:
其中:[J]是由传播矩阵表示的连乘形式,△n为纵波体积膨胀系数,Wn为横波角位移,C为相速度,μ0为介质的拉梅常数,w0为法向位移;
因为[J]中的元素是相速度C和常数R的函数,所以式(1)表示了C与R的关系,从而导出了相速度C的频散函数
F(C,f)=0;
通过二分法求解得到不同频率下的速度值,最终得到面波速度随频率的变化——频散曲线;
面波反演的主要包括频散曲线C-f的求取和对频散曲线的解释,
反演过程中的目标函数为
Φ(X)=∑[di-Ci(xj)]2 (2)
其中:Ci为相速度的理论值,di(i=1,2,...,M)为面波相速度的实测值,xj(j=1,2,...,N)为近地表模型参数;
目标函数待到最小的充要条件是
根据式(3)可以求出
ATAΔX=AT·ΔC (4)
式(4)中的A为Jacobi矩阵,C为相速度,其具体形式为:
面波反演是对介质各层厚度和横波速度同时反演的过程,具体方法是:
(a)给定初值:根据实测频散曲线的变化规律,做出定性解释,即确定出可能的层数及各层的横波速度和厚度的变化范围,选择反演开始时各层横波速度和厚度的初值,也就是建立一个初始地层模型;
(b)将这个初始地层模型的理论频散曲线与实测瑞波频散曲线相比较,求得理论频散曲线与实测频散曲线之间的偏差ΔC(f);
(c)采用阻尼最小二乘法进行误差修正,调整初始地层模型的参数,重新计算调整后的地层模型的理论频散曲线,再与实测频散曲线比较,得到新的偏差ΔC(f);
(d)重复上述过程,反复调整地层模型参数,是理论频散曲线不断逼近实测频散曲线,直至偏差ΔC(f)小于给定的量值标准;
第四步:静校正量的计算
近地表静校正,使用替换速度静校正方法,利用直接反演的横波速度,对炮点和检波点分开求取静校正量。
2.如权利要求1所述的一种基于表层横波速度反演的纵横波联合静校正方法,其特征在于:
在反演过程中,还包括将实验室测出的横波速度VS、纵波速度VP、密度ρ数据,根据公式
VP=F(VS)
ρ=G(VP)
反求出的函数关系引入到反演过程中的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310520027.9A CN104730579B (zh) | 2013-12-18 | 2013-12-18 | 一种基于表层横波速度反演的纵横波联合静校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310520027.9A CN104730579B (zh) | 2013-12-18 | 2013-12-18 | 一种基于表层横波速度反演的纵横波联合静校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104730579A true CN104730579A (zh) | 2015-06-24 |
CN104730579B CN104730579B (zh) | 2018-02-13 |
Family
ID=53454650
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310520027.9A Active CN104730579B (zh) | 2013-12-18 | 2013-12-18 | 一种基于表层横波速度反演的纵横波联合静校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104730579B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106353805A (zh) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | 一种直观的微地震事件分类方法 |
CN106547024A (zh) * | 2015-09-16 | 2017-03-29 | 中国石油化工股份有限公司 | 针对微地震射孔数据的剩余静校正量估计方法和装置 |
CN106772575A (zh) * | 2016-11-28 | 2017-05-31 | 安徽理工大学 | 一种基于折射波与面波联合反演剩余煤层厚度的方法 |
CN107942387A (zh) * | 2017-10-23 | 2018-04-20 | 中国石油天然气集团公司 | 一种续至波衰减处理方法及装置 |
CN108363099A (zh) * | 2018-02-11 | 2018-08-03 | 中国石油化工股份有限公司 | 基于模拟退火和多尺度脊波分析的纵横波相关方法 |
WO2020029015A1 (zh) * | 2018-08-06 | 2020-02-13 | 南方科技大学 | 一种人工源面波勘探方法、面波勘探装置及终端设备 |
CN110879412A (zh) * | 2019-10-31 | 2020-03-13 | 南方科技大学 | 地下横波速度反演方法、装置、计算设备及存储介质 |
CN111856555A (zh) * | 2020-06-19 | 2020-10-30 | 同济大学 | 一种基于面波多尺度窗分析的地下探测方法 |
CN112394412A (zh) * | 2020-10-30 | 2021-02-23 | 中国石油天然气集团有限公司 | 一种基于浅层叠加建模的横波静校正方法及装置 |
CN112859172A (zh) * | 2021-03-05 | 2021-05-28 | 中国石油天然气集团有限公司 | 纵横波微测井资料处理方法及装置 |
CN112925025A (zh) * | 2021-01-26 | 2021-06-08 | 中国地质科学院地球物理地球化学勘查研究所 | 基于面波的地震数据静校正方法、装置及电子设备 |
CN113093279A (zh) * | 2020-01-08 | 2021-07-09 | 中国石油天然气集团有限公司 | 转换波静校正方法及装置 |
CN114185093A (zh) * | 2021-12-07 | 2022-03-15 | 中国石油大学(北京) | 一种基于瑞雷面波反演的近地表速度模型建立方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020075759A1 (en) * | 2001-01-21 | 2002-06-20 | Xixiang Zhou | Static correction method for exploration seismic data using first arrivals of seismic waves |
CN101907727A (zh) * | 2010-08-17 | 2010-12-08 | 中国科学院地质与地球物理研究所 | 一种面波多分量转换波静校正方法 |
US20110120724A1 (en) * | 2008-08-11 | 2011-05-26 | Krohn Christine E | Estimation of Soil Properties Using Waveforms of Seismic Surface Waves |
CN102749643A (zh) * | 2011-04-22 | 2012-10-24 | 中国石油天然气股份有限公司 | 一种波动方程正演的瑞利面波频散响应计算方法及其装置 |
-
2013
- 2013-12-18 CN CN201310520027.9A patent/CN104730579B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020075759A1 (en) * | 2001-01-21 | 2002-06-20 | Xixiang Zhou | Static correction method for exploration seismic data using first arrivals of seismic waves |
CN1128373C (zh) * | 2001-01-21 | 2003-11-19 | 周熙襄 | 用初至波对地震勘探资料进行短波长静校正处理的方法 |
US20110120724A1 (en) * | 2008-08-11 | 2011-05-26 | Krohn Christine E | Estimation of Soil Properties Using Waveforms of Seismic Surface Waves |
CN101907727A (zh) * | 2010-08-17 | 2010-12-08 | 中国科学院地质与地球物理研究所 | 一种面波多分量转换波静校正方法 |
CN102749643A (zh) * | 2011-04-22 | 2012-10-24 | 中国石油天然气股份有限公司 | 一种波动方程正演的瑞利面波频散响应计算方法及其装置 |
Non-Patent Citations (3)
Title |
---|
彭文 等: ""利用瑞雷面波调查表层结构的可行性"", 《天然气工业》 * |
李晶: ""面波在地震波场中的特性研究及其应用"", 《万方数据博士学位论文库》 * |
王东才: ""综合利用地质雷达和瞬态瑞利波法进行地表静校正"", 《物探装备》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106353805B (zh) * | 2015-07-17 | 2019-12-20 | 中国石油化工股份有限公司 | 一种直观的微地震事件分类方法 |
CN106353805A (zh) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | 一种直观的微地震事件分类方法 |
CN106547024A (zh) * | 2015-09-16 | 2017-03-29 | 中国石油化工股份有限公司 | 针对微地震射孔数据的剩余静校正量估计方法和装置 |
CN106547024B (zh) * | 2015-09-16 | 2018-08-07 | 中国石油化工股份有限公司 | 针对微地震射孔数据的剩余静校正量估计方法和装置 |
CN106772575A (zh) * | 2016-11-28 | 2017-05-31 | 安徽理工大学 | 一种基于折射波与面波联合反演剩余煤层厚度的方法 |
CN107942387A (zh) * | 2017-10-23 | 2018-04-20 | 中国石油天然气集团公司 | 一种续至波衰减处理方法及装置 |
CN108363099A (zh) * | 2018-02-11 | 2018-08-03 | 中国石油化工股份有限公司 | 基于模拟退火和多尺度脊波分析的纵横波相关方法 |
US11137511B2 (en) | 2018-08-06 | 2021-10-05 | Southern University Of Science And Technology | Active source surface wave prospecting method, surface wave exploration device and computer-readable storage medium |
WO2020029015A1 (zh) * | 2018-08-06 | 2020-02-13 | 南方科技大学 | 一种人工源面波勘探方法、面波勘探装置及终端设备 |
CN110879412A (zh) * | 2019-10-31 | 2020-03-13 | 南方科技大学 | 地下横波速度反演方法、装置、计算设备及存储介质 |
CN113093279A (zh) * | 2020-01-08 | 2021-07-09 | 中国石油天然气集团有限公司 | 转换波静校正方法及装置 |
CN113093279B (zh) * | 2020-01-08 | 2024-05-28 | 中国石油天然气集团有限公司 | 转换波静校正方法及装置 |
CN111856555A (zh) * | 2020-06-19 | 2020-10-30 | 同济大学 | 一种基于面波多尺度窗分析的地下探测方法 |
CN112394412A (zh) * | 2020-10-30 | 2021-02-23 | 中国石油天然气集团有限公司 | 一种基于浅层叠加建模的横波静校正方法及装置 |
CN112925025A (zh) * | 2021-01-26 | 2021-06-08 | 中国地质科学院地球物理地球化学勘查研究所 | 基于面波的地震数据静校正方法、装置及电子设备 |
CN112859172A (zh) * | 2021-03-05 | 2021-05-28 | 中国石油天然气集团有限公司 | 纵横波微测井资料处理方法及装置 |
CN114185093A (zh) * | 2021-12-07 | 2022-03-15 | 中国石油大学(北京) | 一种基于瑞雷面波反演的近地表速度模型建立方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104730579B (zh) | 2018-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104730579A (zh) | 一种基于表层横波速度反演的纵横波联合静校正方法 | |
US9128205B2 (en) | Process for creating image gathers | |
CN104570125B (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
CN101630014B (zh) | 一种利用垂直地震剖面数据对各向异性介质成像的方法 | |
CN101329405B (zh) | 一种简单的多参数地震反演方法 | |
CN102053261B (zh) | 一种地震数据处理方法 | |
CN103645503B (zh) | 一种三维时间域照明分析及振幅补偿方法 | |
CN102841379B (zh) | 一种基于共散射点道集的叠前时间偏移与速度分析方法 | |
CN104570079A (zh) | 一种纵波、转换横波地震资料的时间匹配方法 | |
CN104656142A (zh) | 一种利用垂直地震剖面与测井联合的地震层位标定方法 | |
CN103064115B (zh) | 一种射线参数域纵波与转换波匹配方法 | |
CN101598803B (zh) | 一种直接得到转换波叠加剖面的方法 | |
CN102073064B (zh) | 一种利用相位信息提高速度谱分辨率的方法 | |
CN104570110A (zh) | 一种基于纵横波匹配的多分量资料联合速度分析方法 | |
CN105093281A (zh) | 一种反演框架下的地震多波建模方法 | |
CN102213769A (zh) | 一种利用三维垂直地震剖面资料确定各向异性参数的方法 | |
CN104820242B (zh) | 一种面向叠前反演的道集振幅分频补偿方法 | |
CN104280777A (zh) | 一种压制陆上地震资料多次波干扰的方法 | |
CN103149588B (zh) | 一种利用井震标定计算vti各向异性参数的方法及系统 | |
CN103576200A (zh) | 低信噪比区浅层波阻抗界面静校正方法 | |
CN103645497A (zh) | 一种基于出射角的多分量波场分离方法 | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
CN105093301A (zh) | 共成像点反射角角道集的生成方法及装置 | |
CN104316966A (zh) | 一种流体识别方法及系统 | |
CN101609165A (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 |