[0042] 在步骤102,输入增强地震数据,设定搜索参数,m为蚂蚁的数目,n表示地震采样点,a为采样点的属性值。在算法初始时刻,将m只蚂蚁随机放到n个采样点上,蚂蚁同时从一个采样点到另一个采样点,逐步完成搜索。Nc为算法循环迭代次数,1
[0043]
[0044]
[0045] 其中Sk(i)是k允许集,Sk(i)=C-禁忌集为蚂蚁k下一步允许选择的采样点,C是所有采样点的集合{1,2,3,...,n},α为信息素因子,表示信息素的重要程度,数值越大就体现信息素在蚂蚁行进过程中越重要。ηij表示采样点i,j之间的能见度即属性变化值倒数,ηij=1/(a(j)-a(i)),属性变化越小,能见度越大,被选择的期望越高,β为期望启发式因子,表示能见度的重要程度,数值越大就体现具体问题的先验信息在蚂蚁行进过程中越重要。当全部的蚂蚁完成了一次循环后,信息素有增有减,因此信息素强度需要调整,设信息素挥发系数 每次循环后采样点i,j之间的信息素强度如式(3、4)所示:
[0046] τi,j(t+n)=(1-ρ)*τi,j(t)+Δτi,j(t) (3)
[0047] 式中
[0048]
[0049] Δτkij(t)表示第k只蚂蚁在此次循环时留在采样点i,j之间的信息素量,Δτkij(t)有蚁量、蚁密、蚁周等模型获得。Δτi,j(t)表示此次循环中采样点i,j之间的信息素量增量。
[0050] 基本的蚁群算法实现如下:
[0051] 初始化阶段。令当前循环次数Nc=0,时间t=0;设各个路径上的初始化信息量τi,j(0)为一常量。信息素的增量Δτi,j(t)=0;设定最大循环次数NCmax地方,将m只蚂蚁放置于n个采样点上,初始化参数ηij、α,β。
[0052] 循环阶段。如果NcNCmax,则循环结束。
[0053] 每只蚂蚁k按式(1)移动下一个顶点j,j∈k允许集.
[0054] 将每一只蚂蚁上一步走过的采样点放到该蚂蚁的禁忌表中,更新禁忌表。
[0055] 根据(2)式更新每条路径上信息量。
[0056] Nc=Nc+1,返回步骤(2)。
[0057] 循环结束,算法终止输出处理数据体。
[0058] 流程进入到104。
[0059] 在步骤104,判断波组特征变化是否满足搜索参数的要求。当满足搜索参数要求时,流程进入到步骤105。当不满足搜索参数要求时,流程返回到步骤103,重新设定参数。
[0060] 在步骤105,将满足搜索参数的波组特征变化处作为超剥点,并将其记录下来。流程进入到步骤106。
[0061] 在步骤106,通过步骤104生成的数据体,追踪解释该数据体得到超剥点,将超剥点连成超剥线,输出超剥线,从而得到地层超剥线在平面上的分布。
[0062] 以下为应用本发明的两个验算正演模型和两个具体实施例:
[0063] 正演模型1,超覆型地质模型及超覆点识别
[0064] 建立图2所示的超覆型砂泥岩间互地质模型,模型包括6层层状砂体,厚度5~15m,不整合面倾角5°~45°,砂岩速度2800m/s,泥岩速度2300m/s,不整合面以下地层速度4300m/s。
[0065] 从主频30Hz合成地震记录看,越靠近超覆点层状砂体振幅越弱,识别的超覆点和实际超覆点相距较远,二者相距大小和不整合面倾角、砂体厚度、砂泥岩速度差及主频有关。当不整合面较缓时,合成地震记录得到的超剥点和实际超剥点距离相差较大,如不整合面倾角在5°左右时,二者相距300m;当不整合面较陡时,其合成地震记录反映的超剥点和实际超剥点相差就小,如不整合面倾角在45°左右时,二者相差70m左右。
[0066] 而蚁群算法识别的超覆点位置和实际超覆点位置一般相差较小。当不整合面越缓时,蚁群算法识别的超覆点和实际超覆点距离相差越大,不整合面倾角5°左右时,相距100m。当不整合面越陡时,蚁群算法识别的超覆点和实际超覆点相差越小;不整合面倾角
45°左右时,相差10m左右,二者基本吻合。
[0067] 正演模型2,剥蚀型地质模型及剥蚀点识别
[0068] 建立图3所示的剥蚀型砂泥岩间互地质模型,模型包括5层层状砂体,地层倾角10°,厚度5~50m,砂岩速度2800m/s,泥岩速度2300m/s。上覆地层倾角5°,地层为泥岩和砂岩,砂岩速度2400m/s,泥岩速度2100m/s。
[0069] 从主频30Hz合成地震记录看地震响应,上覆地层和剥蚀地层之间形成中强振幅,砂体地震响应可识别的剥蚀点和实际剥蚀点相距较远。地震响应可识别剥蚀点距离最前端剥蚀点随砂体厚度增大而增大,砂体厚度15m时砂体地震响应可识别的剥蚀点和实际剥蚀点最前端相距约400m;砂体厚度45m时相距约700m。砂体地震响应可识别的剥蚀点和实际最后端剥蚀点距离约200m,和砂体厚度相关性较小。砂体地震响应可识别的剥蚀点与实际最后端剥蚀点垂直高度和地层倾角及上覆地层倾角有关,砂体地层倾角10°,上覆地层倾角5°时,二者垂直距离约20m,砂体地震响应可识别的剥蚀点与实际最前端剥蚀点垂直高度不但和地层倾角有关,还和砂体厚度有关。
[0070] 蚁群算法识别的剥蚀点位置和实际剥蚀点位置相差较小,相差20m以内,二者基本吻合。蚂蚁追踪算法可以准确识别剥蚀砂体的最前端剥蚀点,而合成地震记录不能识别。
[0071] 实施例1,草桥北坡沙河街组地层超剥带识别
[0072] 草桥北坡位于济阳坳陷东营凹陷南部斜坡带,自南向北沙河街组地层逐渐遭受剥蚀,被馆陶组地层覆盖,发育以地层圈闭为主的多类型圈闭,具有较大勘探潜力。草桥北坡沙河街组Ⅰ、Ⅱ、Ⅲ级不整合面附近超覆、剥蚀现象明显,但地层超剥点常规地震资料较难识别,尤其是Ⅱ、Ⅲ级不整合面。该地区Ⅲ级不整合面由沉积间断面(或侵蚀面)和相转换面(由湖平面突变造成的沉积体系域转换)组成。分别位于孔店组-沙四下内部(T8—Ek2/Ek1、T8’—Ek1/Es4下)、沙四上与东营组之间(T7—沙四上内部纯上/纯下、T6—沙三段底界面、T4—Es3中/Es3上、T2’—Es1/Ed)。从地震剖面上看这些不整合面有的表现为强振幅、高连续反射特征,但由于受到地震资料分辨率限制(频率在25Hz左右),尖灭点附近表现为复合地震波反射,实际地层尖灭位置不能准确刻画;有的不整合面在强反射屏蔽下表现为空白弱反射,追踪困难。使用蚁群算法对这些不整合面进行处理,由于蚁群追踪最大限度保留了地震横向连续性信息,处理后地层超剥点较地震剖面有明显提升,原来复合波有效分离,超覆地层尖灭点表现清晰,剥蚀点与不整合面接触关系更为明确,尖灭点位置非常清楚,相对于直接在地震剖面上识别,在生物启发计算剖面上向前推50米左右,更易于划分地层结构,如图4所示。生物启发三维显示图上可以从不同方向进行超剥线的识别,可以看出各层系地层遭受剥蚀现象清晰,超剥点明确,如图5所示。
[0073] 实施例2,义和庄凸起中古生界地层超剥带识别
[0074] 义和庄凸起位于济阳坳陷车镇凹陷南部沾化凹陷西部,发现中生界和上下古生界多套含油层系,中生界和上古生界石炭-二叠系为首构造控制的层状油藏,以断块和地层不整合油藏为主。由于多期构造运动影响,中-古生界潜山内幕构造复杂,断层刻画难度大,地震反射特征比较杂乱,与上覆第三系地层接触关系不明确,地层结构不明确,构造和地层超剥带落实困难。利用生物启发计算处理后,可以明显看到中-古生界发育高角度层状地层,地层呈平行或亚平行反射特征,北厚南薄,发育两条较大的南倾断层,中-古生界地层与上覆地层接触关系清楚,各层剥蚀点(尖灭点)位置明确,使用生物启发计算参数剖面更易于刻画地层结构,如图6所示。