CN106154329B - 利用vsp速度校正声波均方根速度的方法及其系统 - Google Patents
利用vsp速度校正声波均方根速度的方法及其系统 Download PDFInfo
- Publication number
- CN106154329B CN106154329B CN201510181167.7A CN201510181167A CN106154329B CN 106154329 B CN106154329 B CN 106154329B CN 201510181167 A CN201510181167 A CN 201510181167A CN 106154329 B CN106154329 B CN 106154329B
- Authority
- CN
- China
- Prior art keywords
- sound wave
- velocity
- root mean
- mean sequare
- vsp
- 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.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种利用VSP速度校正声波均方根速度的方法及其系统,该方法包括步骤:S1、对声波时差曲线进行井眼环境影响校正后,将声波时差转换成声波层速度;S2、对声波层速度进行深时转换,得到时间域的声波层速度,计算得到声波均方根速度并计算VSP均方根速度S3、计算每个时间点的声波均方根速度与VSP均方根速度的误差Δvrms(ti):S4、设计一个多项式:f(t)=a0+a1t+a2t2+…+antn,使其与Δvrms(ti)曲线最佳匹配,从而确定系数a0,a1,a2…an;S5、计算校正后的声波均方根速度本发明减少了拾取的随机性,从而提高了均方根速度的精度。
Description
技术领域
本发明涉及声波测井处理的领域,尤其涉及一种利用VSP速度校正声波均方根速度的方法及其系统。
背景技术
众所周知,由于声波测井和地震勘探的激发震源不同,接受方式也不同,导致对于同一地层,由声波测井得到的速度(均方根速度和层速度)与地震资料速度分析得到的速度数值不相等,因此,声波测井速度不能直接应用到地震资料速度分析中,必须经过校正。
均方根速度是地震速度谱拾取的速度,由它通过DIX公式转换得到地层的层速度,建立层速度模型,应用到地震资料的叠前深度偏移处理中,因此,均方根速度正确与否直接影响层速度模型的正确性。由于VSP(Vetical Seismic Profile,垂直地震剖面)测井得到的速度不仅准确,而且与地面地震得到的速度数值相当,遗憾的是VSP测井数量很少,资料有限,而声波测井数量多,资料丰富,因此,目前多数研究者采用VSP层速度校正声波测井层速度,但实现方法比较繁琐,需要分层段进行。
发明内容
本发明的特征和优点在下文的描述中部分地陈述,或者可从该描述显而易见,或者可通过实践本发明而学习。
为克服现有技术的问题,本发明提供一种利用VSP速度校正声波均方根速度的方法及其系统,采用均方根速度误差曲线的多项式拟合,从而减少拾取的随机性,从而提高均方根速度的精度。
本发明解决上述技术问题所采用的技术方案如下:
根据本发明的一个方面,提供一种利用VSP速度校正声波均方根速度的方法,其特征在于,包括:
S1、对声波时差曲线进行井眼环境影响校正后,将声波时差转换成声波层速度;
S2、对声波层速度进行深时转换,得到时间域的声波层速度,计算得到声波均方根速度并计算VSP均方根速度
S3、计算每个时间点的声波均方根速度与VSP均方根速度的误差:
S4、设计一个多项式:
使其与Δvrms(ti)曲线最佳匹配,从而确定系数a0,a1,a2…
S5、计算校正后的声波均方根速度
根据本发明的一个实施例,在步骤S1中,声波层速度为声波时差的倒数。
根据本发明的一个实施例,在步骤S2中,在对声波层速度进行深时转换时,采取的公式为:
其中,d是深度采样间隔,vj是声波层速度,ti是第i个点的时间。
根据本发明的一个实施例,在步骤S2中,在计算声波均方根速度时,采取的公式为:
其中,v(i)和分别是第i个点的时间域的声波层速度和时间值。
根据本发明的一个实施例,在步骤S2中,在计算VSP均方根速度时,采取的公式为:
其中,v(i)和分别是第i个点的VSP层速度和时间值。
根据本发明的另一个方面,还提供一种利用VSP速度校正声波均方根速度的系统,其特征在于,包括:
声波层速度计算单元,用于对声波时差曲线进行井眼环境影响校正后,将声波时差转换成声波层速度;
均方根速度计算单元,用于对声波层速度进行深时转换,得到时间域的声波层速度,计算得到声波均方根速度还用于计算VSP均方根速度
误差计算单元,用于计算每个时间点的声波均方根速度与VSP均方根速度的误差:
多项式拟合系数确定单元,用于设计一个多项式:f(t)=a0+a1t+a2t2+...+antn
使其与Δvrms(ti)曲线最佳匹配,从而确定系数a0,a1,a2…an
声波均方根速度校正计算单元,用于计算校正后的声波均方根速度
根据本发明的一个实施例,声波层速度为声波时差的倒数。
根据本发明的一个实施例,均方根速度计算单元,用于采取以下公式对声波层速度进行深时转换:
其中,d是深度采样间隔,vj是声波层速度,ti是第i个点的时间。
根据本发明的一个实施例,均方根速度计算单元,还用于采取以下公式计算声波均方根速度:
其中,v(i)和分别是第i个点的时间域的声波层速度和时间值。
根据本发明的一个实施例,均方根速度计算单元,还用于采取以下公式计算与VSP均方根速度:
其中,v(i)和分别是第i个点的VSP层速度和时间值。
通过阅读说明书,本领域普通技术人员将更好地了解这些技术方案的特征和内容。
附图说明
下面通过参考附图并结合实例具体地描述本发明,本发明的优点和实现方式将会更加明显,其中附图所示内容仅用于对本发明的解释说明,而不构成对本发明的任何意义上的限制,在附图中:
图1为本发明实施例的利用VSP速度校正声波均方根速度的方法的流程示意图。
图2为本发明实施例的利用VSP速度校正声波均方根速度的系统的结构示意图。
图3为本发明实施例的某井未经校正的声波均方根速度和VSP均方根速度的曲线示意图。
图4为本发明实施例的某井经校正后的声波均方根速度和VSP均方根速度的曲线示意图。
具体实施方式
如图1所示,本发明提供一种利用VSP速度校正声波均方根速度的方法,其特征在于,包括:
S1、对声波时差曲线进行井眼环境影响校正后,将声波时差转换成声波层速度;
S2、对声波层速度进行深时转换,得到时间域的声波层速度,计算得到声波均方根速度并计算VSP均方根速度
S3、计算每个时间点的声波均方根速度与VSP均方根速度的误差:
S4、设计一个多项式:f(t)=a0+a1t+a2t2+...+antn
使其与Δvrms(ti)曲线最佳匹配,从而确定系数a0,a1,a2…an
S5、计算校正后的声波均方根速度
在步骤S1中,对声波时差曲线进行井眼环境影响校正后,还可以进行编辑消除奇异值,再将声波时差转换成声波层速度,其中,声波层速度为声波时差的倒数。
在步骤S2中,在对声波层速度进行深时转换,采取的公式为:
其中,j=1,2,3,...i,i=1,2,3,...N,N为总点数,d是深度采样间隔,vj是声波层速度,第i个点的时间ti是前i个点的层间时间的累加求和。
通过上述公式计算深度域的层速度对应的时间值,就得到了时间域的层速度,在本实施例中,时间不是等间隔采样,但这不会影响后续均方根速度的计算。
而在计算声波均方根速度时,采取的公式为:
其中,v(i)和分别是第i个点的时间域的声波层速度和时间值。更具体来说,v(i)是时间域的声波测井速度,就是由深度域的层速度经过深时转换得到的,其数值大小没有改变,即vj与v(i)在数值上是相等的,只是变到了时间域;而就是ti,即第i个点的时间。
在计算VSP均方根速度时,采取的公式同样为:
此时,v(i)和分别是第i个点的VSP层速度和时间值。
本发明还提供一种利用VSP速度校正声波均方根速度的系统,如图2所示,包括:
声波层速度计算单元11,用于对声波时差曲线进行井眼环境影响校正后,将声波时差转换成声波层速度;
均方根速度计算单元12,用于对声波层速度进行深时转换,得到时间域的声波层速度,计算得到声波均方根速度还用于计算VSP均方根速度
误差计算单元13,用于计算每个时间点的声波均方根速度与VSP均方根速度的误差:
多项式拟合系数确定单元14,用于设计一个多项式:f(t)=a0+a1t+a2t2+...+antn
使其与Δvrms(ti)曲线最佳匹配,从而确定系数a0,a1,a2…an
声波均方根速度校正计算单元15,用于计算校正后的声波均方根速度
声波层速度计算单元11还用于对声波时差曲线进行井眼环境影响校正后,进行编辑消除奇异值。声波层速度计算单元11计算得出的声波层速度为声波时差的倒数。
均方根速度计算单元12,用于采取以下公式对声波层速度进行深时转换:
其中,j=1,2,3,...i,i=1,2,3,...N,N为总点数,d是深度采样间隔,vj是声波层速度,第i个点的时间ti是前i个点的层间时间的累加求和。
均方根速度计算单元12还用于采取以下公式计算声波均方根速度:
其中,v(i)和分别是第i个点的时间域的声波层速度和时间值。
均方根速度计算单元12还用于采取同样的公式计算VSP均方根速度:
此时,v(i)和分别是第i个点的VSP层速度和时间值。
请参照图3,图3为本发明实施例的某井未经校正的声波均方根速度和VSP均方根速度的曲线示意图,其中黑线表示的是未经校正的声波均方根速度,虚线表示的是VSP均方根速度,可以看到,声波均方根速度和VSP均方根速度是不一致的,尤其是浅层,差别最大,在700ms附近,达到500m/s。请再参照图4,图4为本发明实施例的某井经校正后的声波均方根速度和VSP均方根速度的曲线示意图,其中黑线表示的是经校正后的声波均方根速度,虚线表示的是VSP均方根速度,可以看到,经过校正后的声波均方根速度和VSP均方根速度基本一致。由此可见,本发明提供的利用VSP速度校正声波均方根速度的方法及其系统是有效的。
本发明提出了一种提供的利用VSP速度校正声波均方根速度的方法及其系统,该方法包括有四个处理步骤:①声波时差曲线的预处理;②深时转换和均方根速度计算;③均方根速度误差曲线的多项式拟合;④根据拟合公式对声波均方根速度校正。其中关键是第③步。该方法较传统的层速度校正方法简单、快速,很容易实现,通过实际验证了该方法的有效性。
以上参照附图说明了本发明的优选实施例,本领域技术人员不脱离本发明的范围和实质,可以有多种变型方案实现本发明。举例而言,作为一个实施例的部分示出或描述的特征可用于另一实施例以得到又一实施例。以上仅为本发明较佳可行的实施例而已,并非因此局限本发明的权利范围,凡运用本发明说明书及附图内容所作的等效变化,均包含于本发明的权利范围之内。
Claims (10)
1.一种利用VSP速度校正声波均方根速度的方法,其特征在于,包括:
S1、对声波时差曲线进行井眼环境影响校正后,将声波时差转换成声波层速度;
S2、对所述声波层速度进行深时转换,得到时间域的声波层速度,计算得到声波均方根速度并计算VSP均方根速度
S3、计算每个时间点的声波均方根速度与VSP均方根速度的误差:
S4、设计一个多项式:f(t)=a0+a1t+a2t2+…+antn
使其与Δvrms(ti)曲线最佳匹配,从而确定系数a0,a1,a2…an
S5、计算校正后的声波均方根速度
2.根据权利要求1所述利用VSP速度校正声波均方根速度的方法,其特征在于,在步骤S1中,所述声波层速度为所述声波时差的倒数。
3.根据权利要求1所述利用VSP速度校正声波均方根速度的方法,其特征在于,在所述步骤S2中,在对所述声波层速度进行深时转换时,采取的公式为:
其中,d是深度采样间隔,vj是声波层速度,ti是第i个点的时间。
4.根据权利要求1所述利用VSP速度校正声波均方根速度的方法,其特征在于,在所述步骤S2中,在计算所述声波均方根速度时,采取的公式为:
其中,v(i)和分别是第i个点的时间域的声波层速度和时间值。
5.根据权利要求1所述利用VSP速度校正声波均方根速度的方法,其特征在于,在所述步骤S2中,在计算所述VSP均方根速度时,采取的公式为:
其中,v'(i)和分别是第i个点的VSP层速度和时间值。
6.一种利用VSP速度校正声波均方根速度的系统,其特征在于,包括:
声波层速度计算单元,用于对声波时差曲线进行井眼环境影响校正后,将声波时差转换成声波层速度;
均方根速度计算单元,用于对所述声波层速度进行深时转换,得到时间域的声波层速度,计算得到声波均方根速度还用于计算VSP均方根速度
误差计算单元,用于计算每个时间点的声波均方根速度与VSP均方根速度的误差:
多项式拟合系数确定单元,用于设计一个多项式:f(t)=a0+a1t+a2t2+…+antn
使其与Δvrms(ti)曲线最佳匹配,从而确定系数a0,a1,a2…an
声波均方根速度校正计算单元,用于计算校正后的声波均方根速度
7.根据权利要求6所述利用VSP速度校正声波均方根速度的系统,其特征在于,所述声波层速度为所述声波时差的倒数。
8.根据权利要求6所述利用VSP速度校正声波均方根速度的系统,其特征在于,所述均方根速度计算单元,用于采取以下公式对所述声波层速度进行深时转换:
其中,d是深度采样间隔,vj是声波层速度,ti是第i个点的时间。
9.根据权利要求6所述利用VSP速度校正声波均方根速度的系统,其特征在于,所述均方根速度计算单元,还用于采取以下公式计算所述声波均方根速度:
其中,v(i)和分别是第i个点的时间域的声波层速度和时间值。
10.根据权利要求6所述利用VSP速度校正声波均方根速度的系统,其特征在于,所述均方根速度计算单元,还用于采取以下公式计算所述与VSP均方根速度:
其中,v'(i)和分别是第i个点的VSP层速度和时间值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510181167.7A CN106154329B (zh) | 2015-04-16 | 2015-04-16 | 利用vsp速度校正声波均方根速度的方法及其系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510181167.7A CN106154329B (zh) | 2015-04-16 | 2015-04-16 | 利用vsp速度校正声波均方根速度的方法及其系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106154329A CN106154329A (zh) | 2016-11-23 |
CN106154329B true CN106154329B (zh) | 2018-08-31 |
Family
ID=58058057
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510181167.7A Active CN106154329B (zh) | 2015-04-16 | 2015-04-16 | 利用vsp速度校正声波均方根速度的方法及其系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106154329B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111060969B (zh) * | 2019-12-25 | 2022-03-18 | 恒泰艾普(北京)能源科技研究院有限公司 | 一种井控q补偿方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7525874B2 (en) * | 2006-04-21 | 2009-04-28 | Prism Seismic Inc. | Method for converting seismic data from the time domain to the depth domain |
CN101533103B (zh) * | 2009-04-13 | 2012-03-21 | 中国石油天然气集团公司 | 地震波速度场构建中井约束横向可变h-v曲线构建方法 |
CN103790569B (zh) * | 2012-10-26 | 2016-12-21 | 中国石油化工股份有限公司 | 一种利用vsp时深关系对声波测井资料进行校正的方法 |
CN103837893B (zh) * | 2014-03-10 | 2017-06-30 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于vsp资料的声波测井曲线校正方法 |
-
2015
- 2015-04-16 CN CN201510181167.7A patent/CN106154329B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106154329A (zh) | 2016-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Harri et al. | Pressure observations by the Curiosity rover: Initial results | |
Adams et al. | Temporal variability of near‐bottom dissolved oxygen during upwelling off central Oregon | |
CN111273318B (zh) | 一种基于抛物线的区域对流层湿延迟计算方法 | |
CN109581501B (zh) | 用于沙漠区深度域速度建模的方法 | |
CN103837895A (zh) | 拟合初至波获取短波长静校正量的方法 | |
CN104570102B (zh) | 一种近地表速度模型与中深层速度模型的融合方法 | |
CN104118578B (zh) | 一种微小卫星平台多个敏感器数据动态融合系统及方法 | |
CN105588883B (zh) | 三维岩石力学参数获取方法和系统 | |
CN107817516B (zh) | 基于初至波信息的近地表建模方法及系统 | |
CN113552625B (zh) | 一种用于常规陆域地震数据的多尺度全波形反演方法 | |
CN104459782A (zh) | 使用薄层的沿层层析反演地震速度建模方法以及建模单元 | |
CN104330827B (zh) | 一种地表模型静校正的处理方法 | |
CN104678440A (zh) | 井约束二维地震变速度场非线性误差校正方法 | |
CN104865597A (zh) | 一种深度域层速度初始模型的建模方法 | |
CN103412333A (zh) | 一种静校正基准面确定方法 | |
CN104765055A (zh) | Gps测站坐标时间序列周期性探测方法及系统 | |
CN103576197B (zh) | 一种转换波角道集抽取方法 | |
NO20101734A1 (no) | Hastighetsmodeller for en enkeltbronner og for et sett av bronner | |
CN109931952A (zh) | 未知纬度条件下捷联惯导直接解析式粗对准方法 | |
CN106154329B (zh) | 利用vsp速度校正声波均方根速度的方法及其系统 | |
CN104914469B (zh) | 一种转换横波中的静校正方法及装置 | |
CN103487834A (zh) | 转换波共检波点叠加静校正方法 | |
CN103790569B (zh) | 一种利用vsp时深关系对声波测井资料进行校正的方法 | |
CN101950034B (zh) | 一种用于地震反演的自然电位曲线基线漂移校正方法 | |
CN102353992B (zh) | 基于中心点规则化的起伏地表动校正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |