泌阳凹陷复杂断裂带地震勘探采集处理方法研究与应用
|
摘要
王集—新庄地区位于泌阳凹陷北部斜坡的中东部的王集—新庄鼻状构造带上,具有良好的构造背景,是油气运移的有利指向区。该区主要发育一系列北东向断层和少量北西向断层,断层的相互交切并与构造背景配置形成了一系列断鼻、断块,为油气的聚集提供良好的圈闭条件。北东方向的候庄三角洲砂体和北西方向的王集三角洲砂体在该区广泛分布,具有良好的储层条件,总体来说该区是形成油气聚集的有利场所。
目前,在该区已取得了较大发现,探明石油地质储量近4000万吨,主要分布在核三上段,核三下段,目前所发现的储量比例较小,并且主要分布在王集西部地区,在王集东部和新庄地区基本没发现核三下段地质储量。综合分析研究认为,该区核三下段的成藏条件与核三上段类似,部分井在核三下段也发现了不同级别的油气显示,还有个别井见到了工业油流,认为核三下段也是形成油气聚集的有利场所,且核三下段分布面积广,在斜坡地带核三上段大部分被剥蚀的地区,核三下段仍有较好保存。
但是,王集-新庄地区目前的地震资料仅能满足浅层勘探的需要,由于中深层该区地震资料品质较差,难以满足进一步勘探的需要。由于受当时技术条件限制,且以前三维地震勘探以核三上地震勘探目的层,已有的地震资料从核三段Ⅲ砂组以下地震反射同相轴连续性变差,信噪比低,反射品质较差。受地表地震地质条件和深部地震地质条件横向不均匀性的影响,造成了资料品质的横向不均匀性。小断块群成像不清,横向和垂向分辨率较低,层位追踪及一些小断层的解释仍有一定的模糊性、多解性,很难对核三下段构造特征进行准确、精细的描述。因此,针对王集-新庄地区部署了重新采集任务,针对泌阳凹陷核三下段地震反射信号弱、信噪比和分辨率低的问题,开展技术攻关,提高中深层地震资料的信噪比、连续性和对比精度,提高小断块群成像精度,搞清复杂断块群的构造断裂特征。
本文针对复杂断裂带地震反射信号弱、信噪比和分辨率低的问题,基于复杂断裂带的地震地质特点,开展高精度地震采集,采用基于波动方程正演模拟和照明分析技术进行野外观测系统的设计和优化,达到对复杂断裂带地下地质目标的最佳观测,最大限度获取完整丰富的波场信息。在处理过程中,充分利用高精度采集的优势,开展等效模型层析静校正方法研究,精确计算近地表模型和静校正量,提高分辨率和小断层成像精度；采用正演与共轭反演算子叠前深度偏移算法和速度分析方法,使断层准确归位,提高复杂断裂带的成像精度,获得准确的构造信息和储层信息。形成了一套适应复杂断裂带油藏勘探的采集和处理方法,为合理落实复杂断裂带的构造面貌,提高钻探效果提供技术支撑。
本文的主要工作及结论如下：
(1)本研究利用高精度的交错网格高阶有限差分方法,利用双程波方程模拟了泌阳凹陷复杂断裂带模型中地震波的传播,进而对王集断鼻构造带的多级断裂构造进行了地震波照明分析,实现野外观测系统的设计和优化。确定最佳分辨率面元的大小、地面最佳采集位置、最大最小炮检距、排列长度、排列滚动等,直观显示能量在地下介质中的分布情况,分析地震波传播规律和地层能量响应,达到对复杂断裂带地下地质目标的最佳观测,最大限度获取完整丰富的波场信息。在有效目的层范围内本次地震采集的观测系统在偏移距和方位角的分布上基本接近于高精度宽方位采集,可对地下进行全方位照明,既有利于浅层不整合面的完整显现和追踪,也有利于复杂断裂的成像,有利于提高地震资料的信噪比,满足构造解释和岩性勘探需要。
(2)根据模型照明分析技术,综合判断后采纳了方案1。在南部采用24L8S136R204F观测系统,北部采用24L8S144R108F观测系统,能保证220m以下的目的层信噪比,并且炮密度小成本较低,炮检距集中度较好,施工组织方便。根据方案1的偏移距分布显示,近偏移距和远偏移距覆盖都比较均匀,,200米以内偏移距基本全部覆盖。近偏移距的均匀分布保证了剖面无缺口,有利于超浅层及不整合面成像。
(3)在泌阳凹复杂陷断裂带,层析静校正比折射静校正和高程静校正有更好的效果。它可使弯曲和连续性差的同相轴形态变得更加自然顺畅,信噪比增高,层析静校正后浅部同相轴细节丰富,频率增高,说明达到了同相叠加的效果。
(4)引入了正演与共轭反演算子,实现了炮域波动方程延拓法叠前深偏移。该算法稳定,适用于复杂速度场,并且计算精度高。在泌阳凹陷复杂断裂带开展提高炮域波动方程延拓法叠前深度偏移的成像精度和计算效率的方法研究。偏移结果信号保真度较好,地震反射层特征明显,各个主要层段波组特征清楚,地质现象丰富,利于地层追踪对比,小断层断面成像清楚,断点清晰,偏移归位合理,断裂系统成像得到明显改善,构造形态改善明显,同相轴自然,信噪比和连续性都得到提高。
(5)在浮动基准面处理和精细选择切除函数,有效利用高精度采集的资料,保留浅层有效反射信息,浅层不整合面保留清楚；
(6)处理结果分辨率得到有效提高。资料频带较宽,有效反射信号频带基本分布在5～100Hz,核二和核三上段(剖面500ms以上)主频达到45-50Hz,核三下段(剖面700ms～1700ms)主频达到30～40Hz。基底成像准确,形态清楚,有利于区域构造研究、沉积演化研究和应力分析。
本论文创新点如下：
(1)引入了正演与共轭反演算子,实现了炮域波动方程延拓法叠前深偏移。该算法稳定,适用于复杂速度场,并且计算精度高。在泌阳凹陷复杂断裂带开展提高炮域波动方程延拓法叠前深度偏移的成像精度和计算效率的方法研究。偏移结果信号保真度较好,地震反射层特征明显,各个主要层段波组特征清楚,地质现象丰富,利于地层追踪对比,小断层断面成像清楚,断点清晰,偏移归位合理,断裂系统成像得到明显改善,构造形态改善明显,同相轴自然,信噪比和连续性都得到提高。
(2)首次将波动方程正演和照明技术应用于泌阳凹陷复杂断裂带观测系统设计,并取得了较好的应用效果。波动方程正演和照明技术解决了泌阳凹陷复杂断裂带数据采集要求的高精度、宽方位、偏移距均匀、覆盖次数均匀等问题,实现了对地下地质目标的均匀和全息观测。
Wangji-Xinzhuang area, located in the Wangji-Xinzhuang nose structure zone, which is in the central and eastern part of the northern slopes of the Biyang Depression, has a good tectonic setting, and is a favorable directional zone of hydrocarbon migration. The area mainly consists of a series of NE-trending faults and a small amount of NW-trending faults. Configuration of intersection faults and tectonic setting lead to the formation of a series of fault blocks and faulted noses, providing trap conditions for accumulation of oil and gas. The sandstone body of Houzhuang delta in the northeast area and of Wangji delta in the northwest area is distribute widely in this area. So the area is favorable for the formation of oil and gas accumulation.
Currently, a great discovery has been made in this area that proven oil reserves here is nearly 40 million tons, which mainly distribute in upper member of Eh3. While small proportion of reservoir has been found in underpart member of Eh3, which mainly distribute in western area of Wangji, and none of the geologic reserve in underpart member of Eh3 was found in eastern Wangji and Xinzhuang area. Comprehensive analysis results indicate that reservoir forming condition in upper member of Eh3 is similar to that in underpart member of Eh3. In some wells have been displayed different level evidences of oil and gas in underpart member of Eh3, with a few wells obtaining commercial oil. Underpart member of Eh3 is regarded as the favorable place for hydrocarbon migration and accumulation. The underpart member of Eh3 distribute widely, which was still preserved well after large-scale erosions in upper member of Eh3.
However, the current seismic data in Wangji-Xinzhuang region is only able to meet the needs of the shallow seismic exploration and it is difficult to meet the needs of future exploration because of the poor seismic data quality. Due to technological constraints at that time, and in the previous seismic three-dimensional exploration, the upper member of the Eh3 being the exploration target, the existing seismic data from the formation below the member H3 III Sand shows the worse continuity of seismic reflection, low SNR, and poor quality reflection. The imaging of small fault blocks is unclear with low lateral resolution and vertical resolution. Horizon tracking and interpretation of small faults is going with a certain degree of ambiguity and multi-resolution. It is difficult to describe exactly and accurately about the structure of underpart member of Eh3.Therefore, for Wangji-Xinzhuang region, there is a deployment of new acquisition task. Be aimed at the problems such as weak earthquake reflection, low SNR and low-resolution of underpart member of Eh3 in Biyang depression, it is necessary to strive to make technological breakthrough to increase the SNR, continuity and Comparison of accuracy of the seismic data in underpart member of Eh3 in Biyang depression, so as to increase the imaging precision of small fault blocks and beating out the fracture feature of the structure of the small fault blocks.
This topic launched a feasibility studies about technical design in acquisition area, we analyzed the main technology and construction difficulty in the survey areas, and systematically analyzed the previous information according to geological tasks and technical requirements. In the base of Reconnaissance of the work area, data collection and analysis of the original three-dimensional data, we proposed feasibility studies and geometry design. Based on fine structure analysis, fault isolation and recognition of translocation system, the pre-stack migration is applied to the high-precision 3D seismic data in this area. Preserving shallow unconformable surface as possible, we increased the SNR of the seismic data and imaging precision of the complex fault, to provide accurate information for structural and lithology interpretation.
The main contribution and conclusion are as follows:
(1)High precession staggering grid high-order finite-difference algorithm was used in this research. Seismic wave propagation of complex fault belt in Biyang Depression was simulated with two-way wave equation. Then, seismic illumination analysis of multilevel tectonic system of Wangji fault nose structure was complicated, and the geometry system was designed and optimized. The bin size with optimum resolution, the best acquisition position, the maximum and minimum offsets, the arrangement length, and the arrangement rolling are all determined. Thus, it is possible to show the energy distribution in the underground medium, analyze the wave propagating pattern and stratum energy responds, and optimally observe the geological target in complex fault belt, in order to obtain the wave field information as complete and abundant as possible. In the effective range of target strata, this acquisition observation system was similar to high-density with wide-azimuth acquisition in terms of distribution of the offset and the azimuth. Therefore, it is possible to illuminate the underground medium omnidirectionally. In another word, it is favorable to not only the tracing of the shallow layers, but also the imagination of complex fault. Thus, signal to noise ratio of the seismic data was increased, and the tectonic interpretation and lithological exploration can be carried out.
(2)After comprehensive judgment based on illumination analysis, the first scheme is commended.24L8S136R204F observation system is used in the south area, and 24L8S144R108F observation system is used in the north area. The signal to noise ratios of target strata below-220m could be accepted, density of the resources were small, the cost was relatively low, the concentration ratio were high, and it was convenient to operate. According to the offset distribution in the first scheme, both the near offset and the far offset were uniformly covered, and the offsets within 200m were almost full-covered. The near offset uniformly coverage was beneficial to the imagination of ultra-shallow layers and unconformities.
(3)In the complex fault belt of Biyang, tomographic static correction is more effective than refraction static correction and elevation static correction. The bending and discontinuous events became more natural and fluent, and the signal to noise ratio were increased. After tomographic static correction the shallow events details were abundant, and the frequency were improved. Thus, the effect of in-phase stacking was achieved.
(4)The forward operator and conjugate inverse was introduced in this paper. Pre-stack depth migration with shot domain wave equation continuation was completed. This algorithm is stable, and is available for complex velocity field, and the calculation precision is high. The methods were studied, to improve the imaging precision and calculating efficiency of the pre-stack depth migration with shot domain wave equation continuation. Signal fidelity of the migration results were relatively high, seismic reflection characters were obvious, and wave group characteristics of the main layer sections were clear. The geological phenomenon was abundant, and this made the strata tracing simple. The section imagination of small faults was clear, the breakpoint was determined, and the migration was reasonable. In the end, imagination of the fault system was improved obviously, the structural configuration was improved, the events were natural, and both the signal to noise ratio and the continuation were increased.
(5)Chopping functions were processed and optimized on floating datum, and high-density acquisition data were used effectively. The effective reflection information of shallow layers was maintained, and the shallow unconformities were clear.
(6)The resolution of the results was improved effectively. The frequency band were relatively wide, and the effective reflection frequencies were 5～100Hz. Dominant frequencies of the upper second and third member of Hetaoyuan formation (500ms and the above of the section) reached 45～50Hz. Dominant frequencies of the lower third member of Hetaoyuan formation (700ms～1700ms of the section) reached 30～40Hz. Imagination of the basement was accurate, and the morphology was clear, which were available to do research on regional tectonics, sedimentary evolution, and stress analysis.
The innovative points of this paper are as follows:
(1)The forward operator and conjugate inverse was introduced in this paper. Pre-stack depth migration with shot domain wave equation continuation was completed. This algorithm is stable, and is available for complex velocity field, and the calculation precision is high. The methods were studied, to improve the imaging precision and calculating efficiency of the pre-stack depth migration with shot domain wave equation continuation. Signal fidelity of the migration results were relatively high, seismic reflection characters were obvious, and wave group characteristics of the main layer sections were clear. The geological phenomenon was abundant, and this made the strata tracing simple. The section imagination of small faults was clear, the breakpoint was determined, and the migration was reasonable. In the end, imagination of the fault system was improved obviously, the structural configuration was improved, the events were natural, and both the signal to noise ratio and the continuation were increased.
It is the first time that the wave equation modeling and seismic illumination method applied in the geometry design, the result is very well. The technology meets the requirement of high density, wide azimuth, equal offset and fold number in acquisition parameters. It makes the equal and all azimuth observation of the geological target become reality.
目前,在该区已取得了较大发现,探明石油地质储量近4000万吨,主要分布在核三上段,核三下段,目前所发现的储量比例较小,并且主要分布在王集西部地区,在王集东部和新庄地区基本没发现核三下段地质储量。综合分析研究认为,该区核三下段的成藏条件与核三上段类似,部分井在核三下段也发现了不同级别的油气显示,还有个别井见到了工业油流,认为核三下段也是形成油气聚集的有利场所,且核三下段分布面积广,在斜坡地带核三上段大部分被剥蚀的地区,核三下段仍有较好保存。
但是,王集-新庄地区目前的地震资料仅能满足浅层勘探的需要,由于中深层该区地震资料品质较差,难以满足进一步勘探的需要。由于受当时技术条件限制,且以前三维地震勘探以核三上地震勘探目的层,已有的地震资料从核三段Ⅲ砂组以下地震反射同相轴连续性变差,信噪比低,反射品质较差。受地表地震地质条件和深部地震地质条件横向不均匀性的影响,造成了资料品质的横向不均匀性。小断块群成像不清,横向和垂向分辨率较低,层位追踪及一些小断层的解释仍有一定的模糊性、多解性,很难对核三下段构造特征进行准确、精细的描述。因此,针对王集-新庄地区部署了重新采集任务,针对泌阳凹陷核三下段地震反射信号弱、信噪比和分辨率低的问题,开展技术攻关,提高中深层地震资料的信噪比、连续性和对比精度,提高小断块群成像精度,搞清复杂断块群的构造断裂特征。
本文针对复杂断裂带地震反射信号弱、信噪比和分辨率低的问题,基于复杂断裂带的地震地质特点,开展高精度地震采集,采用基于波动方程正演模拟和照明分析技术进行野外观测系统的设计和优化,达到对复杂断裂带地下地质目标的最佳观测,最大限度获取完整丰富的波场信息。在处理过程中,充分利用高精度采集的优势,开展等效模型层析静校正方法研究,精确计算近地表模型和静校正量,提高分辨率和小断层成像精度；采用正演与共轭反演算子叠前深度偏移算法和速度分析方法,使断层准确归位,提高复杂断裂带的成像精度,获得准确的构造信息和储层信息。形成了一套适应复杂断裂带油藏勘探的采集和处理方法,为合理落实复杂断裂带的构造面貌,提高钻探效果提供技术支撑。
本文的主要工作及结论如下：
(1)本研究利用高精度的交错网格高阶有限差分方法,利用双程波方程模拟了泌阳凹陷复杂断裂带模型中地震波的传播,进而对王集断鼻构造带的多级断裂构造进行了地震波照明分析,实现野外观测系统的设计和优化。确定最佳分辨率面元的大小、地面最佳采集位置、最大最小炮检距、排列长度、排列滚动等,直观显示能量在地下介质中的分布情况,分析地震波传播规律和地层能量响应,达到对复杂断裂带地下地质目标的最佳观测,最大限度获取完整丰富的波场信息。在有效目的层范围内本次地震采集的观测系统在偏移距和方位角的分布上基本接近于高精度宽方位采集,可对地下进行全方位照明,既有利于浅层不整合面的完整显现和追踪,也有利于复杂断裂的成像,有利于提高地震资料的信噪比,满足构造解释和岩性勘探需要。
(2)根据模型照明分析技术,综合判断后采纳了方案1。在南部采用24L8S136R204F观测系统,北部采用24L8S144R108F观测系统,能保证220m以下的目的层信噪比,并且炮密度小成本较低,炮检距集中度较好,施工组织方便。根据方案1的偏移距分布显示,近偏移距和远偏移距覆盖都比较均匀,,200米以内偏移距基本全部覆盖。近偏移距的均匀分布保证了剖面无缺口,有利于超浅层及不整合面成像。
(3)在泌阳凹复杂陷断裂带,层析静校正比折射静校正和高程静校正有更好的效果。它可使弯曲和连续性差的同相轴形态变得更加自然顺畅,信噪比增高,层析静校正后浅部同相轴细节丰富,频率增高,说明达到了同相叠加的效果。
(4)引入了正演与共轭反演算子,实现了炮域波动方程延拓法叠前深偏移。该算法稳定,适用于复杂速度场,并且计算精度高。在泌阳凹陷复杂断裂带开展提高炮域波动方程延拓法叠前深度偏移的成像精度和计算效率的方法研究。偏移结果信号保真度较好,地震反射层特征明显,各个主要层段波组特征清楚,地质现象丰富,利于地层追踪对比,小断层断面成像清楚,断点清晰,偏移归位合理,断裂系统成像得到明显改善,构造形态改善明显,同相轴自然,信噪比和连续性都得到提高。
(5)在浮动基准面处理和精细选择切除函数,有效利用高精度采集的资料,保留浅层有效反射信息,浅层不整合面保留清楚；
(6)处理结果分辨率得到有效提高。资料频带较宽,有效反射信号频带基本分布在5～100Hz,核二和核三上段(剖面500ms以上)主频达到45-50Hz,核三下段(剖面700ms～1700ms)主频达到30～40Hz。基底成像准确,形态清楚,有利于区域构造研究、沉积演化研究和应力分析。
本论文创新点如下：
(1)引入了正演与共轭反演算子,实现了炮域波动方程延拓法叠前深偏移。该算法稳定,适用于复杂速度场,并且计算精度高。在泌阳凹陷复杂断裂带开展提高炮域波动方程延拓法叠前深度偏移的成像精度和计算效率的方法研究。偏移结果信号保真度较好,地震反射层特征明显,各个主要层段波组特征清楚,地质现象丰富,利于地层追踪对比,小断层断面成像清楚,断点清晰,偏移归位合理,断裂系统成像得到明显改善,构造形态改善明显,同相轴自然,信噪比和连续性都得到提高。
(2)首次将波动方程正演和照明技术应用于泌阳凹陷复杂断裂带观测系统设计,并取得了较好的应用效果。波动方程正演和照明技术解决了泌阳凹陷复杂断裂带数据采集要求的高精度、宽方位、偏移距均匀、覆盖次数均匀等问题,实现了对地下地质目标的均匀和全息观测。
Wangji-Xinzhuang area, located in the Wangji-Xinzhuang nose structure zone, which is in the central and eastern part of the northern slopes of the Biyang Depression, has a good tectonic setting, and is a favorable directional zone of hydrocarbon migration. The area mainly consists of a series of NE-trending faults and a small amount of NW-trending faults. Configuration of intersection faults and tectonic setting lead to the formation of a series of fault blocks and faulted noses, providing trap conditions for accumulation of oil and gas. The sandstone body of Houzhuang delta in the northeast area and of Wangji delta in the northwest area is distribute widely in this area. So the area is favorable for the formation of oil and gas accumulation.
Currently, a great discovery has been made in this area that proven oil reserves here is nearly 40 million tons, which mainly distribute in upper member of Eh3. While small proportion of reservoir has been found in underpart member of Eh3, which mainly distribute in western area of Wangji, and none of the geologic reserve in underpart member of Eh3 was found in eastern Wangji and Xinzhuang area. Comprehensive analysis results indicate that reservoir forming condition in upper member of Eh3 is similar to that in underpart member of Eh3. In some wells have been displayed different level evidences of oil and gas in underpart member of Eh3, with a few wells obtaining commercial oil. Underpart member of Eh3 is regarded as the favorable place for hydrocarbon migration and accumulation. The underpart member of Eh3 distribute widely, which was still preserved well after large-scale erosions in upper member of Eh3.
However, the current seismic data in Wangji-Xinzhuang region is only able to meet the needs of the shallow seismic exploration and it is difficult to meet the needs of future exploration because of the poor seismic data quality. Due to technological constraints at that time, and in the previous seismic three-dimensional exploration, the upper member of the Eh3 being the exploration target, the existing seismic data from the formation below the member H3 III Sand shows the worse continuity of seismic reflection, low SNR, and poor quality reflection. The imaging of small fault blocks is unclear with low lateral resolution and vertical resolution. Horizon tracking and interpretation of small faults is going with a certain degree of ambiguity and multi-resolution. It is difficult to describe exactly and accurately about the structure of underpart member of Eh3.Therefore, for Wangji-Xinzhuang region, there is a deployment of new acquisition task. Be aimed at the problems such as weak earthquake reflection, low SNR and low-resolution of underpart member of Eh3 in Biyang depression, it is necessary to strive to make technological breakthrough to increase the SNR, continuity and Comparison of accuracy of the seismic data in underpart member of Eh3 in Biyang depression, so as to increase the imaging precision of small fault blocks and beating out the fracture feature of the structure of the small fault blocks.
This topic launched a feasibility studies about technical design in acquisition area, we analyzed the main technology and construction difficulty in the survey areas, and systematically analyzed the previous information according to geological tasks and technical requirements. In the base of Reconnaissance of the work area, data collection and analysis of the original three-dimensional data, we proposed feasibility studies and geometry design. Based on fine structure analysis, fault isolation and recognition of translocation system, the pre-stack migration is applied to the high-precision 3D seismic data in this area. Preserving shallow unconformable surface as possible, we increased the SNR of the seismic data and imaging precision of the complex fault, to provide accurate information for structural and lithology interpretation.
The main contribution and conclusion are as follows:
(1)High precession staggering grid high-order finite-difference algorithm was used in this research. Seismic wave propagation of complex fault belt in Biyang Depression was simulated with two-way wave equation. Then, seismic illumination analysis of multilevel tectonic system of Wangji fault nose structure was complicated, and the geometry system was designed and optimized. The bin size with optimum resolution, the best acquisition position, the maximum and minimum offsets, the arrangement length, and the arrangement rolling are all determined. Thus, it is possible to show the energy distribution in the underground medium, analyze the wave propagating pattern and stratum energy responds, and optimally observe the geological target in complex fault belt, in order to obtain the wave field information as complete and abundant as possible. In the effective range of target strata, this acquisition observation system was similar to high-density with wide-azimuth acquisition in terms of distribution of the offset and the azimuth. Therefore, it is possible to illuminate the underground medium omnidirectionally. In another word, it is favorable to not only the tracing of the shallow layers, but also the imagination of complex fault. Thus, signal to noise ratio of the seismic data was increased, and the tectonic interpretation and lithological exploration can be carried out.
(2)After comprehensive judgment based on illumination analysis, the first scheme is commended.24L8S136R204F observation system is used in the south area, and 24L8S144R108F observation system is used in the north area. The signal to noise ratios of target strata below-220m could be accepted, density of the resources were small, the cost was relatively low, the concentration ratio were high, and it was convenient to operate. According to the offset distribution in the first scheme, both the near offset and the far offset were uniformly covered, and the offsets within 200m were almost full-covered. The near offset uniformly coverage was beneficial to the imagination of ultra-shallow layers and unconformities.
(3)In the complex fault belt of Biyang, tomographic static correction is more effective than refraction static correction and elevation static correction. The bending and discontinuous events became more natural and fluent, and the signal to noise ratio were increased. After tomographic static correction the shallow events details were abundant, and the frequency were improved. Thus, the effect of in-phase stacking was achieved.
(4)The forward operator and conjugate inverse was introduced in this paper. Pre-stack depth migration with shot domain wave equation continuation was completed. This algorithm is stable, and is available for complex velocity field, and the calculation precision is high. The methods were studied, to improve the imaging precision and calculating efficiency of the pre-stack depth migration with shot domain wave equation continuation. Signal fidelity of the migration results were relatively high, seismic reflection characters were obvious, and wave group characteristics of the main layer sections were clear. The geological phenomenon was abundant, and this made the strata tracing simple. The section imagination of small faults was clear, the breakpoint was determined, and the migration was reasonable. In the end, imagination of the fault system was improved obviously, the structural configuration was improved, the events were natural, and both the signal to noise ratio and the continuation were increased.
(5)Chopping functions were processed and optimized on floating datum, and high-density acquisition data were used effectively. The effective reflection information of shallow layers was maintained, and the shallow unconformities were clear.
(6)The resolution of the results was improved effectively. The frequency band were relatively wide, and the effective reflection frequencies were 5～100Hz. Dominant frequencies of the upper second and third member of Hetaoyuan formation (500ms and the above of the section) reached 45～50Hz. Dominant frequencies of the lower third member of Hetaoyuan formation (700ms～1700ms of the section) reached 30～40Hz. Imagination of the basement was accurate, and the morphology was clear, which were available to do research on regional tectonics, sedimentary evolution, and stress analysis.
The innovative points of this paper are as follows:
(1)The forward operator and conjugate inverse was introduced in this paper. Pre-stack depth migration with shot domain wave equation continuation was completed. This algorithm is stable, and is available for complex velocity field, and the calculation precision is high. The methods were studied, to improve the imaging precision and calculating efficiency of the pre-stack depth migration with shot domain wave equation continuation. Signal fidelity of the migration results were relatively high, seismic reflection characters were obvious, and wave group characteristics of the main layer sections were clear. The geological phenomenon was abundant, and this made the strata tracing simple. The section imagination of small faults was clear, the breakpoint was determined, and the migration was reasonable. In the end, imagination of the fault system was improved obviously, the structural configuration was improved, the events were natural, and both the signal to noise ratio and the continuation were increased.
It is the first time that the wave equation modeling and seismic illumination method applied in the geometry design, the result is very well. The technology meets the requirement of high density, wide azimuth, equal offset and fold number in acquisition parameters. It makes the equal and all azimuth observation of the geological target become reality.
引文
[1]陈生昌,马在田,吴如山.波动方程双程地下方向照明分析,同济大学学报(自然科学版),2007,35(5)：681-684
[2]冯伟,王华忠,吴如山等,面向目标的小束源照明和成像,地球物理学进展,2006,21(3)：802-808
[3]董良国,吴晓丰,唐海忠,等.逆掩推覆构造的地震波照明与观测系统优化,石油物探,2006,45(1)：40-47
[4]Liu Shouwei, Geng Jianhua, Feng Wei. Controlled illumination and seismic acquisition geometry for target-oriented imaging, Applied Geophysics,2005, V1.2(4):230-234
[5]Dong Liangguo, Cai Xiwei, et al. Seismic wave illumination for Yumen overthrust structure and seismic survey optimal design, SEG Annual Meeting,2005,48-51
[6]冯伟,吴如山,马在田.单向传播子的局部平面波分解以及在定向照明分析和目标特征成像中的应用,地球物理学进展,2007,22(1)：171-178
[7]裴正林,波动方程地震定向照明分析,石油地球物理勘探,2008,43(6)：645-651
[8]Luo M, Cao J, Xie X B, et al,2004, Comparison of illumination analysis using one-way and full-wave propagators:74th Annual International Meeting, SEG, Expanded Abstracts, 67-70
[9]吕公河,尹成,周星合,等.基于采集目标的地震照明度的精确模拟,石油地球物理勘探,2006,41(3)：258-261
[10]熊晓军,贺振华,文晓涛,等.基于波动方程的观测系统设计方法研究,中国矿业大学学报,2007,36(1),132-136
[11]赵殿栋,郭建,王咸彬,等.基于模型面向目标的观测系统优化设计技术[J].2006,2(2)：119-122
[12]陈秀梅,王华忠,程玖兵等.地表旋转控制照明叠前深度偏移[J].地球物理学报,2004,47(2)：306-311
[13]陈秀梅,王华忠,马在田.多地质目标地表控制照明叠前深度偏移[J].石油地球物理勘探,2003,38(5)：507-511
[14]吕林,李斐,李振春等.聚焦点控制照明波动方程叠前深度偏移[J].岩性油气藏,2007,19(1)：113-116
[15]林伯香,孙晶梅,徐颖,等.几种常用静校正方法的讨论[J].石油物探,2006,23(4)：368-372.
[16]杨红霞,高林.复杂地表地区地震勘探静校正方法应用研究[J].世界地质,2002,21(3)：294-298.
[17]方伍宝,李满树,孙爱萍,等.基于Born近似的波动方程静校正技术[J].石油物探,2004,43(1)：27-29.
[18]何光明,贺振华,黄德济,等.几种静校正方法的比较研究[J].物探化探计算技术,2006,28(4)：311-314.
[19]王西文,徐凤银,冯云发,等.三湖地区地震资料处理中的静校正方法研究[J].石油地球物理勘探,2007,42(5)：572-580.
[20]张海新,丁伟,刘秀玉,等.永新地区近地表校正方法[J].石油地球勘探,2008,43(增刊2)：63-66.
[21]薛为平,宋阳,刘居文,等.准噶尔盆地石南地区静校正方法应用研究[J].石油物探,2006,45(6)：615-618.
[22]李录明,罗省贤.三维波场延拓复杂表层模型校正方法及应用[J].地球物理学进展,2005,20(4)：1027-1034.
[23]朱生旺.近地表速度菲涅耳带初至层析反演[J].新疆石油地质,2007,28(6)：753-755.
[24]陈立新,王振华.SMODE1层析反演近地表模型静校正[J].天然气工业,2007,27(增刊A)：128-129.
[25]何光明,贺振华黄德济王翠华.非线性层析静校正在川西地区资料处理中的应用[J].石油物探,2006,45(1)：89-92.
[26]徐凌,崔兴福,齐莉,等.复杂近地表层析反演与波场延拓联合基准面校正[J].石油物探,2007,46(3)：227-230.
[27]韩晓丽,杨长春,麻三怀,等.复杂山区初至波层析反演静校正[J].地球物理学进展,2008,23(2)：475-483.
[28]杨兰锁,王世煜,刘财,等.射线层析静校正[J].吉林大学学报(地球科学版),2006,36(专辑)：50-75.
[29]苏贵仕,黄莉莉,赵明秋.伪三维层析反演静校正技术及其应用[J].石油地球物理堪探,2007,42(增刊)：111-114.
[30]时磊,顾汉明,刘洪雷.无射线追踪法层析静校正及其应用[J].工程地球物理学报,2008,5(5)：575-578.
[31]秦亚玲,王彦春,王瑞祥,等.用多项式拟合求取沙漠地区表层速度的静校正方法[J].天然气工业,2009,29(6)：36-39.
[32]赵传雪,王丽,吴靖,等.有限差分波动方程基准面校正方法及其在丘陵地区的应用[J].石油物探,2009,48(5)：505-509.
[33]张印堂.折射层析反演静校正及其在伊朗卡山地区的应用[J].油气地球物理,2005,3(3)：11-13.
[34]林伯香,孙晶梅,刘清林.层析成像低速带速度反演和静校正方法[J].石油物探,2002,41(2)：137-140.
[35]林伯香,肖万富,李博.层析静校正在黄土塬弯宽线资料处理中的应用[J].石油物探,2007,46(4)：418-420.
[36]林伯香,孙晶梅,徐颖,等.复杂地表条件静校正中的3D表层速度层析反演研究[J].石油 物探,2005,22(5)：455-457.
[37]陈世军,张建中.初至波射线层析成像在复杂区静校正中的应用[J].石油物探,2006,45(1)：35-39.
[38]何光明,贺振华,黄德济,等.无射线追踪层析静校技术在川西复杂地区资料处理中的应用[J].成都理工大学学报(自然科学版),2006,33(6)：599-602.
[39]左海,宋小平,魏庚雨.折射和层析静校正量计算方法的综合应用[J].石油地球物理勘探,2007(增刊)：35-39.
[40]宋涛,工秀云,杨世文,等.层析反演静校正技术在永新地区的应用[J].石油天然气学报(江汉石油学院学报),2008,30(5)：228-230.
[41]张继国,刘连升.复杂区初至层析反演静校正[J].石油地球物理勘探,2006,41(4)：383-395.
[42]Cui Xingfu, Li Hongbing, Hu Ying, Liang Hong, and Qi Li. Wavefield continuation datuming using a near surface model[J]. APPLIED GEOPHYSICS,2007,4(2):94-100
[43]冯泽元,李培明,唐海忠,等.利用层析反演技术解决山地复杂区静校正问题[J].石油物探,2005,44(3)：284-287.
[44]马在田.高阶方程偏移的分裂算法[J].地球物理学报,1983,26：377～388.
[45]马在田.论反射地震偏移成像[J].勘探地球物理进展,2002,25(3)：1-5.
[46]陈生昌,马在田.地震成像中的Born近似与Rytov近似比较,天然气工业,2007,27(1),37-40
[47]陈生昌,马在田,Wu Rushan,波动方程深度偏移的局部裂步Fourier传播算子,计算物理,2006,23(5)：604-608
[48]Rietveld W E A, Berkhout A J, Pre-stack depth migration by means of controlled illumination, Geophysics,1994,59(5):801-809
[49]Wu R S, Wang Y, and Gao J H, Beamlet migration based on local perturbation theory,70th Annual meeting, SEG Expanded abstracts.1008-1011
[50]Ye Yueming, Li Zhenchun, et al., Beamlet prestack depth migration and illumination:A test based on the Marmousi model, Applied Geophysics,2006,3(4):203-209
[51]王勇,程金星,庞全康等.弹性波正演技术在地震采集中的应用[J].小型油气藏,2007,12：21-25
[52]程玖兵,王华忠,马在田.频率空间域有限差分法叠前深度偏移[J].地球物理学报,2001,44(3)：389-395
[53]De Hoop, M, Le R—J, and Wu R. S., Generalization of the phase-screen approximation for the scattering of acoustic waves. Wave motion,2000,31:43-70
[54]李录明,罗省贤.相移加校正变范围三维叠前深度偏移[J].成都理工学院学报,1999,26(4)：329-334
[55]张叔伦,孙沛勇.基于平面波合成的傅里叶有限差分叠前深度偏移[J].石油地球物理勘探,1999,34(1)：1-7
[56]李锋,顾汉明.泌阳凹陷南部陡坡带地震资料的炮域波动方程叠前深度偏移[J].地质 科技情报,2007,26(3)：104-108
[57]曹文俊,李振春,韩文功.‘波动方程叠前深度偏移对速度模型的敏感性分析[J].天然气工业,2007,27(增刊A)：120-122
[58]马淑芳,李振春.波动方程叠前深度偏移方法综述[J].勘探地球物理进展,2007,30(3)：153-161
[59]李振春,姚云霞,马在田,等.波动方程共成像点道集偏移速度建模[J].地球物理学报,2003,46(1)：86-93
[60]徐琴,李振春,鲍伟.波动方程偏移中的成像条件研究进展[J].勘探地球物理进展,2008,31(4)：247-252
[61]杨敬磊,李振春,叶月明,等.地震照明叠前深度偏移方法综述[J].地球物理学进展,2008,23(1)：146-152
[62]刘定进,印兴耀.傅里叶有限差分法保幅叠前深度偏移方法[J].地球物理学报,2007,50(1)：268-27
[63]曹文俊,李振春,陈业全.高精度三维叠前深度偏移在DP地区的应用[J].天然气工业,2007,27(增刊A)：123-124
[64]刘定进,印兴耀.共炮检距道集波动方程保幅叠前深度偏移方法[J].地球物理学进展,2007,22(2)：492-501
[65]许卓,韩立国,廉玉广.基于等效偏移距的偏移成像方法[J].吉林大学学报,2007,37(增刊)：57-60
[66]李振春,叶月明,仝兆岐.起伏地表下基于高阶广义屏算子的叠前深度偏移[J].勘探地球物理进展,2007,30(5)：377-381
[67]胡振远,杨文采,程振炎.苏北大陆科学钻探靶区深反射地震的叠前深度偏移[J].地球物理学报,2000,43(4)：488-497
[68]Audebert, F, D. Nichols, T. Rekdal, B. Biondi, H. Iamaley, and H. Unianeta,1997, Imaging complex geologic structure with single-arrival Kirchhoff prestack depth migration. Geophysics,62,1533-1543
[69]Berkhout, A. J.,2001, Seismic imaging beyond depth migration:Geophysics,66,1895-1912
[70]Berryhill, J. R.,1979, Wave-equation datuming:Geophysics,44,1329-1344
[71]Bourgeois, A., M. Bourget, P. Lailly, M. Poulet, P. Ricarte, and R. Versteeg,1991, Marmousi, model and data:in The Marmousi experience, EAGE,5-9
[72]Chen Xiumei,Wang Huazhong, Cheng Jiubing and Ma Zaitian. Prestack depth migration by surface controlled illumination in a target oriented way. Expanded Abstracts of the 72nd S EG Mtg,2002
[73]Claerbout, J. F.,1971, Toward a unified theory of reflector imaging:Geophysics,36,3, 467-481
[74]Dave Hale,1984, Dip moveout by Fourier transform, Geophysics Vol.49, No.:P 741-757
[75]Dave Hale,1991,A nonaliases integral method for dip moveout, Geophysics Vol.56, No.6,P 795-805
[76]David Forel and Gerald H. F Gardner,1988, Three-dimension perspective on two dimension dip moveout, Geophysics Vol.53, No.5:P 604-610
[77]Stoffa, P. L., Fokkema, J. T., de Luna Freire, R. M., and Kessinger, W. P.,1990, Split-step Fourier migration:Geophysics,55,410-421.12 CREWES Research Report—Volume 18 (2006)
[78]Etgen, J. T.,1994, Stability analysis of explicit depth extrapolation through laterally varying media:64th Annual International Meeting, SEG, Expanded Abstracts, 1266-1269.183
[79]Ferguson, R.,1999, Seismic Imaging in Heterogeneous Anisotropic Media by Nonstationary Phase Shift, PhD thesis, University of Calgary, Canada
[80]Ferguson, R., G. Margrave,2002, Prestack depth migration by symmetric nonstationary phase shift:Geophysics,67,594-603
[81]Fishman, L., J. McCoy,1985, A new class of propagation models based on a factorization of the Helmholtz equation, Geophys. J. R. Astr. Soc.,80,439-461
[82]Fishman, L.,2005, Phase Space and Path Integral Methods in Seismic Wave Propagation and Imaging, Distinguished Chair Lecture Notes, PIMS, University of Calgary
[83]Francois Audebert.Integrated approach to subsalt depth imaging:Synthetic case study,TLE,2002
[84]Francois Audebert.Prestack depth migration of an Alberta Foothills data set—The Husky experience,Geophysics,1998,63 (2),pp.392-398
[85]Gazdag, J., and Sguazzero, P,1984, Migration of seismic data by phase shift plus interpolation, Geophysics,49,124-231
[86]Gazdag, J.,1978, Wave equation migration with the phase shift method:Geophysics:43, 1342-1351.
[87]Geiger, H.,2001, Relative-Amplitude-Preserving Prestack Time Migration by theEquivalent Offset Method, PhD thesis, University of Calgary.
[88]Gray, S. H., and W. May,1994, Kirchhoff migration using eikonal equation traveltimes: Geophysics,59,810-817.
[89]Gray, S. H.,1997, Where is the zero-velocity layer?:Geophysics,62,266-269.
[90]Gray, S.,2001, Seismic Imaging:Geophysics,66,15-17.
[91]Gray, S., H., J. Etgen, J. Dellinger, and D. Whitmore,2001, Seismic migration problems and solutions, Geophysics,66,1622-1640.
[92]Grossman, J. P., Margrave, G. F., and Lamoureux,M. P.,2002, Fast wavefield extrapolation by phase-shift in the nonuniform Gabor domain:CREWES Research Report, 14.
[93]Hale, D.,1991, Stable explicit depth extrapolation of seismic wavefields:Geophysics,56, 1770-1777.
[94]Holberg, O.,1988, Towards optimum one-way wave propagation, Geophys. Prosp.,36, 99-114.
[95]Huang, L. Y., M. C. Fehler, and R. S. Wu,1999, Extended local Born Fourier migration method:Geophysics,64,1524-1534.
[96]Hubral, P, Schneider, J., and Tgyel, M.,1996, A unified approach to 3—seismic reflection imaging. Part 1:Basic concepts. Geophysics,49,741-757.
[97]Kosloff, D., J. Sherwood, Z. Koren, E. Machet, and Y. Falkovitz,1996, Velocity and interface determination by tomography of depth migration gathers:Geophysics,61, 1511-1523.
[98]Lafond, C., and A. Levander,1993, Migration moveout analysis and depth focusing: Geophysics,58,91-100.
[99]Le Rousseau, J. H., and M. V. de Hoop,2001, Modeling and imaging with the scalar generalized-screen algorithms in isotropic media:Geophysics,66,1551-1568.
[100]Lee, W. B., and L. Zhang,1992, Residual shot profile migration:Geophysics,57, 812-822.
[101]MacKay, S., and R. Abma,1993, Depth-focusing analysis using a wavefront-curvature criterion:Geophysics,58,1148-1156.185
[102]MacKay, S.,1994, Efficient wavefield extrapolation to irregular surfaces using finite differences:zero-velocity datuming:64th Annual International Meeting, SEQExpanded Abstracts,1564-1567.
[103]Margrave, G. F.,1998, Direct Fourier migration for vertical velocity variations:CREWES Research Report,10.
[104]Margrave, G. F., and R. J. Ferguson,1999, Wavefield extrapolation by nonstationary phase shift:Geophysics,64,1067-1078.
[105]Margrave, G. F., and Z. Yao,2000, Downward continuation from topography with a laterally variable depth step,70th Annual International Meeting, SEG, Expanded Abstracts, 481-484.
[106]Margrave, G, F. Lamoureux, P. Gibson, R. Bale, and J. Grossman,2002, Exact wavefield extrapolation in 2D for v(x):CREWES Research Report,14.
[107]Margrave, G. F.,2003, Methods of Seismic Data Processing, Geophysics 657 Course lecture notes, University of Calgary.
[108]Margrave, G. F., H. D. Geiger, S. M. Al-Saleh, and M. P. Lamoureux,2006, Improving explicit depth migration with a stabilizing Wiener filter and spatial resampling: Geophysics,71,111-120.
[109]McClellan. J. H., and T. W. Parks,1972, Equiripple approximation of fan filters: Geophysics,37. S73-583.
[110]Margrave, G. F., Al-Saleh, S. M., Geiger, H. D., and Lamoureux, M. P.,2004, The FOCI algorithm for seismic depth migration:CREWES Research Report,16.
[111]Margrave, G. F., Geiger, H. D., Al-Saleh, S. M., and Lamoureux, M. P.,2006, Improving explicit seismic depth migration with a stabilizing wiener filter and spatial resampling: Geophysics,71, No.3, S111-S120.
[112]Nautiyal, A., Gray, S. H., and Whitmore, N. D.,1993, Stability versus accuracy for an explicit wave-field eatrapolation operatoer. Geophysics,58,277-283.
[113]O'Brien, M., and Gray, S.,1996, Can we image beneath salt? The Leading Edge,15, 17-24.
[114]Popovici, A. M., Prestack migration by split-step DSR, Geopthysics,1996,61,1412-1416.
[115]Ren, J., C. Gerrard, J. McClean, M. Orlovich, and A. Long,2005, An efficient 3D explicit finite-difference prestack depth migration:CSEG Recorder,44-46.
[116]Ristow, D., and T. Ruhl,1994, Fourier finite-difference migration:Geophysics,59,12, 1882-1893.
[117]Reshef, M.,1991,Depth migration from irregular surfaces with depth extrapolation methods,Geophysics,56th,119-122.
[118]Rietveld W E A and BerkhoutA J. Prestack depth migration by means of controlled illumination. Geophysics,1994,59 (5):801-809
[119]Ruhl, T.,1996, Finite-difference migration derived from the Kirchhoff-Helmholtz integral, Geophysics,61,1394-1399.
[120]Samuel H. Gray,What are you interpreting? Is your seismic image flawed?
[121]Samuel H. Gray, Combining geology and prestack seismic for structural model building in the Gas belt,CSEG,2006
[122]Samuel H. Gray.Seismic imaging and velocity analysis for an Alberta Foothills seismic survey,Geophysics,2001,66 (3),pp.721-732
[123]Samuel H. Gray.Gaussian beam migration of common-shot records,2005,geophysics,70 (4),pp. S71-S77,
[124]Samuel H. Gray, Angle gathers for Gaussian beam depth migration, CSEG,2006
[125]Schneider W A. Integral formulation for migration in two and three dimensions. Geophysics,1978,43:49-76
[126]Selesnick, I.W., M. Lang, and C. S. Burrus,1996, Constrained Least Squares Design of FIR Filters with Specified Transition Bands, IEEE,44,1879-1892.
[127]Shtivelman, V and Canning, A.1988, Datum correction by wave-equation extrapolation,Geophysics,53(10),1311-1322.
[128]Soubaras, R.,1996, Explicit 3-D migration using equiripple polynomial expansion and Laplacian synthesis, Geophysics,61,1386-1393.
[129]Stole, R. H.,1978, Migration by Fourier transform; Geophysics,43,23-08.
[130]Stolk,C.,and W.Symes,2004, Kinematic artifacts in prestack depth migration:Geophysics, 69,562-575.
[131]Stork, C.,1992, Reflection tomography in the post migrated domain:Geophysics, 57,680-692.
[132]Tarantola, A.,1984, Inversion of seismic reflection data in the acoustic approximation: Geophysics,49,1259-1266.
[133]Thorbecke, J., K.Wapenaar, and G Swinnen,2004, Design of one-way wavefield extrapolation operators, using smooth functions in WLSQ optimization, Geophysics,69, 1037-1045.
[134]Valenciano, A., and B. Biondi,2003,2-D deconvolution imaging condition for shotprofile migration,73rd Annual International Meeting, SEG, Expanded Abstracts,1059-1062.
[135]Varela, C.L.,P.L. Stoffa,and M.K.Sen,1998, Background velocity estimation using non-linear optimization for reflection tomography and migration misfit:Geophys. Prosp., 46,51-78.
[136]Versteeg, R., and Grau, G,1991,The Marmousi experience. Proc.1990 EAEG workshop on pratical aspects of seismic data inversion, Eur. Assoc. Expl. Geopgys.
[137]Wang Chengxiang, Zhao Bo and Zhang Guanquan.Pre2stack depth migration experiment with multi-controlled illumination. Expanded Abstracts of the 71st SEG Mtg,2001
[138]Wang, B., and K. Pann,1998, Model-based Interpretation of Focusing Panels for Depth Focusing Analysis,68th Annual International Meeting, SEG, Expanded Abstracts, 1596-1599.
[139]Wang B., F. Qin, V. Dirks, P. Guillaume, F. Audebert, and D. Epili,2005,3D finite angle tomography based on focusing analysis,75th Annual International Meeting, SEG, Expanded Abstracts,2546-2549.
[140]Withmore, N. D.,1983, Iterative depth igration, by backwards time Propagating.53rd Ann. Internat. Mtg., Soc. Expl. Geophys.
[141]Yang Kai, Wang Huaxhong and Ma Zaitian.1999, Wave equation damming from irregulatar surfaces using finite difference scheme. Expanded Abstracts of 69th SEG Mtg, 1456-1468.
[142]Zhang Guanquan, ZhangWensheng and Hao Xianjun. Prestack depth migration with common shot and synthesis shot records. Expanded Abstracts of the 69th SEG Mtg,1999
[2]冯伟,王华忠,吴如山等,面向目标的小束源照明和成像,地球物理学进展,2006,21(3)：802-808
[3]董良国,吴晓丰,唐海忠,等.逆掩推覆构造的地震波照明与观测系统优化,石油物探,2006,45(1)：40-47
[4]Liu Shouwei, Geng Jianhua, Feng Wei. Controlled illumination and seismic acquisition geometry for target-oriented imaging, Applied Geophysics,2005, V1.2(4):230-234
[5]Dong Liangguo, Cai Xiwei, et al. Seismic wave illumination for Yumen overthrust structure and seismic survey optimal design, SEG Annual Meeting,2005,48-51
[6]冯伟,吴如山,马在田.单向传播子的局部平面波分解以及在定向照明分析和目标特征成像中的应用,地球物理学进展,2007,22(1)：171-178
[7]裴正林,波动方程地震定向照明分析,石油地球物理勘探,2008,43(6)：645-651
[8]Luo M, Cao J, Xie X B, et al,2004, Comparison of illumination analysis using one-way and full-wave propagators:74th Annual International Meeting, SEG, Expanded Abstracts, 67-70
[9]吕公河,尹成,周星合,等.基于采集目标的地震照明度的精确模拟,石油地球物理勘探,2006,41(3)：258-261
[10]熊晓军,贺振华,文晓涛,等.基于波动方程的观测系统设计方法研究,中国矿业大学学报,2007,36(1),132-136
[11]赵殿栋,郭建,王咸彬,等.基于模型面向目标的观测系统优化设计技术[J].2006,2(2)：119-122
[12]陈秀梅,王华忠,程玖兵等.地表旋转控制照明叠前深度偏移[J].地球物理学报,2004,47(2)：306-311
[13]陈秀梅,王华忠,马在田.多地质目标地表控制照明叠前深度偏移[J].石油地球物理勘探,2003,38(5)：507-511
[14]吕林,李斐,李振春等.聚焦点控制照明波动方程叠前深度偏移[J].岩性油气藏,2007,19(1)：113-116
[15]林伯香,孙晶梅,徐颖,等.几种常用静校正方法的讨论[J].石油物探,2006,23(4)：368-372.
[16]杨红霞,高林.复杂地表地区地震勘探静校正方法应用研究[J].世界地质,2002,21(3)：294-298.
[17]方伍宝,李满树,孙爱萍,等.基于Born近似的波动方程静校正技术[J].石油物探,2004,43(1)：27-29.
[18]何光明,贺振华,黄德济,等.几种静校正方法的比较研究[J].物探化探计算技术,2006,28(4)：311-314.
[19]王西文,徐凤银,冯云发,等.三湖地区地震资料处理中的静校正方法研究[J].石油地球物理勘探,2007,42(5)：572-580.
[20]张海新,丁伟,刘秀玉,等.永新地区近地表校正方法[J].石油地球勘探,2008,43(增刊2)：63-66.
[21]薛为平,宋阳,刘居文,等.准噶尔盆地石南地区静校正方法应用研究[J].石油物探,2006,45(6)：615-618.
[22]李录明,罗省贤.三维波场延拓复杂表层模型校正方法及应用[J].地球物理学进展,2005,20(4)：1027-1034.
[23]朱生旺.近地表速度菲涅耳带初至层析反演[J].新疆石油地质,2007,28(6)：753-755.
[24]陈立新,王振华.SMODE1层析反演近地表模型静校正[J].天然气工业,2007,27(增刊A)：128-129.
[25]何光明,贺振华黄德济王翠华.非线性层析静校正在川西地区资料处理中的应用[J].石油物探,2006,45(1)：89-92.
[26]徐凌,崔兴福,齐莉,等.复杂近地表层析反演与波场延拓联合基准面校正[J].石油物探,2007,46(3)：227-230.
[27]韩晓丽,杨长春,麻三怀,等.复杂山区初至波层析反演静校正[J].地球物理学进展,2008,23(2)：475-483.
[28]杨兰锁,王世煜,刘财,等.射线层析静校正[J].吉林大学学报(地球科学版),2006,36(专辑)：50-75.
[29]苏贵仕,黄莉莉,赵明秋.伪三维层析反演静校正技术及其应用[J].石油地球物理堪探,2007,42(增刊)：111-114.
[30]时磊,顾汉明,刘洪雷.无射线追踪法层析静校正及其应用[J].工程地球物理学报,2008,5(5)：575-578.
[31]秦亚玲,王彦春,王瑞祥,等.用多项式拟合求取沙漠地区表层速度的静校正方法[J].天然气工业,2009,29(6)：36-39.
[32]赵传雪,王丽,吴靖,等.有限差分波动方程基准面校正方法及其在丘陵地区的应用[J].石油物探,2009,48(5)：505-509.
[33]张印堂.折射层析反演静校正及其在伊朗卡山地区的应用[J].油气地球物理,2005,3(3)：11-13.
[34]林伯香,孙晶梅,刘清林.层析成像低速带速度反演和静校正方法[J].石油物探,2002,41(2)：137-140.
[35]林伯香,肖万富,李博.层析静校正在黄土塬弯宽线资料处理中的应用[J].石油物探,2007,46(4)：418-420.
[36]林伯香,孙晶梅,徐颖,等.复杂地表条件静校正中的3D表层速度层析反演研究[J].石油 物探,2005,22(5)：455-457.
[37]陈世军,张建中.初至波射线层析成像在复杂区静校正中的应用[J].石油物探,2006,45(1)：35-39.
[38]何光明,贺振华,黄德济,等.无射线追踪层析静校技术在川西复杂地区资料处理中的应用[J].成都理工大学学报(自然科学版),2006,33(6)：599-602.
[39]左海,宋小平,魏庚雨.折射和层析静校正量计算方法的综合应用[J].石油地球物理勘探,2007(增刊)：35-39.
[40]宋涛,工秀云,杨世文,等.层析反演静校正技术在永新地区的应用[J].石油天然气学报(江汉石油学院学报),2008,30(5)：228-230.
[41]张继国,刘连升.复杂区初至层析反演静校正[J].石油地球物理勘探,2006,41(4)：383-395.
[42]Cui Xingfu, Li Hongbing, Hu Ying, Liang Hong, and Qi Li. Wavefield continuation datuming using a near surface model[J]. APPLIED GEOPHYSICS,2007,4(2):94-100
[43]冯泽元,李培明,唐海忠,等.利用层析反演技术解决山地复杂区静校正问题[J].石油物探,2005,44(3)：284-287.
[44]马在田.高阶方程偏移的分裂算法[J].地球物理学报,1983,26：377～388.
[45]马在田.论反射地震偏移成像[J].勘探地球物理进展,2002,25(3)：1-5.
[46]陈生昌,马在田.地震成像中的Born近似与Rytov近似比较,天然气工业,2007,27(1),37-40
[47]陈生昌,马在田,Wu Rushan,波动方程深度偏移的局部裂步Fourier传播算子,计算物理,2006,23(5)：604-608
[48]Rietveld W E A, Berkhout A J, Pre-stack depth migration by means of controlled illumination, Geophysics,1994,59(5):801-809
[49]Wu R S, Wang Y, and Gao J H, Beamlet migration based on local perturbation theory,70th Annual meeting, SEG Expanded abstracts.1008-1011
[50]Ye Yueming, Li Zhenchun, et al., Beamlet prestack depth migration and illumination:A test based on the Marmousi model, Applied Geophysics,2006,3(4):203-209
[51]王勇,程金星,庞全康等.弹性波正演技术在地震采集中的应用[J].小型油气藏,2007,12：21-25
[52]程玖兵,王华忠,马在田.频率空间域有限差分法叠前深度偏移[J].地球物理学报,2001,44(3)：389-395
[53]De Hoop, M, Le R—J, and Wu R. S., Generalization of the phase-screen approximation for the scattering of acoustic waves. Wave motion,2000,31:43-70
[54]李录明,罗省贤.相移加校正变范围三维叠前深度偏移[J].成都理工学院学报,1999,26(4)：329-334
[55]张叔伦,孙沛勇.基于平面波合成的傅里叶有限差分叠前深度偏移[J].石油地球物理勘探,1999,34(1)：1-7
[56]李锋,顾汉明.泌阳凹陷南部陡坡带地震资料的炮域波动方程叠前深度偏移[J].地质 科技情报,2007,26(3)：104-108
[57]曹文俊,李振春,韩文功.‘波动方程叠前深度偏移对速度模型的敏感性分析[J].天然气工业,2007,27(增刊A)：120-122
[58]马淑芳,李振春.波动方程叠前深度偏移方法综述[J].勘探地球物理进展,2007,30(3)：153-161
[59]李振春,姚云霞,马在田,等.波动方程共成像点道集偏移速度建模[J].地球物理学报,2003,46(1)：86-93
[60]徐琴,李振春,鲍伟.波动方程偏移中的成像条件研究进展[J].勘探地球物理进展,2008,31(4)：247-252
[61]杨敬磊,李振春,叶月明,等.地震照明叠前深度偏移方法综述[J].地球物理学进展,2008,23(1)：146-152
[62]刘定进,印兴耀.傅里叶有限差分法保幅叠前深度偏移方法[J].地球物理学报,2007,50(1)：268-27
[63]曹文俊,李振春,陈业全.高精度三维叠前深度偏移在DP地区的应用[J].天然气工业,2007,27(增刊A)：123-124
[64]刘定进,印兴耀.共炮检距道集波动方程保幅叠前深度偏移方法[J].地球物理学进展,2007,22(2)：492-501
[65]许卓,韩立国,廉玉广.基于等效偏移距的偏移成像方法[J].吉林大学学报,2007,37(增刊)：57-60
[66]李振春,叶月明,仝兆岐.起伏地表下基于高阶广义屏算子的叠前深度偏移[J].勘探地球物理进展,2007,30(5)：377-381
[67]胡振远,杨文采,程振炎.苏北大陆科学钻探靶区深反射地震的叠前深度偏移[J].地球物理学报,2000,43(4)：488-497
[68]Audebert, F, D. Nichols, T. Rekdal, B. Biondi, H. Iamaley, and H. Unianeta,1997, Imaging complex geologic structure with single-arrival Kirchhoff prestack depth migration. Geophysics,62,1533-1543
[69]Berkhout, A. J.,2001, Seismic imaging beyond depth migration:Geophysics,66,1895-1912
[70]Berryhill, J. R.,1979, Wave-equation datuming:Geophysics,44,1329-1344
[71]Bourgeois, A., M. Bourget, P. Lailly, M. Poulet, P. Ricarte, and R. Versteeg,1991, Marmousi, model and data:in The Marmousi experience, EAGE,5-9
[72]Chen Xiumei,Wang Huazhong, Cheng Jiubing and Ma Zaitian. Prestack depth migration by surface controlled illumination in a target oriented way. Expanded Abstracts of the 72nd S EG Mtg,2002
[73]Claerbout, J. F.,1971, Toward a unified theory of reflector imaging:Geophysics,36,3, 467-481
[74]Dave Hale,1984, Dip moveout by Fourier transform, Geophysics Vol.49, No.:P 741-757
[75]Dave Hale,1991,A nonaliases integral method for dip moveout, Geophysics Vol.56, No.6,P 795-805
[76]David Forel and Gerald H. F Gardner,1988, Three-dimension perspective on two dimension dip moveout, Geophysics Vol.53, No.5:P 604-610
[77]Stoffa, P. L., Fokkema, J. T., de Luna Freire, R. M., and Kessinger, W. P.,1990, Split-step Fourier migration:Geophysics,55,410-421.12 CREWES Research Report—Volume 18 (2006)
[78]Etgen, J. T.,1994, Stability analysis of explicit depth extrapolation through laterally varying media:64th Annual International Meeting, SEG, Expanded Abstracts, 1266-1269.183
[79]Ferguson, R.,1999, Seismic Imaging in Heterogeneous Anisotropic Media by Nonstationary Phase Shift, PhD thesis, University of Calgary, Canada
[80]Ferguson, R., G. Margrave,2002, Prestack depth migration by symmetric nonstationary phase shift:Geophysics,67,594-603
[81]Fishman, L., J. McCoy,1985, A new class of propagation models based on a factorization of the Helmholtz equation, Geophys. J. R. Astr. Soc.,80,439-461
[82]Fishman, L.,2005, Phase Space and Path Integral Methods in Seismic Wave Propagation and Imaging, Distinguished Chair Lecture Notes, PIMS, University of Calgary
[83]Francois Audebert.Integrated approach to subsalt depth imaging:Synthetic case study,TLE,2002
[84]Francois Audebert.Prestack depth migration of an Alberta Foothills data set—The Husky experience,Geophysics,1998,63 (2),pp.392-398
[85]Gazdag, J., and Sguazzero, P,1984, Migration of seismic data by phase shift plus interpolation, Geophysics,49,124-231
[86]Gazdag, J.,1978, Wave equation migration with the phase shift method:Geophysics:43, 1342-1351.
[87]Geiger, H.,2001, Relative-Amplitude-Preserving Prestack Time Migration by theEquivalent Offset Method, PhD thesis, University of Calgary.
[88]Gray, S. H., and W. May,1994, Kirchhoff migration using eikonal equation traveltimes: Geophysics,59,810-817.
[89]Gray, S. H.,1997, Where is the zero-velocity layer?:Geophysics,62,266-269.
[90]Gray, S.,2001, Seismic Imaging:Geophysics,66,15-17.
[91]Gray, S., H., J. Etgen, J. Dellinger, and D. Whitmore,2001, Seismic migration problems and solutions, Geophysics,66,1622-1640.
[92]Grossman, J. P., Margrave, G. F., and Lamoureux,M. P.,2002, Fast wavefield extrapolation by phase-shift in the nonuniform Gabor domain:CREWES Research Report, 14.
[93]Hale, D.,1991, Stable explicit depth extrapolation of seismic wavefields:Geophysics,56, 1770-1777.
[94]Holberg, O.,1988, Towards optimum one-way wave propagation, Geophys. Prosp.,36, 99-114.
[95]Huang, L. Y., M. C. Fehler, and R. S. Wu,1999, Extended local Born Fourier migration method:Geophysics,64,1524-1534.
[96]Hubral, P, Schneider, J., and Tgyel, M.,1996, A unified approach to 3—seismic reflection imaging. Part 1:Basic concepts. Geophysics,49,741-757.
[97]Kosloff, D., J. Sherwood, Z. Koren, E. Machet, and Y. Falkovitz,1996, Velocity and interface determination by tomography of depth migration gathers:Geophysics,61, 1511-1523.
[98]Lafond, C., and A. Levander,1993, Migration moveout analysis and depth focusing: Geophysics,58,91-100.
[99]Le Rousseau, J. H., and M. V. de Hoop,2001, Modeling and imaging with the scalar generalized-screen algorithms in isotropic media:Geophysics,66,1551-1568.
[100]Lee, W. B., and L. Zhang,1992, Residual shot profile migration:Geophysics,57, 812-822.
[101]MacKay, S., and R. Abma,1993, Depth-focusing analysis using a wavefront-curvature criterion:Geophysics,58,1148-1156.185
[102]MacKay, S.,1994, Efficient wavefield extrapolation to irregular surfaces using finite differences:zero-velocity datuming:64th Annual International Meeting, SEQExpanded Abstracts,1564-1567.
[103]Margrave, G. F.,1998, Direct Fourier migration for vertical velocity variations:CREWES Research Report,10.
[104]Margrave, G. F., and R. J. Ferguson,1999, Wavefield extrapolation by nonstationary phase shift:Geophysics,64,1067-1078.
[105]Margrave, G. F., and Z. Yao,2000, Downward continuation from topography with a laterally variable depth step,70th Annual International Meeting, SEG, Expanded Abstracts, 481-484.
[106]Margrave, G, F. Lamoureux, P. Gibson, R. Bale, and J. Grossman,2002, Exact wavefield extrapolation in 2D for v(x):CREWES Research Report,14.
[107]Margrave, G. F.,2003, Methods of Seismic Data Processing, Geophysics 657 Course lecture notes, University of Calgary.
[108]Margrave, G. F., H. D. Geiger, S. M. Al-Saleh, and M. P. Lamoureux,2006, Improving explicit depth migration with a stabilizing Wiener filter and spatial resampling: Geophysics,71,111-120.
[109]McClellan. J. H., and T. W. Parks,1972, Equiripple approximation of fan filters: Geophysics,37. S73-583.
[110]Margrave, G. F., Al-Saleh, S. M., Geiger, H. D., and Lamoureux, M. P.,2004, The FOCI algorithm for seismic depth migration:CREWES Research Report,16.
[111]Margrave, G. F., Geiger, H. D., Al-Saleh, S. M., and Lamoureux, M. P.,2006, Improving explicit seismic depth migration with a stabilizing wiener filter and spatial resampling: Geophysics,71, No.3, S111-S120.
[112]Nautiyal, A., Gray, S. H., and Whitmore, N. D.,1993, Stability versus accuracy for an explicit wave-field eatrapolation operatoer. Geophysics,58,277-283.
[113]O'Brien, M., and Gray, S.,1996, Can we image beneath salt? The Leading Edge,15, 17-24.
[114]Popovici, A. M., Prestack migration by split-step DSR, Geopthysics,1996,61,1412-1416.
[115]Ren, J., C. Gerrard, J. McClean, M. Orlovich, and A. Long,2005, An efficient 3D explicit finite-difference prestack depth migration:CSEG Recorder,44-46.
[116]Ristow, D., and T. Ruhl,1994, Fourier finite-difference migration:Geophysics,59,12, 1882-1893.
[117]Reshef, M.,1991,Depth migration from irregular surfaces with depth extrapolation methods,Geophysics,56th,119-122.
[118]Rietveld W E A and BerkhoutA J. Prestack depth migration by means of controlled illumination. Geophysics,1994,59 (5):801-809
[119]Ruhl, T.,1996, Finite-difference migration derived from the Kirchhoff-Helmholtz integral, Geophysics,61,1394-1399.
[120]Samuel H. Gray,What are you interpreting? Is your seismic image flawed?
[121]Samuel H. Gray, Combining geology and prestack seismic for structural model building in the Gas belt,CSEG,2006
[122]Samuel H. Gray.Seismic imaging and velocity analysis for an Alberta Foothills seismic survey,Geophysics,2001,66 (3),pp.721-732
[123]Samuel H. Gray.Gaussian beam migration of common-shot records,2005,geophysics,70 (4),pp. S71-S77,
[124]Samuel H. Gray, Angle gathers for Gaussian beam depth migration, CSEG,2006
[125]Schneider W A. Integral formulation for migration in two and three dimensions. Geophysics,1978,43:49-76
[126]Selesnick, I.W., M. Lang, and C. S. Burrus,1996, Constrained Least Squares Design of FIR Filters with Specified Transition Bands, IEEE,44,1879-1892.
[127]Shtivelman, V and Canning, A.1988, Datum correction by wave-equation extrapolation,Geophysics,53(10),1311-1322.
[128]Soubaras, R.,1996, Explicit 3-D migration using equiripple polynomial expansion and Laplacian synthesis, Geophysics,61,1386-1393.
[129]Stole, R. H.,1978, Migration by Fourier transform; Geophysics,43,23-08.
[130]Stolk,C.,and W.Symes,2004, Kinematic artifacts in prestack depth migration:Geophysics, 69,562-575.
[131]Stork, C.,1992, Reflection tomography in the post migrated domain:Geophysics, 57,680-692.
[132]Tarantola, A.,1984, Inversion of seismic reflection data in the acoustic approximation: Geophysics,49,1259-1266.
[133]Thorbecke, J., K.Wapenaar, and G Swinnen,2004, Design of one-way wavefield extrapolation operators, using smooth functions in WLSQ optimization, Geophysics,69, 1037-1045.
[134]Valenciano, A., and B. Biondi,2003,2-D deconvolution imaging condition for shotprofile migration,73rd Annual International Meeting, SEG, Expanded Abstracts,1059-1062.
[135]Varela, C.L.,P.L. Stoffa,and M.K.Sen,1998, Background velocity estimation using non-linear optimization for reflection tomography and migration misfit:Geophys. Prosp., 46,51-78.
[136]Versteeg, R., and Grau, G,1991,The Marmousi experience. Proc.1990 EAEG workshop on pratical aspects of seismic data inversion, Eur. Assoc. Expl. Geopgys.
[137]Wang Chengxiang, Zhao Bo and Zhang Guanquan.Pre2stack depth migration experiment with multi-controlled illumination. Expanded Abstracts of the 71st SEG Mtg,2001
[138]Wang, B., and K. Pann,1998, Model-based Interpretation of Focusing Panels for Depth Focusing Analysis,68th Annual International Meeting, SEG, Expanded Abstracts, 1596-1599.
[139]Wang B., F. Qin, V. Dirks, P. Guillaume, F. Audebert, and D. Epili,2005,3D finite angle tomography based on focusing analysis,75th Annual International Meeting, SEG, Expanded Abstracts,2546-2549.
[140]Withmore, N. D.,1983, Iterative depth igration, by backwards time Propagating.53rd Ann. Internat. Mtg., Soc. Expl. Geophys.
[141]Yang Kai, Wang Huaxhong and Ma Zaitian.1999, Wave equation damming from irregulatar surfaces using finite difference scheme. Expanded Abstracts of 69th SEG Mtg, 1456-1468.
[142]Zhang Guanquan, ZhangWensheng and Hao Xianjun. Prestack depth migration with common shot and synthesis shot records. Expanded Abstracts of the 69th SEG Mtg,1999