CN111965707A - 一种复杂构造含逆掩断裂的地震反演储层预测方法 - Google Patents

一种复杂构造含逆掩断裂的地震反演储层预测方法 Download PDF

Info

Publication number
CN111965707A
CN111965707A CN202010907578.0A CN202010907578A CN111965707A CN 111965707 A CN111965707 A CN 111965707A CN 202010907578 A CN202010907578 A CN 202010907578A CN 111965707 A CN111965707 A CN 111965707A
Authority
CN
China
Prior art keywords
seismic
module
geological
inversion
detection
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.)
Pending
Application number
CN202010907578.0A
Other languages
English (en)
Inventor
方光建
王锦喜
郭晓中
余川
曾春林
焦伟伟
余忠樯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing Institute of Geology and Mineral Resources
Original Assignee
Chongqing Institute of Geology and Mineral Resources
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing Institute of Geology and Mineral Resources filed Critical Chongqing Institute of Geology and Mineral Resources
Priority to CN202010907578.0A priority Critical patent/CN111965707A/zh
Publication of CN111965707A publication Critical patent/CN111965707A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48Processing data
    • G01V1/50Analysing data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616Data from specific type of measurement
    • G01V2210/6169Data from specific type of measurement using well-logging

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明属于地震预测技术领域,公开了一种复杂构造含逆掩断裂的地震反演储层预测方法,所述复杂构造含逆掩断裂的地震反演储层预测系统包括:地层探测模块、岩性探测模块、地质构造图绘制模块、中央控制模块、地质模型构建程序、历史信息获取模块、地震子波采集模块、地震信息合成模块、地震反演模块、分析模块、地震预测模块。本发明进行地层和岩性的探测,绘制地质构造图以及依照地形构造图建立地形三维模型,能实现对预测区域的地形、地质的真实展现,方便进行地震反演模拟,实现含逆掩断层上下盘反演计算;通过对地震历史数据的获取实现对地震预测准确性的提升;预测方法简单,能够实现对不同地质区域的地震预测,预测结果准确。

Description

