CN112147684B - 表征同沉积断层活动强度的方法和装置 - Google Patents
表征同沉积断层活动强度的方法和装置 Download PDFInfo
- Publication number
- CN112147684B CN112147684B CN201910572342.3A CN201910572342A CN112147684B CN 112147684 B CN112147684 B CN 112147684B CN 201910572342 A CN201910572342 A CN 201910572342A CN 112147684 B CN112147684 B CN 112147684B
- Authority
- CN
- China
- Prior art keywords
- fault
- stratum
- seismic
- determining
- well
- 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
- 230000000694 effects Effects 0.000 title claims abstract description 87
- 238000000034 method Methods 0.000 title claims abstract description 59
- 230000008859 change Effects 0.000 claims abstract description 12
- 238000012512 characterization method Methods 0.000 claims abstract description 8
- 238000013517 stratification Methods 0.000 claims description 21
- 238000002679 ablation Methods 0.000 claims description 20
- 238000005056 compaction Methods 0.000 claims description 8
- 238000011084 recovery Methods 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 6
- 238000011160 research Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 5
- 238000011835 investigation Methods 0.000 claims description 5
- 239000011435 rock Substances 0.000 claims description 4
- 238000000926 separation method Methods 0.000 claims description 2
- 238000011161 development Methods 0.000 description 17
- 230000008021 deposition Effects 0.000 description 15
- 230000007480 spreading Effects 0.000 description 6
- 230000007797 corrosion Effects 0.000 description 5
- 238000005260 corrosion Methods 0.000 description 5
- 230000001174 ascending effect Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005553 drilling Methods 0.000 description 4
- 230000002441 reversible effect Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000000737 periodic effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 230000003628 erosive effect Effects 0.000 description 2
- 239000013049 sediment Substances 0.000 description 2
- 238000005266 casting Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 230000003313 weakening 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. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- 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. analysis, for interpretation, for correction
- G01V1/30—Analysis
Abstract
本发明提供了一种表征同沉积断层活动强度的方法和装置,该方法包括:根据断层的三维地震数据体,确定垂直于断层走向且经过断层附近预定距离的井的时间域地震剖面;将垂直于断层走向且经过断层的井的时间域地震剖面,转换为深度域地震剖面;根据深度域地震剖面的海拔深度和地层格架中断层的铅直断距,确定各目标层的斜率;利用斜率,表征同沉积断层活动强度。本发明采用实际海拔深度和断层断距相结合的方法,利用线段斜率来反映断层活动强度,实现定量化精细表征同沉积断层多期次活动强度差异,并且表达地质含义更加充分,对纵向上不同时期断层活动强度的周期性变化体现更加直观。
Description
技术领域
本发明涉及油气田勘探开发阶段的分析领域,尤其涉及一种表征同沉积断层活动强度的方法和装置。
背景技术
当前,大量油气田勘探开发实践表明,在早期油气勘探领域,断层活动对圈闭的形成、可容空间的大小、宏观沉积体系和储层的展布、油气的运移、油气藏的保存和破坏具有重要的影响;在后期油气开发阶段,断层活动的重要意义体现在其对单个储集体的规模、叠置样式和展布方向的控制、井间连通关系的影响、断层伴生裂缝发育程度及其对开发效果和酸化压裂效果的影响。
因此,对断层的活动性开展研究,特别是同沉积断层(在地层沉积过程中形成发育的断层)活动强度大小以及不同时期活动强度的周期性变化的分析就显得格外重要。目前对在断层活动与否以及同沉积断层活动强度方面主要有以下研究方法。①生长指数法:生长指数虽然能够定量分析同沉积断层的活动强弱,但是其计算结果受上升盘地层剥蚀厚度影响较大,且无量纲,地质含义不够明确。②断层落差法:方法简单可行,不足之处在于没有体现时间概念。③断层活动速率法:断层活动速率法能够克服上升盘地层剥蚀或缺失的影响,且引入了时间概念,但前提是准确确定地层单元的沉积时间,而在油气田开发阶段较小地层单元条件下,准确确定不同地层单元的沉积时间难度很大。
此外,目前断层活动强度主要以数据表格或者直方图的形式来表现,单纯以数字表示地质含义不充分且不够直观。故,现有表征同沉积断层多期次活动强度的方法不够理想,针对其中存在的问题需建立更为合适的表征方法。
发明内容
为了充分且直观地表征断层活动强度的周期性变化,本发明提供一种表征同沉积断层活动强度的方法和装置。
第一方面,本发明提供一种表征同沉积断层活动强度的方法,所述方法包括:
根据断层的三维地震数据体,确定垂直于断层走向且经过断层附近预定距离的井的时间域地震剖面;
将所述垂直于断层走向且经过断层的井的时间域地震剖面,转换为深度域地震剖面;
根据所述深度域地震剖面的海拔深度和地层格架中断层的铅直断距,确定各目标层的斜率;
利用所述斜率,表征同沉积断层活动强度。
第二方面,本发明提供一种表征同沉积断层活动强度的装置,所述装置包括:
根据断层的三维地震数据体,确定垂直于断层走向且经过断层附近预定距离的井的时间域地震剖面;
将所述垂直于断层走向且经过断层的井的时间域地震剖面,转换为深度域地震剖面;
根据所述深度域地震剖面的海拔深度和地层格架中断层的铅直断距,确定各目标层的斜率;
活动强度表征模块,用于利用所述斜率,表征同沉积断层活动强度。
第三方面,本发明提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现第一方面提供的表征同沉积断层活动强度的方法的步骤。
第四方面,本发明提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现第一方面提供的表征同沉积断层活动强度的方法的步骤。
本发明采用实际海拔深度和断层断距相结合的方法,利用线段斜率来反映断层活动强度,实现定量化精细表征同沉积断层多期次活动强度差异,并且表达地质含义更加充分,对纵向上不同时期断层活动强度的周期性变化体现更加直观。
为让本发明的上述和其他目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附图式,作详细说明如下。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的表征同沉积断层活动强度的方法流程示意图;
图2为本发明实施例提供的井震标定合成地震记录;
图3为本发明实施例提供的时深关系拟合曲线;
图4为本发明实施例提供的断层活动强度定量表征折线图;
图5为本发明实施例提供的征同沉积断层活动强度的装置的框图;
图6为本发明实施例提供的电子设备框图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
目前在断层活动以及同沉积断层活动强度方面主要有以下研究方法。①生长指数法:Thorsen在1963最早提出断层生长指数这一概念,并用于判断断层活动强度。他将生长指数定义为断层两侧同一地层单元下降盘地层厚度与上升盘地层厚度之比。生长指数大于1代表断层活动,等于1代表断层不活动。生长指数虽然能够定量分析同沉积断层的活动强弱,但是其计算结果受上升盘地层剥蚀厚度影响较大,且无量纲,地质含义不够明确。②断层落差法:断层落差指的是同沉积断层上盘、下盘相同层位顶面或底面的铅直高度。断层落差不仅不受上升盘地层缺失的影响,且地质含义明确,方法简单可行,不足之处在于没有体现时间概念。③断层活动速率法:所谓断层活动速率指的是某段时间内断层两盘特定地层单元的厚度差与时间跨度的比值,断层活动速率越大,断层活动强度越大。断层活动速率法能够克服上升盘地层剥蚀或缺失的影响,且引入了时间概念,但前提是准确确定地层单元的沉积时间,而在油气田开发阶段较小地层单元条件下,准确确定不同地层单元的沉积时间难度很大。
此外,目前断层活动强度主要以数据表格或者直方图的形式来表现,如采用图示方法展示了断层的活动性,但主要从断层平面不同部位活动强度差异角度考虑,对纵向上不同时期断层活动强度的周期性变化考虑不够充分且表现不够直观。为解决上述问题,本发明实施例提供一种表征同沉积断层活动强度的方法,如图1所示,该方法包括:
步骤S101,根据断层的三维地震数据体,确定垂直于断层走向且经过断层附近预定距离的井的时间域地震剖面;
步骤S102,将垂直于断层走向且经过断层的井的时间域地震剖面,转换为深度域地震剖面;
步骤S103,根据深度域地震剖面的海拔深度和地层格架中断层的铅直断距,确定各目标层的斜率;
步骤S104,利用斜率,表征同沉积断层活动强度。
具体为,利用三维地震数据体对断层进行三维空间解释,明确断层的走向、倾向,根据断层的走向建立垂直于断层走向并且经过断层附近预定距离的井的时间域地震剖面。在这里预定距离可为断层附近的50m,100m等,可根据实际情况进行设置,不作具体数值限定。
将垂直于断层走向的时间域地震剖面和断层附近井的时间域地震剖面导出,并加载至2Dmove软件中,设定初速度和速度随深度的变化系数r,进行时深转换,将时间域地震剖面转换为深度域地震剖面,并定义该深度域地震剖面顶界面海拔深度值为0,向下海拔深度为正值且逐渐增大。
根据深度域地震剖面的海拔深度和研究区原始的地层格架中断层的铅直断距,确定各目标层的斜率,也即是在深度域地震剖面图一侧,以深度域地震剖面海拔深度0值为原点,以海拔深度建立横坐标轴x,将断层下降盘各目的层顶面海拔深度标注在深度坐标轴x之上,从原点开始,以铅直断距大小建立纵坐标轴y,将各目的层顶面铅直断距大小标注在坐标轴y之上。
将各目的层顶面海拔深度,如:a,以及铅直断距,如:b,在坐标系中投点,得到相应的(ai,bi)点,得到相应的(ai,bi)点,即{(a1,b1)、(a2,b2)、(a3,b3)…(an+1,bn+1)},其中n为目的层个数,按照纵坐标海拔深度值由小到大的顺序将各(ai,bi)点两两相连,得到一条由若干条线段组成的折线。计算各两点间组成线段的斜率Ki,Ki=(bi+1-bi/ai+1-ai),将各线段的斜率Ki标注在线段上,当Ki越小,即线段越平行x坐标,则代表目的层i沉积时期断层活动强度越小,反之则代表断层在目的层i沉积时期活动强度越大。
本发明实施例采用实际海拔深度和断层断距相结合的方法,利用线段斜率来反映断层活动强度,实现定量化精细表征同沉积断层多期次活动强度差异,并且表达地质含义更加充分,对纵向上不同时期断层活动强度的周期性变化体现更加直观。
基于上述实施例的内容,作为一种可选实施例:该方法还包括:
根据地震资料,建立典型井的地质分层与地震同相轴的对应关系;
根据地质分层与地震同相轴的对应关系,对三维地震区块的地层进行追踪对比,得到三维精细等时地层格架,确定断层位置。
具体为,根据地震资料,建立地质分层与地震同相轴的对应关系,在地震时间剖面上对地质分层开展侧向追踪对比,井震结合完成整个三维地震区块的地层追踪对比,从而建立三维精细等时地层格架,明确地层发育、剥蚀特征,从而确定断层位置。在这里,典型井是指取心长度大、取心较为完整、测井曲线系列较为完整且曲线特征明显的井,本发明实施例以典型井为计算标准,得到的地层发育更准确。
本发明实施例根据地质分层与地震同相轴的对应关系,对三维地震区块的地层进行追踪对比,使得到的地层发育更准确,能快速找到断层位置。
基于上述实施例的内容,作为一种可选实施例:根据断层的三维地震数据体,确定垂直于断层走向且经过断层附近预定距离的井的时间域地震剖面之前还包括:
根据岩心、录井和测井曲线资料,建立典型井地层划分标准;
根据典型井地层划分标准,对研究区所有井进行地层划分。
具体为,对研究区岩心、测井曲线资料较为丰富的井开展研究,建立典型井的地层划分标准,然后以典型井为中心,建立南北向和东西向两个方向的十字剖面,开展过十字剖面井的地层划分对比,之后过十字剖面上各井建立垂直于十字剖面的连井剖面,逐步完成全区所有井的地层划分对比。
基于上述实施例的内容,作为一种可选实施例:根据地震资料,建立典型井的地质分层与地震同相轴的对应关系包括:
根据地震资料,建立典型井的时间域地震剖面,进行合成地震记录;
根据合成地震记录和实际地震剖面,得到典型井的地质分层与地震同相轴的对应关系。
具体为,建立典型井的时间域地震剖面,利用典型井和剖面上所有具有声波时差曲线的井进行合成地震记录,如图2所示,确保合成记录与实际地震剖面对应关系良好,从而建立时间域地震剖面和深度域井剖面的时深关系,将典型井上的地质分层对应到相应的地震同相轴,如图3所示,可以看到地震反射时间与井上深度具有良好的正相关性,进而得到典型井的地质分层与地震同相轴的对应关系。
基于上述实施例的内容,作为一种可选实施例:根据地质分层与地震同相轴的对应关系,对三维地震区块的地层进行追踪对比,得到三维精细等时地层格架,确定断层位置包括:
利用井震结合完成对三维地震区块的地层进行追踪对比,得到三维精细等时地层格架,确定断层位置。
具体为,地质分层与地震同相轴的对应关系,在地震时间剖面上对地质分层开展侧向追踪对比,井震结合完成整个三维地震区块的地层追踪对比,从而建立三维精细等时地层格架,明确地层发育、剥蚀特征,从而确定断层位置。
基于上述实施例的内容,作为一种可选实施例:根据深度域地震剖面的海拔深度和研究区原始的地层格架中断层的铅直断距,确定各目标层的斜率包括:
根据邻层厚度比值法和参考层厚度变化率进行地层剥蚀厚度恢复,还原剥蚀前的地层格架;
对还原剥蚀前的地层格架中的地层压实量进行校正,得到原始的地层格架;
根据原始的地层格架,确定目标层的断层铅直断距;
根据深度域地震剖面中各目标层的海拔深度和目标层的断层铅直断距,确定各目标层的斜率。
具体为,根据地层格架发育特征,确定地层剥蚀情况,进行地层剥蚀厚度恢复,在油气田勘探阶段以及评价阶段,钻井数较少的情况下,主要利用地震资料进行地层剥蚀厚度的恢复。在油气田开发阶段,特别是开发中后期,此时钻井数量较多,可以在地震剥蚀厚度恢复的基础上,利用丰富的钻井资料选取邻层厚度比值法及参考层厚度变化率法进行地层剥蚀厚度恢复。
例如,在油气田属于开发后期,钻井数较多,平均井距约为200-400m,因此在地震剥蚀厚度恢复的基础上重点采用邻层厚度比值法及参考层厚度变化率法进行地层剥蚀厚度恢复,然后根据现今井上和地震地层厚度,以及计算出的地层剥蚀厚度,还原剥蚀前的地层格架展布特征。
在2Dmove软件中对地层压实量进行校正,得到压实前的地层厚度,即沉积刚结束时的沉积物厚度,从而得到沉积刚结束时原始的地层格架。例如,图4中,研究区原始沉积发育S1、S2、S3、S4、S5、S6、S7 1、S7 2、S7 3、S7 4共10个目的层单元,此时同沉积逆断层上盘地层厚度薄,下盘地层厚度厚。
研究区S1、S2、S3、S4、S5、S6、S7 1、S7 2、S7 3、S7 4每个目的层,分别统计其断层两盘同一地层顶面的铅直断距大小,然后根据深度域地震剖面的海拔深度和目标层的断层的铅直断距,确定各目标层的斜率。如图4所示,将S1、S2、S3、S4、S5、S6、S7 1、S7 2、S7 3、S7 4十个目标层的海拔深度和铅直断距表示在坐标系中,并两两连接,计算各两点间组成线段的斜率,线段越平行x坐标,斜率值越小,则代表目的层i沉积时期断层活动强度越小,反之则代表断层在目的层i沉积时期活动强度越大。
图中,S5顶到S6顶,即S5层,这两点所连线段斜率为2.8,代表在S5层沉积时期,同沉积逆断层活动强度很大。而S6顶到S7 1顶,即S6层,这两点所连线段斜率为0,代表在S6层沉积时期,断层不活动。总体来说,这10个目的层可以分为两个构造活动阶段,第一阶段为S1、S2、S3、S4、S5层沉积时期,同沉积逆断层总体活动强度较大,并随时间呈阶梯式减弱的特征;第二阶段为S6、S7 1、S7 2、S7 3、S7 4层沉积时期,同沉积逆断层总体活动强度较小,且早期(S7 2、S7 3、S7 4层)活动强度略大于晚期(S6、S7 1层)。
需要说明的是,在上述各实施例中表征同沉积断层活动强度的方法,均适用于同沉积正断层和同沉积逆断层。
本发明实施例通过还原剥蚀前的地层格架展布,以及进行地层压实量进行校正,得到原始的地层格架,进而准确确定目标层的断层铅直断距。
根据本发明的另一个方面,本发明实施例还提供一种表征同沉积断层活动强度的装置,参见图5,图5为本发明实施例提供的表征同沉积断层活动强度的装置的框图。该装置用于在前述各实施例中进行表征同沉积断层活动强度。因此,在前述各实施例中的表征同沉积断层活动强度的方法中的描述和定义,可以用于本发明实施例中各执行模块的理解。
如图所示,该装置包括:
时间域地震剖面确定模块501,用于根据断层的三维地震数据体,确定垂直于断层走向且经过断层附近预定距离的井的时间域地震剖面;
深度域地震剖面确定模块502,用于将垂直于断层走向且经过断层的井的时间域地震剖面,转换为深度域地震剖面;
斜率确定模块503,用于根据深度域地震剖面的海拔深度和地层格架中断层的铅直断距,确定各目标层的斜率;
活动强度表征模块504,用于利用斜率,表征同沉积断层活动强度。
本发明实施例采用实际海拔深度和断层断距相结合的方法,利用线段斜率来反映断层活动强度,实现定量化精细表征同沉积断层多期次活动强度差异,并且表达地质含义更加充分,对纵向上不同时期断层活动强度的周期性变化体现更加直观。
基于上述实施例的内容,作为一种可选实施例:该装置还包括:
对应关系建立模块,用于根据地震资料,建立典型井的地质分层与地震同相轴的对应关系;
断层位置确定模块,用于根据地质分层与地震同相轴的对应关系,对三维地震区块的地层进行追踪对比,得到三维精细等时地层格架,确定断层位置。
本发明实施例根据地质分层与地震同相轴的对应关系,对三维地震区块的地层进行追踪对比,使得到的地层发育更准确,能快速找到断层位置。
基于上述实施例的内容,作为一种可选实施例:时间域地震剖面确定模块之前还包括:
典型井地层划分模块,用于根据岩心、录井和测井曲线资料,建立典型井地层划分标准;
研究区所有井地层划分模块,用于根据典型井地层划分标准,对研究区所有井进行地层划分。
基于上述实施例的内容,作为一种可选实施例:对应关系建立模块包括:
合成地震记录确定单元,用于根据地震资料,建立典型井的时间域地震剖面,进行合成地震记录;
对应关系确定单元,用于根据合成地震记录和实际地震剖面,得到典型井的地质分层与地震同相轴的对应关系。
基于上述实施例的内容,作为一种可选实施例:断层位置确定模块包括:
断层位置确定单元,用于利用井震结合完成对三维地震区块的地层进行追踪对比,得到三维精细等时地层格架,确定断层位置。
基于上述实施例的内容,作为一种可选实施例:斜率确定模块包括:
剥蚀前的地层格架还原单元,用于根据邻层厚度比值法和参考层厚度变化率进行地层剥蚀厚度恢复,还原剥蚀前的地层格架展布;
原始的地层格架确定单元,用于对还原剥蚀前的地层格架中的地层压实量进行校正,得到原始的地层格架;
断层铅直断距确定单元,用于根据原始的地层格架,确定目标层的断层铅直断距;
各目标层的斜率确定单元,用于根据深度域地震剖面中各目标层的海拔深度和目标层的断层铅直断距,确定各目标层的斜率。
本发明实施例通过还原剥蚀前的地层格架展布,以及进行地层压实量进行校正,得到原始的地层格架,进而准确确定目标层的断层铅直断距。
图6为本发明实施例提供的电子设备框图,如图6所示,该设备包括:处理器601、存储器602和总线603;
其中,处理器601及存储602分别通过总线603完成相互间的通信;处理器601用于调用存储器602中的程序指令,以执行上述实施例所提供的表征同沉积断层活动强度的方法,例如包括:根据断层的三维地震数据体,确定垂直于断层走向且经过断层附近预定距离的井的时间域地震剖面;将垂直于断层走向且经过断层的井的时间域地震剖面,转换为深度域地震剖面;根据深度域地震剖面的海拔深度和地层格架中断层的铅直断距,确定各目标层的斜率;利用斜率,表征同沉积断层活动强度。
本发明实施例提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现表征同沉积断层活动强度的方法的步骤。例如包括:根据断层的三维地震数据体,确定垂直于断层走向且经过断层附近预定距离的井的时间域地震剖面;将垂直于断层走向且经过断层的井的时间域地震剖面,转换为深度域地震剖面;根据深度域地震剖面的海拔深度和地层格架中断层的铅直断距,确定各目标层的斜率;利用斜率,表征同沉积断层活动强度。
以上所描述的装置实施例仅仅是示意性的,其中作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分的方法。
最后,本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种表征同沉积断层活动强度的方法,其特征在于,所述方法包括:
根据断层的三维地震数据体,确定垂直于断层走向且经过断层附近预定距离的井的时间域地震剖面;
将所述垂直于断层走向且经过断层的井的时间域地震剖面,转换为深度域地震剖面;
根据所述深度域地震剖面的海拔深度和地层格架中断层的铅直断距,确定各目标层的斜率;
利用所述斜率,表征同沉积断层活动强度;
所述方法还包括:
根据地震资料,建立典型井的地质分层与地震同相轴的对应关系;典型井是指取心长度大、取心较为完整、测井曲线系列较为完整且曲线特征明显的井;
根据所述地质分层与地震同相轴的对应关系,对三维地震区块的地层进行追踪对比,得到三维精细等时地层格架,确定断层位置;
所述根据地震资料,建立典型井的地质分层与地震同相轴的对应关系包括:
根据地震资料,建立典型井的时间域地震剖面,进行合成地震记录;
根据所述合成地震记录和实际地震剖面,得到典型井的地质分层与地震同相轴的对应关系;
根据地震资料,建立典型井的时间域地震剖面,进行合成地震记录,包括:
建立典型井的时间域地震剖面,利用典型井和剖面上所有具有声波时差曲线的井进行合成地震记录。
2.根据权利要求1所述的方法,其特征在于,所述根据断层的三维地震数据体,确定垂直于断层走向且经过断层附近预定距离的井的时间域地震剖面之前还包括:
根据岩心、录井和测井曲线资料,建立典型井地层划分标准;
根据典型井地层划分标准,对研究区所有井进行地层划分。
3.根据权利要求1所述的方法,其特征在于,所述根据地质分层与地震同相轴的对应关系,对三维地震区块的地层进行追踪对比,得到三维精细等时地层格架,确定断层位置包括:
利用井震结合完成对三维地震区块的地层进行追踪对比,得到三维精细等时地层格架,确定断层位置。
4.根据权利要求1所述的方法,其特征在于,所述根据所述深度域地震剖面的海拔深度和研究区原始的地层格架中断层的铅直断距,确定各目标层的斜率包括:
根据邻层厚度比值法和参考层厚度变化率进行地层剥蚀厚度恢复,还原剥蚀前的地层格架;
对还原剥蚀前的地层格架中的地层压实量进行校正,得到原始的地层格架;
根据所述原始的地层格架,确定目标层的断层铅直断距;
根据深度域地震剖面中各目标层的海拔深度和目标层的断层铅直断距,确定各目标层的斜率。
5.一种表征同沉积断层活动强度的装置,其特征在于,所述装置包括:
时间域地震剖面确定模块,用于根据断层的三维地震数据体,确定垂直于断层走向且经过断层附近预定距离的井的时间域地震剖面;
深度域地震剖面确定模块,用于将所述垂直于断层走向且经过断层的井的时间域地震剖面,转换为深度域地震剖面;
斜率确定模块,用于根据所述深度域地震剖面的海拔深度和地层格架中断层的铅直断距,确定各目标层的斜率;
活动强度表征模块,用于利用所述斜率,表征同沉积断层活动强度;
所述装置还包括:
对应关系建立模块,用于根据地震资料,建立典型井的地质分层与地震同相轴的对应关系;典型井是指取心长度大、取心较为完整、测井曲线系列较为完整且曲线特征明显的井;
断层位置确定模块,用于根据所述地质分层与地震同相轴的对应关系,对三维地震区块的地层进行追踪对比,得到三维精细等时地层格架,确定断层位置;
所述对应关系建立模块包括:
合成地震记录确定单元,用于根据地震资料,建立典型井的时间域地震剖面,进行合成地震记录;
对应关系确定单元,用于根据所述合成地震记录和实际地震剖面,得到典型井的地质分层与地震同相轴的对应关系;
合成地震记录确定单元具体用于:
建立典型井的时间域地震剖面,利用典型井和剖面上所有具有声波时差曲线的井进行合成地震记录。
6.根据权利要求5所述的装置,其特征在于,所述时间域地震剖面确定模块之前还包括:
典型井地层划分模块,用于根据岩心、录井和测井曲线资料,建立典型井地层划分标准;
研究区所有井地层划分模块,用于根据典型井地层划分标准,对研究区所有井进行地层划分。
7.根据权利要求5所述的装置,其特征在于,所述断层位置确定模块包括:
断层位置确定单元,用于利用井震结合完成对三维地震区块的地层进行追踪对比,得到三维精细等时地层格架,确定断层位置。
8.根据权利要求5所述的装置,其特征在于,斜率确定模块包括:
剥蚀前的地层格架还原单元,用于根据邻层厚度比值法和参考层厚度变化率进行地层剥蚀厚度恢复,还原剥蚀前的地层格架;
原始的地层格架确定单元,用于对还原剥蚀前的地层格架中的地层压实量进行校正,得到原始的地层格架;
断层铅直断距确定单元,用于根据所述原始的地层格架,确定目标层的断层铅直断距;
各目标层的斜率确定单元,用于根据深度域地震剖面中各目标层的海拔深度和目标层的断层铅直断距,确定各目标层的斜率。
9.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1至4任一项所述表征同沉积断层活动强度的方法的步骤。
10.一种非暂态计算机可读存储介质,其上存储有计算机程序,其特征在于,该计算机程序被处理器执行时实现如权利要求1至4任一项所述表征同沉积断层活动强度的方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910572342.3A CN112147684B (zh) | 2019-06-28 | 2019-06-28 | 表征同沉积断层活动强度的方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910572342.3A CN112147684B (zh) | 2019-06-28 | 2019-06-28 | 表征同沉积断层活动强度的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112147684A CN112147684A (zh) | 2020-12-29 |
CN112147684B true CN112147684B (zh) | 2023-12-26 |
Family
ID=73869077
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910572342.3A Active CN112147684B (zh) | 2019-06-28 | 2019-06-28 | 表征同沉积断层活动强度的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112147684B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114764150A (zh) * | 2021-01-14 | 2022-07-19 | 中国石油天然气股份有限公司 | 一种断层活动强度和活动速率的计算方法及系统 |
CN112965111B (zh) * | 2021-02-05 | 2023-05-12 | 中国石油化工股份有限公司 | 一种断层解释方法 |
CN113376694B (zh) * | 2021-06-07 | 2024-01-23 | 中海石油(中国)有限公司深圳分公司 | 断陷盆地的活动性判断方法、装置、电子设备和存储介质 |
CN114415237B (zh) * | 2022-01-24 | 2022-08-23 | 东北石油大学 | 基于三维地震资料的砂岩型铀矿控矿断裂识别方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5537320A (en) * | 1994-10-31 | 1996-07-16 | Landmark Graphics Corporation | Method and apparatus for identifying fault curves in seismic data |
CN101046516A (zh) * | 2007-04-22 | 2007-10-03 | 罗仁泽 | 利用垂直地震剖面和微测井进行地震信号补偿方法 |
CN106199709A (zh) * | 2016-06-29 | 2016-12-07 | 中国石油化工股份有限公司 | 扭张断裂活动性的评价方法 |
CN106772602A (zh) * | 2016-12-22 | 2017-05-31 | 中国科学院地质与地球物理研究所 | 生长断层活动强度定量表征方法 |
CN108333627A (zh) * | 2018-01-16 | 2018-07-27 | 成都理工大学 | 火成岩地区断裂真伪的识别方法及装置 |
-
2019
- 2019-06-28 CN CN201910572342.3A patent/CN112147684B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5537320A (en) * | 1994-10-31 | 1996-07-16 | Landmark Graphics Corporation | Method and apparatus for identifying fault curves in seismic data |
CN101046516A (zh) * | 2007-04-22 | 2007-10-03 | 罗仁泽 | 利用垂直地震剖面和微测井进行地震信号补偿方法 |
CN106199709A (zh) * | 2016-06-29 | 2016-12-07 | 中国石油化工股份有限公司 | 扭张断裂活动性的评价方法 |
CN106772602A (zh) * | 2016-12-22 | 2017-05-31 | 中国科学院地质与地球物理研究所 | 生长断层活动强度定量表征方法 |
CN108333627A (zh) * | 2018-01-16 | 2018-07-27 | 成都理工大学 | 火成岩地区断裂真伪的识别方法及装置 |
Non-Patent Citations (1)
Title |
---|
断距系数在南阳凹陷张店油田中的应用;王冶;《石油地质与工程》;第31卷(第5期);第47-49页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112147684A (zh) | 2020-12-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112147684B (zh) | 表征同沉积断层活动强度的方法和装置 | |
CN110361781B (zh) | 一种基于地震反射振幅分析的沉积期古地貌恢复方法 | |
CN103472484A (zh) | 基于rs三维敏感地震属性分析的水平井轨迹优化方法 | |
CN104914465A (zh) | 火山岩裂缝定量预测方法及装置 | |
CN105319585B (zh) | 一种利用薄层干涉振幅恢复识别油气藏的方法 | |
CN104914470B (zh) | 碳酸盐岩缝洞储层储量校正方法 | |
CN107885893B (zh) | 描述二氧化碳驱储层非均质性的构建方法 | |
CN109543915B (zh) | 基于测井信息对全井段烃源岩总有机碳含量的识别方法 | |
CN105240006A (zh) | 一种适用于火山岩储层的油水层识别方法 | |
CN110632655B (zh) | 非均质型断层的封闭性分析方法及装置 | |
CN110632657B (zh) | 泥岩涂抹型断层封闭性分析方法及装置 | |
CN112489208A (zh) | 基于蚂蚁算法的裂缝片提取方法和三维地质模型构建方法 | |
CN111025391B (zh) | 一种断层活动性的定量评价方法 | |
CN110244358B (zh) | 一种构造破坏导致的油气逸散区的判识方法 | |
CN113742962B (zh) | 一种基于水平井的页岩储层属性三维建模方法 | |
CN110821484A (zh) | 一种陆相基质型页岩油地质储量的计算方法和存储介质 | |
CN106353813A (zh) | 基于阵列声波测井的流体性质识别方法 | |
CN110568149B (zh) | 沉积盆地烃源岩生排烃史精细快速定量模拟方法 | |
CN105528732B (zh) | 一种试气井产能预测方法 | |
CN112394416A (zh) | 非均质断裂控制储层预测方法及装置 | |
CN113971334A (zh) | 基于三维应力场数值模拟预测裂缝有效性的方法及装置 | |
CN105445797A (zh) | 一种获取圈闭含油气概率的方法 | |
God'spower et al. | Effects of erosional scours on reservoir properties of heterolithic, distal lower-shoreface sandstones | |
CN113495293B (zh) | 油藏流体预测方法及装置 | |
CN112835098B (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 |