一种地震勘探相关／叠加数据处理方法技术

本发明专利技术涉及一种地震勘探相关／叠加数据处理方法，（１）快速傅里叶变换算法与哈特利变换性质的结合混合算法；（２）ＦＦＴ－Ｈａｒｔｌｙ混合算法中分段与求和；（３）多核、多ＣＰＵ并行处理的线程安排；通过在ＩＢＭ?３７５５服务器（４核）上进行测试可以得到如下指标：在不加载噪声控制情况下，２毫秒采样率，２４秒扫描信号，３２秒地震采集，相关／叠加运算后得到８秒记录长度的情况下，处理４０００道时间不大于３秒，处理８０００道不大于４秒水平，对于加载噪声控制部分时（分段加权叠加）运算８０００道（每道记录长度与前面相同）不大于７秒。

【技术实现步骤摘要】

本专利技术涉及一种大型地震仪器野外地震勘探数据采集快速相关/叠加数据处理 方法。
技术介绍
利用可控震源作为地震采集信号激发源，进行野外地震勘探数据采集施工已经非 常普遍。对采集的地震数据进行相关/叠加运算与处理的最常使用的方法是通过硬件的手 段制作相关运算器与叠加运算器，进行实时的运算处理。以国外某种大型地震仪器为例，叠 加运算与相关运算分别由插在计算机总线上的噪声编辑处理板卡(即，NEP板)和快速傅 立叶处理板卡(即，FTP板)来完成。其中FTP板拥有2片摩托罗拉96002DSP处理器，利 用快速傅里叶变换算法(FFT)在频域执行相关运算。这两块处理板配合使用在两毫秒采样 率二十六秒采集长度的条件下，最大处理能力为1800道。使用硬件进行相关/叠加运算的优点是非常明显的，那就是实时性好、反应迅速。 例如摩托罗拉96002DSP处理器为精简指令集，具有每秒一千三百五十万次的处理能力。但使用硬件进行相关/叠加运算也存在一些难以避免的问题(1)随着地震勘探施工方法的不断发展与丰富，使得同一时刻参与数据采集的地 震道数不断扩大，同一时刻需要进行相关/叠加运算的数据海量增加，使得使用硬件相关/ 叠加运算器的数据吞吐能力出现了瓶颈。(2)由于受到硬件设计目标局限性的影响，以及选用的电子器件技术指标的影响， 这种“瓶颈”很难在后期使用阶段通过某种扩展手段进行根本性的解决。(3)相关/叠加运算器的硬件电路设计复杂，研发周期长，难以在短时间保持研发 与施工需求变化同步。(4)某些行业进行地震采集施工时(如煤田地质行业)对同一时刻进行的地震 采集道数要求比较小，如果配置处理能力过强的相关/叠加运算器硬件设备，从资源成本 的角度上讲，显得过于浪费。但为此专门开发一种相应的相关/叠加运算器的硬件设备又 得不偿失。(5)研发相关/叠加运算器硬件设备，资金投入大，风险高。随着通用计算机CPU运算能力的不断提高，硬件相关/叠加运算器有逐步转为全 软件方式实现。使用软件进行相关运算常用的方法为利用傅立叶变换性质和快速傅立叶算法，对 一组数字信号利用两次正变换和一次逆变换，进行运算求解。同时域算法相比，提高了计算 效率。由于在地震勘探数据采集过程中进行相关运算的信号均属于实序列，而傅立叶变 换是复序列，在运算过程中存在很大的冗余，限制了存储效率与运算速度的进一步提高。利 用实序列傅立叶变换的对称性进行一步简化运算方法，通过将两组实序列构造成复序列结 构进行运算，可通过一次变换得到两组实序列的傅立叶变换。这比直接采用傅里叶变换减少了一次傅立叶变换的处理过程，从而节省了近三分之一的运算量，大大缩短了运算时间。以国外某种大型地震仪器为例，使用SUN Blade 2500工作站(8G内存)在2毫秒 采样率，26秒采集长度的条件下，最大处理能力为1800道。使用软件进行相关/叠加运算的优点是(1)可扩展性好，用户可以根据施工处理需要，对运行相关/叠加运算的工作站配 置指标进行适当的调整。(2)相对于研发硬件相关/叠加运算器而言，开发软件相关/叠加运算器研发周期短。(3)利用软件手段实现通用的相关/叠加算法较为简单，因此与研发硬件相关/叠 加运算器相比风险小。(4)研发所需资金投入较小。现今使用的软件相关/叠加运算急需解决的问题(1)同硬件相关/叠加运算器相比，实时性差，处理同等数量级的数据时软件相关 /叠加运算器需要更长的时间。(2)难以适应地震勘探施工中大道数地震数据采集的发展趋势。
技术实现思路
本专利技术的目的是提供一种对FFT-Hartly混合算法的计算流程和运算步骤进行重 新的编排与设计，增加相关运算器的运算效率，提高运算速度，减少运算时间的地震勘探相 关/叠加数据处理方法。从地震采集施工数据流回收工作流程出发，结合多核、多CPU工作 站(服务器)运算与并行处理机制，设计、编排相关/叠加运算的工作处理顺序。在最大 限度上复用数据流回收工作流程中，步骤间工作状态相互转换的空白时间，予以实现相关/ 叠加运算器拥有快速处理大道数地震采集施工中地震采集数据的处理能力。相关/叠加运算器处理能力的预期目标(1)在同一时刻进行8000道数据采集，地震采集施工技术指标为，2毫秒采样率， 可控震源扫描长度为24秒，进行32秒地震数据采集，相关/叠加运算后输出8秒记录长度 时，相关/叠加运算的处理时间应在8秒内完成。(2)在同一时刻进行4000道数据采集，地震采集施工技术指标为，2毫秒采样率， 可控震源扫描长度为24秒，进行32秒地震数据采集，相关/叠加运算后输出8秒记录长度 时，相关/叠加运算的处理时间应在4秒内完成。本专利技术所述的地震勘探相关/叠加数据处理方法是由(1)快速傅里叶变换算法与哈特利变换性质的结合(FFT-Hartly混合算法)地震勘探数据采集中的相关运算在从时域到频域，频域到时域的变换过程中展 开，正变换后对频域中的数据进行如下构造 其反变换过程则有 以FFT正变换的形式实现了从频域向时域的反变换过程。(2) FFT-Hartly混合算法中的分段与求和将一组实序列分成几段，分别进行FFT运算并保存中间结果及kx.y. 0 )，下 标j为正整数，代表第几段，当这组实序列所有分段全部运算结束后，通过将保存的中间结 果进行累加求出及hxy Ck )再进行后续的运算；设需要输出的相关rxy(n)样点长度为L，扫描记录y (n)样点长度为M，采集记录 x (n)样点长度为N，将y (n)分为每段长为p的J段，则有M = JXp ；对x (n)以R为长度分 段，且满足K = 2kl彡L+p，kl代表某一满足条件的正整数，当& = 2kl > L+p时，要通过补 零的方法使得补零后的L+p段与&相等才能进行运算；如果x(n)某一分段，通常是最后一 段，由于样点数的不足，没有达到的要求，也应该通过补零的手段填满剩余的部分再进行运 算； 求出i hxy {k )可以根据FFT-Hartly混合算法中讨论过的结论求出rxy(n)；(3)多核、多CPU并行处理的线程安排将一个地震道采集数据多个分段中的一个数据分段作为一个单独的处理线程，以 多核CPU中的一个物理内核作为一个独立处理单元，由这个处理单元完成对这个数据分段的运算与处理，并以中间结果的形式存入临时缓冲区，(即，灭^。(k))。随后会有其它数据分段处理线程分配到这个处理单进行运算，周而复始直至所有 数据分段的处理线程全部计算完毕，再到缓冲区中取出中间结果，以数据分段的应对关系进行相加求和 这样每一个地震道采集数据都会求得一个及Hxy (k )，然后再将每一个地震道采集数据的及t (A;)作为一个处理线程，分配到多核cpu中的一个物理内核中进行最后 的运算，求出最后的计算结果(lf，rxy(n))。这样软件相关/叠加运算器在多核、多CPU的条件下，可以同一时刻完成对多个地 震道采集数据进行并行计算与处理。综上所述，通过对算法结构的有机调整、计算过程中最大程度的时间复用以及大 量采用多线程技术的软件工程手段，使得相关/叠加运算器的计算效率、处理能力得到了 巨大的提高。专利技术的效果通过在IBM 3755服务器(4核)上进行测试可以得到如下指标在不加载噪声控 制情况下，2本文档来自技高网...

