CN1165561A

CN1165561A - 反射剪切波地震方法

Info

Publication number: CN1165561A
Application number: CN96191049A
Authority: CN
Inventors: 迈克尔·C·米勒; 利昂·汤姆森; 伊利娅·茨万科林
Original assignee: BP Corp North America Inc
Current assignee: BP Corp North America Inc
Priority date: 1995-10-06
Filing date: 1996-10-03
Publication date: 1997-11-19
Anticipated expiration: 2016-10-03
Also published as: CA2202168A1; AU699854B2; MX9704217A; WO1997013165A3; US5835452A; CN1148585C; EP0796441A2; NO971184D0; DE69617902D1; NO971184L; NO315767B1; EP0796441B1; PL321156A1; WO1997013165A2; AU7389996A; RU2181493C2; ATE210831T1

Abstract

一种在煤层构造层中测量及采用各向异性的方法。该方法包括以下步骤：通过对与慢剪切波极化方向重合的数据分量作一静平移以及对既不与慢剪切波极化方向重合又不与快剪切波极化方向重合的数据分量做上述静平移量一半的平移来分离指明了极化变化的各层；获得煤层的快极化轨迹包络的幅值量度；将煤层数据组旋转偏离主方向45℃；获得煤层至少一个混合/不匹配极化轨迹包络的幅值量度；以及获得煤层中快极化轨迹包络幅值与一个混合/不匹配极化轨迹包络幅值的比值。

Description

反射剪切波地震方法

本发明是属于一般的石油和天燃气勘探领域，特别是关于利用各向异性层上的表面反射剪切波进行地震勘探的方法。

识别和开发煤床的甲浣资源变得越来越重要。为达到这一目的，就需要有助于煤层及其厚度，几何构造，岩性学，结构和纹的预钻探识别方法。在特殊的地点，先后的钻探，开采和其它地质研究可以提供这些信息，但是，许可开采的区域常常位于已有矿区和芯孔位置之外。除此之外，这些煤层参数中横向迅速变化的蕴藏量也意味着这些传统的数据太稀疏以致不能提供出许可开采的煤资源的恰当面貌(见Henson，Jr.，H.和J.L.Sexton的“采用高解析地震反射方法预矿研究浅煤层”，地球物理学，56，(9)，第1494-1503页，1991)。

地震覆盖在揭示局部细节中可以起着重要作用。压缩波地震数据已经被采矿业用来提供有关煤层埋藏位置、厚度，几何构造和结构的信息，从而指导开采活动(见Gochioco，L.M.和S.A.Cotten的“采用高解析地震反射技术定位地下煤矿中的断层”，地球物理学，54，(12)，第1521-1527页，1989；Gochioco，L.M.的“来自薄床和煤层的转动效应和干涉反涉”，地球物理学，56，(8)，第1288-1295页，1991；以及Gochioco，L.M.的“在由煤层包围的薄层介质中的干涉反射的模型研究”，地球物理学，57，(9)，第1209-1216页，1992；poole，G.G.Duncan、L.Leung的“用于煤勘探的水平断错垂直地震断面勘查”，“SEG 1991 Exp.Abs.，第562-563页；以及urosevic，M.、B.J.Evans、以及P.J.Hatherly的“三维地震方法在煤层精细断层探测中的应用”，SEG 1992 Exp.Abs.，第254-256页)。通常采矿业的目的是避开完全的割理和完全断裂的煤层以便减少煤层甲浣渗透进煤矿中从而增加长墙操作的安全性。相反，高割理和裂纹区却对钻探孔的煤层甲烷生产有益(见Ting，F.T.C.的“在煤床中的割理的起源和隔开”，压力容器技术杂志，第624-626页，1977；Ayoub，J.、L.Colson、D.Johnston、J.Levine的“学会生产煤床沼气”，油田概览，第27-40页，1991)；Gash，B.W.的“用于煤床沼气生产的煤中‘岩性’的测量”，SPE 22909_，1991，第221-230页；以及Gash，B.W.、R.F.Volz和J.M.Corgan的“割理取向和约束压力在煤的割理多孔性、渗透性以及相对渗透性方面的影响”，1993年国际煤床沼气发展研讨会，1993年5月17日至21日，第247-255页)。

既然割理和裂纹的存在是所希望的，因此那些在钻探前可从远处识别这些裂纹和割理的方法将对煤层甲烷资源的优化经济开采提供潜力。

曾经发展了一种识别裂纹和裂纹区集气点的剪切波方法(见Thomsen，L.的“在方位角各向异性介质方面的反射地震学”，地球地理学，53，(3)，第304-313页，1988；以及Alford，R.M.，的“在方位角各向异性存在时的剪力数据”，SEG公约扩展的摘要，第476-479页，1986)，该方法并被应用于有裂纹的碳质碎岩中(见Mueller，M.C.的“采用多分量剪切波表面地震作为水平钻探的前驱波预测断裂强度的横向变异性”，地球地理学国际杂志，107，(3)，第387-408页，1991)。

在推广这些原始技术过程中，面临的主要挑战是来自煤层的三个特性，即：

煤层的薄层特性(相对于典型的地震波长度)

煤层极低的声阻(该声阻导致反系数比典型的沉积层大五倍)

对于对裂纹敏感的波，这一阻抗在强裂纹区域甚至减少的更低

这一事实。

当煤层在一个煤层系(CBS)中垂直堆叠时，这三个特性就同时存在，因此层间的多次复合可导致极有害的煤层反射并在CBS中产生不可解释的地震波震幅。(见Pietsch，K.和R。Slusarczyk的“高解析地震探测法在波兰煤矿中的应用”，地球地理学，57，(1)，第171-180页，1992；Gochioco，L.M.的“来自薄床和煤层的转动效应和干涉反射”，地球地理学，56，(8)，第1288-1295，1991；以及Gochioco，L.M.的“在由煤层围绕着的薄层介质中的干涉反射的模型研究”，地球地理学，57，(9)，第1209-1216页，1992)。因此，传统剪切波旋转和震幅技术的任何推广都应适应于这些地震波的挑战。

此外，裂纹引起的各向异性的理论，只对于单组环形裂纹的特殊情况研究得比较透彻。通常文献中讨论的单组环形裂纹引起的各向异性是由Hudson给出(见Hudsonm，J.A.的“断裂固体的整体特性”，MathProc.Camb.phil.Soc.，88，第371页，1980年)。Thomsen(见Thomsen，L.的“在多孔岩石中因裂纹对齐而产生的弹性各向异性”，欧洲协会勘探地球地理学公约摘要，53，第244-245页，(1991)；Geoph.Prospctg.，第＿-＿页，1995)指出如果裂纹是存在于一多孔岩石中(相对于实体而言)，那么多孔性对各向异性起作用，即使(按照定义)对各向异性没有更优先的方位。这一现象的发生是由于在地震频率点，裂纹和孔中流体间的流体压力的平衡。这一现象的普遍性由Rathore等的试验所证实(见Rathore，J.S.、Fjaer，E.、Holct，R.M.、Renlie，L.的“具有受控制断裂几何形状的人工合成沙夺的声学各向异性”，Geoph.Prospctng，，第＿-＿页，1995)。

因此，当前对裂纹引起的各向异性的认识明显不充分，至少需要推广到水平长、垂向短的正交裂纹组。相信带裂纹的集气点的直接探测要比间接探测更可取。Thomsen通过分析地震剪切波的极化方向来直接探测方位各向异性这一课题作过一个介绍(见Thomsen，L.的在方位角各向异性介质的反射地震学，53，(3)，第304-313页，1988，以及所引用的参考资料)。裂纹的间接探测(例如，采用超常的P波方法，探测近地表下构造并假设这些构造一定引起裂纹)，由于其在特点，很难令人满意。当然，这种直接探测依整于裂纹引起的方位各向异性的探测。然而，采用P-波进行方位各向导性的探测在历史上其它裂纹工作中(例如，Austin Chalk)还没有成功过；因此，对在CBS范围内的应用能否成功还不容乐观。

由于大部分岩石是方位各向导性的，因此，一般来说，垂直传播的剪切波并不按照波源建立的极化来传。而是，向量分解成两个垂直传播的波，两个波(在相当近似的程度上)相互正交极化，且传播速度不同。一般来说，这两个波极化方向既不是横线方向也不是顺线方向；而是由岩石不是由源确定的特殊极化方向。(熟悉这门学科的人会知道自然界很少存在“SH”和“SV”波；这两种波主要存在于教课书中！)

在最简单情况下(其中方位各向异性是由另外的各向同性的基质中的单组垂直裂纹或割理引起)，垂直传播波的两个特殊的极化方向是沿着裂纹并垂直于裂纹的。当然，两个都与传播方向成横向并且传播速度不同，因此剪切波分裂开。各向异性(即，在这种情况中根据方向的不同速度的不同以及极化方向的不同)可定义为：

γ = \frac{V^{| |} - V^{&perp;}}{V^{&perp;}} - - - - (1)

(即垂直传播的剪切波其与裂纹平行(‖)的极化波和与裂纹垂直(⊥)的极化波速度的相对差)由于V¹比V^⊥大，γ为正。选择的记号与Thomsen一致(见Thomsen，L.的“弱弹性各向异性”，地球地理学，51，(10)，第1954-1966页，1986)。

两个波向下独立传播，在交界面(这里假设为水平)反射开并回到表面，在表面记录下两个波。一般来说，横线接收器可接收到两个来波；顺线接收器也同样可做到。一般来说，无论源方向是顺线还是横线，都有上述结果；源的方位仅影响依赖于向量三角形分解的两个模式的相对激励。

常规的多源/多接收器(MS/MR)采集(具有顺线和横线，顺线和横线接收器)可得到四个轨迹，这些轨迹构成带有依赖于时间的分量的2×2张量S。选择“x”作为顺线方向(由于通常把X-Z作为横截面方向，这样显得更自然，)选择“y”作为横线方向，其指标为通常的方式(x，y)＝(1，2)，那么顺线源，顺线接收器的轨迹记作s₁₁(t)。那么横线/横线轨迹就为S₂₂(t)，不匹配的源和接收器的轨迹，为s₁₂(t)和s₂₁(t)，位于非角线元。

一般说，如果两种剪切模式间的延迟不可忽略但比波列持续时间小，那么(两种模式间)复杂的干涉就会出现在四个轨迹的每一个上，导致对轨迹解释的困难。但是在源与岩石的两个特殊方向重合的特殊情况，只有一个波被激起。例如，如果顺线源碰巧与“快”方向平行，那么只有快模式被激起，它被顺线接收器记录而不匹配接收器接收不到信号。另一模式只能被横线源激起并且只能由横线接收器接收。

由于四个观测到的轨迹sij(通常不可解释)是一张量，它可以通过张量旋转操作；旋转到一与岩石特殊方向重合的新座标系中。Alford(见Alford，R.M.的“在方位角各向异性存在时的剪力数据”，SEG公约扩展摘要，第476-479页，(1986)和美国专利4,817,061；4,803,666；5,029,146以及4,903,244)，指出，通过MS/MR旋转技术(该技术以他的名字命名)，可以(从四个不可解释的轨迹中)计算出：

^*一个“快”轨迹(标记为22，图1)，记录了一个不受裂纹影响的剪切波(在一组垂直平行裂纹的简单情况下)

^*一个“慢”轨迹(标记为11，图1)，记录了一个受裂纹强烈影响的剪切波；

^*两个近似为零的轨迹(标记为12和21，图1)。

Thomsen(见Thomsen.L.的在方位角各向异性介质的反射地震学，地球地理学，53，(3)，第304-313页，1988年)揭示了一种向量旋转方法。Thomsen给出基本Alford方法和一个只采用一个极化源的备选方法的推导(即一个单元/多接收器，SS/MR技术，也见美国专利4,888,743和4,933,913)

