CN104237943A - 无表层约束层析反演静校正方法 - Google Patents
无表层约束层析反演静校正方法 Download PDFInfo
- Publication number
- CN104237943A CN104237943A CN201310240759.2A CN201310240759A CN104237943A CN 104237943 A CN104237943 A CN 104237943A CN 201310240759 A CN201310240759 A CN 201310240759A CN 104237943 A CN104237943 A CN 104237943A
- Authority
- CN
- China
- Prior art keywords
- cmp
- pile
- surface layer
- section
- static correction
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明是提出了完全不依赖表层调查资料的无表层约束层析反演静校正方法,直接对初至数据在CMP域分段,建立CMP段作为控制点,对每一段内所有CMP点初至按偏移距排列形成时距图;然后对时距图采用折射法进行解释,得到该CMP段的解释结果;将所有的CMP段解释结果进行内插就得到了整个测线表层初始模型;最后以此初始模型约束通过层析反演表层最终模型,计算静校正量。本发明在复杂山地区提高地表的约束控制密度,同时利用远偏移距折射初至信息提高对地表纵向的解释深度,克服了复杂山地微测井密度和深度难以满足建模要求的难点,解决复杂山地表层建模与静校正问题。
Description
技术领域
本发明涉及地球物理勘探方法,是一种不依赖表层调查资料约束的无表层约束层析反演静校正方法。
背景技术
在复杂山地区,地表起伏剧烈,表层风化严重,表层速度横向变化快,低降速带巨厚,并且折射界面不稳定或者难以识别的地区勘探,对地震数据常规的静校正方法很难解决复杂山地区的静校正问题。目前公开的基于表层调查内插的模型静校正和初至折射静校正方法视图解决这些问题,它们充分利用了不同偏移距各类波(直达波、折射波,透射波、反射波等)的初至,不需要稳定折射层,对速度、厚度变化剧烈能较好适应,是解决复杂山地静校正问题的手段。
上述的从初至层析反演静校正看,要减少反演多解性,提高复杂山地层析反演精度,必须引入先验信息进行约束反演,这样,建立比较合理的初始近地表模型是层析反演静校正的基础。对于复杂山地,目前主要依赖于微测井调查建立初始近地表模型。但是利用微测井建模存在以下问题:1、复杂山地地形起伏和表层的剧烈变化,需要有足够密度的微测井表层调查控制点来精确描述表层速度结构,必将消耗巨大的人力、物力资源成本,无形中增加了勘探成本和HSE(健康、安全、环保标准)风险,也会制约野外勘探生产效率;2、受山地钻机能力的制约,在低降速带巨厚区,微测井无法打穿低降速层,调查深度不够,难以满足近地表建模精度要求。基于上述理由,实际中依靠表层调查资料进行约束层析反演会遇到困难,结果也不够精确。
发明内容
本发明针对现有技术存在的问题,提出了一种降低勘探成本、提高效率,不依赖野外表层调查资料,直接用CMP域初至构建近地表的初始模型的无表层约束层析反演静校正方法。
本发明通过以下步骤实现:
1)将探区整条测线按CMP(共中心点)桩号进行分段;
所述的分段是:将测线按CMP桩号均匀分开;
所述的分段是:对测线CMP桩号为1-100,分为10段,桩号1-10为一段,桩号11-20为一段,其余类推;
所述的分段是:对采用三维勘探数据,在inline沿测线方向将测线按CMP桩号均匀分开,crossline垂直测线方向按CMP线号均匀分开,形成CMP块;
所述的分段是:对采用三维勘探数据,inline方向桩号1-10为一段,与crossline垂直测线方向CMP线号101-105为一段围成一个CMP块;
2)对落在每一段或块内的所有CMP的初至按绝对偏移距排序,形成CMP段或块的时距图;
所述的排序是:按每个炮检对的检波点桩号或炮点桩号减去CMP桩号的绝对值排序。
所述的时距图横轴为偏移距,纵轴为该炮检对的初至时间。
3)通过线性拟合时距曲线的斜率求出低降速带不同地层的速度,再利用折射时距曲线与时间轴的交点求出交叉时t,进而通过折射波时距方程求出各层厚度,从而求出该CMP段对应的近地表速度、厚度信息;
4)将所有CMP段的速度、厚度信息内插出整条测线的表层速度、厚度信息,得到表层模型;
5)以上述表层模型为约束条件,进行初至层析反演近地表最终模型,计算得到静校正量。
本发明采用CMP域初至,通过CMP灵活分段,在复杂山地区通过划分较密的CMP段(块)提高对复杂地表的约束控制密度,同时利用远偏移距折射初至信息提高对地表纵向的解释深度,克服了复杂山地微测井密度和深度难以满足建模要求的难点,充分发挥其在密度和深度上的优势,可构建更合理的初始地表模型,结合层析反演,建立近地表更深的速度场,解决复杂山地表层建模与静校正问题。
附图说明
图1CMP域初至分段示意图。
图2CMP段内初至时距示意图。
图3无表层约束层析反演静校正方法流程图。
具体实施方式
以下结合附图详细说明本发明。
步骤一:将探区整条测线按CMP(共中心点)桩号进行分段。
所述的分段是:将测线按CMP桩号均匀分开;
所述的分段是:对测线CMP桩号为1-100,分为10段,桩号1-10为一段,桩号11-20为一段,其余类推;如图1所示,CMP桩号1-3为一CMP段。
所述的分段是:对采用三维勘探数据,在inline沿测线方向将测线按CMP桩号均匀分开,crossline垂直测线方向按CMP线号均匀分开,形成CMP块;
所述的分段是:对采用三维勘探数据,inline方向桩号1-10为一段,与crossline垂直测线方向CMP线号101-105为一段围成一个CMP块;
步骤二:对落在每一段或块内的所有CMP的初至按绝对偏移距排序,形成CMP段或块的时距图;如示意图2,该CMP段含3个CMP点的初至。
所述的排序是:按每个炮检对的检波点桩号或炮点桩号减去CMP桩号的绝对值排序。
所述的时距图横轴为偏移距,纵轴为该炮检对的初至时间。
步骤三:CMP段(块)速度分层。利用时距曲线与时间轴的交点求出交叉时t,通过公式(1)、(2)求出各层厚度,从而建立每个CMP段近地表初始模型。
第一层厚度: 公式(1)
第一层厚度: 公式(2)
其中:Tj为第j层的交时,为第j层的速度,均可通过线拟合求出;为地层厚度,可通过公式(1)、(2)求出。
步骤四:将所有CMP段的速度、厚度信息内插出整条测线的表层速度、厚度信息,得到表层模型;
步骤五:以上述表层模型为约束条件,进行初至层析反演近地表最终模型,计算得到静校正量。图3为该方法流程图。
Claims (7)
1.一种无表层约束层析反演静校正方法,特点是通过以下步骤实现:
1)将探区整条测线按CMP桩号进行分段;
2)对落在每一段或块内的所有CMP的初至按绝对偏移距排序,形成CMP段或块的时距图;
3)通过线性拟合时距曲线的斜率求出低降速带不同地层的速度,再利用折射时距曲线与时间轴的交点求出交叉时t,进而通过折射波时距方程求出各层厚度,从而求出该CMP段对应的近地表速度、厚度信息;
4)将所有CMP段的速度、厚度信息内插出整条测线的表层速度、厚度信息,得到表层模型;
5)以上述表层模型为约束条件,进行初至层析反演近地表最终模型,计算得到静校正量。
2.根据权利要求1所述的方法,特点是步骤1)所述的分段是:将测线按CMP桩号均匀分开。
3.根据权利要求1所述的方法,特点是步骤1)所述的分段是:对测线CMP桩号为1-100,分为10段,桩号1-10为一段,桩号11-20为一段,其余类推。
4.根据权利要求1所述的方法,特点是步骤1)所述的分段是:对采用三维勘探数据,在inline沿测线方向将测线按CMP桩号均匀分开,crossline垂直测线方向按CMP线号均匀分开,形成CMP块。
5.根据权利要求1所述的方法,特点是步骤1)所述的分段是:对采用三维勘探数据,inline方向桩号1-10为一段,与crossline垂直测线方向 CMP线号101-105为一段围成一个CMP块。
6.根据权利要求1所述的方法,特点是步骤2)所述的排序是:按每个炮检对的检波点桩号或炮点桩号减去CMP桩号的绝对值排序。
7.根据权利要求1所述的方法,特点是步骤2)所述的时距图横轴为偏移距,纵轴为该炮检对的初至时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310240759.2A CN104237943B (zh) | 2013-06-18 | 2013-06-18 | 无表层约束层析反演静校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310240759.2A CN104237943B (zh) | 2013-06-18 | 2013-06-18 | 无表层约束层析反演静校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104237943A true CN104237943A (zh) | 2014-12-24 |
CN104237943B CN104237943B (zh) | 2017-09-08 |
Family
ID=52226387
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310240759.2A Active CN104237943B (zh) | 2013-06-18 | 2013-06-18 | 无表层约束层析反演静校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104237943B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104898165A (zh) * | 2015-03-20 | 2015-09-09 | 中国石油天然气股份有限公司 | 一种黄土塬近道约束层析静校正方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1480134A2 (en) * | 2003-05-22 | 2004-11-24 | Cognos Incorporated | System and method for logging actions for transformating a data model |
CN1948999A (zh) * | 2005-10-10 | 2007-04-18 | 中国石油天然气股份有限公司 | 近似层替换静校正方法 |
CN101598803A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种直接得到转换波叠加剖面的方法 |
US20100256916A1 (en) * | 2009-04-03 | 2010-10-07 | Chevron U.S.A. Inc. | Method for target-oriented reverse time migration for prestack depth imaging |
CN102073064A (zh) * | 2009-11-25 | 2011-05-25 | 中国石油天然气集团公司 | 一种利用相位信息提高速度谱分辨率的方法 |
-
2013
- 2013-06-18 CN CN201310240759.2A patent/CN104237943B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1480134A2 (en) * | 2003-05-22 | 2004-11-24 | Cognos Incorporated | System and method for logging actions for transformating a data model |
CN1948999A (zh) * | 2005-10-10 | 2007-04-18 | 中国石油天然气股份有限公司 | 近似层替换静校正方法 |
CN101598803A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种直接得到转换波叠加剖面的方法 |
US20100256916A1 (en) * | 2009-04-03 | 2010-10-07 | Chevron U.S.A. Inc. | Method for target-oriented reverse time migration for prestack depth imaging |
CN102073064A (zh) * | 2009-11-25 | 2011-05-25 | 中国石油天然气集团公司 | 一种利用相位信息提高速度谱分辨率的方法 |
Non-Patent Citations (4)
Title |
---|
熊定钰: "复杂地表静校正方法研究及应用", 《中国博士学位论文全文数据库 基础科学辑》 * |
王咸彬: "西部复杂山地地震勘探关键技术应用研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
马云东: "叠加法静校正的实现与研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
齐春艳,等: "大庆长垣油田近地表模型约束折射波层析静校正", 《石油地球物理勘探》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104898165A (zh) * | 2015-03-20 | 2015-09-09 | 中国石油天然气股份有限公司 | 一种黄土塬近道约束层析静校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104237943B (zh) | 2017-09-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106054248B (zh) | 一种基于大面积致密储层地震岩石物理反演方法 | |
CN104267429B (zh) | 确定地层压力的方法及装置 | |
CN103258091B (zh) | 非常规油气藏水平井段三维岩体力学模型建立的方法及装置 | |
CN103573250B (zh) | 一种计算水平井井眼到地层上下界面距离的方法 | |
CN103410504B (zh) | 一种水平井/大斜度井确定地层真电阻率的方法及装置 | |
CN104297787B (zh) | 河流相低渗致密砂岩储层的三维岩相数据处理方法及装置 | |
CN105510993A (zh) | 前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法 | |
CN102636812B (zh) | 一种获得碳酸盐岩储层储集空间体积的方法 | |
CN104635261B (zh) | 山前带重磁电震联合分体建模方法 | |
CN105549082B (zh) | 超深碳酸盐储层三维地质力学场建立方法和系统 | |
CN107329171A (zh) | 深度域储层地震反演方法及装置 | |
CN104849762B (zh) | 利用水平井入层点信息校正油藏顶面微构造的方法 | |
CN104375172B (zh) | 火山岩下伏地层的构造形态正确成像的方法 | |
CN103293562A (zh) | 一种碳酸盐岩储层地质储量的确定方法及设备 | |
CN108121009B (zh) | 开发后期复杂断块变速构造成图方法 | |
CN103149588B (zh) | 一种利用井震标定计算vti各向异性参数的方法及系统 | |
CN103838936A (zh) | 一种适用于浊积砂低渗透储层的高精度构造应力场模拟方法 | |
CN107728205B (zh) | 一种地层压力预测方法 | |
CN103064118A (zh) | 一种利用声波测井定量洞穴充填程度的方法 | |
CN109164487A (zh) | 一种基于模型建立平均速度场的方法及精细成图的方法 | |
CN105738952A (zh) | 一种水平井区储层岩石相建模方法 | |
CN103454676B (zh) | 一种非线性等时地层切片的确定方法及系统 | |
CN105626056A (zh) | 一种确定目标区三维地层孔隙压力的方法和装置 | |
CN103399345A (zh) | 一种潜山裂缝分布的勘测方法与装置 | |
CN104422955A (zh) | 一种利用旅行时变化量进行各向异性参数提取的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |