CN108897046B - 远探测声波逆时偏移成像方法及装置 - Google Patents
远探测声波逆时偏移成像方法及装置 Download PDFInfo
- Publication number
- CN108897046B CN108897046B CN201810966059.4A CN201810966059A CN108897046B CN 108897046 B CN108897046 B CN 108897046B CN 201810966059 A CN201810966059 A CN 201810966059A CN 108897046 B CN108897046 B CN 108897046B
- Authority
- CN
- China
- Prior art keywords
- time
- wave
- wave field
- reverse
- excitation
- 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
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 102
- 238000003384 imaging method Methods 0.000 title claims abstract description 94
- 230000005012 migration Effects 0.000 title claims abstract description 82
- 238000013508 migration Methods 0.000 title claims abstract description 82
- 230000002441 reversible effect Effects 0.000 title claims abstract description 69
- 230000005284 excitation Effects 0.000 claims abstract description 71
- 238000000034 method Methods 0.000 claims abstract description 49
- 230000008569 process Effects 0.000 claims abstract description 21
- 238000007781 pre-processing Methods 0.000 claims abstract description 10
- 238000004590 computer program Methods 0.000 claims description 13
- 238000001914 filtration Methods 0.000 claims description 8
- 238000005381 potential energy Methods 0.000 claims description 5
- 238000013075 data extraction Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 14
- 238000013500 data storage Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 17
- 238000005516 engineering process Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 5
- 238000005070 sampling Methods 0.000 description 5
- 239000000523 sample Substances 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000005669 field effect Effects 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000005404 monopole Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本申请实施例提供了一种远探测声波逆时偏移成像方法及装置，该方法包括：获取远探测声波测井数据并进行预处理；从预处理后的远探测声波测井数据中提取纵横波时差及反射波；根据纵横波时差构建偏移速度模型；基于偏移速度模型进行第一次震源波场正向延拓，以在正向延拓过程中利用激发振幅成像条件获得激发时间；基于偏移速度模型进行第二次震源波场正向延拓，并保存以激发时间为中心的指定时窗内的震源波场；基于反射波进行检波点波场逆向延拓，并保存以激发时间为中心的指定时窗内的检波点波场；对指定时窗内的震源波场及检波点波场进行互相关成像。本申请实施例可以减少远探测声波逆时偏移成像的数据存储，从而提高处理效率。
Description
技术领域
本申请涉及远探测声波测井技术领域，尤其是涉及一种远探测声波逆时偏移成像方法及装置。
背景技术
常规声波测井利用井筒滑行波获取井旁地层信息，无法探知井壁纵深地带缝洞体发育情况。远探测声波测井是近年来发展起来的新一代声波测井技术，该技术借鉴地面地震勘探思想，打破了常规声波测井技术利用井筒波的模式思维，创造性地利用反射波探测井旁缝洞等隐蔽储层构造，将探测距离从井周1-3米扩展至井旁几十米，填补了测井和地震间的勘探空白，对于复杂储层油气勘探具有重要意义。
大量研究实践表明，远探测声波测井资料核心处理方法主要分为两部分：反射波提取和偏移成像。偏移成像源自地震勘探，是将地面接收到的发射信号恢复成地下几何空间构造的一类技术。先后有学者将基于射线理论的Kirchhoff偏移和基于单程波理论的FK偏移引入到远探测声波测井偏移成像中来，实现了对井壁外围缝洞异常体的成像，但其成像精度差且多解性较强，油田现场效果较差。逆时偏移基于双程波动方程，较其它偏移方法其成像精度更高。但由于逆时偏移需要对应时刻的正传和反传波场互相关成像，而正传和反传在时间上反向，其不可避免要将所有时刻的正传波场存储下来，巨大的数据存储和读写(尤其对于三维数据)成为制约逆时偏移在实际生产中高效应用的重要瓶颈。
发明内容
本申请实施例的目的在于提供一种远探测声波逆时偏移成像方法及装置，以提高远探测声波逆时偏移成像的处理效率。
为达到上述目的，一方面，本申请实施例提供了一种远探测声波逆时偏移成像方法，包括：
获取远探测声波测井数据并进行预处理；
从预处理后的远探测声波测井数据中提取纵横波时差及反射波；
根据所述纵横波时差构建偏移速度模型；
基于所述偏移速度模型进行第一次震源波场正向延拓，以在正向延拓过程中利用激发振幅成像条件获得激发时间；
基于所述偏移速度模型进行第二次震源波场正向延拓，并保存以所述激发时间为中心的指定时窗内的震源波场；
基于所述反射波进行检波点波场逆向延拓，并保存以所述激发时间为中心的指定时窗内的检波点波场；
对指定时窗内的震源波场及检波点波场进行互相关成像。
本申请实施例的远探测声波逆时偏移成像方法，所述震源波场正向延拓基于以下时间正向延拓有限差分公式实现：
其中，n表示变量在时间网格的位置；lx和lz分别表示变量在空间网格x方向及z方向的位置；Δt表示时间间隔；Δx和Δz分别表示空间网格x方向及z方向的间隔；δx和δz分别表示空间x方向及z方向的差分；2N表示空间差分阶数；vx和vz分别表示震源波场速度变量的x分量及z分量；σxx和σzz分别表示震源波场应力变量的xx分量及zz分量；ρ表示能量密度；λ表示拉梅参数。
本申请实施例的远探测声波逆时偏移成像方法，所述在正向延拓过程中利用激发振幅成像条件获得激发时间，包括：
根据公式
确定激发振幅，并确定所述激发振幅对应的激发时间；
其中，E、K和W分别表示震源波场中的总能量密度、动能密度和势能密度；ρ表示能量密度；vi表示速度变量，σij表示应力变量，εij表示应变变量。
本申请实施例的远探测声波逆时偏移成像方法，所述检波点波场逆向延拓基于以下时间逆向延拓有限差分公式实现：
其中，n表示变量在时间网格的位置；lx和lz分别表示变量在空间网格x方向及z方向的位置；Δt表示时间间隔；Δx和Δz分别表示空间网格x方向及z方向的间隔；δx和δz分别表示空间x方向及z方向的差分；2N表示空间差分阶数；vx和vz分别表示震源波场速度变量的x分量及z分量；σxx和σzz分别表示震源波场应力变量的xx分量及zz分量；ρ表示能量密度；λ表示拉梅参数。
本申请实施例的远探测声波逆时偏移成像方法，所述从预处理后的远探测声波测井数据中提取反射波，包括：
根据中值滤波、FK滤波或多尺度相关对比法从预处理后的远探测声波测井数据中提取反射波；
根据倾斜中值滤波滤除所述反射波中的地层无效反射波。
另一方面，本申请实施例还提供了一种远探测声波逆时偏移成像装置，包括：
数据获取模块，用于获取远探测声波测井数据并进行预处理；
数据提取模块，用于从预处理后的远探测声波测井数据中提取纵横波时差及反射波；
模型构建模块，用于根据所述纵横波时差构建偏移速度模型；
第一正向延拓模块，用于基于所述偏移速度模型进行第一次震源波场正向延拓，以在正向延拓过程中利用激发振幅成像条件获得激发时间；
第二正向延拓模块，用于基于所述偏移速度模型进行第二次震源波场正向延拓，并保存以所述激发时间为中心的指定时窗内的震源波场；
逆向延拓模块，用于基于所述反射波进行检波点波场逆向延拓，并保存以所述激发时间为中心的指定时窗内的检波点波场；
互相关成像模块，用于对指定时窗内的震源波场及检波点波场进行互相关成像。
本申请实施例的远探测声波逆时偏移成像装置，所述震源波场正向延拓基于以下时间正向延拓有限差分公式实现：
其中，n表示变量在时间网格的位置；lx和lz分别表示变量在空间网格x方向及z方向的位置；Δt表示时间间隔；Δx和Δz分别表示空间网格x方向及z方向的间隔；δx和δz分别表示空间x方向及z方向的差分；2N表示空间差分阶数；vx和vz分别表示震源波场速度变量的x分量及z分量；σxx和σzz分别表示震源波场应力变量的xx分量及zz分量；ρ表示能量密度；λ表示拉梅参数。
本申请实施例的远探测声波逆时偏移成像装置，所述在正向延拓过程中利用激发振幅成像条件获得激发时间，包括：
根据公式
确定激发振幅，并确定所述激发振幅对应的激发时间；
其中，E、K和W分别表示震源波场中的总能量密度、动能密度和势能密度；ρ表示能量密度；vi表示速度变量，σij表示应力变量，εij表示应变变量。
本申请实施例的远探测声波逆时偏移成像装置，所述检波点波场逆向延拓基于以下时间逆向延拓有限差分公式实现：
其中，n表示变量在时间网格的位置；lx和lz分别表示变量在空间网格x方向及z方向的位置；Δt表示时间间隔；Δx和Δz分别表示空间网格x方向及z方向的间隔；δx和δz分别表示空间x方向及z方向的差分；2N表示空间差分阶数；vx和vz分别表示震源波场速度变量的x分量及z分量；σxx和σzz分别表示震源波场应力变量的xx分量及zz分量；ρ表示能量密度；λ表示拉梅参数。
另一方面，本申请实施例还提供了另一种远探测声波逆时偏移成像装置，包括存储器、处理器、以及存储在所述存储器上的计算机程序，所述计算机程序被所述处理器运行时执行如下步骤：
获取远探测声波测井数据并进行预处理；
从预处理后的远探测声波测井数据中提取纵横波时差及反射波；
根据所述纵横波时差构建偏移速度模型；
基于所述偏移速度模型进行第一次震源波场正向延拓，以在正向延拓过程中利用激发振幅成像条件获得激发时间；
基于所述偏移速度模型进行第二次震源波场正向延拓，并保存以所述激发时间为中心的指定时窗内的震源波场；
基于所述反射波进行检波点波场逆向延拓，并保存以所述激发时间为中心的指定时窗内的检波点波场；
对指定时窗内的震源波场及检波点波场进行互相关成像。
由以上本申请实施例提供的技术方案可见，本申请实施例中，先基于偏移速度模型进行一次震源波场正向延拓，以在该正向延拓过程中利用激发振幅成像条件获得激发时间；然后在后续基于偏移速度模型再次震源波场正向延拓时，可以仅记录以激发时间为中心的指定时窗内的震源波场；并且，在后续基于反射波进行检波点波场逆向延拓时，可以仅记录以激发时间为中心的指定时窗内的检波点波场；在此基础上，可仅对指定时窗内的震源波场及检波点波场进行互相关成像，从而不仅可以大大节省存储空间，还大大减少了后续互相关成像的计算量，从而极大的提高了远探测声波逆时偏移成像的处理效率。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本申请一实施例中远探测声波逆时偏移成像方法的流程图；
图2a为常规逆时偏移成像技术的互相关示意图；
图2b为本申请一实施例中远探测声波逆时偏移成像方法的互相关示意图；
图3为本申请一实施例中单炮激发振幅的示意图；
图4为本申请一实施例中单炮激发时间的示意图；
图5a为采用常规逆时偏移成像技术获得的逆时偏移结果示意图；
图5b为采用本申请实施例的远探测声波逆时偏移成像方法获得的逆时偏移结果示意图；
图6为本申请一实施例中远探测声波逆时偏移成像装置的结构框图；
图7为本申请另一实施例中远探测声波逆时偏移成像装置的结构框图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。例如在下面描述中，在第一部件上方形成第二部件，可以包括第一部件和第二部件以直接接触方式形成的实施例，还可以包括第一部件和第二部件以非直接接触方式(即第一部件和第二部件之间还可以包括额外的部件)形成的实施例等。
而且，为了便于描述，本申请一些实施例可以使用诸如“在…上方”、“在…之下”、“顶部”、“下方”等空间相对术语，以描述如实施例各附图所示的一个元件或部件与另一个(或另一些)元件或部件之间的关系。应当理解的是，除了附图中描述的方位之外，空间相对术语还旨在包括装置在使用或操作中的不同方位。例如若附图中的装置被翻转，则被描述为“在”其他元件或部件“下方”或“之下”的元件或部件，随后将被定位为“在”其他元件或部件“上方”或“之上”。
参考图1所示，本申请实施例的远探测声波逆时偏移成像方法可以包括以下步骤：
S101、获取远探测声波测井数据并进行预处理。
本申请一实施例中，可利用远探测声波仪器进行测井以获取单极纵波或四分量偶极横波远探测数据，同时还可以记录仪器类型、采样间隔和采样点等信息，以便用于后期偏移处理的观测系统设计。在一示例性实施例中，某井例如可选择型号为XMAC-F1的远探测声波仪器进行远探测声波测井。其中，四分量偶极反射波资料时间采样间隔可以为36微秒，采样点数可以为400，深度采样间隔可以为0.1524米。
本申请一实施例中，所述的预处理例如可以是根据波形增益曲线和延迟曲线进行资料波形恢复，以利于后续处理。
S102、从预处理后的远探测声波测井数据中提取纵横波时差及反射波。
本申请一实施例中，可基于本领域常规纵横波时差技术从测井数据中提取纵横波时差，对此不再赘述。本申请一实施例中，所述从预处理后的远探测声波测井数据中提取反射波例如可以包括以下步骤：
首先，可根据中值滤波、FK滤波或多尺度相关对比法从预处理后的远探测声波测井数据中提取反射波。
然后，可根据倾斜中值滤波滤除所述反射波中的地层无效反射波。
S103、根据所述纵横波时差构建偏移速度模型。
本申请一实施例中，可基于本领域常规偏移速度建模技术构建偏移速度模型。
S104、基于所述偏移速度模型进行第一次震源波场正向延拓，以在正向延拓过程中利用激发振幅成像条件获得激发时间。
本申请一实施例中，可以利用有限差分法进行震源波场正向延拓。例如在一示例性实施例中，可基于以下时间正向延拓有限差分公式实现震源波场正向延拓：
其中，n表示变量在时间网格的位置；lx和lz分别表示变量在空间网格x方向及z方向的位置；Δt表示时间间隔；Δx和Δz分别表示空间网格x方向及z方向的间隔；δx和δz分别表示空间x方向及z方向的差分；2N表示空间差分阶数；vx和vz分别表示震源波场速度变量的x分量及z分量；σxx和σzz分别表示震源波场应力变量的xx分量及zz分量；ρ表示能量密度；λ表示拉梅参数。
本申请一实施例中，在震源波场正向延拓过程中，可计算所有网格上的最大能量密度(即激发振幅)(例如图3所示)，并确定最大能量出现的时间，即激发时间(例如图4所示)。在一示例性实施例中，可根据公式
确定激发振幅，并确定所述激发振幅对应的激发时间。其中，E、K和W分别表示震源波场中的总能量密度、动能密度和势能密度；ρ表示能量密度；vi表示速度变量，σij表示应力变量，εij表示应变变量。
S105、基于所述偏移速度模型进行第二次震源波场正向延拓，并保存以所述激发时间为中心的指定时窗内的震源波场。
本步骤中的震源波场正向延拓可参考步骤S104。在本申请一实施例中，可以步骤S104中的激发时间为中心来指定时窗，在此基础上，当基于所述偏移速度模型进行第二次震源波场正向延拓时，可以仅记录指定时窗内的震源波场。而常规逆时偏移成像技术需存储所有时间点的正传波场。因此，与常规逆时偏移成像技术相比，本申请实施例仅记录指定时窗内的震源波场，可以大大节省存储空间，且大大减少了后续互相关成像的数据量，从而极大的提高了远探测声波逆时偏移成像的处理效率。
S106、基于所述反射波进行检波点波场逆向延拓，并保存以所述激发时间为中心的指定时窗内的检波点波场。
本申请一实施例中，可以利用有限差分法进行检波点波场正向延拓。例如在一示例性实施例中，可基于以下时间逆向延拓有限差分公式实现检波点波场逆向延拓：
其中，n表示变量在时间网格的位置；lx和lz分别表示变量在空间网格x方向及z方向的位置；Δt表示时间间隔；Δx和Δz分别表示空间网格x方向及z方向的间隔；δx和δz分别表示空间x方向及z方向的差分；2N表示空间差分阶数；vx和vz分别表示震源波场速度变量的x分量及z分量；σxx和σzz分别表示震源波场应力变量的xx分量及zz分量；ρ表示能量密度；λ表示拉梅参数。
同样，在本申请一实施例中，可以步骤S104中的激发时间为中心来指定时窗，在此基础上，当基于所述反射波进行检波点波场逆向延拓时，可以仅记录指定时窗内的检波点波场。而常规逆时偏移成像技术需存储所有时间点的检波点波场。因此，与常规逆时偏移成像技术相比，本申请实施例仅记录指定时窗内的检波点波场，可以大大节省存储空间，且大大减少了后续互相关成像的数据量，从而极大的提高了远探测声波逆时偏移成像的处理效率。
S107、对指定时窗内的震源波场及检波点波场进行互相关成像。
本申请一实施例中，可以基于常规的互相关成像技术对指定时窗内的震源波场及检波点波场进行互相关成像，从而获得成像剖面。如图2a所示，为常规逆时偏移成像技术的互相关示意图；从图2a中可以看出，所有时间点的震源波场及所有时间点的检波点波场都参与了互相关成像。然而，实际上，只有部分时间点的震源波场及检波点波场是有效数据。因此，在本申请实施例中，可仅让指定时窗内的震源波场及检波点波场参与了互相关成像(例如图2b所示)，从而极大的降低了计算量，提高了远探测声波逆时偏移成像的处理效率。
在一示例性实施例中，图5a和图5b分别示出了采用常规逆时偏移成像技术获得的逆时偏移结果示意图；以及采用本申请实施例的远探测声波逆时偏移成像方法获得的逆时偏移结果示意图。对比图5a和图5b可以看出，采用本申请实施例的远探测声波逆时偏移成像方法获得的逆时偏移结果，基本可以达到采用常规逆时偏移成像技术获得的逆时偏移结果的效果。但是，相比于常规逆时偏移成像技术，本申请实施例的远探测声波逆时偏移成像方法的处理效率要高得多。
参考图6所示，本申请实施例的一种远探测声波逆时偏移成像装置可以包括：
数据获取模块61，可以用于获取远探测声波测井数据并进行预处理；
数据提取模块62，可以用于从预处理后的远探测声波测井数据中提取纵横波时差及反射波；
模型构建模块63，可以用于根据所述纵横波时差构建偏移速度模型；
第一正向延拓模块64，可以用于基于所述偏移速度模型进行第一次震源波场正向延拓，以在正向延拓过程中利用激发振幅成像条件获得激发时间；
第二正向延拓模块65，可以用于基于所述偏移速度模型进行第二次震源波场正向延拓，并保存以所述激发时间为中心的指定时窗内的震源波场；
逆向延拓模块66，可以用于基于所述反射波进行检波点波场逆向延拓，并保存以所述激发时间为中心的指定时窗内的检波点波场；
互相关成像模块67，可以用于对指定时窗内的震源波场及检波点波场进行互相关成像。
参考图7所示，本申请实施例的另一种远探测声波逆时偏移成像装置可以包括存储器、处理器、以及存储在所述存储器上的计算机程序，所述计算机程序被所述处理器运行时执行如下步骤：
获取远探测声波测井数据并进行预处理；
从预处理后的远探测声波测井数据中提取纵横波时差及反射波；
根据所述纵横波时差构建偏移速度模型；
基于所述偏移速度模型进行第一次震源波场正向延拓，以在正向延拓过程中利用激发振幅成像条件获得激发时间；
基于所述偏移速度模型进行第二次震源波场正向延拓，并保存以所述激发时间为中心的指定时窗内的震源波场；
基于所述反射波进行检波点波场逆向延拓，并保存以所述激发时间为中心的指定时窗内的检波点波场；
对指定时窗内的震源波场及检波点波场进行互相关成像。
虽然上文描述的过程流程包括以特定顺序出现的多个操作，但是，应当清楚了解，这些过程可以包括更多或更少的操作，这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
Claims (10)
1.一种远探测声波逆时偏移成像方法，其特征在于，包括：
获取远探测声波测井数据并进行预处理；
从预处理后的远探测声波测井数据中提取纵横波时差及反射波；
根据所述纵横波时差构建偏移速度模型；
基于所述偏移速度模型进行第一次震源波场正向延拓，以在正向延拓过程中利用激发振幅成像条件获得激发时间；
基于所述偏移速度模型进行第二次震源波场正向延拓，并保存以所述激发时间为中心的指定时窗内的震源波场；
基于所述反射波进行检波点波场逆向延拓，并保存以所述激发时间为中心的指定时窗内的检波点波场；
对指定时窗内的震源波场及检波点波场进行互相关成像。
2.如权利要求1所述的远探测声波逆时偏移成像方法，其特征在于，所述震源波场正向延拓基于以下时间正向延拓有限差分公式实现：
其中，n表示变量在时间网格的位置；lx和lz分别表示变量在空间网格x方向及z方向的位置；Δt表示时间间隔；δx和δz分别表示空间x方向及z方向的差分；vx和vz分别表示震源波场速度变量的x分量及z分量；σxx和σzz分别表示震源波场应力变量的xx分量及zz分量；ρ表示能量密度；λ表示拉梅参数。
3.如权利要求1所述的远探测声波逆时偏移成像方法，其特征在于，所述在正向延拓过程中利用激发振幅成像条件获得激发时间，包括：
根据公式
确定激发振幅，并确定所述激发振幅对应的激发时间；
其中，E、K和W分别表示震源波场中的总能量密度、动能密度和势能密度；ρ表示能量密度；vi表示速度变量，σij表示应力变量，εij表示应变变量。
4.如权利要求1所述的远探测声波逆时偏移成像方法，其特征在于，所述检波点波场逆向延拓基于以下时间逆向延拓有限差分公式实现：
其中，n表示变量在时间网格的位置；lx和lz分别表示变量在空间网格x方向及z方向的位置；Δt表示时间间隔；δx和δz分别表示空间x方向及z方向的差分；vx和vz分别表示震源波场速度变量的x分量及z分量；σxx和σzz分别表示震源波场应力变量的xx分量及zz分量；ρ表示能量密度；λ表示拉梅参数。
5.如权利要求1所述的远探测声波逆时偏移成像方法，其特征在于，所述从预处理后的远探测声波测井数据中提取反射波，包括：
根据中值滤波、FK滤波或多尺度相关对比法从预处理后的远探测声波测井数据中提取反射波；
根据倾斜中值滤波滤除所述反射波中的地层无效反射波。
6.一种远探测声波逆时偏移成像装置，其特征在于，包括：
数据获取模块，用于获取远探测声波测井数据并进行预处理；
数据提取模块，用于从预处理后的远探测声波测井数据中提取纵横波时差及反射波；
模型构建模块，用于根据所述纵横波时差构建偏移速度模型；
第一正向延拓模块，用于基于所述偏移速度模型进行第一次震源波场正向延拓，以在正向延拓过程中利用激发振幅成像条件获得激发时间；
第二正向延拓模块，用于基于所述偏移速度模型进行第二次震源波场正向延拓，并保存以所述激发时间为中心的指定时窗内的震源波场；
逆向延拓模块，用于基于所述反射波进行检波点波场逆向延拓，并保存以所述激发时间为中心的指定时窗内的检波点波场；
互相关成像模块，用于对指定时窗内的震源波场及检波点波场进行互相关成像。
7.如权利要求6所述的远探测声波逆时偏移成像装置，其特征在于，所述震源波场正向延拓基于以下时间正向延拓有限差分公式实现：
其中，n表示变量在时间网格的位置；lx和lz分别表示变量在空间网格x方向及z方向的位置；Δt表示时间间隔；δx和δz分别表示空间x方向及z方向的差分；vx和vz分别表示震源波场速度变量的x分量及z分量；σxx和σzz分别表示震源波场应力变量的xx分量及zz分量；ρ表示能量密度；λ表示拉梅参数。
8.如权利要求6所述的远探测声波逆时偏移成像装置，其特征在于，所述在正向延拓过程中利用激发振幅成像条件获得激发时间，包括：
根据公式
确定激发振幅，并确定所述激发振幅对应的激发时间；
其中，E、K和W分别表示震源波场中的总能量密度、动能密度和势能密度；ρ表示能量密度；vi表示速度变量，σij表示应力变量，εij表示应变变量。
9.如权利要求6所述的远探测声波逆时偏移成像装置，其特征在于，所述检波点波场逆向延拓基于以下时间逆向延拓有限差分公式实现：
其中，n表示变量在时间网格的位置；lx和lz分别表示变量在空间网格x方向及z方向的位置；Δt表示时间间隔；δx和δz分别表示空间x方向及z方向的差分；vx和vz分别表示震源波场速度变量的x分量及z分量；σxx和σzz分别表示震源波场应力变量的xx分量及zz分量；ρ表示能量密度；λ表示拉梅参数。
10.一种远探测声波逆时偏移成像装置，包括存储器、处理器、以及存储在所述存储器上的计算机程序，其特征在于，所述计算机程序被所述处理器运行时执行如下步骤：
获取远探测声波测井数据并进行预处理；
从预处理后的远探测声波测井数据中提取纵横波时差及反射波；
根据所述纵横波时差构建偏移速度模型；
基于所述偏移速度模型进行第一次震源波场正向延拓，以在正向延拓过程中利用激发振幅成像条件获得激发时间；
基于所述偏移速度模型进行第二次震源波场正向延拓，并保存以所述激发时间为中心的指定时窗内的震源波场；
基于所述反射波进行检波点波场逆向延拓，并保存以所述激发时间为中心的指定时窗内的检波点波场；
对指定时窗内的震源波场及检波点波场进行互相关成像。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810966059.4A CN108897046B (zh) | 2018-08-23 | 2018-08-23 | 远探测声波逆时偏移成像方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810966059.4A CN108897046B (zh) | 2018-08-23 | 2018-08-23 | 远探测声波逆时偏移成像方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108897046A CN108897046A (zh) | 2018-11-27 |
| CN108897046B true CN108897046B (zh) | 2020-05-08 |
Family
ID=64358286
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810966059.4A Active CN108897046B (zh) | 2018-08-23 | 2018-08-23 | 远探测声波逆时偏移成像方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108897046B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110018517B (zh) * | 2019-05-07 | 2021-04-02 | 西安石油大学 | 一种多尺度地面微地震逆时干涉定位方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6956790B2 (en) * | 2003-04-25 | 2005-10-18 | Schlumberger Technology Corporation | Borehole sonic data processing method |
| CN103758511B (zh) * | 2013-11-25 | 2016-08-17 | 中国石油天然气股份有限公司 | 一种井下逆时偏移成像识别隐蔽储层的方法及装置 |
| CN104533396A (zh) * | 2014-12-31 | 2015-04-22 | 中国石油天然气集团公司 | 一种远探测声波的处理方法 |
| CN106199721B (zh) * | 2016-07-04 | 2018-06-15 | 中国石油集团川庆钻探工程有限公司 | 从阵列声波测井资料中提取反射波的方法 |
| CN106526678B (zh) * | 2016-10-26 | 2019-08-23 | 中国科学院声学研究所 | 一种反射声波测井的波场分离方法及装置 |
| CN108035708B (zh) * | 2017-11-20 | 2021-04-30 | 中国石油天然气股份有限公司 | 一种去除地层界面反射波的方法及装置 |
-
2018
- 2018-08-23 CN CN201810966059.4A patent/CN108897046B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108897046A (zh) | 2018-11-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Claerbout | Synthesis of a layered medium from its acoustic transmission response | |
| CN110133715B (zh) | 一种基于初至时差和波形叠加的微地震震源定位方法 | |
| Sava | Micro-earthquake monitoring with sparsely sampled data | |
| CN102176054B (zh) | 近地表综合信息处理解释方法 | |
| CN103926623B (zh) | 一种压制逆时偏移低频噪音的方法 | |
| CN103257361A (zh) | 基于Zoeppritz方程近似式的油气预测方法及系统 | |
| CN102053263B (zh) | 调查表层结构的方法 | |
| CN103758511B (zh) | 一种井下逆时偏移成像识别隐蔽储层的方法及装置 | |
| US8867308B2 (en) | System, method and computer program product for stacking seismic noise data to analyze seismic events | |
| CN103792576A (zh) | 一种基于梯度结构张量的储层非均质检测方法及设备 | |
| RU2737846C2 (ru) | Система установки наземных сейсмических датчиков с парами смежных многокомпонентных сейсмических датчиков на расстоянии в среднем по меньшей мере двадцать метров | |
| US9952341B2 (en) | Systems and methods for aligning a monitor seismic survey with a baseline seismic survey | |
| CN104459777A (zh) | 基于流体体积模量avo反演的流体识别的方法及系统 | |
| CN109001813A (zh) | 一种压制多次波的方法、装置及系统 | |
| CN112305591B (zh) | 隧道超前地质预报方法、计算机可读存储介质 | |
| CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
| CN117546051A (zh) | 使用s波速度模型和机器学习进行地震成像的方法和系统 | |
| CN114924315A (zh) | 一种多态式地震勘探方法和系统 | |
| Pafeng et al. | Prestack waveform inversion of three-dimensional seismic data—An example from the Rock Springs Uplift, Wyoming, USA | |
| CN109709602B (zh) | 一种远探测声波偏移成像方法、装置及系统 | |
| CN108897046B (zh) | 远探测声波逆时偏移成像方法及装置 | |
| CN104199088A (zh) | 一种提取入射角道集的方法及系统 | |
| Li et al. | Relative elastic interferometric imaging for microseismic source location | |
| CN103543467B (zh) | 地震弱信号增强处理方法 | |
| RU2011148308A (ru) | Способ комплексной обработки геофизических данных и технологическая система "литоскан" для его осуществления |
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 |