这些旋转方法假设在向下达到反射层整个深度上(可能除去几个薄层)只有单一的各向异性轴方向。这意味着剪切模式只分裂一次。假设了垂直传播，然而它是典型地应用于堆叠轨迹的。这需要进一步假设这样的堆叠轨迹是减噪的垂直入射轨迹的合理近似。使旋转方法可能正确的条件现在还不甚了解。

然而，Willis(见Willis，H.A.、Rethford，G.L.、Bielanski，E.的“方位角各向异性一剪切波数据质量上的发生和影响”，SEG 1986 Exp.Abs.，第479页，1986)指出，找出一个单一的(不随深度变)不确定性在±10°左右，并满足上述结果的旋转角θ通常是可能的。典型地，将采集的四分量数据组(即顺线，横线/顺线，横线接收器)，用不同的角度(例如θ＝15°，30°，45°，60°，75°等)旋转，选择出旋转后使非对向线分量极小化的角度。自动处理过程是可取的，但并不真真需要。

得到(即旋转得到)的两个对角轨迹(“主轨迹”)通常都象P波轨迹一样容易解释，(即，它们通常表现出好的反射层连续性，即使是大深度情况下。)与好的P-波数据一样，这些轨迹可用传统的剪切波处理方法(例如：岩性学识别，亮点解释等)来处理。此外，两种轨迹间的比较通常可显示出一小的相应反射层间的时间变化的延迟。在粗糙化的间隔上，延迟的增加作出了这一粗糙间隔上平均各向异性的一个量度，因此给出该间隔上裂纹的一个量度。

因此，可以说，旋转以后，每个对角轨迹可让人透过裂纹来勾划出地质构造并利用剪切数据达到“传统目的”。另一方面，两个对角轨迹的差别也让人看到了裂纹。这一差别(两个相应反射层间延迟时间的增加)的提供的一种各向异性(从而是裂纹)的量度的空间分辨度必然较差。

图l显示了这种方法在非常厚(70m)的煤层中应用的可能情况主岩层是各向同性的，γ＝10％，勘察线垂直于于煤层裂纹。这当然不是一个真实的构造，但是这一例子说明了前面用于识别粗糙化层的各向异性的方法的原理(见Martin，M.和T.Davis的“剪切波双折射：估算断裂矿藏的新工具”，前沿，6(10)，第22-28页，1987)。第一个反射波(在0.25秒到达)标明了煤层的顶部位置；由于这些波路经完全在碎覆盖层(这里假设为各向同性)中，所以两个主轨迹线11和22同时到达。第二个反射波(从厚煤层底部反射的波)中，“快”22轨迹到达的略早，表明煤层是方位各向异性。当然，以毫秒计的延迟时间Δt取决于平均各向异性与层厚度(用毫秒表示)的积；对于下面讨论的薄煤层系也有同样的结果。在许多应用中，各向异性要比图1计算中采用的要小得多，因此相应的垂直分辨度非常差。

Thomsen(见Thomsen，L.的方位角各向异性介质的反射地震学，地球地理学，53，(3)，第304-313页，1988)提出分裂剪切波正入射幅值的差异给出了反射层处各向异性(从而，裂纹强度差)突变的高可解量度。Thomsen的建议是一进步，尽管其幅值(相对于到达时间)的不可靠性是众所周知的，这是因为传播时间方法的总的成功(通过旋转)证明各种分量的幅值在相互比较时的确有意义，虽然其绝对幅值本身没有意义。Thomsen的建议由Mueller(见Mueller，M.C.的“采用多分量剪切波表面地震作为水平钻探的前驱波预测断裂强度的横向可变性”，Geoph.J.Intlo，107，(3)，第387-408页，1991)所证实，Muller发现在Texas中部的Austin Chalk，裂纹的高强度区可通过慢极化段的模糊点识别。

Thomsen的建议，在其最早的情形中，是基于正入射剪切波的反射。考虑一位于无裂纹的入射岩层(下标为1)下的具有垂直向裂纹的反射岩层(下标为2)，快模式(极化的‖于裂纹)的反射系数为1。

R_{| |} = \frac{P_{2} {V_{2}}^{| |} - P_{1} V_{1}}{P_{2} {V_{2}}^{| |} + P_{1} V_{1}} - - - - (2)

剪切速度V₂ ¹(为快剪切波，极化方向平行于裂纹)不受裂纹的影响，因此，对于Austin Chalk例子，V₂ ¹是比V₁大的数，为上层页岩的剪切速度。此外，在这一范围内，P2＞P100因此，在Mueller的证实文献(见Mueller，M.C.的“采用多分量剪切波表面地震作为水平钻探的前驱波预测断裂强度的横向可变性”，Geoph.J.Intl.，107，(3)，第387-408页，1991)中，R₁是一大的(负)数。方程2中前面的负号是由于剪切位移的习惯符号引起(也参见Aki，K.和P.G.Richards的量子地震学：理论和方法，W.H.Freeman & Co.，旧金山，1980)。

同样，慢模式(极化的⊥于裂纹)的反射系数为：

R_{&perp;} = - \frac{P_{2} {V_{2}}^{&perp;} - P_{1} V_{1}}{P_{2} {V_{2}}^{&perp;} + P_{1} V_{1}} = R_{| |} + γ / 2 - - - - (3)

其中γ为剪切波各向异性。剪切速度V₂ ^⊥由于裂纹小于V₂ ¹，因此V₂ ^⊥接近页岩速度V₁，(例如，前面引用的Mueller的Austin Chalk情况)。因此，在裂纹强的那些地方，R_⊥的绝对值较小(即较小的负数)。如果γ在10％的量级，R₁大约为5％，那么R_⊥约为零(即，大约100％的衰减)。Mueller在慢(⊥)段中得到的“模糊点”，识别出强裂纹的集气点，该点被以后的水平钻探证实。

此外，封闭在碎沉积岩中的薄煤层是前面没有提到的一种特殊情况。传播时间方法和反射幅值方法均不完全令人满意。传播时间方法由于薄层几何构造常常导致煤层系内短延迟回响间的复杂干涉而不能令人满意，因此，除一些非规则的CBS外，顶部和底部的反射均不清晰。因此，测量到的煤层系上积累的时间延迟不可靠。反射幅值方法的失败是由于薄的煤层构造，也是由于煤层的非裂纹蓄积速度小于顶层岩石的剪切速度等其它原因。

因此，需要找到一种用于探测封闭在或镶嵌在快沉积碎岩中的薄煤层示的增强的方位各向异性(从而增强的裂纹)区的地震特征。这一特征应不依赖于薄煤层的细部几何构造，也不应依赖于从原煤层上限定好的反射。此外，还应当认识到煤层系内内部多复合Coda的复杂性。

本发明总的目的是提出一种利用煤层系(CBS)表面反射剪切波特征的过程和方法；

然而本发明的另一目的是提出一种利用煤层系各向异性特性来识别自然裂纹区的过程和方法。

本发明的一个目的是识别一个可用来表证CBS各向异性特征的参数，该数对煤层中各向异性的量值敏感，对煤层系内各向异性材料的组成敏感，对各向异性层的厚度敏感，并且实际上不依赖于几何构造的细节。

本发明的另一目的(从广义上说)是提供一种识别和勾画地表层裂纹的过程和方法，从狭义上说也就是识别和勾画煤层甲烷勘探中的集气点的过程和方法。

本发明再一个目的是提供一种补尝CBS的覆盖层的各向异性的方法。

本发明的一个特殊目的是揭示煤的特性和其在煤层系中的几何构造是如何影响剪切波的特征以及这些特征的模型是如何导出补尝煤层系覆盖层各向异性和探测及评估割理方位的过程的。

按照本发明，实现了一种处理方位各向异性随深度变化的地层上的表面反射数据的方法。在本发明的一个实施例中，这种方法包括下列步骤：从一组多源和多接收器采集一剪切波数据堆积；将每个剪切波数据堆积的轨迹旋转一相应于数据堆积方位与最浅层方位各向异性方位之差的角度；确定最浅层各向异性方位的垂直变化方位。测量线层方位各向异性程度；完成数据堆积的层分离以获得最浅层的主时间序列，其中层分离包括考虑方位各向异性而修正数据堆积的慢极化轨迹和数据堆积的混合极化轨迹的步骤；对最浅层下面的一层重复上述步骤；将上面形成的数据旋转一相应于下层的方位各向异性方向与最浅层方位各向异性方向之差度以识别下面层的方位各向异性主方向。

在本发明的一个特殊实施例中，分离方法包括下述步骤：识别与慢剪切波极化方向重合的数据分量，识别与快剪切波极化方向一重合的数据分量以及识别与慢剪切波极化方向和快剪切波极化方向均不重合的不匹配数据分量，确定快剪切波轨迹和慢剪切波轨迹间的明显的延迟，确定快剪切波和慢剪切波的自然极化方向，对与慢剪切波极化方向重合的数据分量做一静平移；对既不与慢剪切波极地方向重合也不与快剪切波极化方向重合的分量做上述静平移量一半的平移；对时间小于层底部的数据组分量静噪。

本发明还揭示了一种量度和评估象煤层系这样的地表层中各向异性的方法。在本发明的一特殊实施例中，这一方法包括步骤：采集煤层系上的一个多分量表面反射地震剪切波数据组；将指明了极化变化的最浅层深度上的上述数据组的所有数据旋转一由一煤层上指明了极化变化的最浅层所确定的方位角以建立该煤层中各向异性主方向；将指明了极化变化的最浅层深度下的数据组的所有数据从该煤层各向异性主方向旋转开大约45°；获取旋转数据的混合不匹配极化轨迹中至少一个的包迹幅值的一个平均量度。

本发明的一个特殊实施例，揭示一种量度碳质层系中各向异性的方法。该方法包括下列步骤：通过对与慢剪切波极化方向一致的分量做一静平移以对与慢剪切波极化方向和快剪切波极化方向均不一致的分量做上述静平移量的一半的平移来分离极化变化指明了的层；将主时间序列数据组旋转一相应于最浅层方位各向异性方向与后一层方位各向异性方向之差的角度。其后，本发明的方法还可包括一附加的获得快极化轨迹包迹幅值量度的步骤；将该数据组从主方向上旋转开45°；获取至少混合/不匹配极化轨迹中之一个的包迹幅值的量度；获取快极化轨迹包迹平均幅值与混合/不匹配极化轨迹之一的包迹幅值的比值。

本发明所讨论的情形(即封闭在碎沉积岩中的薄煤层)不同于前面所述两种众所周知的方法，(即传播时间方法和反射幅值方法)，而是两种方法的改进。传播时间方法不完全令人满意，因为薄层几何构造常导致煤层系内短延迟回响(“桩一腿”复合)间的复杂干涉，使得顶部和底部的反射(一般情况下)都不清晰，(虽然在一些CBS的特列下可能清晰)。因此，量测到的煤层系中传播所积累的时间延迟不可靠。反射幅值方法失败的原因是由于薄层几何构造，也还由于在这种情况下煤层的非裂纹蓄积速度小于覆盖层岩石的该速度。并且，方程(2)中P₂＜P₁，使得R₁很大且为正。对于⊥-极化，V₂ ^⊥甚至更小，使得R_⊥甚至比R₁还大。因此，裂纹区的集气点导致慢⊥一段的亮点而不是模糊点。此外，由于无裂纹的(即固有的)反射系数如此大(R₁约为5％，图1)，使得小的γ值不能显著改变幅值只能以一个定量方式改变。这一效应在图1中很明显，其中顶部反射层的幅值对两个极化方向都非常小。

本发明的许多其它优点，从下面本发明的评细说明，实施例，权利要求书，以及附图中就可明显看到。

图1为从各向同性基质岩中的有裂纹煤床上采集的四分量多源/多接收器(MS/MR)数据的例子。

