技術(shù)特征：

1.一種采集腳印壓制方法，其特征在于，包括：

將待處理的地震數(shù)據(jù)按照偏移距和方位角分選到由炮線和檢波線組成的十字排列域內(nèi)，得到N個(gè)共偏移距共方位角域道集，所述共偏移距共方位角域道集是指具有相同偏移距和相同方位角的地震道的集合，N為相同偏移距和相同方位角區(qū)域被覆蓋的次數(shù)，N為大于或等于1的整數(shù)；

分別將所述炮線之間的距離和所述檢波線之間的距離進(jìn)行減半處理，根據(jù)每個(gè)所述共偏移距共方位角域道集中的每個(gè)地震道的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，計(jì)算出新增加的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，每個(gè)所述共偏移距共方位角域道集生成四個(gè)子共偏移距共方位角域道集；

對(duì)每個(gè)所述子共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值，得到加密數(shù)據(jù)，所述加密數(shù)據(jù)包括原始地震數(shù)據(jù)和插值后的數(shù)據(jù)；

將所述加密數(shù)據(jù)按照偏移距和方位角重新分選到由減半后的炮線距和檢波線距組成的十字排列域內(nèi)，得到4N個(gè)子共偏移距共方位角域道集；

對(duì)所述4N個(gè)子共偏移距共方位角域道集進(jìn)行疊加前偏移處理，得到壓制采集腳印后的地震數(shù)據(jù)分析結(jié)果。

2.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述根據(jù)每個(gè)所述共偏移距共方位角域道集中的每個(gè)地震道的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，計(jì)算出新增加的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，具體包括：

通過改變所述共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)坐標(biāo)，而不改變所述共偏移距共方位角域道集中每個(gè)地震道的檢波點(diǎn)坐標(biāo)，計(jì)算出新增加的炮點(diǎn)的坐標(biāo)；

或者

通過改變所述共偏移距共方位角域道集中每個(gè)地震道的檢波點(diǎn)坐標(biāo)，而不改變所述共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)坐標(biāo)，計(jì)算出新增加的檢波點(diǎn)的坐標(biāo)；

或者

通過改變所述共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)坐標(biāo)和檢波點(diǎn)坐標(biāo)，計(jì)算出新增加的炮點(diǎn)的坐標(biāo)和檢波點(diǎn)的坐標(biāo)。

3.根據(jù)權(quán)利要求2所述的方法，其特征在于，所述通過改變所述共偏移 距共方位角域道集中每個(gè)地震道的炮點(diǎn)坐標(biāo)，而不改變所述共偏移距共方位角域道集中每個(gè)地震道的檢波點(diǎn)坐標(biāo)，計(jì)算出新增加的炮點(diǎn)的坐標(biāo)，具體包括：

采用下述坐標(biāo)變換公式(1)得到所述新增加的炮點(diǎn)的坐標(biāo)：

$\left(\begin{array}{c}x\\ y\end{array}\right)=\left(\begin{array}{cc}1& 0\\ 0& 1\end{array}\right)\left(\begin{array}{c}{x}_s\\ y_s\end{array}\right)+\left(\begin{array}{c}\frac{-{\mathrm{dl}}_s}{2}cos\mathrm{\θ}\\ \frac{{\mathrm{dl}}_s}{2}sin\mathrm{\θ}\end{array}\right)---\left(1\right)$

其中，x_s和y_s分別為所述共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)的橫坐標(biāo)和縱坐標(biāo)，x和y分別為所述新增加的炮點(diǎn)的橫坐標(biāo)和縱坐標(biāo)，θ為所述共偏移距共方位角域道集中炮線的方位角，dl_s為所述共偏移距共方位角域道集中的炮線距。

4.根據(jù)權(quán)利要求2所述的方法，其特征在于，所述通過改變所述共偏移距共方位角域道集中每個(gè)地震道的檢波點(diǎn)坐標(biāo)，而不改變所述共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)坐標(biāo)，計(jì)算出新增加的檢波點(diǎn)的坐標(biāo)，具體包括：

采用下述坐標(biāo)變換公式(2)得到所述新增加的檢波點(diǎn)的坐標(biāo)：

$\left(\begin{array}{c}{x}^{\′}\\ y^{\′}\end{array}\right)=\left(\begin{array}{cc}1& 0\\ 0& 1\end{array}\right)\left(\begin{array}{c}{x}_r\\ y_r\end{array}\right)+\left(\begin{array}{c}\frac{{\mathrm{dl}}_r}{2}cos\mathrm{\α}\\ \frac{{\mathrm{dl}}_r}{2}\sin \mathrm{\α}\end{array}\right)---\left(2\right)$

其中，x_r和y_r分別為所述共偏移距共方位角域道集中每個(gè)地震道的檢波點(diǎn)的橫坐標(biāo)和縱坐標(biāo)，x′和y′分別為所述新增加的檢波點(diǎn)的橫坐標(biāo)和縱坐標(biāo)，α為所述共偏移距共方位角域道集中檢波線的方位角，dl_r為所述共偏移距共方位角域道集中的檢波線距。

5.根據(jù)權(quán)利要求1～4任一項(xiàng)所述的方法，其特征在于，所述對(duì)每個(gè)所述子共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值，得到加密數(shù)據(jù)，具體包括：

利用五維插值法對(duì)每個(gè)所述子共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值處理，得到加密數(shù)據(jù)。

6.一種采集腳印壓制系統(tǒng)，其特征在于，包括：

獲取模塊，用于將待處理的地震數(shù)據(jù)按照偏移距和方位角分選到由炮線和檢波線組成的十字排列域內(nèi)，得到N個(gè)共偏移距共方位角域道集，所述共 偏移距共方位角域道集是指具有相同偏移距和相同方位角的地震道的集合，N為相同偏移距和相同方位角區(qū)域被覆蓋的次數(shù)，N為大于或等于1的整數(shù)；

處理模塊，用于分別將所述炮線之間的距離和所述檢波線之間的距離進(jìn)行減半處理，根據(jù)每個(gè)所述共偏移距共方位角域道集中的每個(gè)地震道的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，計(jì)算出新增加的炮點(diǎn)和/或檢波點(diǎn)的坐標(biāo)，每個(gè)所述共偏移距共方位角域道集生成四個(gè)子共偏移距共方位角域道集；

插值模塊，用于對(duì)每個(gè)所述子共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值，得到加密數(shù)據(jù)，所述加密數(shù)據(jù)包括原始地震數(shù)據(jù)和插值后的數(shù)據(jù)；

所述獲取模塊，還用于將所述加密數(shù)據(jù)按照偏移距和方位角重新分選到由減半后的炮線距和檢波線距組成的十字排列域內(nèi)，得到4N個(gè)子共偏移距共方位角域道集；

所述處理模塊，還用于對(duì)所述4N個(gè)子共偏移距共方位角域道集進(jìn)行疊加前偏移處理，得到壓制采集腳印后的地震數(shù)據(jù)分析結(jié)果。

7.根據(jù)權(quán)利要求6所述的系統(tǒng)，其特征在于，所述處理模塊，包括：第一處理單元、第二處理單元和第三處理單元；

所述第一處理單元，用于通過改變所述共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)坐標(biāo)，而不改變所述共偏移距共方位角域道集中每個(gè)地震道的檢波點(diǎn)坐標(biāo)，計(jì)算出新增加的炮點(diǎn)的坐標(biāo)；

所述第二處理單元，用于通過改變所述共偏移距共方位角域道集中每個(gè)地震道的檢波點(diǎn)坐標(biāo)，而不改變所述共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)坐標(biāo)，計(jì)算出新增加的檢波點(diǎn)的坐標(biāo)；

所述第三處理單元，用于通過改變所述共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)坐標(biāo)和檢波點(diǎn)坐標(biāo)，計(jì)算出新增加的炮點(diǎn)的坐標(biāo)和檢波點(diǎn)的坐標(biāo)。

8.根據(jù)權(quán)利要求7所述的系統(tǒng)，其特征在于，所述第一處理單元，具體用于采用下述坐標(biāo)變換公式(1)得到所述新增加的炮點(diǎn)的坐標(biāo)：

$\left(\begin{array}{c}x\\ y\end{array}\right)=\left(\begin{array}{cc}1& 0\\ 0& 1\end{array}\right)\left(\begin{array}{c}{x}_s\\ y_s\end{array}\right)+\left(\begin{array}{c}\frac{-{\mathrm{dl}}_s}{2}cos\mathrm{\θ}\\ \frac{{\mathrm{dl}}_s}{2}sin\mathrm{\θ}\end{array}\right)---\left(1\right)$

其中，x_s和y_s分別為所述共偏移距共方位角域道集中每個(gè)地震道的炮點(diǎn)的橫坐標(biāo)和縱坐標(biāo)，x和y分別為所述新增加的炮點(diǎn)的橫坐標(biāo)和縱坐標(biāo)，θ為 所述共偏移距共方位角域道集中炮線的方位角，dl_s為所述共偏移距共方位角域道集中的炮線距。

9.根據(jù)權(quán)利要求7所述的系統(tǒng)，其特征在于，所述第二處理單元，具體用于采用下述坐標(biāo)變換公式(2)得到所述新增加的檢波點(diǎn)的坐標(biāo)：

$\left(\begin{array}{c}{x}^{\′}\\ y^{\′}\end{array}\right)=\left(\begin{array}{cc}1& 0\\ 0& 1\end{array}\right)\left(\begin{array}{c}{x}_r\\ y_r\end{array}\right)+\left(\begin{array}{c}\frac{{\mathrm{dl}}_r}{2}cos\mathrm{\α}\\ \frac{{\mathrm{dl}}_r}{2}\sin \mathrm{\α}\end{array}\right)---\left(2\right)$

其中，x_r和y_r分別為所述共偏移距共方位角域道集中每個(gè)地震道的檢波點(diǎn)的橫坐標(biāo)和縱坐標(biāo)，x′和y′分別為所述新增加的檢波點(diǎn)的橫坐標(biāo)和縱坐標(biāo)，α為所述共偏移距共方位角域道集中檢波線的方位角，dl_r為所述共偏移距共方位角域道集中的檢波線距。

10.根據(jù)權(quán)利要求6～9任一項(xiàng)所述的系統(tǒng)，其特征在于，所述插值模塊，具體用于利用五維插值法對(duì)每個(gè)所述子共偏移距共方位角域道集進(jìn)行數(shù)據(jù)插值處理，得到加密數(shù)據(jù)。