静校正(讨论版)

1、四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司静校正处理静校正处理（黄大云）（黄大云）2008.4四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司内容内容 静校正概述静校正概述 高程静校正高程静校正 折射静校正折射静校正 层析静校正层析静校正 剩余静校正剩余静校正四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司地震资料常规处理一般流程地震资料常规处理一般流程数据输入数据输入观测系统加载观测系统加载坏炮坏道去除坏炮坏道去除噪声压制噪声压制振幅补偿振幅补偿静校正静校正反褶积反褶积叠加叠加剩余静校正剩余

2、静校正速度分析速度分析动校正动校正偏移偏移四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 问题问题:在地震资料处理中为什么要做静校正在地震资料处理中为什么要做静校正?静校正概述静校正概述 1、 为使地下构造形态不受地表起伏的影响为使地下构造形态不受地表起伏的影响; 2、一个探区需使用统一的基准面、一个探区需使用统一的基准面,以便测线之间接点以便测线之间接点; 3、地震资料处理的一些重要步骤是在反射波时距曲线为双、地震资料处理的一些重要步骤是在反射波时距曲线为双 曲线的前提下进行的（速度分析、动校正等），但反射曲线的前提下进行的（速度分析、动校正等），但反射 波时距曲线为双曲线的条

3、件是：地表水平、上覆介质速度波时距曲线为双曲线的条件是：地表水平、上覆介质速度 为常数。为了后续处理的需要，应把反射时间校正到炮点、为常数。为了后续处理的需要，应把反射时间校正到炮点、 检波点均在一个水平面上一样的情况。检波点均在一个水平面上一样的情况。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司22202vxttP PV VS SM MxR R021txt21（双曲线方程）（双曲线方程）v地表水平、上覆介质速度为常数时地表水平、上覆介质速度为常数时,地震波传播路径如下图地震波传播路径如下图:四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 在地震勘探中，由于受地形

4、起伏、爆炸井深不一、低在地震勘探中，由于受地形起伏、爆炸井深不一、低降速带的厚度和速度变化等因素的影响，地震波的反射时降速带的厚度和速度变化等因素的影响，地震波的反射时距曲线已不是理论上的双曲线分布。如何对记录的反射时距曲线已不是理论上的双曲线分布。如何对记录的反射时间进行校正，使反射时距曲线还原为双曲线，这就是静校间进行校正，使反射时距曲线还原为双曲线，这就是静校正要解决的主要问题。静校正通常包括三个方面：正要解决的主要问题。静校正通常包括三个方面：1 1、野外高程变化引起的时差校正；、野外高程变化引起的时差校正；2 2、低、降速带引起的时差校正；、低、降速带引起的时差校正；3 3、以上两种

5、校正留下的残余和其它因素引起的剩余时差校正。、以上两种校正留下的残余和其它因素引起的剩余时差校正。 四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司目前我们常用的静校正方法目前我们常用的静校正方法高程静校正高程静校正折射静校正折射静校正层析静校正层析静校正剩余静校正剩余静校正四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 解决地形起伏、爆炸井深不一引起的静校正问解决地形起伏、爆炸井深不一引起的静校正问题可用野外高程校正方法解决。该方法利用野外测题可用野外高程校正方法解决。该方法利用野外测量成果和预定的基准面高程及其上覆层速度来计算量成果和预定的基准面高程及其上覆层速度

6、来计算校正量。校正量。高程静校正高程静校正四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司P PV VS1S1R1R1M M基准面地面S SR R炮点校正量检波点校正量高程校正示意图高程校正示意图误差来自于把三角形的一个边误差来自于把三角形的一个边的长度看作另两边之和，其误的长度看作另两边之和，其误差随基准面离地面的距离和炮差随基准面离地面的距离和炮检距的增加而变大。检距的增加而变大。这种计算校正量的这种计算校正量的方法有时误差太大方法有时误差太大四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 通常不直接校正到标准基准面而是先校正到通常不直接校正到标准基准面而是先校正

7、到浮动浮动基准面，以后再校正到标准基准面。基准面，以后再校正到标准基准面。问题问题: :1 1、为什么要使用浮动基准面、为什么要使用浮动基准面? ?2 2、在什么情况下须使用、在什么情况下须使用浮动浮动基准面基准面? ?3 3、数据在什么时候校正到固定基准面、数据在什么时候校正到固定基准面? ?四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司P PV VS1S1R1R1M M基准面地面S SR R炮点校正量检波点校正量高程校正示意图高程校正示意图浮动基准面S2S2R2R2 浮动基准面与地面浮动基准面与地面高差较小高差较小, ,可降低静可降低静校正计算误差校正计算误差问题问题1:为什

8、么要使用浮动基准面为什么要使用浮动基准面?四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 浮动基准面在局部范围内可看作是平面，它可能是水平面，也可能是倾浮动基准面在局部范围内可看作是平面，它可能是水平面，也可能是倾斜面。如果是水平面，静校正的目的已经达到。如果是倾斜面，则反射波斜面。如果是水平面，静校正的目的已经达到。如果是倾斜面，则反射波走时图如左图。我们将该图做旋转，使倾斜的浮动基准面变为水平，原本走时图如左图。我们将该图做旋转，使倾斜的浮动基准面变为水平，原本水平的反射界面就变为倾斜。这样一来，我们就把倾斜基准面的问题变成水平的反射界面就变为倾斜。这样一来，我们就把倾斜基准

