基于地震干涉法的波场成像技术及应用研究
|
摘要
地震干涉法的提出及其在野外实际数据上的应用，极大的丰富了对地震波传播规律的认识和研究，被认为是地球物理发展中的一个重大进展。地震干涉法的适用范围从一维的层状均匀声波介质推广到三维弹性介质、任意非均匀介质或各向异性介质，甚至可被推广到运动的介质以及衰减介质。地震干涉法的基本思想是从混沌无序的地震信号中提取出有用信号，通过对记录到的地震信号进行干涉，得到新的地震信号，新形成的地震信号不仅包含了原始地震信号的特性，而且能反映出原始信号所不具有的某些重要特征。地震波干涉方法分为主动源地震干涉和被动源地震干涉，其震源可以是确定性震源也可以是噪声震源。主动源地震干涉法为地震勘探提供了一种新的处理方法，可以用于多种数据类型。以垂直地震剖面（VSP）数据类型的应用为代表，通过地震干涉法可以将VSP数据转换成其他类型的数据，如单井地震剖面（SWP）数据，地表地震剖面（SSP）数据等。被动源地震干涉法用于从野外地震噪声中提取出有用信号来获取地震波穿过介质的地质特性并推断地下地质构造的特征。本文中，我们分别就地震干涉法在主动源地震数据和被动源地震数据上的应用进行深入的研究和讨论。
首先阐述了地震干涉法的基本原理，以互易定理为基础，按照无损任意非均匀介质的假设，推导准确的声学和弹性动力学格林函数的表述。而重构的格林函数中不仅包含两个检波点之间的直达波场，也包含所有一次反射和多次散射。对于地震干涉提取格林函数的表述是基于用于相关处理中的检波器是被一个封闭的震源面所包围的这样一个假设。一般情况下，往往不能满足地震干涉应用中所需的震源分布条件，因此文中通过数值模拟实验调查当理论条件不能被满足，且震源位于地下深部位置处时，震源分布特性对重构的反射响应的影响进行讨论。在地震勘探的实际应用中，通常不能满足理论上所需的连续震源分布，但在干涉计算过程中，震源间隔应尽量布设的密集一些，以防止空间假频的产生。当震源只集中分布于地下局部区域内时，重构的道集中只有部分反射同相轴能被恢复出来，且大大降低了重构响应的信噪比。一般而言，可以利用地震干涉法从地下的透射响应中重构出有效的反射信息，并对重构的反射信息进一步应用偏移处理，偏移成像可以改善重构共炮点道集中反射信息的不完整性，获得令人满意的结果。
文中对地震干涉法在主动源数据上的应用进行了深入系统的研究，主要针对的是VSP型数据。当常规地震勘探数据处理方法不能对上覆盖层复杂变化的地震资料获得满意的成像效果时，可以通过虚拟震源方法将VSP数据转换成SWP数据，变换之后，原来位于地表的震源被重建到井下位置处，震源-检波器的井下排列离地质目标体更近，且可以有效的避免近地表变化和复杂上覆盖对数据的影响，不需要知道上覆盖层的复杂速度模型，为处理这种复杂地震资料提供了一种有效途径。本文主要利用虚拟震源方法探测地下陡构造地质体，通过数值模拟发现地质体越陡峭，获得的成像效果越好，并且虚拟震源方法可以自动校正静校正量。虚拟震源方法对于弹性波数据的应用也是有效的，能够获得合理的成像结果。此外，模拟实验表明，虚拟震源方法应用于稀疏观测系统仍然是非常有效的，尽管随着震源及检波器间隔的增大，信噪比随之下降，但是仍能构建出正确的成像结果。
标准VSP数据的成像精度虽然很高，但是成像范围却很有限，仅能对从第一个检波器开始向下的一个三角锥形区域内的地质体进行成像。为了扩展地下地质构造的成像范围，可以利用地震干涉中的相关算法将VSP数据重构为虚拟SSP数据，并对虚拟SSP数据进行成像处理。由此得到的结果中，地下地质构造的成像范围等同于具有相同震源覆盖范围的常规地表地震勘探中获得的成像范围。此外，该变换校正了井中的静校正量，并且不需要知道井中震源激发时间，震源以及检波器的位置信息。文中利用VSP-SSP相关变换将正演模拟的VSP数据转换成虚拟SSP数据，并对重构的地表地震数据进行常规处理，获得地下地质构造的成像。从结果中看出，和标准VSP数据的成像范围相比，地下的成像范围被大大扩展了。通过数值模拟测试了观测系统对VSP-SSP相关变换方法所得结果的影响。当采用井下不同分布范围上局部检波器所获得的响应时，应用VSP-SSP相关变换方法重构的SSP道集受到不同程度的影响，部分反射同相轴不能被正确的重构出来，道集的信噪比也降低了。但是采用稀疏的观测系统时，VSP-SSP相关变换方法依然能够重构出质量较好的SSP道集。当需要对采集自同一地区的VSP和SSP数据进行整合时，可以利用VSP-SSP褶积变换将VSP数据转换成虚拟SSP数据。形成的虚拟SSP数据可以和实际SSP数据进行对比分析，用于确定地下反射体的岩性。而且在该变换过程中也不需要知道速度模型和VSP测量中检波器的位置信息。我们利用VSP-SSP褶积变换将正演模拟的VSP数据转换成虚拟SSP数据，并对重构的地表地震数据进行常规处理，获得地下地质构造的成像。与通过相关算法从VSP数据中重构出来的虚拟SSP数据相比，褶积算法重构的虚拟SSP道集的信噪比更高，对地下较深处反射体的分辨率更好，但是成像范围却很有限。同样通过数值模拟测试了观测系统对VSP-SSP褶积变换方法所得结果的影响。采用仅分布于井下局部范围上检波器所获得的响应时，应用VSP-SSP褶积变换方法重构的SSP道集受到不同程度的影响，位于浅层的反射同相轴不能被重构出来，但依然保持了较高的信噪比。对于稀疏的观测系统，VSP-SSP褶积变换方法依然能够重构出信噪比较高的虚拟SSP道集。将VSP-SSP相关变换和VSP-SSP褶积变换相联合，能够将两种变换方法的优点结合到一起，对构建出的虚拟SSP道集进行常规成像处理，所得结果的成像范围和相关算法中的一样，大大扩展了地下地质构造的成像范围，而成像精度和褶积算法中的一样，对来自于深部的地质体也能分辨出来，且具有较高的信噪比。
被动源地震数据主要采集自地下的微震，记录的背景噪声信号通常能量较弱，信噪比较低。直接对被动源数据进行成像处理，研究储层的地质构造有一定的技术难度。利用地震干涉法从背景噪声数据中提取出新的地震数据，能够把无序的噪声变成有用的信号，将重构的面波响应和反射波信息用于推断地下浅层和深部的地质构造特征。本文通过数值模拟和野外实际数据深入研究了地震干涉法在被动源数据上的应用。首先专门针对被动源地震数据的特点，设计了适合的正演模拟方案，并模拟生成了被动源地震数据。然后对采集自德国Ketzin地区的CO2冲注试验区记录时长将近25小时的背景噪声进行被动源地震干涉技术的应用研究。根据该工区内两条地震测线上所有可用检波器接收到的背景噪声记录，重构出共炮点道集。在重构的道集中可以观察到面波、折射波和反射波信息，与该地区主动源地震测量产生的共炮点道集相比较，验证了虚拟共炮点道集中重构的面波响应和反射波响应的有效性。通过产生共偏移距叠加道集提高虚拟共炮点道集的信噪比。利用重构的面波响应的分散特性估计出面波分散曲线，通过分散曲线的反演获得该地区地下的1D和2D剪切波速度剖面，获得的速度模型与根据该地区主动源数据的旅行时层析反演得到的速度模型有很好的一致性，但是从被动源数据中获得的速度模型对地表附近的浅层具有更高的分辨率。此外，与该地区其他主动源资料中获得的速度剖面相比，也有很好的匹配性。对重构的反射信息进行处理，产生的成像结果能够反映地下地质构造的特征，与对相同测线记录的主动源数据进行处理得到的成像结果有很好的一致性，为被动源地震勘探的实际应用提供了有力的技术支持。
Seismic interferometry method and its application in the field data, greatly enrichedthe understanding and study of seismic wave propagation, it is considered to be a majoradvance in the development of geophysics. The scope of application of seismicinterferometry method is from the one-dimensinal layered acoustic medium to thethree-dimensional elastic medium, arbitrary inhomogeneous media or anisotropic media,and even can be extended to the moving media and attenuation media. The basic idea ofseismic interferometry method is retrieval of useful signal from the chaotic disorderseismic signal, new seismic signal is generated by the interference of recorded seismicsignal, the new seismic signal not only include the properties of original seismic signal, butalso retrieve some important features that the original signal does not have. Seismicinterferometry method can be used in both passive seismic measurements and controlledsource seismic measurements, the source can be certain source or noise source. Activesource seismic interferometry provide a new approach for seismic exploration, it can beused for a variety of data types. Take application of vertical seismic profile (VSP) data asrepresentative, VSP data can be transferred into the other data type by seismicinterferometry method, such as single seismic profile (SWP) data, surface seismic profile(SSP). Passive source seismic interferometry method applied to retrieval of useful signalfrom ambient seismic noise to obtain geological properties of medium, and then infer thecharacteristics of subsurface structure. In this paper, we respectively study and discussionthe application of seismic interferometry method for active source seismic data and passivesource seismic data.
First of all, describe the basic principle of seismic interferometry methods, based onreciprocity theorems, underlying assumption of inhomogeneous, lossless medium, acousticand elastic Green’s function representations can be derived. The reconstructed Green’sfunction not only include the direct wave between two receivers, but also include primaryreflection and multiple scattering. Retrieved Green’s function by seismic interferometry isbased on the assumption that receivers are surrounded by a closed source surface. Undernormao circumstances, the source distribution conditions required by seismicinterferometry can not be met. So in this paper, we investigate that when the circumstancethat theory condition can not be met, the sources distributed in the subsurface, how the properties of source distribution effect on the retrieval of reflected waves. In the practicalapplication of seismic exploration, it always hard to meet the theoretically required for acontinuois source distribution, but during the calculation, in order to avoid aliasing, thesource interval should be set intensively. When the source only distributed in the localsubsurface area, only part of reflected event can be reconstructed in the retrieved gathers,the signal to noise ration is decreased seriously. In general, retrieved reflection informationfrom subsurface transmission response by seismic interferometry method is very effective.Process the retrieved reflection information and then used to migration process, theresulted image can be improve the incomplete reflected waves in the retrieved commonsource gathers and obtain satisfactory results.
In the paper, we study on application of seismic interferometry method for activesource seismic data, here mainly target at VSP data. When the conventional seismic dataprocessing methods failed to acquire the satisfactory imaging results of seismic datainclude overburden complexity, we could apply the virtual source method to convert thetraditional vertical seismic profile gathers into virtual single well profile gathers. After that,redatumed geometry of source-receiver arrays is closer to the structure target and belowthe complex overburden so that we can get the better image resolution of the target afterconventional process and avoid the effected by unknown velocities of complex overburden,provide an efficient way to handle this complex seismic data. In the paper, the virtualsource method applied to detect steep structure such as faults present in the subsurface.Though the synthetic data found that the steeper structure the better image results, andautomatically correct the statics. The virtual source method is still work well for elasticdata sets. Moreover, modeling experiment show that virtual source method is still effectivefor the sparse geometry, we could get the good imaging results. Although the signal tonoise ratio decreased as source interval and receiver interval increased, the image result iscorrect.
Although the imaging accuracy of the standard VSP data is high, the imaging range isvery limited that is defined by a small triangle, with the triangle apex at the shallowestreceiver. In order expand the imaging range of subsurface structure, one can transformVSP data into virtual SSP data by correlation algorithm of seismic interferometry, and thenthe imaging process used to virtual SSP data. The imaging range of subsurface structurefrom the results has same area as that of the original VSP source distribution along thesurface. In addition, this transform eliminates well statics and the need to know the sourceor receiver position in the well and excitation time of source in the well. In this paper,VSP-SSP correlation transformation effectively transform synthetic VSP data into virtual SSP data, and then migrated by a standard surface seismic data processing to give ainterferometric image of subsurface structure. From the results, compared with standardVSP image, VSP-SSP correlation transformation is a huge increase in the subsurfaceillumination. Test on the effect of geometry on the VSP-SSP correlation transformation bynumerical simulation. When the receivers in the well distributed on local area,reconstructed SSP gathers by VSP-SSP correlation transformation with various degrees ofeffect, part of reflected event can not be reconstructed correctly, the signal to noise rationof gathers decreased. When the VSP-SSP correlation transformation applied to sparsegeometry, the reconstructed SSP gathers still has good quality. When integrate VSP andSSP data that recorded over the same area, one can transform VSP data into virtual SSPdata by VSP-SSP convolution transformation. Resulted virtual SSP data can be comparedwith actual SSP data, and used to determine the lithology of subsurface reflector body.During the transformation process, it is no need to know the location information ofreceivers in the well. We transform VSP data into virtual SSP data by VSP-SSPconvolution transformation, and then the conventional seismic processing used to virtualsurface gathers to obtain the image of subsurface structure. Compared with virtual SSPdata generated by correlation transformation, the signal to noise ratio of virtual SSP datagenerated by convolution transformation is higher, the resolution of deeper reflector isbetter, but the imaging range is limited. Test on the effect of sparse geometry on theVSP-SSP convolution transformation by numerical simulation. When the receivers in thewell distributed on local area, reconstructed SSP gathers by VSP-SSP convolutiontransformation with various degrees of effect, reflected event of shallow reflector can notbe reconstructed correctly, but the signal to noise ration of gathers is still high. When theVSP-SSP correlation transformation applied to sparse geometry, the reconstructed SSPgathers still has good quality. Combine VSP-SSP correlation transformation and VSP-SSPconvolution transformation, the advantages from two transformations combined together.Process the virtual SSP gathers by conventional image processing, the resulted imagingrange is same with it from correlation transformation and result in a huge increase in thesubsurface illumination, at the same time, the imaging accuracy is same with it fromconvolution transformation and deeper subsurface structure can be detected.
Passive source seismic data collected from the underground microseismic, generally,recorded ambient noise signal has relatively weak energy and low signal to noise ratio. It ishard to directly process the passive seismic data and use it to study of geological structurein the reservoir. Retrieved new seismic data set from ambient noise by seismicinterferometry, the noise into useful signal. The retrieved surface waves and reflected waves can be used to infer geological properties in the shallow subsurface structure anddeeper subsurface sreucture. This paper study application of seismic interferometrymethod for passive source seismic data by synthetic data and ambient seismic noise datasets. First of all, design forward modeling scheme suit for passive source seismic data,and then generate synthetic passive seismic data. Then we apply the seismic interferometrytechnique to approximately25hours of recordings of ambient seismic noise at the Ketzinexperimental CO2sequestration site, Germany. Common source gathers were generatedfrom the ambient noise for all available receivers along two seismic lines bycross-correlation of noise records. The retrieved response includes surface waves, refractedwaves and reflected waves. Comparison with surface waves and reflected waves fromactive source seismic survey, confirming the validity of the passive data processing.Stacking of common offset gathers increases the signal to noise ratio of retrieved commonsource gathers. A surface wave dispersion curve could be estimated from the retrievedsource gathers. Inversion of the dispersion curve results in an average1D and2D S-waveprofile below the site that is consistent with velocity models from traveltime tomographyof active source data from the area, but has higher resolution for the shallow subsurface. Inaddition, compared with velocity profile from other previous active source data, velocityprofile from surface wave dispersion curves also show similarities. Data processing of theretrieved common source gathers results in a stacked section that is in reasonableagreement with that of an active source section from same survey lines. The stackedsection from passive seismic survey have the potential to provide information on thedeeper subsurface structure， provide a effective technical support for the practicalapplication of passive source seismic exploration.
首先阐述了地震干涉法的基本原理，以互易定理为基础，按照无损任意非均匀介质的假设，推导准确的声学和弹性动力学格林函数的表述。而重构的格林函数中不仅包含两个检波点之间的直达波场，也包含所有一次反射和多次散射。对于地震干涉提取格林函数的表述是基于用于相关处理中的检波器是被一个封闭的震源面所包围的这样一个假设。一般情况下，往往不能满足地震干涉应用中所需的震源分布条件，因此文中通过数值模拟实验调查当理论条件不能被满足，且震源位于地下深部位置处时，震源分布特性对重构的反射响应的影响进行讨论。在地震勘探的实际应用中，通常不能满足理论上所需的连续震源分布，但在干涉计算过程中，震源间隔应尽量布设的密集一些，以防止空间假频的产生。当震源只集中分布于地下局部区域内时，重构的道集中只有部分反射同相轴能被恢复出来，且大大降低了重构响应的信噪比。一般而言，可以利用地震干涉法从地下的透射响应中重构出有效的反射信息，并对重构的反射信息进一步应用偏移处理，偏移成像可以改善重构共炮点道集中反射信息的不完整性，获得令人满意的结果。
文中对地震干涉法在主动源数据上的应用进行了深入系统的研究，主要针对的是VSP型数据。当常规地震勘探数据处理方法不能对上覆盖层复杂变化的地震资料获得满意的成像效果时，可以通过虚拟震源方法将VSP数据转换成SWP数据，变换之后，原来位于地表的震源被重建到井下位置处，震源-检波器的井下排列离地质目标体更近，且可以有效的避免近地表变化和复杂上覆盖对数据的影响，不需要知道上覆盖层的复杂速度模型，为处理这种复杂地震资料提供了一种有效途径。本文主要利用虚拟震源方法探测地下陡构造地质体，通过数值模拟发现地质体越陡峭，获得的成像效果越好，并且虚拟震源方法可以自动校正静校正量。虚拟震源方法对于弹性波数据的应用也是有效的，能够获得合理的成像结果。此外，模拟实验表明，虚拟震源方法应用于稀疏观测系统仍然是非常有效的，尽管随着震源及检波器间隔的增大，信噪比随之下降，但是仍能构建出正确的成像结果。
标准VSP数据的成像精度虽然很高，但是成像范围却很有限，仅能对从第一个检波器开始向下的一个三角锥形区域内的地质体进行成像。为了扩展地下地质构造的成像范围，可以利用地震干涉中的相关算法将VSP数据重构为虚拟SSP数据，并对虚拟SSP数据进行成像处理。由此得到的结果中，地下地质构造的成像范围等同于具有相同震源覆盖范围的常规地表地震勘探中获得的成像范围。此外，该变换校正了井中的静校正量，并且不需要知道井中震源激发时间，震源以及检波器的位置信息。文中利用VSP-SSP相关变换将正演模拟的VSP数据转换成虚拟SSP数据，并对重构的地表地震数据进行常规处理，获得地下地质构造的成像。从结果中看出，和标准VSP数据的成像范围相比，地下的成像范围被大大扩展了。通过数值模拟测试了观测系统对VSP-SSP相关变换方法所得结果的影响。当采用井下不同分布范围上局部检波器所获得的响应时，应用VSP-SSP相关变换方法重构的SSP道集受到不同程度的影响，部分反射同相轴不能被正确的重构出来，道集的信噪比也降低了。但是采用稀疏的观测系统时，VSP-SSP相关变换方法依然能够重构出质量较好的SSP道集。当需要对采集自同一地区的VSP和SSP数据进行整合时，可以利用VSP-SSP褶积变换将VSP数据转换成虚拟SSP数据。形成的虚拟SSP数据可以和实际SSP数据进行对比分析，用于确定地下反射体的岩性。而且在该变换过程中也不需要知道速度模型和VSP测量中检波器的位置信息。我们利用VSP-SSP褶积变换将正演模拟的VSP数据转换成虚拟SSP数据，并对重构的地表地震数据进行常规处理，获得地下地质构造的成像。与通过相关算法从VSP数据中重构出来的虚拟SSP数据相比，褶积算法重构的虚拟SSP道集的信噪比更高，对地下较深处反射体的分辨率更好，但是成像范围却很有限。同样通过数值模拟测试了观测系统对VSP-SSP褶积变换方法所得结果的影响。采用仅分布于井下局部范围上检波器所获得的响应时，应用VSP-SSP褶积变换方法重构的SSP道集受到不同程度的影响，位于浅层的反射同相轴不能被重构出来，但依然保持了较高的信噪比。对于稀疏的观测系统，VSP-SSP褶积变换方法依然能够重构出信噪比较高的虚拟SSP道集。将VSP-SSP相关变换和VSP-SSP褶积变换相联合，能够将两种变换方法的优点结合到一起，对构建出的虚拟SSP道集进行常规成像处理，所得结果的成像范围和相关算法中的一样，大大扩展了地下地质构造的成像范围，而成像精度和褶积算法中的一样，对来自于深部的地质体也能分辨出来，且具有较高的信噪比。
被动源地震数据主要采集自地下的微震，记录的背景噪声信号通常能量较弱，信噪比较低。直接对被动源数据进行成像处理，研究储层的地质构造有一定的技术难度。利用地震干涉法从背景噪声数据中提取出新的地震数据，能够把无序的噪声变成有用的信号，将重构的面波响应和反射波信息用于推断地下浅层和深部的地质构造特征。本文通过数值模拟和野外实际数据深入研究了地震干涉法在被动源数据上的应用。首先专门针对被动源地震数据的特点，设计了适合的正演模拟方案，并模拟生成了被动源地震数据。然后对采集自德国Ketzin地区的CO2冲注试验区记录时长将近25小时的背景噪声进行被动源地震干涉技术的应用研究。根据该工区内两条地震测线上所有可用检波器接收到的背景噪声记录，重构出共炮点道集。在重构的道集中可以观察到面波、折射波和反射波信息，与该地区主动源地震测量产生的共炮点道集相比较，验证了虚拟共炮点道集中重构的面波响应和反射波响应的有效性。通过产生共偏移距叠加道集提高虚拟共炮点道集的信噪比。利用重构的面波响应的分散特性估计出面波分散曲线，通过分散曲线的反演获得该地区地下的1D和2D剪切波速度剖面，获得的速度模型与根据该地区主动源数据的旅行时层析反演得到的速度模型有很好的一致性，但是从被动源数据中获得的速度模型对地表附近的浅层具有更高的分辨率。此外，与该地区其他主动源资料中获得的速度剖面相比，也有很好的匹配性。对重构的反射信息进行处理，产生的成像结果能够反映地下地质构造的特征，与对相同测线记录的主动源数据进行处理得到的成像结果有很好的一致性，为被动源地震勘探的实际应用提供了有力的技术支持。
Seismic interferometry method and its application in the field data, greatly enrichedthe understanding and study of seismic wave propagation, it is considered to be a majoradvance in the development of geophysics. The scope of application of seismicinterferometry method is from the one-dimensinal layered acoustic medium to thethree-dimensional elastic medium, arbitrary inhomogeneous media or anisotropic media,and even can be extended to the moving media and attenuation media. The basic idea ofseismic interferometry method is retrieval of useful signal from the chaotic disorderseismic signal, new seismic signal is generated by the interference of recorded seismicsignal, the new seismic signal not only include the properties of original seismic signal, butalso retrieve some important features that the original signal does not have. Seismicinterferometry method can be used in both passive seismic measurements and controlledsource seismic measurements, the source can be certain source or noise source. Activesource seismic interferometry provide a new approach for seismic exploration, it can beused for a variety of data types. Take application of vertical seismic profile (VSP) data asrepresentative, VSP data can be transferred into the other data type by seismicinterferometry method, such as single seismic profile (SWP) data, surface seismic profile(SSP). Passive source seismic interferometry method applied to retrieval of useful signalfrom ambient seismic noise to obtain geological properties of medium, and then infer thecharacteristics of subsurface structure. In this paper, we respectively study and discussionthe application of seismic interferometry method for active source seismic data and passivesource seismic data.
First of all, describe the basic principle of seismic interferometry methods, based onreciprocity theorems, underlying assumption of inhomogeneous, lossless medium, acousticand elastic Green’s function representations can be derived. The reconstructed Green’sfunction not only include the direct wave between two receivers, but also include primaryreflection and multiple scattering. Retrieved Green’s function by seismic interferometry isbased on the assumption that receivers are surrounded by a closed source surface. Undernormao circumstances, the source distribution conditions required by seismicinterferometry can not be met. So in this paper, we investigate that when the circumstancethat theory condition can not be met, the sources distributed in the subsurface, how the properties of source distribution effect on the retrieval of reflected waves. In the practicalapplication of seismic exploration, it always hard to meet the theoretically required for acontinuois source distribution, but during the calculation, in order to avoid aliasing, thesource interval should be set intensively. When the source only distributed in the localsubsurface area, only part of reflected event can be reconstructed in the retrieved gathers,the signal to noise ration is decreased seriously. In general, retrieved reflection informationfrom subsurface transmission response by seismic interferometry method is very effective.Process the retrieved reflection information and then used to migration process, theresulted image can be improve the incomplete reflected waves in the retrieved commonsource gathers and obtain satisfactory results.
In the paper, we study on application of seismic interferometry method for activesource seismic data, here mainly target at VSP data. When the conventional seismic dataprocessing methods failed to acquire the satisfactory imaging results of seismic datainclude overburden complexity, we could apply the virtual source method to convert thetraditional vertical seismic profile gathers into virtual single well profile gathers. After that,redatumed geometry of source-receiver arrays is closer to the structure target and belowthe complex overburden so that we can get the better image resolution of the target afterconventional process and avoid the effected by unknown velocities of complex overburden,provide an efficient way to handle this complex seismic data. In the paper, the virtualsource method applied to detect steep structure such as faults present in the subsurface.Though the synthetic data found that the steeper structure the better image results, andautomatically correct the statics. The virtual source method is still work well for elasticdata sets. Moreover, modeling experiment show that virtual source method is still effectivefor the sparse geometry, we could get the good imaging results. Although the signal tonoise ratio decreased as source interval and receiver interval increased, the image result iscorrect.
Although the imaging accuracy of the standard VSP data is high, the imaging range isvery limited that is defined by a small triangle, with the triangle apex at the shallowestreceiver. In order expand the imaging range of subsurface structure, one can transformVSP data into virtual SSP data by correlation algorithm of seismic interferometry, and thenthe imaging process used to virtual SSP data. The imaging range of subsurface structurefrom the results has same area as that of the original VSP source distribution along thesurface. In addition, this transform eliminates well statics and the need to know the sourceor receiver position in the well and excitation time of source in the well. In this paper,VSP-SSP correlation transformation effectively transform synthetic VSP data into virtual SSP data, and then migrated by a standard surface seismic data processing to give ainterferometric image of subsurface structure. From the results, compared with standardVSP image, VSP-SSP correlation transformation is a huge increase in the subsurfaceillumination. Test on the effect of geometry on the VSP-SSP correlation transformation bynumerical simulation. When the receivers in the well distributed on local area,reconstructed SSP gathers by VSP-SSP correlation transformation with various degrees ofeffect, part of reflected event can not be reconstructed correctly, the signal to noise rationof gathers decreased. When the VSP-SSP correlation transformation applied to sparsegeometry, the reconstructed SSP gathers still has good quality. When integrate VSP andSSP data that recorded over the same area, one can transform VSP data into virtual SSPdata by VSP-SSP convolution transformation. Resulted virtual SSP data can be comparedwith actual SSP data, and used to determine the lithology of subsurface reflector body.During the transformation process, it is no need to know the location information ofreceivers in the well. We transform VSP data into virtual SSP data by VSP-SSPconvolution transformation, and then the conventional seismic processing used to virtualsurface gathers to obtain the image of subsurface structure. Compared with virtual SSPdata generated by correlation transformation, the signal to noise ratio of virtual SSP datagenerated by convolution transformation is higher, the resolution of deeper reflector isbetter, but the imaging range is limited. Test on the effect of sparse geometry on theVSP-SSP convolution transformation by numerical simulation. When the receivers in thewell distributed on local area, reconstructed SSP gathers by VSP-SSP convolutiontransformation with various degrees of effect, reflected event of shallow reflector can notbe reconstructed correctly, but the signal to noise ration of gathers is still high. When theVSP-SSP correlation transformation applied to sparse geometry, the reconstructed SSPgathers still has good quality. Combine VSP-SSP correlation transformation and VSP-SSPconvolution transformation, the advantages from two transformations combined together.Process the virtual SSP gathers by conventional image processing, the resulted imagingrange is same with it from correlation transformation and result in a huge increase in thesubsurface illumination, at the same time, the imaging accuracy is same with it fromconvolution transformation and deeper subsurface structure can be detected.
Passive source seismic data collected from the underground microseismic, generally,recorded ambient noise signal has relatively weak energy and low signal to noise ratio. It ishard to directly process the passive seismic data and use it to study of geological structurein the reservoir. Retrieved new seismic data set from ambient noise by seismicinterferometry, the noise into useful signal. The retrieved surface waves and reflected waves can be used to infer geological properties in the shallow subsurface structure anddeeper subsurface sreucture. This paper study application of seismic interferometrymethod for passive source seismic data by synthetic data and ambient seismic noise datasets. First of all, design forward modeling scheme suit for passive source seismic data,and then generate synthetic passive seismic data. Then we apply the seismic interferometrytechnique to approximately25hours of recordings of ambient seismic noise at the Ketzinexperimental CO2sequestration site, Germany. Common source gathers were generatedfrom the ambient noise for all available receivers along two seismic lines bycross-correlation of noise records. The retrieved response includes surface waves, refractedwaves and reflected waves. Comparison with surface waves and reflected waves fromactive source seismic survey, confirming the validity of the passive data processing.Stacking of common offset gathers increases the signal to noise ratio of retrieved commonsource gathers. A surface wave dispersion curve could be estimated from the retrievedsource gathers. Inversion of the dispersion curve results in an average1D and2D S-waveprofile below the site that is consistent with velocity models from traveltime tomographyof active source data from the area, but has higher resolution for the shallow subsurface. Inaddition, compared with velocity profile from other previous active source data, velocityprofile from surface wave dispersion curves also show similarities. Data processing of theretrieved common source gathers results in a stacked section that is in reasonableagreement with that of an active source section from same survey lines. The stackedsection from passive seismic survey have the potential to provide information on thedeeper subsurface structure， provide a effective technical support for the practicalapplication of passive source seismic exploration.
引文
[1] Wapenaar K, Snieder R. From order to disorder to order: A philosophical view on seismicinterferometry,2007[C].77th Annual International Metting, SEG, Expanded Abstracts,2683–2687.
[2]王伟涛,倪四道,王宝善.地球背景噪声干涉应用研究的新进展[J].中国地震,2011,27(1):1-13.
[3] Wang B, Zhu P, Chen Y, et al. Continuous subsurface velocity measurement with coda waveinterferometry [J]. J. Geophys. Res.,2006,113: B12313, doi:10.1029/2007JB005023.
[4] Yao H, R D van der Hilst, de Hoop M V. Surface-wave array tomography in SE Tibet from ambientseismic noise and two-station analysis: I-phase velocity maps [J]. Geophys. J. Int.,2006,166:732–744.
[5] Ritzwoller M, Shapiro N, Pasyanos M, Bensen G, Tang Y. Short-period surface wave dispersionmeasurements from ambient seismic noise in North Africa, the Middle East, and Central Asia,161–170[R].27th Seismic Research Review: Ground-Based Nuclear Explosion MonitoringTechnologies,2005.
[6] Florent B, Nikolai M S, Michel C, Slexandre N, Valérie F.3-D surfacewave tomography of thePiton de la Fournaise volcano using seismic noise correlations [J]. Geophysical Research Letters,2007,34: L02305, doi:10.1029/2006GL028585.
[7] Claerbout J F. Synthesis of a layered medium from its acoustic transmission response [J].Geophysics,1968,33:264–269.
[8] Rickett J, Claerbout J F. Acoustic daylight imaging via spectral factorization: helioseismology andreservoir monitoring [J]. The Leading Edge,1999,18:957-960.
[9] Wapenaar C P A, Draganov D, Thorbecke J, Fokkema J T. Theory of acoustic daylight imagingrevisited,2002[C].72nd SEG annual meeting, Expanded Abstracts,2269-2272.
[10]Scherbaum F. Seismic imaging of the site response using microearthquake recordings, Part I:Method [J]. Bull. Seis. Soc. Am.,1987a,77:1905–1923.
[11]Scherbaum F. Seismic imaging of the site response using microearthquake recordings, Part II:Application to the Swabian Jura, southwest Germany, seismic network[J]. Bull. Seis. Soc. Am.,1987b,77:1924–1944.
[12]Cole S. Passive seismic and drill-bit experiments using2-D arrays [D]. Ph.D. dissertation, StanfordUniversity,1995.
[13]Daneshvar M, Clay C, Savage M. Passive seismic imaging using microearthquakes [J].Geophysics,1995,60:1178–86.
[14]Baskir C, Wells C. Sourceless reflection seismic exploration [J]. Geophysics,1975,40,158.
[15]Duvall T L, S M Jefferies, J W Harvey, M A Pomerantz. Time-distance helioseismology [J].Nature,1993,362:430–432.
[16]Gilles P M, Duvall T L, Scherrer P H, Bogart R S., A subsurface flow of material from the Sunequator’s to its poles [J]. Nature,1997,390:52–54.
[17]Cassereau D, Fink M. Focusing with plane time reversal mirrors: an efficient alternative to closedcavities [J]. Journal of the Acoustical Society of America,1993,94(4):2373-2386.
[18]Fink M. Time reversal mirrors [J]. J. Phys. Rev.: Appl. Phys.,1993,26,1330-1350.
[19]Fink M. Time reversed acoustics [J]. Physics Today,1997,50:34-40.
[20]Roux P, Fink M. Green′s function estimation using secondary sources in a shallow waterenvironment [J]. Journal of the Acoustical Society of America,2003,113:1406-1416.
[21]Snieder R, Grêt A, Douma H, Scales J. Coda wave interferometry for estimating nonlinearbehavior in seismic velocity [J]. Science,2002,295:2253-2255.
[22]Snieder R. Coda wave interferometry and the equilibration of energy in elastic media [J]. Phys.Rev. E,2002,66:046615-1-046615-8.
[23]Snieder R, Hagerty M. Monitoring change in Volcanic Interiors using Coda WaveInterferometry:Application to Arenal Volcano, Costa Rica [J]. Geophys. Res. Lett.,2004,31(9): L09608.1-L09608.5.
[24]Grêt A R, Snieder R, Aster C, Kyle P R. Monitoring rapid temporal changes in a volcano with codawave interferometry [J]. Geophys. Res. Lett.,2005,32(6): L06304.1-L06304.
[25]Snieder R. The theory of coda wave interferometry [J]. Pureand Appl. Geophys,2006,163:455-473.
[26]Fink M.1992, Time-reversal of ultrasonic fields: Basic principles [J]. IEEE Transactions onUltrasonics, Ferroelectrics, and Frequency Control,39,555–566.
[27]Fink M. Time-reversal acoustics in complex environments [J]. Geophysics,2006,71(4):SI151–SI164.
[28]Lobkis O I, R L Weaver. On the emergence of the Green’s function in the correlations of a diffusefield [J]. Journal of the Acoustical Society of America,2001,110:3011–3017.
[29]Lerosey G, de Rosny J, Tourin A, Fink M. Focusing beyond the diffraction limit with far-fieldtime reversal [J]. Science,2007,315,1120–1122.
[30]Derode A, Larose E, Campillo M, Fink M. How to estimate the Green’s function of aheterogeneous medium between two passive sensors-Application to acoustic waves [J]. AppliedPhysics Letters,2003a,83:3054–3056.
[31]Derode A, Larose E, Tanter M, de Rosny J, Tourin A, Campillo M, Fink M. Recovering theGreen’s function from field-field correlations in an open scattering medium (L)[J]. Journal of theAcoustical Society of America,2003b,113:2973–2976.
[32]Blomgren P, Papanicolaou G, Zhao H. Super-resolution in time-reversal Acoustics [J]. J. Acoust.Soc. Am.,2002,111(1),230–48.
[33]Borcea L, Papanicolaou G, Tsogka C, Berryman J. Imaging and time reversal in random media [J].Inverse Problems,2002,18,1247–79.
[34]Borcea L, Papanicolaou G, Tsogka C. Coherent interferometric imaging in clutter [J], Geophysics,2006,71, SI165–SI175.
[35]Schuster G T. Seismic interferometry [M]. New York: Cambridge University Press,2009.
[36]Weaver R L, Lobkis O I. Ultrasonics without a source: thermal fluctuation correlations at MHzfrequencies [J]. Physical Review Letters,2001,87:134301-1–134301-4.
[37]Weaver R L, Lobkis O I.2002, On the emergence of the Green’s function in the correlations of adiffuse field: Pulse-echo using thermal phonons [J]. Ultrasonics,40,435–439.
[38]Campillo M, Paul A. Long-range correlations in the diffuse seismic coda [J]. Science,2003,299:547–549.
[39]Weaver R L, Lobkis O I. Elastic wave thermal fluctuations, ultrasonic waveforms by correlation ofthermal phonons [J]. J. Acoust. Soc. Am.,2003,113:2611–2621.
[40]Weaver R L, Lobkis O I. Diffuse fields in open systems and the emergence of the Green’s function[J]. J. Acoust. Soc. Am.,2004,116(5):2731–2734.
[41]Larose E, Derode A, Campillo M, Fink M. Imaging from one-bit correlations of wide-band diffusewavefield [J]. J. Appl. Phys.,2004,95:8393–8399.
[42]Larose E, Lobkis O I, Weaver R L. Passive correlation imaging of a buried scatterer[J]. J. Acoust.Soc. Am.,2006,119:3549–3552.
[43]Malcolm A, Scales J, Tiggelen B. Extracting the Green function from diffuse, equipartitionedwaves [J]. Physical Review E,2004,70:015601(R)-1–015601(R)-4.
[44]Roux P, Kuperman W A, the NPAL Group. Extracting coherent wavefronts from acoustic ambientnoise in the ocean [J]. J. Acoust. Soc. Am.,2004,116:1995–2003.
[45]Sabra K G, Roux P, Kuperman W A. Arrival-time structure of the time-averaged ambient noisecross-correlation function in an oceanic waveguide [J]. J. Acoust. Soc. Am.,2005,117:164–174.
[46]Sabra K G, P Roux, W A Kuperman.2005, Arrival-time structure of the time-averaged ambientnoise cross-correlation function in an oceanic waveguide [J]. Journal of the Acoustical Society ofAmerica,117,1:164–174.
[47]Shapiro N, Campillo M. Emergence of broadband Rayleigh waves from correlations of theambient seismic noise [J]. Geophys. Res. Lett.,2004,31: L07614.
[48]Shapiro N, Campillo M, Stehly L, Ritzwoller M. High-resolution surface-wave tomography fromambient seismic noise [J]. Science,2005,307:1615–18.
[49]Sabra K, Gerstoft P, Roux P, Kuperman W, Fehler M. Extracting time-domain Green’s functionsfrom ambient seismic noise [J]. Geophys. Res. Lett.,2005b,32: L03310.
[50]Gerstoft P, Sabra K, Roux P, Kuperman W A, Fehler M. Green’s functions extraction andsurface-wave tomography from microseisms in southern California [J]. Geophysics,2006,71:SI23–SI32.
[51]Larose E, Margerin L, Derode A, van Tiggelen B, Campillo M, Shapiro N, Paul A, Stehly L, TanterM. Correlation of random fields: An interdisciplinary review [J]. Geophysics,2006,71: SI11–SI21.
[52]Gudmundsson O, Khan A, Voss P. Rayleigh-wave group-velocity of the Icelandic crust fromcorrelation of ambient seismic noise [J]. Geophys. Res. Lett.,2007,34: L14314.
[53]Gouedard P, Stehly L, Brenguier F, Campillo M, de Verdiere Y, Larose E, Margerin L, Roux P,Sanchez-Sesma F, Shapiro N, Weaver R. Cross-correlation of random fields: Mathematicalapproach and applications [J]. Geophysical Prosp.2008,56:375–393.
[54]Nunziata C, De Nisco G, Panza G F. S-wave profiles from noise cross correlation at small scale [J].Engineering Geology,2009,105:3–4.
[55]Shragge J, Artman B, Wilson C. Teleseismic shot-profile migration [J]. Geophysics,2006,71:SI221–SI229.
[56]Schuster G T, Rickett J. Daylight imaging in V(x,y,z) media [R]. Utah Tomography andModeling-Migration Project Midyear Report and SEP Report at Stanford University,2000.
[57]Yu J, Schuster G T, Crosscorrelogram migration of IVSPWD data [J].71st Ann. Internat. Mtg.,SEG Expanded Abstracts,2001,456–459.
[58]Yu J, Schuster G T, Enhancing illumination coverage of VSP data by cross-correlogram migration,74th Ann. Internat. Mtg., SEG Expanded Abstracts,2004,2501–2504.
[59]Yu J, Schuster G T, Crosscorrelogram migration of inverse vertical seismic profile data [J].Geophysics,2006,71, S1–S11.
[60]Hornby B E, Yu J. Interferometric imaging of a salt flank using walkaway VSP data [J]. TheLeading Edge,2007,26:760–763.
[61]Sheng J, Schuster G T, Nowack R. Imaging of crustal layers by teleseismic ghosts [C]. AmericanGeophysical Union, Fall Mtg.,2002, S32C–0658.
[62]Sheng J, Schuster G T, Pankow K, Pechmann J, Nowack R. Coherence-weighted wavepathmigration of teleseismic data [C]. American Geophysical Union, Fall Mtg.,2003, S11E–0344.
[63]Sheng J. Migrating multiples and primaries in CDP data by crosscorrelation Migration [J].71stAnn. Internat. Mtg., SEG Expanded Abstracts,2001,1297–1300.
[64]He R, Hornby B, Schuster G.3D wave-equation interferometric migration of VSP free-surfacemultiples [J]. Geophysics,2007,72: S195–S203.
[65]Schuster G T, Yu J, Sheng J, Rickett J. Interferometric/daylight seismic imaging [J]. GeophysicalJournal International,2004,157:838–852.
[66]Wapenaar K, Draganov D, Thorbecke J, Fokkema J. Theory of acoustic daylight imaging revisited[C].72nd Annual International Meeting, SEG, Expanded Abstracts,2002,2269–2272.
[67]Wapenaar K. Synthesis of an inhomogeneous medium from its acoustic transmission response [J].Geophysics,2003,68:1756-1759.
[68]Wapenaar K, Thorbecke J, Draganov D. Relations between reflection and transmission responsesof three-dimensional inhomogeneous media [J]. Geophysical Journal International,2004,156:179–194.
[69]Wapenaar K. Retrieving the elastodynamic Green’s function of an arbitrary inhomogeneousmedium by crosscorrelations [J]. Physical Review Letters,2004,93:254301-1–254301-4.
[70]Wapenaar K, Fokkema J, Snieder R. Retrieving the Green’s function in an open system bycross-correlation: A comparison of approaches(L)[J]. Journal of the Acoustical Society of America,2005,118:2783–2786.
[71]Wapenaar K. Green’s function retrieval by cross-correlation in case of one-sided illumination [J].Geophysical Research Letters,2006,33: L19304-1-L19304-6.
[72]Willis M E, R Lu, X Campman, M N Toks z, Y Zhang, M V de Hoop. A novel application oftime-reversed acoustics: Salt-dome flank imaging using walkaway VSP surveys [J]. Geophysics,2006,71(2): A7–A11.
[73]Artman B. Imaging passive seismic data [J]. Geophysics,2006,71:177-187.
[74]Calvert R W, Bakulin A, Jones T C. Virtual sources, a new way to remove overburden problems[J].66th Mtg., European Association of Geoscientists and Engineers, Extended Abstracts,2004,P234.
[75]Bakulin A, Calvert R. Virtual source: New method for imaging and4D below complex overburden[C].74th Ann. Internat. Mtg., SEG Expanded Abstracts,2004,2477–2480.
[76]Bakulin A, Calvert R. The virtual source method: Theory and case study [J]. Geophysics,2006,71:SI139–SI150.
[77]Berkhout A, Verschuur D. Transformation of multiples into primary reflections,73rd Ann. Internat.Mtg., SEG Expanded Abstracts,2003,1925–1928.
[78]Berkhout A, Verschuur D. Imaging of multiple reflections, Geophysics,2006,71, SI209–SI220.
[79]Jiang Z, Yu J, Hornby B, Schuster G T. Migration of multiples [J]. The Leading Edge,2005,315–318.
[80]Jiang Z, Sheng J, Yu J, Schuster G T, Hornby B. Migration methods for imaging different-ordermultiples [J]. Geoph. Prosp.,2007,55,1–19.
[81]Guitton A. Shot-profile migration of multiple reflections [J].72nd Ann. Internat. Mtg., SEGExpanded Abstracts,2002,1296–1299.
[82]Shan G, Guitton A. Migration of surface-related multiples: Tests on Sigsbee2B dataset [J].74thAnn. Internat. Mtg., SEG Expanded Abstracts,2004,1285–1288.
[83]Meng F, Schuster G T. Autocorrelogram imaging of Friendswood RVSP data [R] Annual UtahReport,320-338.
[84]Xiao X, Schuster G T. Redatuming CDP data below salt with VSP Green’s function [J].76th Ann.Internat. Mtg., SEG Expanded Abstracts,2006,3511–3515.
[85]Xiao X, Zhou M, Schuster G T. Salt-flank delineation by interferometric imaging of transmitted Pto S waves, Geophysics,2006,71, SI197–SI207.
[86]Draganov D, Wapenaar K, Thorbecke J. Passive seismic imaging in the presence of white noisesources [J]. The Leading Edge,2004,23,889–892.
[87]Draganov D, Wapenaar K, Mulder W, Singer J. Seismic interferometry on background-noise fielddata [C].76th Annual International Meeting, SEG, Expanded Abstracts,2006,590–593.
[88]Draganov D, Wapenaar K, Thorbecke J. Seismic interferometry: Reconstructing the earth’sreflection response [J]. Geophysics,2006,71: SI61–SI70.
[89]Draganov D, Wapenaar K, Mulder W, Singer J, Verdel A. Retrieval of reflections from seismicbackground-noise measurements [J]. Geophysical Research Letters,2007,34, L04305-1–L04305-4.
[90]Draganov D, Campman X, Thorbecke J, Verdel A, Wapenaar K. Reflection images from seismicnoise [J]. Geophysics,2009,74: A63–A67.
[91]Artman B. Imaging passive seismic data [J]. Geophysics,2006,71, SI177–SI187.
[92]Wapenaar K, Draganov D, Snieder R, Campman X, Verdel A. Tutorial on seismic interferometry:Part1-Basic principles and applications [J]. Geophysics,2010,75(5):75A195-75A209.
[93]Dorman J, Ewing M. Numerical inversion of surface-wave dispersion data and crust-mantlestructure in the New York-Pennsylvania area [J]. J. geophys. Res.,1962,67:5227–5241.
[94]Lu R, Venkataraman A, Payne M, Zhang J. Application of noise interferometry to obtaintime-lapse velocity variations during a steam stimulation cycle at Cold Lake [C].79th AnnualInternational Metting, SEG, Expanded Abstracts,2009,1673–1677.
[95]Snieder R, Safak E. Extracting the building response using seismic interferometry: Theory andapplication to the Millikan Library in Pasadena, California [J]. Bulletin of the SeismologicalSociety of America,2006,96:586–598.
[96]Curtis A P, Gerstoft H, Sato R, Snieder R, Wapenaar K. Seismic interferometry—Turning noiseinto signal [J]. The Leading Edge,2006,25:1082–1092.
[97]Wapenaar K, van der Neut J, Ruigrok E. Passive seismic interferometry by multidimensionaldeconvolution [J]. Geophysics,2008,73(6):51-56.
[98]Vasconcelos I, Snieder R. Interferometry by deconvolution, Part1—Theory for acoustic waves andnumerical examples [J]. Geophysics,2008,73(1): S115-S128.
[99]Vasconcelos I, Snieder R. Interferometry by deconvolution, Part2—Theory for elastic waves andapplication to drill-bit seismic imaging [J]. Geophysics,2008,73(3):129-141.
[100]Schuster G T. Seismic interferometric imaging with waveforms [R]. Utah Tomography andModeling-Migration Project Mideyear Report,1999,121-130
[101]Schuster G T, Followill F, Lewis K J, Yu J, Liu Z. Autocorrelogram migration: Theory [J].Geophysics,2003,68(5):1685-1694.
[102]Mercier J, Bostock M, Baig A. Improved Green’s functions for passive-source structural studies[J]. Geophysics,2006,71: SI95-SI102
[103]Wang B, Zhu P, Chen Y, et al. Continuous subsurface velocity measurement with coda waveinterferometry [J]. J. Geophys. Res.,2006,113: B12313, doi:10.1029/2007JB005023
[104]房立华,吴建平,吕作勇.华北地区基于噪声的瑞利面波群速度层析成像[J].地球物理学报,2009,52(3):663-671.
[105]齐诚,陈棋福,陈颙.利用背景噪声进行地震成像的新方法[J].地球物理学进展,2007,22(3):771-777.
[106]汪承灏,魏炜.改进的时间反转法用于有界面时超声目标探测的鉴别[J].声学学报,2002,27(3):193-197.
[107]吴昊,张碧星,汪承灏.声波在固体板中的多径传播及其时间反转声场[J].声学学报,2005,30(3):215-221.
[108]陆铭慧,张碧星,汪承灏.时间反转法在水下通信中的应用[J].声学学报,2005,30(4):215-221.
[109]Wang B, Zhu G, Gao J. Joint interferometry imaging of walkaway VSP data [J]. AppliedGeophysics,2010,7(1):41-48.
[110]Hill S J. Geophysics bright spots [J]. The Leading Edge,2006,25:1082-1092.
[111]单新建,马瑾,王长林,柳稼航.利用差分干涉雷达测量技术(D-InSAR)提取同震形变场[J].地震学报,2007,24(4),413-420.
[112]朱新慧,孙付平.用甚长基线干涉测量数据检测冰期后地壳回弹[J].地球物理学报,2005,48(2):308-313.
[113]Wapenaar K, Fokkema J. Green’s function representations for seismic interferometry [J].Geophysics,2006,71(4): SI33–SI46.
[114]de Hoop A T. An elastodynamic reciprocity theorem for linear, viscoelastic media [J]. AppliedScientific Research,1966,16:39–45.
[115]de Hoop A T. Time-domain reciprocity theorems for acoustic wave fields in fluids with relaxation[J]. Journal of the Acoustic Society of America,1988,84,1877–1882.
[116]Aki K, Richards P G. Quantitative Seismology: Theory and Methods [M]. San Francisco:Freenman and Coopration,1980,146-177.
[117]何樵登.地震波理论[M].北京:地质出版社,1988.
[118]de Hoop A T. Handbook of radiation and scattering of waves [M]. Academic Press, Amsterdam,1995.
[119]Rayleigh J W S. The theory of sound, Volume II [M]. Macmillan and Co.1878.
[120]Fokkema J T, van den Berg P M. Seismic applications of acoustic reciprocity [M]. ElsevierScience,1993.
[121]Bojarski N N. Generalized reaction principles and reciprocity theorems for the wave equations,and the relationship between the time advanced and the time-retarded fields [J]. Journal of theAcoustic Society of America,1983,74,281–285.
[122]Bleistein N. Mathematical methods for wave phenomena: Academic Press Inc,1984.
[123]Knopoff L, Gangi A F. Seismic reciprocity [J]. Geophysics,1959,24,681–691.
[124]Aki K, Richards P G, Quantitative seismology, vol.1[M]. W. H. Freeman and Company,1980.
[125]Mehta K, Snieder R. Acquisition geometry requirements for generating virtual-source data [J].The Leading Edge,2008,27(5):620-629.
[126]Hardage B. Vertical seismic profiling: Principles [M]. Elsevier Science Publications Co., Inc.,2000.
[127]Blakeslee S, Chen S, Krebs J, Srnka L. TVSC: Twin VSP simulation of cross-well data, Astrategy for low-cost monitoring of EOR processes [J].63rd Ann. Internat. Mtg., SEG ExpandedAbstracts,1993,12:9–12.
[128]Cao W, Schuster G T. Natural redatuming of VSP data to virtual CDP data [R].2005UTAM Ann.Report, University of Utah,2005,231–248.
[129]Krebs J, Fara D, Berlin A. Accurate migration using offset-checkshot surveys [J]65th Ann.Internat. Mtg., SEG Expanded Abstracts,1995,14:1186–1189.
[130]Schuster G T. Reverse time migration=generalized diffraction stack migration [J].72nd Ann.Internat. Mtg., SEG Expanded Abstracts,2002,21:1280–1283.
[131]Brandsberg-Dahl S, Hornby B, Xiao X. Migration of surface seismic data with VSP Green’sfunctions [J]. The Leading Edge,2007,26:778–781.
[132]Xiao X. Local reverse-time migration with VSP Green’s functions [D]. Salt Lake City, Utah:University of Utah,2008.
[133]Minato S, Onishi K, Matsuoka T, Okajima Y, Tsuchiyama J, Nobuoka D, Azuma H, Iwamoto T.Cross-well seismic survey without borehole source [J].77th Ann. Internat. Mtg., SEG ExpandedAbstracts,2007,26:1357–1361.
[134]Hornby B, Yu J. Single-well imaging of a salt flank using walkaway VSP data [J].76th Ann.Internat. Mtg., SEG Expanded Abstracts,2006,25:3492–3496.
[135]Hornby B, Yu J, Sharp J, Ray A, Quist Y, Regone C. VSP: Beyond time-to-depth [J]. The LeadingEdge,2006,25:446–452.
[136]He R, Hornby B, Schuster G T.3D wave-equation interferometric migration of VSP multiples [J].76th Ann. Internat. Mtg., SEG Expanded Abstracts,2006,25:3442–3445.
[137]Willis M, Lu R, Campman X, Toksoz N, Zhang Y, de Hoop M. A novel application oftime-reversed acoustics: Salt-dome flank imaging using walkaway VSP surveys [J]. Geophysics,2006,71: O43–O51.
[138]Korneev V, Bakulin A. On the fundamentals of the virtual source method [J]. Geophysics,2006,71: A13-A17.
[139]Chakroun N, Fink M, Wu F. Time reversal processing in ultrasonic nondestructive testing [J].IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control,1995,42:1087–1098.
[140]Draeger C, Cassereau D, Fink M. Acoustic time reversal with mode conversion at a solid-fluidinterface [J]. Applied Physics Letters,1998,72:1567–1569.
[141]Fink M, Cassereau D, Derode A, Prada C, Roux P, Tanter M, Thomas J L, Wu F. Time-reversedacoustics [J]. Reports on Progress in Physics,2000,63:1933–1995.
[142]Fink M, Prada C. Acoustic time-reversal mirrors [J]. Inverse Problems,2001,17: R1–R38.
[143]Tsogka C, Papanicolaou G. reversal through a solid-liquid interface and super-resolution [J].Inverse Problems,2002,18:1639–1657.
[144]Morse P, Feshbach H. Methods of Theoretical Physics (Part I)[M]. New York: McGraw-Hill,1953.
[145]F rster A, Norden B, Zinck-J rgensen K, Frykman P, Kulenkampff J, Spangenberg E, Erzinger J,Zimmer M, Kopp J, Borm G, Juhlin C, Cosma C, Hurter S. Baseline characterization of theCO2SINK geological storage site at Ketzin, Germany [J]. Environmental Geosciences,2006,13:145–160.
[146]Martens S, Liebscher A, M ller F, Würdemann H, Schilling F, Kühn, M, Ketzin Group. Progressreport on the first European on-shore CO2storage site at Ketzin (Germany)-Second year ofinjection [R]. Energy Procedia,2011,4:3246–3253.
[147]Juhlin C, Giese R, Zinck-J rgensen K, Cosma C, Kazemeini H, Juhojuntti N, Lüth S, Norden B,F rster A.3D baseline seismics at Ketzin, Germany: The CO2SINK project [J]. Geophysics,2007,72: B121–B132.
[148]Derode A, Tourin A, Fink M. Ultrasonic pulse compression with one-bit time reversal throughmultiple scattering [J]. J. Appl. Phys.,1999,85:6343–6352.
[149]Dong S, Sheng J, Schuster J T. Theory and practice of refraction interferometry [J].76th AnnualInternational Metting, SEG, Expanded Abstracts,2006,26:3021–3025.
[150]Mikesell D, van Wijk K, Calvert A, Haney M. The virtual refraction: Useful spurious energy inseismic interferometry [J]. Geophysics,2009,74: A13–A17.
[151]Nichols J, Mikesell D, van Wijk K. Application of the virtual refraction to near-surfacecharacterization at the Boise Hydrogeophysical Research Site [J]. Geophysical Prospecting,2010,58:1011–1021.
[152]Rayleigh L. On waves propagated along the plane surface of an elastic solid [J]. Proc. Lond.Math. Soc.,1885,17:4-11.
[153]Haskell N A. The dispersion of surface waves on multilayered media [J]. Bulletin of theSeismological Society of America,1953,43,17-34.
[154]Babuska V, Cara M. Seismic anisotropy in the earth [M]. Academic Publishers,1991.
[155]Sheriff R E, Geldart L P. Exploration seismology (volume II): History, theory, and dataacquisition [M]. New York: Cambridge University Press,1983.
[156]Stokoe K H, Wright S G, Bay J A, Ro sset J M. Characterization of geotechnical sites by SASWmethod [M]. In: R.D. Woods, Editor, Geophysical Characterization of Sites: Oxford Publishers,1994.
[157]Aki K. Space and time spectra of stationary stochastic waves, with special reference tomicrotremors [J]. Bull. Earthq. Res. Inst.,1957,35:415–456.
[158]Matthews M C, Hope V S, Clayton C R I. The use of surface waves in the determination ofground stiffness profiles [J]. Geotechnical Eng., Proc.Inst.Civ.Eng,1996,119:84-95.
[159]Aki K, Richards P G. Quantitative Seismology [M]. San Francisco: W. H. Freeman and Company,1980.
[160]Dorman J, Ewing M. Numerical inversion of surface-wave dispersion data and crust-mantlestructure in the New York-Pennsylvania area [J]. Journal of Geophysical Research,1962,67:5227-5241.
[161]Song Y Y, Castagna J P, Black R A, Knapp R W. Sensitivity of near-surface shear-wave velocitydetermination from Rayleigh and Love waves [J]. Technical Program with Biographies, SEG,59th Annual Meeting, Dallas, TX,1989,8:509-512.
[162]Turner G. Aliasing in the tau-p transform and the removal of spatially aliased coherent noise [J].Geophysics,1990,55:1496-1503.
[163]Rix G J, Leipski E A. Accuracy and resolution of surface wave inversion [J]. In: Bhatia S K,Blaney G W,(Eds), Recent Advances in Instrumentation, Data Acquisition and Testing in SoilDynamics. Am. Soc. Civil Eng.,1991, pp.17-32.
[164]Nazarian S, Stokoe II K H, Hudson W R. Use of spectral analysis of surface waves method fordetermination of moduli and thicknesses of pavement systems [J]. Transport Research Record,1983,930:38-45.
[165]Sanchez-Salinero I, Rosset J M, Shao K Y, Stokoe II K H, Rix G J. Analytical evaluation ofvariables affecting surface wave testing of pavements [J]. Transport Res Rec,1987,1136:86-95.
[166]Sheu J C, Stokoe II K H, Roesset J M. Effect of reflected waves in SASW testing of pavements[J]. Transportation Res. Record,1988,1196:51-61.
[167]Stokoe II K H, Rix G J, Nazarian S. In situ seismic testing with surface waves [J].12nd Internat.Conf. on Soil Mechanics and Found. Engin. Proceedings,1989,331-334.
[168]Gucunski N, Wood R D. Instrumentation for SASW testing [J]. In: Bhatia S K, Blaney G W,(Eds), Recent advances in instrumentation, data acquisition and testing in soil dynamics: Am. Soc.Of Civil Engin., Geotechnical special publication,1991,29:1-16.
[169]Hiltunen D R. Nondestructive evaluation of pavement systems by the SASW method [J].Geotechnical News,1991,9:22-25.
[170]Ali Z, Osama H. Inversion of SASW dispersion curves based on maximum flexibility coefficientsin the wave number domain [J]. Soil Dynamics and Earthquake Engineering,2005,26(8):735-752.
[171]Yordkayhun S, Giese R, Juhlin C, Cosma C. Shallow velocity-depth model using first arrivaltraveltime inversion at the CO2SINK site, Ketzin, Germany [J]. J. Applied Geophysics,2007,63:68–79.
[172]Yordkayhun S, Juhlin C, Norden B.3D seismic reflection surveying at the CO2SINK project site,Ketzin, Germany: A study for extracting shallow subsurface information [J]. Near SurfaceGeophysics,2009a,7:75–91.
[173]Yordkayhun S, Tryggvason A, Norden B, Juhlin C, Bergman B.3D seismic traveltimetomography imaging of the shallow subsurface at the CO2SINK project site, Ketzin, Germany [J].Geophysics,2009b,74: G1–G15
[174]Draganov D, Campman X, Thorbecke J, Verdel A, Wapenaar K. Event-Driven SeismicInterferometry With Ambient Seismic Noise [J].72nd EAGE Conference and Exhibition,,2010，G013.
[175]Bergmann P, Yang C, Luth S, Juhlin C, Cosma C. Time-lapse processing of2D seismic profileswith testing of static correction methods at the CO2injection site Ketzin (Germany)[J]. J.Applied Geophysics,2011,75:124–139.
[2]王伟涛,倪四道,王宝善.地球背景噪声干涉应用研究的新进展[J].中国地震,2011,27(1):1-13.
[3] Wang B, Zhu P, Chen Y, et al. Continuous subsurface velocity measurement with coda waveinterferometry [J]. J. Geophys. Res.,2006,113: B12313, doi:10.1029/2007JB005023.
[4] Yao H, R D van der Hilst, de Hoop M V. Surface-wave array tomography in SE Tibet from ambientseismic noise and two-station analysis: I-phase velocity maps [J]. Geophys. J. Int.,2006,166:732–744.
[5] Ritzwoller M, Shapiro N, Pasyanos M, Bensen G, Tang Y. Short-period surface wave dispersionmeasurements from ambient seismic noise in North Africa, the Middle East, and Central Asia,161–170[R].27th Seismic Research Review: Ground-Based Nuclear Explosion MonitoringTechnologies,2005.
[6] Florent B, Nikolai M S, Michel C, Slexandre N, Valérie F.3-D surfacewave tomography of thePiton de la Fournaise volcano using seismic noise correlations [J]. Geophysical Research Letters,2007,34: L02305, doi:10.1029/2006GL028585.
[7] Claerbout J F. Synthesis of a layered medium from its acoustic transmission response [J].Geophysics,1968,33:264–269.
[8] Rickett J, Claerbout J F. Acoustic daylight imaging via spectral factorization: helioseismology andreservoir monitoring [J]. The Leading Edge,1999,18:957-960.
[9] Wapenaar C P A, Draganov D, Thorbecke J, Fokkema J T. Theory of acoustic daylight imagingrevisited,2002[C].72nd SEG annual meeting, Expanded Abstracts,2269-2272.
[10]Scherbaum F. Seismic imaging of the site response using microearthquake recordings, Part I:Method [J]. Bull. Seis. Soc. Am.,1987a,77:1905–1923.
[11]Scherbaum F. Seismic imaging of the site response using microearthquake recordings, Part II:Application to the Swabian Jura, southwest Germany, seismic network[J]. Bull. Seis. Soc. Am.,1987b,77:1924–1944.
[12]Cole S. Passive seismic and drill-bit experiments using2-D arrays [D]. Ph.D. dissertation, StanfordUniversity,1995.
[13]Daneshvar M, Clay C, Savage M. Passive seismic imaging using microearthquakes [J].Geophysics,1995,60:1178–86.
[14]Baskir C, Wells C. Sourceless reflection seismic exploration [J]. Geophysics,1975,40,158.
[15]Duvall T L, S M Jefferies, J W Harvey, M A Pomerantz. Time-distance helioseismology [J].Nature,1993,362:430–432.
[16]Gilles P M, Duvall T L, Scherrer P H, Bogart R S., A subsurface flow of material from the Sunequator’s to its poles [J]. Nature,1997,390:52–54.
[17]Cassereau D, Fink M. Focusing with plane time reversal mirrors: an efficient alternative to closedcavities [J]. Journal of the Acoustical Society of America,1993,94(4):2373-2386.
[18]Fink M. Time reversal mirrors [J]. J. Phys. Rev.: Appl. Phys.,1993,26,1330-1350.
[19]Fink M. Time reversed acoustics [J]. Physics Today,1997,50:34-40.
[20]Roux P, Fink M. Green′s function estimation using secondary sources in a shallow waterenvironment [J]. Journal of the Acoustical Society of America,2003,113:1406-1416.
[21]Snieder R, Grêt A, Douma H, Scales J. Coda wave interferometry for estimating nonlinearbehavior in seismic velocity [J]. Science,2002,295:2253-2255.
[22]Snieder R. Coda wave interferometry and the equilibration of energy in elastic media [J]. Phys.Rev. E,2002,66:046615-1-046615-8.
[23]Snieder R, Hagerty M. Monitoring change in Volcanic Interiors using Coda WaveInterferometry:Application to Arenal Volcano, Costa Rica [J]. Geophys. Res. Lett.,2004,31(9): L09608.1-L09608.5.
[24]Grêt A R, Snieder R, Aster C, Kyle P R. Monitoring rapid temporal changes in a volcano with codawave interferometry [J]. Geophys. Res. Lett.,2005,32(6): L06304.1-L06304.
[25]Snieder R. The theory of coda wave interferometry [J]. Pureand Appl. Geophys,2006,163:455-473.
[26]Fink M.1992, Time-reversal of ultrasonic fields: Basic principles [J]. IEEE Transactions onUltrasonics, Ferroelectrics, and Frequency Control,39,555–566.
[27]Fink M. Time-reversal acoustics in complex environments [J]. Geophysics,2006,71(4):SI151–SI164.
[28]Lobkis O I, R L Weaver. On the emergence of the Green’s function in the correlations of a diffusefield [J]. Journal of the Acoustical Society of America,2001,110:3011–3017.
[29]Lerosey G, de Rosny J, Tourin A, Fink M. Focusing beyond the diffraction limit with far-fieldtime reversal [J]. Science,2007,315,1120–1122.
[30]Derode A, Larose E, Campillo M, Fink M. How to estimate the Green’s function of aheterogeneous medium between two passive sensors-Application to acoustic waves [J]. AppliedPhysics Letters,2003a,83:3054–3056.
[31]Derode A, Larose E, Tanter M, de Rosny J, Tourin A, Campillo M, Fink M. Recovering theGreen’s function from field-field correlations in an open scattering medium (L)[J]. Journal of theAcoustical Society of America,2003b,113:2973–2976.
[32]Blomgren P, Papanicolaou G, Zhao H. Super-resolution in time-reversal Acoustics [J]. J. Acoust.Soc. Am.,2002,111(1),230–48.
[33]Borcea L, Papanicolaou G, Tsogka C, Berryman J. Imaging and time reversal in random media [J].Inverse Problems,2002,18,1247–79.
[34]Borcea L, Papanicolaou G, Tsogka C. Coherent interferometric imaging in clutter [J], Geophysics,2006,71, SI165–SI175.
[35]Schuster G T. Seismic interferometry [M]. New York: Cambridge University Press,2009.
[36]Weaver R L, Lobkis O I. Ultrasonics without a source: thermal fluctuation correlations at MHzfrequencies [J]. Physical Review Letters,2001,87:134301-1–134301-4.
[37]Weaver R L, Lobkis O I.2002, On the emergence of the Green’s function in the correlations of adiffuse field: Pulse-echo using thermal phonons [J]. Ultrasonics,40,435–439.
[38]Campillo M, Paul A. Long-range correlations in the diffuse seismic coda [J]. Science,2003,299:547–549.
[39]Weaver R L, Lobkis O I. Elastic wave thermal fluctuations, ultrasonic waveforms by correlation ofthermal phonons [J]. J. Acoust. Soc. Am.,2003,113:2611–2621.
[40]Weaver R L, Lobkis O I. Diffuse fields in open systems and the emergence of the Green’s function[J]. J. Acoust. Soc. Am.,2004,116(5):2731–2734.
[41]Larose E, Derode A, Campillo M, Fink M. Imaging from one-bit correlations of wide-band diffusewavefield [J]. J. Appl. Phys.,2004,95:8393–8399.
[42]Larose E, Lobkis O I, Weaver R L. Passive correlation imaging of a buried scatterer[J]. J. Acoust.Soc. Am.,2006,119:3549–3552.
[43]Malcolm A, Scales J, Tiggelen B. Extracting the Green function from diffuse, equipartitionedwaves [J]. Physical Review E,2004,70:015601(R)-1–015601(R)-4.
[44]Roux P, Kuperman W A, the NPAL Group. Extracting coherent wavefronts from acoustic ambientnoise in the ocean [J]. J. Acoust. Soc. Am.,2004,116:1995–2003.
[45]Sabra K G, Roux P, Kuperman W A. Arrival-time structure of the time-averaged ambient noisecross-correlation function in an oceanic waveguide [J]. J. Acoust. Soc. Am.,2005,117:164–174.
[46]Sabra K G, P Roux, W A Kuperman.2005, Arrival-time structure of the time-averaged ambientnoise cross-correlation function in an oceanic waveguide [J]. Journal of the Acoustical Society ofAmerica,117,1:164–174.
[47]Shapiro N, Campillo M. Emergence of broadband Rayleigh waves from correlations of theambient seismic noise [J]. Geophys. Res. Lett.,2004,31: L07614.
[48]Shapiro N, Campillo M, Stehly L, Ritzwoller M. High-resolution surface-wave tomography fromambient seismic noise [J]. Science,2005,307:1615–18.
[49]Sabra K, Gerstoft P, Roux P, Kuperman W, Fehler M. Extracting time-domain Green’s functionsfrom ambient seismic noise [J]. Geophys. Res. Lett.,2005b,32: L03310.
[50]Gerstoft P, Sabra K, Roux P, Kuperman W A, Fehler M. Green’s functions extraction andsurface-wave tomography from microseisms in southern California [J]. Geophysics,2006,71:SI23–SI32.
[51]Larose E, Margerin L, Derode A, van Tiggelen B, Campillo M, Shapiro N, Paul A, Stehly L, TanterM. Correlation of random fields: An interdisciplinary review [J]. Geophysics,2006,71: SI11–SI21.
[52]Gudmundsson O, Khan A, Voss P. Rayleigh-wave group-velocity of the Icelandic crust fromcorrelation of ambient seismic noise [J]. Geophys. Res. Lett.,2007,34: L14314.
[53]Gouedard P, Stehly L, Brenguier F, Campillo M, de Verdiere Y, Larose E, Margerin L, Roux P,Sanchez-Sesma F, Shapiro N, Weaver R. Cross-correlation of random fields: Mathematicalapproach and applications [J]. Geophysical Prosp.2008,56:375–393.
[54]Nunziata C, De Nisco G, Panza G F. S-wave profiles from noise cross correlation at small scale [J].Engineering Geology,2009,105:3–4.
[55]Shragge J, Artman B, Wilson C. Teleseismic shot-profile migration [J]. Geophysics,2006,71:SI221–SI229.
[56]Schuster G T, Rickett J. Daylight imaging in V(x,y,z) media [R]. Utah Tomography andModeling-Migration Project Midyear Report and SEP Report at Stanford University,2000.
[57]Yu J, Schuster G T, Crosscorrelogram migration of IVSPWD data [J].71st Ann. Internat. Mtg.,SEG Expanded Abstracts,2001,456–459.
[58]Yu J, Schuster G T, Enhancing illumination coverage of VSP data by cross-correlogram migration,74th Ann. Internat. Mtg., SEG Expanded Abstracts,2004,2501–2504.
[59]Yu J, Schuster G T, Crosscorrelogram migration of inverse vertical seismic profile data [J].Geophysics,2006,71, S1–S11.
[60]Hornby B E, Yu J. Interferometric imaging of a salt flank using walkaway VSP data [J]. TheLeading Edge,2007,26:760–763.
[61]Sheng J, Schuster G T, Nowack R. Imaging of crustal layers by teleseismic ghosts [C]. AmericanGeophysical Union, Fall Mtg.,2002, S32C–0658.
[62]Sheng J, Schuster G T, Pankow K, Pechmann J, Nowack R. Coherence-weighted wavepathmigration of teleseismic data [C]. American Geophysical Union, Fall Mtg.,2003, S11E–0344.
[63]Sheng J. Migrating multiples and primaries in CDP data by crosscorrelation Migration [J].71stAnn. Internat. Mtg., SEG Expanded Abstracts,2001,1297–1300.
[64]He R, Hornby B, Schuster G.3D wave-equation interferometric migration of VSP free-surfacemultiples [J]. Geophysics,2007,72: S195–S203.
[65]Schuster G T, Yu J, Sheng J, Rickett J. Interferometric/daylight seismic imaging [J]. GeophysicalJournal International,2004,157:838–852.
[66]Wapenaar K, Draganov D, Thorbecke J, Fokkema J. Theory of acoustic daylight imaging revisited[C].72nd Annual International Meeting, SEG, Expanded Abstracts,2002,2269–2272.
[67]Wapenaar K. Synthesis of an inhomogeneous medium from its acoustic transmission response [J].Geophysics,2003,68:1756-1759.
[68]Wapenaar K, Thorbecke J, Draganov D. Relations between reflection and transmission responsesof three-dimensional inhomogeneous media [J]. Geophysical Journal International,2004,156:179–194.
[69]Wapenaar K. Retrieving the elastodynamic Green’s function of an arbitrary inhomogeneousmedium by crosscorrelations [J]. Physical Review Letters,2004,93:254301-1–254301-4.
[70]Wapenaar K, Fokkema J, Snieder R. Retrieving the Green’s function in an open system bycross-correlation: A comparison of approaches(L)[J]. Journal of the Acoustical Society of America,2005,118:2783–2786.
[71]Wapenaar K. Green’s function retrieval by cross-correlation in case of one-sided illumination [J].Geophysical Research Letters,2006,33: L19304-1-L19304-6.
[72]Willis M E, R Lu, X Campman, M N Toks z, Y Zhang, M V de Hoop. A novel application oftime-reversed acoustics: Salt-dome flank imaging using walkaway VSP surveys [J]. Geophysics,2006,71(2): A7–A11.
[73]Artman B. Imaging passive seismic data [J]. Geophysics,2006,71:177-187.
[74]Calvert R W, Bakulin A, Jones T C. Virtual sources, a new way to remove overburden problems[J].66th Mtg., European Association of Geoscientists and Engineers, Extended Abstracts,2004,P234.
[75]Bakulin A, Calvert R. Virtual source: New method for imaging and4D below complex overburden[C].74th Ann. Internat. Mtg., SEG Expanded Abstracts,2004,2477–2480.
[76]Bakulin A, Calvert R. The virtual source method: Theory and case study [J]. Geophysics,2006,71:SI139–SI150.
[77]Berkhout A, Verschuur D. Transformation of multiples into primary reflections,73rd Ann. Internat.Mtg., SEG Expanded Abstracts,2003,1925–1928.
[78]Berkhout A, Verschuur D. Imaging of multiple reflections, Geophysics,2006,71, SI209–SI220.
[79]Jiang Z, Yu J, Hornby B, Schuster G T. Migration of multiples [J]. The Leading Edge,2005,315–318.
[80]Jiang Z, Sheng J, Yu J, Schuster G T, Hornby B. Migration methods for imaging different-ordermultiples [J]. Geoph. Prosp.,2007,55,1–19.
[81]Guitton A. Shot-profile migration of multiple reflections [J].72nd Ann. Internat. Mtg., SEGExpanded Abstracts,2002,1296–1299.
[82]Shan G, Guitton A. Migration of surface-related multiples: Tests on Sigsbee2B dataset [J].74thAnn. Internat. Mtg., SEG Expanded Abstracts,2004,1285–1288.
[83]Meng F, Schuster G T. Autocorrelogram imaging of Friendswood RVSP data [R] Annual UtahReport,320-338.
[84]Xiao X, Schuster G T. Redatuming CDP data below salt with VSP Green’s function [J].76th Ann.Internat. Mtg., SEG Expanded Abstracts,2006,3511–3515.
[85]Xiao X, Zhou M, Schuster G T. Salt-flank delineation by interferometric imaging of transmitted Pto S waves, Geophysics,2006,71, SI197–SI207.
[86]Draganov D, Wapenaar K, Thorbecke J. Passive seismic imaging in the presence of white noisesources [J]. The Leading Edge,2004,23,889–892.
[87]Draganov D, Wapenaar K, Mulder W, Singer J. Seismic interferometry on background-noise fielddata [C].76th Annual International Meeting, SEG, Expanded Abstracts,2006,590–593.
[88]Draganov D, Wapenaar K, Thorbecke J. Seismic interferometry: Reconstructing the earth’sreflection response [J]. Geophysics,2006,71: SI61–SI70.
[89]Draganov D, Wapenaar K, Mulder W, Singer J, Verdel A. Retrieval of reflections from seismicbackground-noise measurements [J]. Geophysical Research Letters,2007,34, L04305-1–L04305-4.
[90]Draganov D, Campman X, Thorbecke J, Verdel A, Wapenaar K. Reflection images from seismicnoise [J]. Geophysics,2009,74: A63–A67.
[91]Artman B. Imaging passive seismic data [J]. Geophysics,2006,71, SI177–SI187.
[92]Wapenaar K, Draganov D, Snieder R, Campman X, Verdel A. Tutorial on seismic interferometry:Part1-Basic principles and applications [J]. Geophysics,2010,75(5):75A195-75A209.
[93]Dorman J, Ewing M. Numerical inversion of surface-wave dispersion data and crust-mantlestructure in the New York-Pennsylvania area [J]. J. geophys. Res.,1962,67:5227–5241.
[94]Lu R, Venkataraman A, Payne M, Zhang J. Application of noise interferometry to obtaintime-lapse velocity variations during a steam stimulation cycle at Cold Lake [C].79th AnnualInternational Metting, SEG, Expanded Abstracts,2009,1673–1677.
[95]Snieder R, Safak E. Extracting the building response using seismic interferometry: Theory andapplication to the Millikan Library in Pasadena, California [J]. Bulletin of the SeismologicalSociety of America,2006,96:586–598.
[96]Curtis A P, Gerstoft H, Sato R, Snieder R, Wapenaar K. Seismic interferometry—Turning noiseinto signal [J]. The Leading Edge,2006,25:1082–1092.
[97]Wapenaar K, van der Neut J, Ruigrok E. Passive seismic interferometry by multidimensionaldeconvolution [J]. Geophysics,2008,73(6):51-56.
[98]Vasconcelos I, Snieder R. Interferometry by deconvolution, Part1—Theory for acoustic waves andnumerical examples [J]. Geophysics,2008,73(1): S115-S128.
[99]Vasconcelos I, Snieder R. Interferometry by deconvolution, Part2—Theory for elastic waves andapplication to drill-bit seismic imaging [J]. Geophysics,2008,73(3):129-141.
[100]Schuster G T. Seismic interferometric imaging with waveforms [R]. Utah Tomography andModeling-Migration Project Mideyear Report,1999,121-130
[101]Schuster G T, Followill F, Lewis K J, Yu J, Liu Z. Autocorrelogram migration: Theory [J].Geophysics,2003,68(5):1685-1694.
[102]Mercier J, Bostock M, Baig A. Improved Green’s functions for passive-source structural studies[J]. Geophysics,2006,71: SI95-SI102
[103]Wang B, Zhu P, Chen Y, et al. Continuous subsurface velocity measurement with coda waveinterferometry [J]. J. Geophys. Res.,2006,113: B12313, doi:10.1029/2007JB005023
[104]房立华,吴建平,吕作勇.华北地区基于噪声的瑞利面波群速度层析成像[J].地球物理学报,2009,52(3):663-671.
[105]齐诚,陈棋福,陈颙.利用背景噪声进行地震成像的新方法[J].地球物理学进展,2007,22(3):771-777.
[106]汪承灏,魏炜.改进的时间反转法用于有界面时超声目标探测的鉴别[J].声学学报,2002,27(3):193-197.
[107]吴昊,张碧星,汪承灏.声波在固体板中的多径传播及其时间反转声场[J].声学学报,2005,30(3):215-221.
[108]陆铭慧,张碧星,汪承灏.时间反转法在水下通信中的应用[J].声学学报,2005,30(4):215-221.
[109]Wang B, Zhu G, Gao J. Joint interferometry imaging of walkaway VSP data [J]. AppliedGeophysics,2010,7(1):41-48.
[110]Hill S J. Geophysics bright spots [J]. The Leading Edge,2006,25:1082-1092.
[111]单新建,马瑾,王长林,柳稼航.利用差分干涉雷达测量技术(D-InSAR)提取同震形变场[J].地震学报,2007,24(4),413-420.
[112]朱新慧,孙付平.用甚长基线干涉测量数据检测冰期后地壳回弹[J].地球物理学报,2005,48(2):308-313.
[113]Wapenaar K, Fokkema J. Green’s function representations for seismic interferometry [J].Geophysics,2006,71(4): SI33–SI46.
[114]de Hoop A T. An elastodynamic reciprocity theorem for linear, viscoelastic media [J]. AppliedScientific Research,1966,16:39–45.
[115]de Hoop A T. Time-domain reciprocity theorems for acoustic wave fields in fluids with relaxation[J]. Journal of the Acoustic Society of America,1988,84,1877–1882.
[116]Aki K, Richards P G. Quantitative Seismology: Theory and Methods [M]. San Francisco:Freenman and Coopration,1980,146-177.
[117]何樵登.地震波理论[M].北京:地质出版社,1988.
[118]de Hoop A T. Handbook of radiation and scattering of waves [M]. Academic Press, Amsterdam,1995.
[119]Rayleigh J W S. The theory of sound, Volume II [M]. Macmillan and Co.1878.
[120]Fokkema J T, van den Berg P M. Seismic applications of acoustic reciprocity [M]. ElsevierScience,1993.
[121]Bojarski N N. Generalized reaction principles and reciprocity theorems for the wave equations,and the relationship between the time advanced and the time-retarded fields [J]. Journal of theAcoustic Society of America,1983,74,281–285.
[122]Bleistein N. Mathematical methods for wave phenomena: Academic Press Inc,1984.
[123]Knopoff L, Gangi A F. Seismic reciprocity [J]. Geophysics,1959,24,681–691.
[124]Aki K, Richards P G, Quantitative seismology, vol.1[M]. W. H. Freeman and Company,1980.
[125]Mehta K, Snieder R. Acquisition geometry requirements for generating virtual-source data [J].The Leading Edge,2008,27(5):620-629.
[126]Hardage B. Vertical seismic profiling: Principles [M]. Elsevier Science Publications Co., Inc.,2000.
[127]Blakeslee S, Chen S, Krebs J, Srnka L. TVSC: Twin VSP simulation of cross-well data, Astrategy for low-cost monitoring of EOR processes [J].63rd Ann. Internat. Mtg., SEG ExpandedAbstracts,1993,12:9–12.
[128]Cao W, Schuster G T. Natural redatuming of VSP data to virtual CDP data [R].2005UTAM Ann.Report, University of Utah,2005,231–248.
[129]Krebs J, Fara D, Berlin A. Accurate migration using offset-checkshot surveys [J]65th Ann.Internat. Mtg., SEG Expanded Abstracts,1995,14:1186–1189.
[130]Schuster G T. Reverse time migration=generalized diffraction stack migration [J].72nd Ann.Internat. Mtg., SEG Expanded Abstracts,2002,21:1280–1283.
[131]Brandsberg-Dahl S, Hornby B, Xiao X. Migration of surface seismic data with VSP Green’sfunctions [J]. The Leading Edge,2007,26:778–781.
[132]Xiao X. Local reverse-time migration with VSP Green’s functions [D]. Salt Lake City, Utah:University of Utah,2008.
[133]Minato S, Onishi K, Matsuoka T, Okajima Y, Tsuchiyama J, Nobuoka D, Azuma H, Iwamoto T.Cross-well seismic survey without borehole source [J].77th Ann. Internat. Mtg., SEG ExpandedAbstracts,2007,26:1357–1361.
[134]Hornby B, Yu J. Single-well imaging of a salt flank using walkaway VSP data [J].76th Ann.Internat. Mtg., SEG Expanded Abstracts,2006,25:3492–3496.
[135]Hornby B, Yu J, Sharp J, Ray A, Quist Y, Regone C. VSP: Beyond time-to-depth [J]. The LeadingEdge,2006,25:446–452.
[136]He R, Hornby B, Schuster G T.3D wave-equation interferometric migration of VSP multiples [J].76th Ann. Internat. Mtg., SEG Expanded Abstracts,2006,25:3442–3445.
[137]Willis M, Lu R, Campman X, Toksoz N, Zhang Y, de Hoop M. A novel application oftime-reversed acoustics: Salt-dome flank imaging using walkaway VSP surveys [J]. Geophysics,2006,71: O43–O51.
[138]Korneev V, Bakulin A. On the fundamentals of the virtual source method [J]. Geophysics,2006,71: A13-A17.
[139]Chakroun N, Fink M, Wu F. Time reversal processing in ultrasonic nondestructive testing [J].IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control,1995,42:1087–1098.
[140]Draeger C, Cassereau D, Fink M. Acoustic time reversal with mode conversion at a solid-fluidinterface [J]. Applied Physics Letters,1998,72:1567–1569.
[141]Fink M, Cassereau D, Derode A, Prada C, Roux P, Tanter M, Thomas J L, Wu F. Time-reversedacoustics [J]. Reports on Progress in Physics,2000,63:1933–1995.
[142]Fink M, Prada C. Acoustic time-reversal mirrors [J]. Inverse Problems,2001,17: R1–R38.
[143]Tsogka C, Papanicolaou G. reversal through a solid-liquid interface and super-resolution [J].Inverse Problems,2002,18:1639–1657.
[144]Morse P, Feshbach H. Methods of Theoretical Physics (Part I)[M]. New York: McGraw-Hill,1953.
[145]F rster A, Norden B, Zinck-J rgensen K, Frykman P, Kulenkampff J, Spangenberg E, Erzinger J,Zimmer M, Kopp J, Borm G, Juhlin C, Cosma C, Hurter S. Baseline characterization of theCO2SINK geological storage site at Ketzin, Germany [J]. Environmental Geosciences,2006,13:145–160.
[146]Martens S, Liebscher A, M ller F, Würdemann H, Schilling F, Kühn, M, Ketzin Group. Progressreport on the first European on-shore CO2storage site at Ketzin (Germany)-Second year ofinjection [R]. Energy Procedia,2011,4:3246–3253.
[147]Juhlin C, Giese R, Zinck-J rgensen K, Cosma C, Kazemeini H, Juhojuntti N, Lüth S, Norden B,F rster A.3D baseline seismics at Ketzin, Germany: The CO2SINK project [J]. Geophysics,2007,72: B121–B132.
[148]Derode A, Tourin A, Fink M. Ultrasonic pulse compression with one-bit time reversal throughmultiple scattering [J]. J. Appl. Phys.,1999,85:6343–6352.
[149]Dong S, Sheng J, Schuster J T. Theory and practice of refraction interferometry [J].76th AnnualInternational Metting, SEG, Expanded Abstracts,2006,26:3021–3025.
[150]Mikesell D, van Wijk K, Calvert A, Haney M. The virtual refraction: Useful spurious energy inseismic interferometry [J]. Geophysics,2009,74: A13–A17.
[151]Nichols J, Mikesell D, van Wijk K. Application of the virtual refraction to near-surfacecharacterization at the Boise Hydrogeophysical Research Site [J]. Geophysical Prospecting,2010,58:1011–1021.
[152]Rayleigh L. On waves propagated along the plane surface of an elastic solid [J]. Proc. Lond.Math. Soc.,1885,17:4-11.
[153]Haskell N A. The dispersion of surface waves on multilayered media [J]. Bulletin of theSeismological Society of America,1953,43,17-34.
[154]Babuska V, Cara M. Seismic anisotropy in the earth [M]. Academic Publishers,1991.
[155]Sheriff R E, Geldart L P. Exploration seismology (volume II): History, theory, and dataacquisition [M]. New York: Cambridge University Press,1983.
[156]Stokoe K H, Wright S G, Bay J A, Ro sset J M. Characterization of geotechnical sites by SASWmethod [M]. In: R.D. Woods, Editor, Geophysical Characterization of Sites: Oxford Publishers,1994.
[157]Aki K. Space and time spectra of stationary stochastic waves, with special reference tomicrotremors [J]. Bull. Earthq. Res. Inst.,1957,35:415–456.
[158]Matthews M C, Hope V S, Clayton C R I. The use of surface waves in the determination ofground stiffness profiles [J]. Geotechnical Eng., Proc.Inst.Civ.Eng,1996,119:84-95.
[159]Aki K, Richards P G. Quantitative Seismology [M]. San Francisco: W. H. Freeman and Company,1980.
[160]Dorman J, Ewing M. Numerical inversion of surface-wave dispersion data and crust-mantlestructure in the New York-Pennsylvania area [J]. Journal of Geophysical Research,1962,67:5227-5241.
[161]Song Y Y, Castagna J P, Black R A, Knapp R W. Sensitivity of near-surface shear-wave velocitydetermination from Rayleigh and Love waves [J]. Technical Program with Biographies, SEG,59th Annual Meeting, Dallas, TX,1989,8:509-512.
[162]Turner G. Aliasing in the tau-p transform and the removal of spatially aliased coherent noise [J].Geophysics,1990,55:1496-1503.
[163]Rix G J, Leipski E A. Accuracy and resolution of surface wave inversion [J]. In: Bhatia S K,Blaney G W,(Eds), Recent Advances in Instrumentation, Data Acquisition and Testing in SoilDynamics. Am. Soc. Civil Eng.,1991, pp.17-32.
[164]Nazarian S, Stokoe II K H, Hudson W R. Use of spectral analysis of surface waves method fordetermination of moduli and thicknesses of pavement systems [J]. Transport Research Record,1983,930:38-45.
[165]Sanchez-Salinero I, Rosset J M, Shao K Y, Stokoe II K H, Rix G J. Analytical evaluation ofvariables affecting surface wave testing of pavements [J]. Transport Res Rec,1987,1136:86-95.
[166]Sheu J C, Stokoe II K H, Roesset J M. Effect of reflected waves in SASW testing of pavements[J]. Transportation Res. Record,1988,1196:51-61.
[167]Stokoe II K H, Rix G J, Nazarian S. In situ seismic testing with surface waves [J].12nd Internat.Conf. on Soil Mechanics and Found. Engin. Proceedings,1989,331-334.
[168]Gucunski N, Wood R D. Instrumentation for SASW testing [J]. In: Bhatia S K, Blaney G W,(Eds), Recent advances in instrumentation, data acquisition and testing in soil dynamics: Am. Soc.Of Civil Engin., Geotechnical special publication,1991,29:1-16.
[169]Hiltunen D R. Nondestructive evaluation of pavement systems by the SASW method [J].Geotechnical News,1991,9:22-25.
[170]Ali Z, Osama H. Inversion of SASW dispersion curves based on maximum flexibility coefficientsin the wave number domain [J]. Soil Dynamics and Earthquake Engineering,2005,26(8):735-752.
[171]Yordkayhun S, Giese R, Juhlin C, Cosma C. Shallow velocity-depth model using first arrivaltraveltime inversion at the CO2SINK site, Ketzin, Germany [J]. J. Applied Geophysics,2007,63:68–79.
[172]Yordkayhun S, Juhlin C, Norden B.3D seismic reflection surveying at the CO2SINK project site,Ketzin, Germany: A study for extracting shallow subsurface information [J]. Near SurfaceGeophysics,2009a,7:75–91.
[173]Yordkayhun S, Tryggvason A, Norden B, Juhlin C, Bergman B.3D seismic traveltimetomography imaging of the shallow subsurface at the CO2SINK project site, Ketzin, Germany [J].Geophysics,2009b,74: G1–G15
[174]Draganov D, Campman X, Thorbecke J, Verdel A, Wapenaar K. Event-Driven SeismicInterferometry With Ambient Seismic Noise [J].72nd EAGE Conference and Exhibition,,2010，G013.
[175]Bergmann P, Yang C, Luth S, Juhlin C, Cosma C. Time-lapse processing of2D seismic profileswith testing of static correction methods at the CO2injection site Ketzin (Germany)[J]. J.Applied Geophysics,2011,75:124–139.
© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号 地址：北京市海淀区学院路29号 邮编：100083 电话：办公室：(+86 10)66554848；文献借阅、咨询服务、科技查新：66554700 |