CN109932749B - 一种井震标定方法、装置 - Google Patents
一种井震标定方法、装置 Download PDFInfo
- Publication number
- CN109932749B CN109932749B CN201910164555.2A CN201910164555A CN109932749B CN 109932749 B CN109932749 B CN 109932749B CN 201910164555 A CN201910164555 A CN 201910164555A CN 109932749 B CN109932749 B CN 109932749B
- Authority
- CN
- China
- Prior art keywords
- sampling point
- well
- calculated
- depth
- wells
- 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
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000005070 sampling Methods 0.000 claims abstract description 380
- 239000003550 marker Substances 0.000 claims abstract description 59
- 238000011160 research Methods 0.000 claims abstract description 34
- 238000007689 inspection Methods 0.000 claims description 37
- 230000002159 abnormal effect Effects 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000011835 investigation Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect 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
- 239000011435 rock Substances 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种井震标定方法、装置，属于油藏地球物理领域。所述方法包括：确定研究区域内的地震标志层，以及地震标志层对应的地质层位；确定研究区域内的多个待计算井；对于每个待计算井，根据地震解释成果确定第一采样点，根据第一采样点，确定待计算井内每个采样点的采样点时深关系；根据每个采样点的采样点时深关系，汇总得到待计算井的计算井时深关系；根据每个待计算井的计算井时深关系，对研究区域内多个待计算井进行井震标定。本发明通过上述方法可以快速完成井震标定工作，不需要进行人工分析，且可以同时对多个井进行井震标定，在保证准确性的前提下，节约了时间，提高了效率，缩短了工作周期。
Description
技术领域
本发明涉及油藏地球物理领域，特别涉及一种井震标定方法、装置。
背景技术
井震标定是联系地震与地质的桥梁，精细的地震标定是进行层位解释和油藏描述的基础和前提。井震标定不准确可能会导致储层构造深度误差、储层预测不准，以及油藏认识偏移。
现有井震标定方法是通过声波、密度测井资料计算反射系数，利用发射系数与子波的褶积构建与地震道相似的合成地震记录，地震解释人员通过对合成地震记录与井旁地震道作对比与调整来进行井震标定。
上述井震标定方法需要人工分析，且工作量大，对于技术人员的地震解释经验具有依赖性，导致分析结果具有较强的主观性，效率低，工作周期长。
发明内容
本发明实施例提供了一种井震标定方法、装置,能够解决现有技术中分析结果主观性强、效率低以及工作周期长的问题。所述技术方案如下：
一方面，提供了一种井震标定方法，所述方法包括：
确定研究区域内的地震标志层，以及所述地震标志层对应的地质层位；
确定所述研究区域内的多个待计算井，所述多个待计算井为包含声波曲线的多个井；
对于每个待计算井，根据地震解释成果确定第一采样点，所述第一采样点为所述地震标志层与所述待计算井的交叉点；
根据所述第一采样点，确定所述待计算井内每个采样点的采样点时深关系；
根据所述待计算井内每个采样点的采样点时深关系，汇总得到所述待计算井的计算井时深关系；
根据所述每个待计算井的计算井时深关系，对所述研究区域内多个待计算井进行井震标定。
在一种可能的实现方式中，所述确定所述研究区域内的多个待计算井，所述多个待计算井为包含声波曲线的多个井，包括：
在所述研究区域内选取具有声波曲线的多个井；
对所述具有声波曲线的多个井进行声波曲线质量检查，得到所述多个井的检查结果；
根据所述多个井的检查结果，将所述多个井中检查结果为合格的井作为所述待计算井。
在另一种可能的实现方式中，所述对所述具有声波曲线的多个井进行声波曲线质量检查，得到所述多个井的检查结果，包括：
对于每个井，确定所述井的声波曲线的数值与单位是否相匹配，以及确定所述声波曲线中是否具有异常值，以及确定所述声波曲线是否覆盖所述井的全井段；
当所述井的声波曲线的数值与单位相匹配、所述声波曲线中无异常值且所述声波曲线覆盖所述井的全井段时，确定所述井的检查结果为合格。
在另一种可能的实现方式中，所述根据所述第一采样点，确定所述待计算井内每个采样点的采样点时深关系，包括：
根据所述第一采样点，确定所述第一采样点的第一深度以及所述第一采样点的第一双程旅行时间；
根据所述第一采样点，确定所述第一采样点的第一深度以及所述第一采样点的第一双程旅行时间；
根据所述待计算井的声波曲线确定所述待计算井内每个采样点的层速度；
根据所述第一双程旅行时间以及所述第一深度，确定所述待计算井内其他采样点的双程旅行时间；
根据所述待计算井内每个采样点的双程旅行时间、每个采样点的深度以及每个采样点的层速度，确定每个采样点的采样点时深关系。
在另一种可能的实现方式中，所述根据所述待计算井的声波曲线确定所述待计算井内每个采样点的层速度，包括：
根据所述声波曲线确定每个采样点对应的声波曲线值；
将所述第一采样点对应的第一声波曲线值与第二采样点对应的第二声波曲线值的平均值作为第一采样点的层速度，所述第二采样点为在计算方向上距离所述第一采样点最近的采样点；
根据所述第二采样点对应的第二声波曲线值与第三采样点对应的第三声波曲线值确定所述第二采样点的层速度，直到确定所述待计算井内每个采样点的层速度，所述第三采样点为在计算方向上距离所述第二采样点最近的采样点。
在另一种可能的实现方式中，所述根据所述第一双程旅行时间、所述第一深度以及所述每个采样点的层速度，确定所述待计算井内其他采样点的双程旅行时间，包括：
确定第二采样点和所述第一采样点之间的第一测井采样深度差，所述第一测井采样深度差为所述第一深度和所述第二采样点的第二深度之间的差值，所述第二采样点为在计算方向上距离所述第一采样点最近的采样点；
根据所述第一双程旅行时间、所述第一测井采样深度差和所述第二采样点的第二层速度，通过以下公式一，确定所述第二采样点的第二双程旅行时间；公式一：T(2)＝T(1)+2×△d(1)/VINT (2)；
其中，T(2)表示第二采样点的第二双程旅行时间，T(1)表示第一采样点的双程旅行时间，△d(1)表示第一测井采样深度差，VINT(2)表示第二层速度；
根据所述第二双程旅行时间以及所述第二深度，确定所述待计算井内第三采样点的第三双程旅行时间，直到确定出所述待计算井内每个采样点的双程旅行时间为止，所述第三采样点为在计算方向上距离所述第二采样点最近的采样点。
另一方面，提供了一种井震标定装置，所述装置包括：
第一确定模块，被配置为确定研究区域内的地震标志层，以及所述地震标志层对应的地质层位；
第二确定模块，被配置为确定所述研究区域内的多个待计算井，所述多个待计算井为包含声波曲线的多个井；
第三确定模块，被配置为对于每个待计算井，根据地震解释成果确定第一采样点，所述第一采样点为所述地震标志层与所述待计算井的交叉点；
第四确定模块，被配置为根据所述第一采样点，确定所述待计算井内每个采样点的采样点时深关系；
汇总模块，被配置为根据所述待计算井内每个采样点的采样点时深关系，汇总得到所述待计算井的计算井时深关系；
标定模块，被配置为根据所述每个待计算井的计算井时深关系，对所述研究区域内多个待计算井进行井震标定。
在一种可能的实现方式中，所述第二确定模块，还被配置为在所述研究区域内选取具有声波曲线的多个井；对所述具有声波曲线的多个井进行声波曲线质量检查，得到所述多个井的检查结果；根据所述多个井的检查结果，将所述多个井中检查结果为合格的井作为所述待计算井。
在另一种可能的实现方式中，所述第二确定模块，还被配置为对于每个井，确定所述井的声波曲线的数值与单位是否相匹配，以及确定所述声波曲线中是否具有异常值，以及确定所述声波曲线是否覆盖所述井的全井段；当所述井的声波曲线的数值与单位相匹配、所述声波曲线中无异常值且所述声波曲线覆盖所述井的全井段时，确定所述井的检查结果为合格。
在另一种可能的实现方式中，所述第四确定模块，还被配置为根据所述第一采样点，确定所述第一采样点的第一深度以及所述第一采样点的第一双程旅行时间；根据所述待计算井的声波曲线确定所述待计算井内每个采样点的层速度；根据所述第一双程旅行时间以及所述第一深度，确定所述待计算井内其他采样点的双程旅行时间；根据所述待计算井内每个采样点的双程旅行时间、每个采样点的深度以及每个采样点的层速度，确定每个采样点的采样点时深关系。
在另一种可能的实现方式中，所述第四确定模块，还被配置为根据所述声波曲线确定每个采样点对应的声波曲线值；将所述第一采样点对应的第一声波曲线值与第二采样点对应的第二声波曲线值的平均值作为第一采样点的层速度，所述第二采样点为在计算方向上距离所述第一采样点最近的采样点；根据所述第二采样点对应的第二声波曲线值与第三采样点对应的第三声波曲线值确定第二采样点的层速度，直到确定所述待计算井内每个采样点的层速度，所述第三采样点为在计算方向上距离所述第二采样点最近的采样点。
在另一种可能的实现方式中，所述第四确定模块，还被配置为确定第二采样点和所述第一采样点之间的第一测井采样深度差，所述第一测井采样深度差为所述第一深度和所述第二采样点的第二深度之间的差值，所述第二采样点为在计算方向上距离所述第一采样点最近的采样点；根据所述第一双程旅行时间、所述第一测井采样深度差和所述第二采样点的第二层速度，通过以下公式一，确定所述第二采样点的第二双程旅行时间；
公式一：T(2)＝T(1)+2×△d(1)/VINT(2)；其中，T(2)表示第二采样点的第二双程旅行时间，T(1)表示第一采样点的双程旅行时间，△d(1)表示第一测井采样深度差，VINT(2)表示第二层速度；根据所述第二双程旅行时间以及所述第二深度，确定所述待计算井内第三采样点的第三双程旅行时间，直到确定出所述待计算井内每个采样点的双程旅行时间为止，所述第三采样点为在计算方向上距离所述第二采样点最近的采样点。
本发明实施例提供的技术方案带来的有益效果是：
通过确定研究区域内的地震标志层，以及地震标志层对应的地质层位，对于每个待计算井，根据地震解释成果确定第一采样点，根据第一采样点，确定待计算井内每个采样点的采样点时深关系，汇总得到待计算井的计算井时深关系，对研究区域内多个待计算井进行井震标定，不需要进行人工分析，且可以同时对多个井进行井震标定，在保证准确性的前提下，节约了时间，提高了效率，缩短了工作周期。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种井震标定方法的流程图；
图2是本发明实施例提供的另一种井震标定方法的流程图；
图3是本发明实施例提供的另一种井震标定方法的流程图；
图4是本发明实施例提供的一种地震标志层与地质层位对应关系图；
图5是本发明实施例提供的一种井震标定装置的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
图1是本发明实施例提供的一种井震标定方法的流程图。参见图1，本发明提供了一种井震标定方法，该方法为基于标志层控制的快速井震标定方法，适用于井网密度大于10口/km2的研究区域内的所有地层。从研究区域内选择一个地震标志层，并确定与该地震标志层相对应的地质层位，选取测井曲线中包含声波曲线的多个待计算井，对于每个待计算井，根据地震解释成果确定第一采样点，根据第一采样点得到其他采样点的时深关系，得到每个待计算井的时深关系，从而对研究区域内的多个待计算井进行井震标定。
在本发明实施例中，通过确定研究区域内的地震标志层以及地震标志层对应的地质层位，对于每个待计算井，根据地震解释成果确定第一采样点，根据第一采样点，确定待计算井内每个采样点的采样点时深关系，汇总得到待计算井的计算井时深关系，对研究区域内多个待计算井进行井震标定，在保证准确性的前提下，节约了时间，提高了效率，缩短了工作周期。
图2是本发明实施例提供的另一种井震标定方法的流程图。参见图2，该实施例包括：
在步骤201中，终端确定研究区域内的地震标志层，以及该地震标志层对应的地质层位。
在步骤202中，终端确定研究区域内的多个待计算井，该多个待计算井为包含声波曲线的多个井。
在步骤203中，对于每个待计算井，终端根据地震解释确定第一采样点，第一采样点为地震标志层与待计算井的交叉点。
在步骤204中，终端根据第一采样点，确定待计算井内每个采样点的采样点时深关系。
在步骤205中，终端根据待计算井内每个采样点的采样点时深关系，汇总得到待计算井的计算井时深关系。
在步骤206中，终端根据每个待计算井的计算井时深关系，对研究区域内多个待计算井进行井震标定。
在本发明实施例中，通过确定研究区域内的地震标志层以及地震标志层对应的地质层位，对于每个待计算井，根据地震解释成果确定第一采样点，根据第一采样点，确定待计算井内每个采样点的采样点时深关系，汇总得到待计算井的计算井时深关系，对研究区域内多个待计算井进行井震标定，在保证准确性的前提下，节约了时间，提高了效率，缩短了工作周期。
图3是本发明实施例提供的另一种井震标定方法的流程图。参见图3，该实施例包括：
在步骤301中，终端确定研究区域内的地震标志层，以及该地震标志层对应的地质层位。
终端首先确定研究区域，该研究区域为井网密度大于10口/km2的任一研究区域，在本发明实施例中，以大庆长恒油田密井网区A区块作为研究区域为例进行说明。终端确定研究区域后，在该研究区域内选取一个地震标志层，其中，该研究区域中包含多个地震标志层，可以随机选取多个地震标志层中的一个地震标志层，也可以选择该研究区域目的层附近的标志层，该目的层可以为该研究区域内需要进行研究的地质层。确定地震标志层后，根据地震解释成果中包含的地震标志层与地质层位的关系，获取与该地震标志层相对应的地质层位。
例如，参见图4所示，选取研究区域内地震标志层特征最明显的层位T1，根据地震解释成果获取该T1层对应的地质层位为SⅡ1层。
需要说明的是，上述地震标志层，是指一层或一组具有明显特征可作为地层对比标志的岩层，具有化石和岩性特征明显、层位稳定、分布范围广、易于鉴别的特点，通常情况下该地震标志层为地震解释成果中已存在的资料。
在步骤302中，终端确定研究区域内的多个待计算井。
其中，多个待计算井为包含声波曲线的多个井，相应的，本步骤可以通过以下步骤(1)至(3)来实现：
(1)终端在研究区域内选取具有声波曲线的多个井。
终端在研究区域内的井中，选取测井曲线中包含声波曲线的多个井，之所以选取带有声波曲线的井，是因为在进行时深转换时需要用到速度这一信息，声波曲线中包含进行时深转换时需要用到的速度信息，而选取其他曲线不能得到进行时深转换时需要用到的速度信息。
(2)终端对具有声波曲线的多个井进行声波曲线质量检查，得到多个井的检查结果。
对于每个井，终端确定该井的声波曲线的数值与单位是否相匹配，以及确定该声波曲线中是否具有异常值，以及确定该声波曲线是否覆盖井的全井段。当终端确定该井的声波曲线的数值与单位相匹配、声波曲线中无异常值且声波曲线覆盖井的全井段时，则确定该井的检查结果为合格；当终端确定该井的声波曲线的数值与单位不相匹配或声波曲线中包含异常值或声波曲线未覆盖井的全井段时，则确定该井的检查结果为不合格。终端对多个井进行声波质量检查后的得到多个井的检查结果。
(3)终端根据多个井的检查结果，将多个井中检查结果为合格的井作为待计算井。
在步骤303中，对于每个待计算井，终端根据地震解释成果确定第一采样点。
终端根据地震解释成果确定第一采样点，该第一采样点为地震标志层与待计算井的交叉点。
对于每个待计算井，该待计算井具有垂直方向的井轨迹，地震标志层可以认为是接近水平方向的二维面，则待计算井与地震标志层相交，在待计算井与地震标志层相交的位置选择一点作为第一采样点，该第一采样点可以为在相交位置随机选择的点，也可以为相交位置的中心点，还可以选择待计算井与地震标志层以及地质层位相交的一点作为第一采样点。本发明对此不进行具体限制。
例如参见图4所示，地质层位SII1层与地震标志层T1相对应，井1为待计算井，井1与SII1层以及T1层交于一点，将该点作为第一采样点。
在步骤304中，对于每个待计算井，终端根据第一采样点，确定第一采样点的第一深度以及第一采样点的第一双程旅行时间。
第一采样点的第一深度为地震标志层对应的地质层位的深度，终端中存储有地震标志层与地质层位的对应关系，根据该对应关系，获取地震标志层包含的第一采样点在对应的地质层位中的第一地质层位点，获取该第一地质层位点的深度，该深度可以用DEP(deep，深度)来表示，将该第一地质层位点的深度作为第一采样点的第一深度，该第一深度可以用DEP(1)来表示。
终端获取地震解释结果，该地震解释结果中包含待计算井中第一采样点的双程旅行时间，该双程旅行时间可以用T(t ime，时间)来表示，终端根据该地震解释结果获取第一采样点的第一双程旅行时间，该第一双程旅行时间可以用T(1)来表示。
在步骤305中，对于每个待计算井，终端根据第一采样点的第一深度、第一双程旅行时间以及每个采样点的层速度，确定待计算井内每个采样点的采样点时深关系。
本步骤可以通过以下步骤(1)至(3)来实现：
(1)终端根据待计算井的声波曲线确定该待计算井内每个采样点的层速度。
终端获取待计算井的声波曲线，该声波曲线中包含该待计算井内每个采样点的声波曲线值，终端根据该声波曲线获取每个采样点对应的声波曲线值。将第一采样点对应的第一声波曲线值与第二采样点对应的第二声波曲线值的平均值作为第一采样点的层速度，其中，该第二采样点为在计算方向上距离第一采样点最近的采样点。根据第二采样点对应的第二声波曲线值与第三采样点对应的第三声波曲线值确定第二采样点的层速度，即将第二采样点对应的第二声波曲线值以及第三采样点对应的声波曲线值的平均值，作为第二采样点的第二层速度，第三采样点为在计算方向上距离第二采样点最近的采样点。按照同样的方式获取第三采样点的第三层速度，直到确定待计算井内每个采样点的层速度，层速度可以用VINT(velocity，速度)来表示，则每个采样点的层速度可以用VINT(n)表示，n表示第n个采样点。
(2)终端根据第一双程旅行时间、第一深度以及每个采样点的层速度，确定待计算井内其他采样点的双程旅行时间。
终端根据第一深度，得到该第一深度和第二采样点的第二深度之间的差值，将该差值作为第二采样点和第一采样点之间的第一测井采样深度差，该第二采样点为在计算方向上与第一采样点距离最近的采样点。
在步骤(1)中，终端已经获取了每个采样点的层速度VINT(n)，n表示第n个采样点，终端根据第一双程旅行时间、第一测井采样深度差和第二采样点的第二层速度，通过以下公式一，确定第二采样点的第二双程旅行时间；公式一：T(2)＝T(1)+2×△d(1)/VINT(2)。其中，T(2)表示第二采样点的第二双程旅行时间，T(1)表示第一采样点的双程旅行时间，△d(1)表示第一测井采样深度差，VINT(2)表示第二层速度。终端根据第二双程旅行时间以及第二深度，确定待计算井内第三采样点的第三双程旅行时间，直到确定出待计算井内计算方向上每个采样点的双程旅行时间为止，该第三采样点为在计算方向上距离所述第二采样点最近的采样点。
其中，终端确定待计算井内计算方向上每个采样点的双程旅行时间时，可以通过公式二来进行计算，该公式二用于根据已知的第n个采样点的双程旅行时间、与已知的第n个采样点相邻的第n+1个采样点的深测井采样深度差以及与已知第n个采样点相邻的第n+1个采样点的层速度，获取与已知第n个采样点相邻的第n+1个采样点的双程旅行时间。
该公式二如下：
T(n+1)＝T(n)+2×△d(n)/VINT(n+1)
其中，T(n+1)表示第n+1个采样点的双程旅行时间，T(n)表示第n个采样点的双程旅行时间，△d(n)表示第n个采样点和第n+1个采样点的测井采样深度差，VINT(n+1)表示根据所述声波曲线计算的第n+1个采样点的层速度。
终端获取地震解释成果，该地震解释成果中包含每个采样点对应的深度，终端获取每个采样点对应的深度，依次计算两个采样点之间的测井采样深度差△d。第二采样点和第三采样点的第二测井采样深度差可以表示为△d(2)，第三采样点和第四采样点的第三测井采样深度差可以表示为△d(3)，以此类推。
例如，终端通过公式二获取第三采样点的第三双程旅行时间T(3)时，将第二双程旅行时间T(2)、第二测井采样深度差△d(2)以及第三采样点的第三层速度VINT(3)代入上述公式二中，得到T(3)＝T(2)+2×△d(2)/VINT(3)；终端将获取到的第三采样点的第三双程旅行时间T(3)、第三测井采样深度差△d(3)以及第四采样点的第四层速度VINT(4)代入上述公式二中，得到第四采样点的第四双程旅行时间T(4)，T(4)＝T(3)+2×△d(3)/VINT(4)。依次进行计算，直到确定出待计算井内计算方向上每个采样点的双程旅行时间为止。
需要说明的是，由于标志层与待计算井相交的位置不在待计算井的两端，即第一采样点不在待计算井的两端，终端需要先确定计算方向，该计算方向分为由标志层出发指向地面(即由第一采样点出发指向井口)的第一计算方向，以及由标志层出发指向地底(即由第一采样点出发指向井底)的第二计算方向。其中第一计算方向和第二计算方向仅是为了方便区分进行的划分，不存在先后顺序。终端在第一计算方向上进行的计算与在第二计算方向上进行的计算遵循同样的计算方式，在本发明中，不再重复介绍两个计算方向的计算，以第二计算方向上采样点的计算为例进行说明。
还需要说明的是，第一计算方向上每个采样点的双程旅行时间和第二计算方向上每个采样点的双程旅行时间是可以同时进行计算的，终端最终获取到待计算井内每个采样点的双程旅行时间。
(3)终端根据待计算井内每个采样点的双程旅行时间、每个采样点的深度以及每个采样点的层速度，确定每个采样点的采样点时深关系。
对于每个采样点，该采样点的双程旅行时间与该采样点的深度成正比，该采样点的双程旅行时间与该采样点的层速度成反比，双程旅行时间与深度和层速度之间的关系即为采样点时深关系。T(n)＝DEP(n)/VINT(n)，其中，n表示第n个采样点。
例如，对于第四采样点，该第四采样点的第四双程旅行时间可以表示为T(4)，该第四采样点的深度可以表示为DEP(4)，该第四采样点的层速度的为VINT(4)，则该第四采样点的采样点时深关系为T(4)＝DEP(4)/VINT(4)。
在步骤306中，终端根据待计算井内每个采样点的采样点时深关系，汇总得到待计算井的计算井时深关系。
对于每一个待计算井，终端将上述步骤获取到的待计算井内的每个采样点的采样点时深关系进行汇总，汇总得到该待计算井的计算井时深关系。
在步骤307中，终端根据每个待计算井的计算井时深关系，对研究区域内多个待计算井进行井震标定。
终端可以同时对研究区域内的多个待计算井同时进行计算井时深关系的建立，多个待计算井建立完成计算井时深关系后，终端根据每个待计算井的计算井时深关系对研究区域内多个待计算井进行井震标定。
在本发明实施例中，通过确定研究区域内的地震标志层以及地震标志层对应的地质层位，对于每个待计算井，根据地震解释成果确定第一采样点，根据第一采样点，结合测井曲线中的声波曲线，确定待计算井内每个采样点的采样点时深关系，汇总得到待计算井的计算井时深关系，对研究区域内多个待计算井进行井震标定，在保证准确性的前提下，节约了时间，提高了效率，缩短了工作周期。
图5是本发明实施例提供的一种井震标定装置结构示意图，参见图5所示，该装置包括第一确定模块501、第二确定模块502、第三确定模块503、第四确定模块504、汇总模块505以及标定模块506。
第一确定模块501，被配置为确定研究区域内的地震标志层，以及地震标志层对应的地质层位；
第二确定模块502，被配置为确定研究区域内的多个待计算井，多个待计算井为包含声波曲线的多个井；
第三确定模块503，被配置为对于每个待计算井，根据地震解释成果确定第一采样点，第一采样点为地震标志层与待计算井的交叉点；
第四确定模块504，被配置为根据第一采样点，确定待计算井内每个采样点的采样点时深关系；
汇总模块505，被配置为根据每个采样点的采样点时深关系，汇总得到待计算井的计算井时深关系；
标定模块506，被配置为根据每个待计算井的计算井时深关系，对研究区域内多个待计算井进行井震标定。
在一种可能的实现方式中，第二确定模块502，还被配置为在研究区域内选取具有声波曲线的多个井；对具有声波曲线的多个井进行声波曲线质量检查，得到多个井的检查结果；根据多个井的检查结果，将多个井中检查结果为合格的井作为待计算井。
在另一种可能的实现方式中，第二确定模块502，还被配置为对于每个井，确定井的声波曲线的数值与单位是否相匹配，以及确定声波曲线中是否具有异常值，以及确定声波曲线是否覆盖井的全井段；当井的声波曲线的数值与单位相匹配、声波曲线中无异常值且声波曲线覆盖井的全井段时，确定井的检查结果为合格。
在另一种可能的实现方式中，第四确定模块504，还被配置为根据第一采样点，确定第一采样点的第一深度以及第一采样点的第一双程旅行时间；根据待计算井的声波曲线确定待计算井内每个采样点的层速度；根据第一双程旅行时间、第一深度以及每个采样点的层速度，确定待计算井内其他采样点的双程旅行时间；根据每个采样点的双程旅行时间、每个采样点的深度以及每个采样点的层速度，确定每个采样点的采样点时深关系。
在另一种可能的实现方式中，第四确定模块504，还被配置为根据声波曲线确定每个采样点对应的声波曲线值；将第一采样点对应的第一声波曲线值与第二采样点对应的第二声波曲线值的平均值作为第一采样点的层速度，第二采样点为在计算方向上距离第一采样点最近的采样点；根据第二采样点对应的第二声波曲线值与第三采样点对应的第三声波曲线值确定第二采样点的层速度，直到确定待计算井内每个采样点的层速度，第三采样点为在计算方向上距离第二采样点最近的采样点。
在另一种可能的实现方式中，第四确定模块504，还被配置为确定第二采样点和所述第一采样点之间的第一测井采样深度差，所述第一测井采样深度差为第一深度和第二采样点的第二深度之间的差值，第二采样点为在计算方向上与第一采样点距离最近的采样点；根据第一双程旅行时间、第一测井采样深度差和第二采样点的第二层速度，通过以下公式一，确定第二采样点的第二双程旅行时间；公式一：T(2)＝T(1)+2×△d(1)/VINT(2)；其中，T(2)表示第二采样点的第二双程旅行时间，T(1)表示第一采样点的双程旅行时间，△d(1)表示第一测井采样深度差，VINT(2)表示第二层速度；根据第二双程旅行时间以及第二深度，确定待计算井内第三采样点的第三双程旅行时间，直到确定出待计算井内每个采样点的双程旅行时间为止，第三采样点为在计算方向上距离第二采样点最近的采样点。
在本发明实施例中，通过确定研究区域内的地震标志层以及地震标志层对应的地质层位，对于每个待计算井，根据地震解释成果确定第一采样点，根据第一采样点，结合声波曲线，确定待计算井内每个采样点的采样点时深关系，汇总得到待计算井的计算井时深关系，对研究区域内多个待计算井进行井震标定，在保证准确性的前提下，节约了时间，提高了效率，缩短了工作周期。
上述所有可选技术方案，可以采用任意结合形成本发明的可选实施例，在此不再一一赘述。
本领域技术人员在考虑说明书及实践这里公开的发明后，将容易想到本发明的其它实施方案。本发明旨在涵盖本发明的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的，本发明的真正范围和精神由下面的权利要求指出。
应当理解的是，本发明并不局限于上面已经描述并在附图中示出的精确结构，并且可以在不脱离其范围进行各种修改和改变。本发明的范围仅由所附的权利要求来限制。
Claims (10)
1.一种井震标定方法，其特征在于，所述方法包括：
确定研究区域内的地震标志层，以及所述地震标志层对应的地质层位；
确定所述研究区域内的多个待计算井，所述多个待计算井为包含声波曲线的多个井；
对于每个待计算井，根据地震解释成果确定第一采样点，所述第一采样点为所述地震标志层与所述待计算井的交叉点；
根据所述第一采样点，确定所述第一采样点的第一深度以及所述第一采样点的第一双程旅行时间；
根据所述待计算井的声波曲线确定所述待计算井内每个采样点的层速度；
根据所述第一双程旅行时间、所述第一深度以及所述每个采样点的层速度，确定所述待计算井内其他采样点的双程旅行时间；
根据所述待计算井内每个采样点的双程旅行时间、每个采样点的深度以及每个采样点的层速度，确定每个采样点的采样点时深关系；
根据所述待计算井内每个采样点的采样点时深关系，汇总得到所述待计算井的计算井时深关系；
根据所述每个待计算井的计算井时深关系，对所述研究区域内多个待计算井进行井震标定。
2.根据权利要求1所述的方法，其特征在于，所述确定所述研究区域内的多个待计算井，所述多个待计算井为包含声波曲线的多个井，包括：
在所述研究区域内选取具有声波曲线的多个井；
对所述具有声波曲线的多个井进行声波曲线质量检查，得到所述多个井的检查结果；
根据所述多个井的检查结果，将所述多个井中检查结果为合格的井作为所述待计算井。
3.根据权利要求2所述的方法，其特征在于，所述对所述具有声波曲线的多个井进行声波曲线质量检查，得到所述多个井的检查结果，包括：
对于每个井，确定所述井的声波曲线的数值与单位是否相匹配，以及确定所述声波曲线中是否具有异常值，以及确定所述声波曲线是否覆盖所述井的全井段；
当所述井的声波曲线的数值与单位相匹配、所述声波曲线中无异常值且所述声波曲线覆盖所述井的全井段时，确定所述井的检查结果为合格。
4.根据权利要求1所述的方法，其特征在于，所述根据所述待计算井的声波曲线确定所述待计算井内每个采样点的层速度，包括：
根据所述声波曲线确定每个采样点对应的声波曲线值；
将所述第一采样点对应的第一声波曲线值与第二采样点对应的第二声波曲线值的平均值作为第一采样点的层速度，所述第二采样点为在计算方向上距离所述第一采样点最近的采样点；
根据所述第二采样点对应的第二声波曲线值与第三采样点对应的第三声波曲线值确定所述第二采样点的层速度，直到确定所述待计算井内每个采样点的层速度，所述第三采样点为在计算方向上距离所述第二采样点最近的采样点。
5.根据权利要求1所述的方法，其特征在于，所述根据所述第一双程旅行时间、所述第一深度以及所述每个采样点的层速度，确定所述待计算井内其他采样点的双程旅行时间，包括：
确定第二采样点和所述第一采样点之间的第一测井采样深度差，所述第一测井采样深度差为所述第一深度和所述第二采样点的第二深度之间的差值，所述第二采样点为在计算方向上距离所述第一采样点最近的采样点；
根据所述第一双程旅行时间、所述第一测井采样深度差和所述第二采样点的第二层速度，通过以下公式一，确定所述第二采样点的第二双程旅行时间；
公式一：T(2)＝T(1)+2×△d(1)/VINT (2)；
其中，T(2)表示第二采样点的第二双程旅行时间，T(1)表示第一采样点的双程旅行时间，△d(1)表示第一测井采样深度差，VINT(2)表示第二层速度；
根据所述第二双程旅行时间以及所述第二深度，确定所述待计算井内第三采样点的第三双程旅行时间，直到确定出所述待计算井内每个采样点的双程旅行时间为止，所述第三采样点为在计算方向上距离所述第二采样点最近的采样点。
6.一种井震标定装置，其特征在于，所述装置包括：
第一确定模块，被配置为确定研究区域内的地震标志层，以及所述地震标志层对应的地质层位；
第二确定模块，被配置为确定所述研究区域内的多个待计算井，所述多个待计算井为包含声波曲线的多个井；
第三确定模块，被配置为对于每个待计算井，根据地震解释成果确定第一采样点，所述第一采样点为所述地震标志层与所述待计算井的交叉点；
第四确定模块，被配置为根据所述第一采样点，确定所述第一采样点的第一深度以及所述第一采样点的第一双程旅行时间；根据所述待计算井的声波曲线确定所述待计算井内每个采样点的层速度；根据所述第一双程旅行时间、所述第一深度以及所述每个采样点的层速度，确定所述待计算井内其他采样点的双程旅行时间；根据所述待计算井内每个采样点的双程旅行时间、每个采样点的深度以及每个采样点的层速度，确定每个采样点的采样点时深关系；
汇总模块，被配置为根据所述待计算井内每个采样点的采样点时深关系，汇总得到所述待计算井的计算井时深关系；
标定模块，被配置为根据所述每个待计算井的计算井时深关系，对所述研究区域内多个待计算井进行井震标定。
7.根据权利要求6所述的装置，其特征在于，所述第二确定模块，还被配置为在所述研究区域内选取具有声波曲线的多个井；对所述具有声波曲线的多个井进行声波曲线质量检查，得到所述多个井的检查结果；根据所述多个井的检查结果，将所述多个井中检查结果为合格的井作为所述待计算井。
8.根据权利要求7所述的装置，其特征在于，所述第二确定模块，还被配置为对于每个井，确定所述井的声波曲线的数值与单位是否相匹配，以及确定所述声波曲线中是否具有异常值，以及确定所述声波曲线是否覆盖所述井的全井段；当所述井的声波曲线的数值与单位相匹配、所述声波曲线中无异常值且所述声波曲线覆盖所述井的全井段时，确定所述井的检查结果为合格。
9.根据权利要求6所述的装置，其特征在于，所述第四确定模块，还被配置为根据所述声波曲线确定每个采样点对应的声波曲线值；将所述第一采样点对应的第一声波曲线值与第二采样点对应的第二声波曲线值的平均值作为第一采样点的层速度，所述第二采样点为在计算方向上距离所述第一采样点最近的采样点；根据所述第二采样点对应的第二声波曲线值与第三采样点对应的第三声波曲线值确定第二采样点的层速度，直到确定所述待计算井内每个采样点的层速度，所述第三采样点为在计算方向上距离所述第二采样点最近的采样点。
10.根据权利要求6所述的装置，其特征在于，所述第四确定模块，还被配置为确定第二采样点和所述第一采样点之间的第一测井采样深度差，所述第一测井采样深度差为所述第一深度和所述第二采样点的第二深度之间的差值，所述第二采样点为在计算方向上距离所述第一采样点最近的采样点；根据所述第一双程旅行时间、所述第一测井采样深度差和所述第二采样点的第二层速度，通过以下公式一，确定所述第二采样点的第二双程旅行时间；公式一：T(2)＝T(1)+2×△d(1)/VINT(2)；其中，T(2)表示第二采样点的第二双程旅行时间，T(1)表示第一采样点的双程旅行时间，△d(1)表示第一测井采样深度差，VINT(2)表示第二层速度；根据所述第二双程旅行时间以及所述第二深度，确定所述待计算井内第三采样点的第三双程旅行时间，直到确定出所述待计算井内每个采样点的双程旅行时间为止，所述第三采样点为在计算方向上距离所述第二采样点最近的采样点。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910164555.2A CN109932749B (zh) | 2019-03-05 | 2019-03-05 | 一种井震标定方法、装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910164555.2A CN109932749B (zh) | 2019-03-05 | 2019-03-05 | 一种井震标定方法、装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109932749A CN109932749A (zh) | 2019-06-25 |
| CN109932749B true CN109932749B (zh) | 2020-12-01 |
Family
ID=66986305
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910164555.2A Active CN109932749B (zh) | 2019-03-05 | 2019-03-05 | 一种井震标定方法、装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109932749B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110579806B (zh) * | 2019-09-11 | 2022-03-22 | 大庆油田有限责任公司 | 一种密井网条件下快速的井震标定方法 |
| CN110673209B (zh) * | 2019-10-13 | 2021-06-04 | 东北石油大学 | 一种井震标定方法 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103576199B (zh) * | 2012-07-19 | 2016-03-09 | 中国石油天然气集团公司 | 一种地震和测井时深关系自动转换的方法 |
| CN102967881B (zh) * | 2012-11-12 | 2015-07-08 | 中国石油天然气股份有限公司 | 用于地震资料解释的地质小层数据深时转换方法及装置 |
| CN103245971B (zh) * | 2013-04-25 | 2015-12-02 | 中国石油天然气股份有限公司 | 井中断点引导的地震小断层解释方法及装置 |
| CN104656142B (zh) * | 2013-11-19 | 2017-06-06 | 中国石油天然气集团公司 | 一种利用垂直地震剖面与测井联合的地震层位标定方法 |
| EP3092515B1 (en) * | 2014-03-11 | 2019-04-24 | Landmark Graphics Corporation | Validation of depth-depth curves using time-seismic depth interval velocity |
| CN105044773B (zh) * | 2015-08-24 | 2017-08-25 | 中国石油集团东方地球物理勘探有限责任公司 | 转换波井震标定方法及装置 |
| CN105759312B (zh) * | 2016-03-01 | 2018-10-16 | 中国海洋石油集团有限公司 | 一种储层特征法井震标定方法 |
| CN109283577B (zh) * | 2017-07-20 | 2023-02-28 | 中国石油化工股份有限公司 | 一种地震层位标定方法 |
| CN107831542B (zh) * | 2017-07-24 | 2019-05-31 | 中国石油化工股份有限公司 | Ddw高精度深度域井震匹配方法 |
-
2019
- 2019-03-05 CN CN201910164555.2A patent/CN109932749B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109932749A (zh) | 2019-06-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5835882A (en) | Method for determining barriers to reservoir flow | |
| US5835883A (en) | Method for determining distribution of reservoir permeability, porosity and pseudo relative permeability | |
| EP2846175B1 (en) | Seismic survey analysis | |
| US5838634A (en) | Method of generating 3-D geologic models incorporating geologic and geophysical constraints | |
| US4969130A (en) | System for monitoring the changes in fluid content of a petroleum reservoir | |
| WO1998034190A9 (en) | Method for determining distribution of reservoir permeability, porosity and pseudo relative permeability | |
| CN111596364B (zh) | 基于高精度层序地层格架的地震沉积微相组合分析方法 | |
| CN109870719B (zh) | 一种碳酸盐岩致密薄储层的井位布设方法、装置及系统 | |
| CN1882852A (zh) | 储层区的延时地震勘测 | |
| US7280952B2 (en) | Well planning using seismic coherence | |
| CN109932749B (zh) | 一种井震标定方法、装置 | |
| CN112363226A (zh) | 一种非常规油气有利区地球物理预测方法 | |
| US6754587B1 (en) | Method for processing seismic data | |
| CN112505754B (zh) | 基于高精度层序格架模型的井震协同划分沉积微相的方法 | |
| CN112578475B (zh) | 基于数据挖掘的致密储层双甜点识别方法 | |
| CN111624649B (zh) | 一种零偏移距vsp建立横向变速层速度模型的方法和装置 | |
| CN114861515A (zh) | 层速度数据体的计算方法、装置、设备及介质 | |
| Aziz et al. | The study of OpenDtect seismic data interpretation and visualization package in relation to seismic interpretation and visualization models | |
| CN115453627A (zh) | 一种无需模式识别的小生境粒子群面波反演方法 | |
| CN110795513B (zh) | 河流相源储异位型致密油气甜点区分布的预测方法 | |
| CN117236172B (zh) | 一种针对火成岩的深度学习低频模型构建方法及系统 | |
| CN117250658B (zh) | 建立研究区的地震数据集的方法 | |
| CN113550734B (zh) | 页岩气水平井地震随钻跟踪方法及装置 | |
| Smith et al. | The Chinguetti deepwater turbidite field, Mauritania: reserve estimation and field development using uncertainty management and experimental designs for multiple scenario 3D models | |
| CN113064202A (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 |