CN1337583A - 基于地质地球物理模型的三维地震观测系统优化设计方法 - Google Patents
基于地质地球物理模型的三维地震观测系统优化设计方法 Download PDFInfo
- Publication number
- CN1337583A CN1337583A CN 01120445 CN01120445A CN1337583A CN 1337583 A CN1337583 A CN 1337583A CN 01120445 CN01120445 CN 01120445 CN 01120445 A CN01120445 A CN 01120445A CN 1337583 A CN1337583 A CN 1337583A
- Authority
- CN
- China
- Prior art keywords
- seismic
- imaging results
- geophysical model
- data
- design
- 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.)
- Granted
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明是一种基于地质地球物理模型的三维地震观测系统优化设计方法,包括以下步骤:A)建立先验的地质地球物理模型;B)设计两种或多种三维地震数据采集的观测系统;C)采用虚谱法三维波动方程数值模拟与三维地震物理模型模拟两种方法联合正演来完成数据采集工作;D)以共反射点为基础,对采集的数据进行叠前深度偏移处理,得到地震处理几种成像结果;E)对所述地震处理几种成像结果进行成像质量的综合对比评价,经过对比,确定其中一个最佳地震成像结果。本方法适用范围广泛,应用效果良好,对指导石油天然气地震勘探实际生产、提高油气勘探成功率具有现实意义。
Description
技术领域
本发明涉及石油天然气和其它矿产资源及地质工程中三维地震勘探的一种技术设计优化方法,特别是一种基于地质地球物理模型的三维地震观测系统优化设计方法。
背景技术
地震勘探三大环节是采集、处理与解释。采集是基础;地面上震源激发点位组成的激发系统和检波器接收点位组成的接收系统,必须通过正确的几何分布,即优化的地震观测系统,才能采集到有效的完整的地下地质构造信息。因此,观测系统在地震勘探中起重要作用,它的设计优劣直接关系到勘探任务的成败。又因为三维地震工作量大、投资高,三维工区往往布置在复杂构造地区,所以三维观测系统优化设计就倍为重要。
至今,国内外广为采用的三维观测系统设计软件多是以水平层状地质模型和共中心点(CMP)叠加原理为基础的,它们不适用于复杂构造、大倾角地层的三维地震技术设计。在复杂构造情况下,应用这类设计方案,不仅不能得到高质量的地震成像,甚至可能得出错误的地震成像结果。
发明内容
因此,本发明的目的使提供一种基于地质地球物理模型的三维地震观测系统优化设计方法。
该方法以勘探地区先验的三维地质和地球物理模型为基础,采用数值模拟与物理模拟联合正演技术,根据共反射点(CRP)原理,以地质目标高质量地震成像为评价标准,将三维地震采集技术设计方案与地震数据处理的成像质量直接联系起来,以此实现优化三维观测系统设计这一采集技术目标。
本发明提供的基于地质地球物理模型的三维地震观测系统优化设计方法包括以下步骤:
A)根据勘探地区现有的勘探成果,即复杂地质构造和复杂地质体形态及其地球物理参数建立先验的地质地球物理模型;B)以上述地质地球物理模型为勘探目标,设计两种或多种三维地震数据采集的观测系统;C)以该模型作为三维地震观测对象,按照观测系统设计所要求的技术参数和施工图形在实验室模拟野外地震数据采集过程,获得模拟的野外原始地震记录,具体作法是采用虚谱法三维波动方程数值模拟与三维地震物理模型模拟两种方法联合正演来完成这一数据采集工作;D)对采集的数据进行三维地震数据处理,其中以共反射点为基础,对数据进行叠前深度偏移处理,得到地震处理几种成像结果;E)以先验的地质地球物理模型作为评价与对比标准,以信噪比、分辨率和构造形态特征作为技术分析主要内容,对地震处理几种成像结果进行成像质量的综合对比评价,经过对比,确定其中一个最佳地震成像结果。
此外,上述步骤E)还包括:如果这个成像结果还没有达到地震勘探的目标要求,则修改原先的三维观测系统设计方案,随后重新进行三维地震数据采集处理工作;如果这个成像结果已经达到地震勘探的目标要求,则输出三维地震观测系统优化设计流程图。
下面结合附图对本发明进行详细说明,以便于进一步理解本发明的上述目的、特点以及由此带来的优点。
附图说明
图1是本发明的基于地质地球物理模型的三维地震观测系统优化设计方法的流程图;
图2(a)是根据先验资料制作的千米桥潜山构造物理模型;
图2(b)是南北方向穿过潜山顶部的二维337测线地质剖面图;
图3是在337测线上分别采用虚谱法三维波动方程正演模拟和采用三维地震物理模拟所得到的地震记录;
图4是三种观测系统采集的地震数据水平叠加成像结果(337线),其中图4(a)是8线2炮纵24次横1次观测系统;图4(b)是8线4炮纵24次横2次观测系统;图4(c)是8线8炮纵24次横2次观测系统;
图5是不同宽窄度三维采集的地震数据偏移成像结果(圆丘顶部成像的水平切片图),其中,图5(a)是窄三维采集的地震数据偏移,图5(b)是宽三维采集的地震数据偏移。
具体实施方式
中国石油天然气集团公司物探重点实验室从1996年开始探索一种针对勘探地区具体的复杂地质构造为目标的,以CRP地震成像为基础的,新的三维采集技术设计方法,该方法工作流程如图1所示。
首先根据勘探地区现有的勘探成果(复杂地质构造和复杂地质体形态及其地球物理参数)建立先验的地质地球物理模型,通常建立地下介质的地震速度模型;第二步以上述地质地球物理模型为勘探目标,设计两种或多种三维地震数据采集的观测系统;第三步以该模型作为三维地震观测对象,按照观测系统设计所要求的技术参数和施工图形在实验室模拟野外地震数据采集过程,本方法特点是采用数学计算和物理实验两种方法模拟野外地震施工,获得模拟的野外原始地震记录(参见图3),具体是采用虚谱法三维波动方程数值模拟与三维地震物理模型模拟两种方法联合正演来完成这一数据采集工作;第四步是对采集的数据进行三维地震数据处理,处理原理以共反射点(CRP)为基础,核心技术为叠前深度偏移处理,即一种当前流行的地震成像技术,从而得到地震处理几种成像结果;第五步以先验的地质地球物理模型作为评价与对比标准,以信噪比、分辩率和构造形态特征作为技术分析主要内容,对地震处理几种成像结果进行成像质量的综合对比评价,经过对比,确定其中一个最佳地震成像结果;如果这个成像结果还没有达到地震勘探的目标要求,则要修改原先的三维观测系统设计方案,随后重新进行三维地震数据采集处理工作;如果这个成像结果已经达到地震勘探的目标要求,则它所对应的三维观测系统就是该地区三维地震勘探优化设计的最后结果,随后输出三维地震观测系统优化设计工作流程图。此时,整个优化设计工作也同时宣告完成。
三维地震观测系统设计的主要内容是三维观测模板(几线、几炮制)选择和三维观测系统参数设计。三维观测系统参数主要有三维地震观测宽窄度、面元尺寸、覆盖或照明次数、照明度、地震道总数、最大炮检距等。在确定三维地震观测宽窄度N时,利用申请人提出的下列公式: 式中:Lc/LI-排列片横纵比;
nc/nI-覆盖次数横纵比;
C1-与排列片横纵比有关的系数;
C2-与覆盖次数横纵比有关的系数,且C2=1-C1,一般可
取C1=C2=0.5;
θ-子区中心面元方位角分布范围(弧度)。当:N<0.5时为窄方位观测系统;
N≥0.5时为宽方位观测系统;
N≥0.85时为全方位观测系统;
结合我国东西部复杂地区的油气勘探任务,中国石油天然气集团公司(CNPC)物探重点实验室多次将这一方法应用到实际的三维地震技术设计中,例如大港千米桥潜山模型,新疆车排子地区逆掩断裂带模型,新疆伊奇克里克褶皱冲断带模型,江苏句容和洪泽的复杂断块模型等,均取得了很好的效果。下面以大港千米桥潜山构造三维地震观测系统优化设计为例,对这一方法作一简单介绍。
图2(a)为根据先验资料制作的千米桥潜山构造物理模型,图2(b)为南北方向穿过潜山顶部的二维337测线地质剖面图。图3是在337测线上分别采用虚谱法三维波动方程正演模拟和采用三维地震物理模拟所得到的地震记录,比较这两张记录,可以认为两者是相似的。对该地区采用了几种不同的观测系统进行了三维地震数据采集工作,通过处理,得到了不同的地震成像结果(见图4和图5)。随后以地质目标作为标准,对各种成像结果进行评价,优化观测系统设计。例如图4(c)的图像信噪比比图4(b)和图4(a)的高,则应该认为图4(c)相对应的观测系统是其中最佳的;又例如图5(b)圆丘顶部成像比较好,聚焦质量比图5(b)高,说明应该选择宽方位三维观测系统进行地震施工。上述的一些实际应用表明:本方法适用范围广泛,应用效果良好,对指导石油天然气地震勘探实际生产、提高油气勘探成功率具有现实意义。
Claims (2)
1、一种基于地质地球物理模型的三维地震观测系统优化设计方法,包括以下步骤:
A)根据勘探地区现有的勘探成果,即复杂地质构造和复杂地质体形态及其地球物理参数建立先验的地质地球物理模型;
B)以上述地质地球物理模型为勘探目标,设计两种或多种三维地震数据采集的观测系统;
C)以该模型作为三维地震观测对象,按照观测系统设计所要求的技术参数和施工图形在实验室模拟野外地震数据采集过程,获得模拟的野外原始地震记录,具体作法是采用虚谱法三维波动方程数值模拟与三维地震物理模型模拟两种方法联合正演来完成这一数据采集工作;
D)对上述采集的数据进行三维地震数据处理,其中以共反射点为基础,对采集的数据进行叠前深度偏移处理,得到地震处理几种成像结果;
E)以先验的地质地球物理模型作为评价与对比标准,以信噪比、分辨率和构造形态特征作为技术分析主要内容,对所述地震处理几种成像结果进行成像质量的综合对比评价,经过对比,确定其中一个最佳地震成像结果。
2、根据权利要求1所述的方法,其中所述步骤E)还包括:
如果成像结果还没有达到地震勘探的目标要求,则修改原先的三维观测系统设计方案,随后重新进行三维地震数据采集处理工作;
如果成像结果已经达到地震勘探的目标要求,则输出三维地震观测系统优化设计流程图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB011204451A CN1158539C (zh) | 2001-07-13 | 2001-07-13 | 基于地质地球物理模型的三维地震观测系统优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB011204451A CN1158539C (zh) | 2001-07-13 | 2001-07-13 | 基于地质地球物理模型的三维地震观测系统优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1337583A true CN1337583A (zh) | 2002-02-27 |
CN1158539C CN1158539C (zh) | 2004-07-21 |
Family
ID=4664139
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB011204451A Expired - Fee Related CN1158539C (zh) | 2001-07-13 | 2001-07-13 | 基于地质地球物理模型的三维地震观测系统优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1158539C (zh) |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1299127C (zh) * | 2004-03-26 | 2007-02-07 | 中国石油天然气集团公司 | 地震观测系统优化设计的层状介质双聚焦方法及其应用 |
CN100344992C (zh) * | 2003-10-23 | 2007-10-24 | Pgs美洲公司 | 用于对叠前地震数据成像的方法 |
CN100429527C (zh) * | 2005-03-25 | 2008-10-29 | 中国石油天然气集团公司 | 三维地震资料处理质量监控技术 |
CN102053271A (zh) * | 2010-11-17 | 2011-05-11 | 中国石油天然气股份有限公司 | 地震采集质量分析方法及装置 |
CN102169187A (zh) * | 2010-12-01 | 2011-08-31 | 中国科学院地质与地球物理研究所 | 三维地震观测系统整体均衡性定量分析方法 |
CN101625417B (zh) * | 2008-07-08 | 2011-09-07 | 中国石油集团东方地球物理勘探有限责任公司 | 一种优化垂直地震剖面观测系统设计的方法 |
CN102495943A (zh) * | 2011-10-31 | 2012-06-13 | 国家卫星海洋应用中心 | 地球物理模型建模方法 |
CN102023310B (zh) * | 2009-09-19 | 2013-05-15 | 中国石油化工集团公司 | 三维观测系统均匀性分析方法 |
CN103116183A (zh) * | 2011-11-16 | 2013-05-22 | 中国石油天然气集团公司 | 一种石油地震采集面元覆盖次数属性体切片成图方法 |
CN103149590A (zh) * | 2013-02-26 | 2013-06-12 | 佟小龙 | 地球物理成像方法及装置 |
CN103777229A (zh) * | 2012-10-23 | 2014-05-07 | 中国石油化工股份有限公司 | 一种面向目的层的vsp观测系统设计方法 |
CN103901465A (zh) * | 2014-03-21 | 2014-07-02 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 全息三维地震勘探观测系统的设计方法 |
CN104050359A (zh) * | 2014-05-30 | 2014-09-17 | 中国石油大学(华东) | 一种基于三维观测系统排列片数据分割的正演模拟方法 |
CN104090297A (zh) * | 2014-06-14 | 2014-10-08 | 吉林大学 | 一种优化地震采集观测系统的逆向照明方法 |
CN104142518A (zh) * | 2013-10-29 | 2014-11-12 | 中国石油化工股份有限公司 | 一种地震观测系统叠前时间偏移响应的分析方法 |
CN104536041A (zh) * | 2014-12-17 | 2015-04-22 | 中国石油天然气集团公司 | 一种地震观测系统参数的优化方法 |
CN104635259A (zh) * | 2013-11-06 | 2015-05-20 | 中国石油天然气集团公司 | 一种确定观测系统排列片范围的方法及装置 |
CN109031462A (zh) * | 2017-06-08 | 2018-12-18 | 道达尔公司 | 用于评估感兴趣的区域上的地球物理勘测采集几何结构的方法 |
CN111723517A (zh) * | 2020-05-26 | 2020-09-29 | 河海大学 | 一种融合机器学习原理和逐渐变形法的地球物理模型优化方法 |
CN114722590A (zh) * | 2022-03-23 | 2022-07-08 | 中石化石油工程技术服务有限公司 | 基于地球物理模型的随机采集观测系统设计优选方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102147480A (zh) * | 2011-01-12 | 2011-08-10 | 中国石油集团川庆钻探工程有限公司 | 一种三维观测系统测试方法 |
-
2001
- 2001-07-13 CN CNB011204451A patent/CN1158539C/zh not_active Expired - Fee Related
Cited By (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100344992C (zh) * | 2003-10-23 | 2007-10-24 | Pgs美洲公司 | 用于对叠前地震数据成像的方法 |
CN1299127C (zh) * | 2004-03-26 | 2007-02-07 | 中国石油天然气集团公司 | 地震观测系统优化设计的层状介质双聚焦方法及其应用 |
CN100429527C (zh) * | 2005-03-25 | 2008-10-29 | 中国石油天然气集团公司 | 三维地震资料处理质量监控技术 |
CN101625417B (zh) * | 2008-07-08 | 2011-09-07 | 中国石油集团东方地球物理勘探有限责任公司 | 一种优化垂直地震剖面观测系统设计的方法 |
CN102023310B (zh) * | 2009-09-19 | 2013-05-15 | 中国石油化工集团公司 | 三维观测系统均匀性分析方法 |
CN102053271A (zh) * | 2010-11-17 | 2011-05-11 | 中国石油天然气股份有限公司 | 地震采集质量分析方法及装置 |
CN102053271B (zh) * | 2010-11-17 | 2012-09-26 | 中国石油天然气股份有限公司 | 地震采集质量分析方法及装置 |
CN102169187B (zh) * | 2010-12-01 | 2013-01-23 | 中国科学院地质与地球物理研究所 | 三维地震观测系统整体均衡性定量分析方法 |
CN102169187A (zh) * | 2010-12-01 | 2011-08-31 | 中国科学院地质与地球物理研究所 | 三维地震观测系统整体均衡性定量分析方法 |
CN102495943A (zh) * | 2011-10-31 | 2012-06-13 | 国家卫星海洋应用中心 | 地球物理模型建模方法 |
CN102495943B (zh) * | 2011-10-31 | 2015-04-22 | 国家卫星海洋应用中心 | 地球物理模型建模方法 |
CN103116183A (zh) * | 2011-11-16 | 2013-05-22 | 中国石油天然气集团公司 | 一种石油地震采集面元覆盖次数属性体切片成图方法 |
CN103116183B (zh) * | 2011-11-16 | 2015-05-27 | 中国石油天然气集团公司 | 一种石油地震采集面元覆盖次数属性体切片成图方法 |
CN103777229A (zh) * | 2012-10-23 | 2014-05-07 | 中国石油化工股份有限公司 | 一种面向目的层的vsp观测系统设计方法 |
CN103149590B (zh) * | 2013-02-26 | 2016-01-27 | 佟小龙 | 地球物理成像方法及装置 |
CN103149590A (zh) * | 2013-02-26 | 2013-06-12 | 佟小龙 | 地球物理成像方法及装置 |
CN104142518A (zh) * | 2013-10-29 | 2014-11-12 | 中国石油化工股份有限公司 | 一种地震观测系统叠前时间偏移响应的分析方法 |
CN104142518B (zh) * | 2013-10-29 | 2017-05-03 | 中国石油化工股份有限公司 | 一种地震观测系统叠前时间偏移响应的分析方法 |
CN104635259B (zh) * | 2013-11-06 | 2017-09-01 | 中国石油天然气集团公司 | 一种确定观测系统排列片范围的方法及装置 |
CN104635259A (zh) * | 2013-11-06 | 2015-05-20 | 中国石油天然气集团公司 | 一种确定观测系统排列片范围的方法及装置 |
CN103901465A (zh) * | 2014-03-21 | 2014-07-02 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 全息三维地震勘探观测系统的设计方法 |
CN104050359A (zh) * | 2014-05-30 | 2014-09-17 | 中国石油大学(华东) | 一种基于三维观测系统排列片数据分割的正演模拟方法 |
CN104050359B (zh) * | 2014-05-30 | 2016-10-05 | 中国石油大学(华东) | 一种基于三维观测系统排列片数据分割的正演模拟方法 |
CN104090297A (zh) * | 2014-06-14 | 2014-10-08 | 吉林大学 | 一种优化地震采集观测系统的逆向照明方法 |
CN104536041A (zh) * | 2014-12-17 | 2015-04-22 | 中国石油天然气集团公司 | 一种地震观测系统参数的优化方法 |
CN109031462A (zh) * | 2017-06-08 | 2018-12-18 | 道达尔公司 | 用于评估感兴趣的区域上的地球物理勘测采集几何结构的方法 |
CN111723517A (zh) * | 2020-05-26 | 2020-09-29 | 河海大学 | 一种融合机器学习原理和逐渐变形法的地球物理模型优化方法 |
CN111723517B (zh) * | 2020-05-26 | 2022-09-02 | 河海大学 | 一种融合机器学习原理和逐渐变形法的地球物理模型优化方法 |
CN114722590A (zh) * | 2022-03-23 | 2022-07-08 | 中石化石油工程技术服务有限公司 | 基于地球物理模型的随机采集观测系统设计优选方法 |
CN114722590B (zh) * | 2022-03-23 | 2024-03-22 | 中石化石油工程技术服务有限公司 | 基于地球物理模型的随机采集观测系统设计优选方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1158539C (zh) | 2004-07-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1158539C (zh) | 基于地质地球物理模型的三维地震观测系统优化设计方法 | |
CN108957549B (zh) | 一种辫状河沉积非均质致密砂岩气藏地质建模方法 | |
CN102053270B (zh) | 一种基于沉积地层单元的地震相分析方法 | |
US8694261B1 (en) | 3D-well log invention | |
CN102759745B (zh) | 一种基于数字地质露头模型正演的碳酸盐岩储层预测方法 | |
Wilson et al. | From outcrop to flow simulation: Constructing discrete fracture models from a LIDAR survey | |
CN104570125B (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
CN100557464C (zh) | 基于叠前波场模拟的地震勘探层位标定方法 | |
CN101329407B (zh) | 一种快速p-sv转换波直接模拟以确定地层岩性、岩相变化的方法 | |
CN111596978A (zh) | 用人工智能进行岩相分类的网页显示方法、模块和系统 | |
CN104142518B (zh) | 一种地震观测系统叠前时间偏移响应的分析方法 | |
CN104678434A (zh) | 一种预测储层裂缝发育参数的方法 | |
CN102967883A (zh) | 通过页岩气叠前弹性参数反演预测岩石脆性概率的方法 | |
CN111722284B (zh) | 一种基于道集数据建立速度深度模型的方法 | |
CN106777707B (zh) | 一种利用改进的蜘蛛网图进行测井岩性定量识别的方法 | |
US6249746B1 (en) | Automated seismic isochron analysis | |
CN112505754B (zh) | 基于高精度层序格架模型的井震协同划分沉积微相的方法 | |
Tebo et al. | Use of volume-based 3-D seismic attribute analysis to characterize physical-property distribution: a case study to delineate sedimentologic heterogeneity at the Appleton field, southwestern Alabama, USA | |
Eppinger et al. | 2d near‐surface full‐waveform tomography reveals bedrock controls on critical zone architecture | |
CN109143398A (zh) | 一种自动网格层析深度域速度的建模方法 | |
CN102147480A (zh) | 一种三维观测系统测试方法 | |
Nahm et al. | Interpretation and practical applications of 4C-3D seismic data, East Cameron gas fields, Gulf of Mexico | |
Bryant et al. | Integrated 3-D geological modeling of the C1 Sands reservoir, Maui Field, offshore New Zealand | |
Wang et al. | Lithological Reservoir Exploration Technology | |
Adlakha | Fracture analysis and modelling of the South Arne field |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |