地震三维观测系统共中心点面元属性的并行计算方法.pdfVIP

第39卷第5期 物 探 与 化 探 V0l_39．No．5 2015年 10月 GEOPHYSICAL ＆ GEOCHEMICALEXPLORATION Oct．，2015 doi：10．11720／wtyht．2015．5．32 李逢春，王润秋，赵薇薇，等．地震三维观测系统共中心点面元属性的并行计算方法[J]．物探与化探，2015，39(5)：1074～1078．http：／／doi．org／10 1I720／wtyht．2015．5．32 “FC，WangRQ，ZhaoWW，eta1．Commonmid-pointbinattributenaalysisofthree-dimensionalseismicgeometryusingparallelcomputationmethod[J] GeophysicalandGeochemicalExploration，2015，39(5)：1074—1078．http：／／doi．org／lO．11720／wtyht．2015．5．32 地震三维观测系统共中心点面元属性的并行计算方法 李逢春 ，王润秋 ，赵薇薇 ，杜清波 ，郭武 ，王汉钧 (1．中国石油大学(北京)地球物理与信息IY-程学院，北京 102249；2．东方地球物理公司采集技术 支持部，河北 涿州 072751) 摘 要：针对高密度地震勘探下的三维观测系统海量面元属性计算需求 ，基于OpenCL并行编程模型．研究并实现了 适用于异构环境下面元属性的并行计算方法。根据面元分析算法特点将整个网格分成多个子区域．子区域可由不 同设备以粗粒度方式并行，而子区域内部以炮检对方式细粒度并行。采用HPZ820工作站，在CPU+GPU异构混合 平台下的试算结果表明，异构多核处理速度是CPU(单核)的3O倍以上 ，数据生成速率高于300MB／s。 关键词：地震观测系统 ；面元分析；并行计算；异构计算 ；OpenCL 中图分类号：P631．4 文献标识码：A 文章编号：1000—8918(2015)05—1074—05 为了合理设计野外观测系统方案，需要对不同 个关系片计算后的统计和累加结果，CMP面元属性 观测系统的共 中心点(commonmid．point，CMP)面元 计算的工作量大、耗时长，对存储容量和I／0的要求 属．I生进行对比分析。常规 CMP面元属性分析有覆 更高，成为生产应用的瓶颈。目前，普遍应用于野外 盖次数、炮检距和炮检方位角。对覆盖次数、炮检距 施工设计 的软件主要有美 国 ION公 司的 Green 和炮检方位角的直观和精确描述，是优化观测系统 Mountain软件 和 GEDCO公司的 OMNI软件_9]。 设计的关键因素之一，而CMP面元属性计算则是实 均面临计算时间长、结果数据量太大导致计算失败 现观测系统面元属性分析的基本技术手段[1-3]。 等问题 。这些软件面对超大观测系统和高覆盖次数 近年来，随着地球物理探测技术的发展，特别是 (千次以上)面元分析任务的挑战，也显得无能为 万道以上地震仪的推广应用，高密度地震勘探技术 力。在对 CMP面元属性分析的过程中，由于传统的 正逐渐成为地震勘探技术发展的主流Ⅲ4]。高密度 面元属性计算方法是串行执行的，即依次计算每个 地震勘探通过减小面元尺度和提高空间采样率来增 面元的面元属性 ，因此 ，在多核硬件环境下 ．不能充 加采样密度 。目前，常规地震勘探面元一般为25m 分利用计算资源 ；由于串行执行CMP面元属性计算 x25131，而为了适应未来的石油勘探开发的精度要 的效率低，而不能满足高密度地震勘探带来海量 求，可能采用5m的空问采样密度，这对于山地勘 CMP面元属性计算的需求．从而导致了分析 CMP 探和小断块、薄储层、小砂体、小尺度孔洞的识别以 面元属性时，效率低、资源利用率低的问题。 及精细油藏描述具有重要意