CN112505773A - 一种断裂继承性量化计算和判定方法 - Google Patents
一种断裂继承性量化计算和判定方法 Download PDFInfo
- Publication number
- CN112505773A CN112505773A CN202011473545.6A CN202011473545A CN112505773A CN 112505773 A CN112505773 A CN 112505773A CN 202011473545 A CN202011473545 A CN 202011473545A CN 112505773 A CN112505773 A CN 112505773A
- Authority
- CN
- China
- Prior art keywords
- fracture
- target
- inheritance
- distance
- horizon
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 37
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000035945 sensitivity Effects 0.000 claims abstract description 5
- 206010017076 Fracture Diseases 0.000 claims description 111
- 208000010392 Bone Fractures Diseases 0.000 claims description 105
- 238000005070 sampling Methods 0.000 claims description 15
- 238000010586 diagram Methods 0.000 claims description 12
- 238000007476 Maximum Likelihood Methods 0.000 claims description 6
- 238000000547 structure data Methods 0.000 claims description 6
- 238000005553 drilling Methods 0.000 abstract description 2
- 230000003872 anastomosis Effects 0.000 abstract 1
- 238000011002 quantification Methods 0.000 description 3
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
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. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种断裂继承性量化计算和判定方法，其步骤为：步骤一，在目标层精细解释层位基础上进行不同属性分析，确定目标断裂的敏感属性；步骤二，采用断裂敏感属性和剖面结合三维空间综合解析进行目标断裂精细解释；步骤三，根据解释后的断裂在剖面上的特征，量化统计目标层断裂不同层位的垂向断距大小，并进行数据整理；步骤四，绘制目标层位断裂垂向断距大小趋势图，确定目标层位断裂的继承性，并利用相关系数计算公式来量化两者之间的继承率R。本发明能直观有效量化判定断裂继承性，能够大幅度提高储层预测吻合率、断溶体单井能力以及断溶体钻遇率，为勘探目标的优选，井位优化部署，提供强有力的技术支撑。
Description
技术领域
本发明涉及一种断裂继承性量化计算和判定方法，属于地震勘探领域。
背景技术
断控碳酸盐岩油气藏是西部主要油藏类型之一，西部顺北地区通源断裂、缝洞体发育，成藏条件有利，但该区地质目标体识别难度大，一方面是碳酸盐岩储层上奥陶统与中奥陶统波组抗差异大，对下伏目标层有很大的屏蔽作用；另一方面断裂成像不清晰，常规方法难以识别，导致断裂性质、规模以及期次演化认识不清，规模储层地震识别模式还不完善，制约了该区进一步的发展。正确认识断裂的性质、规模以及期次演化直接影响着断溶体油气藏的发现、储层预测以及后期开发，而断裂继承性是确定断裂性质、规模以及期次演化主要因素之一，以往断裂继承性通常是地质学家根据区域地质情况进行定性判断，方法单一，主观因素比重大，缺少有效的物探和量化计算方法，而且该方法仅限于经验丰富，对该区地质情况非常熟悉的地质专家，局限性大，不能在技术人员之间进行大范围内推广。因此，如何有效和直观判定断裂继承性尤为重要，为此，我们发明了一种断裂继承性量化计算和判定方法，解决了以上技术问题，能够直观有效直观量化判定断裂继承性，前景广阔，具有非常重要的意义。
发明内容
本发明的目的在于提供一种断裂继承性量化计算和判定方法，以克服上述现有技术存在的问题，实现有效直观量化判定断裂继承性。
为实现上述目的，本发明采取的技术方案是：一种断裂继承性量化计算和判定方法，其包括以下步骤：
步骤一，在目标层精细解释层位基础上进行不同属性分析，确定目标断裂的敏感属性；
步骤二，采用断裂敏感属性和剖面结合三维空间综合解析进行目标断裂精细解释；
步骤三，根据解释后的断裂在剖面上的特征，量化统计目标层断裂不同层位的垂向断距大小，并进行数据整理；
步骤四，绘制目标层位断裂垂向断距大小趋势图，确定目标层位断裂的继承性，并利用相关系数计算公式来量化两者之间的继承率R；其中，R越大说明上下目标层断裂相关系数越大，继承性越好，R越小说明上下目标层断裂相关系数越小，继承性越差。
进一步的，所述步骤一中，不同属性分析包括针对主干断裂和次级断裂进行AFE、相干、张量、不连续性、曲波多尺度相干、曲率、蚂蚁体、最大似然性的分析。
进一步的，所述步骤三中，量化统计目标层断裂不同层位的垂向断距大小，具体做法如下：
(1)在目标层精细解释层位基础上获取地质层面构造数据，所述地质层面构造数据包含多个取样点的位置信息以及取样点对应的时间值；
(2)以断裂最近稳定取样点作为计算点为标准，分别确定目标断裂两侧与断距计算点相对应的取样点为断距标准计算点；
(3)用断裂两侧标准计算点相对应的时间值进行相减得出断距大小。
进一步的，利用EXCEL对计算得出的断距大小取绝对值则得出最终不同层位的垂向断距大小，将不同目标层断裂垂向断距数据整理在一起，保证采样点数一致。
进一步的，所述步骤四中，利用EXCEL对步骤三整理好的目标层断裂垂向断距数据进行插入带平滑线的散点图，并根据断距大小确定图表标题、纵横向坐标轴值、图例格式以及图例位置，绘制目标层位断裂垂向断距大小趋势图，确定目标层位断裂的继承性，并利用相关系数计算公式来量化两者之间的继承率R。
进一步的，所述继承率R的计算公式如下：
其中，x和y分别为上下目标层断裂垂向断距离散值，xave和yave分别为上下目标层断裂垂向断距平均值，R为上下目标层断裂继承率，即上下断裂的相似系数。
本发明的有益效果是：
本发明克服了以前断裂继承性往往是地质学家根据区域地质情况进行定性判断，方法单一，主观因素比重大，缺少有效的物探和量化计算方法，而且该方法仅限于经验丰富，对该区地质情况非常熟悉的地质专家，局限性大，不能在技术人员之间进行大范围内推广等问题。
本发明能够直观有效直观量化判定断裂继承性，前景广阔，具有非常重要的意义。该方法的完善和推广能够大幅度提高储层预测吻合率、断溶体单井能力以及断溶体钻遇率，为勘探目标的优选，井位优化部署，提供强有力的技术支撑。
附图说明
图1是针对目标大尺度断裂多种平面属性图；其中(A)为AFE属性图、(B)为不连续性属性图、(C)为相干属性图、(D)为曲波多尺度相干属性图。
图2是针对目标中小尺度断裂多种平面属性图；其中(A)为最大似然性属性图、(B)为曲率属性图、(C)为蚂蚁体属性图、(D)为曲波中尺度属性图。
图3是典型断裂剖面、平面属性以及目标断裂三维可视化断层解释图.
图4是目标层断裂典型剖面。
图5是目标层断裂T7 4、T7 0层位的垂向断距大小数据整理统计图。
图6是T7 4、T7 0断距量化统计图。
图7是新疆YL地区T7 0、T7 4、T8 0、T9 0断距数据整理图。
图8是新疆YL地区T7 0、T7 4、T8 0、T9 0断距量化统计图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
本实施例公开了一种断裂继承性量化计算和判定方法，其包括以下步骤：
步骤1，在目标层精细解释层位基础上进行AFE、相干、张量、不连续性、曲波多尺度相干、曲率、蚂蚁体、最大似然性等不同属性分析，针对主干断裂和次级断裂进行了AFE、相干、张量、不连续性、曲波多尺度相干、曲率、蚂蚁体、最大似然性等不同属性分析，从测试效果分析可知，AFE、相干、不连续性、曲波多尺度相干能够清晰地识别和刻画出大尺度断裂特征，T74主干断裂平面特征清晰明显；曲率、蚂蚁体、最大似然性、曲波中尺度相干能够精细识别和刻画次级断裂特征，T74主干断裂内部分支断裂平面特征较明显，信息丰富。确定目标断裂的敏感属性，如图1和图2。
步骤2，在第一步敏感属性优选基础上，采用断裂敏感属性和剖面结合三维空间综合解析进行目标断裂精细解释，如图3和图4。
步骤3在精细解释目标断裂之后，根据解释后的断裂在剖面上的特征，量化统计目标层断裂不同层位的垂向断距大小，具体做法如下：
(1)在目标层精细解释层位基础上获取地质层面构造数据，所述地质层面构造数据包含多个取样点的位置信息以及取样点对应的时间值。
(2)以断裂最近稳定取样点作为计算点为标准，分别确定目标断裂两侧与所述断距计算点相对应的所述取样点为断距标准计算点。
(3)用断裂两侧标准计算点相对应的时间值进行相减得出断距大小。
利用EXCEL对计算得出的断距大小取绝对值则得出最终不同层位的垂向断距大小，将不同目标层断裂垂向断距数据整理在一起，保证采样点数一致，如图5。
步骤4，利用EXCEL对第三步整理好的目标层断裂垂向断距数据进行插入带平滑线的散点图，并根据断距大小确定图表标题、纵横向坐标轴值、图例格式以及图例位置等，绘制目标层位断裂垂向断距大小趋势图，确定目标层位断裂的继承性，并利用相关系数计算公式来量化两者之间的继承率R。
其中，x和y分别为上下目标层断裂垂向断距离散值，xave和yave分别为上下目标层断裂垂向断距平均值，R为上下目标层断裂继承率，即上下断裂的相似系数。R越大说明上下目标层断裂相关系数越大，继承性越好，R越小说明上下目标层断裂相关系数越小，继承性越差，如图6。
图8是新疆YL地区目标层断裂断距量化统计图表，采用本次发明的断裂继承性量化计算和判定方法对不同时段目标层进行断裂继承性量化计算，为了更科学地对比分析不同地质目标的断裂继承性，选取T70、T74、T80、T90四个目标层进行对比，分别对应的是奥陶系顶界面反射、中奥陶统一间房组顶界面反射、寒武系顶界反射以及震旦系顶界反射，目标层位时窗范围分别为3915-1075ms、4525-4645ms、5020-5150ms、5810-5935ms、从图7和图8断距大小统计图表中可以看出四者具有明显继承性，通过量化计算T90-T80、T80-T74、T74-T70三者继承性的相似度分别为0.96、0.76和0.78，由此可知，寒武系和震旦系继承性最好，奥陶系顶界面和中奥陶继承性次之，奥陶系及寒武系继承性相对最差，为该区后面断裂性质、规模以及期次演化的确定和分析提供了有力的指导。
以上显示和描述了本发明的基本原理、主要特征和优点。本领域的普通技术人员应该了解，上述实施例不以任何形式限制本发明的保护范围，凡采用等同替换等方式所获得的技术方案，均落于本发明的保护范围内。
本发明未涉及部分均与现有技术相同或可采用现有技术加以实现。
Claims (6)
1.一种断裂继承性量化计算和判定方法，其特征在于包括以下步骤：
步骤一，在目标层精细解释层位基础上进行不同属性分析，确定目标断裂的敏感属性；
步骤二，采用断裂敏感属性和剖面结合三维空间综合解析进行目标断裂精细解释；
步骤三，根据解释后的断裂在剖面上的特征，量化统计目标层断裂不同层位的垂向断距大小，并进行数据整理；
步骤四，绘制目标层位断裂垂向断距大小趋势图，确定目标层位断裂的继承性，并利用相关系数计算公式来量化两者之间的继承率R；其中，R越大说明上下目标层断裂相关系数越大，继承性越好，R越小说明上下目标层断裂相关系数越小，继承性越差。
2.根据权利要求1所述的一种断裂继承性量化计算和判定方法，其特征在于，所述步骤一中，不同属性分析包括针对主干断裂和次级断裂进行AFE、相干、张量、不连续性、曲波多尺度相干、曲率、蚂蚁体、最大似然性的分析。
3.根据权利要求1所述的一种断裂继承性量化计算和判定方法，其特征在于，所述步骤三中，量化统计目标层断裂不同层位的垂向断距大小，具体做法如下：
(1)在目标层精细解释层位基础上获取地质层面构造数据，所述地质层面构造数据包含多个取样点的位置信息以及取样点对应的时间值；
(2)以断裂最近稳定取样点作为计算点为标准，分别确定目标断裂两侧与断距计算点相对应的取样点为断距标准计算点；
(3)用断裂两侧标准计算点相对应的时间值进行相减得出断距大小。
4.根据权利要求3所述的一种断裂继承性量化计算和判定方法，其特征在于，利用EXCEL对计算得出的断距大小取绝对值则得出最终不同层位的垂向断距大小，将不同目标层断裂垂向断距数据整理在一起，保证采样点数一致。
5.根据权利要求1所述的一种断裂继承性量化计算和判定方法，其特征在于，所述步骤四中，利用EXCEL对步骤三整理好的目标层断裂垂向断距数据进行插入带平滑线的散点图，并根据断距大小确定图表标题、纵横向坐标轴值、图例格式以及图例位置，绘制目标层位断裂垂向断距大小趋势图，确定目标层位断裂的继承性，并利用相关系数计算公式来量化两者之间的继承率R。
6.根据权利要求1或5所述的一种断裂继承性量化计算和判定方法，其特征在于，所述继承率R的计算公式如下：
其中，x和y分别为上下目标层断裂垂向断距离散值，xave和yave分别为上下目标层断裂垂向断距平均值，R为上下目标层断裂继承率，即上下断裂的相似系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011473545.6A CN112505773B (zh) | 2020-12-15 | 2020-12-15 | 一种断裂继承性量化计算和判定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011473545.6A CN112505773B (zh) | 2020-12-15 | 2020-12-15 | 一种断裂继承性量化计算和判定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112505773A true CN112505773A (zh) | 2021-03-16 |
| CN112505773B CN112505773B (zh) | 2022-10-28 |
Family
ID=74973321
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011473545.6A Active CN112505773B (zh) | 2020-12-15 | 2020-12-15 | 一种断裂继承性量化计算和判定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112505773B (zh) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4698759A (en) * | 1985-03-28 | 1987-10-06 | The United States Of America As Represented By The United States Department Of Energy | Process for structural geologic analysis of topography and point data |
| CN104636980A (zh) * | 2013-11-13 | 2015-05-20 | 中国石油化工股份有限公司 | 针对河道砂油藏类型油气汇集条件的地球物理表征方法 |
| CN104914481A (zh) * | 2015-05-13 | 2015-09-16 | 中国石油大学(华东) | 一种火山岩区断层封闭性的综合评价方法 |
| CN105093300A (zh) * | 2015-07-27 | 2015-11-25 | 中国石油天然气股份有限公司 | 一种地质体边界识别方法及装置 |
| CN105317427A (zh) * | 2014-06-05 | 2016-02-10 | 中国石油化工股份有限公司 | 描述老断层顶面构造图的表示方法 |
| US20160209544A1 (en) * | 2015-01-15 | 2016-07-21 | Chevron U.S.A. Inc. | Quantitative assessment of plate tectonic models |
| CN108680952A (zh) * | 2018-04-03 | 2018-10-19 | 刘敬寿 | 一种走滑断层构造演化解析方法 |
| CN110133736A (zh) * | 2019-05-30 | 2019-08-16 | 安徽省勘查技术院(安徽省地质矿产勘查局能源勘查中心) | 一种覆盖区继承性断裂构造的金矿识别方法及系统 |
| CN110858001A (zh) * | 2018-08-22 | 2020-03-03 | 中国石油化工股份有限公司 | 一种针对深层碳酸盐岩走滑断裂带的解析方法 |
-
2020
- 2020-12-15 CN CN202011473545.6A patent/CN112505773B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4698759A (en) * | 1985-03-28 | 1987-10-06 | The United States Of America As Represented By The United States Department Of Energy | Process for structural geologic analysis of topography and point data |
| CN104636980A (zh) * | 2013-11-13 | 2015-05-20 | 中国石油化工股份有限公司 | 针对河道砂油藏类型油气汇集条件的地球物理表征方法 |
| CN105317427A (zh) * | 2014-06-05 | 2016-02-10 | 中国石油化工股份有限公司 | 描述老断层顶面构造图的表示方法 |
| US20160209544A1 (en) * | 2015-01-15 | 2016-07-21 | Chevron U.S.A. Inc. | Quantitative assessment of plate tectonic models |
| CN104914481A (zh) * | 2015-05-13 | 2015-09-16 | 中国石油大学(华东) | 一种火山岩区断层封闭性的综合评价方法 |
| CN105093300A (zh) * | 2015-07-27 | 2015-11-25 | 中国石油天然气股份有限公司 | 一种地质体边界识别方法及装置 |
| CN108680952A (zh) * | 2018-04-03 | 2018-10-19 | 刘敬寿 | 一种走滑断层构造演化解析方法 |
| CN110858001A (zh) * | 2018-08-22 | 2020-03-03 | 中国石油化工股份有限公司 | 一种针对深层碳酸盐岩走滑断裂带的解析方法 |
| CN110133736A (zh) * | 2019-05-30 | 2019-08-16 | 安徽省勘查技术院(安徽省地质矿产勘查局能源勘查中心) | 一种覆盖区继承性断裂构造的金矿识别方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 张伟健: "芳22-芳36区块扶余油层精细构造解释及储层预测", 《中国优秀硕士学位论文全文数据库》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112505773B (zh) | 2022-10-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104747185B (zh) | 非均质油藏储层综合分类评价方法 | |
| CN110674841B (zh) | 一种基于聚类算法的测井曲线识别方法 | |
| CN111458767B (zh) | 一种基于交汇图法识别岩性的方法及系统 | |
| CN103135131B (zh) | 一种针对裂缝性储层预测的解释装置 | |
| Zhang et al. | Architecture characterization of Ordovician fault-controlled paleokarst carbonate reservoirs in Tuoputai, Tahe oilfield, Tarim Basin, NW China | |
| CN103245971A (zh) | 井中断点引导的地震小断层解释方法及装置 | |
| CN103777243A (zh) | 砂泥岩薄互层储层厚度预测方法 | |
| CN110308444B (zh) | 道路层位智能识别及干扰源排除方法 | |
| CN105240006B (zh) | 一种适用于火山岩储层的油水层识别方法 | |
| CN107356965A (zh) | 基于加权叠加噪音压制策略的反射系数反演储层预测方法 | |
| CN111487691B (zh) | 一种致密砂岩储层岩性及砂体结构定量识别方法 | |
| CN108663714A (zh) | 一种沉积微相刻画方法 | |
| CN113254497A (zh) | 一种基于地球化学成分数据的综合分析及异常提取方法 | |
| CN113050168B (zh) | 基于阵列声波测井和声波远探测测井资料的裂缝有效性评价方法 | |
| CN110847887B (zh) | 一种细粒沉积陆相页岩裂缝识别评价方法 | |
| CN103777244B (zh) | 一种地震裂缝属性体的定量分析方法 | |
| CN112505773B (zh) | 一种断裂继承性量化计算和判定方法 | |
| CN116201535B (zh) | 油气藏目标井标志地层自动划分方法、装置及设备 | |
| CN110469327B (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 | |
| CN116027454A (zh) | 一种砂岩厚度等值线图的成图方法及智能绘制装置 | |
| AU2019373096B2 (en) | A method and system for sample classification | |
| CN109598049B (zh) | 钻孔岩石裂隙发育程度及区域岩石裂隙发育规律的方法 | |
| CN106291712B (zh) | 一种针对钻后数据的数据处理方法 | |
| CN118378101B (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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240904 Address after: 100027 Chaoyangmen North Street, Chaoyang District, Chaoyang District, Beijing Patentee after: SINOPEC Group Country or region after: China Patentee after: SINOPEC OILFIELD SERVICE Corp. Patentee after: SINOPEC PETROLEUM ENGINEERING GEOPHYSICS Co.,Ltd. Address before: 100029 Chaoyang District, Beijing Hui Xin Street six, Twelfth level. Patentee before: SINOPEC OILFIELD SERVICE Corp. Country or region before: China Patentee before: SINOPEC PETROLEUM ENGINEERING GEOPHYSICS Co.,Ltd. |
|
| TR01 | Transfer of patent right |