CN108106582A - 基于多期构造活动背景下走滑断裂的走滑位移计算方法 - Google Patents
基于多期构造活动背景下走滑断裂的走滑位移计算方法 Download PDFInfo
- Publication number
- CN108106582A CN108106582A CN201711340849.3A CN201711340849A CN108106582A CN 108106582 A CN108106582 A CN 108106582A CN 201711340849 A CN201711340849 A CN 201711340849A CN 108106582 A CN108106582 A CN 108106582A
- Authority
- CN
- China
- Prior art keywords
- fault
- fracture
- slip
- disc
- sliding
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/32—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring the deformation in a solid
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种基于多期构造活动背景下走滑断裂的走滑位移计算方法,其中,方法包括:确定主构造轴线;测量走滑断裂第i组断层左右两盘基底卷入逆冲构造上、下盘的变形量Lupi、Lupi′、Ldowni及Ldowni′;测量走滑断裂第i组断层左右两盘滑脱断层的变形量Mi及Mi′;计算走滑断裂第i组断层左右两盘基底构造上、下盘构造线与主滑移带的夹角Φsi、Φsi、Φxi及Φsi′,上覆走滑断裂第i组断层左右两盘上覆滑脱层构造线与主滑移带的夹角ΦMi以及ΦMi′,分别计算Lupi、Lupi′、Ldowni、Ldowni′、Mi、以及Mi′在主滑移带的分量;根据上述分量计算走滑断裂的走滑位移量。采用本方案,可实现多期构造活动背景下走滑断层走滑位移量的计算,为后续规模开发井钻探及储量规模评价提供参考。
Description
技术领域
本申请实施例涉及地质勘探技术领域,具体涉及一种基于多期构造活动背景下走滑断裂的走滑位移计算方法。
背景技术
近年来,随着塔里木盆地内沿着走滑断裂带钻探成果的不断提升,尤其盆地内顺北地区沿北东向走滑断裂钻探的顺北系列井位的不断突破,以及不同强度级别上钻井的产能差异,进一步证实了区域内走滑断裂对深层岩溶缝洞型储层形成改造及油气运聚具有重要控制作用,从而加大对走滑断裂的识别及强度划分显得更加重要,而盆地内走滑断层走滑位移量是衡量走滑断裂活动强度及影响范围的一条重要指标。
目前,走滑断层平面位移计算的方法主要为以下几种:计算方法一:通过测量盖层扭动构造的压缩或拉伸量,估算出基底走滑断层的平移距离;计算方法二:在走滑拉分盆地中盆地沉降(或抬升)速率与边界断层走滑速率之间的数值关系,间接求取走滑位移;计算方法三:将走滑断层两侧同一地质体(或同时代地层)的位错距离,作为走滑运动指标来计算走滑断层平面位移大小;计算方法四,当沿着断层两侧沉积单元具有显著差异时,说明该断层可能经历了走滑运动,则可将不同沉积单元重叠的宽度作为走滑断层最小的平面走滑位移;计算方法五:在滑脱层与撕裂断层发育区域,利用剖面上逆冲双重构造的几何学形态恢复,来计算撕裂断层平面上走滑距离(挤压的缩短距离),忽略构造变形横向扩展性,根据构造变形前后剖面面积守恒原理,推算出撕裂断层最终平面走滑位移大小受控于逆冲断片的数量、双重构造面积以及双重构造的高度;计算方法六:撕裂断层和逆冲推覆构造共同发育的区域,推覆体或大型逆掩断层的上盘岩席,在侵蚀作用下蚀穿断面,使下盘岩石局部出露地表,形成“构造窗”,“构造窗”边缘到逆冲断层最大距离即为平面上撕裂断层的最小走滑距离(平面推覆的最小距离)。
然而,发明人在实现本发明的过程中发现,现有技术中的上述方式至少存在下述缺陷:
1、基于测量盖层扭动构造的压缩或拉伸量,计算基底走滑断层的平移距离的方法,一般假设断面直立无旋转,未考虑不同旋转角度下的走滑位移计算方法,旋转情况下走滑断裂断面复杂,在地震剖面较难找准同一条断裂分量,计算选点与拉伸收缩量数值读取难度大,并且未考虑旋转同时对计算精度的影响。
2、走滑拉分盆地中根据盆地沉降(或抬升)速率与边界断层走滑速率之间的数值关系,间接求取走滑位移的方法,主要适用于计算走滑拉分盆地中主动边界断层,并不不适用于盆地内部规模较小的对储层及油气运聚具有重要控制作用的走滑断裂位移量计算。
3、利用“构造窗”边缘到逆冲断层边界距离推算走滑位移的方法一般仅适用平面上露头发育区域,且推算出仅为最小的平面走滑位移。此外,根据剖面上逆冲双重构造的面积守恒原理,推算撕裂断层平面上走滑距离(挤压的缩短距离)的方法,不仅忽略了构造变形横向扩展性,同时也未考虑盖层滑脱量、剥蚀情况及对盘逆冲分量的影响。
4、利用特殊地质体(地震异常反射体、被走滑切割断层等)在走滑断裂两盘位移的变化来计算走滑量,但走滑断裂内标志性特殊地质体发现概率低,常规地震剖面上很难寻找到能作为走滑断裂位移量识别标志的特殊地质体,并且即使存在特殊地质体,也仅能对某一盖层内部走滑量进行表征,无法对走滑断裂纵横向整体活动性进行表征,并且选取不同的地质体对比参考往往导致不同的结果,达不到精细刻画的要求。
发明内容
鉴于上述问题,本发明提供了一种克服上述问题或者至少部分地解决上述问题的基于多期构造活动背景下走滑断裂的走滑位移计算方法。
根据本发明的一个方面提供了一种基于多期构造活动背景下走滑断裂的走滑位移计算方法,包括:
确定主构造轴线;
垂直于主构造轴线分别测量走滑断裂第i组断层左盘基底卷入逆冲构造上盘的变形量Lupi、走滑断裂第i组断层右盘基底卷入逆冲构造上盘的变形量Lupi′、走滑断裂第i组断层左盘基底卷入逆冲构造下盘的变形量Ldowni、以及走滑断裂第i组断层右盘基底卷入逆冲构造下盘的变形量Ldowni′;
垂直于主构造轴线分别测量出走滑断裂第i组断层左盘滑脱断层的变形量Mi、以及走滑断裂第i组断层右盘滑脱断层的变形量Mi′;
计算走滑断裂第i组断层左盘基底构造上盘构造线与主滑移带的夹角Φsi,走滑断裂第i组断层右盘基底构造上盘构造线与主滑移带的夹角Φsi′,走滑断裂第i组断层左盘基底构造下盘构造线与主滑移带的夹角Φxi,走滑断裂第i组断层右盘滑脱层主滑移带下盘构造线与主滑移带的夹角Φsi′,上覆走滑断裂第i组断层左盘上覆滑脱层构造线与主滑移带的夹角ΦMi,以及上覆走滑断裂第i组断层右盘上覆滑脱层构造线与主滑移带的夹角ΦMi,
分别计算Lupi、Lupi′、Ldowni、Ldowni′、Mi、以及Mi′在所述主滑移带的分量以及
根据所述以及计算走滑断裂的走滑位移量。
根据本发明提供的基于多期构造活动背景下走滑断裂的走滑位移计算方法,首先确定主构造轴线;并测量走滑断裂第i组断层左右两盘基底卷入逆冲构造上、下盘的变形量Lupi、Lupi′、Ldowni及Ldowni′;以及测量走滑断裂第i组断层左右两盘滑脱断层的变形量Mi及Mi′;计算走滑断裂第i组断层左右两盘基底构造上、下盘构造线与主滑移带的夹角Φsi、Φsi、Φxi及Φsi′,上覆走滑断裂第i组断层左右两盘上覆滑脱层构造线与主滑移带的夹角ΦMi以及ΦMi′,分别计算Lupi、Lupi′、Ldowni、Ldowni′、Mi、以及Mi′在主滑移带的分量;根据上述分量计算走滑断裂的走滑位移量。采用本方案,可实现多期构造活动背景下走滑断层走滑位移量的计算,为后续开发井钻探及储量规模评价提供参考,并可进一步为走滑断层相关的缝洞型储层新区块的经济效益评价提供评价基础,从而提高钻井成功率,提升勘探开发经济效益。
上述说明仅是本申请实施例技术方案的概述,为了能够更清楚了解本申请实施例的技术手段,而可依照说明书的内容予以实施,并且为了让本申请实施例的上述和其它目的、特征和优点能够更明显易懂,以下特举本申请的具体实施方式。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1示出了根据本发明一个实施例提供的基于多期构造活动背景下走滑断裂的走滑位移计算方法的流程示意图;
图2示出了多期构造活动背景下走滑断裂构造的示意图;
图3示出了根据本发明另一个实施例提供的基于多期构造活动背景下走滑断裂的走滑位移计算方法的流程示意图;
图4示出了一种采样点选取方式。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例,然而应当理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本公开,并且能够将本公开的范围完整的传达给本领域的技术人员。
图1示出了根据本发明一个实施例提供的基于多期构造活动背景下走滑断裂的走滑位移量计算方法。如图1所示,该方法包括:
步骤S110,确定主构造轴线。
其中,多期构造活动背景下走滑断裂构造可如图2所示。图2所示的多期构造活动背景下走滑断裂构造中包括基底构造21、第一滑脱层22、第二滑脱层23、逆冲构造24以及走滑断层25。其中,逆冲构造24是指岩石沿可观察到的或可推断出的低角(<30°)位移面产生大于5km水平位移的断层。在此,本领域人员应该理解的是,本发明中的多期构造活动背景下走滑断裂的构造包括但不限于图2所示的走滑断裂构造,本领域技术人员可根据实际情况获得不同的多期构造活动背景下的走滑断裂结构。
具体地,在计算走滑位移量过程中,首先选取研究区域,对研究区域中的走滑断裂进行断层构造分析,如断层解释,断层分区等,在对走滑断裂进行断层构造分析后进一步确定走滑断裂的主构造轴线(构造枢纽线)。其中,确定主构造轴线的具体方法本领域技术人员可自行选取,本发明在此不做限定。例如,可根据标志层的构造等值线图来确定主构造轴线等。
步骤S120,垂直于主构造轴线分别测量走滑断裂第i组断层左盘基底卷入逆冲构造上盘的变形量Lupi、走滑断裂第i组断层右盘基底卷入逆冲构造上盘的变形量Lupi′、走滑断裂第i组断层左盘基底卷入逆冲构造下盘的变形量Ldowni以及走滑断裂第i组断层右盘基底卷入逆冲构造下盘的变形量Ldowni′。
其中,i为大于或等于1的整数,且i小于或等于n,n为走滑断裂中划分的断层组数。在具体的实施过程中,可根据实际的走滑断裂构造确定走滑断裂构造中包含的断层组数。
以i=1为例,在垂直于主构造轴线的方向分别测量走滑断裂第一组断层左盘基底卷入逆冲构造上盘的变形量Lup1、走滑断裂断层右盘基底卷入逆冲构造上盘的变形量Lup1′、走滑断裂断层左盘基底卷入逆冲构造下盘的变形量Ldown1以及走滑断裂断层右盘基底卷入逆冲构造下盘的变形量Ldown1′。
步骤S130,垂直于主构造轴线分别测量出走滑断裂第i组断层左盘滑脱断层的变形量Mi、以及走滑断裂第i组断层右盘滑脱断层的变形量Mi′。
具体地,沿着滑脱断层主应变方向提取滑脱标志层对应的构造线数据,根据提取出的滑脱标志层对应的构造线数据,确定走滑断裂第i组断层左盘滑脱断层的变形量Mi、以及走滑断裂第i组断层右盘滑脱断层变形量Mi′。
步骤S140,计算走滑断裂第i组断层左盘基底构造上盘构造线与主滑移带的夹角Φsi,走滑断裂第i组断层右盘基底构造上盘构造线与主滑移带的夹角Φsi′,走滑断裂第i组断层左盘基底构造下盘构造线与主滑移带的夹角Φxi,走滑断裂第i组断层右盘滑脱层主滑移带下盘构造线与主滑移带的夹角Φsi′,上覆走滑断裂第i组断层左盘上覆滑脱层构造线与主滑移带的夹角ΦMi,以及上覆走滑断裂第i组断层右盘上覆滑脱层构造线与主滑移带的夹角ΦMi。
其中,本发明对Φsi、Φsi′、Φxi、Φsi′、ΦMi、和/或ΦMi的具体测量方法及测量精度等不做限定,本领域技术人员可根据实际的测量环境选择相应的测量方法及测量精度。
步骤S150,分别计算Lupi、Lupi′、Ldowni、Ldowni′、Mi、以及Mi′在主滑移带的分量LupΦsi、LupΦsi′、LdownΦxi、LdownΦxi′、MΦMi及MΦMi′。
其中,在计算以及过程中,首先计算Φsi、Φsi′、Φxi、Φsi′、ΦMi、以及ΦMi的余弦值,根据Φsi、Φsi′、Φxi、Φsi′、ΦMi、及ΦMi的余弦值与其对应的Lupi、Lupi′、Ldowni、Ldowni′、Mi、以及Mi′的乘法运算即可获得Lupi、Lupi′、Ldowni、Ldowni′、Mi、以及Mi′在主滑移带的分量以及
步骤S160,根据LupΦsi、LupΦsi′、LdownΦxi、LdownΦxi′、MΦMi及MΦMi′计算走滑断裂的走滑位移量。
具体地,根据以下公式1-1计算走滑断裂的走滑位移量:
其中,公式1-1中w为走滑断裂的走滑位移量,i为走滑断裂第i组断层,i为大于或等于1的整数,n为走滑断层中的断层组数。
根据本实施例提供的基于多期构造活动背景下走滑断裂的走滑位移计算方法,首先确定主构造轴线;并测量走滑断裂第i组断层左右两盘基底卷入逆冲构造上、下盘的变形量Lupi、Lupi′、Ldowni及Ldowni′;以及测量走滑断裂第i组断层左右两盘滑脱断层的变形量Mi及Mi′;计算走滑断裂第i组断层左右两盘基底构造上、下盘构造线与主滑移带的夹角Φsi、Φsi、Φxi、以及Φsi′,上覆走滑断裂第i组断层左右两盘上覆滑脱层构造线与主滑移带的夹角ΦMi以及ΦMi′,分别计算Lupi、Lupi′、Ldowni、Ldowni′、Mi、以及Mi′在主滑移带的分量;根据上述分量计算走滑断裂的走滑位移量。采用本方案,可实现多期构造活动背景下走滑断层走滑位移量的计算,为后续开发井钻探及储量规模评价提供参考,并可进一步为走滑断层相关的缝洞型储层新区块的经济效益评价提供评价基础,从而提高钻井成功率,提升勘探开发经济效益。
图3示出了根据本发明另一个实施例提供的基于多期构造活动背景下走滑断裂的走滑位移量计算方法。如图3所示,该方法包括:
步骤S310,追踪走滑断裂结构中的标志层,根据追踪的标志层确定主构造轴线。
具体地,首先选取走滑断裂结构中的标志层。在选取标志层的过程中,可根据获得地震数据,确定能够连续追踪的稳定界面,如泥岩、和/或膏盐岩与其他岩性界面形成的产生滑脱构造的强波阻抗界面(在地震剖面上变现为连续强反射特性),将该能够连续追踪的稳定界面确定为标志层。
进一步地,对走滑断裂上盘主要标志层以及走滑断裂下盘主要标志层进行追踪,根据追踪的走滑断裂上盘主要标志层以及走滑断裂下盘主要标志层的层位以及走滑断裂的断层构造特征确定主构造轴线。例如,根据走滑断裂上盘主要标志层的构造等值线图或者三维立体确定出的构造枢纽线即为主构造轴线。
步骤S320,垂直于主构造轴线分别测量走滑断裂第i组断层左盘基底卷入逆冲构造上盘的变形量Lupi、走滑断裂第i组断层右盘基底卷入逆冲构造上盘的变形量Lupi′、走滑断裂第i组断层左盘基底卷入逆冲构造下盘的变形量Ldowni以及走滑断裂第i组断层右盘基底卷入逆冲构造下盘的变形量Ldowni′。
其中,i为大于或等于1的整数,且i小于或等于n,n为走滑断裂中划分的断层组数。在具体的实施过程中,可根据实际的走滑断裂构造确定走滑断裂构造中包含的断层组数。
具体地,根据公式3-1计算走滑断裂第i组断层左盘基底卷入逆冲构造上盘的变形量Lupi。
其中,(xk,yk,zk)为沿着走滑断裂第i组断层左盘基底构造上盘构造线的第k个取样点的三维坐标值,(xk+1,yk+1,zk+1)为沿着滑断裂第i组断层左盘基底构造上盘构造线的第k+1个取样点的三维坐标值,其中,k为大于等于1的自然数,n1为沿着第i组断层左盘基底构造上盘构造线取样点的个数。
具体地,如图4所示,Oup为上盘零点,Odown为下盘零点,Fup为上盘终点,Fdown为下盘终点,则A1为第k个取样点,A2为第k+1个取样点,即A1及A2取样点位于上盘零点及上盘终点之间。
根据公式3-2计算走滑断裂第i组断层右盘基底卷入逆冲构造上盘的变形量Lupi′:
其中,(xk′,yk′,zk′)为沿着走滑断裂第i组断层右盘基底构造上盘构造线的第k个取样点的三维坐标值,(xk+1′,yk+1′,zk+1′)为沿着滑断裂第i组断层右盘基底构造上盘构造线的第k+1个取样点的三维坐标值,其中,k为大于等于1的自然数,n2为沿着第i组断层右盘基底构造上盘构造线取样点的个数。
根据公式3-3计算走滑断裂第i组断层左盘基底卷入逆冲构造下盘的变形量Ldowni::
其中,(xk",yk",zk")为沿着走滑断裂第i组断层左盘基底构造下盘构造线的第k个取样点的三维坐标值,(xk+1",yk+1",zk+1")为沿着滑断裂第i组断层左盘基底构造下盘构造线的第k+1个取样点的三维坐标值,其中,k为大于等于1的自然数,n3为沿着第i组断层左盘基底构造下盘构造线取样点的个数。
根据公式3-4计算走滑断裂第i组断层右盘基底卷入逆冲构造下盘的变形量Ldowni′:
其中,(xk”’,yk”',zk”')为沿着走滑断裂第i组断层左盘基底构造下盘构造线的第k个取样点的三维坐标值,(xk+1”',yk+1”',zk+1”')为沿着滑断裂第i组断层右盘基底构造下盘构造线的第k+1个取样点的三维坐标值,其中,k为大于等于1的自然数,n4为沿着第i组断层右盘基底构造下盘构造线取样点的个数。
在一种可选的实施方式中,n1=n2=n3=n4。进一步可选的,在选取取样点的过程中,可选取横向采用率为25米,纵向采用率为25米的采样方式。以根据公式3-1计算走滑断裂第i组断层左盘基底卷入逆冲构造上盘的变形量Lupi为例,当采用横向采用率为25米,纵向采用率为25米的采样方式时,xk+1=xk+25,yk+1=yk+25。其中,本领域技术人员可根据实际的地震数据选取合适的采样方式,本发明对具体的采样方式不做限定。
步骤S330,垂直于主构造轴线分别测量出走滑断裂第i组断层左盘滑脱断层的变形量Mi、以及走滑断裂第i组断层右盘滑脱断层的变形量Mi′。
具体地,根据公式3-5计算走滑断裂第i组断层左盘上覆滑脱层构造线的变形量Mi:
其中,(xk””,yk””,zk””)为沿着走滑断裂第i组断层左盘上覆滑脱层构造线的第k个取样点的三维坐标值,(xk+1””,yk+1””,zk+1””)为沿着滑断裂第i组断层左盘上覆滑脱层构造线的第k+1个取样点的三维坐标值,其中,k为大于等于1的自然数,m1为沿着第i组断层左盘上覆滑脱层构造线取样点的个数。
根据公式3-6计算走滑断裂第i组断层右盘上覆滑脱层构造线的变形量Mi':
其中,(xk””',yk””',zk””')为沿着走滑断裂第i组断层右盘上覆滑脱层构造线的第k个取样点的三维坐标值,(xk+1””',yk+1””',zk+1””’)为沿着滑断裂第i组断层右盘上覆滑脱层构造线的第k+1个取样点的三维坐标值,其中,k为大于等于1的自然数,m2为沿着第i组断层右盘上覆滑脱层构造线取样点的个数。
在一种可选的实施方式中,n1=n2=n3=n4=m1=m2。即在计算Lupi、Lupi′、Ldowni、Ldowni′、Mi、以及Mi′过程中可选取相同个数的取样点。
步骤S340,对Lupi、Lupi′、Ldowni、Ldowni′、Mi、和/或Mi′进行剥蚀量修正。
具体地,由于在多期构造活动背景下走滑断裂中的有关标志层中会有冲出地表而遭受剥蚀的情形,从而本步骤中可采用趋势分析法、地层对比法等方法对Lupi、Lupi′、Ldowni、Ldowni′、Mi、和/或Mi′进行剥蚀量修正。
步骤S350,计算走滑断裂第i组断层左盘基底构造上盘构造线与主滑移带的夹角Φsi,走滑断裂第i组断层右盘基底构造上盘构造线与主滑移带的夹角Φsi′,走滑断裂第i组断层左盘基底构造下盘构造线与主滑移带的夹角Φxi,走滑断裂第i组断层右盘滑脱层主滑移带下盘构造线与主滑移带的夹角Φsi′,上覆走滑断裂第i组断层左盘上覆滑脱层构造线与主滑移带的夹角ΦMi,以及上覆走滑断裂第i组断层右盘上覆滑脱层构造线与主滑移带的夹角ΦMi。
其中,本发明对Φsi、Φsi′、Φxi、Φsi′、ΦMi、和/或ΦMi的具体测量方法及测量精度等不做限定,本领域技术人员可根据实际的测量环境选择相应的测量方法及测量精度。
步骤S360,分别计算Lupi、Lupi′、Ldowni、Ldowni′、Mi、以及Mi′在主滑移带的分量LupΦsi、LupΦsi′、LdownΦxi、LdownΦxi′、MΦMi及MΦMi′。
具体地,根据计算
根据计算
根据计算
根据计算
根据计算以及
根据计算
步骤S370,根据LupΦsi、LupΦsi′、LdownΦxi、LdownΦxi′、MΦMi及MΦMi′计算走滑断裂的走滑位移量。
具体地,可参照图1所示基于多期构造活动背景下走滑断裂的走滑位移量计算方法实施例中的公式1-1计算走滑断裂的走滑位移量的具体阐述,在此不做赘述。
根据本实施例提供的基于多期构造活动背景下走滑断裂的走滑位移计算方法,首先追踪走滑断裂结构中的标志层,根据追踪的标志层确定主构造轴线;并测量走滑断裂第i组断层左右两盘基底卷入逆冲构造上、下盘的变形量Lupi、Lupi′、Ldowni及Ldowni′;以及测量走滑断裂第i组断层左右两盘滑脱断层的变形量Mi及Mi′;并对对Lupi、Lupi′、Ldowni、Ldowni′、Mi、和/或Mi′进行剥蚀量修正,从而进一步提高走滑位移准确性的提高;计算走滑断裂第i组断层左右两盘基底构造上、下盘构造线与主滑移带的夹角Φsi、Φsi、Φxi及Φsi′,上覆走滑断裂第i组断层左右两盘上覆滑脱层构造线与主滑移带的夹角ΦMi以及ΦMi′,分别计算Lupi、Lupi′、Ldowni、Ldowni′、Mi、以及Mi′在主滑移带的分量;根据上述分量计算走滑断裂的走滑位移量。采用本方案,可实现多期构造活动背景下走滑断层走滑位移量的计算,为后续开发井钻探及储量规模评价提供参考,并可进一步为走滑断层相关的缝洞型储层新区块的经济效益评价提供评价基础,从而提高钻井成功率,提升勘探开发经济效益。
最后应说明的是:以上实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。
Claims (8)
1.一种基于多期构造活动背景下走滑断裂的走滑位移量计算方法,其包括:
确定主构造轴线;
垂直于主构造轴线分别测量走滑断裂第i组断层左盘基底卷入逆冲构造上盘的变形量Lupi、走滑断裂第i组断层右盘基底卷入逆冲构造上盘的变形量Lupi′、走滑断裂第i组断层左盘基底卷入逆冲构造下盘的变形量Ldowni、以及走滑断裂第i组断层右盘基底卷入逆冲构造下盘的变形量Ldowni′;
垂直于主构造轴线分别测量出走滑断裂第i组断层左盘滑脱断层的变形量Mi、以及走滑断裂第i组断层右盘滑脱断层的变形量Mi′;
计算走滑断裂第i组断层左盘基底构造上盘构造线与主滑移带的夹角Φsi,走滑断裂第i组断层右盘基底构造上盘构造线与主滑移带的夹角Φsi′,走滑断裂第i组断层左盘基底构造下盘构造线与主滑移带的夹角Φxi,走滑断裂第i组断层右盘滑脱层主滑移带下盘构造线与主滑移带的夹角Φsi′,上覆走滑断裂第i组断层左盘上覆滑脱层构造线与主滑移带的夹角ΦMi,以及上覆走滑断裂第i组断层右盘上覆滑脱层构造线与主滑移带的夹角ΦMi;
分别计算Lupi、Lupi′、Ldowni、Ldowni′、Mi、以及Mi′在所述主滑移带的分量LupΦsi、LupΦsi′、LdownΦxi、LdownΦxi′、MΦMi及MΦMi′;
根据LupΦsi、LupΦsi′、LdownΦxi、LdownΦxi′、MΦMi及MΦMi′计算走滑断裂的走滑位移量。
2.根据权利要求1所述的方法,其特征在于,所述根据LupΦsi、LupΦsi′、LdownΦxi、LdownΦxi′、MΦMi及MΦMi′计算走滑断裂的走滑位移量具体包括:
根据以下公式计算走滑断裂的走滑位移量:
<mrow>
<mi>w</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>L</mi>
<mrow>
<msub>
<mi>up&Phi;</mi>
<mrow>
<mi>s</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</msub>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>L</mi>
<mrow>
<msub>
<mi>down&Phi;</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</msub>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>M</mi>
<msub>
<mi>&Phi;</mi>
<mrow>
<mi>M</mi>
<mi>i</mi>
</mrow>
</msub>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msubsup>
<mi>L</mi>
<mrow>
<msub>
<mi>up&Phi;</mi>
<mrow>
<mi>s</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msubsup>
<mi>L</mi>
<mrow>
<msub>
<mi>down&Phi;</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msubsup>
<mi>M</mi>
<msub>
<mi>&Phi;</mi>
<mrow>
<mi>M</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&prime;</mo>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
其中,w为走滑断裂的走滑位移量,i为走滑断裂第i组断层,i为大于或等于1的整数,n为大于或等于1的整数,且n为走滑断裂断层组数。
3.根据权利要求1或2所述的方法,其特征在于,所述确定主构造轴线进一步包括:
对走滑断裂上盘主要标志层以及走滑断裂下盘主要标志层进行追踪;
根据追踪的所述走滑断裂上盘主要标志层以及走滑断裂下盘主要标志层的层位以及走滑断裂的断层构造特征确定主构造轴线。
4.根据权利要求1或2所述的方法,其特征在于,所述垂直于主构造轴线分别测量走滑断裂第i组断层左盘基底卷入逆冲构造的变形量Lupi、走滑断裂第i组断层右盘基底卷入逆冲构造的变形量Lupi′、走滑断裂第i组断层左盘基底卷入逆冲构造下盘的变形量Ldowni、以及走滑断裂第i组断层右盘基底卷入逆冲构造下盘的变形量Ldowni′进一步包括:
根据公式:计算走滑断裂第i组断层左盘基底卷入逆冲构造上盘的变形量Lupi,其中,(xk,yk,zk)为沿着走滑断裂第i组断层左盘基底构造上盘构造线的第k个取样点的三维坐标值,(xk+1,yk+1,zk+1)为沿着滑断裂第i组断层左盘基底构造上盘构造线的第k+1个取样点的三维坐标值,其中,k为大于等于1的自然数,n1为沿着第i组断层左盘基底构造上盘构造线取样点的个数;
根据公式:计算走滑断裂第i组断层右盘基底卷入逆冲构造上盘的变形量Lupi′,其中,(xk′,yk′,zk′)为沿着走滑断裂第i组断层右盘基底构造上盘构造线的第k个取样点的三维坐标值,(xk+1′,yk+1′,zk+1′)为沿着滑断裂第i组断层右盘基底构造上盘构造线的第k+1个取样点的三维坐标值,其中,k为大于等于1的自然数,n2为沿着第i组断层右盘基底构造上盘构造线取样点的个数;
根据公式:计算走滑断裂第i组断层左盘基底卷入逆冲构造下盘的变形量Ldowni,其中,(xk",yk",zk")为沿着走滑断裂第i组断层左盘基底构造下盘构造线的第k个取样点的三维坐标值,(xk+1",yk+1",zk+1")为沿着滑断裂第i组断层左盘基底构造下盘构造线的第k+1个取样点的三维坐标值,其中,k为大于等于1的自然数,n3为沿着第i组断层左盘基底构造下盘构造线取样点的个数;
根据公式:计算走滑断裂第i组断层右盘基底卷入逆冲构造下盘的变形量Ldowni′,其中,(xk"',yk"',zk"')为沿着走滑断裂第i组断层左盘基底构造下盘构造线的第k个取样点的三维坐标值,(xk+1"',yk+1"',zk+1"')为沿着滑断裂第i组断层右盘基底构造下盘构造线的第k+1个取样点的三维坐标值,其中,k为大于等于1的自然数,n4为沿着第i组断层右盘基底构造下盘构造线取样点的个数。
5.根据权利要求4所述的方法,其特征在于,所述垂直于主构造轴线分别测量出走滑断裂第i组断层左盘滑脱断层的变形量Mi、以及走滑断裂第i组断层右盘滑脱断层的变形量Mi′具体包括:
根据公式:计算走滑断裂第i组断层左盘上覆滑脱层构造线的变形量Mi,其中,(xk"”,yk"”,zk"”)为沿着走滑断裂第i组断层左盘上覆滑脱层构造线的第k个取样点的三维坐标值,(xk+1"”,yk+1"”,zk+1"”)为沿着滑断裂第i组断层左盘上覆滑脱层构造线的第k+1个取样点的三维坐标值,其中,k为大于等于1的自然数,m1为沿着第i组断层左盘上覆滑脱层构造线取样点的个数;
根据公式:计算走滑断裂第i组断层右盘上覆滑脱层构造线的变形量Mi',其中,(xk"”',yk"”',zk"”')为沿着走滑断裂第i组断层右盘上覆滑脱层构造线的第k个取样点的三维坐标值,(xk+1"”',yk+1"”',zk+1'"”)为沿着滑断裂第i组断层右盘上覆滑脱层构造线的第k+1个取样点的三维坐标值,其中,k为大于等于1的自然数,m2为沿着第i组断层右盘上覆滑脱层构造线取样点的个数。
6.根据权利要求5所述的方法,其特征在于,n1=n2=n3=n4=m1=m2。
7.根据权利要求1或2所述的方法,其特征在于,所述分别计算Lupi、Lupi′、Ldowni、Ldowni′、Mi、以及Mi′在所述主滑移带的分量 以及具体包括:
根据计算
根据计算
根据计算
根据计算
根据计算以及
根据计算
8.根据权利要求1或2所述的方法,其特征在于,在所述分别计算Lupi、Lupi′、Ldowni、Ldowni′、Mi、以及Mi′在所述主滑移带的分量 以及之前,所述方法还包括:
对所述Lupi、Lupi′、Ldowni、Ldowni′、Mi、和/或Mi′进行剥蚀量修正。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711340849.3A CN108106582B (zh) | 2017-12-14 | 2017-12-14 | 基于多期构造活动背景下走滑断裂的走滑位移计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711340849.3A CN108106582B (zh) | 2017-12-14 | 2017-12-14 | 基于多期构造活动背景下走滑断裂的走滑位移计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108106582A true CN108106582A (zh) | 2018-06-01 |
CN108106582B CN108106582B (zh) | 2019-12-24 |
Family
ID=62216022
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711340849.3A Active CN108106582B (zh) | 2017-12-14 | 2017-12-14 | 基于多期构造活动背景下走滑断裂的走滑位移计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108106582B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109147531A (zh) * | 2018-07-26 | 2019-01-04 | 中国石油大学(华东) | 一种用于模拟复杂基底卷入型构造发育的变形装置 |
CN109211179A (zh) * | 2018-09-11 | 2019-01-15 | 福建工程学院 | 采矿滑裂面变形监测方法 |
CN109839661A (zh) * | 2019-01-18 | 2019-06-04 | 长江大学 | 确定走滑断层活动时间段和走滑位移的方法及设备 |
CN112859159A (zh) * | 2021-01-13 | 2021-05-28 | 中国石油大学(北京) | 一种雁列式正断层走滑带走滑位移量的计算方法 |
CN114647004A (zh) * | 2022-02-25 | 2022-06-21 | 中国地质大学(北京) | 一种用于确认地下走滑断裂不同时期的滑动方向的方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004197386A (ja) * | 2002-12-18 | 2004-07-15 | Railway Technical Res Inst | シールドトンネル通過時の地盤変位の予測方法 |
US8589080B2 (en) * | 2011-03-14 | 2013-11-19 | Southwest Research Institute | Geological stress inversion using fault displacement and slip tendency |
CN103955007A (zh) * | 2014-05-20 | 2014-07-30 | 中国石油化工股份有限公司胜利油田分公司西部新区研究院 | 复杂山前构造带的综合建模方法及建立的地质结构模型 |
CN104181595A (zh) * | 2014-08-24 | 2014-12-03 | 西南石油大学 | 用于复杂拉张构造体系断层伴生裂缝定量识别的新方法 |
CN105911589A (zh) * | 2016-06-27 | 2016-08-31 | 中国石油大学(华东) | 运用地震资料确定三维空间内断层滑动参数的定量化方法 |
CN106646605A (zh) * | 2016-10-16 | 2017-05-10 | 中国地质大学(北京) | 一种复杂断裂区岩层曲率优化计算方法 |
CN106842290A (zh) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | 一种扭张构造的平面表征方法 |
-
2017
- 2017-12-14 CN CN201711340849.3A patent/CN108106582B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004197386A (ja) * | 2002-12-18 | 2004-07-15 | Railway Technical Res Inst | シールドトンネル通過時の地盤変位の予測方法 |
US8589080B2 (en) * | 2011-03-14 | 2013-11-19 | Southwest Research Institute | Geological stress inversion using fault displacement and slip tendency |
CN103955007A (zh) * | 2014-05-20 | 2014-07-30 | 中国石油化工股份有限公司胜利油田分公司西部新区研究院 | 复杂山前构造带的综合建模方法及建立的地质结构模型 |
CN104181595A (zh) * | 2014-08-24 | 2014-12-03 | 西南石油大学 | 用于复杂拉张构造体系断层伴生裂缝定量识别的新方法 |
CN106842290A (zh) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | 一种扭张构造的平面表征方法 |
CN105911589A (zh) * | 2016-06-27 | 2016-08-31 | 中国石油大学(华东) | 运用地震资料确定三维空间内断层滑动参数的定量化方法 |
CN106646605A (zh) * | 2016-10-16 | 2017-05-10 | 中国地质大学(北京) | 一种复杂断裂区岩层曲率优化计算方法 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109147531A (zh) * | 2018-07-26 | 2019-01-04 | 中国石油大学(华东) | 一种用于模拟复杂基底卷入型构造发育的变形装置 |
CN109211179A (zh) * | 2018-09-11 | 2019-01-15 | 福建工程学院 | 采矿滑裂面变形监测方法 |
CN109839661A (zh) * | 2019-01-18 | 2019-06-04 | 长江大学 | 确定走滑断层活动时间段和走滑位移的方法及设备 |
CN112859159A (zh) * | 2021-01-13 | 2021-05-28 | 中国石油大学(北京) | 一种雁列式正断层走滑带走滑位移量的计算方法 |
CN114647004A (zh) * | 2022-02-25 | 2022-06-21 | 中国地质大学(北京) | 一种用于确认地下走滑断裂不同时期的滑动方向的方法 |
CN114647004B (zh) * | 2022-02-25 | 2024-03-15 | 中国地质大学(北京) | 一种用于确认地下走滑断裂不同时期的滑动方向的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108106582B (zh) | 2019-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108106582B (zh) | 基于多期构造活动背景下走滑断裂的走滑位移计算方法 | |
CN108301825B (zh) | 一种高温高压储层孔隙压力分析方法 | |
CN110424945B (zh) | 页岩油水平井部署方法及装置 | |
CN104881547B (zh) | 一种用于定向井井眼轨迹的误差分析方法 | |
CN108005646B (zh) | 基于随钻电磁波测井资料的地层各向异性电阻率提取方法 | |
CN110275204B (zh) | 地下走滑断裂水平滑移距定量分析方法 | |
CN110073246B (zh) | 与质量控制有关的改进的方法 | |
CN101647023A (zh) | 用于远程表征在钻孔附近的断层的方法和装置 | |
NO20140524A1 (no) | Inversjonsbasert arbeidsflyt for å behandle kjernetetthetsbilder i steilvinklede og horisontale brønner | |
Wang et al. | Calculating detachment depth and dip angle in sedimentary wedges using the area–depth graph | |
Zhou et al. | Improved wellbore quality using a novel real-time tortuosity index | |
AU2011374333B2 (en) | Method and system of correlating a measured log to a predicted log | |
CN110007344A (zh) | 一种断溶体储层连通性的地震识别方法及装置 | |
CN107103165A (zh) | 一种计算油田储层陡倾角剪切裂缝线密度的方法 | |
US11041381B2 (en) | Systems and methods for measuring rate of penetration | |
CN114063154B (zh) | 走滑断层位移量计算方法、装置、电子设备及介质 | |
Barrash et al. | Hydraulic conductivity distribution from multi-level slug tests and multivariate facies associations in a conglomeratic fluvial aquifer, Boise Hydrogeophysical Research Site | |
CN110618455B (zh) | 一种沉积盆地构造隆升的定量评价方法 | |
CN110764161B (zh) | 一种通过钻探资料判识碳酸盐岩断裂破碎带的综合方法 | |
CN116148929A (zh) | 井震结合上覆地层等效密度估算方法 | |
CN114296135A (zh) | 确定断层空间分段性、生长历史和真实长度的方法 | |
CN113050161A (zh) | 一种基于断层交面线的断裂走向断距图解计算方法 | |
CN111123362A (zh) | 裂缝性砂砾岩储层品质分类方法、装置及可读存储介质 | |
CN113514876B (zh) | 一种基于地层速度校正的设计井深度计算方法 | |
CN102944897B (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 |