CN113031073A - 叠前深度偏移速度确定方法、装置及设备 - Google Patents
叠前深度偏移速度确定方法、装置及设备 Download PDFInfo
- Publication number
- CN113031073A CN113031073A CN201911347263.9A CN201911347263A CN113031073A CN 113031073 A CN113031073 A CN 113031073A CN 201911347263 A CN201911347263 A CN 201911347263A CN 113031073 A CN113031073 A CN 113031073A
- Authority
- CN
- China
- Prior art keywords
- sampling
- velocity
- longitudinal wave
- location
- speed
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 64
- 238000013508 migration Methods 0.000 title claims abstract description 59
- 230000005012 migration Effects 0.000 title claims abstract description 59
- 238000005070 sampling Methods 0.000 claims abstract description 180
- 230000002159 abnormal effect Effects 0.000 claims abstract description 56
- 238000001514 detection method Methods 0.000 claims description 22
- 238000012545 processing Methods 0.000 claims description 13
- 238000003384 imaging method Methods 0.000 abstract description 8
- 238000010586 diagram Methods 0.000 description 19
- 230000008569 process Effects 0.000 description 13
- 230000008859 change Effects 0.000 description 9
- 239000011435 rock Substances 0.000 description 8
- 230000000694 effects Effects 0.000 description 5
- 238000005457 optimization Methods 0.000 description 5
- 238000004880 explosion Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 239000002360 explosive Substances 0.000 description 2
- 230000008021 deposition Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000006467 substitution reaction 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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
- G01V2210/512—Pre-stack
-
- 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/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明实施例提供一种叠前深度偏移速度确定方法、装置及设备，根据多个勘探井中的多个采样地点的位置和预设网格，确定第一地点，多个第一地点和多个采样地点形成预设网格；分别获取多个采样地点和第一地点的横波速度和纵波速度；根据多个采样地点和第一地点的横波速度和纵波速度，在多个采样地点和多个第一地点中确定异常地点，异常地点对应的横波速度和/或纵波速度异常，提高了叠前深度偏移速度准确率，在横波速度剧烈变化的情况下可以得到更加精确的成像。
Description
技术领域
本发明实施例涉及地质勘探技术领域，尤其涉及一种叠前深度偏移速度确定方法、装置及设备。
背景技术
地震勘探是利用人工的方法引起地壳振动(如炸药爆炸、可控震源振动)，再用速度检测仪检测爆炸后各接收点的速度信息，对检测到的速度信息加工处理后得到地震资料，根据地震资料对地下岩层的构造进行成像，从而推断出地下岩层的性质和形态。
现有技术中，通常采用叠前深度偏移速度的方法对地下岩层的构造进行成像。在实际应用过程中，根据叠前深度偏移速度剖面速度的变化梯度和趋势判断是否有速度异常点，从而判断速度剖面是否合理。
然而，上述方法只考虑了一维空间速度的变化趋势，即只考虑了纵波速度，实际上，叠前深度偏移速度由于不同剖面横向变化而导致地震成像横向的差异比较大，因此，叠前深度偏移速度准确率较低。
发明内容
本发明实施例提供一种叠前深度偏移速度确定方法、装置及设备，用于解决叠前深度偏移速度准确率不高的问题。
第一方面，本发明实施例提供一种叠前深度偏移速度确定方法，包括：
根据所述多个勘探井中的多个采样地点的位置和预设网格，确定第一地点，所述多个第一地点和所述多个采样地点形成所述预设网格；
分别获取所述多个采样地点和所述第一地点的横波速度和纵波速度；
根据所述多个采样地点和所述第一地点的横波速度和纵波速度，在所述多个采样地点和所述多个第一地点中确定异常地点，所述异常地点对应的横波速度和/或纵波速度异常。
在一种可能的实施方式中，获取所述多个采样地点的纵波速度，包括：
针对任意一个勘探井，获取所述勘探井对应的检测速度，所述检测速度为所述勘探井中的速度检测器采集得到的速度；
根据所述检测速度和预设速度模型，确定所述勘探井中的多个采样地点的纵波速度。
在一种可能的实施方式中，所述采样地点包括第一采样地点和第二采样地点，所述第一采样地点位于所述勘探井内，所述第二采样点位于所述勘探井口；根据所述检测速度和预设速度模型，确定所述勘探井中的多个采样地点的纵波速度，包括：
根据所述检测速度和预设速度模型，确定所述第一采样地点的纵波速度；
将所述勘探井中的第一采样地点的纵波速度的平均值，确定为所述勘探井对应的第二采样点的纵波速度。
在一种可能的实施方式中，针对任意一个采样地点，获取所述采样地点的横波速度，包括：
获取所述采样地点对应的至少一个第一参考地点，所述第一参考地点与所述采样地点位于同一水平面；
按照预设规则对所述至少一个第一参考地点的纵波速度进行处理，得到所述采样地点的横波速度。
在一种可能的实施方式中，获取所述第一地点的纵波速度，包括：
确定所述第一地点对应的参考勘探井；
根据所述参考勘探井中的第一采样地点的纵波速度，确定所述第一地点的纵波速度，所述第一采样地点与所述第一地点位于同一水平面。
在一种可能的实施方式中，获取所述第一地点的横波速度，包括：
获取所述第一地点对应的至少一个第二参考地点，所述第二参考地点与所述第一地点位于同一水平面；
按照预设规则对所述至少一个第二参考地点的纵波速度进行处理，得到所述第一地点的横波速度。
在一种可能的实施方式中，根据所述多个采样地点和所述第一地点的横波速度和纵波速度，在所述多个采样地点和所述多个第一地点中确定异常地点，所述异常地点对应的横波速度和/或纵波速度异常，包括：
若所述多个采样地点和所述第一地点中的第二地点的横波速度与纵波速度的差值大于预设阈值，则确定所述第二地点为异常地点。
第二方面，本发明实施例提供一种叠前深度偏移速度确定装置，包括：
第一确定模块，用于根据所述多个勘探井中的多个采样地点的位置和预设网格，确定第一地点，所述多个第一地点和所述多个采样地点形成所述预设网格；
获取模块，用于分别获取所述多个采样地点和所述第一地点的横波速度和纵波速度；
第二确定模块，用于根据所述多个采样地点和所述第一地点的横波速度和纵波速度，在所述多个采样地点和所述多个第一地点中确定异常地点，所述异常地点对应的横波速度和/或纵波速度异常。
在一种可能的实施方式中，所述获取模块包括第一获取单元和第一确定单元，其中，
所述第一获取单元用于，针对任意一个勘探井，获取所述勘探井对应的检测速度，所述检测速度为所述勘探井中的速度检测器采集得到的速度；
所述第一确定单元用于，根据所述检测速度和预设速度模型，确定所述勘探井中的多个采样地点的纵波速度。
在一种可能的实施方式中，所述采样地点包括第一采样地点和第二采样地点，所述第一采样地点位于所述勘探井内，所述第二采样点位于所述勘探井口；所述第一确定单元具体用于：
根据所述检测速度和预设速度模型，确定所述第一采样地点的纵波速度；
将所述勘探井中的第一采样地点的纵波速度的平均值，确定为所述勘探井对应的第二采样点的纵波速度。
在一种可能的实施方式中，所述获取模块具体用于：
获取所述采样地点对应的至少一个第一参考地点，所述第一参考地点与所述采样地点位于同一水平面；
按照预设规则对所述至少一个第一参考地点的纵波速度进行处理，得到所述采样地点的横波速度。
在一种可能的实施方式中，所述获取模块具体用于：
确定所述第一地点对应的参考勘探井；
根据所述参考勘探井中的第一采样地点的纵波速度，确定所述第一地点的纵波速度，所述第一采样地点与所述第一地点位于同一水平面。
在一种可能的实施方式中，所述获取模块具体用于：
获取所述第一地点对应的至少一个第二参考地点，所述第二参考地点与所述第一地点位于同一水平面；
按照预设规则对所述至少一个第二参考地点的纵波速度进行处理，得到所述第一地点的横波速度。
在一种可能的实施方式中，所述第二确定模块具体用于：
若所述多个采样地点和所述第一地点中的第二地点的横波速度与纵波速度的差值大于预设阈值，则确定所述第二地点为异常地点。
第三方面，本发明实施例提供一种叠前深度偏移速度确定装置，包括：至少一个处理器和存储器；
所述存储器存储计算机执行指令；
所述至少一个处理器执行所述存储器存储的计算机执行指令，使得所述至少一个处理器执行如权利要求1至7任一项所述的叠前深度偏移速度确定方法。
第四方面，本发明实施例提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机执行指令，当处理器执行所述计算机执行指令时，实现如第一方面任一项所述的叠前深度偏移速度确定方法。
本发明实施例提供的一种叠前深度偏移速度确定方法、装置及设备，根据多个勘探井中的多个采样地点的位置和预设网格，确定第一地点，多个第一地点和多个采样地点形成预设网格；分别获取多个采样地点和第一地点的横波速度和纵波速度；根据多个采样地点和第一地点的横波速度和纵波速度，在多个采样地点和多个第一地点中确定异常地点，异常地点对应的横波速度和/或纵波速度异常，因此，通过上述方法不仅考虑到纵波速度的变化，同时也考虑了横波速度的变化趋势，提高了叠前深度偏移速度准确率，在横波速度剧烈变化的情况下可以得到更加精确的成像。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种叠前深度偏移速度确定方法的流程示意图；
图2为本发明实施例提供的预设网格的结构示意图；
图3为本发明实施例提供的纵波速度标定方法的示意图；
图4为本发明实施例提供的速度拟三维插值的网格示意图；
图5为本发明实施例提供的异常地点的三维空间速度场中二维速度剖面图；
图6为本发明实施例提供的优化前后V4平面速度对比图；
图7为本发明实施例提供的优化速度前后地震剖面成像效果对比图；
图8为本发明实施例提供的另一种叠前深度偏移速度确定方法的流程示意图；
图9为任意一个勘探井中待获取纵波速度的多个采样地点的示意图；
图10为本发明实施例提供的任意水平面内采样地点和第一参考地点示意图；
图11为本发明实施例提供的二维地震速度剖面示意图；
图12为本发明实施例提供的一种叠前深度偏移速度确定装置的结构示意图；
图13为本发明实施例提供的另一种叠前深度偏移速度确定装置的结构示意图；
图14为本发明实施例提供的叠前深度偏移速度确定装置的硬件结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优地点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
目前，通常采用地震勘探的方式进行油气勘探，地震勘探是利用人工的方法引起地壳振动(如炸药爆炸、可控震源振动)，再用速度检测仪检测爆炸后各接收点的速度信息，对检测到的速度信息加工处理后得到地震资料，根据地震资料对地下岩层的构造进行成像，从而推断出地下岩层的性质和形态。
现有技术中，通常采用叠前深度偏移速度的方法对地下岩层的构造进行成像。在实际应用过程中，根据叠前深度偏移速度剖面速度的变化梯度和趋势判断是否有速度异常点，从而判断速度剖面是否合理。
上述方法只考虑了一维空间速度的变化趋势，即只考虑了纵波速度，然而油气勘探面临越来越复杂的地表条件和地下地质构造特征，叠前深度偏移速度由于不同剖面横向变化而导致地震成像横向的差异比较大，当介质存在明显横波速度变化时，由于地震波传播遇到速度变化明显界面时会产生绕射的影响，现有的叠前深度偏移速度不能使反射层正确归位，叠前深度偏移速度准确率不高，在横波速度剧烈变化的情况下得不到精确的成像。因此，对复杂条件下，构造成像的方法提出了更高的要求。
本发明通过根据多个勘探井中的多个采样地点的位置和预设网格，确定第一地点，多个第一地点和多个采样地点形成预设网格；分别获取多个采样地点和第一地点的横波速度和纵波速度；根据多个采样地点和第一地点的横波速度和纵波速度，在多个采样地点和多个第一地点中确定异常地点，异常地点对应的横波速度和/或纵波速度异常，并对异常地点进行质控，通过上述方法不仅考虑到纵波速度的变化，同时也考虑了横波速度的变化趋势，提高了叠前深度偏移速度准确率，在横波速度剧烈变化的情况下可以得到更加精确的成像。
下面，通过具体实施例对本申请所示的技术方案进行详细说明。需要说明的是，下面几个具体实施例可以相互结合，对于相同或相似的内容，在不同的实施例中不再进行重复说明。
图1为本发明实施例提供的一种叠前深度偏移速度确定方法的流程示意图。请参见图1，该方法可以包括：
S101、根据多个勘探井中的多个采样地点的位置和预设网格，确定第一地点，多个第一地点和多个采样地点形成预设网格。
本发明实施例的执行主体可以为终端设备或者服务器，也可以为设置在终端设备或者服务器中的叠前深度偏移速度确定装置。可选的，终端设备可以为手机、电脑等。叠前深度偏移速度确定装置可以通过软件实现，也可以通过软件和硬件的结合实现。为了便于描述，下面为执行主体为终端设备为例进行说明。
其中，勘探井是一种地质调查过程中用于勘探而钻探的探井，当基岩埋藏较深，探槽无法达到地质目的或受地面条件影响，无法施工探槽时采用的一种占地面积较小的浅部地质工程，可给调查者提供关于地下条件的重要资料，例如，可以包括岩层描述、地下水位深度、水质和水井的性能。本发明实施例中，勘探井的设置主要为了在油气勘探的圈闭预探阶段，确定岩层构造、形态和查明地层层序及接触关系。
采样地点可以为勘探井中多个可以放置速度检测仪的用来检测速度的地点。
预设网格为叠前深度偏移速度确定装置重速度拟三维差值的栅格点所组成的网格，图2为本发明实施例提供的预设网格的结构示意图，请参见图2，该预设网格可以通过获取到的各栅格点的速度，得到各速度层平面速度图，实现对地下岩层的构造进行成像。
第一地点可以为通过采样地点确定的除采样地点之外其他的地点。多个第一地点和多个采样地点可以形成预设网格。
S102、分别获取多个采样地点和第一地点的横波速度和纵波速度。
纵波速度为只考虑一维空间速度的变化趋势而近似得到的地震速度，大多数沉积地区将地震界面的上覆层近似看成连续介质而无线束差异从而使得地震速度的研究得以简化。
横波速度为地震速度的横波速度分量，用于指示地震速度的横向变化。在地质情况复杂、速度多变的地区中，只考虑一维空间速度的变化趋势会存在较大的误差，精确获取地震速度的横波速度可以提高叠前深度偏移速度的精度。
可选的，在S102之前，可以通过井速度标定方法验证现有的叠前深度偏移速度的合理性，再通过叠前深度偏移速度获取多个采样地点和第一地点的纵波速度。下面，结合图3，对通过井速度标定方法验证现有的叠前深度偏移速度的合理性进行说明。
图3为本发明实施例提供的纵波速度标定方法的示意图，请参见图3，该方法可以验证现有的叠前深度偏移速度模型是否合理。利用工区已完钻井S1、S2井声波和垂直地震测深(Vertical Seismic Profiling，VSP)资料，与利用叠前深度偏移速度模型获取的钻井S1、S2的纵波速度进行对比拟合，判断现有的叠前深度偏移速度模型的合理性。从图3所示的两口井声波速度与纵波速度的对比来看，两者整体趋势基本一致，因此认为利用该速度模型获取得到的纵波速度是合理的。即，通过现有的叠前深度偏移速度模型可以得到待测区多个采样地点和第一地点的纵波速度，该纵波速度是合理的。
因此，利用现有的叠前深度偏移速度模型可以获取多个采样地点和第一地点的纵波速度。
S103、根据多个采样地点和第一地点的横波速度和纵波速度，在多个采样地点和多个第一地点中确定异常地点，异常地点对应的横波速度和/或纵波速度异常。
若多个采样地点和第一地点中的第二地点的横波速度与纵波速度的差值大于预设阈值，则确定第二地点为异常地点。
可选的，在S103之前，可以先进行速度拟三维插值。图4为本发明实施例提供的速度拟三维插值的网格示意图，请参见图4，H1、H2为通过S101、S102获取得到纵波速度和横波速度的平面层，Hn为待求取速度的拟三维插值的平面层。具体的，可以通过如下方式进行速度拟三维插值：
式中，Vn为待求取速度，VH1、VH2为上下两层速度，ZH1、ZH2为上下两层深度值。
可选的，速度拟三维插值之后，可以通过散点网格成图得到三维空间速度场。通过三维空间速度场可以更加直观地确定异常地点。
图5为本发明实施例提供的异常地点的三维空间速度场中二维速度剖面图。请参见图5，其中，(A)为三维空间速度场中V4平面速度图，(B)为三维空间速度场中V5平面速度图，椭圆形框中所示区域即指示该区域中的地点为异常地点。从图5可以看出，通过三维空间速度场可以更加直观的在多个采样地点和多个第一地点中确定异常地点。
可选的，在确定异常地点之后，可以对三维空间速度场中异常地点的横波速度与纵波速度进行优化。
可选的，可以通过如下方式将找到的异常地点的横波速度与纵波速度进行优化：
将找到的异常地点的纵波速度与横波速度做第一次平均，得到平面速度图，对该平面速度图进行质控，如果横波速度变化趋势不符合地质规律，则继续对第一次平均得到的速度与横波速度再做一次平均，不断迭代优化直至合理。当迭代优化后的误差小于第二预设阈值时，则认为优化后的速度场为合理的。下面，结合图6，对优化前后的速度场进行对比说明。
图6为本发明实施例提供的优化前后V4平面速度对比图，其中，(A)为本发明实施例提供的优化前V4平面速度图，(B)为本发明实施例提供的优化后V4平面速度图。请参见图6，可以看出，对异常地点的横波速度与纵波速度进行优化后，该地点的地震速度趋于正常。
图7为本发明实施例提供的优化速度前后地震剖面成像效果对比图。请参见图7，其中，(A)为优化前地震剖面成像效果图，(B)为优化后地震剖面成像效果图。通过本发明实施例提供的一种叠前深度偏移速度确定方法，可以在地质情况复杂、速度多变的地区，得到更加精确的成像。
本发明实施例提供的一种叠前深度偏移速度确定方法、装置及设备，根据多个勘探井中的多个采样地点的位置和预设网格，确定第一地点，多个第一地点和多个采样地点形成预设网格；分别获取多个采样地点和第一地点的横波速度和纵波速度；根据多个采样地点和第一地点的横波速度和纵波速度，在多个采样地点和多个第一地点中确定异常地点，异常地点对应的横波速度和/或纵波速度异常，因此，通过上述方法不仅考虑到纵波速度的变化，同时也考虑了横波速度的变化趋势，提高了叠前深度偏移速度准确率，在横波速度剧烈变化的情况下可以得到更加精确的成像。
图8为本发明实施例提供的另一种叠前深度偏移速度确定方法的流程示意图。请参见图8，该方法可以包括：
S201、根据多个勘探井中的多个采样地点的位置和预设网格，确定第一地点，多个第一地点和多个采样地点形成预设网格。
需要说明的是，S201的执行过程可以参见S101的执行过程，此处不再进行赘述。
S202、获取多个采样地点的纵波速度。
采样地点包括第一采样地点和第二采样地点，第一采样地点位于勘探井内，第二采样点位于勘探井口。第一采样地点的纵波速度可以通过速度检测器以及预设速度模型得到，第二采样地点的纵波速度可以通过对第一采样地点的纵波速度求取平均速度得到。下面，对该过程进行详细说明：
图9为任意一个勘探井中待获取纵波速度的多个采样地点的示意图，请参见图9，V2-Vi分别表示第一采样地点的纵波速度，采样地点V1表示第二采样地点的纵波速度。
首先，针对任意一个勘探井，获取勘探井对应的检测速度，检测速度为勘探井中的速度检测器采集得到的地震波的纵波速度。
其次，根据检测速度和预设速度模型，确定第一采样地点V2-Vi的纵波速度。预设速度模型即为现有的叠前深度偏移速度模型，该模型的合理性已通过井速度标定方法得以验证。将检测速度输入该预设速度模型，可以得到第一采样地点V1的纵波速度。
最后，将勘探井中的第一采样地点的纵波速度V2-Vi的平均值，确定为勘探井对应的第二采样点的纵波速度V1。例如，可以通过如下公式确定第二采样点的纵波速度：
式中，VRMS为待求取的第二采样点的纵波速度，N为采样点个数，Vi为第一采样点的纵波速度；
采用上述公式依次求取多个勘探井中多个采样地点的纵波速度。
通过上述公式可以计算得到勘探井中多个采样地点的纵波速度。通过勘探井中速度检测器采集得到检测速度以及预设叠前深度偏移速度模型，可以确定第一采样地点的纵波速度，根据第一采样地点的纵波速度求取平均值可以获取第二采样地点的纵波速度，因此，可以确定勘探井中多个采样地点的纵波速度。
S203、针对任意一个采样地点，获取采样地点的横波速度。
获取采样地点对应的至少一个第一参考地点，第一参考地点与采样地点位于同一水平面；
按照预设规则对至少一个第一参考地点的纵波速度进行处理，得到采样地点的横波速度。
下面，结合图10，对通过采样地点的纵波速度获取采样地点横波速度的过程进行说明。
图10为本发明实施例提供的任意水平面内采样地点和第一参考地点示意图。请参见图10，V(x，y)所在地点为待求取横波速度的采样地点，V(xi,yi)所在地点为已知纵波速度第一参考地点，其中，i可以为任意整数，图中以4示出，第一参考地点可以为不同勘探井中同一平面的采样地点。
可选的，按照通过如下预设规则对至少一个第一参考地点的纵波速度进行处理，得到采样地点的横波速度：
式中V(x，y)为待求取的采样地点的横波速度，V(xi,yi)为第一参考地点的纵波速度，λi为距离加权因子，k为比例系数，N为采样点数。
通过上述公式可以根据第一参考地点的纵波速度，计算得到采样地点的横波速度。
S204、获取第一地点的纵波速度。
确定第一地点对应的参考勘探井，参考勘探井的数量可以为一个或者多个，第一地点为工区同一平面内，需求取地震速度(包括纵波速度和横波速度)的除第一采样地点之外其他的地点。
确定好第一地点对应的参考勘探井之后，将参考勘探井中第一采样地点的纵波速度输入叠前深度偏移速度模型，得到第一地点的纵波速度。
S205、获取第一地点的横波速度。
获取第一地点对应的至少一个第二参考地点，第二参考地点与第一地点位于同一水平面，第二参考地点可以为已获取纵波速度的第一地点中的任意地点，第二参考地点的数量可以为1个或者多个；
按照预设规则对至少一个第二参考地点的纵波速度进行处理，得到第一地点的横波速度。
下面，结合图11，以水平面为6层、第二参考地点为4个，对按照预设规则得到第一地点的横波速度的过程进行详细说明。
图11为本发明实施例提供的二维地震速度剖面示意图。请参见图11：
首先，在工区内建立若干平面层，例如，平面层可以为V1-V6；
其次，对工区内二维地震速度剖面每个平面层进行划线标记，每个平面层中第二参考地点的纵波速度已知，具体计算过程请参见上述步骤；
最后，通过各平面层的第二参考地点的纵波速度，得到V1-V6共6个速度层的第二参考地点的横波速度平面图。
可选的，可以通过如下公式对第二参考地点的纵波速度进行处理，得到第二参考地点的横波速度：
式中V(x，y)为待求取的第二参考地点的横波速度，V(xi,yi)为第二参考地点的纵波速度，λi为距离加权因子，k为比例系数，N为参考地点的点数。
通过上述公式可以根据第二参考地点的纵波速度，计算得到第二参考地点的横波速度。
S206、根据多个采样地点和第一地点的横波速度和纵波速度，在多个采样地点和多个第一地点中确定异常地点。
需要说明的是，S206的执行过程可以参见S103的执行过程，此处不再进行赘述。
本发明实施例提供的一种叠前深度偏移速度确定方法、装置及设备，根据多个勘探井中的多个采样地点的位置和预设网格，确定第一地点，多个第一地点和多个采样地点形成预设网格；分别获取多个采样地点和第一地点的横波速度和纵波速度；根据多个采样地点和第一地点的横波速度和纵波速度，在多个采样地点和多个第一地点中确定异常地点，异常地点对应的横波速度和/或纵波速度异常，因此，通过上述方法不仅考虑到纵波速度的变化，同时也考虑了横波速度的变化趋势，提高了叠前深度偏移速度准确率，在横波速度剧烈变化的情况下可以得到更加精确的成像。
图12为本发明实施例提供的一种叠前深度偏移速度确定装置的结构示意图。请参见图12，该装置包括：
第一确定模块11，用于根据所述多个勘探井中的多个采样地点的位置和预设网格，确定第一地点，所述多个第一地点和所述多个采样地点形成所述预设网格；
获取模块12，用于分别获取所述多个采样地点和所述第一地点的横波速度和纵波速度；
第二确定模块13，用于根据所述多个采样地点和所述第一地点的横波速度和纵波速度，在所述多个采样地点和所述多个第一地点中确定异常地点，所述异常地点对应的横波速度和/或纵波速度异常。
本发明实施例提供的叠前深度偏移速度确定装置10可以执行上述方法实施例所示的技术方案，其实现原理以及有益效果类似，此处不再进行赘述。
图13为本发明实施例提供的另一种叠前深度偏移速度确定装置的结构示意图，在图12所示实施例的基础上，请参见图13，所述获取模块12包括获取单元121和确定单元122，其中，
所述获取单元121用于，针对任意一个勘探井，获取所述勘探井对应的检测速度，所述检测速度为所述勘探井中的速度检测器采集得到的速度；
所述确定单元122用于，根据所述检测速度和预设速度模型，确定所述勘探井中的多个采样地点的纵波速度。
在另一种可能的实施方式中，所述采样地点包括第一采样地点和第二采样地点，所述第一采样地点位于所述勘探井内，所述第二采样点位于所述勘探井口；所述确定单元122具体用于：
根据所述检测速度和预设速度模型，确定所述第一采样地点的纵波速度；
将所述勘探井中的第一采样地点的纵波速度的平均值，确定为所述勘探井对应的第二采样点的纵波速度。
在另一种可能的实施方式中，所述获取模块12具体用于：
获取所述采样地点对应的至少一个第一参考地点，所述第一参考地点与所述采样地点位于同一水平面；
按照预设规则对所述至少一个第一参考地点的纵波速度进行处理，得到所述采样地点的横波速度。
在另一种可能的实施方式中，所述获取模块12具体用于：
确定所述第一地点对应的参考勘探井；
根据所述参考勘探井中的第一采样地点的纵波速度，确定所述第一地点的纵波速度，所述第一采样地点与所述第一地点位于同一水平面。
在另一种可能的实施方式中，所述获取模块12具体用于：
获取所述第一地点对应的至少一个第二参考地点，所述第二参考地点与所述第一地点位于同一水平面；
按照预设规则对所述至少一个第二参考地点的纵波速度进行处理，得到所述第一地点的横波速度。
在另一种可能的实施方式中，所述第二确定模块13具体用于：
若所述多个采样地点和所述第一地点中的第二地点的横波速度与纵波速度的差值大于预设阈值，则确定所述第二地点为异常地点。
本发明实施例提供的叠前深度偏移速度确定装置10可以执行上述方法实施例所示的技术方案，其实现原理以及有益效果类似，此处不再进行赘述。
图14为本发明实施例提供的叠前深度偏移速度确定装置的硬件结构示意图，如图14所示，该叠前深度偏移速度确定装置20包括：至少一个处理器21和存储器22。其中，处理器21和存储器22通过总线23连接。
在具体实现过程中，至少一个处理器21执行所述存储器22存储的计算机执行指令，使得至少一个处理器21执行如上的叠前深度偏移速度确定方法。
处理器21的具体实现过程可参见上述方法实施例，其实现原理和技术效果类似，本实施例此处不再赘述。
在上述图14所示的实施例中，应理解，处理器可以是中央处理单元(英文：CentralProcessing Unit，简称：CPU)，还可以是其他通用处理器、数字信号处理器(英文：DigitalSignal Processor，简称：DSP)、专用集成电路(英文：Application Specific IntegratedCircuit，简称：ASIC)等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合发明所公开的方法的步骤可以直接体现为硬件处理器执行完成，或者用处理器中的硬件及软件模块组合执行完成。
存储器可能包含高速RAM存储器，也可能还包括非易失性存储NVM，例如至少一个磁盘存储器。
总线可以是工业标准体系结构(Industry Standard Architecture，ISA)总线、外部设备互连(Peripheral Component，PCI)总线或扩展工业标准体系结构(ExtendedIndustry Standard Architecture，EISA)总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示，本申请附图中的总线并不限定仅有一根总线或一种类型的总线。
本申请还提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机执行指令，当处理器执行所述计算机执行指令时，实现如上所述的叠前深度偏移速度确定方法。
上述的计算机可读存储介质，上述可读存储介质可以是由任何类型的易失性或非易失性存储设备或者它们的组合实现，如静态随机存取存储器(SRAM)，电可擦除可编程只读存储器(EEPROM)，可擦除可编程只读存储器(EPROM)，可编程只读存储器(PROM)，只读存储器(ROM)，磁存储器，快闪存储器，磁盘或光盘。可读存储介质可以是通用或专用计算机能够存取的任何可用介质。
一种示例性的可读存储介质耦合至处理器，从而使处理器能够从该可读存储介质读取信息，且可向该可读存储介质写入信息。当然，可读存储介质也可以是处理器的组成部分。处理器和可读存储介质可以位于专用集成电路(Application Specific IntegratedCircuits，简称：ASIC)中。当然，处理器和可读存储介质也可以作为分立组件存在于设备中。
所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
本领域普通技术人员可以理解：实现上述各方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程序在执行时，执行包括上述各方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (10)
1.一种叠前深度偏移速度确定方法，其特征在于，包括：
根据所述多个勘探井中的多个采样地点的位置和预设网格，确定第一地点，所述多个第一地点和所述多个采样地点形成所述预设网格；
分别获取所述多个采样地点和所述第一地点的横波速度和纵波速度；
根据所述多个采样地点和所述第一地点的横波速度和纵波速度，在所述多个采样地点和所述多个第一地点中确定异常地点，所述异常地点对应的横波速度和/或纵波速度异常。
2.根据权利要求1所述的方法，其特征在于，获取所述多个采样地点的纵波速度，包括：
针对任意一个勘探井，获取所述勘探井对应的检测速度，所述检测速度为所述勘探井中的速度检测器采集得到的速度；
根据所述检测速度和预设速度模型，确定所述勘探井中的多个采样地点的纵波速度。
3.根据权利要求2所述的方法，其特征在于，所述采样地点包括第一采样地点和第二采样地点，所述第一采样地点位于所述勘探井内，所述第二采样点位于所述勘探井口；根据所述检测速度和预设速度模型，确定所述勘探井中的多个采样地点的纵波速度，包括：
根据所述检测速度和预设速度模型，确定所述第一采样地点的纵波速度；
将所述勘探井中的第一采样地点的纵波速度的平均值，确定为所述勘探井对应的第二采样点的纵波速度。
4.根据权利要求1-3任一项所述的方法，其特征在于，针对任意一个采样地点，获取所述采样地点的横波速度，包括：
获取所述采样地点对应的至少一个第一参考地点，所述第一参考地点与所述采样地点位于同一水平面；
按照预设规则对所述至少一个第一参考地点的纵波速度进行处理，得到所述采样地点的横波速度。
5.根据权利要求1-3任一项所述的方法，其特征在于，获取所述第一地点的纵波速度，包括：
确定所述第一地点对应的参考勘探井；
根据所述参考勘探井中的第一采样地点的纵波速度，确定所述第一地点的纵波速度，所述第一采样地点与所述第一地点位于同一水平面。
6.根据权利要求1-3任一项所述的方法，其特征在于，获取所述第一地点的横波速度，包括：
获取所述第一地点对应的至少一个第二参考地点，所述第二参考地点与所述第一地点位于同一水平面；
按照预设规则对所述至少一个第二参考地点的纵波速度进行处理，得到所述第一地点的横波速度。
7.根据权利要求1-3任一项所述的方法，其特征在于，根据所述多个采样地点和所述第一地点的横波速度和纵波速度，在所述多个采样地点和所述多个第一地点中确定异常地点，所述异常地点对应的横波速度和/或纵波速度异常，包括：
若所述多个采样地点和所述第一地点中的第二地点的横波速度与纵波速度的差值大于预设阈值，则确定所述第二地点为异常地点。
8.一种叠前深度偏移速度确定装置，其特征在于，包括：
第一确定模块，用于根据所述多个勘探井中的多个采样地点的位置和预设网格，确定第一地点，所述多个第一地点和所述多个采样地点形成所述预设网格；
获取模块，用于分别获取所述多个采样地点和所述第一地点的横波速度和纵波速度；
第二确定模块，用于根据所述多个采样地点和所述第一地点的横波速度和纵波速度，在所述多个采样地点和所述多个第一地点中确定异常地点，所述异常地点对应的横波速度和/或纵波速度异常。
9.一种叠前深度偏移速度确定装置，其特征在于，包括：至少一个处理器和存储器；
所述存储器存储计算机执行指令；
所述至少一个处理器执行所述存储器存储的计算机执行指令，使得所述至少一个处理器执行如权利要求1至7任一项所述的叠前深度偏移速度确定方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有计算机执行指令，当处理器执行所述计算机执行指令时，实现如权利要求1至7任一项所述的叠前深度偏移速度确定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911347263.9A CN113031073B (zh) | 2019-12-24 | 2019-12-24 | 叠前深度偏移速度确定方法、装置及设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911347263.9A CN113031073B (zh) | 2019-12-24 | 2019-12-24 | 叠前深度偏移速度确定方法、装置及设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113031073A true CN113031073A (zh) | 2021-06-25 |
| CN113031073B CN113031073B (zh) | 2024-03-01 |
Family
ID=76451743
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911347263.9A Active CN113031073B (zh) | 2019-12-24 | 2019-12-24 | 叠前深度偏移速度确定方法、装置及设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113031073B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4736349A (en) * | 1987-04-24 | 1988-04-05 | Mobil Oil Corporation | Method for estimating shear wave reflection data from acquired compressional wave reflection data |
| US4870580A (en) * | 1983-12-30 | 1989-09-26 | Schlumberger Technology Corporation | Compressional/shear wave separation in vertical seismic profiling |
| CN101419292A (zh) * | 2007-10-25 | 2009-04-29 | 中国石油天然气集团公司 | 采用纵波源多分量地震数据生成横波地震剖面的方法 |
| CN103076629A (zh) * | 2012-09-21 | 2013-05-01 | 中国石油天然气集团公司 | 基于纵横地震波速度比进行油气勘探的方法及设备 |
| CN103760603A (zh) * | 2014-01-28 | 2014-04-30 | 中国石油大学(北京) | 转换波地震数据的叠前时间偏移方法及装置 |
| CN106959467A (zh) * | 2017-03-20 | 2017-07-18 | 中国科学院地质与地球物理研究所 | 地震波速度反演方法和装置 |
-
2019
- 2019-12-24 CN CN201911347263.9A patent/CN113031073B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4870580A (en) * | 1983-12-30 | 1989-09-26 | Schlumberger Technology Corporation | Compressional/shear wave separation in vertical seismic profiling |
| US4736349A (en) * | 1987-04-24 | 1988-04-05 | Mobil Oil Corporation | Method for estimating shear wave reflection data from acquired compressional wave reflection data |
| CN101419292A (zh) * | 2007-10-25 | 2009-04-29 | 中国石油天然气集团公司 | 采用纵波源多分量地震数据生成横波地震剖面的方法 |
| CN103076629A (zh) * | 2012-09-21 | 2013-05-01 | 中国石油天然气集团公司 | 基于纵横地震波速度比进行油气勘探的方法及设备 |
| CN103760603A (zh) * | 2014-01-28 | 2014-04-30 | 中国石油大学(北京) | 转换波地震数据的叠前时间偏移方法及装置 |
| CN106959467A (zh) * | 2017-03-20 | 2017-07-18 | 中国科学院地质与地球物理研究所 | 地震波速度反演方法和装置 |
Non-Patent Citations (1)
| Title |
|---|
| 宁斌 等: "转换波叠前深度偏移成像应用研究", 2015年物探技术研讨会, pages 717 - 720 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113031073B (zh) | 2024-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2683618C (en) | Inverse-vector method for smoothing dips and azimuths | |
| US6721694B1 (en) | Method and system for representing the depths of the floors of the oceans | |
| US8873333B2 (en) | Geophysical signal processing | |
| US20160018543A1 (en) | Quality check of compressed data sampling interpolation for seismic information | |
| EA032186B1 (ru) | Сейсмическая адаптивная фокусировка | |
| US10534100B2 (en) | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on time-lapse seismic data | |
| CN113267822B (zh) | 基于地形约束因子权重优化提高海底地形反演精度的方法 | |
| AU2009311479A1 (en) | Systems and methods for computing and validating a variogram model | |
| US10884149B2 (en) | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic data | |
| CN107831536B (zh) | 沙漠地表地震资料处理浮动基准面建立方法 | |
| CN111208558B (zh) | 超深低幅度三维地质构造的建立方法及装置 | |
| CN113031073A (zh) | 叠前深度偏移速度确定方法、装置及设备 | |
| US10718876B2 (en) | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic inversions | |
| US10067254B2 (en) | Removal of an estimated acquisition effect from a marine survey measurement | |
| US20220221605A1 (en) | Decompaction of subsurface region using decompaction velocity | |
| WO2018207090A1 (en) | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic data | |
| US11573347B2 (en) | Computing program product and method that interpolates wavelets coefficients and estimates spatial varying wavelets using the covariance interpolation method in the data space over a survey region having multiple well locations | |
| US20240168189A1 (en) | Systems, devices, and methods for generating average velocity maps of subsurface formations | |
| CN114355453A (zh) | 一种地质勘探方法及终端设备 | |
| CN115184987A (zh) | 地层倾角信息的提取方法、装置、服务器及存储介质 | |
| US10466379B2 (en) | Method for enhancing the resolution of a seismic image | |
| CN116400412A (zh) | 地震层位提取方法、装置、电子设备及计算机存储介质 | |
| CN111142160A (zh) | 时间推移地震观测数据的分析方法及装置 | |
| CN115113273A (zh) | 裂缝参数的确定方法、装置、设备及存储介质 | |
| CN115576023A (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 |