CN110441821A - 一种基于变解释测网密度的断层快速解释方法 - Google Patents

一种基于变解释测网密度的断层快速解释方法 Download PDF

Info

Publication number
CN110441821A
CN110441821A CN201910829780.3A CN201910829780A CN110441821A CN 110441821 A CN110441821 A CN 110441821A CN 201910829780 A CN201910829780 A CN 201910829780A CN 110441821 A CN110441821 A CN 110441821A
Authority
CN
China
Prior art keywords
tomography
explained
fault
explanation
body attribute
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
Application number
CN201910829780.3A
Other languages
English (en)
Other versions
CN110441821B (zh
Inventor
范廷恩
杜昕
张显文
张晶玉
高云峰
董建华
王海峰
马良涛
范洪军
何明薇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Research Center of CNOOC China Ltd
CNOOC China Ltd
Original Assignee
Beijing Research Center of CNOOC China Ltd
CNOOC China Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Research Center of CNOOC China Ltd, CNOOC China Ltd filed Critical Beijing Research Center of CNOOC China Ltd
Priority to CN201910829780.3A priority Critical patent/CN110441821B/zh
Publication of CN110441821A publication Critical patent/CN110441821A/zh
Application granted granted Critical
Publication of CN110441821B publication Critical patent/CN110441821B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/301Analysis for determining seismic cross-sections or geostructures

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)
  • Apparatus For Radiation Diagnosis (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种基于变解释测网密度的断层快速解释方法,包括以下步骤:1)对待解释的地震数据进行解释性处理;2)利用方差体技术提取解释目的层附近的方差体属性切片;3)在提取的属性切片上解释一条断层,确定该条断层在平面上具有的拐点数目及位置;4)垂直于该条断层走向方向,分别提取平面上过断层线的首、尾两个端点,以及过该条断层各个拐点的地震剖面,在该地震剖面上完成该条断层的解释,形成变解释测网密度断层解释;5)在该地震剖面之间进行的断层内部插值性解释,完成该条断层的解释工作;6)重复步骤3)‑5)直到完成工区内所有断层的解释。本发明不仅具有更高的工作效率与解释精度,而且能够保证断层的空间闭合。

Description

