CN112558144B - 基于零偏vsp数据确定裂缝方位的方法、装置及电子设备 - Google Patents

基于零偏vsp数据确定裂缝方位的方法、装置及电子设备 Download PDF

Info

Publication number
CN112558144B
CN112558144B CN202011371093.0A CN202011371093A CN112558144B CN 112558144 B CN112558144 B CN 112558144B CN 202011371093 A CN202011371093 A CN 202011371093A CN 112558144 B CN112558144 B CN 112558144B
Authority
CN
China
Prior art keywords
component
wave
target
transverse wave
fast
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
Application number
CN202011371093.0A
Other languages
English (en)
Other versions
CN112558144A (zh
Inventor
李梦琦
芦俊
王赟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Geosciences Beijing
Original Assignee
China University of Geosciences Beijing
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Geosciences Beijing filed Critical China University of Geosciences Beijing
Priority to CN202011371093.0A priority Critical patent/CN112558144B/zh
Publication of CN112558144A publication Critical patent/CN112558144A/zh
Application granted granted Critical
Publication of CN112558144B publication Critical patent/CN112558144B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
    • G01V1/18—Receiving elements, e.g. seismometer, geophone or torque detectors, for localised single point measurements
    • G01V1/181—Geophones
    • G01V1/184—Multi-component geophones
    • 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
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures
    • 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
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • 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
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/51—Migration
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/65—Source localisation, e.g. faults, hypocenters or reservoirs

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地震数据进行处理，得到上行波在预设三维正交坐标系的水平方向上的R分量和T分量；基于预设波场分离方法，对R分量和T分量进行波矢分离处理，得到R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波；基于R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，确定目标夹角；基于目标夹角、三分量检波器在预设三维正交坐标系的水平方向上的R分量方位，确定目标裂缝的方位信息。

Description

基于零偏VSP数据确定裂缝方位的方法、装置及电子设备
技术领域
本发明涉及计算机技术领域，尤其涉及一种基于零偏VSP数据确定裂缝方位的方法、装置及电子设备。
背景技术
随着石油天然气勘探与开发进程的推进，常规的孔隙型油气藏储量日益减少，勘探方向已经逐步转向裂缝型油气藏，因而，利用地震数据预测地层中裂缝的方位成为油气勘探领域的重要研究方向。
目前，可以通过对横波分裂现象的分析，预测地层中裂缝的方位信息，例如，可以对多方位等径源距(Walkaround)的垂直地震剖面(Vertical Seismic Profiling，VSP)数据，进行矢量分解与合成，以确定地层中裂缝的方位信息。
Walkaround VSP数据携带了横波的方位各向异性信息，但是，地层中存在多个裂缝地层时，地震波在裂缝层处产生的分裂横波在上行过程中会受到上覆地层的影响，同时，分裂横波的层匹配也存在着误差，因而，通过Walkaround VSP数据确定裂缝的方位信息，存在方位信息确定准确性差的问题。
发明内容
本发明实施例的目的是提供一种基于零偏VSP数据确定裂缝方位的方法、装置及电子设备，以解决现有技术中在确定地层中裂缝的方位信息时，存在的方位信息确定准确性差的问题。
为解决上述技术问题，本发明实施例是这样实现的：
第一方面，本发明实施例提供的一种基于零偏VSP数据确定裂缝方位的方法，所述方法包括：
获取目标井的三分量检波器检测到的零偏移距VSP地震数据；
基于预设下行波去除算法，对所述零偏移距VSP地震数据进行处理，得到上行波在预设三维正交坐标系的水平方向上的R分量和T分量；
基于预设波场分离方法，对所述R分量和T分量进行波矢分离处理，得到R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波；
基于所述R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，确定目标夹角，所述目标夹角为由所述快横波在所述预设三维正交坐标系的水平方向上的方位和所述慢横波在所述预设三维正交坐标系的水平方向上的方位确定的夹角；
基于所述目标夹角、所述三分量检波器在所述预设三维正交坐标系的水平方向上的R分量方位，确定目标裂缝的方位信息。
可选地，所述基于所述R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，确定目标夹角，包括：
基于所述零偏移距VSP地震数据，确定纵波初至位置；
在所述R分量剖面和所述T分量剖面上，确定所述快横波和所述慢横波的目标相同起始点；
在所述R分量剖面上，基于所述纵波初至位置以及所述目标相同起始点，确定第一数据范围，并在所述T分量剖面上，基于所述纵波初至位置以及所述目标相同起始点，确定第二数据范围；
基于所述第一数据范围内的所述R分量上的快横波和R分量上的慢横波，以及所述第二数据范围内的所述T分量上的快横波和T分量上的快横波，确定所述目标夹角。
可选地，所述基于所述第一数据范围内的所述R分量上的快横波和R分量上的慢横波，以及所述第二数据范围内的所述T分量上的快横波和T分量上的快横波，确定所述目标夹角，包括：
基于预设第一滑动时窗，在所述R分量剖面上，对所述第一数据范围内的所述R分量上的快横波进行扫描处理，得到第一目标快横波；
基于所述预设第一滑动时窗，在所述R分量剖面上，对所述第一数据范围内的所述R分量上的慢横波进行扫描处理，得到第一目标慢横波；
基于预设第二滑动时窗，在所述T分量剖面上，对所述第二数据范围内的所述T分量上的快横波进行扫描处理，得到第二目标快横波；
基于所述预设第二滑动时窗，在所述T分量剖面上，对所述第二数据范围内的所述T分量上的慢横波进行扫描处理，得到第二目标慢横波；
基于所述第一目标快横波、所述第一目标慢横波、所述第二目标快横波、所述第二目标慢横波，确定所述目标夹角。
可选地，所述基于预设第一滑动时窗，在所述R分量剖面上，对所述第一数据范围内的所述R分量上的快横波进行扫描处理，得到第一目标快横波，包括：
将每个所述第一滑动时窗在所述R分量剖面内包含的R分量上的快横波波矢振幅的均值，作为每个所述第一滑动时窗对应的第一快横波波矢在R分量上的投影，其中，每两个相邻所述第一滑动时窗包含预设数量的相同的检波器；
基于每个所述第一滑动时窗对应的第一快横波波矢在R分量上的投影，确定所述第一目标快横波。
可选地，所述预设第一滑动时窗和所述第二滑动时窗为平行四边形滑动时窗，所述预设第一滑动时窗在所述R分量剖面上的斜边的斜率与所述R分量上的快横波的同相轴的斜率相同，所述预设第二滑动时窗在所述R分量剖面上的斜边的斜率与所述R分量上的慢横波的同相轴的斜率相同。
可选地，所述基于所述第一目标快横波、所述第一目标慢横波、所述第二目标快横波、所述第二目标慢横波，确定所述目标夹角，包括：
将所述第一目标快横波和所述第二目标快横波代入公式
Figure BDA0002806686830000031
得到第一夹角，其中，θ1为所述第一夹角，PS1(R)为所述第一目标快横波，PS1(T)为所述第二目标快横波；
将所述第一目标慢横波和所述第二目标慢横波代入公式
Figure BDA0002806686830000032
得到第二夹角，其中，θ2为所述第二夹角，PS2(R)为所述第一目标慢横波，PS2(T)为所述第二目标慢横波；
将所述第一夹角和所述第二夹角代入公式
θ＝θ1*f+θ2*(1-f)
得到所述目标夹角，其中，θ为所述目标夹角，θ1为所述第一夹角，θ2为所述第二夹角，f为预设加权系数。
可选地，所述基于所述目标夹角、所述三分量检波器在所述预设三维正交坐标系的水平方向上的R分量方位，确定目标裂缝的方位信息，包括：
将所述目标夹角和所述三分量检波器在所述预设三维正交坐标系的水平方向上的R分量方位的和，作为所述目标裂缝的方位角；
基于所述纵波初至位置，确定纵波单程旅行时；
基于所述目标裂缝的方位角和所述纵波单程旅行时，确定所述目标裂缝的方位信息。
第二方面，本发明实施例提供了一种基于零偏VSP数据确定裂缝方位的装置，所述装置包括：
数据获取模块，用于获取目标井的三分量检波器检测到的零偏移距VSP地震数据；
数据处理模块，用于基于预设下行波去除算法，对所述零偏移距VSP地震数据进行处理，得到上行波在预设三维正交坐标系的水平方向上的R分量和T分量；
波场分离模块，用于基于预设波场分离方法，对所述R分量和T分量进行波矢分离处理，得到R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波；
夹角确定模块，用于基于所述R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，确定目标夹角，所述目标夹角为由所述快横波在所述预设三维正交坐标系的水平方向上的方位和所述慢横波在所述预设三维正交坐标系的水平方向上的方位确定的夹角；
方位确定模块，用于基于所述目标夹角、所述三分量检波器在所述预设三维正交坐标系的水平方向上的R分量方位，确定目标裂缝的方位信息。
可选地，所述夹角确定模块，用于：
基于所述零偏移距VSP地震数据，确定纵波初至位置；
在所述R分量剖面和所述T分量剖面上，确定所述快横波和所述慢横波的目标相同起始点；
在所述R分量剖面上，基于所述纵波初至位置以及所述目标相同起始点，确定第一数据范围，并在所述T分量剖面上，基于所述纵波初至位置以及所述目标相同起始点，确定第二数据范围；
基于所述第一数据范围内的所述R分量上的快横波和R分量上的慢横波，以及所述第二数据范围内的所述T分量上的快横波和T分量上的快横波，确定所述目标夹角。
可选地，所述夹角确定模块，用于：
基于预设第一滑动时窗，在所述R分量剖面上，对所述第一数据范围内的所述R分量上的快横波进行扫描处理，得到第一目标快横波；
基于所述预设第一滑动时窗，在所述R分量剖面上，对所述第一数据范围内的所述R分量上的慢横波进行扫描处理，得到第一目标慢横波；
基于预设第二滑动时窗，在所述T分量剖面上，对所述第二数据范围内的所述T分量上的快横波进行扫描处理，得到第二目标快横波；
基于所述预设第二滑动时窗，在所述T分量剖面上，对所述第二数据范围内的所述T分量上的慢横波进行扫描处理，得到第二目标慢横波；
基于所述第一目标快横波、所述第一目标慢横波、所述第二目标快横波、所述第二目标慢横波，确定所述目标夹角。
可选地，所述夹角确定模块，用于：
将每个所述第一滑动时窗在所述R分量剖面内包含的R分量上的快横波波矢振幅的均值，作为每个所述第一滑动时窗对应的第一快横波波矢在R分量上的投影，其中，每两个相邻所述第一滑动时窗包含预设数量的相同的检波器；
基于每个所述第一滑动时窗对应的第一快横波波矢在R分量上的投影，确定所述第一目标快横波。
可选地，所述预设第一滑动时窗和所述第二滑动时窗为平行四边形滑动时窗，所述预设第一滑动时窗在所述R分量剖面上的斜边的斜率与所述R分量上的快横波的同相轴的斜率相同，所述预设第二滑动时窗在所述R分量剖面上的斜边的斜率与所述R分量上的慢横波的同相轴的斜率相同。
可选地，所述夹角确定模块，用于：
将所述第一目标快横波和所述第二目标快横波代入公式
Figure BDA0002806686830000051
得到第一夹角，其中，θ1为所述第一夹角，PS1(R)为所述第一目标快横波，PS1(T)为所述第二目标快横波；
将所述第一目标慢横波和所述第二目标慢横波代入公式
Figure BDA0002806686830000052
得到第二夹角，其中，θ2为所述第二夹角，PS2(R)为所述第一目标慢横波，PS2(T)为所述第二目标慢横波；
将所述第一夹角和所述第二夹角代入公式
θ＝θ1*f+θ2*(1-f)
得到所述目标夹角，其中，θ为所述目标夹角，θ1为所述第一夹角，θ2为所述第二夹角，f为预设加权系数。
可选地，所述方位确定模块，用于：
将所述目标夹角和所述三分量检波器在所述预设三维正交坐标系的水平方向上的R分量方位的和，作为所述目标裂缝的方位角；
基于所述纵波初至位置，确定纵波单程旅行时；
基于所述目标裂缝的方位角和所述纵波单程旅行时，确定所述目标裂缝的方位信息。
第三方面，本发明实施例提供一种电子设备，包括处理器、存储器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述计算机程序被所述处理器执行时实现上述第一方面提供的基于零偏VSP数据确定裂缝方位的方法的步骤。
第四方面，本发明实施例提供一种计算机可读存储介质，所述计算机可读存储介质上存储计算机程序，所述计算机程序被处理器执行时实现上述第一方面提供的基于零偏VSP数据确定裂缝方位的方法的步骤。
由以上本发明实施例提供的技术方案可见，本发明实施例通过获取目标井的三分量检波器检测到的零偏移距VSP地震数据，基于预设下行波去除算法，对零偏移距VSP地震数据进行处理，得到上行波在预设三维正交坐标系的水平方向上的R分量和T分量，基于预设波场分离方法，对R分量和T分量进行波矢分离处理，得到R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，基于R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，确定目标夹角，目标夹角为由快横波在预设三维正交坐标系的水平方向上的方位和慢横波在预设三维正交坐标系的水平方向上的方位确定的夹角，基于目标夹角、三分量检波器在预设三维正交坐标系的水平方向上的R分量方位，确定目标裂缝的方位信息。这样，由于通过预设下行波去除算法，对零偏移距VSP地震数据进行了处理，所以，可以避免上覆地层的影响以及层匹配的需求，在通过波场分离方法，将上行波中视速度不同的快横波和慢横波分离为R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，这样，基于分离后的数据，就可以准确的确定目标裂缝的方位信息，提高目标裂缝的方位信息的确定准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明一种基于零偏VSP数据确定裂缝方位的方法的流程示意图；
图2为本发明一种零偏VSP数据获取方法的示意图；
图3为本发明一种三维正交坐标系的示意图；
图4为本发明一种获取上行波在预设三维正交坐标系的水平方向上的R分量和T分量的示意图；
图5为本发明一种第一数据范围和第二数据范围的示意图；
图6为本发明一种目标夹角的确定方法的示意图；
图7为本发明一种基于零偏VSP数据确定裂缝方位的方法的流程示意图；
图8为本发明一种第一滑动时窗和第二滑动时窗的示意图；
图9为本发明一种上行波在预设三维正交坐标系的水平方向上的R分量和T分量的示意图；
图10为本发明一种目标裂缝的方位信息的示意图；
图11为本发明一种基于零偏VSP数据确定裂缝方位的装置的结构示意图；
图12为本发明一种电子设备的结构示意图。
具体实施方式
本发明实施例提供一种基于零偏VSP数据确定裂缝方位的方法、装置及电子设备。
为了使本技术领域的人员更好地理解本发明中的技术方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
实施例一
如图1所示，本说明书实施例提供一种基于零偏VSP数据确定裂缝方位的方法，该方法的执行主体可以为服务器，该服务器可以是独立的服务器，也可以是由多个服务器组成的服务器集群。该方法具体可以包括以下步骤：
在S102中，获取目标井的三分量检波器检测到的零偏移距VSP地震数据。
其中，目标井可以是通过相关的电子器件等进行野外采集的区域，三分量检波器可以是多波勘探时使用的特种检波器，三分量检波器与单分量的常规地震检波器不同，每个三分量检波器内装有三个互相垂直的传感器，以记录质点振动速度向量的三个分量(垂直Z分量、水平X以及水平Y分量)，用于同时记录纵波、横波、转换波，零偏移距VSP地震数据可以为通过相关的勘探器件进行实地勘探时记录的勘探资料中获取的已知数据，零偏移距VSP地震数据采集，是指在离目标井最近的一个炮点激发震源时，进行相应的VSP数据采集。
在实施中，随着石油天然气勘探与开发进程的推进，常规的孔隙型油气藏储量日益减少，勘探方向已经逐步转向裂缝型油气藏，因而，利用地震数据预测地层中裂缝的方位成为油气勘探领域的重要研究方向。
目前，可以通过对横波分裂现象的分析，预测地层中裂缝的方位信息，例如，可以对多方位等径源距(Walkaround)的垂直地震剖面(Vertical Seismic Profiling，VSP)数据，进行矢量分解与合成，以确定地层中裂缝的方位信息。
Walkaround VSP数据携带了横波的方位各向异性信息，但是，地层中存在多个裂缝地层时，地震波在裂缝层处产生的分裂横波在上行过程中会受到上覆地层的影响，同时，分裂横波的层匹配也存在着误差，因而，通过Walkaround VSP数据确定裂缝的方位信息，存在方位信息确定准确性差的问题。为此，本发明实施例提供一种能够解决上述问题的技术方案，具体可以参见下述内容：
如图2所示，在离目标井最近的一个震源(即与目标井井口之间的偏移距约等于零的震源)激发时，由震源激发的下行P波在地下反射界面处形成上行转换横波(即PSV波)，并在经过裂缝地层时会分裂为偏振方位平行于裂缝平面的快横波(PS1波)和垂直于裂缝平面的慢横波(PS2波)，即快横波和慢横波的偏振相互垂直。此时，目标井内的多个三分量检波器可以接收到PS1波和PS2波，即三分量检波器可以进行相应的VSP数据采集，并存储采集到的的目标井的零偏移距VSP地震数据，服务器可以从预先存储的实际勘探资料中获取该目标井的三分量检波器检测到的零偏移距VSP地震数据。
在S104中，基于预设下行波去除算法，对零偏移距VSP地震数据进行处理，得到上行波在预设三维正交坐标系的水平方向上的R分量和T分量。
其中，预设下行波去除算法可以是任意能够在保证地震波矢量特征不变、保持三分量相对振幅的前提下，去除下行波的算法，如图3所示，预设三维正交坐标系可以是由三分量检波器检测到的垂直方向上的Z分量，以及水平方向上的X分量、Y分量构成，上行波在预设三维正交坐标系的水平方向上的R分量和T分量可以是由上行波在三维正交坐标系的水平方向上X分量、Y分量和Z分量确定的。
在实施中，例如，如图4所示，服务器可以在Z分量上拾取纵波初至的时窗，并应用于X分量和Y分量，得到X、Y分量上的纵波初至。
在按照X分量上纵波初至能量最大，Y分量上纵波初至能量最小的原则，服务器可以对X分量和Y分量进行水平旋转，以得到水平方向上的R分量和T分量，此时，R分量可以在炮检平面内，T分量可以在水平面内与R分量正交。
利用旋转后的R分量的纵波初至和Z分量的纵波初至，在炮检平面内进行旋转。服务器可以将R分量的纵波旋转到Z分量上，将Z分量上的横波旋转到R分量上，恢复纵波和横波的波形泄露，此时，反射波矢量只投影在R分量上。
服务器可以将处理后的R分量、T分量和Z分量进行三分量一体化振幅补偿，在补偿振幅衰减的同时，保持各分量振幅的相对关系，服务器在利用线性Radon变换约束的中值滤波进行上下行波分离，以得到上行波在水平方向上的R分量和T分量。
上述下行波去除算法是一种可选地、可实现的去除算法，在实际应用场景中，还可以有多种不同的下行波去除算法，可以根据实际应用场景的不同而有所不同，本发明实施例对此不做具体限定。
在S106中，基于预设波场分离方法，对R分量和T分量进行波矢分离处理，得到R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波。
其中，预设波场分离方法可以是任意能够将上行波(包括快横波和慢横波)的波矢从上行波在水平方向上的R分量和T分量上分离出来的方法。
在实施中，由于快横波(PS1波)与慢横波(PS2波)的波速不同，所以，如图5(a)所示，在R分量上的PS1(R)与PS2(R)波会呈现出不同的视速度，如图5(b)所示，在T分量上的PS1(T)与PS2(T)波也会呈现出不同的视速度。
因此，可以通过预设波场分离方法(如基于视速度约束的f-k滤波方法)对R分量和T分量进行波矢分离处理，得到R分量上的快横波(即PS1(R))、T分量上的快横波(即PS1(T))、R分量上的慢横波(PS2(R))以及T分量上的慢横波(即PS2(T))。
在S108中，基于R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，确定目标夹角。
其中，目标夹角可以为由快横波在预设三维正交坐标系的水平方向上的方位和慢横波在预设三维正交坐标系的水平方向上的方位确定的夹角。
在实施中，如图6所示，服务器可以通过R分量上的快横波和T分量上的快横波，可以确定与快横波对应的夹角，通过R分量上的慢横波和T分量上的慢横波，可以确定与慢横波对应的夹角。服务器可以基于与快横波对应的夹角，以及与慢横波对应的夹角，确定目标夹角。
在S110中，基于目标夹角、三分量检波器在预设三维正交坐标系的水平方向上的R分量方位，确定目标裂缝的方位信息。
在实施中，服务器可以将目标夹角与三分量检波器R分量方位的和，确定为目标裂缝的方位角，基于目标裂缝的方位角，可以确定目标裂缝的方位信息。
本发明实施例提供一种基于零偏VSP数据确定裂缝方位的方法，本发明实施例通过获取目标井的三分量检波器检测到的零偏移距VSP地震数据，基于预设下行波去除算法，对零偏移距VSP地震数据进行处理，得到上行波在预设三维正交坐标系的水平方向上的R分量和T分量，基于预设波场分离方法，对R分量和T分量进行波矢分离处理，得到R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，基于R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，确定目标夹角，目标夹角为由快横波在预设三维正交坐标系的水平方向上的方位和慢横波在预设三维正交坐标系的水平方向上的方位确定的夹角，基于目标夹角、三分量检波器在预设三维正交坐标系的水平方向上的R分量方位，确定目标裂缝的方位信息。这样，由于通过预设下行波去除算法，对零偏移距VSP地震数据进行了处理，所以，可以避免上覆地层的影响以及层匹配的需求，在通过波场分离方法，将视速度不同的快横波和慢横波分离为R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，这样，基于分离后的数据，就可以准确的确定目标裂缝的方位信息，提高目标裂缝的方位信息的确定准确性。
实施例二
如图7所示，本发明实施例提供一种基于零偏VSP数据确定裂缝方位的方法，该方法的执行主体可以为服务器，该服务器可以是独立的服务器，也可以是由多个服务器组成的服务器集群。该方法具体可以包括以下步骤：
在S702中，获取目标井的三分量检波器检测到的零偏移距VSP地震数据。
在S704中，基于预设下行波去除算法，对零偏移距VSP地震数据进行处理，得到上行波在预设三维正交坐标系的水平方向上的R分量和T分量。
在S706中，基于预设波场分离方法，对R分量和T分量进行波矢分离处理，得到R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波。
上述S702～S706的具体处理过程可以参见上述实施例一中S102～S106的相关内容，在此不再赘述。
在S708中，基于零偏移距VSP地震数据，确定纵波初至位置。
其中，纵波初至位置可以是由三分量检波器最先接收到的纵波初至波确定的位置。
在实施中，如图5所示，服务器可以基于零偏移距VSP地震数据，分别在R分量剖面和T分量剖面上确定由纵波初至波确定的纵波初至位置。
在S710中，在R分量剖面和T分量剖面上，确定快横波和慢横波的目标相同起始点。
其中，如图5所示，目标相同起始点可以是由快横波和慢横波发生分裂的瞬时位置，在R分量剖面和T分量剖面上确定的起始点。
在S712中，在R分量剖面上，基于纵波初至位置以及目标相同起始点，确定第一数据范围，并在T分量剖面上，基于纵波初至位置以及目标相同起始点，确定第二数据范围。
在实施中，如图5(a)所示，第一数据范围可以是在R分量剖面上，由目标相同起始点对应的三分量检波器检测到的数据与纵波初至位置确定的数据范围。
如图5(b)所示，第二数据范围可以是在T分量剖面上，由目标相同起始点对应的三分量检波器检测到的数据与纵波初至位置确定的数据范围。
在S714中，基于第一数据范围内的R分量上的快横波和R分量上的慢横波，以及第二数据范围内的T分量上的快横波和T分量上的快横波，确定目标夹角。
在实施中，在实际应用中，上述S714的处理方式可以多种多样，以下提供一种可选的实现方式，具体可以参见下述步骤一～步骤五处理。
步骤一，基于预设第一滑动时窗，在R分量剖面上，对第一数据范围内的R分量上的快横波进行扫描处理，得到第一目标快横波。
在实施中，如图8(a)所示，预设第一滑动时窗可以为平行四边形滑动时窗，预设第一滑动时窗在R分量剖面上的斜边的斜率可以与R分量上的快横波的同相轴的斜率相同。
第一滑动时窗可以在第一数据范围内，从R分量上的快横波的初始位置滑动至目标相同起始点与纵波初至位置相交的位置。服务器可以将每个第一滑动时窗在R分量剖面内包含的R分量上的快横波波矢振幅的均值，作为每个第一滑动时窗对应的第一快横波波矢在R分量上的投影，其中，每两个相邻第一滑动时窗包含预设数量的相同的检波器。
服务器可以基于每个第一滑动时窗对应的第一快横波波矢在R分量上的投影，确定第一目标快横波。例如，服务器可以将每个第一滑动时窗对应的第一快横波波矢振幅的均值，作为第一目标快横波。
步骤二，基于预设第一滑动时窗，在R分量剖面上，对第一数据范围内的R分量上的慢横波进行扫描处理，得到第一目标慢横波。
在实施中，如图8(c)所示，服务器也可以在R分量剖面上，对第一数据范围内的R分量上的慢横波进行扫描处理，得到每个第一滑动时窗对应的第一慢横波，在基于每个第一滑动时窗对应的第一慢横波，确定第一目标慢横波。
步骤三，基于预设第二滑动时窗，在T分量剖面上，对第二数据范围内的T分量上的快横波进行扫描处理，得到第二目标快横波。
步骤四，基于预设第二滑动时窗，在T分量剖面上，对第二数据范围内的T分量上的慢横波进行扫描处理，得到第二目标慢横波。
其中，如图8(b)和图8(d)所示，预设第二滑动时窗可以为平行四边形滑动时窗，预设第二滑动时窗在R分量剖面上的斜边的斜率可以与R分量上的慢横波的同相轴的斜率相同。
步骤五，基于第一目标快横波、第一目标慢横波、第二目标快横波、第二目标慢横波，确定目标夹角。
在实施中，服务器可以将第一目标快横波和第二目标快横波代入公式
Figure BDA0002806686830000131
得到第一夹角，其中，θ1为第一夹角，PS1(R)为第一目标快横波，PS1(T)为第二目标快横波。
服务器可以将第一目标慢横波和第二目标慢横波代入公式
Figure BDA0002806686830000132
得到第二夹角，其中，θ2为第二夹角，PS2(R)为第一目标慢横波，PS2(T)为第二目标慢横波。
由于用于确定第一夹角所需要的R分量上的快横波、T分量上的快横波，与用于确定第二夹角所需要的R分量上的慢横波、T分量上的慢横波的视速度不同，因此可以采用两个斜率不同的平行四边形滑动时窗(即第一滑动时窗和第二滑动时窗)。这两个平行四边形滑动时窗在水平方向跨越的三分量检波器的数量相同，在竖直方向上的采样点数也相同。
服务器可以将第一夹角和第二夹角代入公式
θ＝θ1*f+θ2*(1-f)
得到目标夹角，其中，θ为目标夹角，θ1为第一夹角，θ2为第二夹角，f为预设加权系数，f可以随着快横波较慢横波的信噪比的提高而提高，且f的取值范围可以为预设范围(如0–1之间)。
在S716中，将目标夹角和三分量检波器在预设三维正交坐标系的水平方向上的R分量方位的和，作为目标裂缝的方位角。
在S718中，基于纵波初至位置，确定纵波单程旅行时。
在S720中，基于目标裂缝的方位角和纵波单程旅行时，确定目标裂缝的方位信息。
在实施中，目前，裂缝的方位信息的预测方法大多是使用动校正后的R分量和T分量，即将地震记录校正为图射线路径的自激自收的模式。但是，由于距离目标井较近的三分量检波器(如图2中的1号三分量检波器)检测到的快横波和慢横波之间存在较大时差，难以进行矢量计算，且PS波的射线路径经过多个地层，由此确定的目标裂缝的方位信息会受到上覆地层的影响。
因此，可以使用动校正前的数据，来确定目标裂缝的方位信息。4个波(即R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波)在纵波初至位置处具有相同的起始点(即目标相同起始点)，其传播模式可以如图2中的射线路径②所示，下行P波在界面处反射后形成的上行PS波在经过极小的传播路径后就可以由204号检波器接收。此时，PS波虽然发生了分裂，但因为传播路径极小而无法在地震剖面处显示出时差。因此，可以沿纵波初至位置对上述4个波进行处理，以得到目标裂缝的方位信息。
例如，服务器可以获取目标井的三分量检波器检测到的零偏移距VSP地震数据，并对零偏移距VSP地震数据进行处理，得到如图9所示的上行波在预设三维正交坐标系的水平方向上的R分量和T分量，图9中标示的R分量上的一组分裂横波和T分量上的一组分裂横波均在纵波初至位置处没有时差，但是随着检波号由大到小、上行分裂横波旅行时的增加，分裂横波之间的时差也随之增加。
因此，服务器可以使用波场分离方法，对R分量和T分量进行波矢分离处理，得到R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波。服务器可以利用平行四边形滑动时窗(即第一滑动时窗和第二滑动时窗)分别确定第一夹角和第二夹角，再由第一夹角和第二夹角确定目标夹角，并基于纵波初至位置，确定纵波单程旅行时，最后，由目标裂缝的方位角和纵波单程旅行时，可以确定如图10所示的目标裂缝的方位信息。该段目标裂缝的方位角主要分布在50–80°之间。
本发明实施例提供一种基于零偏VSP数据确定裂缝方位的方法，本发明实施例通过获取目标井的三分量检波器检测到的零偏移距VSP地震数据，基于预设下行波去除算法，对零偏移距VSP地震数据进行处理，得到上行波在预设三维正交坐标系的水平方向上的R分量和T分量，基于预设波场分离方法，对R分量和T分量进行波矢分离处理，得到R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，基于R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，确定目标夹角，目标夹角为由快横波在预设三维正交坐标系的水平方向上的方位和慢横波在预设三维正交坐标系的水平方向上的方位确定的夹角，基于目标夹角、三分量检波器在预设三维正交坐标系的水平方向上的R分量方位，确定目标裂缝的方位信息。这样，由于通过预设下行波去除算法，对零偏移距VSP地震数据进行了处理，所以，可以避免上覆地层的影响以及层匹配的需求，在通过波场分离方法，将视速度不同的快横波和慢横波分离为R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，这样，基于分离后的数据，就可以准确的确定目标裂缝的方位信息，提高目标裂缝的方位信息的确定准确性。
实施例三
以上为本说明书实施例提供的基于零偏VSP数据确定裂缝方位的方法，基于同样的思路，本说明书实施例还提供一种基于零偏VSP数据确定裂缝方位的装置，如图11所示。
该基于零偏VSP数据确定裂缝方位的装置包括：数据获取模块1101、数据处理模块1102、波场分离模块1103、夹角确定模块1104和方位确定模块1105，其中：
数据获取模块1101，用于获取目标井的三分量检波器检测到的零偏移距VSP地震数据；
数据处理模块1102，用于基于预设下行波去除算法，对所述零偏移距VSP地震数据进行处理，得到上行波在预设三维正交坐标系的水平方向上的R分量和T分量；
波场分离模块1103，用于基于预设波场分离方法，对所述R分量和T分量进行波矢分离处理，得到R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波；
夹角确定模块1104，用于基于所述R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，确定目标夹角，所述目标夹角为由所述快横波在所述预设三维正交坐标系的水平方向上的方位和所述慢横波在所述预设三维正交坐标系的水平方向上的方位确定的夹角；
方位确定模块1105，用于基于所述目标夹角、所述三分量检波器在所述预设三维正交坐标系的水平方向上的R分量方位，确定目标裂缝的方位信息。
本发明实施例中，所述夹角确定模块1104，用于：
基于所述零偏移距VSP地震数据，确定纵波初至位置；
在所述R分量剖面和所述T分量剖面上，确定所述快横波和所述慢横波的目标相同起始点；
在所述R分量剖面上，基于所述纵波初至位置以及所述目标相同起始点，确定第一数据范围，并在所述T分量剖面上，基于所述纵波初至位置以及所述目标相同起始点，确定第二数据范围；
基于所述第一数据范围内的所述R分量上的快横波和R分量上的慢横波，以及所述第二数据范围内的所述T分量上的快横波和T分量上的快横波，确定所述目标夹角。
本发明实施例中，所述夹角确定模块1104，用于：
基于预设第一滑动时窗，在所述R分量剖面上，对所述第一数据范围内的所述R分量上的快横波进行扫描处理，得到第一目标快横波；
基于所述预设第一滑动时窗，在所述R分量剖面上，对所述第一数据范围内的所述R分量上的慢横波进行扫描处理，得到第一目标慢横波；
基于预设第二滑动时窗，在所述T分量剖面上，对所述第二数据范围内的所述T分量上的快横波进行扫描处理，得到第二目标快横波；
基于所述预设第二滑动时窗，在所述T分量剖面上，对所述第二数据范围内的所述T分量上的慢横波进行扫描处理，得到第二目标慢横波；
基于所述第一目标快横波、所述第一目标慢横波、所述第二目标快横波、所述第二目标慢横波，确定所述目标夹角。
本发明实施例中，所述夹角确定模块1104，用于：
将每个所述第一滑动时窗在所述R分量剖面内包含的R分量上的快横波波矢振幅的均值，作为每个所述第一滑动时窗对应的第一快横波波矢在R分量上的投影，其中，每两个相邻所述第一滑动时窗包含预设数量的相同的检波器；
基于每个所述第一滑动时窗对应的第一快横波波矢在R分量上的投影，确定所述第一目标快横波。
本发明实施例中，所述预设第一滑动时窗和所述第二滑动时窗为平行四边形滑动时窗，所述预设第一滑动时窗在所述R分量剖面上的斜边的斜率与所述R分量上的快横波的同相轴的斜率相同，所述预设第二滑动时窗在所述R分量剖面上的斜边的斜率与所述R分量上的慢横波的同相轴的斜率相同。
本发明实施例中，所述夹角确定模块1104，用于：
将所述第一目标快横波和所述第二目标快横波代入公式
Figure BDA0002806686830000161
得到第一夹角，其中，θ1为所述第一夹角，PS1(R)为所述第一目标快横波，PS1(T)为所述第二目标快横波；
将所述第一目标慢横波和所述第二目标慢横波代入公式
Figure BDA0002806686830000171
得到第二夹角，其中，θ2为所述第二夹角，PS2(R)为所述第一目标慢横波，PS2(T)为所述第二目标慢横波；
将所述第一夹角和所述第二夹角代入公式
θ＝θ1*f+θ2*(1-f)
得到所述目标夹角，其中，θ为所述目标夹角，θ1为所述第一夹角，θ2为所述第二夹角，f为预设加权系数。
本发明实施例中，所述方位确定模块1105，用于：
将所述目标夹角和所述三分量检波器在所述预设三维正交坐标系的水平方向上的R分量方位的和，作为所述目标裂缝的方位角；
基于所述纵波初至位置，确定纵波单程旅行时；
基于所述目标裂缝的方位角和所述纵波单程旅行时，确定所述目标裂缝的方位信息。
本发明实施例提供一种基于零偏VSP数据确定裂缝方位的方法，本发明实施例通过获取目标井的三分量检波器检测到的零偏移距VSP地震数据，基于预设下行波去除算法，对零偏移距VSP地震数据进行处理，得到上行波在预设三维正交坐标系的水平方向上的R分量和T分量，基于预设波场分离方法，对R分量和T分量进行波矢分离处理，得到R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，基于R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，确定目标夹角，目标夹角为由快横波在预设三维正交坐标系的水平方向上的方位和慢横波在所述预设三维正交坐标系的水平方向上的方位确定的夹角，基于目标夹角、三分量检波器在预设三维正交坐标系的水平方向上的R分量方位，确定目标裂缝的方位信息。这样，由于通过预设下行波去除算法，对零偏移距VSP地震数据进行了处理，所以，可以避免上覆地层的影响以及层匹配的需求，在通过波场分离方法，将视速度不同的快横波和慢横波分离为R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，这样，基于分离后的数据，就可以准确的确定目标裂缝的方位信息，提高目标裂缝的方位信息的确定准确性。
实施例四
图11为实现本发明各个实施例的一种电子设备的硬件结构示意图，
该电子设备1200包括但不限于：射频单元1201、网络模块1202、音频输出单元1203、输入单元1204、传感器1205、显示单元1206、用户输入单元1204、接口单元1208、存储器1209、处理器1210、以及电源1211等部件。本领域技术人员可以理解，图11中示出的电子设备结构并不构成对电子设备的限定，电子设备可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。在本发明实施例中，电子设备包括但不限于手机、平板电脑、笔记本电脑、掌上电脑、车载终端、可穿戴设备、以及计步器等。
其中，处理器1210，用于：获取目标井的三分量检波器检测到的零偏移距VSP地震数据；基于预设下行波去除算法，对所述零偏移距VSP地震数据进行处理，得到上行波在预设三维正交坐标系的水平方向上的R分量和T分量；基于预设波场分离方法，对所述R分量和T分量进行波矢分离处理，得到R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波；基于所述R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，确定目标夹角，所述目标夹角为由所述快横波在所述预设三维正交坐标系的水平方向上的方位和所述慢横波在所述预设三维正交坐标系的水平方向上的方位确定的夹角；基于所述目标夹角、所述三分量检波器在所述预设三维正交坐标系的水平方向上的R分量方位，确定目标裂缝的方位信息。
所述处理器1210，还用于：基于所述零偏移距VSP地震数据，确定纵波初至位置；在所述R分量剖面和所述T分量剖面上，确定所述快横波和所述慢横波的目标相同起始点；在所述R分量剖面上，基于所述纵波初至位置以及所述目标相同起始点，确定第一数据范围，并在所述T分量剖面上，基于所述纵波初至位置以及所述目标相同起始点，确定第二数据范围；基于所述第一数据范围内的所述R分量上的快横波和R分量上的慢横波，以及所述第二数据范围内的所述T分量上的快横波和T分量上的快横波，确定所述目标夹角。
所述处理器1210，还用于：基于预设第一滑动时窗，在所述R分量剖面上，对所述第一数据范围内的所述R分量上的快横波进行扫描处理，得到第一目标快横波；基于所述预设第一滑动时窗，在所述R分量剖面上，对所述第一数据范围内的所述R分量上的慢横波进行扫描处理，得到第一目标慢横波；基于预设第二滑动时窗，在所述T分量剖面上，对所述第二数据范围内的所述T分量上的快横波进行扫描处理，得到第二目标快横波；基于所述预设第二滑动时窗，在所述T分量剖面上，对所述第二数据范围内的所述T分量上的慢横波进行扫描处理，得到第二目标慢横波；基于所述第一目标快横波、所述第一目标慢横波、所述第二目标快横波、所述第二目标慢横波，确定所述目标夹角。
所述处理器1210，还用于：将每个所述第一滑动时窗在所述R分量剖面内包含的R分量上的快横波波矢振幅的均值，作为每个所述第一滑动时窗对应的第一快横波波矢在R分量上的投影，其中，每两个相邻所述第一滑动时窗包含预设数量的相同的检波器；基于每个所述第一滑动时窗对应的第一快横波波矢在R分量上的投影，确定所述第一目标快横波。
此外，所述预设第一滑动时窗和所述第二滑动时窗为平行四边形滑动时窗，所述预设第一滑动时窗在所述R分量剖面上的斜边的斜率与所述R分量上的快横波的同相轴的斜率相同，所述预设第二滑动时窗在所述R分量剖面上的斜边的斜率与所述R分量上的慢横波的同相轴的斜率相同。
所述处理器1210，还用于：将所述第一目标快横波和所述第二目标快横波代入公式
Figure BDA0002806686830000191
得到第一夹角，其中，θ1为所述第一夹角，PS1(R)为所述第一目标快横波，PS1(T)为所述第二目标快横波；将所述第一目标慢横波和所述第二目标慢横波代入公式
Figure BDA0002806686830000192
得到第二夹角，其中，θ2为所述第二夹角，PS2(R)为所述第一目标慢横波，PS2(T)为所述第二目标慢横波；将所述第一夹角和所述第二夹角代入公式θ＝θ1*f+θ2*(1-f)得到所述目标夹角，其中，θ为所述目标夹角，θ1为所述第一夹角，θ2为所述第二夹角，f为预设加权系数。
所述处理器1210，还用于：将所述目标夹角和所述三分量检波器在所述预设三维正交坐标系的水平方向上的R分量方位的和，作为所述目标裂缝的方位角；基于所述纵波初至位置，确定纵波单程旅行时；基于所述目标裂缝的方位角和所述纵波单程旅行时，确定所述目标裂缝的方位信息。
本发明实施例提供一种电子设备，本发明实施例通过获取目标井的三分量检波器检测到的零偏移距VSP地震数据，基于预设下行波去除算法，对零偏移距VSP地震数据进行处理，得到上行波在预设三维正交坐标系的水平方向上的R分量和T分量，基于预设波场分离方法，对R分量和T分量进行波矢分离处理，得到R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，基于R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，确定目标夹角，目标夹角为由快横波在预设三维正交坐标系的水平方向上的方位和慢横波在所述预设三维正交坐标系的水平方向上的方位确定的夹角，基于目标夹角、三分量检波器在预设三维正交坐标系的水平方向上的R分量方位，确定目标裂缝的方位信息。这样，由于通过预设下行波去除算法，对零偏移距VSP地震数据进行了处理，所以，可以避免上覆地层的影响以及层匹配的需求，在通过波场分离方法，将视速度不同的快横波和慢横波分离为R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，这样，基于分离后的数据，就可以准确的确定目标裂缝的方位信息，提高目标裂缝的方位信息的确定准确性。
应理解的是，本发明实施例中，射频单元1201可用于收发信息或通话过程中，信号的接收和发送，具体的，将来自基站的下行数据接收后，给处理器1210处理；另外，将上行的数据发送给基站。通常，射频单元1201包括但不限于天线、至少一个放大器、收发信机、耦合器、低噪声放大器、双工器等。此外，射频单元1201还可以通过无线通信系统与网络和其他设备通信。
电子设备通过网络模块1202为用户提供了无线的宽带互联网访问，如帮助用户收发电子邮件、浏览网页和访问流式媒体等。
音频输出单元1203可以将射频单元1201或网络模块1202接收的或者在存储器1209中存储的音频数据转换成音频信号并且输出为声音。而且，音频输出单元1203还可以提供与电子设备1200执行的特定功能相关的音频输出(例如，呼叫信号接收声音、消息接收声音等等)。音频输出单元1203包括扬声器、蜂鸣器以及受话器等。
输入单元1204用于接收音频或视频信号。输入单元1204可以包括图形处理器(Graphics Processing Unit，GPU)12041和麦克风12042，图形处理器12041对在视频捕获模式或图像捕获模式中由图像捕获装置(如摄像头)获得的静态图片或视频的图像数据进行处理。处理后的图像帧可以显示在显示单元1206上。经图形处理器12041处理后的图像帧可以存储在存储器1209(或其它存储介质)中或者经由射频单元1201或网络模块1202进行发送。麦克风12042可以接收声音，并且能够将这样的声音处理为音频数据。处理后的音频数据可以在电话通话模式的情况下转换为可经由射频单元1201发送到移动通信基站的格式输出。
电子设备1200还包括至少一种传感器1205，比如光传感器、运动传感器以及其他传感器。具体地，光传感器包括环境光传感器及接近传感器，其中，环境光传感器可根据环境光线的明暗来调节显示面板12061的亮度，接近传感器可在电子设备1200移动到耳边时，关闭显示面板12061和/或背光。作为运动传感器的一种，加速计传感器可检测各个方向上(一般为三轴)加速度的大小，静止时可检测出重力的大小及方向，可用于识别电子设备姿态(比如横竖屏切换、相关游戏、磁力计姿态校准)、振动识别相关功能(比如计步器、敲击)等；传感器1205还可以包括指纹传感器、压力传感器、虹膜传感器、分子传感器、陀螺仪、气压计、湿度计、温度计、红外线传感器等，在此不再赘述。
显示单元1206用于显示由用户输入的信息或提供给用户的信息。显示单元1206可包括显示面板12061，可以采用液晶显示器(Liquid Crystal Display，LCD)、有机发光二极管(Organic Light-Emitting Diode,OLED)等形式来配置显示面板12061。
用户输入单元1204可用于接收输入的数字或字符信息，以及产生与电子设备的用户设置以及功能控制有关的键信号输入。具体地，用户输入单元1204包括触控面板12071以及其他输入设备12072。触控面板12071，也称为触摸屏，可收集用户在其上或附近的触摸操作(比如用户使用手指、触笔等任何适合的物体或附件在触控面板12071上或在触控面板12071附近的操作)。触控面板12071可包括触摸检测装置和触摸控制器两个部分。其中，触摸检测装置检测用户的触摸方位，并检测触摸操作带来的信号，将信号传送给触摸控制器；触摸控制器从触摸检测装置上接收触摸信息，并将它转换成触点坐标，再送给处理器1210，接收处理器1210发来的命令并加以执行。此外，可以采用电阻式、电容式、红外线以及表面声波等多种类型实现触控面板12071。除了触控面板12071，用户输入单元1204还可以包括其他输入设备12072。具体地，其他输入设备12072可以包括但不限于物理键盘、功能键(比如音量控制按键、开关按键等)、轨迹球、鼠标、操作杆，在此不再赘述。
进一步的，触控面板12071可覆盖在显示面板12061上，当触控面板12071检测到在其上或附近的触摸操作后，传送给处理器1210以确定触摸事件的类型，随后处理器1210根据触摸事件的类型在显示面板12061上提供相应的视觉输出。虽然在图11中，触控面板12071与显示面板12061是作为两个独立的部件来实现电子设备的输入和输出功能，但是在某些实施例中，可以将触控面板12071与显示面板12061集成而实现电子设备的输入和输出功能，具体此处不做限定。
接口单元1208为外部装置与电子设备1200连接的接口。例如，外部装置可以包括有线或无线头戴式耳机端口、外部电源(或电池充电器)端口、有线或无线数据端口、存储卡端口、用于连接具有识别模块的装置的端口、音频输入/输出(I/O)端口、视频I/O端口、耳机端口等等。接口单元1208可以用于接收来自外部装置的输入(例如，数据信息、电力等等)并且将接收到的输入传输到电子设备1200内的一个或多个元件或者可以用于在电子设备1200和外部装置之间传输数据。
存储器1209可用于存储软件程序以及各种数据。存储器1209可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序(比如声音播放功能、图像播放功能等)等；存储数据区可存储根据手机的使用所创建的数据(比如音频数据、电话本等)等。此外，存储器409可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。
处理器1210是电子设备的控制中心，利用各种接口和线路连接整个电子设备的各个部分，通过运行或执行存储在存储器1209内的软件程序和/或模块，以及调用存储在存储器1209内的数据，执行电子设备的各种功能和处理数据，从而对电子设备进行整体监控。处理器1210可包括一个或多个处理单元；优选的，处理器1210可集成应用处理器和调制解调处理器，其中，应用处理器主要处理操作系统、用户界面和应用程序等，调制解调处理器主要处理无线通信。可以理解的是，上述调制解调处理器也可以不集成到处理器1210中。
电子设备1200还可以包括给各个部件供电的电源1211(比如电池)，优选的，电源1211可以通过电源管理系统与处理器1210逻辑相连，从而通过电源管理系统实现管理充电、放电、以及功耗管理等功能。
优选的，本发明实施例还提供一种电子设备，包括处理器1210，存储器1209，存储在存储器1209上并可在所述处理器1210上运行的计算机程序，该计算机程序被处理器1210执行时实现上述基于零偏VSP数据确定裂缝方位的方法实施例的各个过程，且能达到相同的技术效果，为避免重复，这里不再赘述。
实施例五
本发明实施例还提供一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述基于零偏VSP数据确定裂缝方位的方法实施例的各个过程，且能达到相同的技术效果，为避免重复，这里不再赘述。其中，所述的计算机可读存储介质，如只读存储器(Read-Only Memory，简称ROM)、随机存取存储器(RandomAccess Memory，简称RAM)、磁碟或者光盘等。
本发明实施例提供一种计算机可读存储介质，通过获取目标井的三分量检波器检测到的零偏移距VSP地震数据，基于预设下行波去除算法，对零偏移距VSP地震数据进行处理，得到上行波在预设三维正交坐标系的水平方向上的R分量和T分量，基于预设波场分离方法，对R分量和T分量进行波矢分离处理，得到R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，基于R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，确定目标夹角，目标夹角为由快横波在预设三维正交坐标系的水平方向上的方位和慢横波在所述预设三维正交坐标系的水平方向上的方位确定的夹角，基于目标夹角、三分量检波器在预设三维正交坐标系的水平方向上的R分量方位，确定目标裂缝的方位信息。这样，由于通过预设下行波去除算法，对零偏移距VSP地震数据进行了处理，所以，可以避免上覆地层的影响以及层匹配的需求，在通过波场分离方法，将视速度不同的快横波和慢横波分离为R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，这样，基于分离后的数据，就可以准确的确定目标裂缝的方位信息，提高目标裂缝的方位信息的确定准确性。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本发明的实施例可提供为方法、系统或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
以上所述仅为本发明的实施例而已，并不用于限制本发明。对于本领域技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。

Claims (10)

1.一种基于零偏VSP数据确定裂缝方位的方法，包括：
获取目标井的三分量检波器检测到的零偏移距VSP地震数据；
基于预设下行波去除算法，对所述零偏移距VSP地震数据进行处理，得到上行波在预设三维正交坐标系的水平方向上的R分量和T分量；其中，所述预设下行波去除算法是能够在保证地震波矢量特征不变、保持三分量相对振幅的前提下，去除下行波的算法；
基于预设波场分离方法，对所述R分量和T分量进行波矢分离处理，得到R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波；
基于所述R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，确定目标夹角，所述目标夹角为由所述快横波在所述预设三维正交坐标系的水平方向上的方位和所述慢横波在所述预设三维正交坐标系的水平方向上的方位确定的夹角；
基于所述目标夹角、所述三分量检波器在所述预设三维正交坐标系的水平方向上的R分量方位，确定目标裂缝的方位信息。
2.根据权利要求1所述的方法，其特征在于，所述基于所述R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，确定目标夹角，包括：
基于所述零偏移距VSP地震数据，确定纵波初至位置；
在所述R分量剖面和所述T分量剖面上，确定所述快横波和所述慢横波的目标相同起始点；
在所述R分量剖面上，基于所述纵波初至位置以及所述目标相同起始点，确定第一数据范围，并在所述T分量剖面上，基于所述纵波初至位置以及所述目标相同起始点，确定第二数据范围；
基于所述第一数据范围内的所述R分量上的快横波和R分量上的慢横波，以及所述第二数据范围内的所述T分量上的快横波和T分量上的快横波，确定所述目标夹角。
3.根据权利要求2所述的方法，其特征在于，所述基于所述第一数据范围内的所述R分量上的快横波和R分量上的慢横波，以及所述第二数据范围内的所述T分量上的快横波和T分量上的快横波，确定所述目标夹角，包括：
基于预设第一滑动时窗，在所述R分量剖面上，对所述第一数据范围内的所述R分量上的快横波进行扫描处理，得到第一目标快横波；
基于所述预设第一滑动时窗，在所述R分量剖面上，对所述第一数据范围内的所述R分量上的慢横波进行扫描处理，得到第一目标慢横波；
基于预设第二滑动时窗，在所述T分量剖面上，对所述第二数据范围内的所述T分量上的快横波进行扫描处理，得到第二目标快横波；
基于所述预设第二滑动时窗，在所述T分量剖面上，对所述第二数据范围内的所述T分量上的慢横波进行扫描处理，得到第二目标慢横波；
基于所述第一目标快横波、所述第一目标慢横波、所述第二目标快横波、所述第二目标慢横波，确定所述目标夹角。
4.根据权利要求3所述的方法，其特征在于，所述基于预设第一滑动时窗，在所述R分量剖面上，对所述第一数据范围内的所述R分量上的快横波进行扫描处理，得到第一目标快横波，包括：
将每个所述第一滑动时窗在所述R分量剖面内包含的R分量上的快横波波矢振幅的均值，作为每个所述第一滑动时窗对应的第一快横波波矢在R分量上的投影，其中，每两个相邻所述第一滑动时窗包含预设数量的相同的检波器；
基于每个所述第一滑动时窗对应的第一快横波波矢在R分量上的投影，确定所述第一目标快横波。
5.根据权利要求3所述的方法，其特征在于，所述预设第一滑动时窗和所述第二滑动时窗为平行四边形滑动时窗，所述预设第一滑动时窗在所述R分量剖面上的斜边的斜率与所述R分量上的快横波的同相轴的斜率相同，所述预设第二滑动时窗在所述R分量剖面上的斜边的斜率与所述R分量上的慢横波的同相轴的斜率相同。
6.根据权利要求4所述的方法，其特征在于，所述基于所述第一目标快横波、所述第一目标慢横波、所述第二目标快横波、所述第二目标慢横波，确定所述目标夹角，包括：
将所述第一目标快横波和所述第二目标快横波代入公式
Figure FDA0003056888450000031
得到第一夹角，其中，θ1为所述第一夹角，PS1(R)为所述第一目标快横波，PS1(T)为所述第二目标快横波；
将所述第一目标慢横波和所述第二目标慢横波代入公式
Figure FDA0003056888450000032
得到第二夹角，其中，θ2为所述第二夹角，PS2(R)为所述第一目标慢横波，PS2(T)为所述第二目标慢横波；
将所述第一夹角和所述第二夹角代入公式
θ＝θ1*f+θ2*(1-f)
得到所述目标夹角，其中，θ为所述目标夹角，θ1为所述第一夹角，θ2为所述第二夹角，f为预设加权系数。
7.根据权利要求6所述的方法，其特征在于，所述基于所述目标夹角、所述三分量检波器在所述预设三维正交坐标系的水平方向上的R分量方位，确定目标裂缝的方位信息，包括：
将所述目标夹角和所述三分量检波器在所述预设三维正交坐标系的水平方向上的R分量方位的和，作为所述目标裂缝的方位角；
基于所述纵波初至位置，确定纵波单程旅行时；
基于所述目标裂缝的方位角和所述纵波单程旅行时，确定所述目标裂缝的方位信息。
8.一种基于零偏VSP数据确定裂缝方位的装置，包括：
数据获取模块，用于获取目标井的三分量检波器检测到的零偏移距VSP地震数据；
数据处理模块，用于基于预设下行波去除算法，对所述零偏移距VSP地震数据进行处理，得到上行波在预设三维正交坐标系的水平方向上的R分量和T分量；其中，所述预设下行波去除算法是能够在保证地震波矢量特征不变、保持三分量相对振幅的前提下，去除下行波的算法；
波场分离模块，用于基于预设波场分离方法，对所述R分量和T分量进行波矢分离处理，得到R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波；
夹角确定模块，用于基于所述R分量上的快横波、T分量上的快横波、R分量上的慢横波以及T分量上的慢横波，确定目标夹角，所述目标夹角为由所述快横波在所述预设三维正交坐标系的水平方向上的方位和所述慢横波在所述预设三维正交坐标系的水平方向上的方位确定的夹角；
方位确定模块，用于基于所述目标夹角、所述三分量检波器在所述预设三维正交坐标系的水平方向上的R分量方位，确定目标裂缝的方位信息。
9.根据权利要求8所述的装置，其特征在于，所述夹角确定模块，用于：
基于所述零偏移距VSP地震数据，确定纵波初至位置；
在所述R分量剖面和所述T分量剖面上，确定所述快横波和所述慢横波的目标相同起始点；
在所述R分量剖面上，基于所述纵波初至位置以及所述目标相同起始点，确定第一数据范围，并在所述T分量剖面上，基于所述纵波初至位置以及所述目标相同起始点，确定第二数据范围；
基于所述第一数据范围内的所述R分量上的快横波和R分量上的慢横波，以及所述第二数据范围内的所述T分量上的快横波和T分量上的快横波，确定所述目标夹角。
10.一种电子设备，其特征在于，包括处理器、存储器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述计算机程序被所述处理器执行时实现如权利要求1至7中任一项所述的基于零偏VSP数据确定裂缝方位的方法的步骤。
CN202011371093.0A 2020-11-30 2020-11-30 基于零偏vsp数据确定裂缝方位的方法、装置及电子设备 Active CN112558144B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011371093.0A CN112558144B (zh) 2020-11-30 2020-11-30 基于零偏vsp数据确定裂缝方位的方法、装置及电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011371093.0A CN112558144B (zh) 2020-11-30 2020-11-30 基于零偏vsp数据确定裂缝方位的方法、装置及电子设备

Publications (2)

Publication Number Publication Date
CN112558144A CN112558144A (zh) 2021-03-26
CN112558144B true CN112558144B (zh) 2021-07-06

Family

ID=75046690

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011371093.0A Active CN112558144B (zh) 2020-11-30 2020-11-30 基于零偏vsp数据确定裂缝方位的方法、装置及电子设备

Country Status (1)

Country Link
CN (1) CN112558144B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115061194B (zh) * 2022-06-02 2023-04-28 中国地质大学（北京） 一种三分量vsp预测两组裂缝的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106468781A (zh) * 2015-08-19 2017-03-01 中国石油化工股份有限公司 一种基于最小熵旋转法的裂缝预测方法
CN107132575A (zh) * 2017-05-12 2017-09-05 中国地质大学（北京） 基于横波极化分析预测裂缝方位角的方法
CN109799529A (zh) * 2019-01-31 2019-05-24 河海大学 一种基于互相关的横波分裂vsp裂缝预测方法
US20200209418A1 (en) * 2018-12-27 2020-07-02 Saexploration, Inc. Method and Apparatus for Active Seismic Shear Wave Monitoring of Hydro-Fracturing of Oil and Gas Reservoirs Using Arrays of Multi-Component Sensors and Controlled Seismic Sources

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9784863B2 (en) * 2011-09-26 2017-10-10 Exxonmobil Upstream Research Company S-wave anisotropy estimate by automated image registration
US20140334262A1 (en) * 2013-04-22 2014-11-13 Robert H. Brune Method and Apparatus for Active Seismic Shear Wave Monitoring of Hydro-Fracturing of Oil and Gas Reservoirs Using Arrays of Multi-Component Sensors and Controlled Seismic Sources

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106468781A (zh) * 2015-08-19 2017-03-01 中国石油化工股份有限公司 一种基于最小熵旋转法的裂缝预测方法
CN107132575A (zh) * 2017-05-12 2017-09-05 中国地质大学（北京） 基于横波极化分析预测裂缝方位角的方法
US20200209418A1 (en) * 2018-12-27 2020-07-02 Saexploration, Inc. Method and Apparatus for Active Seismic Shear Wave Monitoring of Hydro-Fracturing of Oil and Gas Reservoirs Using Arrays of Multi-Component Sensors and Controlled Seismic Sources
CN109799529A (zh) * 2019-01-31 2019-05-24 河海大学 一种基于互相关的横波分裂vsp裂缝预测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Optimizing the design of vertical seismic profiling(VSP) for imaging fracture zones over hardrock basement geothermal environments;Fabienne Reiser 等;《Journal of Applied Geophysics》;20170221;第139卷;第25-35页 *
三分量VSP资料在裂缝检测中的应用;曹立斌 等;《石油物探》;20080131;第47卷(第1期);第67-71页 *

Also Published As

Publication number Publication date
CN112558144A (zh) 2021-03-26

Similar Documents

Publication Publication Date Title
CN107576986B (zh) 一种确定海底地震反射波旅行时的方法和装置
CN107505653B (zh) 一种确定叠前偏移时间结果的方法和装置
CN110174699B (zh) 确定断层封闭性的方法、装置及存储介质
CN110389380B (zh) 地震剖面的同相轴自动拾取方法、装置及存储介质
CN112558144B (zh) 基于零偏vsp数据确定裂缝方位的方法、装置及电子设备
CN110988993B (zh) 一种偏移成像方法、装置及电子设备
EP2693301A1 (en) Method and apparatus for identifying the posture of an object using ultrasonic waves
CN111427082A (zh) 地震响应模板的确定方法、装置及存储介质
CN111077570A (zh) 确定油藏中辉绿岩的侵入位置的方法、装置和存储介质
CN113009573B (zh) 地震数据采集方法、装置、计算机设备及存储介质
CN111897001B (zh) 波场确定方法与装置
CN114721046A (zh) 缝洞检测方法、装置、计算机存储介质
CN111505717B (zh) 一种三维vsp共面元成像方法与装置
US20140286125A1 (en) Seismic acquisition method and apparatus
CN111551987A (zh) 利用变偏移距vsp数据确定岩石弹性性质的方法与装置
CN113341456B (zh) 一种基于非规则网格的地震偏移方法、装置及电子设备
CN112698391B (zh) 砂体刻画方法、装置、设备及存储介质
CN117518269A (zh) 海洋勘探中压制多次波的方法、装置及计算机存储介质
CN110058306B (zh) 确定三维速度体的方法、装置及计算机存储介质
CN115508894B (zh) 地面网格参数的确定方法、装置、设备和存储介质
CN118151244A (zh) 对地震数据进行剩余静校正的方法和装置
CN112987105B (zh) 定量预测地下岩盐层分布的方法、装置、终端和存储介质
CN115542389A (zh) 地层温度的预测方法、装置及计算机存储介质
CN110630254B (zh) 油水层的区分方法、装置及计算机可读存储介质
CN115561813A (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
GR01 Patent grant
GR01 Patent grant