技術特征：

1.一種雙對雙差地震定位方法，其特征在于，所述方法包括：

基于地震的初至波絕對到時數據或者地震波形的互相關數據，按預設規則，提取觀測的地震事件對(i,j)到同一臺站的到時差數據以及同一地震到臺站對(k,l)的到時差數據，進而提取觀測的雙對到時差數據；所述雙對為所述地震事件對(i,j)和所述臺站對(k,l)的組合，其中，i和j分別表示地震事件對中的兩個地震事件，地震事件i和地震事件j的位置相鄰，l和k分別表示臺站對中的兩個臺站；

計算所述地震事件對(i,j)到臺站k的事件對雙差和所述地震事件對(i,j)到臺站l的事件對雙差dr_l^ij；所述事件對雙差為所述地震事件對(i,j)的到時差數據的觀測值與預測值的差；

或計算所述地震事件i到臺站對(k,l)的臺站對雙差和所述地震事件j到所述臺站對(k,l)的臺站對雙差所述臺站對雙差為所述臺站對(k,l)的到時差數據的觀測值與預測值的差；

計算與dr_l^ij的差值，或計算與的差值，得到所述地震事件對(i,j)到所述臺站對(k,l)的雙對雙差所述雙對雙差為所述雙對的到時差數據的觀測值與預測值的差；

根據臺站對雙差與地震位置相對于地震初始位置的變化量的線性關系，以及雙對雙差與地震位置相對于地震初始位置的變化量的線性關系，構建反演系統；

采取分級加權的方式來反演得到各個地震事件的地震位置。

2.根據權利要求1所述的方法，其特征在于，所述計算所述地震事件對(i,j)到臺站k的事件對雙差和所述地震事件對(i,j)到臺站l的事件對雙差dr_l^ij，包括：

根據所述地震事件對(i,j)中地震事件i和地震事件j到臺站k的到時和的觀測值與預測值的差，得到所述地震事件i和所述地震事件j到臺站k的到時殘差和并根據所述地震事件i和所述地震事件j到所述臺站k的到時殘差和得到所述地震事件對(i,j)到所述臺站k的事件對雙差所述地震事件對(i,j)為多個地震事件對中的任意一個地震事件對，

其中，dx為地震位置的變化量，dτ為發震時刻的變化量，為觀測的事件對(i,j)到時差數據，為預測的事件對(i,j)到時差數據；

根據所述地震事件i和所述地震事件j到臺站l的到時T_lⁱ和T_l^j的觀測值與預測值的差，得到所述地震事件i和所述地震事件j到臺站l的到時殘差r_lⁱ和r_l^j，并根據所述地震事件i和所述地震事件j到所述臺站l的到時殘差r_lⁱ和r_l^j，得到所述地震事件對(i,j)到所述臺站l的事件對雙差dr_l^ij，

其中，dx為地震位置的變化量，dτ為發震時刻的變化量，(T_lⁱ-T_l^j)^obs為觀測的事件對(i,j)到時差數據，(T_lⁱ-T_l^j)^cal為預測的事件對(i,j)到時差數據。

3.根據權利要求1所述的方法，其特征在于，所述計算所述地震事件i到臺站對(k,l)的臺站對雙差和所述地震事件j到所述臺站對(k,l)的臺站對雙差包括：

根據所述地震事件i到所述臺站k和所述臺站l的到時和T_lⁱ的觀測值與預測值的差，得到所述地震事件i到所述臺站k和所述臺站l的到時殘差和r_lⁱ，并根據所述地震事件i到所述臺站k和所述臺站l的到時殘差和r_lⁱ，得到地震i到所述臺站對(k,l)的臺站對雙差所述臺站對(k,l)為多個臺站對中的任意一個臺站對；

其中，s_k為所述臺站k的校正項，s_l為所述臺站l的校正項，為觀測的臺站對(k,l)到時差數據，為預測的臺站對(k,l)到時差數據；