一种复杂构造含逆掩断裂的地震反演储层预测方法
技术领域
本发明属于地震预测技术领域,尤其涉及一种复杂构造含逆掩断裂的地震反演储层预测方法。
背景技术
目前:地震记录是震源产生的地震子波传播到地下遇到反射界面反射回地震后被记录的结果。通常认为地震波在地下介质中传播可以简化成Wilson模型方程,即地震记录是地震子波与反射系数序列的褶积结果。地震波在地下介质中传播的速度受介质的密度、孔隙度、流体类型、流体含量等物性参数影响,地下介质的作用就相当于一个滤波器,因此,地面接收的地震信号是地下介质的综合响应。
地震反演是根据反射地震记录求取地下介质物性参数的反问题求解的统称。目前地震反演方法主要有道积分、稀疏脉冲反演、基于模型的反演、随机反演等算法。但这些反演计算都是在正断层或无断层的模型下计算的,无法完成含逆断层的反演计算。然而,随着地震勘探开发的发展,现在已经过渡到复杂构造勘探阶段,含逆掩断裂的地质构造已经成为勘探开发中非常普遍的地质现象。但现有技术中的反演方法无法实现含逆掩断层上下盘反演计算,且进行多次反演的结果不准确,无法满足科研生产的需要。
通过上述分析,现有技术存在的问题及缺陷为:现有技术中的反演方法无法实现含逆掩断层上下盘反演计算,且进行多次反演的结果不准确,无法满足科研生产的需要。
发明内容
针对现有技术存在的问题,本发明提供了一种复杂构造含逆掩断裂的地震反演储层预测方法。
本发明是这样实现的,一种复杂构造含逆掩断裂的地震反演储层预测系统,所述复杂构造含逆掩断裂的地震反演储层预测系统包括:
地层探测模块、岩性探测模块、地质构造图绘制模块、中央控制模块、地质模型构建程序、历史信息获取模块、地震子波采集模块、地震信息合成模块、地震反演模块、分析模块、地震预测模块;
地层探测模块,与中央控制模块连接,用于通过地层探测程序进行地层探测,得到不同探测点不同深度对应的地层信息;
岩性探测模块,与中央控制模块连接,用于通过岩性探测程序进行岩性探测,得到不同地层对应的岩性;
地质构造图绘制模块,与中央控制模块连接,用于通过地质构造图绘制程序进行地质构造剖面图的绘制,展示地层、岩性信息;
中央控制模块,与地层探测模块、岩性探测模块、地质构造图绘制模块、地质模型构建模块、历史信息获取模块、地震子波采集模块、地震信息合成模块、地震反演模块、分析模块、地震预测模块连接,用于通过主控机控制各个模块正常运行;
地质模型构建模块,与中央控制模块连接,用于通过地质模型构建程序依照绘制的地质构造剖面图进行地质模型的构建,得到三维地质模型;
历史信息获取模块,与中央控制模块连接,用于通过历史信息获取模块仓数据库中获取预测区域地震的历史信息,包括地震的时间、震源、震级;
地震子波采集模块,与中央控制模块连接,用于通过地震子波采集程序采集地震子波;
地震信息合成模块,与中央控制模块连接,用于通过地震信息合成模块依据地震波传播的褶积数学物理模型,用测井曲线制作校正合成地震信息;
地震反演模块,与中央控制模块连接,用于通过地震反演程序在建立的地质模型中进行地震反演;
分析模块,与中央控制模块连接,用于通过分析程序进行地震反演结果的分析,得到不同地质对地震的影响;
地震预测模块,与中央控制模块连接,用于通过地震预测程序预测地震概率、时间以及地震强度。
进一步,所述地质模型构建模块包括:
地质数据输入单元,用于输入地层探测模块和岩性探测模块采集的数据并确定地质三维建模的有效数据;
参数定义单元,用于为地质建模对象的地质属性参数进行定义;
边界线构建单元,用于根据地质数据和定义参数建立地质对象的边界线;
综合模型构建单元,用于通过对地质对象的边界线进行耦合和运算构建地质综合模型。
进一步,所述地震子波采集模块包括:
参数采集单元,用于获取时间域的地震子波、最大频率、最小频率、采样点数以及采样间隔;
转换单元,用于将时间域地震子波转化至频率域地震子波;
校正单元,用于基于频率域地震子波,获取单频校正因子,并进行频率域地震子波的校正;
输出单元,用于基于校正后的频率域地震子波获取估计地震子波。
本发明的另一目的在于提供一种复杂构造含逆掩断裂的地震反演储层预测方法,所述复杂构造含逆掩断裂的地震反演储层预测方法包括以下步骤:
步骤一,通过地层探测模块利用地层探测程序进行地层探测,得到不同探测点不同深度对应的地层信息;通过岩性探测模块利用岩性探测程序进行岩性探测,得到不同地层对应的岩性;
步骤二,通过地质构造图绘制模块利用地质构造图绘制程序进行地质构造剖面图的绘制,展示地层、岩性信息;通过地质模型构建模块利用地质模型构建程序依照绘制的地质构造剖面图进行地质模型的构建,得到三维地质模型;
步骤三,通过历史信息获取模块利用历史信息获取模块仓数据库中获取预测区域历史地震的时间、震源、震级;
步骤四,通过地震子波采集模块利用地震子波采集程序采集地震子波;通过地震信息合成模块利用地震信息合成模块依据地震波传播的褶积数学物理模型,用测井曲线制作合成地震信息;
步骤五,通过地震反演模块利用地震反演程序在建立的地质模型中进行地震反演;
步骤六,通过分析模块利用分析程序进行地震反演结果的分析,得到不同地质对地震的影响;通过地震预测模块利用地震预测程序预测地震概率、时间以及地震强度。
进一步,步骤一中,所述通过地层探测模块利用地层探测程序进行地层探测,得到不同探测点不同深度对应的地层信息,具体包括:
步骤A,在探测器上设置由多个发射线圈和接收线圈组成的线圈系,所述多个发射线圈发射不同频率的探测信号;
步骤B,调节所述多个发射线圈与所述接收线圈的距离,使预设深度的地层对所述每个频率的探测信号的响应信号具有相同的预设响应参数值,得到所述预设深度对应的线圈系的优化位置参数;
步骤C,根据要探测的目标深度和对应的所述优化位置参数设置所述线圈系,并向所述目标深度的地层发射所述不同频率的探测信号;
步骤D,测量所述目标深度的地层对所述探测信号的响应信号,得到目标深度的地层信息;
步骤E,进行目标深度的变换,得到不同深度的地层信息。
进一步,步骤B中,所述调节所述多个发射线圈与所述接收线圈的距离,使预设深度的地层对所述每个频率的探测信号的响应信号具有相同的预设响应参数值,包括:
获取不同频率的探测信号的伪几何因子曲线;确定所述预设深度和所述预设响应参数值对应的目标点;调节所述多个发射线圈与所述接收线圈的距离,使得所述不同频率的探测信号的伪几何因子曲线都经过所述目标点;将所述发射线圈与所述接收线圈的距离记为相应频率下所述预设深度对应的线圈系的优化位置参数。
进一步,步骤二中,所述通过地质构造图绘制模块利用地质构造图绘制程序进行地质构造剖面图的绘制,具体为:
基于探测的地层和岩性信息建立协同克立格方程组矩阵;根据所述协同克立格方程组矩阵建立二维地质剖面图。
进一步,所述根据所述协同克立格方程组矩阵建立二维地质剖面图,包括以下步骤:
1)基于采集的地质、岩性数据建立地质属性数据库;
2)定义地质对象之间的时空关系;
3)根据所述地质属性数据库和所述时空关系采用协同克立格方程组矩阵构建第一势场;
4)根据二维地质剖面图的位置和所述第一势场构建第一地质边界线;
5)根据所述第一地质边界线的地质属性确定二维地质剖面图。
进一步,步骤二中,所述通过地质模型构建模块利用地质模型构建程序依照绘制的地质构造剖面图进行地质模型的构建,得到三维地质模型,具体包括以下步骤:
(1)采用网格索引与包围盒碰撞检测算法优化欧式距离;
(2)根据优化后的所述欧式距离,采用协同克立格方程组矩阵构建第二势场;
(3)根据所述第二势场和三维地质模型剖面的位置构建第二地质边界线;
(4)根据所述第二地质边界线的地质属性确定三维地质模型剖面图。
进一步,步骤(3)中,所述根据所述第二势场和三维地质模型剖面的位置构建第二地质边界线,具体为:
先建三角网格,然后在所述三角网格上根据所述第二势场和三维地质剖面的位置追踪获得第二地质边界线。
进一步,步骤五中,所述通过地震反演模块利用地震反演程序在建立的地质模型中进行地震反演,具体为:
进行地震资料多尺度分解;针对每个尺度或频率内的地震反射记录,开展地震反演。
进一步,步骤六中,所述通过分析模块利用分析程序进行地震反演结果的分析,得到不同地质对地震的影响,包括:
确定频变加权系数,并在此基础上提取频变粘弹性流体因子;基于频变粘弹性流体因子,对流体的分布状态进行预测,得到地震预测信息。
结合上述的所有技术方案,本发明所具备的优点及积极效果为:本发明进行地层和岩性的探测,绘制地质构造图以及依照地形构造图建立地形三维地质模型,能实现对预测区域的地形、地质的真实展现,方便进行地震反演模拟,实现含逆掩断层上下盘反演计算;并且本申请构建的地质三维地质模型能够准确体现该预测区域的地质情况,进行地震反演得到的结构更准确;通过对地震历史数据的获取以及融合,能够实现对地震预测准确性的提升。本发明提供的预测方法简单,能够实现对不同地质区域的地震预测,预测结果准确。
附图说明
为了更清楚地说明本申请实施例的技术方案,下面将对本申请实施例中所需要使用的附图做简单的介绍,显而易见地,下面所描述的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的复杂构造含逆掩断裂的地震反演储层预测系统的结构框图。
图中:1、地层探测模块;2、岩性探测模块;3、地质构造图绘制模块;4、中央控制模块;5、地质模型构建程序;6、历史信息获取模块;7、地震子波采集模块;8、地震信息合成模块;9、地震反演模块;10、分析模块;11、地震预测模块。
图2是本发明实施例提供的复杂构造含逆掩断裂的地震反演储层预测方法的流程图。
图3是本发明实施例提供的通过地层探测模块利用地层探测程序进行地层探测,得到不同探测点不同深度对应的地层信息的流程图。
图4是本发明实施例提供的根据所述协同克立格方程组矩阵建立二维地质剖面图的流程图。
图5是本发明实施例提供的通过地质模型构建模块利用地质模型构建程序依照绘制的地质构造剖面图进行地质模型的构建,得到三维地质模型的流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
针对现有技术存在的问题,本发明提供了一种复杂构造含逆掩断裂的地震反演储层预测方法,下面结合附图对本发明作详细的描述。
如图1所示,本发明实施例提供的复杂构造含逆掩断裂的地震反演储层预测系统包括:
地层探测模块1、岩性探测模块2、地质构造图绘制模块3、中央控制模块4、地质模型构建程序5、历史信息获取模块6、地震子波采集模块7、地震信息合成模块8、地震反演模块9、分析模块10、地震预测模块11;
地层探测模块1,与中央控制模块4连接,用于通过地层探测程序进行地层探测,得到不同探测点不同深度对应的地层信息;
岩性探测模块2,与中央控制模块4连接,用于通过岩性探测程序进行岩性探测,得到不同地层对应的岩性;
地质构造图绘制模块3,与中央控制模块4连接,用于通过地质构造图绘制程序进行地质构造剖面图的绘制,展示地层、岩性信息;
中央控制模块4,与地层探测模块1、岩性探测模块2、地质构造图绘制模块3、地质模型构建程序5、历史信息获取模块6、地震子波采集模块7、地震信息合成模块8、地震反演模块9、分析模块10、地震预测模块11连接,用于通过主控机控制各个模块正常运行;
地质模型构建模块5,与中央控制模块4连接,用于通过地质模型构建程序依照绘制的地质构造剖面图进行地质模型的构建,得到三维地质模型;
历史信息获取模块6,与中央控制模块4连接,用于通过历史信息获取模块仓数据库中获取预测区域地震的历史信息,包括地震的时间、震源、震级;
地震子波采集模块7,与中央控制模块4连接,用于通过地震子波采集程序采集地震子波;
地震信息合成模块8,与中央控制模块4连接,用于通过地震信息合成模块依据地震波传播的褶积数学物理模型,用测井曲线制作合成地震信息;
地震反演模块9,与中央控制模块4连接,用于通过地震反演程序在建立的地质模型中进行地震反演;
分析模块10,与中央控制模块4连接,用于通过分析程序进行地震反演结果的分析,得到不同地质对地震的影响;
地震预测模块11,与中央控制模块4连接,用于通过地震预测程序预测地震概率、时间以及地震强度。
本发明实施例中的地质模型构建模块包括:
地质数据输入单元,用于输入地层探测模块和岩性探测模块采集的数据并确定地质三维建模的有效数据;
参数定义单元,用于为地质建模对象的地质属性参数进行定义;
边界线构建单元,用于根据地质数据和定义参数建立地质对象的边界线;
综合模型构建单元,用于通过对地质对象的边界线进行耦合和运算构建地质综合模型。
本发明实施例中的地震子波采集模块包括:
参数采集单元,用于获取时间域的地震子波、最大频率、最小频率、采样点数以及采样间隔;
转换单元,用于将时间域地震子波转化至频率域地震子波;
校正单元,用于基于频率域地震子波,获取单频校正因子,并进行频率域地震子波的校正;
输出单元,用于基于校正后的频率域地震子波获取估计地震子波。
如图2所示,本发明实施例提供的复杂构造含逆掩断裂的地震反演储层预测方法包括以下步骤:
S101,通过地层探测模块利用地层探测程序进行地层探测,得到不同探测点不同深度对应的地层信息;通过岩性探测模块利用岩性探测程序进行岩性探测,得到不同地层对应的岩性;
S102,通过地质构造图绘制模块利用地质构造图绘制程序进行地质构造剖面图的绘制,展示地层、岩性信息;通过地质模型构建模块利用地质模型构建程序依照绘制的地质构造剖面图进行地质模型的构建,得到三维地质模型;
S103,通过历史信息获取模块利用历史信息获取模块仓数据库中获取预测区域历史地震的时间、震源、震级;
S104,通过地震子波采集模块利用地震子波采集程序采集地震子波;通过地震信息合成模块利用地震信息合成模块依据地震波传播的褶积数学物理模型,用测井曲线制作合成地震信息;
S105,通过地震反演模块利用地震反演程序在建立的地质模型中进行地震反演;
S106,通过分析模块利用分析程序进行地震反演结果的分析,得到不同地质对地震的影响;通过地震预测模块利用地震预测程序预测地震概率、时间以及地震强度。
如图3所示,步骤S101中,本发明实施例提供的通过地层探测模块利用地层探测程序进行地层探测,得到不同探测点不同深度对应的地层信息,具体包括:
S201,在探测器上设置由多个发射线圈和接收线圈组成的线圈系,所述多个发射线圈发射不同频率的探测信号;
S202,调节所述多个发射线圈与所述接收线圈的距离,使预设深度的地层对所述每个频率的探测信号的响应信号具有相同的预设响应参数值,得到所述预设深度对应的线圈系的优化位置参数;
S203,根据要探测的目标深度和对应的所述优化位置参数设置所述线圈系,并向所述目标深度的地层发射所述不同频率的探测信号;
S204,测量所述目标深度的地层对所述探测信号的响应信号,得到目标深度的地层信息;
S205,进行目标深度的变换,得到不同深度的地层信息。
步骤S202中,本发明实施例提供的调节所述多个发射线圈与所述接收线圈的距离,使预设深度的地层对所述每个频率的探测信号的响应信号具有相同的预设响应参数值,包括:
获取不同频率的探测信号的伪几何因子曲线;确定所述预设深度和所述预设响应参数值对应的目标点;调节所述多个发射线圈与所述接收线圈的距离,使得所述不同频率的探测信号的伪几何因子曲线都经过所述目标点;将所述发射线圈与所述接收线圈的距离记为相应频率下所述预设深度对应的线圈系的优化位置参数。
步骤S102中,本发明实施例提供的通过地质构造图绘制模块利用地质构造图绘制程序进行地质构造剖面图的绘制,具体为:
基于探测的地层和岩性信息建立协同克立格方程组矩阵;根据所述协同克立格方程组矩阵建立二维地质剖面图。
如图4所示,本发明实施例提供的根据所述协同克立格方程组矩阵建立二维地质剖面图,包括以下步骤:
S301,基于采集的地质、岩性数据建立地质属性数据库;
S302,定义地质对象之间的时空关系;
S303,根据所述地质属性数据库和所述时空关系采用协同克立格方程组矩阵构建第一势场;
S304,根据二维地质剖面图的位置和所述第一势场构建第一地质边界线;
S305,根据所述第一地质边界线的地质属性确定二维地质剖面图。
如图5所示,步骤S102中,本发明实施例提供的通过地质模型构建模块利用地质模型构建程序依照绘制的地质构造剖面图进行地质模型的构建,得到三维地质模型,具体包括以下步骤:
S401,采用网格索引与包围盒碰撞检测算法优化欧式距离;
S402,根据优化后的所述欧式距离,采用协同克立格方程组矩阵构建第二势场;
S403,根据所述第二势场和三维地质模型剖面的位置构建第二地质边界线;
S404,根据所述第二地质边界线的地质属性确定三维地质模型剖面图。
步骤S403中,本发明实施例提供的根据所述第二势场和三维地质模型剖面的位置构建第二地质边界线,具体为:
先建三角网格,然后在所述三角网格上根据所述第二势场和三维地质剖面的位置追踪获得第二地质边界线。
步骤S105中,本发明实施例提供的通过地震反演模块利用地震反演程序在建立的地质模型中进行地震反演,具体为:
进行地震资料多尺度分解;针对每个尺度或频率内的地震反射记录,开展地震反演。
步骤S106中,本发明实施例提供的通过分析模块利用分析程序进行地震反演结果的分析,得到不同地质对地震的影响,包括:
确定频变加权系数,并在此基础上提取频变粘弹性流体因子;基于频变粘弹性流体因子,对流体的分布状态进行预测,得到地震预测信息。
以上所述,仅为本发明较优的具体的实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,都应涵盖在本发明的保护范围之内。

