主权项 |
一种消除速度误差影响的震源定位方法,其特征在于:所述方法通过已知震源的位置,求取一个与速度相关的速度误差模型,把速度误差模型应用到目标函数中,然后扫描所有可能的震源位置,找到使目标函数的值最小的点的坐标就是待定位震源的坐标;所述求取一个与速度相关的速度误差模型是这样实现的:设有n个检波点,取已知震源A(x<sub>a</sub>,y<sub>a</sub>,z<sub>a</sub>),通过公式(3)计算震源A到每个检波点的计算旅行时为T<sub>i</sub><sup>A</sup>,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><msub><mi>T</mi><mi>i</mi></msub><mi>A</mi></msup><mo>=</mo><mfrac><msubsup><mi>l</mi><mi>i</mi><mi>A</mi></msubsup><mi>v</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001005762010000011.GIF" wi="1189" he="199" /></maths>其中,T<sub>i</sub><sup>A</sup>为震源A到第i个检波点的计算旅行时,<img file="FDA0001005762010000013.GIF" wi="69" he="103" />为震源A到第i检波点的距离,v为传播速度,那么i检波点处的速度误差模型为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>ϵ</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msup><msub><mi>T</mi><mi>i</mi></msub><mi>A</mi></msup></mrow><msubsup><mi>t</mi><mi>i</mi><mi>A</mi></msubsup></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001005762010000012.GIF" wi="963" he="253" /></maths>其中<img file="FDA0001005762010000014.GIF" wi="54" he="78" />为震源A传播到第i个检波点的观测旅行时,ε<sub>i</sub>为i检波点处的速度误差比例,所述把速度误差模型应用到目标函数中,然后扫描所有可能的震源位置,找到使目标函数的值最小的点的坐标就是待定位震源的坐标是这样实现的:对于任意一个待定位震源进行如下处理:(1)从该震源到各个检波点的观测旅行时中找到最小观测旅行时,该最小观测旅行时所对应的检波点设为检波点K,然后计算出该震源到每个检波点的观测旅行时与该震源到检波点K的观测旅行时之间的时间差,该时间差为观测旅行时差;(2)在给定的范围内,根据所述速度误差模型和公式(5),计算出每个震源到各个检波点的计算旅行时,并求取每个震源到各个检波点的计算旅行时与每个震源到检波点K的计算旅行时之间的时间差,该时间差为计算旅行时差,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>T</mi><mi>i</mi></msub><mo>=</mo><mfrac><msub><mi>l</mi><mi>i</mi></msub><mrow><msub><mi>vϵ</mi><mi>i</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001005762010000021.GIF" wi="990" he="175" /></maths>其中,T<sub>i</sub>表示震源到第i个检波点的计算旅行时,l<sub>i</sub>表示震源到第i个检波点的距离;(3)以步骤(1)得到的观测旅行时差和步骤(2)得到的计算旅行时差之间的差的平方和为目标函数,遍历给定范围内的所有点,使所述目标函数的值最小的点的坐标即为待定位震源的坐标。 |