2017, Vol. 32
2. 中国石化胜利油田分公司物探研究院, 东营 257022
2. Geophysical Research Institute of Shengli Oil Field Branch, SINOPEC, Dongying 257022, China
地震勘探算法是各类采集、处理和解释软件系统的核心，是油气地球物理勘探开发的灵魂(魏福吉和徐维秀，2013).各类改进与创新算法的出现往往会推动油气勘探开发行业的发展和进步，如：时频分析(高静怀等，2004；刁瑞和冯玉苹，2012)、井间地震(冯玉苹等，2013)、保幅反褶积(刁瑞等，2013)、三维波动方程正演模拟(张金海等，2007)、叠前时间偏移、叠前深度偏移、逆时偏移(吴国忱和秦海旭，2014；李振春等，2014)、全波形反演(杨午阳等，2013)、微地震监测(姜福兴等，2006；刁瑞等，2017)等.随着高效勘探的进展，可控震源的使用越来越普遍，国内外学者对可控震源进行了深入研究.由于可控震源的机械装置和震动装置的非线性以及与大地耦合等特点，可控震源记录中存在谐波干扰(佟训乾等，2012)，尤其为提高生产效率，采用滑动扫描等高效采集技术施工时，谐波干扰更为明显，相应的谐波干扰压制技术随之不断发展.Ras等(1999)通过正演模拟滑动扫描记录中的谐波干扰，并选择合适的滑动扫描时间在频率域分离谐波干扰.曹务祥(2004)提出了组合扫描压制谐波，谐波干扰随扫描频宽的降低而减少，当扫描频宽达到一个倍频程时谐波干扰即消失.Cao(2010)利用地面力信号代替参考信号，采用褶积代替相关来压制谐波干扰，消除复杂近地表影响，提高数据品质.黄建平等(2012)提出了一种可以适用于滑动扫描技术的相移滤波法.韩文功等(2015)基于谐波干扰是谐波畸变信号与反射系数的褶积，提出了基于反射系数和地面力信号的谐波干扰压制方法.谐波干扰压制技术的不断进步，保障了可控震源高效采集的实施，降低了生产成本.
地球物理行业中各类编程语言很多，目前较为常用和流行的是C/C++、Fortran、Matlab、Object Pascal、Java、C#、Visual Basic等(任慧等，2007；李勤等，2010).不同的编程语言各有优缺点，刻意区分语言的好坏并没有实际意义，重要的是如何在实际编程中发挥他们各自的优点，实现地震勘探算法的高效集成与应用(潘大夫等，2009；王鹏，2014).在油气地球物理算法研究中，常常采用C/C++、Fortran、Matlab语言等进行算法编程，形成较为独立的地震勘探算法和实用程序，较好地完成了地震数据采集、处理和解释等工作，在油气地球物理勘探中发挥了重要作用(王宏琳和陈继红，2010；徐维秀等，2016).在综合前人研究的基础上，笔者将地震勘探算法特点归纳为以下几个方面：(1) 各类算法程序独立运行，没有统一的软件平台和专门的软件界面；(2) 通常在算法程序内部进行读写操作及参数修改，没有数据输入、输出界面；(3) 各算法相互独立，由不同的科研人员采用不同的编程语言进行开发；(4) 输入、输出数据格式复杂多样，数据格式不统一，不同科研人员之间无法有效交流；(5) 数据显示和交互操作功能欠缺，需要借助第三方软件平台等.
地震勘探算法的特点制约了算法的成果转化和推广应用，主要包括以下几个方面的问题：(1) 由于单一化的处理功能，不能形成较为完善的处理流程，限制了技术的推广与应用；(2) 复杂、繁琐的数据格式转换与过程处理，制约了算法的应用与完善开发；(3) 没有商业化软件界面封装，不能提供友好的图形化显示功能，难以吸引用户使用，科研成果转化和软件系统产业化比较困难.针对地震勘探算法的特点和存在的问题，开展了插件模式集成及相关技术的研究，将不同编程语言实现的地震勘探算法，通过插件模式进行算法模块集成，设计形成了有效的功能技术模块，并保证软件模块的可重用性和稳定性.基于插件模式集成技术形成了完整、有效的软件系统，以可控震源正演模拟与谐波压制技术为例，验证了算法集成的正确性和有效性，在实际应用中取得了良好的效果.
1 插件模式集成技术研究 1.1 插件开发机制插件是可独立完成某一个或一系列功能的模块，通常插件由宿主程序加载，不能独立运行.宿主程序是承载插件运行的环境，为插件提供基本服务，通常情况下插件由宿主程序管理和控制.宿主程序由核心开发人员编写，并预先定义好系统所提供的基本服务接口和插件接口，然后由其他开发人员根据系统插件接口实现编写插件功能.插件的机制如图 1所示.
|
图 1 插件机制示意图 Figure 1 The plug-in mechanism schematic diagram |
插件服务引擎提供构建软件系统需要的抽象插件处理机制，为应用开发提供插件处理接口规范和基类，实现基于插件的地震数据采集应用程序框架和地震勘探算法管理框架，支持各种地震勘探算法应用扩展.主框架定义插件的加载、卸载、通讯接口规则等，实现对插件库的管理操作，插件方则根据接口的规则实现一套插件接口，以供主框架集成.
地震数据采集工程软件平台支持功能性插件和方法性插件，功能性插件允许具备自己的功能操作和显示，集成过程相对复杂；而方法性插件只提供算法的输入和输出，集成过程较简单、快捷.针对具体地震勘探算法的特点，采用方法性插件的方式进行地震勘探算法模块的集成.图 2为基于插件的应用集成框架结构示意图，观测系统设计、近地表资料分析、正演模拟与照明分析、地质/地震建模等各类地震勘探算法，均可通过插件的形式集成到主框架软件系统中.
|
图 2 基于插件的应用集成框架结构 Figure 2 Application integration framework based on plug-in structure diagram |
算法插件需要继承于ISgmMathRegister类、ISgmMathShow类和ISgmMathEvaluate类.ISgmMathRegister类是算法插件注册接口类，用于获取算法插件的注册信息. ISgmMathShow类是算法插件的显示接口类，用于调用算法插件中的显示界面，展示算法的效果. ISgmMathEvaluate类是算法插件的计算接口类，用于调用算法插件的计算方法.插件的接口包括插件信息接口、加载插件信息接口和卸载插件信息接口.
由于编程人员的习惯和地震勘探算法需求的不同，使用的编程语言也各不相同，目前常用的编程语言有C/C++、Fortran、Matlab语言等.如果将这些编程语言都转换成C/C++语言，则重新编码工作量巨大、集成周期较长，不利于算法的应用和推广.在此采用的针对性策略是：采用插件模式进行算法模块集成，由于Fortran语言和C/C++语言可以跟QT进行混合编程，所以对于Fortran和C/C++语言编写的程序，生成可执行文件供QT程序调用即可；由于Matlab跟QT无法进行混合编程，所以对于Matlab语言编写的程序，生成动态库文件，然后新建C++项目调用该动态库，生成可执行文件供QT程序调用.插件模式集成技术实现快捷、应用方便，不会增加额外的工作量，集成周期比较短，具有较高的实用价值，有利于各类地震勘探算法的应用和推广.
2 SEGY数据服务地震勘探算法中常用到的数据格式为SEGY格式，为此我们对SEGY格式进行了解析，以此实现在各种不同格式地震数据的统一管理.标准的SEGY格式包括由3200字节的字符串和400字节的二进制组成的卷头字段，以及240字节的道头字和地震道数据(表 1).3200字节的字符串包括ASCII码和EBCDIC码两种方式.SEGY数据格式包括：4字节IBM浮点SEGY、4字节整型SEGY、2字节整型SEGY、4字节IEEE浮点SEGY和1字节整型SEGY.
|
表 1 SEGY文件结构 Table 1 Tht file structure of SEGY |
通过对SEGY格式的解析，设计了SEGY格式访问接口，如图 3所示，该访问接口提供了标准SEGY格式的输入与输出.由于地震勘探算法的差异，以及编程的不规范等问题，有的算法内部只支持无卷头、无道头的二进制SEGY格式数据，有的算法只支持有道头格式的SEGY格式数据，有的算法支持的是不标准的SEGY格式数据.由于数据格式的差异，造成了复杂的格式转换难题，不利于算法的应用和数据的统一管理.
|
图 3 SEGY数据访问接口类图 Figure 3 SEGY data access interface class diagram |
针对数据格式差异的难题，在不修改原有地震勘探算法的基础上，提供了一个标准转换接口，使得地震勘探算法的输入、输出数据格式标准化，集成后的地震勘探算法使用方便，数据统一管理，大大提高了工作效率.以地震数据采集工程软件提供的地震数据管理服务为基础，以地震图形组件库为支撑，开发了“SEGY地震数据浏览器”，支持无卷头和有卷头各类地震数据的浏览显示.
3 模块集成测试及应用效果 3.1 可控震源正演算法集成及应用可控震源扫描技术分为线性扫描、非线性扫描和伪随机扫描，通过可控震源仿真模拟可优选最佳参数，从而提高地震勘探分辨率(王忠仁等，2006).基于算法插件集成模式完成了地震勘探算法的集成，以可控震源滑动扫描正演模拟算法为例，进行算法插件集成模式的测试与应用.地震数值模拟是地震勘探和地震学的重要基础(吴国忱等，2005；李美娟等，2016).可控震源滑动扫描正演算法是由Fortran语言编写的，首先提取算法中的输入参数，并以参数文件的形式进行读写，然后将该算法编译成可执行文件，在制作算法插件时调用即可.滑动扫描正演模拟算法插件的实现类图如图 4所示.
|
图 4 滑动扫描正演模拟算法插件实现类图 Figure 4 Sliding sweep forward modeling algorithm implementation class diagram |
可控震源扫描正演算法集成的具体流程为：首先，实现算法插件的接口；其次，直接调用已制作好的可执行文件，完成算法插件的制作；再次，在插件制作完成后，进行插件的注册；最后，在注册完成后在算法导航树上，即可使用加载的可控震源滑动扫描正演模拟算法插件.目前，SeisWay 3.0软件系统已完成可控震源单炮正演模拟、可控震源滑动扫描正演模拟、可控震源混叠数据成像算法的集成工作.可控震源滑动扫描正演模拟集成效果如图 5所示，软件系统提供对算法插件的计算、信息查看和删除等操作.点击“计算”按钮，打开参数输入界面，可以手动填写参数信息，也可以导入已有的参数文件，在此基础上进行修改完善，填写完成后保存参数文件即可.点击“确定”按钮即开始运算，运算完成后，在运行窗口中可显示运行信息.
|
图 5 可控震源滑动扫描正演模拟集成效果 Figure 5 Slide vibroseis scanning forward modeling integration module |
插件模块运行完成后，通过SeisWay 3.0软件系统提供的SEGY显示功能，可以对计算结果进行查看，图 6a是可控震源滑动扫描正演模拟相关前炮记录，图 6b是可控震源滑动扫描正演模拟相关后炮记录.通过室内可控震源滑动扫描正演模拟方法，能够模拟多种可控震源高效采集技术的数据，结合野外实际采集地震数据，分析反射波对应的地形特征和波场特性.另外，通过模拟的炮记录和波场快照记录，研究高效采集数据谐波形成的机理，可为野外谐波压制和分析提供指导.
|
图 6 可控震源滑动扫描正演模拟相关前(a)、后(b)单炮记录 Figure 6 The correlation shot befor (a) and after (b) of slide vibroseis scanning forward modeling |
高密度宽方位三维地震勘探技术是获取高品质地震资料的有效途径，可控震源施工是最有效、经济实用的方法，已成为地震勘探中重要的激发震源之一.在国际地震勘探领域，可控震源已成为最主要的激发震源，并被广泛应用于油气勘探开发中(Sicking et al., 2009；Martin and Munoz, 2010).根据可控震源的激发原理，由于相关过程而形成谐波虚反射，谐波干扰严重影响地震资料的品质，必须进行有效的压制(刘世文等，1998；Wei and Phillips, 2010).王忠仁等(2010)利用匹配伪随机序列偶调制正弦载波信号作为可控震源驱动信号进行扫描激发，得到两个匹配的原始地震记录，对两个匹配的原始地震记录分别解码后再进行叠加形成最终的综合解码地震剖面.胡江涛等(2013)提出了谐波应该被看作有效信号的观点，然后扩展了纯相位滤波器方法，并用于谐波信号的分离，最后对分离后的谐波记录分别进行相关并叠加得到最终记录.林娟等(2014)尝试应用同步扫描邻炮干扰压制方法滤除滑动扫描资料上的谐波干扰，取得了较好的效果.曲英铭等(2016)提出了一种基于最小二乘思想的自适应匹配预测滤波的谐波压制方法.目前，SeisWay 3.0软件系统已完成可控震源同步交涉干扰压制、自适应谐波压制、可控震源相移法谐波压制、可控震源分频中值滤波谐波压制等算法的插件集成工作.可控震源自适应谐波压制模块的集成效果如图 7所示.
|
图 7 可控震源谐波压制算法模块集成效果 Figure 7 Vibroseis harmonic suppression algorithm module integration module |
图 8a是可控震源自适应谐波压制前单炮记录，图 8b是可控震源自适应谐波压制后单炮记录.对比噪声压制前、后的单炮记录，可控震源谐波压制后，单炮记录的信噪比明显提升，有效提高了地震资料的品质.通过算法插件模式集成后，可控震源谐波压制技术系列能够方便应用和有效推广.
|
图 8 可控震源谐波压制前(a)、后(b)单炮记录对比 Figure 8 The shot record before (a) and after (b) of vibroseis harmonic suppression |
通过对地震勘探算法的现状研究和分析，研究了地震勘探算法插件模块集成中的关键技术，研发了较为广泛的地震勘探算法开发方式，形成了地震勘探算法的模块集成思路，实现了C++、Fortran、MATLAB等不同编程语言算法的插件模块集成，以可控震源正演模拟与谐波压制技术为例，验证了算法集成的正确性和有效性.
4.2同步震源激发和滑动扫描是可控震源高效采集技术，随着可控震源高效采集市场的不断扩大，在提高施工效率的同时产生谐波干扰和邻炮干扰，如何在提高采集效率的同时获取高品质地震资料是一个重要研究方向，可控震源谐波干扰压制、同步震源有效分离等技术将越来越重要.
4.3在地震勘探算法的开发与模块集成过程中，要遵循一定的标准和规则，采用统一的模块集成模式，使得集成模块具有操作简单、使用方便、风格一致的特点，有利于科研成果的转化和推广应用.另外，开放自由的插件模块加载、卸载方式，有利于自主研发软件系统的算法扩展，能够最大程度的促进科研成果的推广应用，具有一定的实用价值和推广前景.
致谢 感谢中国石油大学(华东)李振春教授、黄建平教授、曲英铭博士、黄金强博士的鼎力帮助和大力支持！| [] | Cao W X. 2004. Using array sweeping to suppress harmonic distortion[J]. Oil Geophysical Prospecting (in Chinese), 39(6): 711–715. |
| [] | Cao W X. 2010. To attenuate harmonic distortion by the force signal of vibrator[C].//2010 SEG Annual Meeting. Denver, Colorado:Expanded Abstracts, 157-161. |
| [] | Diao R, Feng Y P. 2012. Hydrocarbon recognition based on spectrum absorption characteristics of reservoir[J]. Oil Geophysical Prospecting (in Chinese), 47(5): 766–772. |
| [] | Diao R, Shang X M, Rui Y J, et al. 2013. Study on spectral modeling deconvolution in time-frequency domain[J]. Lithologic Reservoirs (in Chinese), 25(1): 116–121. |
| [] | Diao R, Wu G C, Shang X M, et al. 2017. Joint correction method based on 3D seismic and surface microseismic data[J]. Chinese Journal of Geophysics(in Chinese), 60(1): 283–292. DOI:10.6038/cjg20170123 |
| [] | Feng Y P, Diao R, Xu W X, et al. 2013. Study on P and S wave separation method in cross-well seismic and its application[J]. Lithologic Reservoirs (in Chinese), 25(5): 94–99. |
| [] | Gao J H, Man W S, Chen S M, et al. 2004. Recognition of signals from colored noise background in generalized S transformation domain[J]. Chinese Journal of Geophysics(in Chinese), 47(5): 869–875. DOI:10.1002/cjg2.v47.5 |
| [] | Han W G, Qu Y M, Hu L X, et al. 2015. Method of suppressing harmonic noise based on reflection coefficient and ground force signal[J]. Progress in Geophysics (in Chinese), 30(6): 2647–2659. DOI:10.6038/pg20150624 |
| [] | Hu J T, Wang H Z, Wang Z L. 2013. Broadband seismic data with vibroseis harmonics[J]. Oil Geophysical Prospecting (in Chinese), 48(6): 876–881. |
| [] | Huang J P, Zhou X F, Guo J, et al. 2012. Method of harmonic noise elimination in slip sweep data[J]. Journal of China University of Petroleum (in Chinese), 36(2): 81–85. |
| [] | Jiang F X, Yang S H, Cheng Y H, et al. 2006. A study on microseismic monitoring of rock burst in coal mine[J]. Chinese Journal of Geophysics (in Chinese), 49(5): 1511–1516. DOI:10.3321/j.issn:0001-5733.2006.05.032 |
| [] | Li M J, Wu G C, Chen H. 2016. Frequency domain scalar wave forward modeling method based on an average-derivative method[J]. Progress in Geophysics (in Chinese), 31(6): 2564–2573. DOI:10.6038/pg20160628 |
| [] | Li Q, Li Q C, Feng H. 2010. Application of mixed-language programming of C# and FORTRAN to seismic velocity analysis[J]. Progress in Geophysics (in Chinese), 25(4): 1503–1507. DOI:10.3969/j.issn.1004-2903.2010.04.044 |
| [] | Li Z C, Guo Z B, Tian K. 2014. Least-squares reverse time migration in visco-acoustic medium[J]. Chinese Journal of Geophysics (in Chinese), 57(1): 214–228. DOI:10.6038/cjg20140118 |
| [] | Lin J, Luo Y, Liu Y W, et al. 2014. Harmonic interference suppression on vibroseis slip sweep data[J]. Oil Geophysical Prospecting (in Chinese), 49(5): 852–856. |
| [] | Liu S W, Chen Z S. 1998. The suppression of harmonic interference of vibroseis[J]. Oil Geophysical Prospecting (in Chinese), 33(S1): 25–30. |
| [] | Martin F D, Munoz P A. 2010. Deharmonics:A method for harmonic noise removal on vibroseis data[C].//72nd EAGE Conference and Exhibition Incorporating SPE EUROPEC 2010. SPE. |
| [] | Pan D F, Wang B, Zhou Z Q. 2009. Research on mixed programming technology of Matlab and C/C+ +[J]. Computer Engineering and Design (in Chinese), 30(2): 465–468. |
| [] | Qu Y M, Li Z C, Huang J P, et al. 2016. Vibroseis harmonics suppression based on adaptive matched predictive filtering[J]. Oil Geophysical Prospecting (in Chinese), 51(6): 1075–1083. |
| [] | Ras P, Daly M, Baeten G. 1999. Harmonic distortion in slip sweep records[C].//1999 SEG Annual Meeting. Houston, Texas:Expanded Abstracts, 609-612. |
| [] | Ren H, Zhou Z H, Zhang C C. 2007. Mixed-language compiling between Fortran and C/C+ +[J]. Computer Engineering and Design (in Chinese), 28(17): 4096–4098, 4111. |
| [] | Sicking C, Fleure T, Nelan S, et al. 2009. Slip sweep harmonic noise rejection on correlated shot data[C].//2009 SEG Annual Meeting. Houston, Texas:Expanded Abstracts. |
| [] | Tong X Q, Lin J, Jiang T, et al. 2012. Summary of development of land vibrator[J]. Progress in Geophysics (in Chinese), 27(5): 1912–1921. DOI:10.6038/j.issn.1004-2903.2012.05.011 |
| [] | Wang H L, Chen J H. 2010. Studies on geophysical software integration environment[J]. Oil Geophysical Prospecting (in Chinese), 45(2): 299–305. |
| [] | Wang P. 2014. Three-dimensional visualization of colliery TEM data base on MATLAB[J]. Progress in Geophysics (in Chinese), 29(3): 1277–1283. DOI:10.6038/pg20140338 |
| [] | Wang Z R, Chen Z B, Zhang L H, et al. 2006. Numerical simulation for seismic response to nonlinear sweeping signal from a vibrator[J]. Progress in Geophysics (in Chinese), 21(3): 759–761. DOI:10.3969/j.issn.1004-2903.2006.03.011 |
| [] | Wang Z R, Gao J, Lin J. 2010. Matched sweep method for seismic vibrators[J]. Chinese Journal of Geophysics (in Chinese), 53(11): 2754–2759. DOI:10.3969/j.issn.0001-5733.2010.11.024 |
| [] | Wei F J, Xu W X. 2013. Applied integration framework technology for seismic acquisition software[J]. Oil Geophysical Prospecting (in Chinese), 48(5): 809–815. |
| [] | Wei Z H, Phillips T F. 2010. Harmonic distortion reduction on seismic vibrators[J]. The Leading Edge, 29(3): 256–261. DOI:10.1190/1.3353719 |
| [] | Wu G C, Qin H X. 2014. Elastic reverse time migration in isotropic medium based on random boundary[J]. Progress in Geophysics (in Chinese), 29(4): 1815–1821. DOI:10.6038/pg20140444 |
| [] | Wu G C, Wang H Z. 2005. Analysis of numerical dispersion in wave-field simulation[J]. Progress in Geophysics (in Chinese), 20(1): 58–65. DOI:10.3969/j.issn.1004-2903.2005.01.012 |
| [] | Xu W X, Duan W X, Yang J, et al. 2016. Technology discussion on the integration and transplantation of seismic data acquisition engineering software[J]. Geophysical Prospecting for Petroleum (in Chinese), 55(3): 350–356, 387. |
| [] | Yang W Y, Wang X W, Yong X S, et al. 2013. The review of seismic full waveform inversion method[J]. Progress in Geophysics (in Chinese), 28(2): 766–776. DOI:10.6038/pg20130225 |
| [] | Zhang J H, Wang W M, Zhao L F, et al. 2007. Modeling 3-D scalar waves using the Fourier finite-difference method[J]. Chinese Journal of Geophysics (in Chinese), 50(6): 1854–1862. DOI:10.3321/j.issn:0001-5733.2007.06.028 |
| [] | 曹务祥. 2004. 组合扫描压制谐波畸变[J]. 石油地球物理勘探, 39(6): 711–715. |
| [] | 刁瑞, 冯玉苹. 2012. 利用储层的频谱吸收特性识别油气藏[J]. 石油地球物理勘探, 47(5): 766–772. |
| [] | 刁瑞, 尚新民, 芮拥军, 等. 2013. 时频域谱模拟反褶积方法研究[J]. 岩性油气藏, 25(1): 116–121. |
| [] | 刁瑞, 吴国忱, 尚新民, 等. 2017. 三维地震与地面微地震联合校正方法[J]. 地球物理学报, 60(1): 283–292. DOI:10.6038/cjg20170123 |
| [] | 冯玉苹, 刁瑞, 徐维秀, 等. 2013. 井间地震纵横波分离方法研究及应用[J]. 岩性油气藏, 25(5): 94–99. |
| [] | 高静怀, 满蔚仕, 陈树民. 2004. 广义S变换域有色噪声与信号识别方法[J]. 地球物理学报, 47(5): 869–875. |
| [] | 韩文功, 曲英铭, 胡立新, 等. 2015. 基于反射系数和地面力信号的谐波干扰消除法[J]. 地球物理学进展, 30(6): 2647–2659. DOI:10.6038/pg20150624 |
| [] | 胡江涛, 王华忠, 王兆磊. 2013. 利用可控震源谐波提高地震数据分辨率[J]. 石油地球物理勘探, 48(6): 876–881. |
| [] | 黄建平, 周学锋, 郭军, 等. 2012. 滑动扫描记录中压制谐波干扰方法[J]. 中国石油大学学报(自然科学版), 36(2): 81–85. |
| [] | 姜福兴, 杨淑华, 成云海, 等. 2006. 煤矿冲击地压的微地震监测研究[J]. 地球物理学报, 49(5): 1511–1516. DOI:10.3321/j.issn:0001-5733.2006.05.032 |
| [] | 李美娟, 吴国忱, 陈浩. 2016. 基于平均导数的标量波频率域正演数值模拟方法[J]. 地球物理学进展, 31(6): 2564–2573. DOI:10.6038/pg20160628 |
| [] | 李勤, 李庆春, 冯宏. 2010. C#和FORTRAN混合编程在地震速度分析中的应用[J]. 地球物理学进展, 25(4): 1503–1507. DOI:10.3969/j.issn.1004-2903.2010.04.044 |
| [] | 李振春, 郭振波, 田坤. 2014. 黏声介质最小平方逆时偏移[J]. 地球物理学报, 57(1): 214–228. DOI:10.6038/cjg20140118 |
| [] | 林娟, 罗勇, 刘宜文, 等. 2014. 可控震源滑动扫描谐波干扰压制方法[J]. 石油地球物理勘探, 49(5): 852–856. |
| [] | 刘世文, 陈振声. 1998. 可控震源谐振干扰的压制[J]. 石油地球物理勘探, 33(S1): 25–30. |
| [] | 潘大夫, 汪渤, 周志强. 2009. Matlab与C/C++混合编程技术研究[J]. 计算机工程与设计, 30(2): 465–468. |
| [] | 曲英铭, 李振春, 黄建平, 等. 2016. 自适应匹配预测滤波压制可控震源谐波[J]. 石油地球物理勘探, 51(6): 1075–1083. |
| [] | 任慧, 周振红, 张成才. 2007. Fortran与C/C++的混合编译[J]. 计算机工程与设计, 28(17): 4096–4098, 4111. DOI:10.3969/j.issn.1000-7024.2007.17.011 |
| [] | 佟训乾, 林君, 姜弢, 等. 2012. 陆地可控震源发展综述[J]. 地球物理学进展, 27(5): 1912–1921. DOI:10.6038/j.issn.1004-2903.2012.05.011 |
| [] | 王宏琳, 陈继红. 2010. 地球物理软件集成环境研究[J]. 石油地球物理勘探, 45(2): 299–305. |
| [] | 王鹏. 2014. 基于MATLAB的煤矿TEM数据体三维可视化技术[J]. 地球物理学进展, 29(3): 1277–1283. DOI:10.6038/pg20140338 |
| [] | 王忠仁, 陈祖斌, 张林行, 等. 2006. 可控震源非线性扫描地震响应的数值模拟[J]. 地球物理学进展, 21(3): 759–761. DOI:10.3969/j.issn.1004-2903.2006.03.011 |
| [] | 王忠仁, 高健, 林君. 2010. 可控震源匹配扫描方法研究[J]. 地球物理学报, 53(11): 2754–2759. DOI:10.3969/j.issn.0001-5733.2010.11.024 |
| [] | 魏福吉, 徐维秀. 2013. 面向地震数据采集工程软件的应用集成框架技术[J]. 石油地球物理勘探, 48(5): 809–815. |
| [] | 吴国忱, 秦海旭. 2014. 各向同性弹性波随机边界逆时偏移[J]. 地球物理学进展, 29(4): 1815–1821. DOI:10.6038/pg20140444 |
| [] | 吴国忱, 王华忠. 2005. 波场模拟中的数值频散分析与校正策略[J]. 地球物理学进展, 20(1): 58–65. DOI:10.3969/j.issn.1004-2903.2005.01.012 |
| [] | 徐维秀, 段卫星, 杨晶, 等. 2016. 地震勘探采集工程软件集成与移植技术探讨[J]. 石油物探, 55(3): 350–356, 387. |
| [] | 杨午阳, 王西文, 雍学善, 等. 2013. 地震全波形反演方法研究综述[J]. 地球物理学进展, 28(2): 766–776. DOI:10.6038/pg20130225 |
| [] | 张金海, 王卫民, 赵连锋, 等. 2007. 傅里叶有限差分法三维波动方程正演模拟[J]. 地球物理学报, 50(6): 1854–1862. DOI:10.3321/j.issn:0001-5733.2007.06.028 |