CN111077575A - 一种深度域速度建模方法与装置 - Google Patents
一种深度域速度建模方法与装置 Download PDFInfo
- Publication number
- CN111077575A CN111077575A CN201811217379.6A CN201811217379A CN111077575A CN 111077575 A CN111077575 A CN 111077575A CN 201811217379 A CN201811217379 A CN 201811217379A CN 111077575 A CN111077575 A CN 111077575A
- Authority
- CN
- China
- Prior art keywords
- depth
- velocity
- domain
- seismic
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/308—Time lapse or 4D effects, e.g. production related effects to the formation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
Abstract
本发明涉及一种深度域速度建模方法与装置,属于地震数据处理技术领域,本发明将井资料和地震层位信息用于深度域速度模型的整个过程中,利用三维叠前时间偏移速度得到深度域层速度模型,将声波测井数据代替工区范围内有井区域的地震速度,作为有井区域的速度模型,在地震解释层位的约束下,井速度和地震速度进行融合得到无井区域的速度模型,最后将有井区域的速度模型与无井区域的速度模型合并,得到初始的深度域速度模型,对初始的深度域速度模型进行网格层析反演优化,得到优化的深度域速度模型。本发明提高了深度域速度模型的精度,化简了网格层析反演的迭代次数,减小了进行叠前深度偏移的井震误差,满足勘探需求。
Description
技术领域
本发明属于地震数据处理技术领域,具体涉及一种深度域速度建模方法与装置。
背景技术
叠前深度偏移对速度模型的精度要求极高,速度不准确会造成复杂地区的偏移结果存在较大的井震误差,严重影响井位部署的成功率。目前,常规的深度速度建模是以时间偏移速度为初始模型,进行基于偏移道集剩余曲率的网格层析反演。这种仅靠地震数据的建模方法,只能反演出平滑、低频的背景速度模型,对复杂构造和特殊地质体,难以反演出其构造形态和细节特征。
近年来,利用声波测井速度和地震层位约束提高深度域速度模型精度的研究持续进行,取得了一些成果。例如,公告号为CN106597533B的中国专利提出了“一种用于山前带地震资料处理的深度域速度建模方法”,该方法利用速度约束反演对叠前时间偏移速度进行DIX公式转换为深度域层速度体,得到初始层速度模型,再进行网格层析反演优化速度模型。但该方法只在网格层析反演中利用井约束进行优化,并没有将井资料和地震层位信息用于初始速度建模以及应用于整个速度建模过程中,造成建立的深度域速度模型的精度不高,增加了网格层析反演的迭代次数,以致利用该深度域速度模型进行叠前深度偏移存在较大的井震误差。
发明内容
本发明的目的是提供一种深度域速度建模方法与装置,用于解决采用现有深度域速度建模方法建立深度域速度模型精度低、增加了网格层析反演的迭代次数,以致进行叠前深度偏移存在较大的井震误差的问题。
为解决上述技术问题,本发明提出一种深度域速度建模方法,包括以下步骤:
1)获取工区范围内的地震数据、声波测井数据和地震解释层位数据,其中,地震数据包括三维叠前时间偏移速度,声波测井数据包括深度、与深度对应的声波速度,地震解释层位数据为设定区域的深度解释层位;
2)利用三维叠前时间偏移速度进行建模,得到初始时间域的速度模型,对初始时间域的速度模型进行两次DIX转换,得到深度域层速度模型;
3)将声波测井数据代替深度域层速度模型中有井区域的地震速度,得到有井区域的速度模型;对于深度域层速度模型中的无井区域,按照滑动加权计算方法,采用克里金差值方法沿着地震层位对邻近区域的井速度和地震速度进行不同权重的插值计算得到速度值,得到无井区域的速度模型;将有井区域的速度模型与无井区域的速度模型合并,得到初始深度域速度模型;
4)利用地震解释层位数据,对初始深度域速度模型进行网格层析反演优化,得到优化的深度域速度模型。
本发明将井资料和地震层位信息用于深度域速度模型的整个过程中,利用地震数据包括的三维叠前时间偏移速度得到深度域层速度模型,将声波测井数据代替工区范围内有井区域的地震速度,作为有井区域的速度模型,在地震解释层位的约束下,井速度和地震速度进行融合得到无井区域的速度模型,最后将有井区域的速度模型与无井区域的速度模型合并,得到初始的深度域速度模型,对初始的深度域速度模型进行网格层析反演优化,得到优化的深度域速度模型。本发明的深度域速度建模方法提高了深度域速度模型的精度,化简了网格层析反演的迭代次数,减小了进行叠前深度偏移的井震误差,满足勘探需求。
为了解决声波测井数据中井深度采样不均匀的问题,步骤3)中在声波测井数据代替深度域层速度模型中有井区域的地震速度之前,还包括:根据目标层纵向分辨率的要求,对声波测井数据进行重采样,且采样点之间深度间隔一致,得到重采样处理后的声波测井数据。
为了去除声波测井数据的噪声,步骤3)中在声波测井数据代替深度域层速度模型中有井区域的地震速度之前,还包括:采用滑动平均数字滤波法去除声波测井数据的噪声,得到去噪处理后的声波测井数据。
为了实现声波测井数据对深度域层速度模型中有井区域的地震速度的代替,得到深度域层速度模型之后,还包括:根据地震解释层位数据将深度域层速度模型在三维空间划分成不同深度区域的深度域层速度模型,根据划分的各个深度区域,将声波测井数据代替相应深度区域的深度域层速度模型中有井区域的地震速度。
为解决上述技术问题,本发明还提出一种深度域速度建模装置,包括存储器和处理器,以及存储在存储器上并在处理器上运行的计算机程序,处理器与存储器相耦合,处理器执行计算机程序时实现以下步骤:
1)获取工区范围内的地震数据、声波测井数据和地震解释层位数据,其中,地震数据包括三维叠前时间偏移速度,声波测井数据包括深度、与深度对应的声波速度,地震解释层位数据为设定区域的深度解释层位;
2)利用三维叠前时间偏移速度进行建模,得到初始时间域的速度模型,对初始时间域的速度模型进行DIX转换,得到深度域层速度模型;
3)将声波测井数据代替深度域层速度模型中有井区域的地震速度,得到有井区域的速度模型;对于深度域层速度模型中的无井区域,按照滑动加权计算方法,采用克里金差值方法沿着地震层位对邻近区域的井速度和地震速度进行不同权重的插值计算得到速度值,得到无井区域的速度模型;将有井区域的速度模型与无井区域的速度模型合并,得到初始深度域速度模型;
4)利用地震解释层位数据,对初始深度域速度模型进行网格层析反演优化,得到优化的深度域速度模型。
本发明将井资料和地震层位信息用于深度域速度模型的整个过程中,利用地震数据包括的三维叠前时间偏移速度得到深度域层速度模型,将声波测井数据代替工区范围内有井区域的地震速度,作为有井区域的速度模型,在地震解释层位的约束下,井速度和地震速度进行融合得到无井区域的速度模型,最后将有井区域的速度模型与无井区域的速度模型合并,得到初始的深度域速度模型,对初始的深度域速度模型进行网格层析反演优化,得到优化的深度域速度模型。本发明的深度域速度建模装置提高了深度域速度模型的精度,化简了网格层析反演的迭代次数,减小了进行叠前深度偏移的井震误差,满足勘探需求。
为了解决声波测井数据中井深度采样不均匀的问题,步骤3)中在声波测井数据代替深度域层速度模型中有井区域的地震速度之前,还包括:根据目标层纵向分辨率的要求,对声波测井数据进行重采样,且采样点之间深度间隔一致,得到重采样处理后的声波测井数据。
为了去除声波测井数据的噪声,步骤3)中在声波测井数据代替深度域层速度模型中有井区域的地震速度之前,还包括:采用滑动平均数字滤波法去除声波测井数据的噪声,得到去噪处理后的声波测井数据。
为了实现声波测井数据对深度域层速度模型中有井区域的地震速度的代替,得到深度域层速度模型之后,还包括:根据地震解释层位数据将深度域层速度模型在三维空间划分成不同深度区域的深度域层速度模型,根据划分的各个深度区域,将声波测井数据代替相应深度区域的深度域层速度模型中有井区域的地震速度。
附图说明
图1是本发明的深度域速度建模方法流程图;
图2是本发明的初始的深度域速度模型示意图;
图3是本发明的处理后的测井声波速度的示意图;
图4是本发明在测井约束下的初始深度域速度模型示意图;
图5是本发明的优化深度域速度模型示意图;
图6-1是本发明的设定区域的叠前时间偏移成像示意图;
图6-2是本发明的设定区域的速度模型叠前深度偏移成像示意图;
图7-1是本发明的另一设定区域的叠前时间偏移成像示意图;
图7-2是本发明的另一设定区域的速度模型叠前深度偏移成像示意图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的说明。
本发明以四川盆地某探区三维地震数据为例,提出一种深度域速度建模方法,如图1所示,包括以下步骤:
获取该四川盆地某探区范围内的地震数据、声波测井数据和地震解释层位数据。其中,地震数据包括三维叠前时间偏移速度,声波测井数据包括深度信息、与深度对应的声波速度信息,地震解释层位数据为设定区域的深度解释层位信息。
获取的三维地震数据应为经过了叠前静校正处理、叠前去噪处理、地表一致性振幅处理、叠前反褶积处理和剩余静校正处理后的共中心点道集数据,排列方式应按照主测线、联络线、偏移距增序排列。
获取深度解释层位信息的设定区域按照深度从浅至深依次包括:千佛崖组底部层位、须家河组须五段底部层位、须家河组须三段底部层位、嘉陵江组嘉四五段底部层位、飞仙关组飞一二段底部层位。
本发明获取的三维叠前时间偏移速度是共中心点道集数据中在主测线的方向间隔250m的稀疏离散速度点组成的离散的速度线,和在联络线的方向间隔500m的稀疏离散速度点组成的离散的速度线。利用已知的所有离散的速度线在四川盆地某探区的范围采用克里金插值方法计算剩余的空点的速度值。克里金插值方法为:令i,k表示速度模型(x,z)方向的离散点,i=1,2…,n;zi,k为对应离散点的深度,νi,k为对应离散点的速度,则i′,k处的速度值ν′i,k可采用一个线性组合得到:
其中,λi为权值,λi可选取估计方差最小的方法计算得到,估计方差最小的方法的计算式如下:
其中,C(xi-xi′)为νi,k和νi′,k的协方差函数。
然后,利用三维叠前时间偏移速度进行建模,得到初始时间域的速度模型,对初始时间域的速度模型进行一次DIX转换,得到时间域层速度模型,对时间域层速度模型进行时深转换,得到如图2所示的深度域层速度模型,从图2中可以看出深度域层速度模型与地下构造吻合度不高,纵向分辨率低。
对初始时间域的速度模型通过DIX转换,得到时间域层速度模型,对应的转换公式为:
其中,νint为时间域层速度,t1为第一反射层的旅行时,t2为第二反射层的旅行时,νrms1为第一反射层的速度,νrms2为第二反射层的速度。
根据目标层纵向分辨率的要求(识别小于10m的断层),对声波测井数据进行重采样,重采样间隔为5m,且采样点之间深度间隔一致,得到重采样处理后的声波测井数据。以p9井的声波测井数据为例,如图3所示,深度-速度对中深度采样不均匀,速度纵向分辨率高,包括高频成分(>25Hz)和中高频成分(15-25Hz),且高频成分中含有异常值和高频振荡噪声。
对p9井深度采样不均匀的问题,以纵向5m的深度间隔重新采样。假设p9井的井速度深度采样点1,2,3,…,i,…对应的深度为d(1),d(2),d(3),...,d(i)...,深度采样间隔为Δd(i),对应的测井速度为ν(1),ν(2),ν(3),...,ν(i)...,其中,由于深度采样间隔为不规则采样,导致每两个深度点之间对应的采样间隔Δd(i)不同。以5m规则间隔对原深度点重新采样,原深度值对应的新采样点为为四舍五入取整,得到深度间隔一致的声波测井数据。
以p9井为例,原深度点d(1),d(2),d(3),d(4),d(5),d(6),...,深度数值分别为164.875,165,165.125,165.25,165.625,165.75,……。则重采样后深度点深度值分别为: 即32,32,33,33,33,33,…,当相邻的两个或多个点重采样后具有相同数值时,选择原数值最接近的点为重采样后的数值。165与165的差值小于162.125,165.25,165.625与165的差值,因此选择165重采样得到33及对应的速度值。
然后,采用滑动平均数字滤波法去除重采样处理后的声波测井数据的噪声,具体的,利用滑动平均数字滤波法去除p9井重采样处理后的声波测井数据的异常值干扰(即噪声):
在当前深度采样点前、当前深度采样点后分别连续地取m个深度采样点,用(2m+1)个深度采样点的测井速度值(包括当前深度采样点的测井速度值在内),依次地计算出当前深度采样点的测井速度值的滑动平均值;再取下一深度采样点作为当前深度采样点,按照上述方法计算下一深度采样点的测井速度值的滑动平均值;然后,再取下一深度采样点,直到计算出所有深度采样点的测井速度值的滑动平均值,就此,便可消除异常值干扰。例如,以m=2为例,当前深度采样点的测井速度值的滑动平均值的计算公式如下:
其中,νi为当前深度采样点的测井速度值,为当前深度采样点的测井速度值的滑动平均值(即去噪处理后的声波测井数据),νi-1为当前深度采样点前一个深度采样点的测井速度值,νi-2为νi-1前一个深度采样点的测井速度值,νi+1当前深度采样点后一个深度采样点的测井速度值,νi+2为νi+1后一个深度采样点的测井速度值。
利用地震解释层位数据中的深度解释层位信息将深度域层速度模型在三维空间切分成不同深度区域的深度域层速度模型,根据划分的各个深度区域的深度域层速度模型中有井区域的地震速度,替换为重采样处理和去噪处理后的声波测井数据,作为新的有井区域的速度模型。
对于不同深度区域的深度域层速度模型中的无井区域,按照滑动加权计算,采用克里金差值方法沿着地震层位对邻近区域的井速度和地震速度进行不同权重的插值计算得到速度值,得到无井区域的速度模型;将有井区域的速度模型与无井区域的速度模型合并,得到初始的深度域速度模型,如图4所示。
利用地震解释层位数据,按照解释层位从浅至深,对初始深度域速度模型进行网格层析反演优化,得到优化的深度域速度模型,如图5所示。图6-1为设定区域的工区叠前时间偏移剖面,图6-2为应用本发明的方法在设定区域的叠前深度偏移剖面,叠前深度偏移结果中偏移噪声更小,地层清晰可连续追踪,断层成像更好。图7-1为另一设定区域的工区叠前时间偏移剖面,图7-2为应用本发明的方法在另一设定区域的叠前深度偏移剖面,叠前深度偏移结果中负向区断裂更清楚,断层关系更明确。
另外,作为其他实施方式,还可以不对声波测井数据进行重采样处理和去噪处理,直接将声波测井数据代替工区范围内有井区域的地震速度,作为新的有井区域的速度模型;或者,只将声波测井数据进行重采样处理或去噪处理,将重采样处理或去噪处理后的声波测井数据代替工区范围内有井区域的地震速度,作为新的有井区域的速度模型。
以上所述仅为本发明的优选实施例,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (8)
1.一种深度域速度建模方法,其特征在于,包括以下步骤:
1)获取工区范围内的地震数据、声波测井数据和地震解释层位数据,其中,地震数据包括三维叠前时间偏移速度,声波测井数据包括深度、与深度对应的声波速度,地震解释层位数据为设定区域的深度解释层位;
2)利用三维叠前时间偏移速度进行建模,得到初始时间域的速度模型,对初始时间域的速度模型进行两次DIX转换,得到深度域层速度模型;
3)将声波测井数据代替深度域层速度模型中有井区域的地震速度,得到有井区域的速度模型;对于深度域层速度模型中的无井区域,按照滑动加权计算方法,采用克里金差值方法沿着地震层位对邻近区域的井速度和地震速度进行不同权重的插值计算得到速度值,得到无井区域的速度模型;将有井区域的速度模型与无井区域的速度模型合并,得到初始深度域速度模型;
4)利用地震解释层位数据,对初始深度域速度模型进行网格层析反演优化,得到优化的深度域速度模型。
2.根据权利要求1所述的深度域速度建模方法,其特征在于,步骤3)中在声波测井数据代替深度域层速度模型中有井区域的地震速度之前,还包括:根据目标层纵向分辨率的要求,对声波测井数据进行重采样,且采样点之间深度间隔一致,得到重采样处理后的声波测井数据。
3.根据权利要求1所述的深度域速度建模方法,其特征在于,步骤3)中在声波测井数据代替深度域层速度模型中有井区域的地震速度之前,还包括:采用滑动平均数字滤波法去除声波测井数据的噪声,得到去噪处理后的声波测井数据。
4.根据权利要求1所述的深度域速度建模方法,其特征在于,得到深度域层速度模型之后,还包括:根据所述地震解释层位数据将深度域层速度模型在三维空间划分成不同深度区域的深度域层速度模型,根据划分的各个深度区域,将声波测井数据代替相应深度区域的深度域层速度模型中有井区域的地震速度。
5.一种深度域速度建模装置,其特征在于,包括存储器和处理器,以及存储在所述存储器上并在所述处理器上运行的计算机程序,所述处理器与所述存储器相耦合,所述处理器执行所述计算机程序时实现以下步骤:
1)获取工区范围内的地震数据、声波测井数据和地震解释层位数据,其中,地震数据包括三维叠前时间偏移速度,声波测井数据包括深度、与深度对应的声波速度,地震解释层位数据为设定区域的深度解释层位;
2)利用三维叠前时间偏移速度进行建模,得到初始时间域的速度模型,对初始时间域的速度模型进行两次DIX转换,得到深度域层速度模型;
3)将声波测井数据代替深度域层速度模型中有井区域的地震速度,得到有井区域的速度模型;对于深度域层速度模型中的无井区域,按照滑动加权计算方法,采用克里金差值方法沿着地震层位对邻近区域的井速度和地震速度进行不同权重的插值计算得到速度值,得到无井区域的速度模型;将有井区域的速度模型与无井区域的速度模型合并,得到初始深度域速度模型;
4)利用地震解释层位数据,对初始深度域速度模型进行网格层析反演优化,得到优化的深度域速度模型。
6.根据权利要求5所述的深度域速度建模装置,其特征在于,步骤3)中在声波测井数据代替深度域层速度模型中有井区域的地震速度之前,还包括:根据目标层纵向分辨率的要求,对声波测井数据进行重采样,且采样点之间深度间隔一致,得到重采样处理后的声波测井数据。
7.根据权利要求5所述的深度域速度建模装置,其特征在于,步骤3)中在声波测井数据代替深度域层速度模型中有井区域的地震速度之前,还包括:采用滑动平均数字滤波法去除声波测井数据的噪声,得到去噪处理后的声波测井数据。
8.根据权利要求5所述的深度域速度建模装置,其特征在于,得到深度域层速度模型之后,还包括:根据所述地震解释层位数据将深度域层速度模型在三维空间划分成不同深度区域的深度域层速度模型,根据划分的各个深度区域,将声波测井数据代替相应深度区域的深度域层速度模型中有井区域的地震速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811217379.6A CN111077575B (zh) | 2018-10-18 | 2018-10-18 | 一种深度域速度建模方法与装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811217379.6A CN111077575B (zh) | 2018-10-18 | 2018-10-18 | 一种深度域速度建模方法与装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111077575A true CN111077575A (zh) | 2020-04-28 |
CN111077575B CN111077575B (zh) | 2021-11-19 |
Family
ID=70308904
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811217379.6A Active CN111077575B (zh) | 2018-10-18 | 2018-10-18 | 一种深度域速度建模方法与装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111077575B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111538080A (zh) * | 2020-05-28 | 2020-08-14 | 中国石油天然气集团有限公司 | 地震成像的方法 |
CN111624649A (zh) * | 2020-06-05 | 2020-09-04 | 中油奥博(成都)科技有限公司 | 一种零偏移距vsp建立横向变速层速度模型的方法和装置 |
CN111736213A (zh) * | 2020-07-07 | 2020-10-02 | 中油奥博(成都)科技有限公司 | 一种变偏移距VSP Kirchhoff偏移速度分析方法和装置 |
CN112529981A (zh) * | 2020-11-24 | 2021-03-19 | 中国煤炭地质总局地球物理勘探研究院 | 一种深度域煤层底板等高线成图的方法 |
CN113945978A (zh) * | 2020-07-15 | 2022-01-18 | 中国石油天然气股份有限公司 | 一种深度域地震数据优化方法及装置 |
CN113960658A (zh) * | 2020-07-20 | 2022-01-21 | 中国石油天然气股份有限公司 | 一种基于地质地震模型的测井约束速度建模方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090010102A1 (en) * | 2007-06-29 | 2009-01-08 | Karine Broto | Method for adjusting a seismic wave velocity model according to information recorded in wells |
CN106443774A (zh) * | 2016-11-16 | 2017-02-22 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种提高起伏地表叠前深度偏移成像精度的方法 |
CN107942379A (zh) * | 2017-10-12 | 2018-04-20 | 中国石油化工股份有限公司 | 一种提高复杂断块速度模型精度的方法 |
-
2018
- 2018-10-18 CN CN201811217379.6A patent/CN111077575B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090010102A1 (en) * | 2007-06-29 | 2009-01-08 | Karine Broto | Method for adjusting a seismic wave velocity model according to information recorded in wells |
CN106443774A (zh) * | 2016-11-16 | 2017-02-22 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种提高起伏地表叠前深度偏移成像精度的方法 |
CN107942379A (zh) * | 2017-10-12 | 2018-04-20 | 中国石油化工股份有限公司 | 一种提高复杂断块速度模型精度的方法 |
Non-Patent Citations (3)
Title |
---|
WOOK B. LEE 等: "三维地质统计速度模型的建立:地压环境下的盐成像", 《国外油气勘探》 * |
YANG CHONGYANG ET AL.: "Application and Discussion of Well logging control Velocity Depth Domain Modeling Technology in the East of Tarim basin Desert Arear", 《SEG FULL WAVEFORM INVERSION: DEPTH MODEL BUILDING WORKSHOP》 * |
罗勇等: "复杂构造地震叠前深度偏移速度模型构建及效果", 《新疆石油地质》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111538080A (zh) * | 2020-05-28 | 2020-08-14 | 中国石油天然气集团有限公司 | 地震成像的方法 |
CN111624649A (zh) * | 2020-06-05 | 2020-09-04 | 中油奥博(成都)科技有限公司 | 一种零偏移距vsp建立横向变速层速度模型的方法和装置 |
CN111624649B (zh) * | 2020-06-05 | 2022-05-20 | 中油奥博(成都)科技有限公司 | 一种零偏移距vsp建立横向变速层速度模型的方法和装置 |
CN111736213A (zh) * | 2020-07-07 | 2020-10-02 | 中油奥博(成都)科技有限公司 | 一种变偏移距VSP Kirchhoff偏移速度分析方法和装置 |
CN111736213B (zh) * | 2020-07-07 | 2022-05-20 | 中油奥博(成都)科技有限公司 | 一种变偏移距VSP Kirchhoff偏移速度分析方法和装置 |
CN113945978A (zh) * | 2020-07-15 | 2022-01-18 | 中国石油天然气股份有限公司 | 一种深度域地震数据优化方法及装置 |
CN113945978B (zh) * | 2020-07-15 | 2023-09-26 | 中国石油天然气股份有限公司 | 一种深度域地震数据优化方法及装置 |
CN113960658A (zh) * | 2020-07-20 | 2022-01-21 | 中国石油天然气股份有限公司 | 一种基于地质地震模型的测井约束速度建模方法及装置 |
CN113960658B (zh) * | 2020-07-20 | 2024-01-30 | 中国石油天然气股份有限公司 | 一种基于地质地震模型的测井约束速度建模方法及装置 |
CN112529981A (zh) * | 2020-11-24 | 2021-03-19 | 中国煤炭地质总局地球物理勘探研究院 | 一种深度域煤层底板等高线成图的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111077575B (zh) | 2021-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111077575B (zh) | 一种深度域速度建模方法与装置 | |
CN108802812B (zh) | 一种井震融合的地层岩性反演方法 | |
EP2062071B1 (en) | Method for building velocity models for pre-stack depth migration via the simultaneous joint inversion of seismic, gravity and magnetotelluric data | |
CN107783187B (zh) | 一种将测井速度和地震速度结合建立三维速度场的方法 | |
CN113552625B (zh) | 一种用于常规陆域地震数据的多尺度全波形反演方法 | |
CN106646613A (zh) | 深度域多尺度井控建模与成像联合处理方法 | |
CN113740901A (zh) | 基于复杂起伏地表的陆上地震数据全波形反演方法及装置 | |
CN102338887B (zh) | 不规则尺寸空变网格层析成像静校正方法 | |
WO2005026776A1 (en) | Wide-offset-range pre-stack depth migration method for seismic exploration | |
CN111123359B (zh) | 随钻测井与地层格架约束的井周地震成像探测方法及装置 | |
CN106842291B (zh) | 一种基于叠前地震射线阻抗反演的不整合圈闭储层岩性预测方法 | |
CN103217715B (zh) | 多尺度规则网格层析反演静校正方法 | |
AU739128B2 (en) | A method of seismic processing, and in particular a 3D seismic prospection method implementing seismic data migration | |
CN111538080B (zh) | 地震成像的方法 | |
CN112147700A (zh) | 速度异常区的低频模型构建方法及系统 | |
CN111060967A (zh) | 一种井控速度场建模方法 | |
CN110967757A (zh) | 随钻实时叠前地震反演方法及系统 | |
CN113960658B (zh) | 一种基于地质地震模型的测井约束速度建模方法及装置 | |
CN113433588B (zh) | 一种分偏移距扫描叠加的近地表速度分析方法 | |
CN111175825B (zh) | 深度域速度建模方法 | |
CN113945978B (zh) | 一种深度域地震数据优化方法及装置 | |
CN112014875B (zh) | 叠前地震反演方法及装置 | |
CN112444848B (zh) | 一种全波形反演的方法及系统 | |
CN114814949B (zh) | 一种浅层逆vsp初至层析及地层预测方法 | |
CN114966826B (zh) | 基于野外近地表调查资料的高低频融合近地表建模方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |