买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：中国石油化工集团有限公司;中石化石油工程技术服务有限公司;中石化石油工程地球物理有限公司;中石化石油工程地球物理有限公司胜利分公司

摘要：本发明提供一种海量地震源数据SEGD分段分类混合压缩方法，包括：解编SEGD数据文件，对获取的质控相关参数做霍夫曼编码；记录辅助道特定时窗内起跳时间和幅值信息；依次求取排列上每个地震道的时差校正因子，利用时差校正因子对该地震道采样做时差校正，计算该排列与上一个排列时差校正后采样的变异系数；如果变异系数大于阈值，对采样数据做整型小波变换，利用动态阈值对小波系数量化，并将量化后的数据做算术编码；若变异系数不大于阈值，则对两排列数据差值做无损编码；将以上分段、分类处理的数据按一定规则进行组码，形成字符流。该方法满足海量地震采集现场质控时效要求，又不因数据压缩造成失真。

主权项：1.海量地震源数据SEGD分段分类混合压缩方法，其特征在于，该海量地震源数据SEGD分段分类混合压缩方法包括：步骤1：解编SEGD数据文件，从各类头块和道头中获取所有质控相关的参数，对其做霍夫曼编码；步骤2：对于辅助道数据，记录特定时窗内的起跳时间、幅值信息；步骤3：根据工区勘探目标要求，按时空范围剔除排列中与质控无关的采样数据；步骤4：通过求取排列中地震道的时差校正因子，对地震道做时差校正，在该排列所有道时差校正后，计算该排列与上个排列所有采样的变异系数；步骤5：如果该变异系数大于设定的阈值，则对该排列采样做整型归一化处理；步骤6：采用Daubechies9小波基函数对归一化数据进行整数小波变换；步骤7：由小波系数计算动态阈值，采用多级树集合分裂算法量化小波系数，精度控制由受控因子限定；步骤8：对量化后的数据做算术编码，并统计编码摘要信息；步骤9：如果步骤4求取的变异系数不大于设定的阈值，则对该排列与上一排列采样的差值做差分编码并归一化；步骤10：对归一化后的数据再做霍夫曼编码，并统计编码摘要信息；步骤11：将以上分段、分类处理的数据按一定规则进行组码，形成字符流；在步骤1中，先解编总头块，计算通道组位置及数据道位置，识别辅助道和数据道，然后解编各类头块和道头，剔除标记道头字与地震道索引信息以及其它与现场资料分析及质控无关的参数，然后，对质控参数做霍夫曼编码；在步骤2中，扫描辅助道数据，记录特定时窗内的参考信号和验证信号的典型数据，包括起跳时间、幅值，直接用于组码；在步骤3中，时空范围的数据剔除是基于监控目标确定，剔除的内容是远偏排列或地震道以及深层地震数据；在步骤4中，时差校正因子用于该排列地震道数据的时差校正，其计算公式为： 式中：Δts,i为当前排列s中第i道的时差校正因子，单位为s，px、py为炮点横、纵坐标，单位为m，ps,ix、ps.iy为当前排列s中第i道的横、纵坐标，单位为m，v为低降速带速度，单位为ms；变异系数反映两个相邻地震道采样的相似程度，其计算公式为： 式中：cv为变异系数，M为该排列地震道数，σ为标准偏差，定义如下： 式中：T为该排列数据剔除后地震道的最大采样数，xs,itj、xs-1,itj分别为当前排列s和上一个排列s-1中第i道第j个离散点采样值；在步骤5中，阈值的设定是与震源类型、地表与地下地质条件密切相关的，由施工前的单炮试验资料分析确定；所述的归一化处理计算公式为： 式中：zs,itj为当前排列s中第i道第j个采样的归一化后数据；int[·]为求整函数；k为压缩位，它是不大于32的正整型值，一般取8、16或24；xs,itj为当前排列s中第i道第j个离散点采样值；minxs,i、maxxs,i分别为当前排列s中第i道的最小采样值和最大采样值；在步骤6中，整数小波变换采用Daubechies9小波基函数，分解层数为3-5层；在步骤7中，多级树集合分裂算法的动态阈值是由小波系数计算得到的，该动态阈值的计算公式为： 式中：Tp为动态阈值，p为循环次数，p＝0,1,2,…,k-1；为向下取整函数；lg·为以10为底的对数函数；max·为最大值函数；Wfa,b表示小波系数，a为时间尺度，b为频率尺度；在步骤7中，受控因子是由压缩位、压缩时间及压缩前后数据的均方误差求得的，其计算公式如下： 式中：CF为受控因子，单位：db；ξ为压缩时间因子，用于调控压缩时间与失真度，取0ξ100，k仍为压缩位，MSE为均方误差，定义如下： 式中：M为当前排列数据剔除后的地震道数，T为当前排列数据剔除后地震道的最大采样数，zs,itj、分别为当前排列s第i道第j个采样点压缩前后的数据；在步骤8中，统计的摘要信息内容包括数据剔除的时空范围M与T、压缩位k、压缩时间因子ξ、归一化因子、地震道采样极值minxs,i和maxxs,i；在步骤9中，对该排列与上一排列采样的差值做差分编码并整型归一化，归一化计算公式为： 式中：zs,itj为当前排列s中第i道第j个采样的归一化后数据；int[·]为求整函数；k为压缩位；为当前排列s中第i道第j个采样值与前一个排列相应采样的差值；为当前排列s中第i道差值的最小值和最大值；在步骤10中，对归一化后的数据做霍夫曼编码，并统计编码摘要信息，内容包括阈值、压缩位、采样差值的最小值和最大值在步骤11中，组码是按照步骤1、步骤2、步骤8、步骤10中数据的顺序组织，其中，对每一类信息制定相应的序列化规则。

全文数据：

权利要求：

百度查询： 中国石油化工集团有限公司 中石化石油工程技术服务有限公司 中石化石油工程地球物理有限公司 中石化石油工程地球物理有限公司胜利分公司 海量地震源数据SEG D分段分类混合压缩方法