[发明专利]基于遗传算法的剩余静校正求取方法

说明书

技术领域

本发明涉及一种基于遗传算法的剩余静校正求取方法，该方法可自动、 科学确定遗传算法计算参数，应用于石油地震勘探，属于地震勘探资料处理 领域，是一种获取精确的高频静校正量的有效方法。

背景技术

我国年石油天然气对外依存度极高，而我国剩余油气资源超过一半聚集 在山地复杂构造中，及时发现探明山地复杂构造油气资源迫在眉睫。这些地 区地震地表条件复杂，静校正问题严重，导致多次覆盖资料同一反射点因为 地震道之间存在延迟时差，叠加成像困难，其精度难以满足后期的构造、储 层解释，由这类品质的资料提供的建议井位可靠性不高，将会增加油气勘探 的成本和风险。

为降低这类因素带来的成本和风险，物探人员尝试采用遗传算法求剩余 静校正量，减少共反射点道集的延迟时差，其基本思想是从模型群体开始搜 索，把模型参数编码，将模型空间的点映射到染色体空间的染色体，然后通 过选择、交换和变异等遗传操作对模型群体进行繁殖，逐次迭代，最后求取 非线性反演问题整体极值所对应的最优解或近似最优解。

然而，常规的遗传算法求取剩余静校正的实际性能受控制参数(如种群数 目、种群大小、变异概率等)的影响，如何确定这些控制参数成为一直困扰物 探人员的瓶颈问题，也直接制约了该技术在生产上的推广应用。

发明内容

为解决常规遗传算法参数难以确定的问题，申请人提出一种基于遗传算 法的剩余静校正求取方法，可以根据允许的最大时移量计算出种群大小、种 群代数和个体变异率，自动计算出剩余静校正量。

本发明提供了一种基于遗传算法的剩余静校正求取方法。所述方法用于 处理二维地震资料，其包括以下步骤：

A、将不同的炮点或检波点的真实值分别作为遗传算法的不同个体中的 单个基因，将多个个体构成初始种群，种群大小popSize通过式1确定，式1 为：popSize＝(maxSta/sampInv)×f_popSize，式1中maxSta为最大时移量， 在10～60ms之间选择；sampInv为地震数据的采样间隔，单位为ms；f_popSize为 种群大小的调控因子，在0.5～1.5之间取值；

B、计算种群中每个个体的适应度值，随后对种群中的所有个体按适应度 值由大到小排序，选择保留适应度值最高的个体作为交换对象；

C、在种群的第k个个体和所述交换对象上各随机产生两个交换位置，将 所述第k个个体的两个交换位置之间的基因段与所述交换对象的两个交换位 置之间的基因段进行交换，产生子代新个体，其中，k＝popSize×0.4，k须 取整；

D、对步骤C中所述进行交换后获得的种群，保留适应度值最高的个体作为精英进入子代，对于其余个体，若随机数r大于该个体的变异概率，则该个体随机产生两个变异位置，并将该个体的两个变异位置处的基因随机变为允许移动的最大时移量maxSta范围内的随机值，从而产生子代新个体，若随机数r不大于该个体的变异概率，则不变异，同时该个体直接进入子代作为子代个体，其中，随机数r在[0，1]间随机产生，个体的变异概率通过以下方式获得，先对该个体所在种群中的所有个体按适应度值由大到小排序，排序后，第i个个体的变异概率根据式2得到，式2为：