一种基于变解释测网密度的断层快速解释方法
技术领域
本发明涉及一种快速断层解释方法,特别是关于一种针对复杂断块油气田的基于变解释测网密度的断层快速解释方法。
背景技术
由于复杂断块油气田断层多期次发育,断裂系统复杂,因此这类油田进入开发中后期后,伴随着剩余油分布预测及井位优化的需求,对地震构造解释尤其断层解释的精度亦有更为严苛的要求。近年来随着三维地震资料采集与处理技术的进步,以三维空间全自动解释为核心的全三维解释技术得到迅速发展。全三维解释技术根据地质认识设置空间种子点,通过自动追踪实现层位、断层的解释与砂体顶底的雕刻。在资料品质较好的情况下,该技术能取得令人满意的效果,然而该技术适用条件过于理想化,在实际生产应用中,针对各类地质条件,尤其是复杂断块型地质条件,构造形态复杂,地震资料品质差,该技术的效果不佳,缺乏实用性,存在断层不闭合,需反复修改的问题。
在此背景下,地震构造解释重新回到三维资料二维解释的轨道,即逐剖面进行解释然后在空间中寻求层位与断层的闭合,通过线—面—体思想,以三维工区抽稀后的主测线、联络线为框架,逐步加密测线解释来实现构造解释。但是,对于大型三维数据体,按照二维解释方法进行层位追踪、断层解释及断层平面组合,会存在层位、断层不闭合、断层平面组合不合理等问题,严重影响工作效率。因此,如何有效挖掘三维资料蕴含的内在信息,大幅度提高断层解释工作效率是目前地震断层解释工作亟需解决的问题。
发明内容
针对上述问题,本发明的目的是提供一种基于变解释测网密度的断层快速解释方法,该方法既能在复杂断块油田快速开展断层解释,又能保障断层空间闭合。
为实现上述目的,本发明采取以下技术方案:一种基于变解释测网密度的断层快速解释方法,包括以下步骤:
1)对待解释的地震数据进行提频、降噪、滤波等解释性处理,以提高地震数据品质,增强地震反射同相轴横向连续性,突出断层边界显示,得到能够充分突出断层信息的地震数据体;
2)利用方差体技术,提取解释目的层附近的方差体属性切片,作为确定断层拐点的数据基础;其中,解释目的层指在待解释的地震数据上确定的需要解释地层所在的时间范围,例如解释一个处于1200毫秒左右的地层,则目的层可以大致设置为1000-1500毫秒;
3)在提取的方差体属性切片上解释一条断层,基于该方差体属性切片上解释的断层形态,确定该条断层在方差体属性平面上具有的拐点数目及位置;
4)垂直于该条断层走向方向,分别提取步骤3)中方差体属性平面上过断层线的首、尾两个端点,以及过该条断层各个拐点的地震剖面,在该地震剖面上完成该条断层的解释,从而形成变解释测网密度断层解释;
5)利用断面网格化技术指导在该地震剖面之间进行的断层内部插值性解释,至此完成该条断层的解释工作;
6)重复步骤3)-5)直到完成工区内所有断层的解释,同时也完成该条断层的平面组合。
所述的断层快速解释方法,优选的,在所述步骤3)中,基于方差体属性切片上解释的断层形态,确定该条断层在方差体属性平面上具有的拐点数目及位置的方法如下:
①将提取的方差体属性切片上解释的某一条断层看做一条任意形状曲线,并定义该任意形状曲线的两端分别为该条断层的首端与尾端;
②从该任意形状曲线的首端开始,使用多个直线段近似拟合该任意形状曲线,假设能够完全拟合该任意形状曲线的形状时所用的直线段数目为n,则该条任意形状曲线存在的拐点数目为n-1,由此确定该条断层在方差体属性平面上具有的拐点数目;
③过相邻两条直线段的交点,与该任意形状曲线做连线,其中距离最短的线段与该任意形状曲线的交点可确定为该条断层拐点位置。
所述的断层快速解释方法,优选的,在所述步骤4)中,形成变解释测网密度断层解释的方法如下:
各个地震剖面之间的距离是由断层在方差体属性平面上的复杂程度所决定的,当该条断层在方差体属性平面上极限简单时,解释的断层线为直线,拐点数目为0,则只需要在该条断层的首尾两个端点对应地震剖面上开展解释;当该条断层在方差体属性平面上具有复杂形态时,则根据断层复杂程度自适应增加地震剖面断层解释数目,由此形成变解释测网密度断层解释。
所述的断层快速解释方法,优选的,在所述步骤5)中,利用断面网格化技术进行地震剖面之间的断层内部插值性解释的方法如下:
①预先定义插值数目为i,则在方差体属性平面上解释的断层线的相邻拐点之间,以及该条断层首、尾端点与邻近拐点之间再全面等分寻找i个点;
②沿着垂直于断层走向方向,提取上述i个点所对应得地震剖面,然后在该地震剖面上根据断面网格化技术指导开展断层解释,从而完成断层内部插值性解释。
本发明由于采取以上技术方案,其具有以下优点:本发明以变解释测网密度解释为核心,通过锁定断层拐点,以断层复杂程度为导向,快速完成断层剖面解释的同时也实现了断层的平面组合。本发明较传统二维断层解释技术具有更高的工作效率与解释精度,较传统三维断层解释技术更能够保证断层的空间闭合。
附图说明
图1是本发明的流程示意图;
图2是本发明实施例中某油田断层经过解释性处理后的地震数据体示意图;
图3是本发明实施例中某油田断层解释时从经过解释性处理后的地震数据体中提取的方差体属性切片示意图;
图4是本发明实施例中某油田断层解释时在上述方差体平面属性上解释的三条断层示意图;
图5是本发明实施例中三条断层的解释结果示意图;
图6是本发明实施例中对三条断层进行插值解释的结果示意图。
具体实施方式
以下将结合附图对本发明的较佳实施例进行详细说明,以便更清楚理解本发明的目的、特点和优点。应理解的是,附图所示的实施例并不是对本发明范围的限制,而只是为了说明本发明技术方案的实质精神。
下面以某进入开发中后期的复杂断块油田断层解释为例对本发明作进一步说明,如图1所示,本实施例提出的基于变解释测网密度的断层快速解释方法,包括以下步骤:
1)对该油田用于断层解释的地震数据进行提高分辨率、噪声压制、结构滤波等解释性处理,得到能够充分突出断层信息的地震数据体(如图2所示)。
2)利用方差体技术,提取解释目的层附近的属性切片,作为确定断层拐点的数据基础(如图3所示)。
3)在提取的方差体属性切片上解释三条断层(该三条断层在图4中从左至右显示,分别定义为断层1、断层2和断层3)。基于方差体属性切片上解释的这三条断层的形态,确定断层1具有5个拐点,断层2具有0个拐点(在平面上呈直线),断层3具有3个拐点。如图4所示,其中圆点代表断层的端点或拐点。
4)垂直于断层走向方向,分别提取过断层1、断层2和断层3的首、尾两个端点,以及过三条断层各个拐点的地震剖面,在上述地震剖面上完成这三条断层的解释,从而形成变解释测网密度断层解释。其中断层1解释2(首尾端点)+5(拐点数)=7个地震剖面,断层2解释2(首尾端点)+0(拐点数)=2个地震剖面,断层3解释2(首尾端点)+3(拐点数)=5个地震剖面,最终解释结果如图5所示。
5)利用断面网格化技术指导在上述地震剖面之间进行的断层内部插值性解释,相邻两个地震剖面间再选一个地震剖面进行断层解释,则断层1再解释2+5-1=6(间隔数)个地震剖面,断层2再解释2+0-1=1(间隔数)个地震剖面,断层3再解释2+3-1=4(间隔数)个地震剖面,最终解释结果如图6所示,其中圆点代表断层插值解释时提取的地震剖面位置,虚线代表断层插值解释的结果。
6)重复步骤3)-5)直到完成工区内所有断层的解释,同时也完成了断层的平面组合。
上述各实施例仅用于说明本发明,其中各部件的结构、连接方式和制作工艺等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。

