[发明专利]一种地震纵、横波波场分离与去噪的方法

摘要

本发明提供一种地震纵、横波波场分离与去噪的方法，包括：根据含M0道，每道N0个采样点的地震数据生成二维数组X1；将X1中的地震数据进行正常时差校正，生成校正数据L1；对L1进行奇异值分解，生成初始奇异值矩阵∑1；根据∑1的奇异值谱获取∑1的畸变位置；将畸变位置后的奇异值构成高通滤波奇异值矩阵∑3；奇异值分解∑3，生成高通滤波重构信号LHp；将LHp进行正常时差反校正处理，生成高通反校正数据X3；根据M0、N0对X3进行扩边处理，生成扩边矩阵X4；对X4进行二维傅里叶正变换得到f-k谱Y1；根据Y1提取目标信号的f-k谱，生成目标f-k谱Y2；对Y2进行二维傅里叶反变换，生成反变换数据X5；对X5进行缩边处理，生成缩边后数据X6；将X6转换为与采样点的地震数据相同的格式并输出。

主权项

一种地震纵、横波波场分离与去噪的方法，其特征在于，所述方法包括：根据含有M0道，每道有N0个采样点的地震数据生成二维数组X1；依据纵波叠加速度和P‑P反射波时距规律，将所述二维数组X1中的地震数据进行正常时差校正，生成校正数据L1；对所述校正数据L1进行奇异值分解，生成初始奇异值矩阵∑1；根据所述初始奇异值矩阵∑1的奇异值谱获取所述初始奇异值矩阵∑1的畸变位置；将所述畸变位置前的奇异值构成低通滤波奇异值矩阵∑2，将所述畸变位置后的奇异值构成高通滤波奇异值矩阵∑3；奇异值分解所述高通滤波奇异值矩阵∑3，生成高通滤波重构信号LHp；依据纵波叠加速度和P‑P反射波时距规律，将所述高通滤波重构信号LHp进行正常时差反校正处理，生成高通反校正数据X3；根据所述M0、N0的数值判断是否需要对所述高通反校正数据X3进行扩边处理，如果是，对所述高通反校正数据X3进行扩边处理，生成扩边矩阵X4；对所述扩边矩阵X4进行二维傅里叶正变换得到f‑k谱Y1；根据所述f‑k谱Y1提取目标信号的f‑k谱，生成目标f‑k谱Y2；对所述目标f‑k谱Y2进行二维傅里叶反变换，生成反变换数据X5；对所述反变换数据X5进行缩边处理，生成缩边后数据X6；将所述缩边后数据X6转换为与所述采样点的地震数据相同的格式并输出。