图2为Narellan模型煤层的几何说明。

图3表示的为沿图2所示Narellan模型煤层各向异性主轴所采集的四分量数据组。

图4表示的为沿图2所示Narellan模型煤层主断裂走向的45°方向上所采集的四分量数据组。

图5描述的是一个特殊的有兴趣的区域的CBS特性较正曲线

图6是Silesian模型煤层的几何描述；

图7表示的为沿图6 silesian模型煤层裂纹方向的45°轴方向采集的四分量数据组。

图8是San Juan模型煤层的几何描述。

图9表示产的为沿图8 san Juan模型煤层断裂方向的45°轴方向采集的四分量数据集。

图10表示的是从带有轻向异性并具有与煤层裂纹走向成30°角的有裂纹的覆盖层的图2 Narellan模型煤层中计算出的零偏移反射假想数据。

图11和12表示的是图10的原始数据相对于初始座标系分别旋转20°和30°度的结果。

图13描述的是进行了覆盖层补偿后的图12的数据。

图14表示了图13轨迹的旋转结果以寻找各向异性煤层主轴。

图15表示的是为得到参数R_a，将图14数据正则旋转45°的结果。

图16A，16B，16C表示的是原始数据，优化旋转数据和旋转45°的分离数据。

附图所表示以及下面将要详细描述的是本发明的一个特殊实施例，但本发明可以有许多不同形式的实施例。应当认识到，这里揭示的内容是从例举本发明原理的角度考虑的，并不想使本发明仅限于所例举说明的具体实施例。

如图2所示，煤层几何构造近似地模拟了澳大利亚新南威尔士(NSW)Narellan芯孔处最厚的40层煤床结构。最厚的煤层位于地下大约3000英尺的页岩和沙石处。在此盆地其它地方以及其它有兴趣的盆地(例如，San Juan，Silesia)具有不同的煤层结构(见图6和图8)。

沿图2煤层主裂纹带垂直方向，采集的一个MS/MR剪切线的四分量数据组结果(仅垂直传播)表示在图3中。假设各向同性覆盖层γ＝20％并且x轴垂直于煤层的断裂带。最快的极化轨迹22由煤层的顶层反射(接近0.5秒)，接着是从煤层内部反射的复杂的混合回响，而最底下煤层的反射则不清晰。本文中，虽然在各煤层内部有复杂的回响，但它们都有同样的极化。对应于慢极化轨迹11的反射到达时间是相同的(因为在本模型中覆盖层是各向同性的)，并且也跟随着一个复杂的回响(细节不同)。不谐调轨迹12和21在这种特殊的采集环境中如上所示均为零。如果没有一点裂纹，这些轨迹将类似，而不管线的方向如何，并且所显示的两个轨迹是一致的而不只是类似。

沿与主裂纹走向相交45°的线上，同样的采集结果表示于图4。顺线/顺线轨迹11含有快模式和慢模式二者的成份，信号比图3更复杂。

横线/横线轨迹22也有与上相同的特点。直观来说，图3与图4最明显的区别是不匹配轨迹的强信号，该信号总体幅度可与对角线轨迹的信号相比。这一现象的出现是因为快模式和慢模式都在煤层内回响，但它们不同的速度使这些回响相互不匹配。由于这些不匹配的回响，不匹配向量分量不能自动的相互抵消，结果造成一个强的合成信号。对于这一方位的线，两种模式在顺线和横线两个方向都有向量分量，因此都显示出不匹配轨迹。

在另一文章中，Alford(见Alford，R.M.的在方位角各向异性存在时的剪切数据，SEG公约扩展摘要，(1986)第476-479页)指出MS/MR数据集中不匹配轨迹的高能量是方位各向异性的灵敏的指示。在那里所指的各向异性位于覆盖层中。这里，覆盖层是各向同性；所指的各向异性位于煤层内。因此这些不匹配轨迹的能量值是煤层各向异性的量度，因此也就是断裂强度的量度。

