[发明专利]一种地震和测井时深关系自动转换的方法

摘要

一种地震和测井时深关系自动转换的方法，采集井A位置的地震子波，在井A周围选取已知井B位置地震子波，将井B的地震子波作为A井的地震子波，利用井A声波时差计算初始的时深关系，把井A测井数据转换到时间域，计算合成地震数据和整体时间移动大小，合成地震数据和对应时窗内的合成地震数据之间的相关系数。本发明可得到更加准确、合理的时深关系。

主权项

1.一种地震和测井时深关系自动转换的方法，特点是提供以下具体实施步骤：1）采集探区地震勘探数据，获得井A位置的地震子波；在井A周围选取已知测井资料的井B，提取井B位置地震子波；将井B的地震子波作为A井的地震子波，长度为WL；2）利用井A声波时差测井数据按照下式计算第i个采样点处深度d_i所对应的时间t_i，得到初始的时深关系；d_i＝d₁+(i-1)×Δd    i＝2，…，nti=t1+Σj=2iΔd×DTi]]>其中：t₁是第一个采样点d₁对应的时间，由下式得到：t₁=d₁×DT₁；d₁是在第一个采样点的深度；DT₁是第一个采样点的声波时差，单位为秒/米；i是采样点数；d_i是在第i个采样点对应的深度；DT_i是第i个采样点的声波时差；n是测井数据中总的采样点数；△d是测井数据的采样间隔；式中深度、时间和声波时差的单位均为米制标准单位；3）根据时深关系把井A测井数据转换到时间域，由下式计算反射系数R_i：Ri=ρi+1vi+1-ρiviρi+1vi+1+ρivi]]>其中：ρ和v分别为速度和密度，单位为米制标准单位；i为采样点数；4）由下式计算合成地震数据A_s：A_s＝R_i*W其中W是步骤1）提取的井B地震子波，R_i是步骤3）得到的反射系数；5）计算整体时间移动大小；取一个时窗长度WIN₁,这个时窗长度小于测井数据的实际长度，并且满足：Wsh≦WIN₁≦TS₂-TS₁+Wsh；其中：Wsh是时间值；TS₁和TS₂是测井中根据地质分层得到的目标储层时间范围；TS₁-Wsh至TS₂+Wsh是合成地震数据的时间范围；然后，①从合成地震数据A_S中取以TS₁-Wsh+WIN₁/2为中心，长度为WIN₁时窗内的数据，记为As₁¹；把As₁¹以△t为采样间隔，等于地震数据的采样间隔，在地震数据的范围内逐渐上、下移动，用下式计算移到位置处的观测地震数据Ao₁¹和As₁¹的相关系数c：c=Cor(Ao11,As11)]]>记录下相关系数c取得最大值对应的移动时间T_shift¹；②将合成地震数据的时窗中点下移△t，重复步骤①，得到对应的T_shift²；③重复步骤②，直到时窗的中点移动到TS₂+Wsh-WIN₁/2，且时窗在有测井数据的范围以内，得到各个位置处的T_shiftⁱ；④计算T_shift的概率密度函数；⑤最大概率密度函数对应的T_shift^max，确定为需要计算的整体时间移动大小；⑥将合成地震数据整体移动得到新的时深关系；6）用下式计算合成地震数据A_s2和对应时窗内的合成地震数据A_o2之间的相关系数c_old：c_old=Cor(A_o2,A_s2)式中：A_s2是TS₁-Wsh+T_shift^max至TS₂+Wsh+T_shift^max时间内的合成地震数据；然后，利用下面的步骤局部调整时深关系：①从合成地震数据A_s2中随机地选取一个时间点，把该时间点对应的采样点上移或者下移一个采样间隔的距离，得到另一合成地震数据A_s2’，计算新相关系数c_new；把迭代次数k的值增加1；②若新相关系数c_new大于相关系数c_old，则用另一合成地震数据A_s2’替代合成地震数据A_s2，并且让相关系数c_old等于新相关系数c_new；跳到步骤④；否则，进入下一步；③随机产生一个0到1之间的数；如果这个数大于那么用另一合成地震数据A_s2’替代合成地震数据A_s2，并且让相关系数c_old等于新相关系数c_new；否则，保留移动前的合成地震数据A_s2；其中e为自然对数的底，w为引入数学变量；④如果迭代次数k是w递减迭代次数k_L的整数倍，那么用w的递减因子l和w的积替代w；⑤如果从新相关系数c_new小于相关系数的最大值c_m并且k小于迭代次数的最大值k_m，那么返回到步骤①；否则，进入下一步；⑥完成地震和测井时深关系自动转换。