[发明专利]DAS海量数据频率波数域数据压缩存储方法

权利要求书

1.DAS海量数据频率波数域数据压缩存储方法，其特征在于，包括以下步骤：根据DAS采集数据的特点，将空间时间域数据，通过快速傅里叶变换为空间频率域、设计滤波器高频截断切除、快速傅里叶逆变换为空间时间域、数据转置为时间空间域、快速傅里叶变换为时间波数域、设计滤波器高波数截断切除、快速傅里叶逆变换为时间空间域、数据转置为空间时间域；具体包括以下步骤：

S1、输入分布式光纤声波传感即DAS数据InDas_M×N，输入频率压缩参数f_com频率压缩斜坡参数f_comB、波数压缩参数k_com、波数压缩斜坡参数k_comB；

其中，InDas_M×N是空间时间域DAS数据，M是时间采样点数，N是空间采样点数；

S2、采用快速傅里叶变换将步骤S1的空间时间域数据InDas_M×N变换为空间频率域数据FDas_M×N；

FDas_M×N＝FFT(InDas_M×N)

其中，InDas_M×N是步骤S1的空间时间域数据，FFT是快速傅里叶变换函数，FDas_M×N是空间频率域数据，M是时间采样点数和频率采样点数，N是空间采样点数；

S3、用步骤S1的频率压缩参数f_com、频率压缩斜坡参数f_comB设计滤波器，高频截断切除步骤S2的空间频率域数据FDas_M×N，得到空间频率域数据FDas_MC×N；

其中，i∈[1,M]，M是时间采样点数和频率采样点数，N是空间采样点数，f_i是第i个频率，f_com是频率压缩参数，f_comB是频率压缩斜坡参数，empty是清空函数，MC是高频截断后的时间采样点数和频率采样点数，FDas_MC×N是高频截断切除的空间频率域数据；

S4、采用快速傅里叶逆变换将步骤S3的空间频率域数据FDas_MC×N变换为空间时间域IFDas_MC×N；

IFDas_MC×N＝IFFT(FDas_MC×N)

其中，FDas_MC×N是空间频率域数据，MC是高频截断后的时间采样点数和频率采样点数，N是空间采样点数，IFFT是快速傅里叶逆变换函数，IFDas_MC×N是空间时间域数据；

S5、步骤S4的空间时间域数据IFDas_MC×N转置为时间空间域数据TIFDas_N×MC；

TIFDas_N×MC＝Trans(IFDas_MC×N)

其中，IFDas_MC×N是空间时间域数据，MC是高频截断后的时间采样点数和频率采样点数，N是空间采样点数，Trans是转置函数，TIFDas_N×MC是时间空间域数据；

S6、采用快速傅里叶变换将步骤S5的时间空间域数据TIFDas_N×MC变换为时间波数域数据KTIFDas_N×MC；

KTIFDas_N×MC＝FFT(TIFDas_N×MC)

其中，TIFDas_N×MC是时间空间域数据，MC是高频截断后的时间采样点数和频率采样点数，N是空间采样点数和波数采样点数，FFT是快速傅里叶变换函数，KTIFDas_N×MC是时间波数域数据；

S7、用步骤S1的波数压缩参数k_com、波数压缩斜坡参数k_comB设计滤波器，高波数截断切除步骤S6的时间波数域数据KTIFDas_N×MC，得到时间波数域数据KTIFDas_NC×MC；

其中，j∈[1,N]，NC是高波数截断后的空间采样点数和波数采样点数，k_j是第j个频率，k_com是波数压缩参数，k_comB是波数压缩斜坡参数，empty是清空函数，MC是高频截断后的时间采样点数和频率采样点数，KTIFDas_NC×MC是高波数截断切除的时间波数域数据；

S8、采用快速傅里叶逆变换将步骤S7的时间波数域数据KTIFDas_NC×MC变换为时间空间域IKTIFDas_NC×MC；

IKTIFDas_NC×MC＝IFFT(KTIFDas_NC×MC)

其中，KTIFDas_NC×MC是时间波数域数据，MC是高频截断后的时间采样点数和频率采样点数，NC是高波数截断后的空间采样点数和波数采样点数，IFFT是快速傅里叶逆变换函数，IKTIFDas_NC×MC是时间空间域数据；

S9、步骤S8的时间空间域数据IKTIFDas_NC×MC转置为空间时间域TIKTIFDas_MC×NC，即完成了DAS海量数据频率波数域数据压缩；

TIKTIFDas_MC×NC＝Trans(IKTIFDas_NC×MC)

其中，IKTIFDas_NC×MC是时间空间域数据，MC是高频截断后的时间采样点数和频率采样点数，NC是高波数截断后的空间采样点数和波数采样点数，Trans是转置函数，TIKTIFDas_MC×NC是空间时间域数据。