2. 中国石油东方地球物理公司研究院, 河北涿州072750
2. Research Institute of BGP, CNPC, Zhuozhou 072750, China
中国西部地区地表起伏剧烈, 低降速带巨厚, 高速地层和低速覆盖层分布无规律, 地震波射线的入射角和出射角变化范围大, 静校正会造成地震波射线在走时计算时出现误差, 进而造成波场失真[1-2]。
常规时间域叠前道集的静校正处理是将静校正和剩余静校正的总量在CMP道集域进行高频、低频分离, 分离出的低频分量是以CMP为单元统计出的时间域地形平均面, 即CMP参考面。早期的叠前深度偏移处理通常通过替换速度将时间域的CMP参考面换算到深度域, 作为深度偏移基准面。复杂山地地表高程变化大, 造成基于真实地表与CMP参考面计算得到的深度偏移基准面差异大, 因而射线路径的走时误差大, 波场畸变严重。为了解决波场畸变问题, 目前深度偏移成像时通常采用对地表高程进行小尺度平滑, 选取地表圆滑面作为深度域偏移基准面的处理方法。地表圆滑面反映地表高程的变化形态, 应尽可能减少其与地表的高差。该处理思路是为了消除地表圆滑面与CMP参考面之间的低频时差, 将激发点和接收点的真实位置垂直校正到地表圆滑面上, 对偏移前的道集数据应用高频静校正, 以减小真实激发点、接收点与地表圆滑面之间的射线旅行时误差。
以上方法技术虽然削弱了道集时差对偏移成像的影响, 但仍然存在以下问题:①虽然地震数据中引入了高频静校时差, 但道集时差对地震波场走时的影响并未完全消除; ②深度偏移基准面是根据时间域基准面静校正量求得的, 它与后续建立的深度域近地表速度模型不一定匹配。
1.2 浅层速度建模问题分析叠前深度偏移的核心是获得准确的速度模型。中、深层速度模型的准确程度不仅取决于道集数据信噪比, 还依赖于浅表层速度模型的精度。浅层速度误差会传递至下伏地层, 进而影响深层速度模型的精度, 因此建立准确的浅层速度模型对深度偏移成像至关重要。利用反射波进行速度层析反演, 要求反射波具有高信噪比、大排列长度、可准确拾取共反射点(common reflection point, CRP)道集的剩余延迟(residual moveout, RMO)。但因三维地震浅层资料存在覆盖道数少、信噪比低的缺点, 使得共成像点道集浅层的反射波同相轴很难被识别, 剩余延迟拾取不准确。因此用于评价偏移速度正确与否的判别准则(同相轴校平)失去其依据, 利用反射波层析成像速度优化方法无法求取准确的近地表及浅层地层速度。图 1是我国西部某复杂山地的三维地震速度谱及CRP道集, 可以看出中深层的反射信号具有较高信噪比以及较大排列长度, 浅层地震资料信噪比低、覆盖道数少, 难以识别有效反射信息, 因此利用反射波层析建立准确的浅层速度模型存在困难。
利用初至波层析反演出高精度的浅表层速度模型, 并与偏移基准面配套应用, 可建立合理的深度偏移浅层速度模型, 实现真地表叠前深度偏移。本文研究并提出了面向复杂山地低信噪比的真地表浅层速度建模技术, 核心技术思路包括如下内容:
1) 在微测井[3]、钻测井等速度信息约束的基础上, 利用回折波层析反演[4-7]方法建立近地表及浅层速度模型, 以提高浅层速度模型的精度。
2) 静校正直接面向深度域成像, 即将时间域大炮初至数据的回折波反演得到的近地表模型直接作为深度域浅层速度模型, 将时间域静校正仅作为叠前去噪等处理方法的辅助手段, 因此需要在叠前深度偏移前对应用了基准面静校正量的地震数据进行静校正“反应用”, 时间域静校正引起的时差问题可采用真地表偏移面+浅层建模+深度偏移的处理方法解决, 根据深度偏移成像质量决定静校正方法及优选参数。
3) 在表层速度准确的情况下, 选取包含高精度旅行时算法的层析成像和偏移方法, 完成整体速度模型的建立以实现真地表叠前深度偏移。
2.2 微测井约束回折波分层层析反演复杂山地采集的浅层地震资料信噪比低、覆盖次数少、缺乏远炮检距波场信息, 反射波层析成像无法有效求取准确的近地表速度, 因此大炮初至数据的回折波层析反演技术[8-11]成为时间域静校正处理和深度偏移近地表速度建模的有效技术手段。
回折波层析反演考虑了地震波在近地表介质中传播的各种现象, 对回折波或连续折射直达波通过高精度的交互反演得到近地表介质速度变化情况, 从而建立复杂的近地表模型。但受采集资料的限制, 当最小炮检距、低降速层厚度较小时, 利用初至波信息无法准确刻画表层极低速度的信息, 需要联合微测井和近炮检距初至波信息反演低速层速度, 再以反演速度为约束, 联合微测井和中、远炮检距初至波信息反演高精度的降速层速度, 从而建立精确的近地表速度模型。
大炮初至数据的回折波层析反演技术对近地表结构进行高密度速度单元划分, 将回折波的走时描述为对介质慢度函数沿射线路径的线积分, 回折波射线追踪利用初始速度模型和炮检点位置计算旅行时, 初至时间与正演旅行时之差可被反演成射线经过的每个速度单元慢度的扰动。
为了使反演结果与实际模型更吻合并使反演结果稳定合理, 本文采用了微测井约束回折波层析反演方法。该方法的目标函数包含了速度约束项, 先将微测井解释的近地表模型作为约束条件[11], 再与旅行时残差共同构成目标函数, 最后求解拉格朗日约束条件下的目标函数最小二乘解。在层析反演中, 由M条射线和N个未知数建立的层析矩阵方程为[12-15]:
$ \boldsymbol{A}·Δ \boldsymbol{S}=Δ \boldsymbol{T} $ | (1) |
式中:A为射线路径集合的Jacobi矩阵; ΔS为慢度修正向量; ΔT为旅行时残差矩阵。可分别表示为:
$ \begin{array}{c} \boldsymbol{A}=\left[\begin{array}{cccc} \frac{\partial t_{1}}{\partial S_{1}} & \frac{\partial t_{1}}{\partial S_{2}} & \cdots & \frac{\partial t_{1}}{\partial S_{N}} \\ \frac{\partial t_{2}}{\partial S_{1}} & \frac{\partial t_{2}}{\partial S_{2}} & \cdots & \frac{\partial t_{2}}{\partial S_{N}} \\ \vdots & \vdots & & \vdots \\ \frac{\partial t_{M}}{\partial S_{1}} & \frac{\partial t_{M}}{\partial S_{2}} & \cdots & \frac{\partial t_{M}}{\partial S_{N}} \end{array}\right]\\ \Delta \boldsymbol{S}=\left[\begin{array}{c} \Delta S_{1} \\ \Delta S_{2} \\ \vdots \\ \Delta S_{N} \end{array}\right] \quad \Delta \boldsymbol{T}=\left[\begin{array}{c} \Delta t_{1} \\ \Delta t_{2} \\ \vdots \\ \Delta t_{N} \end{array}\right] \end{array} $ | (2) |
假设通过微测井获得了L个准确的速度信息, 即方程(1)中的N个速度单元中有L个已知量, 用L维方程可表示为:
$ \boldsymbol{C}·Δ \boldsymbol{S}= \boldsymbol{F} $ | (3) |
式中:C为已知的L个约束条件的Jacobi矩阵; F为约束旅行时残差矩阵。
其中,
$ \begin{array}{c} \boldsymbol{C}=\left[\begin{array}{cccc} C_{11} & C_{12} & \cdots & C_{1 N} \\ C_{21} & C_{22} & \cdots & C_{2 N} \\ \cdots & \cdots & & \cdots \\ C_{L 1} & C_{L 2} & \cdots & C_{L N} \end{array}\right]\\ \Delta \boldsymbol{S}=\left[\begin{array}{c} \Delta S_{1} \\ \Delta S_{2} \\ \vdots \\ \Delta S_{N} \end{array}\right] \quad \Delta \boldsymbol{F}=\left[\begin{array}{c} f_{1} \\ f_{2} \\ \vdots \\ f_{L} \end{array}\right] \end{array} $ | (4) |
在初至波层析中, 假如
$ C_{i j}=\left\{\begin{array}{ll} 1 & i=j \\ 0 & i \neq j \end{array} \quad f_{i}=\tilde{S}_{i}-S_{i}\right. $ | (5) |
式中:Si为第i个单元上一次迭代层析反演慢度值;
联合方程(1)和方程(3), 构建成的方程可表示为:
$ \boldsymbol{B}·Δ \boldsymbol{S}= \boldsymbol{H} $ | (6) |
其中,
$ \boldsymbol{B}=\left[\begin{array}{c} \boldsymbol{A} \\ \sqrt{\lambda} \boldsymbol{C} \end{array}\right] \quad \boldsymbol{H}=\left[\begin{array}{c} \Delta \boldsymbol{T} \\ \sqrt{\lambda} \boldsymbol{F} \end{array}\right] $ | (7) |
式中:B是(M+L)×N维矩阵; λ为约束权重系数。
拉格朗日最优约束的目标函数可表示为:
$ \begin{array}{c} ξ(Δ \boldsymbol{S})=( \boldsymbol{B}Δ \boldsymbol{S}- \boldsymbol{H})^{\rm T}( \boldsymbol{B}Δ \boldsymbol{S}- \boldsymbol{H})\\ \;\;\;\;\;\;\;\;\;\;\;\;\;=( \boldsymbol{A}Δ \boldsymbol{S}-Δ \boldsymbol{T})^{\rm T}( \boldsymbol{A}Δ \boldsymbol{S}-Δ \boldsymbol{T})+\\ \;\;\;\;\;\;\;\;λ( \boldsymbol{C}Δ \boldsymbol{S}- \boldsymbol{F})^{\rm T}( \boldsymbol{C}Δ \boldsymbol{S}- \boldsymbol{F}) \end{array} $ | (8) |
式中:ξ(ΔS)为总目标函数; (AΔS-ΔT)T(AΔS-ΔT)为观测数据目标函数; λ(CΔS-F)T(CΔS-F)为先验约束条件的目标函数。当目标函数最小时, 即当▽ξ(ΔS)=0, 利用(8)式计算得到的ΔS为最优解。如果输入部分将已知的真实速度作为约束条件, 利用方程(6)可以获得更准确的解。从目标函数中可看出约束权重系数λ=0时, 相当于没有约束, λ越大, 相当于约束方程在总目标函数中所占权重越大, 即约束越强。
利用迭代层析反演中层析反演速度与初始模型速度的比值R, 可以确定系数λ值:
$ R= \frac{{V_{1}}}{{V_{m}}}\\ λ=1- \frac{{R-R_{\min}}}{{R_{\max}-R_{\min}}} $ | (9) |
式中:V1是初始模型速度与层析迭代反演速度中的较小者; Vm是两者中的较大者; Rmin和Rmax分别为R的最小值和最大值。由(9)式可知R越大, 则λ越小; R越小, 则λ值越大。根据(9)式, 利用迭代过程中的反演速度与初始模型速度差, 可实现约束权重系数的自适应计算。
在上述计算过程中, 可通过自适应算法确定λ值。对于极浅层速度模型, 由于缺少小炮检距地震数据, 微测井数据的λ可赋值为1。随着深度的增加, 微测井约束的λ值逐渐下降, 在接近高速层顶面时, 初始模型速度的贡献越来越少, 此时利用大炮初至数据就能充分解决微测井资料空间采样不均匀和探测深度不足的问题, 此时应赋予大炮初至数据较大的λ值。
在实际应用中, 要对收集的微测井、小折射、钻测井等资料的可靠性进行认真分析和筛选[16-17], 确保近地表调查资料的准确性, 从而提高近地表速度模型的反演精度。图 2是微测井约束前、后反演的浅层速度模型, 可以看到采用微测井约束后, 浅层速度模型与实际结果更吻合。
初始速度模型的建立对最终速度模型的质量以及速度模型的迭代修正效率起着至关重要的作用, 由于叠前深度偏移的速度模型是逐层反演的, 因此在初始模型建立的过程中, 浅层速度模型的精度会严重影响中、深层速度模型迭代修正过程。根据复杂山地地区的资料特点, 初始速度模型建模时, 将近地表速度模型与中、深层速度模型分开建立, 然后再融合, 具体思路如下:
1) 利用微测井约束回折波分层层析反演求取高精度近地表速度模型。
2) 中深层的速度模型的建立:①利用时间域均方根速度场, 通过约束速度反演得到时间域层速度; ②在构造模型约束条件下, 通过时深转换获取深度域层速度; ③利用VSP等钻、测井资料对速度场进行修正得到中、深层速度模型。
3) 在近地表速度模型中选取一个平滑的高速层顶界面作为速度融合界面, 在设定的时窗内将近地表速度镶嵌到中、深层速度模型中, 从而获得用于深度偏移的初始速度模型, 在后续速度模型迭代修正过程中, 保持近地表速度模型不变, 仅对高速层顶界面以下的地层进行速度迭代更新。
解决起伏地表的波场畸变问题, 不仅需要高精度的近地表速度模型, 还必须将偏移基准面选为真实地表面, 从而实现真地表叠前深度偏移。为此需要正确处理时间域静校正与偏移基准面的关系。既然深度域浅层速度模型与时间域静校正是一致的, 即已经在深度域利用近地表速度解决了时间域静校正问题, 那么应用了静校正的时间域偏移前道集数据应该在深度域偏移前对应用的静校正展开“反应用”, 使深度偏移前的数据无需进行基准面静校正。
对柴达木盆地西部某工区回折波进行层析反演迭代, 得到的浅层速度模型及叠前深度偏移剖面如图 3所示, 本工区共进行了4轮迭代, 可以看出随着迭代次数的增加, 速度场细节逐渐清晰, 最终得到了精确的浅层速度模型, 叠前深度偏移剖面成像质量也得到了大幅改善。对该工区应用微测井等资料约束前、后的回折波层析反演结果如图 4所示, 图 4c和图 4d中蓝线代表VSP测井速度, 红线代表层析成像反演的层速度。可以看出, 经微测井资料约束后, 回折波层析反演的表层速度横向刻画更精细, 整体速度与VSP测井速度在浅、中层一致性更好, 浅、中层速度很好地反映出逆冲断层速度上冲的趋势, 符合目前的地质认识, 叠前深度偏移剖面上反射波同相轴连续性明显增强, 断面波特征更清楚, 断层下部成像质量也得到明显改善。
目前制约复杂山地真地表叠前深度偏移实际应用的因素主要包括高精度浅表层速度建模和高精度偏移算法。提高旅行时计算精度, 使旅行时计算结果能够描述真实近地表地层速度横向变化特点及规律, 即让深度偏移解决时间域的静校正问题, 有利于深度域高精度浅层速度模型建立。在建立浅层速度模型的基础上, 实现真地表叠前深度偏移还需要选取合适的层析成像和偏移方法, 但无论如何, 旅行时计算都是重要步骤之一, 在条件允许的情况下, 须尽可能采用波前构建法和频率域波动方程数值解法进行旅行时计算。建立高精度的整体速度模型后, 可采用波场外推的偏移方法实现真地表叠前深度偏移。
采用多方位各向异性层析成像或者更为先进的绕射层析成像方法可提高速度建模精度。通常采用由浅及深、逐步递进的建模方法, 在获得高精度的浅层速度模型后, 采用基于井控地质导向约束的TTI网格层析成像迭代技术完成整体速度建模。该技术是在数字化露头、地层倾角测井、三维地震构造解释模型等多信息的约束下, 通过联合反演减少层析解的不确定性, 并确保速度模型在符合地质及构造认识规律的条件下迭代更新[18-19]。实际资料处理结果表明, 网格尺度影响走时计算精度, 进而对速度模型的细节刻画造成影响, 因此在层析反演过程中, 可通过逐步细分层析反演的网格, 提高反演结果的分辨率。
3 应用效果分析图 5和图 6是西部某工区近地表圆滑面叠前深度偏移剖面和真地表叠前深度偏移剖面, 该区开展过多轮次的处理攻关, 以往速度建模都是在近地表圆滑面上开展的, 存在断层刻画不清楚, 构造细节未落实, 部分地区井震深度误差大的问题。此次采用了真地表叠前深度偏移技术, 首先通过微测井约束初至回转波分层层析反演建立深度0~1500 m的速度模型, 然后采用基于井控地质导向约束的多方位网格层析[20]成像及TTI网格层析成像建立并优化深度大于1 500 m的速度模型, 最后利用TTI真地表逆时叠前深度偏移实现高精度偏移成像。由图 5可见, 真地表叠前深度偏移剖面上逆冲断层下盘盐内断背斜构造形态清晰可靠, 成像效果较近地表圆滑面的叠前深度偏移成果明显改善, 主要表现在以下方面:
1) 真地表叠前深度偏移剖面, 反射波特征突出, 波组特征清晰, 偏移归位准确, 断点、断面位置更清楚, 有利于精细构造解释。
2) 真地表叠前深度偏移较好地恢复了地下界面的真实形态。
如图 6所示, 与以往叠前深度偏移的结果相比, 真地表叠前深度偏移得到的剖面构造形态变化明显, 利用本工区已经完钻的两口井资料验证可知, 真地表叠前深度偏移的结果与井吻合程度高, 证明真地表叠前深度偏移的结果真实可靠。
4 结论与认识总结真地表浅层速度建模的研究及其实际应用成果得出以下结论:
1) 真地表叠前深度偏移浅层建模技术在很大程度上解决了近地表圆滑偏移面与时间域静校正引起的波场畸变问题。
2) 微测井约束回折波分层层析反演是一种可靠地建立高精度浅层速度模型的方法, 在多种信息的约束下, 只有在获得高精度浅层速度模型的基础上, 才能够建立准确的中、深层模型。
3) 复杂山地的真地表叠前深度偏移成像精度受多种因素的影响, 只有从旅行时算法、层析成像、偏移算法等多方面深入开展理论研究与实际资料试验, 提高其中每一个环节的计算精度, 才能够解决山前带的复杂构造成像问题。
[1] |
刘定进, 蒋波, 李博, 等. 起伏地表逆时偏移在复杂山前带地震成像中的应用[J]. 石油地球物理勘探, 2016, 51(2): 315-323. LIU D J, JIANG B, LI B, et al. Application of reverse time migration of undulated surface in seismic imaging of complex piedmont area[J]. Oil Geophysical Prospecting, 2016, 51(2): 315-323. |
[2] |
王东, 李文卉, 熊登, 等. 静校正对波动方程叠前深度偏移成像效果的影响[J]. 石油地球物理勘探, 2016, 51(2): 76-80. WANG D, LI W H, XIONG D, et al. Aeffect of static correction on wave equation prestack depth migration imaging[J]. Oil Geophysical Prospecting, 2016, 51(2): 76-80. |
[3] |
李明海. 地面微测井在山地表层结构调查中的应用[J]. 勘探地球物理进展, 2008, 31(5): 378-382. LI M H. Application of uphole methods in the investigation of near surface structurein Mountainous Areas[J]. Progress in Exploration Geophysics, 2008, 31(5): 378-382. |
[4] |
李录明, 罗省贤. 复杂三维表层模型层析反演和静校正[J]. 石油地球物理勘探, 2003, 38(6): 636-642. LI L M, LUO X X. Tomographic inversion for 3D complex near-surface model and static correction[J]. Oil Geophysical Prospecting, 2003, 38(6): 636-642. DOI:10.3321/j.issn:1000-7210.2003.06.010 |
[5] |
张凯, 李振春. 双复杂条件下层析速度反演方法研究[J]. 地球物理学进展, 2013, 28(6): 3001-3006. ZHANG K, LI Z C. Tomographic velocity inversion method with dual-complexity[J]. Progress in Geophysics, 2013, 28(6): 3001-3006. |
[6] |
李录明, 罗省贤, 赵波, 等. 初至波表层模型层析反演[J]. 石油地球物理勘探, 2000, 35(5): 559-564. LI L M, LUO X X, ZHAO B, et al. Tomographic inversion of first break in surface model[J]. Oil Geophysical Prospecting, 2000, 35(5): 559-564. DOI:10.3321/j.issn:1000-7210.2000.05.002 |
[7] |
季天愚, 李澈, 李景业, 等. 无射线追踪层析静校正技术及其应用研究[J]. 物探化探计算技术, 2014, 36(4): 475-480. JI T Y, LI C, LI J Y, et al. Non-ray tracing tomography static correction technology and its application[J]. Computing Techniques for Geophysical and Geochemical Exploration, 2014, 36(4): 475-480. DOI:10.3969/j.issn.1001-1749.2014.04.18 |
[8] |
刘小民, 邬达理, 梁硕博, 等. 潜水波胖射线走时层析速度反演及其在深度偏移速度建模中的应用[J]. 石油物探, 2017, 56(5): 718-726. LIU X M, WU D L, LIANG S B, et al. Diving wave tomography velocity inversion using fat ray in prestack depth migration[J]. Geophysical Prospecting for Petroleum, 2017, 56(5): 718-726. DOI:10.3969/j.issn.1000-1441.2017.05.012 |
[9] |
郎玉泉. 初至波层析反演表层速度的应用实践[J]. 中国煤炭地质, 2018, 30(3): 76-81. LANG Y Q. Application practices of preliminary wave tomographic inversion near-surface velocity[J]. Coal Geology of China, 2018, 30(3): 76-81. DOI:10.3969/j.issn.1674-1803.2018.03.13 |
[10] |
胡自多, 贺振华, 王西文, 等. 无射线追踪层析静校正技术在黄土塬区的应用[J]. 石油地球物理勘探, 2010, 45(1): 94-100. HU Z D, HE Z H, WANG X W, et al. Application of non-ray tracing tomography static correction technology in loess plateau area[J]. Oil Geophysical Prospecting, 2010, 45(1): 94-100. |
[11] |
袁刚, 冯心远, 蒋波, 等. 约束层析反演及其在地震速度计算中的应用[J]. 石油物探, 2013, 52(1): 55-59. YUAN G, FENG X Y, JIANG B, et al. Constrained tomographic inversion and its application in seismic velocity calculation[J]. Geophysical Prospecting for Petroleum, 2013, 52(1): 55-59. DOI:10.3969/j.issn.1000-1441.2013.01.009 |
[12] |
王孝, 曾华会, 刘文卿, 等. 基于微测井分步约束的近地表速度层析反演[J]. 石油地球物理勘探, 2018, 53(增刊1): 69-74. WANG X, ZENG H H, LIU W Q, et al. Tomographic inversion of near surface velocity based on the substep restriction of microlog[J]. Oil Geophysical Prospecting, 2018, 53(S1): 69-74. |
[13] |
高利东, 刘连升, 张志禹. 基于约束初至波层析反演的近地表建模技术[J]. 石油化工应用, 2013, 32(6): 31-33. GAO L D, LIU L S, ZHANG Z Y. Constrained first arrival tomography inversion for near surface model building[J]. Petrochemical Industry Application, 2013, 32(6): 31-33. DOI:10.3969/j.issn.1673-5285.2013.06.010 |
[14] |
杨国权, 丁鹏程, 李振春, 等. 应用线性程函方程和整形正则化的三维初至波旅行时层析[J]. 石油地球物理勘探, 2017, 52(2): 264-272. YANG G Q, DING P C, LI Z C, et al. Travel time tomography of three-dimensional first break using linear equation and shaping regularization[J]. Oil Geophysical Prospecting, 2017, 52(2): 264-272. |
[15] |
赵玲芝, 谷跃民, 张建中. 多信息融合的近地表速度建模技术及应用[J]. 石油地球物理勘探, 2017, 52(1): 34-41. ZHAO L Z, GU Y M, ZHANG J Z. Multi-information fusion near-surface velocity modeling technology and its application[J]. Oil Geophysical Prospecting, 2017, 52(1): 34-41. |
[16] |
王振宇, 杨勤勇, 李振春, 等. 近地表速度建模研究现状及发展趋势[J]. 地球科学进展, 2014, 29(10): 1138-1148. WANG Z Y, YANG Q Y, LI Z C, et al. Research status and development trend of near-surface velocity modeling[J]. Advances in Earth Science, 2014, 29(10): 1138-1148. DOI:10.11867/j.issn.1001-8166.2014.10.1138 |
[17] |
赵玲芝, 冯许魁, 祖云飞, 等. 复杂地表区叠前深度偏移近地表速度应用技术研究[J]. 地球科学进展, 2018, 33(3): 1176-1183. ZHAO L Z, FENG X K, ZU Y F, et al. Applied research of the complex near-surface velocity in the depth migration[J]. Advances in Earth Science, 2018, 33(3): 1176-1183. |
[18] |
罗宾.地震资料叠前偏移成像——方法、原理和优缺点分析[M].王克斌, 译.北京: 石油工业出版社, 2012: 61-71 ROBEIN E.Seismic image—A review of the techniques, their principles, merits and limitations[M].WANG K B translators.Beijing: Petroleum Industry Press, 2012: 61-71 |
[19] |
陈云峰, 杨淑卿. 利用层析成像方法反演近地表速度模型[J]. 海洋地质前沿, 2014, 30(5): 47-51. CHENG Y F, YANG S Q. Building of near surface velocity model using tomography inversion method[J]. Marine Geology Frontiers, 2014, 30(5): 47-51. |
[20] |
方勇, 温铁民, 李虹, 等. 多方位网格层析成像技术及应用效果[J]. 物探化探计算技术, 2016, 38(5): 677-680. FANG Y, WEN T M, LI H, et al. Multidirectional grid tomography technology and its application effect[J]. Computing Techniques for Geophysical and Geochemical Exploration, 2016, 38(5): 677-680. DOI:10.3969/j.issn.1001-1749.2016.05.17 |