当然，实际上不必沿与(初始未知)的裂纹走向成45°的线上采集数据。可以沿任一方便的线采集数据并对数据做数值旋转以找到能产生最近似图3结果的角度。这一方法建立了裂纹走向相对地图坐标系的方向。接着旋转45°就产生图4的等效结果，这里不匹配轨迹的幅度成为断裂强度的一个量度。我们可以预想断裂强度将各处不相同，因此这为直接探测断裂点提供了一种量度，一种地震特征或一种诊断方法。特别是，不匹配轨迹12或21(位于45°)的包迹幅度可看作是煤层平均各向异性(和厚度)的一个好的量度。这一不匹配轨迹是沿各向异性主轴采集的主轨迹11和22(见图3)之差的直接函数。在没有方位各向异性情况下，主轨迹将是一致的，并且45°方向的不谐调轨迹是零。为定义一个无量纲的平均CBS各向异性的量，将完全依赖于方位各向异性(见图4)的不匹配轨迹的包迹幅值A₁₂(45°)或A₂₁(45°)用最佳方位的快轨迹的包迹幅值A₂₂(0°)(该幅值完全独立于这样的各向异性，而依赖于煤层几何结构来正则化：

当然，幅值是CBS区域的平均值。只要Ra值大就意味着是一块高各向异性(因此高断裂强度)或大煤层厚度或二者皆有的区域。

模拟结果

剪切波地震在煤层断裂探测上的应用标志着分离剪切波技术的一种新的适用领域。虽然上述这种方法是基于现场证明了的物理原理，但数值模拟可以深入了解这种效果的强度，它与煤层断裂强度间的关系以及可能由覆盖层引起的失真。

下面讨论的人工剪切波地震反射是从Sydeny，silesian和San Juan盆地的典型煤层采集获得。模型是从已有的声波和强度曲线建立。虽然这些曲线提供了煤层的几何构造以及密度和垂直的P波和S波的速度，但煤层断裂系统的参数仍很难量化。对煤层的穿透性起主要作用的贯穿裂纹(大的割理)太大不能通芯样的测量来描述，因此模型中采用的断裂特性必须假定。

在模拟中假设煤层的各向异性是由在所有煤层中具有相同方位的单个平行垂直裂纹系引起。在煤层中通常存在的第二个割理系所引起的复杂性将被忽。这一假设的结论将被检验。

模型中另一忽略的效应是薄的水平层叠，这一层叠(在没有裂纹情况下)在地震频率上导致一种有效的具有垂直对称轴的横向各向同性介质(垂直横向各向同性即VTI)。对于垂直传递的剪切波VTI的存在可以合理地忽，但对于斜向入射波VTI会变得更重要。

在大部分例子中，覆盖层模拟成均匀的各同同性层。但是对于具有方位各向异性的覆盖层(或许方位各不相同)也做了考虑，并提出了这种覆盖层对剪切波反射影响的一种修正方法。

众所周知，纯各向同性基岩中的单组平行裂纹使介质横向各向同性，其对称轴垂直于裂纹平面。在现在的垂直裂纹情况下，对称轴为水平向(水平横向各向同性，即HTI)。均匀横向各向同性的介质可完全由其密度以及由Thomsen导出的五个弹性参数来描述(见Thomsen，L.的“弱弹性各向异性”，地球地理学，51，(10)，第1954-1966页，1986)：

^*速度V_P0(对于P波)和V_S0(对于S波)方向为对称轴方向(在这里为水平向)

^*三个无量纲各向异性系数(ε，δ和γ)

本文中垂直传播的剪切波的速度和幅值主要依赖于两个参数：S波速度(V_S0)和各向异性参数γ，这些参数确定了垂直入射波剪切分高度(即方程1)。应当注意在弹性平面波的近似中，其它三个参数对垂直传播波S波没有影响。对于更真实的点源辐射情况，比率V_P0/V_S0，系数ε和δ通过缓慢表面的形状影响剪切波幅值。然而，这些影响可认为是次等重要的。因此，煤层各向层性的程度可由参数γ来描述(即，说20％各向异性“意指γ＝0.2等等)。附录A证明了γ正比的断裂强度的值。

在所有模型中，下列参数用于基质岩和煤：

基质岩(沙石)：V_P＝4.5km/s，V_S＝2.65km/s，ρ＝2.37g/cm³

各向同性煤：V_P＝2.3km./s，V_S＝1.045km/s，ρ＝1.39g/cm³

各向异性煤：V_P0＝2.3km/s，V_S0＝883km/s，ρ＝1.39g/cm³

γ＝0.2；ε＝0；δ＝-0.16

模拟采用基于反射方法的全波形程序来完成(参见Garmany，J.的“在分层介质中弹性动力学特征解析的一些特性”，Geoph.J.R.A.S.，78，565-569，1983)。这一程序计算位于自由表面的实际点源的反射响应(包括所有可能的复合与转换)下面讨论的结果代表一个来自具有不同方向的点水平力的零偏移反射。采用的是拐角频率为10，20，90和120Hz的Ormsby弱波。图1中的以形状接近于波源脉冲。采用有限元程序对其中一个模形的计算检验了反射程序的精确性。两种方法的结果表明是近似的，因此证明了我们的方法的精确性。

下面的讨论主要集中在Narellan模(图2)，该模型表示了Sydney盆地一个芯孔处煤床的几何结构。该模拟得出了下列结果：

1.参数R_a对煤层的各向异性程度非常灵敏。参照图2的Narellan模型，保持模型的几何结构和基岩参数不变，煤层各向异性程度(相应地，断裂强度)改变了。对于小的和中等的各向异性，参数R_a实际上正比于各向异性γ，对于γ＝0.15，R_a达到了23％。虽然，对于高各向异性，曲线开始平缓，但做为可操作的有用的各向异性的量度，可以认为R_a(在最重要的范围γ＝0到25％内对各向异性足够灵敏。无疑，对于γ＞0，1到0.15，R_a的值超过了噪声水平，可以可靠地在现场测量到。

2.参数R_a对于煤的各向异性量，煤层内各向异性材料的裂纹以及各向异性煤层的厚度相当灵敏，但实际上独立于煤层的几何细节(即煤层的位置)。复合干涉特性不仅由煤层的弹性参数确定，而且由煤层的几何结构确定。如果煤层与主波长相比非常薄，并且互相非常接近，各层特性就会象一单个的均匀的层一样，其有效参数由各层的弹性参数的平均值确定，而不依赖于各层的位置细节。然而，至今所考虑的模型还没有超过这一长波长限制。

例I

为测试参数R_a对层结构细节的灵敏度，将Narellan模型(图2)每个煤层向上或向下(随机)移动1m，而不改变层的厚度以及其弹性参数。虽然每个波列的形状改变了，但参数R_a保持相同(4％以内)。因此R_a作为各向同性的最度实际上不依赖于CBS的任何几何细节。

例II

在另一测试中，所有煤层的厚度(在γ＝0.2的Narellan模型中)都增加20％。新模型的参数R_a相对原始模型的28.3％增加到40.5％。这并不感到奇怪，因为各向异性材料的裂纹的增加导致CBS内快和慢模式间回响延迟的增大。因此，45°线上的不匹配轨迹的幅值变高(虽然各向异性煤层的裂纹对R_a的影响不是光滑的)。

当然，很有可能对于CBS内不同的煤层其各向异性的量也在变化。这些变化的影响可通过假定Narellano模型中所有奇数层(由上数1，3，5，……)γ＝5％，所有偶数层(2.4.5……)γ＝20％来检验。这样模型的R_a值计算结果为11％，该结果更接近于均匀的5％的各向异性的结果(R_a＝8.7％)，而不是均匀的20％的各向异性结果(R_a＝28.3％)。这一现象可用前三层煤层(特别是第一层，Bulli层)对不匹配轨迹的最大振幅的影响来解释。从原理上讲，通过在移动的时间窗来计算参数R_a以获得原始的依赖于深度的各向异性量度来改进所提出的方法是可行的。

3.这一技术可以应用于其它领域。曾经对从波兰的Silesian盆地以及科罗拉多和新墨西哥的San Juan盆地分别地的芯孔记录所确定的煤层几何构造所表示的模型做过研究。在两种模型中覆盖层都假设为各向同性；为了比较，煤层中S-波的各向异性假设为20％。

Silesian模型(图6和7)的剪切波特征非常近似于Narellan模型的剪切波特征(见图2)，其参数R_a接受30％。对不匹配分量起主要作用的是靠近上部的几个煤层。

San Juan模型(图8和9)有点不同于前两个模型。它含有几个相对厚的煤层，该煤层引起“coda”波列中的低频反射。由于在这些厚层中快和慢模式间时间延迟的增大，其对不匹配分量的作用也相对增大。与煤层各向异性程度相同的Narellan模型相比，R_a参数由28.3％达到了47％。不匹配轨迹上的最大能量是由最厚的几个煤层产生。这一结果进一步证实了R_a值不仅依赖于各向异性的程度而且依赖于CBS内各向异性材料的裂纹的观察结论。在定量分析中，对任何特殊区域的CBS特征建立类似于图5的校正曲线可能是最好的方法。

4.覆盖层的各向异性的影响可以修正。在前面所有例子中覆盖层都假设为纯各向同性。然而，由于上表层中普遍存在着方位各向异性(见Willis，H.A.、Rethford，G.L.，Bielanski，E.的“方位角各向异性一剪切波数据质量的发生和影响”，SEG 1986 Exp.Abs.，第479页，1986)，很可能CBS上的碎屑层也具有轻度裂纹，并且这一裂纹的方向很可能不同于煤层的割理方向。覆盖层方位各向异性的存在显著地改变了CBS的剪切波特征，这需要分开讨论。

开始，计算了图2但具有弱各向异性覆盖层的Narellan模型的零偏移反射(图10)。覆盖层的γ参数是0.02，比煤层的γ参数(γ＝0.2)小10倍。虽然覆盖层各向异性很弱，但由于其厚度，积累的延迟时间(在煤层顶部)却很显著(即大约10毫秒)。

在这一例子中，覆盖层的裂纹方向与煤层的割理方向成30°。如图3所示，图10的合成数据是沿煤层各向异性主轴(x轴垂直于割理)采集。然而，在CBS时间窗中不匹配轨迹信号非常强，这在CBS上没有各向异性覆盖层的模型计算结果图3中是不存在的。明显这是前面分析方法的一个问题。

如果覆盖层是各向异性的，入射S波方解成与覆盖层裂纹(分别)成水平和垂直极化的快和慢模式。在最上层煤层顶部，两种模式又分解成与煤层裂纹(割理)平行和垂直极化的两个分量。这四种模式在CBS内回响，引起柱-腿复合，并且(在它们返回表面时)在覆盖层底部又分解。这形成了非常复杂的波场，这一波场不仅反映了煤层，也反映了覆盖层裂纹的影响。

如果煤层的裂纹与覆盖层相同，这一场面将大大减化。这种情况下，每个在覆盖层分离的S波在煤层中只产生一个S波(具有相同极化的S波)，并且覆盖层的影响被限制在一个附加的时间延迟内。然而，这附加的延迟引起不匹配轨迹的非零信号，即使CBS是各向异性的，这是因为两个波互不同相。因此，不分离出上层的各向异性上述简单的方法就不能应用。这一分析的最后结论是必须考虑覆盖层的方位各向异性的效应。下面的讨论将表明覆盖层和CBS裂纹方位不同的更一般的情况并不特别复杂。因此，这一层分离方法有广泛的应用，甚至超出CBM的范围，这将在下一段分开讨论。

总体上说，作为一种裂纹区集气点的诊断方法，这一模拟方法表明地震特征/参数R_a确定是煤层系各向异性好的整体量度。虽然这一特征对CBS几何构造的细节不灵敏，但却不仅依赖于煤层各向异性程度，而且还依赖于CBS内各向异性物质的裂纹。

各向异性方向沿垂直方向变化的模拟

Winterstein和Meadows(见Winterstein，D.F.和M.A.Meadows的“在Lost Hills矿场的剪切波极化和表面下应力方向，地球地理学，56，第1331-1348页，1991；以及美国专利5,060,203及5,060,204)曾提出一种层分离技术作为解决零偏移VSP中不同深度的裂纹方位的方法。他们认为他们的VSP技术同样可用于表面反射的范畴；但他们没有给出例子。他们解释为：对于表面反射数据信噪比太低；并且如果各向异性改变方向的地层与强反射层不精确重合，这一技术的假设就失败，且导致精度的损失。

虽然这两个观察结果是正确的，因为两个分离模式之间的时间延迟随传播路径逐渐发展(即，延迟并不在各向异性改变的地层立刻出现)，但(各向异性改变的地层和反射层之间)精确的重合并不需要。除此之外，最近关于二极剪切各向异性记录的经验(见M.C.Mueller的“Interpretation of Dipole Shear Anisotropy Log in a Thrust BeltSetting”，SEG扩展摘要，第＿＿页，1995)表明各向异性变化的地层用地层边界很难修正，即使是采用与煤层尺度相同或更低的分辨率的研究方法。这样的经验使人有信心认为不精确的重合通常不是一个问题。同样地，表面反射数据的信噪比也没发现有严重问题(参见T.A.Chaimov、G.J.Beaudoin、W.W.Haggard和M.C.Mueller的“剪切波各向异性和煤床沼气预测性”，SEG扩展摘要，第＿＿页，(1995)，以及其人中的现场数据例子)。

Winterstein和Meadows注意到Alford旋转要求两个不匹配轨迹(在噪声限制的范围内)是一致的，因此进一步指出垂直变化的各向异性方向在VSP范围内破坏了这一对称性。然而，在(具有真实的垂直传播目的)反射情况内，由于互逆原理(考虑到源与接收器位置的互换)，具有或不具有垂直变化的各向异性方向，不匹配轨迹都是一致的。在实践中这已被普遍观测到，即使是采用堆垒轨迹而不是直实垂直轨迹的情况也是如此。

回到本发明的方法，象Winterstel和Meadows一样，各向异性的主方向假设成是粗糙化层基深度的分段常数。对于各向异性的量设有如下的限制，该限制定义为(见Thomsen，L.的“方位角各向异性介质的反射地震学”，地球地理学，53，(3)，第304-313页，1988)：

γ (z) = \frac{{V_{| |}}^{2} - {V_{&perp;}}^{2}}{2 {V_{| |}}^{2}} = \frac{V_{| |} - V_{&perp;}}{V_{&perp;}} - - - - (5)

γ(z)的变化可以是缓慢的或剧烈的或随V₁(z)本身变化。这里，V₁是平行于“快”主方向极化的垂直传播剪切模式的速度，V_⊥是垂直于“快”主方向极化的垂直传播剪切模式的速度。应当知道，如果γ太大(＞10％)，并且岩石的弹性对称元素不包括一个水平对称平面，那么下述方法的应用是近似的，而不是精确的。

按照早期的工作，堆叠波被认为是垂直入射轨迹的精确替代，尽管这种情况有效的条件不十分清楚。下面的模拟结果是从垂直变化的HT1(水平横向各向同性)对称性介质中真实的垂直传播得到的。

解决方法是从反射剪切波数据中分离出上层各向异性的影响。我们将发现。处理覆盖层和cbs中裂纹方向不同的更一般的模型并不比上面所讨论的简单模型更困难。

多分量数据反射层分离

层分离的基本目标是将所记录的反射轨迹变换为相应于CBS上面的纯各向同性层的新的一组轨迹。即，各向异性效应只需从数据中分离出。从地震特征/参数Ra的高度讲，基本目标就是去除覆盖层的影响，获得表征CBS内各向异性特征的参数Ra(即方程4)。应当知道，本发明的反射层分离方法可以用于其它方面，而不能仅限于Ra的求解。它还可用于许多沉积层的方面，而不只是CBS模型；如上说过，煤层只是增强了覆盖层各向异性的效应。从某种意义上说，(即从反射的角度上说)它是对Winterstein和Meadows所提出的一种解决零偏VSP中随深度变化的裂纹方位问题即所谓的“极化层分离”技术的改进。假定认为各向异性主方向沿垂直方向变化是普遍的，就可推断出这项技术将可用于探测有裂纹区的集气点(即煤床的甲烷和密闭天燃气)。

这一方法是在源层的基础上完成的(即，按照Winsterstein和Meadows所做假定)，认为各向异性方向(虽然不需要其数值)在粗糙化的层内是分段常数。这一似乎真实的假设一般与各向异性最终的原因也即局域应力场有关。下面所述方法的成功肯定了这一假设在实际中的有效性。与Winsterstein和Meadows一样采用了相同的准则来确定这些粗糙化的层的顶层和底层，即层的厚度足以带来明显的(两种剪切模式之间的)时间延迟。在实践中，这并不是一个困难问题，虽然如前面Winsterstein和Meadows所建议的需对判据做出解释。

第一步是通过对图10中数据组做MS/MR旋转找出覆盖层的各向异性主方向。通常，四分量场数据的旋转每隔10°或15°完成一次，以便选择出使不匹配轨迹的能量最小化的角度。图11和12显示了原始数据相对于初始座标系分别旋转20°和30°的结果。从煤层最顶层反射即第一个和最突出的到达波)的不匹配分量在旋转30°后几乎消失(图12)。20°的旋转(图11)明显不充分。因此，30°的旋转将数据变换到覆盖层的主轴方向，从而确定了覆盖层的各向同性方向(当然，初始的采集可以是任何角度)。应当想起煤层最快的方向是南北向而覆盖层最快的方向是北偏东30°。

由于来自第一个煤层底部的反射的影响以及CBS顶部快和慢模式间的偶合效应，图12中第一个来波的不匹配分量并没完全消失，然而按照本方法覆盖层的各向异性主方向通过数据旋转很容易确定。图12中(慢)11轨迹和(快)22轨迹到达时间间的差异给出了覆盖层快和慢反射模式间的时间延迟(10毫秒)。这一差别用肉眼无法观测到但通过横向相关的方法很容易确定。图12中不匹配轨迹上CBS时间窗内回响/信号的存在是CBS和覆盖层两者裂纹方向不同引起的各向异性效应。为了确定它，必须消除覆盖层各向异性对图12中数据的影响。

层分离方法的目的是将从覆盖层各向异性主轴获得的反射轨迹(图12)变换成对应于CBS上面纯各向同性层的一组不同的数据。因此，在建立起最上一层粗糙层(具有均匀的各向异性方向，但任意的各向异性和阻抗变换)的底限后，下一步骤是通过将慢模式11踪迹上移10毫移来对从底部(也即CBS顶部)反射的慢模式11轨迹和快模式22轨迹做时间补偿。然而，不匹配轨迹12和21只上移该量值的一半(5毫秒)，因为这两种轨迹传播的一个路程(例如，向下)只作为慢模式，而另一路程是快模式。快模式轨迹22一点也没平移；最后轨迹结果如图13所示。这一方法消除了由覆盖层各向异性引起的时间延迟。覆盖层中两种剪切模式间的衰减和几何发散的差异假设小得可以忽略。这一假设在层分离异法的现场应用中已得到证明(后面将对一个例子做详细说明)。

熟悉这一技术领域的人将会意识到对不匹配轨迹的特殊处理只是本反射层分离方法区别于Winsterstein和Meadows的VSP层分离方法的一个方面(见Winterstein，D.F.和M.A.Meadows的“在Lost Hills矿场的剪切波极化和表面下应力方向”，地球地理学，56，第1331-1348页，1991；以及美国专利5,060,203和5,060,204)。按Winsterstein和Meadows的方法对三个延迟轨迹做同样的处理不会给出正确的结果。

做了这样的修正以后，可以完成下一步3。图13的数据组完全可以按照仿佛它们是从一纯各向同性的覆盖层采集而来的方式处理。通过对图13的轨迹做张量旋转(这种情况是30°旋转)，就找到了煤层的各向异性主轴(见图14)。不匹配分量的零能量证明层分离方法是正确的。不匹配轨迹的零值也告诉我们在CBS顶层的反射/传播过程中，两种剪切模式的耦合是非常弱的。

图14中的轨迹作为图10中的原始数对应于同样的座标系；两种图间的差别完全是由于为使轨迹从煤层各向异性角度更易解释的层分离方法造成的。当然，除这一教学方法以外，关于图10的方法没有什么特别；尤其是，本方法不需要以前的关于各向异性，煤层几何构造等的知识。

最后，进行一个45°的正则的旋转以获得参数Ra，即CBS各向异性的总体量度。将图15与图14做一比较是有趣的，图4是对同样模型做的同样的采集，但却具有各向同性的覆盖层。这些图(和Ra的值)实际是一致的，这一进一步证实了这一层分离方法的可靠性。

应当强调本发明的极化层分离方法可以重复多次，只有手头拥有具有几个裂纹方向不同的层的介质资料。

然而，层分离方法的多次应用可能会导致误差随深度的积累，因此这一方法应区别应用。

详细实施例

上述煤层中裂纹探测的剪力波方法曾在新南威尔士的Sydney盆地做过现场测试。采用的是通常的多源/多接收器(MS/MR)采集技术；但是处理方法和各向异性的评估设计成适合于由薄的多层煤层系引起的特殊情况。

在Sydney盆地感兴趣的区域，上部煤层(Bulli煤层)出现在400到900的深度范围(见图12所示的一个典型煤层系的代表)。近表面由三叠纪时期的Hawkesbury砂石构成(即这些单元的外露部分形成了悉尼壮观的海岸礁石。表面状况为土壤层厚度变化的起伏的牧场和乡村丛林。由于目标相对浅，并且CBS裂纹群实质上相当小(即，100米的量级或更小)，减小了采集深度点的间距。勘察区域大到足以探测到可能的集气点，而不管可能的变化的数据质量。

记录设备包括两个Geosource Dss-10 24通道系统，3-分量地音探测器和2 Hertz型M-13水平振动器。记录电子仪器是带有地力相位和地力控制的Pelton Model5 Advanc I。振动器产生50,000磅的力，其地力信号与参考信号同相位。拥有的记录能力(即48个通道)，均分给三分量地音探测器，也即对探测器的每次使用只有16个通道/模式。因此在采集计划的安排上充分利用与地音探测器相对来说更扣动的振动器。也即，对地音探测器的每次使用，振动点都要占用多次。计划扫描长度为16秒范围为8到80Hz，记录成20秒的记录文件。横相相关性由当地的Vendor完成。剪切模式(即顺线/横线)分离和记录设备是常规的，数据通过一个ufh程序来整理(即多分量数据一边进行顺线，横线记录的多路传输，一边进行多分量数据的整理)，以后采用PROMAX来完成这项工作。

由于使用的是有限的列(即接收器和源)，表面波主导着原始现场普通探测采集。看到由顺限源相对于横线源的运动激发的不同表面模式被记录进顺线、横线和垂直地音探测器是可能的。例如，纵向波的第一个折线很容易从顺线源对垂直接收器的记录认出，然而各种剪切和表面模式常相互关联，并且影响振个记录，使所有反射都模糊不清。各种波的速度为：

直接纵波 400m/sec

直接剪切波 1900m/sec

表面模式 750m/sec到1500m/sec

目标剪切堆叠速度 2000m/sec

目标间隔Vs 1000m/sec

目标间隔Vp 2500m/sec

覆盖层间隔Vs 2100m/sec

覆盖层间隔Vp 3600m/sec

多分量方法一个可能发生混乱的方面是模式的分离和组合。每个现场记录相应于一特定的源模式，顺线或模线。每个多分量记录的轨迹将被多通道传输：顺线，横线和垂直。一旦完成正确的分离，数据就被汇集成3-，4-，6-或9-分量的图线。这些记录可推导出一个6分量的数据组(即，横线振动器和顺线振动器对顺线接收器，横线接收器和垂直接收器)。前面所述分析需四个水平分量。一旦这些分量被分离处理它们是简单的，因为所有分量都是同一处理流程。

处理需要折射静力学(即，从横线/模线分量计算)，然后，通常需要一个简单的放大或者减噪技术来保持(各不同分量)相对的幅值。这意味着，如果例如需要自动增益(AGC)定标量必须计算一次并储存起来对所有分量级统一的应用。这一计算方法可以通过堆积由每个单个分量派生出的标定量文件来提供/给“多分量”。类似地，减噪技术不应当涉及难以追踪的放大技术。扩散速度计算应当进行一次，然后一致应用于所有四个分量。计算以横线/横线分量做为典型进行完成，虽然任何具有好的(预-旋转)反射层强度的分量都可用来计算。通常，预旋转剪切波数据质量并不适合于外观速度的分析。在这种情况下，常速度堆迭是非常有用的。

各向异性评估可通过预旋转或后旋转完成。但是通常不鼓励采用预旋转，因为可得到的旋转算法是用于正入射路径。标准的首先通过的各向异性分析是由强堆叠完成。四个分离的堆积可认为是一个2×2的矩阵。旋转分析通过保持记录为这一矩阵格式，并对每一旋转产生一个分离的矩阵来大大地简化。

首先通过的旋转的标准角度为0°，15°30°，45°，60°，75°和90°(对这方面知识熟悉的人会知道90°旋转只是转换了主元和非对角元的位置，即矩阵变换操作)。附录B给出的是完成这一多次旋转和适当显示工作的一个C语言程序。旋转完以后，顺线和横线术语不再合适。需要用相对于罗盘方向的粒子运动来规定分量。在完成了一系列的旋转后，就容易通过研究非对角分量的反射器强度来观察到最佳的旋转。非对角元素反射器强度最小化的地方，就达到了最佳的旋转(即，图16B)。在这首次通过的角度附近用较小的角度差，重复这一步骤，就可保证得到最佳角。附录C是四分量(4C)向反射和VSP剪切波数据层分离的数学描述。

非对角区上相关能量的存在表明了方位各向异性。通过旋转使这一能量较小化，就修正了勘察(线)方位与覆盖层主轴之间的不同轴性。由于在覆盖层时间窗内任何时间非对角区都不存在有相于能量，可以给出这一最后结论。

按照本发明，在补偿了覆盖层各向异性后，将最优旋转矩阵完成一45°的旋转，展现出CBS时间窗内非对角区的相干能量。这指出和量度了覆盖层下CBS内的各向异性区域。由于指明了高裂纹渗透性，其方向是沿着CBS内各向异性的快方向，这些区域被解释为高裂纹强度的集气点的指示(见T.A.Chaimov、G.J.Beaudoin、W.W.Haggard和M.C.Mueller的“剪切波各向异性和煤床沼气预测性”，SEG扩展摘要，第＿＿页，1995)。根据本发明，由于在此没有看到显著的相关能量(见图16C)，可以得出结论，有很少或没有CBS方位各向异性，因此在此地点有很少或没有CBS裂纹。这一试验结论由下面得到证实：

^*芯棒分析(芯中无裂纹)，以及

^*VSP分析(即，采用4分量VSP数据组的下孔旋转)。

对于那些熟悉此行的人，从前面的描述，可以看到，众多的变化，选择和修改是显而易见的。因此，本描述只能解释为说明，其目的是传授那些精通此行的人完成本发明的方法。可以进行各种改变，材料可代替，本发明的特点可以利用。例如，虽然本发明是从煤层的范围来描述的，但它的方法可以推广到其它碳质系例如antrium和泥盆纪页岩。除此这外，在覆盖层实质上是方位各向同性的情况，不必进行层分离。同样，本发明的层分离方法不限制在碳质层素的范围内，它可用于具有垂直变化的各向异性的地表下几何构造的广阔范围内。最后，应当认识到如果地面下构造事实上是由分段常数的多向异性方向来描述，那么在每个连续的粗糙化层内，数据矩阵无论是从VSP的角度还是反射的角度讲都是对称的。根据试验构造按照上述方法或者Winsterstein和Meadows的方法，就可成功地分离出各个层的各向异性。因此，应当承认在不脱离所附权利要求所界定的本发明的精神实质和范畴的情况下可以做出各种修改、选择、变化等。当然。我们要把所有这些在权利要求范围内的修正包容在所附的权利要求中。

附录A：各向异性与割理

煤床中的割理通常为张拉裂纹(即，位移为割理平面法向，而不是象剪切裂纹在平面内方向)。通常有两组割理，均为垂直向(或为煤层的法向，如果煤层是倾斜的)，两组割理互折正交。这一正交特性并非偶然，而是应力张量正交性的结果。从宏观的角度，这是一基本物理规律的表现。

尤其是，割理方位受局部水平应力场控制。如一些人所指出的，如果割理是由于煤(在煤化过程中)处在垂直压缩场中产生体积收缩而引起，那么割理方位将是多边形的方位，正如泥土裂纹或玄武岩中柱状连接一样。这样的体积收缩无疑存在，但割理方位由水平应力场确定，并且这一方位在与一个地震Fresnel区相比要大的水平距离上是稳定的。

这两个正交的裂纹组通常是不等同的。一组(面割理)较长，另一组(碎割理)通常在碰到前一组割理时终断开。液压传导性明显依赖于割理的长度和宽度两者，而地震各向异性只依赖于割理的长度(见后面的讨论，这一特性导出了方程Sa)。因此，用单组裂纹(即，面裂纹)来近似地震各向异性是合理的。

文献中通常讨论的单组环形裂纹引起的各向异性是由Hudson给出(见Hudson，J.A.的“断裂固体的整体特性”，Math Proc.Camb.Phil.Soc.，88，第371页，1980)。Thomsen(见Thomsen，L.的“多孔岩石中因对齐的裂纹而产生的弹性各向异性”，欧洲协会勘探地球地理学公约摘要，53，第244-245页，1991；Geoph.Prospctg.，1995)表明，如果裂纹发生在多孔岩(相对于土壤)，那么多孔性对各向异性起作用，即使(按定义)对各向异性没有优先的方向性。这一现象的发表是由于在地震频率点，裂缝中流体与孔中流体处于压力平衡。这一判断为Rathore等的实验所证实(见Rathore，J.S.、Fjaer，E.、Holct，R.M.、Renlie，L.的“具有受控制的裂纹几何图形的人工合成沙岩的声学各向异性”，Geoph.Prospctng.，第＿-＿页，1995)。

在煤层中，煤基质的渗透性太慢使这种效应可以忽略。方程减少为Hudson方程(见Hudson，J.A.的“断裂固体的整体特性”，MathProc.Camb.Phil.Soc.88，第371页，1980)的表示法，这些方程为V_p(θ)＝V_p(90°)[1+δcos²θsin²θ+εcos⁴θ] (1a)

这里，V_p(90°)和V_s(90°)是波沿对称方向传播的P和S速度(即，由于裂纹是垂直的，水平波方向与垂直方向成90°)；θ是波的入射角，以垂直于裂纹的平面的垂直方向开始量度。各向异性参数，Thomsen是以弹性模量给出，还可通过下面方程以裂纹密度给出(参见Thomsen L.的“多孔岩石中因对齐的裂纹所产生的弹性各向异性”，欧洲协会勘探地球地理学公约摘要，53，第244-245页，1991；Geoph.Prospctg.，第＿-＿页，1995)方程为：

ε＝0 (2a)

δ = - \frac{16}{3} (\frac{1 - 2 v}{2 - v}) η - - - - (2 b)

γ = \frac{8}{3} (\frac{1 - v}{2 - v}) η - - - - (2 c)

其中V是煤基质的泊松比，无量纲裂纹密度是

η = N_{v} < \frac{d^{3}}{8} > - - - - (3 a)

其中N_v是单位体积的裂纹数目，d是裂纹直径。括号表示体积平均。

对于垂直入射速度：V_p(0°)＝V_p(90°) (4a)

方程(4b)和(4c)与方程3明显一致。

在这种特殊情况下(即垂直传播无孔)，可进一步推导出关于煤层范围内的更满意的理论。考虑到裂纹对着1-方向(即x-方向)，并将相关的(面)裂纹密度记作η₁。下一步，作为替代，考虑沿2-方向的一组(碎)裂纹，裂纹密度为η₂。最后，考虑这两组裂纹的结合；在线性低裂纹密度范围内，很明显裂纹效应可假设为线性，因此剪切各向异性为：

对于煤，泊松比v＝0.33似乎是可取的(见Greenhalgh，S.A.和D.W.Emerson的“来自新威尔士的悉尼盆地的煤测量岩石的弹性特性”，勘探地球地理学，17，第157-163页，1986)，所以方程(5a)的系数是16/3。

从方程(5a)可清楚看到裂纹的大小强烈地影响着各向异性。如果碎割理为面割理的一半长，并且数量一样多，那么碎裂纹密度η₂为面裂纹密度η₁值的八分之一。两个正交的裂纹组通常不等价。一组(面割理)较长；另一组(碎割理)在第一组处终结。液体传导性明显依赖于割理的长度和宽度二者，而地震各向异性只依赖于割理的长度。因此，用单组裂纹(即面割理)来近似地震的各向异性是合理的。因此，至少在现在在方程(5a)以及以后的讨论中忽略碎割理，相信是合理的。

附录B：旋转程序

下面是一C语言程序，该程序可允许用户采用四数据矩阵文件(以一个“根名”后跟.11，.12，.21，.22扩展名命名)并且可使其按要求的“起始角度”：“终结角度”和“增量角度”来旋转。该程序然后将四输入文件转换成一输出文件中的单个记录，如图16a所示，其中完成的每一旋转作为含有4-分量数据矩阵的分开的记录附在输出文件上。

#！ /bin/csh
# this is a C-shell script for ITERATIVE CONSTANT ANGLE
# INCREMENT ROTATIONS
#
#   data set name(root) NM
#   start angle，       sa
#   end angle，         ea
#   angle increment，   da
#
# the″@″sign assigns a numeric value to the name
# immediately following it and the number in ″[]″denotes
# the position of the nth argument on the# command line
unset noclobber filec
if(S#argv＜3)then
   echo″ ″
   echo″ ″
   echo″ ″
   echo″rotzs2.x：Iterative constant angle rotations″
   echo″ ″
   echo″″
   echo″Input：″
   echo″″
   echo″  File Name(root)    -N[NM]″
   echo″  Start angle        -sa[sa]″
   echo″  End angle          -ea[ea]″
   echo″  Angle Increment    -da[da]″
   echo″″
   echo″Usage：″
   echo″″
   echo″rotzs2.x -N[]-sa[]-ea[]-da[]″
   exit
endif
set NM＝′usparg-N″ ″S*′
@ sa＝′usparg -sa″O″S*′
@ ea＝′usparg -ea″O″S*′
@ da＝′usparg -da″O″S*′
				
				<dp n="d27"/>
@ nsamp＝′getval-NSNM.11-LNumSmp′
@ nsi＝′getval-NSNM.11-LSmpInt′
@ nrecs＝′getval-NSNM.11-LNumRec′
@ tnrecs＝2 * Snrecs
@ te＝Snsamp * Snsi
@ dte＝2 * Ste
@ mte＝-1 * Ste
if(Ssa＝＝Sea)then

    echo″You must provide reasonable rotation parameters(-sa[]，-ea[]，-da[])″

    exit 1
endif
echo″Data Set Name(root)＝″SNM
echo″Start Angle＝        ″Ssa
echo″End Angle＝          ″Sea
echo″Angle Increment＝    ″Sda
# initialize the internal parameter vv to start at sa
@ vv＝Ssa
# initialize a iteration counter
@ numiter＝1
# initialize a temp parm
@ temp＝-1
@ vw＝0
# loop for as long as vv is less than or equal to ea
###################################   start while loop
#################################
while(Svv＜＝Sea)
# print to stdout then current value of vv

        echo″starting on angle＝″Svv

        echo″at″′date′

        echo″this is iteration number＝″Snumiter
# for the angle vv do const.angle rotations(on cdp stacked
# data)

    rotzs2-NSNM-RSvv
# each const.angle rotation is a record appended to the
# final output file for each component 11，12，21，22
# Assemble a rotation matrix for each angle output from
# rotzs2.
#                 11，12    
#                 21，22
# give each component twice the input samples(time)for
# vertical stacking
if(Svv＞-1 &amp;&amp; Svv＜10)then

     wind-eSdte-NSNM.R″+″OSvv.11-OSNM.11.wind

     wind-eSdte-NSNM.R″+″OSvv.12-OSNM.12.wind

     wind-sSmte-NSNM.R″+″OSvv.21-OSNM.21.wind

     wind-sSmte-NSNM.R″+″OSvv.22-OSNM.22.wind

    /usr/bin/rm SNM.R″+″OSvv.11
				
				<dp n="d28"/>
       /usr/bin/rm SNM.R″+″OSvv.12

       /usr/bin/rm SNM.R″+″OSvv.21

       /usr/bin/rm SNM.R″+″OSvv.22
endif
if(Svv＞-10 &amp;&amp; Svv＜-1)then

     @vw＝Svv * Stemp

     wind-eSdte-NSNM.R″-″OSvw.11-OSNM.11.wind

　　 wind-eSdte-NSNM.R″-″OSvw.12-OSNM.12.wind

　　 wind-sSmte-NSNM.R″-″OSvw.21-OSNM.21.wind

　　 wind-sSmte-NSNM.R″-″OSvw.22-OSNM.22.wind

　　/usr/bin/rm SNM.R″-″OSvw.11

　　/usr/bin/rm SNM.R″-″OSvw.12

　　/usr/bin/rm SNM.R″-″OSvw.21

　　/usr/bin/rm SNM.R″-″OSvw.22
endif
if(Svv＜-9)then

　　 wind-eSdte-NSNM.RSvv.11-OSNM.11.wind    

　　 wind-eSdte-NSNM.RSvv.12-OSNM.12.wind

　　 wind-sSmte-NSNM.RSvv.21-OSNM.21.wind

　　 wind-sSmte-NSNM.RSvv.22-OSNM.22.wind

　　/usr/bin/rm SNM.RSvv.11

　　/usr/bin/rm SNM.RSvv.12

　　/usr/bin/rm SNM.RSvv.21

　　/usr/bin/rm SNM.RSvv.22
endif
if(Svv＞9)then

　　 wind-eSdte-NSNM.R″+″Svv.11-OSNM.11.wind

　　 wind-eSdte-NSNM.R″+″Svv.12-OSNM.12.wind

　　 wind-sSmte-NSNM.R″+″Svv.21-OSNM.21.wind

　　 wind-sSmte-NSNM.R″+″Svv.22-OSNM.22.wind

　　/usr/bin/rm SNM.R″+″Svv.11

　　/usr/bin/rm SNM.R″+″Svv.12

　　/usr/bin/rm SNM.R″+″Svv.21

　　/usr/bin/rm SNM.R″+″Svv.22
endif
# vertically stack the 11 on top of the 21 and the 12 ontop
# of the 22

　　 vstak-N1SNM.11.wind-N2SNM.21.wind-OSNM.col1

　　 vstak-N1SNM.12.wind-N2SNM.22.wind-OSNM.col2

　　/usr/bin/rm SNM.11.wind

　　/usr/bin/rm SNM.12.wind

　　/usr/bin/rm SNM.21.wind

　　/usr/bin/rm SNM.22.wind
# append col2 next to(after)col1

　　append-NSNM.col2-OSNM.col1
				
				<dp n="d29"/>
# get record numbering and traces cleaned up before assembly
# in output file
# note that when modifying trace headers you must use a
# different-O from-N

　　 utop-R1-LStnrecs-kORecNum＝Snumiter-NSNM.col1
-OSNM.temp
if(Snumiter＝＝1)mv SNM.temp SNM.matrix
if(Snumiter＞1)append-KSNM.temp-OSNM.matrix
# now increase the current(numeric)value of vv by the
# value da

　　    @vv＝Svv+Sda
# increase the iteration counter by 1

　　    @numiter＝Snumiter+1

　　    rmprint
end
##################################    end ″while″loop
################################
# clean up your mess buster...
/usr/bin/rm SNM.col1
/usr/bin/rm SNM.col2
/usr/bin/rm SNM.temp
exit 0

附录C

垂直变化的方位各向异性的剪切波数据

粗糙层分离

本发明的原理表示在图A和B中。

图A大略地表示出地表的两个层10和20，每层中都散布着一组平行排列的垂直裂纹12和14(该裂纹形成了每层的方位各向异性)，但每层的方位各向异性又不同。箭头表示的是表面18上产生的固定极化方向的剪切波的传播路径。当波向下传播时，波在层10中分成两个波16a和16b，然后在第二层中分成4个波，在通过第一层反射回来时又分成6个波，等等。记录的信号(含有所有这些到达波)是非常复杂的，没有处理是很难解释的。

图B大略表示了本发明是如何处理图A的数据使数据仿佛是出自图B简单的构造。本发明以上层10中“分离”出各向异性，使其成为有效的各向同性，并且波16在该层中向下传播或向上返回过程中没有分裂。这将问题减化为单层中波的分裂分析。反射层分离(4C)

4分量(4C)反射地震学装置或试验的接收信号(即源的两横相方位，接收器的两个横向方位)可写为：S₀(2t₁+2t₂)＝ (方程1)2R(-θ₁)P₁_T²¹R(θ₁-θ₂)P₂_RP₂_R(θ₂-θ₁)T¹² P₁_R(θ₁)s₀w(t)S₀

这里，除非另外注明，黑体字的项都为矩阵。

S₀(2t₁+2t₂)是一个2×2的地震数据轨迹矩阵，每一列代表由相应源分量激起的接收器的两个分量。按照通常的实际作法，通过对向外移进行修正，每一个这种轨迹由一个代替的法向入射轨迹构成，这一入射轨迹通过几个模糊的轨迹堆叠而成。下标“o”表示该向量是相对于勘察座标系。自变量(2t₁+2t₂)表示通过近似层1和层2穿过两个路径(即向下和向上回来)的波，如图A和图B所示。

目标是求解方程(1)中的“主时间序”(即，解出穿过层系，相应于纯模式，快或慢的束波)。从这些束波中，就推导出各个层各向异性的方位和值。变量注释顺序从右边开始：

没有自变量的S₀是源的方位矩阵。假设一个顺线源和一个强度相等的横线，S₀可记作

S_{0} = [\begin{matrix} 1 & 0 \\ 0 & 1 \end{matrix}] = 1;

w(t)是地震扰波。

S₀是初始幅值；

R(θ)是旋转矩阵，使向下传播波主坐标系中的矩阵调整到层1中更优化的方向(即在图A和B的例子中，层1的裂纹走向)，并且

P₁是顶层的传播矩阵，

“_”指卷积；

A₁ ^f是考虑了几何传播、衰减、发散等用于快分量的传播“滤波器”；

t₁ ^f＝z₁/v₁ ^f是单向路程的传播时间；

A₁ ^s和t₁ ^s是慢分量的对应量。

T¹²是层1和层2间的比例传递系数。

R(θ₂-θ₁)进一步旋转到层2的角度。

P₂是层2的传播算子为：

R是层2的底层反射系数矩阵(假设为对角矩阵)

注意以下的旋转矩阵(例如R(θ₂-θ₁)，R(-θ₁)，等)的自变为前面旋转自变量的负号，表示相反旋转的意思。因子“2”是考虑到与自由表而相互作用。

通过对S₀(2t₁+2t₂)左乘R(θ₁)右乘R(-θ₁)，数据可旋转到与层1各向异性主方向相重合的座标系。方程(1)左侧就变成以层1座标系表示的4C数据：

S₁(2t₁+2t₂)

＝R(θ₁)S₀(2t₁+2t₂)R(-θ₁) (方程6)

＝P₁R(θ₁-θ₂)P₂_2T²¹Rs₀w(t)T¹²_P₂_R(θ₂-θ₁)P₁ (方程7)

在两个角度相同的特殊情况，表达式减化成各向异性方位均匀的情况，即，Alford(见Alford，R.M.的“方位角各向异性存在时的剪切数据”，第56届Ann.Internat.Mtg.Soc.Exploc.Geophys.，扩展摘要，476-479，(1986))以及Thomsen(见Thomsen，L.的“方位角各向异性介质的反射地震学”，地球地理学，53(3)，304-313，(1988))所考虑的情况，方程(7)减化到他们的方程。在小于t₁ ^f的时间内，初始未知θ₁角的选择应使方程(2)旋转后的数据的非对角元极小。在大于t₁ ^f的时间上非对角元上大的能量可以用来识别t₁ ^f。

定义模式之间在层1底部的时间延迟为

/Δt₁≡t₁ ^s-t^f。

层1各向异性的分离可利用一模式促进算子完成：

方程(7)左右乘D₁并假设传播算子A₁ ^f和A₁ ^s相同，那么：

因此D₁的应用，促进了慢模式(即分离了层1的各向异性)分离层方程(7)左边项为：

S₁(2t₁+2t₂)≡D₁_S₁(2t₁+2t₂)_D₁ (方程10)S₁(2t₁+2t₂)＝ (方程11)

A_{1} &CircleTimes; δ (t - t_{1}^{'}) R (θ_{1} - θ_{2}) &CircleTimes; P_{2} &CircleTimes; {2 T}^{21} R s_{o} w (t) T^{12} &CircleTimes; P_{2} &CircleTimes; R (θ_{2} - θ_{1}) δ (t - t_{1}^{'}) &CircleTimes; A_{1}

方程(12)方括号中的量是“层分离的主时间序”：

注意22分量较大的延迟。

按照前述方法，层2各向异性的方位和值可以确定。解方程(14)，S(t)＝R(θ₂-θ₁)S₁(2t₁+2t₂)R(θ₁-θ₂) (方程15)

上式完全类似于方程(6)。从表面看，分离层主时间序只是分离层再旋转的数据。

层分离方法由下式给出(参见方程(10))：S₁(2t₁+2t₂)＝

按方程(6)，上面矩阵元素中的上标“1”表示该数据已旋转到层1座标系。自变量表明非对角线轨迹作了向后的单程模式延迟平移，而22踪迹作了向后的双程延迟平移(即两倍)。后面的层的分离可以用此方法推广完成。VSP层分离(4C)

从VSP角度，那些熟悉此行的人会知道来自4-分量(4C)VSP装置或试验(即接收器与源座标系重合)的接收信号可写为：S₀(t₁+t₂)＝R(-θ₂)P₂_R(θ₂-θ₁)T¹²P₁_R(θ₁)S₀W(t)S₀ (方程18)假设顺线和横线源相等，并旋转到层-1的主轴，上式成为：S₁(t₁+t₂)＝ (方程19)R(θ₁-θ₂)P₂_S₀W(t)T¹²R(θ₂-θ₁)_P₁

仅采用一个模式促进算子D₁(因为仅向下传播)进行层分离，方程

(19)成为：S₁(t₁+t₂)≡S₁(t₁+t₂)_D₁ (方程20)

方括号中的量是分离层主时间序列：

象前面一样通过解方程(2 1)可以得到：S(t)＝R(θ₂-θ₁)S₁(t₁+t₂)R(θ₁-θ₂) (方程23)

表面看，分离层主时间-序列只是分离层再旋转数据。因此，层分离方法为：

按照方程(6)在这些矩阵元素中，下标“1”表示该数据已旋转到层-1座标系。自变量表明只有第二列的(旋转的)数据作了时间平移。这一方法不同于Winterstein和Meadow的方法(见Winterstein，D.F.和M.A.Meadows的“Los Hills矿场中的剪切波极化和表面下应力方向”，地球地理学，56(9)，1331-1348，以及“在Cymric和Aairoad空隙油田具有深度的热波极化方位角的变化”，地球地理学，56(9)，1331-1348(1991)，在其中对其它非对角轨迹也作了时间平移。他们给出了启发性的自变量，而没有推导和例子。

Claims

1.一种评估地表随深度变化的方位各向异性属的表面反射数据的方法包括如下步骤：

I)采用多源和多接收器，单源和多接收器以多源和单接收器中之一种方法，采集一剪切波数据堆积。

II)将每个上述剪切波数据堆积的轨迹旋转一个相应于上述数据堆积方位与最浅层方位各向异性方向之差的角度。

III)确定上述层中各向异性方向垂直变化的方位。

IV)量测上述层方位各向异性的程度。

V)完成上述数据堆积的层分离以获得所谓的主时间序列，上述层分离包括上述数据堆积的慢极化轨迹和上述方位各向异性数据堆积的混合极化轨迹的修正。

VI)将步骤(V)形成的数据组旋转一个相应于下面层的方位各向异性方向与上面层方位各向异性方向之的角度以识别上述下面层的方位各向异性主方向。

VII)对上述最浅层下面的至少一个层重复步骤(III)到(VI)的步骤。

2.如权利要求1所述的方法，其特征在于步骤(VI)的数据组是四分量数据组，包括：一个顺线源和顺线接收器分量；一个顺线源和横线接收器分量；一个横线源和顺线接收器分量；一个横线源和横线源接收顺分量。

3.如权利要求1所述方法，其特征在于步骤IV)包括上述层慢方向极化轨迹相对于快方向极化轨迹总延迟Δt的确定步骤。

4.如权利要求3所述方法，其特征在于包括步骤：将上述慢方向极化的上述数据的轨迹在时间上向后平移上述延迟Δt；并且将混合方向极化的轨迹平移上述延迟的一半。

5.如权利要求4所述方法，其特征在于步骤(V)还包括对小于t_bot的时间，将上述数据组轨迹静噪的步骤，这里t_bot对应于上述层的底部。

6.一种从具有一最浅层的方位各向异性沿层深变化的地质构造中采集的表面反射剪切波轨迹的处理方法包括步骤：

I)识别快主时间序列，慢主时间序列和最浅层方位各向异性方向；

II)识别层的底部；

III)采用上述步骤(I)快主时间序列和慢主时间序列来确定层的方位各向异性的程度。

IV)采用层的地震分量并形成一多分量数据组；

V)由于层的方位各向异性的效应修正上述数据组的慢方向极化轨迹。

VI)考虑层的方位各向异性的效应修正上述数据组的混合方向极化轨迹。

VII)对于时间小于层底的时间，对上述数据组的原始轨迹静噪。

VHI)将步骤(VII)形成的数据组旋转一个对应于上述后一层的方位各向异性方向与最浅层方位各向异性方向之差的角度以识别上述后一层中的主裂纹系。

IX)对后面的另一层重复步骤(I)到(VIII)。

7.如权利要求6所述的方法，其特征在于对每层的处理步骤(III)包括确定上述层的慢主时间序列相对于快主时间序列的总的延迟Δt的步骤。

8.如权利要求7所述的方法，其特征在于对每层的处理步骤(V)包括将上述数据组的所述慢方向极化的轨迹在时间上向后平移上述延迟Δt的步骤；并且步骤(VI)包括将混合方向极化的轨迹平移上述延迟的一半的步骤。