9、面的问题变成了倾斜反射面的问题。而了倾斜反射面的问题。而倾斜反射面的时差校正将由倾斜反射面的动校正去解决。倾斜反射面的时差校正将由倾斜反射面的动校正去解决。 如果地下反射界面原本就是倾斜的，旋转后变成什么样呢？如果地下反射界面原本就是倾斜的，旋转后变成什么样呢？SRRS 静校正的目不是要将炮点、检波点校正到同一水平面上吗？静校正的目不是要将炮点、检波点校正到同一水平面上吗？只校正到浮动基准面，动校正的条件满足吗？只校正到浮动基准面，动校正的条件满足吗？四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 通常用炮点、检波点的平滑面作为浮动基准面，平滑长度通常用炮点、检波点的平滑面作为浮

10、动基准面，平滑长度一般为排列长度。浮动基准面静校正只计算炮点、检波点到浮一般为排列长度。浮动基准面静校正只计算炮点、检波点到浮动基准面的校正量，以后再计算浮动基准面到标准基准面的校动基准面的校正量，以后再计算浮动基准面到标准基准面的校正量。炮点、检波点实际静校正量等于两部分相加。正量。炮点、检波点实际静校正量等于两部分相加。 当地表起伏太大时须使用浮动基准面当地表起伏太大时须使用浮动基准面!问题问题3:数据在什么时候校正到固定基准面数据在什么时候校正到固定基准面?问题问题2:在什么情况下须使用浮动基准面在什么情况下须使用浮动基准面?四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司三

11、种可能的选择三种可能的选择: 1、动校前、动校前 2、动校后、动校后 3、叠加后、叠加后1 1、动校前、动校前: :与直接使用固定基准面差别不大与直接使用固定基准面差别不大2 2、动校后、动校后, ,即速度分析和动校正在浮动基准面上进行即速度分析和动校正在浮动基准面上进行: :浮动基准面是平滑的浮动基准面是平滑的, ,在局部范围内可看成直线在局部范围内可看成直线, ,这时的静这时的静校校正相当于倾斜地层的动校正正相当于倾斜地层的动校正,倾角为浮动基准面倾角和地层倾角为浮动基准面倾角和地层倾角的综合倾角的综合。动校后动校后CMP道集内同相轴拉平道集内同相轴拉平,此时校正到固此时校正到固定基准面后

12、同相轴会有时差定基准面后同相轴会有时差,但因浮动基准面比地表平滑但因浮动基准面比地表平滑,时时差不及差不及直接使用固定基准面那么大直接使用固定基准面那么大, ,有时有时或许或许可用剩余静校可用剩余静校正解决。正解决。3 3、叠加后、叠加后: :理论上叠加成像最好理论上叠加成像最好, ,但这时已没有炮点但这时已没有炮点、检波点的检波点的 概念概念, ,校正须按校正须按CMPCMP点进行。点进行。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 野外高程校正中的替换速度为常数，当探野外高程校正中的替换速度为常数，当探区内无低降速带变化或变化很小可忽略不计时，区内无低降速带变化或变化很小

13、可忽略不计时，使用野外高程校正即可解决静校正问题。但对使用野外高程校正即可解决静校正问题。但对于复杂地区而言，情况就不那么简单，因为在于复杂地区而言，情况就不那么简单，因为在这些地区，低降速带的厚度和速度通常都存在这些地区，低降速带的厚度和速度通常都存在变化，这时就不能使用高程静校正了。变化，这时就不能使用高程静校正了。 四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司折射静校正折射静校正 如果低降速带存在，则它与高速层的分界面如果低降速带存在，则它与高速层的分界面是一个很好的反射界面。当入射波以临界角入射是一个很好的反射界面。当入射波以临界角入射到该界面时，将产生折射波。折射波传

14、到地面为到该界面时，将产生折射波。折射波传到地面为检波器所接收，根据折射波的初至时间可估算低检波器所接收，根据折射波的初至时间可估算低降速带的速度和炮点及检波点的延迟时间，这就降速带的速度和炮点及检波点的延迟时间，这就是折射静校正的基本思路。是折射静校正的基本思路。 四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司A点激发，点激发，B点接收的折射波到达时间点接收的折射波到达时间TAB为：为： TAB = AC/VR+CD/VS+DB/VR 或或 TAB = AH/（VR cos）+AB-（AH+BG）tan/VS+BG/（VR cos） 考虑到考虑到sin=VR/VS ,且且VR

15、VS , 上式变为上式变为 TAB = TA +AB/ VS + TB （折射波旅行时方程）（折射波旅行时方程） 其中，其中，TA = AHcos/VR 为炮点的延迟时间为炮点的延迟时间 (时间深度时间深度)， TB = BGcos/VR 为检波点的延迟时间为检波点的延迟时间(时间深度时间深度) 。 对以上对以上折射波旅行时方程可用互换法或迭代法求解。折射波旅行时方程可用互换法或迭代法求解。 折射示意图折射示意图 临界角：sin=VR/ VS 地层速度：VSABHGCD低降速带速度：VR地面折射面四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司因为临界角因为临界角：当低降速带速度当

16、低降速带速度 与地层速度与地层速度 差别不大时差别不大时,临界角会临界角会很大很大,在有限的排列内可能收不到折射或折射波到达时在有限的排列内可能收不到折射或折射波到达时,其其它波已先行到达它波已先行到达, 所以低降速带速度与地层速度的差别必所以低降速带速度与地层速度的差别必须足够大才能在记录上拾取到折射波初至。须足够大才能在记录上拾取到折射波初至。 SRVV /sinRVSV问题问题1:在什么情况下能在记录上拾取到折射波初至在什么情况下能在记录上拾取到折射波初至?四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 折射静校正的两个必要条件：折射静校正的两个必要条件： 1 1、存在稳定

17、的折射层；、存在稳定的折射层； 2 2、能够从炮记录的初至中拾取到折射波初至。、能够从炮记录的初至中拾取到折射波初至。 如果不存在稳定的折射层或者或不能拾取到折射波初至，如果不存在稳定的折射层或者或不能拾取到折射波初至，则无法实现折射静校正。则无法实现折射静校正。 四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司问题问题2:折射波初至在记录上有什么特征折射波初至在记录上有什么特征?为什么为什么? 与直达波初至比较与直达波初至比较, ,相邻道时差小相邻道时差小, ,所以显得平直所以显得平直, ,两端与直达波相接形成拐点。两端与直达波相接形成拐点。 因为相邻道折射波到达时间之间的差由它

18、们传因为相邻道折射波到达时间之间的差由它们传播路径不同引起的播路径不同引起的, ,而不同的这段路径为高速。而不同的这段路径为高速。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司Omega系统的折射静校正模块系统的折射静校正模块 EGRM EGRMEGRM是是Extnded Generalized Reciprocal Method(Extnded Generalized Reciprocal Method(扩展扩展常规互换法常规互换法) )的简称的简称, ,该模块根据折射初至拾取时间使用互该模块根据折射初至拾取时间使用互换法或换法或( (和和) )高斯高斯- -赛德尔迭代法求解获

19、得炮点项和检波点项赛德尔迭代法求解获得炮点项和检波点项的折射时间深度和折射速度层速度信息，并由此构建出折的折射时间深度和折射速度层速度信息，并由此构建出折射层的速度深度模型，再并根据速度深度模型计算出静校射层的速度深度模型，再并根据速度深度模型计算出静校正量。正量。 总体静校正量总体静校正量= 地面到中间基准面时间地面到中间基准面时间+ +中间基准中间基准面到标准基准面的时间面到标准基准面的时间以折射层速度计算以折射层速度计算以替换速度计算以替换速度计算四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司三种提取折射信息的方法三种提取折射信息的方法: :1 1、扩展常规互换法扩展常规互

20、换法,该技术是常规互换法技术的扩展该技术是常规互换法技术的扩展 能获得准确有效的折射延迟时估算能获得准确有效的折射延迟时估算,且与折射速度场的误差相对无关且与折射速度场的误差相对无关。2、高斯高斯- -赛德尔迭代分解法赛德尔迭代分解法( (最小平方法最小平方法) ) 有很强的外推特性有很强的外推特性, ,因该方法要同时估算炮点因该方法要同时估算炮点、检波点折射延迟时和射检波点折射延迟时和射速度速度, ,所以与折射速度场的误差有关所以与折射速度场的误差有关。3 3、混合、混合EGRM/EGRM/迭代算法迭代算法 以上两种算法组合以上两种算法组合。 三种方法既可单独使用三种方法既可单独使用, ,也

21、可按某种组合使用也可按某种组合使用。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司扩展常规互换法扩展常规互换法四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司互换法折射速度计算互换法折射速度计算AG1G2G3G4G5V1VnVn+1V1VnBG1G2G3G4G5Vn+1DT1=TAG1-TBG1= 负值负值DO1=AG1-BG1= 负值负值DT3=TAG3-TBG3= 0DO3=AG3-BG3= 0DT5=TAG5-TBG5= 正值正值DO5=AG5-BG5= 正值正值DT= 初至时差初至时差DO= 偏移距差偏移距差A、B：炮点：炮点G： 接收点接收点四川石油管理局地

22、球物理勘探公司四川石油管理局地球物理勘探公司DTDOG1G2G3G4G5+_Vn+1= 拟合直线斜率拟合直线斜率/2DT：时间差：时间差DO：偏移距差：偏移距差四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司互换法延迟时间互换法延迟时间(以时间表示的深度以时间表示的深度)计算计算XYABV1V2DADBDR121VDBVDRVDATAB当当V2比比V1大得多时大得多时,上式可近似地改为上式可近似地改为: （1）TAB=折射波到达时间；折射波到达时间； ZA=A点的折射层厚度；点的折射层厚度；ZA= B点的折射层厚度；点的折射层厚度； = A到到B点的地面距离；点的地面距离； =临界

