一种水下垂直缆地震波干涉成像方法及装置
阅读说明：本技术 一种水下垂直缆地震波干涉成像方法及装置 (Underwater vertical cable seismic wave interference imaging method and device ) 是由 陈学磊 王祥春 于 2021-10-19 设计创作，主要内容包括：本发明公开了一种水下垂直缆地震波干涉成像方法,所述方法包括：建立地质模型；选取若干个炮点中的两个炮点位置A和位置B；水下设置若干个检波器,将同一检波器接收到的位置A和位置B的两炮数据进行互相关运算；重复进行每个检波器的互相关运算,直至完成检波器域循环；获得一组干涉处理后的数据记录,进行叠加后得到以位置A模拟作为虚源,位置B作为模拟检波器的单道记录；完成炮域计算和共炮点道集,最后得到地质成像。本发明利用直达波与多次波,而不是将多次波作为噪声在处理过程中被滤除,而作为有效信息参与处理过程；通过构建虚震源,可以扩大数据量,因此能够获得更高质量的剖面成像。(The invention discloses an underwater vertical cable seismic wave interference imaging method, which comprises the following steps: establishing a geological model; selecting two shot point positions A and B in the plurality of shot points; arranging a plurality of detectors underwater, and performing cross-correlation operation on two shot data of a position A and a position B received by the same detector; repeatedly performing cross-correlation operation on each detector until the detector domain cycle is completed; obtaining a group of data records after interference processing, and obtaining a single-channel record with a position A simulation as a virtual source and a position B as a simulation detector after superposition; and completing shot domain calculation and shot point gather sharing, and finally obtaining geological imaging. The invention utilizes direct waves and multiples, rather than filtering the multiples as noise in the processing process, and participates in the processing process as effective information; by constructing virtual sources, the amount of data can be enlarged, and therefore higher quality sectional imaging can be obtained.)
技术领域
本发明是关于地质勘探技术，特别是关于一种水下垂直缆地震波干涉成像方法及装置。
背景技术
随着勘探技术的发展，油气勘探的重点逐渐从陆上转向对海洋，众所周知，对于天然气水合物最为有效精确的勘探方法就是地震勘探。在海洋地震勘探中，拖缆、海底电缆、海底地震仪等多种海上地震观测系统的应用及各种海洋地震数据处理方法也因此而日臻成熟，然而这些成熟的观测系统对于一些特殊的区域探测能力有限，成像效果并不能达到最优，也存在一些不足。
近些年来，海洋垂直缆地震法Vertical Cable Seismic，即VCS，因其宽频，广角，易识别多次波等特征而受到了国内外地球物理学者们的广泛关注。这一技术通常用于局部构造探测及大倾角地层的成像，是其他地震探测方法的一种有效的补充，同时也是高精度海洋地震勘探中不可或缺的一个环节。
垂直缆最早被应用于美国海军反潜战略装备系统中，世界上最早的垂直缆地震勘探发生在1987年的墨西哥湾，共投放了三条垂直缆。首次垂直缆地震勘探的目的，是为了解决该海域下地质构造的岩性储层成像难题。1989年，垂直缆地震数据采集系统被用于美国路易斯安娜州水深约600米的近海中，其中投放了12条垂直缆，海面布设了180条炮线，每一炮都通过调查船上的遥测系统进行实时监控。1992年，为了对墨西哥湾水深超过1000m的地下岩层进行成像，在水中投放了6条垂直缆，与同地区的拖缆采集数据做了同样的叠前深度偏移成像处理技术，发现垂直缆数据的成像精度明显大于拖缆数据的成像精度。1997年，Andersonden等人同样在墨西哥湾分别用拖缆和垂直缆进行地震数据的采集，并对采集后的拖缆和垂直缆数据进行了同等条件下的处理，偏移技术均采用的是叠前深度偏移技术，通过控制变量的研究方法得知，垂直缆数据比拖缆数据能呈现出更复杂的微小地质构造信息。1998年，Guimaraes等人对盐丘地质构造模型进行了正演，并模拟出垂直缆地震采集系统对其接收。同时，在处理过程中联合了模拟出的垂直缆上下行波场，最终证明了垂直缆观测系统的一些优缺点。
随着垂直缆这项技术在国外地震勘探中不断取得进展，国内诸多学者近年也逐步开展了垂直缆地震勘探这项技术的试采及研究，我国在垂直缆仪器的研制方面也是刚刚起步。
1998年，黄必铭翻译了Ward的与垂直缆相关的文章，主要介绍了地震采集可行性的Fuji实验。2015年，何勇等人对垂直缆观测系统的射线照明度进行了分析，为最佳的垂直缆观测系统设计提供了依据。2015年，Bian等人通过实验测试，表明采用多尺度的全波形反演技术，低频数据和多垂直缆采集系统的组合，可以提供高分辨率和高保真度的地下构造模型。2015年，杨晓辉等人使用Radon变换的方法对垂直缆数据进行了上下行波场的分离。2018年，安振芳等人分别对单垂直缆和多垂直缆设计了相应的观测系统，并从采集角度出发为改善垂直缆成像效果提供了一系列有益的结论。目前，我国在垂直缆地震勘探方面虽仍处于起步阶段，广州海洋地质调查局在我国南海区域已经开展了垂直缆地震勘探的试采及相关垂直缆系统的研制。2019年，胡守旺使用镜像偏移方法，将实际VCS勘探数据中多次波作为有效波进行常规方法成像，旨在处理出较为精细的垂直缆成像剖面。
地震干涉法发展至今日，其主要有三种具体形式，相关型地震干涉法，卷积型地震干涉法和反卷积型地震干涉法。
1968年，Claerbout利用水平层状介质地下放炮得到的地震记录，将从底部来的透射信息进行自相关运算，得到了描述场源关系的格林函数，后来这个方法被命名为“日光成像”。Claerbout打开地震干涉法的大门，成为地球物理行业里地震干涉技术的先驱者。该结论不仅适用于自激自收的情况，也可以应用在两个不同的检波器记录到的地震信号做互相关运算，这样可得到一个新的地震信号。Claerbout教授推断，这一结论同样可以推广到三维各向异性的介质中，然而并没有给出证明。2000年，Schuster教授通过对偏移后的地震数据进行固定相位分析，对地下构造信息成像，2001年在其发表的一篇会议文章中正式使用地震波干涉法这一名词。2002年，为了证明了Claerbout的猜想，Wapenaar等人结合格林定理进行论证，Wapenaar的证明为地震干涉法的发展提供了数学与物理基础，此后，地震干涉技术未来的发展十分光明。2004年，Calvert等人提出了虚拟震源法，推动了相干成像技术在实际上的应用，这种方法可以将基准面VSP数据转化成单井剖面记录数据。随后一段时间，该方法逐渐用到各个领域，如被动源干涉成像，体波干涉成像和面波干涉成像等等。2009年，Schuster教授结合前人对地震干涉法成像的总结，完成Seismic Interferometry一书，书中对地震干涉法在VSP以及各种地震勘探观测系统有着一定的描写。
同时期，国内也有很多学者对地震干涉法进行深入研究，其中在天然地震学上的应用研究较早，也取得了较好的效果，吴世萍教授团队和陈国金教授团队等对虚震源法进行了研究学习。2009年，常旭在主动震源方面对地震相干偏移与数据自参照偏移的关系进行了讨论。2011年，吴世萍研究了基于虚源估计的复杂上覆地层下地震相干成像。同年，乔宝平讨论了波场分离对虚源成像质量的影响。陶毅等人在文章中建立实际地质模型，并总结出用于干涉处理的一个流程。2012年倪雪灿在其硕士学位论文中分析了地震干涉法在VSP中的应用原理，同时具体描述了三种干涉方法的数学物理原理，证明了应用地震干涉法可以实现波场的重构。由于在干涉法处理的过程中，多次波作为有效成分，因而与常规VSP一次波成像对比，地震干涉法处理扩大了地下的照明范围。常旭等人也在VSP领域中对地震相干偏移有所研究。
在地震干涉法成像的研究方面，2004年，Schuster将Claerbout的日光照明方法统一为干涉成像方法，并从实际应用的角度展示了该方法的优势所在。2006年，Fink研究了复杂环境条件下的逆时声波成像，探讨了逆时窗的长度以及逆时信号的频带宽度的作用。同年，Willis和Homby分别利用逆时声波原理和干涉成像原理，实现了Walkaway VSP数据的盐丘侧翼成像。2006年，Xiao则利用地震干涉法实现了转换横波的盐丘侧翼成像；Hohl和Mateeva通过地震干涉法实现了海底电缆数据的被动震源反射波的成像。Shragge利用实际监测到的远震数据体，实现了炮集记录下的干涉偏移；Berkhout和Vershuur利用地震干涉法将多次波噪音转换成一次波进行了偏移成像，取得了不错的效果。在2006年，Borcea研究了随机介质条件下的相干成像技术；Yu和Schuster研究了逆垂直地震剖面数据体IVSP的互相关道集的偏移成像；Zhou对共深度域地震数据体的逆时偏移结果和干涉偏移结果进行了对比分析，发现两者都能很好地避开上覆复杂地层的影响。2007年，He实现了垂直地震剖面VSP数据体中自由表面多次波的3D波动方程干涉偏移。Draganov利用给定排列下检波器记录到的环境噪音，提取出反射波信息，并对其进行偏移成像，取得较好的效果。
2016年，Olivier Carrière对墨西哥湾处某海域的实际数据分别进行了常规地震处理与地震干涉技术处理，在文章中，主要使用的是相关型地震干涉。通过对比两种处理方法得到的地下PP反射系数图对比发现，在海洋地震勘探中，地震干涉技术依然适用。
在海洋地震勘探中，由于检波器成本昂贵，数量相对于陆地勘探较少，接收的海底有用信息有限，因此，对于接收的信号来说要尽可能的充分利用。
在常规处理方法中，如图1为海洋垂直缆地震法勘探示意图，多次波在成像的过程中被认作为噪声而被滤除，然而，多次波中仍然含有大量有用信息，处理这部分数据仍然可以得到地下构造信息；同时，海洋垂直缆勘探在设计之初，检波器需要悬挂于一根缆上，其采集方式为：电缆底部悬挂重物，沉放到海底；顶部悬挂浮球，浮于海面；检波器等间隔或不等间隔挂载于电缆上。这种挂载方式在采集的过程中，受洋流的影响，检波器在水中的位置也是不精确的，因此后期处理也需要做二次定位等工作。
发明内容
本发明的目的在于提供一种水下垂直缆地震波干涉成像方法，其能够利用直达波与多次波，通过构建虚震源以扩大数据量，获得更高质量的剖面成像；成像过程简单快速，并且结果更为准确。
为实现上述目的，本发明提供了一种水下垂直缆地震波干涉成像方法，所述方法包括：
建立地质模型，正演得到共炮点道集；
选取若干个炮点中的两个炮点位置A和位置B；
水下设置若干个检波器，将同一检波器接收到的位置A和位置B的两炮数据进行互相关运算；重复进行每个检波器的互相关运算，直至完成检波器域循环；
获得一组干涉处理后的数据记录，进行叠加后得到以位置A模拟作为虚源，位置B作为模拟检波器的单道记录；
完成炮域计算；若未完成炮域计算则再次进行“垂直缆上设置若干个检波器，将同一检波器接收到的位置A和位置B的两炮数据进行互相关运算”步骤；
计算处理得到以位置A为虚源的共炮点道集；
计算处理后得到地质成像。
在本发明的一实施方式中，炮点设置在垂直缆竖直方向以外的水面上，检波器位于垂直缆上。
在本发明的一实施方式中，位置A和位置B的两炮数据包括直达波与多次波信号。
在本发明的一实施方式中，所述位置A模拟作为虚源，位置B作为模拟检波器的单道记录，即将位置B炮点到检波器的直达波信号检出，以位置A炮点信号经地层反射到位置B后再次反射到检波器的地震波信号消去直达波信号后作为单道记录的信号数据。
在本发明的一实施方式中，检波器等间隔或不等间隔挂载于垂直缆上；垂直缆在水中位置为水面100m以下。
在本发明的一实施方式中，炮点为10-350个；炮点间距为6-50m。
在本发明的一实施方式中，检波器为10-95个；相邻检波器的间距为5-30米。
水中检波器数量越多，出来的结果越精确，但是，实际勘探过程中，检波器数量有限，几个或者几十个为佳；本发明方法不仅可以处理检波器数量多的情况，对于检波器数量少依然可以成像。实际海洋勘探中，炮间距一般是12.5米，大的炮间距也依然适用；检波器这个在实际沉放的过程中可能会很深，根据实际水体深度而定。
本发明还提供了一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述程序时实现上述任一水下垂直缆地震波干涉成像方法的步骤。
本发明还提供了一种非暂态计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述任一水下垂直缆地震波干涉成像方法的步骤。
本发明还提供了一种计算机程序产品，包括计算机程序，所述计算机程序被处理器执行时实现上述任一水下垂直缆地震波干涉成像方法的步骤。
与现有技术相比，根据本发明的水下垂直缆地震波干涉成像方法将地震干涉法应用到海洋垂直缆地震勘探中，利用直达波与多次波，而不是将多次波作为噪声在处理过程中被滤除，而作为有效信息参与处理过程；通过构建虚震源，可以扩大数据量，因此能够获得更高质量的剖面成像；垂直缆上检波器挂于缆之上，受水流影响其位置不准确，常规方法在后期需要二次定位等工作，而本发明在处理数据时，不需要精确的观测系统即可进行计算成像，成像结果显示，具有本发明方法具有较高的准确性；与常规处理方法相比，本发明方法成像简单快速，并且结果相对准确，可以为后续处理解释等提供可靠参考。
附图说明
图1是根据本发明背景技术中海洋垂直缆地震法勘探示意图；
图2是根据本发明水下垂直缆地震波干涉成像方法信号路径示意图；
图3是根据本发明水下垂直缆地震波干涉成像方法的流程图；
图4是根据本发明一实施方式的垂直缆地震勘探模型示意图；
图5是根据本发明一实施方式的正演地震记录图；
图6是根据本发明一实施方式的波场分离后得到直达波地震记录图；
图7是根据本发明一实施方式的波场分离后得到直达波与多次拨地震记录图；
图8是根据本发明一实施方式的计算处理后得到的虚源记录图；
图9是根据本发明一实施方式的虚源记录处理后得到的时间域叠加剖面图；
图10是根据本发明一实施方式的一次原始地震波信号记录图；
图11是根据本发明一实施方式的波场分离后直达波信号图；
图12是根据本发明一实施方式的波场分离后多次波信号图；
图13是根据本发明一实施方式的水下垂直缆地震波干涉成像方法处理实际数据得到的虚源记录；
图14是根据本发明一实施方式的数据处理后虚源记录叠加剖面图；
图15是根据本发明一实施方式的CMP104拾取速度谱图；
图16是根据本发明一实施方式的虚源记录速度分析后建立的速度场图；
图17是根据本发明一实施方式的电子设备的结构示意图；
主要附图标记说明：
1701-处理器；1702-存储器；1703-通信总线；1704-通信接口；
其中：图10-图13中纵坐标为时间，单位为ms，横坐标为trace；图9和图14的横坐标为cmp，纵坐标为时间，单位为s。
具体实施方式
下面结合附图，对本发明的具体实施方式进行详细描述，但应当理解本发明的保护范围并不受具体实施方式的限制。
除非另有其它明确表示，否则在整个说明书和权利要求书中，术语“包括”或其变换如“包含”或“包括有”等等将被理解为包括所陈述的元件或组成部分，而并未排除其它元件或其它组成部分。
如图2至图17所示，根据本发明优选实施方式的水下垂直缆地震波干涉成像方法。
本发明水下垂直缆地震波干涉成像方法处理的过程中，通过构建虚震源，从而实现对数据的处理。图2是本发明垂直缆地震波干涉成像法在VCS勘探中的处理原理。首先根据互易定理，将VCS勘探中的炮检点进行互换，实际VCS勘探中，PQ是位于海水面下的炮点，S是位于垂直缆上的检波点，在图2中，首先将炮检点互换，S点激发后，P点最先接收到直达波信号，随后在海水—空气界面产生反射，波继续向下传播，遇见界面后发生发射，Q点的检波点接收到反射波，将炮检互换，PQ接收到的波场信息可认为是PQ分别激发，S点接收。将PQ接收到的直达波和多次波进行干涉法处理后，消除了共同路径PS段，即可得到P为炮点Q为检波点的虚源记录。
通过上述处理原理可知，水中检波器的位置在处理时并不影响虚源记录的生成，因此该方法在处理数据时，并不需要精确的观测系统信息，即检波器的位置的轻微的波动不影响数据的准确性；同时还可以简单快速成像，为后续精确处理解释提供参考。
在其他常规地震数据处理方法中，多次波作为噪声会被滤除，然而实际中，多次波中仍然包含有反映地下构造的有用信息。在图2中可以看出，地震干涉法在VCS勘探处理中，直达波与多次波进行处理后，得到新的虚源记录，在这一过程中，多次波作为主要对象参与其中的处理，实现了多次波中的信息充分利用。
根据图2中生成的虚源记录中，炮检点均位于同一基准面，在获取的原始数据中，炮点位于海水面，检波器是悬挂于水中的垂直缆上。同时，新生成的虚源记录中，检波点的位置与炮点位置相同，相对于原始数据来说，能够增加检波点的数据，增大数据量，对于后续的成像来提供更多的地层反射信息和数据。
本发明垂直缆地震干涉法处理流程：在处理的过程中，根据上述的原理方法，对相关处理进行编程处理，制定了以下的处理流程，实现对数据的处理，并对得到的虚源记录进行处理成像。水下垂直缆地震波干涉成像方法的流程图如图3所示：建立地质模型，正演得到共炮点道集；选取若干个炮点中的两个炮点位置A和位置B；水下设置若干个检波器，将同一检波器接收到的位置A和位置B的两炮数据进行互相关运算；重复进行每个检波器的互相关运算，直至完成检波器域循环；获得一组干涉处理后的数据记录，进行叠加后得到以位置A模拟作为虚源，位置B作为模拟检波器的单道记录；完成炮域计算；若未完成炮域计算则再次进行“垂直缆上设置若干个检波器，将同一检波器接收到的位置A和位置B的两炮数据进行互相关运算”步骤；计算处理得到以位置A为虚源的共炮点道集；计算处理后得到地质成像。
对共炮点道集进行地震波干涉法的处理得到虚源记录的处理流程如图3所示。实际VCS勘探数据为共接收点道集，因此前期需要对数据进行预处理工作。如果要得到一个以位置A为虚震源和位置B为检波点的新的地震记录，需要提取位置A为检波器对应于每一炮的地震道，以及对应于同一个炮记录的位置B为检波器所记录的地震信号，并对他们进行相关，直到完成所有炮对B检波器和A检波器的地震信号的提取和相关。对相关图谱所有的地震道的信号进行叠加将得到一个新的地震信号，该信号可解释为以位置A为震源位置B为检波器的地震信号。然后再进行以位置A为虚震源与其他炮点处检波器接收到的地震信号进行互相关运算随后叠加，直至完成所有炮点处的计算，即可得到一个以A为虚源点其他炮点处为虚检波点的新的虚源地震记录，最后，对所有炮点完成虚源点的计算，得到所有炮点对应的虚源记录。得到虚源记录后，即可使用相关的处理软件完成动校正、叠加等常规处理，实现成像。
本发明的实施例中，首先需要构建地质模型，是为了更好地验证地震干涉法在VCS勘探中的应用，首先建立了一个实际勘探地质模型，垂直缆地震勘探模型示意图如图4所示。在模型中，共计模拟四种层位，水平的海底和两个倾斜界面，速度依次是1500m/s、1800m/s、2100m/s和2400m/s，海面共放置炮点250个，水中模拟一条垂直缆，垂直缆在水中位置为水下200-600m，缆上共放置51个检波器，正演时记录长度3.5s，采样间隔1ms，正演过程中用到的子波为50Hz雷克子波。
经过正演后得到250个单炮记录，如图5正演地震记录图，为选取炮点位置为125时的单炮记录，从单炮记录中可以看出，接收到的信号包含有上下行波，因为干涉法用的主要为直达波和多次波，多为下行波，因此需要对数据进行上下行波分离。经过分离后得到直达波和多次波。
将所有的正演地震数据均进行分离后，然后按照图3所示的处理流程进行处理，经过地震干涉法处理后，垂直缆上的检波点重建到海面炮点处，得到炮检点均为于海面的虚源地震记录，如图8计算处理后得到的虚源记录图所示。
在本实施例中，经过地震干涉法处理后得到的虚源记录中，存在有部分计算噪声，因为正演过程中并未加入噪声，因此可以断定这些干扰是计算过程中产生的，因此在后续的处理中首先需要完成去噪工作。得到的虚源记录后，使用商业化处理软件，完成速度分析，动校正和叠加等常工作，处理得到该模型的叠加剖面。图9为虚源记录处理后得到的时间域叠加剖面图。在叠加剖面中可以看出，从上到下依次有三条强能量同相轴，依次为海底和两个倾斜界面。根据原始模型中各层速度和深度的关系，计算得到的双程旅行时与叠加剖面上时间相吻合，证明该地震干涉法在处理垂直缆地震勘探模型中的准确性，为实际数据处理奠定了理论试验基础。
本发明的一个具体实施方式为，在模型数据处理结果的基础上，完成实际VCS勘探数据的处理。本次实际数据采集于琼东南海域，该地震调查工区距离三亚市180千米，调查区水深约为1400-1600米，海底地形变化相对平缓。海洋垂直缆电缆长300米，缆上共有12个水听器，间距25米，采样率2ms，记录长度为8s。本工区共布设采集炮线58条，每条炮线共计1134炮，炮间距为12.5米。
由于实际数据中包含有多种噪音干扰，因此首先对实际数据做了去噪等预处理工作，随后将直达波和多次波进行分离。图10为垂直缆上一个检波器的一次原始地震波信号记录图。图11为波场分离后直达波信号图，12为波场分离后多次波信号图。
在实际处理的过程中，将垂直缆上的12个检波器数据按照图3所示的处理流程进行地震干涉法处理，并得到实际数据对应的虚源记录，图13为处理实际数据得到的虚源记录。对模型数据处理得到的虚源记录后，将得到的数据加载到GeoEast软件里，完成速度分析，动校正和叠加等常工作，处理得到该模型的叠加剖面。图15为CMP104拾取速度谱图，是速度分析过程中共中心点道集(CMP)的速度谱。图16为虚源记录速度分析后建立的速度场图，是速度分析过程中拾取速度谱建立的速度场，在小号CMP内，由于覆盖次数较少，因而拾得的速度不准确，随着覆盖次数的增大，从CMP120开始，速度逐渐趋于稳定，同时速度在时间的变化与模型中地层变化相吻合。在此速度分析的基础上，通过拾取的速度完成数据的动校正工作。
从图14数据处理后虚源记录叠加剖面图可以看出，直达波对应的时间在2000ms左右，海水中纵波速度大约为1500m/s，经过计算，海底界面约为1500m，与工区实际海水深度大致相同，成像结果较为准确。同时，该叠加剖面中，地层同相轴清晰且变化平稳，浅部地层成像效果较好，与实际情况相吻合。因此，在实际数据处理的过程中，该方法能够获得一个准确的结果，可以为后续的精确处理以及解释提供参考。
本发明实施例同时提供了一种电子设备，如图17是本发明实施例提供的一种电子设备的结构示意图，该电子设备可以包括：处理器(processor)1701、通信接口(Communications Interface)1704、存储器(memory)1702和通信总线1703，其中，处理器，通信接口，存储器通过通信总线完成相互间的通信。处理器可以调用存储器中的逻辑指令，以执行上述各实施例中提供的水下垂直缆地震波干涉成像方法，该方法包括：建立地质模型，正演得到共炮点道集；选取若干个炮点中的两个炮点位置A和位置B；水下设置若干个检波器，将同一检波器接收到的位置A和位置B的两炮数据进行互相关运算；重复进行每个检波器的互相关运算，直至完成检波器域循环；获得一组干涉处理后的数据记录，进行叠加后得到以位置A模拟作为虚源，位置B作为模拟检波器的单道记录；完成炮域计算；若未完成炮域计算则再次进行“垂直缆上设置若干个检波器，将同一检波器接收到的位置A和位置B的两炮数据进行互相关运算”步骤；计算处理得到以位置A为虚源的共炮点道集；计算处理后得到地质成像。
此外，上述的存储器1702中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
另一方面，本发明实施例还提供一种计算机程序产品，所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所述计算机程序包括程序指令，当所述程序指令被计算机执行时，计算机能够执行上述各实施例中提供的水下垂直缆地震波干涉成像方法，该方法包括：建立地质模型，正演得到共炮点道集；选取若干个炮点中的两个炮点位置A和位置B；水下设置若干个检波器，将同一检波器接收到的位置A和位置B的两炮数据进行互相关运算；重复进行每个检波器的互相关运算，直至完成检波器域循环；获得一组干涉处理后的数据记录，进行叠加后得到以位置A模拟作为虚源，位置B作为模拟检波器的单道记录；完成炮域计算；若未完成炮域计算则再次进行“垂直缆上设置若干个检波器，将同一检波器接收到的位置A和位置B的两炮数据进行互相关运算”步骤；计算处理得到以位置A为虚源的共炮点道集；计算处理后得到地质成像。
又一方面，本发明实施例还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各实施例提供的水下垂直缆地震波干涉成像方法，该方法包括：建立地质模型，正演得到共炮点道集；选取若干个炮点中的两个炮点位置A和位置B；水下设置若干个检波器，将同一检波器接收到的位置A和位置B的两炮数据进行互相关运算；重复进行每个检波器的互相关运算，直至完成检波器域循环；获得一组干涉处理后的数据记录，进行叠加后得到以位置A模拟作为虚源，位置B作为模拟检波器的单道记录；完成炮域计算；若未完成炮域计算则再次进行“垂直缆上设置若干个检波器，将同一检波器接收到的位置A和位置B的两炮数据进行互相关运算”步骤；计算处理得到以位置A为虚源的共炮点道集；计算处理后得到地质成像。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
前述对本发明的具体示例性实施方案的描述是为了说明和例证的目的。这些描述并非想将本发明限定为所公开的精确形式，并且很显然，根据上述教导，可以进行很多改变和变化。对示例性实施例进行选择和描述的目的在于解释本发明的特定原理及其实际应用，从而使得本领域的技术人员能够实现并利用本发明的各种不同的示例性实施方案以及各种不同的选择和改变。本发明的范围意在由权利要求书及其等同形式所限定。
- 上一篇：一种医用注射器针头装配设备
- 下一篇：一种考虑倾角时差的反射点位置确定方法