根據所述地震事件j到所述臺站k和所述臺站l的到時和T_l^j的觀測值與預測值的差，得到所述地震事件j到所述臺站k和所述臺站l的到時殘差和r_l^j并根據所述地震事件j到所述臺站k和所述臺站l的到時殘差和r_l^j，得到地震j到所述臺站對(k,l)的臺站對雙差

其中，s_k為所述臺站k的校正項，s_l為所述臺站l的校正項，為觀測的臺站對(k,l)到時差數據，為預測的臺站對(k,l)到時差數據。

4.根據權利要求1所述的方法，其特征在于，計算與dr_l^ij的差值，或計算與的差值，得到所述地震事件對(i,j)到所述臺站對(k,l)的雙對雙差包括：

計算地震事件對(i,j)到臺站k的事件對雙差與地震事件對(i,j)到臺站l的事件對雙差dr_l^ij的差值，或計算地震事件i到臺站對(k,l)的臺站對雙差與地震事件j到臺站對(k,l)的臺站對雙差的差值，得到所述地震事件對(i,j)到所述臺站對(k,l)的雙對雙差所述地震事件對(i,j)與所述臺站對(k,l)的雙對組合為多個雙對中的任意一個雙對；

$\begin{array}{c}{\mathrm{dr}}_{kl}^{ij}={\{(T_k^i-T_k^j)-(T_l^i-T_l^j)\}}^{obs}-{\{(T_k^i-T_k^j)-(T_l^i-T_l^j)\}}^{cal}\\ ={\{(T_k^i-T_l^i)-(T_k^j-T_l^j)\}}^{obs}-{\{(T_k^i-T_l^i)-(T_k^j-T_l^j)\}}^{cal}\end{array},$

其中，和為雙對到時差數據。

5.根據權利要求1所述的方法，其特征在于，所述根據臺站對雙差與地震位置相對于地震初始位置的變化量的線性關系，以及雙對雙差與地震位置相對于地震初始位置的變化量的線性關系，構建反演系統，包括：

所述臺站對雙差與地震位置相對于地震初始位置的變化量的線性關系公式為：

${\mathrm{dr}}_{kl}^i=r_k^i-r_l^i={(T_k^i-T_l^i)}^{obs}-{(T_k^i-T_l^i)}^{cal}={\Σ}_{m=1}^3(\frac{\∂T_k^i}{\∂x_m^i}-\frac{\∂T_l^i}{\∂x_m^i}){\mathrm{dx}}_m^i+s_k-s_l;$

所述雙對雙差與地震位置相對于地震初始位置的變化量的線性關系公式為：

其中，為地震事件i相對于地震初始位置在三個方向上的變化量；

將上述兩個線性關系公式以矩陣的形式表示，構成反演系統：

$\left[\begin{array}{c}WG\\ \mathrm{\λ}I\end{array}\right]m=\left[\begin{array}{c}Wd\\ 0\end{array}\right]$

其中，G為臺站對到時差與雙對到時差數據對地震位置的偏導數矩陣，所述偏導數矩陣根據預先設定的各個所述地震事件的初始位置計算得到，所述偏導數矩陣的維度是M*3N，M為所述臺站對到時差數據數量與雙對到時差數據數量的總和，N為所述地震事件數量，m是N個地震事件的位置在x,y,z三個方向上的變化量向量[dx,dy,dz]，d為各個所述臺站對雙差數據與各個所述雙對雙差數據，W為權重矩陣，λ為阻力因子，I為單位陣；

為所述臺站對雙差設置第一權重，為所述雙對雙差設置第二權重。

6.根據權利要求1所述的方法，其特征在于，所述采取分級加權的方式來反演得到各個地震事件的地震位置，包括：

根據LSQR算法求解所述反演系統，得到各個所述地震事件的地震位置，并將各個所述地震事件的地震位置分別作為各個所述地震事件的新的初始位置，開始迭代求解各個所述地震事件的地震位置，具體為：在初期迭代中，使所述臺站對雙差的第一權重大于所述雙對雙差的第二權重，直到迭代次數達到第一預設次數；在后期迭代中，使所述雙對雙差的第二權重大于所述臺站對雙差的第一權重，直到迭代次數達到第二預設次數，迭代結束，得到最終各個所述地震事件的地震位置。

