CN105353412B - 一种井震联合平均速度场的计算方法及系统 - Google Patents
一种井震联合平均速度场的计算方法及系统 Download PDFInfo
- Publication number
- CN105353412B CN105353412B CN201510933562.6A CN201510933562A CN105353412B CN 105353412 B CN105353412 B CN 105353412B CN 201510933562 A CN201510933562 A CN 201510933562A CN 105353412 B CN105353412 B CN 105353412B
- Authority
- CN
- China
- Prior art keywords
- seismic
- average
- logging
- well
- velocity
- 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
Links
- 238000000205 computational method Methods 0.000 title abstract 3
- 208000035126 Facies Diseases 0.000 claims abstract description 168
- 238000000034 method Methods 0.000 claims abstract description 61
- 238000005553 drilling Methods 0.000 claims abstract description 19
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
- 238000004364 calculation method Methods 0.000 claims description 21
- 230000008021 deposition Effects 0.000 claims description 10
- 239000002131 composite material Substances 0.000 claims description 6
- 238000007621 cluster analysis Methods 0.000 claims description 2
- 238000005192 partition Methods 0.000 claims 1
- 230000006870 function Effects 0.000 description 80
- 238000006243 chemical reaction Methods 0.000 description 12
- 230000008859 change Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 8
- 238000010276 construction Methods 0.000 description 8
- 230000008901 benefit Effects 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000013213 extrapolation Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000002948 stochastic simulation Methods 0.000 description 1
- 238000012916 structural analysis Methods 0.000 description 1
- 230000009466 transformation 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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- 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
-
- 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/66—Subsurface modeling
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种井震联合平均速度场的计算方法及计算系统，包括：利用目的层的测井相类型及目的层的地震相类型，划分目的层的沉积相，并对所述沉积相进行编码，得到沉积相编码；利用地震叠加速度场，通过Dix公式计算得到地震平均速度；利用所述测井资料，计算得到目的层在各个钻井位置处的测井平均速度；将所述沉积相编码作为约束项，所述地震平均速度作为次级变量，所述测井平均速度作为主变量，利用协克里金估计算法进行协克里金插值，得到井震联合平均速度场；该方法利用沉积相约束下的高精度井震联合平均速度场计算方法，确保速度场的横向变化趋势与实际地质特征相吻合，从而提高构造图时深转换的精度。
Description
技术领域
本发明涉及石油地球物理勘探领域，特别涉及一种井震联合平均速度场的计算方法及系统。
背景技术
随着油气勘探程度的不断提高，勘探区域的地质条件也变得越来越复杂，地下介质非均质性强、速度纵横向变化大等问题对高精度速度场建模技术提出了挑战，同时也严重影响着构造图等地震解释成果图件的时深转换精度。在构造图时深转换中通常采用平均速度场，而目前大部分平均速度场建模方法都是以地震速度转换或测井速度内插外推为主。其中，以地震为主的平均速度场建模方法主要有叠加速度分析方法、层析反演方法等。常规叠加速度分析法具有较高的计算效率，但不适用于界面倾斜和速度横向变化的情况，通过与射线追踪方法相结合能够在一定程度上解决层速度横向变化的问题，但依然很难解决复杂构造区域的速度场建模问题。
而层析速度反演法则以观测数据与合成数据的误差最小为原则，反演出地下介质的结构与速度分布等重要信息，该方法的精度高于叠加速度分析所建立的平均速度场，但要求具有较高的资料信噪比和可靠的初始速度模型。基于地震资料的平均速度场建模方法在横向上具有较好的连续性，但纵向精度不高，速度场的准确性和可靠性难以评价，且该速度场在应用中往往缺乏明确的地质意义。
以测井资料为主的平均速度场建模方法主要包括克里金估计法、随机模拟以及随机反演等方法。克立金估计通过开展变差函数理论和结构分析，在有限区域内对区域化变量的取值进行最优无偏估计，该方法能够有效反映变量的空间结构性，而且能给出速度的估计精度，因此在实际中得到了广泛的应用。但克里金估计是一种局部估计方法，对估计值的整体空间相关性考虑不够，加上该方法属于光滑内插，因此很容易导致一些有意义的异常带被“光滑”掉。随机模拟与随机反演法则是为了克服克里金估计的光滑效应而产生的，在实际中也有较为广泛的应用。显然，以测井资料为主的速度场建模方法具有纵向精度高的优势，但建立的速度场横向连续性往往不够可靠，并且在钻井资料较少时作为“硬数据”的测井数据无法提供充足的先验信息，特别是在构造复杂、横向变化快的河流相储层中，直接利用测井资料进行内插外推建立速度场很容易产生构造假象。因此，亟需发展一种能够将地震横向连续性与测井纵向高分辨率特征有效结合的方法。
发明内容
本发明的目的是提供一种井震联合平均速度场的计算方法及计算系统，利用沉积相约束下的高精度井震联合平均速度场计算方法，确保速度场的横向变化趋势与实际地质特征相吻合，从而提高构造图时深转换的精度。
为解决上述技术问题，本发明提供一种井震联合平均速度场的计算方法，包括：
利用目的层的测井相类型及目的层的地震相类型，划分目的层的沉积相，并对所述沉积相进行编码，得到沉积相编码；
利用地震叠加速度场，通过Dix公式计算得到地震平均速度；
利用所述测井资料，计算得到目的层在各个钻井位置处的测井平均速度；
将所述沉积相编码作为约束项，所述地震平均速度作为次级变量，所述测井平均速度作为主变量，利用协克里金估计算法进行协克里金插值，得到井震联合平均速度场。
其中，对所述沉积相进行编码，得到沉积相编码包括：
根据划分后的沉积相的空间范围的大小，按照从大到小的顺序，从整数1开始赋予沉积相对应的数字标识，得到沉积相编码。
其中，所述将所述沉积相编码作为约束项，所述地震平均速度作为次级变量，所述测井平均速度作为主变量，利用协克里金估计算法进行协克里金插值，得到井震联合平均速度场，包括：
根据所述地震平均速度，所述测井平均速度以及沉积相编码，建立工区所对应的变差函数；
将每个沉积相编码作为约束项，并利用工区所对应的变差函数，根据协克里金估计的无偏性及最优性原则，计算得到各个待估点搜索范围内的各地震平均速度对应的协克里金加权系数值，每个钻井位置处测井平均速度对应的协克里金加权系数值，以及各沉积相编码对应的协克里金加权系数值，进行协克里金估计，得到井震联合平均速度场。
其中，根据所述地震平均速度，所述测井平均速度以及沉积相编码，建立工区所对应的变差函数，包括：
根据所述地震平均速度，计算所述地震平均速度的离散变差函数，根据所述地震平均速度的离散变差函数选择对应的理论变差函数模型进行拟合，建立地震平均速度的理论变差函数；
根据所述测井平均速度，计算所述测井平均速度的离散变差函数，根据所述测井平均速度的离散变差函数选择对应的理论变差函数模型进行拟合，建立井点处平均速度的理论变差函数；
根据所述地震平均速度及所述测井平均速度，计算所述地震平均速度及所述测井平均速度的离散变差函数，根据所述地震平均速度及所述测井平均速度的离散变差函数选择对应的理论变差函数模型进行拟合，建立井点平均速度与地震平均速度的理论互变差函数；
根据所述沉积相编码，计算所述沉积相编码的离散变差函数，根据所述沉积相编码的离散变差函数选择对应的理论变差函数模型进行拟合，建立沉积相编码的理论变差函数。
其中，所述将每个沉积相编码作为约束项，并利用工区所对应的变差函数，根据协克里金估计的无偏性及最优性原则，计算得到各个待估点搜索范围内的各地震平均速度对应的协克里金加权系数值，每个钻井位置处测井平均速度对应的协克里金加权系数值，以及各沉积相编码对应的协克里金加权系数值，进行协克里金估计，得到井震联合平均速度场，包括：
利用进行协克里金插值，得到井震联合平均速度场；
其中，X*(u0)是u0位置的速度估计值；是在位置上的测井平均速度值，是位置上的测井平均速度值对应的加权系数值；是在位置上的地震平均速度采样值，是位置上的地震平均速度采样值的加权系数值，是在位置上的沉积相编码，是位置上的沉积相编码的加权系数值，N为沉积相的数量。
其中，所述利用目的层的测井相类型及目的层的地震相类型，划分目的层的沉积相，包括：
利用测井资料与地震层位特征进行精细井震标定，并根据精细井震标定结果进行层位解释；
根据层位解释后的地震层位，提取有效地震属性集；
通过对所述有效地震属性集进行聚类分析，确定目的层的地震相类型；
以所述目的层的测井相类型作为约束，利用所述目的层的地震相类型，划分所述目的层的沉积相。
其中，还包括：
利用所述井震联合平均速度场将时间域构造图转换成对应的深度域构造图。
本发明提供一种井震联合平均速度场的计算系统，包括：
沉积相模块，用于利用目的层的测井相类型及目的层的地震相类型，划分目的层的沉积相，并对所述沉积相进行编码，得到沉积相编码；
地震平均速度面模块，用于利用地震叠加速度场，通过Dix公式计算得到地震平均速度；
测井平均速度模块，用于利用所述测井资料，计算得到目的层在各个钻井位置处的测井平均速度；
井震联合平均速度场模块，用于将所述沉积相编码作为约束项，所述地震平均速度作为次级变量，所述测井平均速度作为主变量，利用协克里金估计算法进行协克里金插值，得到井震联合平均速度场。
其中，所述井震联合平均速度场模块包括：
变差函数建立单元，用于根据所述地震平均速度，所述测井平均速度以及沉积相编码，建立工区所对应的变差函数；
协克里金估计单元，用于将每个沉积相编码作为约束项，并利用工区所对应的变差函数，根据协克里金估计的无偏性及最优性原则，计算得到各个待估点搜索范围内的各地震平均速度对应的协克里金加权系数值，每个钻井位置处测井平均速度对应的协克里金加权系数值，以及各沉积相编码对应的协克里金加权系数值，进行协克里金估计，得到井震联合平均速度场。
其中，所述协克里金估计单元具体用于：
利用进行协克里金插值，得到井震联合平均速度场；
其中，X*(u0)是在u0位置的速度估计值；是在位置上的测井平均速度值，是位置上的测井平均速度值对应的加权系数值；是在位置上的地震平均速度采样值，是位置上的地震平均速度采样值的加权系数值，是在位置上的沉积相编码，是位置上的沉积相编码的加权系数值，N为沉积相的数量。
本发明所提供的井震联合平均速度场的计算方法及计算系统，包括：利用目的层的测井相类型及目的层的地震相类型，划分目的层的沉积相，并对所述沉积相进行编码，得到沉积相编码；利用地震叠加速度场，通过Dix公式计算得到地震平均速度；利用所述测井资料，计算得到目的层在各个钻井位置处的测井平均速度；将所述沉积相编码作为约束项，所述地震平均速度作为次级变量，所述测井平均速度作为主变量，利用协克里金估计算法进行协克里金插值，得到井震联合平均速度场；
利用沉积相特征来更为合理地控制速度场横向变化特征，以测井速度为主变量，沉积相和地震平均速度为次级变量，采用协克里金估计算法将测井速度的高纵向分辨率特征与沉积相和地震速度的横向连续性进行有机结合；从而解决了现有技术中单纯以地震资料为主的平均速度场计算方法的纵向精度较低，而以测井资料为主的平均速度场计算方法往往难以合理控制速度面的横向趋势变化；该方法利用沉积相约束下的高精度井震联合平均速度场计算方法，确保速度场的横向变化趋势与实际地质特征相吻合，从而提高构造图时深转换的精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
图1为本发明实施例所提供的井震联合平均速度场的计算方法的流程图；
图2为本发明实施例所提供的沉积相获取的方法的流程图；
图3为本发明实施例所提供的地震平均速度的理论变差函数的示意图；
图4为本发明实施例所提供的井点处平均速度的理论变差函数的示意图；
图5为本发明实施例所提供的井点平均速度与地震平均速度的理论互变差函数的示意图；
图6为本发明实施例所提供的井震联合平均速度场的计算系统的结构框图。
具体实施方式
本发明的核心是提供一种井震联合平均速度场的计算方法及计算方法，利用沉积相约束下的高精度井震联合平均速度场计算方法，确保速度场的横向变化趋势与实际地质特征相吻合，从而提高构造图时深转换的精度。
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
请参考图1，图1为本发明实施例所提供的井震联合平均速度场的计算方法的流程图；该方法可以包括：
s100、利用目的层的测井相类型及目的层的地震相类型，划分目的层的沉积相，并对所述沉积相进行编码，得到沉积相编码；
其中，获取沉积相编码的具体过程如下：请参考图2，优选的，所述利用目的层的测井相类型及目的层的地震相类型，划分目的层的沉积相可以包括：
s200、利用测井资料与地震层位特征进行精细井震标定，并根据精细井震标定结果进行层位解释；
其中，精细井震标定是地震资料层位解释的基本依据，可以通过对测井资料进行环境校正并去除野值，利用声波测井速度将曲线从深度域转换到时间域，逐层计算出反射系数，从地震剖面中提取地震子波并与反射系数进行褶积得到合成地震记录，根据测井分层与地震层位特征开展精细时深关系调整，确保井旁地震道与合成地震记录之间达到最佳匹配。基于整个地震工区所有井位的精细井震标定，可以建立起地质目标解释的“种子点集”，然后由点到线、由线到面实现目标地震层位的精细解释，并绘制出目的层等T0构造图。
s210、根据层位解释后的地震层位，提取有效地震属性集；
其中，地震属性是地震相划分的基础，基于上述步骤解释的精细地震层位沿层选择合适的时窗，逐一提取瞬时振幅、均方根振幅、瞬时频率等各种地震属性信息，并对提取出的所有地震属性做K-L变换进行降维映射，获得有效的主成分分量，并进一步结合目的层位置处的地质信息精确计算地震属性与测井数据之间的相关性，根据相关性的大小获得与目的层空间变化特征密切相关的有效地震属性作为地震属性集。
s220、通过对所述有效地震属性集进行聚类分析，确定目的层的地震相类型；
其中，结合上述过程确定的有效地震属性集，采用神经网络开展属性聚类分析例如采用Kohonen神经网络进行属性聚类分析，参考全区整体地质特征和地震资料波形特征确定分类数、精确度等最佳基本聚类参数，并将地震属性聚类结果与原始地震资料波形特征开展对比分析，根据地震资料的外部形态、内部反射结构、顶底接触关系以及地震资料的运动学与动力学特征，确定地震相类型及其代表的地质意义，为下一步沉积相的划分做好准备。
s230、以所述目的层的测井相类型作为约束，利用所述目的层的地震相类型，划分所述目的层的沉积相。
其中，以目的层测井相类别作为约束，在地震相的基础上结合工区地质认识划分目的层沉积相，并进一步对沉积相进行编码以参与后续计算，得到沉积相编码。可选的，对所述沉积相进行编码，得到沉积相编码可以包括：
根据划分后的沉积相的空间范围的大小，按照从大到小的顺序，从整数1开始赋予沉积相对应的数字标识，得到沉积相编码。即编码时以沉积相的空间范围大小为顺序，分别用整数1、2、3、……、N代表相应的沉积相(其中N为沉积相的数量)，并将这些数字作为沉积相指示变量f引入到后续的多变量数据整合中。
s110、利用地震叠加速度场，通过Dix公式计算得到地震平均速度；
其中，地震平均速度面的计算是基于地震叠加速度的平均速度计算。具体过程可以是：
地震速度是井震联合速度场建模的重要约束条件，该方法充分采用地震资料处理过程中的叠加速度场，通过Dix公式(即公式1)将地震叠加速度体进行转换得到层速度体：
其中：t0,n，t0,n-1分别为基准层到第n层和基准层到第n-1层的地震波旅行时，VR,n，VR,n-1分别为第n层和第n-1层的地震叠加速度，Vn为第n层的层速度。
进一步可以通过公式(2)将层速度体转换成地震平均速度体，由此可以沿着目的层层位提取出地震平均速度面。
其中：vi为第i层的层速度，ti为第i层的旅行时，Vav为地震平均速度面。
s120、利用所述测井资料，计算得到目的层在各个钻井位置处的测井平均速度；
其中，测井资料中的平均速度是井震联合速度场计算的基础数据，基于步骤100井震精细标定中确定的时深关系确定目的层的钻遇深度h，根据目的层等T0构造图提取各钻井井位所对应的旅行时t，并采用公式(3)计算得到目的层在各钻井分层处的平均速度v：
v＝h/t (3)
s130、将所述沉积相编码作为约束项，所述地震平均速度作为次级变量，所述测井平均速度作为主变量，利用协克里金估计算法进行协克里金插值，得到井震联合平均速度场。
其中，沉积相编码作为约束相，使得不同沉积相中计算的速度场更加符合各沉积相的地质特征。即利用沉积相特征来更为合理地控制速度场横向变化特征。以测井速度为主变量，沉积相和地震平均速度为次级变量，采用协克里金估计算法将测井速度的高纵向分辨率特征与沉积相和地震速度的横向连续性进行有机结合的高精度井震联合平均速度场计算方法，确保速度场的横向变化趋势与实际地质特征相吻合，从而提高构造图时深转换的精度。优选的，具体过程可以如下：将所述沉积相编码作为约束项，所述地震平均速度作为次级变量，所述测井平均速度作为主变量，利用协克里金估计算法进行协克里金插值，得到井震联合平均速度场可以包括：
根据所述地震平均速度，所述测井平均速度以及沉积相编码，建立工区所对应的变差函数；
其中，变差函数的建立是采用协克里金方法开展速度场建模的核心步骤，变差函数的精确可以使得后续进行协克里金插值时，更加精确，得到的加权系数值更加合理。具体的，这里的工区所对应的变差函数的具体建立过程可以如下：
根据所述地震平均速度，计算所述地震平均速度的离散变差函数，根据所述地震平均速度的离散变差函数选择对应的理论变差函数模型进行拟合，建立地震平均速度的理论变差函数；请参考图3，通过该变差函数可以确定基台值为1，变程约为8000。
根据所述测井平均速度，计算所述测井平均速度的离散变差函数，根据所述测井平均速度的离散变差函数选择对应的理论变差函数模型进行拟合，建立井点处平均速度的理论变差函数；请参考图4，通过该变差函数可以确定基台值为1.15，变程约为7050。
根据所述地震平均速度及所述测井平均速度，计算所述地震平均速度及所述测井平均速度的离散变差函数，根据所述地震平均速度及所述测井平均速度的离散变差函数选择对应的理论变差函数模型进行拟合，建立井点平均速度与地震平均速度的理论互变差函数；请参考图5，通过该变差函数可以确定基台值为1.05，变程约为6500。
根据所述沉积相编码，计算所述沉积相编码的离散变差函数，根据所述沉积相编码的离散变差函数选择对应的理论变差函数模型进行拟合，建立沉积相编码的理论变差函数。
其中，变差方程的计算如公式(4)所示：
其中，Z(u)为位置u处的值，Z(u+h)为位置u+h处的值。实际上变差方程γ(h)就是指距离为h的任意两点的方差的一半。
其中，计算出的变差方程是些随h变化的离散值，但在对协克里金权系数进行计算时需要的是一个随h值变化的函数关系，因此需要找一个函数关系来对这些离散值进行拟合，也就是理论模型。同时不同的理论模型可以反映工区的多样的结构特征。上述的拟合过程是根据公式(4)计算得到的变差方程计算得到的离散值来选择合适工区目的层速度变化特征的理论模型，其中变差方程的理论模型按变差函数在原点处的性质可分为球状模型、指数模型、高斯模型等选择合适的模型进行拟合得到变差函数。图3-图5三个图中圆点为求得的实验变差函数，可以发现三个实验变差函数均没有块金效应，表明该数据在短距离范围内的平均速度连续性较好。实线为变差函数拟合曲线，可以看出三个图中变差函数可以用球状模型进行较好的拟合，从而为后续的协克里金估计提供了可靠的理论支撑。
将每个沉积相编码作为约束项，并利用工区所对应的变差函数，根据协克里金估计的无偏性及最优性原则，计算得到各个待估点搜索范围内的各地震平均速度对应的协克里金加权系数值，每个钻井位置处测井平均速度对应的协克里金加权系数值，以及各沉积相编码对应的协克里金加权系数值，进行协克里金估计，得到井震联合平均速度场。
其中，协克里金估计算法能够实现测井速度、地震速度和沉积相特征的整合，通过将沉积相指示变量引入到协克里金估计。优选的，具体过程可以如下：
利用进行协克里金插值，得到井震联合平均速度场；
其中，X*(u0)是u0位置的速度估计值；是在位置上的测井平均速度值，是位置上的测井平均速度值对应的加权系数值；是在位置上的地震平均速度采样值，是位置上的地震平均速度采样值的加权系数值，是在位置上的沉积相编码，是位置上的沉积相编码的加权系数值，N为沉积相的数量。
其中，以沉积相作为约束，选取井点平均速度和地震平均速度对应的变差函数，二者的互变差函数以及沉积相编码的变差函数，以井点平均速度作为硬数据，地震平均速度作为软数据进行协克里金插值，得到最终的平均速度面。
基于上述技术方案，本发明实施例提供的井震联合平均速度场的计算方法，利用沉积相特征来更为合理地控制速度场横向变化特征，以测井速度为主变量，沉积相和地震平均速度为次级变量，采用协克里金估计算法将测井速度的高纵向分辨率特征与沉积相和地震速度的横向连续性进行有机结合；从而解决了现有技术中单纯以地震资料为主的平均速度场计算方法的纵向精度较低，而以测井资料为主的平均速度场计算方法往往难以合理控制速度面的横向趋势变化；因此该方法利用沉积相约束下的高精度井震联合平均速度场计算方法，确保速度场的横向变化趋势与实际地质特征相吻合，从而提高构造图时深转换的精度。
即以沉积相和地震速度来约束平均速度场的横向变化，利用沉积相和地震速度与测井速度之间的空间相关性对平均速度场进行空间估计，提高速度场建模的精度和合理性。
基于上述技术方案，该方法还可以包括：
利用所述井震联合平均速度场将时间域构造图转换成对应的深度域构造图。
其中，井震联合平均速度场计算旨在获得具有地质意义的高精度速度场，从而为构造图时深转换提供合理可靠的速度场信息，利用上述计算方法所得到的平均速度场可以将步骤100精细解释所确定的目的层等T0构造图转换成对应的深度域构造图，实现构造图的高精度时深转换。
即具体过程为首先需要开展精细井震标定和层位解释，提取并深入分析地震属性特征，通过聚类分析建立地震相模式并划分沉积相，并进一步采用协克里金估计来整合井点平均速度、地震平均速度和沉积相特征，建立起具有地质意义的高精度平均速度场，实现构造图的精确时深转换。
基于上述技术方案，本发明实施例提供的井震联合平均速度场的计算方法，引入沉积相信息作为约束，采用协克里金估计算法实现了地震速度计算和测井速度计算优势的有机结合，充分发挥了测井速度纵向精度高和地震叠加速度横向连续性好的优势，通过沉积相特征的引入赋予平均速度场更为合理的地质意义，确保了速度面的横向变化趋势更为准确可靠，从而保证了构造图时深转换的精度，为石油地震勘探提供更为准确的构造信息。
本发明实施例提供了井震联合平均速度场的计算方法，通过上述方法能够确保速度场的横向变化趋势与实际地质特征相吻合。
下面对本发明实施例提供的井震联合平均速度场的计算系统进行介绍，下文描述的井震联合平均速度场的计算系统与上文描述的井震联合平均速度场的计算方法可相互对应参照。
请参考图6，图6为本发明实施例所提供的井震联合平均速度场的计算系统的结构框图，该系统可以包括：
沉积相模块100，用于利用目的层的测井相类型及目的层的地震相类型，划分目的层的沉积相，并对所述沉积相进行编码，得到沉积相编码；
地震平均速度面模块200，用于利用地震叠加速度场，通过Dix公式计算得到地震平均速度；
测井平均速度模块300，用于利用所述测井资料，计算得到目的层在各个钻井位置处的测井平均速度；
井震联合平均速度场模块400，用于将所述沉积相编码作为约束项，所述地震平均速度作为次级变量，所述测井平均速度作为主变量，利用协克里金估计算法进行协克里金插值，得到井震联合平均速度场。
可选的，沉积相模块100包括：
标定解释单元，用于利用测井资料与地震层位特征进行精细井震标定，并根据精细井震标定结果进行层位解释；
根地震属性集单元，用于据层位解释后的地震层位，提取有效地震属性集；
聚类单元，用于通过对所述有效地震属性集进行聚类分析，确定目的层的地震相类型；
沉积相单元，用于以所述目的层的测井相类型作为约束，利用所述目的层的地震相类型，划分所述目的层的沉积相；
编码单元，根据划分后的沉积相的空间范围的大小，按照从大到小的顺序，从整数1开始赋予沉积相对应的数字标识，得到沉积相编码。
可选的，所述井震联合平均速度场模块400包括：
变差函数建立单元，用于根据所述地震平均速度，所述测井平均速度以及沉积相编码，建立工区所对应的变差函数；
协克里金估计单元，用于将每个沉积相编码作为约束项，并利用工区所对应的变差函数，根据协克里金估计的无偏性及最优性原则，计算得到各个待估点搜索范围内的各地震平均速度对应的协克里金加权系数值，每个钻井位置处测井平均速度对应的协克里金加权系数值，以及各沉积相编码对应的协克里金加权系数值，进行协克里金估计，得到井震联合平均速度场。
可选的，所述变差函数建立单元具体包括：
第一建立子单元，用于根据所述地震平均速度，计算所述地震平均速度的离散变差函数，根据所述地震平均速度的离散变差函数选择对应的理论变差函数模型进行拟合，建立地震平均速度的理论变差函数；
第二建立子单元，用于根据所述测井平均速度，计算所述测井平均速度的离散变差函数，根据所述测井平均速度的离散变差函数选择对应的理论变差函数模型进行拟合，建立井点处平均速度的理论变差函数；
第三建立子单元，用于根据所述地震平均速度及所述测井平均速度，计算所述地震平均速度及所述测井平均速度的离散变差函数，根据所述地震平均速度及所述测井平均速度的离散变差函数选择对应的理论变差函数模型进行拟合，建立井点平均速度与地震平均速度的理论互变差函数；
第四建立子单元，用于根据所述沉积相编码，计算所述沉积相编码的离散变差函数，根据所述沉积相编码的离散变差函数选择对应的理论变差函数模型进行拟合，建立沉积相编码的理论变差函数。
可选的，所述协克里金估计单元具体用于：
利用进行协克里金插值，得到井震联合平均速度场；
其中，是在u0位置的速度估计值；是在位置上的测井平均速度值，是位置上的测井平均速度值对应的加权系数值；是在位置上的地震平均速度采样值，是位置上的地震平均速度采样值的加权系数值，是在位置上的沉积相编码，是位置上的沉积相编码的加权系数值，N为沉积相的数量。
基于上述任意技术方案，该系统还可以包括：
转换模块，用于利用所述井震联合平均速度场将时间域构造图转换成对应的深度域构造图。
说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
专业人员还可以进一步意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
以上对本发明所提供的井震联合平均速度场的计算方法及计算系统进行了详细介绍。本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰也落入本发明权利要求的保护范围内。
Claims (8)
1.一种井震联合平均速度场的计算方法，其特征在于，包括：
利用目的层的测井相类型及目的层的地震相类型，划分目的层的沉积相，并对所述沉积相进行编码，得到沉积相编码；
利用地震叠加速度场，通过Dix公式计算得到地震平均速度；
利用所述测井资料，计算得到目的层在各个钻井位置处的测井平均速度；
将所述沉积相编码作为约束项，所述地震平均速度作为次级变量，所述测井平均速度作为主变量，利用进行协克里金插值，得到井震联合平均速度场；
其中，X*(u0)是u0位置的速度估计值；是在位置上的测井平均速度值，是位置上的测井平均速度值对应的加权系数值；是在位置上的地震平均速度采样值，是位置上的地震平均速度采样值的加权系数值，是在位置上的沉积相编码，是位置上的沉积相编码的加权系数值，N为沉积相的数量。
2.如权利要求1所述的计算方法，其特征在于，对所述沉积相进行编码，得到沉积相编码包括：
根据划分后的沉积相的空间范围的大小，按照从大到小的顺序，从整数1开始赋予沉积相对应的数字标识，得到沉积相编码。
3.如权利要求2所述的计算方法，其特征在于，所述将所述沉积相编码作为约束项，所述地震平均速度作为次级变量，所述测井平均速度作为主变量，利用协克里金估计算法进行协克里金插值，得到井震联合平均速度场，包括：
根据所述地震平均速度，所述测井平均速度以及沉积相编码，建立工区所对应的变差函数；
将每个沉积相编码作为约束项，并利用工区所对应的变差函数，根据协克里金估计的无偏性及最优性原则，计算得到各个待估点搜索范围内的各地震平均速度对应的协克里金加权系数值，每个钻井位置处测井平均速度对应的协克里金加权系数值，以及各沉积相编码对应的协克里金加权系数值，进行协克里金估计，得到井震联合平均速度场。
4.如权利要求3所述的计算方法，其特征在于，根据所述地震平均速度，所述测井平均速度以及沉积相编码，建立工区所对应的变差函数，包括：
根据所述地震平均速度，计算所述地震平均速度的离散变差函数，根据所述地震平均速度的离散变差函数选择对应的理论变差函数模型进行拟合，建立地震平均速度的理论变差函数；
根据所述测井平均速度，计算所述测井平均速度的离散变差函数，根据所述测井平均速度的离散变差函数选择对应的理论变差函数模型进行拟合，建立井点处平均速度的理论变差函数；
根据所述地震平均速度及所述测井平均速度，计算所述地震平均速度及所述测井平均速度的离散变差函数，根据所述地震平均速度及所述测井平均速度的离散变差函数选择对应的理论变差函数模型进行拟合，建立井点平均速度与地震平均速度的理论互变差函数；
根据所述沉积相编码，计算所述沉积相编码的离散变差函数，根据所述沉积相编码的离散变差函数选择对应的理论变差函数模型进行拟合，建立沉积相编码的理论变差函数。
5.如权利要求1所述的计算方法，其特征在于，所述利用目的层的测井相类型及目的层的地震相类型，划分目的层的沉积相，包括：
利用测井资料与地震层位特征进行精细井震标定，并根据精细井震标定结果进行层位解释；
根据层位解释后的地震层位，提取有效地震属性集；
通过对所述有效地震属性集进行聚类分析，确定目的层的地震相类型；
以所述目的层的测井相类型作为约束，利用所述目的层的地震相类型，划分所述目的层的沉积相。
6.如权利要求1至5任一项所述的计算方法，其特征在于，还包括：
利用所述井震联合平均速度场将时间域构造图转换成对应的深度域构造图。
7.一种井震联合平均速度场的计算系统，其特征在于，包括：
沉积相模块，用于利用目的层的测井相类型及目的层的地震相类型，划分目的层的沉积相，并对所述沉积相进行编码，得到沉积相编码；
地震平均速度面模块，用于利用地震叠加速度场，通过Dix公式计算得到地震平均速度；
测井平均速度模块，用于利用所述测井资料，计算得到目的层在各个钻井位置处的测井平均速度；
井震联合平均速度场模块，用于将所述沉积相编码作为约束项，所述地震平均速度作为次级变量，所述测井平均速度作为主变量，利用进行协克里金插值，得到井震联合平均速度场；
其中，X*(u0)是在u0位置的速度估计值；是在位置上的测井平均速度值，是位置上的测井平均速度值对应的加权系数值；是在位置上的地震平均速度采样值，是位置上的地震平均速度采样值的加权系数值，是在位置上的沉积相编码，是位置上的沉积相编码的加权系数值，N为沉积相的数量。
8.如权利要求7所述的计算系统，其特征在于，所述井震联合平均速度场模块包括：
变差函数建立单元，用于根据所述地震平均速度，所述测井平均速度以及沉积相编码，建立工区所对应的变差函数；
协克里金估计单元，用于将每个沉积相编码作为约束项，并利用工区所对应的变差函数，根据协克里金估计的无偏性及最优性原则，计算得到各个待估点搜索范围内的各地震平均速度对应的协克里金加权系数值，每个钻井位置处测井平均速度对应的协克里金加权系数值，以及各沉积相编码对应的协克里金加权系数值，进行协克里金估计，得到井震联合平均速度场。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510933562.6A CN105353412B (zh) | 2015-12-14 | 2015-12-14 | 一种井震联合平均速度场的计算方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510933562.6A CN105353412B (zh) | 2015-12-14 | 2015-12-14 | 一种井震联合平均速度场的计算方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105353412A CN105353412A (zh) | 2016-02-24 |
| CN105353412B true CN105353412B (zh) | 2017-08-29 |
Family
ID=55329413
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510933562.6A Active CN105353412B (zh) | 2015-12-14 | 2015-12-14 | 一种井震联合平均速度场的计算方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105353412B (zh) |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106094028A (zh) * | 2016-06-22 | 2016-11-09 | 中国石油化工股份有限公司 | 基于位移场的深度域地震属性提取方法 |
| CN105974475B (zh) * | 2016-06-24 | 2018-02-06 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种时深转换方法 |
| CN106125133B (zh) * | 2016-07-04 | 2018-08-17 | 中国海洋石油总公司 | 一种基于气云区约束下的精细速度建模方法 |
| CN107884827A (zh) * | 2016-09-29 | 2018-04-06 | 中国石油化工股份有限公司 | 地震叠加速度转层速度的方法及装置 |
| CN106772587A (zh) * | 2017-02-23 | 2017-05-31 | 河海大学 | 基于同位多相协同克里金的地震弹性参数相控建模方法 |
| CN108663713B (zh) * | 2017-03-27 | 2021-05-25 | 中国石油化工股份有限公司 | 一种建立深度域构造模型的方法 |
| CN107688200B (zh) * | 2017-04-12 | 2019-09-03 | 恒泰艾普集团股份有限公司 | 针对低幅度构造畸变的速度校正方法 |
| CN107607993B (zh) * | 2017-09-07 | 2019-05-31 | 中国石油大学(北京) | 一种确定叠加速度的方法、装置及系统 |
| CN109655884B (zh) * | 2017-10-10 | 2021-01-12 | 中国石油化工股份有限公司 | 一种速度拾取异常点去除方法及系统 |
| CN110297269B (zh) * | 2018-03-23 | 2020-12-04 | 中国石油化工股份有限公司 | 一种基于改进型bp神经网络的地震数据双向预测插值方法 |
| CN111257945B (zh) * | 2018-11-30 | 2021-03-16 | 中国石油化工股份有限公司 | 一种更新已钻开井段井周地层地震速度的方法 |
| CN112394396A (zh) * | 2019-08-14 | 2021-02-23 | 中国石油天然气股份有限公司 | 一种地震随机反演方法、设备及系统 |
| CN112782756B (zh) * | 2019-11-08 | 2024-06-25 | 中国石油天然气集团有限公司 | 一种基于自适应构造约束的约束层速度反演方法及系统 |
| CN111611714B (zh) * | 2020-05-26 | 2021-08-20 | 中国人民解放军国防科技大学 | 一种超声速混合层降阶方法及系统 |
| CN114167496B (zh) * | 2020-09-11 | 2024-09-03 | 中国石油化工股份有限公司 | 一种复杂构造区建立变速平均速度场的方法 |
| CN116381782A (zh) * | 2023-02-17 | 2023-07-04 | 中海石油(中国)有限公司深圳分公司 | 一种低幅度构造成图方法、设备及存储设备 |
| CN116522688B (zh) * | 2023-06-29 | 2023-09-15 | 北京城建勘测设计研究院有限责任公司 | 一种井控多信息融合工程地质建模方法和装置 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101533103A (zh) * | 2009-04-13 | 2009-09-16 | 中国石油集团川庆钻探工程有限公司 | 地震波速度场构建中井约束横向可变h-v曲线构建方法 |
| CN102590864A (zh) * | 2011-12-31 | 2012-07-18 | 中国石油集团西北地质研究所 | 两步法层析反演近地表建模方法 |
| CN103543478A (zh) * | 2012-07-11 | 2014-01-29 | 中国石油化工股份有限公司 | 地质形态插值的km方法 |
| CN104155687A (zh) * | 2013-05-15 | 2014-11-19 | 中国石油化工股份有限公司 | 一种相控叠后声波阻抗反演方法 |
| CN104200115A (zh) * | 2014-09-12 | 2014-12-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于地质统计学模拟的全地层速度建模方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7525874B2 (en) * | 2006-04-21 | 2009-04-28 | Prism Seismic Inc. | Method for converting seismic data from the time domain to the depth domain |
-
2015
- 2015-12-14 CN CN201510933562.6A patent/CN105353412B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101533103A (zh) * | 2009-04-13 | 2009-09-16 | 中国石油集团川庆钻探工程有限公司 | 地震波速度场构建中井约束横向可变h-v曲线构建方法 |
| CN102590864A (zh) * | 2011-12-31 | 2012-07-18 | 中国石油集团西北地质研究所 | 两步法层析反演近地表建模方法 |
| CN103543478A (zh) * | 2012-07-11 | 2014-01-29 | 中国石油化工股份有限公司 | 地质形态插值的km方法 |
| CN104155687A (zh) * | 2013-05-15 | 2014-11-19 | 中国石油化工股份有限公司 | 一种相控叠后声波阻抗反演方法 |
| CN104200115A (zh) * | 2014-09-12 | 2014-12-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于地质统计学模拟的全地层速度建模方法 |
Non-Patent Citations (5)
| Title |
|---|
| "A Geostatistical Approach for Two Dimensional Seismic Velocity Modelling";Claprood.M 等;《Geological Survey of Canada》;20111231;第1-21页 * |
| "储层建模中地质统计学整合地震数据的方法及研究进展";印兴耀 等;《石油地球物理勘探》;20020831;第37卷(第4期);第423-430页 * |
| "地震数据约束相控建模方法研究";周单;《中国优秀硕士学位论文全文数据库（基础科学辑）》;20110415(第04（2011）期);第1-64页 * |
| "应用序贯指示模拟方法模拟沉积微相";冯国庆 等;《西南石油学院学报》;20010430;第23卷(第2期);第1-4页 * |
| 谢清惠 等."贝叶斯-克里金估计方法在变速成图中的应用".《石油物探》.2012,第51卷(第3期),第251-256页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105353412A (zh) | 2016-02-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105353412B (zh) | 一种井震联合平均速度场的计算方法及系统 | |
| CN113759425B (zh) | 井震联合评价深层古岩溶储层充填特征的方法与系统 | |
| CN113759424B (zh) | 基于频谱分解和机器学习的岩溶储层充填分析方法和系统 | |
| CN108802812B (zh) | 一种井震融合的地层岩性反演方法 | |
| CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
| US8670288B2 (en) | Velocity model for well time-depth conversion | |
| CN109541685B (zh) | 一种河道砂体识别方法 | |
| CN107065013B (zh) | 一种地震尺度下的层速度确定方法及装置 | |
| CN108005646B (zh) | 基于随钻电磁波测井资料的地层各向异性电阻率提取方法 | |
| EP2160633A1 (en) | Creating an absorption parameter model | |
| CN105588883B (zh) | 三维岩石力学参数获取方法和系统 | |
| CN102937720A (zh) | 井控提高地震资料分辨率的方法 | |
| US8599644B2 (en) | Velocity models for a single well and for a set of wells | |
| CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
| CN114994758B (zh) | 碳酸盐岩断控储层的波阻抗提取与结构表征方法和系统 | |
| CN113267822B (zh) | 基于地形约束因子权重优化提高海底地形反演精度的方法 | |
| CN111856566A (zh) | 湖相滩坝砂体中薄储层预测方法及装置 | |
| CN109471163A (zh) | 基于地质体连续性建模的高精度井拓反演方法 | |
| CN111983683A (zh) | 一种少井条件下盐下湖相灰岩储层预测方法和系统 | |
| CN109100791B (zh) | 基于纵横向空间约束的速度反演方法 | |
| CN108508481B (zh) | 一种纵波转换波地震数据时间匹配的方法、装置及系统 | |
| CN106869919B (zh) | 一种三角洲前缘薄油藏识别方法 | |
| CN110703329B (zh) | 基于弱振幅地震反射形成机制的岩性油藏边界确定方法 | |
| CN112147700A (zh) | 速度异常区的低频模型构建方法及系统 | |
| CN113534263B (zh) | 一种不依赖测井资料的含油气饱和度预测方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |