1.一种在被动源地震数据中自动提取体波的处理方法，其特征在于，包括：
对原始数据进行数据截取；
对截取到的数据进行数据选择；
对选择的数据进行时间域强振幅控制；
对进行了时间域强振幅控制后的数据再进行频率域强振幅控制；以及
数据输出。
2.如权利要求1所述的在被动源地震数据中自动提取体波的处理方法，其特征在于，对原始数据进行数据截取包括：预设第一阈值，并将所述原始数据中的高振幅记录降低到所述第一阈值以下。
3.如权利要求1所述的在被动源地震数据中自动提取体波的处理方法，其特征在于，对截取到的数据进行数据选择包括：
预设短时间窗口，并计算整个时间记录上每个所述短时间窗口内的均方根振幅值；
根据整个时间记录上均方根振幅值能量的分布，选取一个全局阈值；以及
将超过全局阈值的所对应的短时间窗口内的数据置为零，得到第一数据。
4.如权利要求1所述的在被动源地震数据中自动提取体波的处理方法，其特征在于，所述短时间窗口为1秒，且所述全局阈值为最大均方根振幅值的50％。
5.如权利要求1所述的在被动源地震数据中自动提取体波的处理方法，其特征在于，对选择的数据进行时间域强振幅控制包括：对所述第一数据使用基于l1迭代加权的最小二乘方法进行计算得到第二数据，公式(1)为，
0≈w(h-d)
其中，d为原始数据，h为选择的数据，w为对角矩阵，
6.如权利要求1所述的在被动源地震数据中自动提取体波的处理方法，其特征在于，对进行了时间域强振幅控制后的数据再进行频率域强振幅控制包括：
对所述第二数据进行傅立叶变换，并分别用公式(2)和公式(3)计算信号的谱幅值m和相位
其中，f为傅立叶变换算子，r为实部，i为虚部，arctan为反正切函数；
定义第二阈值，并应用滑动窗口对窗口内数据的谱幅值m进行筛选；
将高于第二阀值的谱幅值m降低到第二阀值，得到处理后的谱幅值md；
根据公式(4)使用谱幅值md重建信号，公式(4)为，
其中，f-1为傅立叶反变换算子，dr为重建的信号，r为实部。