Claims (5)

1.一种基于变解释测网密度的断层快速解释方法,其特征在于,包括以下步骤:
1)对待解释的地震数据进行解释性处理,以提高地震数据品质,增强地震反射同相轴横向连续性,突出断层边界显示,得到能够充分突出断层信息的地震数据体;
2)利用方差体技术,提取解释目的层附近的方差体属性切片,作为确定断层拐点的数据基础;
3)在提取的方差体属性切片上解释一条断层,基于该方差体属性切片上解释的断层形态,确定该条断层在方差体属性平面上具有的拐点数目及位置;
4)垂直于该条断层走向方向,分别提取步骤3)中方差体属性平面上过断层线的首、尾两个端点,以及过该条断层各个拐点的地震剖面,在该地震剖面上完成该条断层的解释,从而形成变解释测网密度断层解释;
5)利用断面网格化技术指导在该地震剖面之间进行的断层内部插值性解释,至此完成该条断层的解释工作;
6)重复步骤3)-5)直到完成工区内所有断层的解释,同时完成该条断层的平面组合。
2.根据权利要求1所述的断层快速解释方法,其特征在于,在所述步骤3)中,基于方差体属性切片上解释的断层形态,确定该条断层在方差体属性平面上具有的拐点数目及位置的方法如下:
①将提取的方差体属性切片上解释的某一条断层看做一条任意形状曲线,并定义该任意形状曲线的两端分别为该条断层的首端与尾端;
②从该任意形状曲线的首端开始,使用多个直线段近似拟合该任意形状曲线,假设能够完全拟合该任意形状曲线的形状时所用的直线段数目为n,则该条任意形状曲线存在的拐点数目为n-1,由此确定该条断层在方差体属性平面上具有的拐点数目;
③过相邻两条直线段的交点,与该任意形状曲线做连线,其中距离最短的线段与该任意形状曲线的交点可确定为该条断层拐点位置。
3.根据权利要求1所述的断层快速解释方法,其特征在于,在所述步骤4)中,形成变解释测网密度断层解释的方法如下:
各个地震剖面之间的距离是由断层在方差体属性平面上的复杂程度所决定的,当该条断层在方差体属性平面上极限简单时,解释的断层线为直线,拐点数目为0,则只需要在该条断层的首尾两个端点对应地震剖面上开展解释;当该条断层在方差体属性平面上具有复杂形态时,则根据断层复杂程度自适应增加地震剖面断层解释数目,由此形成变解释测网密度断层解释。
4.根据权利要求1所述的断层快速解释方法,其特征在于,在所述步骤5)中,利用断面网格化技术进行地震剖面之间的断层内部插值性解释的方法如下:
①预先定义插值数目为i,则在方差体属性平面上解释的断层线的相邻拐点之间,以及该条断层首、尾端点与邻近拐点之间再全面等分寻找i个点;
②沿着垂直于断层走向方向,提取上述i个点所对应得地震剖面,然后在该地震剖面上根据断面网格化技术指导开展断层解释,从而完成断层内部插值性解释。
5.根据权利要求1所述的断层快速解释方法,其特征在于,在所述步骤1)中,解释性处理包括对待解释的地震数据进行提频、降噪和滤波处理。
CN201910829780.3A 2019-09-03 2019-09-03 一种基于变解释测网密度的断层快速解释方法 Active CN110441821B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910829780.3A CN110441821B (zh) 2019-09-03 2019-09-03 一种基于变解释测网密度的断层快速解释方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910829780.3A CN110441821B (zh) 2019-09-03 2019-09-03 一种基于变解释测网密度的断层快速解释方法

Publications (2)

Publication Number Publication Date
CN110441821A true CN110441821A (zh) 2019-11-12
CN110441821B CN110441821B (zh) 2020-10-09

Family

ID=68438949

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910829780.3A Active CN110441821B (zh) 2019-09-03 2019-09-03 一种基于变解释测网密度的断层快速解释方法

Country Status (1)

Country Link
CN (1) CN110441821B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111856572A (zh) * 2020-07-06 2020-10-30 中国石油大学(北京) 一种断层破碎带宽度的确定方法及装置

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5229976A (en) * 1991-11-06 1993-07-20 Conoco Inc. Method for creating a numerical model of the physical properties within the earth
WO2002073240A1 (en) * 2001-03-13 2002-09-19 Conoco Phillips Company Method and process for prediction of subsurface fluid and rock pressures in the earth
EP1653411A1 (fr) * 2004-10-26 2006-05-03 Total S.A. Procédé et programme de reconstruction de plan de fracture
WO2008111037A2 (en) * 2007-03-12 2008-09-18 Geomage 2003 Ltd A method for identifying and analyzing faults/fractures using reflected and diffracted waves
CN103135136A (zh) * 2011-11-25 2013-06-05 中国石油化工股份有限公司 一种三维地震数据体的断层自动解释装置
CN103760597A (zh) * 2013-12-27 2014-04-30 淮南万泰电子股份有限公司 一种矿山断层自动识别方法
CN103941286A (zh) * 2014-04-16 2014-07-23 中国海洋石油总公司 一种基于平面导航的快速三维层位解释方法
WO2015070022A1 (en) * 2013-11-08 2015-05-14 Schlumberger Canada Limited Stratigraphic and structural interpretation of deviated and horizontal wellbores
CN106556869A (zh) * 2017-01-19 2017-04-05 中国石油大学(华东) 一种拉张环境下定量刻画断层走滑位移的方法
CN109975872A (zh) * 2019-04-26 2019-07-05 河南理工大学 一种利用Croel DRAW软件直接编制平衡剖面的方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5229976A (en) * 1991-11-06 1993-07-20 Conoco Inc. Method for creating a numerical model of the physical properties within the earth
WO2002073240A1 (en) * 2001-03-13 2002-09-19 Conoco Phillips Company Method and process for prediction of subsurface fluid and rock pressures in the earth
EP1653411A1 (fr) * 2004-10-26 2006-05-03 Total S.A. Procédé et programme de reconstruction de plan de fracture
WO2008111037A2 (en) * 2007-03-12 2008-09-18 Geomage 2003 Ltd A method for identifying and analyzing faults/fractures using reflected and diffracted waves
CN103135136A (zh) * 2011-11-25 2013-06-05 中国石油化工股份有限公司 一种三维地震数据体的断层自动解释装置
WO2015070022A1 (en) * 2013-11-08 2015-05-14 Schlumberger Canada Limited Stratigraphic and structural interpretation of deviated and horizontal wellbores
CN103760597A (zh) * 2013-12-27 2014-04-30 淮南万泰电子股份有限公司 一种矿山断层自动识别方法
CN103941286A (zh) * 2014-04-16 2014-07-23 中国海洋石油总公司 一种基于平面导航的快速三维层位解释方法
CN106556869A (zh) * 2017-01-19 2017-04-05 中国石油大学(华东) 一种拉张环境下定量刻画断层走滑位移的方法
CN109975872A (zh) * 2019-04-26 2019-07-05 河南理工大学 一种利用Croel DRAW软件直接编制平衡剖面的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
方红萍 等: ""断层识别与定量解释方法进展"", 《工程地球物理学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111856572A (zh) * 2020-07-06 2020-10-30 中国石油大学(北京) 一种断层破碎带宽度的确定方法及装置
CN111856572B (zh) * 2020-07-06 2021-07-20 中国石油大学(北京) 一种断层破碎带宽度的确定方法及装置

