主权项 |
一种介质品质因子的逐次差分进化估计方法,其特征在于,包括如下步骤,步骤1,用WEPIF方法估计实际VSP道集的品质因子和速度前,先对实际VSP道集进行预处理,去除不同炮子波的不一致性;步骤2,构建一个含M个个体的初始群体X,表示为,X=[x<sub>1</sub>,x<sub>2</sub>,…,x<sub>M</sub>];群体中的第k个个体表示为,x<sub>k</sub>=[x<sub>k,1</sub>,x<sub>k,2</sub>,…,x<sub>k,N</sub>],k=1,…,M;其中,N为总道数,M为地层数量,每个地层对应一组要进行估计的品质因子和速度;步骤3,按相邻两个检波器的距离将地下介质划分为不同地层,将每一道作为一个子成分,用L<sub>1</sub>范数定义第k个个体的局部适应度函数为,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>f</mi><mrow><mi>k</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mi>P</mi><mo>-</mo><msub><mi>P</mi><mn>1</mn></msub></mrow><mrow><mi>P</mi><mo>+</mo><msub><mi>P</mi><mn>2</mn></msub></mrow></munderover><mo>|</mo><mo>[</mo><msub><mi>d</mi><mrow><mi>o</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>d</mi><mrow><mi>k</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>]</mo><mo>·</mo><mi>g</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>|</mo><mo>;</mo></mrow>]]></math><img file="FDA0000687293200000011.GIF" wi="789" he="169" /></maths>从而得到每道的局部适应度值;式中,d<sub>o,j</sub>(·)为第j道观测数据,d<sub>k,j</sub>(·)为由第k个个体正演计算的第j道数据,P为观测数据直达波瞬时振幅谱的峰值位置,P<sub>1</sub>、P<sub>2</sub>分别为子波瞬时振幅的包络峰值左右两侧的第一个波谷位置,g(n)是高斯窗函数;步骤4,根据实际VSP道集和每道的局部适应度,对初始群体X中的每个子成分进行优化更新;进行正演模拟逐道的优化每个子成分,当所有个体中当前层和下一层的参数均选择最优值后,再进行层位循环,逐层反演参数;反演第j层的参数时,将群体X中第j列的所有变量进行升序排列,根据变异策略求取变异后的待估品质因子和速度;步骤5,重复步骤4,将整个群体更新完后根据全局适应度的值选择最优个体,从而得到地下介质的品质因子和速度。 |