CN116931089B - 基于厚度分区的时深标定方法、装置、电子设备及介质 - Google Patents
基于厚度分区的时深标定方法、装置、电子设备及介质 Download PDFInfo
- Publication number
- CN116931089B CN116931089B CN202210325161.2A CN202210325161A CN116931089B CN 116931089 B CN116931089 B CN 116931089B CN 202210325161 A CN202210325161 A CN 202210325161A CN 116931089 B CN116931089 B CN 116931089B
- Authority
- CN
- China
- Prior art keywords
- time
- well
- depth
- thickness
- nonstandard
- 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
- 238000005192 partition Methods 0.000 title claims abstract description 106
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000006870 function Effects 0.000 claims description 31
- 238000000638 solvent extraction Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 8
- 239000002344 surface layer Substances 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 5
- 238000013316 zoning Methods 0.000 claims description 3
- 230000008021 deposition Effects 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 12
- 239000010410 layer Substances 0.000 description 10
- 238000013517 stratification Methods 0.000 description 8
- 230000035939 shock Effects 0.000 description 6
- 238000012937 correction Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000003786 synthesis reaction Methods 0.000 description 3
- 208000019888 Circadian rhythm sleep disease Diseases 0.000 description 2
- 208000001456 Jet Lag Syndrome Diseases 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 208000033915 jet lag type circadian rhythm sleep disease Diseases 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000002224 dissection 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
- 238000011160 research Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000004062 sedimentation Methods 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)针对全井段的局部储层高精度特殊标定,同时运用了多时窗变子波标定方法和声波时差进行均匀补偿方法。
以上这些方法随滚动勘探开发程度的提高,单位面积的钻井密度也越来越大,随钻井密度的提高,对井震闭合、解释与反演的精度提出更高的要求;同时,油公司为节约测井成本,往往仅对部分关键井测量密度与声波曲线;另外,目前的时深关系相互拷贝的方法,因受沉积厚度以及近地表高程的影响无法确保目的层地质分层与地震反射的一一对应关系。总之,据目前“地震-测井-地质”一体化综合研究精细化的形势发展与预测效率要求,存在无法提高单井标定效率、无法利用无密度与声波的井对地震波进行精细解释与地质解剖、无法克服地表高程变化与沉积变化对直接时深拷贝的影响。
因此,有必要开发一种基于厚度分区的时深标定方法、装置、电子设备及介质。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种基于厚度分区的时深标定方法、装置、电子设备及介质,其能够通过不同厚度分区内的标准井,克服沉积厚度对时深关系的影响,实现对有声波或无声波曲线井的高精度、批量高效标定。
第一方面,本公开实施例提供了一种基于厚度分区的时深标定方法,包括:
根据目标区域中标准井的井震标定结果进行自动追踪解释,获得地震层位顶、底的解释结果;
根据地震层位顶、底的解释结果计算时间厚度,进而根据所述时间厚度在目标区域进行区域划分,获得多个厚度分区;
确定每一个厚度分区的标准井,进而获得该标准井的时深关系;
在每一个厚度分区内,根据所述标准井的时深关系与所述非标准井的深度,计算所述非标准井的时间,实现基于该厚度分区的时深标定。
优选地,所述地震层位顶的解释结果为:
Htop(x,y,t)=(xi,yj,ti) (1)
所述地震层位底的解释结果为:
Hbot(x,y,t)=(xi,yj,ti) (2)
其中,Htop(x,y,t)为地震层位顶的解释结果,Hbot(x,y,t)为地震层位底的解释结果,xi,yj,ti分别代表层位点在不同位置横坐标、纵坐标以及时间值。
优选地,通过公式(3)计算时间厚度:
HD(xi,yi,Δti)=Hbot(xi,yi,ti)-Htop(xi,yi,ti) (3)
其中,HD(xi,yi,Δti)为时间厚度。
优选地,在每一个厚度分区内,根据所述标准井的时深关系与所述非标准井的深度,计算所述非标准井的时间:
在每一个厚度分区内,根据所述标准井的时深关系与所述非标准井的深度,确定所述非标准井的深度对应的所述标准井的时间;
以所述标准井的时间作为所述非标准井的深度对应的初始时间,获得所述非标准井的伪时间函数;
根据地震层位的解释结果与所述非标准井的深度,获得所述非标准井的实际时间;
计算所述非标准井的初始时间与实际时间的时差;
根据所述非标准井的深度与初始时间以及所述时差,计算所有非标准井的时间,实现基于厚度分区的时深标定。
优选地,所述非标准井的伪时间函数为:
其中,k代表分区序号,T为非标准井的初始时间,d为非标准井的深度,xi、yi为非标准井的纵横坐标。
优选地,根据地震层位的解释结果与所述非标准井的深度,获得所述非标准井的实际时间包括:
根据地震层位的解释结果获得地震层位的解释结果时间域的空间函数;
将所述非标准井的深度代入地震层位的解释结果时间域的空间函数,获得所述非标准井的实际时间。
优选地,通过公式(5)计算所有非标准井的时间:
其中,为非标准井的时间。
作为本公开实施例的一种具体实现方式,
第二方面,本公开实施例还提供了一种基于厚度分区的时深标定装置,包括:
追踪解释模块,根据目标区域中标准井的井震标定结果进行自动追踪解释,获得地震层位顶、底的解释结果;
厚度分区模块,根据地震层位顶、底的解释结果计算时间厚度,进而根据所述时间厚度在目标区域进行区域划分,获得多个厚度分区;
标准井计算模块,确定每一个厚度分区的标准井,进而获得该标准井的时深关系;
标定模块,在每一个厚度分区内,根据所述标准井的时深关系与所述非标准井的深度,计算所述非标准井的时间,实现基于该厚度分区的时深标定。
优选地,所述地震层位顶的解释结果为:
Htop(x,y,t)=(xi,yj,ti) (1)
所述地震层位底的解释结果为:
Hbot(x,y,t)=(xi,yj,ti) (2)
其中,Htop(x,y,t)为地震层位顶的解释结果,Hbot(x,y,t)为地震层位底的解释结果,xi,yj,ti分别代表层位点在不同位置横坐标、纵坐标以及时间值。
优选地,通过公式(3)计算时间厚度:
HD(xi,yi,Δti)=Hbot(xi,yi,ti)-Htop(xi,yi,ti) (3)
其中,HD(xi,yi,Δti)为时间厚度。
优选地,在每一个厚度分区内,根据所述标准井的时深关系与所述非标准井的深度,计算所述非标准井的时间:
在每一个厚度分区内,根据所述标准井的时深关系与所述非标准井的深度,确定所述非标准井的深度对应的所述标准井的时间;
以所述标准井的时间作为所述非标准井的深度对应的初始时间,获得所述非标准井的伪时间函数;
根据地震层位的解释结果与所述非标准井的深度,获得所述非标准井的实际时间;
计算所述非标准井的初始时间与实际时间的时差;
根据所述非标准井的深度与初始时间以及所述时差,计算所有非标准井的时间,实现基于厚度分区的时深标定。
优选地,所述非标准井的伪时间函数为:
其中,k代表分区序号,T为非标准井的初始时间,d为非标准井的深度,xi、yi为非标准井的纵横坐标。
优选地,根据地震层位的解释结果与所述非标准井的深度,获得所述非标准井的实际时间包括:
根据地震层位的解释结果获得地震层位的解释结果时间域的空间函数;
将所述非标准井的深度代入地震层位的解释结果时间域的空间函数,获得所述非标准井的实际时间。
优选地,通过公式(5)计算所有非标准井的时间:
其中,为非标准井的时间。
第三方面,本公开实施例还提供了一种电子设备,该电子设备包括:
存储器,存储有可执行指令;
处理器,所述处理器运行所述存储器中的所述可执行指令,以实现所述的基于厚度分区的时深标定方法。
第四方面,本公开实施例还提供了一种计算机可读存储介质,该计算机可读存储介质存储有计算机程序,该计算机程序被处理器执行时实现所述的基于厚度分区的时深标定方法。
其有益效果在于:本发明为克服标定层段内因沉积厚度对时深关系的影响,通过计算不同分区内标准井分层的时深与层位反射时间的时间差异,实现快速批量的井震标定,获取储层段更加精确的时深关系。
本发明的方法和装置具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的,或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述,这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的一个实施例的基于厚度分区的时深标定方法的步骤的流程图。
图2示出了根据本发明的一个实施例的不同分区标准井位置以及声波时差井分布的示意图。
图3示出了根据本发明的一个实施例的无声波时差井分布的示意图。
图4示出了根据本发明的一个实施例的标定结果对比的示意图。
图5示出了根据本发明的一个实施例的A区标准地质分层与解释层位时差直方图。
图6示出了根据本发明的一个实施例的B区标准地质分层与解释层位时差直方图。
图7示出了根据本发明的一个实施例的C区标准地质分层与解释层位时差直方图。
图8示出了根据本发明的一个实施例的D区标准地质分层与解释层位时差直方图。
图9示出了根据本发明的一个实施例的B分区伪时深关系井震标定的示意图。
图10示出了根据本发明的一个实施例的B分区校正后深关系井震标定的示意图。
图11示出了根据本发明的一个实施例的一种基于厚度分区的时深标定装置的框图。
附图标记说明:
201、追踪解释模块;202、厚度分区模块;203、标准井计算模块;204、标定模块。
具体实施方式
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
本发明提供一种基于厚度分区的时深标定方法,包括:
根据目标区域中标准井的井震标定结果进行自动追踪解释,获得地震层位顶、底的解释结果;
根据地震层位顶、底的解释结果计算时间厚度,进而根据时间厚度在目标区域进行区域划分,获得多个厚度分区;
确定每一个厚度分区的标准井,进而获得该标准井的时深关系;
在每一个厚度分区内,根据标准井的时深关系与非标准井的深度,计算非标准井的时间,实现基于该厚度分区的时深标定。
在一个示例中,地震层位顶的解释结果为:
Htop(x,y,t)=(xi,yj,ti) (1)
地震层位底的解释结果为:
Hbot(x,y,t)=(xi,yj,ti) (2)
其中,Htop(x,y,t)为地震层位顶的解释结果,Hbot(x,y,t)为地震层位底的解释结果,xi,yj,ti分别代表层位点在不同位置横坐标、纵坐标以及时间值。
在一个示例中,通过公式(3)计算时间厚度:
HD(xi,yi,Δti)=Hbot(xi,yi,ti)-Htop(xi,yi,ti) (3)
其中,HD(xi,yi,Δti)为时间厚度。
在一个示例中,在每一个厚度分区内,根据标准井的时深关系与非标准井的深度,计算非标准井的时间:
在每一个厚度分区内,根据标准井的时深关系与非标准井的深度,确定非标准井的深度对应的标准井的时间;
以标准井的时间作为非标准井的深度对应的初始时间,获得非标准井的伪时间函数;
根据地震层位的解释结果与非标准井的深度,获得非标准井的实际时间;
计算非标准井的初始时间与实际时间的时差;
根据非标准井的深度与初始时间以及时差,计算所有非标准井的时间,实现基于厚度分区的时深标定。
在一个示例中,非标准井的伪时间函数为:
其中,k代表分区序号,T为非标准井的初始时间,d为非标准井的深度,xi、yi为非标准井的纵横坐标。
在一个示例中,根据地震层位的解释结果与非标准井的深度,获得非标准井的实际时间包括:
根据地震层位的解释结果获得地震层位的解释结果时间域的空间函数;
将非标准井的深度代入地震层位的解释结果时间域的空间函数,获得非标准井的实际时间。
在一个示例中,通过公式(5)计算所有非标准井的时间:
其中,为非标准井的时间。
具体地,在三维地震覆盖区选取无扩径影响的密度与声波曲线的标准井,或对测井数据进行环境校正。
利用密度与声波曲线计算反射系数,确定子波主频类型,计算反射系数,进行VTI垂向各项异性标定,获取高精度的时深关系。
根据某一口标准井的井震标定结果,选取储层附近上下的强反射层或不整合面的地震反射同相轴,并对其自动追踪解释,高效率地获取反射界面的层位解释结果。设地震层位顶、底的解释结果时间域的空间函数分别为公式(1)、公式(2),进而利用公式(3)求差后获取时间厚度。
利用公式(6)、公式(7)将沉积厚度减掉最小值,然后除以常量后取整,对分区编号k=[Δti-MinΔti/50],如50ms的级差进行分区,200-250、250-300、350-400、400-450等。
其中,xi,yj,ti分别代表层位点在不同位置横坐标、纵坐标以及时间值,i,j,k为自然数。
根据厚度分区的结果,在每个分区内选取一口标准井,分别做合成记录标定。每个分区获取的时深结果如下:
将分区序号对应的标准井时深,分别将时深拷贝至各分区内无时深关系的井,获得伪时间函数为公式(4)。
在同一坐标点下,分别获取地质分层时间值与地震层位时间的差值:
即当xi=xo且yi=yq时:
则每口井的时差为:
SpickiΔti=WellPicki(Ti,di,xi,yi)-H(xo,yp,tq)
其中,i=(1,2,3....m),m为全工区井的总数量,o,p,q为自然数。
将不同井的时间差异量SpickiΔti,各自以常量的方式应用公式(4)的时间序列中,同时保持各自测井深度不变。亦即,仅更新所有井各自的时间序列。则批量计算表达式可表示为公式(5)。
如上可完成同一厚度分区的标定,针对每一个厚度分区均完成标定,实现对目标区域的基于厚度分区的时深标定。
本发明还提供一种基于厚度分区的时深标定装置,包括:
追踪解释模块,根据目标区域中标准井的井震标定结果进行自动追踪解释,获得地震层位顶、底的解释结果;
厚度分区模块,根据地震层位顶、底的解释结果计算时间厚度,进而根据时间厚度在目标区域进行区域划分,获得多个厚度分区;
标准井计算模块,确定每一个厚度分区的标准井,进而获得该标准井的时深关系;
标定模块,在每一个厚度分区内,根据标准井的时深关系与非标准井的深度,计算非标准井的时间,实现基于该厚度分区的时深标定。
在一个示例中,地震层位顶的解释结果为:
Htop(x,y,t)=(xi,yj,ti) (1)
地震层位底的解释结果为:
Hbot(x,y,t)=(xi,yj,ti) (2)
其中,Htop(x,y,t)为地震层位顶的解释结果,Hbot(x,y,t)为地震层位底的解释结果,xi,yj,ti分别代表层位点在不同位置横坐标、纵坐标以及时间值。
在一个示例中,通过公式(3)计算时间厚度:
HD(xi,yi,Δti)=Hbot(xi,yi,ti)-Htop(xi,yi,ti) (3)
其中,HD(xi,yi,Δti)为时间厚度。
在一个示例中,在每一个厚度分区内,根据标准井的时深关系与非标准井的深度,计算非标准井的时间:
在每一个厚度分区内,根据标准井的时深关系与非标准井的深度,确定非标准井的深度对应的标准井的时间;
以标准井的时间作为非标准井的深度对应的初始时间,获得非标准井的伪时间函数;
根据地震层位的解释结果与非标准井的深度,获得非标准井的实际时间;
计算非标准井的初始时间与实际时间的时差;
根据非标准井的深度与初始时间以及时差,计算所有非标准井的时间,实现基于厚度分区的时深标定。
在一个示例中,非标准井的伪时间函数为:
其中,k代表分区序号,T为非标准井的初始时间,d为非标准井的深度,xi、yi为非标准井的纵横坐标。
在一个示例中,根据地震层位的解释结果与非标准井的深度,获得非标准井的实际时间包括:
根据地震层位的解释结果获得地震层位的解释结果时间域的空间函数;
将非标准井的深度代入地震层位的解释结果时间域的空间函数,获得非标准井的实际时间。
在一个示例中,通过公式(5)计算所有非标准井的时间:
其中,为非标准井的时间。
具体地,在三维地震覆盖区选取无扩径影响的密度与声波曲线的标准井,或对测井数据进行环境校正。
利用密度与声波曲线计算反射系数,确定子波主频类型,计算反射系数,进行VTI垂向各项异性标定,获取高精度的时深关系。
根据某一口标准井的井震标定结果,选取储层附近上下的强反射层或不整合面的地震反射同相轴,并对其自动追踪解释,高效率地获取反射界面的层位解释结果。设地震层位顶、底的解释结果时间域的空间函数分别为公式(1)、公式(2),进而利用公式(3)求差后获取时间厚度。
利用公式(6)、公式(7)将沉积厚度减掉最小值,然后除以常量后取整,对分区编号k=[Δti-MinΔti/50],如50ms的级差进行分区,200-250、250-300、350-400、400-450等。
其中,xi,yj,ti分别代表层位点在不同位置横坐标、纵坐标以及时间值,i,j,k为自然数。
根据厚度分区的结果,在每个分区内选取一口标准井,分别做合成记录标定。每个分区获取的时深结果如下:
将分区序号对应的标准井时深,分别将时深拷贝至各分区内无时深关系的井,获得伪时间函数为公式(4)。
在同一坐标点下,分别获取地质分层时间值与地震层位时间的差值:
即当xi=xo且yi=yq时:
则每口井的时差为:
SpickiΔti=WellPicki(Ti,di,xi,yi)-H(xo,yp,tq)
其中,i=(1,2,3....m),m为全工区井的总数量,o,p,q为自然数。
将不同井的时间差异量SpickiΔti,各自以常量的方式应用公式(4)的时间序列中,同时保持各自测井深度不变。亦即,仅更新所有井各自的时间序列。则批量计算表达式可表示为公式(5)。
如上可完成同一厚度分区的标定,针对每一个厚度分区均完成标定,实现对目标区域的基于厚度分区的时深标定。
本发明还提供一种电子设备,电子设备包括:存储器,存储有可执行指令;处理器,处理器运行存储器中的可执行指令,以实现上述的基于厚度分区的时深标定方法。
本发明还提供一种计算机可读存储介质,该计算机可读存储介质存储有计算机程序,该计算机程序被处理器执行时实现上述的基于厚度分区的时深标定方法。
为便于理解本发明实施例的方案及其效果,以下给出四个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
实施例1
图1示出了根据本发明的基于厚度分区的时深标定方法的步骤的流程图。
如图1所示,该基于厚度分区的时深标定方法包括:步骤101,根据目标区域中标准井的井震标定结果进行自动追踪解释,获得地震层位顶、底的解释结果;步骤102,根据地震层位顶、底的解释结果计算时间厚度,进而根据时间厚度在目标区域进行区域划分,获得多个厚度分区;步骤103,确定每一个厚度分区的标准井,进而获得该标准井的时深关系;步骤104,在每一个厚度分区内,根据标准井的时深关系与非标准井的深度,计算非标准井的时间,实现基于该厚度分区的时深标定。
确定目的层层段并解释顶底反射层位,据解释结果计算厚度分区,选取确定不同分区经过环境校正的含有密度与声波曲线的标准井后,根据地震资料选取相应主频的子波,对声波时差曲线进行时域内的积分,求取时间序列后,据密度与声波求取时域内的反射系数并于子波进行卷积,获取初始的合成记录后,再对其进行精细的VTI波组标定。其次,将不同分区标准井的时深拷贝至不同分区内无时深关系的井,获取不同分区的伪时深关系,同时根据伪时深关系获取不同分区内所有井的标准地质分层时间,即标准地质分层(如煤层等)的伪时间值;同时根据不同分区标准井的时深关系,获取标准分层所对应得反射同相轴时间值、纵横坐标。最后,以具有相同的纵横坐标为必要条件,计算不同分区所有井标准地质层的伪时间与地震标准层反射时间的时间差。最后将求取的时间差,应用至不同分区所有井伪时深的时间序列,批量完成目标层段内所有井的时深关系求取。
图2示出了根据本发明的一个实施例的不同分区标准井位置以及声波时差井分布的示意图。
图3示出了根据本发明的一个实施例的无声波时差井分布的示意图。
图4示出了根据本发明的一个实施例的标定结果对比的示意图。
选取如图2所示的不同分区标准井W393、W248、W277、W403,据其中的密度与声波求取的时间域反射系数如图3所示。与子波褶积的传统标定与VTI标定如图4中a与b所示,两者对比可明显看出图4中b的合成记录与井旁的地震道的相关关系更高。
图5示出了根据本发明的一个实施例的A区标准地质分层与解释层位时差直方图。
图6示出了根据本发明的一个实施例的B区标准地质分层与解释层位时差直方图。
图7示出了根据本发明的一个实施例的C区标准地质分层与解释层位时差直方图。
图8示出了根据本发明的一个实施例的D区标准地质分层与解释层位时差直方图。
表1为部分井标准地质分层的地震反射时间与其伪时深分层时间的时间差值表,在时域地质分层的时深关系中,时间误差在-30ms~30ms之间,其中图5-8为A区-D区的误差表的直方图显示,即采用统一的时深关系,不能确保标定的准确性。
表1
井名 | 分层与层位时间差 | 分层点对应的X坐标 | 分层点对应的X坐标 |
W378 | 2.48 | 595223.30 | 5877775.30 |
W379 | -10.94 | 605074.20 | 5870782.60 |
W380 | -11.93 | 604653.60 | 5870827.00 |
W381 | 8.64 | 601452.90 | 5873404.00 |
W382 | -13.95 | 604528.80 | 5870776.00 |
W383 | 14.62 | 599226.90 | 5875579.00 |
W384 | 4.22 | 598604.60 | 5873803.50 |
W385 | 7.41 | 599427.50 | 5873813.50 |
W386 | 2.78 | 597048.40 | 5874345.50 |
W387 | -16.62 | 605145.10 | 5869670.00 |
W388 | -15.31 | 603633.10 | 5870910.00 |
W389 | -24.25 | 599402.00 | 5872481.00 |
W390 | -22.81 | 598622.70 | 5873205.00 |
W391 | -10.71 | 604657.50 | 5869941.50 |
W392 | 10.69 | 597093.80 | 5875500.00 |
W393 | -16.00 | 591393.40 | 5879873.00 |
W394 | 5.94 | 594013.80 | 5877050.00 |
W395 | -19.37 | 598793.00 | 5871460.50 |
W396 | -22.77 | 598743.00 | 5871457.50 |
表2为W150伪时深校正前后的对比表,其中更新前后的时差为一常量15.6ms,同时保持TVD,SSTVD或MD的值不变,也就是仅改变时间序列。
表2
井名 | 初始合成记录时间 | 更新后合成记录时间 | 时间差 | TVDSS | TVD |
W150 | 569.00 | 584.60 | 15.60 | 404.36 | 542.64 |
W150 | 571.00 | 586.60 | 15.60 | 402.42 | 544.58 |
W150 | 573.00 | 588.60 | 15.60 | 398.60 | 548.40 |
W150 | 575.00 | 590.60 | 15.60 | 393.14 | 553.86 |
W150 | 577.00 | 592.60 | 15.60 | 387.53 | 559.47 |
W150 | 579.00 | 594.60 | 15.60 | 381.87 | 565.13 |
W150 | 581.00 | 596.60 | 15.60 | 376.41 | 570.59 |
W150 | 583.00 | 598.60 | 15.60 | 371.11 | 575.89 |
W150 | 585.00 | 600.60 | 15.60 | 365.64 | 581.36 |
W150 | 587.00 | 602.60 | 15.60 | 360.58 | 586.42 |
W150 | 589.00 | 604.60 | 15.60 | 356.94 | 590.06 |
W150 | 591.00 | 606.60 | 15.60 | 353.42 | 593.58 |
W150 | 593.00 | 608.60 | 15.60 | 350.12 | 596.88 |
W150 | 595.00 | 610.60 | 15.60 | 346.84 | 600.17 |
W150 | 597.00 | 612.60 | 15.60 | 343.52 | 603.48 |
W150 | 599.00 | 614.60 | 15.60 | 340.23 | 606.77 |
图9示出了根据本发明的一个实施例的B分区伪时深关系井震标定的示意图。
图9为应用伪时深校正前后的标定,直接的时深拷贝导致W133,W150~W154井井震之间存在很大标定差,导致对储层反射位置认识不清,无法精确地进行地震反射界面的解释。
图10示出了根据本发明的一个实施例的B分区校正后深关系井震标定的示意图。
图10为批量校正后的无声波时差井的时深标定结果,明显可见砂体的反射位置基本处于黄色虚线的解释层位上。但W154、153、152以及133分层位置依然与解释层位存在标定误差。
实施例2
图11示出了根据本发明的一个实施例的一种基于厚度分区的时深标定装置的框图。
如图11所示,该基于厚度分区的时深标定装置,包括:
追踪解释模块201,根据目标区域中标准井的井震标定结果进行自动追踪解释,获得地震层位顶、底的解释结果;
厚度分区模块202,根据地震层位顶、底的解释结果计算时间厚度,进而根据时间厚度在目标区域进行区域划分,获得多个厚度分区;
标准井计算模块203,确定每一个厚度分区的标准井,进而获得该标准井的时深关系;
标定模块204,在每一个厚度分区内,根据标准井的时深关系与非标准井的深度,计算非标准井的时间,实现基于该厚度分区的时深标定。
作为可选方案,地震层位顶的解释结果为:
Htop(x,y,t)=(xi,yj,ti) (1)
地震层位底的解释结果为:
Hbot(x,y,t)=(xi,yj,ti) (2)
其中,Htop(x,y,t)为地震层位顶的解释结果,Hbot(x,y,t)为地震层位底的解释结果,xi,yj,ti分别代表层位点在不同位置横坐标、纵坐标以及时间值。
作为可选方案,通过公式(3)计算时间厚度:
HD(xi,yi,Δti)=Hbot(xi,yi,ti)-Htop(xi,yi,ti) (3)
其中,HD(xi,yi,Δti)为时间厚度。
作为可选方案,在每一个厚度分区内,根据标准井的时深关系与非标准井的深度,计算非标准井的时间:
在每一个厚度分区内,根据标准井的时深关系与非标准井的深度,确定非标准井的深度对应的标准井的时间;
以标准井的时间作为非标准井的深度对应的初始时间,获得非标准井的伪时间函数;
根据地震层位的解释结果与非标准井的深度,获得非标准井的实际时间;
计算非标准井的初始时间与实际时间的时差;
根据非标准井的深度与初始时间以及时差,计算所有非标准井的时间,实现基于厚度分区的时深标定。
作为可选方案,非标准井的伪时间函数为:
/>
其中,k代表分区序号,T为非标准井的初始时间,d为非标准井的深度,xi、yi为非标准井的纵横坐标。
作为可选方案,根据地震层位的解释结果与非标准井的深度,获得非标准井的实际时间包括:
根据地震层位的解释结果获得地震层位的解释结果时间域的空间函数;
将非标准井的深度代入地震层位的解释结果时间域的空间函数,获得非标准井的实际时间。
作为可选方案,通过公式(5)计算所有非标准井的时间:
其中,为非标准井的时间。
实施例3
本公开提供一种电子设备包括,该电子设备包括:存储器,存储有可执行指令;处理器,处理器运行存储器中的可执行指令,以实现上述基于厚度分区的时深标定方法。
根据本公开实施例的电子设备包括存储器和处理器。
该存储器用于存储非暂时性计算机可读指令。具体地,存储器可以包括一个或多个计算机程序产品,该计算机程序产品可以包括各种形式的计算机可读存储介质,例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(RAM)和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(ROM)、硬盘、闪存等。
该处理器可以是中央处理单元(CPU)或者具有数据处理能力和/或指令执行能力的其它形式的处理单元,并且可以控制电子设备中的其它组件以执行期望的功能。在本公开的一个实施例中,该处理器用于运行该存储器中存储的该计算机可读指令。
本领域技术人员应能理解,为了解决如何获得良好用户体验效果的技术问题,本实施例中也可以包括诸如通信总线、接口等公知的结构,这些公知的结构也应包含在本公开的保护范围之内。
有关本实施例的详细说明可以参考前述各实施例中的相应说明,在此不再赘述。
实施例4
本公开实施例提供一种计算机可读存储介质,该计算机可读存储介质存储有计算机程序,该计算机程序被处理器执行时实现所述的基于厚度分区的时深标定方法。
根据本公开实施例的计算机可读存储介质,其上存储有非暂时性计算机可读指令。当该非暂时性计算机可读指令由处理器运行时,执行前述的本公开各实施例方法的全部或部分步骤。
上述计算机可读存储介质包括但不限于:光存储介质(例如:CD-ROM和DVD)、磁光存储介质(例如:MO)、磁存储介质(例如:磁带或移动硬盘)、具有内置的可重写非易失性存储器的媒体(例如:存储卡)和具有内置ROM的媒体(例如:ROM盒)。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (7)
1.一种基于厚度分区的时深标定方法,其特征在于,包括:
根据目标区域中标准井的井震标定结果进行自动追踪解释,获得地震层位顶、底的解释结果;
根据地震层位顶、底的解释结果计算时间厚度,进而根据所述时间厚度在目标区域进行区域划分,获得多个厚度分区;
确定每一个厚度分区的标准井,进而获得该标准井的时深关系;
在每一个厚度分区内,根据所述标准井的时深关系与非标准井的深度,计算所述非标准井的时间,实现基于该厚度分区的时深标定;
其中,在每一个厚度分区内,根据所述标准井的时深关系与所述非标准井的深度,计算所述非标准井的时间:
在每一个厚度分区内,根据所述标准井的时深关系与所述非标准井的深度,确定所述非标准井的深度对应的所述标准井的时间;
以所述标准井的时间作为所述非标准井的深度对应的初始时间,获得所述非标准井的伪时间函数;
根据地震层位的解释结果与所述非标准井的深度,获得所述非标准井的实际时间;
计算所述非标准井的初始时间与实际时间的时差;
根据所述非标准井的深度与初始时间以及所述时差,计算所有非标准井的时间,实现基于厚度分区的时深标定;
其中,所述非标准井的伪时间函数为:
其中,k代表分区序号,T为非标准井的初始时间,d为非标准井的深度,x、y为非标准井的横纵坐标;
其中,通过公式(5)计算所有非标准井的时间:
其中,为非标准井的时间,SpickiΔti为时间差异量。
2.根据权利要求1所述的基于厚度分区的时深标定方法,其中,所述地震层位顶的解释结果为:
Htop(x,y,t)=(xi,yj,ti) (1)
所述地震层位底的解释结果为:
Hbot(x,y,t)=(xi,yj,ti) (2)
其中,Htop(x,y,t)为地震层位顶的解释结果,Hbot(x,y,t)为地震层位底的解释结果,xi,yj,ti分别代表层位点在不同位置横坐标、纵坐标以及时间值。
3.根据权利要求2所述的基于厚度分区的时深标定方法,其中,通过公式(3)计算时间厚度:
HD(xi,yi,Δti)=Hbot(xi,yi,ti)-Htop(xi,yi,ti) (3)
其中,HD(xi,yi,Δti)为时间厚度。
4.根据权利要求1所述的基于厚度分区的时深标定方法,其中,根据地震层位的解释结果与所述非标准井的深度,获得所述非标准井的实际时间包括:
根据地震层位的解释结果获得地震层位的解释结果时间域的空间函数;
将所述非标准井的深度代入地震层位的解释结果时间域的空间函数,获得所述非标准井的实际时间。
5.一种基于厚度分区的时深标定装置,其特征在于,包括:
追踪解释模块,根据目标区域中标准井的井震标定结果进行自动追踪解释,获得地震层位顶、底的解释结果;
厚度分区模块,根据地震层位顶、底的解释结果计算时间厚度,进而根据所述时间厚度在目标区域进行区域划分,获得多个厚度分区;
标准井计算模块,确定每一个厚度分区的标准井,进而获得该标准井的时深关系;
标定模块,在每一个厚度分区内,根据所述标准井的时深关系与非标准井的深度,计算所述非标准井的时间,实现基于该厚度分区的时深标定;
其中,在每一个厚度分区内,根据所述标准井的时深关系与所述非标准井的深度,计算所述非标准井的时间:
在每一个厚度分区内,根据所述标准井的时深关系与所述非标准井的深度,确定所述非标准井的深度对应的所述标准井的时间;
以所述标准井的时间作为所述非标准井的深度对应的初始时间,获得所述非标准井的伪时间函数;
根据地震层位的解释结果与所述非标准井的深度,获得所述非标准井的实际时间;
计算所述非标准井的初始时间与实际时间的时差;
根据所述非标准井的深度与初始时间以及所述时差,计算所有非标准井的时间,实现基于厚度分区的时深标定;
其中,所述非标准井的伪时间函数为:
其中,k代表分区序号,T为非标准井的初始时间,d为非标准井的深度,x、y为非标准井的横纵坐标;
其中,通过公式(5)计算所有非标准井的时间:
其中,为非标准井的时间,SpickiΔti为时间差异量。
6.一种电子设备,其特征在于,所述电子设备包括:
存储器,存储有可执行指令;
处理器,所述处理器运行所述存储器中的所述可执行指令,以实现权利要求1-4中任一项所述的基于厚度分区的时深标定方法。
7.一种计算机可读存储介质,其特征在于,该计算机可读存储介质存储有计算机程序,该计算机程序被处理器执行时实现权利要求1-4中任一项所述的基于厚度分区的时深标定方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210325161.2A CN116931089B (zh) | 2022-03-30 | 2022-03-30 | 基于厚度分区的时深标定方法、装置、电子设备及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210325161.2A CN116931089B (zh) | 2022-03-30 | 2022-03-30 | 基于厚度分区的时深标定方法、装置、电子设备及介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116931089A CN116931089A (zh) | 2023-10-24 |
CN116931089B true CN116931089B (zh) | 2024-04-05 |
Family
ID=88375865
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210325161.2A Active CN116931089B (zh) | 2022-03-30 | 2022-03-30 | 基于厚度分区的时深标定方法、装置、电子设备及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116931089B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102520447A (zh) * | 2012-01-06 | 2012-06-27 | 中国石油化工股份有限公司胜利油田分公司西部新区研究中心 | 油亮点识别系统 |
CN106772570A (zh) * | 2015-11-24 | 2017-05-31 | 中国石油化工股份有限公司 | 火成岩开启缝地震预测方法 |
CN110246138A (zh) * | 2018-03-09 | 2019-09-17 | 中国石油化工股份有限公司 | 一种数字岩心图像孔喉分割的方法 |
CN114154283A (zh) * | 2020-09-07 | 2022-03-08 | 中国石油化工股份有限公司 | 储层大孔径孔隙度计算方法、装置、电子设备及介质 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8332152B2 (en) * | 2008-03-19 | 2012-12-11 | Baker Hughes Incorporated | Method and apparatus for eliminating drill effect in pulse induction measurements |
-
2022
- 2022-03-30 CN CN202210325161.2A patent/CN116931089B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102520447A (zh) * | 2012-01-06 | 2012-06-27 | 中国石油化工股份有限公司胜利油田分公司西部新区研究中心 | 油亮点识别系统 |
CN106772570A (zh) * | 2015-11-24 | 2017-05-31 | 中国石油化工股份有限公司 | 火成岩开启缝地震预测方法 |
CN110246138A (zh) * | 2018-03-09 | 2019-09-17 | 中国石油化工股份有限公司 | 一种数字岩心图像孔喉分割的方法 |
CN114154283A (zh) * | 2020-09-07 | 2022-03-08 | 中国石油化工股份有限公司 | 储层大孔径孔隙度计算方法、装置、电子设备及介质 |
Non-Patent Citations (1)
Title |
---|
二连盆地巴音都兰凹陷隐蔽油藏勘探;李正文,等;中国石油勘探;20020630;第7卷(第2期);14-30 * |
Also Published As
Publication number | Publication date |
---|---|
CN116931089A (zh) | 2023-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10802171B2 (en) | High resolution seismic data derived from pre-stack inversion and machine learning | |
CN108802812B (zh) | 一种井震融合的地层岩性反演方法 | |
EP3204799B1 (en) | Conditioning of object or event based reservoir models using local multiple-point statistics simulations | |
US7337069B2 (en) | Method of simulating the sedimentary deposition in a basin respecting the thicknesses of the sedimentary sequences | |
US20090119076A1 (en) | Method for Generating a 3D Earth Model | |
US7065004B2 (en) | Horizon-based residual depth migration velocity analysis | |
US20130289958A1 (en) | Method and System for Modeling Anomalous Density Zones in Geophysical Exploration | |
CN109541685B (zh) | 一种河道砂体识别方法 | |
CN110058316B (zh) | 一种基于电阻率等值性原理的电磁测深约束反演方法 | |
CN105223617A (zh) | 利用地震属性控制初始波阻抗的建模方法及系统 | |
CN102937720A (zh) | 井控提高地震资料分辨率的方法 | |
WO2011056421A2 (en) | Sifting models of a subsurface structure | |
CN105425299A (zh) | 确定地层裂缝分布的方法和装置 | |
CN112114374B (zh) | 复杂地质体的多密度界面反演方法 | |
CN111983683B (zh) | 一种少井条件下盐下湖相灰岩储层预测方法和系统 | |
EP4018226B1 (en) | System and method for seismic imaging of complex subsurface volumes | |
CN116931089B (zh) | 基于厚度分区的时深标定方法、装置、电子设备及介质 | |
CN116856845A (zh) | 一种致密砂岩气藏水平井地质导向方法 | |
CN116859466B (zh) | 井震时深批量标定方法、装置、电子设备及介质 | |
US20210255346A1 (en) | System and Method for Seismic Imaging with Amplitude Recovery | |
CN112180466B (zh) | 一种河道砂体识别方法及装置 | |
CN112394391B (zh) | 雷克子波分数阶导数的地震相定量表征方法 | |
Dong | Integration of time-lapse seismic data into automatic history matching | |
CN111198401B (zh) | 基于多地球物理参数相关性的测井曲线校正方法与装置 | |
CN114488290A (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 |