23、角。临界角。这说明折射波到达时间可表示为这说明折射波到达时间可表示为A、B两点的时间深度之和再加上两点的时间深度之和再加上AB的距离的距离/V2。121VCOSZBVABVCOSZATABAB四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 某一站点某一站点G的时间深度的时间深度TD(G)需用需用A、B两个炮点的三个折射旅行时两个炮点的三个折射旅行时: （2） TAB= A点激发点激发B点接收的折射时间；点接收的折射时间； TAG= A点激发点激发G点接收的折射时间；点接收的折射时间； TBG= B点激发点激发G点接收的折射时间。点接收的折射时间。2)()(TABTGBTAGGTD

24、AGBV1临界角临界角V2四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 并非所有点都能使用互换法求取其时间深度，如测线两端、炮点不规则并非所有点都能使用互换法求取其时间深度，如测线两端、炮点不规则段等。对这些点必须使用其它办法段等。对这些点必须使用其它办法: :1 1、Time-depth ProjectTime-depth Project 利用相邻点的时间深度做最小平方（线性）拟合利用相邻点的时间深度做最小平方（线性）拟合2 2、Refractor Elevation ProjectRefractor Elevation Project 先将先将相临点的时间深度转换为折射层高

25、程，再做最小平方（线性）拟合相临点的时间深度转换为折射层高程，再做最小平方（线性）拟合3 3、Overlay Extrapolation of Time-DepthOverlay Extrapolation of Time-Depth 将拾取的初至时间依次连续扩展到整条测线，这些扩展时间就叫做重复时将拾取的初至时间依次连续扩展到整条测线，这些扩展时间就叫做重复时间（间（overlay-times），重复时间与要估算时间深度的区域的临近站点的时），重复时间与要估算时间深度的区域的临近站点的时间进行最小平方意义下的比较，并按最小平方拟合线扩展到估算区域。这种间进行最小平方意义下的比较，并按最小平方

26、拟合线扩展到估算区域。这种方法等同于外插和内插地结合。方法等同于外插和内插地结合。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司折射模型的建立折射模型的建立 为了建立折射模型，先将时间为了建立折射模型，先将时间-深度表达式做一下变换：深度表达式做一下变换： 将式（将式（1）代入式（）代入式（2），经一系列计算并考虑到），经一系列计算并考虑到 最后得到：最后得到： （3） 其中，其中，ZG为折射层厚度为折射层厚度,在式（在式（3）中，）中，ZG、V1和临界角和临界角 都是未知的，但临界角可由都是未知的，但临界角可由V1、V2求得。前求得。前面我们用直线拟合方式求得了面我们用直线拟合

27、方式求得了V2,所以实际上只有所以实际上只有ZG和和V1两个未知数，而折射层两个未知数，而折射层的速度的速度-深度正是由它们确定。为了得到深度正是由它们确定。为了得到ZG和和V1，我们再来看另一个观测值：临，我们再来看另一个观测值：临界距离，所谓临界距离即折射波最先道达那一道的炮检距。临界距离界距离，所谓临界距离即折射波最先道达那一道的炮检距。临界距离XC与临界角与临界角的关系是的关系是:GBAGAB1cos)(VZGGTD21sin1VV)4(tan2ZGXC四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司式（式（4）说明，在）说明，在V2已知的前提下，临界距离已知的前提下，临界

28、距离XC也是也是ZG和和V1的函数。于是由的函数。于是由 和和两个方程可最终求得两个方程可最终求得ZG和和V1。有了各站点（。有了各站点（G）的厚度和速度）的厚度和速度V1，折射层的速度，折射层的速度深度模型也就确定了。深度模型也就确定了。1cos)(VZGGTDtan2ZGXC四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司高斯高斯-赛德尔迭代旅行时分解赛德尔迭代旅行时分解 该方法就是用迭代分解的办法计算各站点折射层的时间深度，有了时间深度，即可该方法就是用迭代分解的办法计算各站点折射层的时间深度，有了时间深度，即可按照前面的方法得到折射层的速度深度模型。式（按照前面的方法得到折

29、射层的速度深度模型。式（1）表示折射旅行时是炮点、检波点时间深度之和再加上偏移距表示折射旅行时是炮点、检波点时间深度之和再加上偏移距/V2，为便于叙述，可写成：，为便于叙述，可写成：显然迭代计算公式就是显然迭代计算公式就是 其中，其中，TA =炮点炮点A的时间深度，的时间深度，TB=检波点检波点B的时间深度，的时间深度， 称为构造项或偏移称为构造项或偏移距距-速度项速度项,由于由于V2可通过前面的互换法求得，这一项可当作已知，这时迭代只需前两可通过前面的互换法求得，这一项可当作已知，这时迭代只需前两个公式即可，若使用三个公式迭代，则用互换法求得的构造项用作迭代初值。个公式即可，若使用三个公式迭

30、代，则用互换法求得的构造项用作迭代初值。121VCOSZBVABVCOSZATABTBVABTATAB2TBVABTABTA2TAVABTABTB2TBTATABVAB22VAB四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司TA和和TB的初值估算的初值估算 为了获得某个站点的初始时间深度，模块要把以该点为中点的所有折射旅行时收为了获得某个站点的初始时间深度，模块要把以该点为中点的所有折射旅行时收集起来：集起来： 如图所示，要得到站点如图所示，要得到站点120的初始时间深度，就把炮点的初始时间深度，就把炮点110到检波点到检波点130、炮点、炮点134到检波点到检波点106、炮点、

31、炮点101到检波点到检波点139等等折射旅行时收集起来，它们显然都是等等折射旅行时收集起来，它们显然都是偏移距偏移距-速度项的函数（左图）。将偏移速度项的函数（左图）。将偏移-速度项减去，再求平均，然后除以速度项减去，再求平均，然后除以2，即，即为该站点的初始时间深度（右图）。为该站点的初始时间深度（右图）。120101 106110130 134139TBVABTATAB2四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司混合混合EGRM/迭代算法迭代算法 用互换法得到的各站点的时间深度和（或）偏移用互换法得到的各站点的时间深度和（或）偏移距距- -速度值作为初始值迭代计算出各站点

32、最终的时间速度值作为初始值迭代计算出各站点最终的时间深度和（或）偏移距深度和（或）偏移距- -速度值速度值。 迭代法和混合法的折射层速度深度模型建立与前迭代法和混合法的折射层速度深度模型建立与前述互换法相同。述互换法相同。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司静校正量计算静校正量计算 计算静校正量需用两个基准面，一个是中间基准面计算静校正量需用两个基准面，一个是中间基准面, ,一个是最终基准一个是最终基准面。任一站点的静校正量等于从该点到中间基准面的时间延迟和从中间面。任一站点的静校正量等于从该点到中间基准面的时间延迟和从中间基准面到最终基准面时间延迟的代数和。从站点到中

33、间基准面的时间延基准面到最终基准面时间延迟的代数和。从站点到中间基准面的时间延迟计算使用前面所建立的速度深度模型的速度，从中间基准面到最终基迟计算使用前面所建立的速度深度模型的速度，从中间基准面到最终基准面时间延迟的计算使用替换速度。准面时间延迟的计算使用替换速度。 EGRMEGRM将中间基准面定为折射层底界的一个平滑面。因此，若折射层底将中间基准面定为折射层底界的一个平滑面。因此，若折射层底界高程变化不大，则中间基准面为水平，若折射层底界高程变化较大，界高程变化不大，则中间基准面为水平，若折射层底界高程变化较大，则中间基准面就不是水平的。另外用户还可以用参数对中间基准面进行则中间基准面就不是

34、水平的。另外用户还可以用参数对中间基准面进行修改或控制。修改或控制。 折射层的速度深度模型可以是单层也可以是多层，但即使我们只要折射层的速度深度模型可以是单层也可以是多层，但即使我们只要求单层模型，模块也把井底以上的地表当作一层即所谓的虚拟折射层，求单层模型，模块也把井底以上的地表当作一层即所谓的虚拟折射层，所以实际上仍是两层。多层模型的静值计算公式如下：所以实际上仍是两层。多层模型的静值计算公式如下：四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司多层模型的静值计算公式多层模型的静值计算公式其中，其中， TS = 总的静校正量总的静校正量 THICKi = 第第i 层的厚度层的厚

35、度 Vi = 第第 i层的速度层的速度 ERn = 第第 n层的高程层的高程 EI = 中间基准面高程中间基准面高程 Vn+1 = 第第 n层以下的速度层以下的速度 ED = 最终基准高程最终基准高程 VR = 替换速度替换速度 TUP = 井口时间井口时间UPRDInInniiiTVEEVEERVTHICKTS11四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司主要参数主要参数一一 处理参数组处理参数组(*PROCESS_OPTIONS)1 1、算法类型、算法类型(SOLUTION_TYPE)(SOLUTION_TYPE)三个选项三个选项: :(1)(1) EGRM:EGRM:所

36、有折射层都用扩展常规互换算法所有折射层都用扩展常规互换算法(2)(2) ITERATIVE:ITERATIVE:指定层使用迭代法指定层使用迭代法, ,其余用扩展常规互换算法其余用扩展常规互换算法(3)(3) HYBRID:HYBRID:使用扩展常规互换算法使用扩展常规互换算法/ /迭代法混合算法迭代法混合算法, ,即用扩展常规即用扩展常规互换算法求得的时间深度值作为初值进行迭代互换算法求得的时间深度值作为初值进行迭代 缺省缺省:EGRM:EGRM四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司2、最终基准面高程、最终基准面高程(DATUM_ELEV) 缺省缺省: #DATUM_EL

37、EV3、替换速度、替换速度(REPL_VEL) 缺省缺省: #REPLACE_VEL4、初至拾取时间调整值、初至拾取时间调整值(BULK_SHIFT) 所有初至拾取时间都加上该值所有初至拾取时间都加上该值 缺省缺省: 0.05、时间深度偏差门槛值、时间深度偏差门槛值(VARIANCE_THRESHOLDL) 计算偏差超过该门槛值的时间深度将做加权平滑计算偏差超过该门槛值的时间深度将做加权平滑,并在迭并在迭 代计算中不用作初值。代计算中不用作初值。 缺省缺省: 16ms四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司二二 折射层参数组折射层参数组(*LAYER_CONTROL) 该参

38、数组以参数表的形式给出该参数组以参数表的形式给出1 1、重复站点最小个数、重复站点最小个数(MIN_OVERLAP)(MIN_OVERLAP) 在测线两端或炮点缺失区域不能获得时间深度或折射高程在测线两端或炮点缺失区域不能获得时间深度或折射高程, ,需用最小平需用最小平方拟合外插获得。本参数指定做方拟合外插获得。本参数指定做最小平方拟合的站点个数。最小平方拟合的站点个数。 缺省缺省:25:252 2、折射模型的层号、折射模型的层号(LAYER)(LAYER) 缺省缺省: NOT_USED(: NOT_USED(无缺省无缺省) ) 3 3、用于该折射层处理的起始站点号、用于该折射层处理的起始站点

39、号 缺省缺省:GENERAL_FIRST_STN:GENERAL_FIRST_STN4 4、用于、用于该折射层处理该折射层处理的最后站点号的最后站点号 缺省缺省:GENERAL_LAST_STN:GENERAL_LAST_STN四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司5 5、非互换估算类型、非互换估算类型(ESTIMATION_TYPE)(ESTIMATION_TYPE) 指定生成测线两端或其它不能使用互换法区域的时间深度的估算类型指定生成测线两端或其它不能使用互换法区域的时间深度的估算类型 TIME_DEPTH TIME_DEPTH 用相邻站点的时间深度做最小平方拟合外插

40、用相邻站点的时间深度做最小平方拟合外插 ELEVATION ELEVATION 用相邻站点的折射高程做最小平方拟合外插用相邻站点的折射高程做最小平方拟合外插 EXTRAP_TIME_DEPTH EXTRAP_TIME_DEPTH 用相邻站点的拾取时间做最小平方拟合外插用相邻站点的拾取时间做最小平方拟合外插, ,然后然后 计算时间深度。计算时间深度。 EXTRAP_ELEV EXTRAP_ELEV 用相邻站点的折射高程做最小平方拟合外插用相邻站点的折射高程做最小平方拟合外插, , 然后计算时然后计算时 间深度。间深度。6 6、中间站点号、中间站点号(MID_STN)(MID_STN) 指定该层内

41、的站点号指定该层内的站点号, ,它的最小偏移距、它的最小偏移距、最大偏移距和内插方式由后面的最大偏移距和内插方式由后面的参数指定。参数指定。 FIRST FIRST 第一个站点第一个站点 LASTLAST 最后一个站点最后一个站点 其它数组其它数组 ( (无缺省无缺省) )四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司7 7、最小偏移距、最小偏移距(MIN_OFFSET)(MIN_OFFSET) 该中间站点号的最小偏移距该中间站点号的最小偏移距 大于等于大于等于0 0 缺省缺省:0.0:0.08 8、最大偏移距、最大偏移距(MIN_OFFSET)(MIN_OFFSET) 该中间站

42、点号的最大偏移距该中间站点号的最大偏移距 大于等于最小偏移距大于等于最小偏移距 缺省缺省: :最小偏移距最小偏移距 由最小偏移距和最大偏移距确定该偏移距范围内的初至拾取时间将由最小偏移距和最大偏移距确定该偏移距范围内的初至拾取时间将 用作时间深度和折射速度计算。用作时间深度和折射速度计算。9 9、内插偏移距、内插偏移距 是否对下一个指定站点的偏移距做内插是否对下一个指定站点的偏移距做内插 YES YES 内插内插 NO NO 不内插不内插 插缺省插缺省: :YES 四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司三三 迭代求解参数组迭代求解参数组(*ITERATIONS) 当算法类

43、型选择了当算法类型选择了ITERATIVEITERATIVE时时, ,必须使用本参数组必须使用本参数组1 1、迭代次数、迭代次数(MUN_ITER)(MUN_ITER) 重复估算炮点和检波点项的次数。重复估算炮点和检波点项的次数。 缺省缺省:6:62 2、分解分量项顺序、分解分量项顺序 SD SD 炮点项炮点项, ,检波点项检波点项 DS DS 检波点项检波点项, ,炮点项炮点项 DSC DSC 检波点项检波点项, ,炮点项炮点项, ,构造项构造项 SDC SDC 炮点项炮点项, , 检波点项检波点项, ,构造项构造项 CDS CDS 构造项构造项, ,检波点项检波点项, ,炮点项炮点项 CS

44、D CSD 构造项构造项, ,炮点项炮点项, ,检波点项检波点项 缺省缺省: SD : SD 四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司3 3、使用迭代法的折射层号、使用迭代法的折射层号 无缺省无缺省四四 混合迭代方法参数组混合迭代方法参数组(HYBRID_METHOD)(HYBRID_METHOD)1 1、删除出现过大值的道、删除出现过大值的道(DELETE_OUTLIERS)(DELETE_OUTLIERS) 删除或保留那些在统计中出现过大值的道删除或保留那些在统计中出现过大值的道 NO NO 不删除不删除( (建议建议) ) YES YES 删除删除 缺省缺省:NO:N

45、O四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司2 2、过大值删除门槛、过大值删除门槛(OUTLIER_THRESHOLD)(OUTLIER_THRESHOLD) 任何小于该值的拾取时间都不当作过大值任何小于该值的拾取时间都不当作过大值 缺省缺省:5ms:5ms3 3、每个站点的炮点、检波点时间深度分别计算、每个站点的炮点、检波点时间深度分别计算(STN_DEP)(STN_DEP) YES YES 分别计算分别计算 NO NO 同一站点的炮点、检波点时间深度同为它们的平均值同一站点的炮点、检波点时间深度同为它们的平均值 缺省缺省:YES:YES4 4、混合算法的构造项、混合算法的

46、构造项(SLACK_TERM)(SLACK_TERM) 构造项是偏移距构造项是偏移距/ /折射速度折射速度, ,而偏移距是确定的。本参数指定构而偏移距是确定的。本参数指定构 造项是速度还是中间站点一致的构造项造项是速度还是中间站点一致的构造项 VELOCITYS VELOCITYS 折射速度折射速度 STRUCTURAL STRUCTURAL 构造项构造项 缺省缺省: VELOCITYS: VELOCITYS 四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 五五 速度控制参数组速度控制参数组(VEL） 每个站点的地表速度必须指定，可用本参数组指定也可用井口控制参数每个站点的地表速

47、度必须指定，可用本参数组指定也可用井口控制参数组指定。若用本参数组指定，则以下表的形式：组指定。若用本参数组指定，则以下表的形式： LAYER SIN VEL INTERPOLATE LAYER 101 1400 NO参数参数1 1、层号（、层号（LAYERLAYER） 定义速度控制的层号定义速度控制的层号四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司2 2、站号（、站号（STNSTN） 速度控制点的站点号速度控制点的站点号 FIRST FIRST 第一个站点第一个站点 LAST LAST 最后一个站点最后一个站点 其它数组其它数组 ( (无缺省无缺省) )3 3、该点的层速度（

48、、该点的层速度（VELVEL）4 4、速度内插（、速度内插（INTERPOLATEINTERPOLATE） 从该点到下一个指定站点之间要不要速度内插从该点到下一个指定站点之间要不要速度内插 NO NO 不内插不内插 YES YES 要内插要内插 缺省：缺省：NONO四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 六六 时间时间-深度控制参数组深度控制参数组(TIME_DEPTHS） 用户提供时间用户提供时间- -深度。深度。若指定本参数组，则以下表的形式：若指定本参数组，则以下表的形式： LAYER SIN TIME_DEPTHL INTERPOLATE LAYER 101 30

49、1 NO参数参数1 1、层号（、层号（LAYERLAYER） 指定时间指定时间- -深度控制的层号深度控制的层号2 2、站号（、站号（STNSTN） 时间时间- -深度控制点的站点号深度控制点的站点号 FIRST FIRST 第一个站点第一个站点 LAST LAST 最后一个站点最后一个站点 其它数组其它数组 ( (无缺省无缺省) ) 四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司3 3、该点的时间、该点的时间- -深度（深度（TIME_DEPTHTIME_DEPTH）4 4、时间、时间- -深度内插（深度内插（INTERPOLATEINTERPOLATE） 从该点到下一个指定站

50、点之间要不要时间深度内插从该点到下一个指定站点之间要不要时间深度内插 NO NO 不内插不内插 YES YES 要内插要内插 缺省：缺省：NONO说明说明: :使用本参数组就意味着由用户给出指定折射层的时间深度，而不希望使用本参数组就意味着由用户给出指定折射层的时间深度，而不希望 模块去估算。如果要这样做就不应为该折射层提供偏移距范围。模块去估算。如果要这样做就不应为该折射层提供偏移距范围。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司 七七 折射高程控制参数组折射高程控制参数组(ELEV） 用户提供折射高程。若指定本参数组，则以下表的形式：用户提供折射高程。若指定本参数组，则以

51、下表的形式： LAYER ADJUSTMET SIN ELEV INTERPOLATE LAYER ADJUSTMET 101 111 NO参数参数1 1、层号（、层号（LAYERLAYER） 指定折射高程控制的层号指定折射高程控制的层号2 2、更改方法（、更改方法（ADJUSTMETADJUSTMET） TIME_DEPTH TIME_DEPTH 将高程更改为时间将高程更改为时间- -深度（速度无效）深度（速度无效） VELOCITY VELOCITY 将高程更改为速度（时间将高程更改为速度（时间- -深度无效）深度无效） 缺省：缺省： VELOCITYVELOCITY 该参数是指在绘图时将

52、高程值变换为时间深度或速度。该参数是指在绘图时将高程值变换为时间深度或速度。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司2 2、站号（、站号（STNSTN） 折射高程控制点的站点号折射高程控制点的站点号 FIRST FIRST 第一个站点第一个站点 LAST LAST 最后一个站点最后一个站点 其它数组其它数组 ( (无缺省无缺省) )3 3、该点的折射高程（、该点的折射高程（ELEVELEV）4 4、时间、时间- -深度内插（深度内插（INTERPOLATEINTERPOLATE） 从该点到下一个指定站点之间的折射高程要不要内插从该点到下一个指定站点之间的折射高程要不要内插

53、NO NO 不内插不内插 YES YES 要内插要内插 缺省：缺省：NONO说明：本参数组通常是在做了第一轮说明：本参数组通常是在做了第一轮EGRMEGRM之后，分析人员发现在那些外拓的之后，分析人员发现在那些外拓的 区域折射模型不合常理，这时用户就可迭代地用时间区域折射模型不合常理，这时用户就可迭代地用时间- -深度控制参数组深度控制参数组 和折射高程控制参数组给出一个指定的外拓。和折射高程控制参数组给出一个指定的外拓。四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司八八 速度平滑参数组（速度平滑参数组（VEL_SMOOTHINGVEL_SMOOTHING）例：例：LAYER S

54、IN SMOOTH_DIAMETER INTERPOLATELAYER 100 55 NO参数参数1 1、层号（、层号（LAYERLAYER） 指定折射高程控制的层号指定折射高程控制的层号2 2、站号（、站号（STNSTN） 折射高程控制点的站点号折射高程控制点的站点号 FIRST FIRST 第一个站点第一个站点 LAST LAST 最后一个站点最后一个站点 其它数组其它数组 ( (无缺省无缺省) )3 3、平滑直径（、平滑直径（ SMOOTH_DIAMETERSMOOTH_DIAMETER） 该站点速度的平滑直径，以站点个数表示。该站点速度的平滑直径，以站点个数表示。 缺省：缺省：0 0

55、（不平滑）（不平滑）四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司4 4、平滑直径内插（、平滑直径内插（INTERPOLATEINTERPOLATE） 从该点到下一个指定站点之间的平滑直径要不要内插从该点到下一个指定站点之间的平滑直径要不要内插 NO NO 不内插不内插 （缺省）（缺省） YES YES 要内插要内插四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司9 9时间时间- -深度平滑参数组（深度平滑参数组（TIME_SMOOTHINGTIME_SMOOTHING）十十 中间基准面平滑参数组（中间基准面平滑参数组（DATUM_SMOOTHINGDATUM_SMO

56、OTHING）十一十一 折射高程平滑参数组（折射高程平滑参数组（ELEV_SMOOTHINGELEV_SMOOTHING） 以上三个平滑参数组的格式与速度平滑参数组相同。以上三个平滑参数组的格式与速度平滑参数组相同。 时间时间- -深度平滑参数组平滑折射层的时间深度，由于它相当于平滑所计深度平滑参数组平滑折射层的时间深度，由于它相当于平滑所计算的静校正量，所以要小心使用，大的平滑直径将影响长波长和短波长静算的静校正量，所以要小心使用，大的平滑直径将影响长波长和短波长静校正量的解。校正量的解。 中间基准面平滑参数组用以产生中间基准面，各站点最低折射深度的平中间基准面平滑参数组用以产生中间基准面，

57、各站点最低折射深度的平滑结果就是中间基准面。滑结果就是中间基准面。 当折射高程出现明显的异常，且这种异常很可能是速度引起的时，可用当折射高程出现明显的异常，且这种异常很可能是速度引起的时，可用折射高程平滑参数组平滑所计算的折射高程。折射高程的平滑将使该折射折射高程平滑参数组平滑所计算的折射高程。折射高程的平滑将使该折射层以上的速度得到修正。（折射高程显的异常根据质量控制图判断）层以上的速度得到修正。（折射高程显的异常根据质量控制图判断）四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司十二十二 浮动中间基准面参数组（浮动中间基准面参数组（INTERMEDIATE_DATUM） 如前所

58、述如前所述, , ,EGRMEGRM模块根据折射深度自动生成中间基准面。用户可用本参数组修模块根据折射深度自动生成中间基准面。用户可用本参数组修改这个这个基准面或重新指定一个平滑中间基准面。改这个这个基准面或重新指定一个平滑中间基准面。 修改中间基准面就是将该中间基准面做上下平移，所用参数是：修改中间基准面就是将该中间基准面做上下平移，所用参数是：平滑中间基准面平移（平滑中间基准面平移（DATUM_SHIFTDATUM_SHIFT）缺省：缺省：0.00.0 重新指定平滑中间基准面则用以下参数表重新指定平滑中间基准面则用以下参数表: : STN ELEV 60 11001 1、站点号（、站点号（

59、STNSTN） 折射高程控制点的站点号折射高程控制点的站点号 FIRST FIRST 第一个站点第一个站点 LAST LAST 最后一个站点最后一个站点 其它数组其它数组 ( (无缺省无缺省) )2 2、该站点的基准面高程（、该站点的基准面高程（ELEVELEV） 无缺省无缺省四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司13井口控制参数组（井口控制参数组（UPHOLE） 该参数组将井口信息用于模块的计算中。该参数组将井口信息用于模块的计算中。 参数参数1、井口时间应用方向（、井口时间应用方向（DIRECTION） 本参数指定井口时间应用方向是沿倾斜路径或是垂直路径。本参数指定井

60、口时间应用方向是沿倾斜路径或是垂直路径。 RAYPATH 沿倾斜路径沿倾斜路径 VERTICAL 沿垂直路径沿垂直路径 缺省：缺省： 垂直（井较浅时）垂直（井较浅时）BCAD井底井底四川石油管理局地球物理勘探公司四川石油管理局地球物理勘探公司2 2、虚拟折射层的应用（、虚拟折射层的应用（PSEUDO_REFRACTOR_APPPSEUDO_REFRACTOR_APP） 本参数指定如何应用井口信息本参数指定如何应用井口信息 VELOCITIES VELOCITIES 井口信息用于修改用户指定的近地表速度场井口信息用于修改用户指定的近地表速度场 ELEVATIONS ELEVATIONS 将井口深