地球物理学进展  2014, Vol. 29 Issue (6): 2738-2744   PDF    
地震资料处理中的整体化研究思路分析
邱铁成1, 庄祖垠1, 徐强1, 温华华2, 俞小林1    
1. 中海油服物探事业部, 天津 300451;
2. 中海石油(中国)有限公司深圳分公司, 广州 510240
摘要:通过对大量地震资料处理过程的研究,结合实例分析,总结了整体化研究思路对地震资料处理的重要性,并从四个方面对整体化研究思路进行了分析,即处理项目的整体计划,技术、经验、理论的整体匹配,不同专业人才的整体合作,处理软硬件资源的整体协调.最终认为,可以将地震资料处理看成是一个整体系统,在处理工作中灵活运用整体化研究思路可以提高工作效率并有效解决物探难题.
关键词地震资料处理     研究思路     整体化     处理解释一体化     系统    
The analysis on integrity research method of seismic data processing
QIU Tie-cheng1, ZHUANG Zu-yin1, XU Qiang1, WEN Hua-hua2, YU Xiao-lin1    
1. Geophysical-China Oilfield Services Limited, TianJin 300451, China;
2. CNOOC China Limited-Shenzhen, Guangzhou 510240, China
Abstract: Through the investigation on the previous study of seismic data processing,combined with case analysis, emphasized the importance on the integrity research method of seismic data processing, analyzed the integrity research ideas from four aspects, namely the entire project of seismic data processing, the overall matching about technique, experience and theory, the overall cooperation between various professionals, the overall coordination between processing software and hardware resources. Finally, seismic data processing can be seen as a system, in the seismic data processing work, by flexibly using integrity research method, we can improve work efficiency and quickly solve geophysical difficult problems.
Key words: seismic data processing     research method     integrity     integration of processing and interpretation     system    
0 引 言

随着油气勘探活动的不断深入,物探、地质等多学科也在逐步交叉、融合发展,这些学科不仅是人们实践经验的总结,很多都已上升为勘探理论.一位有名的西方石油地质学家说过,石油归根结底是在人们头脑中找到的,这说明了勘探工作者的研究思路影响甚至决定了找油的方向.

在地震资料处理方面,研究思路同样非常重要,近些年,前人对各种处理方法及技术均进行了深入研究,有很多研究的思路很新颖.例如,在叠前深度偏移方面,段鹏飞等(2013)针对典型的横向各向同性(TI)介质,根据深度域构造成像与偏移速度分析的需要,研究了基于射线理论的局部角度域叠前深度偏移成像方法;吴帮玉等(2013)将局部余弦基小波束波场分解、传播与观测系统沉降法叠前深度偏移相结合,推导了源-检波器观测系统沉降法传播算子;杜启振等(2011)为提高地震成像结果的准确性并真实反映实际地震波场在介质中的传播特性,提出了直接利用多分量地震数据进行弹性波角度域偏移速度分析的方法;刘礼农和张剑锋(2011)将各向同性复杂构造介质中的最优分裂Fourier方法发展为各向异性介质中的波场延拓算法,并就其工业应用中遇到的提高计算效率及反假频成像策略问题进行了讨论.在去噪方面,袁艳华等(2013)分析了curvelet变换在地震数据处理中的优势及不足,提出了改进的非二次幂curvelet变换,提升了curvelet变换对地震数据有效频带、方向和低频信息的精细处理能力;刘洋等(2011)提出了一种新的加权中值滤波方法,属于数据驱动方法,其主要特点是更加适合地震数据的非稳态特点.

以上实例对地震勘探数据处理的研究思路有很好的启示意义.作为新技术的叠前深度偏移方面,有很多前沿问题需要解决,研究思路较宽泛而不受常规思路的制约,出现了很多有益的探索,可以不断开拓出新思路;作为常规技术的去噪方面,去噪方法的一些改进对实际资料能产生明显效果,说明研究思路的稍作改变看似简单,处理结果有时却会发生质变.这里,笔者不具体研究某一种处理方法的优劣,而是结合以上实例和在处理项目中的体会,明确一种整体化研究思路并对其进行相应分析. 1 地震资料处理与研究思路

地震资料处理是地震勘探工作中的重要一环,处理工作中,善于思考和运用合理的研究思路不仅可以高效地解决处理难题,而且可以促进技术、方法的创新.要做到真正地熟练掌握处理工作,不仅是善于运用处理软件解决常规生产任务,而且要善于灵活运用各种知识解决一些要求思路新、见解独特的地质或物探难题,而探索处理工作中的研究思路、提高处理思维能力就成为处理工作的重中之重.好的研究思路是能够将掌握的物探知识在处理工作中灵活应用,在解决各种处理问题时都觉得思路开阔、得心应手,图 1是珠江口盆地A区块原始资料不同类型多次波衰减思路,可见,即使在面对多次波衰减等常规处理问题时,也可以尝试各种研究思路并综合优选出最佳方案.

图 1 珠江口盆地A区块原始资料不同类型多次波衰减思路 Fig. 1 Method of different multiple attenuation in Pearl River Mouth Basin A block

地震资料处理是一个复杂的动态过程,面对不同的处理目的,研究思路就会相应改变,而有些处理员在处理工作中形成了思维定势,陷入了经验主义的误区,认为所谓的处理工作就是一个基本不变的处理流程,只需要在进行不同目的处理时,进行少量的模块和参数替换即可,这种想法在面对日益复杂的处理项目时是行不通的,而随着叠前深度偏移、逆时偏移、海上OBC等(Chandola et al., 2003倪成洲等,2008谢飞等,2008杨仁虎等,2010刘红伟等, 20102011杜启振等,2011刘礼农和张剑锋,2011马继涛等,2011段鹏飞等,2013吴帮玉等,2013严红勇和刘洋,2013刘文卿等,2013)新技术新方法的不断成熟,对处理工作提出了更高的要求,必须具有发散性思维、逆向思维等思考能力,尝试各种研究思路,进而快速有效地解决问题. 2 整体化研究思路分析

整体化研究思路就是要在资料处理中充分运用整体思维观念思考和解决问题,避免进入片面思维的误区.下面就整体化研究思路在处理工作中的应用进行相应分析. 2.1 处理项目的整体计划 2.1.1 整体计划的目的性和可执行性

(1)计划要强调目的性

处理项目是一项复杂的系统工程,从最初的数据解编、置观测系统、去噪等,到后期的偏移和偏移后处理,每一步都需要投入大量的时间进行模块、参数试验和前后资料对比分析等,而一般的三维处理项目的周期通常只有两三个月,这么短的时间要想出好的结果,就必须制定严密的计划并严格按计划执行,否则,当遇到原始资料品质较差、地质条件复杂(李振春等,2006岳玉波等,2012)、甚至个别参数出现调试失误等一些问题时,就会延缓处理进度或使后续处理时间紧张,进而影响最终处理效果.

项目计划要合理分配资源,重要的处理步骤要多投入时间和人力、物力,这样才能解决主要矛盾,处理结果也不会比预期有太大偏差.严密的项目计划也保证了正确的工作方向和有序的工作进度,防止浪费大量资源做无用功.例如,在处理过程中的速度分析方面,需要投入更多的时间进行细致研究,不同的地质认识对速度场的解释影响较大,与地质人员进行有效沟通以强化对研究区的地质认识,条件允许情况下与地质人员共同解释速度场,这在地质复杂区和重要目标区会取得明显效果.图 2是珠江口盆地A区块不同百分比偏移速度场叠加剖面对比图,三个对比剖面的处理模块和参数完全一致,从图中看出,100%速度场剖面总体较合理但局部需要调整,96%速度场剖面的断面不清晰,104%速度场剖面的断层位置不合理,将经过细致调整后的速度场进行偏移,然后输出速度线,经过多次反复的速度迭代分析,就能较好地把握住工区速度场特征,更精确的偏移速度场使得新处理资料(图 3)在花状断层区比老资料成像更清楚、合理,地层的接触关系更可靠.

图 2 珠江口盆地A区块不同百分比偏移速度场叠加剖面对比图 Fig. 2 Comparison of different percentage migration velocity field profiles in Pearl River Mouth Basin A block

图 3 珠江口盆地A区块相同测线新老偏移速度场叠加剖面对比图 Fig. 3 Comparison of new and old migration velocity field profiles in Pearl River Mouth Basin A block

(2)计划的可执行性

制定项目处理计划的目的是为了使目标变得更加具有可执行性,严密的计划必然使每天的工作内容都安排得充实,避免忘掉重要的事情和漫无目的地试验,每周、每天都要首先把项目最重要的事情列出来,只要保证计划的有效执行,一定会得到令人满意的结果.

计划也要强调执行标准,例如处理流程的标准化,处理质量控制标准化(郝振江和陈小宏,2007赵玉莲等,2013)等等,这样在执行计划时才能有的放矢、防止个别环节出现差错影响后续处理结果甚至造成生产事故.对项目组成员来说,工作不仅要符合基本的处理标准,更应当对自我有更高的要求,向行业专家学习工作作风、良好习惯,这样不仅提高工作技能、完成项目要求,更重要的是不断开阔眼界、提高综合素质. 2.1.2 生产与试验处理的整体协调

(1)生产时间与试验时间相协调

资料处理的生产任务是按时按质完成项目,而试验是为了对参数、模块甚至流程进行对比和优选,以保证最好的处理结果,两者的目的是一致的,但很多时候由于处理时间的限制,试验时间与生产进度往往有冲突,这时就要对两者有合理的协调,生产时间方面,要对处理软件的使用、生产作业的编制(王秀闽等,2007)进行合理的安排,以提高效率;在试验时间方面,应利用等待生产作业的间隙来思索试验方案、处理模块和参数的选择、试验作业的编制.

(2)生产的严谨与试验的灵活相协调

生产作业中的模块、参数的选择相对严谨,要求能够保证处理结果的稳定性和合理性.生产的严谨是保证每一步处理后的资料在完全准确的情况下再进行下一步生产处理,这样才能从根本上保障整体计划的顺利执行.一旦某一步出现问题,而后续未及时发现,则会浪费掉大量时间去查找错误,解决办法是在生产处理时要认真细心、多总结经验、多与人沟通、力争少出错.试验作业中则应尽可能灵活选取更多的模块、参数甚至流程进行尝试,这样才能从中优选最佳方案进行正式生产. 2.2 技术、经验、理论的整体匹配

资料处理人员要具有处理方面的技术、经验、理论等综合素质,从技术方面来说,要求处理员对处理软件、运算设备等能够熟练掌握,还要求有大量的项目研究和管理经验,无论是二维、三维或多波多分量项目,还是复杂山地、海陆过渡带、深水等(吴琼等,2008常旭等,2008李中元等,2011张英德等,2012张山等,2012)处理项目,只有承担更多的处理项目,才能积累足够的实践经验,最终做到触类旁通,在经验方面达到专业人才的要求.理论比技术和经验的要求更进一步,当我们做完一个处理项目后,要善于从理论角度与以往项目进行对比总结,才能达到实践与理论相结合、不断提升整体技术水平. 2.2.1 技术、经验、理论的匹配程度决定处理效果

同一处理项目,可以选择不同的处理思路、模块,处理水平的高低取决于技术、经验、理论的匹配程度,三者的合理匹配才能花费更少的时间、精力并得到更好的处理效果.

经验、技术与理论的提升是一个由量变到质变的过程,这些只有在不断的资料处理实践中得到逐步地提升.成熟的处理员除了要不断学习专业理论,更要有大量的生产、科研项目经验的积累,没有量的积累难以达到质变的效果.处理工作在某种程度上讲是单调枯燥的,在面对几百种处理模块、大量的处理和试验参数、海量的原始数据、紧张的工作节奏时,极易导致处理思路单一、进度缓慢,造成处理效果不佳.

这里强调的技术、经验和理论,不仅指专业方面,更包括项目管理、工作兴趣等各方面,例如,在进行长周期复杂处理项目时,应当把项目目标分解为多个个小的目标,同时分清楚主次,做好时间和精力分配,在完成这些小目标后,最终处理效果一定不会与处理目的相差太远.因为再复杂的问题也是从一些小问题堆积而成的,认真执行项目计划、不放过可能影响处理效果的小问题、多做试验对比,处理效果则会逐步显现,有时甚至能找到解决一些处理难题的更好方法.即使在进行相对枯燥的速度谱解释时,也可以针对偏移、叠加等不同需要选取不同方法,从而激发兴趣、提高效率. 2.2.2 处理新技术与常规技术的协调

新技术与常规技术两者是相辅相成的关系,新技术除了有其技术优势外,也有相应理论相支撑,开发应用好会成为主流技术,现今广泛应用的三维PSTM处理就是由新技术发展成熟的,新技术应用得合理可以提高工作效率和处理质量,例如,图 4是渤海湾盆地B区块常规资料与OBC采集处理资料对比图,通过对比可见,采用新的采集处理技术后使得地震资料的成像质量有了较大的提高,以前的成像模糊区变得更清晰,地层内部信息更丰富.

图 4 渤海湾盆地B区块常规资料与OBC采集处理资料对比图 Fig. 4 Comparison of conventional data and OBC data in Bohai Bay Basin B block

新技术也有理论不完善、技术缺陷或成本高昂等劣势,如果过分强调新技术,也会忽视常规技术理论基础完备、技术成熟、成本合理等优势.无论是新技术还是常规技术,都应重视其各自的优势,在解决不同的处理问题时,如果时间允许,要多进行新、老技术的试验对比,合理选择处理技术,达到二者优势互补的目的. 2.3 不同专业人才的整体合作

处理人员要不断地提升自身的专业能力、更快地掌握处理复杂问题的综合技能,必然要求物探、地质、计算机等各类人才之间加强合作沟通、互相学习才能更快地掌握处理复杂问题的综合技能.而在面对日益更新的技术及理论发展潮流和复杂的地质条件时,也要求各专业人才互相合作、共同攻坚克难.
2.3.1 处理解释一体化

交叉学科人才交流合作,有利于开拓研究思路、打破固有的思维模式、找到更为有效的解决问题的方法,例如,在面对同一复杂地质问题时,处理人员与解释人员思考的角度不同,思考的结论也有差别,各自偏重于某一方面,一体化的目的就是要接受彼此的观点,最终得出更全面正确的认识.图 5是珠江口盆地A区块相同测线新老处理结果对比图,图中显示两条边界大断层,经过处理解释人员的深入交流,认为此处为两条断层的应力转换区,新处理结果的阶梯状断层比老资料的铲式断层在地质认识上更合理,统一认识后,新解释的速度场更准确,偏移叠加剖面成像更清晰、波组特征更清楚.

图 5 珠江口盆地A区块相同测线新老处理结果对比图 Fig. 5 Comparison of new and old data processing results in Pearl River Mouth Basin A block

从另一方面来说,一体化合作可以防止出现信息不对称带来的认识偏差,同样的地震资料,处理人员看到的是物探信息且擅长解决物探问题,解释人员看到的则是地质信息而擅长解决地质问题,而处理工作要解决的是综合问题,必然要求资料处理时掌握更多的物探、地质综合信息,处理解释一体化给处理和解释人员提供了信息交流平台,解决更多综合难题.同时,一体化也可以整合人才资源、实现企业效益最大化. 2.3.2 项目成员相互协作

根据项目成员每个人的性格特点、专长来合理分配工作任务,发挥各自的优点.例如,新老处理员结合是经验与技术的合理优化,老员工的优势是技术经验丰富,劣势是精力和时间不够用,而年轻员工则正好相反,两者的合作刚好取长补短,在面对复杂处理问题时,老员工提供经验帮助,而新员工则充分运用时间多做试验和生产处理,最终保证资料处理质量和进度. 2.4 处理软硬件资源的整体协调

