层析反演静校正手册第一期

1、第 1 期回折波层析反演静校正方法绿山地球物理公司北京代表处使用手册北京市海淀区中关村南大街 1号友谊宾馆苏园61132-33室 邮编 100873电话 (10) 6845.8445 传真 (10) 6849.8375Email : 方法介绍Fathom Tomography(FathTomo)回折波层析反演方法是绿山公司最新开发的静校正应用模块，是对现有绿山静校正软件包的补充和扩展。与现有的延迟时方法相比，层析方法提供了一种不同于折射模型的静校正计算方法。Fathom将折射信息分解到炮点和检波点的延迟时间和折射层的速度，而FathTomo方法则不同。差别重要在于对地球物理和地学观点的基本模型假设不同，层析将地球看作更复杂的模型。层析方法可以在延迟时方法应用困难的地区得到更好的近地表模型和静校正量。Fathom方法针对近地表模型进行逐步的分析和解释。软件包中的每一个模块都起到决定作用。FathTomo不规定用户对算法的细节和地球物理关系有具体的了解，只要有了初始模型和初至时间，即能求解得到最终模型。整个过程不需要进行人为干预，同时，用户需要输入或控制的参数较少。使用建议Fathom T

2、omography可以对任何观测系统求解模型和计算静校正量。但是，对下列观测系统类型，我们建议在GeoScribe II模块中为观测系统提供附加的测量信息，这样才干求得更好的静校正值。当Fathom Tomography模块打开数据库后，这些附加的信息将自动参与计算。炮点为地面炸药激发的观测系统。在G2中，用户需提供True Z Coordinate of the Shot Point Location和True Z Coordinate of the Receiver Point Location。炮点在井下炸药激发的观测系统。在G2中，用户需要提供True Z Coordinate of the Shot Point Location、True Z Coordinate of the Receiver Point Location、Uphole Time、Shot Depth - Depth of the “Down Hole” Explosive。（OBDC）海底电缆观测系统，在近水面激发。用户需提供Z Coordinate Value of “0” for the Shot P

3、oint Location、Shot Depth - Depth below the Water Surface for Source、Water Depth at each Shot Point Location、Z Coordinate at Water Bottom for each Receiver Location。海上拖缆、近水面激发观测系统。在G2中，用户需提供Z Coordinate value of “0” for the Shot Point Location、Z Coordinate value of “0” for the Receiver Point Location、Water Depth at each Shot/Receiver Point Location、Shot Depth(Depth below the Water Surface) for the Source、Receiver Depth(Depth below the Water Surface) for the Receiver。回折波层析反演原理FathTomo是一个全三维回折波层析反演

4、方法，通过建立近地表速度模型，计算静校正量，合用于二维和三维资料。该方法是一个速度反演过程。对任何观测系统，用回折波或连续折射直达波，交互反演近地表的速度变化。层析速度算法涉及多个环节：（a）初至时间拾取；（b）模型离散化和初始速度；（c）射线追踪和分割；（d）剩余时间（误差）计算；（e）速度更新，减小误差。简朴流程如下：初始模型射线追踪层 析计算剩余时间输出模型初至时间层析运算涉及一个正演过程即计算每个炮检对的旅行时间，和一个反演过程即速度得到交互更新，生成速度模型。通过若干次迭代（每次迭代都是一个射线追踪、剩余时间计算和速度更新的周期）得到一个平滑的速度模型。一般来讲，根据不同的初始模型和不同噪音水平的初至时间，通过8-10次迭代后，误差都可以得到收敛。层析反演流程和参数（一）大炮初至拾取层析反演方法需要炮点和检波点的坐标信息和大炮初至时间。用GeoScribe II模块建立观测系统数据库，用Picker模块拾取地震道的初至时间，或者从其他解决系统输入拾取的初至时间。拾取的初至时间代表观测到的通过模型的旅行时间，用于计算剩余（误差）时间。（二）模型空间离散化和初始速度模型定义模型空

5、间离散化涉及，第一步将三维模型空间提成一系列的网格（Voxel）。网格大小可以人为选择，一般而言，网格应当足够小才干分辨感爱好的最小地质体；同时，网格也应当足够大，才干获得充足的射线，增长成像的可信度。第二步需要定义初始模型速度。在抱负情况下，可以根据该区域已知的地质和地球物理信息来定义初始模型速度。模型成像域网格的数量是反演中未知参数，炮检对的个数拟定了已知的参数（旅行时间）这个过程需要在射线追踪之前完毕。图(1)是层析中定义模型网格的界面。缺省值是根据观测系统和Bin（3D）或CMP（2D）大小计算出来的。在这里给出合适的参数非常很重要，保证模型可以对的离散化。模型网格大小在空间上不一定规定是正方体，在Inline、Crossline和垂直方向都可以不同样，这是绿山层析方法很突出的优点。由于对于某些地震观测系统，假如网格在Inline和Crossline方向的大小相同时，计算结果会产生假像。本地震观测系统的接受线距与道距之比较大时，长方形网格比较好，由于射线穿透的深度与排列长度相比较小，通常网格的垂直大小应当小于Inline和Crossline方向的大小比较有利。模型网格定义完毕之

