CN109471162A - 层间多次波处理方法、系统、电子设备及可读介质 - Google Patents
层间多次波处理方法、系统、电子设备及可读介质 Download PDFInfo
- Publication number
- CN109471162A CN109471162A CN201811166416.5A CN201811166416A CN109471162A CN 109471162 A CN109471162 A CN 109471162A CN 201811166416 A CN201811166416 A CN 201811166416A CN 109471162 A CN109471162 A CN 109471162A
- Authority
- CN
- China
- Prior art keywords
- wave field
- wave
- data
- vsp
- separation
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 16
- 238000000926 separation method Methods 0.000 claims abstract description 74
- 238000000034 method Methods 0.000 claims abstract description 45
- 238000012545 processing Methods 0.000 claims description 44
- 238000011144 upstream manufacturing Methods 0.000 claims description 31
- 238000012937 correction Methods 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 9
- 230000005284 excitation Effects 0.000 claims description 7
- 238000003860 storage Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 6
- 239000000284 extract Substances 0.000 claims description 5
- 239000011229 interlayer Substances 0.000 claims description 5
- 230000001419 dependent effect Effects 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 abstract description 7
- 238000007906 compression Methods 0.000 abstract description 4
- 238000003825 pressing Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 18
- 239000003245 coal Substances 0.000 description 8
- 238000001914 filtration Methods 0.000 description 5
- 230000006854 communication Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000012790 confirmation Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000010410 layer Substances 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
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
-
- 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
- G01V1/30—Analysis
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
本发明提供了层间多次波处理方法、系统、电子设备及可读介质。该方法包括：包括：分离VSP波场；根据地面地震数据和VSP波场的分离数据识别层间多次波；根据VSP波场的分离数据求取反褶积算子；利用该反褶积算子压制所述多次波。本发明通过井中地震和地面地震的地震数据结合，通过对VSP波场进行分离，从而结合地面数据比对识别层间多次波，并通过求取反褶积算子的方式，利用反褶积算子对多次波进行压制处理，从而使得多次波压制效果好，深部地层波阻特征更加清晰明显，使得地震数据更加真实反映地下地质的变化。
Description
技术领域
本发明涉及石油勘探领域，更具体的，涉及层间多次波处理方法、系统、电子设备及可读介质。
背景技术
多次波处理一直是地震资料处理的难题之一，多次波分为自由界面多次波(全程多次波)和层间多次波，自由界面多次波的识别与处理技术已经相当成熟，从开始的基于速度差异变化转化到不同域进行处理，到目前基于波动方程的多次波压制方法较少的依赖速度和地质信息，应用效果理想，已经成为地震处理的常规环节和主要方法；而层间多次波与一次波的时差、频率差异小识别困难，同时对地震信号的影响更加严重、隐蔽，对储层识别产生不可忽视的影响，随着地震勘探的不断深入，要求地震资料品质越来越高，目前缺少解决层间多次波的识别和压制问题的方法。
发明内容
为了解决层间多次波的识别和压制问题，本发明实施例提供一种层间多次波处理方法、系统、电子设备及可读介质。
在某些实施例中，一种层间多次波处理方法，包括：
分离VSP波场；
根据地面地震数据和VSP波场的分离数据识别层间多次波；
根据VSP波场的分离数据求取反褶积算子；
利用该反褶积算子压制所述多次波。
在某些实施例中，一种层间多次波处理系统，包括：
分离模块，分离VSP波场；
识别模块，根据地面地震数据和VSP波场的分离数据识别层间多次波；
求取模块，根据VSP波场的分离数据求取反褶积算子；
处理模块，利用该反褶积算子压制所述多次波。
在某些实施例中，一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上层间多次波处理方法的步骤。
在某些实施例中，一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如上层间多次波处理方法的步骤。
本发明通过井中地震(VSP)和地面地震的地震数据结合，通过对VSP波场进行分离，从而结合地面数据比对识别层间多次波，并通过求取反褶积算子的方式，利用反褶积算子对多次波进行压制处理，从而使得多次波压制效果好，深部地层波阻特征更加清晰明显，使得地震数据更加真实反映地下地质的变化。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1示出本发明实施例中一种层间多次波处理方法的流程示意图。
图2示出本发明实施例中VSP波场的观测图。
图3示出本发明实施例中VSP波场的示意图。
图4示出本发明实施例中VSP波场的模型建立图。
图5示出本发明实施例中VSP波场观测的系统定义的示意图。
图6示出图1中步骤S1的流程示意图。
图7示出本发明实施例中VSP波场的记录示意图。
图8示出本发明实施例中VSP波场下行波坡面图。
图9示出本发明实施例中VSP下行波提取算子图。
图10示出本发明实施例中正演记录去多次波前叠加剖面图。
图11示出本发明实施例中经过多次压制的剖面图。
图12示出图1中步骤S4的流程示意图。
图13示出现有技术地震资料速度谱图。
图14示出现有技术地震资料道集图。
图15示出现有技术地震资料叠加剖面图。
图16示出本发明实施例中VSP波场动校正后的剖面图。
图17示出本发明实施例中VSP波场去下行波后的剖面图。
图18示出本发明实施例中VSP波场反褶积后的剖面图。
图19示出本发明实施例中VSP波场标定示意图。
图20示出本发明实施例中VSP波场与地震道集进行匹配识别图。
图21示出进行层间多次波处理方法前的VSP波场剖面。
图22示出进行层间多次波处理方法后的VSP波场剖面。
图23示出本发明实施例中一种层间多次波处理系统结构示意图。
图24示出图23中分离模块的结构示意图。
图25示出图23中处理模块的结构示意图。
图26示出适用于用来实现本发明实施例的终端设备或服务器的计算机设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
随着地震勘探的不断深入，要求地震资料品质越来越高，目前缺少解决层间多次波的识别和压制问题的方法。
多次波压制方法通常分为两大类：一类是基于信号分析处理的滤波方法，如τ-р变换，f-k变换和预测反褶积等；另一类是基于波动方程的预测减去法，如波场外推法、逆散射级数方法等等。第一类方法由于考虑多次波和一次波速度差异性，相似的动校正速度减弱了在τ-р域变换的有效性，多次波和一次波具有相似的频率和振幅，减弱预测滤波的有效性；多次波和一次波具有相似的斜率，减弱了f-k滤波的有效性。第二类方法利用多次波产生机理来进行预测和压制多次波，考虑了多次波的传播运动学和动力学特征，目前方法还处于理论研究阶段，工业生产上还未看到可操作的流程方法。
以下论述涉及的符号表示：VSP表示井中地震；x(k)表示观测信号，s(k)表示信号；n(k)表示噪声；s'(k)表示输出信号；e(k)表示s(k)与s'(k)的误差；E表示均分误差；J表示h(k)的函数，当J最小时，h(k)为反褶积算子。
图1示出本发明实施例中的层间多次波处理方法，结合图1可知，包括：
S1:分离VSP波场；
S2:根据地面地震数据和VSP波场的分离数据识别层间多次波；
S3:根据VSP波场的分离数据求取反褶积算子；
S4:利用该反褶积算子压制所述多次波。
本发明通过井中地震(VSP)和地面地震的地震数据结合，通过对VSP波场进行分离，从而结合地面数据比对识别层间多次波，并通过求取反褶积算子的方式，利用反褶积算子对多次波进行压制处理，从而使得多次波压制效果好，深部地层波阻特征更加清晰明显，使得地震数据更加真实反映地下地质的变化。
结合图2和图3，VSP在井中观测可以避免或减少地面上的干扰，VSP观测时检波器在井中不同深度接收，为波场赋予了深度域信息，尤其是在井眼周围，反射波的深度和对应的地质层位可以准确标定，地层信息更易于估计，这也是地面地震所不具备的。地面的各种面波、散射和人文噪声很难避免。所以，前者更易于记录和识别各种地震波场(例如，通过图3中的采集区的波场特性，可以更加精确地预测预测区的波场特性)。
在具体实施例中，所述分离VSP波场，包括：将VSP波场分离为上行波场和下行波场，获取上行波场的波场分离数据和下行波场的波场分离数据，其中，所述上行波场的波场分离数据和下行波场的波场分离数据共同组成VSP波场的分离数据。
图4和图5为具体实施例中模型和观测系统定义的示意图，激发点置于地表自-1000m至1000m，激发点间隔100m，共201个激发点。接收点分为两部分，地面部分与激发点位置完全重合，井下部分自0m至1500m，接收点间隔100m，共151个接收点。需要说明的是，图中未示出所有激发点和接收点，图中的激发点和/或接收点的数目仅仅是示意性的。
可以通过非线性方式对VSP波场数据进行分离，生成上行波场和下行波场，本发明对分离VSP波场所用的方法不做限制，例如，可以采用常规的众志滤波、F-K滤波的方法，但上述方法存在混波、低频化、振幅畸变等缺陷，作为改进，本发明实施例中可以采用以下方法对VSP波场分离，结合图6，包括：
S11:对VSP波场的数据进行预处理；
S12:对预处理后的VSP波场的数据进行振幅衰减补偿；
S13:对补偿后的VSP波场的数据进行非线性分离，生成上行波场和下行波场；
S14:对分离后的上行波场和下行波场进行动校正NMO处理，得到上行波场的波场分离数据和下行波场的波场分离数据。
图7示出了VSP波场记录图，图8示出了VSP波场分离后的下行波剖面图，显然，通过上述分离方法，能够达到高保真波场分离(图8中的下行波部分最大程度地与图7的原下行波部分一致)，从而提高了数据的准确性。
在某些实施例中，对地面地震数据、所述上行波场的波场分离数据和下行波场的波场分离数据进行时间和振幅标定，判断层间多次波出现的位置。从而将VSP波场数据和地面地震数据结合识别出层间多层波出现的位置，使得数据更准确。
进一步的，如图9所示，根据下行波场的波场分离数据，提取对应不同炮点的统计子波，将其中一个子波作为期望输出，建立反褶积计算方程组，求取反褶积算子。
优选地，作为期望输出的子波为雷克子波，雷克子波作为地震子波中的一种，具有广泛的代表性和通用性。
进一步的，所述反褶积计算方程组为Toeplitz方程组。其处理过程如下：
设观测信号x(k)是由信号s(k)与噪声n(k)相加组成,二者的自相关与互相关都是已知的：
x(k)＝s(k)+n(k) (1)
设输出信号s'(k)：
记s(k)与s'(k)的误差为：
e(k)＝s(k)-s'(k) (3)
均方误差最小,即：
J＝E[e2(k)]＝min (4)
显然J是h(k)的函数：
上述问题就是求h(k)＝？时,J最小，h(k)即求取的反褶积算子。
基于VSP下行波反褶积方法，将反褶积算子应用到VSP波场数据过程中，得到如图10和图11所示的最终反演结果。与图7、图8相对比，可见多次波压制彻底，确认达到最佳反褶积处理效果。
结合图12，上述步骤中，利用该反褶积算子压制所述多次波，包括：
S21:抽取地面地震共中心点道集，并与VSP波场的数据进行一致性处理；
S22:将VSP波场的不同炮点的记录与求取的反褶积算子作反褶积处理，得到VSP波场反褶积后压制至少一次的数据；
S23:对地面地震数据不同的炮点的记录与求取的反褶积算子作相关处理，以消除炮点子波的非一致性，并压制所述多次波。
优选地，可以通过对反褶积后的道集叠加并输出数据，与叠加前的数据和VSP波场的走廊剖面比对，检测所述压制所述多次波的效果。
进一步的，为了提高压制精度，实现高保真压制，可以对地面地震叠前道集进行振幅一致性处理、速度解释和动校正NMO处理。
结合国内某区块三维数据处理实例，结合图13-图15，在该区块所在工区位置沙漠腹地，地表环境以起伏沙丘为主，工区内最大高程差接近100m，且大部分区域沙层厚度超过200m，深部地层中存一套强反射煤层，从以往资料来看，该区域地震资料存在较强多次波，其多次波来源于低速度沙层底界面，同时由于煤层的屏蔽作用，深度特征比较不清晰，从速度谱和道集分析，煤下目的层多次波发育，影响速度分析和深层精确成像。
通过上述实施例所述的处理方法，结合图16-图18，VSP数据经过动校正去掉下行波后可以明显观察到在煤层下面强烈发育一套层间多次波，通过下行波提取反褶积算子，经过反褶积运算处理过的数据可以看到层间多次波得到明显的衰减和消除。
如图19和图20所示，通过反褶积压制多次波后的VSP数据，可以看到在含有多次波污染的区域和正常地层反射区域存在明显的波组特征差异，与地面地震道集进行匹配标定识别，在煤层反射动校正拉平后，在煤层下面存在与正常反射波存在时差的低速反转现象，这组低速反射与VSP多次污染区相吻合，说明地震反射存在明显的多次波干扰。
结合图21和图22，对比可知，在基于井地联合识别和处理多次波的剖面中，多次波被压制效果较好，煤层以下的在深部地层波阻特征更加清晰明显，衰减地震假象给油气藏认识带来的影响，使得地震数据更加真实反映地下地质的变化。
图23示出了本发明实施例中层间多次波处理系统，包括：分离模块101，分离VSP波场；识别模块102，根据地面地震数据和VSP波场的分离数据识别层间多次波；求取模块103，根据VSP波场的分离数据求取反褶积算子；处理模块104，利用该反褶积算子压制所述多次波。本发明通过井中地震(VSP)和地面地震的地震数据结合，通过对VSP波场进行分离，从而结合地面数据比对识别层间多次波，并通过求取反褶积算子的方式，利用反褶积算子对多次波进行压制处理，从而使得多次波压制效果好，深部地层波阻特征更加清晰明显，使得地震数据更加真实反映地下地质的变化。
所述分离模块101配置为：将VSP波场分离为上行波场和下行波场，获取上行波场的波场分离数据和下行波场的波场分离数据，其中，所述上行波场的波场分离数据和下行波场的波场分离数据共同组成VSP波场的分离数据。
可以通过非线性方式对VSP波场数据进行分离，生成上行波场和下行波场，本发明对分离VSP波场所用的方法不做限制，例如，可以采用常规的众志滤波、F-K滤波的方法，但上述方法存在混波、低频化、振幅畸变等缺陷，作为改进，本发明实施例中可以采用以下方法对VSP波场分离，结合图24可知，所述分离单元111包括：
预处理单元1111，对VSP波场的数据进行预处理；
衰减补偿单元1112，对预处理后的VSP波场的数据进行振幅衰减补偿；
波场分离单元1113，对补偿后的VSP波场的数据进行非线性分离，生成上行波场和下行波场；
动校正单元1114，对分离后的上行波场和下行波场进行动校正NMO处理，得到上行波场的波场分离数据和下行波场的波场分离数据。
通过上述分离方法，能够达到高保真波场分离，从而提高了数据的准确性。
在某些实施例中，所述识别模块进一步包括：标定单元112，对地面地震数据、所述上行波场的波场分离数据和下行波场的波场分离数据进行时间和振幅标定，判断层间多次波出现的位置。从而将VSP波场数据和地面地震数据结合识别出层间多层波出现的位置，使得数据更准确。
进一步的，所述求取模块配置为：
根据下行波场的波场分离数据，提取对应不同炮点的统计子波，将其中一个子波作为期望输出，建立反褶积计算方程组，求取反褶积算子。
优选地，作为期望输出的子波为雷克子波，雷克子波作为地震子波中的一种，具有广泛的代表性和通用性。
进一步的，所述反褶积计算方程组为Toeplitz方程组。其处理过程如下：
设观测信号x(k)是由信号s(k)与噪声n(k)相加组成,二者的自相关与互相关都是已知的：
x(k)＝s(k)+n(k) (1)
设输出信号s'(k)：
记s(k)与s'(k)的误差为：
e(k)＝s(k)-s'(k) (3)
均方误差最小,即：
J＝E[e2(k)]＝min (4)
显然J是h(k)的函数：
上述问题就是求h(k)＝？时,J最小，h(k)即求取的反褶积算子。
基于VSP下行波反褶积方法，将反褶积算子应用到VSP波场数据过程中，得到如图10和图11所示的最终反演结果。与图7、图8相对比，可见多次波压制彻底，确认达到最佳反褶积处理效果。
在某些实施例中，结合图25，所述处理模块包括：
一致性处理单元121，抽取地震共中心点道集，并与VSP波场数据，进行一致性处理；
反褶积处理单元122，将VSP波场的不同炮点记录与求取的反褶积算子作反褶积处理，得到VSP波场反褶积后压制至少一次的数据；
相关处理单元123，对地面地震数据不同的炮点记录与求取的反褶积算子作相关处理，以消除炮点子波非一致性，并压制所述多次波。
优选地，处理模块可以包括比对单元124，从而实现对反褶积后的道集叠加并输出数据，与叠加前的数据和VSP波场的走廊剖面比对，检测所述压制所述多次波的效果。
进一步的，为了提高压制精度，实现高保真压制，所述处理模块进一步包括：
地面地震处理单元125，对地面地震叠前道集进行振幅一致性处理、速度解释和动校正NMO处理。
在上述国内某区块三维数据处理的实施例中，可以通过对比图21和图22知晓，在基于井地联合识别和处理多次波的剖面中，多次波被压制效果较好，煤层以下的在深部地层波阻特征更加清晰明显，衰减地震假象给油气藏认识带来的影响，使得地震数据更加真实反映地下地质的变化，本发明不予赘述。
进一步的，本发明的一些具体实施例提供一种计算机设备，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上所述的由客户端执行的方法，或者，所述处理器执行所述程序时实现如上所述的由服务器执行的方法。
下面参考图26，其示出了适于用来实现本申请实施例的终端设备或服务器的计算机设备600的结构示意图。
如图26所示，计算机设备600包括中央处理单元(CPU)601，其可以根据存储在只读存储器(ROM)602中的程序或者从存储部分608加载到随机访问存储器(RAM))603中的程序而执行各种适当的工作和处理。在RAM603中，还存储有系统600操作所需的各种程序和数据。CPU601、ROM602、以及RAM603通过总线604彼此相连。输入/输出(I/O)接口605也连接至总线604。
以下部件连接至I/O接口605：包括键盘、鼠标等的输入部分606；包括诸如阴极射线管(CRT)、液晶显示器(LCD)等以及扬声器等的输出部分607；包括硬盘等的存储部分608；以及包括诸如LAN卡，调制解调器等的网络接口卡的通信部分609。通信部分609经由诸如因特网的网络执行通信处理。驱动器610也根据需要连接至I/O接口606。可拆卸介质611，诸如磁盘、光盘、磁光盘、半导体存储器等等，根据需要安装在驱动器610上，以便于从其上读出的计算机程序根据需要被安装如存储部分608。
特别地，根据本发明的实施例，上文参考流程图描述的过程可以被实现为计算机软件程序。例如，本发明的实施例包括一种计算机程序产品，其包括有形地包含在机器可读介质上的计算机程序，所述计算机程序包括用于执行流程图所示的方法的程序代码。在这样的实施例中，该计算机程序可以通过通信部分609从网络上被下载和安装，和/或从可拆卸介质611被安装。
附图中的流程图和框图，图示了按照本发明各种实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分，所述模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发送。例如两个接连地表示的方框实际上可以基本并行地执行，他们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或操作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
Claims (22)
1.一种层间多次波处理方法，其特征在于，包括：
分离VSP波场；
根据地面地震数据和VSP波场的分离数据识别层间多次波；
根据VSP波场的分离数据求取反褶积算子；
利用该反褶积算子压制所述多次波。
2.根据权利要求1所述的方法，其特征在于，所述分离VSP波场，包括：
将VSP波场分离为上行波场和下行波场，获取上行波场的波场分离数据和下行波场的波场分离数据，其中，所述上行波场的波场分离数据和下行波场的波场分离数据共同组成VSP波场的分离数据。
3.根据权利要求2所述的方法，其特征在于，所述将VSP波场分离为上行波场和下行波场，获取上行波场的波场分离数据和下行波场的波场分离数据，包括：
对VSP波场的数据进行预处理；
对预处理后的VSP波场的数据进行振幅衰减补偿；
对补偿后的VSP波场的数据进行非线性波场分离，生成上行波场和下行波场；
对分离后的上行波场和下行波场进行动校正NMO处理，得到上行波场的波场分离数据和下行波场的波场分离数据。
4.根据权利要求2所述的方法，其特征在于，所述根据地面地震数据和VSP波场的分离数据识别层间多次波，包括：
对地面地震数据、所述上行波场的波场分离数据和下行波场的波场分离数据进行时间和振幅标定，判断层间多次波出现的位置。
5.根据权利要求2所述的方法，其特征在于，所述根据VSP波场的分离数据求取反褶积算子，包括：
根据下行波场的波场分离数据，提取对应不同炮点的统计子波，将其中一个子波作为期望输出，建立反褶积计算方程组，求取反褶积算子。
6.根据权利要求5所述的方法，其特征在于，作为期望输出的子波为雷克子波。
7.根据权利要求5所述的方法，其特征在于，所述反褶积计算方程组包括：
x(k)＝s(k)+n(k)
e(k)＝s(k)-s'(k)
J＝E[e2(k)]＝min
其中，x(k)表示观测信号，s(k)表示信号，n(k)表示噪声，s'(k)表示输出信号，e(k)表示s(k)与s'(k)的误差，E表示均分误差，J表示h(k)的函数，当J最小时，h(k)为反褶积算子。
8.根据权利要求1所述的方法，其特征在于，所述利用该反褶积算子压制所述多次波，包括：
抽取地面地震共中心点道集，并与VSP波场的数据进行一致性处理；
将VSP波场的不同震源激发点的记录与求取的反褶积算子作反褶积处理，得到VSP波场反褶积后压制至少一次的数据；
对地面地震数据不同的炮点的记录与求取的反褶积算子作相关处理，以消除炮点子波的非一致性，并压制所述多次波。
9.根据权利要求8所述的方法，其特征在于，所述利用该反褶积算子压制所述多次波，进一步包括：
对反褶积后的共面元中心点集叠加并输出数据，与叠加前的数据和VSP波场的走廊剖面比对，检测所述压制所述多次波的效果。
10.根据权利要求8所述的方法，其特征在于，所述利用该反褶积算子压制所述多次波，进一步包括：
对地面地震叠前道集进行振幅一致性处理、速度解释和动校正NMO处理。
11.一种层间多次波处理系统，其特征在于，包括：
分离模块，分离VSP波场；
识别模块，根据地面地震数据和VSP波场的分离数据识别层间多次波；
求取模块，根据VSP波场的分离数据求取反褶积算子；
处理模块，利用该反褶积算子压制所述多次波。
12.根据权利要求11所述的系统，其特征在于，所述分离模块配置为：
将VSP波场分离为上行波场和下行波场，获取上行波场的波场分离数据和下行波场的波场分离数据，其中，所述上行波场的波场分离数据和下行波场的波场分离数据共同组成VSP波场的分离数据。
13.根据权利要求12所述的系统，其特征在于，所述分离单元包括：
预处理单元，对VSP波场的数据进行预处理；
衰减补偿单元，对预处理后的VSP波场的数据进行振幅衰减补偿；
波场分离单元，对补偿后的VSP波场的数据进行非线性分离，生成上行波场和下行波场；
动校正单元，对分离后的上行波场和下行波场进行动校正NMO处理，得到上行波场的波场分离数据和下行波场的波场分离数据。
14.根据权利要求12所述的系统，其特征在于，所述识别模块包括：
标定单元，对地面地震数据、所述上行波场的波场分离数据和下行波场的波场分离数据进行时间和振幅标定，判断层间多次波出现的位置。
15.根据权利要求12所述的系统，其特征在于，所述求取模块配置为：
根据下行波场的波场分离数据，提取对应不同炮点的统计子波，将其中一个子波作为期望输出，建立反褶积计算方程组，求取反褶积算子。
16.根据权利要求15所述的系统，其特征在于，作为期望输出的子波为雷克子波。
17.根据权利要求15所述的系统，其特征在于，所述反褶积计算方程组包括：
x(k)＝s(k)+n(k)
e(k)＝s(k)-s'(k)
J＝E[e2(k)]＝min
其中，x(k)表示观测信号，s(k)表示信号，n(k)表示噪声，s'(k)表示输出信号，e(k)表示s(k)与s'(k)的误差，E表示均分误差，J表示h(k)的函数，当J最小时，h(k)为反褶积算子。
18.根据权利要求11所述的系统，其特征在于，所述处理模块包括：
一致性处理单元，抽取地震共中心点道集，并与VSP波场数据，进行一致性处理；
反褶积处理单元，将VSP波场的不同炮点记录与求取的反褶积算子作反褶积处理，得到VSP波场反褶积后消除多次波后的数据；
相关处理单元，对地面地震数据不同的炮点的记录与求取的反褶积算子作相关处理，以消除炮点子波的非一致性，并压制所述多次波。
19.根据权利要求18所述的系统，其特征在于，所述处理模块进一步包括：
比对单元，对反褶积后的道集叠加并输出数据，与叠加前的数据和VSP波场的走廊剖面比对，检测所述压制所述多次波的效果。
20.根据权利要求18所述的系统，其特征在于，所述处理模块进一步包括：
地面地震处理单元，对地面地震叠前道集进行振幅一致性处理、速度解释和动校正NMO处理。
21.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至10任一项所述的层间多次波处理方法的步骤。
22.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至10任一项所述的层间多次波处理方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811166416.5A CN109471162A (zh) | 2018-10-08 | 2018-10-08 | 层间多次波处理方法、系统、电子设备及可读介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811166416.5A CN109471162A (zh) | 2018-10-08 | 2018-10-08 | 层间多次波处理方法、系统、电子设备及可读介质 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109471162A true CN109471162A (zh) | 2019-03-15 |
Family
ID=65663358
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811166416.5A Pending CN109471162A (zh) | 2018-10-08 | 2018-10-08 | 层间多次波处理方法、系统、电子设备及可读介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109471162A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111352156A (zh) * | 2020-03-09 | 2020-06-30 | 吉林大学 | 多次波无关的地震资料自动速度分析方法 |
| CN117214950A (zh) * | 2023-09-25 | 2023-12-12 | 中海石油(中国)有限公司上海分公司 | 一种多次波压制方法、装置、设备及存储介质 |
| US11880012B2 (en) | 2020-01-20 | 2024-01-23 | China National Petroleum Corporation | Method and apparatus for extracting downgoing wavelet and attenuation parameters by using vertical seismic data |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2013152468A1 (zh) * | 2012-04-13 | 2013-10-17 | 中国石油天然气集团公司 | 一种地层品质因子反演方法 |
| CN103376466A (zh) * | 2012-04-13 | 2013-10-30 | 中国石油化工股份有限公司 | 一种多次波压制方法 |
| CN104502974A (zh) * | 2014-12-31 | 2015-04-08 | 中国石油天然气集团公司 | 一种压制多次波的组合方法及装置 |
| CN106094027A (zh) * | 2016-06-01 | 2016-11-09 | 中国海洋石油总公司 | 一种垂直地震剖面vsp钻前压力预测方法和系统 |
| CN107678062A (zh) * | 2017-09-15 | 2018-02-09 | 上海海洋大学 | 双曲Radon域综合预测反褶积和反馈循环方法压制多次波模型构建方法 |
| CN108107485A (zh) * | 2017-11-17 | 2018-06-01 | 中国石油天然气集团公司 | 地震数据多次波污染程度的评估方法及装置 |
-
2018
- 2018-10-08 CN CN201811166416.5A patent/CN109471162A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2013152468A1 (zh) * | 2012-04-13 | 2013-10-17 | 中国石油天然气集团公司 | 一种地层品质因子反演方法 |
| CN103376466A (zh) * | 2012-04-13 | 2013-10-30 | 中国石油化工股份有限公司 | 一种多次波压制方法 |
| CN104502974A (zh) * | 2014-12-31 | 2015-04-08 | 中国石油天然气集团公司 | 一种压制多次波的组合方法及装置 |
| CN106094027A (zh) * | 2016-06-01 | 2016-11-09 | 中国海洋石油总公司 | 一种垂直地震剖面vsp钻前压力预测方法和系统 |
| CN107678062A (zh) * | 2017-09-15 | 2018-02-09 | 上海海洋大学 | 双曲Radon域综合预测反褶积和反馈循环方法压制多次波模型构建方法 |
| CN108107485A (zh) * | 2017-11-17 | 2018-06-01 | 中国石油天然气集团公司 | 地震数据多次波污染程度的评估方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 蔡志东 等: "井地联合地震数据反褶积", 《石油地球物理勘探》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11880012B2 (en) | 2020-01-20 | 2024-01-23 | China National Petroleum Corporation | Method and apparatus for extracting downgoing wavelet and attenuation parameters by using vertical seismic data |
| CN111352156A (zh) * | 2020-03-09 | 2020-06-30 | 吉林大学 | 多次波无关的地震资料自动速度分析方法 |
| CN111352156B (zh) * | 2020-03-09 | 2021-09-17 | 吉林大学 | 多次波无关的地震资料自动速度分析方法 |
| CN117214950A (zh) * | 2023-09-25 | 2023-12-12 | 中海石油(中国)有限公司上海分公司 | 一种多次波压制方法、装置、设备及存储介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101334483B (zh) | 一种在地震数据处理中衰减瑞雷波散射噪声的方法 | |
| CN105334535B (zh) | 一种识别薄储层隐蔽岩性油气藏的方法 | |
| US6041018A (en) | Method for correcting amplitude and phase differences between time-lapse seismic surveys | |
| US6987706B2 (en) | Removal of noise from seismic data using high resolution radon transformations | |
| US7616524B1 (en) | Wavelet based intercept attribute for seismic exploration | |
| CN107356965B (zh) | 基于加权叠加噪音压制策略的反射系数反演储层预测方法 | |
| WO2020033465A1 (en) | Surface wave estimation and removal from seismic data | |
| CN105116449B (zh) | 一种弱反射储层的识别方法 | |
| CN109471162A (zh) | 层间多次波处理方法、系统、电子设备及可读介质 | |
| Witten et al. | Microseismic image-domain velocity inversion: Marcellus Shale case study | |
| CN105938203A (zh) | 一种储层特性的检测方法及装置 | |
| Malagnini et al. | Estimating absolute site effects | |
| Xie et al. | An improved ant-tracking workflow based on divided-frequency data for fracture detection | |
| CN111077577B (zh) | 井地联合储层描述方法及装置 | |
| CN107561578B (zh) | 一种针对宽方位地震资料的面波压制方法 | |
| CN107678065B (zh) | 提高地震分辨率的保构造井控空间反褶积方法和装置 | |
| Piazza et al. | M-factorial kriging-an efficient aid to noisy seismic data interpretation | |
| Perico et al. | Fault characterization in a postsalt reservoir interval, Jubarte Field (Campos Basin), using seismic attributes and machine learning | |
| Kneller et al. | Integration of broadband seismic data into reservoir characterization workflows: A case study from the Campos Basin, Brazil | |
| CN105989240A (zh) | 砂层预测方法及装置 | |
| Durrani et al. | A pragmatic approach for prestack seismic data conditioning and quality inspection for quantitative seismic reservoir characterization | |
| CN106125138B (zh) | 一种叠前crp道集拓频处理的方法 | |
| CN104459788B (zh) | 一种确定地震数据采集参数的方法及系统 | |
| Johann et al. | Spectral decomposition reveals geological hidden features in the amplitude maps from a deep water reservoir in the campos basin | |
| CN103760599B (zh) | 一种微小断层检测方法及断层检测装置 |
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 |
Application publication date: 20190315 |
|
| WD01 | Invention patent application deemed withdrawn after publication |