【技术保护点】
一种地震勘探相关／叠加数据处理方法，其特征在于：是由（１）快速傅里叶变换算法与哈特利变换性质的结合混合算法：地震勘探数据采集中的相关运算在从时域到频域，频域到时域的变换过程中展开，正变换后对频域中的数据进行如下构造：Ｒ↓［ｘｙ］↑［ｈ］（ｋ）＝Ｘ↓［Ｒ］（ｋ）［Ｙ↓［Ｒ］（ｋ）＋Ｙ↓［Ｉ］（ｋ）］－Ｘ↓［Ｉ］（ｋ）［Ｙ↓［Ｒ］（ｋ）－Ｙ↓［Ｉ］（ｋ）］ｋ＝０，……，Ｎ／２Ｒ↓［ｘｙ］↑［ｈ］（Ｎ－ｋ）＝Ｘ↓［Ｒ］（ｋ）［Ｙ↓［Ｒ］（ｋ）－Ｙ↓［Ｉ］（ｋ）］＋Ｘ↓［Ｉ］（ｋ）［理线程，分配到多核ＣＰＵ中的一个物理内核中进行最后的运算，求出最后的计算结果ｒ↓［ｘｙ］（ｎ）。Ｙ↓［Ｒ］（ｋ）＋Ｙ↓［Ｉ］（ｋ）］其反变换过程则有：ｒ↓［ｘｙ］（ｎ）＝ＩＤＨＴ［Ｒ↓［ｘｙ］↑［ｈ］（ｋ）］＝１／ＮＤＨＴ［Ｒ↓［ｘｙ］↑［ｈ］（ｋ）］＝１／Ｎ｛Ｒｅ（ＤＦＴ［Ｒ↓［ｘｙ］↑［ｈ］（ｋ）］）－Ｉｍ（ＤＦＴ［Ｒ↓［ｘｙ］↑［ｈ］（ｋ）］）｝以ＦＦＴ正变换的形式实现了从频域向时域的反变换过程；（２）ＦＦＴ－Ｈａｒｔｌｙ混合算法中分段与求和将一组实序列分成几段，分别进行ＦＦＴ运算并保存中间结果Ｒ↓［ｘ↓［ｊ］ｙ↓［ｊ］］↑［ｈ］（ｋ），下标ｊ为正整数，代表第几段，当这组实序列所有分段全部运算结束后，通过将保存的中间结果进行累加求出Ｒ↓［ｘｙ］↑［ｈ］（ｋ）再进行后续的运算；设需要输出的相关ｒ↓［ｘｙ］（ｎ）样点长度为Ｌ，扫描记录ｙ（ｎ）样点长度为Ｍ，采集记录ｘ（ｎ）样点长度为Ｎ，将ｙ（ｎ）分为每段长为ｐ的Ｊ段，则有Ｍ＝Ｊ×ｐ；对ｘ（ｎ）以Ｎ１为长度分段，且满足Ｎ↓［１］＝２↑［ｋ１］≥Ｌ＋ｐ，ｋ１代表某一满足条件的正整数，当Ｎ↓［１］＝２↑［ｋ１］＞Ｌ＋ｐ时，要通过补零的方法使得补零后的Ｌ＋ｐ段与Ｎ↓［１］相等才能进行运算；如果ｘ（ｎ）某一分段，通常是最后一段，由于样点数的不足，没有达到的要求，也应该通过补零的手段填满剩余的部分再进行运算；Ｒ↓［ｘｙ］↑［ｈ］（ｋ）＝＊Ｒ↓［ｘ↓［ｊ］ｙ↓［ｊ］］↑［ｈ］（ｋ）０≤ｋ≤Ｎ－１０≤ｊ≤Ｊ－１求出Ｒ↓［ｘｙ］↑［ｈ］（ｋ）可以根据ＦＦＴ－Ｈａｒｔｌｙ混合算法中讨论过的结论求出ｒ↓［ｘｙ］（ｎ）；（３）多核、多ＣＰＵ并行处理的线程安排将一个地震道采集数据多个分段中的一个数据分段作为一个单独的处理线程，以多核ＣＰＵ中的一个物理内核作为一个独立处理单元，由这个处理单元完成对这个数据分段的运算与...

【技术特征摘要】

【专利技术属性】
技术研发人员：刘益成，穆群英，沈孝科，赵培根，魏启，李永权，王浩，张红震，张洁，
申请(专利权)人：西安英诺瓦物探装备有限公司，
类型：发明
国别省市：87[中国|西安]