6、后，要为每一个网格的中心赋予速度值，射线追踪需要的网格中心之间的速度值可以通过不均匀内插获得。注意：模型速度的剧烈不连续性会导致射线分散，不适合层析方法中所用的射线追踪算法。初始速度模型可以根据现有的地质信息，如测井资料来建立，或运用梯度变化的速度模型。通常而言，运用合适的初始速度模型开始迭代可以提高最终成像的可信度。此外，初始速度模型应当保证地震能量可以从炮点传播到检波点。例如，常速模型不适合地表地震采集观测系统，由于从炮点到检波点没有回折射线。图(1). 模型空间定义图(2)表达层析中定义速度模型的界面。假如选择梯度速度模型，用户需要输入开始的速度和梯度变化因子。初始速度模型也可以通过输入GeoScribe II数据库中的折射速度信息，或输入ASCII格式的速度信息以及手工定义来完毕。图(2). 定义初始速度模型（三）射线追踪和分割射线追踪和分割涉及计算每个炮检对的旅行时间和射线途径，以及根据射线途径获得通过每一个网格的射线段。射线段的长度在反演问题中用于计算加权和更新速度。根据不同的物理模型和数字模型，有多种射线追踪算法。在绿山的层析方法中，我们采用Um和Thurber于1987

7、年提出的最大速度梯度射线追踪三维算法，这是一种两点射线追踪方法。我们之所以采用这种算法，是由于它的计算效率以及射线追踪方法的优点，可以避免内插。这种方法根据费马原理（Fermats Principle），在炮点和检波点之间通过计算最小的旅行时间，找到两点之间的射线途径，而不是严格地验证Snell定律。对每一条射线途径，都是一个迭代的计算过程。网格中射线段的数量根据射线的数量成倍增长，直到用户设定的精度规定。射线途径的数量按照3，5，9，17，33急剧增长，假如用户设定了最大值，最终将达成65。增长射线数量会减少计算效率，但可以在炮点和检波点之间得到一个较平滑的射线途径。这种算法不规定有岩性边界或水平连续层面。计算出炮点和检波点之间一条射线途径之后，接着对射线进行分割，这个过程涉及计算射线途径长度，和该射线通过的网格在Inline、Crossline 及深度方向的相应的坐标。在层析中，这个信息为每个网格的速度更新提供加权系数。图(3)表达层析中设立批量解决参数的界面。所有层析迭代的参数都在这里定义。图(3). 设立批量解决参数用户可以更改射线追踪的缺省参数，如图(4)：图(4). 射线追

8、踪参数Minimum Ray Path Points对每一个计算途径，定义最小的射线途径反射点数，缺省为33。假如该值较小，射线追踪速度会更快；假如该值较大，射线途径将会更平滑。Maximum Difference on Travel Time射线追踪采用Perturbation Scheme方案。射线段数量增长一倍，即跳到下一个位置点。计算每个炮、检对总的旅行时间，和前一个计算时间对比，假如精度达成用户的设定，那么，该炮、检对的射线追踪停止。用户也许会注意到，缺省的精度值是0.001毫秒，与我们通常的初至时间相比，这已经是一个很高的精度了。但这个参数在这里是规定射线追踪达成一个合理的收敛，而不是层析所要的精度。假如用户指定这个值为1，那么就无法得到对的的回折波途径。经验表白，选0.001或更小的值，可以得到较好的效果。Maximum Perturbation Allowed(ft or m)如上所述，射线段数量增长一倍，即跳到下一个位置点。在本次迭代中，当前坐标与前一次迭代相比，最大误差不能超过用户指定的这个值。缺省为网格深度方向大小的2.5倍。假如起始速度模型的梯度较大，则选用网格

9、深度大小的1.5倍较合适。Shot Decimation Increment对于一块较大的三维资料来讲，用户也许会减少参与层析运算的炮点的数量，缺省为所有参与。假如选3，则代表数据库中每隔三炮有一炮参与层析计算。用户可以选择数据库中所有的初至都参与计算，也可以指定由Branch模块定义的某折射层的初至参与计算。两种选择各有优点。假如选择用所有的初至，算法将会把初至按照偏移距的大小顺序采用。这与射线途径多次覆盖才干产生收敛速度模型的假设一致。假如选择某折射层的初至参与计算，算法只采用相应偏移距范围的初至，然而，层与层之间交叉区域的初至质量会较差，这些较差的初至会干扰最终的层析效果。此外，对某些观测系统，在层析方法中限制偏移距范围也许会得到较好的结果。这种做法就是只采用指定层的初至，通常是最浅的折射层。用户也可以用此外一种方法对射线追踪进行限制，如图(5)。图(5). 定义参与射线追踪的初至射线追踪可以根据偏移距范围、初至时间范围和炮、检点的方位角范围来进行限制。缺省情况下，对此无限制。用户可以根据观测系统和指定的偏移距范围内初至的质量，进行实验。（四）剩余时间（误差）计算剩余时间是指观测到的初至时间（拾取的初至）与计算旅行时间的差异。计算旅行时间是根据射线途径通过的每个网格的旅行时间求和得出的。（五）速度更新层析过程的最后一步是更新每个网格的速度，根据迭代背投图表，采用同步迭代重建技术（SIRT）。在批量解决参数设立中，可以改变缺省的反演参数。如图(6)。Inversion Loops速度更新依据迭代背投图表。一次反演循环指所有的炮、检对的单次过程。射线从炮点到检波点，要通过很多网格。通过每个网格的射线段的数量在单个网格速度更新时起到加权作用。射线途径的误差被分派到射线通过的每个网格中。对所有的炮、检对进行连续解决，计算结果是所有射线段通过每个网格的平均网格速度。这个过程反复循环，直到

《层析反演静校正手册第一期》由会员大米分享，可在线阅读，更多相关《层析反演静校正手册第一期》请在金锄头文库上搜索。