CN112526610B - 一种约束表层建模的三维地震采集激发井深设计方法 - Google Patents
一种约束表层建模的三维地震采集激发井深设计方法 Download PDFInfo
- Publication number
- CN112526610B CN112526610B CN201910877008.9A CN201910877008A CN112526610B CN 112526610 B CN112526610 B CN 112526610B CN 201910877008 A CN201910877008 A CN 201910877008A CN 112526610 B CN112526610 B CN 112526610B
- Authority
- CN
- China
- Prior art keywords
- model
- well depth
- velocity
- data
- excitation
- 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.)
- Active
Links
- 230000005284 excitation Effects 0.000 title claims abstract description 61
- 238000013461 design Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000002344 surface layer Substances 0.000 title claims abstract description 27
- 239000010410 layer Substances 0.000 claims abstract description 28
- 238000005259 measurement Methods 0.000 claims abstract description 15
- 238000011160 research Methods 0.000 claims description 25
- 238000009499 grossing Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 abstract description 8
- 238000012937 correction Methods 0.000 description 8
- 230000003068 static effect Effects 0.000 description 7
- 238000011835 investigation Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 5
- 238000010276 construction Methods 0.000 description 3
- 230000009189 diving Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 238000004883 computer application Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种约束表层建模的三维地震采集激发井深设计方法,属于地震资料采集处理技术领域,通过更新初始速度模型能够得到用于设计激发井深的高精度近地表速度模型,即先建立初始速度模型,根据该初始速度模型反演建立近地表模型,结合地震采集的新测量成果、表层调查点资料的约束条件,在近地表模型的基础上反演得到高精度近地表速度模型。最后将采集区的设计炮点数据和不同地表类型区域对应的激发层速度界面代入该高精度近地表速度模型,从而得到高精度的激发井深。本发明的三维地震采集激发井深设计方法提高了速度模型的精度,优化了井深设计,且优化后的激发井深精度较高,能够保证地震资料品质、低幅构造和属性预测精度。
Description
技术领域
本发明属于地震资料采集处理技术领域,具体涉及一种约束表层建模的三维地震采集激发井深设计方法。
背景技术
地震采集是油气地震勘探中的第一道工序,地震采集是指按照地震采集设计,在地面激发地震波,通过地层传播后的地震信号经过放大、滤波和数字转换被地面检波器接收的过程。在井炮地震勘探野外采集中,需通过合理精确的分线设计保证在有利激发层中进行激发,因此,在地震勘探施工中首先需要通过激发试验和近地表调查获得采集区的最佳激发岩性及其空间展布,以合理设计激发参数。
常规激发井深设计方法是:首先进行小折射、微测井和潜水面等表层结构调查,确定低降速带厚度,然后通过经验和各类内插算法计算得到全区的最佳激发岩性界面分布,结合激发试验论证,进而确定最佳激发井深,具体参见作者葛利华等人于2017年3月在期刊《吉林大学学报》第47卷第2期公开的文献《辽西葫芦岛东部表层调查方法比对实验》、作者吕芙瑶等人于2013年6月在期刊《石油工业计算机应用》第2期公开的文献《基于局部矩形网格构建的激发井深设计》,和魏明阳在于2014年在期刊《地球物理学进展》第29卷第6期公开的文献《巴楚地区基于多重条件下约束下激发参数设计》。常规激发井深设计方法在表层结构变化小、高速层速度横向差异不大的地区,能取得较好的效果,但在地表高程起伏大、表层结构横向变化较大的区域,低测资料难以控制空间速度和厚度的变化,导致设计井深出现内插精度不足和随地形起伏变化的问题,降低了根据上述设计井深得到的地震资料品质,甚至会出现把近地表异常变化引起的振幅不一致性当成地质现象的错误解释。
发明内容
本发明的目的是提供一种约束表层建模的三维地震采集激发井深设计方法,用于解决采用常规激发井深设计方法设计井深出现内插精度不足和随地形起伏变化的问题。
基于上述目的,三维地震采集激发井深设计方法的技术方案如下:
获取研究区的前期地震数据和前期表层资料,根据所述前期地震数据和确定的近地表模型反演深度和反演偏移距优化大炮初至;根据优化后的大炮初至,结合所述近地表模型反演深度,建立初始速度模型;将所述前期表层资料作为约束条件,在所述初始速度模型的基础上进行反演,建立近地表模型;
获取研究区地震采集的新测量成果和新表层调查点资料,该新测量成果包括设计炮点数据,新表层调查点资料为微测井得到的数据;
将所述新测量成果和新表层调查点资料作为新的约束条件,在近地表模型的基础上进行反演,得到高精度近地表速度模型;
根据获取采集区的设计炮点数据和不同地表类型区域对应的激发层速度界面,结合所述高精度近地表速度模型,得到相应设计炮点的激发井深。
上述技术方案的有益效果为:
本发明的一种约束表层建模的三维地震采集激发井深设计方法,通过更新初始速度模型能够得到用于得到激发井深的高精度近地表速度模型,即结合研究区的地表和低降速带变化特点,依据大炮初至先建立初始速度模型,根据该初始速度模型反演建立近地表模型,再将结合新测量成果和新表层调查点资料作为约束条件,在近地表模型的基础上反演得到高精度近地表速度模型。最后将采集区的设计炮点数据和不同地表类型区域对应的激发层速度界面代入该高精度近地表速度模型,从而得到高精度的激发井深。本发明的三维地震采集激发井深设计方法提高了速度模型的精度,优化了井深设计,且优化后的激发井深精度较高,能够保证地震资料品质、低幅构造和属性预测精度,不存在内插精度不足和随地形起伏变化的问题。
进一步的,所述近地表模型反演深度和反演偏移距是根据研究区的地表和低降速带变化特点确定的。
进一步的,所述激发层速度界面是依据研究区的地表类型、试验论证参数和表层岩性调查确定的。
进一步的,在将所述激发层速度界面作为输入加载至所述高精度近地表模型前,还包括对激发层速度界面进行大半径平滑,得到统一平滑半径的激发层速度界面。
附图说明
图1是本发明的三维地震采集激发井深设计方法流程图;
图2是本发明的表层资料转化为适合静校正处理软件的文本图;
图3是本发明的设计炮点激发井深文本图;
图4是现有技术中采用常规激发井深方法得到的地震单炮记录图;
图5是采用本发明方法得到的地震单炮记录图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的说明。
本实施例提出一种约束表层建模的三维地震采集激发井深设计方法,如图1所示,包括以下步骤:
收集研究区的前期地震数据和前期表层资料(小折射和微测井得到的数据);依据研究区的地表和低降速带变化特点,确定近地表模型反演深度和反演偏移距,优化大炮初至;依据优化后的大炮初至,结合近地表模型反演深度,建立近地表初始速度模型;将所述前期表层资料作为约束条件,在近地表初始速度模型的基础上进行反演建模,建立近地表模型。
收集研究区地震采集的新测量成果(包括设计炮点的X、Y坐标、高程)和新表层调查点资料(微测井得到的数据),将新测量成果和表层调查点资料作为新的约束条件,在近地表模型的基础上进行反演,得到高精度近地表速度模型。
依据研究区地震采集施工范围和技术设计,获取采集区的设计炮点数据;依据研究区的地表类型、试验论证参数和表层岩性调查,确定不同地表类型区域对应的激发层速度界面;将设计炮点数据和激发层速度界面作为输入,加载到高精度近地表速度模型中,输出地震采集区设计炮点的激发井深。
以某盆地某工区的三维地震资料采集处理为例,已知该研究区范围为685km2,属于沙漠草原地貌,地表高程在1340-1500米之间,地形整体趋势相对平坦(分为高地区、低洼区、断裂区和斜坡区),局部起伏高差达近百米,沙土覆盖厚,近地表厚度和速度横向变化剧烈;研究区的地质背景相对简单,地震成果显示中生界和古生界地层极为平缓,整体呈东北高西南低的斜坡,东北部位发育大型应力构造,断裂出露地表;研究区内已有前期两个年度的二维地震采集资料,共计50条测线,密度为4×2km,表层低测资料点(小折射、微测井)近千个,有已完钻井约8口(包括探井、直井、水平井),井点分布相对均匀。研究区的地表类型多样、低降速带的厚度和速度横向变化剧烈增加了地震采集井深设计的难度,前期地震资料(即前期地震数据)和低测资料丰富,具备多信息约束建模的基本条件,且常规激发井深设计方法存在内插精度不足和随地形起伏变化的问题,因此,采用本发明的三维地震采集激发井深设计方法,设计研究区地震采集的激发井深,具体操作步骤如下:
(1)收集研究区的前期地震数据和前期表层资料;前期地震数据为地震采集技术领域标准segy或segd格式的炮集数据和炮集数据对应的地震采集技术领域标准sps文件;将所述前期地震数据在静校正处理软件中打开,通过线性动校正后逐点拾取初至旅行时,得到dat或txt格式的大炮初至,作为建立近地表速度模型和本发明进行反演的资料基础;前期表层资料为研究区前期实施的各个低测调查点的小折射和微测井解释成果,包括表层调查点的XY坐标、岩性特征和解释成果分层数据表格。
(2)依据研究区的地表和低降速带变化特点,包括地表类型、低降速带厚度和速度的纵向、横向变化,确定近地表模型反演深度和反演偏移距,优化大炮初至;针对反演参数,包括反演深度、反演偏移距、反演次数、反演面元等,通过迭代误差和射线密度的指标,测试近地表模型反演深度和反演偏移距范围,进而确定大炮初至;依据大炮初至,结合反演深度,在静校正处理相关软件中点取分层数据的速度和厚度,建立近地表初始速度模型。
(3)将前期表层资料作为约束条件,在初始速度模型的基础上进行层析反演,建立近地表模型。
该步骤隐含一个简化处理,由于收集到的所有表层资料是解释成果分层数据表格,需要将此表格的数据转化为适合静校正处理软件的文本文件,该文本文件是包含X坐标、Y坐标、低速层厚度H0、低速层速度V0、降速层厚度H1、降速层速度V1、高速层厚度H2、高速层速度V2等八列的prn或txt文件,如图2所示;将此文件加载到处理软件中,作为约束条件,采用约束层析反演算法,建立近地表模型。
(4)收集和加载研究区地震采集的新测量成果和新表层调查点资料,将此作为新的约束条件,在步骤(3)得到的近地表模型的基础上进行反演建模,获得高精度近地表速度模型。
(5)依据研究区地震采集施工范围和技术设计,获取采集区的设计炮点数据。
该步骤隐含一个简化处理,需要将采集区的设计炮点数据转化为适合静校正处理软件的文本文件,该文本文件是包含线点号、X坐标、Y坐标、地表高程、插值井深等五列的prn或txt文件。
(6)获取试验论证参数和表层岩性调查,其中,试验论证参数主要包括单炮记录的一些定量化指标,比如:能量、频率、信噪比、子波等;表层岩性调查主要是指微测井岩性解释,按照表层调查(表层调查包括地表露头调查、潜水面调查和微测井解释),确定微测井的深度,以及每一段深度对应的岩性,比如:表层沙、含泥砂岩、含砾砂岩及对应的深度和速度。
依据研究区的地表类型、试验论证参数和表层岩性调查,确定不同地表类型区域对应的激发层速度界面,以此适应不同地表类型区域的井深设计。比如:地表类型为高地平坦区,结合该区域的表层岩性调查中的速度,用单炮记录的试验论证参数,确定该区域的最佳激发层速度面。
由于激发层速度界面是按照速度固定值确定的,在速度模型上会突变,为消除速度突变,保证激发层速度界面的速度趋势,激发层速度界面需要进行大半径平滑,平滑效果要足以体现层位速度的变化趋势;不同地表类型的激发层速度界面采用统一的平滑半径;
(7)将步骤(5)得到的设计炮点数据和步骤(6)得到的激发层速度界面作为输入,加载到步骤(4)得到的高精度速度模型中,输出地震采集区设计炮点的激发井深;计算得到的激发井深以文本格式导出,如图3所示,各列依次为线点号、X坐标、Y坐标、地表高程、基底厚度、设计井深、原始插值井深、设计误差等八列的prn或txt文件。
本发明主要利用研究区内的高精度近地表速度模型的速度趋势来设计地震采集激发井深,并不是将前期表层资料通过内插算法进行插值,不存在内插精度不足和随地形起伏变化的情况。依据本发明的三维地震采集激发井深设计方法得到的激发井深,通过单炮质量对比分析看出,地震资料品质得到大幅提高,如图4和图5所示,两炮的位置相距20m,药量相同,图4井深为16m,图5井深为37.2m。本发明的三维地震采集激发井深设计方法的主要优势是利用了研究区的多种信息融合,获得的高精度近地表速度模型进一步削弱常规静校正对处理的影响,降低处理成果的解释误差,提高静校正处理精度和地震成果属性预测精度。
Claims (4)
1.一种约束表层建模的三维地震采集激发井深设计方法,其特征在于,包括以下步骤:
获取研究区的前期地震数据和前期表层资料,根据所述前期地震数据和确定的近地表模型反演深度和反演偏移距优化大炮初至;根据优化后的大炮初至,结合所述近地表模型反演深度,建立初始速度模型;将所述前期表层资料作为约束条件,在所述初始速度模型的基础上进行反演,建立近地表模型;
获取研究区地震采集的新测量成果和新表层调查点资料,该新测量成果包括设计炮点数据,新表层调查点资料为微测井得到的数据;
将所述新测量成果和新表层调查点资料作为新的约束条件,在近地表模型的基础上进行反演,得到高精度近地表速度模型;
根据获取采集区的设计炮点数据和不同地表类型区域对应的激发层速度界面,结合所述高精度近地表速度模型,得到相应设计炮点的激发井深。
2.根据权利要求1所述的约束表层建模的三维地震采集激发井深设计方法,其特征在于,所述近地表模型反演深度和反演偏移距是根据研究区的地表和低降速带变化特点确定的。
3.根据权利要求1所述的约束表层建模的三维地震采集激发井深设计方法,其特征在于,所述激发层速度界面是依据研究区的地表类型、试验论证参数和表层岩性调查确定的。
4.根据权利要求1所述的约束表层建模的三维地震采集激发井深设计方法,其特征在于,在将所述激发层速度界面作为输入加载至所述高精度近地表模型前,还包括对激发层速度界面进行大半径平滑,得到统一平滑半径的激发层速度界面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910877008.9A CN112526610B (zh) | 2019-09-17 | 2019-09-17 | 一种约束表层建模的三维地震采集激发井深设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910877008.9A CN112526610B (zh) | 2019-09-17 | 2019-09-17 | 一种约束表层建模的三维地震采集激发井深设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112526610A CN112526610A (zh) | 2021-03-19 |
CN112526610B true CN112526610B (zh) | 2023-03-21 |
Family
ID=74974673
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910877008.9A Active CN112526610B (zh) | 2019-09-17 | 2019-09-17 | 一种约束表层建模的三维地震采集激发井深设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112526610B (zh) |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101788686B (zh) * | 2010-01-06 | 2012-03-14 | 中国石油化工集团公司 | 基于一致近地表模型的逐点激发井深设计方法 |
US9091786B2 (en) * | 2010-10-15 | 2015-07-28 | Oz YILMAZ | Image based effective medium modeling of the near surface earth formation |
CN103869368B (zh) * | 2012-12-10 | 2016-05-11 | 中国石油天然气集团公司 | 一种无表层调查资料约束的大炮初至综合建模静校正方法 |
WO2015196779A1 (zh) * | 2014-06-24 | 2015-12-30 | 夏正元 | 利用叠加速度求取高精度地震波速度的方法 |
CN104870067A (zh) * | 2014-09-19 | 2015-08-26 | 杨顺伟 | 一种确定均匀介质中最佳激发井深的方法 |
CN105277978B (zh) * | 2015-10-15 | 2018-01-05 | 中国石油天然气集团公司 | 一种确定近地表速度模型的方法及装置 |
CN106569282B (zh) * | 2016-10-17 | 2019-03-12 | 中国石油化工股份有限公司 | 一种地震采集激发井深设计方法 |
CN109143366A (zh) * | 2017-06-27 | 2019-01-04 | 中国石油化工股份有限公司 | 近地表初至层析静校正方法及计算机可读存储介质 |
CN108064348B (zh) * | 2017-10-12 | 2020-05-05 | 南方科技大学 | 一种基于两点射线追踪的地震走时层析反演方法 |
-
2019
- 2019-09-17 CN CN201910877008.9A patent/CN112526610B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112526610A (zh) | 2021-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dix | Seismic velocities from surface measurements | |
Regone et al. | Geologic model building in SEAM Phase II—Land seismic challenges | |
CN110133715B (zh) | 一种基于初至时差和波形叠加的微地震震源定位方法 | |
Catchings et al. | High-resolution seismic velocities and shallow structure of the San Andreas fault zone at Middle Mountain, Parkfield, California | |
Schaeffer et al. | A low‐velocity zone atop the transition zone in northwestern Canada | |
CN109884710B (zh) | 针对激发井深设计的微测井层析成像方法 | |
CN110501757B (zh) | 点线面三维逐点井深设计方法 | |
CN109799540B (zh) | 基于地质信息约束下的火山岩型铀矿床磁化率反演方法 | |
CN113031068B (zh) | 一种基于反射系数精确式的基追踪叠前地震反演方法 | |
Youngmann | Spatial data structures for modeling subsurface features | |
CN105607119B (zh) | 近地表模型构建方法与静校正量求取方法 | |
CN109655918B (zh) | 地面浅井微地震监测观测台站位置确定方法及系统 | |
CN110737018B (zh) | Vsp地震资料各向异性建模方法 | |
CN111308558B (zh) | 页岩气水平井纵波时差校正方法 | |
Hatayama et al. | Control factors of spatial variation of long-period strong ground motions in the Yufutsu sedimentary basin, Hokkaido, during the M w 8.0 2003 Tokachi-oki, Japan, earthquake | |
CN112526610B (zh) | 一种约束表层建模的三维地震采集激发井深设计方法 | |
CN107179553A (zh) | 基于双压实规律的沙漠区表层静校正方法 | |
CN109188522A (zh) | 速度场构建方法及装置 | |
CN105259577B (zh) | 一种确定地层界面的角度信息的方法及装置 | |
Jin et al. | Comparison of Rayleigh wave dispersion relations from three surface wave measurements in a complex-layered system | |
CN102809761A (zh) | 一种获取煤岩的速度特征的控制方法 | |
CN114200538A (zh) | 地应力方向确定方法、装置、存储介质和计算机设备 | |
Maufroy et al. | Travel time inversion from ground level to gallery: protocol for the characterization of P‐wave seismic signature in a fractured‐porous Urgonian platform at hectometric scale | |
CN112305594A (zh) | 非均质储层的油气分布确定方法及系统 | |
Yang et al. | Crustal and upper mantle structure beneath the Southeastern United States from joint inversion of receiver functions and rayleigh wave dispersion |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |