CN115993650A - 一种基于棱柱波的地震干涉成像方法 - Google Patents
一种基于棱柱波的地震干涉成像方法 Download PDFInfo
- Publication number
- CN115993650A CN115993650A CN202310279095.4A CN202310279095A CN115993650A CN 115993650 A CN115993650 A CN 115993650A CN 202310279095 A CN202310279095 A CN 202310279095A CN 115993650 A CN115993650 A CN 115993650A
- Authority
- CN
- China
- Prior art keywords
- wave
- prismatic
- excitation
- interference
- green
- 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
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于棱柱波的地震干涉成像方法，涉及石油勘探技术领域，包括如下步骤：输入速度场、偏移结果及基准面，建立观测系统；生成干涉需要的三种格林函数；推导双重格林函数干涉方程；输出校准后的格林函数并作逆傅里叶变换，得到炮记录；用得到的炮记录做逆时偏移。本发明采用了基于互相关原理的地震干涉方程，实现两个格林函数之间相同部分的消去处理，首次将地震干涉技术与棱柱波相结合，能够实现声波介质下的观测系统下移，观测系统下移后的炮记录只包含基准面下方的介质信息，不受基准面上覆介质的影响，因此相较于原始的炮记录而言，校准后得到的炮记录的有效长度小，进而后期逆时偏移成像的运算成本降低。
Description
技术领域
本发明涉及石油勘探技术领域，尤其涉及一种基于棱柱波的地震干涉成像方法。
背景技术
传统的地震干涉技术是针对垂直井间系统设计的，而实际操作中这种采集系统占比很少，因此有必要提出适用于地表观测系统的地震干涉技术。现有的适用于地表观测系统的地震干涉技术具有众多约束，例如靶区面积远小于实际工具面积、上覆介质水平均匀等，因此有必要克服上述束缚，研发出适用于复杂模型的地震干涉技术。
发明内容
为解决上述技术问题，本发明公开了一种基于棱柱波的地震干涉成像方法。
为实现上述目的，本发明采用下述技术方案：
一种基于棱柱波的地震干涉成像方法，包括如下步骤：
步骤S1、输入速度场、偏移结果及基准面，建立观测系统；
步骤S2、生成干涉需要的三种格林函数；
步骤S3、推导双重格林函数干涉方程，得到校准后的格林函数；
步骤S4、输出校准后的格林函数并作逆傅里叶变换，得到炮记录；
步骤S5、用步骤S4得到的炮记录做逆时偏移。
可选地，步骤S2中，生成干涉需要的三种格林函数的步骤，包括：
不同观测系统下地震波的一次波炮记录，如式（1）~式（3）所示：
（1）
（2）
（3）
不同观测系统下地震波的棱柱波炮记录，如式（4）~式（6）所示：
（4）
（5）
（6）
不同观测系统下地震波的棱柱波炮记录用格林函数G表达，如式（7）~式（9）所示：
（7）
（8）
（9）
式中，
G表示格林函数，
X，
Y，
A，
B表示坐标位置，
top表示计算范围是基准面上侧的模型区域，
all表示计算范围是全模型区域，
S 0 表示地表处激发震源，
S 1 表示地下基准面处激发震源，表示一次波，表示棱柱波，
v表示速度，
s表示震源，表示梯度算子，
ω表示角频率，表示全模型偏移结果，表示基准面上覆介质偏移结果。
可选地，步骤S3中，推导双重格林函数干涉方程，得到校准后的格林函数的步骤包括：
步骤S31、X激发、Y接收的棱柱波炮记录与X激发、B接收的棱柱波炮记录互相关处理，前者传播路径为X-B-Y，后者传播路径为X-B，两者互相关处理消去X-B，得到B激发、Y接收的棱柱波炮记录，如式（10）所示：
（10）
步骤S32、B激发、Y接收的棱柱波与A激发、Y接收的棱柱波炮记录互相关处理，前者传播路径为B-A-Y，后者传播路径为A-Y，两者互相关处理消去A-Y，得到B激发、A接收的棱柱波炮记录，如式（11）所示：
（11）
式中，
G表示格林函数，
X，
Y，
A，
B表示坐标位置，表示取共轭，
m表示地表处检波点个数，
n表示地下基准面处检波点个数，
i表示虚部单位，
ω表示角频率。
可选地，步骤S4中，输出校准后的格林函数并作逆傅里叶变换，得到炮记录的步骤，包括：
将棱柱波炮记录从频率域中变换到时间域中，如式（12）所示：
（12）
式中，
G表示格林函数，
A，
B表示坐标位置，
t 0表示激发时间，
t max表示人为选择的有效信号时长，
ifft表示逆傅里叶变换，
output表示输出。
可选地，步骤S5中，用步骤S4得到的炮记录做逆时偏移的步骤，包括：
成像采用互相关条件，如式（13）所示：
（13）
式中，
I表示成像结果，
S表示正传波场，
R表示反传波场，
x，
z表示坐标位置，
t表示时间，
T表示总接收时间。
本发明的有益效果是，
1、本发明的方法采用了基于互相关原理的地震干涉方程，可以实现两个格林函数之间相同部分的消去处理。本发明首次将地震干涉技术与棱柱波相结合，能够实现声波介质下的观测系统下移，观测系统下移后的炮记录只包含基准面下方的介质信息，输出的地下基准面激发、地下基准面接收的炮记录，不受基准面上覆介质的影响，因此相较于原始的炮记录而言，校准后得到的炮记录的有效长度小，进而后期逆时偏移成像的运算成本降低，其逆时偏移成像结果优于常规全局模型的对应部分。
2、本发明克服了传统干涉技术对靶区面积的要求，突破了传统干涉只适用于垂直井间观测系统的束缚，且无需做炮记录的道集转换，进一步提升运算效率。本发明不但适用于简单模型，也适用于复杂模型。
附图说明
图1为本发明一种基于棱柱波的地震干涉成像方法的流程图；
图2为本发明一实施例示出的地表激发地表接收的观测系统；
图3为本发明一实施例示出的地表激发地下基准面接收的观测系统；
图4为本发明一实施例示出的地下基准面激发地表接收的观测系统；
图5为本发明一实施例示出的地下基准面激发地下基准面接收的观测系统；
图6为本发明一实施例示出的基于互相关原理的双重格林函数干涉技术示意图；
图7为本发明一实施例示出的两层水平模型；
图8为本发明图7的输入炮记录，其中，（a）为图2所示观测系统下得到的炮记录；（b）为图3所示观测系统下得到的炮记录；（c）为图4所示观测系统下得到的炮记录；
图9为本发明图7的干涉效果炮记录对比，其中，（a）为重校准得到的炮记录；（b）为真实情况下的炮记录；
图10为本发明图7的干涉效果逆时偏移结果对比，（a）为图9（a）的干涉效果逆时偏移结果；（b）图8（a）的成像结果只取基准面下的部分；
图11为本发明一实施例示出的盐丘模型；
图12为本发明图11的输入炮记录，其中，（a）为图2所示观测系统下得到的炮记录，（b）为图3所示观测系统下得到的炮记录，（c）为图4所示观测系统下得到的炮记录；
图13为本发明图11的干涉效果炮记录对比，其中，（a）为重校准得到的炮记录，（b）为真实情况下的炮记录；
图14为本发明图11的干涉效果逆时偏移结果对比，（a）为图13（a）的干涉效果逆时偏移结果；（b）为图12（a）的成像结果只取基准面下的部分。
实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
一种基于棱柱波的地震干涉成像方法，如图1所示，包括如下步骤：
步骤S1、输入速度场、偏移结果及基准面，建立观测系统；
步骤S2、生成干涉需要的三种格林函数；
步骤S3、推导双重格林函数干涉方程，得到校准后的格林函数；
步骤S4、输出校准后的格林函数并作逆傅里叶变换，得到炮记录；
步骤S5、用步骤S4得到的炮记录做逆时偏移。
在该可选实施例中，炮记录采用交错网格、吸收边界条件，精度为时间二阶空间十阶的有限差分算法计算得到。
可选地，步骤S2中，生成干涉需要的三种格林函数的步骤，包括：
不同观测系统下地震波的一次波炮记录，如式（1）~式（3）所示：
（1）
（2）
（3）
式（1）表达了图2所示的观测系统，该系统是在整个速度模型中计算，震源和检波点都放置在地表，得到的是整个模型的一次波炮记录；式（2）表达了图3所示的观测系统，该系统是在基准面以上的速度模型中计算，震源放置在地表，检波点放置在地下基准面，得到的是基准面以上的模型的一次波炮记录；式（3）表达了图4所示的观测系统，与式（2）相同的是，该系统也是在基准面以上的速度模型中计算，得到的是基准面以上的模型的另一种传播路径下的一次波炮记录，与式（2）不同的是，该系统震源放置在地下基准面，检波点放置在地表；
不同观测系统下地震波的棱柱波炮记录，如式（4）~式（6）所示：
（4）
（5）
（6）
式（4）表达了图2所示的观测系统，该系统是在整个速度模型中计算，震源和检波点都放置在地表，得到的是整个模型的棱柱波炮记录；式（5）表达了图3所示的观测系统，该系统是在基准面以上的速度模型中计算，震源放置在地表，检波点放置在地下基准面，得到的是基准面以上的模型的棱柱波炮记录；式（6）表达了图4所示的观测系统，与式（5）相同的是，该系统也是在基准面以上的速度模型中计算，得到的是基准面以上的模型的另一种传播路径下的棱柱波炮记录，与式（5）不同的是，该系统震源放置在地下基准面，检波点放置在地表；
为了方便表达，不同观测系统下地震波的棱柱波炮记录用格林函数G表达，如式（7）~式（9）所示：
（7）
（8）
（9）
式中，
G表示格林函数，
X，
Y，
A，
B表示坐标位置，
top表示计算范围是基准面上侧的模型区域，
all表示计算范围是全模型区域，
S 0 表示地表处激发震源，
S 1 表示地下基准面处激发震源，表示一次波，表示棱柱波，
v表示速度，
s表示震源，表示梯度算子，
ω表示角频率，表示全模型偏移结果，表示基准面上覆介质偏移结果。
在该可选实施例中，式（7）~式（9）分别对应了式（4）~式（6）。
通过地震波有限差分正演模拟可以得到式（7）~式（9）所示的棱柱波炮记录，接下来便是不同观测系统下的棱柱波炮记录的互相关处理。
可选地，步骤S3中，推导双重格林函数干涉方程，得到校准后的格林函数的步骤包括：
步骤S31、X激发、Y接收的棱柱波炮记录与X激发、B接收的棱柱波炮记录互相关处理，前者传播路径为X-B-Y，后者传播路径为X-B，两者互相关处理消去X-B，得到B激发、Y接收的棱柱波炮记录，如式（10）所示：
（10）
步骤S32、B激发、Y接收的棱柱波与A激发、Y接收的棱柱波炮记录互相关处理，前者传播路径为B-A-Y，后者传播路径为A-Y，两者互相关处理消去A-Y，得到B激发、A接收的棱柱波炮记录，如式（11）所示：
（11）
式中，
G表示格林函数，
X，
Y，
A，
B表示坐标位置，表示取共轭，
m表示地表处检波点个数，
n表示地下基准面处检波点个数，
i表示虚部单位，
ω表示角频率。
上述过程均在频率域中完成，后续处理需要把棱柱波炮记录从频率域中变换到时间域中。
可选地，步骤S4中，输出校准后的格林函数并作逆傅里叶变换，得到炮记录的步骤，包括：
将棱柱波炮记录从频率域中变换到时间域中，如式（12）所示：
（12）
式中，
G表示格林函数，
A，
B表示坐标位置，
t 0表示激发时间，
t max表示人为选择的有效信号时长，
ifft表示逆傅里叶变换，
output表示输出；
需要注意的是，经过两次互相关处理得到的棱柱波炮记录依然与原始的棱柱波炮记录长度一致，但是校准后的炮记录只含基准面下侧模型的信息，所以校准后的炮记录下侧零值（无效值）很多，该实施例只截取上侧有效信号进行后期的成像处理。
基于式（1）~式（12）得到的棱柱波炮记录是来自基准面下侧模型，因此研究区域被锁定在基准面下侧模型，基准面上侧模型已不再被研究，因此，该实施例逆时偏移只在基准面下侧模型中完成，炮记录变短和计算模型的纵向尺度变小均减少了运算时间。
可选地，步骤S5中，用步骤S4得到的炮记录做逆时偏移的步骤，包括：
成像采用互相关条件，如式（13）所示：
（13）
式中，
I表示成像结果，
S表示正传波场，
R表示反传波场，
x，
z表示坐标位置，
t表示时间，
T表示总接收时间。
应用实验
本发明一种基于棱柱波的地震干涉成像方法，应用于两个模型，取得了理想的成像效果。图2示出了地表激发、地表接收的观测系统，这是常规的观测系统，也是本实施例需要处理的观测系统。图3示出了地表激发、地下基准面接收的观测系统，它生成了本实施例第一步想要消去的地震波传播路径的炮记录。图4示出了地下基准面激发、地表接收的观测系统，它生成了本实施例第二步想要消去的地震波传播路径的炮记录。图5示出了地下基准面激发、地下基准面接收的观测系统，这也就是本实施例想要得到的观测系统，即通过本发明实施例提出的算法将图2所示的观测系统产生的炮记录转换为图5所示的观测系统产生的炮记录。
图6示出了基于互相关原理的双重格林函数干涉方法示意图，即：本方法的目标是输入图2~图4所示观测系统下得到炮记录，输出图5所示观测系统下得到炮记录。第一步在X-B-A-Y这一传播路径下消去X-B，得到B-A-Y这一传播路径；第二步在B-A-Y这一传播路径消去A-Y，得到最终想要的B-A这一传播路径。
先对简单模型进行测试，如图7所示的两层水平模型，模型横向网格点数为251，纵向网格点数为251，基准面选在纵向深度1500米处。
图8示出了图7模型下的输入炮记录，其中（a）为图2所示观测系统下得到的炮记录，（b）为图3所示观测系统下得到的炮记录，（c）为图4所示观测系统下得到的炮记录。
图9示出了干涉效果炮记录对比图，其中，（a）为重校准得到炮记录，（b）为真实情况下炮记录，通过对比可以看出本发明方法得到炮记录具有真实炮记录的特征。
图10示出了干涉效果逆时偏移对比图，其中，（a）为图9（a）的干涉效果逆时偏移结果，（b）为图8（a）的成像结果只取基准面下的部分，通过对比可以看出本发明方法得到炮记录的逆时偏移结果与全局速度模型逆时偏移结果的基准面下的部分一样。
再对复杂模型进行测试，如图11所示的盐丘模型，模型横向网格点数为641，纵向网格点数为170，基准面选在纵向深度550米处。
图12示出了图11模型下的输入炮记录，其中（a）为图2所示观测系统下得到的炮记录，（b）为图3所示观测系统下得到的炮记录，（c）为图4所示观测系统下得到的炮记录。
图13示出了干涉效果炮记录对比图，其中，（a）为重校准得到炮记录，（b）为真实情况下炮记录，通过对比可以看出本方法得到炮记录具有真实炮记录的特征。
图14示出了干涉效果逆时偏移对比图，其中，（a）为图13（a）的干涉效果逆时偏移结果，（b）为图12（a）的成像结果只取基准面下的部分，通过对比可以看出本方法得到炮记录的逆时偏移结果与全局速度模型逆时偏移结果的基准面下的部分一样，且盐丘下侧（箭头指示）成像结果更好。
本发明提出的一种基于棱柱波的地震干涉成像方法可以实现声波介质下的观测系统下移处理，即将观测系统从地表下移到地下基准面，对于复杂模型可以提高基准面下侧的成像精度。本发明提出方法的目标是实现不同观测系统之间的炮记录转换，这一过程是在频率域中完成，即为重校准，在得到重校准的炮记录后，在时间域里完成后续的成像工作。
当然，上述说明并非是对本发明的限制，本发明也并不仅限于上述举例，本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换，也应属于本发明的保护范围。
Claims (5)
1.一种基于棱柱波的地震干涉成像方法，其特征在于，包括如下步骤：
步骤S1、输入速度场、偏移结果及基准面，建立观测系统；
步骤S2、生成干涉需要的三种格林函数；
步骤S3、推导双重格林函数干涉方程，得到校准后的格林函数；
步骤S4、输出校准后的格林函数并作逆傅里叶变换，得到炮记录；
步骤S5、用步骤S4得到的炮记录做逆时偏移。
2.如权利要求1所述的一种基于棱柱波的地震干涉成像方法，其特征在于，步骤S2中，生成干涉需要的三种格林函数的步骤，包括：
不同观测系统下地震波的一次波炮记录，如式（1）~式（3）所示：
（1）
（2）
（3）
不同观测系统下地震波的棱柱波炮记录，如式（4）~式（6）所示：
（4）
（5）
（6）
不同观测系统下地震波的棱柱波炮记录用格林函数G表达，如式（7）~式（9）所示：
（7）
（8）
（9）
式中，G表示格林函数，X，Y，A，B表示坐标位置，top表示计算范围是基准面上侧的模型区域，all表示计算范围是全模型区域，S 0 表示地表处激发震源，S 1 表示地下基准面处激发震源，表示一次波，表示棱柱波，v表示速度，s表示震源，表示梯度算子，ω表示角频率，表示全模型偏移结果，表示基准面上覆介质偏移结果。
3.如权利要求1所述的一种基于棱柱波的地震干涉成像方法，其特征在于，步骤S3中，推导双重格林函数干涉方程，得到校准后的格林函数的步骤包括：
步骤S31、X激发、Y接收的棱柱波炮记录与X激发、B接收的棱柱波炮记录互相关处理，前者传播路径为X-B-Y，后者传播路径为X-B，两者互相关处理消去X-B，得到B激发、Y接收的棱柱波炮记录，如式（10）所示：
（10）
步骤S32、B激发、Y接收的棱柱波与A激发、Y接收的棱柱波炮记录互相关处理，前者传播路径为B-A-Y，后者传播路径为A-Y，两者互相关处理消去A-Y，得到B激发、A接收的棱柱波炮记录，如式（11）所示：
（11）
式中，G表示格林函数，X，Y，A，B表示坐标位置，表示取共轭，m表示地表处检波点个数，n表示地下基准面处检波点个数，i表示虚部单位，ω表示角频率。
4.如权利要求1所述的一种基于棱柱波的地震干涉成像方法，其特征在于，步骤S4中，输出校准后的格林函数并作逆傅里叶变换，得到炮记录的步骤，包括：
将棱柱波炮记录从频率域中变换到时间域中，如式（12）所示：
（12）
式中，G表示格林函数，A，B表示坐标位置，t 0表示激发时间，t max表示人为选择的有效信号时长，ifft表示逆傅里叶变换，output表示输出。
5.如权利要求1所述的一种基于棱柱波的地震干涉成像方法，其特征在于，步骤S5中，用步骤S4得到的炮记录做逆时偏移的步骤，包括：
成像采用互相关条件，如式（13）所示：
（13）
式中，I表示成像结果，S表示正传波场，R表示反传波场，x，z表示坐标位置，t表示时间，T表示总接收时间。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310279095.4A CN115993650B (zh) | 2023-03-22 | 2023-03-22 | 一种基于棱柱波的地震干涉成像方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310279095.4A CN115993650B (zh) | 2023-03-22 | 2023-03-22 | 一种基于棱柱波的地震干涉成像方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115993650A true CN115993650A (zh) | 2023-04-21 |
| CN115993650B CN115993650B (zh) | 2023-06-06 |
Family
ID=85993676
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310279095.4A Active CN115993650B (zh) | 2023-03-22 | 2023-03-22 | 一种基于棱柱波的地震干涉成像方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115993650B (zh) |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100256916A1 (en) * | 2009-04-03 | 2010-10-07 | Chevron U.S.A. Inc. | Method for target-oriented reverse time migration for prestack depth imaging |
| CN102879816A (zh) * | 2012-07-17 | 2013-01-16 | 中国科学院地质与地球物理研究所 | 一种地震多次波偏移方法 |
| KR101549388B1 (ko) * | 2014-10-17 | 2015-09-02 | 한국지질자원연구원 | 탄성파 다성분 자료에 대한 중합전 ｅｇｓ 구조보정 방법 |
| CN105676277A (zh) * | 2015-12-30 | 2016-06-15 | 中国石油大学(华东) | 一种提高高陡构造速度反演效率的全波形联合反演方法 |
| CN106526677A (zh) * | 2016-10-26 | 2017-03-22 | 中海石油(中国)有限公司 | 一种海上自适应压制鬼波的宽频逆时偏移成像方法 |
| US20170363759A1 (en) * | 2016-06-17 | 2017-12-21 | Cgg Services Sa | System and method for seismic interferometry optimized data acquisition |
| CN108037526A (zh) * | 2017-11-23 | 2018-05-15 | 中国石油大学(华东) | 基于全波波场vsp/rvsp地震资料的逆时偏移方法 |
| CN110058303A (zh) * | 2019-05-06 | 2019-07-26 | 吉林大学 | 声波各向异性逆时偏移混合方法 |
| CN110376646A (zh) * | 2019-07-19 | 2019-10-25 | 中国石油大学(华东) | 一种基于曲坐标系纵横波解构方程的弹性棱柱波逆时偏移成像方法 |
| CN110488354A (zh) * | 2019-07-19 | 2019-11-22 | 中国石油大学(华东) | 一种q补偿的起伏地表棱柱波与一次波联合最小二乘逆时偏移成像方法 |
| CN111158049A (zh) * | 2019-12-27 | 2020-05-15 | 同济大学 | 一种基于散射积分法的地震逆时偏移成像方法 |
| CN111505718A (zh) * | 2020-04-28 | 2020-08-07 | 西安交通大学 | 一种高分辨率地下结构保幅成像方法 |
| CN112083493A (zh) * | 2020-08-19 | 2020-12-15 | 中国石油大学(华东) | 一种三维c-τ坐标系的圆锥波编码多震源最小二乘逆时偏移成像方法 |
| CN114942472A (zh) * | 2022-03-29 | 2022-08-26 | 中国石油大学(华东) | 一种基于上行射线追踪策略的偏移成像方法及其设备 |
-
2023
- 2023-03-22 CN CN202310279095.4A patent/CN115993650B/zh active Active
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100256916A1 (en) * | 2009-04-03 | 2010-10-07 | Chevron U.S.A. Inc. | Method for target-oriented reverse time migration for prestack depth imaging |
| CN102879816A (zh) * | 2012-07-17 | 2013-01-16 | 中国科学院地质与地球物理研究所 | 一种地震多次波偏移方法 |
| KR101549388B1 (ko) * | 2014-10-17 | 2015-09-02 | 한국지질자원연구원 | 탄성파 다성분 자료에 대한 중합전 ｅｇｓ 구조보정 방법 |
| CN105676277A (zh) * | 2015-12-30 | 2016-06-15 | 中国石油大学(华东) | 一种提高高陡构造速度反演效率的全波形联合反演方法 |
| US20170363759A1 (en) * | 2016-06-17 | 2017-12-21 | Cgg Services Sa | System and method for seismic interferometry optimized data acquisition |
| CN106526677A (zh) * | 2016-10-26 | 2017-03-22 | 中海石油(中国)有限公司 | 一种海上自适应压制鬼波的宽频逆时偏移成像方法 |
| CN108037526A (zh) * | 2017-11-23 | 2018-05-15 | 中国石油大学(华东) | 基于全波波场vsp/rvsp地震资料的逆时偏移方法 |
| CN110058303A (zh) * | 2019-05-06 | 2019-07-26 | 吉林大学 | 声波各向异性逆时偏移混合方法 |
| CN110376646A (zh) * | 2019-07-19 | 2019-10-25 | 中国石油大学(华东) | 一种基于曲坐标系纵横波解构方程的弹性棱柱波逆时偏移成像方法 |
| CN110488354A (zh) * | 2019-07-19 | 2019-11-22 | 中国石油大学(华东) | 一种q补偿的起伏地表棱柱波与一次波联合最小二乘逆时偏移成像方法 |
| CN111158049A (zh) * | 2019-12-27 | 2020-05-15 | 同济大学 | 一种基于散射积分法的地震逆时偏移成像方法 |
| CN111505718A (zh) * | 2020-04-28 | 2020-08-07 | 西安交通大学 | 一种高分辨率地下结构保幅成像方法 |
| CN112083493A (zh) * | 2020-08-19 | 2020-12-15 | 中国石油大学(华东) | 一种三维c-τ坐标系的圆锥波编码多震源最小二乘逆时偏移成像方法 |
| CN114942472A (zh) * | 2022-03-29 | 2022-08-26 | 中国石油大学(华东) | 一种基于上行射线追踪策略的偏移成像方法及其设备 |
Non-Patent Citations (3)
| Title |
|---|
| 张宇;: "从成像到反演:叠前深度偏移的理论、实践与发展", 石油物探, no. 01 * |
| 王娟;李振春;陶丽;: "逆时偏移成像条件研究", 地球物理学进展, no. 03 * |
| 谷丙洛;韩建光;陈军;黄建平;李振春;: "碳酸盐岩油气储层高精度成像方法研究", 地质学报, no. 08 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115993650B (zh) | 2023-06-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102395902B (zh) | 使用快速面向目标照明计算的地震成像系统及方法 | |
| KR101548976B1 (ko) | 지진 표면파들의 파형들을 사용하는 토양 특성들의 추정 | |
| CN111158049B (zh) | 一种基于散射积分法的地震逆时偏移成像方法 | |
| CN102695970B (zh) | 用于表征石油或者天然气储集层随时间演变的改进方法 | |
| CN108241173B (zh) | 一种地震资料偏移成像方法及系统 | |
| CN112593922B (zh) | 一种阵列声波测井评价固井二界面胶结质量的方法及装置 | |
| CN109738952B (zh) | 基于全波形反演驱动的被动源直接偏移成像方法 | |
| CN111290017B (zh) | 一种震电波场联合提取瑞雷波频散特征的面波勘探方法 | |
| CN112285767B (zh) | 海底地震仪四分量海洋面波多阶频散能量成像装置及方法 | |
| EA032186B1 (ru) | Сейсмическая адаптивная фокусировка | |
| CN112327362B (zh) | 速度域的海底多次波预测与追踪衰减方法 | |
| CN109655918B (zh) | 地面浅井微地震监测观测台站位置确定方法及系统 | |
| CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
| CN110579795A (zh) | 基于被动源地震波形及其逆时成像的联合速度反演方法 | |
| CN109239776B (zh) | 一种地震波传播正演模拟方法和装置 | |
| CN111665563B (zh) | 基于聚焦分析的叠前偏移垂向分辨率评价方法 | |
| CN115993650B (zh) | 一种基于棱柱波的地震干涉成像方法 | |
| CN110907989A (zh) | 重建拟地面地震反射波成像方法及系统 | |
| CN112630830B (zh) | 一种基于高斯加权的反射波全波形反演方法及系统 | |
| CN110579799B (zh) | 一种等旅行时间间隔的地震采集观测方法及系统 | |
| CN111175822B (zh) | 改进直接包络反演与扰动分解的强散射介质反演方法 | |
| CN113447981B (zh) | 一种基于共成像点道集的反射全波形反演方法 | |
| CN112305615B (zh) | 一种地震资料角度域共成像点道集提取方法及系统 | |
| CN114942472A (zh) | 一种基于上行射线追踪策略的偏移成像方法及其设备 | |
| WO2002023222A1 (en) | Illumination weighted imaging condition for migrated seismic data |
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 |