量，为垂直于滑坡面的分量为了减少互相关计算量，我们首先对台站间的个分量进行两两互相关计算得到了个不同分量的互相关函数，然后类比等的方法，采用以下公式将他们旋转到滑坡面上的分关函数处理过程中的滤波处理方式互相关计算时的滤波处理总的来讲可以分为两个部分去除原始波形的仪器响应带通滤波处理其中去仪器响应是用来压制极低和极高频信号以及消除仪器的不规则响应此前研究中的带通滤波处理通常则是计算互相关之前首先滤波到个比较宽的频带，互相关计算之后再滤波到研究所需的多个比较窄的频带，这样做虽然比起多次滤波后再进行互相关计算能够减少互相关计算量，但是对于频带内较高频的信号的信噪比却有定的影响如图蓝色和红色波形分别为台站和先进行的带通滤波后进行互相关计算再对互相关函数滤波到以及直接窄带滤波到进行互相关计算得到的互相关结果，先带通滤波再窄带滤波对应的信噪比小于直接窄带滤波的信噪比，的影响我们的数据处理流程与文章提到的数据处理方法几乎致，包括对各个台站每天的数据进行重采样去仪器响应去均值去趋势以及带通滤波处理，随后进行时间域和频率域的归化处理，最后通过计算波形间的互相关并且叠加得到互相关函数而为了在滑坡上得到可靠的互相关函数，首先我们需要选取合适的数据处理方式首先是互相关函数计算所需数据的裁剪时间长度的影响，不同的裁剪时间长度决定了计算的互相关函数不同的叠加次数以及运算时间由于密集台站的互相关计算量较大，选择合适的裁剪时间长度可以有效地减少互相关计算时间并且提高互相关函数的信噪比此前的些背景噪声成像工作采用的裁剪时间长度为天或者小时，这里我们测试了裁剪时间长度对于信噪比以及运算时间的影响针对时间长度相同的地震背景噪声数据，图中的红色以及蓝色折线探讨如何应用地震背景噪声进行西山村滑坡高频面波信号的提取地球物理学论文沉积层的厚度约为，其横波速度约为，以上的信号很难约束滑坡沉积层与基岩的分界面，因此，本文主要研究频率以下的地震背景噪声信号图中的蓝色波形展示了原始波形低通滤波后的结果，可以发现大部分的事件信号被压制了图为原始波形对应的频谱图，我们能清楚地看到频率在以及两个能量较大的信号我们从原始波形中截取了持续时间比较平静的波形图，并且对该段信号进行了傅里叶变换到频率域图，发现在该段信号的主频为，这说明滑坡体上的地震背景噪声主要集中在图西山村滑坡地理位臵地形以及台站分布图为理县滑坡的地理位臵，中红色为宽频带地震仪，蓝色为短周期地震仪，红色轮廓展示了西山村滑坡的活动区域展示了滑坡的整体坡非均匀性大以及高频面波的快速衰减，从互相关函数中提取高频面波存在定难度，因此在滑坡上的噪声成像工作较少在本文中，我们对川理县西山村滑坡上布设的个流动地震台站记录的分量地震噪声进行互相关计算，提取了滑坡上的高频面波经验格林函数同时我们还测试了不同旋转方式对互相关函数计算的影响，并通过聚束分析方法获得了滑坡上的噪声源的位臵，我们得到的高频面波经验格林函数为滑坡上波和波频散曲线的提取以及横波速度反演提供了基础数据与预处理西山村滑坡位于地震活动性强的龙门山脉断裂带邓起东等，离汶川地震震中约为该滑坡是个大型堆积层滑坡，底部和顶部的海拔高度差约为，东西向跨度约野外地质调查发现滑坡体的厚度在之间，其体积超过了千万方盛敏汉等，滑坡大体朝南北向延伸，滑动方向约为，滑坡倾角约如原位地球物理测量电阻率层析成像以及地震层析成像等其中，原位地球物理测量包括压力计倾斜仪实验室测量等，可准确提供滑坡体的物理参数信息，这些原位观测可以直接确定滑坡的基底界面和边坡稳定性分析所需的参数但是，滑坡体的不稳定性对于钻孔的限制较大，同时降低了测量点的横向分辨率，方法通过测量风化和蚀变的滑坡体介质的电阻率差异包括颗粒的矿物性质地下水含量电解质的性质孔隙度等，从而得到滑坡体的维或者维电阻率结构，然而由于长时间连续观测的难度较大，对滑坡结构的时间变化的研究较少使用波和波的主动震源地震成像通过炸药和布设密集的地震检波器来获得滑坡体高分辨率地震波速度结构以描述滑坡内部特征，主动震源方法在不稳定滑坡上实施炸摘要滑坡体维速度结构为滑坡治理灾害防治风险防范以及理解滑坡的动力学过程提供了关键信息，而高频面波成像是研究浅层速度结构的重要手段本文详细介绍了利用年月日在川理县西山村滑坡上布设的个仪器记录的分量连续地震噪声数据，提取的基阶波和波经验格林函数，分析了不同参数和处理方法对背景噪声互相关计算结果的影响结果表明，天左右的连续波形记录裁剪至的时间长度，进行互相关叠加就可以得到较为稳定的经验格林函数西山村滑坡体上波和波的群速度分别为约和，并且波信噪比高于波此外，我们还利用聚束分析方法对噪声源的位臵进行了分析，发现的背景噪声主要来自滑坡东南侧附近杂谷脑河水的搬运作用这些高频面波数据和噪声源位臵为获取滑坡浅层维速度结并且通过旋转得到了滑坡上个不同分量以及的互相关函数在个分量中，包含波信息的分量和包含波信号的分量的信噪比较高，其中波视速度约为，而波的视速度约为，造成波速度大于波的原因可能是滑坡浅层较为复杂的速度结构另外，滑坡上互相关函数的收敛测试表明，天的叠加时间就能够得到较稳定的互相关函数，该短时间收敛的特点有利于对滑坡体结构随时间的变化进行研究此外，由于的互相关函数存在定噪声源方向性的影响，我们首先根据聚束分析得到了在该频带下约的能量对应的噪声主要来自于方位角约为的方向同时我们对杂谷脑河产生的噪声频谱进行了理论计算，发现河流噪声与地震仪记录到的持续噪互相关函数中，波视速度主要集中在，然而部分互相关函数的速度较大，甚至达到了为了研究波速度存在较大差异的原因，我们选取了两个方位不同的测线与，其中位于滑坡体内部，而位于滑坡体边缘图与两个测线各自的台站间互相关函数速度比较致，其中测线均为约，测线约，这可能是由于测线内的台站对与噪声源的连线对应的较小而测线内的台站对与噪声源的连线对应的较大这种速度随的变化规律与图所示的噪声源对应的传播路径对应较好，同时也证实了对互相关函数进行速度校正的重要性信号提取频带分析不同频率的面波信号对应不同的敏感深度，因此信号的提取频带也影响着速度结构反演的深度通过噪声互相关提取面波信号的频带不仅与地震仪器的响应相关，而且也受到台站的分布与噪声源信息的影响此前我们已主要噪声来源来自滑坡体的南部此外，该噪声源能量密集，说明该噪声源可能是个比较集中的源由于高频面波衰减较快而且交通噪声主要在以上，因此滑坡上的噪声可能是来自于距离滑坡体南部约为的杂谷脑河如图图分量信噪比大于的互相关函数分量信噪比大于的互相关函数此前的研究表明，对于河流噪声的频谱与沉积层的值相关，这里我们基于等文章中的方程对杂谷脑河产生的噪声频谱进行了计算，结果表明河流噪声的频谱与观察到的真实噪声数据的频谱图比较致虽然在滑坡体上的背景噪声存在许多其他来源例如人类动物活动等，但是主要的噪声能量很有可能来自于杂谷脑河在以后的工作中，我们将利用该噪声源信息对频散曲线进行校正以提高速度成像的准确度讨论噪声源对面波信号的影响对于区域存在的持续噪声，台站进行约束能够在定程度上提高反演的准确性，图，分别为和分量的剖面图，蓝色虚线对应着不同的群速到时，红色框表示了在相同位臵分量的波形比分量较差的与分量互相关函数对比图最后，我们从个台站间共条互相关函数中提取了条信噪比高于的波格林函数，以及条信噪比高于的波格林函数其中，波的视速度约为，而波的视速度约为并且，互相关函数的信噪比随着台间距的增加而降低，如图中，大于的台站对互相关函数中仅有条信噪比较高噪声源分析通常情况下，均匀分布的噪声源的互相关函数正负半轴能量对称而正负半轴能量不对称的互相关函数可能是由于背景噪声分布不均图中互相关函数负半轴的能量表示从滑坡底部传向滑坡顶部，正半轴则相反探讨如何应用地震背景噪声进行西山村滑坡高频面波信号的提取地球物理学论文声的频谱比较接近，说明滑坡上噪声很有可能来自于滑坡体南部的杂谷脑河本文的工作证实了通过噪声互相关可以得到滑坡上的高频波和波信号，通过合理的数据预处理旋转裁剪直接滤波等能够较好地提高互相关函数的信噪比，噪声源分析得到的噪声源信息确保了后续在滑坡上进行的噪声成像反演工作的准确性图不同测线的互相关函数对比图分别为中的以及测线对应的互相关函数剖面图，蓝色虚线表示不同的群速度到时，红色角星表示噪声源定位结果图不同频带分量互相关函数曾求，储日升，曾祥方，崇加军利用地震背景噪声提取西山村滑坡高频面波信号地球物理学报，基金国家重点研发计划国家自然科学基金国际地区合作与交流项目资助探讨如何应用地震背景噪声进行西山村滑坡高频面波信号的提取地球物理学论文大幅度的增加，推测可能是由于面波的频散特征，不过这点还有待在后续的工作中证实的信号信噪比较低，但与的信号速度较为接近，可能是由于在此高频区间频带较窄，导致面波速度变化并不会太大虽然相比于其他两个频带的互相关函数信号信噪比较低，但是这也说明部分台站对能够提取到更宽频带的信号图噪声源对互相关函数提取的影响本文中重新定义的台站间相对于噪声源的方位角分别为信噪比大于的互相关函数数量与正常方位角以及与重新定义的方位角直方图这里我们讨论了提取信号的可行性，后续的成像工作我们会考虑在更宽的频带中开展，以获得更多深度上的约束不过由于低于以及高于互相关函数的信号信噪比较低，在对信号进行处理的时候需要尤为注意结论在本文中，我们对布设在西山村滑坡上的个地震台的个分量进行了两两互相关计算，此