曙光4000A上三维叠前深度偏移并行计算的应用设计

万方数据万方数据万方数据・208・物探化探计算技术28卷图2曙光4000A的存储系统Fig.2ThestoragesystemofDawning4000A围3曙光4000A4一wayOpteron访存机制Fig.34-wayopteroncallmemorymodeofDawnin94000A(1公平共享式调度(FairshareScheduling可以规定用户或用户组对计算资源的使用分额,从而实现针对不同用户实行不同的管理策略,以保证计算资源能被公平合理的使用。(2抢占式调度(PreemptionScheduling保证高优先级的作业能够抢占低优先级作业的资源,优先获得执行权限,从而保证资源使用的优先级。(3预约式调度(ReservationScheduling是根据整个机器的负载情况,然后由用户预约某个资源指标,当系统指标满足用户需求的资源时允许此作业运行。在进行炮域波动方程叠前深度偏移计算时,一般采用预约式调度和抢占式调度的方式提交作业。提交炮域波动方程叠前深度偏移作业时,采用预约式调度方式,指定必须将4个CPU全部分配给作业运行,调度级别为normal。虽然一个结点少于四个CPU并不影响作业的提交,但在此结点上,如有其它作业运行,就可能出现内存分配失衡,从而影响作业运行效率。接收程序与炮域波动方程叠前深度偏移运行结点绑定,采用抢占式调度方式,将优先级提高到VIP级别提交运行,收集深度偏移运行的中间结果。4三维地震数据波动方程叠前深度偏移计算某三维工区属高陡构造褶皱体系,存在不同方向的构造叠加与干扰,并造成深、浅构造的分异,形成多层次的构造变形样式。数据采集是中间放炮观测系统,12线48炮,接收道数为12木216=2592道,炮点距50m,炮线距450m,道距50m,偏移距25m,采样间隔2ms,记录长度6s,每道3001个样点,单炮记录见下页图4。实际成像面积为451.10km2,计算4252炮,三维地震资料数据约200GB。由于此三维数据是按四束排列采集,因此计算时根据束线数据分别处理。叠前深度偏移延拓总深度为12000m,延拓步长10m,纵、横向扩边均为30道。|||II万方数据万方数据・２１０・物探化探计算技术２８卷图４单炮记录Ｆｉｇ．４Ｓｅｉｓｍｉｃｓｈｏｔｒｅｃｏｒｄ万方数据图５各结点计算一炮后收集的数据Ｆｉｇ～Ｔｈｅｍｉｇｒａｔｉｏｎｓｅｃｔｉｏｎｆｒ。ｏｍａ图６一束线的计算结果Ｆｉｇ．６Ｔｈｅａｍｉｇｒａｔｉｏｎｓｅｃｔｉｏｎｆｒｏｍｓｈｏｔｅａｃｈｃｏｍｐｕｔｅｒｎｏｄｅｓｔｒｅａｍｓｅｉｓｍｉｃｄａｔａ参考文献：¨１Ｊ［Ｊ］．勘探地球物理进展，２００２．［７］孔祥宁，赵改善，王华忠，等．叠前偏移计算中多文件并行Ｉ／Ｏ技术［Ｊ］．物探化探计算技术，２００５，２７（２）：９７．李建江．一种基于共享存储的叠前深度偏移并行算法［Ｊ］．软件学报，２００２．［８］罗省贤．基于ＭＰＩ的快速富里叶变换并行算法［Ｊ］，物探化探计算技术，２０００，２２（２）：１５４．［９］心１ｊ赵改善，孔祥宁，王于静，等．“位集群计算平台波动方程叠前深度偏移的性能优化［Ｊ］．勘探地球物理进展，２００５．程玖兵，于富文．波动方程共炮检距道集叠前深度偏移［Ｊ］．石油地球物理勘探，２００１．Ｂ１ｊＨ１Ｊ孙建国，方伍宝，赵改善，等．复杂条件下深度成像方法技术研究［ｃ］．中石化科研项目成果报告，２００３．上海超级计算中心．高性能计算发展与应用［ｃ］，曙光４０００Ａ超级计算机专辑，２００５．赵改善．叠前深度偏移成像技术［Ｊ］．石油物探译丛，１９９８．庞世明．高性能计算技术及其在油气勘探中的应用［１０］潘宏勋，方伍宝．波动方程叠前深度偏移中的并行策略（油气