CN100349010C - 三维地质层位自动追踪方法 - Google Patents

三维地质层位自动追踪方法 Download PDF

Info

Publication number
CN100349010C
CN100349010C CNB2004101026486A CN200410102648A CN100349010C CN 100349010 C CN100349010 C CN 100349010C CN B2004101026486 A CNB2004101026486 A CN B2004101026486A CN 200410102648 A CN200410102648 A CN 200410102648A CN 100349010 C CN100349010 C CN 100349010C
Authority
CN
China
Prior art keywords
seed points
tomography
impact point
tracking
bounding box
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.)
Expired - Fee Related
Application number
CNB2004101026486A
Other languages
English (en)
Other versions
CN1797039A (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.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
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 China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CNB2004101026486A priority Critical patent/CN100349010C/zh
Publication of CN1797039A publication Critical patent/CN1797039A/zh
Application granted granted Critical
Publication of CN100349010C publication Critical patent/CN100349010C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

一种三维地质层位自动追踪方法,在三维地震数据体上定义种子点;定义追踪范围;定义层位自动追踪拾取位置参数;选择拾取法;选择追踪法;选择断层,并拟合成断层面;确定第一个当前工作种子点,在其相邻的线和道上寻找目标点;目标点与种子点对比,得到符合层位自动追踪拾取方法的点,并确定新的种子点,用替代手工拾取层位,提高追踪质量,能在整个地震数据上进行的三维地质层位自动追踪。

Description

三维地质层位自动追踪方法
技术领域
本发明涉及地球物理勘探地震数据处理解释技术,是一种三维地质层位自动追踪方法。
背景技术
在得到地震数据后,需要在地震数据生成的图中搜寻出地质层位的位置与走向,然后显示出三维断层的形态与地质体结构。
在以往的工作中,通常是由解释员在剖面上手工逐一标注点位、层位,对于大量的三维地震数据,有时也采用空间层位自动追踪的方法,这样可以快速处理大量地震数据、搜寻地质层位。
实际的地震数据通常存在大量断层。当层位遇到断层时,同一地质层位被断层切割、错断。以往在没有断层约束的情况下使用自动追踪方式,当遇到断层时,往往会穿越断层,错误地追踪到其他的层位。因此空间层位自动追踪的使用范围受到限制,追踪质量大为降低,使之难以广泛使用。
发明内容
本发明所要解决的技术问题是提供一种替代手工拾取层位、提高追踪质量、可在整个地震数据中进行三维地质层位自动追踪的方法。
本发明的三维地质层位自动追踪方法,对采集的地震数据进行包含以下步骤的层位追踪:
(1)在三维地震数据体中定义种子点;
(2)定义追踪范围;
(3)定义层位自动追踪拾取位置参数;
(4)选择拾取方法;
(5)选择追踪方法;
(6)选择断层并拟合断层面;
(7)确定第一个当前工作种子点,在其相邻的线与道中寻找目标点;
(8)目标点与种子点对比,得到符合层位自动追踪拾取方法的点,并确定新的种子点;
(9)选用未使用的种子点作为新的种子出发点,重复步骤(7)、(8)进行新一轮追踪,直至追踪范围内的所有种子点皆被追踪,得到地震地质层位。
本发明步骤(1)、(7)、(8)中定义或确定的种子点在同一时间方向只使用一次
本发明确定种子点的方法为用户使用鼠标在剖面上直接定义或利用已经解释过的层位数据作为种子点。
本发明所述的拾取位置为波峰或波谷或过零点值(Zero(-/+),Zero(+/-))。
所述的拾取方法为相关比较算法与振幅比较算法的组合。
所述的追踪方法为简单方法、两点校验方法、四点校验方法:
其中简单方法:判断得到的符合层位自动追踪拾取方法的目标点是否落在断层上,不在断层上的目标点成为新的种子点,放入种子队列中。
两点校验方法:判断得到的符合层位自动追踪拾取方法的目标点是否落在断层上,不在断层上的目标点与原种子点再进行比较,判断是否符合层位自动追踪拾取方法,若符合,则目标点成为新的种子点,放入种子队列中。
四点校验方法:判断得到的符合层位自动追踪拾取方法的目标点是否落在断层上,不在断层上的目标点,在其相邻的线与和道中寻找第二个目标点,同上判断,直至第三个目标点,然后第三个目标点与原种子点再进行比较,判断是否符合层位自动追踪拾取方法,若符合,则目标点成为新的种子点,放入种子队列中,否则三个目标点都不能成为种子点。
本发明所述的目标点符合未被追踪过并符合层位自动追踪条件。
本发明选择断层并自动拟合断层面的方法包括:
(1)用户选取作为断层阻挡的断层名。
(2)筛选断层数据,去除只有断层名而没有解释数据的断层,去除无法构建断面的断层,如一个断层段无法构建成面进行约束;
(3)把筛选后的每一个断层自动拟合成断层面。
一个断层由若干断层段(折线段)组成,使用断层做阻挡,必须把断层用三角形拼接成面。
本发明判断目标点是否在断层面上时,首先计算种子点与目标点构成的包围盒与一个断层面构成的包围盒是否相交,如果不相交,继续与下一个断层面构成的包围盒判断比较;如果相交进入下一步判断,判断种子点与目标点的包围盒是否与断层面上的一个三角形构成的包围盒相交,如果不相交,继续与下一个三角形构成的包围盒判断,如果相交,进入第三步判断交点是否在种子点与目标点的连线间并同时在三角形内。如果是,则判断目标点与该断层面相交,自动追踪停止,从种子队列中选取下一个种子点进行追踪;如果不是,则目标点不在断层面上,目标点被拾取成为新的种子点。
所述种子点与目标点构成的包围盒,是种子点与目标点构成连线段,比较这两点的坐标,得到最大与最小的x,y,z坐标值,由此构建一个可把种子点与目标点构成的连线段完全包含在内的长方体。
本发明的用于地震数据处理的空间层位自动追踪方法用软件完成快速、精确的空间层位追踪,提高追踪质量,从而提高工作效率。
附图说明
图1是种子点与目标点构成的包围盒示意图,
1-种子点,2-目标点;
图2是断层面构成的包围盒示意图;
图3是三角形构成的包围盒示意图;
图4是使用断层约束的效果示意图,3-层位,4-断层,由于受到断层4的阻挡,两个层位3在断层处断裂、错位;
图5a、图5b是未使用断层约束的地震数据层位追踪效果示意图,可以清楚地观察到,遇到断层的地方由于没有断层约束,层位没有被准确拾取到,而是被错误地拾取到黑色箭头所示的位置;
图6是断层用三角形拼接的方法拟和成断层面的效果示意图;
图7是在三维地震数据的一张剖面中拾取种子点的效果示意图;
图8是种子点与相邻目标点的示意图,三角形代表目标点,圆点代表种子点,箭头代表种子点与目标点比较;
图9是种子点与相邻目标点使用两点校验的示意图,三角形代表目标点,圆点代表种子点,细箭头代表种子点与目标点比较,粗箭头代表目标点与种子点校验;
图10是种子点与相邻目标点使用四点校验的示意图,三角形代表目标点,圆点代表种子点,细箭头代表种子点与目标点比较,粗箭头代表目标点与种子点校验。
具体实施方式
实施例1
(1)在采集的地震数据中用鼠标沿着Line方向的第720线剖面中定义了3个种子点,如图7所示5、6、7三点,存贮于种子点队列中,存贮的顺序与拾取种子点的顺序一致。种子点队列是线性队列,如下图说明种子点的存放。
第一个拾取的种子点 第二个拾取的种子点 第三个拾取的种子点
(2)定义追踪范围,采用的实际体数据是Line方向:从450线到790线;CDP方向:从350道到500道;Time方向:从1350ms到1600ms,定义的追踪范围只能在Line,CDP范围之内。根据实际数据,选择追踪全部Line,CDP范围。
(3)定义层位自动追踪参数,其中拾取位置为波峰,拾取方法为相关比较算法与振幅比较算法,追踪方法为简单方法。
(4)选取作为断层yt1作为阻挡的断层,筛选断层数据,去除没有实际数据的断层,去除无法构建断面的断层,如一个断层段无法构建成面进行约束,把筛选后的每一个断层自动拟合成断层面,拟合效果如图6所示。
(5)把种子队列中序号为1的种子点作为第一个当前工作种子点,把其在相邻的线与道中的4个点作为目标点,如图8所示。
(6)取4个目标点中的一个,首先判断此目标点是否被拾取过,如果被拾取过,此点被放弃,不能成为种子点;如果没有拾取过,接着判断使用相关算法与振幅比较算法计算得到的相关系数与振幅比较系数是否满足(即大于等于)用户给定的相关系数与振幅比较系数。如果有一个条件不满足,此目标点不能被成为种子点,同时记录此点已经被拾取过,防止再次计算。接着取剩下三个目标点中的一个重复上述操作;如果此目标点同时满足给定的条件,采用简单追踪方法,判断得到的符合层位自动追踪拾取方法的目标点是否落在断层上。
首先计算种子点与目标点构成的包围盒与一个断层面构成的包围盒是否相交,如果不相交,继续与下一个断层面构成的包围盒判断比较;如果相交进入下一步判断,判断种子点与目标点的包围盒是否与断层面上的一个三角形构成的包围盒相交,如果不相交,继续与下一个三角形构成的包围盒判断,如果相交,进入第三步,判断交点是否在种子点与目标点的连线间并同时在三角形内,如果是,则判断目标点与该断层面相交,否则目标点不在断层面上。如果目标点与断层面相交,这个目标点便不能成为种子点,同时将已经被拾取过不在断层上的目标点变成新的种子点,放入种子队列中。取剩余两个目标点中的一个重复上述步骤,直至相邻的4个目标点都如上述步骤判断过。第一个种子点与相邻目标点判断过后,从种子点队列剔除刚才工作的种子点,种子点队列的情况如下:
第二个种子点 第三个种子点 第四个种子点 .......
(7)选用第二个种子点作为新的种子出发点,重复步骤(5)、(6)进行新一轮追踪,直至种子点队列中的种子点全部追踪完毕,停止追踪,得到地震数据中的地质层位。
实施例2
(1)在采集的地震数据中利用已经解释过的层位数据作为种子点,存贮于种子点队列中。
(2)定义追踪范围。
(3)定义层位自动追踪参数,其中拾取位置为波谷,拾取方法为相关比较算法与振幅比较算法,追踪方法为两点校验方法。
(4)用户选取作为断层阻挡的断层名,筛选断层数据,去除没有实际数据的断层,去除无法构建断面的断层,把筛选后的每一个断层自动拟合成断层面。
(5)把种子队列中序号为1的种子点作为第一个当前工作种子点,在其相邻的线与道中寻找目标点;
(6)未被追踪过并符合层位追踪条件的目标点与种子点对比,得到符合层位自动追踪拾取方法的点,接着判断目标点能否变成新的种子点。
采用两点校验的追踪方法,判断得到的符合层位自动追踪拾取方法的目标点是否落在断层上同实施例1,不在断层上的目标点与原种子点再进行比较,如图9所示,判断是否符合层位自动追踪拾取方法,若符合,则目标点变成新的种子点,放入种子队列中。
(7)选用未使用的种子点作为新的种子出发点,重复步骤(5)、(6)进行新一轮追踪,直至追踪范围内的所有种子点皆被追踪过,得到地震数据中含有断层的地质层位。
实施例3
(1)在采集的地震数据中利用已经解释过的层位数据作为种子点,存贮于种子点队列中。
(2)定义追踪范围。
(3)定义层位自动追踪参数,其中拾取位置为过零点值,拾取方法为相关比较算法与振幅比较算法,追踪方法为四点校验方法。
(4)用户选取作为断层阻挡的断层名,筛选断层数据,去除没有实际数据的断层,去除无法构建断面的断层,把筛选后的每一个断层自动拟合成断层面。
(5)把种子队列中序号为1的种子点作为第一个当前工作种子点,在其相邻的线与道中寻找目标点;
(6)未被追踪过并符合层位追踪条件的目标点与种子点对比,得到符合层位自动追踪拾取方法的点,接着判断目标点能否变成新的种子点。
采用四点校验的追踪方法,判断得到的符合层位自动追踪拾取方法的目标点是否落在断层上同实施例1,不在断层上的目标点,在其相邻的线与道中寻找第二个目标点,同上判断,直至第三个目标点,然后第三个目标点与原种子点再进行比较,如图10所示,判断是否符合层位自动追踪拾取方法,若符合,则目标点变成新的种子点,放入种子队列中,否则三个目标点都不能变成种子点。
(7)选用未使用的种子点作为新的种子出发点,重复步骤(5)、(6)进行新一轮追踪,直至追踪范围内的所有种子点皆被追踪过,得到地震数据中需要的地质层位。

Claims (10)

1、一种三维地质层位自动追踪方法,对采集到的地震数据进行包含以下步骤的层位追踪:
(1)在三维地震数据体上定义种子点,在采集到的地震数据中利用已经解释过的层位数据作为种子点,存贮于种子里队列中;
(2)定义追踪范围;
(3)定义层位自动追踪拾取位置参数;
(4)选择拾取法;
(5)选择追踪法;
(6)选择断层,并拟合成断层面;
(7)确定第一个当前工作种子点,在其相邻的线和道上寻找目标点;
(8)目标点与种子点对比,得到符合层位自动追踪拾取方法的点,并确定新的种子点;
(9)选用未使用的种子点作为新的种子出发点,重复步骤(7),(8)进行新一轮的追踪,直至追踪范围内的所有种子点皆被追踪,得到地震地质层位。
2、根据权利要求1所述的三维地质层位自动追踪方法,其特征在于步骤(1)、(7)、(8)中定义或确定的种子点在同一时间方向只使用一次
3、根据权利要求1所述的三维地质层位自动追踪方法,其特征在于步骤(1)定义种子点是根据已知层位数据在剖面上直接定义。
4、根据权利要求1所述的三维地质层位自动追踪方法,其特征在于所述的拾取位置参数为波峰或波谷或过零点值。
5、根据权利要求1所述的三维地质层位自动追踪方法,其特征在于所述的拾取法为相关比较算法和振幅比较算法的组合。
6、根据权利要求1所述的三维地质层位自动追踪方法,其特征在于所述的追踪法为简单法或两点校验法或四点校验法。
7、根据权利要求1所述的三维地质层位自动追踪方法,其特征在于步骤(6)选择断层并拟合成断层面为:
选取作为断层阻挡的断层;
筛选断层数据,去除没有解释数据的的断层,去除无法构建断面的断层;
把筛选后的每一个断层拟和成断层面。
8、根据权利要求1或6所述的三维地质层位自动追踪方法,其特征在于步骤(5)追踪法在判断目标点是否在断层面上,采用计算包围盒与一个断层面构成的包围盒是否相交,不相交继续和下一个断层面构成的包围盒比较判断。
9、根据权利要求8所述的三维地质层位自动追踪方法,其特征在与断层面构成的包围盒相交的种子点和目标点的包围盒与该断层面上的一个三角形构成的包围盒是否相交,不相交继续和下一个三角形的包围盒比较。
10、根据权利要求9所述的三维地质层位自动追踪方法,其特征在与三角形构成的包围盒相交的种子点和目标点的包围盒,判断交点是否在种子点和目标点的连线间并在三角形内,如果是,则判断目标点和某一个断层面相交,否则目标点不在断层面上。
CNB2004101026486A 2004-12-29 2004-12-29 三维地质层位自动追踪方法 Expired - Fee Related CN100349010C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2004101026486A CN100349010C (zh) 2004-12-29 2004-12-29 三维地质层位自动追踪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2004101026486A CN100349010C (zh) 2004-12-29 2004-12-29 三维地质层位自动追踪方法

Publications (2)

Publication Number Publication Date
CN1797039A CN1797039A (zh) 2006-07-05
CN100349010C true CN100349010C (zh) 2007-11-14

Family

ID=36818242

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2004101026486A Expired - Fee Related CN100349010C (zh) 2004-12-29 2004-12-29 三维地质层位自动追踪方法

Country Status (1)

Country Link
CN (1) CN100349010C (zh)

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101702237B (zh) * 2009-11-10 2011-11-23 南京大学 一种基于二次误差测度的层位面拟合方法
US8849574B2 (en) * 2011-07-28 2014-09-30 Chevron U.S.A. Inc. System and method for perturbing an initial horizon-picking solution to follow local features of a volume
CN102819040B (zh) * 2012-08-22 2014-12-10 电子科技大学 基于中心扩散加倾角属性的三维地震层位自动追踪方法
CN102819688A (zh) * 2012-08-29 2012-12-12 电子科技大学 基于半监督分类的二维地震数据全层位追踪方法
US10577895B2 (en) 2012-11-20 2020-03-03 Drilling Info, Inc. Energy deposit discovery system and method
CN103869356B (zh) * 2012-12-10 2016-10-26 中国石油天然气集团公司 一种地震反射界面闭合方法
CN104062681B (zh) * 2013-03-22 2016-12-28 中国石油化工股份有限公司 一种基于分数阶导数的地震层位追踪预处理方法
US10853893B2 (en) 2013-04-17 2020-12-01 Drilling Info, Inc. System and method for automatically correlating geologic tops
US10459098B2 (en) * 2013-04-17 2019-10-29 Drilling Info, Inc. System and method for automatically correlating geologic tops
CN104375175B (zh) * 2013-08-15 2017-08-04 中国石油天然气集团公司 倾角传播法层位自动追踪方法
CN103412331B (zh) * 2013-08-30 2015-10-28 电子科技大学 一种三维地震断层自动提取方法
CN103592681B (zh) * 2013-09-16 2016-05-04 电子科技大学 一种基于信号分类的地震图像层位追踪方法
CN104597494B (zh) * 2013-10-31 2017-06-13 中国石油天然气集团公司 地震地层体分析方法及装置
CN104656131A (zh) * 2013-11-22 2015-05-27 中国石油天然气集团公司 一种基于波形恢复技术的提高地震层位自动追踪精度的方法
CN103901467A (zh) * 2014-03-18 2014-07-02 中国石油集团川庆钻探工程有限公司地球物理勘探公司 三维地震数据层位追踪的方法
CN104181596B (zh) * 2014-08-27 2017-01-11 中国石油集团东方地球物理勘探有限责任公司 地质层位自动追踪方法及装置
CN104199092A (zh) * 2014-08-31 2014-12-10 电子科技大学 基于多层次框架的三维全层位自动追踪方法
CN105022092B (zh) * 2014-09-15 2017-11-28 中国矿业大学(北京) 用于确定断层展布格局的地震资料解释方法和设备
US10908316B2 (en) 2015-10-15 2021-02-02 Drilling Info, Inc. Raster log digitization system and method
CN111796324B (zh) 2019-04-09 2023-02-10 中国石油天然气股份有限公司 地震全层位追踪方法及装置
CN112305596A (zh) * 2019-07-26 2021-02-02 中国石油天然气集团有限公司 断层控制下的层位自动追踪方法及装置
CN112068196B (zh) * 2020-09-07 2021-11-16 北京中恒利华石油技术研究所 一种网格化层位追踪方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5537365A (en) * 1993-03-30 1996-07-16 Landmark Graphics Corporation Apparatus and method for evaluation of picking horizons in 3-D seismic data

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5537365A (en) * 1993-03-30 1996-07-16 Landmark Graphics Corporation Apparatus and method for evaluation of picking horizons in 3-D seismic data

Also Published As

Publication number Publication date
CN1797039A (zh) 2006-07-05

Similar Documents

Publication Publication Date Title
CN100349010C (zh) 三维地质层位自动追踪方法
US6823266B2 (en) Method for performing object-based connectivity analysis in 3-D seismic data volumes
US6516274B2 (en) Method for imaging discontinuities in seismic data using dip-steering
EP0652447B1 (en) Method and apparatus for finding horizons in 3D seismic data
US5056066A (en) Method for attribute tracking in seismic data
CN105787281B (zh) 一种海浪有效波高场的融合反演方法和装置
US5432751A (en) Method and apparatus for finding horizons in 3D seismic data
CA3092055C (en) Seismic data acquisition using designed non-uniform receiver spacing
CN102395902B (zh) 使用快速面向目标照明计算的地震成像系统及方法
US5757663A (en) Hydrocarbon reservoir connectivity tool using cells and pay indicators
CN104375175B (zh) 倾角传播法层位自动追踪方法
CN103969682B (zh) 一种钻测井资料深度匹配方法及系统
WO2010096759A1 (en) Dynamic extraction of extrema-based geometric primitives in 3d seismic voxel volume
CN104597494B (zh) 地震地层体分析方法及装置
CN105319576B (zh) 一种凹多边形区域的观测系统自动满覆盖快速布设方法
CN110727030B (zh) 一种基于子波多极值特征联合分析的储层预测方法
CN110967737B (zh) 一种构造约束的初始模型构建方法
AU2002312550B2 (en) Method for performing object-based connectivity analysis in 3-D seismic data volumes
GB2583910A (en) Method of analysing seismic data
CN108051858A (zh) 优化的地震勘探方法及装置
CN109991663A (zh) 工区地震速度体校正方法和装置
AU2002312550A1 (en) Method for performing object-based connectivity analysis in 3-D seismic data volumes
CN112363221B (zh) Walkaway VSP测线的布设方法及装置
US20220196863A1 (en) Automated horizon layer extraction from seismic data for wellbore operation control
CN117270042A (zh) 走滑断层的确定方法及确定装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20071114