Claims (10)

1.一种复杂构造含逆掩断裂的地震反演储层预测方法,其特征在于,所述复杂构造含逆掩断裂的地震反演储层预测方法包括以下步骤:
步骤一,通过地层探测模块利用地层探测程序进行地层探测,得到不同探测点不同深度对应的地层信息;通过岩性探测模块利用岩性探测程序进行岩性探测,得到不同地层对应的岩性;
所述通过地层探测模块利用地层探测程序进行地层探测,得到不同探测点不同深度对应的地层信息,具体包括:
步骤A,在探测器上设置由多个发射线圈和接收线圈组成的线圈系,所述多个发射线圈发射不同频率的探测信号;
步骤B,调节所述多个发射线圈与所述接收线圈的距离,使预设深度的地层对所述每个频率的探测信号的响应信号具有相同的预设响应参数值,得到所述预设深度对应的线圈系的优化位置参数;
步骤C,根据要探测的目标深度和对应的所述优化位置参数设置所述线圈系,并向所述目标深度的地层发射所述不同频率的探测信号;
步骤D,测量所述目标深度的地层对所述探测信号的响应信号,得到目标深度的地层信息;
步骤E,进行目标深度的变换,得到不同深度的地层信息;
步骤二,通过地质构造图绘制模块利用地质构造图绘制程序进行地质构造剖面图的绘制,展示地层、岩性信息;通过地质模型构建模块利用地质模型构建程序依照绘制的地质构造剖面图进行地质模型的构建,得到三维地质模型;
步骤三,通过历史信息获取模块利用历史信息获取模块仓数据库中获取预测区域历史地震的时间、震源、震级;
步骤四,通过地震子波采集模块利用地震子波采集程序采集地震子波;通过地震信息合成模块利用地震信息合成模块依据地震波传播的褶积数学物理模型,用测井曲线制作校正合成地震信息;
步骤五,通过地震反演模块利用地震反演程序在建立的地质模型中进行地震反演;
步骤六,通过分析模块利用分析程序进行地震反演结果的分析,得到不同地质对地震的影响;通过地震预测模块利用地震预测程序预测地震概率、时间以及地震强度。
2.如权利要求1所述的复杂构造含逆掩断裂的地震反演储层预测方法,其特征在于,步骤B中,所述调节所述多个发射线圈与所述接收线圈的距离,使预设深度的地层对所述每个频率的探测信号的响应信号具有相同的预设响应参数值,包括:
获取不同频率的探测信号的伪几何因子曲线;确定所述预设深度和所述预设响应参数值对应的目标点;调节所述多个发射线圈与所述接收线圈的距离,使得所述不同频率的探测信号的伪几何因子曲线都经过所述目标点;将所述发射线圈与所述接收线圈的距离记为相应频率下所述预设深度对应的线圈系的优化位置参数。
3.如权利要求1所述的复杂构造含逆掩断裂的地震反演储层预测方法,其特征在于,步骤二中,所述通过地质构造图绘制模块利用地质构造图绘制程序进行地质构造剖面图的绘制,具体为:
基于探测的地层和岩性信息建立协同克立格方程组矩阵;根据所述协同克立格方程组矩阵建立二维地质剖面图;
所述根据所述协同克立格方程组矩阵建立二维地质剖面图,包括以下步骤:
1)基于采集的地质、岩性数据建立地质属性数据库;
2)定义地质对象之间的时空关系;
3)根据所述地质属性数据库和所述时空关系采用协同克立格方程组矩阵构建第一势场;
4)根据二维地质剖面图的位置和所述第一势场构建第一地质边界线;
5)根据所述第一地质边界线的地质属性确定二维地质剖面图。
4.如权利要求1所述的复杂构造含逆掩断裂的地震反演储层预测方法,其特征在于,步骤二中,所述通过地质模型构建模块利用地质模型构建程序依照绘制的地质构造剖面图进行地质模型的构建,得到三维地质模型,具体包括以下步骤:
(1)采用网格索引与包围盒碰撞检测算法优化欧式距离;
(2)根据优化后的所述欧式距离,采用协同克立格方程组矩阵构建第二势场;
(3)根据所述第二势场和三维地质模型剖面的位置构建第二地质边界线;
(4)根据所述第二地质边界线的地质属性确定三维地质模型剖面图。
5.如权利要求4所述的复杂构造含逆掩断裂的地震反演储层预测方法,其特征在于,步骤(3)中,所述根据所述第二势场和三维地质模型剖面的位置构建第二地质边界线,具体为:
先建三角网格,然后在所述三角网格上根据所述第二势场和三维地质剖面的位置追踪获得第二地质边界线。
6.如权利要求1所述的复杂构造含逆掩断裂的地震反演储层预测方法,其特征在于,步骤五中,所述通过地震反演模块利用地震反演程序在建立的地质模型中进行地震反演,具体为:
进行地震资料多尺度分解;针对每个尺度或频率内的地震反射记录,开展地震反演。
7.如权利要求1所述的复杂构造含逆掩断裂的地震反演储层预测方法,其特征在于,步骤六中,所述通过分析模块利用分析程序进行地震反演结果的分析,得到不同地质对地震的影响,包括:
确定频变加权系数,并在此基础上提取频变粘弹性流体因子;基于频变粘弹性流体因子,对流体的分布状态进行预测,得到地震预测信息。
8.一种实施权利要求1~7任意一项所述的复杂构造含逆掩断裂的地震反演储层预测方法的复杂构造含逆掩断裂的地震反演储层预测系统,其特征在于,所述复杂构造含逆掩断裂的地震反演储层预测系统包括:
地层探测模块、岩性探测模块、地质构造图绘制模块、中央控制模块、地质模型构建程序、历史信息获取模块、地震子波采集模块、地震信息合成模块、地震反演模块、分析模块、地震预测模块;
地层探测模块,与中央控制模块连接,用于通过地层探测程序进行地层探测,得到不同探测点不同深度对应的地层信息;
岩性探测模块,与中央控制模块连接,用于通过岩性探测程序进行岩性探测,得到不同地层对应的岩性;
地质构造图绘制模块,与中央控制模块连接,用于通过地质构造图绘制程序进行地质构造剖面图的绘制,展示地层、岩性信息;
中央控制模块,与地层探测模块、岩性探测模块、地质构造图绘制模块、地质模型构建程序、历史信息获取模块、地震子波采集模块、地震信息合成模块、地震反演模块、分析模块、地震预测模块连接,用于通过主控机控制各个模块正常运行;
地质模型构建程序,与中央控制模块连接,用于通过地质模型构建程序依照绘制的地质构造剖面图进行地质模型的构建,得到三维地质模型;
历史信息获取模块,与中央控制模块连接,用于通过历史信息获取模块仓数据库中获取预测区域地震的历史信息,包括地震的时间、震源、震级;
地震子波采集模块,与中央控制模块连接,用于通过地震子波采集程序采集地震子波;
地震信息合成模块,与中央控制模块连接,用于通过地震信息合成模块依据地震波传播的褶积数学物理模型,用测井曲线制作校正合成地震信息;
地震反演模块,与中央控制模块连接,用于通过地震反演程序在建立的地质模型中进行地震反演;
分析模块,与中央控制模块连接,用于通过分析程序进行地震反演结果的分析,得到不同地质对地震的影响;
地震预测模块,与中央控制模块连接,用于通过地震预测程序预测地震概率、时间以及地震强度。
9.如权利要求8所述的复杂构造含逆掩断裂的地震反演储层预测系统,其特征在于,所述地质模型构建模块包括:
地质数据输入单元,用于输入地层探测模块和岩性探测模块采集的数据并确定地质三维建模的有效数据;
参数定义单元,用于为地质建模对象的地质属性参数进行定义;
边界线构建单元,用于根据地质数据和定义参数建立地质对象的边界线;
综合模型构建单元,用于通过对地质对象的边界线进行耦合和运算构建地质综合模型。
10.如权利要求8所述的复杂构造含逆掩断裂的地震反演储层预测系统,其特征在于,所述地震子波采集模块包括:
参数采集单元,用于获取时间域的地震子波、最大频率、最小频率、采样点数以及采样间隔;
转换单元,用于将时间域地震子波转化至频率域地震子波;
校正单元,用于基于频率域地震子波,获取单频校正因子,并进行频率域地震子波的校正;
输出单元,用于基于校正后的频率域地震子波获取估计地震子波。
CN202010907578.0A 2020-09-02 2020-09-02 一种复杂构造含逆掩断裂的地震反演储层预测方法 Pending CN111965707A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010907578.0A CN111965707A (zh) 2020-09-02 2020-09-02 一种复杂构造含逆掩断裂的地震反演储层预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010907578.0A CN111965707A (zh) 2020-09-02 2020-09-02 一种复杂构造含逆掩断裂的地震反演储层预测方法

