System.ArgumentOutOfRangeException: 索引和长度必须引用该字符串内的位置。 参数名: length 在 System.String.Substring(Int32 startIndex, Int32 length) 在 zhuanliShow.Bind()
【技术实现步骤摘要】
本专利技术属于石油天然气勘探领域,特别涉及一种时差搜索范围动态调整的地层纵波时差提取方法、系统。
技术介绍
1、利用阵列声波测井数据提取地层纵波时差是测井资料处理解释的重要环节。目前,地层纵波时差提取主要采用时间-时差相关分析法(stc),该方法基于阵列声波测井中不同接收器接收的同一模式波波形具有相似性这一特点来计算地层纵波时差。传统stc方法计算时,首先在时间域和时差域内,采用一定长度的窗长求取不同时间和不同时差下不同接收器上波形的相关系数;其次,根据实际地层情况,设置纵波时差范围,在该时差范围内搜索时间-时差对应的二维相关系数谱的极大值;最后,利用相关系数极大值对应的时差来确定地层纵波时差。申请号为cn201611236689.3的中国专利公开了一种lwf存储式声波测井慢度提取方法,该方法通过先获取测井声波波形数据;再确定声波测井深度域上声波接收器r1接收声波的首波负峰序列ti;在深度域上分别对调整后的声波接收器r1接收的声波数据、声波接收器r2接收的声波数据和声波接收器r3接收的声波数据进行基线归零以及样条拟合;然后确定调整后的声波接收器r1的纵波窗长以及慢度s搜索范围;五、慢度序列si的提取:501、确定慢度分割步长δs,502、获取共源法相关匹配中的慢度序列sfi,503、获取共收法相关匹配中的慢度sri,504、计算深度域上声波测井的慢度序列si。该方法通过涉及到多个接收器获取声波,涉及的计算量较大。
2、但随着勘探目标的日益复杂,同一井中不同井段地层岩性多变,纵波时差往往差别较大,整个井段采用同样的搜索
技术实现思路
1、针对上述问题,一方面本专利技术提出了一种时差搜索范围动态调整的地层纵波时差提取方法,所述提取方法包括:
2、分析1号接收器波形的能量比值随波形记录时间的变化,确定最大能量比值对应的时间作为1号接收器的纵波到时;
3、计算根据所述1号接收器的纵波到时,计算时差域内不同时差对应的波形相关系数;
4、动态调整纵波时差的搜索范围,在搜索范围内拾取所述波形相关系数的最大值对应的时差作为地层的纵波时差;
5、重复上述步骤,直到完成整个测量井段的纵波时差提取。
6、进一步地,分析1号接收器波形的能量比值随波形记录时间的变化采用长短时窗能量比值法。
7、进一步地,在分析1号接收器波形的能量比值随波形记录时间的变化之前还包括:
8、在目的层段利用多极子阵列声波测井仪器进行井下测量,收集全波测量数据,并对所述全波测量数据进行预处理。
9、进一步地,对所述全波测量数据进行预处理包括:对全波测量数据进行深度偏移校正、波形开始时间校正、波形去增益和数字滤波。
10、进一步地,确定最大能量比值对应的时间采用以下公式:
11、
12、式中,r(i)为第i个时间采样点处的能量比值;wk为第k个时间采样点的波形值;n1为前时窗包含的时间采样点数,n2为后时窗包含的时间采样点数;最大能量比值对应的波形记录时间即为1号接收器的纵波到时。
13、进一步地,计算时差域内不同时差对应的波形相关系数采用以下公式:
14、
15、式中,ρ(s)为不同时差对应的波形相关系数;s为地层模式波时差;tw为相关分析窗长;m为仪器接收器个数;δz为接收器间距;m为接收器序号;t1为时窗的起始时间;wm为m号接收器上对应采样时间的波形值。
16、进一步地,所述动态调整纵波时差的搜索范围包括以下步骤:
17、根据实际地层情况确定纵波时差的初始搜索范围[s1,s2],进行某一深度点的纵波时差提取;
18、利用已处理井段的纵波时差计算结果动态调整纵波时差的搜索范围;
19、当已处理井段大于等于设定深度段d时,利用当前深度点前d深度段范围内纵波时差的计算结果计算当前深度点的纵波时差预测值;
20、将纵波时差初始搜索范围中心点移动到当前深度点的纵波时差预测值,得到调整后的纵波时差搜索范围[s′1,s′2]。
21、进一步地,当已处理井段小于设定深度段d时,所述地层纵波时差初始搜索范围[s1,s2]内拾取相关系数的最大值对应的时差。
22、进一步地,所述计算当前深度点的纵波时差预测值采用如下公式:
23、
24、式中,为当前深度点的纵波时差预测值;s(i)为第i个深度点的纵波时差值;n为设定深度段d范围内波形的深度采样点数,s(n)为设定深度段d范围内最后一个深度点的时差,s(n-1)为设定深度段d范围内倒数第二个深度点的时差。
25、进一步地,当已处理井段的深度大于等于设定深度段时,确定所述地层纵波时差包括以下步骤:
26、确定调整后的纵波时差搜索范围[s′1,s′2]内相关系数ρ(s)的最大值ρ(sm)和对应的模式波时差sm、[s′1,s′2]内相关系数ρ(s)的极大值ρ(se)和对应的模式波时差se,分别计算sm、se与的相关误差re1和re2;
27、设置相关系数ρ(s)的阈值为e1,纵波时差计算结果和纵波时差预测值的相关误差阈值为e2,若re2≤e2且ρ(se)≥e1时,取se作为当前深度的地层纵波时差;反之,取sm作为当前深度的地层纵波时差。
28、另一方面本专利技术还提出了一种时差搜索范围动态调整的地层纵波时差提取系统,所述提取系统包括:
29、纵波到时确定单元,用于分析1号接收器波形的能量比值随波形记录时间的变化,确定最大能量比值对应的时间作为1号接收器的纵波到时;
30、波形相关系数计算单元,用于计算根据所述1号接收器的纵波到时,计算时差域内不同时差对应的波形相关系数;
31、纵波时差调整单元,用于根据所述不同时差对应的波形相关系数动态调整纵波时差的搜索范围,在搜索范围内拾取波形相关系数的最大值对应的时差作为地层的纵波时差;
32、循环单元,用于重复上述步骤,直到完成整个测量井段的纵波时差提取。
33、进一步地,所述提取系统还包括:
34、预处理单元,用于在分析1号接收器波形的能量比值随波形记录时间的变化之前,在目的层段利用多极子阵列声波测井仪器进行井下测量,并收集全波测量数据对全波测量数据进行预处理;
35、对所述全波测量数据进行预处理包括对全波测量数据进行深度偏移校正、波形开始时间校正、波形去增益和数字滤波。
36、进一步地,所述纵波到时确定单元中调用以下公式确定最大能量比值对应的时间:
37、
38、式中,r(i)为第i个时间采样点处的能量比值;wk为第k个时间采样点的波形值;n1为前时本文档来自技高网...
【技术保护点】
1.一种时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,所述提取方法包括:
2.根据权利要求1所述的时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,
3.根据权利要求1或2所述的时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,
4.根据权利要求3所述的时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,
5.根据权利要求1所述时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,
6.根据权利要求1所述时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,
7.根据权利要求1所述的时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,
8.根据权利要求7所述的时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,
9.根据权利要求7所述的时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,
10.根据权利要求7所述的时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,
11.一种时差搜索范围动态调整的地层纵波时差提取系
12.根据权利要求11所述的时差搜索范围动态调整的地层纵波时差提取系统,其特征在于,所述提取系统还包括:
13.根据权利要求11或12所述的时差搜索范围动态调整的地层纵波时差提取系统,其特征在于,
14.根据权利要求11或12所述的时差搜索范围动态调整的地层纵波时差提取系统,其特征在于,
...【技术特征摘要】
1.一种时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,所述提取方法包括:
2.根据权利要求1所述的时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,
3.根据权利要求1或2所述的时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,
4.根据权利要求3所述的时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,
5.根据权利要求1所述时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,
6.根据权利要求1所述时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,
7.根据权利要求1所述的时差搜索范围动态调整的地层纵波时差提取方法,其特征在于,
8.根据权利要求7...
【专利技术属性】
技术研发人员:王浩,王才志,刘英明,夏守姬,原野,傅海成,
申请(专利权)人:中国石油天然气集团有限公司,
类型:发明
国别省市:
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。