CN116859466B - 井震时深批量标定方法、装置、电子设备及介质 - Google Patents
井震时深批量标定方法、装置、电子设备及介质 Download PDFInfo
- Publication number
- CN116859466B CN116859466B CN202210311472.3A CN202210311472A CN116859466B CN 116859466 B CN116859466 B CN 116859466B CN 202210311472 A CN202210311472 A CN 202210311472A CN 116859466 B CN116859466 B CN 116859466B
- Authority
- CN
- China
- Prior art keywords
- time
- well
- depth
- nonstandard
- standard
- 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 34
- 230000006870 function Effects 0.000 claims description 45
- 238000004364 calculation method Methods 0.000 claims description 22
- 238000005070 sampling Methods 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 6
- 230000008021 deposition Effects 0.000 claims description 4
- 230000035939 shock Effects 0.000 claims description 4
- 230000002194 synthesizing effect Effects 0.000 claims description 4
- 238000013517 stratification Methods 0.000 abstract description 13
- 238000010586 diagram Methods 0.000 description 14
- 239000004576 sand Substances 0.000 description 8
- 238000012937 correction Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 4
- 239000003245 coal Substances 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000011161 development Methods 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
- 208000019888 Circadian rhythm sleep disease Diseases 0.000 description 2
- 208000001456 Jet Lag Syndrome Diseases 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000032798 delamination Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 208000033915 jet lag type circadian rhythm sleep disease Diseases 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 239000007787 solid Substances 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
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
本申请公开了一种井震时深批量标定方法、装置、电子设备及介质。该方法可以包括：确定标准井，进而获得标准井的时深关系；根据标准井的时深关系与非标准井的深度，获得非标准井的初始时间；根据地震层位的解释结果与非标准井的深度，获得非标准井的实际时间；计算非标准井的初始时间与实际时间的时差；根据非标准井的深度与初始时间以及时差，计算所有非标准井的时间，实现井震时深批量标定。本发明不仅能够实现有声波曲线井的高效率批量标定，以及对无声波与密度曲线的井进行批量标定，同时还可以快速检查地质分层的合理性。
Description
技术领域
本发明涉及油气地球物理技术领域，更具体地，涉及一种井震时深批量标定方法、装置、电子设备及介质。
背景技术
钻前储层预测工作正逐步从传统意义上的前期勘探向开发目标延伸，“地震-测井-地质”一体化综合研究的形势发展要求储层预测走向精细化。在目标解释与预测中，合成地震记录是将井震相结合的一条纽带，它最终将抽象地震数据与实际地质模型连接起来，为地震资料解释与储层预测的必要环节。
针对目的层的地震信息提取均离不开对测井信息的标定。通过高精度的合成地震记录，赋予地震剖面上反射波同相轴确切的地质意义，可为精细储层描述打下坚实的基础。为了获取高品质的合成记录，不同学者分别从理论和资料处理等不同角度分析了影响合成地震记录的主要因素，并在原有基础上加以改进，从而形成相应的方法。总体来讲，可概括为三种：(1)对于目标层，通常是根据地震波组特征进行相对标定，采用人机交互式的拉伸或压缩方法进行；(2)对于大尺度层或整个地层标定，为获取完整大间距的时深关系，一般采用多时窗变子波的标定方法；(3)针对全井段的局部储层高精度特殊标定，同时运用了多时窗变子波标定方法和声波时差进行均匀补偿方法。
以上这些依靠密度与声波求取反射系数，并通过子波褶积的得到井旁合成地震记录的单井标定方法，在客观上都提高了合成记录标定的精度，基本满足了岩性圈闭勘探的需求，在油气勘探中发挥着重要作用。但据目前“地震-测井-地质”一体化综合研究精细化的形势发展与预测效率要求，存在无法提高单井标定效率、无法利用无密度与声波的井对地震波进行精细解释与地质解剖、无法克服地表高程变化与沉积变化对直接时深拷贝的影响。
因此，有必要开发一种井震时深批量标定方法、装置、电子设备及介质。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种井震时深批量标定方法、装置、电子设备及介质，其不仅能够实现有声波曲线井的高效率批量标定，以及对无声波与密度曲线的井进行批量标定，同时还可以快速检查地质分层的合理性。
第一方面，本公开实施例提供了一种井震时深批量标定方法，包括：
确定标准井，进而获得标准井的时深关系；
根据所述标准井的时深关系与所述非标准井的深度，获得所述非标准井的初始时间；
根据地震层位的解释结果与所述非标准井的深度，获得所述非标准井的实际时间；
计算所述非标准井的初始时间与实际时间的时差；
根据所述非标准井的深度与初始时间以及所述时差，计算所有非标准井的时间，实现井震时深批量标定。
优选地，根据所述标准井的时深关系与所述非标准井的深度，获得所述非标准井的初始时间包括：
根据所述标准井的时深关系与所述非标准井的深度，确定所述非标准井的深度对应的所述标准井的时间；
以所述标准井的时间作为所述非标准井的深度对应的初始时间，获得所述非标准井的伪时间函数。
优选地，所述非标准井的伪时间函数为：
WellPicki(t,d)＝(Ti,di,xi,yi) (1)
其中，T为非标准井的初始时间，d为非标准井的深度，xi、yi为非标准井的纵横坐标。
优选地，根据地震层位的解释结果与所述非标准井的深度，获得所述非标准井的实际时间包括：
根据地震层位的解释结果获得地震层位的解释结果时间域的空间函数；
将所述非标准井的深度代入地震层位的解释结果时间域的空间函数，获得所述非标准井的实际时间。
优选地，地震层位的解释结果时间域的空间函数为：
H(x,y,t)＝(xi,yj,tk) (2)
其中，xi,yj,tk分别代表层位点在不同位置横坐标、纵坐标以及时间值。
优选地，通过公式(3)计算所述非标准井的初始时间与实际时间的时差：
SpickiΔti＝WellPicki(Ti，di，xi，yi)-H(xo，yp，tq) (3)
其中，SpickiΔti为时差。
优选地，通过公式(4)计算所有非标准井的时间：
其中，为非标准井的时间。
作为本公开实施例的一种具体实现方式，
第二方面，本公开实施例还提供了一种井震时深批量标定装置，包括：
标准井计算模块，确定标准井，进而获得标准井的时深关系；
初始时间计算模块，根据所述标准井的时深关系与所述非标准井的深度，获得所述非标准井的初始时间；
实际时间计算模块，根据地震层位的解释结果与所述非标准井的深度，获得所述非标准井的实际时间；
时差计算模块，计算所述非标准井的初始时间与实际时间的时差；
批量标定模块，根据所述非标准井的深度与初始时间以及所述时差，计算所有非标准井的时间，实现井震时深批量标定。
优选地，根据所述标准井的时深关系与所述非标准井的深度，获得所述非标准井的初始时间包括：
根据所述标准井的时深关系与所述非标准井的深度，确定所述非标准井的深度对应的所述标准井的时间；
以所述标准井的时间作为所述非标准井的深度对应的初始时间，获得所述非标准井的伪时间函数。
优选地，所述非标准井的伪时间函数为：
WellPicki(t,d)＝(Ti,di,xi,yi) (1)
其中，T为非标准井的初始时间，d为非标准井的深度，xi、yi为非标准井的纵横坐标。
优选地，根据地震层位的解释结果与所述非标准井的深度，获得所述非标准井的实际时间包括：
根据地震层位的解释结果获得地震层位的解释结果时间域的空间函数；
将所述非标准井的深度代入地震层位的解释结果时间域的空间函数，获得所述非标准井的实际时间。
优选地，地震层位的解释结果时间域的空间函数为：
H(x,y,t)＝(xi,yj,tk) (2)
其中，xi,yj,tk分别代表层位点在不同位置横坐标、纵坐标以及时间值。
优选地，通过公式(3)计算所述非标准井的初始时间与实际时间的时差：
SpickiΔti＝WellPicki(Ti，di，xi，yi)-H(xo，yp，tq) (3)
其中，SpickiΔti为时差。
优选地，通过公式(4)计算所有非标准井的时间：
其中，为非标准井的时间。
第三方面，本公开实施例还提供了一种电子设备，该电子设备包括：
存储器，存储有可执行指令；
处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的井震时深批量标定方法。
第四方面，本公开实施例还提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的井震时深批量标定方法。
其有益效果在于：针对目前在高密度井区单井标定效率低，以及在地震解释与储层预测中无法有效利用无声波曲线井标定的缺陷，本发明不仅能够实现有声波曲线井的高效率批量标定，以及对无声波与密度曲线的井进行批量标定，同时还可以快速检查地质分层的合理性。
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1示出了根据本发明的一个实施例的井震时深批量标定方法的步骤的流程图。
图2示出了根据本发明的一个实施例的部分井无声波曲线的初始时深关系标定的示意图。
图3示出了根据本发明的一个实施例的部分井无声波曲线批量校正后时深关系标定的示意图。
图4示出了根据本发明的一个实施例的标准井位置以及声波时差井分布的示意图。
图5示出了根据本发明的一个实施例的无声波时差井分布的示意图。
图6示出了根据本发明的一个实施例的标定结果对比的示意图。
图7示出了根据本发明的一个实施例的标准地质分层与解释层位时差直方图。
图8a和图8b分别示出了根据本发明的一个实施例的伪时深关系与批量校正后的时深关系图。
图9示出了根据本发明的一个实施例的地质分层可靠性分析的示意图。
图10示出了根据本发明的一个实施例的一种井震时深批量标定装置的框图。
附图标记说明：
201、标准井计算模块；202、初始时间计算模块；203、实际时间计算模块；204、时差计算模块；205、批量标定模块。
具体实施方式
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
本发明提供一种井震时深批量标定方法，包括：
确定标准井，进而获得标准井的时深关系；
根据标准井的时深关系与非标准井的深度，获得非标准井的初始时间；
根据地震层位的解释结果与非标准井的深度，获得非标准井的实际时间；
计算非标准井的初始时间与实际时间的时差；
根据非标准井的深度与初始时间以及时差，计算所有非标准井的时间，实现井震时深批量标定。
在一个示例中，根据标准井的时深关系与非标准井的深度，获得非标准井的初始时间包括：
根据标准井的时深关系与非标准井的深度，确定非标准井的深度对应的标准井的时间；
以标准井的时间作为非标准井的深度对应的初始时间，获得非标准井的伪时间函数。
在一个示例中，非标准井的伪时间函数为：
WellPicki(t,d)＝(Ti,di,xi,yi) (1)
其中，T为非标准井的初始时间，d为非标准井的深度，xi、yi为非标准井的纵横坐标。
在一个示例中，根据地震层位的解释结果与非标准井的深度，获得非标准井的实际时间包括：
根据地震层位的解释结果获得地震层位的解释结果时间域的空间函数；
将非标准井的深度代入地震层位的解释结果时间域的空间函数，获得非标准井的实际时间。
在一个示例中，地震层位的解释结果时间域的空间函数为：
H(x,y,t)＝(xi,yj,tk) (2)
其中，xi,yj,tk分别代表层位点在不同位置横坐标、纵坐标以及时间值。
在一个示例中，通过公式(3)计算非标准井的初始时间与实际时间的时差：
SpickiΔti＝WellPicki(Ti，di，xi，yi)-H(xo，yp，tq) (3)
其中，SpickiΔti为时差。
在一个示例中，通过公式(4)计算所有非标准井的时间：
其中，为非标准井的时间。
具体地，在三维地震覆盖区选取无扩径影响的密度与声波曲线的标准井，或对测井数据进行环境校正。
利用密度与声波曲线计算反射系数，确定子波主频类型，计算反射系数，进行VTI垂向各项异性标定，获取高精度的时深关系。若沉积厚度变化较大，可对选不同位置的标准井对储层进行分区标定，并获取每个分区储层段上下的时深关系；
设子波为w(t)：
w(t)＝[1-2(πf0t)2]exp[-(πf0t)2] (5)
则标准井的反射系数序列函数分别为：
其中，为标准井在深度i点的双程旅行时，ρi、Δti分别为在深度i点密度与声波时差瞬时值，而Δh为测井深度的采样率,i＝(1,2,3....m)。
合成地震记录卷积公式：
S(t)＝R(t)*W(t) (7)
将公式(5)与(6)代入公式(7)则有：
S(t)＝Rw(ti)*[1-2(πf0ti)2]exp[-(πf0ti)2] (8)
进行VTI标定后的反射系数序列为：
其中，为标准井在深度i点的VTI垂向各项异性衰减双程旅行时，ρi、Δti分别为在深度j点密度与声波时差瞬时值，而Δh为测井深度的采样率，β为VTI垂向各项异性衰减系数。
据合成地震记录卷积公式S(t)＝R(t)*W(t)，将公式(9)与(4)代入公式(7)则有：
S(Ti)＝Rw(Ti)*[1-2(πf0Ti)2]exp[-(πf0Ti)2] (10)
根据公式(10)与地震采集记录标定的结果，获取准确的时深对应关系，并将标准井时深拷贝至工区内所有井，对应的时间与深度序列、及相应的纵横坐标的时深函数，即标准井的时深关系为：
根据标准井的时深关系可直接获取所有井的标准地质分层的时间T，深度d，纵横坐标xi与yi的函数关系表示为公式(1)。
根据标准井的井震标定结果，选取储层附近强反射层或不整合面的地震反射同相轴，并对其自动追踪解释，高效率地获取反射界面的层位解释结果。设地震层位的解释结果时间域的空间函数为公式(2)。
据每口井的标准地质分层函数(1)与层位函数(2)，在同一坐标点下，分别获取地质分层时间值与地震层位时间的差值：
WellPicki(t,d)∩H(x,y,t)＝(Ti,di,xi,yi)∩(xo,yp,tq) (12)
即当xi＝xo且yi＝yq时，每口井的时差为公式(3)。
将不同井的时间差异量SpickiΔti，各自以常量的方式应用标准井的时深关系的时间序列中，同时保持各自测井深度不变。仅更新所有井各自的时间序列，则批量计算表达式可表示为公式(4)。
本发明还提供一种井震时深批量标定装置，包括：
标准井计算模块，确定标准井，进而获得标准井的时深关系；
初始时间计算模块，根据标准井的时深关系与非标准井的深度，获得非标准井的初始时间；
实际时间计算模块，根据地震层位的解释结果与非标准井的深度，获得非标准井的实际时间；
时差计算模块，计算非标准井的初始时间与实际时间的时差；
批量标定模块，根据非标准井的深度与初始时间以及时差，计算所有非标准井的时间，实现井震时深批量标定。
在一个示例中，根据标准井的时深关系与非标准井的深度，获得非标准井的初始时间包括：
根据标准井的时深关系与非标准井的深度，确定非标准井的深度对应的标准井的时间；
以标准井的时间作为非标准井的深度对应的初始时间，获得非标准井的伪时间函数。
在一个示例中，非标准井的伪时间函数为：
WellPicki(t,d)＝(Ti,di,xi,yi) (1)
其中，T为非标准井的初始时间，d为非标准井的深度，xi、yi为非标准井的纵横坐标。
在一个示例中，根据地震层位的解释结果与非标准井的深度，获得非标准井的实际时间包括：
根据地震层位的解释结果获得地震层位的解释结果时间域的空间函数；
将非标准井的深度代入地震层位的解释结果时间域的空间函数，获得非标准井的实际时间。
在一个示例中，地震层位的解释结果时间域的空间函数为：
H(x,y,t)＝(xi,yj,tk) (2)
其中，xi,yj,tk分别代表层位点在不同位置横坐标、纵坐标以及时间值。
在一个示例中，通过公式(3)计算非标准井的初始时间与实际时间的时差：
SpickiΔti＝WellPicki(Ti，di，xi，yi)-H(xo，yp，tq) (3)
其中，SpickiΔti为时差。
在一个示例中，通过公式(4)计算所有非标准井的时间：
其中，为非标准井的时间。
具体地，在三维地震覆盖区选取无扩径影响的密度与声波曲线的标准井，或对测井数据进行环境校正。
利用密度与声波曲线计算反射系数，确定子波主频类型，计算反射系数，进行VTI垂向各项异性标定，获取高精度的时深关系。若沉积厚度变化较大，可对选不同位置的标准井对储层进行分区标定，并获取每个分区储层段上下的时深关系；
设子波为w(t)：
w(t)＝[1-2(πf0t)2]exp[-(πf0t)2] (5)
则标准井的反射系数序列函数分别为：
其中，为标准井在深度i点的双程旅行时，ρi、Δti分别为在深度i点密度与声波时差瞬时值，而Δh为测井深度的采样率,i＝(1,2,3....m)。
合成地震记录卷积公式：
S(t)＝R(t)*W(t) (7)
将公式(5)与(6)代入公式(7)则有：
S(t)＝Rw(ti)*[1-2(πf0ti)2]exp[-(πf0ti)2] (8)
进行VTI标定后的反射系数序列为：
其中，为标准井在深度i点的VTI垂向各项异性衰减双程旅行时，ρi、Δti分别为在深度j点密度与声波时差瞬时值，而Δh为测井深度的采样率，β为VTI垂向各项异性衰减系数。
据合成地震记录卷积公式S(t)＝R(t)*W(t)，将公式(9)与(4)代入公式(7)则有：
S(Ti)＝Rw(Ti)*[1-2(πf0Ti)2]exp[-(πf0Ti)2] (10)
根据公式(10)与地震采集记录标定的结果，获取准确的时深对应关系，并将标准井时深拷贝至工区内所有井，对应的时间与深度序列、及相应的纵横坐标的时深函数，即标准井的时深关系为：
根据标准井的时深关系可直接获取所有井的标准地质分层的时间T，深度d，纵横坐标xi与yi的函数关系表示为公式(1)。
根据标准井的井震标定结果，选取储层附近强反射层或不整合面的地震反射同相轴，并对其自动追踪解释，高效率地获取反射界面的层位解释结果。设地震层位的解释结果时间域的空间函数为公式(2)。
据每口井的标准地质分层函数(1)与层位函数(2)，在同一坐标点下，分别获取地质分层时间值与地震层位时间的差值：
WellPicki(t,d)∩H(x,y,t)＝(Ti,di,xi,yi)∩(xo,yp,tq) (12)
即当xi＝xo且yi＝yq时，每口井的时差为公式(3)。
将不同井的时间差异量SpickiΔti，各自以常量的方式应用标准井的时深关系的时间序列中，同时保持各自测井深度不变。仅更新所有井各自的时间序列，则批量计算表达式可表示为公式(4)。
本发明还提供一种电子设备，电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述的井震时深批量标定方法。
本发明还提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述的井震时深批量标定方法。
为便于理解本发明实施例的方案及其效果，以下给出四个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
实施例1
图1示出了根据本发明的井震时深批量标定方法的步骤的流程图。
如图1所示，该井震时深批量标定方法包括：步骤101，确定标准井，进而获得标准井的时深关系；步骤102，根据标准井的时深关系与非标准井的深度，获得非标准井的初始时间；步骤103，根据地震层位的解释结果与非标准井的深度，获得非标准井的实际时间；步骤104，计算非标准井的初始时间与实际时间的时差；步骤105，根据非标准井的深度与初始时间以及时差，计算所有非标准井的时间，实现井震时深批量标定。
图2示出了根据本发明的一个实施例的部分井无声波曲线的初始时深关系标定的示意图。
利用优选的标准井W393获取的时深关系，将时深关系直接批量应用至其他无声波与密度曲线的井，获取的结果如图1所示，连井剖面中仅有W392的地质分层(Coal)与其对应的反射层位相吻合，误差较小。其他无声波密度曲线的井地质反射层与煤层反射层位的标定匹配都有非常大误差。其次，储层sand的标定误差也通常会引起储层地震反射的层位解释问题。
图3示出了根据本发明的一个实施例的部分井无声波曲线批量校正后时深关系标定的示意图。
利用优选的标准井W393获取的时深关系，将时深关系直接批量应用至其他无声波与密度曲线的井后，利用本文批量校正的时深算法，利用新的时深关系得到的标定结果如图2所示，Coal煤层反射界面与其对应的分层界面一一对应，无声波密度曲线的井地质反射层与解释层位的标定匹配误差为零。其次，下部的储层分层与厚度(sand)也与其地震反射的强弱一一对应，在砂体解释过程中可有效提高反射界面的解释精度。
图4示出了根据本发明的一个实施例的标准井位置以及声波时差井分布的示意图。
图5示出了根据本发明的一个实施例的无声波时差井分布的示意图。
以南美某工区为例，工区面积900平方公里，三维覆盖区内总计417口井，如图4与图5所示，井控密度约2口/km2；其中，具有声波时差与密度的井共计27口，地震解释与储层预测的井控密度为3口/100km2，远不能满足精细解释与预测的要求。为提高砂体解释精度，采用了本发明所提出的方法。
图6示出了根据本发明的一个实施例的标定结果对比的示意图。
选取如图4所示的标准井W393，据其中的密度与声波求取的时间域反射系数如图6所示。子波褶积的传统标定与VTI标定如图6中a与b所示，两者对比可明显看出图6中b的VTI标定合成记录比传统标定更好，与井旁的地震道的相关关系更高。
图7示出了根据本发明的一个实施例的标准地质分层与解释层位时差直方图。
表1为部分井标准地质分层的地震反射时间与其伪时深分层时间的时间差值表，在时域地质分层的时深关系中，时间误差在-30ms～30ms之间，其中图7为误差表的直方图显示，即采用统一的时深关系，不能确保标定的准确性。
表1
图8a和图8b分别示出了根据本发明的一个实施例的伪时深关系与批量校正后的时深关系图。
表2为W150伪时深校正前后的对比表，其中更新前后的时差为一常量15.6ms，同时保持TVD，SSTVD或MD的值不变，也就是仅改变时间序列。图8a和图8b为部分井校正前后的时深关系对比图。
表2
| 井名 | 初始合成记录时间 | 更新后合成记录时间 | 时间差 | TVDSS | TVD |
| W150 | 601.00 | 616.60 | 15.60 | 336.87 | 610.13 |
| W150 | 603.00 | 618.60 | 15.60 | 333.49 | 613.51 |
| W150 | 605.00 | 620.60 | 15.60 | 330.04 | 616.96 |
| W150 | 607.00 | 622.60 | 15.60 | 326.63 | 620.37 |
| W150 | 609.00 | 624.60 | 15.60 | 323.35 | 623.65 |
| W150 | 611.00 | 626.60 | 15.60 | 319.93 | 627.07 |
| W150 | 613.00 | 628.60 | 15.60 | 316.52 | 630.48 |
| W150 | 615.00 | 630.60 | 15.60 | 313.25 | 633.75 |
| W150 | 617.00 | 632.60 | 15.60 | 310.09 | 636.91 |
| W150 | 619.00 | 634.60 | 15.60 | 306.74 | 640.26 |
| W150 | 621.00 | 636.60 | 15.60 | 303.74 | 643.26 |
| W150 | 623.00 | 638.60 | 15.60 | 300.72 | 646.28 |
| W150 | 625.00 | 640.60 | 15.60 | 297.91 | 649.09 |
| W150 | 627.00 | 642.60 | 15.60 | 295.10 | 651.91 |
| W150 | 629.00 | 644.60 | 15.60 | 292.28 | 654.72 |
| W150 | 631.00 | 646.60 | 15.60 | 289.16 | 657.84 |
图9示出了根据本发明的一个实施例的地质分层可靠性分析的示意图。
由图9的GR放射性特征可见W150在RCCK分层附近为单个砂体的反射，而由W151至W154在RCCK分层附近为一砂组反射，顶部为一薄砂，底部为一厚砂；从岩性反射界面而言，由此确定真实分层的位置应该位于薄砂的顶部，所以根据标定的结果可快速有效地检验地质分层的可靠性，并及时进行分层调整，同时深化地质认识。
实施例2
图10示出了根据本发明的一个实施例的一种井震时深批量标定装置的框图。
如图10所示，该井震时深批量标定装置，包括：
标准井计算模块201，确定标准井，进而获得标准井的时深关系；
初始时间计算模块202，根据标准井的时深关系与非标准井的深度，获得非标准井的初始时间；
实际时间计算模块203，根据地震层位的解释结果与非标准井的深度，获得非标准井的实际时间；
时差计算模块204，计算非标准井的初始时间与实际时间的时差；
批量标定模块205，根据非标准井的深度与初始时间以及时差，计算所有非标准井的时间，实现井震时深批量标定。
作为可选方案，根据标准井的时深关系与非标准井的深度，获得非标准井的初始时间包括：
根据标准井的时深关系与非标准井的深度，确定非标准井的深度对应的标准井的时间；
以标准井的时间作为非标准井的深度对应的初始时间，获得非标准井的伪时间函数。
作为可选方案，非标准井的伪时间函数为：
WellPicki(t,d)＝(Ti,di,xi,yi) (1)
其中，T为非标准井的初始时间，d为非标准井的深度，xi、yi为非标准井的纵横坐标。
作为可选方案，根据地震层位的解释结果与非标准井的深度，获得非标准井的实际时间包括：
根据地震层位的解释结果获得地震层位的解释结果时间域的空间函数；
将非标准井的深度代入地震层位的解释结果时间域的空间函数，获得非标准井的实际时间。
作为可选方案，地震层位的解释结果时间域的空间函数为：
H(x,y,t)＝(xi,yj,tk) (2)
其中，xi,yj,tk分别代表层位点在不同位置横坐标、纵坐标以及时间值。
作为可选方案，通过公式(3)计算非标准井的初始时间与实际时间的时差：
SpickiΔti＝WellPicki(Ti，di，xi，yi)-H(xo，yp，tq) (3)
其中，SpickiΔti为时差。
作为可选方案，通过公式(4)计算所有非标准井的时间：
其中，为非标准井的时间。
实施例3
本公开提供一种电子设备包括，该电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述井震时深批量标定方法。
根据本公开实施例的电子设备包括存储器和处理器。
该存储器用于存储非暂时性计算机可读指令。具体地，存储器可以包括一个或多个计算机程序产品，该计算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(RAM)和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(ROM)、硬盘、闪存等。
该处理器可以是中央处理单元(CPU)或者具有数据处理能力和/或指令执行能力的其它形式的处理单元，并且可以控制电子设备中的其它组件以执行期望的功能。在本公开的一个实施例中，该处理器用于运行该存储器中存储的该计算机可读指令。
本领域技术人员应能理解，为了解决如何获得良好用户体验效果的技术问题，本实施例中也可以包括诸如通信总线、接口等公知的结构，这些公知的结构也应包含在本公开的保护范围之内。
有关本实施例的详细说明可以参考前述各实施例中的相应说明，在此不再赘述。
实施例4
本公开实施例提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的井震时深批量标定方法。
根据本公开实施例的计算机可读存储介质，其上存储有非暂时性计算机可读指令。当该非暂时性计算机可读指令由处理器运行时，执行前述的本公开各实施例方法的全部或部分步骤。
上述计算机可读存储介质包括但不限于：光存储介质(例如：CD－ROM和DVD)、磁光存储介质(例如：MO)、磁存储介质(例如：磁带或移动硬盘)、具有内置的可重写非易失性存储器的媒体(例如：存储卡)和具有内置ROM的媒体(例如：ROM盒)。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (7)
1.一种井震时深批量标定方法，其特征在于，包括：
确定标准井，进而获得标准井的时深关系；
根据所述标准井的时深关系与非标准井的深度，获得所述非标准井的初始时间；
根据地震层位的解释结果与所述非标准井的深度，获得所述非标准井的实际时间；
计算所述非标准井的初始时间与实际时间的时差；
根据所述非标准井的深度与初始时间以及所述时差，计算所有非标准井的时间，实现井震时深批量标定；
其中，根据所述标准井的时深关系与所述非标准井的深度，获得所述非标准井的初始时间包括：
根据所述标准井的时深关系与所述非标准井的深度，确定所述非标准井的深度对应的所述标准井的时间；
以所述标准井的时间作为所述非标准井的深度对应的初始时间，获得所述非标准井的伪时间函数；
其中，所述非标准井的伪时间函数为：
WellPicki(t,d)＝(Ti,di,xi,yi) (1)
其中，T为非标准井的初始时间，d为非标准井的深度，xi、yi为非标准井的纵横坐标；
通过公式(4)计算所有非标准井的时间：
其中，为非标准井的时间，SpickiΔti为时差；
其中，确定标准井，进而获得标准井的时深关系包括：利用密度与声波曲线计算反射系数，确定子波主频类型，计算反射系数，进行VTI垂向各项异性标定，获取时深关系；若沉积厚度变化大，对不同位置的标准井对储层进行分区标定，并获取每个分区储层段上下的时深关系；
设子波为w(t)：
w(t)＝[1-2(πf0t)2]exp[-(πf0t)2] (5)
则标准井的反射系数序列函数分别为：
其中，为标准井在深度i点的双程旅行时，ρi、Δti分别为在深度i点密度与声波时差瞬时值，而Δh为测井深度的采样率，i＝(1,2,3....m)；
合成地震记录卷积公式：
S(t)＝R(t)*W(t) (7)
将公式(5)与(6)代入公式(7)则有：
S(t)＝Rw(ti)*[1-2(πf0ti)2]exp[-(πf0ti)2] (8)
进行VTI标定后的反射系数序列为：
其中，为标准井在深度i点的VTI垂向各项异性衰减双程旅行时，ρi、Δti分别为在深度j点密度与声波时差瞬时值，而Δh为测井深度的采样率，β为VTI垂向各项异性衰减系数。
2.根据权利要求1所述的井震时深批量标定方法，其中，根据地震层位的解释结果与所述非标准井的深度，获得所述非标准井的实际时间包括：
根据地震层位的解释结果获得地震层位的解释结果时间域的空间函数；
将所述非标准井的深度代入地震层位的解释结果时间域的空间函数，获得所述非标准井的实际时间。
3.根据权利要求2所述的井震时深批量标定方法，其中，地震层位的解释结果时间域的空间函数为：
H(x,y,t)＝(xi,yj,tk) (2)
其中，xi,yj,tk分别代表层位点在不同位置横坐标、纵坐标以及时间值。
4.根据权利要求1所述的井震时深批量标定方法，其中，通过公式(3)计算所述非标准井的初始时间与实际时间的时差：
SpickiΔti＝WellPicki(Ti,di,xi,yi)-H(xo,yp,tq) (3)。
5.一种井震时深批量标定装置，其特征在于，包括：
标准井计算模块，确定标准井，进而获得标准井的时深关系；
初始时间计算模块，根据所述标准井的时深关系与非标准井的深度，获得所述非标准井的初始时间；
实际时间计算模块，根据地震层位的解释结果与所述非标准井的深度，获得所述非标准井的实际时间；
时差计算模块，计算所述非标准井的初始时间与实际时间的时差；
批量标定模块，根据所述非标准井的深度与初始时间以及所述时差，计算所有非标准井的时间，实现井震时深批量标定；
其中，根据所述标准井的时深关系与所述非标准井的深度，获得所述非标准井的初始时间包括：
根据所述标准井的时深关系与所述非标准井的深度，确定所述非标准井的深度对应的所述标准井的时间；
以所述标准井的时间作为所述非标准井的深度对应的初始时间，获得所述非标准井的伪时间函数；
其中，所述非标准井的伪时间函数为：
WellPicki(t,d)＝(Ti,di,xi,yi) (1)
其中，T为非标准井的初始时间，d为非标准井的深度，xi、yi为非标准井的纵横坐标；
通过公式(4)计算所有非标准井的时间：
其中，为非标准井的时间，SpickiΔti为时差；
其中，确定标准井，进而获得标准井的时深关系包括：利用密度与声波曲线计算反射系数，确定子波主频类型，计算反射系数，进行VTI垂向各项异性标定，获取时深关系；若沉积厚度变化大，对不同位置的标准井对储层进行分区标定，并获取每个分区储层段上下的时深关系；
设子波为w(t)：
w(t)＝[1-2(πf0t)2]exp[-(πf0t)2] (5)
则标准井的反射系数序列函数分别为：
其中，为标准井在深度i点的双程旅行时，ρi、Δti分别为在深度i点密度与声波时差瞬时值，而Δh为测井深度的采样率，i＝(1,2,3....m)；
合成地震记录卷积公式：
S(t)＝R(t)*W(t) (7)
将公式(5)与(6)代入公式(7)则有：
S(t)＝Rw(ti)*[1-2(πf0ti)2]exp[-(πf0ti)2] (8)
进行VTI标定后的反射系数序列为：
其中，为标准井在深度i点的VTI垂向各项异性衰减双程旅行时，ρi、Δti分别为在深度j点密度与声波时差瞬时值，而Δh为测井深度的采样率，β为VTI垂向各项异性衰减系数。
6.一种电子设备，其特征在于，所述电子设备包括：
存储器，存储有可执行指令；
处理器，所述处理器运行所述存储器中的所述可执行指令，以实现权利要求1-4中任一项所述的井震时深批量标定方法。
7.一种计算机可读存储介质，其特征在于，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现权利要求1-4中任一项所述的井震时深批量标定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210311472.3A CN116859466B (zh) | 2022-03-28 | 2022-03-28 | 井震时深批量标定方法、装置、电子设备及介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210311472.3A CN116859466B (zh) | 2022-03-28 | 2022-03-28 | 井震时深批量标定方法、装置、电子设备及介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116859466A CN116859466A (zh) | 2023-10-10 |
| CN116859466B true CN116859466B (zh) | 2024-05-31 |
Family
ID=88229133
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210311472.3A Active CN116859466B (zh) | 2022-03-28 | 2022-03-28 | 井震时深批量标定方法、装置、电子设备及介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116859466B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102520448A (zh) * | 2012-01-06 | 2012-06-27 | 中国石油化工股份有限公司胜利油田分公司西部新区研究中心 | 年代地层旋回域井震标定系统 |
| CN103089232A (zh) * | 2011-10-27 | 2013-05-08 | 中国石油集团长城钻探工程有限公司 | 测井设备的检测和一致性校正方法 |
| CN107130955A (zh) * | 2016-02-26 | 2017-09-05 | 中国石油化工股份有限公司 | 井底流压的确定方法及储集体天然能量的确定方法 |
| CN113625360A (zh) * | 2020-05-08 | 2021-11-09 | 中国石油化工股份有限公司 | 微裂缝地层产量预测方法、系统、电子设备及介质 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011055359A1 (en) * | 2009-11-04 | 2011-05-12 | Acousticeye Ltd | Apparatus and method for measuring liquid level in a well |
-
2022
- 2022-03-28 CN CN202210311472.3A patent/CN116859466B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103089232A (zh) * | 2011-10-27 | 2013-05-08 | 中国石油集团长城钻探工程有限公司 | 测井设备的检测和一致性校正方法 |
| CN102520448A (zh) * | 2012-01-06 | 2012-06-27 | 中国石油化工股份有限公司胜利油田分公司西部新区研究中心 | 年代地层旋回域井震标定系统 |
| CN107130955A (zh) * | 2016-02-26 | 2017-09-05 | 中国石油化工股份有限公司 | 井底流压的确定方法及储集体天然能量的确定方法 |
| CN113625360A (zh) * | 2020-05-08 | 2021-11-09 | 中国石油化工股份有限公司 | 微裂缝地层产量预测方法、系统、电子设备及介质 |
Non-Patent Citations (1)
| Title |
|---|
| 鄂尔多斯盆地黄土塬三维地震技术突破及勘探开发效果;付锁堂，等;中国石油勘探;20200131;第25卷(第1期);67-77 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116859466A (zh) | 2023-10-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103454685B (zh) | 利用测井约束波阻抗反演预测砂体厚度的方法和装置 | |
| CN102937720B (zh) | 井控提高地震资料分辨率的方法 | |
| CN104516018B (zh) | 一种地球物理勘探中岩性约束下的孔隙度反演方法 | |
| CN109738945B (zh) | 一种利用叠前深度偏移成果直接生成构造图的方法 | |
| CN108931814A (zh) | 一种基于多属性融合的基岩裂缝预测的方法 | |
| CN102176054B (zh) | 近地表综合信息处理解释方法 | |
| CN111257926A (zh) | 一种利用老旧地震资料进行古河谷铀储层预测的方法 | |
| WO2008070596A1 (en) | Identification of fracture clusters in rock formations | |
| CN111123359B (zh) | 随钻测井与地层格架约束的井周地震成像探测方法及装置 | |
| CN111722284B (zh) | 一种基于道集数据建立速度深度模型的方法 | |
| CN111983683B (zh) | 一种少井条件下盐下湖相灰岩储层预测方法和系统 | |
| WO2018227034A1 (en) | Anisotropic nmo correction and its application to attenuate noises in vsp data | |
| CN105425299A (zh) | 确定地层裂缝分布的方法和装置 | |
| CN106842299A (zh) | 一种基于地震属性的裂缝定量化预测的方法 | |
| CN116856845A (zh) | 一种致密砂岩气藏水平井地质导向方法 | |
| CN116859466B (zh) | 井震时深批量标定方法、装置、电子设备及介质 | |
| Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
| CN112147700A (zh) | 速度异常区的低频模型构建方法及系统 | |
| CN116931089B (zh) | 基于厚度分区的时深标定方法、装置、电子设备及介质 | |
| CN107255832A (zh) | 一种地下构造的反演方法 | |
| CN109839659B (zh) | 一种对叠前深度偏移剖面进行迭代优化的方法 | |
| CN112180466B (zh) | 一种河道砂体识别方法及装置 | |
| CN113495293B (zh) | 油藏流体预测方法及装置 | |
| Nourollah et al. | Seal potential of shale sequences through seismic anisotropy: Case study from Exmouth Sub-basin, Australia | |
| CN110941029B (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 |