显小断裂和窄河道等小规模的地质体边界。
[0067] 本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤可以通 过程序来指令相关的硬件来完成，该程序可以存储于一计算机可读取存储介质中，比如 R0M/RAM、磁碟、光盘等。
[0068] 以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详 细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保 护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本 发明的保护范围之内。
【主权项】
1. 一种地质体边界识别方法，其特征在于，所述的地质体边界识别方法包括： 设置一长度L为2M+1的滑动窗口，并设置长度L的阶值为N，获取2M+1个采样点，并根 据所述2M+1个采样点计算高斯窗函数G(n)，M为自然数； 利用倾角扫描方法确定每个所述采样点所在的局部反射界面的最佳倾角对； 分别将每一所述采样点作为分析样点，并以所述分析样点为中心，根据所述分析样点 的最佳倾角对从多个方向上分别读取2M+1个点的振幅值，生成多个地震信号序列Sk(n)，k =1，2, ? ? ?，i，n= 0, 1，? ? ?，2M，i= 2M+1 ; 根据所述地震信号序列&(1〇及高斯窗函数G(n)计算新地震序列Gsk(n); 计算所述新地震序列Gsk(n)的广义希尔伯特变换值GHk; 根据所述广义希尔伯特变换值GHk生成所述分析样点的边界检测值； 根据多个所述的边界检测值生成一数据体，并对所述数据体进行沿层切片分析，识别 地质体边界。2. 根据权利要求1所述的地质体边界识别方法，其特征在于，利用倾角扫描方法确定 每个所述采样点所在的局部反射界面的最佳倾角对，包括： 设置一个以米样点（x，y)为中心点、长短轴分别为a、b的横向窗口； 通过所述横向窗口获取所述采样点的地震参数； 将所述地震参数代入采样点相干计算公式，获取所述采样点的相干值； 分别设置线方向、道方向地层视倾角的极值P_min，P_max，Q_min，Q_max; 按一定的步长AP、AQ确定地质体反射界面的视倾角P、Q; 根据奈奎斯特采样定律获取np*nq个倾角对（P，Q)及对应的相干值c(t，P，Q)，其中，np =(P_max_P_min)/AP，nq= (Q_max-Q_min)/AQ； 以视倾角对（PbQM)为中心点选取C*d个视倾角（Pi，QJ及其对应的相干值c(t，Pi，QJ拟合成一个曲面G(P，Q)，其中，P1=PL±c*AP，Qn=QM±d*AP; 所述曲面G(P，Q)的最大值对应的倾角对（P，Q)为所述采样点的最佳倾角对。3. 根据权利要求1所述的地质体边界识别方法，其特征在于，所述的高斯窗函数为：其中，a为标准偏差的倒数。4. 根据权利要求1所述的地质体边界识别方法，其特征在于，所述的新地震序列 Gsk(n) =Sk (n)*G(n),k= 1, 2,. . . ,i〇5. 根据权利要求1所述的地质体边界识别方法，其特征在于，根据所述广义希尔伯特 变换值GHk生成所述分析样点的边界检测值，包括： 计算所述广义希尔伯特变换值6&的虚部的绝对值IGHk; 判断多个所述绝对值IGHk中的最大值，并将所述最大值作为所述分析样点的边界检测 值。6. -种地质体边界识别装置，其特征在于，所述的地质体边界识别装置包括： 高斯窗函数计算单元，用于设置一长度L为2M+1的滑动窗口，并设置长度L的阶值N， 获取2M+1个采样点，并根据所述2M+1个采样点计算高斯窗函数G(n)，M为自然数； 最佳倾角对确定单元，用于利用倾角扫描装置确定每个所述采样点所在的局部反射界 面的最佳倾角对； 地震信号序列生成单元，用于分别将每一所述采样点作为分析样点，并以所述分析样 点为中心，根据所述分析样点的最佳倾角对从多个方向上分别读取2M+1个点的振幅值，生 成多个地震f曰可序列Sk (n)，k= 1，2，? ? ?，i，n= 0，1，? ? ?，2M，i= 2M+1 ; 新地震序列计算单元，用于根据所述地震信号序列&(1〇及高斯窗函数G(n)计算新地 震序列Gsk(n); 广义希尔伯特变换值计算单元，用于计算所述新地震序列Gsk(n)的广义希尔伯特变换 值GHk; 边界检测值生成单元，用于根据所述广义希尔伯特变换值6氏生成所述分析样点的边 界检测值； 地质体边界识别单元，用于根据多个所述的边界检测值生成一数据体，并对所述数据 体进行沿层切片分析，识别地质体边界。7. 根据权利要求6所述的地质体边界识别装置，其特征在于，所述最佳倾角对确定单 元具体用于： 设置一个以米样点（x，y)为中心点、长短轴分别为a、b的横向窗口； 通过所述横向窗口获取所述采样点的地震参数； 将所述地震参数代入采样点相干计算公式，获取所述采样点的相干值； 分别设置线方向、道方向地层视倾角的极值P_min，P_max，Q_min，Q_max; 按一定的步长AP、AQ确定地质体反射界面的视倾角P、Q; 根据Nyquist采样定律获取np*nq个倾角对（P, Q)及对应的相干值c (t, P, Q)，其中，np=(P_max_P_min)/AP，nq= (Q_max-Q_min)/A Q； 以视倾角对（PbQM)为中心点选取C*d个视倾角（Pi，QJ及其对应的相干值c(t，Pi，QJ拟合成一个曲面G(P，Q)，其中，P1=PL±c*AP，Qn=QM±d*AP; 所述曲面G(P，Q)的最大值对应的倾角对（P，Q)为所述采样点的最佳倾角对。8. 根据权利要求6所述的地质体边界识别装置，其特征在于，所述的高斯窗函数为：其中，a为标准偏差的倒数。9. 根据权利要求6所述的地质体边界识别装置，其特征在于，所述的新地震序列 Gsk(n) =Sk (n)*G(n),k= 1, 2,. . . ,i〇10. 根据权利要求6所述的地质体边界识别装置，其特征在于，所述边界检测值生成单 元具体用于： 计算所述广义希尔伯特变换值6&的虚部的绝对值IGHk; 判断多个所述绝对值IGHk中的最大值，并将所述最大值作为所述分析样点的边界检测 值。
【专利摘要】本发明提供一种地质体边界识别方法及装置，该地质体边界识别方法包括：设置一长度L为2M+1的滑动窗口，并设置长度L的阶值为N，获取2M+1个采样点，并根据2M+1个采样点计算高斯窗函数G(n)；利用倾角扫描方法确定每个采样点所在的局部反射界面的最佳倾角对；分别将每一采样点作为分析样点，并以分析样点为中心，根据分析样点的视倾角从多个方向上分别读取2M+1个点的振幅值，生成多个地震信号序列Sk(n),k＝1,2,...,i，n＝0,1,...,2M，i＝2M+1；根据地震信号序列Sk(n)及高斯窗函数G(n)计算新地震序列GSk(n)；计算新地震序列GSk(n)的广义希尔伯特变换值GHk；根据广义希尔伯特变换值GHk生成分析样点的边界检测值；根据多个边界检测值生成一数据体，并对数据体进行沿层切片分析，识别地质体边界。
【IPC分类】G01V1/30, G01V1/28
【公开号】CN105093300
【申请号】CN201510445467
【发明人】孙勤华, 潘建国, 张虎权, 王宏斌, 孙东, 刘晓梅
【申请人】中国石油天然气股份有限公司
【公开日】2015年11月25日
【申请日】2015年7月27日