CN108508481A - 一种纵波转换波地震数据时间匹配的方法、装置及系统 - Google Patents
一种纵波转换波地震数据时间匹配的方法、装置及系统 Download PDFInfo
- Publication number
- CN108508481A CN108508481A CN201810471934.1A CN201810471934A CN108508481A CN 108508481 A CN108508481 A CN 108508481A CN 201810471934 A CN201810471934 A CN 201810471934A CN 108508481 A CN108508481 A CN 108508481A
- Authority
- CN
- China
- Prior art keywords
- wave
- depth
- time
- converted
- seismic data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 74
- 238000006243 chemical reaction Methods 0.000 claims abstract description 37
- 238000012545 processing Methods 0.000 claims abstract description 23
- 238000012937 correction Methods 0.000 claims description 14
- 238000013213 extrapolation Methods 0.000 claims description 12
- 238000003860 storage Methods 0.000 claims description 12
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 241001269238 Data Species 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 14
- 230000006870 function Effects 0.000 description 13
- 230000005012 migration Effects 0.000 description 9
- 238000013508 migration Methods 0.000 description 9
- 230000015572 biosynthetic process Effects 0.000 description 7
- 238000004590 computer program Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 238000005755 formation reaction Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 7
- 238000004891 communication Methods 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000000750 progressive effect Effects 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000011017 operating method Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001131 transforming effect Effects 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
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
本申请实施例公开了一种纵波转换波地震数据时间匹配的方法、装置及系统，所述方法包括对纵波以及转换波的原始时间域地震数据分别进行时间深度转换，获得纵波以及转换波的深度域地震数据；根据所述纵波以及转换波的深度域地震数据分别计算纵波以及转换波在深度域的地震能量数据；基于预设的纵波与转换波在相同深度处的地震能量数据体最大相关性规则对转换波的深度域地震数据进行深度修正处理；将深度修正处理后的所述转换波的深度域地震数据转换到纵波时间域，获得时间匹配后的转换波时间域地震数据。利用本申请各个实施例，可以提高纵波转换波地震数据时间匹配的效率与准确性。
Description
技术领域
本发明涉及石油天然气地震勘探技术领域，特别地，涉及一种纵波转换波地震数据时间匹配的方法、装置及系统。
背景技术
近年来，转换波地震勘探技术在全球范围内呈现出较快的发展趋势，该技术在岩性油气藏、隐蔽性油气藏、非常规油气藏等领域有很好的应用前景。然而，要充分利用转换波地震勘探资料提取纵横波时差比、振幅比、纵横波速度比、泊松比、体积模量、剪切模量、拉梅系数等油气储层的岩性与物性参数，实现储层预测、裂缝检测及流体识别等勘探目标，首先需要解决纵波转换波地震剖面的时间匹配问题。
纵波转换波地震剖面的时间匹配问题是由于转换波传播路径的非对称性引起的，这导致纵波与转换波地震资料的联合解释困难。由于地震纵波与地震横波的传播速度不同，造成同一个地下地质层位在纵波地震资料和转换波地震资料上的传播时间相差很大，因此阻碍了技术人员对地下同一层位利用纵波与转换波地震资料进行联合解释。因此，如何将纵波地震剖面与转换波地震剖面进行时间匹配(使相同地下地质层位在纵波和转换波地震资料上的传播时间一致)，一直是困扰油气勘探工程师的技术难题。传统的解决方案是在时间域通过人工交互解释的方法，将纵波与转换波地震剖面进行时间匹配，但这种方法基于对勘探区域地下地层的经验认识，误差较大。
发明内容
本申请实施例的目的在于提供一种纵波转换波地震数据时间匹配的方法、装置及系统，可以提高纵波转换波地震数据时间匹配的效率与准确性。
本申请提供的一种纵波转换波地震数据时间匹配的方法、装置及系统是通过包括以下方式实现的：
一种纵波转换波地震数据时间匹配的方法，包括：
对纵波以及转换波的原始时间域地震数据分别进行时间深度转换，获得纵波以及转换波的深度域地震数据；
根据所述纵波以及转换波的深度域地震数据分别计算纵波以及转换波在深度域的地震能量数据；
基于预设的纵波与转换波在相同深度处的地震能量数据体最大相关性规则对转换波的深度域地震数据进行深度修正处理；
将深度修正处理后的所述转换波的深度域地震数据转换到纵波时间域，获得时间匹配后的转换波时间域地震数据。
本申请实施例的纵波转换波地震数据时间匹配的方法，所述对纵波的原始时间域地震数据进行时间深度转换，包括：
对纵波的原始时间域地震数据进行层位分析，获得纵波时间层位数据；
根据纵波均方根速度对纵波时间层位数据进行时间深度转换，获得纵波深度层位数据；
对所述纵波深度层位数据利用二维内插以及外推算法确定纵波的深度域地震数据。
本申请实施例的纵波转换波地震数据时间匹配的方法，所述对转换波的原始时间域地震数据进行时间深度转换，包括：
对转换波的原始时间域地震数据进行层位分析，获得转换波时间层位数据；
根据转换波均方根速度以及纵横波垂直速度比值数据对转换波时间层位数据进行时间深度转换，获得转换波深度层位数据；
对所述转换波深度层位数据利用二维内插以及外推算法确定转换波的深度域地震数据。
本申请实施例的纵波转换波地震数据时间匹配的方法，所述根据纵波均方根速度对纵波时间层位数据进行时间深度转换，包括：
根据如下公式对纵波时间层位数据进行时间深度转换：
其中，表示纵波深度层位数据，表示纵波均方根速度，表示纵波时间层位数据，i表示时间层位序号，j表示每一时间层位的控制点序号。
本申请实施例的纵波转换波地震数据时间匹配的方法，所述根据转换波均方根速度以及纵横波垂直速度比值数据对转换波时间层位数据进行时间深度转换，包括：
根据如下公式对转换波时间层位数据进行时间深度转换：
其中，表示纵波深度层位数据，表示纵波均方根速度，表示纵波时间层位数据，表示纵横波垂直速度比值，i表示时间层位序号，j表示每一时间层位的控制点序号。
本申请实施例的纵波转换波地震数据时间匹配的方法，所述分别计算纵波以及转换波在深度域的地震能量数据，包括：
对所述纵波以及转换波的深度域地震数据分别进行希尔伯特变换；
根据纵波以及转换波希尔伯特变换后的深度域地震数据确定纵波以及转换波在深度域的地震能量数据。
本申请实施例的纵波转换波地震数据时间匹配的方法，所述对转换波的深度域地震数据进行深度修正处理，包括：
计算转换波深度层位的能量数据体与纵波各深度层位的能量数据体的相关值；
根据最大相关值对应的纵波深度层位数据确定所述转换波深度层位的最优深度修正值，利用所述最优深度修正值对转换波的深度域地震数据的深度进行修正。
另一方面，本申请实施例还提供一种纵波转换波地震数据时间匹配的装置，包括：
时深转换模块，用于对纵波以及转换波的原始时间域地震数据分别进行时间深度转换，获得纵波以及转换波的深度域地震数据；
能量计算模块，用于根据所述纵波以及转换波的深度域地震数据分别计算纵波以及转换波在深度域的地震能量数据；
深度修正模块，用于基于预设的纵波与转换波在相同深度处的地震能量数据体最大相关性规则对转换波的深度域地震数据进行深度修正处理；
时间匹配模块，用于将深度修正处理后的所述转换波的深度域地震数据转换到纵波时间域，获得时间匹配后的转换波时间域地震数据。
本申请实施例的一种纵波转换波地震数据时间匹配的装置，包括处理器及用于存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：
对纵波以及转换波的原始时间域地震数据分别进行时间深度转换，获得纵波以及转换波的深度域地震数据；
根据所述纵波以及转换波的深度域地震数据分别计算纵波以及转换波在深度域的地震能量数据；
基于预设的纵波与转换波在相同深度处的地震能量数据体最大相关性规则对转换波的深度域地震数据进行深度修正处理；
将深度修正处理后的所述转换波的深度域地震数据转换到纵波时间域，获得时间匹配后的转换波时间域地震数据。
另一方面，本申请实施例还提供一种纵波转换波地震数据时间匹配的系统，包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现上述任意一个实施例中所述方法的步骤。
本说明书一个或多个实施例提供的一种纵波转换波地震数据时间匹配的方法、装置及系统，可以通过首先将纵波、转换波地震数据从时间域转换到深度域，并在深度域，根据纵波与转换波在相同深度处的地震能量数据体最大相关性规则，对转换波地震数据进行深度层位修正，完成纵波与转换波在深度域的匹配。然后再将深度匹配后的转换波地震数据从深度域转换到纵波时间域，最终实现纵波转换波地震数据的时间匹配。利用本申请各个实施例，可以提高纵波转换波时间匹配的效率以及准确度，缩短纵波转换波地震资料联合解释的周期并提升联合解释的效果。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本说明书提供的一种纵波转换波地震数据时间匹配的方法实施例的流程示意图；
图2为本说明书提供的一个实施例中转换波地震勘探示意图；
图3为本说明书提供的一个实施例中纵波的时间层位剖面以及均方根速度剖面示意图；
图4为本说明书提供的一个实例中纵波以及转换波的原始时间域地震剖面示意图；
图5为本说明书提供的一个实例中纵波的时间层位剖面以及深度层位剖面示意图；
图6为本说明书提供的一个实例中转换波的时间层位剖面以及深度层位剖面示意图；
图7为本说明书提供的一个实例中纵波转换波深度域能量值相关值扫描曲线示意图；
图8为本说明书提供的一个实例中时间匹配后的纵波以及转换波的时间域地震剖面示意图；
图9为本说明书提供的一种纵波转换波地震数据时间匹配的装置实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本说明书中的技术方案，下面将结合本说明书一个或多个实施例中的附图，对本说明书一个或多个实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是说明书一部分实施例，而不是全部的实施例。基于说明书一个或多个实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本说明书实施例方案保护的范围。
纵波与转换波地震剖面的时间匹配问题是由于转换波传播路径的非对称性(参考图2)引起的，这导致利用纵波与转换波的地震数据联合进行油气藏解释困难。由于地震纵波与地震转换波的传播速度不同，造成同一地下地质层位在纵波地震资料和转换波地震资料上对应的传播时间相差很大。因此，利用纵波与转换波资料进行联合反演与解释，需要把它们校正到同一时间域(即相同地下地层层位在纵波与转换波地震资料上的传播时间一致)。而如何将纵波地震剖面与转换波地震剖面进行时间匹配，一致是困扰油气勘探的技术难题。
本申请通过首先将纵波、转换波地震数据从时间域转换到深度域，并在深度域，根据同一地下层位在纵波与转换波的深度域地震剖面上能量值相似性，对转换波地震数据进行深度层位校正，完成纵波与转换波在深度域的层位匹配。然后再将深度域层位匹配后的转换波地震数据从深度域转换到纵波时间域，最终实现纵波转换波地震数据的时间匹配。与传统在时间域进行纵波转换波地震数据的时间匹配的方法技术相比，本申请提出的方法可以在降低纵波转换波时间匹配难度的同时，进一步提高匹配效率，最终缩短纵波转换波地震资料联合解释的周期并提升效果。
本申请各个实施例中，所述转换波可以指地震纵波在地下地层界面处发生反射或者透射时所产生的地震横波。转换波地震勘探技术是在油气勘探的目标区域的地表利用炸药或者机械装置激发下行传播入地的地震纵波(简称P波)，并利用置于地表的检波器在一定的时间范围内(通常不大于10秒)按照规定的时间间隔记录地下地层反射的上行传播的地震横波(简称S波)在地表引发的质点振动形态，并将质点振动形态的量化值按照离散时间顺序记录下来并保存在磁盘介质中，供油气勘探技术人员在后期利用计算机软件对记录下来的转换波地震数据进行处理与解释。
图1是本说明书提供的所述一种纵波转换波地震数据时间匹配的方法实施例流程示意图。虽然本说明书提供了如下述实施例或附图所示的方法操作步骤或装置结构，但基于常规或者无需创造性的劳动在所述方法或装置中可以包括更多或者部分合并后更少的操作步骤或模块单元。在逻辑性上不存在必要因果关系的步骤或结构中，这些步骤的执行顺序或装置的模块结构不限于本说明书实施例或附图所示的执行顺序或模块结构。所述的方法或模块结构的在实际中的装置、服务器或终端产品应用时，可以按照实施例或者附图所示的方法或模块结构进行顺序执行或者并行执行(例如并行处理器或者多线程处理的环境、甚至包括分布式处理、服务器集群的实施环境)。
具体的一个实施例如图1所示，本说明书提供的纵波转换波地震数据时间匹配的方法的一个实施例中，所述方法可以包括：
S2、对纵波以及转换波的原始时间域地震数据分别进行时间深度转换，获得纵波以及转换波的深度域地震数据。
本实施例中，可以分别获取纵波以及转换波的原始时间域地震数据，然后根据地层层位速度数据或者构建的空间速度场获取时深转换速度数据，之后，根据所述时深转换速度数据对所述原始时间域地震数据进行时间深度转换(简称时深转换)，从而将采集的时间域的纵波以及转换波的地震数据从时间域转换到深度域，获得纵波的深度域地震数据以及转换波的深度域地震数据。
本说明书的一个实施例中，可以首先对纵波的原始时间域地震数据进行层位分析，获得纵波时间层位数据。然后，根据纵波均方根速度对纵波时间层位数据进行时间深度转换，获得纵波深度层位数据，再对所述纵波深度层位数据利用二维内插以及外推算法，获得纵波的深度域地震数据。
首先，可以在纵波的原始时间域地震数据的显示图上，根据对纵波地震剖面上的层位关系进行分析交互解释层位，建立纵波地震剖面的时间层位模型。然后，进一步根据地震数据计算得到相应的均方根速度数据，获得相应的纵波均方根速度剖面。之后，从均方根速度剖面中读取对应时间层位模型上每一个控制点处的纵波均方根速度值其中，t表示时间，x表示地震道，v表示均方根速度，pp表示纵波，i表示层位序号，j表示每一层位对应的控制点序号。如图3所示，图3左表示纵波地震剖面时间层位模型示意图，图3右表示相应的纵波均方根速度剖面示意图。
然后，根据纵波均方根速度对纵波时间层位数据进行时间深度转换，获得纵波深度层位数据，再以此为基础，通过二维内插以及外推算法，确定纵波在深度域的地震数据，从而提高获得的深度域地震数据的准确性。
本说明书的一个实施例中，对于纵波地震时间层位模型可以根据纵波均方根速度计算每一层的等效均方根速度具体如下：
然后，对于任意一控制点j，计算相应的地震纵波的层速度可以为：
利用上述层速度根据下述公式可以进一步计算出每一控制点的深度位置为：
从而，基于纵波均方根速度，利用上述公式(1)、(2)、(3)实现纵波的时间深度转换，将纵波时间层位数据转换到深度域，获得纵波深度层位数据。
本说明书的另一个实施例中，以类似的方法，在转换波的原始时间域地震数据的显示图上，建立相应的转换波时间层位模型，并获取转换波均方根速度数据以及纵横波速度比(即纵波与转换波的速度比值)数据。然后，确定转换波时间层位模型中每一个控制点处的转换波均方根速度值以及纵横波垂直速度比值其中，ps表示转换波。然后，根据转换波均方根速度以及纵横波垂直速度比值数据对转换波时间层位数据进行时间深度转换，获得转换波深度层位数据。再对所述转换波深度层位数据利用二维内插以及外推算法更加准确的确定出转换波的深度域地震数据。
具体的，对于转换波地震时间层位模型，转换波等效均方根速度可以表示为：
以及等效纵横波垂直速度比值：
则，对于任意一控制点j，转换波相应的层速度可以通过如下公式计算确定：
相应的，转换波每一控制点对应的深度位置可以表示为：
从而，利用上述公式(3)以及(7)可以实现纵波以及转换波时间层位模型的时深转换。然后，以此为基础，通过二维内插以及外推算法，计算出纵波以及转换波地震数据中每一采样点所对应的深度位置，从而将时间域地震数据转换到深度域，获得纵波以及转换波的深度域地震数据。分别记为和其中，表示纵波的深度域地震数据，表示转换波的深度域地震数据，m为采样点序号，l为地震道序号。
S4、基于预设的纵波与转换波在相同深度处的地震能量数据体最大相关性规则对转换波的深度域地震数据进行深度修正处理。
本实施例中，可以首先根据纵波以及转换波的深度域地震数据，分别计算获得纵波以及转换波在深度域的地震能量数据。然后，基于预设的纵波与转换波在相同深度处的地震能量数据体最大相关性规则对转换波的深度域地震数据进行深度修正处理。其中，所述能量数据体最大相关性规则可以包括：纵波与转换波在相同深度处对应的能量数据体的相关程度最大，或者纵波与转换波在相同深度处对应的能量数据体具有最大相关系数。然后，如可以对比转换波某一深度h1对应的地震能量数据体与纵波各深度处的地震能量数据体之间的相关程度，从而确定出最大相关程度对应的纵波深度值h2。之后，可以根据h2与h1差值对转换波深度h1的深度域地震数据进行深度修正。
例如，本说明书的一个实施例中，可以对纵波的深度域地震数据进行希尔伯特变换得到变换后的数据然后进一步的利用下述公式确定纵波的深度域地震能量数据
同理，可以计算得到转换波的深度域地震能量数据
然后，本说明书的一个实施例中，可以计算转换波深度层位的能量数据体与纵波各深度层位的能量数据体的相关值，根据最大相关值对应的纵波深度层位数据确定所述转换波深度层位的最优深度修正量，利用所述最优深度修正量对转换波的深度域地震数据的深度进行修正。
例如，针对深度域的某一层位可以根据所述地震能量数据体最大相关性规则构建如下目标函数：
其中，wi(z)表示转换波的第i深度层位对应深度修正量，表示转换波某深度层位的能量数据体与纵波深度各层位的能量数据体之间的相关值。
然后，可以在预先给定的范围中自动计算并搜索的最大值，从而可以得到转换波在层位对应的最优深度修正值{wi(z)optimal,i＝1,2,...,I}。依据上述公式(10)，对所有层位分析完毕后，获得转换波在深度域每一层位对应的最优深度修正值。然后，利用所述最优深度修正值对转换波的相应深度层位进行修正处理。相应的，修正后的转换波深度层位可以表示为：
然后，将每个层位的最优深度修正值进行平滑，获得所有层位的最优深度修正值之后，对其进行内插以及外推，得到转换波的深度域地震数据的所有采样点的最优深度修正值。利用每个采样点的最优深度修正值对转换波的深度域地震数据进行深度修正处理，获得深度修正处理后的转换波的深度域地震数据。从而完成纵波与转换波在深度域的匹配
S6、将深度修正处理后的所述转换波的深度域地震数据转换到纵波时间域，获得时间匹配后的转换波时间域地震数据。
本实施例中，可以获取深度修正处理后的转换波深度域地震数据，将深度修正处理后的转换波深度域地震数据的深度代入公式(3)的左边，计算出相应的时间。从而将深度修正处理后的转换波深度域地震数据转换到纵波时间域，获得时间匹配后的转换波时间域地震数据。从而实现了纵波与转换波地震数据在同一时间坐标中显示与分析的目的。
为了使得进一步来验证本方法的实用性和可行性，本说明书还提供了应用上述方案的另一个具体实例。本实例依托中国西南某油田的地震勘探数据进行时间匹配的过程为例进行说明，具体实施方案如下：
(1)资料收集
收集勘探工区中同一条测线的纵波叠前时间偏移数据体与转换波叠前时间偏移数据体，要求两个数据体的面元划分参数保持一致；收集勘探工区中地震纵波均方根速度与转换波均方根速度资料；收集勘探工区中纵横波速度比值资料。图4表示油田同一测线的纵波地震叠前时间偏移数据(左)与转换波地震叠前时间偏移数据(右)的对比显示，可见相同测线的地下地层在纵波数据和转换波数据上有明显的时间差异。
(2)建立时间层位模型
显示该测线的纵波叠前时间偏移数据体，人工交互解释建立时间层位模型；显示该测线的转换波叠前时间偏移数据体，人工交互解释建立时间层位模型。其中，图5所示是按照本说明书提供的方案对图4中的纵波地震数据建立时间层位模型(左图)与时间层位模型的时深转换后的深度层位模型(右图)示意图。图6所示是按照本说明书提供的方案对图4中的转换波地震数据建立时间层位模型(左图)及其时间层位模型的时深转换后的深度层位模型(右图)示意图。
(3)时间深度转换
根据纵波时间层位模型上的控制点的位置读取对应的纵波均方根速度值，根据公式(3)计算纵波时间层位模型上控制点所对应的深度位置，并通过内插与外推的方法计算纵波地震资料所有采样点对应的深度位置，将纵波叠前时间偏移数据体转换到深度域；根据转换波时间层位模型上的控制点位置读取对应的转换波均方根速度值与纵横波速度比值，根据公式(7)计算转换波时间层位模型上控制点所对应的深度位置，并通过内插与外推的方法计算转换波地震资料所有采样点对应的深度位置，将转换波叠前时间偏移数据体转换到深度域。
(4)深度域能量值计算
根据公式(8)和(9)计算纵波与转换波深度域数据体的能量值。其中，图7是图4中的纵波与转换波数据在深度域的能量数据体的某个层位控制点处的相关值扫描曲线，其中虚线圆圈中的圆点为相关值最大位置，圆点所对应的横坐标的读数为层位深度修正值。
(5)深度域纵波转换波层位匹配
根据公式(10)逐层计算纵波与转换波的能量值数据体的相关值，并将最大相关值对应的深度修正量记录下来，并对每个层位的修正量进行平滑。得到所有层位的深度修正量之后，对其进行内插外推，得到转换波深度域数据体的所有采样点的深度修正量，并对转换波深度域数据体进行层位深度修正。
(6)匹配后的深度域转换波数据体转换到时间域
根据公式(3)计算深度域转换波数据体对应的时间值，将深度域转换波数据体转换到时间域，得到了时间匹配后的转换波地震数据。其中，图8是该油田同一条测线的纵波地震叠前时间偏移数据(左)与匹配后的转换波地震数据(右)显示，对比可见，左右两图中的主要反射层的层位时间基本一致。
从而，利用本说明书上述一个或者多个实施例提供的方案，相对于传统在时间域进行纵波转换波时间匹配的方案，效率及准确性更高。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。具体的可以参照前述相关处理相关实施例的描述，在此不做一一赘述。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
本说明书一个或多个实施例提供的一种纵波转换波地震数据时间匹配的方法，可以通过首先将纵波、转换波地震数据从时间域转换到深度域，并在深度域，根据纵波与转换波在相同深度处的地震能量数据体最大相关性规则，对转换波地震数据进行深度层位修正，完成纵波与转换波在深度域的匹配。然后再将深度匹配后的转换波地震数据从深度域转换到纵波时间域，最终实现纵波转换波地震数据的时间匹配。利用本申请各个实施例，可以提高纵波转换波时间匹配的效率以及准确度，缩短纵波转换波地震资料联合解释的周期并提升联合解释的效果。
基于上述所述的纵波转换波地震数据时间匹配的方法，本说明书一个或多个实施例还提供一种纵波转换波地震数据时间匹配的装置。所述的装置可以包括使用了本说明书实施例所述方法的系统、软件(应用)、模块、组件、服务器等并结合必要的实施硬件的装置。基于同一创新构思，本说明书实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似，因此本说明书实施例具体的装置的实施可以参见前述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。具体的，图9说明书提供的一种纵波转换波地震数据时间匹配的装置实施例的模块结构示意图，如图9，所述装置可以包括：
时深转换模块102，可以用于对纵波以及转换波的原始时间域地震数据分别进行时间深度转换，获得纵波以及转换波的深度域地震数据；
能量计算模块104，可以用于根据所述纵波以及转换波的深度域地震数据分别计算纵波以及转换波在深度域的地震能量数据；
深度修正模块106，可以用于基于预设的纵波与转换波在相同深度处的地震能量数据体最大相关性规则对转换波的深度域地震数据进行深度修正处理；
时间匹配模块108，可以用于将深度修正处理后的所述转换波的深度域地震数据转换到纵波时间域，获得时间匹配后的转换波时间域地震数据。
需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
本说明书一个或多个实施例提供的一种纵波转换波地震数据时间匹配的装置，可以通过首先将纵波、转换波地震数据从时间域转换到深度域，并在深度域，根据纵波与转换波在相同深度处的地震能量数据体最大相关性规则，对转换波地震数据进行深度层位修正，完成纵波与转换波在深度域的匹配。然后再将深度匹配后的转换波地震数据从深度域转换到纵波时间域，最终实现纵波转换波地震数据的时间匹配。利用本申请各个实施例，可以提高纵波转换波时间匹配的效率以及准确度，缩短纵波转换波地震资料联合解释的周期并提升联合解释的效果。
本说明书提供的上述实施例所述的方法或装置可以通过计算机程序实现业务逻辑并记录在存储介质上，所述的存储介质可以计算机读取并执行，实现本说明书实施例所描述方案的效果。因此，本说明书还提供一种纵波转换波地震数据时间匹配的装置，包括处理器及存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：
对纵波以及转换波的原始时间域地震数据分别进行时间深度转换，获得纵波以及转换波的深度域地震数据；
根据所述纵波以及转换波的深度域地震数据分别计算纵波以及转换波在深度域的地震能量数据；
基于预设的纵波与转换波在相同深度处的地震能量数据体最大相关性规则对转换波的深度域地震数据进行深度修正处理；
将深度修正处理后的所述转换波的深度域地震数据转换到纵波时间域，获得时间匹配后的转换波时间域地震数据。
所述存储介质可以包括用于存储信息的物理装置，通常是将信息数字化后再以利用电、磁或者光学等方式的媒体加以存储。所述存储介质有可以包括：利用电能方式存储信息的装置如，各式存储器，如RAM、ROM等；利用磁能方式存储信息的装置如，硬盘、软盘、磁带、磁芯存储器、磁泡存储器、U盘；利用光学方式存储信息的装置如，CD或DVD。当然，还有其他方式的可读存储介质，例如量子存储器、石墨烯存储器等等。
需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
上述实施例所述的一种纵波转换波地震数据时间匹配的装置，可以通过首先将纵波、转换波地震数据从时间域转换到深度域，并在深度域，根据纵波与转换波在相同深度处的地震能量数据体最大相关性规则，对转换波地震数据进行深度层位修正，完成纵波与转换波在深度域的匹配。然后再将深度匹配后的转换波地震数据从深度域转换到纵波时间域，最终实现纵波转换波地震数据的时间匹配。利用本申请各个实施例，可以提高纵波转换波时间匹配的效率以及准确度，缩短纵波转换波地震资料联合解释的周期并提升联合解释的效果。
本说明书还提供一种纵波转换波地震数据时间匹配的系统，所述系统可以为单独的纵波转换波地震数据时间匹配的系统，也可以应用在多种类型的地震勘探采集或者处理、解释系统中。所述的系统可以为单独的服务器，也可以包括使用了本说明书的一个或多个所述方法或一个或多个实施例装置的服务器集群、系统(包括分布式系统)、软件(应用)、实际操作装置、逻辑门电路装置、量子计算机等并结合必要的实施硬件的终端装置。所述纵波转换波地震数据时间匹配的系统可以包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现上述任意一个或者多个实施例中所述方法的步骤。
需要说明的，上述所述的系统根据方法或者装置实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
上述实施例所述的一种纵波转换波地震数据时间匹配的系统，可以通过首先将纵波、转换波地震数据从时间域转换到深度域，并在深度域，根据纵波与转换波在相同深度处的地震能量数据体最大相关性规则，对转换波地震数据进行深度层位修正，完成纵波与转换波在深度域的匹配。然后再将深度匹配后的转换波地震数据从深度域转换到纵波时间域，最终实现纵波转换波地震数据的时间匹配。利用本申请各个实施例，可以提高纵波转换波时间匹配的效率以及准确度，缩短纵波转换波地震资料联合解释的周期并提升联合解释的效果。
需要说明的是，本说明书上述所述的装置或者系统根据相关方法实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照方法实施例的描述，在此不作一一赘述。本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于硬件+程序类、存储介质+程序实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
上述实施例阐明的系统、装置、模块或单元，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。一种典型的实现设备为计算机。具体的，计算机例如可以为个人计算机、膝上型计算机、车载人机交互设备、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本说明书一个或多个时可以把各模块的功能在同一个或多个软件和/或硬件中实现，也可以将实现同一功能的模块由多个子模块或子单元的组合实现等。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法或者设备中还存在另外的相同要素。
本领域技术人员应明白，本说明书一个或多个实施例可提供为方法、系统或计算机程序产品。因此，本说明书一个或多个实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书一个或多个实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书一个或多个实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本本说明书一个或多个实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述并不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
Claims (10)
1.一种纵波转换波地震数据时间匹配的方法，其特征在于，包括：
对纵波以及转换波的原始时间域地震数据分别进行时间深度转换，获得纵波以及转换波的深度域地震数据；
根据所述纵波以及转换波的深度域地震数据分别计算纵波以及转换波在深度域的地震能量数据；
基于预设的纵波与转换波在相同深度处的地震能量数据体最大相关性规则对转换波的深度域地震数据进行深度修正处理；
将深度修正处理后的所述转换波的深度域地震数据转换到纵波时间域，获得时间匹配后的转换波时间域地震数据。
2.根据权利要求1所述的纵波转换波地震数据时间匹配的方法，其特征在于，所述对纵波的原始时间域地震数据进行时间深度转换，包括：
对纵波的原始时间域地震数据进行层位分析，获得纵波时间层位数据；
根据纵波均方根速度对纵波时间层位数据进行时间深度转换，获得纵波深度层位数据；
对所述纵波深度层位数据利用二维内插以及外推算法确定纵波的深度域地震数据。
3.根据权利要求1或2所述的纵波转换波地震数据时间匹配的方法，其特征在于，所述对转换波的原始时间域地震数据进行时间深度转换，包括：
对转换波的原始时间域地震数据进行层位分析，获得转换波时间层位数据；
根据转换波均方根速度以及纵横波垂直速度比值数据对转换波时间层位数据进行时间深度转换，获得转换波深度层位数据；
对所述转换波深度层位数据利用二维内插以及外推算法确定转换波的深度域地震数据。
4.根据权利要求2所述的纵波转换波地震数据时间匹配的方法，其特征在于，所述根据纵波均方根速度对纵波时间层位数据进行时间深度转换，包括：
根据如下公式对纵波时间层位数据进行时间深度转换：
其中，表示纵波深度层位数据，表示纵波均方根速度，表示纵波时间层位数据，i表示时间层位序号，j表示每一时间层位的控制点序号。
5.根据权利要求3所述的纵波转换波地震数据时间匹配的方法，其特征在于，所述根据转换波均方根速度以及纵横波垂直速度比值数据对转换波时间层位数据进行时间深度转换，包括：
根据如下公式对转换波时间层位数据进行时间深度转换：
其中，表示纵波深度层位数据，表示纵波均方根速度，表示纵波时间层位数据，表示纵横波垂直速度比值，i表示时间层位序号，j表示每一时间层位的控制点序号。
6.根据权利要求1所述的纵波转换波地震数据时间匹配的方法，其特征在于，所述分别计算纵波以及转换波在深度域的地震能量数据，包括：
对所述纵波以及转换波的深度域地震数据分别进行希尔伯特变换；
根据纵波以及转换波希尔伯特变换后的深度域地震数据确定纵波以及转换波在深度域的地震能量数据。
7.根据权利要求1所述的纵波转换波地震数据时间匹配的方法，其特征在于，所述对转换波的深度域地震数据进行深度修正处理，包括：
计算转换波深度层位的能量数据体与纵波各深度层位的能量数据体的相关值；
根据最大相关值对应的纵波深度层位数据确定所述转换波深度层位的最优深度修正值，利用所述最优深度修正值对转换波的深度域地震数据的深度进行修正。
8.一种纵波转换波地震数据时间匹配的装置，其特征在于，包括：
时深转换模块，用于对纵波以及转换波的原始时间域地震数据分别进行时间深度转换，获得纵波以及转换波的深度域地震数据；
能量计算模块，用于根据所述纵波以及转换波的深度域地震数据分别计算纵波以及转换波在深度域的地震能量数据；
深度修正模块，用于基于预设的纵波与转换波在相同深度处的地震能量数据体最大相关性规则对转换波的深度域地震数据进行深度修正处理；
时间匹配模块，用于将深度修正处理后的所述转换波的深度域地震数据转换到纵波时间域，获得时间匹配后的转换波时间域地震数据。
9.一种纵波转换波地震数据时间匹配的装置，其特征在于，包括处理器及用于存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括以下步骤：
对纵波以及转换波的原始时间域地震数据分别进行时间深度转换，获得纵波以及转换波的深度域地震数据；
根据所述纵波以及转换波的深度域地震数据分别计算纵波以及转换波在深度域的地震能量数据；
基于预设的纵波与转换波在相同深度处的地震能量数据体最大相关性规则对转换波的深度域地震数据进行深度修正处理；
将深度修正处理后的所述转换波的深度域地震数据转换到纵波时间域，获得时间匹配后的转换波时间域地震数据。
10.一种纵波转换波地震数据时间匹配的系统，其特征在于，包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现权利要求1-7中任意一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810471934.1A CN108508481B (zh) | 2018-05-17 | 2018-05-17 | 一种纵波转换波地震数据时间匹配的方法、装置及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810471934.1A CN108508481B (zh) | 2018-05-17 | 2018-05-17 | 一种纵波转换波地震数据时间匹配的方法、装置及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108508481A true CN108508481A (zh) | 2018-09-07 |
| CN108508481B CN108508481B (zh) | 2019-11-08 |
Family
ID=63401033
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810471934.1A Active CN108508481B (zh) | 2018-05-17 | 2018-05-17 | 一种纵波转换波地震数据时间匹配的方法、装置及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108508481B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111665558A (zh) * | 2019-03-07 | 2020-09-15 | 中国石油天然气集团有限公司 | 一种油气分布检测方法和系统 |
| CN112305607A (zh) * | 2020-10-28 | 2021-02-02 | 中国石油天然气集团有限公司 | 地震数据纵剖面图形绘制方法及装置 |
| CN112362148A (zh) * | 2020-09-21 | 2021-02-12 | 中铁第四勘察设计院集团有限公司 | 基于斜拉索索力影响面加载动态识别车辆轴重方法和系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030021184A1 (en) * | 2001-06-08 | 2003-01-30 | Yaohui Zhang | Extraction of P-wave and S-wave velocities from multi-component seismic data by joint velocity inversion |
| CN103091709A (zh) * | 2012-12-25 | 2013-05-08 | 中国石油天然气集团公司 | 获得纵波、转换波地震数据时间匹配关系的方法和装置 |
| CN104570079A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 一种纵波、转换横波地震资料的时间匹配方法 |
| CN106226811A (zh) * | 2016-07-05 | 2016-12-14 | 中国石油天然气集团公司 | 一种基于属性的纵波数据和转换波数据匹配的方法及装置 |
-
2018
- 2018-05-17 CN CN201810471934.1A patent/CN108508481B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030021184A1 (en) * | 2001-06-08 | 2003-01-30 | Yaohui Zhang | Extraction of P-wave and S-wave velocities from multi-component seismic data by joint velocity inversion |
| CN103091709A (zh) * | 2012-12-25 | 2013-05-08 | 中国石油天然气集团公司 | 获得纵波、转换波地震数据时间匹配关系的方法和装置 |
| CN104570079A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 一种纵波、转换横波地震资料的时间匹配方法 |
| CN106226811A (zh) * | 2016-07-05 | 2016-12-14 | 中国石油天然气集团公司 | 一种基于属性的纵波数据和转换波数据匹配的方法及装置 |
Non-Patent Citations (2)
| Title |
|---|
| 饶静: "《基于有理插值的纵横波匹配方法应用研究》", 《成都理工大学硕士学位论文》 * |
| 龚雪萍等: "《改进纵波与转换波时间匹配方法研究》", 《石油地球物理勘探》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111665558A (zh) * | 2019-03-07 | 2020-09-15 | 中国石油天然气集团有限公司 | 一种油气分布检测方法和系统 |
| CN111665558B (zh) * | 2019-03-07 | 2023-02-28 | 中国石油天然气集团有限公司 | 一种油气分布检测方法和系统 |
| CN112362148A (zh) * | 2020-09-21 | 2021-02-12 | 中铁第四勘察设计院集团有限公司 | 基于斜拉索索力影响面加载动态识别车辆轴重方法和系统 |
| CN112362148B (zh) * | 2020-09-21 | 2022-01-14 | 中铁第四勘察设计院集团有限公司 | 基于斜拉索索力影响面加载动态识别车辆轴重方法和系统 |
| CN112305607A (zh) * | 2020-10-28 | 2021-02-02 | 中国石油天然气集团有限公司 | 地震数据纵剖面图形绘制方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108508481B (zh) | 2019-11-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107817535B (zh) | 薄夹层的确定方法和装置 | |
| US10995592B2 (en) | Method and system for analyzing the uncertainty of subsurface model | |
| CN107976713B (zh) | 一种高维地震数据输入下去除沉积背景的方法及装置 | |
| CN106772579B (zh) | 一种薄煤层中地震叠前反演方法和装置 | |
| CN107783185B (zh) | 一种层析静校正的处理方法及装置 | |
| CN110031896A (zh) | 基于多点地质统计学先验信息的地震随机反演方法及装置 | |
| CN108508481B (zh) | 一种纵波转换波地震数据时间匹配的方法、装置及系统 | |
| CN109870719A (zh) | 一种碳酸盐岩致密薄储层分布确定方法、装置及系统 | |
| CN109001813A (zh) | 一种压制多次波的方法、装置及系统 | |
| CN109143358B (zh) | 一种获得深层、超深层碎屑岩地层压力结构的方法及装置 | |
| CN106707344B (zh) | 一种地层层序划分方法及装置 | |
| CN109188520A (zh) | 薄储层厚度预测方法及装置 | |
| CN104199092A (zh) | 基于多层次框架的三维全层位自动追踪方法 | |
| Lu et al. | Improving the image quality of elastic reverse-time migration in the dip-angle domain using deep learning | |
| CN110837115A (zh) | 陆相混积岩致密储层岩性的地震识别方法及装置 | |
| Wei et al. | Small-data-driven fast seismic simulations for complex media using physics-informed Fourier neural operators | |
| Artun et al. | Intelligent seismic inversion workflow for high-resolution reservoir characterization | |
| RU2490677C2 (ru) | Способ комплексной обработки геофизических данных и технологическая система "литоскан" для его осуществления | |
| CN106338764A (zh) | 生物启发计算地层圈闭油气藏超剥线识别方法 | |
| CN103543478A (zh) | 地质形态插值的km方法 | |
| CN106842316B (zh) | 裂缝确定方法和装置 | |
| Maity et al. | Reservoir characterization of an unconventional reservoir by integrating microseismic, seismic, and well log data | |
| CN108761533B (zh) | 一种确定纵横波速度比的方法、装置及系统 | |
| Jika et al. | Application of geostatistics in facies modeling of Reservoir-E,“Hatch Field” offshore Niger Delta Basin, Nigeria | |
| CN108492014A (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 |