CN111948707A - 基于小层波形校正的多波地震资料匹配方法及装置 - Google Patents
基于小层波形校正的多波地震资料匹配方法及装置 Download PDFInfo
- Publication number
- CN111948707A CN111948707A CN201910408871.XA CN201910408871A CN111948707A CN 111948707 A CN111948707 A CN 111948707A CN 201910408871 A CN201910408871 A CN 201910408871A CN 111948707 A CN111948707 A CN 111948707A
- Authority
- CN
- China
- Prior art keywords
- wave
- data
- seismic
- waveform
- wave data
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 66
- 238000012937 correction Methods 0.000 title claims abstract description 55
- 238000012545 processing Methods 0.000 claims abstract description 56
- 238000006073 displacement reaction Methods 0.000 claims abstract description 39
- 238000009499 grossing Methods 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 15
- 230000006870 function Effects 0.000 claims description 13
- 238000005070 sampling Methods 0.000 claims description 13
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 238000003491 array Methods 0.000 claims 1
- 239000003208 petroleum Substances 0.000 description 11
- 238000004891 communication Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 10
- 238000007405 data analysis Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012546 transfer Methods 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/28—Processing seismic data, e.g. for interpretation or for event detection
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
本申请提供一种基于小层波形校正的多波地震资料匹配方法及装置，方法包括：针对同一目标地质区域的多波地震资料中的PP波数据和PS波数据进行初始匹配；若相对于PP波数据对应的地震剖面，经初始匹配后的PS波数据对应的地震剖面出现局部畸变，则以PP波数据和PS波数据的纵波道数据为基准，拉平PP波数据和PS波数据的横波同一位置的波形数据；依次获取被拉平的波形数据的位移量，并根据各个位移量对纵波和横波之间进行小层拉平处理，以消除PS波数据对应的地震剖面的局部畸变。本申请能够有效消除多波地震资料的初始匹配后出现的PS波数据的局部畸变现象，实现多波地震资料的小层多波数据精细匹配。
Description
技术领域
本申请涉及石油勘探技术领域，具体涉及一种基于小层波形校正的多波地震资料匹配方法及装置。
背景技术
石油与天然气地震勘探中，除了应用常规的地震反射纵波(以下称为PP波)数据来进行地质构造和储层研究外，还需要结合转换横波(以下称为PS波)地震数据来进行综合研究，特别是在气云成像、裂缝预测、岩性和流体识别等领域。这种应用PP波与PS波地震数据，进行综合地质分析的技术，称为多波地震解释。由于PP波与PS波在地下传播速度有差异，造成PP波与PS波两套地震数据在时间上产生了明显的差异(时差)，只有消除这一差异，才能使这两类数据在同一尺度上进行研究分析，最常用的消除差异方法是以PP波地震数据为基准，将PS波地震数据压缩到与PP波地震数据一致的时间范围内，这项技术被称为多波数据匹配，参见图1和图2。
多波匹配技术主要包括解释层位匹配(层位是分别在PP波地震剖面上和PS波地震剖面上各自解释的层位)和逐点交互匹配(分别在PP波地震剖面上和PS波地震剖面上确立属于一个地质反射界面的对比点)，用这些方法对PP波、PS波数据进行初始匹配，初始匹配后PP波与PS波数据大套的地质反射界面能够基本对应一致，但由于进行数据匹配的控制点，在纵向上和横向上都是离散分布的，这就造成匹配后的PS波数据，某些区域会不同程度地发生地震同相轴(地震记录上各道振动相位相同的波峰或波谷的连线)的异常扭曲和错断现象，参见图3，因此需要进一步更为准确的多波数据匹配手段。
发明内容
针对现有技术中的问题，本申请提供一种基于小层波形校正的多波地震资料匹配方法及装置，能够有效消除多波地震资料的初始匹配后出现的PS波数据的局部畸变现象，实现所述多波地震资料的小层多波数据精细匹配。
为解决上述技术问题，本申请提供以下技术方案：
第一方面，本申请提供一种基于小层波形校正的多波地震资料匹配方法，包括：
针对同一目标地质区域的多波地震资料中的PP波数据和PS波数据进行初始匹配；
若相对于所述PP波数据对应的地震剖面，经初始匹配后的所述PS波数据对应的地震剖面出现局部畸变，则以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据；
依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理，以消除所述PS波数据对应的地震剖面的所述局部畸变。
进一步地，所述以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据，包括：
分别在针对同一目标地质区域的多波地震资料中抽取一道PP波和一道PS波；
应用二维数组表示抽取的所述PP波和PS波；
根据预先获取的同相轴精细拉平处理参数中的波形对比窗口的大小及波形对比窗口的移动量，依次判断是否在某个所述波形对比窗口的中心点设置拉平种子点，若是，则通过设置该拉平种子点拉平所述PP波数据和PS波数据的横波同一位置的波形数据。
进一步地，所述应用二维数组表示抽取的所述PP波和PS波，包括：
将抽取的一道所述PP波和一道所述PS波存储成一个用于进行同相轴精细拉平处理的道集，并应用I行及J列的二维数组D来表示该道集；
其中，I用于表示每个地震道的采样点数，J用于表示道集的地震道数，且J＝2，在所述二维数组D中，第i行及第j列的值D(i，j)用于表示该道集的第j个地震道的第i个采样点的振幅值。
进一步地，所述同相轴精细拉平处理参数中还包含有所述波形对比窗口的搜索半径。
进一步地，所述依次获取被拉平的所述波形数据的位移量，包括：
依次获取所述波形对比窗口的一个面上的所述拉平种子点的移动量；
依次应用平滑函数对所述波形对比窗口的一个面上的所述拉平种子点的移动量进行横向平滑处理。
进一步地，所述根据各个所述位移量对所述纵波和横波之间进行小层拉平处理，包括：
根据各个所述拉平种子点的移动量建立拉平后坐标，并将原始的PS数据根据该拉平后坐标进行差值，而后返回重新抽取一道PP波和一道PS波后，重复执行所述以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据，以及依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理的迭代过程，直至消除所述PS波数据对应的地震剖面的全部的所述局部畸变。
第二方面，本申请提供一种基于小层波形校正的多波地震资料匹配装置，包括：
初始匹配模块，用于针对同一目标地质区域的多波地震资料中的PP波数据和PS波数据进行初始匹配；
波形数据拉平模块，用于若相对于所述PP波数据对应的地震剖面，经初始匹配后的所述PS波数据对应的地震剖面出现局部畸变，则以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据；
局部畸变消除模块，用于依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理，以消除所述PS波数据对应的地震剖面的所述局部畸变。
进一步地，所述波形数据拉平模块包括：
抽取单元，用于分别在针对同一目标地质区域的多波地震资料中抽取一道PP波和一道PS波；
二维数组表示单元，用于应用二维数组表示抽取的所述PP波和PS波；
拉平种子点设置单元，用于根据预先获取的同相轴精细拉平处理参数中的波形对比窗口的大小及波形对比窗口的移动量，依次判断是否在某个所述波形对比窗口的中心点设置拉平种子点，若是，则通过设置该拉平种子点拉平所述PP波数据和PS波数据的横波同一位置的波形数据。
进一步地，所述二维数组表示单元包括：
道集表示子单元，用于将抽取的一道所述PP波和一道所述PS波存储成一个用于进行同相轴精细拉平处理的道集，并应用I行及J列的二维数组D来表示该道集；
其中，I用于表示每个地震道的采样点数，J用于表示道集的地震道数，且J＝2，在所述二维数组D中，第i行及第j列的值D(i，j)用于表示该道集的第j个地震道的第i个采样点的振幅值。
进一步地，所述同相轴精细拉平处理参数中还包含有所述波形对比窗口的搜索半径。
进一步地，所述局部畸变消除模块包括：
移动量获取单元，用于依次获取所述波形对比窗口的一个面上的所述拉平种子点的移动量；
横向平滑处理单元，用于依次应用平滑函数对所述波形对比窗口的一个面上的所述拉平种子点的移动量进行横向平滑处理。
进一步地，所述局部畸变消除模块还包括：
迭代处理单元，用于根据各个所述拉平种子点的移动量建立拉平后坐标，并将原始的PS数据根据该拉平后坐标进行差值，而后返回重新抽取一道PP波和一道PS波后，重复执行所述以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据，以及依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理的迭代过程，直至消除所述PS波数据对应的地震剖面的全部的所述局部畸变。
第三方面，本申请提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述基于小层波形校正的多波地震资料匹配方法的步骤。
第四方面，本申请提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现所述的基于小层波形校正的多波地震资料匹配方法的步骤。
由上述技术方案可知，本申请提供一种基于小层波形校正的多波地震资料匹配方法及装置，通过针对同一目标地质区域的多波地震资料中的PP波数据和PS波数据进行初始匹配；若相对于所述PP波数据对应的地震剖面，经初始匹配后的所述PS波数据对应的地震剖面出现局部畸变，则以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据；依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理，以消除所述PS波数据对应的地震剖面的所述局部畸变，能够有效消除多波地震资料的初始匹配后出现的PS波数据的局部畸变现象，实现所述多波地震资料的小层多波数据精细匹配，进而能够根据该精细匹配结果提高石油数据分析的准确性，以进一步提高石油勘探的准确性和可靠性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为初始匹配前的PP-PS多波数据比较示意图。
图2为初始匹配后的PP-PS多波数据匹配示意图。
图3为初始匹配后PS波地震剖面示意图。
图4为本申请实施例中的基于小层波形校正的多波地震资料匹配系统的结构示意图。
图5为本申请实施例中的基于小层波形校正的多波地震资料匹配方法的流程示意图。
图6为本申请应用实例中的基于小层波形校正的多波地震资料匹配方法的结构示意图。
图7为应用所述基于小层波形校正的多波地震资料匹配方法之前与之后的PS波地震剖面的比较示意图。
图8为本申请实施例中的基于小层波形校正的多波地震资料匹配装置的结构示意图。
图9为本申请实施例中的电子设备的结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
考虑到现有的多波地震资料匹配的方式中存在的匹配结果不准确，无法实现对多波地震资料的精细匹配的问题，本申请提供一种基于小层波形校正的多波地震资料匹配方法及装置，通过针对同一目标地质区域的多波地震资料中的PP波数据和PS波数据进行初始匹配，得到所述PS波数据对应的地震剖面和所述PP波数据对应的地震剖面；若相对于所述PP波数据对应的地震剖面，所述PS波数据对应的地震剖面出现局部畸变，则以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据；依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理，以消除所述PS波数据对应的地震剖面的所述局部畸变，能够有效消除多波地震资料的初始匹配后出现的PS波数据的局部畸变现象，实现所述多波地震资料的小层多波数据精细匹配，进而能够根据该精细匹配结果提高石油数据分析的准确性，以进一步提高石油勘探的准确性和可靠性。
基于上述内容，本申请还提供一种基于小层波形校正的多波地震资料匹配系统，基于小层波形校正的多波地震资料匹配系统包含有一基于小层波形校正的多波地震资料匹配装置，参见图4，该基于小层波形校正的多波地震资料匹配装置可以为一种服务器1，该服务器1可以与至少一个客户端设备2通信连接，服务器1还可以与石油勘探现场设备通信连接。其中，用户可以通过登录客户端设备2的网站或者APP等，与服务器1建立连接。服务器1可以根据用户在线自客户端设备2发出的操作指令在线或离线执行基于小层波形校正的多波地震资料匹配过程，具体可以针对同一目标地质区域的多波地震资料中的PP波数据和PS波数据进行初始匹配，得到所述PS波数据对应的地震剖面和所述PP波数据对应的地震剖面；若相对于所述PP波数据对应的地震剖面，所述PS波数据对应的地震剖面出现局部畸变，则以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据；依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理，以消除所述PS波数据对应的地震剖面的所述局部畸变，并将小层拉平处理结果进行显示，或将小层拉平处理结果在线发送至客户端设备2，使得客户端设备2对小层拉平处理结果进行显示。
可以理解的是，所述客户端设备2可以包括智能手机、平板电子设备、网络机顶盒、便携式计算机、台式电脑、个人数字助理(PDA)、车载设备、智能穿戴设备等。其中，所述智能穿戴设备可以包括智能眼镜、智能手表、智能手环等。
在实际应用中，进行基于小层波形校正的多波地震资料匹配的部分可以在如上述内容所述的服务器1侧执行，即，如图4所示的架构，也可以所有的操作都在所述客户端设备2中完成。具体可以根据所述客户端设备2的处理能力，以及用户使用场景的限制等进行选择。本申请对此不作限定。若所有的操作都在所述客户端设备2中完成，所述客户端设备2还可以包括处理器。
上述的客户端设备可以具有通信模块(即通信单元)，可以与远程的服务器进行通信连接，实现与所述服务器的数据传输。所述服务器可以包括任务调度中心一侧的服务器，其他的实施场景中也可以包括中间平台的服务器，例如与任务调度中心服务器有通信链接的第三方服务器平台的服务器。所述的服务器可以包括单台计算机设备，也可以包括多个服务器组成的服务器集群，或者分布式装置的服务器结构。
所述服务器与所述客户端设备之间可以使用任何合适的网络协议进行通信，包括在本申请提交日尚未开发出的网络协议。所述网络协议例如可以包括TCP/IP协议、UDP/IP协议、HTTP协议、HTTPS协议等。当然，所述网络协议例如还可以包括在上述协议之上使用的RPC协议(Remote Procedure Call Protocol，远程过程调用协议)、REST协议(Representational State Transfer，表述性状态转移协议)等。
为了能够有效消除多波地震资料的初始匹配后出现的PS波数据的局部畸变现象，实现所述多波地震资料的小层多波数据精细匹配，本申请实施例提供一种执行主体为基于小层波形校正的多波地震资料匹配装置的基于小层波形校正的多波地震资料匹配方法，参见图5，所述基于小层波形校正的多波地震资料匹配方法具体包含有如下内容：
步骤100：针对同一目标地质区域的多波地震资料中的PP波数据和PS波数据进行初始匹配。
可以理解的是，如果把测量点作为地球的一极，穿过地心和测量点相对应的为另一极，中间的距离按照和地心与测量点之间的张角划分为180度，也可以说震中距在10度以内的称为近震，10度——105度的为远震，105度以上的称为极远震。对于远震和极远震，不能再象近震一样认为地震波是直线传播的，必须考虑地震波在传播时的折射。从震源发出，以纵波传播到中途，在地表反射一次后仍以纵波形式传播到测量点的称为PP波。以及，从震源发出，以纵波传播到中途，在地表反射后以横波形式传播到测量点的称为PS波。
步骤200：若相对于所述PP波数据对应的地震剖面，经初始匹配后的所述PS波数据对应的地震剖面出现局部畸变，则以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据。
步骤300：依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理，以消除所述PS波数据对应的地震剖面的所述局部畸变。
在本申请的基于小层波形校正的多波地震资料匹配方法的一个实施例中，所述步骤200具体包含有如下内容：
步骤201：分别在针对同一目标地质区域的多波地震资料中抽取一道PP波和一道PS波。
步骤202：应用二维数组表示抽取的所述PP波和PS波。
步骤203：根据预先获取的同相轴精细拉平处理参数中的波形对比窗口的大小及波形对比窗口的移动量，依次判断是否在某个所述波形对比窗口的中心点设置拉平种子点，若是，则通过设置该拉平种子点拉平所述PP波数据和PS波数据的横波同一位置的波形数据。
在本申请的基于小层波形校正的多波地震资料匹配方法的一个实施例中，所述步骤202具体包含有如下内容：
步骤202A：将抽取的一道所述PP波和一道所述PS波存储成一个用于进行同相轴精细拉平处理的道集，并应用I行及J列的二维数组D来表示该道集；其中，I用于表示每个地震道的采样点数，J用于表示道集的地震道数，且J＝2，在所述二维数组D中，第i行及第j列的值D(i，j)用于表示该道集的第j个地震道的第i个采样点的振幅值。
可以理解的是，所述同相轴精细拉平处理参数中还包含有所述波形对比窗口的搜索半径。
在本申请的基于小层波形校正的多波地震资料匹配方法的一个实施例中，所述步骤300具体包含有如下内容：
步骤301：依次获取所述波形对比窗口的一个面上的所述拉平种子点的移动量；
步骤302：依次应用平滑函数对所述波形对比窗口的一个面上的所述拉平种子点的移动量进行横向平滑处理。
步骤302：根据各个所述拉平种子点的移动量建立拉平后坐标，并将原始的PS数据根据该拉平后坐标进行差值，而后返回重新抽取一道PP波和一道PS波后，重复执行所述以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据，以及依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理的迭代过程，直至消除所述PS波数据对应的地震剖面的全部的所述局部畸变。
从上述描述可知，本申请实施例提供的基于小层波形校正的多波地震资料匹配方法，能够有效消除多波地震资料的初始匹配后出现的PS波数据的局部畸变现象，实现所述多波地震资料的小层多波数据精细匹配，进而能够根据该精细匹配结果提高石油数据分析的准确性，以进一步提高石油勘探的准确性和可靠性。
为进一步说明本方案，在纵、横波地震数据初始匹配的基础上，进行小层多波数据的精细匹配，本申请还提供一种基于小层波形校正的多波地震资料匹配方法的具体实施方式，以纵波道数据为基准，拉平经过初始匹配的横波同一位置的数据，计算拉平的位移量，根据位移量实现纵、横波之间的小层拉平，从而达到小层多波数据精细匹配的目的，以消除初始匹配后PS波数据出现的局部畸变现象，参见图6，所述基于小层波形校正的多波地震资料匹配方法具体包含有如下内容：
S1：分别抽取针对同一目标地质区域的地震资料的一道PP波和一道PS波，具体为抽取PP波数据的一道数据trace1和PS波对应的一道数据trace2。
S2：用I行2列的二维数组D表示，I表示一道数据的点数，2表示道数，具体为：
将两个数据道存储成一个将要进行同相轴精细拉平处理的道集，其中trace1为第一道，trace2为第二道。用I行、J列的二维数组D来表示。其中I为每个地震道的采样点数，J为道集的地震道数，这里取J为2。二维数组D，其第i行，第j列的值D(i,j)表示道集的第j个地震道的第i个采样点的振幅值。
S3：设定同相轴精细拉平处理参数，具体为：
设定同相轴精细拉平处理参数：Nw、Nm、Ns、Toli、Tolg，其中，Nw为波形对比窗口的大小，Nw为整数；Nm为波形对比窗口的移动量，Nm为整数；Ns为波形对比窗口的搜索半径，Ns为整数。
S4：以Nw为窗口大小，以Nm为窗口移动量，逐窗口判定是否在某个波形对比窗口的中心点设置拉平种子点，具体为：
以Nw为窗口，以Nw为窗口移动量，逐窗口判定是否在某个波形对比窗口的中心点设置拉平种子点，并计算拉平种子点移动量，具体计算方法为：对于第k个波形对比窗口，计算地震道与地震道之间的最优移动量及最大相关系数。定义J行J列的二维数组C用来存储地震道与地震道之间的最大相关系数。所述的二维数组C，其第j1行、第j2列的值C(j1,j2)表示道集的第j1个地震道和第j2个地震道的最大相关系数。所述的二维数组S，其第j1行、第j2列的值S(j1,j2)表示道集的第j1个地震道和第j2个地震道的最优移动量。其具体计算方法为：以计算C(j1,j2)和S(j1,j2)为例，设is＝(k-1)Nm，ie＝is+Nw，定义集合
定义长度为(Nw+1)的向量f，f的第i个元素
定义2Ns+1个长度为(Nw+1)的向量gl，gl的第i个元素
则最大相关系数C(j1,j2)由以下公式给出：
最优移动量S(j1,j2)，即取得C(j1,j2)时的l值，有以下公式给出：
对于第k个波形对比窗口,计算拉平种子点移动量,定义长度为J的向量m用来存储每个地震道的拉平种子点移动量,其中,所述的计算拉平种子点移动量为m(j2)＝S(j1,j2)-S(j1,j1)。
S5：计算拉平种子点移动量，具体为：
求出一个面上的种子点的移动量，利用平滑函数，对该面上种子点的移动量进行横向平滑处理。即对每一道同一坐标点上的种子点的移动量进行平滑处理。
S6：计算拉伸采样坐标，具体为：
在计算拉伸采样坐标中，定义I行，J列的二维数组Xnew来存储拉伸采样坐标，所述的计算拉伸采样坐标的方法为：假设需要在Q个对比窗口的中心点设置拉平种子点，设其中第q个中心点以及它所对应的拉平种子点移动量分别为iq和
定义集合
则对于任意一个定值
Xnew(i,j)＝i+mi(j)；
S7：根据拉伸采样坐标实施同相轴精细拉平，具体为：
根据Xnew里种子点的坐标线性插出每一道上每个点的采样坐标。
定义I行J列的二维数组X用来存储原始采样坐标，定义I行J列的二维数组Dnew用来存储同相轴精细拉平处理后的道集；
设置原始采样坐标，具体地，X(i,j)＝i,i＝1,2,...,I；j＝1,2,...,J；
对于任意一个定值
以X的第j列Xj为自变量、以D的第j列Dj为函数值，构建三次样条差值函数
则
用Dnew的第一列代替trace1，用Dnew的第二列代替trace2。
然后循环的拉平其它道，并且用拉平后的地震道数据替换初始的地震道数据。
将计算好的PS波再次代入算法开始S1，作为PS波进行输入，不断对效果进行迭代修正。
应用本应用实例的基于小层波形校正的多波地震资料匹配方法进行精细匹配后，PS波地震剖面消除不正常时差，同相轴横向连续性得到了明显改善，参见图7。
为了能够有效消除多波地震资料的初始匹配后出现的PS波数据的局部畸变现象，实现所述多波地震资料的小层多波数据精细匹配，本申请实施例提供一种用于实现所述基于小层波形校正的多波地震资料匹配方法中的全部或部分内容的基于小层波形校正的多波地震资料匹配装置，参见图8，所述基于小层波形校正的多波地震资料匹配装置具体包含有如下内容：
初始匹配模块10，用于针对同一目标地质区域的多波地震资料中的PP波数据和PS波数据进行初始匹配。
波形数据拉平模块20，用于若相对于所述PP波数据对应的地震剖面，经初始匹配后的所述PS波数据对应的地震剖面出现局部畸变，则以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据。
局部畸变消除模块30，用于依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理，以消除所述PS波数据对应的地震剖面的所述局部畸变。
在本申请提供的所述基于小层波形校正的多波地震资料匹配装置的实施例中，所述波形数据拉平模块20具体包含有如下内容：
抽取单元，用于分别在针对同一目标地质区域的多波地震资料中抽取一道PP波和一道PS波；
二维数组表示单元，用于应用二维数组表示抽取的所述PP波和PS波；
拉平种子点设置单元，用于根据预先获取的同相轴精细拉平处理参数中的波形对比窗口的大小及波形对比窗口的移动量，依次判断是否在某个所述波形对比窗口的中心点设置拉平种子点，若是，则通过设置该拉平种子点拉平所述PP波数据和PS波数据的横波同一位置的波形数据。
在本申请提供的所述基于小层波形校正的多波地震资料匹配装置的实施例中，所述二维数组表示单元具体包含有如下内容：
道集表示子单元，用于将抽取的一道所述PP波和一道所述PS波存储成一个用于进行同相轴精细拉平处理的道集，并应用I行及J列的二维数组D来表示该道集；
其中，I用于表示每个地震道的采样点数，J用于表示道集的地震道数，且J＝2，在所述二维数组D中，第i行及第j列的值D(i，j)用于表示该道集的第j个地震道的第i个采样点的振幅值。
其中，所述同相轴精细拉平处理参数中还包含有所述波形对比窗口的搜索半径。
在本申请提供的所述基于小层波形校正的多波地震资料匹配装置的实施例中，所述局部畸变消除模块30具体包含有如下内容：
移动量获取单元，用于依次获取所述波形对比窗口的一个面上的所述拉平种子点的移动量；
横向平滑处理单元，用于依次应用平滑函数对所述波形对比窗口的一个面上的所述拉平种子点的移动量进行横向平滑处理。
迭代处理单元，用于根据各个所述拉平种子点的移动量建立拉平后坐标，并将原始的PS数据根据该拉平后坐标进行差值，而后返回重新抽取一道PP波和一道PS波后，重复执行所述以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据，以及依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理的迭代过程，直至消除所述PS波数据对应的地震剖面的全部的所述局部畸变。
从上述描述可知，本申请实施例提供的基于小层波形校正的多波地震资料匹配装置，能够有效消除多波地震资料的初始匹配后出现的PS波数据的局部畸变现象，实现所述多波地震资料的小层多波数据精细匹配，进而能够根据该精细匹配结果提高石油数据分析的准确性，以进一步提高石油勘探的准确性和可靠性。
从硬件层面来看，本申请的实施例还提供能够实现上述实施例中的基于小层波形校正的多波地震资料匹配方法中全部步骤的一种电子设备的具体实施方式，参见图9，所述电子设备具体包括如下内容：
处理器(processor)601、存储器(memory)602、通信接口(CommunicationsInterface)603和总线604；
其中，所述处理器601、存储器602、通信接口603通过所述总线604完成相互间的通信；所述通信接口603用于实现基于小层波形校正的多波地震资料匹配装置、客户终端、现场设备以及其他参与机构之间的信息传输；
所述处理器601用于调用所述存储器602中的计算机程序，所述处理器执行所述计算机程序时实现上述实施例中的基于小层波形校正的多波地震资料匹配方法中的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：
步骤100：针对同一目标地质区域的多波地震资料中的PP波数据和PS波数据进行初始匹配。
步骤200：若相对于所述PP波数据对应的地震剖面，经初始匹配后的所述PS波数据对应的地震剖面出现局部畸变，则以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据。
步骤300：依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理，以消除所述PS波数据对应的地震剖面的所述局部畸变。
从上述描述可知，本申请实施例提供的电子设备，能够有效消除多波地震资料的初始匹配后出现的PS波数据的局部畸变现象，实现所述多波地震资料的小层多波数据精细匹配，进而能够根据该精细匹配结果提高石油数据分析的准确性，以进一步提高石油勘探的准确性和可靠性。
本申请的实施例还提供能够实现上述实施例中的基于小层波形校正的多波地震资料匹配方法中全部步骤的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的基于小层波形校正的多波地震资料匹配方法的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：
步骤100：针对同一目标地质区域的多波地震资料中的PP波数据和PS波数据进行初始匹配。
步骤200：若相对于所述PP波数据对应的地震剖面，经初始匹配后的所述PS波数据对应的地震剖面出现局部畸变，则以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据。
步骤300：依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理，以消除所述PS波数据对应的地震剖面的所述局部畸变。
从上述描述可知，本申请实施例提供的计算机可读存储介质，能够有效消除多波地震资料的初始匹配后出现的PS波数据的局部畸变现象，实现所述多波地震资料的小层多波数据精细匹配，进而能够根据该精细匹配结果提高石油数据分析的准确性，以进一步提高石油勘探的准确性和可靠性。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。
本说明书实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书的实施例而已，并不用于限制本说明书实施例。对于本领域技术人员来说，本说明书实施例可以有各种更改和变化。凡在本说明书实施例的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书实施例的权利要求范围之内。
Claims (14)
1.一种基于小层波形校正的多波地震资料匹配方法，其特征在于，包括：
针对同一目标地质区域的多波地震资料中的PP波数据和PS波数据进行初始匹配；
若相对于所述PP波数据对应的地震剖面，经初始匹配后的所述PS波数据对应的地震剖面出现局部畸变，则以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据；
依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理，以消除所述PS波数据对应的地震剖面的所述局部畸变。
2.根据权利要求1所述的基于小层波形校正的多波地震资料匹配方法，其特征在于，所述以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据，包括：
分别在针对同一目标地质区域的多波地震资料中抽取一道PP波和一道PS波；
应用二维数组表示抽取的所述PP波和PS波；
根据预先获取的同相轴精细拉平处理参数中的波形对比窗口的大小及波形对比窗口的移动量，依次判断是否在某个所述波形对比窗口的中心点设置拉平种子点，若是，则通过设置该拉平种子点拉平所述PP波数据和PS波数据的横波同一位置的波形数据。
3.根据权利要求2所述的基于小层波形校正的多波地震资料匹配方法，其特征在于，所述应用二维数组表示抽取的所述PP波和PS波，包括：
将抽取的一道所述PP波和一道所述PS波存储成一个用于进行同相轴精细拉平处理的道集，并应用I行及J列的二维数组D来表示该道集；
其中，I用于表示每个地震道的采样点数，J用于表示道集的地震道数，且J＝2，在所述二维数组D中，第i行及第j列的值D(i，j)用于表示该道集的第j个地震道的第i个采样点的振幅值。
4.根据权利要求2所述的基于小层波形校正的多波地震资料匹配方法，其特征在于，所述同相轴精细拉平处理参数中还包含有所述波形对比窗口的搜索半径。
5.根据权利要求2至4任一项所述的基于小层波形校正的多波地震资料匹配方法，其特征在于，所述依次获取被拉平的所述波形数据的位移量，包括：
依次获取所述波形对比窗口的一个面上的所述拉平种子点的移动量；
依次应用平滑函数对所述波形对比窗口的一个面上的所述拉平种子点的移动量进行横向平滑处理。
6.根据权利要求5所述的基于小层波形校正的多波地震资料匹配方法，其特征在于，所述根据各个所述位移量对所述纵波和横波之间进行小层拉平处理，包括：
根据各个所述拉平种子点的移动量建立拉平后坐标，并将原始的PS数据根据该拉平后坐标进行差值，而后返回重新抽取一道PP波和一道PS波后，重复执行所述以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据，以及依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理的迭代过程，直至消除所述PS波数据对应的地震剖面的全部的所述局部畸变。
7.一种基于小层波形校正的多波地震资料匹配装置，其特征在于，包括：
初始匹配模块，用于针对同一目标地质区域的多波地震资料中的PP波数据和PS波数据进行初始匹配；
波形数据拉平模块，用于若相对于所述PP波数据对应的地震剖面，经初始匹配后的所述PS波数据对应的地震剖面出现局部畸变，则以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据；
局部畸变消除模块，用于依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理，以消除所述PS波数据对应的地震剖面的所述局部畸变。
8.根据权利要求7所述的基于小层波形校正的多波地震资料匹配装置，其特征在于，所述波形数据拉平模块包括：
抽取单元，用于分别在针对同一目标地质区域的多波地震资料中抽取一道PP波和一道PS波；
二维数组表示单元，用于应用二维数组表示抽取的所述PP波和PS波；
拉平种子点设置单元，用于根据预先获取的同相轴精细拉平处理参数中的波形对比窗口的大小及波形对比窗口的移动量，依次判断是否在某个所述波形对比窗口的中心点设置拉平种子点，若是，则通过设置该拉平种子点拉平所述PP波数据和PS波数据的横波同一位置的波形数据。
9.根据权利要求8所述的基于小层波形校正的多波地震资料匹配装置，其特征在于，所述二维数组表示单元包括：
道集表示子单元，用于将抽取的一道所述PP波和一道所述PS波存储成一个用于进行同相轴精细拉平处理的道集，并应用I行及J列的二维数组D来表示该道集；
其中，I用于表示每个地震道的采样点数，J用于表示道集的地震道数，且J＝2，在所述二维数组D中，第i行及第j列的值D(i，j)用于表示该道集的第j个地震道的第i个采样点的振幅值。
10.根据权利要求8所述的基于小层波形校正的多波地震资料匹配装置，其特征在于，所述同相轴精细拉平处理参数中还包含有所述波形对比窗口的搜索半径。
11.根据权利要求8至10任一项所述的基于小层波形校正的多波地震资料匹配装置，其特征在于，所述局部畸变消除模块包括：
移动量获取单元，用于依次获取所述波形对比窗口的一个面上的所述拉平种子点的移动量；
横向平滑处理单元，用于依次应用平滑函数对所述波形对比窗口的一个面上的所述拉平种子点的移动量进行横向平滑处理。
12.根据权利要求11所述的基于小层波形校正的多波地震资料匹配装置，其特征在于，所述局部畸变消除模块还包括：
迭代处理单元，用于根据各个所述拉平种子点的移动量建立拉平后坐标，并将原始的PS数据根据该拉平后坐标进行差值，而后返回重新抽取一道PP波和一道PS波后，重复执行所述以所述PP波数据和PS波数据的纵波道数据为基准，拉平所述PP波数据和PS波数据的横波同一位置的波形数据，以及依次获取被拉平的所述波形数据的位移量，并根据各个所述位移量对所述纵波和横波之间进行小层拉平处理的迭代过程，直至消除所述PS波数据对应的地震剖面的全部的所述局部畸变。
13.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至6任一项所述基于小层波形校正的多波地震资料匹配方法的步骤。
14.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至6任一项所述的基于小层波形校正的多波地震资料匹配方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910408871.XA CN111948707A (zh) | 2019-05-16 | 2019-05-16 | 基于小层波形校正的多波地震资料匹配方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910408871.XA CN111948707A (zh) | 2019-05-16 | 2019-05-16 | 基于小层波形校正的多波地震资料匹配方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111948707A true CN111948707A (zh) | 2020-11-17 |
Family
ID=73336685
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910408871.XA Pending CN111948707A (zh) | 2019-05-16 | 2019-05-16 | 基于小层波形校正的多波地震资料匹配方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111948707A (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4686657A (en) * | 1984-09-24 | 1987-08-11 | Mobil Oil Corporation | Matching P wave and shear wave sections in seismic exploration |
| US20100177595A1 (en) * | 2009-01-13 | 2010-07-15 | Vijay Khare | Using Seismic Attributes for Data Alignment and Seismic Inversion In Joint PP/PS Seismic Analysis |
| CN102879821A (zh) * | 2012-09-26 | 2013-01-16 | 中国石油天然气股份有限公司 | 一种针对地震叠前道集的同相轴精细拉平处理方法 |
| CN104865600A (zh) * | 2015-06-10 | 2015-08-26 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 将地震记录中的横波数据和纵波数据进行匹配的方法 |
| CN108897039A (zh) * | 2018-08-27 | 2018-11-27 | 电子科技大学 | 基于层位控制的dtw多波匹配方法 |
-
2019
- 2019-05-16 CN CN201910408871.XA patent/CN111948707A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4686657A (en) * | 1984-09-24 | 1987-08-11 | Mobil Oil Corporation | Matching P wave and shear wave sections in seismic exploration |
| US20100177595A1 (en) * | 2009-01-13 | 2010-07-15 | Vijay Khare | Using Seismic Attributes for Data Alignment and Seismic Inversion In Joint PP/PS Seismic Analysis |
| CN102879821A (zh) * | 2012-09-26 | 2013-01-16 | 中国石油天然气股份有限公司 | 一种针对地震叠前道集的同相轴精细拉平处理方法 |
| CN104865600A (zh) * | 2015-06-10 | 2015-08-26 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 将地震记录中的横波数据和纵波数据进行匹配的方法 |
| CN108897039A (zh) * | 2018-08-27 | 2018-11-27 | 电子科技大学 | 基于层位控制的dtw多波匹配方法 |
Non-Patent Citations (2)
| Title |
|---|
| NECATI GULUNAY ET AL.: "Gather flattening", 《THE LEADING EDGE》 * |
| 李涛: "多波高精度匹配方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106646611B (zh) | 一种纵波和转换波地震数据时间匹配的方法及装置 | |
| EA025919B1 (ru) | Системы и способы моделирования трехмерных геологических структур | |
| Haney et al. | Causal instrument corrections for short‐period and broadband seismometers | |
| US11768912B2 (en) | Performing multivariate time series prediction with three-dimensional transformations | |
| CN110082823B (zh) | 一种地震数据插值方法及装置 | |
| CN110295892B (zh) | 多极子阵列声波测井中横波衰减因子的确定方法及装置 | |
| US20240045093A1 (en) | Denoising methods and systems for preserving amplitude variation with offset (avo) features of pre-stack seismic data | |
| CN106324665A (zh) | 反演裂缝密度的方法及系统 | |
| CN110799856B (zh) | 使用与地下方位角和反射角相关的全波形反演产生速度模型 | |
| Piras et al. | Fast and realistic large-scale structure from machine-learning-augmented random field simulations | |
| CN104122584A (zh) | 根据地震数据确定方向性的方法及装置 | |
| CN111830559A (zh) | 一种地震资料拓频处理的方法及装置 | |
| CN107561588B (zh) | 一种地震数据噪声压制方法及装置 | |
| US11397273B2 (en) | Full waveform inversion in the midpoint-offset domain | |
| CN109100803A (zh) | 微断裂的确定方法和装置 | |
| CN111948707A (zh) | 基于小层波形校正的多波地震资料匹配方法及装置 | |
| Cui et al. | A tour of Marchenko redatuming: Focusing the subsurface wavefield | |
| US20230140656A1 (en) | Method and system for determining seismic processing parameters using machine learning | |
| CA2693247A1 (en) | Seismic attributes for structural analysis | |
| US10082588B2 (en) | Adaptive structure-oriented operator | |
| CN111948706B (zh) | 正交各向异性介质地震成像方法及装置 | |
| CN113721296A (zh) | 远震数据处理方法及装置 | |
| CN112083490B (zh) | 地震数据噪声衰减方法及装置 | |
| CN110824555B (zh) | 地震能量均衡方法及装置、计算机可读存储介质 | |
| Taufik et al. | A neural network based global traveltime function (GlobeNN) |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20201117 |