Also Published As

Publication number Publication date
CN110441821B (zh) 2020-10-09

Similar Documents

Publication Publication Date Title
CN104850682B (zh) 基于位置的多点地质统计学建模方法
CA2737008C (en) Smoothing of stair-stepped geometry in grids
CN108415077A (zh) 新的边缘检测低序级断层识别方法
CN105701319B (zh) 一种水平井规则开发井网下的沉积微相建模方法
CN104574511A (zh) 一种快速递进的三维地质建模方法
CN101738638B (zh) 石油地震构造图等值线的填充方法
CN102867332B (zh) 基于复杂边界约束的多级细分网格曲面拟合方法
CN104166163A (zh) 基于三维大数据量地震数据体的断层曲面自动提取方法
CN106709987B (zh) 一种三维地质剖面模型动态构建方法
CN103941286A (zh) 一种基于平面导航的快速三维层位解释方法
CN103941288A (zh) 一种基于平面导航的三维目标地质体解释方法
CN103941287A (zh) 一种基于平面导航的快速三维断层解释方法
CN106338761A (zh) 三维镶嵌式实体建模方法
CN102353986B (zh) 含逆断层的复杂地质构造的全三维联动地震解释方法
CN110441821A (zh) 一种基于变解释测网密度的断层快速解释方法
CN102903149B (zh) 地质模型的成块成体方法以及装置
CN112489208A (zh) 基于蚂蚁算法的裂缝片提取方法和三维地质模型构建方法
Taixian et al. Techniques for high-efficient development of offshore fluvial oilfields
CN106652033B (zh) 一种地质剖面自然网格的剖分方法
CN110703334B (zh) 一种复杂断块油田层位快速解释方法
CN110211232A (zh) 一种基于钻孔数据的三棱柱分解建模方法及系统
CN109782344A (zh) 沉积层序边界识别方法及装置
CN108508486A (zh) 一种基于构造演化地层应变的裂缝预测分析方法
CN115035258A (zh) 一种高效的基于cad钻孔柱状图的城市三维地质建模方法
Shao et al. 3D Geological Modeling under Extremely Complex Geological Conditions.

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