资料处理的对象往往是海量数据,因此对处理软件及硬件的要求很高,处理人员必须熟悉计算机系统及硬件等性能,根据不同的处理需求应用不同的软、硬件资源,达到资源的合理利用. 2.4.1 软件与硬件的合理配置

目前处理行业中选用的处理软件种类繁多,工作站版的处理软件较为普遍,优点是性能稳定、运算量大、效率高,常见处理软件系统有OMEGA、CGG、FOCUS等;微机版处理软件具有操作简便且硬件要求不高的优点,但性能不够稳定、运算量小,更适合现场处理或简单的处理试验.因此应根据研究需要来选择不同的处理软件,同时,不同的处理软件则应配置合适的硬件资源.例如,根据需求在进行逆时偏移等大量数据处理时,运用CPU/GPU联合作为计算核心(李博等,2010匡斌等,2011),而非单一使用CPU作计算,既节省成本又提高了效率.

随着资料处理难度的加大,对处理软件的要求也在不断提高,每一次处理软件升级后,对硬件的配置也提高了要求.处理员、系统人员、设备管理员之间应及时沟通,使得处理软件与硬件之间合理配置,而目前处理硬件的升级赶不上软件升级的现象较为普遍,这点应当逐渐被重视起来,避免仅由于软硬件配置不合理造成的处理进度缓慢、处理效果不佳等问题. 2.4.2 各处理项目间的软硬件资源协调

根据处理需求,普遍存在多个处理项目并行的问题,但一个处理单位软硬件资源是有限的,这时应协调好每个项目的处理进度,尽量将进度错开,避免多个项目同时等待使用少量偏移节点情况,造成软硬件资源分配不合理现象出现.如果发生计算节点等资源紧张问题时,应根据项目轻重缓急,依次排序处理,同时,应将计算量大、周期较长的处理作业安排在晚上或者周末进行运算,避免产生工作时段计算节点不够用,而非工作时段大量节点闲置的资源浪费现象. 3 结 论

(1)处理项目的周期相对较短,短时间内要想得到好的处理效果,就必须制定严密的整体计划,计划要强调目的性和可执行性,这可以保障处理工作质量和进度,也保障了资料处理中的试验与生产时间的合理协调.

(2)处理人员要具有技术、经验、理论等综合素质,技术上,要对处理软件、运算设备等熟练掌握;要有大量的项目研究和管理经验,善于钻研物探及地质理论,提高理论联系实践的能力.

(3)资料处理中,物探、地质、计算机等专业人才之间加强沟通交流,可以更好更快地解决物探地质综合难题,同时也可以整合人才资源、实现企业效益最大化.

(4)处理人员要熟悉各种处理软件、硬件等方面的性能,根据不同的处理需求应用不同的软、硬件资源,达到资源的整体有效利用.

(5)资料处理可以被看成是一个整体系统,该系统由处理员、处理算法、处理硬件等各要素组成.若处理目的简单,即组成系统的要素的性质和要素间的关系相对简单,这是一个稳定系统,这种情况下对资料处理的研究思路要求相对简单.随着处理目的难度加大,该系统中各要素的性质、要素间的关系更加不稳定,即为不稳定系统,这时就需要运用整体化研究思路将系统中的各要素变得更加稳定,进而有效解决资料处理难题. 致 谢 本研究得到中海油服物探事业部处理解释中心的陈继宗经理、张善刚高工等人的帮助,在此表示衷心感谢!

参考文献
[1] Chandola S K, Ramakrishna K, Sharma A K, et al. 2003. Dynamic fold solutions for improved sub-surface imaging: A transition zone 3D-OBC survey in East Coast of India[Z]. SEG Technical Program Expanded Abstracts, 22.
[2] Chang X, Liu Y K, Du X D, et al. 2008. Seismic imaging under the irregular deep water bottom[J]. Chinese Journal of Geophysics (in Chinese), 51(1): 228-234.
[3] Du Q Z, Li F, Hou B, et al. 2011. Angle-domain migration velocity analysis based on elastic-wave Kirchhoff prestack depth migration [J]. Chinese Journal of Geophysics (in Chinese), 54(5): 1327-1339, doi: 10.3969/j.issn.0001-5733.2011.05.022.
[4] Duan P F, Cheng J B, Chen S P, et al. 2013. Local angle-domain ray tracing and prestack depth migration in TI medium[J]. Chinese Journal of Geophysics (in Chinese), 56(1): 269-279, doi: 10.6038/cjg20130128.
[5] Hao Z J, Chen X H. 2007. A study on time-lapse seismic cross-equalization QC method based on correlation parameters[J]. Progress in Geophysics (in Chinese), 22(3): 929-933, doi: 10.3969/j.issn.1004-2903.2007.03.040.
[6] Kuang B, Du J X, Wang H Z, et al. 2011. 3D wave equation prestack depth migration on Graphic Processing Unit[J]. Oil Geophysical Prospecting (in Chinese), 46(5): 705-709.
[7] Li B, Liu H W, Liu G F, et al. 2010. Computational strategy of seismic pre-stack reverse time migration on CPU/GPU[J]. Chinese Journal of Geophysics (in Chinese), 53(12): 2938-2943, doi: 10.3969/j.issn.0001-5733.2010.12.017.
[8] Li Z C, Sun X D, Liu H. 2006. Common reflection surface stack for rugged surface topography[J]. Chinese Journal of Geophysics (in Chinese), 49(6): 1794-1801.
[9] Li Z Y, Wu D L, Cai L, et al. 2011. Processing techniques for the low S/N seismic data of the complex mountainous area in MD[J]. Geophysical Prospecting for Petroleum (in Chinese), 50(6): 595-599.
[10] Liu H W, Liu H, Zou Z, et al. 2010. The problems of denoise and storage in seismic reverse time migration[J]. Chinese Journal of Geophysics (in Chinese), 53(9): 2171-2180, doi: 10.3969/j.issn.0001-5733.2010.09.017.
[11] Liu H W, Liu H, Li B, et al. 2011. Pre-stack reverse time migration for rugged topography and GPU acceleration technology[J]. Chinese Journal of Geophysics (in Chinese), 54(7): 1883-1892, doi: 10.3969/j.issn.0001-5733.2011.07.022.
[12] Liu L N, Zhang J F. 2011. Wave equation prestack depth migration method in 3D VTI media[J]. Chinese Journal of Geophysics (in Chinese), 54(11): 2906-2915, doi: 10.3969/j.issn.0001-5733.2011.11.020.
[13] Liu W Q, Wang X W, Liu H, et al. 2013. Application of velocity modeling and reverse time migration to subsalt structure[J]. Chinese Journal of Geophysics (in Chinese), 56(2): 616-625, doi: 10.6038/cjg20130225.
[14] Liu Y, Wang D, Liu C, et al. 2011. Weighted median filter based on local correlation and its application to poststack random noise attenuation[J]. Chinese Journal of Geophysics (in Chinese), 54(2): 358-367, doi: 10.3969/j.issn.0001-5733.2011.02.012.
[15] Ma J T, Mrival S K, Chen X H, et al. 2011. OBC multiple attenuation technique using SRME theory[J]. Chinese Journal of Geophysics (in Chinese), 54(11): 2960-2966, doi: 10.3969/j.issn.0001-5733.2011.11.026.
[16] Ni C Z, Quan H Y, Chen G, et al. 2008. A new method for high-precision first break secondary positioning—Search method[J]. Oil Geophysical Prospecting (in Chinese), 43(2): 131-133, 157.
[17] Wang X M, Liu H, Liu Y J, et al. 2007. Editing and running interactive job in seismic data processing system[J]. Progress in Geophysics (in Chinese), 22(3): 860-864, doi: 10.3969/j.issn.1004-2903.2007.03.030.
[18] Wu B Y, Wu R S, Gao J H. 2013. Survey sinking prestack depth migration using local cosine basis beamlets[J]. Chinese Journal of Geophysics (in Chinese), 56(2): 635-643, doi: 10.6038/cjg20130227.
[19] Wu Q, Yang C C, Zhang W Z, et al. 2008. Research on processing technology for matching seismic data acquired in seashore and marine prospect[J]. Progress in Geophysics (in Chinese), 23(3): 761-767.
[20] Xie F, Chang X, Liu Y K. 2008. Converted-wave pre-stack time migration in volcanic area[J]. Chinese Journal of Geophysics (in Chinese), 51(6): 1899-1908.
[21] Yan H Y, Liu Y. 2013. Acoustic prestack reverse time migration using the adaptive high-order finite-difference method in time-space domain[J]. Chinese Journal of Geophysics (in Chinese), 56(3): 971-984, doi: 10.6038/cjg20130325.
[22] Yang R H, Chang X, Liu Y K. 2010. The influence factors analyses of imaging precision in pre-stack reverse time migration[J]. Chinese Journal of Geophysics (in Chinese), 53(8): 1902-1913, doi: 10.3969/j.issn.0001-5733.2010.08.016.
[23] Yuan Y H, Wang Y B, Liu Y K, et al. 2013. Non-dyadic Curvelet transform and its application in seismic noise elimination[J]. Chinese Journal of Geophysics (in Chinese), 56(3): 1023-1032, doi: 10.6038/cjg20130330.
[24] Yue Y B, Li Z C, Qian Z P, et al. 2012. Amplitude-preserved Gaussian beam migration under complex topographic conditions[J]. Chinese Journal of Geophysics (in Chinese), 55(4): 1376-1383, doi: 10.6038/j.issn.0001-5733.2012.04.033.
[25] Zhang S, Bai J H, Li Y P, et al. 2012. offset VSP data processing method for complex piedmont[J]. Geophysical Prospecting for Petroleum (in Chinese), 51(2): 164-171.
[26] Zhang Y D, Peng J Y, Hao L Y, et al. 2012. Time-depth conversion difficulties and technical countermeasures for overseas deepwater and complex geological conditions[J]. Progress in Geophysics (in Chinese), 27(4): 1484-1492, doi: 10.6038/j.issn.1004-2903.2012.04.023.
[27] Zhao Y L, Wang Y C, Wang Z Q, et al. 2013. Application of seismic attributes QC technology in seismic data processing[J]. Progress in Geophysics (in Chinese), 28(2): 860-868, doi: 10.6038/pg20130237.
[28] 常旭, 刘伊克, 杜向东,等. 2008. 深水崎岖海底地震数据成像方法与应用[J]. 地球物理学报, 51(1): 228-234.
[29] 杜启振, 李芳, 侯波,等. 2011. 角度域弹性波Kirchhoff叠前深度偏移速度分析方法[J]. 地球物理学报, 54(5): 1327-1339, doi: 10.3969/j.issn.0001-5733.2011.05.022.
[30] 段鹏飞, 程玖兵, 陈三平,等. 2013. TI介质局部角度域射线追踪与叠前深度偏移成像[J]. 地球物理学报, 56(1): 269-279, doi: 10.6038/cjg20130128.
[31] 郝振江, 陈小宏. 2007. 基于相关参数的时移地震互均化质量监控方法研究[J]. 地球物理学进展, 22(3): 929-933, doi: 10.3969/j.issn.1004-2903.2007.03.040.
[32] 匡斌, 杜继修, 王华忠,等. 2011. 基于GPU计算平台的三维波动方程叠前深度偏移[J]. 石油地球物理勘探, 46(5): 705-709.
[33] 李博. 刘红伟. 刘国峰,等. 2010. 地震叠前逆时偏移算法的CPU/GPU实施对策[J]. 地球物理学报, 53(12): 2938-2943, doi: 10.3969/j.issn.0001-5733.2010.12.017.
[34] 李振春, 孙小东, 刘洪. 2006. 复杂地表条件下共反射面元(CRS)叠加方法研究[J]. 地球物理学报, 49(6): 1794-1801.
[35] 李中元, 邬达理, 蔡俩,等. 2011. MD复杂山地低信噪比地震资料处理技术[J]. 石油物探, 50(6): 595-599.
[36] 刘红伟, 刘洪, 邹振,等. 2010. 地震叠前逆时偏移中的去噪与存储[J]. 地球物理学报, 53(9): 2171-2180, doi: 10.3969/j.issn.0001-5733.2010.09.017.
[37] 刘红伟, 刘洪, 李博,等. 2011. 起伏地表叠前逆时偏移理论及GPU加速技术[J]. 地球物理学报, 54(7): 1883-1892, doi: 10.3969/j.issn.0001-5733.2011.07.022.
[38] 刘礼农, 张剑锋. 2011. 三维各向异性介质中的波动方程叠前深度偏移方法[J]. 地球物理学报, 54(11): 2906-2915, doi: 10.3969/j.issn.0001-5733.2011.11.020.
[39] 刘文卿, 王西文, 刘洪,等. 2013. 盐下构造速度建模与逆时偏移成像研究及应用[J]. 地球物理学报, 56(2): 616-625, doi: 10.6038/cjg20130225.
[40] 刘洋, 王典, 刘财,等. 2011. 局部相关加权中值滤波技术及其在叠后随机噪声衰减中的应用[J]. 地球物理学报, 54(2): 358-367, doi: 10.3969/j.issn.0001-5733.2011.02.012.
[41] 马继涛, Mrinal S K, 陈小宏, 等. 2011. 海底电缆多次波压制方法研究[J]. 地球物理学报, 54(11): 2960-2966, doi: 10.3969/j.issn.0001-5733.2011.11.026.
[42] 倪成洲, 全海燕, 陈刚,等. 2008. 一种高精度初至波二次定位新方法—搜索法[J]. 石油地球物理勘探, 43(2): 131-133, 157.
[43] 王秀闽, 刘洪, 刘永江,等. 2007. 地震数据处理系统中交互作业编辑与运行[J]. 地球物理学进展, 22(3): 860-864, doi: 10.3969/j.issn.1004-2903.2007.03.030.
[44] 吴帮玉, 吴如山, 高静怀. 2013. 基于局部余弦基小波束的观测系统沉降法叠前深度偏移[J]. 地球物理学报, 56(2): 635-643, doi: 10.6038/cjg20130227.
[45] 吴琼, 杨长春, 张文忠,等. 2008. 过渡带拼接地震资料处理方法研究[J]. 地球物理学进展, 23(3): 761-767.
[46] 谢飞, 常旭, 刘伊克. 2008. 火山岩地震屏蔽层的转换波叠前时间偏移成像[J]. 地球物理学报, 51(6): 1899-1908.
[47] 严红勇, 刘洋. 2013. 基于时空域自适应高阶有限差分的声波叠前逆时偏移[J]. 地球物理学报, 56(3): 971-984, doi: 10.6038/cjg20130325.
[48] 杨仁虎, 常旭, 刘伊克. 2010. 叠前逆时偏移影响因素分析[J]. 地球物理学报, 53(8): 1902-1913, doi: 10.3969/j.issn.0001-5733.2010.08.016.
[49] 袁艳华, 王一博, 刘伊克,等. 2013. 非二次幂Curvelet变换及其在地震噪声压制中的应用[J]. 地球物理学报, 56(3): 1023-1032, doi: 10.6038/cjg20130330.
[50] 岳玉波, 李振春, 钱忠平,等. 2012. 复杂地表条件下保幅高斯束偏移[J]. 地球物理学报, 55(4): 1376-1383, doi: 10.6038/j.issn.0001-5733.2012.04.033.
[51] 张山, 白俊辉, 李彦鹏,等. 2012. 复杂山前带有偏VSP数据处理方法研究[J]. 石油物探, 51(2): 164-171.
[52] 张英德, 彭佳勇, 郝立业,等. 2012. 海外深水复杂地质条件下时深转换难点及技术对策[J]. 地球物理学进展, 27(4): 1484-1492, doi: 10.6038/j.issn.1004-2903.2012.04.023.
[53] 赵玉莲, 王宇超, 王振强,等. 2013. 地震属性质量控制技术在地震数据处理中的应用[J]. 地球物理学进展, 28(2): 860-868, doi: 10.6038/pg20130237.