庐枞盆地金属矿地震波场精细模拟及属性应用研究
|
摘要
传统的金属矿勘查主要依赖重、磁、电等方法,一般适应于中浅层的矿产勘查,当探测目标埋深逐渐增大时,由于方法本身的限制,异常强度和分辨率急剧下降,难于适应深部矿产勘查的需求。地震方法具有探测深度大,分辨率高和探测结果准确可靠等特点,成为深部矿产勘探最有前景的技术。
为了验证金属矿区开展地震勘探的可行性,在庐枞地区罗河、泥河和大包庄矿区开展了金属矿区的地震试验,包括地震的数据采集和处理解释以及地质-地球物理模型的建立和地震波场正演模拟等研究。
与石油地震勘探相比,金属矿地震勘探有其自身的技术难点：比如金属矿区往往地质构造非常复杂,地层连续性差,探测目标小,岩石蚀变严重,地下介质不均匀性非常强,矿区通常背景干扰严重。采集的地震资料信噪比较低,缺乏标志层和测井资料,金属矿区的地震解释十分困难。
金属矿地震勘探一般需要在地质构造复杂的火山岩和变质岩有关的山区或丘陵地带进行数据采集。野外数据采集时应根据具体地质条件和参数试验合理布设激发炮点和选择采集参数,并保证检波点与地面的良好耦合,确保采集质量。
金属矿区地震数据处理时,经过对资料信噪比、频率、干扰波、静校正、动校正和速度分析等问题的具体分析及应用针对性处理技术措施后,资料品质得到明显提高,构造特征也较为清晰。合理的应用叠前和叠后去噪,提高了资料的信噪比；应用“无射线层析反演静校正”技术,较好地解决了研究区因高差和低降速带厚度变化等而引起的静校正问题,同相轴连续性有所提高；高精度交互速度分析及分频剩余静校正的应用,更加直观也更加准确的拾取叠加速度,较好地解决了低信噪比区域的速度拾取及短波长静校正问题；应用“无拉伸NMO”技术较好地解决了因拉伸畸变切除影响浅层成像的问题,浅层资料得到较大改善。
由于研究区有较多钻孔信息和岩石物性测量数据,有利于开展矿区模型的地震正演模拟。选取地震测线附近的钻孔,将钻孔投影到地震测线上,并根据这些钻孔信息,向两侧外推岩性分布,绘制出地质模型。在矿区进行物性参数测量,并统计出各种岩性的密度和速度转换关系,建立合理的物性参数库,将物性参数赋予地质模型的各个地质体,使之转换为复杂的地球物理模型；为研究地震方法对矿体和围岩的分辨能力,将复杂模型进行适当简化变为简化的地球物理模型。
然而,地下介质往往并不是均匀的,尤其是在构造活动剧烈,岩石蚀变严重的金属矿区,这种不均匀性更为显著,通常建立的均匀介质模型并不符合实际地下情况。随机介质就是非均匀介质中较为典型的一种,它主要考虑了小尺度条件下的介质非均匀性。自相关函数、尺度、方差(或者标准差)、粗糙系数是随机算法的四个要素,另外为了描述随机介质的方向性,加入了倾角参数。粗糙系数从小到大,局部由粗糙变得逐渐平滑；方差越大,介质速度扰动越大；横、纵尺度的对比关系体现了介质整体的分布特征,当横纵尺度比为几倍时,介质扰动呈现断续的似层状分布特征,当横向尺度远远大于纵向尺度时,介质呈现层状分布特征；不同倾角的随机介质能更准确的体现地下介质扰动的方向特性。本文通过建立矿区的多种随机介质模型,讨论了随机参数对正演模拟地震波场的影响。
地震波场正演模拟技术,因其直观、形象地反映地下介质中波场传播特征和规律,对人们理解矿体波场特征、解释实际地震资料,表征地下介质结构与岩性特征等,均具有重要的理论和实际意义。金属矿区地质构造复杂,矿体形态多变,地震资料的信噪比一般较低,波组特征不太明显,地震解释十分困难,对比研究正演模拟结果和实际资料,有助于认识和识别各种地震波场特征。
本文用波动方程交错网格有限差分法对各种理论模型进行正演模拟,并对模拟结果做纵波叠后、叠前时间偏移和转换波叠前时间偏移处理。结果发现：复杂模型中很多同相轴出现了一些间断,但整体形态能看到矿体和主要目标层位的反射信息；简化模型同相轴比较连续并且更清晰,突出了矿体和主要目标层位的反射信息,易于识别和分辨矿体和围岩的反射；随机介质模型波场复杂,产生了很多小的反射,形态上和实际资料更相似,偏移剖面仍能看到一些波阻抗差较大层位的反射,能够看到矿体反射。CMP叠加剖面能量在倾斜层位和矿体上没有正确归位,绕射现象明显,需要进行偏移处理；叠后克希霍夫时间偏移使同相轴归位,能够大体确定矿体的分布；叠前克希霍夫时间偏移使能量聚焦效果更好,同相轴更清晰；转换波叠前时间偏移分辨率明显提高,同时矿体边界反射更连续,矿体下边界及陡倾角层位的反射更清晰,在纵波能量不足的地方转换波能量相对更强；处理时应结合叠前克希霍夫时间偏移和转换波叠前时间偏移共同进行解释。
通过改变地震子波主频和随机介质模型中的各种随机参数,讨论频率和随机参数对地震波场的影响。发现：频率越高,地震分辨率越高,地震子波主频大于70Hz时能分辨一些复杂矿体内部的反射：随机方差大于5%时,火山岩层的反射多数淹没于随机介质产生的杂乱波场中,但由于矿体和围岩有较大波阻抗差,依然能见到矿体反射；横向相关长度较大时,产生断续的似层状地震波场；几组不同相关长度随机介质构成的组合模型产生的波场包含了这几种相关长度的组合特征；粗糙系数越大,随机介质越平滑,地震波场杂乱的现象有所减弱。另外通过建立立陡状矿体随机介质模型,并进行正演模拟,认为当矿体的随机方差远大于围岩介质随机方差时,能够看到矿体区域内部的反射和散射特征,因此随机方差相差很大时用地震方法可以探测到陡倾矿体。
根据区域地质情况和测线附近的钻孔信息,结合正演模拟和实际采集资料的地震波场特征进行地质解释,发现了沉积红盆的“箕状”结构和内部分层结构及火山岩层的3层结构。对偏移剖面进行属性分析,发现不同性质的岩石其地震动力学特征存在明显差异。从频率上看,沉积岩低频成分丰富,有第四纪覆盖的火山岩区也含低频成分,而第四纪覆盖较少的火山岩区则缺少低频信息,火山岩区高频信息相对丰富。沉积岩反射同相轴连续性较好,能量强,火山岩层和矿体反射连续性较差,多见一些凌乱的强能量团。小波时频分析的分频剖面能更好的观测每个频段各个不同反射区域的能量特征。
总体看来,矿体规模较大,产状相对平缓时,在地震剖面上可能直接观测到矿体反射信息。另外火山构造对成矿具有重要的控制作用,如果能够对控矿构造进行有效探测,则可以间接地为找矿提供有用信息。由于金属矿区地下反射的复杂性,直接探测矿体还存在一定问题,核心是如何区分矿体反射与其它反射,因此,如果作为一种勘探方法在未知区应用时,更可信的是追踪和探测已知控矿地质体,或控矿构造的深部延伸。
Traditional metal mineral exploration depends on gravity, magnetic, electric and other methods, generally adapting in the shallow. When the probe target depth increases, abnormal strength and resolution sharp decline because of the limitation of the method itself. So it is difficult to adapt to the needs of deep mineral exploration. Seismic exploration method has great depth, high resolution, and is accurate and reliable, is the most promising technology in deep mineral exploration.
To verify the feasibility of seismic exploration in metal mineral region, we carried out seismic experiment on Luohe, Nihe and Dabaozhuang deposits in Luzong area, including seismic data acquisition, processing, interpretation, geological geophysical model construction, seismic wave field modeling and other research.
Compared with the oil seismic exploration, the seismic exploration in metal mining area has its own technical difficulties:Such as geological structure in metal mining is often very complex, stratigraphic continuity is poor, the detection targets is small, rock alteration is serious, subsurface heterogeneity is very strong, and mine has usually a serious background interference. Seismic interpretation of metal mining is very difficult for the low signal to noise ratio of seismic data, the lack of marker beds and log data.
Seismic data acquisition generally occurs in mountainous or hilly areas where volcanic metamorphic rocks have complex geological structure. Field data acquisition should select reasonable shot position and acquisition parameters based on circumstances and parameters test, and ensure receivers coupled well with the ground, ensure acquisition quality.
Metal mining seismic data processing, through the specific analysis of SNR, frequency, interference wave, static correction. NMO correction, velocity et al., and the application of specific measures, while processing seismic data of metal mining, data quality has been significantly improved, structural features is also more clear. Reasonable application of prestack and post-stack noise removal, the SNR of seismic data is improved; application of "no-ray tomography inversion static correction" technology, static correction problem due to elevation changes and low-velocity layer thickness is better resolved in the study area, and the continuity of phase axis is improved; The application of high-precision interactive velocity analysis and residual static correction of frequency splitting, more intuitive and more accurate stacking velocity picking, can solve the speed picked up and short wavelength static problem in the low SNR region; Application of "no stretch NMO" technology can solve the shallow imaging problems affected due to removal of the stretch distortion, shallow data improved greatly.
Since the study area have more information of boreholes and petrophysical measurement data, mining model is conducive to carry out seismic modeling. Select boreholes nearby seismic lines, and project onto the seismic, and extrapolate to both sides of the rock distribution based on these drilling information, then draw a geological model. Take the petrophysical parameters measurement, and establish reasonable petrophysical parameters library using the statistics results of the relationship of density and velocity, and then convert the geological model to the complex geophysical model by put the petrophysical parameters to the geological bodies; To research the resolution of Seismic methods to the ore body and surrounding rock, complex model will be transformed to a simplified geophysical model.
However, the underground media are often not uniform, especially in the metal mining with intense tectonic activity and serious alteration of rock, this heterogeneity is more significant. The homogeneous medium mode established is often not meet with the actual ground conditions. Random media is a typical heterogeneous medium, it is mainly considered under the conditions of small-scale. Autocorrelation function, scale, variance (or standard deviation), roughness coefficient are the four elements for the random algorithm. Roughness coefficient from small to large, local media change from rough to smoothing; Variance more larger, the disturbance of velocity is more larger; The relation of horizontal and vertical scale reflects the overall distribution of the media, when the ratio of horizontal and vertical scale is several times, the media shows like layered distribution with intermittent disturbances, when the horizontal scale is much larger than the vertical scale, the media presentation layer distribution; Different angle of the random medium can more accurately reflects the direction of subsurface media disturbance characteristics. Discuss the influence of random parameters to the forward modeling of seismic wave field by establishing a variety of mining random medium models.
Seismic wave field modeling technique have important theoretical and practical significance, for its visual image to reflect the characteristics and the law of subsurface wave field propagation, people understand the wave field characteristics of the ore bodies, to explain the seismic data, characterization of subsurface structure and lithology. Geological structures are complex in metal mining, ore body shape is changeable, SNR is generally low, characteristics of composite waves is less obvious, seismic interpretation is very difficult. Compared to the results of forward modeling and practical information to help understand and identify characteristics of various seismic wave field.
The various theoretical models are modeled using wave equation with staggered-grid finite difference method, and then do wave post-stack, the pre-stack time migration of p-wave and pre-stack time migration of converted wave. The results showed that:phase axis is intermittent partly in complex model, but the reflection information of overall shape of ore bodies and the main objective layers can be seen; The phase axis is more continuous and clearly in simplified model, highlight the reflection information of ore bodies and the main objective layers; The wave field is very complex in random medium model, resulting in many small reflections, is more similar to practical information in form, Some reflection of ore bodies and the layers with larger impedance difference can be still seen in the migration profile. Energy of CMP stack section in an tilted layer and ore bodies does not have the correct homing, diffraction phenomenon is obvious, so needed for migration; Kirchhoff post-stack time migration can generally replace phase axis and determine the distribution of ore bodies; Kirchhoff pre-stack time migration is better to focus the energy, has clearer phase axis; Converted wave pre-stack time migration significantly improve the resolution, while the reflection of ore boundary is more continuous, the reflections of the lower boundary of ore bodies and the steep angle layers are more clearly, converter Wave energy is relatively stronger where in the lack of p-wave reflection energy place; So p-wave Kirchhoff pre-stack time migration and converted wave pre-stack time migration should be combined to do interpretation.
Discuss the influence of frequencies and random parameters to seismic wave field, through changing the seismic wave dominant frequencies and random parameters. Found:the higher the frequency, the seismic resolution is better, when the seismic wave frequency is higher than 70Hz, some internal reflection of the orebody can be distinguished; When random variance is greater than 5%, the majority reflection of volcanic layer submerge in the random medium wave field, but still able to see the ore body reflection for the large impedance difference between ore body and surrounding; When the horizontal correlation length is large, generate intermittent liked layered seismic wave field; when roughness coefficient is greater, the random medium on smooth, the messy phenomenon of seismic wave field has been declined. Through the establishment of steep-shaped ore bodies random medium model, and the forward modeling, the reflection and scattering characteristics within the ore bodies can be seen when the random variance of ore body is much larger than the surrounding rock. So when the random variance is large difference, steeply dipping ore bodies can be detected using seismic methods.
Carry out geological interpretation based on regional geology and drillhole information near the seismic line, combined with forward simulation and actual data characteristics of the seismic wave field, and found that the deposition red basin's "Kei-like" structure and the internal hierarchical structure and volcanic 3 layers structure. Attribute analysis of migration profiles can found that the different nature of the rock significantly different dynamic characteristics. From the frequency, the sedimentary rock is rich in low frequency, there are also low-frequency components in the volcanic region covered with quaternary, and volcanic region without quaternary coverage is lack of low frequency, high frequency information is relatively rich in volcanic areas. Sedimentary areas is better continuity of phase axis, the energy intensity, volcanic areas is poor continuity, and has some messy strong energy group. Wavelet time-frequency analysis can better observe the energy characteristics of different reflection region of the each frequency band profile.
Generally speaking, the large and relatively flat ore bodies may generate effective information in the seismic reflection profile. And volcanic structures is very important to control the mineralization. If it can be effectively detected, then will provide useful information for prospecting. There are still some problems of direct detection of ore bodies, because of the complexity of the subsurface reflection, the core of the problem is how to distinguish between ore bodies reflection and other reflections. So when the seismic exploration method is used in unknown area, the tracking and detection of known ore-controlling geologic body or structures in the deep extension is more credible.
为了验证金属矿区开展地震勘探的可行性,在庐枞地区罗河、泥河和大包庄矿区开展了金属矿区的地震试验,包括地震的数据采集和处理解释以及地质-地球物理模型的建立和地震波场正演模拟等研究。
与石油地震勘探相比,金属矿地震勘探有其自身的技术难点：比如金属矿区往往地质构造非常复杂,地层连续性差,探测目标小,岩石蚀变严重,地下介质不均匀性非常强,矿区通常背景干扰严重。采集的地震资料信噪比较低,缺乏标志层和测井资料,金属矿区的地震解释十分困难。
金属矿地震勘探一般需要在地质构造复杂的火山岩和变质岩有关的山区或丘陵地带进行数据采集。野外数据采集时应根据具体地质条件和参数试验合理布设激发炮点和选择采集参数,并保证检波点与地面的良好耦合,确保采集质量。
金属矿区地震数据处理时,经过对资料信噪比、频率、干扰波、静校正、动校正和速度分析等问题的具体分析及应用针对性处理技术措施后,资料品质得到明显提高,构造特征也较为清晰。合理的应用叠前和叠后去噪,提高了资料的信噪比；应用“无射线层析反演静校正”技术,较好地解决了研究区因高差和低降速带厚度变化等而引起的静校正问题,同相轴连续性有所提高；高精度交互速度分析及分频剩余静校正的应用,更加直观也更加准确的拾取叠加速度,较好地解决了低信噪比区域的速度拾取及短波长静校正问题；应用“无拉伸NMO”技术较好地解决了因拉伸畸变切除影响浅层成像的问题,浅层资料得到较大改善。
由于研究区有较多钻孔信息和岩石物性测量数据,有利于开展矿区模型的地震正演模拟。选取地震测线附近的钻孔,将钻孔投影到地震测线上,并根据这些钻孔信息,向两侧外推岩性分布,绘制出地质模型。在矿区进行物性参数测量,并统计出各种岩性的密度和速度转换关系,建立合理的物性参数库,将物性参数赋予地质模型的各个地质体,使之转换为复杂的地球物理模型；为研究地震方法对矿体和围岩的分辨能力,将复杂模型进行适当简化变为简化的地球物理模型。
然而,地下介质往往并不是均匀的,尤其是在构造活动剧烈,岩石蚀变严重的金属矿区,这种不均匀性更为显著,通常建立的均匀介质模型并不符合实际地下情况。随机介质就是非均匀介质中较为典型的一种,它主要考虑了小尺度条件下的介质非均匀性。自相关函数、尺度、方差(或者标准差)、粗糙系数是随机算法的四个要素,另外为了描述随机介质的方向性,加入了倾角参数。粗糙系数从小到大,局部由粗糙变得逐渐平滑；方差越大,介质速度扰动越大；横、纵尺度的对比关系体现了介质整体的分布特征,当横纵尺度比为几倍时,介质扰动呈现断续的似层状分布特征,当横向尺度远远大于纵向尺度时,介质呈现层状分布特征；不同倾角的随机介质能更准确的体现地下介质扰动的方向特性。本文通过建立矿区的多种随机介质模型,讨论了随机参数对正演模拟地震波场的影响。
地震波场正演模拟技术,因其直观、形象地反映地下介质中波场传播特征和规律,对人们理解矿体波场特征、解释实际地震资料,表征地下介质结构与岩性特征等,均具有重要的理论和实际意义。金属矿区地质构造复杂,矿体形态多变,地震资料的信噪比一般较低,波组特征不太明显,地震解释十分困难,对比研究正演模拟结果和实际资料,有助于认识和识别各种地震波场特征。
本文用波动方程交错网格有限差分法对各种理论模型进行正演模拟,并对模拟结果做纵波叠后、叠前时间偏移和转换波叠前时间偏移处理。结果发现：复杂模型中很多同相轴出现了一些间断,但整体形态能看到矿体和主要目标层位的反射信息；简化模型同相轴比较连续并且更清晰,突出了矿体和主要目标层位的反射信息,易于识别和分辨矿体和围岩的反射；随机介质模型波场复杂,产生了很多小的反射,形态上和实际资料更相似,偏移剖面仍能看到一些波阻抗差较大层位的反射,能够看到矿体反射。CMP叠加剖面能量在倾斜层位和矿体上没有正确归位,绕射现象明显,需要进行偏移处理；叠后克希霍夫时间偏移使同相轴归位,能够大体确定矿体的分布；叠前克希霍夫时间偏移使能量聚焦效果更好,同相轴更清晰；转换波叠前时间偏移分辨率明显提高,同时矿体边界反射更连续,矿体下边界及陡倾角层位的反射更清晰,在纵波能量不足的地方转换波能量相对更强；处理时应结合叠前克希霍夫时间偏移和转换波叠前时间偏移共同进行解释。
通过改变地震子波主频和随机介质模型中的各种随机参数,讨论频率和随机参数对地震波场的影响。发现：频率越高,地震分辨率越高,地震子波主频大于70Hz时能分辨一些复杂矿体内部的反射：随机方差大于5%时,火山岩层的反射多数淹没于随机介质产生的杂乱波场中,但由于矿体和围岩有较大波阻抗差,依然能见到矿体反射；横向相关长度较大时,产生断续的似层状地震波场；几组不同相关长度随机介质构成的组合模型产生的波场包含了这几种相关长度的组合特征；粗糙系数越大,随机介质越平滑,地震波场杂乱的现象有所减弱。另外通过建立立陡状矿体随机介质模型,并进行正演模拟,认为当矿体的随机方差远大于围岩介质随机方差时,能够看到矿体区域内部的反射和散射特征,因此随机方差相差很大时用地震方法可以探测到陡倾矿体。
根据区域地质情况和测线附近的钻孔信息,结合正演模拟和实际采集资料的地震波场特征进行地质解释,发现了沉积红盆的“箕状”结构和内部分层结构及火山岩层的3层结构。对偏移剖面进行属性分析,发现不同性质的岩石其地震动力学特征存在明显差异。从频率上看,沉积岩低频成分丰富,有第四纪覆盖的火山岩区也含低频成分,而第四纪覆盖较少的火山岩区则缺少低频信息,火山岩区高频信息相对丰富。沉积岩反射同相轴连续性较好,能量强,火山岩层和矿体反射连续性较差,多见一些凌乱的强能量团。小波时频分析的分频剖面能更好的观测每个频段各个不同反射区域的能量特征。
总体看来,矿体规模较大,产状相对平缓时,在地震剖面上可能直接观测到矿体反射信息。另外火山构造对成矿具有重要的控制作用,如果能够对控矿构造进行有效探测,则可以间接地为找矿提供有用信息。由于金属矿区地下反射的复杂性,直接探测矿体还存在一定问题,核心是如何区分矿体反射与其它反射,因此,如果作为一种勘探方法在未知区应用时,更可信的是追踪和探测已知控矿地质体,或控矿构造的深部延伸。
Traditional metal mineral exploration depends on gravity, magnetic, electric and other methods, generally adapting in the shallow. When the probe target depth increases, abnormal strength and resolution sharp decline because of the limitation of the method itself. So it is difficult to adapt to the needs of deep mineral exploration. Seismic exploration method has great depth, high resolution, and is accurate and reliable, is the most promising technology in deep mineral exploration.
To verify the feasibility of seismic exploration in metal mineral region, we carried out seismic experiment on Luohe, Nihe and Dabaozhuang deposits in Luzong area, including seismic data acquisition, processing, interpretation, geological geophysical model construction, seismic wave field modeling and other research.
Compared with the oil seismic exploration, the seismic exploration in metal mining area has its own technical difficulties:Such as geological structure in metal mining is often very complex, stratigraphic continuity is poor, the detection targets is small, rock alteration is serious, subsurface heterogeneity is very strong, and mine has usually a serious background interference. Seismic interpretation of metal mining is very difficult for the low signal to noise ratio of seismic data, the lack of marker beds and log data.
Seismic data acquisition generally occurs in mountainous or hilly areas where volcanic metamorphic rocks have complex geological structure. Field data acquisition should select reasonable shot position and acquisition parameters based on circumstances and parameters test, and ensure receivers coupled well with the ground, ensure acquisition quality.
Metal mining seismic data processing, through the specific analysis of SNR, frequency, interference wave, static correction. NMO correction, velocity et al., and the application of specific measures, while processing seismic data of metal mining, data quality has been significantly improved, structural features is also more clear. Reasonable application of prestack and post-stack noise removal, the SNR of seismic data is improved; application of "no-ray tomography inversion static correction" technology, static correction problem due to elevation changes and low-velocity layer thickness is better resolved in the study area, and the continuity of phase axis is improved; The application of high-precision interactive velocity analysis and residual static correction of frequency splitting, more intuitive and more accurate stacking velocity picking, can solve the speed picked up and short wavelength static problem in the low SNR region; Application of "no stretch NMO" technology can solve the shallow imaging problems affected due to removal of the stretch distortion, shallow data improved greatly.
Since the study area have more information of boreholes and petrophysical measurement data, mining model is conducive to carry out seismic modeling. Select boreholes nearby seismic lines, and project onto the seismic, and extrapolate to both sides of the rock distribution based on these drilling information, then draw a geological model. Take the petrophysical parameters measurement, and establish reasonable petrophysical parameters library using the statistics results of the relationship of density and velocity, and then convert the geological model to the complex geophysical model by put the petrophysical parameters to the geological bodies; To research the resolution of Seismic methods to the ore body and surrounding rock, complex model will be transformed to a simplified geophysical model.
However, the underground media are often not uniform, especially in the metal mining with intense tectonic activity and serious alteration of rock, this heterogeneity is more significant. The homogeneous medium mode established is often not meet with the actual ground conditions. Random media is a typical heterogeneous medium, it is mainly considered under the conditions of small-scale. Autocorrelation function, scale, variance (or standard deviation), roughness coefficient are the four elements for the random algorithm. Roughness coefficient from small to large, local media change from rough to smoothing; Variance more larger, the disturbance of velocity is more larger; The relation of horizontal and vertical scale reflects the overall distribution of the media, when the ratio of horizontal and vertical scale is several times, the media shows like layered distribution with intermittent disturbances, when the horizontal scale is much larger than the vertical scale, the media presentation layer distribution; Different angle of the random medium can more accurately reflects the direction of subsurface media disturbance characteristics. Discuss the influence of random parameters to the forward modeling of seismic wave field by establishing a variety of mining random medium models.
Seismic wave field modeling technique have important theoretical and practical significance, for its visual image to reflect the characteristics and the law of subsurface wave field propagation, people understand the wave field characteristics of the ore bodies, to explain the seismic data, characterization of subsurface structure and lithology. Geological structures are complex in metal mining, ore body shape is changeable, SNR is generally low, characteristics of composite waves is less obvious, seismic interpretation is very difficult. Compared to the results of forward modeling and practical information to help understand and identify characteristics of various seismic wave field.
The various theoretical models are modeled using wave equation with staggered-grid finite difference method, and then do wave post-stack, the pre-stack time migration of p-wave and pre-stack time migration of converted wave. The results showed that:phase axis is intermittent partly in complex model, but the reflection information of overall shape of ore bodies and the main objective layers can be seen; The phase axis is more continuous and clearly in simplified model, highlight the reflection information of ore bodies and the main objective layers; The wave field is very complex in random medium model, resulting in many small reflections, is more similar to practical information in form, Some reflection of ore bodies and the layers with larger impedance difference can be still seen in the migration profile. Energy of CMP stack section in an tilted layer and ore bodies does not have the correct homing, diffraction phenomenon is obvious, so needed for migration; Kirchhoff post-stack time migration can generally replace phase axis and determine the distribution of ore bodies; Kirchhoff pre-stack time migration is better to focus the energy, has clearer phase axis; Converted wave pre-stack time migration significantly improve the resolution, while the reflection of ore boundary is more continuous, the reflections of the lower boundary of ore bodies and the steep angle layers are more clearly, converter Wave energy is relatively stronger where in the lack of p-wave reflection energy place; So p-wave Kirchhoff pre-stack time migration and converted wave pre-stack time migration should be combined to do interpretation.
Discuss the influence of frequencies and random parameters to seismic wave field, through changing the seismic wave dominant frequencies and random parameters. Found:the higher the frequency, the seismic resolution is better, when the seismic wave frequency is higher than 70Hz, some internal reflection of the orebody can be distinguished; When random variance is greater than 5%, the majority reflection of volcanic layer submerge in the random medium wave field, but still able to see the ore body reflection for the large impedance difference between ore body and surrounding; When the horizontal correlation length is large, generate intermittent liked layered seismic wave field; when roughness coefficient is greater, the random medium on smooth, the messy phenomenon of seismic wave field has been declined. Through the establishment of steep-shaped ore bodies random medium model, and the forward modeling, the reflection and scattering characteristics within the ore bodies can be seen when the random variance of ore body is much larger than the surrounding rock. So when the random variance is large difference, steeply dipping ore bodies can be detected using seismic methods.
Carry out geological interpretation based on regional geology and drillhole information near the seismic line, combined with forward simulation and actual data characteristics of the seismic wave field, and found that the deposition red basin's "Kei-like" structure and the internal hierarchical structure and volcanic 3 layers structure. Attribute analysis of migration profiles can found that the different nature of the rock significantly different dynamic characteristics. From the frequency, the sedimentary rock is rich in low frequency, there are also low-frequency components in the volcanic region covered with quaternary, and volcanic region without quaternary coverage is lack of low frequency, high frequency information is relatively rich in volcanic areas. Sedimentary areas is better continuity of phase axis, the energy intensity, volcanic areas is poor continuity, and has some messy strong energy group. Wavelet time-frequency analysis can better observe the energy characteristics of different reflection region of the each frequency band profile.
Generally speaking, the large and relatively flat ore bodies may generate effective information in the seismic reflection profile. And volcanic structures is very important to control the mineralization. If it can be effectively detected, then will provide useful information for prospecting. There are still some problems of direct detection of ore bodies, because of the complexity of the subsurface reflection, the core of the problem is how to distinguish between ore bodies reflection and other reflections. So when the seismic exploration method is used in unknown area, the tracking and detection of known ore-controlling geologic body or structures in the deep extension is more credible.
引文
[1]徐明才,高景华等.金属矿地震勘探(M).北京：地质出版社,2009：1-263.
[2]Eaton D W, Milkereit B, Salisbury M. Seismic methods for deep mineral exploration:mature technologies adapted to new targets[J]. The Leading Edge, 2003,22(6):580-585.
[3]Aminzadeh F. Future geophysical technology trends [M]. The Leading Edge, 1996:739-742.
[4]吕庆田,廉玉广,赵金花.反射地震技术在成矿地质背景与深部矿产勘查中的应用：现状与前景[J].地质学报,2010,84(6)：771-787.
[5]Brown L, Barazangi M, Kaufman S, et al. The first decade of COCORP:1974-1984. in Reflection Profiling:A Global Perspective[M]// Barazangi M and Brown L, eds, Geodynamics Series,13:107-120, American Geophysical Union, Washington, D C.1986.
[6]Matthews D H, Smith C., Deep seismic reflection profiling of the continental lithosphere[J]. Geophys. J. Roy. Astr. Soc.,1987,89:209-216
[7]Klemperer S L, Hobbs R. The BIRPS Atlas[M]. edited by Cambridge University Press, Cambridge, UK, hardcover,1992.
[8]John Wilson. LITHOPROBE:Dancing Elephants & Floating Continents-The story of Canada beneath your feet[M]., Published by Key Porter Books,2003.
[9]Australia Geodynamics Cooperative Research Center. Annual Report 1998/1999. 1999.
[10]Cook F, Brown L, Kaufman S, et al. COCORP seismic reflection profiling of the Appalachian orogen beneath the Coastal Plain of Georgia[J]. Geol. Soc. Amer. Bull.,92:738-748.
[11]Brun J P, Wenzel F, ECORS-DEKORP team. Crustal-scale structure of the southern Rhinegraben from ECORS-DEKORP seismic reflection data[J].. Geology,1981,19:758-762
[12]Zhao W J, Nelson K D. Project INDEPTH Team. Deep seismic reflection evidence for continental underthrusting beneath southern Tibet[J]. Nature,1993, 366:557-559.
[13]Nelson K D, Zhao W J. Project INDEPTH Team. Partially molten middle crust beneath Southern Tibet:Synthesis of Project INDEPTH results[J]. Science,1996, 274:1684-1688.
[14]董树文,李廷栋SinoProbe—中国深部探测实验[J].地质学报,2009,83(9)：895-909.
[15]Perron G and Calvert A J. Shallow, high-resolution seismic imaging at the ansil mining camp in the Abitibi greenstone belt[J]. Geophysics,1998, 63(2):379-391.
[16]Milkereit B, Green A, Sudbury Working Group.1992. Deep geometry of the Sudbury structure from seismic reflection profiling[J]. Geology,1992, 20:807-811.
[17]Roy B, Clowes R M. Seismic and potential-field imaging of the Guichon Creek batholith, Brithsh Columia, Canada, to delineate structures hosting porphyry copper deposits[J]. Geophysics,2000,65(5):1418-1434
[18]Pretorius C C, Trewick W F, Fourie A, Irons C. Application of 3-D seismic to mine planning at Vaal Reefs gold mine, number 10 shaft, Republic of South Africa[J]. Geophysics,2000,65 (6):1862-1870
[19]Stuart G W, Jolley S J, Polome L G. Application of 3-D seismic attributes analysis to mine planning:Target gold deposit, South Africa[J]. The Leading Edge 2000,19(7):736-742.
[20]White D, Boerner D, Wu J J, et al. Mineral exploration in the Thompson nickel belt, Manitoba, Canada, using seismic and controlled-source EM methods[J]. Geophysics,2000,65(6):1871-1881.
[21]Milkereit B, Berrer E K, King A R, et al. Development of 3-D seismic exploration technology for deep nickel-copper deposits-A case history from the Sudbury basin, Canada[J]. Geophysics,2000,65(6):1890-1899.
[22]吕庆田,侯增谦,史大年等.铜陵狮子山金属矿地震探测结果及对区域找矿的意义[J].矿床地质,2004,23(3)：390-398.
[23]吕庆田,史大年,赵金花,严加永,徐明才.深部矿产勘查的地震学方法：问题与前景一铜陵矿集区的应用实例[J].地质通报,2005,24(3)：211-218.
[24]吕庆田,韩立国,严加永,廉玉广,史大年,颜庭杰.庐枞矿集区火山气液型铁硫矿床及控矿构造的反射地震成像[J].岩石学报,2010,26(9)：2598-2612.
[25]徐明才,高景华,荣立新,等.地面地震层析技术在金属矿勘查中的试验研究[J].物探与化探,2005,29(4)：299-303.
[26]史大年,吕庆田,徐明才.铜陵矿集区地壳浅表结构的地震层析研究[J].矿床地质,2004,23(3)：383-389.
[27]Li T L, Eaton D W. Delineating the Tuwu porphyry copper deposit at Xinjiang, China, with seismic-reflection profiling[J]. Geophysics.2005,70(6):B53-B60.
[28]梁光河,蔡新平.地震勘探在山东蓬家夼金矿深部预测中的应用[J].矿产与地质,2001,15(6)：743-748.
[29]徐明才,高景华,柴铭涛,等.金属矿地震勘查的方法技术[J].岩土工程界,1997,6(4)：41-46.
[30]徐明才,高景华,柴铭涛,等.寻找隐伏金属矿的地震方法技术研究[J].物探与化探,1997,21(6)：69-75.
[31]徐明才,高景华.从金属矿地震方法的试验效果探讨其应用前景[J].中国地质,2004,31(1)：110-114.
[32]徐明才,高景华,荣立新,等.散射波地震方法在蔡家营多金属矿区的试验研究[J].物探与化探,2003,27(1)：52-57.
[33]高景华,徐明才,荣立新,等.小热泉子铜矿区地震方法试验研究[J].地质与勘探，2004,40(6)：49-54.
[34]勾丽敏,刘学伟,雷鹏,等.金属矿地震勘探技术方法研究综述—散射波地震勘探方法[J].勘探地球物理进展,2007，30(2)：65-90.
[35]Mitchell, A H G, Garson M S. Mineral deposits and global tectonic settings[M]. New York:Academic Press,1981:1-405.
[36]Sawkins F J. Metal deposits in relation to plate tectonics[M]. Berlin:Springer Verlag,1984:1-325.
[37]Kerrich R, Goldfarb R J, Richards J P. Metallogenic Provinces in an Evolving Geodynamic Framework[J]. Economic Geology 100th anniversary Volume, 2005:1097-1136
[38]Beaumont C, Quinlan G. A geodynamic framework for interpreting crustal-scale seismic-reflectivity patterns in compressional orogens[J]. Geophysics Journal International,1994,116:754-783.
[39]Drummond B J, Goleby B R, Owen A J, Yeates A N, Swager C, Zhang Y, and Jackson J K. Seismic reflection imaging of mineral systems:Three case histories[J]. Geophysics,2000,65(6):1852-1861.
[40]Oliver N H S. Linking of regional and local hydrothermal systems in the mid-crust by shearing and faulting[J]. Tectonophysics,2001,335:147-161.
[41]Goleby B R, Blewett R S, Korsch R J, et al. Deep seismic reflection profiling in the Archaean northeastern Yilgarn Craton, Western Australia:Implication for crustal architecture and mineral potential [J]. Tectonophysics,2004,388: 119-133.
[42]Drummond B J, Goleby B R. Seismic reflection images of major ore-controlling structure in the Eastern Goldfields,Western Australia[J]. Expl. Geophys.,1993, 24:473-478.
[43]Swager C P, Goleby B R, Drummond B J, Rattenbury M S, Williams P R. Crustal structure of granite-greenstone terranes in the Eastern Goldfields, Yilgarn Craton, as revealed by seismic reflection profiling[J]. Precambrian Research,1997, 83(1-3):43-56.
[44]Goleby B R, Rattenbury M S, Swager, Swager, C P, et al. Archaean crustal structure from seismic reflection profiling, Eastern Goldfields, Western Australia[C]. Austral. Geol. Surv. Org. record,1993/15.
[45]Upton P, Hobbs B, Ord A, Zhang Y, Drummond B, Archibald N. Thermal and deformation modelling of the Yilgarn deep seismic transect[C]. Geodynamics and Ore Deposits Conference, Australia Geodynamics Cooperative Research Centre, Ballarat, Victoria,1997:22-25.
[46]Sorjonen-Ward P, Zhang Y, Zhao C. Numerical modelling of orogenic processes and gold mineralisation in the southeastern part of the Yilgarn Craton, Western Australia[J]. Australian Journal of Earth Sciences,2002,49:935-964.
[47]Jones T, Nur A. The Nature of seismic reflection from deep crustal fault zones[J]. Geophys. Res.,1984,89:3153-3171
[48]Malehmir A, Tryggvason A, Lickorish H, Weihed P. Regional structural profiles in the western part of the Palaeoproterozoic Skellefte ore district, northern Sweden[J]. Precambrin Research,2007,159:1-18.
[49]Tryggvason A, Malehmir A, Rodriguez-Tablante J, et al. Reflection seismic investigation in the western part of the palaeoproterozoic VHMS bearing Skellefte Districts, northern Sweden[J]. Economic Geology,2006, 101:1039-1054.
[50]Salisbury M H, Milkereit B, and Bleeker W. Seismic imaging of massive sulfide deposits:Part I. rock properies[J]. Economic Geology,1996,91(5):821-828.
[51]Salibury M H, Harvey C W, and Matthews. The acoustic properties of ores and host rocks in hardrock terranes[M]//In D. W. Eaton., B. Milkereit, M. H. Salisbury, Eds., Hardrock Seismic Exploration,2003.
[52]Bohlen T, Muller C, Milkereit B.2003. Elastic seismic-wave scattering from massive sulfide orebodies:on the role of composition and shape[M]//In D. W. Eaton., B. Milkereit, M. H. Salisbury, Eds., Hardrock Seismic Exploration, Geophysical developments.2003,10:70-90.
[53]Clarke G J, Eaton D W. Influence of morphology and surface roughness on the seismic response of massive sulfides, based on elastic-wave kirchhoff modeling[M]. In D. W. Eaton., B. Milkereit, M. H. Salisbury, Eds., Hardrock Seismic Exploration, Geophysical developments,2003,10:45-58.
[54]Eaton D E. Weak elastic scatting from massive sulfide ore bodies[J]. Geophysics. 1999,64(1):289-299.
[55]Hobbs R W.2003.3D modelling of seismic-wave propagation using complex elastic screen, with application to mineral exploration[M]//In D. W. Eaton., B. Milkereit, M. H. Salisbury, Eds., Hardrock Seismic Exploration, Geophysical developments,2003,10:59-69.
[56]Adame E, Milkereit B, and Salmon B.2004.3D seismic exploration in the Val d'Or mining camp, Quebec[C]. SEG Int'l Exposition and 74th Annual Meeting, Denver, Colorado,10-15 October,2004.
[57]Calvert A J, and Li Y.Seismic reflection imaging over a massive sulfide deposit at the Matagami mining camp, Quebec[J]. Geophysics,1999,64(l):24-32.
[58]Adam E, Perron G, Arnold G, et al.3D seismic imaging for VMS deposit exploration, Matagami, Quebec[M]//In D. W. Eaton., B. Milkereit, M. H. Salisbury, Eds., Hardrock Seismic Exploration, Geophysical developments,2003, 10:229-246.
[59]Milkereit B, Berrer E K, King A R, et al. Development of 3D seismic exploration technology for deep nickel-copper deposits-A case history from the Sudbury basin, Canada[J]. Geophysics,2000,65(6):1890-1899.
[60]Salisbury M H, Milkereit B, Ascough G, et al. Physical properties and seismic imaging of massive sulfides[J]. Geophysics,2000,65(6):1882-1889.
[61]Matthews L. Base metal exploration:Looking deeper and adding value with seismic data[C]. Canadian Society of Exploration Geophysicists Recorder, 2002:37-43.
[62]Bellefleur G, Muller C, Snyder D, Matthews L. Downhole seismic imaging of a massive sulfide orebody with mode-converted waves, Halfmile Lake, New Brunswick, Canada[J]. Geophysics,2004,69(2):318-329.
[63]Fountain D, Salisbury M H, Percival J. Seismic structure of the continental crust based on rock velocity measurements from the Kapuskasing uplift[J]. J. Geophys. Res.,1990,95:1167-1186.
[64]Pretorius C C, Jamison A A, Irons C. Seismic exploration in the Witwatersrand basin, Republic of South Africa[J]. Exploration 87, Proceedings,1987,3: 241-253.
[65]Roberts B. Zaleski E, Perron G et al. Seismic Exploration of the Manitouwadge Greenstone Belt, Ontario:A case history [M]//In D.W. Eaton.,B.Milkereit, M. H. Salisbury, Eds.,Hardrock Seismic Exploration, Geophysical developments,2003, 10:110-126.
[66]Aki K. Analysis of the Seismic Coda of Local Earthquakes as Scattered Waves[J]. Geophysics,1969,74(4):615-631.
[67]Godfrey R, Muir F, Rocca F. Modeling seismic impedance with Markovchains[J]. Geophysics,1980,45:1351-1372.
[68]Frankel A, Clayton R W, A finite-difference simulation of wave propagation in two-dimensional random media[M].Bulletin of the Seismological Society of America,1984,74(6):2167-2186.
[69]Frankel A, Clayton R W. Finite-difference simulation of seismic scattering: Implications for the propagation of short-period seismic waves in the crust and models of crustal heterogeneity [J]. Journal of Geophysical Research,1986,91: 6465-6489.
[70]Stanke F E, Burridge R. Comparison of spatial and ensemble averaging methods applied to wave propagation in finely layered media[C].60th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts,1990:1062-1065.
[71]De Hoop M V, Bumidge R, Chang H W. Wave propagation with tunnelingin a highly discontinuous layered medium[J]. Wave Motion,1991,13:307-329.
[72]Kerner C. Anisotropy in sedimentary rocks modeled as random media[J]. Geophysics,1992,57:564-576.
[73]Ikelle L T, Yung S K.2-D random media with ellipsoidal autocorrelation functions[J]. Geophysics,1993,58(9):1359-1372.
[74]Kneib G, Kerner C. A ccurate and efficient seismic modeling in random media[J]. Geophysics,1993,58(4):576-588.
[75]Levander A R, Gibson B S. Wide-angle seismic reflections from two dimensional target zones [J]. Geophys. Res.1991.96:10251-10260.
[76]Leary P C. Quantifying crustal fracture heterogeneity by seismic scattering[J]. Geophys. J. Int.1995,122:125-142.
[77]Larkin S P, Levander A. A deterministic and stochastic velocity model for the Salton Trough/Basin and Range transition zone and constraints on magmatism during rifting[J]. Geophys. Res.1996,101:27883-27897.
[78]Holliger K. Upper-crustal velocity heterogeneities as derived from a variety of P-wave sonic logs[J]. Geophys. J. Int.1997,125:813-829.
[79]Wu R, Xu Z, Li X. Heterogeneity spectrum and scale-anisotropy in the upper crust revealed by the German Continental Deep-Drilling (KTB) Holes[J]. Geophys. Res. Lett.,1994,21:911-914.
[80]Levander A, Hobbs R W. The crust as a heterogeneous'optical'medium, or 'crocodiles in the mist'[J]. Tectonophysics,1994,232:281-297.
[81]Kneib G. The statistical nature of the upper continental crystalline crust derived from in seismic measurements[J]. Geophys. J. Int,1995,122:594-616.
[82]Holliger K, Green A G, Juhlin C. Stochastic analysis of sonic logs from the upper crystalline crust methodology[J]. Tectonophysics,1996,264:341-356.
[83]Frenje L, Juhlin C. Scattering of seismic waves simulated by finite difference modeling in random media:Application to the Gravberg-1 well, Sweden[J]. Tectonophysics,1998,293:61-68.
[84]Thomas B. Parallel 3-D viscoelastic finite difference seismic modeling[J]. Computers & Geosciences,2001,28:887-899.
[85]奚先,姚姚.二维随机介质及波动方程正演模拟[J].石油地球物理勘探,2001,36(5)：546-552.
[86]奚先,姚姚.随机介质模型的模拟与混合型随机介质[J].地球科学—中国地质大学学报,2002,27(1)：67-71.
[87]奚先,姚姚.二维粘弹性随机介质中的波场特征分析[J].地球物理学进展,2004,19(3),608-615.
[88]奚先,姚姚.二维横向各向同性弹性随机介质中的波场特征[J].地球物理学进展,2004,19(4)：924-932.
[89]奚先,姚姚.非平稳随机介质模型[J].石油地球物理勘探,2005,40(1)：71-75.
[90]姚姚,奚先.随机介质模型正演模拟及其地震波场分析[J].石油物探,2002,41(1)：31-36.
[91]Cerveny V, Molotkov I A. Ray method in seismology[M]. Karlova University: Praha,1977.
[92]张中杰,何樵登.含裂隙介质中地震波运动学问题的正演模拟[J].石油地球物理勘探,1989,24(3),290-300.
[93]Pao Y H, Varatharajulu V. Huygens' principle, radiation conditions, and integral formulas for the scattering of elastic waves[J]. J Acoust Soc Am,1976,59: 1361-1371.
[94]Bennett C L, Mieras H.J. Acoust Soc.Am.,1981,69(5):1261-1265.
[95]Bouchon M. Diffraction of elastic waves by cracks or cavities using the discrete wavenumber method[J]. Journal of the acoustical society of America.1987, 81:1671-167.
[96]Bouchon M, Campillo M, Gaffet S. A boundary integral equation-discrete wavenumber representation method to study wave propagation in multilayered media having irregular interfaces, Geophysics,1989,54,1134-1140.
[97]Bakamjian B. Boundary integrals:An efficient methodfor modeling seismic wave propagation in 3-D[C].62th Ann.Internat.Mtg.Soc. Expl.Geophys. Tulsa:Univ.of Tulsa,1992:1232-1235.
[98]符力耘,牟永光.弹性波边界元法正演模拟[J].地球物理学报,1994,37(04)：521-529.
[99]闫贫,何樵登.二维横向各向同性介质中的格林函数积分法合成记录[J].物化探计算技术,1995,17(2)：1-6.
[100]Fu L Y,Mu Y QYang H J. Forward problem of nonlinear Fredholm integral equation in reference medium via velocity-weighted wavefield function[J]. Geophysics,1997,62:650-656.
[101]Fu L Y. Numerical study of generalized Lippmann-Schwinger integral equation including surface topography[J]. Geophysics,2003,68:665-671.
[102]Lysmer J, Drake L A. A finite element method for seismology[J]. Methods in Computational Physics,1972:11.
[103]杨宝俊,何樵登.有限元素法的一种频率域计算方法及其应用[J].石油物探,1982,21(1)：56-64.
[104]Marfurt K J. Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations[J]. Geophysics,1984,49(5):533-549.
[105]Seron F J, Sanz F J, Kindelan M. Finite-element method for elastic wave propagation[M]. Comm.Appl.Numerical Methods,1990.
[106]Seron F J, Sabadell F J, Sabadell F J. A numerical laboratory for simulation and visualization of seismic wavefields[J]. Geophysical Prospocting,1996, 44:603-642.
[107]Padovani E, Seriani G, Seriani G. Low-and high-order finite element method:experience in seismic modeling[J]. Journal of Computational Acoustics, 1994,2:371-422.
[108]Sarma G S, Mallick K, Gadhinglajkar V R. Nonreflecting boundary condition in finite-element formulation for an elastic wave equation[J]. Geophysics,1998, 63(3):1006-1023.
[109]牟永光.三维复杂介质地震物理模拟[M].北京：石油工业出版社,2003.
[110]Gazdag J. Modeling of the acoustic wave equation with transform method[J]. Geophysics.1981.46(6):854-859.
[111]Kosloff D, Baysal E. Forward modeling by a Fourier Method[J]. Geophysics, 1982,47(10):1402-1412.
[112]Fornberg B. The pseudospectral method:Comparisons with finite-differences for the elastic wave equations[J]. Geophysics,1987,52(4):483--501.
[113]Reshef M, Kosloff D. Three-dimensional acoustic modeling by the Fourier method[J]. Geophysics,1988,3(9):1175-1183.
[114]刘洋,董敏煜.各向异性介质中的方位AVO[J]石油地球物理勘探,1999,34(3)：260-268.
[115]候安宁,何樵登.各向异性介质中弹性波特征的伪谱法模拟研究[J].石油物探,1995,34(2)：36-45.
[116]张文生,何樵登.二维横向各向同性介质的伪谱法正演模拟[J].石油地球物理勘探,1998,33(3)：310-319.
[117]傅旦丹,何樵登.正交各向异性介质地震弹性波场的伪谱法正演模拟[J].石油物探,2001,40：8-14.
[118]Alterman, Z, Karal F C. Propagation of elastic waves in layered media by finite-difference methods[J]. Bull. Seism. Soc. Am.1968,58:367-398.
[119]Boore D M. Finite-difference methods for seismic wave propagation in heterogeneous materials in Methods in computational physics[M]. B.A.Bolt, ed.,Academic Press,ine.1972:11.
[120]Alford R M, Kelly K R, Boore D M. Accuracy of Finite-difference Modeling of the Acoustic Wave Equation[J]. Geophysics,1974,39:834-842.
[121]Kelly K R, Ward R W, Treitel S, et al. Synthetic seismograms, a finite difference approach[J]. Geophysics,1976,41(1):2-27.
[122]Dablain M A. The application of high2order differencing to the scalar wave equation[J]. Geophysics,1986,51(1):54-66.
[123]Virieux J. SH-wave propagation in heterogeneous media:velocity-stress finite difference method[J]. Geophysics,1984,49(11):1933-1957.
[124]Virieux J. P-SV wave propagation in heterogeneous media, Velocity-stress finite-difference method[J]. Geophysics,1986,51(4):889-901.
[125]Igel H, Mora P, Riollet B. Anisotropic wave propagation through finite difference grids[J]. Geophysics,1995,60(4):1203-1216.
[126]Dong Z, McMechan G A.3-D viscoelasitc anisotropic modeling ofdata from a multi-component,multi-azimuth experiment in northeast Texas[J]. Geophysics, 1995,60(4):1128-1150.
[127]Ramos-Martinez J, Ortega A A. Shear-wave splitting at vertical incidence in media containing intersecting fracture systems[C].68th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts,1998.
[128]Ramos-Martinez J, Ortega A A, McMechan G A.3-D seismic modeling for cracked media:shear-wave splitting at zero offset[J]. Geophysics,2000,65(1): 211-221.
[129]Robertsson Johan O A, Blanch J O, Symes W W. Viscoelastic finite-difference modeling[J]. Geophysics,1994,59(9):1444-1456.
[130]Robertsson J O A. A numerical free-surface condition for elastic/viscoelastic finite-difference modeling in the presence of topography[J]. Geophysics,1996, 61(6):1921-1934.
[131]Graves R W. Simulating seismic wave propagation in 3D elastic media using staggered-grid finite-difference[J].BSSA,1996,86:1091-1106.
[132]侯安宁,何樵登.各向异性介质中弹性波动方程高阶差分及其稳定性研究[J].地球物理学报,1995,38(4)：519-527.
[133]裴正林.三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟[J].石油大学学报(自然科学版),2004,28(5)：23-29.
[134]裴正林.三维双相各向异性介质弹性波方程交错网格高阶有限差分法模拟[J].中国石油大学学报(自然科学版)，2006,30(2)：608-615.
[135]裴正林.层状各向异性介质中横波分裂和再分裂数值模拟[J].石油地球物理勘探,2006,41(1)：17-25.
[136]Jastram C, Tessemer E. Elastic modeling on a grid of vertically varying spacing[J]. Geophysical Prospecting,1994,42(4):357-370.
[137]Falk J, Tessmer E, Gajevski D. Tube wave modeling by the finite-differences method with varying grid spacing[J]. Pure and Applied Geophysics,1996,148(3): 77-93.
[138]Oprsal I, Zahradnik J. Elastic finite-Difference method for irregular grids[J]. Geophysics,1999,64(1):240-250.
[139]Pitarka A.3D Elastic Finite-Difference Modeling of Seismic Motion Using Staggered Grids with Nonuniform Spacing[J]. Bulletin of Seismological Society of America,1999,89 (1):54-68.
[140]Nordstrom J. Carpenter M H. High-order finite difference methods, multidimensional linear problems and curvilinear coordinates[J]. Journal of Computational Physics,2001,173(1):149-174.
[141]常印佛,刘湘培,吴言昌.长江中下游铜铁成矿带[M].北京：地质出版社,1991：71-76.
[142]翟裕生,姚书振,林新多.长江中下游地区铁铜矿床.北京：地质出版社,1992：1-120.
[143]任启江,刘孝善,徐兆文,胡受奚,胡文瑄.安徽庐枞中生代火山构造洼地及其成矿作用.北京：地质出版社,1991：1-206.
[144]周涛发,范裕,袁峰,陆三明,尚世贵,David R Cooke, Sebastien Meffre,赵国春.安徽庐枞(庐江-枞阳)盆地火山岩的年代学及其意义[J].中国科学,2008,38(11)：1342-1353.
[145]安徽省地质局庐枞地区铁矿会战指挥部.安徽省庐江罗河铁矿详细地质勘探报告[R].1980.
[146]黄清涛,尹恭沛.安徽庐江罗河铁矿[R]//中华人民共和国地质矿产部地质专报-矿床与矿产(第14号),安徽省地质矿产局,北京：地质出版社,1989.
[147]安徽省地质局三二七地质队.安徽省庐江县大包庄硫铁矿床详细普查地质报告[R].1981.
[148]安徽省地质调查院.安徽省庐江县泥河铁矿勘探设计书[R].2008.
[149]Holbrook W S, Mooney W D, Christensen N I.1992, The seismic velocity structure of the deep continental crust[M]. In Fountain D M, Arculus R, Kay R. eds., The lower continental crust:Amsterdam, Elsevier,1992:1-43.
[150]Yilmaz O.1987. Seismic data processing, Series:Investigation in Geophysics. Volume 2 [M]. Tulsa:Society of Exploration Geophysicists,1987:1-526.
[151]马中高,解吉高.2005.岩石的纵、横波速度与密度的规律研究[J].地球物理学进展,2005,20(4)：905-910.
[152]吴明安,等.泥河矿区勘探报告[R].2010.
[153]Ergintav S, Canitez N. Modeling of multi-scale media in discrete form[J]. Journal of Seismic Exploration,1997,6:77-96.
[154]雍运动.三维粘弹随机介质地震波场并行模拟研究[D].长春：吉林大学硕士学位论文,2007.
[155]Flatte S M, Wu R S. Small-scale structure in the lithosphere and asthenosphere deduced from arrival time and amplitude fluctuations at NORSAR[J]. Geophys Res,1988,93:6601-6614.
[156]王德利,何樵登,韩立国.裂隙型单斜介质中多方位地面三分量记录模拟[J].地球物理学报,2005,48(2)：386-393.
[157]Levander A R. Fourth-order finite-difference p-sv seismograms [J]. Geophysics. 1988,53(11):1425-143.
[158]王者江.基于BISQ机制的三维双相介质正演模拟及传播特性研究[D].长春：吉林大学博士学位论文，2008.
[159]王恩利.基于各向异性的复合介质弹性波场模拟与特征分析[D].长春：吉林大学博士学位论文,2008.
[160]Fornberg. Pseudospectral approximation of the elastic wave equation on a staggered grid[C].59th Ann. Intemat. Mtg. Soc. Expl.Geophys. Expanded Abstracts,1989:1047-1049.
[161]董良国.一阶弹性波动方程交错网格高阶差分解法[J].地球物理学报,2000，43(3)：411-419.
[162]侯安宁.各向异性弹性波及其波动方程正反演研究[D].长春：长春地质学院博士学位论文,1994.
[163]Cerjan C, Kosloff D. A non-reflection boundary condition for discrete acoustic and elastic wave equation[J]. Geophysics,1985,50:705-708.
[164]伊尔马兹著,刘怀山等译,地震资料分析—地震资料处理、反演和解释[M].北京：石油工业出版社,2006：482-487.
[165]邓志文.复杂山地地震勘探[M].北京：石油工业出版社，2006,70-71.Tessmer G and Behle A. Common reflection point data-stacking technique for converted waves[J]. Geophysical prospecting,1988,36:671-688.
[166]Taner m t, Koehler F. Velocity Spectra-Digital computer derivation and applications of velocity function[J].Geophysics,1969,34960:859-881.
[167]黄绪德,杨文霞.转换波地震勘探[M].北京：石油工来出版社,2008：111-113.
[2]Eaton D W, Milkereit B, Salisbury M. Seismic methods for deep mineral exploration:mature technologies adapted to new targets[J]. The Leading Edge, 2003,22(6):580-585.
[3]Aminzadeh F. Future geophysical technology trends [M]. The Leading Edge, 1996:739-742.
[4]吕庆田,廉玉广,赵金花.反射地震技术在成矿地质背景与深部矿产勘查中的应用：现状与前景[J].地质学报,2010,84(6)：771-787.
[5]Brown L, Barazangi M, Kaufman S, et al. The first decade of COCORP:1974-1984. in Reflection Profiling:A Global Perspective[M]// Barazangi M and Brown L, eds, Geodynamics Series,13:107-120, American Geophysical Union, Washington, D C.1986.
[6]Matthews D H, Smith C., Deep seismic reflection profiling of the continental lithosphere[J]. Geophys. J. Roy. Astr. Soc.,1987,89:209-216
[7]Klemperer S L, Hobbs R. The BIRPS Atlas[M]. edited by Cambridge University Press, Cambridge, UK, hardcover,1992.
[8]John Wilson. LITHOPROBE:Dancing Elephants & Floating Continents-The story of Canada beneath your feet[M]., Published by Key Porter Books,2003.
[9]Australia Geodynamics Cooperative Research Center. Annual Report 1998/1999. 1999.
[10]Cook F, Brown L, Kaufman S, et al. COCORP seismic reflection profiling of the Appalachian orogen beneath the Coastal Plain of Georgia[J]. Geol. Soc. Amer. Bull.,92:738-748.
[11]Brun J P, Wenzel F, ECORS-DEKORP team. Crustal-scale structure of the southern Rhinegraben from ECORS-DEKORP seismic reflection data[J].. Geology,1981,19:758-762
[12]Zhao W J, Nelson K D. Project INDEPTH Team. Deep seismic reflection evidence for continental underthrusting beneath southern Tibet[J]. Nature,1993, 366:557-559.
[13]Nelson K D, Zhao W J. Project INDEPTH Team. Partially molten middle crust beneath Southern Tibet:Synthesis of Project INDEPTH results[J]. Science,1996, 274:1684-1688.
[14]董树文,李廷栋SinoProbe—中国深部探测实验[J].地质学报,2009,83(9)：895-909.
[15]Perron G and Calvert A J. Shallow, high-resolution seismic imaging at the ansil mining camp in the Abitibi greenstone belt[J]. Geophysics,1998, 63(2):379-391.
[16]Milkereit B, Green A, Sudbury Working Group.1992. Deep geometry of the Sudbury structure from seismic reflection profiling[J]. Geology,1992, 20:807-811.
[17]Roy B, Clowes R M. Seismic and potential-field imaging of the Guichon Creek batholith, Brithsh Columia, Canada, to delineate structures hosting porphyry copper deposits[J]. Geophysics,2000,65(5):1418-1434
[18]Pretorius C C, Trewick W F, Fourie A, Irons C. Application of 3-D seismic to mine planning at Vaal Reefs gold mine, number 10 shaft, Republic of South Africa[J]. Geophysics,2000,65 (6):1862-1870
[19]Stuart G W, Jolley S J, Polome L G. Application of 3-D seismic attributes analysis to mine planning:Target gold deposit, South Africa[J]. The Leading Edge 2000,19(7):736-742.
[20]White D, Boerner D, Wu J J, et al. Mineral exploration in the Thompson nickel belt, Manitoba, Canada, using seismic and controlled-source EM methods[J]. Geophysics,2000,65(6):1871-1881.
[21]Milkereit B, Berrer E K, King A R, et al. Development of 3-D seismic exploration technology for deep nickel-copper deposits-A case history from the Sudbury basin, Canada[J]. Geophysics,2000,65(6):1890-1899.
[22]吕庆田,侯增谦,史大年等.铜陵狮子山金属矿地震探测结果及对区域找矿的意义[J].矿床地质,2004,23(3)：390-398.
[23]吕庆田,史大年,赵金花,严加永,徐明才.深部矿产勘查的地震学方法：问题与前景一铜陵矿集区的应用实例[J].地质通报,2005,24(3)：211-218.
[24]吕庆田,韩立国,严加永,廉玉广,史大年,颜庭杰.庐枞矿集区火山气液型铁硫矿床及控矿构造的反射地震成像[J].岩石学报,2010,26(9)：2598-2612.
[25]徐明才,高景华,荣立新,等.地面地震层析技术在金属矿勘查中的试验研究[J].物探与化探,2005,29(4)：299-303.
[26]史大年,吕庆田,徐明才.铜陵矿集区地壳浅表结构的地震层析研究[J].矿床地质,2004,23(3)：383-389.
[27]Li T L, Eaton D W. Delineating the Tuwu porphyry copper deposit at Xinjiang, China, with seismic-reflection profiling[J]. Geophysics.2005,70(6):B53-B60.
[28]梁光河,蔡新平.地震勘探在山东蓬家夼金矿深部预测中的应用[J].矿产与地质,2001,15(6)：743-748.
[29]徐明才,高景华,柴铭涛,等.金属矿地震勘查的方法技术[J].岩土工程界,1997,6(4)：41-46.
[30]徐明才,高景华,柴铭涛,等.寻找隐伏金属矿的地震方法技术研究[J].物探与化探,1997,21(6)：69-75.
[31]徐明才,高景华.从金属矿地震方法的试验效果探讨其应用前景[J].中国地质,2004,31(1)：110-114.
[32]徐明才,高景华,荣立新,等.散射波地震方法在蔡家营多金属矿区的试验研究[J].物探与化探,2003,27(1)：52-57.
[33]高景华,徐明才,荣立新,等.小热泉子铜矿区地震方法试验研究[J].地质与勘探，2004,40(6)：49-54.
[34]勾丽敏,刘学伟,雷鹏,等.金属矿地震勘探技术方法研究综述—散射波地震勘探方法[J].勘探地球物理进展,2007，30(2)：65-90.
[35]Mitchell, A H G, Garson M S. Mineral deposits and global tectonic settings[M]. New York:Academic Press,1981:1-405.
[36]Sawkins F J. Metal deposits in relation to plate tectonics[M]. Berlin:Springer Verlag,1984:1-325.
[37]Kerrich R, Goldfarb R J, Richards J P. Metallogenic Provinces in an Evolving Geodynamic Framework[J]. Economic Geology 100th anniversary Volume, 2005:1097-1136
[38]Beaumont C, Quinlan G. A geodynamic framework for interpreting crustal-scale seismic-reflectivity patterns in compressional orogens[J]. Geophysics Journal International,1994,116:754-783.
[39]Drummond B J, Goleby B R, Owen A J, Yeates A N, Swager C, Zhang Y, and Jackson J K. Seismic reflection imaging of mineral systems:Three case histories[J]. Geophysics,2000,65(6):1852-1861.
[40]Oliver N H S. Linking of regional and local hydrothermal systems in the mid-crust by shearing and faulting[J]. Tectonophysics,2001,335:147-161.
[41]Goleby B R, Blewett R S, Korsch R J, et al. Deep seismic reflection profiling in the Archaean northeastern Yilgarn Craton, Western Australia:Implication for crustal architecture and mineral potential [J]. Tectonophysics,2004,388: 119-133.
[42]Drummond B J, Goleby B R. Seismic reflection images of major ore-controlling structure in the Eastern Goldfields,Western Australia[J]. Expl. Geophys.,1993, 24:473-478.
[43]Swager C P, Goleby B R, Drummond B J, Rattenbury M S, Williams P R. Crustal structure of granite-greenstone terranes in the Eastern Goldfields, Yilgarn Craton, as revealed by seismic reflection profiling[J]. Precambrian Research,1997, 83(1-3):43-56.
[44]Goleby B R, Rattenbury M S, Swager, Swager, C P, et al. Archaean crustal structure from seismic reflection profiling, Eastern Goldfields, Western Australia[C]. Austral. Geol. Surv. Org. record,1993/15.
[45]Upton P, Hobbs B, Ord A, Zhang Y, Drummond B, Archibald N. Thermal and deformation modelling of the Yilgarn deep seismic transect[C]. Geodynamics and Ore Deposits Conference, Australia Geodynamics Cooperative Research Centre, Ballarat, Victoria,1997:22-25.
[46]Sorjonen-Ward P, Zhang Y, Zhao C. Numerical modelling of orogenic processes and gold mineralisation in the southeastern part of the Yilgarn Craton, Western Australia[J]. Australian Journal of Earth Sciences,2002,49:935-964.
[47]Jones T, Nur A. The Nature of seismic reflection from deep crustal fault zones[J]. Geophys. Res.,1984,89:3153-3171
[48]Malehmir A, Tryggvason A, Lickorish H, Weihed P. Regional structural profiles in the western part of the Palaeoproterozoic Skellefte ore district, northern Sweden[J]. Precambrin Research,2007,159:1-18.
[49]Tryggvason A, Malehmir A, Rodriguez-Tablante J, et al. Reflection seismic investigation in the western part of the palaeoproterozoic VHMS bearing Skellefte Districts, northern Sweden[J]. Economic Geology,2006, 101:1039-1054.
[50]Salisbury M H, Milkereit B, and Bleeker W. Seismic imaging of massive sulfide deposits:Part I. rock properies[J]. Economic Geology,1996,91(5):821-828.
[51]Salibury M H, Harvey C W, and Matthews. The acoustic properties of ores and host rocks in hardrock terranes[M]//In D. W. Eaton., B. Milkereit, M. H. Salisbury, Eds., Hardrock Seismic Exploration,2003.
[52]Bohlen T, Muller C, Milkereit B.2003. Elastic seismic-wave scattering from massive sulfide orebodies:on the role of composition and shape[M]//In D. W. Eaton., B. Milkereit, M. H. Salisbury, Eds., Hardrock Seismic Exploration, Geophysical developments.2003,10:70-90.
[53]Clarke G J, Eaton D W. Influence of morphology and surface roughness on the seismic response of massive sulfides, based on elastic-wave kirchhoff modeling[M]. In D. W. Eaton., B. Milkereit, M. H. Salisbury, Eds., Hardrock Seismic Exploration, Geophysical developments,2003,10:45-58.
[54]Eaton D E. Weak elastic scatting from massive sulfide ore bodies[J]. Geophysics. 1999,64(1):289-299.
[55]Hobbs R W.2003.3D modelling of seismic-wave propagation using complex elastic screen, with application to mineral exploration[M]//In D. W. Eaton., B. Milkereit, M. H. Salisbury, Eds., Hardrock Seismic Exploration, Geophysical developments,2003,10:59-69.
[56]Adame E, Milkereit B, and Salmon B.2004.3D seismic exploration in the Val d'Or mining camp, Quebec[C]. SEG Int'l Exposition and 74th Annual Meeting, Denver, Colorado,10-15 October,2004.
[57]Calvert A J, and Li Y.Seismic reflection imaging over a massive sulfide deposit at the Matagami mining camp, Quebec[J]. Geophysics,1999,64(l):24-32.
[58]Adam E, Perron G, Arnold G, et al.3D seismic imaging for VMS deposit exploration, Matagami, Quebec[M]//In D. W. Eaton., B. Milkereit, M. H. Salisbury, Eds., Hardrock Seismic Exploration, Geophysical developments,2003, 10:229-246.
[59]Milkereit B, Berrer E K, King A R, et al. Development of 3D seismic exploration technology for deep nickel-copper deposits-A case history from the Sudbury basin, Canada[J]. Geophysics,2000,65(6):1890-1899.
[60]Salisbury M H, Milkereit B, Ascough G, et al. Physical properties and seismic imaging of massive sulfides[J]. Geophysics,2000,65(6):1882-1889.
[61]Matthews L. Base metal exploration:Looking deeper and adding value with seismic data[C]. Canadian Society of Exploration Geophysicists Recorder, 2002:37-43.
[62]Bellefleur G, Muller C, Snyder D, Matthews L. Downhole seismic imaging of a massive sulfide orebody with mode-converted waves, Halfmile Lake, New Brunswick, Canada[J]. Geophysics,2004,69(2):318-329.
[63]Fountain D, Salisbury M H, Percival J. Seismic structure of the continental crust based on rock velocity measurements from the Kapuskasing uplift[J]. J. Geophys. Res.,1990,95:1167-1186.
[64]Pretorius C C, Jamison A A, Irons C. Seismic exploration in the Witwatersrand basin, Republic of South Africa[J]. Exploration 87, Proceedings,1987,3: 241-253.
[65]Roberts B. Zaleski E, Perron G et al. Seismic Exploration of the Manitouwadge Greenstone Belt, Ontario:A case history [M]//In D.W. Eaton.,B.Milkereit, M. H. Salisbury, Eds.,Hardrock Seismic Exploration, Geophysical developments,2003, 10:110-126.
[66]Aki K. Analysis of the Seismic Coda of Local Earthquakes as Scattered Waves[J]. Geophysics,1969,74(4):615-631.
[67]Godfrey R, Muir F, Rocca F. Modeling seismic impedance with Markovchains[J]. Geophysics,1980,45:1351-1372.
[68]Frankel A, Clayton R W, A finite-difference simulation of wave propagation in two-dimensional random media[M].Bulletin of the Seismological Society of America,1984,74(6):2167-2186.
[69]Frankel A, Clayton R W. Finite-difference simulation of seismic scattering: Implications for the propagation of short-period seismic waves in the crust and models of crustal heterogeneity [J]. Journal of Geophysical Research,1986,91: 6465-6489.
[70]Stanke F E, Burridge R. Comparison of spatial and ensemble averaging methods applied to wave propagation in finely layered media[C].60th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts,1990:1062-1065.
[71]De Hoop M V, Bumidge R, Chang H W. Wave propagation with tunnelingin a highly discontinuous layered medium[J]. Wave Motion,1991,13:307-329.
[72]Kerner C. Anisotropy in sedimentary rocks modeled as random media[J]. Geophysics,1992,57:564-576.
[73]Ikelle L T, Yung S K.2-D random media with ellipsoidal autocorrelation functions[J]. Geophysics,1993,58(9):1359-1372.
[74]Kneib G, Kerner C. A ccurate and efficient seismic modeling in random media[J]. Geophysics,1993,58(4):576-588.
[75]Levander A R, Gibson B S. Wide-angle seismic reflections from two dimensional target zones [J]. Geophys. Res.1991.96:10251-10260.
[76]Leary P C. Quantifying crustal fracture heterogeneity by seismic scattering[J]. Geophys. J. Int.1995,122:125-142.
[77]Larkin S P, Levander A. A deterministic and stochastic velocity model for the Salton Trough/Basin and Range transition zone and constraints on magmatism during rifting[J]. Geophys. Res.1996,101:27883-27897.
[78]Holliger K. Upper-crustal velocity heterogeneities as derived from a variety of P-wave sonic logs[J]. Geophys. J. Int.1997,125:813-829.
[79]Wu R, Xu Z, Li X. Heterogeneity spectrum and scale-anisotropy in the upper crust revealed by the German Continental Deep-Drilling (KTB) Holes[J]. Geophys. Res. Lett.,1994,21:911-914.
[80]Levander A, Hobbs R W. The crust as a heterogeneous'optical'medium, or 'crocodiles in the mist'[J]. Tectonophysics,1994,232:281-297.
[81]Kneib G. The statistical nature of the upper continental crystalline crust derived from in seismic measurements[J]. Geophys. J. Int,1995,122:594-616.
[82]Holliger K, Green A G, Juhlin C. Stochastic analysis of sonic logs from the upper crystalline crust methodology[J]. Tectonophysics,1996,264:341-356.
[83]Frenje L, Juhlin C. Scattering of seismic waves simulated by finite difference modeling in random media:Application to the Gravberg-1 well, Sweden[J]. Tectonophysics,1998,293:61-68.
[84]Thomas B. Parallel 3-D viscoelastic finite difference seismic modeling[J]. Computers & Geosciences,2001,28:887-899.
[85]奚先,姚姚.二维随机介质及波动方程正演模拟[J].石油地球物理勘探,2001,36(5)：546-552.
[86]奚先,姚姚.随机介质模型的模拟与混合型随机介质[J].地球科学—中国地质大学学报,2002,27(1)：67-71.
[87]奚先,姚姚.二维粘弹性随机介质中的波场特征分析[J].地球物理学进展,2004,19(3),608-615.
[88]奚先,姚姚.二维横向各向同性弹性随机介质中的波场特征[J].地球物理学进展,2004,19(4)：924-932.
[89]奚先,姚姚.非平稳随机介质模型[J].石油地球物理勘探,2005,40(1)：71-75.
[90]姚姚,奚先.随机介质模型正演模拟及其地震波场分析[J].石油物探,2002,41(1)：31-36.
[91]Cerveny V, Molotkov I A. Ray method in seismology[M]. Karlova University: Praha,1977.
[92]张中杰,何樵登.含裂隙介质中地震波运动学问题的正演模拟[J].石油地球物理勘探,1989,24(3),290-300.
[93]Pao Y H, Varatharajulu V. Huygens' principle, radiation conditions, and integral formulas for the scattering of elastic waves[J]. J Acoust Soc Am,1976,59: 1361-1371.
[94]Bennett C L, Mieras H.J. Acoust Soc.Am.,1981,69(5):1261-1265.
[95]Bouchon M. Diffraction of elastic waves by cracks or cavities using the discrete wavenumber method[J]. Journal of the acoustical society of America.1987, 81:1671-167.
[96]Bouchon M, Campillo M, Gaffet S. A boundary integral equation-discrete wavenumber representation method to study wave propagation in multilayered media having irregular interfaces, Geophysics,1989,54,1134-1140.
[97]Bakamjian B. Boundary integrals:An efficient methodfor modeling seismic wave propagation in 3-D[C].62th Ann.Internat.Mtg.Soc. Expl.Geophys. Tulsa:Univ.of Tulsa,1992:1232-1235.
[98]符力耘,牟永光.弹性波边界元法正演模拟[J].地球物理学报,1994,37(04)：521-529.
[99]闫贫,何樵登.二维横向各向同性介质中的格林函数积分法合成记录[J].物化探计算技术,1995,17(2)：1-6.
[100]Fu L Y,Mu Y QYang H J. Forward problem of nonlinear Fredholm integral equation in reference medium via velocity-weighted wavefield function[J]. Geophysics,1997,62:650-656.
[101]Fu L Y. Numerical study of generalized Lippmann-Schwinger integral equation including surface topography[J]. Geophysics,2003,68:665-671.
[102]Lysmer J, Drake L A. A finite element method for seismology[J]. Methods in Computational Physics,1972:11.
[103]杨宝俊,何樵登.有限元素法的一种频率域计算方法及其应用[J].石油物探,1982,21(1)：56-64.
[104]Marfurt K J. Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations[J]. Geophysics,1984,49(5):533-549.
[105]Seron F J, Sanz F J, Kindelan M. Finite-element method for elastic wave propagation[M]. Comm.Appl.Numerical Methods,1990.
[106]Seron F J, Sabadell F J, Sabadell F J. A numerical laboratory for simulation and visualization of seismic wavefields[J]. Geophysical Prospocting,1996, 44:603-642.
[107]Padovani E, Seriani G, Seriani G. Low-and high-order finite element method:experience in seismic modeling[J]. Journal of Computational Acoustics, 1994,2:371-422.
[108]Sarma G S, Mallick K, Gadhinglajkar V R. Nonreflecting boundary condition in finite-element formulation for an elastic wave equation[J]. Geophysics,1998, 63(3):1006-1023.
[109]牟永光.三维复杂介质地震物理模拟[M].北京：石油工业出版社,2003.
[110]Gazdag J. Modeling of the acoustic wave equation with transform method[J]. Geophysics.1981.46(6):854-859.
[111]Kosloff D, Baysal E. Forward modeling by a Fourier Method[J]. Geophysics, 1982,47(10):1402-1412.
[112]Fornberg B. The pseudospectral method:Comparisons with finite-differences for the elastic wave equations[J]. Geophysics,1987,52(4):483--501.
[113]Reshef M, Kosloff D. Three-dimensional acoustic modeling by the Fourier method[J]. Geophysics,1988,3(9):1175-1183.
[114]刘洋,董敏煜.各向异性介质中的方位AVO[J]石油地球物理勘探,1999,34(3)：260-268.
[115]候安宁,何樵登.各向异性介质中弹性波特征的伪谱法模拟研究[J].石油物探,1995,34(2)：36-45.
[116]张文生,何樵登.二维横向各向同性介质的伪谱法正演模拟[J].石油地球物理勘探,1998,33(3)：310-319.
[117]傅旦丹,何樵登.正交各向异性介质地震弹性波场的伪谱法正演模拟[J].石油物探,2001,40：8-14.
[118]Alterman, Z, Karal F C. Propagation of elastic waves in layered media by finite-difference methods[J]. Bull. Seism. Soc. Am.1968,58:367-398.
[119]Boore D M. Finite-difference methods for seismic wave propagation in heterogeneous materials in Methods in computational physics[M]. B.A.Bolt, ed.,Academic Press,ine.1972:11.
[120]Alford R M, Kelly K R, Boore D M. Accuracy of Finite-difference Modeling of the Acoustic Wave Equation[J]. Geophysics,1974,39:834-842.
[121]Kelly K R, Ward R W, Treitel S, et al. Synthetic seismograms, a finite difference approach[J]. Geophysics,1976,41(1):2-27.
[122]Dablain M A. The application of high2order differencing to the scalar wave equation[J]. Geophysics,1986,51(1):54-66.
[123]Virieux J. SH-wave propagation in heterogeneous media:velocity-stress finite difference method[J]. Geophysics,1984,49(11):1933-1957.
[124]Virieux J. P-SV wave propagation in heterogeneous media, Velocity-stress finite-difference method[J]. Geophysics,1986,51(4):889-901.
[125]Igel H, Mora P, Riollet B. Anisotropic wave propagation through finite difference grids[J]. Geophysics,1995,60(4):1203-1216.
[126]Dong Z, McMechan G A.3-D viscoelasitc anisotropic modeling ofdata from a multi-component,multi-azimuth experiment in northeast Texas[J]. Geophysics, 1995,60(4):1128-1150.
[127]Ramos-Martinez J, Ortega A A. Shear-wave splitting at vertical incidence in media containing intersecting fracture systems[C].68th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts,1998.
[128]Ramos-Martinez J, Ortega A A, McMechan G A.3-D seismic modeling for cracked media:shear-wave splitting at zero offset[J]. Geophysics,2000,65(1): 211-221.
[129]Robertsson Johan O A, Blanch J O, Symes W W. Viscoelastic finite-difference modeling[J]. Geophysics,1994,59(9):1444-1456.
[130]Robertsson J O A. A numerical free-surface condition for elastic/viscoelastic finite-difference modeling in the presence of topography[J]. Geophysics,1996, 61(6):1921-1934.
[131]Graves R W. Simulating seismic wave propagation in 3D elastic media using staggered-grid finite-difference[J].BSSA,1996,86:1091-1106.
[132]侯安宁,何樵登.各向异性介质中弹性波动方程高阶差分及其稳定性研究[J].地球物理学报,1995,38(4)：519-527.
[133]裴正林.三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟[J].石油大学学报(自然科学版),2004,28(5)：23-29.
[134]裴正林.三维双相各向异性介质弹性波方程交错网格高阶有限差分法模拟[J].中国石油大学学报(自然科学版)，2006,30(2)：608-615.
[135]裴正林.层状各向异性介质中横波分裂和再分裂数值模拟[J].石油地球物理勘探,2006,41(1)：17-25.
[136]Jastram C, Tessemer E. Elastic modeling on a grid of vertically varying spacing[J]. Geophysical Prospecting,1994,42(4):357-370.
[137]Falk J, Tessmer E, Gajevski D. Tube wave modeling by the finite-differences method with varying grid spacing[J]. Pure and Applied Geophysics,1996,148(3): 77-93.
[138]Oprsal I, Zahradnik J. Elastic finite-Difference method for irregular grids[J]. Geophysics,1999,64(1):240-250.
[139]Pitarka A.3D Elastic Finite-Difference Modeling of Seismic Motion Using Staggered Grids with Nonuniform Spacing[J]. Bulletin of Seismological Society of America,1999,89 (1):54-68.
[140]Nordstrom J. Carpenter M H. High-order finite difference methods, multidimensional linear problems and curvilinear coordinates[J]. Journal of Computational Physics,2001,173(1):149-174.
[141]常印佛,刘湘培,吴言昌.长江中下游铜铁成矿带[M].北京：地质出版社,1991：71-76.
[142]翟裕生,姚书振,林新多.长江中下游地区铁铜矿床.北京：地质出版社,1992：1-120.
[143]任启江,刘孝善,徐兆文,胡受奚,胡文瑄.安徽庐枞中生代火山构造洼地及其成矿作用.北京：地质出版社,1991：1-206.
[144]周涛发,范裕,袁峰,陆三明,尚世贵,David R Cooke, Sebastien Meffre,赵国春.安徽庐枞(庐江-枞阳)盆地火山岩的年代学及其意义[J].中国科学,2008,38(11)：1342-1353.
[145]安徽省地质局庐枞地区铁矿会战指挥部.安徽省庐江罗河铁矿详细地质勘探报告[R].1980.
[146]黄清涛,尹恭沛.安徽庐江罗河铁矿[R]//中华人民共和国地质矿产部地质专报-矿床与矿产(第14号),安徽省地质矿产局,北京：地质出版社,1989.
[147]安徽省地质局三二七地质队.安徽省庐江县大包庄硫铁矿床详细普查地质报告[R].1981.
[148]安徽省地质调查院.安徽省庐江县泥河铁矿勘探设计书[R].2008.
[149]Holbrook W S, Mooney W D, Christensen N I.1992, The seismic velocity structure of the deep continental crust[M]. In Fountain D M, Arculus R, Kay R. eds., The lower continental crust:Amsterdam, Elsevier,1992:1-43.
[150]Yilmaz O.1987. Seismic data processing, Series:Investigation in Geophysics. Volume 2 [M]. Tulsa:Society of Exploration Geophysicists,1987:1-526.
[151]马中高,解吉高.2005.岩石的纵、横波速度与密度的规律研究[J].地球物理学进展,2005,20(4)：905-910.
[152]吴明安,等.泥河矿区勘探报告[R].2010.
[153]Ergintav S, Canitez N. Modeling of multi-scale media in discrete form[J]. Journal of Seismic Exploration,1997,6:77-96.
[154]雍运动.三维粘弹随机介质地震波场并行模拟研究[D].长春：吉林大学硕士学位论文,2007.
[155]Flatte S M, Wu R S. Small-scale structure in the lithosphere and asthenosphere deduced from arrival time and amplitude fluctuations at NORSAR[J]. Geophys Res,1988,93:6601-6614.
[156]王德利,何樵登,韩立国.裂隙型单斜介质中多方位地面三分量记录模拟[J].地球物理学报,2005,48(2)：386-393.
[157]Levander A R. Fourth-order finite-difference p-sv seismograms [J]. Geophysics. 1988,53(11):1425-143.
[158]王者江.基于BISQ机制的三维双相介质正演模拟及传播特性研究[D].长春：吉林大学博士学位论文，2008.
[159]王恩利.基于各向异性的复合介质弹性波场模拟与特征分析[D].长春：吉林大学博士学位论文,2008.
[160]Fornberg. Pseudospectral approximation of the elastic wave equation on a staggered grid[C].59th Ann. Intemat. Mtg. Soc. Expl.Geophys. Expanded Abstracts,1989:1047-1049.
[161]董良国.一阶弹性波动方程交错网格高阶差分解法[J].地球物理学报,2000，43(3)：411-419.
[162]侯安宁.各向异性弹性波及其波动方程正反演研究[D].长春：长春地质学院博士学位论文,1994.
[163]Cerjan C, Kosloff D. A non-reflection boundary condition for discrete acoustic and elastic wave equation[J]. Geophysics,1985,50:705-708.
[164]伊尔马兹著,刘怀山等译,地震资料分析—地震资料处理、反演和解释[M].北京：石油工业出版社,2006：482-487.
[165]邓志文.复杂山地地震勘探[M].北京：石油工业出版社，2006,70-71.Tessmer G and Behle A. Common reflection point data-stacking technique for converted waves[J]. Geophysical prospecting,1988,36:671-688.
[166]Taner m t, Koehler F. Velocity Spectra-Digital computer derivation and applications of velocity function[J].Geophysics,1969,34960:859-881.
[167]黄绪德,杨文霞.转换波地震勘探[M].北京：石油工来出版社,2008：111-113.