Publications (1)

Publication Number Publication Date
CN111965707A true CN111965707A (zh) 2020-11-20

Family

ID=73392243

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010907578.0A Pending CN111965707A (zh) 2020-09-02 2020-09-02 一种复杂构造含逆掩断裂的地震反演储层预测方法

Country Status (1)

Country Link
CN (1) CN111965707A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116821999A (zh) * 2023-05-11 2023-09-29 青海省地质调查院(青海省地质矿产研究院、青海省地质遥感中心) 一种板块尺度三维地质模型构建方法及系统

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116821999A (zh) * 2023-05-11 2023-09-29 青海省地质调查院(青海省地质矿产研究院、青海省地质遥感中心) 一种板块尺度三维地质模型构建方法及系统
CN116821999B (zh) * 2023-05-11 2024-01-02 青海省地质调查院(青海省地质矿产研究院、青海省地质遥感中心) 一种板块尺度三维地质模型构建方法及系统

Similar Documents

Publication Publication Date Title
CA2920499C (en) Stratigraphic function
AU2010236999B2 (en) Interferometric seismic data processing
US6928367B2 (en) Reservoir fracture characterization
EP2171499B1 (en) Method for determining seismic data quality
CN111596364B (zh) 基于高精度层序地层格架的地震沉积微相组合分析方法
NO327326B1 (no) Fremgangsmate og system for a bringe data i overensstemmelse mellom seismisk og bronnloggeskalaer i 3-D geologisk modellering
AU2015280679B2 (en) Instantaneous isochron attribute-based geobody identification for reservoir modeling
US7280952B2 (en) Well planning using seismic coherence
EA031663B1 (ru) Способ вычислительной геологии
EP3978961B1 (en) System and method for quantitative seismic integration modeling workflow
EA030770B1 (ru) Система и способ адаптивной сейсмической оптики
US20150109887A1 (en) Sonic adaptor for converting sonic or ultrasonic waveform data for use with a seismic-based computer program
CN111965707A (zh) 一种复杂构造含逆掩断裂的地震反演储层预测方法
Barison et al. Processing and interpretation of seismic reflection data from the Los Humeros super-hot geothermal system
CN111077578B (zh) 岩层分布预测方法和装置
CN111722284A (zh) 一种基于道集数据建立速度深度模型的方法
CN104111476A (zh) 构建地层速度场的方法及装置
Kruiver et al. Capturing spatial variability in the regional Ground Motion Model of Groningen, the Netherlands
CN113848593A (zh) 一种定量预测含煤地层中岩浆岩侵蚀区的方法
CN115880455A (zh) 基于深度学习的三维智能插值方法
US11022710B2 (en) Identifying anomalies in a subterranean formation based on seismic attributes
CN116520394A (zh) 基于地震测井双驱动融合多尺度信息直接预测孔隙度方法
CN113806674A (zh) 古河道纵向尺度的量化方法、装置、电子设备及存储介质
US11852768B2 (en) Multimodal approach to target stratigraphic plays through seismic sequence stratigraphy, rock physics, seismic inversion and machine learning
EP3669210B1 (en) Method for obtaining estimates of a model parameter so as to characterise the evolution of a subsurface volume over a time period using time-lapse seismic

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