7.一種雙對雙差地震定位裝置，其特征在于，所述裝置包括：

提取單元，用于基于地震的初至波絕對到時數據或者地震波形的互相關數據，按預設規則，提取觀測的地震事件對(i,j)到同一臺站的到時差數據以及同一地震到臺站對(k,l)的到時差數據，進而提取觀測的雙對到時差數據；所述雙對為所述地震事件對(i,j)和所述臺站對(k,l)的組合，其中，i和j分別表示地震事件對中的兩個地震事件，地震事件i和地震事件j的位置相鄰，l和k分別表示臺站對中的兩個臺站；

第一計算單元，用于計算所述地震事件對(i,j)到臺站k的事件對雙差和所述地震事件對(i,j)到臺站l的事件對雙差dr_l^ij；所述事件對雙差為所述地震事件對(i,j)的到時差數據的觀測值與預測值的差；

第二計算單元，用于計算所述地震事件i到臺站對(k,l)的臺站對雙差和所述地震事件j到所述臺站對(k,l)的臺站對雙差所述臺站對雙差為所述臺站對(k,l)的到時差數據的觀測值與預測值的差；

第三計算單元，用于計算與dr_l^ij的差值，或計算與的差值，得到所述地震事件對(i,j)到所述臺站對(k,l)的雙對雙差所述雙對雙差為所述雙對的到時差數據的觀測值與預測值的差；

構建單元，用于根據臺站對雙差與地震位置相對于地震初始位置的變化量的線性關系，以及雙對雙差與地震位置相對于地震初始位置的變化量的線性關系，構建反演系統；

第四計算單元，用于采取分級加權的方式來反演得到各個地震事件的地震位置。

8.根據權利要求7所述的裝置，其特征在于，所述第一計算單元包括：

第一計算子單元，用于根據所述地震事件對(i,j)中地震事件i和地震事件j到臺站k的到時和的觀測值與預測值的差，得到所述地震事件i和所述地震事件j到臺站k的到時殘差和并根據所述地震事件i和所述地震事件j到所述臺站k的到時殘差和得到所述地震事件對(i,j)到所述臺站k的事件對雙差所述地震事件對(i,j)為多個地震事件對中的任意一個地震事件對，

其中，dx為地震位置的變化量，dτ為發震時刻的變化量，為觀測的事件對(i,j)到時差數據，為預測的事件對(i,j)到時差數據；

第二計算子單元，用于根據所述地震事件i和所述地震事件j到臺站l的到時T_lⁱ和T_l^j的觀測值與預測值的差，得到所述地震事件i和所述地震事件j到臺站l的到時殘差r_lⁱ和r_l^j，并根據所述地震事件i和所述地震事件j到所述臺站l的到時殘差r_lⁱ和r_l^j，得到所述地震事件對(i,j)到所述臺站l的事件對雙差dr_l^ij，

其中，dx為地震位置的變化量，dτ為發震時刻的變化量，(T_lⁱ-T_l^j)^obs為觀測的事件對(i,j)到時差數據，(T_lⁱ-T_l^j)^cal為預測的事件對(i,j)到時差數據。

9.根據權利要求7所述的裝置，其特征在于，所述第二計算單元包括：

第三計算子單元，用于根據所述地震事件i到所述臺站k和所述臺站l的到時和T_lⁱ的觀測值與預測值的差，得到所述地震事件i到所述臺站k和所述臺站l的到時殘差和r_lⁱ，并根據所述地震事件i到所述臺站k和所述臺站l的到時殘差和r_lⁱ，得到地震i到所述臺站對(k,l)的臺站對雙差所述臺站對(k,l)為多個臺站對中的任意一個臺站對；

其中，s_k為所述臺站k的校正項，s_l為所述臺站l的校正項，為觀測的臺站對(k,l)到時差數據，為預測的臺站對(k,l)到時差數據；

第四計算子單元，用于根據所述地震事件j到所述臺站k和所述臺站l的到時和T_l^j的觀測值與預測值的差，得到所述地震事件j到所述臺站k和所述臺站l的到時殘差和r_l^j并根據所述地震事件j到所述臺站k和所述臺站l的到時殘差和r_l^j，得到地震j到所述臺站對(k,l)的臺站對雙差

其中，s_k為所述臺站k的校正項，s_l為所述臺站l的校正項，為觀測的臺站對(k,l)到時差數據，為預測的臺站對(k,l)到時差數據。

10.根據權利要求7所述的裝置，其特征在于，所述第三計算單元具體用于：計算地震事件對(i,j)到臺站k的事件對雙差與地震事件對(i,j)到臺站l的事件對雙差dr_l^ij的差值，或計算地震事件i到臺站對(k,l)的臺站對雙差與地震事件j到臺站對(k,l)的臺站對雙差的差值，得到所述地震事件對(i,j)到所述臺站對(k,l)的雙對雙差所述地震事件對(i,j)與所述臺站對(k,l)的雙對組合為多個雙對中的任意一個雙對；

$\begin{array}{c}{\mathrm{dr}}_{kl}^{ij}={\{(T_k^i-T_k^j)-(T_l^i-T_l^j)\}}^{obs}-{\{(T_k^i-T_k^j)-(T_l^i-T_l^j)\}}^{cal}\\ ={\{(T_k^i-T_l^i)-(T_k^j-T_l^j)\}}^{obs}-{\{(T_k^i-T_l^i)-(T_k^j-T_l^j)\}}^{cal}\end{array},$

其中，和為雙對到時差數據。

11.根據權利要求7所述的裝置，其特征在于，所述構建單元包括：

構建子單元，用于將臺站對雙差與地震位置相對于地震初始位置的變化量的線性關系公式，以及雙對雙差與地震位置相對于地震初始位置的變化量的線性關系公式以矩陣的形式進行表示，構成反演系統：

$\left[\begin{array}{c}WG\\ \mathrm{\λ}I\end{array}\right]m=\left[\begin{array}{c}Wd\\ 0\end{array}\right]$

其中，G為臺站對到時差與雙對到時差數據對地震位置的偏導數矩陣，所述偏導數矩陣根據預先設定的各個所述地震事件的初始位置計算得到，所述偏導數矩陣的維度是M*3N，M為所述臺站對到時差數據數量與雙對到時差數據數量的總和，N為所述地震事件數量，m是N個地震事件的位置在x,y,z三個方向上的變化量向量[dx,dy,dz]，d為各個所述臺站對雙差數據與各個所述雙對雙差數據，W為權重矩陣，λ為阻力因子，I為單位陣；

所述臺站對雙差與地震位置相對于地震初始位置的變化量的線性關系公式為：

${\mathrm{dr}}_{kl}^i=r_k^i-r_l^i={(T_k^i-T_l^i)}^{obs}-{(T_k^i-T_l^i)}^{cal}={\mathrm{\Σ}}_{m=1}^3(\frac{\∂T_k^i}{\∂x_m^i}-\frac{\∂T_l^i}{\∂x_m^i}){\mathrm{dx}}_m^i+s_k-s_l;$

所述雙對雙差與地震位置相對于地震初始位置的變化量的線性關系公式為：

其中，為地震事件i相對于地震初始位置在三個方向上的變化量；

設置子單元，用于為所述臺站對雙差設置第一權重，為所述雙對雙差設置第二權重。

12.根據權利要求7所述的裝置，其特征在于，所述第四計算單元具體用于：根據LSQR算法求解所述反演系統，得到各個所述地震事件的地震位置，并將各個所述地震事件的地震位置分別作為各個所述地震事件的新的初始位置，開始迭代求解各個所述地震事件的地震位置，具體為：在初期迭代中，使所述臺站對雙差的第一權重大于所述雙對雙差的第二權重，直到迭代次數達到第一預設次數；在后期迭代中，使所述雙對雙差的第二權重大于所述臺站對雙差的第一權重，直到迭代次數達到第二預設次數，迭代結束，得到最終各個所述地震事件的地震位置。