9.如权利要求6所述的方法，其特征在于步骤(IV)所述数据组是一四分量数据组，该数据组包括：一个顺线源和顺线接收器的分量，一个顺线源与一个横线接收器的分量；一个横线源和一个顺线接收器的分量；一个横线源与一个横线接收器的分量。

10.一种评估方位各向异性随深度变化的地质层的表面反射数据的方法，包括下列步骤：

I)沿勘察线采集一来自地质层的剪切波堆积，上述数据堆积包括多项地震轨迹；

II)利用步骤(I)的数据堆积识别所述勘察线与第一层主裂纹系方位间的夹角以及所述第一层的底部。

III)比较上述主裂纹系快地震轨迹与慢地震轨迹以获得慢方向极化轨迹相对于快方向极化轨迹的总延迟Δt的量度。

IV)采用上述第一层的地震分量形成一多分量数据组；

V)将上述数据组慢方向极化轨迹在时间上向后平移上述总时间延迟Δt；

VI)将上述上述数据组不谐调方向极化的轨迹在时间上向后平移上述总延迟Δt的一半。

VII)将步骤(VI)的形成的上述多分量数据组旋转一相对于所述下一层方位各向异性方向与上述第一层方位各向异性方向之差的角度以认别上述下一层方位各向异性的方位。

VIII)对上述第一层下面的一层重复步骤(III)到(VII)的操作。

11.如权利要求10所述的方法，其特征在于步骤(I)的数据堆积是从多源和多接收器，单源和多接收器以及多源和单接收器中的一种采集得来；其中步骤(IV)的多分量数据组包括：一个顺线源和顺线接收器分量，一个顺线源和横线接收器分量，一个横线源和顺线接收器分量和一个横线源和横线接收器分量。

12.如权利要求11所述的方法，其特征在于步骤(II)所述主裂纹系由一快主时间序列或一慢主时间序列来表征；每个时间序列是通过将上述剪切波数据堆积轨迹旋转一相对于上述堆积的方向与上述层的方位各向异性方向之差的角度来获得。

13.如权利要求10所述的方法，其特征在于步骤(VII)包括将时间t_bot小于上述层底部时间的上述数据组的原始轨迹静噪的步骤。

14.如权利要求10所述的方法，其特征还在于包括下述步骤：在上述层形成的窗上获得快极化轨迹的包迹的幅值A。

15.如权利要求14所述的方法，特征在于还包括下述步骤：将步骤(VIII)形成的数据组旋转到偏离主方向45°的方向上；在上述窗上，获得至少一个旋转数据组的/不匹配极化轨迹的包迹幅值A₄₅的量度。

16.如权利要求15所述的方法，还包括形成混合/不匹配极化轨迹的包迹幅值A₄₅的量度与快极化轨迹的包迹的幅值A₀的量度的比值以识别裂纹的方位的步骤。

17.如权利要求15所述方法，还包括步骤：用快极化轨迹的幅值A₀的量度对至少一个混合/不匹配极化轨迹的幅值A₄₅的量度进行正则化以识别裂纹方位。

18.如权利要求14所述的方法，其特征在于快极化轨迹的幅值A₀的量度代表其在上述下一层形成的窗上的平均值。

19.一种分析地震表面反射剪切波数据以评估地表下裂纹的方法包括下列步骤：

(a)采集一包含多个地震轨迹的表面反射剪切波数据。

(b)通过对与慢剪切波极化方向重合的数据分量做一静平移以及对既不与慢剪切波极化方向重合又不与快剪切波极化方向重合的数据分量做上述静平移量一半的平移来分离指明了极化变化的最浅层。

(c)将最浅层下面的主时间序列数据单元旋转一相应于最浅层的方位各向异性方向与下一层方位各向异性方向之差的角度。

20.一种处理相对薄，随深度变化的方位各向异性层上的表面反射数据的地震方法包括步骤：

I)从多源和多接收器，单源和多接收器以及多源和单接收器中之一种采集一剪切波数据堆积，上述剪切波数据堆积，特征在于至少包含一种四分量数据组，该数据组包括一顺线源分量，一顺线接收器分量，一横线源分量，一横线接收器分量。

II)将上述剪切波数据堆积的每个分量旋转一相应于上述数据堆积方位与最浅层的方位各向异性方向之差的角度。

III)确定在上述最浅层中各向异性方向垂直变化的方位。

IV)通过确定上述层慢方向极化轨迹相对于快方向极化轨迹的总时间延迟Δt来量测上述最浅层的方位各向异性的程度。

V)完成上述数据堆积的层分离以获得最浅层的主时间序列，上述考虑方位各向异性的层分离方法包括通过将上述数据组慢方向极化轨迹在时间上向后平移上述延迟Δt，以及将混合方向极化轨迹平移上述延迟一半来修正上述数据堆积的慢极化轨迹和上述数据混合极化轨迹的步骤。

VI)对最浅层下面的层重复步骤(III)到(V)。

VII)将步骤(VI)形成的数据旋转一相应于上述下层的方位各向异性方向与上述最浅层方位各向异性方向之差的角度以识别上述下层的方位各向异性的主方向。

21.一种分析表面反射地震剪切波数据以评估剪切波极化随深度的变化的方法，其特征在于包括步骤：

a)从指明了极化变化的层的数据中进行分离，其步骤为：

I)识别每层与慢剪切波极化方向把一致的数据分量，与快剪切波极化方向相一致的数据分量以及与上述快和慢剪切波极化方向均不一致的数据分量；

H)确定快剪切波和慢剪切波在各层深度上的明显延迟。

III)确定上述快剪切波和上速慢剪切波的自然极化方向。

IV)对与慢剪切波极化方向一致的数据分量作一静平移。

V)对既不与慢剪切波极化方向一致又不与快剪切波极化方向一致的分量作上述静平移量一半的平移以修正方位各向异性。

VI)将最浅层下面层的所有上述数据旋转一个从指明了极化变化的上述最浅层中确定的方位角度，以识别上述最浅层中的方位各向异性主方向。

b)获得上述最浅层下面至少一层的快极化轨迹包迹幅值的一个量度。

c)获得混合/不匹配极化轨迹中至少一个的包迹幅值的量度以识别裂纹方位。

d)比较步骤(b)和(c)所获得量度。

22.如权利要求21所述的方法，其特征在于步骤(b)和(c)所述量度是幅值的平均值；步骤(d)是通过形成至少一个混合/不匹配极化轨迹的上述幅值的平均值与快极化轨迹包迹的幅值平均值的比来完成。

23.如权利要求21所述的方法，其特征在于步骤(d)包括用快极化轨迹的包迹幅值的最度对至少一个混合/不匹配极化轨迹的幅值量度进行正则化的步骤。

24.如权利要求21所述的方法，其特征在于上述静平移由上述层中慢方向极化轨迹相对于快方向极化轨迹的总延迟Δt来代表。

25.一种分析表面反射地震剪切波数据以评估剪切波极化随深度变化的方法包括分离层的步骤，该步骤中通过对与剪切波极化方向相一致的数据分量作一静平移以修正方位各向异性来指明极化的变化，一个改进的方法包括步骤：

a)对与慢剪切波极化方向一致的数据分量作一静平移。

b)对既不与慢剪切波极化方向一致又不与快剪切波方向一致的数据分量作上述静平移一半的平移，其中上述静平移代表着上述层中慢方向极化轨迹相对于快方向极化轨迹的总时延迟。

26.权利要求25方法，还包括步骤：

c)将主时间序列数据组旋转一个相应于最浅层方位各向异性方向与后一层方位各向异性方向之差的角度。

27.如权利要求26所述的方法，其特征在于剪切波数据组是从一包括至少一层相对薄的煤层的多个随深度变化的方位各向异性层中采集的，其特征还在于包括步骤：

d)对于向下至到上述煤层的各层重复(a)到(c)的步骤。

e)获得上述煤层快极化数据分量包迹幅值的一个量度。

f)将上述煤层的上述数据分量从主方向上旋转开45°。

g)获得上述煤层混合/不匹配极化数据分量中之一的包迹幅值的最度。

h)采用一个混合/不匹配极化数据分量的包迹幅值的量度与快极化数据分量包迹幅值的量度的比值来评估所述煤层的裂纹方向。

28.一种分析从一般的各向同性覆盖层和至少一层相对薄的碳质矿床中采集的表面反射地震剪切波数据组的方法，包括下面步骤：

a)获得碳质矿床快极化数据分量的主方向。

b)将上述碳质矿床的数据分量旋转到偏离主方赂大约45°的方向。

c)获得步骤(b)混合/不匹配极化数据分量中之一的包迹幅值的量度A₄₅。

29.如权利要求28所述的方法，还包括步骤：在完成(b)步骤前，获取碳质矿床快极化数据分量包迹幅值的量度A₀；并且通过形成步骤(c)所述量度与碳质矿床快极化数据包迹幅值的量度A₀的比值来评估碳质矿床的裂纹。

30.如权利要求28所述的方法，其特征在于碳质矿床是一包含多个割理的煤层；并且还包括通过上述煤层的泊松比和割理密度来校正上述比值的步骤。

31.如权利要求30所述的方法，其特征在于所述割理由环形裂缝表征；并且割理密度为环形裂缝直径的体积平均值。

32.如权利要求30所述的方法，其特征在于所述割理包括碎割理和面割理，并且割理密度代表的是面割理的容量。

33.一种评估包含有各向异性覆盖层和至少一个碳质层的碳质系中各向异性的方法，包括步骤：

a)采集碳质系上一多分量表面反射地震剪切波数据。

b)通过下述步骤分离覆盖层效应：

识别上述数据组与慢剪切波极化方向一致的数据分量，与快剪切波极化方向一致的数据分量，既不与慢剪切波极化方向又不与快剪切波极化方向一致的不匹配数据分量。

确定快剪切波轨迹与慢剪切波轨迹间的明显时间延迟。

确定上述快剪切波轨迹和上述慢剪切波轨迹的自然极化方向，

对与慢剪切波极化方向一致的数据分量做一静平移。

对即不与慢剪切波极化方向一致又不与快剪切波极化方向一致的分量作一上述静平移量的部分的平移以修正覆盖层中的方位各向异性。

将指明了极化变化的最浅层下面的上述数据组的所有数据旋转一由极化变化指明了的上述碳质矿床上面的最浅层确定的方位角度以建立碳质矿床的各向异性主方向。

c)获取碳质矿床快极化轨迹包迹幅值A₀的一个量度。

d)将指明了极化变化的上述最浅层下的所述数据组的所有数据从碳质层中各向异性主方向旋转开一锐角。

e)获取碳质矿床中至少一个混合/不匹配极化轨迹的包迹幅值A₀的量度。

f)获取至少一个混合/不匹配极化轨迹包迹值A₀的量度与快极化轨迹的包迹幅值A₀的量度的比值。

34.如权利要求33所述的方法，其特征在于碳质矿床是一个包含多个割理的煤床；并且包括通过测量煤层的泊松比和割理密度来校正上述步骤(f)的比值。

35.如权利要求34所述的方法，其特征在于所述割理由环形裂纹来表征；其中所述割理密度是环形裂纹的容量表示。

36.如权利要求34所述的方法，其特征在于所述割理包括碎割理和面割理，其中割理密度是面割理的容量表示。

37.如要求33所述的方法，其特征在于步骤(d)所述锐角在30°到60°之间。

38.如权利要求37所述的方法，其特征在于步骤(d)所述锐角是45°。

39.如权利要求33所述的方法，其特征在于步骤(c)所述碳质矿床的快极化轨迹包迹的幅值A₀的量度是其平均值的表示。

40.如权利要求33所述的方法，其特征在于步骤(e)所述碳质矿床的至少一个混合/不匹配极化轨迹的包迹幅值A₀的量度是其平均值的表示。

41.如权利要求33所述的方法，其特征在于步骤(b)包括对时间小于覆盖层底部时间的分量进行静噪。

42.如权利要求33所述的方法，其特征在于步骤(b)中所述静平移动一部分是在3/8和5/8之间。

43.如权利要求42所述的方法，其特征在于所述静平移量的一部分量一半。

44.一种采用从煤层系上采集的多分量表面反射地震剪切波数据组来评估煤层系裂纹的方法，包括步骤：

a)将指明了极化变化的上述最浅层下的所述数据组的所有数据，旋转一由极化变化指明了的上述煤层系上面最浅层确定的方位角度以建立煤层系的各向异性主方向。

b)获取步骤(a)所述煤层快极化轨迹包迹的幅值平均值的量度。

c)将指明了极化变化的最浅层深度下面的数据组的所有数据从所述煤层的各向异性主方向旋转开大约45°。

d)获取步骤(c)煤层中至少一个混合/不匹配极化轨迹包迹的幅值平均值的量度。

e)获取上述至少一个混合/不匹配极化轨迹包迹的幅值的量度对快极化轨迹包迹幅值的量度的比值。

45.如权利要求44所述的方法，其特征在于煤层系含有一各向异性的覆盖层，并且在完成步骤(a)之前，从上述向下至到煤层的数据组中分离出覆盖层的效应。

46.如权利要求45所述方法，其特征在于上述覆盖层效应从数据中的分离是如下进行：

识别与慢剪切波极化方向一致的数据分量，与快剪切波极化方向一致的数据分量以及与慢剪切波极化方向和快剪切波极化方向均不一致的不匹配数据分量。

确定快剪切波轨迹和慢剪切波轨迹之间的明显的延迟。

确定快剪切波轨迹和慢剪切波轨迹的自然极化方向。

对与慢剪切波极化方向一致的数据分量做一静平移。

对与慢剪切波极化方向和快剪切波极化方向均不一致的分量进行上述静平移量一半的平移。

对时间小于覆盖层底部的数据组分量静噪。

47.如权利要求44所述的方法，其特征在于步骤(a)的煤层包括多个碎割理和面割理；并且包括将上述比值作为至少一个上述煤层的泊松比和煤层内面割理体积密度的函数来校正。

48.一种评估地表随深度变化的方位各向异性层上表面反射地震数据的方法，包括步骤：

I)从层中采集一个至少两个分量的剪切波数据堆积，上述数据堆积包括多个地震轨迹；

II)利用上述步骤(I)的数据堆积识别包括最浅层在内的各层的快主时间序列，慢主时间序列和层的方位各向异性方向。

III)识别步骤(II)层的底部时间。

IV)利用步骤(II)的快主时间序列和慢主时间序列确定步骤(III)层的方位各向异性的程度。

V)复制该层地震分量以形成一多分量数据组。

VI)考虑该层的方位各向异性的效应对上述数据组慢方向极化轨迹进行修正。

VII)对于时间小于步骤(III)确定的该层底部时间的上述数据组的原始轨迹静噪。

VIII)利用(VII)步骤的数据组识别后一层的快主时间序列，慢主时间序列和方位各向异性。

IX)对后一层重复步骤(II)到(VIII)。

49.如权利要求48所述的方法，其特征在于步骤(I)的数据堆积沿一勘察线采集；步骤(II)是通过将上述至少两分量数据组旋转一个相应于上述勘察线方向与上述层方位各向异性方位0之差的角度来完成。

50.如权利要求49所述的方法，其特征在于步骤(IV)包括确定该层慢主时间序相对于快主时间序的总延迟Δt的量度的步骤。

51.如权利要求50所述的方法，其特征在于步骤(VI)包括将上述慢主时间序列在时间上向后平移上述延迟Δt以及用依赖于角度的比例来修正上述轨迹的步骤。